Next Article in Journal
Gene Expression Programming Model for Tribological Behavior of Novel SiC–ZrO2–Al Hybrid Composites
Next Article in Special Issue
Temperature during Repetitive Short-Term Operation of a Brake with Functionally Graded Friction Element
Previous Article in Journal
Tuning the Electronic and Charge Transport Properties of Schiff Base Compounds by Electron Donor and/or Acceptor Groups
Previous Article in Special Issue
Influence of the Additive of Ceramic and Intermetallic Powders on the Friction Properties and Temperature of the Wet Clutch Disc
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effect of the Transverse Functional Gradient of the Thin Interfacial Inclusion Material on the Stress Distribution of the Bimaterial under Longitudinal Shear

1
Department of Applied Mathematics, Faculty of Computer Science and Telecommunications, Cracow University of Technology, Warszawska Str. 24, 31-155 Cracow, Poland
2
Department of Applied Mathematics and Physics, Ukrainian Academy of Printing, Pidgolosko Str. 19, 79020 Lviv, Ukraine
3
Department of Mechanics and Applied Computer Science, Faculty of Mechanical Engineering, Bialystok University of Technology, Wiejska 45C, 15-351 Bialystok, Poland
*
Author to whom correspondence should be addressed.
Submission received: 3 November 2022 / Revised: 21 November 2022 / Accepted: 28 November 2022 / Published: 2 December 2022

Abstract

:
The effect of a functional gradient in the cross-section material (FGM) of a thin ribbon-like interfacial deformable inclusion on the stress–strain state of a piecewise homogeneous linear–elastic matrix under longitudinal shear conditions is considered. Based on the equations of elasticity theory, a mathematical model of such an FGM inclusion is constructed. An analytic–numerical analysis of the stress fields for some typical cases of the continuous functional gradient dependence of the mechanical properties of the inclusion material is performed. It is proposed to apply the constructed solutions to select the functional gradient properties of the inclusion material to optimize the stress–strain state in its vicinity under the given stresses. The derived equations are suitable with minor modifications for the description of micro-, meso- and nanoscale inclusions. Moreover, the conclusions and calculation results are easily transferable to similar problems of thermal conductivity and thermoelasticity with possible frictional heat dissipation.

1. Introduction

Many structural materials contain numerous thin inhomogeneities in the form of inclusions of different origins [1,2,3,4,5,6,7]. Quite often, these inclusions are used as elements to reinforce the structural parts of machines and structures or as fillers in composite materials or coatings [8,9,10,11,12,13,14,15]. The use of nanocomposites with specific properties in engineering and technology has significantly shifted the interest from the study of objects at the macro level (100–10−1 m) and micro level (10−3–10−6 m) to the nano level (10−9 m) [16,17,18,19,20].
One of the typical examples of composite materials is the structure with thin ribbon inclusions. Structural elements made with the use of FGM have proven to be rather effective in practice [21,22,23]. In this way, it is possible to achieve a significant improvement in their mechanical, rheological, thermal, or other properties or the formation of protective thin layers [18,24,25,26,27].
The mathematical modeling of nanostructures requires the construction of more complex constitutional laws in comparison with the macro level [28,29]. Therefore, it is important to construct methods for studying the stress–strain state of such structures. To model a thin inclusion, there are mainly two basic approaches using analytical methods. The first one is based on the use of Eshelby’s analytical solution [21,30,31] for an ellipsoidal inclusion, in which a limiting transition with a decrease in one of the characteristic dimensions of the inclusion is performed. However, its application to thin interphase inclusions is impossible. The second one is based on the principle of the conjugation of continua of different measurability [32,33,34,35,36,37] and the method of jump functions [28,29,38,39,40,41]. According to this method, the inclusions are replaced by a certain surface (in the two-dimensional case, a line) of the discontinuity of the physical–mechanical fields, which describes the perturbing effect of a thin inclusion. Successful attempts were made in [17,33,42,43,44,45] to apply it to consider the influence of various physical and contact nonlinearities in the antiplane problem of elasticity theory for two compressed half-spaces with interfacial defects. The frictional slip with possible heat generation for contacting bodies [34,44,46,47,48,49] and the boundary element approach [50,51,52] were also considered here. Inhomogeneity of the mechanical properties of structural materials can be both designed for a specific purpose (FGM) and a consequence of technological processes of obtaining new materials and their processing (FSW, ball-burnishing process, etc.) [53,54]. Such factors cause additional complexity in the constitutive relations for the mathematical modeling of the behavior of such materials. However, the case of a thin inclusion of an inhomogeneous material has not been practically studied.
The process of improving the mathematical models of FGM [11,25,26,55,56,57,58] is complicated by the complex geometry of structural elements and the consideration of imperfections in the contact of their components. This is especially important to ensure their qualitative design, both in terms of mechanical strength [11,59,60,61,62,63] and in terms of the consideration of thermal, magnetic, and piezoelectric load factors [13,21,24,44,53,54,58,64,65].
The works [17,31,35,39] have been devoted to the consideration of surface energy and stresses in nanocomposites. Consideration of the heterogeneity of inclusions’ properties at the micro and nano scale is particularly important because the heterogeneity density of matter (discrepancy variation in mechanical and other properties) of matter bodies with a decrease in their scale usually increases, and the impact of this heterogeneity increases even further.
This publication aims to develop the method of jump functions and to construct a convenient structured and highly versatile approach to study the longitudinal displacement and thermal heating of bodies with thin inclusions of arbitrary physical nature, including those made from FGM.

