无网格GALERKIN方法及裂纹扩展数值模拟研究

【摘要】： 本文以裂隙岩体为研究背景,利用无网格Galerkin方法和偶应力理论下的无网格Galerkin方法,根据本文提出的最大周向拉应力裂纹扩展准则,对二维裂纹扩展和三维裂纹扩展问题进行了数值模拟研究。 首先推导了二维无网格Galerkin方法、偶应力理论下的无网格Galerkin方法及三维无网格Galerkin方法的基本公式。然后根据最小势能原理,导出了弹性体的离散方程。以经典的二维和三维弹性力学问题为算例研究了无网格方法的程序实施和几个重要参数的合理配置。采用罚参数法来释放位移约束,引入边界条件;基于节点控制的思想,确定动态影响半径,那么在裂纹前沿节点加密区域适当减少了最小二乘的计算量,可以适当提高无网格方法的计算效率。 利用空间解析几何理论,建立了裂纹前沿局部坐标系,计算了局部坐标和全局坐标的变换矩阵,导出了局部坐标中周向正应力与全局坐标中应力分量的变换公式,为裂纹曲面(线)的描述和裂纹扩展数值模拟做好了充分的数学准备工作。 在研究裂纹扩展时,采用可视准则,体现裂纹对节点之间、节点与高斯积分点之间的隔离作用。 选择合适的裂纹扩展准则,确定裂纹扩展条件、扩展方向和扩展步长是研究裂纹扩展的关键。考虑到目前常用的二维裂纹扩展准则存在缺陷和没有三维裂纹扩展准则,在其它学者的思想基础上,本文提出了一个新的裂纹扩展准则——最大周向拉应力裂纹扩展准则。该准则不仅给出了裂纹扩展方向、扩展条件,还给出了确定扩展步长的依据,摒弃了计算应力强度因子带来的误差和麻烦,尤其是三维裂纹问题,克服了靠经验给定扩展步长的弊端。 结合本文提出的最大周向拉应力裂纹扩展准则,考虑到三维裂纹的扩展曲面可能出现弯折、扭曲现象,把裂纹扩展曲面分成若干个片断面,并将每个裂纹片断面划分成三角形表示,实现了裂纹曲面的数学描述。 在裂纹扩展的数值模拟中,为了识别裂纹形态和提高模拟精度,在裂纹面两侧增加了附加节点,并对处在裂纹扩展段或扩展片断面上的节点实施了劈裂技术,分别置于两侧的裂纹表面上。 研究无水压平面裂纹扩展数值模拟,并与其它数值计算方法的结果和相关的试验结果进行了对比,得到了一致的结论,确保了最大周向拉应力裂纹扩展准则的有效性。在此基础上,分别研究了裂隙水压和偶应力对裂纹扩展的影响,得到了一些新结论。 由于三维裂纹扩展问题的复杂性,目前对三维裂纹扩展数值模拟的研究成果较少。用传统的计算方法仅对币形裂纹的Ⅰ型扩展进行了数值模拟。本文以三维无网格Galerkin方法为基础,采用最大周向拉应力裂纹扩展准则,首先对椭圆形裂纹的Ⅰ型扩展进行了数值模拟,数值模拟结果与试验结果一致,再一次证实了三维无网格Galerkin方法程序实施技术的正确性、最大周向拉应力裂纹扩展准则的正确性和三维裂纹扩展曲面数学描述的有效性。在此基础上,对三维表面斜置半圆形裂纹的扩展进行了数值模拟,了解了复杂三维裂纹不同扩展阶段的扩展方式、扩展曲面的形态,以及空间方位对三维裂纹扩展模式的影响,填补了国内在三维复杂裂纹扩展数值模拟方面的空白。 通过三维裂纹的数值模拟可以清楚地了解三维裂纹扩展的全貌,弥补了三维裂隙断裂试验与理论研究的不足,为进一步研究三维裂隙断裂的本质提供一定的参考价值。