Bridge Health Monitoring via Displacement Reconstruction-Based NB-IoT Technology

: An aged bridge’s performance is a ﬀ ected by degradation and becomes one of the major concerns in maintenance. A preliminary, simple and workable procedure of bridge damage detection is required to minimize maintenance costs. In the past, frequency is one of the most common indicators to detect damage occurrence. Recent research found that using frequency as a health indicator still has room to improve. Alternatively, dynamic displacement is used as an indicator in the current study. These dynamic displacements are reconstructed based on measured acceleration records from micro electro mechanical system (MEMS) sensors. The Newmark-beta method with Windows is proposed to acquire the reconstructed displacements of considered bridges. To demonstrate the accuracy and applicability of the proposed approach, three di ﬀ erent experiments are carried out; (i) A small scale bridge with the implementation of MEMS acceleration sensors; (ii) a numerical complex ﬁnite element method (FEM) bridge model; (iii) an actual bridge with the implementation of MEMS acceleration sensors and narrow bandwidth Internet of things (NB-IoT) technology. The ﬁrst experiment shows that the proposed method can successfully identify the di ﬀ erence between damaged / undamaged bridges and determine damage location. The second experiment indicates that the proposed method is able to identify the di ﬀ erence between sti ﬀ ened / unsti ﬀ ened bridges. The last experiment shows the applicability of the proposed method on an actual bridge health monitoring project.


Introduction
Over a decade, many researchers have attempted to address the issue of bridge health monitoring, in which one of the main focuses is to deliver a simple and sustainable monitoring procedure. Often, the monitoring procedure is required not only accurately to detect the damaged status but also to minimize the cost, including no interruption in transportation. Rytter [1] classify four damage identification levels, in which Level 1 detects damage occurrence; Level 2 localizes bridge damage; Level 3 quantifies the severity of the damage, and Level 4 predicts the remaining service life of a bridge structure. Various approaches are developed to accomplish their own goals of either only Level 1 or up to Level 4 identification. Most approaches could be classified into two groups: model update-based and featured-based approaches.
Model update-based methods require users to construct a numerical structural model either using simplified multiple degrees of freedom (MDOF) or a complex finite element method (FEM) model. The numerical model is built using the initial undamaged bridge property, and its dynamic vehicles and transforms these records into dynamic displacements as the feature for our study. Based on Lee et al. [22] and the Newmark-beta method, a new formulation of the reconstructed displacement algorithm is proposed here. The proposed algorithm can deliver more accurate results without sacrificing cost, and details are introduced in the section of Methodology. Continuous records of reconstructed displacements are used to diagnose the bridge conditions. If any significant difference is observed, damage may occur. Since such a process is repeatedly performed at different locations, a platform or network should be built to fulfill the needs: a long record at several places. In light of this, the narrow bandwidth Internet of things (NB-IoT) technology is adopted and implemented with MEMS sensors in our study. With NB-IoT implementation, the recorded data could be transported directly to servers that are usually located at offices and not connected to the sensors. Table 1 provides a summary of the aforementioned approaches. Users can select a suitable method based on their own needs. Detailed comparison of accuracy, noise resistance or sensitivity to different loading conditions is worth investigation; however, it is beyond the scope of the current study. Table 1. Summary of methods in the literature.

Methods Descriptions
Time-series damage index [13,14] Damage index can be extracted from measured accelerations for two different states under identical loading conditions.
Static influence line deflection [15][16][17] Dynamic properties are not needed. The influencing factors are loading, stiffness, and mass of bridges.
Acceleration sensor attached on passing vehicle [20] Six damage features are extracted from measured accelerations on the passing vehicle data without using equipment on bridges.
Strain fiber optic sensor [19] Damage could be distinguished from the increasing recorded strain embedded in concrete.
Tilt-angle sensor [12] Only accelerations in the gravity direction are measured, and damage status can be distinguished from the constructed rotational influence line.
Global positioning system [20] Measurement without any contact with bridges.
The proposed method Reference-free displacements are constructed from measured accelerations for Two different states without interrupting bridge traffic.
To demonstrate the accuracy and applicability of the proposed method, three different experiments are carried out. (i) Small scale bridge testing with passing vehicles under two different bridge conditions (original and damaged) is considered. In this first experiment, the accuracy of the proposed reconstruction algorithm from accelerations will be validated by comparing the calculated displacements with those from linear variable differential transformers (LVDT). The proposed Level 1 damage detection is also performed and compared to that of using frequency indicator. (ii) A numerical model built by FEM with various artificial damage locations is considered. The proposed approach of Level 2 damage detection and the effect of the number of sensors on a bridge damage detection are examined. (iii) Full-scale bridge testing with the implementation of MEMS accelerometers with NB-IoT technology under daily transportation is considered. The third experiment is intended to demonstrate the applicability of the proposed method on a bridge.

Methodology
Three subsections are included: Section 2.1 derives the formula for displacements reconstructed from accelerations. The idea behind the formula is also provided. Section 2.2 provides MEMS and NB-IoT information that is needed in our approach. Section 2.3 describes the proposed damage identification method of Levels 1 and Level 2. Basically, the proposed method is not able to distinguish different types of bridge damages. The proposed Levels 1 and 2 indicators only reveal the information of vertical deflections for a damaged bridge due to varied causes, such as foundation settlement or shear failure. To be specific, the Level 1 indicator detects damage occurrence and the Level 2 indicator localizes bridge damage. However, the performance of both indicators depends on the number of sensors mounted in the bridge. The more sensors, the more detailed the damaged location.

