Search

In this paper, the contact problem between a cylinder and a half-space with a crystalline anisotropic behavior is solved. The model is based on semi-analytical methods to solve a three dimensional contact problem. A numerical technique based on a Voronoi tessellation is implemented using the Eshelby's equivalent inclusion method to account for the effect of the material microstructure on the contact pressure distribution and subsurface stresses. Fast Fourier Transforms (3D and 2D) are used to reduce the computation cost of the simulations. An application of this method to compute the scatter in fatigue life of rolling element bearings is also presented. Three different critical stresses are used in the Lundberg-Palmgren equation and results are compared in term of Weibull plot slope.

This paper presents a damage modeling approach to simulate the formation of butterfly wings in rolling contact fatigue coupled with a fast method for solving 3D contact problems with the presence of a spherical heterogeneity. The damage model is implemented in a semi-analytical model using the Eshelby’ s equivalent inclusion method in the contact solver. The proposed technique can easily include the effects of heterogeneous inclusions to the homogeneous solution in the contact algorithm. Contact pressure and subsurface stress field computation time is kept small due to the use of 3D and 2D Fast Fourier Transforms. Cuboidal inclusions with the same size as the discretization of the half-space are superimposed to represent the microstructural alterations during butterfly wings formation.

This paper presents a fast method of solving 3D contact problems when one of the mating bodies has an elastic-damageable behavior. The damage model is implemented in a semi-analytical model using Eshelby’ s equivalent inclusion method in the contact solver. The proposed technique can be seen as an enrichment technique for which the effect of heterogeneous inclusions is surimposed on the homogeneous solution in the contact algorithm. Contact pressure and subsurface stress field computation time is kept small due to a massive use of 3D and 2D Fast Fourier Transforms. Cuboidal inclusions with the same size as the discretization of the half-space and with the same elastic properties are surimposed. The damage model affects the elastic properties of the cuboidal inclusions. The emphasis is put on the effects of the fretting regimes on the contact pressure and damage evolution.