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

Omid Ghaffaripour  Golnaz A. Esgandani  Arman Khoshghalb  Babak Shahbodaghkhan 《国际地质力学数值与分析法杂志》2019,43(11):1919-1955

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     

Élise G. Devoie  Ryan F. Connon  James R. Craig  William L. Quinton 《水文研究》2020,34(23):4541-4546

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          下载免费PDF全文

Jérôme Duriez  Richard Wan  Mehdi Pouragha  Félix Darve 《国际地质力学数值与分析法杂志》2018,42(8):959-978

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     

Daichao Sheng  David W. Smith  Scott W. Sloan  Antonio Gens 《国际地质力学数值与分析法杂志》2003,27(9):767-790

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     

Tung-Lin Tsai  Hung-En Chen  Jinn-Chuang Yang 《Environmental Geology》2008,55(6):1269-1277

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.  相似文献   

强降雨条件下土质边坡瞬态稳定性分析   总被引：9，自引：1，他引：8  

吴长富  朱向荣  尹小涛  刘雪梅 《岩土力学》2008,29(2):386-391

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

非饱和膨胀土边坡的稳定性分析   总被引：2，自引：0，他引：2  

吴珺华  YUAN Jun-ping  袁俊平  卢廷浩 《岩土力学》2008,29(Z1):363-367

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

Design of Inclined Covers with Capillary Barrier Effect     

Serge-Étienne Parent  Alexandre Cabral 《Geotechnical and Geological Engineering》2006,24(3):689-710

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.  相似文献   

地下非均匀非饱和带中地下洞室的渗流问题数值模拟--介质参数的灵敏度分析   总被引：3，自引：0，他引：3  

李国敏  Chin-FuTsang 《地球科学》2003,28(5):497-504

影响任何一个地下废物处置场所长期行为的一个重要因素是进入地下废物贮置洞室内的渗流量的大小，预测地下洞室中的渗流量是困难的，特别是当地下废物贮置洞室位于非均质非饱和带中，三维数值模型用于研究地下非均质中非饱和流动及流入地下洞室中的渗流量，讨论了非均质模型与均质模型的比较以及数值剖分尺度对计算结果的影响，进入洞室中的渗漏率随着引入系统中入渗量的增加而增大，选用4个参数来衡量渗流场的非均质程度：(1)介质的平均渗透率K0；(2)Van Genuchten参数α；(3)渗透率空间分布相关尺度；(4)渗透率空间分布变化的标准差σ，根据一个随机实现的渗透率分布，通过改变平均渗透率来研究其对流入洞室中渗透量的影响，对一个固定的入渗率而言，流入洞室中的渗漏率将随着Van Genuchten参数α的减小而减小，模拟结果表明流入洞室的渗漏率与介质的平均渗透率相关，即随平均渗透率的增大，流入洞室的渗漏率亦增大，流入洞室的渗漏率还高度依赖于非均质渗透率场的空间分布相关长度与标准差，一个大的相关长度或高的标准差均能导致流入洞室的渗漏率增大。  相似文献   

Mechanisms of riverbank failure along the Arno River,central Italy     

Stefano Dapporto  Massimo Rinaldi  Nicola Casagli  Pietro Vannocci 《地球表面变化过程与地形》2003,28(12):1303-1323

Riverbanks along the Arno River have been investigated with the aims of de?ning the main mechanisms of failure and retreat, their spatial distribution, and their causes. Geomorphological aspects were investigated by a reconnaissance of riverbank processes, for a number (26) of representative sites. Laboratory and in situ tests were then performed on a selected number of riverbanks (15). Based on the material characteristics, six main typologies of riverbanks have been de?ned, with homogeneous ?ne‐grained and composite banks representing the most frequent types. Slab‐type failures are the most frequent mechanism observed on ?ne‐grained banks, while cantilever failures prevail on composite banks. The role of river stage and related pore water pressure distributions in triggering the main observed mechanisms of failure has been investigated using two different types of stability analysis. The ?rst was conducted for 15 riverbanks, using the limit equilibrium method and considering simpli?ed hypotheses for pore water pressure distribution (annulment of negative pore pressures in the portion of the bank between low water stage and peak stage). Stability conditions and predicted mechanisms of failure are shown to be in reasonably good agreement with ?eld observations. Three riverbanks, representative of the main alluvial reaches of the river, were then selected for a more detailed bank stability analysis, consisting of: (a) de?nition of characteristic hydrographs of the reach with different return periods; (b) modelling of saturated and unsaturated ?ow using ?nite element seepage analysis; and (c) stability analysis with the limit equilibrium method, by adopting pore water pressure values derived from the seepage analysis. The results are compared to those obtained from the previous simpli?ed analysis, and are used to investigate the different responses, in terms of stability, to different hydrological and riverbank conditions. Copyright © 2003 John Wiley & Sons, Ltd.  相似文献