Analytical Bending and Stress Analysis of Variable Thickness FGM Auxetic Conical/Cylindrical Shells with General Tractions

Abstract

In the present research, a unified semi-analytical solution that incorporates influence of the auxeticity (negative Poisson ratio) of the material into elastic responses of the variable thickness functionally graded conical and cylindrical shells and circular/annular plates is developed. The top or bottom layers of the shell/plate may be subjected to general non-uniform normal and shear tractions. The mentioned material and loading complexities have not been investigated before and consequently, the presented comprehensive results are quite new. The proposed unified formulation is developed using the principle of minimum total potential energy and solved using Taylor’s transform, for some combinations of the simply supported, clamped and free edge conditions. Accuracy of results of the proposed unified solution is verified by results of the three-dimensional theory of elasticity extracted from the ABAQUS finite element analysis code. Finally, a comprehensive parametric study including evaluation of individual/simultaneous effects of the auxeticity, structure configuration, shear/normal traction, thickness variability, and boundary conditions on the resulting lateral deflection and in-plane stress distributions of the considered shell and plate structures is accomplished.

Keywords:
Taylor’s transform; Bending and stress analysis; Auxetic FGM conical and cylindrical shells; Non-uniform shear and normal tractions

1 INTRODUCTION

Simpler manufacturing requirements and more geometrical consistency with the nature of the applied loads have led to extensive use of the shells of revolution in many mechanical, civil, aerospace, vehicular, and other engineering structures; especially as rotating, pressure containing, or aerodynamic components. These structures may sometimes be fabricated from functionally graded materials (FGMs) to benefit from continuously varying volume fractions of the different material phases of the single construction instead of using the multilayered structures that exhibit jumps in the material properties at the interfaces between layers. Circular/annular plates and cylindrical/conical shells are among the most common axisymmetric structures. Comprehensive descriptions of the theories presented so far on isotropic, composite, and FGM plates and doubly-curved shells may be found in books published by Reddy (2003), Shen (2009), Carrera et al. (2011), and Tornabene and Fantuzzi (2014). The cylindrical or conical shells (e.g., CNG pressure vessels and conical caps) produced by deep drawing or spinning of the billets or designed to meet specified criteria, are generally variable-thickness. On the other hand, the conical caps or thin cylindrical components may experience normal and shear tractions; e.g., when they are subjected to aerodynamic, hydrodynamic, or internal/external fluid flow forces.

In contrast to the extensive researches accomplished on displacement analysis, vibration (Nie and Zhong, 2007; Tornabene, 2009; Setoodeh et al., 2011; Shariyat and Alipour, 2011; Shariyat, 2012; Tahani et al., 2013; Duc, 2013; Mercan et al., 2016) or buckling (Sofiyev, 2007; Shariyat, 2008; Alipour and Shariyat, 2011, 2013; Shariyat and Asgari, 2013; Sun et al., 2014) of the FGM circular/annular plates and cylindrical/conical shells, limited works have been published on the stress analysis of the mentioned structures. A semi-inverse series elasticity solution was presented by Pelletier and Vel (2006) for thermoelastic analysis of a transversely graded orthotropic cylindrical shell. An mth-order shear-deformation theory was presented by Matsunaga (2009) for vibration and buckling analyses of the transversely graded cylindrical shells. Viola et al. (2012) conducted a parametric investigation coupled with a stress recovery procedure for the FG cylindrical shells. A high-order theory was presented by Zozulya and Zhang (2012) for the axisymmetric FGM cylindrical shells. Asymmetric bending of the FGM circular plates was investigated by Nosier and Fallah (2009), using the first-order shear-deformation theory (FSDT). Wang et al. (2010) presented a three-dimensional solution for axisymmetric bending of FGM circular plates. Sburlati and Bardella (2011) developed three-dimensional elastic solutions for the axisymmetric FGM circular plates. Behravan Rad and Shariyat (2013) presented a 3D elasticity analysis for bidirectional FGM annular plates subjected to normal and shear tractions, on elastic foundations. Golmakani and Kadkhodayan (2014) investigated axisymmetric bending and stretching of FGM circular/annular plates with variable thicknesses under thermomechanical loads, using the FSDT and the finite difference technique. Shariyat and Mohammadjani (2013,2014) conducted 3D elasticity stress and displacement analyses of rotating bidirectional FGM annular plates subjected to non-uniform loads and non-uniform elastic foundations. Alipour and Shariyat (2012) and Shariyat Alipour (2014) incorporated 3D elasticity corrections into results of the single-layer/sandwich plate theories of circular/annular plates under non-uniform shear and normal tractions. Recently, Behravan Rad and Shariyat (2015) conducted a three-dimensional magneto-elastic analysis for asymmetric variable thickness porous FGM circular plates with non-uniform shear and normal tractions and elastic foundations.

In the stress analysis field of the conical shells, Asemi et al. (2011) furnished a finite element elastic analysis for a bidirectional functionally graded thick truncated cone under combined internal, external, and axial pressures. Aghdam et al. (2011) conducted a bending analysis for FGM conical panels subjected to non-uniform distributed loadings, using the FSDT and solving the governing equations by the extended Kantorovich method. Setoodeh et al. (2012) studied transient dynamic behavior of transversely graded axisymmetric truncated variable thickness conical shells subjected to internal/external pressure shocks. Abediokhchi et al. (2013) presented a bending analysis for the FGM conical panels based on the FSDT and generalized differential quadrature (GDQ) solution method. The FSDT and the perturbation theory were employed by Ghannad et al. (2013) for deformation and stress analysis of axisymmetric clamped variable thickness FGM cylindrical shells under internal pressures. Tornabene et al. (2013) applied the GDQ Method, FSDT, and stress recovery approach for stress analysis of transversely graded and laminated composite shells of revolution. Sofiyev (2014) investigated large deformation dynamic buckling of truncated conical shells with functionally graded coatings under axial loads. Static analysis of FGM truncated conical shells subjected to meridian, circumferential and normal uniform loadings was conducted by Viola et al. (2014), using a third-order shear-deformation theory, GDQ method, and the stress recovery approach. A hybrid method composed of the layerwise theory, DQM, and Fourier series expansion was developed by Selahi et al. (2014) for transient analysis of FGM truncated conical shells with variable thickness, under asymmetric dynamic pressures. Tornabene et al. (2015) investigated recovery of the static transverse normal and shear stresses in FGM doubly-curved sandwich shells of revolution using Carrera’s unified formulation and the DQM. Using the FSDT and dividing the FGM shell into homogeneous disks, and using the finite element method, Zamani Nejad (2015) performed elastic analysis of rotating axially graded truncated conical shells under non-uniform pressures.

Very rare articles may be found in literature on elastic analysis of the auxetic structures. Khan and Hilton (2009) demonstrated that for the non-homogeneous linear elastic materials, the Poisson's ratio becomes both space and time dependent under time dependent stresses. Hilton et al. (2009) presented the idea of optimum 3-D anisotropic designer materials, including auxetic viscoelastic functionally graded ones. Brighenti (2014) studied effects of the auxetic materials on lateral deflections of isotropic and homogeneous layered circular plates under uniform transverse loads, using the classical bending plate theory. Performance of sandwich panels with auxetic lattice cores confined between metallic facets were studied by Imbalzano et al. (2015) under localised impacts at high strain rates and a remarkable reduction was observed in the maximum back facet displacements, due to the localised densification of the auxetic core. Based on Reddy’s first-order shear deformation plate theory, Duc and Cong (2016) analyzed nonlinear dynamic response and vibration of sandwich plates with auxetic honeycomb cores on elastic foundations subjected to blast and mechanical loads. The only published paper on stress analysis of the auxetic structures was due to Alipour to Shariyat (2015) who developed analytical zigzag solutions with 3D elasticity corrections for bending and stress analysis of circular/annular composite sandwich plates with auxetic cores.

In the present research, a unified analytical formulation for stress and displacement analysis of conical and cylindrical shells and annular plates with variable thickness is developed, using the FSDT with 3D elasticity correction. It is the first time that effects of using auxetic functionally graded materials on distributions of the displacements and stresses of the conical and cylindrical shells is investigated. Therefore, the reported comprehensive results are quite new. Both the top and bottom layers of the shell/plate may be subjected to non-uniform normal and shear tractions. Edges of the shell/plate may be simply supported, clamped or free. The proposed unified formulation is derived based on the principle of minimum total potential energy and solved analytically. Accuracy of results of the proposed Taylor transform solution is verified by those of the three-dimensional theory of elasticity obtained from the ABAQUS finite element analysis code.

2 THE UNIFIED FORMULATION

2.1 Description of Variations of the Thickness, Material Properties, and Normal and Shear Tractions

Geometric and load parameters of the considered general shell are illustrated in Fig. 1. The two limit values of the semi vertex angle (α), i.e., α = 0 and α = π/2 (which correspond respectively, to the φ = π/2 and φ = 0 meridian angles) are associated respectively, with the cylindrical shell and annular plate configurations.

Figure 1:
Geometric parameters and coordinate system of the considered variable thickness functionally graded auxetic shell with arbitrary distributions of the normal and shear tractions.

As may be noted from Fig. 1, thickness of the shell may vary in a general manner; so that distances of the top and bottom layers of thickness may vary arbitrary, e.g., as follows, with respect to the reference surface of the shell, in the meridian direction:

(1)

where the subscripts t and b denote the top and bottom surfaces of the shell, respectively, and the coordinate z is measured from the reference surface of the shell (Fig. 1):

(2)

So that the local thickness (h) of the shell may be defined as:

(3)

Moreover, the top and bottom layers of the thickness may be prone to normal and shear tractions whose intensity vary arbitrary along the meridian direction, e.g., according to the following parabolic forms (Fig. 1):

(4)

where the symbols q and τ stand for the normal and shear tractions and. ,q are amplitudes of the shear and normal tractions, respectively, and (αt, βt, γt), (αb, βb, γb)( η, λ, ξ) are the relevant loading coefficients.,

It is assumed that the conical shell (or equivalently, the cylindrical shell or annular/circular plate) are fabricated from a mixture of ceramic and metal constituent materials whose material compositions vary continuously in the thickness direction. Some researchers have shown that results of Voigt’s rule of mixtures and Mori-Tanaka micromechanical homogenization techniques are generally close (Ferreira et al., 2005; Alipour and Shariyat, 2015), unless constituent materials with quite different Poisson ratios are employed (Shariyat and Asemi, 2014). On the other hand, one of motivations of the present research is investigation of effects of the auxicity nature of the material properties through comparative analyses. Therefore, both models may be used in the mentioned comparative analyses. According to the rule of mixture, transverse variations of the effective elastic moduli of the shell/plate may be expressed as follows:

(5)

where g is the positive definite volume fraction index and the subscripts c and m refer to the ceramic and metallic constituent materials, respectively. The ceramic volume fraction, _^Vc , is so defined that, the top layer (z = _^ht (x)) becomes a metal-rich (_^Vc = 0) and the bottom layer (wherein, z = -_^hb (x) = _^ht (x) -h) becomes a ceramic-rich (_^Vc = 1) one.

2.2 The Governing Equations

According to the first-order shear-deformation theory, the displacement field of the axisymmetric shell may be described as follows:

(6)

where _^u0 and w are the meridian and transverse displacement components of the middle surface, respectively, and _^ψx is rotation of the meridian section of the shell/plate.

The strain-displacement relations of the conical shell are:

(7)

where _^εx and ε _{θ ?}are the normal strain components and _^γxz is the transverse shear strain. _^R0 is the radius of the cone at a plane perpendicular to cone axis:

(8)

where _^Rb corresponds to the starting edge of the truncated conical shell. The corresponding stresses may be obtained based on the generalized Hooke’s law as:

(9)

The governing equations of the shell/plate may be derived by using the principle of minimum total potential energy:

(10)

where δ∏, δU, and δW are increments of the total potential energy, strain energy, and work of the externally applied loads, respectively:

(11)

where A and V stand for the area and volume and κ is the shear correction factor, e.g., κ - 5/6, and:

(12)

In the first-order shell/plate theories, the assumption of zero or constant transverse shear stresses is amended through considering work of the transversely distributed loads (normal traction) of the shell/plate. The solution may be enhanced further by enforcing the traction conditions _{^{τ = τt}} and _{^{τ = τb}} at the top and bottom surfaces of the shell/plate. In contrast to the elasticity and equilibrium equations approaches, in the energy version of the three-dimensional theory of elasticity (that is the base of almost all the commercial finite element analysis codes), the traction or force boundary conditions may be imposed through including of the relevant works. Therefore, distributions of the stresses and displacements so arranged that they finally become compatible with the mentioned traction conditions, through the principle of minimum potential energy. This fact is valid even for the stress fields formed around the concentrated loads of the structures.

Substituting Eq. (7) into Eq. (9) and substituting the resulting equations into Eq. (11) enables expressing of Eq. (10) in terms of the displacement parameters and their increments (δu ₀. ^{_{δψx, δw}} ):

(13)

where, e.g., A ≡ A (_{^{u0, ψx, w, φ, q,}}, , η, λ, ξ, αt, βt, γt, αb, βb, γb, g, Ec, Em, at, bt, ct, ab, bb, cb, x0). Since the displacement increments (δu0. δψx, δw) are generally non-zero values for the interior points of the shell/plate, it may be deduced that their coefficients appeared in Eq. (13), i.e., 𝒜, ℬ, and 𝒞 are zero for those points. Therefore, the governing equations of the shell become:

(14)

The stress resultants M, N, Q are:

(15)

where based on Eqs. (6), (7), and (9), one has:

(16)

or in expanded form:

(17)

Definition of the coefficients appeared in the right hand side of Eq. (17) is not included here, to save space. In this regard, Eq. (1) is employed as well.

The governing equations (14) may be rewritten based on Eqs. (15), as follows:

(18-a)

(18-b)

(18-c)

2.3 Boundary Conditions

Some more common edge conditions are considered here to develop the semi-analytical solution:

(19-a)

(19-b)

(19-c)

3 THE UNIFIED TAYLOR TRANSFORM SOLUTION

According to the Taylor transform for the finite domains (Shariyat and Alipour, 2011; Alipour and Shariyat, 2011; Shariyat and Asemi, 2014), the displacement components may be expanded in series forms about the initial value of the meridian coordinate:

(20)

Therefore, one may define the following expressions that have been appeared in Eq. (18):

(21)

where x = x ₀ is associated with a proper center of expansion. In the present research, x ₀ = x is chosen. By substituting Eqs. (20) and (21) into the governing Eq. (14) and performing some manipulations, the transformed form of Eq. (14) becomes:

(22-a)

(22-b)

(22-c)

The expanded form of the displacement components introduced in Eq. (20) may be used to express the stress resultants appeared in the boundary conditions (19) by Taylor’s expansions in terms of the displacement series coefficients:

(23)

By solving Eqs. (22-a) to (22-c) for i=0...N, an algebraic system of 3(N+1) equations in terms of the coefficients of Taylor’s series of the unknown displacement parameters, _^Ui , Ψ_i and _^Wi (i=2...N+2) is resulted and may be solved in terms of the 6 unknown displacement parameters U ₀, U ₁, Ψ₀, Ψ₁, W ₀ and W ₁. These six displacement parameters may be determined by substituting the unknown displacement parameters _^Ui , Ψ_i and _^Wi (i=2...N+2) in terms of U ₀, U ₁, Ψ₀, Ψ₁, W ₀ and W ₁, into the six transformed edge conditions of the shell/plate [3 edge conditions for each individual edge, according to Eq. (19)]; so that the resulting system of equation becomes:

(24)

where

(25)

4 RESULTS AND DISCUSSIONS

In the present research, the results are verified by the three-dimensional elasticity results extracted from ABAQUS finite element analysis code, whenever required. In this regard, number of the terms of the Taylor transformation series is so chosen that almost no change is observed by increasing this number (50 terms are chosen for each of the displacement components). Similarly, number of the 8-node biquadratic CAX8R axisymmetric solid elements used in ABAQUS code for verification of the results is adopted based on the same criteria. It was noticed that adopting about 3000 elements (where three elements are used in the transverse direction) leads to convergent results.

Moduli of materials of the top and bottom layers of the shell/plate are chosen as follows:

_^Et = 70GPa, _^Eb = GPa

But the Poisson ratio is varied to evaluate effects of the auxeticity of the material. The top and bottom layers of the shell/plate are respectively, metallic- and ceramic-rich layers. The ceramic face of the structure is subjected to normal and shear traction. Since effects of the material index have been investigated by numerous researchers and its qualitative effects may be deduced without the need to perform numerical analyses, it is chosen as g=1 which corresponds to linear variations of the volume fractions of the constituent materials (a more practical case). _^Rb and L are chosen equal to 1m. Furthermore three values are considered for the semi-vertex angle, α = π / 2, π / 4,0 which correspond to annular plate and conical and cylindrical shells, respectively.

4.1 Influence of the Auxeticity on the Stress and Deflections of Uniform Thickness FGM Shells/Plates Under Shear Loading and Various Boundary Conditions

Less researches have been published on plates or shells subjected to shear tractions, especially, conical shells. On the other hand, the stresses and displacements induced by shear tractions are significantly dependent on the shear modulus and consequently, Poisson ratio of the structures, particularly those fabricated from heterogeneous materials, e.g., FGMs.

The considered structures are prone to a shear traction of magnitude h/L=0.1. Denoting the simply supported, clamped, and free edge conditions by the abbreviations S, C, and F, respectively, three types of boundary conditions are considered here: = 1MPa. Thickness of the structure is assumed here to be uniform and equal to

1. C-C boundary condition,

2. S-C boundary condition,

3. F-C boundary condition.

where the first and second letters are associated with the edges located at x=0 and L, respectively.

In the first stage, an annular plate is considered. Influences of the auxeticity of the material are investigated in Figs. 2 to 4 on radial distributions of lateral deflection of plates having the mentioned three boundary conditions, respectively. Magnitude of auxeticity of the material depends on the manufacturing process. In the present research, a practical range of the Poisson ratio is scanned with equal steps, specifically, behavior of the structure is evaluated for ν = 0.3,0.-0.3, -0.6 and -0.9. Figs. 2 to 4 reveal that auxeticity of the material significantly affects overall stiffness and consequently, lateral deflection of the plate. Results are verified by those obtained from ABAQUS code, for ν = 0.3 and -0.9. As Figs. 2 to 4 show, the differences between results of the successive values of the Poisson ratio grow as the Poisson ratio becomes more negative; so that the for ν = -0.9, the structure becomes extremely rigid.

Figure 2:
Influence of the auxeticity of the material on the radial distribution of the lateral deflection of the C-C annular plate, under shear traction.

Figure 3:
Effects of the auxeticity of the material on the radial distribution of the lateral deflection of the S-C annular plate, under shear traction.

Figure 4:
Effects of the auxeticity of the material on the radial distribution of the lateral deflection of the F-C annular plate, under shear traction.

