Analysis of pressuremeter geometry effects in clay using critical state models

The conventional interpretation methods of pressuremeter testing effectively approximate pressuremeter membranes as infinitely long. As a result, the effects of the two‐dimensional geometry of pressuremeters are ignored, leading to an overestimation of soil shear strength by pressuremeter testing, as demonstrated in several previous studies. This paper presents results of a numerical study of two‐dimensional geometry effects on self‐boring pressuremeter tests in undrained clay. The results are obtained using critical state soil models with an effective stress formulation. This is in contrast to most (if not all) existing studies on pressuremeter geometry effects, which were based on perfectly plastic soil models (e.g. Yu (Cavity expansion theory and its application to the analysis of pressuremeters. DPhil Thesis, The University of Oxford, 1990), Yeung and Carter (Proc. 3rd Int. Symp. on Pressuremeters, 1990), and Houlsby and Carter (Géotechnique, 1993; 43 (4):567–576)). The present study suggests that the overestimation of soil strength due to the neglect of finite pressuremeter length is significantly affected by the soil model used in the calculations. It is found that for clays with a high overconsolidation ratio (OCR) the strength overestimation predicted using critical state soil models could be considerably smaller than that predicted using perfectly plastic soil models. The main conclusion of this numerical study is that care must be exercised before directly applying any numerically determined pressuremeter geometry correction factors in practice. Copyright © 2005 John Wiley & Sons, Ltd.     

Genealogy of hypoplasticity and barodesy

This article is an attempt at providing an insight into the development of hypoplasticity (including barodesy, which is a recent development of hypoplasticity) as a theory elaborated since 1977, when the first version was published by the first author, until present. The multiplicity of the many versions published since then is hard to overlook. This article presents a review and insight into the evolution of a theory and the struggle to formulate a satisfactory constitutive law. Among the many proposed versions, we focus on those ones that can be seen as changes of paradigm. Copyright © 2016 John Wiley & Sons, Ltd.     

Analytical modeling on consolidation of stiffened deep mixed column-reinforced soft soil under embankment

Stiffened deep mixed (SDM) column is a new ground improvement technique to improve soft soil, which can be used to increase bearing capacity, reduce deformation, and enhance stability of soft soil. This technique has been successfully adopted to support the highway and railway embankments over soft soils in China and other countries. However, there have been limited investigations on its consolidation under embankment loading. This paper developed an analytical solution for the consolidation of embankment over soft soil with SDM column in which core pile is equal to or shorter than outer DM column. The consolidation problem was simplified as a consolidation of composite soil considering the load shear effect of core pile. The developed solution was verified by a comparison with the results computed by three-dimensional (3-D) finite element analysis. A parametric study based on the derived solution was conducted to investigate influence factors—length of core pile, diameter of core pile, diameter of SDM column, modulus of DM column, and permeability coefficient of DM column—on the consolidation behavior of SDM column-supported embankment over soft soil. The developed solution was applied to a case history of SDM column-supported embankment, and a good agreement was found between the predictions and the field measurements.     

淄博市西宝山耐火黏土矿地质特征及找矿方向

西宝山耐火黏土矿为沉积型矿床,目前揭露1个耐火黏土矿矿层,赋存于石盒子群万山组底部“B”层黏土岩的下部,厚度较稳定,勘查开发潜力较大。通过对其地质特征、矿层特征、矿石结构构造及矿床成因分析,提出了找矿方向。认为该耐火黏土矿受沉积地层控制,呈层状分布,分布较稳定。依据该区地层产出情况,在青龙山断裂与禹王山断裂之间及两侧均具有良好的找矿前景,这为下一步找矿工作奠定了基础。     

Integration of aquifer geology,groundwater flow and arsenic distribution in deltaic aquifers – A unifying concept

