UH模型在有限元分析中的应用

由姚仰平等提出的超固结土统一硬化模型(UH模型)能够描述超固结土的硬化、软化、剪缩和剪胀等应力应变特性以及应力路径对它的影响。通过改进,直线型的Hvorslev线由抛物线型的Hvorslev线代替,UH模型的参数与剑桥模型(Cam-clay模型)的参数相同。应用该模型,对三轴试验进行有限元模拟,分析了与应用Cam-clay模型所得应力应变关系曲线的异同点;针对刚性承压板试验,进行了地基承载力的三维有限元分析,比较了与应用Cam-clay模型和摩尔库仑模型(Mohr-Coulomb模型)等所得p-s曲线的不同,分析了地基中不同点处的应力-应变关系曲线,并与Boussinesq弹性理论计算的应力做了比较;进行了地基固结沉降数值分析,比较了孔压、沉降的变化规律,体现了UH模型的优越性。     

考虑渐近状态特性的超固结土本构模型

基于超固结土的UH模型,结合渐近状态概念,将适用于饱和砂土的渐近状态本构模型推广为考虑渐近状态特性的超固结土本构模型,用于描述土的超固结性和渐近状态特性。通过采用变换应力方法实现了模型的三维化。新提出的本构模型简单、参数较少,与剑桥模型相比,仅增加了一个材料参数,即伏斯列夫面斜率。该模型能够合理反映超固结土的硬化、软化、临界状态、剪缩、剪胀、应力路径依赖性、渐近状态等特性。     

K0固结饱和土柱孔扩张问题弹塑性分析

采用K0固结各向异性土体本构模型,将柱孔扩张后周围土体分为弹性区和塑性区,根据柱孔扩张理论和边界条件,推导出K0固结状态下饱和天然土体柱孔扩张问题弹塑性区的应力、塑性区半径以及超孔隙水压力的理论解答。同时,通过算例与修正剑桥模型解答进行对比分析,结果表明,土的不同初始固结状态对柱孔扩张后孔周围的应力和超孔隙水压力产生很大影响,采用考虑K0固结诱发各向异性土体本构模型所得到的应力和超孔隙水压力解答大于修正剑桥模型的解答,但塑性区影响半径却明显小于后者。     

基于SMP准则柱孔扩张问题相似解

以柱孔扩张理论为基础,将柱孔扩张后塑性区半径b作为时间刻度,引入扩孔速度概念,将塑性区内任意一点的应力、应变描述成塑性区半径b和柱孔扩张后扩张体单元距孔心距离r的函数,采用平面应变条件下无黏性土的SMP准则和不相关联流动法则,考虑到无黏性土的剪胀性,推导出柱孔扩张后极限扩张力 的相似解。通过参数分析表明,土体强度和剪胀性对塑性区半径和极限扩张力产生很大的影响,即随着土体强度和剪胀性的增加,极限扩张力明显增加,相应塑性区半径却随之减小。通过与基于摩尔-库仑准则解答进行对比分析表明,考虑中间主应力影响的SMP准则所得到的极限扩张力明显高于基于摩尔-库仑准则的解答。该解答为无黏性土中进行静力触探（CPT）试验和静压桩等岩土工程问题提供必要的理论依据。     

统一硬化模型与下加载面模型的理论关系

Hashiguchi下加载面本构模型是国际上具有较大影响力的超固结土本构模型之一。因此,选取该模型进行理论和构建方法分析,并与UH模型进行对比。对比结果表明,Hashiguchi下加载面模型通过定义单一的数学公式来拟合超固结土内变量的变化规律,从而得到了超固结土应力−应变关系的公式;而统一硬化（unified hardening,简称UH）模型掌握了超固结土体的强度特性,因而能够更准确、合理地描述超固结土的应力−应变关系。相较于Hashiguchi下加载面本构模型,UH模型在理论上更为先进。同时,两个模型的试验预测与试验数据的对比也验证了UH模型在数值计算上的正确性与优越性。     

特殊应力路径下岸坡饱和土体变形特性模拟

库岸边坡土体受到库水位涨落的周期性影响,其受力过程可描述为在前期固结压力范围内有效应力的多次循环加、卸载过程。通过室内应力路径三轴试验模拟土体孔隙水压力的循环加、卸载,并通过对统一硬化模型的扩展和加、卸载准则的修正,来预测土体的变形发展规律,研究结果表明:循环加、卸载过程中,土体的体积变形呈周期性弹性变化,剪应变呈螺旋上升趋势且增幅随循环次数逐渐减小;卸围压屈服阶段,土体变形呈现体积膨胀的剪胀特征和剪应变增幅加快的脆变趋势,且土体剪切变形在6%左右时即出现实测孔压骤减和峰值强度点;建立压缩曲线系数与加、卸载次数的关系,并修正超固结土的加、卸载准则,扩展后的统一硬化模型能够较好地预测特殊应力路径下岸坡饱和土体的变形特性。     

不同应力路径下超固结黏土试样变形局部化分析

基于改进伏斯列夫面超固结黏土三维本构模型,利用有限元软件ABAQUS材料子程序接口,采用回映应力更新算法,实现了该模型在有限元分析中的应用。通过该模型与比奥固结理论的耦合,对超固结比为8的超固结黏土在三轴压缩、三轴伸长及平面应变应力条件下的变形局部化问题,进行了水-土耦合弹塑性有限元分析。分析结果表明：剪切带带内、带外点经历不同应力路径;剪切带带外单元经历了体缩、剪胀及被吸水体缩过程,而剪切带带内单元一直保持剪胀趋势;剪切带的形成伴随着剪胀,剪切带内、外出现了负的孔压,且孔压的分布也具有局部化特性。关于剪切带带内、外的孔隙水压及体变变化趋势与剪切速率有关,而平面应变介于三轴压缩与三轴伸长之间,但平面应变较早出现剪切带。孔隙水的迁移速度影响剪切带带内单元的剪胀,进而影响剪切带的形成及发展;而围压和弱单元位置也对剪切带的形成也有影响。     

考虑Hansbo渗流的饱和黏土一维弹黏塑性固结分析

为进一步深入探讨饱和黏土的固结机制,引入考虑时间效应的统一硬化(UH)本构模型描述饱和黏土变形的弹黏塑性,以及Hansbo渗流方程描述固结过程中的非Darcy渗流,修正了饱和黏土一维固结方程,并利用有限差分法进行了求解。通过与文献中一维流变固结试验结果和理论计算结果的对比,验证了所提计算方法和UH模型的适用性。然后探讨了Hansbo渗流参数及UH模型参数对流变固结过程的影响。计算结果表明:土体的黏滞效应是引起流变固结初期地基不排水面附近出现孔压升高现象的主要原因,而渗流的非Darcy特性增强了这一效应,同时土体的这两个特性还延缓了固结中后期饱和黏土地基中孔压的整体消散和地基的沉降。另外,随着回弹指数和超固结参数的增大,或渗透指数的减小,上述流变固结特性更加明显。但是Hansbo渗流参数和渗透指数对地基最终沉降量均无影响。     

饱和土体小孔扩张问题的弹塑性解析解

将土体在圆孔扩张过程中的应力分布分为三个区域，基于修正剑桥模型，推导了圆孔扩张过程中土体在三个区域的超孔隙水压力的解析表达式，对比分析了固结比对土体中应力分布的影响。分析推导可为沉桩、旁压试验、静力触探等岩土工程问题提供理论依据。     

用快剪试验推求先期固结压力

土体先期固结压力是判断土体应力历史的关键性指标。依据Wroth在一维固结条件下提出的不排水强度与超固结比（OCR）的关系式,利用快速直剪的剪切强度计算临界状态孔隙水压力系数Λ₀,从而计算出OCR以及先期固结压力。对南京某黏土的重塑饱和黏土的试验分析结果表明,快剪试验确定的超固结比与设计超固结比基本一致,因此可作为求解土体先期固结压力的一个途径。      

胶结结构性土统一硬化模型

结构性土颗粒间的胶结使试样剪切破坏最终应力比高于相应重塑土,也限制了试样剪切时体积应变的自由发挥。在考虑结构垮塌为主的结构性土统一硬化(UH)模型基础上,将应力空间中静止的临界状态线扩展为动态的移动临界状态线。据此,通过建立新的屈服面方程并修正剪胀方程,将结构性土统一硬化(UH)模型扩展为胶结结构性土统一硬化(UH)模型。相对于原模型,新模型增加了1个模型参数,即初始胶结应力,反映土颗粒之间的初始胶结作用。通过4种结构性土试验数据与模型预测对照表明：所提模型能够较合理地描述结构性土等向压缩、常规三轴排水与不排水剪切等特性。     

An analytical study of the effect of penetration rate on piezocone tests in clay

Finite element cavity expansion analysis investigating the effect of penetration rate on piezocone tests in clay is presented. A coupled analysis was performed, in which the rate of cavity expansion was linked to the penetration rate of the cone and the cone angle, using the assumption that the deformation was wholly radial, and took place only between the cone tip and the cone shoulder. The soil was modelled using modified cam clay with two sets of parameters and varying values of overconsolidation ratio (OCR). The influence of penetration rate on the stress and pore pressure distributions was examined. For slower penetration rates, the excess pore pressure at the cone shoulder is lower since consolidation is permitted coincident with penetration. The radial profiles of post‐penetration voids ratio demonstrate that partially drained penetration is permitted by volume change in the near field, in addition to radial movement in the far field. The radial distribution of excess pore pressure after slow penetration differs from the undrained case, with a relatively low radial gradient existing at the cone face. As a result, the dissipation curves after slow penetration lag behind those following fast penetration. The cone velocity is made dimensionless by normalizing with the coefficient of consolidation and the cone diameter. ‘Backbone’ curves of normalized velocity against normalized tip resistance and excess pore pressure capturing the transition from undrained to drained penetration are derived. The normalized pore pressure backbone curve is unique, whilst the normalized tip resistance shows a small dependency on OCR. These backbone penetration curves are compared with centrifuge model piezocone tests conducted at varying rates, and subsequent dissipation tests. The numerical and experimental results suggest that the value of consolidation coefficient operative during the dissipation phase is 2–4 times higher than the virgin compression value due to changes in the operative soil stiffness, as demonstrated from the stress paths of individual soil elements. The use of multi‐rate penetration tests to deduce values of consolidation coefficient is discussed, in light of these differences. Copyright © 2005 John Wiley & Sons, Ltd.     