As Fig. 1 shows, a positive shear traction on the top layer induces a bending moment in the clockwise direction with respect to the lower layers; for this reason, an in-ward lateral deflection may be induced in the conical shell (or equivalently, a downward lateral deflection in the annular plate), as Fig. 2 (and later, Fig. 14) implies. In addition to the induced moment, the shear traction tends to move particles of the top layer of the shell/plate. Due the kinematic of deformation of the plate/shell and the constituted curvature, slopes of the lateral deflection curve becomes negative and positive at the inner/beginning and outer/end regions of the plate/shell, respectively. As the shear traction tends to move the top layer along the tangent to the deformed plate/shell; it tends to move the plate/shell downward/inwards in the vicinity of the smaller edge of the plate/shell (due to the negative slope) whereas it pushes the top layer upward/outwards (due to the positive slope), in the neighborhood of the larger edge of the annular plate or conical shell. When the Poisson ratio is positive, the plate/shell is more compliant; so that the plate/shell may experience a small upward/outward deflection in the neighborhood of the larger edge, as present and ABAQUS results of Fig. 2 (and later, Fig. 14) confirm. When material of the plate/shell is more auxetic, the plate/shell becomes stiffer; so that, the force effect of the shear traction diminishes in comparison to influence of the bending moment of the sections.

To present a more general imagination regarding effects of the auxeticity of the material on the transverse distributions of the resulting in-plane stresses, asymmetric sections located at x=0.3 and 0.8m (where x denotes the radial distance from the inner edge of the plate) are considered. Results are plotted in Figs. 5 to 7, for the three types of boundary conditions, respectively. These Figures show that the shear traction has produced bending moments about the reference surface of the plate. However, since transverse distribution of the modulus of elasticity of the section is not uniform, the resulting transverse distributions of the in-plane stress are non-linear. Results show an excellent concordance with the ABAQUS 3D-elasticity results.

Figure 5:
Effects of the auxeticity of the material on transverse distribution of the in-plane stress of the C-C annular plate, under shear traction.

Figure 6:
Effects of the auxeticity of the material on transverse distribution of the in-plane stress of the S-C annular plate, under shear traction.

Figure 7:
Effects of the auxeticity of the material on transverse distribution of the in-plane stress of the F-C annular plate, under shear traction.

The displacement results obtained either by present research or from ABAQUS code (some of them are not reported here due to their coincidence with the presented analytical results) show that trends of variations of results are different for the ν = - 0.9, ..., 0 and ν = 0, ..., 0.3 ranges of the Poisson ratio. On the other hand, the bending mechanism is quite different for the three types of boundary conditions, e.g., for the clamped-clamped and free-clamped plates. The stress results show that as the movability and therefore, rigidity of the plate increase, effects of the auxeticity of the material become more pronounced; so that results of the C-C and F-C plates exhibit respectively, the most and the least sensitivity to the auxeticity of the material. Furthermore, in contrast to the somewhat linear transverse distribution of the in-plane stress of the C-C and S-C plates in the neighborhood of the clamped edge, this distribution is remarkably distorted in the F-C plate.

Axial distributions of the lateral deflection of the auxetic FGM cylindrical shell are plotted for the three types of the boundary conditions in Figs. 8 to 10, respectively. Since the resulting moments due to the imposed shear traction, geometry, and boundary conditions are symmetric with respect to the mid-length section of the C-C shell, the inflection point has almost formed at that section of the cylindrical shell. As may be noted from Fig. 8, trends of variations of the axial distribution of the lateral deflection are different for the ν = - 0.9, ..., 0 and ν = 0, ..., 0.3 ranges of the Poisson ratio; so that even direction of the lateral deflection has changed for the negative Poisson ratios. As before, results are verified by ABAQUS code for the lower and higher limits of the Poisson ratio range and a good agreement is noticed. The same phenomenon has occurred for the S-C shell (Fig. 9). As Fig. 9 shows, location of the maximum lateral deflection has also changed and shifted toward the mid-length section of the S-C shell for more negative Poisson ratios; i.e., the difference between the edge conditions of the S-C shell has been compensated to a great extent by the overall rigidity induced by auxeticity of the material. As may be readily seen in Fig. 10, for the F-C plate, the shear traction has brought about outwards movements in the neighborhood of the clamped edge but inwards movements of the shell particles in the vicinity of the free edge, for low Poisson ratios. The auxeticity of the material mixture has significantly affected the resulting curvatures and the trend of the lateral deflections.

Figure 8:
Influence of the auxeticity of the material on axial distribution of lateral deflection of a C-C cylindrical shell subjected to a shear traction.

Figure 9:
Influence of the auxeticity of the material on axial distribution of lateral deflection of a S-C cylindrical shell subjected to a shear traction.

Figure 10:
Influence of the auxeticity of the material on axial distribution of lateral deflection of a F-C cylindrical shell subjected to a shear traction.

Transverse distributions of the in-plane stress of the shell are illustrated for various Poisson ratios in Figs. 11 to 13, for the three types of the boundary conditions, respectively. It is worth noting that due to presence of the axial shear traction, no axial plane of symmetry may be defined here for the stress distribution. It is expected that higher lateral deflections and stresses form in a shell with more movability (the F-C shell) and the opposite is also true (C-C shell has the smallest deflections and stresses). Furthermore, effects of the auxeticity of the material have higher effects on regions located at the neighborhood of the first edge (with generally more movability); although the resulting complicated deformation and stress fields and their interactions with the negative Poisson ratio effects may not enable exactly predictions of trends of the variations. Some of the ABAQUS results are also included to confirm that the results are accurate enough.

Figure 11:
Influence of the auxeticity of the material on transverse distribution of the in-plane stress of a C-C cylindrical shell subjected to a shear traction.

Figure 12:
Influence of the auxeticity of the material on transverse distribution of the in-plane stress of a S-C cylindrical shell subjected to a shear traction.

Figure 13:
Influence of the auxeticity of the material on transverse distribution of the in-plane stress of a F-C cylindrical shell subjected to a shear traction.

