Numerical analysis on seepage failures of dike due to water level-up and rainfall using a water–soil-coupled smoothed particle hydrodynamics model

For seepage failures of dike due to water level-up and rainfall, surface infiltration and strength change induced by suction reduction are important factors; thus, numerical analysis should consider the coupling of water and soil, as well as the effect of saturation to obtain more precise failure mechanism. Based on the advanced smoothed particle hydrodynamics (SPH) method, this work proposed a two-phase-coupled SPH model in coordination with a novel constitutive model for unsaturated soils. Then, a triaxial compression test is simulated to check the applicability of the SPH method on the soil phase. After that, the failure test of a dike due to water level-up is discretized and simulated, from which the seepage process, the distribution of maximum shear strain, the slip surface, and pore water pressure are obtained. The two-phase-coupled SPH model is also applied to a slope failure test of heavy rainfall, and the results are compared to the model test. Finally, a dike failure test due to rainfall is analyzed using the proposed SPH model to reproduce the surface infiltration and suction reduction. The proposed SPH model provides several insights of seepage failures and can be a helpful tool for the analysis of dike failures induced by water level-up and rainfall.     

地震动加速度幅值对土坡失稳后滑坡体大小及致灾范围的影响

地震滑坡的致灾范围是判断滑坡能否会对已有建构筑物造成损失、确定预警疏散范围的重要依据,因此对地震土坡破坏后的滑坡体大小和致灾范围进行研究具有重要的意义.本研究基于SPH动力分析方法,结合弹塑性本构模型和固体力学控制方程建立了地震土坡破坏的动力分析模型;通过设置振动边界粒子和自由场边界粒子,实现了地震动加速度的施加以及自...     

考虑循环软化的非线性动力本构模型在FLAC3D中的实现

已有震害研究表明,震后边坡会因持续变形而破坏,且伴随着土体强度逐渐降低的现象,即土的循环软化行为。因此,有必要研究考虑循环软化的非线性动力本构模型以用于复杂条件下地震边坡稳定性分析。在已有的非线性动力本构模型基础上,提出了考虑循环软化的处理方法。同时,在FLAC3D平台上实现了本构模型二次开发,并通过了理论公式与已有文献中试验数据的验证。结果表明：计算出的骨干曲线与理论公式一致,且计算出的动剪切模量比及阻尼比与试验数据吻合较好,能够克服Hardin-Drnevich模型和Davidenkov模型在较大应变处（>0.01%）过高地估算阻尼比的缺陷;考虑了循环软化后,计算出的剪切强度有明显降低,且当遇到骨干曲线剪应力可以连续地过渡到软化后的主干曲线上,模型的收敛性较好。所开发的本构模型可为大应变条件下软土场地及边坡地震灾害评估提供支持。     

地震动作用下均匀土质边坡的临界转角分析

均质土坡变形的整体旋转角度（即转角）与其弧形滑动特征相协调,是比永久位移更为优越的土坡变形与稳定性表征量。地震动瞬态弹性变形应变能相当于土坡系统势能的动态随机性“涨落”,临界点（临界转角）作为系统本身的一个固有特性,不受地震动随机性“涨落”的影响,土坡滑动失稳的临界转角实质上对应着系统势能的突变点。基于均质土坡地震动失稳的特点,结合突变理论,给出了确定失稳临界转角的计算公式,发现其仅与滑带介质本构模型参数 、 和滑弧半径R有关,但在本质上体现了诸因素的影响而具有广泛适用性。     

修正有限元强度折减法与失稳判据在边坡稳定分析中的应用

现有的有限元强度折减法用于边坡稳定性分析,在土体变形特性的认识和失稳判据上有缺陷。基于土体的应力-应变曲线和邓肯-张非线性弹性模型进行有限元强度折减法计算分析,在现有方法仅对强度参数--黏聚力c和内摩擦角φ进行折减的基础上,提出了切线弹性模量Et的折减算法,更全面地体现折减的内涵;提出了以潜在滑动区域的能量积分变化作为判据,根据能量积分与折减系数关系曲线的特征,将其拟合成两段直线,以其交点所对应的折减系数作为边坡的安全系数。比较分析了修正Duncan有限元强度折减法、失稳判据与文献的计算结果,同时分析了边坡滑动区能量积分计算范围对结果的影响,并用于南水北调东线实际工程计算分析。结果表明：修正Duncan有限元强度折减法能较全面考虑土体的变形和强度特性,新的边坡滑动失稳判据可行且具有较好的计算稳定性,能保证工程建设的安全。     

