基于残余强度修正的岩石损伤软化统计本构模型研究

岩石随着围压的增大,残余强度的增加幅度比峰值强度大,残余强度逐渐成为影响岩石全应力-应变曲线峰后段的主要因素,因此,建立岩石损伤软化统计本构模型时,对残余强度进行修正是非常必要的。基于岩石应变强度理论以及岩石微元强度服从Weibull随机分布的假设,考虑岩石峰后残余强度对损伤变量进行修正,在微元破坏符合Hoek-Brown屈服准则条件下,建立了能够反映岩石峰后软化特征的三维损伤统计本构模型;依据岩石试验全应力-应变曲线的几何特征,推导出本构模型参数的数学表达式,并基于花岗岩室内试验数据对模型参数进行探讨,研究了Weibull分布参数与本构模型的关系,探讨了损伤修正系数的取值和岩石累积损伤的扩展过程;将建立的本构模型理论曲线与4种岩石（斑状二长花岗岩、细晶大理岩、砂岩和粉砂质泥岩）不同围压下的常规三轴压缩试验曲线进行对比分析,结果表明,该模型能很好地描述岩石破裂过程的全应力-应变关系和表征岩石残余强度特征。这对于岩石损伤软化问题以及岩石加固处理措施的研究均具有重要的意义。     

岩石全应力-应变过程关键特征点损伤特征分析

基于Weibull分布统计损伤力学理论,针对传统本构模型对初始状态损伤和峰后特征表征不足的局限进行了优化,建立能够表征初始损伤对力学行为响应的本构方程,并基于优化的本构方程确定岩石全应力-应变过程中关键特征点的力学特性和损伤状态关系,将岩石受荷损伤劣化过程划分为6个损伤阶段,阐释了岩石全应力-应变过程中的损伤演化机制,同时结合RFPA-2D数值仿真分析确定岩体在峰后曲线上的结构失稳点应变是峰值应变的(1.00～1.32)倍;损伤阈值约为0.4～0.6,损伤起始强度约为峰值强度15%;结构面贯通破坏强度约为峰值强度的41%。从理论上确定了残余强度点应变约为峰值应变的4倍。     

砂岩单轴压缩与干湿循环耦合损伤试验研究

设计并进行了不同干湿循环次数下的砂岩单轴压缩试验,得到试验应力-应变曲线。观察对比各组曲线的各阶段特征,分析干湿循环作用对砂岩的变形破坏和力学性能的影响规律,发现干湿循环交替过程导致砂岩试件单轴抗压强度和弹性模量均显著降低并逐渐趋于一个定值,破坏特征则由脆性向脆—延性转变,砂岩受到了不可逆的渐进损伤作用。建立了单轴压缩与干湿循环耦合作用的随机损伤力学模型,推导出随机损伤力学本构关系表达式,拟合出理论曲线并与试验曲线进行对比分析,证明该模型能较好地反映干湿循环作用下砂岩的力学特性曲线。     

不同含水率与围压下伊犁高温冻土三轴力学试验特性研究

为了探究含水率与围压变化对高温冻土物理力学性质的影响，以新疆伊犁河谷高温冻结黄土为研究对象，开展了黄土的矿物成分、物理性质，以及不同含水率和围压条件下冻土的三轴压缩试验。结果表明：伊犁黄土的粉粒与黏粒粒组含量占比较高，对冻融作用的反应敏感。低含水率时表现为应变软化现象，破坏形态以脆性剪切破坏为主，饱和含水率时表现为应变硬化现象，破坏形态以塑性鼓胀变形破坏为主，软化系数随含水率增大而逐渐减小。随着含水率增大，峰残内摩擦角逐渐降低，峰残黏聚力逐渐增大，变形模量逐渐增大。随着围压增大，弹性模量和损伤演化特征参数均逐渐降低，引入的损伤力学本构模型能够较好地描述高温冻土在不同含水率和围压影响下的应力应变全过程。研究成果可为伊犁河谷冻融滑坡成灾机理研究提供力学参数与理论依据支撑。     

基于断裂及高温损伤的岩石蠕变模型研究

为了反映热-力耦合作用下岩石蠕变变形的全过程,依据断裂力学原理提出了岩石裂纹扩展的临界损伤应力和一个新的可描述岩石在稳态蠕变阶段与临界损伤应力相关的非线性黏性分量,在传统西原模型和Burgers模型的基础上,将指数形式的损伤变量、临界损伤应力以及与其有关的非线性黏性分量引入到流变微分方程,通过叠加原理推导了考虑温度效应的单轴和三轴压缩条件下岩石的流变本构关系,建立了岩石的热-力耦合损伤蠕变本构模型。利用不同温度、不同应力条件下花岗岩的三轴蠕变试验曲线和本文蠕变模型的计算曲线进行比较,结果表明本文蠕变模型能够较好地模拟岩石在初始瞬态、稳态和加速蠕变阶段全过程的变形规律,验证了所建模型的有效性和合理性。该模型为分析高温、高应力环境下岩石工程的长期变形和稳定情况提供了理论依据。     

考虑损伤阀值影响的钢纤维混凝土损伤本构模型研究

从室内试验得到的不同体积率下的钢纤维混凝土全应力-应变曲线的特性,探讨了损伤变量的变化规律及建立钢纤维混凝土损伤本构模型时考虑损伤阀值影响的必要性。基于连续损伤力学理论和统计强度理论,在传统的两参数Weibull分布函数基础上引入损伤阀值参数（位置参数）,建立可考虑损伤阀值影响的钢纤维混凝土损伤本构模型,该模型可以反映钢纤维混凝土在低应力水平或变形较小时的线弹性变形特性;根据经验确定损伤阀值参数后,通过混凝土应力-应变全曲线的几何条件求解Weibull分布函数另外两个参数的表达式。最后,基于MATLAB的分析计算结果,并通过与钢纤维混凝土单轴压缩试验实测结果对比,证明模型可以很好的反映单轴受压状态下钢纤维混凝土的应力-应变关系,探讨了损伤阀值取值的大小。     

陇东Q3结构性黄土压剪损伤本构模型试验研究

Q3黄土在我国西北地区分布广泛,且大部分都具有显著的结构性。随着国家“一带一路”倡议的不断向前推进,作为丝路沿线的西北黄土分布地区将迎来新的建设大潮。黄土的结构变形特性非常复杂,深入研究黄土在压缩、剪切条件下的结构损伤变形特性,并依此构建黄土的本构关系在理论研究及现场工程应用中具有非常重要的意义。通过对黄土进行均等压缩试验及三轴剪切试验,基于损伤力学思想,提出黄土结构在均等压缩条件下的平均正应力损伤比,在剪切条件下的平均正应力损伤比及偏应力损伤比。根据弹性、塑性应变确定塑性势线,进而确定其屈服函数;将确定的黄土结构损伤比引入到屈服函数中,得到一定结构损伤时黄土的屈服函数表达式;验证了选取塑性体应变作为本构模型硬化参量的合理性;根据硬化参量与相关试验参数的联系,推导出结构性黄土在压剪条件下的损伤本构模型。经过实测应力-应变曲线与本构模型推算得到的应力-应变曲线对比可知,所建立的本构模型可以较好地反映黄土在压剪条件下结构损伤演化变形过程,具有较好的工程应用前景。     

不同应力路径下剪切带的数值模拟

采用回映应力更新算法,编写了基于伏斯列夫面的超固结黏土本构关系模型子程序,嵌入非线性有限元软件ABAQUS。通过对单元试验进行三轴压缩、三轴伸长及平面应变等问题的模型预测,再现了超固结黏土在不同初始超固结比和应力路径时的变形和强度特性,从而验证了子程序的正确性。借助该本构模型,对三轴压缩、三轴伸长及平面应变应力路径下超固结黏土体变形局部化问题,进行了三维数值模拟。分析结果表明：超固结黏土在三轴压缩及伸长状态时,土体变形局部化在应力-应变关系软化时出现,而平面应变状态时,在应力-应变关系硬化阶段出现,其超固结黏土的剪胀特性在剪切带的形成过程中起重要作用。     

季节冻土区路基土的损伤本构模型研究北大核心CSCD

为探索不同冻融循环次数和围压作用对季节冻土区路基土在不同含水率条件下的应力-应变关系及损伤机理,以辽宁阜新市某公路区间为试验路段,采用环刀法选取试验路段路基土样,通过冻融循环及三轴压缩试验,得到在不同冻融循环次数及围压作用下路基土的应力-应变曲线变化规律;根据损伤力学及统计学原理,将Weibull分布与Lemaitre有效应力原理相结合,建立了季节冻土区路基土损伤本构模型。结果表明:冻融循环作用可导致路基土试件强度降低,且随着试件含水率的增加,应力-应变曲线呈现由应变软化逐渐表现为应变硬化特征,试件最优含水率为应变软化与应变硬化的分界点;随着围压的增大试件强度增加,在冻融循环次数较少且围压较低时,试件的应力-应变曲线出现峰值应力,曲线表现出应变软化特征,在冻融循环次数较多且围压较大时容易出现应变硬化现象;经对比分析,所建立的损伤本构模型与试验应力-应变曲线吻合度较好,且模型所需参数均可通过三轴试验获得,说明该模型能够较好地描述季节冻土区路基土的应力-应变关系,具有实用性。另由试验结果可知,为降低东北地区冻融循环作用对季节冻土区路基强度的影响,提前做好路基的防排水工作对于路基冻融病害防治是非常重要的。     