Groundwater arsenic (As) presents a public health risk of great magnitude in densely populated Asian delta regions, most acutely in the Bengal Basin (West Bengal, India and Bangladesh). Research has focused on the sources, mobilisation, and heterogeneity of groundwater As, but a consistent explanation of As distribution from local to basin scale remains elusive. We show for the Bengal Aquifer System that the numerous, discontinuous silt‐clay layers together with surface topography impose a hierarchical pattern of groundwater flow, which constrains As penetration into the aquifer and controls its redistribution towards discharge zones, where it is re‐sequestered to solid phases. This is particularly so for the discrete periods of As release to groundwater in the shallow subsurface associated with sea level high‐stand conditions of Quaternary inter‐glacial periods. We propose a hypothesis concerning groundwater flow ( S ilt‐clay layers I mpose H ierarchical groundwater flow patterns constraining A rsenic progression [SIHA]), which links consensus views on the As source and history of sedimentation in the basin to the variety of spatial and depth distributions of groundwater As reported in the literature. SIHA reconciles apparent inconsistencies between independent, in some cases contrasting, field observations. We infer that lithological and topographic controls on groundwater flow, inherent to SIHA, apply more generally to deltaic aquifers elsewhere. The analysis suggests that groundwater As may persist in the aquifers of Asian deltas over thousands of years, but in certain regions, particularly at deeper levels, As will not exceed low background concentrations unless groundwater flow systems are short‐circuited by excessive pumping.     

Dynamic analysis of strain localization in water‐saturated clay using a cyclic elasto‐viscoplastic model

A computational framework is presented for dynamic strain localization and deformation analyses of water‐saturated clay by using a cyclic elasto‐viscoplastic constitutive model. In the model, the nonlinear kinematic hardening rule and softening due to the structural degradation of soil particles are considered. In order to appropriately simulate the large deformation phenomenon in strain localization analysis, the dynamic finite element formulation for a two‐phase mixture is derived in the updated Lagrangian framework. The shear band development is shown through the distributions of viscoplastic shear strain, the axial strain, the mean effective stress, and the pore water pressure in a normally consolidated clay specimen. From the local stress–strain relations, more brittleness is found inside the shear bands than outside of them. The effects of partially drained conditions and mesh‐size dependency on the shear banding are also investigated. The effect of a partially drained boundary is found to be insignificant on the dynamic shear band propagation because of the rapid rate of applied loading and low permeability of the clay. Using the finer mesh results in slightly narrower shear bands; nonetheless, the results manifest convergency through the mesh refinement in terms of the overall shape of shear banding and stress–strain relations. Copyright © 2013 John Wiley & Sons, Ltd.     

A bounding surface model for anisotropically overconsolidated clay incorporating thermodynamics admissible rotational hardening rule

A bounding surface model is formulated to simulate the behavior of clays that are subject to an anisotropic consolidation stress history. Conventional rotational hardening is revisited from the perspective of thermodynamics. As the free energy cannot be accumulated infinitely upon critical state failure, the deviatoric back stress must vanish. This requires the rotated yield surface to be turned back to eventually align on the hydrostatic axis in the stress plane. Noting that most of the previous propositions violate this restriction, an innovative rotational hardening rule is formulated that is thermodynamically admissible. The bounding surface framework that employs the modified yield surface is applied to simulate elastoplastic deformations for overconsolidated clays, with which the overprediction of strength on the “dry” side can be greatly improved with reasonable results. Other important features, including contractive or dilative response and hardening or softening behavior, can also be well-captured. It has been shown that the model can simulate three types of reconstituted clays that are sheared with initial conditions over a wide range of anisotropic consolidation stress ratios and overconsolidation ratios under both triaxial undrained and drained conditions. Limitations and potential improvement of the model regarding the fabric anisotropy at critical state have been discussed.     

Bound and mobile soil water isotope ratios are affected by soil texture and mineralogy,whereas extraction method influences their measurement

