Next Article in Journal
Fatigue Performance Prediction of RC Beams Based on Optimized Machine Learning Technology
Next Article in Special Issue
Residual Stress Redistribution Analysis in the Repair Welding of AA6082-T6 Aluminum Alloy Joints: Experiment and Simulation
Previous Article in Journal
A Fukui Analysis of an Arginine-Modified Carbon Surface for the Electrochemical Sensing of Dopamine
Previous Article in Special Issue
The Effect of Nucleating Agents on Enzyme-Induced Carbonate Precipitation and Corresponding Microscopic Mechanisms
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Magneto-Thermo-Elastic Theoretical Solution for Functionally Graded Thick-Walled Tube under Magnetic, Thermal and Mechanical Loads Based on Voigt Method

1
College of Civil Engineering and Architecture, Zhejiang University, Hangzhou 310058, China
2
College of Civil Engineering, Taiyuan University of Technology, Taiyuan 030024, China
3
Zhejiang Provincial Key Laboratory of Space Structures, Hangzhou 310058, China
*
Author to whom correspondence should be addressed.
Materials 2022, 15(18), 6345; https://doi.org/10.3390/ma15186345
Submission received: 19 July 2022 / Revised: 15 August 2022 / Accepted: 7 September 2022 / Published: 13 September 2022
(This article belongs to the Special Issue Fatigue Behavior, Lifetime Prediction and Modeling of Welding Process)

Abstract

:
In this study, the mechanical responses of a functionally graded thick-walled tube simultaneously under magnetic, thermal and mechanical loads are studied. Based on the assumption that the volume fraction of each phase material is distributed as a power function, the Voigt method is used to obtain the stress–strain relationship of the functionally graded materials (FGMs). The influences of the relevant material parameters including volume fraction, thermal expansion coefficient, and Poisson’s ratio on the magneto-thermo-elastic theoretical solution are deeply studied and discussed. Furthermore, when some of the parameters are set as special values, the research results can be degenerated to two coupled loads which are consistent with the existing researches. The results of this paper provide theoretical support for the practical design and application of the FGM tube under the combined action of magnetic, thermal and mechanical loads.

1. Introduction

Functionally graded materials (FGMs) [1,2,3] are a type of newly developed composite material with continuous material and mechanical properties, which effectively eliminate stress concentrations in common laminates. FGMs can be designed in various shapes as needed, such as disc [4], plate [5,6,7], beam [8] and cylinder [9], etc. As a commonly used structure in engineering, functionally graded thick-walled tubes have attracted extensive attention and research in the materials industry due to their unique physical and mechanical properties, and have begun to be applied in many different practical engineering fields, such as aerospace, magnetic storage components, magnetic structural components, etc.
The magneto-thermo-elastic environment is a fairly common phenomenon in industrial applications, such as in nuclear devices [10], development of a highly sensitive superconducting magnetometer, electrical equipment [11] and optics, etc. [12]. In the last couple of decades, a number of researchers have investigated this problem.
The researches on theoretical solutions for the uniform circular tube with mechanical loads, thermal loads and magnetic fields are relatively mature [13,14,15]. In recent years, some scholars have studied theoretical solutions for the graded thick-walled tube, based on certain assumptions (see Table 1).
When the responses of the functionally graded tube under the combination of the above three loads are investigated, it is necessary to make assumptions about all the above parameters [16,17,18,19,20,21,22,23,24,25,26,27,28]. In past studies, Poisson’s ratio is set as a constant and most material parameters are assumed as power functions, such as elastic modulus E(r) = E0rn, thermal expansion coefficient α(r) = α0rn, thermal conductivity k(r) = k0rn, and magnetic permeability μ(r) = μ0rn. Some researchers also assume these material parameters as exponential functions [27], such as E(r) = E0ern, α(r) = α0em, k(r) = k0ern, μ(r) = μ0 ern (E0, α0, k0, μ0, m and n are the material constants, and r is the radial coordinate).
A summary of the existing literature finds that: (1) the indexes m or n are basically the same, but this is inconsistent with objective reality if all material parameters have the same form. Some authors [28] have proposed that the indexes should not be the same. However, it is also considered the same in ultimate numerical analysis because the problems are complicated by different indexes. (2) As is well known, the elastic modulus and Poisson’s ratio are basic material parameters of the uniform isotropic materials. Therefore, it is obviously more reasonable to simultaneously make a hypothesis of elastic modulus and Poisson’s ratio when assuming the entire material parameters for a functionally graded tube. However, most scholars think that Poisson’s ratio of each layer of the tube has a small effect on the responses of the whole tube, and set Poisson’s ratio as a constant to avoid the complexity of solving equilibrium equations. Only a few scholars think that Poisson’s ratio should also be considered as another material parameter [29,30,31]. If Poisson’s ratio is also assumed as a specific function, often along with complex equilibrium equations, it is difficult to obtain theoretical solutions. (3) From the aspect of the production of engineering materials, it is difficult to change material parameters according to a certain rule. The FGMs are achieved by combining two or more homogeneous materials in a controlled proportion; this does not conform to the material’s practical manufacturing method if one simply assumes these material parameters.
In order to overcome the above problems, the equivalent material parameters of the thick-walled tube are given out using the Voigt method, based on the volume fractions and respective material parameters of the two-phase materials that compose the tube instead of assuming the material parameters. Based on this method, the existing literature [32,33,34] has given theoretical solutions of thick-walled tubes under a single load or two loads acting together. In this topic, thick-walled tubes simultaneously under magnetic, thermal and mechanical loads are studied, and the results are compared with previous literature [32,33,34] when some parameters are set as special values.

2. Material Models and Properties

