Abstract

Structural pounding can cause considerable damage and even lead to collapse of structures. Most research focuses on modeling, parameter investigation, and mitigation approaches. With the development of structural health monitoring, the on-line detection of pounding becomes possible. The detection of pounding can provide useful information of potential damage of structures. This paper proposed using wavelet scalograms of dynamic response to detect pounding and examined the feasibility of this method. Numerical investigations were performed on a pounding system that consisted of a damped single-degree-of-freedom (SDOF) structure and a rigid barrier. Hertz contact model was used to simulate pounding behavior. The responses and pounding forces of the system under harmonic and earthquake excitations were numerically solved. The wavelet scalograms of acceleration responses were used to identify poundings. It was found that the scalograms can indicate the occurrence of pounding and occurrence time very well. The severity of the poundings was also approximately estimated. Experimental studies were carried out, in which shake table tests were conducted on a bridge model that underwent pounding between its different components during ground motion excitation. The wavelet scalograms of the bridge responses indicated pounding occurrence quite well. Hence the conclusions from the numerical studies were verified experimentally.

1. Introduction

The different phase vibrations of neighboring buildings or adjacent parts of the same building or bridge can result in pounding under earthquake excitation if the separation distance between them is not sufficient. The pounding can cause considerable damage and even lead to collapse of structures, for example, the bridges with columns of unequal heights [1, 2]. Damage of buildings and bridges due to pounding has been documented in the reports of many earthquakes by researchers. The existed pounding research concentrated on modeling of pounding systems, parameter investigations, and mitigation approaches.

Most mathematical pounding modeling methods fall into two categories: (a) stereomechanical impact approach, which is primarily based on impulse-momentum law; (b) contact element approach, a force-based approach. The first approach was used in [3] and other literature. Hertz contact model, a contact element approach, has been extensively used to model impact [47]. The numerical study portion of this paper used this model.

With the development of structural health monitoring systems in bridges and other civil structures, the on-line detection of the presence of pounding and parameter determination of the pounding system becomes possible. Valuable information about pounding can be extracted from the response measurements to provide early warning about damage to the structures. There is rare research on this useful topic based on the authors’ knowledge. Wavelets allow for powerful and flexible methods for dealing with fundamental engineering problems and have been used within a broader scope [8]. This paper proposed to use wavelet scalograms to detect the occurrence of pounding.

Wavelet scalogram is square of the modulus of the wavelet transform. It can uncover previously concealed information about the nature of nonstationary processes. It has been widely used in vibration signal analysis, which includes denoising, fault diagnostics, structural performance testing, ground motion analysis, bridge response due to vortex shedding, and many other forms of analysis [9, 10].

The primary objective of this paper is to investigate the applicability of wavelet scalograms in seismic pounding detection. Numerical and experimental studies were conducted to achieve this aim. In the numerical study, the Hertz contact model was introduced into the pounding modeling, and the properties of a large-scale bridge segment were used for an SDOF concrete structure while the appropriate value of an impact stiffness parameter was calculated. The responses of the pounding system between the SDOF structure subjected to harmonic and earthquake excitations and a stationary rigid structure were computed and processed to obtain wavelet scalograms. In the experimental studies, shake table tests were conducted on a steel bridge model. The pounding between the different units of the bridge was detected by analyzing the dynamic response measurements. The effectiveness of the wavelet scalogram in identifying seismic pounding was studied.

2. Investigation of the Pounding Model

2.1. Pounding Model

In the numerical study, pounding was simulated by an impact oscillator, which was proposed by Davis [5] and later adopted in the literature [4, 6, 7]. An SDOF damped oscillator subjected to ground motion as illustrated in Figure 1 was studied. The pounding occurred between the oscillator and a neighboring stationary rigid structure. The neighboring structure can also be flexible structures; however, for the sake of simplicity, this paper investigated the rigid case only.

The equation of motion for the nonlinear pounding model of the SDOF system subjected to ground motion in Figure 1 is where m denotes the system mass, c denotes the damping coefficient, , , denote the relative acceleration, velocity and displacement relative to the ground, the superimposed dot denotes differentiation with respect to time, and is the restoring force.

This study focused on an elastic structure. For this case, the damping coefficient and restoring force can be expressed as where and denote the damping ratio and initial stiffness, respectively.

