Document Type: Research Paper
Authors
^{1} Mechanical Engineering Department, Faculty of Engineering, University of Kurdistan, Sanandaj, Iran
^{2} Mechanical Engineering Department, Khajeh Nasir Toosi University of Technology, Tehran, Iran
^{3} Mechanical Engineering Department, Urmia University, Urmia, Iran
Abstract
Keywords

Mechanics of Advanced Composite Structures 6 (2019) 81 – 94


Semnan University 
Mechanics of Advanced Composite Structures journal homepage: http://MACS.journals.semnan.ac.ir 
Transient Thermal Stresses Analysis in a FPGM Cylinder
N. Habibi ^{a}^{*}, S. Samawati ^{b}, O. Ahmadi ^{c}
^{a }Mechanical Engineering Department, University of Kurdistan, Sanandaj, 6617715175, Iran,
^{b }Mechanical Engineering Department, Khajeh Nasir Toosi University of Technology, Tehran, Iran,
^{c }Mechanical Engineering Department, Urmia University, Urmia, Iran
Paper INFO 

ABSTRACT 
Paper history: Received 20180529 Received in revised form 20180722 Accepted 20190425 
The present study aims to investigate the analysis of stress, strain, displacement, and electrical potential of a thick hollow cylinder made of FGPM under mechanical and thermal loads. Distribution of mechanical property of material is considered along the shell stick through the power distribution function. Thermal loads have been taken to signify the difference of temperature between outer and inner surfaces for each type of mechanical property. After extracting and solving the differential equations in transient state and the observation of mechanical and thermal boundary conditions, governing functions are obtained through the following parameters: thermal conduction nonhomogeneous parameters, thermal linear distribution coefficient, elastic stiffness constant, piezoelectric coefficient, and dielectric constants.




Keywords: Piezoelectric FGPM Hollow cylinder Transient Thermal Stress 


