基于扩展有限元模型的动态应力强度因子计算

中国机械工程

基于扩展有限元模型的动态应力强度因子计算

汪必升1;李毅波1,2;廖雅诗3;李剑1

1.中南大学机电工程学院，长沙，410083
2.中南大学高性能复杂制造国家重点实验室，长沙，410083
3.中南大学轻合金研究院，长沙，410083

出版日期:2019-06-10 发布日期:2019-06-11
基金资助:
国家自然科学基金资助项目(51575535);
中南大学研究生科研创新项目(2018zzts448);
柳州市重点研发计划资助项目(2018BB30501)

Calculation of Dynamic Stress Intensity Factors Based on XFEM Model

WANG Bisheng1;LI Yibo1,2;LIAO Yashi3;LI Jian1

1.College of Mechanical and Electrical Engineering, Central South University,Changsha,410083
2.State Key Laboratory of High Performance and Complex Manufacturing,Central South University,Changsha,410083
3.Light Alloy Research Institute,Central South University,Changsha,410083

Online:2019-06-10 Published:2019-06-11

摘要/Abstract

摘要： 为准确计算基于扩展有限元法（XFEM）的裂纹扩展模型中的应力强度因子，在ABAQUS软件中建立中心裂纹平板和三点弯曲的XFEM模型，采用相互作用积分法，通过用户子程序接口分别实现了Ⅰ型、Ⅱ型断裂模式下裂纹扩展过程中应力强度因子的计算；研究了网格密度与积分半径对XFEM模型应力强度因子计算精度的影响规律，研究结果表明：当网格密度因子为0.012~0.016、相对积分半径为3时，应力强度因子收敛至稳定值，计算误差不超过3%。利用所提方法与程序计算了单边带孔疲劳裂纹扩展试样的动态应力强度因子，试验结果表明：基于Paris理论预测的剩余寿命与疲劳试验结果误差为5.3%，进一步验证了所提方法与程序的正确性。

关键词: 扩展有限元法, 应力强度因子, 相互作用积分法, 网格密度因子, 相对积分半径

Abstract: In order to accurately calculate the stress intensity factors in the crack propagation model based on XFEM, the XFEM model of central crack plate and three-point bending were established in ABAQUS software. The stress intensity factors in the crack propagation processes were realized under mode Ⅰ and Ⅱ fracture modes respectively by using interactive integration method and user subroutine interface. The influences of mesh density and integral radius on the accuracy of stress intensity factor calculation for XFEM were studied. The research results show that when the mesh density factor is as 0.012~0.016 and the relative integral radius is as 3, the stress intensity factors converge to a stable value and the calculation errors are less than 3%. The dynamic stress intensity factors of fatigue crack growth specimens with single side hole were calculated by using the proposed method and programs. The testing results show that the errors between the residual life predicted by Paris theory and the fatigue test results are as 5.3%, which further verifies the correctness of the proposed method and programs.

Key words:
extended finite element method(XFEM), stress intensity factor, interaction integral method, mesh density factor, relative integral radius

中图分类号:

O346

汪必升1;李毅波1,2;廖雅诗3;李剑1. 基于扩展有限元模型的动态应力强度因子计算[J]. 中国机械工程.

WANG Bisheng1;LI Yibo1,2;LIAO Yashi3;LI Jian1. Calculation of Dynamic Stress Intensity Factors Based on XFEM Model[J]. China Mechanical Engineering.

参考文献

