0
Journal of Tribology | Volume 131 | Issue 2 | Research Paper
Research Papers: Other (Seals, Manufacturing)

A Voronoi Finite Element Study of Fatigue Life Scatter in Rolling Contacts

[+] Author and Article Information
Behrooz Jalalahmadi

School of Mechanical Engineering, Purdue University, West Lafayette, IN 47907bjalalah@purdue.edu

Farshid Sadeghi

School of Mechanical Engineering, Purdue University, West Lafayette, IN 47907sadeghi@ecn.purdue.edu

J. Tribol 131(2), 022203 (Mar 04, 2009) (15 pages) doi:10.1115/1.3063818 History: Received June 19, 2008; Revised September 30, 2008; Published March 04, 2009
Article
References
Figures
Tables
Citing Articles

Microlevel material failure has been recognized as one of the main modes of failure for rolling contact fatigue (RCF) of bearing. Therefore, microlevel features of materials will be of significant importance to RCF investigation. At the microlevel, materials consist of randomly shaped and sized grains, which cannot be properly analyzed using the classical and commercially available finite element method. Hence, in this investigation, a Voronoi finite element method (VFEM) was developed to simulate the microstructure of bearing materials. The VFEM was then used to investigate the effects of microstructure randomness on rolling contact fatigue. Here two different types of randomness are considered: (i) randomness in the microstructure due to random shapes and sizes of the material grains, and (ii) the randomness in the material properties considering a normally (Gaussian) distributed elastic modulus. In this investigation, in order to determine the fatigue life, the model proposed by Raje (“A Numerical Model for Life Scatter in Rolling Element Bearings,” ASME J. Tribol., 130, pp. 011011-1–011011-10), which is based on the Lundberg–Palmgren theory (“Dynamic Capacity of Rolling Bearings,” Acta Polytech. Scand., Mech. Eng. Ser., 1(3), pp. 7–53), is used. This model relates fatigue life to a critical stress quantity and its corresponding depth, but instead of explicitly assuming a Weibull distribution of fatigue lives, the life distribution is obtained as an outcome of numerical simulations. We consider the maximum range of orthogonal shear stress and the maximum shear stress as the critical stress quantities. Forty domains are considered to study the effects of microstructure on the fatigue life of bearings. It is observed that the Weibull slope calculated for the obtained fatigue lives is in good agreement with previous experimental studies and analytical results. Introduction of inhomogeneous elastic modulus and initial flaws within the material domain increases the average critical stresses and decreases the Weibull slope.
Copyright © 2009 by American Society of Mechanical Engineers
Your Session has timed out. Please sign back in to continue.

References

