应变局部化弹黏塑性分析的内嵌局部软化带模型

针对岩土介质结构在破坏过程中局部化变形的问题,结合位移不连续的思想,提出内嵌局部软化带模型来捕捉结构中的局部化带。通过虚功原理建立了含局部化带影响的弹黏塑性的有限元计算模式,其中分叉理论作为局部化判断条件。模型将局部化带的形成视为一个黏塑性屈服流动过程,从而能够连续地描述局部化变形前后的力学性质。特点是计算量小、物理意义明确,可以方便地整合到传统有限元分析程序中。算例表明,计算模型是合理和有效的。     

路堤侧向变形性状及预测模型研究

结合大量实测资料和已有的研究成果,分析了路堤荷载作用下水平向应变随深度的变化规律,得出了侧向变形沿深度呈“弓”型分布的变化规律;提出了符合该变化规律的预测模型及参数求解方法,模型中a、b和c等3个待定参数能简单地求出。工程实践表明,该预测模型不仅反映了侧向变形随深度的变化规律,且有一定的精度,能反映出最大侧向变形发生的位置,应用效果良好。     

Crust restoration for the western Lachlan Orogen using the strain-reversal,area-balancing technique: implications for crustal components and original thicknesses

Strain reversal of structural/stratigraphic profiles at different scales in the western Lachlan Orogen provides a perspective on original crustal thickness estimates, the former depositional basin width of the proto-western Lachlan Orogen, the original sedimentary-fan thickness, and the possible length extent of lower crust lost by subduction. Retrodeformation using strain-reversal techniques allows basin reconstruction giving an original width of the western Lachlan Orogen basin receptor of between 800 km (minimum) and ～1150 km (maximum), depending on the amount of stratal duplication allowed in the turbidites. Crude area balancing of the regional cross-section, adding in sectional volume lost by erosion and assuming strain compatibility between the upper and lower crust, suggests that the predeformation crustal thickness ranges between 15 km and ～21 km, with a lower crustal thickness (oceanic lithosphere) of ～9 km and a turbidite fan thickness of ～6 km (minimum) and ～12 km (maximum allowable), respectively. Disparity between the calculated fan thickness and that derived from measured stratigraphic sections adjusted for strain (～6 km) indicates that some form of crustal stacking must be important in structural thickening of the turbidite crustal component. By varying shortening due to fault stacking, mass balance dictates the mismatch of the upper crustal (uc) and lower crustal (lc) retrodeformed lengths, and therefore provides an estimate of lower crustal loss by subduction. End members range from: (i) a 12 km-thick fan without fault duplication, a basin width of ～800 km where uc = lc giving no lower crustal loss by subduction; to (ii) a ～6 km fan, requiring duplication by faulting, a basin of ～1150 km where uc > lc, and ～360 km of lower crust length (～30%) lost by subduction. This suggests that the total thickness of underplated igneous material in the western Lachlan Orogen is low, probably < ～2 km.     

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应变软化模型适合应用于隧道中岩质相近区段围岩的研究;采用收敛-约束法对隧道支护结构的稳定性进行分析,可为隧道工程的初期支护优化设计及安全性评价提供参考。     

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.