Next Article in Journal
Field Applications of Ultra High Frequency Techniques for Defect Detection in GIS
Previous Article in Journal
A Distributed and Context-Aware Task Assignment Mechanism for Collaborative Mobile Edge Computing
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Optimal Sensor Placement Based on Eigenvalues Analysis for Sensing Deformation of Wing Frame Using iFEM

1
Key Laboratory of Electronic Equipment Structure Design of Ministry of Education, Xidian University, Xi’an 710071, China
2
Xinjiang Observatory, National Astronomical Observatories, Chinese Academy of Sciences, Urumqi 830011, China
*
Author to whom correspondence should be addressed.
Sensors 2018, 18(8), 2424; https://doi.org/10.3390/s18082424
Submission received: 13 June 2018 / Revised: 20 July 2018 / Accepted: 21 July 2018 / Published: 25 July 2018
(This article belongs to the Section Physical Sensors)

Abstract

:
For real time monitoring of the wing state, in this paper, the inverse Finite Element Method (iFEM) is applied, which describes the displacement field of beam according to the Timoshenko theory, to sense the wing frame deformation. In order to maintain the accuracy and stability of frame deformation sensing with iFEM, an optimal placement model of strain sensors based on eigenvalue analysis is constructed. Through the model solution with the Particle Swarm Optimization (PSO) algorithm, two different optimal placement schemes of sensors are obtained. Finally, a simulation is performed on a simple cantilever beam and a static load experiment is conducted on an aluminum alloy wing frame. The results demonstrate that the iFEM is able to accurately sense the deformation of the wing frame, when the two optimal placement schemes of sensors are used.

1. Introduction

Deformation sensing in real-time is an essential technology for providing feedback to the actuation and control systems of smart structures, especially of next-generation aircrafts. Moreover, when the detailed state of structural deformation is known, other necessary response quantities, such as stress and failure criteria, can also be assessed [1,2,3]. A research of deformation sensing focuses on the utilization of in situ strain measurements, captured from a network of strain sensors, to estimate the structural deformation, which is commonly referred to as shape sensing. Specifically, fiber Bragg grating (FBG) sensors have been extensively studied for shape sensing due to the corresponding lightness, accuracy and ease of embedding [4,5,6,7].
Nevertheless, the structural deformation sensing from in situ strain data always represents an inverse problem. The inverse problem is frequently ill-posed, since it does not satisfy conditions of existence, uniqueness and stability [8,9]. Although many types of inverse problems and their applications have been proposed, few researchers have dealt with the deformation sensing of the wing shape. The current methods used to sense the wing shape can be divided into two categories [10]: one category focuses on the deformation sensing of the skin comprised of plate/shell structures to reflect the deformation of the wing shape and the other category emphasizes on the deformation sensing of a wing frame made of beam structures.
The effectiveness of the modal transformation scheme was analyzed and verified by Bogert et al. [11], in which the in-situ surface strain is used to sense the bending plate deformation. Despite the advantage of this method, accurate mode shapes and extensive eigenvalue analysis are required; for eigenvalue analysis, the elastic and inertial material properties need to be described in detail for application in high-fidelity finite element methods. On the basis of the least-squares and first-order plate theory, which includes the membrane and the bending and transverse shear deformations, a mathematical framework was developed by Tessler and Spangler to conduct the deformation sensing of plate and shell—called the inverse Finite Element Method (iFEM) [8,12]. The advantage of iFEM is that the structural deformations are sensed only from strain measurements, without the finite element model and a priori knowledge, such as damping, loading and material properties [13,14]. Through the iFEM framework and the Refined Zigzag Theory (RZT), Cerracchio and Gherlone et al. explored an effective way to sense the shape and stress of multilayered composites and sandwich structures [15,16,17]. Nevertheless, the kinematic variables in the above iFEM framework are derived from the three-node inverse-shell (i3-RZT) element or three-node inverse-plate (iMIN3) element, without including the hierarchical drilling rotation degrees-of-freedom (DOF). In order to extend the practical application of iFEM for the shape-sensing analysis of large-scale structures (such as the ship), a four-node quadrilateral inverse-shell element named iQS4 was formulated by Kefal et al. that utilized the kinematic assumptions of the first-order and transverse-shear deformation theory, to add the hierarchical drilling rotation DOF into the kinematic variables [18,19]. In a series of works [14,15,16,17,18,19,20,21,22,23], the capability of iFEM framework is verified through simulation or experimentation for 3D deformation sensing of plate/shell structures of aerospace vehicles and ships. However, for the aforementioned iFEM schemes, the strain sensors need to be attached on the top and bottom of the plate/shell structures. When the sensors are placed on the exterior surface of the structure, the protection mechanism are performed on the FBG or resistance strain sensors, which affects the aerodynamic performance of the wing. To avoid this problem, the strain sensors can be placed on the wing frame covered by the wing skin. If so, the wing frame deformation can be calculated to reflect the wing shape change.
In order to conduct the deformation sensing of beam structures, Ko et al. developed a load-independent scheme to approximate the beam curvature, named Ko’s Displacement Theory [24]. The proposed scheme employs Euler-Bernoulli beam theory to account for the beam deformation, in which the discrete surface strain measurements are integrated by using piece-wise continuous polynomials to sense the beam deformation. The experimental tests show that Ko’s Displacement Theory is sufficiently accurate for the deflection sensing of aircraft wings [25,26]. Regardless of its advantage, the scheme is only suitable for one-dimensional deformation of a beam structure, requiring a high number of strain sensors. According to Timoshenko beam theory, Gherlone et al. analyzed the displacement field of a constant cross-section beam, consequently constructing the relationship between the displacement and the surface strain data through iFEM [27,28]. Furthermore, Gherlone et al. conducted a comparative research on the aforementioned three shape sensing approaches: the iFEM, the modal transformation method and the Ko’s Displacement Theory [10]. It is found that the iFEM method is slightly more accurate and attractive than the other two methods for shape sensing. The experimental verifications demonstrate that the accuracy of iFEM for the beam deformation sensing is severely affected by the locations where the strain sensors are placed but the optimal placement criterion of the sensors is not mentioned [29]. Meanwhile, the stability and existence of iFEM are also severely affected by the locations where the strain sensors are placed. For instance, all sensors are placed parallel to the centroidal axes of the beam element.
The existing optimal placements of the strain sensors are based on the mode shape analysis of the structure, such as the modal kinetic energy (MKE) method, the effective independence (EI), the modal assurance criteria (MinMAC) and the drive point residue (DPR) [30,31,32,33]. Based on a mode approach, Geng et al. proposed an optimization scheme and improved genetic algorithms (GA) for optimal FBG sensor placement; consequently, the best reconstruction effects are obtained for the reconstruction of flexible plate structures [34]. Yang et al. developed a robust optimal sensor placement for uncertain structures [35]. In their schemes, the fitness function based on the Fisher information matrix (FIM) is derived and extended to the interval parameters. Genetic algorithm is employed to obtain the optimal solution. Being different to the traditional sensors, FBG sensors possess different detecting probe length at different angles. In view of the above difference, Yi et al. developed a detective model for FBG sensor which fit its feature based on probability model and the optimal result of the FBG sensor placement is solved by Particle Swarm Optimization (PSO) algorithm [36]. In Reference [37], Soman et al. proposed an optimal placement model based on Modal Identification (MI) and Accurate Mode Shape Expansion (AMSE) to place the strain sensors and accelerometers simultaneously. The results indicate that the use of a multi-type sensor system can improve the quality of Structure Health Monitoring (SHM). Given the fact that the iFEM reconstructs the structure deformation without the mode shape, all previous optimal placement schemes of the strain sensors are not suitable for iFEM.
When the strain sensors are placed on the ill-suited locations, the coefficient matrix of the equations may be ill-posed since the essence of deformation sensing with iFEM focuses on the solution of a linear system of equations. In [38,39,40], the optimal placement of sensors is regarded as the parameter identification of linear system of equations, using the condition number as the optimal object to construct the optimal model of sensors. When the condition number in a matrix is higher than 1016, the solution algorithm can return results with no accuracy at all; thus, such a matrix is numerically singular and linear systems with this matrix are not solved [39]. It is considered that a matrix is morbid when the condition number of the matrix is higher than 103. Nevertheless, in our research, it is found that the condition number of the coefficient matrix in iFEM tends to fall into [103, 1016] when the coefficient matrix is nonsingular; consequently, the bigger condition number in this range does not mean that the coefficient matrix is morbid. Therefore, the condition number is not suitable for the optimal sensor placement for iFEM.
This paper focuses on the construction of an optimal sensor placement model, based on eigenvalue analysis of the relationship between the sensor placement and the stability of sensing of the beam deformation with iFEM. This is conducted with an aim to maintain the sensing stability of iFEM. The contents of this paper are: firstly, the beam deformation sensing process through iFEM based on Timoshenko beam theory was reviewed. Following, the instability of sensing is discussed; consequently, the optimal placement model of strain sensors based on eigenvalue analysis is constructed. The two optimal results are obtained through the model solution with PSO algorithm. Moreover, a high-fidelity finite element model of the beam is constructed with ANSYS, which generate the discrete strain data to replace the experimental strain measurements captured from strain sensors and which produce the deformation data to replace the experimental deformation. This is executed to verify the robustness of iFEM under the optimal placement schemes. Finally, to examine the iFEM sensing capability for the two optimal placement schemes, deformation sensing test is conducted on an aluminum alloy wing frame model.

