Next Article in Journal
PCEP: Few-Shot Model-Based Source Camera Identification
Previous Article in Journal
Flexible Metamaterial-Based Frequency Selective Surface with Square and Circular Split Ring Resonators Combinations for X-Band Applications
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Interfacial Stresses for a Coated Irregularly Shaped Hole Embedded in an Infinite Solid under Point Heat Singularity

Department of Mechanical Engineering, National Taiwan University of Science and Technology, Taipei 106335, Taiwan
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(4), 802; https://doi.org/10.3390/math11040802
Submission received: 5 January 2023 / Revised: 2 February 2023 / Accepted: 2 February 2023 / Published: 4 February 2023

Abstract

:
Particle-reinforced metals are being developed for advanced heat dissipation applications. However, an irregularly shaped void develops during eutectic solidification and enhances interfacial stress induced by visco-plastic deformation in temperature gradient conditions. An analytical solution to an irregularly shaped coated hole embedded in an infinite substrate under an arbitrarily located heat source or sink is presented. For a coated polygonal hole with any number of edges, a rapidly convergent series solution of the temperature and stress functions is expressed in an elegant form using conformal mapping, the analytic continuation theorem, and the alternation method. The iterations of the trial-and-error method are utilized to obtain the solution for the correction terms. First, temperature contours are obtained to provide an optimal suggestion that a larger thermal conductivity of the coating layer exhibits better heat absorption capacity. Furthermore, interfacial stresses between a coating layer and substrate increase if the strength of a point thermal singularity and thermal mismatch increases. This study provides a detailed explanation for the growth of an irregular void at an ambient temperature gradient.

1. Introduction

Metal matrix composites (MMCs) are used in the automotive industry or can power electronic devices owing to their heterogeneous microstructure or remarkably high thermal fatigue resistance [1,2], respectively. To exhibit higher energy absorption capacity and reduced weight, dual-phase steel [3,4,5] is integrated with hard martensite lath and softer ferrite matrix and is responsible for reinforcing particles; whereas ductility is incorporated by the matrix. Regarding the heat produced on ceramic chips in electronic devices, an excellent thermal management solution is required. To avoid delamination during temperature gradients, advanced MMCs [6,7] are being developed by combining the high thermal conductivity of the matrix and low thermal expansion of ceramic. Although the mentioned MMC exhibits excellent performance when subjected to thermal and deformation conditions, its drawback includes low wettability, weak interface bonding strength, and thermal mismatch, which are critical for void formation. Consequently, MMCs suffer thermal fatigue damage under temperature gradient conditions [2,8].
Observing from the experiment, the silicon content [9] or carbon diamond [10] in the matrix performs irregularly shaped voids during eutectic solidification. Subsequently, the pores grow further owing to larger plastic deformation and thermal mismatch during cooling to room temperature solidification. To realize the porosity behavior in MMC by the finite element method, Rivera-Salinas et al. [11] studied the effect of porosity in SiC particles embedded in the Al matrix and concluded that a stress concentrator appears around the contact points among the particle and matrix. Both strain and stress in MMC material during thermal loading have been investigated in detail. As is known, interfacial stress plays a dominant role in increasing the driving force [12] to cause damage initiation and crack propagation. However, a comprehensive quantitative analysis of interfacial stress under temperature loading in MMC is challenging. The difficulty in discovering interfacial stress lies in the fact that the detected interfacial stresses are disturbed by neighbor material in the experimental method; moreover, the interfacial stresses are interfered with by heterogeneous microstructure, thus displaying an unstable solution in the simulation method. Hence, the purpose of this investigation was to theoretically and quantitatively determine the evolution of interfacial stress and temperature distribution under a temperature gradient environment in composite material, including void, inclusion, and matrix.
For a micromechanical analytical model of porous solids with particles, Ramakrishnan and Arunachalam [13] predicted the effective elastic moduli of composite material as a function of the volume of pores using the continuum mechanics method. Additional theoretical models for porous solids and validation with experiments can be discovered in this article [14]. As mentioned above, deriving the analytical model of the interaction of multiple material systems is difficult because the stress distribution around pores, particles, and matrix disturbs each other [11]. Nonetheless, the quantitative behavior of the interface can be investigated under the assumption of continuation conditions across the interfaces from complex variable theory [15]. An analytical solution to a tri-material of circularly cylindrical layered [16], fibrous composite [17], and an elliptical coating hole [18] were investigated under a point thermal singularity condition. The aforementioned studies are all static situations, although many researchers have remarked on the dynamic analysis subject to moving thermal loads [19,20]. Nevertheless, the investigation of simple geometry as a circle boundary cannot comprehensively achieve a realistic microstructure material. In experimental observation, different shape inclusions can be microstructurally identified in MMC as transformation induced by plasticity-steel (TRIP) with 5% Mg-YSZ particles. In Figure 1, the ceramic particles of Mg-YSZ can be recognized as a triangle (Figure 1a), quadrilateral (Figure 1b), and polygon (Figure 1c) shaped inclusions [21,22]. Hence, numerous researchers studied the irregular shape inclusion problem theoretically [23,24,25,26,27].
The aforementioned studies primarily investigated the problems involving only two phases, namely the inclusion and matrix. Moreover, microscopic observation of the Al matrix with SiC [9] or carbon diamond [10] particles and 20NiCr3H steel [28] showed an irregular pore, inclusion, and matrix. Furthermore, the microscopic observation revealed an irregular void, iron inclusion, and steel matrix, as shown in Figure 1d. In addition, an irregular void, inclusion (SiC), and matrix (AlSi7) were detected by Schöbel et al. [9,10]. According to these geometric suggestions, an irregularly shaped hole with a coating layer embedded in the matrix is established, as shown in Figure 2a. A hole can be observed to exist inside and attach to the coating layer embedded in the matrix. In the theoretical model limitation, an “equiangular polygon” is investigated for realizing irregularly shaped inclusion. This composite system, including irregular holes, coating layers, and matrix has been discussed under different conditions [29,30,31,32,33]. However, the previous study [33] can only be utilized to deal with certain triangular or square-coated hole problems under thermal loading. Hence, the general solution of temperature and stress function and singular terms of an irregularly shaped pore with an arbitrary edge and a coating layer are investigated for the first time in this study. Many studies have demonstrated that the majority of materials experience thermal stress (tension or compression) and strain when subjected to different thermal loadings [34,35]. A qualitative and quantitative investigation of thermal stress on the interface is presented in this study.
In this study, the temperature distribution was solved first with a visualization of the temperature contour. The next step is to determine the proper interfacial stress of a coated pentagonal hole subjected to point thermal singularity. Additionally, parameter analysis is performed, including the shape factor, shear modulus, and thermal expansion coefficient. The remaining of this study is divided into the following sections. Section 2.1 provides the flow diagram of this problem. The general formulation of isotropic thermo-elasticity in two dimensions, as well as mapping techniques, are introduced in Section 2.2. Analytical solutions in the general form of the temperature and interfacial stress of the coated polygonal layer are provided in Section 3.1 and Section 3.2, respectively. The correction terms required to rectify the singular terms are determined using the trial-and-error method in MATLAB, which is described in Section 3.3. In Section 4, practical applications are discussed through a quantitative analysis. Section 5 presents the conclusions.

2. Mathematical Modeling

2.1. Flow Diagram

A flow diagram of the study is shown in Figure 3. The theoretical approach, which comprises the mapping method, analytical continuation, and alternative methodology, is described first. The temperature field is analyzed utilizing theoretical techniques and iterations. To derive the stress function and correction terms, the iterative method and the trial-and-error method are then suggested. The numerical computations, including those for temperature contours and interfacial stresses, are therefore carried out using MATLAB. Finally, output results are acquired and thoroughly analyzed.
An analytical solution to the temperature function of the coated polygonal hole can be obtained by using the iterative technique in the temperature field. The homogenous solution of the temperature field is initially introduced. The continuity conditions (temperature and resultant heat flux) are then considered to be satisfied by the outer interface for the first iteration of the temperature field. The inner interface must fulfill the adiabatic condition in order for the second iteration to proceed. It is possible to achieve an analytical solution for the temperature field by repeating the aforementioned iterative techniques until n iterations.
Determining an analytical solution to the stress function in the stress field can be aided by the iterative technique and the trial-and-error method. The homogeneous solution for the stress field is given initially and the unknown solution for the stress field is determined using the analytical continuation theorem. After that, it is ensured that the outer interface should satisfy the continuity conditions (displacement and resultant force) for the first iteration of the stress field. Meanwhile, MATLAB is used to implement a while-loop and the trial-and-error method for determining the correction terms of the outer interface. The traction-free condition is then satisfied by the inner interface and correction terms are applied to derive the stress function of the inner interface. The analytical solution to the stress field can be obtained by repeating the abovementioned iterative steps n times.
The numerical computations are implemented by MATLAB, such as temperature contours and interfacial stresses, after obtaining the analytical solutions for the temperature and stress fields. The final step is undertaking a thorough quantitative and qualitative analysis to compare the output findings and offer reasonable explanations for pertinent practical situations.