Reconstructed Dynamic Displacement Algorithm
The objective of the reconstructed dynamic displacement algorithm is to translate measured accelerations into dynamic displacements, initially introduced by Lee et al. [23]. The main concept is to minimize the difference between two accelerations: one is from the records, and the other one is transferred from the assumed displacements. Figure 1 illustrates the procedure of the algorithm. As shown in Figure 1, noise-contaminated acceleration data are the input, and the reconstructed displacement is the output. To be specific, at each time of reconstruction, a window length (N w ) of acceleration are chosen, and only a single displacement is reconstructed. The window slides from left to right with the same increment of the acceleration data. Taking window i = 1 as an example, as shown in Figure 1, acceleration data inside window i = 1 are the input for the reconstruction algorithm, and the output is a single displacement termed as d i=1 . Such a process is repeatedly performed until window i = j. The length of the window (N w ) is determined by an optimization process; details are explained later.

Methodology
Three subsections are included: Section 2.1 derives the formula for displacements reconstructed from accelerations. The idea behind the formula is also provided. Section 2.2 provides MEMS and NB-IoT information that is needed in our approach. Section 2.3 describes the proposed damage identification method of Levels 1 and Level 2. Basically, the proposed method is not able to distinguish different types of bridge damages. The proposed Levels 1 and 2 indicators only reveal the information of vertical deflections for a damaged bridge due to varied causes, such as foundation settlement or shear failure. To be specific, the Level 1 indicator detects damage occurrence and the Level 2 indicator localizes bridge damage. However, the performance of both indicators depends on the number of sensors mounted in the bridge. The more sensors, the more detailed the damaged location.

Reconstructed Dynamic Displacement Algorithm
The objective of the reconstructed dynamic displacement algorithm is to translate measured accelerations into dynamic displacements, initially introduced by Lee et al. [23]. The main concept is to minimize the difference between two accelerations: one is from the records, and the other one is transferred from the assumed displacements. Figure 1 illustrates the procedure of the algorithm. As shown in Figure 1, noise-contaminated acceleration data are the input, and the reconstructed displacement is the output. To be specific, at each time of reconstruction, a window length (Nw) of acceleration are chosen, and only a single displacement is reconstructed. The window slides from left to right with the same increment of the acceleration data. Taking window i = 1 as an example, as shown in Figure 1, acceleration data inside window i = 1 are the input for the reconstruction algorithm, and the output is a single displacement termed as di = 1. Such a process is repeatedly performed until window i = j. The length of the window (Nw) is determined by an optimization process; details are explained later. Lee et al. [22] adopted the central difference method to transfer the assumed displacement to acceleration, in which the reconstructed displacement does not contain noise. Instead of using the Lee et al. [22] adopted the central difference method to transfer the assumed displacement to acceleration, in which the reconstructed displacement does not contain noise. Instead of using the central difference method, the Newmark-beta method is utilized in this study to explore the possibility of enhancing the accuracy in displacement reconstruction. Equations (1) and (2) show the Newmark-beta equations, where x i , .
x i , and ..
x i are the displacement, velocity, and acceleration, Appl. Sci. 2020, 10, 8878 5 of 26 respectively. ∆t refers to the time interval, γ and β are the Newmark constants. In this study, γ and β are taken as 1/2 and 1/6, respectively. Such assumption represents a linear acceleration process: .. ..
Equation (1) could be rewritten in a matrix form, as indicated in Equation (3). As shown in Equation (3), matrix a is a set of acceleration data, and in the proposed algorithm, matrix a is the measured acceleration while L d and L 2 are the constant. The dimension of the matrix a is adjusted to be the same length as each window (from i = 1 to N w ). Thus, the size of matrix L d and matrix L 2 are N w × (N w + 1) and N w × N w , respectively. Matrix v in Equation (3) can be considered as the reconstructed velocity. However, it is understood that the exact solution for v is not possible. where: ..

T
To find the reconstructed velocity, Equation (3) could be transformed into a minimization problem as described in Equation (4). Setting the first derivation with respect to v equals to zero, leading to a new equation, as indicated in Equation (5).
Similarly, Equation (3) can be translated to a matrix form, as shown in Equation (6), in which the size of L 3 is N w × N w . where: where L 3 is a constant and u is defined as the reconstructed displacements, details of obtaining u are provided below. Again, Equation (6) can be considered as a minimization problem as well, in which the objective is formulated, as shown in Equation (7). Equation (7) is differentiated with respect to u to have Equation (2). Using v in Equation (5) and combining it with Equation (8), one can have Equation (9), which is the formula of reconstructed displacement with measured acceleration (a) data as input.
Min(u) = L d u − (v∆t + L 3 a) 2 (7) The Thikonov regularization is often used to tackle the ill condition of the matrix inverse problem encountered from the reconstruction process. In this study, instead of using Thikonov regularization, Moore-Penrose pseudo inverse is utilized to handle the inverse of matrix L d . Please note that although Equation (9) yields a set of displacements, only the middle point will be taken as the reconstructed displacement, then the window will continue to slide to obtain another reconstructed displacement. For a given window, the constructed displacement at the middle is the most accurate one among others. To ensure accuracy, a window-sliding technique is utilized wherein each window, only the outcome at the mid-point data is used. That is, for each window operation, a single point is obtained as the reconstructed displacement.
Besides acceleration data, another parameter needed for Equation (9) is the window size (N w ). N w is affected by the frequency and interval time of acceleration records. To determine the size of N w , this study utilizes a metaheuristic optimization technique to determine the appropriate N w value for each frequency and interval. A single degree of freedom (SDOF) structure with various frequencies from 1-10 Hz with an increment of 1 Hz is generated, in which the accelerations are recorded in different time intervals such as 0.005 s, 0.01 s, and 0.02 s. As a result, 30 different cases are generated, and each SDOF structure is excited by the El Centro earthquake, as shown in Figure 2. Figure 3 shows the ground accelerations of El Centro. Since the structure and excited vibration are given, the exact displacement and acceleration of the structure can be calculated using the numerical integration method. The calculated structural accelerations are then added with noise and transformed to displacements using Equation (9). The added noise is an arbitrary random wave with frequencies between 1 to 200 Hz. In addition, the maximum amplitude of the added noise is twice the absolute mean of structural accelerations.  Metaheuristic optimization is performed to find the most suitable Nw value for each case. As illustrated in Figure 2, the objective of the optimization is to minimize root-mean-square error (RMSE) between exact and reconstructed displacements. Symbiotic organism search (SOS) [23] is used in this study. Although metaheuristic optimization probably needs more computational time compared to that of the conventional gradient-based optimization, metaheuristic optimization could avoid to trap at a local minimum point. Table 2 summarizes the metaheuristic optimization used in this study, including the values of parameters adopted. Figure 4 shows the optimization results of using the SOS algorithm. Based on the power regression lines, the formula for the Nw value could be formulated as shown in Equation (10), where f is the frequency of acceleration and Δt is interval time. Metaheuristic optimization is performed to find the most suitable N w value for each case. As illustrated in Figure 2, the objective of the optimization is to minimize root-mean-square error (RMSE) between exact and reconstructed displacements. Symbiotic organism search (SOS) [23] is used in this study. Although metaheuristic optimization probably needs more computational time Appl. Sci. 2020, 10, 8878 7 of 26 compared to that of the conventional gradient-based optimization, metaheuristic optimization could avoid to trap at a local minimum point. Table 2 summarizes the metaheuristic optimization used in this study, including the values of parameters adopted. Figure 4 shows the optimization results of using the SOS algorithm. Based on the power regression lines, the formula for the N w value could be formulated as shown in Equation (10), where f is the frequency of acceleration and ∆t is interval time. Metaheuristic optimization is performed to find the most suitable Nw value for each case. As illustrated in Figure 2, the objective of the optimization is to minimize root-mean-square error (RMSE) between exact and reconstructed displacements. Symbiotic organism search (SOS) [23] is used in this study. Although metaheuristic optimization probably needs more computational time compared to that of the conventional gradient-based optimization, metaheuristic optimization could avoid to trap at a local minimum point. Table 2 summarizes the metaheuristic optimization used in this study, including the values of parameters adopted. Figure 4 shows the optimization results of using the SOS algorithm. Based on the power regression lines, the formula for the Nw value could be formulated as shown in Equation (10), where f is the frequency of acceleration and Δt is interval time.
Please note that the value of f in Equation (10) only has a moderate influence on the reconstructed displacement. Therefore, an approximate value of the fundamental frequency is often enough. To be specific, it should be acceptable to use a value of f with a 10% error. Using the generated 30 cases, the proposed problem is compared to the previous work by H.S. Lee et al. Table 3 shows the required Nw values for the proposed and the literature methods, respectively. It is seen that the proposed method requires less Nw that lessens the computational cost, and reduces the effect of time lag. In addition to the cost, Table 4 shows that the proposed method also delivers displacements with higher accuracy. It is seen that the performance of the current algorithm is able to reduce the required window size and enhance the accuracy of constructed displacement compared to literature. The current algorithm incorporates the Newmark-Beta as the main component to construct the displacement and utilizes the Moore-Penrose pseudo inverse to tackle the non-symmetrical inverse problem. While the window size is optimized using symbiotic optimization search (SOS). However, it is difficult to Please note that the value of f in Equation (10) only has a moderate influence on the reconstructed displacement. Therefore, an approximate value of the fundamental frequency is often enough. To be specific, it should be acceptable to use a value of f with a 10% error. Using the generated 30 cases, the proposed problem is compared to the previous work by H.S. Lee et al. Table 3 shows the required N w values for the proposed and the literature methods, respectively. It is seen that the proposed method requires less N w that lessens the computational cost, and reduces the effect of time lag. In addition to the cost, Table 4 shows that the proposed method also delivers displacements with higher accuracy. It is seen that the performance of the current algorithm is able to reduce the required window size Appl. Sci. 2020, 10, 8878 8 of 26 and enhance the accuracy of constructed displacement compared to literature. The current algorithm incorporates the Newmark-Beta as the main component to construct the displacement and utilizes the Moore-Penrose pseudo inverse to tackle the non-symmetrical inverse problem. While the window size is optimized using symbiotic optimization search (SOS). However, it is difficult to distinguish the contribution of each component to the increase in performance. Table 2. Summary of metaheuristic optimization used in this study.

