双轴应力下非贯通节理岩体压缩损伤本构模型

针对目前节理岩体损伤变量定义中大多仅考虑节理长度、倾角等几何性质,而未考虑节理抗剪强度等力学性质的不足,基于断裂力学中的由于单个节理存在引起的附加应变能增量与损伤力学中的损伤应变能释放量相关联的观点,推导出了在双轴应力下含单条非贯通闭合节理岩体的损伤变量计算公式,并根据断裂力学理论对双轴压缩荷载下的单个节理尖端应力强度因子计算方法进行了研究,得出了应力强度因子KⅠ、KⅡ的计算公式;同时考虑多节理间的相互作用,给出了单组单排及单组多排非贯通节理尖端应力强度因子计算公式,由此建立了相应的节理岩体双轴压缩损伤本构模型,并利用该模型进行了相应的算例分析。结果表明：对含单条非贯通闭合节理的岩体而言,当节理倾角小于其内摩擦角时,岩体强度与完整岩石相同,岩体损伤为0,而后随着节理倾角增加,岩体强度、损伤随节理倾角的变化分别呈开口向上及向下的抛物线,当节理倾角约为60°时,岩体损伤最大,强度最低。随着节理长度增加,岩体损伤增加,而随着节理内摩擦角的增加,岩体损伤则减小;对含单组单排非贯通闭合节理的岩体而言,当节理总长度一定时,随着单条节理长度的减小及节理条数的增加,岩体损伤则逐渐减小,但其减小幅度与节理条数并不呈线性关系。     

非贯通节理岩体损伤变量计算方法研究

作为岩体组成部分的非贯通节理对岩体力学特性有着重要影响,然而几乎目前所有的岩体损伤变量计算方法都仅考虑节理几何参数对岩体力学特性的影响。对含单组非贯通节理的岩体力学特性进行研究,提出一个能够同时考虑节理几何及强度参数对岩体力学特性影响的新的岩体损伤变量计算方法。首先,采用弹性余能等效假设代替Lemaitre应变等效假设研究由节理引起的岩体各向异性损伤,并基于断裂力学中单个节理引起的附加应变能增量与损伤力学应变能释放量相关联的观点,推导出由单条节理引起的损伤变量计算公式。其次,根据断裂力学理论获得了单轴压缩下单条节理尖端应力强度因子（SIF）KⅠ、KⅡ的计算公式。最后,通过考虑节理间的相互作用给出了单组单排或多排节理尖端应力强度因子KⅠ、KⅡ的计算公式,得到了单组节理引起的岩体损伤变量计算公式,并与已有试验结果的对比分析证明了该公式的合理性。     

基于宏细观损伤耦合的非贯通裂隙岩体本构模型

针对非贯通裂隙岩体工程结构中的受荷岩体,提出受荷细观损伤与裂隙宏观损伤的概念。以完整岩石的初始损伤状态作为基准损伤状态,综合考虑裂隙宏观缺陷的存在、微裂纹细观缺陷在受荷下的损伤扩展以及宏细观缺陷在受荷过程中的耦合,基于Lemaitre应变等效假设,推导了考虑宏细观缺陷耦合的复合损伤变量,并给出同时考虑试件尺寸、裂隙几何与力学特性的宏观损伤变量的计算公式,从而建立了基于宏细观缺陷耦合的非贯通裂隙岩体在荷载作用下的损伤本构模型。用宏细观损伤耦合的本构模型来描述非贯通裂隙岩体在受荷过程中的细观损伤演化与宏观损伤行为,与非贯通裂隙岩体实际受荷情况符合较好。研究结果表明：（1）完整岩样和裂隙岩样的应力-应变行为在峰值强度之前差异较大,峰值强度以后差异逐渐减小,最后趋于一致,二者具有相近的残余强度;（2）裂隙岩体强度随裂隙贯通率的增加而增大,随裂隙倾角的变化具有明显的各向异性,同时还与裂隙面的内摩擦角有关;（3）裂隙倾角为90°时,裂隙岩样的峰值强度最高;张开型裂隙岩样的裂隙倾角为45°时,峰值强度最低;（4）非贯通裂隙岩体工程结构中的受荷岩体,其力学性能由受荷细观损伤与裂隙宏观损伤及其耦合效应所决定,基于宏细观损伤耦合的复合损伤变量可以较好地反映非贯通裂隙岩样的力学特性。     

节理岩体本构模型与强度理论的统计断裂力学分析

本文提出了节理岩体的三维本构模型、强度理论及破坏概率的表达式。这里把节理岩体视作等效的各向异性连续介质,其具有裂隙张量S_(ifst)和柔度张量C_(ifst)。为了获得本构模型、强度判据及破坏概率,这里引用了能量平衡方程ω=ω_0+ω_c、断裂力学能量判据、岩体节理统计理论及可靠性理论,并引入了剩余剪应力τ_n,有效应力σ_(if)及最可能出现的最大节理半径a_m等概念。这套理论的初步应用是令人满意的。     

综合考虑宏微观复合损伤的节理岩体本构模型

针对目前节理岩体损伤本构模型中仅考虑节理等宏观缺陷造成的损伤,而没有考虑岩块中微裂纹等微观损伤的不足,提出了综合考虑宏微观损伤的节理岩体本构模型。其中微观损伤模型采用基于应变强度理论和岩石微元强度服从Weibull分布假定的统计损伤模型,把其应用于被节理切割而成的岩块。而宏观损伤模型主要基于连续介质力学原理,把节理对岩体性质的劣化作为损伤来考虑。在考虑节理传压及传剪系数的基础上,提出了综合考虑宏微观复合损伤的节理岩体本构模型。算例表明宏观节理的存在大大削弱了岩体的强度,降低了其刚度,增大了其柔性。所提出的模型能够较好地反映宏微观两类损伤对岩体应力-应变关系及强度的影响,较为合理。     

考虑损伤的节理本构模型

本文在弹塑性损伤的理论框架内,讨论了节理等地质间断面的本构模型。这个模型能够反映节理面的损伤弱化,扩容和弹性刚度劣化等复杂特性。这个模型的另一优点是,塑性变形增量与屈服面是非正交的,但本构矩阵具有对称性。这种对称性在岩石力学的理论研究和数值分析中是至关紧要的。     

单轴受压状态下钢纤维混凝土损伤本构模型研究

针对钢纤维混凝土在单调载荷条件下损伤行为的复杂性,根据Weibull统计分布理论和等效应变假定原理,得出了钢纤维混凝土损伤本构模型,基于试验数据,确定了模型参数。分析计算结果显示,该模型可以有效地描述单轴受压状态下钢纤维混凝土的损伤过程。     

考虑裂隙变形参数的岩体单轴压缩损伤模型

为了更准确地预测裂隙对岩体压缩力学特性(如强度和刚度)的影响,需要建立更为合理的裂隙岩体压缩损伤模型。为此,基于相关试验数据和裂隙岩体单轴压缩力学行为,采用损伤及断裂理论对目前断续裂隙岩体压缩损伤模型中存在的不足进行了深入分析,并对其进行了3方面的改进,即：不再将裂隙变形参数视为定值、考虑了裂隙面上法向正应力产生的负第一应力强度因子(K_I)、考虑了裂隙面上有效剪应力产生的K_I,由此最终提出了考虑裂隙变形参数的岩体单轴压缩损伤本构模型。最后采用试验数据对该模型的合理性进行了验证,发现与现有模型相比,该模型明显提高了岩体单轴压缩弹性模量和损伤值的预测精度,尤其是当裂隙倾角为0°时,该模型计算得到的弹性模量为4.306 MPa,与实测弹性模量4.310 MPa几乎相同。因此,该模型能够很好地刻画岩体单轴压缩力学行为,这也说明考虑裂隙变形参数对岩体单轴压缩力学特性的影响是十分必要的。该研究可为准确预测裂隙岩体的单轴压缩力学行为提供参考。     

基于复合损伤的节理岩体动态本构模型研究

基于运用霍普金森压杆（SHPB）装置对节理岩体动载试验得出的数据,从节理面倾角、贯通度、厚度、组数、填充物及应变率等不同方面分析各因素对节理岩体力学特性的影响,通过对节理岩体在高应变率下的损伤机制和破坏形式进行分析的基础上,基于复合损伤理论,对广义Bingham模型进行改进,构造了节理岩体材料在不同应变率下动态响应的本构模型,对模型计算和试验结果的比较表明,该模型能够很好地描述节理岩体动荷载下初始弹性变形阶段、稳态塑性变形阶段和加速变形破坏阶段的应力-应变关系,并且理论与试验结果吻合较好,从而证明了该模型的正确性和合理性。     

基于摩尔-库仑准则的断续节理岩体复合损伤本构模型