2.2. Problem Construction

Figure 1 shows several irregularly shaped inclusions, which must be discussed in detail. To further explore the problem of a coated polygonal hole, we consider an equiangular polygonal hole with a coating layer inlaid with an infinite plane that is subjected to point thermal singularity. The general solution of correction terms and temperature and stress fields are described in this study. Consider that the substrate Ω a and coating layer Ω b are assumed to be perfectly bonded across the interface Γ 1 . The interface Γ 2 shown in Figure 2a between the coating layer Ω b and the hole must satisfy the adiabatic and traction-free conditions in temperature and stress fields, respectively. The two undetermined stress functions ( ϕ z ,   ψ z ) are utilized to determine the displacement and the resultant force. Assume that θ z = g z is the undetermined temperature function. The relevant formula of isotropic thermo-elasticity in two-dimension is given as follows [15,36].
2 G u x + i u y = κ ϕ z z ϕ z ¯ ψ z ¯ + 2 G β g z d z
F y + i F x = ϕ z + z ϕ z ¯ + ψ z ¯
For plane strain deformation, κ = 3 4 ν and β = 1 + ν α is considered. The parameters G , ν , and α denote the shear modulus, Poisson’s ratio and the coefficient of thermal expansion, respectively. Additionally, the function (   ) represents the complex conjugate, and the prime sign implies the derivative associated with z = x + i y . To address the current problem, the mapping function is implemented as follows [15]:
z = m ζ = R ζ + w ζ λ , 0 w < 1 λ
After using the mapping function, the coated polygonal hole in the physical plane (see Figure 2a) with two interfaces Γ 1 and Γ 2 maps onto the concentric circles consisting of two interfaces L   1 and L   2 in mathematical plane (see Figure 2b). The inside and outside radii of the concentric circles are set to be unity and r , respectively in the ζ p l a n e , as shown in Figure 2b. Assume that S a and S b respectively represents the substrate the coating layer with the interface L   1 after the transformation. The interface L   2 shown in Figure 2b is between the coating layer S b and a unit circular hole. Moreover, the parameter λ + 1 describes the number of edges on the coated polygonal hole in Equation (3). For example, λ = 2 denotes a triangle; λ = 3 denotes a square, and λ = 4 denotes a pentagon.
Additionally, the magnification factor R in Equation (3) is related to particle size. Although “particle size” is being conducted by many researchers [37,38], the parameter R was assumed to be unity in this study to simplify the complicated geometry. Another parameter in Equation (3) that might affect the sharpness of polygon corners is the shape factor w . The different shape factors w of a pentagonal hole ( λ = 4 ) with a coating layer are illustrated in Figure 4. Generally, the edge of the pentagonal hole becomes sharper when w is larger. However, the range of the shape factors w is defined in Equation (3). The parameter w  ranges between 0 to 1 / λ . When w = 1 / λ , the polygonal problem will develop into a hypocycloidal crack. Hence, this study only discusses the range of 0 w < 1 / λ to avoid singularity with w = 1 / λ .
After applying the mapping function, Equations (1) and (2) can be transferred to mathematical plane form, rewritten as follows:
2 G u x + i u y = κ ϕ ζ m ζ m ζ ¯ ϕ ζ ¯ ψ ζ ¯ + 2 G β g ζ m ζ d ζ
F y + i F x = ϕ ζ + m ζ m ζ ¯ ϕ ζ ¯ + ψ ζ ¯

3. Mathematical Formulation

3.1. Temperature Field

Consider that the temperature function fulfills the homogeneous differential equation (Laplace equation) of the steady-state heat conduction problem [36]. The complex potential θ z = g z is related to temperature T and heat flux h as
T = Re θ ζ
h = q x d y q y d x = k Im θ ζ
where Re and Im respectively represent the real and imaginary parts in Equations (6) and (7). The components q x and q y in Equation (7) are the x and y- directions of the heat flux, respectively. The parameter k is the coefficient of the heat conductivity.
The homogeneous solution of the temperature field subject to point thermal singularity in an infinite substrate is written as follows [36]
θ 0 ζ = Q l o g ζ ζ 0
where Q = q 0 / 2 π k a (or Q = q 0 / 2 π k a ) is defined as a point heat source (or a point heat sink); q 0 and ζ 0 are the strength and position of the point thermal singularity, respectively. The temperature functions of a coated hole under a point thermal singularity are assumed as
θ ζ = θ 0 ζ + n = 1 θ a n ζ   ζ S a n = 1 θ b n ζ + n = 1 θ n ζ ζ S b
where θ a n ζ is holomorphic in S a region ( ζ r ); θ n ζ and θ b n ζ are both holomorphic in S b region ( ζ r ).
Using the analytical continuation and alternation method, the temperature functions in Equation (9) can be obtained in terms of θ 0 ζ as
θ ζ = θ 0 ζ + V b a θ 0 ¯ r 2 ζ + U a b U b a n = 1 V a b n 1 θ 0 ¯ 1 r 2 n 1 1 ζ   , ζ S a U b a θ 0 ζ + V a b U b a n = 1 V a b n 1 θ 0 1 r 2 n 1 ζ r 2 + U b a n = 1 V a b n 1 θ 0 ¯ 1 r 2 n 1 1 ζ         , ζ S b
where
V a b = k b k a / k a + k b
U b a = 2 k a / k a + k b
Integration of Equation (10) with ζ yields:
g ζ = Q R λ 1 w ζ 0 λ V b a r 2 ζ 0 ¯ U a b U b a r 2 ζ 0 ¯ ( r 2 V a b ) log ζ + g a i ζ + g a o ζ   , ζ S a U b a Q R λ 1 w ζ 0 λ + V a b w ζ 0 λ ( r 2 λ V a b ) r 2 ζ 0 ¯ ( r 2 V a b ) log ζ + g b i ζ + g b o ζ   , ζ S b
where the temperature functions g a i ζ and g b i ζ are holomorphic in the region S b ; and g a o ζ , g b o ζ are holomorphic in region S a . It should be noted that the solutions in Equation (17) are the series solution for arbitrarily shaped coated polygonal holes. The temperature functions g a i ζ , g a o ζ , g b i ζ , and g b o ζ in Equations (13) are respectively derived in Appendix A.

3.2. Stress Field

