Pile foundation in alternate layered liquefiable and non-liquefiable soil deposits subjected to earthquake loading

Pile foundations are still the preferred foundation system for high-rise structures in earthquake-prone regions.Pile foundations have experienced failures in past earthquakes due to liquefaction.Research on pile foundations in liquefiable soils has primarily focused on the pile foundation behavior in two or three-layered soil profiles.However,in natural occurrence,it may occur in alternative layers of liquefiable and non-liquefiable soil.However,the experimental and/or numerical studies on the layered effect on pile foundations have not been widely addressed in the literature.Most of the design codes across the world do not explicitly mention the effect of sandwiched non-liquefiable soil layers on the pile response.In the present study,the behavior of an end-bearing pile in layered liquefiable and non-liquefiable soil deposit is studied numerically.This study found that the kinematic bending moment is higher and governs the design when the effect of the sandwiched non-liquefied layer is considered in the analysis as opposed to when its effect is ignored.Therefore,ignoring the effect of the sandwiched non-liquefied layer in a liquefiable soil deposit might be a nonconservative design approach.

Pile Running in Layered Soils

This paper presents a case study on incidents of offshore pile running in layered soils.The study provides a detailed description of the seabed soil data,pile driving records,and field surveillance video observations.Three-dimensional large deformation finite element(LDFE)analyses were conducted to retrospectively analyze the incidents,considering the remoulding of seabed soil and degradation of the pile-soil interface in the LDFE modeling.By comparing the field observations with the LDFE analysis,the mechanism of pile running was discussed,with a focus on investigating the pile penetration resistance in each layer.The study revealed that pile running in layered soils primarily resulted from a significant reduction in pile base resistance when transitioning from a strong layer to an adjacent weak layer.To further investigate the pile running mechanism in layered soils,a parametric study on the strength variation of adjacent soil layers and its influence on pile base resistance was conducted.Lastly,a simplified prediction model of pile base resistance,suitable for assessing the risk of pile running in layered soils,was proposed.

Mechanisms to explain soil liquefaction triggering,development,and persistence during an earthquake

Mechanisms have been proposed to explain the triggering,development,and persistence of soil liquefaction.The mechanism explaining the horizontal failure plane(triggering)and its depth below the phreatic surface is governed by the flux properties and effective stress at that plane.At the failure plane,the pore water pressure was higher than the effective stress,and the volume change was the highest.The pore water pressure is a function of the soil profile features(particularly the phreatic zone width)and bedrock motion(horizontal acceleration).The volume change at the failure plane is a function of the intrinsic permeability of the soil and bedrock displacement.The failure plane was predicted to occur during the oscillation with the highest amplitude,disregarding further bedrock motion,which was consistent with low seismic energy densities.Two mechanisms were proposed to explain the persistence of soil liquefaction.The first is the existence of low-permeability layers in the depth range in which the failure planes are predicted to occur.The other allows for the persistence and development of soil liquefaction;it is consistent with homogeneous soils and requires water inflow from bedrock water springs.The latter explains many of the features of soil liquefaction observed during earthquakes,namely,surficial effects,“instant”liquefaction,and the occurrence of short-and long-term changes in the level of the phreatic surfaces.This model(hypothesis),the relationship between the flux characteristics and loss of soil shear strength,provides self-consistent constraints on the depth below the phreatic surfaces where the failure planes are observed(expected to occur).It requires further experimental and observational evidence.Similar reasoning can be used to explain other saturated soil phenomena.

Longitudinal dynamic response of pile in layered soil based on virtual soil pile model

Taking the effect of finite soil layers below pile end into account,the longitudinal dynamic response of pile undergoing dynamic loading in layered soil was theoretically investigated.Firstly,finite soil layers below pile end are modeled as virtual soil pile whose cross-section area is the same as that of the pile and the soil layers surrounding the pile are described by the plane strain model.Then,by virtue of Laplace transform and impedance function transfer method,the analytical solution of longitudinal dynamic response at the pile head in frequency domain is yielded.Also,the semi-analytical solution in time domain undergoing half-cycle sine pulse at the pile head is obtained by means of inverse Laplace transform.Based on these solutions,a parametric study is conducted to analyze emphatically the effects of parameters of soil below pile end on velocity admittance and reflected wave signals at the pile head.Additionally,a comparison with other models with different supporting conditions from soil below pile end is performed to verify the model presented.

