考虑变形模量劣化的应变软化模型在FLAC3D中的开发与验证

FLAC3D软件中的应变软化本构模型是通过抗剪强度参数劣化来实现的,没有考虑变形模量劣化造成的影响。针对这一问题,采用了变形模量劣化的损伤本构模型,详细推导了劣化损伤本构模型的差分格式,研究该模型在FLAC3D中二次开发的流程,在VC++环境下实现了劣化损伤本构模型在FLAC3D中的二次开发;通过与深埋圆形洞室应变软化数值解算例对比,进一步验证了该模型弹塑性和损伤的力学性质,分析了变形模量劣化对围岩变形特性的影响。由于该模型能够反映围岩的劣化损伤特性和大变形特点,因此在工程应用中有一定的实用价值。     

一种求解应变软化岩体中球腔问题的数值方法

均匀地应力场中球形洞室开挖问题可以简化为空间轴对称问题,为分析和评价球形洞室开挖的稳定性,采用有限差分方法对平衡方程与变形协调方程进行离散,将洞室围岩潜在塑性区划分为一系列的球壳,并考虑岩体的应变软化行为,可以得到满足线性Mohr-Coulomb（H-C）准则和非线性Hoek-Brown（H-B）屈服准则的围岩体塑性区范围及相应的应力与变形分布。计算结果表明,随着离散步的增加,数值解逐渐收敛,文中解与已有文献的结果完全吻合,且计算效率更高。对一组强度参数遵循CWFS（黏聚力弱化与摩擦角强化）规律的围岩体内塑性区范围、变形与应力的分布进行了分析,并与常规的应变软化模型的计算结果进行了对比,CWFS模型的结果更合理一些。     

基于Hoek-Brown准则的非等比强度折减方法

针对节理岩质边坡稳定性分析问题,开展了基于Hoek-Brown破坏准则的强度折减法研究。边坡的失稳过程是一个渐进累积破坏过程,破坏过程中各强度参数的衰减程度不同,在强度折减法中对应的折减系数也不同,不同折减系数间关系的确定及其综合安全系数的定义值得深入探讨。首先从岩体材料软化（硬化）特性出发,根据岩体强度参数从峰值强度到残余强度的变化规律,推导了Hoek-Brown破坏准则中3个强度参数折减系数间的数学关系式,然后以折减前后滑面上抗滑力之比定义综合安全系数来评价边坡的稳定性,最后结合算例验证了该方法的合理性,可有效应用于节理岩体边坡的稳定性分析。     

模拟应变软化岩石三轴试验过程曲线

在经典弹塑性理论框架下,认为岩石材料在应变软化过程中遵守摩尔-库仑强度准则,塑性变形服从非关联流动法则。对两组岩石在三轴试验中呈现的应变软化行为进行了分析,获取了峰值前后强度参数的大小。采用应变软化模拟方法,将应力-应变曲线峰后应变软化阶段简化为一系列应力跌落和塑性流动的脆塑性过程,在给定的强度参数演化规律基础上,得到了岩石的应力-应变曲线、侧向应变-轴向应变曲线、体应变-轴向应变曲线和塑性体应变-轴向塑性应变曲线。模拟得到的应力-应变曲线与试验结果吻合较好,其他曲线与相关试验测试规律基本一致。研究成果对确定岩石强度参数与认识岩石峰后强度演化规律有一定的指导意义。     

基于统一强度理论深埋圆形隧道围岩的剪胀分析

基于连续介质理论中岩体的剪胀角与围压和塑性剪切应变密切相关,隧洞周边岩体的应力状态因开挖卸荷而发生应力重分布,迫使其围压由原地应力逐渐衰减,塑性剪切应变不断增加,引起剪胀效应呈非线性变化。首先,基于统一强度理论和非关联流动法则,将潜在塑性区围岩按等围压释放划分为若干同心圆,提出了考虑中间主应力和非线性剪胀性的有限差分法,计算应变软化围岩的力学问题,并以实例验证其正确性。其次,通过参数分析,研究塑性区内岩体的剪胀角受中间主应力、临界软化系数和支护力的影响规律。研究结果表明,中间主应力主要影响剪胀角的峰值,随着中间主应力效应增加,剪胀峰值增加;临界软化系数主要影响剪胀角的变化率,随着临界软化系数的增加,剪胀角变化缓慢;中间主应力和临界软化系数共同影响塑性区剪胀角的变化;随着支护力的增加,洞壁处的剪胀角增加;双剪强度理论计算的位移值较小,应谨慎采用,同时采用Mohr-Coulomb强度准则时可以适当考虑围岩的承载潜力。     

