应用蒙特卡洛法确定节理岩体的连通率和综合抗剪强度指标

本文提出了一种确定节理岩体连通率和岩体综合抗剪强度指标的新方法。该方法对传统的连通率的概念作了扩充,将连通率的确定归结为在岩体结构面网络图中寻找节理面和岩桥组合破坏路径的问题。在连通率的确定过程中,考虑了节理和岩桥的组合破坏机制。并应用动态规划原理,基本解决了岩体中存在多组节理面情况下的连通率确定问题,所得成果在漫湾电站左岸边坡稳定分析研究中获得了成功的应用。     

三维结构面连通率的随机模拟计算

通过对裂隙岩体中大量随机分布的结构面调查资料的统计分析，得出结构面分布的统计规律，运用随机模拟的方法，建立了裂隙岩体的三维结构面网络。通过三维结构面网络图像，在大量实际工程岩体破坏面特征调查的基础上，提出了三维结构面连通率的概念，并应用岩桥破坏机理，选取通过破坏面的结构面，采用三维曲面拟合的方法，建立了三维结构面连通率的计算方法。     

岩桥破坏的简化模型及在节理岩体模拟网络中的应用

提出一个简化的岩桥破坏 /结构面扩展模型 ,该模型符合现有试验观测成果 ,在力学机理上基于Lajtai对岩桥剪切破坏的研究 ,能方便地计算实际问题中各种结构面组合情况下的抗剪力。利用先进的优化方法 ,该模型用于在结构面模拟网络中搜索抗剪力最小的结构面-岩桥组合 ,从而研究裂隙岩体的结构面连通特性和综合抗剪强度。算例和分析表明 ,在不同的法向应力作用下 ,岩体沿某一组结构面受剪时 ,破坏路径上的连通率为常数 ,这时可以使用连通率求取岩体的综合抗剪强度指标。在剪切方向与某一组结构面的平均延展方向不相同时 ,破坏路径上的连通率随法向应力的改变而改变 ,这种情况下岩体的综合抗剪强度只能通过模拟试验直接求得。     

节理岩体三维综合抗剪强度数值模拟研究

采用不确定性理论建立的三维节理网络模拟技术是计算机技术在岩体力学领域应尉的一个重要方面。该技术利用野外不连续面现场地质调查所获得的原始数据，建立结构面几何参数分布的概率统计模型，进而应用蒙特卡罗模拟方法构筑不连续面在三维空间的组合形态，借助计算机生成三维网络模型。在建立三维网络模型的基础上，通过搜索不连续面和岩桥组合破坏的临界路径来确定连通率并提供三维综合抗剪强度指标，在搜索过程中应用了具有鲜明生物背景的遗传算法。     

节理岩体结构面连通性研究及其应用

岩体结构面通过几何学和力学两方面影响岩体稳定性。因此，可将结构面连通性研究分解为2类问题：（1）反映空间方位效应及岩体水力学特性的几何学连通性；（2）反映岩体强度的力学连通性：从两方面能更全面的反映岩体结构面的工程性质。首先采用单测线法对节理岩体结构面进行统计测量，对获取的结构面的产状和迹长进行取样偏差校正和蒙特卡洛模拟，从而形成结构面网络模型。然后，对网络模型进行分析与计算，得到结构面在几何学意义上的连通性：并由网络模型寻求真正的优势控制性结构面。最后，以此类优势结构面为基础，与其他随机结构面进行耦合分析，结合岩桥破坏机制，运用动态规划方法搜索抗剪力最小路径，得到更符合工程实际的结构面在力学意义上的连通率。     

岩体结构面连通率估计的概率模型及其工程应用

H-H连通率估计的概率模型是计算岩体结构面连通率的有效方法之一，但是由于它要求结构面均匀分布而且锁固段长度不能超过概率模型方法的预测极限，因此在实际应用中受到了较大限制。为了克服这一局限性，详细分析了产生这一现象的原因，建立了相应的等效模型。在此基础上，提出了广义H—H连通率估计的概率模型，扩大了H—H连通率估计的概率模型的应用范围，并通过对某水电工程侧裂结构面连通率的预测，进一步证实了该模型的正确性和可靠性。     

基于平均迹长修正新方法的岩体三维连通率及抗剪强度参数计算

 岩体结构面连通率是反映结构面延伸程度和连通状况的一个重要参数,是确定裂隙岩体抗剪强度指标的关键和分析边坡稳定性的前提。对三维连通率的求解方法及其在估算抗剪强度参数中的应用进行探讨,在考虑结构面迹长和迹线与统计窗边线的交角均遵循一般概率分布类型的基础上,从概率统计角度推导出窗口法中结构面平均迹长的估算公式。利用基于平均迹长修正公式的三维结构面网络模拟方法对贵阳市鱼简河水利工程导流洞进口的岩体结构面进行模拟,在其三维可视化模型的基础上求得模拟区内不同截面的三维连通率,并利用三维连通率结果估算出岩体的综合抗剪强度参数。计算结果表明：该模拟区内岩体总体连通性较好,各向异性明显,70°～80°∠20°～40°范围内的截面为危险截面,最大连通率为0.914 3,最小黏聚力和内摩擦角分别为1.135 MPa和25.7°,工程建设中应予以重视。     

