Open access

Finite Element Analysis of Functionally Graded Piezoelectric Spheres

Written By

A. Ghorbanpour Arani, R. Kolahchi, A. A. Mosalaei Barzoki, A. Loghman and F. Ebrahimi

Published: 10 October 2012

DOI: 10.5772/52932

From the Edited Volume

Finite Element Analysis - Applications in Mechanical Engineering

Edited by Farzad Ebrahimi

Chapter metrics overview

2,971 Chapter Downloads

View Full Metrics

1. Introduction

A smart structure typically comprises of one or more active (or functional) materials. These active materials act in a unique way in which couple at least two of the following fields to provide the required functionality: mechanical, electrical, magnetic, thermal, chemical and optical. Through this coupling, these materials have the ability to change their shape, respond to external stimuli and vary their physical, geometrical and rheological properties. In modern technologies there has been an intense interest in FGPMs which are used in smart structures. It is well known that piezoelectric materials produce an electric field when deformed, and undergo deformation when subjected to an electric field. The coupling nature of piezoelectric materials has conducted wide applications in electro-mechanical and electric devices, such as electro-mechanical actuators, sensors and transducers. For example, piezoelectric actuators can be used to modify the shape of an airfoil, thereby reducing transverse vortices [1], or to maintain proper tension with overhead electrical wires on a locomotive pantograph [2].

For homogeneous piezoelectric media, problems of radially-polarized piezoelectric bodies were considered and solved analytically by Chen [3]. Sinha [4] obtained the solution of the problem of static radial deformation of a piezoelectric spherical shell and under a given voltage difference between these surfaces, coupled with a radial distribution of temperature from the inner to the outer surface. Ghorbanpour et al. [5] investigated the stress and electric potential fields in piezoelectric hollow spheres. Stress in piezoelectric hollow sphere under thermal environment was developed by Saadatfar and Rastgoo [6]. Dai and Wang [7] presented the thermo-electro-elastic transient responses in piezoelectric hollow structures. Dai and Fu [8] studied the electromagneto transient stress and perturbation of magnetic field vector in transversely isotropic piezoelectric solid spheres.

In-homogenity was considered in a number of studies. Elastic analysis of internally pressurized thick-walled spherical pressure vessels of functionally graded materials (FGMs) investigated by You et al. [9]. Analytical solution for a non-homogeneous isotropic piezoelectric hollow sphere was presented by Ding et al. [10]. Effect of material in-homogeneity on electro-thermo-mechanical behaviors of functionally graded piezoelectric rotating cylinder was considered by Ghorbanpour et al. [11]. Wang and Xu [12] studied the effect of material inhomogeneity on electromechanical behaviors of functionally graded piezoelectric spherical structures. Magnetothermoelastic problems of FGM spheres are studied by Ghorbanpour et al. [13].

Sladek et al. [14] derived Local integral equations for numerical solution of 3-D problems in linear elasticity of FGMs viewed as 2-D axisymmetric problems while the meshless local Petrov-Galerkin method was applied to transient dynamic problems in 3D axisymmetric piezoelectric solids with continuously non-homogeneous material properties subjected to mechanical and thermal loads by Sladek et al.[15]. They concluded that this method is promising for numerical analysis of multi-field problems like piezoelectric or thermoelastic problems, which cannot be solved efficiently by the conventional boundary element method.

Motivated by these ideas, new applications of piezoelectric sensors and actuators are being introduced and expanded for a number of geometric configurations. In this chapter, a hollow sphere composed of a radially polarized transversely isotropic piezoelectric material, e.g., PZT-4, which is subjected to mechanical and thermal loads, together with a potential difference induced by electrodes attached to the inner and outer surfaces of the annular sphere is considered. All mechanical, thermal and piezoelectric properties of the FGPM hollow sphere, except for the Poisson’s ratio, are assumed to depend on the radius r and expressed in terms of its power functions. Hence, the equation of equilibrium in the radially polarized form is reduced to a system of second–order ordinary differential equation and is solved analytically for four different sets of boundary conditions. Finally, the thermal stresses, electric potential and displacement distributions are shown for different material in-homogeneity Also, a three-dimensional finite element analysis of asymmetric closed and open spheres with different boundary conditions subjected to an internal pressure and a uniform temperature field has also been carried out using ANSYS software.

Advertisement

2. Electromechanical coupling

