Next Article in Journal
Iterants, Majorana Fermions and the Majorana-Dirac Equation
Next Article in Special Issue
Multifrequency Topological Derivative Approach to Inverse Scattering Problems in Attenuating Media
Previous Article in Journal
Advantages of the Discrete Stochastic Arithmetic to Validate the Results of the Taylor Expansion Method to Solve the Generalized Abel’s Integral Equation
Previous Article in Special Issue
Application of Anti-Diagonal Averaging in Response Reconstruction
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Online Dynamic Load Identification Based on Extended Kalman Filter for Structures with Varying Parameters

1
Mechatronic Engineering College, Jinling Institute of Technology, Nanjing 211169, China
2
State Key Laboratory of Mechanics and Control of Mechanical Structures, Nanjing University of Aeronautics and Astronautics, Nanjing 210016, China
3
Institute for Infrastructure and Environment, Heriot-Watt University, Edinburgh EH14 4AS, UK
*
Authors to whom correspondence should be addressed.
Symmetry 2021, 13(8), 1372; https://doi.org/10.3390/sym13081372
Submission received: 15 June 2021 / Revised: 12 July 2021 / Accepted: 25 July 2021 / Published: 28 July 2021
(This article belongs to the Special Issue Advanced Mathematical and Simulation Methods for Inverse Problems)

Abstract

:
Dynamic load identification is an inverse problem concerned with finding the load applied on a structure when the dynamic characteristics and the response of the structure are known. In engineering applications, some of the structure parameters such as the mass or the stiffness may be unknown and/or may change in time. In this paper, an online dynamic load identification algorithm based on an extended Kalman filter is proposed. The algorithm not only identifies the load by measuring the structural response but also identifies the unknown structure parameters and tracks their changes. We discuss the proposed algorithm for the cases when the unknown parameters are the stiffness or the mass coefficients. Furthermore, for a system with many degrees of freedom and to achieve online computations, we implement the model reduction theory. Thus, we reduce the number of degrees of freedom in the resulting symmetric system before applying the proposed extended Kalman filter algorithm. The algorithm is used to recover the dynamic loads in three numerical examples. It is also used to identify the dynamic load in a lab experiment for a structure with varying parameters. The simulations and the experimental results show that the proposed algorithm is effective and can simultaneously identify the parameters and any changes in them as well as the applied dynamic load.

1. Introduction

In many applications related to structural safety and reliability, it is necessary to obtain the dynamic loads active on a structure. Furthermore, it is possible that the dynamic parameters of such a structure vary with time. In such cases, especially for gradually varying parameters, the dynamic characteristics of a structure cannot be known in advance. Therefore, in order to identify the load acting on such structures, it is important to identify the structural parameters before or simultaneously as we identify the load by solving an inverse problem starting from the structural responses. Such information will not only be useful for monitoring the loads but also to monitor the characteristics of a structure as it changes in real time.
The methods of parameters’ identification can be divided into time-domain and frequency-domain methods [1,2,3,4]. In the frequency-domain, the methods are in general concerned with identifying the modal parameters of a structure, and then using the modal coordinate transformation relationship to obtain the physical parameters [5,6,7]. As for the time-domain, several scholars have proposed a variety of identification methods for the modal parameters as well as the physical parameters. The proposed methods include: Random Decrement Technique (RDT) [8], Ibrahim Time Domain Method (ITD) [9], Stochastic Subspace Identification (SSI) [10], Eigensystem Realization Algorithm (ERA) [11] and Natural Excitation Technique (NExT) [12,13]. However, using these methods requires the excitation to be white noise, or the response to be a free decay process. There are also extended Kalman filter (EKF) [14], least Square Estimation (LSE) [15], H filter [16] and Monte Carlo filter [17] to identify the physical parameters for a structure in the time-domain.
Similarly, the dynamic load identification can also be classified into frequency- and time-domain methods. The frequency-domain methods are relatively more mature, but they have certain limitations. The identification accuracy strongly depends on the time length of the response data acquisition. Their performance is negatively impacted if the response data are only available for a short time period. The time-domain methods are relatively recent compared to the frequency-domain. In general, time-domain methods express the system characteristics as a Duhamel integral, and then obtain the load acting on a structure using a convolution response. Eills [18] is the first to identify the aerodynamic load acting on a structure using the frequency-domain. Bartlett [19] identified the load on the hub by measuring the acceleration response of a helicopter. Inoue [20] used deconvolution method to identify the load in the frequency-domain. Milana [21] proposed a method of load identification by coherent analysis of a structural response. Liu et al. [22] used an interpolation function to fit the load to be identified. Zhang [23] deduced the method of identifying two dimensional distributed dynamic loads in generalized orthogonal domains. By introducing the Tikhonov regularization algorithm, Jiang [24,25] improved the accuracy of identifying distributed random dynamic loads in the generalized orthogonal polynomial domain. Hory [26] first used the method of time discretization to identify aircraft flight loads. Wu et al. [27] identified the load in time domain based on Green’s function. Gunawan [28] improved the recognition effect of the impact load with the help of a L-curve and the truncated singular value method. Song [29] proposed a dynamic load identification method based on the Kalman filter, which took the strain signal as the observation signal and identified the unknown load on a cantilever beam through recursive iteration. The above methods can only identify the excitation when the dynamic characteristics of the structure are known.
Starting in the past decade, a lot of efforts were invested in solving the inverse problem of identifying the loads when the structural parameters are unknown while the dynamic model is given. Based on the classical Kalman filter (KF) [30], Yang [31] proposed an adaptive Kalman filtering algorithm. T. Lenkutis [32] modified the sinusoidal approximation using a windowing function, which can remove sharp jumps in a generated road profile. Y. Lei et al. [33] used an extended Kalman filter to identify the load and the state vectors through multiple stages. The algorithm can identify the dynamic load when the structural parameters are unknown and track changes in the structural parameters. Naets et al. [34] reconstituted the structural state vector parameters and excitation into an augmented vector, and then used the extended Kalman filter to identify the load and the unknown parameters. Du et al. [35] studied the identification of the blade tip vortex and the motion characteristics of the vortex. However, the “Drift phenomenon” is observed in the displacement and in the load identification if the measured acceleration is contaminated with noise. To solve this problem, Huang et al. [36] proposed a data fusion method, which used partial data fusion of displacement and acceleration responses in the observation equation. Lei et al. [37] proposed a dynamic load identification algorithm based on the recursive least square method. Zhang and He [38] introduced the projection matrix to identify jointly structural parameters and unknown loads. Yu [39] combined the theory of the Hilbert transform and Kalman filter to identify the parameters and the dynamic load of a shear frame structure.
In this paper, we are concerned with the dynamic load identification where structural parameters are unknown and change gradually in time. An example of such a system is a launched rocket where the mass is reduced in time due to the fuel consumption as well as moving through stages of a multistage rocket. Standard load identification methods cannot identify the load acting on such systems. We propose a load identification algorithm based on an extended Kalman filter. We validate the algorithm using numerical and experimental studies where we also consider data contaminated with noise. The algorithm is used to identify properties of structures with gradually varying stiffness or with gradually varying mass. To deal with continuous structures, they are first discretized using the finite element method. The discretization produces a symmetric linear system of equations. If the system includes many degrees of freedom, we implement the model reduction theory to reduce the number of degrees of freedom before applying the proposed extended Kalman filter algorithm. The rest of this paper is organized as follows: the load identification method is presented in Section 2; in Section 3, we discuss implementing the model reduction theory to reduce the number of degrees of freedom; in Section 4, three numerical examples are presented to study the efficiency and the accuracy of the proposed algorithm using different structural parameters and different noise conditions, while in Section 5, the proposed algorithm is used to identify the dynamic load of a simply supported beam with gradually varying mass characteristic in an experimental setting; we finish in Section 6 with some concluding remarks.

