Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
A Polarization Independent Quasi-TEM Metamaterial Absorber for X and Ku Band Sensing Applications
Next Article in Special Issue
Damage Localization of Beam Bridges Using Quasi-Static Strain Influence Lines Based on the BOTDA Technique
Previous Article in Journal
Design and Implementation of Respiration Rate Measurement System Using an Information Filter on an Embedded Device
Previous Article in Special Issue
Model Updating for Nam O Bridge Using Particle Swarm Optimization Algorithm and Genetic Algorithm
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Drive-By Bridge Frequency Identification under Operational Roadway Speeds Employing Frequency Independent Underdamped Pinning Stochastic Resonance (FI-UPSR)

1
Department of Civil, Construction, and Environmental Engineering, The University of Alabama at Birmingham, 1075 13th St S, Birmingham, AL 35205, USA
2
School of Civil Engineering, University College Dublin, Newstead Block B, Belfield, Dublin D04V1W8, Ireland
*
Author to whom correspondence should be addressed.
Sensors 2018, 18(12), 4207; https://doi.org/10.3390/s18124207
Submission received: 2 November 2018 / Revised: 22 November 2018 / Accepted: 24 November 2018 / Published: 30 November 2018
(This article belongs to the Special Issue Bridge Structural Health Monitoring and Damage Identification)

Abstract

:
Recently, drive-by bridge inspection has attracted increasing attention in the bridge monitoring field. A number of studies have given confidence in the feasibility of the approach to detect, quantify, and localize damages. However, the speed of the inspection truck represents a major obstacle to the success of this method. High speeds are essential to induce a significant amount of kinetic energy to stimulate the bridge modes of vibration. On the other hand, low speeds are necessary to collect more data and to attenuate the vibration of the vehicle due to the roughness of the road and, hence, magnify the bridge influence on the vehicle responses. This article introduces Frequency Independent Underdamped Pinning Stochastic Resonance (FI-UPSR) as a new technique, which possesses the ability to extract bridge dynamic properties from the responses of a vehicle that passes over the bridge at high speed. Stochastic Resonance (SR) is a phenomenon where feeble information such as weak signals can be amplified through the assistance of background noise. In this study, bridge vibrations that are present in the vehicle responses when it passes over the bridge are the feeble information while the noise counts for the effect of the road roughness on the vehicle vibration. UPSR is one of the SR models that has been chosen in this study for its suitability to extract the bridge vibration. The main contributions of this article are: (1) introducing a Frequency Independent-Stochastic Resonance model known as the FI-UPSR and (2) implementing this model to extract the bridge vibration from the responses of a fast passing vehicle.

1. Introduction

The degradation in the structural integrity of highway bridges is attributed to aging and significant increases in freight volumes. In the United States, approximately 11% of bridges on the transport network has been classified as structurally defective [1]. This percentage is expected to spike at 25% within a ten-year period [1]. Notably, the most recent federal estimate for the cost of bridge rehabilitation projects is found to be $123 billion [2]. Therefore, the structural safety assessment of bridges on the road network has become an essential area of research. To this end, the detection of incipient defects has become a crucial issue where early detection and repair of damage could allow for prioritization of resources, reduce the retrofitting budget, and, simultaneously, preserve the functionality of the network. This has led the field of Structural Health Monitoring (SHM) identifying the condition of existing structures using sensor readings [3,4,5,6,7,8,9,10]. Conventional Structural Health Monitoring techniques require onsite instrumentation, which has disadvantages (e.g., necessity for a continuous source of power, data acquisition, and transmission electronics as well as sensor installation and maintenance) that have slowed their uptakes [11,12]. This is in addition to the impact on the capacity of the roadway due to the associated closure of bridge lanes [13,14,15]. Recently, researchers have investigated the feasibility of moving the instrumentation from the bridge structure to a passing vehicle to assess the bridge health condition, which has been referred to as “Drive-by Bridge Inspection” [16]. Yang et al. [17,18] were the first to publish a preliminary study investigating the feasibility of the approach. The researchers developed a mathematical expression for the Vehicle-Bridge Interaction (VBI) problem to explore the extraction of the fundamental bridge frequency from the vehicle’s response. They found that the response of the vehicle includes four main frequencies: (1) the driving frequency v/L (Hz) (where ‘v’ is the speed of the vehicle and ‘L’ is the bridge span), (2) the vehicle frequency fv (Hz) and (3) two frequencies, which are shifted by a value equal to half of the driving frequency from the bridge first natural frequency, fb, i.e., f b ± v 2 L (Hz). Furthermore, they declared that higher speeds have a positive impact on the approach since they increase the amplitude of the bridge frequency in the spectrum and, hence, increase the visibility of the bridge frequency components. In the same vein, they investigated utilizing the vehicle responses in detecting the change in the bridge damping ratio. They observed a significant drop in the spectrum power as the bridge damping ratio increased. The findings of this study suggested that the approach is feasible.
Afterwards, a number of studies were carried out studying the variation in the acceleration spectra of the vehicle response as a measure of bridge deterioration in the presence of the road roughness [16,19,20,21,22,23,24]. These found that road roughness induces a considerable vibration in the vehicle suspension system, which results in the signal being overwhelmed by the dynamic characteristics of the vehicle. These works concluded that the use of lower speeds for the inspection vehicle provides more data, attenuates the road roughness effect on the recorded signal, and, hence, the bridge responses become evident in the signal. In a study shown in Reference [6], the effect of road roughness on the vehicle response was substantially removed theoretically using a half car model by subtracting the acceleration signals of two consecutive axles. Under high speed, the bridge frequency has dominated the response spectrum of the subtracted signal. Elhattab et al. [25,26] have investigated the concept using an explicit VBI solver, which is included in a commercial Finite Element Program (LS-Dyna) and similar findings were observed. The approach was confined to numerical simulations and no field demonstration has been carried out. To probe the fidelity of the drive-by approach, Lin and Yang [27] utilized field test data to obtain the natural frequency of the Da-Wu-Lun Bridge in Taiwan from passing vehicle responses. They reported a successful identification of the bridge frequency for speeds below 40 km/h (24.85 mph). For higher speeds, the bridge frequency was swamped by vibrations arising from the roughness of the road. Oshima, Yamaguchi [28] showed the effectiveness of a vehicle equipped with a vibrator in extracting bridge frequencies. However, they recommended using lower speeds to attain higher resolution.
In summary, previous studies have identified the speed of the inspection vehicle as a key challenge to accurately identify the bridge frequencies. When the speed is increased, the vehicle’s excitation is increased due to the road roughness. Thus, the bridge response is less prominent due to excessive vehicle vibration. This article is devoted to the magnification of the feeble response of the bridge in the vehicle signal utilizing the Stochastic Resonance (SR) phenomenon.
Stochastic Resonance (SR) is a process of weak response detection in a signal overwhelmed by a global system response. It is one of the established mathematical models for nonlinear systems where generally feeble input information (such as a weak signal) can be amplified and optimized by exploiting the noise in the signal. The approach was proposed in 1981 by Benzi, Sutera [29] as an explanation of the observed periodicity in the ice ages on earth. Since then, the approach has spurred interest and many researchers have investigated the application of this hypothesis in different disciplines [30]. SR was originally built to work for bi-stable/overdamped systems. Bi-stable systems are systems in which the feeble information is oscillating between two stable positions. The overdamped term implies that the observed system’s stiffness is negligible in comparison with the surrounding domain. To illustrate the essence of these two underlying terminologies, a common example of SR implementation by crickets will be discussed. A cricket’s hair cells receive a signal from the surrounding noisy environment. Then they apply the SR technique to explore the possibility of a bird attack [31]. In this example, the feeble information is the periodic air pressure change caused by the predator’s wings, which varies between two stable positions (positive and negative as the wings go up and down). Since the influence of the pressure perturbation around the bird’s wing on the surrounding domain is negligible, the system is an overdamped one. Several studies have investigated the cooperative effect of noise in bi-stable/overdamped systems in detecting weak signals corrupted by a heavy background noise [32,33,34,35,36,37,38,39]. Following the publication of these work, researchers have investigated employing SR for tri-stable [40,41,42] and multi-stable [43] systems and underdamped/bi-stable systems [44,45,46,47,48]. These methods require a priori definition of the targeted frequency to extract the weak signal. Wang, He [49] introduced the Adaptive Multiscale Noise Tuning SR (AMSTSR) technique, which does not require prior knowledge of the targeted frequency. However, the approach is built for bi-stable systems. Recently, Zhang, He [50] introduced an SR model with the capability to detect mono-stable (i.e., the weak property is stable at one particular point) and bi-stable systems known as Underdamped Pinning Stochastic Resonance (UPSR). However, the approach requires identifying the frequency of interest in advance.
This paper is intended to pave the way for the exploitation of SR in the Civil Structural Health Monitoring domain with an emphasis on Drive-by Bridge Inspection. Presuming that most Civil Engineering structures are lightly dampened in nature and are mostly mono-stable systems (i.e., the structure tends to stabilize at a certain mode of vibration), the authors herein adopt the UPSR model for its suitability to the problem under study. The article will look at potential synergies of incorporating visual demarcation with the UPSR technique to obviate the need to define the frequency of weak signals. The proposed visualization approach will facilitate the observation of the SNR attribute of signals that lie in a specified frequency band where SNR (i.e., high SNR value) is exploited as a metric to discriminate between feeble signals and the background noise. The model stemming from this step will be introduced as Frequency Independent-UPSR (FI-UPSR). The proposed technique will be used to detect weak signals in a record disrupted by heavy noise. Afterwards, the approach will be implemented to extract bridge frequency from the response of a fast vehicle. The approach will be explored numerically using a VBI MATLAB model and experimentally using a full-scale field test.
The paper is organized as follows: Section 2 introduces Frequency Independent-UPSR (FI-UPSR) and the proposed visual demarcation approach. Section 3 demonstrates numerically the implementation of FI-UPSR in extracting bridge frequencies from the responses of a quickly passing vehicle. Lastly, a full-scale field test data is used to examine the fidelity of the proposed approach in the field. This will be presented in Section 4.

