不同含水率软弱结构面剪切蠕变试验及模型研究

软弱结构面的蠕变特性是岩体重要力学特性之一,其蠕变特性常控制着岩体的蠕变变形和长期强度。软弱结构面中充填物的含水率是影响其蠕变特性的重要因素之一,通过开展不同含水率砂岩软弱结构面剪切蠕变试验,得到了具不同含水率的砂岩软弱结构面剪切蠕变试验曲线。基于得到的剪切蠕变试验曲线,分析了含水率对软弱结构面剪切蠕变特性及长期抗剪强度的影响规律。最后通过模型辨识,采用改进的粘弹塑性模型来描述考虑含水率变化的砂岩软弱结构面蠕变特性。     

非规则砂岩节理剪切变形本构关系试验研究

岩石节理剪切变形对岩体工程的安全性和稳定性具有重要影响。为研究常法向应力下岩石节理剪切变形本构关系,采用RDS-200型岩石直剪仪对非规则砂岩节理进行了不同法向应力下的直剪试验。根据岩石节理剪切应力?位移全程曲线形状特征,将其依次划分为峰前压密阶段、线性阶段、屈服阶段和峰后软化阶段;根据剪切应力在峰后软化阶段降低幅度和速率大小,将岩石节理剪切应力?位移曲线划分为3种类型：峰后平台型、峰后缓降型和峰后跌落型。基于岩石节理剪切应力?位移曲线各阶段剪切变形特征,采用分段函数建立了岩石节理剪切变形本构模型。与其他模型相比,新提出的岩石节理剪切变形本构模型对试验数据拟合精度更高,更好地描述了岩石节理剪切应力?位移全程曲线。另外,在通过岩石节理直剪试验由经验公式确定模型参数之后,所提出本构模型可在不同法向应力下实现对不同粗糙度岩石节理剪切应力?位移曲线的预测。研究结果对岩石节理剪切变形的数值模拟和工程估计具有一定的实用价值。     

海相沉积软黏土的弹塑性本构模型研究

在深入探讨海相沉积原状软黏土压缩、变形等力学特性和详细分析加载屈服面随荷载情况变化的基础上,确认了海相沉积原状软黏土的强度、变形特性与结构屈服应力密切相关。即当固结压力小于结构屈服应力时,其力学特性与超固结重塑土的力学特性类似;当固结压力大于结构屈服应力时,其力学特性与正常固结重塑土的力学特性类似。为描述海相沉积原状软黏土的上述力学特性,将姚仰平等提出的超固结重塑土本构模型引入到海相沉积软黏土弹塑性本构模型的构建中。在本构模型构建过程中,考虑了海相沉积原状软黏土具有的抗拉强度及其演化规律,软黏土强度包线的特点及其进一步修正的表达式,使模型更符合海相原状软黏土的强度、变形特性。最后,将3种不同海相沉积软黏土固结排水剪切试验得到的应力-应变-体变曲线与模型预测结果进行对比。比较结果显示,本文提出的弹塑性本构模型能很好地描述海相沉积原状软黏土的剪缩硬化、剪胀软化以及变形的应力水平依存性等力学特性。     

基于Plesha本构的岩石节理多层结构模型研究

在Dong等工作的基础上[1],研究了基于Plesha本构的岩石节理多层结构模型。模型将节理面分解为不同层次细观结构面,每层结构面代表一级粗糙度,力学响应发生在最底层（基本面）,上层结构面受力性状由其下层结构面平均化得到,依次进行得到节理宏观力学性状。节理破坏由下而上分层次进行,基本面破坏后,其上层结构面转化为基本面,基本面采用Plesha节理模型。模型能模拟弹性变形、滑动变形、磨损、剪断、压碎、分离等作用机理。定义的接触因子与绕流因子能考虑节理接触面积及粗糙度变化对渗流的影响。通过将粗糙度定义为等效起伏角,可以模拟节理循环剪切性状。采用ABAQUS的用户子程序UEL进行了模型验证与参数分析。     

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

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

千枚岩的剪切流变特性研究