2. Formulation of the Problem

Consider an unbounded isotropic structure consisting of two half-spaces with elastic moduli G k   ( k = 1 , 2 ) , which is subjected to an external longitudinal shear load determined by uniformly distributed infinity stresses σ y z and σ y z , concentrated forces of intensity Q k and screw dislocations with vector Burgers component b k at points ς k S k ( k = 1 ,   2 ) . The stresses must satisfy the conditions σ x z 2 G 1 = σ x z 1 G 2 at infinity to ensure the straightness of the interface.
Let us investigate the stress–strain state (SSS) of the body section with a plane x O y perpendicular to the direction O z of its longitudinal displacement (external problem). The plane sections of half-spaces perpendicular to this axis form two half-planes S k   ( k = 1 , 2 ) and the abscissa axis corresponds to the interface L ~ x between them (Figure 1). At the interface of half-spaces (plane xOz), there is a tunnel section L = [ a ; a ] in the direction of the shear axis z, in which a certain object of general thickness 2 h   ( h a ) is inserted.
According to the paradigm of the method of jump functions [36], the presence of a thin inclusion in the bulk is modeled by jumps in the components of stress and displacement vectors in [38,40,41]:
[ σ y z ] h σ y z σ y z + = f 3 ( x ) ,   x L [ w x ] h w x w + x = [ σ x z G ] h σ x z G 1 σ x z + G 2 = f 6 ( x ) ;
f 3 ( x ) = f 6 ( x ) = 0 ,   if   x L .
It is hereinafter marked [ ] h = ( x , h ) ( x , + h ) ,   h = ( x , h ) + ( x , + h ) ; superscripts “+” and “−” correspond to the boundary values of the functions on the upper and lower banks of the line L .
The mathematical model of a thin inclusion is given as complicated conditions of imperfect contact between opposite matrix surfaces adjacent to the inclusion (internal problem) [28,38,39,40,65]. The general model of a thin, physically nonlinear inclusion is presented in [28,29,38], where the methods of modeling thin objects involve the integration of the defining relations describing the physical and mechanical properties of the material of the inclusion, with the subsequent consideration of the smallness of one of the linear dimensions of the inclusion.
Let us consider a similar model for a thin inclusion, assuming that the mechanical properties of the inclusion material are coordinate-dependent. This will allow us to model the inclusions of a functionally graded material:
{ G x i n ( x ) w i n x h ( x ) 2 σ x z i n ( a ) 1 h a x [ σ y z i n ] h ( ξ ) d ξ = 0 , G y i n ( x ) [ w i n ] h ( x ) + h σ y z i n h ( x ) = 0 .
Here, G x i n ( x ) , G y i n ( x ) are the variable shear moduli of the inclusion’s material. As a special case, considering their values to be constant, we obtain Hooke’s law. The upper index “in” denotes the terms describing the inclusion material’s SSS components.
Contact between matrix components and the inclusion at L and between the bimaterial structure components along a line L \ L is supposed to be mechanically perfect,
w ( x , ± h ) = w i n ( x , ± h ) , σ y z i n ( x , ± h ) = σ y z k ( x , ± h ) ,   x L , w ( x , + 0 ) = w ( x , 0 ) , σ y z 2 ( x , + 0 ) = σ y z 1 ( x , 0 ) ,    x L \ L ,
or frictional in some areas x L f L , as was considered in works [41,46,47,48],
σ y z i n ( x , ± h ) = σ y z k ( x , ± h ) = sgn [ w i n ] h τ y z max .
Here, τ y z K max is the limit value of shear stresses, at which the slippage begins. In this case, however, additional iterative methods should be applied to determine the area of the slip zones depending on the specific types of external loading of the composite [41].

3. Materials and Methods

Expressions for the components of the stress tensor and the derivatives of displacements on the line L of the infinite plane S = S 1 S 2 , as well as inside the latter, can be obtained by applying the results of [37] to the solution of the external problem
σ y z ( z ) + i σ x z ( z ) = σ y z 0 ( z ) + i σ x z 0 ( z ) + + i p k g 3 ( z ) C g 6 ( z ) ( z S k ;   r = 3 , 6 ;   k = 1 , 2 ) ;
σ y z k ± ( x ) = p k f 3 ( x ) C g 6 ( x ) + σ y z 0 ± ( x ) ,   σ x z k ± ( x ) = C f 6 ( x ) + p k g 3 ( x ) + σ x z 0 ± ( x ) ,  
where the notation [28,41] is introduced:
g r ( z ) 1 π L f r ( x ) d x x z , s r ( x ) a x f r ( x ) d x , p = 1 G 1 + G 2 ,   p k = G k p ,   C = G 1 G 2 p .
The superscript “+” corresponds to k = 2 and “–” corresponds to k = 1. Values marked with superscript “0” correspond to values in a continuous medium without inclusions under the same external load (homogeneous solution) [28,41].
Using (7), (8) and boundary conditions (4), it is easy to obtain from model (3) a system of singular integral equations:
{ ( p 2 p 1 ) f 6 ( x ) + 2 p g 3 ( x ) s 3 ( x ) h G x i n ( x ) = F 3 ( x ) , ( p 2 p 1 ) f 3 ( x ) + 2 C g 6 ( x ) G y i n ( x ) s 6 ( x ) h = F 6 ( x ) , F 3 ( x ) = 2 G x i n σ x z i n ( a ) ( σ x z 2 0 ( x ) / G 2 + σ x z 1 0 ( x ) / G 1 ) , F 6 ( x ) = σ y z 0 ( x ) G y i n σ y z k 0 ( x ) G k G y i n h [ w 0 ] ( a ) .
Balance conditions on the power balance and unambiguity of displacements while moving around the thin defect must be added to the solution of the external problem:
a a f 3 ( ξ ) d ξ = 2 h ( σ x z i n ( a ) σ x z i n ( a ) ) , a a f 6 ( ξ ) d ξ = [ w ] ( a ) [ w ] ( a ) .
Solving (9) and (10) using the methods in [29,38,41], it is easy to obtain a system of linear algebraic equations with unknown coefficients of the decomposition of the jump functions f r ( x ) into a series by orthogonal Jacobi or Chebyshev polynomials.
An important aspect of the study of the strength of such structures is the improvement of their strength criteria. In fracture mechanics, it is acceptable to use the stress intensity factor to describe the behavior of the SSS in the vicinity of the crack tip [42,45,61,62,63,66]. This is not sufficient for the case of a thin deformable inclusion. In [45], the authors obtained the two-term asymptotical expressions for the distribution of SSS in the vicinity of the thin inclusion tips using the introduced generalized stress intensity factors (GSIF):
K 31 + i K 32 = lim r 0 ( θ = 0 , π ) 2 π r ( σ y z + i σ x z ) .
Here, ( r , θ ) is a system of polar coordinates with the origin near the right or the left tip of the inclusion z = ± r exp ( i θ ) ± a .
Considering the well-known mathematical analogy [67], the obtained solutions to the antiplane problem can be regarded as solutions to the accordant heat conduction problem, if we take into account the correspondence of the values
w T ,   w x T x ,   w y T y ,   q x σ x z ,   q y σ y z ,   G x λ x ,   G y λ y ,   K 31 k q y ,   K 32 k q x .
The terms are as follows: T —temperature, q x ,   q y —heat flows, λ x ,   λ y —thermal conductivity coefficients, k q y ,   k q x —heat flow intensity factors [40].

4. Numerical Results and Discussion

Since the main focus of this article is to investigate the effect of the functional gradient on the mechanical properties of the inclusion material, we will limit ourselves to one of the most representative variants of the structure loading: homogeneous longitudinal shear σ y z = τ and σ x z k = τ k   ( k = 1 , 2 ) at infinity. However, the calculations for loading by concentrated force factors or dislocations do not make any fundamental difference except for the necessity to consider the locality of their application [29,41].
The dependence G x ( x ) , G y ( x ) on coordinate x for mathematical modeling can be defined as an arbitrary function (linear, exponential, power, periodic [58], etc.), which adequately reflects the desired practical properties of the material. To illustrate the method, let us consider one of the illustrative variants of the functional gradient of the inclusion material—the piecewise linear one:
G x ( x ) = G y ( x ) = { ( G 01 G 0 ) x a + G 01 , x [ a , 0 ] ; ( G 02 G 01 ) x a + G 01 , x [ 0 , a ] ,
where G 0 , G 01 , G 02 —some given constants.
To significantly reduce the number of calculations without loss of generality, it is convenient to use the following dimensionless quantities, marked by symbol “~” (tilde) on top:
x ˜ = x / a ,   h ˜ = h / a ,   y ˜ = y / a , G ˜ x i n ( x ˜ ) = G x i n ( x ) / G g a v ,   G ˜ y i n ( x ˜ ) = G y i n ( x ) / G g a v , τ ˜ k = τ k / G g a v ,   τ ˜ = τ / G g a v ,   G g a v = G 1 G 2 , G ˜ 0 = G 0 / G g a v , G ˜ 01 = G 01 / G g a v , G ˜ 02 = G 02 / G g a v , σ ˜ x z ( x ˜ ) = σ x z ( x ) / G g a v ,   σ ˜ y z ( x ˜ ) = σ y z ( x ) / G g a v . K ˜ 31 = K 31 + 2 C ˜ G g a v π a ,   K ˜ 32 = K 32 + 2 p 2 G g a v π a ,
where K 31 +   K 32 + are the GSIFs near the tip x = + a of the inclusion.
Figure 2, Figure 3, Figure 4, Figure 5, Figure 6, Figure 7, Figure 8, Figure 9, Figure 10 and Figure 11 illustrate the dependence of the stress–strain behavior of the matrix in the inclusion vicinity on the variation in the parameters G ˜ 0 , G ˜ 01 , G ˜ 02 , the values of which were chosen to reveal a qualitative picture of the FGM effect on the stress–strain parameters. It can be immediately concluded from Figure 2 and Figure 3 that under the load τ ˜ , the dimensionless K ˜ 31 + are expected to decrease with the increasing shear moduli of any part of the inclusion, while at K ˜ 32 + they appear to increase with increasing load τ ˜ k .
The effect of changes in the moduli G ˜ x ( x ) , G ˜ y ( x ) on the stresses σ ˜ y z , σ ˜ x z on the inclusion surface is more obvious if we choose a linear growth law for them along the inclusion axis (Figure 4, Figure 5 and Figure 6). The magnitude of the surface stresses increases significantly in the stiffer part of the inclusion. The larger the stiffness gradient, the more significant the increase.
The choice of the piecewise linear law of moduli G ˜ x ( x ) , G ˜ y ( x ) change in the Formulae (14) as G ˜ 01 = G ˜ 02 (variant 1) or G ˜ 0 = G ˜ 01 (variant 2) has a more contrasting effect on the surface stresses σ ˜ y z , σ ˜ x z , especially in the vicinity of the gradient breaking point x = 0 (Figure 7 and Figure 8). Moreover, variant 2 of the functional dependence of the inclusion material moduli leads to partial unloading in the softer part of the inclusion near the breaking point x = 0 (Figure 8).
Figure 9, Figure 10 and Figure 11 illustrate the changes in the stress field in the matrix in the inclusion vicinity under different variants of the law of functional change of the inclusion material moduli. The trends towards a decrease in the stress magnitudes in the vicinity of the stiffer parts of the inclusion are visible.

5. Conclusions

The proposed sufficiently simple and mathematically correct methodology made it possible for us to construct, for the first time, a mathematical model of a deformable thin linear interfacial inclusion with essentially inhomogeneous linear mechanical properties. Such a model can be used to simulate a thin inclusion from a functionally graded material and to solve the corresponding problems of defining the stress–strain field of the corresponding micro- or nanostructures by efficient analytical–numerical methods (the jump function method and its modifications), without the need to involve purely numerical approaches (in particular, FEM).
The calculations of the stress–strain field components for simple test cases of the functional dependence of the shear moduli of inclusion material have demonstrated the expected qualitative picture of their effect on the variation in the FGM parameters. In particular, (1) the stress magnitude increases significantly in the vicinity of the inclusion regions with increased stiffness; (2) the combination of the inclusion materials from parts with piecewise linear mechanical characteristics may lead to partial unloading of the inclusion and matrix in their softer part in the vicinity of the breaking point of the gradient dependence of the inclusion material parameters; (3) the contrast of the stress field changes of the inclusion and matrix is proportional to the increase in the gradient dependence.
The conclusions and calculation results are easily transferable to analogous problems of thermal conductivity and thermoelasticity with possible frictional heat generation and can be used for recommendations on the optimal operating parameters of structures.
The discussed conclusions can be useful in designing the functionally gradient mechanical properties of the material of inclusions and in the optimization of engineering structures to increase their strength and service life. The proposed method is effective for solving a wide class of problems of deformation of solids with thin deformable inclusions of finite length and can be used for SSS calculation for different FGM inclusions.

Author Contributions

Conceptualization, Y.P. and H.S.; methodology, Y.P. and H.S.; software, Y.P. and L.P.; validation, Y.P. and H.S.; formal analysis, Y.P.; investigation, Y.P. and H.S.; resources, H.S.; data curation, Y.P.; writing—original draft preparation, Y.P.; writing—review and editing, Y.P.; visualization, Y.P. and L.P.; supervision, H.S.; project administration, H.S.; funding acquisition, H.S. All authors have read and agreed to the published version of the manuscript.

Funding

This investigation was performed within the framework of research project No. 2017/27/B/ST8/01249, funded by the National Science Centre, Poland, and with project financing through the program of the Minister of Education and Science of Poland named “Regional Initiative of Excellence” in 2019–2022, project No. 011/RID/2018/19; amount of financing: 12,000,000 PLN.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are openly available at https://ua1lib.org/book/665574/5c937e (accessed on 9 January 2022), reference number [37]; https://0-doi-org.brum.beds.ac.uk/10.3390/ma15041435 (accessed on 9 January 2022), reference number [38]; and https://0-doi-org.brum.beds.ac.uk/10.3390/ma14174928 (accessed on 9 January 2022), reference number [41].

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

FGMfunctionally gradient material;
GSIFgeneralized stress intensity factor;
SSSstress–strain state
x , y , z Cartesian coordinates;
f r jump functions;
E k ,   ν k , G k ,   G x i n K ,   G y i n K elastic properties of the materials;
S k half-planes (sections of the solid);
w ,   σ x z ,   σ y z ,   σ x x ,   σ y y displacement, stresses (components of SSS);
L = [ a ; a ] line, modeling the presence of thin inclusion;
Q k ,   b k magnitudes of concentrated forces and screw dislocations;
σ y z , σ x z k uniformly distributed in infinity shear stresses.

References

  1. Mura, T.; Barnett, D.M. Micromechanics of Defects in Solids. J. Appl. Mech. 1983, 50, 477. [Google Scholar] [CrossRef]
  2. Wang, Y.; Huang, Z. Analytical Micromechanics Models for Elastoplastic Behavior of Long Fibrous Composites: A Critical Review and Comparative Study. Materials 2018, 11, 1919. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Wang, J.; Karihaloo, B.L.E.; Duan, H.L. Nano-mechanics or how to extend continuum mechanics to nanoscale. Bull. Pol. Acad. Sci. Tech. Sci. 2007, 55, 133–140. [Google Scholar]
  4. Zhou, K.; Hoh, H.J.; Wang, X.; Keer, L.M.; Pang, J.H.; Song, B.; Wang, Q.J. A review of recent works on inclusions. Mech. Mater. 2013, 60, 144–158. [Google Scholar] [CrossRef]
  5. Mencik, J. Mechanics of Components with Treated or Coated Solids; Kluwer Academic Publishing: Dordrecht, The Netherlands, 1996; 360p. [Google Scholar]
  6. Nemat-Nasser, S.; Hori, M. Micromechanics: Overall Properties of Heterogeneous Materials; North-Holland Series in Applied Mathematics and Mechanics; Elsevier: Amsterdam, The Netherlands, 1993; Volume 37, 687p. [Google Scholar]
  7. Pasternak, I.; Sulym, H.; Ilchuk, N. Interaction of Physicomechanical Fields in Bodies with Thin Structural Inhomogeneities: A Survey. J. Math. Sci. 2021, 253, 63–83. [Google Scholar] [CrossRef]
  8. Chen, J. Determination of thermal stress intensity factors for an interface crack in a graded orthotropic coating-substrate structure. Int. J. Fract. 2005, 133, 303–328. [Google Scholar] [CrossRef]
  9. Johnson, W.S.; Masters, J.; O’Brien, T.K.; Naik, R. Simplified Micromechanical Equations for Thermal Residual Stress Analysis of Coated Fiber Composites. J. Compos. Technol. Res. 1992, 14, 182–186. [Google Scholar] [CrossRef] [Green Version]
  10. Shevchuk, V.A. Modeling and Computation of Heat Transfer in a System “Body-Multilayer Coating”. Heat Transf. Res. 2006, 37, 412–423. [Google Scholar] [CrossRef]
  11. Kashtalyan, M.; Menshykova, M. Three-dimensional analysis of a functionally graded coating/substrate system of finite thickness. Philos. Trans. R. Soc. Lond. Ser. A Math. Phys. Eng. Sci. 2008, 366, 1821–1826. [Google Scholar] [CrossRef]
  12. Bao, G.; Wang, L. Multiple cracking in functionally graded ceramic/metal coatings. Int. J. Solids Struct. 1995, 32, 2853–2871. [Google Scholar] [CrossRef]
  13. Zhao, J.; Silberschmidt, V.V. Microstructure-based damage and fracture modelling of alumina coatings. Comput. Mater. Sci. 2005, 32, 620–628. [Google Scholar] [CrossRef]
  14. Hashin, Z. Thin interphase/imperfect interface in elasticity with application to coated fiber composites. J. Mech. Phys. Solids 2002, 50, 2509–2537. [Google Scholar] [CrossRef]
  15. Kulchytsky-Zhyhailo, R.; Matysiak, S.J.; Bajkowski, A.S. Semi-analytical solution of three-dimensional thermoelastic problem for half-space with gradient coating. J. Therm. Stress. 2018, 41, 1169–1181. [Google Scholar] [CrossRef]
  16. Davim, P.; Constantinos, A. (Eds.) Nanocomposites. Materials, Manufacturing and Engineering; Charitidis, Walter de Gruyter GmbH: Berlin, Germany; Boston, MA, USA, 2013; 211p. [Google Scholar]
  17. Kim, C.I.; Schiavone, P.; Ru, C.-Q. The Effects of Surface Elasticity on an Elastic Solid with Mode-III Crack: Complete Solution. Trans. ASME J. Appl. Mech. 2009, 77, 021011. [Google Scholar] [CrossRef]
  18. Uchida, Y. Properties of functionally graded materials, Manufactured by progressive lamination method for applications. Aichi Inst. Technol. Res. Rep. 2004, 39–B, 39–51. [Google Scholar]
  19. Ichikawa, K. (Ed.) Functionally Graded Materials in the 21ST Century, A Workshop on Trends and Forecasts; Kluver Academic Publishers: Boston, MA, USA; Dordrecht, The Netherlands; London, UK, 2000; 264p. [Google Scholar]
  20. Miyamoto, Y.; Kaysser, W.A.; Rabin, B.H.; Kawasaki, A.; Ford, R.G. Functionally Graded Materials: Design, Processing and Applications; Materials Technology Series; Springer: Berlin/Heidelberg, Germany, 1999; pp. 247–313. [Google Scholar] [CrossRef]
  21. Duan, H.; Wang, J.; Huang, Z.; Karihaloo, B. Eshelby formalism for nano-inhomogeneities. Proc. R. Soc. A Math. Phys. Eng. Sci. 2005, 461, 3335–3353. [Google Scholar] [CrossRef]
  22. Elperin, T.; Rudin, G. Thermal stresses in a coating–substrate assembly caused by internal heat source. J. Therm. Stress. 2016, 39, 90–102. [Google Scholar] [CrossRef]
  23. Zhang, N.; Khan, T.; Guo, H.; Shi, S.; Zhong, W.; Zhang, W. Functionally Graded Materials: An Overview of Stability, Buckling, and Free Vibration Analysis. Adv. Mater. Sci. Eng. 2019, 2019, 1354150. [Google Scholar] [CrossRef] [Green Version]
  24. Li, Y.; Feng, Z.; Hao, L.; Huang, L.; Xin, C.; Wang, Y.; Bilotti, E.; Essa, K.; Zhang, H.; Li, Z.; et al. A Review on Functionally Graded Materials and Structures via Additive Manufacturing: From Multi-Scale Design to Versatile Functional Properties. Adv. Mater. Technol. 2020, 5, 1900981. [Google Scholar] [CrossRef]
  25. Bishop, A.; Lin, C.Y.; Navaratnam, M.; Rawlings, R.D.; McShane, H.B. A functionally gradient material produced by a powder metallurgical process. J. Mater. Sci. Lett. 1993, 12, 1516–1518. [Google Scholar] [CrossRef]
  26. Kieback, B.; Neubrand, A.; Riedel, H. Processing techniques for functionally graded materials. Mater. Sci. Eng. A 2003, 362, 81–105. [Google Scholar] [CrossRef]
  27. Mistler, R.E. High strength alumina substrates produced by a multiple-layer casting technique. Am. Ceram. Soc. Bull. 1973, 52, 850–854. [Google Scholar]
  28. Sulym, H.T.; Piskozub, I.Z. Nonlinear Deformation of a Thin Interface Inclusion. Mater. Sci. 2018, 53, 600–608. [Google Scholar] [CrossRef]
  29. Sulym, H.; Piskozub, Y.; Polanski, J. Antiplane Deformation of a Bimaterial with Thin Interfacial Nonlinear Elastic Inclusion. Acta Mech. Autom. 2018, 12, 190–195. [Google Scholar] [CrossRef] [Green Version]
  30. Eshelby, J.D. The determination of the elastic field of an ellipsoidal inclusion, and related problems. Proc. Soc. Lond. A Math. Phys. Eng. Sci. 1957, 241, 376–396. [Google Scholar] [CrossRef]
  31. Sharma, P.; Ganti, S. Size-Dependent Eshelby’s Tensor for Embedded Nano-Inclusions Incorporating Surface/Interface Energies. ASME J. Appl. Mech. 2004, 71, 663–671. [Google Scholar] [CrossRef] [Green Version]
  32. Benveniste, Y.; Baum, G. An interface model of a graded three-dimensional anisotropic curved interphase. Proc. R. Soc. A Math. Phys. Eng. Sci. 2006, 463, 419–434. [Google Scholar] [CrossRef]
  33. Benveniste, Y.; Miloh, T. Imperfect soft and stiff interfaces in two-dimensional elasticity. Mech. Mater. 2001, 33, 309–323. [Google Scholar] [CrossRef]
  34. Martynyak, R.M.; Serednytska, K.I. Contact Problems of Thermoelasticity for Interface Cracks in Bimaterials; Rastr-7: Lviv, Ukraine, 2017; 168p. [Google Scholar]
  35. Gurtin, M.E.; Murdoch, A.I. A continuum theory of elastic material surfaces. Arch. Ration. Mech. Anal. 1975, 57, 291–323. [Google Scholar] [CrossRef]
  36. Steigmann, D.J.; Ogden, R.W. Elastic surface—Substrate interactions. Proc. R. Soc. A Math. Phys. Eng. Sci. 1999, 455, 437–474. [Google Scholar] [CrossRef]
  37. Sulym, H.T. Bases of Mathematical Theory of Thermo-elastic Equilibrium of Solids Containing Thin Inclusions; Research and Publishing Center of NTSh: L’viv, Ukraine, 2007; 716p, Available online: https://ua1lib.org/book/665574/5c937e (accessed on 19 January 2022). (In Ukrainian)
  38. Hutsaylyuk, V.; Piskozub, Y.; Piskozub, L.; Sulym, H. Deformation and Strength Parameters of a Composite Structure with a Thin Multilayer Ribbon-like Inclusion. Materials 2022, 15, 1435. [Google Scholar] [CrossRef] [PubMed]
  39. Piskozub, J.Z. Effect of surface stresses on the tensely deformed state of thin interface microinclusion. Math. Model. Comput. 2021, 8, 69–77. [Google Scholar] [CrossRef]
  40. Sulim, G.T.; Piskozub, J.Z. Thermoelastic equilibrium of piecewise homogeneous solids with thin inclusions. J. Eng. Math. Spec. Issue Thermomechanics 2008, 61, 315–337. [Google Scholar] [CrossRef]
  41. Piskozub, Y.; Sulym, H. Effect of Frictional Slipping on the Strength of Ribbon-Reinforced Composite. Materials 2021, 14, 4928. [Google Scholar] [CrossRef] [PubMed]
  42. Evtushenko, A.A.; Sulim, G.T. Stress concentration near a cavity filled with a liquid. Sov. Mater. Sci. 1981, 16, 546–549. [Google Scholar] [CrossRef]
  43. Lee, Y.-D.; Erdogan, F. Interface cracking of FGM coatings under steady-state heat flow. Eng. Fract. Mech. 1998, 59, 361–380. [Google Scholar] [CrossRef]
  44. Martynyak, R.M.; Kryshtafovych, A. Friction contact of two elastic half-planes with local recesses in boundary. J. Frict. Wear 2000, 21, 6–15. [Google Scholar]
  45. Piskozub, Y.; Sulym, H. Asymptotics of stresses in the vicinity of a thin elastic interphase inclusion. Mater. Sci. 1996, 32, 421–432. [Google Scholar] [CrossRef]
  46. Sulym, H.; Piskozub, L.; Piskozub, Y.; Pasternak, I. Antiplane Deformation of a Bimaterial Containing an Interfacial Crack with the Account of Friction I. Single Loading. Acta Mech. Autom. 2015, 9, 115–121. [Google Scholar] [CrossRef] [Green Version]
  47. Sulym, H.; Piskozub, L.; Piskozub, Y.; Pasternak, I. Antiplane Deformation of a Bimaterial Containing an Interfacial Crack with the Account of Friction 2. Repeating and Cyclic Loading. Acta Mech. Autom. 2015, 9, 178–184. [Google Scholar] [CrossRef] [Green Version]
  48. Sulym, H.; Pasternak, I.; Piskozub, L. Longitudinal shear of a bi-material with fractional sliding contact in the interfacial crack. J. Theor. Appl. Mech. 2016, 54, 529–539. [Google Scholar] [CrossRef] [Green Version]
  49. Yevtushenko, A.; Rozniakowska, M.; Kuciej, M. Transient temperature processes in composite strip and homogeneous foundation. Int. Commun. Heat Mass Transf. 2007, 34, 1108–1118. [Google Scholar] [CrossRef]
  50. Pasternak, I.M.; Pasternak, R.M.; Sulym, H.T. 2D boundary element analysis of defective thermoelectroelastic bimaterial with thermally imperfect but mechanically and electrically perfect interface. Eng. Anal. Bound. Elem. 2015, 61, 194–206. [Google Scholar] [CrossRef]
  51. Pasternak, Y.M.; Sulym, H.T.; Piskozub, L.G. Model tonkogo vklyucheniya v usloviyah yego idealnogo i neidealnogo kontaktnogo vzaimodeystviya s okruzhayuschim materialom. In Proceedings of the Trudy VI Mezhdunar. Simp. po Tribofatike MSTF 2010, Minsk, Belarus, 25 October–1 November 2010; BGU: Minsk, Belarus, 2010; Volume 2, Chapter 2. pp. 399–404. [Google Scholar]
  52. Chen, X.; Liu, Q. Thermal stress analysis of multi-layer thin films and coatings by an advanced boundary element method. Comput. Model. Eng. Sci. 2001, 2, 337–349. [Google Scholar]
  53. Rodríguez, A.; Calleja, A.; de Lacalle, L.N.L.; Pereira, O.; González, H.; Urbikain, G.; Laye, J. Burnishing of FSW Aluminum Al–Cu–Li Components. Metals 2019, 9, 260. [Google Scholar] [CrossRef] [Green Version]
  54. Sánchez, A.J.; Egea, A.; Rodríguez, A.; Celentano, D.; Calleja, A.; López de Lacalle, L.N. Joining metrics enhancement when combining FSW and ball-burnishing in a 2050 aluminum alloy. Surf. Coat. Technol. 2019, 367, 327–335. [Google Scholar] [CrossRef] [Green Version]
  55. Bouhamed, A.; Jrad, H.; Mars, J.; Wali, M.; Gamaoun, F.; Dammak, F. Homogenization of elasto-plastic functionally graded material based on representative volume element: Application to incremental forming process. Int. J. Mech. Sci. 2019, 160, 412–420. [Google Scholar] [CrossRef]
  56. Chen, Y.Z. Study of multiply-layered cylinders made of functionally graded materials using the transfer matrix method. J. Mech. Mater. Struct. 2011, 6, 641–657. [Google Scholar] [CrossRef] [Green Version]
  57. Jin, Z.-H. An asymptotic solution of temperature field in a strip of a functionally graded material. Int. Commun. Heat Mass Transf. 2002, 29, 887–895. [Google Scholar] [CrossRef]
  58. Majak, J.; Mikola, M.; Pohlak, M.; Eerme, M.; Karunanidhi, R. Modelling FGM materials. An accurate function approxi-mation algorithms. In Proceedings of the Modern Materials and Manufacturing (MMM 2021), Tallinn, Estonia, 27–29 April 2021; IOP Conference Series: Materials Science and Engineering. 2021; Volume 1140. [Google Scholar]
  59. Takahashi, M.; Itoh, Y.; Kashiwaya, H. Fabrication and evaluation of W/Cu gradient material by sintering and infiltration technique. In Proceedings of the First International Symposium on Functionally Gradient Materials FGM, Tokyo, Japan, 8–9 October 1990; Functionally Gradient Materials Forum, The Society of Non-Traditional Technology. 1990; pp. 129–134. [Google Scholar]
  60. El-Borgi, S.; Erdogan, F.; Ben Hatira, F. Stress intensity factors for an interface crack between a functionally graded coating and a homogeneous substrate. Int. J. Fract. 2003, 123, 139–162. [Google Scholar] [CrossRef]
  61. Erdogan, F.; Wu, B.H. Crack problem in FGM layers under thermal stresses. J. Therm. Stress. 1996, 19, 237–265. [Google Scholar] [CrossRef]
  62. Wang, B.L.; Mai, Y.W.; Noda, N. Fracture mechanics analysis model for functionally graded materials with arbitrarily dis-tributed properties. Int. J. Fract. 2002, 116, 161–177. [Google Scholar] [CrossRef]
  63. Yildirim, B.; Dag, S.; Erdogan, F. Three dimensional fracture analysis of FGM coatings under thermomechanical loading. Int. J. Fract. 2005, 132, 371–397. [Google Scholar] [CrossRef]
  64. Guo, L.-C.; Noda, N.; Ishihara, M. Thermal Stress Intensity Factors for a Normal Surface Crack in a Functionally Graded Coating Structure. J. Therm. Stress. 2007, 31, 149–164. [Google Scholar] [CrossRef]
  65. Sulym, H.T.; Piskozub, Y.Z. Conditions of contact interaction (a survey). Mat. Metod. Fiz. Mekhanichni Polya 2004, 47, 110–125. [Google Scholar]
  66. Ding, S.-H.; Li, X. Thermal stress intensity factors for an interface crack in a functionally graded layered structures. Arch. Appl. Mech. 2010, 81, 943–955. [Google Scholar] [CrossRef]
  67. Kit, G.S. On Analogy between Longitudinal Shear and Stationary Thermal Conductivity of Bodies with Inclusions and Cracks. Dokl. USSR Ser. A 1977, 4, 336–340. [Google Scholar]
Figure 1. The loading and geometric scheme of the problem.
Figure 1. The loading and geometric scheme of the problem.
Materials 15 08591 g001
Figure 2. Influence of the parameters G ˜ 0 , G ˜ 01 , G ˜ 02 on the GSIF K ˜ 31 + under the load, uniformly distributed on infinity stress σ y z = τ .
Figure 2. Influence of the parameters G ˜ 0 , G ˜ 01 , G ˜ 02 on the GSIF K ˜ 31 + under the load, uniformly distributed on infinity stress σ y z = τ .
Materials 15 08591 g002
Figure 3. Influence of the parameters G ˜ 0 , G ˜ 01 , G ˜ 02 on the GSIF K ˜ 32 + under the load, uniformly distributed on infinity stress σ x z k = τ k ( k = 1 , 2 ) .
Figure 3. Influence of the parameters G ˜ 0 , G ˜ 01 , G ˜ 02 on the GSIF K ˜ 32 + under the load, uniformly distributed on infinity stress σ x z k = τ k ( k = 1 , 2 ) .
Materials 15 08591 g003
Figure 4. Stress distribution along with the upper interface (inclusion–matrix half-space S 2 ) with a linear distribution of material stiffness.
Figure 4. Stress distribution along with the upper interface (inclusion–matrix half-space S 2 ) with a linear distribution of material stiffness.
Materials 15 08591 g004
Figure 5. Stress distribution along with the upper interface (inclusion–matrix half-space S 2 ) with a linear distribution of material stiffness.
Figure 5. Stress distribution along with the upper interface (inclusion–matrix half-space S 2 ) with a linear distribution of material stiffness.
Materials 15 08591 g005
Figure 6. Stress distribution along with the upper interface (inclusion–matrix half-space S 2 ) with a linear distribution of material stiffness.
Figure 6. Stress distribution along with the upper interface (inclusion–matrix half-space S 2 ) with a linear distribution of material stiffness.
Materials 15 08591 g006
Figure 7. Stress distribution along with the upper interface (inclusion—matrix half-space S 2 ) with a piecewise linear distribution of material stiffness.
Figure 7. Stress distribution along with the upper interface (inclusion—matrix half-space S 2 ) with a piecewise linear distribution of material stiffness.
Materials 15 08591 g007
Figure 8. Stress distribution along with the upper interface (inclusion—matrix half-space S 2 ) with a piecewise linear distribution of material stiffness.
Figure 8. Stress distribution along with the upper interface (inclusion—matrix half-space S 2 ) with a piecewise linear distribution of material stiffness.
Materials 15 08591 g008
Figure 9. Stress distribution in the matrix at the vicinity of the inclusion with a linear distribution of material stiffness.
Figure 9. Stress distribution in the matrix at the vicinity of the inclusion with a linear distribution of material stiffness.
Materials 15 08591 g009
Figure 10. Stress distribution in the matrix at the vicinity of the inclusion with the piecewise linear distribution of material stiffness.
Figure 10. Stress distribution in the matrix at the vicinity of the inclusion with the piecewise linear distribution of material stiffness.
Materials 15 08591 g010
Figure 11. Stress distribution in the matrix at the vicinity of the inclusion with the piecewise linear distribution of material stiffness.
Figure 11. Stress distribution in the matrix at the vicinity of the inclusion with the piecewise linear distribution of material stiffness.
Materials 15 08591 g011
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Piskozub, Y.; Piskozub, L.; Sulym, H. Effect of the Transverse Functional Gradient of the Thin Interfacial Inclusion Material on the Stress Distribution of the Bimaterial under Longitudinal Shear. Materials 2022, 15, 8591. https://0-doi-org.brum.beds.ac.uk/10.3390/ma15238591

AMA Style

Piskozub Y, Piskozub L, Sulym H. Effect of the Transverse Functional Gradient of the Thin Interfacial Inclusion Material on the Stress Distribution of the Bimaterial under Longitudinal Shear. Materials. 2022; 15(23):8591. https://0-doi-org.brum.beds.ac.uk/10.3390/ma15238591

Chicago/Turabian Style

Piskozub, Yosyf, Liubov Piskozub, and Heorhiy Sulym. 2022. "Effect of the Transverse Functional Gradient of the Thin Interfacial Inclusion Material on the Stress Distribution of the Bimaterial under Longitudinal Shear" Materials 15, no. 23: 8591. https://0-doi-org.brum.beds.ac.uk/10.3390/ma15238591

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop