粘土干湿循环中裂缝演变过程的数值模拟

在非饱和土简化固结理论基础上,对粘土在干湿过程中,表面裂缝从形成、发展到闭合的发展过程进行了数值模拟,模拟结果说明,本文采用的基于非饱和简化固结理论的计算方法是可行、实用的.     

Modeling 3-D desiccation soil crack networks using a mesh fragmentation technique

The problem of desiccation cracks in soils has received increasing attention in the last few years, in both experimental investigations and modeling. Experimental research has been mainly focused on the behavior of slurries subjected to drying in plates of different shapes, sizes and thickness. The main objectives of these studies were to learn about the process of crack formation under controlled environmental conditions, and also to understand better the impact of different factors (e.g. soil type, boundary conditions, soil thickness) on the morphology of the crack network. As for the numerical modeling, different approaches have been proposed to describe the behavior of drying cracks in soils. One aspect that it is still difficult to simulate properly is the 3-D crack pattern typically observed in desiccated soils. In this work we present a numerical technique to model the behavior of drying soils. The proposed approach inserts high aspect ratio elements in-between standard elements of a finite element mesh. This mesh fragmentation technique can be easily adapted to standard finite element programs. We used this technique to analyze multiple case studies related to soil desiccation cracks developed under laboratory and field conditions. We focused our attention in some key factors that control the 3-D morphology of the drying cracks network in soils. We show that the proposed technique is able to simulate very satisfactorily the main patterns typically observed in cracked soils.     

Experimental Set-up for Determining Soil Water Retention Curves for Granular Soils During Drying

Abstract: Soil water retention curves (SWRCs) provide an important means of describing the response of unsaturated soils during drying / wetting processes in terms of variations of degree of saturation, water content or void ratio with suction. A key consideration in generating these curves is how to measure the suction. Frequently the filter paper technique is adopted, especially when high suctions are developed, e.g. with plastic clays. As each measurement takes at least a week with this technique, it can take months or years to generate a full SWRC in drying and wetting. Developments in laboratory tensiometers now allow matrix suctions up to about 1.5 MPa to be measured. With such a device it is possible to develop SWRCs for granular soils such as silts and clays in hours or days by continuous measurement. This paper describes an experimental set-up that was developed to measure changes in volume, water content and matrix suction during drying of three granular soils. Limitations of the apparatus and usefulness of the curves are discussed.     

Numerical analysis of seepage–deformation in unsaturated soils

A coupled elastic–plastic finite element analysis based on simplified consolidation theory for unsaturated soils is used to investigate the coupling processes of water infiltration and deformation. By introducing a reduced suction and an elastic–plastic constitutive equation for the soil skeleton, the simplified consolidation theory for unsaturated soils is incorporated into an in-house finite element code. Using the proposed numerical method, the generation of pore water pressure and development of deformation can be simulated under evaporation or rainfall infiltration conditions. Through a parametric study and comparison with the test results, the proposed method is found to describe well the characteristics during water evaporation/infiltration into unsaturated soils. Finally, an unsaturated soil slope with water infiltration is analyzed in detail to investigate the development of the displacement and generation of pore water pressure.     

A simple analytical solution to one‐dimensional consolidation for unsaturated soils

This paper presents a simple analytical solution to Fredlund and Hasan's one‐dimensional (1‐D) consolidation theory for unsaturated soils. The coefficients of permeability and volume change for unsaturated soils are assumed to remain constant throughout the consolidation process. The mathematical expression of the present solution is much simpler compared with the previous available solutions in the literature. Two new variables are introduced to transform the two coupled governing equations of pore‐water and pore‐air pressures into an equivalent set of partial differential equations, which are easily solved with standard mathematical formulas. It is shown that the present analytical solution can be degenerated into that of Terzaghi consolidation for fully saturated condition. The analytical solutions to 1‐D consolidation of an unsaturated soil subjected to instantaneous loading, ramp loading, and exponential loading, for different drainage conditions and initial pore pressure conditions, are summarized in tables for ease of use by practical engineers. In the case studies, the analytical results show good agreement with the available analytical solution in the literature. The consolidation behaviors of unsaturated soils are investigated. The average degree of consolidation at different loading patterns and drainage conditions is presented. The pore‐water pressure isochrones for two different drainage conditions and three initial pore pressure distributions are presented and discussed. Copyright © 2013 John Wiley & Sons, Ltd.     

Study of desiccation crack initiation and development at ground surface

Cracks are common in clayey or expansive soils and provide preferential pathways for water infiltration into soils. A field study was conducted to investigate the mechanisms of initiation and development of desiccation cracks on two soils. Based on results of the field study and measurements of soil-water retention curves and soil shrinkage properties in the laboratory, the conditions of crack initiation and the development of crack geometric parameters (i.e., crack porosity, crack aperture, and crack density) with water content or drying time were quantified. The results show that desiccation cracks developed in three stages: initial stage, primary stage, and steady state stage. In the initial stage, few cracks developed with gradually decreasing water content. When the water content reached a critical value for crack initiation, cracks developed quickly and this was the start of the primary stage. The critical suction at crack initiation was calculated using a stress criterion, which is in the range of 5.3–21.3 times the preconsolidation pressure of the soil. As the water content approached the shrinkage limit of the soil, cracks developed slowly and approached a steady state. The cracks were found to be repeatable during three drying–wetting cycles.     

Zarka shakedown modelling of expansive soils subjected to wetting and drying cycles

Unsaturated expansive soils subjected to wetting and drying cycles result in huge differential settlements of structures built on these materials. The existed models for these materials present large number of parameters that lead to time-consuming procedure to characterise their mechanical behaviour during wetting–drying cycles. In this context, Zarka shakedown theory previously applied to the mechanical loading of granular materials has been used for expansive soils subjected to suction cycles. The parameters of this shakedown-based model were calibrated for two different expansive soils. The comparisons between the experimental results and the calculations for the different tests, demonstrate the capacity of Zarka shakedown theory to simulate the mechanical behaviour of unsaturated expansive soils.     

Semi-analytical solution to one-dimensional consolidation for viscoelastic unsaturated soils

This paper presents a semi-analytical solution to one-dimensional consolidation of viscoelastic unsaturated soils with a finite thickness under oedometric conditions and subjected to a sudden loading. The solution is obtained by using Lee’s correspondence principle based on the semi-analytical solution to one-dimensional consolidation of elastic unsaturated soils. The boundary contains the top surface permeable to water and air and the bottom impermeable to water and air. A typical example is given to show the evolution of excess pore-air and pore-water pressures as well as the total degree of consolidation of the soil layer with time for different ratios of air–water permeability coefficient, elastic modulus and viscoelastic coefficient. The one-dimensional consolidation behavior of viscoelastic unsaturated soil is discussed according to the semi-analytical solution. These results contribute to a better understanding of the consolidation behavior of viscoelastic unsaturated soils.     

FEM validation of a double porosity elastic model for consolidation of structurally complex clayey soils

Laboratory consolidation of structured clayey soils is analysed in this paper. The research is carried out by two different methods. The first one treats the soil as an isotropic homogeneous equivalent Double Porosity (DP) medium. The second method rests on the extensive application of the Finite Element Method (FEM) to combinations of different soils, composing 2D or fully 3D ordered structured media that schematically discretize the complex material. Two reference problems, representing typical situations of 1D laboratory consolidation of structured soils, are considered. For each problem, solution is obtained through integration of the equations governing the consolidation of the DP medium as well as via FEM applied to the ordered schemes composed of different materials. The presence of conventional experimental devices to ensure the drainage of the sample is taken into account through appropriate boundary conditions. Comparison of FEM results with theoretical results clearly points out the ability of the DP model to represent consolidation processes of structurally complex soils. Limits of applicability of the DP model may arise when the rate of fluid exchange between the two porous systems is represented through oversimplified relations. Results of computations, obtained having assigned reasonable values to the meso‐structural and to the experimental apparatus parameters, point out that a partially efficient drainage apparatus strongly influences the distribution along the sample and the time evolution of the interstitial water pressure acting in both systems of pores. Data of consolidation tests in a Rowe's cell on samples of artificially fissured clays reported in the literature are compared with the analytical and numerical results showing a significant agreement. Copyright © 2000 John Wiley & Sons, Ltd.     

