论文部分内容阅读
在工程领域中,常发生起源于断裂或终止于裂纹扩展的灾难性破坏事故。从离散裂缝断裂力学角度加以数值模拟研究结构的断裂破坏过程,对于提高工程设计水平有着重要的意义。扩展有限元法虽然在模拟裂纹扩展时无需进行网格重构,但需预先获得裂尖解析形式的渐进解去构造复杂的增强函数,在形成单元刚度阵时常需采用特殊的数值积分技术；比例边界有限元法作为一种半解析的方法在求解无限域和裂尖奇异性问题方面优势明显,两者衔接于有限元法理论内,可建立一种结合二者优势的断裂耦合数值模型。本文首先从虚功原理出发,利用边界的位移协调与力平衡,将两种方法结合,建立了扩展比例边界有限元法,该方法在裂纹贯穿单元采用Heaviside增强函数描述裂纹面两侧的不连续位移,在裂尖单元则采用半解析的比例边界有限元描述奇异应力场。其次,基于扩展比例边界有限元法推导了裂纹面有任意分布荷载、正交各向异性材料以及双材料界面断裂等复杂问题的位移场和应力强度因子的求解公式。随后,将水平集的思想集入扩展比例边界有限元法中,进一步提升了该方法处理裂纹扩展问题的能力。最后,将扩展比例边界有限元法运用于实际结构,对带预置裂缝的Koyna重力坝进行了超载水头作用下的线弹性断裂分析,得到了坝体裂缝扩展过程和位移响应。在数值算例中,通过编写完整的Matlab分析计算程序,求解了混合型裂纹、裂纹面受任意分布荷载作用的半无限大板、单边缺口的三点弯曲梁、四点剪切梁以及带预置裂纹的Koyna重力坝等断裂问题。通过算例可发现扩展有限元法具有如下特点：1)无需预先知道裂尖位移场或应力场的渐进展开形式；2)在裂尖不需要引入增强函数,无需采用特殊的数值积分技术,直接生成裂尖刚度阵；3)对多种应力奇异类型均可根据定义由统一的表达式直接求解广义应力强度因子；4)在裂纹扩展时,拓扑关系简单,不需增加新的节点和单元,超单元重建的工作量很小。为增强算法的工程适用性,通过参数敏感性分析还可得出如下结论：1)裂尖超单元的层数会影响裂尖断裂特征的计算精度。随着超单元边界节点的增加,计算误差会迅速收敛在工程允许的范围,通常只需采用3层以上的网格,就能获得很好的计算效果。2)算法的网格依赖性较低,特别是当网格尺寸小于扩展步长时,网格大小对裂纹扩展形态的影响就可以忽略；3)扩展步长会影响裂纹的扩展过程,过大的步长会导致计算结果的不稳定,并使裂纹路径出现明显的“锯齿”现象,为了获得稳定收敛的解,可将裂纹扩展步长取为裂缝处构件宽度的1/10～1/25。