ISRN Civil Engineering

Volume 2012, Article ID 989785, 14 pages

http://dx.doi.org/10.5402/2012/989785

## Behaviour and Optimization Aids of Composite Stiffened Hypar Shell Roofs with Cutout under Free Vibration

Department of Civil Engineering, Meghnad Saha Institute of Technology, Kolkata 700107, India

Received 22 February 2012; Accepted 19 April 2012

Academic Editor: I. Smith

Copyright © 2012 Sarmila Sahoo. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

A scrutiny of the literature reveals that the free vibration characteristics of stiffened composite hypar shell with cutout are missing. So a generalized finite element formulation for the stiffened hyperbolic paraboloidal shells bounded by straight edges (commonly called as hypar shells) is attempted using an eight-noded curved quadratic isoparametric element for shell with a three-noded beam element for stiffener. Numerical problems of earlier investigators are solved as benchmark problems to validate the approach. A number of problems are further solved by varying the size of the cutouts and their positions with respect to the shell centre for different edge constraints. The results are presented in the form of figures and tables. The results are further analysed to suggest guidelines to select optimum size and position of the cutout with respect to shell centre considering the different practical constraints.

#### 1. Introduction

The advent of laminated composites in civil engineering applications has provided a new impetus to the researchers to explore the different aspects of composite structural elements including different forms of shells. A skewed hypar shell is aesthetically appealing and being doubly ruled is easy to cast. Moreover, this configuration can allow entry of north light and due to this advantage it finds use as roofing units in practical civil engineering situations. Shell structures that are normally thin-walled exhibit improved performances with stiffeners, particularly when the shell surface is provided with cutouts. Cutout is sometimes necessary in roof structure for the passage of light, to provide accessibility to other parts of the structure, for venting and also sometimes for alteration of resonant frequency. Basic knowledge of the free vibration characteristics of stiffened composite skewed hypar shell with cutout is essential for using these forms confidently.

As early as in 1982, Reddy [1] carried out the finite element analysis of composite plate with cutout and presented the effects of parametric variations on linear and nonlinear frequencies. Later in 1989, Malhotra et al. [2] presented the effect of fibre orientation and size of cutout on natural frequency on orthotropic square plates with square cutouts for different boundary conditions using the Rayleigh-Ritz method. Sivasubramonian et al. [3] reported free vibration of curved panels with cutout. They analysed the effect of cutouts on the natural frequencies of plates with some classical boundary conditions. Later Sivakumar et al. [4], Rossi [5], Huang and Sakiyama [6], and Hota and Padhi [7] studied free vibration of plate with various cutout geometries. Chakravorty et al. [8] reported some results in order to study the effect of concentric cutout on natural frequency of different shell options. In 1999, Sivasubramonian et al. [9] studied the free vibration characteristics of longitudinally stiffened square panels with symmetrical square cutouts by using the finite element method. The size of the cutout (symmetrically located) as well as curvature of the panels is varied. Hota and Chakravorty [10] published useful information about free vibration of stiffened conoidal shell roofs with cutout. Later, Nanda and Bandyopadhyay [11] investigated the effect of different parametric variation on nonlinear free vibration characteristics of cylindrical shell with cutout. The finite element model using an eight-noded C^{0} continuity, isoparametric quadrilateral element is used to study the dynamic behaviour. In a recent paper, Sahoo [12] carried out detailed free vibration analysis of composite unstiffened hypar shells with different practical boundary conditions.

Thus, it is found that results of free vibration characteristics of isotropic as well as composite stiffened composite shell panels with cutout are scanty in the existing body of literature. Only Chakravorty et al. [8] deals with some results of free vibration of hypar shell with concentric cutout, but to the best of author knowledge information regarding the free vibration behaviour of composite stiffened hypar shell with cutout is missing in the literature. In the present paper the free vibration of stiffened hypar shell with cutouts (Figure 1) is studied considering different boundary conditions. The variation of fundamental frequency due to change in eccentricity of cutout along and direction is also considered.

#### 2. Mathematical Formulation

##### 2.1. Finite Element Formulation for Shell

A laminated composite hypar shell of uniform thickness (Figure 2) and twist radius of curvature is considered. Keeping the total thickness the same, the thickness may consist of any number of thin laminae each of which may be arbitrarily oriented at an angle with reference to the *X*-axis of the coordinate system. An eight-noded curved quadratic isoparametric finite element (Figure 3) is used. The five degrees of freedom taken into consideration at each node include two in-plane and one transverse displacement and two rotations about the - and -axes. Sahoo and Chakravorty [13] reported in an earlier paper the strain displacement and constitutive relationships together with the systematic development of stiffness matrix for the shell element.

##### 2.2. Finite Element Formulation for Stiffener of the Shell

Three noded curved isoparametric beam elements (Figure 3) are used to model the stiffeners, which are taken to run only along the boundaries of the shell elements. In the stiffener element, each node has four degrees of freedom that is , , , and for -stiffener and , , , and for -stiffener. The generalized force-displacement relation of stiffeners can be expressed as (the notations have been defined in the nomenclature): where ; , ; .

Elasticity matrices are as follows: where and , , , and are explained in an earlier paper by Sahoo and Chakravorty [13].

Here the shear correction factor is taken as 5/6. The sectional parameters are calculated with respect to the midsurface of the shell by which the effect of eccentricities of stiffeners is automatically included. The element stiffness matrices are of the following forms: The integrals are converted to isoparametric coordinates and are carried out by 2-point Gauss quadrature. Finally, the element stiffness matrix of the stiffened shell is obtained by appropriate matching of the nodes of the stiffener and shell elements through the connectivity matrix and is given as The element stiffness matrices are assembled to get the global matrices.

##### 2.3. Element Mass Matrix

The element mass matrix for shell is obtained from the integral where in which Element mass matrix for stiffener element Here, is a diagonal matrix.

The mass matrix of the stiffened shell element is the sum of the matrices of the shell and the stiffeners matched at the appropriate nodes. The element mass matrices are assembled to get the global matrices.

##### 2.4. Modeling the Cutout

The code developed can take the position and size of cutout as input. The program is capable of generating nonuniform finite element mesh all over the shell surface. So the element size is gradually decreased near the cutout margins. One such typical mesh arrangement is shown in Figure 4. Such finite element mesh is redefined in steps and a particular grid is chosen to obtain the fundamental frequency when the result does not improve by more than one percent on further refining. Convergence of results is ensured in all the problems taken up here.

##### 2.5. Solution Procedure for Free Vibration Analysis

The free vibration analysis involves determination of natural frequencies from the condition

This is a generalized eigenvalue problem and is solved by the subspace iteration algorithm.

#### 3. Validation Study

The results of Table 1 show that the agreement of present results with the earlier ones is excellent and the correctness of the stiffener formulation is established. Free vibration of simply supported and clamped hypar shell with (0/90)_{4} lamination with cutouts is also considered. The fundamental frequencies of hypar shell with cutout obtained by the present method agree well with those reported by Chakravorty et al. [8] as evident from Table 2, establishing the correctness of the cutout formulation. Thus, it is evident that the finite element model proposed here can successfully analyse vibration problems of stiffened skewed hypar composite shells with cutout which is reflected by close agreement of present results with benchmark ones.

The present approach uses the improved first order approximation theory for thin shells [16] considering the radius of cross curvature. For this class of thin shells, a shear correction factor of unity is found to yield good results. It is observed that the results remain the same when analysis is repeated with the commonly used shear correction factor of .

#### 4. Results and Discussion

In order to study the effect of cutout size and position on the free vibration response, additional problems for hypar shells with 0/90/0/90 and +45/−45/+45/−45 lamination and different boundary conditions have been solved. The selection of the 0/90/0/90 and +45/−45/+45/−45 lamination is based on an earlier study by Sahoo and Chakravorty [13] which revealed that repeating 0/90 unit and +45/−45 unit more than once and keeping the total shell thickness constant does not improve the fundamental frequency to an appreciable extent. The positions of the cutouts are varied along both of the plan directions of the shell for different practical boundary conditions to study the effect of eccentricity of cutout on the fundamental frequency.

##### 4.1. Free Vibration Behaviour of Shells with Concentric Cutouts

Figures 5 and 6 furnish the results of nondimensional frequency () of 0/90/0/90 and +45/−45/+45/−45 stiffened hypar shells with cutout. The shells considered are of square plan form () and the cutouts are also taken to be square in plan (). The cutouts placed concentrically on the shell surface. The cutout sizes (i.e., ) are varied from 0 to 0.4 and boundary conditions are varied along the four edges. The stiffeners are placed along the cutout periphery and extended up to the edge of the shell. The boundary conditions are designated by describing the support clamped or simply supported as C or S taken in an anticlockwise order from the edge . This means a shell with CSCS boundary is clamped along and simply supported along and clamped along and simply supported along . The material and geometric properties of shells and cutouts are mentioned along with the figures.

###### 4.1.1. Effect of Cutout Sizes

From the figures it is seen that when a cutout is introduced to a stiffened shell the fundamental frequency increases in all the cases. This increasing trend is noticed for both cross-ply and angle-ply shells. This initial increase in frequency is due to the fact that with the introduction of cutout, numbers of stiffeners are increased from two to four in the present study. In order to study the effect of cutout size, in more details, the ratio of the fundamental frequency of a concentric punctured shell to that of a shell without cutout are expressed in percentage. The increase or decrease in percentage of fundamental frequency from the full shell is denoted by . Tables 3 and 4 contain such values for 0/90/0/90 and +45/−45/+45/−45 shells, respectively. Negative sign indicates decrease in frequency. It is evident from Tables 3 and 4 that in all the cases with the introduction of cutout with the frequencies increase. But further increase in cutout size, that is, when , fundamental frequency may increase or decrease. When the cutout size is further increased, but the number and dimensions of the stiffeners do not change, the shell surface undergoes loss of both mass and stiffness. It is evident from Tables 3 and 4 that when number of boundary constrains are more (shells with more clamped edges), loss of stiffness is more significant than loss of mass. Hence, fundamental frequency decreases except in case of clamped angle ply shell. For this shell when , fundamental frequency increases but further increase of cutout size decreases the fundamental frequency. But for shells with less number of boundary constraints, the fundamental frequency increases with the increase in cutout sizes. This trend is true for both cross- and angle-ply shells.

In such cases, as the cutout grows in size, the loss of mass is more significant than loss of stiffness, and hence, the frequency increases. From Tables 3 and 4, it is found that with introduction of a cutout of , in shell surface, the decrease in frequency is not more than 5%. This leads to the engineering conclusion that concentric cutouts with stiffened margins may be provided safely on shell surfaces for functional requirements up to .

###### 4.1.2. Effect of Boundary Conditions

The boundary conditions have been divided into four groups, so that the combinations in a particular group have equal number of boundary reactions. The groups are of the following forms: Group I contains CCCC shells. Group II contains CSCC, CCSC, and CCCS shells. Group III contains CSSC, CCSS, CSCS, and SCSC shells. Group IV contains SSSS shells.

As evident from Figures 5 and 6, fundamental frequencies of members belonging to different boundary combinations may be regrouped according to performance.

According to the values of (), Group III may be subdivided into Group IIIa and Group IIIb for both 0/90/0/90 and +45/−45/+45/−45 shells. Group I contains CCCC shells. Group II contains CSCC, CCSC, and CCCS shells. Group IIIa contains CSSC, CCSS shells. Group IIIb contains CSCS and SCSC shells. Group IV contains SSSS shells.

This observation indicates that the impact of number of boundary constraints is far more important than their arrangement. But the impact of arrangement of boundary constraints is not negligible in case of the shells which have two clamped edges. It is found that when two adjacent edges are clamped, the frequency attains greater value than when two alternate edges are clamped.

The frequencies are further studied and marks are given to the options of clamping the edges of a simply supported shell in order to gradually improve performances. Tables 5 and 6 furnish such clamping options for cross-ply and angle-ply shells, respectively. The scale is chosen like this, 0 is assigned to a simply supported shell and 100 to a clamped shell. These marks are furnished for cutouts with . These tables will help a practicing engineer. If one takes the frequency of a clamped shell as upper limit and that of the simply supported as lower limit, one can easily realize the efficiency of a particular boundary condition.

##### 4.2. Effect of Eccentricity of Cutout Position on Fundamental Frequency

To study the effect of cutout positions on fundamental frequencies, results are obtained for different locations of a cutout with . Each of the non-dimensional coordinates of the cutout centre (, ) is varied from 0.2 to 0.8 along both the plan directions so that the distance of a cutout margin from the shell boundary is not less than one tenth of the plan dimension of the shell. The study is carried out for all the twelve boundary conditions for both 0/90/0/90 and +45/−45/+45/−45 hypar shells. The ratio of the fundamental frequency of a shell with an eccentric puncture to that of a shell with concentric puncture (obtainable from Figures 5 and 6) expressed in percentage is denoted by . Tables 7 and 8 contain the value of for 0/90/0/90 and +45/−45/+45/−45 hypar shells.

It may be seen that the fundamental frequency is maximum when cutout is along the centre line of the shell. This is true for cross-ply shells with four edges simply supported or clamped. In case of clamped angle ply shells, similar trend is observed, but in case of a simply supported angle ply shells, the fundamental frequency is maximum when the cutout centre is along the diagonal. For the shells which have one edge simply supported and others clamped the fundamental frequency increases towards the simply supported edge. This observation is true for cross ply as well as angle ply shells. For the cross ply shells which have two opposite edges clamped and others opposite edges simply supported, the shell become stiffer when the cutout is placed along the centre line of the shell. But in case of angle ply shell with same boundary condition, the fundamental frequency of shell increases when the cutout centre is towards the clamped edges. It is further noticed that the shells with three edges simply supported and other clamped fundamental frequency is maximum when cutout centre is along the line which is equidistant from two simply supported edges. Along the other direction, when the cutout centre shifts towards the simply supported edge, which is opposite to be clamped edge, the shell becomes stiffer. But reverse is the case for angle ply shells. Here, the shifting of cutout centre towards the simply supported edges makes the shell flexible.

Tables 9 and 10 provide the maximum values of together with the corresponding values of and indicating the positions of the cutouts. These tables further specify rectangular zones within which the centre of the cutout may be varied so that the value of is always greater than or equal to 95 and less than or equal to 90 (resp.). It is to be noted that the cutout centre may be placed at some points beyond the zones indicated in Tables 9 and 10 to obtain similar values of , but only zones rectangular in plan are identified in the tables. These tables will enable practicing engineers to get an idea of the maximum eccentricity of a cutout which can be permitted if the fundamental frequency of a concentrically punctured shell is not to suffer a drastic reduction in value.

#### 5. Conclusions

From the present study, the following conclusions are drawn.(1)The finite element code used here is suitable for analysing free vibration problems of stiffened hypar shell roofs with cutouts as this approach produces results in close agreement with those of the benchmark problems.(2)Concentric cutouts may be provided safely on stiffened hypar shell surfaces for functional requirements up to .(3)The arrangement of boundary constraints along the four edges is far more important than their actual number as long as the free vibration stiffness is concerned.(4)Fundamental frequency undergoes marked improvement when the edge is converted to clamped from simply supported condition.(5)The relative free vibration performances of stiffened hypar shells for different combinations of edge conditions along the four sides are expected to be very useful in decision making for practicing engineers.(6)The information regarding the behaviour of stiffened hypar shell with eccentric cutouts for wide range of eccentricity and boundary conditions may be used as design aids by structural engineers.(7)As furnished, the specific zones within which the cutout centre may be varied such that the loss of fundamental frequency is within 5% with respect to a shell with concentric cutout also the specific zones which should be avoided (where the frequency decrease is greater than equal to 10%) will help to make a decision regarding the eccentricity of the cutout centre that can be allowed.

#### Notations

, : | Length and width of shell in plan |

, : | Length and width of cutout in plan |

: | Width of stiffener in general |

, : | Width of and stiffeners, respectively |

: | Strain displacement matrix of stiffener element |

: | Rise of hypar shell |

: | Depth of stiffener in general |

, : | Depth of and stiffeners, respectively |

: | Eccentricity of stiffeners with respect to midsurface of shell |

, : | Eccentricities of and stiffeners with respect to mid surface of shell |

, : | Elastic moduli |

, , : | Shear moduli of a lamina with respect to 1, 2, and 3 axes of fibre |

: | Shell thickness |

, : | Moment resultants of stiffeners |

: | Number of plies in a laminate |

, : | Number of stiffeners along and directions, respectively |

, : | Axial force resultants of stiffeners |

: | Ratio of the fundamental frequency of a concentric punctured shell to that of a shell without cutout expressed in percentage |

, : | Transverse shear resultants of stiffeners |

: | Ratio of the fundamental frequency of a shell with an eccentric puncture to that of a shell with concentric puncture expressed in percentage |

: | Radii of cross curvature of hypar shell |

, : | Torsion resultants of stiffeners |

, : | Axial and transverse translational degrees of freedom at each node of -stiffener element |

, : | Axial and transverse translational degrees of freedom at each node of -stiffener element |

, , : | Local co-ordinate axes |

, , : | Global co-ordinate axes |

: | Distance of bottom of the kth ply from midsurface of a laminate |

, : | Rotational degrees of freedom at each node of -stiffener element |