The infinitely long FGM thick-walled tube under uniform magnetic fields Hz is subjected to mechanical loads P and thermal loads T, as shown in Figure 1. Cylindrical polar coordinates (r, θ, z) are used, and the inner and outer radii are designated as a and b, respectively. Pa and Ta are the loads acting on the inner surface, and Pb and Tb are the loads acting on the outer surface.
The FGM thick-walled tube is composed of two distinct materials that are isotropic linear elastic Material A and linear elastic Material B. The interaction between the two materials is not considered. The volume fraction c(r) of Material A is described according to the power law [34].
c ( r ) = c 0 1 k ( r / b ) n
where, c 0 , k and n are the material parameters, r is the radius. Different distribution rules can be obtained by adjusting the material parameters, and the value of c(r) stays between 0 and 1.

3. Magneto-Thermo-Elastic Theoretical Solution

The relevant literature has presented the thermo-elastic theoretical solution [32] and magneto-elastic theoretical solution [33] of the thick-wall tube. In this paper, the magneto-thermo-elastic theoretical solution is further deduced on the basis of the above researches using the Voigt method, which considers that the strain of each phase material is equal and the stress is not equal when different phase materials come into contact with each other. The Voigt method can avoid making assumptions about parameters such as elastic modulus, Poisson’s ratio, thermal conductivity, thermal expansion coefficient, and magnetic permeability, and can obtain equivalent material parameters of FGMs that are more in accordance with the actual situation.
The stresses of the FGM thick-walled tube under thermal and mechanical loads are given as [32]
σ r = λ ¯ u r + λ ¯ + 2 G ¯ d u d r c ( r ) α 0 ( 3 λ 0 + 2 μ 0 ) + 1 c ( r ) α 1 ( 3 λ 1 + 2 μ 1 ) T ( r ) σ θ = λ ¯ + 2 G ¯ u r + λ ¯ d u d r c ( r ) α 0 ( 3 λ 0 + 2 μ 0 ) + 1 c ( r ) α 1 ( 3 λ 1 + 2 μ 1 ) T ( r ) σ z = λ ¯ u r + d u d r c ( r ) α 0 ( 3 λ 0 + 2 μ 0 ) + 1 c ( r ) α 1 ( 3 λ 1 + 2 μ 1 ) T ( r )
where, the superscript i = 0, 1, respectively, correspond to Materials A and B; λ i and G i are Lamé constants; the subscript r, θ, z represent radial, circumferential and axial direction, respectively; α i is the thermal expansion coefficient, u is the displacement in radial direction; λ ¯ and G ¯ are defined as
λ ¯ = c ( r ) λ 0 + 1 c ( r ) λ 1 G ¯ = c ( r ) G 0 + 1 c ( r ) G 1
The following assumptions are made: (a) each material component of the FGM tube is non-ferromagnetic and non-ferroelectric; (b) the Thompson effects are ignored; (c) the displacement currents are ignored; the radial Lorentz’s stress of the thick-walled tube subjected to magnetic field Hz is obtained based on consideration of Maxwell’s equations [33], described as
f r = H z 2 μ ( r ) 2 u r 2 + 1 r u r u r 2
where
μ ( r ) = c ( r ) μ 1 + 1 c ( r ) μ 2
Substituting Equations (2) and (4) into the equilibrium equation of tube subjected to a magnetic field in the cylindrical coordinate system d σ r d r + σ r σ θ r + f r = 0 [33], the ordinary differential equation for displacement can be induced as
r ϕ 1 ϕ 2 r n d 2 u d r 2 + ϕ 1 ϕ 3 r n d u d r ϕ 1 + ϕ 4 r n u r = f ( r )
where
ϕ 1 = c 0 λ 1 + 2 G 1 + μ 1 H z 2 + ( 1 c 0 ) λ 2 + 2 G 2 + μ 2 H z 2 ϕ 2 = c 0 k λ 1 + 2 G 1 + μ 1 H z 2 λ 2 2 G 2 μ 2 H z 2 / b n ϕ 3 = ( n + 1 ) ϕ 2 c 0 k n H z 2 ( μ 1 μ 2 ) / b n ϕ 4 = c 0 k n ( λ 1 λ 2 ) / b n ϕ 2 f ( r ) = A n r n ln r + B + C r 2 n + D r n
and constants in f ( r ) are
A = c 0 C 1 k b n α 1 3 λ 1 + 2 μ 1 α 0 3 λ 0 + 2 μ 0 B = C 1 α 1 1 c 0 3 λ 1 + 2 μ 1 + c 0 α 0 3 λ 0 + 2 μ 0 C = 2 k 0 k 1 c 0 k ( c 0 k 1 + k 0 c 0 k 0 ) b n A D = ( n C 2 + C 1 ) A C 1 + C B 2 A
As is well known, ϕ 1 includes parameters for Material A and Material B, which are generally not equal to zero, so Equation (6) can be transformed into
r 1 ϕ 2 ϕ 1 r n d 2 u d r 2 + 1 ϕ 3 ϕ 1 r n d u d r 1 + ϕ 4 ϕ 1 r n u r = C ϕ 1 r 2 n + n A ϕ 1 r n ln r + D ϕ 1 r n + B ϕ 1
Substituting x = χ ( r ) = ϕ 2 ϕ 1 r n into Equation (9)
x 2 1 x d 2 u d x 2 + x 1 n 1 + ϕ 3 / ϕ 2 n x d u d x 1 n 2 1 + ϕ 4 ϕ 2 x u = x 2 1 x g x
where
g ( x ) = ϕ 1 ϕ 2 x 1 / n C ϕ 1 ϕ 2 x 2 + A ϕ 2 x ln ϕ 1 ϕ 2 x + D ϕ 2 x + B ϕ 1 n 2 x 2 ( 1 x )
Equation (10) is a second-order inhomogeneous ordinary differential equations, and its solution includes two parts: the general solution u c ( x ) of the homogeneous equation and the specific solution u p ( x ) of the inhomogeneous equation. The general solution in the interval −1 < x < 1 can be expressed as
u c ( x ) = B 1 u 1 ( x ) + B 2 u 2 ( x ) = B 1 r F ( α , β , δ ; x ) + B 2 1 r F ( α δ + 1 , β δ + 1 , 2 δ ; x )
where, B1 and B2 are integration constant, F is the hypergeometric series, which has been defined in reference [32]. For common functionally graded material parameters, x is mainly concentrated in −1 < x < 1. For all the other x values, it can be converted to −1 < x < 1 by simple variable substitution, which will not be repeated here.
The coefficients in the hypergeometric series are
δ = 1 + 2 n ,   α = ϕ 3 / ϕ 2 1 2 4 ϕ 4 / ϕ 2 + ϕ 3 / ϕ 2 + 1 2 n ,   β = ϕ 3 / ϕ 2 + 1 n α
A specific solution of Equation (10) is easily obtained by the constant variation method.
u p ( x ) = u 1 ( x ) x a x u 2 ( t ) g ( t ) W ( t ) d t + u 2 ( x ) x a x u 1 ( t ) g ( t ) W ( t ) d t
where, x a = χ ( a ) , and Wronskian is
W ( x ) = u 1 ( x ) u 2 ( x ) u 2 ( x ) u 1 ( x )
Then, the solution of Equation (10) can be written as
u ( x ) = u 1 ( x ) B 1 x a x u 2 ( t ) g ( t ) W ( t ) d t + u 2 ( x ) B 2 + x a x u 1 ( t ) g ( t ) W ( t ) d t
Substituting Equation (16) into the Equation (2), stresses in all directions and the magnetic field of the FGM tube can be induced as
σ r = λ ¯ u ( x ) r + λ ¯ + 2 G ¯ u ( x ) c ( r ) α 0 ( 3 λ 0 + 2 μ 0 ) + 1 c ( r ) α 1 ( 3 λ 1 + 2 μ 1 ) T ( r ) σ θ = λ ¯ + 2 G ¯ u ( x ) r + λ ¯ u ( x ) c ( r ) α 0 ( 3 λ 0 + 2 μ 0 ) + 1 c ( r ) α 1 ( 3 λ 1 + 2 μ 1 ) T ( r ) σ z = λ ¯ u ( x ) r + u ( x ) c ( r ) α 0 ( 3 λ 0 + 2 μ 0 ) + 1 c ( r ) α 1 ( 3 λ 1 + 2 μ 1 ) T ( r ) h z = H z u ( x ) r + u ( x )
With stress boundary conditions σ r r = a = p a and σ r r = b = p b , the specific values of the constants B1 and B2 can be written as
B 1 = p ( λ ¯ ( b ) + 2 G ¯ ( b ) ) Q ( b ) + λ ¯ ( b ) Q ( b ) / b / B 0 B 2 = p ( λ ¯ ( b ) + 2 G ¯ ( b ) ) P ( b ) + λ ¯ ( b ) P ( b ) / b / B 0
where
B 0 = ( λ ¯ ( b ) + 2 G ¯ ( b ) ) Q ( b ) + λ ¯ ( b ) Q ( b ) / b ( λ ¯ ( a ) + 2 G ¯ ( a ) ) P ( a ) + λ ¯ ( a ) P ( a ) / a ( λ ¯ ( a ) + 2 G ¯ ( a ) ) Q ( a ) + λ ¯ ( a ) Q ( a ) / a ( λ ¯ ( b ) + 2 G ¯ ( b ) ) P ( b ) + λ ¯ ( b ) P ( b ) / b

4. Results and Discussion

The dimensionless representation of the radial coordinate, the inner radius, the stresses, the radial displacements and the magnetic field are r ¯ = r / b , a ¯ = a / b , σ ¯ i j = σ i j E 1 α 1 f 2 , u ¯ = u E 1 p a α 1 f 2 b , h ¯ z = h z E 1 H z p a α 1 f 2 . The parameters values are taken as a ¯ = 0.7 , c 0 = 1 , k = 1, E 1 = 70   Gpa , E 0 = 3 E 1 , μ 1 = μ 2 = 4 π × 10 7 H / m .
In the following simulations, the mechanical responses of the thick-walled tube within magnetic loads (Hz = 2.23 × 109 A/m), thermal loads (Ta = 0 °C, Tb = 100 °C) and mechanical loads (Pa = 1 GPa, Pb = 0 GPa) are discussed.

4.1. Effects of Parameter n

The results of different parameters n (n = 1.5, 3, 5, 10) are presented to discuss the influence of the volume fraction. In this section, v 0 = v 1 = 0.3 , α 1 / α 0 = 2 .
According to Equation (1), the volume fraction of Material A at key locations such as internal, middle and external positions are given, as shown in Figure 2. With the increase of parameter n, the volume fraction of Material A displays a nonlinear increasing trend. For n = 10, the volume fraction approaches 1.0 at the inner radius r ¯ = 0.7 . For different parameter n, the volume fraction of Material A decreases from the inner radius to outer radius, and equals zero at the outer radius r ¯ = 1.0 .
To make the results distinct and understandable, Table 2 lists the extremum values at the inner and outer radii corresponding to different parameter n.
Figure 3, Figure 4, Figure 5, Figure 6 and Figure 7 show the effect of parameter n on the radial displacement u ¯ , the radial stress σ ¯ r , the circumferential stress σ ¯ θ , the axial stress σ ¯ z and the perturbation of magnetic field h ¯ z , respectively.
Figure 3 shows that parameter n has a marked impact on the value of radial displacement u ¯ , but almost has no impact on the curve law. It is clear that the radial displacement u ¯ decreases by increasing the values of parameter n. The variable value of u ¯ from n = 1.5 to n = 3 is greater than that from n = 5 to n = 10, which indicates that the influence is relatively larger with smaller n, as listed in Table 2.
Comparing Figure 4, Figure 5 and Figure 6, it is found that the parameter n has significant influence on the circumferential stress σ ¯ θ and the axial stress σ ¯ z , but the influence on the radial stress σ ¯ r can be negligible. From the inner surface to the outer surface, the radial stress σ ¯ r shows the same nonlinear increasing curve shape with different parameter n, and all the curves obey the stress boundary condition that is σ r ¯ ( r ¯ ) = 1 at the inner radius ( r ¯ = 0.7 ) and σ r ¯ ( r ¯ ) = 0 at the outer radius ( r ¯ = 1.0 ), as shown in Figure 4 and listed in Table 2. The stresses in other directions ( σ ¯ θ or σ ¯ z ) show a downward trend from the inner radius to the outer radius, and the values exhibit fluctuations corresponding to different parameter n at the same radius. The curvature of curves are different and the curves intersect near the outer radius, as shown in Figure 5 and Figure 6.
Figure 7 reveals that the perturbation of magnetic field h ¯ z increases by increasing the parameter n, and all the values are negative with horizontal distribution law. The value of h ¯ z increases by about 0.11 from 1.5 to 3, and about 0.09 from 3 to 10, which indicates that the increase rate of h ¯ z shows a downward trend as the increase of parameter n.

