岩石三维表面裂纹扩展机理数值模拟研究

岩石破坏的本质原因是由于内部裂隙的萌生、扩展与贯通过程。从三维的角度出发,采用细观损伤数值模拟方法,模拟单轴压缩下含预制三维表面裂纹的岩石试样的破坏过程。数值模拟得到了表面裂隙内部扩展、贯通过程,动态再现翼型裂纹、壳体裂纹的形态,探讨三维裂纹内部的受力机制,推测可能发生的断裂类型,进一步探讨三维裂纹扩展规律。研究结果表明：①反翼型裂纹并不一定萌生于预制裂纹端部,是由于翼型裂纹扩展后应力释放后的拉应力引起;②壳体裂纹的萌生与扩展阶段是由Ⅲ型加载断裂主导,而翼型裂纹扩展至一定长度之后停滞不前;③除了反翼型裂纹之外,还新发现了一种由壳体裂纹萌生出的次生裂纹,这种裂纹的扩展引起试样整体失稳崩溃;④岩石Ⅲ型加载（反平面剪切）难以获得Ⅲ型断裂破坏,壳体裂纹是由于Ⅲ型加载下的拉应力引起,实际上属于Ⅰ型与Ⅱ型复合裂纹;⑤非均匀性对岩石表面裂纹扩展影响很大,相对均匀岩石中难以出现曲线翼型裂纹或反翼裂纹。研究结果对于岩石三维裂隙扩展机理的物理力学实验与理论分析都具有参考意义。     

等效岩体三维随机节理网络模型构建及其在两河口水电站中的应用

节理的存在对水电高陡岩质边坡的力学性质有重要影响,如何构建反映三维节理分布特征的等效岩体计算模型,是分析与评价岩体力学特性的关键。本文基于损伤力学和统计强度理论,在三维岩石破裂过程分析（RFPA3D）软件的基础上,提出了一种计算等效岩体三维随机节理网络模型的新方法。首先,基于Baecher模型和Monte-Carlo方法,在RFPA3D软件中实现了三维随机离散节理网络（Discrete fracture network,DFN）模型的重构。然后,利用RFPA3D软件内嵌DFN模型,赋予节理和岩石不同的力学参数,构建了工程尺度等效岩体三维随机节理网络模型,实现了三维随机节理岩体破裂过程、变形和强度等力学性质的分析。最后,以两河口水电站左岸边坡坝址区下游节理岩体为研究对象,验证了三维随机DFN模型的准确性,开展了研究区内节理岩体尺寸效应研究,并获得了研究区内岩体的表征单元体（Representative elementary volume,REV）和等效力学参数。该研究成果为等效岩体力学行为分析提供一种新方法。      

应力波形对岩石爆生裂纹扩展机制影响的数值模拟

运用RFPA3D动力分析软件模拟了冲击动力作用下含预制裂纹岩石的裂纹扩展过程,探究了应力波峰值、能量、上升及下降速率对岩石裂纹扩展过程的影响。研究表明动载下岩石裂纹扩展形态受应力波上升速率影响,应力波上升速率越快,孔周边岩石越破碎;应力波能量影响裂纹扩展长度,能量越大裂纹扩展越长,而相同能量条件下,应力波上升速率越小,裂纹扩展距离越远,但孔边破碎程度越弱;上升速率和应力波上升沿能量共同影响着炮孔粉碎区半径。数值模拟结果很好地揭示了不同应力波峰值、能量与上升/下降速率对岩石的破碎机制,在实际爆破作业中可以通过水炮泥封口或者采用空气柱间隔装药结构来延长应力波作用时间,以达到扩大爆破影响范围的目的,而通过选取合适类型与配比的炸药来提升应力波上升速率从而增强孔边破碎效果。      

软弱夹层对深部地下洞室围岩损伤模式的影响

软弱夹层等不良结构面对深部地下洞室的围岩稳定有着重要影响,应用有限元数值分析系统RFPA对软弱夹层影响下的深部地下洞室围岩损伤演化过程进行了模拟分析,重点研究了软弱夹层的展布位置、地应力大小对洞室围岩损伤区模式的影响.分析表明,有软弱夹层时,地下洞室开挖后围岩应力分布恶化,围岩应力扰动范围增大、相应的损伤区范围增大;地应力越大时,侧压系数的大小更多地主导着洞室围岩的损伤模式.研究成果可望对深部地下洞室选址及施工设计有借鉴意义.     

刀具破岩机制研究现状及发展趋势

在分析大量文献的基础上,从理论模型、物理试验和数值模拟三个方面,对近些年国内外刀具破岩规律研究进展情况进行归纳和总结,并对各种研究方法的适用性、优缺点进行了评述,侧重论述了数值模拟研究进展,探讨了当前研究热点问题,包括节理岩体中的破岩规律、开裂准则、动力破岩和三维数值模拟研究等。总体来看,有关刀具破岩机理的研究取得了长足进展,尤其是数值计算从不同层次上揭示了刀具破岩机制,包括滚刀作用下岩石裂纹扩展规律、滚刀与岩石的相互作用力学模型、材料力学参数、滚刀作用方式、围岩地应力以及刀具几何参数等因素对破岩效率的影响,数值计算方法极大地促进了刀具破岩规律的深入研究;最后总结了当前刀具破岩受力模型、动态加载与三维模拟等面临的挑战问题。     

高性能计算技术在岩石力学课程教学中的应用

岩石力学是一门实践性很强的课程,工程实例、实验模型和理论模型的结合将直接影响课程的教学效果。介绍了把基于高性能计算技术的软、硬件平台应用于岩石力学教学的情况,拓宽了传统实验教学中的实验对象,实现了对工程实例、实验模型和理论模型的精细数值表征;通过鼓励学生自发设计模型进行数值模拟分析,激发了学生的思维能力和创造性,有利于学生理解和掌握教学内容。在高性能计算机日益普及的形势下,对国内高校岩石力学教学提供了有益的参考。     

水压致裂过程的三维数值模拟研究

基于RFPA数值分析方法和并行计算技术,建立了反映岩石细观损伤演化过程的三维渗流–应力–损伤耦合模型。对具有120万单元的方形岩石材料模型,进行了4组不同应力状态下水压致裂过程的三维大规模科学计算分析。计算结果分析表明：起裂压力与失稳压力并不重合,起始裂纹均为张性,裂纹扩展形式、表面平整度、走向、扩展失稳过程以及裂纹的空间分布形态受应力状态影响。当竖直方向为最大主应力方向时裂纹呈空间竖片分布,当水平应力差较大时裂纹表面形态平整,失稳到来较快;当竖直方向为最小主应力方向时裂纹的空间分布呈水平片状;不等的主应力情况下裂纹总是分布在最小主应力面内;当三向主应力相等时,裂纹起裂位置和扩展方向具有竞争趋势,空间分布不具规律,裂缝分支较多。数值模拟结果与物理实验结果有着较好的吻合,该研究对水压致裂工程设计有一定参考价值。     

岩石破裂失稳声发射监测频段信息识别研究

首先对声发射信号进行分频处理,分析岩石损伤破裂演化过程中声发射信号频段分布规律。根据岩石断裂模式与声发射频率存在一定的对应关系,证明了岩石破裂监测中声发射前兆优势频段的存在,进一步构建前兆优势频段的计算方法,确定声发射监测岩石破裂的最优短临预报参数。通过开展花岗岩圆形隧洞模型水平卸荷试验,分析整个过程的声发射信号,寻找前兆优势频段,以此验证方法的可靠性。研究结果表明,水平方向卸荷瞬间,隧洞两壁发生劈裂破坏,频段(31.25~62 k Hz)的小波能量占比达73%以上。临近最终破裂,频段(7.8125~15.625 k Hz)变能系数CD6?由前一刻的0.5突增到15,具有典型的灾变响应特征,可以确定该频段为声发射前兆优势频段。选择CD6??作为宏观破裂的短临预报参数,再结合临近最终破裂前声发射主频在频段(90~105 k Hz)和(13~20 k Hz)所出现的响应突现规律,可以为预测花岗岩圆形隧洞模型的失稳破坏提供预警参量。     

主应力对洞室围岩失稳破坏行为的影响研究

 针对复杂应力环境三向应力状态下地下洞室围岩破坏条件和破裂机制,采用统计损伤理论和数值模拟方法,建立三维非均匀性地质模型,考虑应力的三维效应,引入强度折减法,一方面在保持模型边界条件的同时实现洞室围岩的逐步破坏,另一方面以此定量评价不同应力场中洞室安全稳定状况,探讨不同侧压力系数和轴向应力条件下洞室围岩破坏模式,探究中间主应力对洞室稳定性的影响以及深部岩体分区破裂化现象产生条件和破裂规律等。结果表明,侧压力系数影响洞室围岩初始破裂形成部位和发展趋势;不同的轴向应力使得洞室围岩破裂区域和范围显著不同,在不同的侧压力系数条件下,轴向应力影响洞室稳定的规律存在差异;不同方向中间主应力对洞室围岩安全稳定状况的影响是不同的;当洞室轴线方向与最大水平应力方向平行时,较大的轴向应力会使洞室围岩产生分区破裂化现象,围岩破坏的区域也是拉应变集中的区域等。这些结果对进一步揭示地下洞室围岩非线性变形破坏行为,评价岩土工程安全稳定性,采取合理的支护措施等均具有重要意义。     

应力波作用下三点弯曲梁破坏过程数值模拟

岩石类准脆性材料三点弯曲梁冲击荷载作用下的破坏受应力波加载速率以及应力波峰值的影响。使用新发展的并行RFPA2D-Dynamic数值模拟工具对三点弯曲梁冲击荷载作用下的破坏过程进行了数值模拟,研究了不同加载速率和加载峰值的破坏规律。数值模拟结果与物理实验结果具有较好的一致性。在低加载速率下,三点弯曲梁发生拉伸破坏,且只出现一条裂纹;在较大加载速率时出现多条平行拉伸裂纹;在大加载速率时出现拉伸、剪切复合裂纹。应力波峰值较小时形成单一的拉伸裂纹,而随着加载峰值的提高,形成多条平行拉伸裂纹。