SciELO - Scientific Electronic Library Online

vol.13 issue3Vibration Analysis of a Magnetoelectroelastic Rectangular Plate Based on a Higher-Order Shear Deformation Theory author indexsubject indexarticles search
Home Pagealphabetic serial listing  

Services on Demand




Related links


Latin American Journal of Solids and Structures

Print version ISSN 1679-7817On-line version ISSN 1679-7825

Lat. Am. j. solids struct. vol.13 no.3 Rio de Janeiro Mar. 2016 


Hygrothermoelastic Analysis of Orthotropic Cylindrical Shells

J.S. Mohamed Alia 

Saleh Alsubarib 

Yulfian Aminandac 

a,b,c International Islamic University Malaysia (IIUM), Malaysia. a b c


In this work, the combined effect of moisture and temperature on the bending behavior of simply supported orthotropic cylindrical shells has been investigated. Initially three dimensional equilibrium equations of thermoelasticity, simplified to the case of generalized plane strain deformations in the axial direction are solved analytically for an orthotropic cylindrical shell strip under thermal loading. Based on the realistic variation of displacements from the elasticity approach, a new higher order shear deformation theory was proposed for the analysis of an orthotropic cylindrical shell strip under hygrothermal and mechanical loading. The zigzag form of the displacement is incorporated via the Murakami zigzag function. Results are presented for mechanical and thermal loading for various layups and they are validated against the derived elasticity solution. The significance of retaining various higher-order terms in the present model, in evaluating the stresses and deflection for composite laminates is brought out clearly through parametric study. Useful results for combined hygrothermal loading are presented in tabular and graphical form. It is expected that the numerical results presented herein will serve as bench mark in future.

Keywords: Hygrothermal; thermoelasticity; composite laminates; higher order shear deformation; zigzag function


Applications of the laminated composite structures have been increasing in the field of aerospace, automobile, civil and other engineering, as these structures possess advantages of high specific stiffness and high specific strength with a good tailoring capability. Aerospace structures such as wings or fuselages are mostly assemblies of shell structures, these structures are exposed to environmental conditions such as high temperature and moisture effect that are detrimental to the strength and stiffness of the composite structures. Hence the analysis of such composite shell structures under hygrothermal loading has been of considerable interest to researchers.

Classical shell theories originally developed for analysis of thin elastic shells based on the Love- Kirchhoff assumptions are not applicable for the analysis of laminated composite shells, as these theories neglect the effect of transverse shear deformations. Thus the application of such theories to laminated composite shells where shear deformation is very significant, such as thick laminates, may lead to errors in calculating deflections and stresses. Whereas higher order theories which account for shear deformation are useful to accurately model such laminated composite shells.

In order to establish the accuracy of any shell theory, it is important to compare the results with the benchmark results based on theory of elasticity, because they represent closer approximation to actual behavior of the structure than the commonly used engineering theories of bending. Thus such a solution can be deemed as a standard benchmark for determining the range of applicability of classical as well as other engineering theories. This has been the prime motivation for three dimensional analyses of beam, plate and shell structures. Such elasticity solutions are limited only for specific geometry and boundary conditions; hence one has to resort to approximate classical and higher order theories for complex loading and geometry.

Considerable work on the analysis of laminated composite shells has been reported by Noor et al. (1991) and Noor and Burton (1992). Other notable work on laminated shells under thermal/hygrothermal loading based on theory of elasticity/higher order models are as follows.

Huang and Tauchert (1991) & (1992) have presented elasticity solutions for cross-ply cylindrical panels and doubly curved cross-ply laminates. Burton and Noor (1994) presented three-dimensional thermoelasticity solutions for sandwich panels and cylindrical shells subjected to thermo- mechanical loads. Khdier (1996) presented thermoelastic analysis of cross-ply circular cylindrical shells using the higher order theory of JN Reddy. Ding and Tang (1998) have derived exact thermoelastic solution for an axisymmetric problem of thick closed laminated shells while Xia and Ding (2001) derived a three-dimensional thermoelastic solution for laminated cantilever cylindrical shell. Khare et al. (2003) presented closed-form solution of 2D higher-order shear deformation theories for the thermo-mechanical analysis of simply supported doubly curved cross-ply laminated shells. In addition, Sayman Onur (2005) presented the stress analysis of thick multi-layered composite cylindrical shells under hygrothermal loading using elasticity solutions based on the Lekhnitskii's theory for plane strain case.