2. Inverse Algorithm

The general vibration differential equation of a n-DOF system is written as
M p ¨ t + C p ˙ t + K p t = B u f t .  
Here, the parameters M ,   C and K represent the mass, the damping and the stiffness symmetric matrices of n × n for a given structure, while the variables p ¨ ( t ) , p ˙ t and p t are acceleration, velocity and displacement vectors of n × 1, respectively.
Starting form Equation (1) and employing x t = p t p ˙ t to describe the state vector of the structure, a linear state space of the vibration differential equation can be obtained with
x ˙ t = 0 I M 1 K M 1 C x t + 0 M 1 B u f t   .
where the influence matrix B u is the n × s matrix which is composed of 0 and 1 to represent the location of the applied excitation so that the position of the excitation on the structure is allocated 1 and otherwise is allocated 0. The matrix 0 is a zero matrix. The variable f t is the excitation vector of s × 1.
If some of the parameters in the mass or the stiffness or the damping matrices are unknown, we combine these unknown parameters with the state vector x t to form the extended state vector z ( t ) in Equation (3)
z t = x t α = p t p ˙ t α .
Therefore, the state updating equation of the extended state vector can be rewritten as
z ˙ t = p ˙ t M 1 K p t M 1 C p ˙ t + M 1 B u f t 0 = f c z t , f t .  
Clearly, the state updating equation with the extended state vector is a nonlinear equation.
The form of measurement available in the measurement updating equation is related to the type of observation that can be gathered from the structure. Usually, the acceleration values measured at the degrees of freedom are the easiest to get. The measurement updated equation can be written as
y t = H 0 M 1 K p t C p ˙ t + B u f t = h z t , f t .  
Thus, the state updating equation and the measurement updating equation are nonlinear equations with respect to two variables. Applying the Taylor expansion, the above two nonlinear equations are transformed into a sum of polynomials. We retain only the first-order polynomials which results into the following
f c z t , f t f c z k 1 | k 1 , f ^ k 1 + z f k 1 c · z t z k 1 | k 1 + f f k 1 c · f t f ^ k 1 ,
h z t , f t h z k | k 1 , f ^ k 1 + z h k · z t z k | k 1 + f h k · f t f ^ k 1 .
where z f k 1 c and z h k are the first-order partial derivatives of the state updating equation and the measurement updating equation to the extended state vector z t , respectively. Then, f f k 1 c and f h k are the first-order partial derivatives of the state updating equation and the measurement updating equation to the load vector. The subscript k or k 1 denotes the time step. The parameter z k | k 1 indicates the prior estimation of the extended state vector at the k th time step starting from the previous time step, namely, k 1 .
In order to identify the extended state vector and the unknown load, the extended Kalman filter can be divided into three steps: state updating, excitation identification and measurement updating, which is described below.
This starts with the following equations
z k | k 1 = z k 1 | k 1 + t k 1 t k f c z k 1 | k 1 , f ^ k 1 d t ,
f ^ k = J k y k h k z k | k 1 , f ^ k 1 + f h k · f ^ k 1 ,
z k | k = z k | k 1 + K k y k h z k | k 1 , f ^ k 1 f h k · f ^ k f ^ k 1 .
The matrices J k and K k need to be determined in order to estimate the state z k | k and the excitation f ^ k as the minimum variance unbiased estimation of the system state z k and the real excitation f k .
(1)
Step 1: State updating
The recursive relation of the extended state vector can be written as
z k = z k 1 + t k 1 t k f c z t , f t d t + w = z k 1 + Δ t · [ f c z k 1 | k 1 , f ^ k 1 + z f k 1 c · z k 1 z k 1 | k 1 + f f k 1 c · f k 1 f ^ k 1 ] + w .
Then, the prior estimation of the extended state vector can be written as
z k | k 1 = z k 1 | k 1 + t k 1 t k f c z k 1 | k 1 , f k 1 d t   = z k 1 | k 1 + Δ t · f c z k 1 | k 1 , f k 1 .
Thus, the prior estimation error is
z ˜ k | k 1 = z k z k | k 1 = z k 1 z k 1 | k 1 + Δ t · z f k 1 c · z k 1 z k 1 | k 1 + f f k 1 c · f k 1 f ^ k 1 + w = 1 + Δ t · z f k 1 c · z k 1 z k 1 | k 1 + Δ t · f f k 1 c · f k 1 f ^ k 1 + w = 1 + Δ t · z f k 1 c · z ˜ k 1 | k 1 + Δ t · f f k 1 c · f ˜ k 1 + w .
where z ˜ k 1 | k 1 is the posterior estimation error at time ( k 1 ) Δ t , and f ˜ k 1 is the load estimation error at time ( k 1 ) Δ t .
We can obtain the variance of prior estimation error using
P k | k 1 z = E ( z ˜ k | k 1 · z ˜ k | k 1 ) T = 1 + Δ t · z f k 1 c Δ t · f f k 1 c P k 1 | k 1 z P k 1 z f P k 1 f z P k 1 f ( 1 + Δ t · z f k 1 c ) T ( Δ t · f f k 1 c ) T + Q .
with P k 1 | k 1 z being the variance matrix of a posterior error estimation at the time instant ( k 1 ) Δ t , while P k 1 f is the variance matrix of the load error at the same time instant. The values of P k 1 z f and P k 1 f z are the cross covariance matrix of the extended state vector and load:
P k 1 | k 1 z = E ( z ˜ k 1 | k 1 · z ˜ k 1 | k 1 ) T ,
P k 1 f = E ( f ˜ k 1 · f ˜ k 1 ) T ,
P k 1 z f = ( P k 1 f z ) T = E ( z ˜ k 1 | k 1 · f ˜ k 1 ) T .
(2)
Step 2: Excitation identification
Firstly, we need to establish the equation about the load to be solved:
y ˜ k = y k h k z k | k 1 , f ^ k 1 + f h k · f ^ k 1 .
Then, by substituting Equation (7) into the above Equation (18), we get
y ˜ k = f h k · f k + z h k · z k z k | k 1 + v = f h k · f k + z h k · z ˜ k | k 1 + v = f h k · f k + e ˜
where e ˜ = z h k · z ˜ k | k 1 + v , and the mean value of e ˜ is 0. As the components of e ˜ have different variances, it is necessary to select the appropriate weighting matrix W to optimize the estimation result. Then,
f ^ k = [ ( f h k ) T · W · f h k ] 1 ( f h k ) T W · y ˜ k ,
R ˜ k = E e ˜ k e ˜ k T = ( z h k ) T · P k | k 1 z · z h k + R k .
where W is the inverse matrix of the variance matrix R ˜ k , and R ˜ k is the variance matrix of e ˜ . The matrix J k is
J k = [ ( f h k ) T R ˜ k 1 f h k ] 1 ( f h k ) T R ˜ k 1 .
In this case, the variance matrix of the error of the excitation estimation f ^ k is as follows:
P k f = [ ( f h k ) T · R ˜ k 1 · f h k ] 1 .
(3)
Step 3: Measurement updating
After determining the matrix J k , we also need to determine the Kalman gain matrix K k . Substituting Equation (23) into Equation (11), we obtain the posterior estimation z k | k :
z k | k = z k | k 1 + K k [ y k h z k | k 1 , f ^ k 1   f h k · J k y k h k z k | k 1 , f ^ k 1 + f h k · f ^ k 1 f ^ k 1 ]   = z k | k 1 + K k I f h k · J k y k h z k | k 1 , f ^ k 1 + f h k · f ^ k 1 .
Thus, the error of posterior estimation z ˜ k | k is
z ˜ k | k = z k z k | k = z k z k | k 1 K k I f h k · J k y k h z k | k 1 , f ^ k 1 + f h k · f ^ k 1 = z k z k | k 1 K k I f h k · J k z h k · z k z k | k 1 + f h k · f k + v = I K k I f h k · J k z h k · z ˜ k | k 1 + K k I f h k · J k f h k · f k + K k I f h k · J k · v .
Considering z k | k 1 is an unbiased estimation, in order to ensure that z k | k is also an unbiased estimation, the following equation must be satisfied:
K k I f h k · J k f h k = 0 .
Therefore, the variance of posterior estimation can be written as
P k | k z = E ( z ˜ k | k · z ˜ k | k ) T = I K k I f h k · J k z h k · P k | k 1 z · I K k I f h k · J k z h k ] T + K k I f h k · J k · R · K k I f h k · J k ] T   = I L k z h k · P k | k 1 z · ( I L k z h k ) T + L k · R · L k T .
Among them,
L k = K k I f h k · J k .
In order to make the posterior estimation optimal, the variance matrix P k | k z of the posterior error estimation should be minimized:
P k | k z K k = 0 .
According to the above Equation (29), the Kalman filter gain K k can be obtained as follows:
K k = P k | k 1 z ( z h k ) T R ˜ k 1 .  
By substituting the above Equation (30) into Equation (27), the variance matrix of posterior error estimation is rewritten as
P k | k z = P k | k 1 z K k · ( R ˜ k f h k · P k f · f h k ) T · K k .  
Then, we can get
P k z f = P k | k 1 z ( z h k ) T J k T = K k · f h k · P k f .  
In order to reduce the numbers of the given initial conditions, the order of the load identification algorithm based on the extended Kalman filter can be adjusted including all the above three steps. It should be stressed that the order adjustment does not affect the efficiency of the algorithm. It only involves a different initial value which has little influence on the whole process of load identification.
In this paper, we distinguish between two cases:
  • the unknown parameters are the stiffness or the damping coefficients
  • the unknown parameter are the mass coefficients.
First, we discuss the case of unknown stiffness or damping. We can get the partial derivative matrix of the state update equation and the measurement update equation to the extended state vector based on Equations (4) and (5) such as
z f k 1 c = f c p   f c p ˙   f c α z = z k 1 k 1 = M 1 0 I 0 0 K C K α 1 p k 1 k 1 C α 1 p ˙ k 1 k 1 K α m p k 1 k 1 C α m p ˙ k 1 k 1 0 0 0 0  
z h k = h p h p ˙ h α z = z k k 1 = M 1 H 0 K C K α 1 p k k 1 C α 1 p ˙ k k 1 K α s p k k 1 C α s p ˙ k k 1
The partial derivative matrices of the state update equation and the measurement update equation with respect to the load vector are
f f k 1 c = 0 M 1 B u 0 ,
f h k = H 0 M 1 B u .
In this case, f f k 1 c   and   f h k are constant matrices. The measurement update step and the load identification step can be simplified as
z k | k = z k | k 1 + K k y k h ˉ z k | k 1 f h k · f ^ k ,
f ^ k = J k y k h ˉ k z k | k 1 .
Thus, the algorithm flow of dynamic load identification can be obtained.
Second is the case where the mass coefficients are unknown. We define the symmetric matrix M 1 as
M 1 = 1 m 1 1 m 2 1 m n .
All n mass coefficients of the structure are then inserted in the extended state vector. We can thus get the partial derivatives of the state updating equation and the measurement updating equation as
z f k 1 c = f c p f c p ˙ f c α z = z k k = 0 I 0 M 1 K M 1 C E k 0 0 0 ,
z h k = h p h p ˙ h α z = z k k 1 = H 0 M 1 K H 0 M 1 C H 0 E ˉ k .
The symmetric matrices E k and E ˉ k can be written as
E k = 1 m 1 2 ξ k 1 1 m 2 2 ξ k 2 1 m n 2 ξ k n , E ˉ k = 1 m 1 2 ξ ˉ k 1 1 m 2 2 ξ ˉ k 2 1 m n 2 ξ ˉ k n .
Let ξ k ( i ) , ξ ˉ k ( i ) be the i-th element of the vector ξ k , ξ ˉ k at time ( k 1 ) Δ t
ξ k = K p k | k C p k | k + B u f k ,
ξ ˉ k = K p k | k 1 C p k | k 1 + B u f k 1 .
The partial derivatives of the state updating equation and the measurement updating equation with respect to the load vector are the same as the Equations (35) and (36). However, they are not constant matrices in this case because the mass coefficient will change in the filtering process. The mass coefficient must be replaced back to the mass matrix in each filtering. The algorithm procedure is then the same as that of the unknown parameter stiffness or damping coefficient, but the difference is in calculating the Jacobian matrix.
For simplicity, we summarize the algorithm steps of the load identification in Algorithm 1.
Algorithm 1. Flow of load identification based on extended Kalman filter
  1. Given initial conditions
               z 0 | 1 , P 0 | 1 z a n d   f 0
  2. Excitation identification step
            R ˜ k = ( z h k ) T · P k | k 1 z · z h k + R
           J k = [ ( f h k ) T R ˜ k 1 f h k ] 1 ( f h k ) T R ˜ k 1
          f ^ k = J k y k h k z k | k 1 , f ^ k 1 + f h k · f ^ k 1
             P k f = [ ( f h k ) T · R ˜ k 1 · f h k ] 1
  3. Measurement update step
             K k = P k | k 1 z ( z h k ) T R ˜ k 1
      z k | k = z k | k 1 + K k y k h z k | k 1 , f ^ k 1 f h k · f ^ k f ^ k 1
        P k | k z = P k | k 1 z K k · ( R ˜ k f h k · P k f · f h k ) T · K k
            P k | k z f = ( P k | k f z ) T = K k · f h k · P k f
  4. State update step
              z k + 1 | k = z k | k + t k t k + 1 f c z k | k , f ^ k d t
      P k + 1 | k z = 1 + Δ t · z f k c Δ t · f f k c P k | k z P k z f P k f z P k f ( 1 + Δ t · z f k c ) T ( Δ t · f f k c ) T + Q

3. Model Reduction Strategy

