Research Article  Open Access
Free Vibration Analysis of FibreMetal Laminated Beams via Hierarchical OneDimensional Models
Abstract
This paper presents a free vibration analysis of beams made of fibremetal laminated beans. Due to its attractive properties, this class of composites has gained more and more importance in the aeronautic field. Several higherorder displacementsbased theories as well as classical models (EulerBernoulli’s and Timoshenko’s ones) are derived, assuming Carrera’s Unified Formulation by a priori approximating the displacement field over the cross section in a compact form. The governing differential equations and the boundary conditions are derived in a general form that corresponds to a generic term in the displacement field approximation. The resulting fundamental term, named “nucleus”, does not depend upon the approximation order , which is a free parameter of the formulation. A Naviertype, closed form solution is used. Simply supported beams are, therefore, investigated. Slender and short beams are considered. Three and fivelayer beams are studied. Bending, shear, torsional, and axial modes and frequencies are presented. Results are assessed for threedimensional FEM solutions obtained by a commercial finite element code using threedimensional elements showing that the proposed approach is accurate yet computationally effective.
1. Introduction
In FibreMetal Laminates (FMLs), Figure 1, metallic and FibreReinforced Polymers (FRPs) are synergistically merged such that improved properties with respect to the individual constituents are obtained. Thanks to the presence of the metallic layer, FMLs present higher fracture toughness, yield stress, and compressive strength when compared to conventional FRPs [1]. Resistance to impact and environmental solicitations is improved by placing a metallic layer at top and the bottom of the structure. FMLs also offer several advantages upon monolithic metallic laminates such as high strength and stiffnesstoweight ratios and fibrebridging mechanism against crack growth in the metallic layers. In the case of heating due to fire, monolithic aluminium skin layers last about to seconds while the metallic layers in FMLs do not melt away and maintain a residual load carrying capability because the fibres may act as firewall [2]. Although FMLs per kilogram can cost five to ten times more than aluminium laminates, the cost of a finished product can be very close, even disregarding the positive effect of FMLs on direct operative costs (such as maintenance). Saving in weight of the structure (about ) can be obtained [3]. FMLs idea has been conceived in the late 70s at the Delft University of Technology when trying to lower the cost of full composites while maintaining the weight saving they comported. A breakthrough was obtained by inserting reinforcing fibres to the adhesive of already available bonded metallic sheets. In 1982, aramid reinforced aluminium laminates (ARALL) were commercially available. Unidirectional glass fibre embedded reinforced FMLs (GLARE) followed [3]. Besides aramid or glass fibres, FRPs made of unidirectional carbon (CARALL), vinylon (Virall), and both carbon and Kevlar or jute natural fibres in polymeric matrices have been presented. Woven fibrereinforced epoxy plies have been also considered. The metallic layers (whose thickness range is about to mm) are often made of aluminium alloys such as , , or . Other possible solutions account for magnesium or titanium. Mechanical or adhesive bonding is, then, used to connect the metallic and FRPs layers. FMLs have been conceived to meet the need of reduced weight and high performances in aeronautics. Fokker was one of the first to investigate the use of FMLs as panels for the wing of the F27. Following this example, several companies (e.g., Airbus, Boing, and Embraer) started to use FMLs [4]. For instance, FMLs are used in the fuselage skins of the Airbus A380. Titaniumcarbon FRPs solutions have been conceived for high temperature applications in supersonic transport and space applications. Over the years, different applications have been considered on FMLs tubes for chemical and nuclear industry, electronics packaging and fuel tanks, and blast resistant containers as well as sport goods (e.g., bike frames, tennis rackets).
Botelho et al. [4] presented a review of FMLs material properties, whereas Wu and Yang [5] thoroughly reviewed GLARE properties. Van Rooijen et al. [6] investigated the influence of constituents properties on FMLs mechanical properties by means of a rule of mixture in terms of the metal volume fraction. Cortes and Cantwell [7] studied numerically and experimentally the tensile properties of a lightweight FML based on a titanium alloy and carbon fibrereinforced FMLs. Sinke [8] presented a review of several GLARE structural solutions as well as manufacturing. FMLs offer a wide variety of choices due to the high number of involved parameters and possible configurations [9]. Due to the presence of FRPs layers, both the structure and the laminates have to be accounted for in the design process. Alderliesten [10] suggested that a topdown design approach (from structural requirements down to individual constituents) should be used. Sinke [8] outlined the need for multiscale modelling and design. Sawant and Muliana [11, 12] used this approach for the thermoviscoelastic behaviour of FMLs.
As laminated structures, advanced models developed for laminates can be used [15] provided they are extended to account for the hybrid behaviour deriving from the metallic layers [8]. As outlined in Alderliesten [10], FMLs mechanical response is generally predicted by the classical laminate theory [16]. Teply et al. [17] used a linear layerwise plate element for the analysis of ARALL laminates. A geometrical and material nonlinear solidlike shell element was formulated by Hashagen et al. [18]. Hagenbeek [19] proposed a solidlike shell element for the thermomechanical problem and results were assessed for experimental test. Hausmann et al. [20] presented some analytical and numerical solutions for determining the thermal residual stresses due to processing.
A free vibration analysis of FML beams via hierarchical onedimensional models is addressed in this paper. Models are derived via Carreras Unified Formulation that was previously derived for plates and shells [21–25] and extended to beams [26–31]. Through a concise notation for the displacement field, the governing differential equations and the corresponding boundary conditions are reduced to a “fundamental nucleus” that does not depend upon the approximation order. This latter can be assumed as a formulation free parameter. Displacementbased theories that account for nonclassical effects, such as transverse shear and cross section in and outofplane warping, can be formulated. Governing differential equations are solved via a Navier’s closed form solution. Slender as well as short beams are investigated. Three and fivelayer FML beams are studied. The presented models are validated towards threedimensional FEM solutions obtained by commercial code Ansys.
2. Preliminaries
In presenting the proposed model, the methodology in Giunta et al. [32] has been followed. A beam is a structure whose axial extension () is predominant if compared to any other dimension orthogonal to it. To identify the cross section (), it suffices to intersect the beam with planes that are orthogonal to its axis. A Cartesian reference system is adopted:  and axes are two orthogonal directions lying on . The coordinate is coincident to the axis of the beam and is bounded such that . Cross section geometry and reference system are reported in Figure 1. The displacement field is
in which , , and are the displacement components along , , and axes. Superscript ‘’ represents the transposition operator. The strain, , and stress, , tensors in Voigt notation are grouped into two vectors: , lying on planes orthogonal to where
and , lying on the cross section:
In the case of a linear analysis, the following straindisplacement geometrical relations hold:
Subscripts ‘’, ‘’, and ‘’, when preceded by comma, represent derivation versus the corresponding spatial coordinate. A compact vectorial notation can be adopted for (4):
where , , and are the following differential matrix operators:
and is the unit matrix. Under the hypothesis of linear elastic behaviour of the material, the generalised Hooke law holds:
in which is the material stiffness matrix. According to (2) and (3), it reads
Matrices , , , and in (8) are
where terms are the material stiffness coefficients:
where
and are Young’s moduli, Poisson’s ratios, and the shear moduli, respectively.
3. Hierarchical Beam Approximation
The displacement field is, a priori, assumed over the cross section in the following way:
represents the number of unknowns while is a generic expansion function over the cross section. The compact expression is based on Einstein’s notation: a repeated index indicates summation. This generic kinematic field can be used to formulate several displacementbased theories. Thanks to this notation, problem’s governing differential equations and boundary conditions can be derived in terms of a single “fundamental nucleus”. In this paper, Maclaurin’s polynomials are assumed as approximation functions . and as functions of the expansion order can be obtained via Pascal’s triangle as shown in Table 1; see Giunta et al. [13]. The actual governing differential equations and boundary conditions due to a fixed approximation order and polynomials type are obtained straightforwardly via summation of the nucleus corresponding to each term of the expansion. According to the previous choice of polynomial function, the generic order displacement field is
As far as the firstorder approximation order is concerned, the kinematic field is
Classical models, such as Timoshenko beam theory (TBT),
and EulerBernoulli beam theory (EBT),
are directly derived from the firstorder approximation model. Higherorder models yield a more detailed description of the shear mechanics (no shear correction coefficient is required), of the in and outofsection deformations, of the coupling of the spatial directions due to Poisson’s effect, and of the torsional mechanics than classical models do. EBT theory neglects them all, since it was formulated to describe the bending mechanics. TBT model accounts for constant shear stress and strain components. In the case of classical models, to contrast the phenomenon known in literature as Poisson’s locking, reduced Hook’s law should be used; i.e., the material stiffness coefficients should be corrected. For more details see Giunta et al. [28] and Carrera and Brischetto [33].
4. Governing Equations
The strong form of the governing differential equations and the boundary conditions is obtained via the Principle of Virtual Displacements (PVD):
where stands for a virtual variation, represents the internal work, and stands for the inertial work.
4.1. Virtual Variation of the Strain Energy
According to the grouping of the stress and strain components in (2) and (3), the virtual variation of the strain energy is considered as sum of two contributes:
where represents the volume of the beam. Using the geometrical relations, (5), the material constitutive equations, (8), and the unified hierarchical approximation of the displacements, (12), and integrating by parts, (18) reads
where the differential matrix operator operates on and the differential matrix operator operates on .
Equation (19) in a compact vectorial form reads
where is the differential linear stiffness matrix nucleus, that is, the representative term of the governing equations. is the differential linear matrix nucleus of the boundary conditions. The components of the differential linear stiffness matrix are
The generic term is a cross section moment:
As far as the boundary conditions are concerned, the components of are
4.2. Virtual Variation of the Inertial Work
The virtual variation of the inertial work is
Double dots stand for a second derivative with respect to time (). Accounting for (12), (24) becomes
Posing
and being Kronecker’s delta, then, the components of the inertial matrix are
Thus, the compact vectorial form of (25) is
4.3. Governing Equations’ and Boundary Conditions’ Fundamental Nuclei
Finally, the explicit form of the fundamental nucleus of the governing equations is
The boundary conditions are
where either displacement at one beam end is prescribed (the virtual variation is then nil) or the terms between brackets are nil (corresponding to a zero higherorder resultant coherent with a simply supported configuration).
5. Closed Form Analytical Solution
The resulting boundary value problem is solved via a Naviertype solution. The following displacement field is adopted:
where is
represents the halfwave number along the beam axis, is the angular frequency, and are the maximal amplitudes of the displacement components.
The displacement field in (31) satisfies boundary conditions (30) since
By substituting (31) into linear algebraic system (29), the following generalised eigenvalue problem is obtained:
In a compact vectorial form, (29) reduces to
where is the global mass matrix, is the global stiffness matrix, and is the amplitude of the displacement field.
6. Results and Discussion
Beams with a square cross section are considered. The beam is bounded such that
with being the length, m the width, and m the height of the beam; see Figure 1. A slenderness ratio as high as (slender beams) and as low as five (thick beams) is accounted for. ARALL and CARALL FMLs are considered. Three and fivelayer beams are investigated. Table 2 presents the properties of the FRP layers (see Ravishankar et al. [14]). Subscript “” is used to address the composite material. Mechanical properties of aluminium are: Young’s modulus equal to GPa, Poisson’s ratio equal to , and density equal to kg/ (subscript “” represents the metallic material). Modes with a halfwave () are studied and the natural frequency is put into the following dimensionless form:

In order to validate the models, the results are compared to the threedimensional finite element solution obtained through the commercial code Ansys. Reference solutions are addressed within the paper as “FEM ” where the subscript “” stands for the number of elements along the y or, equivalently, zaxis (the axial to cross section length ratio being equal to ten). The quadratic twentynode solid element “SOLID186” is used. In order to show the convergence of the reference threedimensional FEM solution (up to four significant digits for all the considered results), two meshes are considered: The first has and in the second one . As far as the computational costs are concerned, the degrees of freedom (DOFs) of the threedimensional FEM solution over a beam cross section as function of are
For , the DOFs are then 843. For a fixed approximation order , the total DOFs of the proposed solutions are
In the case of the highest considered expansion order (), they are .
Bending (on and plane), torsional, shear (on and plane), and axial modes are investigated. For slender beams (), the shear modes are much higher than the other considered modes (after all classical models hypotheses are satisfied in this case) and they are not reported. In all the tables, the frequencies are ordered according to the order of appearance in the reference threedimensional FEM solution: the first mode (“Mode 1”) is a bending mode in the plane, the second one (“Mode 2”) corresponds to a bending mode in the plane, the third one (“Mode 3”) is a torsional mode, the forth one (“Mode 4”) is an axial mode, the fifth one (“Mode 5”) is a shear mode in the plane, and, finally, the sixth one (“Mode 6”) is a shear mode in the plane. Higherorder beam model respects this order, whereas this is not always the case for classical and some loworder model. The mode swapping can be easily recognised by the fact that for a given model the frequencies are not in an ascending order. It should be noted that for two types of considered FML beams the order of apparition of the different modes is the same. Tables 3 and 4 present the dimensionless frequencies in the case of a threelayer slender beam. Classical theories accurately predict the bending and axial frequencies. A fourthorder model is required to provide good results for the torsional mode, the maximal difference being about , at worst. The case of a lengthtoheight ratio equal to ten is presented in Tables 5 and 6. For the ARALL beam, a thirdorder model accurately predicts the bending and the shear frequencies. A fourthorder theory is required for the torsional and axial mode, the error being less than compared with the reference solution. Classical theories yield good results for the bending and axial frequencies of CARALL beams. A fourthorder model is required to provide the torsional and shear frequencies when a difference with the reference threedimensional FEM solution of less than is required. The shear mode in the plane for the CARALL beam is the only exception. For this one, a seventhorder model is needed. Tables 7 and 8 show the frequencies for thick beams (). Classical theories do not accurately predict the different modes. For the ARALL beam, loworder models () will accurately predict the bending frequencies in the plane and the torsional and the shear frequencies in the plane. Higherorder theories are required for the other modes. A fifthorder theory predicts accurately the bending and shear modes in the plane, the error being less than compared to the reference solution. For the CARALL beam, a fourthorder approximation is required for the bending and the axial frequencies; the error being less than . For the torsional mode and the shear mode in the plane, a fifthorder is sufficient for the same accuracy. For the shear mode in the plane, a seventhorder model is needed.
 