MEMS Accelerometer and NB-IoT Technology
As described earlier, to have continuous or stable records in the period of the bridge health monitoring, NB-IoT technology is utilized to transmit the acceleration signals from the bridge site to the server. In this study, the ADXL355 acceleration sensor is used, as shown in Figure 5a. ADXL is selected due to its affordable cost and its adequate ability to sense the motion of the bridge. ADXL355 requires power around 3.6 V and has a sensing range of ±3 g that is suitable for this study. The sampling frequency is 125 Hz, and the working temperature is from −40 • C up to 125 • C, which is suitable for bridge monitoring tasks. The acceleration sensor ADXL355 is mounted to Nordic Semiconductor nRF9160, as shown in Figure 5b that will read the input signal from the sensor.
As described earlier, to have continuous or stable records in the period of the bridge health monitoring, NB-IoT technology is utilized to transmit the acceleration signals from the bridge site to the server. In this study, the ADXL355 acceleration sensor is used, as shown in Figure 5a. ADXL is selected due to its affordable cost and its adequate ability to sense the motion of the bridge. ADXL355 requires power around 3.6 V and has a sensing range of ±3 g that is suitable for this study. The sampling frequency is 125 Hz, and the working temperature is from −40 °C up to 125 °C, which is suitable for bridge monitoring tasks. The acceleration sensor ADXL355 is mounted to Nordic Semiconductor nRF9160, as shown in Figure 5b that will read the input signal from the sensor. nRF9160 is equipped with long term evolution (LTE) with bandwidth 700-960 MHz and 1710-2200 MHz that will send the signal to base transceiver station (BTS) and to the off-site server as shown in Figure 6. The low-cost sensor and NB-IoT package with low power requirements enable the bridge monitoring to be performed at an acceptable price at length. nRF9160 is equipped with long term evolution (LTE) with bandwidth 700-960 MHz and 1710-2200 MHz that will send the signal to base transceiver station (BTS) and to the off-site server as shown in Figure 6. The low-cost sensor and NB-IoT package with low power requirements enable the bridge monitoring to be performed at an acceptable price at length.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 10 of 28 ADXL355 + nRF9160 Base transceiver station Server . Figure 6. General scheme of signal transmission from an on-site sensor to an off-site server.

Proposed Level 1 and Level 2 Damage Detection
This section explains the proposed method that intends to detect damage occurrence and location. The proposed method utilizes the displacement reconstruction algorithm explained in an earlier section. In the bridge monitoring task, accelerometers often spread along a bridge span, as illustrated in Figure 7. Assuming that i number of sensor points (SP) spread along a bridge, as shown in Figure 7 and each sensor records acceleration data. If there is damage occurred, the collected accelerations naturally include two sets of records. Namely, the original and observed sets. These two sets of accelerations are then transformed into reconstructed displacements (RDisp) using Equation (10). The root-mean-square (RMS) of RDisp for each sensor point (RMSDi) is calculated as formulated in Equation (11). where n is the quantity of data.
where n is the quantity of data. Figure 6. General scheme of signal transmission from an on-site sensor to an off-site server.

Proposed Level 1 and Level 2 Damage Detection
This section explains the proposed method that intends to detect damage occurrence and location. The proposed method utilizes the displacement reconstruction algorithm explained in an earlier section. In the bridge monitoring task, accelerometers often spread along a bridge span, as illustrated in Figure 7. Assuming that i number of sensor points (SP) spread along a bridge, as shown in Figure 7 and each sensor records acceleration data. If there is damage occurred, the collected accelerations naturally include two sets of records. Namely, the original and observed sets. These two sets of accelerations are then transformed into reconstructed displacements (RDisp) using Equation (10). The root-mean-square (RMS) of RDisp for each sensor point (RMSD i ) is calculated as formulated in Equation (11).
where n is the quantity of data.
where n is the quantity of data. Figure 7. Illustration of sensor points and damage index 1.
Each original and observed sets consist of i number of RMSD, and for either set and is normalized with its own maximum value if shown by a diagram with respect to the sensor position will create a deflection shape as shown in Figure 7. Level 1 damage indicator is the shape difference between the original and observed RMSDs, as shown in Figure 8. Once there is a significant difference between these two RMSDs, the bridge will be diagnosed as a damaged one. The larger the damage indicator, the more severely damaged the bridge. A threshold could be suggested based on available Each original and observed sets consist of i number of RMSD, and for either set and is normalized with its own maximum value if shown by a diagram with respect to the sensor position will create a deflection shape as shown in Figure 7. Level 1 damage indicator is the shape difference between the original and observed RMSDs, as shown in Figure 8. Once there is a significant difference between these two RMSDs, the bridge will be diagnosed as a damaged one. The larger the damage indicator, the more severely damaged the bridge. A threshold could be suggested based on available data. As shown in Figure 8, the difference between two RMSDs at the SP 4 position is zero since two sets of data are normalized with respect to the same midpoint.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 11 of 28 data. As shown in Figure 8, the difference between two RMSDs at the SP4 position is zero since two sets of data are normalized with respect to the same midpoint. For Level 2 damage indicator, the dRMSD, which is the difference between adjacent RMSDs as formulated in Equation (12), is calculated ( . Same as Level 1 damage indicator, original and observed dRMSDs are computed as illustrated in Figure 9. Subsequently, the difference between original and observed dRMSDs are plotted with respect to the sensor point, as shown in Figure 9. Each segment of the bridge is represented by a single value of damage indicator (i.e., dRMSD) and If one of the values is particularly prominent, indicating that a particular segment is in a damaged status. The more sensors, the more accurate in locating the damaged position.  For Level 2 damage indicator, the dRMSD, which is the difference between adjacent RMSDs as formulated in Equation (12), is calculated ( RMSD i is Zero ). Same as Level 1 damage indicator, original and observed dRMSDs are computed as illustrated in Figure 9. Subsequently, the difference between original and observed dRMSDs are plotted with respect to the sensor point, as shown in Figure 9. Each segment of the bridge is represented by a single value of damage indicator (i.e., dRMSD) and If one of the values is particularly prominent, indicating that a particular segment is in a damaged status. The more sensors, the more accurate in locating the damaged position.
dRMSD) and If one of the values is particularly prominent, indicating that a particular segment is in a damaged status. The more sensors, the more accurate in locating the damaged position.

Bridge Health Monitoring Experiment
To validate the aforementioned displacement reconstruction method and damage indicators, three different experiments are introduced in this section.

Bridge Health Monitoring Experiment
To validate the aforementioned displacement reconstruction method and damage indicators, three different experiments are introduced in this section.