The proposed algorithm works online to recover the applied load and the system characteristics as they change in real-time, However, when considering a large system with thousands of degrees of freedom, it becomes impossible to perform the computations in real-time. In this case, the model reduction technique can be used to reduce the involved computations. Here, we adapt the reduction strategy. The main idea in this strategy is to divide the degrees of freedom of the finite element model into master and slave degrees of freedom. The master degrees of freedom are generally those that have forces applied at them, or we need to read the results at them. The other degrees of freedom can be defined as slave degrees of freedom. The slave degrees of freedom are expressed in terms of the master degree of freedom which significantly reduces the required computations.
Next, we explain the reduction strategy using a beam finite element with four degrees of freedom associated to two nodes so that each node has a rotational and a translational degree of freedom. The rotational degrees of freedom are difficult to measure, while the translational ones are easy. Using the reduction method, the dynamic equations of the beam elements can be transformed into condensed dynamic equations with only translational degrees of freedom.
The mass associated to the translational degrees of freedom is considered and is assumed to be concentrated, while the rotational degrees of freedom mass are ignored. In this case, the stiffness matrix and the mass matrix of the beam element can be rewritten as
       u 1        u 2        θ 1            θ 2    u 1    u 2    θ 1    θ 2 K i = E i I i l 3 · 12 12 6 l 6 l 12 12 6 l 6 l 6 l 6 l 4 l 2 2 l 2 6 l 6 l 2 l 2 4 l 2 ,   M i = ρ i A i l 2 · 1 0 0 0 0 1 0 0 0 0 0 0 0 0 0 0 .
The global stiffness and mass matrices are composed of their counterpart elementary matrices. The first half of the global stiffness and mass matrices correspond to the translational degrees of freedom, while the second half corresponds to the rotational degrees of freedom. Then, the global stiffness matrix and mass matrix can be written as
K = K u u K u θ K θ u K θ θ ,   M = M u u 0 0 0 .
Substituting the global stiffness and mass matrices into the differential equation of motion, we can get
M u u 0 0 0 X ¨ u X ¨ θ + α M u u 0 0 0 + β K u u K u θ K θ u K θ θ X ˙ u X ˙ θ + K u u K u θ K θ u K θ θ X u X θ = F t 0 .
Among them, the parameter X u and X θ mean translational displacement and rotational displacement, respectively. We also assume that the external load only acts on the beam structure in the translational direction. Thereafter, Equation (47) can be divided into two equations:
M u u X ¨ u + α M u u + β K u u X ˙ u + β K u θ X ˙ θ + K u u X u + K u θ X θ = F t ,
β K θ u X ˙ u + β K θ θ X ˙ θ + K θ u X u + K θ θ X θ = 0 .
We can get
β X ˙ θ + X θ = K θ θ 1 K θ u β X ˙ u + X u .
Substituting Equation (50) into Equation (48), we can get the motion equation of reduction
M g u · X ¨ u + C g u · X ˙ u + K g u · X u = F t .
where
M g u = M u u ,
C g u = α M u u + β K u u K u θ K θ θ 1 K θ u ,
K g u = K u u K u θ K θ θ 1 K θ u .
Substitute the condensed dynamic equations into the proposed dynamic load identification algorithm. It can be seen that the process of load identification is the same as that without the reduction. Similar to before and after the reduction, the state update equation and the measurement update equation in the algorithm can be rewritten as
d d t X u X ˙ u α = X ˙ u M g u 1 K g u X u C g u X ˙ u + B u f t 0 = f c z t , f t ,
y t = H 0 M g u 1 C g u X ˙ u K g u X u + B u f t = h z t , f t .

4. Numerical Validation

In this section, three numerical examples are used to validate the proposed algorithm. The first two represent systems of multi-degrees of freedom while the third a beam. Two of the examples are dedicated to systems with unknown mass while the third to unknown stiffness. The measurement taken in the examples is assumed to be contaminated with low and high levels of noise. In the last example, the proposed algorithm is combined with the reduction algorithm to identify the system and its load.

4.1. Three-Degrees-of-Freedom with Varying Mass

In the first example, a three-degrees-of-freedom system is selected. A schematic diagram of the structure is shown in Figure 1. The damping of the system is assumed to be Rayleigh damping, with the damping coefficients α = 0.05, β = 0.02.
Here, m1 = m2 = m3 = 1 kg, k1 = k2 = k3 = k4 = 200 N/m. We assume the mass coefficient m2 increases gradually from 1 kg to 3 kg starting 1.5 s and finishing at 3.5 s. The responses of the three masses are selected as the measurement vector, while the three masses are included in the extended state vector as unknown parameters. Assuming that the estimation of the initial response is correct and the estimation of the initial parameters is biased, the initial extended state vector can be set as z 0 | 1 = [ 0 , 0 , 0 , 0 , 0 , 0 , 1 , 3 , 4 ] T , which is a 9 × 1 column vector. We set the covariance matrix of the model noise to Q = d i a g 1 × 10 6 ones 1 , 6 , 1 × 10 4 ones 1 , 3 , and the covariance matrix of the measurement noise to R = 1 × 10 6 eye 3 .
A load of f = s i n ( 5 π t ) + 2 s i n ( 2 π t ) is applied on the mass block m1 and 5%, 10% noise are applied to the measurement vector, respectively. The identification results are shown in the Figure 2, Figure 3, Figure 4 and Figure 5. The left plots are for the case with 5% noise, while the right side is for the case with 10% noise.
The results show that the proposed algorithm can accurately and simultaneously identify the varying masses and the unknown loads on the structure when the noise level is 5% and 10%. The values of the unknown parameters identified by the algorithm are recorded before the mass changes (t = 1 s) and after it changes (t = 5 s), and are compared with the real values in Table 1. The identification errors are shown in the same table. It can be seen that the algorithm has successfully identified the unknown mass values.
In order to illustrate the identification accuracy of the unknown loads, the relative error (RE) and the correlation coefficient (r) methods are used to verify the results. The relative error refers to the ratio of the absolute error between the identification result and the exact value divided by the exact value. The smaller the relative error is, the higher the accuracy of the identification result. The correlation coefficient (r) is the quantity to study the level of the linear correlation between the variables and is expressed as a percentage. The closer the correlation coefficient is to 1, the better the identification. The errors are shown in Table 2. It can be seen that the algorithm can accurately identify the unknown load in general even when the measured data is contaminated with 5% or 10% noise.

4.2. Five-Degree-of-Freedom with Varying Stiffness

