Load‐displacement and bearing capacity of foundations on granular soils using a multi‐surface kinematic constitutive soil model

A finite element approach based on an advanced multi‐surface kinematic constitutive model is used to evaluate the bearing capacity of footings resting on granular soils. Unlike simple elastic‐perfectly plastic models, often applied to granular foundation problems, the present model realistically accounts for stress dependency of the friction angle, strain softening–hardening and non‐associativity. After the model and its implementation into a finite element code are briefly discussed, the numerical difficulty due to the singularity at the footing edge is addressed. The bearing capacity factor N_γ is then calculated for different granular materials. The effect of footing size, shape, relative density and roughness on the ultimate bearing capacity are studied and the computed results compare very favourably with the general experimental trends. In addition, it is shown that the finite element solution can clearly represent counteracting mechanisms of progressive failure which have an important effect on the bearing capacity of granular foundations. Copyright © 2006 John Wiley & Sons, Ltd.     

不饱和土动弹塑性本构模型研究

以饱和度与有效应力为状态变量,通过引入描述不饱和与饱和土孔隙比差的状态变量,将Zhang等提出的饱和土体应力诱导各向异性动弹塑性本构模型推广到不饱和土体中,使其可描述不饱和土在动力循环荷载作用下的力学特性行为。通过对已有不饱和土体在完全不排水条件下的动三轴试验进行理论模拟,验证了所提出不饱和土本构模型的正确性。最后基于所提出本构模型,讨论了在不排水条件下初始饱和度对不饱和土动力特性研究。结果表明,不饱和土在动力荷载作用下,土体的孔隙比将减少,导致饱和度增加;当初始饱和度较高时,不饱和土会转化为饱和土,从而发生液化现象。该研究成果对研究不饱和土在地震等动力荷载作用下的力学特性行为具有重要意义。     

Application of unconstrained optimization and sensitivity analysis to calibration of a soil constitutive model

Large sets of soil experimental data (field and laboratory) are becoming increasingly available for calibration of soil constitutive models. A challenging task is to calibrate a potentially large number of model parameters to satisfactorily match many data sets simultaneously. This calibration effort can be facilitated by optimization techniques. The current study aims to explore systematic approaches for exercising optimization and sensitivity analysis in the area of soil constitutive modelling. Analytical, semi‐analytical and numerical optimization techniques are employed to calibrate a multi‐surface‐plasticity sand model. Calibration is based on results from a number of drained triaxial sample tests and a dynamic centrifuge liquefaction test. The analytical and semi‐analytical approaches and associated sensitivity analysis are applied to calibrate the model non‐linear shear stress–strain response. Thereafter, model parameters controlling shear–volume coupling effects (dilatancy) are calibrated using a solid–fluid fully coupled finite element program in conjunction with an advanced numerical optimization code. A related sensitivity study reveals the challenges often encountered in optimizing highly non‐linear functions. Overall, this study demonstrates applicability and limitations of optimization techniques for constitutive model calibration. Copyright © 2003 John Wiley & Sons, Ltd.     

A constitutive model for unsaturated soil interfaces

A constitutive model based on the disturbed state concept is presented to describe the behavior of interfaces in unsaturated soil. The model is an extension of an existing model developed for a sand–steel interface. As opposed to the original model, the modified model incorporates two independent stress variables, which are the net normal stress and matric suction. The saturated and dry state of the interface can be modeled as a special case using the constitutive model presented in this paper. The modified model is capable of capturing the main features of unsaturated interfaces observed during laboratory testing, including increasing shear strength and strain softening with increasing suction and net normal stress and increasing dilatancy with increasing suction. Laboratory tests were carried out on unsaturated interfaces in a modified direct shear test apparatus. The observed behavior of interfaces between unsaturated soil and steel plates (rough and smooth) is presented in comparison with model predictions. Copyright © 2008 John Wiley & Sons, Ltd.     

在ABAQUS中开发实现Duncan-Chang本构模型

