Application of the adaptive neuro-fuzzy inference system for prediction of soil liquefaction

The determination of liquefaction potential of soils induced by earthquake is a major concern and an essential criterion in the design process of the civil engineering structures. A purely empirical interpretation of the filed case histories relating to liquefaction potential is often not well constrained due to the complication associated with this problem. In this study, an integrated fuzzy neural network model, called Adaptive Neuro-Fuzzy Inference System (ANFIS), is developed for the assessment of liquefaction potential. The model is trained with large databases of liquefaction case histories. Nine parameters such as earthquake magnitude, the water table, the total vertical stress, the effective vertical stress, the depth, the peak acceleration at the ground surface, the cyclic stress ratio, the mean grain size, and the measured cone penetration test tip resistance were used as input parameters. The results revealed that the ANFIS model is a fairly promising approach for the prediction of the soil liquefaction potential and capable of representing the complex relationship between seismic properties of soils and their liquefaction potential.     

基于MPL神经网络的地震作用下砂土液化评估及预测

砂土在地震的作用下会产生剧烈的液化现象,液化引发的地基失稳会对道路、建筑物、堤坝等各类设施造成严重危害。因此,地震作用下的砂土液化判别预测一直是地质灾害领域研究的热点问题。本文使用过去几十年发生在世界各地的166组地震作用下砂土液化实例数据,通过大量数据训练和参数分析建立了基于机器学习的地震作用下砂土液化判别模型。结果表明,当网络结构为6(输入层)-15(隐藏层)-1(输出层)、训练函数为Levenberg-Marquardt时,对地震液化预测效果较好,最大准确率可达96%。参数分析结果表明不同参数对网络预测准确率影响程度不一:锥端阻力、地表归一化峰值水平加速度影响相对较大;地震震级、总垂向应力、有效垂向应力影响中等;贯入深度对其影响较小。因此在不同网络预测精度要求的条件下,可考虑适当简化输入参数。     

纤维加筋砂土抗液化试验与数值模拟

砂土液化是导致重大地震灾害的主要原因之一。本研究探讨了天然纤维加筋砂土在循环荷载作用下的抗液化性能。在不排水条件下,对具有不同纤维含量的加筋砂土试样进行了一系列循环三轴试验,研究了饱和砂土的液化特性以及循环剪应变幅值、纤维含量对饱和砂土抗液化性能的影响。此外,通过模拟已完成的循环三轴试验,建立了二维有限元数值模型,并对具有不同纤维含量的加筋砂土进行了参数标定。研究结果表明:(1)增加循环剪应变幅值将促进超孔压累积,使得滞回曲线斜率和平均有效应力减小速度加快;(2)纤维的存在能够减缓超孔压的累积,随着纤维含量增加,加筋砂土抗液化能力得到明显提高;(3)标定后的本构模型参数能可靠地用于模拟纤维加筋砂土的液化响应。研究结果为饱和砂土抗液化问题与纤维加筋砂土的数值模拟提供了有价值的参考。     

A cyclic elastoplastic constitutive model and effect of non-associativity on the response of liquefiable sandy soils

A large number of constitutive models for geomaterials, such as soils and rocks, have been proposed over the last three decades. Those models have been implemented into computer codes and have been successfully used to solve practical engineering problems particularly under monotonic loading conditions. Compared with the models for monotonic loadings, more improvements for cyclic models are necessary in order to obtain more accurate predictions for the dynamic behavior of geomaterials, e.g., the behavior during earthquakes. A cyclic elastoplastic model has been developed in this study for sandy soils; it is based on the kinematical hardening rule with a yield function that includes the changes in the stress ratio and the mean effective stress considering the degradation of the yield surface. From a simulation with the present model, it has been found that strong non-associativity leads to a large decrease in the mean effective stress during cyclic deformations under undrained conditions, while the model with the associated flow rule does not. This result is quite important because the mean effective stress becomes almost zero at the state of full liquefaction. Compared with the experimental results, the model can accurately reproduce the cyclic behavior of soil.     

小坡度海底土层地震液化诱发滑移分析方法

地震可使海底砂质、粉质土层液化并导致上部土层的滑移。基于有效应力有限元动力分析方法和Newmark刚性滑块理论,提出了一种计算海底小坡度（≤5o）土层地震液化引起侧向滑移的简化方法。该方法将波浪荷载简化为海底恒定的上覆压力和初始孔压,忽略了海水粘性对海底土层地震反应的影响,利用改进的Seed孔压模型进行动力分析和液化判别,用Newmark滑块理论计算了土层侧向滑移。通过算例和对比分析,研究了海水深度和土层坡度对侧向滑移的影响,表明该方法的有效性,可为近海工程场地地震地质灾害评价提供参考数据。     

Working mechanism of cutoff walls in reducing uplift of large underground structures induced by soil liquefaction

The uplift of large underground structures in saturated liquefiable soils under strong earthquake loadings may induce severe damages to the structures. Various mitigation procedures have been proposed to alleviate such damage, among which installation of cutoff walls next to underground structures was found to be effective. However, the working mechanism of cutoff walls in alleviating uplift of underground structures and the corresponding design parameters are still not clear. The liquefaction induced uplift behaviour of a subway tunnel in saturated sandy deposit over a layer of non-liquefiable soils and the working mechanism of cutoff walls for uplift mitigation purpose were investigated using the fully coupled dynamic finite element code DIANA Swandyne-II. A generalized plasticity model capable of simulating both cyclic liquefaction and pressure dependency of soils was used to model the sandy deposit. It is found that the small effective unit weight of underground structures, the development of excess pore pressure and the flow of liquefied soils were the sufficient and necessary conditions for underground structures to uplift during earthquakes. Cutoff walls could restrain the flow or deformation of liquefied soils and inhibit the uplift of underground structures but they could not necessarily prevent the liquefaction of the enclosed soils. After earthquake loadings, underground structures might settle due to the consolidation of soils and cutoff walls could also reduce the magnitude of settlement. The design parameters of cutoff walls, including the acting lateral pressure, the position, the stiffness and the permeability of cutoff walls, were also analyzed, the findings of which, together with the unveiled working mechanism, would be relevant for the design of cutoff walls for uplift mitigation purpose.     

饱和砂土液化判别与放大效应数值模拟研究

为研究地震作用下饱和砂土液化判别及地震放大效应的影响因素,采用边界面塑性模型框架内开发的砂土本构模型,基于开源有限元平台OpenSees建立了一维剪切梁土柱模型。以循环应力比CSR和循环抗力比CRR为控制指标,对比了不同液化判别方法的差异,分析了地震荷载类型和砂土相对密度对液化判别和放大效应的影响。研究表明:与数值模拟结果相比,Seed简化法计算的CSR更大,判断饱和砂土场地发生液化的可能性更高;冲击型地震波较振动型地震波更容易使饱和砂土场地发生液化,砂土相对密度越小场地越容易发生液化;放大系数随埋深的减小而增大,振动型地震波引起的放大效应整体大于冲击型,埋深较大时放大系数随砂土相对密度的增大而减小。     

含细粒珊瑚土抗液化特性试验研究

珊瑚土作为新兴的热带地区岛礁与港口工程的首选土工材料,宽级配是其结构组成的主要特征。由于相关工程面临较高的地震风险,珊瑚土的抗液化能力正逐渐引起重视。为探究含细粒珊瑚土的抗液化能力,以东太平洋某热带港口工程的实际珊瑚土场地为背景,通过大粒径循环三轴液化试验测试了设计相对密实度为0.4～0.8的3组代表性级配试样及剔除细粒的两组试样的饱和不排水动强度。试验结果表明：幂函数可以模拟含细粒珊瑚土的循环应力比与液化所需振次关系;细粒的存在与相对密实度的提高不会显著提高珊瑚土抗液化能力;珊瑚土液化过程的超静孔压发展模式与砂土相近,两参数或三参数的反正弦模型可以较好地模拟含细粒珊瑚土的液化超孔压发展过程。研究表明,含细粒珊瑚土仍然属于可液化土类。以背景工程为例,同类型工程在设计施工及使用阶段都需要考虑对地震液化灾害的设防,该研究为珊瑚土液化防治工作提供了技术支持。     

A rational method for development of limit state for liquefaction evaluation based on shear wave velocity measurements

This papers presents a new approach for developing a limit state for liquefaction evaluation based on field performance data. As an example to illustrate the new approach, a database that consists of, among many other features, in situ shear wave velocity measurements and field observations of liquefaction/non‐liquefaction in historic earthquakes is analysed. This database is first used to train a neural network to classify liquefaction/non‐liquefaction based on soil resistance parameters and load parameters. The successfully trained and tested neural network is then used to establish a limit state, a multiple dimension boundary that separates ‘zone’ of liquefaction from ‘zone’ of non‐liquefaction. The limit state yields cyclic resistance ratio for a given set of soil resistance parameters. Examination of all cases in the database show that the developed limit state has a high degree of accuracy in predicting the occurrence of liquefaction/non‐liquefaction. The developed neural network model can accurately predict the cyclic resistance ratio of soils. Copyright © 2000 John Wiley & Sons, Ltd.     

