NGI‐ADP: Anisotropic shear strength model for clay

Many geotechnical problems involve undrained behavior of clay and the capacity in undrained loading. Most constitutive models used today are effective stress based and only indirectly obtain values for the undrained shear strength. To match the design profiles of undrained shear strengths, in active (A), direct simple shear (D) and passive (P) modes of loading are complicated. This paper presents the elastoplastic constitutive model NGI‐ADP which is based on the undrained shear strength approach with direct input of shear strengths. Consequently, exact match with design undrained shear strengths profiles is obtained and the well‐known anisotropy of undrained shear strength and stiffness is accounted for in the constitutive model. A non‐linear stress path‐dependent hardening relationship is used, defined from direct input of failure strains in the three directions of shearing represented by triaxial compression, direct simple shear and triaxial extension. With its clear input parameters the model has significant advantages for design analysis of undrained problems. The constitutive model is implemented, into finite element codes, with an implicit integration scheme. Its performance is demonstrated by a finite element analysis of a bearing capacity problem. Copyright © 2011 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.     

A rigorous semi‐analytical solution for undrained cylindrical cavity expansion in critical state soils

This paper presents a generalized, rigorous and simple large strain solution for the undrained expansion of a vertical cylindrical cavity in critical state soils using a rate‐based plasticity formulation: the initial stress field is taken as anisotropic, that is with horizontal stresses that differ from the vertical stress, and the soil is assumed to satisfy any two‐invariant constitutive model from the critical state (Cam‐clay) family; no simplifying assumption is made during the mathematical derivation; calculating the effective stresses around the cavity requires the solution of a nonlinear equation by means of the Newton–Raphson method in combination with quadrature. Cavity expansion curves and stress distributions in the soil are then presented for different critical state models (including the modified Cam‐clay model). The solution derived can be useful for estimating the instantaneous response of saturated low‐permeability soils around piles and self‐boring pressuremeters and can serve as trustworthy benchmark for numerical analysis codes. Copyright © 2016 John Wiley & Sons, Ltd.     

Analysis of undrained cavity expansion in elasto‐plastic soils with non‐linear elasticity

A large strain analysis of undrained expansion of a spherical/cylindrical cavity in a soil modelled as non‐linear elastic modified Cam clay material is presented. The stress–strain response of the soil is assumed to obey non‐linear elasticity until yielding. A power‐law characteristic or a hyperbolic stress–strain curve is used to describe the gradual reduction of soil stiffness with shear strain. It is assumed that, after yielding, the elasto‐plastic behaviour of the soil can be described by the modified Cam clay model. Based on a closed‐form stress–strain response in undrained condition, a numerical solution is obtained with the aid of simple numerical integration technique. The results show that the stresses and the pore pressure in the soil around an expanded cavity are significantly affected by the non‐linear elasticity, especially if the soil is overconsolidated. The difference between large strain and small strain solutions in the elastic zone is not significant. The stresses and the pore pressure at the cavity wall can be expressed as an approximate closed‐form solution. Copyright © 2001 John Wiley & Sons, Ltd.     

Numerical predictions of the behaviour of soft clay with two anisotropic elastoplastic models

The numerical predictions obtained with the Melanie and MIT-E3 models are compared. Firstly, the performance of the constitutive models is checked against undrained triaxial tests. The models are then used to replicate the behaviour of an embankment built on soft clay. The numerical results are compared with the field data in terms of settlements, lateral displacements and excess pore pressures. Additionally, the numerical predictions are also analysed in terms of yield area, contours of vertical effective stresses, horizontal effective stresses and shear stresses and in terms of effective stress paths.     

Effects of transport of pore water and material heterogeneity on strain localization of fluid‐saturated gradient‐dependent viscoplastic geomaterial

The purpose of the present paper is to clarify the effects of permeability and initial heterogeneity on the strain localization of fluid‐saturated cohesive soil modelled by a strain gradient‐dependent poro‐viscoplastic constitutive model. The effects of permeability and gradient parameters on the growth rate of the fluctuation were obtained by a linear instability analysis. Deformation behaviour of clay specimens modelled as a viscoplastic model with a second order strain gradient during shear was numerically analysed by a soil–water coupled FEM under both globally undrained and partially drained conditions. It was found that the deformation pattern and the stress–strain curve greatly depend on the permeability, the drainage conditions and the initial non‐homogeneous properties. Copyright © 2005 John Wiley & Sons, Ltd.     

Drained and undrained pullout capacity of a stiff inclusion in a saturated poroelastic matrix

