冻结黏土卸载状态下双屈服面流变本构关系研究

卸载状态下冻土三轴蠕变试验结果表明,冻土具有明显的剪胀性,且体积应变不能忽略。为了全面反映冻土在卸载状态下的蠕变变形规律,以简单的黏弹塑模型为基本流变元件,提出了冻土服从椭圆-DP双屈服面流变本构力学模型。通过有限元程序的二次开发,将该模型嵌入到ADINA有限元程序中。对冻土三轴蠕变试验过程进行了数值模拟表明,椭圆-DP型双屈服面本构力学模型能够较好地反映冻土的剪胀性和体积应变等变形性质,数值模拟和试验结果吻合良好,可为冻土工程数值计算提供参考。     

Forecast of creep settlements of heavy structures using pressuremeter tests

This paper presents a new method proposed for the prediction of long term settlement of very heavy structures based on a numerical interpretation of long term pressuremeter creep tests. The constitutive model used for the modelling of the soil creep is described. The constants of the constitutive model are determined by fitting a simulated curve, given by a simple Finite Element Method, to the results of a long term pressuremeter creep test. The calculation of long term settlement under nuclear power plants, using the soil parameters defined by this method, is presented and compared with the in situ measured values.     

抛物线型屈服面人工冻土蠕变本构模型研究

利用自行研制的W3Z-200型冻土三轴压缩蠕变仪,对淮南地区深表土冻土进行了三轴蠕变试验。根据试验结果,提出以标准黏弹塑性模型为基础,并耦合温度自由度,应用抛物线型屈服准则对模型中的黏塑性项进行改进,建立了新的本构模型来描述高围压复杂应力状态下冻土的蠕变变形特征。基于ADINA有限元软件平台,开发了冻土本构子程序。对淮南某冻结井筒开挖过程进行了数值模拟,计算结果验证了冻土抛物线型黏弹塑本构模型的正确性和合理性。     

人工冻土黏弹塑性本构参数反分析研究

通过大量的人工冻结黏土蠕变试验研究,获得在高偏应力下人工冻结黏土具有显著的蠕变特性：当应力水平较低时,人工冻土发生黏弹性蠕变变性;当应力水平较高时,人工冻土发生黏塑性蠕变变性。因此,以西原模型为基础,增加一个黏弹性元件项和用抛物线屈服函数代替塑性屈服项,构建了高应力下抛物线型黏弹塑性蠕变本构模型。改进的西原模型可以较好描述人工冻土剪切蠕变特性。基于小生境技术的遗传算法和人工冻土黏弹塑性模型有限元理论,在FEPG（有限元程序生成器）自动生成平台上研发了位移反分析程序。利用现场实测数据对深井冻结井筒开挖过程进行了位移反分析,获得了与试验相一致的冻土本构力学模型参数,数值分析结果与实测数据之间的误差在6%以内,反演结果验证了模型的可行性。冻结井筒开挖过程的位移反分析获得的人工冻土本构模型参数,能反映冻结壁施工过程和有效冻结范围内冻土本构参数平均值,因而具有广泛的代表性,对冻结壁应力场、位移场和稳定性预测具有重要工程意义。     

粘弹塑本构模型及用于冻土数值计算的柔度矩阵

随着寒区经济发展及冻结法施工技术在工程中的推广应用,对冻土蠕变的深入研究显得尤为重要.对冻土单轴流变试验结果分析,用非线性牛顿体替代线性牛顿体,从而对西原模型进行了改进,得到了冻土三向应力状态下的本构方程;并利用D-P屈服准则推导出应用与数值计算的柔度矩阵.通过有限元程序的二次开发,将本构模型添加到大型非线性程序ADINA中.用三轴蠕变试验检验了粘弹塑本构模型,获得相应的试验拟合曲线.结果表明:三轴蠕变试验与理论计算结果吻合良好,说明粘弹塑本构模型可以描述冻土在高应力下反映的非衰减性蠕变变形特征.     

冻土旁压试验与单轴试验对比