To facilitate the derivation of boundary conditions, an auxiliary stress function ω ζ is introduced as
ω ζ = m ¯ r 2 ζ m ζ ϕ ζ + ψ ζ
In order to satisfy the single-valued conditions of displacement and resultant force, Equations (4) and (5), the stress functions of the present problem must be formulated as follows:
ϕ ζ = M log ζ + n = 1 ϕ a n ζ       , ζ S a N log ζ + n = 1 ϕ b n ζ + n = 1 ϕ n ζ   , ζ S b
ω ζ =   M ¯ log ζ + n = 1 ω a n ζ       , ζ S a   N ¯ log ζ + n = 1 ω b n ζ + n = 1 ω n ζ   , ζ S b
where ϕ a n ζ and ω a n ζ are holomorphic in S a region ( ζ r ); ϕ n ζ , ω n ζ , ϕ b n ζ , and ω b n ζ are holomorphic in S b regions ( ζ r ).
The general expressions of M and N for different values of λ are given by:
M = 2 G a β a ( κ a + 1 ) Q R λ 1 U a b U b a r 2 ζ 0 ¯ ( r 2 V a b ) + V b a r 2 ζ 0 ¯ w ζ 0 λ
N = 2 G b β b ( κ b + 1 ) U b a Q R λ 1 r 2 ζ 0 ¯ ( r 2 V a b ) V a b w ζ 0 λ ( r 2 λ V a b ) w ζ 0 λ
In Equation (15), the coefficients M and N for the logarithmic terms are derived in Appendix A.
First, the stress functions at the interface L 1 must satisfy the continuity conditions. The general form of the first analytical continuation at interface L 1 yields:
ϕ a 1 ( ζ ) = 2 G a G b κ a G b + G a β b g b o ( ζ ) β a g a o ( ζ ) ω b 1 ( ζ ) = 2 G a G b κ a G b + G a β b g b o ¯ r 2 ζ β a g a o ¯ r 2 ζ m = 1 λ + 1 C 1 m ¯ r 2 ζ r 2 ( λ w ) 1 / λ + 1 e 2 ( m 1 ) i π / λ + 1 ϕ b 1 ( ζ ) = 2 G a G b κ b G a + G b β a g a i ( ζ ) β b g b i ( ζ ) m = 1 λ + 1 Π a b C 1 m ζ r 2 ( λ w ) 1 / λ + 1 e 2 ( m 1 ) i π / λ + 1 ω a 1 ( ζ ) = 2 G a G b κ b G a + G b β a g a i ¯ r 2 ζ β b g b i ¯ r 2 ζ m = 1 λ + 1 ( 1 + Π a b ) C 1 m ¯ r 2 ζ r 2 ( λ w ) 1 / λ + 1 e 2 ( m 1 ) i π / λ + 1
where
Π a b = G a G b / G a κ b + G b
C 1 m = 1 + r 2 λ + 2 λ w 2 λ λ 1 w 2 λ 1 / λ + 1 λ + 1 r 2 λ 2 e 4 ( m 1 ) i π / λ + 1 ϕ b 1 ¯ ( λ w ) 1 / λ + 1 e 2 ( m 1 ) i π / λ + 1 , m = 1 , 2 , , λ + 1
the interface L 2 between the coating layer and the hole must fulfill the traction-free condition. The general form of the second analytical continuation at interface L 2 yields:
ϕ 1 ( ζ ) = ω b 1 * ¯ 1 ζ m = 1 λ + 1 C 1 m * ζ ( λ w ) 1 / λ + 1 e 2 ( m 1 ) i π / λ + 1 ω 1 ( ζ ) = ϕ b 1 ¯ 1 ζ + m = 1 λ + 1 C 1 m * ¯ 1 ζ ( λ w ) 1 / λ + 1 e 2 ( m 1 ) i π / λ + 1
where
ω b 1 * ζ = m ¯ 1 ζ m ¯ r 2 ζ m ζ ϕ b 1 ζ + ω b 1 ζ
C 1 m * = 1 + 1 λ w 2 λ λ 1 w 2 λ 1 / λ + 1 λ + 1 e 4 ( m 1 ) i π / λ + 1 ϕ b 1 ¯ ( λ w ) 1 / λ + 1 e 2 ( m 1 ) i π / λ + 1 , m = 1 , 2 , , λ + 1
Third, returning to the stress functions at the interface L 1 , they must satisfy the continuity conditions at the interface L 1 . The general solution of the third analytical continuation at interface L 1 yields:
ϕ a 2 ( ζ ) = ( 1 + Λ a b ) ϕ 1 ( ζ ) ω b 2 ( ζ ) = Λ a b ϕ 1 ¯ r 2 ζ m = 1 λ + 1 C 2 m ¯ r 2 ζ r 2 ( λ w ) 1 / λ + 1 e 2 ( m 1 ) i π / λ + 1 ϕ b 2 ( ζ ) = Π a b ω 1 * ¯ r 2 ζ m = 1 λ + 1 Π a b C 2 m ζ r 2 ( λ w ) 1 / λ + 1 e 2 ( m 1 ) i π / λ + 1 ω a 2 ( ζ ) = ( 1 + Π a b ) ω 1 * ( ζ ) m = 1 λ + 1 ( 1 + Π a b ) C 2 m ¯ r 2 ζ r 2 ( λ w ) 1 / λ + 1 e 2 ( m 1 ) i π / λ + 1
where
Λ a b = G a κ b G b κ a G b κ a + G a
ω 1 * ζ = m ¯ r 2 ζ m ¯ 1 ζ m ζ ϕ 1 ζ + ω 1 ζ
C 2 m = 1 + r 2 λ + 2 λ w 2 λ λ 1 w 2 λ 1 / λ + 1 λ + 1 r 2 λ 2 e 4 ( m 1 ) i π / λ + 1 ϕ b 2 ¯ ( λ w ) 1 / λ + 1 e 2 ( m 1 ) i π / λ + 1 , m = 1 , 2 , , λ + 1
After repeating the analytical continuation in association with the alternation method, the remaining stress functions are determined as;
ϕ a n ( ζ ) = ( 1 + Λ a b ) ϕ n 1 ( ζ ) ω b n ( ζ ) = Λ a b ϕ n 1 ¯ r 2 ζ m = 1 λ + 1 C n m ¯ r 2 ζ r 2 ( λ w ) 1 / λ + 1 e 2 ( m 1 ) i π / λ + 1 ϕ b n ( ζ ) = Π a b ω n 1 * ¯ r 2 ζ m = 1 λ + 1 Π a b C n m ζ r 2 ( λ w ) 1 / λ + 1 e 2 ( m 1 ) i π / λ + 1 ω a n ( ζ ) = ( 1 + Π a b ) ω n 1 * ( ζ ) m = 1 λ + 1 ( 1 + Π a b ) C n m ¯ r 2 ζ r 2 ( λ w ) 1 / λ + 1 e 2 ( m 1 ) i π / λ + 1 n = 3 , 4 , 5
ϕ n ( ζ ) = ω b n * ¯ 1 ζ m = 1 λ + 1 C n m * ζ ( λ w ) 1 / λ + 1 e 2 ( m 1 ) i π / λ + 1 ω n ( ζ ) = ϕ b n ¯ 1 ζ + m = 1 λ + 1 C n m * ¯ 1 ζ ( λ w ) 1 / λ + 1 e 2 ( m 1 ) i π / λ + 1 n = 2 , 3 , 4
where
C n m = 1 + r 2 λ + 2 λ w 2 λ λ 1 w 2 λ 1 / λ + 1 λ + 1 r 2 λ 2 e 4 ( m 1 ) i π / λ + 1 ϕ b n ¯ ( λ w ) 1 / λ + 1 e 2 ( m 1 ) i π / λ + 1 , m = 1 , 2 , , λ + 1
C n m * = 1 + 1 λ w 2 λ λ 1 w 2 λ 1 / λ + 1 λ + 1 e 4 ( m 1 ) i π / λ + 1 ϕ b n ¯ ( λ w ) 1 / λ + 1 e 2 ( m 1 ) i π / λ + 1 , m = 1 , 2 , , λ + 1

3.3. Trial-and-Error Method

After the mathematical solution is obtained, the numerical calculation is implemented using MATLAB. The most difficult part of exporting numerical results is dealing with correction terms that appear in the stress field. According to the flow diagram (see Figure 3), the zero iteration for the stress field is the homogeneous solution ϕ 0 ζ and ψ 0 ζ . The first iteration happens on interface L 1 , where the stress functions are ϕ a 1 ζ ,   ω a 1 ζ ,   ϕ b 1 ζ , and ω b 1 ζ . The continuity conditions (displacement and resultant force continuous) must satisfy along the interface L 1 . However, C 1 m and ϕ b 1 ζ are mutually dependent functions, which means the correction term C 1 m cannot be determined directly. Hence, the trial-and-error method is used to obtain the correction term. A discriminant equation is introduced initially and expressed as follows;
f 1 = C 1 m   f C 1 m   i , ( 33 ) C 1 m i = C 1 m   f ,   when   f 1 > 0.0001   ( iteration   continues ) , ( 34 ) C 1 m = C 1 m   i ,   when   f 1 < 0.0001   ( iteration   stops ) , ( 35 )
where C 1 m   i and C 1 m f respectively represent the initial and final values of the correction term C 1 m , and f 1 is the absolute difference between C 1 m   i and C 1 m f . The detailed iterative process is as follows:
First, the term ϕ b 1 ζ can be derived by substituting C 1 m   i = 1 into Equation (19). Next, C 1 m f can be obtained by substituting ϕ b 1 ζ into Equation (27). After the terms of C 1 m   i and C 1 m f are determined, the difference between C 1 m f and C 1 m   i could be calculated as displayed in Equation (33). Based on the while-loop, the iteration continued or stopped when f 1 is larger or less than 0.0001, respectively. The value C 1 m f replaces the value C 1 m   i when the difference is larger than 0.0001, as shown in Equation (34). Then, the new term ϕ b 1 ζ could be derived by substituting the replaced term C 1 m   i into Equation (19) again. Simultaneously, C 1 m f could be yielded by substituting ϕ b 1 ζ . After that, the status of the iteration depends on the difference between C 1 m   i and C 1 m f as displayed in Equation (33). Conversely, the iteration could be forcibly terminated because of the convergent value for f 1 when the difference value is less than 0.0001, as shown in Equation (35). Herein, the latest value of C 1 m   i would be defined as the approximate correction term C 1 m . Finally, the actual stress function ϕ b 1 ζ could be solved when the approximate correction term C 1 m is confirmed.
Because of the traction-free condition, the additional terms ϕ 1 ζ and ω 1 ζ  must be satisfied along the interface L 2 . L 2 . The correction term C 1 m * for the interface L 2 could be determined directly because C 1 m * involves ϕ b 1 ζ as shown in Equation (24), where ϕ b 1 ζ is yielded by the differential of the latest obtained term ϕ b 1 ζ in the first iteration for the interface L 1 . The approximate correction term C 1 m * is derived by substituting the term ϕ b 1 ζ into Equation (24). Lastly, the actual stress function ϕ 1 ζ could be solved by substituting the approximate correction term C 1 m * into Equation (22).
The second iteration for the interface L 1 is considered after the stress function ϕ 1 ζ is determined. The continuity conditions of displacement and resultant force must be satisfied by the stress functions ϕ a 2 ζ ,   ω a 2 ζ ,   ϕ b 2 ζ , and ω b 2 ζ along the interface L 1 . A discriminant equation of the second iteration for the interface L 1 is represented as follows:
f 2 = C 2 m f C 2 m i , ( 36 ) C 2 m i = C 2 m f ,   when   f 2 > 0.0001   ( iteration   continues ) , ( 37 ) C 2 m = C 2 m i ,   when   f 2 < 0.0001   ( iteration   stops ) , ( 38 )
where C 2 m   i and C 2 m f respectively represent the initial and final values of the correction term C 2 m , and f 2 is the absolute difference between C 2 m   i and C 2 m f . The detailed iterative process is as follows:
Similar to the first iteration, the term ϕ b 2 ζ could be calculated by substituting C 2 m   i = 1 into Equation (25). Meanwhile, C 2 m f can be yielded by substituting ϕ b 2 ζ into Equation (28). Next, f 2 is calculated after the terms of C 2 m   i and C 2 m f are determined, as displayed in Equation (36). The iteration continues or stops when f 2 is larger or less than 0.0001, respectively, as described in Equations (37) and (38). Then, the approximate correction term C 2 m is derived when the iteration stops. After that, the actual stress function ϕ 2 ζ for the interface L 2 could be determined by substituting the approximate correction term C 2 m * , which is yielded by substituting the latest obtained term of ϕ b 2 ζ . Finally, the general rules of the iteration for the polygonal hole with a coating layer could be summarized based on the above skill until n-times iteration.
A general form of the discriminant equation for the n-times iteration is written as:
f n = C n m f C n m i , ( 39 ) C n m i = C n m f ,   when   f n > 0.0001   ( iteration   continues ) ,   n = 1 , 2 , ( 40 ) C n m = C n m i ,   when   f n < 0.0001   ( iteration   stops ) , ( 41 )
where f n is the difference value; C n m i and C n m f are the initial and final values of the correction terms C n m in the corresponding situation, respectively.
By following the rules mentioned in this section, a series of correction terms could be determined.

