Quantification of Dynamic Properties of Pile Using Ensemble Empirical Mode Decomposition

*is paper investigated dynamical interactions between pile and frozen ground by using the ensemble empirical mode decomposition (EEMD)method. Unlike the conventional empirical mode decomposition (EMD)method, EEMD is found to be able to separate the mode patterns of pile response signals of different scales without causing mode mixing. *e identified dynamic properties using the EEMDmethod are more accurate than those obtained from conventional methods. EEMD-based results can be used to reliably and accurately characterize pile-frozen soil interactions and help designing infrastructure foundations under permafrost condition.


Introduction
To quantify dynamical properties of structures, timefrequency analysis can be used to effectively reveal the constituent frequency components of nonstationary signals and their time variation features, as well as transient events such as impulses.To date, various time-frequency analysis methods have been proposed [1,2].However, the inherent drawbacks of these time-frequency analysis methods limit their effectiveness in analyzing complex vibration signals.
For instance, linear transforms such as the short time Fourier transform (STFT) and wavelet transform are subject to the Heisenberg uncertainty principle, that is, the best time localization and highest frequency resolution cannot be achieved simultaneously.One of them can only be enhanced at the expense of the other, and hence, the time-frequency resolutions of linear transforms are limited.In addition, the basis in either the Fourier or wavelet transform is fixed.erefore, the conventional methods lack adaptability in simultaneously matching the complicated components inherent in complex vibration signals, such as the fundamental frequency and their harmonics, combined harmonics, chaos, system-induced impulses, and other transient vibration.As a typical representative of bilinear time-frequency representations, Wigner-Ville distribution has the best time-frequency resolution, but it has the inevitable cross-term interferences for multiple component signals.Such cross-term interferences complicate the interpretation of signal features in the timefrequency domain and make it unsuitable to analyzing the complex vibration signals.Various modified bilinear timefrequency distributions including Cohen and affine class distributions may suppress the negative effect of cross-terms but will compromise time-frequency resolution and autoterm integrity.Spectrogram better preserves energy distribution (comparing to wavelet transform), but it has issues with time/frequency resolution and gives rather rough estimates.Wavelet transform and scalograms provide good results but add artifacts on frequencies close to zero or close to the edges.Also, the conventional methods distort energy distribution with high-frequency components having less power comparing to lower ones.
e well-known Hilbert-Huang transform has fine timefrequency resolution and is free of cross-term interferences, but it essentially relies on the empirical mode decomposition (EMD) using spline interpolation [2].It is susceptible to singularities in signals, and may produce pseudointrinsic mode functions (IMF), thus masking or interfering the timefrequency structure of true signal components.Oscillations of di erent amplitudes are found in a mode, or similar oscillations are encountered in di erent modes.To avoid this problem, Wu and Huang proposed ensemble EMD (EEMD) [3], a method based on the EMD algorithm.e proposed method follows a study of the statistical characteristics of white noise, involves a noise-assisted analysis, and adds white noise of a uniform frequency distribution into EMD to avoid mode mixing.
In the Northern area, there are some bridges damaged directly due to the occurrence of frozen ground.Recently, some research has identi ed that the e ects of frozen ground on seismic response of bridge foundations are signi cant, which has validated the requirements for design guidelines to include the consideration of the frozen ground e ects [4][5][6][7].Seasonally frozen ground across the Northern area is signi cantly sti er than unfrozen ground.For bridges supported on deep foundations, the overall bridge sti ness increases measurably in winter season [7].As bridge sti ness increases, seismic demand increases.To address signi cant changes in bridge pier boundary conditions, additional knowledge of the subsurface is needed.Extra information is needed to address a bridge's ductile performance in bridge design for earthquake prevention event.e study of the e ects of frozen ground on foundation has been required for bridge designers.Several studies have been performed to evaluate the soil-foundation-structure interaction [4,5,8,9].
ere always has been a need for engineers to reliably and accurately quantify dynamical interactions between foundation and frozen ground by using the latest method.
As a relatively new technique, EEMD has been successfully applied in some engineering problems.For example, varied fault patterns in machine bearings have been identi ed using EEMD, whereas these fault frequencies could not be identi ed using fast Fourier transform (FFT) [10].In [11], EEMD is applied to extract speci c components from a data set of passenger ows of a station in the metro system and illustrates the time variants of short-term passenger ow for these data sets.e EEMD method is used to analyze carbon price pattern and extract carbon price signals from a multiscale viewpoint [12].Also, EEMD is used for trend prediction of market indexes [13].
In this study, the pile dynamic behavior in frozen ground condition has been determined by using the EEMD method.A test pile was constructed and tested for free vibration in winter conditions in Alaska to assess the e ects of frozen soils on the performance of single piles.To accurately quantify the dynamic properties of pile on frozen soil, we use EEMD to decompose the transient response of pile into serial components and then use these components to identify the system dynamic properties including natural frequencies and mode damping.e study provided a basis for accurately and reliably analyzing the frozen ground e ects on a pile structure and for designing bridge foundations.

Test Setup and Measurements
e test site is at mile 0.9 Farmers Loop Road at Fairbanks in Alaska.Some mechanical properties of soil at the test site (Fairbanks silt) are documented in a research report and a journal paper [14,15].Preliminary investigation revealed a thin layer of gravel overlaying silt and a permafrost table located 20 to 22 feet (6.096 m to 6.7056 m) below ground surface.
A single-pile eld test was conducted by using 25-footlong (7.62 m), 16-inch diameter (40.64 cm) steel-jacketed reinforced concrete piles embedded approximate 20 feet (6.096 m) in silt.Eight Number 6 ASTM A706 reinforcing bars provided 1.93% of the cross-sectional area.Shear reinforcement consisted of Number 3 reinforcing bars bent into rings (or hoops) and spaced on 12-inch (30.48 cm) centers outside of the longitudinal reinforcing bars; the outside diameter of each hoop was 1 inch (2.54 cm) less than the inside diameter of the steel pipe jacket.e pipe wall thickness is 0.375 inch (0.9525 cm).
e purpose of this study was to evaluate pile behavior during winter month.
e dynamic measurement was used to determine nature frequency and damping for the test pile.
e eld test was performed in the month of January, and the ground temperature was 1.94 °F (−16.7 °C).e pile tops were a free end, a lateral load was applied to the pile, and the load was quickly released with a quick disconnect device.
e accelerometer was used to measure how the pile vibrated once the load was released (Figure 1).

Response Signal and Spectrum Analysis
Figure 2 shows the measured acceleration signal.Figure 3 shows the FFT spectrum of this acceleration signal, and there are two speci c peaks f 1 101 Hz and f 2 1072 Hz. Figure 4 shows the spectrogram of the measured acceleration signal.From Figure 4(a), we can see the dominant component  Advances in Civil Engineering corresponding to speci c frequency f 1 101 Hz; f 2 only exists a short period decay to zero after 0.1 second.is can be attributed as the initial nonlinear interaction of pile and soil; the other parts have harmonics and the variation of frequency which suggest the existence of nonlinear e ect.Additional, from the 3D plot of Figure 4(b), we can see that the amplitude of the dominant component with speci c frequency f 1 101 Hz is substantially higher than the rest of components.e system's response can be approximated by this component.

