#### Abstract

A thorough review of the existing literature reflects that forced vibration studies of laminated composite conoidal shells with complicated boundary conditions are missing. Hence, the present paper aims to fill the lacuna. A finite element code utilizing eight-noded doubly curved elements together with modified Sanders’ first approximation theory for thin shells is used to study the forced vibration behavior of moderately thin laminated composite conoidal shells subjected to three different uniformly distributed time-dependent forces. Newmark’s direct time integration method is used to solve the dynamic problem. Results obtained using the present code are compared with the values available in the literature, and a good agreement of the results confirms the accuracy of the proposed code. The transient responses of the laminated shell are studied meticulously for parametric variations like boundary conditions and stacking orders of cross and angle-ply laminates and are compared with bending responses of the shell to conclude on the necessity of the dynamic study.

#### 1. Introduction

Laminated composites gained popularity to fabricate plates and shells since second half of the last century. High strength/stiffness to weight ratio, low cost of fabrication, and better durability of the laminated composites popularized them in the weight-sensitive engineering applications. Moreover, the stiffness parameters of this material can be altered by varying the fiber orientations and lamina stacking sequences which made them a lucrative option to the engineers. Naturally, a good number of research reports started publishing on laminated plates and shells. A number of researchers like Reddy [1], Reddy and Chandrashekhara [2], Ribeiro [3], and Nanda and Bandyopadhyay [4, 5] reported bending and dynamic responses of laminated shell configurations.

A conoidal shell is doubly curved, easy to cast, added stiff surface which is used in the civil engineering industry to cover large column free open spaces as one sees in airport terminal buildings, shopping malls, and in car parking lots. Moreover, these shell structures allow entry of daylight and natural air which makes them popular roofing units to medicinal plants and food processing units. In a number of situations like drop hammer condition in the industry, snow loading in low temperature areas, the hit by wind borne debris develops dynamic forces on these shells. To ensure an uninterrupted service life of the shell roofs, their relative static and dynamic performances are needed to be understood in detail.

A number of researchers worked on isotropic and laminated composite conoidal shells. Static bending responses of isotropic conoidal shells were reported by Ghosh and Bandyopadhyay [6, 7]. The authors used finite element method and reported deflections and stress resultants of the conoidal shells. Chakravorty et al. [8] studied fundamental frequencies of composite conoidal shells with cut-outs. They also reported forced vibration responses of laminated conoidal shells with point supported boundaries. Later, Nayak and Bandyopadhyay [9] studied forced vibration responses of stiffened isotropic conoidal shells. Apart from the transient characteristics, the authors [10] also carried out free vibration analysis of composite conoidal shells. Das and Chakravorty [11, 12] worked on bending and free vibration of laminated conoidal shells with and without stiffeners, respectively. Bending characteristics of delaminated conoidal shells were reported by Kumari and Chakravorty [13, 14], while Pradyumna and Bandyopadhyay [15, 16], studied vibration and dynamic instability behaviors of laminated conoidal shells using the higher order shear deformation theory.

It is found from the literature survey that although a number of authors worked on dynamic analysis of laminated conoidal shells, none reported forced vibration responses of these shells with complicated boundary conditions. Hence, the present paper aims to study the deflections and stress resultants of laminated conoidal shells with complicated boundary conditions subjected to three different uniformly distributed time-dependent forces. Moreover, the static deflections and fundamental frequencies of these shells are also reported to understand the overall behavior of the conoidal shells.

#### 2. Mathematical Formulation

##### 2.1. Governing Differential Equation

Let us consider a laminated composite conoidal shell (Figure 1) of uniform thickness “” and radii of curvature “” and “.” Each of the thin laminae may be oriented at an angle “” with reference to the (Figure 3). Lagrange’s equation of motion is derived using Hamilton’s principle to express the governing differential equation of motion for an undamped forced vibration condition of the shell undergoing small displacements, where : work done by conservative forces in a shell element and is given by : kinetic energy of the shell element and is given by where represents time derivative of shell deflections with respect to time.

: work done by external forces and is given by where indicates external transverse load intensity.

##### 2.2. Laminate Constitutive Relationship

