Static Analysis of Cylindrical Shells Akshaya Dhananjay Patil 1, Mayuri Madhukar Jadhav 2 1,2 Assistant Professor, Dept. of Civil Engineering, Padmabhooshan Vasantraodada Patil Institute Of Technology, Shivaji University, Kolhapur, India. ---------------------------------------------------------------------***--------------------------------------------------------------------- Abstract - In this paper, an isotropic cylindrical shell is studied under the external pressure. For and uniformly distributed load and sinusoidal load at simply supported end conditions, cylindrical shell is studied for its deflection and the von-mises stresses are analyzed. Analytical modeling is based on first order shear deformation theory (FOST) and a finite element computational tool ABAQUS [2] is used to model the isotropic cylindrical shell. Key Words: Isotropic Cylindrical shell,fost,abaqus, von- Mises stresses. 1.INTRODUCTION Shells are commonly used in many engineering structures, such as aerospace, automotive and submarine structures. Isometric shells are getting more attention recently. Because of the simplicity of shell theories, it is favorable to analyze shell structures with shell theories instead of three dimensional (3D) theories of elasticity. Shell theories redeem the difficulty of shell analyses by employing certain assumptions on the behavior of displacements along the thickness direction. For instance, First order expansion of inplane displacements leads to First order shear deformation theories (FSDTs). A difficulty rise and appears in both in strain displacements and stress resultants in the equation of derivation of basic equation of shell. The term was neglected by first analyst of composite thin shells which is understandable for thin shells.although the importance of the inclusion of the term has been tested for the thicker shells and proves to be essential. 2. Formulation of FSDT 2.1 Introduction The FSDT developed by Mindlin [12] accounts for the shear deformation effect by the way of a linear variation of the inplane displacements through the thickness. It is noted that the theory developed by Reissner [13] also accounts for the shear deformation effect. However, the Reissner theory is not similar with the Mindlin theory like erroneous perception of many researchers through the use of misleading descriptions such as Reissner-Mindlin plates and FSDT of Reissner.The major difference between two theories was established by Wang et al. by derivating the bending relationships between Mindlin and Reissner quantities for a general plate problems. Since the Reissner theory was based on the assumption of a linear bending stress distribution, its formulation will inevitably lead to the displacement variation being not necessarily linear across the plate thickness. Thus, it is incorrect to refer to the Reissner theory as the FSDT which implies a linear variation of the displacements through the thickness. Another difference between two theories is that the normal stress which was included in the Reissner theory was omitted in the Mindlin. The FSDT was used to model FG shells. Reddy and Chin [14] studied the dynamic response of FG cylinders and plates subjected to two different types of thermal loadings using the FSDT and the finite element method. 2.2 Definition of Displacement field 2017, IRJET Impact Factor value: 6.171 ISO 9001:2008 Certified Journal Page 618 (1a) In the above relations, the terms and are the displacements of a general point ( ) in the domain in x, y and z directions respectively. The parameters the inplane displacements and are is the transverse displacement of a point (x, y) on the element middle plane. The functions are the rotations of the normal to the element middle plane about y- and x-axis respectively. 2.3 Strain displacement relation With the definition of strains from the linear theory of elasticity, assuming, the general straindisplacement relations in the curvilinear co-ordinate system are given as follows:
accomplished through the transformations.the final relations are as follows: Substituting the expressions for displacements at any point within the space given by Eqs. (1a) for the displacement models considered herein, the linear strains in terms of middle surface displacements, for each displacement model can be obtained as follows: (2b) These equations in compacted form may be written as Where, 2.4 Stress-strain relations and stress resultants Assuming the principal material axes (1,2,3) and the shell axes (x,y,z) in the curvilinear co-ordinate system,the three dimensional stress-strain relations for an cylindrical shell with reference to the principal material axes for the theory to be developed based on the displacement are defined as follows: in which the coefficients of the Q matrix, called as reduced elastic constants are defined in Appendix A. The components of the strain vector and the corresponding components of the stress- resultant vector are defined as follows: (2a) Where, These equations in compacted form may be written as As mentioned earlier, the relations given by Eq. (2a) are the stress-strain constitutive relations for the cylindrical shell referred to shell s principal material axes (1,2,3). The principal material axes of shell may not coincide with the reference axes of the shell (x,y,z).it is therefore necessary to transform the constitutive relations from the shell s material axes (1,2,3) to reference axes (x,y,z).this is conveniently 2.5 Equilibrium equations For equilibrium equations, the total potential energy must be stationary and using the definitions of stress-resultants and mid-surface strains stated in above sections principal of virtual work yields Where, U is the strain energy and W represents the work done by external forces. These are evaluated as follows: 2017, IRJET Impact Factor value: 6.171 ISO 9001:2008 Certified Journal Page 619
, in which a and b are the dimensions of shell middle surface along the x and y axes respectively. (3a) Integration through thickness and by substituting in terms of strains and introducing stress resultants, the above relations transform in the following form The exact form of spatial variation of mid-surface displacements is given by (4a) (3b) In the above equation is the distributed transverse load. The governing equations of equilibrium can be derived from eq. (3b) by integrating the displacement gradients in mid surface strains by parts and setting the coefficients of derivatives of mid-surface displacements to zero separately. Thus one obtains the following equilibrium equations. 3. Analysis in ABAQUS To find the dimensionless displacements of an isotropic cylindrical shells having following properties. a/b= 1,, υ=0.25 1. Find displacement for a/h=20 and for different ratios of a) a/r= 0.5 b) a/r=1 c) a/r=2 (3c) 2. Find the displacement for a/h=10, and for different ratios of a) a/r= 0.5 b) a/r=1 c) a/r=2 In addition following line integrals are also obtained (3d) 2.6 Closed form solutions Sinusoidal variation of transverse load is considered as under: Fig-1 Modeling in ABAQUS 2017, IRJET Impact Factor value: 6.171 ISO 9001:2008 Certified Journal Page 620
Transverse Displacement (w) Transverse Displacement (w) International Research Journal of Engineering and Technology (IRJET) e-issn: 2395-0056 Table 2. Non-Dimensional transverse deflection for longitudinal edges simply supported cylindrical shell under uniformly distributed loading with a/h = 10. a/r 3D FSDTQ FSDT ABAQUS (present) 0.5 40.875 40.956 40.699 41.1 1 28.415 28.333 28.009 29 2 12.242 12.108 11.972 11.52 Fig-2 Shows stresses and deflection [ ] % error w.r.t. FSDT results Table 1. Non-Dimensional transverse deflection for longitudinal edges simply supported cylindrical shell under uniformly distributed loading with a/h = 20. [ ] % error w.r.t. FSDT results a/r 3D FSDTQ FSDT ABAQUS (present) 0.5 26.698 26.679 26.596 25.87 1 11.427 11.388 11.337 11.19 2 1.6384 3.0070 2.9994 2.92 45 40 35 30 25 20 15 10 5 FSDT ABAQUS result 30 25 20 15 10 5 0 Fig.3. Non- dimensional transverse displacement vs. FSDT results having a/h = 20. 0.5 1 2 a/r FSDT 0 ABAQUS result 0.5 1 2 a/r Fig.4. Non- dimensional transverse displacement vs. FSDT results having a/h = 10. 3. CONCLUSIONS It is concluded that for a simply supported isotropic cylindrical shell subjected to uniformly distributed load, for a constant ratio of a/central displacement decreases as a/r ratio increases. 4. REFERENCES 1. Schurg M,Wagner W.,Gruttmann F., An enhanced FSDT model for the calculation of interlaminar shear stresses in composite plate structures, Compt Mech (2009)44:765.doi:1007/s00466-009-0410-7 2. ABAQUS User s Manual (2003). ABAQUS Version.6.10 2017, IRJET Impact Factor value: 6.171 ISO 9001:2008 Certified Journal Page 621
3. A.S.Oktem, J.L. Mantari, C. Guedes Soares, Static response of functionally graded plates and doublycurved shells based on a higher order shear deformation theory,european Journal of Mechanics A/Solids36(2012)163-172 4. Erasmo Viola, Francesco Tornabene, Nicholas Fantuzzi, Static analysis of completely doubly-curved laminated shells and panels using general higher-order shear deformation theories, Composite Structures101(2013)59-93 5. Dao Huy Bich, Dao Van Dung, Vu Hoai Nam, Nguyen ThiPhuong, Nonlinear static and dynamic buckling analysis of imperfect eccentrically stiffened functionally graded circular cylindrical thin shells under axial compression,international Journal of Mechanical Sciences 74(2013)190-200 12. R. D. Mindlin, Influence of rotatory inertia and shear on flexural motions of istotropic elastic plates, Journal ofapplied Mechanics. (1951) 31 38 13. E. Reissner, The effect of transverse shear deformation on bending of elastic plates, Journal of Applied mechanics. 12 (1945) 69 77. 14. Reddy,Chin, Thermomechanical analysis of functionally graded cylinders and plates, Journal of Thermal stresses 21(6),593-626 15. Timoshenko SP, Woinowsky-Krieger S. Theory of plates and shells. McGraw-Hill; 1959 6. J.L. Mantari, E.V.Granados, A refined FSDT for the static analysis of functionally graded sandwich plates,thin- Walled Structures90(2015)150-158 S.M.Shiyekar, Prashant Lavate, Flexure of power law governed functionally graded plates using ABAQUS UMAT,Aerospace Science and Technology 46(2015)51-59 7. Ankit Gupta,Mohammad Talha, Recent development in modeling and analysis of functionally graded materials and structures,progress in Aerospace Sciences M. Memar Ardestani, B.Soltani, Sh. Shams, Analysis of functionally graded stiffened plates based on FSDT utilizing reproducing kernel particle method,composite Structures112(2014)231-240 8. Erasmo viola, Luigi Rossetti, Nicholas Fantuzzi, Francesco Tornabene, Static analysis of functionally graded conical shells and panels using the generalized unconstrained third order theory coupled with the stress recovery,composite Structures 112 (2014) 44-65 9. Dao Huy Bich, Dao Van Dung, Vu Hoai Nam, Nonlineardynamical analysis of eccentrically stiffened functionally graded cylindrical panels,journal of CompositeStructure 94(2012)2465-2473 10. Erasmo Viola, Luigi Rossetti, Nicholas Fantuzzi, Numerical investigation of functionally graded cylindrical shells and panels using the generalized unconstrained third order theory coupled with the stress recovery,composite Structures 94 (2012) 3736-3758 11. Ebrahim Asadi, Wencho Wang, Mohamad S. Qatu, Static and vibration analyses of thick deep laminated cylindrical shells using 3D and various shear deformation theories,composite Structures 94 (2012) 494-500 2017, IRJET Impact Factor value: 6.171 ISO 9001:2008 Certified Journal Page 622