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.
KangJHHosseinkhaniBWilliamsCA, et al.Solute redistribution in the nanocrystalline structure formed in bearing steels. Scr Mater2013; 69: 630–633. doi:10.1016/j.scriptamat.2013.07.017
2.
BeyerTSadeghiFChaiseT, et al.A coupled damage model and a semi-analytical contact solver to simulate butterfly wing formation around nonmetallic inclusions. Int J Fatigue2019; 127: 445-460.
3.
HouJDongJX. Characteristics of microscopic inclusions and their influence on fatigue crack propagation in FGH96 superalloy. Mater Sci Forum2017; 898: 505–516.
4.
MiaoJINYuLLiuX, et al.High temperature stress–strain curves of MnS and their applications in finite element simulation. Trans Nonferrous Met Soc China2018; 28: 2082–2093.
5.
LipaS. Numerical analysis of the effect of nonmetallic inclusions location in carburized layer of 16MNCR5 steel on initiation of fatigue cracking. Adv Sci Technol Res J2018; 12(3): 66-74.
6.
Arreola-HerreraRCruz-RamirezARivera-SalinasJE, 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. Theor Appl Fract Mech2018; 94: 134–146.
7.
DziobaILipiecS. Fracture mechanisms of S355 steel—experimental research, FEM simulation and SEM observation. Materials2019; 12: 3959.
8.
HoriMNemat-NasserS. Double-inclusion model and overall moduli of multi-phase composites. Mech Mater1993; 14: 189–206.
9.
ShodjaHSarvestaniA. Elastic fields in double inhomogeneity by the equivalent inclusion method. J Appl Mech2001; 68: 3–10.
10.
LuPLeongYPallathadkaP, et al.Effective moduli of nanoparticle reinforced composites considering interphase effect by extended double-inclusion model—theory and explicit expressions. Int J Eng Sci2013; 73: 33–55.
11.
EshelbyJD. The determination of the elastic field of an ellipsoidal inclusion, and related problems. Proc R Soc Lond A Math Phys Eng Sci1957; 241: 376–396.
12.
JacqCNéliasDLormandG, et al.Development of a three-dimensional semi-analytical elastic–plastic contact code. J Tribol2002; 124: 653.
13.
ChenWWWangQJWangF, et al.Three-dimensional repeated elasto-plastic point contacts, rolling, and sliding. J Appl Mech2008; 75: 021021.
14.
ZhouQJinXWangZ, et al.Numerical implementation of the equivalent inclusion method for 2D arbitrarily shaped inhomogeneities. J Elast2015; 118: 39–61.
15.
ZhouKKeerLWangQ, et al. Interaction of multiple inhomogeneous inclusions beneath a surface. Comput Methods Appl Mech Eng2012; 217–220: 25–33.
16.
AmuzugaKChaiseTDuvalA, et al.Fully coupled resolution of heterogeneous elastic–plastic contact problem. J Tribol2016; 138: 021403.
17.
ZhangMZhaoNGlawsP, et al. Elasto-plastic contact of materials containing double-layered inhomogeneities. Int J Solids Struct2017; 126–127: 208–224.
18.
WangZJinXKeerLM, et al.Novel model for partial-slip contact involving a material with inhomogeneity. J Tribol2013; 135: 1–15. doi:10.1115/1.4024548
19.
DongQZhouKChenW, et al.Partial slip contact modeling of heterogeneous elasto-plastic materials. Int J Mech Sci2016; 114: 98–110.
20.
WangZZhuDWangQ. Elastohydrodynamic lubrication of inhomogeneous materials using the equivalent inclusion method. J Tribol2013; 136: 021501.
21.
DongQLiYWeiJ, et al.Layered structures with rough surfaces and interfaces at contact loading. Int J Mech Sci2020; 178: 105611.
22.
DongQWeiJZhaoB, et al.Contact behaviors of multilayered structures with interfacial cracks. Eng Fract Mech2020; 225: 106816.
23.
YangWZhouQHuangY, et al.A thermoelastic contact model between a sliding ball and a stationary half space distributed with spherical inhomogeneities. Tribol Int2019; 131: 33–44.
24.
WangZJinXKeerLM, et al.A numerical approach for analyzing three-dimensional steady-state rolling contact including creep using a fast semi-analytical method. Tribol Trans2012; 55: 446–457.
25.
ZhangMZhaoNWangZ, et al.Efficient numerical method with a dual-grid scheme for contact of inhomogeneous materials and its applications. Comput Mech2018; 62: 991–1007. doi:10.1007/s00466-018-1543-3
26.
LiuSJinXWangZ, et al.Analytical solution for elastic fields caused by eigenstrains in a half-space and numerical implementation based on FFT. Int J Plast2012; 35: 135–154.
27.
ZhouQJinXWangZ, et al. A mesh differential refinement scheme for solving elastic fields of half-space inclusion problems. Tribol Int2016; 93: 124–136.
28.
ZhouQJinXWangZ, et al.Numerical EIM with 3D FFT for the contact with a smooth or rough surface involving complicated and distributed inhomogeneities. Tribol Int2016; 93: 91–103.
29.
ZhouKChenWKeerL, et al. Multiple 3D inhomogeneous inclusions in a half space under contact loading. Mech Mater2011; 43: 444–457.
30.
OlverPJ. Numerical analysis lecture notes. Published online 2008.
31.
KrylovAN. On the numerical solution of equation by which are determined in technical problems the frequencies of small vibrations of material systems. News Acad Sci USSR1931; 7: 491–539.
BrownPNSaadY. Hybrid Krylov methods for nonlinear systems of equations. SIAM J Sci Stat Comput1990; 11: 450–481.
34.
ChanTFJacksonKR. Nonlinearly preconditioned Krylov subspace methods for discrete Newton algorithms. SIAM J Sci Stat Comput1984; 5: 533–542.
35.
KnollDAKeyesDE. Jacobian-free Newton-Krylov methods: a survey of approaches and applications. J Comput Phys2004; 193: 357–397. doi:10.1016/j.jcp.2003.08.010
36.
WangFGuoHLiuW, et al.Control of MnS inclusions in high-and low-sulfur steel by tellurium treatment. Materials2019; 12: 1034.
37.
BruceTLongHDwyer-JoyceRS. Threshold maps for inclusion-initiated micro-cracks and white etching areas in bearing steel: the role of impact loading and surface sliding. Tribol Lett2018; 66: 111.
38.
BruceTRoundingELongH, et al.Characterisation of white etching crack damage in wind turbine gearbox bearings. Wear2015; 338: 164–177.
39.
Al-TameemiHALongHDwyer-JoyceRS. Initiation of sub-surface micro-cracks and white etching areas from debonding at non-metallic inclusions in wind turbine gearbox bearing. Wear2018; 406–407: 22–32. doi:10.1016/j.wear.2018.03.008
40.
Al-BedhanyJHLongH. Microscopic investigation of subsurface initiated damage of wind turbine gearbox bearings. J Phys Conf Ser IOP Publishing, 2018; 1106(1): 012029. doi:10.1088/1742-6596/1106/1/012029