The laminate stress resultant vector is expressed using the following relationship: where is the laminate constitutive relationship matrix and is adopted from Kumari and Chakravorty [13]. While obtaining the shear stress resultants, proper shear correction factors (as used by Das and Chakravorty [11]) are applied. The laminate stress resultant vector is expressed as (see Figure 2).

##### 2.3. Finite Element Formulation

An eight-noded curved quadratic isoparametric element (Figure 3) with continuity is used in the proposed finite element code. The element displacement field is described as follows: where , , and denote displacements of the shell along , , and axes and and define rotations of the normal to the shell midsurface about and axes, respectively. The element degrees of freedom are expressed in terms of their nodal values by the following relationship: where denotes the shape functions of the element and adopted is from Das and Chakravorty [11].

Following the modified Sanders’ first approximation theory for thin shells, the midsurface strain vector of the shell is expressed as where The strain vector is related to the nodal degrees of freedom by the strain displacement matrix which is the same used here as those were reported by Das and Chakravorty [11].

Using the strain displacement and laminate constitutive relations, (2) can be expressed as Using (8), (3) is given by where Similarly, we get from (4), Using (11), (12), and (14), (1) can be expressed as The coefficient of the element acceleration vector of (15) represents the element mass matrix , and that of the element displacement vector represents the element stiffness matrix . The term on the right-hand side of (15) represents the element load vector.

Thus, (15) results in where The element mass matrix, stiffness matrix, and load vector are transformed to isoparametric coordinates and for numerical integration by Gauss quadrature rule. The element matrices are assembled to get the respective global matrices with due consideration for curvature of the conoidal shell surface [17].

Consider The dynamic equation of motion in the global form is

##### 2.4. The Static Problem

If the inertia force term of (19) is dropped, and the displacement and load vectors are assumed to be time independent, then the following equation of static equilibrium is obtained: The previous equation is solved by the Gauss elimination method [18].

##### 2.5. The Free Vibration Problem

If the load vector of (19) is dropped, the equation of free vibration is obtained as The free vibration analysis involves determination of natural frequencies and is solved by subspace iteration algorithm [18].

##### 2.6. The Forced Vibration Problem

The global force vector in (19) is transient in nature and solved using Newmark’s direct time integration method [18].

#### 3. Numerical Problems

In order to verify the correctness of the static bending and free vibration formulations of the proposed code, the static deflections and fundamental frequencies of laminated spherical shell are evaluated and compared with the values reported by Reddy [1]. The comparisons are shown in Tables 1 and 2 with which the material and geometric properties of the shell are also reported as footnote. A comparison of results for varying mesh size are also reported in the tables.

The accuracy of the proposed code in predicting dynamic deflections of a conoidal shell is verified by comparing the present results with the values reported by Nayak and Bandyopadhyay [9]. The authors worked on isotropic conoidal shell, and the present code is modified accordingly by assigning equal values of elastic and shear moduli along all the directions. The comparison of results is furnished in Figure 4 where the material and geometric properties of the isotropic shell are also reported.

Apart from the verification problems, the present code is further utilized to study forced vibration responses of composite conoidal shells with two different complicated boundary conditions (Figure 5). The shells are subjected to three different time-dependent forces which are depicted in Figure 6. Keeping the total laminate thickness constant, symmetric and antisymmetric two, three, and four-layered laminates of cross and angle-ply laminations are considered. The nondimensional downward static deflections, nondimensional fundamental frequencies, and dynamic magnification of deflections for all three load cases of the conoidal shell are reported in Table 3. The nondimensional dynamic displacements are furnished in Table 4. The dynamic magnification factors (henceforth referred as DMF) are calculated as the ratios of maximum dynamic responses to the corresponding static ones. The material and geometric properties of the conoidal shells are furnished with the tables. In case of bending analysis, the magnitude of the static load is considered equal to the peak step load value of load-time history considered here.

**(a) A step load**

**(b) A rectangular pulse load**

**(c) Half sinusoidal pulse load**

#### 4. Results and Discussion

