Study on natural characteristics of fiber metal laminates thin plates under cantilever boundary
Yunshan Liu^{1} , Feng Shang^{2} , Zhuo Xu^{3} , Bangchun Wen^{4}
^{1}Guidaojiaotong Polytechnic Institute, Shenyang, 110023, China
^{2}Beijing Xingyou Engineering Project Management Co., Ltd, Chengdu, Sichuan, 610017, China
^{3, 4}School of Mechanical Engineering and Automation, Northeastern University, Shenyang, 110004, China
^{1}Corresponding author
Journal of Vibroengineering, Vol. 22, Issue 4, 2020, p. 909922.
https://doi.org/10.21595/jve.2020.21383
Received 18 March 2020; received in revised form 14 May 2020; accepted 3 June 2020; published 30 June 2020
JVE Conferences
Through the combination of theory and experiment, the natural characteristics of the fiber metal laminates thin plates under cantilever boundary are analyzed and verified. Based on the mechanics of composites and classical laminated plate theory, the theoretical model is established. The orthogonal polynomial method and the energy method are used to solve the natural characteristics. Meanwhile the calculation process is proposed. And then, the natural characteristics of a TA2/TC500 fiber metal laminates thin plate are tested. It is found that comparing the calculation results of the frequencies with the test ones, the errors are within the range of 3.4 % to 4.5 %, the trends of modal shapes are consistent as well, thus the effectiveness of above method has been verified.
Keywords: FMLs thin plates, cantilever boundary, orthogonal polynomial method, natural characteristics.
1. Introduction
Fiber metal laminates (FMLs) are new types of composite materials, which are composed of fiber reinforced layers and metal layers alternately [13]. Due to FMLs contained both of the advantages of metal materials, such as good toughness, strong impact resistance, high damage tolerance etc., and composite materials, such as high specific strength and specific stiffness, corrosion resistance and fatigue resistance, they have widely used in the wing, tail and hatch of advanced aircraft [46]. Since the above composite materials and the typical structural parts such as beams, plates and shells are often used in harsh environment, it is easy to produce problems such as excessive vibration, fatigue failure and fatigue wear. Therefore, it is necessary to study the vibration characteristics of FMLs thin plates which have important engineering and academic significances in the field of theoretical analysis, dynamic design and fault diagnosis [7].
The natural characteristics, such as natural frequencies and modal shapes, are the basis of indepth study on the vibration characteristics of structural system. So far, a lot of researches on natural characteristics of FMLs thin plates had been done. Harras [8] et al. established a theoretical model of GLARE 3 fiber metal laminated thin plate under clamped boundary condition and the natural characteristics were obtained. Utilizing the experimental technology, Botelho [9] analyzed the damping characteristics of aluminum plate, carbon fiber/resin composite plate and glass fiber/aluminum alloy FMLs plate. Based on the firstorder shear deformation theory, Shooshtari and Razavi [10] deduced the nonlinear ordinary differential equations of the fiber metal laminated thin plate by employing the Galerkin method. Meanwhile, the linear and nonlinear natural characteristics of the structure were analyzed. Payeganeh [11] et al. studied the dynamic characteristics of FMLs plate under lowspeed impact, and it was found that parameters such as the lamination sequences, aspect ratios, impact speeds, and mass had significant effect on the dynamic characteristics of FMLs plate. Using Ritz method and ABAQUS finite element method, Ghasemi [12] et al. obtained the dimensionless natural frequency of fiberreinforced metal laminates under simply supported boundary condition. Then, the influences of geometric parameters and different distributions of metal layers on the vibration parameters were analyzed. Rahimi [13] et al. proposed a threedimensional elastic analysis theory of FMLs based on the state space differential quadrature method. And the dimensionless frequency of the annular FMLs plate was calculated. Employing Navier method and Hamilton principle, Mahi [14] et al. presented the energy equation of various thin plates by using Ritz method. Then the natural frequencies the corresponding plates under different boundaries were solved. Iriondo [15] et al. studied the damping characteristics of the traditional FMLs plate and the self reinforced polypropylene FMLs plate by utilizing the resonance experiment under forced vibration. Sayyad and Ghugal [16] used the triangle shear method and the theory of normal deformation to establish a dynamic model of multilayer laminates under simple boundary conditions based on the virtual work principle. Li [17] et al. studied natural frequencies of CTi FMLs beams and plates under cantilever boundary condition by employing experiments and ABAQUS finite element method. Meanwhile, the influences of different sizes and different numbers of metal layers on their natural frequencies were analyzed. However, the theoretical derivation had not been done.
It can be found from the available literatures that the researches on the natural characteristics of FMLs thin plate under cantilever boundary were not enough, and most of the literatures had not been verified by experiments. Therefore, it is necessary to conduct further research on its natural characteristics. In this paper, a theoretical model of FMLs under cantilever boundary condition is established based on composite mechanics and classic laminate plate theory. The natural characteristics of the FMLs thin plates are solved by utilizing the orthogonal polynomial method and the energy method. As an example, to demonstrate the feasibility of the proposed model, the experimental test of a TA2/TC500 laminated thin plate is implemented. The calculated and measured results show a good consistency.
2. Theoretical solution of the natural characteristics of the fiber metal laminates thin plates under the cantilever boundary
2.1. Theoretical model
FMLs thin plate is made of metal material and $n$layer orthotropic fiber and matrix composite fiber reinforced material. Its theoretical model is shown in Fig. 1. In Fig. 1, a rectangular coordinate system is established with the middle plane as the reference plane. The length of the thin plate is $a$ and the width is $b$. The sides of the $yoz$ plane are fixed. The thickness of the metal layer is ${h}_{1}$ while the thickness of the fiber layer is ${h}_{2}$. Therefore, the total thickness of the thin plate is $h=3{h}_{1}+2{h}_{2}$, the angle between the longitudinal direction of the fiber and the $x$ direction is $\theta $. Its material parameters are expressed as follows, ${E}_{1}$ is the elastic modulus of the metal, ${\upsilon}_{1}$ is the Poisson’s ratio of the metal, and ${E}_{21}$, ${E}_{22}$ and ${G}_{12}$ are the Young’s modulus of the fiber layer along the fiber direction, vertical fiber direction and inplane shear in the fiber layer respectively, ${\upsilon}_{12}$ and ${\upsilon}_{21}$ are the Poisson’s ratios along the fiber direction and the vertical fiber direction, respectively, and ${\rho}_{1}$ and ${\rho}_{2}$ represent the densities of the metal layer and the fiber layer, respectively.
Fig. 1. Theoretical model of the FMLs thin plate under cantilever boundary
Due to the model is a symmetric structure, based on Kirchhoff hypothesis and classical thin plate theory, its displacement components $u$, $v$, $w$ along the $x$, $y$, and $z$ directions can be expressed as [18]:
where ${u}_{0}$, ${v}_{0}$ and ${w}_{0}$ represent the displacement of the middle plane of the thin plate respectively, and $t$ represents the time.
The strain at any point of the thin plate can be expressed by displacement [18]:
where, ${\epsilon}_{x}$, ${\epsilon}_{y}$, ${\epsilon}_{z}$ and ${\gamma}_{yz}$, ${\gamma}_{xz}$ and ${\gamma}_{xy}$ represent the line strain and shear strain in $x$, $y$ and $z$ directions respectively.
Therefore, the stressstrain relationship should be represented as [18]:
where, for the metal layer, ${Q}_{11}={Q}_{22}=\frac{{E}_{1}}{1{\upsilon}_{1}^{2}}$, ${Q}_{12}={Q}_{21}=\frac{{\upsilon}_{1}{E}_{1}}{1{\upsilon}_{1}^{2}}$, ${Q}_{66}=\frac{{E}_{1}}{2\left(1+{\upsilon}_{1}\right)}$; for fiber layer, ${Q}_{11}=\frac{{E}_{21}}{1{\upsilon}_{12}{\upsilon}_{21}}$, ${Q}_{12}={Q}_{21}=\frac{{\upsilon}_{12}{E}_{22}}{1{\upsilon}_{12}{\upsilon}_{21}}=\frac{{\upsilon}_{21}{E}_{21}}{1{\upsilon}_{12}{\upsilon}_{21}}$, ${Q}_{22}=\frac{{E}_{22}}{1{\upsilon}_{12}{\upsilon}_{21}}$, ${Q}_{66}={G}_{12}$.
In the fiber layer, when the angle between the fiber direction and the $x$axis is $\theta $, by using the rotation axis formula, the stressstrain relationship of the $k$th layer fiberboard is obtained as follows:
where:
where $k$ represents the $k$th layer fiber and ${\theta}_{k}$ represents the angle between the $k$th layer fiber and the $x$axis direction.
2.2. Energy equations
In the bending vibration of thin plates, the kinetic energy and strain energy of the metal layer of the model can be expressed as follows:
The kinetic energy and strain energy of the fiber layer can be expressed as follows:
where $A$ is the area of the thin plate. Eq. (1) and Eq. (2) are introduced into Eq. (6), then Eq. (3) and Eq. (4) are introduced into Eq. (6a)(6d) and (6e)(6f) respectively. The expressions of kinetic energy and strain energy of metal layer and fiber layer can be obtained through the middle displacement ${w}_{0}$:
$\left.+{Q}_{22}{\left(\frac{{\partial}^{2}{w}_{0}}{\partial {y}^{2}}\right)}^{2}+4{Q}_{66}{\left(\frac{{\partial}^{2}{w}_{0}}{\partial x\partial y}\right)}^{2}\right]dzdA,$
$\left.+{Q}_{22}{\left(\frac{{\partial}^{2}{w}_{0}}{\partial {y}^{2}}\right)}^{2}+4{Q}_{66}{\left(\frac{{\partial}^{2}{w}_{0}}{\partial x\partial y}\right)}^{2}\right]dzdA,$
$+{Q}_{22}{\left(\frac{{\partial}^{2}{w}_{0}}{\partial {y}^{2}}\right)}^{2}+\left.4{Q}_{66}{\left(\frac{{\partial}^{2}{w}_{0}}{\partial x\partial y}\right)}^{2}\right]dzdA$
$+{\overline{Q}}_{22}{\left(\frac{{\partial}^{2}{w}_{0}}{\partial {y}^{2}}\right)}^{2}+\left.4{\overline{Q}}_{26}\frac{{\partial}^{2}{w}_{0}}{\partial {y}^{2}}\cdot \frac{{\partial}^{2}{w}_{0}}{\partial x\partial y}+4{\overline{Q}}_{66}{\left(\frac{{\partial}^{2}{w}_{0}}{\partial x\partial y}\right)}^{2}\right]dzdA,$
$\left.+{\overline{Q}}_{22}{\left(\frac{{\partial}^{2}{w}_{0}}{\partial {y}^{2}}\right)}^{2}+4{\overline{Q}}_{26}\frac{{\partial}^{2}{w}_{0}}{\partial {y}^{2}}\cdot \frac{{\partial}^{2}{w}_{0}}{\partial x\partial y}+4{\overline{Q}}_{66}{\left(\frac{{\partial}^{2}{w}_{0}}{\partial x\partial y}\right)}^{2}\right]dzdA.$
By adding the kinetic energy and strain energy of each layer, the total kinetic energy and strain energy can be obtained as:
2.3. Solution by orthogonal polynomial method
As the orthogonal polynomial method has clear solution principle, fast calculation speed, and can be applied to solve the natural characteristics of various boundaries, reference [14], using this method, the mid plane displacement is expressed as:
where $M$ and $N$ are called the truncation coefficients for solution, ${A}_{mn}$ is a pending parameter, $\omega $ is the angular natural frequency, ${P}_{m}\left(\xi \right)$ and ${P}_{n}\left(\eta \right)$ are orthogonal characteristic polynomials, and their expression is as follows:
where ${B}_{k}$ and ${C}_{k}$ are coefficient parameters, $\varphi \left(\xi \right)$ and $\phi \left(\eta \right)$ are polynomial functions for determining boundary conditions, and their expression is as follows:
where $W\left(\zeta \right)$ is the weighting function used in the orthogonalization process, and usually $W\left(\zeta \right)=1$; $p$, $q$, $r$, $s$ represent boundary state parameters at $x=$0, $x=a$, $y=$0, $y=b$ respectively. When the boundary state is free, the value of the corresponding boundary state parameter is 0; When the boundary state is a simply supported boundary, the value of the corresponding boundary state parameter is 1; When the boundary state is a fixed boundary, the corresponding boundary state parameter is 2. According to the model established in this paper, the boundary state parameters $p$, $q$, $r$, $s$ correspond to the boundary state. Therefore, the values of $p$, $q$, $r$, $s$ are 2, 0, 0, 0 [14].
By substituting Eq. (9) into Eq. (8), the expressions of maximum kinetic energy ${E}_{k}^{max}$ and maximum strain energy ${U}_{s}^{max}$ with undetermined parameters can be obtained by making $\mathrm{c}\mathrm{o}\mathrm{s}\left(\omega t\right)=$1 and $\mathrm{s}\mathrm{i}\mathrm{n}\left(\omega t\right)=$1, respectively:
$\left.+{\left(\sum _{m=1}^{M}\sum _{n=1}^{N}{A}_{mn}{P}_{m}{P}_{n}\right)}^{2}\right\}dzdA+\frac{1}{2}{\int}_{A}{\int}_{{h}_{1}/2}^{{h}_{1}/2}{\rho}_{1}{\omega}^{2}\left\{{z}^{2}{\left(\sum _{m=1}^{M}\sum _{n=1}^{N}{A}_{mn}\frac{d{P}_{m}}{dx}{P}_{n}\right)}^{2}\right.$
$+{z}^{2}{\left(\sum _{m=1}^{M}\sum _{n=1}^{N}{A}_{mn}{P}_{m}\frac{d{P}_{n}}{dy}\right)}^{2}\left.+{\left(\sum _{m=1}^{M}\sum _{n=1}^{N}{A}_{mn}{P}_{m}{P}_{n}\right)}^{2}\right\}dzdA$
$+\frac{1}{2}{\int}_{A}{\int}_{h/2}^{{h}_{2}{h}_{1}/2}{\rho}_{1}{\omega}^{2}\left\{{z}^{2}{\left(\sum _{m=1}^{M}\sum _{n=1}^{N}{A}_{mn}\frac{d{P}_{m}}{dx}{P}_{n}\right)}^{2}\right.+{z}^{2}{\left(\sum _{m=1}^{M}\sum _{n=1}^{N}{A}_{mn}{P}_{m}\frac{d{P}_{n}}{dy}\right)}^{2}$
$\left.+{\left(\sum _{m=1}^{M}\sum _{n=1}^{N}{A}_{mn}{P}_{m}{P}_{n}\right)}^{2}\right\}dzdA+\frac{1}{2}{\int}_{A}{\int}_{{h}_{1}/2}^{{h}_{2}+{h}_{1}/2}{\rho}_{2}{\omega}^{2}\left\{{z}^{2}{\left(\sum _{m=1}^{M}\sum _{n=1}^{N}{A}_{mn}\frac{d{P}_{m}}{dx}{P}_{n}\right)}^{2}\right.$
$+\frac{1}{2}{\int}_{A}{\int}_{{h}_{2}{h}_{1}/2}^{{h}_{1}/2}{\rho}_{2}{\omega}^{2}\left\{{z}^{2}{\left(\sum _{m=1}^{M}\sum _{n=1}^{N}{A}_{mn}\frac{d{P}_{m}}{dx}{P}_{n}\right)}^{2}\right.+{z}^{2}{\left(\sum _{m=1}^{M}\sum _{n=1}^{N}{A}_{mn}{P}_{m}\frac{d{P}_{n}}{dy}\right)}^{2}$
$\left.+{\left(\sum _{m=1}^{M}\sum _{n=1}^{N}{A}_{mn}{P}_{m}{P}_{n}\right)}^{2}\right\}dzdA,$
$+2{Q}_{12}\left(\sum _{m=1}^{M}\sum _{n=1}^{N}{A}_{mn}\frac{{d}^{2}{P}_{m}}{d{x}^{2}}{P}_{n}\right)\left(\sum _{m=1}^{M}\sum _{n=1}^{N}{A}_{mn}{P}_{m}\frac{{d}^{2}{P}_{n}}{d{y}^{2}}\right)$
$+{Q}_{22}{\left(\sum _{m=1}^{M}\sum _{n=1}^{N}{A}_{mn}{P}_{m}\frac{{d}^{2}{P}_{n}}{d{y}^{2}}\right)}^{2}\left.+4{Q}_{66}{\left(\sum _{m=1}^{M}\sum _{n=1}^{N}{A}_{mn}\frac{d{P}_{m}}{dx}\frac{d{P}_{n}}{dy}\right)}^{2}\right\}dzdA$
$+\frac{1}{2}{\int}_{A}{\int}_{{h}_{1}/2}^{{h}_{1}/2}{z}^{2}\left\{{Q}_{11}{\left(\sum _{m=1}^{M}\sum _{n=1}^{N}{A}_{mn}\frac{{d}^{2}{P}_{m}}{d{x}^{2}}{P}_{n}\right)}^{2}\right.$
$+2{Q}_{12}\left(\sum _{m=1}^{M}\sum _{n=1}^{N}{A}_{mn}\frac{{d}^{2}{P}_{m}}{d{x}^{2}}{P}_{n}\right)\left(\sum _{m=1}^{M}\sum _{n=1}^{N}{A}_{mn}{P}_{m}\frac{{d}^{2}{P}_{n}}{d{y}^{2}}\right)$
$+{Q}_{22}{\left(\sum _{m=1}^{M}\sum _{n=1}^{N}{A}_{mn}{P}_{m}\frac{{d}^{2}{P}_{n}}{d{y}^{2}}\right)}^{2}\left.+4{Q}_{66}{\left(\sum _{m=1}^{M}\sum _{n=1}^{N}{A}_{mn}\frac{d{P}_{m}}{dx}\frac{d{P}_{n}}{dy}\right)}^{2}\right\}dzdA$
$+\frac{1}{2}{\int}_{A}{\int}_{h/2}^{{h}_{2}{h}_{1}/2}{z}^{2}\left\{{Q}_{11}{\left(\sum _{m=1}^{M}\sum _{n=1}^{N}{A}_{mn}\frac{{d}^{2}{P}_{m}}{d{x}^{2}}{P}_{n}\right)}^{2}\right.$
$+2{Q}_{12}\left(\sum _{m=1}^{M}\sum _{n=1}^{N}{A}_{mn}\frac{{d}^{2}{P}_{m}}{d{x}^{2}}{P}_{n}\right)\left(\sum _{m=1}^{M}\sum _{n=1}^{N}{A}_{mn}{P}_{m}\frac{{d}^{2}{P}_{n}}{d{y}^{2}}\right)$
$+{Q}_{22}{\left(\sum _{m=1}^{M}\sum _{n=1}^{N}{A}_{mn}{P}_{m}\frac{{d}^{2}{P}_{n}}{d{y}^{2}}\right)}^{2}\left.+4{Q}_{66}{\left(\sum _{m=1}^{M}\sum _{n=1}^{N}{A}_{mn}\frac{d{P}_{m}}{dx}\frac{d{P}_{n}}{dy}\right)}^{2}\right\}dzdA$
$+\frac{1}{2}{\int}_{A}{\int}_{{h}_{1}/2}^{{h}_{2}+{h}_{1}/2}{z}^{2}\left\{{\overline{Q}}_{11}{\left(\sum _{m=1}^{M}\sum _{n=1}^{N}{A}_{mn}\frac{{d}^{2}{P}_{m}}{d{x}^{2}}{P}_{n}\right)}^{2}\right.$
$+2{\overline{Q}}_{12}\left(\sum _{m=1}^{M}\sum _{n=1}^{N}{A}_{mn}\frac{{d}^{2}{P}_{m}}{d{x}^{2}}{P}_{n}\right)\left(\sum _{m=1}^{M}\sum _{n=1}^{N}{A}_{mn}{P}_{m}\frac{{d}^{2}{P}_{n}}{d{y}^{2}}\right)$
$+4{\stackrel{}{Q}}_{16}\left(\sum _{m=1}^{M}\sum _{n=1}^{N}{A}_{mn}\frac{{d}^{2}{P}_{m}}{d{x}^{2}}{P}_{n}\right)\left(\sum _{m=1}^{M}\sum _{n=1}^{N}{A}_{mn}\frac{d{P}_{m}}{dx}\frac{d{P}_{n}}{dy}\right)$
$+4{\stackrel{}{Q}}_{26}\left(\sum _{m=1}^{M}\sum _{n=1}^{N}{A}_{mn}{P}_{m}\frac{{d}^{2}{P}_{n}}{d{y}^{2}}\right)\left(\sum _{m=1}^{M}\sum _{n=1}^{N}{A}_{mn}\frac{d{P}_{m}}{dx}\frac{d{P}_{n}}{dy}\right)$
$+\frac{1}{2}{\int}_{A}{\int}_{{h}_{2}{h}_{1}/2}^{{h}_{1}/2}{z}^{2}\left\{{\overline{Q}}_{11}{\left(\sum _{m=1}^{M}\sum _{n=1}^{N}{A}_{mn}\frac{{d}^{2}{P}_{m}}{d{x}^{2}}{P}_{n}\right)}^{2}\right.$
$+2{\overline{Q}}_{12}\left(\sum _{m=1}^{M}\sum _{n=1}^{N}{A}_{mn}\frac{{d}^{2}{P}_{m}}{d{x}^{2}}{P}_{n}\right)\left(\sum _{m=1}^{M}\sum _{n=1}^{N}{A}_{mn}{P}_{m}\frac{{d}^{2}{P}_{n}}{d{y}^{2}}\right)$
$+4{\stackrel{}{Q}}_{16}\left(\sum _{m=1}^{M}\sum _{n=1}^{N}{A}_{mn}\frac{{d}^{2}{P}_{m}}{d{x}^{2}}{P}_{n}\right)\left(\sum _{m=1}^{M}\sum _{n=1}^{N}{A}_{mn}\frac{d{P}_{m}}{dx}\frac{d{P}_{n}}{dy}\right)$
$+4{\stackrel{}{Q}}_{26}\left(\sum _{m=1}^{M}\sum _{n=1}^{N}{A}_{mn}{P}_{m}\frac{{d}^{2}{P}_{n}}{d{y}^{2}}\right)\left(\sum _{m=1}^{M}\sum _{n=1}^{N}{A}_{mn}\frac{d{P}_{m}}{dx}\frac{d{P}_{n}}{dy}\right)$
$+{\stackrel{}{Q}}_{22}{\left(\sum _{m=1}^{M}\sum _{n=1}^{N}{A}_{mn}{P}_{m}\frac{{d}^{2}{P}_{n}}{d{y}^{2}}\right)}^{2}\left.+4{\stackrel{}{Q}}_{66}{\left(\sum _{m=1}^{M}\sum _{n=1}^{N}{A}_{mn}\frac{d{P}_{m}}{dx}\frac{d{P}_{n}}{dy}\right)}^{2}\right\}dzdA.$
So the energy equation can be expressed as:
Take the partial derivative of the energy equation $F$ with respect to the undetermined coefficients ${A}_{mn}$. Then, set it equal to 0, thus:
Solving the generalized eigenvalue problem of the minimum pending parameters of the Eq. (13) can be obtained:
where $K$ and $M$ are the stiffness matrix and mass matrix of the structural system, respectively, and the feature vector $q={\left[{A}_{11},{A}_{12},\dots {A}_{mn}\right]}^{T}$. If Eq. (15) is solved, the determinant of the coefficient matrix of the feature vector $q$ is 0, that is:
The natural frequencies of the structure can be obtained by solving Eq. (16), the accuracy of which is determined by the values of $M$ and $N$. the larger the values of $M$ and $N$, the higher the accuracy of the results. In this study, setting $M=$8 and $N=$8 can satisfy the calculation accuracy [19].
Taking the obtained natural frequencies back into Eq. (15), the undetermined parameters of ${A}_{mn}$ can be obtained, and taking them into Eq. (9), the modal shapes of each natural frequency can be obtained.
3. Analysis process of natural characteristics of fiber metal laminates thin plate
In this section, based on the MATLAB software, the corresponding program is written, and the analysis process of the natural characteristics of the FMLs thin plate under the cantilever boundary condition is proposed. The flowchart is given in Fig. 2, and corresponding specific steps are as follows:
1) Input the geometrical parameters and material parameters of FMLs thin plate.
Firstly, the length, width, thickness of metal layer and fiber layer, fiber angle of each layer and other geometric parameters of FMLs thin plate are given. Then, the elastic modulus of metal and fiber, Poisson’s ratio, material density, metal elastic modulus, elastic modulus and shear modulus of fiber along the fiber direction and perpendicular to the fiber direction are input respectively to prepare for calculation.
2) Obtain the expression of maximum kinetic energy and maximum strain energy.
The displacement component expression Eq. (1) is substituted into the general kinetic energy expression Eq. (8a) and the general strain energy expression Eq. (8b), and the later generations are arranged into the middle plane displacement expression Eq. (14) of the orthogonal polynomial method. According to different boundaries, ${P}_{m}\left(\xi \right)$ and ${P}_{n}\left(\eta \right)$ are determined. After finishing, set $\mathrm{c}\mathrm{o}\mathrm{s}\left(\omega t\right)=1$ to obtain the maximum kinetic energy ${E}_{k}^{max}$, and set $\mathrm{s}\mathrm{i}\mathrm{n}\left(\omega t\right)=$1 to obtain the maximum strain energy ${U}_{s}^{max}$.
3) Obtain the natural frequencies.
Take the maximum kinetic energy and maximum strain energy obtained in step 2 into Eq. (13), and get the specific expression of energy equation $F$. By solving the generalized eigenvalue problem of the minimum undetermined parameter, i.e. let $\mathrm{d}\mathrm{e}\mathrm{t}(K{\omega}^{2}M)=$ 0, the natural frequencies $\omega $ of FMLs can be obtained.
4) Obtain the modal shapes.
Taking the obtained natural frequencies $\omega $ into Eq. (15), the undetermined parameters ${A}_{mn}$ can be obtained, and the obtained parameters can be taken back to the displacement expression Eq. (9) to obtain the modal shape of each order.
Fig. 2. The flowchart of solving the natural characteristics of FMLs plate
4. A case study
Taking a TA2/TC500 FMLs thin plate as the research object, the natural frequencies and modal shapes of the thin plates are tested under cantilever boundary condition. The length, width and height of the object are 200×300×2.65 mm. Among them, the thickness of the metal layer is 0.3 mm, the density is 4.15×10^{3} kg/m^{3}, the elastic modulus is 108 GPa, and the Poisson’s ratio is 0.3; The total thickness of a fiber layer is 0.875 mm, the elastic model of the fiber along the fiber direction is 136 GPa, the elastic modulus of the vertical fiber direction is 7.92 GPa, the shear modulus is 3.39 GPa, Poisson’s ratio 0.3, density 1780 kg/m^{3}, the laminate configuration is [(0°/90°)_{s}/0°/(90°/0°)_{s}]. 14 layers are laid in total, and each layer has the same thickness and fiber volume fraction. The short side is clamped during installation, and the clamping length is 30 mm.
Fig. 3. The natural characteristics experiment system of FMLs thin plates
In order to verify the correctness of the calculation method proposed in this paper, a test system is built as shown in Fig. 3. The experiment adopts the method of multipoint excitation and singlepoint response, which is carried out by a modal hammer. Before the test, the tested plate is equally divided into 9 parts in the length direction and 10 parts in the width direction. The exciting points are numbered from 1 to 110. Furthermore, the response point is set at point 33. Meanwhile, the test bandwidth is set as 1600 Hz, the number of spectrum lines is 4096, and the frequency resolution is 0.39 Hz. In order to improve the test accuracy, a force index window function is also added to the excitation signal, and an exponential window function is added to the response signal.
During the test, PCB086C01 modal hammer is used to excit the FMLs plate three times at each exciting point, and LMS mobile frontend is used to collect the excitation signal and response signal. Finally, the first 7 natural frequencies and modal shapes can be obtained through the storage and analysis of notebook computer equipped with LMS Test. lab 14A analysis software. The obtained natural frequencies and modal shapes are shown in Table 1. For comparison purposes, the frequencies and modal shapes calculated by MATLAB software and corresponding errors are listed in Table 1.
By comparing the calculation results with the experimental ones, it can be found that the errors between the calculation results of the natural frequency of the FMLs thin plate and the experimental results are 3.4 % and 4.5 %, which is within an allowable range. Meanwhile, the corresponding modal shapes are completely consistent. Thus, the correctness of the theoretical analysis method can be verified.
5. Parameters studies
The theoretical model has been verified in section 4. Therefore, in this section, the relationship between the material and geometric parameters of FMLs plate and its frequencies will be discussed. For the convenience of calculation, the thickness of the fiber layer is set to be 0.1 mm, and the laminate configuration is [0°/90°]_{s}.
5.1. The influence of material parameters on the natural frequency of FMLs plate
In this section, in order to compare the relationship between different material parameters and the natural frequencies of the structures studied, the densities and elastic moduli of aluminum, titanium and stainless steel materials are respectively selected for calculation and discussion. The corresponding material parameters are shown in Table 2. The length × width × thickness of the plate is set to be 300 × 200 × 2.9 mm. The thickness of the matal layer is 0.3mm, therefore, the numbers of the fiber in the fiber layers are 10. Meanwhile, the calculated results of the first 7 orders of the natural frequencies of the FMLs plates with different material parameters are shown in Table 3.
Table 1. The first 7 natural frequencies and modal shapes of FMLs thin plate by calculation and experiment
Order