4. Numerical Results

4.1. Temperature Distribution

The general form of temperature potential of a coated polygonal hole subjected to point thermal singularity is expressed in Equation (13). The temperature field for a coated triangle, square, and pentagonal hole with r = 1.1 under a point thermal singularity located at (3,0) and (−3,0) is considered in this study. Figure 5, Figure 6 and Figure 7a,b show the isothermal contours for a coated equiangular polygonal hole with different ratios of thermal conductivity. Furthermore, Figure 5, Figure 6 and Figure 7c show the quantitative analysis of a coated polygonal hole with different thermal conductivities.
Figure 5, Figure 6 and Figure 7a,b show the temperature contour distribution for different ratios of thermal conductivity. It can be pointed out that the isothermal contours have different changes in different ratios of thermal conductivity with k b / k a = 10 and k b / k a = 100 . When the thermal conductivity of the coating layer is increased, the contour lines around the coated hole are more dispersed. On the contrary, when the thermal conductivity of the coating layer is decreased, the contour lines around the coated hole are concentrated. It means that the higher thermal conductivity of the coating layer has a better heat dissipation effect. That is, a higher thermal conductivity can avoid a larger temperature gradient and prevent damage to the material at high temperatures. This phenomenon can be observed in the case of k b / k a = 100 , which can be seen from Figure 5, Figure 6 and Figure 7b.
The previously proposed isothermal contours belong to qualitative analysis that is based on personal visual perception. To obtain more accurate temperature results, the following arguments are explored with appropriate quantitative analysis. To discuss the temperature evolution on the x-axis, Figure 5, Figure 6 and Figure 7c provide a quantitative analysis of a coated equiangular triangle, square, and pentagonal hole with different thermal conductivities. For the limiting case of a polygonal hole, k b / k a = 1 , the heat flow is always along the adiabatic hole. It can be discovered that there exists a larger gradient due to which the temperature increases and decreases rapidly around the coating layer. This is because the theoretically adiabatic condition is assumed around the hole. Hence, the temperature is concentrated on the coating layer ((1.1,0) and (−1.1,0)), indicating that temperature approaches zero close to the hole due to adiabatic conditions. Furthermore, it can be evidently observed from the enlarged view in the upper left corner of Figure 5, Figure 6 and Figure 7c that when k b / k a = 1 , the temperature around the hole is 0.24; when k b / k a = 10 , the temperature drops to 0.12, and when k b / k a = 100 , the temperature drops further to near zero. It can be realized that, when the thermal conductivity rises by a factor of 10, the temperature can be reduced by half. The reason for this phenomenon is that the coating layer with a high thermal conductivity has the effect of heat absorption. This shows that when thermal conductivity is larger, the effect of endothermic is more obvious. The same trend can be seen from the different shaped coated holes in Figure 5, Figure 6 and Figure 7c.
When the entire material containing a coating layer is in a high-temperature environment, the coating layer with a higher thermal conductivity could display a lower temperature gradient environment. Herein, it will not be damaged in harsh temperature conditions, and the coating layer can achieve the effect of protecting the irregularly shaped holes and substrate. Furthermore, materials with higher thermal conductivity should be chosen because they can resist thermal fatigue and also have good heat absorption [10].

4.2. Interfacial Stresses

The summarized general stress functions of a coated equiangular polygonal hole applied to a point thermal singularity are given in Equations (29) and (30). In this study, all the results of stress fields are based on a coated pentagonal hole; that is, λ = 4 in Equation (5). In this section, we discuss a parametric analysis for point thermal singularity, shear modulus, thermal expansion coefficient, and shape factor.
Figure 8 shows a parametric analysis of different strengths of point thermal singularity for the coated pentagonal hole. A point heat source or sink was assumed to be located at point (3,0) along the x-axis. As predicted, when the strength of a point thermal singularity is high, the interfacial stress is increased. In Figure 8, the two curves at the bottom represent point heat source. The point heat source causes a compressive normal stress at the interface because the coating layer expands at high temperatures and the substrate resists it. Conversely, the two curves at the top represent the point heat sink. Tensile normal stress appears at the interface because the coating layer shrinks at low temperatures and the substrate invades it. The interfacial normal stress of each corner of the pentagonal hole under different magnitudes of heat loading is shown in Table 1. The interfacial stress has a positive and negative value under heat sink and source, respectively. The interfacial stress reduces by half as the magnitude of the thermal loading decreases by half. This can provide a reasonable explanation of the void evolution in MMC during the cooling stage. Both Al matrix with carbon diamond and AlSi7 matrix with SiC particle composites [10] exhibit the formation of irregularly shaped voids along the interface at eutectic solidification. Small rounded pores were observed to have formed within the SiC particles owing to the shrinkage of aluminum at the eutectic solidification. Subsequently, the small pores coalesced into larger irregular pores during the cooling stage [10]. Larger voids are generated because higher interfacial tensile stress occurs between the coating and adjacent materials when the temperature drops.
The interfacial normal and shear stresses for a pentagonal coated hole subjected to a point heat source is respectively symmetric and anti-symmetric with respect to the x-axis, as shown in Figure 9, Figure 10 and Figure 11, and this trend is similar to those reported in a previous study [33]. The extreme values of interfacial normal stress occurred at the corners of a pentagonal coating layer (0°, 72°, 144°, 216°, and 288°). In contrast, the interfacial shear stress vanishes at the corners of a pentagonal coating layer (0°, 72°, 144°, 216°, and 288°). Table 2 presents the interfacial normal stress of each corner of the pentagonal hole with different parameters. Table 2 depicts that the interfacial normal stress exhibited a maximum value as the corner was 0°. The corner 0° of a coated pentagonal hole was closer to the point heat source, which implies that a larger temperature gradient exists. Therefore, the interfacial normal stress was larger. Conversely, interfacial normal stress was smaller when the corner was 144° and 216°, which are far away from the point heat source.
Figure 9 expresses the interfacial stresses with different shear modulus ratios of a coated pentagonal hole subjected to a point heat source. Evidently, the interfacial stresses increased with a larger shear modulus of the coating layer. Previous research [39] reported that different modulus ratios of the coatings cause “elastic mismatch” and act as an important factor in controlling failure modes. In fracture mechanics, the interfacial delamination results from an increase in driving force when the elastic mismatch increases. The effect of different thermal expansion coefficient ratios around a coated pentagonal hole is shown in Figure 10. Due to different thermal expansion coefficients of the coating layer and the substrate, a phenomenon known as “thermal mismatch” occurs during the cooling phase. In addition, the interface between the substrate and coating layer may easily fracture due to the increased tensile stress caused by the thermal mismatch. Furthermore, Table 2 and Figure 9 and Figure 10 present that the interfacial stress is dominated by the thermal expansion coefficient compared to the shear modulus. Figure 11 presents the interfacial stress around a coated pentagonal hole with different shape factors. The interfacial stress around the coating layer was observed to be more intensified when the shape factor was larger. When the shape factor is large, the coating layer becomes sharper, and serious stress concentration is likely to occur at the corners of the interface. For practical engineering applications, the debonding mechanism is likely to result in interfacial failure due to interfacial stress concentration at the corners of the composites [40]. In addition, for TRIP-steel of MMC, the point of ductile damage begins from the sharper corners of the interface of the austenite matrix and zirconia particle when the global strain increases. However, most of the zirconia particles were cracked, void coalescence began to appear, and the damage progressed faster when the stress reached a critical value [41].

