Mathematical Problems in Engineering

Volume 2013 (2013), Article ID 947693, 13 pages

http://dx.doi.org/10.1155/2013/947693

Research Article

Three-Dimensional Piezothermoelastic Stress of a Finite Functionally Graded Cylindrical Shell with Piezoelectric Layer

Yong-gang Kang, Zhong-qi Wang, and Gongnan Xie

The Key Laboratory of Contemporary Design and Integrated Manufacturing Technology, Ministry of Education, Northwestern Polytechnical University, P.O. Box 409, Xi’An, Shaanxi 710072, China

Received 27 September 2013; Accepted 27 October 2013

Academic Editor: Manyu Xiao

Copyright © 2013 Yong-gang Kang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

Three-dimensional piezothermoelastic solutions for a finite functionally graded cylindrical shell with piezoelectric layer are carried out in this paper. The cylindrical shell is simply supported at four end edges and is subjected to axisymmetric thermomechanical loads. The piezoelectric layers are polarized along radial direction as a sensor. The material properties are assumed to be temperature independent and radially dependent but are assumed to be homogeneous in each layer; the variables are expanded in Fourier series to satisfy the boundary conditions and multilayer approach is used. Numerical results of mullite/molybdenum functionally graded cylindrical shell are presented; the temperature change, stresses, electric potential, and electric displacement distributions are given and briefly discussed.

1. Introduction

Functionally graded materials (FGMs) are new kind of inhomogeneous composite materials with continuously varying properties. The spatial grading in composition and microstructure of FGMs can be adjusted in order to obtain specific properties or to reduce the magnitude of mechanical and thermal stresses. Piezoelectric materials are probably the most popular active materials which act as sensors and actuators and have been used widely in structural health monitoring, vibration, and noise control as well as many other areas. Among these, piezoelectric bimorph and multimorph are commonly employed. Usually, these structures are made of two or more layers of piezoelectric sheets and are jointed by bonding agents. However, such laminated piezoelectric structures suffer from high stress concentration near the interface due to the abrupt changes in both material composition and thermoelectroelastic properties, which can cause severe deterioration of the bonding layer strength and reduce the lifetime of the structures. To overcome the drawbacks and meet some particular requirements for performance and reliability, the functionally graded materials (FGMs) and piezoelectric material are used together, which called functionally graded piezoelectric materials (FGPMs). A new type of structure was developed.

The fabrication and property investigation of FGPMs have attracted great attention from the research community. Among the early investigators, PintoCarreia et al. [1] derived and solved the coupled displacement and electrical field equations for a piezoelectric cylindrical shell based on the third-order shear deformation theory and the finite element method. Using classical laminate theory, the piezothermoelastic behavior of a piezoelectric composite actuator with functionally graded microstructure was analyzed [2]. Ootao and Tanigawa [3–7] exactly analyzed the transient piezothermoelastic problem of a functionally graded hollow sphere, a simply supported cylindrical composite panel which composed of crossply or angleply and piezoelectric layer under the state of generalized plane deformation, a rectangular composite plate composed of cross-ply laminate and piezoelectric material, and a functionally graded rectangular plate bonded to a piezoelectric plate due to partial heat supply. Using the power series expansion method, Chen and his coworkers [8, 9] successfully derived two-dimensional exact elasticity solution and studied both the direct and inverse piezoelectric effects of the piezothermoelastic behavior of circular laminated piezoelectric cylindrical shells and functionally graded piezoelectric cylindrical shells subjected to axisymmetric loading. Kapuria and Achary [10] presented a new coupled consistent third-order theory for static thermoelectromechanical response of hybrid piezoelectric laminated plates. Shao [11] presented solutions for temperature, displacement and stress fields in an FG hollow cylinder using a multilayered approach based on the laminated composite theory. Huang et al. [12] derived the analytical solutions for FGP beams under both mechanical and electrical loads from the two-dimensional equations of piezoelectricity, in which the elastic, piezoelectric and dielectric coefficients of the piezoelectric beams were assumed to vary along the beam thickness direction only. Yang and Xiang [13] used the Timoshenko beam theory to investigate the static bending and dynamic response of FGP actuators under combined thermalelectromechanical loading. Ying and Wang [14] employed the separation of variable technique to carry out two-dimensional elastic-dynamic analysis of simply supported hollow cylinders of finite length excited by nonuniform thermal shock. Reference [15] proposed improved layered shell finite element solution for the coupled thermoelectromechanical responses of smart fiber reinforced composite shell panels under piezothermoelastic loading and showed that piezoelectric effect has a significant influence on the response of such shells. The bending behavior of a circularly curved FGP cantilever actuator under an applied electrical load and heat conduction was investigated [16]. Infinitesimal axisymmetric deformation of a functionally graded shell with piezoelectric layers perfectly bonded to its inner and outer surfaces subjected to thermo-electro-mechanical loads was studied by Alibeigloo [17]. Alashti and Khorsand [18] carried out three-dimensional static analyses of FG cylindrical shells with piezoelectric layers under the effect of thermoelectromechanical loads, using the differential quadrature method.