Furthermore, Naidu and Sinha (2005) presented a finite element formulation based on FSDT for large deformation behavior of laminated composite cylindrical shell panels subjected to hygrothermal environments. Jacob and Senthil (2006) derived an exact solution for thermoelastic response of functionally graded orthotropic cylindrical shells. In another study, Brischetto and Carrera (2008) study the thermal stress analysis of composite shell by refined multilayered models. A study on the free-edge effects in laminated circular cylindrical shell panels subjected to hygrothermal loading is carried out by Nosier and Miri (2009) by utilizing displacement-based theories. A complete review on the analysis of laminated composite shells under thermal and hygrothermal loading from 2000 to 2010 has been recently reported by Qatu et al. (2012). Yifeng and Yu (2011) constructed a hygrothermoelastic model capable of performing a hygrothermoelastic stress analysis of composite laminates by using variational-asymptotic method. Later on, Lal et al. (2011) investigated the effect of random system properties on transverse nonlinear central deflection of laminated composite spherical shell panel subjected to hygro-thermo-mechanical loading using higher order theory proposed by JN Reddy. Brischetto and Carrera (2012) studied the static analysis of multilayered smart shells subjected to mechanical, thermal and electrical loads by using Carrera Unified Formulation (CUF).

Recently in 2013, Jafari et al. (2013) presented the thermo-elastic response of a 2-D functionally graded open cylindrical shell with temperature-dependent material properties based on third-order shear deformation plate theory of JN Reddy.

Very recently, Brischetto (2013) analyzed multilayered composite and sandwich shells under hygrothermal loads using refined two-dimensional models that has been developed based on Carrera Unified Formulation (CUF).

A careful study of the references on thermally loaded laminated shells based on elasticity reveals that most of them provide just the solutions but not useful numerical results in the form of deflections and stresses for various span to thickness ratio as has been reported by Pagano (1970b) & (1970a) for laminated plates under mechanical loading, Ren (1987), Bhaskar and Varadan (1993) for laminated shells under mechanical loading and Bhaskar et al. (1996) for plates under thermal loading. Taking care to note, that, these works has formed a well-established basis for assessment, as well as development of the methodology for various laminate theories and is most commonly referred to, one can recognize the absence of such important benchmark results for the flexure of thermally loaded orthotropic laminated shells. Thus the objective of this present study is to present benchmark results in a form suitable for the evaluation of various shell theories for orthotropic laminated shells under thermal loading. Having derived the benchmark results, then from the realistic variations of displacements based on the elasticity approach were used to develop an accurate higher order model for analysis of laminates under hygrothermal loading.


2.1 Elasticity Solution

A N-layered cylindrical shell strip that is infinitely long in the axial direction as shown Fig. 1 is considered. The fiber orientation of any layer is either parallel along θ or z direction. The equations governing the cylindrical bending of the shell are:

2.1.1 Equilibrium Equations

where a subscript comma denotes differentiation.

Figure 1: Geometry and coordinate system of shell 

2.1.2 Thermoelastic Stress-Strain Law

where Cij's and (i's are the stiffness coefficients and the coefficients of thermal expansion, respectively, transformed from the material co-ordinates to the geometric co-ordinates. T(r, θ) is the change in temperature with reference to the stress-free isothermal state.

2.1.3 Strain -Displacement Relations

where ur and uθ are the displacements in r and θ directions respectively.

Using Eqns. (2) and (3), the equilibrium equations (Eqn. (1)) can be written in terms of displacements as

where [L] and {t} are function of Cij's and (i's, and {u} is the displacement vector.

