非连续变形分析方法的改进及应用1）

首先,采用位移收敛判据对非连续变形的岩石断裂分析进行改进并程序实现. 其次,将非连续变形分析方法和裂隙网络渗流模型结合起来,建立了渗流应力耦合分析模型,研究了裂隙岩体变形对渗流的影响以及渗流与应力耦合作用下裂隙岩体的变形破坏特征. 研究发现,主干裂隙控制着渗流场分布,考虑耦合作用时的水头值比不考虑耦合作用时大,说明考虑耦合后坡体更容易变形失稳. 改进前后裂纹扩展形态基本不变,表明程序计算的收敛判据是有效的.     

非连续变形分析方法的改进及应用

首先,采用位移收敛判据对非连续变形的岩石断裂分析进行改进并程序实现.其次,将非连续变形分析方法和裂隙网络渗流模型结合起来,建立了渗流应力耦合分析模型,研究了裂隙岩体变形对渗流的影响以及渗流与应力耦合作用下裂隙岩体的变形破坏特征.研究发现,主干裂隙控制着渗流场分布,考虑耦合作用时的水头值比不考虑耦合作用时大,说明考虑耦合后坡体更容易变形失稳.改进前后裂纹扩展形态基本不变,表明程序计算的收敛判据是有效的.     

岩体水力劈裂的应力-渗流-损伤耦合模型研究

利用无单元法追踪裂纹扩展的优势研究模拟岩体水力劈裂的数值分析模型。定义损伤变量以描述岩体的损伤程度对岩体渗透性和强度的影响,建立了岩体水力劈裂的应力-渗流-损伤耦合分析模型,并编制了基于无单元法的计算程序。综合考虑应力、渗流及损伤之间的相互作用影响,分析了具有初始裂纹的岩体平面应力模型的裂纹扩展过程,指出考虑渗流作用时的裂纹扩展角大于不考虑渗流作用时的裂纹扩展角。同时,裂纹内水压力和渗流的作用对于裂纹的扩展方向和扩展长度具有较大的影响。     

脆性岩石破裂过程损伤与渗流耦合数值模型研究

大量的实验结果表明，脆性岩石的渗透性不是一个常量，而是应力和应力诱发损伤破裂的函数．建立了一个描述非均匀岩石渗流-应力-损伤耦合数学模型(FSD Model)，开发出岩石破裂过程渗流-应力-损伤耦合分析计算系统(F-RFPA^2D)．在该系统中，单元的力学、水力学性质根据统计分布而变化，以体现材料的随机不均质性，材料在开裂破坏过程中流体压力传递通过单元渗流-损伤耦合迭代来实现．该系统能够对岩石试件在孔隙水压力和双轴荷载作用下裂纹的萌生、扩展过程中渗透率演化规律及其渗流-应力耦合机制进行模拟分析．最后，给出两个算例：算例1模拟载荷作用下岩石应力应变-渗透率全过程．模拟结果表明，非均匀性对岩石的应力峰值强度、峰值前后其渗透性演化规律及其破裂机制影响十分明显，模拟结果和实验结果较为一致；算例2模拟孔隙水压力作用下岩石拉伸断裂过程，通过和物理实验对比验证，验证了模型计算结果的可靠性。     

水力压裂中裂纹扩展的数值模拟

水力压裂是在高压粘滞流体或清水作用下地层内裂缝起裂与扩展的过程。由于包含岩石断裂和流-固耦合等复杂问题,对该过程的数值模拟具有相当大的挑战性。本文建立基于有限元与离散元混合方法的裂纹模型,模拟岩石裂纹扩展,实现了连续向非连续的转化;建立双重介质流动模型,裂隙流作为孔隙渗流的压力边界,孔隙渗流反作用裂隙的压力求解,处理了流体在基岩与人工裂缝中的协调流动;将裂纹模型与流体流动模式进行结合,建立断裂-应力-渗流耦合形式的力学模型,进一步分析了水力压裂的基本过程,综合多种数值计算方法,编写程序,在验证岩体裂纹模型与双重介质流动模型有效性的基础上,对压裂过程进行复现,将模拟结果与文献结果进行了对比,并讨论了所构建模型的优缺点。     

岩石破裂与渗流耦合过程细观力学模型

从岩石细观非均匀性的特点出发,应用弹性损伤力学、B iot渗流力学理论,对岩石破裂过程渗流-应力-损伤耦合模型进行了描述.算例表明,该模型能较好的模拟出岩石材料在水压和载荷作用下裂纹扩展路径不规则发展及稳态、瞬态渗流过程,裂纹扩展对渗流路径和流动过程起到明显的控制作用.     