2. Frequency Independent Underdamped Pinning Stochastic Resonance (FI-UPSR)

2.1. Background to Underdamped Pinning Stochastic Resonance (UPSR)

The principle underlying the Stochastic Resonance phenomenon can be explained as follows: consider a heavily damped mass moving in a symmetric double potential well, V(x), as shown in Figure 1a. The mass is moving due to a fluctuating force f(t) contaminated by an additive noise having intensity, D. The application of a weak periodic force f(t) will tilt the potential well up and down. However, with a low intensity of force, it is not big enough to let the particle surmount the potential barrier and move from one potential minima to the other (Figure 1b). In contrast, the noise induces a more intense hopping for the particle around the potential minima (Figure 1c). By tuning the noise periodicity (or changing the potential geometry, V(x)), the hopping induced by the noise can be synchronized with the weak periodic force (Figure 1d) [30]. Since noise in most cases is the uncontrollable term, the SR adjusts the shape of its potential well (V(x)) to synchronize the hopping between potential minima by changing the parameters of the potential well. Within this framework, the weak periodic force f(t) refers to the targeted feeble information, which is the input weak signal (s(t)). This signal is swamped by a severe background noise (n(t)). The extraction of this feeble signal will be done by resonating the hopping between the two potential minima through changes in the parameters of the potential well (V(x)). The position of the particle at resonance (x(t)) represents the extracted feeble feature (herein the extracted weak signal), which is amplified when the hopping between potential minima is synchronized. The dynamic equation of this type of SR is presented in Equation (1).
d 2 x d t 2 = V ( x ) γ d x d t + { s ( t ) + n ( t ) }
where x ( t ) is the extracted weak signal, V ( x ) is the first derivative of the potential well (V(x)), γ is the system damping, and { s ( t ) + n ( t ) } is the input signal added with noise. Equation (1) implies that the Stochastic Resonance phenomenon is a single degree of freedom filtering process where feeble signals are amplified ( x ( t ) ) by approaching resonance. There exists a plethora of models to describe different types of potential shapes. Herein, we adopt the UPSR model for its suitability to the problem understudy.
The potential function of the UPSR model [50] is shown below.
V ( x ) = V 0 V d ( exp ( ( x + x 0 ) 2 L 2 ) + exp ( ( x x 0 ) 2 L 2 ) )
where V 0 is a constant and it will not be considered since it will have no effect after the function is differentiated, as shown in Equation (1), V d is the depth of the pinning, L is the length between the two pinnings, and ± x 0 is the center of each pinning. The potential function of the UPSR is governed by three parameters (i.e., V d , L, and x 0 ) that provide effective representation for mono-stable and bi-stable systems. Zhang, He [50] provide a numerical discretization for Equation (1) considering the potential function given in Equation (2), which is presented in Appendix A. As previously noted, to extract the weak signal ( x ( t ) ), the potential function is tuned until the noise is synchronized with the weak signal. In other words, when Equation (2) approaches resonance. At this point (i.e., when the signal and the noise are tuned), the extracted signal will have the highest Signal to Noise Ratio (SNR). Therefore, the discretized version of Equation (1) (i.e., Equation (A1) presented in Appendix A) will be solved several times for different potential parameters (i.e., Vd, L, and x 0 ) until this condition is achieved. The selection for the parameter ranges is done arbitrarily. Therefore, the ranges are changed randomly until the signal of the highest SNR is extracted. This framework is followed due to the lack of knowledge about the mathematical relationship between the values of the selected parameters (i.e., Vd, L, and x 0 ) and the characteristics of the extracted signal (e.g., frequency and amplitude). This is considered to be one of the major drawbacks of the SR filtering process. Section 2.2 of this article addresses this issue by deriving a direct relationship between the potential parameters of the UPSR model and the frequency of the extracted signal. This relationship is then used to establish a criterion for parameter initialization (i.e., Vd, L, and x 0 ).
The SNR of the extracted signal is calculated below [50].
S N R = 10 log 10 A s 0 N / 2 A i
where A s is the power of the targeted signal (i.e., the original one, s(t), rather than the extracted one, x(t)), and 0 N / 2 A i is the summation of the power of the computed signal x(t). The UPSR process is summarized in Figure 2. More details of the UPSR method are provided in Reference [50].
The preceding discussion deals with the equations, procedures, and terminologies of the UPSR technique. The following section presents the changes made to make the UPSR approach frequency independent.

2.2. Frequency Independent Underdamped Pinning Stochastic Resonance FI-UPSR

The core of the proposed FI-UPSR approach is to move from searching for the optimum potential parameters (i.e., Vd, L, and x 0 ) that maximize the SNR for a targeted signal to observe the SNR attributable to a specified frequency band where the signal that exhibits a relatively high SNR ratio in this band is the weak signal. To this end, the SNR values are explored in a three-dimensional plot to facilitate the identification of the highest SNR visually. The SNRs are computed by first generating a set of different combinations for the potential parameters. Specifically, at a certain potential depth (Vd), a range of potential shapes are generated (by varying L and x0). For each combination (i.e., Vd, L, and x 0 ), the corresponding signal is extracted by using Equation (A1) (Appendix A). Afterwards, the SNR of this signal is computed. The calculation of the SNR needs to be independent of the targeted signal power ( A s , which corresponds to the weak signal ( s ( t ) )). Therefore, Equation (3) has been modified to be as follows where A x is the maximum power of the computed signal ( x ( t ) ).
S N R = 10 log 10 A x 0 N / 2 A i
The selection of the potential parameter values and ranges is a substantial contributor to the success of the FI-UPSR approach. As highlighted, the FI-UPSR user will specify the upper and lower limits for the frequency band of interest [fminfmax]. Consequently, the corresponding potential parameter values for the selected frequency limits must be computed [ { V d ,   L ,   and   x 0 } f o r   f m i n { V d ,   L ,   and   x 0 } f o r   f m a x ] to generate a set of parameter combinations, which will be utilized to generate the SNR surface plot. However, there is no direct relation between the values of the potential parameters (i.e., Vd, L, and x0) and the frequency of the extracted signal ( x ( t ) ), which is one of the main drawbacks of conventional SR filters. Rather than relying on expert judgment on the selection of the proper potential values, we have derived a direct approximate relationship between the potential parameters and the frequency of the extracted signal.
The relationship is presented in Equation (5) of which the full derivation is provided in Appendix B.
f a p p r o x i m a t e R = V d π 2 L 2
where R is a rescaling factor that is used to rescale the time step in Equation (A1) and fapproximate is the approximate frequency of the extracted signal x(t).
The value of the potential depth plays a prominent role in the filtering process. Small potential depth (Vd) provides a small potential barrier. Thus, the particle jumps through the barrier without resistance. As Vd increases, surmounting the potential barrier becomes harder for the particle unless the noise and the force are tuned. As a result, high Vd values block unfavorable noise and provide better extraction for the weak signal. At a certain potential depth Vd, the values of L corresponding to a specific frequency range ( f 1 ,   f 2 , , f n ) can be computed using Equation (5). Lastly, the values of the pinning centers ( ± x 0 ) are related to the values of L where: L 2 x 0 . This relation will be explained in the following section after visualizing the SNR surface plot.
The rescaling factor R is a function of the time step size or the available sampling frequency and it impacts the accuracy of Equation (A1). The damping term in Equation (1) is used to smoothen the oscillation between the potentials to better recognize resonance when approached. The effect of the values of the rescaling factor R and the damping term, γ, on the filtering process will be explored later in the paper. Table 1 provides recommended values for V d , R, and γ that can be used as an initial trial. The values of L are calculated using Equation (5) and the values of x 0 are calculated using the L 2 x 0 relation.
The FI-UPSR procedure is summarized in the following Figure 3.

2.3. Investigating the FI-UPSR Approach

The feasibility of the proposed FI-UPSR method will be examined using a sinusoidal signal ( s = A sin ( 2 π f t ) where A = 0.5 and f = 10 Hz) contaminated with a Gaussian white noise of −15 dB in power. The signal record and the signal Power Spectral Density (PSD) are shown in Figure 4. It is evident that the signal has been totally disrupted by the noise. The contaminated signal has been processed using the proposed FI-UPSR approach. In this case, the damping coefficient was chosen as γ = 0.5, Vd was 100, the scaling factor, R was set to 50, and the sampling frequency, fs was 2000 Hz. Determining the frequency range is incumbent on the user where it is presumed that the frequency of interest should lie within a range that is known a priori. In this example, the scanned frequency window ranges from 5 to 80 Hz in increments of 0.1 Hz. The limits for the scanned window are chosen so the frequency of the pure signal (f = 10 Hz) lies between these limits. The parameter, x 0 , varies from −22.5 to 22.5 in increments of 0.05. The selection of this range will be discussed in detail later in the paper. The SNR surface is plotted by using a perceptual color map that includes hue, illumination, and saturation dimensions to provide a better visualization surface topography, as shown in Figure 5 [51,52]. The map has a luminance dimension, which provides brightness for the crests and darkness for troughs. Visual examination of the figure reveals a number of important characteristics. First, it is evident that the system is symmetric around x 0 = 0. Another point of note is that the left and the right dark red-colored regions of the figure (the triangles) have SNR values consistently close to zero. This attribute is associated with the geometry of the UPSR potential, which can be reshaped to provide a mono-stable or a bi-stable potential well. Zhang, He [50] showed that a bi-stable model is obtained when L 2 x 0 (boundary depicted in the figure with a dashed yellow line) while a mono-stable model is obtained otherwise. This area is in the first category (i.e., bi-stable model where L 2 x 0 ). This intriguing feature is essential to implement SR in health monitoring applications since most of the civil structures are mono-stable systems. We will focus on the mono-stable zone rather than visualizing the whole x0 domain. Figure 6a presents a surface plot of the mono-stable region of Figure 4 after eliminating the bi-stable zones. In the same vein, the SNR plot exhibits a significant peak wave within the boundaries of this mono-stable system where for x 0 = 0, fapproximate = 10.0 Hz. fapproximate is an approximate frequency for the extracted signal computed using Equation (5), but it does not represent the actual signal frequency. The actual signal frequency is identified by transferring the extracted signal, x(t), from the time domain to the frequency domain. To extract the weak signal, we manually select the parameter values from the peak wave and apply Equation (A1). Any point on the wave crest can be used. In this scenario, we selected the point where x 0 = 0 and the associated L is 15.92 (at fapproximate = 10 Hz in the figure). The extracted signal and the corresponding spectrum are shown in Figure 6b–c. From the results, it is evident that the proposed approach has an excellent capacity for detecting feeble signals.
Returning back to the damping term, γ, it smooths the oscillation of the particle between the potential minima. Low damping coefficients will result in a stronger oscillation that develops spurious peaks in the SNR-Surface plot and, therefore, observing the weak signal in the SNR-surface plot becomes challenging. On the other hand, high damping coefficients severely damp the particle movement, which might affect the properties of the signal. Accordingly, a probable selection for the damping term is important. This can be performed by testing the SNR plots for different damping values. For the previous problem, the damping is set equal to 0.05, 0.5, and 0.95. The results are shown in Figure 7. The scaling factor R has a similar effect. Figure 8 presents the SNR plots for γ = 0.5 when R = 1, 50, and 200. For bridge structures, the recommended value for γ is between 0.08–0.2. For other infrastructures, further investigation is needed to determine the recommended ranges for γ.

3. Extract Bridge Frequency from a Fast Passing Vehicle Signal Using FI-UPSR

