1. Introduction
Currently, scholars have conducted extensive research on the topic of damage identification in bridge structures both domestically and internationally. The analysis methods for damage identification in bridge structures include the following: (1) identification methods based on dynamic fingerprints; (2) identification methods based on model correction; (3) identification methods based on static characteristics; and (4) identification methods based on time-frequency signals [
1,
2,
3,
4,
5]. Dynamic fingerprints, indicative of structural dynamic characteristics, involve discerning the structural damage state by observing changes in dynamic fingerprints before and after damage occurrence [
6,
7]. Huang et al. [
8] proposed a vibration-based non-destructive global damage identification method using a genetic algorithm. By combining frequency and modal shapes, this method can be applied to identify the location and severity of structural damage under the influence of temperature variations and noise. Model correction involves iteratively adjusting simulated parameters towards responses that closely resemble the actual static–dynamic responses of real structures. Innocenzi et al. [
9] integrated three methods—static analysis, dynamic analysis, and model correction—to conduct static and dynamic tests on a cable-stayed steel–concrete composite bridge. They developed a reliable structural digital twin model, demonstrating that expanded measurements and improved models can effectively explain the behavior of complex structures. This work lays a solid foundation for future bridge structural health monitoring.
The damage identification method based on static characteristics is typically employed when the structure is in a static state. Zhou et al. [
10] considered the uncertainty of cross-sectional parameters and proposed a method to locate damage in simply supported beam bridges using the influence line interpolation of arbitrary cross-sectional rotation angles. However, this method is characterized by a relatively simple model and narrow applicability. Son et al. [
11] utilized static displacement, slope, and curvature to detect damage in bridges. They calculated static displacement through finite element analysis, determined slope and curvature using the central difference method, and thereby assessed the sensitivity of damaged bridges. Wang et al. [
12] introduced a damage identification analysis method based on the influence line interpolation index of mid-span displacement, addressing the issue of uncertain bending stiffness of the main beam and effectively locating the damage position. To obtain comprehensive data and reduce dependence on initial data, Sun et al. [
13] measured the vertical displacement influence line of a bent bridge in a quasi-static state. By obtaining the curvature of the influence line through second-order differencing and employing gap smoothing techniques, they constructed a damage index for the identification of damage in such structures.
The method of identification based on time-frequency signals refers to the analysis of vibration response signals under external excitation when a structure undergoes damage. These signals often exhibit non-stationary and non-linear characteristics. Analyzing the time-frequency characteristics of vibration responses allows for a deeper exploration of hidden features within the signals [
14]. Wu et al. [
15] employed wavelet packet analysis to decompose acceleration data into wavelet packets. The energy values of the decomposed frequency bands (wavelet packet energy (WPE)) were used as different dimensions of MDS. They constructed a damage index combining WPE and MDS values, determining the damage location based on whether the MDS values of units exceeded a threshold. Yen et al. [
16] proposed the concept of nodal energy based on the principles of wavelet packet analysis. Their research indicated that nodal energy exhibited higher robustness compared to decomposition coefficients in representing signal features. Sun et al. [
17] introduced a statistical pattern classification method based on wavelet packet transform (WPT). The core of this method lies in the ability of wavelet packet transform to extract subtle anomalies from vibration signals, enabling structural damage identification. Soleymani et al. [
18] utilized time-domain modal testing and wavelet analysis to identify damage in reinforced concrete beams. They generated various damage scenarios of different severity and locations using a numerical model of an RC beam, recorded acceleration time histories of damaged and undamaged structures, and ultimately determined the location and severity of damage through wavelet analysis.
In previous studies, damage identification of bridge structures often relied on single indicators such as displacement influence lines and wavelet packet energy. However, there has been limited research on the combined analysis of different indicators, and there is a lack of comparative analysis among different damage indicators. In current structural damage identification research, the majority of methods are concentrated on theoretical analysis and finite element numerical simulations. The reliability and applicability of theoretical analysis methods lack corresponding validation through experimental model testing.
Therefore, in response to the aforementioned challenges, this study proposes a bridge structural damage identification method that combines two indicators: Quasi-Static Displacement-Induced Linearity Error Curve (QSDIL) and Relative Energy Rate of Wavelet Packet Energy Spectrum (RES). Firstly, the QSDIL, a damage feature indicator, is constructed to preliminarily locate the position of damage in bridge structures based on a quasi-static displacement influence linearity error. Subsequently, following the principles of wavelet packet decomposition and the denoising theory of wavelet packet transform, the RES, another damage feature indicator, is constructed for accurate damage localization in bridge structures.
Simultaneously, in the laboratory, bridge structural test models are manufactured, and a bridge structural damage identification experimental system is designed to conduct experiments on bridge structural damage identification. Using the finite element software ANSYS 19.2, corresponding bridge models are established, and experimental values are compared and analyzed against theoretical calculations under various test conditions. Furthermore, this study explores the influence of factors such as single-point damage, multi-point damage, signal noise, lane position, and vehicle weight on the identification results. This method requires only a small number of sensors to achieve precise localization of bridge structural damage positions, providing valuable insights and references for experimental research on bridge structural damage identification.
2. Bridge Structural Damage Identification Method Based on Quasi-Static Displacement Influence Line and Wavelet Packet Analysis
2.1. Analysis of Differences in Quasi-Static Displacement Influence Line Theory
The displacement influence line is introduced based on the concept of force influence lines [
19]. It represents the curve of the displacement at an observation point due to the movement of a unit load along the span of the structure. This curve varies as the position of the unit load changes. It is used to determine displacement values under the combined action of multiple loads and to identify the most critical location for moving loads on a structure. The quasi-static displacement influence line is an extension of the displacement influence line. It involves subjecting the structure to quasi-static loading to obtain the quasi-static displacement influence line. This line is then used for damage identification based on the structure’s response to quasi-static loading. In the context of bridge structures, a point’s displacement influence line represents the displacement curve generated at that point as a unit of concentrated force moves along the bridge direction due to load movement.
Consider a unit of concentrated force Fp = 1 with a constant direction and magnitude applied to a simply supported single-span beam. This force is moving at an extremely low velocity from support A towards support B. Let point C denote the location of a localized damage in the beam, where the flexural stiffness of the damaged region within the range [
c−ε,
c+
ε] is EI′, while the remaining undamaged segments of the beam have a flexural stiffness of EI, as illustrated in
Figure 1.
Damage information of a bridge structure is analyzed by the difference in quasi-static displacement influence lines, denoted as Δω(x), before and after damage at an arbitrary section D of a simply supported beam. The expression for the function of Δω(x) with respect to the distance from the moving load to point A can be derived as follows:
- (1)
- (2)
When
, then
- (3)
When
, then
- (4)
When
, then
From the above equation, it is evident that the minimum value of the QSDIL curve must occur within the damaged region. Therefore, by analyzing the location of the peak on the QSDIL curve, it is possible to identify the damaged area, assess whether structural damage has occurred, and thus make an initial determination of the damaged region.
2.2. The Theory of Wavelet Packet Analysis
Wavelet packet analysis naturally possesses denoising capabilities. Its algorithm core inherits the fundamental characteristics of wavelet denoising while enhancing time-frequency resolution. It can simultaneously distinguish between the high- and low-frequency components of a signal at the same frequency [
20]. In digital signal analysis, signal energy refers to the sum of the squared amplitudes of the signal at various points within a certain time frame. Typically, a signal
x(
t) is decomposed into high-frequency and low-frequency components at different scales, and then the energy of the signal within each wavelet packet frequency band is summed, thereby obtaining the total energy of the entire signal.
The process of calculating the energy of a signal using wavelet packet analysis is as follows.
The signal
x(
t) collected by the sensor is decomposed using wavelet packet analysis, resulting in the following:
cji,r represents the wavelet packet coefficient, where
i is the scale index,
r is the position index, and
j represents the frequency index.
ψi,r,,j(
t) represents the wavelet packet function, which is a set of standard orthogonal bases.
By performing a wavelet packet decomposition on the original signal, the obtained wavelet coefficients can be utilized for calculating the energy components of the signal at various frequency bands. Within each wavelet packet frequency band, the computation of the sum of squared wavelet coefficients yields the energy component for that specific frequency band. Summing up the energy components from all frequency bands results in the total wavelet packet energy. At this stage, the total energy of the signal is given by the following:
From the orthogonality of wavelet packets, it follows that
In the equation, the energy of wavelet packet components is denoted as Et, and can be considered as the energy contained within the component signal. Ei represents the energy of the i-th frequency band within the structure, while Et represents the total energy of the wavelet packet decomposition. Here, i corresponds to the number of frequency bands and j corresponds to the number of wavelet packet coefficients within each frequency band.
RES is a damage characteristic indicator based on the total energy change rate of wavelet packet coefficients. It is expressed as follows:
RES represents the relative energy change rate of the wavelet packet, E(t,undam) represents the total energy of undamaged wavelet packets in the structure, and E(t,dam) represents the total energy of wavelet packets in the structure with damage.
In the presence of damage in a bridge structure, the structural response signals display oscillations within a particular frequency range. Within this frequency range, specific segments will experience alterations, and the energy of signal components within particular frequency bands, as determined through wavelet packet analysis, will likewise change. Consequently, by comparing the energy of signal components within specific frequency bands before and after structural damage, it becomes feasible to characterize the intrinsic properties of the structure and make assessments regarding its structural integrity.
2.3. Process for Structural Damage Identification
The damage identification procedure proposed in this paper is grounded in a method for damage identification and analysis, which includes an initial damage assessment and precise damage localization. It involves the analysis and processing of response signals from bridge structures to facilitate the completion of identification and analysis. This method demands only a limited number of sensors for accurately localizing damage within the bridge structure, as illustrated in
Figure 2. The specific steps are outlined as follows:
- (1)
Bridge test model fabrication: Determine the bridge type, structural design, and materials for the experimental bridge. Create detailed drawings for the Bridge test model and proceed to construct the truss bridge structure model.
- (2)
Bridge test model validation and adjustment: Establish finite element models of the bridge structure using software like Ansys 19.2 and Midas Civil 2021 (v1.1). Conduct static and dynamic analysis on both the bridge test model and the finite element model to validate the accuracy of the bridge test model.
- (3)
Presetting damage scenarios: Assign numerical identifiers to truss members based on their structural forms, considering the characteristics of the truss structure. Predefine multiple damage scenarios for the bridge structure.
- (4)
Acquisition of displacement response data: Utilize moving loads as the loading method and position sensors at the mid-span of the main beam. In the laboratory setting, apply simulated vehicle loads to the bridge structure. Utilize displacement sensors to capture the displacement-time responses of the bridge structure under both undamaged conditions and the preset damage scenarios. Acquire displacement response data of the bridge structure and apply wavelet packet analysis for noise reduction.
- (5)
Based on the measurement results from step (4), the displacement influence lines under the assumed damage condition should be subtracted from the displacement influence lines under the undamaged state. This subtraction yields the QSDIL damage index, which is used for the preliminary assessment of bridge structure damage and the rough localization of the damaged area.
- (6)
In step (5), after approximately localizing the damaged area, an increased number of sensors should be strategically placed in that region. These additional sensors are used to measure the dynamic response signals of acceleration time history at various points within the initially identified damage location.
- (7)
Localization of bridge structural damage: Utilize continuous wavelet transform for wavelet packet decomposition and calculate the RES damage indicator. This should be followed by the plotting of the RES damage curve. Precise spatial localization of structural damage in the bridge is achieved by analyzing the peaks and abrupt changes in the RES curve.
Figure 2.
Damage identification process flowchart.
Figure 2.
Damage identification process flowchart.
3. Experimental Design for Identifying Structural Damage in Bridges
3.1. Fabrication and Validation of Bridge Model
The selection of materials for bridge models should consider the experimental objectives and the applicable range of materials. Commonly used model materials include metal, gypsum, plastic, wood, and micro-concrete. Wood, compared to previous materials, possesses advantages such as simple material selection, easy model fabrication, high strength, and stability. In bridge damage identification experiments, structures made of wood are easily disassembled, and damaged beams can be readily replaced. Therefore, this paper adopts an experimental model with a wooden truss benchmark structure in its design.
Therefore, in this study, a wooden truss benchmark structure was chosen as the bridge test model. The structural dimensions are as follows: 1.8 m × 0.4 m × 0.3 m, with a main truss spacing of 0.3 m and a bay length of 0.3 m. The truss bridge has a span (L) of 1.8 m, a height (H) of 0.3 m, and a cross-sectional area of 100 mm
2. The members are bonded using hot-melt adhesive. The experimental test beam has a total length of 1.8 m, with an effective test segment of 1.6 m. There are 0.1 m of bridge deck on each side as a buffer zone for accelerating and decelerating moving loads. The bridge deck is constructed from a single piece of Japanese white pine, measuring 1.8 m in length, 0.4 m in width, and 0.005 m in thickness. The boundary conditions on both ends are mobile supports. By changing the cross-sectional shape of the truss beam at the damage location, the degree of damage is quantified. The bridge test model is shown in
Figure 3, and the relevant wooden material parameters are listed in
Table 1.
To validate the applicability of the bridge test model, a combined methodology involving experiments and numerical analysis was employed to comprehensively assess its static and dynamic characteristics. Distinct finite element models for the bridge structure were meticulously formulated using ANSYS and Midas software. In the ANSYS model, the upper and lower chord members, as well as the diagonal bracing elements, were discretized using Link8 elements, while the bridge deck was simulated using Shell181 elements. The comprehensive 3D finite element model of the bridge in ANSYS comprises a total of 470 elements and 375 nodes, as visually represented in
Figure 4a.
In the Midas model, the truss bridge components were defined utilizing custom materials fabricated from pine wood. The upper chords, lower chords, vertical members, and diagonal bracing elements have been represented as rod elements with a diameter of 10 × 10 mm. The entire 3D finite element model of the bridge within the Midas software encompasses 77 elements and 46 nodes, as elucidated in
Figure 4b.
To assess the behavior of the bridge, both a static load test and a dynamic resonance test were conducted. In the static load test, a 3 kg test weight was positioned at the mid-span of the bridge, and the deflection at this location was quantified using a laser displacement sensor (model LK-G3000, KEYENCE, Osaka, Japan). The dynamic resonance test aimed to determine the bridge’s fundamental frequency, and it employed a single-point excitation technique. The bridge underwent free vibration excitation by being struck with a rubber hammer. Vibration signals were then collected via piezoelectric accelerometers (model DH105E, Jiangsu Donghua Testing, Taizhou, China) positioned at critical locations, including the bridge’s supports, 1/4 span, mid-span along the central axis, and at the deck’s edge.
Subsequently, the collected acceleration data underwent spectral analysis to extract the bridge’s first-order vertical resonance frequency, which is graphically depicted in
Figure 5. The measured values of mid-span deflection and the first-order vertical resonance frequency are compared with their corresponding theoretical values, as presented in
Table 2.
The data provided in
Table 2 illustrate a close alignment between the measured values of mid-span deflection and the first-order vertical resonance frequency for the bridge test model when compared to the theoretical values. The variations observed are all within a margin of less than 5%. This strongly suggests that both the experimental and numerical models of the bridge are dependable and can be confidently employed for subsequent structural damage identification and analysis.
3.2. Experimental Trolley Model
The experimental vehicle is comprised of a metal body equipped with an FS-GR3E drive system. The dimensions of the vehicle’s body measure 0.2 m × 0.09 m × 0.05 m, with a total weight of 10 N. The rear cargo compartment of the vehicle is designed with dimensions of 7.5 cm × 3.5 cm × 3.5 cm, allowing for the placement of a specific number of counterweights to adapt to varying cargo loads. It is equipped with the FS-GT3C remote control system, which operates on a 3-channel 2.4 G transmission signal, providing controlled speed within a designated range. The visual representation of both the vehicle and the accompanying remote control system is provided in
Figure 6.
3.3. Instrumentation and Essential Equipment
The truss structure primarily exhibits vertical vibrations. Therefore, when performing signal extraction, only vertical (
Z-axis) accelerations need to be considered. Consequently, it is sufficient to deploy acceleration and displacement sensors in the vertical direction. To initially locate damage in the bridge structure, displacement sensing is required, necessitating the placement of a displacement sensor at the midspan of the bridge structure. After the preliminary damage localization, additional acceleration sensors are introduced to achieve precise localization of structural damage and identify the spatial damage locations. The arrangement of acceleration and displacement sensors is illustrated in
Figure 7.
The chosen data acquisition system for testing is the TZT3828E dynamic–static signal test and analysis signal acquisition instrument, along with its accompanying dynamic acquisition analysis software. The dynamic–static signal test and analysis signal acquisition instrument are depicted in
Figure 8. The operational block diagram of the multi-channel data acquisition system used is illustrated in
Figure 9.
3.4. Experimental System for Identifying Structural Damage in Bridges
The bridge structure damage identification test system consists of vehicle test models and bridge test models. The bridge test model is divided into an acceleration runway, a bridge test segment, and a deceleration runway. During the experiment, the vehicle can be controlled remotely to accelerate, decelerate, maintain a constant speed, and apply brakes. The vehicle accelerates from rest to the set speed in the acceleration runway, travels at a constant speed through the bridge test segment, and then brakes to a stop in the deceleration runway. The test system based on the bridge test model for bridge structure damage identification is depicted in
Figure 10.