Next Article in Journal
Low-Energy Transfer Design of Heliocentric Formation Using Lunar Swingby on the Example of LISA
Previous Article in Journal
Numerical Study on Tandem-Rotor Autorotation in Forward Flight
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Random Dynamic Load Identification with Noise for Aircraft via Attention Based 1D-CNN

1
Key Laboratory of Civil Aviation Aircraft Airworthiness Certification Technology, Civil Aviation University of China, Tianjin 300300, China
2
School of Safety Science and Engineering, Civil Aviation University of China, Tianjin 300300, China
3
School of Information Engineering, Southwest University of Science and Technology, Mianyang 621010, China
*
Author to whom correspondence should be addressed.
Aerospace 2023, 10(1), 16; https://doi.org/10.3390/aerospace10010016
Submission received: 24 October 2022 / Revised: 10 December 2022 / Accepted: 20 December 2022 / Published: 25 December 2022
(This article belongs to the Section Aeronautics)

Abstract

:
Dynamic load identification plays an important role in the field of fault diagnosis and structural modification design for aircraft. In conventional dynamic load identification approaches, accurate structural modeling is usually needed, which is difficult to obtain for highly nonlinear or unknown structures. In this paper, a one-dimensional convolution neural network with multiple modules is proposed for random dynamic load identification of aircraft. Firstly, the convolution module is designed for temporal feature extraction. Secondly, the extracted features are linearly weighted based on the contributions to the final output. The contributions are learned in a data driven manner via the designed attention module. Lastly, the dynamic load of a certain time stamp is predicted from the learned and weighted features. The proposed model is trained and tested using the real data from a GARTEUR aircraft model. Extensive experimental results with qualitative and quantitative evaluations have demonstrated the identification performance with satisfactory accuracy of the proposed approach under different strengths of load noises.

1. Introduction

With the development of the aviation industry in recent years and the wide application of civil aircraft in public transportation, the research on the safety of aircraft operation has also received more and more attention [1,2,3,4,5]. Usually, the aircraft is working in complex mechanical environments and is subjected to different types of dynamic forces or loads, e.g., the dynamic loads from engine rotor, and the random dynamic loads from the airflow disturbance. It is important to analyze the influence of dynamic load on the aircraft body. The accurate measured external dynamic loads of the aircraft plays an important role in the overall design of the aircraft structure, the strength check and the environmental prediction. However, in practice, limited by factors such as shape design and sensor installation arrangements, it is difficult, if not impossible, to directly measure the external dynamic loads on the aircraft structures. In some cases, the vibration response at certain positions can be measured, and it can be utilized for identifying the external dynamic load information [6]. The load identification can also be used in fault diagnosis. For example, there might be excessive clearance of the aircraft actuator. As the actuator is usually connected with the frame of the aircraft structure, the load can be regarded as a load at some point on the aircraft structure. The identification of such loads would be useful for fault diagnosis and localization. However, in real situations, there could be multiple excessive clearances at different positions, the unknown loads would interference each other, resulting in the increase of the challenge of accurate load identification.
In early research of dynamic load identification, the direct inversion method is often adopted [7]. It builds an algebraic equation between the loads and the responses, and then estimates the dynamic loads from the response by inverting the coefficient matrix in the algebraic equation. While the direct inversion method is simple in principle and easy to apply in practice, it is often difficult to build the equation. For unknown or complicated structures, the problem is often ill-posed, which limits the accuracy of dynamic loads identification. Regularization-based methods improve the ill-posed problem of direct inversion. By adding reasonable boundary conditions to the ill-posed problem, the regularized optimal solution of the equation is sought. Jacquelin et al. [8] provides the L-curve and the general cross-validation (GCV) criterion for determining regularization parameters in truncated singular value decomposition (TSVD) or Tikhonov regularization approaches. Wang et al. [9] propose a regularization-based method with different regularization operators for identification tasks of complex structures. A Bayesian principle-based augmented Tikhonov regularization method is proposed in [10] to improve the load identification accuracy under severe ill-posed conditions. Bianchi et al. [11] proposes a fractional In addition to the TSVD and Tikhonov regularizations, the L 1 norm-based regularization methods for load identification [12,13,14] are also proposed. Li et al. [15] propose a L q norm-based regularization method, in which the q can be learned from a posteriori knowledge.
For the above-mentioned conventional load identification approaches, an accurate dynamic model of the structure is also required. For complex structures, the establishment of an accurate structural dynamic model is often difficult. It can be seen that conventional methods rely too much on the transfer function of the structures. The measurement accuracy of the transfer function, the ill-posedness caused by the rank deficiency of transfer function matrix and the selection of regularization parameters would all affect the accuracy of the identification results.
With the recent rapid developments in artificial intelligence, neural networks (NNs) are now applied in the field of load identification. As a data-driven approach, NN-based methods avoid the need for accurate structural modeling, which could be suitable for load identification of unknown complex structures. In the literature, approaches with different types of NNs [16,17,18,19,20,21] are proposed to identify the static or dynamic load for different structures. Trivailo et al. [16] address the load identification problem of the empennage of a F/A-18 fighter aircraft via a shallow NN. Chen et al. [17] develop a neural network with two hidden layers to solve the impact load conditions for an elastoplastic hemispherical metal shell. Data generated from the finite element method are utilized for network training and testing. Wang et al. [18] propose to use a feed-forward NN to solve the nonlinear identification problem for a one-stage spur gearbox, an extended Kalman filter based algorithm is applied for NN parameters training. Zhou et al. [19] propose a long short-term memory (LSTM)-based recurrent NN model for impact load identification. Hossain et al. [20] compare the performance of radial basis function based NN and multilayer perceptron (MLP) for impact identification of a rectangular plate. A comparison between the NN and LSTM-based identification methods on different types of loads for aircraft is provided in [1].
The main problem of the above-mentioned feed-forward NN or the MLP is that the amount of parameters in the NN is relatively large. To reduce the number of parameters in the NN model, the convolutional neural network (CNN) can be applied. The CNN would process only part of the input data via convolution with a kernel. The parameters of CNN, or the weight of the kernel, are less in quantity compared to those of the full connect feed-forward NNs or the MLP, in which the entire input data are processed. In the literature, CNNs [22] are applied for fault diagnosis [23,24], damage detection [25] based on vibration data in both the time domain [23,25] and the frequency domain [24]. The problems of fault diagnosis or damage detection are to classify the vibration in time or frequency domain into a certain number of classes, e.g., different types of faults. As for the problem of load identification, the goal is to estimate the continuous, time-varying load signal based on the vibration data. Therefore, these problem are different. Xia et al. [26] applied a time delay neural network (TDNN) for dynamic load identification for an aircraft rudder model. The TDNN can be regarded as a special CNN with one layer. Yang et al. [27] proposes to use a two-layer CNN and two-layer MLP to regress the dynamic loads. Yang et al. [28] utilize the hand-crafted wavelet features of signals as an input of a LSTM network and identify the dynamic on a three-degrees-of-freedom vibration system. Models in [26,27] do share some similarities with the proposed model. However, on one hand, Ref. [26] tends to use one layer but large kernels (size of 50) while multiple, small kernels (size of 2) in a multiple layer fashion are utilized in the proposed model. On the other hand, in the proposed model, we tend to further obtain the contributions for the CNN learned features, i.e., the outputs of the convolutions, to the final prediction via a designed attention mechanism and to fuse the output from small convolutions with the learned contributions. These would be the main differences of the proposed approach with methods in [26,27].
To address the load identification problem for unknown structures, especially for random dynamic loads identification problem on the aircraft structure with the existence of the interference of another noise load, in this paper, a novel one-dimensional convolution neural network (1D-CNN) with attention mechanism based approach is proposed. The network would make use of the convolution to extract high-level temporal features of the input multiple synchronized vibration response signals. Considering the fact that different features would contribute differently to the final load prediction, an attention module is designed to learn the contributions for the features from the training data. The dimensions of the input(s) and the output(s) of the proposed approach can be flexible to fit single input single out (SISO), multiple inputs single output (MISO) and multiple inputs multiple outputs (MIMO) load identification problems. Data obtained from a real self-build GARTEUR aircraft model is used for model parameters training and testing. Extensive experimental results show that the proposed approach can accurately identify the random dynamic loads in the presence of noises. Quantitative comparisons shows that the proposed approach outperforms conventional regularization-based approaches and shallow NN based-approaches.
The rest of the paper is structured as follows. In Section 2, we describe the problem of load identification, as well as the proposed deep learning model. Section 3 presents the details of the experimental results. We give the limitations and conclude the paper in Section 4.

2. Proposed Method

In this section, we describe the proposed approach for random dynamic load identification via attention-based 1D-CNN. The problem of load identification is first analyzed. The proposed 1D-CNN model is then described, as well as the implementation details for the network.

2.1. Load Identification Problem

Considering a linear vibration system with p dynamic loads, the corresponding acceleration response y ¨ ( t ) in time domain at a certain measuring point can be calculated from the dynamic loads according to the principle of linear superposition,
y ¨ ( t ) = i = 1 p 0 t h i ( t τ ) f i ( τ ) d τ ,
where h i ( · ) is the unknown response function for the i-th dynamic load at the measuring point, f i the i-th dynamic load. The measurements for acceleration obtained in real experiments are in discrete form. Equation (1) can be discretized by
y ¨ ( t ) = i = 1 p n = 0 N l o a d h i ( n Δ t ) f i ( t n Δ t ) ,
where Δ t is a time interval for discretization, N l o a d = t Δ t and N l o a d N .
As for the problem of load identification, it is to estimate the dynamic load f i in Equation (2) in time domain from the acceleration measurements y ¨ at different times. From Equation (2), it can be seen that the dynamic load at a certain time t, would affect the response y ¨ of one measuring point at time t, t + Δ t , t + 2 Δ t , ⋯. If the vibration system is not linear, the integral equation of Equation (1) would be invalid. The loads can be seen as a function of responses at time t, t + Δ t , t + 2 Δ t , ⋯. More generally, if there are K different measuring points in the system, the dynamic load f i ( t ) can be represented by a function g ( · ) of acceleration measurements from K measuring points,
f i ( t ) = g y ¨ 1 ( t ) , y ¨ 1 ( t + Δ t ) , , y ¨ 2 ( t ) , y ¨ 2 ( t + Δ t ) , , y ¨ K ( t ) , y ¨ K ( t + Δ t ) , ,
where y ¨ k ( t + n Δ t ) , k [ 1 , K ] , n [ 0 , N ] is the response of measuring point k at time t + n Δ t , respectively. The problem of load identification is therefore to estimate the value for f i ( t ) given the response measurements of y ¨ 1 ( t ) , y ¨ 1 ( t + Δ t ) , , y ¨ 2 ( t ) , y ¨ 2 ( t + Δ t ) , , y ¨ K ( t ) , y ¨ K ( t + Δ t ) , .

2.2. The Proposed 1D-CNN-Based Model

Usually, it is difficult to explicitly model the function of g ( · ) in Equation (3) for unknown or complicated structures. Therefore, in this paper, we would utilize the deep NN to model g ( · ) in an implicit and data-driven manner. More specifically, an attention-based 1D-CNN is proposed. Figure 1 shows the framework of the propose approach. The network consists of the convolution modules (blue dotted box in Figure 1), the attention modules (purple dotted box color in Figure 1), a prediction module for final load prediction. The network is designed in an end-to-end manner, the inputs of the network are the vibration signals from K measuring points, and the outputs are the estimated dynamic load. Details of network modules are given below.

2.2.1. Convolution Module