Longitudinal vibration of pile in layered soil based on Rayleigh-Love rod theory and fictitious soil-pile model

The dynamic response of pile in layered soil is theoretically investigated when considering the transverse inertia effect.Firstly, the fictitious soil-pile model is employed to simulate the dynamic interaction between the pile and the soil layers beneath pile toe. The dynamic interactions of adjacent soil layers along the vertical direction are simplified as distributed Voigt models.Meanwhile, the pile and fictitious soil-pile are assumed to be viscoelastic Rayleigh-Love rods, and both the radial and vertical displacement continuity conditions at the soil-pile interface are taken into consideration. On this basis, the analytical solution for dynamic response at the pile head is derived in the frequency domain and the corresponding quasi-analytical solution in the time domain is then obtained by means of the convolution theorem. Following this, the accuracy and parameter value of the hypothetical boundaries for soil-layer interfaces are discussed. Comparisons with published solution and measured data are carried out to verify the rationality of the present solution. Parametric analyses are further conducted by using the present solution to investigate the relationships between the transverse inertia effects and soil-pile parameters.

Deep Penetration of Spudcan Foundation into Double Layered Soils

The spudcan foundation has been widely used in offshore engineering for jack-up rigs. However, “punch through' failure often occurs where a stronger soil layer overlays a softer soil layer. In this study, spudcan penetration into double layered soils is investigated numerically. The soil profile is set up as a stronger soil layer overlaying a softer soil layer, with the soil strength ratio (bottom soil strength / top soil strength) varied from 0.1 to 1.0 (1 means uniform soil). The bearing behaviour is discussed and the bearing capacity factors are given for various cases involving different layer thicknesses and different strength ratios of the two clay layers. The development of the plastic zones and the effect of soil self-weight on the bearing capacity are also discussed. From this study, it is found that, when a spudcan is distant from the soil layer boundary, the spudcan can be analysed with single soil layer data. However, when a spudcan becomes closer to the soil boundary layer, the influence of the lower soft soil layer is significant, and the bearing capacity of the spudcan decreases. The critical distance is an indication of the occurrence of “punch through' failure. The critical distance between the spudcan and the layer boundary is larger for a rough spudcan than the one for a smooth one, and the critical distance decreases with increasing soil strength ratio. The depth of cavity formed during initial spudcan penetration depends on the top layer soil strength, soil strength ratio and unit soil self-weight, and the cavity affects the spudcan bearing behaviour as well.

A theoretical analysis of vertical dynamic response of large-diameter pipe piles in layered soil

Considering the viscous damping of the soil and soil-pile vertical coupled vibration,a computational model of large-diameter pipe pile in layered soil was established.The analytical solution in frequency domain was derived by Laplace transformation method.The responses in time domain were obtained by inverse Fourier transformation.The results of the analytical solution proposed agree well with the solutions in homogenous soil.The effects of the shear modulus and damping coefficients of the soil at both outer and inner sides of the pipe pile were researched.The results indicate that the shear modulus of the outer soil has more influence on velocity admittance than the inner soil.The smaller the shear modulus,the larger the amplitude of velocity admittance.The velocity admittance weakened by the damping of the outer soil is more obvious than that weakened by the damping of the inner soil.The displacements of the piles with the same damping coefficients of the outer soil have less difference.Moreover,the effects of the distribution of soil layers are analyzed.The results indicate that the effect of the upper soil layer on dynamic response of the pipe pile is more obvious than that of the bottom soil layer.A larger damping coefficient of the upper layer results in a smaller velocity admittance.The dynamic response of the pipe pile in layered soil is close to that of the pipe pile in homogenous soil when the properties of the upper soil layer are the same.

Mode-based equivalent multi-degree-of-freedom system for one-dimensional viscoelastic response analysis of layered soil deposit

Discrete models such as the lumped parameter model and the finite element model are widely used in the solution of soil amplification of earthquakes. However, neither of the models will accurately estimate the natural frequencies of soil deposit, nor simulate a damping of frequency independence. This research develops a new discrete model for onedimensional viscoelastic response analysis of layered soil deposit based on the mode equivalence method. The new discrete model is a one-dimensional equivalent multi-degree-of-freedom（MDOF） system characterized by a series of concentrated masses, springs and dashpots with a special configuration. The dynamic response of the equivalent MDOF system is analytically derived and the physical parameters are formulated in terms of modal properties. The equivalent MDOF system is verified through a comparison of amplification functions with the available theoretical solutions. The appropriate number of degrees of freedom（DOFs） in the equivalent MDOF system is estimated. A comparative study of the equivalent MDOF system with the existing discrete models is performed. It is shown that the proposed equivalent MDOF system can exactly present the natural frequencies and the hysteretic damping of soil deposits and provide more accurate results with fewer DOFs.

Stochastic seismic response of multi-layered soil with random layer heights

This paper deals with the effect of layer height randomness on the seismic response of a layered soil. These parameters are assumed to be lognormal random variables. The analysis is carried out via Monte Carlo simulations coupled with the stiffness matrix method. A parametric study is conducted to derive the stochastic behavior of the peak ground acceleration and its response spectrum,the transfer function and the amplification factors. The input soil characteristics correspond to a site in Mexico City and the input seismic accelerations correspond to the Loma Prieta earthquake. It is found that the layer height heterogeneity causes a widening of the frequency content and a slight increase in the fundamental frequency of the soil profile,indicating that the resonance phenomenon is a concern for a large number of structures. Variation of the layer height randomness acts as a variation of the incident angle,i.e.,a decrease of the amplitude and a shift of the resonant frequencies.

Three-dimensional consolidation deformation analysis of porous layered soft soils considering asymmetric effects

Long-term settlements for underground structures, such as tunnels and pipelines, are generally observed after the completion of construction in soft clay. The soil consolidation characteristic has great influences on the long-term deformation for underground structures. A three-dimensional consolidation analysis method under the asymmetric loads is developed for porous layered soil based on Biot's classical theory. Time-displacement effects can be fully considered in this work and the analytical solutions are obtained by the state space approach in the Cartesian coordinate. The Laplace and double Fourier integral transform are applied to the state variables in order to reduce the partial differential equations into algebraic differential equations and easily obtain the state space solution. Starting from the governing equations of saturated porous soil, the basic relationship of state space variables is established between the ground surface and the arbitrary depth in the integral transform domain. Based on the continuity conditions and boundary conditions of the multi-layered pore soil model, the multi-layered pore half-space solutions are obtained by means of the transfer matrix method and the inverse integral transforms. The accuracy of proposed method is demonstrated with existing classical solutions. The results indicate that the porous homogenous soils as well as the porous non-homogenous layered soils can be considered in this proposed method. When the consolidation time factor is 0.01, the value of immediate consolidation settlement coefficient calculated by the weighted homogenous solution is 27.4% bigger than the one calculated by the non-homogeneity solution. When the consolidation time factor is 0.05, the value of excess pore water pressure for the weighted homogenous solution is 27.2% bigger than the one for the non-homogeneity solution. It is shown that the material non-homogeneity has a great influence on the long-term settlements and the dissipation process of excess pore water pressure.

Axisymmetric Finite Element Method for Analysis of Plate on Layered Soil

To obtain the fundamental solution of soil has become the key problem for the semi-analytical and semi-numerical (SASN) method in analyzing plate on layered soil. By applying axisymmetric finite element method (FEM),an expression relating the surface settlement and the reaction of the layered soil can be obtained. Such a reaction can be treated as load acting on the applied external load. Having the plate modelled by four-node elements,the governing equation of the plate can be formed and solved. In this case, the fundamental solution can be introduced into the global soil stiffness matrix and five-node or nine-node element soil stiffness matrix.The existing commercial FEM software can be used to solve the fundamental solution of soil, which can bypass the complicated formula derivation and boasts high computational efficiency as well.

Nonlinear analysis of pile load-settlement behavior in layered soil

A simplified approach is presented to analyze the single pile settlement in multilayered soil. First, a fictitious soil-pile model is employed to consider the effect of layered soil beneath pile toe on pile settlement behavior. Two approximation methods are proposed to simplify the nonlinear load transfer function and simulate the nonlinear compression of fictitious soil-pile, respectively. On this basis, an efficient program is developed. The procedures for determining the main parameters of mathematical model are discussed. Comparisons with two well-documented field experimental pile loading tests are conducted to verify the rationality of the present method. Further studies are also made to evaluate the practicability of the proposed approach when a soft substratum exists, and the results suggest that the proposed method can provide a constructive means for assessing the settlement of a single pile for use in engineering design.

ONE-DIMENSIONAL CONSOLIDATION OF LAYERED SOILS WITH IMPEDED BOUNDARIES UNDER TIME-DEPENDENT LOADINGS

On the basis of Terzaghi's one-dimensional consolidation theory, the variation of effective stress ratio in layered saturated soils with impeded boundaries under time-dependent loading was studied. By the method of Laplace transform, the solution was presented. Influences of different kinds of cyclic loadings and impeded boundaries conditions were discussed. Through numerical inversion of Laplace transform, useful illustrations were given considering several common time-dependent loadings. Pervious or impervious boundary condition is just the special case of the problem considered here. Compared with average index method,the results from the method illustrated are more accurate.

In-Situ Evaluation of Evaporation in Layered Soils （Tukulu, Sepane and Swartland）

Knowledge about the influence of soil layers on evaporation is essential for the optimization of infield rainwater harvesting （IRWH） in the semi-arid areas of the Free State province of South Africa. Among the soils earmarked for 1RWH development include the Tukulu, Sepane and Swartland soil types that have contrasting soil layers. These soils have to capture and store rainwater within the soil profile layers away from the evaporation zone. To determine how the three soils release and deliver soil water at the evaporating site, a 21-day evaporation experiment was conducted on pre-drained monoliths. Instantaneous soil water content （SWC） from in-situ and soil water characteristic curve （SWCC） from laboratory was measured. Separate soil samples of 15 mm thickness were also evaporated under the same conditions to establish the extent of drying and hydraulic gradient at the soil surface. The Darcian evaporative flux and unsaturated hydraulic conductivity （K-coefficient） were also determined. At the surface suctions of magnitude greater than 1,500 kPa were observed from all monoliths. Total contributions to evaporation from the Tukulu, Sepane and Swartland were 43, 51 and 70 mm, respectively. The low contributions were explained by the presence of the prismacutanic C-horizon in the Tukulu and Sepane at respective depths of 600 and 700 mm. This layer was associated with the steepest suction gradient that restrained further upward fluxes by subsequent lowering for the K-coefficient with more than two orders of magnitudes within a narrow range of SWC. However, the presence of the pedocutanic B-horizon at depths of 300 mm undermined this restrictive function through the appreciable capillary activity demonstrated by clays at near evaporating surfaces. The shallowness and deficiency in structure of the Swartland was consistent with the high contribution to evaporation that gave this soil a dry soil water regime. It was therefore concluded that the Tukulu offered soil profile layers that could reasonably satisfy the soil water conservation requirements for IRWH.

Analysis of Interaction Factors between Two Piles in Layered Soils

A rigorous analytical method is presented, which takes into account the pile stiffening effects, using the theory of the transfer matrix-bottom rigidity for calculating the interaction factor between two identical piles in multi-layered soils subjected to vertical loads. Following the technique proposed by Muki and Sternberg, the problem is decomposed into extended soil layers and two fictitious piles characterized respectively by Young's moduli of the layered soils and those of the differences between the piles and the layered soils. The unknown axial forces along fictitious piles are determined by solving a Fredholm integral equation of the second kind, which imposes the compatibility condition that the axial strains of the fictitious piles are equal to those corresponding to the centroidal axes of the extended layered soils. The real pile displacements can be calculated based on the determined fictitious pile forces, and finally, the desired pile interaction factors may be obtained. Selected results from parametrical studies are presented to confirm the validity of the proposed approach and portray the influence of the governing parameters on the pile interaction.

Theoretical Formula for Calculating Negative Skin Friction of Pile in Layered Soil and Its Application

Based on Zeevaert＇s method, a theoretical formula was developed to calculate the negative skin friction of pile in layered soil. For practical purpose,a cut-and-try method was proposed to determine neutral point. Case studies indicate that the total calculated negative skin friction was in agreement with the measured one, which verifies the feasibility and practicability of theoretical formula. Furthermore, the methods for calculating efficiency factor of drag load and settlement were also given.

Soil Moisture Variation under Different Irrigation Quota

[Objective] The aim was to study variation of soil moisture under different irrigation quota.[Method] By using Trime-TDR apparatus,soil moisture with different irrigation quota infiltration was measured;combining the characteristics of soil texture,curve characteristics of soil moisture variation with soil depth under different irrigation quota were analyzed.[Result] Different irrigation quota has resulted in variation of soil moisture in different layer depth.Soil moisture is 9.88%,17%,25% and 24.45% in so...

A Simplified Method for the Stress Analysis of Underground Transfer Structures Crossing Multiple Subway Tunnels

According to the design specifications,the construction of extended piles involves traversing the tunnel’s upper region and extending to the underlying rock layer.To address this challenge,a subterranean transfer structure spanning multiple subway tunnels was proposed.Deliberating on the function of piles in the transfer structure as springs with axial and bending stiffness,and taking into account the force balance and deformation coordination conditions of beams and plates within the transfer structure,we established a simplified mechanical model that incorporates soil stratification by combining it with the Winkler elastic foundation beam model.The resolved established simplifiedmechanicalmodel employed finite difference technology and the Newton-Simpsonmethod,elucidating the mechanical mechanism of the transfer structure.The research findings suggest that the load carried by the upper structural columns can be transferred to the pile foundation beneath the beams through the transfer structure,subsequently reaching the deep soil layer and ensuring minimal impact on adjacent tunnels.The established simplified analysis method can be used for stress analysis of the transfer structure,concurrently considering soil stratification,pile foundation behavior,and plate action.The pile length,pile section size,and beam section size within the transfer structure should account for the characteristics of the upper load,ensuring an even distribution of the beam bending moment.

Preliminary Analysis on Abnormal Granularity Layers of Soil Profile and the Response of Relative Factors in Loess Plateau

[Objective] The paper was to analyze organic carbon content （SOC）, granularity, total nitrogen content （TN）, carbon-nitrogen ratio （C/N）, calcium carbonate content （CaCO3） of 1cm soil profiles in returning forest in Zhifanggou watershed of Ansai County in Loess Plateau, so as to study the changes of physical and chemical properties in abnormal layer of soil reflected with granularity, as well as the physical and chemical responses of soil. [Method] Three quadrats with the size of 10 m×10 m were randomly selected in three sampling plots in Loess Plateau, three profiles in upper, middle and lower slope were excavated, and the samples were collected with interval of 10 cm; the surface layer with the depth of 0-10 cm was divided into two layers of 0-5 and 5-10 cm for sampling, respectively. Eleven samples were collected in each profile with a total of 99 samples. Its organic carbon content, granularity, total nitrogen content, carbon-nitrogen ratio and CaCO3 content were analyzed. [Result] The soil profiles in three sampling sites contained five characteristic layers, including a1, b1, b2, c1 and c2, the content of soil granule with particle size less than 0.02 mm decreased, and those with particle size 0.02 mm increased, the organic carbon content and C/N value （a1, b1, b2, c2） increased, but the increase trend of CaCO3 content was not obvious. [Conclusion] The study shows that the characteristic soil layer is commonly existed in loess region, especially the eroded loess region, which should be paid attention in the research fields of modern soil science and ecology.

Degradation of Soil Properties due to Erosion on Sloping Land in Southern Jiangsu Province,China

Soil erosion accelerates soil degradation. Some natural soils and cultivated soils on sloping land in southern Jiangsu Province, China were chosen to study soil degradation associated with erosion. Soil erosion intensity was investigated using the 137Cs tracer method. Soil particle-size distribution, soil organic matter (OM), total nitrogen (TN) and total phosphorus (TP) were measured, and the effects of erosion on soil physical and chemical properties were analyzed statistically using SYSTAT8.0. Results indicated that erosion intensity of cultivated soils was greater than that of the natural soils, suggesting that cultivation increased soil loss. Erosion also led to an increase of coarser soil particle proportion, especially in natural soils. In addition, silt was the primary soil particle lost due to erosion. However, in cultivated fields, coarser soil particles over time were attributed not only to soil erosion but also to mechanical eluviation as a result of farming activities. Moreover, erosion caused a decrease in soil OM, TN and TP as well as thinning of the soil layer.