Document Type : Research Article
Authors
^{1} Faculty of Mechanical and Mechatronics Engineering, Shahrood University of Technology, Shahrood, Iran
^{2} Faculty of Mechanical and Mechatronics Engineering, Shahrood University of Technology, Semnan, Iran
Abstract
Keywords
ElectroElastic Analysis of Finite Length FGPM Cylinders Subjected to Electromechanical Loading Using FirstOrder Electric Potential Theory
M. Parhizkar Yaghoobi, M. Ghannad^{*}
Faculty of Mechanical and Mechatronics Engineering, Shahrood University of Technology, Semnan, Iran
KEYWORDS 

ABSTRACT 
2D ElectroElastic FGPM Cylinders Piezoelectric Shear Deformation Theory FirstOrder Electric Potential Theory

In this study, an analytical solution is presented based on the voltage feedback control method for the twodimensional electroelastic static response of functionally graded piezoelectric material (FGPM) cylinders. Using firstorder shear deformation theory as well as firstorder electric potential theory and applying the energy method, a differential equations system is extracted, which is solved as a classical eigenvalue problem. The results show the significant impact of heterogeneity on the electromechanical behavior of the cylinders. Furthermore, control gain affects the electric potential and electromechanical behavior of the head where the voltage is applied. The present research also introduces an analytical solution with no limitation to specific conditions in cylinder heads and without any need for convergence check. Moreover, the results show that any changes in cylinder head conditions affect the behavior of FGPM cylinders. The results were compared with those from the finite element method (FEM), leading us to a reasonable agreement. 
Nowadays, researchers are looking for materials that, in addition to being endowed with such features as constructional suitableness like low weight and high strength among others, would have other desirable properties with a view to a new revolution in structures, mechanical devices and equipment. Piezoelectric materials are substances in which electric fields are deformed or their deformation makes electric field in them. Piezoelectric materials can be considered as smart materials because the electromechanical behavior of these materials enables them to be used as sensors, actuators or active controllers of the structure. In order to improve the structures, the idea of functionally graded materials (FGMs) or heterogeneous materials was developed, which led to improved distribution properties changes. Such materials can prove to have better continuity and integrity properties than conventional composite materials, and their stresses and changes can have a continuous state which leads to an increase in the material strength. Therefore, researchers have focused on the behavior of structures with these materials by studying their behavior from different perspectives. Among the various structures, shells are of particular importance mainly due to the high efficiency of their behavior against load bearing. Many examples of these types of structures, such as dome roofs, tankers, fluid storage tanks, nuclear reactors, projectiles, missiles, etc., are found in artificial systems made by humans. These samples are also found in natural systems, for example, in skulls, shells and bones among others. These samples represent merely a limited portion of the extensive use of shells in natural and artificial systems [1, 2]. Among diverse shell geometries, cylinders have attracted researchers the most due to their high rate of usage and easy construction. That is why a section is dedicated to them in most curriculums and we review it in the relevant literature that follows.
Dai et al. in 2010 [3] obtained a heterogeneous electrothermomagnetoelastic analysis using the plane elasticity theory (PET). They considered the heterogeneity of the distribution of properties in the cylinder in an exponential state. In their study, the heat transfer and electrostatic equilibrium were considered as onedimensional and in radial direction. It was also assumed that the cylinder is in the uniform magnetic field located in the direction of its axis. That led to the creation of heterogeneous differential equations. By solving these equations and applying boundary conditions, they studied the electrothermoelastic behavior of the cylinder under electrothermoelastic loading. In the same year, Ghannad and ZamaniNejad [4] analyzed the homogeneous doubleheaded clamped cylinder under internal pressure using FirstOrder Shear Deformation Theory (FSDT) and studied the shear stress in the cylinder. They showed that in areas far from the boundaries of the two cylinder’s heads the shear stress tends to zero, causing stresses and displacements to be only a function of cylinder radius, while at the boundary it is a function of the radius and length of the cylinder. They also observed that at points far from the boundary, there is a good agreement between the theory of shear deformation and the plane elasticity. Ghorbanpour Arani et al. [5] presented an exact solution for an axisymmetric functionally graded piezoelectric (FGP) rotating disk with constant thickness subjected to an electric field and thermal gradient. Using PET, a nonhomogeneous second order differential equation is derived and solved.
In 2012, Ghannad and ZamaniNejad [6], using the theory of plane elasticity and the definition of structural relationships in general state (plane stress and plane strain), obtained a complete solution to the heterogeneous thick wall cylinders that were under internal and external pressure. They considered the heterogeneity of properties for the elasticity modulus as the exponential state and assumed that the Poisson coefficient was constant. They also showed that in order to increase or decrease stress and displacement, positive or negative values should be taken for the heterogeneity constant. Also, they [7] analyzed the heterogeneous thick wall cylinders that were subjected to internal and external pressure, using the shear deformation theory.
In 2013, Dai and Jiang [8] investigated the electromagnetothermoelastic behavior of a FGP solid cylinder. The solid cylinder was located in a uniform magnetic field and subjected to a thermal load, electric excitation, and pressure in the outer layer. Using long cylinder assumption (plane strain condition) and regardless of the end effects of the FGPM solid cylinder, they presented an analytical solution for the problem and studied the heterogeneity of properties effects on electromagnetothermoelastic behavior. Kargarnovin et al. [9] presented the exact planar solution for electroelastic analysis of FGP structures under thermoelectromechanical loadings. The material properties were assumed to vary exponentially along the thickness. Exploiting the potential functions for stress and induction, the solution was obtained in a closed form manner. The influence of material heterogeneity was examined against the electroelastic reaction of the FGP media.
In 2014 and 2015, Ghannad and Gharooni [1011], using the firstorder shear deformation theory and higherorder, achieved the displacements field and stresses in the heterogeneous cylinder whose properties changed as an exponential function under internal and external pressure. Khorshidi and Bakhsheshy analyzed the Free Natural Frequency of an FG Composite Rectangular Plate Coupled with Fluid using Rayleigh–Ritz Method [12]. In this research, the influence of the main parameters of the problem on wet natural frequencies is discussed. Dynamic elasticity solution for a clamped, laminated cylindrical shell with two orthotropic layers bounded with a piezoelectric layer and subjected to impulse load distributed on inner surface was presented by Saviz [13].
Jabbari and Aghdam in 2015 [14] conducted the thermoelastic analysis of threelayered cylinders, consisting of two piezoelectric homogeneous cylinders and a heterogeneous cylinder. Using PET, the governing differential equations of the problem at plane strain state and based on onedimensional thermal hypothesis are solved. The impact of heterogeneity of properties and control of the mechanical behavior of the cylinder by piezoelectric layers was observed. A free vibration analysis of the FG rectangular nanoplates is investigated by Khorshidi et al. [15]. Nonlocal elasticity theory is employed and the effect of transverse shear deformation and rotary inertia is investigated by the exponential shear deformation theory. The effect of various parameters including nonlocal parameters and power law indexes on natural frequencies is studied.
Atrian et al., in 2015 [16] obtained a thermoelastic analysis of piezoelectric heterogeneous cylinders that were axialasymmetrically loaded. The electrothermoelastic behavior of the heterogeneous thick wall cylinder was investigated in two dimensions and in the radial and circumferential directions. The heterogeneity of the properties was modeled as an exponential function that was changed in the radial direction. In this study, the effect of heterogeneity and the effect of electrothermoelastic loading on heterogeneous cylinder behavior was investigated. The active forced vibration control of circular plates coupled with piezoelectric layers and laminated composite rectangular plate resting on linear line support is presented by [17] and [18], respectively. The results of this research show that by applying control force at the plates, the amplitude of transverse deformation will be controlled and reduced.
Loghman and Parsa [19] using PET, extracted a closed form solution for the electromagnetothermoelastic response to a thick doublelayered cylinder made from a homogeneous interlayer and an FGPM outer layer. Ghannad and Parhizkar Yaghoobi in 2015 and 2017 presented the twodimensional thermoelasticity solution of homogeneous and heterogeneous thick wall cylinders that were placed under their thermomechanical loading on the inner and outer layers for different boundary conditions in two cylindrical heads [2021]. Using the energy method and the first order shear deformation theory and assuming linear variations of the thermal field along the thickness, they obtained the system of governing equations. In their study, the effects of different boundary conditions on thermoelastic behavior of homogeneous and heterogeneous cylinders were studied. Arefi et al. [22] presented an analytical method for the threedimensional free vibration analysis of an FG cylindrical shell integrated by two thin FGP layers. The FSDT and quadratic forms formulation are used for the displacement field and electric potential, respectively. Using the energy method and considering the von Karman nonlinear straindisplacement relations, the equations of motion are derived. In this research the effects of the nonhomogeneity of materials and geometry on natural frequency are studied. The sizedependent free vibration analysis of a sandwich nanoplate is studied by Arefi and Zenkour [23]. In this paper, the governing equation are derived using FSDT and a cosine form for electric potential distributions. The influence of important parameters of the problem such as the applied electric, magnetic potentials, etc. on the vibration characteristics of the sandwich nanoplate are discussed.
The exact closedform solution for sound radiation of vibrating circular plates coupled with piezoelectric layers is presented by [24]. In this paper, the effect of boundary conditions, piezoelectric thickness, and the piezoelectric layer on the acoustical parameters was examined and discussed in details. Using differential quadrature method (DQM), static analysis of FG solid circular and annular plates integrated with piezoelectric layers was performed by Alibeigloo [25]. It was assumed that the plate has various edge boundary conditions and whose material properties vary in an exponential law across the thickness. The effects of material heterogeneity, geometry parameter, and various boundary conditions on the response of FG structures were studied. Mehditabar et al. [26] conducted research on the thermoselasticity analysis FGP rotating hollow cylindrical shell subjected to dynamic loads. The material properties were varied in radial direction of cylinder by a power function. Using DQM and finite difference methods, the solution of boundary/initial value equations of the problem are provided. Free vibrations of composite rectangular piezoelectric nanoplate using modified shear deformation theories are analyzed by Khorshidi et al. [27]. The governing equations are derived based on the energy method and using the exponential shear deformation and trigonometric shear deformation theories. The effect of small scale and geometry parameters on the natural frequency of composite rectangular piezoelectric nanoplate is studied.
Parhizkar Yaghoobi et al. [28] in 2018 performed a longitudinal piezoelectric heterogeneous cylinder (plane strain state) and thin piezoelectric heterogeneous discs (plane stress states) analysis of stress and active control analysis. In their research, they assumed that the internal surface of the structure is like a sensor and its external surface is like an actuator, and is based on the method of solving the Lame. They finally provided an exact onedimensional solution (in radial direction). The effect of heterogeneity constant and control coefficient factor on the electrothermomechanical behavior are studied. Arefi and Bidgoli [29] presented electroelastic analysis of the piezoelectric doubly curved shells resting on Winkler's foundation subjected to applied voltage. Using the energy method and the first order shear deformation theory, they derived the system of governing equations. The influence of important parameters such as applied electric potential and Winkler’s parameter were investigated on the electroelastic behaviors of the structure.
Creep stress and strain in Laminated Composite Pressure Vessels are analyzed by Ghasemi and Hosseinpour [30]. The classical lamination theory is used to derive the governing equations as a secondorder equation to determine the radial, circumferential, axial, and effective stresses in the cylinder wall. The effect of orientation of fibers on creep strain distribution in the wall of a cylinder is studied. The twodimensional electroelastic analysis of FG carbon nanotubes reinforced composite cylindrical pressure vessels integrated with piezoelectric facesheets as sensor and actuator is extracted by Arefi et al. [31]. The governing equations of the problem are derived based on the principle of virtual work. Using the classical eigenvalueeigenvector method the problem is solved and the electroelastic behavior of the structure is studied. Also, the effect of characteristics and distribution of porosity on the electroelastic behavior of laminated vessels with piezoelectric facesheets based on higher order modeling is analyzed by Arefi et al. [32].
The electroelastic analysis of FG carbon nanotubes reinforced composite cylindrical shells with piezoelectric layers based on thirdorder shear deformation theory is done by Mohammadi et al. [33]. Using the classical method, the analytical solution for the problem is obtained and the results are discussed. The thermomagnetomechanical longterm creep behavior of threephase nanocomposite cylinder is investigated by Ghasemi and Hosseinpour [34]. The results show that the values of creep strain and radial displacement in the magnetic field are lower than without a magnetic field, for two layups. Also, the effects of the magnetic field in the creep behavior of threephase laminated composite cylindrical shells are analyzed by [35]. It was found that temperature increase in the magnetic field is less effective on the increased values of creep strain and radial displacement. The fluidstructure interaction of vibrating composite piezoelectric plates is investigated by Khorshidi and Karimi [36]. The governing equations are derived based on Hamilton’s principle and using the exponential shear deformation theory and a cosine form for electric potential distributions. The natural frequencies of the fluidstructure system are computed using the Galerkin method and the influence of the main variables on natural frequencies are discussed.
Considering the literature background of the subject, it seems that the problem of twodimensional solution of the electroelasticity static response of FG thick wall cylinders in the general case of boundary conditions with a voltage feedback control method has not been taken into account, and electroelastic analysis has been achieved only in limited cases with simplification in the boundary conditions. For this reason, using the firstorder shear deformation theory, firstorder electric potential theory and applying the energy method the governing equations to heterogeneous cylinders that are under electromechanical loading are obtained under arbitrary electromechanical boundary conditions with a voltage feedback control method. The governing equilibrium equations are arranged in the form of a set of ordinary differential equations, which are solved by analytical method. By applying arbitrary electromechanical boundary conditions with a voltage feedback control method in the two cylinder’s heads, the analytical constants of ordinary differential equations system are obtained and the numerical results for the cylinder are shown under two different boundary conditions. These results are compared with the results of FEM.
The analyzed cylindrical shell is axial symmetry in terms of geometry, properties, loading, and boundary conditions, and it is possible to neglect changes in the circumferential direction ( ); therefore, the functions representing the mechanical and electrical behavior are only a function of the radius r and the x axis of the cylinder.
According to what was stated, the displacement field in this case is as follows.
(1) 
, and are the displacement components in radial, circumferential, and axial directions, and and are a function of cylinder radius and axis.
The studied cylindrical shell has a continuous wall and is made of heterogeneous and transversely isotropic materials, whose material property in the radial direction is different from the two other directions. The heterogeneity of properties in this research is considered as a power state similar to [7], and using the relationship (2) the heterogeneous modeling of properties for analysis is presented.
(2) 
In the above relations is the internal radius of the cylinder and is the property of the material in the inner layer of the cylinder, which can have mechanical and electrical properties such as modulus of elasticity, density, Poisson ratio, piezoelectric constants and dielectric stability constants. in Eq. (2) above is the heterogeneity constant for property; the heterogeneity constant can take positive and negative real values. It is noteworthy that denotes the homogeneous materials.
To use the firstorder shear deformation theory, as shown in Fig. 1, the distance between each point of the cylindrical shell to its axis of symmetry r ranges is expressed as its distance to the midline of the cylinder z plus the distance from the middle layer to the cylinder R symmetry axis. This work causes the coordinate change the differential component from r to z (dr,(r,θ,x)→dz,(r,θ,x)) [7].
Fig. 1. Schematic of geometry and cross section of FGPM cylinder.
Given what has been said, we have.
(3) 
Also, for z and x variations we have:
(4) 
The displacement field for axial symmetric cylinder is described using the first order shear deformation theory as follows [7].
(5) 
In the above relation, , and are the components of the displacement field along radial, circumferential and axial directions; also, according to the shear deformation theory, and have zeroorder and oneorder components of radial displacement, respectively and and are zeroorder and oneorder components of axial displacement, respectively which are only a function of axial coordination.
Also, using the firstorder electric potential we have [37].
(6) 
Where and are respectively zeroorder and first order electric potentials, which are the only functions of the xaxis coordination.
The electric field and strain field are extracted using gradient relationships in the cylindrical coordinates [37] and the relations (5) and (6).
(7) 
, and are radial, circumferential, and axial normal strains respectively. Also, is the radialaxial shear strain.
(8) 
, and are the components of the electric field in the radial, circumferential and axial directions respectively. The constitutive equations are expressed as follows [38].
(9) 
In Eq. (9), , , and are radial, circumferential and axial normal stresses respectively, is the radialaxial shear stress, and D_{z} and D_{x} are the components of the electric displacement vector in radial direction and axial directions respectively. Also, c_{11}, c_{12}, c_{22}, c_{23}, and c55 are the elastic constants e_{11}, e_{12}, and e_{35} are piezoelectricstress constants and and are the dielectric constants in the abovementioned relationship.
For heterogeneity modeling in properties, the relationship (2) or power distribution is used, and the properties used in the structural relation (9) are modeled in an exponential state. So we will have.
(10) 
In the relation (10) c_{11i}, c_{12i}, c_{22i}, c_{23i}, and c_{55i} are the elastic constants values, e_{11i}, e_{12i}, e_{35i} are the piezoelectricstress constants and and are the values of dielectric constants in the heterogeneous inner cylinder and n is the heterogeneity constant.
Fig. 2. FGPM cylinder cross section under loading on inner and outer surfaces.
Figure 2 shows the crosssection of a heterogeneous cylinder that is electromechanically loaded at its internal and external surfaces. As can be seen, the pressure P_{i} and the electric charge density Q_{i} are applied to the cylinder inner surface and the pressure P_{o} and the electric charge density Q_{o} are applied to the outer surface of the cylinder.
Mechanical and electrical resultants are defined as follows.
(11) 