In this paper, a functionally graded cylindrical shell with piezoelectric layer is investigated. The thermal and mechanical loads are applied on the cylinder asymmetrically. In order to obtain analytical solutions for temperature, displacements, and stresses for the three-dimensional piezothermoelastic problem, it is assumed that the piezoelectric functionally graded cylindrical shell is composed of fictitious layers and one piezoelectric layer, and each layer of the functionally graded material is homogeneous and isotropic.

2. Basic Equations

A piezoelectric functionally graded cylindrical shell with finite length , internal radius , external radius , and circumferential angle is considered (as shown in Figure 1), and it is assumed that the piezoelectric layer is perfectly bonded, the inner radius of piezoelectric layer is , and the radius of the interface between the th and th layer is . Cylindrical coordinates , , and are used in analysis. This functionally graded cylindrical shell is simply supported at its four end edges and subjected to nonuniform steady-state thermal loads and on the inner and outer surface. The temperatures at end edges keep zero. Nonuniform pressure loadings and are applied to the inner and outer surfaces of the panel, respectively. Details of piezoelectric functionally graded cylindrical shell are shown in Figure 1.

Figure 1: A piezoelectric functionally graded cylindrical shell.

It is assumed that the piezoelectric functionally graded cylindrical shell is composed of fictitious layers and one piezoelectric layer, and each layer of the functionally graded material is homogeneous and isotropic. The law of mixtures proposed by Kerner, which was derived based on the assumption that there is a granular phase embedded in a matrix phase, is applied to describe the variations of material properties in the radial direction of the functionally graded cylinder shell. Young’s modulus , Poisson’s ratio , thermal conductivity coefficient , and thermal expansion coefficient can be expressed as

where the subscripts and represent the metal and ceramic, respectively. is the volume fraction in the functionally graded material:

Here, is a constant material. The effects of the constant on the temperature, displacements, and thermal stresses are fully discussed in [6, 19]. For the sake of brevity, discussion of these effects is omitted here.

The basic piezothermoelastic equations for the th layer can be expressed as follows.

2.1. Conduction Equations and Equilibrium Equations

The heat conduction equation can be expressed as

where is the temperature in the th layer. While the equilibrium equation of the piezoelectric functionally graded cylindrical shell is

the electrostatics equation can be written as

The Constitutive equation can be

where

In the case of the stress-strain relations of functionally graded cylindrical shell ( to ), the terms and on the right side in (5) vanished. And , are elastic constants and thermoelastic coefficients.

The constitutive equations for the electric field are

2.2. Boundary and Continuity Conditions

Consider the following:

where and are the thermal dimensionless conductivity parameter in th and th layer, respectively. and are temperatures in the inner and outer surfaces.

For the simply supported panel, we assumed that the edges of the piezoelectric layer are electrically grounded, so the boundary conditions and continuous conditions between layers of displacements and stresses can be expressed as

