- About this Journal ·
- Abstracting and Indexing ·
- Aims and Scope ·
- Article Processing Charges ·
- Author Guidelines ·
- Bibliographic Information ·
- Citations to this Journal ·
- Contact Information ·
- Editorial Board ·
- Editorial Workflow ·
- Free eTOC Alerts ·
- Publication Ethics ·
- Recently Accepted Articles ·
- Reviewers Acknowledgment ·
- Submit a Manuscript ·
- Subscription Information ·
- Table of Contents

Journal of Engineering

Volume 2013 (2013), Article ID 230120, 18 pages

http://dx.doi.org/10.1155/2013/230120

## Dynamic Characters of Stiffened Composite Conoidal Shell Roofs with Cutouts: Design Aids and Selection Guidelines

Department of Civil Engineering, Heritage Institute of Technology, Kolkata 700107, India

Received 11 December 2012; Revised 3 May 2013; Accepted 14 May 2013

Academic Editor: Jun Li

Copyright © 2013 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

Dynamic characteristics of stiffened composite conoidal shells with cutout are analyzed in terms of the natural frequency and mode shapes. A finite element code is developed for the purpose by combining an eight-noded curved shell element with a three-noded curved beam element. The code is validated by solving benchmark problems available in the literature and comparing the results. The size of the cutouts and their positions with respect to the shell centre are varied for different edge constraints of cross-ply and angle-ply laminated composite conoids. The effects of these parametric variations on the fundamental frequencies and mode shapes are considered in details. The results furnished here may be readily used by practicing engineers dealing with stiffened composite conoids with cutouts central or eccentric.

#### 1. Introduction

Laminated composite structures are gaining wide importance in various fields of aerospace and civil engineering. Shell roof structures can be conveniently built with composite materials that have many attributes, besides high specific strength and stiffness. Among the different shell panels which are commonly used as roofing units in civil engineering practice, the conoidal shell has a special position due to a number of advantages it offers. Conoidal shells are often used to cover large column-free areas. Being ruled surfaces, they provide ease of casting and also allow north light in. Hence, this shell is preferred in many places, particularly in medical, chemical, and food processing industries where entry of north light is desirable. Application of conoids in these industries often necessitates cutouts for the passage of light, service lines, and also sometimes for alteration of resonant frequency. In practice, the margin of the cutouts must be stiffened to take account of stress concentration effects. An in-depth study including bending, buckling, vibration, and impact is required to exploit the possibilities of these curved forms. The present investigation is, however, restricted only to the free vibration behaviour. A generalized formulation for the doubly curved laminated composite shell has been presented using the eight-noded curved quadratic isoparametric finite element including three radii of curvature. Some of the important contributions on the investigation of conoidal shells are briefly reviewed here.

The research on conoidal shell started about four decades ago. In 1964, Hadid [1] analysed static characteristics of conoidal shells using the variational method. The research was carried forward and improved by researchers like Brebbia and Hadid [2], Choi [3], Ghosh and Bandyopadhyay [4, 5], Dey et al. [6], and Das and Bandyopadhyay [7]. Dey et al. [6] provided a significant contribution on static analysis of conoidal shell. Chakravorty et al. [8] applied the finite element technique to explore the free vibration characteristics of shallow isotropic conoids and also observed the effects of excluding some of the inertia terms from the mass matrix on the first four natural frequencies. Chakravorty et al. [9–11] published a series of papers where they reported on free and forced vibration characteristics of graphite-epoxy composite conoidal shells with regular boundary conditions. Later, Nayak and Bandyopadhyay [12–15] reported free vibration of stiffened isotropic and composite conoidal shells. Das and Chakravorty [16, 17] considered bending and free vibration characteristics of unpunctured and unstiffened composite conoids. Hota and Chakravorty [18] studied isotropic punctured conoidal shells with complicated boundary conditions along the four edges, but no such study about composite conoidal shells is available in the literature. Also, they did not furnish any information on vibration mode shapes. It is also seen from the recent reviews [19, 20] that dynamic characteristics of stiffened conoidal shells with cutout are still missing in the literature. The present study thus focuses on the free vibrations of graphite-epoxy laminated composite stiffened conoids with cutout both in terms of the natural frequencies and mode shapes. The results so obtained may be readily used by practicing engineers dealing with stiffened composite conoids with cutouts. The novelty of the present study lies in the consideration of vibration mode shapes of stiffened composite conoids in presence of cutouts.

#### 2. Mathematical Formulation

A laminated composite conoidal shell of uniform thickness (Figure 1), and radius of curvature , and radius of cross 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 -axis of the coordinate system. The constitutive equations for the shell are given by (a list of notations is separately given) where The force and moment resultants are expressed as The submatrices , , , and of the elasticity matrix are functions of Young’s moduli, shear moduli, and Poisson’s ratio of the laminates. They also depend on the angle which the individual lamina of a laminate makes with the global -axis. The detailed expressions of the elements of the elasticity matrix are available in several references including Vasiliev et al. [21] and Qatu [22].

The strain-displacement relations on the basis of improved first-order approximation theory for thin shell (Dey et al. [6]) are established as where the first vector is the midsurface strain for a conoidal shell and the second vector is the curvature.

#### 3. Finite Element Formulation

##### 3.1. Finite Element Formulation for Shell

An eight-noded curved quadratic isoparametric finite element is used for conoidal shell analysis. The five degrees of freedom taken into consideration at each node are , , , , *.* The following expressions establish the relations between the displacement at any point with respect to the coordinates and and the nodal degrees of freedom
where the shape functions derived from a cubic interpolation polynomial [6] are
The generalized displacement vector of an element is expressed in terms of the shape functions and nodal degrees of freedom as
that is,

###### 3.1.1. Element Stiffness Matrix

The strain-displacement relation is given by where The element stiffness matrix is

###### 3.1.2. Element Mass Matrix

*The element mass matrix is obtained from the integral*
where
in which

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

Three-noded curved isoparametric beam element (Figure 2) 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 where The generalized displacements of the -stiffener and the shell are related by the transformation matrix , where This transformation is required due to curvature of -stiffener, and is the appropriate portion of the displacement vector of the shell excluding the displacement component along the -axis.

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

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.

###### 3.2.1. 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 3 × 3 diagonal matrix.

Consider 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.

##### 3.3. 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 3. 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.

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

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

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

#### 4. Numerical Examples

The validity of the present approach is checked through solution of benchmark problems. The first problem, free vibration of stiffened clamped conoid, was solved earlier by Nayak and Bandyopadhyay [13]. The second is the free vibration of composite conoid with cutouts solved by Chakravorty et al. [11]. The results obtained by the present method, along with the published results, are presented in Tables 1 and 2, respectively.

Additional problems for conoids with cutouts are solved, varying the size and position of cutout along both of the plan directions of the shell for different practical boundary conditions. In order to study the effect of cutout size on the free vibration response, results for unpunctured conoids are also included in the study.

#### 5. Results and Discussion

It is found from Table 1 that the fundamental frequencies of stiffened conoids obtained by the present method agree well with those reported by Nayak and Bandyopadhyay [13]. Here, monotonic convergence is noted as the mesh is made progressively finer. Thus, the correctness of the stiffened shell element used here is established. It is evident from Table 2 that the present results agree with those of Chakravorty et al. [11], and the fact that the cutouts are properly modeled in the present formulation is thus established.

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

Tables 3 and 4 show the results for the non-dimensional fundamental frequency of composite cross-ply and angle-ply stiffened conoidal shells for different cutout sizes and various combinations of boundary conditions along the four edges. The shells considered are of square planform , and the cutouts are also taken to be square in plan . The cutout sizes (i.e., ) are varied from 0 to 0.4, and boundary conditions are varied along the four edges. Cutouts are concentric on shell surface. 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 that a shell with CSCS boundary is clamped along , simply supported along , clamped along , and simply supported along . The material and geometric properties of shells and cutouts are mentioned along with the figures.

###### 5.1.1. Effect of Cutout Size on Fundamental Frequency

From Tables 3 and 4, it is seen that when a cutout is introduced to a stiffened shell, the fundamental frequencies increase. This increasing trend continues up to for both cross- and angle ply shells except some angle ply shells with . The initial increase in frequency may be explained by the fact that when a cutout is introduced to an unpunctured surface, the number of stiffeners increases from two to four in the present study. When the cutout size is further increased, the number and dimensions of the stiffeners do not change, but the shell surface undergoes loss of both mass and stiffness. As the cutout grows in size, the loss of mass is more significant than that of stiffness, and hence the frequency increases. But for some angle ply shells with further increase in the size of the cutout, the loss of stiffness gradually becomes more important than that of mass, resulting in decrease in fundamental frequency. This leads to the engineering conclusion that cutouts with stiffened margins may always safely be provided on shell surfaces for functional requirements.

###### 5.1.2. Effect of Boundary Conditions on Fundamental Frequency

The boundary conditions may be divided into six groups, considering number of boundary constraints. The combinations in a particular group have equal number of boundary reactions. These groups are Group I: CCCC shells, Group II: CSCC, CCSC, and SCCC shells, Group III: CSSC, SSCC, CSCS, and SCSC shells, Group IV: CSSS, SSSC, and SSCS shells, Group V: SSSS shells, Group VI: Corner point supported shell.

It is seen from Tables 3 and 4 that fundamental frequencies of members belonging to the same groups of boundary combinations may not have close values. So, the different boundary conditions may be regrouped according to performance. According to the values of , the following groups may be identified.

For cross ply shells Group 1: Contains CCCC, CCSC, and SCSC boundaries which exhibit relatively high frequencies. Group 2: Contains CSCC, CCCS, CSSC, CCSS, CSCS, SSSC, CSSS, and SSCS which exhibit intermediate values of frequencies. Group 3: Contains SSSS and corner point supported boundaries which exhibit relatively low values of frequencies.

Similarly for angle ply shells: Group 1: Contains CCCC, CCSC, CSCC, CCCS, CSSC, CCSS, SCSC, and CSCS boundaries which exhibit relatively high frequencies. Group 2: Contains SSSC, CSSS, SSCS, and SSSS boundaries which exhibit intermediate values of frequencies. Group 3: Contains corner point supported shells which exhibit relatively low values of frequencies.

It is evident from the present study that the free vibration characteristics mostly depend on the arrangement of boundary constraints rather than their actual number. It can be seen from the present study that if the higher parabolic edge along is released from clamped to simply supported, there is hardly any change of frequency for cross ply shell. But for angle ply shells, if the edge along higher parabolic edge is released, fundamental frequency even increases more than that of a clamped shell. For cross ply shells, if the edge along or is released, that is, along the straight edges, frequency values undergo marked decrease. The results indicate that the edge along or should preferably be clamped in order to achieve higher frequency values, and if the edge has to be released for functional reason, the edge along and of a conoid must be clamped to make up for the loss of frequency. But for angle ply shells, if any two edges are released, the change in fundamental frequency is not so significant.

Tables 5 and 6 show the efficiency of a particular clamping option in improving the fundamental frequency of a shell with minimum number of boundary constraints relative to that of a clamped shell. Marks are assigned to each boundary combination in a scale assigning a value of 0 to the frequency of a corner point supported shell and 100 to that of a fully clamped shell. These marks are furnished for cutouts with These tables will enable a practicing engineer to realize at a glance the efficiency of a particular boundary condition in improving the frequency of a shell, taking that of clamped shell as the upper limit.

###### 5.1.3. Mode Shapes

The mode shapes corresponding to the fundamental modes of vibration are plotted in Figures 4 and 5 for cross-ply and angle ply shells, respectively. The normalized displacements are drawn with the shell midsurface as the reference for all the support conditions and for all the laminations used here. The fundamental mode is clearly a bending mode for all the boundary conditions for cross-ply and angle-ply shells, except for corner point supported shell. For corner point supported shells, the fundamental mode shapes are complicated. With the introduction of cutout, mode shapes remain almost similar. When the size of the cutout is increased from 0.2 to 0.4, the fundamental modes of vibration do not change to an appreciable amount.

##### 5.2. Effect of Eccentricity of Cutout Position

###### 5.2.1. Fundamental Frequency

The effect of eccentricity of cutout positions on fundamental frequencies is studied from the results obtained for different locations of a cutout with . The non-dimensional coordinates of the cutout centre were varied from 0.2 to 0.8 along each direction, so that the distance of a cutout margin from the shell boundary was not less than one tenth of the plan dimension of the shell. The margins of cutouts were stiffened with four stiffeners. The study was carried out for all the thirteen boundary conditions for both cross ply and angle ply shells. The fundamental frequency of a shell with an eccentric cutout is expressed as a percentage of fundamental frequency of a shell with a concentric cutout. This percentage is denoted by . In Tables 7 and 8, such results are furnished.

It can be seen that eccentricity of the cutout along the length of the shell towards the parabolic edges makes it more flexible. It is also seen that towards the lower parabolic edge, value is greater than that of the higher parabolic edge. This means that if a designer has to provide an eccentric cutout along the length, he should preferably place it towards the lower height boundary. The exception is there in some cases of cross ply shells. For cross ply shells with three edges simply supported when cutout shifts towards the lower parabolic edge (including the lower parabolic edge) the shell becomes stiffer. Again for corner point supported shells, values increase towards the parabolic edges and are maximum along lower parabolic edge. For cross ply shells, four out of thirteen boundary conditions yield the maximum value of along , four yield maximum values of along , and others show maximum values within to 0.5.