考虑滑坡渐进破坏的改进简布条分法

为考虑位移因素对于滑坡渐进性发育的影响,基于简布条分法基本原理,提出一个新的能够考虑滑带土抗剪能力与剪切位移关系的滑坡稳定性简化评价方法。该方法从以下两个方面对简布条分法进行了改进：（1）将滑动面剪应力-剪切位移本构模型与沿滑动面的剪切位移模型相结合,以考虑滑坡随着变形发展渐进破坏的过程;（2）提出局部安全系数的概念并对整体安全系数进行重新定义,物理意义更加明确。推导了该方法的计算公式,给出了其迭代计算思路,并编制Matlab计算程序。通过计算实例表明,应变软化边坡的安全系数不仅取决于土的抗剪强度指标,还与其剪应力-剪切位移关系密切相关,基于峰值强度的极限平衡条分法应用于实际工程是偏于危险的。同时该方法还可以有效地模拟边坡从开始破坏到最终滑动面贯通的发育过程,实现了坡体位移与安全系数之间一一对应的定量关系,进而实现边坡的稳定性的预测预报。     

边坡变形破坏过程的大变形有限元分析

在用有限元法分析边坡稳定性时,引入计算大变形问题的更新的拉格朗日方法,推导了边坡大变形弹塑性有限元分析的方程式。采用边坡某一幅值的等效塑性剪应变区,从坡脚到坡顶贯通前的折减系数作为边坡安全系数。在此基础上,采用弹塑性大变形有限元分析软件计算了均质土坡不同坡角的安全系数,将其与小变形分析的结果进行了对比分析,结果表明：用弹塑性大变形有限元分析边坡失稳破坏的过程中,既考虑了岩土材料的非线性,又考虑了边坡的几何非线性,使计算结果更趋合理。并结合东深供水改造工程BIII2边坡进行了大变形有限元分析,计算结果与勘查到的实际边坡的滑动面分布位置比较接近。研究表明：该方法尤其适宜于软土类边坡或基坑的稳定性分析。     

Failure mechanism of a slope with a thin soft band triggered by intensive rainfall

Rainfall is considered as one of the paramount factors for slope failures in many regions around the world, particularly in tropical and subtropical regions. To study the effect of rainfall storm and its duration on the stability of slopes with a thin soft band layer, a 2D seepage numerical analysis and experimental investigation were implemented on an unsaturated model, consisting of clayey sand soils with a thin soft layer inclined to the horizontal level by 30° at a slope angle of 50°. It was subjected to intensive rainfall 20 mm/h for 8 h. Positive pore water pressures and horizontal earth pressures were monitored during the rainstorm using sensors distributed inside the experimental model. Both the experimental and numerical simulation results showed that the stability of the slope decreased during the time of the rainfall storm. Infiltration of rainwater resulted in reduction of soil shear strength, due to the loss occurring in soil suction after 1 h of rainstorm; the tension cracks appeared at the top of the slope and a certain displacement was observed in the sliding blocks. During the time of rainstorm, the infiltrated water flowed out from the slope through the weak interlayer near to the toe causing piping and local failure, so the formed cracks at the top of slope grew and expanded due to sliding of the failed soil blocks. Moreover, the ground water table rose and the positive pore water pressures increased, resulting in a reduction of effective stress, which is considered as a main factor in soil shear strength. A surface runoff was also present following the full saturation of the slope, leading to dragging the fine particles with water flow causing erosion. The combination of piping and erosion effects led to a quick local failure at the toe, as well as sliding of the failed blocks and spreading of the cracks.     

Constitutive modeling of earthquake-induced slides on clays along slip surfaces

The paper first proposes and validates a constitutive model simulating the change of resistance along a slip surface of clay for both the undrained and drained cases. The proposed model is based on (a) the critical state theory and (b) the assumption that the critical state changes once failure is reached, in terms of the further shear displacement. Under undrained conditions, the proposed model simulates the excess pore pressure generation and, subsequently, the continuous change of resistance along the slip surface from its initial value to the peak strength and then, at large displacement, the residual value. The latter can be measured in constant-volume ring shear tests. Then, the developed constitutive model is implemented in the multi-block sliding system model for the prediction of the triggering and deformation of slides. The improved model is applied at the well-documented Fourth Avenue landslide of the 1964 Alaska earthquake.     

地震作用前后极限平衡法和强度折减法的边坡稳定分析对比

