Research Article | Open Access

# Research on Lateral Force of Pile Based on Liquefaction Site Effect

**Academic Editor:**Filippo Santucci de Magistris

#### Abstract

This paper presents a theoretical investigation on the lateral force of pile in liquefaction site. Regarded liquefied soil as fluid, the vector method can be used to analyze the liquefaction velocity field and solve the analytical solution of the dynamic field by making use of the principle of fluid mechanics. In addition, by solving the velocity field with vector symbol operation method, the analytical expression of the lateral force in the liquefied flow field is obtained, and the sensitivity of the parameters in the analytical expression is analyzed. The results show that the stress field of the pile contains both the pressure resistance caused by surface pressure and the friction resistance caused by shear stress, when the liquefied soil flows laterally. The lateral forces on the pile are mainly composed of inertial forces and damping forces and are related to density, fluid viscosity, pile radius, and vibration frequency. With the increase of density, fluid viscosity, and pile radius, the added mass and added damping increase gradually. In a certain range, added mass and added damping are sensitive to vibration frequency.

#### 1. Introduction

Sand liquefaction caused by almost every earthquake is accompanied by large deformation, and the lateral displacement caused by seismic liquefaction causes great damage to the pile foundation [1]. Hamada [2, 3] conducted in-depth field investigations and studied liquefaction sites and their effects in the light of the 1964 Kobe earthquake in Japan and the 1983 Chubu earthquake, in which most pile foundation damage in liquefaction sites was recorded. In 1994, the Northridge earthquake in the United States and the Hanshin earthquake in Japan in 1995 caused damage to the pile foundation of numerous bridges, buildings, storage tanks, wharves, and other structures in the liquefaction area. In addition, pile foundation damage in liquefaction sites was found among many earthquakes in Mexico, Costa Rica, India, Chile, and other countries. Cubrinovski et al. [4] reported in detail the seismic damage of bridge pile foundations in liquefaction sites in several earthquakes in New Zealand during 2010-2011. The Magnitude 9.0 earthquake that struck in the East Japan Sea on March 11, 2011, caused a large-scale liquefaction. A large number of piles on the banks of the Tone River are destroyed by site liquefaction [5]. In China, the Tangshan earthquake in 1976 caused the collapse of the pile foundation bridges such as the Yuehe Bridge, Shengli Bridge, and Luanhe Highway Bridge due to sand liquefaction and lateral displacement of the soil. In the 2008 Wenchuan earthquake, more than 10 bridges were subjected to different degrees of liquefaction damage [6].

The researchers mainly used the solid mechanics concept of pile-soil interaction to explain the failure of pile foundation in liquefiable site. It is believed that the inertial force of the upper structure is the main reason for the strength failure of the pile foundation. Ting [7] through the horizontal load test of the full-size single pile in the liquefaction site found that the soil pressure between the pile and soil under dynamic load is lower than the corresponding value under static load. Tokimatsu and Suzuki [8] used a large-scale shaking table model test to study the pore pressure response of the soil around the pile and its influence on the p-y curve. Yao et al. [9] used the shaking table model test to study the dynamic interaction behavior of pile-soil-superstructure system during liquefaction and used the centralized mass model to analyze the dynamic response. Choobbasti et al. [10] used the finite difference method to establish a model and studied the response of the pile-soil interaction system after soil liquefaction. Samui et al. [11] analyzed the lateral buckling of pile groups caused by the axial loads on pile foundation under the condition of soil liquefaction around the pile. Li and Motamed [12] established a two-dimensional nonlinear dynamic finite element model, studied the interaction of pile-soil system under the consideration of two-dimensional effective stress, and compared the analysis results of finite element model and test results. Wang et al. [13, 14] conducted a series of shaking table test and found that sand liquefaction will cause the foundation acceleration response to change and the maximum bending moment of the pile appears in the pile head position during the vibration process. Hussien et al. [15] carried out a series of centrifuge experiments, mainly studied the dynamic response characteristics of pile group in sand, and compared the bending moment distribution of pile foundations.