The subsequent characterization of electromechanical coupling covers the various classes of piezoelectric materials. Details with respect to definition and determination of the constants describing these materials have been standardized by the Institute of Electrical and Electronics Engineers [16]. Stresses σ and strains ε on the mechanical side, as well as flux density D and field strength E on the electrostatic side, may be arbitrarily combined as follows [17,18]

{σD}=[CEeeTε]{εE},E1

whereCE, ε, eand eT are the fourth-order elasticity tensor, the dielectric permittivity tensor, third order tensor of piezoelectric coefficient and transpose of it, respectively.

Assuming total strain tensor to be the sum of mechanical (M) and thermal (T) strains [19, 20]

ε=εM+εT,E2

where

ε={εrrεθθεζζεrζεζθεrθ},E3
εM=(u,rw,θrsinζ+cotζrv+urv,ζr+uru,ζr+v,rvrv,θrsinζ+w,ζrcotζrwu,θrsinζ+w,rwr), εT=(αrTαθTαζT000),E4

It is also noted that the electric field tensor E can be written in terms of electric potential ϕ as [21]

E=gradϕ.E5
Advertisement

3. Formulation for electrothermoelastic FGPM spheres

A hollow FGPM sphere with an inner radius ri and outer radiusro is considered. The sphere is subjected to an internal and external pressures Pi andPo, an electric potential φand a distributed temperature field T(r) (Fig. 1). It is assumed that, only the radial displacement Ur is nonzero and electric potential is the functions of radial coordinater, Thus

Ur=u(r),Uζ=Uθ=0,ϕ=ϕ(r).E6

Figure 1.

Hollow FGPM sphere subject to uniform temperature field, uniform internal pressure, uniform external pressure and applied voltage V.

The equilibrium equation of the FGPM sphere in the absence of body force and the Maxwell's equation for free electric charge density are [18, 22]

σrr,r+2(σrrσθθ)r=0,E7
Drr,r+2rDrr=0,E8

where σii(i=r,θ) is the stress tensor and Drris the radial electric displacement.

Also, the radial and circumferential strain and the relation between electric field and electric potential are reduced to

εrr=u,r,E9
εθθ=εζζ=ur,E10
Err=ϕ,r.E11

The constitutive relations of spherically radially polarized piezoelectric media and the component of radial electric displacement vector also can be written as [23, 24]

{σrrσθθσζζDrr}=[C11C12C13e11C21C22C23e12C31C32C33e13e11e12e1311]{εrrαrT(r)εθθαθT(r)εζζαζT(r)Err}E12

For transversely isotropic properties, when the concerned axis of rotation is oriented in the radial direction, the elasticity and piezoelectric coefficient tensors are summarized to [25]

C12=C13=C21,C22=C33,C32=C23, e12=e13.E13

It is appropriate to introduce the following dimensionless quantities as

σi=σiiC22(i=r,θ), ci=C1iC22(i=1,2,3),Ei=e1iE0(i=1,2,3), E0=C2211,U=urri, ξ=rri, η=rori, β=β1E0, Φ=ϕϕ0, ϕ0=riC2211, Dr=DrrE0.E14

Using the above dimensionless variables, Eqs. (7) and (8) can be expressed as

σr,ξ+2(σrσθ)ξ=0,E15
Dr,ξ+2Drξ=0.E16

Before substituting the component of the electric field in Maxwell’s equation, appropriate power functions for all properties are assumed as [26]

Γr=Γ0(ξ)γ,E17

in which Γrrepresents the general properties of the sphere such as the elastic, piezoelectric, dielectric coefficients and thermal conductivity, and Γ0 corresponds to the value of the coefficients at the outer surface. Substituting Eqs. (14) and (17) into Eq. (12), yields

{σrσθσζDr}=ξγ([C1C2C2E1C21C3E2C2C31E2E1E2E21]{U,ξξγαrT(ξ)UξξγαθT(ξ)UξξγαθT(ξ)Φ,ξ})E18
(18)
Advertisement

4. Electrothermoelastic analysis of FGPM spheres

The solution of Eq. (16) is

Dr=A1ξ2,E19

where A1is a constant. Substituting Eq. (19) into Eq. (18), we obtain

{σrσθ}=ξγ(([C1C2C2C21C3]+[E12E1E2E1E2E1E2E22E22])({U,ξξγαrT(ξ)UξξγαθT(ξ)UξξγαθT(ξ)})[E1E2](A1ξγ2))E20

