基于颗粒离散元法的盘形刀具破碎性能研究.docxVIP

基于颗粒离散元法的盘形刀具破碎性能研究 隧道全段（tm）广泛应用于地下空间隧道开挖工程中。在sdm波形刀破岩的情况下，挖掘过程的核心内容是波形刀破岩的勘探。因此，对波形刀破岩的研究非常重要。盘形滚刀破岩过程中, 地质条件是极其复杂的, 如围压、节理、岩溶等, 这些地质因素对盘形滚刀破岩必然有较大的影响。Chen L H等通过实验研究得出了当围压达到一个临界值时, 裂纹的扩展就会受到限制, 且围压的改变对刀具的贯入力影响不大的结论。张魁等利用UDEC软件模拟了在不同围压的作用下盘形滚刀破岩过程, 研究了围压和最优刀间距之间的关系。Ma H S等利用RFPA软件研究了不同围压对盘形滚刀破岩的影响, 得到了围压的改变会改变裂纹扩展方向、有效裂纹长度等结论。孙金山等利用PFC2D软件研究了节理强度等对TBM滚刀破岩的影响, 得到了较好的结论。Bejari H等利用UDEC软件研究了节理间距和节理倾角对单刀破岩的影响, 研究表明节理间距增大会减小刀具的切入率, 节理倾角在25°~30°时最有利于刀具破岩。刘红岩等通过实验研究了节理对岩石的破坏模式、单轴抗压强度、弹性模量的影响。Ghazvinian A等通过实验研究表明岩石的强度随围压增大而增大, 岩石强度随节理倾角的增大先增大后减小, 其中节理倾角为45°时岩石强度最小。这些研究成果对于研究不同围压、节理特征下的盘形滚刀破岩具有很大的参考价值。在此基础上, 本文借助于PFC2D离散元软件对在不同围压、节理特征下的盘形滚刀破岩进行了数值模拟, 研究围压和节理特征对盘形滚刀破岩的影响, 为数值模拟技术在刀具破岩研究中的应用提供参考。 1 滚滚轴破岩的数值模型 1.1 宏观力学特性 锦屏二级水电站地层中主要以大理岩和花岗岩等硬岩为主, 这类硬岩石单轴抗压强度普遍在100MPa左右, 属于较高强度岩石地层, 选择该地层的某一硬岩为研究对象, 其岩石试样的宏观力学参数见表1。 岩石试样的细观力学参数一般是通过模拟单轴压缩、巴西劈裂和直剪等数值实验进行反复标定, 不断去匹配岩石试样的宏观力学参数得到的。设定岩石试样数值模型的最大和最小颗粒半径之比为1.66, 颗粒的最小半径取0.8 mm;颗粒密度ρ=3 375 kg/m3, 由于岩石试样存在空隙率, 所以颗粒密度大于岩石试样的宏观密度;选用接触黏结模型来模拟颗粒之间的接触方式。单轴压缩数值试验岩石试样尺寸为100 mm×50 mm, 如图1 (a) 所示, 上下两道墙体为加载压盘模拟压缩实验时用的加载板, 赋予上下墙一定的速度来给试样进行加载;巴西劈裂数值试验岩石试样的直径为50 mm, 如图1 (b) 所示, 通过给两侧墙一个缓慢的速度来给圆盘加载, 直至圆盘发生破裂;直剪数值试验尺寸为100 mm×100 mm, 如图1 (c) 所示, 通过伺服加压程序给上下墙体施加不同恒定的正应力, 下剪切盒固定, 上剪切盒以恒定的速度进行剪切。经过上述不断标定选取合适的细观参数可得到岩石试样在PFC2D中的细观参数见表2。 1.2 刚性墙体剪切试验 在PFC2D软件中, 节理是在先生成岩石试样以后再通过jset命令设置的。通过定位设置指定某一平面两侧的颗粒间的接触处, 将此接触处两侧的颗粒接触强度以及摩擦因数重新设定, 即可在岩石试样中形成节理。为了研究节理力学特性, 通过PFC2D内置FISH语言建立直剪数值试验模型。节理直剪试验模型如图2所示, 由刚性墙体组成上下2个剪切盒, 其中位于下方的剪切盒固定, 位于上方的剪切盒可以按照设定的速度向左移动, 对试样进行剪切。模型试样的上下墙通过伺服加压程序, 可以设定不同的正应力, 如图2所示的上下墙上的箭头。试验时记录剪切过程中的剪切应力和剪切位移。模型所剪岩石试样尺寸为100 mm×100 mm, 在其中间位置水平方向设置一连通节理。 岩石试样采用的细观参数见表2, 通过jset命令生成节理后, 节理处颗粒的细观参数被重新赋值, 其中节理处颗粒的黏结强度取50 k Pa, 摩擦因数取0.15。在不同的法向应力作用下进行数值直剪数值试验, 得到剪切应力与剪切位移的关系, 如图3所示。 提取数值直剪试验中得到的峰值剪切应力, 可得其与法向应力的关系, 当法向应力为0.1, 1.0, 1.5, 2.0, 2.5 MPa时, 峰值剪切应力分别为9.05, 9.62, 10.2, 10.8, 11.4 MPa, 并根据上面的数据拟合得其关系如图4所示。从图中可知该节理的黏聚力和内摩擦角分别为8.45 MPa和50°。 1.3 盘形1盘形小程序 盘形滚刀在破岩过程中, 一方面由于切向摩擦力的作用, 盘形滚刀绕自身中心自转, 另一方面盘形滚刀受法向推力侵入岩体, 沿法向方向不断贯入。而就裂纹的产生和扩展而言, 主要是