超固结非饱和土本构关系的三维化方法

本构模型的三维化对反映材料在三维应力状态下的力学特性具有重要意义。在π平面上,广义Mises三维化方法不能反映应力洛德角对材料变形和强度的影响,g(θ)等方法不能合理地反映应力水平对变形和强度的影响。变换应力三维化方法考虑了应力洛德角影响,数学表达简洁,不增加模型参数,便于计算应用。将作者等人所提的超固结非饱和土本构模型与基于SMP（spatially mobilized plane）的变换应力方法结合,实现了超固结非饱和土的三维化,达到了从剪切屈服到剪切破坏的连续过渡。该模型能够描述不同应力洛德角下超固结非饱和土的硬化、软化、剪胀特性,并与试验数据进行了对比     

强降雨条件下基岩型层状边坡入渗模型及稳定性研究

强降雨作用下基岩型层状边坡易发生失稳破坏,给人们生命财产造成重大损失。为探索基岩型层状边坡的降雨入渗过程,以Green-Ampt入渗模型为基础,考虑了边坡几何特征及饱和带渗流作用,得到了边坡不同入渗阶段入渗率和湿润峰深度的计算式,建立了适用于基岩型层状边坡的降雨入渗计算模型;在分析基岩型层状边坡稳定性时,考虑了饱和带的渗流作用,联合降雨入渗模型与极限平衡法,分析了基岩面与湿润峰面稳定性的变化规律,得到了计算基岩型层状边坡安全系数的解析表达式;结果表明,该模型的计算与试验结果具有一致性。利用该计算模型与传统计算方法对张家湾滑坡进行了强度为30 mm/h的降雨入渗与稳定性分析,得到了不同计算方法下张家湾滑坡湿润峰深度及安全系数随降雨历时增加而发生的变化规律,结果表明该计算模型在分析基岩型层状边坡稳定性时优于传统分析方法。     

A density‐dependent elastoplastic hydro‐mechanical model for unsaturated compacted soils

This paper presents a three‐dimensional elastoplastic constitutive model for predicting the hydraulic and mechanical behaviour of unsaturated soils. It is based on experimental results obtained from a series of controlled‐suction triaxial tests on unsaturated compacted clay with different initial densities. Hydraulic hysteresis in the water‐retention behaviour is modelled as an elastoplastic process, with the elastic part modelled by a series of scanning curves and the elastoplastic part modelled by the main drying and wetting curves. The effect of void ratio on the water‐retention behaviour is studied using data obtained from controlled‐suction wetting–drying cyclic tests on unsaturated compacted clay with different initial densities. The effect of the degree of saturation on the stress–strain‐strength behaviour and the effect of void ratio on the water‐retention behaviour are considered in the model, as is the effect of suction on the hydraulic and mechanical behaviour. The initial density dependency of the compacted soil behaviour is modelled by experimental relationships between the initial density and the corresponding yield stress and, thereby, between the initial density and the normal compression line. The model is generalized to three‐dimensional stress states by assuming that the shapes of the failure and yield surfaces in the deviatoric stress plane are given by the Matsuoka–Nakai criterion. Model predictions of the stress–strain and water‐retention behaviour are compared with those obtained from triaxial tests with different initial densities under isotropic compression, triaxial compression and triaxial extension, with or without variation in suction. The comparisons indicate that the model accurately predicts the hydraulic and mechanical behaviour of unsaturated compacted soils with different initial densities using the same material constant. Copyright © 2006 John Wiley & Sons, Ltd.     

Pool-fills: a window to palaeoflood history and response in bedrock-confined rivers