为了认识冻土旁压试验结果与常规试验结果之间的关系,在室内分别开展了冻结重塑黏土的旁压试验和单轴压缩试验,并对试验结果进行了对比分析。结果表明,在各级荷载作用下旁压曲线一般都呈现应变速率衰减的趋势,而单轴曲线在冻土破坏时会出现渐进流动阶段。旁压试验的应力-应变曲线呈现应变硬化型,而且出现初始拟弹性阶段;而单轴试验的应力-应变曲线则属于应变软化型,并在轴向应变大约为10%时达到剪应力峰值。温度相同时,旁压试验的剪切强度以及初始弹性模量都要大于单轴试验,且温度越低差值越大。     

Undrained and partially drained behaviour of end bearing piles and bells founded in untreated calcarenite

The subject paper¹ details finite element analyses of underreamed footings installed at the tips of deeply embedded pipe piles. The analysis model incorporated a double-hardening constitutive law and was calibrated using extensive field and laboratory test results. The results of the study were used to judge the relative effects of various aspects of the design, such as partial drainage and layering. Ultimately, this analysis was used to assess the adequacy of the footings as a remedial strengthening measure for a major offshore platform.     

冻粘土蠕变损伤耦合本构关系研究

以冻土试验为基础,提出了冻粘土服从DP屈服准则的粘塑损伤变量,通过推导得到了相关联流动法则下的冻土粘弹塑性损伤耦合本构方程.为了研究粘弹塑性和损伤的耦合作用,将建立的冻土材料粘弹塑性损伤耦合本构模型,通过用户子程序嵌入到ADINA有限元中.结果表明:粘弹塑损伤耦合本构模型能较好的模拟冻粘土加速蠕变过程,与实测蠕变吻合较好.     

高温冻土的蠕变特性试验及蠕变模型研究

为了研究高温冻土蠕变变形特征以及各影响因素对蠕变的作用, 分别在含水量15%、 25%及35%, 荷载100 kPa、 200 kPa及300 kPa, 温度-1.5 ℃、 -0.7 ℃及-0.3 ℃的条件下开展了室内单轴蠕变试验, 分析在无侧限条件下高温冻土在不同温度、 荷载及含冰量条件下的蠕变变形特性。结果表明： 在当前试验条件下, 冻土蠕变变形非常可观, 且蠕变曲线都没有出现渐进流阶段; 温度是影响冻土蠕变的最重要的外在因素, 而含冰量是影响冻土蠕变的关键内在因素; 在高含冰量条件下温度及荷载的改变对蠕变速率的影响非常显著, 甚至引起量级上的差别。在现有试验条件下, 高温冻土蠕变过程可利用Burgers黏弹性模型来较好地描述。     

Analysis of liquefaction susceptibility of nearly saturated sands

A new constitutive formulation for simulating the behaviour of nearly saturated sands under seismic loads is presented. The formulation is based on combining the Henry's law for dissolution of gas in water, the ideal or perfect gas law and the law of conservation of mass. The effects of transient air dissolution in water on the compressibility of partially saturated soils are also taken into account. The model was calibrated based on numerical simulations of isotropically consolidated cyclic triaxial tests conducted on partially saturated samples of Toyoura sand. A multi‐yield plasticity soil constitutive model implemented in the finite element code DYNAFLOW was used for these numerical simulations. It is shown that the formulation proposed here is able to reasonably predict the soil cyclic undrained behaviour at various degrees of saturation (95% and higher). Copyright © 2006 John Wiley & Sons, Ltd.     

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.     

青藏直流输电线路冻土长期抗剪强度预测及影响因素分析

