0引言花岗岩是钨矿围岩中常见的一种岩石。
目前许多钨矿矿山为了延长矿山服务年限,急需加强矿山外围和深部勘探,这就需要进行大量的矿岩破碎工程。
冲击载荷作用下破碎岩石是目前运用较为广泛的一种破岩方法,但其破岩机理还有待进一步研究。
近年来,非线性科学中的分形理论、系统科学中的突变理论和数值模拟方法等现代科学理论开始渗入岩石破碎研究领域[1-5],取得了一些成果。
但是,分形理论在岩石破碎中主要被用于研究破碎块度,无法建立破碎结果与破碎的物理机制之间的相互关系;突变理论主要研究系统参数发生微小变化后结构的稳定与否,它在岩石破碎领域的应用目前主要局限于地震与岩爆等由缓慢变形到突然破裂的现象,尚没有用于研究凿岩等由于冲击载荷产生的破坏现象;数值模拟方法是用来沟通理论模型和实验研究的桥梁,但岩石的本构模型准确性、客观性以及岩石性质参数和边界条件还需要加以研究。
因此冲击载荷破岩的数值模拟与试验研究相接合对丰富冲击破岩理论具有重要意义。
1冲击载荷破岩特点1.1岩石断裂形态如图1所示,压头侵入硬岩时一般产生径向、中间和侧向裂纹,同时在压头下方还会形成一个密实核,在钠玻璃上进行侵入试验表明,形成的密实核近似于半球形,其主要特征是发生了剪切变形[6]。
由此可看出侧向裂纹是从剪切变形区底部起裂的。
侧向裂纹一般在卸载过程产生并扩展,中间裂纹产生于加载过程,并在卸载过程有部分弹性恢复;径向裂纹既可产生于加载过程,又可出现在卸载期间,但不论何时产生都在卸载过程继续发展。
1.2冲击破岩特点岩石在冲击载荷作用下将引起应力波在岩石中的传播。
对应力波的描述包括频谱(即频率分布、能量分布、波速、波长等,同时岩石中裂纹或缺陷也有尺寸分布或谱图,这样便能根据载荷作用特点与具体岩石对象确定微观与细观以及宏观尺度的划分标准,并找出分析重点,从而采用相应的方法来对其动态破坏过程进行研究。
典型的冲击载荷有机械冲击载荷与爆炸冲击载荷。
机械冲击是动能传播到一个系统,其发生传递的时间比该系统的自由振动周期要短。
冲击载荷的这种短暂特征,意味着在材料中会产生很高的应变文章编号:1009-0622(200705-0005-05冲击载荷作用下岩石破碎数值模拟及试验研究谢世勇1,王艳霞2,赵伏军3(1.赣州有色冶金研究所,江西赣州341000;2.河北工业大学电气与自动化学院,天津300130;3.湖南科技大学能源与安全工程学院,湖南湘潭411201摘要:运用数值模拟方法研究了冲击载荷作用下的岩石微观结构特性变化,通过对岩石内部应力分布及裂纹分布图可以得出冲击载荷破岩时岩石内部在垂直方向产生一个拉应力区和一个压应力区,刀具下方的拉应力最大;中间裂纹扩展速度比径向和侧向裂纹都要快,裂纹面呈三角形状。
在多功能试验装置上对花岗岩进行冲击载荷破岩试验,试验结果表明:单纯提高冲击载荷对破岩效果不会有很大的提高,存在一个最佳破岩冲击能(63J。
关键词:冲击载荷;岩石破碎;数值模拟;试验研究中图分类号:TU457文献标识码:A收稿日期:2007-08-21作者简介:谢世勇(1982-,男,江西南康人,硕士,从事矿山地压控制、矿岩破碎等方向研究。
压头密实核径向裂纹侧向裂纹中间裂纹图1侵入断裂形态第22卷第5期2007年10月Vol.22,No.5Oct.2007China Tungsten Industry第22卷速率,对材料的变形破碎行为和机制产生重大影响。
爆炸冲击载荷比机械冲击载荷的强度要高得多。
爆炸发生以后,气体产生的压力瞬间可达到一万多个兆帕,温度达到几千度,在高温高压的作用下,临近岩石受到强烈压缩,颗粒被压碎,甚至发生液化和汽化。
岩石因受爆炸产物的挤压而向外运动,形成一个空腔,同时在岩石中产生一种强冲击波向外传播。
随着传播距离的增加,波阵面压力下降,不再引起岩石结构的破坏,而只产生质点的振动,这就是爆炸引起的地震效应[7]。
M.P.马弗柳托夫[8]采用直径1.95mm的压头做了自由落体冲击岩石的试验,他研究了动载荷冲击岩石时的破碎问题,其结果如图2所示。
曲线1是当冲击功不大时,在岩石表面只能见到压头冲击的痕迹, 围绕压头边缘出现裂纹带;曲线2是当冲击功增加到一定值时,载荷达到第一极大值,压头底下的岩石发生与静载荷相似的破碎,称为第一次大体积破碎(剪切体崩离,主压力体压碎;曲线3是当冲击功大于产生第一次大体积破碎的能量,使得压头侵入深度增加;曲线4是当冲击功相当大时,载荷达到第二极大值,岩石产生第二次大体积破碎。
依次循环。
由上述分析可知,动载荷破岩同样具有跃进式的特点。
冲击破碎主要适用于脆性和弹脆性岩石,岩石的脆性愈明显、愈坚硬,破碎时的载荷极大值就愈大。
脆性岩石在外载荷达到极大值后产生破裂或大体积破碎,而弹塑性岩石在破坏后仍有部分强度,在一定的外部作用下方可产生体积破碎。
2数值模拟2.1岩石及PDC刀具模型岩石及PDC刀具均采用三维实体建模,动静组合破岩试验的花岗岩试件尺寸为1000mm×500mm×500mm,PDC刀具的规格是"13.4mm×15mm。
刀具的模型采用实际尺寸建模,由于岩石受刀具影响范围较小,所以取刀具与岩面接触正下方的100mm×50mm×50mm范围建模。
花岗岩试件和PDC刀具的力学参数见表1和2,划分网格后的岩石及刀具模型见图3。
表1花岩石试件的力学参数表2刀具材料的物理力学性能(a模型整体三维视图网格(b刀具与岩面接触点放大网格(c模型整体二维视图网格(d刀具与岩面接触点放大网格图3岩石与PDC刀具模型网格2.2加载与求解在试验中,刀具刀刃的锐边圆弧(即锋利的边缘与岩面接触,在模拟中为了简化运算,现假设刀具与岩面接触为点面接触,且在加载过程中刀具只有垂直方向的位移,无水平移动,即不存在切削现象。
由于只考虑施加单轴冲击载荷(48、63、78J,所以分析类型选择瞬态求解(Transient,加载时间为0.0002s(中等应变速率。
只需将岩石模型的底面施加位移约束(即x、y、z3个方向都约束,模型前后左右4个面为自由面。
冲击载荷直接施加在刀具与岩面相接触的节点上,设定100个子步,最小子步为50,最大子步为100,为了便于收敛,打开自动调整时间积分步长。
2.3结果分析为了分析在单轴冲击载荷作用下试件的垂直应力和裂纹分布随子载荷步变化的情况,选取具有代表性的子载荷步计算结果来分析。
图4为冲击能48J作用下各子载荷步在垂直方向上的应力等值线分布云图,图5为各子载荷步产生的裂纹分布图,图6为不同冲击能(48、63、78J时裂纹分布图。
体积密度g/mm3抗压强度MPa抗拉强度MPa弹性模量GPa泊松比2.640167.518.6670.21 密度g/mm3努普硬度GPa抗弯强度GPa抗压强度GPa断裂韧性弹性模量GPa泊松比耐磨性3.5250￣800.6￣1.17￣87￣9560￣8000.08250 24600.51.01.52.02.51234P/kNh/min图2岩石冲击力与浸深关系6第5期(1在冲击载荷作用的初期,即在第1载荷步的时候,此时在垂直方向产生比较大范围的压应力,拉应力集中出现在刀具的下方,但还不足以使岩石产生裂纹,如图4所示。
(2在第2载荷步时,此时垂直方向在刀具与岩面接触点的下方出现一个较大的拉应力区,如图4所示,裂纹最先在刀具与岩面接触点下方开始产生,但裂纹分布面较小,如图4和图5所示。
(3从第8载荷步开始,刀具下方的压应力进一步增大,同时拉应力也相应增大,岩石裂纹面进一步扩大,裂纹扩展迅速,形成较大面积的破碎,并且裂纹往深部比往两侧扩展要快,如图4和图5所示。
(4到第27载荷步时,应力分布较为稳定,在垂直方向产生一个拉应力区和一个压应力区,刀具下方的拉应力最大,裂纹扩展形成一个较稳定的“三角形”裂纹面,如图5所示。
第1步第2步第10步第12步第21步第27步图448J 冲击能垂直方向应力等值线云图图548J 冲击能作用下各子载荷步产生的裂纹分布图6不同冲击能时裂纹分布第2步第8步第10步第12步第21步第27步48J 63J 78J谢世勇,等:冲击载荷作用下岩石破碎数值模拟及试验研究7第22卷(5从图5和图6中可以看出,裂纹首先产生在刀具与岩面接触点上,然后中间裂纹扩展最快,径向裂纹比侧向裂纹扩展要快,岩石表面破碎较快;冲击载荷作用下岩石内部裂纹形成的破碎角较小;随着冲击载荷的增加,裂纹分布的深度增加较快,两侧的扩展则较缓慢,裂纹分布的面积也随即增加,但裂纹面积在63J时增加较快。
3试验研究3.1试验装置该装置主要由轴向静压加载装置、落锤冲击加载装置、刀具及夹具装置和测试系统等部分组成[9-10], 如图7所示。
图7多功能破岩试验台结构图1-机架;2-轴向加压油缸;3-冲击杆传动装置;4-冲击杆;5-升降横梁;6-刀具夹;7-冲击调速电机;8-皮带传动装置; 9-轨轮车;10-岩样;11-导轨冲击加载机构或冲击杆传动装置主要组成部件为机架、调速电机、皮带传动机构、齿轮传动机构、链轮传动机构、托销、冲击杆(锤、配重等。
刀具安装机构通过法兰盘及螺栓固定在升降横梁上,主要由刀头、碟形弹簧、刀具杆、传感器(分别测量静载荷和冲击应力波等组成。
在升降横梁板上可安装多把刀具同时破碎岩石,也可变换刀具夹安装滚刀或其他刀具破岩。
选用湖南飞碟新材料有限公司生产的圆柱形聚晶金刚石复合片(PDC压头,规格见表2。
3.2试验方法和步骤本试验中冲击高度为300mm,冲击能分别为48、63、78J。
每次破碎岩石之后,清除岩粉,测量破碎坑体积,每一试验重复3次,取其平均值。
试验完成后,采用DH3840应变放大器(动态应变仪和DH5932数据采集仪记录和分析试验数据。
实验步骤如下:(1将试验前制作好的花岗岩试件吊放在轨轮车平台上,在台面与试样之间放置不同规格的钢架,以随时调节试样高度。
(2试验参数设置。
花岗岩的冲击高度为300mm,加载波形为似正弦波,幅值设置超过试件的估计强度值。
因为存在动态实验,所以选择动载频率和数据采样频率高的数据采集/处理系统。
(3启动各加载设备,以不同能级的冲击能来实现冲击破岩的试验,获得各个对应的侵入深度、破岩体积等,以侵入深度、破岩体积和比能耗来评价破岩效果。
(4设置极限保护以保障设备安全。
3.3试验结果及分析岩石的破碎效果以破碎体积和破碎比能(即破碎单位体积岩石所消耗的能量来度量,由表3的试验结果可计算出破碎比能。
破碎比能[11]由公式(1确定: E S=W/V(1式中:ES为破碎比能,J/mm3;V为破碎岩石体积,mm3;W为冲击能量,J。
而W=P u/2式中:P为作用在刀具上的力,用WOB表示,N;u为刀具切削深度,m。