力学学报 ›› 2018, Vol. 50 ›› Issue (3): 599-610.doi: 10.6052/0459-1879-18-030

• 固体力学 • 上一篇    下一篇

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

文龙飞, 王理想, 田荣*   

  1. 中物院高性能数值模拟软件中心, 北京 100088;
    北京应用物理与计算数学研究所, 北京 100088
  • 收稿日期:2018-01-30 出版日期:2018-05-15 发布日期:2018-06-11
  • 通讯作者: 田荣, 研究员, 主要研究方向: 计算力学与高性能计算. E-mail: tian_rong@iapcm.ac.cn
  • 基金资助:
    国家重点研发计划(2016YFB0201002, 2016YFB0201004),国家自然科学基金(11472274, 91530319),科学挑战专题(JCKY2016212A502)资助项目.

ACCURATE COMPUTATION ON DYNAMIC SIFS USING IMPROVED XFEM

Wen Longfei, Wang Lixiang, Tian Rong*   

  1. CAEP Software Center for High Performance Numerical Simulation, Beijing 100088, China;
    Institute of Applied Physics and Computational Mathematics, Beijing 100088, China
  • Received:2018-01-30 Online:2018-05-15 Published:2018-06-11

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

关键词: 改进型XFEM|无额外自由度单位分解|扩展裂纹|动态应力强度因子|线弹性断裂力学

Abstract: Compared to the standard extended finite element method (XFEM), the improved XFEM overcomes, in theory, the linear dependence and the ill-conditioning issues, and improves in magnitude of orders the efficiency of linear system solve. In particular, in dynamic crack propagation problems, thanks to the exclusion of the extra dofs on crack tip enriched nodes, the new method eliminates the issue of energy inconsistency or the inconservitive energy transfer caused by dof dynamics, and provides more accurate dynamic stress intensity factors (DSIFs) with much less numerical oscillation. To the best of our knowledge, numerical solution on DSIFs for crack propagation under dynamic loading remains engineeringly unsatisfied. In this paper, the extra-dof free XFEM is extended to implicitly dynamic crack propagation problems-still a remaining difficulty for the current XFEMs. The implicit Newmark algorithm is used for time discretization. A dynamic interaction integral method is employed for DSIFs for both stationary and moving cracks under dynamic loading. Compared with the interaction integral method for static cracks, the dynamic method considers the effects of crack growth speed and inertia terms. The paper investigates in detail the influences of element size, mass matrix formulation, time step size, crack tip enriched zone size, inertia term, crack growth speed, and J-domain mesh/cell sizes of the interaction integral. Numerical tests show satisified accuracy and efficiency of the new method for dynamic crack problems. In particular on the challenging benchmark problem "Mode I semi-infinite stationary and then moving crack", the new improved XFEM provides the best DSIFs up to the time of publication.

Key words: improved XFEM|extra-dof free partition of unity approximation|crack growth|dynamic stress intensity factor|linear elastic fracture mechanics

中图分类号: 

  • TB115