In this study a distributed temperature field due to steady-state heat conduction has been considered. Using Eq. (17) for the thermal conductivity property, the heat conduction equation without any heat source is written in spherical coordinate as [22, 27]

1ξ2(K0ξγ+2T(ξ),ξ),ξ=0,E21
atξ=1T(ξ)=Ta,atξ=ηT(ξ),ξ+hT(ξ)=0,E22

where h is the ratio of the convective heat-transfer coefficient and K0 is the nominal heat conductivity coefficient. Integrating Eq. (21) twice yields

T(ξ)=B1γ+1ξγ1+B2,E23

Constants B1andB2are obtained using thermal boundary conditions which shown in Eq. (22).

Finally, substituting Eq. (20) and (24) into Eq. (15) yields the following non-homogeneous Cauchy differential equation

ξ22Uξ2+D1ξUξ+D2U=D4B1+D5B2ξ(1+γ)+D3A1ξ(1+γ)E24

where Di(i=1,...8)are defined in Appendix A.

The exact solution for Eq. (24) is written as follows

ug=K1exp(q1ξ)ug1+K2exp(q2ξ)ug2,E25
q1,q2=(1D1)±(D11)24D22.E26

The particular solution of the differential Eq. (24) may be obtained as

up=ξq1u1+ξq2u2,E27

where

u1=ξq2R(ξ)W(q1,q2),u2=ξq1R(ξ)W(q1,q2),E28

in which R(ξ) is the expression on the right hand side of Eq. (24) and W(ξ) is defined as

W(q1,q2)=|ug1ug2(ug1)(ug2)|.E29

Combining Eqs. (25)-(29) one can obtain the particular solution as

up=D3ξ1γ(q2+γ1)(q1+γ1)A1+D5ξ3+γ(3q2+γ)(3q1+γ)B2+D4ξ3(q23)(q13)B1.E30

The complete solution for Um in terms of the non-dimensional radial coordinate is written as

U=ug+up,E31

whereK1, K2and A1 are unknown constants. Substituting the displacement from Eq. (31) into Eq. (20) the radial and circumferential stresses are obtained as

σr(ξ)=((2C2+C1q1+2E1E2+E12q1)ξq11+γ)K1+((2C2+C1q2+2E1E2+E12q2)ξq21+γ)K2+((2C2+2E1E2+(3+γ)(E12+C1))D4ξ2(q23)(q13)+(2(C2+E1E2)αθ+(E12+C1)αr)ξγ1γ+1)B1+((2C2+2E1E2+(3+γ)(E12+C1))D5ξ2+2γ(3q2+γ)(3q1+γ)(2(C2+E1E2)αθ+(E12+C1)αr)ξ2γ)B2+((2C2+2E1E2+(1γ)(E12+C1))D3(q2+γ1)(q1+γ1)2(E1)ξ2)A1E32
σθ(ξ)=((C3+C2q1+E1E2q1+2E22+1)ξq11+γ)K1+((C3+2E22+E1E2q2+C2q2+1)ξq21+γ)K2+((C3+3C2+2E22+3E1E2+1)D4ξ2+γ(q23)(q13)+((C3+2E22+1)αθ+(C2+E1E2)αr)ξγ1γ+1)B1+((C3+2E22+1+(3+γ)(E1E2+C2))D5ξ2+2γ(3q2m+γ)(3q1m+γ)((1+C3+2E22)αθ+(C2+E1E2)αr)ξ2γ)B2+((C3+2E1E2+2E22+1+(1γ)(E1E2+C2))D3(q2+γ1)(q1+γ1)E2ξ2)A1E33

Substituting U from above into the last term of Eq. (18), (Φ,ξ) and combining with Eq. (19) and performing the integrating, electric potential is obtained as

Φ(ξ)=((2E2q1+E1)ξq1)K1+((2E2q2+E1)ξq2)K2+((γE1E12E2)D3ξ1γ(q2+γ1)(q1+γ1)(1γ)+ξγ1(1+γ))A1+((E1+23E2)D4ξ3(q23)(q13)+(2E2αθ+E1αr)ln(ξ)γ+1)B1+((3E1+γE1+2E2)D5ξγ+3(q2+γ1)(q1+γ1)(3+γ)(2E2αθ+E1αr)ξγ+1γ+1)B2+A2,E34

whereσr,σθ and Φ are radial stress, hoop stress and electric potential, respectively. Two sets of mechanical and electrical loading boundary conditions are considered in this investigation which in normalized form are written as

