Study on the Time Domain Semi Analytical Method for Horizontal Vibration of Pile in Saturated Clay

: Under the framework of Biot porous media theory, a fractional order Kelvin model is used to describe the rheological effects of soil skeletons, and a coupled vibration model of saturated clay and a pile foundation is constructed. The Laplace transform is used to derive the analytical solution of the control equation in the transformation domain, and then the time–domain solution is obtained through numerical inversion. By analyzing numerical examples, the displacement and internal force response of pile foundations under horizontal vibration loads, as well as the inﬂuence of parameters, are studied. The results show that the displacement and internal force response of pile foundation vibrations in saturated clay foundations have a delayed effect. The stronger the rheological properties of the foundation soil, the more obvious the delay, the lower the load frequency, and the more signiﬁcant the inﬂuence of the rheological properties on the delayed effect. The stronger the rheological properties of the soil, the smaller the displacement amplitude of the pile foundation vibration, and the higher the load frequency, the greater the decrease in displacement amplitude. The stronger the rheological properties of the soil, the smaller the positive bending moment of the pile body, while the negative bending moment increases. Both positive and negative shear forces increase, but the shear force at the top of the pile is not affected. Therefore, when designing pile foundations in saturated clay foundations, it is necessary to appropriately increase the pile foundation or increase the reinforcement to meet the shear resistance of the pile foundation. The results of this study can provide a valuable reference for geotechnical and seismic engineers in pile foundation design.


Introduction
Saturated clay is a typical weak soil that cannot be directly used as a bearing foundation for structures.Therefore, in practical engineering, pile foundations are usually used to pass through saturated clay layers to directly transmit loads to the actuating layer, while pile foundations are buried in saturated clay.During the service process, pile foundations not only need to bear the vertical loads transmitted by the upper structure but also bear various horizontal dynamic loads, such as environmental vibration, earthquakes, and impacts [1][2][3][4][5].Therefore, the horizontal vibration response of pile foundations has always been a key issue of concern in engineering.At present, the design and analysis methods for the bearing capacity of pile foundations are very mature, while dynamic analysis is still in the development stage.Amin [6] established a coupled three-dimensional dynamic analysis model for piles in a liquefied foundation, studied the dynamic characteristics of piles, and analyzed the reliability of the model under dynamic loads using centrifugal tests.J.M. [7] studied the influence of seismic body wave types and their incidence angles on the dynamic response of pile foundations.Nguyen [8] used numerical simulation methods to simulate the dynamic behavior of soil, pile foundations, and structures under seismic excitation while considering the seismic performance of buildings located on soft soil.Kaustav [9] used finite difference software to analyze the deflection and bending moment behavior along the pile length in non liquefied and liquefied soils under seismic loading conditions, and studied the horizontal and vertical inertial interactions as well as the kinematic interactions caused by free field motion.Kamel [10] used numerical simulation methods to study the response of adjacent pile foundations, strong grouting protection, and soil interactions induced by tunnel construction.It can be seen that further research on the horizontal vibration characteristics of pile foundations in saturated clay is of great significance for engineering construction, seismic resistance, and the safe service of saturated clay sites.
Based on the Biot porous medium theory, scholars, such as Jin [11], have conducted many studies on the horizontal vibration characteristics of pile foundations in saturated soil.According to the mixed boundary value conditions, the dual integral equations of the horizontal vibration of a rigid disk on a saturated poroelastic half-space are established, analytically examining the horizontal vibration of a rigid disk on a saturated poroelastic half-space.The dynamic response of pile groups partially buried in layered saturated soil under horizontal harmonic loads was studied [12].The Winkler model was used to derive the soil-pile dynamic interaction factor.Based on the dynamic interaction factor, the horizontal impedance of the pile group was obtained using the superposition principle, providing a valuable reference for geotechnical and seismic engineers in pile foundation design.Zhang [13] developed a simplified analytical solution for the impedance function of viscoelastic pile groups partially embedded in a layered, transversely isotropic, liquid-saturated, viscoelastic soil medium under a vertical time-harmonic load.Liang [14] provided a closed-form solution for the horizontal vibration of pipe piles in saturated soil, considering radial heterogeneity effects.Liu [15] made beneficial extensions to the horizontal vibration of piles in saturated soil by considering the variable cross-section and defects of the pile body.In addition, Hu [16] studied the horizontal vibration of partially buried piles for offshore structures.In the latest reported study, Zou [17] derived the dynamic impedance frequency domain analytical formula for horizontal vibrations of pile foundations in unsaturated soil using integral equation theory.Cui [18] studied the frequency domain response of horizontal vibrations of a single pile, considering the non-uniformity of the soil around the pile.Xu [19] conducted vibration table experiments to test the vibration response of pile foundations in clay foundations.The above research has made important contributions toward revealing the dynamic impedance and frequency domain characteristics of horizontal vibrations of pile foundations in saturated soil.However, it can also be seen that the above studies all use frequency domain analysis methods and obtain steady-state solutions without considering the viscous rheology of saturated soil.
In practical engineering, soil particles undergo rearrangement and skeleton dislocation under load, resulting in significant rheological properties (relaxation, creep) [20][21][22].Rheology leads to a time-dependent stress-strain relationship in the soil and also alters the time-domain response of pile foundation vibrations.To investigate the time-domain response of pile foundation horizontal vibrations under the influence of saturated clay rheological effects, this paper uses a fractional order Kelvin rheological model to describe the stress-strain relationship of the soil skeleton, establishes a coupled vibration model of saturated clay and pile foundation, and uses Laplace transform to solve the control equation to obtain a semi-analytical time-domain solution.Based on case analysis, the displacement and internal force response of pile foundations in a saturated clay foundation under horizontal vibration loads were studied, and the influence of soil physical parameters was discussed.

Fractional Order Kelvin Model
The classic Kelvin model is composed of a spring element K and a Newton dashpot N in parallel, as shown in Figure 1.The fractional Kelvin rheological model is obtained by replacing the Newton dashpot with the Abel dashpot in the classical Kelvin model [23], and the constitutive equation is described as: where σ and ε are stress and strain, as functions of the time coordinate t; E s is the elastic modulus of spring element K; and η and a are the viscosity coefficient and viscous order of Abel dashpot N, respectively.Equation (1) shows that the dashpot N will degenerate into a linear elastic solid in the case of a = 0 and degenerate into an ideal Newtonian fluid in the case of a = 1.In the general case of 0 < a < 1, dashpot N represents a viscoelastic solid, and the value of a represents the proportion of solid properties to fluid properties [24].

Fractional Order Kelvin Model
The classic Kelvin model is composed of a spring element K and a Newton dashpot N in parallel, as shown in Figure 1.The fractional Kelvin rheological model is obtained by replacing the Newton dashpot with the Abel dashpot in the classical Kelvin model [23], and the constitutive equation is described as: where σ and ε are stress and strain, as functions of the time coordinate t; Es is the elastic modulus of spring element K; and η and a are the viscosity coefficient and viscous order of Abel dashpot N, respectively.Equation (1) shows that the dashpot N will degenerate into a linear elastic solid in the case of a = 0 and degenerate into an ideal Newtonian fluid in the case of a = 1.In the general case of 0 < a < 1, dashpot N represents a viscoelastic solid, and the value of a represents the proportion of solid properties to fluid properties [24].The Laplace transform is introduced for a time-domain function f(t) as follows: where ( ) f s is the Laplace transform function corresponding to the function ( ) where σ and ε are stress and strain tensors, θ is volumetric strain, I is a three-order identity matrix, and λ and µ are the Lamé constants.According to the definition, the expression of λ and µ can be obtained as The Laplace transform is introduced for a time-domain function f (t) as follows: where f (s) is the Laplace transform function corresponding to the function f (t), s is the transform parameter of the time coordinate t, ζ is any real number, and i = √ −1 is the imaginary unit.
Using Equation (2) to perform the Laplace transform on Equation (1), we obtain Equation (3) reflects the stress-strain relationship in the Laplace transform domain with a complex modulus Ê = E s + ηs a .Equation (3) is expanded to a three-dimensional stress state and written in tensor form as where σ and ε are stress and strain tensors, θ is volumetric strain, I is a three-order identity matrix, and λ and µ are the Lamé constants.According to the definition, the expression of λ and µ can be obtained as where v is the Poisson's ratio.