5. Conclusions

A general solution for a polygonal hole with a coating layer located at an infinite substrate under point thermal singularity is investigated in this study. A rapidly convergent series solution of both temperature and stress functions is obtained in an elegant form for the composites problem of an irregularly shaped void, inclusion, and substrate on the premise of using the conformal mapping function, analytical continuation theorem, and alternation method. The correction terms, expressed as a general solution, are solved using the trial-and-error method in MATLAB. The general solution proposed in this study is applied to the case with any number of edges of the polygonal hole with the coating layer. As numerical examples, the distribution of temperature contours provides the best optimal design with high thermal conductivity of the coating layer, thereby resulting in better heat absorption capacity through quantitative analysis. The interfacial stress demonstrated a larger stress concentration that occurred at the corners of the coated polygonal pores under thermal loading. Moreover, a reasonable explanation is provided for the growth of irregular voids in MMC during the cooling stage.

Author Contributions

Conceptualization: C.-K.C. and Y.-L.L.; methodology, Y.-L.L. and C.-K.C.; software, Y.-L.L. and S.-C.T.; formal analysis, Y.-L.L. and S.-C.T.; investigation, Y.-L.L.; resources, S.-C.T.; data curation, Y.-L.L. and S.-C.T.; writing—original draft preparation, Y.-L.L. and S.-C.T.; writing—review and editing, C.-K.C. and S.-C.T.; visualization, Y.-L.L. and S.-C.T.; project administration, C.-K.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by the Ministry of Science and Technology of Taiwan under the grant MOST 106-2221-E-011-065-MY3.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available on request from the corresponding author after obtaining permission from the authorized individual.

Acknowledgments

The authors acknowledge the support of Markus Kirschner and Heiko Winderlich (IMF, Technische Universität Bergakademie Freiberg) help regarding the capture microstructure of TRIP-steel with Mg-PSZ.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Determination of M and N in Equations (19) and (20)

Appendix A.1. For the Case of a Triangular Hole Where λ = 2

The temperature function g a ζ can be divided into three parts:
g a ζ = g a i ζ + g a o ζ + A log ζ
where g a o ζ and g a i ζ are defined as
g a o ( ζ ) = Q R w ζ 2 l o g ζ 0 2 w ζ ζ 0 + V b a Q R ζ l o g r 2 ζ ζ 0 ¯ ζ l o g ζ 0 ¯ r 2 ζ 0 ¯ l o g r 2 ζ ζ 0 ¯ + V b a Q R w ζ 2 l o g r 2 ζ ζ 0 ¯ w ζ 0 2 ¯ r 4 l o g r 2 ζ ζ 0 ¯ + U a b U b a Q R n = 1 V a b n 1 ζ l o g r 2 2 n ζ ζ 0 ¯ + w ζ 2 l o g r 2 2 n ζ ζ 0 ¯ U a b U b a Q R n = 1 V a b n 1 r 2 2 n ζ 0 ¯ + w ζ 0 2 ¯ r 4 4 n l o g r 2 2 n ζ ζ 0 ¯ U a b U b a Q R ζ 1 V a b l o g ζ 0 ¯ V b a Q R w ζ 0 ¯ r 2 ζ + w 2 ζ 2 U a b U b a Q R n = 1 V a b n 1 w ζ 0 ¯ r 2 2 n ζ + w 2 ζ 2
g a i ( ζ ) = Q R ζ ζ 0 l o g ζ ζ 0 1 + w ζ 2 w ζ 0 2 l o g ζ ζ 0 w ζ 2 l o g ( ζ 0 ) + w ζ ζ 0 + V b a Q R ζ l o g ζ 0 ¯ + U a b U b a Q R ζ l o g ζ 0 ¯ 1 V a b
The constant A in (A1) is defined as
A = Q R w ζ 0 2 V b a r 2 ζ 0 ¯ U a b U b a r 2 ζ 0 ¯ r 2 V a b
The displacement can be expressed as
u x + i u y = 1 2 G a κ a ϕ ζ ω ζ ¯ + β a g a ζ
This leads to
u x + i u y c = 1 2 G a κ a ϕ ζ ω ζ ¯ + β a g a ζ c = 0
where f ζ c = f r , 2 π f r , 0 . The coefficient of the logarithmic term must vanish, which yields
1 + κ a 2 G a M + β a A log ζ = 0
Solving for M in Equation (A5) yields
M = 2 G a β a ( κ a + 1 ) Q R U a b U b a r 2 ζ 0 ¯ ( r 2 V a b ) + V b a r 2 ζ 0 ¯ w ζ 0 2
The temperature potential function g b ζ must be split into three parts:
g b ζ = g b i ζ + g b o ζ + B log ζ
where g b o ζ and g b i ζ are defined as
g b o ( ζ ) = U b a Q R w ζ 2 l o g ( ζ 0 ) 2 w ζ ζ 0 + V a b w l o g ζ 0 ζ 2 1 V a b 2 V a b w ζ ζ 0 r 2 V a b + U b a Q R n = 1 V a b n 1 ζ + w ζ 2 l o g r 2 2 n ζ ζ 0 ¯ U b a Q R n = 1 V a b n 1 r 2 2 n ζ 0 ¯ + w ζ 0 2 ¯ r 4 4 n l o g r 2 2 n ζ ζ 0 ¯ U b a Q R ζ l o g ζ 0 ¯ 1 V a b + w ζ 0 ¯ ζ 1 V a b r 2 + w 2 ζ 2 1 V a b
g b i ( ζ ) = U b a Q R ζ ζ 0 l o g ζ ζ 0 1 + U b a Q R w ζ 2 w ζ 0 2 l o g ζ ζ 0 w ζ 2 l o g ζ 0 + w ζ ζ 0 + U b a Q R ζ n = 1 V a b n l o g r 2 n ζ ζ 0 + U b a Q R w ζ 2 n = 1 V a b n l o g 1 ζ r 2 n ζ 0 + U b a V a b Q R w ζ ζ 0 r 2 V a b U b a Q R n = 1 V a b n r 2 n ζ 0 + w ζ 0 2 r 4 n l o g ζ r 2 n ζ 0 + U b a Q R ζ l o g ζ 0 ¯ 1 V a b V a b 1 V a b
And the constant B in (A7) is defined as
B = U b a Q R w ζ 0 2 + V a b w ζ 0 2 r 4 V a b r 2 ζ 0 ¯ r 2 V a b
The displacement can be expressed as
u x + i u y = 1 2 G b κ b ϕ ζ ω ζ ¯ + β b g b ζ
which leads to
1 + κ b 2 G b N + β b B log ζ = 0
Substituting (A8) into (A10), we obtain
N = 2 G b β b ( κ b + 1 ) U b a Q R 2 r 2 ζ 0 ¯ ( r 2 V a b ) V a b w ζ 0 3 ( r 6 V a b ) w ζ 0 3

Appendix A.2. For the Case of a Square Hole Where λ = 3

