For the purpose of demonstrating the applicability of the previously derived theories, the problem of a laminated circular cylindrical shell under internal pressure and edge loadings will be examined. The cylinder is assumed to consist of boron/epoxy composite layers. Each layer is taken to be homogeneous but anisotropic with an arbitrary orientation of th e elastic axes. We need not consider the restriction of the symmetry of the layering due to the non -homogeneity considered in the original development of the theory expressed by the constitutive equations. Thus, each layer can possess a different thickne ss. Keywords: Shell Theory, circular cylindrical shell, applicability INTRODUCTION We assume here that the contact between layers is such that the strains are continuous functions of the thickness coordinate. As the C’s are piecewise continuous functions, the in -plane stresses are also. We would expect them to be discontinuous at the juncture of layers of dissimilar materials. The transverse stresses are continuous functions of the thick ness coordinate. Although as mentioned above the theory developed can consider, random layering, numerical results are to be carried out for a four layer symme tric angle ply configuration (see Appendix C for asymmetric cross -ply lay ering). For this configuration the angle of elastic axes y is oriented at +γ , -γ, -γ, +γ , with the shell axis and the layers are of equal thickness. The coefficients Aij Bij Dij Eij and Fij for this configuration can be determined by making use of the results for layered shells developed in Appendix B. Let the cylinder be subjected to an internal pressure p, an axial force per unit circumferential length N and a torque T. The axial force is taken to be applied at r=a+H such that a moment N(H -d) is produced about the reference surface r=a+d. We introduce dimensionless external force and moments as follows: Equation (1) Asian Institute of Research Engineering and Technology Quarterly Reviews Vol.4, No. 1, 202 1 11 The cylinder is taken to be clamped at both ends but free to rotate and extend axially at one end. The edge conditions can thus be written as Equation (2) Here, l is the dimensionless length of the cylinder. The above edge conditions are assumed to represent a close approximation to the test set -up used for obtaining the mechanical properties of the composites. In the theories developed in the previous chapters, the distance d at which the stress resultants were defined was left arbitrary. We now choose it to be such that there exists no coupling between N z and K 1 and Mz and ε1d This can be achieved by setting the first component of sub - matrix [ B ] equal to zero, B ij = 0 This yield Equation (3) For a homogeneous, isotropic material we obtain, for example, As the loading applied at the end of the shell is axi -symmetric, all the stresses and strains are also taken to be axi -symmetric. We thus can set all the φ derivatives in the expressions for the stresses and strains and in the equations for the displaceme nts equal to zero. The resulting simplifications and the general solution corresponding to each of the der ived theories is now given. As the theory corresponding to axial length scale a(a/h)1/2 and circumferential length scale a has been used in the analysis of long cylindrical shells with edges of the form φ = constant , it will not be considered. a) Theory associated with Length Scales a Strain -Displacement Relations Asian Institute of Research Engineering and Technology Quarterly Reviews Vol.4, No. 1, 202 1 12 Governing Equations Equation (4) By com bi ni ng th e fi r st t wo equ at i on s, we o bt ain ex pr essio n s rel at in g Vz and Vr and V θ and Vr , respectively, Equation (5) where the coefficients R and £ are defined as follows: Equation (6) Integrating equations with respect to x yields where C 1* and C 2* are constants of integrations. On substituting the equations into the last equation of (4) we obtain, Wh er e Equation (1 0) sh ows that the r adial deform ation dep end s on t he pr essure P* and the material properties of the individual layers. These effective properties i n turn dep end o n th e an gle y bet ween th e elastic axi s of t he material of each layer and the shell axis. a : Inside Radius of Cylindrical Shell h : Total Thickness of the Shell Wall d; Distance (thickness) from Inside Radius and to mechanical neutral surface Si : Radius of Each Layer of Wall (I = 1, 2, 3 -- - to the number of layer) L : Longitudinal Length Scale to be defined, Also Actual Length of the Cyl indrical Shell Π : Circumferential Length Scale of cylindrical shell to be defined Ei : Young’s Moduli in I Direction Gij : Shear Moduli in i -j Face Sij : Compliance Matrix of Materials of Each Layer r : Radial Coordinate Π : Circumferential Length Scale to be defined Ɣ : Angle of Fiber Orientation σ : Normal Stresses ε : Normal Strains z, θ, r : Generalized Coordinates in Longitudinal, Circumferential and Radial Asian Institute of Research Engineering and Technology Quarterly Reviews Vol.4, No. 1, 202 1 13 Directions Respectively τ : Shear Stresses εij : Shear Strains in i -j Face λ ; Shell Thickness / Inside Radius (h/a) Cij : Elastic Moduli in General X, φ, Y : Non Dimensional Coordinate System in Longitudinal, Circumferential and Radial Directions Respectively Figure 1: A cylindrical shell showing dimensions, deformations and stresses Figure 2: Anisotropic materials Figure 3: Non -homogenous materials Asian Institute of Research Engineering and Technology Quarterly Reviews Vol.4, No. 1, 202 1 14 Note that since we have only four unknown constants, Ci*. ( i = 1, 4 ), the solution of this theory can accommodate four boundary conditions. We therefore must abandon some of the boundary conditions specified in equations (9.3). The chosen two boundary conditions are: Equation (7) where the external forces, >N" and T are as defined by equation (2) The equation (10) thus has a particular solution only Equation (8) From now on the superscript p wm indicate the particular solution. b) Theory Associated with Axial Length Scale (ah)1/2 and Circumferential Length Scale a Strain -Displacement Relations C*_ and C| arise from the integration Equation (9) Governing Equations Equation (10) From the first two equations of above we obtain the following relations for V and V Q: Equation (11) Asian Institute of Research Engineering and Technology Quarterly Reviews Vol.4, No. 1, 202 1 15 where C* i's (i = l,3) are the integration constants to be determined later and the coefficients D and E are defined as follows: We can now express the third equation of (15) in terms of a single variable V by substituting (16) into it. This yield Equation (12 ) where the constants N i (i = 1.2.3) are as follows: The homogeneous solution of equation (9.18) can be expressed as, Equation (13) Here, the A i' s are constants to be determined by the edge conditions speci fied earlier, N^ and N, are given by Governing Equations Equation (14) From the first two equations above, we obtain the following relations for V and V Q: where C* i's (i = l,3) are the integration constants to be determined later and the coefficients D and E are defined as follows: Asian Institute of Research Engineering and Technology Quarterly Reviews Vol.4, No. 1, 202 1 16 We can now express the third equation of (9.15) in terms of a single variable V by substituting (9.16) into it. This yield, Equation (15) where the constants N i (i = 1.2.3) are as follows: Equation (16) The homogeneous solution of above equation can be expressed as, Equation (17) Here, the A i' s are constants to be determined by the edge conditions speci fied. N^ and N, are given by Equation (18) where and £, is the dimensionless coordinate originating from the far edge of the shell defined as, The particular solution of Thus the complete solution is given by c) Theory Associated with Length Scales (ah)1/2 Asian Institute of Research Engineering and Technology Quarterly Reviews Vol.4, No. 1, 202 1 17 Strain -Displacement Relations Equation (19) Governing Equations The first two equations of the above will be as follows when solved for V z and Vθ yield, Equation (20) where the coefficients F, G, P. Q are defined as follows: and Q is given by On substituting the above equations into the third equation of constitutive equations we obtain differential equation for V only, where C 1* and C 2* are constants of integrat ion and Asian Institute of Research Engineering and Technology Quarterly Reviews Vol.4, No. 1, 202 1 18 Equation (21) The homogeneous solution of equation is representable in a similar fashion, where A 1's, and N 5 and N 6 are as defined previously, The particular solution is given by The complete solution is then Equation (22) d) Uniformly Valid Theory Strain -Displacement Relations Governing Equations Asian Institute of Research Engineering and Technology Quarterly Reviews Vol.4, No. 1, 202 1 19 Equation (23) The first two equations of above give us the following relations for Uz and Uθ: where C 1* and C 2* are the constants of integrations and the coefficients R, S, T, W are defined as follows: and On substituting (9.38) into the third equation of (9.36) we obtain a differential equation for U only, Equation (24) where APPPLICATION The homogeneous solution can be written as Asian Institute of Research Engineering and Technology Quarterly Reviews Vol.4, No. 1, 202 1 20 where A i ( i = l, 2, 3, 4) are the constants of integration to be determined and N5 , N6 are as defined previously. As we are now dealing with actual coordinates (except for the thickness coordinate), we define n as follows: where L is the actual length of the cylinder. The constants of integration are to be determined from the following edge conditions: The particular solution is and the complete solution of equation is then given by Equation (25) Having obtained the above solutions for each of the theories, numerical calculations are now carried out for a shell of the following dimensions: We thus have a thickness to radius ratio of Each of the layers is taken to be .025 in. thick and thus the dimensionless distances from the bottom of the first layer are given by Equation (26) As mentioned previously, each layer of the symmetric angle ply configu ration (elastic symmetry axes y are oriented at ( +γ , - γ, - γ, + γ ) is taken to be orthotropic with engineering elastic coefficients representing those for a boron/epoxy material system, Equation (27) Asian Institute of Research Engineering and Technology Quarterly Reviews Vol.4, No. 1, 202 1 21 Here direction 1 signifies the direction parall el to the fibers while 2 is the transverse direction. Angles chosen were = 0 , 15, 30, 45 and 60. Use of the results of Appendices A and B and the transformation equations (6) then yields the mechanical properties for the different symmetric angle ply configurations. We next apply the following edge loads: and take Shown in Figs. 3 - 5 is the variation of the dimensionless radial dis placement with the actual distance along the axis for the different theories. The reference surface for the chosen configuration is given by As mentioned above, the theory associated with length scales a is a membrane type theory and it's radial displacement is a function of the dimension less pressure and the two integration constants determined from the edge conditions (12). As Fig.3 demonst rates the radial displacement of this theory is constant over the entire length of the shell. The variation of the magnitude of radial displacement due to the change of cross -ply angle is almost identical compared to the other theories except for the fact that the theory cannot describe the deformation pattern due to the boundary cond itions while the other theories showing the radial deformations of the so-called edge effect zone. The theory associated with longitudinal length scale (ah)1/2 and circum ferential length scale a is similar to the axi -symmetric version of the theory of length scales (ah)1/2 in the following aspects: a) Expressions for strains and curvatures are identical as they were shown in (14) and (26). This is due to the fact that the theory of length scales (ah)1/2 are much simplified by the axi -symmetric property while the other theory is closer in fashion to the axi -symm etric behavior by it ’s nature because of the larger circumferential length scale we used for the theory, i.e. a. b) Although the expressions for the particular solutions are different as indicated in (24) and (33), the combined form of governing equations and the homogeneous solutions, as shown in (20) and (32), are identical in form. This is due to the same length scales being used in longitudinal direction, (ah)1/2 , for both theories and again, the axial symmetry. In obtaining the homogeneous solutions for both theories, it was assumed that the cylinder has such material properties and geometric dimensions so as to justify the decay type solutions (20) and (32). In order to have these decay type solutions we first must have that the value of term as shown in (22) must be real. Secondly, the dimensionless shell length (must be sufficiently larger compared to the axial length scale used in the basic formulation of the theories , (ah)1/2, so that interaction effects from the opposite edges may be neglected. The condition for satisfying this can be obtained by comparing the two decay terms in equation (20), exp ( -N5 x) and exp ( -N55 ζ). This leads to where Asian Institute of Research Engineering and Technology Quarterly Reviews Vol.4, No. 1, 202 1 22 The restriction of the cylinder length l to be larger than r is important in the analysis of cylindrical shells due to the difference, of nature of the solution. For the cylinder shorter than r, the edge conditions have an effect on each other and the solution is no longer of the decay type. Edge conditions in this case govern the deformation pattern as well as the magnitude. A short cylinder under external pressure and closed at both ends deform axi -symmetrically and ca n be considered a typical example of a problem where the solution has a decay length shorter than r. It must be noted here that unlike for isotropic homogeneous shells, the decay length r depends not only the shell geometry, h and a, but also on the mater ial pro perties of each laminate. CONCLUSION As stated previously, Figs. 4 and 5 show that the radial displacement of the shell at distances from the edge greater than r, from now on called the edge effect boundary layer, is nearly identical for both theories and close in magnitude to that of the solution which is obtained for length scales a. This is because, in the regions away from, the particular part of the solutions of gov erning equations dominate while the homogeneous solutions are more importan t within the boundary layer regions. Because the results shown in the figures are nearly identical for the problem considered, no numer ical calculations of the uniformly valid solution given above is carried out. It is also seen that wide variations in the magnitude of radial displace ment take place with change in the cross -ply angle. The maximum displacement occurs at γ =30 degree while the minimum displacement is at γ = 60 degree. In each case, the displacements increase with increas e in γ up to γ = 30 degree and thereafter decrease. The attached Figures show that the edge effect is sharper for small angle G than for larger ones. Similarly, deeper penetration of the edge effect is shown for small angles γ while weak and smooth edge effects are the case for large cross -ply angles. Shown in Figs. 6 and 7 are the dimensionless displacement of an isotropic material of elastic coefficient 30x106 psi and in Figs s h o w n are of single layer boron/epoxy composite we used for four layers case. Circumferential component of stress resultant is also shown in Figures