Fig. 14 that illustrates meridian distribution of the lateral deflection of the C-C conical shells, show that rigidity of the auxetic shell is larger in the vicinity of the smaller edge; as effects of the initial clamped edge (retaining the zero slope condition) on the interior regions have decayed with a smaller rate, so that these effect last further and location of the maximum lateral deflection shifts toward the larger edge. However, for positive Poisson ratios, the larger lateral deflections have occurred in regions that are closer to the smaller edge. On the other hand, as Fig. 15 that has been plotted for the S-C conical shell implies, due to the difference in direction of the bending moment caused by the shear traction at x=0.3 and x=0.8, direction of the curvature of the shell is different in the mentioned sections. Fig. 16 reveals that the resulting curvatures become remarkably smaller at ν = - 0.9, as the previous cases.

Figure 14:
Effects of the auxeticity of the material on meridian distribution of the lateral deflection of a C-C truncated conical shell subjected to a shear traction.

Figure 15:
Effects of the auxeticity of the material on meridian distribution of the lateral deflection of a S-C truncated conical shell subjected to a shear traction.

Figure 16:
Effects of the auxeticity of the material on meridian distribution of the lateral deflection of a F-C truncated conical shell subjected to a shear traction.

Distributions of the lateral deflection and in-plane stress of the truncated conical shell have to resemble distributions of it special configurations, i.e., the cylindrical shell and annular plate. Comparing of Figs. 14 to 16 and the transverse distributions of the in-plane stress shown in Figs. 17 to 19 for the three types of the boundary conditions with those obtained for the annular plate and cylindrical shell, reveals that trends of present results are more similar to those of the annular plate (whose rigidity is different at its two edges, as the present case). Curvature of the stress distribution curves is due to transverse variations of the material properties. Moreover, since in addition to the bending moment, the shear traction exerts compressive stresses on one of the considered sections and tensile stresses on the other, location of the neutral layer is not identical for the two sections, as may be inferred from Figs. 17 to 19. Since in contrast to the conventional materials, the auxeticity causes the material to shrink in all directions in compressive stresses and expand in all directions in tensile stresses (so that the cross section increases remarkably), the stress field varies in a complicated manner; so that the stresses have become less in the considered sections of the C-C conical shell in more negative Poisson ratios but increased significantly in the neighborhood of the simply supported and free edges of the S-C and F-C shells, respectively.

Figure 17:
Effects of the auxeticity of the material on transverse distribution of the in-plane stress of a C-C truncated conical shell subjected to a shear traction.

Figure 18:
Effects of the auxeticity of the material on transverse distribution of the in-plane stress of a S-C truncated conical shell subjected to a shear traction.

Figure 19:
Effects of the auxeticity of the material on transverse distribution of the in-plane stress of a F-C truncated conical shell subjected to a shear traction.

4.2 Uniform Thickness FGM Shells/Plates Under Combined Shear and Normal Tractions

In the present section, responses of functionally graded auxecit shells/plates subjected to simultaneous normal and shear tractions are studied. Each of these distinct tractions may have different interactions with the negative Poisson ratios. As the previous section, h/L ratio is equal to 0.1. To present a more useful discussion, transverse distribution of the in-plane displacement component of the plate is also investigated for the x/L=0.3 and 0.8 sections.

Radial distributions of the lateral deflection of the annular plate are shown in Figs. 20 to 22, for the three types of the boundary conditions, respectively. As before, trends of variations of the distributions of the lateral deflection are different for the ν = - 0.9, ..., 0 and ν = 0, ..., 0.3 intervals of the Poisson ratio. However, the trends are more robust, since the transverse load has the dominant effects. Results related to the lower and higher limits of the Poisson ratio are also verified by the 3D elasticity results extracted from ABAQUS.

Figure 20:
Effects of the auxeticity of the material on radial distribution of lateral deflection of a C-C annular plate subjected to combined shear and normal tractions.

Figure 21:
Effects of the auxeticity of the material on radial distribution of lateral deflection of a S-C annular plate subjected to combined shear and normal tractions.

Figure 22:
Effects of the auxeticity of the material on radial distribution of lateral deflection of a F-C annular plate subjected to combined shear and normal tractions.

The transverse distribution of the in-plane stress of the C-C plate plotted in Fig. 23, show that location of the neutral layer is affected not only by the shear traction but also especially by the transverse load. Moreover, gradients of the distribution are higher in the neighborhood of the internal edge.

Figure 23:
Effects of the auxeticity of the material on transverse distribution of in-plane stress of a C-C annular plate subjected to combined shear and normal tractions.

Transverse distribution of the in-plane displacement component is depicted in Fig. 24. From this figure, one may deduce that both the extensional and bending rigidities increase drastically by employing higher negative Poisson ratios.

Figure 24:
Effects of the auxeticity of the material on transverse distribution of in-plane displacement of a C-C annular plate subjected to combined shear and normal tractions.

Axial distributions of the lateral deflection of the cylindrical shell for various Poisson ratios are illustrated in Figs. 25 to 27, for the three types of the boundary condition. As may be noted, location of the maximum lateral deflection shifts to the initial edge, as auxeticity of the shell increases. Moreover, in this case (as the proceeding section), in contrast to the annular plate, a noticeable maximum is achieved in the F-C shell. However, since the auxeticity increases rigidity of the shell and extends influence of the edge condition on farther interior regions, this apparent maximum has vanished for higher negative Poisson ratios.

Figure 25:
Axial distribution of lateral deflection of a C-C cylindrical shell subjected to combined shear and normal tractions, for various Poisson ratios.

Figure 26:
Axial distribution of lateral deflection of a S-C cylindrical shell subjected to combined shear and normal tractions, for various Poisson ratios.

Figure 27:
Axial distribution of lateral deflection of a F-C cylindrical shell subjected to combined shear and normal tractions, for various Poisson ratios.

Fig. 28 reveals that in this case, transverse distribution of the in-plane stress manifests severe gradients in the neighborhood of the second edge of the shell that exaggerates at larger negative Poisson ratios. As the previous case, Fig. 29 that presents effects of the Poisson ration on transverse distribution of the in-plane displacement of the C-C shell, shows the significant growth in the bending and extensive rigidities of the sections by increasing the auxeticity of the material mixture.

Figure 28:
Transverse distribution of in-plane stress of a C-C cylindrical shell subjected to combined shear and normal tractions, for various Poisson ratios.

Figure 29:
Transverse distribution of in-plane displacement component of a C-C cylindrical shell subjected to combined shear and normal tractions, for various Poisson ratios.