The temperature function of g a i ζ ,   g a o ζ ,   g b i ζ ,   g b o ζ are shown below
g a o ( ζ ) = Q R 2 w ζ 3 l o g ζ 0 w ζ ζ 0 3 2 ζ 0 + 3 2 ζ + V b a Q R 2 ζ l o g r 2 ζ ζ 0 ¯ ζ l o g ζ 0 ¯ r 2 ζ 0 ¯ l o g r 2 ζ ζ 0 ¯ + V b a Q R 2 w ζ 3 w ζ 0 3 ¯ r 6 l o g r 2 ζ ζ 0 ¯ + U a b U b a Q R 2 n = 1 V a b n 1 ζ + w ζ 3 l o g r 2 2 n ζ ζ 0 ¯ U a b U b a Q R 2 n = 1 V a b n 1 r 2 2 n ζ 0 ¯ + w ζ 0 3 ¯ r 6 6 n l o g r 2 2 n ζ ζ 0 ¯ U a b U b a Q R 2 ζ 1 V a b l o g ζ 0 ¯ V b a Q R 2 w ζ 0 ¯ 2 ζ 2 r 2 + w 3 ζ 3 + w ζ 0 2 ¯ r 4 ζ U a b U b a Q R 2 n = 1 V a b n 1 w ζ 0 ¯ 2 r 2 2 n ζ 2 + w 3 ζ 3 + w ζ 0 2 ¯ r 4 4 n ζ
g a i ( ζ ) = Q R 2 ζ ζ 0 l o g ζ ζ 0 1 + w ζ 3 w ζ 0 3 l o g ζ ζ 0 w ζ 3 l o g ( ζ 0 ) + w ζ ζ 0 1 2 ζ 0 + 1 ζ + V b a Q R 2 ζ l o g ζ 0 ¯ + U a b U b a Q R 2 ζ l o g ζ 0 ¯ 1 V a b
g b o ( ζ ) = U b a Q R 2 w ζ 3 l o g ( ζ 0 ) 3 w ζ ζ 0 1 2 ζ 0 + 1 2 ζ + V a b w l o g ζ 0 ζ 3 1 V a b 3 V a b w 2 ζ 2 ζ 0 r 2 V a b 3 V a b w 2 ζ ζ 0 2 r 4 V a b + U b a Q R 2 n = 1 V a b n 1 ζ + w ζ 3 l o g r 2 2 n ζ ζ 0 ¯ U b a Q R 2 n = 1 V a b n 1 r 2 2 n ζ 0 ¯ + w ζ 0 3 ¯ r 6 6 n l o g r 2 2 n ζ ζ 0 ¯ U b a Q R 2 ζ l o g ζ 0 ¯ 1 V a b + w ζ 0 ¯ 2 ζ 2 1 V a b r 2 + w ζ 0 2 ¯ ζ 1 V a b r 4 + w 3 ζ 3 1 V a b
g b i ( ζ ) = U b a Q R 2 ζ ζ 0 l o g ζ ζ 0 1 + U b a Q R 2 w ζ 3 w ζ 0 3 l o g ζ ζ 0 w ζ 3 l o g ζ 0 + w ζ ζ 0 1 ζ + 1 2 ζ 0 + U b a Q R 2 ζ n = 1 V a b n l o g ζ r 2 n ζ 0 + U b a Q R 2 w ζ 3 n = 1 V a b n l o g 1 ζ r 2 n ζ 0 + U b a V a b Q R 2 w ζ 2 ζ 0 r 2 V a b + w 2 ζ ζ 0 2 r 4 V a b U b a Q R 2 n = 1 V a b n r 2 n ζ 0 + w ζ 0 3 r 6 n l o g ζ r 2 n ζ 0 + U b a Q R 2 ζ l o g ζ 0 ¯ 1 V a b V a b 1 V a b
The constants M and N for λ = 3 yield
M = 2 G a β a κ a + 1 Q R 2 U a b U b a r 2 ζ 0 ¯ r 2 V a b + V b a r 2 ζ 0 ¯ w ζ 0 3
N = 2 G b β b κ b + 1 U b a Q R 2 r 2 ζ 0 ¯ r 2 V a b V a b w ζ 0 3 r 6 V a b w ζ 0 3

Appendix A.3. For the Case of a Pentagonal Hole Where λ = 4

The temperature function of g a i ζ ,   g a o ζ ,   g b i ζ ,   g b o ζ are shown below
g a o ( ζ ) = Q R 3 w ζ 4 l o g ζ 0 w ζ ζ 0 4 3 ζ 0 2 + 1 ζ ζ 0 + 4 3 ζ 2 + V b a Q R 3 ζ l o g r 2 ζ ζ 0 ¯ ζ l o g ζ 0 ¯ r 2 ζ 0 ¯ l o g r 2 ζ ζ 0 ¯ + V b a Q R 3 w ζ 4 w ζ 0 4 ¯ r 8 l o g r 2 ζ ζ 0 ¯ + U a b U b a Q R 3 n = 1 V a b n 1 ζ + w ζ 4 l o g r 2 2 n ζ ζ 0 ¯ U a b U b a Q R 3 n = 1 V a b n 1 r 2 2 n ζ 0 ¯ + w ζ 0 4 ¯ r 8 8 n l o g r 2 2 n ζ ζ 0 ¯ U a b U b a Q R 3 ζ 1 V a b l o g ζ 0 ¯ V b a Q R 3 w 4 ζ 4 + w ζ 0 ¯ 3 r 2 ζ 3 + w ζ 0 2 ¯ 2 r 4 ζ 2 + w ζ 0 3 ¯ r 6 ζ U a b U b a Q R 3 n = 1 V a b n 1 w 4 ζ 4 + w ζ 0 ¯ 3 r 2 2 n ζ 3 + w ζ 0 2 ¯ 2 r 4 4 n ζ 2 + w ζ 0 3 ¯ r 6 6 n ζ
g a i ( ζ ) = Q R 3 ζ ζ 0 l o g ζ ζ 0 1 + w ζ 4 w ζ 0 4 l o g ζ ζ 0 w ζ 4 l o g ( ζ 0 ) + w ζ ζ 0 1 ζ 2 + 1 2 ζ ζ 0 + 1 3 ζ 0 2 + V b a Q R 3 ζ l o g ζ 0 ¯ + U a b U b a Q R 3 ζ l o g ζ 0 ¯ 1 V a b
g b o ( ζ ) = U b a Q R 3 w ζ 4 l o g ( ζ 0 ) w ζ ζ 0 4 3 ζ 2 + 1 ζ ζ 0 4 3 ζ 0 2 + V a b w l o g ζ 0 ζ 4 1 V a b 4 V a b w 3 ζ 3 ζ 0 r 2 V a b V a b w ζ 2 ζ 0 2 r 4 V a b 4 V a b w 3 ζ ζ 0 3 r 6 V a b + U b a Q R 3 n = 1 V a b n 1 ζ + w ζ 4 l o g r 2 2 n ζ ζ 0 ¯ U b a Q R 3 n = 1 V a b n 1 r 2 2 n ζ 0 ¯ + w ζ 0 4 ¯ r 8 8 n l o g r 2 2 n ζ ζ 0 ¯ U b a Q R 3 ζ l o g ζ 0 ¯ 1 V a b + w 4 ζ 4 1 V a b + w ζ 0 ¯ 3 ζ 3 1 V a b r 2 + w ζ 0 2 ¯ 2 ζ 2 1 V a b r 4 + w ζ 0 3 ¯ ζ 1 V a b r 6
g b i ( ζ ) = U b a Q R 3 ζ ζ 0 l o g ζ ζ 0 1 + U b a Q R 3 w ζ 4 w ζ 0 4 l o g ζ ζ 0 w ζ 4 l o g ζ 0 + w ζ ζ 0 1 ζ 2 + 1 2 ζ ζ 0 + 1 3 ζ 0 2 + U b a Q R 3 ζ n = 1 V a b n l o g ζ r 2 n ζ 0 + U b a Q R 3 w ζ 4 n = 1 V a b n l o g 1 ζ r 2 n ζ 0 + U b a V a b Q R 3 w ζ 3 ζ 0 r 2 V a b + w 2 ζ 2 ζ 0 2 r 4 V a b + w 3 ζ ζ 0 3 r 6 V a b U b a Q R 3 n = 1 V a b n r 2 n ζ 0 + w ζ 0 4 r 8 n l o g ζ r 2 n ζ 0 + U b a Q R 3 ζ l o g ζ 0 ¯ 1 V a b V a b 1 V a b
The constants M and N for λ = 4 yield
M = 2 G a β a κ a + 1 Q R 3 U a b U b a r 2 ζ 0 ¯ r 2 V a b + V b a r 2 ζ 0 ¯ w ζ 0 4
N = 2 G b β b κ b + 1 U b a Q R 3 r 2 ζ 0 ¯ r 2 V a b V a b w ζ 0 4 r 8 V a b w ζ 0 4