不同耦合介质爆破裂纹动态扩展特性研究

不耦合装药结构能有效改善爆破效果。为探究不同耦合介质对爆破效果的影响，利用透射式数字激光焦散线实验系统，对有机玻璃板(polymethyl methacrylate, PMMA)试件开展水、空气、沙土3种耦合介质的不耦合爆破试验，研究不同耦合介质爆后裂纹形貌，计算裂纹分形维数，定量评价岩石破坏程度。分析裂纹扩展速度与动态应力强度因子变化规律，对比不同耦合介质下爆生裂纹的动态力学行为。结合数值模拟进行辅助分析，解释不同耦合介质对数值模拟中爆破损伤范围、孔壁压力时间曲线和岩石模型的总能量曲线的影响规律。结果表明：水耦合爆后裂纹分形维数最大，炮孔周围裂纹数量众多，主裂纹扩展长度最大，爆破效果最佳。沙土耦合爆破需要更多的能量使裂尖起裂，其动态应力强度因子最大。模拟得到的损伤云图与实验结果一致，水耦合损伤范围最大，空气耦合次之，沙土耦合损伤范围集中在炮孔周围。而孔壁压力时程曲线、能量时程曲线则验证爆破实验结果，表明水耦合爆破能更好地利用爆炸能量，提高试件的破坏效果。     

拉伸荷载下分支裂隙破坏机理研究

为研究拉伸荷载下分支裂隙对破坏模式的影响,保持主裂隙参数不变,改变分支裂隙倾角和长度,利用扩展有限元方法模拟了弯折裂隙的动态扩展,总结了分支裂隙参数变化对破坏模式的影响。利用ABAQUS中的轮廓积分计算了分支裂隙尖端应力强度因子,并根据最大周向应力准则计算起裂角。结果表明:拉伸荷载下分支裂隙出现三种破坏模式;分支裂隙倾角和长度均对破坏模式有一定的影响。I型应力强度因子与分支裂隙倾角关系曲线呈斜"S"型,相应II型应力强度因子曲线呈上凸型;由于分支裂隙存在非尖端破坏,利用裂隙尖端应力强度因子判断开裂应结合相应的破坏模式。     

考虑闭合效应和三维应力约束的表面裂纹扩展模拟

假定承受Ⅰ型常幅载荷下的表面裂纹在扩展中的形状保持为半椭圆,利用Newman半椭圆表面裂纹应力强度因子公式计算应力强度因子。提出了等效厚度的概念,利用穿透直裂纹的研究结果,考虑表面裂纹扩展中塑性致闭和三维应力约束效应。基于Elber模型建立了三维表面裂纹扩展模型。数值模拟了表面裂纹扩展过程,研究了裂纹形状变化及规律,计算了裂纹扩展寿命。将计算结果与有关试验结果进行了对比,吻合较好。     

不同节理角度和地应力条件下岩石双孔爆破的数值模拟

在采矿工程、地下交通工程和水利水电工程等岩体开挖工程中,爆破法仍是一种主要的破岩方法. 实际天然岩体中存在的裂隙、节理等不连续面和所处的地应力场环境,都会对爆破荷载传播过程和最终的破碎效果产生重要影响. 本文把岩石爆破视为爆炸应力波动态作用和爆生气体压力准静态作用两个独立的先后作用过程,同时考虑初始应力场的静态作用,建立岩石爆破的力学模型. 基于渗流方程描述爆生气体在爆生裂纹中的传播过程,进而基于流固耦合理论实现爆生气体准静态压力对裂纹尖端的力学作用. 爆炸应力波主要在孔周引起压碎区和径向微裂纹区,随后爆生气体压力楔入径向裂纹尖端,促使裂纹进一步扩展,形成径向主裂纹. 数值模型可以再现孔周压碎区、径向微裂纹区、径向主裂纹区萌生、扩展的完整演化过程. 针对不同节理角度和地应力条件下岩石双孔爆破过程的数值模拟结果表明,初始地应力场的压应力作用不利于爆生裂纹的萌生与扩展,但节理的存在对裂纹的扩展具有明显的导向和促进作用,有利于爆生裂纹沿节理面方向的扩展.      

冲击载荷作用下岩体破坏规律的数值流形方法模拟研究

为了更好地模拟岩体这种典型的连续与非连续介质共存的复杂结构体在冲击载荷作用下的破坏规律,在近年来新出现的数值流形方法基本理论的基础上,利用其算法程序对岩体中存在的不连续面如节理、裂隙等对岩体破裂效果的影响进行了数值模拟分析,计算结果充分显示了数值流形方法在模拟这类问题中的准确性和有效性。     

