动载下缝端应力强度因子计算的扩展有限元法-应用数学和力学.PDFVIP

摇 , 34 6 Applied Mathematics and Mechanics摇 摇 2013年6月15 日出版 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 摇 Vol.34,No.6,Jun.15,2013 文章编号:1000鄄0887(2013)06鄄0586鄄12 訫 应用数学和力学编委会,ISSN 1000鄄0887 动动载下缝端应力强度因子 计算的扩展有限元法* 1,2 2 江守燕 杜成斌 ,摇 (1. 河海大学 水利水电学院,南京210098; 2. 河海大学 力学与材料学院,南京210098) 摘要:摇 在扩展有限元法(extendedfinite element methods,XFEM)的理论框架下,重点研究了动荷 载作用下稳定裂纹尖端动态应力强度因子(dynamic stress intensity factors,DSIFs)的求解方法.根 据XFEM 的位移模式,推导了动力XFEM 的支配方程,采用 Newmark 隐式算法进行时间积分. 同 时,提出一种XFEM质量矩阵的集中策略,给出了求解DSIFs 的相互作用积分方法,与静态问题的 相互作用积分方法相比,增加了惯性项的贡献.最后,若干典型算例的计算结果表明:XFEM可以准 确评价稳定裂纹尖端的DSIFs,建议的质量矩阵集中策略是有效的,为得到正确的DSIFs,惯性项的 贡献不可忽略. 关摇 键摇 词:摇 扩展有限元法;摇 动载;摇 裂纹;摇 动应力强度因子;摇 相互作用积分 中图分类号:摇 TB115;O346.1摇 摇 摇 文献标志码:摇 A DOI:10.3879/ j.issn.1000鄄0887.2013.06.005 引摇 摇 言 扩展有限元法(extended finite element methods,XFEM)及其基本思想已被广泛地用于分 [1] [2] 析工程领域中的各种不连续问题,如结构的断裂分析 、非均质材料的热分析 等.在结构的 断裂分析中,很多实际结构多是受动载条件,因此,动力XFEM 的研究显得十分必要,而动应 力强度因子(dynamic stress intensity factors,DSIFs)是判断裂纹在动载下是否发展以及发展方 向的重要参数. [3] [4] Belytschko等 和他的学生 Chen 最早利用XFEM对动力开裂问题进行了研究;随后, [5鄄6] R佴thor佴等 发展了时空扩展有限元法(STX鄄FEM)去处理动载下裂纹生长问题;Menouillard [7鄄9] 等 基于显式时间积分算法研究了改进函数的不同形式对DSIFs及动力裂纹生长的影响,他 们提出通过纠正力(correction force)的方法释放裂尖单元,从而避免裂纹生长过程由于裂尖单 元的突然释放产生的伪应力波;针对显式积分算法的稳定性问题,Menouillard等[10鄄11]研究了积 分步长的取值,同时指出常规有限元处理动力问题时的质量矩阵集中策略不适用于基于显式 算法的动力XFEM,并给出了一种质量矩阵的集中方法;基于谱单元法和显式时间积分算法, * 收稿日期:摇 2012鄄09鄄11;修订日期:摇 2013鄄05鄄16 基金项目:摇 国家自然科学基金资助项目;中国博士后科学基金资助项目(2013M530233) 作者简介:摇 江守燕(1986—),女,安徽舒城人,博士(E鄄mail:syjiang@hhu.edu.cn); 杜成斌(1965—),男,