Experiment Setup
The experiment setup is shown in Figure 10, where the bridge has a 5.6 m span. During the testing, the excitation of the bridge is induced by a moving cart on the bridge span, as indicated in Figure 10. MEMS accelerometer with a sampling rate of 125 Hz and LVDT with a sampling rate of 200 Hz is equipped at the first quarter, middle, and the third quarter of the bridge as shown in Figure 10.
For bridge frequency, the span length of small scale bridge is set so that its first three mode frequencies are 3.7 Hz, 15 Hz, and 30 Hz, respectively. As shown, these three frequencies are reasonable compared to those of a real bridge. For the mass ratio, two different types of ratios are used: light and heavy. The bridge weight is 291 kg, and the weights of vehicles are 30 kg and 47 kg, respectively. The mass ratio for the light case is around 10%, and the mass ratio for the heavy case is 16%. The heavy case is intended to model two consecutive trucks or a trailer truck with full loading. For the length ratio, the light case is 7%, and the heavy case is 16%. For the speed ratio (α), the similarity ratio adopted is described in Equation (13) To simulate various conditions of bridge and passing vehicle, 12 different scenarios are designed as indicated in Table 5. Two different bridge conditions are used, namely the original and the stiffened ones. The stiffened bridge refers to a stiffener is attached to one of the bridge spans to increase its stiffness, as shown in Figure 11b. Figure 10a indicates the location of the stiffened span. Two vehicles with different weights are used, namely heavy and light carts. The weights are 30.4 kgf and 47.3 kgf for the light and heavy carts, respectively. Three different speeds are used, and they are 0.4, 0.8, and 1.2 m/s, respectively. In total, 12 different scenarios are considered as listed in Table 5. For each scenario, the test will be repeated 10 times; each time, the cart will run back and forth once, as shown in Figure 10a.
Two vehicles with different weights are used, namely heavy and light carts. The weights are 30.4 kgf and 47.3 kgf for the light and heavy carts, respectively. Three different speeds are used, and they are 0.4, 0.8, and 1.2 m/s, respectively. In total, 12 different scenarios are considered as listed in Table 5. For each scenario, the test will be repeated 10 times; each time, the cart will run back and forth once, as shown in Figure 10a.      Table 6 shows the dynamic bridge properties, such as the first three frequencies for the original and stiffened bridges. It is seen that although this is a small scale bridge, its dynamic property is adjusted so that it is similar to an actual bridge condition. In Table 6, frequency is measured under free vibration 5 times. It is shown that the stiffened bridge only slightly increases its frequency from 3.733 to 3.867 Hz. Table 6 also displays the coefficient of variation (cov) among 5 testing data. It is shown that covs are very small, indicating the measurement is very accurate. Table 7 shows the  Table 6 shows the dynamic bridge properties, such as the first three frequencies for the original and stiffened bridges. It is seen that although this is a small scale bridge, its dynamic property is adjusted so that it is similar to an actual bridge condition. In Table 6, frequency is measured under free vibration 5 times. It is shown that the stiffened bridge only slightly increases its frequency from 3.733 to 3.867 Hz. Table 6 also displays the coefficient of variation (cov) among 5 testing data. It is shown that covs are very small, indicating the measurement is very accurate. Table 7 shows the damping ratio for the original and stiffened bridges under free vibration. Similar to frequency measurement, the number of measurements is 5.

Level 1 Detection Using Frequency
From 10 experiments of each scenario, statistics of the fundamental frequency are shown in Table 8. The fundamental frequency is calculated using the direct fast Fourier transform (FFT). Mean, Maximum, and Minimum among 10 experiments are presented. It is shown that frequencies are changed accordingly due to varied vehicle speeds and masses. The increasing vehicle can reduce the fundamental frequency while increasing vehicle mass tends to decrease fundamental frequency. In addition, the position of the vehicle also could affect the mass matrix of the entire bridge, resulting in a different frequency. For the current state, this study only reports the phenomena found in our small scale experiment; no further statement is concluded for other small scale or real bridge conditions. Instead of using direct FFT, FFT with window is performed. Results of windowed FFT are normalized to 1. Figures 12 and 13 show the displacements and FFT spectrogram for scenario LVB0S1 and LVB0S3, respectively. Figures 12 and 13 contain three plots. The first plot shows the dynamic displacement at mid-span under a moving load where the static displacement has been filtered out, and the second plot shows the FFT spectrogram, where each window contains 1000 data and 999 of them are overlapped with the next windows. The third plot presents two FFTs at different times indicated in the second plot. Figures 12 and 13 show that different cart positions could yield different fundamental frequencies. Using frequency as a bridge health indicator with passing vehicles should be cautious.

Level 1 Detection Using Reconstructed Displacements
This subsection introduces the results of the Level 1 detection using reconstructed displacements explained in Section 2.1. As described, both MEMS accelerometers and LVDT are used for measurement, as shown in Figure 10a; comparisons between two devices are conducted to validate the accuracy of the proposed approach. Figures 14 and 15 (top) shows the measured accelerations using MEMS for scenarios of LVB0S1 and LVB0S3, respectively. Figures 14 and 15 (bottom) shows the displacements from LVDT and reconstructed displacement using MEMS for scenarios of LVB0S1 and LVB0S3, respectively. It is shown that reconstructed displacement is almost identical to those of LVDT. Note that noise has been filtered out in the reconstructed displacements.

