随机分布贯穿裂隙岩体峰后应力-应变关系模型

基于峰后软化阶段强度参数逐渐演化这一行为,采用Mohr-Coulomb强度准则,分别以岩石的最大主应变和裂隙的切向位移作为软化参数,假设强度参数为软化参数分段线性函数的条件下,给出岩石的峰后应力-应变关系和裂隙的峰后应力-切向位移关系的求法,在此基础上,利用概率方法给出了随机分布贯穿裂隙岩体峰后应力-应变关系式的求法.最后结合算例讨论裂隙的平均间距、法向刚度和切向刚度对峰后应力-应变曲线的影响.研究结果表明:裂隙的平均间距、法向刚度和剪切刚度越小,裂隙岩体的应变越大.     

围压作用下岩石峰后应力-应变关系模型

假定峰后软化阶段任意一点皆满足摩尔库仑极限破坏条件,建立岩石峰后后继屈服面模型。利用岩石峰值处和残余阶段初始处应力状态绘制几组莫尔应力圆,获得不同围压下峰值处和残余阶段初始处的黏聚力c和内摩擦角φ,然后分别建立黏聚力c和内摩擦角φ同围压的函数关系。基于莫尔-库仑准则,以最大主应变ε作为软化参数,以分段线性函数形式演化黏聚力c和内摩擦角φ的正切,反求内摩擦角φ,从而建立岩石峰后应力-应变关系模型。数值算例研究Tennessee大理岩三轴试验的峰后应变软化过程,模拟结果与试验结果较吻合。研究结果表明:建立的峰后应力-应变关系模型合理而可靠,可以较好地描述不同围压下岩石的峰后力学行为。     

岩石软化时效变形的弹黏塑性模型表征

为同时描述岩石应变软化以及时间相依的变形特征,基于超固结土的下负荷面模型,依据相对过应力思路,以Hashiguchi的下负荷面为参考屈服面,建立一种能表征岩石软化时效变形特征的弹黏塑性本构模型。相对于原下负荷面模型,新模型增加了2个能通过不同应变速率三轴压缩试验测定的率敏性参数c0和m',并结合岩石的特点用弹性模量E0和变形模量Ed代替原模型中的回弹指数κ与压缩指数λ。新模型采用以当前应力、黏塑性应变以及黏塑性应变速率为状态变量的动屈服准则函数,给出了基于Newton-Raphson迭代的应力积分算法,且成功地将其嵌入到大型有限元软件ABAQUS中。利用新模型对岩石率敏性三轴试验进行模拟验证,结果表明:预测结果与试验结果吻合良好;模型参数物理意义清晰、简单可测,能正确描述岩石的率敏性及应变软化等特征,可用于复杂边值问题的有限元计算。     

峰后脆性对非均质岩石试样破坏及全部变形的影响

在平面应变压缩条件下,采用FLAC模拟峰后脆性对含初始随机材料缺陷的岩样的破坏过程、前兆、声发射及全部变形特征的影响.利用若干FISH函数预置初始随机缺陷,计算全部变形特征,并统计每10个时间步内的破坏单元数.密实的岩石服从莫尔库仑剪破坏与拉破坏复合的破坏准则,破坏之后呈现应变软化-理想塑性行为.缺陷在破坏之后经历理想塑性行为.不同脆性岩石的峰后的应力-轴向应变曲线、应力-侧向应变曲线、侧向应变-轴向应变曲线、体积应变-轴向应变曲线及由侧向应变及轴向应变计算得到的泊松比-轴向应变曲线于应力峰值之前发生分离.由于缺陷的依次破坏,在初始加载阶段之后,计算得到的泊松比随着轴向应变的增加而直线上升,这使侧向应变增加的速度超过轴向应变增加的速度.随着峰后脆性的降低,岩样失稳破坏的前兆变得明显,破坏变得不突然.脆性越强的材料,体积膨胀之后发生应变软化越早;剪切应变集中的位置越少;在硬化阶段及软化阶段,声发射持续的时间越短.     

基于强度参数演化的裂隙岩体峰后应力-应变关系求法

完整岩块峰后应力应变关系的表达不能较好反映裂隙岩体峰后特性,为得到适合裂隙岩体峰后应力-应变关系,在试验数据基础上,基于库伦强度准则,分析裂隙岩体峰后强度参数的演化行为。根据裂隙岩体在不同围压和裂隙倾角时的不同峰后特性,将峰后应力-应变关系简化为A和B两种类型。以应变作为应变软化参数,假设黏聚力、内摩擦角为应变的分段线性函数,给出裂隙岩体峰后两种类型应力-应变关系式的求法。将A和B两种类型应力-应变关系式算出结果与试验结果进行比较,得出用裂隙岩体的应力-应变关系求得峰后应力-应变曲线与试验数据基本一致。     

片状岩石全应力应变压缩过程及力学机制试验研究

为了解片状岩石的全应力应变压缩过程及力学机制,利用RMT-150C岩石力学刚性伺服试验系统,对绿泥石片岩和云母石英片岩进行压缩力学试验.基于试验结果,得到岩石全应力应变压缩过程,提出峰前压密、线弹性、微裂隙稳态扩展和微裂隙不稳态扩展4个变形阶段以及峰后应变软化、残余应力稳定两个变形阶段,分析岩石在各阶段的强度变形特性以及体应变、横向应变和轴向应变之间的关系;详细探讨压缩条件下水对试样物理力学性质的影响规律,重点对比峰值强度、弹性(变形)模量等指标的敏感程度和围压效应;并揭示片状岩石在压缩状态下的破裂模式,将整个破坏过程分为单一剪切破裂、鼓状破裂、Y字形破裂及平行劈裂4种形式.     

堆石颗粒浸水破碎引起堆石料湿化变形研究

为研究堆石料的湿化变形,首先通过室内颗粒破碎试验分析了粒径对岩石颗粒强度和软化系数的影响;然后依据散粒集合体应力与应变张量表达式和岩石颗粒的软化系数,建立了饱和与干燥试样的应力-应变转换关系,结合双线法预测了堆石料的湿化变形;最后通过算例分析了该方法的有效性.结果表明：浸水湿化会降低岩石颗粒的强度,饱和与干燥玄武岩颗粒的强度均服从Weibull分布;颗粒特征强度随粒径的增加而降低,但是Weibull模量和颗粒的软化系数没有明显的粒径相关性.堆石料饱和与干燥试样的应力-应变关系可通过岩石颗粒软化系数进行变换;由饱和试样的应力-应变关系和岩石颗粒软化系数可估算堆石料浸水后因颗粒破碎引起的湿化变形.算例中预测结果与试验结果整体比较接近,表明岩石颗粒浸水后强度降低导致的颗粒破碎是堆石料产生湿化变形的主要原因.     

单轴压缩条件下普通混凝土柱的峰后非线性尺寸效应

将单轴压缩条件下遭受到剪切带(峰值应力之后呈现线性应变软化行为)形式的单一剪切破坏普通混凝土试样的峰后应力-应变曲线斜率的解析解推广为非线性情形.在峰值应力之前,采用Scott模型描述非线性的本构关系.剪切带的非线性应变软化本构关系由导出的最短普通混凝土试样峰后斜率反算.利用得到的峰后本构关系,对其他较长的普通混凝土试样的应力-应变曲线进行了预测.预测的峰后应力-应变曲线依赖于试样的高度,且与实验结果吻合.估算的剪切带内部平均塑性剪切应变远大于在单轴压缩条件下测得的轴向应变的极限值.若测得的峰后应力-应变曲线被视为本构关系,则普通混凝土柱的峰后延性将被极大地低估.     

应变空间的软岩统一弹塑性软化本构模型

为使以拉为正的统一强度理论表达式适用于岩土工程界以压为正的习惯表达，首先将以压为正的单向压缩和单向拉伸试验条件代入统一强度理论表达式，确定了相关参数，建立了以压为正的统一强度理论表达式，该表达式考虑了材料拉压不等及中间主应力的影响．将所建立的统一强度理论表达式引入应变空间，针对一种硅藻质软岩材料，利用基于应变空间的弹塑性本构方程对软岩的软化模型进行了分析和探讨，建立了应变空间表述的三轴状态下的弹塑性本构关系式，并结合硅藻质软岩的三轴试验结果进行了数值模拟分析，结果表明应变空间的统一弹塑性软化本构方程的数值模拟结果非常接近试验结果，可以很好地模拟软岩的应变软化特征．     

基于FLAC3D模拟的岩石单轴 压缩试验应力应变统计规律

在三维数值模拟的基础上,利用FLAC~(3D)模拟了岩石单轴压缩试验过程.根据全应力应变曲线的4个破坏阶段(弹性阶段、微破裂稳定发展阶段、非稳定破裂发展阶段和破裂后阶段),采用内嵌的FISH语言编程提取出模型中相应破坏阶段的全部应力应变数据,并利用统计分析软件@Risk对大量的应力和应变数据进行统计分析,以研究岩石变形各阶段中试样整体的应力和应变概率分布规律.本文不仅验证了大量岩石微单元的受力服从Weibull分布的规律,还得出了径向和轴向应变分别服从Logistic分布和BetaGeneral分布的结论.     

脆性颗粒材料的超弹性软化模型及应变局部化模拟

在Volokh软化模型的基础上,考虑颗粒材料的变形和破坏机制,认为颗粒材料以剪切变形和破坏为主,提出了表示软化程度的软化因子概念,发展了两种能够部分表征颗粒材料力学行为特征的软化模式,并基于ABAQUS用户材料子程序(UMAT)开发相应程序.数值算例考察了上述模型模拟颗粒材料应变局部化现象和软化现象的能力.结果表明了所建议模型的有效性和可行性.     

混凝土软化特性的数值模拟

采用微裂纹区分段线性应变软化律，并考虑由微裂纹演化而造成的弹性刚度劣化对裂缝扩展的影响，在给定σ－ω软化律的基础上，用应变增量迭代法，数值模拟了混凝土三点弯曲梁的Ｐ－△全曲线。     

考虑粗粒土应变软化特性和剪胀性的本构模型

在传统细粒土本构模型的基础上进行了改进,建立了一个能较好地描述粗粒土应变软化特性和剪胀性的弹塑性本构模型.该模型采用双屈服面形式,克服了单屈服面模型存在的一些不足,可同时反映剪切变形和压缩变形机理.模型在应变软化特性描述方面,提出了一个利用残余状态应力比和峰值应力比的应变软化公式,较为合理地反映了粗粒土的应变软化现象.在剪胀性描述方面,考虑了状态转换应力比与初始有效围压的相关性.数值模拟结果与试验结果较为接近,表明该模型描述粗粒土在低围压和相对中高围压下的应变软化特性和剪胀性方面具有一定的优越性.     

一种受剪细观损伤单元模型及其应用

基于剪切损伤破坏机制,发展了相应的细观损伤单元模型,它由弹脆性的微弹簧和模拟裂缝滑移的塑性元件组成.混凝土单轴受压时的宏观性能可用受剪细观损伤单元的并联体来解释,由此建立了混凝土单轴受压的损伤本构关系模型,并建议了塑性变形的简化计算方法.通过引入混凝土剪切单元破坏应变分布随机场,推导了应力与损伤的均值、标准差的表达式.根据混凝土单轴受压时的应力-应变全曲线试验结果,结合随机建模原理和优化算法确定随机场参数和材料参数,将分析得出的受压随机应力-应变关系与试验结果进行比较,二者吻合良好.     

软岩应变软化特性的数值解析探讨

基于实验研究,提出了一种新的考虑时间因素的增量形应力－应变－时间关系,用来描述硅藻软岩的应变软化特性．分别对正常固结和超固结领域软岩试样的应力－应变关系,特别是应变软化特性进行了数值模拟计算,得出了与实验相近的结果．通过对计算结果和实验结果的比较分析,探讨了该计算模型的合理性和可行性,验证了软岩的应变软化特性受应力应变时间依存性的影响．     

颗粒材料变形破坏与影响因素细宏观分析

考虑到颗粒流(PFC)表示的颗粒速度场和位移场并不能很直观地表示出颗粒材料在加载破坏中的变形特征,从宏观的角度引入易于理解和接受的颗粒集合体的有效应变定义,建立了考虑不同周围压力、颗粒粒径、粗糙程度、接触刚度等因素在内的颗粒材料数值模型,并结合一定的初边值条件,编制了相应的Fish语言程序.将细观数值分析与宏观表示相结合,分析了双轴压缩试验中颗粒材料变形破坏的规律和各种基本的细、宏观影响因素.数值结果表明,基于宏观连续体概念的有效应变充分显示了颗粒材料变形破坏的发展形成过程及相关特征:破坏区域完全形成于应力-应变曲线峰值之后的软化阶段,域内有较大的孔隙,域中的颗粒显示出明显的颗粒旋转;颗粒的法向与切向刚度、粗糙程度、颗粒集合体试样所受的压力水平,以及颗粒粒径等对颗粒材料的承载能力,破坏区域的形成、型式与厚度等有一定的影响.     

关于非线性非局部应变梯度模型的探讨

考虑线性软化模型不能反映材料软化变形复杂性的不足,选取高斯正态分布函数作为非局部理论的权函数,同时采用指数型应力衰减模式考虑软化的非线性特征,进一步地将这种非线性非局部理论与采用Laplace算子考虑塑性应变梯度效应的塑性梯度理论相结合,建议了一种非线性非局部应变梯度模型,并据此针对各向同性单向拉伸杆的力学响应进行了分析,将计算结果与线性软化模型进行了对比.结果表明:在非线性软化模型中,尽管塑性应变分布形式与线性软化模型的相同,但塑性应变同时依赖于应力降与破坏极限应变,两者最大应变是不同的.     

混凝土应变软化关系及其确定

在对混凝土应变软化关系物理意义及试验方法进行分析的基础上，提出一种新方法确定混凝土应变软化关系；这种方法根据比较容易获得的试验结果，通过数值计算，间接确定混凝土的应变软化关系．实际应用结果表明这种方法是比较简便有效的