case 1:σr(1)=1, σr(η)=0,Φ(1)=0,Φ(η)=0E35
case 2:σr(1)=0, σr(η)=0,Φ(1)=1,Φ(η)=0E36

In case 1, the FGPM hollow sphere is subjected to an internal uniform pressure without any imposed electric potential and external pressure. However in this case the induced electric potential existed across the thickness. In this case, the sphere acts as an sensor. In the second case, an electrical potential difference is applied between the inner and outer surfaces of the sphere without any internal and external pressures. In this case, the sphere acts as an actuator.

For the above mentioned cases 1, and 2 the system of linear algebraic equations for the constantsK1, K2, A1and A2 of the Eqs. (32), (33) and (34) can be written in the following from

[K1K2A1A2]=inv([m11m12m13m14m21m22m23m24m31m32m33m34m41m42m43m44])×[b1b2b3b4],E37

where the mij and bi (i,j=1,...4) are defined in Appendix B

Advertisement

5. Numerical results and discussion

5.1. Analytical solution

The numerical results are showing the variation of stresses, electric potential and displacement across the thickness of the FGPM sphere for different material inhomogenity parameterγ. Presented results are for the two cases of different boundary conditions with aspect ratioη=1.3. The plots in these figures correspond toTi=323K andTo=298K. The piezoelectric material PZT-4 has been selected because of its technical applications. Mechanical and electrical properties of piezoelectric material, PZT-4 are tabulated in Table 1[28].

PropertyC11C12C22C23e11e12e11αr0ρ
PZT-4115 Gpa74.3 Gpa139 Gpa77.8 Gpa15.1 C/m2-5.2 C/m23.87e-9 F/m2 e - 5 1/K2 e -6 1/K7500 kg/m3

Table 1.

Mechanical, electrical and thermal properties for PZT-4

5.1.1. Case 1

Results of the first case are illustrated in Figs. 2 to 5. Radial stresses for different material in-homogeneity parameters γare shown in Fig. 2. Radial stresses satisfy the mechanical boundary conditions at the inner and outer surfaces of the FGPM sphere. The maximum absolute values of radial stresses belong to a material identified by in-homogeneity parameter γ=1.5 the minimum absolute values of which belong toγ=1.5. In this case there is no imposed electric potential however, the induced electric potentials for different material in-homogeneity parameters γ are shown in Fig. 3. Electric potentials satisfy the electrical boundary conditions at the inner and outer surfaces of the FGPM sphere. It is also obvious that higher induced electric potentials belong to higher absolute values of compressive radial stresses. In this case circumferential stresses shown in Fig. 4 are highly tensile throughout thickness and except for the material γ=1.5 their maximum values located at the inner and their minimum values located at the outer surfaces of the sphere. Displacements are illustrated in Fig. 5 for all material properties. Displacements are positive throughout the thickness and they smoothly decrease from their maximum value at the inner surface to their minimum value at the outer surface of the FGPM sphere. Maximum values of displacements belong toγ=1.5 and minimum values belong toγ=1.5

Figure 2.

Case 1: Distributions of the radial stress for different values of γ.

Figure 3.

Case 1: Distributions of the circumferential stress for different values of γ.

Figure 4.

Case 1: Distributions of the electric potential for different values of γ.

Figure 5.

Case 1: Distributions of the radial displacement for different values of γ.

5.1.2. Case 2

Results of this case are illustrated in Figs. 6 to 9. In this case the imposed electric potential satisfies the electrical boundary conditions at the inner and outer surfaces of the sphere. The maximum electric potentials belong to γ=1.5 the minimum values of which belong toγ=1.5. In this case there is no applied pressure at the inner and outer surfaces of the sphere however the induced compressive radial stresses satisfy the free mechanical boundary conditions. The maximum absolute values of the induced compressive radial stresses belong to the same maximum value of electric potential. Circumferential induced stresses are compressive throughout thickness for different material in-homogeneity parametersγ. However, for negative parameters γ the minimum values of circumferential stresses located at the inner surface while for positive parametersγ their minimum values located at the outer surface of the FGPM sphere. The induced displacement is negative across the thickness for all material parameters. Their minimum values located at the inner and their maximum values located at the outer surfaces of the FGPM sphere. It is interesting to compare the induced radial and circumferential stresses in this case with the residual stresses locked in the sphere during the autofrettage process of spheres made of uniform material. One might conclude that by easily imposing an electric potential there is no need to autofrettage these vessels.

