New Developments in Geotechnical Earthquake Engineering

Geotechnical seismic engineering is an important area in geotechnical engineering, and its topic is to resolve the seismic problems related to geotechnical engineering. With development for more than 50 years, the time domain analysis theory and the frequency analysis theory have been generally established in the geotechnical seismic engineering [1–6]. However, a common problem still exists in the time domain analysis theory and the frequency domain theory, which is that the research methods for geotechnical seismic engineering are still the separate application of time domain or frequency domain, and no combination of above two is considered. [7–11]. For example, time domain analysis theory can only consider the time histories of acceleration, velocity, and displacement but not the frequency contents of ground motions. Frequency domain analysis theory can only consider the frequency contents of ground motions, while time histories of acceleration, velocity, and displacement are not included. Seismic wave is however a very complex nonstationary signal, whose amplitudes and frequencies change with time. Present outcomes in the time domain theory and the frequency domain theory cannot reflect synthetically the characteristics of geotechnical seismic engineering, very close to an oldChinese saying “Theblindman feels an elephant—to take a part for the whole.” Therefore, research on the timefrequency-amplitude analysis theory considering the time, frequency, and amplitude can be an important task, which will be a new direction of geotechnical seismic engineering.


Review on the Present Main Defect in the Geotechnical Earthquake Engineering
Geotechnical seismic engineering is an important area in geotechnical engineering, and its topic is to resolve the seismic problems related to geotechnical engineering.With development for more than 50 years, the time domain analysis theory and the frequency analysis theory have been generally established in the geotechnical seismic engineering [1][2][3][4][5][6].However, a common problem still exists in the time domain analysis theory and the frequency domain theory, which is that the research methods for geotechnical seismic engineering are still the separate application of time domain or frequency domain, and no combination of above two is considered.[7][8][9][10][11].For example, time domain analysis theory can only consider the time histories of acceleration, velocity, and displacement but not the frequency contents of ground motions.Frequency domain analysis theory can only consider the frequency contents of ground motions, while time histories of acceleration, velocity, and displacement are not included.Seismic wave is however a very complex nonstationary signal, whose amplitudes and frequencies change with time.Present outcomes in the time domain theory and the frequency domain theory cannot reflect synthetically the characteristics of geotechnical seismic engineering, very close to an old Chinese saying "The blind man feels an elephant-to take a part for the whole." Therefore, research on the timefrequency-amplitude analysis theory considering the time, frequency, and amplitude can be an important task, which will be a new direction of geotechnical seismic engineering.