In recent years, lateral displacement caused by seismic liquefaction has caused great losses to pile. Due to the significant difference between the mechanical properties of the liquefaction soil and the solid soil, such as the low shear resistance, high fluidity, and other fluid properties, it is a novel research idea to study the liquefaction of saturated sand soil on the basis of fluid mechanics. The basic idea of adopting fluid mechanics to study the liquefaction effect is to view the liquefied sand as a kind of fluid and then to study the liquefaction effects such as flow deformation or soil-structure interaction by using the computational methods of mechanics.

Relevant scholars have done some researches and analysis and have obtained some useful research conclusions. Towhata et al. [16] took the liquefied soil as fluid for the first time in the analysis of large lateral deformation of liquefied soil. Based on the results of the shaking table test, the lateral deformation curve of liquefied soil was assumed to be the sinusoidal function curve with the maximum at the top and minimum at the bottom of the layer. After that, Hamada and Wakamatsu [17] carried out extensive surveys on the earthquakes that occurred in Japan in the 1980s and 1990s and conducted a series of shaking table tests to study the flow characteristics of liquefied soil and obtained two important conclusions: (1) Liquefied soil is highly close to a kind of pseudoplastic fluid. With the dissipation of pore water pressure after the earthquake, the stiffness of liquefied soil is recovered and gradually converted to solid. (2) The flow of liquefied soil meets the Reynolds similarity rate, and its flow velocity is proportional to the square root of the thickness of the liquefaction layer. These conclusions have played an important role in promoting the flow deformation analysis of liquefied soil, and their correctness has been proved by subsequent studies. Thereafter, Tamate and Towhata [18] regarded the liquefied soil as viscous fluid and analyzed the flow characteristic of liquefied sand by using the principle of minimum potential energy. Hadush et al. [19] treated liquefied sand as Bingham fluid with yield strength, proposed and proved a numerical method based on fluid dynamics theory to predict deformation caused by liquefied soil, and used this method to reproduce liquefaction flow by the shaking table test. Sawicki and Mierczyński [20] regarded liquefied sand as a fluid and studied the settlement of buildings in liquefied sand and the flow slip of the inclined slope by using the Navier–stokes formula, a classical equilibrium equation in fluid mechanics. Huang et al. [21] regarded liquefied sand as viscous fluid and analyzed the flow deformation of liquefied sand by using the VOF model. Monstassar and de Buhan [22] regarded liquefied sand as Bingham fluid with yield strength and developed a simplified one-dimensional numerical model to analyze the horizontal displacement of liquefied sand. Chen et al. [23] studied the hydromechanical properties of liquefied sand, and the results showed that the behavior of postliquefaction fine sand under monotonic axial compression loading can be described as the combination fluid model of a power-law time-independent fluid and a power-law time-dependent fluid.

In this paper, the liquefied soil is regarded as fluid, and the dynamic field analysis of liquefaction site is carried out based on the principle of fluid mechanics and the vector symbol operation method, and the analytical solution of the kinetic force field is solved. After obtaining the analytical expression of lateral force of pile considering liquefaction effect, it is known that the lateral force of the liquefaction site on pile is mainly composed of inertia force and damping force and is related to fluid density, fluid viscosity, pile radius, and vibration frequency. Finally, the sensitivity analysis of the parameters in the analytical solution is carried out to obtain the sensitivity degree of different parameters.

#### 2. Study on Dynamic Distribution Field of Liquefaction Site

##### 2.1. Analytical Solution of Velocity Field in Liquefaction Site

###### 2.1.1. Velocity Field Analysis Based on Vector Method

In the process of soil liquefaction, the soil changes from solid state to liquid state due to the increase of pore water pressure and the decrease of effective stress [24]. Therefore, by introducing the viewpoint of fluid mechanics, the velocity field of the liquefaction site considering the pile effect can be obtained, as shown in Figure 1.

Assume the radius of cross section of the pile is *R*. The density and viscosity of soil, when the liquefaction flow state occurs, are and , respectively.

The velocity of the pile in the liquefied flow field is , where is the vibration frequency of the pile. represents the velocity of the pile at the initial moment of soil liquefaction, which is a constant vector. Assume the liquefaction flow of soil is infinite. Thus, when its velocity is , the lateral velocity field around the pile is . The flow velocity field of the liquefaction site caused by the pile movement is , and then .