To simplify the solving process of the 3D piezothermomechanical problem, the following dimensionless variables are introduced:

where is reference value of radius and , , , and are reference values of Young’s modulus, thermal expansion coefficient, thermal diffusion coefficient, thermal conductivity coefficient, and temperature, respectively, where are the stress components, are the strain components, (, , ) are the displacement components, are the elastic stiffness constants, are the electric displacement components, are the piezoelectric coefficients, are the dielectric constants, is the pyroelectric constant, and is the piezoelectric modulus.

The dimensionless displacement-expressed equilibrium equations can be expressed as

The dimensionless displacement-expressed electrostatic equation can be expressed as

Using Navier trigonometric series, solutions of (3), (13a) to (13c) and (14) which satisfy the ends displacements boundary conditions (11a) to (11b), (14) can be written as

where , ; , , , , and are unknown functions, Substituting trigonometric series (15a) to (15e) into (3) and (13a) to (13c) and (14), we can obtain

Furthermore, substituting the trigonometric series (15a) to (15e) into boundary and continuous conditions (10c) to (10e) and (11c) to (11f), we can obtain

Using the orthogonality of trigonometric functions, we can obtain

3. Theoretical Analysis

Based on the series solving method of ordinary differential equations, if the coefficient items of the ordinary differential equations were analytical at point and could be expressed as Taylor’s series of , the solutions of (17a)–(17c) and (18) can also be expressed as the following Taylor’s series [17]:

Substituting series (22) into (16) to (18), comparing the coefficient of , we can obtain the following recurrence equation:

where . From the recurrence equations (23) and (24a) to (24d), we can see that all coefficients , , , , and in Taylor’s series can be briefly expressed by , , , ; , ; , ; , and :

where the coefficients , ; ; ; and can be derived from the recurrence equations (23) and (24a) to (24d). Substituting the solved series coefficients (25a) to (25d) into Taylor series (22), we can obtain , , , and .

Where , ; , ; , ; , ; and are unknown constants and they can be determined by substituting solutions (24a)–(24d) into the boundary and continuous conditions (19a) to (19c) and (20a) to (20c).

Substituting solutions , , , , and into trigonometric series (15a)–(15e), we can obtain the steady-state piezothermomechanical solution for the finite cylindrical composite panel. Furthermore, using geometric equation and constitutive equation, the analytical solutions of thermo-mechanical stresses for the laminated cylinder can be obtained.

4. Numerical Results and Discussion

To illustrate the above analysis, we investigate the piezothermoelastic stresses in a functionally graded cylindrical shell with a piezoelectric layer on outside. The dimensions , , , and , the thickness of piezoelectric layer is 0.002, the thickness of each layer is assumed to be equivalent. The reference values of temperature, Young’s modulus, and thermal expansion coefficient are  K,  GPa, and  K⁻¹, respectively.

In this paper, the functionally graded finite cylinder made of mullite and molybdenum is analyzed. The thermoelastic properties of mullite and molybdenum are given as follows [11, 19]:

The material properties of piezoelectric layer (cadmium selenite) are summarized as

And it is assumed that the cylindrical shell is subjected to the following nonumiform thermal and mechanical loadings on its outer surfaces:

Figures 2 to 8 illustrate the numerical results of dimensionless temperature, displacement, and stresses of functionally graded cylindrical shell with piezoelectric layer.

(a)

(b)

Figure 2: Dimensionless temperature distribution: (a) at section , (b) .

Figure 2 shows the numerical results of dimensionless temperature distributions. Due to the graded variation of the material properties along the radial direction, the temperature reduced in a nonlinear form through thickness direction. For the thermal load is a function of sin, the maximum temperature in each layer is located along the axial .

Figure 3 shows the axial displacement; the magnitude of axial displacement becomes larger near the outer surface. Figures 4 to 12 show the numerical results of dimensionless radial, hoop, axial, and shear stresses, respectively. The radial stress increases from the inner surface to outer surface nonlinearly; the radial stress is compressive for the inner heating and pressure loads.