1
Furuma, K., Shirota, S., and Hirakawa, K., 1975, “The Subsurface-Initiated and the Surface-Initiated Rolling Fatigue Life of Bearing Steels,” "Proceedings of the JSLE-ASLE International Conference on Lubrication", Tokyo, pp. 475–483.
2
Ioannides, E., and Harris, T. A., 1985, “A New Fatigue Life Model for Rolling Bearings,” ASME J. Tribol., 107 (3), pp. 367–378.
3
Zhou, R. S., 1993, “Surface Topography and Fatigue Life of Rolling Contact Bearings,” Tribol. Trans.  [CrossRef], 36 , pp. 329–340.
4
Zaretsky, E. V., 1994, “Design for Life, Plan for Death,” Mach. Des., 66 (15), pp. 55–59.
5
Cheng, W., and Cheng, H. S., 1995, “Semi-Analytical Modeling of Crack Initiation Dominant Contact Fatigue for Roller Bearings,” "Proceedings of the 1995 Joint ASME/STLE Tribology Conference", Orlando, FL.
6
Ringsberg, J. W., 2001, “Life Prediction of Rolling Contact Fatigue Crack Initiation,” Int. J. Fatigue  [CrossRef], 23 (7), pp. 575–586.
7
Guo, Y. B., and Barkey, M. E., 2004, “Modeling of Rolling Contact Fatigue for Hard Machined Components With Process-Induced Residual Stress,” Int. J. Fatigue  [CrossRef], 26 (6), pp. 605–613.
8
Li, Y. G., Kang, G. Z., Wang, C. G., Dou, P., and Wang, J., 2006, “Vertical Short-Crack Behavior and Its Application in Rolling Contact Fatigue,” Int. J. Fatigue, 28 (7), pp. 804–811.
9
Lee, H. Y., and Lee, D. Y., 2007, “The Effects of Indentation on Rolling Contact Fatigue Under Line Contact,” J. Korean Inst. Metals Mater., 45 (3), pp. 163–168.
10
Girodin, D., Dudragne, G., Courbon, J., and Vincent, A., 2006, “Statistical Analysis of Nonmetallic Inclusions for the Estimation of Rolling Contact Fatigue Range and Quality Control of Bearing Steel,” J. ASTM Int., 3 (7), pp. 85–100.
11
Andersson, J., 2005, “The Influence of Grain Size Variation on Metal Fatigue,” Int. J. Fatigue, 27 (8), pp. 847–852.
12
Miller, K. J., 1999, “A Historical Perspective of the Important Parameters of Metal Fatigue and Problems for the Next Century,” "Proceedings of the Seventh International Fatigue Congress, Fatigue ‘99", Beijing, pp. 15–39.
13
Bogdanski, S., and Trajer, M., 2005, “A Dimensionless Multi-Size Finite Element Model of a Rolling Contact Fatigue Crack,” Wear, 258 , pp. 1265–1272.
14
Melander, A., 1997, “A Finite Element Study of Short Cracks With Different Inclusion Types Under Rolling Contact Fatigue Load,” Int. J. Fatigue  [CrossRef], 19 (1), pp. 13–24.
15
Liu, Y., Liu, L., and Mahadevan, S., 2007, “Analysis of Subsurface Crack Propagation Under Rolling Contact Loading in Railroad Wheels Using FEM,” Eng. Fract. Mech., 74 , pp. 2659–2674.
16
Guo, Y. B., and Barkey, M. E., 2004, “FE-Simulation of the Effects of Machining-Induced Residual Stress Profile on Rolling Contact of Hard Machined Components,” Int. J. Mech. Sci.  [CrossRef], 46 , pp. 371–388.
17
Sraml, M., Flasker, J., and Potrc, I., 2003, “Numerical Procedure for Predicting the Rolling Contact Fatigue Crack Initiation,” Int. J. Fatigue, 25 , pp. 585–595.
18
Ringsberg, J. W., Loo-Morrey, M., Josefson, B. L., Kapoor, A., and Beynon, J. H., 2000, “Prediction of Fatigue Crack Initiation for Rolling Contact Fatigue,” Int. J. Fatigue  [CrossRef], 22 , pp. 205–215.
19
Ghosh, S., and Mallett, R. L., 1994, “Voronoi Cell Finite Elements,” Compos. Struct., 50 , pp. 33–46.
20
Ghosh, S., and Moorthy, S., 1995, “Elastic-Plastic Analysis of Arbitrary Heterogeneous Materials With the Voronoi Finite Element Method,” Comput. Methods Appl. Mech. Eng., 121 , pp. 373–409.
21
Moorthy, S., and Ghosh, S., 1996, “A Model for Analysis of Arbitrary Composite and Porous Microstructures With Voronoi Cell Finite Elements,” Int. J. Numer. Methods Eng.  [CrossRef], 39 , pp. 2363–2398.
22
Grujicic, M., and Zhang, Y., 1998, “Determination of Effective Elastic Properties of Functionally Graded Materials Using Voronoi Cell Finite Element Method,” Mater. Sci. Eng., A, 251 , pp. 64–76.
23
Guo, R., Shi, H. J., and Yao, Z. H., 2003, “Modeling of Interfacial Debonding Crack in Particle Reinforced Composites Using Voronoi Cell Finite Element Method,” Comput. Mech., 32 , pp. 52–59.
24
Vena, P., and Gastaldi, D., 2005, “A Voronoi Cell Finite Element Model for the Indentation of Graded Ceramic Composites,” Composites, Part B, 36 , pp. 115–126.
25
Okabe, A., and Boots, B., 1992, "Spatial Tessellations: Concepts and Applications of Voronoi Diagrams", Wiley, New York.
26
Moller, J., 1994, "Lectures Notes on Random Voronoi Tessellations", Springer, Berlin.
27
Lundberg, G., and Palmgren, A., 1947, “Dynamic Capacity of Rolling Bearings,” Acta Polytech. Scand., Mech. Eng. Ser., 1 (3), pp. 7–53.
28
Raje, N., Sadeghi, F., Rateick, R. G., and Hoeprich, M. R., 2008, “A Numerical Model for Life Scatter in Rolling Element Bearings,” ASME J. Tribol.  [CrossRef], 130 , pp. 011011-1–011011-10.
29
Zaretsky, E. V., Parker, R. J., and Anderson, W. J., 1969, “A Study of Residual Stress Induced During Rolling,” ASME J. Lubr. Technol., 91 , pp. 314–319.
30
Chen, L., Chen, Q., and Shao, E., 1989, “Study on Initiation and Propagation Angles of Sub-Surface Cracks in GCr15 Bearing Steel Under Rolling Contact,” Wear  [CrossRef], 133 , pp. 205–218.
31
Chen, Q., Shao, E., Zhao, D., Gue, J., and Fan, Z., 1991, “Measurement of the Critical Size of Inclusions Initiating Contact Fatigue Cracks and Its Application in Bearing Steels,” Wear  [CrossRef], 147 , pp. 285–294.
32
Yoshioka, T., 1993, “Detection of Rolling Contact Sub-Surface Fatigue Cracks Using Acoustic Emissions Technique,” Lubr. Eng., 94 (4), pp. 303–308.
33
Harris, T. A., 2001, "Rolling Bearing Analysis", Wiley, New York, p. 696.