in (1) is the impact force between the colliding structures. Using the Hertz contact model, the impact force is represented by a nonlinear spring of stiffness and can be expressed as where a is the separation distance between the colliding structures. The Hertz exponent 3/2 is a typical value to render (1) as nonlinear. While pounding in the real world may have different values, the alternative values have little effect on the system’s response [5, 6]. is the impact stiffness parameter, which is dependent on the material properties and geometric surfaces of the contact bodies [11]. The calculation of values of for some special cases can be found in [12]. For the pounding between concrete bodies, the values of typically range from 2~80 kN/mm3/2 as demonstrated in [11]. However, these values are based on small-scale experiments. They are not suitable for large-scale situations [7].

In the numerical examples used in this paper, the impact between a sphere of radius and massive plane surface of = was considered. The appropriate impact stiffness parameter is calculated as where , is the radius of the flexible structure in Figure 1, , , , and are the elastic modulus and Poisson’s ratios of the contact structures, respectively. For the neighboring rigid barrier, assuming , then , (4) becomes Equation (5) was used in Section 4 to calculate .

2.2. Solution of the Equation of Motion

The equation of motion (1) was solved by using the Dormand-Prince 4th-order-runge-kutta method with adaptive size as well as the Dormand-Prince 8th-order-runge-kutta method. When the results of these two methods were compared, it was found that the 8th order method was no better than the former method. Therefore, only the results of the 4th-order method were presented in this paper.

The solution method for the pounding system was verified by using a numerical model from [4]. The SDOF elastic system with mass = 350 Mg, elastic stiffness = 10.5 MN/m, damping ratio = 0.02, separation gap = 0.025 m, and impact stiffness parameter β = 80 kN/mm3/2 was used as a verification example. Two cases of sine excitations with 2 HZ and 0.67 HZ sharing a common magnitude of 0.1 g were each applied separately. The responses of displacement, acceleration, and pounding force calculated in this study were compared against those presented in [4]. They matched perfectly, so the method in this paper was verified to be accurate. This comparison study was not presented herein.

3. Wavelet Transform and Scalogram

The wavelets generated by dilation and translation from the analyzing wavelet are defined as where a is a scale parameter and b is a translation parameter.

and are usually called mother and son wavelets, respectively. It is supposed that and the analyzing wavelet satisfy the admissibility condition where is the Fourier transform. It can be derived from the previous conditions that,

Various wavelets have been developed, which include Haar, Morlet, Mexican Hat, Hermitian, and Daubechies wavelets. To select a proper mother wavelet is important. The Morlet wavelet is a complex function that can be considered as modulated Gaussian function by sinusoidal wave; it is a cosine signal decaying exponentially on both sides. It has the advantage of allowing trade between time and frequency resolution. It was proved to be an effective tool for feature extraction for mechanical vibration signals diagnosis [13, 14]. It has also been applied in seismic signal detection [15]. Morlet wavelet is similar to impulse component, which is symptom of damage in many dynamic signals. Because of these properties, the complex Morlet wavelet was chosen for this research, which is defined as where is a bandwidth parameter and is a wavelet center frequency [16]. = 1.5, = 1 are used for the wavelet scalograms in Sections 4 and 5 of this paper.

The wavelet transform is a linear transform that correlates a finite-energy function with a series of oscillating functions it is defined as where the asterisk denotes a complex conjugate. The modulus is defined as

The phase is defined as

The modulus and phase were used in [17] to investigate nonlinear and chaotic behavior of structural systems.

The wavelet scalogram is defined as the square of the modulus of the wavelet transform and represents energy density that measures the energy of in the Heisenberg box of each wavelet centered in [18]:

It is used for analysis of nonstationary signals and can be understood as a spectrum with a constant relative bandwidth. The scalogram is employed in this study to determine the occurrence of pounding for both numerical and experimental data.

4. Numerical Investigations

The results of pounding of an impact oscillator on its neighboring rigid barrier, as illustrated in Figure 1, were investigated. The properties of the shorter concrete bridge segment in [3] were used for the SDOF structure: mass m = 1.2 × 106 kg, elastic stiffness = 107 MN/m, damping ratio = 0.05, and a natural period T = 0.67 seconds. The other parameters were = 4.92 m from [7], Poisson ratio = 0.15, compressive strength (27.58 MPa), and consequently concrete strength = 3.605 × 106 psi = 24.856 GPa. Substitute the above values into (5) to calculate the impact stiffness parameter, and it was obtained that β = 757 kN/mm3/2. The separation distance was a = 12.5 mm. The structure was excited by two single frequency sinusoidal ground motions with frequencies = 4 Hz and = 0.8 Hz, respectively: = −0.1 g × sin; = −0.1g × sin (). 1940 El Centro North-South component scaled to peak ground motion (PGA) of 0.1 g and 0.05 g, respectively, were also employed as ground excitations. The duration of 31.14 seconds was used for all the cases.

The responses of the pounding system were computed by the procedures discussed in Section 2. The wavelet scalograms of the acceleration responses were calculated by the method expanded upon in Section 3. The frequency axes of the scalograms were in log scale. The objective was to find a method to detect the presence of pounding from the time history signals that may include the displacement/velocity/acceleration. Both the amplitude spectra and scalogram were calculated from acceleration responses. The pounding forces that were computed numerically addressed the actual condition but were treated as unknown throughout the detection process.

The analyses of responses due to stationary sinusoidal excitation are discussed first. Figure 2 shows the responses, single-sided amplitude spectra and scalogram for the case = 4 Hz, at which pounding forces were all computed to zero during this session, thus indicating no pounding presence. The displacement and acceleration were able to maintain the approximately sinusoidal shapes after the first few seconds of oscillation. The single-sided amplitude spectrum shows that the predominant frequency is approximately same to the excitation frequency of 4 Hz. The scalogram provdies a clear picture of energy distribution with respect to time and frequency.

Figure 3 shows the displacement, acceleration, pounding forces, single-sided amplitude spectra, and scalogram for the case = 0.8 Hz, which is closer to the natural frequency of the SDOF structure, and consequently generates larger displacement and acceleration responses. In Figures 3(a)3(b), the displacement and acceleration change to periodical signals after about 5 seconds. Although the responses were not regular as the nonpounding responses, there were no acceleration spikes and it was not certain whether pounding has occurred. The pounding forces as shown in Figure 3(c) indicate the exact time and duration of the poundings. The spectrum in Figure 3(d) shows many prominent frequencies, and additionally the predominant frequency is 1.6 Hz rather than the excitation frequency 0.8 Hz, which indicates that pounding may have occurred at this time. But this judgment needs a priori knowledge of the structure and the excitation. Compared with the nonpounding scalogram shown in Figure 2(d), the scalogram in Figure 3(e) has many abrupt changes of energy with respect to the frequency. When Figures 3(e) and 3(c) are compared, it is obvious that these changes occur at the time of pounding. Further, the changes of frequency contents indicate the magnitude of the pounding force. Therefore the scalogram of acceleration response can indicate the presence of pounding, the time when pounding occurs, and relative intensity of the pounding.

The dynamic responses due to nonstationary earthquake excitation were analyzed. Figures 4 and 5 show the analysis results due to earthquake excitation with different PGAs. The earthquake record of 1940 El Centro north-south component scaled to PGA of 0.1g is illustrated in Figure 4(a), which contains its most energy approximately in the frequency range of 0.2~10 Hz. It was observed that there are spikes in the acceleration histories shown in Figure 4(c), which indicates that the pounding may have occurred. The pounding forces shown in Figure 4(d) illustrate the presence of pounding due to this excitation. Unfortunately, the pounding cannot be detected by amplitude spectra very well. The wavelet scalogram for PGA = 0.1 g is illustrated in Figure 4(f), and the scalogram for PGA =0.05 g that does not generate pounding is shown in Figure 5(f). When Figures 4(f) and 4(d) are compared, it is clear that the higher frequencies determined by the scalogram are associated with the occurrence of pounding. The scalogram in Figure 5(f) for the nonpounding case does not show these types of changes, this additionally proves that the scalogram is capable of detecting pounding. It can be concluded that acceleration spikes and wavelet scalogram can detect pounding effectively.

5. Experimental Data Analysis

A series of shaking table tests were conducted on seismic pounding of a high-pier bridge by Meng et al. [19] at Institute of Engineering Mechanics, China Earthquake Administration in 2009. The measured accelerations were analyzed in this paper by using a wavelet scalogram to verify the ability of scalograms to detect pounding.

5.1. Shake Table Tests on Seismic Pounding of High-Pier Bridge Models

A high-pier bridge model with one main span and two approach spans was constructed as illustrated in Figure 6; the prototype of it was a high-pier bridge in Yunan, China. The superstructures and the 1.50 m tall piers were built with steel plates and pipes, respectively. The main span and approach spans had different sizes. Different ballast masses were added to satisfy the similitude laws. Consequently they had different natural frequencies. The gaps between the main span and the approach spans were 0.01 m.

The bridge model was placed on the three-dimensional shake table. Seven accelerometers, 11 displacement sensors, and 20 strain gauges were attached to the bridge model. Two sets of dynamic tests were conducted on this model. The first tests were performed on the previous described bridge model, and the second tests were implemented on the previous model installed with shape memory alloy metal rubber dampers (SMAMRDs). The SMAMRD reduced the pounding effects significantly. The applied ground motions simulated the El Centro, Taft earthquake and artificial ground motions for both test sets with same PGAs. In this study, the accelerations a4 from the El Centro ground motion tests were analyzed. The PGA of the El Centro table motion was 0.22 g for both the test sets with and without SMAMRD. A typical table acceleration is shown in Figure 7.

5.2. Pounding Detection by Wavelet Scalogram

The accelerations recorded by accelerometer a4 for both test sets with and without SMAMRD were analyzed by wavelet scalograms. The results of the two tests are shown in Figures 8 and 9 respectively. The analysis methods were the same as those used in nummrical investigations.

From the acceleration response, single-sided amplitude spectra, and the scalogram in Figure 8 for the test set with SMAMRD, it is observed that in the time ranges around 3 and 18 seconds, the peak values of acceleration become 2 times as those of other time ranges and the scalograms extended to higher frequencies 3 times as those of other time ranges. The changes in the scalogram together with the acceleration indicate clearly that the pounding occurs at the time ranges around 3 and 18 seconds. However, no clear pounding information can be drawn from the spectra.

From the response results shown in Figure 9 for the test set without SMAMRD, it was observed that in the time ranges of 3.5~9.2, 15.0, 18.5~22, and 24.0 seconds, the peak values of acceleration become about 3 times as those of other time ranges and the scalograms extended to higher frequencies about 5 times those of other time ranges. Therefore, it can be concluded that the pounding occurs at the above time ranges by analyzing the acceleration spikes and wavelet scalograms. However, the pounding information from the spectrum is inconclusive.

Both test sets in Figures 8 and 9 refer to the same El Centro ground motion with the same PGA. Comparing the times of pounding and the scalogram that represent the energy density in Figures 8(c) and 9(c), the times of pounding occurrence and the scalogram magnitudes in Figure 9(c) are more frequent and larger than those of Figure 8(c). The scalogram in Figure 9 is extended to much higher frequencies of 100 Hz compared to about 37 Hz in Figure 8. These results are consistent with the conclusion in [19] that the SMAMRD can reduce pounding substantially.

6. Conclusions

In this paper, the pounding detection by wavelet scalogram was examined by numerical studies as well as experimental studies.

In the numerical studies, the Hertz contact model was carefully integrated into the one-sided pounding system between an elastic SDOF concrete structure and a rigid barrier. The pounding forces, displacement, and acceleration of the pounding system were solved from the numerical models for both stationary sinusoidal and non-stationary earthquake excitations. Acceleration spikes were present in the earthquake excitation case, but not in the harmonic one. The wavelet scalogram proved its ability to identify pounding occurrences in both cases. It can also be used to indicate the time of pounding occurrence and the approximate severity of each pounding session.

In the experimental studies, the acceleration data from two sets of shake table tests on a steel bridge model with and without SMAMRD were analyzed. Pounding occurred between the different components of the bridge for both sets of tests. The wavelet scalogram of the accelerations detected the pounding effectively, which verified the conclusions from the numerical studies. The acceleration spikes were also able to indicate the pounding for this case.

It can be concluded that the wavelet scalogram is capable of indicating the presence of structural pounding quite well. The acceleration spikes can be used together with the wavelet scalogram for pounding detection. Further analytical and experimental studies on more practical cases are needed.