2. Inverse Finite Element Method Algorithm for Beam Deformation

In the iFEM algorithm, the displacement fields of an isotropic, straight and circular cross section beam element can be described from Timoshenko beam theory [27,28,29]:
{ u x ( x , y , z ) = u ( x ) + z θ y ( x ) y θ z ( x ) v y ( x , y , z ) = v ( x ) z θ x ( x ) u z ( x , y , z ) = w ( x ) + y θ x ( x ) ,
where, u x ,   v y , w z are the displacements along the x , y   a n d   z axes; u ( x ) , v ( x ) and w ( x ) denote the displacements along the center axis ( x [ 0 , l ] ,   y = z = 0 ); θ x ( x ) ,   θ y ( x )   a n d   θ z ( x ) are the rotations around the three coordinate axes (see Figure 1). The six kinematic variables are grouped in vector form as u { u ( x ) ,   v ( x ) ,   w ( x ) ,   θ x ,   θ y ,   θ z } T .
According to the small-strain hypothesis, the linear strains are given by Equation (1):
{ ε x ( x , y , z ) = e 1 ( x ) + z e 2 ( x ) + y e 3 ( x ) γ x z ( x , y ) = e 4 ( x ) + y e 6 ( x ) γ x y ( x , z ) = e 5 ( x ) z e 6 ( x ) ,  
where, the section strains e ( u ) = [ e 1 ,   e 2 , , e 6 ] T are theoretically defined as:
  e 1 ( x ) = u x , x ( x ) e 3 ( x ) = θ z , x ( x ) e 5 ( x ) = v y , x ( x ) θ z ( x ) e 2 ( x ) = θ y , x ( x ) e 4 ( x ) = w z , x ( x ) + θ y ( x ) e 6 ( x ) = θ x , x ( x ) .  
The iFEM senses the beam deformation shape by minimizing a weighted least-squares function φ , containing the section strains vector e ε computed from in situ surface strain data and the theoretical section strains vector e ( u ) defined through Equation (3):
φ ( u ) = e ( u ) e ε 2 .
In the finite element framework, the kinematic variables vector u(x) is obtained through interpolation from certain shape functions N ( x ) and the nodal degrees of freedom u e as following:
u ( x ) = N ( x ) u e .
With regard to different loading cases, the orders of shape functions N ( x ) differ. For the end-node loads, the order of N ( x ) is c 0 ; whereas for the uniformly distributed loads, the order is c 1 [28]. Thus, the total least-square function is a sum of N individual element contributions φ j ( u ( x ) , e ε ) :
  φ = j = 1 N φ j ( u ( x ) , e ε ) .  
Taking into account the effects of axial stretching, bending, twisting and transverse shearing, the element contributions φ j ( u ( x ) , e ε ) are given by the dot product as:
  φ j ( u ( x ) , e ε ) = k = 1 6 w k e φ k e ,  
with { w k e } = { w 1 0 , w 2 0 ( I y e A e ) , w 3 0 ( I z e A e ) , w 4 0 , w 5 0 , w 6 0 ( I p e A e ) }
  φ k e = l n i = 1 n [ e k ( x i ) e k ε ( x i ) ] 2 ( k = 1 , 2 , , 6 ) ,  
where, w k 0 ( k = 1 , 2 , , 6 ) are identically set as 1; A e , I y e , I z e and I p e are the cross-section area, the second moments of the area according to the y- and z-axis and the polar moment of the area of the beam element, respectively. l is the length of the beam element; x i ( 0 x i l ) and n are the axial coordinate of the locations where the section strains are evaluated and the number of locations, respectively. For the end-node loads, n m i n = 2 and for the uniformly distributed loads, n m i n = 3 [27,28,29,41].
The invoking of Equations (3) and (5) results in the theoretical section strains as following:
  e ( x ) = B ( x ) u e ,  
where, the matrix B ( x ) contains the derivatives of the shape functions N ( x ) .
Substituting Equation (9) into Equation (8), Equation (4) yields the following quadratic form:
  φ ( u ) = 1 2 ( u e ) T k e u e ( u e ) T f e + C e ,  
where, C e is a constant vector; k e and f e are defined as follows:
  k k e = L n i = 1 n [ B k T ( x i ) B k ( x i ) ]         f k e = L n i = 1 n [ B k T ( x i ) e k ε ( x i ) ] .
Finally, the minimization of function φ ( u ) in Equation (10) in terms of u e yields the sensing equation for the beam element deformation:
  k e u e = f e .
The two key steps in iFEM are: (1) the selection of suitable shape functions for the beam deformation sensing with iFEM; (2) the calculation of section strains from the measured surface strain data. For step (1), literatures [27,28,29,41,42] provide the detailed derivation process and results. For step (2), the section strains are computed through the following equation [27,28]:
  ε ( x i , θ i , β i ) = e 1 ε ( x i ) ( c β 2 v s β 2 ) + e 2 ε ( x i ) ( c β 2 v s β 2 ) s θ R + e 3 ε ( x i ) ( c β 2 v s β 2 ) c θ R  
+ e 4 ε ( x i ) c β s β c θ e 5 ε ( x i ) c β s β s θ + e 6 ε ( x i ) c β s β R              = [ c β 2 v s β 2 , ( c β 2 v s β 2 ) s θ R , ( c β 2 v s β 2 ) c θ R , c β s β c θ , c β s β s θ , c β s β R ] × e ε ( x i )
  = T ( x i , θ i , β i ) × e ε ( x i )  
with   c β   cos β i , s β s i n β i , c θ c o s θ i , s θ s i n θ i ,
where, e ε ( x i ) = { e 1 ε ( x i ) , e 2 ε ( x i ) , , e 6 ε ( x i ) | ( i = 1 , 2 , , m ) } is the in-situ section strains vector at location x i , along the x-axis. R is the external radius of the beam element. ε ( x i , θ i , β i ) denotes the measured surface strain at the location ( x i , θ i , β i ) , which is signed by the cylindrical coordinate system (see Figure 2). T ( x i , θ i , β i ) indicates the transformation relationship between the surface strain measurements and the section strains.

3. Construction of Optimal Placement of Sensors

When the section strains vector e ε is calculated from the surface strain measurements only with Equation (13), the strain measurements ε ( x i ,   θ i ,   β i ) (i = 1, 2…6) and the section strains vector e ε must be distributed in one cross-section. Consequently, 6 × n strain measurements are required to calculate the n section strain vectors. In the case of end-node loads, 12 strain measurements are required for the calculation of Section 2 strain vectors; in the case of uniformly distributed loads, 18 strain measurements are required for the calculation of Section 3 strain vectors. As aforementioned in the literatures [27,28,29,41], the number of strain measurements is reduced through the constitutive equations of Equation (14) (refer to Figure 3) and the equilibrium equations of Equation (15).
N = A x e 1 Q y = G y e 5 Q z = G z e 4 M x = J x e 6 M y = D y e 2 M z = D z e 3 ,
where, the section forces ( N ,   Q y   a n d   Q z ) and moments ( M x ,   M y   a n d   M z ) are related to the section strains e i ( x ) . A x = E A is the axial rigidity, where A is the area of the cross-section of the beam element. G y = k y 2 G A , G z = k z 2 G A are the shear rigidities with k y 2 and k z 2 denoting the shear correction factors, where G is the shear modulus. For the thin-walled section, k y 2 = k z 2 = 0.531; for the thick-walled section, k y 2 = k z 2 = 0.62; and for the solid section, k y 2 = k z 2 = 0.887 [28]. J x = G I p is the torsional rigidity. D y = E I y and D z = E I z denote the bending rigidities.
N x + q x = 0 Q y x + q y = 0 Q z x + q z = 0 M x x = 0 M y x Q z = 0 M z x Q y = 0 .  
When the forms of distribution loads q x , q y and q z are known, the forms of the section strains e ε can be estimated. For the end-node loads, the section stains e 1 ε ( x i ) , e 4 ε ( x i ) , e 5 ε ( x i ) and e 6 ε ( x i ) are constant and   e 2 ε ( x i )   and   e 3 ε ( x i ) are linear. For the uniformly distributed transverse loads, the section strains e 1 ε ( x i )   and   e 6 ε ( x i ) are constant, e 4 ε ( x i )   and   e 5 ε ( x i ) are linear and e 2 ε ( x i )   and   e 3 ε ( x i ) are parabolic. e 4 ε ( x i ) = E I z G k y 2 A e 2 ,   x ε ( x i ) = m 1 e 2 , x ε ( x i ) , e 5 ε ( x i ) = E I y G k z 2 A e 3 , x ε ( x i ) = m 2 e 3 , x ε ( x i )   [28,41]. Through the aforementioned results, the section strains are calculated as:
For the end-node loads, the section strains are expressed as:
e 1 ε ( x i ) = a 1 e 2 ε ( x i ) = a 2 x i + a 4 e 4 ε ( x i ) = m 1 a 2 e 6 ε ( x i ) = a 6 e 3 ε ( x i ) = a 3 x i + a 5 e 5 ε ( x i ) = m 2 a 3 ,
or in a matrix equation:
e ε ( x i ) = { e 1 ε ( x i ) , e 2 ε ( x i ) , , e 6 ε ( x i ) } T                   = [ 1 0 0   0 x i 0   0 0 x i   0 0 0   m 1 0 0   0 m 2 0 0 1 0   0 0 1   0 0 0   0 0 0   0 0 0   0 0 1 ] × [ a 1 , a 2 , a 3 , a 4 , a 5 , a 6 ] T = T 1 ( x i ) × p 1 ,
where, p 1 = [ a 1 , a 2 , a 3 , a 4 , a 5 , a 6 ] T is a constant parameters vector; T 1 x i indicates the transfer matrix between p 1 and the section strains vector e ε ( x i ) .
For the uniformly distributed transverse loads, the section strains are expressed as:
e 1 ε ( x i ) = b 1 e 2 ε ( x i ) = b 2 x i 2 + b 3 x i + b 4 e 4 ε ( x i ) = 2 m 1 b 2 x i + m 1 b 3 e 6 ε ( x i ) = b 8 e 3 ε ( x i ) = b 5 x i 2 + b 6 x i + b 7 e 5 ε ( x i ) = 2 m 2 b 5 x i + m 2 b 6 ,
or in a matrix equation:
e ε ( x i ) = { e 1 ε ( x i ) , e 2 ε ( x i ) , , e 6 ε ( x i ) } T  
= [ 1 0 0 0 0 0 0 0 0 x i 2 x i 1 0 0 0 0 0 0 0 0 x i 2 x i 1 0 0 2 m 1 x i m 1 0 0 0 0 0 0 0 0 0 2 m 2 x i m 2 0 0 0 0 0 0 0 0 0 1 ] × [ b 1 , b 2 , b 3 , b 4 , b 5 , b 6 , b 7 , b 8 ] T
= T 2 ( x i ) × p 2 ,
where, p 2 = [ b 1 , b 2 , b 3 , b 4 , b 5 , b 6 , b 7 , b 8 ] T is a constant parameters vector; T 2 x i indicates the transfer matrix between p 2 and the section strains vector e ε ( x i ) . Following, the relationship between the undetermined parameters vector p 1   o r   p 2 and any measured strain data ε ( x i ,   θ i , β i ) is expressed as:
  ε ( x i , θ i , β i ) = T ( x i , θ i , β i ) T 1 ( x i ) × ε ( x i , θ i , β i ) = T ( x i , θ i , β i ) T 2 ( x i ) × p 2 .
When the parameters vector p 1   o r   p 2 is solved, the arbitrary section strains vector is calculated through Equations (16) and (18).
  e ε ( x j ) = T 1 ( x j ) × ( T ( x i , θ i , β i ) T 1 ( x i ) ) 1 × ε ( x i , θ i , β i ) or   e ε ( x j ) = T 2 ( x j ) × ( T ( x i , θ i , β i ) T 2 ( x i ) ) 1 × ε ( x i , θ i , β i ) ,
where, i = 1 ,   2 , n   is the location where the section strains vector is calculated; n is the number of sections stated in Equation (8); j = 1 ,   2 , , m and m is the minimum number of the strain sensors used to capture the surface strains. The value of m is different under different loading cases analyzed through Equations (16) and (18). For the end-node loads, m = 6 , whereas for the uniformly distributed loads, m = 8 [28]. When the section strains vector e ε is obtained, the kinematic variables u e can be directly calculated from Equation (12).
However, through in-depth research, it is found that the transformation matrices T ( x i ,   θ i ,   β i ) T 1 ( x i ) and T ( x i ,   θ i , β i ) T 2 ( x i ) are ill-conditioned or even singular, when the strain sensors are placed at the inappropriate locations along the beam surface, such as β i in Equation (13) set to the same value 0 o . Thus, the deformation reconstructed from Equation (12) does not necessarily satisfy the conditions of existence, uniqueness and stability [8].
The constant parameter vectors p 1 and p 2 determination is regarded as the solution of a linear system of equations. The stability of Equation (20) solution depends on whether the product matrices T ( x i , θ i ,   β i ) T 1 ( x i ) and T ( x i ,   θ i ,   β i ) T 2 ( x i ) are well-conditioned or ill-conditioned. In actual measurements, the sensor placement is different from the pre-set placement. This means that one sensor is set at the node ( x i ,   θ i ,   β i ) ideally but the actual sensor location might be at the node ( x i ,   θ i ,   β i ) and the strain measurement contains the noise. Therefore, the strain input errors are added to Equation (20):
  ε ( x i ,   θ i , β i ) + Δ ε ( x i , θ i , β i ) = T ( x i , θ i , β i ) × H ( x i ) × ( p + Δ p ) Δ ε ( x i , θ i , β i ) = ε ( x i , θ i , β i ) ε ( x i , θ i , β i )   ( i = 1 ,   2 , , m ) ,
where, Δ ε ( x i , θ i , β i ) indicate the strain input errors, which couple the errors of sensor placement with the noise of the strain measurement system. ε ( x i , θ i , β i ) and ε ( x i , θ i , β i ) are the ideal strain input at the node ( x i , θ i , β i ) and the actual strain input, respectively. H = T 1   o r   T 2 , Δ p is the errors vector, resulting from the above errors. Following, the relative error between Δ p and p is estimated as:
Δ p p A × A 1 Δ ε ε   = c o n d ( A ) Δ ε ε w i t h   A = T ( x i , θ i , β i ) H ( x i ) ,   c o n d ( A ) = A × A 1 ,
where, · indicates the matrix norm and c o n d ( A ) signifies the condition number of matrix A . To a certain extent, the condition number can reflect the morbidity degree of a matrix. From Equation (23), it is observed that the relative error is high and the matrix is ill-posed, when the condition number is high. Kunsoo Huh and J.L. Stein proposed that the reason for the morbidity of a matrix is that a good distribution of eigenvalues for the matrix does not exist. For well distributed eigenvalues, large differences among the matrix eigenvalues exist [43], such as the high value in m i n | λ j λ k | ( j k , λ j and λ k are the eigenvalues of the matrix). Therefore, the optimal placement model of sensors is constructed as:
f i n d   f ( A ( x , θ , β ) ) = m a x ( min | λ j λ k | ) ( j k ,     j ,     k = 1 ,   2 , , m ) ,   ( ( x , θ , β ) = [ x i , θ i , β i , , x m , θ m , β m ] ) s . t .   x i [ L / 5 ,     4 L / 5 ] ,   θ i ϵ [ 180 ° , 180 ° ]   β i = 0 °   o r   45 ° ,    i = 1 ,   2 , m
where, λ i (i = 1, 2,…,m) are the eigenvalues of the matrix A . ( x ,   θ ,   β ) represents the locations where the sensors are placed and m is the number of the sensors used to capture the strain data of the structure. For the end-node loads, m = 6 , whereas for the uniformly distributed loads, m = 8 [28]. In view of application environment in engineering, it is difficult to set FBG strain sensors on the clamped node and on the free end node of the beam element. Moreover, the curve of the beam surface affects the strain measurement accuracy when a non-zero angle exists between the sensor and the generatrix of the beam surface [29]. Thus, it is difficult to set two or more FBG sensors at one node along the surface of the beam element, as the fiber sensors are not stacked (see Figure 4). Consequently, x i [ L / 5 ,   4 L / 5 ] is set, where only one sensor is placed at β i = 45 ° and the other sensors are placed at β i = 0 ° .
The optimal placement model of Equation (24) is a multi-parameter optimization problem, which is solved by using the PSO algorithm with good convergence speed. The optimal result is efficiently obtained in hyperspace.
In a population that includes N particles, the status of the i-th particle is described by the current location and the current velocity. a i = ( x i 1 , θ i 1 , β i 1 , , x i m , θ i m , β i m ) 1 3 m is the current location of the i-th particle and v i = ( v x 1 , v θ 1 , v β 1 , v x m , v θ m , v β m ) 1 3 m is the current velocity of the i-th particle ( i = 1 ,   2 , ,   N ; m   is   the   number   of   strain   sensors ) . The fitness function of PSO algorithm is defined as:
f ( a i ) = f ( A ( x , θ , β ) ) = m a x ( min | λ j λ k | ) ( j k , j , k = 1 , 2 , , m ) ,
( ( x , θ , β ) = ( x i 1 , θ i 1 , β i 1 , , x i m , θ i m , β i m ) 1 3 m ) ,
and the updated velocity and location are expressed as:
  v i k + 1 = α v i k + c 1 r 1 × ( p i k a i k ) + c 2 r 2 × ( p g k a i k ) a i k + 1 = a i k + v i k + 1 ,
where, α ,   c 1   a n d   c 2 are constants; r 1   a n d   r 2 are uniformly distributed random numbers between 0 and 1; a i k and v i k arethe current location and the current velocity of the i-th particle in the k-th iteration, respectively; p i k is the individual best location of the i-th particle at step k; p g k is the local best location of the population at step k. p i k and p g k are determined as:
  p i k = { p i k 1 ,   f ( a i k ) f ( p i k ) a i k ,   f ( a i k ) > f ( p i k )   p g k { p 0 k , p 1 k , , p N k } | g l o c k = f ( p g k ) = max { f ( p 0 k ) , f ( p 1 k ) , , f ( p N k ) } g g l o k = m a x { g l o c 0 , g l o c 1 , , g l o c k } ,
where, g l o c k and g g l o k are defined as the local extremum and the global extremum of the population, respectively. The key steps of the optimal model Equation (24) solution with PSO algorithm are described as:
Step 1.
PSO algorithm initializing. The size of the particle swarm is set to N = 50 and the maximum iteration is set to k m a x = 1000 ; the initial location a 0 and the initial velocity v 0 are set to the random values included in the constraint of Equation (24). The corresponding fitness value of each particle, which is calculated from the fitness function Equation (25), is set as the initial individual extremum f ( p i 0 ) , (i = 1,2,…,N). Following, the initial local extremum g l o c 0 and the global extremum g b e s t 0 of the population at the initial phase are selected for Equation (27). The threshold of the algorithm breaking is set to M = 1000 . When the global extremum g b e s t 0 is equal to M , the algorithm breaks and the current particle location is the optimal placement of the sensors. Else Step 2 will follow.
Step 2.
Velocity and location update of each particle with Equation (26). The fitness value of each particle f ( a i 1 ) , which is calculated from Equation (25), is compared to the initial individual extremum f ( p i 0 ) . If f ( a i 1 ) > f ( p i 0 ) , the individual extremum is replaced with f ( a i 1 ), or else, the individual extremum is invariant. Subsequently, the local extremum ( g l o c 1 ) and the global extremum g b e s t 1 of the population at step 1 are selected with Equation (27). When the global extremum g b e s t 1 is equal to M , the algorithm breaks; or else, the iteration continues and k = k + 1.
Step 3.
PSO algorithm termination. The iteration will not break until the iteration number k = k m a x or the current global extremum g b e s t k is equal to M .
The flow chat of above optimization procedure is shown in Figure 5.

4. Verifications through Simulations and Experimentation

In order to assess the optimal placement effectiveness of the sensors proposed in this paper, simulations and model experimentation are performed. In Section 4.1, a finite element model of a thin-walled beam with a circular cross-section is constructed through high-fidelity direct Finite Element analyses software (ANSYS 14.5, ANSYS, Southpointe, PA, USA), used to assess the robustness of the optimal placement schemes of the strain sensors. In Section 4.2, a wing frame model is tested under different static loads and presented, in order to assess the feasibility of iFEM for sensing wing frame deformation under the condition that the sensors are optimally placed.

4.1. Simulation Verification of Beam

To model the thin-walled beam structure, we made use of the Beam188 elements module with two nodes, based on the Timoshenko beam theory. The Young’s modulus is E = 75 , 000   MPa , the Poisson ratio is v = 0.3 and the density is ρ = 2557   kg / m 3 . The beam length is 660 mm, the external radius is R e x t = 13   mm and the thickness is s = 1.5   mm . The finite element model of the beam is divided into 200 elements and the cross-section of the beam is divided into 120 sectors (see Figure 6).
For comparison, two normal sensors placements (C1 and C2) are quoted from the literature [28], while two optimal results (C3 and C4) are obtained from the Equation (24) model with PSO. The four placements of the strain sensors are presented in Table 1, where C1 and C3 are suited to the end-node loading, C2 and C4 are suited to the uniformly distributed loading and f ( A ( x , θ , β ) ) indicates the function value in the optimal model of Equation (24).
Three different static loading tests are performed on the beam model. For the first test, the two end-node forces are loaded on the model, −100 N in direction Y and 80 N in direction Z (Loading A, Figure 7a); for the second test, two uniform distribution forces are loaded along the beam: −1 N in direction Y and 1.5 N in direction Z (Loading B, Figure 7b); for the third test, complex loads (Loading C) producing uniform distribution of loads (1 N in direction Y) with four equal node forces are loaded on the beam (150 N in direction Z, Figure 7c).
In lieu of the experimentally measured surface strains and the deformation results, the corresponding data are extracted from ANSYS. The deformations calculated from iFEM, the deformations extracted from ANSYS, the absolute error (AE) and the percentage errors (PE) for the free-end displacements along with the rotations of the beam model are presented in Table 2, Table 3 and Table 4. The displacements and rotations are expressed in mm and rad, respectively.
  P E = | δ i F E M δ A N S Y S | | δ A N S Y S | % A E = | δ i F E M δ A N S Y S |  
where, δ = ( u ( x ) ,   v ( x ) ,   w ( x ) ,   θ ( x ) ,   θ ( y ) ,   θ ( z ) ) ; the superscript ‘iFEM’ refers to the calculated value; and ‘ANSYS’ refers to the extracted value from ANSYS.
The comparisons of Table 2, Table 3 and Table 4 demonstrate that the two optimal placements (C3 and C4) present higher accuracy (maximum percent error is below 1%) for the iFEM beam deformation calculation compared to the schemes (C1 and C2) in the literature [28]. As discussed in section 3, strain input errors Δ ε ( x i , θ i , β i ) exist in Equation (22). It is assumed that the theoretical locations, where the strain sensors are placed, are unaltered. This means that the transformation matrix A = T ( x i , θ i , β i ) H ( x i ) is invariable but the actual locations, where the strain sensors are placed, will change. Thus, the strain inputs in Equation (22) do not correspond to the theoretical locations. The difference between the theoretical locations and the actual locations of strain sensors is set:
Δ x i [ 0.05 L , 0.05 L ] , Δ θ i [ 9 ° , 9 ° ]   and   Δ β i [ 9 ° , 9 ° ]   ( i = 1 ,   2 , , m ) .
The strain measurement system errors are assumed to obey the Gaussian error distributions, which have zero mean value and three-standard deviations equal to 5% of ANSYS simulation strain value, such as Δ ε i [ 5 % ε i ,   5 % ε i ] . These disturbances are added 1000 times and the highest errors are selected as the worst sensing results (presented in Table 5, Table 6 and Table 7, the displacements and the rotations are expressed in mm and rad, respectively).
The comparisons from Table 5, Table 6 and Table 7 show that the two optimal sensor placement schemes of C3 and C4 are more robust than the two quoted schemes of C1 and C2. For the schemes C1 and C2, the absolute and relative errors drastically increase; the maximum relative error is 1225.16% in direction Y and the maximum absolute error is 298.55 mm in direction Z in the case of complex loads. In contrast, for the two optimal sensor placement schemes, the maximum relative error is 103.95% in direction Y in the case of complex loads, whereas the maximum relative error is below 50% in the other two loading cases. The reason why difference increases sharply between the computed results and the extracted data is that the small changes in strain inputs cause large changes in the p 1 and p 2 parameter calculation through Equation (22). From Table 1, it is observed that the large differences among the matrix eigenvalues for the schemes of C3 and C4 (both exceed 10) are higher compared to the schemes of C1 and C2 (both equal to approximately 1). Consequently, the distributions of eigenvalues for C3 and C4 are better compared to the C1 and C2; and the changes of the calculated results for the C3 and C4 schemes are more stable compared to C1 and C2.

4.2. Experimental Verification of Wing Frame

The experimental tests under different static loadings are performed on a wing frame model, in order to verify the feasibility and effectiveness of iFEM structural deformation sensing for two optimal placements of sensors. The frame model is made of an aluminum alloy, combining two thin-walled beams with several thin-walled plates. The length and thickness of each beam are 2 m and 1.5 mm, respectively but the external radiuses differ for the two beams: one radius is 13 mm as the main spar, whereas the other radius is 11.5 mm. The entire frame is divided into three sections and the length of each section is L = 666 mm (see Figure 8A). Displacement measurements are conducted at different locations along the main spar with position sensors (see Figure 8A,C), which sends the infrared lights to the CCD cameras of the 3D optical measurement instruction (see Figure 8B NDI Optrotrak Certus, NDI, Canada) to reflect the structure deformation. The accuracy of NDI Optrotrak Certus is 0.1 mm in its measurement range. The displacements captured from NDI are used to assess the accuracy of deformation sensing from the strain data with iFEM.
The FBG strain sensor is a novel strain measurement device, based on the light wavelength shift which is caused by the FBG grating deformation generated by the tension/compressive force or change of temperature. The deformation for the unit length of the grating is labeled as the strain. Without considering the effect of the temperature change, the strain is calculated through Equation (30) [44,45].
ε i = K × ( λ e n d ( i ) λ i n i ( i ) λ i n i )   w i t h   K = 1 P e  
where, λ e n d ( i ) and λ i n i ( i ) are the wavelength shift and initial wavelength of i-th FBG sensor; P e is the photo-optical coefficient of the fiber; the strain measurement ε i is expressed as micro strain.
In the model test, the experimental strain data are obtained from the strain measurement system composed of FBG strain sensors (Fiber Bragg Grating|os1100, Micron Optics, Shelburne, VT, USA) and the FBG interrogator (Optical Sensing Instrument|Si 155, Micron Optics, Atlanta, GA, USA). Twenty-four FBG strain sensors (the range of initial wavelength is [1527 nm, 1564 nm]) are placed at different locations along the main spar and used to capture the surface strain. Since the locations in placement scheme C3 are contained in scheme C4 (refer to Table 1), eight FBG sensors are placed on the surface of every section according to scheme C4. The entire experiment system is presented in Figure 8D.
For the static tests, two different loading cases are considered (see Figure 8E,F):
I. the end-node of the main spar is loaded 6 times (for the total weight of F = 5.57 KG);
II. the main spar is loaded 5 times with uniform-distributed loads (for the total weight of F = 6.02 KG).
The details of the above loads, the maximum deformations calculated with iFEM and the maximum deformations captured by NDI are presented in Table 8. The optical measurement instruction only captures the displacement of position sensors stuck to the structure surface. Therefore, the comparisons among iFEM calculations and NDI measurements focuses on the displacements only. To visually describe the structural deformation, the comparisons of the discrete point displacements among iFEM calculations and NDI measurements are shown in Figure 9 (for the end-node load v and the uniformly distributed load iv).
The root-mean-square difference (RMSD) is used to assess the sensing accuracy of the entire structure through iFEM [46].
  R M S D = i = 1 j ( d i s p N D I ( x i ) d i s p i F E M ( x i ) ) 2 / j 2  
where, d i s p ( x i ) is the displacement of one node along the beam centroidal axis in one direction; the superscript ‘NDI’ refers to the deformation values captured from NDI; ‘iFEM’ refers to the predicted values computed from in situ strain measurements with iFEM; and j is the number of the nodes used to describe the beam deformation. In the test, the number of nodes used to describe the beam deformation is 16. The comparisons are shown in Table 9.
The frame is divided into three sections and every section is regarded as an iFEM element, so the entire frame is calculated three times through iFEM. The calculations in the second and third sections of the frame are affected by the computation in the first section. Consequently, accumulative errors occur during the deformation sensing of the entire frame structure. In Table 8, it is observed that the maximum differences of displacements among the iFEM calculations and the NDI captures are below 12%; for the main deformation v(x) under all loading cases, the errors remain at approximately 5%, when the deformations increase. Although the percentage errors are relatively high, the RMSD is acceptable for the C3 and C4 schemes: (1) under the end-node loads, the maximum RMSD is 0.644 mm for the main deformation v ( x ) and 0.234 mm for the minor deformation w ( x ) ; (2) under the uniformly distributed loads, the RMSD is 0.32 mm for the main deformation v ( x ) and 0.219 mm for the minor deformation w ( x ) (refer to Table 9). The corresponding figures also demonstrate that the deformations calculated with iFEM for the sensor placements of C3 and C4 are quite similar to the deformations measured with the optical measurement instruction (NDI). Nevertheless, the number of measurement strains for C4 is higher compared to the C3 number. Therefore, more disturbances exist for the C4 than for the C3, which results in the sensing accuracy of C4 being slightly worse than the C3. Finally, it is discovered that the accuracies in the model experiment are worse compared to the simulation, which also demonstrate that the practical location errors of the strain sensors affect the deformation sensing of the entire structure.

5. Conclusions

The study successfully validates the feasibility of iFEM used for deformation sensing of a wing frame. As aforementioned in the introduction, the locations where the strain sensors are placed affect the accuracy of iFEM for sensing. Therefore, in the paper, an optimal placement model of sensors is proposed to confirm the suitable placement of sensors for the accuracy and stability of iFEM of sensing frame deformation to be maintained. Following experimental test, the results demonstrate that the iFEM is able to precisely sense the deformation of the wing frame for the two optimal placement schemes of C3 and C4. Nevertheless, it is found that certain limits exist for the iFEM application: (1) the cross-section of the beam must be constant along the central axis of the beam and the section must be un-deformed when the loads are applied on the beam; (2) the boundary condition of the sensing equation must be known. Moreover, the test results also demonstrate that the strain measurement system errors (combining location errors of sensors with measurement errors of sensors) affect the accuracy of sensing results. Further works will be focused on the location calibrations of the strain sensors, as well as the relationship between the structural strain and the wavelength shift of the FBG sensor.

Author Contributions

Y.Z. carried out the majority of work on the paper, analyzing the singularity of the sensing model, constructing the optimal placement model of sensors is constructed, verifying the sensor optimal placement and writing the article; J.D. constructed the finite element model of the wing-like framework; H.B. and Q.X. gave advice regarding the writing of the article.

Funding

This work is financially supported by the National Natural Science Foundation of China through Grants 51775401 and 51490660, the National Key Basic Research Program of China through Grant 2015CB857100, the CAS “Light of West China” Program. No. 2016-QNXZ-A-7 and the joint Research Fund in Astronomy (No. U1731135) under cooperative between the National Science Foundation of China (NSFC) and Chinese Academy of Sciences (CAS).

Acknowledgments

Special appreciations are given to the editor and the anonymous reviewers for their careful reviews and valuable comments.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Esper, I.D.M.; Rosa, P.F.F. Heading controller for a fixed wing UAV with reduced control surfaces based on ANFIS. In Proceedings of the 45th Annual IEE/IFIP International Conference on Dependable Systems and Networks Workshops (DSN-W’15), Rio de Janeiro, Brazil, 22–25 June 2015; pp. 118–123. [Google Scholar]
  2. Kayacan, E.; Khanesar, M.A.; Rubio-Hervas, J.; Reyhanoglu, M. Learning Control of Fixed-Wing Unmanned Aerial Vehicles Using Fuzzy Neural Networks. Int. J. Aerosp. Mech., 2017, 2017, 1–12. [Google Scholar] [CrossRef]
  3. Pak, C. Wing shape sensing from measured strain. AIAA J. 2016, 54, 1064–1073. [Google Scholar] [CrossRef]
  4. Hong, C.Y.; Zhang, Y.F.; Zhang, M.X.; Leung, L.M.; Liu, L.Q. Application of FBG sensors for geotechnical monitoring, a review of sensor design, implementation methods and packaging techniques. Sens. Actuators A Phys. 2016, 244, 184–197. [Google Scholar] [CrossRef]
  5. Di Sante, R. Fiber optic sensors for structural health monitoring of aircraft composite structures: Recent advances and applications. Sensors 2015, 15, 18666–18713. [Google Scholar] [CrossRef] [PubMed]
  6. Vorathin, E.; Hafizi, Z.M.; Ghani, S.C.; Lim, K.S. Real-time monitoring system of composite aircraft wings utilizing Fibre Bragg Grating Sensor. Innov. Aerosp. Eng. Technol. 2016, 152, 012024. [Google Scholar] [CrossRef]
  7. Amanzadeh, M.; Aminossadati, S.M.; Kizil, M.S.; Rakić, A.D. Recent Developments in Fiber Optic Shape Sensing. Measurement 2018, 128, 119–137. [Google Scholar] [CrossRef]
  8. Tessler, A.; Spangler, J.L. A variational principle for reconstruction of elastic deformations in shear deformable plates and shells. In NASA Langley Research Center TM-212445; NASA Research Center: Hampton, VA, USA, 2003. [Google Scholar]
  9. Tikhonov, A.N.; Arsenin, V.Y. Solutions of Ill-Posed Problems; Winston and Sons: Washington, DC, USA, 1977. [Google Scholar]
  10. Gherlone, M.; Cerracchio, P.; Mattone, M. Shape sensing methods: Review and experimental comparison on a wing-shaped plate. Prog. Aerosp. Sci. 2018, 99, 14–26. [Google Scholar] [CrossRef]
  11. Bogert, P.; Haugse, E.; Gehrki, R. Structural shape identification from experimental strains using a modal transformation technique. In Proceedings of the 44th AIAA/ASME/ASCE/AHS Structures, Structural Dynamics, and Materials Confere Norfolk, Norfolk, VA, USA, 7–10 April 2003. [Google Scholar]
  12. Tessler, A.; Spangler, J.L. A least-squares variational method for full-field reconstruction of elastic deformations in shear-deformable plates and shells. Comput. Methods Appl. Mech. Eng. 2005, 194, 327–339. [Google Scholar] [CrossRef]
  13. Quach, C.; Vazquez, S.; Tessler, A.; Moore, J.; Cooper, E.; Spangler, J. Structural anomaly detection using fiber optic sensors and inverse finite element method. In Proceedings of the AIAA Guidance, Navigation, and Control Conference and Exhibit, San Francisco, CA, USA, 15–18 August 2005. [Google Scholar]
  14. Song, X.; Liang, D. Dynamic displacement prediction of beam structures using fiber bragg grating sensors. Optik 2018, 158, 1410–1416. [Google Scholar] [CrossRef]
  15. Cerracchio, P.; Gherlone, M.; Tessler, A. Real-time displacement monitoring of a composite stiffened panel subjected to mechanical and thermal loads. Meccanica 2015, 50, 2487–2496. [Google Scholar] [CrossRef]
  16. Cerracchio, P.; Gherlone, M.; Di Sciuva, M.; Tessler, A. A novel approach for displacement and stress monitoring of sandwich structures based on the inverse finite element method. Compos. Struct. 2015, 127, 69–76. [Google Scholar] [CrossRef]
  17. Cerracchio, P.; Gherlone, M.; Di Sciuva, M.; Tessler, A. Shape and stress sensing of multilayered composite and sandwich structures using an inverse finite element method. In Proceedings of the International Conference on Computational Methods for Coupled Problems in Science and Engineering-COUPLED PROBLEMS, Lbiza, Spain, 17–19 June 2013. [Google Scholar]
  18. Kefal, A.; Oterkus, E. Displacement and stress monitoring of a Panamax containership using inverse finite element method. Ocean Eng. 2016, 119, 16–29. [Google Scholar] [CrossRef] [Green Version]
  19. Kefal, A.; Tessler, A.; Oterkus, E. An enhanced inverse finite element method for displacement and stress monitoring of multilayered composite and sandwich structures. Compos. Struct. 2017, 179, 514–540. [Google Scholar] [CrossRef]
  20. Papa, U.; Russo, S.; Lamboglia, A.; Del Core, G.; Iannuzzo, G. Health structure monitoring for the design of an innovative UAS fixed wing through inverse finite element method (iFEM). Aerosp. Sci. Technol. 2017, 69, 439–448. [Google Scholar] [CrossRef]
  21. Kefal, A.; Oterkus, E. Displacement and stress monitoring of a chemical tanker based on inverse finite element method. Ocean Eng. 2016, 112, 33–46. [Google Scholar] [CrossRef] [Green Version]
  22. Kefal, A.; Yildiz, M. Modeling of sensor placement strategy for shape sensing and structural health monitoring of a wing-shaped sandwich panel using inverse finite element method. Sensors 2017, 17, 2775. [Google Scholar] [CrossRef] [PubMed]
  23. Kefal, A.; Oterkus, E.; Tessler, A.; Spangler, J.L. A quadrilateral inverse-shell element with drilling degrees of freedom for shape sensing and structural health monitoring. Eng. Sci. Technol. Int. J. 2016, 19, 1299–1313. [Google Scholar] [CrossRef] [Green Version]
  24. Ko, W.L.; Richards, W.L.; Fleischer, V.T. Applications of the Ko displacement theory to the deformed shape predictions of the doubly-tapered Ikhana wing. In NASA/TP-2009-214652; NASN Dryden Flight Research Center: Edwards, CA, USA, 2009. [Google Scholar]
  25. Jutte, C.V.; Ko, W.L.; Stephens, C.A.; Bakalyar, J.A.; Richards, W.L. Deformed Shape Calculation of a Full-Scale Wing Using Fiber Optic Strain Data from a Ground Loads Test. In NASA Langley Research Center, Rept, TP-215975; NASA Research Center: Hampton, VA, USA, 2011. [Google Scholar]
  26. Derkevorkian, A.; Masri, S.F.; Alvarenga, J.; Boussalis, H.; Bakalyar, J.; Richards, W.L. Strain-based deformation estimating algorithm for control and monitoring applications. AIAA J. 2013, 51, 2231–2240. [Google Scholar] [CrossRef]
  27. Gherlone, M.; Cerracchio, P.; Mattone, M.; Di Sciuva, M.; Tessler, A. Dynamic shape reconstruction of three-dimensional frame structures using the inverse finite element method. In Proceedings of the 3rd ECCOMAS Thematic Conference on Computational Methods in Structural Dynamics and Earthquake Engineering, Corfu, Greece, 26–28 May 2011. [Google Scholar]
  28. Gherlone, M.; Cerracchio, P.; Mattone, M.; Di Sciuva, M.; Tessler, A. Shape sensing of 3D frame structures using an inverse finite element method. Int. J. Solids Struct. 2012, 49, 3100–3112. [Google Scholar] [CrossRef]
  29. Gherlone, M.; Cerracchio, P.; Mattone, M.; Di Sciuva, M.; Tessler, A. An inverse finite element method for beam shape sensing: Theoretical framework and experimental validation. Smart Mater. Struct. 2014, 23, 1–13. [Google Scholar] [CrossRef]
  30. Papadopoulos, M.; Garcia, E. Sensor placement methodologies for dynamic testing. AIAA J. 1998, 36, 256–263. [Google Scholar] [CrossRef]
  31. Li, D.S.; Li, H.N.; Fritzen, C.P. Load dependent sensor placement method: Theory and experimental validation. Mech. Syst. Signal Process. 2012, 31, 217–227. [Google Scholar] [CrossRef]
  32. Rao, A.R.M.; Lakshmi, K.; Krishnakumar, S. A generalized optimal sensor placement technique for structural health monitoring and system identification. Procedia Eng. 2014, 86, 529–538. [Google Scholar]
  33. Kammer, D.C. Optimal sensor placement for modal identification using system-realization methods. Adapt. Struct. Forum 2012, 19, 729–731. [Google Scholar]
  34. Lu, G.; Zhu, X.; Zhang, H.; Gao, Z.; Liu, K. Optimal placement of FBG sensors for reconstruction of flexible plate structures using modal approach. In Proceedings of the 2015 34th Chinese Control Conference (CCC), Hangzhou, China, 28–30 July 2015; pp. 4587–4592. [Google Scholar]
  35. Yang, C.; Lu, Z.; Yang, Z. Robust Optimal Sensor Placement for Uncertain Structures with Interval Parameters. IEEE Sens. J. 2018, 99, 2031–2041. [Google Scholar] [CrossRef]
  36. Yi, X.; Liu, J.; Ye, X.; Lin, S. An optimal placement of FBG sensor network based on probability model. Opt. Int. J. Light Electron Opt. 2013, 124, 1045–1048. [Google Scholar] [CrossRef]
  37. Soman, R.N.; Onoufrioua, T.; Kyriakidesb, M.A.; Votsisc, R.A.; Chrysostomou, C.Z. Multi-type, multi-sensor placement optimization for structural health monitoring of long span bridges. Smart Struct. Syst. 2014, 14, 55–70. [Google Scholar] [CrossRef]
  38. Du, X.; Sun, L.; Xiao, F.; Wang, R. Localization algorithm based on minimum condition number for wireless sensor networks. J. Electron. 2013, 30, 25–32. [Google Scholar] [CrossRef]
  39. Kannan, R.; Hendry, S.; Higham, N.J.; Tisseur, F. Detecting the cause of ill-conditioning in structural finite element models. Comput. Struc. 2014, 133, 79–89. [Google Scholar] [CrossRef]
  40. Shadi, K.; Dehghani, H.L.; Gholampour, I. Well-conditioned sensor placement for range-only localization. In Proceedings of the 2011 4th IFIP International Conference on New Technologies, Mobility and Security, Paris, France, 7–10 February 2011; pp. 1–5. [Google Scholar]
  41. Zhao, Y.; Bao, H.; Duan, X.; Fang, H. The Application Research of Inverse Finite Element Method for Frame Deformation Estimating. Int. J. Aerosp. Eng. 2017, 2017, 1326309. [Google Scholar] [CrossRef]
  42. Tessler, A.; Dong, S.B. On a hierarchy of conforming Timoshenko beam elements. Comput. Struct. 1981, 14, 335–344. [Google Scholar] [CrossRef]
  43. Huh, K.; Stein, J.L. Well-conditioned observer design for observer-based monitoring systems. J. Dyn. Syst. Meas. Control 1995, 117, 592–599. [Google Scholar] [CrossRef]
  44. Werneck, M.M.; Allil, R.C.; Ribeiro, B.A.; de Nazaré, F.V. A Guide to Fiber Bragg Grating Sensors. In Current Trends in Short-and Long-Period Fiber Gratings; INTECH: Vienna, Austria, 2013; pp. 1–24. [Google Scholar]
  45. Opoka, S.; Soman, R.; Mieloszyk, M.; Ostachowicz, W. Damage detection and localization method based on a frequency spectrum change in a scaled tripod model with strain rosettes. Mar. Struct. 2016, 49, 163–179. [Google Scholar] [CrossRef]
  46. Soman, R.; Mieloszyk, M.; Ostachowicz, W. A two-step damage assessment method based on frequency spectrum change in a scaled wind turbine tripod with strain rosettes. Mar. Struct. 2018, 61, 419–433. [Google Scholar] [CrossRef]
Figure 1. Beam kinematic variable fields.
Figure 1. Beam kinematic variable fields.
Sensors 18 02424 g001
Figure 2. Location and coordinate of strain sensor placed on the external beam surface.
Figure 2. Location and coordinate of strain sensor placed on the external beam surface.
Sensors 18 02424 g002
Figure 3. Beam section forces and moments.
Figure 3. Beam section forces and moments.
Sensors 18 02424 g003
Figure 4. Fiber Bragg grating (FBG) sensor (a) and FBG strain rosette (b).
Figure 4. Fiber Bragg grating (FBG) sensor (a) and FBG strain rosette (b).
Sensors 18 02424 g004
Figure 5. Flow chat of solving the optimal model with Particle Swarm Optimization (PSO) algorithm.
Figure 5. Flow chat of solving the optimal model with Particle Swarm Optimization (PSO) algorithm.
Sensors 18 02424 g005
Figure 6. Finite element model of beam and its cross-section division.
Figure 6. Finite element model of beam and its cross-section division.
Sensors 18 02424 g006
Figure 7. Loading cases on cantilever beam: (a) end-node forces: FY = −100 N, FZ = 80 N; (b) uniform distribution loads: QZ = 1.5 N, QY = −1 N; (c) complex loads: uniform loads (QY = 1 N) in direction Y and four equal node forces in direction Z (FZ1 = FZ2 = FZ3 = FZ4 = 150 N).
Figure 7. Loading cases on cantilever beam: (a) end-node forces: FY = −100 N, FZ = 80 N; (b) uniform distribution loads: QZ = 1.5 N, QY = −1 N; (c) complex loads: uniform loads (QY = 1 N) in direction Y and four equal node forces in direction Z (FZ1 = FZ2 = FZ3 = FZ4 = 150 N).
Sensors 18 02424 g007
Figure 8. Wing frame model tests.
Figure 8. Wing frame model tests.
Sensors 18 02424 g008
Figure 9. Comparisons among iFEM calculations for different placements of sensors and NDI measurements for two different loads.
Figure 9. Comparisons among iFEM calculations for different placements of sensors and NDI measurements for two different loads.
Sensors 18 02424 g009aSensors 18 02424 g009b
Table 1. Description of strain sensor configurations. For each configuration, axial locations and orientations of sensors reported as ( x i , θ i , β i ) and angles expressed in degrees.
Table 1. Description of strain sensor configurations. For each configuration, axial locations and orientations of sensors reported as ( x i , θ i , β i ) and angles expressed in degrees.
End-Node LoadingUniformly Distributed Loading
C1C3C2C4
ε 1 (0.5L,−120,0)(0.2L,−120,0)(0.33L,−120,45)(0.2L,−120,0)
ε 2 (0.5L,−120,45)(0.2L,0,0)(0.5L,−120,0)(0.2L,0,0)
ε 3 (0.5L,0,0)(0.2L,120,0)(0.5L,−120,45)(0.2L,120,0)
ε 4 (0.5L,0,45)(0.8L,−120,0)(0.5L,0,0)(0.5L,−120,0)
ε 5 (0.5L,120,0)(0.8L,0,45)(0.5L,0,45)(0.5L,0,0)
ε 6 (0.5L,120,45)(0.8L,120,0)(0.5L,120,0)(0.8L,−120,0)
ε 7 --(0.5L,120,45)(0.8L,0,45)
ε 8 --(0.66L,120,45)(0.8L,120,0)
f ( A ( x , θ , β ) ) 1.0711.451.0910.53
Table 2. Comparisons among deformations computed with inverse Finite Element Method (iFEM) and deformations extracted from ANSYS for end-node load (Loading A).
Table 2. Comparisons among deformations computed with inverse Finite Element Method (iFEM) and deformations extracted from ANSYS for end-node load (Loading A).
Deformation u ( x ) v ( x ) w ( x ) θ ( x ) θ ( y ) θ ( z )
ANSYS0−14.7411.790−0.03−0.03
iFEMC10−14.8811.900−0.03−0.03
AEC100.140.11000
PEC10%0.95%0.93%0%0%0%
iFEMC30−14.6811.740−0.03−0.03
AEC300.060.05000
PEC30%0.41%0.42%0%0%0%
iFEMC20−17.2514.800−0.03−0.03
AEC202.513.01000
PEC20%17%25.5%0%0%0%
iFEMC40−14.6711.740−0.03−0.03
AEC400.070.05000
PEC40%0.47%0.42%0%0%0%
AE = absolute error, PE = percentage error.
Table 3. Comparisons among deformations computed with iFEM and deformations extracted from ANSYS for uniform distribution loads (Loading B).
Table 3. Comparisons among deformations computed with iFEM and deformations extracted from ANSYS for uniform distribution loads (Loading B).
Deformation u ( x ) v ( x )   w ( x )   θ ( x )   θ ( y )   θ ( z )  
ANSYS0−11.1416.710−0.03−0.02
iFEMC10−9.4614.190−0.02−0.02
AEC101.682.5200.010
PEC10%15.08%15.08%0%33.33%0%
iFEMC30−11.216.810−0.03−0.02
AEC300.060.1000
PEC30%0.54%0.6%0%0%0%
iFEMC20−15.5620.070−0.04−0.03
AEC204.423.3600.010.01
PEC20%39.68%20.11%0%33.33%50%
iFEMC40−11.0516.580−0.03−0.02
AEC400.090.13000
PEC40%0.81%0.78%0%0%0%
AE = absolute error, PE = percentage error.
Table 4. Comparisons among deformations computed with iFEM and deformations extracted from ANSYS for complex loads (Loading C).
Table 4. Comparisons among deformations computed with iFEM and deformations extracted from ANSYS for complex loads (Loading C).
Deformation u ( x )   v ( x )   w ( x )   θ ( x )   θ ( y )   θ ( z )  
ANSYS011.1344.350−0.090.02
iFEMC109.4636.040−0.070.02
AEC101.678.3100.020
PEC10%15%18.74%0%22.22%0%
iFEMC3011.2144.410−0.090.02
AEC300.080.06000
PEC30%0.72%0.14%0%0%0%
iFEMC2010.7745.130−0.10.02
AEC200.360.7800.010
PEC20%3.23%1.76%0%11.11%0%
iFEMC4011.0543.910−0.090.02
AEC400.080.44000
PEC40%0.72%0.99%0%0%0%
AE = absolute error, PE = percentage error.
Table 5. Comparisons among deformations computed with iFEM and deformations extracted from ANSYS in case of end-node loads (Loading A). Strain inputs contained errors.
Table 5. Comparisons among deformations computed with iFEM and deformations extracted from ANSYS in case of end-node loads (Loading A). Strain inputs contained errors.
Deformation u ( x )   v ( x )   w ( x )   θ ( x )   θ ( y )   θ ( z )  
ANSYS0−14.7411.790−0.03−0.03
iFEMC10.046.750.370−0.03−0.03
AEC10.0421.4438.58000
PEC1-145.45%327.23%0%0%0%
iFEMC30−17.617.910−0.02−0.04
AEC302.873.8800.010.01
PEC30%19.47%32.91%0%33.33%33.33%
iFEMC20.04−113.7154.580−0.38−0.25
AEC20.0498.96142.7900.350.22
PEC2-671.37%1211.11%0%1166.67%733.33%
iFEMC40.11−16.712.850−0.02−0.04
AEC40.111.961.0600.010.01
PEC4-13.30%8.99%0%33.33%33.33%
AE = absolute error, PE = percentage error.
Table 6. Comparisons among deformations computed with iFEM and deformations extracted from ANSYS in case of uniform distribution loads (Loading B). Strain inputs contained errors.
Table 6. Comparisons among deformations computed with iFEM and deformations extracted from ANSYS in case of uniform distribution loads (Loading B). Strain inputs contained errors.
Deformation u ( x )   v ( x )   w ( x )   θ ( x )   θ ( y )   θ ( z )  
ANSYS0−11.116.60−0.03−0.02
iFEMC105.8940.250−0.03−0.01
AEC1016.9923.65000.01
PEC10%153.06%142.47%0%0%50%
iFEMC30.05−15.6413.560−0.03−0.03
AEC30.054.543.04000.01
PEC3-40.9%18.31%0%0%50%
iFEMC20−93.11118.450−0.29−0.20
AEC2082.01101.8500.260.18
PEC20%738.83%613.55%0%866.67%900%
iFEMC40.16−13.4817.580−0.03−0.03
AEC40.162.380.98000.01
PEC4-21.44%5.9%0%0%50%
AE = absolute error, PE = percentage error.
Table 7. Comparisons among deformations computed with iFEM and deformations extracted from ANSYS in case of complex loads (Loading C). Strain inputs contained errors.
Table 7. Comparisons among deformations computed with iFEM and deformations extracted from ANSYS in case of complex loads (Loading C). Strain inputs contained errors.
Deformation u ( x )   v ( x )   w ( x )   θ ( x )   θ ( y )   θ ( z )  
ANSYS0.011.1344.350−0.090.02
iFEMC10.0654.35−12.940−0.060.01
AEC10.0643.2257.2900.030.01
PEC1-388.32%129.18%0%33.33%50%
iFEMC30.261.0236.47−0.05−0.08−0.01
AEC30.2610.117.880.050.010.03
PEC3-90.84%17.77%-11.11%150%
iFEMC20−125.23342.90−0.87−0.23
AEC20136.36298.5500.780.25
PEC20%1225.16%673.17%0%866.67%1250%
iFEMC40.04−0.4440.950−0.090
AEC40.0411.573.4000.02
PEC4-103.95%7.67%0%0%100%
AE = absolute error, PE = percentage error.
Table 8. Loads for static tests (loads expressed in kg and displacements expressed in mm).
Table 8. Loads for static tests (loads expressed in kg and displacements expressed in mm).
End-Node LoadsUniformly Distributed Loads
Timesiiiiiiivvviiiiiiiivv
loads1.532.223.114.054.985.571.192.663.644.626.02
u ( x ) N D I −0.1−0.2−0.3−0.3−0.4−0.4−0.1−0.1−0.1−0.2−0.2
u ( x ) i F E M C 3 0.0−0.1−0.1−0.2−0.2−0.300−0.1−0.1−0.1
u ( x ) PEC 3 100%50%66.6%33.3%50%25%100%100%0%50%50%
u ( x ) i F E M C 4 0.00.0−0.1−0.2−0.3−0.3000−0.1−0.1
u ( x ) PEC 4 100%100%66.6%33.3%25%25%100%100%100%50%50%
v ( x ) N D I 3.86.19.513.216.118.10.92.53.95.06.4
v ( x ) i F E M C 3 3.75.79.213.917.419.21.02.74.15.26.7
v ( x ) PEC 3 2.6%6.6%3.2%5.3%4.3%5.0%11.1%8%5.1%4%4.7%
v ( x ) i F E M C 4 4.06.49.912.615.317.10.82.23.64.76.0
v ( x ) PEC 4 5.3%4.9%4.2%4.5%4.9%5.5%11.1%12%7.7%6%6.3%
w ( x ) N D I 0.81.32.12.93.74.1−0.10.10.20.30.4
w ( x ) i F E M C 3 0.91.32.23.13.94.300.10.20.20.4
w ( x ) PEC 3 11.3%0%4.8%6.9%5.4%4.9%100%0%0%33.3%0%
w ( x ) i F E M C 4 0.71.22.02.73.43.800.10.10.30.4
w ( x ) PEC 4 11.3%7.7%4.8%6.9%8.1%7.3%100%0%50%0%0%
Table 9. Root-mean-square difference (RMSD) of iFEM calculation for different loads (RMSD expressed in mm).
Table 9. Root-mean-square difference (RMSD) of iFEM calculation for different loads (RMSD expressed in mm).
End-Node LoadsUniformly Distributed Loads
Timesiiiiiiivvviiiiiiiivv
u ( x ) N D I −0.1−0.2−0.3−0.3−0.4−0.4−0.1−0.1−0.1−0.2−0.2
R M S D i F E M C 3 for u(x)0.0680.0970.1510.2050.2770.2960.0180.0370.0620.0790.102
R M S D i F E M C 4 for u(x)0.0690.0980.1520.2070.2800.2980.0180.0380.0620.0800.103
v ( x ) N D I 3.86.19.513.216.118.10.92.53.95.06.4
R M S D i F E M C 3 for v(x)0.1660.1760.2920.5180.3810.4160.1320.1690.2300.2610.308
R M S D i F E M C 4 for v(x)0.1960.2460.3570.5410.6410.6440.1340.1740.2360.2690.320
w ( x ) N D I 0.81.32.12.93.74.1−0.10.10.20.30.4
R M S D i F E M C 3 for w(x)0.0490.0550.0560.1690.1420.1150.0610.1310.1070.1190.155
R M S D i F E M C 4 for w(x)0.090.090.1640.3170.2080.2340.0700.1560.1460.1690.219

Share and Cite

MDPI and ACS Style

Zhao, Y.; Du, J.; Bao, H.; Xu, Q. Optimal Sensor Placement Based on Eigenvalues Analysis for Sensing Deformation of Wing Frame Using iFEM. Sensors 2018, 18, 2424. https://doi.org/10.3390/s18082424

AMA Style

Zhao Y, Du J, Bao H, Xu Q. Optimal Sensor Placement Based on Eigenvalues Analysis for Sensing Deformation of Wing Frame Using iFEM. Sensors. 2018; 18(8):2424. https://doi.org/10.3390/s18082424

Chicago/Turabian Style

Zhao, Yong, Jingli Du, Hong Bao, and Qian Xu. 2018. "Optimal Sensor Placement Based on Eigenvalues Analysis for Sensing Deformation of Wing Frame Using iFEM" Sensors 18, no. 8: 2424. https://doi.org/10.3390/s18082424

APA Style

Zhao, Y., Du, J., Bao, H., & Xu, Q. (2018). Optimal Sensor Placement Based on Eigenvalues Analysis for Sensing Deformation of Wing Frame Using iFEM. Sensors, 18(8), 2424. https://doi.org/10.3390/s18082424

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop