Next Article in Journal
Therapy of Neurophysiological Changes after Oral and Maxillofacial Surgery—A Systematic Review
Next Article in Special Issue
Prediction of Pile Bearing Capacity Using XGBoost Algorithm: Modeling and Performance Evaluation
Previous Article in Journal
Analytical Model for Enhancing the Adoptability of Continuous Descent Approach at Airports
Previous Article in Special Issue
Responses of Laterally Loaded Single Piles Subjected to Various Loading Rates in a Poroelastic Soil
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Dynamic Lateral Response of the Partially-Embedded Single Piles in Layered Soil

1
School of Civil Engineering, Henan University of Science and Technology, Luoyang 471023, China
2
Engineering Technology Research Center of Safety and Protection of Buildings of Henan Province, Luoyang 471023, China
3
College of Civil Engineering, Hunan University, Changsha 410082, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2022, 12(3), 1504; https://doi.org/10.3390/app12031504
Submission received: 25 November 2021 / Revised: 7 January 2022 / Accepted: 28 January 2022 / Published: 30 January 2022
(This article belongs to the Special Issue Recent Progress on Advanced Foundation Engineering)

Abstract

:
Scouring can reduce the strength and rigidity of the pile–soil system and become one of the major causes for the failure of the structure of the partially-embedded single pile. The stratification of the soil fields has a significant influence on the internal force and deformation of laterally-loaded piles. A dynamic model of the laterally-loaded single pile in layered soil is established employing Hamilton’s principle based on the modified Vlasov foundation model. Then, the finite difference method is used to obtain the numerical matrix containing the control equations of the single pile to achieve accurate modeling of the soil–structure interaction (SSI) system affected by scouring, so as to solve the natural frequencies of the single pile. Green’s function method obtains the analytical solution of the forced vibration of the single pile. The effects of scouring and the layered soil on the dynamic response of the single pile are studied by numerical calculation and parameter analysis. It is shown that the dynamic model of the partially-embedded single pile in layered soil based on the modified Vlasov foundation model can accurately predict the dynamic characteristics of pile foundation affected by scouring. As the scouring degree intensifies, the first-order natural frequencies of the single pile in layered soil decrease significantly. The subgrade reaction coefficient of each layer of soil in the modified Vlasov foundation model decreases, and the shear coefficient increases. The first-order natural frequencies of the single pile at each scour level increase with the increase in the thickness of the underlying soil. When the elastic modulus of the first layer of soil is increased by one time, the first-order natural frequencies of the single pile are increased by about 20%.

1. Introduction

As a common form of deep foundation, pile foundations are usually used in bridges, pile foundations of the abutment, and hydraulic structures. They are often subjected to horizontal loads such as wind, waves, traffic, and earthquakes. Existing studies have shown that scouring reduces the strength and stiffness of the pile–soil interaction system and becomes one of the significant causes of structural failure [1,2,3]. Lin et al. proposed a simplified method to evaluate the effect of scour-hole dimensions on the response of laterally-loaded piles in the sand [4]. Zhang et al. studied the effects of stress history, scour depth, scour width, and scour-hole slope angle on the response of laterally-loaded piles in soft clay [5].
Due to the requirement of lateral load-bearing of pile foundations, the dynamic research of laterally-loaded piles has become a hot issue [6]. Scholars at home and abroad have conducted extensive studies on the dynamic characteristics of single piles under horizontal loads. Their theoretical models can be classified into the Continuum model, finite element or boundary element model, and Winkler foundation beam model. Using the continuum model, Poulos analyzed the deflection of piles in clay under a quasi-static cyclic load [7]. Based on the coupling of finite element method and boundary element method, Millán et al. studied the dynamic response of pile support structure under harmonic excitation [8]. Using the linear elastic Winkler foundation model, Novak analyzed the dynamic response of a single pile in the frequency domain [9].
For the study of the dynamic response of the partially-embedded single pile, it is critical to simulate the effect of the soil–structure interaction accurately. The Vlasov foundation model not only overcomes the shortcomings of the Winkler model that cannot show the shear characteristics of the foundation, but is also developed from the elastic half-space model, reducing the complexity of the half-space model, and has a complete theoretical basis [10]. Vallabhan correlated the elastic foundation parameters with the attenuation parameters and proposed a modified Vlasov foundation model [11]. Based on the modified Vlasov foundation model, Das and Sargand obtained the equation of displacement of laterally-loaded piles under dynamic load [12]. Since most of the soil fields in the project are layered, the physical properties of the soil fields that change with the buried depth have a significant impact on the mechanical properties of the pile foundation. Using the Mindlin integral equation, Pise numerically solved the horizontally-loaded pile in the two-layer homogeneous isotropic elastic foundation [13]. Basu et al. analyzed the elastic solution of horizontally-loaded piles in layered soil [14].
The numerical method Is the simplest way to solve the dynamic response of pile foundations in complex sites, and the results are accurate [15,16]. For the evaluation of the seismic behavior of the foundation system, M. Zucca et al. have conducted different types of non-linear numerical analysis, taking into account the soil–structure interaction effects [17]. Using the finite difference method, Bao et al. analyzed the effect of the soil–structure interaction on the natural frequencies of partially-embedded single piles [18]. Compared with fully-embedded single piles, the boundary conditions of partially-embedded single piles in layered soil are more complicated. It is more difficult to use the method of undetermined coefficients to solve the function of dynamic response. The above-mentioned finite element method or boundary element method usually requires much work and often requires professional software and complex models. Green’s function method has a clear concept and can obtain an accurate solution of the closed-form of the system, which is especially important for determining the dynamic response of the pile foundation. Using Green’s function method, Abu-Hilal obtained the steady-state solution of the Euler–Bernoulli beam [19]. Foda and Abdul-Jabbar obtained an analytical solution for the forced vibration of a supported beam [20]. Liang et al. analyzed the dynamic impedance of pile groups in soft clay and the effect of scouring [21]. Therefore, the combination of the finite difference method and Green’s function method will effectively solve the dynamic response of the partially-embedded single pile affected by scouring in the layered soil.
In order to accurately reveal the effect of scouring on the dynamic response of the partially-embedded single pile in layered soil, a dynamic model of the laterally-loaded single pile in layered soil is established employing Hamilton’s principle based on the modified Vlasov foundation model. The numerical matrix containing the control equation of a partially-embedded single pile is obtained by the finite difference method. The soil–structure interaction system is modeled accurately by changing the type and number of units in the process of scouring. Thereby, we can obtain a more accurate natural frequency of the single pile under common boundary conditions (free-fixed piles and free-free piles) and analyze the effect of scouring on its dynamic characteristics. Green’s function method obtains the analytical solution of the forced vibration of the single pile. The effects of elastic modulus and the thickness of layered soil on the dynamic response of a single pile are analyzed.

2. Solving Equations of Motion

2.1. Dynamic Model