Figure 6.

Case 2: Distributions of the radial stress for different values of γ.

5.2. Validation

The results of this investigation are validated with the recently published paper by Wang and Xu [12] which is shown in Figs. 10 and 11. There are a very good agreement among the results and the only small differences are due to thermal stresses which are not considered by Wang and Xu.

Figure 7.

Case 2: Distributions of the circumferential stress for different values of γ.

Figure 8.

Case 2: Distributions of the electric potential for different values of γ.

Figure 9.

Case 2: Distributions of the radial displacement for different values of γ.

Figure 10.

Case 2: Comparison of the radial stress distributions with Ref. [12] for homogeneous piezoelectric hollow sphere.

Figure 11.

Case 2: Comparison of the electric potential distributions with Ref. [12] for homogeneous piezoelectric hollow sphere.

5.3. Finite element solution

In order to develop the one-dimensional solution to a three-dimensional approach finite element analysis of a sphere subjected to an internal pressure and a uniform temperature field has been carried out using ANSYS finite element software. A three-dimensional element identified by solid 191 is selected because it is an appropriate element for the FGPM structures. Sphere has been divide into eight layers by a controlled mesh system along radius and the mechanical, electrical and thermal properties are functionally defined according to power law Eq. (17) forγ=1.5. A controlled mesh in which very fine elements are located at the supports where stress concentration existed is employed in this method. However farther from the supports a coarser mesh is dominated. In this work, two cases for sphere are considered as follows:

5.3.1. Three- dimensional sphere

In this case, consider a sphere with two asymmetric simply supported boundary conditions on the outer surface of the sphere as shown in Fig. 12. For this boundary condition dimensionless effective stresses versus normalized radius at two cross sections (i.e. A-A and B-B) are depicted in Fig. 13. Section A-A is selected to pass through supported point on the outer surface of the sphere and section B-B is an arbitrary section as shown in Fig. 12. It can be seen from this figure that the maximum effective stress for the above mentioned sections occur at the inner surface of the sphere and the effective stresses are decreasing with increasing radius forη=1.3. Total dimensionless displacement versus dimensionless radius for η=1.3 at two cross sections of A-A and B-B are demonstrated in Fig. 14. As can be seen from this figure the maximum displacement occur at the inner surface of the sphere and displacement value is decreasing with increasing dimensionless radius so that for section A-A, the zero value of displacement at the outer surface satisfies the boundary condition at this point.

Figure 12.

A schematic of asymmetric thick-walled sphere with simply-simply supported.

Figure 13.

Effective stress distribution along the radius of asymmetric closed sphere.

Figure 14.

Total displacement distribution along the radius of asymmetric closed sphere.

5.3.2. Three- dimensional open sphere

The geometry and loading condition as well as its boundary conditions are shown in Fig. 15. Three different boundary conditions are considered in this case. These boundary conditions are clamped-clamped, clamped-simply and simply-simply supported respectively.

The solution obtained by the software clearly indicates the most critical region of the sphere. In the most critical region normalized effective stress distribution and the total displacements are plotted in Figs. 16 and 17 along normalized radius at all node points for the above mentioned three boundary conditions. Fig 16 shows that in general the effective stresses are decreasing along radius to an absolute minimum and then increasing to their maximum values located at the outer surface of the vessel. For simply-simply supported boundary condition this absolute minimum is located near the outer surface of the vessel, however for the clamped-clamped condition it is nearly at the middle surface of the vessel. For the clamped-simply supported condition this minimum is somewhere between the previous two cases.

It has been found that the magnitude of effective stresses at all node points are higher for the clamped-clamped condition and are lower for the simply-simply supported condition. It can be observed from Fig. 17 that the maximum displacements for the three boundary conditions are located at the inner surface and they are decreasing to zero value at the outer surface of the sphere. It is also found that the displacement curve for simply-simply support condition is higher than other boundary conditions.

Figure 15.

The three-dimensional finite element model for open sphere subjected to internal pressure with clamped-clamped boundary conditions.

Figure 16.

Effective stresses distribution along the radius of the open sphere with different boundary conditions.

Figure 17.

Total displacement distribution along the radius of the open sphere with different boundary conditions.

Advertisement

6. Conclusions