References

  1. Zweben, C. Metal-matrix composites for electronic packaging. J. Miner. Met. Mater. Soc. 1992, 44, 15–23. [Google Scholar] [CrossRef]
  2. Beffort, O.; Vaucher, S.; Khalid, F. On the thermal and chemical stability of diamond during processing of Al/diamond composites by liquid metal infiltration (squeeze casting). Diam. Relat. Mater. 2004, 13, 1834–1843. [Google Scholar] [CrossRef]
  3. Jagadeesh, G.V.; Setti, S.G. A review on micromechanical methods for evaluation of mechanical behavior of particulate reinforced metal matrix composites. J. Mater. Sci. 2020, 55, 9848–9882. [Google Scholar] [CrossRef]
  4. Hussein, T.; Umar, M.; Qayyum, F.; Guk, S.; Prahl, U. Micromechanical Effect of Martensite Attributes on Forming Limits of Dual-Phase Steels Investigated by Crystal Plasticity-Based Numerical Simulations. Crystals 2022, 12, 155. [Google Scholar] [CrossRef]
  5. Qayyum, F.; Umar, M.; Guk, S.; Schmidtchen, M.; Kawalla, R.; Prahl, U. Effect of the 3rd Dimension within the Representative Volume Element (RVE) on Damage Initiation and Propagation during Full-Phase Numerical Simulations of Single and Multi-Phase Steels. Materials 2020, 14, 42. [Google Scholar] [CrossRef]
  6. Schöbel, M.; Degischer, H.; Vaucher, S.; Hofmann, M.; Cloetens, P. Reinforcement architectures and thermal fatigue in diamond particle-reinforced aluminum. Acta Mater. 2010, 58, 6421–6430. [Google Scholar] [CrossRef]
  7. Nam, T.H.; Requena, G.; Degischer, P. Thermal expansion behaviour of aluminum matrix composites with densely packed SiC particles. Compos. Part A Appl. Sci. Manuf. 2008, 39, 856–865. [Google Scholar] [CrossRef]
  8. Beffort, O.; Khalid, F.; Weber, L.; Ruch, P.; Klotz, U.; Meier, S.; Kleiner, S. Interface formation in infiltrated Al(Si)/diamond composites. Diam. Relat. Mater. 2006, 15, 1250–1260. [Google Scholar] [CrossRef]
  9. Schöbel, M.; Altendorfer, W.; Degischer, H.; Vaucher, S.; Buslaps, T.; Di Michiel, M.; Hofmann, M. Internal stresses and voids in SiC particle reinforced aluminum composites for heat sink applications. Compos. Sci. Technol. 2011, 71, 724–733. [Google Scholar] [CrossRef]
  10. Schöbel, M.; Requena, G.; Fiedler, G.; Tolnai, D.; Vaucher, S.; Degischer, H. Void formation in metal matrix composites by solidification and shrinkage of an AlSi7 matrix between densely packed particles. Compos. Part A Appl. Sci. Manuf. 2014, 66, 103–108. [Google Scholar] [CrossRef]
  11. Rivera-Salinas, J.E.; Gregorio-Jáuregui, K.M.; Romero-Serrano, J.A.; Cruz-Ramírez, A.; Hernández-Hernández, E.; Miranda-Pérez, A.; Gutierréz-Pérez, V.H. Simulation on the Effect of Porosity in the Elastic Modulus of SiC Particle Reinforced Al Matrix Composites. Metals 2020, 10, 391. [Google Scholar] [CrossRef]
  12. Evans, A.G.; Mumm, D.R.; Hutchinson, J.W.; Meier, G.H.; Pettit, F.S. Mechanisms controlling the durability of thermal barrier coatings. Prog. Mater. Sci. 2001, 46, 505–553. [Google Scholar] [CrossRef]
  13. Ramakrishnan, N.; Arunachalam, V.S. Effective elastic moduli of porous solids. J. Mater. Sci. 1990, 25, 3930–3937. [Google Scholar] [CrossRef]
  14. Ramakrishnan, N.; Arunachalam, V.S. Effective Elastic Moduli of Porous Ceramic Materials. J. Am. Ceram. Soc. 1993, 76, 2745–2752. [Google Scholar] [CrossRef]
  15. Muskhelishvili, N.I. Some Basic Problems of the Mathematical Theory of Elasticity; Noordhoff: Groningen, The Netherlands, 1953; Volume 15. [Google Scholar] [CrossRef]
  16. Chao, C.K.; Chen, F.M.; Shen, M.H. Green’s Functions for a Point Heat Source in Circularly Cylindrical Layered Media. J. Therm. Stress. 2006, 29, 809–847. [Google Scholar] [CrossRef]
  17. Chao, C.; Chen, F.; Lin, T. Green’s function for a point heat source embedded in an infinite body with two circular elastic inclusions. Appl. Math. Model. 2018, 56, 254–274. [Google Scholar] [CrossRef]
  18. Chao, C.K.; Chen, C.K.; Chen, F.M. Interfacial stresses induced by a point heat source in an isotropic plate with a reinforced elliptical hole. Comput. Model. Eng. Sci. (CMES) 2010, 63, 1. [Google Scholar] [CrossRef]
  19. Zenkour, A.M.; Mashat, D.S.; Allehaibi, A.M. Thermoelastic Coupling Response of an Unbounded Solid with a Cylindrical Cavity Due to a Moving Heat Source. Mathematics 2022, 10, 9. [Google Scholar] [CrossRef]
  20. Attia, M.A.; Melaibari, A.; Shanab, R.A.; Eltaher, M.A. Dynamic Analysis of Sigmoid Bidirectional FG Microbeams under Moving Load and Thermal Load: Analytical Laplace Solution. Mathematics 2022, 10, 4797. [Google Scholar] [CrossRef]
  21. Chiu, C.; Tseng, S.; Qayyum, F.; Guk, S.; Chao, C.; Prahl, U. Local deformation and interfacial damage behavior of partially stabilized zirconia-reinforced metastable austenitic steel composites: Numerical simulation and validation. Mater. Des. 2022, 225, 111515. [Google Scholar] [CrossRef]
  22. Tseng, S.-C.; Chiu, C.-C.; Qayyum, F.; Guk, S.; Chao, C.-K.; Prahl, U. The Effect of the Energy Release Rate on the Local Damage Evolution in TRIP Steel Composite Reinforced with Zirconia Particles. Materials 2023, 16, 134. [Google Scholar] [CrossRef] [PubMed]
  23. Xie, C.; Wu, Y.; Liu, Z. Stress fields and effective modulus of piezoelectric fiber composite with arbitrary shaped inclusion under in-plane mechanical and anti-plane electric loadings. Math. Mech. Solids 2019, 24, 3180–3199. [Google Scholar] [CrossRef]
  24. Luo, J.-C.; Gao, C.-F. Stress field of a coated arbitrary shape inclusion. Meccanica 2011, 46, 1055–1071. [Google Scholar] [CrossRef]
  25. Yoshikawa, K.; Hasebe, N. Green’s function of the displacement boundary value problem for a heat source in an infinite plane with an arbitrary shaped rigid inclusion. Arch. Appl. Mech. 1999, 69, 227–239. [Google Scholar] [CrossRef]
  26. Tang, J.-Y.; Yang, H.-B. An alternative numerical scheme for calculating the thermal stresses around an inclusion of arbitrary shape in an elastic plane under uniform remote in-plane heat flux. Acta Mech. 2019, 230, 2399–2412. [Google Scholar] [CrossRef]
  27. Yang, Q.; Zhu, W.; Li, Y.; Zhang, H. Stress field of a functionally graded coated inclusion of arbitrary shape. Acta Mech. 2018, 229, 1687–1701. [Google Scholar] [CrossRef]
  28. Zhou, X.; Shao, Z.; Tian, F.; Hopper, C.; Jiang, J. Microstructural effects on central crack formation in hot cross-wedge-rolled high-strength steel parts. J. Mater. Sci. 2020, 55, 9608–9622. [Google Scholar] [CrossRef]
  29. Tseng, S.C.; Chao, C.K.; Chen, F.M. Interfacial Stresses of a Coated Square Hole Induced by a Remote Uniform Heat Flow. Int. J. Appl. Mech. 2020, 12, 2050063. [Google Scholar] [CrossRef]
  30. Tseng, S.C.; Chao, C.K.; Guo, J.Y. Failure Analysis of a Polygonal Void with an Oxide Layer in a Cracked Matrix. Int. J. Appl. Mech. 2021, 13, 2150099. [Google Scholar] [CrossRef]
  31. Chiu, C.; Tseng, S.; Chao, C.; Guo, J. Stress Intensity Factors for a Non-Circular Hole with Inclusion Layer Embedded in a Cracked Matrix. Aerospace 2021, 9, 17. [Google Scholar] [CrossRef]
  32. Liao, Y.L.; Tseng, S.C.; Chao, C.K. Stress analysis of an inclusion layer bonded to an irregularly shaped pore under an edge dislocation or a concentrated load. J. Mech. 2022, 38, 397–409. [Google Scholar] [CrossRef]
  33. Tseng, S.C.; Chao, C.K.; Chen, F.M.; Chiu, W.C. Interfacial stresses of a coated polygonal hole subject to a point heat source. J. Therm. Stress. 2020, 43, 1487–1512. [Google Scholar] [CrossRef]
  34. Xue, X.-Y.; Wen, S.-R.; Sun, J.-Y.; He, X.-T. One- and Two-Dimensional Analytical Solutions of Thermal Stress for Bimodular Functionally Graded Beams under Arbitrary Temperature Rise Modes. Mathematics 2022, 10, 1756. [Google Scholar] [CrossRef]
  35. He, X.-T.; Zhang, M.-Q.; Pang, B.; Sun, J.-Y. Solution of the Thermoelastic Problem for a Two-Dimensional Curved Beam with Bimodular Effects. Mathematics 2022, 10, 3002. [Google Scholar] [CrossRef]
  36. Bogdanoff, J.L. Note on Thermal Stresses. J. Appl. Mech. 1954, 21, 88. [Google Scholar] [CrossRef]
  37. Qayyum, F.; Chaudhry, A.A.; Guk, S.; Schmidtchen, M.; Kawalla, R.; Prahl, U. Effect of 3D Representative Volume Element (RVE) Thickness on Stress and Strain Partitioning in Crystal Plasticity Simulations of Multi-Phase Materials. Crystals 2020, 10, 944. [Google Scholar] [CrossRef]
  38. Wu, C.; Wang, D.; Guo, Y. An immersed particle modeling technique for the three-dimensional large strain simulation of particulate-reinforced metal-matrix composites. Appl. Math. Model. 2016, 40, 2500–2513. [Google Scholar] [CrossRef]
  39. Jiang, P.; Fan, X.; Sun, Y.; Li, D.; Li, B.; Wang, T. Competition mechanism of interfacial cracks in thermal barrier coating system. Mater. Des. 2017, 132, 559–566. [Google Scholar] [CrossRef]
  40. Evans, A.G.; He, M.Y.; Hutchinson, J.W. Mechanics-based scaling laws for the durability of thermal barrier coatings. Prog. Mater. Sci. 2001, 46, 249–271. [Google Scholar] [CrossRef]
  41. Qayyum, F.; Guk, S.; Prahl, U. Studying the Damage Evolution and the Micro-Mechanical Response of X8CrMnNi16-6-6 TRIP Steel Matrix and 10% Zirconia Particle Composite Using a Calibrated Physics and Crystal-Plasticity-Based Numerical Simulation Model. Crystals 2021, 11, 759. [Google Scholar] [CrossRef]