Influence of the negative Poisson ratio on the meridian distribution of the lateral deflection for the three types of the edge conditions is illustrated in Figs. 30 to 32, and effects of the auxeticity on transverse distributions of the in-plane stress and displacement components of the C-C truncated conical shell are plotted in Figs. 33 and 34, respectively. Comparing these results with those of the annular plate and the cylindrical shell reveals that responses of the truncated conical shell resemble those of the annular plate and the cylindrical shell in some aspects; so that a combination of characteristics of the mentioned structures are reproduced here. However, the in-plane stresses are higher than those of the cylindrical shell. Moreover, in the present situations, the lateral deflection of free edge of the F-C shell becomes larger when more negative Poisson ratios are employed (Fig. 32). Indeed, in highly auxetic materials, the work consumed for establishment of the intermediate curvatures in the shell is used for deformation of the free edge, due to the higher bending rigidity of the shell.

Figure 30:
Meridian distribution of lateral deflection of a C-C truncated conical shell subjected to combined shear and normal tractions, for various Poisson ratios.

Figure 31:
Meridian distribution of lateral deflection of a S-C truncated conical shell subjected to combined shear and normal tractions, for various Poisson ratios.

Figure 32:
Meridian distribution of lateral deflection of a F-C truncated conical shell subjected to combined shear and normal tractions, for various Poisson ratios.

Figure 33:
Transverse distribution of in-plane stress of a C-C truncated conical shell subjected to combined shear and normal tractions, for various Poisson ratios.

Figure 34:
Transverse distribution of in-plane displacement of a C-C truncated conical shell subjected to combined shear and normal tractions, for various Poisson ratios.

4.3 Plates and Shells Under Non-Uniform Shear and Normal Tractions