In the second numerical test, we study a five-degrees-of-freedom system. We show in Figure 6 a schematic diagram of the structure. In this test, we consider a gradually varying stiffness. The damping of the system is assumed to be Rayleigh damping of the coefficients α = 0.05, β = 0.02.
The considered masses and stiffnesses are m1 = m2 = m3 = m4 = m5 = m6 = 1 kg and k1 = k2 = k3 = k4 = k5 = k6 = 200 N/m. The stiffness coefficient k4 of the structure decreases gradually from 200 N/m to 120 N/m starting at the time instant t = 1.5 s and finishes at t = 3.5 s. The response of all the mass blocks is recorded in the measurement vector. Assuming that the stiffness coefficients k1 and k2 are accurately estimated, the remaining stiffness coefficients are substituted into the extended state vector as unknown parameters. We start by considering the initial response estimation to be correct and that the initial parameter estimation has a certain deviation. The initial extended state vector can be set to z 0 | 1 = ( z e r o s 1 , 10 , 120 , 220 , 160 , 180 ) T , which is a 14 × 1 column vector.
The mass m1 is excited with a load given by f 1 = s i n ( 5 π t ) + 2 s i n ( 2 π t ) , while the mass m2 by a load f 2 = s i n ( 6 π t ) . We consider the measurement to be polluted with noise so that we first add 1% and then 5% white noise to the measurement vector. Again, we use the proposed algorithm to track the unknown parameters as well as identify the unknown load. Figure 7, Figure 8, Figure 9, Figure 10 and Figure 11 show the identified stiffness while Figure 12 and Figure 13 show the identified loads. The left plots show the identification results when the measurement is contaminated with 1% noise, while the right plots with 5% noise. The values of the unknown parameters identified by the algorithm are recorded before the parameters change at t = 1 s and again after the parameters change at t = 5 s, and the results are compared to the real values in Table 3. The corresponding errors are also shown in the table. The errors in the identified load under different noise conditions are shown in the Table 4.
The results show that the proposed algorithm can successfully identify the unknown load with varying stiffness even when the measurement is polluted with noise. The algorithm is also capable of tracking the changes in the stiffness as it recovers the load. The results suggest that the identification results of the proposed algorithm are accurate.

4.3. Cantilever with Varying Mass

In the third example, we study a cantilever with a varying mass. The considered beam length is L = 1 m with a rectangular cross section of the width b = 0.1 m and the height h = 0.01 m. The material density is ρ = 2770   kg / m 3 . The Young’s modulus is E = 71 GPa. The damping characteristic of the cantilever is assumed to be Rayleigh damping with the coefficients being α = 0.8 , β = 0.004 . We mesh the cantilever using six beam finite elements. Figure 13 shows the cantilever and the considered mesh. In this example, we consider the mass ρ i A i of the beam element to be gradually changing. We assume the mass characteristics ρ 4 A 4 of the element number 4 to decrease linearly between 1.5 s and 3.5 s to 60% of its original value. The masses of the other elements do not change. A sinusoidal force F = 200 s i n ( 5 π t ) is applied at the end of the cantilever. Again, two cases are considered where 1% and then 5% white Gaussian noise is added to the measured response.
The reduction method proposed in Section 3 is used to reduce the dynamic equation of the beam structure to a dynamic equation with only translational degrees of freedom. Then, the proposed algorithm is used to identify the load and the parameters of the condensed beam structure. We assume that the mass coefficient ρ 6 A 6 is accurately estimated, while other elements’ mass characteristics are assumed unknown. Thus, the extended state vector of the proposed algorithm can be written as z = ( p , p ˙ , ρ 1 A 1 , ρ 2 A 2 , ρ 3 A 3 , ρ 4 A 4 , ρ 5 A 5 ) T . Let the initial extended state vector be z 0 | 1 = ( z e r o s 1 , 12 , ρ A , 0.6 ρ A . ones 1 , 4 ) T , which is a 17 × 1 column vector, and all the translational degrees of freedom on the beam are observed. Figure 14, Figure 15, Figure 16, Figure 17 and Figure 18 show the identification results of the unknown mass parameters, while Figure 19 shows the identified load. Again, on the left side, we show the results when 1% noise is applied to the measurement, while to the right are the results correspond to 5% noise pollution.
The results’ accuracy clearly shows that the reduction algorithm works well with the proposed load identification algorithm. The load as well as the unknown mass of each beam element are recovered with good accuracy. The errors of the identification results are shown in Table 5 and Table 6 for the mass and the load, respectively. The results with both levels of noise are accurate with engineering accuracy, i.e., the error is the same level as that of the pollution noise. Obviously, increasing the noise level results in larger errors, but the algorithm remains to provide meaningful results.

5. Experimental Validation

As a final validation, an online dynamic load identification experiment of a simply supported beam with varying mass characteristics is designed to verify the proposed algorithm experimentally. The structural parameters of the simply supported beam are shown in Table 7.
Before we run the experiment, the natural frequencies of the simply supported beam are evaluated both experimentally and numerically using the finite element method. The natural frequencies and natural mode shapes of the beam are shown in Table 8 where the numerical and experimental results are compared. It can be seen that the natural frequencies of the experimental results are very close to that of the simulation model. This is reflected in the first four mode shapes where the relative errors are less than 1%. Thus, we conclude that the finite element model is reliable.
To vary the beam mass, we fix an aluminum box at 0.48 m from the left of the simply supported beam. The aluminum box is filled with steel sand. To change the beam mass, the box has holes at the bottom so that it leaks the steel sand as the experiment progresses. The aluminum box with the iron sand and the holes is shown in Figure 20, while a schematic diagram of the experiment is shown in Figure 21. The experimental setting is shown in Figure 22. All the used instruments are listed in Table 9. The sampling rate is set to 8196 Hz, and the excitation signal on the exciter is measured by a force sensor. The beam is meshed into eight beam elements. The aluminum box is fixed to element 6. The acceleration responses are measured by sensors placed at the nodes as shown in Figure 21. Two loading cases are considered. First, a sinusoidal load is applied, and then the experiment is repeated for a triangular wave load.
Using the acceleration measurements at the nodes, the proposed identification algorithm is used to identify the unknown load and the unknown changing mass. The identification results are evaluated by comparing them to the actual load measured by the force sensor. The mass identification results and evaluation results are shown in Figure 23, Figure 24, Figure 25 and Figure 26. Table 10 shows the error in the identified mass of individual elements. The errors are shown before the box has started leaking the sand and then again immediately after all the sand has leaked.
It can be seen that the proposed algorithm can identify the unknown mass with good accuracy regarding the parameters of the structure as well as the unknown loads acting on it. Other than the identification error of the mass characteristic ρ 1 A 1 , the identification of other mass coefficients is good. The errors between the identification results and the real values are due to the noise in the measurement responses and the inaccuracy caused by the approximation of the reduction method. Since most of the identification errors of these mass characteristic is within a close range of the element mass, it can be considered that the parameters identified by the proposed algorithm are accurate. Figure 27 and Figure 28 show the identification results of the sinusoidal loads before and during the beam mass change. While Figure 29 and Figure 30 show the identification results of the triangular loads before and during the beam mass change. In general, it can be seen that the algorithm has captured the loads with good accuracy. To quantify the identification accuracy further, we show in Table 11 the relative errors and the correlation coefficient for the identified results against the experimental results.
It can be seen in the figures that the accuracy of the load identification is indeed accurate where all the details of the applied load are captured. However, in the initial stage for the first 0.05 s, there are fluctuations in the results. This is caused by the initial errors in estimating the mass coefficients at this initial phase. With the identification progress, the mass parameter converges to accurate values, and the load identification accuracy significantly improves. The errors in Table 11 show the correlation coefficients of the sinusoidal load and the triangular load to be more than 99%, which suggests an excellent overall accuracy in the algorithm. This accuracy in identifying the dynamic load and in capturing the mass parameters has been achieved despite the varying nature of the problem.

6. Conclusions