Calculation frequency (Hz) $A$

Test
frequency
(Hz) $B$

Error (%)
$\left\frac{AB}{B}\right$

Calculation
vibration mode

Test
vibration mode

1

34.4

35.7

3.8



2

71.9

75.1

4.5



3

215.6

222.7

3.3



4

296.9

309

4.1



5

475.0

491.2

3.4



6

603.3

624.1

3.4



7

699.6

726.9

3.9



Table 2. The material parameters of the FMLs plate
Type

Material

Elasticity modulus (GPa)

Density (kg/m^{3})

Case 1

Aluminum

72

2700

Case 2

Titanium

108

4150

Case 3

Stainless steel

180

7800

Table 3. The calculated frequencies of FMLs plate with different material parameters
Order

Natural frequencies (Hz)


Case 1

Case 2

Case 3


1

32.2

32.5

31.8

2

65.7

71.7

95.1

3

195.0

202.9

198.6

4

305.9

322.3

329.7

5

462.5

475.4

474.3

6

565.6

572.7

572.4

7

667.3

693.9

697.0

It can be found from the above results that the natural frequency of the structure is the smallest when the material parameters of aluminum are selected. Furthermore, when the material parameters of titanium and stainless steel are chosen, the difference of the natural frequencies of FMLs plate is very small. It may be caused by the combined effect of the modulus of elasticity and the density.
5.2. The influence of geometric parameters on the natural frequency of FMLs plate
In this section, the influence of different geometric structures on the natural frequency of FMLs plate is discussed by setting different metal layer thickness and fiber layer thickness. The elasticity modulus and density are set to be the same as them in Section 4. And the length × width × thickness of the plate is set to be 300 × 200 × 3.0 mm. Besides, the corresponding geometric structure parameters are shown in Table 4. The calculation results of natural frequencies for different geometric parameters are shown in Table 5.
Table 4. The geometric construction of the FMLs plate
Type