To numerically examine the proposed method, a Vehicle Bridge Interaction model [53] has been adopted to simulate the vehicle crossing over the bridge. The VBI model has been built using the fully-integrated Euler-Bernoulli Beam element, which has two degrees of freedom per node (vertical displacement and out of plane rotation). Axial displacement is omitted to avoid membrane locking. The Implicit Newmark-Beta integration scheme is utilized to solve the VBI equation of motion. A numerical damping is added to the scheme to suppress instabilities of the higher mode responses. The time step is set equal to 0.0005 s to capture the changes in the vehicle acceleration when it traverses the bridge. The vehicle and the bridge properties are listed in Table 2 and Table 3, respectively. Figure 9 illustrates the adopted VBI model. The road roughness has been randomly generated using ISO-8608 [54] specification (Roughness class “A”). The vehicle speed is 25 m/s. The vehicle passes over a 10 m approach distance before crossing the bridge. The recorded acceleration has been contaminated by a 90 dB Gaussian Noise. This value is chosen to match the records of the sensor used in the field experiment [55]. The acceleration of the vehicle’s axle mass and its PSD are illustrated in Figure 10. As expected, the vehicle response overwhelms the spectrum.
The axle acceleration has been processed using the FI-UPSR technique. Herein, Vd = 100, R = 50, and γ = 0.1 while the values of x0 and L have been selected to provide a scanning window bounded by [4 Hz → 25 Hz]. The SNR surface plot is illustrated in Figure 11.
The first point of note is that the road roughness and noise effects have been substantially removed and few peaks were kept in the SNR plot. In Figure 11a, the values of fapproximate corresponding to the highest SNR are shown. The exact signals frequencies are identified by extracting the weak signals that correspond to the SNR peaks. Then plot the signal PSD depicted in Figure 12.
Two values have been determined for the first bridge frequency. Both values relate to the first bridge frequency, shifted by ± half the driving-frequency [17,18] (5.67 ± 25 2 × 15 = 4.84 and 6.5). It is evident in the plot that the 6.5 Hz frequency has more power than the 4.84 Hz frequency. An ambiguous dominant peak is visible in the plot and is associated with a frequency of 13.32 Hz. This is apparently due to the influence of the road roughness. Thus, the approach could not totally attenuate the road roughness effect at higher frequencies.
The results suggest that the FI-UPSR is capable of detecting the lower bridge frequencies from a fast crossing vehicle measurement. However, the existence of an obscure frequency (13.46 Hz) requires a plausible interpretation to firmly root the approach before expanding it to the field scale. The authors presume that this is due to the nature of the representation for the contact between the vehicle and the road. The simplified Quarter-Car model touches the road at a single point. In tandem with that, the used road roughness is defined each 2.5 cm along the vehicle path. Thus, the vehicle will observe a stronger vibration due to the instantaneous change in the profile. This vibration produces multiple local maxima in the SNR surface plot especially at the higher frequencies, as presented in Figure 13a. In the field, on the other hand, the contact area between the vehicle tires and road roughness will attenuate the roughness effect on the vehicle signal. To simulate the condition in the field, the road roughness has been regenerated using a 25-cm increment. Figure 14 presents the PSD plots for the displacement history of the quarter car contact point and for the two scenarios (the 2.5-cm increment in Figure 14a and the 25-cm increment in Figure 14b). It is evident that increasing the length of the contact attenuates the effect of the road roughness on the vehicle vibration. The FI-UPSR surface plot has been regenerated utilizing the new vehicle acceleration. The results are presented in Figure 13b. The figure shows that, for milder profiles, local maxima have less power. This point requires further investigation to clearly depict the effect of the contact area between the vehicle and the road on the Drive-by Bridge Inspection approaches.
Regardless of this point, the influence of road roughness in exciting the vehicle and the dominance of the vehicle frequencies was overcome in this example for the first and the second bridge frequency.

4. FI-UPSR Fidelity for Full Scale Field Test Data