分别运用极限平衡法和ANSYS强度折减法对国内某边坡地震作用前后稳定性进行了分析,考虑地震作用时采用了拟静力法施加水平和垂直两个方向的地震力.根据地震作用前后两种方法得到的边坡稳定系数、滑动面、滑动体的对比分析发现,两种方法得出的安全系数、滑动面相差不大,说明了强度折减法运用到边坡稳定分析中是合理的.强度折减法得出安全系数均比极限平衡法略大,这是由于强度折减法考虑了土体的应力-应变关系,使得计算结果更符合实际.通过地震作用前后边坡的水平位移图的对比,总结了一些规律,为地震区边坡的锚固提出相应的建议.     

Lagrangian meshfree particles method (SPH) for large deformation and failure flows of geomaterial using elastic–plastic soil constitutive model

Simulation of large deformation and post‐failure of geomaterial in the framework of smoothed particle hydrodynamics (SPH) are presented in this study. The Drucker–Prager model with associated and non‐associated plastic flow rules is implemented into the SPH code to describe elastic–plastic soil behavior. In contrast to previous work on SPH for solids, where the hydrostatic pressure is often estimated from density by an equation of state, this study proposes to calculate the hydrostatic pressure of soil directly from constitutive models. Results obtained in this paper show that the original SPH method, which has been successfully applied to a vast range of problems, is unable to directly solve elastic–plastic flows of soil because of the so‐called SPH tensile instability. This numerical instability may result in unrealistic fracture and particles clustering in SPH simulation. For non‐cohesive soil, the instability is not serious and can be completely removed by using a tension cracking treatment from soil constitutive model and thereby give realistic soil behavior. However, the serious tensile instability that is found in SPH application for cohesive soil requires a special treatment to overcome this problem. In this paper, an artificial stress method is applied to remove the SPH numerical instability in cohesive soil. A number of numerical tests are carried out to check the capability of SPH in the current application. Numerical results are then compared with experimental and finite element method solutions. The good agreement obtained from these comparisons suggests that SPH can be extended to general geotechnical problems. Copyright © 2008 John Wiley & Sons, Ltd.     

基于光滑粒子流体动力学的流动法则对土体滑坡大变形的影响探讨

光滑粒子流体动力学（SPH）是近年来发展起来的一种纯拉格朗日无网格法,并因其在大变形领域内的优势而受到广泛关注。在进行滑坡大变形分析时,流动法则及剪胀角的选取对于边坡失稳后的运动特性有重要的影响。采用Fortran自行编写了基于SPH的边坡稳定性及失稳后大变形分析程序,然后通过2个经典的算例,讨论了不同流动法则及剪胀角的选取对滑坡大变形分析精度的影响。结果表明：（1）剪胀角的选取对土体失稳后的滑动距离有显著影响,随着剪胀角的增大,土体的滑动速度与距离呈明显增大趋势;（2）在关联性流动法则及非关联性流动法则? =1/2? （? 为剪胀角,? 为摩擦角）条件下,土体在大变形过程中会产生过度膨胀,且运动速度与距离要大于实际情况,而采用关联性流动法? =0时,对于非膨胀土可以得到比较令人满意的结果,但对于膨胀土体会使得运动速度和运动距离过小,不符合实际情况。建议在计算膨胀土大变形时,采用非关联性流动法则且适当考虑膨胀性（? =(0.1～0.2)? ）,可以得到较好的结果。     

考虑强度退化效应的堤坝抗震稳定性评价方法

考虑土的动强度随振动孔隙压力上升的衰减效应,将拟静力极限平衡分析和滑动体位移分析相结合,提出了堤坝抗震稳定性评价方法.首先,基于土工动力有限元分析,确定坝坡潜在滑动土体的平均加速度时程,进而,基于拟静力概念,采用极限平衡分析,确定坝坡的安全系数随时间的变化历程,其中当安全系数瞬时小于1时,表明坝坡在地震中处于瞬时超载状态,采用Newmark刚体滑块模型估算瞬时超载所产生的滑动位移,将各个超载阶段的滑动位移叠加,求得设计地震动作用下堤坝边坡的累积滑移量,根据这种方法所进行的数值计算与分析表明,考虑强度循环退化效应后所得到的坝坡滑移量更为合理.     

Numerical modeling of seismic triggering,evolution, and deposition of rapid landslides: Application to Higashi–Takezawa (2004)

A mathematical model is developed for the dynamic analysis of earthquake‐triggered rapid landslides, considering two mechanically coupled systems: (a) the accelerating deformable body of the slide and (b) the rapidly deforming shear band at the base of the slide. The main body of the slide is considered as a one‐phase mixture of Newtonian incompressible fluids and Coulomb solids sliding on a plane of variable inclination. The evolution of the landslide is modeled via a depth‐integrated model of the Savage–Hutter type coupled with: (a) a cyclic hysteretic constitutive model of the Bouc–Wen type and (b) Voellmy's rheology for the deformation of the material within the shear band. The original shallow‐water equations that govern the landslide motion are appropriately reformulated to account for inertial forces due to seismic loading, and to allow for a smooth transition between the active and the passive state. The capability of the developed model is tested against the Higashi–Takezawa landslide. Triggered by the 2004 Niigata‐ken Chuetsu earthquake, the slide produced about 100m displacement of a large wedge from an originally rather mild slope. The mechanism of material softening inside the shear band responsible for the surprisingly large run‐out of the landslide is described by a set of equations for grain crushing‐induced pore‐water pressures. The back‐analysis reveals interesting patterns on the flow dynamics, and the numerical results compare well with field observations. It is shown that the mechanism of material softening is a crucial factor for the initiation and evolution of the landslide, while viscoplastic frictional resistance is a key requirement for successfully reproducing the field data. Copyright © 2009 John Wiley & Sons, Ltd.     

Modelling of fibre–cohesive soil mixtures

A new constitutive model for fibre-reinforced cohesive soil is proposed. The model combines a Cam-Clay like bounding surface model with an elastic–plastic one-dimensional fibrous element model. A “smearing procedure”, which can consider any spatial distribution of fibre orientation, is employed to transform discrete tensile forces developed in the fibres into stresses for the composite material. The fibre stress contribution is bounded by both degradation of soil–fibre bonding due to pull-out mechanism and tensile strength of the fibres. Eventual occurrence of fibre breakage is also considered. The model performances are analysed for both consolidation and shearing loading modes, and qualitative comparison is performed with experimental data available in the literature. For consolidation loading, tensile stresses are not developed in the fibres and thus the fibre effect is rather limited. For drained shear loading, addition of fibres can result in a consistent shear strength increase. The beneficial effect of fibres seems to be controlled by two parameters: the fibre tensile stiffness and the fibre/soil strain ratio that accounts for any possible slippage or shear deformation at the fibre/soil matrix interface. For undrained shear loading, the strengthening effect of the fibres appears to be counteracted by the increase in pore water pressure, induced by the additional confining contribution of the fibres. In agreement with published experimental data, the model suggests also that the moisture content is a key factor governing fibre effectiveness for undrained shearing. Finally, analysis of the model predicted critical states for fibre-reinforced cohesive soil is provided.     

地震作用下华丽高速公路金沙江桥华坪岸顺层边坡动力稳定性评价

以华丽金沙江桥华坪岸顺层边坡地震稳定问题为研究对象,在slide下利用Janbu-simplified法计算了其最不利滑动面和安全系数,为模型试验确定了宽900 m×高450 m的试验范围。采用拟静力法,在Phase 2D下利用有限元强度折减法模拟模型边坡,计算了0.00g、0.10g、0.15g、0.20g、0.25g、0.30g、0.35g等7个地震工况的安全系数和位移。对比不同工况下的安全系数、位移与水平地震加速度曲线,结果表明,（1）极限平衡法与强度折减法所得安全系数在1.0附近安全系数离散性最小,远离1.0则安全系数的离散性增大,差异是两类方法极限状态定义的不同造成的。基于安全系数确定的临界水平地震加速度为0.20g;（2）剪应变揭示的边坡破坏模式为桥墩位置拉剪破坏、T2凝灰岩主滑和前缘近于水平剪出破坏,坡体内部的变形响应敏感性为水平位移＞竖向位移＞剪应变。基于变形演化规律确定的临界值为0.15g～0.20g,相应水平位移预警值为10.2 mm,室内模型试验验证了上述分析的合理性;（3）综合确定华坪岸顺层边坡临界水平地震加速度为0.2g,50年超越概率为10%和100年超越概率为2%对应的安全系数分别约为1.2和1.1,华坪岸边坡在现行抗震设防标准下是安全的。当前选择安全系数作为边坡稳定性评价指标,而安全监控实践中一般选择可操作性和预警性强的变形值作为控制指标,所以边坡稳定性评价需要采取安全系数与临界位移相结合的方法。     

海底泥流的流变试验及强度模型

