Prediction of Crack Initiation Direction for Inclined Crack Under Biaxial Loading by Finite Element Method

Document Type : Research Paper


1 Department of Mechanical Engineering, College of Technology, G. B. Pant University of Agriculture and Technology

2 Mechanical Engineering Department, National Institute of Technology, Jamshedpur, Jharkhand

3 Deputy Manager (Design) Air Craft Upgrade Research and Design Centre, Hindustan Aeronautics Limited


This paper presents a simple method based on strain energy density criterion to study the crack initiation angle by finite element method under biaxial loading condition. The crack surface relative displacement method is used to eliminate the calculation of the stress intensity factors which are normally required. The analysis is performed using higher order four node quadrilateral element. The results by finite element method are compared with DET (determinant of stress tensor criterion) and strain energy density criteria. Finite element results are in well agreement with the experimental and analytical results.


[1] Erdogan F., Sih G.C., 1963, On the crack extension in plates under plane loading and transverse shear, Journal of Basic Engineering 85: 19-27.
[2] Sih G.C., 1974, Strain energy density factor applied to mixed mode crack problems, International Journal of Fracture 11: 305-321.
[3] Theocaris P.S., Andrianopoulos N.P., 1982, The T-Criterion applied to ductile fracture, International Journal of Fracture 20: R125-130.
[4] Hussain M.A., Pu S.L., Underwood J.H., 1974, Strain energy release rate for a crack under combined mode I and mode II, Fracture Analysis, ASTM STP 560: 2-28.
[5] Papadopoulos G.A., 1988, Crack initiation under biaxial loading, Engineering Fracture Mechanics 29: 585-598.
[6] Shih C.F., 1974, Small scale yielding analysis of mixed mode plane strain crack problem, Fracture Analysis, ASTM, STP 560: 187-210.
[7] Obta M., 1984, On stress field near a stationery crack tip, Mechanics of Materials 3: 35-243.
[8] Saka M., Abe H., Tanaka S.,1986, Numerical analysis of blunting of crack tip in a ductile material under small scale yielding and mixed mode loading, Computational Mechanics 1: 11-19.
[9] Dong P., Pan J., 1990, Plane strain mixed mode near tip fields in elastic perfectly plastic solids under small scale yielding condition, International Journal of Fracture 45: 243-262.
[10] Dong P., Pan J., 1990, Plane stress mixed mode near tip fields in elastic perfectly solids, Engineering Fracture Mechanics 37: 43-57.
[11] Sedmak A., 1984, Finite element evaluation of fracture mechanics parameter using rapid mesh refinement, Advance in Fracture Research, 1095-1106.
[12] Guydish Jacob J., Fleming J.F., 1978, Optimization of the finite element mesh for the solution of fracture problems, Engineering Fracture Mechanics 10:31-42.
[13] Mahanty D.K., Maiti S.K., 1990, Experimental and finite element studies on mode I and mixed mode (I and II) stable crack growth-I, Engineering Fracture Mechanics 37: 1237-1250.
[14] Mahanty D.K., Maiti S.K.,1990, Experimental and finite element studies on mode I and mixed mode (I and II) stable crack growth-II, Engineering Fracture Mechanics 37: 1251-1275.
[15] Ju S.H., 2010, Finite element calculation of stress intensity factors for interface notches, Computer Methods in Applied Mechanics and Engineering 199 (33-36): 2273-2280.
[16] Réthoré J., Roux S., Hild F., 2010, Mixed-mode crack propagation using a Hybrid Analytical and extended Finite Element Method, Comptes Rendus Mécanique 338(3): 121-126.
[17] Benrahou K.H., Benguediab M., Belhouari M., Nait-Abdelaziz M., Imad A., 2007, Estimation of the plastic zone by finite element method under mixed mode (I and II) loading, Computational Materials Science 38(4): 595-601.
[18] Ghorbanpoor A., Zhang J., 1990, Boundary element analysis of crack growth for mixed mode center slant crack problems, Engineering Fracture Mechanics 36(5): 661-668.
[19] Wang J., Chow C.L., 1989, Mixed mode ductile fracture studies with non proportional loading based on continuum damage mechanics, Journal of Engineering material and Technology III, 204-209.
[20] Sun Y.J., Xu L.M., 1985, Further studies on crack tip plasticity under mixed mode loading, in: Proceedings of 1985 Spring Conference, Experimental Mechanics Publication, Las Vegas, 20-25.
[21] Lee K.Y., Lee J.D., Liebowitz H., 1997, Finite element analysis of slow crack growth process in mixed mode fracture, Engineering Fracture Mechanics 56: 551-577.
[22] Seibi A.C., Zamrik S.Y., 2003, Prediction of crack initiation direction for surface flaws under biaxial loading, Journal of Pressure Vessel and Technology 125: 65-75.
[23] Ling L. H., Woo C.W., 1984, On angle crack initiation under biaxial loading, Journal of Strain Analysis 19(1): 51-59.
[24] Shlyannikov V.N., Kislovaa S.Y., Tumanova A.V., 2010, Inclined semi-elliptical crack for predicting crack growth direction based on apparent stress intensity factors, Theoretical and Applied Fracture Mechanics 53(3): 185-193.
[25] Kibler J.J., Roberts R.., 1970, The effect of biaxial stresses on fatigue and fracture, Journal of Engineering for Industry 92: 727-734.
[26] Hilton P.D., 1973, Plastic intensity factors for cracked plates subjected to biaxial loading, International Journal of Fracture 9: 149-156.
[27] Liebowitz H., Lee J.D., Eftis J., 1978, Biaxial load effects in Fracture mechanics, Engineering Fracture Mechanics 10: 315-335.
[28] Hafeele P. M., Lee J.D., 1995, The constant stress term, Engineering Fracture Mechanics 50: 869-882.
[29] Yan X.., 2006, Multiple crack fatigue growth modeling by displacement discontinuity method with crack-tip elements, Applied Mathematical Modeling 30: 489-508.
[30] Kamaya M.., 2003, A crack growth evaluation method for interacting multiple cracks. JSME International Journal, Series A 46(1): 15-23.
[31] Masayuki Kamaya, 2008, Growth evaluation of multiple interacting surface cracks. Part I: Experiments and simulation of coalesced crack, Engineering Fracture Mechanics 75:1336-1349.
[32] Masayuki Kamaya, 2008, Growth evaluation of multiple interacting surface cracks. Part II: Growth evaluation of parallel cracks, Engineering Fracture Mechanics 75: 1350-1366.
[33] Chow C.L., Lau K.J., 1976, Finite element analysis of cracked bodies to determine stress intensity factors, Journal of Strain Analysis 11, 18-25.
[34] Woo C.W., Kuruppu M.D., 1982, Use of finite element method for determining stress intensity factors with a conic-section simulation model of crack surface, 20: 163-178.
[35] Sneddon I.N., Lowengrub M., 1970, Crack Problems in Classical Theory of Elasticity, John Wiley, New York.
[36] Gope P.C., 2002, Investigation into the crack propagation under static and fatigue loading, PhD Thesis, Ranchi University, Ranchi.
[37] Bahuguna S., Kothari A., Awasthi P., 2001, Study on Mixed Mode Crack Propagation under Uniaxial Tensile Loading, B. Tech Project report, Govind Ballabh Pant University of Agriculture and Technology.
[38] Golos K., Wasiluk B., 2000, Role of plastic zone in crack growth direction under mixed mode loading, International Journal of Fatigue 102(4): 341-353.
[39] Seibi A.C., Zamrik S.Y., 1997, Prediction of crack initiation direction for surface flaws under biaxial loading, in: 5th International Conference on Biaxial/Multiaxial Fatigue and Fracture Cracow ‘97 Poland 2, 611-622.
[40] Ewing P.D., Williams J.G., 1974, The fracture of spherical shells under pressure and circular tubes with angled cracks in torsion, International Journal of Fracture 10 (4): 537-544.