Channel‐scale sedimentary units associated with bedrock‐controlled riffle‐pool morphology are examined in detail along Sandy Creek gorge, an ephemeral stream in arid south‐eastern central Australia. Pool‐fills comprise cut‐and‐fill assemblages of poorly sorted sediments ranging in texture from muds to boulders. Five unit types are defined based on particle size, sedimentary structures, geometry and bounding surface character: (1) coarse‐grained bar platform; (2) fine‐grained bar supraplatform; (3) fine‐grained pool‐fill; (4) fine‐grained bench; and (5) modern pool‐fill. The last coarse‐grained unit currently lining the pools suggests an altered sedimentation style over the post‐settlement period (post‐ad 1860s). Situated at bedrock valley constrictions, pool‐fills are compared with other sedimentary units associated with recirculating currents: eddy bars and slackwater deposits. But only the fine‐grained bench units reflect eddy recirculation; the pool‐fills are principally forced‐bars associated with bedrock‐controlled or ‘forced’ riffle‐pool morphology. A late Holocene palaeoflood history is proposed based on radiocarbon ages from the pool‐fills: multiple phases of cut‐and‐fill activity were preceded by a superflood 3400–1900 years ago that eroded the pool‐fills to bedrock. The resilience of the pool‐fills was illustrated by the passage of a 1‐in‐100‐year flood in 1992, which caused only minor erosion. The presence of pool‐fills may provide a window to past phases of river activity that cannot be extracted from either historical records/observations or palaeoflood slackwater sediment analyses. The formation and sedimentary preservation potential of these landforms reflect a combination of hydraulic and structural influences, but the occurrence of high‐magnitude floods exerts the dominant control.     

A kinematic hardening based model for unsaturated soils considering different hydraulic conditions

At present, several of the existing elastoplastic constitutive models are adapted for describing the stress–strain behavior of unsaturated soils. However, most of them present certain limitations in this field. These limitations can be related to the basic model and/or added unsaturated state variables and formulations. In this regard, inability to model the hydro‐mechanical behavior in constant water (CW) conditions is an example of these limitations. In this paper, an advanced version of CJS model is selected for adaptation to the unsaturated states. Adaptation to unsaturated states is achieved in the framework of effective stress approach. Effective stress equation and unsaturated state variables are selected based on the recent research existing in the literature. The developed model is capable of describing the complex behavior of unsaturated soil in the CW condition in addition to predicting the behavior at failure and post–failure, nonlinear elastoplastic behavior at low levels of stress and strain (by selecting a very small elastic domain), as well as wetting and collapse behaviors. In order to validate the model, results of triaxial tests in CD and CW conditions are used. The validation results indicate the good capability of the proposed model. Behavior of the unsaturated soils during wetting is an important issue. For this reason, the model is also evaluated based on the results of wetting and collapse triaxial tests. A comparison between the tests and simulation results shows that the model is able to predict the soil behavior under the wetting path. Copyright © 2016 John Wiley & Sons, Ltd.     

Using modified state surface approach to select parameter values in the Barcelona basic model

The Barcelona basic model (BBM) successfully explained many key features of unsaturated soils and received extensive acceptance. It is also one of the few elastoplastic constitutive models for unsaturated soils that have been implemented within finite element codes and applied to the analysis of real boundary value problems. The BBM was proposed in incremental forms according to theories of soil plasticity in which individual aspects of the isotropic virgin behavior are controlled by multiple parameters, whereas at the same time, a single parameter controls more than one aspect of soil behavior. Although a variety of methods have been recently developed for calibrating model parameters for elastoplastic soil models, at present, there are no well‐established, simple, and objective methods for selecting parameter values in the BBM from laboratory tests. This has been one of the major obstacles to the dissemination of this constitutive model beyond the research context. This article presents an optimization approach especially developed for simple and objective identification of material parameters in the BBM. This is achieved by combining a modified state surface approach, recently proposed to model the elastoplastic behavior of unsaturated soils under isotropic stress conditions, with the Newton or quasi‐Newton method to simultaneously determine the five parameters governing isotropic virgin behavior in the BBM. The comparison between results using the proposed method and an existing method for the same laboratory tests was discussed from which the simplicity and objectivity of the proposed method were evaluated. Copyright © 2012 John Wiley & Sons, Ltd.     

Diffuse instabilities with transition to localization in loose granular materials