美国HKS公司开发的软件产品ABAQUS是目前世界上最强大的非线性有限元分析工具之一。该软件包括了众多材料本构模型,但尚缺少国内外土工数值分析中广泛采用的Duncan-Chang本构模型。这一缺憾影响了ABAQUS软件在土工分析中的应用范围。本文旨在介绍ABAQUS中开发Duncan-Chang材料本构模型实现方法,给出了开发过程的概要,完成了两个典型常规三轴压缩模型问题数值测试。结果表明：在ABAQUS中增加Duncan-Chang材料本构模型后,不仅可以充分地利用该软件强大的非线性求解平台,而且还有可能完成复杂土工应力应变的有限元数值分析问题,同时,具有计算速度快、计算精度高和前后处理快捷方便的优点,极大地降低了土工分析程序开发的难度,并减少了维护工作量。     

一种基于修正西原模型的冻土蠕变本构关系

冻土在高应力水平下通常会表现出加速蠕变特性。经典西原（Nishihara）模型可以反映衰减蠕变和稳态蠕变,但不适用于描述加速蠕变。通过对Nishihara模型中黏性元件进行修正,将黏弹性元件中的定常的黏滞系数修正为时间的函数,并在黏塑性元件中采用非定常、非牛顿黏性元件,得到的蠕变柔量可以考虑应力水平和时间的影响,使其能够反映其加速蠕变的过程。用提出的模型分别对一维和三维应力、不同应力水平下的冻土蠕变试验数据进行了拟合,拟合结果能够反映冻土从衰减蠕变和稳定蠕变状态进入加速蠕变的过程,说明了该模型对于描述冻土在不同应力下不同蠕变状态的适用性。     

土体非线弹性-塑性本构模型

把本质上属于亚弹性本构模型,在岩土工程界广为应用的Duncan-Chang模型与服从Drucker-PragerMohr-Column屈服准则的弹塑性本构模型相结合,推出了非线弹性-塑性的组合本构模型,以克服经典的弹塑性模型不能考虑岩土材料在塑性屈服前的非线性行为以及一般的Duncan-Chang模型不适用于应力水平接近于屈服或破坏状态等缺点。所建议的本构模型简单实用,能较好地综合利用工程地质勘察资料和常规的土工试验数据。针对实际边坡具体地层条件的分析表明,该模型及其非线性析算法的正确性和可行性,而且在岩土工程方面具有广泛的实际应用价值。     

考虑试验阻尼效应的一种土体动力双型抛物线本构模型

针对土动力试验中土体反向卸载（或加载）曲线形状有时表现为先变化平缓后变化剧烈这一特征,基于土体实验结果的 和 曲线,建立了不规则荷载作用下可调的双型抛物线形式的本构曲线规则,阐述了A型和B型抛物线的函数表达式及实现过程。采用9种方法对试验阻尼比的模拟对比,表明所建立的加卸载准则尤其能准确拟合试验阻尼值,是对广义Masing准则的一种简化及实用化。由于调制因子 的可变性,它具有更强地模拟试验本构曲线形状的能力。针对台湾Lotung DHB钻井台阵试验场地,以等效线性化方法SHAKE91、LSSRLI-1,时域非线性方法DESRA-2、Pyke方法、非Masing准则模型1、“隐式应力阻尼等效”模型、“阻尼比退化系数”模型和非Masing准则模型2以及双型抛物线本构模型等9种方法（或模型）分别进行了非线性地震反应分析,并将结果与强震观测记录进行比较,结果表明,双型抛物线本构模型计算的复合加速度地震动的峰值大小、波形大小相对关系及后续波形,与实际地震记录较相一致,说明了模型的合理性以及实际工程应用的可行性。     

水泥土细观破裂过程的损伤本构模型

根据对饱水环境中水泥土力学特性的宏观试验和细观破裂过程试验研究,建立了水泥土细观孔隙损伤变量和损伤本构模型,该水泥土的损伤变量不仅与水泥土的初始损伤有关,而且与水泥土微小孔隙中的孔隙水的密度有关,所建立的损伤本构模型与饱水情况下水泥土的单轴压缩试验曲线进行了对比,分析结果表明,所提出的水泥土细观孔隙损伤本构模型与试验值基本吻合,很好地描述了水泥土的损伤演化特性。     

非饱和土本构模型在预测地表沉降中的应用