Static deflections and fundamental frequencies obtained using the present formulation are in close agreement [refer to Tables 1 and 2] with the values reported by Reddy [1]. Close agreement of static deflections as reported in Table 1 ensures accurate bending formulation for laminated shells in the current code. The deflection values monotonically increase as the finite element mesh is made finer in steps because a coarse mesh indicates a rigid modelling of the plate stiffness. For the laminate, showing orthropy about the midplane, interestingly the deflection improves greatly when the mesh changes from a to a grid. For other laminates, the same trend is observed but the mesh does not yield results which are so widely different from the converged ones. Correct incorporation of mass matrix of laminated shells is proved by the results of the second benchmark problem indicated in Table 2. The frequency values decrease gradually as the mesh is made finer for the reasons explained previously. As it is found that a mesh grossly overestimates the plate stiffness for lamination, as a natural outcome the frequency value for this grid is much higher than the converged result. The response curves furnished in Figure 4 show that the natures of published response curves and those obtained by present approach are in close agreement and thus validates correct incorporation of Newmark’s method in the present code.

The shells are ranked from 1 to 16 in terms of static and dynamic deflections with the first rank being given to the shell option showing the least value. In cases of fundamental frequency, the first rank goes to the shell exhibiting the highest value. Such ranks are provided in parentheses in Tables 3 and 4 and are very helpful to understand the relative behaviors of the shell combinations comprehensively.

The complicated boundary conditions considered here (Figure 5) have equal numbers of support degrees of freedom locked but arranged in a different manner. Naturally, a practicing engineer may be curious to know the relative performances of the conoidal shells in terms of these two-edge conditions. It is interesting to note from Table 3 that for any given lamination the CCSS boundary condition shows less static deflections when they are compared with the values obtained for the SSCC one. The conoidal shell is curved along one of its plan directions and has straight edges along the other plan direction. The arch direction of the shell is relatively stronger than the beam direction as it combines the bending and axial rigidity together to sustain external loads. But unlike the cylindrical shell, the conoidal shell has its curved boundaries at different elevations. This is why the bending stiffness of the shell is enhanced substantially on clamping the lower parabolic arch of the shell, and the static deflections decrease. The fundamental frequencies of the conoidal shell listed in the same table also vary in a similar manner as it was observed in case of the static deflections, although the symmetric angle-ply laminations show marginally higher values for the SSCC shell. The and laminates yield 3.71% and 2.18% higher frequency values, respectively, for the SSCC boundary condition. This observation leads to conclude that clamping the lower arch of the conoid is advantageous to achieve higher bending and dynamic stiffnesses.

The dynamic deflections of all the shell options considered here are studied for three different load cases: step load of infinite duration (load case I), step load of finite duration (load case II), and half sinusoidal pulse load (dynamic load case III). The results are furnished in Table 4. The results show that the dynamic deflections follow the same pattern as it was noted in cases of static deflections and fundamental frequencies where the CCSS shell performed better than the SSCC one. Although the magnification factors for all the load cases listed in Table 3 vary in a different manner where in a number of cases the SSCC shell shows lesser DMF values, it is found from the table that the DMF values for dynamic load cases I and II are all greater than 2. The load-time history of load cases I and II is identical up to 1 second, and after that the dynamic load corresponding to load case I continues to act with the same intensity, and for load case II, the load is withdrawn. Interestingly, the magnification factors are same for both the load cases or marginally higher for load case II. In load case III, both the application and withdrawn of pressure are gradual and this is why the DMF values are marginally greater than unity.

It is also interesting to observe that for both the boundary conditions considered here the angle-ply shells yield lesser values of static and dynamic deflections when they are compared with their counterparts among cross-ply shells. It is further noted form Table 3 that among the symmetric and antisymmetric angle-ply laminates, the four-layered antisymmetric one () shows the best overall performance for the CCSS boundary condition. This means that in order to get the least possible static and dynamic deflections of a conoidal shell, an engineer’s choice will be the CCSS boundary condition with the four-layered antisymmetric angle ply laminate (). The above conclusion is further reinforced by the fact that in terms of fundamental frequency also this particular laminate with the aforesaid boundary condition is the best available choice.