(a)

(b)

Figure 3: Dimensionless axial displacement distribution: (a) at section , (b) .

(a)

(b)

Figure 4: Distribution of dimensionless radial stress: (a) , , (b) .

Figure 5(a) shows the variation of circumferential stress on the thickness direction at point , . Due to the assumed thermomechanical loadings, the circumferential stress is compressive stress at the inner surface and is tensile stress at the outer surface. Because of the mismatch material properties, the circumferential stress is not continuous across the interfaces. Due to the assumed thermomechanical loading, the distribution of circumferential stress is nonuniform. It can be seen that maximum values of the circumferential stress occur at the same surface where the temperature loadings are applied. Similar distributions of circumferential stress can be obtained at other sections in the thickness direction. For the sake of brevity, the distributions in other thickness section are omitted here. Figure 6(a) shows the distributions of axial stress on the variation of thickness direction at and . in Figure 6(b), due to the assumed boundary condition and thermomechanical loading, the axial stress is tensile stress at the inner surface and varies from compressive to tensile stress, which is similar to circumferential stress.

(a)

(b)

Figure 5: Distribution of dimensionless circumferential stress: (a) , , (b) .

(a)

(b)

Figure 6: Distribution of dimensionless axial stress: (a) at point , . (b) .

Figures 7 to 9 show the distributions of three shear stresses at different sections. It is seen that the maximum values of all three shear stresses are smaller than those of normal stresses. Due to the assumed thermo-mechanical loadings, distributions of the three shear stresses show different distribution. According to the boundary condition and material properties, the maximum shear stresses and appear at the middle radius, while the shear stress varies from compressive stress to tensile stress along its length. The maximum shear stress appears at the middle along the length and the middle thickness. Figures 10, 11, and Figure 12 show the variation of electric potential and displacement along the thickness; the electric potential distribution is similar to the temperature distribution. The electric displacement varied nonlinearly.

Figure 7: Distribution of dimensionless shear stress ().

(a)

(b)

Figure 8: Distribution of dimensionless shear stress : (a) , (b) at point .0, .

Figure 9: Distribution of dimensionless shear stress ().

(a)

(b)

Figure 10: Distribution of dimensionless electric potential: (a) at section , (b) .

(a)

(b)

Figure 11: Distribution of dimensionless electric displacement : (a) at section , (b) .

(a)

(b)

Figure 12: Distribution of dimensionless electric displacement : (a) at section , (b) .

5. Conclusions

The piezothermoelastic analysis for a functionally graded cylindrical shell with piezoelectric layer is presented in this paper. The cylindrical shell is simply supported at four end edges and subjected to thermomechanical loadings on its outer surfaces. The piezoelectric layers are polarized along radial direction as a sensor. An exact 3D analytical solution is derived by using multilayered approximate approach and series solving method. Referred to the presented solution, other temperature boundary conditions can be solved by suitable forms of series solutions to satisfy both temperature and simply supported boundary conditions. It is also worth noting that the trigonometric series, which are used to separate variables, are only suitable for simply supported at four end edges. For other supported boundary conditions, another suitable series form must be considered to satisfy the related boundary conditions. For different functionally graded material, this analytical method can also be used.

Acknowledgment

This project is supported by the National Natural Science Foundation of China (Grant nos. 51375396 and 11202164).

References

I. F. PintoCarreia, C. M. MotaSoares, and C. A. MotaSoares, “Analysis of piezo-laminated axisymmetric shell: a semi analytical higher order model,” in Computational Methods For Shell and Spatial Structures, pp. 1–19, IASS-IACM, 2000.

M. Taya, A. A. Almajid, M. Dunn, and H. Takahashi, “Design of bimorph piezo-composite actuators with functionally graded microstructure,” Sensors and Actuators A, vol. 107, no. 3, pp. 248–260, 2003.

Y. Ootao and Y. Tanigawa, “Transient piezothermoelastic analysis for a functionally graded thermopiezoelectric hollow sphere,” Composite Structures, vol. 81, no. 4, pp. 540–549, 2007.