The temperature variation is taken as

The longitudinal edges are assumed to be simply-supported such that

2.1.4 Solution

A solution for Eqn. (4) satisfying the conditions of Eqn. (6), can be assumed as

where ϕr and ϕθ are functions of r. Substitution of this solution in Eqn. (4) yields a set of coupled fourth order system of non-homogeneous ordinary differential equations with variable coefficients, which can be further reduced to a set of ordinary differential equations with constant coefficients by a simple transformation [Wylie and Barrett (2012)]. Thus the final solution in terms of complementary solution and particular integral results in evaluation of four arbitrary constants per layer. Thus for an N-layer laminate there will be 4N unknown constants. These can be determined by enforcing the continuity of the displacements ur and uθ and the transverse stresses τ and σrr at the (N-1) interfaces of the laminate, along with the following lateral surface conditions

2.2 Higher Order Model

2.2.1 Definition of Displacement Field

Based on the results from elasticity solution, it can be proposed that the following complete model with 13 terms are required to accurately predict the behavior of a finite cross-ply shell under thermal loading

To demonstrate the capability of the proposed model, the formulation was made for the case of cylindrical bending similar to the elasticity solution with only v and w (uθand ur).

Consider an infinitely long laminated cylindrical shell panel. The panel coordinate system is such that 0 ≤ x ≤ ∞, 0 ≤ θ ≤ θm and -h/2 ≤ z ≤ h/2. The mean radius of the shell strip is considered as R and the load is uniform along the x axis so that the shell undergoes cylindrical bending (εx = 0).

The other displacement model assumed here is HSDT6 by neglecting the higher order terms in w which is described as:

where v and w are the displacements at any point in the laminate. The parameter vo is the in-plane displacement, and wo is the transverse displacement of a point on the middle plane. The functions Φθ is the rotation of the normal to the middle plane about θ. The other parameters λθ, ζθ, w1 and Γ are unknown higher-order terms which are function of θ only. ψk is the zigzag function term as defined by Murakami (Murakami, 1986)

where zk is a local transverse coordinate with its origin at the center of the kth layer and hk is the corresponding layer thickness.

The stress-stain law for the kth layer with respect to the (x, θ) coordinate is given by:

where the are the transformed elastic coefficients and ΔT = T - T0, ΔC = C - C0, in which, T(x, θ, z) and C(x, θ, z) are the temperature and moisture loads variations through the thickness respectively, whereas T0 and C0 are the reference temperature and moisture concentration respectively. αx, αθ, and αz are the thermal expansion coefficients and βx, βθ and βz are the moisture expansion coefficients respectively.

Expression for the principle of virtual work is

where q is the load per unit area on the mid-surface of the shell.

Integrating the displacement gradient given in Eq. (14) for the case of cylindrical bending by parts and by setting the coefficient of δv0, δw0, δØθ, δλθ, δζθ, δSθ, δw1 and δΓ to zero separately, we obtain the following equilibrium equations with associated boundary conditions.

The equilibrium equations for displacement models HSDT6 are the first six equations of displacement model HSDT8 given in Eq. (15)

Boundary conditions are given as

At θ = 0 constant: one from each of the following bracketed quantities should be specified:

The generalized stress resultants occurring in Eq. (15) are defined as follows:

The integration being from -hk/2 to hk/2 and the summation being carried out over the N layers of shells.

From (Eqns. 10-13), the stress resultants can be expressed in terms of displacements, hence the equilibrium equations can be written in terms of generalized displacements in operator form as:

where [L] is a matrix of differential operators with respect to θ, {δ} represent the displacement vector [v0 Øθ λθ ζθ Sθ w0 w1 Γ] and {f} is the load vector that depends on the mechanical and hygrothermal load.

The solution assumed for Eq. (18) expressed in terms of displacements which satisfies the above boundary conditions Eq. (16) is:

