Fully coupled elastoplastic hydro-mechanical analysis of unsaturated porous media using a meshfree method

This paper presents the first application of an advanced meshfree method, ie, the edge-based smoothed point interpolation method (ESPIM), in simulation of the coupled hydro-mechanical behaviour of unsaturated porous media. In the proposed technique, the problem domain is spatially discretised using a triangular background mesh, and the polynomial point interpolation method combined with a simple node selection scheme is adopted for creating nodal shape functions. Smoothing domains are formed on top of the background mesh, and a constant smoothed strain, created by applying the smoothing operation over the smoothing domains, is assigned to each smoothing domain. The deformation and flow models are developed based on the equilibrium equation of the mixture, and linear momentum and mass balance equations of the fluid phases, respectively. The effective stress approach is followed to account for the coupling between the flow and deformation models. Further coupling among the phases is captured through a hysteretic soil water retention model that evolves with changes in void ratio. An advanced elastoplastic constitutive model within the context of the bounding surface plasticity theory is employed for predicting the nonlinear behaviour of soil skeleton. Time discretisation is performed by adopting a three-point discretisation method with growing time steps to avoid temporal instabilities. A modified Newton-Raphson framework is designed for dealing with nonlinearities of the discretised system of equations. The performance of the numerical model is examined through a number of numerical examples. The state-of-the-art computational scheme developed is useful for simulation of geotechnical engineering problems involving unsaturated soils.     

Subsurface flow measurements using passive flux meters in variably-saturated cold-regions landscapes

To date, passive flux meters have predominantly been applied in temperate environments for tracking the movement of contaminants in groundwater. This study applies these instruments to reduce uncertainty in (typically instantaneous) flux measurements made in a low-gradient, wetland dominated, discontinuous permafrost environment. This method supports improved estimation of unsaturated and over-winter subsurface flows which are very difficult to quantify using hydraulic gradient-based approaches. Improved subsurface flow estimates can play a key role in understanding the water budget of this landscape.     

Revisiting the existence of an effective stress for wet granular soils with micromechanics

A possible effective stress variable for wet granular materials is numerically investigated based on an adapted discrete element method (DEM) model for an ideal three‐phase system. The DEM simulations consider granular materials made of nearly monodisperse spherical particles, in the pendular regime with the pore fluid mixture consisting of distinct water menisci bridging particle pairs. The contact force‐related stress contribution to the total stresses is isolated and tested as the effective stress candidate for dense or loose systems. It is first recalled that this contact stress tensor is indeed an adequate effective stress that describes stress limit states of wet samples with the same Mohr‐Coulomb criterion associated with their dry counterparts. As for constitutive relationships, it is demonstrated that the contact stress tensor used in conjunction with dry constitutive relations does describe the strains of wet samples during an initial strain regime but not beyond. Outside this so‐called quasi‐static strain regime, whose extent is much greater for dense than loose materials, dramatic changes in the contact network prevent macroscale contact stress‐strain relationships to apply in the same manner to dry and unsaturated conditions. The presented numerical results also reveal unexpected constitutive bifurcations for the loose material, related to stick‐slip macrobehavior.     

Finite element formulation and algorithms for unsaturated soils. Part II: Verification and application

The finite‐element formulation and integration algorithms developed in Part I are used to analyse a number of practical problems involving unsaturated and saturated soils. The formulation and algorithms perform well for all the cases analysed, with the robustness of the latter being largely insensitive to user‐defined parameters such as the number of coarse time steps and error control tolerances. The efficiency of the algorithms, as measured by the CPU time consumed, does not depend on the number of coarse time steps, but may be influenced by the error control tolerances. Based on the analyses presented here, typical values for the error control tolerances are suggested. It is also shown that the constitutive modelling framework presented in Part I can, by adjusting one constitutive equation and one or two material parameters, be used to simulate soils that expand or collapse upon wetting. Treating the suction as a strain variable instead of a stress variable proves to be an efficient and robust way of solving suction‐dependent plastic yielding. Moreover, the concept of the constitutive stress is a particularly convenient way of handling the transition between saturation and unsaturation. Copyright © 2003 John Wiley & Sons, Ltd.     

Numerical modeling of rainstorm-induced shallow landslides in saturated and unsaturated soils

