平台中心角对巴西劈裂过程及破坏模式的影响

平台中心角对巴西劈裂过程及破坏模式的影响

岩石的抗应力是岩石力学和工程的重要参数之一。目前,测定岩石抗拉强度的方法有:直接拉伸试验,巴西劈裂法,正方形板对轴压裂试验,对轴压模拉伸试验等。由于巴西劈裂容易操作和对试样要求较低,所以通常采用巴西劈裂法测定岩石的抗拉强度。1978年,国际岩石力学学会(ISRM)推荐该方法为测定岩石抗拉强度的方法。国内,各种岩石力学试验标准也推荐使用巴西劈裂法测定岩石抗拉强度。按照规程要求,试验机对巴西圆盘施加径向荷载使其破坏。假定其为线性荷载,基于平面应力弹性方法进行求解,可得到圆盘各处的应力解析解。从而,通过计算可得到岩石的抗拉强度。然而,这一经典的测试方法还有很多不足。事实上,巴西劈裂和直接拉伸所测的抗拉强度之间有较大的差别。进行巴西劈裂试验时,加载处会因为应力集中首先脆断,巴西劈裂的中心起裂条件不能得到保证。为了改善加载处的应力集中现象,可以加工2个平台面作为加载面,这就是平台巴西劈裂试验。那么,就需要确定平台中心角的大小和抗拉强度计算式的修正系数。尤明庆等利用ANSYS软件,结合格里菲斯强度准则对平台圆盘中心的受力状态和岩石抗拉强度进行了分析。喻勇勇基于Mohr强度理论,提出了一个利用平台圆盘测试岩石抗拉强度的计算公式。尤明庆等还通过试验得出抗拉强度随平台中心角2α增大的规律,并推荐平台中心角为20°~30°。于庆磊等利用RFPA程序,观察并分析了平台中心角对试样破坏模式和应力分布的影响。已有的研究大都利用有限元方法和某种强度理论进行分析,然后确定平台中心角及建立抗拉强度计算公式。然而岩石材料自身的非均质、不确定性往往被忽略而且岩石内部破坏过程也没有被重视。本文作者基于颗粒流数值模拟程序建立了平台巴西劈裂试验的细观离散元模型,研究了劈裂过程中裂纹发展规律。通过对不同平台中心角下的巴西平台试验的分析,确定了合理的平台中心角的值和该平台角下岩石抗拉强度的计算式。1岩石模型和测试过程的p1.1基于fish函数的土体模型PFC作为一种离散元数值分析方法,通过分析圆形颗粒的运动及其相互作用来求解岩石力学问题。使用BPM(黏结颗粒模型)时,模型内的颗粒随机生成,其粒径服从可以服从均匀分布或正态分布;通过内置的FISH函数,黏结模型的强度参数服从正态分布,从而模型具有一定的非均质性和随机性。这就能很好的表示岩土体的非均值性和参数的离散性。PFC通过定义简单的颗粒接触模型,就可以表示材料复杂的破坏机制,避免了使用岩土体难以获得的本构关系。1.1.1最大作用力pfc接触刚度模型提供了接触力和位移的弹性变化关系,可表示为:式中:Fin和△Fis分别为法向和切向力;Kn和Ks分别是法向刚度,切向刚度;Un和ΔUs分别为法向重叠量和切向位移增量。滑动模型通过定义摩擦系数得到最大摩擦力,在模拟中根据判断最大摩擦力Fsmax和静摩擦力Fis的大小关系来判断是否会滑动。最大摩擦力的计算公式如下:式中:μ为摩擦因数。PFC提供了2种基本的黏结模型:平行黏结模型和接触黏结模型。因为平行黏结模型可以很好的表示黏结断裂后,材料细观结构刚度的降低,所以对于岩石材料而言更为适合。平行黏结模型可以假想为胶结聚合物,颗粒在接触点处胶结在一起。同时黏结物自身受力时,允许产生一定的变形,当接触力超过黏结强度时,黏结破坏。1.1.2实体接触力的形成过程PFC是按照时步迭代算法不断计算,直至得到问题的解。计算过程的每一步包含2个基本部分(如图1所示)。每一计算步开始时,首先更新接触关系和实体(球和墙)位置的数据。然后对接触应用力和位移定律,从而得到实体间的接触力。接着,就可以根据已得到的接触力,应用牛顿第二定律决定每个颗粒的运动,如此不断循环。为了得到稳定,可靠的计算结果,计算时步必须足够的小。1.2pbc试验过程1.2.1初始应力条件的设置本次模型的建立可分为如下4步。(2)通过半径扩大的方法得到颗粒集合体初始的各向同性应力条件。(3)去除漂浮粒子。球的配位数小于4的颗粒定义为漂浮粒子,并予去除。(4)设置平行黏结模型。这样平台巴西圆盘的岩石模型最终生成。1.2.2宏观参数标定PFC的输入变量并不是人们所熟悉的宏观参数,而必须经过多次标定后才能找到的细观参数。所以使用该数值方法时,首先要假定其输入参数,然后对模拟结果和真实材料的力学性质进行比较;不断调整输入参数,比较其结果,最终两者结果相一致的一组参数就是所需的模型参数。本次标定过程使用了某矿山的大理岩试样,所取的宏观参数具有一定的代表性。标定过程选择了表1所示的宏观参数,进行了巴西劈裂和单轴压缩虚拟试验。标定过程后,颗粒半径在0.16~0.26mm范围之间,模型主要微观参数如表2所示。1.2.3墙与颗粒间的运动在圆盘的2个平台面上分别设置一个墙,设墙的速度为1×10-2m/s(在物理试验中速度很大;但对于PFC来说,约4×10-10m/步已经足够小)进行加载。为便于问题的研究,墙的摩擦因数设为0,即加载过程中墙与颗粒之间无摩擦作用。当某点黏结破坏时,FISH程序会在两球球心连线的法线方向上用一条线段标志破裂时所产生的裂纹。2结果分析2.1抗裂特性分析本文通过观察不同平台中心角圆盘内裂纹的扩展过程,发现了岩石破坏的一般过程。图2所示为平台中心角为15°时,所做的平台巴西圆盘劈裂过程示意图。在图2中,黑色线段表示为裂纹。由于岩石内随机分布的缺陷,当其受力时,在微孔或微裂隙末端易造成应力集中,从而在低应力状态下岩石内部就会有微裂纹的产生和扩展。所以,在加载的初始阶段,产生了随机分布的裂纹(如图2(a)和(b)所示),这和基于强度理论分析的结论相异。可以看到:端部应力集中效应依然存在,端部附近分布的裂纹较其他地方要集中。当平台中心角大于20°时,这种效应已不明显(如图3所示)。随着荷载的不断增加,裂纹的出现变得有序。裂纹集中出现在已有的裂纹周围并且裂纹之间相互搭接,沿径向不断扩展。最终,产生的裂纹沿径向贯通,试样破坏。实际上,岩石的破裂过程是破坏单元从无序分布到集中有序的自组织过程。在加载过程的初始阶段,荷载较小,产生于岩石内的裂纹带有很大的随机性。所以,本文并不赞同在岩石均匀性假设的前提下,通过强度理论来判断岩石这种特殊材料的断裂起始位置。根据尤明庆等的分析,平台中心角2α≥20°时,平台圆盘满足中心起裂的条件。然而在这种情况下,试样受压断裂时,初始的裂纹依然是随机分布,这一点是由具体一块岩石中随机分布的缺陷决定的。随着荷载的增加,岩石的破裂过程逐渐变得有序。通过应力分析和强度理论,可以确定相对应力最大的位置在圆盘的受压直径上,也就是圆盘会对径劈裂为相等的2部分。试验结果表明:不断增加的裂纹沿着圆盘的受压直径,集中分布在已有的裂纹周围。随着裂纹的不断增加,裂纹之间不断相互作用逐渐贯通,最终形成一个破坏面。观察到的破坏的过程和理论上的分析还有一定的差异,其主要原因是岩石的非均质性以及破裂过程的非线性。2.2拉裂裂试验现象整个劈裂过程中,裂纹的扩展和应力变化具有很大的相关性(如图4所示)。在加载初始阶段应力较小,裂纹很少而且扩展缓慢。在临近破坏荷载时,裂纹发展迅速,弹性应变能迅速释放,驱动裂纹迅速扩展、贯通。峰值应力后产生的大量裂纹主要是断裂面拉开时的拉伸裂纹。这个阶段拉裂纹变化最为明显,试件由体积压缩转为扩容,随后宏观断裂面形成。平台中心角会影响到试件内的应力分布,进而对整个劈裂过程产生影响。如图5所示,平台中心角越大,劈裂试验过程中裂纹产生得越多、增长得越快,特别当平台中心角大于25°时,裂纹数量增加得更加迅速。裂纹数量随着应变增加而增加,可以作为岩石内损伤的演化和裂纹的扩展程度的表征。这意味着中心角越大,岩石损伤越大,破坏模式变得更加复杂。2.3平台中心角抗拉能力当平台中心角不大于20°时,可以看到巴西圆盘可以沿直径劈裂为两半。当平台中心角大于20°时,巴西盘的断面变得复杂,而且平台圆盘中心裂纹扩展方向和受力方向有一定的夹角(如图6所示)。随着平台中心角增大,拉应力降低,而压、拉应力比增大,圆盘破裂方式从单纯的拉断到应力组合的复杂断裂方式。当平台中心角等于26.7°时,这种组合破坏方式表现的更加明显。所以,当平台中心角大于25°以后,试验测量所得到的名义抗拉强度能否反映岩石的抗拉性能值得怀疑。然而,平台中心角越大,端部的集中应力影响就越小和加工可操作性和精度越好。目前已有的研究成果推荐平台中心角一般为20°~30°。综合以上分析,本文推荐平台中心角为20°。2.4模型建立和修正确定平台中心角为20°后,就可以对岩石抗拉强度进行测定。对于α=20°的平台圆盘巴西劈裂其名义抗拉强度σt可以表示为:式中:k为选定的平台中心角为20°的抗拉强度计算式的修正系数。P为试件劈裂破坏时的最大荷载(N);d为岩石圆盘的直径(m);t为圆盘试件的厚度(m)。st′为直接拉伸时的岩石抗拉强度。所以k(28)st/st′。数值模拟实验操作方便、可重复性强,因此本文基于PFC2D进行了虚拟试验仿真模拟了平台巴西劈裂和直接拉伸试验。采用如上标定过程中的参数,经过平台巴西劈裂数值实验,得到了平台中心角为20°时的名义抗拉强度st′=29.5MPa;采用同样的参数,建立模型(如图7所示,试件长60mm,直径50mm),对两侧球体赋恒定速度为0.01m/s(约4×10-10m/步)进行加载。最终试件破坏,测得直接拉伸的抗拉强度为σt=27.48MPa。最后,算得修正系数,k=0.9315。文献所给结果为0.9644,2种方法得到的修正系数差别并不大。那么,平台中心角20°的平台巴西劈裂抗拉强度计算公式为:3不同平台中心角的巴西使用情况(1)对平台圆盘劈裂过程进行观察,发现了岩石内裂纹的形成、发展的一般过程。由于岩石的非均质性,岩石破坏单元的产生和发展是一个从无序到有序的自组织过程。(2)平台中心角越大,裂纹的数量越多和发展速度更快,