It is observed from Table 7 that if the eccentricity of a cutout is varied along the width, the shell becomes stiffer when the cutout shifts towards clamped edges. So, for functional purposes, if a shift of central cutout is required, eccentricity of a cutout along the width should preferably be towards the clamped straight edge. For shells having two straight edges of identical boundary condition, the maximum fundamental frequency occurs along . For corner point supported shells, the maximum fundamental frequency always occurs along the boundary of the shell. All these are true for cross ply shells only. For an angle ply shell, such unified trend is not observed, and the boundary conditions and the fundamental frequency behave in a complex manner as evident from Table 8. But for corner point supported angle-ply shells also, the maximum values of are along the boundary.

Tables 9 and 10 provide the maximum values of together with the position of the cutout. These tables also show the rectangular zones within which is always greater than or equal to 95 and 90. It is to be noted that at some other points, values may have similar values, but only the zone rectangular in plan has been identified. This study identifies the specific zones within which the cutout centre may be moved so that the loss of frequency is less than 5% or 10%, respectively, with respect to a shell with a central cutout. This will help a practicing engineer to make a decision regarding the eccentricity of the cutout centre that can be allowed.

###### 5.2.2. Mode Shapes

The mode shapes corresponding to the fundamental modes of vibration are plotted in Figures 6, 7, 8, and 9 for cross-ply and angle-ply shells of CCCC and CCSC shells for different eccentric positions of the cutout. As CCCC and CCSC shells are most efficient with respect to number of restraints, the mode shapes of these shells are shown as typical results. All the mode shapes are bending modes. It is found that for different position of cutout, mode shapes are somewhat similar, only the crest and trough positions change.

The present study considers the dynamic characteristics of stiffened composite conoidal shells with square cutout in terms of the natural frequency and mode shapes. The size of the cutouts and their positions with respect to the shell centre are varied for different edge constraints of cross-ply and angle-ply laminated composite conoids. The effects of these parametric variations on the fundamental frequencies and mode shapes are considered in details. However, the effect of the shape and orientation of the cutout on the dynamic characters of the conoid has not been considered in the present study. Future studies will evaluate these aspects.

#### 6. Conclusions

The following conclusions are drawn from the present study.(1) As this approach produces results in close agreement with those of the benchmark problems, the finite element code used here is suitable for analyzing free vibration problems of stiffened conoidal roof panels with cutouts. The present study reveals that cutouts with stiffened margins may always safely be provided on shell surfaces for functional requirements.(2) The arrangement of boundary constraints along the four edges is far more important than their actual number; so far the free vibration is concerned. The relative-free vibration performances of shells for different combinations of edge conditions along the four sides are expected to be very useful in decision making for practicing engineers.(3) The information regarding the behaviour of stiffened conoids with eccentric cutouts for a wide spectrum of eccentricity and boundary conditions for cross ply and angle ply shells may also be used as design aids for structural engineers.

#### 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 elements |

: | Depth of stiffener in general |

: | Depth of and stiffeners,respectively |

: | Element displacement |

Eccentricities of both and -direction stiffeners with respect to shell midsurface | |

: | Elastic moduli |

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

Shell thickness | |

Higher height of conoid | |

Lower height of conoid | |

: | Moment resultants |

: | Torsion resultant |

Number of plies in a laminate | |

: | Shape functions |

: | Inplane force resultants |

: | Inplane shear resultant |

: | Transverse shear resultant |

: | Radii of curvature and cross curvature of shell, respectively |

Translational degrees of freedom | |

Local coordinate axes | |

Global coordinate axes | |

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

Rotational degrees of freedom | |

: | Inplane strain component |

Angle of twist | |

: | Shearing strain components |

: | Poisson’s ratios |

Isoparametric coordinates | |

Density of material | |

: | Inplane stress components |

: | Shearing stress components |

Natural frequency | |

: | Nondimensional natural frequency. |

#### References

