降雨诱发浅层滑坡稳定性的计算模型研究

我国是一个滑坡灾害频发的国家,众多事实表明：降雨是影响边坡稳定性,导致边坡失稳的最主要和最普遍的环境因素,是浅层滑坡的触发因素。为了更好地对降雨诱发浅层滑坡进行研究,采用非饱和土VG模型与改进的Green-Ampt入渗模型对Mein-Larson降雨入渗模型进行改进,并结合无限边坡提出了一个降雨诱发浅层滑坡的简化计算模型。与以往提出的简化计算模型相比,该模型既考虑了坡面倾斜的影响,又考虑了非饱和土的特性,并可用于两种降雨形式下的边坡浅层稳定性估算,具有更广的应用范围。通过与有限元得到的结果进行比较可得：在不同降雨条件下,该计算模型得到的各项结果与数值解是接近的,安全系数计算结果是偏于安全的,因此,可将该计算模型用于降雨诱发浅层滑坡的近似估算;该计算模型公式简单,便于计算,计算效率较高。     

Generalized plasticity state parameter‐based model for saturated and unsaturated soils. Part II: Unsaturated soil modeling

The aim of this paper is to extend the generalized plasticity state parameter‐based model presented in part 1 to reproduce the hydro‐mechanical behavior of unsaturated soils. The proposed model is based on two pairs of stress–strain variables and a suitable hardening law taking into account the bonding—debonding effect of suction and degree of saturation. A generalized state parameter for unsaturated state is proposed to reproduce soil behavior using a single set of material parameters. Generalized plasticity gives a suitable framework to reproduce not only monotonic stress path but also cyclic behavior. The hydraulic hysteresis during a drying—wetting cycle and the void ratio effect on the hydraulic behavior is introduced. Comparison between model simulations and a series of experimental data available, both cohesive and granular, are given to illustrate the accuracy of the enhanced generalized plasticity equation. Copyright © 2010 John Wiley & Sons, Ltd.     

Experimental analysis and modelling of shallow landslides

The paper presents the results of some experimental tests reproducing the triggering mechanism of a special kind of shallow landslides induced by rainfalls (soil slip) in a physical 1-g model. The experimental data have been employed to verify the capability of a simplified stability model to describe the phenomenon and to back-analyse its occurrence in a case history (Pizzo d’Alvano, Campania Region 1998). The method enables a direct correlation between the safety factor of a slope and rainfall intensity, as well as antecedent rainfalls.     

The influence of the geological and geomorphological settings on shallow landslides. An example in a temperate climate environment: the June 19, 1996 event in northwestern Tuscany (Italy)

On June 19, 1996, an extremely heavy rainstorm hit a restricted area in the Apuan Alps (northwestern Tuscany, Italy). Its max intensity concentrated over an area of about 150 km² astride the Apuan chain, where 474 mm was recorded in about 12 h (21% of the mean annual precipitation, with an intensity up to 158 mm/h). The storm caused floods and hundreds of landslides and debris flows, which produced huge damage (hundreds of millions of Euros), partially destroyed villages and killed 14 people. This paper reports the results obtained from a detailed field survey and aerial view interpretation. In the most severely involved area, 647 main landslides were investigated, mapped and related to the geologic, geomorphic and vegetational factors of the source areas. This was in order to define the influence of these factors and contribute to an evaluation of the landslide hazard in the study area. An assessment was also made of the total area and volume of material mobilised by landsliding. The study area, about 46 km² wide, includes three typically mountainous basins, characterised by narrow, deep cut valleys and steep slopes, where many rock types outcrop. Most of the landslides were shallow and linear, referable to complex, earth and debris translational slide, which quickly developed into flow (soil slip–debris flow). Usually, they involved colluvium and started in hollows underlain by metamorphic rock (metasandstone and phyllite), often dipping downslope. Therefore, bedrock lithology and impermeability appeared to be important factors in the localisation of the landslide phenomena. The investigation of the geomorphic and land use features in the source areas also frequently highlighted a rectilinear profile of the slope, a high slope gradient (31–45°) and dense chestnut wood cover. In the area, about 985,000 m² (2.1% of 46 km²) was affected by landsliding and about 700,000 m² of this area was covered by chestnut forest. The landslides removed about 7000 trees. The volume of mobilised material was about 1,360,000 m³; about 220,000 m³ remained on the slopes, while the rest poured into the streams. In addition, about 945,000 m³ was mobilised by the torrential erosion in the riverbeds.     