To study the transient responses of the composite conoidal shells after removal of the dynamic load, a more detailed transient analysis is carried out, and the results are furnished graphically in Figures 7, 8, 9, 10, and 11. The shell option showing the best overall performance is chosen for that study. The nondimensional dynamic deflections and force and moment resultants of the shell are studied. The shell responses are reported for all three load cases, and for each of them the study is carried out upto 2 seconds, that is, for a period of 1 second even after the load is withdrawn to see the after effect of the transient load. Dynamic responses of the shell are studied at the location of the maximum static displacement, and the respective and coordinates are presented with the figures. The loading conditions for load cases I and II are in fact identical up to time = 1 second. This is why the response curves up to 1 second for these two load cases coincides. In both of these cases since the load is suddenly applied, the elastic rebound of the structure is quite strong, and sharp closely spaced crests and troughs are noted. For load case III, however, as the load is very gradually applied and withdrawn and load is of lesser magnitude, the deflection response curve almost resembles the shape of the loading curve.

The dynamic deflections for load cases I and II are reported along with the static deflection in Figure 7, and they are observed to oscillate around the static deflection upto 1 second. After the sudden removal of the dynamic load in case of load case II, the shell response magnifies drastically and even results in upward movements. This is why the DMF values for load cases I and II are almost the same or marginally higher for load case II. This observation also brings out the importance of forced vibration analysis as static simplification of dynamic problems using a suitable factor cannot account for such reversal of dynamic displacements. It is also important to note here that in case of the dynamic loads, its time variation is more significant than the magnitude. The dynamic deflections for load case I continue its same vibrating pattern for 2 seconds. On the other hands for load case III when the load is applied and withdrawn gradually, the shell vibrates with positive and negative peaks of equal magnitude and in a frequency which is nearly equal to its fundamental frequency.

The force and moment resultants of the composite shell (Figures 8, 9, 10, and 11) follow the same pattern as it was observed in case of the dynamic deflections. Reversals of all the force and moment resultants are noted in the figures although it is more prominent for the moment resultants compared to than that of the force resultants in case of load case I where the dynamic load continued to act for 2 seconds.

#### 5. Conclusion

The following conclusions can be drawn from this study.(1)The present code accurately incorporates the finite element method to solve static, free, and forced vibration problems of laminated composite conoidal shells as it is indicated by the solutions of the benchmark problems.(2)It is concluded that among the boundary conditions considered in the present study, the lower parabolic arch of the conoidal shell should be clamped to achieve higher bending and dynamic stiffnesses.(3)It is further indicated that among the considered laminations, the laminate shows the best overall performance by yielding the least static and dynamic deflections and the highest fundamental frequency.(4)When the dynamic responses of the composite shell are studied for varying time, it is found that the deflections and stress resultants show reversal in nature. Hence a full dynamics analysis is needed to be carried out to get the magnitudes of the reversed stresses, imposing a factor of safety on the static values cannot take care of this phenomenon.

#### Notations

Area of the shell | |

: | Length and width of shell in plan |

: | Flexural rigidity matrix of the shell |

: | Global displacement vector |

: | Global acceleration vector |

: | Young’s moduli |

and 3: | Local coordinates of a lamina |

: | Shear moduli |

: | Higher and lower heights of the conoidal shell, respectively |

: | Shell thickness |

: | Moment resultants per unit length |

: | Nondimensional moment resultants |

: | Nondimensional torsion resultant |

: | Force resultants per unit length |

: | Nondimensional inplane force resultants |

: | Nondimensional inplane shear resultant |

: | Number of elements |

: | Transverse shear resultants per unit length |

: | Peak value of the transient loads |

: | Radius of cross curvature of conoidal shell |

: | Radius of curvature of conoidal shell along |

: | Radius of the spherical shell |

: | Translational degree of freedoms along , , and directions, respectively |

: | Volume of the shell |

: | Nondimensional transverse displacement |

: | Global coordinates axes |

: | Rotations of the shell about and , respectively |

: | Shear strains |

: | Inplane strains along and , respectively |

: | Angle of lamination with respect to of the conoidal shell |

: | Poisson’s ratio |

: | Natural coordinates of isoparametric elements |

: | Mass density |

: | Inplane stresses along and , respectively |

: | Shear stresses |

: | Curvatures of conoidal shell due to load |

: | Fundamental frequency in radian/sec |

: | Nondimensional fundamental frequency . |

#### Acknowledgment

The first author gratefully acknowledges the financial assistance of the Council of Scientific and Industrial Research, India, through the Senior Research Fellowship vide Grant no. 09/096 (0686) 2k11-EMR-I.