1
无网格法理论及MATLAB程序
1.9.4.3 7.4.3 受压裂纹应力强度因子
7.4.3 受压裂纹应力强度因子

如图7.9(a)所示,裂纹试样尺寸为50 mm×50 mm,,裂纹半长为a=5mm,裂纹中心位于试样中心。材料参数取为E=5 GPa,泊松比v=0.35,摩擦系数μ=0。试样顶部和底部施加均匀的单轴压缩荷载σ=10000 Pa。节点离散见图7.9(a)。

在裂纹面光滑情况下,该问题的应力强度因子的理论解[134]

为正则化应力强度因子,采用相互作用积分计算了不同裂纹倾角θ 对应的应力强度因子。图7.9(b)给出了本文方法计算的应力强度因子与θ 的变化关系,并与理论结果进行了比较。

图7.9 受压裂纹试样及其应力强度因子计算结果