Numerical methods for solving the equivalent inclusion equation in semi-analytical models

Equivalent inclusion method is the basis for semi-analytical models in tackling inhomogeneity problems. Equivalent eigenstrains are obtained by solving the consistency equation system of the equivalent inclusion method and then stress disturbances caused by inhomogeneities are determined. The equivalent inclusion method equation system can only be solved numerically, but the current fixed-point iteration method may not be able to achieve deep convergence when the Young's modulus of inhomogeneity is lower than that of the matrix material. The most significant innovation of this paper is to reveal the non-convergence mechanism of the current method. Considering the limitation, the Jacobian-free Newton Krylov algorithm is selected to solve the equivalent inclusion method equation. Results indicate that the new algorithm has significant advantages of computing accuracy and efficiency compared with the classic method.

[1]  Qin Fan,et al.  Partial slip contact modeling of heterogeneous elasto-plastic materials , 2016 .

[2]  Bin Zhao,et al.  Contact behaviors of multilayered structures with interfacial cracks , 2020 .

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

[4]  Microscopic investigation of subsurface initiated damage of wind turbine gearbox bearings , 2018, Journal of Physics: Conference Series.

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

[6]  Qian Wang,et al.  Efficient numerical method with a dual-grid scheme for contact of inhomogeneous materials and its applications , 2018 .

[7]  Jiaxu Wang,et al.  A mesh differential refinement scheme for solving elastic fields of half-space inclusion problems , 2016 .

[8]  Pin Lu,et al.  Effective moduli of nanoparticle reinforced composites considering interphase effect by extended double-inclusion model – Theory and explicit expressions , 2013 .

[9]  A thermoelastic contact model between a sliding ball and a stationary half space distributed with spherical inhomogeneities , 2019, Tribology International.

[10]  I. Dzioba,et al.  Fracture Mechanisms of S355 Steel—Experimental Research, FEM Simulation and SEM Observation , 2019, Materials.

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

[12]  F. Sadeghi,et al.  A coupled damage model and a semi-analytical contact solver to simulate butterfly wing formation around nonmetallic inclusions , 2019, International Journal of Fatigue.

[13]  Leon M Keer,et al.  Analytical solution for elastic fields caused by eigenstrains in a half-space and numerical implementation based on FFT , 2012 .

[14]  A. Cruz-Ramírez,et al.  The effect of non-metallic inclusions on the mechanical properties of 32 CDV 13 steel and their mechanical stress analysis by numerical simulation , 2018 .

[15]  Yousef Saad,et al.  Hybrid Krylov Methods for Nonlinear Systems of Equations , 1990, SIAM J. Sci. Comput..

[16]  Sebastian Lipa,et al.  NUMERICAL ANALYSIS OF THE EFFECT OF NONMETALLIC INCLUSIONS LOCATION IN CARBURIZED LAYER OF 16MnCr5 STEEL ON INITIATION OF FATIGUE CRACKING , 2018, Advances in Science and Technology Research Journal.

[17]  Daniel Nelias,et al.  Fully Coupled Resolution of Heterogeneous Elastic–Plastic Contact Problem , 2016 .

[18]  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.

[19]  Jiaxu Wang,et al.  Numerical EIM with 3D FFT for the contact with a smooth or rough surface involving complicated and distributed inhomogeneities , 2016 .

[20]  H. A. Al-Tameemi,et al.  Initiation of sub-surface micro-cracks and white etching areas from debonding at non-metallic inclusions in wind turbine gearbox bearing , 2018, Wear.

[21]  R. S. Dwyer-Joyce,et al.  Threshold Maps for Inclusion-Initiated Micro-Cracks and White Etching Areas in Bearing Steel: The Role of Impact Loading and Surface Sliding , 2018, Tribology Letters.

[22]  L. Keer,et al.  Elasto-plastic contact of materials containing double-layered inhomogeneities , 2017 .

[23]  T. Chan,et al.  Nonlinearly Preconditioned Krylov Subspace Methods for Discrete Newton Algorithms , 1984 .

[24]  D. Keyes,et al.  Jacobian-free Newton-Krylov methods: a survey of approaches and applications , 2004 .

[25]  Leon M Keer,et al.  Novel Model for Partial-Slip Contact Involving a Material With Inhomogeneity , 2013 .

[26]  Yousef Saad,et al.  A Flexible Inner-Outer Preconditioned GMRES Algorithm , 1993, SIAM J. Sci. Comput..

[27]  Yan Li,et al.  Layered structures with rough surfaces and interfaces at contact loading , 2020 .

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

[29]  Hui Long,et al.  Characterisation of white etching crack damage in wind turbine gearbox bearings , 2015 .

[30]  Leon M Keer,et al.  Numerical Implementation of the Equivalent Inclusion Method for 2D Arbitrarily Shaped Inhomogeneities , 2015 .

[31]  Dong Zhu,et al.  Elastohydrodynamic Lubrication of Inhomogeneous Materials Using the Equivalent Inclusion Method , 2014 .

[32]  P. Bagot,et al.  Solute redistribution in the nanocrystalline structure formed in bearing steels , 2013 .

[33]  Wei Liu,et al.  Control of MnS Inclusions in High- and Low-Sulfur Steel by Tellurium Treatment , 2019, Materials.

[34]  M. Jin,et al.  High temperature stress–strain curves of MnS and their applications in finite element simulation , 2018, Transactions of Nonferrous Metals Society of China.

[35]  Leon M Keer,et al.  A Numerical Approach for Analyzing Three-Dimensional Steady-State Rolling Contact Including Creep Using a Fast Semi-Analytical Method , 2012 .

[36]  Hossein M. Shodja,et al.  Elastic Fields in Double Inhomogeneity by the Equivalent Inclusion Method , 2001 .