蠕变性和时效特征是岩石的固有力学属性。岩石的流变特性是解释和分析一些地质构造现象、岩体变形破坏以及岩体长期稳定性预测的重要依据。通过对千枚岩沿片理面的剪切流变试验研究,建立千枚岩流变本构模型;并在此基础上,采用伯格斯体模型对该岩石力学参数进行了拟合分析,获取千枚岩沿片理面的长期强度及剪切蠕变本构方程,为该类岩体流变分析与计算提供了理论依据。     

自然与饱水条件下异性结构面的剪切特性研究

利用室内剪切试验装置,开展不同正应力、自然与饱水条件下岩体异性结构面的剪切试验,分析异性结构面的剪切变形和强度特性,利用JRC-JCS抗剪强度力学模型将试验结果与计算结果进行对比分析。结果表明,（1）异性结构面的剪切变形曲线属于弹塑性残余变形类型。相比自然状态试件,饱水试件的弹性区减小,较早进入峰值区或塑性区;（2）异性结构面的剪切刚度均随法向力的增大而增大,但同一正应力下饱水试件的剪切刚度要小于自然状态的剪切刚度,在低法向力下尤为显著;（3）水对岩体异性结构面剪切强度参数有不同程度的弱化,尤其对内摩擦角 的劣化明显;（4）只要计算参数选取得当,并考虑实际岩体的受力及水等情况,JRC-JCS计算模型可以为实际工程岩体的结构面强度提供较为准确的参数。     

重复剪切作用下泥质充填结构面剪切力学特性试验研究

结构面的力学特性直接关系工程岩体的变形稳定,而充填物对结构面的力学特性存在较大的影响.为了准确地量化充填前后结构面剪切性能的变化,以三峡库区典型库岸边坡节理岩体为研究对象,采用劈裂法制备了人工节理,按照充填度100％进行考虑,设计进行了充填前、后结构面的重复剪切试验.结果表明:(1)基于重复剪切过程中结构面抗剪性能劣化...     

结构面的剪切蠕变特性研究

岩体结构面的蠕变特性研究对解决岩石力学实际问题具有很重要的意义。本文利用水泥砂浆制作的规则齿形结构面试样进行了剪切蠕变室内试验,针对不同角度结构面试件在不同法向应力条件下的剪切蠕变特性进行了研究,并根据蠕变试验结果,提出了结构面剪切蠕变经验公式;对结构面的长期强度特性进行研究,并选取Burgers模型来描述蠕变的剪切变形特性,文章最后讨论了Burgers模型对于描述结构面剪切蠕变特性的参数特征。     

岩体锚固效应及锚杆的解析本构模型研究

利用Mohr-Coulomb、Hoek-Brown以及Duncan-Chang理论分别分析了块状和碎块状岩体锚固后的物理效应,通过分析认为,岩体锚固后可以有效地提高岩体的凝聚力和软弱结构面的抗剪强度,增强岩体的弹性模量,改善岩体的力学性质。通过对拉拔试验测试结果分析,在一定的假设条件下,推导了锚杆与注浆体或岩体耦合情况下的解析本构方程,并对其进行了参数分析,认为提高拉拔力、增大锚杆直径和锚固段长度可以有效地改善其锚固效果,并提出了锚固临界值的概念。在前人工作的基础上,利用提出的耦合解析本构模型建立了非耦合状态下锚杆的解析本构模型,通过计算认为,该解析模型是合理的。     

碎屑砂岩三轴压缩下强度和变形特性试验研究

基于碎屑岩组织结构疏松、含水率较高、物理力学性能较差、呈孔隙式胶结接触等特点,对碎屑砂岩首先开展了物理特性试验分析,认为其微、细观结构复杂、内部破坏严重,矿物成分为石英、长石、绢云母等,化学成分以SiO2为主,属微透水、小孔隙率砂岩,且渗水化学侵蚀并不显著。其次,开展了静水压力、单轴压缩和三轴压缩试验,研究了碎屑砂岩的强度和变形破坏特性。最后,初步探索了物理特性与强度变形特性的关系。结果表明,静水压力为2.6 MPa时,岩样内部微缺陷压密完成;单轴压缩曲线呈明显6阶段特征,峰值应力达0.98 MPa,属脆-延性破坏;三轴压缩条件下,岩样呈压缩为主的延性扩容破坏,轴向压缩和环向体积扩容达6%和4%;曲线无明显破坏荷载,呈现非线性、塑性硬化、存在屈服平台和体积由压缩向扩容过渡等特性。且体积扩容破损应力与屈服应力基本相同,扩容转折点随围压增加而增大,围压可增强岩样抵抗变形破坏的能力。试验结果旨在为岩石工程稳定分析及本构模型构建提供可靠的依据。     

