利用有限元线法奇异线映射技巧计算应力强度因子.pdfVIP

《工程力学》增刊1999年 ·085· 利用有限元线法奇异线映射技巧 计算应力强度因子· 徐永君 袁驷 中国科学院力学研究所．北京．100080) (滴华大学土木系，北京．100084) 提要本文讨论了有限元线法的奇异线映射技巧在具有任意切口、任意切口边界条件的多材料反平 面断裂问题和平面断裂问题的应用，不仅使得计算速度和计算精度有显著提高，同时使得位移的导数项易 于较精确地给出。利用位移的导数项与特征展开项的低阶特征系数问的关系，可直接将低阶特征系数给 出，并确定应力强度因子。 关键词：育限元线法，奇异线映射，应力强度因子。 1引言 利用奇异线映射技巧求解奇异性问题时可以改变其求解性态，特别是使奇异性问题的 求解精度和求解效率具有较大幅度的提高【1．2】。结合ODE求解器口。5】的自适应等优势使得有 Element of 限元线法(FiniIe MetIlod LineS，简称FEMOL)计算奇异问题具有更多的突出特 点，突破了其它一些方法(有限元法)的仅对裂纹问题才能采用四分之一布点的局限性。区 间映射技巧和FEMoL布线的灵活性使得可以直接采用更为灵活通用的变量代换映射技 巧，对任意奇异性问题都能方便地采用奇异线映射技巧求解，计算量大为减少和求解精度 明显提高。因此可以在奇异点附近的范围内更出色地描述奇异问题对应的奇异场。经过奇 异线映射变换后，可将其含有奇异项的级数(奇异场，如应力场)转换为不含奇异项的级数， 使计算奇异问题转换为计算非奇异问题，可直接得到低阶特征根对应的特征展开系数．从 而利用特征项系数与应力强度因子的关系直接有效地计算应力强度因子。 一’～一、 N ，b，力、 ≥≮专产一1- l ●y J 、一二哆‘ 圈IN种材料二雏切口问题 图2FEMOL的径向布线同格 本文利用FEMoL奇异线映射技巧计算了反平面断裂问题和平面断裂问题。不仅对采用 国家杰出基金和博上后基金资助项目 馀永君．男．1964lO出生．博士 ‘086‘ ‘工程力学》增刊1999年 奇异线映射技巧和不采用奇异线映射技巧的求解效率进行了比较，同时也给出了应力强度 冈子的计算算例。 2断裂问题的用旺奇异线映射 考虑如图l所示的多(Ⅳ)种材料二维切口问题，为分析讨论方便，除特别需要我们一 般略去表示不同材料域的角标。对于反平面切口问题的位移“，，西和平面切口问题的应力 势函数烈r，研，熟知的willi锄s特征展开可分别表示为 (-) 反平面：以r，功= ∑吼，‘Z(印 JtU．1．ij．一 (ii) 平面：船，功= ∑8，n1E@ (1) j=m1．2^一 其中Z(印和E(研均为^。对应的经过归一化的特征向量，其表达形式分别为 Z(研=一，cosz，口+雪，sin五．口 f(印=』．sin(五．+1)口+B，co《五，+1)口阜c，sjn(A，一1)口+D，co《2。一1)一 对于平面切口问题特征报五，对应的位移部分也为，^阶的项 ”rt唑业掣 f2 巡 0，l，2，3…… (2) ¨97=r-1 篙掣 一 特征根五，相对应的应力部分均为，^“阶的项 I ∽姗：E：i艺≯圳^z^…．． 仃厅 =，‘。口』【E’(印+以