Bending mode in the plane . Bending mode in the plane . Torsional mode. Axial mode. Mode not provided by the theory. 
 
Bending mode in the plane . Bending mode in the plane . Torsional mode. Axial mode. Mode not provided by the theory. 
 
Bending mode in the plane . Bending mode in the plane . Torsional mode. Axial mode. Shear mode in the plane . Shear mode in the plane . Mode not provided by the theory. 
 
Bending mode in the plane . Bending mode in the plane . Torsional mode. Axial mode. Shear mode in the plane . Shear mode in the plane . Mode not provided by the theory. 
 
Bending mode in the plane . Bending mode in the plane . Torsional mode. Axial mode. Shear mode in the plane . Shear mode in the plane . Mode not provided by the theory. 
 
Bending mode in the plane . Bending mode in the plane . Torsional mode. Axial mode. Shear mode in the plane . Shear mode in the plane . Mode not provided by the theory. 
Figures 2 and 3 present the axial and the shear mode in the plane of threelayer ARALL FML beams with lengthtoheight ratio equal to five. Those for CARALL are very similar and are not reported for the sake of brevity. From these figures, the difference in the material properties of the aluminium and the fibres is clearly visible since the cross section is no longer rigid.
Finally, results for fivelayer FML beams with a lengthtoheight ratio equal to ten are reported in Tables 9 and 10. Classical theories do not accurately predict the different modes. For an ARALL beam, a thirdorder theory is needed for the bending modes and the shear mode in the plane. Fourthorder and secondorder models predict accurately the torsional and the axial frequencies, respectively. For the shear mode in the plane, a seventhorder theory is required for a percentage error less than . For the CARALL beam, classical theories are accurate only for the axial mode, higherorder models being required for the other modes. For the bending modes, fifth and sixthorder models are needed. An eighthorder theory yields good results for the torsional mode, whereas shear modes call for a ninthorder theory.
 
Bending mode in the plane . Bending mode in the plane . Torsional mode. Axial mode. Shear mode in the plane . Shear mode in the plane . Mode not provided by the theory. 
 
Bending mode in the plane . Bending mode in the plane . Torsional mode. Axial mode. Shear mode in the plane . Shear mode in the plane . 