The liquefied soil can be viewed as incompressible fluid and can be obtained according to the following continuous equation:

According to the vector theory, the polar vector corresponds to the axial vector , to guarantee that comes true, where is the gradient, so we can obtainwhere is the axial vector, which can be expressed as the vector product of two polar vectors. Since is only related to the radius and velocity of the polar vector, must be composed of , and can be expressed by the gradient . For the slender structure of pile, it can be considered that is independent of the axial position coordinate *z*; that is, has the vector product form of . In addition, since and is a constant vector, then

Therefore, the velocity field caused by the vibration of the pile at the velocity in the fluid field is

###### 2.1.2. Solution of the Velocity Field with Vector Symbol Operation Method

Since the liquefied soil is an incompressible fluid, the momentum equation can be obtained as follows:where is the fluid viscosity of liquefied soil, is the density of liquefied soil, and is the stress field around the pile, is Laplacian.

Since the fluid motion caused by soil liquefaction belongs to a uniform flow with low Reynolds number, the inertia term in the formula can be ignored, and the linearized N–S equation can be obtained:

Take the curl from left to right of the above equation, and because , we can obtain

Substitute equation (4) into equation (7) and use the vector symbol operation to obtain

Then

It is known from the boundary conditions at infinity that ; therefore,and the first derivative can be solved: , where .

Velocity field can be obtained by using the vector symbol operation method in the cylindrical coordinate system:

According to the boundary condition , ; thus, and .

According to equation (11), the components of in the normal, tangential, and axial directions are, respectively,

##### 2.2. Analytical Solution of Stress Field in Liquefaction Site

Since flowing liquefied soil is an incompressible fluid, according to the Navier–Stokes equation and stress tensor formula, equations (3), (4), and (12)–(14) can be substituted into it to obtain

Therefore, the pressure distribution around the pile is

Thus, the stress field around the pile can be obtained as

According to the calculation results, , , and are nonzero values. It can be found that the stress field around the pile includes both the pressure resistance caused by the surface pressure and the friction resistance caused by the shear stress.

#### 3. Research Studies on Lateral Force of Pile under Liquefaction Site Flow Condition

##### 3.1. Analytical Solution of Lateral Force of Pile

The lateral force applied to the pile is parallel to the velocity of soil liquefaction flow. Therefore, all surface forces acting on the surface of the pile are projected to the direction of fluid velocity, and integrate along the whole cylinder of the pile. The lateral force acting on the unit length pile is

On the pile cylinder , according to formula (16)–(19),

According to formula (23) and according to the definition of fluid mechanics principle and , by integration,where is the fluid viscosity of liquefied soil, is the density of liquefied soil, is the radius of pile, and is the vibration frequency of pile.

Formula (25) shows that the lateral force of the liquefaction site on pile is mainly composed of inertia force and damping force, and its dynamic characteristics are related to density, fluid viscosity, pile radius, and vibration frequency of pile. Supposewhere is the dynamically added mass of pile and is the dynamic added damping of pile. Both of these factors reflect the dynamic characteristics of pile under transverse vibration in the liquefaction flow field.

Finally, the lateral force of liquefaction site on the pile can be expressed as

##### 3.2. Parametric Sensitivity Analysis of Analytical Solutions of Lateral Force

###### 3.2.1. Method of Sensitivity Analysis

Sensitivity is defined as the derivative of a function with respect to a parameter; at the parameter , sensitivity of dynamic added mass function or dynamic added damping function with respect to iswhere reflects the monotonicity of function or on parameter . The value of indicates the sensitivity of each function to the parameter . The larger the value, the more sensitive the function is to .

The sensitivity of density, fluid viscosity, pile radius, and vibration frequency is analyzed below, and take the reference value of each parameter as density of liquefied soil , fluid viscosity , radius of pile , and vibration frequency of pile . During the sensitivity analysis of each parameter, the other parameters remain unchanged based on the change of a single parameter.

###### 3.2.2. Density of Liquefied Soil