裂缝对膨胀土抗剪强度指标影响的试验研究

基于室内直剪试验,测定了膨胀土经历干湿循环后的抗剪强度指标。试验基本反映了膨胀土在干湿循环过程中裂缝的开展现象,揭示了正是由于裂缝的开展使膨胀土强度降低。结果表明,无论是黏聚力还是内摩擦角随裂缝的开展,都会产生一定程度的衰减,相比之下,裂缝对黏聚力的影响更加突出,对其取值也应更为谨慎。膨胀土的抗剪强度随裂缝的开展而衰减的规律符合双曲线模式,提出反映膨胀土抗剪强度随干湿循环次数增加,实际上是随裂缝的开展,而降低规律的经验公式。工程实际中膨胀土的抗剪强度反映的是包含各种裂隙的土体复合强度,因此在设计中,强度参数的取用应以考虑裂隙作用的土体强度为基本依据。所建议的方法可以使今后研究膨胀土裂隙的发生、发展对抗剪强度的影响变得相对简单,具有可操作性。     

Influencing factors of geometrical structure of surface shrinkage cracks in clayey soils

In order to investigate the effects of temperature, thickness of soil layer, wetting and drying cycles and soil types on geometrical structure of surface shrinkage cracks in clayey soils, special software Crack Image Analysis System (CIAS) for analyzing shrinkage crack patterns was developed. Eight groups of soil samples were prepared and subjected to drying to crack in laboratory. The number of crack segments and intersections, average crack length, width and aggregate area, crack intensity factor (CIF), and the corresponding probability density functions (PDF) of these parameters were determined by analyzing several crack patterns derived from different experimental conditions. The results show that the soil cracking behavior and the geometrical structure of crack patterns are significantly influenced by these considered factors. There is a tendency of crack length, width, aggregate area and their most probable value (MPV) related to the PDF increases with temperature increase. With thicker soil layers, the average crack length, width, aggregate area and CIF are increased, and the main distribution ranges of crack length, width and aggregate area are increased also. When the soil is subjected to multiple wetting–drying cycles, the soil surface generates more irregular and coarse cracks. The number of short and narrow crack segments increases significantly, and the CIF decreases with an increase in wetting–drying cycles. It is also observed that the extent of cracking is directly related to the soil fines fraction and its plasticity index (I_P). The greatest CIF and crack width are observed in the soils with the largest fines fraction and highest I_P. In addition, the ratio of numbers of crack segments to intersections ranges from 1.5 to 2, and cracking mainly takes place in three stages: main-cracks initiation stage; sub-cracks initiation stage; terminal stable stage.     

A bounding surface model for unsaturated soils considering the microscopic pore structure and interparticle bonding effect due to water menisci

The interparticle bonding effect due to water menisci plays an important role in the hydromechanical coupling properties of unsaturated soils. This paper presents an unsaturated hydromechanical coupling model that considers the influence of matric suction, degree of saturation, and microscopic pore structure on the interparticle bonding effect. The enhanced effective stress and bonding variable are selected as constitutive variables. The bonding variable is correlated with the ratio between unsaturated void ratio and saturated void ratio. The deformation characteristics of unsaturated soils are described based on the bounding surface plasticity theory. A soil–water characteristic model that considers deformation and hydraulic hysteresis is integrated into the constitutive model to achieve hydromechanical coupling. The proposed model can effectively describe the hydromechanical coupling characteristics and the meniscus bonding force of unsaturated bimodal structure soils; the model parameters can be easily obtained through routine experiments. The experimental results of unsaturated isotropic compression, the wetting/drying cycle, and unsaturated triaxial shear tests are used to validate the capability of the proposed model.

     