Similarity solution for cavity expansion in thermoplastic soil

This paper presents an analytical solution for cavity expansion in thermoplastic soil considering non‐isothermal conditions. The constitutive relationship of thermoplasticity is described by Laloui's advanced and unified constitutive model for environmental geomechanical thermal effect (ACMEG‐T), which is based on multi‐mechanism plasticity and bounding surface theory. The problem is formulated by incorporating ACMEG‐T into the theoretical framework of cavity expansion, yielding a series of partial differential equations (PDEs). Subsequently, the PDEs are transformed into a system of first‐order ordinary differential equations (ODEs) using a similarity solution technique. Solutions to the response parameters of cavity expansion (stress, excess pore pressure, and displacement) can then be obtained by solving the ODEs numerically using mathematical software. The results suggest that soil temperature has a significant influence on the pressure‐expansion relationships and distributions of stress and excess pore pressure around the cavity wall. The proposed solution quantifies the influence of temperature on cavity expansion for the first time and provides a theoretical framework for predicting thermoplastic soil behavior around the cavity wall. The solution found in this paper can be used as a theoretical tool that can potentially be employed in geotechnical engineering problems, such as thermal cone penetration tests, and nuclear waste disposal problems.     

Analytical solution of undrained cylindrical cavity expansion in saturated soil under anisotropic initial stress

This paper presents an analytical solution for undrained elasto-plastic cylindrical cavity expansion in saturated soil under anisotropic initial stress. The problem is formulated by assuming small-strain deformation in the elastic zone and large-strain deformation in the plastic zone. Plastic yielding is determined by the Tresca failure criterion and an associated flow rule. Two stress functions are used to describe the stress state in the two zones around the cavity. The elasto-plastic boundary can subsequently be determined by solving the two stress functions under the stress boundary conditions. Additionally, the cavity pressure-expansion relationship, the total stress and the excess pore pressure around the cavity wall under anisotropic initial stress can be easily obtained by application of a numerical integration. The results show that the cavity pressure and excess pore pressure under the isotropic initial total stress (K = 1) are larger than those under the anisotropic initial stress (K > 1 and K < 1), which is true at all states of the expansion. The higher value of K develops lower stress and pore pressure around the cavity wall at the ultimate states. However, the stress and excess pore pressure are not sensitive to the value of K. The present solution may be used for analyzing the uplift capacity of plate anchors in soils and Horizontal Directional Drilling (HDD) problems such as the tunneling, and pipeline installation.     

XCC pile installation effect in soft soil ground: A simplified analytical model

This paper presents a simplified analytical model: the X-shaped cylindrical cavity expansion method (XCEM), which can be used to interpret and predict the displacement, stress, and excess pore pressure caused by the X-section cast in place concrete (XCC) pile installation in soft soil. Analytical solutions for the displacement and strain field are obtained with a streamlined solution, base on the strain path method (SPM). The stress and excess pore pressure can then be subsequently determined with the assumption of the elastic-perfectly plastic soil behavior. A positive agreement was found when the theoretical prediction of the displacement, stress, and excess pores pressure was compared against the field measurements. The proposed XCEM improves the conventional cavity expansion method (CEM) and offers a framework for understanding the non-circular cross-section penetrator problem, which is different from the conventional circular penetrator problem.     

高低压下不同力学特性的砂土统一模型

为构建能够反映砂土高低压下不同力学特性的统一模型,分析了砂土在较大的压力范围内的力学试验结果并获取其强度、等向压缩以及临界状态特性参数。通过引入应力路径相关因子来修正塑性应变增量中与应力路径相关的部分,从而使得模型硬化参量能够反映密实砂土在常压下的剪胀特性。同时,基于砂土的临界状态特性提出能够与砂土内部状态相对应的潜在状态面概念,由屈服面与潜在状态面间的动态关系确定加载过程中的动态密实参数及潜在强度,进而使得硬化参量也能够反映密实砂土在常压下的软化特性及高压下的剪缩、硬化特性。分析模型屈服面及潜在状态面间的演化规律并对不同压力等级下的砂土受荷力学行为进行模拟预测,证实了该模型能够反映密实砂土常压下剪胀软化及高压下剪缩硬化的特性。