(12) 

(13) 

(14) 

(15) 
Based on the principle of virtual work, energy changes in the electromechanical potential of the structure are equal to the variation of work due to external forces applied to the structure δw=δU [37, 3940] where U is the total electromechanical energy of the whole body and W is the whole work due to the applied of the pressure and electrical charge density on the inner and outer surfaces of the cylinder. The energy of electromechanical potential of the whole body and the work resulting from these forces are calculated as follows.
(16) 

(17) 
In Eq. (17) the physical interpretation of the terms are the electrical work done after applying the electrical surface charge density on the inner and outer surfaces of the cylinder respectively.
For work and energy changes, we have:
(18) 
By substituting the Eqs. (79) in (18), the application of the principle of virtual work and performing mathematical operations and simplifying, the system of governing equations is obtained as follows.
(19) 
Also, for boundary conditions we have:
(20) 
In this section an analytical solution for the equation system (19) is presented. As shown in the system of Equations (19), mechanical and electrical resultants are more numerous than the equations. Therefore, it is not possible to solve the system of equations (19) on the basis of resultants, and the equations (19) must be expressed on the basis of displacement field (4 unknown components) and electric potential (2 unknown components) so that they would be solved. Thus, using the relations (11) to (15) (definition of mechanical and electrical resultants) and the relations (5) to (10) and applying mathematical operations, the governing equations system is rewritten as equations (21) based on the displacement field and electric potential components. So we will have:
(21) 
Where , , and are the coefficients matrices and and are respectively electromechanical force quasivector and electromechanical quasivector, whose nonzero components are defined and presented in Appendix A. Moreover, the "T" superscript shows the transpose of quasivector. K_{9} and K_{10} are solved constants created by integrating the first and fifth equations of governing Eq. (19), and F_{3} to F_{6} are presented in Appendix A. Solving this equation’s system consists of two parts; particular solving section and general solving section .
(22) 
As can be seen from Eq. (21), the heterogeneous part of the differential equation system is a firstorder polynomial. Therefore, solving a particular part is considered as a firstorder polynomial with nondefinite coefficients. Therefore, we have:
(23) 
These indefinite (unknown) coefficients become definite (known) as follows.
(24) 
To solve the general problem, the equation system should be solved in the following form.
(25) 
The general solution is considered in the following form.
(26) 
The general solution form of (26) is valid for Equation (21) when the coefficients matrix C does not have zero rows or columns. For this purpose, in deriving equations (21), the first and fifth equations of the system of equations (19) have been integrated. Furthermore, in defining the electromechanical quasivector, instead of zeroorder axial displacement component and zeroorder component of electric potential, their derivatives have been used.
By substituting the relation 26 in equation 25, the eigenvalue problem is obtained.
(27) 
To this end, the determinant of the Eq. (27) must be set to zero.
(28) 
By solving the above eigenvalue problem whose eigenvalue equation is an eightorder one and by finding eigenvalue vectors, the following general solution is obtained.
(29) 
Therefore, the solution of the system of equations (21) is obtained.
(30) 
Given that (30) yields the obtained solution for the derivatives of the functions and , in order to obtain these functions, it is necessary to integrate the proposed solution for the derivation of these functions, in which case the constants K_{11} and K_{12} are obtained. Now these boundary constants are obtained by applying boundary conditions in the two cylinder heads.
In this section, an adaptive voltage feedback control method is applied to supply control voltage for activating the voltage in the upper cylinder head. It is assumed that the voltage of lower cylinder head ( ) that is without electrical surface charge is measured by the embedded voltage sensor and fed back to the upper cylinder head as a control voltage by the controller. The control voltage used for operating the upper cylinder head is proportional to the voltage of lower cylinder head with a control gain (control factor) C_{f} and may be expressed as follows [14, 28, 41].
(31) 
The mechanical boundary conditions in both cylinder’s heads could be of essential conditions. For example, the equation for the cylinder with two clamped heads in these conditions is as follows:
(32) 
The mechanical boundary conditions on the two cylindrical heads can be a combination of essential and natural conditions, for example, for a cylinder with lower clamped head and upper free head, the mechanical boundary conditions in both heads of the cylinders are expressed as follows:
(33) 
Therefore, by applying the boundary conditions on the two cylindrical heads – 6 conditions in one head (4 mechanical conditions and 2 electrical conditions) and 6 conditions on the other head (4 mechanical conditions and 2 electrical conditions) – K_{1}K_{12} constants can be achieved. It is noteworthy that the process of mathematical solution using coding in MAPLE 13 software is achieved.
For a numerical study, a heterogeneous cylinder with an inner and outer diameter of 40 mm and 60 mm respectively and a length of 800 mm is considered. The internal layer of this heterogeneous cylinder is made of PZT4 and has the mechanical and electrical properties mentioned in the inner layer in accordance with Table 1 [42].
In order to study the effect of heterogeneity of properties and control gain on the electromechanical behavior of the cylinder, the heterogeneity constant (n) for the values n = 2, 1, 0, 1 and 2, and the control gain (C_{f}) for the values C_{f} = 10, 5, 0, 5, 10 and 10, are studied. It is noteworthy that the inner layer of this cylinder under the surface charge density is and its outer layer is under pressure of 6 bar. This heterogeneous cylinder is studied in two cases: (a) two clamped heads; and (b) one lower clamped head and the upper free head. The static electric field that causes depolarization (E_{de}) often varies between and , depending on the temperature of the applied electrical field [43]. The maximum amount of electric field imposed by C_{f} = 10 inside the cylinder is about that it is less than the lower limit of the assumed depolarizing electric field. The select ranges for the control factor must be chosen based on the amount of electrical surface charge density at the inner and outer surfaces of the cylinder.
In order to validate the results, the voltage of the lower end of cylinder (b) was assumed at zero. Figs. 3 and 4 show the distribution of radial displacement and electric potential in the FG cylinder (b) for validation. The results are predicted based on two analytical methods (present research and PET [28]) and FEM.
Table 1. Mechanical and electrical properties of the inner layer of FGPM cylinder.
Properties of PZT4 