Horizontal Coupled Vibration Model of Pile Soil
As shown in Figure 2, end-bearing piles are buried in a saturated clay foundation using a cylindrical coordinate system (r-θ-z) to build a pile-soil coupling system model.The bottom of the pile is a fixed end, with a pile length of H and a pile diameter of r 0 .The pile top is subjected to a sinusoidal horizontal vibration load P(t), which is defined as: ( ) ( ) where v is the Poisson's ratio.

Horizontal Coupled Vibration Model of Pile Soil
As shown in Figure 2, end-bearing piles are buried in a saturated clay foundation using a cylindrical coordinate system (r-θ-z) to build a pile-soil coupling system model.The bottom of the pile is a fixed end, with a pile length of H and a pile diameter of r0.The pile top is subjected to a sinusoidal horizontal vibration load P(t), which is defined as: In the formula, P0 is the load amplitude, ω is the load frequency.Introducing the Laplace transform of time t into Equation ( 6), the expression of cyclic load in the Laplace transform domain is obtained as: where ( ) is the Laplace transformation of P(t), and s is the transformation parameter.
For the problem of horizontal vibration, ignoring the vertical displacement of the soil, the motion equation of saturated clay in the cylindrical coordinate system is: where u is the displacement tensor of the soil skeleton, w is the relative displacement between the pore water and the soil skeleton, ρ denotes the bulk density of the porous material, and ρ = (1 − n)ρs + nρf.n is the porosity, and ρs and ρf are the density of the soil skeleton and pore water, respectively.In addition, kf is the permeability coefficient of the soil.γw = ρf g is the gravity density of the pore water, where g denotes the gravitational acceleration, ▽ represents the Hamiltonian operator, and in the cylin- In the formula, P 0 is the load amplitude, ω is the load frequency.Introducing the Laplace transform of time t into Equation ( 6), the expression of cyclic load in the Laplace transform domain is obtained as: where P(s) is the Laplace transformation of P(t), and s is the transformation parameter.
For the problem of horizontal vibration, ignoring the vertical displacement of the soil, the motion equation of saturated clay in the cylindrical coordinate system is: where u is the displacement tensor of the soil skeleton, w is the relative displacement between the pore water and the soil skeleton, ρ denotes the bulk density of the porous material, and ρ = (1 − n)ρ s + nρ f .n is the porosity, and ρ s and ρ f are the density of the soil skeleton and pore water, respectively.In addition, k f is the permeability coefficient of the soil.γ w = ρ f g is the gravity density of the pore water, where g denotes the gravitational acceleration, represents the Hamiltonian operator, and ∇ = ∂ ∂r + 1 r ∂ ∂ϕ + ∂ ∂z in the cylindrical coordinate system.σ a is the total stress tensor on the porous material unit, and p is the pore water pressure.Taking compressive stress as the positive direction, the principle of effective stress is expressed as where σ is the effective stress tensor of the soil skeleton, and α s is the Biot-Willis coefficient.
The seepage continuity equation is expressed as where M is a Biot coefficient to describe the compressibility of the pore water; the superscript "•" represents the first-order differential of time.
The fractional order Kelvin model is introduced into the soil skeleton; that is, the effective stress σ satisfies Equation (4).Substituting Equations ( 4) and (10) into Equation ( 8) and then performing Laplace transform, the following equation can be derived: where Due to the study of horizontal vibrations of pile foundations in this article, the pile top pressure and constraint stiffness can be ignored.Considering the influence of the pile shear effect, the Timoshenko beam model is used to establish the pile foundation motion equation as follows: where u p and ϕ p are the horizontal and angular displacements; E p , G p , I p , S, ρ P are the elastic modulus, shear modulus, interface moment of inertia, cross-sectional area, and density, respectively; and q is the load concentration of pile-soil interaction.
The initial condition for the pile-soil system is: Considering that the surface is free and completely permeable, the boundary condition of the soil on the surface is: The boundary condition of the soil at radial infinity is: The boundary condition of the soil on the top surface of the bedrock is: This article considers the problem of small deformation vibrations.Assuming that the pile-soil is in complete contact and the contact surface is impermeable, the boundary conditions of the pile-soil interface are: The boundary conditions for the free end of the pile top are: The boundary condition for the fixed end of the pile bottom is:

Solving the Control Equation of Saturated Clay Movement
By introducing potential functions φ s , ψ s , φ f , ψ f and Helmholtz decomposition of the soil skeleton displacement and water displacement, we obtain Substituting Equation (23) into Equations ( 11) and ( 12) and organizing it, we obtain: where Using the method of separating variables to solve Equation (24), Substituting Equation ( 25) into Equation ( 24) and using boundary condition Equations ( 16)- (18), we obtain the general solutions of T 1 (z) and T 2 (z) as follows [16]: In the equation, γ k = iπ(2k − 1)/(2H), k = 1, 2, 3; K 1 () represents the Bessel function of the second type of first-order imaginary argument, with A 1 ~A5 as undetermined coefficients; β 1 , β 2 , β 3 and ε 1 , ε 2 , ε 3 are all known functions related to s, and the specific form is shown in Appendix A. Substituting Equation (26) into Equation ( 25) yields: where ch() represents the hyperbolic cosine function, and B 1 , B 2 , and B 3 are the reorganized undetermined coefficients.
Substituting Equation (27) into Equation ( 23) yields the displacement solution of the soil skeleton and water in the Laplace transform domain in saturated soil: Using Equation ( 11), the pore pressure in the Laplace domain is obtained as: Substituting boundary condition Equations ( 16) and ( 20) into Equation ( 30), we obtain: where κ 1 , κ 2 are known coefficients, and the specific expression is shown in Appendix A.
According to the stress balance conditions at the pile-soil interface, there are: Substituting Equations ( 28)- (30) into Equation ( 32) and taking r = r 0 , the load concentration of the pile-soil interaction is obtained as:

Solution of Pile Foundation Vibration Equation
Performing Laplace transformation on Equations ( 13) and ( 14), we obtain: By using Equation (34), it can be obtained that: Taking the partial derivative of z on both sides of Equation (34) and then substituting Equation (35) to compile the decoupled pile foundation motion equation: The general solution of the fourth-order ordinary differential Equation ( 36) is [25]: where sh() represents a hyperbolic sine function, and D 1 ~D8 are undetermined coefficients; Q 1 and Q 2 are a set of special solutions to Equation (34).In Equations ( 36) and (37), the expressions for χ 1 , χ 2 , C 1 , C 2 , Q 1 , and Q 2 can be found in Appendix A.
From the relationship neutralized by Equation (37), it can be obtained that: where η 1 ~η4 are known coefficients, and the specific form is shown in Appendix A.
Using the pile-soil contact coupling condition Equation ( 19) and the pile end boundary condition Equations ( 21) and ( 22), five boundary equations were constructed to determine the undetermined coefficients B 1 and D 1 -D 4 .At this point, the displacement field, stress field, and pore pressure of the pile-soil system have been determined.Combining the initial condition Equation ( 15) for Laplace inversion, a time-domain numerical solution is obtained.To ensure the accuracy of the long-term inversion, the Crump method [26,27] is used for Laplace inversion in this paper.

