Thermo-Elastic Analysis of Clamped-Clamped Thick FGM Cylinders by Using Third-Order Shear Deformation Theory

Using the third-order shear deformation theory (TSDT), an analytical solution for deformations and stresses of axisymmetric clamped-clamped thick cylindrical shells made of functionally graded material (FGM) subjected to internal pressure and thermal loading are presented. The material properties are graded along the radial direction according to power functions of the radial direction. It is assumed that Poisson’s ratio is constant across the cylinder thickness. The differential equations governing were generally derived, making use of TSDT. Following that, the set of non-homogenous linear differential equations for the cylinder with clamped-clamped ends was solved, and the effect of loading and supports on the stresses and displacements was investigated. The problem was also solved, using the finite element method (FEM), and the results of which were compared with those of the analytical method. Furthermore, the effect of increases in the temperature gradient on displacement and stress values has been studied. Finally, in order to investigate the effect of third-order approximations on displacements and stresses, a comparison between the results of firstand third-order shear deformation theory has been made.

Latin American Journal of Solids and Structures 13 (2016) 750-774 ties significantly improves the mechanical strength and fracture toughness of FGMs (Cao et al. 2012).Since thick cylindrical shells such as hollow cylinders are used in various industrial fields as structure materials, thermo-elastic analysis of these structures under mechanical and thermal loads, which may cause undesirable stress and deformation, has become important.In order to reduce these stresses and deformations, hollow cylindrical shells made up of FGMs have been widely used.Most of the existing literature deals with thermo-elastic analysis of FG shells.Reddy and Liu (1985) developed a simple third-order shear deformation shell theory, in which the transverse shear strains are assumed to be parabolically distributed across the shell thickness and which contains the same number of dependent unknowns as in the first-order shear deformation theory.Moreover, this theory requires no shear correction factors.Fukui and Yamanaka (1992) the Navier solution to derive the governing equation for a thick-walled FGM tube under internal pressure and solved the equations obtained numerically by means of the Runge-Kutta method.Eipakchi et al. (2003) investigated the governing equations of homogeneous cylinders with variable thickness, using FSDT and presented the solution of the equations, using perturbation theory.They further extended their previous work by considering homogenous and isotropic conical shells with variable thickness, using FSDT and SSDT (second-order shear deformation theory) and solved the conducted equations by perturbation theory (Eipakchi et al. 2008).Xiang et al. (2006) obtained the exact solution of FGM hollow cylinders in the state of plane strain with exponential function of elasticity modulus along the radius.Shi et al. (2007) analyzed heterogeneous cylindrical shells with power function of elasticity modulus by using multilayer method with homogeneous layers.Thickwalled FGM cylinders in plane strain state with exponentially-varying material properties were solved by Tutuncu (2007) using Frobenius method.Shao and Ma (2008) solved the governing difential equations of functionally graded circular hollow cylinders with exponential variation of material properties under mechanical loads and linearly increasing boundary temperature, using series expansion.Nejad, Rahimi, and Ghannad (2009) developed a 3-D set of field equations of FGM thick shells of revolution in curvilinear coordinate system by tensor calculus.Ghannad and Nejad (2010) presented the general method of derivation and the analysis of internally pressurized thick-walled cylinders with clamped-clamped ends.Eipakchi (2010) calculated stresses and displacements of a thick conical shell with varying thickness under non-uniform internal pressure analytically, using third-order shear deformation theory.Arefi and Rahimi (2010) conducted thermo-elastic analysis of a functionally graded cylinder under internal pressure, using first-order shear deformation theory for plane strain condition.Arani et al. (2011) investigated the effect of material inhomogeneity on electro-thermo-mechanical behaviors of FG piezoelectric rotating cylinders under internal and external pressure with varying properties of power function.Keles and Conker (2011) indicated transient hyperbolic heat condition in thick-walled FGM cylinders and spheres with exponentially-varying properties.Jabbari et al. (2011) indicated mechanical and thermal stresses in an functionally graded porous material (FGPM) hollow cylinder subjected to radially symmetric thermal loads and solved the governing equations by using Fourier series.Ghannad et al. (2012) presented a closed-form analytical solution for clamped-clamped thick cylindrical shells with variable thickness subjected to constant internal pressure based on the first-order shear deformation theory (FSDT).Ghannad et al. (2013) investigated the elastic analysis of pressurized thick cylindrical shells with variable thickness made of functionally graded materials.They used the matched asymptotic method (MAM) of Latin American Journal of Solids and Structures 13 (2016) 750-774 the perturbation theory in order to convert the resultant equations into a system of algebraic equations and two systems of differential equations with constant coefficients.Three dimensional free vibration and transient response of a cylindrical panel made of two directional functionally graded materials based on three dimensional equations of elasticity and subjected to internal impact loading is considered by Zafarmand et al. (2015).Nejad et al. (2015) presented a semi-analytical solution for the purpose of elastic analysis of rotating thick truncated conical shells made of FGM under non-uniform pressure by using the first-order shear deformation theory (FSDT) and multilayer method (MLM).Ghannad and Gharooni (2015) investigated displacements and stresses in pressurized thick FGM cylinders with exponential variation of material properties based on TSDT.
Most of the studies reviewed have investigated the heat transfer and gradient temperature of FGM cylinders and have presented no accurate solution especially the ones which use shear deformation theory.The results of previous analyses of the cylindrical shells show that FSDT is not an appropriate theory for the purpose of stress analysis, suggesting that the accuracy of the stresses (especially radial stresses), unlike displacements which have been calculated directly by SDT, are not acceptable.Furthermore, it could be seen that the thermo-elastic analysis of pressurized thick FGM cylinders based on third-order shear deformation theory has not been studied in the literature.On the other hand, other theories such as classical theories or plane elasticity theories (PET) have not acceptable results because of failing to take into account shear stresses.
In order to improve the approximation of shear deformation theory (especially under thermal loading) and to take into account the effect of shear stresses and strains, the general method of derivation and thermo-elastic analysis of pressurized thick-walled cylindrical shells under the effect of gradient temperature has been presented in this paper by using HSDT.The material of the cylinder is assumed to be isotropic and heterogeneous with radially varying mechanical and physical properties continuously along the thickness with a power function and constant Poisson's ratio.The coefficient matrices of the governing equations have been derived in the unique abbreviated form.The obtained equations are solved under the generalized clamped-clamped conditions.Furthermore, the effects of loading, temperature gradient and inhomogeneity of FG materials on the stresses and displacements have been investigated.The analytical results of shear deformation theory are compared with the numerical results of FEM.Finally, the higher efficiency and accuracy of the thirdorder shear deformation theory, compared with the lower one, have been shown.