[1]文龙飞, 王理想, 田荣. 动载下裂纹应力强度因子计算的改进型扩展有限元法[J]. 力学学报, 2018, 50(3):599-610.
WEN Longfei, WANG Lixiang, TIAN Rong. Accurate Computation on Dynamic SIFs Using Improved XFEM[J]. Chinese Journal of Theoretical and Applied Mechanics, 2018, 50(3):599-610.
[2]BRANCO R, ANTUNES F V. Finite Element Modeling and Analysis of Crack Shape Evolution in Mode-Ⅰ Fatigue Middle Cracked Tension Specimens[J]. Engineering Fracture Mechanics, 2008,75(10):3020-3037.
[3]SAGERESAN N, DRATHI R. Crack Propagation in Concrete Using Meshless Method[J]. CMES Computer Modeling in Engineering & Sciences, 2008, 32(2):103-112.
[4]WEN P H, ALIABADI M H. Dual Boundary Element Method for Modelling Curved Crack Paths[J]. International Journal of Fracture, 2012, 176(1):127-133.
[5]BELYTSCHKO T, BLACK T. Elastic Crack Growth in Finite Elements with Minimal Remeshing[J]. International Journal for Numerical Methods in Engineering, 1999, 45(5):601-620.
[6]NEHAR K C, HACHI B E, CAZES F, et al. Evaluation of Stress Intensity Factors for Bi-material Interface Cracks Using Displacement Jump Methods[J]. Acta Mechanica Sinica, 2017，33(6):1-14.
[7]周博，薛世峰. 基于扩展有限元的应力强度因子的位移外推法[J]. 力学与实践,2017(4):371-378.
ZHOU Bo, XUE Shifeng. Displacement Extrapolation Method for Stress Intensity Factor Based on Extended Finite Method[J]. Mechanics in Engineering, 2017(4):371-378.
[8]GUO Licheng, GUO Fengnan, YU Hongjun, et al. An Interaction Energy Integral Method for Nonhomogeneous Materials with Interfaces under Thermal Loading[J]. International Journal of Solids & Structures, 2012, 49(2):355-365.
[9]ZHENG Anxing, LUO Xianqi. Numerical Study of Quasi-static Crack Growth Problems Based on Extended Finite Element Method[J]. Journal of Shanghai Jiaotong University (Science), 2014, 19(6):736-746.
[10]艾书民,于明,成晓鸣，等. 基于Franc3D软件的三维裂纹扩展分析与应用[J]. 机械强度, 2018,40(1):251-254.
AI Shumin, YU ming, CHENG Xiaoming, et al.Analysis and Application of Three-dimensional Crack Growth Based on Franc3D[J]. Journal of Mechanical Strength, 2018,40(1):251-254.
[11]龚迪光, 曲占庆, 李建雄, 等. 基于ABAQUS平台的水力裂缝扩展有限元模拟研究[J]. 岩土力学, 2016, 37(5):1512-1520.
GONG Diguang, QU Zhanqing, LI Jianxiong, et al. Extended Finite Element Simulation of Hydraulic Fracture Based on ABAQUS Platform[J]. Rock and Soil Mechanics, 2016, 37(5):1512-1520.
[12]ZHU Zhenyang, CHEN Weiming, ZHANG Guoxing. Cracking Propagation of Hardening Concrete Based on the Extended Finite Element Method[J]. Journal of Wuhan University of Technology（Materials Science), 2017(5):1132-1139.
[13]WEI Mingdong, DAI Feng, XU Nuwen, et al. Fracture Prediction of Rocks under Mode Ⅰ and Mode Ⅱ Loading Using the Generalized Maximum Tangential Strain Criterion[J]. Engineering Fracture Mechanics,2017,186:21-38.
[14]BENVENUTI E. An Effective XFEM with Equi-valent Eigenstrain for Stress Intensity Factors of Homogeneous Plates[J]. Computer Methods in Applied Mechanics and Engineering, 2017,4(5):427-454.
[15]王自东，陈少华. 高等断裂力学[M]. 北京：科学出版社,2009.
WANG Zidong, CHEN Shaohua. Advanced Fracture Mechanics[M]. Beijing:China Science Publication, 2009.
[16]申杰斌,唐东林. 一种考虑应力梯度的疲劳寿命预测方法[J]. 中国机械工程, 2017, 28(1):40-44.
SHEN Jiebin, TANG Donglin. Predicting Method for Fatigue Life with Stress Gradient[J]. China Mechanical Engineering, 2017, 28(1):40-44.
[17]LIAO Yashi, LI Yibo, PAN Qing, et al. Residual Fatigue Life Analysis and Comparison of an Aluminum Lithium Alloy Structural Repair for Aviation Applications[J]. Engineering Fracture Mechanics, 2018,194:262-280.