土坝坝体底部劈裂灌浆加固效果研究

介绍了传统堤坝劈裂灌浆工存在的主要问题,通过对灌浆起始劈裂点的力学求解,分析了劈裂灌浆浆液扩散规律;提出了坝体底部首先劈裂的灌浆工艺。以山东某水库堤坝劈裂灌浆工程为研究对象,利用三维有限差分数值方法分析了坝体底部劈裂灌浆引起坝体的应力和变形特征。结果表明,采用应变软化模型能有效的模拟土体材料的弹塑性力学行为,反映灌浆完成后浆脉的固结和周围土体的湿陷作用,即“浆坝互压”的作用机制,分析结果与实测结果吻合较好;同时得到了一些有价值的结论,为劈裂灌浆防渗加固工程的施工和设计提供理论指导。     

考虑界面软化特性的垃圾填埋场斜坡上土工膜内力分析

为分析垃圾填埋场斜坡上防渗系统中的土工膜在上覆垃圾堆体重力作用下产生的拉应力,考虑土工膜-黏土界面的应变软化特性,通过建立土工膜受力的平衡方程获得了土工膜拉应力的解析解。与传统方法相比,该方法可分析土工膜-黏土界面不同位置所处的3种应力状态：弹性状态、软化状态及残余状态,从而使得垃圾填埋场斜坡上土工膜的受力分析更为合理。参数分析显示,当软化阶段的强度损失及残余位移与峰值位移的差值较小时,土工膜-黏土界面的峰值强度对土工膜最大拉应力的影响较小,而其残余强度是主要控制因素;而当土工膜-黏土界面的软化特征较明显时,则峰值强度和残余位移的影响也更为明显。另外,土工膜上覆堆体高度、斜坡坡度及土工膜锚固长度对土工膜拉应力也有较大影响。     

土与结构物接触面损伤本构模型

根据粗糙接触面变形机理,建立了基于损伤力学基本原理的接触面本构模型来描述其力学特性。所建议的损伤模型可以用一个统一的表达式来表述,能够较好地反映土与结构物接触面剪切过程中的应变软化和剪胀等力学特性;模型参数较少,物理意义明确。对接触面直剪试验和单剪试验成果进行了验证,模型计算与试验结果吻合较好,表明建议的土与结构物粗糙接触面损伤模型是合理可行的。     

煤岩两体模型变形破坏数值模拟

采用拉格朗日元法,在弹性岩石与弹性-应变软化煤体所构成的平面应变两体模型的上、下端面上不存在水平方向摩擦力条件下,模拟了模型的破坏过程、岩石高度对模型及煤体全程应力-应变曲线、煤体变形速率、煤体破坏模式及剪切应变增量分布的影响。结果表明,当模型的全程应力-应变曲线达到峰值时煤体内部的剪切带图案已经十分明显,在模型的应变硬化阶段,煤体中的应变局部化可视为模型失稳破坏的前兆,随岩石高度的增加,模型应力-应变曲线的软化段变得陡峭,这与单轴压缩条件下的解析解在定性上是一致的;煤体应力-应变曲线的软化段变得平缓,煤体消耗能量的能力增强;弹性阶段煤体的变形速率降低;煤体内部的剪切应变增量增加。煤体应力-应变曲线的软化段的斜率、弹性阶段煤体的变形速率、煤体内部的剪切应变增量及塑性耗散能都受岩石高度的影响,说明了岩石几何尺寸对煤体的影响(煤岩相互作用)是不容忽视的。     

考虑剪胀效应的深埋圆形隧道围岩应变软化弹塑性解

隧道围岩非线性体积膨胀影响施工安全,如何正确评价隧道围岩塑性区域内扩容机制非常重要。多数研究从理论上推导深埋圆形隧洞应变软化围岩应力-应变场时仅采用恒定或线性变化剪胀模型,针对此不足,提出了一种基于有限差分法的分析方法,能够合理地考虑围岩非线性剪胀效应及其应变软化特性。利用已有的研究成果,验证了该方法的合理性,并进一步探讨了不同质量石灰岩与支护压力下剪胀系数在围岩塑性区域内的影响因素,比较了恒定与非线性变化剪胀模型下围岩变形的不同。结果表明:对于地质强度指标(GSI)较小、质量较差的岩体,塑性区域内主要由围压控制剪胀效应程度;恒定与非线性变化剪胀模型下围岩洞壁变形差别显著。