This module would utilize the convolution to extract temporal features of the inputs and generate feature maps for further processing. Considering the inputs of the network are signals of acceleration measurements from K measurement points in the time domain, the 1D convolution is adopted, which is frequently used in signal processing problems. The module includes a 1D convolution layer, a batch normalization layer and an activation function.
1D convolution layer. This layer would take the vibration signals (for the first module) or the feature maps (for the second module) as input. For a discrete input signal f, the output of the 1D convolution is given by
o ( n ) = i = 1 M k ( i ) f ( n + i ) ,
where k is the convolution kernel, M the size of k. The operator ∗ is the cross-correlation. In the network, there would be multiple corresponding kernels for a single 1D input. The number of the kernels is the filter number N f n . The parameters of all kernels would be learned from data. Details for the input and output dimensions are discussed in Section 2.3. To enlarge the receptive field of the convolution, the dilated convolution, which is also called the convolution with holes [29,30,31], is used. The spacing between the values in a kernel is the dilation rate d.
Batch normalization. Considering the data of load identification is in small quantity, to reduce the possible over-fitting, the batch normalization (BN) is applied in this module. In this layer, the output of the dilated convolution layer would be normalized to a standard Gaussian distribution. Given the mean and variance estimated of the output x from the previous layer, i.e., m e a n ( x ) and v a r ( x ) , respectively, the output for the BN is given by
y = B N ( x ) = x m e a n ( x ) v a r ( x ) + ϵ ,
where ϵ is a pre-defined small constant.
The Activation function. The output of the BN layer is then processed via the activation function. In the proposed approach, the Gaussian error linear unit (GELU) function [32] is used. Suppose x is the input, the GELU is given by [32],
G E L U ( x ) = x Φ ( x ) ,
where Φ ( x ) = P ( X x ) , X N ( 0 , 1 ) is the standard Gaussian cumulative distribution function. The GELU is a modified rectified linear unit (ReLU) activation function with nonlinearity for the inputs, and show good performances in speech and vision tasks [32]. The outputs of the activation function are passed to the next module.

2.2.2. Attention Module

The attention mechanism is a special designed module in the NN models to learn and compute the contribution of input data to output data. Inspired by the SE-Net [33], an attention module is designed in the proposed approach. Figure 2 shows the sketch of this module.
Suppose that the input for the attention module is a feature map with the dimension of H × C . The feature maps can be regarded as C different H × 1 features (see Figure 2 left). The attention module tends to weigh the C features using a learned weight, or, the contributions. The contributions of all C features are learnt in a data-driven manner. First, the feature map H × C is down-sampled to 1 × C via two pooling operation. In the “MaxPool” branch, the maximum value of each H × 1 features is selected, while in the “AvgPool” branch, the mean value of each H × 1 features is stored. Secondly, the two 1 × C vector are transposed to C × 1 and connected to two separate fully connected NNs, i.e., the MLPs. The MLP can be regraded as an implicit function. Here, the MLPs are used to map the input into two contribution vectors. The outputs of the two MLPs are summed. The summed vector, which is shown in Figure 2 as a colored vector, is the learned contribution of each H × 1 vector. The vector is called the attention and is used for weighting each H × 1 vector. The weighted H × C feature map is passed to the next module. By weighting the 1D-CNN features using the 1 × C attention vectors, the features could be better fused.

2.2.3. Prediction Module

In this module, the learned feature maps are projected to a scalar, which is the predicted load at a certain time stamp. Suppose that the dimension of the input feature map is H × C . The input is first pooled to 1 × C by average pooling. The 1 × C features are then weighted and summed to 1 × 1 via a fully connected (FC) network, which would be flattened to a scalar. Additionally, the weights in the FC layer are learnable.

2.3. Model Implementation Details

We discuss the details for model implementation in this subsection, including the dimensions of input and output, hyper-parameters of each modules and layers. Suppose that there are K measurement points in the system. The discrete sampled accelerations from a time window, e.g., from t to t + N Δ t , are inputted to the network. We regard the input as N different 1D vectors of K × 1 dimensions. In practice, we set N = 256 , and K = 6 . For the two 1D convolution layers in the network, the size of kernel M, the number of kernels (or the filter number) N f n , the dilation rate d are given in Table 1. The stride s defines the step size of the kernel when traversing the inputs. Considering K is small in our system, we choose to use a small kernel size M. As for the activation function used in the two 1D convolution layers, although we choose the GELU function, the choice for activation function is not limited. The ReLU or Sigmoid function can also be utilized. It is noteworthy that the impact on algorithm performances and metrics with the three activation functions is limited.
As for the attention module, recall that two branches are designed, in which the MLPs are used for learning the attention, respectively. The MLP is a fully connected NN with one hidden layer. The number of nodes, which are also known as the neurons, in the hidden layer are set to be identical for the two branches. Table 2 shows the number of the neurons in the input, hidden, and output layer of the MLPs in both attention modules in the network. The ReLU function is used for activation function in the MLPs.
With the above-mentioned input and hyper-parameters of the network layers, the dimensions of the feature maps of all layers can be determined, which is shown in Table 3. The output of the network would be a scalar, which is the predicted load at time t. It can be seen from Table 3 that two-convolution and two-attention modules are used in the proposed network. It is noteworthy that more copies of the two modules could be added subsequently between the input layer and prediction module, and the total amount of parameters, computational resources for model training and inference are thus increased.
It is noteworthy that the dimension of input and output is not strictly limited. By changing the channels, i.e., K, of inputs, the proposed method can also be utilized for load identification problems with less or more than six channels. As well, the dimensions of the output can also be altered, which can be used in MIMO identification problems. The flexibility of the proposed method in terms of input and output dimensions is verified in the experimental results.

3. Experimental Results

In this section, we describe the details of the experimental results for the proposed approach. The equipment and procedure for data collection are given in Section 3.1. Details of model training are described in Section 3.2. The quantitative evaluations and comparisons are presented in Section 3.3. An ablation study of the proposed approach is also presented in Section 3.4.

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 3 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 # 1 ), the left wing (R # 2 ), the right wing (R # 3 ), the middle of the fuselage (R # 4 ), the rear of the fuselage ((R # 5 )), and the horizontal tail (R # 6 ). 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 0 % , 5 % , 10 % , 20 % , 40 % , 60 % , 80 % , 100 % 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 40 % .
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.2. Model Training