针对工程实际中断续节理裂隙岩体的损伤本构模型,假设岩石微元强度服从Weibull随机分布,以摩尔-库仑破坏准则作为描述微元强度的表示方法,推导出细观损伤变量。利用能量和断裂力学理论,综合考虑节理几何特征及力学特性,推导宏观损伤变量计算公式。基于Lemaitre应变等效假设,考虑宏细观缺陷耦合作用,推导出复合损伤变量,建立了基于摩尔-库仑准则的宏细观缺陷耦合作用的断续裂隙岩体损伤本构模型。研究结果表明：（1）采用摩尔-库仑准则作为描述微元强度的统计分布变量建立的损伤模型能够较好地反映岩石内部缺陷分布和变形特征,该模型真实地反映岩石微元强度受应力状态的影响。（2）该模型建立的理论曲线与断续节理岩体的试验曲线吻合较好。（3）节理裂隙岩体宏观损伤变量及峰值强度随节理倾角的变化规律与综合考虑宏细观耦合作用下的损伤变量及裂隙岩体峰值强度随节理倾角的变化规律基本一致。（4）宏细观耦合作用下的等效弹性模量与节理贯通率呈非线性负相关;在节理倾角一定的情况下,损伤变量与节理长度呈非线性正相关;在贯通率较小时,岩体的宏观损伤变量与内摩擦角的关系呈线性负相关变化,贯通率达到一定程度,线性关系变成非线性关系。     

基于动态CT识别的冻土单轴压缩损伤本构模型

采用连续介质力学与热力学方法,建立了冻土单轴压缩损伤本构模型。以兰州冻结黄土为例,基于冻土单轴压缩动态CT试验,对冻土附加损伤的2个阶段,即塑性损伤和微裂纹扩展损伤分别采用硬化曲线法和动态CT识别方法进行了损伤计算,在此基础上给出了冻土附加损伤与弹性应变的对应关系,建立了基于试验的冻土损伤演化方程,并进行了有效应力的计算。     

岩石声速与其损伤及声发射关系研究

首先,建立了单轴压缩过程岩石损伤参量、应变与声速之间的定量关系式,分析了不同均质度对单轴压缩过程岩石声速的影响,结果表明,随着均匀度的增加,单轴压缩过程中声速由平缓变化到急剧变化,这与已有的岩石声发射数值模拟分析结果是一致的;其次,根据建立的声速与应变的关系公式,通过单轴压缩过程岩石声速与应变实测结果的回归分析,得到了具有较高精度的回归方程,从而通过试验验证了所建立的关系式的正确性;最后,从损伤的角度讨论了单轴压缩过程岩石声速与声发射的关系,得出了Kaiser点应位于声速初始下降点附近的结论,为岩石声发射测量地应力试验中Kaiser点的确定提供了新的方法。     

拉应力条件下岩石细观力学本构模型和渗透系数张量研究(I):各向异性损伤本构模型

在拉应力作用下的考虑渗流的岩石本构模型研究尚且很少。采用细观力学研究了岩石损伤变形过程,考虑了渗流过程。由于选用了新的微裂纹扩展准则,提高了解方程效率;采用细观力学成果建立了考虑渗流的微裂纹面不连续变形对岩石变形影响规律,并采用断裂力学成果建立了考虑渗流的微裂纹应力强度因子;建立了同时考虑微裂纹扩展过程和渗流过程的岩石各向异性损伤本构模型,该模型也适用于岩体。     

基于考虑动力本构的连续−非连续方法的单轴压缩岩样开裂过程模拟

动载作用下岩石的破坏规律研究对于众多地质灾害的机制分析和预防具有重要的理论及实际意义。鉴于数值模拟研究的优势,应大力发展适于岩石动力断裂过程模拟的数值方法。在自主开发的拉格朗日元与离散元耦合连续−非连续方法的基础上,采用朱−王−唐本构模型取代了广义胡克定律,发展了考虑动力本构的连续−非连续方法,其正确性通过模拟不同加载速度时砂岩试样的单轴压缩试验进行了验证。通过统计裂缝区段数目随着岩样的纵向应变的演化规律,并监测岩样左、右对称线上多个测点的最小主应力的演化规律,开展了不同加载速度时单轴压缩花岗岩试样的变形−开裂过程研究,阐明了岩样的开裂机制。研究发现,剪裂缝以雁列式展布,整体上形成剪切带。随着时步数目的增加,各测点的最小主应力均呈波动下降−震荡上升的变化趋势。震荡上升阶段对应岩样的应变软化阶段。测点分离后最小主应力的震荡幅度较大,这是由于节点分离和单元接触激发了较大的应力波。剪切带尖端的最小主应力集中会使测点发生剪切分离。当岩样的三角块向下楔入时,下方测点的应力状态类似于紧凑拉伸试验进而发生拉伸分离。