In this research, the electro-thermo-mechanical behavior of radially polarized FGPM hollow sphere was investigated. An analytic solution technique was developed for the electro-thermo-mechanical problem, where stresses were produced under combined thermomechanical and electrical loading conditions. Variation of normalized stresses, electric potential and displacement of four sets of boundary conditions for different material in-homogeneity parameters γ were plotted against dimensionless radius. In general, radial stresses and electric potentials satisfy the mechanical and electrical boundary conditions at the inner and outer surfaces of the FGPM sphere. It was concluded that higher absolute values of compressive radial stresses are associated with the higher induced electric potentials throughout the thickness in all cases. It was found that the induced radial and circumferential stresses of an imposed electric potential is similar to the residual stresses locked in the sphere during the autofrettage process of these vessels. Therefore, one might concluded that by easily imposing an electric potential there is no need to autofrettage these vessels. It was interesting to see that the compressive circumferential stresses due to an external pressure were very similar to the induced circumferential stresses resulted from imposing an electric potential. Moreover a three-dimensional finite element analysis of an asymmetric sphere subjected to an internal pressure and a uniform temperature field has been carried out using ANSYS software. In this study closed and opened spheres with different boundary conditions were considered. The finite element analysis indicated that the values of effective stress and total displacement at all node points along the thickness of the open sphere were the highest and lowest for the clamped-clamped condition, respectively.

Appendix A

D1=γ(C1+E12)+2C1+2E12C1+E12,D2=γ(2C2+2E1E2)+2C2+2E1E222C34E22C1+E12,D3=2E2C1+E12 D4=(γ+1)(2C2αθ+C1αr+E1(2E2αθ+E1αr))+2(2C2αθ+C1αr)+2(E1E2)(2E2αθ+E1αr)((1+C3)αθ+C2αr)(C1+E12)(γ+1)D5=γ(4C2αθ2C1αr2E1(2E2αθ+E1αr))2(2C2αθ+C1αr)(2E1E2)(2E2αθ+E1αr)+(2(1+C3)αθ+2C2αr)(C1+E12)

Appendix B

mij(i,j=1,...4) are the coefficient ofK1, K2, A1 and A2 which defined as
m11=2C2+C1q1+2E1E2+E12q, m12=2C2+C1q2+2E1E2+E12q2, m13=(2C2+2E1E2+(1γ)(E12+C1))D3(q2+γ1)(q1+γ1)2(E1)
m21=(2C2+C1q1+2E1E2+E12q1)ηq11+γ, m22=(2C2+C1q2+2E1E2+E12q2)ηq21+γ, m23=(2C2+2E1E2+(1γ)(E12+C1))D3(q2+γ1)(q1+γ1)2(E1)η2, m31=(2E2q1+E1), m33=(γE1E12E2)D3(q2+γ1)(q1+γ1)(1γ)+1(1+γ), m41=(2E2q1+E1)ηq1 , m42=(2E2q2+E1)ηq2
m43=(γE1E12E2)D3η1γ(q2+γ1)(q1+γ1)(1γ)+ηγ1(1+γ), m41=0 , m24=0, m34=1 , m44=1
bi(i=1,2)are correspond to cases 1 and 24 which denoted as follows:b1=σr(1), b1=σr(ξ)