© 2019 Published by Semnan University Press. All rights reserved. 
The cylinders are widely used under applied loads and in recent years, they have experienced the emergence of governing equations on cylinders and methods sought to solve such equations. Jabbari et al. [1] investigated the effect of transient state thermal stress on a FGPM cylinder. They used direct method to solve the problem and compared the results with the results of the power function method.
The presented method is very extensible even to solve the mixed boundary condition problems. They studied the distribution of stress in cylinder for different boundary conditions. Jabbari et al. [2] in another research investigated the effect of mechanical and transient state thermal stresses in FGM cylinder analytically through the Bessel functions. The analysis shows that as the power law indices increase, the mechanical stresses increased as well until they reach a constant value, where the mechanical stresses remain constant from the metal surface to the ceramic surface of the cylinder; this value of mechanical effective stress yields the optimum pressure vessel. Soufiye’s study [3] scrutinized FGM cone cylindrical shell under widespread pressed loads and press of hydrostatic, where some formulas were obtained. GingHua Zhang et al. [4] surveyed the thermodynamical behavior and the effects of type and load on one FGM, applying Van Karmen Theory (VKT) through numerical solution method under the change of the surface shape. HuiShen et al. [5] investigated the postbuckling FGM hybrid shells under pressed loads in thermal hoops; then, they analyzed the postbuckling for an FGM piezoelectric cylindrical shell under hydrostatic pressure and electrical load in thermal hoops. Shao [6] presented threedimensional solution to obtain the stress fields in an FGM cylindrical plate with a finite length; mechanical and thermal loads were tested on the plate and the results were presented in graphs.
Piezoelectric materials show coupled effects between elastic and electric fields; these materials have been widely used as the actuators or sensors in smart composite material systems. Many analytical studies concerning piezoelastic or piezothermoelastic problems have been reported and the results have been published in book forms [7]. Recently, a new type of piezoelectric material named functionally graded piezoelectric materials (FGPM), with material constants varying continuously in terms of thickness direction has been developed [810]. Electromagnetothermoelastic treatment for a FGPM hollow cylinder has been measured [11]. Obata et al. [12] presented the solution for thermal stresses of a hollow thick cylinder made of FGM under 2D transient temperature distribution. Shariyat et al. [13] presented the nonlinear transient thermal stress and elastic wave propagation analyses of thick temperaturedependent FGM cylinders, through a secondorder pointcollocation method. Lü Chen et al. [14] investigated the elastic mechanical behavior of Nanoscaled FGM films incorporating surface energies. Afsar et al. [15] studied the inverse problems of material distributions for prescribed apparent fracture toughness in FGM coatings around a circular hole in infinite elastic media. Farid et al. [16] investigated the 3D temperature dependent free vibration analysis of FGM curved panels resting on twoparameter elastic foundation using a hybrid semianalytic, differential quadrature method. Bagri and Eslami [17] showed the general coupled thermoelasticity of FGM annular disk, considering the Lord–Shulman Theory. Samsam et al. [18] studied the buckling of thick FG plates under mechanical and thermal loads. Jabbari et al. [19] studied an axisymmetric mechanical and thermal stresses in thick, short length FGM cylinder. They applied separation of variables and Complex Fourier series to solve the heat conduction and Navier equation. Thieme et al. [20] presented the Titanium powder sintering for the preparation of a porous FGM destined as a skeletal replacement implant. Jabbari et al. [21] used the generalized Bessel function to investigate the axisymmetric mechanical and thermal stresses in thick short length FGM cylinders. Asghari and Ghafoori [22] investigated the 3D elasticity solution for functionally graded rotating disks. Khoshgoftar et al. [23] presented the thermoelastic analysis of a thick walled FGPM cylinder.
Jabbari et al. [24] investigated the transient solution of asymmetric mechanical and thermal stresses for hollow cylinders made of functionally graded material. They analytically obtained the temperature distribution, as function of radial and circumferential directions and time, using the the method of separation of variables and generalized Bessel function, and a direct method were used to solve the Navier equations, using the Euler equation and complex Fourier series. The difference between reference [24] and the present work is in cylindermaterial, which at present work is used piezoelectric material in FGM cylinder.
In another research [25], an analytical method of a thermoelastic problem for a medium with FG material properties was developed in a theoretical manner for the ellipticcylindrical coordinate system subjected to the presumption that the material properties except for Poisson's ratio and density are assumed to vary arbitrarily with the exponential law in the radial direction. The Stress Functions Extraction in a hollow cylinder under heating and cooling was investigated by Lamba et al. [26].
Thermal elastic stress distribution occurred on FGM long hollow cylinders was analytically defined under thermal, mechanical and thermomechanical loads. In the closed form, solutions for elastic stresses and displacements are obtained analytically using the infinitesimal deformation theory of elasticity [27].
Khobragade et al. [2836] investigated the temperature distribution, displacement function, and stresses of a thin as well as thick hollow cylinder and also have established the displacement function, temperature distribution and stresses of a semiinfinite cylinder. In addition, Cursun et al. [37], has performed an elastic stress analysis of annular FGMs discs under both uniform pressures on the inner surface and a linearly decreasing temperature distribution.
The study of temperature distribution and thermal stresses of a FG thick hollow cylinder with temperaturedependent material properties by Manthena and Kedar [38] was investigated. All the material properties except Poisson’s ratio are assumed to be dependent on temperature and spatial coordinate z. The twodimensional transient heat conduction equation was solved under convective heat transfer condition with varying point heat source.
The application of this structure is mostly in military structures, aerospace, and also the medical engineering. The main difference between this article and other similar studies is in the solution method and how to achieve the answer. In addition, in similar previous studies, more than FGM was used while FPGM was used in this work.
2. Equilibrium Equations
The main objective of the present study is surveying a hollow thick FGPM cylinder with inner and outer radii, a and b, which are shown in Fig. 1. The condition is in transient state; therefore, the component of time is observed in temperature equations.
This picture is the designed FGM Cylinder model in ABAQUS and aims to show a general schema. Strain–displacement equations are as follow:
(1) 
where u and v are the displacement components along r and θ directions, respectively. The stressstrain equations could be expressed as follow [24]:
(21) 

(22) 

(23) 
The material used for this sample is orthotropic.
Fig. 1.Schematic view of a FGPM cylinder
where:
(3) 
In the above relationship, coefficients C_{ij, }e_{ij}, ε_{ii, }g_{ij} (i=1,2), are the constants of elastic, piezoelectric, dielectric and temperature coefficient.
Furthermore, C_{ij}, e_{ij} (i, j = 1, 2, 4) are elastic, piezoelectric constants. Also, Z_{r}, Z_{q}; the material parameters, (E_{r}, E_{θ}); the electric potential, T(r,θ); temperature distribution and α_{r}, α_{θ} are thermal expansion coefficients, respectively. And governing equations on potential electrical field are [24]:
(41) 

(42) 

(43) 
Here, C_{ij}, e_{ij} and χ_{ii} are elastic, piezoelectric, and dielectric constants, g_{ij} (i=1, 2) signify the thermal coefficients. In addition, electrostatic and equilibrium equations are as follow:
(51) 

(52) 

(53) 
Using weightpower equations, the modelling of nonhomogeneous property can be defined as follow [24]:
(6) 

where are the material parameters and m is the powerlow indices of the material. The equilibrium equations are obtained as follow in terms of the displacement components:


(7) 

(8) 

(9) 

where:
(10) 
In the next section, first temperature distribution is obtained in transient state for cylinder, the application of which will solve the differential equations.
2.1. Temperature Distribution
It is assumed that in hollow cylinder, and represent the temperatures in internal and external radii, respectively. Therefore, heat conduction equation for cylindrical coordinates will be as follows:
(11) 
where, R is the inner thermal source, k is the coefficient of thermal conduction, is density and is the specific heat capacity, the components depending on the property of cylindrical source are as follow:
(12) 
where the material parameters and m are is the powerlow indices of the material. Substituting Eq. (12) into Eq. (11), heat conduction equation in FGM cylinder can yield the following equation:
(13) 
Initial and boundary conditions are shown as:
(14) 
The solution of heat conduction equation (11) is assumed to be:
(15) 
is defined in such a way that the boundary conditions for are zero Such that the definition of guarantees the vanishing of boundary conditions for :
(16) 
A and B are the unknown functions obtained from Eqs. [37] and [38]. Therefore:
(17) 

(18) 

(19) 
To solve the Eq. (17), using the method of separation of variables, onecan write:
(20) 
In addition, i is a complex number. where:
(21) 
and are the Bessel function of the first kind:
(221) 

(222) 
On the other hand, substituting equation (20) in relation (17), equation (23) can be obtained:
(23) 
. (Dot) is time derivative, is displacement derivative. By using (20):
(24) 

(25) 

Substituting (20) in previous relations (15, 17&19) gives the following system of differential equations:
(261) 

(262) 

By putting the obtained determinant of coefficients in this differential equation, equal to zero the positive roots will emerge as follow:
(27) 
Substituting Eq (24) into Eq (17) gives:
(28) 
where the following relation is extracted:
(29) 
The previous equation is an ordinary differential equation which will be solved through the application of:
(30) 
The solution of the differential equation (31) is as follows:
dt) 
(31) 
where:
(32) 
therefore:
(33) 
Using the initial condition, indeed, equation (17) gives:
(34) 
So, the equilibrium equation is extracted as follows:
(35) 
where:
(36) 
and as unknown functions are obtained by solving the following 2equations and 2unknown system:
(37) 

(38) 
2.2. Equilibrium Equations Solution
Equations (7), (8), and (9) could be solved using the power series as follow [24]:
(39) 
By substituting (36) and (31) into (7), (8) and (9):

(40) 

(41) 

2k+1) 
(42) 

Solution of homogeneous part of three upper equations is as follow:

(43) 
where , and are constants obtained using boundary conditions and g means the general solution, through substituting (43) into (42), (41) and (40):

(44) 

(45) 

(46) 
In order to provide the evident solution, the above equations must be equal to zero. Therefore, Eigen vectors are as follow:

(47) 
where, is the more general case for η,Also, η as a parameter in equation (43), only one case n is used, Whereas, it is observed in equation (39) that displacements are a set of special cases for to ; therefore, for determining the general solution, is used. Consequently, the general solution is the linear combination of eigen vectors as follows:

(48) 
where, by applying (47), and are obtained. Particular solution is assumed as follows [24]:
(491) 

(492) 

(493) 
By substituting upper equations and their derivations in (40), (41) and (42) and definition of Bessel function of the first kind and the value of temperature distribution, four 3equations and 3unknowns systems would be resulted .The coefficients are obtained by solving four systems as follow:
(501) 

(502) 