地下水位降低将导致非饱和土体中净平均应力和基质吸力发生变化,从而引起土体体积变化。利用GDS非饱和土三轴仪,研究了重塑非饱和黏土在干燥过程中试样含水率变化和总体积变化。针对某一典型基坑,运用大型有限元软件ABAQUS模拟了井点降水而形成的二维饱和-非饱和稳态渗流。饱和土区域和非饱和土区域沉降变形有着本质的区别,基于试验和数值模拟结果,根据Fredlund非饱和土弹性体变本构模型和饱和土有效应力原理,分别估算了饱和土区和非饱和土区沉降变形量。结果表明,在总沉降量中非饱和土区的沉降量是不可忽略的。随着降水深度的增大,非饱和土区域沉降量在总沉降量中所占比例将逐渐增大。     

Numerical study of soil heterogeneity effects on contaminant transport in unsaturated soil (model development and validation)

The movement of chemicals through soil to groundwater is a major cause of degradation of water resources. In many cases, serious human and stock health implications are associated with this form of pollution. The study of the effects of different factors involved in transport phenomena can provide valuable information to find the best remediation approaches. Numerical models are increasingly being used for predicting or analyzing solute transport processes in soils and groundwater. This article presents the development of a stochastic finite element model for the simulation of contaminant transport through soils with the main focus being on the incorporation of the effects of soil heterogeneity in the model. The governing equations of contaminant transport are presented. The mathematical framework and the numerical implementation of the model are described. The comparison of the results obtained from the developed stochastic model with those obtained from a deterministic method and some experimental results shows that the stochastic model is capable of predicting the transport of solutes in unsaturated soil with higher accuracy than deterministic one. The importance of the consideration of the effects of soil heterogeneity on contaminant fate is highlighted through a sensitivity analysis regarding the variance of saturated hydraulic conductivity as an index of soil heterogeneity. Copyright © 2011 John Wiley & Sons, Ltd.     

基于判别分析法的地震砂土液化预测研究

将距离判别分析方法应用于砂土液化的预测问题中,建立了砂土液化预测的距离判别模型。选用震级、研究深度、震中距、标贯击数、地下水位及地震持续时间等6项指标作为判别因子,以大量的工程实例数据作为学习样本进行训练,建立了线性判别函数对待评样本进行了评价。研究结果表明,距离判别分析模型判别砂土液化效果良好,预测准确度高,回判估计误判率低,可望成为砂土液化预测的有效手段。     

结构性土二元介质本构模型及破损率影响因素的研究

岩土类材料是非均质的具有微缺陷的天然材料,岩土破损力学将岩土材料抽象成由结构块和结构带组成的非均质二元结构体。破损率演化方程反映了岩土材料在外部荷载作用下,结构块向结构带动态变化程度的过程。从细观尺度上分析结构性土体变形、破坏过程发展机制,揭示土体在荷载作用下结构块向结构带转化的过程。最后,针对杨陵黄土,通过三轴剪切试验,研究围压、含水率对破损率演化规律的影响。研究结果表明,单参数模型的破损率与围压和含水率没有关系,使从实验上证明参数模型的破损率是大主应变的函数。     

软粘土的各向异性和小应变条件下的本构模型(ASM)

以能量方程为基础推导出各向异性和小应变条件下弹塑性的本构方程,并将屈服面与修正桥模型和试验测得的屈服面进行了对比,表明新的屈服面比修正剑桥模型更接近试验结果,采用从初始应力到状态边界面的距离为参数的系数来修正硬化准则以模拟土体小应谱条件下的应力－应变特征。将新的本构模型编入有限元程序,对小应变试验进行了计算分析,计算结果表明,新的模型比修正剑桥模型能更好地反映土体的主要力学特性。     

Constitutive modelling of a reinforced soil using hierarchical model