Algorithm Verification
To achieve comparative verification, the soil model in this paper was degraded to linear elastic saturated soil, and the pile foundation model was degraded to the Euler beam.According to beam theory, when the shear modulus of Timoshenko beams is sufficiently large, the cross-section does not undergo shear deformation, and Timoshenko beams degenerate into Euler beams [28].Reference [29] used the Euler beam and linear elastic saturated soil models to provide the displacement of the pile top in saturated soil under a triangular impact load.The soil model is degraded into elastic saturated soil by taking the viscosity order a = 0, and the Timoshenko beam is degraded into the Euler beam by taking the pile foundation shear modulus G p = 10E p .The other calculation parameters are the same as those in reference [29] (see Table 1).Meanwhile, the load P(t) is defined as a triangular impact load, and its time-domain and Laplace transform domain expressions are: where t 0 is the duration of the load action, and H() represents the Heaviside step function.Substitute Equation (39) into the system equation and calculate the time-domain solution of the pile top displacement through MATLAB(2018) programming.Figure 3 shows a comparison of the calculation results of the pile top displacement, and the results show that the degenerate solution in this paper is in good agreement with the literature solution (dimensionless displacement u * p = u p /r 0 and dimensionless time t * = t/t 0 in the figure).From Figure 3, it can also be seen that when considering the viscous rheological behavior of the soil (a = 0.5), the peak displacement of the pile top significantly decreases, and the deformation recovery rate after loading slows down, which is consistent with the conclusion obtained by Yin [30]   From Figure 3, it can also be seen that when considering the viscous rheologic behavior of the soil (a = 0.5), the peak displacement of the pile top significantly decreas and the deformation recovery rate after loading slows down, which is consistent with t conclusion obtained by Yin [30] through clay consolidation experiments.To further verify the reliability of the proposed method, a comparison was ma between the proposed method and the finite element method for the horizontal vibrati model of pile foundations in fractional-order saturated clay.Based on COM SOL(6.0)multi-physicsfield analysis software, a finite element model of soil was co structed under the solid mechanics porous media seepage module (as shown in Figure The stress-strain relationship of the soil skeleton was described using the fractional ord Kelvin Voigt constitutive model, and an infinite element domain was used at the rad boundary to make the soil model consistent with the theoretical model in this pap Under the structural mechanics module, Timoshenko beams are used to simulate p foundations, and the pile-soil interface is a coupled contact.To further verify the reliability of the proposed method, a comparison was made between the proposed method and the finite element method for the horizontal vibration model of pile foundations in fractional-order saturated clay.Based on COMSOL(6.0)multiphysicsfield analysis software, a finite element model of soil was constructed under the solid mechanics porous media seepage module (as shown in Figure 4).The stress-strain relationship of the soil skeleton was described using the fractional order Kelvin Voigt constitutive model, and an infinite element domain was used at the radial boundary to make the soil model consistent with the theoretical model in this paper.Under the structural mechanics module, Timoshenko beams are used to simulate pile foundations, and the pile-soil interface is a coupled contact.Figure 5 shows the comparison between the solution proposed in this paper and the finite element solution.The figure shows that the results of the two can be well matched when the load frequency is small, and when the load frequency is high (ω = 4π), the finite element calculation results show instability.Due to the use of the time integration method in finite element calculations, a very small integration step size is required to ensure calculation accuracy when the frequency is high.However, a small step size will result in very low calculation efficiency (when the step size is on the order of 10 -4 , the calculation time is about 3 h).Compared to other methods, the semi analytical method Figure 5 shows the comparison between the solution proposed in this paper and the finite element solution.The figure shows that the results of the two can be well matched when the load frequency is small, and when the load frequency is high (ω = 4π), the finite element calculation results show instability.Due to the use of the time integration Appl.Sci.2024, 14, 778 10 of 16 method in finite element calculations, a very small integration step size is required to ensure calculation accuracy when the frequency is high.However, a small step size will result in very low calculation efficiency (when the step size is on the order of 10 −4 , the calculation time is about 3 h).Compared to other methods, the semi analytical method proposed in this paper has higher computational efficiency and stability.Figure 5 shows the comparison between the solution proposed in this paper an finite element solution.The figure shows that the results of the two can be well mat when the load frequency is small, and when the load frequency is high (ω = 4π), the element calculation results show instability.Due to the use of the time integr method in finite element calculations, a very small integration step size is requir ensure calculation accuracy when the frequency is high.However, a small step size result in very low calculation efficiency (when the step size is on the order of 10 - calculation time is about 3 h).Compared to other methods, the semi analytical me proposed in this paper has higher computational efficiency and stability.

Pile Foundation Vibration Analysis
Based on Table 1, add viscosity coefficient η = 1 × 10 8 Pa•s, and use load fun Equation ( 6) to calculate the horizontal vibration response of pile foundation in satu clay foundation.Figure 6 shows the time-domain response of pile top displacemen der different shear modulus Gp values.It can be seen from the figure that as the s modulus of the pile foundation increases, the amplitude of pile top displacemen creases.When Gp increases to 6Ep, the displacement amplitude no longer changes, cating that the Timoshenko beam has converged to the Euler beam. Figure 7 show vibration response of pile top displacement under different load frequencies.The fi shows that the larger the load frequency, the smaller the amplitude of pile top disp ment.This result is consistent with basic vibration theory [31,32].

Pile Foundation Vibration Analysis
Based on Table 1, add viscosity coefficient η = 1 × 10 8 Pa•s, and use load function Equation ( 6) to calculate the horizontal vibration response of pile foundation in saturated clay foundation.Figure 6 shows the time-domain response of pile top displacement under different shear modulus G p values.It can be seen from the figure that as the shear modulus of the pile foundation increases, the amplitude of pile top displacement decreases.When G p increases to 6E p , the displacement amplitude no longer changes, indicating that the Timoshenko beam has converged to the Euler beam. Figure 7 shows the vibration response of pile top displacement under different load frequencies.The figure shows that the larger the load frequency, the smaller the amplitude of pile top displacement.This result is consistent with basic vibration theory [31,32].Figure 8 reflects the influence of soil viscosity order on the response of pile top di placement.The figure shows that the larger the viscosity order, the smaller the amplitud of pile top displacement, and the peak point of displacement moves backward on th time axis, reflecting the delayed effect of soil creep.According to the definition of th Kelvin rheological model, the higher the viscosity order, the stronger the rheologic properties.Therefore, the results in Figure 6 indicate that the stronger the rheologic properties of the soil, the smaller the displacement response amplitude of the pile fou dation vibration and the more obvious the delay effect.Figure 9 further indicates that a the order of soil viscosity increases, the amplitude of pile top displacement decreas almost linearly, and the higher the load frequency, the greater the decreasing rate of di placement amplitude.Figure 8 reflects the influence of soil viscosity order on the response of pile top displacement.The figure shows that the larger the viscosity order, the smaller the amplitude of pile top displacement, and the peak point of displacement moves backward on the time axis, reflecting the delayed effect of soil creep.According to the definition of the Kelvin rheological model, the higher the viscosity order, the stronger the rheological properties.Therefore, the results in Figure 6 indicate that the stronger the rheological properties of the soil, the smaller the displacement response amplitude of the pile foundation vibration and the more obvious the delay effect.Figure 9 further indicates that as the order of soil viscosity increases, the amplitude of pile top displacement decreases almost linearly, and the higher the load frequency, the greater the decreasing rate of displacement amplitude.To further explore the parameter correlation of the displacement delay effect, this article defines delay time t a = t 1 − T/4, where t 1 is the time to the first displacement peak and T = 2π/ω is the load period.Figure 10 shows the variation of delay time ta with the viscosity order.It can be seen from the figure that as the viscosity order increases, the delay time correspondingly increases.It is worth noting that the lower the load frequency, the more significant the change in delay time with the viscosity order.To further explore the parameter correlation of the displacement delay effect, this article defines delay time T π ω = is the load period.Figure 10 shows the variation of delay time ta with the viscosity order.It can be seen from the figure that as the viscosity order increases, the delay time correspondingly increases.It is worth noting that the lower the load frequency, the more significant the change in delay time with the viscosity order.T π ω = is the load period.Figure 10 shows the variation delay time ta with the viscosity order.It can be seen from the figure that as the viscosi order increases, the delay time correspondingly increases.It is worth noting that th lower the load frequency, the more significant the change in delay time with the viscosi order.Figures 11 and 12 show the distribution of pile bending moment M p = ∂ 2 u p /∂z 2 and shear force Q p = ∂ 3 u p /∂z 3 with depth at t = T/4, respectively.Figure 11 shows that as the viscosity order increases, the positive bending moment of the pile decreases while the negative bending moment increases.Figure 12 shows that as the viscosity order increases, both the positive and negative shear forces will increase.However, it should be considered that due to the delay effect, the bending moment and shear force values at t = T/4 cannot represent the peak values at each viscosity order.Therefore, Figures 13-15 show the time-domain response of bending moment and shear force at the characteristic depth.Figure 13 shows that there is also a delay effect in the bending moment of the pile body, and as the viscosity order increases, the amplitude of the bending moment decreases slightly.Figure 14 shows that the shear force at the top of the pile is almost unaffected by the viscosity order, and there is no significant delay effect or amplitude change phenomenon.Figure 15 shows that at larger depths, there is a significant delay in shear force, and the larger the viscosity order, the greater the shear amplitude.
Figure 13 shows that there is also a delay effect in the bending moment of the pile bod and as the viscosity order increases, the amplitude of the bending moment decreas slightly.Figure 14 shows that the shear force at the top of the pile is almost unaffected the viscosity order, and there is no significant delay effect or amplitude change ph nomenon.Figure 15 shows that at larger depths, there is a significant delay in shear for and the larger the viscosity order, the greater the shear amplitude.Figure 13 shows that there is also a delay effect in the bending moment of the pile bod and as the viscosity order increases, the amplitude of the bending moment decrease slightly.Figure 14 shows that the shear force at the top of the pile is almost unaffected b the viscosity order, and there is no significant delay effect or amplitude change ph nomenon.Figure 15 shows that at larger depths, there is a significant delay in shear forc and the larger the viscosity order, the greater the shear amplitude.

Conclusions
This article uses the fractional order Kelvin model to describe the rheological prop erties of soil, considers small deformation to establish a coupled vibration model of sat urated clay and pile foundation, and uses Laplace transform to solve the control equation to obtain a semi analytical solution in the time domain.Based on case analysis, the hori

Conclusions
This article uses the fractional order Kelvin model to describe the rheological properties of soil, considers small deformation to establish a coupled vibration model of saturated clay and pile foundation, and uses Laplace transform to solve the control equation to obtain a semi analytical solution in the time domain.Based on case analysis, the hori-

Conclusions
This article uses the fractional order Kelvin model to describe the rheological properties of soil, considers small deformation to establish a coupled vibration model of saturated clay and pile foundation, and uses Laplace transform to solve the control equation to obtain a semi analytical solution in the time domain.Based on case analysis, the horizontal vibration characteristics of pile foundations in a saturated clay foundation were studied, and the following conclusions were summarized: (1) The displacement and internal force response of pile vibrations in saturated clay foundations have a delayed effect.The stronger the rheological properties of the foundation soil, the more obvious the delay, and the lower the load frequency, the more significant the influence of rheological properties on the delay effect.(2) The stronger the rheological properties of the soil, the smaller the displacement amplitude of the pile foundation, and the higher the load frequency, the greater the decrease in displacement amplitude.The positive bending moment of the pile decreases with the increase in soil rheological properties, while the negative bending moment increases accordingly.

Figure 1 .
Figure 1.Schematic diagram of the fractional Kelvin model.

Figure 1 .
Figure 1.Schematic diagram of the fractional Kelvin model.

Figure 2 .
Figure 2. Horizontal coupled vibration model of pile-soil system.

Figure 2 .
Figure 2. Horizontal coupled vibration model of pile-soil system.

Figure 3 .
Figure 3.Comparison of pile top displacement results (imensionless displacement * p u and di-

Figure 3 .
Figure 3.Comparison of pile top displacement results (imensionless displacement u * p and dimensionless time t * ).

Figure 4 .
Figure 4. Model of finite element analysis.

Figure 4 .
Figure 4. Model of finite element analysis.

Figure 5 .
Figure 5.Comparison of pile top displacement with FEM.

Figure 5 .
Figure 5.Comparison of pile top displacement with FEM.

Figure 6 .Figure 6 .
Figure 6.Influence of shear modulus of pile on the vibration displacement of pile top (ω = 2π, a = 0) Figure 6.Influence of shear modulus of pile on the vibration displacement of pile top (ω = 2π, a = 0).

Figure 6 .Figure 7 .
Figure 6.Influence of shear modulus of pile on the vibration displacement of pile top (ω = 2π, a = 0

Figure 9 .
Figure 9.The variation of pile top displacement amplitude with the order of soil viscosity (Gp = Ep).

1 t
is the time corresponding to the first displacement peak and 2

Figure 10 .
Figure 10.The variation of displacement delay time with the order of soil viscosity (Gp = Ep).

Figure 9 .
Figure 9.The variation of pile top displacement amplitude with the order of soil viscosity (G p = E p ).

Figure 9 . 1 t
Figure 9.The variation of pile top displacement amplitude with the order of soil viscosity (Gp = Ep

Figure 10 .
Figure 10.The variation of displacement delay time with the order of soil viscosity (Gp = Ep).

Figure 10 .
Figure 10.The variation of displacement delay time with the order of soil viscosity (G p = E p ).