断裂问题的扩展有限元法研究  

Study of fracture problem with extended finite element method

作　　者：茹忠亮[1] 朱传锐[1] 张友良[2] 赵洪波[1] 

机构地区：[1]河南理工大学土木工程学院,河南焦作454010 [2]中国科学院武汉岩土力学研究所岩土力学与工程国家重点实验室,武汉430071

出　　处：《岩土力学》2011年第7期2171-2176,共6页Rock and Soil Mechanics

基　　金：国家自然科学基金(No.50804014;No.41072224);教育部新世纪优秀人才计划(No.NCET-08-0662);河南省教育厅自然科学研究计划(No.2008A44005);河南省高等学校青年骨干教师资助计划(No.2009GGJS-037)

摘　　要：扩展有限元(extended finite element method,XFEM)是近年来发展起来的、在常规有限元框架内求解不连续问题的有效数值计算方法,其基于单位分解的思想,在常规有限元位移模式中加入能够反映裂纹面不连续性的跳跃函数及裂尖渐进位移场函数,避免了采用常规有限元计算断裂问题时需要对裂纹尖端重新加密网格造成的不便。在推导扩展有限元算法的基础上,分析了应力强度因子的J积分计算方法及积分区域的选取。采用XFEM对I型裂纹进行了计算,有限元网格独立于裂纹面,无需在裂纹尖端加密网格;分析了积分区域、网格密度对应力强度因子计算精度的影响,指出了计算应力强度因子的合适参数,验证了此方法的可靠性和准确性。The extended finite element method （XFEM） is a numerical method for modeling discontinuities within the classical finite element framework developed in recent years. In XFEM, jump functions and asymptotic crack-tip fields were added to enrich the finite element displacement using the framework of partition of unity to model the discontinuities of cracks. A key advantage of XFEM is that in fracture problems the finite element mesh does not need to be updated and remeshed. Based on the algorithm of XFEM, the J integral methods which calculate the stress intensity factor（SIF） are introduced. A numerical case with model I fracture is analyzed by XFEM. The results illustrate that the FEM mesh is independent of the crack interface. The major factors which influence the accuracy of SIF are discussed, such as integral domain and mesh density. Proper parameters to calculate the SIF are given; and the feasibility and accuracy of XFEM are validated.

关 键 词：扩展有限元 裂纹 形函数 数值方法 

分 类 号：TU443[建筑科学—岩土工程]