Bridges are key components of transportation infrastructure, albeit one in eleven bridges in the U.S. were structurally defi-cient [1]. The high cost and time usually required to inspect aging bridges desperately calls for advanced sensing and data analysis techniques.
The use of vibration signals collected from travelling vehicles to monitor structures (Figure 1) has recently become a viable alternative to traditional structural health monitoring approaches [2, 3, 4, 5, 6, 7]. This approach does not require intensive deployment and maintenance. Previous work on vehicle-vibration-based bridge health monitoring (or indirect BHM) mainly falls into two categories: modal analysis and data-driven approaches. Modal analysis focuses on identifying modal parameters of a bridge, such as natural frequencies
[8, 9, 10], mode shapes [11, 10] and damping [12]. Datadriven approaches use signal processing and machine learning techniques to extract informative features for diagnosing damage [7, 13, 14].
Fig. 1. Vehicle-vibration-based BHM. Features extracted from vehicle vibration signals are used to diagnose damages.
However, to make the indirect BHM (IBHM) approaches practical, there are three main challenges to address. First, since the vibration signals are indirect measurements of structure’s vibrations, modal parameters identified using modal analysis, e.g.,[8, 9, 11, 10, 12], are sensitive to the vehicles properties, environmental factors and noise. Second, purely data-driven methods, e.g., [7, 13, 14, 15], can suffer from overfitting to the available data (i.e., a set of bridges with known damage labels) and achieve significantly worse performance when applied to other bridges. Third, labeled data are limited in quantity. Especially for full-scale bridges, it is expensive, time-consuming, and impractical to obtain vehicle vibrations with corresponding damage labels. It is also unrealistic to damage bridges for sourcing damage labels.
We introduce a physics-guided signal-processing algorithm to extract an informative feature, which can estimate and localize damage in a bridge. To handle the first challenge, the extracted feature should be sensitive to damage instead of to uncertainties. In this work, we use the synchrosqueezed wavelet transform (SWT) [16, 17] to represent the vehicle acceleration signal in the time-frequency plane and reconstruct a DS component using the inverse SWT (ISWT) within a frequency band determined by pre-identified system properties. We use the SWT and ISWT because SWT can represent the non-stationary and time-varying vehicle acceleration as a superposition of intrinsic mode functions (IMF)-type components (our desired feature has the same type), and ISWT
To appear in Proc. ICASSP2020, May 04-08, 2020, Barcelona, Spain cIEEE 2020
can reconstruct our non-stationary component within a frequency range without having the mode mixing effect [18]. To addresses the second and the third challenges, a DS & DI feature is obtained by multiplying the reconstructed component by a DI factor. It is obtained from the solution of a vehicle bridge interaction system (VBIS). Because this feature is DI and not extracted by a trainable model with the supervision of damage states, diagnosing damage using this feature does not suffer from overfitting and can across multiple bridges. We verify the DS & DI properties of our feature by visualization and by comparing it to other time- and frequency-domain features for estimating stiffness reductions and locations.
The main contributions of this paper are 1) We cast the IBHM problem as a signal decomposition, thus affording us the tools from the signal processing community; 2) We use the synchrosqueezed wavelet transform to extract a DS & DI feature for IBHM; and 3) We validate the DS & DI properties of it through extensive experiments.
To model the VBIS, we consider a commonly used model, a sprung mass (representing vehicle) traveling with a constant speed on a simply supported beam (representing bridge). This model provides physical foundations of the IBHM problem, which help us to formulate this problem as a signaldecomposition problem.
Vehicle-bridge interaction system. The derivation of the theoretical formulation of the VBIS follows the same assumptions and system geometry as presented in [7], albeit a local damage is considered in this work.
Let x be the coordinate of the beam with the origin at the left support; the density, cross section area, damping coefficient and stiffness of the beam, respectively;
the Dirac delta function; g the gravity constant; v the moving speed of the vehicle; and y(t) the vertical displacement of the vehicle chassis;
the weight, stiffness and damping coefficient of the vehicle, respectively. The vertical acceleration of the vehicle chassis,
is our measurements. The equations of motion for the VBIS are
where is the central location of the damage;
and
are the damage length and the percentage reduction of the stiffness, respectively;
is the stiffness of the undamaged beam segment.
and
are the parameters we want to infer. Note that many previous works [8, 11, 9, 7] consider the stiffness as a constant. In this paper, we solve the VBIS with the stiffness-reduction-type local damage.
Damage-sensitive and domain-invariant features. We now solve Eq. (2), look for a DS & DI feature and formulate the IBHM problem as a signal decomposition problem. By omitting the damping of the beam and the vehicle, we write the acceleration of the moving vehicle as
where , and
are constants depending on properties of the bridge and the vehicle;
are the n-th mode’s equivalent resonant frequency, mass, and stiffness, respectively; and is the natural frequency of the vehicle. Note that when we solve Eq. (2), because
, we have
.
For our VBIS, local damage has relatively small influence on the bridge mode shape. Thus, we can approximate
celeration by analyzing Eq. (4). First, Eq. (4) has the following three components within three frequency bands:
• may encoder damage information in
. The dominant frequency of this term is
, which changes as vehicle properties change;
• may encode damage information in
and
. The dominant frequency of this term is
, which changes as bridge properties change;
• may encode damage information in both
and
, and the multiplication of the derivatives of
amplifies the damage information. Also, the dominant frequency of this term is
, which does not change as the bridge and vehicle properties change.
For each n, the third component is thus DS & DI once we multiply it by factor . Second, because of the local influence of the damage,
is non-stationary and its instantaneous frequency, which is defined as
, is time-varying and
. Third, for VBIS,
is very small and generally does not overlap with the vehicle frequency
and the shifting bridge frequency
. This property of the vehicle acceleration signal also indicates that the third component is a good candidate for our desired feature because it can be extracted from the original signal.
To extract the DS & DI feature from the vehicle acceleration signal, we can cast the IBHM problem as a signal decomposition problem. We decompose the non-stationary vehicle acceleration and extract the desired feature by reconstructing the non-stationary component within a frequency band that includes . We consider approximating the DS & DI feature for the first mode that is
Note
has time-varying instantaneous frequencies. Also,
is the sum of harmonic functions, hyperbolic functions and multiplication of harmonic and hyperbolic functions, so that the desired feature can be expressed as a superposition of IMFs.
Our proposed method uses SWT [16, 17] to represent the vehicle acceleration signal in the time-frequency plane and reconstruct the desired feature using ISWT. The reconstructed feature is used to estimate and localize damage. SWT has three advantages in solving our problem. First, it assumes that signals are approximately harmonic locally and have a slowly time-varying instantaneous frequency. This transform has the ability to decompose a non-stationary and time-varying signal as a superposition of IMF-type components. Second, comparing with the conventional time-frequency methods, such as short-time Fourier transform (STFT) and continuous wavelet transform (CWT), this empirical model decomposition (EMD)-like approach can further sharpen the time-frequency representation and enhance frequency localization [17]. Third, to localize the damage, we need to reconstruct the time-domain signal in the damage-related frequency band as the vehicle moves. ISWT can directly reconstruct a component within a selected frequency band and avoid the mode-mixing effect encountered by the EMD method [18]. The SWT has three steps:
1) Calculate the wavelet coefficients of the signal
where a is the scale, b is the time offset, and is the complex conjugate of wavelet. We use analytic Morlet wavelet; 2) Estimate the instantaneous frequencies for the signal
3) Reallocate the energy of CWT coefficients to enhance frequency localization. At discrete scales , the SWT of
is only calculated at the centers
of the frequency range
, where
. The SWT is
where . Component of the original signal in band
is estimated by ISWT:
where returns the real part of the function;
represents the Fourier transform of the complex conjugate of the wavelet. Further details about the implementation of SWT and ISWT can be found in [16, 17]. Our feature extraction algorithm is shown in Algorithm 1.
Finite element simulations and the dataset. Finite element models (FEMs) of the VBIS are employed to create a dataset. We have five bridges (Bridge 1, 2, 3, 4 & 5) simulated as Euler-Bernoulli beams of lengths 25 m, 25 m, 25 m, 20 m and 30 m, respectively; natural frequency of and
(for
), respectively; and they all have the same undamaged stiffness,
. We use the same properties (with reference to [9]) for the simulated vehicle:
, and v = 3 m/s.
Damage proxy is introduced by reducing the stiffness of one element at different locations. We use the same beam element with the length of 0.6 m for all bridges. The stiffness reduction () ranges from 70% to 30% with an interval of 10%. For each stiffness reduction level, the simulation is run at seven different damage locations (
is every eighth of the span), and for each damage scenario, the simulation is run ten times. For each simulation run, random forces, which follows a Gaussian distribution with zero mean and 0.1 N variance, are applied on each node of the FEMs. This added disturbing force is a process noise that propagates through time and varies in space. In total, we generate 5 (bridges)
5 (reduction levels)
(damage locations)+1 (undamaged bridge)]
(trials) = 2000 vehicle acceleration records.
Visualization of the proposed feature. We verify the DS & DI properties of our feature by visualization. The data for visualization is created by FEMs without adding the disturbing force. In our experiments, we use the proposed method to extract the desired feature within . Figures 2 (a) and (b) show vehicle accelerations and our proposed feature with different damage locations, respectively. Figures 2 (c) and (d) show the signals and the feature with different stiffness reductions. By visualization, we can easily localize and compare the damage simulated by reducing stiffness. To verify if the proposed feature is DI, in Figure 3 we visualize the ISWT of accelerations collected from the simulated vehicle traveling on Bridges 1 to 5 with different simulated damage before and after multiplying by
. We can observe that, after the multiplication, the features for different bridges having the same damage match each other, which indicates that our feature is DI. Though the results also contain boundary effects of the transform that require further investigation, it is easy to see that at the normalized locations of damage the proposed feature exhibit high sensitivity.
Fig. 2. Raw signals and our proposed feature for vehicle ac- celerations with different damage locations (a, and b) and different stiffness reduction levels (c and d). The blue marks indicate the damages. We can visually localize and compare the stiffness reductions using our proposed feature, which verifies that our feature is DS.
Fig. 3. ISWT within of vehicle accelerations traveling on Bridges 1, 2, 3, 4 & 5 with different damage (a) before and (b) after multiplying by
. Our features for bridges with different properties and the same damage match each other around damage locations. However, the features before applying the proposed multiplication do not match each other. This visualization verifies that our feature is DI.
Damage localization and quantification. To further verify if our proposed feature is DS & DI, we use a multi-task learning model proposed in [7] to estimate and localize stiffness reductions in a supervised and a semi-supervised fashion. The model’s input includes raw accelerations (raw data); band-pass filtered signals within , and
(Bandpass 1, 2&3); inverse CWT within the above three bands (ICWT 1, 2 & 3); the first three IMFs calculated using EMD (IMF 1, 2 & 3); ISWT within the three bands (Ours, ISWT 2 & 3); and spectrograms calculated using STFT, CWT and SWT.
Table 1. Damage localization and quantification results in terms of RMSE. Lower error means better performance.
The supervised task examines how sensitive each input is to damage. The second and third columns in Table 1 present this task’s results (30% for testing) in terms of root mean squared error (RMSE). Using the proposed feature, we obtain the best stiffness reduction estimation and localization results.
Table 1 also shows results for the semi-supervised regressions, where we have two sub-tasks: test if the feature is DI across bridges with different lengths (Different L) and with different natural frequencies (Different ). As shown in columns 4 to 7 of the table, using our feature, we obtain the best stiffness reduction estimations for the two sub-tasks, and the best damage localization results for bridges having different lengths and the same frequency. For damage localization in the second sub-task, IMF 1 provides the best result.
We introduce a physics-guided signal decomposition method to extract a DS & DI feature from vehicle accelerations for IBHM. The SWT is used to represent the data in the time-frequency plane, and the desired feature is reconstructed using ISWT within a damage-related frequency band. We verify and evaluate the DS & DI properties of the extracted feature for IBHM using the simulated data generated from FEMs. Among the six experiments we conducted, five of them exhibited the best damage quantification and localization results across different bridges using our proposed feature.
[1] ASCE, “Asce’s 2017 infrastructure report card,” 2017.
[2] George Lederman, Siheng Chen, James H Garrett Jr, Je- lena Kovaˇcevi´c, Hae Young Noh, and Jacobo Bielak, “Track-monitoring from the dynamic response of an operational train,” Mechanical Systems and Signal Processing, vol. 87, pp. 1–16, 2017.
[3] George Lederman, Siheng Chen, James H Garrett Jr, Je- lena Kovaˇcevi´c, Hae Young Noh, and Jacobo Bielak, “Track monitoring from the dynamic response of a passing train: a sparse approach,” Mechanical Systems and Signal Processing, vol. 90, pp. 141–153, 2017.
[4] George Lederman, Siheng Chen, James H Garrett Jr, Je- lena Kovaˇcevi´c, Hae Young Noh, and Jacobo Bielak, “A data fusion approach for track monitoring from multiple in-service trains,” Mechanical Systems and Signal Processing, vol. 95, pp. 363–379, 2017.
[5] Jingxiao Liu, Siheng Chen, George Lederman, David B Kramer, Hae Young Noh, Jacobo Bielak, James H Garrett Jr, Jelena Kovaˇcevi´c, and Mario Berg´es, “Dynamic responses, gps positions and environmental conditions of two light rail vehicles in pittsburgh,” Scientific data, vol. 6, no. 1, pp. 1–11, 2019.
[6] Jingxiao Liu, Yujie Wei, Mario Berg´es, Jacobo Bielak, James H Garrett Jr, and HaeYoung Noh, “Detecting anomalies in longitudinal elevation of track geometry using train dynamic responses via a variational autoencoder,” in Sensors and Smart Structures Technologies for Civil, Mechanical, and Aerospace Systems 2019. International Society for Optics and Photonics, 2019, vol. 10970, p. 109701B.
[7] Jingxiao Liu, Mario Berg´es, Jacobo Bielak, James H Garrett Jr, Jelena Kovaˇcevi´c, and Hae Young Noh, “A damage localization and quantification algorithm for indirect structural health monitoring of bridges using multi-task learning,” in AIP Conference Proceedings. AIP Publishing, 2019, vol. 2102, p. 090003.
[8] YB Yang and CW Lin, “Vehicle–bridge interaction dy- namics and potential applications,” Journal of sound and vibration, vol. 284, no. 1, pp. 205–226, 2005.
[9] YB Yang and Wei-Fan Chen, “Extraction of bridge fre- quencies from a moving test vehicle by stochastic subspace identification,” Journal of Bridge Engineering, vol. 21, no. 3, pp. 04015053, 2015.
[10] Thomas J Matarazzo and Shamim N Pakzad, “Scalable structural modal identification using dynamic sensor network data with stridex,” Computer-Aided Civil
and Infrastructure Engineering, vol. 33, no. 1, pp. 4–20, 2018.
[11] YB Yang, Yue-ju Li, and Kai Chun Chang, “Construct- ing the mode shapes of a bridge from a passing vehicle: a theoretical study,” Smart Structures and Systems, vol. 13, no. 5, pp. 797–819, 2014.
[12] YB Yang, Bin Zhang, Yanan Chen, Yao Qian, and Yun- tian Wu, “Bridge damping identification by vehicle scanning method,” Engineering Structures, vol. 183, pp. 637–645, 2019.
[13] George Lederman, Zihao Wang, Jacobo Bielak, Hae Young Noh, James H Garrett, Siheng Chen, Jelena Kovaˇcevi´c, Fernando Cerda, and Piervincenzo Rizzo, “Damage quantification and localization algorithms for indirect shm of bridges,” in Proc. Int. Conf. Bridge Maint., Safety Manag., Shanghai, China, 2014.
[14] Abdollah Malekjafarian, Fatemeh Golpayegani, Callum Moloney, and Siobh´an Clarke, “A machine learning approach to bridge-damage detection using responses measured on a passing vehicle,” Sensors, vol. 19, no. 18, pp. 4035, 2019.
[15] Jingxiao Liu, Siheng Chen, Mario Berg´es, Jacobo Bielak, James H Garrett, Jelena Kovaˇcevi´c, and Hae Young Noh, “Diagnosis algorithms for indirect structural health monitoring of a bridge model via dimensionality reduction,” Mechanical Systems and Signal Processing, vol. 136, pp. 106454, 2020.
[16] Ingrid Daubechies, Jianfeng Lu, and Hau-Tieng Wu, “Synchrosqueezed wavelet transforms: A tool for empirical mode decomposition,” arXiv preprint arXiv:0912.2437, 2009.
[17] Ingrid Daubechies, Jianfeng Lu, and Hau-Tieng Wu, “Synchrosqueezed wavelet transforms: An empirical mode decomposition-like tool,” Applied and computational harmonic analysis, vol. 30, no. 2, pp. 243–261, 2011.
[18] Danilo P Mandic, Naveed ur Rehman, Zhaohua Wu, and Norden E Huang, “Empirical mode decompositionbased time-frequency analysis of multivariate signals: The power of adaptive data analysis,” IEEE signal processing magazine, vol. 30, no. 6, pp. 74–86, 2013.
[19] Jingxiao Liu, Susu Xu, Mario Berg´es, Jacobo Bielak, James H Garrett Jr, and Hae Young Noh, “An expectation-maximization algorithm-based framework for vehicle-vibration-based indirect structural health monitoring of bridges,” Structural Health Monitoring 2019, 2019.