Simulating Fatigue Crack Growth in Spiral Bevel Pinion

Abstract : New results for predicting crack trajectory and fatigue life for a spiral bevel pinion using the Finite Element Method (FEM) are reported. The predictions presented are based on linear elastic fracture mechanics combined with the FEM, incorporating plasticity induced fatigue crack closure and moving gear tooth loads. The analyses were carried out using a parallel FEM solver, which calculates stress intensity factors using equivalent domain J-integral method. Fatigue life predictions were made based on a modified Paris model incorporating crack closure. To obtain a more detailed understanding of the contact between a cracked pinion tooth in mesh with an uncracked gear tooth, three-dimensional contact analyses were performed on a spiral bevel gear set incorporating a crack. The goal in carrying out these analyses was to capture the redistribution of contact loads due to crack growth. Results of these analyses showed the expected trend of decreasing tooth loads carried by the cracked tooth with increasing crack length. It was also showed that this decrease in contact loads bad an impact on the stress intensity factor values and therefore would also affect the crack trajectory and fatigue life predictions.

[1]  Gerd Heber,et al.  Parallel FEM Simulation of Crack Propagation - Challenges, Status, and Perspectives , 2000, IPDPS Workshops.

[2]  A. Needleman,et al.  A COMPARISON OF METHODS FOR CALCULATING ENERGY RELEASE RATES , 1985 .

[3]  Y. Murakami Stress Intensity Factors Handbook , 2006 .

[4]  H. G. deLorenzi,et al.  3-D elastic-plastic fracture mechanics with ADINA , 1981 .

[5]  Brian Moran,et al.  Energy release rate along a three-dimensional crack front in a thermally stressed body , 1986, International Journal of Fracture.

[6]  Paul A. Wawrzynek,et al.  Simulating Fatigue Crack Growth in Spiral Bevel Gears , 2001 .

[7]  J. Rice A path-independent integral and the approximate analysis of strain , 1968 .

[8]  S. Chan,et al.  On the Finite Element Method in Linear Fracture Mechanics , 1970 .

[9]  Brian Moran,et al.  Crack tip and associated domain integrals from momentum and energy balance , 1987 .

[10]  H. D. Bui,et al.  Associated path independent J-integrals for separating mixed modes , 1983 .

[11]  Robert H. Dodds,et al.  Numerical Evaluation of Domain and Contour Integrals for Nonlinear Fracture Mechanics: Formulation and Implementation Aspects , 1988 .

[12]  Satya N. Atluri,et al.  Comparison of different methods of evaluation of weight functions for 2-D mixed-mode fracture analyses , 1989 .

[13]  M. German,et al.  Crack extension modeling with singular quadratic isoparametric elements , 1976 .

[14]  A. Ingraffea,et al.  Stress‐intensity factor computation in three dimensions with quarter‐point elements , 1980 .

[15]  G. D. Bibel,et al.  Meshing of a Spiral Bevel Gearset With 3D Finite Element Analysis. , 1996 .

[16]  S. Vijayakar A combined surface integral and finite element solution for a three‐dimensional contact problem , 1991 .

[17]  S. Atluri,et al.  CALCULATION OF FRACTURE MECHANICS PARAMETERS FOR AN ARBITRARY THREE-DIMENSIONAL CRACK, BY THE ‘EQUIVALENT DOMAIN INTEGRAL’ METHOD , 1987 .

[18]  Yi Zhang,et al.  Local Synthesis and Tooth Contact Analysis of Face-Milled Spiral Bevel Gears , 1991 .

[19]  J. Barbera,et al.  Contact mechanics , 1999 .

[20]  B. Moran,et al.  A general treatment of crack tip contour integrals , 1987 .

[21]  G. D. Bibel,et al.  Comparison of Experimental and Analytical Tooth Bending Stress of Aerospace Spiral Bevel Gears , 1999 .

[22]  G. D. Bibel,et al.  Contact Stress Analysis of Spiral Bevel Gears Using Finite Element Analysis , 1995 .

[23]  R. Barsoum On the use of isoparametric finite elements in linear fracture mechanics , 1976 .