In this paper, an online dynamic load identification method based on an extended Kalman filter is proposed. The proposed algorithm is also combined with the reduction method which is used to reduce the number of degrees of freedom in the system. The reduced system leads to improved efficiency by reducing the rank of the matrix to be solved, hence, the amount of necessary computation. The algorithm is implemented for cases when the unknown parameters are either the stiffness or the mass of a structure. It can be seen that the different unknown parameter types only affect the Jacobian matrix, but the algorithm flow in a similar way for both cases.
Several examples are used to test and validate the proposed approach. The studied examples include multi-degrees of freedom systems as well as a cantilever where the measured data are contaminated with different levels of noise. The algorithm was also validated in a lab experiment designed to validate the approach in a physical test on a simply supported beam with a varying mass. The reduction method is successfully implemented with the proposed algorithm to eliminate the rotational degrees of freedom in the beam elements. The numerical as well as the experimental results prove that the proposed algorithm can simultaneously identify the change in the structural parameters and capture them as well as identify the unknown dynamic load acting on the structure. All the presented results show the accuracy of the proposed algorithm and its efficiency in capturing the dynamic loads as well as the changing parameters in the system even when the measurement data are polluted with noise. Finally, for future works, we aim to develop a theoretical proof for the method converge. However, developing such a theoretical proof can be especially difficult given the varying nature of the identified parameters.

Author Contributions

Conceptualization, H.L. and J.J.; methodology, H.L., J.J. and M.S.M.; software, H.L.; validation, J.J., and M.S.M.; data curation, H.L., J.J. and M.S.M.; writing—original draft preparation, H.L.; writing—review and editing, J.J. and M.S.M.; funding acquisition, H.L., J.J. and M.S.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China, No. 51775270, Qing Lan Project, the Research Fund Incubation Project provided by Jinling Institute of Technology, No. jit-fhxm-201914. The financial support provided by the project of Qatar National Research Fund under the contract NPRP11S-1220-170112.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

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

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Bernal, D.; Gunes, B. Flexibility Based Approach for Damage Characterization: Benchmark Application. J. Eng. Mech. 2004, 130, 61–70. [Google Scholar] [CrossRef]
  2. Cheng, W.; Bo, Y. A Time Domain Identification Method for Distributed Dynamic Loads in Spatial Structures. Appl. Math. Mech. 2017, 38, 967–978. [Google Scholar]
  3. Koh, C.G.; See, L.M.; Balendra, T. Estimation of structural Parameters in time domain: A substructure approach. Earthq. Eng. Struct. Dyn. 1991, 20, 787–801. [Google Scholar] [CrossRef]
  4. Brincker, R.; Zhang, L.; Andersen, P. Modal identification of output-only systems using frequency domain decomposition. Smart Mater. Struct. 2001, 10, 441–445. [Google Scholar] [CrossRef] [Green Version]
  5. Weixin, R.; Harik, I.E. Modal analysis of the Cumberland river bridge on I-24 highway in West Kentucky. In Proceedings of the 19th International Modal Analysis Conference, Kissimmee, FL, USA, 5–8 February 2001; pp. 21–27. [Google Scholar]
  6. Huokun, L.; Jijian, L. Operational Modal Parameter Identification of High Arch Dam under Flood Discharge Excitation Based on Frequency Domain Decomposition Method. J. Vib. Shock 2008, 27, 149–153. [Google Scholar]
  7. Xu, B.; He, J.; Rovekamp, R.; Dyked, S.J. Structural parameters and dynamic loading identification from incomplete measurements: Approach and validation. Mech. Syst. Signal Process. 2012, 28, 244–257. [Google Scholar] [CrossRef]
  8. Cole, H.A. On Line Failure Detection and Damping Measurement of Aerospace Structure by Random Decrement Signature. AIAA 1968, 68, 288–319. [Google Scholar]
  9. Ibrahim, S.R. Efficient Random Decrement Computation for identification of Ambient Responses. In Proceedings of the 19th International Modal Analysis Conference, Kissimmee, FL, USA, 5–8 February 2001; pp. 1–6. [Google Scholar]
  10. Peeters, B.; De Roeck, G. Reference-based stochastic subspace identification for out-put-only modal analysis. Mech. Syst. Signal Process. 1999, 13, 855–878. [Google Scholar] [CrossRef] [Green Version]
  11. Juang, J.N.; Pappa, R.S. An eigensystem realization algorithm for modal parameter identification and model reduction. J. Guid. Control Dyn. 1985, 8, 620–627. [Google Scholar] [CrossRef]
  12. James, G.H.; Carne, T.G.; Lauffer, J.P. The natural excitation technique(NEXT)for modal parameter extraction from operating structures. Int. J. Anal. Exp. Modal Anal. 1995, 10, 260–277. [Google Scholar]
  13. Arrar, C.R.; James, G.H. System identification from ambient vibration measurements on a Bridge. J. Sound Vib. 1997, 205, 1–18. [Google Scholar]
  14. Hoshiya, M.; Saito, E. Structural Identification by Extended Kalman Filter. J. Eng. Mech. 1984, 110, 1757–1770. [Google Scholar] [CrossRef]
  15. Loh, C.H.; Lin, C.Y.; Huang, C.C. Time Domain Identification of Frames under Earthquake Loading. J. Eng. Mech. 2000, 126, 693–703. [Google Scholar] [CrossRef]
  16. Sato, T.; Qi, K. Adaptive H∞ Filter: Its Application to Structural Identification. J. Eng. Mech. 1998, 124, 1233–1240. [Google Scholar] [CrossRef]
  17. Yoshida, I. Damage detection using Monte Carlo filter based on non-Gaussian noise. In Proceedings of the Structural Safety and Reliability, ICOSSAR’01, Newport Beach, CA, USA, 17–22 June 2001; Swet Zeitinger: Lisse, The Netherlands, 2001. [Google Scholar]
  18. Ellis, N. A technique for evaluating the fluctuating aerodynamic forces on a flexible building. In Proceedings of the International Symposium on Vibrating Problems in Industry, Keswick, UK, 10–12 April 1973. [Google Scholar]
  19. Bartlett, F.D.; Flannelly, W.D. Model verification of force determination for measuring vibration loads. J. Am. Helicopter Soc. 1979, 19, 10–18. [Google Scholar] [CrossRef]
  20. Hirotsugu, I.; Hiroyuki, I.; Kikuo, K.; Toshikazu, S. Measurement of Impact Load by Using an Inverse Analysis Technique. JSME Int. J. Ser. 1 Solid Mech. Strength Mater. 1991, 34, 453–458. [Google Scholar]
  21. Milana, S.; Sinibaldi, G.; Marino, L.; Culla, A. Load identification by coherence analysis of structural response. Acoust. Soc. Am. J. 2017, 141, 3833. [Google Scholar] [CrossRef]
  22. Yang, J.N.; Huang, H.; Lin, S. Sequential non-linear least-square estimation for damage identification of structures. Int. J. Non-Linear Mech. 2006, 41, 124–140. [Google Scholar] [CrossRef]
  23. Zhang, F.; Tang, X.; Qin, Y. The Identification Method Research of Dynamic Random Load Distributing on Structure. J. Vibrat. Shock 2012, 31, 1–4, 15. [Google Scholar]
  24. Jiang, J.; Mohamed, M.S.; Seaid, M.; Li, H. Identifying the wavenumber for the inverse Helmholtz problem using an enriched finite element formulation. Comput. Methods Appl. Mech. Eng. 2018, 340, 615–629. [Google Scholar] [CrossRef] [Green Version]
  25. Jiang, J.; Mohamed, M.S.; Seaid, M.; Li, H. Inverse algorithm for real-time road roughness estimation for autonomous vehicles. Arch. Appl. Mech. 2020, 90, 1333–1348. [Google Scholar] [CrossRef]
  26. Ory, H.; Glaser, H.; Holzdeppe, D. Reconstruction of forcing functions based on measured structural responses. In Proceedings of the 2nd International Symposium on Aeroelasticity and Structural Dynamics, Aachen, Germany, 1–3 April 1985; pp. 656–668. [Google Scholar]
  27. Wu, E.; Tsai, T.D.; Yen, C.S. Two methods for determining impact-force history on elastic plates. Exp. Mech. 1995, 35, 11–18. [Google Scholar] [CrossRef]
  28. Gunawan, F.E.; Homma, H.; Morisawa, Y. Impact force estimation by quadratic spline approximation. J. Solid Mech. Mater. Eng. 2008, 2, 1092–1103. [Google Scholar] [CrossRef]
  29. Xuegang, S.; Peng, L.; Zhuming, C. An Algorithm of Dynamic Load Identification Based on FGB Sensor and Kalman Filter. Acta Opt. Sin. 2018, 38, 0328012. [Google Scholar] [CrossRef]
  30. Kalman, R.E. A New Approach to Linear Filtering and Prediction Problems. J. Fluids Eng. 1960, 82, 35–45. [Google Scholar] [CrossRef] [Green Version]
  31. Yang, J.N.; Pan, S.; Huang, H. An adaptive extended Kalman filter for structural damage identifications II: Unknown inputs. Struct. Control Health Monit. 2007, 14, 497–521. [Google Scholar] [CrossRef]
  32. Lenkutis, T.; Cerškus, A.; Šešok, N. Road Surface Profile Synthesis: Assessment of Suitability for Simulation. Symmetry 2021, 13, 68. [Google Scholar] [CrossRef]
  33. Lei, Y.; Jiang, Y.; Xu, Z. Structural damage detection with limited input and output measurement signals. Mech. Syst. Signal Process. 2012, 28, 229–243. [Google Scholar] [CrossRef]
  34. Naets, F.; Croes, J.; Desmet, W. An online coupled state/input/parameter estimation approach for structural dynamics. Comput. Methods Appl. Mech. Eng. 2015, 283, 1167–1188. [Google Scholar] [CrossRef]
  35. Du, H.; Kong, W.; Wang, Y.; Liu, W.; Huang, M.; Zhang, W.; Tang, M. Research on Rotorcraft Blade Tip Vortex Identification and Motion Characteristics in Hovering State. Symmetry 2020, 12, 196. [Google Scholar] [CrossRef] [Green Version]
  36. Huang, Q.; Xu, Y.L.; Liu, H.J. An efficient algorithm for simultaneous identification of time-varying structural parameters and unknown excitations of a building structure. Eng. Struct. 2015, 98, 29–37. [Google Scholar] [CrossRef]
  37. Lei, Y.; Xia, D.; Erazo, K.; Nagarajaiah, S. A novel unscented Kalman filter for recursive state-input-system identification of nonlinear systems. Mech. Syst. Signal Process. 2019, 127, 120–135. [Google Scholar] [CrossRef]
  38. Xiaoxiong, Z.; Jia, H. Identification of Structural Parameters and Unknown Excitations Basend on the Extended Kalman Filter. Eng. Mech. 2019, 36, 221–230. [Google Scholar]
  39. Anhua, Y. Parameter Identification of Shear-Type Frame Structures Based on Kalman Filter; Fujian Agriculture and Forestry University: Fuzhou, China, 2019. [Google Scholar]