用投影法求算岩体结构面三维连通率

岩体结构面三维连通率是一个内涵空间岩体稳定性评价的重要参数,它完全不同于常规的与水文地质或岩体水力学相关的连通率。由于岩体中结构面发育的随机性,采用一般的方法很难得到岩体结构面在空间的三维连通率。介绍的岩体结构面三维连通率是基于现场工程岩体结构面三维网络数值模拟的基础之上,采用理想截面投影法来求取工程岩体关键部位、特定方向上的三维连通率,并通过计算机编程来实现计算,直接应用于工程中。     

寻优方法在岩体结构面连通特性和综合抗剪强度研究中的应用

首先，通过一简单的算例，详细介绍了遗传算法和动态规划在复杂的岩体结构面网络中搜索抗剪力最小的结构面组合的方法，进而将两种方法进行对比，并通过实例说明两种方法的寻优结果是一致的。最后，讨论了剪切边界条件的选取对计算结果的影响。实例分析显示，固定剪切路径端点情况下的寻优计算的结果与不固定剪切路径端点情况下的寻优计算的结果是不同的，前者得出的连通率较小，综合抗剪强度较大。     

岩体结构面连通率研究进展及应用

岩体结构面连通率的研究一直是岩石力学研究领域的难点之一。一些学者已经相继提出了多种确定连通率的方法，但许多方法可能引起的误差较大，实用性也较差。广义H—H连通率公式和基于网络模拟的计算方法是近年发展起来的比较成熟且实用的方法。文章评述了各种方法的适用性和有效性，并采用基于结构面网络模拟的“带宽投影法”进行了实例分析，最后讨论了岩体结构面连通率的研究趋向。     

基于3D打印技术的模拟柱状节理岩体试样制备方法

以规则正六棱柱型柱状节理岩体为对象,介绍了基于3D打印技术的模拟柱状节理岩体试样制备方法,即在SolidWork中建立柱状节理网络3D打印数字模型,以光敏树脂为打印材料,采用3D打印技术打印具有不同倾角的圆柱形柱状节理网络模型,并以此为模具,采用白水泥浆为模型材料进行浇筑,经过拆模、黏结等程序制备具有不同倾角的模拟柱状节理岩体试件;再通过对试件进行单轴抗压强度试验,分析模拟柱状节理岩体试件的强度、变形和破坏特征,初步验证了上述基于3D打印技术的制样方法能够较好考虑柱状节理岩体的结构特性,反映柱状节理对岩体强度的弱化效应,所制备的柱状节理岩体试件基本满足柱状节理岩体室内试验的要求。     

应用单元劈裂法模拟三维内嵌裂纹扩展

 对于岩石内嵌裂纹很难通过物理试验进行研究,在应用有限单元法模拟内嵌裂纹扩展时,三维节理面建模、网格划分以及节理扩展后的网格处理一直是难点问题。为简单、高效地分析含有内嵌型不连续面的节理岩体的破坏模式,利用三维单元劈裂法来解决这一问题。单元劈裂法是利用四面体单元的几何性质来构造的一种特殊的四结点节理单元。当有裂纹穿过四面体单元时,位于裂纹两侧的结点可形成3或4个接触点对,由此推导出劈裂单元的刚度矩阵。采用该方法模拟内嵌裂纹扩展时,可以不用考虑节理岩体的几何完整性,而直接对整个岩体进行网格划分。因而一些单元则被内嵌裂纹相切,对于与裂纹相切的单元采用劈裂单元刚度矩阵,而其余单元则采用常规单元刚度矩阵。由此免去了单独设置三维节理单元的麻烦以及为设置节理单元而不断改变网格划分方案的麻烦。这种优势在裂纹数目较多时尤为突出,非常适合节理岩体计算。通过对内嵌椭圆裂纹扩展模拟并与相关物理试验做对比,得出该方法能够再现三维内嵌裂纹扩展的基本特征,说明该方法是有效的。     

三维条件下应用遗传算法与Monte-Carlo法确定节理岩体的综合抗剪强度

使用遗传算法并充分利用二维研究成果，提出了在结构面三维模拟网络中搜寻结构面与岩桥的最小抗剪力组合的断层扫描方法，从而可以应用Monte-Carlo法的原理，在计算机上多次模拟岩体的剪切试验，最终确定节理岩体的综合抗剪强度。该方法反映了结构面在三维岩体空间的几何分布特征，抓住了控制岩体力学性质的主要因素，并通过使用二维的结构面扩展模型反映了结构面与岩桥的破坏机理及结构面与岩桥的相互作用机理。实例分析的计算结果与反演结果相符，说明了该方法的合理性。     

节理岩体渗流的无网格与有限元耦合方法

节理岩体的渗透破坏、节理岩体中石油渗流、地下工程的防渗设计等无不与渗流计算有关,而节理岩体渗流介质的非均质性和各向异性特性、边界条件和几何形状的复杂性、初始条件的不确定性等因素制约着解析法在渗流计算中的应用。应用无网格法来求解节理岩体的渗流问题,建立节理岩体渗流的无网格与有限元耦合方法。通过变分原理详细推导无网格与有限元耦合方法求解节理岩体渗流的离散方程及相关计算公式,并利用耦合方法对边界条件进行处理。最后,给出的数值算例说明该方法的正确性和有效性。     

基于广义Hoek-Brown非线性破坏准则的节理岩体地基承载力研究

 针对节理岩体非线性破坏特征,发展基于Hoek-Brown破坏准则的临界滑动场理论进行节理岩体地基承载力的计算。首先将Hoek-Brown准则的剪切强度逐点等效到Mohr-Coulomb强度线上,求得每点的瞬时内摩擦角和瞬时黏聚力;在此基础上,改进基于Mohr-Coulomb准则的临界滑动场理论,建立新的迭代算法,将Hoek-Brown强度准则与临界滑动场理论结合起来,求解被动土压力和地基承载力。同时,综合分析GSI和mi值对地基承载力的影响。算例对比分析结果表明,该方法能迅速准确地确定节理岩体地基最危险破坏滑面并得到相应的地基极限承载力值。     

考虑岩体节理分布特性的各向异性强度准则研究

 针对岩体节理分布的特性,以Mohr-Coulomb和Rankine准则为基础,引入结构面网络模拟计算所得的岩体连通率和强度参数,研究岩体在各向异性强度准则下的理论方法。在此基础上,对岩体各向异性强度准则及其数值实现方法进行研究,并编制相应的计算程序对坝体–岩基系统的深层滑动面进行动力分析。分析结果表明,岩体在各向异性强度准则下的破坏面不一定与最大剪应力面重合;综合考虑Mohr-Coulomb和Rankine准则可改进Mohr-Coulomb准则对岩体抗拉强度过高估计的缺点,更符合工程实际;所建立的各向异性强度准则能真实地模拟各向抗力与响应的相互关系所反映的节理裂隙岩体的破坏特征。     

柱状节理岩体压缩破坏过程模拟及机制分析

 将柱状节理岩体整体考虑为一种广义宏观复合材料是一种行之有效的研究方法,为此建立复合型多弱面软化模型,采用VC++语言成功开发并嵌入FLAC3D程序。以白鹤滩水电站的地质条件为研究背景,对柱状节理岩体进行单轴压缩和三轴压缩的模拟计算,分析不同节理面倾角和倾向的试块破坏过程,得出各种情况下破裂面分布和破坏形态,并总结具有普遍意义的破坏规律。通过三轴试验还得到在处于合适的地应力条件时,柱状节理岩体整体强度及节理面之间的咬合力均较高的结论,与现场勘测结果一致。同时,研究压缩破坏过程中外载和变形的规律、塑性破坏区的产生、破裂区的分布和发展以及能量的耗散等,从多方面对柱状节理岩体的破坏机制进行探讨,以便更好了解该类岩体的力学特性。     

非连续变形分析(DDA)中断续节理扩展的模拟方法研究和试验验证

在非连续变形分析方法中,可用虚节理表示岩石破坏过程中潜在的破坏路径,因此断续节理中节理面开裂、扩展、贯通过程认为是虚节理(岩桥)向实节理转化的过程,这一过程伴随着虚节理(岩桥)力学性质参数的弱化,且弱化程度受到虚节理向实节理转化程度的影响。因此,利用表征虚节理向实节理转化程度的虚节理连通率对虚节理力学性质参数的弱化规律进行研究,得出力学性质参数弱化函数;在此基础上对Jennings强度准则和最大抗拉强度准则进行处理,提出应用于非连续变形分析计算的断续节理强度形式;最后,将虚节理力学性质参数弱化规律和处理后的断续节理强度表达式运用于非连续变形分析计算程序中,实现对断续节理扩展过程的模拟算法。并通过剪切试验计算结果与室内试验结果的比较,验证该算法的正确性。     

拉剪条件下节理岩体中锚杆的力学作用分析

在分析穿过节理的锚杆与岩体相互作用的机理后，采用线弹性断裂力学的方法，分析在拉剪综合作用下锚杆对裂纹尖端应力强度因子产生的贡献，并揭示了各种应力作用情况以及锚杆与节理面之间不同的夹角下锚杆的作用规律。计算分析结果表明，锚杆的作用使节理端部的应力强度因子发生转换，从而明显降低了对岩体破坏产生主要作用的应力强度因子。这是锚杆能够加固节理岩体的重要原因。     

断裂破坏强度模型在大型油库山体边坡支护中的应用

岩体中存在大量孔洞、节理和裂隙，它们作为岩体的初始损伤而存在其中。正是由于这些节理裂隙的初始损伤、损伤演化、扩展断裂等缺陷导致岩体产生了各种形式的破坏。根据节理裂隙的扩展断裂机理，深入研究了节理扩展及扩展相互作用对岩体强度特性的影响，建立了节理岩体断裂破坏强度模型，并将该模型用于指导广东大湾油库山体高边坡支护设计，取得了令人满意的效果，创造了显著的经济效益。