For the assessment of shallow landslides triggered by rainfall, the physically based model coupling the infinite slope stability analysis with the hydrological modeling in nearly saturated soil has commonly been used due to its simplicity. However, in that model the rainfall infiltration in unsaturated soil could not be reliably simulated because a linear diffusion-type Richards’ equation rather than the complete Richards’ equation was used. In addition, the effect of matric suction on the shear strength of soil was not actually considered. Therefore, except the shallow landslide in saturated soil due to groundwater table rise, the shallow landslide induced by the loss in unsaturated shear strength due to the dissipation of matric suction could not be reliably assessed. In this study, a physically based model capable of assessing shallow landslides in variably saturated soils is developed by adopting the complete Richards’ equation with the effect of slope angle in the rainfall infiltration modeling and using the extended Mohr–Coulomb failure criterion to describe the unsaturated shear strength in the soil failure modeling. The influence of rainfall intensity and duration on shallow landslide is investigated using the developed model. The result shows that the rainfall intensity and duration seem to have similar influence on shallow landslides respectively triggered by the increase of positive pore water pressure in saturated soil and induced by the dissipation of matric suction in unsaturated soil. The rainfall duration threshold decreases with the increase in rainfall intensity, but remains constant for large rainfall intensity.     

强降雨条件下土质边坡瞬态稳定性分析

基于实际降雨气象资料,设计了单峰降雨和8个不同时间间隔的双峰降雨计算方案,利用非饱和土力学理论,对边坡的瞬态稳定性进行了计算和分析,研究了水分在坡体内的运移对边坡稳定性的时间空间影响效应,同时考察了降雨入渗造成的土性渗透特性的改变。分析发现：一次降雨的影响历时约12 d,降雨后约0.5 d该类土质边坡最危险;不同时间间隔的两次降雨对边坡稳定的影响比单峰降雨的最小安全系数滞后约0.3 d～0.8 d,影响历时基本保持不变;土体物理力学性质、边坡最危险滑动面及其对应的安全系数随水分在坡体内渗透运移而变化。     

新疆巴里坤膨润土的工艺应用矿物学特性研究

从工艺矿物学角度系统研究了新疆巴里坤膨润土矿的矿石特征、化学成分、物理性质、矿物组成与结构、工艺应用属性等,认为巴里坤膨润土矿蒙脱石相对含量高,原土适用于石油钻井、铸造和冶金工业,将膨润土深加工后,可制得符合市场需求的有机膨润土、无机凝胶和高纯度蒙脱石等深加工产品.     

非饱和膨胀土边坡的稳定性分析

为了提高非饱和膨胀土边坡的稳定,采用了土工膜覆盖边坡表面的方法.基于非饱和渗流理论,采用有限元程序对用膜覆盖和不覆盖的边坡分别进行了降雨蒸发干湿循环条件下的数值分析.在计算结果基础上,采用极限平衡法分析了边坡的稳定性.结果表明,用膜覆盖后边坡能长期保持非饱和状态,土体吸力较高,其稳定安全系数明显高于无膜保护的边坡,用复合土工膜覆盖边坡有利于提高非饱和膨胀土边坡的稳定.     

刘房子钠基膨润土矿石质量研究及应用

刘房子钠基膨润土矿床是我国膨润土矿石质量优异的矿区,经过三年的详查与勘探,探明钠基膨润土矿石2 300万t,钠蒙脱石含量高可达98%,平均含量达63.74%。本文通过X射线衍射、红外吸收光谱分析、矿物的粒度分析,并结合电子扫描对矿石质量进行了研究,应用热分析方法,确定刘房子钠基膨润土矿石具有的优良特征。     

Design of Inclined Covers with Capillary Barrier Effect

A design procedure is proposed to minimize water infiltration into landfills by optimizing the water diversion length of inclined covers with capillary barrier effect (CCBE). This design procedure is based on a conceptual, mathematical and numerical approach and aims at selecting materials and optimizing layer thickness. Selection among candidate materials is made based on their hydraulic conductivity functions and on a threshold infiltration rate imposed on the designer. The capillary break layer (CBL; bottom layer) is characterized by a weak capillarity, while the moisture retention layer (MRL; upper layer) is characterized by a compromise between strong capillarity and high hydraulic conductivity. The thickness of the CBL corresponds to the height where suction reaches its maximum value for a given infiltration rate. This height can be calculated using the Kisch [Géotechnique 9 (1959)] model. The optimal thickness of the MRL is determined by applying an adaptation of the Ross [Water Resources Research 26 (1990)] model. The results obtained using the proposed design procedure were compared to those obtained from numerical simulations performed using a finite element unsaturated seepage software. The procedure was applied for two cover systems; one where deinking by-products (DBP) were used as MRL and sand as CBL and another where sand was used as MRL and gravel as CBL. Using this procedure, it has been shown that an infiltration control system composed of thin layers of sand over gravel is highly efficient in terms of diversion length and that its efficiency can be enhanced by placing a hydraulic barrier – such as a layer of DBP – above the MRL.