Figure 1. Schematic diagram of the three-degrees-of-freedom system.
Figure 1. Schematic diagram of the three-degrees-of-freedom system.
Symmetry 13 01372 g001
Figure 2. Identification results of the mass coefficient m1.
Figure 2. Identification results of the mass coefficient m1.
Symmetry 13 01372 g002
Figure 3. Identification result of the mass coefficient m2.
Figure 3. Identification result of the mass coefficient m2.
Symmetry 13 01372 g003
Figure 4. Identification result of the mass coefficient m3.
Figure 4. Identification result of the mass coefficient m3.
Symmetry 13 01372 g004
Figure 5. Identification result of the unknown load.
Figure 5. Identification result of the unknown load.
Symmetry 13 01372 g005
Figure 6. Schematic diagram of the five-degree-of-freedom system.
Figure 6. Schematic diagram of the five-degree-of-freedom system.
Symmetry 13 01372 g006
Figure 7. Identification results of the unknown stiffness k3.
Figure 7. Identification results of the unknown stiffness k3.
Symmetry 13 01372 g007
Figure 8. Identification results of the unknown stiffness k4.
Figure 8. Identification results of the unknown stiffness k4.
Symmetry 13 01372 g008
Figure 9. Identification results of the unknown stiffness k5.
Figure 9. Identification results of the unknown stiffness k5.
Symmetry 13 01372 g009
Figure 10. Identification results of the unknown stiffness k6.
Figure 10. Identification results of the unknown stiffness k6.
Symmetry 13 01372 g010
Figure 11. Identification result of the unknown load f1.
Figure 11. Identification result of the unknown load f1.
Symmetry 13 01372 g011
Figure 12. Identification result of the unknown load f2.
Figure 12. Identification result of the unknown load f2.
Symmetry 13 01372 g012
Figure 13. Schematic diagram of cantilever structure.
Figure 13. Schematic diagram of cantilever structure.
Symmetry 13 01372 g013
Figure 14. Identification results of the mass ρ 1 A 1 .
Figure 14. Identification results of the mass ρ 1 A 1 .
Symmetry 13 01372 g014
Figure 15. Identification results of the mass ρ 2 A 2 .
Figure 15. Identification results of the mass ρ 2 A 2 .
Symmetry 13 01372 g015
Figure 16. Identification results of the mass ρ 3 A 3 .
Figure 16. Identification results of the mass ρ 3 A 3 .
Symmetry 13 01372 g016
Figure 17. Identification results of the mass ρ 4 A 4 .
Figure 17. Identification results of the mass ρ 4 A 4 .
Symmetry 13 01372 g017
Figure 18. Identification results of the mass ρ 5 A 5 .
Figure 18. Identification results of the mass ρ 5 A 5 .
Symmetry 13 01372 g018
Figure 19. Identification result of the unknown load.
Figure 19. Identification result of the unknown load.
Symmetry 13 01372 g019
Figure 20. The aluminum box with iron sand.
Figure 20. The aluminum box with iron sand.
Symmetry 13 01372 g020
Figure 21. The test scheme of the dynamic load identification test.
Figure 21. The test scheme of the dynamic load identification test.
Symmetry 13 01372 g021
Figure 22. The test site of the dynamic load identification test.
Figure 22. The test site of the dynamic load identification test.
Symmetry 13 01372 g022
Figure 23. Mass identification results ( ρ 1 A 1 on the left and ρ 2 A 2 on the right).
Figure 23. Mass identification results ( ρ 1 A 1 on the left and ρ 2 A 2 on the right).
Symmetry 13 01372 g023
Figure 24. Mass identification results ( ρ 3 A 3 on the left and ρ 4 A 4 on the right).
Figure 24. Mass identification results ( ρ 3 A 3 on the left and ρ 4 A 4 on the right).
Symmetry 13 01372 g024
Figure 25. Mass identification results ( ρ 5 A 5 on the left and ρ 6 A 6 on the right).
Figure 25. Mass identification results ( ρ 5 A 5 on the left and ρ 6 A 6 on the right).
Symmetry 13 01372 g025
Figure 26. Mass identification results ( ρ 7 A 7 on the left and ρ 8 A 8 on the right).
Figure 26. Mass identification results ( ρ 7 A 7 on the left and ρ 8 A 8 on the right).
Symmetry 13 01372 g026
Figure 27. Identification result of the sinusoidal wave load in the initial phase before the mass change.
Figure 27. Identification result of the sinusoidal wave load in the initial phase before the mass change.
Symmetry 13 01372 g027
Figure 28. Identification result of the sinusoidal wave load during the mass change phase.
Figure 28. Identification result of the sinusoidal wave load during the mass change phase.
Symmetry 13 01372 g028
Figure 29. Identification result of the triangular load in the initial phase before the mass change.
Figure 29. Identification result of the triangular load in the initial phase before the mass change.
Symmetry 13 01372 g029
Figure 30. Identification result of the triangular load during the mass change phase.
Figure 30. Identification result of the triangular load during the mass change phase.
Symmetry 13 01372 g030
Table 1. Parameter identification error.
Table 1. Parameter identification error.
Mass Coefficientt = 1 st = 5 s
Identification
(kg)
Exact
(kg)
Error
(%)
Identification
(kg)
Exact
(kg)
Error
(%)
5%
noise
m11.01311.30.986711.3
m21.08718.73.02230.7
m30.985611.41.00610.6
10%
noise
m11.00210.21.142114.2
m21.02512.53.04934.9
m31.05115.11.139113.9
Table 2. Load identification error.
Table 2. Load identification error.
Noise LevelRelative Error (RE) (%)Correlation Coefficient (r) (%)
5% noise9.5599.48
10% noise17.7598.44
Table 3. Parameter identification error.
Table 3. Parameter identification error.
Stiffness Coefficientt = 1 st = 5 s
Identification
(N/m)
Exact
(N/m)
Error
(%)
Identification
(N/m)
Exact
(N/m)
Error
(%)
1%
noise
k3198.82000.6199.32000.35
k4199.22000.4118.41201.33
k5199.62000.2198.92000.55
k6199.82000.1199.22000.4
5%
noise
k3187.22006.4194.42002.8
k4191.82004.1118.81201
k5198.72000.65197.52001.25
k62012000.5195.62002.2
Table 4. Load identification error.
Table 4. Load identification error.
Noise LevelRelative Error (RE) (%)Correlation Coefficient (r) (%)
1% noisef12.6499.97
f24.9899.88
5% noisef110.2499.32
f219.8598.05
Table 5. Parameter identification error.
Table 5. Parameter identification error.
Mass Characteristict = 1 st = 5 s
Identification
(kg/m)
Exact
(kg/m)
Error
(%)
Identification
(kg/m)
Exact
(kg/m)
Error
(%)
1%
noise
ρ 1 A 1 2.762.773.612.9032.774.8
ρ 2 A 2 2.8182.771.732.5062.779.53
ρ 3 A 3 2.7652.770.182.9442.776.28
ρ 4 A 4 2.7732.770.111.5541.6626.49
ρ 5 A 5 2.7852.770.542.9082.774.98
5%
noise
ρ 1 A 1 2.7732.770.112.8292.772.13
ρ 2 A 2 2.2852.7717.512.3492.7715.20
ρ 3 A 3 3.1722.7714.513.1552.7713.89
ρ 4 A 4 2.5182.779.091.8631.6627.26
ρ 5 A 5 2.9242.775.562.9492.776.46
Table 6. Load identification error.
Table 6. Load identification error.
Noise LevelRelative Error (RE) (%)Correlation Coefficient (r) (%)
1% noise2.7499.96
5% noise13.6699.08
Table 7. Parameter table of simply supported beam.
Table 7. Parameter table of simply supported beam.
Parameters of Simply Supported BeamValue
Length a 0.7 m
Width b 0.04 m
Thickness h 0.008 m
Density ρ 7800 kg/m3
Elastic Modulus E 209 GPa
Poisson’s ratio μ 0.30
Table 8. Natural frequency comparison.
Table 8. Natural frequency comparison.
OrderExperimental Results (Hz)Simulation Results (Hz)Error (%)
137.86437.80510.156
2151.74151.2640.314
3337.58340.86370.963
4603.5609.24730.943
Table 9. Force identification experiment equipment.
Table 9. Force identification experiment equipment.
Equipment ClassificationName
Vibration ExciterJZT-2 Permanent magnet exciter
Power amplifierHEAS-50 Power amplifier
Dynamic signal acquisition boardNI PXIe-4499 Capture card
Signal acquisition instrumentNI PXI
SensorPCB 356A33 accelerometer
SoftwareNI Signal Express
Table 10. Parameter identification errors.
Table 10. Parameter identification errors.
Mass CharacteristicBefore Mass Changes After Mass Changed
Identification
(kg/m)
Exact
(kg/m)
Error
(%)
Identification
(kg/m)
Exact
(kg/m)
Error
(%)
ρ 1 A 1 2.0972.52817.052.1682.52814.24
ρ 2 A 2 2.3952.5285.262.5402.5280.47
ρ 3 A 3 2.2462.52811.152.4512.5283.04
ρ 4 A 4 2.5572.5281.142.4772.5282.02
ρ 5 A 5 2.5962.5282.692.6222.5283.72
ρ 6 A 6 5.5825.3853.662.7282.5287.91
ρ 7 A 7 2.5882.5282.372.5622.5281.34
ρ 8 A 8 2.4572.5282.812.5032.5280.98
Table 11. Load identification errors.
Table 11. Load identification errors.
TypeRelative Error (RE) (%)Correlation Coefficient (r) (%)
sine9.3199.57
triangular wave12.9799.16
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Li, H.; Jiang, J.; Mohamed, M.S. Online Dynamic Load Identification Based on Extended Kalman Filter for Structures with Varying Parameters. Symmetry 2021, 13, 1372. https://doi.org/10.3390/sym13081372

AMA Style

Li H, Jiang J, Mohamed MS. Online Dynamic Load Identification Based on Extended Kalman Filter for Structures with Varying Parameters. Symmetry. 2021; 13(8):1372. https://doi.org/10.3390/sym13081372

Chicago/Turabian Style

Li, Hongqiu, Jinhui Jiang, and M Shadi Mohamed. 2021. "Online Dynamic Load Identification Based on Extended Kalman Filter for Structures with Varying Parameters" Symmetry 13, no. 8: 1372. https://doi.org/10.3390/sym13081372

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