Modelling the behaviour of an embankment on soft clay with different constitutive models

The paper investigates the effect of constitutive models on the predicted response of a simplified benchmark problem, an embankment on soft soil. The soft soil is assumed to have the properties of POKO clay from Finland and five different constitutive models are used to model the deposit. Two of the models are isotropic models, i.e. the Modified Cam Clay model and the Soft‐Soil model. The other models are recently proposed constitutive models that account for plastic anisotropy. The S‐CLAY1 and S‐CLAY1S models are embedded in a standard elasto‐plastic framework and account for anisotropy via a rotational hardening law. In addition, the S‐CLAY1S model accounts for bonding and destructuration. In contrast, the Multilaminate Model for Clay (MMC) accounts for plastic anisotropy by utilizing so‐called multilaminate framework. The results of numerical simulations show that accounting for anisotropy results in notable differences in the predicted settlements and horizontal movements compared to the predictions using the isotropic models. There are also significant differences in the K₀ predictions by the different constitutive models and this has a significant impact on the results. Copyright © 2006 John Wiley & Sons, Ltd.     

Failure envelopes and plastic potentials for eccentrically loaded surface footings on undrained soil

The failure envelope in V–M space for surface foundations on undrained material under eccentric loading can be determined using an extended version of the scaling (or effective area) concept of Meyerhof. A similar displacement transformation allows production of the plastic potential. The two‐dimensional finite element analyses of fully attached foundations subject to combined vertical (V) and moment (M) loading have been used to calculate appropriate scaling points for deduction of the failure envelope and plastic potential. Failure envelopes and plastic potentials are presented for footings on uniform and non‐uniform undrained material and it is seen that the equivalent ‘critical state’ or ‘parallel point’ lies slightly beneath the peak moment capacity. For accurate prediction of failure envelopes for footings on non‐uniform strength soil, consideration must be made of the apparent reduction of the soil heterogeneity as the area of footing in contact with the soil decreases. Copyright © 2001 John Wiley & Sons, Ltd.     

A non‐coaxial critical state soil model and its application to simple shear simulations

The yield vertex non‐coaxial theory is implemented into a critical state soil model, CASM (Int. J. Numer. Anal. Meth. Geomech. 1998; 22 :621–653) to investigate the non‐coaxial influences on the stress–strain simulations of real soil behaviour in the presence of principal stress rotations. The CASM is a unified clay and sand model, developed based on the soil critical state concept and the state parameter concept. Without loss of simplicity, it is capable of simulating the behaviour of sands and clays within a wide range of densities. The non‐coaxial CASM is employed to simulate the simple shear responses of Erksak sand and Weald clay under different densities and initial stress states. Dependence of the soil behaviour on the Lode angle and different plastic flow rules in the deviatoric plane are also considered in the study of non‐coaxial influences. All the predictions indicate that the use of the non‐coaxial model makes the orientations of the principal stress and the principal strain rate different during the early stage of shearing, and they approach the same ultimate values with an increase in loading. These ultimate orientations are dependent on the density of soils, and independent of their initial stress states. The use of the non‐coaxial model also softens the shear stress evolutions, compared with the coaxial model. It is also found that the ultimate shear strengths by using the coaxial and non‐coaxial models are dependent on the plastic flow rules in the deviatoric plane. Copyright © 2006 John Wiley & Sons, Ltd.     

重塑硅藻土抗剪强度的环剪试验研究

硅藻土是一种强结构性土,具有密度小、孔隙大、干燥土块吸水率高等特点。天然硅藻土的抗剪强度大,但在外界扰动或风化后,结构性强度丧失,易发生滑坡等斜坡变形破坏现象,因此,研究其在大变形下的力学性质很有必要。利用环剪仪可以研究土体在较大剪切位移条件下抗剪强度变化规律,对重塑硅藻土进行环剪试验研究,得到了不同竖向压力下的应力-位移曲线,认为重塑硅藻土为典型的应变软化型土;根据测试结果得到了重塑硅藻土的抗剪强度和残余强度,并从物质组成和微观结构特征方面对重塑硅藻土的力学性质进行了分析,为硅藻土滑坡的成因分析和灾害防治提供了依据。     

Anisotropic viscoplastic modeling of rate‐dependent behavior of clay

The objective of this study is to derive an effective stress‐based constitutive law capable of predicting rate‐dependent stress–strain, stress path and undrained shear strength and creep behavior. The flow rule used in the MIT‐E3 model and viscoplasticity theory is employed in the derivation. The model adopts the yield surface capable of representing the yield behavior of the Taipei silty clay and assumes that it is initially symmetric about the K₀‐line. A method is then developed to compute the gyration and expansion of the loading surface to simulate the anisotropic behavior due to the principal stress rotation after shear. There are 11 parameters required for the model to describe the soil behavior and six of them are exactly the same as those used in the Modified Cam‐clay model. The five additional parameters can be obtained by parametric studies or conventional soil tests, such as consolidation tests, triaxial compression and extension tests. Finally, verification of the model for the anisotropic behavior, creep behavior and the rate‐dependent undrained stress–strain and shear strength of the Taipei silty clay is conducted. Copyright © 2010 John Wiley & Sons, Ltd.     