Drained triaxial tests are conducted on natural and reinforced sand under various stress paths. Direct shear tests and pull-out tests are conducted on soil–reinforcement interface and on reinforcement, respectively. The effects of two types of reinforcement, viz, woven and non-woven geotextile and number of layers of reinforcement are investigated. Hierarchical single surface model is used to depict the behaviour of natural and reinforced soil by treating the soil as a single composite material and by considering soil, reinforcement and interface as independent elements. It is shown that the material parameters are very much affected by the type and the number of layers of reinforcement. The hierarchical model provides satisfactory prediction for both natural and reinforced soil. Copyright © 1999 John Wiley & Sons, Ltd.     

膨胀土卸荷蠕变特性及其非线性蠕变模型

鉴于膨胀土滑坡往往表现为长期性、渐进性等与时间相关的特性,利用GDS应力路径三轴仪对膨胀土进行了三轴卸荷蠕变试验。试验结果表明：当偏应力较小时,膨胀土的蠕变曲线仅出现瞬时变形和衰减蠕变;当偏应力达到一定值时,其蠕变曲线也呈现衰减蠕变、稳态蠕变和加速蠕变3个阶段,但其加速蠕变阶段的特征与一般岩土体不同,其蠕变速度近乎常数。膨胀土的应力?应变等时曲线显示,膨胀土卸荷蠕变具有非线性特征,且其非线性程度与蠕变时间和应力水平相关,蠕变时间越长、应力水平越高,其非线性程度越高。基于非线性流变力学理论,提出了一种非线性四元件蠕变模型,将标准线性体与一个非线性黏壶串联,该模型可描述等围压三轴压缩应力状态下膨胀土轴向应变随时间的演变规律。根据膨胀土卸荷蠕变试验结果,采用曲线拟合法对三维非线性模型的参数进行反演识别。拟合曲线和试验曲线对比显示,两者吻合良好,说明该模型可以很好地描述膨胀土的蠕变特性。此外,基于该蠕变模型获取了膨胀土的临界破坏应力,其与常规剪切破坏应力的比值随着固结压力的减小而减小,表明越接近坡面的土层越容易发生蠕变破坏。     

非饱和土弹塑性模型参数的试验确定及有限元法

简要介绍广泛应用的非饱和土Alonso模型的原理与发展,并通过3组共7个试样的非饱和土三轴试验,得到模型的11个参数。在饱和土Cam本构关系的基础上,推导出一个基于非饱和土Alonso模型的应力-应变增量方程的计算公式,该方程与饱和土的本构方程形式相同。编制了能够考虑净应力和吸力二者对土体硬化规律影响的有限元程序,从而为非饱和土弹塑性计算提供一条途径,便于进一步应用于实际工程。     

Generalized trapezoidal numerical integration of an advanced soil model

This paper investigates the numerical performance of the generalized trapezoidal integration rule by using an advanced soil model. The generalized trapezoidal integration rule can include many other integration algorithms by adjusting a single parameter α ranging from 1 to 0. The soil model used is the recently developed middle surface concept (MSC) sand model which simulates different soil response characteristics by using different pseudo‐yield functions. The generalized trapezoidal rule and MSC sand model are used to simulate the responses of sand samples with different relative densities under various initial and loading conditions. Instead of a single step, multiple loading steps bring the sample to the vicinity of failure. These comprehensive investigations examine and compare the influences of various values of α on the numerical solution of integrated constitutive equations, the convergence of Newton's iterative scheme, and the integration accuracy. Copyright © 2007 John Wiley & Sons, Ltd.     

静三轴试验中水泥土力学特性及本构模型研究

通过不排水三轴压缩试验对水泥土的力学特性进行研究,重点分析了围压与水泥掺入比对水泥土体强度、孔压及刚度的影响。试验结果表明,随着围压与水泥掺入比的增加,水泥土体的强度与刚度均增加,而孔压随之减少;随着轴应变的增加,刚度逐渐减少。而在加载初期,刚度衰减较快,随着轴应变的继续增加曲线渐趋平缓,衰减现象并不明显;随着无量纲化偏应力的增加,无量纲化刚度逐渐减少,并且两者为曲线关系。通过对水泥土的刚度软化规律进行描述,得到了修正的邓肯-张模型。通过对试验曲线进行回归分析确定了模型参数。与实测数据对比发现,该修正模型不仅可以较好地模拟水泥土的应力-应变关系,还体现了水泥土应变软化型曲线的特点。