Shear‐lag analysis is used to obtain closed‐form solutions for the problem of a stiff inclusion embedded in a poroelastic soil matrix. The following assumptions are made: the soil matrix and the inclusion are elastic; plane strain conditions apply; and shear stresses at the soil‐inclusion interface follow Coulomb's friction law. Two solutions are obtained, the first one for drained conditions where no excess pore pressures are generated, and the second one for undrained conditions where excess pore pressures are produced and the soil does not change volume during pullout. The solutions are verified by comparing analytical predictions with numerical results obtained using a finite element method. Predictions from the analytical solutions are also compared with results from experiments conducted in a large‐scale pullout box. Both comparisons show good agreement. The analytical solution shows that the pullout capacity in drained and undrained conditions is overall independent of the relative stiffness of the soil and the inclusion. The most important factor controlling the pullout capacity is the coefficient of friction between the soil and the inclusion. Both drained and undrained pullout capacities increase with the coefficient of friction; although the drained capacity shows a proportional increase, it is not so for the undrained capacity. The ratio of undrained to drained pullout capacity is about 0.9 for friction coefficients smaller than 0.2, but can be as small as 0.6 for a coefficient of friction of 1.0. Copyright © 2006 John Wiley & Sons, Ltd.     

Modelling large plastic deformations of cohesive soils using sequential limit analysis

This paper introduces sequential limit analysis (SLA) as a method for modelling large plastic deformations of purely cohesive materials such as undrained clay. The method involves solving a series of consecutive small‐deformation plastic collapse problems using finite element limit analysis, thus ensuring high levels of accuracy, efficiency, and robustness. The techniques needed to develop an SLA implementation for two‐dimensional (plane strain) problems are described in detail, including model geometry updating routines, treatment of rigid bodies, interfaces and boundaries, and periodic remeshing and interpolation of field variables. A simple total stress‐based constitutive model is used to account for strain softening and strain rate effects. Extensive verifications and validations are performed using analytical solutions and physical model test results, comparing both collapse loads and failure mechanisms, to demonstrate the effectiveness of the SLA approach. Additional solution quality checks on the bracketing discrepancy between lower‐bound and upper‐bound limit analysis solutions, and on the incompressibility of the rigid‐plastic material, are also presented.     

A combined neural network/gradient‐based approach for the identification of constitutive model parameters using self‐boring pressuremeter tests

This paper presents a numerical procedure of material parameter identification for the coupled hydro‐mechanical boundary value problem (BVP) of the self‐boring pressuremeter test (SBPT) in clay. First, the neural network (NN) technique is applied to obtain an initial estimate of model parameters, taking into account the possible drainage conditions during the expansion test. This technique is used to avoid potential pitfalls related to the conventional gradient‐based optimization techniques, considered here as a corrector that improves predicted parameters. Parameter identification based on measurements obtained through the pressuremeter expansion test and two types of holding tests is illustrated on the Modified Cam clay model. NNs are trained using a set of test samples, which are generated by means of finite element simulations of SBPT. The measurements obtained through expansion and consolidation tests are normalized so that NN predictors operate independently of the testing depth. Examples of parameter determination are demonstrated on both numerical and field data. The efficiency of the combined parameter identification in terms of accuracy, effectiveness and computational effort is also discussed. Copyright © 2008 John Wiley & Sons, Ltd.     

正常固结黏土中平板锚基础的吸力和抗拉力

平板锚是新近出现的一种系泊深海浮式结构的基础型式。当黏土地基中的平板锚承受上拔力时,平板上、下表面超静孔压差形成的吸力使其抗拉承载力显著增加,对于风浪等快速加载条件尤其如此。利用有限元软件ADINA建立有效应力形式的轴对称动力有限元模型,研究圆形平板锚在缓慢加载与快速加载时的超静孔压分布与地基破坏型式。快速加载时重黏土和高岭土两组典型正常固结土样所得极限承载力系数与塑性极限分析解一致。进而通过变动参数分析,讨论加载速率和埋深对吸力和总抗拉力极限值的影响,并给出排水和不排水加载条件对应的临界加载速率。结果表明,不排水加载条件的总抗拉力可能达到排水总抗拉力的3倍。     

Stress-strain behaviour of reconstituted illitic clay at different temperatures