软土中T型触探仪贯入阻力的数值模拟

T型触探仪近年来在国外的海洋软土地基原位测试中得到了广泛应用,但贯入速率、地基土强度的各向异性与渐进软化等因素对于贯入阻力的影响没有得到系统的研究。在大型有限元软件ABAQUS平台上进行二次开发,针对基于Tresca屈服准则的理想弹塑性模型进行了相应改进,使之可以比较合理地反映上述因素的影响,进而对软土中T型触探仪的贯入机制进行了比较系统的数值分析。计算结果表明,贯入速率、地基土强度的各向异性与强度软化效应对于T型触探仪的贯入阻力系数影响较大。通过与有关极限分析上限解的对比分析,在一定程度上验证了有限元分析结果的合理性。     

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.     

Discrete element method modeling of inherently anisotropic rocks under uniaxial compression loading

A new numerical approach is proposed in this study to model the mechanical behaviors of inherently anisotropic rocks in which the rock matrix is represented as bonded particle model, and the intrinsic anisotropy is imposed by replacing any parallel bonds dipping within a certain angle range with smooth‐joint contacts. A series of numerical models with β = 0°, 15°, 30°, 45°, 60°, 75°, and 90° are constructed and tested (β is defined as the angle between the normal of weak layers and the maximum principal stress direction). The effect of smooth‐joint parameters on the uniaxial compression strength and Young's modulus is investigated systematically. The simulation results reveal that the normal strength of smooth‐joint mainly affects the behaviors at high anisotropy angles (β > 45°), while the shear strength plays an important role at medium anisotropy angles (30°–75°). The normal stiffness controls the mechanical behaviors at low anisotropy angles. The angle range of parallel bonds being replaced plays an important role on defining the degree of anisotropy. Step‐by‐step procedures for the calibration of micro parameters are recommended. The numerical model is calibrated to reproduce the behaviors of different anisotropic rocks. Detailed analyses are conducted to investigate the brittle failure process by looking at stress‐strain behaviors, increment of micro cracks, initiation and propagation of fractures. Most of these responses agree well with previous experimental findings and can provide new insights into the micro mechanisms related to the anisotropic deformation and failure behaviors. The numerical approach is then applied to simulate the stress‐induced borehole breakouts in anisotropic rock formations at reduced scale. The effect of rock anisotropy and stress anisotropy can be captured. Copyright © 2015 John Wiley & Sons, Ltd.     

Elastic models for nonlinear response of rigid passive piles

Recent study indicates that the response of rigid passive piles is dominated by elastic pile–soil interaction and may be estimated using theory for lateral piles. The difference lies in that passive piles normally are associated with a large scatter of the ratio of maximum bending moment over maximum shear force and induce a limiting pressure that is ~1/3 that on laterally loaded piles. This disparity prompts this study. This paper proposes pressure‐based pile–soil models and develops their associated solutions to capture response of rigid piles subjected to soil movement. The impact of soil movement was encapsulated into a power‐law distributed loading over a sliding depth, and load transfer model was adopted to mimic the pile–soil interaction. The solutions are presented in explicit expressions and can be readily obtained. They are capable of capturing responses of model piles in a sliding soil owing to the impact of sliding depth and relative strength between sliding and stable layer on limiting force prior to ultimate state. In comparison with available solutions for ultimate state, this study reveals the 1/3 limiting pressure (of the active piles) on passive piles was induced by elastic interaction. The current models employing distributed pressure for moving soil are more pertinent to passive piles (rather than plastic soil flow). An example calculation against instrumented model piles is provided, which demonstrates the accuracy of the current solutions for design slope stabilising piles. Copyright © 2014 John Wiley & Sons, Ltd.     

Modelling chemo‐hydro‐mechanical behaviour of unsaturated clays: a feasibility study