基于改进遗传算法的岩体结构面幂函数模型参数辨识

沿结构面的剪切滑移是工程岩体主要的破坏模式之一,选择恰当的岩体结构面本构模型是分析岩体结构面破坏的重要一环。当前常用的结构面本构模型有指数型模型、曲线型模型和幂函数型模型,其模型参数都是根据工程实际经验或室内试验所得,不能很好的模拟实际变形情况。指数型模型得出的模拟值偏大,曲线型模型得出的模拟值偏小,幂函数型模型在较低的法向应力下能较好的模拟剪切变形过程,但在较高的法向应力下得出的模拟值偏小。为此,通过改进的遗传算法对岩体结构面幂函数模型参数进行辨识,得出由"等效参数"构成的幂函数模型。实例验证结果说明,通过改进的遗传算法辨识出的由"等效参数"构成的幂函数模型能够很好的模拟实际变形情况,能够很好的适应于不同的工程实际应用。      

一种基于弹性能释放率的岩石新型统计损伤本构模型

为了揭示岩石变形的破坏机理以及岩石材料产生损伤的本质原因,文章对岩石材料变形规律和力学特性进行分析后,再以损伤变量作为影响岩石变形和力学性能变化的内变量,采用能量原理、有效应力原理和统计损伤理论构建了一种基于弹性能释放率的新型岩石统计损伤本构模型。该损伤模型进一步完善了岩石损伤本构模型的理论体系,弥补了传统损伤模型无法合理解释引发岩石破坏原因的不足。利用岩石试验数据对损伤模型的参数进行确定,并将损伤演化模型代入弹性能-应变模型中,分析在加载过程中岩石弹性能变化的规律。结果表明:模型曲线与试验曲线在峰前变形阶段几乎重合,说明损伤模型可以很好地反映岩石的变形特性;在初始加载阶段,岩石的损伤变量随着轴向应变的增大而增大,说明在荷载作用下,岩石内部裂隙逐渐发展发育,使得岩石材料的损伤逐步积累;在围压达到10 MPa以上时,损伤-应变曲线基本重合。同时,在初始加载时刻,损伤-应变曲线增长率急剧上升,大约在岩石应变为0.01%时,损伤-应变曲线趋于平稳变化状态;且由于岩石在峰值应力点附近损伤迅速累积,进而使得损伤变量在数值上快速增大到1,这说明了围压的增大使得岩石破坏极限得到显著的提升。     

多参数岩体结构面优势分组方法研究

在岩体力学和水力学分析中,搞清岩体结构面的分布规律是一项重要的基础工作。当岩体结构面的产状一致,而其他性质不一致时,岩体结构面的力学和水力学性质是不同的,而传统的方法仅根据岩体结构面的产状进行分组,这势必存在缺陷。提出了一种基于变尺度混沌优化算法的多参数结构面数据的优势组划分方法,通过度量结构面之间的相似性建立目标函数,运用变尺度混沌优化算法搜索目标函数的最优解来确定聚类中心,可同时考虑结构面的多个参数,把具有相似性质的结构面归并成组。通过计算机模拟的结构面数据验证了此算法的正确性,最后将此算法应用于实际工程中测量的多参数结构面数据的优势组划分,得到了清晰可靠的分组结果。     

岩体裂隙模拟的扩展有限元法应用研究

岩体中普遍存在着断层﹑节理和裂隙等结构面,这些结构面的存在和发展对岩体的整体强度﹑变形及稳定性有极大的影响。因此,研究岩体中原生结构面的萌生﹑发展以及贯通演化过程对评估岩体工程安全性和可靠性具有非常重要的理论与现实意义。扩展有限元法（XFEM）作为一种求解不连续问题的有效数值方法,模拟裂隙时独立于网格,因此,在模拟岩体裂隙扩展﹑水力劈裂等方面具有独特优势。针对扩展有限元法的基本理论及其在岩体裂隙扩展模拟中的应用展开了研究,建立了扩展有限元法求解岩体裂隙摩擦接触、岩体裂隙破坏等问题的数值模型,并将计算模型应用于岩质边坡稳定性分析和重力坝坝基断裂破坏等工程问题。     