Figures

Grahic Jump Location
+(a) The values and (b) the critical depths of the maximum range of orthogonal shear stress (Δτxy) and the maximum shear stress (τmax) for the considered 40 domains with different microstructures
View Large  |  Save Figure  |  Download Slide (.ppt)
(a) The values and (b) the critical depths of the maximum range of orthogonal shear stress (Δτxy) and the maximum shear stress (τmax) for the considered 40 domains with different microstructures
Figure 10

(a) The values and (b) the critical depths of the maximum range of orthogonal shear stress (Δτxy) and the maximum shear stress (τmax) for the considered 40 domains with different microstructures

Grahic Jump Location
+Two-parameter Weibull plot for the obtained lives of the 40 considered domains using Δτxy and τmax as the critical stress
View Large  |  Save Figure  |  Download Slide (.ppt)
Two-parameter Weibull plot for the obtained lives of the 40 considered domains using Δτxy and τmax as the critical stress
Figure 12

Two-parameter Weibull plot for the obtained lives of the 40 considered domains using Δτxy and τmax as the critical stress

Grahic Jump Location
+The two different normal Gaussian distributions of the elastic modulus employed for Domain 1 with a 10% variation in the average value of 200 GPa
View Large  |  Save Figure  |  Download Slide (.ppt)
The two different normal Gaussian distributions of the elastic modulus employed for Domain 1 with a 10% variation in the average value of 200 GPa
Figure 13

The two different normal Gaussian distributions of the elastic modulus employed for Domain 1 with a 10% variation in the average value of 200 GPa

Grahic Jump Location
+(a) The values and (b) the critical depths of Δτxy and τmax obtained for the 40 domains using the homogeneous elastic modulus and its normal distribution
View Large  |  Save Figure  |  Download Slide (.ppt)
(a) The values and (b) the critical depths of Δτxy and τmax obtained for the 40 domains using the homogeneous elastic modulus and its normal distribution
Figure 14

(a) The values and (b) the critical depths of Δτxy and τmax obtained for the 40 domains using the homogeneous elastic modulus and its normal distribution

Grahic Jump Location
+The comparison of the Weibull life distribution for cases of the uniform elastic modulus and its normal distribution using (a) Δτxy as the critical stress and (b) τmax as the critical stress
View Large  |  Save Figure  |  Download Slide (.ppt)
The comparison of the Weibull life distribution for cases of the uniform elastic modulus and its normal distribution using (a) Δτxy as the critical stress and (b) τmax as the critical stress
Figure 15