For network training, the mean square error ( M S E ) between the ground truth loads A and identified loads B is used for the loss function, which is defined as
M S E = 1 N i = 1 N ( A i B i ) 2 ,
where N is the amount of data used for training.
We split each 16-second vibration data point into training and testing sets with a proportion of 70 % for training and that of 30 % for testing. The training and testing is done on a platform of a CPU (Intel i5-10500), and a single GPU (Nvidia RTX 3060 Ti with 8 GB memory). For each experiment, the model is trained with 10 epochs, and the total time for model training on the above mentioned platform is about 4 min. The batch size for training is 1024. The Adam optimization is used for network training with a learning rate of 1.0 × 10 4 .

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 M S E , i.e., the R M S E , the correlation coefficient (R), the relative error ( R E ), the mean absolute error ( M A E ), and the mean absolute percentage error ( M A P E ) are used. Suppose that A denotes the ground truth of the dynamic load, and B is the identified dynamic load. N denotes the length of vectors A and B . The definitions of the above metrics are given below.
R = 1 N 1 i = 1 N A i μ A σ A B i μ B σ B ,
where μ A and σ A are the mean and standard deviation for the ground truth of the dynamic load, respectively, μ B and σ B those of the identified dynamic load.
R E = A B 2 A 2 ,
where · denotes the L 2 norm of a vector.
M A E = 1 N i = 1 N A i B i ,
M A P E = 1 N i = 1 N A i B i A i .
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 0 % to 100 % . 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 0 % and 40 % 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., R = 0.7628 when 100 % 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 N = 256 time window are inputted to the network, which would correspond to 0.125 s. In this experiment, we verify the performances of the proposed approach under different input time windows, ranging from N = 64 to N = 384 . Table 5 gives the quantitative results on different metrics under the noise strengths of 0 % , 40 % , and 100 % , 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 N = 256 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 0.125 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 L 2 -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 [ t , t + N Δ t ] 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 [ t , t + 50 Δ t ] , 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 ( 1.0 × 10 4 ) 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 100 % 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 p % of the amplitude of the accelerometer measurement. In the experiment, p % varies from 0 % to 10 % . 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.

3.4. Ablation Study

We perform ablations on the proposed approach to determine the quantitative contribution for load identification of the feature weighting mechanism. The ablations are performed as follows. For the case of “no-att” (no-attention), no attention module would be used in the network, while the dimensions of input, convolution modules, prediction modules would remain. The ablation is performed on the data with 0%, 40%, and 100% noise strength, representing low, mid, and high levels of noises. Table 12 shows the quantitative comparisons on the metrics.
From Table 12, it can be seen that, with the designed feature weighting mechanism, the performance on all metrics would be better than those of the “no-attention” situation. Additionally, the improvement of the metrics would be larger when the noise strengths are larger. The learned attention would help to increase the contribution of more important learned high level features. As for the time for training and testing, for the case of “no-attention”, the time for training is 190 s, and the testing would cost 10 s. For the proposed approach, the time for training is 260 s, and the testing time is 20 s.

4. Conclusions

In this paper, we have presented a novel attention based 1D-CNN model to identify the dynamic loads for aircraft models. As a data-driven approach, the proposed model would avoid the accurate structure modeling. The proposed approach takes advantage of the 1D-CNN to extract temporal features. The learned features are then weighted using learned attentions. The main limitations for the proposed approach is that response data of a certain period of time are needed for model training. In the future, we plan to further analyze the sensibility of sensor positions, while suppressing the potential high sensor measurement noise. A test of the proposed model on real civil aircrafts, while performing a further exploration on the attention mechanism would be other future directions.

Author Contributions

Conceptualization, X.Z.; Data curation, Q.L. and B.X.; Formal analysis, X.Z.; Investigation, B.X.; Methodology, W.H. and Z.F.; Project administration, X.L.; Software, W.H., X.Z. and Q.L.; Supervision, Z.F.; Visualization, Q.L.; Writing—original draft, W.H. and X.Z.; Writing—review and editing, Z.F. and X.L. All authors have read and agreed to the published version of the manuscript.

Funding

This work is supported in part by the Open Fund of Key Laboratory of Civil Aircraft Airworthiness Technology under Grant SH2020112706, in part by the Doctoral Fund of Southwest University of Science and Technology under Grant 19zx7123, in part by the Fundamental Research Funds for the Central Universities (special project of the Civil Aviation University of China) under Grant 3122019163.

Data Availability Statement

The data presented in this study are available on request from the corresponding author.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
GCVgeneral cross-validation
TSVDtruncated singular value decomposition
NNneural network
LSTMlong short term memory
MLPmultilayer perceptron
CNNconventional neural network
TDNNtime delay neural network
1D-CNNone-dimensional convolution neural network
BNbatch normalization
GELUGaussian error linear unit
ReLUrectified linear unit
FCfully connected
MSEmean square error
RMSEroot of mean square error
RErelative error
MAEmean absolute error
MAPEmean absolute percentage error
SISOsingle input single out
MISOmultiple inputs single output
MIMOmultiple inputs multiple outputs