where P1 to P8 are the displacement amplitudes. Using Eq. (19) in Eq. (18) results in a system of linear algebraic equations and by solving these equations, unknown coefficient {δ} can be obtained readily and subsequently in-plane stresses and transverse stresses can be obtained. Although the transverse shear stresses can be calculated from the constitutive relations, these stresses may not satisfy the continuity conditions at the interface between layers, hence transverse shear stresses are obtained by integrating three dimensional equilibrium equations of elasticity.


In this section, initially, benchmark results based on elasticity are generated for cylindrical shell strips under mechanical loading to validate the analytical formulation of both the elasticity approach and the new higher order model developed. Then the results are presented for both mechanical and hygrothermal loads for various cross-ply layups and through the thickness variations of displacements and stresses are given for various R/h ratios.

All laminas are assumed to be of the same thickness and made of the same orthotropic materials with the following lamina properties (Bhaskar et al. (1996)).

Here L and T refer to the properties parallel and perpendicular respectively to the fiber.

All the results are generated for laminates with simply supported boundary condition such that at θ =constant

All reported results for thermal and hygrothermal loads are normalized using the following dimensionless expression except that for the mechanical loading which is stated in the corresponding Table.

3.1 Validation of the Laminate Model

For validating the formulation of the new higher order model HSDT8 of the orthotropic laminated shell, initially the laminate is considered to be under mechanical loading for which the elasticity results are regenerated using Ren's formulation (Ren (1987)). A simply supported laminated shell made of symmetric three-layer cross-ply (90/0/90) under mechanical loading was considered. The results for both the elasticity and higher order model HSDT8 are presented in Table 1 for comparison. The results for both stresses and displacements are found to agree very well even for thick laminates (R/h=4). It should be noted that the reported values of displacements and stresses based on elasticity has been generated using the elasticity formulation derived in this paper.

Table 1: Dimensionless maximum displacements and stresses of cross-ply (90/0/90) cylindrical shell strips subjected to mechanical loading 

* z values are given in parentheses; 𝜃 values are such that the sin or cos function become a maximum

As there are no benchmark results available in the literature for validation of higher order models for thermal loading; firstly, benchmark results based on the developed elasticity solution are generated for simply supported laminated shell strip under thermal loading. The results are presented in Table 2-4 for laminated shells made of (a) a unidirectional shell with the fibers oriented in the θ direction (90); (b) bidirectional laminated cylindrical shell with the T and L directions parallel to the θ-direction in the bottom and top layers respectively, the layers being of equal thickness (90/0); (c) a symmetric 3-ply cylindrical shell with layers of equal thickness-the L direction coincides with the θ-direction in the outer layers, while T is parallel to the θ-direction in the central layer (90/0/90).

Table 2: Dimensionless maximum displacements and stresses of (90) cylindrical shell strips subjected to thermal loading 

* z values are given in parentheses; 𝜃 values are such that the sin or cos function become a maximum

Table 3: Dimensionless maximum displacements and stresses of cross-ply (90/0) cylindrical shell strips subjected to thermal loading 

* z values are given in parentheses; 𝜃 values are such that the sin or cos function become a maximum.

Table 4: Dimensionless maximum displacements and stresses of cross -ply (90/0/90) cylindrical shell strips subjected to thermal loading 

* z values are given in parentheses; 𝜃 values are such that the sin or cos function become a maximum

For these problems the shell strip was considered to be subjected to temperature field linearly varying through the thickness (i.e.. T = 2(z/h) sin(mπθ/θm)). For to understand the capability of the HSDT8 in predicting the stresses and displacements accurately, results based on HSDT8 are also included in Table 2 and 3 for comparison and hence the validation of this higher order model under thermal loading. The results for both stresses and displacements are found to agree well even for thick laminates (R/h= 4) for all cases as shown in the following tables.