Numerical analysis of wetting-induced instabilities in partially saturated soils under plane strain conditions

This paper concerns the hydromechanical stability of partially saturated soils during wetting. The response of a loose silty sand is numerically/theoretically investigated with the main aim of identifying both triggering mechanisms and predisposing factors to instability. This latter is testified by both a rapid increase in the pore water pressure and an unexpected loss of numerical convergence. The study has been conducted at different scales from laboratory tests to boundary value problems, and in both cases, the controllability theory has been employed to both interpret and justify the numerical results.     

Modelling of shallow landslides with machine learning algorithms

This paper introduces three machine learning(ML)algorithms,the‘ensemble'Random Forest(RF),the‘ensemble'Gradient Boosted Regression Tree(GBRT)and the Multi Layer Perceptron neural network(MLP)and applies them to the spatial modelling of shallow landslides near Kvam in Norway.In the development of the ML models,a total of 11 significant landslide controlling factors were selected.The controlling factors relate to the geomorphology,geology,geo-environment and anthropogenic effects:slope angle,aspect,plan curvature,profile curvature,flow accumulation,flow direction,distance to rivers,water content,saturation,rainfall and distance to roads.It is observed that slope angle was the most significant controlling factor in the ML analyses.The performance of the three ML models was evaluated quantitatively based on the Receiver Operating Characteristic(ROC)analysis.The results show that the‘ensemble'GBRT machine learning model yielded the most promising results for the spatial prediction of shallow landslides,with a 95%probability of landslide detection and 87%prediction efficiency.     

非饱和土中波的传播特性

基于已建立的非饱和土中的波动方程,第1次从理论上推证出非饱和土中存在4种体波,即1种S波和3种P波（P1波、P2波和P3波）,并导出了4种体波的波速以及衰减的解析表达式,分析了4种体波的波速、衰减与频率以及饱和度之间的关系。结果表明：3种P波中P3波波速最小,衰减最大;P1波波速最大,衰减最小;P2波则介于二者之间。4种体波都具有弥散性,一般而言,波速随频率的提高而有所减小,而衰减则有所增大。随含气量的提高,P1波的波速减小,P2波和S波的波速增大,而P3波则受含气量的影响很小;P1波和P3波的衰减都有所减小,但变化不是很显著,P2波的衰减基本不随含气量的变化而变化,S波几乎不存在衰减现象。     

Equivalent stress approach in modelling unsaturated soils

This article presents an equivalent stress approach that can be used in many elastoplastic constitutive models for unsaturated soils. The use of the equivalent stress leads to a modified yield locus that is independent of the suction. In addition, the equivalent stress becomes the major stress variable, with suction required only as an additional variable in calculations. The model on the basis of equivalent stress predicts exactly the same soil behaviour, with the sole difference being the use of equivalent stress instead of original stress variables. This article also presents the equivalent stress formulations of several constitutive models for unsaturated soils, including the Barcelona Basic Model. The predictions from these models remain unchanged, with the only difference being in their implementation. Finally, the equivalent stress approach and the net stress approach are compared for the Barcelona Basic Model. Copyright © 2011 John Wiley & Sons, Ltd.     

非饱和土抗剪强度理论的研究进展

非饱和土的工程性质是20世纪90年代以来国际岩土工程界研究的热点,非饱和土的抗剪强度是非饱和土的重要工程性质之一,在参阅了中、外文相关资料的基础上,总结了非饱和土抗剪强度理论的研究成果,提出了非饱和土抗剪强度有待进一步研究的问题。     

非饱和土固结理论新进展

总结了国内外多年来在非饱和土固结理论领域研究的主要理论成果,内容包括国内外非饱和土固结理论研究的概况,涉及线性、非线性、弹塑性和结构性等方面,非饱和土固结理论的研究具有非常好的应用前景。但非饱和土固结理论还很不成熟,有待广大科研工作者的进一步研究。     

铜离子在非饱和土中迁移的离心模型试验研究

利用土工离心机考察了铜离子泄漏后6个月至3年的时间里非饱和土中铜离子的变化和污染物的迁移机理。通过4种黏性上的离心模型试验研究，总结出铜离子在这类土体中的扩散规律。铜离子在一年内，主要受重力影响和水力梯度的作用，随渗流运移，存在一个舌状竖向扩散峰，污染范围近似呈三角形；一年后扩散范围逐渐转变为近似的正方形。竖向扩散速度高于水半向。土体密实度对铜离予的扩散有较为明显的影响，有效的击实或者碾压有助下防止污染物迁移。黏粒含量对铜离子的迁移阻碍非常迅速而且有效，黏粒含量很小的变化都会对迁移范围和速度有较大的影响。     

吸力对非饱和土抗剪强度的贡献

回顾了非饱和土抗剪强度理论，并进行了评价；系统地分析和总结了吸力对抗剪强度的贡献。     

Work input for unsaturated soils considering interfacial effects

In this research, the interfacial energy is taken into account in the deformation work for unsaturated soils. Based on porous media theory, the thermodynamic balance equations for each phase and the interface are used to derive the work input for unsaturated soils. The work input equation serves as the basis and starting point for the choice of stress state variables, based on which the conjugate stresses and strain increments are derived. The influences of the interfaces on the effective stress and the constitutive law for the liquid phase are then discussed based on the work input equation. The effective stress can be expressed as Bishop's type, and the effective stress parameter is shown to be a function of both the degree of saturation and the interfacial area. The constitutive law for the liquid phase under dynamic condition is also presented. The relationship among interfacial area, saturation, and capillary pressure is proposed to calculate the value of the effective stress. Experimental data obtained from literature are used to validate the proposed model equations. Results show that our findings are in accordance with the existing research. Unlike the phenomenal study, our research has a rigorous theoretical basis, which lays a foundation for further research of unsaturated soils considering the interfacial effects.     

An evolutionary approach to modelling the thermomechanical behaviour of unsaturated soils

A new data‐mining approach is presented for modelling of the stress–strain and volume change behaviour of unsaturated soils considering temperature effects. The proposed approach is based on the evolutionary polynomial regression (EPR), which unlike some other data‐mining techniques, generates a transparent and structured representation of the behaviour of systems directly from raw experimental (or field) data. The proposed methodology can operate on large quantities of data in order to capture nonlinear and complex relationships between contributing variables. The developed models allow the user to gain a clear insight into the behaviour of the system. Unsaturated triaxial test data from the literature were used for development and verification of EPR models. The developed models were also used (in a coupled manner) to produce the entire stress path of triaxial tests. Comparison of the EPR model predictions with the experimental data revealed the robustness and capability of the proposed methodology in capturing and reproducing the constitutive thermomechanical behaviour of unsaturated soils. More importantly, the capability of the developed models in accurately generalizing the predictions to unseen data cases was illustrated. The results of a sensitivity analysis showed that the models developed from data are able to capture and represent the physical aspects of the unsaturated soil behaviour accurately. The merits and advantages of the proposed methodology are also discussed. Copyright © 2014 John Wiley & Sons, Ltd.     

非饱和土土-水特征曲线的温度效应

基于热力学的理论,利用van Genuchten的土-水特征曲线表达式,建立了一种能考虑温度影响的土-水特征曲线方程,并给出了推导过程。该表达式综合考虑了温度对表面张力和浸润系数的影响。相对于根据试验数据拟合的表达式,该表达式具有更加坚实的理论基础,更具适用性和一般性。选用MX-80斑脱土和黄土土样的试验结果,将该方程预测结果和试验结果进行了比较;可以看到,所提出的方程能较好地反映非饱和土土-水特征曲线随温度而变化的情况。     

不同土壤转换函数预测砂土非饱和导水率的对比分析