Tests on specimens of reconstituted illitic clay have examined the influence of temperature on the mechanical behaviour of clay soils. The program involved consolidation to effective confining pressures up to 1.5 MPa, heating to 100°C, and tests on normally consolidated and overconsolidated specimens with OCR = 2. The tests included isotropic consolidation, undrained triaxial compression with pore water pressure measurement, drained tests along controlled stress paths to investigate yielding behaviour, and undrained tests which involved heating and measurement of the resulting induced pore water pressures. The large strain strength envelope is independent of temperature. However, peak undrained strengths increase with temperature because smaller pore water pressures are generated during shearing. An important contribution from the study is a series of results for the yielding of illitic clay at three different temperatures. For the first time, there is clear evidence of yield loci decreasing in size with increasing temperature. An associated flow rule can be assumed without serious error. The results contribute to the confirmation of a thermal elastic-plastic soil model developed by the authors from cam clay following the addition of a small number of extra assumptions. Depending on the initial stress state, heating under undrained conditions may produce shear failure.     

Sensitive bounding surface constitutive model for structured clay

The main purpose of the paper is to present a relatively simple, yet realistic, constitutive model for simulations of structured sensitive clays. The proposed constitutive model can simulate 1‐D and isotropic consolidation, and drained and undrained shear response of sensitive structured clay. The proposed sensitive bounding surface model is based on concepts from the modified Cam clay model 8 and bounding surface plasticity 27 , with the addition of a simple degradation law. The key material parameters are M, λ, κ, and ν from the modified Cam clay framework, h from the bounding surface framework to model a smoothed elasto‐plastic transition, and ω_v, ω_q, and S_sr to model softening associated with destructuration. The model has separate parameters to model destructuration caused by volumetric strain and deviatoric strain. The model is capable of modeling unusual behavior of strain softening during 1‐D compression (i.e., a reduction of effective stress as void ratio decreases). A good match between test results and the model simulation is demonstrated. Copyright © 2016 John Wiley & Sons, Ltd.     

Deformation analysis of shallow penetration in clay

A new method of analysis is described for estimating the deformations and strains caused by shallow undrained penetration of piles and caissons in clay. The formulation combines previous analyses for steady, deep penetration, with methods used to compute soil deformations due to near-surface ground loss, and is referred to as the Shallow Strain Path Method (SSPM). Complete analytical solutions for the velocity and strain rates are given for a planar wall, an axisymmetric, closed-ended pile and unplugged, open-ended pile geometries. In these examples, the analyses consider a single source penetrating through the soil at a constant rate, generating a family of penetrometers with rounded tips, referred to as simple wall, pile and tube geometries. Soil deformations and strains are obtained by integrating the velocity and strain rates along the particle paths. The transition from shallow to deep penetration is analysed in detail. Shallow penetration causes heave at the ground surface, while settlements occur only in a thin veneer of material adjacent to the shaft and in a bulb-shaped region around the tip. The size of this region increases with the embedment depth. Deformations inside an open-ended pile/caisson are affected significantly by details of the simple tube wall geometry. © 1997 by John Wiley & Sons, Ltd.     

上海第4层黏土力学特性的弹塑性本构模拟

上海第4层黏土是典型的结构性海相软土,用一个本构模型统一地模拟不同应力路径下的力学特性对数值计算具有重要意义。对UNIFIED模型的结构性及超固结发展函数进行了改进,并提出了一种确定原状土材料参数和初始状态的方法。为了验证修正模型的正确性,用块状取土法取得上海第4层原状土样,进行了固结试验和三轴排水、不排水剪切试验。通过比较试验结果和本构模拟结果,明确了新的本构模型仅用一组材料参数就能统一地模拟上海第4层黏土在固结、排水及不排水三轴试验得到的应力-应变关系。模拟结果揭示了上海第4层黏土的结构比较稳定,即使在经历三轴剪切发生35%轴应变后仍能保持较高位的结构性。     

冲击作用下饱和土性状的试验研究

在三轴条件下,对饱和土（砂土和黏土）进行排水与不排水条件下的冲击试验及冲击后再固结试验,对比研究了不同渗透性土在不同排水条件下的冲击动力响应和冲击后再固结性状。结果表明：饱和黏土不排水冲击时的孔隙水压力随冲击击数增加而升高并逐渐稳定,排水冲击时的孔隙水压力则是先达到峰值然后有所下降;砂土不排水冲击时的冲击能量对孔隙水压力影响最明显;饱和砂土不排水冲击时的轴向应变与冲击击数呈近似线性关系,饱和黏土冲击及饱和砂土排水冲击则呈近二次曲线关系;饱和砂土不排水冲击后再固结阶段的孔隙水压力立即消散为0,同时体变迅速增大到一定值;饱和黏土在冲击后再固结阶段的孔隙水压力在一定时间内逐渐消散完毕,同时体变逐渐增大;饱和黏土排水冲击时,冲击阶段产生的体变占冲击引起总体变的39%～49%,冲击后再固结阶段产生的体变占51%～61%;砂土和黏土的总体变均表现为排水冲击明显大于不排水冲击,改善冲击时的排水条件有利于提高加固效果。     