(503) 

(504) 
In the last two 2equations and 2unknown systems, since coefficients determinant matrix does not vanish, the possible solution of these systems is the evident solution. Solutions of first 2 systems yield other coefficients. Complete solution of displacement is as follows:
(511) 

(512) 

(513) 
For n = 0, the system of Navier equations would lead to the following single differential equation:

(52) 

2k+1) 
(53) 

The upper differential equations are Euler differential equations, where it could be assumed that the solution of the homogeneous part is as follows:

(54) 
Where , and are obtained using boundary conditions. Substituting the upper equation in homogeneous part of Eqs (52) and (53) would give the following result:
(551) 

(552) 
In order to obtain the nontrivial solution of the above equation, the determinant of coefficients of constants B, C and D must vanish, which leads to the evaluation of the eigenvector η, as presented in the following equation:

(56) 
Thus, the general solution, utilizing the linearity lemma, is a linear combination of all values of eigenvalues, which is obtained from:

(57) 
The particular solution of Eqs. (52) and (53) is assumed as:
(581) 

(582) 
By substituting the upper equations in (52) and (53), four 2equations and 2unknown systems would be found as follow, where by solving it, D will be derived:
(591) 

(592) 

(593) 

(594) 
In the Eq. (59), regarding the 2equations and 2unknown systems, since coefficient determinant matrix does not vanish, possible solution of these systems is the evident solution, and the solutions of first twosystems will yield other coefficients. Thus, the complete solution, using the Eqs. (52), (53) and (58) is:

(601) 
(602) 

(603) 
The integration of the straindisplacement relationship given in Eq. (1) strains, and the results of Eq. (2), will yield the electrical field, radial, shear and hoop stresses, and the electrical displacements.
3. Results and Discussions
This section discusses the results of the analysis of FGPM transient state thermal of cylinders and compares the results with the stated references.In this part, the results obtained from this research are compared with the results of quoted references of [24] and [2]. The property of FGM and piezo are represented in Table 1. Inner radius is a=70 cm and outer radius is b=100 cm. Inner thermal source would be
.
Fig. 2 represents the thermal distribution related to radius and Fig. 3 shows the radial displacement in proportion to radius at θ = π/3; then, the obtained results were compared with the findings of [24] and [2] references. The boundary conditions are shown in Table 2.
The obtained results show good compatibility with the results of corresponding references in [24] and [2]. Maximum of difference between present work and the references [2] and [24] are shown in Table 3:
Table 1. The property of materials used in the analysis of transient state thermal of cylinders
FGM 
piezoelectric 
property 

Metal:
Ceramic:

Refs [24] , [2] 

Table 2. The boundary conditions used in the analysis of transient state thermal of cylinders 

electrical 
boundary conditions mechanical 
Thermal 




