The Effect of Fiber Orientation and Laminate Stacking Sequences on the Torsional Natural Frequencies of Laminated Composite Beams
1Mechanical Design and Production Department Faculty of Engineering, Cairo University, Giza Egypt
2Industrial and Manufacturing Engineering Department Faculty of Engineering, Fayoum University, Fayoum Egypt
The composite materials are well known by their excellent combination of high structural stiffness and low weight. The main feature of these anisotropic materials is their ability to be tailored for specific applications by optimizing design parameters such as stacking sequence, ply orientation and performance targets. Finding free torsional vibrations characteristics of laminated composite beams is one of the bases for designing and modeling of industrial products. With these requirements, this work considers the free torsional vibrations for laminated composite beams of doubly symmetrical cross sections. The torsional vibrations of the laminated beams are analyzed analytically based on the classical lamination theory, and accounts for the coupling of flexural and torsional modes due to fiber orientation of the laminated beams are neglected. Also, the torsional vibrations of the laminated beams analyzed by shear deformation theory in which the shear deformation effects are considered. Numerical analysis has been carried out using finite element method (FEM). The finite element software package ANSYS 10.0 is used to perform the numerical analyses using an eight-node layered shell element to describe the torsional vibration of the laminated beams. The rotary inertia and shear deformation effects of the element are taken. The influence of fiber directions and stacking sequences of laminates on torsional natural frequencies were investigated. Also, the effects of boundary conditions are demonstrated. Numerical results, obtained by the ANSYS 10.0, classical lamination theory, and shear deformation theory are presented to highlight the effects of fibers orientation and layers stacking sequence on torsional frequencies of the beams. The results obtained by ANSYS are compared against the classical lamination theory, as well as shear deformation theory.
At a glance: Figures
Keywords: composite materials, laminated composite beams, torsional vibrations, shear deformation, finite element analysis
Journal of Mechanical Design and Vibration, 2013 1 (1),
Received October 19, 2013; Revised November 25, 2013; Accepted December 22, 2013Copyright © 2013 Science and Education Publishing. All Rights Reserved.
Cite this article:
- Hassan, Galal. A., Mohammed. A. Fahmy, and Ibrahim. M. Goda. "The Effect of Fiber Orientation and Laminate Stacking Sequences on the Torsional Natural Frequencies of Laminated Composite Beams." Journal of Mechanical Design and Vibration 1.1 (2013): 20-26.
- Hassan, G. A. , Fahmy, M. A. , & Goda, I. M. (2013). The Effect of Fiber Orientation and Laminate Stacking Sequences on the Torsional Natural Frequencies of Laminated Composite Beams. Journal of Mechanical Design and Vibration, 1(1), 20-26.
- Hassan, Galal. A., Mohammed. A. Fahmy, and Ibrahim. M. Goda. "The Effect of Fiber Orientation and Laminate Stacking Sequences on the Torsional Natural Frequencies of Laminated Composite Beams." Journal of Mechanical Design and Vibration 1, no. 1 (2013): 20-26.
|Import into BibTeX||Import into EndNote||Import into RefMan||Import into RefWorks|
The composite beam members have been increasingly used over the past few decades in the fields of aerospace, civil and mechanical engineering due to their excellent engineering features. A variety of structural components made of composite materials such as turbine blades, vehicle axles, aircraft wing, and helicopter blade can be approximated as laminated composite beams, which requires a deeper understanding of the vibration characteristics of the composite beams as mentioned by . The composite beams are generally used as structural components of light-weight heavy load bearing elements because of the high strength-to-weight and stiffness-to-weight ratios, the ability of being different strengths in different directions and the nature of being tailored to satisfy the design requirements of strength and stiffness in practical designs. The increased use of laminated composite beams requires a better understanding of vibration characteristics of these beams; it is quite essential in the design of composite beams subjected to dynamic loads. Due to the composite beams widely used in a variety of structures as well as their substantial benefits and great promise for future application, the dynamic behaviors of the laminated composite beams have received widespread attention and have been investigated extensively by many researchers. A number of researchers have been developed numerous solution methods to analysis the dynamic behaviors of the laminated composite beams [14, 16, 19].
In designing structures, it is critical to know the natural frequencies of the structure. If a natural frequency of the structure is close to an excitation frequency, then severe vibration of the structure could occur. This condition is called resonance and to avoid resonance, the natural frequencies of the structure must be altered by making suitable adjustments in the design. The study of such free vibrations (free because the structure vibrates with no external forces) is very important in finding the dynamic response of the elastic structure. Thus, in the dynamic analyses, it is quite essential to consider an overview of the free vibration characteristics, including the natural frequencies of these composite structures.
Free vibration analysis of laminated beams has been conducted by significant amount of research. Yõldõrõm and Kõral  studied the out-of-plane free vibration problem of symmetric cross-ply laminated beams using the transfer matrix method. Also, the effects of the rotary inertia and shear deformation are investigated under various boundary conditions. Banerjee, J.  investigated the free vibrations of axially loaded composite Timoshenko beams using the dynamic stiffness matrix method by developing an exact dynamic stiffness matrix of composite beams taking into account the effects of an axial force, shear deformation, and rotatory inertia. The theory includes the material coupling between the bending and torsional modes of deformations. Jun et al.  investigated the free vibration and buckling behaviors of axially loaded laminated composite beams having arbitrary lay-up using the dynamic stiffness method taking into account the influences of axial forces, Poisson effect, axial deformation, shear deformation, and rotary inertia. They developed the exact dynamic stiffness matrix by directly solving the governing differential equations of an axially loaded laminated beam. Abramovich and Livshits  studied the free vibration of non symmetric Cross-ply laminated Composite Beams based on Timoshenko type equations. The effect of coupled longitudinal and transversal displacements, shear deformation and rotary inertia are included in the analysis . Eisenberger, M. et al.  used the dynamic stiffness analysis and the first-order shear deformation theory to study the free vibration of laminated beams. Calim, F.  make study intended to analyze free and forced vibrations of non-uniform composite beams in the Laplace domain. Khdeir and Reddy  have been studied free vibrations of cross-ply laminated beams with arbitrary boundary conditions. Krishnaswamy, S. et al.  gave analytical solutions for the free vibration problem of laminated composite beams. Song and Waas  have been studied both buckling and free vibration analyses of laminated composite beams. They Song and Waas  also investigated the shear deformation effects. Yildirim, V.  used the stiffness method for the solution of the purely in-plane free vibration problem of symmetric cross-ply laminated beams with the rotary inertia, axial and transverse shear deformation effects included by the first-order shear deformation theory. Chandrashekhara and Bangera  investigated the free vibration of angle-ply composite beams by a higher-order shear deformation theory using the shear flexible FEM. Teh and Huang  presented two finite element models based on a first-order theory for the free vibration analysis of fixed-free beams of general orthotropy. Qiao Pizhong and Zou Guiping  presented an analytical study for dynamic behavior of pultruded fiber-reinforced plastic (FRP) composite cantilever I-beams based on a Vlasov-type linear hypothesis.
Several researchers have carried out studies on experimental and theoretical evaluations of flexural- torsional vibration analysis for FRP structural members. Lee and Kim  studied free vibration of a thin-walled laminated composite beam, where a general analytical model applicable to the dynamic behavior of a thin-walled channel section composite is developed. This model is based on the classical lamination theory, and accounts for the coupling of flexural and torsional modes for arbitrary laminate stacking sequence. Shadmehri, F. et al.  studied the flexural–torsional behavior of thin-walled composite beams with closed cross-section and a number of nonclassical effects, such as material anisotropy, transverse shear, are considered in the study. Kollar, LP.  investigated the analysis of Flexural-torsional vibration of open section composite beam with including shear deformation. Qiao et al.  presented a combined analytical and experimental approach to characterize the vibration behavior of pultruded Fiber-Reinforced Plastic (FRP) composite cantilever I-beams.
In engineering practice, we often come across the analysis of structures subjected to vibratory twisting loading, such as aerodynamic or asymmetric traffic forces. Also, composite structural elements consisting of a relatively weak matrix reinforced by stronger inclusions or of different materials in contact are of increasing technological importance in engineering. Steel beams or columns totally encased in concrete are most common examples, while construction using steel beams as stiffeners of concrete plates is a quick, familiar and economical method for long bridge decks or for long span slabs.
The extensive use of the aforementioned structural elements necessitates a rigorous dynamic analysis. Several researchers have dealt with torsional vibration of beams. Exact torsional vibration frequencies were presented by  for the case of circular cross-section shafts subjected to classical boundary conditions avoiding in this way warping effects. Eisenbeger, M.  studied the torsional vibration of open and variable cross section bars by derive analytical method is to form the dynamic stiffness matrix of the bar, including the effect of warping. Kameswara and Mirza  studied the problem of free torsional vibration and buckling of doubly symmetric thin-walled beams of open section, subjected to an axial compressive static load and resting on continuous elastic foundation. Kameswara RC.  investigated the torsional frequencies and mode shapes of generally constrained shafts and piping, where the exact frequency and normal mode expressions are derived for the free torsional vibrations of circular shafts or piping systems constrained by unsymmetrical torsional springs and carrying unequal rotational masses at either end. Evangelos and Sapountzakis  studied the Torsional vibrations of composite bars by (BEM) boundary element method which is developed for the nonuniform torsional vibration problem of doubly symmetric composite bars of arbitrary constant cross-section.
In the present study, the torsional vibration behaviors of symmetrical laminated composite beams are studied. The laminated beam is modeled and analyzed by the FEM. The commercial finite element program ANSYS 10.0 is used to perform a dynamic modelling to the laminated beams by performing an eigenvalue analysis. Mindlin eight-node isoparametric layered shell elements (SHELL 99) are employed in the modeling for describing the torsional vibrations of these beams. Also, analytical models are developed by classical lamination theory and shear deformation theory to study the torsional vibrations of the beams. In the analytical models, the flexural-torsional coupling effects are ignored and pure torsional vibrations are taken. The effects of fiber direction and laminate stacking sequence on the frequencies of torsional vibrations were investigated. Also, the effects of boundary conditions on the torsional frequencies of the laminated beams are demonstrated.
2. Material and Geometry
A generally laminated composite beam with a solid rectangular cross-section of doubly symmetrical cross sections, as shown in Figure 1, is considered to be studied. The laminated beam is made of many plies of orthotropic materials, and the principal material axes of a ply may be oriented at an arbitrary angle with respect to the x-axis. In the right-handed Cartesian coordinate system, the x-axis is coincident with the beam axis and its origin is on the mid-plane of the beam. The length, breadth and thickness of the beam are represented by L, b and h, respectively.
Glass fiber (E-Glass) is used as reinforcement in the form of unidirectional fibers with epoxy resin as matrix for the laminated composite beams. The mechanical properties for fiber and matrix are presented in Table 1 . For all finite element and analytical models, their associate material elastic properties were calculated analytically using the simple rule-of-mixtures as given in . More accurate values can be further obtained with some mechanical testing.
The constituent laminae were considered to be linear elastic and generally orthotropic therefore the concept of engineering constants was used to describe the laminae elastically. A certain set of elastic properties is required as input parameters for the finite element code and for the analytical models. The set of properties required as an input parameter at a material level were E1, E2, E3, G12, G13, G23, v12, v13 and v23 as shown in Table 1; Where 1, 2, and 3 are principal material directions.
3. Modal Analysis by Finite-Element Method, ANSYS
The beams were discretized using (type shell 99) finite element as shown in Figure 2, available in the commercial package ANSYS 10.0. This element has 8 nodes and is constituted by layers that are designated by numbers (LN - Layer Number), increasing from the bottom to the top of the laminate; the last number quantifies the existent total number of layers in the laminate (NL - Total Number of Layers).The element has six degrees of freedom at each node: translations in the nodal x, y, and z directions and rotations about the nodal x, y, and z-axes. The choice of shell 99 element type is based on layered applications of a structural shell model, and the type of results that need to be calculated.
xIJ = Element x-axis if ESYS is not supplied.
x = Element x-axis if ESYS is supplied.
A modal analysis will be carried out using ANSYS 10.0 finite element software to study the frequencies of free torsional vibrations of the mentioned composite laminated beam.
4. Dynamic Modeling by Analytical Methods
In the present paper, the free torsional vibrations of symmetric laminated beams are studied by the classical lamination theory and shear deformation theory.
When the cross-section has one plane of symmetry, one of the circular frequencies belongs to a flexural mode and the other two circular frequencies to flexural-torsional modes; while when the cross-section has two planes of symmetry, the three circular frequencies belong respectively to the flexural modes in the two planes of symmetry and to the pure torsion mode (when the axis of the beam does not bend). The bending–torsion coupling due to stiffness coupling presented in composite beams due to fiber orientation and stacking sequence is neglected.4.1. Classical Beam Theory
The oldest and the well-known beam theory is the Euler–Bernoulli beam theory (or classical beam theory—CBT), in which the shear deformation not included. Although this theory is useful for slender beams, it does not give accurate solutions for thick beams. The beams to be studied are orthotropic and its cross section has two axes of symmetry y and z. The mass is also symmetrical with respect to these axes, and, accordingly, the center of mass coincides with the origin of the y-z coordinate system, so that the flexural-torsional coupling not occurs.
A beam with two cross-sectional planes of symmetry may undergo flexural vibration in either of the two planes of symmetry and torsional vibration . Pure torsional vibrations are focused in this study.
Expressions for the torsional vibration ωBψi of long (>>) and short (<<) orthotropic beams are:
Torsional vibration of long beam is given by:
Torsional vibration of short beam is given by:
Where is the torsional stiffness of the beam; in N.m2, is the warping stiffness of the beam; in N.m4, θ is the polar moment of mass per unit length about the shear center, ρcomp is mass per unit volume, A is the area of the cross section, and µ4Bi and µ4Gi are parameters in the calculation of natural frequencies, which are given in Table 2.
The torsional frequencies of a beam of arbitrary length can be approximated by:
By using the previous two equations of torsional vibration, the torsional frequencies will be:
For symmetric orthotropic laminated beam previously mentioned; the torsional stiffness of the beam can be obtained by this relation,
and Warping stiffness of the beam can be obtained by this relation,
a11: element 1–1 of the laminate extensional compliance matrix (m/N).
d66: element 6–6 of the laminate bending compliance matrix (1/N. m).
The theory, based on the assumption that cross sections remain plane but not perpendicular to the axis is frequently called first-order shear theory. A beam, in which shear deformation is taken into account, is called a Timoshenko beam. In shear deformation theory the effect of the shear deformation is considered in torsional frequencies calculation as given by (La's zlo' and George, 2003).
Torsional vibration with shear deformation ωψi of short (<<) and long (>>and>>) orthotropic beams are:
Torsional vibration of short beam is given by,
Torsional vibration of long beam is given by,
The torsional circular frequency of a beam of arbitrary length can be approximated by:
By using Esq. (8) and (9) then,
Where the torsional shear stiffness is given by,
Where a66 is element 6–6 of the laminate extensional compliance matrix (m/N)
5. Numerical Results and Discussion5.1. Influence of Fiber Angle on Torsional Natural Frequencies
The influences of fiber orientation are investigated by modeling laminated beams of different lay-up construction of clamped – free boundary condition as shown in Figure 1. The analysis was performed to 8-layered symmetrically laminated beam with length 400 mm, width 40 mm and thickness 3.2 mm and the lamination scheme of beams is ranging from θ= 0° to 90°, in increments of 5°.
The results obtained after modeling the beams are presented in Figure 3. It presents the variation of the lowest three torsional natural frequencies of the beam with respect to fiber angle. From the results, It is seen that the torsional frequencies increase with increasing the fiber angle of the laminated beams until reach to significant values in the range from about θ=35° to θ=45°, then the torsional frequencies decrease gradually with increasing the fiber angle until reach to minimum value at θ=90°.
To investigate the influence of laminate stacking sequence, dynamic modeling is performed to 3 set of symmetrical laminates with a total of 8 layers and dimension of 400 mm length, 40 mm width and total thickness 3.2 mm. Each layer in the laminate has the same thickness. The lamination schemes of the beams are as follow: (0/90) 2S, (45/-45) 2S, and (45/ 45/0/90) S. The torsional frequencies are obtained by the analytical methods and by FE package, Table 3. [More Explanation of Table 3 is required regarding Observations and Comparison].
Figure 4, give the variation of the torsional natural frequencies of the laminates with respect to mode number for clamped free end condition, from the results it is already possible to verify the influence of the stacking sequence of the laminate on torsional vibration: the laminate with fibers at +/-45° has a larger torsional natural frequencies than the laminate of +/-45/0/90° and of 0/90° fibers. This was expected, since the natural frequencies are related to the stiffness of the structure and the (+/-45°) is much stiffer on torsion than the +/-45/0/90° and than 0/900 laminate.
The laminate of 0/90° fibers have the lowest torsional natural frequencies than the other lamination schemes and thus, because 50% of the fibers are oriented at 0° direction for 0/90° laminate, and thus appropriate for bending (Flexural Modes). This can be explained by the fact that the fibers oriented at 0° are more appropriate to flexural loads and the fibers oriented at 45° are more appropriate to torsional loads.
The mode shapes associated with the torsional natural frequencies of (0/90) 2S laminated beam are illustrated in Figure 5. They are deduced by FEM ANSYS for the first six torsional frequencies.
The FE analysis by using ANSYS is performed to investigate the influences of boundary conditions on torsional frequencies of the laminated beams. The analysis can be applied to the laminated beam previously mentioned with the same dimensions and geometry. The boundary conditions to be investigated in this study are as follow: C-F Clamped–Free, C-C Clamped–Clamped, C-S Clamped–Simply-Supported, and S-S Simply Supported - Simply Supported.
The influence of boundary conditions on torsional natural frequencies is investigated for the 1st torsional frequency as shown in Figure 6. From the results, It is seen that the clamped-clamped condition has a larger torsional frequency than other boundary fixations, and thus for all fiber angles. The torsional frequency for clamped-simply supported condition comes to be lower than clamped-clamped condition, then simply supported-simply supported comes to be lower than clamped-simply supported, and finally clamped-free condition comes to be lower than other supports.
This study is the first attempt to deal with the torsional natural frequencies of laminated composite beam with different lamination. For this, the influences of ply angle as well as laminate stacking sequences on the torsional natural frequencies of the laminated beams with doubly symmetrical cross- sections are investigated. Also, effects of boundary conditions on torsional natural frequencies are demonstrated.
It is observed that the results obtained by FEM using ANSYS are in good agreement with the analytical solutions, classical lamination theory and shear deformation theory. The analytical results by both classical lamination theory and shear deformation theory; the coupling effects due to bending–twisting coupling stiffness result from fiber orientation and laminate stacking sequences are excluded.
From the results, it is clear that changes in fiber angle as well as laminate stacking sequences yield to different dynamic behavior of the component, that is, different torsional natural frequencies for the same geometry, mass and boundary conditions. This gives the designer one additional degree of freedom to design the laminate - the possibility to change fiber orientations in order to get more (or less) structure stiffness in torsion. Different layers (plies) have different contributions to the overall stiffness of the beam depending on their locations from the mid-plane.
This possibility makes once more these materials very attractive since it makes possible to obtain the desired torsional natural frequencies without increasing mass or changing geometry. In practical applications, it means that if a torsional natural frequency excites the structure, the designer can change the material properties by changing the laminate stacking sequence, instead of re-design the complete structure. It is seen that the torsional frequencies increase with increasing the fiber angle of the laminated beams until reach to significant value at about θ=45°, then the torsional frequencies decrease gradually with increasing the fiber angle until reach to minimum value at θ=90°.
The finite element software package ANSYS is an efficient vibration prediction tool, because of its ability to model the laminated composite beam and reveal fundamental modal frequencies and modal shapes. Finally this study is useful for the designer in order to select the fiber orientation angle to shift the torsional natural frequencies as desired or to control the vibration level.
|||Abramovich A, Livshits A., (1994), “Free vibration of non-symmetric cross-ply laminated composite beams”, Journal of Sound and Vibration, Vol. 38, pp. 597-612.|
|||Banerjee, J.R., (1998), “Free vibration of axially loaded composite timoshenko beams using the dynamic stiffness matrix method”, Computers and Structures, Vol. 69, pp. 197-208.|
|||Calım FF, (2008), “Free and forced vibrations of non-uniform composite beams”, Composite Structures.|
|||Chandrashekhara K., Bangera K.M., (1992), “Free vibration of composite beams using a refined shear flexible beam element”, Computer and Structures, Vol. 43, pp. 719-727.|
|||Danial Gay, Hoa SV., Tsai SW., (2003), Composite materials design and applications, Boca Raton, London New York, Washington.|
|||Eisenberger M, Abramovich H, Shulepov O., (1995), “Dynamic stiffness analysis of laminated beams using a first-order shear deformation theory”, Composite Structure, Vol. 31, pp. 265-271.|
|||Eisenbeger M. (1997) “Torsional vibrations of open and variable cross section bars”, Thin-Walled Structures, Vol. 28, pp. 269-278.|
|||Evangelos J., Sapountzakis, (2005), “Torsional vibrations of composite bars by BEM”, Composite Structures, Vol. 70, pp. 229-239.|
|||Gorman DJ., (1975), “Free vibration analysis of beams and shafts”, New York.|
|||Jun L, Hongxing H, Rongying S., (2008), “Dynamic stiffness analysis for free vibrations of axially loaded laminated composite beams”, Composite Structures, Vol. 84, pp. 87-98.|
|||Kameswara Rao C., Mirza S., (1989), “Torsional vibrations and buckling of thin-walled beams on elastic foundation”, Thin-Walled Structures, Vol. 7, PP. 73-82.|
|||Kameswara RC., (1988), “Torsional frequencies and mode shapes of generally constrained shafts and piping” Journal of Sound and Vibration, Vol. 125, pp. 115-121.|
|||Kapuria S, Alam N., (2006), “Efficient layerwise finite element model for dynamic analysis of laminated piezoelectric beams”, Computer Methods in Applied Mechanics and engineering, Vol. 195, pp. 2742-2760.|
|||Khdeir A.A., Reddy J.N., (1994), “Free vibration of cross-ply laminated beams with arbitrary boundary conditions”, International Journal of Engineering Science, Vol. 32, pp. 1971-1980.|
|||Kollar LP., (2001), “Flexural-torsional vibration of open section composite beam with shear deformation”, International Journal of Solids and Structures, Vol. 38, pp. 7543-7558.|
|||Krishnaswamy S., Chandrashekhara K., Wu W.Z.B., (1992), “Analytical solutions to vibration of generally layered composite beams”. Journal of Sound and Vibration, Vol. 159, pp. 85-99.|
|||La'szlo' P., George S., (2003), Mechanics of composite structures, Cambridge University, New York, United States of America.|
|||Lee J., Kim S.E., (2002), “Flexural-torsional coupled vibration of thin-walled”, Composite Structures, Vol. 80, pp. 133-144.|
|||Matsunaga H., (2001), “Vibration and buckling of multilayered composite beams according to higher order deformation theories”, Journal of Sound and Vibration, Vol. 246, pp. 47-62.|
|||Qiao Pizhong, Zou Guiping., (2002), “Free vibration analysis of fiber-reinforced plastic composite cantilever I-beams”, Mechanics of Advanced Materials and Structure, Vol. 9, pp. 359-373.|
|||Qiao P., Zou G., Song G., (2002) “Analytical and experimental study of vibration behavior of FRP composite I-beams”, 15th ASCE Engineering Mechanics Conference, New York.|
|||Shadmehri F., Haddadpour H., Kouchakzadeh, M.A., (2007), “Flexural-torsional behavior of thin-walled composite beams with closed cross-section composite beams with channel sections”, Thin-Walled Structures, Vol. 45, pp. 699-705.|
|||Song SJ, Waas A., (1997), “Effects of shear deformation on buckling and free vibration of laminated composite beams”, Composite Structures, Vol. 37, pp.3 3-43.|
|||Teh KK, Huang CC., (1979), “The vibration of generally orthotropic beams-A finite element approach”, Journal of Sound and Vibration, Vol. 62, pp. 195-206.|
|||Vinson JR, Sierakowski RL., (2004), the behaviour of structures composed of composite materials, New York, Boston, Dordrecht, and Moscow.|
|||Yildirim V., (2000), “Effect of the longitudinal to transverse moduli ratio on the in-plane natural frequencies of symmetric cross-ply laminated beams by the stiffness method” Composite Structures, Vol. 50, pp. 319-26.|
|||Yõldõrõm V, Kõral E., (2000), “Investigation of the rotary inertia and shear deformation effects on the out-of-plane bending and torsional natural frequencies of laminated beams”, Composite Structures,Vol. 49, pp. 313-320.|