Elastic constants (GPa) 
c_{11} 
115 
c_{12} 
74 

c_{22} 
139 

c_{23} 
78 

c_{55} 
25.6 

Piezoelectricstress constants (C/m^{2}) 
e_{11} 
15.6 
e_{12} 
5.2 

e_{35} 
12.7 

Dielectric Constants (10^{8}C/Vm) 
∊_{11} 
0.562 
∊_{22} 
0.646 
Figure 5 shows the distribution of electric potential, radial and axial displacement in FG cylinders. Because of the existence of electrical and mechanical boundary conditions at the two ends of the cylinder and the interaction of these conditions caused by the piezoelectricity of the cylinder, the distribution field of the displacement field and the electric potential near the boundaries are nonuniform.
Furthermore, in cylinder (a) or a cylinder with two clamped heads, where electrical boundary conditions in two cylinder heads (Eq. (31)) are nonuniform, the contour of the radial displacement behavior distribution is not symmetric with respect to the middle of the cylinder length; It is rather slightly asymmetric.
As can be seen in Fig. 5, cylinder (b), unlike cylinder (a), has only one mechanical constrained head and the upper head is free; this causes the maximum amount of displacements values to increase substantially in the cylinder. The greatest increase is related to the axial displacement of the cylinder, where the amount of maximum axial displacement has increased about 40 times. The radial and axial displacement changes along the thickness (except for axial displacement near the boundaries) are negligible, but the electric potential along the thickness has considerable variation.
Fig. 3. Distribution of radial displacement in the middle length of FGP cylinder (n=1).
Fig. 4. Distribution of electric potential in the middle length of FGP cylinder (n =1).
Fig. 5. Distribution of radial and axial displacement and electric potential in FG cylinders (n=2, C_{f} =10).
In fact, the electric potential sensitivity to the cylinder's thickness is greater than the sensitivity of the displacement field to the thickness. Maximum displacement values occur in or near the boundary areas (cylinder’s head). Therefore, when considering the requirements of the design, engineers should pay attention to the boundary areas of the structure.
Figures 6 and 7 show the radial and axial displacement distribution in the inner and middle layer of heterogeneous cylinders respectively. Reduction of heterogeneity causes a significant increase in the distribution of displacement values. But it has no effect on the location of the emergence of maximum displacement values, and the maximum displacement values occur in the border regions or near the boundaries. As can be seen, the nonconstrained state of one cylinder head causes the maximum displacement to increase, which is a larger increase in axial displacement and is more significant.
Although in the initial review it seems that Fig. 7 is not true for cylinder (b), by considering the mechanical conditions related to the two ends of the cylinder (b) compared with cylinder (a), the axial displacement behavior in the longitudinal direction is valid. Cylinder (a) is a clampedclamped ends one while cylinder (b) is a clampedfree ends one. As can be seen in Fig. 5, a little away from the lower end of the cylinder (clamped condition), the supportive effects of the lower boundary on the axial displacement behavior vanishes and the axial displacement along the longitudinal distribution is almost linear.
Fig. 6. Distribution of radial displacement in the inner layers of FG cylinders (C_{f} =0). 
Fig. 7. Distribution of axial displacement in the middle layer of FG cylinders (C_{f} =0).
Figures 8 and 9 show the effect of control gain on the radial and axial displacement in the inner and middle layer of heterogeneous cylinders respectively. Generally, the effect of control gain on the radial and axial displacement can be ignored and it doesn’t cause a remarkable effect on radial and axial displacement; however, this effect is more visible in the heads of cylinder where the control voltage is applied.
Figure 10 shows the distribution of electric potential in the outer layer of the heterogeneous cylinder (em and ϵm are the average of piezoelectricstress and dielectric constants in the inner layer of FG cylinder respectively). An increase in the heterogeneity constant causes the distribution of the electric potential in the heterogeneous cylinder decrease; in other words, the increase in heterogeneity constant leads to a decrease in the intensity of electric potential changes in the heterogeneous cylinder. The maximum and minimum values of electric potential occurred near the cylinder heads (boundary areas). Therefore, designers and engineers should pay attention to this issue when using and designing piezoelectric heterogeneous cylinders.

Fig. 8. Result of control gain on radial displacement in the inner layer of FG cylinders (n = 2).
Fig. 9. Result of control gain on axial displacement in the middle layer of FG cylinders (n = 2). 
Fig. 10. Distribution of electric potential in the outer layer of FGPM cylinders (C_{f}=0).
Figure 11 shows the effect of control gain on the electric potential in the outer layers of FG cylinders. The control gain has significant effects on the electric potential behavior profile especially in the upper head of cylinder where the voltage is applied. In order to better consider and evaluate the control gain effect on electric potential in the middle of the middle layer, the ζ is defined based on Eq. (33),which φ(0,L/2) and φ^{1.25}(0,L/2) are the electric potential in the middle of the middle layer with arbitrary control gain value and C_{f}=1.25, respectively.
Figure 12 illustrates the effect of control gain on electric potential behavior of FG cylinders in z=0 and x=L/2. Fig. 10 predicts the same behavior for controlability FG cylinder structure like [35] by highlighting the effect of cylinder heads or cylinder boundary conditions. It shows changes in the boundary conditions from clampedclamped ends to clampedfree ends (mechanical boundary conditions) cause changes at about 10 percent in the controlability of FG cylinder.
(34) 
Figure 13 shows the stress distribution inside the FGPM cylinder along the thickness in the middle of the cylinder length. Shear stress is negligible compared to other stresses in distant areas from boundaries and can be neglected.
Fig. 11. Result of control gain on electric potential in the outer layer of FG cylinders (n = 2).
Fig. 12. Result of control gain on the electric potential behavior of FG cylinders in z=0 (n = 2).
Circumferential stress has a dominant effect among the stresses and has great effect on determining the effective stress behavior calculated based on the von Mises failure theory [21]. Also, the axial stress caused by the existence of the free cylinder head in the cylinder (b) has tension state and compression state distribution.
(35) 
Figure 14 shows the effective stress distribution. The maximum amount of effective stress occurs near the boundary areas, which increases with any increase in the heterogeneity constant. As can be seen, the intensity of the heterogeneity constant effect on the cylinder of type (a) or of the cylinder with two constrained heads is greater and causes more changes.
Fig. 13. Distribution of stresses in the middle of FGPM cylinders length (C_{f}=0, n=2).
Fig. 14. Distribution of effective stress in the middle layers of FGPM cylinders (C_{f}=0).
Figure 15 shows the effect of control gain on the effective stress in the middle layer of heterogeneous cylinders. Generally, the effect of control gain on the effective stress is negligible except in the upper cylinder head that the voltage is applied on it. Special attention should be considered for the head of cylinder with applied voltage when it was designed.
Figure 16 shows the distribution of the electric field along the axial direction in heterogeneous cylinders. The heterogeneity constant has no observable effect on the distribution of the electric field along the axial direction. As it is seen, the largest distribution of the electric field occurs in the vicinity of the boundaries or within the boundaries themselves, requiring more attention to these areas.
(36) 
Figure 17 shows the of control gain on electric field in the middle layers of heterogeneous cylinders, which is calculated according to the Eq. (35). As it is seen, the effect of control gain on electric field.is remarkable at the upper head or applied voltage head. Therefore, the designers should pay special attention to this head.
Table 2 shows the results of solving analytical and finite element problems using the ANSYS 12.0 software. As can be seen, there is a good consensus between the predicted results of the two methods that showed in Figs. 6, 7, 10 and 14.
Fig. 15. Result of control gain on effective stress in the middle layer of FG cylinders (n = 2).
Fig. 16. Distribution of axial component of electric field in the middle layers of FGPM cylinders (C_{f} =0).
Fig. 17. Result of control gain on electric field in the middle layers of FG cylinders (n = 2).
Table 2. The numerical results achieved for analytical solution and FE solution of FGPM cylinders (x=L/2, C_{f} =0).
Cylinder (b) 
Cylinder (a) 