Modelling of shear stiffness of unsaturated fine grained soils at very small strains

The shear modulus at very small strains (less than 0.001%) is an important parameter in the design of geotechnical structures subjected to static and cyclic loadings. Although numerous soil models are available for predicting shear modulus of saturated and dry soils, only a few ones can predict shear stiffness at very small strains of unsaturated soils correctly. In this study, a few unsaturated soil models are evaluated critically and compared with a newly developed model. This newly proposed model is verified by using measured shear modulus at very small strains for three different low plasticity fine grained soils available in the literature. It is found that this new model can predict shear modulus at very small strain resulting from an increase and a decrease in mean net stress at constant matric suction for low plasticity fine grained soils. Moreover, this model is able to give a reasonably good prediction on shear stiffness at very small strain during wetting of a collapsible unsaturated soil. In addition, the newly proposed model is illustrated to capture a consistent trend with experimental data of shear stiffness at very small strain for non-collapsible soils obtained during drying–wetting cycles. This evaluation revealed that the newly proposed model has better predictive capabilities than some earlier formulations of the same simplicity. In addition, the proposed model with fewer parameters has similar predictive capability as compared with a more complex model.     

膨胀土干湿循环胀缩裂隙演化的CT试验研究

对南阳重塑膨胀土在干湿循环过程中裂隙的演化进行了CT试验研究，试样的CT图像显示了原有裂隙开展，新裂隙产生，裂缝数量增加并连通，最后形成裂隙网络的过程。定义了基于CT数据的裂隙损伤变量，分析了裂隙损伤变量随累计干缩体变的变化规律。     

分段循环荷载作用下二维非饱和土固结特性分析

基于Dakshanamurthy和Fredlund提出的二维非饱和土固结理论,利用Fourier正弦级数展开、Laplace变换,分别给出了分段循环荷载作用下二维非饱和土固结问题的超孔隙气压力、超孔隙水压力和沉降的半解析解,并应用退化法验证了本文所得半解析解的正确性。然后,结合3种具体的荷载形式,分析了分段循环荷载作用下气相与液相渗透系数之比(k_a/k_w)、水平方向与竖直方向渗透系数之比(k_x/k_z)和荷载特征参数(a)对二维非饱和土固结特性的影响。结果表明：k_a/k_w和k_x/k_z的增大均会加速固结沉降进程;荷载特征参数越大,沉降发展越早,沉降值越小;二维非饱和土固结特性受分段循环荷载作用影响明显。因此,在实际施工过程中改变施工速度和设置径向排水装置可有效控制二维条件下非饱和土体的固结过程,该研究成果可为非饱和土地基的设计和施工提供重要理论依据。     

干湿交替作用后砂岩破裂过程实时观测与分析

采用岩石破裂全过程的细观力学试验系统进行了天然状态及干湿交替作用后的完整和含预制圆孔砂岩试件的单轴压缩破裂试验。利用显微镜观察了砂岩裂纹萌生、扩展以及试件的失稳破坏过程。经过干湿交替作用后岩样更容易在较低的应力状态下发生裂纹的扩展,并且裂纹扩展的初始阶段一般表现为拉破裂。完整试件的裂纹是随机性的三维扩展,而有预制圆孔试件一般是从预制圆孔周边开始扩展。天然状态下的预制圆孔试件在起裂时从圆孔的压应力集中区产生剪切裂纹;而干湿交替作用后,起裂时可能先从拉应力集中区产生拉裂纹,并且裂纹演化过程更为复杂一些,破裂形式多样化。与天然状态相比,经过干湿循环作用后的砂岩破坏脆性降低,强度等力学参数也随干湿交替次数增加而逐渐降低。     

A hypoplastic model for mechanical response of unsaturated soils