Questions persist about interpreting isotope ratios of bound and mobile soil water pools, particularly relative to clay content and extraction conditions. Interactions between pools and resulting extracted water isotope composition are presumably related to soil texture, yet few studies have manipulated the bound pool to understand its influence on soil water processes. Using a series of drying and spiking experiments, we effectively labelled bound and mobile water pools in soils with varying clay content. Soils were first vacuum dried to remove residual water, which was then replaced with heavy isotope-enriched water prior to oven drying and spiking with heavy isotope-depleted water. Water was extracted via centrifugation or cryogenic vacuum distillation (at four temperatures) and analysed for oxygen and hydrogen isotope ratios via isotope ratio mass spectrometry. Water from centrifuged samples fell along a mixing line between the two added waters but was more enriched in heavy isotopes than the depleted label, demonstrating that despite oven drying, a residual pool remains and mixes with the mobile water. Soils with higher clay + silt content appeared to have a larger bound pool. Water from vacuum distillation samples have a significant temperature effect, with high temperature extractions yielding progressively more heavy isotope-enriched values, suggesting that Rayleigh fractionation occurred at low temperatures in the vacuum line. By distinctly labelling bound and mobile soil water pools, we detected interactions between the two that were dependent on soil texture. Although neither extraction method appeared to completely extract the combined bound and mobile (total water) pool, centrifugation and high temperature cryogenic vacuum distillations were comparable for both δ²H and δ¹⁸O of soil water isotope ratios.     

Micromechanical approach to swelling behavior of capillary-porous media with coupled physics

A detailed multiscale analysis is presented of the swelling phenomenon in unsaturated clay-rich materials in the linear regime through homogenization. Herein, the structural complexity of the material is formulated as a three-scale, triple porosity medium within which microstructural information is transmitted across the various scales, leading ultimately to an enriched stress-deformation relation at the macroscopic scale. As a side note, such derived relationship leads to a tensorial stress partitioning that is reminiscent of a Terzaghi-like effective stress measure. Otherwise, a major result that stands out from previous works is the explicit expression of swelling stress and capillary stress in terms of micromechanical interactions at the very fine scale down to the clay platelet level, along with capillary stress emerging due to interactions between fluid phases at the different scales, including surface tension, pore size, and morphology. More importantly, the swelling stress is correlated with the disjoining forces due to electrochemical effects of charged ions on clay minerals and van der Waals forces at the nanoscale. The resulting analytical expressions also elucidate the role of the various physics in the deformational behavior of clayey material. Finally, the capability of the proposed formulation in capturing salient behaviors of unsaturated expansive clays is illustrated through some numerical examples.     

Diffusive mass exchange of non‐reactive substances in dual‐porosity porous systems – column experiments under saturated conditions

Diffusive mass exchange into immobile water regions within heterogeneous porous aquifers influences the fate of solutes. The percentage of immobile water is often unidentified in natural aquifers though. Hence, the mathematical prediction of solute transport in such heterogeneous aquifers remains challenging. The objective of this study was to find a simple analytical model approach that allows quantifying properties of mobile and immobile water regions and the portion of immobile water in a porous system. Therefore, the Single Fissure Dispersion Model (SFDM), which takes into account diffusive mass exchange between mobile and immobile water zones, was applied to model transport in well‐defined saturated dual‐porosity column experiments. Direct and indirect model validation was performed by running experiments at different flow velocities and using conservative tracer with different molecular diffusion coefficients. In another column setup, immobile water regions were randomly distributed to test the model applicability and to determine the portion of immobile water. In all setups, the tracer concentration curves showed differences in normalized maximum peak concentration, tailing and mass recovery according to their diffusion coefficients. These findings were more pronounced at lower flow rates (larger flow times) indicating the dependency of diffusive mass exchange into immobile water regions on tracers' molecular diffusion coefficients. The SFDM simulated all data with high model efficiency. Successful model validation supported the physical meaning of fitted model parameters. This study showed that the SFDM, developed for fissured aquifers, is applicable in porous media and can be used to determine porosity and volume of regions with immobile water. Copyright © 2015 John Wiley & Sons, Ltd.