References

  1. Candon, M.; Esposito, M.; Fayek, H.; Levinski, O.; Koschel, S.; Joseph, N.; Carrese, R.; Marzocca, P. Advanced multi-input system identification for next generation aircraft loads monitoring using linear regression, neural networks and deep learning. Mech. Syst. Signal Process. 2022, 171, 108809. [Google Scholar] [CrossRef]
  2. Morse, L.; Cartabia, L.; Mallardo, V. Reliability-based bottom-up manufacturing cost optimisation for composite aircraft structures. Struct. Multidiscip. Optim. 2022, 65, 1–27. [Google Scholar] [CrossRef]
  3. Dangut, M.D.; Jennions, I.K.; King, S.; Skaf, Z. Application of deep reinforcement learning for extremely rare failure prediction in aircraft maintenance. Mech. Syst. Signal Process. 2022, 171, 108873. [Google Scholar] [CrossRef]
  4. Cartocci, N.; Napolitano, M.R.; Costante, G.; Valigi, P.; Fravolini, M.L. Aircraft robust data-driven multiple sensor fault diagnosis based on optimality criteria. Mech. Syst. Signal Process. 2022, 170, 108668. [Google Scholar] [CrossRef]
  5. Agrapart, Q.; Nyssen, F.; Lavazec, D.; Dufrénoy, P.; Batailly, A. Multi-physics numerical simulation of an experimentally predicted rubbing event in aircraft engines. J. Sound Vib. 2019, 460, 114869. [Google Scholar] [CrossRef]
  6. Liu, R.; Dobriban, E.; Hou, Z.; Qian, K. Dynamic load identification for mechanical systems: A review. Arch. Comput. Methods Eng. 2022, 29, 831–863. [Google Scholar] [CrossRef]
  7. Uhl, T. The inverse identification problem and its technical application. Arch. Appl. Mech. 2007, 77, 325–337. [Google Scholar] [CrossRef]
  8. Jacquelin, E.; Bennani, A.; Hamelin, P. Force reconstruction: Analysis and regularization of a deconvolution problem. J. Sound Vib. 2003, 265, 81–107. [Google Scholar] [CrossRef]
  9. Wang, L.; Huang, Y.; Xie, Y.; Du, Y. A new regularization method for dynamic load identification. Sci. Prog. 2020, 103, 0036850420931283. [Google Scholar] [CrossRef]
  10. Jiang, J.; Tang, H.; Mohamed, M.S.; Luo, S.; Chen, J. Augmented tikhonov regularization method for dynamic load identification. Appl. Sci. 2020, 10, 6348. [Google Scholar] [CrossRef]
  11. Bianchi, D.; Buccini, A.; Donatelli, M.; Serra-Capizzano, S. Iterated fractional Tikhonov regularization. Inverse Probl. 2015, 31, 055005. [Google Scholar] [CrossRef] [Green Version]
  12. Bao, Y.; Li, H.; Chen, Z.; Zhang, F.; Guo, A. Sparse l1 optimization-based identification approach for the distribution of moving heavy vehicle loads on cable-stayed bridges. Struct. Control Health Monit. 2016, 23, 144–155. [Google Scholar] [CrossRef]
  13. Pan, C.D.; Yu, L.; Liu, H.L.; Chen, Z.P.; Luo, W.F. Moving force identification based on redundant concatenated dictionary and weighted l1-norm regularization. Mech. Syst. Signal Process. 2018, 98, 32–49. [Google Scholar] [CrossRef]
  14. Liu, J.; Qiao, B.; Chen, Y.; Zhu, Y.; He, W.; Chen, X. Impact force reconstruction and localization using nonconvex overlapping group sparsity. Mech. Syst. Signal Process. 2022, 162, 107983. [Google Scholar] [CrossRef]
  15. Li, Q.; Lu, Q. Time domain force identification based on adaptive Lq regularization. J. Vib. Control 2018, 24, 5610–5626. [Google Scholar] [CrossRef]
  16. Trivailo, P.; Carn, C. The inverse determination of aerodynamic loading from structural response data using neural networks. Inverse Probl. Sci. Eng. 2006, 14, 379–395. [Google Scholar] [CrossRef]
  17. Chen, G.; Li, T.; Chen, Q.; Ren, S.; Wang, C.; Li, S. Application of deep learning neural network to identify collision load conditions based on permanent plastic deformation of shell structures. Comput. Mech. 2019, 64, 435–449. [Google Scholar] [CrossRef]
  18. Wang, J.; Zhang, J.; Wang, Y.; Fang, X.; Zhao, Y. Nonlinear identification of one–stage spur gearbox based on pseudo–linear neural network. Neurocomputing 2018, 308, 75–86. [Google Scholar] [CrossRef]
  19. Zhou, J.; Dong, L.; Guan, W.; Yan, J. Impact load identification of nonlinear structures using deep Recurrent Neural Network. Mech. Syst. Signal Process. 2019, 133, 106292. [Google Scholar] [CrossRef]
  20. Hossain, M.S.; Ong, Z.C.; Ismail, Z.; Khoo, S.Y. A comparative study of vibrational response based impact force localization and quantification using radial basis function network and multilayer perceptron. Expert Syst. Appl. 2017, 85, 87–98. [Google Scholar] [CrossRef]
  21. Gao, Y.; Yu, X.; Chen, L.; Huang, D. Impact Load Identification Algorithm of Helicopter Weapon Pylon Based on Time-Domain Response Signal. Aerospace 2022, 9, 388. [Google Scholar] [CrossRef]
  22. Kiranyaz, S.; Avci, O.; Abdeljaber, O.; Ince, T.; Gabbouj, M.; Inman, D.J. 1D convolutional neural networks and applications: A survey. Mech. Syst. Signal Process. 2021, 151, 107398. [Google Scholar] [CrossRef]
  23. Du, C.; Zhang, X.; Zhong, R.; Li, F.; Yu, F.; Rong, Y.; Gong, Y. Unmanned aerial vehicle rotor fault diagnosis based on interval sampling reconstruction of vibration signals and a one-dimensional convolutional neural network deep learning method. Meas. Sci. Technol. 2022, 33, 065003. [Google Scholar] [CrossRef]
  24. Guo, D.; Zhong, M.; Ji, H.; Liu, Y.; Yang, R. A hybrid feature model and deep learning based fault diagnosis for unmanned aerial vehicle sensors. Neurocomputing 2018, 319, 155–163. [Google Scholar] [CrossRef]
  25. Dang, V.H.; Vu, T.C.; Nguyen, B.D.; Nguyen, Q.H.; Nguyen, T.D. Structural damage detection framework based on graph convolutional network directly using vibration data. Structures 2022, 38, 40–51. [Google Scholar] [CrossRef]
  26. Xia, P.; Yang, T.; Xu, J.; Wang, L.; Yang, Z. Reversed time sequence dynamic load identification method using time delay neural network. Acta Aeronaut. Astronaut. Sin. 2021, 42, 224452. [Google Scholar]
  27. Yang, H.; Jiang, J.; Chen, G.; Zhao, J. Dynamic load identification based on deep convolution neural network. Mech. Syst. Signal Process. 2023, 185, 109757. [Google Scholar] [CrossRef]
  28. Yang, T.; Yang, Z.; Liang, S.; Kang, Z.; You, J. Feature extraction and indentification of stationary random dynamic load using deep neural network. Acta Aeronaut. Astronaut. Sin. 2022, 43, 225952. [Google Scholar]
  29. van den Oord, A.; Dieleman, S.; Zen, H.; Simonyan, K.; Vinyals, O.; Graves, A.; Kalchbrenner, N.; Senior, A.; Kavukcuoglu, K. Wavenet: A generative model for raw audio. arXiv 2016, arXiv:1609.03499. [Google Scholar]
  30. Yu, F.; Koltun, V. Multi-Scale Context Aggregation by Dilated Convolutions. In Proceedings of the International Conference on Learning Representations (ICLR), San Juan, Puerto Rico, 2–4 May 2016. [Google Scholar]
  31. Ran, L.; Zhang, Y.; Hua, G. CANNET: Context aware nonlocal convolutional networks for semantic image segmentation. In Proceedings of the 2015 IEEE International Conference on Image Processing (ICIP), Quebec City, QC, Canada, 27–30 September 2015; pp. 4669–4673. [Google Scholar]
  32. Hendrycks, D.; Gimpel, K. Gaussian Error Linear Units (GELUs). arXiv 2020, arXiv:cs.LG/1606.08415. [Google Scholar]
  33. Hu, J.; Shen, L.; Albanie, S.; Sun, G.; Wu, E. Squeeze-and-Excitation Networks. IEEE Trans. Pattern Anal. Mach. Intell. 2020, 42, 2011–2023. [Google Scholar] [CrossRef] [Green Version]
  34. Balmes, E.; Wright, J.R. GARTEUR group on ground vibration testing: Results from the test of a single structure by 12 laboratories in Europe. In Proceedings of the International Design Engineering Technical Conferences and Computers and Information in Engineering Conference, Sacramento, CA, USA, 14–17 September 1997; American Society of Mechanical Engineers: New York, NY, USA, 1997; Volume 80401, p. V01AT03A004. [Google Scholar]
  35. Chen, J.; Jiang, J.; Wang, K.; Zhang, F. Optimal Placement of Actuators for Active Vibration Control Using EER and Genetic Algorithm. In Proceedings of the 2019 IEEE 10th International Conference on Mechanical and Aerospace Engineering (ICMAE), Brussels, Belgium, 22–25 July 2019; pp. 449–453. [Google Scholar] [CrossRef]
  36. Bao, N.; Wang, C. A Monte Carlo simulation based inverse propagation method for stochastic model updating. Mech. Syst. Signal Process. 2015, 60–61, 928–944. [Google Scholar] [CrossRef]
  37. Rett, S.; Nabarrete, A.; Arbelo, M.; Góes, L.C.; Guimarães, G. Results of the GVT of the Unmodified GARTEUR SM-AG-19 Testbed in South America. In Proceedings of the 51st AIAA/ASME/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference 18th AIAA/ASME/AHS Adaptive Structures Conference 12th, Orlando, FL, USA, 12–15 April 2010; p. 2808. [Google Scholar]
  38. Siemens-Simcenter. Simcenter Testlab, formerly called LMS Test.Lab. 2022. Available online: https://community.sw.siemens.com/s/article/simcenter-testlab (accessed on 10 December 2022).
  39. Peeters, B.; Van der Auweraer, H.; Guillaume, P.; Leuridan, J. The PolyMAX frequency-domain method: A new standard for modal parameter estimation? Shock Vib. 2004, 11, 395–409. [Google Scholar] [CrossRef]
  40. Link, M.; Friswell, M. Working Group 1: Generation of validated structural dynamic models—results of a benchmark study utilising the GARTEUR SM-AG19 test-bed. Mech. Syst. Signal Process. 2003, 17, 9–20. [Google Scholar] [CrossRef]
  41. Hansen, P.C. Regularization tools version 4.0 for Matlab 7.3. Numer. Algorithms 2007, 46, 189–194. [Google Scholar] [CrossRef]