- H. A. Hadid,
*An analytical and experimental investigation into the bending theory of elastic conoidal shells [Ph.D. thesis]*, University of Southampton, 1964. - C. Brebbia and H. Hadid, Analysis of plates and shells using rectangular curved elements, CE/5/71 Civil engineering department, University of Southampton, 1971.
- C. K. Choi, “A conoidal shell analysis by modified isoparametric element,”
*Computers and Structures*, vol. 18, no. 5, pp. 921–924, 1984. View at Scopus - B. Ghosh and J. N. Bandyopadhyay, “Bending analysis of conoidal shells using curved quadratic isoparametric element,”
*Computers and Structures*, vol. 33, no. 3, pp. 717–728, 1989. View at Scopus - B. Ghosh and J. N. Bandyopadhyay, “Approximate bending analysis of conoidal shells using the galerkin method,”
*Computers and Structures*, vol. 36, no. 5, pp. 801–805, 1990. View at Scopus - A. Dey, J. N. Bandyopadhyay, and P. K. Sinha, “Finite element analysis of laminated composite conoidal shell structures,”
*Computers and Structures*, vol. 43, no. 3, pp. 469–476, 1992. View at Scopus - A. K. Das and J. N. Bandyopadhyay, “Theoretical and experimental studies on conoidal shells,”
*Computers and Structures*, vol. 49, no. 3, pp. 531–536, 1993. View at Scopus - D. Chakravorty, J. N. Bandyopadhyay, and P. K. Sinha, “Free vibration analysis of point-supported laminated composite doubly curved shells-A finite element approach,”
*Computers and Structures*, vol. 54, no. 2, pp. 191–198, 1995. View at Scopus - D. Chakravorty, J. N. Bandyopadhyay, and P. K. Sinha, “Finite element free vibration analysis of point supported laminated composite cylindrical shells,”
*Journal of Sound and Vibration*, vol. 181, no. 1, pp. 43–52, 1995. View at Publisher · View at Google Scholar · View at Scopus - D. Chakravorty, J. N. Bandyopadhyay, and P. K. Sinha, “Finite element free vibration analysis of doubly curved laminated composite shells,”
*Journal of Sound and Vibration*, vol. 191, no. 4, pp. 491–504, 1996. 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 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. N. Nayak and J. N. Bandyopadhyay, “Free vibration analysis and design aids of stiffened conoidal shells,”
*Journal of Engineering Mechanics*, vol. 128, no. 4, pp. 419–427, 2002. View at Publisher · View at Google Scholar · View at Scopus - A. N. Nayak and J. N. Bandyopadhyay, “Free vibration analysis of laminated stiffened shells,”
*Journal of Engineering Mechanics*, vol. 131, no. 1, pp. 100–105, 2005. View at Publisher · View at Google Scholar · View at Scopus - A. N. Nayak and J. N. Bandyopadhyay, “Dynamic response analysis of stiffened conoidal shells,”
*Journal of Sound and Vibration*, vol. 291, no. 3–5, pp. 1288–1297, 2006. View at Publisher · View at Google Scholar · View at Scopus - H. S. Das and D. Chakravorty, “Design aids and selection guidelines for composite conoidal shell roofs—a finite element application,”
*Journal of Reinforced Plastics and Composites*, vol. 26, no. 17, pp. 1793–1819, 2007. View at Publisher · View at Google Scholar · View at Scopus - H. S. Das and D. Chakravorty, “Natural frequencies and mode shapes of composite conoids with complicated boundary conditions,”
*Journal of Reinforced Plastics and Composites*, vol. 27, no. 13, pp. 1397–1415, 2008. View at Publisher · View at Google Scholar · View at Scopus - S. S. Hota and D. Chakravorty, “Free vibration of stiffened conoidal shell roofs with cutouts,”
*JVC/Journal of Vibration and Control*, vol. 13, no. 3, pp. 221–240, 2007. View at Publisher · View at Google Scholar · View at Scopus - M. S. Qatu, E. Asadi, and W. Wang, “Review of recent literature on static analyses of composite shells: 2000–2010,”
*Open Journal of Composite Materials*, vol. 2, pp. 61–86, 2012. - M. S. Qatu, R. W. Sullivan, and W. Wang, “Recent research advances on the dynamic analysis of composite shells: 2000–2009,”
*Composite Structures*, vol. 93, no. 1, pp. 14–31, 2010. View at Publisher · View at Google Scholar · View at Scopus - V. V. Vasiliev, R. M. Jones, and L. L. Man,
*Mechanics of Composite Structures*, Taylor and Francis, New York, NY, USA, 1993. - M. S. Qatu,
*Vibration of Laminated Shells and Plates*, Elsevier, London, UK, 2004. - S. Sahoo and D. Chakravorty, “Finite element bending behaviour of composite hyperbolic paraboloidal shells with various edge conditions,”
*Journal of Strain Analysis for Engineering Design*, vol. 39, no. 5, pp. 499–513, 2004. View at Publisher · View at Google Scholar · View at Scopus