加卸载应力作用下煤岩渗透率演化模型研究

为研究开采过程中煤岩力学行为及渗透率演化规律,运用含瓦斯煤热-流-固耦合渗流伺服试验系统,进行了不同加卸载条件下原煤力学及渗流试验,分析了加卸载应力作用对煤岩变形及渗流规律的影响,得到了剪胀角随塑性剪切应变的变化关系,发现了塑性剪切应变在1.6%左右会出现剪胀角急剧变化的现象。根据试验现象和结果,考虑煤岩结构对渗透率的影响,对煤岩在弹性阶段和屈服损伤后的结构进行简化,基于火柴棍模型及渗流理论分析,从应变的角度出发,探讨了加卸载应力对煤岩渗透率的影响,建立了两个不同阶段的渗透率模型（即弹性阶段和损伤阶段）,基于不同阶段的渗透率模型和剪胀角规律,构建了煤岩全过程渗透率模型。所构建的渗透率模型与试验结果对比吻合效果较好,验证了该模型的适用性,可以为实现煤与瓦斯共采提供一定的理论依据。     

岩盐弹塑性损伤耦合模型研究

岩盐力学模型是进行能源岩盐储存工程稳定性分析的基础,而损伤和塑性机制并存且相互耦合是岩盐力学行为的基本特点。采用云应岩盐,进行了多组围压条件下的三轴压缩试验,分析了不同围压下岩盐的变形特征。在试验分析的基础上,提出了一种能够描述岩盐特性的弹塑性损伤耦合的模型,该模型描述了岩盐损伤的演化和塑性变形的耦合关系,并引入了一种非关联的塑性流动法则来描述岩盐从塑性体积压缩到膨胀的转化。采用该模型对在三轴压缩下的岩盐应力-应变关系进行了模拟分析,并与试验数据进行了对比,结果表明该模型能够较好地描述岩盐的主要力学和变形特性。     

考虑降解效应的城市固体废弃物非线性本构模型

基于城市固体废弃物(MSW)的降解及土力学特征,在应力降解模型、等向压缩特性和邓肯-张模型的理论基础上,建立了一个能够描述MSW应力和降解耦合作用的全量非线性本构模型。给出了不同应力状态下单元体体积和轴向变形随时间变化的数学表达,同时给出了所有参数的确定方法。采用该模型对不同龄期和不同围压条件下填埋场现场试样的三轴CD应力-应变曲线以及室内长期一维降解压缩试验进行了模拟。结果表明,模型概念合理、参数明确,能够较全面地描述不同降解龄期垃圾的本构关系,可用于分析填埋场的长期变形。     

不同渗透力的类充填体力学特性及损伤软化本构模型研究

通过离心试验产生的渗透压力对试样进行模拟,开展单轴压缩力学试验来研究渗透力对充填体力学特性的影响,且对充填体试样的变形特性随渗透压力变化的规律进行讨论。研究结果表明,随着渗透力增大,试样应力-应变曲线压密阶段的区间先减小后增大,弹性阶段的区间缩小,屈服阶段不明显。试验过程从低渗透压力到高渗透压力,试样的破坏模式依次表现为拉伸破坏、剪切破坏,且产生的裂纹数目增多,形态趋于复杂。在力学试验的基础上,考虑到试样压密阶段的应力-应变关系,建立不同渗透力的充填体试样损伤软化本构模型。验证结果显示,理论曲线和试验曲线高度吻合,该本构模型适用于分析不同渗透力的充填体单轴压缩力学问题。该研究为超重力离心模拟和地下渗流试验开展提供参考。     

珊瑚砂微生物固化体三轴压缩试验及损伤本构模型研究

利用微生物诱导碳酸钙沉积（MICP）技术固化南海某岛礁的陆域吹填珊瑚砂,对珊瑚砂微生物固化体进行了三轴压缩试验,基于损伤力学理论建立了珊瑚砂微生物固化体的损伤本构模型。结果表明,利用MICP技术固化珊瑚砂效果好,强度高;固化体的三轴压缩应力–应变曲线可分为近似线弹性阶段、屈服阶段与延性流动阶段。将固化体划分为匀质微元进行损伤演化分析,根据连续介质损伤力学的有效应力理论与应变等效假说,定义了损伤变量,假定固化体强度服从双参数的Weibull分布及Druker-Prager准则,建立了损伤本构模型。模型参数包括固化体力学参数和Weibull分布参数,由三轴试验和线性回归法确定,并用试验资料初步验证了模型的合理性。     

膨胀土弹塑性本构关系数值建模方法研究

按照数值建模理论,建立了不同含水率条件下膨胀土的应力-应变关系。通过可视化绘出了在整个应力场（p,q）中的三维应力-应变关系曲面和屈服轨迹。经过对比可见含水率的变化对膨胀土本构关系的影响是明显的。未使用传统的塑性势理论,而是采用数值建模方法,既克服了寻找塑性势解析表达式的困难,又可以任意选择含水率等初始条件,从而使建立的本构数值模型能够真实地反映膨胀土的变形过程     

冻融循环作用下泥质白云岩力学特性及损伤演化规律研究

针对饱水状态下的泥质白云岩,采用BCD–218 C低温数控恒温箱、WDW3100型微机控制电子万能试验机进行不同冻融循环状态下的单轴压缩试验,对其在单轴状态下应力-应变曲线变化特征以及抗压强度、峰值应变、弹性模量、泊松比的变化规律进行了分析。在现有的损伤力学理论基础上,引入斜率增大趋势系数,建立了以冻融循环次数和应变为控制变量的本构模型。结果表明,在冻融循环过程中,单轴抗压强度、峰值应变、弹性模量参数呈指数下降趋势,泊松比呈线性增加的趋势;受冻融循环的影响,应力-应变全过程曲线在第一阶段后出现短暂的线形平缓阶段,具有岩石局部滑移现象;岩石破坏形态从脆性向延性转化,具有峰前塑性硬化和峰后应变软化等行为特征。理论本构模型与试验结果相符合。     

热-水-力耦合作用下非饱和土变形特性的弹塑性模拟

在已有的非饱和土水-力耦合模型基础上,耦合考虑温度影响,建立了一个热-水-力耦合作用下非饱和土弹塑性本构模型。该模型以土骨架平均应力、修正吸力和温度为应力状态变量,以土骨架应变、饱和度和熵为应变状态变量。通过引入与温度相关的屈服面(LY、TY)以及相应的硬化规律来考虑温度对土体变形的影响。利用建立的模型,对文献中不同吸力和温度条件下的等向压缩和三轴排水剪切试验进行预测,预测结果表明,该模型能够较好地定量描述热-水-力耦合作用下非饱和土的变形特性。     

岩石应变软化变形特性及损伤本构模型研究

对滇藏铁路三江段的花岗岩进行了单轴和常规三轴压缩试验,研究了岩石在不同围压条件下的变形特性及其损伤演化规律。试验结果表明:该岩石压缩变形过程中,内部裂隙的演化经历了无起裂、起裂并稳定扩展、加速扩展、扩展变缓以及扩展停止5个阶段,且围压对岩石变形特性有显著影响。在此基础上,提出了基于有损微元生成的损伤演化模型,建立了三轴压缩条件下岩石应变软化损伤本构模型。该模型能够统一描述岩石峰前及峰后全部变形特性及其损伤演化规律,且损伤方程形式简单、参数物理意义明确。理论结果和试验数据具有很好的一致性,且参数分析结果也表明了该模型的合理性。     

考虑多晶冰软化特性的弹塑性损伤本构模型开发及应用

随着全球气候逐渐变暖,冰川、大型冰盖或较厚的冰层发生大规模滑动、崩塌灾害的可能性和危险性不断增加,研究多晶冰的变形特性并建立其本构模型对极地灾害预报具有重要的意义。基于已有的人造多晶冰常规三轴压缩试验数据,发现人造多晶冰具有明显的应变软化特征,构建以等效塑性应变为变量的损伤变量,以摩尔-库仑准则为初始屈服函数,采用相适应流动法则,建立了能够反映多晶冰应变软化特征的弹塑性损伤本构模型。利用不同围压、不同温度条件下的常规三轴压缩试验数据对提出的本构模型进行了验证,并对部分模型参数进行了参数敏感性分析。将建立的弹塑性损伤本构模型写入FLAC^3D软件中,利用单个单元的三轴压缩数值模拟对写入的本构模型进行了验证;然后使用写入的弹塑性损伤本构模型对圆柱状试样开展三轴压缩数值模拟,对数值模拟结果进行了详细的分析。研究结果以期为冰川稳定性分析以及冰川多场耦合研究提供理论及数值基础。     

饱和度对非饱和土力学性质的影响

现在被广泛公认的由Fredlund提出的非饱和土力学的双参数理论,即净应力和吸力为非饱和土的应力状态量,不能直接考虑饱和度或含水率对非饱和土的应力-应变关系和强度的影响。在非饱和土三轴试验结果表明,即使在净应力和吸力路径相同的条件下,具有不同饱和度试样的应力-应变关系和强度也是不同的。其他条件相同时,试样饱和度越高,其应力比-应变关系曲线越高,强度越大。最新的水力-力学特性耦合的弹塑性本构模型可以定量地表示上述非饱和土的性质