On the natural frequency and vibration mode of composite beam with nonuniform crosssection
Botong Li^{1} , Longlei Dong^{2} , Liangliang Zhu^{3} , Xi Chen^{4}
^{1}College of Mechanical Engineering, Beijing University of Technology, Beijing 100124, China
^{1, 3, 4}International Center for Applied Mechanics, SV Lab, School of Aerospace, Xi’an Jiaotong University, Xi’an 710049, China
^{2}State Key Laboratory for Strength and Vibration of Mechanical Structures, School of Aerospace, Xi’an Jiaotong University, Xi’an 710049, China
^{4}Department of Earth and Environmental Engineering, Columbia University, New York, NY 10027, USA
^{4}Corresponding author
Journal of Vibroengineering, Vol. 17, Issue 5, 2015, p. 24912502.
Received 20 January 2015; received in revised form 25 March 2015; accepted 5 May 2015; published 15 August 2015
JVE Conferences
In this paper, the vibratory properties and expression of natural modes of laminated composite beam with variable crosssection ratios of elastic modulus and density along the axis of the beam have been investigated via theoretical analysis. Based on the generalized Hamilton principle, the longitudinal and transverse vibration equations have been deduced by the means of variational method. Then, the natural frequencies of longitudinal and transverse vibration modes have been obtained using the method of power series, which agree well with finite element simulations. The firstorder natural frequencies of longitudinal and transverse of composite beams are plotted as a function of the elastic modulus or densities difference of two components. With distinct material characteristics, the effect of shape factor on the first and second order lateral modes of composite beam is also revealed. In addition, the study shows that the boundary conditions impose a strong effect on the shape factor. The method presented in this paper is not only suitable for the laminated composite beam with variable crosssection, but will also be applicable to more general cases of composite beams of complex geometry and component in vibration mechanics. This controllable vibration performance achieved in this paper may shed some light on and stimulate new architectural design of composite engineering structures.
Keywords: BernoulliEuler composite beams, natural frequencies, vibration mode, method of power series.
1. Introduction
In mechanical, aerospace and civil structures and vehicles, laminated structural components have received wide applications thanks to their distinct advantages of high strength, stiffnesstoweight ratios and high bending rigidity [1, 2]. An important example is the widely used steelconcrete components in bridges [3, 4]. For structural applications, the eigenmodes and vibration characteristics of composite beams are of great importance [5]. Extensive theoretical and numerical approaches were carried out to explore the natural frequencies of laminated structures. Tseng et al. [6] applied the stiffness analysis method involving shear deformation and rotation inertia to determine the natural frequencies of laminated beams with arbitrary curvatures. Banerjee [7] presented the frequency equation and mode shapes of composite Timoshenko beams by symbolic computing. Rao et al. [8] proposed a higherorder mixed theory for determining natural frequencies of several laminated simply supported beams. Chen et al. [9] combined the state space method and the differential quadrature method for freely vibrating laminated beams. An excellent overview of recent advances in straight and curved composite beam models may be found in [10]. Wu and Chen [11] and Qu et al. [12] compared a few free vibration modes of laminated beams using the shear deformationbased theories. By using the RayleighRitz (RR) method, Gunda et al. [13] studied large amplitude vibration of laminated composite beam with symmetric and asymmetric layup orientations.
Most previous studies focused on composite beams with uniform axial properties, while composite beams with variable crosssections may be needed to better fit the stress distribution or to achieve a better structural performance [1417]. Present study emphasizes the normal mode analysis of composite beam with nonuniform constituent crosssection properties. The attention is limited to composite beams with a constant overall geometrical section dimension, but the height ratio between the components varies along the axial direction (Fig. 1), the configuration of which is similar to the widely used and studied scarf joint technique for fabricating composite structures [1820]. It is noted that the algorithms presented here are also applicable to study other composite beams with more complex crosssection properties, for example, the jagged joint interface between two components.
Any general structural deformation can be regarded as a superposition of its fundamental normal modes, therefore, the present study emphasizes the normal mode analysis of composite beam with nonuniform constituent crosssection properties. To focus on the effect of constituent ratio, the attention is limited to composite beams with a constant overall geometrical section dimension, but the height ratio between the components varies along the axial direction (Fig. 1).
Many previous analytical and computational studies focused on the nonlinear vibration equations, using differential quadrature method (DQM) [2122], dynamic stiffness method [23], finite element method (FEM) [24] and finite difference method [25], etc. In this paper, we adopt the variational method to obtain the kinetic equations of composite beams with variable crosssection ratios of elastic modulus, and to deduce the analytical expression of natural modes of longitudinal and transverse vibrations. Numerical validations are carried out using finite element simulations.
Though this study focuses on a rather simple composite beam with two components attached on their inclined surface (Fig. 1), the algorithms presented here can be easily extended to be applicable to studying composite beams of complex geometry and multiple components. One example is to study the vibration mode and frequency of the helicopter rotor blades which consist of multiple different material components and each component’s crosssection may be varying along its length direction. This study is also intended to stimulate new design of composite structures. A widely used composite beam structure is the composite bridges usually composed of a uniform layer of steel, a uniform layer of concrete and probably a uniform layer of fiber reinforced polymer deck. We believe that a more appropriate design of the bridge with varying crosssection of each material component (for example, increasing the thickness of the concrete layer in the middle region of the bridge) may be able to spread the load on the bridge, thus exploiting the inherent advantages of each of its material. Besides, it would be routine to analyze the natural frequency and vibration mode of a number of composite structures with various joint interface geometries [14, 17] based on the formulas presented here.
2. Modeling and mathematics formulation
Consider a beam of rectangular cross section, composed of two materials with their ratio changing along the axial direction. This is equivalent to bonding two rectangular beams with variable cross section, (Fig. 1). Such a configuration may be applicable to the composite components in bridges to allow a better fitting of the bending moment distribution of the bridge due to external loading or used as a scarf joint for composite structures. It is worth noting that our normal mode analysis algorithms presented below are not limited to this configuration, but can be easily adapted to study other composite beams with more complex $x$$y$ crosssection properties. All beams considered here are BernoulliEuler beams.
The elastic modulus, crosssectional area, moment of inertia, length and density of the upper beam are ${E}_{1}$, ${A}_{1}\left(x\right)$, ${I}_{1}\left(x\right)$, $l$, ${\rho}_{1}$, respectively, while the ones of lower beam are ${E}_{2}$, ${A}_{2}\left(x\right)$, ${I}_{2}\left(x\right)$, $l$, ${\rho}_{2}$. The width of the beam is $b$. Thus, the crosssectional areas of upper and lower beams are:
${A}_{2}\left(x\right)=\left({a}_{20}kx\right)b,$
while ${a}_{10}+{a}_{20}=a$ and $k=({a}_{20}{a}_{10})/l$.
The moment of inertia of beams can be derived as:
According to the neutral plane assumption, the longitudinal (axial) deformation $u$ of upper beam and lower beam are the same. We can also obtain the similar conclusion with the transverse deformation $w$. Shear deformation is neglected.
The kinetic energy of the system is:
The energy of the system is:
Using variational principle, one obtains:
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}={\int}_{{t}_{0}}^{{t}_{1}}{\int}_{0}^{1}\left[{\rho}_{1}{A}_{1}\left(x\right)+{\rho}_{2}{A}_{2}\left(x\right)\right]\left(\ddot{u}\delta u\ddot{w}\delta w\right)dxdt,$
$+\delta {\int}_{{t}_{0}}^{{t}_{1}}\frac{1}{2}{\int}_{0}^{1}\left[{E}_{1}{I}_{1}\right(x)+{E}_{2}{I}_{2}(x\left)\right]{\left(\frac{{\partial}^{2}w}{\partial {x}^{2}}\right)}^{2}dxdt$
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}={\int}_{{t}_{0}}^{{t}_{1}}{\int}_{0}^{1}\frac{\partial}{\partial x}\left\{\left[{E}_{1}{A}_{1}\left(x\right)+{E}_{2}{A}_{2}\left(x\right)\right]\frac{\partial u}{\partial x}\right\}\delta udxdt$
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}+{\int}_{{t}_{0}}^{{t}_{1}}{\int}_{0}^{1}\left[{E}_{1}{I}_{1}\left(x\right)+{E}_{2}{I}_{2}\left(x\right)\right]\left(\frac{{\partial}^{2}w}{\partial {x}^{2}}\right)\delta \left(\frac{{\partial}^{2}w}{\partial {x}^{2}}\right)dxdt$
$={\int}_{{t}_{0}}^{{t}_{1}}\left\{{\int}_{0}^{1}\frac{\partial}{\partial x}\left\{\left[{E}_{1}{A}_{1}\left(x\right)+{E}_{2}{A}_{2}\left(x\right)\right]\frac{\partial u}{\partial x}\right\}\delta udx\right.$
$\left.\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}+{\int}_{0}^{1}\frac{{\partial}^{2}}{\partial {x}^{2}}\left\{\left[{E}_{1}{I}_{1}\right(x)+{E}_{2}{I}_{2}(x\left)\right]\left(\frac{{\partial}^{2}w}{\partial {x}^{2}}\right)\right\}\delta wdx\right\}dt.$
Based on the generalized Hamilton principle, the following equations with independent $\delta u$ and $\delta w$ are obtained:
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\left[{E}_{1}{A}_{1}\right(x)+{E}_{2}{A}_{2}(x\left)\right]\frac{{\partial}^{2}u}{\partial {x}^{2}}=0,$
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}+2\left[{E}_{1}\frac{\partial {I}_{1}\left(x\right)}{\partial x}+{E}_{2}\frac{\partial {I}_{2}\left(x\right)}{\partial x}\right]\frac{{\partial}^{3}w}{\partial {x}^{3}}+\left[{E}_{1}\frac{{\partial}^{2}{I}_{1}\left(x\right)}{\partial {x}^{2}}+{E}_{2}\frac{{\partial}^{2}{I}_{2}\left(x\right)}{\partial {x}^{2}}\right]\frac{{\partial}^{2}w}{\partial {x}^{2}}=0.$
Fig. 1. Beams with different crosssectional areas
3. Mathematical transformation and calculation
3.1. The longitudinal deformation equation
The longitudinal deformation equation is considered first. From Eq. (8), it is assumed that the main vibration mode of beam is:
Combining Eqs. (8) and (10):
where ${\alpha}_{1}=({E}_{1}{a}_{10}+{E}_{2}{a}_{20})b$, ${\alpha}_{2}=({E}_{1}{E}_{2})kb$, ${\alpha}_{3}=({\rho}_{1}{a}_{10}+{\rho}_{2}{a}_{20})b$, ${\alpha}_{4}=({\rho}_{1}{\rho}_{2})kb$. Obviously, ${\alpha}_{2}\ne \text{0}$. Assume that:
It follows that ${\alpha}_{2}>0$.
In Eq. (11), ${\alpha}_{2}/({\alpha}_{1}+{\alpha}_{2}x)$ and $({\alpha}_{3}+{\alpha}_{4}x){\omega}^{2}/({\alpha}_{1}+{\alpha}_{2}x)$ are analytical as $x\in \left[0,l\right]\text{,}$ which may be solved using power series. Suppose that:
Substitute Eq. (13) into Eq. (11):
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}+\sum _{n=0}^{\mathrm{\infty}}{\alpha}_{3}{\omega}^{2}{a}_{n}{x}^{n}+\sum _{n=0}^{\mathrm{\infty}}{\alpha}_{4}{\omega}^{2}{a}_{n}{x}^{n+1}=0.$
To obtain the solution, make each power coefficient of $x$ equal to zero, leads to:
${a}_{0}$ and ${a}_{1}$ are different constants.
Assuming that the beam is stationary initially:
Substitute Eq. (16) into Eq. (10), one deduces $\phi =0$.
Suppose the one end of the beam is fixed while the other is free. This kind of boundary condition is taken as an example in the calculating process:
The fundamental modes of the system are:
The frequency equation is:
Thus:
3.2. The transverse deformation equation
The transverse deformation is explored next. The main vibration equation is:
Substitute it into Eq. (9):
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}+(6{\beta}_{1}x+2{\beta}_{2})Y\mathrm{\text{'}}\mathrm{\text{'}}({\alpha}_{3}+{\alpha}_{4}x){\omega}^{2}Y=0,$
where:
${\beta}_{1}=\left({E}_{1}{E}_{2}\right)\frac{b{k}^{3}}{3},\mathrm{}\mathrm{}\mathrm{}{\beta}_{2}=b{k}^{2}\left[{E}_{1}\left({a}_{10}\frac{a}{2}\right){E}_{2}\left(\frac{a}{2}{a}_{20}\right)\right],$
${\beta}_{3}=bk\left[{E}_{1}{\left({a}_{10}\frac{a}{2}\right)}^{2}{E}_{2}{\left(\frac{a}{2}{a}_{20}\right)}^{2}\right],$
${\beta}_{4}=b\left\{\left[{E}_{1}{\left({a}_{10}\frac{a}{2}\right)}^{3}{E}_{2}{\left(\frac{a}{2}{a}_{20}\right)}^{3}\right]+\frac{({E}_{1}+{E}_{2}){a}^{3}}{8}\right\}/3.$
Suppose Eq. (22) has the solution in the form of (see Appendix for justification):
Substitute the above equation into Eq. (22), yields:
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}+\sum _{n=4}^{\mathrm{\infty}}{\beta}_{3}n\left(n1\right)\left(n2\right)\left(n3\right){a}_{n}{x}^{n3}+\sum _{n=4}^{\mathrm{\infty}}{\beta}_{4}n\left(n1\right)\left(n2\right)\left(n3\right){a}_{n}{x}^{n4}$
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}+\sum _{n=3}^{\mathrm{\infty}}6{\beta}_{1}n\left(n1\right)\left(n2\right){a}_{n}{x}^{n1}+\sum _{n=3}^{\mathrm{\infty}}4{\beta}_{2}n\left(n1\right)\left(n2\right){a}_{n}{x}^{n2}$
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}+\sum _{n=3}^{\mathrm{\infty}}2{\beta}_{3}n(n1)(n2){a}_{n}{x}^{n3}+\sum _{n=2}^{\mathrm{\infty}}6{\beta}_{1}n(n1){a}_{n}{x}^{n1}$
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}+\sum _{n=2}^{\mathrm{\infty}}2{\beta}_{2}n(n1){a}_{n}{x}^{n2}\sum _{n=0}^{\mathrm{\infty}}{\alpha}_{3}{\omega}^{2}{a}_{n}{x}^{n}\sum _{n=0}^{\mathrm{\infty}}{\alpha}_{4}{\omega}^{2}{a}_{n}{x}^{n+1}=0.$
To obtain the solution, set each power coefficient of $x$ equal to zero, thus:
$\dots ,$
${a}_{r+4}=\frac{{I}_{r+4}}{{\beta}_{4}\left(r+4\right)\left(r+3\right)\left(r+2\right)\left(r+1\right)},$
${I}_{r+4}={\beta}_{3}(r+3)(r+2{)}^{2}(r+1){a}_{r+3}{\beta}_{2}(r+2{)}^{2}(r+1{)}^{2}{a}_{r+2}$
$\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}\mathrm{}{\beta}_{1}(r+2)(r+1{)}^{2}r{a}_{r+1}+{\alpha}_{3}{\omega}^{2}{a}_{r}+{\alpha}_{4}{\omega}^{2}{a}_{r1},$
$r\ge 1.$
${a}_{0}$, ${a}_{1}$, ${a}_{2}$ and ${a}_{3}$ are different constant. Some boundary conditions are listed in order to discuss the problem further.
Assuming that the beam is stationary initially:
Substitute Eq. (27) into Eq. (21), yields $\phi =0$. Suppose the one end of the beam is fixed while the other is free:
The main modes for the systems are obtained:
The frequency equation is:
Thus:
4. Results and discussion
To compare with the results obtained by methods mentioned above, finite element simulations using ABAQUS [26] are employed to calculate both the firstorder longitudinal and transverse natural modes of beams in six cases. Data in Table 1 reveals the relative difference between the analytical results and the finite element results. In the analytical procedure, a seventhorder polynomial is adopted in the equations.
Table 1. Comparison between ABAQUS results and ones obtained by analytical methods
${a}_{10}={a}_{20}=$ 0.05 m

Longitudinal frequency

Transverse frequency


Analysis (rad/s)

Abaqus (rad/s)

Difference

Analysis (rad/s)

Abaqus (rad/s)

Difference


${\rho}_{1}={\rho}_{2}=$ 7.8×10^{3} kg/m^{3}${E}_{1}={E}_{2}=$ 200×10^{9} Pa

8063.5

7953.7

1.36 %

512.7

509.7

0.59 %

${\rho}_{1}={\rho}_{2}=$ 7.8×10^{3}^{}kg/m^{3}
${E}_{1}=$200×10^{9} Pa
${E}_{2}=$ 150×10^{9} Pa

7542.7

7439.7

1.37 %

480.3

473.0

1.52 %

${\rho}_{1}={\rho}_{2}=$ 7.8×10^{3}^{}kg/m^{3}
${E}_{1}=$200×10^{9} Pa
${E}_{2}=$ 100×10^{9} Pa

6983.2

6886.2

1.39 %

436.0

422.4

3.12 %

${\rho}_{1}={\rho}_{2}=$ 8.0×10^{3}^{}kg/m^{3}
${E}_{1}={E}_{2}=$ 200×10^{9} Pa

7962.3

7853.7

1.36 %

506.2

503.2

0.59 %

${\rho}_{1}=$8.0×10^{3}^{}kg/m^{3}
${\rho}_{2}=$ 6.0×10^{3}^{}kg/m^{3}
${E}_{1}={E}_{2}=$ 200×10^{9} Pa