The curves of the added mass and added damping with the density of liquefied soil are shown in Figure 2. The added mass and added damping increase with fluid density, and it indicates that with the increase of fluid density, the lateral force of liquefied soil on pile increases gradually. The sensitivity decreases with the increase of density, indicating that fluid will become more sensitive to lateral force when the density of fluid becomes smaller. The variation range of sensitivity corresponding to added damping is obviously larger than that corresponding to added mass, which proves that the influence of fluid density on added damping is greater.

**(a)**

**(b)**

###### 3.2.3. Fluid Viscosity of Liquefied Soil

The curves of added mass and added damping with the viscosity of liquefied soil are shown in Figure 3. The added mass and added damping increase with the increase of fluid viscosity, and the sensitivity decreases with the increase of fluid viscosity. It implies that the smaller the fluid viscosity, the more sensitive it is to lateral forces. And when the viscosity is between and , and the corresponding sensitivity drops sharply. This indicates that the change of fluid viscosity has a great influence on the lateral force. The variation range of sensitivity of added damping is larger, which indicates that the influence of fluid viscosity on added damping is greater.

**(a)**

**(b)**

###### 3.2.4. Pile Radius

The curves of added mass and added damping with pile radius are shown in Figure 4. It can be found that as the pile radius increases, the added mass and added damping increase gradually. The sensitivity corresponding to added mass increases accordingly, while the sensitivity corresponding to added damping remains unchanged, indicating that the pile radius is sensitive to added mass, but not to added damping.

**(a)**

**(b)**

###### 3.2.5. Vibration Frequency of Pile

The curves of the added mass and added damping with the vibration frequency of the pile are shown in Figure 5. With the increase of vibration frequency, the added mass decreases sharply and tends to be stable gradually. The added damping increases gradually, and the sensitivity corresponding to the added mass and added damping decreases rapidly at the low frequency stage and tends to a smaller value, indicating that when the frequency is less than 1 Hz, the vibration frequency is particularly sensitive to the lateral force.

**(a)**

**(b)**

##### 3.3. Comparative Analysis of Analytical and Numerical Solution of the Distribution Field around the Pile

Using the CFD fluid module of ABAQUS finite element platform, Abaqus/CFD provides advanced computational fluid dynamics capabilities with extensive support for preprocessing and postprocessing provided in Abaqus/CAE. These scalable parallel CFD simulation capabilities address a broad range of nonlinear coupled fluid-structural problems. Abaqus/CFD can solve the following types of incompressible flow problems: the internal or external flows that are steady-state or transient, span a broad Reynolds number range, and involve complex geometry that may be simulated with Abaqus/CFD. This includes flow problems induced by spatially varying distributed body forces.

A fluid finite element model is established to carry out numerical calculation of liquefied dynamic field, based on the principle of minimum potential energy and basic control equation of flow. A site is adopted to simulate the flow effect of the liquefaction site. The schematic diagram of the site is shown in Figure 6.

**(a)**

**(b)**

For the boundary conditions, in order to better simulate the infinite dynamic flow field, in the whole liquefied soil flow area, all the boundaries except the fluid outlet are set as the velocity inlet. Some parameters used in numerical calculation are as follows: the CFD model is adopted for liquefied soil mass, and thus, (the soil density) and (dynamic viscosity). The elastic model is adopted for pile, and thus, (concrete pile density), (elasticity modulus), and (Poisson’s ratio). In the calculation, the initial stress setting and balance of the soil are firstly carried out to obtain the stress distribution in the preliquefaction state, and then the liquefied soil is set as the fluid material for the dynamic field calculation.

The comparison of the analytical and numerical solutions of the stress distribution field around the pile is shown in Figure 7. It can be seen that the distribution trend of the analytical solution of the stress field is basically the same as that of the finite element solution around the pile and the stress on the pile structure is mainly distributed in the head-on position of liquefied soil flow. This result is consistent with the stresses around the pile in the process of seabed consolidation movement in Literature [25].

#### 4. Conclusion