The comparison of the Weibull life distribution for cases of the uniform elastic modulus and its normal distribution using (a) Δτxy as the critical stress and (b) τmax as the critical stress

Grahic Jump Location
+A sample simulated RVE and the five initial flaws introduced in it
View Large  |  Save Figure  |  Download Slide (.ppt)
A sample simulated RVE and the five initial flaws introduced in it
Figure 16

A sample simulated RVE and the five initial flaws introduced in it

Grahic Jump Location
+Variation in (a) the values and (b) the depths of Δτxy and τmax with the introduction of five initial flaws in the domains
View Large  |  Save Figure  |  Download Slide (.ppt)
Variation in (a) the values and (b) the depths of Δτxy and τmax with the introduction of five initial flaws in the domains
Figure 17

Variation in (a) the values and (b) the depths of Δτxy and τmax with the introduction of five initial flaws in the domains

Grahic Jump Location
+The comparison of the Weibull life distribution between the domains without initial flaw and the domains including five initial flaws using (a) Δτxy as the critical stress and (b) τmax as the critical stress
View Large  |  Save Figure  |  Download Slide (.ppt)
The comparison of the Weibull life distribution between the domains without initial flaw and the domains including five initial flaws using (a) Δτxy as the critical stress and (b) τmax as the critical stress
Figure 18

The comparison of the Weibull life distribution between the domains without initial flaw and the domains including five initial flaws using (a) Δτxy as the critical stress and (b) τmax as the critical stress

Grahic Jump Location
+The comparison of the Weibull life distribution between the domains having different numbers of initial flaws using (a) Δτxy as the critical stress and (b) τmax as the critical stress
View Large  |  Save Figure  |  Download Slide (.ppt)
The comparison of the Weibull life distribution between the domains having different numbers of initial flaws using (a) Δτxy as the critical stress and (b) τmax as the critical stress
Figure 19

The comparison of the Weibull life distribution between the domains having different numbers of initial flaws using (a) Δτxy as the critical stress and (b) τmax as the critical stress

Grahic Jump Location
+The variation in the Weibull slope in terms of numbers of initial flaws in the domains using (a) Δτxy as the critical stress and (b) τmax as the critical stress
View Large  |  Save Figure  |  Download Slide (.ppt)
The variation in the Weibull slope in terms of numbers of initial flaws in the domains using (a) Δτxy as the critical stress and (b) τmax as the critical stress
Figure 20

The variation in the Weibull slope in terms of numbers of initial flaws in the domains using (a) Δτxy as the critical stress and (b) τmax as the critical stress

Grahic Jump Location
+Two-parameter Weibull plot for the critical depths of the 40 considered domains using Δτxy and τmax as the critical stress
View Large  |  Save Figure  |  Download Slide (.ppt)
Two-parameter Weibull plot for the critical depths of the 40 considered domains using Δτxy and τmax as the critical stress
Figure 11

Two-parameter Weibull plot for the critical depths of the 40 considered domains using Δτxy and τmax as the critical stress

Grahic Jump Location
+(a) Discretizing a domain with the Voronoi elements, (b) dividing a Voronoi element into finer triangular elements, and (c) linear triangular element
View Large  |  Save Figure  |  Download Slide (.ppt)
(a) Discretizing a domain with the Voronoi elements, (b) dividing a Voronoi element into finer triangular elements, and (c) linear triangular element
Figure 1

(a) Discretizing a domain with the Voronoi elements, (b) dividing a Voronoi element into finer triangular elements, and (c) linear triangular element

Grahic Jump Location
+The two approaches for storing the global stiffness matrix (K): (a) a single sparse matrix and (b) two much smaller condensed matrices
View Large  |  Save Figure  |  Download Slide (.ppt)
The two approaches for storing the global stiffness matrix (K): (a) a single sparse matrix and (b) two much smaller condensed matrices
Figure 2

The two approaches for storing the global stiffness matrix (K): (a) a single sparse matrix and (b) two much smaller condensed matrices

