Shape reconstruction of aerospace plate structure is an important issue for safe operation of aerospace vehicles. One way to achieve such reconstruction is by constructing smart fiber Bragg grating (FBG) plate structure with discrete distributed FBG sensor arrays using reconstruction algorithms in which error analysis of reconstruction algorithm is a key link. Considering that traditional error analysis methods can only deal with static data, a new dynamic data error analysis method are proposed based on LMS algorithm for shape reconstruction of smart FBG plate structure. Firstly, smart FBG structure and orthogonal curved network based reconstruction method is introduced. Then, a dynamic error analysis model is proposed for dynamic reconstruction error analysis. Thirdly, the parameter identification is done for the proposed dynamic error analysis model based on least mean square (LMS) algorithm. Finally, an experimental verification platform is constructed and experimental dynamic reconstruction analysis is done. Experimental results show that the dynamic characteristics of the reconstruction performance for plate structure can be obtained accurately based on the proposed dynamic error analysis method. The proposed method can also be used for other data acquisition systems and data processing systems as a general error analysis method.
1. Introduction
Plate structure is a common structure form of aircrafts and aerospace vehicles [1–4]. Shape reconstruction of such structure can be achieved by nonvisual methods using distributed FBG sensors, which generally uses a FBG sensor network or a FBG sensor array [5–7]. Compared with visual perception method [8–11], nonvisual perception method needs less amount of collection data, with high acquisition accuracy, high real-time performance, and high anti-interference capability. Currently there are two types of methods: modal method [5, 12–17] and geometrical method [18–22].
Displacement mode superposition method mentioned in [5, 17, 23, 24] is only suitable for shape reconstruction with small amplitude distortion as the strain direction on the measuring surface is uncertain. Geometrical iteration method mentioned in [19, 25, 26] uses curvature data along a single direction and plane curve fitting algorithm to achieve shape reconstruction. Although these methods have good reconstruction performance for pure bending deformation, they are not suited for torsional deformation. The authors of [26] have achieved curved surface reconstruction, but the method is not suited for real-time reconstruction as it needs to solve complex nonlinear equations. The authors of [19] proposed a space curve shape reconstruction algorithm, but it is not suited for bend torsion coupling shape reconstruction. To overcome the shortcomings of the above algorithms, orthogonal curve net based space surface reconstruction algorithm can be employed by using biorthogonal strain data and orthogonal curve net to achieve shape reconstruction of complex vibration shape. Recursive computation can be employed for real-time implementation instead of solving nonlinear equations.
While performance comparison is mentioned in [5, 12–22, 24], static error analysis such as mean square error between real deformation and theoretically calculated deformation is used [19, 23]. However, there is no dynamic performance comparison method to value these proposed methods to provide transmission performance. While structure vibration reconstruction can be understood as a transformation process of the original strain data, it can be treated as a kind of generalized signal transmission. So dynamic data transmission features can be defined including signal amplitude attenuation, phase delay, and noise mixing [27, 28]. So accurate dynamic data error analysis can be done using system identification methods by taking the vibration response as the input of the reconstruction system and the reconstruction results as the output. Least mean square (LMS) is a famous algorithm for adaptive control [29–31], and many applications for system identification can be found in [32–34].
Taking a one end fixed square flexible plate as experimental plate, a new dynamic data error analysis method is proposed based on LMS algorithm. The proposed method can also be used for other data acquisition systems and data processing systems as a general error analysis method. Section 2 introduces the orthogonal curved network based reconstruction method. Section 3 introduces the construction of the experimental smart FBG plate and the fundamental principles of FBG sensing. Section 4 proposes a dynamic error analysis model for dynamic reconstruction error analysis and gives the detailed modelling method including signal amplitude attenuation, phase delay, and noise mixing. Section 5 illustrates the parameter identification for the proposed dynamic error analysis model and how to design filter F. Section 6 introduces the simulation and experimental platform. Section 7 is the simulation and experiment results. Section 8 comes to the conclusion.
2. Orthogonal Curved Network Based Reconstruction Method
The basic process of curved surface reconstruction algorithm based on orthogonal curved network is as follows: constructing the orthogonal curved network by using continuous curvature and establishing the orthogonal moving coordinate system on two directions for orthogonal curved network; calculating node coordinates of orthogonal curvilinear net using the moving coordinate system and using the internal relations of different nodes to implement the coupling relationship between orthogonal coordinate systems. The iterative calculation of node coordinates and the transformation of moving coordinate system to realize all nodes’ coordinate of orthogonal curved network constitute the two steps used to reconstruct the whole shape of curved surface.
2.1. Establishing Orthogonal Curved Network and Moving Coordinate System
According to differential geometry, a regular parametric curved surface S is a continuous mapping from an area D on space E2 to space E3. Cartesian coordinate system can be established on E2 and E3 separately. Using (u,v) to denote the coordinate in E2 and using (x,y,z) to denote the coordinate in E3,(1)S=fxu,v,yu,v,zu,v.
On the basis of curvature continuity, those measuring points which have equal space (denoted by Δs) between each other were chosen. The measuring points were connecting from u direction and v direction separately. Then the orthogonal equal arc length grid is shown in Figure 1.
Schematic diagram of partial equal arc length grid.
The orthogonal curvature of measuring points is divided into the curvature along u direction and the curvature along v direction, which are denoted by ρu(n,m) and ρv(n,m) separately, where (n,m) is the point ordinal. For convenient description, p(n,m) is used to denote discrete points on curved surface, where n is the ordinal along u direction and m is the ordinal along v direction.
A length-preserving correspondence can be established between deformed curved surface and original curved surface because there is a one-to-one correspondence of every point between pre- and postdeformation, and the orthogonal relationship between two curves through nodes will not change under a micro deformation condition. So, the equal arc length grids mentioned above became an orthogonal curve net as shown in Figure 2.
Orthogonal curve net.
If the coordinate of each node on orthogonal curve net shown in Figure 2 can be solved separately by a certain computation method, the whole smooth curved surface can be obtained by traditional curved surface fitting algorithm. Considering that the research object is one-side fixed and there is no large deformation, there exist two boundary conditions:
There is no deformation on the fixed side of curved surface.
The central line along u direction of orthogonal curve net is plane curve.
In order to describe the algorithm conveniently, the moving coordinate system is established along the u direction and v direction of orthogonal curve net separately as shown in Figure 3.
Moving coordinate system.
To ensure the curve continuity, there is a coordinate transformation relationship between adjacent moving coordinate system on the curve. The relationship contains translation, rotation around z-axis, and rotation around x-axis, as shown in Figure 4.
Moving coordinate system transformation.
In Figure 4, θun,m=ρun,m·Δs, where φu(n,m) is torsion angle. The transformation process is moving the coordinate origin from point p(n,m) to point p(n+1,m) firstly, then rotating θu(n,m) degrees around z-axis, and finally rotating φ(n,m) degrees around x-axis. The torsion angle φu(n,m) is obtained by the curvature coupling relationship between two orthogonal directions.
2.2. Nodal Coordinate Computation
Nodal coordinate recursive computation is the key part of the space curved surface reconstruction algorithm. Because the central line along u direction is plane curve which can be obtained precisely by plate curved fitting algorithm, on the condition of using the central line as the boundary, nodal coordinate of the other curve can be obtained by iterative calculation. The micro arc between point p(n,m) and point p(n+1,m) is circle arc. pn,m·x, p(n,m)·y and p(n,m)·z is used to denote the coordinates of point p(n,m) separately. It is easy to know that, in the moving coordinate system that includes point p(n,m), the coordinate p′(n+1,m) of point p(n+1,m) is(2)p′n+1,m·x=sinθun,mρun,m,p′n+1,m·y=1-cosθun,mρun,m,p′n+1,m·z=0.
Here, p(n+1,m) is the coordinate of point n in absolute coordinate system and p′(n+1,m) is the coordinate of point n in the moving coordinate system. Let point p(1,m) be the first point of the curve. The global coordinate system (absolute coordinate system) is established by using point p(1,m) as the original point. So, the absolute coordinate of point p(2,m) can be obtained by formula (2):(3)p2,m·x=sinθu1,mρu1,m,p2,m·y=1-cosθu1,mρu1,m,p2,m·z=0.
Because the moving coordinate system whose original point is point p(2,m) has been transformed,
rotate θu(1,m) degrees around z-axis,
rotate φu(1,m) degrees around x-axis,
move the original point from point p(1,m) to point p(2,m).
So, if Az(θ,m) denotes the rotation transformation matrix that shows that one point rotates θ degrees around z-axis, the coordinate of point p′(3,m) in absolute coordinate system is(4)p3,m=Azθu1,m·Axφu1,m·p′3,m+p2,m.
Similarly, the moving coordinate system where point p(4,m) is located has been transformed relative to absolute coordinate system as the following steps:
Rotate θu(1,m) degrees around z-axis.
Rotate φu(1,m) degrees around x-axis.
Move the original point from point p(1,m) to point p(2,m).
Continue to rotate θu(2,m) degrees around z-axis.
Continue to rotate φu(2,m) degrees around x-axis.
Move the original point from point p(2,m) to point p(3,m).
From the above transformation, the absolute coordinate of p(4,m) can be obtained:(5)p4,m=Azθu1,m·Axφu1,m·Azθu2,m·Axφu2,m·p′4,m+p′3,m+p2,m.
The absolute coordinate (when n>2) of point p(n) can be obtained:(6)pn,m=∏i=1n-2Azθui,m·Axφui,m·p′n,m+pn-1,m.
In formula (6), Ax(φu(n,m)), Az(θu(n,m)) are(7)Axφn,m=1000cosφun,m-sinφun,m0sinφun,mcosφun,m,Azθn,m=cosθun,m-sinθun,m0sinθun,mcosθun,m0001.
One nodal coordinate can be computed according to formula (3) and formulas (6)~(7). All nodal coordinates can be obtained by iterative recursion combined with coupling transformation of moving coordinate system.
2.3. Coupling Transformation of the Moving Coordinate Systems
The moving coordinate systems along u direction and along v direction are moving along the curve on which they are located. And they converge on the next node of orthogonal curve net. On the basis of a nodal coordinate used to compute the coordinate of next node, coupling transformation between two moving coordinate systems along two directions must be performed. If point pun,m and pvn,m are the points in the moving coordinate system, the midpoint between points pun,m and pvn,m is used as the valid value of node pn,m, as shown in Figure 5. After the node pn,m is obtained, the moving coordinate system along u direction and v direction converges on the node pn,m. Because two curves at node pn,m are orthogonal, corresponding axes should converge. The torsion angle of two moving coordinate systems can be obtained.
Coupling relationship between orthogonal coordinate directions.
Firstly, x-axis and y-axis of the two coordinate systems rotate θu(n-1,m) degrees and θv(n,m-1) degrees around their own z-axis separately. (xu′n-1,m,yu′n-1,m,zu′n-1,m) and (xv′n,m-1,yv′n,m-1,zv′n,m-1) are used to denote the obtained coordinate systems. The vector xu′n-1,m denotes x-axis of the moving coordinate system along u direction and so on. If T(xu′n-1,m,yu′n-1,m,zu′n-1,m)·A is used to denote relative coordinate of vector A in the coordinate system of (xu′n-1,m,yu′n-1,m,zu′n-1,m), the torsion angles of the two moving coordinate system can be obtained as(8)φun-1,m=arctgTxu′n-1,m,yu′n-1,m,zu′n-1,m·xv′n,m-1·yTxu′n-1,m,yu′n-1,m,zu′n-1,m·xv′n,m-1·z,φvn,m-1=arctgTxv′n,m-1,yv′n,m-1,zv′n,m-1·zu′n-1,m·yTxv′n,m-1,yv′n,m-1,zv′n,m-1·zu′n-1,m·x.
The two coordinate systems should rotate γ degrees around their y-axis separately to ensure superposition of corresponding coordinate axis:(9)γun-1,m=-γvn,m-1=arcsinTxu′n-1,m,yu′n-1,m,zu′n-1,m·xv′n,m-1·x2.
Thus, new moving coordinate system along u direction can be obtained:(10)xun,m=Ayγun-1,m·Axφun-1,m·Azθun-1,m·xun-1,m,yun,m=Ayγun-1,m·Axφun-1,m·Azθun-1,m·yun-1,m,zun,m=Ayγun-1,m·Axφun-1,m·Azθun-1,m·zun-1,m.
In the same way, new moving coordinate system along v direction can also be obtained:(11)xvn,m=Ayγvn,m-1·Axφvn,m-1·Azθvn,m-1·xvn,m-1,yvn,m=Ayγvn,m-1·Axφvn,m-1·Azθvn,m-1·yvn,m-1,zvn,m=Ayγvn,m-1·Axφvn,m-1·Azθvn,m-1·zvn,m-1.
In (33) and (34), the matrix of rotation angle α around y-axis is(12)Ayα=cosα0sinα010-sinα0cosα.
Therefore, the coupling transformation of the two moving coordinate systems at the node Puv can be obtained on the basis of formula (10) and (11). The coordinate of next node can be obtained based on the transformed moving coordinate system. All nodal coordinates on orthogonal curve net can be obtained by continuous iterative computation which contains node coordinate recursion and coupling transformation of moving coordinate system.
3. Experimental Smart FBG Plate
When the flexible plate structure deformation occurs, the surface would generate a strain. If the impact of environmental temperature changes is ignored, the strain detected by FBG sensors bonded on the structure surface can be considered as strain of the flexible plate structure on the measuring point:(13)ε=ΔλBKε.
Here, ε is the detected strain, ΔλB is the change of the fiber grating center wavelength, and Kε is the strain sensitivity coefficient. For particular fiber grating, Kε is constant.
While the FBG sensor is employed to detect the curvature in the plane of the plate surface,(14)ρ=2hKεΔλB.
Here, ρ is the curvature of measuring point and h is the thickness of the plate. For a particular plate structure and FBG sensor, h and Kε are constants.
The schematic diagram of a smart FBG plate is shown in Figure 6. A compensated FBG sensor is placed in the bottom of the plate to compensate the temperature influence of the FBG sensors.
Schematic diagram of a smart FBG plate.
Using the detected curvature, orthogonal curved network can be constructed by using continuous curvature. Orthogonal moving coordinate system on two directions for orthogonal curved network can be established. Node coordinates of orthogonal curvilinear net can be calculated using the moving coordinate system, and the internal relations of different nodes can be used to implement the coupling relationship between orthogonal coordinate systems. The iterative calculation of node coordinates and the transformation of moving coordinate system to realize all nodes’ coordinate of orthogonal curved network constitute the two steps used to reconstruct the whole shape of curved surface.
To compare the vibration shape reconstruction performance, laser displacement sensor is employed in most researches. The diagram of the error analysis principle is shown in Figure 7.
Error analysis principle for vibration shape reconstruction of smart FBG structure.
As shown in Figure 7, the reconstruction accuracy validation process is as follows: low-frequency excitation of FBG is applied to structure by the exciter. The distributed FBG sensing network could get the changes of the wavelength of FBG grating measuring point. Then, on the basis of FBG wavelength and the structure strain relationship, these changes can be transferred into strain data. Then orthogonal curved network based reconstruction method can be employed for the structure vibration shape reconstruction. In the same time, a laser displacement sensor is employed to measure the actual displacement of the structure testing point. Then the reconstruction error can be obtained by estimating the difference between the measured and estimated displacement.
4. Dynamic Error Analysis Model
As the essence of all these shape reconstruction algorithms is doing data transformation from the original strain data acquired from the distributed FBG sensors, the dynamic analysis method proposed in signal processing area can be employed. An example of the sampled dynamic data of geometrical reconstruction method can be shown in Figure 8. The solid line is the signal measured by the laser displacement sensor. The dashed line is the estimated signal by geometrical shape reconstruction algorithm. The lag of the two curves is the sampling delay of the two data acquisition systems. The strain data acquisition system is in advance of the laser displacement sensor data acquisition system.
An example of the sampled dynamic data.
Dynamic error analysis can be done by modelling of the data transmission process. The model can be treated as a combination of a FIR filter model with a noise. The diagram of the dynamic error analysis model is shown in Figure 9.
Diagram of the dynamic error analysis model.
For a specific measured point on the plate, xj is the actual displacement of the measured point in the plane of the plate surface. zj is the theoretical reconstructed or estimated displacement without reconstruction noise. nj is the reconstruction noise signal and yj is the sum of zj and nj. The length of the FIR model H is M. zj is the linear weighted sum of all past M inputs:(15)zj=∑i=1Mwixij.
Here, wi is the weighting coefficients, xij is the input sequence, and x1j=xj, x2j=xj-1,…, xMj=xj-M+1. Meanwhile, let W=w1,w2,…,wMT and Xj=x1j,x2j,…,xMj; then formula (15) can be written in vector form:(16)zj=XjW.
LMS algorithm can be employed for model identification of the above dynamic error analysis model. The dynamic modelling principle is shown in Figure 10.
Dynamic modelling principle based on LMS algorithm.
As shown in Figure 10, yj filtered by F gets sj. The output of the identification adaptive filter H∗ is dj. sj minus dj is ej. LMS1 is used to adjust W∗ by minimizing the expectation P of the mean square error of ej. While the adaptive filter converges, W∗ is closer to the series of H and F. Since Wf is known, H2∗ can be built. xj filtered through H2∗and F gets s′j. s′j minus dj gets e′j. LMS2 is used to adjust W2∗ by minimizing the expectation P of the mean square error of e′j. While the adaptive filter converges, W2∗ is almost the same with W. As the main purpose of LMS2 is to get the parameter of H, the exact parameter of H∗ is needed. ej can be used as the starting condition of LMS2 algorithm. Only when the expectation P of the mean square error of ej is less than a certain value can LMS2 start.
With W, the amplitude attenuation, phase delay, and noise mixing can be calculated. In Figure 10,(17)yj=zj+nj,(18)ej=sj-dj,(19)e′j=s′j-dj,(20)sj=YjWf,(21)s′j=XjW2∗Wf,(22)P=Ee2j,(23)P′=Ee′j.
Using σ to describe unit amplitude noise and assuming the amplitude of noise N is S, the expectation PN=E[nj];(24)N=σFj+PN.
With the model parameter W2∗ of H2∗, the amplitude attenuation and phase delay can be obtained. The noise mixing can be obtained using (17) and (24). Therefore, the parameter identification of W2∗ is the key of reconstruction error analysis.
The nomenclature is shown in Nomenclature.
5. Parameter Identification and Reconstruction Error Analysis
To analyze the dynamic reconstruction error, identification of the model parameter W2∗ of H2∗ must be done.
5.1. Parameter Identification
The identification of W2∗ can be done using LMS algorithm. From (18) and (20), one can get(25)P=Ezfj-z∗j2=W∗τEXjTXjW∗-2EzfjXjW∗+Ezfj2.
Let(26)Q=EzfjXj,(27)R=EXjTXj,(28)G=Ezfj2.
Equation (28) can be written as(29)P=W∗TRW∗-2QW∗+G.
Equation (29) is a quadratic function of the weight vector W∗. P~W∗ is a concave hyperparaboloid with a unique minimum. Using the steepest descent method,(30)W∗j+1=W∗j+2uejXj.
Similarly,(31)W2∗j+1=W2∗j+2u′e′jXj.
Here, u, u′ are the adjusting step.
5.2. Reconstruction Error Analysis
Using the obtained model parameter W2∗, the reconstruction error including the amplitude attenuation, phase delay, and noise mixing can be analyzed.
The input signal frequency is f. The sampling period is T, τ=2πfT. The sinusoidal response is(32)yj=∑i=1Nwisinj-i+1τ.
Here ψ is the amplitude of the reconstruction data and Δφ is the phase delay of reconstruction results respected to the actual vibration displacement. Equation (17) shows(39)nj=yj-zj.
One can get(40)nj=yj-XjW.
Using (40), nj can be calculated. The expectation and magnitude of the mixing signal can be obtained.
5.3. Design of Filter F
The mixing noise in the algorithm process could be filtered by filter F. As the mixing noise is high frequency, the filter F should be a low pass filter. The upper limit frequency of the filter F should be greater than main natural frequency of the structure. Filter F can be designed with window function method. Assuming that the amplitude-frequency characteristic of this low-pass filter is Hdejω=1, the phase-frequency characteristic is φ(ω)=0, and then the unit sample response of the filter is(41)hdn=12π∫-ππHdejωejωndω=sinωcnπn.
Here ωc is the upper limit of the low-frequency filter. As hd(n) is an infinite function which is symmetric around hd(0), it is a noncausal system. In order to make hd(n) physically realizable, hd(n) should be truncated with the length M+1: hd(-M/2),…,hd(0),…,hd(M/2). Shifting the truncated hd(n), the following is obtained:(42)hn=hdn-M2,n=0,1,…,M.
System (42) is causal system and has a similar frequency response with hd(n). While the phase frequency response of Hd(ejω) is φ(ω)=-Mω/2, φ(ω) has linear phase. The low-pass filter can be obtained by (42):(43)hdn=sinωcn-M/nπn-M/2.
Using (43) and the actual vibration frequency in the experimental procedure, low-pass filter F can be constructed.
6. Simulation and Experiment Platform6.1. Simulation Platform
The diagram of the simulation software is shown in Figure 11.
The diagram of software structure.
The vibration signal simulation module generates sinusoidal vibration signal. Data processing module is responsible for processing the amplitude, phase, and mixing noise of the data. The processing parameters are set by the parameter module. LMS algorithm module is responsible for dynamic model parameter identification. Data error analysis module is responsible for the analysis and calculation of the amplitude attenuation, phase delay, and noise mixing. The feasibility and effectiveness of dynamic error analysis method proposed in this paper could be verified by comparing the difference of the set parameter and the calculated error parameter.
In order to improve the searching efficiency and accuracy, initial u and u′ are set to 0.01. And the step reduces 10% in every iteration.
Define recognition accuracy as η=(1-|setvalue-identificationvalue|/|setvalue|)·100%. A high recognition accuracy could verify a good dynamic error analysis capability.
6.2. Experiment Platform
A square PC plate is used to construct the smart plate. The elastic modulus is 6.9×109Pa, Poisson’s ratio is 0.33, the density is 1200kg/m3, and the size of the plate is 800×800×5 mm. 25 orthogonal pairs of FBG sensors are bonded on the front and back of the plate. Each pair of the orthogonal FBG sensors can measure the structural strain in horizontal direction and vertical direction, respectively.
With the constructed smart FBG plate, a vibration shape reconstruction experiment platform is built as shown in Figure 12. The experimental platform contains a support base, the smart FBG plate, an excitation system, strain measurement system, displacement measurement system, and shape reconstruction software system. The support base of the experimental platform is an optical experiment table. The table type is GZ103PTB. The excitation system contains a signal generator, a power amplifier, and an exciter. The strain measurement system consisted of the FBG sensors and the fiber grating network analyzer. It can collect the grating wavelength signal of the FBG sensors precisely for strain and curvature calculation. Displacement measurement system consisted of three-dimensional guide, motion controller, motor driver, laser displacement sensor, and data acquisition software. It can achieve real-time acquisition of the vibration displacement of the smart FBG plate. The software system contains a server and a client. The server is used for data transmission of the original data. The client is used for real-time shape reconstruction. The reconstruction result can display dynamically and visually on the screen.
Schematic diagram of the experimental platform.
The photo of the experimental platform is shown in Figure 13. The range of FBG sensor central wavelength is 1532 nm~1568 nm. The type of the fiber grating network analyzer is FONA-2008C. Its precision is 1 pm. The type of the signal generator is SFG-2110. The type of the power amplifier is YE5872. The type of the exciter is JZK-10 with maximum exciting force of 200 N. The model of laser displacement sensor is LK-GD500. The accuracy of the three-dimensional guide is 0.0025 mm.
Photo of the experimental platform.
7. Simulation and Experiment Results7.1. Shape Reconstruction Experiment
The proposed dynamic error analysis method can be used to value the shape reconstruction performance. To test the effectiveness of the proposed error analysis method, geometrical reconstruction method proposed in [19] is adopted in our simulation and experiment analysis. Sweeping frequency technique is used to find the natural frequency of the smart plate. While the vibration frequency is higher than the sixth natural frequency, the maximum vibration displacement of the smart plate is relatively small. So the shape reconstruction in our experiment is only done for the first six natural frequencies. The reconstructed vibration shape at the second natural frequency is shown in Figure 14.
Shape reconstruction for the second modal.
The actual vibration shape
The reconstructed shape
7.2. Simulation Verification
To test the effectiveness of proposed dynamic error analysis algorithm, six different excitations are applied. 2.3 Hz, 4.2 Hz, and 5.9 Hz are the first three resonance frequencies and 1.5 Hz, 3.5 Hz, and 5.0 Hz are nonresonance frequencies. The simulation results are shown in Table 1.
Result of simulation analysis.
Freq.
SA
DA
ADP
SDT
DDT
DDP
SN
DN
NDP
1.5 Hz
96.00%
95.84%
99.83%
0.10T
0.10T
100.00%
0.50 cm
0.57 cm
86.00%
2.3 Hz
85.00%
85.05%
99.94%
0.20T
0.10T
100.00%
0.40 cm
0.46 cm
85.00%
3.5 Hz
99.00%
99.45%
99.55%
0.35T
0.35T
100.00%
0.60 cm
0.62 cm
96.67%
4.2 Hz
90.00%
90.73%
99.19%
0.30T
0.30T
100.00%
0.80 cm
0.92 cm
85.00%
5.0 Hz
80.00%
80.16%
99.80%
0.05T
0.05T
100.00%
0.30 cm
0.34 cm
86.67%
5.9 Hz
98.00%
98.83%
99.15%
0.15T
0.15T
100.00%
0.20 cm
0.23 cm
85.00%
In Table 1, SA represents set attenuation of the output signal. DA represents distinguished attenuation with the error analysis system. ADP represents attenuation distinguished precision of the error analysis method. SDT represents set delay time of the signal phase. DDT represents attenuation delay time of the signal phase. DDP represents delay distinguished precision of the signal phase with the error analysis method. SN represents set noise of the actual output signal. DN represents distinguished noise with the error analysis system. NDP represents noise distinguished precision of the error analysis method. The frequencies of 2.3 Hz, 4.2 Hz, and 5.9 Hz are the first three modal frequencies of FBG smart plate structure. As shown in Table 1, for different vibration frequencies, the identification accuracy of signal attenuation value is not less than 99.15%. The identification accuracy of the signal delay is 100.00%; the identification accuracy of the noise mixing is not less than 85.00%. It can be seen that the algorithm can be set to accurately identify the errors of different amplitude attenuation, phase delay, and noise mixing. The identification accuracy of the amplitude attenuation and phase delay is higher, and noise identification is lower. Figure 15 shows the convergence process of the proposed algorithm. Figure 15(a) is the convergence process of LMS1 algorithm, and Figure 15(b) is the convergence process of LMS2 algorithm. It is shown that the proposed algorithm could rapidly converge in the calculation process, indicating that the algorithm has a good searching capability.
LMS algorithm convergence process.
Convergence of ej
Convergence of e′j
7.3. Experimental Verification
In the vibration experiments, the excitation frequencies are the same as in Table 1. 2.3 Hz, 4.2 Hz, and 5.9 Hz are the first three resonance frequencies and 1.5 Hz, 3.5 Hz, and 5.0 Hz are nonresonance frequencies. The reconstruction experiment results for the first three modals are shown in Figure 16. The performance of the proposed reconstruction algorithm is very good. The fitting curve and the actual curve are nearly the same.
The comparison between actual vibration waveform and reconstruction waveform.
Experimental result at 2.3 Hz
Experimental result at 4.2 Hz
Experimental result at 5.9 Hz
Using traditional static error analysis method, the calculation equation for reconstruction accuracy is shown as follows:(44)η1=∑t=1T1-Et/E′tT·100%.
Here, T is the number of the measurement points and E(t) is the difference between the measuring displacement and reconstruction displacement for a measurement point at the time t. E′(t) is the measuring displacement of the measurement point measured at the time t.
Our proposed method can measure the ratio ρ between the actual vibration amplitude and reconstructed vibration amplitude. The calculation equation of the reconstruction accuracy is(45)η2=1-1-ρ·100%.
As shown in Table 1, η1 and η2 are similar, indicating that our proposed algorithm could also value the static reconstruction performance. But amplitude only gives limited information for the reconstruction performance.
Using the obtained actual displacement data and the reconstruction displacement data, the dynamic error analysis results are shown in Table 2.
Data error analysis of the reconstruction result.
Vibration frequency
η1
η2
Phase delay
Mixing noise amplitude
Mixing noise expectation
1.5 Hz
93.41%
96.37%
20.05 ms
0.54 cm
−0.18 cm
2.3 Hz
87.05%
89.05%
21.07 ms
0.40 cm
−0.10 cm
3.5 Hz
79.94%
79.95%
17.11 ms
0.28 cm
−0.06 cm
4.2 Hz
83.14%
85.49%
16.67 ms
0.19 cm
−0.02 cm
5.0 Hz
78.64%
75.80%
18.04 ms
0.14 cm
0.00 cm
5.9 Hz
73.53%
72.14%
20.30 ms
0.09 cm
0.03 cm
As shown in Table 2, the phase delay of each experiment varies in the range 16.67 ms~20.30 ms, indicating that the phase delay is a fixed delay time value. The delay time is much less than the minimum vibration period, 169.5 ms (at 5.9 Hz), indicating that the reconstruction algorithm has a good following ability. The algorithm can also identify the amplitude and expectation of the mixing noise. As the vibration frequency increases, the vibration amplitude decreases. Meanwhile the mixing noise introduced by the algorithm is also decreased.
8. Conclusion
A new dynamic data error analysis method is proposed in this paper. The detailed analysis process is illustrated, including the construction of the dynamic error analysis model and the LMS algorithm based parameter identification method. Simulation and experiment results show that the proposed dynamic error analysis method can accurately obtain the dynamic characteristics of the FBG smart plate structure reconstruction algorithm including signal amplitude attenuation, phase delay, and noise mixing. It can provide a theoretical basis and technical support for the improvement and engineering realization of various shape reconstruction algorithms. It can also apply to some other applications that need dynamic data analysis.
Nomenclaturexj:
Vibration signal at time j of a testing point on the structure surface
yj:
Reconstructed output
H:
System model
zj:
Theoretical output
nj:
The introduced error
F:
Filter
sj:
Filtered yj
H∗:
The adaptive filter to identify H+F
dj:
Output of H∗
ej:
Difference of sj and dj
W∗:
The identified model of H+F
H2∗:
The estimated model of H
s′j:
Filter output
e′j:
Difference of s′j and dj
W2∗:
Model for H.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This paper is sponsored by program of National Natural Science Foundation of China (no. 51175319), Innovation program of Shanghai Municipal Education Commission (no. 13ZZ075), and Shanghai Key Laboratory of Power Station Automation Technology.
BaneenU.GuivantJ. E.A 2D Bayesian approach for damage detection in plate-type structures201557314415210.1784/insi.2014.57.3.144ZhongY.YuanS.QiuL.Multi-impact source localisation on aircraft composite structure using uniform linear PZT sensors array201511331032010.1080/15732479.2013.8787322-s2.0-84914666558ThambiratnamD. P.ShihH. H.ChanT. H.Structural health monitoring—use of mathematics and computersProceedings of the International Conference on Mathematics and Computers in Sciences and in Industry (MCSI '14)September 2014Varna, Bulgaria555910.1109/mcsi.2014.35WeinzettelJ.ReenaasM.SolliC.HertwichE. G.Life cycle assessment of a floating offshore wind turbine200934374274710.1016/j.renene.2008.04.0042-s2.0-54149111707LiP.LiuY.LengJ.A new deformation monitoring method for a flexible variable camber wing based on fiber Bragg grating sensors201425131644165310.1177/1045389x135102202-s2.0-84906672979WangC.ChengL.Use of fiber Bragg grating sensors for monitoring concrete structures with prestressed near-surface mounted carbon fiber-reinforced polymer strips201425216417310.1177/1045389x134893642-s2.0-84891122874ShapiroY.KosaG.WolfA.Shape tracking of planar hyper-flexible beams via embedded PVDF deflection sensors20141941260126710.1109/tmech.2013.22782512-s2.0-84900507881LiuT.BurnerA. W.JonesT. W.BarrowsD. A.Photogrammetric techniques for aerospace applications20125415810.1016/j.paerosci.2012.03.0022-s2.0-84865463596QiuZ.-C.ZhangX.-T.ZhangX.-M.HanJ.-D.A vision-based vibration sensing and active control for a piezoelectric flexible cantilever plate201410.1177/1077546314536752KiangC. T.SpowageA.YoongC. K.Review of control and sensor system of flexible manipulator201577118721310.1007/s10846-014-0071-42-s2.0-84916888532ParkS. W.ParkH. S.KimJ. H.AdeliH.3D displacement measurement model for health monitoring of structures using a motion capture system20155935236210.1016/j.measurement.2014.09.0632-s2.0-84908433310KimH.-I.KangL.-H.HanJ.-H.Shape estimation with distributed fiber Bragg grating sensors for rotating structures201120303501110.1088/0964-1726/20/3/0350112-s2.0-79952856332LiL.LiW.DingP.ZhuX.SunW.Structural shape reconstruction through modal approach using strain gages2014462Berlin, GermanySpringer273281Communications in Computer and Information Science10.1007/978-3-662-45261-5_29BangH.-J.KimH.-I.LeeK.-S.Measurement of strain and bending deflection of a wind turbine tower using arrayed FBG sensors201213122121212610.1007/s12541-012-0281-22-s2.0-84876566874DerkevorkianA.MasriS. F.AlvarengaJ.BoussalisH.BakalyarJ.RichardsW. L.Strain-based deformation shape-estimation algorithm for control and monitoring applications20135192231224010.2514/1.j0522152-s2.0-84883014497BangH.-J.KoS.-W.JangM.-S.KimH.-I.Shape estimation and health monitoring of wind turbine tower using a FBG sensor arrayProceedings of the IEEE International Instrumentation and Measurement Technology Conference (I2MTC '12)May 2012Graz, AustriaIEEE49650010.1109/i2mtc.2012.62294072-s2.0-84864194187WangZ.-C.GengD.RenW.-X.LiuH.-T.Strain modes based dynamic displacement estimation of beam structures with strain sensors2014231212504510.1088/0964-1726/23/12/1250452-s2.0-84911426494XuH.RenW. X.WangZ. C.Deflection estimation of bending beam structures using fiber bragg grating strain sensors201518339540410.1260/1369-4332.18.3.395YiJ.ZhuX.ZhangH.ShenL.QiaoX.Spatial shape reconstruction using orthogonal fiber Bragg grating sensor array201222667968710.1016/j.mechatronics.2011.10.0052-s2.0-84866083136QiaoX.ZhangH.XuJ.ZhuX.Research and simulation of surface fitting algorithm based on surface patches splicing2012Berlin, GermanySpringer286295Communications in Computer and Information Science10.1007/978-3-642-34390-2_33ChoiK.-S.HuhY.-H.KwonI.-B.YoonD.-J.A tip deflection calculation method for a wind turbine blade using temperature compensated FBG sensors201221202500810.1088/0964-1726/21/2/0250082-s2.0-84856378017LiM.ZhangH.LiuK.ZhuX.Experimental platform design and implementation for plate structure shape perception and reconstruction algorithm2014462Berlin, GermanySpringer199209Communications in Computer and Information Science10.1007/978-3-662-45261-5_21RappS.KangL.-H.HanJ.-H.MuellerU. C.BaierH.Displacement field estimation for a two-dimensional structure using fiber Bragg grating sensors200918202500610.1088/0964-1726/18/2/0250062-s2.0-65349173020KimH.-I.HanJ.-H.BangH.-J.Real-time deformed shape estimation of a wind turbine blade using distributed fiber Bragg grating sensors20141791455146710.1002/we.16442-s2.0-84905823688YiJ.ZhangH.QiaoX.ZhuX.Shape monitoring for wing structure using fiber Bragg grating sensorsProceedings of the IEEE 5th International Conference on Advanced Computational Intelligence (ICACI '12)October 2012Nanjing, China1032103610.1109/icaci.2012.64633292-s2.0-84874617672XuJ.ZhangH.ZhuX.LiL.DingP.Curve surface fitting based on an improved genetic algorithmProceedings of the 6th International Congress on Image and Signal Processing (CISP '13)December 2013Hangzhou, ChinaIEEE74775210.1109/cisp.2013.67452642-s2.0-84897805973SenguptaS.GhoraiS. K.An efficient method for determining transmission characteristics of superstructure fiber Bragg grating and its use for multiparameter sensing20152114615310.1016/j.yofte.2014.11.0092-s2.0-84919436789IslamM. S.MohammadA. B.Transmission characteristics of Mach-Zehnder interferometer comb filter based on thermal operation20141171596510.1007/s00340-014-5798-x2-s2.0-84919877975MohapatraS.KarA.ChandraM.Advanced adaptive mechanisms for active noise control: a technical comparisonProceedings of the International Conference on Medical Imaging, m-Health and Emerging Communication Systems (MedCom '14)November 2014Greater Noida, IndiaIEEE37037510.1109/medcom.2014.7006035XiangM.WeiT.Autobalancing of high-speed rotors suspended by magnetic bearings using LMS adaptive feedforward compensation20142091428143610.1177/1077546313479990MR32257302-s2.0-84903551211RoopaS.NarasimhanS. V.Transform domain Variable Step-size Griffiths Least Mean Square adaptive algorithm and its applications20144041028104110.1016/j.compeleceng.2013.11.0252-s2.0-84904785444SalmanM. S.Sparse leaky-LMS algorithm for system identification and its convergence analysis201428101065107210.1002/acs.2428MR32698612-s2.0-84907937049SuG.JinJ.GuY.WangJ.Performance analysis of l0 norm constraint least mean square algorithm20126052223223510.1109/tsp.2012.2184537MR29542042-s2.0-84859976621RajaM. A. Z.ChaudharyN. I.Two-stage fractional least mean square identification algorithm for parameter estimation of CARMA systems201510732733910.1016/j.sigpro.2014.06.0152-s2.0-84922437835