Regarding the liquefied soil as fluid, the dynamic field analysis of liquefaction site is carried out based on the principle of fluid mechanics and the vector symbol operation method, and the analytical solution of the kinetic force field is solved. The analytical expression of the lateral force of pile is obtained considering the liquefaction effect, the sensitivity analysis of the parameters in the analytical solution is carried out, and then the influence of different parameters on the lateral force of pile is obtained. The main conclusions are as follows:(1)According to the analytical solution calculated by the liquefaction dynamic field theory, it can be known that when the liquefied soil flows laterally, the stress field around the pile includes the pressure resistance caused by the surface pressure and the friction resistance caused by the shear stress.(2)The lateral force of liquefied soil on pile is mainly composed of added inertia and added damping, which are related to density, fluid viscosity, pile radius, and vibration frequency.(3)With the increase of density, fluid viscosity, and pile radius, the added mass and added damping are gradually increasing. The influence of density and fluid viscosity on the added damping is greater, and the pile radius is only sensitive to the additional mass. Vibration frequency is particularly sensitive to added mass and added damping within a certain range.

#### Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

#### Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

#### Acknowledgments

The authors gratefully acknowledge the financial support for this study from the project of the National Natural Science Foundation of China (Grant no. 51408281), the Natural Science Foundation of Jiangsu Province (Grant no. BK20140108), and Qing Lan Project (2016).

#### References

