Electro elastic analysis of a pressurized thick-walled functionally graded piezoelectric cylinder using the first order shear deformation theory and energy method/Hermetisku storasieniu funkciskai kokybisku pjezoelektriniu cilindru elektriskai tampri analize remiantis pirmos eiles slyties deformacijos teorija ir energetiniu metodu.

Rahimi, G.H. ; Arefi, M. ; Khoshgoftar, M.J. 等

1. Introduction

Piezoelectrics are new groups of material which can be used as a sensor or actuator in electromechanical systems. These materials can exchange the mechanical deformations into electric potential. Conversely, the electric potential can be exchanged to the mechanical deformation. The piezoelectric sensors or actuators may be designed as many structural elements such as beam, plate or cylindrical shell. The piezoelectric analysis of a functionally graded piezoelectric (FGP) cylindrical pressure vessel is studied in the present paper. A brief review of functionally graded material (FGM) and functionally graded piezoelectric material (FGPM) are performed in introduction.

One of the most applicable structures in the mechanical engineering is the shells. In this study, the cylindrical shell structure is considered. Lame [1] studied the exact solution of a thick walled cylinder under inner and outer pressures. It was supposed the cylinder to be axisymmetric and isotropic. Piezoelectric property has been discovered by Pierre and Jacques Curie in Paris (1888). Shear deformation theory has been proposed by Naghdi and Cooper [2]. The application of first order shear deformation theory for an isotropic cylinder has been proposed by Mirsky and Hermann [3]. In the 1980's one Japanese group of material scientists created new class of materials. Properties of this material are varying continuously and gradually in terms of coordinate system components.

The researches on the thermal and vibration analysis of functionally graded materials have been started in the first years of decade 1990 [4]. Displacement and stress analysis of a functionally graded cylinder under the thermal and mechanical loads is performed analytically by Jabbari et al [5]. It was supposed the material properties are varying as a power function in terms of radial coordinate system. Chen et al [6] investigated the mechanical and electrical analyses of a spherical shell. Liu et al [7] proposed an analytical model for free vibration analysis of a cylindrical shell under mechanical and electrical loads. Mindlin's theory is investigated for this analysis and a sinusoidal function is used for simulation of the electric potential distribution. Wu, Jiang and Liu [8] investigated the elastic stability of a FG cylinder. They employed the shell Donnell's theory to derive the strain-deformation relations. Exact solution of a FGP clamped beam is investigated by Shi and Chen [9]. Peng-Fei and Andrew [10] studied the piezoelectric analysis of a cylindrical shell. Lu et al [11] studied the exact solution of a FGP cylinder under bending. Dai et al [12] analyzed the electromagnetoelastic behavior of FGP cylindrical and spherical pressure vessels. Babaei and Chen [13] presented exact solution of an infinitely long magneto elastic hollow cylinder and solid rotating cylinder that is polarized and magnetized radially. They supposed the cylinder to be orthotropic and investigated the effect of angular velocity on the hoop and radial stresses. Jabbari et al [14] investigated the thermoelastic behavior of a FG cylinder under the thermal and mechanical loads. Khoshgoftar et al [15] investigated the thermoelastic analysis of a FGP cylindrical pressure vessel. They supposed that all mechanical and electrical properties are varying as a power function. This mentioned work was the last comprehensive thermoelastic analysis of a FGP cylindrical shell using the plane elasticity theory. The present paper develops the previous paper significantly using the shear deformation theory and proposes an analytical formulation for a comprehensive analysis of a FGP cylinder. The proposed formulation is validated at the regions that are adequate far from two ends of the cylinder with the previous plane elasticity theory. Thermo-elastic analysis of a functionally graded cylinder is investigated analytically by Arefi and Rahimi [16]. They used the first order shear deformation theory (FSDT) for thermoelastic analysis of a FG structure. The achieved results are compared with those results that have been derived using the plane elasticity theory. Thermoelastic vibration and buckling characteristics of a functionally graded piezoelectric cylindrical shell is investigated analytically by Sheng and Wang [17, 18]. First order shear deformation theory is investigated for simulation of the deformations in structure. Electric potential is considered as a quadratic function along the thickness. The Hamilton's principle and Maxwell's equation are considered for solution of the problem. The critical values of axial load, temperature and voltage are investigated for different boundary conditions. Analytical solution for electromagneto thermoelastic behaviors of a functionally graded piezoelectric hollow cylinder under a uniform magnetic field and subjected to thermoelectromechanical loads is investigated by Dai et al [19]. They presented advantageous of material non homogeneity for design optimization of electro mechanical structures and systems.