Level 1 Detection Using Reconstructed Displacements
This subsection introduces the results of the Level 1 detection using reconstructed displacements explained in Section 2.1. As described, both MEMS accelerometers and LVDT are used for measurement, as shown in Figure 10a; comparisons between two devices are conducted to validate the accuracy of the proposed approach. Figures 14 and 15 (top) shows the measured accelerations using MEMS for scenarios of LVB0S1 and LVB0S3, respectively. Figures 14 and 15 (bottom) shows the displacements from LVDT and reconstructed displacement using MEMS for scenarios of LVB0S1 and LVB0S3, respectively. It is shown that reconstructed displacement is almost identical to those of LVDT. Note that noise has been filtered out in the reconstructed displacements.  As shown in Figure 10, three accelerometers are attached to the bridge. To demonstrate the proposed approach, the mean value of RMSDs between original and stiffened bridges in the H1 direction is compared with those of the H2 direction. Figure 16 shows the mean RMSDs in H1 and H2 directions for (a) original (Cases 1 to 6) and (b) stiffened bridge (Cases 7 to 12). Points No. 1 and No. 5 represent the two ends of the test bridge, while points No. 2, No. 3, and No. 4 represent the sensor location at the first, second, and third quarters of the bridge. Figure 16 (top) clearly shows that the shapes created by the RMSD value have no significant change between H1 and H2 directions since both plots have identical conditions (only the cart moving direction is opposite). Level 1 damage indicator between H1 and H2 directions is less than 5E-3. However, if the Level 1 damage indicator is compared between the original and stiffened bridges, the shapes are quite different from those of  As shown in Figure 10, three accelerometers are attached to the bridge. To demonstrate the proposed approach, the mean value of RMSDs between original and stiffened bridges in the H1 direction is compared with those of the H2 direction. Figure 16 shows the mean RMSDs in H1 and H2 directions for (a) original (Cases 1 to 6) and (b) stiffened bridge (Cases 7 to 12). Points No. 1 and No. 5 represent the two ends of the test bridge, while points No. 2, No. 3, and No. 4 represent the sensor location at the first, second, and third quarters of the bridge. Figure 16 (top) clearly shows that the shapes created by the RMSD value have no significant change between H1 and H2 directions since both plots have identical conditions (only the cart moving direction is opposite). Level 1 damage indicator between H1 and H2 directions is less than 5E-3. However, if the Level 1 damage indicator is compared between the original and stiffened bridges, the shapes are quite different from those of As shown in Figure 10, three accelerometers are attached to the bridge. To demonstrate the proposed approach, the mean value of RMSDs between original and stiffened bridges in the H1 direction is compared with those of the H2 direction. Figure 16 shows the mean RMSDs in H1 and H2 directions for (a) original (Cases 1 to 6) and (b) stiffened bridge (Cases 7 to 12). Points No. 1 and No. 5 represent the two ends of the test bridge, while points No. 2, No. 3, and No. 4 represent the sensor location at the first, second, and third quarters of the bridge. Figure 16 (top) clearly shows that the shapes created by the RMSD value have no significant change between H1 and H2 directions since both plots have identical conditions (only the cart moving direction is opposite). Level 1 damage indicator between H1 and H2 directions is less than 5 × 10 −3 . However, if the Level 1 damage indicator is compared between the original and stiffened bridges, the shapes are quite different from those of Figure 16. As shown in Figure 17, the Level 1 damage indicator is around 1.8 × 10 −2 up to 3.5 × 10 −2 .
The proposed indicator could recognize the occurrence of the stiffened plate while the frequency indicator could not. Need to be noted, the comparison is performed between Cases 1 to 6 and Cases 7 to 12 to reproduce the actual bridge conditions, in which a bridge is often passed by various vehicle types and vehicle speeds.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 18 of 28 Figure 16. As shown in Figure 17, the Level 1 damage indicator is around 1.8E-2 up to 3.5E-2. The proposed indicator could recognize the occurrence of the stiffened plate while the frequency indicator could not. Need to be noted, the comparison is performed between Cases 1 to 6 and Cases 7 to 12 to reproduce the actual bridge conditions, in which a bridge is often passed by various vehicle types and vehicle speeds. Appl. Sci. 2020, 10, x FOR PEER REVIEW 18 of 28 Figure 16. As shown in Figure 17, the Level 1 damage indicator is around 1.8E-2 up to 3.5E-2. The proposed indicator could recognize the occurrence of the stiffened plate while the frequency indicator could not. Need to be noted, the comparison is performed between Cases 1 to 6 and Cases 7 to 12 to reproduce the actual bridge conditions, in which a bridge is often passed by various vehicle types and vehicle speeds.

Analysis of a Numerical Bridge
In this section, a numerical model is created to demonstrate the proposed method in recognizing the damage position for a given bridge, that is, the Level 2 damage detection. Because the accuracy of Level 2 damage detection depends on the number of sensors, a numerical model is more convenient to fulfill this requirement and is therefore adopted here. Only the vertical direction is considered in the numerical model. Both ends are fixed. Mass and stiffness are divided into 200 increments, and 9 sensors are uniformly distributed along the span. Figure 18 shows the bridge dimension and accelerometer positions. As shown, on a 40-m bridge, sensors are arranged every 4 m. The bridge consists of 200 finite elements with assigned mass and stiffness. For the given mass and stiffness, the mass and stiffness matrix could be constructed, and together with the damping matrix, the governing equation is displayed in Equation (13). The bridge property is listed in Table 9.
Appl. Sci. 2020, 10, x FOR PEER REVIEW 19 of 28 Figure 17. Level 1 damage indicator of (a) original and stiffened bridges in H1 direction; (b) original and stiffened bridges in H2 direction.

