3.1. Data Collection
To obtain the data for model training and evaluation, in this paper, a self-built GARTEUR aircraft model is utilized for data collection. The GARTEUR aircraft model is a typical and widely used [
34,
35,
36,
37] standard aircraft model with high compliance, low frequency, and dense frequency characteristics. It is designed to evaluate the accuracy and validity for modal testing methods. The GARTEUR aircraft model consists of six aluminum beams with rectangular section, which are the fuselage, wings, vertical tail, horizontal tail, and wing end counterweight plates. The fuselage is 1.5 m long and the wingspan is 2.0 m. The GARTEUR model is made of 2024-T3 aluminum alloy, with the elastic modulus of 73 GPa, the density of 2780 kg/m
and the Poisson’s ratio of 0.33. The total weight of the model is 43.34 kg. Each component is connected by screw connection, and an adapter plate is added to connect the fuselage and the wing.
Figure 3 shows the self-built GARTEUR aircraft model when collecting experimental data. During all experiments, deformation of the GARTEUR aircraft model is within the range of elastic deformation. The system would be a elastic system satisfying Hooke’s law.
To capture the dynamic load and the corresponding response signals, a testing system is built. The entire system is shown in
Figure 3a, which includes the vibration exciters, the power amplifiers, the accelerometers, the force sensors, the data acquisition hardware and software. The data acquisition hardware is LMS SCADAS III from Siemens, which is also utilized for signal synchronization.
For the dynamic load signals, the amplified electrical signals are sent to the exciters to excite the GARTEUR aircraft model and vibrate through the exciter rod. In the system, two exciters are used. Exciter A (
Figure 3a) provides the random dynamic loads. Random dynamic loads from exciter B (
Figure 3a) are considered as the noise load signals for the system. Two high-precision force sensors are placed at the location of the two exciter rods for measuring the loads, respectively. For the corresponding response signals, six high-precision accelerometers are attached to the GARTEUR aircraft model at different locations.
Figure 3b shows the locations of sensors for the dynamic loads, the noise loads, and the responses. The accelerometers are placed at the blue points, which are the nose (R
), the left wing (R
), the right wing (R
), the middle of the fuselage (R
), the rear of the fuselage ((R
)), and the horizontal tail (R
). Acceleration measurements from the sensors are captured and stored via data acquisition hardware. The dynamic load signals, the noise load signals, and the six response signals are temporally synchronized via the data acquisition hardware. The exciters, the power amplifier, the force sensors, and the accelerometers utilized in the system are shown in
Figure 4.
Throughout the experiment, the random signals of both exciters are white noise. The signal generation module of “LMS Test.lab” from Siemens [
38] is utilized to generate the random excitation. The maximum excitation voltage is 0.4 V.
To verify the performance of the proposed method under different noise strengths, the amplitude of load from the exciter B is set to be
,
,
,
,
,
,
,
of that from the exciter A, respectively. For each configuration of noise strength, 16 s of vibration signals is captured. Since the sampling frequency is 2048 Hz, there would be 32,768 data points for each device (exciter or accelerometer) in one noise strength configuration.
Figure 5 shows an example for the dynamic loads, the noise loads, and the responses, with the noise strength of
.
In order to determine the sampling frequency, we conducted a modal experiment. By utilizing PolyMAX [
39], the modal data between 6.062 Hz (mode 1) to 140.024 Hz (mode 14) of the GARTEUR aircraft model are obtained, which is consistent with that of [
40].
Figure 6 shows the modal assurance criterion plot. We choose to set the sampling frequency to 2048 Hz, which is more than 10 times that of the mode 14.
3.3. Experimental Results
To quantitatively evaluate the accuracy and effectiveness of proposed approach, several different metrics are utilized for evaluation. The square root of
, i.e., the
, the correlation coefficient (
R), the relative error (
), the mean absolute error (
), and the mean absolute percentage error (
) are used. Suppose that
denotes the ground truth of the dynamic load, and
is the identified dynamic load.
N denotes the length of vectors
and
. The definitions of the above metrics are given below.
where
and
are the mean and standard deviation for the ground truth of the dynamic load, respectively,
and
those of the identified dynamic load.
where
denotes the
norm of a vector.
where
denotes the absolute value.
The experiments are designed as follows. Firstly, the proposed approaches is quantitatively analyzed for inputs with different strength of load noise. Secondly, inputs with different time windows are quantitatively analyzed. Finally, we compare the propose method with conventional regularization-based methods and NN-based approach with quantitative evaluations.
Experiment 1. Different strengths of load noises.
In this experiment, we verify the performance of the proposed approach for different strengths of noises. Recall that the dynamic loads from exciter B would be regarded as noises. The amplitude of loads from exciter B would thus be the strength of the noise. In this experiment, the noise strengths range varies from
to
.
Table 4 presents the quantitative evaluations on testing data (30% of 16-second data) on the above mentioned metrics. In the table, the ↑ sign next to the metric indicates that the larger the metric, the better the performance of the method, while the ↓ sign means the opposite. From
Table 4, it can be seen that an increase of noise strength would result in a decrease in algorithm performance. However, the correlation coefficient (
R) is relatively high in all rows in
Table 4, which means the identified loads are consistent with those of the ground truth.
Figure 7 gives the comparisons between the identified dynamic loads and the ground truth values under
and
noise strengths for 0.5 s (from 1.5 s to 2 s). Since the noise strengths are different, the ground truth values in the two sub-figures are different. It can be seen from
Table 4 and
Figure 7 that the proposed approach can accurately identify the dynamic loads when the noise strengths are low. When the noises are strong, the proposed approach can estimate a reasonable identification results, e.g.,
when
noise is added. The accuracy of the results in
Table 4 would show the potential application of the proposed approach in real aircraft fault diagnoses.
Experiment 2. Different time windows of input signals.
Recall that the load signals and measurements signals are all captured at 2048 Hz. For the proposed method, measurements from a
time window are inputted to the network, which would correspond to
s. In this experiment, we verify the performances of the proposed approach under different input time windows, ranging from
to
.
Table 5 gives the quantitative results on different metrics under the noise strengths of
,
, and
, representing three cases of no noise, medium strength of noise, and high strength of noise.
It can be seen from
Table 5 that, the size of time window is important to the performances of the algorithm. A larger input time window, which contains more temporal information for the responses, would improve the performances for load identification. However, when the time window is larger than 256, the improvements would be limited. Therefore, in the other experiments, we choose
as the network input. The results of the experiment would show that, for the proposed algorithm and the experimental setup in the paper, information from an input time window of
s would be enough for estimating the load for one time stamp.
Experiment 3. Quantitative comparisons of different approaches.
In this experiment, we quantitatively compare the proposed approach with other methods. We choose conventional load identification method and a data driven method for the comparisons.
For the conventional methods, we adopt the TSVD based and Tikhonov regularization approaches. The two methods are both regularization based methods. The TSVD method would truncate small singular values in the transfer function matrix, avoiding disturbance terms generated from small singular values. The Tikhonov regularization method regards the load identification problem as a minimization problem with a residual term and a weighted
-norm regularization term, which can be explicitly solved. The regularization parameters for both the two approaches are obtained from the GCV method [
41]. The two methods are implemented by Matlab.
For the data driven method, the TDNN [
26] and a wavelet-LSTM [
28] based method are adopted. The TDNN would take the response signals of six measurement points in a time window
as the input. The data is processed by a MLP with one hidden layer. The output would be the dynamic load at time
t. The TDNN is applied for dynamic load identification for aircraft vertical tail model [
26]. In [
26], a 3-layer NN is designed with the input of time window of
, the hidden layer includes 64 neurons. For a reasonable comparison,
N is set to be 256, which is the same as that of the proposed method. The learning rate (
) and the number of epochs (10) are also the same with those of the proposed method. The TDNN is implemented by Pytorch. For the wavelet-LSTM [
28], the input of the network and the output is the eighth-order Meyer wavelet of the signals. The network is with 2 LSTM layer with 128 neurons and 1 FC layer.
Table 6 gives the quantitative comparisons among the four methods. It can be seen that, our methods would outperform in most of the metrics. A higher correlation coefficient (
R) from the propose method shows that a more consistent identified load with the ground truths. From the comparisons with the TDNN method, it can be seen that the designed network architecture with small but deep convolutional kernel and attention modules in the paper would make more use and fuse better the information of the multiple channel inputs. The proposed method can also outperform LSTM based method. In the meantime, our method avoids the explicit estimation of transferring function and manually tuning of the regularization parameters in the two traditional methods, but resulting in a better identification result.
Figure 8 gives the comparisons for different methods with respect to the ground truth values under noise strengths of
for 0.5 s (from 1.5 s to 2 s). In all sub-figures of
Figure 8, the blue curves are the ground truth values for the random loads, where the identified loads using the TSVD method, those of the Tikhonov method, those of the TDNN method and those of the proposed methods, are shown in magenta, green, yellow, and red, respectively.
Figure 8 shows visually that the identified load from the propose method is more consistent with the ground truths, which would be useful for identification problems in real situations.
Experiment 4. Further MISO and MIMO identification performance verification.
In this experiment, we plan to further show the flexibility of the proposed method in terms of dimensions of inputs and outputs. First, by changing the dimensions of input, the proposed method can be used for different MISO configurations. To verify this, we use the 0% noise data for experiment.
Table 7 shows the quantitative evaluation results.
From
Table 7, it can be seen that the increase of input channels can improve the performances, since more information is provide. Moreover, the designed convolutional modules and attention modules in the proposed network can learns the additional information from extra input channels, resulting in an increase of quantitative metrics.
By comparing the first row of
Table 7 with the performance of conventional methods, i.e., the first and second rows of
Table 6, it can be seen that the proposed network with one-channel input would perform similar with the conventional methods. However, the explicitly estimation of transfer function matrix can be avoided in our method.
Second, by changing the dimensions of the network output, the proposed method can be used for MIMO load identification. To verify the flexibility, the six accelerometers’ signals are used for multiple inputs and the load of both Exciter A and B are estimated. Three sets of data are utilized for evaluation, in which the ratios between the amplitude of Exciter A and that of B are 60%, 80%, and 100%, respectively.
Table 8 presents the quantitative analysis results in this experiment. By comparing results in
Table 8 with those in
Table 4, a performance drop is noticed in estimation of load from Exciter A.
The results in this experiment showed that the proposed method would have the flexibility to fit different configuration of SISO, MISO and MIMO load identification problems. Once again, the explicit modeling of the transferring function from the input to the output can be avoided, while a reasonable identification results can be obtained, which would be meaningful in real-situation applications.
Experiment 5. Experiments on additional sensor measurement noise.
In addition to the interference from another load, another potential factor that might affect the performance of the algorithm would be the measurement noise from the accelerometers. It should be noted that the data used in all experiments in the paper are from real sensors, therefore, there are already measurement noises in the accelerometers output. The goal of this experiment is to further analyze the algorithm performance under sensor measurement noise.
The data when only Exciter A is working are used for evaluation to avoid other interferences. A time series Gaussian noise with the mean value of 0 to is added to each of the accelerometer measurements. The standard deviation is set to be
of the amplitude of the accelerometer measurement. In the experiment,
varies from
to
.
Table 9 show the performances of the proposed method on different additional sensor noises. The quantitative results in
Table 9 shows that the performances of the algorithm would decrease with the increase of the additional sensor noise.
To further the evaluate the performances of the proposed method under sensor noise, the above metrics are compared with the conventional methods, namely, the TSVD method and Tikhonov regularization method.
Table 10 shows the quantitative comparison. It can be seen that the drop rate in metric
R is larger for our method than for the conventional methods. The proposed method could outperform at a low noise strength. It should be noted that the main motivation of the proposed method is not to suppress the sensor noise. A further design of the network and fine-tuning of the parameters would be useful for a better performance under high sensor noise.
Experiment 6. Sensibility experiment on different sensor positions.
As the responses of the load would differ at different positions on the aircraft structure, the performance of the proposed method may be affected by the positions of the sensors. In this experiment, our goal is to test the potential sensibility of the algorithm in different sensor positions. To avoid other potential interferences, only one exciter is used in the experiment.
Figure 9 shows the sensor positions in symmetrical and asymmetrical distributions in the experiment.
Table 11 shows the comparisons of the algorithm performance with inputs from different positions of accelerometers. It can be seen that the performance is different. The case of symmetrical distribution would outperform. From the comparisons, it can be seen that the proposed algorithm is sensitive to the sensor position. Since the complexity of the transfer functions for different positions on the aircraft model is different, the accuracy of the implicit estimation from the neural networks of such “transfer functions” would be different given the same amount of input signals. The differences would make the sensibility of the proposed method.