海底泥流是海底斜坡发生失稳滑动后,经复杂的水土交换作用演变成的稀式海底滑坡体,表现出土体和流体的双重特征,而现有的测试方法难以获得低剪切应变率下连续变化的强度值,不能很好地揭示其综合强度特性。利用新型全流动贯入仪和RST流变仪,对模拟海底泥流在不同剪切应变率下的流变和强度特性进行了多组试验研究,并基于试验结果分析了不排水剪切强度、屈服应力和表观黏度与含水率的相关性。基于剪切稀化理论,提出分阶段拟合模式来描述海底泥流从低剪切应变率到中高剪切应变率的流变关系;通过多种常规流变模式拟合结果的对比分析,显示出新流变模式的适用性和优势。考虑强度软化的影响,建立了全剪切应变率范围内海底泥流的不排水剪切强度模型,为海底泥流运动过程的数值模拟和灾害评价提供科学依据。     

A coupled two‐phase fluid flow and elastoplastic deformation model for unsaturated soils: theory,implementation, and application

Although numerous numerical models have been proposed for simulating the coupled hydromechanical behaviors in unsaturated soils, few studies satisfactorily reproduced the soil–water–air three‐phase coupling processes. Particularly, the impacts of deformation dependence of water retention curve, bonding stress, and gas flow on the coupled processes were less examined within a coupled soil–water–air model. Based on our newly developed constitutive models (Hu et al., 2013, 2014, 2015) in which the soil–water–air couplings have been appropriately captured, this study develops a computer code named F²Mus3D to investigate the coupled processes with a focus on the above impacts. In the numerical implementation, the generalized‐α time integration scheme was adopted to solve the equations, and a return‐mapping implicit stress integration scheme was used to update the state variables. The numerical model was verified by two well‐designed laboratory tests and was applied for modeling the coupled elastoplastic deformation and two‐phase fluid flow processes in a homogenous soil slope induced by rainfall infiltration. The simulation results demonstrated that the numerical model well reproduces the initiation of a sheared zone at the toe of the slope and its propagation toward the crest as the rain infiltration proceeds, which manifests a typical mechanism for rainfall‐induced shallow landslides. The simulated plastic strain and deformation would be remarkably underestimated when the bonding stress and/or the deformation‐dependent nature of hydraulic properties are ignored in the coupled model. But on the contrary, the negligence of gas flow in the slope soil results in an overestimation of the rainfall‐induced deformation. Copyright © 2015 John Wiley & Sons, Ltd.     

液化型路堤边坡动力数值模拟分析

液化型路堤边坡动力稳定性问题涉及岩土工程与工程地震两个学科领域,是边坡工程与砂土液化的交叉课题。采用天然地震记录为输入条件,应用Finn本构关系模型,运用有限差分法,对填土+砂土+卵砾土地层组合的路堤边坡进行了全时程动力分析,探讨了地震作用下路堤边坡的液化初步规律和稳定性。数值模拟结果表明:地震作用引起了路基饱和砂土有效应力急剧减小,并导致路基砂土液化,引起路堤变形破坏。孔隙水压力的积累与消散不仅与地震记录序列存在对应关系,也与砂土所处的位置和深度有密切关系。地表变形破坏主要表现为路堤顶面发生震陷和拉裂破坏,坡底面产生挤压隆起变形。地面以下的变形破坏主要包括土体剪切破坏和深部砂土液化引起的侧向流动破坏。     

膨胀土边坡非饱和渗流及渐进性破坏耦合分析

膨胀土边坡受降雨影响产生膨胀变形,是典型的非饱和土多场耦合问题。为探究降雨入渗对其渐进性破坏的失稳过程,基于饱和-非饱和渗流理论、膨胀土弹塑性本构关系和应变软化理论,利用应变软化模型、FLAC3D二次开发平台和内置FISH语言,提出了一种综合考虑非饱和渗流、膨胀变形和应变软化的多场耦合数值分析法。结合工程实例,通过该方法探讨了降雨入渗条件下膨胀土边坡非饱和渗流、位移响应及渐进性破坏的变化规律。结果表明:膨胀变形和应变软化受控于非饱和渗流的时空分布,对边坡位移响应过程影响显著,也易导致饱和-非饱和分界带形成剪应力集中区。膨胀土边坡渐进性破坏由局部破坏转变为整体性失稳,其塑性破坏区首先随悬挂型暂态饱和区的变化向坡内扩展,雨后逐渐形成第二条由坡脚向坡顶扩展的滑动带,呈现出多重滑动性和后退牵引式的破坏特征。