Figure 1. The network architecture of the propose approach. The inputs are K measurements in the time domain. The convolution modules are utilized for feature extraction, while the attention module are used for feature weighting. The final prediction module would map the learned features to the dynamic load for a certain time stamp.
Figure 1. The network architecture of the propose approach. The inputs are K measurements in the time domain. The convolution modules are utilized for feature extraction, while the attention module are used for feature weighting. The final prediction module would map the learned features to the dynamic load for a certain time stamp.
Aerospace 10 00016 g001
Figure 2. The attention module in the proposed approach. The H × C feature map is firstly reduced to one dimension via two pooling. The two pooling 1D vectors are then inputted to two separate fully connected NNs, respectively. The NNs are used for contribution learning. The outputs of the two NNs are summed to obtained the attention vectors, and are used for feature map weighting.
Figure 2. The attention module in the proposed approach. The H × C feature map is firstly reduced to one dimension via two pooling. The two pooling 1D vectors are then inputted to two separate fully connected NNs, respectively. The NNs are used for contribution learning. The outputs of the two NNs are summed to obtained the attention vectors, and are used for feature map weighting.
Aerospace 10 00016 g002
Figure 3. The self-built GARTEUR aircraft model and equipment for data collection. In the system, two exciters with two force sensors, six response accelerometers are used. The system is synchronized by the data acquisition hardware. The positions of two exciter and six acceleration measurement points are also given. (a) The self-built GARTEUR aircraft model used for data collection. (b) The positions of exciters and measurements points (R in blue).
Figure 3. The self-built GARTEUR aircraft model and equipment for data collection. In the system, two exciters with two force sensors, six response accelerometers are used. The system is synchronized by the data acquisition hardware. The positions of two exciter and six acceleration measurement points are also given. (a) The self-built GARTEUR aircraft model used for data collection. (b) The positions of exciters and measurements points (R in blue).
Aerospace 10 00016 g003
Figure 4. Equipment used in the data collection system. The dynamic loads are from the vibration exciter and amplified using the power amplifier. The dynamic loads are measured by the force sensor. Responses from each measurement points are measured by the accelerometers. (a) The vibration exciter. (b) The power amplifier. (c) The force sensor. (d) The accelerometer.
Figure 4. Equipment used in the data collection system. The dynamic loads are from the vibration exciter and amplified using the power amplifier. The dynamic loads are measured by the force sensor. Responses from each measurement points are measured by the accelerometers. (a) The vibration exciter. (b) The power amplifier. (c) The force sensor. (d) The accelerometer.
Aerospace 10 00016 g004
Figure 5. An example of the dynamic loads from exciter A and B, and the corresponding responses from 6 measurement points. Loads from exciter B are considered as the noises, in this example, the noise strength, which is the ratio between the amplitudes of loads from exciter A and those from exciter B, is 40 % . (a) Dynamic loads from exciter A. (b) Dynamic loads from exciter B. (c) Response from measurements R #1. (d) Response from measurements R #2. (e) Response from measurements R #3. (f) Response from measurements R #4. (g) Response from measurements R #5. (h) Response from measurements R #6.
Figure 5. An example of the dynamic loads from exciter A and B, and the corresponding responses from 6 measurement points. Loads from exciter B are considered as the noises, in this example, the noise strength, which is the ratio between the amplitudes of loads from exciter A and those from exciter B, is 40 % . (a) Dynamic loads from exciter A. (b) Dynamic loads from exciter B. (c) Response from measurements R #1. (d) Response from measurements R #2. (e) Response from measurements R #3. (f) Response from measurements R #4. (g) Response from measurements R #5. (h) Response from measurements R #6.
Aerospace 10 00016 g005
Figure 6. The modal assurance criterion plot of 14 modes of the GARTEUR model in our modal experiment. Modal data between 6.062 Hz (mode 1) to 140.024 Hz (mode 14) of the GARTEUR aircraft model are obtained. The sampling frequency is set to 2048 Hz, which is more than 10 times that of mode 14.
Figure 6. The modal assurance criterion plot of 14 modes of the GARTEUR model in our modal experiment. Modal data between 6.062 Hz (mode 1) to 140.024 Hz (mode 14) of the GARTEUR aircraft model are obtained. The sampling frequency is set to 2048 Hz, which is more than 10 times that of mode 14.
Aerospace 10 00016 g006
Figure 7. Comparisons between the identified random loads and the ground truths under different noise strengths. (a) Comparisons with respect to the ground truths (blue) at noise strength of 0% for the proposed method (red). (b) Comparisons with respect to the ground truths (blue) at noise strength of 40% for the proposed method (red).
Figure 7. Comparisons between the identified random loads and the ground truths under different noise strengths. (a) Comparisons with respect to the ground truths (blue) at noise strength of 0% for the proposed method (red). (b) Comparisons with respect to the ground truths (blue) at noise strength of 40% for the proposed method (red).
Aerospace 10 00016 g007
Figure 8. The identified dynamic loads and the ground truth values for different input time windows. Four different time windows are utilized in this experiment. The noise strength is 100 % in this example. (a) Comparisons of the estimated loads using TSVD method (magenta) with respect to the ground truths (blue). (b) Comparisons of the estimated loads using Tikhonov method (green) with respect to the ground truths (blue). (c) Comparisons of the estimated loads using TDNN method (yellow) with respect to the ground truths (blue). (d) Comparisons of the estimated loads using the proposed method (red) with respect to the ground truths (blue).
Figure 8. The identified dynamic loads and the ground truth values for different input time windows. Four different time windows are utilized in this experiment. The noise strength is 100 % in this example. (a) Comparisons of the estimated loads using TSVD method (magenta) with respect to the ground truths (blue). (b) Comparisons of the estimated loads using Tikhonov method (green) with respect to the ground truths (blue). (c) Comparisons of the estimated loads using TDNN method (yellow) with respect to the ground truths (blue). (d) Comparisons of the estimated loads using the proposed method (red) with respect to the ground truths (blue).
Aerospace 10 00016 g008
Figure 9. Symmetrically and asymmetrically distributed accelerometers in the experiment. In total, 4 accelerometers are used. To avoid other interferences, only 1 exciter is used in exepriment. (a) Symmetrically distributed 4 accelerometers and 1 exciter on the aircraft model. (b) Asymmetrically distributed 4 accelerometers and one exciter on the aircraft model.
Figure 9. Symmetrically and asymmetrically distributed accelerometers in the experiment. In total, 4 accelerometers are used. To avoid other interferences, only 1 exciter is used in exepriment. (a) Symmetrically distributed 4 accelerometers and 1 exciter on the aircraft model. (b) Asymmetrically distributed 4 accelerometers and one exciter on the aircraft model.
Aerospace 10 00016 g009
Table 1. Details for the 1D convolutions in the two convolution modules. In the table, “Conv1D #1” means the first 1d convolution module, and “Conv1D #2” the second.
Table 1. Details for the 1D convolutions in the two convolution modules. In the table, “Conv1D #1” means the first 1d convolution module, and “Conv1D #2” the second.
ModuleKernel Size MFilter Number N fn Dilation Rate dStride s
Conv1D #123221
Conv1D #226421
Table 2. Details on the amount of neurons for the MLPs in the two-attention module.
Table 2. Details on the amount of neurons for the MLPs in the two-attention module.
ModuleInput LayerHidden LayerOutput LayerActivation
ATT #15121024512ReLU
ATT #2102420481024ReLU
Table 3. The structure of the proposed network. The dimensions of the input and the output feature map of each layer are listed.
Table 3. The structure of the proposed network. The dimensions of the input and the output feature map of each layer are listed.
ModuleLayerInput ShapeOutput Shape
INPUTInput(6.256)
CONV1D #1Conv1D(6.256)(4.32)
BN(4.32)(4.32)
GELU(4.32)(4.32)
ATT #1Attention(4.32)(4.32)
CONV1D #2Conv1D(4.32)(2.64)
BN(2.64)(2.64)
GELU(2.64)(2.64)
ATT #2Attention(2.64)(2.64)
PREDICTIONAvgPooling(2.64)(1.64)
FC(1.64)(1.1)
Flatten(1.1)1
OUTPUTOutput1
Table 4. Quantitative analysis on different strengths of load noises.
Table 4. Quantitative analysis on different strengths of load noises.
Noise StrengthRMSE ↓ RE ↓R ↑MAE ↓MAPE ↓
0%0.75170.38670.92320.59392.4077
5%0.70370.37630.92720.55443.2146
10%0.76370.38760.92230.59352.1172
20%0.89500.42350.90740.70281.9737
40%1.02610.51440.85990.80272.2846
60%1.23550.59740.80290.95503.0470
80%1.30870.59820.80221.01552.0965
100%1.25170.64790.76280.97743.3292
Table 5. Quantitative analysis on different time windows as inputs.
Table 5. Quantitative analysis on different time windows as inputs.
Noise StrengthNRMSE ↓RE ↓R ↑MAE ↓MAPE ↓
0%641.22660.63520.77390.94222.6046
1280.97090.50280.86470.75732.9802
2560.75170.38670.92320.59392.4077
3840.75990.39180.92020.59542.4893
40%641.32070.67360.74121.00902.6030
1281.24080.63340.77450.96232.8295
2561.02610.51440.85990.80272.2846
3841.02950.51600.85710.80132.8821
100%641.49450.75700.65571.15973.6410
1281.43900.73350.68591.11823.4713
2561.25170.64790.76280.97743.3292
3841.28380.66980.75031.00143.4256
Table 6. Quantitative comparisons on different methods.
Table 6. Quantitative comparisons on different methods.
Noise StrengthMethodRMSE ↓ RE ↓ R ↑MAE ↓MAPE ↓
0%TSVD1.93921.00880.63851.53735.2440
Tikhonov1.69750.88320.64361.33454.1897
TDNN1.55660.80070.65471.22642.5463
wavelet-LSTM0.86430.44760.89670.68672.4251
ours0.75170.38670.92320.59392.4077
40%TSVD2.42331.23300.47881.91867.4397
Tikhonov2.44871.24580.46951.94347.0369
TDNN1.77960.89220.48801.39692.1664
wavelet-LSTM1.10610.56470.83010.87902.2342
ours1.02610.51440.85990.80272.2846
100%TSVD2.79511.41720.35822.24509.1424
Tikhonov2.60191.32300.35852.07038.4509
TDNN1.83330.94950.32331.46492.5837
wavelet-LSTM1.45040.73930.72111.16894.0700
ours1.25170.64790.76280.97743.3292
Table 7. Quantitative comparisons on different numbers of input channels. By changing the dimension of the inputs of the network, the proposed method can fit different MISO configuration.
Table 7. Quantitative comparisons on different numbers of input channels. By changing the dimension of the inputs of the network, the proposed method can fit different MISO configuration.
Noise StrengthChannel NumberRMSE ↓ RE ↓ R ↑MAE ↓MAPE ↓
0%11.50330.77330.64501.16853.7730
21.50520.77430.64471.16843.8007
31.24810.64200.77960.98143.6928
41.19950.61710.80010.92933.4354
50.77190.39710.91870.61162.3188
60.75170.38670.92320.59392.4077
Table 8. Quantitative analysis on the MIMO identification performances. By changing the output dimensions in the network, the load of Exciter A and B can be simultaneously estimated.
Table 8. Quantitative analysis on the MIMO identification performances. By changing the output dimensions in the network, the load of Exciter A and B can be simultaneously estimated.
Amplitude RatioExciterRMSE ↓ RE ↓ R ↑MAE ↓MAPE ↓
60%A1.21120.58570.81220.94112.9255
B1.21350.77420.64020.95953.6395
80%A1.34960.61680.78731.04932.1058
B1.46710.72920.68441.16956.3053
100%A1.30310.67450.73821.01343.1291
B1.75380.71990.69531.97582.8092
Table 9. Quantitative analysis on different strengths of additional sensor noises.
Table 9. Quantitative analysis on different strengths of additional sensor noises.
Sensor Noise StrengthRMSE ↓ RE ↓R ↑MAE ↓MAPE ↓
0%0.75170.38670.92320.59392.4077
2%0.95060.48900.87420.74343.0471
4%1.17630.60510.80000.90993.4878
6%1.33050.68440.73561.04473.4711
8%1.56920.80720.61531.22753.6680
10%1.56080.80290.60701.22344.1296
Table 10. Quantitative comparisons on different strengths of additional sensor noises.
Table 10. Quantitative comparisons on different strengths of additional sensor noises.
Sensor Noise StrengthMethodRMSE ↓RE ↓R ↑ MAE ↓ MAPE ↓
0%TSVD1.50940.78530.61921.16152.9040
Tikhonov1.54480.80370.62231.20413.4852
ours0.75170.38670.92320.59392.4077
6%TSVD1.57300.81830.57491.20473.0625
Tikhonov1.59200.82820.60501.24133.6896
ours1.33050.68440.73561.04473.4711
10%TSVD1.59600.83030.55751.22652.9476
Tikhonov1.70990.88960.56581.32764.1587
ours1.56080.80290.60701.22344.1296
Table 11. Quantitative comparisons on different distribution of sensors.
Table 11. Quantitative comparisons on different distribution of sensors.
Accelerometers DistributionRMSE ↓RE ↓R ↑MAE ↓ MAPE ↓
Symmetrical1.19950.61710.80010.92933.4354
Asymmetrical2.06090.70900.70791.65352.7237
Table 12. Ablations on different noise strengths.
Table 12. Ablations on different noise strengths.
Noise StrengthMethodRMSE ↓ RE ↓ R ↑MAE ↓ MAPE ↓
0%no-att0.78270.40260.91610.62222.5564
ours0.75170.38670.92320.59392.4077
40%no-att1.06600.53440.85310.83842.2758
ours1.02610.51440.85990.80272.2846
100%no-att1.30400.67500.73801.02022.8646
ours1.25170.64790.76280.97743.3292
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

He, W.; Zhang, X.; Feng, Z.; Leng, Q.; Xu, B.; Li, X. Random Dynamic Load Identification with Noise for Aircraft via Attention Based 1D-CNN. Aerospace 2023, 10, 16. https://doi.org/10.3390/aerospace10010016

AMA Style

He W, Zhang X, Feng Z, Leng Q, Xu B, Li X. Random Dynamic Load Identification with Noise for Aircraft via Attention Based 1D-CNN. Aerospace. 2023; 10(1):16. https://doi.org/10.3390/aerospace10010016

Chicago/Turabian Style

He, Wenbo, Xiaoqiang Zhang, Zhenyu Feng, Qiqi Leng, Bufeng Xu, and Xinmin Li. 2023. "Random Dynamic Load Identification with Noise for Aircraft via Attention Based 1D-CNN" Aerospace 10, no. 1: 16. https://doi.org/10.3390/aerospace10010016

APA Style

He, W., Zhang, X., Feng, Z., Leng, Q., Xu, B., & Li, X. (2023). Random Dynamic Load Identification with Noise for Aircraft via Attention Based 1D-CNN. Aerospace, 10(1), 16. https://doi.org/10.3390/aerospace10010016

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