Numerically observed shear bands in soft sensitive clays

The numerical challenges that arise in modelling shear bands in soft sensitive (SS) clays have not yet been fully resolved. Convincing and well-accepted solutions have yet to be found. This paper presents some novel information related to the shear band phenomenon in SS clays. In this study, the hypothesis is that the generation and dissipation of excess pore pressure from shear bands could regularise the strain softening and result in a mesh independent shear band thickness. The generation and dissipation of excess pressure is modelled by a coupled consolidation process. The simulation aims at modelling two counteracting mechanisms in the SS clay. First, the shear band narrows because of strain softening. Second, the internal pore water pressure drainage reduces the rate of strain softening. This counteracting mechanism provides an inherent regularisation technique for SS clays. This study presents some numerical results involving these two counteracting mechanisms. This study also shows that an inherent internal parameter applicable for SS clays can be defined by the ratio between soil permeability and the applied strain rate. In the case of SS clays, the range of this parameter varies from 0 to 0.0002 mm.     

单轴压缩下岩石应变局部化的应变梯度塑性解

通过在屈服函数中引入材料内部长度建立应变梯度塑性增量本构关系,对岩石材料在单轴压缩条件下应变局部化的二维情况进行了研究。首先,建立了在Mises屈服准则下平面应变情况的岩石材料应变局部化带的带宽与其倾角的关系式,讨论了其带宽的分布范围和最大值与最小值形成的条件,给出岩石材料破坏Ⅰ, Ⅱ类变形的条件,并推出破坏的临界角度;其次,对岩石材料Ⅰ类变形进行了详细地分析,讨论了在不同泊松比 和 条件下材料应变局部化带的倾角的变化范围;最后,利用所得的结论同材料破坏的库仑准则进行对比研究,对内摩擦系数给予了合理的解释。     

岩土体渗流自由面问题的重心插值无网格方法

岩土体的渗透破坏、地下工程的防渗设计等无不与渗流计算有关。针对渗流自由面问题,提出一种重心拉格朗日插值的配点型无网格方法。由于渗流自由面问题的求解区域是不规则区域,该方法通过将不规则求解区域嵌入一个正则矩形区域,在正则区域上采用重心拉格朗日插值近似未知函数,利用配点法离散渗流问题的控制方程,将重心拉格朗日插值的微分矩阵离散成代数方程表达的矩阵形式。将自由面上的边界条件通过重心拉格朗日插值离散,通过置换方程法和附加方程法施加边界条件,利用正则区域上的重心插值配点法,通过迭代确定最终自由面的位置。数值算例表明所提出的无网格方法对于求解渗流自由面问题的正确性和高精度。     

概率性分析在黄土场地动力沉降评价中的应用

天然黄土动残余应变研究主要依赖室内动三轴试验,如何将有限室内试验数据合理地应用于场地动力沉降的定量评价,是岩土地震工程领域中研究涉足尚少的重要科学问题。应用概率统计与蒙特卡洛模拟等非确定性分析方法,借助综合考虑固结应力、结构强度、空间体积特性和地震动荷载等动残余应变关键影响参量的估算模型,提出有效弱化土体物性参量不确定性（离散性与随机性）的场地动力沉降概率性评价的思路与方法。应用结果显示,利用该方法得出的场地动力沉降的概率分布特征,能够反映天然黄土动残余应变的基本认知特点,对明晰地震作用下天然黄土场地精细动力沉降特性以及据此采用适当合理的地基处理方法,具有理论意义与实用价值。     

原状膨胀土剪切力学特性的应变速率效应

为研究应变速率对原状膨胀土力学性状的影响,通过GDS三轴试验系统进行了不同速率和围压下的固结不排水三轴剪切试验,分析了应力−应变曲线、孔隙水压力、剪切强度以及破坏模式随应变速率的变化规律。结果表明：不同应变速率下,膨胀土应力−应变曲线均呈应变硬化型。随着应变速率的增加,不排水剪切强度单调递增,引入应变速率参数ρ_0.9后发现,不排水强度增长率为14.3%～23.2%,平均值为18.4%。低围压下,应变速率对孔隙水压力影响较小,随着围压的增大,孔隙水压力的发展趋势由软化型转变为硬化型,孔隙水压力峰值随应变速率的增大而减小。原状膨胀土应变速率效应与其多裂隙性密切相关,破坏形式表现为小应变速率下主剪切带与次剪切带共存,大应变速率下仅有主剪切带,裂隙或多剪切带的出现强化了膨胀土强度的应变速率效应。     

排水条件下卸荷土体变形特性的真三轴试验研究

在真三轴仪上模拟基坑开挖周边土体的应力路径,进行了排水与不排水两种卸荷试验,通过对试验结果的对比分析可知,土体在不同排水条件下应力-应变关系有着显著的差异,说明排水条件对土体的变形性状有一定的影响。进一步研究了排水条件下卸荷土体广义剪应力与广义剪应变的关系,证明土体广义剪应力与广义剪应变关系曲线呈明显的非线性,并根据土体实际应力路径定义了中主应力系数,推导出了土体卸荷剪切模量计算公式。     

基于Hoek-Brown强度准则的应变软化模型在隧道工程中的应用

基于Hoek-Brown强度准则及量化GSI围岩评级系统,分析岩体强度弱化行为,总结出确定岩体峰值及峰后强度参数的方法,提出应变软化模型的简化计算方法,并以FLAC3D数值模拟软件为工具,采用收敛-约束法对深埋大断面隧道支护结构及围岩稳定性进行分析,计算出支护结构安全系数。通过分析得出：随着围岩应力释放,岩体软化参数随之减小,岩体强度弱化行为突出;高地质强度GSI指标（GSI=75）的岩体虽然强度弱化程度较大,但围岩的变形量仍然较小,而低GSI指标（GSI=25）的岩体表现出了理想弹塑性行为;提出的简化Hoek-Brown应变软化模型适合应用于隧道中岩质相近区段围岩的研究;采用收敛-约束法对隧道支护结构的稳定性进行分析,可为隧道工程的初期支护优化设计及安全性评价提供参考。     

基于Hoek-Brown准则的隧洞围岩变形研究

由于在许多实际条件下,比如节理岩体中,线性的M-C准则不太适用,非线性的Hoek-Brown比较适用,因此,可以尝试使用这一非线性屈服准则对洞室变形进行研究。研究隧洞变形时,将围岩分为弹性区、应变软化区、塑性流动区。采用Hoek-Brown准则和非关联流动法则对洞室变形进行了理论推导;软化区域围岩参数随着塑性变形增加而变化,解析法难以求得应力,采用龙格－库塔方法进行数值计算,求解得到塑性软化区和流动区半径,并最终求得洞室变形。通过算例计算表明,在不考虑软化区和流动区时,方法和Carranza-Torres计算结果相差甚小;随着原岩应力的增加,膨胀角对洞室变形的影响增大。     

Thermo-hydro-mechanical modeling of unsaturated soils using isogeometric analysis: Model development and application to strain localization simulation

This study presents a thermo-hydro-mechanical (THM) model of unsaturated soils using isogeometric analysis (IGA). The framework employs Bézier extraction to connect IGA to the conventional finite element analysis (FEA), featuring the current study as one of the first attempts to develop an IGA-FEA framework for solving THM problems in unsaturated soils. IGA offers higher levels of interelement continuity making it an attractive method for solving highly nonlinear problems. The governing equations of linear momentum, mass, and energy balance are coupled based on the averaging procedure within the hybrid mixture theory. The Drucker-Prager yield surface is used to limit the modified effective stress where the model follows small strain, quasi-static loading conditions. Temperature dependency of the surface tension is implemented in the soil-water retention curve. Nonuniform rational B-splines (NURBS) basis functions are used in the standard Galerkin method and weak formulations of the balance equations. Displacement, capillary pressure, gas pressure, and temperature are four independent quantities that are approximated by NURBS in spatial discretization. The framework is used to simulate strain localization in an undrained dense sand subjected to plane strain biaxial compression under different temperatures and displacement velocities. Results show that an increase in the displacement rate leads to reduction in the equivalent plastic strain while an increase in the temperature leads to an increase in the equivalent plastic strain. The findings suggest that the proposed IGA-based framework offers a viable alternative for solving THM problems in unsaturated soils.     

高应力深部洞室模型试验分区破裂现象机制的初步研究

以经典的芬纳公式为手段,以观察到分区破裂现象的模型试验为背景,通过对简化后的圆形巷道围岩应变场和能量场的分析,发现极限平衡区边界存在径向应变的不连续性和弹性变形能的聚集性。认为与开挖洞室呈同心圆的环状断裂的形成,是极限平衡区在边界上由于应变场的不协调而导致的与弹性区的断裂分离。而分区破裂的产生,则是这种断裂分离形成后极限平衡区向深处发展的结果。解释了分区破裂产生时洞周围岩应变和位移随距离洞室增加呈波浪形变化的非单调性规律,给出了分区破裂形成后各破裂区半径的表达式,关系式表明各破裂区之间存在等比关系,比值与地应力状态及围岩力学参数有关。将分区破裂模型试验的参数代入该式,可得,n =1.1左右,与实验完成后剖开模型的各破裂区分布特征基本相符。