Compute Stress Intensity Factors via Combining Analytical Solutionsaround Crack Tips with Surrounding Numerical Solutions

  • SU Hai-dong ,
  • QI Yong-feng ,
  • GONG Ya-qi
Expand
  • Material and Engineering Structure Department, Yangtze River Scientific Research Institute,Wuhan 430010, China

Received date: 2012-11-13

  Revised date: 2013-06-04

  Online published: 2013-06-04

Abstract

Due to the complex distribution of the displacements and stresses around the crack tip, it is not easy to obtain the Stress Intensity Factor (SIF) with a rapid convergence when using conventional interpolation approaches of numerical methods such as Finite Element Method (FEM). On the basis of Numerical Manifold Method (NMM), a novel method is presented to compute the SIFs via combining analytical solutions with numerical solutions. The Williams expansion is used as the analytical solution, which is formed by applying the constraints of nodal freedoms in the mesh containing the crack tip. High-order polynomial functions are used as the numerical solutions which are connected with the analytical solution via shape functions in the surrounding meshes. Meanwhile, the meshes in NMM need not conform to the physical boundaries including the crack edges, and discontinuous covers are used to allow the cracks arbitrarily align within the meshes, providing the convenience of mesh generation. Numerical example shows the validity of the method. Considering that the Williams expansion is the best approximation for the displacement field around the crack tip, the method has a more rapid convergence than other new methods such as extended Finite Element Method (XFEM).

Cite this article

SU Hai-dong , QI Yong-feng , GONG Ya-qi . Compute Stress Intensity Factors via Combining Analytical Solutionsaround Crack Tips with Surrounding Numerical Solutions[J]. Journal of Changjiang River Scientific Research Institute, 2013 , 30(6) : 83 -89 . DOI: 10.3969/j.issn.1001-5485.2013.06.019

References

[1]范天佑.断裂理论基础[M].北京;科学出版社,2003.(FAN Tian-you. Basic Principle of Fracture [M]. Beijing; Science Press,2003.(in Chinese)
[2] 数值流形方法与非连续变形分析[M].石根华,裴觉民,译.北京;清华大学出版社,1997.(SHI Gen-hua. Numerical Manifold Method (NMM) and Discontinuous Deformation Analysis (DDA) [M]. Translated by PEI Jue-min. Beijing; Tsinghua University Press,1997.(in Chinese)
[3] 李录贤,王铁军. 扩展有限元法(XFEM)及其应用[J]. 力学进展,2005,35(1);5-20.(LI Lu-xian,WANG Tie-jun. The Extended Finite Element Method and Its Applications-A Review [J]. Advances in Mechanics,2005,35(1);5-20.(in Chinese)
[4] 李树忱,程玉民. 考虑裂纹尖端场的数值流形方法[J]. 土木工程学报,2005,38(7);96-101.(LI Shu-chen,CHENG Yu-min. Numerical Manifold Method for Crack Tip Fields[J]. China Civil Engineering Journal,2005,38(7); 96-101.(in Chinese)
[5] WILLIAMS M L. On the Stress Distribution at the Base of a Stationary Crack[J]. Journal of Applied Mechanics,1957,24; 109-114.
[6] KARIHALOO B L,XIAO Q Z. Implementation of HCE on a General FE Mesh for Interacting Multiple Cracks[C]∥Proceedings of the European Congress on Computational Methods in Applied Sciences and Engineering,ECCOMAS,Jyv skyl ; July 24-28,2004.
[7] 彭 俚. 结构断裂分析的Williams广义参数单元[D]. 南宁;广西大学,2008.(PENG Li. Williams Element with Generalized DOF for Structural Fracture Analysis[D]. Nanning; Guangxi University,2008.(in Chinese)
[8] 苏海东,黄玉盈.数值流形方法在流固耦合谐振分析中的应用[J].计算力学学报,2007,24(6);823-828.(SU Hai-dong,HUANG Yu-ying. Application of Numerical Manifold Method in Fluid-Solid Interaction Harmonic Analysis[J]. Chinese Journal of Computational Mechanics,2007,24(6); 823-828.(in Chinese)
[9] 王勖成,邵 敏.有限单元法的基本概念和数值方法[M].北京;清华大学出版社,1988.(WANG Xu-cheng,SHAO Min. Basic Principle of the FEM and Numerical Method[M]. Beijing; Tsinghua University Press,1988.(in Chinese)
[10]中国航空研究院.应力强度因子手册[M]. 北京;科学出版社,1993.(Chinese Institute of Aviation Research. Manual of Stress Intensity Factors[M]. Beijing; Science Press,1993.(in Chinese)
Outlines

/