Through the thickness variations of the displacements and stresses in the (90/0/90) shell strip are presented in Figs 2 and 3. As can be seen from Fig.2, the variation of is small for R/h=10. However, for a thick laminate with R/h=4, varies significantly, thus while developing refined theories for thermal stress analysis of thick laminates, it is necessary to account for the thickness-stretch effect. Based on this observation from elasticity results, higher order terms have been considered for w in the higher order model (HSDT8). The variation of the in-plane displacement for thermal loading are shown in Fig 3, from which it can be noted that there is a considerable warping of the normal and a pronounced zig-zag effect for the case of the thick laminates. To show the importance of zig-zag term in the developed HSDT8 on the accurate prediction of in-plane displacement, through the thickness variation of inplane displacement with (HSDT8) and without zig zag term (HSDTWZ) in the Higher order model for both thermal and mechanical loading are graphed in Fig 4 and Fig 5 respectively. It can be noted from theses figures that zig-zag term is important in v for accurate prediction of displacements.

Figure 2: Varaition of deflection under thermal loading for (90/0/90) shell strip 

Figure 3: Varaition of displacement under thermal loading for (90/0/90) shell strip 

Figure 4: Varaition of displacement under thermal loading for (90/0/90) shell strip 

Figure 5: Varaition of displacement under mechanical loading for (90/0/90) shell strip 

To show the importance of the zig-zag function and the thickness stretch/contraction term, results obtained by dropping these terms from HSDT8 are also included in Figs 4-9 and Table 5. The thickness wise variations of the deflection, in-plane stresses and interlaminar stresses are given in Figs 6-9 for (90/0/90) cylindrical shell strip under thermal loading. From the above Tables and Figures it can be seen that not only the displacements but also the in-plane and the interlaminar stresses are accurate even for thick laminates under thermal loading. In all the above Tables, results and graphs for (or ) have not been presented. It is expected that once the transverse shear stresses are estimated accurately, their use with the integration of three-dimensional equilibrium equation would automatically yield accurate values of (or ).

Figure 6: Varaition of deflection under thermal loading for (90/0/90) shell strip 

Figure 7: In-plane stress under thermal for (90/0/90) shell strip 

Figure 8: In-plane stress under thermal loading for (90/0/90) shell strip 

Figure 9: Transverse shear stress under thermal loading for (90/0/90) shell strip 

Table 5: Dimensionless maximum deflection (-h/2)* of (90/0/90) shell strip under thermal load 

R/h HSDT8 HSDT6 CST Elasticity
4 19.80 2.147 3.145 19.754
10 5.869 2.950 3.145 5.860
20 3.830 3.096 3.145 3.827
50 3.257 3.139 3.145 3.255
100 3.175 3.145 3.145 3.171

* z values are given in parentheses; 𝜃 values are such that the sin or cos function become a maximum

Table 5 compares the results for deflection obtained by HSDT8, HSDT6 (HSDT8 without the higher order terms in w) and CST with results obtained from 3-D elasticity solution. It can be seen that HSDT6 demonstrates a marked accuracy deterioration in the deflection especially at low values of aspect ratio (R/h<10). Thus HSDT6 as well as CST, without the thickness stretching effect being considered, ends up with a very poor estimation of deflection for low values of R/h.

Thus from above Tables 1-5 and Fig 4-9, it is clear that for mechanical/thermal loading, the displacements and stresses obtained using HSDT8 which account for both thickness stretch effect in w and zig-zag function in v agree very well with those obtained by 3D elasticity solution even for very thick laminate with R/h = 4.

3.2 Results for Hygrothermal Loading

Having validated the higher order model HSDT8, a study on the effect of combined loading-thermal and hygroscopic load on the bending response of simply supported (90/0/90) laminated cylindrical shell strip is analyzed.

A linear variation of moisture and thermal distribution across the shell thickness is considered.

where and are the temperature and moisture concentration that are assumed to be 300ºK and 0.01% respectively, where θm is taken to be 1 radian and m=1.

Also the thermal and moisture expansion coefficients are assumed to be (Zenkour et al., 2014):

αT = 22.5 * 10-6/ºK; αL = 0.02 * 10-6/ºK

βL = 0; βT = 0.6 (wt. %H2O)-1

