Three-Dimensional Finite Element Elastic–Plastic Model for Subsurface Initiated Spalling in Rolling Contacts

2013 ◽  
Vol 136 (1) ◽  
Author(s):  
John A. R. Bomidi ◽  
Farshid Sadeghi

In this investigation, a three-dimensional (3D) finite element (FE) model was developed to study subsurface initiated spalling observed in rolling line contact of tribo components such as bearings. An elastic–kinematic hardening–plastic material model is employed to capture the material behavior of bearing steel and is coupled with the continuum damage mechanics (CDM) approach to capture the material degradation due to fatigue. The fatigue damage model employs both stress and accumulated plastic strain based damage evolution laws for fatigue failure initiation and propagation. Failure is modeled by mesh partitioning along unstructured, nonplanar, intergranular paths of the microstructure topology represented by randomly generated Voronoi tessellations. The elastic–plastic model coupled with CDM was used to predict both ratcheting behavior and fatigue damage in heavily loaded contacts. Fatigue damage induced due to the accumulated plastic strains around broken intergranular joints drive the majority of the crack propagation stage, resulting in a lower percentage of life spent in propagation. The 3D FE model was used to determine fatigue life at different contact pressures ranging from 2 to 4.5 GPa for 33 different randomly generated microstructure topology models. The effect of change in contact pressure due to subsurface damage and plastic strain accumulation was also captured by explicitly modeling the rolling contact geometry and the results were compared to those generated assuming a Hertzian pressure profile. The spall shape, fatigue lives, and their dispersion characterized by Weibull slopes obtained from the model correlate well with the previously published experimental results.

2019 ◽  
Vol 141 (4) ◽  
Author(s):  
Zamzam Golmohammadi ◽  
Farshid Sadeghi

A coupled multibody elastic–plastic finite element (FE) model was developed to investigate the effects of surface defects, such as dents on rolling contact fatigue (RCF). The coupled Voronoi FE model was used to determine the contact pressure acting over the surface defect, internal stresses, damage, etc. In order to determine the shape of a dent and material pile up during the over rolling process, a rigid indenter was pressed against an elastic plastic semi-infinite domain. Continuum damage mechanics (CDM) was used to account for material degradation during RCF. Using CDM, spall initiation and propagation in a line contact was modeled and investigated. A parametric study using the model was performed to examine the effects of dent sharpness, pile up ratio, and applied load on the spall formation and fatigue life. The spall patterns were found to be consistent with experimental observations from the open literature. Moreover, the results demonstrated that the dent shape and sharpness had a significant effect on pressure and thus fatigue life. Higher dent sharpness ratios significantly reduced the fatigue life.


1990 ◽  
Vol 112 (3) ◽  
pp. 287-291 ◽  
Author(s):  
F. A. Kolkailah ◽  
A. J. McPhate

In this paper, results from an elastic-plastic finite-element model incorporating the Bodner-Partom model of nonlinear time-dependent material behavior are presented. The parameters in the constitutive model are computed from a leastsquare fit to experimental data obtained from uniaxial stress-strain and creep tests at 650°C. The finite element model of a double-notched specimen is employed to determine the value of the elastic-plastic strain and is compared to experimental data. The constitutive model parameters evaluated in this paper are found to be in good agreement with those obtained by the other investigators. However, the parameters determined by the numerical technique tend to give response that agree with the data better than do graphically determined parameters previously used. The calculated elastic-plastic strain from the model agreed well with the experimental strain.


2014 ◽  
Vol 658 ◽  
pp. 207-212
Author(s):  
Gabriel Popescu

An analytical three-dimensional elastic-plastic over-rolling solution is used to evaluate the plastic strains and residual stresses. Central to this plastic contact formulation is the incremental approach to deal with non-linear material behavior. The Prandtl-Reuss constitutive equations in conjunction with Huber-Mises-Hencky yield criterion and Ramberg-Osgood strain-hardening relationships are applied to describe the plastic behavior of common hardened bearing steel. The model was extended to include the tangential force in the rolling direction, assumed to be proportional to the hertzian contact pressure. Comparisons of three-dimensional pure rolling and rolling/sliding contact results were provided to elucidate the differences in residual stresses and residual profiles in case of kinematic and work-hardening materials.


Author(s):  
Ivan Karin ◽  
Nils Lommatzsch ◽  
Klaus Lipp ◽  
Volker Landersheim ◽  
Holger Hanselka ◽  
...  