As mentioned in the literature, there is not reported a comprehensive analysis about electroelastic analysis of a FGP structure by considering the whole nonzero piezoelectric coefficient and in the general state (last study devoted to Dai et al [19] that the cylinder is analyzed using one dimensional method. The proposed method in that paper has not ability to consider other components of strain and piezoelectric coefficients). Therefore, the present paper employs the comprehensive electroelastic formulation for the analysis of a FGP cylinder under inner pressure as an applied problem. The whole elastic, piezoelectric and dielectric coefficients in constitutive and Maxwell's equations are considered to be nonzero. This subject has been disregarded in the previous papers [15, 19]. Although the mentioned problem at the end of the paper can be considered as a plane strain problem, the present method of solution has enough capability to solve the problem in the general state with considering the whole piezoelectric coefficients. This advantageous condition cannot be understood in the plane strain method.

2. Formulation

In the present paper, the FSDT is employed to simulate the deformations. Based on this theory, deformation of every layer of the cylinder is decomposed into deformation of the middle surface and rotation about outward axis of the middle surface [3]. In order to better understand this theory, it is necessary to expand Lame's solution for a cylindrical pressure vessel. Based on the Lame's theory, symmetrical distribution of the radial displacement u may be obtained as follows [16, 20-24]

u = [c.sub.1]r + [c.sub.2]/r (1)

where r is the radius of every layer of the cylinder. In the general state, this distance can be obtained in terms of the radius of middle surface R and distance of every layer with respect to middle surface [rho]. By substitution of r into Lame's solution (Eq. (1)) and applying the Taylor expansion, Eq. (1) may be obtained as a function of [rho] as follows [16]

