• CN:11-2187/TH
  • ISSN:0577-6686

机械工程学报 ›› 2016, Vol. 52 ›› Issue (10): 91-96.doi: 10.3901/JME.2016.10.091

• 材料科学与工程 • 上一篇    下一篇

扫码分享

基于最大能量释放率原理的裂纹扩展算法的改进*

朱有利, 侯帅, 王燕礼, 孙寒骁   

  1. 装甲兵工程学院装备维修与再制造工程系 北京 100072
  • 出版日期:2016-05-15 发布日期:2016-05-15
  • 作者简介:朱有利,男,1962年出生,教授,博士研究生导师。主要研究方向为表面强化技术和疲劳与断裂。E-mail:zhuyl2011@sina.com;侯帅(通信作者),男,1989年出生,博士研究生。主要研究方向为疲劳断裂与寿命预测。 E-mail:houshuai2012@sina.com

Improvements of the Crack Extension Algorithm Based on the Maximum Energy Release Rate Criterion

ZHU Youli, HOU Shuai, WANG Yanli, SUN Hanxiao   

  1. Faculty of Equipment Maintenance and Remanufacturing Engineering, Academy of Armored Forces Engineering, Beijing 100072)
  • Online:2016-05-15 Published:2016-05-15

摘要: 在基于最大能量释放率判据(Maximum energy release rate, MERR)的复合型裂纹扩展分析中,因需要对多个虚拟裂纹扩展方向进行试算和迭代运算,这会大大增加计算量,特别是为保证计算精度不得不选取很小的试算转角增量时,该问题会很突出。为此,提出分别基于有限差分法(Finite difference method, FDM)和最大周向应力判据(Maximum tangential stress, MTS)的裂纹转角预测方法,以减小计算量或提高计算精度。推导相关计算公式和算法,利用FORTRAN与ANSYS参数化设计语言(ANSYS parameter design language, APDL)语言混合编程的方法实现了基于ANSYS有限元环境下的裂纹扩展自动分析,并通过对带孔薄板的疲劳裂纹扩展分析,验证了所提出方法的有效性。分析发现,FDM的预测误差小于MTS的预测误差。MTS预测的转角误差曲线与转角曲线的一阶导数有相同的变化规律,而FDM预测转角的误差曲线与转角曲线的二阶导数有相同的变化规律。根据MERR和MTS判据得到的KII值趋近于零,这与局部对称判据(Local symmetry, LS)一致,但MERR和LS判据的计算结果更接近。

关键词: 裂纹扩展方向, 有限差分法, 最大能量释放率, 最大周向应力

Abstract: For a mixed-mode crack propagation analysis according to the maximum energy release rate(MERR) criterion, test and iterative calculation in multiple virtual crack extension directions significantly increases computation time, especially, when using smaller trying angle increments to insure accuracy of the analysis. To overcome the inconvenience, algorithms based on the finite difference method(FDM) and the maximum tangential stress(MTS) criterion are proposed for crack kinking angle prediction, so as to reduce computation time or to improve accuracy. Relevant formulation and algorithms are derived and verified for effectiveness via fatigue crack propagation analysis of a thin plate with holes. Automatic crack extension analysis by mixed language programming with FORTRAN and ANSYS parameter design language (APDL) in the ANSYS software is realized. Comparative analysis revealed that kinking angle prediction errors from FDM are smaller than those by MTS. The kinking angle errors resulted from MTS follow the same rule with the second derivatives of the kinking angle curve, while those from FDM follow the same fluctuation with the first derivatives of the kinking angle curve. It is found that KII tends to be zero both for the MERR criterion and for the MTS criterion, which is consistent with the the local symmetry criterion(LS), while the MERR criterion seems more akin to the LS criterion in view of the KII value predicted.

Key words: crack kinking direction, finite difference method, maximum energy release rate, maximum tangential stress

中图分类号: