Transient Thermal Stress Problem of a Functionally Graded Magneto-Electro-Thermoelastic Hollow Sphere.

This article is concerned with the theoretical analysis of the functionally graded magneto-electro-thermoelastic hollow sphere due to uniform surface heating. We analyze the transient thermoelastic problem for a functionally graded hollow sphere constructed of the spherical isotropic and linear magneto-electro-thermoelastic materials using a laminated composite mode as one of theoretical approximation in the spherically symmetric state. As an illustration, we carry out numerical calculations for a functionally graded hollow sphere constructed of piezoelectric and magnetostrictive materials and examine the behaviors in the transient state. The effects of the nonhomogeneity of material on the stresses, electric potential, and magnetic potential are investigated.

1. Introduction

Functionally graded materials (FGMs) are new nonhomogeneous material systems that two or more different material ingredients changes continuously and gradually. The concept of FGMs is applicable to many industrial fields such as aerospace, nuclear energy, chemical plant, electronics and so on. On the other hand, it has recently been found that composites made of piezoelectric and magnetostrictive materials exhibit the magnetoelectric effect, which is not seen in piezoelectric or magnetostrictive materials [1]. These materials are known as multiferroic composites [2]. These composites exhibit a coupling among magnetic, electric, and elastic fields. In the past, various problems in magneto-electro-elastic media that exhibit anisotropic and linear coupling among the magnetic, electric, and elastic fields were analyzed. Examples for the plates and beams were analyzed in the papers [3,4,5]. Examples for the shell type structures were analyzed in the papers [6,7,8]. Examples of functionally graded magneto-electro-elastic media are as follows. Wang and Ding [9] treated spherically symmetric transient responses of a functionally graded magneto-electro-elastic hollow sphere. Ma and Lee [10] analyzed an in-plane problem in functionally graded magneto-electro-elastic bimaterials. Yu and Wu [11] analyzed the propagation of circumferential wave in magneto-electro-elastic functionally graded cylindrical curved plates. Wu and Lu [12] analyzed the 3D dynamics responses of functionally graded magneto-electro-elastic plates. Huang et al. [13] analyzed the static problem of an anisotropic functionally graded magneto-electro-elastic beams subjected to arbitrary loading. Lee and Ma [14] analyzed the two-dimensional problem of two bonded dissimilar half-planes for functionally graded magnetoelectroelastic materials subjected to generalized line forces and screw dislocations.

Examples of the thermal stress problems of electro-magneto-elastic media are as follows, Ganesan et al. [15] analyzed the response of a layered, multiphase magnetoelectroelastic cylinder subjected to an axisymmetric temperature distribution using finite element procedures. Kumaravel et al. [16] analyzed the response of a three-layered magnetoelectroelastic strip subjected to uniform temperature rise and non-uniform temperature distribution using finite element procedures. Hou et al. [17] obtained 2D fundamental solutions of a steady point heat source in infinite and semi-infinite orthotropic electro-magneto-thermo-elastic planes. With regard to transient thermal stress problems of electro-magneto-elastic media, Wang and Niraula [18] analyzed transient thermal fracture in transversely isotropic electro-magneto-elastic cylinders. The exact solution of a transient analysis of multilayered magneto-electro-thermoelastic strip subjected to nonuniform heat supply was obtained in the paper [19]. The exact solution of a transient analysis of multilayered magneto-electro-thermoelastic hollow cylinder subjected to uniform heat supply was obtained in the paper [20]. Though a several transient thermal stress problems of the functionally graded hollow spheres [21,22] using a laminated composite model were analyzed already, theses studies don’t consider a coupling among magnetic, electric, and thermoelastic fields. However, to the author’s knowledge, the transient thermal stress problem for a functionally graded magneto-electro-thermoelastic hollow spheres under unsteady heat supply considering a coupling among magnetic, electric, and thermoelastic fields has not been reported.

In the present article, we have analyzed the transient behavior of a functionally graded magneto-electro-thermoelastic hollow sphere due to uniform surface heating. We assumed that the magneto-electro-thermoelastic materials are polarized and magnetized in the radial direction. We analyze the transient thermal stress problem for a functionally graded hollow sphere constructed of the spherical isotropic and linear magneto-electro-thermoelastic materials using a laminated composite model as one of theoretical approximation. We carried out numerical calculations for a functionally graded hollow sphere composed of piezoelectric and magnetostrictive materials, and examined the effects of the nonhomogeneity of material on the stresses, electric potential, and magnetic potential.

2. Analysis

We consider a functionally graded hollow sphere constructed of the spherical isotropic and linear magneto-electro-thermoelastic materials. We analyze the transient thermal stress problem using a multilayered composite hollow sphere model with a number N of homogeneous layers. The hollow sphere’s inner and outer radii are designated a and b, respectively. r is the outer radius of the ith layer. Throughout this article, the indices i (=1,2,…,N) are associated with the ith layer of a composite hollow sphere from the inner side.

2.1. Heat Conduction Problem

We assumed that the multilayered hollow sphere is initially at zero temperature and its inner and outer surfaces are suddenly heated by surrounding media having constant temperatures T and T with relative heat transfer coefficients h and h, respectively. Then, the temperature distribution is one-dimensional, and the transient heat conduction equation for the ith layer is written in the following form:

The initial and thermal boundary conditions in dimensionless form are

In Equations (1)–(6), we introduced the following dimensionless values: where T is the temperature change; t is time; λ is the thermal conductivity in the radial direction; κ is the thermal diffusivity in the radial direction; and T, λ and κ are typical values of temperature, thermal conductivity, and thermal diffusivity, respectively. To solve the fundamental equation (1), we introduced the Laplace transformation with respect to the variable τ as follows;

Performing the Laplace transformation on Equation (1) under the condition of Equation (2) gives where

The general solution of Equation (9) is where j0() and y0() are zeroth-order spherical Bessel functions of the first and second kind, respectively. Furthermore, A and B are unknown constants determined from the boundary conditions. Substituting Equation (11) into the boundary conditions in the transformed domain from Equations (3)–(6), these equations are represented in matrix form as follows:

Making use of Cramer’s formula, the constants A and B can be determined from Equation (12). Then the temperature solution in the transformed domain is where Δ is the determinant of 2N × 2N matrix [a], and the coefficients and are defined as determinants of a matrix similar to the coefficient matrix [a], in which the (2i − 1)th column or 2ith column is replaced with the constant vector {c}, respectively. Using the residue theorem, we can accomplish the inverse Laplace transformation on Equation (13). Because the single-valued poles of Equation (13) correspond to p = 0 and the roots of Δ = 0, in which the residue for p = 0 gives a solution for the steady state. Accomplishing the inverse Laplace transformation of Equation (13), the solution of Equation (1) is written as follows: where F ie the determinants of 2N × 2N matrix [e], and the coefficients and are defined as determinants of a matrix similar to the coefficient matrix [e], in which the (2i − 1)th column or 2ith column is replaced with the constant vector {c}, respectively. The nonzero elements of the coefficient matrices [a] and [e] and the constant vector {c} are given from the Equations (3)–(6). In Equation (14), is and μ is the jth positive root of the following transcendental equation

2.2. Thermoelastic Problem

We developed the analysis of a multilayered magneto-electro-thermoelastic hollow sphere as a spherically symmetric state. The displacement-strain relations are expressed in dimensionless form as follows: where the comma denotes partial differentiation with respect to the variable that follows. For the spherical isotropic and linear magneto-electro-thermoelastic material, the constitutive relations are expressed in dimensionless form as follows: where

The constitutive equations for the electric and the magnetic fields in dimensionless form are given as

The relation between the electric field intensity and the electric potential in dimensionless form is defined as

The relation between the magnetic field intensity and the magnetic potential in dimensionless form is defined as

The equilibrium equation is expressed in dimensionless form as follows:

If the electric charge density is absent, the equations of electrostatics and magnetostatics are expressed in dimensionless form as follows:

In Equations (17)–(26), the following dimensionless values are introduced: where are the stress components; (, ) are the strain components; u is the displacement in the r direction; α are the coefficients of linear thermal expansion; C are the elastic stiffness constants; D is the electric displacement in the r direction; B is the magnetic flux density in the r direction; e are the piezoelectric coefficients; η1 is the dielectric constant; P1 is the pyroelectric constant; q are the piezomagnetic coefficients; μ1 is the magnetic permeability coefficient; d1 is the magnetoelectric coefficient; m1 is the pyromagnetic constant; and α0, Y0 and d0 are typical values of the coefficient of linear thermal expansion, Young’s modulus, and piezoelectric modulus, respectively.

Substituting Equations (17), (22), and (23) into Equations (18), (20), and (21) and later into Equations (24)–(26), the governing equations of the displacement u, electric potential φ, and magnetic potential ψ in the dimensionless form are written as

If the inner and outer surfaces of the multilayered magneto-electro-thermoelastic hollow sphere are traction free, and the interfaces of each adjoining layer are perfectly bonded, then the boundary conditions of inner and outer surfaces and the conditions of continuity at the interfaces can be represented as follows:

The boundary conditions in the radial direction for the electric and magnetic fields are expressed as

The solutions of Equations (28)–(30) are assumed in the following form:

In Equation (33), the first term on the right-hand side gives the homogeneous solution and the second term gives the particular solution. We now consider the homogeneous solution, and introduce the following equation:

Changing a variable with the use of Equation (34), the homogeneous expression of Equations (28)–(30) are where

By eliminating and between Eqs. (35)-(37), we can obtain an ordinary differential equation about :

The solution of Equation (40) can be expressed as follows when .

From Equations (36), (37) and (41), we can obtain an ordinary differential equation about :

Using Equation (39), the solution of Equation (42) is

From Equations (36), (41) and (43), we can obtain an ordinary differential equation about :

Using Equation (39), the solution of Equation (44) is

In Equations (41)–(45)

In Equations (41), (43) and (45), C () are unknown constants. We have the following relation.

The homogeneous solutions when are omitted here for brevity.

It is difficult to obtain the particular solutions using the temperature solution of Equation (14). In order to obtain the particular solutions, series expansions of Bessel functions given in Equation (14) are used. Equation (14) can be written in the following way: where Here, is the Kronecker delta. The particular solutions , , and are obtained as the function system like Equation (48). Then, the stress components, electric displacement, and magnetic flux density can be evaluated from Equations (41), (43) and (45). Details of the solutions are omitted from here for brevity. The unknown constants in the homogeneous solutions are determined so as to satisfy the boundary conditions in (31) and (32).

3. Numerical Results

To illustrate the foregoing analysis, we consider the functionally graded hollow sphere composed of piezoelectric and magnetostrictive materials. The piezoelectric material is made up of BaTiO3, and the magnetostrictive material is made up of CoFe2O4. Numerical parameters of heat conduction and shape are presented as follows:

The first layer is pure piezoelectric material and the Nth layer is pure magnetostrictive material. It is assumed that the volume fractions of the piezoelectric phase V and the magnetostrictive phase V for other layers are given by the relations

The value of V in ith layer is obtained by calculating the value of V in Equation (34) at the centre point of each layer defined by . To estimate the material properties of FGM, we apply the simplest linear law of mixture. The material constants considered for BaTiO3 and CoFe2O4 are shown in the paper [20]. The typical values of material parameters such as κ, λ, α, Y, and d, used to normalize the numerical data, based on those of BaTiO3 are as follows:

In the numerical calculations, the boundary conditions at the surfaces for the electric and magnetic fields are expressed as

Figure 1, Figure 2, Figure 3, Figure 4 and Figure 5 show the numerical results for M = 1 and N = 10. The variations of temperature change and displacement , along the radial direction are shown in Figure 1 and Figure 2, respectively. From Figure 1 and Figure 2, it is clear that the temperature and displacement increase with time and have the largest values in the steady state. The variations of normal stresses and along the radial direction are shown in Figure 3a and Figure 3b, respectively. Figure 3a reveals that the maximum tensile stress of occurs in the transient state and the maximum compressive stress of occurs in the steady state. From Figure 3b, it is clear that the maximum tensile stress occurs near the outer surface. The variations of electric potential and magnetic potential along the radial direction are shown in Figure 4 and Figure 5, respectively. Figure 4 reveals that the absolute value of the electric potential increases with time, and attains its maximum value in the steady state. The electric potential is almost zero in the tenth layer, i.e. the pure magnetostrictive layer. From Figure 5, it is clear that the absolute value of the magnetic potential increases with time and attains its maximum value in the steady state. The magnetic potential is almost constant in the first layer, i.e. the pure piezoelectric layer.

Figure 1

Variation of the temperature change (M = 1, N = 10).

Figure 2

Variation of the displacement (M = 1, N = 10)

Figure 3

Variation of the thermal stresses (M = 1, N = 10): (a) normal stress ; (b) normal stress .

Figure 4

Variation of the electric potential (M = 1, N = 10).

Figure 5

Variation of the magnetic potential (M = 1, N = 10).

In order to assess the effect of the nonhomogeneous parameter M on the stresses, electric potential, and magnetic potential, the numerical results for N = 10 are shown in Figure 6, Figure 7 and Figure 8. M = 4 shows a piezoelectric material rich, and M = 0.25 shows a magnetostrictive material rich. The variations of stresses and are shown in Figure 6a and Figure 6b, respectively. From Figure 6a, it is clear that the maximum compressive stress of decreases when the parameter M increases in the steady state. From Figure 6b, it is clear that the maximum tensile stress of decreases when the parameter M decreases in the steady state. The variations of electric potential and magnetic potential are shown in Figure 7 and Figure 8, respectively. From Figure 7 and Figure 8, the absolute value of the electric potential in the inner surface is maximum when the parameter M = 1 in the steady state, while that of the magnetic potential is maximum when the parameter M = 0.25 in the steady state.

Figure 6

Variation of the thermal stresses (N = 10): (a) normal stress ; (b) normal stress .

Figure 7

Variation of the electric potential (N = 10).

Figure 8

Variation of the magnetic potential (N = 10).

In order to assess the effect of relaxation of stress values in functionally graded magneto-electro-thermoelastic hollow sphere, the numerical results for the two-layered hollow sphere are shown in Figure 9. Figure 9a, Figure 9b, Figure 9c and Figure 9d show the variations of stresses , , electric potential and magnetic potential, respectively. From Figure 3 and Figure 9, the effect of relaxation of stress distributions for the functionally graded hollow sphere can be clearly seen compared with the two-layered hollow sphere. From Figure 4, Figure 5 and Figure 9, it is clear that the maximum absolute values of the electric potential and magnetic potential for functionally graded hollow sphere are grater than those for the two-layered hollow sphere.

Figure 9

Numerical results for the two-layered hollow sphere (N = 2): (a) normal stress ; (b) normal stress ; (c) electric potential; (d) magnetic potential.

4. Conclusions

In this study, we analyzed the transient thermal stress problem for the functionally graded magneto-electro-thermoelastic hollow sphere due to uniform surface heating using a laminated composite mode by solving the governing equations of the displacement, electric potential and magnetic potential. As an illustration, we carried out numerical calculations for a functionally graded hollow sphere composed of piezoelectric BaTiO3 and magnetostrictive CoFe2O4, and examined the behaviors in the transient state for temperature change, displacement, stress, electric potential and magnetic potential distributions. We investigated the effects of the nonhomogeneity of material on the stresses, electric potential, and magnetic potential. Furthermore, the effect of relaxation of stress values in functionally graded magneto-electro-thermoelastic hollow sphere was investigated. We conclude that we can evaluate not only the thermoelastic response of the functionally graded magneto-electro-thermoelastic hollow sphere, but also the electric and magnetic fields of functionally graded magneto-electro-thermoelastic hollow sphere quantitatively in a transient state.