Effective capabilities of combined chemo‐elasto‐plastic and unsaturated soil models to simulate chemo‐hydro‐mechanical (CHM) behaviour of clays are examined in numerical simulations through selected boundary value problems. The objective is to investigate the feasibility of approaching such complex material behaviour numerically by combining two existing models. The chemo‐mechanical effects are described using the concept of chemical softening consisting of reduction of the pre‐consolidation pressure proposed originally by Hueckel (Can. Geotech. J. 1992; 29 :1071–1086; Int. J. Numer. Anal. Methods Geomech. 1997; 21 :43–72). An additional chemical softening mechanism is considered, consisting in a decrease of cohesion with an increase in contaminant concentration. The influence of partial saturation on the constitutive behaviour is modelled following Barcelona basic model (BBM) formulation (Géotech. 1990; 40 (3):405–430; Can. Geotech. J. 1992; 29 :1013–1032). The equilibrium equations combined with the CHM constitutive relations, and the governing equations for flow of fluids and contaminant transport, are solved numerically using finite element. The emphasis is laid on understanding the role that the individual chemical effects such as chemo‐elastic swelling, or chemo‐plastic consolidation, or finally, chemical loss of cohesion have in the overall response of the soil mass. The numerical problems analysed concern the chemical effects in response to wetting of a clay specimen with an organic liquid in rigid wall consolidometer, during biaxial loading up to failure, and in response to fresh water influx during tunnel excavation in swelling clay. Copyright © 2005 John Wiley & Sons, Ltd.     

A coupled two‐phase fluid flow and elastoplastic deformation model for unsaturated soils: theory,implementation, and application

Although numerous numerical models have been proposed for simulating the coupled hydromechanical behaviors in unsaturated soils, few studies satisfactorily reproduced the soil–water–air three‐phase coupling processes. Particularly, the impacts of deformation dependence of water retention curve, bonding stress, and gas flow on the coupled processes were less examined within a coupled soil–water–air model. Based on our newly developed constitutive models (Hu et al., 2013, 2014, 2015) in which the soil–water–air couplings have been appropriately captured, this study develops a computer code named F²Mus3D to investigate the coupled processes with a focus on the above impacts. In the numerical implementation, the generalized‐α time integration scheme was adopted to solve the equations, and a return‐mapping implicit stress integration scheme was used to update the state variables. The numerical model was verified by two well‐designed laboratory tests and was applied for modeling the coupled elastoplastic deformation and two‐phase fluid flow processes in a homogenous soil slope induced by rainfall infiltration. The simulation results demonstrated that the numerical model well reproduces the initiation of a sheared zone at the toe of the slope and its propagation toward the crest as the rain infiltration proceeds, which manifests a typical mechanism for rainfall‐induced shallow landslides. The simulated plastic strain and deformation would be remarkably underestimated when the bonding stress and/or the deformation‐dependent nature of hydraulic properties are ignored in the coupled model. But on the contrary, the negligence of gas flow in the slope soil results in an overestimation of the rainfall‐induced deformation. Copyright © 2015 John Wiley & Sons, Ltd.     

基于微观力学机制的各向异性结构性砂土的本构模型研究

在岩土破损力学基础上,基于微观破损机制,提出了考虑各向异性的结构性砂土本构理论。采用Lade-Duncan强度准则考虑中主应力对抗剪强度的影响;采用考虑颗粒排列组构的各向异性状态变量A反映各向异性对土体强度和变形的影响;通过相似扩大重塑土的屈服面反映结构性对土性的影响;通过引入非相关联流动法则考虑各向异性和结构性对土体塑性变形的影响。同时,将基于微观力学机制的损伤演化规律引入结构性土的硬化规律;该硬化规律同时考虑了塑性体积应变和剪切应变对各向异性结构性土强度的影响。然后将该模型用于模拟室内三轴压缩试验,初步验证了该模型的合理性和适用性。     

A general method for defining the number of cycles of repeated loading

In geotechnical engineering designs and safety checks, the number of cycles of repeated loading N is frequently used as a parameter for calculating the response of soil structures under the repeated loading. However, there appears to be a lack of rational methods for calculating N for repeated loadings along general stress paths. A new method for defining N is proposed in this paper, based on the assumption that a yield surface may be identified for the material. It is demonstrated that the proposed method provides a rational and useful means for engineering calculation of soil deformations under repeated loadings. Copyright © 2001 John Wiley & Sons, Ltd.     

Influence of initial plastic anisotropy of overconsolidated clays on ground behaviour during tunneling

In order to evaluate the influence of the initial plastic anisotropy in the excavation of a tunnel, a “bubble” bounding surface model for structured soils, formulated by Kavvadas and Belokas (Proc. 10th IACMAG Conf., 2001), was implemented in the explicit finite difference code FLAC. Two different initial stress (K ₀) conditions were considered. The size and shape of the initial bounding surfaces were specified to be consistent with the initial stress field. The distorted and rotated shape of the bounding surface, supported by experimental results, defines the anisotropy of shear strength, which is shown to have a significant influence on the displacements. There is also considerable sensitivity of the soil model to the initial stress field.     

A non‐coaxial and unbalanced plastic flow rule for geomaterials and its application to the shear band orientation prediction