n = 2 
n= 0 
n = 2 
n = 2 
n = 0 
n = 2 



0.000612 
0.000931 
0.001344 
0.000469 
0.000706 
0.001007 
Analytical 


0.000626 
0.000954 
0.001379 
0.000481 
0.000726 
0.001039 
FEM 
U_{r } (mm) 
Inner Layer 
2.24 
2.41 
2.54 
2.49 
2.75 
3.07 
Difference (%) 


0.000598 
0.000919 
0.001340 
0.000430 
0.000654 
0.000942 
Analytical 


0.000604 
0.000935 
0.001379 
0.000435 
0.000667 
0.000977 
FEM 
U_{r } (mm) 
Outer layer 
0.99 
1.71 
2.83 
1.15 
1.95 
3.58 
Difference (%) 


1.2406 
1.8936 
2.7429 
0.9705 
1.4482 
2.0423 
Analytical 


1.4647 
2.1750 
3.0646 
1.2748 
1.8776 
2.6202 
FEM 
σ_{eff } (MPa) 
Inner layer 
15.30 
12.94 
10.50 
23.87 
22.87 
22.06 
Difference (%) 


2.0583 
1.4110 
0.9194 
1.5751 
1.0450 
0.6552 
Analytical 


1.9697 
1.3792 
0.9325 
1.3526 
0.8540 
0.4857 
FEM 
σ_{eff } (MPa) 
Outer layer 
4.50 
2.31 
1.40 
16.45 
23.37 
34.89 
Difference (%) 


644.34 
1079.73 
1703.42 
373.45 
660.44 
1102.63 
Analytical 


587.14 
993.18 
1576.90 
368.53 
660.74 
1114.90 
FEM 
φ (V) 
Inner layer 
9.74 
8.71 
8.02 
1.33 
0.05 
1.10 
Difference (%) 


775.98 
1258.94 
1927.67 
632.10 
1043.48 
1636.92 
Analytical 


768.52 
1227.40 
1833.60 
656.11 
1068.10 
1639.20 
FEM 
φ (V) 
Outer layer 
0.97 
2.57 
5.13 
3.66 
2.30 
0.14 
Difference (%) 


In this study, Active control and the twodimensional analytical solution of piezoelectric FGM cylinders under mechanical and electrical loading in their internal and external layers and different boundary conditions in two cylinder’s heads Using FirstOrder Electric Potential Theory and a voltage feedback control method was presented. A numerical study was carried out for the two states of boundary conditions: (a) two clamped heads; and (b) one lower clamped head and the upper free head, and the effects of heterogeneity, control gain, boundary conditions, and electromechanical loading were studied. Then, the results of this study were compared with the results obtained from FEM, which showed a good agreement between the findings. The results of this study can be summarized as follows.
Appendix
Nonzero components of matrix are calculated as follows:
(A1) 
Which in the above relations K_{s} is the shear correction coefficient and is equal to 5/6 [7]. Also, for the nonzero components of matrix we have:
(A2) 
And for matrix nonzero components, we have:
(A3) 
Quasivector nonzero components are calculated as follows:
(A4) 
References
[1] Ventsel E. and Krauthammer T., 2001. Thin Plates and Shells: Theory, Analysis, and Applications, Marcel Dekker Company.
[2] Tzou H.S., Lee H.J. and Arnold S., 2004. Smart materials, precision sensors/actuators, smart structures, and structronic systems, Mechanics of Advanced Materials and Structures, 11, pp. 367393.
[3] Dai H.L, Hong L., Fu Y.M. and Xiao X., 2010. Analytical Solution for Electromagnetothermoelastic Behaviors of a Functionally Graded Piezoelectric Hollow Cylinder, Applied Mathematical Modelling 34, pp. 343–357.
[4] Ghannad M. and ZamaniNejad M., 2010. Elastic Analysis of Pressurized Thick Hollow Cylindrical Shells with ClampedClamped Ends, Mechanika, 85, pp. 1118.
[5] Ghorbanpour Arani A., Jafarzadeh Jazi A., Abdollahian M., Mozdianfard M.R., Mohammadimehr M. and Amir S., 2011. Exact Solution for Electrothermoelastic Behaviors of a Radially Polarized FGPM Rotating Disk, Journal of Solid Mechanics, 3, pp. 244257.
[6] Ghannad M. and ZamaniNejad M., 2012. Complete Elastic Solution of Pressurized Thick Cylindrical Shells Made of Heterogeneous Functionally Graded Materials, Mechanika, 18, pp. 640649.
[7] Ghannad M. and ZamaniNejad M., 2012. Elastic Analysis of Heterogeneous Thick Cylinders Subjected to Internal or External Pressure Using Shear Deformation Theory, Acta Polytechnica Hungarica, 9, pp. 117136.
[8] Dai H.L. and Jiang H.J., 2013. Analytical study for electromagnetothermoelastic behavior of a functionally graded piezoelectric solid cylinder, Mechanics of Advanced Materials and Structures, 20, pp. 811818.
[9] Kargarnovin M.H., Hashemi R. and Emami A.A., 2013. Electroelastic Analysis of FG Piezoelectric Structures under ThermoElectroMechanical Loadings, Mechanics of Advanced Materials and Structures, 20, pp. 1127.
[10] Ghannad M. and Gharooni H., 2014. Displacements and Stresses in Pressurized Thick FGM Cylinders with Exponentially Varying Properties Based on FSDT, Structural Engineering and Mechanics, 51, pp. 939953.
[11] Ghannad M. and Gharooni H., 2015. ThermoElastic Analysis of ClampedClamped Thick FGM Cylinders by Using ThirdOrder Shear Deformation Theory, Latin American Journal of Solids and Structures, 12, pp. 10241041.
[12] Khorshidi K. and Bakhsheshy A., 2014. Free natural frequency analysis of an FG composite rectangular plate coupled with fluid using Rayleigh–Ritz method, Mechanics of Advanced Composite Structures, 1(2), pp. 131143.
[13] Saviz M.R., 2015. Elasticity Solution of a Laminated Clamped Cylindrical Shell with Piezoelectric Layer under Dynamic Load, Mechanics of Advanced Materials and Structures, 22, pp. 503519.
[14] Jabbari M. and Aghdam M. B., 2015. Asymmetric Thermal Stresses of Hollow FGM Cylinders with Piezoelectric Internal and External Layers, Journal of Solid Mechanics, 7, pp. 327–343.
[15] Khorshidi K., Asgari T. and Fallah A., 2015. Free vibrations analysis of functionally graded rectangular nanoplates based on nonlocal exponential shear deformation theory, Mechanics of Advanced Composite Structures, 2(2), pp. 7993.
[16] Atrian A., JafariFesharaki J. and Nourbakhsh S.H., 2015. ThermoElectromechanical Behavior of Functionally Graded Piezoelectric Hollow Cylinder under NonAxisymmetric Load, Applied Mathematics and Mechanics (English Edition), 36, pp. 939954.
[17] Khorshidi K., Rezaei E., Ghadimi A. A. and Pagoli M., 2015. Active vibration control of circular plates coupled with piezoelectric layers excited by plane sound wave, Applied Mathematical Modelling, 39(34), pp. 12171228.
[18] Khorshidi K., Balali M. and Ghadimi A. A., 2015. Control Forced Vibrations Of Laminated Composite Rectangular Plate Resting On Linear Line Support, Modares Mechanical Engineering, 15(9), pp. 95104.
[19] Loghman A. and Parsa H., 2016. Closed Form Solution for ElectroMagnetoThermoElastic Behaviour of DoubleLayered Composite Cylinder, Journal of Solid Mechanics, 8, pp. 3144.
[20] Ghannad M. and Parhizkar Yaghoobi M., 2015. A Thermoelasticity Solution for Thick Cylinders Subjected to ThermoMechanical Loads Under Various Boundary Conditions, International Journal of Advanced Design and Manufacturing Technology, 8, pp. 112.
[21] Ghannad M. and Parhizkar Yaghoobi M., 2017. 2D Thermo Elastic Behavior of a FG Cylinder under Thermomechanical Loads Using a First Order Temperature Theory, International Journal of Pressure Vessels and Piping, 149, pp. 7592.
[22] Arefi M., Karroubi R. and IraniRahaghiand M., 2016. Free vibration analysis of functionally graded laminated sandwich cylindrical shells integrated with piezoelectric layer, Applied Mathematics and Mechanics (English Edition), 37(7), pp. 821834.
[23] Arefi M. and Zenkour A. M., 2017. Sizedependent free vibration and dynamic analyses of piezoelectromagnetic sandwich nanoplates resting on viscoelastic foundation, Physica B, 521, pp. 188197.
[24] Khorshidi K. and Pagoli M., 2016. Analytical solution for sound radiation of vibrating circular plates coupled with piezoelectric layers, Mechanics of Advanced Composite Structures, 3(2), pp. 8998.
[25] Alibeigloo A., 2018. Thermo Elasticity Solution of Functionally Graded, Solid, Circular, and Annular Plates Integrated with Piezoelectric Layers Using the Differential Quadrature Method, Mechanics of Advanced Materials and Structures, 25, pp. 76684.
[26] Mehditabar A., Rahimi G.H. and Tarahhomi M.H., 2018. ThermoElastic Analysis of a Functionally Graded Piezoelectric Rotating Hollow Cylindrical Shell Subjected to Dynamic Loads, Mechanics of Advanced Materials and Structures, 25, pp. 10681079.
[27] Khorshidi K., Siahpush A. and Fallah A., 2017. Electromechanical free vibrations analysis of composite rectangular piezoelectric nanoplate using modified shear deformation theories, Journal of Science and Technology of Composites, 4(2), pp. 151160.
[28] Parhizkar Yaghoobi M., Ghaffari I. and Ghannad M., 2018. Stress and active control analysis of functionally graded piezoelectric material cylinder and disk under electrothermomechanical loading, Journal of Intelligent Material Systems and Structures, 29(5), pp. 924937.
[29] Arefi M. and MohammadRezaei Bidgoli E., 2018. Electroelastic displacement and stress analysis of the piezoelectric doubly curved shells resting on Winkler's foundation subjected to applied voltage, Mechanics of Advanced Materials and Structures, 0(0), pp. 114.
[30] Ghasemi A. R. and Hosseinpour K., 2018. Creep Strain and Stress Analysis in Laminated Composite Pressure Vessels, Mechanics of Advanced Composite Structures, 5(2), pp. 141147.
[31] Arefi M., Mohammadi M., AmirAhmadi S. and Rabczuk T., 2019. FSDT electroelastic analysis of FGCNTRC cylindrical threelayered pressure vessels with piezoelectric facesheets, ThinWalled Structures, 144, p. 106320.
[32] Arefi M., Mohammadi M. and Rabczuk T., 2019. Effect of Characteristics and Distribution of Porosity on ElectroElastic Analysis of Laminated Vessels with Piezoelectric FaceSheets Based on Higher Order Modeling, Composite Structures, 225, p. 111085.
[33] Mohammadi M., Arefi M and AmirAhmadi S., 2020. TwoDimensional ElectroElastic Analysis of FGCNTRC Cylindrical Laminated Pressure Vessels with Piezoelectric Layers Based on ThirdOrder Shear Deformation Theory, Journal of Pressure Vessel Technology, 142(2), p. 021304.
[34] Ghasemi A. R. and Hosseinpour K., 2018. ThermoMagnetoMechanical LongTerm Creep Behavior of ThreePhase NanoComposite Cylinder, Composites Science and Technology, 167, pp. 7178.
[35] Hosseinpour K. and Ghasemi A. R., 2019. Effects of Magnetic Field in Creep Behavior of ThreePhase Laminated Composite Cylindrical Shells, Mechanics of Advanced Composite Structures 6(1), pp. 5156.
[36] Khorshidi K., and Karimi M., 2020. FluidStructure Interaction of Vibrating Composite Piezoelectric Plates Using Exponential Shear Deformation Theory, Mechanics of Advanced Composite Structures, 7(1), pp. 5969
[37] Benjeddou A., 2000. Advances in Piezoelectric Finite Element Modeling of Adaptive Structural Elements: a Survey, Computers and Structures, 76, pp. 347363.
[38] Wang X. and Zhong Z., 2003. The General Solution of Spherically Isotropic Magnetoelectroelastic Media and Its Applications, European Journal of Mechanics A/Solids, 22, pp. 953969.
[39] Benjeddou A. and Andrianarison O., 2005. A Thermopiezoelectric Mixed Variational Theorem for Smart Multilayered Composites, Computers and Structures, 83, pp. 12661276.
[40] Ghaffari I., Parhizkar Yaghoobi M. and Ghannad M., 2018. Complete Mechanical Behavior Analysis of FG Nano Beam under NonUniform Loading Using NonLocal Theory, Materials Research Express, 5(1), p. 015016.
[41] Li J. and Li F., 2019. Active Control of Thermal Buckling for Plates Using a Temperature Feedback Control Method, Smart Materials and Structures, 28(4), p. 045001.
[42] Yang J., 2006. The Mechanics of Piezoelectric Structures, World Scientific.
[43] Ye Z.G., 2008. Handbook of Dielectric, Piezoelectric and Ferroelectric Materials: Synthesis, Properties and Applications, CRC Press.