The field test was carried out on a skewed pre-stressed bridge consisting of three simply supported spans. The bridge is located on HWY 113 in Bartow County, Atlanta, Georgia between Covered Bridge Rd. and Dry Creek Rd. Each span is 21.3 m from the centers of the two supports. The bridge deck is a reinforced concrete slab integral with five pre-stressed concrete girders. The roadway facility consists of two-lane one-way traffic and one hard shoulder. In addition to vehicle instrumentation, sensors were installed on the first span of the bridge (from the traffic direction) as shown in Figure 15a,b. The inspection vehicle used was instrumented with three Silicon Design 2012-002 accelerometers on its axles. A special accelerometer, which was developed by Dong, Zhu [56], was used to measure the body mass acceleration. This sensor was mounted near the vehicle’s center of gravity, which was illustrated in Figure 15c,d. At the same location, a gyroscope (ADXRS624) was installed to monitor the vehicle pitching motion. All the recorded data during the test were sent wirelessly to a DAQ station, which lies under the bridge. Furthermore, the bridge barriers were equipped with five laser emitters, which were utilized to determine the exact time when the vehicle entered and exited the bridge.
The bridge natural frequencies were first identified by using a vibration test. This used the ELECTRO-SEIS Long Stroke Exciter vibration shaker while the bridge acceleration was recorded using 15 accelerometers installed on the girders (three sensors spaced equally on each of the bridge’s five girders). The sampling frequency was set to 100 Hz while the frequency resolution was 6.25 × 10−3 Hz. The traffic was blocked during the test to halt the traffic’s mass interference with the calculated frequency, which was noted by Yang, Cheng [57]. The test was repeated twice, and the spectra are illustrated in Figure 16. The bridge frequencies are listed in Table 4.
The inspection truck crossed the bridge three times at the speeds listed in Table 5. The truck could not reach the road’s speed limit (90 km/h or 60 mph) since the available approach distance was limited and the truck did not have enough time to accelerate. However, the maximum speed (i.e., for Test 3) is approximately 30% higher than the speed where the bridge frequencies are detectable (i.e., 40 km/h), which was previously introduced in the literature. Each test has two stages. (i) Stationary stage with engine running and (ii) moving stage where the vehicle accelerates to cross the bridge. The recorded accelerations when the truck was stationary and excited only by the vibration of the engine have been used to determine the frequencies of the vehicle mechanical system. The PSD for body mass acceleration (i.e., Accelerometer #4) is presented in Figure 17. The vehicle frequencies are identified using the acceleration spectrum. It is important to identify the vehicle frequencies to distinguish the bridge frequencies in the FI-UPSR plots. The frequency band is set to [5 Hz → 15 Hz]. The proposed FI-UPSR technique was utilized to process the data of the three tests for this frequency band. The SNR plots for Test 3 are presented in Figure 18a.
Three dominant peaks are observed in the surface plot. The peaks at 5 Hz and 10.5 Hz indicate the frequencies of the vehicle. The third peak refers to the first fundamental frequency for the bridge (the peak at fapproximate = 7.48 Hz). To extract the bridge signal, we manually selected the potential parameter values from the crest of the peak-wave and then we applied Equation (A1). Any point on the crest of the peak wave can be utilized. For instance, Point #1 where x0 = 0 and L = 21.28 or Point #2 where x0 = 5.58 and L = 17.51. In this example, we utilized the coordinates of Point #1. The PSD for the extracted signal is presented in Figure 18b. The figure reveals a successful identification for the first bridge frequency. The same surface plot attribute is observed for Test 1 and Test 2. The results for Test 1 and Test 2 are presented in Figure 19a,b.
The results presented in the figure show a successful identification for the first bridge frequency for all three tests. The second frequency has not been observed in any of the surface plots. This may be because the second mode of vibration (i.e., a transverse mode of vibration) interferes less with the body mass vibration. The third frequency has not been successfully extracted from the vehicle responses. The authors presume that the frequency content of the vehicle (10.00 Hz) and the bridge frequency (11.93 Hz) have been fused to provide an average peak between them. This is shown in Figure 19 where the peak takes place between 10.0 Hz to 12.00 Hz.
While the results of the field test demonstrated a successful identification for the first bridge frequency, the SNR-surface plot exhibit multiple local maxima due to the vehicle frequencies that make the bridge frequency identification more challenging. The authors suggest that adopting a truck trailer model where the sensors are installed in the trailer and the trailer suspension properties are chosen carefully to provide low-frequency content may minimize the interference of the vehicle frequencies and positively affect the approach.
Another point of note is the impact of the environmental condition on the accuracy of assessing the bridge frequencies utilizing the drive-by approach. It is often the influence of temperature on the stiffness of the bridge, which results in changes in natural frequency. This point can be addressed by finding or estimating the temperature of the bridge and, hence, applying a correction where the temperature equipment could be mounted on the vehicle to perform the correction instantaneously. Another option is to mitigate the influence of temperature by choosing the best times of day for screening to minimize temperature differences from one test to the other. However, these issues present when utilizing short-time measurements of frequency to monitor the bridge condition.
While the approach holds significant merits and its application can be expanded to other SHM problems, the utilization of a visual peak picking process subjects the results to some uncertainty. The process requires a priori experience about the SNR surface plot shapes to properly identify the peaks. Figure 20b demonstrates the PSD for the extracted signal for Test #3 where the point is taken from the wave trough rather than the wave crest (Point #3 at Figure 20a). As presented in the figure, spurious peaks start to dominate the PSD plot and, therefore, it becomes unviable for identifying the bridge frequency. Further work is needed to implement a suitable peak-picking algorithm to automatically identify the peaks within a certain frequency band.

5. Conclusions

Recently, the drive-by bridge inspection has become a burgeoning field of research in the bridge health-monitoring domain. However, it faces a major obstacle in the need for low speeds of the inspection vehicle. This paper introduces Frequency Independent Underdamped Pinning Stochastic Resonance (FI-UPSR) as a new technique to extract the dynamic characteristics of the bridge from the responses of a passing vehicle operating with a high speed. Stochastic Resonance is a process for the detection of a feeble property in a very rough domain through the addition of noise. This paper invests the UPSR model developed by Zhang, He [50] as a suitable model for Civil structures since it works with mono-stable/underdamped systems, which are features of most of these structures. The UPSR model requires a priori knowledge of the targeted frequency to extract the waveform associated with it. Therefore, the authors have introduced a graphical approach that can be utilized to identify feeble signals in a response disrupted by severe noise. The new approach is based on the plotting of 2D surface plots of SNR values for a specified range of UPSR parameters. By visualizing the plots, the values of the parameters that provide the maximum SNR can be observed and, hence, the weak signals can be extracted. The article has explored numerically the use of FI-UPSR for the extraction of feeble bridge vibrations from a fast passing vehicle and the results have revealed a successful identification of the first and second bridge frequencies. A numerical VBI MATLAB model has been used to simulate the crossing of the vehicle over the bridge at 25 m/s (90 km/h).
To investigate the fidelity of the approach in practice, a full-scale field test has been utilized. The experiment employed an inspection truck equipped with an accelerometer on its center of gravity to extract the fundamental bridge frequencies. First, the data from 15 accelerometers installed on the bridge were used to identify its fundamental frequencies. Afterwards, the inspection truck crossed the bridge at three different speeds. The data of the truck accelerometers have been processed using the FI-UPSR method to identify the bridge frequencies. The results show a successful identification of the first natural frequency for the bridge of all test speeds. The second mode of vibration for the bridge is orthogonal to the traffic direction and, therefore, it was not evident in the vehicle body mass acceleration history.
The proposed FI-UPSR method shows good potential to amplify a weak property in a noisy disrupted signal, which could have significant applications in Civil SHM. In conclusion, the paper presents an alternative methodology to process the data of the vehicle sensors by using rather than suppressing unwanted noise in the recorded signal. Further investigation is required on the impact of environmental effects on the proposed approach.

Author Contributions

All listed authors contributed extensively to the work. A.E. developed the method, wrote the method code, and applied it to the numerical and the experimental data. N.U. is the Principal Investigator of the NSF research project, lead and supervised the research. E.OB. provided advice on drafting and consulted on the technical issues.

Funding

This research was funded by the National Science of Foundation, grant number NSF-CNS-1645863.

Acknowledgments

The authors wish to express their gratitude for the financial support received from the National Science of Foundation (NSF-CNS-1645863) and the Science Foundation Ireland under the US/Ireland project, 16/US/I3277. Any opinions, findings, and conclusions or recommendations expressed in this publication are those of the authors and do not necessarily reflect the view of the sponsors.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

The appendix presents the discretized version of Equation (1) developed by Zhang, He [50].
                     V ( x ) = d V ( x ) d x = V d ( 2 ( x + x 0 ) L 2 exp ( ( x + x 0 ) 2 L 2 ) + 2 ( x x 0 ) L 2 exp ( ( x x 0 ) L 2 ) ) x ( 0 ) = 0 ;   d x d t ( 0 ) = 0 For   i = 2 :   end k x 1 = d x d t ( i 1 ) k y 1 = V ( x ( i 1 ) ) γ d x d t ( i 1 ) + { s ( i 1 ) + n ( i 1 ) }   k x 2 = d x d t ( i 1 ) + h 2 k y 1 k y 2 = V ( x ( i 1 ) + h 2 k x 1 ) γ ( d x d t ( i 1 ) + h 2 k y 1 ) + { s ( i 1 ) + n ( i 1 ) }   k x 3 = d x d t ( i 1 ) + h 2 k y 2 k y 3 = V ( x ( i 1 ) + h 2 k x 2 ) γ ( d x d t ( i 1 ) + h 2 k y 2 ) + { s ( i ) + n ( i ) }   k x 4 = d x d t ( i 1 ) + h k y 3 k y 4 = V ( x ( i 1 ) + h k x 3 ) γ ( d x d t ( i 1 ) + h k y 3 ) + { s ( i ) + n ( i ) }   x ( i ) = x ( i 1 ) + h 6 ( k x 1 + 2 k x 2 + 2 k x 3 + k x 4 ) d x d t ( i ) = d x d t ( i 1 ) + h 6 ( k y 1 + 2 k y 2 + 2 k y 3 + k y 4 )
where h is the calculation step and equals R*dt (dt is the time step and R is a rescaling factor). The parameter, x(t), is the targeted weak property or, in this case, the weak signal.

Appendix B

This section presents the derivation of the approximate relationship between the approximate frequency fapproximate and the potential parameters of the UPSR model. The original SR equation is shown below.
d 2 x d t 2 = V ( x ) γ d x d t + { s ( t ) + n ( t ) }
d 2 x d t 2 + γ d x d t + V ( x ) = { s ( t ) + n ( t ) }
Knowing that:
V ( x ) = V 0 V d ( exp ( ( x + x 0 ) 2 L 2 ) + exp ( ( x x 0 ) 2 L 2 ) )
Then,
V ( x ) = d V ( x ) d x = V d ( 2 ( x + x 0 ) L 2 exp ( ( x + x 0 ) 2 L 2 ) + 2 ( x x 0 ) L 2 exp ( ( x x 0 ) L 2 ) )
The SNR plots in Section 2 and Section 3 illustrate the essential features of the UPSR model. First, for the mono-stable zone, the SNR plot is symmetric around x0 = 0. In addition, the SNR peaks are more evident around x0 = 0. Furthermore, Equation (A5) shows that, for specific values of x0, increasing Vd requires an increase in L to extract the same signal x(t) by balancing the V d / L 2 ratio. Accordingly, use x0 = 0 and a considerably high value for L to simplify the formula. Thus, Equation (A5) yields the following.
V ( x ) 4 V d L 2 x exp ( x 2 L 2 )
For   large   L ,   exp ( x 2 L 2 ) 1
Hence,
V ( x ) 4 V d L 2 x
Substituting in Equation (A3) then gives the equation below.
d 2 x d t 2 + γ d x d t + 4 V d L 2 x = { s ( t ) + n ( t ) }
The equation yields to the form of a forced damped vibration for a single degree of freedom system where the system will resonate when the forcing frequency ( { s ( t ) + n ( t ) } ) matches the system frequency. The system frequency for the left hand side, neglecting the damping term for simplicity, equals the following.
ω = 2 π f a p p r o x i m a t e = 4 V d L 2
Thus:
f a p p r o x i m a t e = V d π 2 L 2
It is worth mentioning that the time step for the discretized equation, as depicted in Equation (A1), is rescaled by a factor R. To consider this, the frequency needs to be re-scaled with 1/R, which is shown below.
f a p p r o x i m a t e R = V d π 2 L 2

References

  1. Davis, S.; DeGood, K.; Donohue, N.; Goldberg, D. The Fix We’re in for: The State of Our Nation’s Busiest Bridges; Transportation for America: Washington, DC, USA, 2013. [Google Scholar]
  2. Federal Highway Administration. 2015 Status of the Nation’s Highways, Bridges, and Transit Conditions & Performance Report to Congress; 016093706X; Federal Highway Administration: Washington, DC, USA, 2017.
  3. Stanbridge, A.; Khan, A.; Ewins, D. Fault identification in vibrating structures using a scanning laser doppler vibrometer. In Proceedings of the International Workshop on Structural Health Monitoring, Stanford, CA, USA, 12–14 September 2017; pp. 56–65. [Google Scholar]
  4. Zhang, L.; Zhao, H.; OBrien, E.J.; Shao, X. Virtual Monitoring of orthotropic steel deck using bridge weigh-in-motion algorithm: Case study. Struct. Health Monit. 2018. [Google Scholar] [CrossRef]
  5. Zang, C.; Friswell, M.I.; Imregun, M. Structural Damage Detection using Independent Component Analysis. Struct. Health Monit. 2004, 3, 69–83. [Google Scholar] [CrossRef]
  6. Keenahan, J.; Obrien, E.J.; McGetrick, P.J.; Gonzalez, A. The use of a dynamic truck–trailer drive-by system to monitor bridge damping. Struct. Health Monit. 2013, 13, 143–157. [Google Scholar] [CrossRef] [Green Version]
  7. Tong, G.; Aiqun, L.; Jianhui, L. Fatigue Life Prediction of Welded Joints in Orthotropic Steel Decks Considering Temperature Effect and Increasing Traffic Flow. Struct. Health Monit. 2008, 7, 189–202. [Google Scholar] [CrossRef]
  8. Fan, W.; Qiao, P. Vibration-based damage identification methods: A review and comparative study. Struct. Health Monit. 2011, 10, 83–111. [Google Scholar] [CrossRef]
  9. Carden, E.P.; Fanning, P. Vibration based condition monitoring: A review. Struct. Health Monit. 2004, 3, 355–377. [Google Scholar] [CrossRef]
  10. Chang, P.C.; Flatau, A.; Liu, S. Health monitoring of civil infrastructure. Struct. Health Monit. 2003, 2, 257–267. [Google Scholar] [CrossRef]
  11. Yang, Y.; Chang, K. Extracting the bridge frequencies indirectly from a passing vehicle: Parametric study. Eng. Struct. 2009, 31, 2448–2459. [Google Scholar] [CrossRef]
  12. Malekjafarian, A.; McGetrick, P.J.; OBrien, E.J. A Review of Indirect Bridge Monitoring Using Passing Vehicles. Shock Vib. 2015, 2015, 286139. [Google Scholar] [CrossRef] [Green Version]
  13. Ramadan, O.; Sisiopiku, V. Impact of Bottleneck Merge Control Strategies on Freeway Level of Service. Transp. Res. Procedia 2016, 15, 583–593. [Google Scholar] [CrossRef]
  14. Ramadan, O.E.; Sisiopiku, V.P. Bottleneck Merge Control Strategies for Work Zones: Available Options and Current Practices. Open J. Civ. Eng. 2015, 5, 428–436. [Google Scholar] [CrossRef]
  15. Ramadan, O.E.; Sisiopiku, V.P. Evaluation of merge control strategies at interstate work zones under peak and off-peak traffic conditions. J. Transp. Technol. 2016, 6, 118–130. [Google Scholar] [CrossRef]
  16. Kim, C.W.; Kawatani, M. Challenge for a drive-by bridge inspection. In Proceedings of the 10th International Conferenceon Structural Safety and Reliability, ICOSSAR, Osaka, Japan, 13–17 September 2009; pp. 758–765. [Google Scholar]
  17. Yang, Y.; Lin, C. Vehicle–bridge interaction dynamics and potential applications. J. Sound Vib. 2005, 284, 205–226. [Google Scholar] [CrossRef]
  18. Yang, Y.-B.; Lin, C.; Yau, J. Extracting bridge frequencies from the dynamic response of a passing vehicle. J. Sound Vib. 2004, 272, 471–493. [Google Scholar] [CrossRef]
  19. McGetrick, P.J.; González, A.; OBrien, E.J. Theoretical investigation of the use of a moving vehicle to identify bridge dynamic parameters. Insight Non-Destr. Test. Cond. Monit. 2009, 51, 433–438. [Google Scholar] [CrossRef] [Green Version]
  20. Toshinami, T.; Kawatani, M.; Kim, C. Feasibility investigation for identifying bridge’s fundamental frequencies from vehicle vibrations. In Proceedings of the Fifth International IABMAS Conference on Bridge Maintenance, Safety, Management and Life-Cycle Optimization, Philadelphia, PA, USA, 11–15 July 2010; p. 108. [Google Scholar]
  21. Keenahan, J.; McGetrick, P.; O’Brien, E.J.; Gonzalez, A. Using instrumented vehicles to detect damage in bridges. In Proceedings of the 15th International Conference on Experimental Mechanics, Porto, Portugal, 22–27 July 2012. Paper No. 2934. [Google Scholar]
  22. O’Brien, E.J.; Keenahan, J.; McGetrick, P.; González, A. Using Instrumented Vehicles to Detect Damage in Bridges. In Proceedings of the BCRI 12-Bridge and Concrete Research in Ireland, Dublin, Ireland, 6–7 September 2012. [Google Scholar]
  23. González, A.; OBrien, E.J.; McGetrick, P.J. Detection of bridge dynamic parameters using an instrumented vehicle. In Proceedings of the 5th World Conference on Structural Control and Monitoring, Tokyo, Japan, 12–14 July 2010. [Google Scholar]
  24. McGetrick, P.J.; González, A.; O’Brien, E.J. Monitoring bridge dynamic behaviour using an instrumented two axle vehicle. In Proceedings of the Bridge and Infrastructure Research in Ireland 2010 (BRI 10), Cork, Ireland, 2–3 September 2010. [Google Scholar]
  25. ElHattab, A.; Uddin, N.; OBrien, E. Drive-by bridge damage detection using non-specialized instrumented vehicle. Bridge Struct. 2017, 12, 73–84. [Google Scholar] [CrossRef]
  26. El-hattab, A.; Uddin, N.; OBrien, E. Drive-By Bridge Damage Detection Using Apparent Profile. In Proceedings of the First International Conference on Advances in Civil Infrastructure and Construction Materials (CISM), Dhaka, Bangladesh, 14–15 December 2015; pp. 33–45. [Google Scholar]
  27. Lin, C.; Yang, Y. Use of a passing vehicle to scan the fundamental bridge frequencies: An experimental verification. Eng. Struct. 2005, 27, 1865–1878. [Google Scholar] [CrossRef]
  28. Oshima, Y.; Yamaguchi, T.; Kobayashi, Y.; Sugiura, K. Eigenfrequency estimation for bridges using the response of a passing vehicle with excitation system. In Proceedings of the Fourth International Conference on Bridge Maintenance, Safety and Management, Seoul, Korea, 13–17 July 2008; pp. 3030–3037. [Google Scholar]
  29. Benzi, R.; Sutera, A.; Vulpiani, A. The mechanism of stochastic resonance. J. Phys. A Math. Gen. 1981, 14, L453. [Google Scholar] [CrossRef]
  30. Gammaitoni, L.; Hänggi, P.; Jung, P.; Marchesoni, F. Stochastic resonance. Rev. Mod. Phys. 1998, 70, 223. [Google Scholar] [CrossRef]
  31. Levin, J.E.; Miller, J.P. Broadband neural encoding in the cricket cereal sensory system enhanced by stochastic resonance. Nature 1996, 380, 165–168. [Google Scholar] [CrossRef] [PubMed]
  32. Chapeau-Blondeau, F. Stochastic resonance and optimal detection of pulse trains by threshold devices. Digit. Signal Process. 1999, 9, 162–177. [Google Scholar] [CrossRef]
  33. Jung, P.; Hänggi, P. Amplification of small signals via stochastic resonance. Phys. Rev. A 1991, 44, 8032–8042. [Google Scholar] [CrossRef] [PubMed]
  34. He, H.-L.; Wang, T.-Y.; Leng, Y.-G.; Zhang, Y.; Li, Q. Study on non-linear filter characteristic and engineering application of cascaded bistable stochastic resonance system. Mech. Syst. Signal Process. 2007, 21, 2740–2749. [Google Scholar] [CrossRef]
  35. He, Q.; Wang, J. Effects of multiscale noise tuning on stochastic resonance for weak signal detection. Digit. Signal Process. 2012, 22, 614–621. [Google Scholar] [CrossRef]
  36. Han, D.; An, S.; Shi, P. Multi-frequency weak signal detection based on wavelet transform and parameter compensation band-pass multi-stable stochastic resonance. Mech. Syst. Signal Process. 2016, 70, 995–1010. [Google Scholar] [CrossRef]
  37. Chen, Z.; Liu, J.; Zhan, C.; He, J.; Wang, W. Reconstructed Order Analysis-Based Vibration Monitoring under Variable Rotation Speed by Using Multiple Blade Tip-Timing Sensors. Sensors 2018, 18, 3235. [Google Scholar] [CrossRef] [PubMed]
  38. Lai, Z.-H.; Leng, Y.-G. Generalized parameter-adjusted stochastic resonance of duffing oscillator and its application to weak-signal detection. Sensors 2015, 15, 21327–21349. [Google Scholar] [CrossRef] [PubMed]
  39. Huang, Q.; Liu, J.; Li, H. A modified adaptive stochastic resonance for detecting faint signal in sensors. Sensors 2007, 7, 157–165. [Google Scholar] [CrossRef]
  40. Lu, S.; He, Q.; Zhang, H.; Zhang, S.; Kong, F. Note: Signal amplification and filtering with a tristable stochastic resonance cantilever. Rev. Sci. Instrum. 2013, 84, 026110. [Google Scholar] [CrossRef] [PubMed]
  41. Pollak, E.; Talkner, P. Kramers turnover theory for a triple well potential. Acta Phys. Pol. B 2001, 32, 361. [Google Scholar]
  42. Zhang, H.; Kong, F.; Lu, S.; He, Q. A Tri-Stable stochastic resonance model and its applying in detection of weak signal. In Proceedings of the 2013 Fourth International Conference on Intelligent Control and Information Processing (ICICIP), Beijing, China, 9–11 June 2013; pp. 199–204. [Google Scholar]
  43. Li, J.; Chen, X.; He, Z. Multi-stable stochastic resonance and its application research on mechanical fault diagnosis. J. Sound Vib. 2013, 332, 5999–6015. [Google Scholar] [CrossRef]
  44. Alfonsi, L.; Gammaitoni, L.; Santucci, S.; Bulsara, A. Intrawell stochastic resonance versus interwell stochastic resonance in underdamped bistable systems. Phys. Rev. E 2000, 62, 299–302. [Google Scholar] [CrossRef]
  45. Almog, R.; Zaitsev, S.; Shtempluck, O.; Buks, E. Signal amplification in a nanomechanical Duffing resonator via stochastic resonance. Appl. Phys. Lett. 2007, 90, 013508. [Google Scholar] [CrossRef] [Green Version]
  46. Kang, Y.-M.; Xu, J.-X.; Xie, Y. Observing stochastic resonance in an underdamped bistable Duffing oscillator by the method of moments. Phys. Rev. E 2003, 68, 036123. [Google Scholar] [CrossRef] [PubMed]
  47. Ray, R.; Sengupta, S. Stochastic resonance in underdamped, bistable systems. Phys. Lett. A 2006, 353, 364–371. [Google Scholar] [CrossRef] [Green Version]
  48. Xu, Y.; Wu, J.; Zhang, H.-Q.; Ma, S.-J. Stochastic resonance phenomenon in an underdamped bistable system driven by weak asymmetric dichotomous noise. Nonlinear Dyn. 2012, 70, 531–539. [Google Scholar] [CrossRef]
  49. Wang, J.; He, Q.; Kong, F. Adaptive multiscale noise tuning stochastic resonance for health diagnosis of rolling element bearings. IEEE Trans. Instrum. Meas. 2015, 64, 564–577. [Google Scholar] [CrossRef]
  50. Zhang, H.; He, Q.; Kong, F. Stochastic resonance in an underdamped system with pinning potential for weak signal detection. Sensors 2015, 15, 21169–21195. [Google Scholar] [CrossRef] [PubMed]
  51. Rogowitz, B.E.; Treinish, L.A. Why Should Engineers and Scientists Be Worried About Color? 1996. Available online: http://www.research.ibm.com/people/l/lloydt/color/color.htm (accessed on 18 March 2010).
  52. Light, A.; Bartlein, P.J. The end of the rainbow? Color schemes for improved data graphics. Eos Trans. Am. Geophys. Union 2004, 85, 385–391. [Google Scholar] [CrossRef]
  53. Elhattab, A.; Uddin, N.; OBrien, E. Drive-by bridge damage monitoring using Bridge Displacement Profile Difference. J. Civ. Struct. Health Monit. 2016, 6, 839–850. [Google Scholar] [CrossRef]
  54. ISO-8608. Mechanical Vibration-Road Surface Profiles-Reporting of Measured Data; International Organization for Standardization (ISO): Geneva, Switzerland, 1995. [Google Scholar]
  55. Wang, Y.; Uddin, N.; Jacobs, L.J.; Kim, J.-Y. Field Validation of a Drive-By Bridge Inspection System with Wireless BWIM+ NDE Devices; National Center for Transportation Systems Productivity and Management: Atlanta, GA, USA, 2016. [Google Scholar]
  56. Dong, X.; Zhu, D.; Wang, Y.; Lynch, J.P.; Swartz, R.A. Design and validation of acceleration measurement using the Martlet wireless sensing system. In Proceedings of the ASME 2014 Conference on Smart Materials, Adaptive Structures and Intelligent Systems, Newport, RI, USA, 8–10 September 2014; p. V001T005A006. [Google Scholar]
  57. Yang, Y.; Cheng, M.; Chang, K. Frequency variation in vehicle–bridge interaction systems. Int. J. Struct. Stab. Dyn. 2013, 13, 1350019. [Google Scholar] [CrossRef]
Figure 1. Conventional Overdamped Bi-stable Stochastic Resonance Model. (a) Model terminology and (b) particle perturbation by the force f (c) particle perturbation by the force f and noise D in the general case. (d) Particle perturbation by the force f and the noise D when resonance is approached.
Figure 1. Conventional Overdamped Bi-stable Stochastic Resonance Model. (a) Model terminology and (b) particle perturbation by the force f (c) particle perturbation by the force f and noise D in the general case. (d) Particle perturbation by the force f and the noise D when resonance is approached.
Sensors 18 04207 g001
Figure 2. UPSR process for extracting the feeble signal at a specific frequency [50].
Figure 2. UPSR process for extracting the feeble signal at a specific frequency [50].
Sensors 18 04207 g002
Figure 3. The FI-UPSR procedure.
Figure 3. The FI-UPSR procedure.
Sensors 18 04207 g003
Figure 4. The raw and contaminated signals and their spectra. (a) Signals record; (b) Signals PSD.
Figure 4. The raw and contaminated signals and their spectra. (a) Signals record; (b) Signals PSD.
Sensors 18 04207 g004
Figure 5. SNR-surface plot for the extracted signals at different potential geometries.
Figure 5. SNR-surface plot for the extracted signals at different potential geometries.
Sensors 18 04207 g005
Figure 6. FI-UPSR results (a) SNR surface plot, (b) extracted weak signal, and (c) PSD plot for the extracted signal.
Figure 6. FI-UPSR results (a) SNR surface plot, (b) extracted weak signal, and (c) PSD plot for the extracted signal.
Sensors 18 04207 g006
Figure 7. Influence of damping on the FI-UPRS method. (a) γ = 0.05 (b) γ = 0.5 (c) γ = 0.95.
Figure 7. Influence of damping on the FI-UPRS method. (a) γ = 0.05 (b) γ = 0.5 (c) γ = 0.95.
Sensors 18 04207 g007
Figure 8. Influence of the scaling factor on the FI-UPRS method. (a) R = 1 (b) R = 50 (c) R = 200.
Figure 8. Influence of the scaling factor on the FI-UPRS method. (a) R = 1 (b) R = 50 (c) R = 200.
Sensors 18 04207 g008aSensors 18 04207 g008b
Figure 9. Theoretical quarter car model.
Figure 9. Theoretical quarter car model.
Sensors 18 04207 g009
Figure 10. (a) Axle mass acceleration and (b) spectrum of the axle mass acceleration.
Figure 10. (a) Axle mass acceleration and (b) spectrum of the axle mass acceleration.
Sensors 18 04207 g010
Figure 11. SNR-surface plots for the vehicle acceleration (a) [4 Hz–25 Hz]. (b) [15 Hz–25 Hz].
Figure 11. SNR-surface plots for the vehicle acceleration (a) [4 Hz–25 Hz]. (b) [15 Hz–25 Hz].
Sensors 18 04207 g011
Figure 12. Extracted Signals PSD (a) first peak at f = 4.99 Hz, (b) second peak at f = 7.12 Hz, (c) third peak at f = 23.41 Hz.
Figure 12. Extracted Signals PSD (a) first peak at f = 4.99 Hz, (b) second peak at f = 7.12 Hz, (c) third peak at f = 23.41 Hz.
Sensors 18 04207 g012
Figure 13. SNR-Surface plot for road profile Class “A” with (a) 2.5-cm increment (b) 25-cm increment.
Figure 13. SNR-Surface plot for road profile Class “A” with (a) 2.5-cm increment (b) 25-cm increment.
Sensors 18 04207 g013
Figure 14. PSD plots for the contact point displacement history of the quarter car: (a) roughness is defined using a 2.5-cm increment (b) roughness is defined using a 25-cm increment.
Figure 14. PSD plots for the contact point displacement history of the quarter car: (a) roughness is defined using a 2.5-cm increment (b) roughness is defined using a 25-cm increment.
Sensors 18 04207 g014
Figure 15. Field instrumentation, (a) bridge overview, (b) bridge instrumentation, (c) truck configuration and axle weights, (d) truck instrumentation, and (e) laser emitter layout.
Figure 15. Field instrumentation, (a) bridge overview, (b) bridge instrumentation, (c) truck configuration and axle weights, (d) truck instrumentation, and (e) laser emitter layout.
Sensors 18 04207 g015
Figure 16. Acceleration PSDs for the Bridge Vibration Tests. (a) Test 1 and (b) Test 2.
Figure 16. Acceleration PSDs for the Bridge Vibration Tests. (a) Test 1 and (b) Test 2.
Sensors 18 04207 g016
Figure 17. Acceleration spectra for the vehicle sensor.
Figure 17. Acceleration spectra for the vehicle sensor.
Sensors 18 04207 g017
Figure 18. Test 3 results (a) FI-UPSR surface plot. (b) PSD for the extracted signal at fapproximate = 7.48 Hz.
Figure 18. Test 3 results (a) FI-UPSR surface plot. (b) PSD for the extracted signal at fapproximate = 7.48 Hz.
Sensors 18 04207 g018
Figure 19. FI-UPSR SNR surface plots and signals PSD for the three tests. (a) Test 1, (b) Test 2, and (c) Test 3.
Figure 19. FI-UPSR SNR surface plots and signals PSD for the three tests. (a) Test 1, (b) Test 2, and (c) Test 3.
Sensors 18 04207 g019
Figure 20. The example illustrates the effect of poor selection of potential parameters on PSD for the extracted signal. (a) FI-UPSR surface plot for Test #3 and (b) PSD plot for extracted signal using the potential parameters for Point #3.
Figure 20. The example illustrates the effect of poor selection of potential parameters on PSD for the extracted signal. (a) FI-UPSR surface plot for Test #3 and (b) PSD plot for extracted signal using the potential parameters for Point #3.
Sensors 18 04207 g020
Table 1. Recommended initial values for Vd, R, and γ.
Table 1. Recommended initial values for Vd, R, and γ.
VdRγ
≥10025–750.1–0.7
For dt = 0.01–0.0005
Table 2. Properties of the quarter car model.
Table 2. Properties of the quarter car model.
PropertyUnitSymbolValue
Body Masskgmb16,600
Axle Masskgms700
Body StiffnessN/mkb2 × 104
Body DampingN.s/mcb10 × 103
Suspension StiffnessN/mks2.75 × 105
Body Bounce FrequencyHzfbounce0.169
Axle Hop FrequencyHzfaxle3.27
Table 3. Properties of the bridge.
Table 3. Properties of the bridge.
PropertyUnitValue
Lengthm15
Mass Per Unit Lengthkg/m28,125
Elastic ModulusMPa35,000
Second Moment of Aream40.5273
1st FrequencyHz5.67
2nd FrequencyHz22.69
3rd FrequencyHz51.05
Table 4. Average bridge frequency (Hz).
Table 4. Average bridge frequency (Hz).
1st2nd3rd4th
7.568.411.9319.04
Table 5. Vehicle speeds km/h.
Table 5. Vehicle speeds km/h.
Test 1Test 2Test 3
375153

Share and Cite

MDPI and ACS Style

Elhattab, A.; Uddin, N.; OBrien, E. Drive-By Bridge Frequency Identification under Operational Roadway Speeds Employing Frequency Independent Underdamped Pinning Stochastic Resonance (FI-UPSR). Sensors 2018, 18, 4207. https://doi.org/10.3390/s18124207

AMA Style

Elhattab A, Uddin N, OBrien E. Drive-By Bridge Frequency Identification under Operational Roadway Speeds Employing Frequency Independent Underdamped Pinning Stochastic Resonance (FI-UPSR). Sensors. 2018; 18(12):4207. https://doi.org/10.3390/s18124207

Chicago/Turabian Style

Elhattab, Ahmed, Nasim Uddin, and Eugene OBrien. 2018. "Drive-By Bridge Frequency Identification under Operational Roadway Speeds Employing Frequency Independent Underdamped Pinning Stochastic Resonance (FI-UPSR)" Sensors 18, no. 12: 4207. https://doi.org/10.3390/s18124207

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