Parameter Identification by Using EEMD Method
To accurately quantify the signal which has been roughly identi ed in the time-frequency spectrogram, this study used the ensemble empirical mode decomposition (EEMD) method which is developed based on the class empirical mode decomposition (EMD) method.Huang et al. [2] established empirical mode decomposition (EMD) which is a time-series demodulation method to extract a signal to a combination of intrinsic simple oscillatory signals (IMF) [8].e original signal x(t) is decomposed into n intrinsic modes c i (t) and a residue r n (t).Assume that the original signal can be represented by e detailed procedure and equations have been included in the published paper [8].However, using this method, there exists a mode mixing problem.Hu man et al. [9] proposed a noise-assisted data analysis method referred to as ensemble EMD (EEMD).In brief, the steps of EEMD method are as follows: rst, add white noise series m i (t) to the initial signal x(t): Second, EMD is applied to decompose the signal to IMFs, c j (t) [9].en, repeat the rst step with di erent white noises for each trail.Finally, estimate the means of corresponding IMF of the decompositions as the desired output: where i is the ensemble number and j is the component number [9].Next, apply the Hilbert transform to all IMFs, to derive model parameters including frequency and damping: After the Hilbert transform, H[c j (t)] and c j (t) form a complex signal Z j (t): en, the envelop of every IMF can be given by in which a j (t) the instantaneous amplitude of x(t) re ects how the energy of x(t) varies with time.e term φ j (t) is the instantaneous phase of x(t).e instantaneous frequency ω(t) is de ned as the time derivative of the instantaneous phase φ(t) as follows: en the original signal x(t) can be expressed as In principle, the measured acceleration response of the pile, € w, can be approximately decomposed by the EMD as follows: Advances in Civil Engineering where € w j (t) is the jth modal acceleration response and c i (t) is the ith IMF: us, the damped natural frequency ω dj can be obtained from the slope of the phase angle plot in ω j (t) versus t, and ς j can be obtained from the slope of the plot in lna ij versus t.
e linear least-squares method can be used to t the plots of lna ij versus time and ω j (t) versus time.e following section illustrates the signal decomposition and parameter identi cation by using the measured acceleration    Advances in Civil Engineering signal.Figure 5 shows the decomposed components of the measured acceleration signal (Figure 2).e spectrum analysis of IMF2 and IMF4 demonstrated that their respective dominant peaks are 1072 Hz and 101 Hz, respectively, which means IMF2 and IMF4 represent the two main natural modes ( f 2 and f 1 ) of the system, respectively.
In order to identify the dynamic properties of the decomposed signal, the obtained plots of phase angle and the logarithmic amplitude versus time for IMF2 and IMF4 are calculated and shown in Figure 6.
e damped natural frequency can also be obtained from the slope of the phase angle plot, and the damping can be obtained from the slope of the amplitude versus t. e linear least-squares method can be used to t the plots.
e damping coe cient corresponding to the mode 101 Hz is identi ed as ς 1 0.015.And the damping coe cient corresponding to the mode 1072 Hz is identi ed as ς 2 0.035.To further consider the complex frequency properties, the Hilbert spectrum is applied to retreat the acceleration signal recorded.Figure 7 shows the Hilbert spectrum of the decomposed signal intrinsic mode functions.From Figure 7, we can see that those respective components are well separated, and the timedependent properties of the frequencies are clearly shown.e dominant mode exhibits modulation centered at about 100 Hz.Advances in Civil Engineering

Conclusion
is paper applied the EEMD method to evaluate the soilpile structure interaction under frozen condition.e method can decompose the vibration signal into different scales and help to identify the natural frequency and damping parameters accurately.
is study identified the system parameter for the frozen soil-pile interaction.e test pile has two predominant natural frequencies: 101 Hz and 1072 Hz, and the corresponding modal damping are 0.015 and 0.035, respectively.e identified parameter can provide bridge engineering the required information for bridge design in a permafrost region.e application of EEMD has potential for other soil-pile vibration signals or other fields.

Figure 1 :
Figure 1: Elevation: test con guration for free-vibration response.e components are (A) soil surface, (B) reaction pile, (C) test pile, (D) chain hoist, (E) quick disconnect, (a) load cell, (b) accelerometer, and (c) displacement measurement.is sketch is not to scale.

Figure 2 :Figure 3 :
Figure 2: Time history of a measured acceleration signal.

Figure 5 :
Figure 5: Decomposed components of a measured acceleration signal.

Figure 6 :Figure 7 :
Figure 6: Characteristic plots of the decomposed signals corresponding to two speci c frequencies: (a) f 2 1072 Hz and (b) f 1 101 Hz.