水射流破碎非均质岩体的极限.docxVIP

水射流破碎非均质岩体的极限 在水和辐射的作用下，岩体破坏的本构关系或失败的第一个问题被机械界广泛接受，并尚未得到解决。目前的研究方法基本上是对实验现象的观察和实验结果的抽象和简化。考虑到问题的不同含义，对问题的简化和假设存在差异，最终结论也不同。 笔者在水射流钻孔的试验中发现,非均质岩体的破坏形式主要是在拉应力作用下的脆性破坏,具体表现为沿着低强度微元产生的径向裂纹、锥状裂纹和横向裂纹及其扩展.岩体在射流的冲击下产生的拉应力和剪应力分别超过了非均质岩体低强度单元的抗拉和抗剪的极限强度,使得岩体中强度低的单元首先破坏,在岩体中形成裂隙.裂隙形成和汇交后,水射流将进入裂隙空间,在水楔作用下,裂隙尖端产生拉应力集中,使裂隙迅速发展和扩大,致使岩体破碎. 1 单元抗压强度的计算 目前关于水射流破岩的理论有:① 把水射流的冲击作用简化为静压力使岩体的某一部位的拉应力或剪应力超过其极限而导致岩体破坏的静态弹性理论;② 在水射流冲击作用下产生的应力波是导致岩体破坏主要原因的应力波破碎理论;③ 认为射流中空穴破裂产生的局部巨大压力导致岩体破坏的气蚀破碎理论;④ 由于裂纹扩展联通导致岩体破坏的裂纹扩展破碎理论. 分析得知,水射流的作用只能使岩体中强度较弱的单元受到的拉应力或剪应力超过其抗拉或抗剪强度极限而直接破坏;而高强度单元则是由于岩体裂隙的发展、连通,以剥离的方式脱离岩体而失效. 对于直接破坏的单元,破坏准则为最大拉应力: |σ3|≥T0, 最大剪应力:τn≥c+σntanφ, 其中,T0为单元的抗拉强度;c为内聚力;σn为受力面上的正应力;φ为内摩擦角. 不能直接破坏的高强度的失效准则,可由岩石的非均质性,通过愈渗理论得到.设非均质岩体的微元强度服从Weibull分布,则岩体单元抗压强度分布函数为 式中,m为岩体的非均质参数;Rc为单元抗压强度;Rc0为岩体单元的平均抗压强度. 设在某一级载荷作用下已破坏的微元体数目为Nf,定义统计损伤变量为已破坏的微元体数目与总微元体数目N之比,即D=Nf/N,这样在任意区间[Rc,Rc+dRc]内已破坏的微元数目为NP(Rc)dRc, 当加载到某一水平Rc时,已破坏的微元数目为 得到用损伤变量D表示的非均质岩体损伤变量演化方程为 根据愈渗理论,当该岩体中破坏单元比例超过一定的值后,则认为这部分单元是游离的,即被剥离岩体而失效.设Mf为愈渗理论值,则得出水射流作用下,岩体的破坏准则为D≥Mf.将式(1)代入得 当水射流冲击压力大于Rc时,岩体中将有比例超过Mf的单元被破坏,由愈渗理论知岩体将被破坏.否则,岩体不能被破坏.因此,Rc为水射流冲击破碎非均质岩体的门槛压力pm,即pm=Rc0[-ln(1-Mf)]1/m. 2 射流冲击煤体的力学模型 笔者用C语言编制了水射流破岩的三维有限元计算模型,并选用空间8节点立方体单元作为有限元计算的基本单元.根据煤层受力的特点,计算中选用垂直对称模型.数值计算模型的大小为200 mm×80 mm×40 mm,比实际钻孔范围小,可定性地分析水射流的钻孔.图1为该模型的网格剖分,该模型共有31 920个单元和35 322个结点.图2为边界条件简化. 由于水射流喷嘴直径很小,其射流对岩体的作用区域很小,但射流压力的变化却很大,因此,在水射流作用区域内,有限元计算网格的划分应该很密.以尽可能反映水射流对岩体的作用机理.水射流冲击作用以外的区域,则选用较疏的剖分网格,以适应计算机的计算能力. 数值模拟中,射流速度在射流截面服从高斯分布,其计算公式为 式中,re为射流截面特征半径,re=0.114x,x为射流冲击距离;r0为喷嘴半径;v0为射流出口速度. 由此得出射流冲击煤体的作用力f=ρqv(1-cosβ), 其中,ρ为流体密度;q为射流体积流量;β为射流方向变化角. 计算中使用的水射流密度为1 000 kg/m3,黏度为0.001 Pa·s, 压力为70 MPa,喷嘴个数为9个,喷嘴直径为0.8 mm.煤岩体弹性模量为1 874.1 MPa,泊松比为0.3,密度为1 400 kg/m3,单轴抗压强度为15.0 MPa,内摩擦角33.7°,孔隙率4.00%,渗透系数14.445 m/s,非均质性系数7.0. 图3为煤岩体在9喷嘴组合喷头水射流钻孔过程中,冲击方向(x方向)对称面的位移等值线,中间部分为已经破碎的部分.由于煤岩体的非均质性和分布喷嘴中心冲击压力高的特点,可以看到形成的凹凸不平的冲击面.并使得煤体强度低的细胞元先破碎,产生一系列大小不等的裂隙.由于位移等值线的特性,煤体产生的微小裂隙被样条曲线平滑掉了.但仍可以看出水射流的冲击破碎和钻孔的形成,首先在冲击面的中心范围内形成裂隙,然后裂纹沿低强度细胞元发展、延伸,逐渐汇聚,发生大块煤岩体的剥落,形成破碎坑.随着射流冲