泥岩巷道围岩弹塑性参数反演分析与长期稳定性预测

以实测的围岩蠕变变形资料为基础,基于现场监测所获得的蠕变本构模型,对某一矿区泥岩体的力学参数进行了弹塑性反演分析,得到的巷道围岩体基本力学参数分别为弹性模量E=2.0 GPa,凝聚力c=1.31 MPa,内摩擦角? =24º。待反演参数水平均值的极差结果显示,通过极差的大小可以判断岩土力学参数的敏感性,对蠕变的影响而言,凝聚力最敏感,其次是内摩擦角,弹性模量再次之。在此基础上,利用反演所得的围岩基本力学参数进行了正演,计算得到了各测点处对应的蠕变位移增量,与实测值相比吻合较好,表明反分析所获得的岩体材料参数综合反映了巷道围岩的力学特性。最后,对巷道围岩变形大小及塑性区范围进行了预测。     

基于人工蜂群算法的岩体结构面多参数优势分组研究

岩体结构面优势分组是研究岩体力学性质与水力特性的基础,通常的分析方法是只根据产状进行划分。鉴于结构面其他特征对岩体力学性质的重要影响,考虑结构面倾向、倾角、迹长、张开度、表面形态5个特征参数,提出了基于人工蜂群算法的岩体结构面多参数优势分组方法。以样本总体离差平方和为目标函数,建立岩体结构面多参数优势分组的数学模型,应用人工蜂群优化算法求解,以目标函数值最小时的解作为聚类中心,并自动确定分组边界。对人工生成的结构面数据的计算结果验证了该方法的正确性,该方法的求解精度是令人满意的。最后,将该方法应用于怒江松塔水电站坝址区岩体结构面多参数优势组的划分,得到了较为合理的分组结果,进一步验证了此方法具有较高的运行效率与工程实用性。     

采动煤岩渗透率演化模型及数值模拟

为描述采动煤岩渗透率演化过程,引入强度退化指数,基于Hoek-Brown强度准则,建立了考虑围压影响的煤岩应变软化力学本构模型。给出了体积应变和渗透率的关系方程,结合应变软化模型建立了采动煤岩渗透率演化模型,并在FLAC下予以实现。通过数值模拟研究了不同围压下圆柱岩样的峰后应变软化力学行为和某煤矿工作面开采过程中煤岩的渗透率演化过程,结果表明：（1）该模型能较好地反映围压对煤岩峰后应变软化行为的影响;（2）随着工作面推进,越来越多的煤岩单元破坏,渗透率也不断增长,逐渐成为瓦斯等流体运移的主要通道。（3）模型能再现采动煤岩渗透率演化的动态过程,从而为煤与瓦斯共采、煤层瓦斯抽放和瓦斯灾害防治提供指导。     

基于细观变形理论的盐岩塑性本构模型研究

唯象学本构模型不能很好解决诸如夹杂含量、温度、应变速率等对盐岩力学性质的影响,更难以解释盐岩变形机制。盐岩为石岩晶体组成,其变形机制主要由多晶结构所控制,故基于固体位错理论研究方法建立的盐岩塑性本构模型更能反映盐岩的变形机制。研究表明,盐岩的塑性-蠕变交互作用机制是（亚）晶粒内部位错的滑移与（亚）晶界及其干涉面内位错的攀移运动之间的耦合。基于此,可确定亚晶（或晶粒）平均尺寸与流动应力之间的关系、（亚）晶内的位错平均密度;建立微观参量（位错、亚晶直径、亚晶界宽度等）演化模式;根据Orowan定律建立盐岩微观-宏观变形联系,从而导出盐岩塑性本构方程。导出的本构方程体现了盐岩塑性-蠕变变形的物理机制,相对于传统的塑性本构方程具有更好的物理意义。