PROBLEM FORMULATION
In shear deformation theory (SDT), the straight lines perpendicular to the central axis of the cylinder do not necessarily remain unchanged after loading and deformation, suggesting that the deformations are axisymmetric and change along the longitudinal direction of cylinder.This means that the elements have rotation, and the shear strain is not zero.
The parameter r is the radius of every layer of cylinder which can be replaced in terms of radi- us of mid-plane ( ) R and distance of every layer with respect to mid-plane ( ) z , as follows (Figure 1 The parameters x and z are the length and the thickness variables which vary in the following intervals: where h and L are the thickness and the length of the cylinder. The displacement field is assumed as a polynomial of a variable ( ) z through the thickness.As the number of terms in the polynomial function increases, the approximate solution is improved as well.The high-order shear deformation theory (HSDT) is employed to simulate the deformation of every layer of the cylinder.Based on HSDT, every component of deformation can be stated by variables that include the displacement and rotation.For an axisymmetric cylindrical shell, axial and radial components of displacement field are assumed to be in the following form: where 0 ( ) u x and 0 ( ) w x are the displacement components of the middle surface.The kinematic relations in the cylindrical coordinates system for an axisymmetric cylinder are: Latin American Journal of Solids and Structures 13 (2016) 750-774 Considering the dimensionless radial coordinate   r as the ratio of radial coordinate   Mechanical and physical properties of inhomogenous and isotropic material (FGM) including modulus of elasticity   E , thermal expansion coefficient    and thermal conductivity coefficient

 
K are supposed to be a power function of dimensionless radial coordinate: Here i E , i  and i K are the modulus of elasticity, thermal expansion coefficient and thermal conductivity coefficient at the inner surface   i r of the cylinder and n is the inhomogeneity con- stant of FG material determined empirically.Generally speaking, the Poisson's ratio    for a thickwalled cylindrical pressure vessel of isotropic FGM varies in a small range.Furthermore, its effects on mechanical stresses are insignificant.For simplicity, the Poisson's ratio is assumed to be constant.
The distribution of material properties from Eqs. ( 6) can be re-written as a function of z by substituting r from Eq. (1) into Eq. ( 5): Latin American Journal of Solids and Structures 13 (2016) 750-774 ( ) Figure 2 shows the distribution of dimensionless material properties   i X X with respect to the dimensionless radius in a heterogeneous cylinder for integer values of n which vary in the range of On the basis of the constitutive equations for inhomogeneous and isotropic materials, the stress components based on components of mechanical and thermal strain are as follows: where  and  are the Lame's constants.Considering variable elasticity modulus for the FGM materials, these two parameters are as follows: is the distribution of temperature gradient along the thickness of the cylinder in the one-dimensional steady state heat transfer condition.Solids and Structures 13 (2016) 750-774 The axial forces based on normal components of stress are as follows:

Latin American Journal of
The bending moments based on normal components of stress are as follows: The higher-order bending moments based on normal components of stress are as follows: The shear force based on shear stress is as follows: The torsional moment based on shear stress is as follows: The higher-order torsional moments based on shear stress are as follows: Latin American Journal of Solids and Structures 13 (2016) 750-774

VIRTUAL WORK PRINCIPLE
Based on the principle of virtual work, the variation of strain energy of the elastic body   U is equal to the variation of external work due to pressure   The strain energy is and the external work consists of internal and external pressure ( i P and o P ) , which is: where sf f  is the surface force of the pressurized cylinder.
The bounds of integrals for x and z variables have been defined in Eq. ( 2).
Variation of the strain energy can be expressed as follows: and the variation of the external work is: Solids and Structures 13 (2016) 750-774 Substituting Eqs. ( 4), ( 7) and ( 8) into Eqs.( 22) and ( 23), using Eq. ( 19) and carrying out the integration by parts, the equilibrium equations for the cylindrical shell with constant thickness under one-dimensional heat transfer and uniform internal and external pressure are obtained in the form of: where P F and T F stand for non-homogeneity of the governing equations which have resulted from the loading of pressure and temperature, respectively.The subscripts x and z in the right terms of each equation show the components of P F and T F along the axial and radial direction, respec- tively.
2 0 0,1, 2, 3 1 The boundary conditions at the two ends of the cylinder are: Eqs. ( 24) express the main governing equations based on the HSDT for the cylindrical shells under mechanical and thermal loading.Eq. ( 27) is the boundary conditions which should be satisfied at the two ends of the cylinder.Structures 13 (2016) 750-774 4 HEAT TRANSFER EQUATION According to Figure 1, heat transfer in the form of conduction along the thickness of the cylinder under the effect of heat sources at the internal and external surface causes thermal strains.For thermo-elastic analysis of the cylinder, the function of temperature gradient distribution along the longitudinal direction should be determined by the solution of heat transfer equation.Fourier's heat transfer law in general form is:

Latin American Journal of Solids and
where T ,   K r , q ,  and C  are temperature distribution, thermal conduction coefficient, heat generation rate, density and specific heat, respectively.In cylindrical coordinate system, the Eq. ( 28) could be re-written in the form of: Axisymmetric thermal load in cylindrical coordinate system causes just radial flow in the cylinder.Therefore, Eq. ( 29) for radial temperature distribution in steady-state condition without any heat generation has been simplified as: Thermal boundary conditions, as shown in Figure 1, consist of temperatures at the internal and external surfaces of the cylinder: Solution of differential Eq. ( 30) by applying thermal boundary conditions according to Eq. ( 31) and considering coefficients of thermal conduction from Eqs. (6-c) yield to the temperature distribution along the thickness of the cylinder as a function of radius: Since the terms of temperature gradient distribution appear in the thermal part of nonhomogeneity in Eq. ( 26)   T F , Eq. ( 32) should be re-written in the terms of ( )

ANALYTICAL SOLUTION
Eqs. ( 24) are the set of differential equations.In order to solve these equations, forces and moments could be written by using Eqs.( 11) to ( 18) in terms of stresses.The stresses could be written in terms of strains by substituting Eqs.(7-a and b) into Eq.( 8) and using Eq. ( 33).By using Eqs.
(4), the strains are converted into the displacement filed components.Finally, a set of linear nonhomogenous differential equations with constant coefficients would result, as follows: where   8 8 A  ,   8 8 B  and   8 8 C  are the coefficient matrices and   F is the force vector, which can be expressed as the set of non-homogeneity of differential equations'.  y is the unknown vector including the components of displacement field as:   in the Eq. ( 35) whose reverse would be needed in the next calculations is irreversible.In order to make   , the first equation in the set of Eqs. ( 24) has been integrated.
In Eqs.(24), it is apparent that 0 u does not exist, but 0 du dx does.In order to calculate de- formations in Eqs.(4), 0 du dx is needed.Therefore, by assuming 0 du dx   as a new parameter which could be indicated in the following terms, we have: Applying the mentioned changes, the unknown vector   y in the set of differential Eqs. ( 24) would be rewritten as follows:     The corresponding coefficient matrices   A ,   B and   C of the new differential Eqs. ( 35) have been defined in appendices.
The solution of Eqs. ( 35) consists of general and particular parts: is substituted in homogeneous Eq. ( 26).
Considering that mx e is not equal to zero, the following determinant which is equal to zero would result.
The above determinant is a sixteen-order polynomial which is a function of m .The determinant's roots are the eigenvalues i m consisting of eight pairs of conjugated roots where a pair of the roots is zero.Substituting the calculated eigenvalues in Eq. ( 42), the corresponding eigenvectors   i V are obtained.Therefore, the general solution has been obtained.
Given that   F in Eq. ( 35) consists of constant parameters, the non-homogenous part of the so- lution for axisymmetric cylinder with constant thickness under uniform pressure is not the function of x .Therefore, the particular solution can be expressed as follows.
Considering Eq. ( 27), clamped-clamped boundary conditions at two ends of the cylinder are as follows: Latin American Journal of Solids and Structures 13 (2016) 750-774 Applying eight boundary conditions at each end of the cylinder, one can calculate sixteen constants comprised of 1 1 4 , , C C  in the general solution and 0 15 , C C in the particular solution.Final- ly, the unknown vector   y which consists of displacement field components would be obtained in terms of x variable based on Eq. ( 41) by determining unknown constants.Using Eq. (3) would yield radial and axial displacements.Stress distribution would be obtained by using Eqs.( 4) and (8).

NUMERICAL RESULTS AND DISCUSSIONS
As a case study, a thick heterogeneous cylinder under uniform internal pressure and heat conduction resulting from temperature gradient of internal and external surfaces and clamped-clamped boundary conditions at the two ends with the following characteristics could be considered: In order to demonstrate the potentials of the presented analytical solution for the purpose of analyzing an FG cylinder, a numerical solution is investigated.The ANSYS 14.5 package is used in the static analysis of thick hollow cylinder with constant thickness.The PLANE82 element in the axisymmetric mode, which is an element with eight nodes and two translational degrees of freedom in the axial and radial directions per each node, has been used to model the mechanical part of the analysis.For the thermal part of the analysis, the appropriate element is PLANE77, which is compatible with initial element for superposing the thermal and mechanical results.The length and the thickness of the cylinder have been divided into 1600 and 40 parts, respectively.The quadratic mapped meshing has been applied for each section of the cylinder.In order to consider power-law varying elastic modulus along the thickness of the cylindrical shell, the thickness of the cylinder has been divided into some homogeneous layers.As applying more than 40 layers along the thickness shows no significant effect in FEM results, the thickness of the cylinder is divided into 40 homogeneous layers.Each layer's properties have been defined as a power function of the distance of layer's middle from the internal layer.Finally, the cylindrical shell consists of some coherent homogeneous layers.In order to merge the layers, the properties of contact location of the layers have been defined the average of left and right limits of the boundaries of the two adjacent layers.For thermal boundary conditions, the temperature gradient values of the inner and outer surfaces with the environment temperature have been applied to the corresponding nodes.Clamped boundary conditions Latin American Journal of Solids and Structures 13 (2016) 750-774 have been exerted by preventing the nodes around the two ends of the cylinder from movement.In the next sections, the numerical and analytical results have been investigated.

Loading and Inhomogeneity Effect
The distribution of the dimensionless radial displacement resulting from the numerical and analytical solution in the middle of a cylinder under o 100 C temperature gradient is depicted in Figures 3  and 4 without and with internal pressure, respectively.It is evident that radial displacement under temperature gradient for negative values of n is higher than the homogenous materials at layers close to the internal surface while at the outer surface, reverse holds true.For positive values of n, the reverse holds true, which means that the heterogonous materials have lower values of radial displacement than the homogenous ones at inner surfaces and the higher values of displacement at outer surfaces.The variation in the displacement of heterogeneous materials is similar to that of homogenous materials close to the middle layer.Under thermal and mechanical loads for negative values of n, the displacements of FGM cylinders are higher than those of the homogeneous cylinder.For positive values of n, the situation is reverse, i.e. the displacement is lower.Figure 4 shows that radial displacements for different values of n increase from internal layer to the middle layer while in the external half of the cylinder's thickness, they decrease around the middle layer to the external layer.Variation of inhomogeneity constants from negative to the positive causes maximum displacement of cylinder's thickness to appear around the external layer.Furthermore, adding internal pressure to the temperature gradient causes uniform displacement along the radial direction.From the viewpoint of low level and small variation of displacements, using materials with negative inhomogeneity constants for thermal loading and positive for mechanical loading is recommended.Figures 7 and 8 show the distribution of the dimensionless circumferential stress resulting from the numerical and analytical solution at 2 x L  for o 100 C temperature gradient under pressurized and non-pressurized cylinder, respectively.For both loadings, the circumferential stress for negative values of n is higher than the homogenous materials at layers close to the internal surface while at the outer surface, reverse holds true.For positive values of n, the reverse holds true, suggesting that the heterogonous materials have less values of circumferential stress than the homogenous ones at inner surfaces and the higher values of stress at outer surfaces.Considering more uniform stress distribution of the layers and less maximum values of stress for 0 n  under combined loading, it is interesting to use FG materials with negative values of n.Internal pressure causes an increase in the values of circumferential stresses, which would result in positive values of stresses around the internal layer.Figure 9 shows the distribution of shear stress resulting from analytical solution along the longitude of cylinder under combined load for 1 n   in different layers.It can be seen that there are shear stresses near the two ends of the cylinder under the effect of clamped boundary conditions, which has been considered in shear deformation theory.Dimensionless radial displacement and von Mises stress distribution under combined load at .

Temperature Gradient Effect
The effect of temperature gradient on radial displacement for 1 n   has been investigated in Figure 12.It could be observed that for 1 n   , the same graph is obtained.It can be seen that increasing temperature gradient leads to increment in values of radial displacements.This increment is equal from the internal layer to the external one for lower gradient temperature while it is higher in the external layers for upper gradient temperature for inhomogenous materials.
The effect of temperature gradient on von Mises stress has been investigated in Figures 13 and  14 for 1 n   .The von Mises stresses, as the radial displacements, for both positive and negative values of n are uniform through the thickness in lower temperature gradient.The increment of temperature gradient yields higher stress values.Upper temperature gradient has brought about higher values of stresses around the internal and external layers and low levels of stresses around the middle layer for 1 n   while it causes reductions in stresses from the internal to the external layer for 1 n   .Therefore, the increment of temperature gradient yields higher displacement and stress along the thickness of the cylinder.

Third-order Approximation Effect
Tables 1 and 2 3 and 4, respectively.The stresses resulting from internal pressure increase from the external layer to the internal one FG materials with positive inhomogeneities yield lower values of the maximum stresses.It has been shown that FSDT method has acceptable results for displacements, but the radial and von Mises stresses resulting from FSDT show a significant difference once compared with the results calculated from HSDT and FEM solution.This difference increases at layers close to the boundaries and the greatest difference occurs in the internal surface   2 z h   while using HSDT yields the values of -80 and 0 for radial stresses at the internal and external layer, respectively, which is compatible with applying pressure.In the cylinder under temperature gradient, the radial stresses at different layers have negative values which are inconsiderable relative to the stresses of pressurized cylinder.
,MPa It is generally observed that modifying initial approximation of displacement field components in comparison with FSDT, HSDT yields error reduction.The reason is that FSDT assumes linear distribution for the radial displacement while according to the Lame's theory (PET), the variation of radial displacement along the thickness of the cylinder has hyperbolic distribution.ly be seen that around the middle layer, lower approximation has acceptable results.However, around the internal and external layers, FSDT yields different distributions resulting from Eq. ( 47).This phenomenon intensifies for thermal loading because displacement distribution along the thickness changes from linear status to nonlinear one.The polynomial distribution of displacement in third-order approximation improves the accuracy of the results.The calculus of stresses from displacements and strains by using constitutive relations intensifies the initial error appearing in displacements.

CONCLUSIONS
In this study, the heterogeneous hollow cylinders with power-law varying elastic modulus along radial direction have been solved by HSDT and FEM, and have been compared with homogenous cylinders.At the boundary areas of a thick-walled cylinder with clamped-clamped ends under thermal and mechanical load, use cannot be made of PET, and SDT must be used instead.The shear stress in boundary areas cannot be ignored, but in areas further away from the boundaries, it could be.Therefore, the non-uniform displacements and stresses at points close to the boundaries are different from the other areas under the effect of shear stresses resulting from boundary conditions and, thus, the PET can be used, provided that the shear strain is zero.It was observed that using negative inhomogeneities in the cylinder causes small decreases under thermal load and considerable increase under mechanical load in radial displacements while the positive inhomogeneities are appropriate under combined loading since less radial displacements appear in the cylinder.In general, the effect of inhomogeneity constants on displacement distribution along the thickness for different temperature gradients is inconsiderable compared with internal pressure.It is observed that the stresses are the function of longitudinal and radial direction.Because of low level and small variation of circumferential stresses, using materials with negative inhomogeneities is more suitable.Unlike the displacements, von Mises stresses are strongly affected by the temperature gradient variations and inhomogeneity constants.Using materials with negative n for the cylinder under thermal or combined loading causes less von Mises stresses.Furthermore, temperature gradient increment Latin American Journal of Solids and Structures 13 (2016) 750-774 causes the external layer to be critical from the viewpoint of displacement and internal layer from the viewpoint of stress.It can be concluded that using FSDT for thermoelastic analysis of cylindrical shells, because of disregarding third approximation terms in displacement field components, is not appropriate.Therefore, third-order approximation yields improved accuracy, especially around the boundary layers.Because the terms resulting from pressure and temperature gradient have been revealed in the non-homogeneity part of the set of governing differential equations, the superposition principle could be utilized for the effect of combined loading on the basis of linear elasticity.In fact, the non-homogeneity vector   F appears as a force vector in the governing equations.

Appendices
where the parameters  and k are as follows

 


Figure 1 :
Figure 1: Geometry of the thick cylindrical shell under internal pressure and temperature gradient.
functions of x which are used to determine the displace- ment field.

Figure 2 :
Figure 2: Distribution of dimensionless material properties in an FGM cylinder.
gradients between the environment of the cylinder and internal and external surfaces, respectively.Considering ref T as the temperature of environment (reference), we have: of the differential Eqs.(24) would be derived as follows:

.
The values of elasticity modulus, thermal expansion coefficient and thermal conductivity coefficient at the internal radius are has been carried out by writing the program in MAPLE 17.

Figure 3 :Figure 4 :
Figure 3: Dimensionless radial displacement distribution in the middle of the cylinder under temperature gradient   o 100 C T  .

Figure 5 :Figure 6 :
Figure 5: Dimensionless circumferential stress distribution in the middle of the cylinder under temperature gradient   o 100 C T  .

Figure 7 :Figure 8 Figure 8 :
Figure 7: Dimensionless von Mises stress distribution in the middle of the cylinder under temperature gradient   o 100 C T  .

Figure 9 :Figure 11 :
Figures 10 and 11, respectively.It could be observed that displacement and stress along the longitudinal direction are uniform except the layers around the clamped boundaries.Therefore, PET results are valid through the length of the cylinder far away from clamped boundaries.

Figure 12 :
Figure 12: Dimensionless radial displacement distribution in the middle of the cylinder under different temperature gradient for 1 n   .

Figure 13 :Figure 14 :
Figure 13: Dimensionless von Mises stress distribution in the middle of the cylinder under different temperature gradient for 1 n   .
present the radial and von Mises stresses of different layers resulting from different solutions in the middle of heterogeneous cylinder ( 2 x L  ) under internal pressure, respectively.The radial and von Mises stresses of different layers resulting from different solutions in the middle of the heterogeneous cylinder (

Figure 15
Figure 15 illustrates the effect of third and first-order shear deformation theories in the middle of the cylinder under temperature gradient T  and internal pressure   i P for 1 n   .It could easi-

Figure 15 :
Figure 15: Dimensionless radial displacement distribution in the middle of the cylinder based on HSDT, FSDT and FEM for 1 n   .

Table 1 :
Radial stress of different layers under internal pressure at

Table 3 :
Radial stress of different layers under

Table 4 :
von Mises stress of different layers under