References

  1. 1. P. Destuynder, A few remarks on the controllability of an aeroacoustic model using piezo-devices, Int. J. Holnicki-Szulc. (1999) 53-62.
  2. 2. H.W. Jiang, F. Schmid, W. Brand, G.R. Tomlinson, Controlling pantograph dynamics using smart technology, Int. J. Holnicki-Szulc. (1999) 125-132.
  3. 3. W.Q. Chen, Problems of radially polarized piezoelastic bodies, Int.J. solids struct. 361998199843174332
  4. 4. D.K. Sinha, Note on the radial deformation of a piezoelectric, polarized spherical shell with a symmetrical distribution, J. Acoust. Soc. 34 (1962) 1073-1075.
  5. 5. A. Ghorbanpour, S. Golabi, M. Saadatfar, Stress and electric potential fields in piezoelectric smart spheres, J. Mech. Sci. Tech. 20 (2006) 1920-1933.
  6. 6. M. Saadatfar, A. Rastgoo, Stress in piezoelectric hollow sphere under thermal environment, J. Mech. Sci. Tech, 22 (2008) 1460-1467.
  7. 7. H.L. Dai, X. Wang, Thermo-electro-elastic transient response, Int. J. solids struct. 42 (2005) 1151-1171.
  8. 8. H.L. Dai, Y.M. Fu, Electromagnetotransient stress and perturbation of magnetic field vector in transversely isotropic piezoelectric solid sphere, Mater.Sci. Eng. B 129200620068692
  9. 9. L.H. You,J.J. Zhang, X.Y. You, Elastic analysis of internally pressurized thick-walled spherical pressure vessels of functionally graded materialsInt. J. Press. Vess. Pip. 8220052005347354
  10. 10. H.J. Ding, H.M. Wang, W.Q. Chen, Analytical solution for a non-homogeneous isotropic piezoelectric hollow sphere, Arch.Appl. Mech. 73200320034962
  11. 11. A. Ghorbanpour Arani, R. Kolahchi, A.A. Mosallaie Barzoki, Effect of material in-homogeneity on electro-thermo-mechanical behaviors of functionally graded piezoelectric rotating cylinder, J. Appl. Math. Model. 35 (2011) 2771-2789.
  12. 12. H.M. Wang, Z.X. Xu, Effect of material inhomogeneity on electromechanical behaviors of functionally graded piezoelectric spherical structures, Comput.Mater. Sci. 4820102010440445
  13. 13. A. Ghorbanpour, M. Salari, H. Khademizadeh, A. Arefmanesh, Magnetothermoelastic problems of FGM spheres, Arch. Appl. Mech. (2010) 189-200.
  14. 14. V. Sladek, J. Sladek, Ch. Zhang, Transient heat conduction analysis in functionally graded materials by the meshless local boundary integral equation method, Comput. Mater. Sci. 28 (2003) 494-504.
  15. 15. J. Sladek, V. Sladek, P. Solek, A. Saez, Dynamic 3D axisymmetric problems in continuously non-homogeneous piezoelectric solids, Int. J. solids struct. 45 (2008) 4523-4542.
  16. 16. Institute of Electrical and Electronics Engineers. Standard on Piezoelectricity, Std (1978) 176-1978 IEEE, New York.
  17. 17. Y.C. Fungn, Foundations of Solid Mechanics, Prentice-Hall, New York, 1965.
  18. 18. H.F. Tiersten, Linear Piezoelectric Plate Vibrations, Plenum Press, New York, 1969
  19. 19. A. Manonukul, F.P.E. Dunne, D. Knowles, S. Williams, Multiaxial creep and cyclic plasticity in nickel-base superalloy C263, Int. J. Plasticity, 21 (2005) 1-20.
  20. 20. H. Martin, ELASTICITY Theory, Applications and Numerics, Elsevier Inc, London, 2005.
  21. 21. H.J. Ding, W.Q. Chen, Three Dimensional Problems of Piezoelasticity, Nova Science, New York, 2001.
  22. 22. M. Sadeghian, H. Ekhteraei Toussi, Axisymmetric yielding of functionally graded spherical vessel under thermo-mechanical loading, Comput. Mater. Sci. 50 (2011) 975-981.
  23. 23. A. Salehi-Khojin, N. Jalili, A comprehensive model for load transfer in nanotube reinforced piezoelectric polymeric composites subjected to electro-thermo-mechanical loadings. Compos: Part B. 39 (2008) 986-998.
  24. 24. Zh. Li, Ch. Wang, Ch. Chen, Effective electromechanical properties of transversely isotropic piezoelectric ceramics with microvoids, Comput.Mater. Sci. 2720032003381392
  25. 25. N. Jalili, Piezoelectric-Based Vibration Control from Macro to Micro/Nano Scale Systems, Boston, 2010.
  26. 26. M.H. Babaei, Z.T. Chen, Analytical solution for the electromechanical behavior of a rotating functionally graded piezoelectric hollow shaft, Arch.Appl. Mech. 7820082008489500
  27. 27. H.L., Dai, L. Hong, Y. M., Fu, X. Xiao, Analytical solution for electro-magneto-thermo-elastic behaviors of a functionally graded piezoelectric hollow cylinder, J. Appl. Math. Model. 34 (2010) 343-357.
  28. 28. M. Saadatfar, A.S. Razavi, Piezoelectric hollow cylinder with thermal gradient, J. Mech. Sci. Tech. 23 (2009) 45-53.

Written By

A. Ghorbanpour Arani, R. Kolahchi, A. A. Mosalaei Barzoki, A. Loghman and F. Ebrahimi

Published: 10 October 2012