This paper is concerned with diffuse and other ensuing failure modes in geomaterials when tested under homogeneous states of shearing in various loading programs and drainage conditions. Material instability is indeed the basic property that accounts for the instability of an initially homogeneous deformation field leading to diffuse failure and strain localization in geomaterials. The former is normally characterized by a runaway type of failure accompanied with a sudden and violent collapse of the material in the absence of any localization phenomena. Against this backdrop, we present a brief overview of material instability in elastoplastic solids where one finds a rich source of theoretical concepts including bifurcation, strain localization, diffuse failure and second‐order work, as well as a considerable body of experiments. Some compelling laboratory experimental studies of material instability with focus to diffuse failure are then presented and interpreted based on the second‐order work. Finally, various material instability analyses using an elastoplastic constitutive and a general finite element analysis of the above‐mentioned laboratory experimental tests are presented as a boundary value problem. It is shown that instability can be captured from otherwise uniform stress, density and hydraulic states, whereas uniform deviatoric loads are being applied on the external boundaries of a specimen. Although the numerical simulations reproduce well the laboratory experimental results, they also highlight the hierarchy of failure modes where localization phenomena emerge in the post‐bifurcation regime as a result of a degradation of homogeneity starting from a diffuse mode signalled by a zero second‐order work. Copyright © 2012 John Wiley & Sons, Ltd.     

Triaxial behaviour of transversely isotropic materials: application to sedimentary rocks

Failure and long‐term behaviour of oriented solids are studied. Transversely isotropic materials are considered and a mathematical formulation that respect the material symmetry is developed and applied to model the triaxial behaviour of sedimentary rocks. Two failure criteria and a viscoplastic constitutive model that describe, respectively, triaxial failure and triaxial creep tests are presented and discussed. The application of the developed models to describe the mechanical behaviour of Tournemire shale shows that theoretical predictions are in good agreement with the experimental data. In the present paper, the developed approach is applied to sedimentary rock materials, nevertheless, it can be generalized to any material that exhibits transverse isotropy. Copyright © 2007 John Wiley & Sons, Ltd.     

A novel computational approach for large deformation and post‐failure analyses of segmental retaining wall systems

Segmental retaining wall (SRW) systems are commonly used in geotechnical practice to stabilize cut and fill slopes. Because of their flexibility, these systems can tolerate minor movements and settlements without incurring damage or crack. Despite these advantages, very few numerical studies of large deformations and post‐failure behavior of SRW systems are found in the current literature. Traditional numerical methods, such as the finite element method, suffer from mesh entanglement, thus are unable to simulate large deformations and flexible behavior of retaining wall blocks in SRW systems. To overcome the above limitations, a novel computational framework based on the smoothed particle hydrodynamics (SPH) method was developed to simulate large deformations and post‐failure behavior of soils and retaining wall blocks in SRW systems. The proposed numerical framework is a hybrid continuum/discontinuum approach that can model soil as an elasto‐plastic material and retaining wall blocks as independent rigid bodies associated with both translational and rotational degrees of freedom. A new contact model is proposed within the SPH framework to simulate the interaction between the soil and the blocks and between the blocks. As an application of the proposed numerical method, a two‐dimensional simulation of an SRW collapse was simulated and compared to experimental results conducted under the same conditions. The results showed that the proposed computational approach provided satisfactory agreement with the experiment. This suggests that the new framework is a promising numerical approach to model SRW systems. Copyright © 2014 John Wiley & Sons, Ltd.     

基岩冲刷的数值模拟

为了探索采用数值模拟的方法来研究基岩冲刷的可行性,对二维散粒体河床的冲刷进行了试验研究,并采用二维k-ε模型,在贴体非正交曲线坐标系中用控制体积法对散粒体河床的冲刷进行了数值模拟,计算的冲坑深度与试验结果符合很好,结果令人满意.选用相同的控制基岩冲刷的参数,对不同下游水垫深度、不同入射流速时基岩冲刷深度的正确模拟说明,在对基岩特性有了充分的了解后,用数值方法来模拟基岩的冲刷过程及估算冲刷深度是合理、可行的.     

Analytical solution for one‐dimensional consolidation of unsaturated soils using eigenfunction expansion method