- H. S. Liu, F. P. Xu, P. C. Li et al., “Lateral spread and large ground displacement caused by earthquake liquefaction,”
*Special Structures*, vol. 14, no. 2, pp. 47–50, 1997. View at: Google Scholar - M. Hamada, “Large ground deformations and their effects on lifelines: 1964 Niigata earthquake: case studies of liquefaction and lifelines performance during past earthquake,” Tech. Rep., Japanese Case Studies, National Centre for Earthquake Engineering Research, Buffalo, NY, USA, 1992, Technical Report NCEER-92-0001. View at: Google Scholar
- M. Hamada, “Large ground deformations and their effects on lifelines: 1983 Nihonkai-Chubu earthquake, case studies of liquefaction and lifelines performance during past earthquake,” Tech. Rep., Japanese Case Studies, National Centre for Earthquake Engineering Research, Buffalo, NY, USA, 1992, Technical Report NCEER-92-0001. View at: Google Scholar
- M. Cubrinovski, J. Haskell, A. Winkley, K. Robinson, and L. Wotherspoon, “Performance of bridges in liquefied deposits during the 2010-2011 Christchurch (New Zealand) earthquakes,”
*Journal of Performance of Constructed Facilities*, vol. 28, no. 1, pp. 24–39, 2014. View at: Publisher Site | Google Scholar - D. Pradel, J. Wartman, and B. Tiwari, “Impact of anthropogenic changes on liquefaction along the Tone River during the 2011 Tohoku earthquake,”
*Natural Hazards Review*, vol. 15, no. 1, pp. 13–26, 2014. View at: Publisher Site | Google Scholar - Z. H. Wu and Z. C. Zhang, “Types and examples of geological hazards in Wenchuan Ms 8.0 earthquake,”
*Acta Geologica Sinica*, vol. 82, no. 12, pp. 1747–1757, 2008. View at: Google Scholar - J. M. Ting, “Full-scale cyclic dynamic lateral pile responses,”
*Journal of Geotechnical Engineering*, vol. 113, no. 1, pp. 30–45, 1987. View at: Publisher Site | Google Scholar - K. Tokimatsu and H. Suzuki, “Power pressure response around pile and its effects on p-y behvaior during soil liquefaction,”
*Soils and Foundations*, vol. 44, no. 6, pp. 101–110, 2004. View at: Publisher Site | Google Scholar - S. Yao, K. Kobayashi, N. Yoshida, and H. Matsuo, “Interactive behavior of soil-pile-superstructure system in transient state to liquefaction by means of large shake table tests,”
*Soil Dynamics and Earthquake Engineering*, vol. 24, no. 5, pp. 397–409, 2004. View at: Publisher Site | Google Scholar - A. J. Choobbasti, M. Saadati, and H. R. Tavakoli, “Seismic response of pile foundations in liquefiable soil: parametric study,”
*Arabian Journal of Geosciences*, vol. 5, no. 6, pp. 1307–1315, 2012. View at: Publisher Site | Google Scholar - P. Samui, S. Bhattacharya, and T. G. Sitharam, “Support vector classifiers for prediction of pile foundation performance in liquefied ground during earthquakes,”
*International Journal of Geotechnical Earthquake Engineering*, vol. 3, no. 2, pp. 42–59, 2012. View at: Publisher Site | Google Scholar - G. Li and R. Motamed, “Finite element modeling of soil-pile response subjected to liquefaction-induced lateral spreading in a large-scale shake table experiment,”
*Soil Dynamics and Earthquake Engineering*, vol. 92, pp. 573–584, 2017. View at: Publisher Site | Google Scholar - R. Wang, J. M. Zhang, and G. Zhang, “Centrifuge shaking table test on single pile lateral spreading soil,”
*Engineering Mechanics*, vol. 29, no. 10, pp. 98–105, 2012. View at: Google Scholar - R. Wang, X. Liu, and J.-M. Zhang, “Numerical analysis of the seismic inertial and kinematic effects on pile bending moment in liquefiable soils,”
*Acta Geotechnica*, vol. 12, no. 4, pp. 773–791, 2016. View at: Publisher Site | Google Scholar - M. N. Hussien, T. Tobita, S. Iai, and M. Karray, “Soil-pile-structure kinematic and inertial interaction observed in geotechnical centrifuge experiments,”
*Soil Dynamics and Earthquake Engineering*, vol. 89, no. 10, pp. 75–84, 2016. View at: Publisher Site | Google Scholar - I. Towhata, Y. Sasaki, K. Tokida, H. Matsumoto, Y. Tamari, and K. Yamada, “Prediction of permanent displacement of liquefied ground by means of minimum energy principle,”
*Soils and Foundations*, vol. 32, no. 3, pp. 97–116, 1992. View at: Publisher Site | Google Scholar - M. Hamada and K. Wakamatsu, “A study on ground displacement caused by soil liquefaction,”
*Doboku Gakkai Ronbunshu*, vol. 43, no. 596, pp. 189–208, 1998. View at: Publisher Site | Google Scholar - S. Tamate and I. Towhata, “Numerical simulation of ground flow caused by seismic liquefaction,”
*Soil Dynamics and Earthquake Engineering*, vol. 18, no. 7, pp. 473–485, 1999. View at: Publisher Site | Google Scholar - S. Hadush, A. Yashima, R. Uzuoka, S. Moriguchi, and K. Sawada, “Liquefaction induced lateral spread analysis using the CIP method,”
*Computers and Geotechnics*, vol. 28, no. 8, pp. 549–574, 2001. View at: Publisher Site | Google Scholar - A. Sawicki and J. Mierczyński, “On the behaviour of liquefied soil,”
*Computers and Geotechnics*, vol. 36, no. 4, pp. 531–536, 2009. View at: Publisher Site | Google Scholar - Y. Huang, W. Mao, H. Zheng, and G. Li, “Computational fluid dynamics modeling of post-liquefaction soil flow using the volume of fluid method,”
*Bulletin of Engineering Geology and the Environment*, vol. 71, no. 2, pp. 359–366, 2012. View at: Publisher Site | Google Scholar - S. Montassar and P. de Buhan, “Numerical prediction of liquefied ground characteristics from back-analysis of lateral spreading centrifuge experiments,”
*Computers and Geotechnics*, vol. 52, no. 32, pp. 7–15, 2013. View at: Publisher Site | Google Scholar - G. X. Chen, E. Q. Zhou, Z. H. Wang et al., “Experimental study on fluid characteristics of medium dense saturated fine sand in pre- and post-liquefaction,”
*Bulletin of Earthquake Engineering*, vol. 14, no. 8, pp. 2185–2212, 2016. View at: Google Scholar - Z. H. Wang, E. Q. Zhou, G. X. Chen et al., “Characteristics of solid-liquid phase change of saturated sand under cyclic loading,”
*Chinese Journal of Geotechnical Engineering*, vol. 34, no. 9, pp. 1604–1610, 2012. View at: Google Scholar - T. Sui, J. Zheng, C. Zhang et al., “Consolidation of unsaturated seabed around an inserted pile foundation and its effects on the wave-induced momentary liquefaction,”
*Ocean Engineering*, vol. 131, pp. 308–321, 2017. View at: Publisher Site | Google Scholar

#### Copyright

Copyright © 2020 Zuo Xi and Zhou Enquan. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.