A new constitutive model is developed for the mechanical behaviour of unsaturated soils based on the theory of hypoplasticity and the effective stress principle. The governing constitutive relations are presented and their application is demonstrated using several experimental data from the literature. Attention is given to the stiffening effect of suction on the mechanical response of unsaturated soils and the phenomenon of wetting‐induced collapse. All model parameters have direct physical interpretation, procedures for their quantification from test data are highlighted. Quantitative predictions of the model are presented for wetting, drying and constant suction tests. Copyright © 2008 John Wiley & Sons, Ltd.     

Hydro-mechanical hypoplastic models for unsaturated soils under isotropic stress conditions

The current work presents two hypoplastic models to predict the hydraulic and mechanical behavior of the unsaturated soils. The hydraulic model is capable to reproduce the observed hysteretic behavior for the drying and wetting processes and incorporates the void ratio dependence. The mechanical model is written in terms of the Bishop effective stress and is limited only to isotropic states. It incorporates a normal consolidation line that depends on the degree of saturation. Both models were developed under the hypoplastic framework and therefore are capable to simulate the observed behavior without defining a yield surface. Simulations of some experiments performed with Pearl clay suggest that the models accurately predict its hydro-mechanical behavior.     

A novel approach for modelling soil–water characteristic curves with hysteresis

Soil–water characteristic curves can be defined as the relationship between the degree of saturation and suction of an unsaturated soil. Geomaterials, such as clays, sands, and geotextiles, usually exhibit hysteresis between drying and wetting curves. In addition, each drying and wetting curve is nonlinear in shape, which may be approximated by sigmoid curves. In geotechnical engineering, it is common to adopt analytical expressions for these curves that must be calibrated iteratively by trying different values for the constitutive parameters. In this paper, a novel approach for modelling the nonlinear saturation–suction response with hysteresis is presented, where a simple differential equation is introduced to describe the shapes of the curves. The great advantage of this new technique is the ease with which the parameters can be determined. In addition, the implementation of the resulting equations into fully hydro-mechanical models for numerical analyses is straightforward. Some features of the behaviour predicted with the new representation are studied and validations against real laboratory curves for soils are presented. The technique is simple, yet versatile due to the rational basis used in the deduction of the equations, which allows for future extensions to soils displaying more complex unsaturated behaviour.     

Numerical modeling of multiphase fluid flow in deforming porous media: A comparison between two- and three-phase models for seismic analysis of earth and rockfill dams

In this paper, a fully coupled numerical model is presented for the finite element analysis of the deforming porous medium interacting with the flow of two immiscible compressible wetting and non-wetting pore fluids. The governing equations involving coupled fluid flow and deformation processes in unsaturated soils are derived within the framework of the generalized Biot theory. The displacements of the solid phase, the pressure of the wetting phase and the capillary pressure are taken as the primary unknowns of the present formulation. The other variables are incorporated into the model using the experimentally determined functions that define the relationship between the hydraulic properties of the porous medium, i.e. saturation, relative permeability and capillary pressure. It is worth mentioning that the imposition of various boundary conditions is feasible notwithstanding the choice of the primary variables. The modified Pastor–Zienkiewicz generalized constitutive model is introduced into the mathematical formulation to simulate the mechanical behavior of the unsaturated soil. The accuracy of the proposed mathematical model for analyzing coupled fluid flows in porous media is verified by the resolution of several numerical examples for which previous solutions are known. Finally, the performance of the computational algorithm in modeling of large-scale porous media problems including the large elasto-plastic deformations is demonstrated through the fully coupled analysis of the failure of two earth and rockfill dams. Furthermore, the three-phase model is compared to its simplified one which simulates the unsaturated porous medium as a two-phase one with static air phase. The paper illustrates the shortcomings of the commonly used simplified approach in the context of seismic analysis of two earth and rockfill dams. It is shown that accounting the pore air as an independent phase significantly influences the unsaturated soil behavior.