r = R + [rho] [right arrow] u = [c.sub.1](R + [rho]) + [c.sub.2]/R+[rho] = [c'.sub.0] + [c'.sub.1][rho] +...+ [c'.sub.m][[rho].sup.m] (2)

This formulation (Eq. (2)) is known as the shear deformation theory (SDT). By setting m = 1, the first order shear deformation theory is employed for the analysis. Based on this theory, every deformation component can be stated by two variables including the rotation and displacement. For a symmetric cylindrical shell, the radial and axial components of deformation may be considered as follows [16]

[u.sub.z] = u + [rho][[psi].sub.z] = w + [rho][[psi].sub.r] (3)

where [u.sub.z], [w.sub.r] are the axial and radial components of deformation, respectively, u,w,[[psi].sub.z],[[psi].sub.r] are only functions of the axial component of coordinate system (z). By considering Eq. (3), the strain components are

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (4)

The previous papers did not consider the piezoelectric structure in comprehensive condition and by considering the whole piezoelectric coefficients. The present paper improves the previous incompleteness and considers a FGP cylindrical shell in complete conditions. Therefore, stress-strain relations are [15]

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (5)

where [C.sub.ijkl], [e.sub.ijk] are elastic stiffness and piezoelectric coefficients, [E.sub.k] is electric field component. Based on Eq. (5), the electric field has no effect on the shear stress. By having the components of the electric field, Eq. (5) can be completed. Electric field is equal to negative divergence of the electric potential. The electric field vector is in accordance with direction of decreasing of the electric potential.

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (6)

where [phi] is the electric potential function. Due to the symmetric condition of the problem (symmetric loading, boundary conditions and material properties), Eq. (6) can be reduced using [partial derivative]/[partial derivative]r = [partial derivative]/[partial derivative][rho] as follows

[E.sub.r] = -[[partial derivative][phi]/[partial derivative][rho]], [E.sub.[theta]] = 0, [E.sub.z] = -[[partial derivative][phi]/[partial derivative]z] (7)

Based on the results of the previous papers [15, 17, 18] the electric potential function may be supposed as a quadratic function in the radial direction and an unknown function in the longitudinal direction

[phi](z, [rho]) = [[phi].sub.0](z) + [rho][[phi].sub.1](z) + [[rho].sup.2][[phi].sub.2](z) (8)

By substitution of Eq. (8) into Eq. (7), the electric field (Eq. (6)) can be reduced to

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (9)

The electric displacement may be obtained as a linear combination of the strain and electric field as follows [17, 18]

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (10)

where [[eta].sub.ik] are dielectric coefficient. By having the components of the stresses, strains, electric field and electric displacements (Eqs. (4), (5), (9) and (10)), the energy equation per unit volume may be obtained. Total energy includes the mechanical and electrical energy. Mechanical energy is equal to one half of multiplying the stress tensor in the corresponding strain tensor. Electrical energy is equal to one half of multiplying the electric displacement tensor in the corresponding electric field tensor. Therefore, energy per unit volume ([bar.u]) may be obtained as follows

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (11)

The total energy must be evaluated by integration of Eq. (11) on the volume of the cylinder. The volume element of the cylinder is 2[pi](R + z)dzdx. Therefore, the total energy of the system is

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (12)

where [F.sub.1](u, w, [[psi].sub.z], [[psi].sub.r], [[phi].sub.0], [[phi].sub.1], [[phi].sub.2], z) is the appropriate functional of the system which can be obtained as follows

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (13)

F(u, w, [[psi].sub.z], [[psi].sub.r], [[phi].sub.0], [[phi].sub.1], [[phi].sub.2], z) can be decomposed to three types of sentences as follows

(u, w, [[psi].sub.z], [[psi].sub.r], [[phi].sub.0], [[phi].sub.1], [[phi].sub.2], z) = [U.sub.S](z) + [U.sub.Piezo] (z) - [U.sub.Die](z) (14)

These sentences include strain energy US (z), piezoelectric energy [U.sub.Piezo](z) and dielectric energy [U.sub.Die](z).

2.1. Calculation of the external works

External works such as pressure [16] or rotational loads [25] can be considered in this section. Energy of internal pressure is equals to multiplying the pressure in the radial deformation of the inner surface of the cylinder. Inner pressure applies in the same direction of the deformation. Eq. (15) indicates work is done by the internal pressure. Fig. 1 shows the schematic figure of the cylindrical pressure vessel.

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (15)

[FIGURE 1 OMITTED]

2.2. Variation of the energy equation

Total energy of the system is obtained by subtraction of Eq. (15) from Eq. (14) as follows

F(u, w, [[psi].sub.z], [[psi].sub.r], [[phi].sub.0], [[phi].sub.1], [[phi].sub.2], z) = [U.sub.S](z) + [U.sub.Piezo](z) - [U.sub.Die](z) - [W'.sub.1] (16)

Every terms of above functional [U.sub.S](z), [U.sub.Piezo](z), [U.sub.Die](z) are demonstrated as follows

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (17)

Eq. (16) includes seven functions. By using Euler equation, variation of Eq. (16) can be expressed as follows

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (18)

Using the Euler equation, final governing differential equation of the system in matrix form is

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (19)

where matrices [G.sub.i], F are functions of [A.sub.i], [C.sub.i], [D.sub.i] that are demonstrated earlier in Eq. (17). The complete set of partial differential equations for a functionally graded piezoelectric shell of revolution with variable thickness and curvature can be studied in future work of authors [26]. The functions of [A.sub.i], [C.sub.i], [D.sub.i] are demonstrated as follows

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]

Matrix F has two nonzero components as follows

[F.sub.3,1] = 2[P.sub.i](R - [h/2]), [F.sub.4,1] = -2[h/2][P.sub.i](R - [h/2])

3. Solution of the problem

For analysis of the problem and comparing the results of the present method with plane elasticity theory (PET), the solution of the problem must be evaluated at the regions that are adequate far from two ends of the cylinder. This solution is evaluated using Eq. (19) as follows

[[G.sub.3]]{X} = {F}, {X} = [{u [[psi].sub.x] w [[psi].sub.z] [[phi].sub.0] [[phi].sub.1] [fa.sub.2]}.sup.T] (20)

The above formulation is used to evaluate the validity of the present method (first order shear deformation theory).

The numerical value of the physical parameters is selected as follows [15]

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (21)

The cylinder is made of functionally graded material that is graded in the radial direction. Therefore the entire properties must be represented as a power function in terms of the radial coordinate.

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (22)

where [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] are the values of the elastic, piezoelectric and dielectric coefficient, respectively, at inner radius of the cylinder.

4. Results

As mentioned earlier, this paper deals with the analysis of a FGP cylinder at regions that are adequate far from two ends of the cylinder. By setting [n.sub.1] = [n.sub.2] = [n.sub.3] = n, Fig. 2 shows the radial distribution of the radial displacement along the thickness for five values of nonhomogenous index (n = 0, [+ or -] 1, [+ or -] 2).

Fig. 3 shows the radial distribution of the electric potential along the thickness direction for five values of nonhomogenous index under 80 MPa internal pressure.

The calculations indicate that the circumferential and axial stresses are two main components of the stress tensor. The circumferential stress is maximum component of stress tensor. The previous formulation has not ability to evaluate the axial stress [15]. Evaluation of the axial stress indicates that the value of the axial stress is significant and must be considered in design calculations.

[FIGURE 2 OMITTED]

Figs. 4 and 5 show the radial distribution of the circumferentialand axial stresses along the thickness direction for five values of nonhomogenous indexes (n = 0, [+ or -] 1, [+ or -] 2).

Fig. 4 indicates that the maximum hoop stress is located at inner radius and the minimum of that at outer radius. The decreasing of the stress from inner to outer radii is maximum for n = -2 and this value decreases with increasing the values of nonhomogenous index. n = 2 presents a uniform distribution of stress along the thickness. Fig. 5 is similar to Fig. 4, behaviorally.

[FIGURE 3 OMITTED]

Figs. 4, 5 indicate that the inner pressure impresses significantly the values of stress at the inner radius rather than the outer radius that pressure is zero. In the other word, the value of stress at the inner radius extremely depends on the value of non homogenous index, while the value of stress at the outer radius weakly depends on the value of nonhomogenous index.

[FIGURE 4 OMITTED]

[FIGURE 5 OMITTED]

For validation, it is appropriate to compare these results with whose results that is obtained using the finite element method. Fig. 6 shows comparison between the obtained results using three methods (PET, FSDT and FEM).

[FIGURE 6 OMITTED]

The main objective of this paper is verification of the FSDT results for electro elastic analysis of a FGP cylinder. Figs. 7, 8 show the radial distribution of the radial displacement and electric potential along the thickness direction for five values of nonhomogenous indexes (n = 0, [+ or -] 1, [+ or -] 2) based on two theories FSDT and PET. Thick lines represent the value of components based on the plane elasticity theory (PET) and thin lines represent the value of component based on FSDT. The numerical difference between two theories is presented in Tables 1, 2 for radial displacement and electric potential, respectively.

[FIGURE 7 OMITTED]

[FIGURE 8 OMITTED]

Minor difference between two theories arises from the first assumption of the problem. The radial displacement based on the FSDT is considered as a linear function of thickness, while the plane elasticity theory solves the problem using the analytical method and calculation of characteristic equation [15]. In spite of plane elasticity theory, the FSDT has an ability to solve the problem in the two dimensional coordinate systems with appropriate boundary conditions. The plane elasticity theory can present two dimensional responses of a FG cylinder only with simply supported end conditions [14].

Fig. 9 shows the radial distribution of percentage of the difference between radial displacements using two theories. It is observed that the maximum difference is located at the surface of applied pressure. This difference decreases uniformly from the inner radius to the middle of the cylinder, approximately. From the middle surface to the outer surface, the difference increases uniformly.

[FIGURE 9 OMITTED]

5. Discussion and conclusion

Thermoelastic analysis of a FGP cylinder was investigated using the FSDT and energy method in this work. The main results that are concluded from the present paper are classified as follows.

1. The distribution of the radial displacement indicates that the maximum value of the radial displacement is located at the inner radius and the minimum value of the radial displacement is located at the outer radius. This result is accordance with the results of the literature [15].

2. The radial distribution of hoop and axial stresses indicates that the inner pressure impresses significantly the stress at inner radius. This distribution indicates that the value of stress at the inner surface of the cylinder depends strongly on the values of nonhomogenous index. This phenomenon is not repeated at outer surface because of zero outer pressure.

3. The comparison between the PET and FSDT indicates that the present results using FSDT have not significant difference with the results using PET [15]. Especially the radial displacement is strongly in accordance with the results of the plane elasticity theory. This accordance indicates that the first order shear deformation theory has sufficient capability to simulate the displacement with well precision. Therefore the first order shear deformation theory (FSDT) can be employed for the analysis of a functionally graded piezoelectric structure as an excellent theory.

4. Finite element modeling has been performed for simulation of the results that is obtained using the first order shear deformation and plane elasticity theories. The obtained results using FEM justified acceptability of the results using the FSDT and PET.

5. The distribution of the axial stress indicates that this component of stress must be regarded in the design calculation. The value of the axial stress is significant in contrast with the order of the circumferential stress.

References

[1.] Timoshenko, S.P. 1976. Strength of Materials: Part II (Advanced Theory and Problems), 3 rd edition, Van Nostrand Reinhold Co, New York.

[2.] Naghdi, P.M.; Cooper, R.M. 1956. Propagation of elastic waves in cylindrical shells including the effects of transverse shear and rotary inertia, Acoustical Sc. America 28(1): 56-63.

[3.] Mirsky, I.; Hermann, G. 1958. Axially motions of thick cylindrical shells, Appl. Mech. 25: 97-102.

[4.] Yamanouchi, M.; Koizumi, M.; Shiota, I. 1990. Proceedings of the first international symposium on functionally gradient materials, Sendai, Japan.

[5.] Jabbari, M.; Sohrabpour, S.; Eslami, M.R. 2002. Mechanical and thermal stresses in a functionally graded hollow cylinder due to radially symmetric loads, Int. J. Pres. Ves. Pip. 79: 493-497. http://dx.doi.org/10.1016/S0308-0161(02)00043-1.

[6.] Chen, W.Q.; Lu, Y.; Ye, J.R.; Cai, J.B. 2002. 3D electroelastic fields in a functionally graded piezoceramic hollow sphere under mechanical and electric loading, Arch. Appl. Mech. 72: 39-51. http://dx.doi.org/10.1007/s004190100184.

[7.] Liu, X.; Wang, Q.; Quek, S.T. 2002. Analytical solution for free vibration of piezoelectric coupled moderately thick circular plates, Int. J. Solid. Struct. 39: 2129-2151. http://dx.doi.org/10.1016/S0020-7683(02)00081-1.

[8.] Wu, L.; Jiang, Z.; Liu, J. 2005. Thermoelastic stability of functionally graded cylindrical shells, Compos. Struct. 70: 60-68. http://dx.doi.org/10.1016/j.compstruct.2004.08.012.

[9.] Shi, Z.; F. Chen, Y. 2004, Functionally graded piezoelectric cantilever beam under load, Arch. Appl. Mech. 74: 237-247.

[10.] Peng-Fei, H.; Andrew, Y.T. 2004. The transient responses of magneto-electro-elastic hollow cylinders, Smart. Mater, Struct. 13: 762-776. http://dx.doi.org/10.1088/0964-1726/13/4/014.

[11.] Lu, P.; Lee, H.P.; Lu, C. 2005. An exact solution for functionally graded piezoelectric laminated in cylindrical bending, Int. J. Mech. Sci. 47: 437-458. http://dx.doi.org/10.1016/j.ijmecsci.2005.01.012.

[12.] Dai, H.L.; Fu, Y.M.; Yang, J.H. 2007. Electromagnetoelastic behaviors of functionally graded piezoelectric solid cylinder and sphere, Acta. Mech. Solida. Sinica. 23: 55-63. http://dx.doi.org/10.1007/s10409-006-0047-0.

[13.] Babaei, M.H.; Chen, Z.T. 2008. Exact solutions for radially polarized and magnetized magnetoelectroelastic rotating cylinders, Smart. Mater. Struct. 17:025035, 11p.

[14.] Jabbari, M.; Bahtui, A.; Eslami, M.R 2009. Axisymmetric mechanical and thermal stresses in thick short length FGM cylinders. Int. J. Pres. Ves. Pip. 86(5): 296-306. http://dx.doi.org/10.1016/j.ijpvp.2008.12.002.

[15.] Khoshgoftar, M.J.; G. Arani, A.; Arefi, M. 2009. Thermoelastic analysis of a thick walled cylinder made of functionally graded piezoelectric material. Smart. Mater. Struct. 18: 115007. http://dx.doi.org/10.1088/0964-1726/18/11/115007.

[16.] Arefi, M.; Rahimi, G.H. 2010, Thermo elastic analysis of a functionally graded cylinder under internal pressure using first order shear deformation theory, Sci. Res. Essays 5(12): 1442-1454.

[17.] Sheng, G.G., Wang, X. 2010, Response and control of functionally graded laminated piezoelectric shells under thermal shock and moving loadings, Compos. Struct. 93: 132-141. http://dx.doi.org/10.1016/j.compstruct.2010.06.007.

[18.] Sheng, G.G.; Wang, X. 2010, Thermoelastic vibration and buckling analysis of functionally graded piezoelectric cylindrical shells. Appl. Math. Modelling 34: 2630-2643. http://dx.doi.org/10.1016/j.apm.2009.11.024.

[19.] Dai, H.L.; Hong, L.; Fu, Y.M.; Xiao, X., 2010. Analytical solution for electromagnetothermoelastic behaviors of a functionally graded piezoelectric hollow cylinder, Appl. Math. Modelling. 34: 343-357. http://dx.doi.org/10.1016/j.apm.2009.04.008.

[20.] Ghannad, K.M.; Rahimi, G.H.; EsmaeIlzadeh, Kh.S. 2008. General plane elasticity solution of axisymmetric functionally graded thick cylindrical shell, Technical and Engineering Journal of Modares (in Persian).

[21.] Ghannad, K.M., Rahimi, G.H.; EsmaeIlzadeh, Kh.S. 2008. General shear deformation of axisymmetric functionally graded thick cylindrical shell, Technical and Engineering Journal of Modares (in Persian).

[22.] Z. Nejad, M.; Rahimi, G.H.; Ghannad, K.M. 2009. Set of field equations for thick shell of revolution made of functionally graded materials in curvilinear coordinate system, Mechanika 3(77): 18-26.

[23.] Ghannad, M.; Z. Nejad, M.; Rahimi, G.H. 2009. Elastic solution of axisymmetric thick truncated conical shells based on first-order shear deformation theory, Mechanika 5(79): 13-20.

[24.] Arefi, M.; Rahimi, G.H. 2011. Comprehensive thermoelastic analysis of a functionally graded cylinder with different boundary conditions under internal pressure using first order shear deformation theory, Mechanika (accepted for publication).

[25.] Rahimi, G.H.; Arefi, M.; Khoshgoftar, M.J. 2011. Application and analysis of functionally graded piezoelectrical rotating cylinder as mechanical sensor subjected to pressure and thermal loads, Appl. Math. Mech. (Engl. Ed) 32(8): 1-12 (In Press).

[26.] Arefi, M.; Rahimi, G.H. 2011. Three dimensional multi field equations of a functionally graded piezoelectric thick shell with variable thickness, curvature and arbitrary nonhomogeneity, Acta Mech, (In Press).

G.H. Rahimi, M. Arefi, M.J. Khoshgoftar

Department of Mechanical Engineering, Tarbiat Modares University, Tehran, Iran, Email: [email protected]

doi: 10.5755/j01.mech.18.3.1875

Table 1

Comparison between the radial displacements based on two theories

r      [rho]   Theories    n = 0     n = 1    n = -1     n = 2

0.6    -0.2      FSDT     0.00203   0.00157   0.00256   0.00119
                 PET      0.00223   0.00176   0.00276   0.00137

0.64   -0.16     FSDT     0.00198   0.00153   0.00249   0.00116
                 PET      0.00213   0.00168   0.00264   0.00130

0.68   -0.12     FSDT     0.00192   0.00149   0.00242   0.00113
                 PET      0.00204   0.00160   0.00254   0.00124

0.72   -0.08     FSDT     0.00187   0.00145   0.00236   0.00110
                 PET      0.00196   0.00154   0.00245   0.00119

0.76   -0.04     FSDT     0.00182   0.00141   0.00229   0.00107
                 PET      0.00190   0.00149   0.00237   0.00114

0.8      0       FSDT     0.00177   0.00137   0.00222   0.00104
                 PET      0.00184   0.00144   0.00230   0.00110

0.84   0.04      FSDT     0.00172   0.00133   0.00216   0.00101
                 PET      0.00179   0.00140   0.00223   0.00107

0.88   0.08      FSDT     0.00166   0.00130   0.00209   0.00098
                 PET      0.00174   0.00136   0.00218   0.00104

0.92   0.12      FSDT     0.00161   0.00126   0.00203   0.00096

                 PET      0.00170   0.00133   0.00213   0.00102

0.96   0.16      FSDT     0.00156   0.00122   0.00196   0.00093
                 PET      0.00166   0.00130   0.00208   0.00100

1       0.2      FSDT     0.00151   0.00118   0.00189   0.00090
                 PET      0.00163   0.00128   0.00204   0.00098

r      [rho]   Theories   n = -2

0.6    -0.2      FSDT     0.00316
                 PET      0.00336

0.64   -0.16     FSDT     0.00307
                 PET      0.00322

0.68   -0.12     FSDT     0.00299
                 PET      0.00310

0.72   -0.08     FSDT     0.00291
                 PET      0.00300

0.76   -0.04     FSDT     0.00282
                 PET      0.00290

0.8      0       FSDT     0.00274
                 PET      0.00281

0.84   0.04      FSDT     0.00266
                 PET      0.00274

0.88   0.08      FSDT     0.00257
                 PET      0.00267

0.92   0.12      FSDT     0.00249
                 PET      0.00261

0.96   0.16      FSDT     0.00240
                 PET      0.00255

1       0.2      FSDT     0.00232
                 PET      0.00250

Table 2

Comparison between the electric potentials based on two theories

r      [rho]   Theories    n = 0      n = 1      n = -1     n = 2

0.6    -0.2      FSDT        0          0          0          0
                 PET         0          0          0          0

0.64   -0.16     FSDT     -5.7E+04   -4.4E+04   -7.2E+04   -3.3E+04
                 PET      -6.2E+04   -5.0E+04   -7.6E+04   -3.8E+04

0.68   -0.12     FSDT     -1.0E+05   -7.8E+04   -1.3E+05   -5.9E+04
                 PET      -1.0E+05   -8.1E+04   -1.3E+05   -6.1E+04

0.72   -0.08     FSDT     -1.3E+05   -1.0E+05   -1.7E+05   -7.7E+04
                 PET      -1.3E+05   -9.8E+04   -1.6E+05   -7.3E+04

0.76   -0.04     FSDT     -1.5E+05   -1.2E+05   -1.9E+05   -8.8E+04
                 PET      -1.4E+05   -1.0E+05   -1.8E+05   -7.6E+04

0.8      0       FSDT     -1.6E+05   -1.2E+05   -2.0E+05   -9.2E+04
                 PET      -1.4E+05   -1.0E+05   -1.8E+05   -7.2E+04

0.84   0.04      FSDT     -1.5E+05   -1.2E+05   -1.9E+05   -8.8E+04
                 PET      -1.2E+05   -9.1E+04   -1.6E+05   -6.4E+04

0.88   0.08      FSDT     -1.3E+05   -1.0E+05   -1.7E+05   -7.7E+04
                 PET      -1.0E+05   -7.5E+04   -1.4E+05   -5.2E+04

0.92   0.12      FSDT     -1.0E+05   -7.8E+04   -1.3E+05   -5.9E+04
                 PET      -7.5E+04   -5.3E+04   -1.0E+05   -3.7E+04

0.96   0.16      FSDT     -5.7E+04   -4.4E+04   -7.2E+04   -3.3E+04
                 PET      -4.0E+04   -2.8E+04   -5.5E+04   -1.9E+04

1       0.2      FSDT        0          0          0          0
                 PET         0          0          0          0

r      [rho]   Theories    n = -2

0.6    -0.2      FSDT        0
                 PET         0

0.64   -0.16     FSDT     -9.0E+04
                 PET      -8.9E+04

0.68   -0.12     FSDT     -1.6E+05
                 PET      -1.5E+05

0.72   -0.08     FSDT     -2.1E+05
                 PET      -2.0E+05

0.76   -0.04     FSDT     -2.4E+05
                 PET      -2.2E+05

0.8      0       FSDT     -2.5E+05
                 PET      -2.2E+05

0.84   0.04      FSDT     -2.4E+05
                 PET      -2.1E+05

0.88   0.08      FSDT     -2.1E+05
                 PET      -1.8E+05

0.92   0.12      FSDT     -1.6E+05
                 PET      -1.3E+05

0.96   0.16      FSDT     -9.0E+04
                 PET      -7.4E+04

1       0.2      FSDT        0
                 PET         0