, : | Rotational degrees of freedom at each node of -stiffener element |

, : | Nodal displacement of stiffener element |

, : | Poisson’s ratios |

: | Density of material |

: | Natural frequency |

: | Non-dimensional natural frequency = . |

#### References

- J. N. Reddy, “Large amplitude flexural vibration of layered composite plates with cutouts,”
*Journal of Sound and Vibration*, vol. 83, no. 1, pp. 1–10, 1982. View at Google Scholar · View at Scopus - S. K. Malhotra, N. Ganesan, and M. A. Veluswami, “Vibration of composite plates with cut-outs,”
*Journal of Aeronautical Society of India*, vol. 41, pp. 61–64, 1989. View at Google Scholar - B. Sivasubramonian, A. M. Kulkarni, G. V. Rao, and A. Krishnan, “Free vibration of curved panels with cutouts,”
*Journal of Sound and Vibration*, vol. 200, no. 2, pp. 227–234, 1997. View at Google Scholar · View at Scopus - K. Sivakumar, N. G. R. Iyengar, and K. Deb, “Free vibration of laminated composite plates with cutout,”
*Journal of Sound and Vibration*, vol. 221, no. 3, pp. 443–470, 1999. View at Google Scholar · View at Scopus - R. E. Rossi, “Transverse vibrations of thin, orthotropic rectangular plates with rectangular cutouts with fixed boundaries,”
*Journal of Sound and Vibration*, vol. 221, no. 4, pp. 733–736, 1999. View at Google Scholar · View at Scopus - M. Huang and T. Sakiyama, “Free vibration analysis of rectangular plates with variously-shaped holes,”
*Journal of Sound and Vibration*, vol. 226, no. 4, pp. 769–786, 1999. View at Google Scholar · View at Scopus - S. S. Hota and P. Padhi, “Vibration of plates with arbitrary shapes of cutouts,”
*Journal of Sound and Vibration*, vol. 302, no. 4-5, pp. 1030–1036, 2007. View at Publisher · View at Google Scholar · View at Scopus - D. Chakravorty, P. K. Sinha, and J. N. Bandyopadhyay, “Applications of FEM on free and forced vibration of laminated shells,”
*Journal of Engineering Mechanics*, vol. 124, no. 1, pp. 1–8, 1998. View at Google Scholar · View at Scopus - B. Sivasubramonian, G. V. Rao, and A. Krishnan, “Free vibration of longitudinally stiffened curved panels with cutout,”
*Journal of Sound and Vibration*, vol. 226, no. 1, pp. 41–55, 1999. View at Google Scholar · View at Scopus - S. S. Hota and D. Chakravorty, “Free vibration of stiffened conoidal shell roofs with cutouts,”
*Journal of Vibration and Control*, vol. 13, no. 3, pp. 221–240, 2007. View at Publisher · View at Google Scholar · View at Scopus - N. Nanda and J. N. Bandyopadhyay, “Nonlinear free vibration analysis of laminated composite cylindrical shells with cutouts,”
*Journal of Reinforced Plastics and Composites*, vol. 26, no. 14, pp. 1413–1427, 2007. View at Publisher · View at Google Scholar · View at Scopus - S. Sahoo, “Free vibration of laminated composite hypar shell roofs with cutouts,”
*Advances in Acoustics and Vibrations*, vol. 2011, Article ID 403684, 13 pages, 2011. View at Publisher · View at Google Scholar - S. Sahoo and D. Chakravorty, “Finite element vibration characteristics of composite hypar shallow shells with various edge supports,”
*Journal of Vibration and Control*, vol. 11, no. 10, pp. 1291–1309, 2005. View at Publisher · View at Google Scholar · View at Scopus - A. Mukherjee and M. Mukhopadhyay, “Finite element free vibration of eccentrically stiffened plates,”
*Computers and Structures*, vol. 30, no. 6, pp. 1303–1317, 1988. View at Google Scholar · View at Scopus - A. N. Nayak and J. N. Bandyopadhyay, “On the free vibration of stiffened shallow shells,”
*Journal of Sound and Vibration*, vol. 255, no. 2, pp. 357–382, 2003. View at Publisher · View at Google Scholar · View at Scopus - A. Dey, J. N. Bandyopadhyay, and P. K. Sinha, “Finite element analysis of laminated composite paraboloid of revolution shells,”
*Computers and Structures*, vol. 44, no. 3, pp. 675–682, 1992. View at Google Scholar · View at Scopus