Table 6 shows the results for deflections of (90/0/90) laminated cylindrical shell strip under combined loads (hygrothermal). When comparing with the results from Table 4, one can see that the hygrothermal load affects the deflection and stresses more than the thermal loading and in some cases, it is almost twice the value of thermal case, as reported for laminated plates (Zenkour et al., 2014).

Table 6: Dimensionless maximum displacements and stresses of cross-ply (90/0/90) cylindrical shell strips subjected to hygrothermal loading 

* z values are given in parentheses; 𝜃 values are such that the sin or cos function become a maximum

Fig 10-13 shows the thickness-wise variation for deflections and stresses of a (90/0/90) laminated cylindrical shell strip under both thermal and hygrothermal loading for a very thick laminate(S=4). When comparing the results for hygrothermal loading with the results from thermal loading, one can see that the hygrothermal load affects the deflection and stresses more than the thermal loading and it can be seen that the stresses increases prominently at the interfaces and in some cases it is almost twice the value of thermal case.

Figure 10: Varaition of deflection under thermal and hygrothermal loading for (90/0/90) shell strip 

Figure 11: In-plane stress under thermal and hygrothermal loading for (90/0/90) shell strip 

Figure 12: In-plane stress under thermal and hygrothermal loading for (90/0/90) shell strip 

Figure 13: Transverse shear stress under thermal and hygrothermal loading for (90/0/90) shell strip 


Thermoelastic solutions for orthotropic cylindrical shells under thermal loading have been presented. Useful benchmark results based on elasticity have been presented for cross-ply laminated cylindrical shell strip. A new accurate higher order theory (HSDT8) based on realistic variations of displacements from the elasticity has been proposed. The results of the proposed model were validated against 3D elasticity solution and found to predict well both the stresses and displacements, even for a thick laminate. Numerical results are presented for various cross-ply layup of cylindrical laminated shell strips under thermal loading. The graphical and tabular results presented bring out clearly the importance of higher order terms in the present theory HSDT8. It was found that for thick laminated shells, inclusion of the thickness stretch/contraction effects are more important only for thermal loading whereas the accounting for the interface slope discontinuity (zig-zag function) of the in-plane displacements is essential for both mechanical and thermal loading. It is believed that the results presented here will serve as benchmark in future.


The authors would like to thank the Ministry of Higher Education, Malaysia for supporting this work through the Research Grant FRGS 12-073-0222.


Ahmed K. Noor, W. Scott Burton and Peter Geame M. (1991). Assessment of computational models for multilayered composite cylinders. International Journal of Solids and Structures, Volume 27(10), 1269-1286. [ Links ]

Ahmed Noor K and Scott Burton W. (1992). Computational Models for High-Temperature Multilayered Composite Plates and Shells. Appl. Mech., Rev. 45(10), 419-446. [ Links ]

Bhaskar, K., & Varadan, T. K. (1993). Elasticity solution for laminated anisotropic cylindrical shells. Journal of Applied Mechanics, 6(March 1993), 1-7. [ Links ]

Bhaskar, K., Varadan, T. K., & Ali, J. S. M. (1996). Thermoelastic solutions for orthotropic and anisotropic composite laminates. Composite :part B, 8368(96), 415-420. [ Links ]

Brischetto, S. (2013). Hygrothermoelastic analysis of multilayered composite and sandwich shells. Journal of Sandwich Structures and Materials, 15(2), 168-202. [ Links ]

Brischetto, S., & Carrera, E. (2008). Thermal Stress Analysis by Refined Multilayered Composite Shell Theories. Journal of Thermal Stresses, 32(1-2), 165-186. [ Links ]

Brischetto, S., & Carrera, E. (2012). Static analysis of multilayered smart shells subjected to mechanical, thermal and electrical loads. Meccanica, 48(5), 1263-1287. [ Links ]

Huang, NN. and Tauchert TR. (1992). Thermal stresses in doubly-curved cross-ply laminates. International Journal of Solids and Structures, 29(8), 991-1000. [ Links ]

Huang, NN, & Tauchert, TR. (1991). Thermoelastic solution for cross-ply cylindrical panels. Journal of Thermal Stresses, Volume 14(2), 37-41. [ Links ]

