Advances in Fault Diagnosis and Defect Detection in Mechanical and Civil Engineering
View this Special IssueResearch Article  Open Access
Ting Wang, Zhimin Wan, Weiguang Zheng, Shuilong He, "Structural Dynamic Load and Parameter Identification Based on Dummy Measurements of Displacement", Shock and Vibration, vol. 2020, Article ID 8886714, 9 pages, 2020. https://doi.org/10.1155/2020/8886714
Structural Dynamic Load and Parameter Identification Based on Dummy Measurements of Displacement
Abstract
Structural load and parameter identification are the essential contents in the field of structural dynamics, and some studies pay attention to the coupled recognition of uncertain structure parameters as well as unacquainted loads. Gillijns and De Moor presented the Kalmantype filter, which was used in the work for coupled identification as the unabbreviated form of GDF. However, it has been demonstrated that it is unstable for the extended GDF (EGDF) method, drifting in the identified unacquainted loads as well as displacement, just like most previous identification methods based on the leastsquare algorithm. In order to deal with this unstable issue, this paper applied the dummy measurements of displacements on a position level and modifies the EGDF algorithm using the information integration method about the accelerated measurements with dummy measurements. Numerical example of a truss is used for validating the applicability of the method in the work, and an influence of covariance matrices in dummy displacements is also considered.
1. Introduction
Structural external load is vital to structural optimization design, failure diagnosis, and health monitoring [1, 2]. Since the external load is difficult to measure by equipment instruments, various deterministic methods have been proposed for load identification. However, structural parameters of the system are uncertain [3], such as uncertain stiffness and damping, which makes the load identification not credible. In recent years, the uncertain methods for structural load and parameter identification have been attracted increasing attention in research.
The simultaneous identification of the parameters of set membership and a load was developed. Meanwhile, it is unusual for the approach needing given states in the applications of engineering [4]. Aiming at the problem of state/input/parameter coupling recognition in structure damage recognition [5], they proposed an EKF (extended Kalman filter algorithm), named EKFUI, based on unacquainted input load. In [6], EKFUI was derived by the leastsquares estimation algorithm, indicating that complicated mathematical reasoning was used to obtain the recursive solution. The above simplified EKFUI in [6] was also applied to identify nonlinear structure parameters. Nevertheless, derivations in the simplified method may bring on the confusion of the Bayesian framework when the prior probability density functions (PDFs) at were regarded as followed PDFs [7]. Naets et al. [8] presented another algorithm called the ADEKF (augmented discrete extended Kalman filter) to recognize coupled parameter/input/state. The ADEKF algorithm looks like the developed KF algorithm [9–11], wherein parameters, uncertain loads, and structural states constitute the extensive state vector in high dimension. More importantly, the ADEKF algorithm requires a matched hypothesis about. A similar method has been used for the offshore wind turbine for the simultaneous identification of the system parameters based on stiffness, states, and the hydrodynamic load in [12]. Lately, the approach of DKFUKF is proposed for the coupling recognition of states, unknown loads, and parameters in [13, 14]. A dual KF (Dual KF) is adopted to estimate the load. Besides, the socalled augment states are estimated by the UKF (unscented KF), consisting of the structural states and the unknown parameters. Nevertheless, the research does not mention the convergence analysis, observability, and stability of the DKFUKF method. Additionally, the DKF algorithm requires the prior supposition about the load covariance, and the assumption value strongly influences the Bayesian filters’ identification quality.
We developed an approach in terms of the GDF algorithm towards the coupling recognition topic, which is called as the extended GDF (EGDF) algorithm [15]. The EGDF algorithm is proposed based on the EKF algorithm’s linearization. Undetermined parameters are added to the structural states forming the extensive ones for recognition. New equations of measurement and state transmission become nonlinear, and the firstorder Taylor expansions are used for linearizing these two equations. The EGDF algorithm and the standard GDF algorithm have a uniform structure, containing updates of time and measurement and load recognition. Their primary difference refers to the statespace formulas’ sensitivity matrices. Nevertheless, the conventional GDF method has been demonstrated to be inherently unstable based solely on limited accelerated measurements. In this situation, it results in the false drift of low frequency in recognizing shift and load [16]. In [15], the information integration of acceleration response and the shift is also used for solving the spurious lowfrequency drift problem. However, the displacement transducer has a much larger volume than the acceleration sensor. The introduction of the displacement transducer requires alteration to the structural characteristic of the system, which is unwanted and unpractical. Also, the price of the displacement transducer is higher than that of the acceleration sensor. Acceleration sensors combined with numerical simulation also can be used to detect mechanical faults of gears and bearings [17, 18].
In this paper, only adopting the acceleration responses is taken as the measure signals to modify the EGDF algorithm based on the dummy measurements of displacements. First, the dummy displacement measurements equation is built. Then, the dummy displacement is combined with the measured acceleration responses to filter the unknown load and parameters for obtaining the stable estimated results. Final, the numerical example of a plane truss structure is conducted for validating the algorithm in the work.
2. Reviewing the EGDF Algorithm
The normal formula of movement for the multidegreesoffreedom’s damped linear structure is expressed aswhere , , and represent the structure’s stiffness matrices, damping, and mass, respectively; is the vector of the nodal displacement; and are the structure’s acceleration and velocity responses; indicates the external load vector; is the impact matrix related to the external load . Usually, the structure’s mass can be obtained accurately. Thus, the mass matrix M is given.
In structure system equation (1), it is assumed that and include the unacquainted parameters to be identified; thus, the initial structural status vector is rewritten as the extensive one.in which represents an uncertain argument vector for estimation. refers to the time step size; so the transmission and measurement equations of the system state are changed into the following discretetime forms taking the sampling frequency of as
In these equations, , and and are the two nonlinear functions associated with the augment state vector. The vectors and represent the vectors of measurement and system noises, which are supposed to be zeromean, irrelevant, and white stochastic signal with the given covariation matrices G_{k} and R_{k} separately.
It is assumed that and represent the posterior estimates of and separately, according to the observed vector . The state variance is denoted as , and and are the initial unbiased state estimate vector and the corresponding variance matrix, respectively, which are supposed to be given. Also, the EGDF formula works requiring the two following conditions: (i) the number of measured signals needs to be more significant compared with that of unacquainted loads. (ii) Since the weighted leastsquares estimation is used to calculate the unacquainted input load inversely, the accelerated response signal at the unacquainted input position should be determined.
Based on the above declaration, the EGDF algorithm can be obtained by the following three steps.(i)Input identification step is(ii)Measurement update step is(iii)Time update step iswhere , , and are the sensitivity matrices as
It is found that the EGDF algorithm is actually dependent on the EKF method. Using the firstorder Taylor division to linearize formula (4) of nonlinear state measurement and transmission, the normal GDF formula is appropriate to identify the augmented states and unacquainted offered burden of an approximately linear power system. In addition, it is proved that edgetype algorithms in terms of the limited accelerated measurements are naturally labile and produce the false displacement of low frequency during identifying input loading and shifts [19]. The accelerated insensitivity to quasistatic components of external loads causes this displacement [16]. The author presented the extension of the EGDF algorithm to recognize the extrinsic loads and structure parameters in terms of fusing the date of shift and acceleration [15]. However, displacement transducers have a large volume and are not suitable for placement in some situations, which also may alter the system properties. It is necessary to develop a novel algorithm based solely on acceleration measurements in engineering applications.
3. The EGDF Algorithm Based on Dummy Measurements
In this paper, data fusion of the accelerometers and dummy measurements of displacements is developed based on the EGDF algorithm for structural dynamics load and parameter identification. The dummy measurements approach is motivated by the one proposed by Chatzi and Fuggin [19, 20] to make the monitoring for civil structure stable. Usually, the deformation of a structural system is bounded to a limited range. Besides, an order of magnitude for deformation is used as a priori estimate that may be implemented in simulation. The work regards the boundary of transformation as the uncertainty of virtual measurement and points out that the assumed transformation is zero. Analog measurements of displacement has the following measurement equation:where H_{dm} is the impact matrix related to the dummy displacements. The formula expresses that, under the unacquainted of the covariation R_{dm}, the location variation of the degrees of freedom is 0. The matrix R_{dm} can be used for the tuning of the KF to obtain the desired results, and it ought to choose a higher order of magnitude than the practical movement of the structural systems. The high covariance based on virtual measurements cannot appropriately limit the drifting in recognition. If the covariance is chosen too small, the estimates will be constrained to generate erroneous results.
Adopting the idea of the dummy measurements, the conventional EGDF algorithm can be extended with (18). Equations (5), (10), and (17) are changed into
And also, the matrix D_{k} in (6) is changed aswhere H_{a} represents the influence matrix related to the extrinsic loading, and, from (20), it can been seen that the matrix D_{k} is a constant matrix in this paper.
They do not benefit quick estimation needed for load estimation because of great uncertainty about the virtual measurement. However, they can stop longterm displacement caused through accelerated measurements and stabilize the approximated filter covariance. If the constant location of the given structural system is not zero, this behavior need to be taken into account for the dummy measurements to obtain a more accurate identification.
4. Numerical Simulation
A numeric example of a plane truss is expressed to estimate the modified EGDF formula in the work based on the dummy measurements. The plane truss structure is supported at two ends, modeling through thirtyone finite elements of the plane truss, the lateral and longitudinal degrees of freedom of free nodes (See Figure 1). All horz and vert members have been with the lengths of 2 and , respectively. Besides, the crosssectional area of all elements is , and assuming Rayleigh damping, the related coefficients are and . The specific mass is kgm^{−3}; moreover, the elastic coefficient is Pa. First, inherent frequencies of the structure are 0.15, 0.41, 0.86, 1.02, 1.39, 1.77, 2.14, and 2.29 Hz, respectively, based on the finite element (FE) computation. At nodes 4 and 9, we apply two external loads, respectively. Load is sin excitation.
(a)
(b)
And the other load is an impact input.
Assume the rigidities for truss elements 17, 15, 14, 10, 7, and 5 are unacquainted parameters for identification, with their original values of 633.0, 759.5, 1,163.5, 1,342.5, 633.0, and 759.5 N/m separately. Thus, the extended status vector refers to . Seven accelerated signals, determined for recognition, are the horizontal ones at node 9, as well as the vertical ones at nodes 10, 9, 7, 5, 4, and 3. At the same time, the measurement signals was used for simulation by adding 4% ambient noise. Original augmented status covariance is composed of three parts which are the covariance vector of the displacement, velocity (, ), and unknown stiffness (), and is considered asin which ndof = 30 represents the amount of trussed structures. The covariant matrices of the system noise and the measured noise have been assumed as and separately. In terms of the proposed formula, two external loads have been recognized with recognized curves (see Figures 2 and 3), comparing the active curves of loads. Both identified loads have the drift problem compared with the real values. Structure shift and speed are recognized together. Figure 4 shows the recognized outcome of vertical displacement and velocity at node 7 and theoretic values. Recognized velocity is precise; however, recognized displacement has the drift problem in low frequency due to the sole acceleration measurements.
(a)
(b)
Therefore, local dummy measurement shifts have been introduced to the measurement accelerations. Vertical shifts of nodes 4 and 9 have been taken as dummy measurements together with the above seven accelerations, and R_{dm} = diag [10^{−4}, 10^{−4}]. Figure 5 shows the comparisons of the identified results of load u_{1}, and Figure 6 shows the amplified segment (14–14.5 s) in Figure 5 for better illustration. Figures 7 and 8 show the identified load u_{2}, shift, and speed. Spurious drift reduces obviously using the proposed EGDF based on dummy measurements of displacement. Table 1 shows the identified stiffness values. Recognized relative error (RE) is defined as follows:
(a)
(b)

The recognized parameters are precise.
As in the above description, the covariance matrix R_{dm} influences the estimations much. The two following cases demonstrate the effect of R_{dm} value. One is a smaller value as R_{dm} = diag [10^{−6}, 10^{−6}], and the other is a larger value as R_{dm} = diag [10^{−2}, 10^{−2}]. The identified results are shown in Figures 9–11. It is found in Figures 9 and 10 that the load is not identified accurately, especially the lowfrequency components, and the identified displacement is smaller compared with the accurate value due to the smaller R_{dm}. From Figures 11 and 12, it is found that the identified load is not accurate as the actual load, and the identified displacement has the lowfrequency drift problem due to the larger R_{dm}.
(a)
(b)
(a)
(b)
5. Conclusions
The work proposed the modified formula by fusing the data about determining acceleration and dummy shift reaction for coupled dynamic load and parameter identification. It reduces the false displacement of low frequency in estimating loads and shifts effectively, which can be seen in the conventional EGDF algorithm using acceleration measurements solely. The algorithm is useful depending on the covariance matrix of dummy displacements. The covariance are selected to be an order of magnitude greater compared with the system’s real movement for a suitable estimate. A numerical example is presented for proving the formula in the work, including the influences of different levels of dummy displacement covariance. Additionally, the algorithm is suitable for the modal domain of identification.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that they have no conflicts of interest.
Acknowledgments
The research was funded by the National Science Foundation of China (Grant No. 52065013), the Technology Foundation of Nantong, China (Grant no. JCZ18056), the Natural Science Foundation of the Jiangsu Higher Education Institutions of China (Grant no. 18KJB580014 and 19KJD14003), Natural Science Foundation of Guangxi (Grant no. 2018GXNSFAA281276 and 2018GXNSFBA281012), InnovationDriven Development Special Fund Project of Guangxi (Guike AA18242033), and Liuzhou Science Research and Planning Development Project (Grant nos. 2020GAAA0403, 2019AD10203, and 2018AA20301).
References
 Y. Gao, X. Liu, and J. Xiang, “FEM simulationbased generative adversarial networks to detect bearing faults,” IEEE Transactions on Industrial Informatics, vol. 16, no. 7, pp. 4961–4971, 2020. View at: Publisher Site  Google Scholar
 Y. Wang, Z. Wei, and J. Yang, “Feature trend extraction and adaptive density peaks search for intelligent fault diagnosis of machines,” IEEE Transactions on Industrial Informatics, vol. 15, no. 1, pp. 105–115, 2019. View at: Publisher Site  Google Scholar
 W. Song, J. Xiang, and Y. Zhong, “Mechanical parameters detection in stepped shafts using the FEM simulation based impulse excitation technique,” Smart Structures and Systems, vol. 20, no. 4, pp. 473–481, 2017. View at: Publisher Site  Google Scholar
 I. Kolmanovsky, I. Sivergina, and J. Sun, “Simultaneous input and parameter estimation with input observers and setmembership parameter bounding: theory and an automotive application,” International Journal of Adaptive Control and Signal Processing, vol. 20, no. 5, pp. 225–246, 2006. View at: Publisher Site  Google Scholar
 J. N. Yang, S. Pan, and H. Huang, “An adaptive extended Kalman filter for structural damage identifications II: unknown inputs,” Structural Control and Health Monitoring, vol. 14, no. 3, pp. 497–521, 2007. View at: Publisher Site  Google Scholar
 Y. Lei, Y. Jiang, and Z. Xu, “Structural damage detection with limited input and output measurement signals,” Mechanical Systems and Signal Processing, vol. 28, pp. 229–243, 2012. View at: Publisher Site  Google Scholar
 Y. Lei, Y. Wu, and T. Li, “Identification of nonlinear structural parameters under limited input and output measurements,” International Journal of Nonlinear Mechanics, vol. 47, no. 10, pp. 1141–1146, 2012. View at: Publisher Site  Google Scholar
 F. Naets, J. Croes, and W. Desmet, “An online coupled state/input/parameter estimation approach for structural dynamics,” Computer Methods in Applied Mechanics and Engineering, vol. 283, pp. 1167–1188, 2015. View at: Publisher Site  Google Scholar
 E. Lourens, E. Reynders, G. De Roeck, G. Degrande, and G. Lombaert, “An augmented Kalman filter for force identification in structural dynamics,” Mechanical Systems and Signal Processing, vol. 27, pp. 446–460, 2012. View at: Publisher Site  Google Scholar
 R. Cumbo, T. Tamarozzi, K. Janssens, and W. Desmet, “Kalmanbased load identification and fullfield estimation analysis on industrial test case,” Mechanical Systems and Signal Processing, vol. 117, pp. 771–785, 2019. View at: Publisher Site  Google Scholar
 R. ÁlvarezBriceño and L. P. R. De Oliveira, “Combining strain and acceleration measurements for random force estimation via Kalman filtering on a cantilevered structure,” Journal of Sound and Vibration, vol. 469, p. 115112, 2020. View at: Publisher Site  Google Scholar
 D. J. M. Fallais, S. Voormeeren, and E. Lourens, “Vibrationbased identification of hydrodynamic loads and system parameters for offshore wind turbine support structures,” Energy Procedia, vol. 94, pp. 191–198, 2016. View at: Publisher Site  Google Scholar
 S. Eftekhar Azam, V. Dertimanis, E. Chatzi, and C. Papadimitriou, “Output only schemes for inputstateparameter estimation of linear systems,” in Proceedings of the 1st International Conference on Uncertainty Quantification in Computational Sciences and Engineering (UNCECOMP 2015), pp. 497–510, Ceete Island, Greece, May 2015. View at: Publisher Site  Google Scholar
 V. Dertimanis, E. Chatzi, S. Eftekhar Azam, and C. Papadimitriou, “Outputonly fatigue prediction of uncertain steel structures,” in Proceedings of the 8th European Workshop on Structural Health Monitoring (EWSHM 2016), Bilbao, Spain, July 2016. View at: Google Scholar
 Z. Wan, T. Wang, L. Li, and Z. Xu, “A novel coupled state/input/parameter identification method for linear structural systems,” Shock and Vibration, vol. 2018, Article ID 7691721, 15 pages, 2018. View at: Publisher Site  Google Scholar
 S. E. Azam, E. Chatzi, and C. Papadimitriou, “A dual Kalman filter approach for state estimation via outputonly acceleration measurements,” Mechanical Systems and Signal Processing, vol. 6061, pp. 866–886, 2015. View at: Publisher Site  Google Scholar
 X. Liu, H. Huang, and J. Xiang, “A personalized diagnosis method to detect faults in gears using numerical simulation and extreme learning machine,” KnowledgeBased Systems, vol. 195, p. 105653, 2020. View at: Publisher Site  Google Scholar
 X. Liu, H. Huang, and J. Xiang, “A personalized diagnosis method to detect faults in a bearing based on acceleration sensors and an FEM simulation driving support vector machine,” Sensors, vol. 20, no. 2, p. 420, 2020. View at: Publisher Site  Google Scholar
 F. Naets, J. Cuadrado, and W. Desmet, “Stable force identification in structural dynamics using Kalman filtering and dummymeasurements,” Mechanical Systems and Signal Processing, vol. 5051, pp. 235–248, 2015. View at: Publisher Site  Google Scholar
 E. N. Chatzi and C. Fuggini, “Structural identification of a supertall tower by gps and accelerometer data fusion using a multirate Kalman filter,” in Proceedings of the 2010 IALCCE Conference, Taipei, Taiwan, October 2010. View at: Google Scholar
Copyright
Copyright © 2020 Ting Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.