This paper introduces an exact analytical solution for governing flow equations for one‐dimensional consolidation in unsaturated soil stratum using the techniques of eigenfunction expansion and Laplace transformation. The homogeneous boundary conditions adopted in this study are as follows: (i) a one‐way drainage system of homogenous soils, in which the top surface is considered as permeable to air and water, whereas the base is an impervious bedrock; and (ii) a two‐way drainage system where both soil ends allow free dissipation of pore‐air and pore‐water pressures. In addition, the analytical development adopts initial conditions capturing both uniform and linear distributions of the initial excess pore pressures within the soil stratum. Eigenfunctions and eigenvalues are parts of the general solution and can be obtained based on the proposed boundary conditions. Besides, the Laplace transform method is adopted to solve the first‐order differential equations. Once equations with transformed domain are all obtained, the final solutions, which are proposed to be functions of time and depth, can be achieved by taking an inverse Laplace transform. To verify the proposed solution, two worked examples are provided to present the consolidation characteristics of unsaturated soils based on the proposed method. The validation of the recent results against other existing analytical solutions is graphically demonstrated. Copyright © 2013 John Wiley & Sons, Ltd.     

Limitations of suction‐controlled triaxial tests in the characterization of unsaturated soils

Behavior of unsaturated soils is influenced by many factors, and the influences of these factors are usually coupled together. Suction‐controlled triaxial (SCTX) tests are considered to allow researchers to investigate influences of individual variables on unsaturated soils under specified stress path with controls of stresses, pore water, and air pressures. In the past 50 years, SCTX testing method has been established as a standard approach to characterize constitutive behavior of unsaturated soils. Most important concepts for modern unsaturated soil mechanics were developed upon results from the SCTX tests. Among these, one of the most important contributions in the constitutive modeling of elasto‐plastic behavior for unsaturated soils is the Barcelona basic model (BBM) proposed by Alonso et al. in 1990. The BBM successfully explained many features of unsaturated soils and received extensive acceptance. However, the SCTX tests are designed based upon the divide‐and‐conquer approach in which an implicit assumption is used: soil behavior is stress‐path independent. However, it is well‐established that unsaturated soil behavior is elasto‐plastic and stress‐path dependent. It is found that the SCTX tests in fact cannot control the stress path of an unsaturated soil during loading. This incapability, in combination with complicated loading/collapse behavior of unsaturated soils, makes the SCTX tests for characterizing unsaturated soil questionable. This paper discusses the limitations of the SCTX tests in the characterization of unsaturated soils. A possible solution to the problem was proposed based on a newly developed modified state surface approach. The discussions are limited for isotropic conditions. Copyright © 2015 John Wiley & Sons, Ltd.     

The lowermost Mississippi River: a mixed bedrock‐alluvial channel

In this study, the distribution of channel‐bed sediment facies in the lowermost Mississippi River is analysed using multibeam data, complemented by sidescan sonar and compressed high‐intensity radar pulse seismic data, as well as grab and core samples of bed material. The channel bed is composed of a discontinuous layer of alluvial sediment and a relict substratum that is exposed on the channel bed and sidewalls. The consolidated substratum is made up of latest Pleistocene and Early Holocene fluvio‐deltaic deposits and is preferentially exposed in the deepest thalweg segments and on channel sidewalls in river bends. The exposed substratum commonly displays a suite of erosional features, including flutes that are quantitatively similar in form to those produced under known laboratory conditions. A total of five bed facies are mapped, three of which include modern alluvial deposits and two facies that are associated with the relict substratum. A radius of curvature analysis applied to the Mississippi River centreline demonstrates that the reach‐scale distribution of channel‐bed facies is related to river planform. From a broader perspective, the distribution of channel‐bed facies is related to channel sinuosity — higher sinuosity promotes greater substratum exposure at the expense of alluvial sediment. For example, the ratio of alluvial cover to substratum is ca 1·5:1 for a 45 km segment of the river that has a sinuosity of 1·76 and this ratio increases to ca 3:1 for a 120 km segment of the river that has a sinuosity of 1·21. The exposed substratum is interpreted as bedrock and, given the relative coverage of alluvial sediment in the channel, the lowermost Mississippi River can be classified as a mixed bedrock‐alluvial channel. The analyses demonstrate that a mixed bedrock‐alluvial channel boundary can be associated with low‐gradient and sand‐bed rivers near their marine outlet.     

浸矿开采对离子型稀土基岩力学特性的影响研究

为研究原地浸矿母液对离子型稀土矿收液底板的侵蚀影响,以离子型稀土矿基岩（半/弱风化花岗岩）为研究对象,开展了室内浸矿模拟试验。对比分析了不同浸矿时间岩样的应力−应变曲线、峰值应力、峰值应变和弹性模量,获取了不同浸矿历时岩样的孔隙度和细观形貌特征,并运用盒维法对岩样的非线性裂纹特征进行了定量化描述。结果表明：浸矿侵蚀作用对离子型稀土矿底板基岩力学性质的影响主要表现在前150 d,破坏形式由贯穿破坏转变为贯穿破坏与片状剥落并存;浸矿母液侵蚀作用下,基岩力学特性呈现先大幅减小后大幅回升,其后在小范围波动现象,浸矿约60 d对基岩力学特性的影响程度最大;浸矿母液侵蚀造成基岩内小孔隙增长,导致表面出现腐蚀及孔洞现象;浸矿时间不同将造成岩石不同的微观结构变化。     

Three‐dimensional elasto‐plastic model for unsaturated compacted soils with different initial densities

This paper presents an elasto‐plastic model for unsaturated compacted soils and experimental results obtained from a series of suction‐controlled triaxial tests on unsaturated compacted clay with different initial densities. The initial density dependency of the compacted soil behaviour is modelled by establishing experimental relationships between the initial density and the corresponding yield stress and thereby between the initial density and the location and slope of normal compression line. The model is generalized to three‐dimensional stress states by assuming that the shapes of the failure surface and the yield surface in the deviatoric plane are given by the extended SMP criterion. A considerable number of the isotropic compression, triaxial compression and extension tests on unsaturated compacted clay with different initial densities were performed using a suction‐controllable triaxial apparatus, to measure the stress–strain–volume change in different stress paths and wetting paths. The model has well‐predicting capabilities to reproduce the mechanical behaviour of specimens compacted under different conditions not only in isotropic compression but also in triaxial compression and triaxial extension. Copyright © 2003 John Wiley & Sons, Ltd.     

大降雨条件下气压力对边坡稳定的影响研究

对于大面积浅层风化土边坡,当下部含有浅水位或不透水基岩层时大降雨将导致下部气体被封闭。随着湿润峰的下移,气压不断增大。封闭气压力不仅降低了雨水在边坡土体的入渗率,而且对边坡的稳定有显著影响。通过分析封闭气压力的形成和相关理论,提出取 大小的气压力头来研究边坡的稳定性（Hc为一水头值,与土体孔隙尺寸分布有关; hd为土体进气值水头）。结合非饱和土的Mohr-Coulomb破坏准则和极限平衡法,将封闭气压力引入到边坡的稳定分析中,建立了考虑气压力影响下的稳定分析模型。与传统的不考虑气压力的稳定分析方法作对比,提出了气压力影响率概念。研究表明,封闭气压力显著降低了边坡的安全系数,传统的无限边坡稳定计算方法偏于危险。研究结果对无限边坡的强降雨安全预报具有较好的指导作用。     

Evaluation of landslide triggering mechanisms in model fill slopes

Hong Kong is particularly susceptible to landslide risk due to the steep natural topography and prolonged periods of high intensity rainfall. Compounding the risk of slope failure is the existence of loose fill slopes which were constructed prior to the 1970s by end-tipping. A clear understanding of the underlying triggering mechanisms of fast landslides in fill slopes is required to analyse landslide risk and to optimise slope stabilisation strategies. The work described here had the objective of evaluating two candidate triggering mechanisms—static liquefaction and the transition from slide to flow due to localised transient pore water pressures—against observations of slope behaviour obtained from highly instrumented centrifuge model tests. These results indicate that static liquefaction is unlikely to occur if the model fill is unsaturated and the depth to bedrock large, as the high compressibility and mobility of air in the unsaturated void spaces allows the model fill slope to accommodate wetting collapse without initiating undrained failure. In contrast, high-speed failures with low-angle run-outs are shown to be easily triggered in model fill slopes from initially slow moving slips driven by localised transient pore water pressures arising from constricted seepage and material layering.     

