Contact analysis in the presence of an ellipsoidal inhomogeneity within a half space

Abstract Many materials contain inhomogeneities or inclusions that may greatly affect their mechanical properties. Such inhomogeneities are for example encountered in the case of composite materials or materials containing precipitates. This paper presents an analysis of contact pressure and subsurface stress field for contact problems in the presence of anisotropic elastic inhomogeneities of ellipsoidal shape. Accounting for any orientation and material properties of the inhomogeneities are the major novelties of this work. The semi-analytical method proposed to solve the contact problem is based on Eshelby’s formalism and uses 2D and 3D Fast Fourier Transforms to speed up the computation. The time and memory necessary are greatly reduced in comparison with the classical finite element method. The model can be seen as an enrichment technique where the enrichment fields from the heterogeneous solution are superimposed to the homogeneous problem. The definition of complex geometries made by combination of inclusions can easily be achieved. A parametric analysis on the effect of elastic properties and geometrical features of the inhomogeneity (size, depth and orientation) is proposed. The model allows to obtain the contact pressure distribution – disturbed by the presence of inhomogeneities – as well as subsurface and matrix/inhomogeneity interface stresses. It is shown that the presence of an inclusion below the contact surface affects significantly the contact pressure and subsurfaces stress distributions when located at a depth lower than 0.7 times the contact radius. The anisotropy directions and material data are also key elements that strongly affect the elastic contact solution. In the case of normal contact between a spherical indenter and an elastic half space containing a single inhomogeneity whose center is located straight below the contact center, the normal stress at the inhomogeneity/matrix interface is mostly compressive. Finally when the axes of the ellipsoidal inclusion do not coincide with the contact problem axes, the pressure distribution is not symmetrical.

[1]  J. D. Eshelby.  The determination of the elastic field of an ellipsoidal inclusion, and related problems , 1957, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.

[2]  Daniel Nelias,et al.  Multiscale computation of fretting wear at the blade/disk interface , 2010 .

[3]  Leon M Keer,et al.  A review of recent works on inclusions , 2013 .

[4]  Krich Sawamiphakdi,et al.  Interaction of multiple inhomogeneous inclusions beneath a surface , 2012 .

[5]  A. Love.  A treatise on the mathematical theory of elasticity , 1892 .

[6]  Daniel Nelias,et al.  A Comprehensive Method to Predict Wear and to Define the Optimum Geometry of Fretting Surfaces , 2006 .

[7]  Ramón Abascal,et al.  Analysis of FRP composites under frictional contact conditions , 2013 .

[8]  Leon M Keer,et al.  A fast method for solving three-dimensional arbitrarily shaped inclusions in a half space , 2009 .

[9]  A. P. Voskamp.  Material Response to Rolling Contact Loading , 1985 .

[10]  I. S. Gradshteyn,et al.  Table of Integrals, Series, and Products , 1976 .

[11]  G. M.,et al.  A Treatise on the Mathematical Theory of Elasticity , 1906, Nature.

[12]  Daniel Nelias,et al.  Contact analysis in presence of spherical inhomogeneities within a half-space , 2010 .

[13]  Toshio Mura,et al.  Two-Ellipsoidal Inhomogeneities by the Equivalent Inclusion Method , 1975 .

[14]  H. Y. Yu,et al.  Elastic inclusions and inhomogeneities in transversely isotropic solids , 1994, Proceedings of the Royal Society of London. Series A: Mathematical and Physical Sciences.

[15]  D. Nélias,et al.  Contact Analyses for Anisotropic Half Space: Effect of the Anisotropy on the Pressure Distribution and Contact Area , 2012 .

[16]  C. Kuo.  Contact stress analysis of an elastic half-plane containing multiple inclusions , 2008 .

[17]  D. Nélias,et al.  Contact Pressure and Residual Strain in 3D Elasto-Plastic Rolling Contact for a Circular or Elliptical Point Contact , 2011 .

[18]  L. Gallego,et al.  Modeling of Fretting Wear Under Gross Slip and Partial Slip Conditions , 2007 .

[19]  Y. P. Chiu,et al.  On the Stress Field and Surface Deformation in a Half Space With a Cuboidal Zone in Which Initial Strains Are Uniform , 1978 .

[20]  J. D. Eshelby,et al.  The elastic field outside an ellipsoidal inclusion , 1959, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.

[21]  Daniel Nelias,et al.  On the Tangential Displacement of a Surface Point Due to a Cuboid of Uniform Plastic Strain in a Half-Space , 2010 .

[22]  V. Vignal,et al.  Local analysis of the mechanical behaviour of inclusions-containing stainless steels under straining conditions , 2003 .

[23]  Qian Wang,et al.  Elastic Fields due to Eigenstrains in a Half-Space , 2005 .