4.2. Effects of Thermal Expansion Coefficient

In this section, four groups of thermal expansion coefficient (α10 = 0.5, 1, 2, 5) are analyzed to discuss the effects of thermal expansion coefficient. α10 represents the thermal expansion coefficient ratio of Materials A and B. In this section, n = 1.5, v 0 = v 1 = 0.3 .
To make the results distinct and understandable, Table 3 lists the extremum values at the inner and outer radii corresponding to different α 1 / α 0 .
Figure 8, Figure 9, Figure 10, Figure 11 and Figure 12 show the effect of the thermal expansion coefficient on the radial displacement u ¯ , the radial stress σ ¯ r , the circumferential stress σ ¯ θ , the axial stress σ ¯ z and the perturbation of magnetic field h ¯ z , respectively.
As depicted in Figure 8, the radial displacement u ¯ shows a clear increasing trend by increasing α 1 / α 0 , and also α 1 / α 0 hardly affect the decline law of u ¯ from inner to outer radius. Different with parameter n, the influence on u ¯ is relatively larger with larger α 1 / α 0 , which is reflected by the data in Table 3.
As shown in Figure 9, similar to the effect of parameter n on the radial stress σ ¯ r , all the curves obey the stress boundary condition. From the inner radius ( r ¯ = 0.7 ) to outer radius ( r ¯ = 1.0 ), the biggest impact occurs at r ¯ = 0.85 , where the radial stress σ ¯ r increases along with the increase of α 1 / α 0 . On the whole, α 1 / α 0 has minor influence on the radial stress σ ¯ r . Unlike the effect on radial stress σ ¯ r , Figure 10 and Figure 11 exhibit clearly that the effects of α 1 / α 0 on the circumferential stress σ ¯ θ and axial stress σ ¯ z are great. In general, the stresses σ ¯ θ (or σ ¯ z ) show an increase trend with the increase of α 1 / α 0 . The influence of α 1 / α 0 on the σ ¯ θ (or σ ¯ z ) is greater for positions close to the inner radius ( r ¯ = 0.7 ) than for those close to the outer radius ( r ¯ = 1.0 ). For the stress σ ¯ θ , the variation is about 0.56 with different α 1 / α 0 at the inner radius ( r ¯ = 0.7 ), and the variation is only 0.15 with different α 1 / α 0 at the outer radius ( r ¯ = 1.0 ), as listed in Table 3. For the stress σ ¯ z , the values with different α 1 / α 0 at the outer radius ( r ¯ = 1.0 ) are basically equal to 0.32. However, at the inner radius ( r ¯ = 0.7 ), the values with different α 1 / α 0 vary greatly and increase from 0.6011 to 0.7694.
The perturbation of magnetic field h ¯ z shown in Figure 12 decreases with the increase of α 1 / α 0 , and all the values are negative with the horizontal distribution law.

4.3. Effects of Poisson’s Ratio

In this subsection, the effects of Poisson’s ratio are dicussed. Poisson’s ratio of Material A ( v 0 ) is set as 0.3, and Poisson’s ratio of Material B ( v 1 ) changes are 0.2, 0.3, 0.4. In this section, n = 1.5, α 1 / α 0 = 2 .
To make the results distinct and understandable, Table 4 lists the extremum values at the inner and outer radii corresponding to different v 1 when v 0 = 0.3 .
Figure 13, Figure 14, Figure 15, Figure 16 and Figure 17 show the effect of Poisson’s ratio on the radial displacement u ¯ , the radial stress σ ¯ r , the circumferential stress σ ¯ θ , the axial stress σ ¯ z and the perturbation of magnetic field h ¯ z , respectively.
The variations of radial displacement u ¯ due to Poisson’s ratio are shown in Figure 13. It is clearly found that the radial displacement u ¯ shows a significant decreasing trend along with the increase of v 1 when v 0 is set as a constant value.
Differently to the influence of the two parameters mentioned above, Poisson’s ratio displays an obvious influence on all the stresses, as shown in Figure 14, Figure 15 and Figure 16. From the curves in Figure 14, the values of radial stress σ ¯ r at the inner radius and outer radius for different Poisson’s ratios are equal to −1.0 and 0.0, respectively, which indicates that all the curves obey the stress boundary condition, as listed in Table 4. However, the values from the inner radius to outer radius show significant fluctuations with the change of Poisson’s ratios, and the largest difference occurs around the middle of tube r ¯ = 0.85 . It is obviously seen from Figure 15 and Figure 16 that Poisson’s ratio has greater effect on circumferential stress σ ¯ θ and axial stress σ ¯ z near the inner radius than that close to the outer radius. Especially for the axial stress σ ¯ z , the same values can be obtained at the external position ( r ¯ = 1.0 ) corresponding to different Poisson’s ratios, but the values change from about 0.43 to about 1.10 at the internal position ( r ¯ = 0.7 ), as shown in Figure 16. Indicated in Figure 17, the increasing Poisson’s ratio v 1 results in significant change of magnetic field h ¯ z at the same radius r ¯ , and also Poisson’s ratio almost has no impact on the horizontal distribution law.
The above discussions indicate that the influence of Poisson’s ratio cannot be ignored when the analysis of functionally graded thick-walled tubes within mechanical, thermal and magnetic loads is conducted.

4.4. Some Special Cases

The characteristic of this paper is that the coupling effect of mechanical load, temperature load, and magnetic field are fully considered, and it has a wide range of applications; that is to say, the present work can degenerate to the same results as previous papers when some of the material parameters take special values. Based on the n = 1.5, α 1 / α 0 = 2 , v 0 = v 1 = 0.3 , the applied load values of three special cases are listed in Table 5.
For case 1, according to the theoretical formulas obtained in this paper and the theoretical formulas in reference [32], the values of material parameters at key locations such as internal ( r ¯ = 0.7 ), middle ( r ¯ = 0.85 ) and external ( r ¯ = 1.0 ) positions are obtained. Similarly, the same operations are performed for case 2 and case 3. All the results are listed in Table 6.
By calculating, it can be found that for different cases, the theoretical solutions obtained in this paper are consistent with the values obtained by the references, as listed in Table 6. This indicates (1) the accuracy of the theoretical research in this paper; (2) the wide application range of this research, which covers the above three working conditions well.

5. Conclusions

The magneto-thermo-elastic theoretical solutions for an infinitely long FGM thick-walled tube composed of two materials are investigated, and the influences of the parameter n, the thermal expansion coefficient and Poisson’s ratio are discussed. The results indicate that all the above parameters have explicit influences on the value of radial displacement, the circumferential stress, the axial stress and the perturbation of the magnetic field. For the radial stress, all the above parameters have no effect at the inner and outer boundaries. In the middle range, the influence of Poisson’s ratio on the radial stress is significant, and the influences of the other two parameters are not significant. The curve changing laws of the radial displacement, radial stress and perturbation of magnetic field are basically not affected by all the parameters. Furthermore, by comparing the research results in this paper with the previous works when some of the parameters take special values, a more extensive scope of application is illustrated. Research in this paper can provide effective guidance for engineering design to make FGM tubes with high reliability in structural performance when subjected to a multi-field environment.

Author Contributions

Conceptualization, T.L. and Y.L.; data curation, T.L. and J.L.; funding acquisition, T.L.; investigation, T.L. and J.L.; methodology, T.L., J.L. and X.L.; project administration, T.L.; supervision, T.L.; validation, J.L. and X.L.; writing—original draft, J.L. and X.L.; writing—review and editing, T.L., J.L. and X.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (No. 51808374) and the Foundation of Zhejiang Provincial Key Laboratory of Space Structures (No. 202108).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Data are contained within the article.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Khan, T.; Ning, Z.; Abdullah, A. State of the art review of functionally graded materials. In Proceedings of the 2nd International Conference on Computing, Mathematics and Engineering Technologies, Sukkur, Pakistan, 30–31 January 2019. [Google Scholar]
  2. Saleh, B.; Jiang, J.; Fathi, R.; Tareq, A.; Xu, Q.; Wang, L.; Song, D.; Ma, A. 30 years of functionally graded materials: An overview of manufacturing methods, applications and future challenges. Compos. Part B Eng. 2020, 201, 108376. [Google Scholar] [CrossRef]
  3. Pasha, A.; Rajaprakash, B.M. Functionally graded materials (FGM) fabrication and its potential challenges and applications. Mater. Today Proc. 2022, 52, 413–418. [Google Scholar] [CrossRef]
  4. Khanna, K.; Gupta, V.K.; Grover, N. Influence of anisotropy on creep in functionally graded variable thickness rotating disc. J. Strain Anal. Eng. Des. 2022, 57, 95–103. [Google Scholar] [CrossRef]
  5. Demirbas, M.D.; Ekici, R.; Apalak, M.K. Thermoelastic analysis of temperature-dependent functionally graded rectangular plates using finite element and finite difference methods. Mech. Adv. Mater. Struct. 2018, 27, 1–18. [Google Scholar] [CrossRef]
  6. Demirbas, M.D.; Apalak, M.K. Thermal Residual Stresses Analyses of Two-Dimensional Functionally Graded Circular Plates with Temperature-Dependent Material Properties. Int. J. Eng. Res. Dev. 2018, 10, 202–213. [Google Scholar]
  7. Neves, A.M.A.; Ferreira, A.J.M.; Carrera, E.; Cinefra, M.; Roque, C.M.C.; Jorge, R.M.N.; Soares, C.M.M. Static, free vibration and buckling analysis of isotropic and sandwich functionally graded plates using a quasi-3D higher-order shear deformation theory and a meshless technique. Compos. Part. B Eng. 2013, 44, 657–674. [Google Scholar] [CrossRef]
  8. Wang, Y.; Yang, C.; Zhang, Y.; Dong, S.; Li, L. Dynamics of a rotating hollow FGM beam in the temperature field. Rev. Adv. Mater. Sci. 2021, 60, 643–662. [Google Scholar] [CrossRef]
  9. Dorduncu, M.; Apalak, M.K.; Reddy, J.N. Stress wave propagation in adhesively bonded functionally graded cylinders: An improved model. J. Adhes. Sci. Technol. 2019, 33, 156–186. [Google Scholar] [CrossRef]
  10. Abd-El-Salam, M.R.; Abd-Alla, A.M.; Hosham, H.A. A numerical solution of magneto-thermoelastic problem in non-homogeneous isotropic cylinder by the finite-difference method. Appl. Math. Model. 2007, 31, 1662–1670. [Google Scholar] [CrossRef]
  11. Chi, L. Research on Vector Hysteresis Model of Electrical Steel Sheets under Coupling Effect of Magneto-Force-Thermal; Shenyang University of Technology: Shenyang, China, 2020. (In Chinese) [Google Scholar]
  12. Mukhopadhyay, S.B.; Roychoudhurj, S.K. Magneto-thermo-elastic interactions in an infinite isotropic elastic cylinder subjected to a periodic loading. Int. J. Eng. Sci. 1997, 35, 437–444. [Google Scholar] [CrossRef]
  13. Akbarzadeh, A.; Chen, Z. Thermo-Magneto-Electro-Elastic responses of rotating hollow cylinders. Mech. Adv. Mater. Struct. 2012, 21, 67–80. [Google Scholar] [CrossRef]
  14. Wu, D.; Yu, L.; Zhao, B.; Yang, L.; Gao, Y. The exact theory of axisymmetric electro-magneto-thermo-elastic circular cylinder. J. Intell. Mater. Syst. Struct. 2015, 26, 770–780. [Google Scholar] [CrossRef]
  15. Saidi, A.; Abouelregal, A. Thermoelastic Model with Higher-order Time-derivatives and Two Phase-lags for an Infinitely Long Cylinder under Initial Stress and Heat Source. J. Appl. Comput. Mech. 2021, 7, 277–291. [Google Scholar]
  16. Dai, H.; Fu, Y. Magnetothermoelastic interactions in hollow structures of functionally graded material subjected to mechanical loads. Int. J. Press. Vessel. Pip. 2007, 84, 132–138. [Google Scholar] [CrossRef]
  17. Arefi, M.; Rahimi, G.H.; Khoshgoftar, M.J. Exact solution of a thick walled functionally gradedpiezoelectric cylinder under mechanical, thermal andelectrical loads in the magnetic field. Smart Struct. Syst. 2012, 9, 427–439. [Google Scholar] [CrossRef]
  18. Heydari, E.; Hashemi, S.J.; Roknizadeh, S. Stress Analysis of Functionally Graded Cylindrical Shells under Magneto-Thermo-Elastic Loads. Am. J. Adv. Sci. Res. 2013, 2, 132–139. [Google Scholar]
  19. Loghman, A.; Parsa, H. Exact solution for magneto-thermo-elastic behaviour of double-walled cylinder made of an inner FGM and an outer homogeneous layer. Int. J. Mech. Sci. 2014, 88, 93–99. [Google Scholar] [CrossRef]
  20. Heydari, E.; Hashemi, S.J.; Roknizadeh, S.A. Magneto-Thermo-Electro-Elastic analysis of functionally graded intelligent cylindrical shells. Am. J. Adv. Sci. Res. 2014, 2, 140–147. [Google Scholar]
  21. Bayat, M.; Rahimi, M.; Saleem, M.; Mohazzab, A.H.; Wudtke, I.; Talebi, H. One-dimensional analysis for magneto-thermo-mechanical response in a functionally graded annular variable-thickness rotating disk. Appl. Math. Model. 2014, 38, 4625–4639. [Google Scholar] [CrossRef]
  22. Atrian, A.; Fesharaki, J.J.; Nourbakhsh, S.H. Thermo-electromechanical behavior of functionally graded piezoelectric hollow cylinder under non-axisymmetric loads. Appl. Math. Mech. 2015, 35, 939–954. [Google Scholar] [CrossRef]
  23. Dini, A.; Abolbashari, M.H. Hygro-thermo-electro-elastic response of a fun-ctionally graded piezoelectric cylinder resting on an elastic foundation subjected to non-axisymmetric loads. Int. J. Press. Vessel. Pip. 2016, 147, 21–40. [Google Scholar] [CrossRef]
  24. Mehditabar, A.; Rahimi, G.H.; Sadrabadi, S.A. Three-dimensional magne to-thermo-elastic analysis of functionally graded cylindrical shell. Appl. Math. Mech. 2017, 38, 479–494. [Google Scholar] [CrossRef]
  25. Bidgoli, M.O.; Loghman, A.; Arefi, M. Three-dimensional thermo-elastic analysis of a rotating cylindrical functionally graded shell subjected to mechanical and thermal loads based on the FSDT formulation. J. Appl. Mech. Tech. Phys. 2019, 60, 899–907. [Google Scholar] [CrossRef]
  26. Fesharaki, J.J.; Roghani, M. Thermo-mechanical behavior of a functionally graded hollow cylinder with an elliptic hole. Braz. Soc. Mech. Sci. 2020, 42, 66–81. [Google Scholar] [CrossRef]
  27. Saadatfar, M.; Aghaie-Khafri, M. Thermoelastic analysis of a rotating functionally graded cylindrical shell with functionally graded sensor and actuator layers on an elastic foundation placed in a constant magnetic field. J. Intell. Mater. Syst. Struct. 2016, 27, 512–527. [Google Scholar] [CrossRef]
  28. Dini, A.; Nematollahi, M.A.; Hosseini, M. Analytical solution for magne-to-thermo-elastic responses of an annular functionally graded sandwich disk by considering internal heat generation and convective boundary condition. J. Sandw. Struct. Mater. 2021, 23, 542–567. [Google Scholar] [CrossRef]
  29. Alavi, F.; Karimi, D.; Bagri, A. An investigation on thermoelastic behaviour of functionally graded thick spherical vessels under combined thermal and mechanical loads. J. Achiev. Mater. Manuf. Eng. 2008, 31, 422–428. [Google Scholar]
  30. Nie, G.; Zhong, Z.; Batra, R.C. Material tailoring for functionally graded hollow cylinders and spheres. Compos. Sci. Technol. 2011, 71, 666–673. [Google Scholar] [CrossRef]
  31. Li, X.; Peng, X. A pressurized functionally graded hollow cylinder with arbitrarily varying material properties. J. Elast. 2009, 96, 81–95. [Google Scholar] [CrossRef]
  32. Xin, L.; Dui, G.; Yang, S.; Zhou, D. Solutions for behavior of a functionally graded thick-walled tube subjected to mechanical and thermal loads. Int. J. Mech. Sci. 2015, 98, 70–79. [Google Scholar] [CrossRef]
  33. Xin, L.; Li, Y.; Pan, D.; Dui, G.; Ju, C. Revisiting the elastic solution for an inner-pressured functionally graded thick-walled tube within a uniform magnetic field. Appl. Math. Mech. 2018, 39, 1485–1498. [Google Scholar] [CrossRef]
  34. Xin, L.; Dui, G.; Yang, S.; Zhang, J. An elasticity solution for functionally graded thick-walled tube subjected to internal pressure. Int. J. Mech. Sci. 2014, 89, 344–349. [Google Scholar] [CrossRef]
Figure 1. A diagram of a long FGM tube subjected to mechanical loads, thermal loads and magnetic field.
Figure 1. A diagram of a long FGM tube subjected to mechanical loads, thermal loads and magnetic field.
Materials 15 06345 g001
Figure 2. Evolution of volume fraction of Material A with different parameter n (c0 = 1, k = 1).
Figure 2. Evolution of volume fraction of Material A with different parameter n (c0 = 1, k = 1).
Materials 15 06345 g002
Figure 3. Effects of n on the radial displacement u ¯ .
Figure 3. Effects of n on the radial displacement u ¯ .
Materials 15 06345 g003
Figure 4. Effects of n on the radial stress σ ¯ r .
Figure 4. Effects of n on the radial stress σ ¯ r .
Materials 15 06345 g004
Figure 5. Effects of n on the circumferential stress σ ¯ θ .
Figure 5. Effects of n on the circumferential stress σ ¯ θ .
Materials 15 06345 g005
Figure 6. Effects of n on the axial stress σ ¯ z .
Figure 6. Effects of n on the axial stress σ ¯ z .
Materials 15 06345 g006
Figure 7. Effects of n on the perturbation of magnetic field h ¯ z .
Figure 7. Effects of n on the perturbation of magnetic field h ¯ z .
Materials 15 06345 g007
Figure 8. Effects of α 1 / α 0 on the radial displacement u ¯ .
Figure 8. Effects of α 1 / α 0 on the radial displacement u ¯ .
Materials 15 06345 g008
Figure 9. Effects of α 1 / α 0 on the radial stress σ ¯ r .
Figure 9. Effects of α 1 / α 0 on the radial stress σ ¯ r .
Materials 15 06345 g009
Figure 10. Effects of α 1 / α 0 on the circumferential stress σ ¯ θ .
Figure 10. Effects of α 1 / α 0 on the circumferential stress σ ¯ θ .
Materials 15 06345 g010
Figure 11. Effects of α 1 / α 0 on the axial stress σ ¯ z .
Figure 11. Effects of α 1 / α 0 on the axial stress σ ¯ z .
Materials 15 06345 g011
Figure 12. Effects of α 1 / α 0 on the perturbation of magnetic field h ¯ z .
Figure 12. Effects of α 1 / α 0 on the perturbation of magnetic field h ¯ z .
Materials 15 06345 g012
Figure 13. Effects of Poisson’s ratio on the radial displacement u ¯ .
Figure 13. Effects of Poisson’s ratio on the radial displacement u ¯ .
Materials 15 06345 g013
Figure 14. Effects of Poisson’s ratio on the radial stress σ ¯ r .
Figure 14. Effects of Poisson’s ratio on the radial stress σ ¯ r .
Materials 15 06345 g014
Figure 15. Effects of Poisson’s ratio on the circumferential stress σ ¯ θ .
Figure 15. Effects of Poisson’s ratio on the circumferential stress σ ¯ θ .
Materials 15 06345 g015
Figure 16. Effects of Poisson’s ratio on the axial stress σ ¯ z .
Figure 16. Effects of Poisson’s ratio on the axial stress σ ¯ z .
Materials 15 06345 g016
Figure 17. Effects of Poisson’s ratio on the perturbation of magnetic field h ¯ z .
Figure 17. Effects of Poisson’s ratio on the perturbation of magnetic field h ¯ z .
Materials 15 06345 g017
Table 1. The assumptions corresponding to different problems.
Table 1. The assumptions corresponding to different problems.
Research ContentsAssumptions
response under mechanical loadselastic modulus, Poisson’s ratio
response under thermal loadsthermal expansion coefficient, thermal conductivity
response within magnetic fieldsmagnetic permeability
Table 2. The extremum values corresponding to different parameter n.
Table 2. The extremum values corresponding to different parameter n.
n1.53510
u ¯ r ¯ = 0.7 1.26081.01630.89460.8150
r ¯ = 1.0 1.04100.84070.74090.6756
σ ¯ r r ¯ = 0.7 −1−1−1−1
r ¯ = 1.0 0000
σ ¯ θ r ¯ = 0.7 3.19083.26353.31253.3376
r ¯ = 1.0 1.12390.90390.79420.7225
σ ¯ z r ¯ = 0.7 0.65720.67900.69370.7012
r ¯ = 1.0 0.32310.25710.22420.2027
h ¯ z r ¯ = 0.7 −0.6229−0.5086−0.4514−0.4129
r ¯ = 1.0 −0.6287−0.5064−0.4494−0.4121
Table 3. The extremum values corresponding to different α 1 / α 0 .
Table 3. The extremum values corresponding to different α 1 / α 0 .
α 1 / α 0 0.5125
u ¯ r ¯ = 0.7 1.19561.21731.26081.3910
r ¯ = 1.0 0.98181.00151.04101.1593
σ ¯ r r ¯ = 0.7 −1−1−1−1
r ¯ = 1.0 0000
σ ¯ θ r ¯ = 0.7 3.00393.06623.19083.5647
r ¯ = 1.0 1.07391.09061.12391.2239
σ ¯ z r ¯ = 0.7 0.60110.61980.65720.7694
r ¯ = 1.0 0.31870.32020.32310.3321
h ¯ z r ¯ = 0.7 −0.5698−0.5875−0.6229−0.7293
r ¯ = 1.0 −0.5675−0.5853−0.6208−0.7274
Table 4. The extremum values corresponding to different v 1 .
Table 4. The extremum values corresponding to different v 1 .
v 1 0.20.30.4
u ¯ r ¯ = 0.7 1.89831.26080.9198
r ¯ = 1.0 1.64901.04100.7133
σ ¯ r r ¯ = 0.7 −1−1−1
r ¯ = 1.0 000
σ ¯ θ r ¯ = 0.7 5.02063.19082.2340
r ¯ = 1.0 1.70021.12390.8259
σ ¯ z r ¯ = 0.7 1.10410.65720.4278
r ¯ = 1.0 0.32600.32310.3163
h ¯ z r ¯ = 0.7 −1.2549−0.6229−0.2753
r ¯ = 1.0 −1.2577−0.6208−0.2704
Table 5. The assumptions corresponding to different problems.
Table 5. The assumptions corresponding to different problems.
Case/ReferParametersLoad Conditions
1/[32]Pa = 1 GPa, Pb = 0 GPa, Ta = 0 °C, Tb = 100 °C, Hz = 0 A/mmechanical and thermal loads
2/[33]Pa = 1 GPa, Pb = 0 GPa, Ta = 0 °C, Tb = 0 °C, Hz = 2.23 × 109 A/mmechanical load and magnetic field
3/[34]Pa = 1 GPa, Pb = 0 GPa, Ta = 0 °C, Tb = 0 °C, Hz = 0 A/mmechanical load
Table 6. The values of material parameters of different cases.
Table 6. The values of material parameters of different cases.
LocationsCaseTheoretical Reference u ¯ σ ¯ r σ ¯ θ σ ¯ z h ¯ z
r ¯ = 0.7 Case 1This paper1.3861−1.00003.55050.76520
[32]
r ¯ = 0.85 This paper1.2364−0.33212.13480.53170
[32]
r ¯ = 1.0 This paper1.1423−0.03501.22020.34160
[32]
r ¯ = 0.7 Case 2This paper1.1740−1.00002.94170.5825−0.5521
[33]
r ¯ = 0.85 This paper1.0422−0.28981.80610.4549−0.5504
[33]
r ¯ = 1.0 This paper0.96220.00001.05730.3172−0.5498
[33]
r ¯ = 0.7 Case 3This paper1.4399−1.00003.70510.81150
[34]
r ¯ = 0.85 This paper1.2860−0.30852.24980.58240
[34]
r ¯ = 1.0 This paper1.18860.00001.30610.39180
[34]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, T.; Li, J.; Liu, X.; Luo, Y. Magneto-Thermo-Elastic Theoretical Solution for Functionally Graded Thick-Walled Tube under Magnetic, Thermal and Mechanical Loads Based on Voigt Method. Materials 2022, 15, 6345. https://doi.org/10.3390/ma15186345

AMA Style

Li T, Li J, Liu X, Luo Y. Magneto-Thermo-Elastic Theoretical Solution for Functionally Graded Thick-Walled Tube under Magnetic, Thermal and Mechanical Loads Based on Voigt Method. Materials. 2022; 15(18):6345. https://doi.org/10.3390/ma15186345

Chicago/Turabian Style

Li, Tiane, Jiabao Li, Xuekang Liu, and Yaozhi Luo. 2022. "Magneto-Thermo-Elastic Theoretical Solution for Functionally Graded Thick-Walled Tube under Magnetic, Thermal and Mechanical Loads Based on Voigt Method" Materials 15, no. 18: 6345. https://doi.org/10.3390/ma15186345

APA Style

Li, T., Li, J., Liu, X., & Luo, Y. (2022). Magneto-Thermo-Elastic Theoretical Solution for Functionally Graded Thick-Walled Tube under Magnetic, Thermal and Mechanical Loads Based on Voigt Method. Materials, 15(18), 6345. https://doi.org/10.3390/ma15186345

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