Thickness of the metal layers (mm)

Thickness of the fiber layers (mm)

Numbers of the fiber layers

Case 4

0

3.0

1

Case 5

0.2

1.2

2

Case 6

0.4

0.9

2

Table 5. The calculated frequencies of FMLs plate with different geometric construction
Order

Natural frequencies (Hz)


Case 4

Case 5

Case 6


1

35.8

34.5

32.9

2

77.4

68.8

62.5

3

224.4

215.5

205.3

4

338.1

320.4

308.4

5

493.0

487.7

465.7

6

605.2

581.7

571.2

7

717.3

701.4

688.6

It can be seen from the above results that with the increase of the thickness of metal layer, the natural frequencies of FMLs plate decrease gradually. The reason for this phenomenon may be that, when the total thickness is certain, the density of the metal layer is higher than that of the fiber layer. With the increase of the thickness of the metal layer, the thickness of the fiber layers gradually decreases, which leads to the increase of the overall mass and then reduces the natural frequency of the structure.
6. Conclusions
In this paper, the natural characteristics of the FMLs thin plate under cantilever boundary condition are analyzed and verified by combining theory with experiment. By comparing the calculation results with the experimental ones of the TA2/TC500 thin plate, it can be seen that the errors of the natural frequencies between the calculated results and the experimental ones are 3.4 % and 4.5 %, which are within an allowable error range. Besides, the corresponding modal shapes are completely consistent, which further verifies the correctness of the theoretical analysis method. By using the method proposed in this paper, the analysis and prediction of the natural characteristics of FMLs thin plate under the cantilever boundary condition can be realized. Meanwhile, by studying the relationship between the parameters, such as material and the geometric ones, and the natural frequencies in the structure, it can be found that the natural frequency of the structure is the smallest when the material parameters of aluminum are selected. And the natural frequencies of FMLs plate decrease gradually with the increase of the thickness of metal layer. Furthermore, the theoretical model is also applicable to the analysis of the natural characteristics of FMLs laminates thin plates under different boundaries conditions such as free, simply supported and clamped. However, the results of this paper have not analyzed and discussed the natural characteristics of FMLs plates with asymmetric and multilayers. And the corresponding results need to be further demonstrated.
Acknowledgements
This study was supported by the National Natural Science Foundation of China 51505070.
References
 Asundi A., Choi A Y N. Fiber metal laminates: An advanced material for future aircraft. Journal of Materials Processing Technology, Vol. 63, Issues 13, 1997, p. 384394. [Publisher]
 Sinmazçelik T., Avcu E., Bora M. Ö., et al. A review: Fibre metal laminates, background, bonding types and applied test methods. Materials and Design, Vol. 32, Issue 7, 2011, p. 36713685. [Publisher]
 Vlot A., Vogelesang L. B., Vries T. J. D. Towards application of fibre metal laminates in large aircraft. Aircraft Engineering and Aerospace Technology, Vol. 71, Issue 6, 1999, p. 558570. [Publisher]
 Marsh G. Airframers exploit composites in battle for supremacy. Reinforced Plastics, Vol. 49, Issue 3, 2005, p. 2632. [Publisher]
 Vogelesang L. B., Vlot A. Development of fibre metal laminates for advanced aerospace structures. Journal of Materials Processing Technology, Vol. 103, Issue 1, 2000, p. 15. [Publisher]
 Wang Y. G., Liang X. Z. Fibermetal laminate technology and large aircraft. Academic Annual Meeting of FRP/Composite Materials, 2010, (in Chinese). [CrossRef]
 Weaver W., Timoshenko W., Young S. P., et al. Vibration Problems in Engineering. Constable and Co. Ltd, 1937. [CrossRef]
 Harras B., Benamar R., White R. G. Experimental and theoretical investigation of the linear and nonlinear dynamic behavior of a glare 3 hybrid composite panel. Journal of Sound and Vibration, Vol. 251, Issue 2, 2002, p. 281315. [Publisher]
 Botelho E. C., Campos A. N., Barros E. D., et al. Damping behavior of continuous fiber/metal composite materials by the free vibration method. Composites Part B, Vol. 37, Issue 2, 2005, p. 255263. [Publisher]
 Shooshtari A., Razavi S. A closed form solution for linear and nonlinear free vibrations of composite and fiber metal laminated rectangular plates. Composite Structures, Vol. 92, Issue 11, 2010, p. 26632675. [Publisher]
 Payeganeh G. H., Ghasemi F. A., Malekzadeh K. Dynamic response of fiber–metal laminates (FMLs) subjected to lowvelocity impact. ThinWalled Structures, Vol. 48, Issue 1, 2010, p. 6270. [Publisher]
 Ghasemi F. A., Paknejad R., Fard K. M. Effects of geometrical and material parameters on free vibration analysis of fiber metal laminated plates. Mechanics and Industry, Vol. 14, Issue 4, 2013, p. 229238. [Publisher]
 Rahimi G. H., Gazor M. S., Hemmatnezhad M., et al. Free vibration analysis of fiber metal laminate annular plate by statespace based differential quadrature method. Advances in Materials Science and Engineering, Vol. 3, Issue 2014, 2014, p. 653659. [CrossRef]
 Mahi A., Bedia E. A. A., Tounsi A. A new hyperbolic shear deformation theory for bending and free vibration analysis of isotropic, functionally graded, sandwich and laminated composite plates. Applied Mathematical Modelling, Vol. 39, Issue 9, 2015, p. 24892508. [Publisher]
 Iriondo J., Aretxabaleta L., Aizpuru A. Characterisation of the elastic and damping properties of traditional FML and FML based on a selfreinforced polypropylene. Composite Structures, Vol. 131, 2015, p. 4754. [Publisher]
 Sayyad A. S., Ghugal Y. M. On the free vibration analysis of laminated composite and sandwich plates: A review of recent literature with some numerical results. Composite Structures, Vol. 129, 2015, p. 177201. [Publisher]
 Li R., Chen X. H., Liu P. Y. Free vibration test and numerical analysis of TiC fiber metal laminate. Acta Materiae Compositae Sinica, Vol. 33, Issue 5, 2016, p. 10641071. [CrossRef]
 Xing Y. F., Liu B. The Precise Solution of Free Vibration of Plates and Shell. Science Press, 2015. [CrossRef]
 Kou H., Yuan H. Rubinduced nonlinear vibrations of a rotating large deflection plate. International Journal of NonLinear Mechanics, Vol. 58, 2014, p. 283294. [Publisher]