Explicit integration of a porosity‐dependent hydro‐mechanical model for unsaturated soils

An adaptive substepping explicit integration scheme is developed for a porosity‐dependent hydro‐mechanical model for unsaturated soils. The model is referred to as the modified σ –Θ model in this paper, which features the employment of the subloading surface plasticity and the stress–saturation approach. On numerical aspects, convex/nonconvex subloading surfaces in the σ –Θ space may result in incorrect loading–unloading decisions during the integration. A new loading–unloading decision method is developed here to solve the problem and then embedded into the explicit integration scheme for the modified σ –Θ model. In addition, to enhance the accuracy of the explicit integration, local errors from both hydraulic and mechanical components are included in the error control for each substep. A drift correction method is also developed to ensure the state point lies on the subloading surface in the σ –Θ space within a set error level. The performance of the loading–unloading decision method for the modified σ –Θ model is discussed through comparing it with the conventional loading–unloading decision method. The importance of involving the hydraulic component in the error control is also demonstrated. The accuracy and efficiency of the proposed adaptive substepping explicit integration scheme for the modified p–Θ model are also studied via several numerical examples. Copyright © 2016 John Wiley & Sons, Ltd.     

Consolidation in spatially random unsaturated soils based on coupled flow‐deformation simulation

This paper integrates random field simulation of soil spatial variability with numerical modeling of coupled flow and deformation to investigate consolidation in spatially random unsaturated soil. The spatial variability of soil properties is simulated using the covariance matrix decomposition method. The random soil properties are imported into an interactive multiphysics software COMSOL to solve the governing partial differential equations. The effects of the spatial variability of Young's modulus and saturated permeability together with unsaturated hydraulic parameters on the dissipation of excess pore water pressure and settlement are investigated using an example of consolidation in a saturated‐unsaturated soil column because of loading. It is found that the surface settlement and the pore water pressure profile during the process of consolidation are significantly affected by the spatially varying Young's modulus. The mean value of the settlement of the spatially random soil is more than 100% greater than that of the deterministic case, and the surface settlement is subject to large uncertainty, which implies that consolidation settlement is difficult to predict accurately based on the conventional deterministic approach. The uncertainty of the settlement increases with the scale of fluctuation because of the averaging effect of spatial variability. The effects of spatial variability of saturated permeability k_sat and air entry parameters are much less significant than that of elastic modulus. The spatial variability of air entry value parameters affects the uncertainties of settlement and excess pore pressure mostly in the unsaturated zone. Copyright © 2016 John Wiley & Sons, Ltd.     

Assessment of slope stability in cohesive soils due to a rainfall

The primary focus in this work is on proposing a methodology for the assessment of stability of natural/engineered slopes in clayey soils subjected to water infiltration. In natural deposits of fine‐grained soils, the presence of water in the vicinity of minerals results in an interparticle bonding. This effect cannot be easily quantified as it involves complex chemical interactions at the micromechanical level. Here, the evolution of strength properties, including the apparent cohesion resulting from initial suction at the irreducible fluid saturation, is described by employing the framework of chemoplasticity. The paper provides first the formulation of the problem; this involves specification of the constitutive relation, development of an implicit return mapping scheme, and the outline of a coupled transient formulation. The framework is then applied to examine the stability of a slope subjected to a prolonged period of intensive rainfall. It is shown that the water infiltration may trigger the loss of stability resulting from the degradation of material properties. Copyright © 2013 John Wiley & Sons, Ltd.     

Influence of pore size and geometry on peat unsaturated hydraulic conductivity computed from 3D computed tomography image analysis

In organic soils, hydraulic conductivity is related to the degree of decomposition and soil compression, which reduce the effective pore diameter and consequently restrict water flow. This study investigates how the size distribution and geometry of air‐filled pores control the unsaturated hydraulic conductivity of peat soils using high‐resolution (45 µm) three‐dimensional (3D) X‐ray computed tomography (CT) and digital image processing of four peat sub‐samples from varying depths under a constant soil water pressure head. Pore structure and configuration in peat were found to be irregular, with volume and cross‐sectional area showing fractal behaviour that suggests pores having smaller values of the fractal dimension in deeper, more decomposed peat, have higher tortuosity and lower connectivity, which influences hydraulic conductivity. The image analysis showed that the large reduction of unsaturated hydraulic conductivity with depth is essentially controlled by air‐filled pore hydraulic radius, tortuosity, air‐filled pore density and the fractal dimension due to degree of decomposition and compression of the organic matter. The comparisons between unsaturated hydraulic conductivity computed from the air‐filled pore size and geometric distribution showed satisfactory agreement with direct measurements using the permeameter method. This understanding is important in characterizing peat properties and its heterogeneity for monitoring the progress of complex flow processes at the field scale in peatlands. Copyright © 2010 John Wiley & Sons, Ltd.     

Dynamics of unsaturated poroelastic solids at finite strain

We derive the governing equations for the dynamic response of unsaturated poroelastic solids at finite strain. We obtain simplified governing equations from the complete coupled formulation by neglecting the material time derivative of the relative velocities and the advection terms of the pore fluids relative to the solid skeleton, leading to a so‐called u^s ? p^w ? p^a formulation. We impose the weak forms of the momentum and mass balance equations at the current configuration and implement the framework numerically using a mixed finite element formulation. We verify the proposed method through comparison with analytical solutions and experiments of quasi‐static processes. We use a neo‐Hookean hyperelastic constitutive model for the solid matrix and demonstrate, through numerical examples, the impact of large deformation on the dynamic response of unsaturated poroelastic solids under a variety of loading conditions. Copyright © 2011 John Wiley & Sons, Ltd.     

Rate and historical change of direct recharge from precipitation constrained by unsaturated zone profiles of chloride and oxygen‐18 in dry river bed of North China Plain

Quantifying of direct recharge derived from precipitation is crucial for assessing sustainability of well‐irrigated agriculture. In the North China Plain, the land use is dominated by groundwater‐irrigated farmland where the direct recharge derived from precipitation and irrigation. To characterize the mean rate and historical variance of direct recharge derived from precipitation, unsaturated zone profiles of chloride and δ¹⁸O in the dry river bed of the Beiyishui River were employed. The results show that archival time scale of the profile covers the duration from 1980 to 2002 (corresponding to depths from 5 to 2 m) which is indicated by matching the δ¹⁸O peaks in the isotope profile with the aridity indexes gained by instrumental records of annual precipitation and annual potential evaporation. Using the chloride mass balance method, the mean rate of the direct recharge corresponding to the archival time scale is estimated to be 3·8 ± 0·8 mm year^?1, which accounts for about 0·7% of the long‐term average annual precipitation. Further, the direct recharge rates vary from 2·1 to 6·8 mm year^?1 since 1980. Despite the subhumid climate, the estimate of recharge rates is in line with other findings in semiarid regions. The low rate of direct recharge is considered as a result of the relative dry climate in recent decades. In dry river bed, unsaturated zone profiles of chloride and δ¹⁸O combined with instrumental records could offer valuable information about the direct recharge derived from precipitation during droughts. Copyright © 2011 John Wiley & Sons, Ltd.     

用高精度小型CT测定骨密度

骨质疏松症是一种常见的代谢性骨病，其诊断的主要依据为骨密度降低。目前的测量骨密度方法在精度、剂量方面不很理想．本文介绍一种外周型高精度小型ＣＴ，精度高、剂量低是一种较理想的测定骨密度的仪器。文中介绍了提高精度和降低剂量的措施，以及仪器在临床应用试验的情况和结果，此仪器并可扩展应用于骨结构分析。     

裂隙概化模型的非饱和渗流试验研究

借鉴多孔介质非饱和渗流试验的研究成果,研制出一套可同时测定单裂隙毛细压力-饱和度以及非饱和渗透系数-毛细压力关系的试验装置.为检验试验装置的可信度和试验原理的正确性,并初步探讨单裂隙非饱和渗流的机理,特制作了一阶梯开度的裂隙概化模型(本文称之为“S-H”裂隙模型),并在上述试验装置上进行了“S-H”裂隙模型的非饱和渗流试验,初步阐明了单裂隙非饱和渗流的一些基本水力特性,同时试验结果也表明本文试验装置和试验原理合理可靠.     

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

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