Within the collaborative research centre 666 “Integral Sheet Metal Design with Higher Order Bifurcations” the innovative manufacturing technologies linear flow-splitting and linear bend-splitting are researched that allow the continuous production of multi-chambered steel profiles in integral style. The massive forming processes create an ultra-fine grained microstructure in the forming area that is characterized by an increased hardness and lower surface roughness compared to as received material. These properties predestine the technology to be used in the production of linear guides. Additionally, the multi-chambered structure of the linear flow-split and -bend components can be used for function integration. To design and evaluate linear guides that use the whole technological potential, the research is focused on a macroscopic and a microscopic point of view. The macroscopic approach is targeting the development of linear flow-split linear guides with integrated functions to provide additional performance values to the established machine parts. Continuously produced guidance systems with innovative functionality can be introduced to a new market with the technology push approach. Preliminary designs of linear flow-split guidance systems and integrated functions are promising. Therefore, an approach to develop new functions for linear flow-split linear guides basing on calculation models and property networks is shown [1]. With this approach, optimized solutions can be created and possible design modifications can be derived. In this contribution, the development and integration of a clamping function for decelerating the slide is presented. Calculation models for analyzing the functionality are presented and validated by finite element models and experiments. The microscopic examination of the profiles aims to investigate the material behavior, particularly of the formed areas. Beside the conventional mechanical and fatigue properties of linear flow-split material ZStE500 [2], the present work focuses on the rolling contact fatigue. This is necessary to evaluate linear flow-split components regarding their eligibility with regard to the rolling contact fatigue behaviour. The Hertz theory for rolling contact fatigue is only valid for homogeneous materials [3]. The flow-split material ZStE500 shows a non-homogeneous behaviour and has to be analyzed with the Finite Element Method in order to determine stresses and strains. In comparison to simulation results with unformed and therefore homogeneous material, the effect of linear flow-split surfaces on the rolling contact behavior is demonstrated. Based on these results, it is possible to start experimental investigations on rolling contact fatigue of linear flow-split components to validate the FE model and determine the performance of linear flow-split flanges for rolling contact fatigue.


1991 ◽  
Vol 113 (3) ◽  
pp. 434-441 ◽  
Author(s):  
S. M. Kulkarni ◽  
G. T. Hahn ◽  
C. A. Rubin ◽  
V. Bhargava

This paper presents an elasto-plastic analysis of the repeated, frictionless, three-dimensional rolling contact similar to the ones produced by the rail-wheel geometry. This paper treats an elliptical contact rolling across a semi-infinite half space. The contact shape and loading: semi-major axis (in the rolling direction), w1 = 8 mm, and semi-minor axis, w2 = 5.88 mm, reflect standard rail and wheel curvatures and a wheel load of 149 KN (33,000 lb). A three-dimensional, elasto-plastic finite element model, developed earlier, is employed together with the elastic-linear-kinematic-hardening-plastic (ELKP) idealization of the cyclic plastic behaviour of a material similar to rail and wheel steels. The calculations present the displacements, the stress-strain distributions, stress-plastic strain histories and the plastic strain ranges in the half-space. The cyclic plasticity approaches a steady state after one contact with further contacts producing open but fully reversed stress-strain hysteresis loops, i.e., plastic shakedown.


2016 ◽  
Vol 138 (3) ◽  
Author(s):  
Tao He ◽  
Dong Zhu ◽  
Jiaxu Wang

Surface plastic deformation due to contact (lubricated or dry) widely exists in many mechanical components, as subsurface stress caused by high-pressure concentrated in the contact zone often exceeds the material yielding limit, and the plastic strain accumulates when the load is increased and/or repeatedly applied to the surface in a rolling contact. However, previous plasto-elastohydrodynamic lubrication (PEHL) studies were mainly for the preliminary case of having a rigid ball (or roller) rotating on a stationary elastic–plastic flat with a fixed contact center, for which the numerical simulation is relatively simple. This paper presents an efficient method for simulating PEHL in a rolling contact. The von Mises yield criteria are used for determining the plastic zone, and the total computation domain is discretized into a number of cuboidal elements underneath the contacting surface, each one is considered as a cuboid with uniform plastic strain inside. The residual stress and surface plastic deformation resulted from the plastic strain can be solved as a half-space eigenstrain–eigenstress problem. A combination of three-dimensional (3D) and two-dimensional (2D) discrete convolution and fast Fourier transform (DC-FFT) techniques is used for accelerating the computation. It is observed that if a rigid ball rolls on an elastic–plastic surface, the characteristics of PEHL lubricant film thickness and pressure distribution are different from those of PEHL in the preliminary cases previously investigated. It is also found that with the increase of rolling cycles, the increment of plastic strain accumulation gradually approaches a stable value or drops down to zero, determined by the applied load and the material hardening properties, eventually causing a groove along the rolling direction. Simulation results for different material hardening properties are also compared to reveal the effect of body materials on the PEHL behaviors.


2002 ◽  
Vol 124 (4) ◽  
pp. 699-708 ◽  
Author(s):  
Yanyao Jiang ◽  
Biqiang Xu ◽  
Huseyin Sehitoglu

Three-dimensional elastic-plastic rolling contact stress analysis was conducted incorporating elastic and plastic shakedown concepts. The Hertzian distribution was assumed for the normal surface contact load over a circular contact area. The tangential forces in both the rolling and lateral directions were considered and were assumed to be proportional to the Hertzian pressure. The elastic and plastic shakedown limits obtained for the three-dimensional contact problem revealed the role of both longitudinal and lateral shear traction on the shakedown results. An advanced cyclic plasticity model was implemented into a finite element code via the material subroutine. Finite element simulations were conducted to study the influences of the tangential surface forces in the two shear directions on residual stresses and residual strains. For all the cases simulated, the p0/k ratio (p0 is the maximum Hertzian pressure and k is the yield stress in shear) was 6.0. The Qx/P ratio, where Qx is the total tangential force on the contact surface in the rolling direction and P is the total normal surface pressure, ranged from 0 to 0.6. The Qy/P ratio (Qy is the total tangential force in the lateral direction) was either zero or 0.25. Residual stresses increase with increasing rolling passes but tend to stabilize. Residual strains also increase but the increase in residual strain per rolling pass (ratchetting rate) decays with rolling cycles. Residual stress levels can be as high as 2k when the Qx/P ratio is 0.6. Local accumulated shear strains can exceed 20 times the yield strain in shear after six rolling passes under extreme conditions. Comparisons of the two-dimensional and three-dimensional rolling contact results were provided to elucidate the differences in residual stresses and ratchetting strain predictions.


Author(s):  
Xin Zhao ◽  
Zili Li ◽  
Rolf Dollevoet

A full-scale 3-D dynamic finite element (FE) model is created to solve the elastic-plastic wheel-rail rolling contact over rough surfaces under different friction forces. Both normal and tangential loads are applied properly. A bi-linear plastic material model is introduced and the real wheel and rail head geometries are simulated. The rolling of a drive wheel with full friction exploitation over a rough contact surface is analyzed in this paper. The stress distributions at the zone with rough surface are derived. From the results at a selected instant, it is found that roughness significantly increases the stress level of the surface layer. Furthermore, plasticity can greatly reduce stress peaks and change stress distributions. The maximum shear stress distribution at the rough surface is also analyzed to assess effects of roughness on fatigue crack initiation.


2015 ◽  
Vol 137 (4) ◽  
Author(s):  
Xiangyun Deng ◽  
Zhiwei Qian ◽  
Rolf Dollevoet

Spin in frictional rolling contact can cause significant stress, which is the key to understanding and predicting the wear and fatigue behavior of contact components, such as wheels, rails, and rolling bearings. The lateral creep force arising from spin influences the kinematics of a wheelset and thus of vehicles. The solution that is currently employed in the field of elasticity and continuum statics was developed by Kalker and uses a boundary element method (BEM). In this paper, a new approach based on Lagrangian explicit finite element (FE) analysis is employed. This approach is able to consider arbitrary geometric profiles of rails and wheels, complex material behavior and dynamic effects, and some other factors. The new approach is demonstrated using a three-dimensional (3D) model of a wheel with a coned profile rolling along a quarter cylinder and can be easily adapted to apply to wheels and rails of arbitrary profiles. The 3D FE model is configured with elastic material properties and is used to obtain both normal and tangential solutions. The results are compared with those of the Hertz theory and the Kalker's model. The 3D FE model is then configured with elastoplastic material properties to study the spin-rolling contact with plasticity. The continuum dynamics phenomenon is captured by the FE model, which enhances the ability of the model to mimic reality. This improvement considerably extends the applicability of the FE model. The model can be applied to fatigue and wear analyses at gauge corners or rails as well as to deep groove bearings, where a large geometrical spin is present and plastic deformation may be of importance.


1995 ◽  
Vol 117 (2) ◽  
pp. 234-243 ◽  
Author(s):  
Maria M.-H. Yu ◽  
Brian Moran ◽  
Leon M. Keer

A novel treatment of a direct procedure for elastic-plastic analysis and shakedown is presented and its application to problems in three-dimensional rolling contact with or without case-hardened layers is demonstrated. The direct approach consists of an operator split technique, which transforms the elastic-plastic problem into a purely elastic problem and a residual problem with prescribed eigenstrains. These eigenstrains are determined using an incremental projection method based on the purely elastic solution and a special representation of the yield condition for a linear-kinematic hardening material. The three-dimensional residual problem is then further split into a plane problem and an anti-plane problem which are readily solved using the finite element method. A significant advantage of the present analysis over the alternative approach of simulating repeated rolling until shakedown occurs is that in the present analysis, the final shakedown solution is obtained directly by solving three elasticity problems. Results are compared with full elastic-plastic finite element calculations available from the literature and good agreement is observed. The effects of surface hardened layers on the distributions of residual stress and displacement are investigated for both two- and three-dimensional contact. The direct approach is shown to be a straightforward and efficient method for obtaining the steady state solution in the analysis of three-dimensional problems in rolling and/or sliding contact.


Sign in / Sign up

Export Citation Format

Share Document