This section is included merely to introduce some of capabilities of the proposed formulation. In this regard, a C-C type shell/plate of thickness ratio h/L=0.2 whose bottom section is subjected to normal and shear traction whose intensity varies according to q = x + x ²) is adopted. As mentioned before, this case may resemble some practical cases, e.g. a cases wherein passage of compressible/incompressible fluids imposes pressures and frictional forces on the internal surfaces of a pipe or nozzle. Since the overall stress variations behavior is similar to those of the previous section, they are not reported here to avoid repetitive results. However, meridian/radial distribution of the lateral deflection of the annular plate as well as the conical and cylindrical shells are illustrated simultaneously in Fig. 35, for various values of the Poisson ratio. Since the higher tractions are imposed in the second half of the meridian coordinate region, location of the maximum lateral deflection has shifted toward that half, in comparison to the previous cases. = (1 +

Figure 35:
Meridian/radial distribution of the lateral deflection of the C-C annular plate and conical and cylindrical shells under variable normal and shear tractions, for various Poisson ratios.

4.4 Plates and Shells with Variable Thickness and Various Boundary Conditions

Finally, simultaneous effects of the auxeticity of the material, thickness variability, and boundary conditions on distributions of the stresses and lateral deflections of the plate and shells are investigated. The shell/plate is subjected to a normal traction whose intensity is q = 1MPa.

First, C-C type plates and shells with a linear variation of the thickness, _^hb = -0.05, _^ht = 0.1x are considered. Meridian/radial distributions of the lateral deflection of the plate and shells are plotted in Fig. 36, for various Poisson ratios. Transverse distribution of the in-plane stress of the annular plate and cylindrical and conical shells are depicted in Figs. 37 to 39, for three sections (x/L=0.3, 0.5, and 0.8) and ABAQUS verification of some of the results is included as well, for various Poisson ratios. As before, lateral deflection of the cylindrical shell is larger; since for an identical _^Rb radius, the cylindrical configuration is more impact (radii of the truncated conical shell are generally larger than _^Rb ) and therefore, its rigidity is higher. Furthermore, since the thickness is greater at the second/exterior edge of the plate/shell, the maximum lateral deflection and the maximum meridian stress have occurred in regions with more flexibility, i.e., those located in the neighborhood of the first/interior edge. Furthermore, due to the occurrence of inflection point whose location may be explored based on Fig. 36, bending direction of the cross section and consequently, sign of the stresses of the plate/shell has changed for the x/L=0.8 section.

Figure 36:
Simultaneous effects of the auxeticity and linear thickness variability on meridian/radial distribution of the lateral deflection of the considered plates and shells.

Figure 37:
Simultaneous effects of the auxeticity and linear thickness variability on transverse distribution of the radial stress of the C-C annular plate.

Figure 38:
Simultaneous effects of the auxeticity and linear thickness variability on transverse distribution of the axial stress of the C-C cylindrical shell.

Figure 39:
Simultaneous effects of the auxeticity and linear thickness variability on transverse distribution of the meridian stress of the C-C truncated conical shell.

In the next stage, a parabolic thickness variation (_^hb = -0.05, _^ht = 0.1x (1+ x)) is adopted and simultaneous effects of the auxeticity and thickness variability are evaluated for the C-C and S-C boundary conditions. The resulting distributions of the lateral deflection are plotted in Figs. 40 and 41, respectively. The previous conclusions are confirmed here once again; since the maximum deflections have occurred in sections that are closer to the first/interior edge of the structure and rigidity of the plate/shell has increased drastically with the auxeticity.

Figure 40:
Simultaneous effects of the auxeticity and parabolic thickness variability on meridian/radial distributions of the lateral deflection of the C-C plates and shells.

Figure 41:
Simultaneous effects of the auxeticity and parabolic thickness variability on meridian/radial distributions of the lateral deflection of the S-C plates and shells.

5 CONCLUSIONS

In the present research, a unified Taylor transform solution is presented for stress and displacement analysis of the variable thickness auxetic functionally graded conical and cylindrical shells and circular/annular plates. The top or bottom layers of the shell/plate may be subjected to general non-uniform normal and shear tractions and the edges may be simply supported, clamped or free.

Some of the obtained practical results are:

• Combination of the auxeticity of the functionally graded material, structure configuration, shear/normal tractions, thickness variability, and edge conditions leads to complicated and sometimes non-predictable behaviors.

• The differences between results of the successive values of the Poisson ratio grow as the Poisson ratio becomes more negative.

• Trends of variations of the results are different for the ν = -0.9, ..., 0 and ν = 0, ..., 0.3 ranges of the Poisson ratio.

• Stress distributions of the C-C and F-C plates exhibit respectively, the most and the least sensitivity to the auxeticity of the material.

• There are similarities between distributions of the lateral deflection and in-plane stress of the truncated conical shells, cylindrical shells, and annular plates, for the three types of the boundary conditions.

• Results reveal that for the considered material properties, the overall lateral rigidity of the plate has grown to up to six to seven times, for ν = -0.9, for all the C-C, S-C, and F-C boundary conditions.

• Since the lateral stiffness is proportional to the inverse of (1 - ν ²), the worst lateral stiffness of the annular plate and conical shell belongs to ν = 0. This hint is confirmed by the results as well.

• In contrast to the lateral deflection results, more regular trends may be noted for the stress results, with varying the Poisson ratio.

• Effects of the auxeticity on the in-plane stresses are more remarkable at sections with more flexibility (in the vicinity of the free or simply supported edges) and regions located at the neighborhood of the first edge.

• For positive Poisson ratios, the larger lateral deflections have occurred in regions that are closer to the smaller edge.

• The maximum deflections have occurred in sections that are closer to the first/interior edge of the structure. However, auxeticity of the material affects the resulting curvatures and location of the maximum lateral deflection.

References

• Abediokhchi, J., Shakouri, M., Kouchakzadeh, M.A. (2013). Bending analysis of moderately thick functionally graded conical panels with various boundary conditions using GDQ method, Composite Structures 103: 68-74.
• Aghdam, M.M., Shahmansouri, N., Bigdeli, K. (2011). Bending analysis of moderately thick functionally graded conical panels, Composite Structures 93: 1376-84.
• Alipour, M.M., Shariyat, M. (2011). Semi-analytical buckling analysis of heterogeneous variable thickness viscoelastic circular plates on elastic foundations, Mechanics Research Communications 38: 594-601.
• Alipour, M.M., Shariyat, M. (2012). An elasticity-equilibrium-based zigzag theory for axisymmetric bending and stress analysis of the functionally graded circular sandwich plates, using a Maclaurin-type series solution, European Journal of Mechanics A/Solids 34: 78-101.
• Alipour, M.M., Shariyat, M. (2013). Semianalytical solution for buckling analysis of variable thickness two-directional functionally graded circular plates with nonuniform elastic foundations, ASCE Journal of Engineering Mechanics 139: 664-676.
• Alipour, M.M., Shariyat, M. (2015). Analytical zigzag formulation with 3D elasticity corrections for bending and stress analysis of circular/annular composite sandwich plates with auxetic cores, Composite Structures 132: 175-197.
• Asemi, K., Salehi, M., Akhlaghi, M. (2011). Elastic solution of a two dimensional functionally graded thick truncated cone with finite length under hydrostatic combined loads, Acta Meccanica 217: 119-34.
• Behravan Rad, A., Shariyat, M. (2013). A three-dimensional elasticity solution for two-directional FGM annular plates with non-uniform elastic foundations subjected to normal and shear tractions, Acta Mechanica Solida Sinica 26: 671-690.
• Behravan Rad, A., Shariyat, M. (2015). Three-dimensional magneto-elastic analysis of asymmetric variable thickness porous FGM circular plates with non-uniform tractions and Kerr elastic foundations, Composite Structures 125: 558-574.
• Brighenti, R. (2014). Smart behavior of layered plates through the use of auxetic materials, Thin-Walled Structures 84: 432-42.
• Carrera, E., Brischetto, S., Nali, P. (2011), Plates and Shells for Smart Structures: Classical and Advanced Theories for Modeling and Analysis, John Wiley & Sons.
• Carrera, E., Brischetto, S. (2008). A Survey with numerical assessment of classical and refined theories for the analysis of sandwich plates, Applied Mechanics Review 62: 010803.
• Duc, N.D. (2013). Nonlinear dynamic response of imperfect eccentrically stiffened FGM double curved shallow shells on elastic foundation, Composite Structures 99: 88-96.
• Duc, N.D., Cong, P.H. (2016). Nonlinear dynamic response and vibration of sandwich composite plates with negative Poisson’s ratio in auxetic honeycombs, Journal of Sandwich Structures and Materials, No. Cic: 1-26, DOI: 10.1177/1099636216674729.
  » https://doi.org/10.1177/1099636216674729
• Ferreira, A.J.M., Batra, R.C., Roque, C.M.C., Qian, L.F., Martins, P.A.L.S. (2005). Static analysis of functionally graded plates using third-order shear deformation theory and a meshless method, Composite Structures 69: 449-57.
• Ghannad, M., Rahimi, G.H., Zamani Nejad, M. (2013). Elastic analysis of pressurized thick cylindrical shells with variable thickness made of functionally graded materials, Compos. Part B 45: 388-396.
• Golmakani, M.E., Kadkhodayan, M. (2014). An investigation into the thermoelastic analysis of circular and annular functionally graded material plates, Mechanics of Advanced Materials and Structures 21: 1-13.
• Hilton, H.H., Lee, D.H., El Fouly A.R.A. (2008). Generalized viscoelastic designer functionally graded auxetic materials engineered/tailored for specific task performances, Mechanics of Time-Dependent Materials 12, No. 2: 151-178.
• Imbalzano, G., Tran, P., Ngo, T.D., Lee, P.V.S. (2015). Three-dimensional modelling of auxetic sandwich panels for localised impact resistance, Journal of Sandwich Structures and Materials, DOI: 10.1177/1099636215618539.
  » https://doi.org/10.1177/1099636215618539
• Khan, K.A., Hilton, H.H. (2009). On Inconstant Poisson's ratios in non-homogeneous elastic media, Journal of Thermal Stresses 33, No. 1: 29-36.
• Matsunaga, H. (2009). Free vibration and stability of functionally graded circular cylindrical shells according to a 2D higher-order deformation theory, Composite Structures 88: 519-531.
• Mercan, K., Demir, Ç., Civalek, Ö. (2016). Vibration analysis of FG cylindrical shells with power-law index using discrete singular convolution technique, Curved and Layered Structures 3, No. 1: 82-90.
• Nie, G., Zhong, Z. (2007). Axisymmetric Bending of Two-Directional Functionally Graded Circular and Annular Plates, Acta Mechanica Solida Sinica 20, No. 4: 289-295.
• Nosier, A., Fallah, F. (2009). Non-linear analysis of functionally graded circular plates under asymmetric transverse loading, International Journal of Non-Linear Mechanics 44: 928 - 942.
• 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: 1131-1158.
• Reddy, J.N. (2003), Mechanics of laminated composite plates and shells: Theory and analysis, 2nd ed. CRC Press (Florida).
• Sburlati, R., Bardella, L. (2011). Three-dimensional elastic solutions for functionally graded circular plates, European Journal of Mechanics A/Solids 30: 219-235.
• Selahi, E., Setoodeh, A.R., Tahani, M. (2014). Three-dimensional transient analysis of functionally graded truncated conical shells with variable thickness subjected to an asymmetric dynamic pressure, International Journal of Pressure Vessels and Piping 119: 29-38.
• Setoodeh, A.R., Tahani, M., Selahi, E. (2011). Hybrid layerwise-differential quadrature transient dynamic analysis of functionally graded axisymmetric cylindrical shells subjected to dynamic pressure, Composite Structures 93: 2663-2670.
• Setoodeh, A.R., Tahani, M., Selahi, E. (2012). Transient dynamic and free vibration analysis of functionally graded truncated conical shells with non-uniform thickness subjected to mechanical shock loading, Composites Part B 43: 2161-2171.
• Shariyat, M. (2008). Dynamic buckling of suddenly loaded imperfect hybrid FGM cylindrical shells with temperature-dependent material properties under thermo-electro-mechanical loads, International Journal of Mechanical Sciences 50: 1561-1571.
• Shariyat, M. (2010a). A generalized high-order global-local plate theory for nonlinear bending and buckling analyses of imperfect sandwich plates subjected to thermo-mechanical loads, Composite Structures 92: 130-143.
• Shariyat, M. (2010b). A generalized global-local high-order theory for bending and vibration analyses of sandwich plates subjected to thermo-mechanical loads, International Journal of Mechanical Sciences 52: 495-514.
• Shariyat, M. (2012). Nonlinear transient stress and wave propagation analyses of the FGM thick cylinders, employing a unified generalized thermoelasticity theory, International Journal of Mechanical Sciences 65: 24-37.
• Shariyat, M., Alipour, M.M. (2011). Differential transform vibration and modal stress analyses of circular plates made of two-directional functionally graded materials resting on elastic foundations, Archive of Applied Mechanics 81: 1289-1306.
• Shariyat, M., Alipour, M.M. (2014). A novel shear correction factor for stress and modal analyses of annular FGM plates with non-uniform inclined tractions and non-uniform elastic foundations, International Journal of Mechanical Sciences 87: 60-71.
• Shariyat, M., Alipour, M.M. (2015). Novel layerwise shear correction factors for zigzag theories of circular sandwich plates with functionally graded layers, Latin American Journal of Solids and Structures, 12: 1362-1396.
• Shariyat, M., Asemi, K. (2014). 3D B-spline finite element nonlinear elasticity buckling analysis of rectangular FGM plates under non-uniform edge loads, using a micromechanical model, Composite Structures 112: 397-408.
• Shariyat, M., Asgari, D. (2013). Nonlinear thermal buckling and postbuckling analyses of imperfect variable thickness temperature-dependent bidirectional functionally graded cylindrical shells, International Journal of Pressure Vessels and Piping 111-112: 310-320.
• Shariyat, M., Mohammadjani, R. (2013). Three-dimensional compatible finite element stress analysis of spinning two-directional FGM annular plates and disks with load and elastic foundation non-uniformities, Latin American Journal of Solids and Structures 10: 859-890.
• Shariyat, M., Mohammadjani, R. (2014). Three-dimensional stress field analysis of rotating thick bidirectional functionally graded axisymmetric annular plates with nonuniform loads and elastic foundations, Journal of Composite Materials 48: 2879-2904.
• Shen, H.-S. (2009), Functionally Graded Materials: Nonlinear Analysis of Plates and Shells, Taylor & Francis Group, LLC.
• Sofiyev, A.H. (2007). Thermoelastic stability of functionally graded truncated conical shells, Composite Structures 77: 56-65.
• Sofiyev, A.H. (2014). On the dynamic buckling of truncated conical shells with functionally graded coatings subject to a time dependent axial load in the large deformation, Composites Part B 58: 524-533.
• Sun, J., Xu, X., Lim, C.W. (2014). Buckling of functionally graded cylindrical shells under combined thermal and compressive loads, J. Therm. Stresses 37: 340-362.
• Tahani, M., Setoodeh, A.R., Selahi, E (2013). Three-dimensional transient analysis of functionally graded cylindrical shells subjected to asymmetric dynamic pressure, Science and Engineering of Composite Materials 20: 75-85.
• Tornabene, F. (2009). Free vibration analysis of functionally graded conical, cylindrical shell and annular plate structures with a four-parameter power-law distribution, Computer Methods in Applied Mechanics and Engineering 198: 2911-2935.
• Tornabene, F., Fantuzzi, N. (2014), Mechanics of Laminated Composite Doubly-Curved Shell Structures, Società Editrice Esculapio (Milan).
• Tornabene, F., Fantuzzi, N., Viola, E., Batra, R.C. (2015). Stress and strain recovery for functionally graded free-form and doubly-curved sandwich shells using higher-order equivalent single layer theory, Composite Structures 119: 67-89.
• Tornabene, F., Viola, E. (2013). Static analysis of functionally graded doubly-curved shells and panels of revolution, Meccanica 48: 901-930.
• Viola, E., Rossetti, L., Fantuzzi, N. (2012). Numerical investigation of functionally graded cylindrical shells and panels using the generalized unconstrained third order theory coupled with the stress recovery, Composite Structures 94: 3736-58.
• Viola, E., Rossetti, L., Fantuzzi, N., Tornabene, F. (2014). Static analysis of functionally graded conical shells and panels using the generalized unconstrained third order theory coupled with the stress recovery, Composite Structures 112: 44-65.
• Wang, Y., Xu, R., Ding, H (2010). Three-dimensional solution of axisymmetric bending of functionally graded circular plates, Composite Structures 92: 1683-1693.
• Wu, Z., Chen, W., Ren, X. (2009). Refined global-local higher-order theory for angle-ply laminated plates under thermo-mechanical loads and finite element model, Composite Structures 88: 643-658.
• Zamani Nejad, M., Jabbari, M., Ghannad, M. (2015). Elastic analysis of FGM rotating thick truncated conical shells with axially-varying properties under non-uniform pressure loading, Composite Structures 122: 561-569.
• Zozulya, V.V., Zhang, Ch. (2012). A high order theory for functionally graded axisymmetric cylindrical shells, International Journal of Mechanical Sciences 60: 12-22.

Publication Dates

• Publication in this collection
  June 2017

History

• Received
  04 Oct 2016
• Reviewed
  07 Mar 2017
• Accepted
  13 Mar 2017