冻土是指0℃以下,含有冰的各种岩石和土壤,由于其复杂的物理力学特性,给青藏直流联网工程的施工和设计带来了很多难题。工程沿线冻土的长期抗剪强度,更是由于土质、含冰量以及密实度分布的复杂特性,影响了其在使用寿命内的安全运营和稳定性。为解决这一工程难题,室内制备重塑试样并在-2℃的温度下对3种土质类型（粉质黏土、粉砂和细砂）、3种含冰量（饱冰、富冰和多冰）的密实和松散冻土进行直剪快剪的蠕变试验,得到了不同土质、含冰量以及密实度类型的冻土剪切蠕变特性,并分析和预测了长期抗剪强度的变化规律和影响因素,结果表明:（1）短时间内冻土承受荷载的能力随含冰量的增加而增大;随着加载时间的延长,冰的流变特性表现出来,使冻土的蠕变变形加大,冻土的强度降低;（2）密实度越大、含冰量越大的冻土,初期强度越高;随着时间的增加,强度开始衰减,且含冰量越大,密实度越小的冻土衰减速度越快;（3）长期强度的变化规律不受土质的影响,说明在工程施工阶段,若单纯考虑长期强度时,只要不同土质冻土保持冻结状态,可同一对待和处理。以上结论为输电线路冻土基础设计、施工和安全运营提供了一定的理论基础,也为在青藏地区开展的其他工程活动提供了可借鉴的试验数据和资料。     

广州某地铁人工冻结法施工热力分析

以广州某地铁采用人工冻结法施工为例,考虑相变和衬砌水泥水化热的生成,根据带相变瞬态温度场的热量控制微分方程,采用伽辽金法推导出的有限元计算公式计算分析了此地铁人工冻结施工过程中的温度场分布规律.并根据温度场的分布特点,确定了开挖时的人工冻结壁厚度.最后,对冻结壁和其周围土体分别应用不同的力学本构模型进行施工过程的力学模拟.结果表明,采用人工冻结法作为一种辅助的施工方法,不但能满足施工过程中结构的强度要求,而且能满足环境沉降量的要求.     

利用旁压试验的砂土初始状态反演分析

讨论了砂土旁压试验的反分析问题。建议利用旁压试验实测数据反演确定土体的初始状态而非材料参数,材料参数可以通过实验室常规试验确定。主要利用了新型的亚塑性本构模型参数不受土体状态影响的特点。主要内容包括：建立有限元数值模型,利用文献中的试验数据确定了特定砂土的材料参数;通过模拟砂土的实验室旁压试验,验证了数值模型;通过改变土体初始应力状态和相对密实度进行数值计算,并利用数值计算成果,建立了可适用于砂土初始状态反演分析的实测旁压荷载和土体中应力与相对密度的数学关系。     

Computer-aided calibration of a soil plasticity model

A non-linear optimization technique based on the quasi-Newton approach is employed to back-calculate certain model parameters of a simple, bounding surface, soil plasticity model from in situ pressuremeter data. The theoretical response corresponding to a given set of parameters is generated by finite element analysis. A semi-analytical procedure is developed for the accurate and efficient evaluation of the gradient of objective function with respect to the model parameters of interest. The BFGS update is used to update the Hessian. Results of a series of numerical experimentation using artificial pressuremeter responses is first reported and discussed. A set of laboratory cavity expansion data is then used to calibrate the constitutive model.     

考虑温度影响的冻土蠕变本构模型

冰作为冻土的基本组成部分,对冻土蠕变的加速蠕变阶段有重要影响。温度通过影响冻土中冰的冻结与融化过程,及其黏塑性流动,引起冻土结构的强化与弱化,从而成为了决定冻土蠕变力学行为的关键因素之一。同时,外部应力也造成冻土的强化与弱化,影响着冻土的蠕变。通过引入硬化因子与损伤因子来考虑温度、应力造成的冻土材料强化与弱化。硬化因子H代表了蠕变过程中强化效应的大小,而损伤因子D则代表了由弱化效应造成的冻土材料相关参数的折减比例,进而提出了适用于冻土的改进西原蠕变本构模型。该模型预测值与试验数据的比较表明：改进的模型不仅能较好地描述初始蠕变阶段、稳定蠕变阶段,而且相比于传统模型,能更好地描述冻土加速蠕变阶段,具有合理性与一定的实用性。     

深部人工冻土黏弹塑损伤耦合本构模型研究及其验证