On integration of a cyclic soil plasticity model

Performance of three classes of explicit and implicit time‐stepping integrators is assessed for a cyclic plasticity constitutive model for sands. The model is representative of an important class of cyclic plasticity models for soils and includes both isotropic and nonlinear kinematic hardening. The implicit algorithm is based on the closest point projection method and the explicit algorithm follows a cutting‐plane integration procedure. A sub‐stepping technique was also implemented. The performance of these algorithms is assessed through a series of numerical simulations ranging from simulations of laboratory tests (such as triaxial and bi‐axial compression, direct shear, and cyclic triaxial tests) to the analysis of a typical boundary value problem of geotechnical earthquake engineering. These simulations show that the closest point projection algorithm remains stable and accurate for relatively large strain increments and for cases where the mean effective stress in a soil element reaches very small values leading to a liquefaction state. It is also shown that while the cutting plane (CP) and sub‐stepping (SS) algorithms provide high efficiency and good accuracy for small to medium size strain increments, their accuracy and efficiency deteriorate faster than the closest point projection method for large strain increments. The CP and SS algorithms also face convergence difficulties in the liquefaction analysis when the soil approaches very small mean effective stresses. Copyright © 2001 John Wiley & Sons, Ltd.     

Numerical simulation of mitigation for liquefaction-induced soil deformations in a sandy ground improved by cement grouting

This paper presents a numerical study of mitigation for liquefaction during earthquake loading. Analyses are carried out using an effective stress based, fully coupled, hybrid, finite element-finite differences approach. The sandy soil behavior is described by means of a cyclic elastoplastic constitutive model, which was developed within the framework of a nonlinear kinematic hardening rule. In theory, the philosophies of mitigation for liquefaction can be summarized as two main concepts, i.e. prevention of excess pore water pressure generation and reduction of liquefaction-induced deformations. This paper is primarily concerned with the latter approach to liquefaction mitigation. Firstly, the numerical method and the analytical procedure are briefly outlined. Subsequently, a case-history study, which includes a liquefaction mitigation technique of cement grouting for ground improvement of a sluice gate, is conducted to illustrate the effectiveness of liquefaction countermeasures. Special emphasis is given to the computed results of excess pore water pressures, displacements, and accelerations during the seismic excitation. Generally, the distinctive patterns of seismic response are accurately reproduced by the numerical simulation. The proposed numerical method is thus considered to capture the fundamental aspects of the problems investigated, and yields results for design purposes. From the results in the case, excess pore water pressures eventually reach fully liquefied state under the input earthquake loading and this cannot be prevented. However, liquefaction-induced lateral spreading of the foundation soils can be effectively reduced by the liquefaction mitigation techniques. An erratum to this article can be found at     

Simulation of liquefaction of unsaturated soil using critical state soil model

Several researchers have reported that the mean effective stress of unsaturated soils having a relatively high degree of saturation gradually decreases under fully undrained cyclic loading conditions, and such soils can be finally liquefied like saturated soils. This paper describes a series of simulations of fully undrained cyclic loading on unsaturated soils, conducted using an elastoplastic model for unsaturated soils. This model is a critical state soil model formulated using effective stress tensor for unsaturated soils, which incorporates the following concepts: (a) the volumetric movement of the state boundary surface containing the critical state line owing to the variation in the degree of saturation; (b) the soil water characteristic curve considering the effects of specific volume and hydraulic hysteresis; and (c) the subloading surface concept for considering the effect of density. Void air is assumed to be an ideal gas obeying Boyle's law. The proposed model is validated through comparisons with past results. The simulation results show that the proposed model properly describes the fully undrained cyclic behavior of unsaturated soils, such as liquefaction, compression, and an increase in the degree of saturation. Finally, the effects of the degree of saturation, void ratio, and confining pressure on the cyclic strength of unsaturated soils are described by the simulation results. The liquefaction resistance of unsaturated soils increases as the degree of saturation and the void ratio decrease, and as the confining pressure increases. Furthermore, the degree of saturation has a greater effect on the liquefaction resistance than the confining pressure and void ratio. Copyright © 2017 John Wiley & Sons, Ltd.     

基于耗散能量的饱和黄土动孔压模型

通过一系列不排水动三轴试验探究了饱和黄土振动液化过程中孔隙水压力和累积耗散能量的演化模式,并讨论了围压、动应力幅值和固结应力比对其演化过程的影响。结果表明：饱和黄土的孔隙水压力和耗散能量随着循环荷载作用逐渐累积。固结围压抑制孔隙水压力增长而消耗更多能量;更大的动应力幅值使得孔隙水压力增长更快而消耗能量更少;等压固结下,孔隙水压力增长至围压从而触发初始液化,而偏压固结下,通常先达到振动液化应变标准而孔隙水压力并没有增长至围压水平,并且固结应力比越大,液化时孔隙水压力越小,消耗能量也更少。归一化孔隙水压力u/σ₀’与累积耗散能量W/W_f之间关系受围压、循环应力比和固结应力比影响较小,可统一用双曲线模型表示。     

黄河口不同强度粉土液化特性的试验研究

研制真空压缩装置以制备不同强度的粉土土样,并进行动三轴试验,测试不同强度粉土的液化特性。建立了土的强度与动剪应力比的关系,确立孔压增长模型及其参数,探讨土的强度和动应力对液化的影响。结果表明,（1）真空压缩装置制备土样过程中,土体超孔压在24 h内基本消散完毕,28 d贯入阻力达300～400 N,不排水抗剪强度达8 kPa,达到黄河口软弱土的强度;（2）土的强度与动剪应力比基本呈线性关系,且土的强度越高,孔压增长曲线越呈现上凸趋势,破坏时的孔压比也越大;（3）指数模型能够较好模拟黄河口粉土孔压增长情况,其中参数a和b分别位于0.77～5.63和0.17～4.65之间;对于孔压比上限,参数a, b分别为0.92和4.65;孔压比下限,参数a, b分别为1.25和0.89。     

土工袋垫层抗液化性能振动台试验研究

土工袋能够有效地约束袋内土体,具有良好的滤水保土作用,可以用于地基抗液化,但土工袋抗液化性能尚未有系统及深入的研究。开展了一系列小振动台试验,验证土工袋垫层的抗液化效果,研究了振动加速度、土工袋层数和排列方式对抗液化性能的影响。结果表明：土工袋具有良好的抗液化效果,袋内土体的超静孔压比小于同深度处周围的土体;土工袋垫层的排水性能是在土工袋本身具有良好透水性的基础上,孔隙水会沿着土工袋与土体的界面以及袋间的空隙排出;相较于不透水刚性垫层,振动过程中土工袋垫层表面基本保持水平,发挥出较好的变形协调性;土工袋层数增加及交错式排列对抗液化效果有利。     

防水措施对堤坝地震液化影响的数值分析

通过基于有效应力的完全耦合动力分析方法,研究了混凝土面板坝在地震作用下的动力响应,详细地讨论了堤坝迎水面斜坡上铺设的防水毯对堤坝地震液化的影响。计算中土体采用Cyclic mobility本构模型,该模型通过在修正剑桥模型上增加应力诱导各向异性、超固结和结构性的概念及其相关发展准则,可以很好地描述可液化土体的动力特性。计算结果表明防水毯的存在可以有效降低坝体的地下水位浸润线,增加坝基的初始有效应力,从而降低坝基地震液化发生的可能性,同时坝体的变形明显减小。     

饱和击实黄土的动力特性研究

通过进行不同固结条件下饱和击实黄土的动三轴试验,研究了饱和击实黄土的动模量、阻尼比、动强度、动孔压及抗液化特性。研究结果表明：饱和击实黄土的动应力-应变关系符合双曲线模型,模型中参数起始动剪切模量和最大动应力与轴向固结应力间均有良好的幂函数关系,且可以对不同固结应力状态归一,固结围压和固结比对阻尼比的影响较小。动剪应力比随固结围压的增大而减小,随固结比的增大而增大。固结围压、固结比以及动应力皆对动孔压比（ ）与振次比关系有显著的影响,而动孔压与破坏时动孔压之比与振次比关系只受固结围压变化的影响,基本上不受固结比和动应力变化的影响,可以用幂函数关系来模拟;在均压固结条件下,当破坏振次小于等于30时,饱和击实黄土不会产生液化,而当破坏振次较大（动应力较小）时可以产生液化;在偏压固结条件下不会产生液化。     

Laboratory Investigation of the Resistance of Tailings and Natural Sediments to Cyclic Loading

