动载下裂纹应力强度因子计算的改进型扩展有限元法

相较于常规扩展有限元法(extended finite element method,XFEM),改进型扩展有限元法(improved XFEM)解决了现有方法线性相关与总体刚度矩阵高度病态问题,在数量级上提升了总体方程的求解效率,克服了现有方法在动力学问题中的能量正确传递、动态应力强度因子数值震荡、精度低下问题.本文基于改进型XFEM,采用Newmark隐式时间积分算法,重点研究了动载荷作用下扩展裂纹尖端应力强度因子的求解方法,与静力学方法相比,增加了裂纹扩展速度项与惯性项的贡献.通过数值算例研究了网格单元尺寸、质量矩阵、时间步长、裂尖加强区域、惯性项、扩展速度项及相互作用积分区域J-domain的网格与单元尺寸对动态应力强度因子求解精度的影响,验证了改进型XFEM计算动态裂纹应力强度因子方法的有效性.针对文献中具有挑战性的"I型半无限长裂纹先稳定后扩展"问题,改进型XFEM给出目前为止精度最好的动态应力强度因子数值解.     

非连续问题中单元分割的模板方法

扩展有限元法 (extended finite element method, XFEM) 因具有裂纹几何独立于模拟网格、裂纹扩展时无需网格重分重映、计算精度高等优点,成为裂纹分析的主流数值方法之一. 但该方法在工程实践中存在单元被裂纹分割的几何困难 —— 现有精确几何分割方法实现复杂、计算量大、鲁棒性差. 为克服这一困难, 本文提出一种基于单元水平集的模板分割方法, 用于非连续单元子剖分和数值积分. 首先, 遍历单元水平集值所有形态并建立标准单元分割模板库; 然后, 根据单元水平集值, 对非标准单元进行形态查询和模板插值; 最后, 套用标准单元分割模板实现单元高效分割和子剖分. 将该方法与常规XFEM、改进型XFEM进行结合,从而应用于孔洞、夹杂、裂纹等非连续问题分析中. 算例分析表明, 本文提出的模板分割方法具有较高计算精度. 由于不引入复杂几何操作, 该模板分割方法同时具有较高计算效率和鲁棒性, 故可为XFEM类方法在实际工程应用中提供有效支撑.      

动载下裂纹应力强度因子计算的改进型扩展有限元法

相较于常规扩展有限元法(extended finite element method, XFEM), 改进型扩展有限元法(improved XFEM) 解决了现有方法线性相关与总体刚度矩阵高度病态问题, 在数量级上提升了总体方程的求解效率, 克服了现有方法在动力学问题中的能量正确传递、动态应力强度因子数值震荡、精度低下问题. 本文基于改进型XFEM, 采用Newmark 隐式时间积分算法, 重点研究了动载荷作用下扩展裂纹尖端应力强度因子的求解方法, 与静力学方法相比, 增加了裂纹扩展速度项与惯性项的贡献. 通过数值算例研究了网格单元尺寸、质量矩阵、时间步长、裂尖加强区域、惯性项、扩展速度项及相互作用积分区域J-domain的网格与单元尺寸对动态应力强度因子求解精度的影响, 验证了改进型XFEM计算动态裂纹应力强度因子方法的有效性. 针对文献中具有挑战性的 "I 型半无限长裂纹先稳定后扩展"问题, 改进型XFEM给出目前为止精度最好的动态应力强度因子数值解.      

改进型XFEM进展

扩展有限元法(XFEM)在诞生后的十几年时间里,引起学术界和工业界的广泛关注,并已经成为目前裂纹分析的主流数值方法。然而,在实际应用中该法一直受到两方面的困扰,(1)总体方程高度病态;刚度阵条件数随网格尺寸呈h-6变化(普通有限元为h-2)。(2)裂尖强化插值由于能量一致性问题无法直接推广应用于动力学计算。前者表现在XFEM稳态问题的迭代求解收敛慢或难以收敛,后者长期以来导致XFEM裂纹扩展动力学计算实施困难。本文认为XFEM目前遇到的种种困难,均与单位分解引入的额外自由度相关。为此,提出了无额外自由度的单位分解插值格式,基于此格式,进一步构造出改进型扩展有限元方法。改进型XFEM具有如下特点,(1)可以消除原有XFEM的线性依赖性和总体方程病态的问题。(2)避免动力学问题中额外自由度引起的质量集中、零临界时间步长问题以及裂纹扩展过程中的能量一致性问题。本文结合静动力学测试问题综述上述改进。