Analysis of a Numerical Bridge
In this section, a numerical model is created to demonstrate the proposed method in recognizing the damage position for a given bridge, that is, the Level 2 damage detection. Because the accuracy of Level 2 damage detection depends on the number of sensors, a numerical model is more convenient to fulfill this requirement and is therefore adopted here. Only the vertical direction is considered in the numerical model. Both ends are fixed. Mass and stiffness are divided into 200 increments, and 9 sensors are uniformly distributed along the span. Figure 18 shows the bridge dimension and accelerometer positions. As shown, on a 40-m bridge, sensors are arranged every 4 m. The bridge consists of 200 finite elements with assigned mass and stiffness. For the given mass and stiffness, the mass and stiffness matrix could be constructed, and together with the damping matrix, the governing equation is displayed in Equation 13. The bridge property is listed in Table 9.
where M, C, and K are the structural mass, damping, and stiffness matrix, respectively. Dynamic analysis is conducted using the state-space method provided by MATLAB SIMULINK, where only displacement in the gravity direction is calculated. The vibration source in this numerical analysis is ambient vibration ( a x  ), which is simulated by white noise with a frequency ranging from 1 to 10 Hz. Figure 19 shows the first three modes of the bridge, where the frequencies of the first three modes are 3.497, 6.994, and 10.490 Hz, respectively. To demonstrate the proposed Level 2 damage detection method, 6 cases of damage scenario are artificially created as listed in Table 10. In cases 1 to 4, a bridge with single damage is considered. On the other hand, a bridge with two damages is considered in cases 5 and 6. The damage is generated by reducing the stiffness of the original bridge by a factor of 0.9 with one-meter-long, as shown in Table 10. For each damage case, the obtained acceleration data from original and damage bridges are used to calculate the Level 2 damage indicator. Figures 20-25 display the Level 2 damage indicators for each case, indicating that the proposed Level 2 indicator is able to provide a promising prediction for the damage location, either with a single or double damages.  ..
where M, C, and K are the structural mass, damping, and stiffness matrix, respectively. Dynamic analysis is conducted using the state-space method provided by MATLAB SIMULINK, where only displacement in the gravity direction is calculated. The vibration source in this numerical analysis is ambient vibration ( .. x a ), which is simulated by white noise with a frequency ranging from 1 to 10 Hz. Figure 19 shows the first three modes of the bridge, where the frequencies of the first three modes are 3.497, 6.994, and 10.490 Hz, respectively. Table 9. Bridge properties.

Bridge Property Value Unit
Young's modulus 20,000 MPa Moment inertia 3.27 × 10 5 mm 4 Total weight 10 tonf To demonstrate the proposed Level 2 damage detection method, 6 cases of damage scenario are artificially created as listed in Table 10. In cases 1 to 4, a bridge with single damage is considered. On the other hand, a bridge with two damages is considered in cases 5 and 6. The damage is generated by reducing the stiffness of the original bridge by a factor of 0.9 with one-meter-long, as shown in Table 10. For each damage case, the obtained acceleration data from original and damage bridges are used to calculate the Level 2 damage indicator. Figures 20-25 display the Level 2 damage indicators for each case, indicating that the proposed Level 2 indicator is able to provide a promising prediction for the damage location, either with a single or double damages. Appl. Sci. 2020, 10, x FOR PEER REVIEW 20 of 28 Figure 19. The first three mode shapes of a numerical bridge.  Figure 19. The first three mode shapes of a numerical bridge.  Figure 19. The first three mode shapes of a numerical bridge.  Figure 19. The first three mode shapes of a numerical bridge.  Figure 19. The first three mode shapes of a numerical bridge.  Figure 19. The first three mode shapes of a numerical bridge.  Figure 19. The first three mode shapes of a numerical bridge.

Full-scale Bridge Experiment
A full-scale bridge is used to demonstrate the implementation of the proposed Level 1 damage detecting method. Wenzhou Bridge, located in Taipei city, is selected as the test bridge because it is an important bridge, and the riverside park under the bridge provides test personnel to facilitate the assembly of test equipment. Wenzhou bridge is constructed in the Year 2000. It is a prestressed girder bridge with four lanes. No damage is visually detected before our testing. The experiments are initiated at the beginning of April 2020, and field data are continuously collected until now. Due to

Full-scale Bridge Experiment
A full-scale bridge is used to demonstrate the implementation of the proposed Level 1 damage detecting method. Wenzhou Bridge, located in Taipei city, is selected as the test bridge because it is an important bridge, and the riverside park under the bridge provides test personnel to facilitate the assembly of test equipment. Wenzhou bridge is constructed in the Year 2000. It is a prestressed girder bridge with four lanes. No damage is visually detected before our testing. The experiments are initiated at the beginning of April 2020, and field data are continuously collected until now. Due to