动载荷作用下裂隙岩体的止裂机理分析

为深层次了解裂隙岩体在动载荷作用下的动态断裂特性及止裂机理,采用TWSRC（tunnel with single radial crack）构型进行中低速冲击实验,选择砂岩作为原材料制作裂隙岩体试样,以落锤冲击试验装置与裂纹扩展计实验系统对裂纹的动态起裂、扩展及止裂过程进行全过程监测,重点研究动态破裂过程的破裂行为及止裂现象。使用有限差分法程序进行数值模拟,验证冲击实验结果的科学性与准确性。研究发现:裂隙岩体的动态断裂过程是由起裂加速-高速扩展-缓慢减速-止裂-再次起裂加速-再次高速扩展等多次循环的过程构成,且止裂区间尺寸为微秒量级;裂隙岩体止裂位置的穿晶断裂比例远小于初始起裂点,青砂岩动态断裂过程的穿晶断裂比例稍大于黑砂岩;裂隙岩体中止裂点再次起裂所需的能量,远小于预制裂纹初始起裂所需要的能量。     

单轴荷载下饱水岩石静态和动态抗压强度的细观力学分析

由于单轴荷载下饱水岩石的动态力学特性与静态力学特性存在很大差异,从宏观上进行力学分析 存在局限性。根据岩石受压全应力应变曲线的细观机制,分析了静态及动态单轴荷载条件下孔隙水影响饱水 岩石裂纹扩展的情况。在静态单轴压缩条件下,初始裂隙受压使自由水产生孔隙水压力,自由水对翼裂纹有 向外挤压的应力,促进裂纹扩展。在动态单轴压缩条件下,自由水会产生粘结力,抑制裂纹扩展。根据翼裂纹 受压扩展原理,推导出饱水单轴条件下动态抗压强度、静态抗压强度的计算公式,在相同断裂韧度下,饱水岩 石静态抗压强度风干岩石静态抗压强度饱水岩石动态抗压强度。对自然风干和饱水砂岩进行单轴静态、 动态压缩实验,结果与理论模型的结果相符。     

不同围压作用下非均匀岩石水压致裂过程的数值模拟

从岩石细观非均匀性的特点出发，提出一个描述非均匀材料渗流和破裂相互作用的数值模型。在这个数值模型中，单元的力学、水力学性质根据统计分布而变化，以体现材料的随机不均质性，材料在开裂破坏过程中流体压力传递通过单元渗流，损伤耦合迭代来实现。算例表明，该模型能较好地模拟出岩石类材料在水力压裂作用下，微结构非均匀分布和不同围压比对破裂模式、失稳压力的影响，非均匀性导致试件的开裂压力、失稳压力明显不同，裂纹扩展路径不规则发展，模拟结果和实验结果较为一致。     

岩石材料动态破坏的近场动力学方法数值模拟

裂纹的聚集、扩展、分叉是岩石等脆性材料破坏失效的关键因素，本文在验证了近场动力学方法在研究岩石类材料裂纹动态扩展方面的有效性之后，采用近场动力学方法分别对冲击载荷作用下含有双裂纹岩石材料和单轴压缩作用下含单斜裂纹的岩石材料进行数值模拟．结果表明，对于双垂直裂纹，其裂纹扩展路径大致与预制裂纹成70°夹角；对于单裂纹，裂纹的扩展路径随裂纹倾角的变化而变化，最终导致构件的整体破坏．数值模拟结果表明近场动力学方法可以很好地模拟岩石等脆性材料的裂纹扩展直至破坏的过程，反映裂纹扩展的物理机理；其作为一种新的基于非局部理论的数值方法，在地下岩体工程方面及页岩气的开采方面会有很好的发展前景．     

边坡岩体结构的三维失稳形式及稳定性分析研究

由两组结构面控制的四面体块体是岩体中最基本的结构体, 其稳定性取决于块体所受荷载及结构面上的摩阻力。正确判定边界面上摩阻力分布是评价块体稳定性的关键。本文提出的块体稳定性分析适用于岩体所受荷载为三维非共点力系。在不同外力作用下块体可产生平动、转动、翘扭及倾倒等各种三维失稳组合形式, 从而得到更为合理的块体稳定性评价。     

渗流对岩体高边坡位移的影响

根据渗流场和位移场耦合机理 ,结合某大型船闸工程 ,采用有限元法进行计算分析和对比 ,结果表明船闸在开挖过程中 ,渗流作用力对开挖边坡的稳定性影响很大。并指出对高边坡进行稳定分析时 ,正确考虑渗流作用力的影响十分重要。