8512.0

8395.5

1.36 %

541.2

540.0

0.22 %

${\rho}_{1}=$8.0×10^{3}^{}kg/m^{3}
${\rho}_{2}=$ 4.0×10^{3}^{}kg/m^{3}
${E}_{1}={E}_{2}=$ 200×10^{9} Pa

9194.0

9066.4

1.38 %

584.6

581.0

0.62 %

4.1. The longitudinal frequency
Table 2 reveals the effect of shape factor ${a}_{10}/{a}_{20}$ on the frequency, as it can be seen from the charts, the boundary conditions imposes a strong effect on the shape factor. If both beam ends are fixed (or free), the longitudinal frequencies of the composite beam have little difference when the shape factors ${a}_{10}/{a}_{20}$ are reciprocal. Only with mixed boundary conditions, can frequencies vary according to the shape factor ${a}_{10}/{a}_{20}$. In addition, with the same material characteristics and shape size, the boundary conditions of both ends fixed and free share no differences.
In Fig. 2(a), the firstorder natural frequency of composite beams is examined. The elastic modulus of the upper beam is fixed at ${E}_{1}=\mathrm{}$200×10^{9} Pa, while the elastic modulus of the lower beam is decreased (that is, ${E}_{1}{E}_{2}$ rises from 0 to 100 (×10^{9} Pa)), which leads to the declination of the firstorder natural frequency of composite beams. For comparison, if a reference beam is made by a uniform material whose elastic modulus is the average of the upper and the lower beam, a set of the firstorder natural frequency of the uniform beam can be obtained (see the blue curve of reference results in Fig. 2(a)). From Fig. 2(a), the firstorder natural frequency of the composite beam is considerably different than that of the reference uniform beam, that is, the difference between blue and black curve is especially significant. The black one composited of two bonded uniform beams, while the blue one is a uniform beam whose elastic modulus is the average of the two prior components. This difference implies nonlinear coupling between the longitudinal modes of the two component beams. The shape gradient also plays an important role in the natural frequency. In this case when the boundary condition is selected as fixed in the left and free in the right, the frequency is higher as ${a}_{10}/{a}_{20}$ rises.
Table 2. Effects of the shape factor on the longitudinal frequency with different boundary conditions
Fixed in both left and right