Grahic Jump Location
+The theoretical internal stress distribution below the surface for a cylinder in contact with a semi-infinite domain
View Large  |  Save Figure  |  Download Slide (.ppt)
The theoretical internal stress distribution below the surface for a cylinder in contact with a semi-infinite domain
Figure 3

The theoretical internal stress distribution below the surface for a cylinder in contact with a semi-infinite domain

Grahic Jump Location
+(a) The dimensions and grain sizes of the domain simulating the semi-infinite domain forming the bearing line contact, and (b) the zoomed view
View Large  |  Save Figure  |  Download Slide (.ppt)
(a) The dimensions and grain sizes of the domain simulating the semi-infinite domain forming the bearing line contact, and (b) the zoomed view
Figure 4

(a) The dimensions and grain sizes of the domain simulating the semi-infinite domain forming the bearing line contact, and (b) the zoomed view

Grahic Jump Location
+Statistical distribution of elements in one of the generated domains: (a) distribution of number of sides of Voronoi elements in the domain and (b) distribution of the elements’ area
View Large  |  Save Figure  |  Download Slide (.ppt)
Statistical distribution of elements in one of the generated domains: (a) distribution of number of sides of Voronoi elements in the domain and (b) distribution of the elements’ area
Figure 5

Statistical distribution of elements in one of the generated domains: (a) distribution of number of sides of Voronoi elements in the domain and (b) distribution of the elements’ area

Grahic Jump Location
+The stress contours obtained for a randomly generated domain under the static Hertzian loading: (a) τxy and (b) σy
View Large  |  Save Figure  |  Download Slide (.ppt)
The stress contours obtained for a randomly generated domain under the static Hertzian loading: (a) τxy and (b) σy
Figure 6

The stress contours obtained for a randomly generated domain under the static Hertzian loading: (a) τxy and (b) σy

Grahic Jump Location
+Variations in normal stress in the y-direction (σy) and orthogonal shear stress (τxy) along the depth for four randomly generated domains
View Large  |  Save Figure  |  Download Slide (.ppt)
Variations in normal stress in the y-direction (σy) and orthogonal shear stress (τxy) along the depth for four randomly generated domains
Figure 7

Variations in normal stress in the y-direction (σy) and orthogonal shear stress (τxy) along the depth for four randomly generated domains

Grahic Jump Location
+The stress contours obtained from ANSYS using the square elements for a domain under the static Hertzian loading
View Large  |  Save Figure  |  Download Slide (.ppt)
The stress contours obtained from ANSYS using the square elements for a domain under the static Hertzian loading
Figure 8

The stress contours obtained from ANSYS using the square elements for a domain under the static Hertzian loading

Grahic Jump Location
+The stress contours obtained from ANSYS using the linear triangular elements for a domain under the static Hertzian loading
View Large  |  Save Figure  |  Download Slide (.ppt)
The stress contours obtained from ANSYS using the linear triangular elements for a domain under the static Hertzian loading
Figure 9

The stress contours obtained from ANSYS using the linear triangular elements for a domain under the static Hertzian loading

Tables

Errata

Discussions

Some tools below are only available to our subscribers or users with an online account.

Related Content

Customize your page view by dragging and repositioning the boxes below.

Related Journal Articles
Related eBook Content
Subsection NB—Class 1 Components
Companion Guide to the ASME Boiler & Pressure Vessel Code, Volume 1, Second Edition> Chapter 6
Data Tabulations
Structural Shear Joints> Chapter 12
Fatigue Analysis in the Connecting Rod of MF285 Tractor by Finite Element Method
International Conference on Advanced Computer Theory and Engineering, 4th (ICACTE 2011)
Analysis of Components in VIII-2
Guidebook for the Design of ASME Section VIII Pressure Vessels, Third Edition> Chapter 8
Openings
Guidebook for the Design of ASME Section VIII Pressure Vessels> Chapter 5
Topic Collections

Sorry! You do not have access to this content. For assistance or to subscribe, please contact us:

  • TELEPHONE: 1-800-843-2763 (Toll-free in the USA)
  • EMAIL: asmedigitalcollection@asme.org
Sign In