Simulation of yielding and stress–stain behavior of shanghai soft clay

In this paper, a simple bounding surface plasticity model is used to reproduce the yielding and stress–strain behavior of the structured soft clay found at Shanghai of China. A series of undrained triaxial tests and drained stress probe tests under isotropic and anisotropic consolidation modes were performed on undisturbed samples of Shanghai soft clay to study the yielding characteristics. The degradation of the clay structure is modeled with an internal variable that allows the size of the bounding surface to decay with accumulated plastic strain. An anisotropic tensor and rotational hardening law are introduced to reflect the initial anisotropy and the evolution of anisotropy. Combined with the isotropic hardening rule, the rotational hardening rule and the degradation law are incorporated into the bounding surface formulation with an associated flow rule. Validity of the model is verified by the undrained isotropic and anisotropic triaxial test and drained stress probe test results for Shanghai soft clay. The effects of stress anisotropy and loss of structure are well captured by the model.     

SBPT测定饱和黏土不排水强度的数值分析

自钻式旁压试验（SBPT）因其扰动小、测试深度大、可以获得应力-应变、超孔隙水压力-时间等数据,在确定地基土性参数和地基承载力上有广阔的应用前景。然而由于用以解释SBPT的柱孔扩张理论（Gibson解）所采用的平面应变假设与实际旁压腔几何特征存在差异,导致试验所确定的黏土不排水剪切强度su与其他原位试验或室内试验结果存在差别。针对旁压腔几何尺寸及应变区间的选择对确定su的影响,基于修正剑桥模型,采用低渗透系数控制加载过程中不排水条件,利用有限元法模拟SBPT,建议了不同应力历史下确定su的应变区间,并给出考虑几何尺寸影响时相应应变区间上su的修正系数。     

饱和软黏土动力学特性循环扭剪试验研究

以饱和软黏土的循环扭剪试验研究为基础,阐明不固结不排水条件下饱和软黏土的动力学特性。在循环扭剪和循环三轴两种不同试验应力状态下,通过研究不同围压和不同静、动应力组合下饱和软黏土的应力等效破坏关系和应变等效破坏关系,提出在循环荷载作用下饱和软黏土的循环破坏同样遵循Mises屈服准则,且在Mises屈服准则下,饱和软黏土的循环强度趋于不变量。循环破坏的过程可以等效为一种拟静力弹塑性循环蠕变,建立了饱和软黏土循环累积变形随静荷载、循环荷载以及循环破坏振次之间的关系式。上述结论分别从应力和变形两个方面阐述饱和软黏土的动力学特性,与试验应力状态和围压无关,可以推广到一般应力状态下。     

循环荷载对粉质黏土力学性质的影响

通过等压固结静、动三轴试验,研究了渤海湾粉质黏土在循环荷载作用下的动力性质和循环荷载后不排水静力性质。试验结果表明,循环应力幅值比越大,平均轴向应变和轴向应变幅值越大;循环应力幅值比达到0.4时,平均轴向应变和轴向应变幅值随着循环周数增加迅速增大;循环应力幅值比相同,固结应力越大,轴向应变幅值越大,而平均轴向应变越小。在较大的循环应力幅值比下,平均孔压比值和孔压幅值比值随着循环周数的增大会达到稳定;循环应力幅值比越大平均孔压比值和孔压幅值比值均越大;相同循环应力幅值,固结应力越大平均孔压比值越小,而固结应力对孔压幅值比值影响较小。循环荷载的作用会导致循环荷载后不排水剪在q-p’平面上有效应力路径和孔压发展表现出超固结土的性质。     

Field measurements of shear strength of an underconsolidated marine clay

This paper presents the observations of cone penetration testing (CPT), in situ vane shear testing and undrained triaxial testing of underconsolidated marine clay in the Craney Island Dredged Material Management Area (CIDMMA), Norfolk, VA. Comparison of vane shear measurements, taken during 1981–1983, to CPT measurements, taken during 1993, indicates a change in undrained shear strength. Both the vane shear and cone penetration resistances are lowest at the mid-heights of the clay layers and the excess pore water pressures are highest at the mid-height of the clay layers, indicating that the clay layer is underconsolidated.