A number of cyclic triaxial tests were carried out on mine tailings and natural sediment samples under undrained conditions to investigate their resistance to cyclic loading. The tests were performed on more than 100 samples with a cyclic shear stress ratio ranging from 0.10 to 0.40 under varying void ratio and the same confining pressure. It was observed that the axial strain and excess pore water pressure increased with the number of loading cycles while the effective stress decreased with increasing number of loading cycles. The liquefaction resistance of the tailings was also observed to be higher than that of natural soils with similar particle size distribution, void ratio and plasticity index. It was observed that the influence of specific gravity on the cyclic strength of mine tailings is significant. The results showed that the cyclic resistance of the tailings was not strongly influenced by plasticity index for low plasticity tailings. A boundary relationship between void ratio and normalized cyclic resistance ratio was established based on the results.     

可液化场地碎石桩复合地基地震动力响应分析

采用砂土液化大变形弹塑性本构模型分析可液化砂土,采用模量随应力与应变变化的等效非线性模型增量形式分析碎石桩,应用FLAC3D有限差分软件对地震动力作用下可液化场地碎石桩复合地基进行三维动力响应分析。模拟分析了在地震作用下碎石桩刚度效应和排水效应对加固处理可液化场地的抗液化效果,从初始小变形到液化后大变形的变形发展,超静孔压累积与消散,及桩与土的变形与应力分配变化等。结果表明,所用模型与方法可合理描述可液化场地碎石桩复合地基在地震作用下场地的动力响应特性和抗液化效果;在地震作用下可液化场地中桩周土体与碎石桩体的竖向应力与水平向剪切应力向碎石桩体集中,竖向有效应力比可降至约1/6～1/3;桩周土体与桩体为非协调变形,剪应变比可达7～10;碎石桩抗液化影响范围约为2.5～3倍桩径,对超过3.5倍桩径范围影响较小;碎石桩与砂土渗透系数比大于100时对降低砂土中超静孔隙水压影响明显;碎石桩对场地的加密效应可显著降低超静孔隙水压力,而碎石桩刚度则对超静孔隙水压力变动影响较小,但有助于减低地面加速度响应峰值。     

Test Investigation on Liquefied Deformation Structure in Saturated Lime-Mud Composites Triggered by Strong Earthquakes

Three identical model boxes were made from transparent plexiglass and angle iron. Using the method of sinking water and according to the sedimentary rhythm of saturated calcium carbonate(lime-mud) intercalated with cohesive soil,calcites with particle sizes diameters of ≤ 5 μm,10–15 μm and 23–30 μm as well as cohesive soil were sunk alternatively in water of three boxes to build three test models,each of which has a specific size of calcite. Pore water pressure gauges were buried in lime-mud layers at different depths in each model,and connected with a computer system to collect pore water pressures. By means of soil tests,physical property parameters and plasticity indices(Ip) were obtained for various grain-sized saturated lime-muds. The lime-muds with Ip ranging from 6.3 to 8.5(lower than 10) are similar to liquid saturated silt in the physical nature,indicating that saturated silt can be liquefied once induced by a strong earthquake. One model cart was pushed quickly along the length direction of the model so that its rigid wheels collided violently with the stone stair,thus generating an artificial earthquake with seismic wave magnitude greater than VI degree. When unidirectional cyclic seismic load of horizontal compression-tension-shear was imposed on the soil layers in the model,enough great pore water pressure has been accumulated within pores of lime-mud,resulting in liquefaction of lime-mud layers. Meanwhile,micro-fractures formed in each soil layer provided channels for liquefaction dewatering,resulting in formation of macroscopic liquefaction deformation,such as liquefied lime-mud volcanoes,liquefied diapir structures,vein-like liquefied structures and liquefied curls,etc. Splendid liquefied lime-mud eruption lasted for two to three hours,which is similar to the sand volcano eruption induced by strong earthquake. However,under the same artificial seismic conditions,development of macroscopic liquefied structures in three experimental models varied in shape,depth and quantity,indicating that excess pore water pressure ratios at initial liquefaction stage and complete liquefaction varied with depth. With size increasing of calcite particle in lime-mud,liquefied depth and deformation extent increase accordingly. The simulation test verifies for the first time that strong earthquakes may cause violent liquefaction of saturated lime-mud composed of micron-size calcite particles,uncovering the puzzled issue whether seafloor lime-mud can be liquefied under strong earthquake. This study not only provides the latest simulation data for explaining the earthquake-induced liquefied deformations of saturated lime-mud and seismic sedimentary events,but also is of great significance for analysis of foundation stability in marine engineering built on the soft calcium carbonate layers in neritic environment.