Measurement of vibrations of composite wings using high-order finite element beam
Ismail Bennamia1 , A. Badereddine2 , T. Zebbiche3
1Institute of Aeronautics and Space Studies, University of Blida 1, B.P. 270, Blida 09000, Algeria
2, 3Aeronautical Sciences Laboratory, Institute of Aeronautics and Space Studies, University of Blida 1, B.P. 270, Blida 09000, Algeria
Journal of Measurements in Engineering, Vol. 6, Issue 3, 2018, p. 143-154.
Received 26 June 2018; received in revised form 9 September 2018; accepted 18 September 2018; published 30 September 2018
The purpose of this study is to show the influence of bending-torsion coupling on natural frequencies and mode shapes of aircraft wings by using two finite element beam formulations. The bending-torsion coupling parameters are the geometric parameter (distance between the mass axis and elastic axis of the cross-section of the beam) and the material coupling due to laminated composites. Cubic and high-order Hermite finite element interpolations are presented in this study, in order to show the influence of geometric and material coupling on natural frequencies and mode shapes. Starting by the governing partial differential equations of motion for the coupled bending-torsion beam with the bending and torsion equations, the Galerkin’s method is used with high-order finite element interpolation to obtain the high-order Hermitian shape functions. The mass and stiffness matrices are obtained using the kinetic energy and potential energy, respectively. The beam finite element has two nodes, the cubic element has three degrees of freedom at each end (transvers displacement, slope and torsion), where the high order element has five degrees of freedom at each end (transvers displacement, slope, curvature, gradient of curvature and torsion). The mass matrix contains geometric coupling terms and the stiffness matrix contains terms of material bending-torsion coupling. The obtained results using cubic and high-order finite element Euler-Bernoulli beam formulations are compared for a free vibration analysis of Goland metallic wing (geometric coupling) and validated with Dynamic Stiffness Method for composite wings.
Keywords: composite beam, Hermite interpolation, high-order finite element beam, natural frequencies and mode shapes.
Many analytical and numerical methods (Finite Element Method, Dynamic Stiffness method DSM, etc.) are adopted to determine approximate natural frequencies and mode shapes of uniform beams . A high-order finite element formulation is developed by Ganesan and Zabehollah  for vibration analysis of tapered composite beams. Vibrations measurements of a wing are investigated experimentally by using the method of time averaged projection moiré analyzed by Maskeliūnas et al. . Omarov et al.  made a dynamic model to determinate the reduced mass and stiffness of flexural vibrating cantilever beams. The measurement of plane vibrations of a two elastic structure are analyzed by Maskeliūnas et al. .
A dynamic stiffness element for free vibration analysis of composite beams and its application to aircraft wings is developed by Banerjee et al. , where a dynamic stiffness matrix of a composite beam is used to investigate it’s free vibration characteristics. A finite element parametric modeling technique of aircraft wing structures is given by Tang and Xi. . The composite beam models are used to study the dynamic response and aeroelasticity for high-aspect-ratio of composite wings [7-10].
A variable-order finite-element method with application to the free-vibration of rotating beams is described by Hodge et Rutkowski , the elasto-dynamic response of a surface stiffened transversely isotropic half-space subjected to an internal time-harmonic normal load is studied in axisymmetric time-harmonic response of a surface-stiffened transversely isotropic half-space by Eskandari et al. . An axisymmetric response of a bi-material full-space reinforced by an interfacial thin film presented by Ahmadi et al.  where Effects of thin film stiffness, material properties, loading depth, and surface/interface effect are studied.
A Kirchhoff thin plate perfectly bonded to the half-space is used to model the surface coating to study the Green’s functions of a surface-stiffened transversely isotropic half-space, Eskandari and Ahmadi .
A new novel collaborative optimization (MGACACO) algorithm based on the GA, ACO, the chaotic optimization method, multi-population strategy, adaptive control parameters and collaborative strategy is proposed to solve the complex optimization problems by Deng et al. [15, 16]. A PSO algorithm is used to optimize the parameters of least squares support vector machines (LS-SVM)  and solve multi-objective gate assignment . A novel two-stage hybrid swarm intelligence optimization algorithm . Other algorithms are used to solve optimization problems with application: a novel parallel hybrid intelligence optimization algorithm , A new feature extraction method based on EEMD and multi-scale fuzzy entropy for motor bearing , a novel collaborative optimization algorithm in solving complex optimization problems , and an improved self-adaptive differential evolution algorithm and its application . A novel vibration suppression method based on fractional order Proportional-Integral-Derivative (PID) controller is proposed by Zhao et al.  in research on vibration suppression method of alternating current motor based on fractional order control strategy.
This paper presents the effect of both geometric and material coupling in free vibration analysis of coupled bending-torsional beams by using a high-order finite element formulation. The bending-torsion coupling coefficient considered here is given without definition of the fiber orientation of laminate beams and the results are validated with those given by the Dynamic Stiffness Method DSM for Euler-Bernoulli beam model .
2. Finite element modeling
The cubic and high-order finite element interpolations are adopted to model the wing structure and deduce its global mass and global stiffness matrices. The kinetic and potential energies are used to obtain the mass matrix and the stiffness matrix of the beam element . The mass matrix contains geometric coupling terms due to (distance between the mass and elastic axis) and the stiffness matrix contains the terms of rigidity bending-torsion coupling (laminated beam or multilayered beams).
Fig. 1. Coupled Euler-Bernoulli beam
2.1. The bending-torsion coupled beam element
The governing partial differential equations of motion for the coupled bending-torsion free natural vibration of the composite beam are given by [1, 25, 26]:
where is the transverse displacement of the beam, is the bending stiffness, is the torsion stiffness and the rigidity of bending-torsion coupling, m, is the mass per unit in length and the moment of inertia about the axis elastic, is the geometric coupling.
2.2. The Euler-Bernoulli beam model
The Euler-Bernoulli equation of beam bending (without externally pressure loading) is :
where is mass density per volume.
2.2.1. Cubic finite element Hermitian beam
The cubic finite element Hermitian beam developed and described by Kown and Bang [27, 31], the cubic element has two nodes at the ends and at each end two degrees of freedom, the transverse displacement and the slope [28, 31]:
The evaluation of the transverse displacement and slope at both nodes gives the followed system: ; ; ; ; is the length of the beam element:
Solving the system Eq. (5) to obtain the coefficients and the transverse displacement assumption gives:
where the Hermitian shape functions of the beam (see Appendix).
2.2.2. High-order finite element beam
The weighted residual, Galerkin’s method is presented, to the beam Eq. (2) to develop the high-order finite element interpolation of the beam and the corresponding matrix equations.
The high-order finite element beam element presented in this study, has two nodes at the ends and at each end node four degrees of freedom , where is the transverse displacement, is the slope, is the curvature and is the gradient of curvature.
The assumption of deflection function is given using high order finite element polynomial function for, the degrees of this polynomial function is seven (eight nodal variables for the beam element):
From the Eq. (8), we can obtained the nodal variables for each node where:
– The slope is computed as:
– The curvature is:
– The gradient of curvature is:
The evaluation of these nodal variables on each node of the beam element gives:
where ; are the shape functions of the high order beam finite element (see Appendix).
2.2.3. The torsion in the beam
A linear shape functions are given for the torsion of the beam . The beam element has two nodes as described in previous sections each node has one degrees of freedom and . The linear function is:
The introducing of boundary conditions in the element , gives:
where are the linear shape functions of the beamin torsion (see Appendix).
2.3. Modal analysis
The free vibration analysis (modal analysis) plays an important role in aircraft design. The energy method is used to determinate the global mass and stiffness matrices of the beam. The natural frequencies and mode shapes are obtained by solving the eigenvalue matrix equation of the wing structure.
2.3.1. Kinetic energy
The kinetic energy of bending-torsion coupled beam is :
The kinetic energy can be written as [28, 29]:
where is the mass matrix of bending-torsion coupled beam [27, 31]:
2.3.2. Potential energy
The potential energy of bending-torsion coupled composite beam is :
The potential energy can be written as [28, 29]:
where is the stiffness matrix of bending-torsion coupled beam :
where the stiffness matrix of bending, is the coupled stiffness matrix and is the stiffness matrix of torsion.
2.3.3. Natural frequencies and mode shapes
The matrix system obtained by the discrete equations of motion for wing structure is [1, 3, 11]:
where and are the global mass and stiffness matrices, respectively.
3. Results and discussions
The free vibration analysis of composite wing is investigated using a cubic and high-order finite element Euler-Bernoulli beam formulations in order to show the effect of geometric coupling shown in Fig. 1 on vibration frequencies and mode shapes for metallic wings and generalized to composite wings where the anisotropy of materials has significant effects on natural frequencies and mode shapes.
3.1. Isotropic beam model (Goland wing)
The Goland’s wing properties are given in Table 1 [1, 7, 30], in order to validate the results obtained by the present work with the Dynamic Stiffness Method results, for a free vibration analysis of metallic wings.
The three first frequencies (rad/s) obtained by the cubic and high-order finite element beam model are compared with the DSM results Banerjee et al.  where the geometric coupling represented by rang of values of are given in Table 2.
Table 1. Goland wing properties
GJ K (Nm2)
Table 2. The three first frequencies of Goland wing
The mode shapes obtained for the uncoupled case and is shown in Fig. 2 (the geometric and material coupling are equal to zero, and this case is similar to Euler-Bernoulli classical beam theory). The obtained results in Fig. 2 show the effect of the geometric coupling on mode shapes and natural frequencies of Goland wing.
It appears clearly that the natural frequencies in Table 2 and the mode shapes in Fig. 2 for the cubic formulation and high-order formulation are very similar to those obtained by Banerjee et al.  for metallic wing using DSM approximation method varying the geometric coupling which affect the naturel frequencies.
Fig. 2. Mode shapes for various geometric coupling of the wing
3.2. Composite wing results
Various values of material coupling are considered to illustrate the effect of bending-torsion coupling due to material anisotropy in mode shapes and natural frequencies.
The material coupling rigidities given in this section are proposed by Banerjee et al. . These rigidities can be obtained by laminated composites.
The mode shapes for range of values of geometric coupling and various material coupling are shown in Figs. 3, 4 and 5. Observe that the obtained mode curves are different, in which the bending and torsion coupled modes is remarkable in each material coupling value.
Fig. 3. Mode shapes for rigidity coupling 1.5×106 and various geometric coupling of cantilever wing
The three first natural frequencies obtained from modal analysis of composite beam using cubic and high-order finite element beam are presented in Tables 3, 4 and 5 varying both geometric and material coupling for the beam, these obtained results are compared to those obtained by Banerjee et al.  using Dynamic Stiffness Method.
Table 3. First frequency results of composite wing
DSM results 
Hermite FEM beam model
High-order FEM Beam model
Table 4. Second frequency results of composite wing
DSM results 
Hermite FEM beam model
High-order FEM Beam model
Fig. 4. Mode shapes for rigidity coupling 2×106 and various geometric coupling of cantilever wing
It appears clearly from Tables 3, 4 and 5 that the three first naturel frequencies , and are almost similar for the Dynamic Stiffness Method  with the cubic and high-order finite element formulations varying both geometric coupling and material coupling.
Figs. 3, 4 and 5 show the bending and torsion mode shapes of the three first natural frequencies of cubic and high-order formulations for the coupled bending-torsional composite wing. The coupled case 0.2 m and 2.0×106 N.m2 curves are almost similar to those given by Banerjee et al. .
Table 5. Third frequency results of composite wing
DSM results 
Hermite FEM beam model
High-order FEM Beam model
Fig. 5. Mode shapes for rigidity coupling 2.5×106 and various geometric coupling of the wing
The free vibration analysis by measurement of natural frequencies and mode shapes of composite wings using cubic and high-order finite element coupled bending-torsion beam is presented and validated in this study. A high-order finite element Euler-Bernoulli beam model has been developed for isotropic materials, in first time, and generalized to laminated composite beams.
The results obtained by modal analysis (natural frequencies and mode shapes) are validated with dynamic stiffness approximation method (DSM). The geometric coupling and material coupling have a very important role in dynamic and aeroelasticity (fluid/structure interaction) problems. The varying of these parameters for Goland metallic wing can be generalized for all metallic and composite wings, in which the choice of aerodynamic or structural design of wing (or tailplane) depends on them (both geometric and material coupling). The results obtained from this study can be used for aeroelasticity studies, dynamic and frequency, response for random excitations obtained from gusting forces and it’s can be generalized to all Aircraft structure.
- Banerjee J. R., Su H., Jayatunga C. A dynamic stiffness element for free Vibration analysis of composite beams and its application to aircraft wings. Computers and Structures, Vol. 86, 2008, p. 573-579. [Publisher]
- Ganesan R., Zabehollah A. Vibration analysis of tapered composite beams using a high-order finite element. Part 1: Formulation, Vol. 77, 2007, p. 306-318. [Search CrossRef]
- Maskeliūnas V., Maskeliūnas R., Paškevičius P., Ragulskis L. Measurement of Vibrations of Wings. Journal of Measurements in Engineering, Vol. 4, Issue 3, 2016, p. 183-187. [Publisher]
- Tamerlan Omarov, Kuralay Tulegenova, Yerulan Bekenov, Gulnara Abdraimova, Algazy Zhauyt, Muslimzhan Ibadullayev Determination of reduced mass and stiffness of flexural vibrating cantilever beam. Journal of Measurements in Engineering, Vol. 6, Issue 1, 2018, p. 1-9. [Publisher]
- Maskeliūnas R., Ragulskis K., Paškevičius P., Patašienė L., Pauliukas A., Ragulskis L. Measurement of plane vibrations of a two dimensional elastic structure. Journal of Measurements in Engineering, Vol. 3, Issue 2, 2015, p. 42-47. [Search CrossRef]
- Tang J., Xi P. A finite element parametric modeling technique of aircraft wing structures. Chinese Journal of Aeronautics, Vol. 26, Issue 5, 2013, p. 1202-1210. [Publisher]
- Cesnik Carlos E. S., Hodgesy Dewey H., Mayuresh Patil J. Aeroelastic analysis of composite wings. Structural Dynamics and Materials Conference, Salt Lake City, Utah, 1996. [Search CrossRef]
- Mayuresh Patil J., Hodges Dewey H., Cesnik Carlos E. S. Nonlinear aeroelastic analysis of aircraft with high-aspect-ratio wings. 39th AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference and Exhibit, Structures, Structural Dynamics, and Materials and Co-located Conferences, 1998. [Search CrossRef]
- Zhao Y. H., Hu H. Y. Structural modeling and aeroelastic analysis of high-aspect-ratio composite wings. Chinese Journal of Aeronautics, Vol. 18, Issue 1, 2005, p. 25-30. [Publisher]
- Haddadpour H., Kouchakzadeh M. A., Shadmehri F. Aeroelastic instability of aircraft composite wings in incompressible flow. Composite Structures, Vol. 83, 2008, p. 93-99. [Publisher]
- Hodges D. Y., Rutkowski M. Y. Free-vibration analysis of rotating beams by a variable-order finite-element method. AIAA Journal, Vol. 19, Issue 11, 1981, p. 1459-1466. [Publisher]
- Eskandari M., Samea P., Ahmadi S. F. Axisymmetric time-harmonic response of a surface-stiffened transversely isotropic half-space. Meccanica, Vol. 52, Issues 1-2, 2017, p. 183-196. [Publisher]
- Ahmadi S. F., Samea P., Eskandari M. Axisymmetric response of a bi-material full-space reinforced by an interfacial thin film. International Journal of Solids and Structures, 90, p. 251-260. [Search CrossRef]
- Eskandari M., Ahmadi S. F. Green’s functions of a surface-stiffened transversely isotropic half-space. International Journal of Solids and Structures, Vol. 49, Issues 23-24, 2012, p. 3282-3290. [Publisher]
- Deng Wu, Zhang Shengjie, Zhao Huimin, Yang Xinhua A novel fault diagnosis method based on integrating empirical wavelet transform and fuzzy entropy for motor bearing. IEEE Access, Vol. 6, 2018, p. 35042-35056. [Publisher]
- Deng Wu, Zhao Huimin, Liu Jingjing, Yan Xiaolin, Li Yuanyuan, Yin Lifeng, Ding Chuanhua An improved CACO algorithm based on adaptive method and multi-variant strategies. Soft Computing, Vol. 19, Issue 3, 2015, p. 701-713. [Publisher]
- Deng Wu, Yao Rui, Zhao Huimin, Yang Xinhua, Li Guangyu A novel intelligent diagnosis method using optimal LS-SVM with improved PSO algorithm. Soft Computing, 2017, https://doi.org/10.1007/s00500-017-2940-9. [Publisher]
- Deng Wu, Zhao Huimin, Yang Xinhua, Xiong Juxia, Li Bo Study on an improved adaptive PSO algorithm for solving multi-objective gate assignment. Applied Soft Computing, Vol. 59, 2017, p. 288-302. [Publisher]
- Deng Wu, Chen Rong, He Bing, Liu Yaqing, Yin Lifeng, Guo Jinghuan A novel two-stage hybrid swarm intelligence optimization algorithm and application. Soft Computing, Vol. 16, Issue 10, 2012, p. 1707-1722. [Publisher]
- Deng Wu, Chen Rong, Gao Jian, Song Yingjie, Xu Junjie A novel parallel hybrid intelligence optimization algorithm for function approximation problem. Computers and Mathematics with Applications, Vol. 63, Issue 1, 2012, p. 325-336. [Publisher]
- Zhao Huimin, Sun Meng, Deng Wu, Yang Xinhua A new feature extraction method based on EEMD and multi-scale fuzzy entropy for motor bearing. Entropy, Vol. 19, Issue 1, 2017, p. 14. [Publisher]
- Deng Wu, Zhao Huimin, Zou Li, Li Guangyu, Yang Xinhua, Wu Daqing A novel collaborative optimization algorithm in solving complex optimization problems. Soft Computing, Vol. 21, Issue 15, 2017, p. 4387-4398. [Publisher]
- Deng Wu, Yang Xinhua, Zou Li, Wang Meng, Liu Yaqing, Li Yuanyuan An improved self-adaptive differential evolution algorithm and its application. Chemometrics and Intelligent Laboratory Systems, Vol. 128, 2013, p. 66-76. [Publisher]
- Zhao H. M., Li D. Y., Deng W., Yang X. H. Research on vibration suppression method of alternating current motor based on fractional order control strategy. Proceedings of the Institution of Mechanical Engineers Part E-Journal of Process Mechanical Engineering, Vol. 231, Issue 3, 2017, p. 786-799. [Publisher]
- Banerjee J. R. Explicite analytical expressions for frequency equation and mode shapes of composite beams. International Journal of Solids and Structures, Vol. 38, 2001, p. 2415-2426. [Publisher]
- Guo S. Aeroelastic optimisation of an aerobatic aircraft wing structure. Aerospace Science and Technology, Vol. 11, 2007, p. 396-404. [Publisher]
- Kwon W. Y., Bang H. The Finite Element Method using MATLAB. CRC Press, New York, 1997. [Search CrossRef]
- Tanka M., Bercin A. N. Finite element bending of a coupled bending and torsional free vibration of uniform beams with an arbitrary cross-section. Journal of Applied Mathematics and Modeling, Vol. 21, 1997, p. 339-344. [Publisher]
- Oz H. R. Calculation of the natural frequencies of a beam mass system using finite element method. Journal of Mathematical and Computational Applications, Vol. 2, 2000, p. 67-75. [Publisher]
- Kirsch Bertrand, Montagnier Olivier, Bénard Emmanuel, Faur Thierry M. Dynamic aeroelastic simulation of composite wing for HALE UAV application. 7th European Conference for Aeronautics and Aerospace Sciences (EUCASS), 2017. [Search CrossRef]
- Wright Jan R., Cooper E. Introduction to Aircraft Aeroelasticity and Loads. John Wiley and Sons, Chichester, 2007. [Publisher]
Multi-mode propagation and diffusion analysis using the three-dimensional second strain gradient elasticity
Mechanical Systems and Signal ProcessingBo Yang, Mohamed Ichchou, Abdelmalek Zine, Christophe Droz
AerospaceZawar Haider, Rafic M. Ajaj, Lakmal Seneviratne
Two-dimensional periodic structures modeling based on second strain gradient elasticity for a beam grid
International Journal of Mechanical SciencesBo Yang, Abdelmalek Zine, Christophe Droz, Mohamed Ichchou
IET Collaborative Intelligent ManufacturingXinhua Wang, Xiulin Zhang, Huajun Gong, Ju Jiang, Hari Mohan Rai