An approach to model the mechanical behavior of transversely isotropic materials

The present work proposes an approach to adapt existing isotropic models to transversely isotropic materials. The main idea is to introduce equivalence relations between the real material and a fictitious isotropic one on which one can take all the advantages of the well‐established isotropic theory. Two applications of this approach are presented here: a failure criterion and a damage model that takes into account the load‐induced anisotropy. In both cases, theoretical predictions are in agreement with the experimental data. In the present paper, the developed approach is applied to sedimentary rock materials; nevertheless, it can be generalized to any material that exhibits transverse isotropy. Copyright © 2015 John Wiley & Sons, Ltd.     

A genetic classification of sinkholes illustrated from evaporite paleokarst exposures in Spain

This contribution analyses the processes involved in the generation of sinkholes from the study of paleokarst features exposed in four Spanish Tertiary basins. Bedrock strata are subhorizontal evaporites, and in three of the basins they include halite and glauberite in the subsurface. Our studies suggest that formation of dolines in these areas results from a wider range of subsidence processes than those included in the most recently published sinkhole classifications; a new genetic classification of sinkholes applicable to both carbonate and evaporite karst areas is thus proposed. With the exception of solution dolines, it defines the main sinkhole types by use of two terms that refer to the material affected by downward gravitational movements (cover, bedrock or caprock) and the main type of process involved (collapse, suffosion or sagging). Sinkholes that result from the combination of several subsidence processes and affect more than one type of material are described by combinations of the different terms with the dominant material or process followed by the secondary one (e.g. bedrock sagging and collapse sinkhole). The mechanism of collapse includes any brittle gravitational deformation of cover and bedrock material, such as upward stoping of cavities by roof failure, development of well-defined failure planes and rock brecciation. Suffosion is the downward migration of cover deposits through dissolutional conduits accompanied with ductile settling. Sagging is the ductile flexure of sediments caused by differential corrosional lowering of the rockhead or interstratal karstification of the soluble bedrock. The paleokarsts we analysed suggest that the sagging mechanism (not included in previous genetic classifications) plays an important role in the generation of sinkholes in evaporites. Moreover, collapse processes are more significant in extent and rate in areas underlain by evaporites than in carbonate karst, primarily due to the greater solubility of the evaporites and the lower mechanical strength and ductile rheology of gypsum and salt rocks.     

Late-glacial and post-glacial deposition in a large, low relief, epicontinental basin: the northern Baltic Sea

This paper presents a model of late‐glacial and post‐glacial deposition for the late‐Neogene sedimentary succession of the Archipelago Sea in the northern Baltic Sea. Four genetically related facies associations are described: (i) an ice‐proximal, acoustically stratified draped unit of glaciolacustrine rhythmites; (ii) an onlapping basin‐fill unit of rotated rhythmite clasts in an acoustically transparent to chaotic matrix interpreted as debris‐flow deposits; (iii) an ice‐distal, acoustically stratified to transparent, draped unit of post‐glacial lacustrine, weakly laminated to homogeneous deposits; and (iv) an acoustically stratified to transparent unit of brackish‐water, organic‐rich sediment drifts. The debris‐flow deposits of the unit 2 pass laterally into slide scars that truncate the unit 1; they are interpreted to result from a time interval of intense seismic activity due to bedrock stress release shortly after deglaciation of the area. Ice‐berg scouring and gravitational failure of oversteepened depositional slopes may also have contributed to the debris‐flow deposition. Comparisons to other late‐Neogene glaciated basins, such as the Hudson Bay or glacial lakes formed along the Laurentide ice sheet, suggest that the Archipelago Sea succession may record development typical for the deglaciation phase of large, low relief, epicontinental basins. The Carboniferous–Permian glacigenic Dwyka Formation in South Africa may provide an ancient analogue for the studied succession. Chronological control for the studied sediments is provided by the independent palaeomagnetic and AMS‐¹⁴C dating methods. In order to facilitate dating of the organic‐poor early post‐glacial deposits of the northern Baltic Sea, the 10 000 year long Lake Nautajärvi palaeomagnetic reference chronology ( Ojala & Saarinen, 2002 ) is extended by 1200 years.