[24]  D. Nélias,et al.  Contact analyses for anisotropic half-space coated with an anisotropic layer: Effect of the anisotropy on the pressure distribution and contact area , 2013 .

[25]  Raymond D. Mindlin,et al.  Thermoelastic Stress in the Semi-Infinite Solid , 1950 .

[26]  J. Willis.  ANISOTROPIC ELASTIC INCLUSION PROBLEMS , 1964 .

[27]  Toshio Mura,et al.  The Elastic Field in a Half Space Due to Ellipsoidal Inclusions With Uniform Dilatational Eigenstrains , 1979 .

[28]  Toshio Mura,et al.  The Elastic Field Outside an Ellipsoidal Inclusion , 1977 .

[29]  Y. P. Chiu,et al.  On the Stress Field Due to Initial Strains in a Cuboid Surrounded by an Infinite Elastic Space , 1977 .

[30]  Gilles Dudragne,et al.  Influence of inclusion pairs, clusters and stringers on the lower bound of the endurance limit of bearing steels , 2003 .

[31]  Daniel Nelias,et al.  A fast and efficient contact algorithm for fretting problems applied to fretting modes I, II and III , 2010 .

[32]  Leon M Keer,et al.  Semi‐analytic solution for multiple interacting three‐dimensional inhomogeneous inclusions of arbitrary shape in an infinite space , 2011 .

[33]  L. Walpole,et al.  The elastic field of an inclusion in an anisotropic medium , 1967, Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences.

[34]  Leon M Keer,et al.  Numerical Simulation for Three Dimensional Elastic-Plastic Contact with Hardening Behavior , 2005 .

[35]  Q. Wang,et al.  A versatile method of discrete convolution and FFT (DC-FFT) for contact analyses , 2000 .

[36]  J. Leroux,et al.  Stick-slip analysis of a circular point contact between a rigid sphere and a flat unidirectional composite with cylindrical fibers , 2011 .

[37]  Wei Chen,et al.  Fast Fourier Transform Based Numerical Methods for Elasto-Plastic Contacts of Nominally Flat Surfaces , 2008 .

[38]  K. Aderogba,et al.  On eigenstresses in a semi-infinite solid , 1976, Mathematical Proceedings of the Cambridge Philosophical Society.

[39]  R. Fougères,et al.  Role of inclusions, surface roughness and operating conditions on rolling contact fatigue , 1999 .

[40]  Daniel Nelias,et al.  Development of a Three-Dimensional Semi-Analytical Elastic-Plastic Contact Code , 2002 .

[41]  R. Asaro,et al.  The non-uniform transformation strain problem for an anisotropic ellipsoidal inclusion , 1975 .

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

[43]  Toshio Mura,et al.  Micromechanics of defects in solids , 1982 .

[44]  Daniel Nelias,et al.  Contact Analyses for Bodies With Frictional Heating and Plastic Behavior , 2005 .

[45]  Wei Chen,et al.  Multiple 3D inhomogeneous inclusions in a half space under contact loading , 2011 .

[46]  Leon M Keer,et al.  Three-Dimensional Repeated Elasto-Plastic Point Contacts, Rolling, and Sliding , 2008 .

[47]  Philippe Gilles,et al.  Modelling of multiple impacts for the prediction of distortions and residual stresses induced by ultrasonic shot peening (USP) , 2012 .

[48]  Daniel Nelias,et al.  Rolling of an Elastic Ellipsoid Upon an Elastic-Plastic Flat , 2007 .

[49]  Elena Kabo,et al.  Fatigue initiation in railway wheels — a numerical study of the influence of defects , 2002 .

[50]  C. H. Kuo,et al.  Stress disturbances caused by the inhomogeneity in an elastic half-space subjected to contact loading , 2007 .

[51]  Hui Wang,et al.  A Simplified Numerical Elastic-Plastic Contact Model for Rough Surfaces , 2009 .

[52]  Daniel Nelias,et al.  On the Effect of Isotropic Hardening on the Coefficient of Restitution for Single or Repeated Impacts Using a Semi-Analytical Method , 2011 .

[53]  A. Mortensen,et al.  EXPERIMENTAL INVESTIGATION OF STRESS AND STRAIN FIELDS IN A DUCTILE MATRIX SURROUNDING AN ELASTIC INCLUSION , 2000 .

[54]  Daniel Nelias,et al.  A Three-Dimensional Semianalytical Model for Elastic-Plastic Sliding Contacts , 2007 .

[55]  Elena Kabo,et al.  Material defects in rolling contact fatigue of railway wheels—the influence of defect size , 2005 .

[56]  Toshio Mura,et al.  The Elastic Inclusion With a Sliding Interface , 1984 .