用边界轮廓法求解断裂力学L积分.pdfVIP

第 19 卷 第 4 期 郑州工业大学学报 Vol119 No14 1998 年 12 月 Journal of Zhengzhou University of Technology Dec1 1998 L 陈颂英 唐委校 冷 霞 郑 雯 ( 山东工业大学化工系, 济南, 250061) 证明断裂力学L 积分方程的被积函数的 度等于零, 将面积分转化为线积分, 使求解 问题的维数降低两维。利用边界轮廓法的结果, 使平面断裂问题L 积分的求解转化为 边界点的位移和面力线性迭加, 避免了求解数值积分。 断裂力学; 边界元法; L 积分 TV 313 L 积分作为断裂力学的不变积分之一, 体现了断裂力学的能量守恒定 , 描述了裂纹尖 端应力奇异性的不变特性, 积分回路中只要包含裂纹尖端的应力奇异性, 其值均与路径无 关。Rice[ 1] 将其解释为绕坐标轴转动时的能量释放率, Hemnann[2] 给出了在特定条件下, L 积分与 J 积分的关系式, Sun Shu un[ 3] 推导了 L 积分的守恒性。目前, 关于 L 积分研究的 2 文献报导还较少, 尚未涉及其计算问题。 本文推证 L 积分方程的被积函数的散度等于零, 将被积函数表达为某一矢量场的旋 度, 再代入积分方程中, 应用 Stokes 公式, 将平面断裂问题的L 积分转化为边界点积分势函 数的点值计算, 从而使求解问题的维数降低两维。利用陈颂英[4] 所建立的边界轮廓法的结 果, 将 L 积分的计算转化为已求得的边界点的位移和面力的代数和, 避免了求解奇异积分。 文中给出平面问题的算例。 1 传统的无体力边界积分方程[ 5] 为: aCc (P) u (P) = Uik(P, Q) Rij ( Q) - (P, Q) ui(Q) e # dS (1) ik i Q E ijk j 5B 式中 Cik 为角张量,P 为源点, Q 为场点, ui 和Rij 分别为位移和应力分量, ej 为 Cartesian 单位矢量,5 B 为弹性体边界, dS 为边界变元矢量, U 和E 为 Kelvin 核函数。 ik ij k 将方程(1) 中的被积函数记为: g U (P, Q) R (Q) - (P, Q) u ( Q) F = ik ij i e (2) k E ij k j 对上式在场点 Q 处求散度, 由文献[4] 知: ¨ # F = 0 (3) Q k 由场论知, 存在矢量 V= 5 e ( 5 称为位移势) , 满足 F = ¨ @ V, 对于二维平面问 k 3 k k Q 题, 将弹性体边界离散为 N 个单元, 原积分方程(1)化为: N