采用改进的Mohr-Coulomb屈服准则(Zienkiewicz-Pande准则), 建立了能描述人工冻土3个阶段蠕变特性的黏弹塑损伤耦合本构模型. 基于该本构模型建立了模型的柔度矩阵, 将本构模型嵌入到ADINA有限元程序中. 计算表明: 数值模拟得到的人工冻土蠕变特性曲线与实测值吻合良好, 最大误差为4.8%, 说明基于改进的Mohr-Coulomb准则建立的黏弹塑损伤耦合本构模型能够准确的描述人工冻土剪切特性和3个阶段的蠕变特性.获得的本构模型具有参数少、容易确定等优点. 深部人工冻土黏弹塑损伤耦合本构模型的建立, 对于人工冻土结构物长期稳定分析和预测具有重大的理论指导意义.     

深厚冲积层人工冻土蠕变参数的模糊遗传反演

冻土的蠕变特性对深井冻结法施工至关重要。针对某矿区人工冻土在-5℃、-10℃、-15℃和-20℃下进行单轴抗压强度试验,发现冻土的抗压强度受冻结温度变化影响,两者间为线性反比例关系。通过小生境原理对传统的遗传算法作模糊随机改进,给出模糊遗传算法的步骤思路,进而运用该算法反演冻土蠕变模型中的参数值,获得各温度下的蠕变模型。试验结果表明：蠕变模型计算值在蠕变各阶段与试验值吻合较好,准确反映了冻土蠕变的整体规律。可见,模糊遗传算法能有效反演蠕变参数,较传统方法更符合工程实际。     

Single‐and multi‐objective genetic algorithm optimization for identifying soil parameters

This paper discusses the quality of the procedure employed in identifying soil parameters by inverse analysis. This procedure includes a FEM‐simulation for which two constitutive models—a linear elastic perfectly plastic Mohr–Coulomb model and a strain‐hardening elasto‐plastic model—are successively considered. Two kinds of optimization algorithms have been used: a deterministic simplex method and a stochastic genetic method. The soil data come from the results of two pressuremeter tests, complemented by triaxial and resonant column testing. First, the inverse analysis has been performed separately on each pressuremeter test. The genetic method presents the advantage of providing a collection of satisfactory solutions, among which a geotechnical engineer has to choose the optimal one based on his scientific background and/or additional analyses based on further experimental test results. This advantage is enhanced when all the constitutive parameters sensitive to the considered problem have to be identified without restrictions in the search space. Second, the experimental values of the two pressuremeter tests have been processed simultaneously, so that the inverse analysis becomes a multi‐objective optimization problem. The genetic method allows the user to choose the most suitable parameter set according to the Pareto frontier and to guarantee the coherence between the tests. The sets of optimized parameters obtained from inverse analyses are then used to calculate the response of a spread footing, which is part of a predictive benchmark. The numerical results with respect to both the constitutive models and the inverse analysis procedure are discussed. Copyright © 2011 John Wiley & Sons, Ltd.     

Undrained Shear Strength of Normally Consolidated and Overconsolidated Clays from Pressuremeter Tests: A Case Study

Undrained shear strength (s_u) of foundation soil of Marquette interchange near Milwaukee, Wisconsin was evaluated from the results of a number of pressuremeter tests conducted on normally consolidated (NC) organic silts and overconsolidated (OC) silty clay. The s_u-values were interpreted from traditional closed-form methods. The pressuremeter geometry and test sequence as well as response of the soil profiles were also simulated using axisymmetric finite element (FE) method with Cam-Clay soil model. The Cam-Clay model parameters were estimated from laboratory tests on undisturbed soil samples. Results show that the s_u estimated from the rate of cavity pressure change with volumetric strain (referred to as direct traditional method) is almost twice the s_u estimated from an indirect traditional method that estimates s_u from shear modulus, in situ horizontal stress, and ultimate cavity pressure obtained from the cavity pressure curves. The s_u-values predicted from the FE models are lower than those estimated from the traditional methods and shows that the assumption of infinite pressuremeter length in traditional methods results in overprediction of undrained shear strength by a factor of 1.5 for NC clay and 2.2 for OC clay. The results of finite element analysis considering Cam-Clay soil model and finite length for pressuremeters suggest the undrained shear strength of 63 ± 7 kPa for NC organic silt and 259 ± 68 kPa for OC silty clay.