Free in both left and right

Fixed in left and free in right


${\rho}_{1}={\rho}_{2}=$ 7.8×10^{3 }kg/m^{3}
${E}_{1}=$200×10^{9} Pa
${E}_{2}=$ 100×10^{9} Pa

${a}_{10}/{a}_{20}=\text{1}/\text{3}$

13710.0

13725.0

6630.6

${a}_{10}/{a}_{20}=\text{1}/\text{1}$

13761.0

13758.0

6886.2


${a}_{10}/{a}_{20}=\text{3}/\text{1}$

13710.0

13725.0

7096.8


${\rho}_{1}=$8.0×10^{3 }kg/m^{3}
${\rho}_{2}=$ 4.0×10^{3 }kg/m^{3}
${E}_{1}={E}_{2}=$ 200×10^{9} Pa

${a}_{10}/{a}_{20}=\text{1}/\text{3}$

18187.0

18187.0

8873.6

${a}_{10}/{a}_{20}=\text{1}/\text{1}$

18119.0

18115.0

9066.4


${a}_{10}/{a}_{20}=\text{3}/\text{1}$

18111.0

18188.0

9386.3

In Fig. 2(b), the firstorder natural frequency of composite beam is plotted as a function of the densities difference of two components. As expected, the firstorder natural frequency of composite beams rises with the reduction of density of lower beam. The blue curve of reference solution in Fig. 2(b) is that of a uniform beam whose density is the average between upper and lower beam. Obviously the reference results have a significant difference not only from the composite beam with constant crosssection but also from varying crosssection composite beam. And the firstorder natural frequency is greater when the size gradient gets bigger.
Fig. 2. The firstorder longitudinal frequency of beams with different crosssectional areas; ${a}_{10}+{a}_{20}=a=\mathrm{}$0.1 m, $l=$ 1 m
a)${E}_{1}=\mathrm{}$200×10^{9} Pa, ${\rho}_{1}={\rho}_{2}=\mathrm{}$7.8×10^{3} kg/m^{3}
b)${E}_{1}={E}_{2}=\mathrm{}$200×10^{9} Pa, ${\rho}_{1}=\mathrm{}$8.0×10^{3} kg/m^{3}
4.2. The transverse frequency and mode of vibration
Table 3 reveals the effects of the shape factor on transverse frequency with four representative boundary conditions combined by fixed, simple supported, and free conditions.
If the composite beam conditions at both ends are the same, the transverse frequency is the same when the shape factors ${a}_{10}/{a}_{20}$ are reciprocal, owing to the mirror symmetry. When the left end is fixed, and the right end sets free, the transverse frequency of composite beam varies obviously according to the shape factor ${a}_{10}/{a}_{20}$. Therefore, the impacts of physical characteristic and shape properties on the transverse frequency with the leftfixed and rightfree boundary conditions are discussed on the following section.
Fig. 3 shows the effect of physical properties and constituent ratio on the transverse frequency of the composite beam. For a composite beam with left end fixed and right end free, if one wants to improve the transverse frequency, a lower beam with bigger modulus or smaller density can be selected without changing the upper beam. Alternatively, one can tune the transverse frequency by changing the shape factor (constituent ratio along the axial direction) of the composite beam. The bigger the shape factor ${a}_{10}/{a}_{20}$ is, the greater the transverse frequencies of the composite beam are.
Table 3. Effects of the shape factor on the transverse frequency with different boundary conditions
Left and right fixed

Left (fixed) right (simplesupported)

Left and right freely supported

Left (fixed) right (free)


${\rho}_{1}={\rho}_{2}=$ 7.8×10^{3 }kg/m^{3}
${E}_{1}=$200×10^{9} Pa
${E}_{2}=$ 100×10^{9} Pa

${a}_{10}/{a}_{20}=\text{1}/\text{3}$

2582.9

1995.6

1646.5

419.8

${a}_{10}/{a}_{20}=\text{1}/\text{1}$

3521.3

2679.8

2191.4

551.0


${a}_{10}/{a}_{20}=\text{3}/\text{1}$

2566.4

2010.9

1650.3

421.1


${\rho}_{1}=$8.0×10^{3 }kg/m^{3}
${\rho}_{2}=$ 4.0×10^{3 }kg/m^{3}
${E}_{1}={E}_{2}=$ 200×10^{9} Pa

${a}_{10}/{a}_{20}=\text{1}/\text{3}$

3527.6

2713.5

2195.0

579.3

${a}_{10}/{a}_{20}=\text{1}/\text{1}$

2582.9

2050.0

1646.5

428.6


${a}_{10}/{a}_{20}=\text{3}/\text{1}$

3521.3

2740.6

2191.4

610.1

Fig. 3. The firstorder transverse frequency with ${a}_{10}+{a}_{20}=a=\mathrm{}$0.1 m, $l=$ 1 m
Fig. 4. The mode of vibration of beams with different crosssectional areas; ${E}_{1}=\mathrm{}$200×10^{9} Pa, ${E}_{2}=\mathrm{}$1×10^{9} Pa, ${\rho}_{1}={\rho}_{2}=\mathrm{}$7.8×10^{3} kg/m^{3}, ${a}_{10}+{a}_{20}=a=\mathrm{}$0.1 m, $l=\mathrm{}$1 m
Comparing with a reference beam whose physical characteristics are the average of the ones of upper and lower beams, when changing the modulus of composite beam, the reference solutions of the single beam can be observed in Fig. 3. As the modulus of the lower beam varies, the reference results are close to the transverse frequency of composite beam with a shape factor of ${a}_{10}/{a}_{20}=\text{1}/\text{1}$, while as the density of the lower beam varies, the reference solution is close to the frequency of the composite beam with a shape factor of ${a}_{10}/{a}_{20}=\text{1}/\text{1}$. Fig. 4 reveals the effect of shape factor on the first and second order lateral modes of composite beam with distinct material characteristics. The left boundary condition is fixed, while the right is free. In this case, the modulus of upper beam is taken as ${E}_{1}=\mathrm{}$200×10^{9} Pa, while the lower is ${E}_{2}=\mathrm{}$1×10^{9} Pa, mimicking metalpolymer composite. For the first mode, the shape factor ${a}_{10}/{a}_{20}$ effect is small; however, for the second mode, the decrease of shape factor ${a}_{10}/{a}_{20}$ can induce the composite beam to bend more severely. At this time, the deformation of the composite beam becomes more exaggerated with a larger curvature.
Authors Botong Li and Longlei Dong made equal contribution.
5. Conclusions
Using the variational method the kinetic equations of composite beams are obtained, and analytical expressions on natural frequencies of longitudinal and transverse vibration of composite beams are deduced by using the method of power series. Finite element simulation is used to verify the theoretical framework. Effects of gradient ratio between components along the axial direction on the natural frequencies of composite beams are revealed, and the first and second order transverse modes of composite beams are drawn. Some of the important findings of the paper are:
1) When the elastic modulus of lower beam is decreased, both the longitudinal and transverse firstorder natural frequencies of composite beams tend to decline, while they rise with the reducing of density of lower beam.
2) The shape gradient plays an important role in the natural frequency. For example, in a case where the left end is fixed and the right end is free, the firstorder longitudinal frequency is higher as the shape factor ${a}_{10}/{a}_{20}$ rises. And the effect of shape factor ${a}_{10}/{a}_{20}$ varies according to the boundary conditions. The coupling between the two bonded constituent beams is nonlinear.
3) With distinct material characteristics, the shape factor ${a}_{10}/{a}_{20}$ effect on the first order lateral mode is small; however，for the second mode, the decrease of shape factor ${a}_{10}/{a}_{20}$ can induce the composite beam to bend more severely.
Through these findings, we emphasize that both the natural frequency and vibration mode are controllable in a significantly wide range. Though the volume of each component of the composite beam is equal, by adjusting the elastic modulus or shape factor, the natural frequency or vibration mode could be changed to avoid or to achieve a specific vibration performance. For example, the natural frequency of a bridge should avoid some specific range in case of resonance.
Acknowledgements
The work of X. C. is supported by the National Natural Science Foundation of China (11172231 and 11372241), AFOSR (FA95501210159), and ARPAE (DEAR0000396). The work of B. L. is supported by the National Natural Science Foundation of China (11402188), the Fundamental Research Funds for the Central Universities (08143047) (2014gjhz16), and the Natural Science Foundation of Shaanxi (2015JQ1018).
References
 Christensen Richard M. Mechanics of Composite Materials. John Wiley and Sons, New York, 1979. [Search CrossRef]
 Jones Robert M. Mechanics of Composite Materials. Hemisphere, New York, 1975. [Search CrossRef]
 Nakamura S.I., Momiyama Y., Hosaka T., Homma K. New technologies of steel/concrete composite bridges. Journal of Constructional Steel Research, Vol. 58, 2002, p. 99130. [Search CrossRef]
 Brozzetti J. Design development of steelconcrete composite bridges in France. Journal of Constructional Steel Research, Vol. 55, 2000, p. 229243. [Search CrossRef]
 JafariTalookolaei Ramazan A., Maryam Abedi, Kargarnovin Mohammad H., Ahmadian Mohammad T. An analytical approach for the free vibration analysis of generally laminated composite beams with shear effect and rotary inertia. International Journal of Mechanical Sciences, Vol. 65, 2012, p. 97104. [Search CrossRef]
 Tseng Y. P., Huang C. S., Kao M. S. Inplane vibration of laminated curved beams with variable curvature by dynamic stiffness analysis. Composite Structures, Vol. 50, 2000, p. 103114. [Search CrossRef]
 Banerjee J. R. Frequency equation and mode shape formulae for composite Timoshenko beams. Composite Structures, Vol. 51, 2001, p. 381388. [Search CrossRef]
 Rao Kamesware M., Desai Y. M., Chitnis M. R. Free vibrations of laminated beams using mixed theory. Composite Structures, Vol. 52, 2001, p. 149160. [Search CrossRef]
 Chen W. Q., Lv C. F., Bian Z. G. Elasticity solution for free vibration of laminated beams. Composite Structures, Vol. 62, 2003, p. 7582. [Search CrossRef]
 Hajianmaleki Mehdi, Qatu Mohamad S. Vibrations of straight and curved composite beams: a review. Composite Structures, Vol. 100, 2013, p. 218232. [Search CrossRef]
 Zhen Wu, Chen Wanji An assessment of several displacementbased theories for the vibration and stability analysis of laminated composite and sandwich beams. Composite Structures, Vol. 84, 2008, p. 337349. [Search CrossRef]
 Qu Yegao, Long Xinhua, Li Hongguang, Meng Guang A variational formulation for dynamic analysis of composite laminated beams based on a general higherorder shear deformation theory. Composite Structures, Vol. 102, 2013, p. 175192. [Search CrossRef]
 Gunda Jagadish Babu, Gupta R. K., Janardhan Ranga G., Rao Venkateswara G. Large amplitude vibration analysis of composite beams: simple closedform solutions. Composite Structures, Vol. 93, Issue 2, 2011, p. 870879. [Search CrossRef]
 Yaning Li, Christine Ortiz, Mary C. Boyce A generalized mechanical model for suture interfaces of arbitrary geometry. Journal of the Mechanics and Physics of Solids, Vol. 61, Issue 4, 2013, p. 11441167. [Search CrossRef]
 Adams R. D., Wake W. C. Structural Adhesive Joints in Engineering. Elsevier, New York, 1984. [Search CrossRef]
 Zeng Q.G., Sun C. T. Novel design of a bonded lap joint. AIAA Journal, Vol. 39, Issue 10, 2001, p. 19911996. [Search CrossRef]
 Haghpanah B., Chiu S., Vaziri A. Adhesively bonded lap joints with extreme interface geometry. International Journal of Adhesion and Adhesives, Vol. 48, 2014, p. 130138. [Search CrossRef]
 Nakano H., Omiya Y., Sekiguchi Y., et al. Threedimensional FEM stress analysis and strength prediction of scarf adhesive joints with similar adherents subjected to static tensile loadings. International Journal of Adhesion and Adhesives, Vol. 54, 2014, p. 4050. [Search CrossRef]
 Kimiaeifar A., Toft H., Lund E., et al. Reliability analysis of adhesive bonded scarf joints. Engineering Structures, Vol. 35, 2012, p. 281287. [Search CrossRef]
 Lubkin J. L. A theory of adhesive scarf joints. Journal of Applied Mechanics, Vol. 24, 1957, p. 255260. [Search CrossRef]
 Atlihan Gökmen, Çallioğlu Hasan, Şahin Conkur E., Topcu Muzaffer, Yücel Uğur Free vibration analysis of the laminated composite beams by using DQM. Journal of Reinforced Plastics and Composites, Vol. 28, Issue 7, 2009, p. 881892. [Search CrossRef]
 Malekzadeh P., Vosoughi A. R. DQM large amplitude vibration of composite beams on nonlinear elastic foundations with restrained edges. Communications in Nonlinear Science and Numerical Simulation, Vol. 14, 2009, p. 906915. [Search CrossRef]
 Damanpack A. R., Khalili S. M. R. Highorder free vibration analysis of sandwich beams with a flexible core using dynamic stiffness method. Composite Structures, Vol. 94, Issue 5, 2012, p. 15031514. [Search CrossRef]
 Vidal P., Polit O. A family of sinus finite elements for the analysis of rectangular laminated beams. Composite Structures, Vol. 84, 2008, p. 5672. [Search CrossRef]
 Özütok A., Madenci E. Free vibration analysis of crossply laminated composite beams by mixed finite element formulation. International Journal of Structural Stability and Dynamics, Vol. 13, Issue 2, 2013, p. 1250056. [Search CrossRef]
 ABAQUS 6.11 User’s Manual. ABAQUS Inc., Providence, RI, 2011. [Search CrossRef]