Y. Ootao and Y. Tanigawa, “Transient piezothermoelasticity for a cylindrical composite panel composed of cross-ply and piezoelectric laminae,” International Journal of Mechanical Sciences, vol. 44, no. 9, pp. 1861–1877, 2002.

Y. Ootao and Y. Tanigawa, “Transient piezothermoelasticity for a cylindrical composite panel composed of angle-ply and piezoelectric laminae,” International Journal of Solids and Structures, vol. 39, no. 23, pp. 5737–5752, 2002.

Y. Ootao and Y. Tanigawa, “Three-dimensional transient piezothermoelasticity in functionally graded rectangular plate bonded to a piezoelectric plate,” International Journal of Solids and Structures, vol. 37, no. 32, pp. 4377–4401, 2000.

Y. Ootao and Y. Tanigawa, “Three-dimensional transient piezothermoelasticity for a rectangular composite plate composed of cross-ply and piezoelectric laminae,” International Journal of Engineering Science, vol. 38, no. 1, pp. 47–71, 2000.

C.-Q. Chen and Y.-P. Shen, “Piezothermoelasticity analysis for a circular cylindrical shell under the state of axisymmetric deformation,” International Journal of Engineering Science, vol. 34, no. 14, pp. 1585–1600, 1996.

X.-H. Wu, Y.-P. Shen, and C. Chen, “An exact solution for functionally graded piezothermoelastic cylindrical shell as sensors or actuators,” Materials Letters, vol. 57, no. 22-23, pp. 3532–3542, 2003.

S. Kapuria and G. G. S. Achary, “A coupled consistent third-order theory for hybrid piezoelectric plates,” Composite Structures, vol. 70, no. 1, pp. 120–133, 2005.

Z. S. Shao, “Mechanical and thermal stresses of a functionally graded circular hollow cylinder with finite length,” International Journal of Pressure Vessels and Piping, vol. 82, no. 3, pp. 155–163, 2005.

D. J. Huang, H. J. Ding, and W. Q. Chen, “Piezoelasticity solutions for functionally graded piezoelectric beams,” Smart Materials and Structures, vol. 16, no. 3, pp. 687–695, 2007.

J. Yang and H. J. Xiang, “Thermo-electro-mechanical characteristics of functionally graded piezoelectric actuators,” Smart Materials and Structures, vol. 16, no. 3, pp. 784–797, 2007.

J. Ying and H. M. Wang, “Axisymmetric thermoelastic analysis in a finite hollow cylinder due to nonuniform thermal shock,” International Journal of Pressure Vessels and Piping, vol. 87, no. 12, pp. 714–720, 2010.

T. Roy, P. Manikandan, and D. Chakraborty, “Improved shell finite element for piezothermoelastic analysis of smart fiber reinforced composite structures,” Finite Elements in Analysis and Design, vol. 87, no. 12, pp. 714–720, 2010.

M. Zaman, Z. Yan, and L. Jiang, “Thermal effect on the bending behavior of curved functionally graded piezoelectric actuators,” International Journal of Applied Mechanics, vol. 2, no. 4, pp. 787–805, 2010.

A. Alibeigloo, “Thermoelastic solution for static deformations of functionally graded cylindrical shell bonded to thin piezoelectric layers,” Composite Structures, vol. 93, no. 2, pp. 961–972, 2011.

R. A. Alashti and M. Khorsand, “Three-dimensional thermo-elastic analysis of a functionally graded cylindrical shell with piezoelectric layers by differential quadrature method,” International Journal of Pressure Vessels and Piping, vol. 88, no. 5-7, pp. 167–180, 2011.

M. Jabbari, S. Sohrabpour, and M. R. Eslami, “Mechanical and thermal stresses in a functionally graded hollow cylinder due to radially symmetric loads,” International Journal of Pressure Vessels and Piping, vol. 79, no. 7, pp. 493–497, 2002.