Full-Scale Bridge Experiment
A full-scale bridge is used to demonstrate the implementation of the proposed Level 1 damage detecting method. Wenzhou Bridge, located in Taipei city, is selected as the test bridge because it is an important bridge, and the riverside park under the bridge provides test personnel to facilitate the assembly of test equipment. Wenzhou bridge is constructed in the Year 2000. It is a prestressed girder bridge with four lanes. No damage is visually detected before our testing. The experiments are initiated at the beginning of April 2020, and field data are continuously collected until now. Due to lack of information and based on the visual inspection, the status of the bridge in April 2020 is assumed to be in a healthy condition. Wenzhou Bridge is an RC prestress girder bridge with dimensions shown in Figure 26. Accelerometers integrated with NB-IoT are attached along the bridge span at the first, second, and third quarters. The power of the equipment is supported by several detached batteries that need to be refilled regularly. The sampling rate is 125 Hz, and acceleration data are recorded every 12 min. Each record consists of 4096 points. That is, each record is approximately 33 s-long. The collected acceleration data are transformed into displacements based on the proposed algorithm. Because the proposed reconstructed displacement algorithm requires the bridge's fundamental frequency, FFT analysis is first performed to obtain the frequency, as shown in Figure 27. It is seen that a technical problem occurred from 6/15 to 6/30 (month/date), causing a gap in collecting data. However, such a gap does not affect the following analysis. The frequency of the bridge is assumed to be 4 Hz that is calculated based on the mean of the calculated FFTs. As explained in Section 3.1.2, different frequencies are observed, as shown in Figure 27, where the frequency is ranging from 3.5 Hz to 4.5 Hz.
When applying the proposed Level 1 damage detection to the Wenzhou Bridge, the condition in April is assumed to be the original condition, which is considered as a reference point for calculating the Level 1 damage indicator. Damage indicator is calculated for each month except June. Figures 28-30 show the Level 1 damage indicator for the Wenzhou Bridge. It is shown that the bridge still remains in the same condition as that in April. Figures 28-31 depict the proposed Level 1 indicator by comparing the normalized RMSD in April. If the difference is large, indicating that there is damage occurrence. As shown, the largest RMSD difference is 8 × 10 −3, which is considered a relatively small value compared to those of the small-scale experiments. It is, therefore, concluded that the bridge condition from May to September does not significantly change from April.
in collecting data. However, such a gap does not affect the following analysis. The frequency of the bridge is assumed to be 4 Hz that is calculated based on the mean of the calculated FFTs. As explained in Section 3.1.2, different frequencies are observed, as shown in Figure 27, where the frequency is ranging from 3.5 Hz to 4.5 Hz.
When applying the proposed Level 1 damage detection to the Wenzhou Bridge, the condition in April is assumed to be the original condition, which is considered as a reference point for calculating the Level 1 damage indicator. Damage indicator is calculated for each month except June. Figures 28  to 30 show the Level 1 damage indicator for the Wenzhou Bridge. It is shown that the bridge still remains in the same condition as that in April. Figures 28-31 depict the proposed Level 1 indicator by comparing the normalized RMSD in April. If the difference is large, indicating that there is damage occurrence. As shown, the largest RMSD difference is 8 × 10 −3, which is considered a relatively small value compared to those of the small-scale experiments. It is, therefore, concluded that the bridge condition from May to September does not significantly change from April.

Conclusions
This study introduces a novel approach for bridge health monitoring utilizing a displacement reconstruction algorithm. The algorithm is developed utilizing the Newmark-beta method, the Moore-Penrose pseudo inverse and the sliding window technique. Level 1 and Level 2 damage detections are considered in this study. Three experiments are carried out to validate and demonstrate the proposed method. Several brief conclusions and discussions could be drawn as described below: 1. The proposed displacement reconstruction algorithm requires less computational time and has higher accuracy compared to the previous study, as provided in Section 2.1. The input parameters of the proposed algorithm are recorded accelerations, fundamental frequency and window size; 2. Level 1 damage detection is performed by first using FFT to identify the fundamental frequency of a bridge. It is shown that various vehicle weights and velocities could affect the dominant

Conclusions
This study introduces a novel approach for bridge health monitoring utilizing a displacement reconstruction algorithm. The algorithm is developed utilizing the Newmark-beta method, the Moore-Penrose pseudo inverse and the sliding window technique. Level 1 and Level 2 damage detections are considered in this study. Three experiments are carried out to validate and demonstrate the proposed method. Several brief conclusions and discussions could be drawn as described below: 1. The proposed displacement reconstruction algorithm requires less computational time and has higher accuracy compared to the previous study, as provided in Section 2.1. The input parameters of the proposed algorithm are recorded accelerations, fundamental frequency and window size; 2. Level 1 damage detection is performed by first using FFT to identify the fundamental frequency of a bridge. It is shown that various vehicle weights and velocities could affect the dominant

Conclusions
This study introduces a novel approach for bridge health monitoring utilizing a displacement reconstruction algorithm. The algorithm is developed utilizing the Newmark-beta method, the Moore-Penrose pseudo inverse and the sliding window technique. Level 1 and Level 2 damage detections are considered in this study. Three experiments are carried out to validate and demonstrate the proposed method. Several brief conclusions and discussions could be drawn as described below: 1.
The proposed displacement reconstruction algorithm requires less computational time and has higher accuracy compared to the previous study, as provided in Section 2.1. The input parameters of the proposed algorithm are recorded accelerations, fundamental frequency and window size; 2.
Level 1 damage detection is performed by first using FFT to identify the fundamental frequency of a bridge. It is shown that various vehicle weights and velocities could affect the dominant frequency. Window-slide FFT is also performed to examine the variation of the dominant frequency per window, and it shows that dominant frequency varies as different vehicle position and speed; 3. Although the stiffness difference between the original and the stiffened bridge in the small-scale bridge experiment is slight, the proposed Level 1 damage detection using reconstructed displacement is able to distinguish the difference between the two; 4.
The proposed Level 2 damage detection is successfully to predict the damage location for either case: single or double damages for the numerical examples; 5.
To apply the proposed damage detection approach to an actual bridge, MEMS accelerometer and NB-IoT are integrated with the proposed methodology. On-site monitoring was conducted from April until November. Based on the proposed Level 1 damage detection, the monitored bridge shows no indication of degradation since April. The proposed method enables on-site monitoring to be conducted without interrupting traffic activity at an acceptable cost.
Please note that the input of the proposed reconstructed displacement algorithm is the measured accelerations. That is, the proposed method cannot consider the effect of each individual factor, including environmental factors, separately. However, the effects of the environmental factors are included in the measured data, and the corresponding influences are considered. In addition, in a real case, the deformation of a bridge is in 3 directions. Although most failure modes are related to the vertical deflection, some failure modes, such as shear failure, is not, and the proposed method may not be suitable in such case. The proposed method is only suitable for a bridge with failure mode related to the gravity direction. The current approach could only provide an early warning for a bridge system. A mechanical-based model is needed to further clarify the possible failure mechanisms for a bridge.