A practical numerical method for large strain liquefaction analysis of saturated soils

Accurate prediction of the liquefaction of saturated soils is based on strong coupling between the pore fluid phase and soil skeleton. A practical numerical method for large strain dynamic analysis of saturated soils is presented. The u–p formulation is used for the governing equations that describe the coupled problem in terms of soil skeleton displacement and excess pore pressure. A mixed finite element and finite difference scheme related to large strain analysis of saturated soils based on the updated Lagrangian method is given. The equilibrium equation of fluid-saturated soils is spatially discretized by the finite element method, whereas terms associated with excess pore pressure in the continuity equation are spatially discretized by the finite difference method. An effective cyclic elasto-plastic constitutive model is adopted to simulate the non-linear behavior of saturated soils under dynamic loading. Several numerical examples that include a saturated soil column and caisson-type quay wall are presented to verify the accuracy of the method and its usefulness and applicability to solutions of large strain liquefaction analysis of saturated soils in practical problems.     

Study of pore pressure variation during liquefaction using two constitutive models for sand

Numerical analyses of liquefiable sand are presented in this paper. Liquefaction phenomenon is an undrained response of saturated sandy soils when they are subjected to static or dynamic loads. A fully coupled dynamic computer code is developed to predict the liquefaction potential of a saturated sandy layer. Coupled dynamic field equations of extended Biot's theory with u–P formulation are used to determine the responses of pore fluid and soil skeleton. Generalized Newmark method is employed for integration in time. The soil behavior is modelled by two constitutive models; a critical state two-surface plasticity model, and a densification model. A class ‘B’ analysis of a centrifuge experiment is performed to simulate the dynamic response of level ground sites. The results of the numerical analyses demonstrate the capability of the critical sate two-surface plasticity model in producing pore pressures that are consistent with observations of the behavior of liquefiable sand in the centrifuge test.     

某长江大桥深度超过15米的砂土层液化势的综合评价

不同抗震设计规范的砂土液化判别方法或国内外其他有代表性的液化判别方法所采用的地震动参数和土性指标及其埋藏条件是不同的,因而采用这些方法对同一工程场地进行液化势预测时其评价结果通常有一些差异,甚至会得到相反的结论。为了给重大工程建设提供较为合理、可信的地基液化势预测结果,采用多种液化判别方法进行场地液化势的综合评价是比较客观的,也是必要的。本文结合某长江大桥桥基工程,采用建筑抗震设计规范的砂土液化判别方法、国内外有代表性的液化判别方法、有限元数值分析法等多种方法逐一对该工程场地砂性土层进行液化判别,并结合室内动三轴液化试验结果,对主桥墩不考虑冲刷条件和考虑一般冲刷深度5m条件时的砂性土层进行了液化势的综合评价,并将各土层的液化势分为液化、可能液化和不液化3个等级,得到了较为合理可靠的判别结果。     

Validation of a new endochronic liquefaction model for granular soil by using centrifuge test data

The application of a new liquefaction constitutive model, based on the endochronic theory applied to densification of sandy soil, to a set of centrifuge tests from the University of British Columbia, is presented in this paper. The model employed herein takes into account, in a unified formulation, contractive and dilative behaviours, and considers the soil collapse as well. First of all, the model is calibrated by means of undrained cyclic simple shear stress test data. The constitutive law of the soil is implemented in the bidimensional coupled finite element code CMLIQ (Cyclic Mobility and LIQuefaction), developed by the authors. Three centrifuge tests are analysed, the seismic loading and the geometry being the same for all of them, namely an improved slope with drain devices or denser soil at some locations. Comparisons between the data provided by the numerical model and the experimental measurements are shown, and, as a result, the accuracy of the model is explored and evaluated.     

Cyclic response of saturated silts

Softening and strength loss of sands with increasing excess pore water pressure under repeated loads is well-known. However, extensive damage to the built environment also occurs at the sites underlain by fine grained soils during seismic shaking. The primary objective of this study is to investigate the factors affecting cyclic behavior of saturated low-plastic silt through laboratory testing. For this purpose, an extensive laboratory testing program including conventional monotonic and cyclic triaxial tests was carried out over reconstituted silt samples. The effects of the inherent soil properties and the effects of loading characteristics on the cyclic response of saturated low-plastic reconstituted silt samples were examined separately. Based on the test results, a model was introduced to estimate the effect of initial shear stress on the cyclic response. Besides, liquefaction susceptibility of the samples was examined via current liquefaction susceptibility criteria.     

单、双向动三轴试验条件下饱和砂动力特性对比

在饱和砂土地震液化评估中,多用单向动三轴试验代替双向动三轴试验获得土体的抗液化强度。由于单向激振和双向激振产生的应力路径不同,土的动力特性受应力路径影响,其合理性尚待证明。本文利用双向振动三轴仪分别采用单、双向振动形式进行饱和砂土振动液化试验,对比了相同动剪应力作用下累积孔压和动强度,分析了较大循环振次作用下单、双向振动方式使试样动力响应产生差异的原因,提出了利用单向三轴试验代替双向振动模拟地震循环剪应力获得砂土动力特性的室内试验方法的适用条件。     

A robust empirical model to estimate earthquake-induced excess pore water pressure in saturated and non-saturated soils

In engineering practice, the liquefaction potential of a sandy soil is usually evaluated with a semi-empirical, stress-based approach computing a factor of safety in free field conditions, defined as the ratio between the liquefaction resistance (capacity) and the seismic demand. By so doing, an estimate of liquefaction potential is obtained, but nothing is known on the pore pressure increments (often expressed in the form of normalized pore pressure ratio r_u) generated by the seismic action when the safety factor is higher than 1. Even though r_u can be estimated using complex numerical analyses, it would be extremely useful to have a simplified procedure to estimate them consistent with the stress-based approach adopted to check the safety conditions. This paper proposes such a procedure with reference to both saturated and unsaturated soils, considering the latter as soils for which partial saturation has been artificially generated with some ground improvement technology to increase cyclic strength and thus tackle liquefaction risk. A simple relationship between the liquefaction free field safety factor FS, and r_u(S_r) is introduced, that generalizes a previous expression proposed by Chiaradonna and Flora (Geotech Lett, 2020. https://doi.org/10.1680/jgele.19.00032) for saturated soils. The new procedure has been successfully verified against some experimental data, coming from laboratory constant amplitude cyclic tests and from centrifuge tests with irregular acceleration time histories for soils having different gradings and densities.

     

Shear wave velocity-based liquefaction evaluation in the great Wenchuan earthquake： a preliminary case study

The great Wenchuan earthquake (Ms = 8.0) in 2008 caused severe damage in the western part of the Chengdu Plain. Soil liquefaction was one of the major causes of damage in the plain areas, and proper evaluation of liquefaction potential is important in the definition of the seismic hazard facing a given region and post-earthquake reconstruction. In this paper, a simplified procedure is proposed for liquefaction assessment of sandy deposits using shear wave velocity (Vs), and soil liquefaction from the Banqiao School site was preliminarily investigated after the earthquake. Boreholes were made at the site and shear wave velocities were measured both by SASW and down-hole methods. Based on the in-situ soil information and Vs profiles, the liquefaction potential of this site was evaluated. The results are reasonably consistent with the actual field behavior observed after the earthquake, indicating that the proposed procedure is effective. The possible effects of gravel and fines contents on liquefaction of sandy soils were also briefly discussed.     

Seismic behavior of breakwaters on complex ground by numerical tests: Liquefaction and post liquefaction ground settlements

A large number of breakwaters have been constructed along coasts to protect humans and infrastructures from tsunamis.There is a risk that foundation soils of these structures may liquefy,or partially liquefy during the earthquake preceding a tsunami,which would greatly reduce the structures’capacity to resist the tsunami.It is necessary to consider not only the soil’s liquefaction behavior due to earthquake motions but also its post-liquefaction behavior because this behavior will affect the breakwater’s capacity to resist an incoming tsunami.In this study,numerical tests based on a sophisticated constitutive model and a soil-water coupled finite element method are used to predict the mechanical behavior of breakwaters and the surrounding soils.Two real breakwaters subjected to two different seismic excitations are examined through numerical simulation.The simulation results show that,earthquakes affect not only the immediate behavior of breakwaters and the surrounding soils but also their long-term settlements due to post-earthquake consolidation.A soil profile with thick clayey layers beneath liquefied soil is more vulnerable to tsunami than a soil profile with only sandy layers.Therefore,quantitatively evaluating the seismic behavior of breakwaters and surrounding soils is important for the design of breakwater structures to resist tsunamis.     

饱和黄土场地原位试验及液化势评价

目前,主要依靠室内动力试验对黄土液化势进行评价。由于黄土特殊的结构性,室内试验对其饱和的过程较为复杂,且与实际场地饱和黄土差异明显,导致室内黄土液化试验结果并不能代表现场饱和黄土的抗液化强度。本文选取兰州市西固区寺儿沟村某饱和黄土场地进行钻孔测试,现场实施了标准贯入试验、静力触探试验以及剪切波速测试。应用Robertson的土类指数分类图对该场地不同含水率黄土的土类进行了界定,确定了饱和黄土属于类砂土,有液化势。应用NCEER推荐方法,计算了3组原位试验数据的饱和黄土循环抗力比（CRR）,通过与1976年唐山地震和1999年集集地震液化土CRR对比,得出了饱和黄土抗液化强度很低的结论。     

饱和原状黄土液化基本特征的振动台试验研究

天然黄土因其较强的结构性,制备大尺寸原状试体非常困难,国内外尚无针对饱和原状黄土实施振动台模拟试验的数据与资料。通过解决大尺寸原状黄土试体现场取备难题,利用振动台模拟试验研究饱和原状黄土液化现象及其基本特征。试验结果表明:饱和原状黄土的液化现象,在超孔隙水压力增长、持续与消散的趋势性上与饱和砂土具有可比性,二者的最大差别在于细节特征方面的不同;饱和度是决定地震作用下天然黄土液化特性的首要条件;试体饱和度约为90.3%的条件下,加载后的最大孔压比约为0.93;饱和度85%、75%和65%可能是天然黄土能否发生液化现象、似液化现象(循环失效)和不考虑循环失效现象的临界值。试验获得的资料与分析结果,对深入理解饱和土体液化物理过程与力学机制意义重大。     

Cyclic resistance and liquefaction behavior of silt and sandy silt soils

The liquefaction behavior and cyclic resistance ratio (CRR) of reconstituted samples of non-plastic silt and sandy silts with 50% and 75% silt content are examined using constant-volume cyclic and monotonic ring shear tests along with bender element shear wave velocity (V_s) measurements. Liquefaction occurred at excess pore water pressure ratios (r_u) between 0.6 and 0.7 associated with cumulative cyclic shear strains (γ) of 4% to 7%, after which cyclic liquefaction ensued with very large shear strains and excess pore water pressure ratio (r_u>0.8). The cyclic ring shear tests demonstrate that cyclic resistance ratio of silt and sandy silts decreases with increasing void ratio, or with decreasing silt content at a certain void ratio. The results also show good agreement with those from cyclic direct simple shear tests on silts and sandy silts. A unique correlation is developed for estimating CRR of silts and sandy silts (with more than 50% silt content) from stress-normalized shear wave velocity measurements (V_s1) with negligible effect of silt content. The results indicate that the existing CRR–V_s1 correlations would underestimate the liquefaction resistance of silts and sandy silt soils.     

基于重塑饱和砂土模型的现场液化试验方法

基于现场开展土体液化问题研究势必成为今后土动力学中的一个重要发展方向。目前人工激振下的现场液化试验方法还不够成熟,尚需进一步探索和发展。本文从试验设备组成、场地地震动激励、试坑布置、饱和砂土模型制备、数据测量与采集等5个方面论述该方法中的主要技术问题。研究表明:动力加载系统激励产生的地震动在0~7m/s2;系统工作频率13~15Hz,饱和砂土模型与基础边缘的距离在0.5~2.5m范围内,更适合进行液化试验;应用水沉法现场制备饱和砂土模型,要重点注意试坑防水和尺寸定位的问题;数据测量与采集中要充分考虑对现场液化问题认识不够这一因素的影响,需对数据测量与采集提出附加要求;试验实例初步表明,该方法可行,适合开展液化问题研究。     

遮帘式板桩码头地基地震液化破坏机理

遮帘式板桩码头作为一种新型的板桩结构型式,其抗震性能研究是设计建造过程中的重要环节。在FEM-FDM水土耦合计算的平台上引入循环弹塑性本构模型,借助FORTRAN编程软件形成饱和砂土动力液化分析的数值方法,可有效模拟饱和砂土在地震动力作用下的非线性及大变形特性,同时也可模拟砂土液化流动对遮帘桩和前墙的动土压力。研究表明:地震作用下可液化土层超孔隙水压力比增长并发生较大的水平流动变形,对前墙的水平破坏大于竖向破坏;前墙剪力最大值位于海床与前墙交界处;遮帘桩剪力最大值位移与前墙底平行的位置;后拉杆拉力逐渐变大,前拉杆拉力逐渐变小。通过对板桩码头地震液化灾害的分析,可为抗震和抗液化设计提供参考依据。     

Probabilistic evaluation of liquefaction potential under earthquake loading

This paper presents a procedure to perform the risk analysis for ground failure by liquefaction. The first part of this study describes the differential equation of a smooth hysteretic model to characterize the behavior of the soil under random loading. The parameters of the proposed model to represent the experimental relationship are discussed. The second part of this study is to develop a method to calculate the probability that a specified volume of soil will liquefy at a given depth in the deposit. The liquefaction is defined as the result of cumulative damage caused by seismic loading. The fatigue life of soil can be determined on the basis of the N---S relationship and Miner's cumulative damage law. The rain-flow method is used to count the number of cycles of stress response of the soil deposit. Finally, the probability of liquefaction is obtained by integration over all the possible ground motion and the fragility curves of liquefaction potential. The sensitivity of the reliability against liquefaction to soil system parameters is also examined.     

关于砂土液化判别的若干意见

简要回顾了我国抗震规范中液化判别方法的发展历史,并对我国若干规范中的液化初步判别的条文进行了比较,提出了修改意见;对我国有代表性的几种液化判别方法及基于Seed简化方法修正的美国NCEER法进行了评述和综合对比,指出了我国现行若干规范中有关液化判别规定方面的错误;此外,还将GB50011－2001《建筑抗震设计规范》液化判别公式转变为类似于NCEER法的表达式,拓宽了该规范方法的应用范围。最后,结合某长江大桥工程的地基液化判别工作,探讨了深层土液化判别的可能性,结果表明,本文建议的方法是可行的。     

Cyclic and post-cyclic shear behavior of low-plasticity silt with varying clay content

This paper investigates the cyclic and post-cyclic shear behavior of low-plasticity silt and the impact of additional clay content. Bentonite clay was added to the low-plasticity Mississippi River Valley (MRV) silt (PI=6) to increase the clay content of the soil. A series of triaxial tests were conducted in the laboratory to examine the shear and pore pressure behavior during and after cyclic loading. As the bentonite content in the reconstituted specimens increased, the excess pore pressure developed at a slower rate and the total excess pore pressure decreased at the end of cyclic loading. In contrast to the MRV silt, the specimens modified with bentonite experienced cyclic softening rather than initial flow liquefaction. The cyclic shear strength increased with an increase in bentonite content. The post-cyclic reconsolidation behavior was a similar to a virgin compression process, and not recompression. Adding bentonite to the MRV silt results in changes in permeability, compressibility, undrained shear strength, and initial stiffness. Additionally, the cyclic loading had a marked effect on the shear behavior of low-plasticity soil with a PI<6, but not noticeable with a PI>6. This study suggests that the behavior of the Mississippi River Valley silt changes from contractive sand-like material to clay-like behavior at a PI≈6 due to the addition of clay.     

Evaluation of temperature and freeze–thaw effects on excess pore pressure generation of fine-grained soils

The November 3, 2002 Denali-Alaska earthquake (M_w=7.9) caused significant liquefaction associated damage to various infrastructure built on fine-grained soils. The seismic response, liquefaction potential, and excess pore pressure generation of soils in cold regions, especially those of fine-grained nature, have not been studied thoroughly and therefore are not well-understood. This paper presents results from an extensive laboratory study on the characteristics of excess pore pressure generation and liquefaction potential of fine-grained soils. Laboratory-constituted soils specimens were tested in four categories: (1) tests on specimens subjected to no thermal conditioning or freeze–thaw cycles; (2) tests on specimens conditioned at 24, 5, 1, 0.5, and −0.2 °C; (3) tests on specimens subjected to 1–4 freeze–thaw cycles; and (4) tests on specimens conditioned at near-freezing temperatures of 0.5 and −0.2 °C through different freeze–thaw paths. Strain-controlled, undrained, cyclic triaxial tests were performed at shear strain levels of 0.005–0.8%. Specimens conditioned at different temperatures were found to generate significantly different pore pressures with cyclic loading. The excess pore pressure generation at near or slightly below freezing was found to change dramatically. A transitional change in the dynamic soil behavior, attributed to unfrozen- or frozen-dominant pore water, was discovered. The threshold shear strain was also found to be influenced by the temperature. Subjecting the soil specimens to 1, 2 and 4 freeze–thaw cycles caused a reduction in excess pore pressure generation and slight change in the threshold shear strain. The temperature conditioning path to reach the target temperature was found to be important on the development of excess pore pressure at near and slightly below-freezing temperatures.