Brief Introduction of the Signal Analysis Technologies
Fourier transform and wavelet transform are two seismic signal analysis technologies [12].Fourier transform is a steady-state analysis technology and suitable for frequency domain analysis.Wavelet transform has some timefrequency domain resolution ability, but it is hard to carry out an accurate time-frequency analysis [13]   intrinsic mode function (known as IMF for short) components by using ensemble empirical mode decomposition (known as EEMD for short) method.The instantaneous frequency of each IMF will be derived by the Hilbert-Huang transform, and then the Hilbert spectrum will be obtained by integrating all the instantaneous frequency spectrums.The analytical signal, amplitude function, phase function, and instantaneous frequency function can be obtained by using the Hilbert-Huang transform.The computational formulas are shown in ()  () . ( From the above formulas, it is known that both amplitude and frequency are the function of time.If the amplitude is calculated by using a combined time and frequency domain approach, the Hilbert spectrum can be obtained, as shown in (2).The marginal spectrum and the instantaneous frequency spectrum can be obtained by using (3) and (4).The Hilbert spectrum can be obtained by integrating amplitude squared, as shown in (5).Consider

New Development in the Geotechnical Earthquake Engineering
Up to now, our team has conducted a field investigation along the road and railway nearly 3000 km in 5.12 Wenchuan earthquake-stricken areas (see Figure 1), a number of field monitoring of typical slopes, more than twenty shaking table tests, a number of numerical simulations and theoretical studies [15,16]; some shaking tabletest models are shown in Figures 2, 3, 4, 5, and 6.On the basis of extensive studies, several important problems in the geotechnical seismic engineering are studied, including analysis methods of seismic responses of the site, analysis methods of seismic responses and seismic stability of the retaining structures, analysis methods of seismic responses of the slope, analysis methods of the landslide mechanism, and the analysis method of seismic stability of the slope.With the help of above-mentioned signal analysis skills in the aerospace area and signal analysis area (e.g., Hilbert-Huang transform), our team has conducted a series of analyses on civil structure depending on the proposed time-frequency-amplitude (known as TFA for short) analysis theories, which are comprised of horizontal layers, inclined  layers, gravity retaining wall, reinforced retaining wall, slopes, and landslide mechanism [17,18].Our team gives the initial prototype of TFA analysis theories in geotechnical seismic engineering, which realizes the organic combination of time and frequency information, and develops a groundbreaking concept to the traditional idea in the geotechnical seismic engineering area.

Essence of the Time-Frequency-Amplitude Analysis Theories
The time-frequency-amplitude analysis theories give reasonable considerations to the time-frequency characteristics of seismic waves, mainly reflected in the elastic displacement magnitude and the frequency of the wave vector.This time, this paper selects Wolong earthquake wave (as shown in Figure 7) to make an explanation; procedures are as follows: first, Wo-Long wave is decomposed into intrinsic mode functions (IMF) using the empirical mode decomposition (EMD) method.The acceleration time histories of each IMF can be seen in Figure 8 and the instant frequency time histories of each IMF can be seen in Figure 9. Based on the acceleration time histories and instant frequency time histories of each IMF, bring them into the related formulas in order to get the stress time histories of each IMF; finally, sum the stress time histories together in order to obtain the total results.

Advantages of the Time-Frequency Analysis Method of Seismic Safety of the Rock Slope
Comparisons with pseudostatic methods using the various countries' earthquake resistant design codes, TFA analysis theory of seismic stability of rock slope is capable of considering the time, frequency, and amplitude of seismic wave, and then a brief introduction of the above-mentioned theory is made as follows.
In the TFA analysis theory of seismic stability of rock slope, the calculated model is shown in Figures 10 and 11.In the model, the sliding mass was divided into seven slices, and an angle of  between the tangent and horizontal directions at any point A was obtained at the middle of a slice.Figure 11 shows the reflection and refraction model at point A. The mechanical parameters in the upper medium (bedrock) are  1 ,  1 ,  1 , and  1 (density,  wave speed,  wave speed, and shear modulus) and  2 ,  2 ,  2 , and  2 in the lower medium (colluvial soils), and the I-I stands for the sliding surface.When transverse wave  1 reaches the I-I surface, some new waves are produced, for example, reflected transverse wave  2 , reflected longitudinal wave  3 , refracted transverse wave  4 , and refracted longitudinal wave  5 . 1 ,   1 ,  1 , and   1 stand for the incident angle of the incident transverse wave, the reflected angle of the reflected longitudinal wave, the refracted angle of the refracted transverse wave, and the refracted angle of the refracted longitudinal wave.Note that the reflected angle of the reflected transverse wave is the same as the incident angle of the incident transverse wave.After the derivation, the stress state at any point can be calculated according to ( 6) and (7).And then, the stress state at whole sliding surface can be calculated according to (8).At  last, the seismic stability of slope can be estimated according to (9). Consider In ( 6) and ( 7),  1 0 stand for the displacement of incident transverse wave,  2 0 stand for the displacement of reflected transverse wave,  3 0 stand for the displacement of reflected longitudinal wave,  4 0 stand for the displacement of refracted transverse wave, and  5 0 stand for the displacement of refracted longitudinal wave;  1  ,  1  stand for the wave vectors in the  and  directions of the incident, reflected, and the refracted waves, respectively;  stand for the cohension of the sliding mass and  stand for the internal friction angle.Consider In ( 8) and ( 9),   represents the sliding shear stress at point ;   = 0.5 ( −1 +  +1 ) × 1 is the slice area,  −1 is the width of the  − 1th slice and  +1 is the width of the  + 1th slice;  0 represents the normal stress at point ;   represents the frictional angle at point ;   represents the cohesion at point .
According to the above-mentioned method, a simple example is used to illustrate the advantage of the TFA analysis theory: model for calculation is shown in Figure 12.Figures 13 and 14 give the adopted sine wave whose instantaneous frequencies gradually increase with time and the calculated parameters can be seen in Table 1.
After the calculation is done, results of TFA analysis method are shown in Figure 15, indicating that the slope would fail after  = 10.0 s.The result of pseudostatic method shows that the surplus sliding force is −47 kN and means that the slope is stable.While the slope is proven to fail in the numerical simulation, which is consistent with the result of TFA analysis method shown in Figure 16, this phenomenon can be elucidated as follows.
(i) At  = 10.0 s, the instantaneous frequency of the seismic wave is 4.9 Hz and the natural frequency calculated according to the literature [17] is 6.8 Hz, with both being very close to each other, leading to a further resonance.However, the instantaneous frequency of the input seismic wave reaches 6.8 Hz at  = 18.2 s, and the resonance shows the strongest value, which leads to the maximum surplus shear force.Then the resonance gradually decreases with difference between the frequencies of the input wave and the natural frequency goes up, and the surplus shear force gradually decreases in response.
(ii) Time-frequency analysis method and numerical method both not only can consider the effect of PGA on the seismic stability of slope, but also consider the effect of frequency on the seismic stability of slope.But pseudo-static method only can consider the effect of PGA on the seismic stability of slope, which leads to the difference between the calculated results.
Synthesizes the above analysis, we can know that, therefore, the TFA analysis method can solve the fundamental flaw existing in the analysis method using the codes.

Application in Practice
Up to now, the TFA analysis theory has been already used to direct the aseismic designs of geotechnical engineering, containing several kilometers of YaLe expressway, and the contribution to several traffic lines opening.

Conclusions and Future Works
The main conclusions of this study are as follows.
(i) Based on the above analysis, this paper gives an initial prototype of time-frequency-amplitude analysis theory in the geotechnical seismic engineering.This theory realizes the organic combination of time information and frequency information and gives reasonable considerations to the time-frequency characteristics of seismic waves, which mainly is reflected in the elastic displacement magnitude and the frequency of the wave vector and that is the core content of TFA analysis theories.
(ii) Comparing to the results of TFA analysis theory of rock slope, pseudostatic method, and numerical analysis method shows that the TFA analysis theory can consider the effect of the time and frequency of the input wave on the seismic stability of slope, which solves the fundamental flaw existing in the analysis method using the codes.At the same time, according to the derivation results, we can know that the time-frequency analysis method of rock slope has some universalization and comprehensiveness, which can consider some factors reasonably as follows: (1) initial geostress field; (2) the inner frictional angle, the cohesion, the reflection coefficients, and the refraction coefficients; (3) the physical and mechanical properties of the bedrock and the regolith; and (4) the instantaneous frequency, the incidence angle, the time histories of the acceleration, the velocity, the displacement, and so forth.
(iii) The time-frequency-amplitude analysis theory is built in compliance with the elastic wave theory and geotechnical seismic engineering in which the elastic state of the system is considered.The next works are making the further studies based on the present outcomes; time-frequency-amplitude analysis theory can be thus proposed to consider the nonlinear of the geotechnical engineering and earthquake wave.

Figure 2 :
Figure 2: Shaking table test of site under nonuniform seismic excitation.

Figure 3 :
Figure 3: Shaking table test of rock slope.

Figure 4 :
Figure 4: Shaking table test of retaining wall.

Table 1 :
Physical and mechanics parameters of bedrock, soil layer, and structural plane.