Table 3. The comparison of temperature and radial displacement distributions along radius at θ = π/3
Parameter 
Present work 
Reference 
Percentage of error 
Temperature 
50 
46 [24] 
8 
Radial displacement 
1.88 
1.85 [2] 
1 
Fig. 2. The comparison of temperature distribution along radius at θ = π/3.
Fig. 3. The comparison of distribution of radial displacement along radius at θ = π/3.
Fig. 4. Oscillation of the temperature with respect to the time at various radii
The following section, investigates the results of parametric and survey studies for the cylinders. Fig. 5 shows the radial displacement along radius at , for different values of (indices in powerlaw). As the Fig. clearly shows, the increase of the radius would result in the reduction of radial displacement. Figs. 4, 6 and 710 show the oscillation of the temperature, radial displacement, hoop displacement, electrical potential, radial stress and shear stress with respect to the component of time at variation radii, respectively.
By increasing the radius, temperature increases and the domain of vibration gets larger. The value of radial displacement increases and the domain of vibration gets larger; also, the value of hoop displacement and electrical potential increase and the domain of vibration gets smaller. The value of radial stress and shear stress decrease and the domain of vibration gets smaller. Note that in the presented diagrams, the variations of parameters are dependent on the variation of ; i.e , when increases temperature, radial displacement, electrical potential (except in outer surface where electrical potential equals the value of boundary conditions, zero), radial stress, and hoop stress decrease. According to the formerly stated points, when gets larger, the resistance of FGM cylinder against thermal conduction increases. Also, the increase of leads to the decrease of stress in source; thus, m is of paramount importance in the selection of the type of the material. Variation of is effective on the variation of confidence coefficient.
Fig. 5. The distribution of radial displacement along radius at θ = π/3.
Fig. 6. Oscillation of the radial displacement with respect to the time at various radii.
Fig. 7. Oscillation of the hoop displacement with respect to the time at various radii.
Fig. 8. Oscillation of the electrical potential with respect to the time at various radii.
Fig. 9. Oscillation of the radial stress with respect to the time at various radii.
Fig. 10. Oscillation of the shear stress with respect to the time at various radii.
4. Conclusion
The results of the present study can be summarized as follow:
References
[1] Jabbari M, Sohrabpour S, Eslami MR, General Solution for mechanical and thermal stresses in a functionally graded hollow cylinder due to nonaxisymmetric steadystate loads, Journal of Applied Mechanics 2003, 70;111118.
[2] Jabbari M, Mohazzab, AH, OneDimensional Moving Heat Source in a Hollow FGM Cylinder, Journal of Pressure Vessel Technology 2009; 131.
[3] Soufiye, AH, buckling analysis of FGM circular shells under combined loads and resting on the Pasternak type elastic foundation, 2011.
[4]ShiRong Li, JingHua Zhang YongGang, Zhao, Nonlinear Thermomechanical PostBuckling of Circular FGM Plate with Geometric Imperfection, 2008.
[5] HuiShen, Shen Noda N, Postbuckling of pressureloaded FGM hybridcylindrical shells in thermal environments, Composite Structures 2007; 77: 546 560.
[6] Shao ZS, Wang TJ, Threedimensional solutions for the stress fields in functionally graded cylindrical panel with finite length and subjected to thermal/mechanical loads. Int J Solids Struct 2006; 43: 3856–3874.
[7] Ding HJ, Chen, WQ, three dimensional problems of piezoelasticity, Nova Science Publishers, New York, 2001.
[8] Wu, CCM, Kahn M, Moy W, Piezoelectric ceramics with functionally gradients: a new application in material design. J Am Ceram Soc 1996; 79: 809–812.
[9] Shelley WF, Wan S, Bowman KJ, functionally graded piezoelectric ceramics. Mater Sci Forum 1999; 308(311): 515520.
[10] Zhu, XH, Zu J, Meng ZY, Zhu JM, Zhou SH Li Q, Micro displacement characteristics and microstructures of functionally graded piezoelectric ceramic actuator, Mater Des 2000; 21: 561–566.
[11] Dai HL, Hong L, Fu, YM, Xiao X, Analytical solution for electro magneto thermoelastic behaviors of a functionally graded piezoelectric hollow cylinder, Appl Math Model 2010; 34: 343–357.
[12] Obata Y, Kanayama K, Ohji T, Noda N, TwoDimensional Unsteady Thermal Stresses in a Partially Heated Circular Cylinder Made of Functionally Graded Material, Proceedings of Thermal Stresses, Pub. Branti Zew, Krakow, Poland 1999; 595–598.
[13] Shariyat M, Lavasani SMH, Khaghani M, Transient thermal stress and elastic wave propagation analyses of thick temperaturedependent FGM cylinders, using a secondorder pointcollocation method, Applied Mathematical Modelling 2009,
[14] Lü CF, Chen WQ, Lim CW, Elastic mechanical behavior of Nanoscaled FGM films incorporating surface energies, Composites Science and Technology 2009;69:11241130.
[15] Afsar AM, Sekine H, Inverse problems of material distributions for prescribed apparent fracture toughness in FGM coatings around a circular hole in infinite elastic media, Composites Science and Technology 2002; 62:10631077.
[16] Farid M, Zahedinejad P, Malekzadeh P, Threedimensional temperature dependent free vibration analysis of functionally graded material curved panels resting on twoparameter elastic foundation using a hybrid semianalytic, differential quadrature method, Materials and Design 2010;31: 213.
[17] Bagri A, Eslami MR, generalized coupled thermoselasticity of functionally graded annular disk considering the Lord–Shulman theory, Composite Structures 2008;83: 168179.
[18] Samsam Shariat, BA, Eslami, MR, buckling of thick functionally graded plates under mechanical and thermal loads, Composite Structures 2007; 78: 433–439.
[19] Jabbari M, Bahtui A, Eslami MR, Axisymmetric mechanical and thermal stresses in thick short length functionally graded material cylinder, International Journal of Pressure Vessels and Piping 2009; 86: 96306.
[20] Thieme M, Wieters KP, Bergner F, Scharnweber D, Worch H, Ndop J, Kim TJ, Grill, Titanium powder sintering for preparation of a porous FGM Destined as a skeletal replacement implant, Materials Science Forum 1999; 308(311): 374382.
[21] Jabbari M, Bahtui A, Eslami MR, Axisymmetric mechanical and thermal stresses in thick short length FGM cylinders, Int. J. Press Vessel Pip 2009; 86: 296–306.
[22] Asghari M, Ghafoori E, A threedimensional elasticity solution for functionally graded rotating disks, Compos. Struct 2010; 92: 10921099.
[23] Khoshgoftar MJ, Ghorbanpour Aranl A, Arefi A, Thermoelastic analysis of a thick walled cylinder made of functionally graded piezoelectric material, Smart Mater. Struct 2009; 18: 115007.
[24] Jabbari M, Vaghari AR, Exact solution for asymmetric transient thermal and mechanical stresses in FGM hollow cylinders with heat source, Structural Engineering and Mechanics 2008; 29 (5): 551565.
[25] Manthena VR, Lamba NK, Kedar GD, Thermal stress analysis in a functionally graded hollow ellipticcylinder subjected to uniform temperature distribution,Applications and Applied Mathematics 2017; 12 (1): 613632.
[26] Lamba NK, Walde RT, Manthena VR. Khobragade NW, Stress Functions in a Hollow Cylinder under Heating and Cooling, Journal of Statistics and Mathematics 2012; 3(3): 118124
[27] Kursun A, Kara E, Cetin E, Aksoy S, Kesimli A, Mechanical and Thermal Stresses in Functionally Graded Cylinders, World Academy of Science, Engineering and Technology International Journal of Aerospace and Mechanical Engineering 2014;8(2): 303308.
[28] Singru SS, Chauthale S, Khobragade NW, Thermal Stress Analysis of a Thick Hollow Cylinder, International Journal of Engineering and Innovative Technology (IJEIT) 2015; 5(3): 157160.
[29] Walde Anjali RC, Pathak C, Khobragade NW, Thermal Stresses of a Solid Cylinder with Internal Heat Source, Int. J. of Engg. And Information Technology (2013); 3(1): 407410.
[30] Jadhav CM, Khobragade NW, An Inverse Thermoelastic Problem of finite length thick hollow cylinder with internal heat sources, Advances in Applied Science Research (2013); 4(3): 302314.
[31] Khobragade, NW, Thermal stresses of a hollow cylinder with radiation type conditions, Int. J. of Engg. And Information Technology (2013); 3(5): 2532.
[32] Khobragade NW, Thermoelastic analysis of a solid circular cylinder, Int. J. of Engg. And Information Technology (2013): 3(5): 155162.
[33] Khobragade NW, Thermoelastic analysis of a thick hollow cylinder with radiation conditions, Int. J. of Engg. And Information Technology (2013): 3(4): 380387.
[34] Gahane TT, Khobragade NW, Transient Thermoelastic Problem of a Semiinfinite Cylinder with Heat Sources, Journal of Statistics and Mathematics (2012); 3(2): 8793.
[35] Hiranwar PC, Khobragade NW, Thermoelastic Problem of a Cylinder with Internal Heat Sources, Journal of Statistics and Mathematics (2012); 3(2): 8793.
[36] Khobragade NW, Khalsa L, Kulkarni A, Thermal Deflection of a Finite Length Hollow Cylinder due to Heat Generation, Int. J. of Engg. And Information Technology (2013); 3(1): 372375.
[37] Kursun a, Topcu M, Tetik T, Stress Analysis of Functionally Graded Disc under Thermal and Mechanical Loads, Physics Engineering(2011); 10: 29492954.
[38] Manthena VR, Kedar GD, Transient thermal stress analysis of a functionally graded thick hollow cylinder with temperaturedependent material properties,Journal of thermal stresses (2018); 41(5): 568582.