Figure 1 shows the analysis diagram of the partially-embedded single pile in the layered soil with the pile top subjected to the lateral dynamic force P 0 ( t ) and the bending moment M 0 ( t ) . The physical parameters of the pile are: length L , radius R , elastic modulus E P , the moment of inertia I P , which l 1 is the length of the non-embedded section pile. As shown in Figure 1a, taking the three-layer soil field as an example, the analysis area is divided into five layers according to the embedded depth: the first layer is the non-embedded section; the second, third, and fourth layers are the embedded section; the 5th layer is the soil section below the pile end; the distance from the bottom of each layer to the pile top is H i ( i = 1 , 2 , 3 , 4 , 5 ), and H 0 = 0 . Lame constants of the elastic foundation are λ s = E s v s ( 1 + v s ) ( 1 2 v s ) and G s = E s 2 ( 1 + v s ) . Taking the undeformed pile top as the origin O , the cylindrical coordinate system O r θ z is established. In Figure 1a, u r , u θ , w z are the radial, tangential, and vertical displacement components of the pile body displacement along with the cylindrical coordinate system with the central axis of the pile as the z-axis. θ is the angle between any point and the r axis. The adopted Vlasov model simplifies the basic equation of the isotropic linear elastic continuum by introducing the constraints of displacement and stress and reduces the complexity of the half-space model. For simplicity, modeling and analysis are based on the following assumptions: the pile is slender, vertical and linear elastic; the soil field is a uniform, isotropic and of linear elastic material; there is no slippage or detachment at the pile–soil interface; the displacement of the pile body outside the direction of external excitation is ignored.
The displacement components of the soil field around the partially-embedded single pile can ignore the vertical displacement. They can be written as:
u r ( r , θ , z , t ) = u ( z , t ) ϕ ( r ) cos θ
u θ ( r , θ , z , t ) = u ( z , t ) ϕ ( r ) sin θ
w z ( r , θ , z , t ) = 0
where u ( z , t ) is the lateral displacement of the pile, w z is the vertical displacement components of the pile body displacement, and ϕ ( r )  is the dimensionless function representing the variation of the soil displacement along the r direction.
Based on the assumption of displacements given in Equations (1)–(3) and the stress–strain relations, the potential energy of the pile and soil system is given by:
U = 1 2 R 0 [ π ( λ s + 3 G s ) u 2 ( d ϕ d r ) 2 + 2 π G s ( d u d z ) 2 ϕ 2 ] r   d r   d z + 1 2 0 L E P I P ( d 2 u d z ) 2 d z + 1 2 L π R 2 G s ( d u d z ) 2 d z
Correspondingly, the kinetic energy of the pile and soil system can be written as:
T = 1 2 0 L m P ( d u d t ) 2 d z + 1 2 L m s ( d u d t ) 2 d z + π 0 R ρ s ( d u d t ) 2 θ 2 r   d r   d z
where m s and m P are the mass of the soil under the pile and the mass per linear meter of the pile, respectively, ρ s is the density of the soil around the pile.
The sum of work of non-conservative forces such as external force and bending moment can be written as:
W = P 0 ( t ) u ( 0 , t ) + M 0 ( t ) d u ( 0 , t ) d z
Substituting Equations (4)–(6) into Hamilton’s principle formulation, we can obtain:
t 1 t 2 i = 1 4 H i 1 H i ( E P I P d 4 u i d z 4 + k i u i 2 t i d 2 u i d z 2 + m i d 2 u i d t 2 ) δ u i d z   dt t 1 t 2 L ( k 4 u 5 2 t 5 d 2 u 5 d z 2 + m 5 d 2 u 5 d t 2 ) δ u 5 d z   dt + t 1 t 2 i = 1 4 H i 1 H i R [ π ( λ s i + 3 G s i ) n 3 ( d 2 ϕ d r 2 r + d ϕ d r ) 2 π ( G s i n 2 ρ s i n 1 ) ϕ r ] δ ϕ d r   dt t 1 t 2 { [ E P I P d 3 u 1 d z 3 2 t 1 d u 1 d z P 0 ] δ u 1 | z = 0 + i = 1 3 [ ( E P I P d 3 u i d z 3 + 2 t i d u i d z ) δ u i + ( E P I P d 3 u i + 1 d z 3 2 t i + 1 d u i + 1 d z ) δ u i + 1 ] z = H i + [ ( E P I P d 3 u 4 d z 3 + 2 t 4 d u 4 d z ) δ u 4 2 t 5 d u 5 d z δ u 5 ] z = L + 2 t 5 d u 5 d z δ u 5 | z = + [ E P I P d 2 u 1 d z 2 + M 0 ( t ) ] δ d u 1 d z | z = 0 + i = 1 3 ( E P I P d 2 u i d z 2 δ d u i d z E P I P d 2 u i + 1 d z 2 δ d u i + 1 d z ) | z = H i + E P I P d 2 u 4 d z 2 δ d u 4 d z | z = L } d t + i = 1 5 2 π H i 1 H i R ρ s i ϕ 2 r d r   dz d u i d t δ u i | t 1 t 2 + L ( λ s i + 3 G s i ) n 3 dz d u 5 d z δ u 5 | t 1 t 2   = 0
where
t i = π G s i R ϕ 2 r d r ;   k i = π ( λ s i + 3 G s i ) R ( d ϕ d r ) 2 r d r ;   m i = m p + 2 π ρ s i R ϕ 2 r d r ; ( i = 1 , 2 , 3 , 4 ) t 5 = π G s 4 R ϕ 2 r d r + 1 2 π G s 4 R 2 ( L z < ) ; n 1 = 0 ( d u i d t ) 2 d z ;   n 2 = 0 ( d u i d z ) 2 d z ; n 3 = 0 u i 2 d z
Therefore, the motion equation of the partially-embedded single pile can be obtained
E P I P d 4 u i d z 4 + k i u i 2 t i d 2 u i d z 2 + m i d 2 u i d t 2 = 0 ( i = 1 , 2 , 3 , 4 )
where k 1 = 0 , 2 t 1 = 0 , m 1 = m p .
With boundary conditions:
E P I P d 3 u 1 d z 3 2 t 1 d u 1 d z P 0 = 0 ( z = 0 )
E P I P d 2 u 1 d z 2 + M 0 ( t ) = 0 ( z = 0 )
( E P I P d 3 u 4 d z 3 + 2 t 4 d u 4 d z ) 2 t 5 d u 5 d z = 0 ( z = L )
k 4 u 5 2 t 5 d 2 u 5 d z 2 + m 5 d 2 u 5 d t 2 = 0 ( L z < )
Collecting the coefficients of δ ϕ for R r in Equation (7), we can obtain:
d 2 ϕ d r 2 r + d ϕ d r ( γ R ) 2 ϕ r = 0
where
( γ R ) 2 = 2 ( G s i n 2 ρ s i n 1 ) ( λ s i + 3 G s i ) n 3 = 2 [ G s i 0 ( d u i d z ) 2 d z ρ s i 0 ( d u i d t ) 2 d z ] ( λ s i + 3 G s i ) 0 u i 2 d z
where i = 1 , 2 , 3 , 4 .
The solution to Equation (13) is:
ϕ ( r ) = K 0 ( r γ R ) K 0 ( γ )
which satisfies the boundary conditions at r = R and r : ϕ ( R ) = 1 and ϕ ( ) = 0 . K 0 (   ) denotes the modified Bessel function of the second kind of order zero.
It is assumed that the pile is undergoing a steady-state harmonic motion. So let
P 0 ( t ) = P 0 e i ω t
M 0 ( t ) = M 0 e i ω t
u ( z , t ) = { u p ( z ) e i ω t , 0 z L u s ( z ) e i ω t , L < z <
where P 0 and M 0 are the amplitudes of lateral load and bending moment; ω is the circular frequency; u p ( z ) and u s ( z ) are the configuration functions of the displacement of the pile and the soil column below the pile end, respectively.
Substituting Equations (16)–(18) into Equation (8) gives:
E P I P d 4 u p i d z 4 2 t i d 2 u p i d z 2 + ( k i m i ω 2 ) u p i = 0 ( i = 1 , 2 , 3 , 4 )
with boundary conditions at z = 0 :
E P I P d 3 u p 1 d z 3 2 t 1 d u p 1 d z P 0 = 0 ( z = 0 )
E P I P d 2 u p 1 d z 2 + M 0 ( t ) = 0 ( free   pile   head )
with boundary conditions at z = H i   ( i = 1 , 2 , 3 ) according to continuity of deformation:
u p i u p i + 1 = 0
d u p i d z d u p i + 1 d z = 0
d 2 u p i d z 2 d 2 u p i + 1 d z 2 = 0
d 3 u p i d z 3 2 t i d u p i d z ( d 3 u p i + 1 d z 3 2 t i + 1 d u p i + 1 d z ) = 0
For free-fixed pile, boundary conditions at z = L are:
u p 4 = 0
d u p 4 d z = 0
u s ( L ) = u p ( L )
For free-free pile, boundary conditions at z = L are:
d 2 u p 4 d z 2 = 0
d 3 u p 4 d z 3 = 0
u s ( L ) = u p ( L )
Solving the differential equation given in Equation (19), with the boundary conditions u s ( ) = 0 and u s ( L ) = u p ( L ) , we get the following solution:
u s = u p ( L ) e a ( z L )
where α = k 4 m 4 ω 2 2 t 5
Meanwhile, the attenuation parameter γ can be written as:
( γ R ) 2 = 2 i = 2 4 [ G s i H i 1 H i ( d u p i d z ) 2 d z + ρ s i ω 2 H i 1 H i u p i 2 d z ] + N i = 2 4 ( λ s i + 3 G s i ) H i 1 H i u p i 2 d z + D
where N = ( G s 4 α + ρ s 4 ω 2 α ) u p 4 2 ( L ) , D = ( λ s 4 + 3 G s 4 2 α ) u p 4 2 ( L )

2.2. Finite Difference Method

By ignoring the effect of external load and damping, Equation (19) can be rewritten as:
u 1 + ω 2 ( m 1 E p I p u 1 ) = 0 , 0 z l 1
u i 2 t i E p I p u i + k i E p I p u i + ω 2 ( m i E p I p u i ) = 0 , ( i = 2 , 3 , 4 ) ,   l 1 < z L
Taking free-free pile as an example, the boundary conditions are:
Pile   top { u 1 ( 0 ) = 0 u 1 ( 0 ) = 0
Continuous { u i ( H i ) = u i + 1 ( H i ) u i ( H i ) = u i + 1 ( H i ) u i ( H i ) = u i + 1 ( H i ) u i ( H i ) = u i + 1 ( H i )
Pile   end { u 4 ( L ) = 0 u 4 ( L ) = 0 ( i = 1 , 2 , 3 )
The discretization of Equations (34) and (35) using the FDM is performed as [18]
δ z 4 u 1 , n + ω 2 ( m 1 E p I p u 1 , n ) = 0 , 0 z l 1
δ z 4 u i , n 2 t i E p I p δ z 2 u i , n + k i E p I p u i , n + ω 2 ( m i E p I p u i , n ) = 0 , ( i = 2 , 3 , 4 ) ,   l 1 < z L
In Equations (39) and (40), the subscript i in u i , n represents the displacements u 1 and u 2 , u 3 , u 4 of the non-embedded section and the embedded section of the pile and the subscript n represents the number of nodes in this section and δ z is the discretization expression of the mode shape function, which is derived using the Taylor series expansion [22].
In the range of l 1 < z L , the second and fourth derivatives of u i ( i = 2 , 3 , 4 ) can be written as:
u i u i , n 2 4 u i , n 1 + 6 u i , n 4 u i , n + 1 + u i , n + 2 Δ z 4
u i u i , n 1 2 u i , n + u i , n + 1 Δ z 2
Considering the virtual nodes u 4 , b + 1 and u 4 , b + 2 , from the boundary condition (36)–(38) we have
z = L : u 4 , b = u 4 , b 1 2 u 4 , b + u 4 , b + 1 Δ z 2 = 0 ,
u 4 , b = u 4 , b + 2 3 u 4 , b + 1 + 3 u 4 , b u 4 , b 1 Δ z 3 = 0
Substituting Equations (41)–(44) into Equation (40), we can obtain
( [ A ] + ω 2 [ B ] ) { U i } = 0
In Equation (45), A and B are diagonal matrices and { U i } = [ u i , a + 1 , u i , a + 2 , u i , b 1 , u i , b ] T ( i = 2 , 3 , 4 ) , where the subscript a is the number of nodes in the non-embedded section, the subscript b is the total number of nodes in a pile. Substituting Equations (36)–(38) into Equation (45), the matrix formulation of Equation (45) can be written as:
( [ B B i N N i 1 0 0 N N i B B i N N i 1 0 1 N N i B B i N N i 0 B B i N N i 1 0 0 1 N N i B B i N N i 0 0 0 4 / 3 N N i 2 B B i + 6 ] + ω 2 [ D D i 0 0 0 0 0 D D i 0 0 0 0 0 D D i 0 0 0 0 0 D D i 0 0 0 0 0 0 0 D D i ] ) [ u 2 , a + 1 u 2 , a + 2 u 2 , a + 3 u 2 , b 1 u 2 , b ] = 0
where B B i = 6 + k i Δ z 4 / E I + 4 t i Δ z 2 / E p I p , N N i = 4 2 t i Δ z 2 / E p I p ( i = 2 , 3 , 4 ) .
From Equations (36)–(38), we can obtain:
z = 0 : u 1 = u 1 , 0 2 u 1 , 1 + u 1 , 2 Δ z 2 = 0
u 1 = u 1 , 2 3 u 1 , 1 + 3 u 1 , 0 u 1 , 1 Δ z 3 = 0
and so
( [ C ] + ω 2 [ D ] ) { U 1 } = 0
where C and D are diagonal matrices and { U 1 } = [ u 1 , 1 , u 1 , 2 , u a 1 , u a ] T . In the same way, substituting boundary condition Equations (36)–(38) into Equation (49), the matrix formulation of Equation (49) can be obtained as:
( [ 1 2 1 0 0 2 5 4 1 0 1 4 6 4 0 0 1 4 6 1 4 0 0 0 1 4 6 ] + ω 2 [ D D 1 0 0 0 0 0 D D 1 0 0 0 0 0 D D 1 0 0 0 0 0 D D 1 0 0 0 0 0 0 0 D D 1 ] ) [ u 1 , 1 u 1 , 2 u 1 , 3 u 1 , a 1 u 1 , a ] = 0
where D D 1 = m 1 Δ z 4 / E p I p .
Considering the continuity condition in Equations (36)–(38), the diagonal matrix of the system is given by
( [ E ] + ω 2 [ F ] ) { U } = ( [ E ] + ω 2 [ F ] ) ( u 1 , 1 u 1 , 2 u 1 , 3 u 1 , a 1 u 1 , a u i , a + 1 u i , a + 2 u i , a + 3 u i , b 1 u i , b ) = 0
where { U } = [ U 1 , U i ] T = [ u 1 , 1 , u 1 , 2 u a 1 , u a , u i , a + 1 , u i , a + 2 , u i , b 1 , u i , b ] T   ( i = 2 , 3 , 4 )
E = ( 1 2 1 0 0 0 2 5 4 1 0 1 4 6 4 0 0 1 4 6 1 0 4 0 1 0 0 0 0 1 N N i B B i 0 N N i 1 0 0 0 0 1 N N i 0 B B i N N i 1 0 0 0 0 1 0 N N i B B i N N i 1 0 1 N N i B B i N N i 1 0 0 1 N N i B B i N N i 0 0 0 0 4 / 3 N N i 2 B B i + 6 ) F = ( D D i 0 0 0 0 0 0 D D i 0 0 0 0 D D i 0 0 0 D D i 0 0 0 0 0 0 0 0 D D i 0 0 0 0 0 0 0 D D i 0 0 0 D D i 0 0 D D i 0 0 D D i 0 0 0 0 0 D D i )
Therefore, we can obtain the natural frequencies of each order of partially-embedded single piles by letting the matrix determinant Equation (51) be zero.

3. Green’s Functions for Dynamic Response

3.1. Embedded Section

We can assume that there is a simple harmonic lateral load at the pile top, which can be written as:
P 0 ( t ) = δ ( z z 0 ) P 0 e i Ω t
where δ (   ) is the Dirac function, Ω is the frequency of the external excitation. We can let Ω = ω when calculating the Green’s functions for the dynamic response of the pile. For the pile top, z 0 = 0 .
Substitution of Equation (52) into Equation (35) gives:
u i + a 1 u i + a 2 u i = b P 0 δ ( z z 0 ) , ( i = 2 , 3 , 4 ) ,   l 1 < z L
where a 1 = 2 t i E p I p , a 2 = k i m i Ω 2 E p I p , b = 1 E p I p .
Using the Green’s function method, the solution of Equation (53) can be written as:
u i ( z ) = l 1 L P 0 ( z 0 ) G ( z ; z 0 ) d z 0
where P 0 ( z 0 ) is the function of the external excitation, G ( z ; z 0 ) is the Green’s function to be determined which is the solution of the following formulation:
u i + a 1 u i + a 2 u i = b δ ( z z 0 )
From Equation (55), it is seen that G ( z ; z 0 ) = U ( z ; z 0 ) . Laplace transformation is applied to the variable z in Equation (55), we can obtain:
u ^ i ( s ; z 0 ) = 1 s 4 + a 1 s 2 + a 2 [ b e s z 0 + ( s 3 + a 1 s ) u i ( 0 ) + ( s 2 + a 1 ) u i ( 0 ) + s u i ( 0 ) + u i ( 0 ) ]
where s = σ + i τ , u i ( 0 ) , u i ( 0 ) , u i ( 0 ) and u i ( 0 ) are constants which can be determined by the boundary conditions. To obtain the inverse transform of u ^ i ( s ; z 0 ) , we assume that s 4 + a 1 s 2 + a 2 = i = 1 4 ( s s i ) . We can get the results of the inverse transform by referring to the literature [23]:
L 1 [ b e s z 0 ( s s 1 ) ( s s 2 ) ( s s 3 ) ( s s 4 ) ] = H ( z z 0 ) [ A 1 ( z z 0 ) b + A 2 ( z z 0 ) b + A 3 ( z z 0 ) b + A 4 ( z z 0 ) b ]
L 1 [ s 3 + a 1 s ( s s 1 ) ( s s 2 ) ( s s 3 ) ( s s 4 ) ] = A 1 ( z ) ( s 1 3 + a 1 s 1 ) + A 2 ( z ) ( s 2 3 + a 1 s 2 ) + A 3 ( z ) ( s 3 3 + a 1 s 3 ) + A 4 ( z ) ( s 4 3 + a 1 s 4 )
L 1 [ s 2 + a 1 ( s s 1 ) ( s s 2 ) ( s s 3 ) ( s s 4 ) ] = A 1 ( z ) ( s 1 2 + a 1 ) + A 2 ( z ) ( s 2 2 + a 1 ) + A 3 ( z ) ( s 3 2 + a 1 ) + A 4 ( z ) ( s 4 2 + a 1 )
L 1 [ s ( s s 1 ) ( s s 2 ) ( s s 3 ) ( s s 4 ) ] = A 1 ( z ) s 1 + A 2 ( z ) s 2 + A 3 ( z ) s 3 + A 4 ( z ) s 4
L 1 [ 1 ( s s 1 ) ( s s 2 ) ( s s 3 ) ( s s 4 ) ] = A 1 ( z ) + A 2 ( z ) + A 3 ( z ) + A 4 ( z )
where H ( ) is the Heaviside function, and A i ( z )   ( i = 1 , 2 , , 4 ) are functions specified by
A 1 ( z ) = e s 1 z ( s 1 s 2 ) ( s 1 s 3 ) ( s 1 s 4 )
A 2 ( z ) = e s 2 z ( s 2 s 1 ) ( s 2 s 3 ) ( s 2 s 4 )
A 3 ( z ) = e s 3 z ( s 3 s 1 ) ( s 3 s 2 ) ( s 3 s 4 )
A 4 ( z ) = e s 4 z ( s 4 s 1 ) ( s 4 s 2 ) ( s 4 s 3 )
From Equation (56), we can obtain:
G ( z ; z 0 ) = L 1 ( b s 4 + a 1 s 2 + a 2 e s z 0 ) + L 1 ( s 3 + a 1 s s 4 + a 1 s 2 + a 2 ) u i ( 0 ) + L 1 ( s 2 + a 1 s 4 + a 1 s 2 + a 2 ) u i ( 0 ) + L 1 ( s s 4 + a 1 s 2 + a 2 ) u i ( 0 ) + L 1 ( 1 s 4 + a 1 s 2 + a 2 ) u i ( 0 )
Substitution of Equations (62)–(65) into Equation (66) gives the undetermined Green function G ( z ; z 0 )
G ( z ; z 0 ) = H ( z z 0 ) ϕ 1 ( z z 0 ) + ϕ 2 ( z ) u i ( 0 ) + ϕ 3 ( z ) u i ( 0 ) + ϕ 4 ( z ) u i ( 0 ) + ϕ 5 ( z ) u i ( 0 )
where ϕ i ( z )   ( i = 1 , 2 , , 5 ) are functions specified by
ϕ 1 ( z ) = i = 1 4 A i ( z ) b
ϕ 2 ( z ) = i = 1 4 A i ( z ) ( s i 3 + a 1 s i )
ϕ 3 ( z ) = i = 1 4 A i ( z ) ( s i 2 + a 1 )
ϕ 4 ( z ) = i = 1 4 A i ( z ) s i
ϕ 5 ( z ) = i = 1 4 A i ( z )
From Equations (62)–(65) and (68)–(72), we can obtain:
ϕ 1 ( k ) ( z ) = i = 1 4 s i k A i ( z ) b
ϕ 2 ( k ) ( z ) = i = 1 4 s i k A i ( z ) b
ϕ 3 ( k ) ( z ) = i = 1 4 s i k A i ( z ) ( s i 2 + a 1 )
ϕ 4 ( k ) ( z ) = i = 1 4 s i k + 1 A i ( z )
ϕ 5 ( k ) ( z ) = i = 1 4 s i k A i ( z )
and so
u i ( z ; z 0 ) = H ( z z 0 ) ϕ 1 ( z z 0 ) + ϕ 2 ( z ) u i ( 0 ) + ϕ 3 ( z ) u i ( 0 ) + ϕ 4 ( z ) u i ( 0 ) + ϕ 5 ( z ) u i ( 0 )
u i ( z ; z 0 ) = ϕ 1 ( z z 0 ) + ϕ 2 ( z ) u i ( 0 ) + ϕ 3 ( z ) u i ( 0 ) + ϕ 4 ( z ) u i ( 0 ) + ϕ 5 ( z ) u i ( 0 )
u i ( z ; z 0 ) = ϕ 1 ( z z 0 ) + ϕ 2 ( z ) u i ( 0 ) + ϕ 3 ( z ) u i ( 0 ) + ϕ 4 ( z ) u i ( 0 ) + ϕ 5 ( z ) u i ( 0 )
u i ( z ; z 0 ) = ϕ 1 ( z z 0 ) + ϕ 2 ( z ) u i ( 0 ) + ϕ 3 ( z ) u i ( 0 ) + ϕ 4 ( z ) u i ( 0 ) + ϕ 5 ( z ) u i ( 0 )
Substitution of Equation (67) into Equation (54) gives the response function of the pile
u i ( z ) = 0 L P 0 ( z 0 ) G ( z ; z 0 ) d z 0 = 0 L P 0 δ ( z z 0 ) [ H ( z z 0 ) ϕ 1 ( z z 0 ) + ϕ 2 ( z ) u i ( 0 ) + ϕ 3 ( z ) u i ( 0 ) + ϕ 4 ( z ) u i ( 0 ) + ϕ 5 ( z ) u i ( 0 ) ] d z 0 = P 0 [ H ( z z 0 ) ϕ 1 ( z z 0 ) + ϕ 2 ( z ) u i ( 0 ) + ϕ 3 ( z ) u i ( 0 ) + ϕ 4 ( z ) u i ( 0 ) + ϕ 5 ( z ) u i ( 0 ) ]
Correspondingly, the dynamic response function of the embedded pile is given by
u i ( z , t ) = P 0 [ H ( z z 0 ) ϕ 1 ( z z 0 ) + ϕ 2 ( z ) u i ( 0 ) + ϕ 3 ( z ) u i ( 0 ) + ϕ 4 ( z ) u i ( 0 ) + ϕ 5 ( z ) u i ( 0 ) ] cos Ω t ,   ( i = 2 , 3 , 4 )

3.2. Non-Embedded Section

As described in Section 3.1, Equation (34) can be written as:
u 1 + κ u 1 = b P 0 δ ( z z 0 )
where κ = m 1 Ω 2 E p I p .
Obviously, the solution of Equation (84) is given by
u 1 ( z ) = 0 L P 0 ( z 0 ) G ( z ; z 0 ) d z 0
Applying the Laplace transform method for the variable z in Equation (84), we can obtain:
u ^ 1 ( r ; z 0 ) = 1 r 4 + κ [ b e r z 0 + r 3 u 1 ( 0 ) + r 2 u 1 ( 0 ) + r u 1 ( 0 ) + u 1 ( 0 ) ]
where κ = m 1 Ω 2 E p I p , r = σ + i τ , and u 1 ( 0 ) , u 1 ( 0 ) , u 1 ( 0 ) , u 1 ( 0 ) are constants which can be determined by the boundary conditions of the pile.
To obtain the inverse transform of u ^ 1 ( r ; z 0 ) , we assume that r 4 + κ = i = 1 4 ( r r i ) .
G ( z ; z 0 ) = L 1 ( b r 4 + κ e r z 0 ) + L 1 ( r 3 r 4 + κ ) u 1 ( 0 ) + L 1 ( r 2 r 4 + κ ) u 1 ( 0 ) + L 1 ( r r 4 + κ ) u 1 ( 0 ) + L 1 ( 1 r 4 + κ ) u 1 ( 0 )
From Equation (87), we can obtain:
G ( z ; z 0 ) = H ( z z 0 ) ψ 1 ( z z 0 ) + ψ 2 ( z ) u 1 ( 0 ) + ψ 3 ( z ) u 1 ( 0 ) + ψ 4 ( z ) u 1 ( 0 ) + ψ 5 ( z ) u 1 ( 0 )
where ψ i ( z ) ( i = 1 , 2 , , 5 ) are functions specified by
ψ 1 ( z ) = i = 1 4 B i ( z ) b
ψ 2 ( z ) = i = 1 4 B i ( z ) r i 3
ψ 3 ( z ) = i = 1 4 B i ( z ) r i 2
ψ 4 ( z ) = i = 1 4 B i ( z ) r i
ψ 5 ( z ) = i = 1 4 B i ( z )
Therefore, the response function of the non-embedded section of the pile is as follows:
u 1 ( z ) = 0 L P 0 ( z 0 ) G ( z ; z 0 ) d z 0 = 0 L P 0 δ ( z z 0 ) [ H ( z z 0 ) ψ 1 ( z z 0 ) + ψ 2 ( z ) u 1 ( 0 ) + ψ 3 ( z ) u 1 ( 0 ) + ψ 4 ( z ) u 1 ( 0 ) + ψ 5 ( z ) u 1 ( 0 ) ] d z 0 = P 0 [ H ( z z 0 ) ψ 1 ( z z 0 ) + ψ 2 ( z ) u 1 ( 0 ) + ψ 3 ( z ) u 1 ( 0 ) + ψ 4 ( z ) u 1 ( 0 ) + ψ 5 ( z ) u 1 ( 0 ) ]
The dynamic response function of the pile of the non-embedded section is given by
u 1 ( z , t ) = P 0 [ H ( z z 0 ) ψ 1 ( z z 0 ) + ψ 2 ( z ) u 1 ( 0 ) + ψ 3 ( z ) u 1 ( 0 ) + ψ 4 ( z ) u 1 ( 0 ) + ψ 5 ( z ) u 1 ( 0 ) ] cos Ω t

4. Numerical Calculations

4.1. Solution Process

According to the modified Vlasov foundation model, the lateral displacement of a partially-embedded single pile in a layered soil field can be obtained by iteratively solving attenuation parameters. The process is shown in Figure 2.

4.2. Validity of the Present Solutions

In order to verify the validity of the finite difference method, the calculation model of Prendergast et al. and the model and calculation method of this paper were used to solve the variation of the first-order natural frequencies of the steel pipe pile with a length of 8.76 m under scouring conditions. The results are compared with the results measured on-site [16].
As shown in Figure 3, when the pile (free-free pile or fixed-free pile) is analyzed under the scour conditions, the scour effect is considered by completely removing the entire soil layer to a certain depth. After scouring, part of the pile body is exposed to the soil surface, and the pile is treated as partially-embedded. In this paper, the pile is divided into 40 pile units. There is no pile–soil interaction above the scour depth so they can be regarded as beam elements. Moreover, they can be regarded as spring-beam elements because of the pile–soil interaction below the scour depth. In the process of numerical calculation, the initial pile’s length of the non-embedded section is 2.19 m (scour level is 0), which can be represented as 10 beam elements and 30 spring-beam elements. The scouring process is simulated by removing spring supports from top to bottom.
The physical parameters of the pile and the foundation are shown in Table 1. The elastic modulus of the foundation is [18]:
E s = ρ s ν c 2 ( 1 + ν s ) ( 1 2 ν s ) 1 ν s
where ρ s and ν s are the soil’s density and Poisson’s ratio, respectively, and ν c is the velocity of the compression wave (213 m/s for sand [18]).
The subgrade reaction coefficient of the model of Prendergast et al. is [24]:
k 0 = 1.0 E s 1 ν s 2 [ E s D p 4 E p I p ] 1 12
where D p is the outer diameter of the pile, and E p I p is the flexural rigidity of the pile.
Figure 4 shows the comparison of the first-order natural frequencies of the partially-embedded single pile under scouring obtained by the model of Prendergast et al., the method in this paper, and the actual measured dates on-site. Referring to Equation (97), the subgrade reaction force in the calculation model of Prendert et al. has nothing to do with the pile diameter. The first-order natural frequencies of the single pile obtained by Prendergast et al.’s calculation model have the same trend as the actual measured dates on-site. The first-order natural frequencies of the single pile calculated by the method in this paper are closer to the measured dates on-site, and they are highly consistent with the measured dates on-site within the range of common scour depth (1–3 m). The results of the calculation model in this paper can more effectively predict the dynamic characteristics of partially-embedded single piles under the scouring action.
The displacement curve of the single pile with pile diameter d = 0.3   m , slenderness ratio λ = 10 , dimensionless frequency a 0 = 0.3 is calculated by using the proposed model and the model presented by the literature [25], which assumes that the pile has a fixed end and a pile top with a constrained angle. The calculation results are compared with the results given by the literature [25] and demonstrated in Figure 5. It can be seen from Figure 5 that the calculation results in this paper are basically consistent with the literature [25], which verifies the validity of Green’s function method in this paper.

5. Parameter Analysis

5.1. Effect of Scour Depth

In this section, the effect of scour depth on the dynamic response of the partially-embedded single pile in layered soil is analyzed. Table 1 shows the physical parameters of the single pile and the layered soil required for the calculation. In the calculation, the partially-embedded single pile passes through the three-layer foundation. The geological condition from top to bottom is: the first layer is the loose sand layer ( E s 1 = 10   MPa ); the second layer is the medium sand layer ( E s 2 = 20   MPa ); the third layer is the tight sand layer ( E s 3 = 50   MPa ). For sand, the Poisson’s ratio is approximately 0.3. The lateral excitation amplitude of the pile top is P 0 = 2 000   N . The setting of the scouring levels is shown in Figure 3, with each scouring level separated by 1.095 m.
Table 2 shows the variation of Vlasov foundation parameters and dynamic characteristics of partially-embedded single piles with the scour levels under the boundary conditions of free-free piles and fixed-free piles. Table 2 shows that the subgrade reaction coefficient k i of each layer of soil decreases with the increase in the scour levels, while the shear coefficient 2 t i increases with the increase in the scour levels.
The frequencies of each order of the partially-embedded single pile decrease rapidly with the increase in scouring levels. The depth of the embedded soil decreases with the scour process, and the pile–soil interaction effect gradually weakens.
l 3 The effect of the degree of scouring on the dynamic response amplitudes of the partially-embedded single pile in the layered soil fields is investigated in Figure 6. As shown in Figure 6, the displacement of the pile top of the single pile increases with the increase in the scour levels. When the scour level reaches 10, that is, the length of the non-embedded section of the pile satisfies l 1 L / 2 ( L is the pile length), the response amplitudes of the pile top of the partially-embedded single pile in the layered soil fields were greater than 0.01 m, and the single pile has presented lateral instability under dynamic load [26]. In addition, as the degree of scouring intensifies, the inverted point of the pile body moves downward, and the displacement at the inverted point increases gradually. The response amplitude of the free-free pile top is slightly larger than that of the free-fixed pile, and the difference between them increases with the increase in scour depth.

5.2. Effect of Soil Thickness

In the process of parameter analysis, the steel pipe piles in Table 1 are still used. The setting of the scouring levels is shown in Figure 3. Partially-embedded single piles pass through the three-layer foundation soil, the elastic modulus is 10, 20, 50 MPa from top to bottom, and the Poisson’s ratio is 0.3. In order to study the effect of the thickness of the bottom layer (third layer) of the layered soil fields on the dynamic response of a single pile, the depth of the first layer of soil is kept constant, and the thickness of the third layer is changed by changing the ratio of the thickness of the second and third layers of soil, in order: h 1 = 2.19   m , h 2 = 3.28   m , h 3 = 1.10   m ; h 1 = 2.19   m , h 2 = 2.19   m , h 3 = 2.19   m ; h 1 = 2.19   m , h 2 = 1.10   m , h 3 = 3.28   m ( h i represents the thickness of the soil layer).
In order to more clearly explore the effect of the thickness of the underlying soil in the layered soil fields on the dynamic response of the partially-embedded single pile, Table 3 shows the effect of h3 on the first-order natural frequencies and amplitudes of the pile top of the pile under varied scour levels. It can be seen from Table 3 that the first-order natural frequencies of the single pile under each scour level increase with the increase in h3. In addition, as the scour degree intensifies, the natural frequencies decrease significantly, and the response amplitudes of the pile top increase significantly. When the scour level reaches 10, the lateral instability of a single pile appears.
The variations of the mode shapes of the partially-embedded single pile under varied scour levels with the change of the thickness of the underlying soil are investigated in Figure 7. It is shown that the response amplitudes of the pile top of the single pile and the displacement amplitudes at the inverted point of the pile body decrease with the increase in the thickness of the underlying soil, regardless of the degree of scour.

5.3. Effect of Elastic Modulus of Soil

This section analyzes the effect of the elastic modulus of the first layer of soil on the dynamic characteristics of the partially-embedded single pile. The first layer of soil around the pile is: mucky soil layer ( E s 1 = 7   MPa ), loose sand layer ( E s 1 = 10   MPa ), silty sand layer ( E s 1 = 15   MPa ) and medium sand layer ( E s 1 = 20   MPa ), successively. The second layer is the medium sand layer ( E s 2 = 20   MPa ). The third layer is the tight sand layer ( E s 3 = 50   MPa ). Table 4 shows the physical parameters of the single pile and layered soil field required for the calculation.
Table 5 shows the variation of the first-order natural frequencies of the partially-embedded single pile and the response amplitudes of the pile top under varied scour levels with the elastic modulus of the first layer of soil. When the scour level is 10, the first layer of soil is completely cleared due to scour, so the dynamic characteristics of the single pile are the same. It can be seen from Table 5 that under each scour level, the response amplitudes of the pile top decrease significantly with the increase in the parameter E s 1 , and the first-order natural frequencies of a single pile increase with the increase in the parameter E s 1 . For example: if the elastic modulus of the soil increases by about 0.5 times, one time, and two times, the first-order natural frequencies of the single pile increase by about 10%, 20%, and 30%, respectively. The subgrade reaction coefficients of the elastic foundation increase with the increase in the elastic modulus of the soil, and accordingly, the lateral constraint of the soil field around the pile is strengthened.
Figure 8 shows the variation of the mode shapes of the partially-embedded single pile under varied scour levels with the change of the elastic modulus of the first layer of soil. Overall, the effect of parameter E s 1 on the mode shapes of the pile body is significant. With the increase in the parameter E s 1 , the lateral displacement amplitudes of the pile top of the single pile under each scour level decrease, while the displacement amplitudes at the inverted point of the pile body increase referring to the detail drawn in Figure 8. Corresponding to the results shown in Table 5, with the increase in the elastic modulus of soil around the pile, the lateral constraint of soil around the pile on the pile foundation increases. The pile-soil interaction effect is enhanced and the lateral deformation of the partially-embedded single pile is reduced.

6. Conclusions

A dynamic model of the laterally-loaded single pile in layered soil is established employing Hamilton’s principle based on the modified Vlasov foundation model. The finite difference method and Green’s function method are combined to obtain the single pile’s first-order natural frequencies and dynamic response configuration. Then, numerical results and discussions are presented to explore the single pile’s dynamic characteristics and examine the effects of scouring, elastic modulus, and the thickness of layered soil on the dynamic response of the partially-embedded single pile in the layered soil fields. The conclusions are as follows:
  • As the scour degree intensifies, the subgrade reaction coefficient of each layer of soil in the modified Vlasov foundation model decreases, and the shear coefficient increases. Furthermore, the first-order natural frequencies of the single pile in layered soil decrease significantly.
  • When scouring to the length of the non-embedded section of the pile satisfies l 1 L / 2 ( L is the pile length), the partially-embedded single pile will demonstrate lateral instability under the action of dynamic load.
  • As the thickness of the underlying soil increases, the first-order natural frequencies of the single pile under varied scour levels increase. During the response, amplitudes of the pile top decrease.
  • As the elastic modulus of the first layer of soil increases, the first-order natural frequencies of the single pile increase, and the response amplitudes of the pile top decrease significantly, while the displacement amplitudes at the inverted point of the pile body increase. The elastic modulus of the first layer of soil is increased by about 0.5 times, 1 time and 2 times, and the first-order natural frequencies of the single pile increase by about 10%, 20%, and 30%.
The model in this paper has certain limitations, such as not considering the slip or separation phenomenon of the pile–soil interface and nonlinear dynamic analysis. The pile–soil dynamic response in most sites is particularly complex. As a weak part of the pile foundation, the pile–soil contact surface is prone to the damage of the soil around the pile and even the detachment of the contact surface, which leads to the lack of binding force of the soil and instability of the pile foundation. Therefore, in the future study of the dynamic response of pile foundations, it is necessary to explore the dynamic characteristics of the laterally-loaded piles, considering the weakening effect of the contact surface.

Author Contributions

Conceptualization, J.M.; methodology, J.M. and S.H.; software, S.H.; validation, S.H., X.G. and D.L.; formal analysis, J.M., Y.G. and Q.L.; investigation, X.G. and D.L.; resources, J.M.; data curation, S.H.; writing—original draft preparation, J.M. and S.H.; writing—review and editing, J.M., S.H. and Q.L.; funding acquisition, J.M. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the National Natural Science Foundation of China (grant number 11502072, 51878265 and 52178331), the Young-backbone Teacher Foundation of Colleges and Universities of Henan Province (grant number 2019GGJS076) and the Key R&D and Promotion Special Project in Henan Province—Science and Technology Research Project (grant number 212102310946).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable. Written informed consent was obtained from the patients to publish this paper.

Data Availability Statement

Not applicable.

Acknowledgments

The authors greatly appreciate the comments from the reviewers, whose comments helped to improve the quality of the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wang, L.H.; Liu, C.L.; Wang, L.J. A simplified method for assessing the dynamic behaviour of a pile in layered saturated soils that considers current scour. Comput. Geotech. 2020, 128, 103831. [Google Scholar] [CrossRef]
  2. Qi, W.G.; Gao, F.P.; Randolph, M.F.; Lehane, B.M. Scour effects on p–y curves for shallowly embedded piles in sand. Géotechnique 2016, 66, 648–660. [Google Scholar] [CrossRef] [Green Version]
  3. Liu, Y.J.; Liu, C. Effects of scour-hole dimensions on lateral behavior of piles in sands. Comput. Geotech. 2019, 111, 30–41. [Google Scholar]
  4. Lin, C.; Han, J.; Caroline, B.; Robert, L.P. Analysis of laterally loaded piles in sand considering scour hole dimensions. J. Geotech. Geoenvironmental Eng. 2014, 140, 04014024. [Google Scholar] [CrossRef]
  5. Zhang, H.; Chen, S.; Liang, F. Effects of scour-hole dimensions and soil stress history on the behavior of laterally loaded piles in soft clay under scour conditions. Comput. Geotech. 2017, 84, 198–209. [Google Scholar] [CrossRef]
  6. Ai, Z.Y.; Li, Z.X. Dynamic analysis of a laterally loaded pile in a transversely isotropic multilayered half-space. Eng. Anal. Bound. Elem. 2015, 54, 68–75. [Google Scholar] [CrossRef]
  7. Poulos, H.G. Single pile response to cyclic lateral load. J. Geotech. Eng. ASCE 1982, 108, 355–375. [Google Scholar] [CrossRef]
  8. Millán, M.; Domínguez, J. Simplified BEM/FEM model for dynamic analysis of structures on piles and pile groups in viscoelastic and poroelastic soil. Eng. Anal. Bound. Elem. 2009, 33, 25–34. [Google Scholar] [CrossRef]
  9. Novak, M. Dynamic stiffness and damping of piles. Can. Geotech. J. 1974, 11, 574–598. [Google Scholar] [CrossRef]
  10. Ma, J.J.; Liu, Q.J.; Zhao, Y.Y. Stability of beams on modified Vlasov foundation subjected to lateral loads acting on the ends. Chin. J. Geotech. Eng. 2008, 30, 850–854. [Google Scholar]
  11. Vallabhan, C.V.G.; Daloglu, A.T. Modified Vlasov model for beams on elastic foundations. J. Geotech. Geoenvironmental Eng. 1991, 117, 956–996. [Google Scholar] [CrossRef]
  12. Das, Y.C.; Sargand, S.M. Forced vibrations of laterally loaded piles. Int. J. Solids Struct. 1999, 36, 4975–4989. [Google Scholar] [CrossRef]
  13. Pise, P.J. Laterally loaded piles in a two-layer soil system. J. Geotech. Eng. Div. 1982, 108, 1177–1181. [Google Scholar] [CrossRef]
  14. Bash, D.; Salgado, R. Elastic analysis of laterally loaded pile in multi-layered soil. Geomech. Geoengin. Int. J. 2007, 2, 183–196. [Google Scholar] [CrossRef]
  15. Avaei, A.; Ghotbi, R.; Aryafar, M. Investigation of pile–soil interaction subjected to lateral loads in layered soils. Am. J. Eng. Appl. Sci. 2008, 1, 76–81. [Google Scholar] [CrossRef] [Green Version]
  16. Prendergast, L.J.; Hester, D.; Gavin, K.; O’Sullivan, J.J. An investigation of the changes in the natural frequency of a pile affected by scour. J. Sound Vib. 2013, 332, 6685–6702. [Google Scholar] [CrossRef] [Green Version]
  17. Zucca, M.; Franchi, A.; Crespi, P.; Longrarini, N.; Ronca, P. The new foundation system for the transept reconstruction of the basilica di collemaggio. In Proceedings of the International Masonry Society Conferences, Milan, Italy, 9–11 July 2018; pp. 2441–2450. [Google Scholar]
  18. Bao, T.; Liu, Z. Evaluation of Winkler model and Pasternak model for dynamic soil-structure interaction analysis of structures partially embedded in soils. Int. J. Geomech. 2020, 20, 4019167. [Google Scholar] [CrossRef]
  19. Abu-hilal, M. Forced vibration of Euler–Bernoulli beams by means of dynamic Green functions. J. Sound Vib. 2003, 267, 191–207. [Google Scholar] [CrossRef]
  20. Foda, M.A.; Abduljabbar, Z. A dynamic Green function formulation for the response of a beam structure to a moving mass. J. Sound Vib. 1998, 210, 295–306. [Google Scholar] [CrossRef]
  21. Liang, F.; Zhang, H.; Huang, M. Influence of flood-induced scour on dynamic impedances of pile groups considering the stress history of undrained soft clay. Soil Dyn. Earthq. Eng. 2017, 96, 76–88. [Google Scholar] [CrossRef]
  22. LeVeque, R. Finite Difference Methods for Ordinary and Partial Differential Equations: Steady-State and Time-Dependent Problems; Society for Industrial and Applied Mathematics: Philadelphia, PA, USA, 2007; pp. 3–10. [Google Scholar]
  23. Asmar, N.H. Partial Differential Equations with Fourier Series and Boundary Value Problems; Dover Publications: New York, NY, USA, 2016; pp. 306–309. [Google Scholar]
  24. Ashford, S.A.; Juirnarongrit, T. Evaluation of pile diameter effect on initial modulus of subgrade reaction. J. Geotech. Geoenvironmental Eng. 2003, 129, 234–242. [Google Scholar] [CrossRef]
  25. Makris, N.; Gazetas, G. Dynamic pile-soil-pile interaction. Part II: Lateral and seismic response. Earthq. Eng. Struct. Dyn. 1992, 21, 145–162. [Google Scholar] [CrossRef]
  26. China Academy of Building Research. Technical Code for Building Pile Foundations: JGJ 94–2008, 1st ed.; China Construction Industry Press: Beijing, China, 2008; pp. 58–59. [Google Scholar]
Figure 1. Schematic diagram of the partially-embedded single pile in layered soil: (a) Laterally-loaded single pile; (b) Analytical model of pile–soil system.
Figure 1. Schematic diagram of the partially-embedded single pile in layered soil: (a) Laterally-loaded single pile; (b) Analytical model of pile–soil system.
Applsci 12 01504 g001
Figure 2. Flow chart of calculation.
Figure 2. Flow chart of calculation.
Applsci 12 01504 g002
Figure 3. Numerical model of the finite difference method and setting of scour levels.
Figure 3. Numerical model of the finite difference method and setting of scour levels.
Applsci 12 01504 g003
Figure 4. Comparison of the first-order natural frequencies calculated by the finite difference method with the measurement result by Prendergast et al. [16,18].
Figure 4. Comparison of the first-order natural frequencies calculated by the finite difference method with the measurement result by Prendergast et al. [16,18].
Applsci 12 01504 g004
Figure 5. Comparison of pile displacement calculated by Green’s functions method and the literature [22,25].
Figure 5. Comparison of pile displacement calculated by Green’s functions method and the literature [22,25].
Applsci 12 01504 g005
Figure 6. Effect of the scour degree on the dynamic response of the single pile in layered foundation: (a) Free-fixed pile; (b) Free-free pile.
Figure 6. Effect of the scour degree on the dynamic response of the single pile in layered foundation: (a) Free-fixed pile; (b) Free-free pile.
Applsci 12 01504 g006
Figure 7. Effect of the thickness of soil on the dynamic response of the single pile under varied scour levels: (a) Free-fixed pile, scour level-5; (b) Free-fixed pile, scour level-10; (c) Free-free pile, scour level-5; (d) Free-free pile, scour level-10.
Figure 7. Effect of the thickness of soil on the dynamic response of the single pile under varied scour levels: (a) Free-fixed pile, scour level-5; (b) Free-fixed pile, scour level-10; (c) Free-free pile, scour level-5; (d) Free-free pile, scour level-10.
Applsci 12 01504 g007aApplsci 12 01504 g007b
Figure 8. Effect of the elastic modulus of soil on the dynamic response of the single pile under varied scour levels: (a) Free-fixed pile, scour level-0; (b) Free-fixed pile, scour level-5 (c) Free-free pile, scour level-0; (d) Free-free pile, scour level-5.
Figure 8. Effect of the elastic modulus of soil on the dynamic response of the single pile under varied scour levels: (a) Free-fixed pile, scour level-0; (b) Free-fixed pile, scour level-5 (c) Free-free pile, scour level-0; (d) Free-free pile, scour level-5.
Applsci 12 01504 g008
Table 1. Physical parameters of the single pile and foundation during the verification.
Table 1. Physical parameters of the single pile and foundation during the verification.
Physical Meaning, Symbols and UnitsNumerical Value
Length of the pile L [m]8.76
Outer diameter of the pile R [m]0.17
Inter diameter of the pile r [m]0.157
Density of the pile ρ [kg/m3]7.8 × 103
Elastic modulus of the pile E P [MPa]2.0 × 105
Poisson’s ratio of the pile ν 0.3
Poisson’s ratio of elastic foundation ν s 0.3
Density of elastic foundation ρ s [kg/m3]2 × 103
Table 2. Effect of the scour on Vlasov foundation parameters and dynamic characteristics of the single pile.
Table 2. Effect of the scour on Vlasov foundation parameters and dynamic characteristics of the single pile.
Type of Single PileScour LevelsDistribution of Soil LayersVlasov Foundation ParametersDynamic Characteristics of Pile
k i   [ Nm 2 ] 2 t i   [ Nm 2 ] Natural Frequency [Hz]Amplitude of Pile Top [m]
Free-fixed pile0 l 1 30,516,730931,85524.750.00270
l 2 61,033,4601,863,709
l 3 152,583,6514,659,273
5 l 1 24,721,8841,411,27715.130.00663
l 2 49,443,7692,822,553
l 3 123,609,4217,056,383
10 l 1 //10.870.01343
l 2 41,639,2433,771,832
l 3 104,098,1069,429,580
15 l 1 //7.720.02424
l 2 40,133,5343,987,401
l 3 100,333,8359,968,503
Free-free pile0 l 1 30,140,280956,03324.110.00330
l 2 60,280,5611,912,067
l 3 150,701,4024,780,167
5 l 1 24,439,5641,441,14914.630.00900
l 2 48,879,1282,882,297
l 3 122,197,8207,205,743
10 l 1 //10.470.02063
l 2 41,411,4223,803,738
l 3 103,528,5559,509,346
15 l 1 //7.340.04313
l 2 40,015,2094,004,815
l 3 100,038,02310,012,037
Table 3. Effect of the thickness of the underlying soil on dynamic characteristics of the single pile under varied scour levels.
Table 3. Effect of the thickness of the underlying soil on dynamic characteristics of the single pile under varied scour levels.
Type of Single Pile The   Thickness   of   the   Soil   Layer   h 3   [ m ] Scour LevelsDynamic Characteristics of Pile
Natural Frequency [Hz]Amplitude of Pile Top [m]
Free-fixed pile1.10024.710.002720
515.030.006851
1010.850.013457
2.19024.750.002701
515.130.006631
1010.870.013429
3.29024.830.002683
515.340.006608
1010.910.012976
Free-free pile1.10024.070.003322
514.5430.009181
1010.4190.020786
2.19024.110.003296
514.6280.009050
1010.4680.020628
3.29024.180.003271
514.6320.008998
1010.5030.019846
Table 4. Physical parameters of the single pile and foundation in the parameter analysis.
Table 4. Physical parameters of the single pile and foundation in the parameter analysis.
Physical Meaning, Symbols and UnitsNumerical Value
Elastic modulus of the first layer of soil E s 1 [MPa]7, 10, 15, 20
Elastic modulus of the second and third layers of soil E s 2 , E s 3 [MPa]20, 50
Poisson’s ratio of elastic foundation ν s 0.3
Thickness of the first, second and third soil layers h i [m] h 1 = h 2 = h 3 = 2.19
Table 5. Effect of the elastic modulus of soil on dynamic characteristics of the single pile under varied scour levels.
Table 5. Effect of the elastic modulus of soil on dynamic characteristics of the single pile under varied scour levels.
Type of Single PileElastic Modulus [MPa]Scour LevelsDynamic Characteristics of Pile
Natural Frequency [Hz]Amplitude of Pile Top [m]
Free-fixed pile7022.690.00311
514.490.00673
1010.870.01343
10024.750.00270
515.130.00663
1010.870.01343
15027.200.00230
515.980.00633
1010.870.01343
20028.970.00208
516.750.00597
1010.870.01343
Free-free pile7022.140.00379
513.990.00904
1010.770.02063
10024.110.00330
514.630.00900
1010.770.02063
15026.460.00283
515.460.00869
1010.770.02063
20028.140.00257
516.150.00828
1010.770.02063
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ma, J.; Han, S.; Gao, X.; Li, D.; Guo, Y.; Liu, Q. Dynamic Lateral Response of the Partially-Embedded Single Piles in Layered Soil. Appl. Sci. 2022, 12, 1504. https://doi.org/10.3390/app12031504

AMA Style

Ma J, Han S, Gao X, Li D, Guo Y, Liu Q. Dynamic Lateral Response of the Partially-Embedded Single Piles in Layered Soil. Applied Sciences. 2022; 12(3):1504. https://doi.org/10.3390/app12031504

Chicago/Turabian Style

Ma, Jianjun, Shujuan Han, Xiaojuan Gao, Da Li, Ying Guo, and Qijian Liu. 2022. "Dynamic Lateral Response of the Partially-Embedded Single Piles in Layered Soil" Applied Sciences 12, no. 3: 1504. https://doi.org/10.3390/app12031504

APA Style

Ma, J., Han, S., Gao, X., Li, D., Guo, Y., & Liu, Q. (2022). Dynamic Lateral Response of the Partially-Embedded Single Piles in Layered Soil. Applied Sciences, 12(3), 1504. https://doi.org/10.3390/app12031504

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop