地质体具有非均匀、非连续、流固耦合等特性。渗流在地质体的失稳破坏中具有两重性，一方面流体可以弱化地质体中节理、断层等结构面从而加剧地质灾害的发生；另一方面，基于渗流的水力压裂技术成功应用于油气开采、防治冲击地压等领域从而提高工程效益。渗流是地质体变形破裂的主要因素，但是流体与地质体的相互作用研究尚有不足，因此建立渗流作用下地质体变形破裂的高精度计算方法是十分有必要的。
本文首先提出了一种任意单元刚体运动及弹性变形一体化计算方法，用于模拟地质体的连续变形阶段；其次，发展了一套岩体破裂的自适应计算方法，实现地质体内部裂缝扩展路径的精确模拟；然后建立了三维岩体流场-破裂耦合计算方法，刻画流体驱动下的岩体破裂过程；最后在实际工程中应用，分析水力压裂作用对煤层坚硬顶板垮落的影响。主要研究内容和结论如下：
（1） 从能量角度出发，选取刚体运动的平动位移、角位移，以及表征弹性变形的应变作为独立的广义变量，基于拉格朗日方程，提出任意单元刚体运动及弹性变形一体化计算方法。通过分解单元，将节点边界转化为广义变量的边界问题，完成了针对该方法的几何边界条件施加。在此基础上发展其高阶形式，实现地质体连续变形阶段应力场的精确模拟。单摆运动、重力场求解以及一维杆弹性波传播等数值算例验证该方法在计算转动问题、静态问题、动态问题方面的正确性和可靠性，最后通过悬臂梁案例验证该方法高阶形式在计算连续变形方面的精确性。
（2） 在上述计算方法的基础上，引入应变强度分布准则，选取破裂度作为裂缝起裂依据，依据最大拉应力-摩尔库仑复合强度准则判断裂缝扩展方向，从而建立岩体破裂的自适应计算方法，并且基于BP神经网络模型发展一种新的网格剖分算法，实现地质体内部裂缝扩展路径的精确、高效计算。通过三点弯曲梁、含孔洞四点弯曲梁和三轴压缩数值案例验证了该方法在拉伸以及剪切破坏形式下的裂缝扩展路径的准确性。
（3） 基于弹簧元模型，发展三维渗流弹簧元计算方法，在此基础上建立三维岩体流场-破裂耦合数值计算方法，解决流体驱动下岩体的破裂问题。该方法采用连续-非连续单元法计算岩体的变形和破裂，采用渗流弹簧元计算流体在裂隙通道中的流动以及流体在岩体内部的孔隙渗流。通过四个理论解和一个实验现象的对比分析，发现该方法可以精确捕捉裂隙通道和岩体内部的流体压力以及精准刻画流体作用下的岩体裂缝形态，最后一个圆盘水力压裂模拟说明该方法不仅可以描述裂缝面内和岩体基质中的流体压力分布，而且同时可以跟踪岩体内部裂缝的萌生、扩展以及贯通过程。
（4） 根据本文发展的计算方法，开展布尔台煤矿4-2煤层顶板水力压裂数值模拟。首先根据4-2煤层建立巴西圆盘试样，数值模拟了巴西圆盘劈裂过程，其裂缝扩展形态以及加载点力和位移关系均与实验解吻合良好；其次从流体压力、垮落步距和竖向应力三个方面数值分析了水力压裂技术对4-2煤层上方坚硬顶板的影响，说明水力压裂技术能够有效弱化顶板岩石强度、缩短垮落步距，达到坚硬煤层顶板下井下工作面的安全高效回采。