Figure 1. (a) Triangle-shaped inclusion, (b) quadrilateral-shaped inclusion, (c) polygonal-shaped inclusion [21,22], and (d) macrostructure of the hole, iron oxide, and steel matrix [28].
Figure 1. (a) Triangle-shaped inclusion, (b) quadrilateral-shaped inclusion, (c) polygonal-shaped inclusion [21,22], and (d) macrostructure of the hole, iron oxide, and steel matrix [28].
Mathematics 11 00802 g001
Figure 2. The (a) physical plane (z-plane) and (b) mathematical plane (ζ-plane) of a polygonal coated hole under thermal singularity.
Figure 2. The (a) physical plane (z-plane) and (b) mathematical plane (ζ-plane) of a polygonal coated hole under thermal singularity.
Mathematics 11 00802 g002
Figure 3. A flow diagram of the current study including theoretical method, temperature field, stress field, and result.
Figure 3. A flow diagram of the current study including theoretical method, temperature field, stress field, and result.
Mathematics 11 00802 g003
Figure 4. The different shape factors w of a pentagonal coated hole in the physical plane, where (a) w = 0.07 and (b) w = 0.1.
Figure 4. The different shape factors w of a pentagonal coated hole in the physical plane, where (a) w = 0.07 and (b) w = 0.1.
Mathematics 11 00802 g004
Figure 5. Contours of temperature distribution for a coated equiangular triangle hole, where (a) kb/ka = 10 and (b) kb/ka = 100. (c) Quantitatively normalized temperature for a coated triangular hole with different ratios of heat conductivity, where w = 0.1.
Figure 5. Contours of temperature distribution for a coated equiangular triangle hole, where (a) kb/ka = 10 and (b) kb/ka = 100. (c) Quantitatively normalized temperature for a coated triangular hole with different ratios of heat conductivity, where w = 0.1.
Mathematics 11 00802 g005
Figure 6. Contours of temperature distribution for a coated equiangular square hole, where (a) kb/ka = 10 and (b) kb/ka = 100. (c) Quantitatively normalized temperature for a coated triangular hole with different ratios of heat conductivity, where w = 0.1.
Figure 6. Contours of temperature distribution for a coated equiangular square hole, where (a) kb/ka = 10 and (b) kb/ka = 100. (c) Quantitatively normalized temperature for a coated triangular hole with different ratios of heat conductivity, where w = 0.1.
Mathematics 11 00802 g006
Figure 7. Contours of temperature distribution for a coated equiangular pentagonal hole, where (a) kb/ka = 10 and (b) kb/ka = 100. (c) Quantitatively normalized temperature for a coated triangular hole with different ratios of heat conductivity, where w = 0.1.
Figure 7. Contours of temperature distribution for a coated equiangular pentagonal hole, where (a) kb/ka = 10 and (b) kb/ka = 100. (c) Quantitatively normalized temperature for a coated triangular hole with different ratios of heat conductivity, where w = 0.1.
Mathematics 11 00802 g007
Figure 8. Interfacial (a) normal and (b) shear stress of a coated pentagonal hole (w = 0.05, Gb/Ga = 3, βbb = 3, kb/ka = 2, νb = 0.3, νa = 0.3).
Figure 8. Interfacial (a) normal and (b) shear stress of a coated pentagonal hole (w = 0.05, Gb/Ga = 3, βbb = 3, kb/ka = 2, νb = 0.3, νa = 0.3).
Mathematics 11 00802 g008
Figure 9. Interfacial (a) normal and (b) shear stress of a coated pentagonal hole (w = 0.05, βb/βb = 3, kb/ka = 2, νb = 0.3, νa = 0.3).
Figure 9. Interfacial (a) normal and (b) shear stress of a coated pentagonal hole (w = 0.05, βb/βb = 3, kb/ka = 2, νb = 0.3, νa = 0.3).
Mathematics 11 00802 g009
Figure 10. Interfacial (a) normal and (b) shear stress of a coated pentagonal hole (w = 0.05, Gb/Ga = 3, kb/ka = 2, νb = 0.3, νa = 0.3).
Figure 10. Interfacial (a) normal and (b) shear stress of a coated pentagonal hole (w = 0.05, Gb/Ga = 3, kb/ka = 2, νb = 0.3, νa = 0.3).
Mathematics 11 00802 g010
Figure 11. Interfacial (a) normal and (b) shear stress of a coated pentagonal hole (Gb/Ga = 3, βb/βb = 3, kb/ka = 2, νb = 0.3, νa = 0.3).
Figure 11. Interfacial (a) normal and (b) shear stress of a coated pentagonal hole (Gb/Ga = 3, βb/βb = 3, kb/ka = 2, νb = 0.3, νa = 0.3).
Mathematics 11 00802 g011
Table 1. The interfacial normal stress of each corner of the pentagonal hole under different magnitudes of heat loading.
Table 1. The interfacial normal stress of each corner of the pentagonal hole under different magnitudes of heat loading.
MagnitudeDegree
72°144°216°288°
Q = −q0/2πka−3.0−2.4−1.8−1.8−2.4
Q = −q0/4πka−1.5−1.2−0.9−0.9−1.2
Q = q0/4πka1.51.20.90.91.2
Q = q0/2πka3.02.41.81.82.4
Table 2. The interfacial normal stress of each corner of the pentagonal hole with different parameters.
Table 2. The interfacial normal stress of each corner of the pentagonal hole with different parameters.
ParametersDegree
72°144°216°288°
Ga/Gb = 3−3.0−2.4−1.8−1.8−2.4
Ga/Gb = 4−4.7−3.7−2.9−2.9−3.7
Ga/Gb = 5−6.5−5.2−4.2−4.2−5.2
βa/βb = 3−3.0−2.4−1.8−1.8−2.4
βa/βb = 4−5.3−4.8−3.5−3.5−4.8
βa/βb = 5−7.6−6.8−5.2−5.2−6.8
w = 0.04−2.5−2.1−1.8−1.8−2.1
w = 0.05−3.0−2.4−2.0−2.0−2.4
w = 0.06−3.8−3.0−2.4−2.4−3.0
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Liao, Y.-L.; Tseng, S.-C.; Chao, C.-K. Interfacial Stresses for a Coated Irregularly Shaped Hole Embedded in an Infinite Solid under Point Heat Singularity. Mathematics 2023, 11, 802. https://doi.org/10.3390/math11040802

AMA Style

Liao Y-L, Tseng S-C, Chao C-K. Interfacial Stresses for a Coated Irregularly Shaped Hole Embedded in an Infinite Solid under Point Heat Singularity. Mathematics. 2023; 11(4):802. https://doi.org/10.3390/math11040802

Chicago/Turabian Style

Liao, Yi-Lun, Shao-Chen Tseng, and Ching-Kong Chao. 2023. "Interfacial Stresses for a Coated Irregularly Shaped Hole Embedded in an Infinite Solid under Point Heat Singularity" Mathematics 11, no. 4: 802. https://doi.org/10.3390/math11040802

APA Style

Liao, Y. -L., Tseng, S. -C., & Chao, C. -K. (2023). Interfacial Stresses for a Coated Irregularly Shaped Hole Embedded in an Infinite Solid under Point Heat Singularity. Mathematics, 11(4), 802. https://doi.org/10.3390/math11040802

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