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

平台中心角对巴西劈裂过程及破坏模式的影响 岩石的抗应力是岩石力学和工程的重要参数之一。目前,测定岩石抗拉强度的方法有:直接拉伸试验,巴西劈裂法,正方形板对轴压裂试验,对轴压模拉伸试验等。由于巴西劈裂容易操作和对试样要求较低,所以通常采用巴西劈裂法测定岩石的抗拉强度。1978年,国际岩石力学学会(ISRM)推荐该方法为测定岩石抗拉强度的方法。国内,各种岩石力学试验标准也推荐使用巴西劈裂法测定岩石抗拉强度。按照规程要求,试验机对巴西圆盘施加径向荷载使其破坏。假定其为线性荷载,基于平面应力弹性方法进行求解,可得到圆盘各处的应力解析解。从而,通过计算可得到岩石的抗拉强度。然而,这一经典的测试方法还有很多不足。事实上,巴西劈裂和直接拉伸所测的抗拉强度之间有较大的差别。进行巴西劈裂试验时,加载处会因为应力集中首先脆断,巴西劈裂的中心起裂条件不能得到保证。为了改善加载处的应力集中现象,可以加工2个平台面作为加载面,这就是平台巴西劈裂试验。那么,就需要确定平台中心角的大小和抗拉强度计算式的修正系数。尤明庆等利用ANSYS软件,结合格里菲斯强度准则对平台圆盘中心的受力状态和岩石抗拉强度进行了分析。喻勇勇基于Mohr强度理论,提出了一个利用平台圆盘测试岩石抗拉强度的计算公式。尤明庆等还通过试验得出抗拉强度随平台中心角2α增大的规律,并推荐平台中心角为20°~30°。于庆磊等利用RFPA程序,观察并分析了平台中心角对试样破坏模式和应力分布的影响。已有的研究大都利用有限元方法和某种强度理论进行分析,然后确定平台中心角及建立抗拉强度计算公式。然而岩石材料自身的非均质、不确定性往往被忽略而且岩石内部破坏过程也没有被重视。本文作者基于颗粒流数值模拟程序建立了平台巴西劈裂试验的细观离散元模型,研究了劈裂过程中裂纹发展规律。通过对不同平台中心角下的巴西平台试验的分析,确定了合理的平台中心角的值和该平台角下岩石抗拉强度的计算式。 1 岩石模型和测试过程的p 1.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.2 pbc试验过程 1.2.1 初始应力条件的设置 本次模型的建立可分为如下4步。 (2)通过半径扩大的方法得到颗粒集合体初始的各向同性应力条件。 (3)去除漂浮粒子。球的配位数小于4的颗粒定义为漂浮粒子,并予去除。 (4)设置平行黏结模型。这样平台巴西圆盘的岩石模型最终生成。 1.2.2 宏观参数标定 PFC的输入变量并不是人们所熟悉的宏观参数,而必须经过多次标定后才能找到的细观参数。所以使用该数值方法时,首先要假定其输入参数,然后对模拟结果和真实材料的力学性质进行比较;不断调整输入参数,比较其结果,最终两者结果相一致的一组参数就是所需的模型参数。 本次标定过程使用了某矿山的大理岩试样,所取的宏观参数具有一定的代表性。标定过程选择了表1所示的宏观参数,进行了巴西劈裂和单轴压缩虚拟试验。标定过程后,颗粒半径在0.16~0.26 mm范围之间,模型主要微观参数如表2所示。 1.2.3 墙与颗粒间的运动 在圆盘的2个平