Free vibration analysis of simply supported beams with solid and thin-walled cross-sections using higher-order theories based on displacement variables
Introduction
The undamped free vibration analysis of structures has always been a major area of activity in structural design. The results of modal analyses are, in fact, of great interest in dynamic response analyses, acoustics, aeroelasticity and also to avoid resonance. Even today, for a certain class of structures, the most convenient way of conducting modal analyses is by means of idealizing the structure by simplified beam models. Beam models are easy to use and important tools for structural analysts. In aircraft structural design, for example, such models are still widely used in the modelling of helicopter rotor blades, aircraft wings, and propeller blades, amongst others.
The classical and oldest one-dimensional (1D) beam theory is that of Euler [1] and Bernoulli [2], hereinafter referred to as EBBM (Euler–Bernoulli beam model), which underwent further developments by Saint-Venant [3], [4] and Timoshenko [5], [6], hereinafter referred to as TBM (Timoshenko beam model). As it is well known, the EBBM does not account for transverse shear deformations, while the TBM incorporates a uniform shear distribution along the cross-section of the beam (see more details in [7]). However, these classical beam models have severe limitations (e.g., the impossibility of dealing with constrained warping and shear-bending couplings). Thus, there are several problems in the engineering practice that cannot be solved with these traditional tools. Deep and thin-walled beams are some examples for which advanced treatment might be necessary.
Many refined beam models can be found in the literature which overcome the shortcomings of classical models. A comprehensive review about existing beam and plate theories was published by Kapania and Raciti [8], [9], who investigated the vibrations, wave propagation, buckling and post-buckling behaviors. Another review about modern theories for beam structures was recently published by Carrera et al. [10]. However, a brief overview about refined 1D models is given here for the sake of completeness. Particular attention should be paid to the pioneering works by Sokolnikoff [11] and Timoshenko and Goodier [12]. Gruttmann et al. [13], [14], [15] computed shear correction factors for torsional and flexural shearing stresses in prismatic beams, arbitrary shaped cross-sections as well as wide- and thin-walled structures. The 3D elasticity equations based on Saint-Venant solution were reduced to beam-like structured by Ladevéze et al. [16], [17], [18] for high aspect ratio beams with thin-walled sections. Yu et al. [19], [20], [21] used the variational asymptotic solution of beams to build an asymptotic series. To enhance the description of the normal and shear stress of the beam, El Fatmi [22], [23] introduced improvements of the displacement models over the beam section by introducing a warping function. With the advent of the finite element method (FEM), various beam models were developed for validation purposes and a considerable overview was provided by Reddy [24], [25], whose works discussed both classical and higher-order 1D elements, together with the problem of shear-locking.
As far as the free vibration analysis is concerned, Eisenberger et al. [26] presented a method to compute the exact vibration frequencies of asymmetrical laminated beams. Three higher-order models to analyze the free vibrations of deep fiber reinforced composite beams were addressed by Marur and Kant [27], and the same authors extended this theory to study vibrations of angle-ply laminated beams by accounting for transverse shear and normal strain effects [28]. A higher order finite element model based on the classical lamination theory was developed by Ganesan and Zabihollah (see [29], [30]), and vibration response from laminated tapered composite beams was subsequently investigated. Kameswara et al. [31] studied a closed form solution with high-order mixed theory for free vibration analysis of composite beams. Numerical examples were computed for beams of various span to height ratios, and the results showed that their theories provide lower natural frequencies than those computed through Timoshenko model in case of thick sandwich beams.
All the publications mentioned above show that refined beam theories and the vibration analysis of slender structures still attract considerable attention of researchers and engineers. The current work presents a new method to deal with the free vibration behavior of beam structures. This method is based on the well-known Carrera unified formulation (CUF), which was introduced by Carrera et al. [32], [33], [34] for plates and shells. CUF was extended to beam structures by Carrera and Giunta [35] in 2010. Since then, various improvements of CUF have taken place and a brief overview is given below. The strength of CUF is that it allows the automatic development and compact formulation of any theory of structures by expressing the 3D displacement field as an expansion series of the generalized unknowns, which lie on the beam axis in the case of 1D models, through certain cross-sectional functions (see [7] for a comprehensive discussion about CUF). Several papers in the literature made use of Taylor series polynomials as cross-sectional functions, and the corresponding models were referred to as TE (Taylor expansion). TE models have demonstrated higher-order capabilities in dealing with various beam problems both in conjunction with FEM methods [36], [37], [38] and exact solutions [39], [40], [41]. More recently, Carrera and Petrolo [42], [43] adopted the Lagrange polynomials to discretize the cross-sectional kinematics and the resulting LE (Lagrange expansion) CUF models have been successfully used for the analysis of both metallic and laminated composite structures. Some of the main advantages of the LE models are that they only involve pure displacement unknowns and allow the component-wise analysis of complex structures [44].
In the previous literature about LE models, FEM was applied to solve the weak form governing equations. In the present work, for the first time, numerically exact solutions of the strong-form equations of motion of LE models for the free vibration analyses of solid and thin-walled structures are presented by assuming simply supported boundary conditions. The present methodology is said to be exact in the sense that it provides exact solution of the equations of motion of a structure once the initial assumptions on the displacement field have been made.
The paper is organized as follows: (i) First, the adopted notation and some preliminary relations are introduced in Section 2. (ii) CUF is then presented in Section 3, along with LE models for beams. (iii) Next, the governing differential equations and natural boundary conditions are derived in Sections 4 and 5. Here, by adopting simply supported boundary conditions, the differential problem is reduced into a linear eigenvalue problem in terms of CUF fundamental nuclei. (iv) Subsequently, a number of significant problems are treated in Section 6. (v) Finally, the main conclusions are outlined.
Section snippets
Preliminaries
The coordinate frame of the generic beam model is shown in Fig. 1. The beam has cross-section Ω and length L. The dimensions along y are . The displacement vector isin which and uz are the displacement components along x-, y- and z-axes, respectively. The superscript “T” represents a transpose. The stress, , and the strain, , components are grouped as follows:In the case of small displacements with respect to a
Unified formulation of beams
Within the framework of CUF (see [7], [46], [47]), the 3D displacement field is expressed as an expansion of the generalized displacements through generic functions where varies within the cross-section; is the generalized displacements vector and M stands for the number of terms in the expansion. According to the Einstein notation, the repeated subscript, τ, indicates summation. In this paper, Lagrange polynomials are used for functions. In
Governing equations
The principle of virtual displacements is used in this paper to derive the equations of motion:where δ stands for a virtual variation operator, represents the strain energy, and is the work of the inertial loadings.
The virtual variation of the strain energy isBy substitution of the geometrical relations (Eq. (3)), the material constitutive equation Eq. (5)), and the unified hierarchical approximation of displacements (Eq. (7)), and after integration by parts,
Analytical solution of the strong-form governing equations
By imposing simply supported boundary conditions, the above differential equations can be solved in analytical form. For doing this, the following generalized displacements, which correspond to a Navier-type solution, are supposed:where α isand m is a positive integer, which represents the half waves number along the beam axis and i is the imaginary unit. The components of vector are the new unknown
Numerical results
The present refined 1D model is compared and evaluated here by analyzing different problems. The results from the LE strong-form solutions are compared with those from classical beam models (EBBM and TBM) and 3D solid models by using the commercial code MSC Nastran. MSC Nastran models are based on HEXA-8 brick elements [48]. All of those models have the same boundary condition that is simply supported at both sides ().
Conclusions
Higher-order models for the analysis of the free vibrations of simply supported beams have been proposed in this paper. By exploiting the Carrera unified formulation (CUF) and Lagrange polynomials to discretize the beam cross-sectional kinematics, refined models with only displacement variables have been developed. The governing differential equations have been formulated in strong form by using the principle of virtual displacement. Subsequently, by assuming simply supported boundary
Acknowledgments
The first author acknowledges China Scholarship Council and the basic scientific research colleges in CAUC (No. 3122014D014).
References (48)
- et al.
A displacement method for the analysis of flexural shear stresses in thin-walled isotropic composite beams
Comput. Struct.
(2002) - et al.
New concepts for linear beam theory with arbitrary geometry and loading
Eur. J. Mech.—A/Solids
(1998) - et al.
Beamlike (Saint-Venant) solutions for fully anisotropic elastic tubes of arbitrary closed cross section
Int. J. Solids Struct.
(2004) Non-uniform warping including the effects of torsion and shear forces. Part Ia general beam theory
Int. J. Solids Struct.
(2007)Non-uniform warping including the effects of torsion and shear forces. Part IIanalytical and numerical applications
Int. J. Solids Struct.
(2007)On locking-free shear deformable beam finite elements
Comput. Methods Appl. Mech. Eng.
(1997)- et al.
Dynamic stiffness analysis of laminated beams using a first order shear deformation theory
Compos. Struct.
(1995) - et al.
Free vibration analysis of fiber reinforced composite beams using higher order theories and finite element modelling
J. Sound Vib.
(1996) - et al.
Vibration analysis of tapered composite beams using a higher order finite element. Part Iformulation
Compos. Struct.
(2007) - et al.
Vibration analysis of tapered composite beams using a higher order finite element. Part Iparametric study
Compos. Struct.
(2007)
Free vibrations of laminated beams using mixed theory
Compos. Struct.
Analysis of thickness locking in classical, refined and mixed multilayered plate theories
Compos. Struct.
Refined beam elements with arbitrary cross-section geometries
Comput. Struct.
Exact dynamic stiffness elements based on one-dimensional higher-order theories for free vibration analysis of solid and thin-walled structures
J. Sound Vib.
Refined dynamic stiffness elements applied to free vibration analysis of generally laminated composite beams with arbitrary boundary conditions
Compos. Struct.
De curvis elasticis
Mémoire sur la flexion des prismes
J. Math. Pures Appl.
Mémoire sur la torsion des prismes
Acad. Sci. l'Inst. Imp. Fr.
On the corrections for shear of the differential equation for transverse vibrations of prismatic bars
Philos. Mag.
On the transverse vibrations of bars of uniform cross section
Philos. Mag.
Beam Structures: Classical and Advanced Theories
Recent advances in analysis of laminated beams and plates. Part Ishear effects and buckling
AIAA J.
Recent advances in analysis of laminated beams and plates. Part IIvibrations and wave propagation
AIAA J.
Cited by (41)
Review of solution methodologies for structural analysis of composites
2024, European Journal of Mechanics, A/SolidsDynamic three-point bending tests under high loading rates
2023, Thin-Walled StructuresThermo-elastic damped nonlinear dynamic response of the initially stressed hybrid GPL/CNT/fiber/polymer composite toroidal shells surrounded by elastic foundation
2022, Composite StructuresCitation Excerpt :Various versions of FEM were developed in the Refs. [21–24] to carry out the thermal nonlinear dynamic analysis of the shells. Various approaches have been employed by many researchers for developing the finite shell elements such as Carrera’s Unified Formulation [25–29]. Some other numerical methods were employed by many researchers around the world, such as point collocation method [30], mesh-free [31–35], isogeometric analysis (IGA) [36,37], and the basis polynomials of the Ritz method [38,39].