线性无关高阶数值流形法在断裂力学中的应用

摘要

数值流形方法(NMM)实现了对连续和非连续问题的统一求解,非常适合应用于求解断裂力学问题。但传统高阶NMM,例如采用一次多项式作为其局部位移函数时,存在线性相关问题,这又在一定程度上限制了NMM的进一步发展和应用。因此,通过在物理片上引入一种新的局部位移函数以及用于模拟裂纹尖端应力场奇异性的位移函数,建立新的NMM求解体系,尝试消除线性相关问题和求解线弹性断裂力学问题。数值算例结果表明:(1)该方法有效地解决了线性相关问题。(2)对于典型的线弹性断裂力学问题,即便在网格密度较小时也能够精确地计算出裂纹尖端的应力强度因子。(3)研究区域内插值点处的应力是连续的。(4)定义在非奇异物理片上自由度具有明确的物理含义,且第3到5个自由度恰好是所对应插值点处的应变分量,可直接获得此处的应力,减少了计算量。最后,所建议方法可以很容易地推广到其他基于单位分解理论的方法中。