以现有的4种土壤转换函数(Pedotransfer Functions, PTFs)(包括ROSETTA、RAWLS、CAMPBELL和VAUCLIN)为例,将其用于土壤水力性质数据库UNSODA中56个典型砂土样本的非饱和导水率预测,对4种PTFs的预测结果进行了误差分析,探讨了饱和导水率对非饱和导水率预测的影响。结果表明,当同种PTFs预测的饱和导水率作为输入时,在4种PTFs中VAUCLIN的预测精度最高,其次为ROSETTA和RAWLS,CAMPBELL的误差最大;当实测饱和导水率作为输入时,RAWLS、CAMPBELL和ROSETTA的预测精度有了不同程度提高,但VAUCLIN的预测误差反而有所增大;饱和导水率对非饱和导水率预测的影响较大,实测饱和导水率作为输入时CAMPBELL预测的非饱和导水率与实测值最为接近。     

Three-dimensional multilayer cylindrical tunnel model for calculating train-induced dynamic stress in saturated soils

This study proposes an improved tunnel model for evaluating train-induced dynamic stress in saturated soils, which can consider multiple moving loads, grouting layer and pore-water pressure. Using Shanghai Metro’s actual parameters for train speed, tunnel, grouting layer and soils, the analysis of the spatial distribution of dynamic stress for soils and stress state of various locations under moving train loads shows that neglecting effects such as pore-water pressure can lead to underestimating dynamic normal stress and overestimating dynamic shear stress in the soils below tunnel. This model can be further extended to investigate principal stress axes rotations and tunnel settlement.     

非饱和土的应力状态变量研究

以热力学和连续介质力学的应力理论为基础,从非饱和土的土骨架平衡方程出发,导出了一系列应力状态变量,目前文献中所使用的应力状态变量都在其中。研究发现：非饱和土的应力状态变量有多种组合形式（5组以上）,不是惟一的;大多数应力状态变量与土的孔隙率、饱和度等物性指标孪生相伴;所有应力状态变量都具有应力的量纲,单纯孔隙率或饱和度及其组合都不构成应力状态变量;简化的Bishop有效应力公式和Fredlund倡导的两个应力状态变量都是简化了的应力状态变量,前者忽略了孔隙率的影响,而后者忽略了与孔隙率及饱和度的关联性。在选择应力状态变量时应考虑理论合理、应用方便、逻辑关系正确的原则,还应做一定的验证工作。     

非饱和土水-力本构模型及其隐式积分算法

在已有工作基础上建立了水力-力学耦合的非饱和土本构模型,在硬化方程中考虑饱和度的影响,同时在土水特征曲线中考虑了塑性体变的影响,从而使模型可以反映非饱和土中的毛细现象与土中弹塑性变形现象的耦合行为。采用隐式积分方法,建立了非饱和土耦合模型的数值模型,并推导了得到了水力-力学耦合的非饱和土的一致切线模量。利用该算法编制了本构模型计算的子程序,使其能向外输出切线刚度矩阵,用于有限元计算。为了验证该算法和程序的正确性,用所编制程序对不同路径下的土体行为进行了预测。通过预测结果与试验结果相对比,表明程序预测结果与试验数据相吻合,模型可以较好地模拟土体的水力-力学耦合行为特性。     

泄滩滑坡碎块石土饱和与非饱和水力学参数的现场试验研究

碎块石土由于块石含量较高,块石粒径较大,其水力学参数的确定具有一定困难。首先,采用双套环法对三峡库区泄滩滑坡的滑体碎块石土饱和渗透系数进行了原位试验,并根据土层孔隙率、颗粒级配等因素采用相关经验公式对试验结果进行了分析。其次,结合使用张力计和体积含水率仪对其土水特征曲线进行了现场模拟试验,并采用Fredlund模型对试验结果进行了拟合分析。最后,根据土水特征曲线和饱和渗透系数,采用经验公式估算了其非饱和渗透系数。试验及分析表明,该滑坡的碎块石土层的饱和渗透系数为(1.78～3.2)×10-2 cm/s,为强渗透性;材料的细颗粒含量越少,有效粒径及控制粒径越大,不均匀系数越小,相应的渗透系数越大。相关研究成果可以为泄滩滑坡非饱和非稳定渗透计算提供参数依据,并对同类型土体非饱和水力学参数的确定具有一定的借鉴意义。