In this paper, the non‐coaxial relation between the principal plastic strain increments and the principal stresses, which results from the internal friction in geomaterials, is analyzed, and the phenomenon of the unbalanced development of plastic flow in two conjugate directions is discussed. A non‐coaxial, unbalanced plastic flow model for Coulomb frictional materials is developed and used to determine the orientation of shear band in geomaterials. It is shown that the unbalanced index r of plastic flow has important effect on the orientation of the shear band, and the orientation determined by the conventional plastic flow theory is only a special case of the proposed model when r=0. This result soundly explains the reason that the geomaterials with the same internal friction angle and dilatancy angle can have very different shear band orientations. In addition, the difference between the intrinsic and apparent dilatancy angles is analyzed, and it is emphasized that the dilatancy angle commonly used in practice is indeed the apparent dilatancy angle. Copyright © 2010 John Wiley & Sons, Ltd.     

Effect of inherent anisotropy on the strength of natural granite residual soil under generalized stress paths

It is universally known that residual soils behave very differently from sedimentary soils. While the latter is widely known as cross-anisotropic, little is known regarding the strength anisotropy of residual soils. This study presents how the inherent anisotropy affects the strength of natural granite residual soils under generalized conditions, where intact specimens were carefully prepared and sheared under triaxial compression, extension, simple shear, and hollow cylinder torsional shear tests. The strength of natural residual soil, in terms of ultimate stress ratio M and undrained shear strength S_u, is found to be significantly anisotropic in a different way from normally consolidated clays with the maximum strength obtained under triaxial compression and the minimum under simple shear or at intermediate principal stress direction. As a result, the existing method failed to measure the anisotropy degree of the studied soil. Two parameters were proposed accordingly to quantify the anisotropic strength under general conditions, taking the special strength anisotropy pattern and cohesive-frictional nature of GRS into account. The proposed parameters enable the direct comparison of strength anisotropy among soils. This study serves as a data set to better understand residual soils regarding their anisotropic behaviors under generalized conditions. Although specific to granite residual soils in China, this study is expected to be more widely applicable to other weathered geomaterials.

     

Quantitative Expression of Heat Flow versus Tectonic Deformation in the China Continent： The Effects of Plastic-Flow Network and Stable Block

1 Introduction It has been understood in the study of terrestrial heat flow that the distribution of heat flow in the interior of continent is influenced by a large number of factors, involving heat sources (e.g. mantle heat flow, heat production of radioactive elements in the crust, magmatic activity, and heat production of tectonic deformation), heat transfercondition (e.g. thermal conductivity and thickness of media), groundwater circulation, etc. On the background of these factors it is pa…     

Development of hyperplasticity models for soil mechanics

Hyperplasticity theory was developed by Collins and Houlsby (Proc. Roy. Soc. Lon. A 1997; 453 :1975–2001) and Houlsby and Puzrin (Int. J. Plasticity 2000; 16 (9):1017–1047). Further research has extended the method to continuous hyperplasticity, in which smooth transitions between elastic and plastic behaviour can be modelled. This paper illustrates a development of a new constitutive model for soils using hyperplasticity theory. The research begins with a simple one‐dimensional elasticity model. This is extended in stages to an elasto‐plastic model with a continuous internal function. The research aims to develop a soil model, which addresses some of the shortcomings of the modified cam‐clay model, specifically the fact that it cannot model small strain stiffness, or the effects of immediate stress history. All expressions used are consistent with critical state soil mechanics terminology. Finally, a numerical implementation of the model using a rate‐dependent algorithm is described. Copyright © 2005 John Wiley & Sons, Ltd.     

An anisotropic elastoplastic model for soft clays based on logarithmic contractancy

A new constitutive model for soft structured clays is developed based on an existing model called S‐CLAY1S, which is a Cam clay type model that accounts for anisotropy and destructuration. The new model (E‐SCLAY1S) uses the framework of logarithmic contractancy to introduce a new parameter that controls the shape of the yield surface as well as the plastic potential (as an assumed associated flow rule is applied). This new parameter can be used to fit the coefficient of earth pressure at rest, the undrained shear strength or the stiffness under shearing stress paths predicted by the model. The improvement to previous constitutive models that account for soil fabric and bonding is formulated within the contractancy framework such that the model predicts the uniqueness of the critical state line and its slope is independent of the contractancy parameter. Good agreement has been found between the model predictions and published laboratory results for triaxial compression tests. An important finding is that the contractancy parameter, and consequently the shape of the yield surface, seems to change with the degree of anisotropy; however, further study is required to investigate this response. From published data, the yield surface for isotropically consolidated clays seems ‘bullet’ or ‘almond’ shaped, similar to that of the Cam clay model; while for anisotropically consolidated clays, the yield surface is more elliptical, like a rotated and distorted modified Cam clay yield surface. © 2015 The Authors. International Journal for Numerical and Analytical Methods in Geomechanics published by John Wiley & Sons Ltd.