Kewei, D. & Limin, T. (1998). Exact Thermoelastic Solution for an Axisymmetric Problem of Thick Closed Laminated Shells. Journal of Thermal Stresses, 21(7), 751-761. [ Links ]

Khare, R. K., Kant, T., & Garg, A. K. (2003). Closed-form thermo-mechanical solutions of higher-order theories of cross-ply laminated shallow shells. Composite Structures, 59(3), 313-340. [ Links ]

Khdier, A. (1996). Thermoeslastic analysis of cross-ply laminated circular cylinderical shells. International Journal of Solids and Structures, 33, 4007-4017. [ Links ]

Lal, A., Singh, B. N., & Anand, S. (2011). Nonlinear bending response of laminated composite spherical shell panel with system randomness subjected to hygro-thermo-mechanical loading. International Journal of Mechanical Sciences, 53(10), 855-866. [ Links ]

Mehrabadi, S. J., & Aragh, B. S. (2013). On the thermal analysis of 2-D temperature-dependent functionally graded open cylindrical shells. Composite Structures, 96, 773-785. [ Links ]

Murakami, H. (1986). Laminated Composite Plate Theory With Improved In-Plane Responses, 53(September 1986), 661-666. [ Links ]

Naidu, N. V. S., & Sinha, P. K. (2005). Nonlinear finite element analysis of laminated composite shells in hygrothermal environments. Composite Structures, 69(4), 387-395. [ Links ]

Nosier, A., & Miri, A. K. (2009). Boundary-layer hygrothermal stresses in laminate composite circular cylindrical shell panels. Archive of Applied Mechanics, 80(4), 413-440. [ Links ]

Pagano, N. J. (1970a). Exact solutions for rectangular bidirectional composites and sandwich plates. Journal of Composite Materials, 4, 20-34. [ Links ]

Pagano, N. J. (1970b). Influence of Shear Coupling in Cylindrical Bending of Anisotropic Laminates. Journal of Composite Materials, 4(3), 330-343. [ Links ]

Pelletier, J. L., & Vel, S. S. (2006). An exact solution for the steady-state thermoelastic response of functionally graded orthotropic cylindrical shells. International Journal of Solids and Structures, 43(5), 1131-1158. [ Links ]

Qatu, M. S. (2012). Review of Recent Literature on Static Analyses of Composite Shells: 2000-2010. Journal of Composite Materials, 02(03), 61-86. [ Links ]

Ren, J. G. (1987). Exact solutions for laminated cylindrical shells in cylindrical bending. Composites Science and Technology, 29(3), 169-187. [ Links ]

Sayman, O. (2005). Analysis of multi-layered composite cylinders under hygrothermal loading. Composites Part A: Applied Science and Manufacturing, 36(7), 923-933. [ Links ]

W. Scott Burton, A. K. N. (1994). Three-dimensional solutions for thermomechanical stresses in sandwich panels and shells. Journal of Engineering Mechanics, Vol. 120, pp. 2044-2071. [ Links ]

Wylie, R. C. and Barrett. L. C. (2012). " Advanced Engineering Mathematics." Tata Mc Graw Hill Education Pvt Ltd, New Delhi. [ Links ]

Xia, L., & Ding, K. (2001). Three-dimensional thermoelastic solution for laminated cantilever cylindrical shell. Aerospace Science and Technology, 5(5), 339-346. [ Links ]

Yifeng, Z., & Yu, W. (2011). A variational asymptotic approach for hygrothermal analysis of composite laminates. Composite Structures, 93(12), 3229-3238. [ Links ]

Zenkour, A. M., Mashat, D. S., & Alghanmi, R. A. (2014). Hygrothermal analysis of antisymmetric cross-ply laminates using a refined plate theory. International Journal of Mechanics and Materials in Design, 10(2), 213-226. [ Links ]

Received: June 26, 2015; Revised: January 02, 2016; Accepted: January 05, 2016

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License