空孔直径对有围压直眼掏槽爆破破碎效果的影响.pdf

空孔直径对有围压直眼掏槽爆破破碎效果的影响张宪堂1,2，马力1,2，余辉1,2，张经双3，朱宝合4，周红敏1,2，王向阳1,2（1.山东科技大学山东省土木工程防灾减灾重点实验室,山东青岛266590；2.山东科技大学土木工程与建筑学院,山东青岛266590；3.安徽理工大学矿山地下工程教育部工程研究中心，安徽淮南232001；4.中国建筑第五工程局有限公司,湖南长沙410004）摘要：掏槽爆破是井巷爆破掘进的关键技术，高围压对掏槽爆破效果有明显的抑制效果，而空孔对于提高掏槽爆破效率具有重要意义，但目前相关研究并不充分。为研究中空孔直径对围压作用下掏槽爆破破碎效果的影响，制作边长为 600mm、砂胶质量比为 21 的立方体水泥砂浆试块，并通过电液伺服加载系统对试块施加双向等围压，进行直眼掏槽爆破水泥砂浆相似模型试验。对爆后爆破块度进行筛分统计分析，分析研究了不同中空孔直径下模型碎块块度分布、平均粒径，通过分形维数表征碎块分形特征，量化分析了碎块破碎程度。采用 ANSYS/LSDYNA 软件基于试验模型参数建立了数值模型，分析了不同空孔直径时不同时刻的应力云图。研究结果表明：炮孔装药量分别为2g 和3g时，中空孔直径为 24mm 和 16mm 的模型碎块分形维数最大，此时碎块粒径主要集中 1975mm的中部区间，且此时碎块总质量最大，分别为 3.128、3.615kg，空孔直径不同时，空孔附近产生的应力集中效应与反射拉伸波效果不同，使得掏槽爆破破碎效果不同，适宜的空孔直径能够有效利用空孔效应，空孔附近应力波峰值和作用范围较大，应力波作用时间较长，因此破碎效果较好；空孔直径过小，空孔附近应力较小，空孔直径过大时，应力波由空孔溢出，应力集中叠加后的压应力较小，且拉伸应力作用区域较小，空孔直径过大或过小都会使大块度碎块占比较高；装药量为 3g 时，空孔直径 32mm 与 16mm 模型碎块的分形维数差值大于 16mm 与无空孔模型碎块分形维数的差值，说明当空孔直径过大时对碎岩块度的影响更大；空孔直径同为 16mm 时，3g 装药量模型的碎块平均粒径小于 2g 装药量模型，分形维数及碎块总质量大于 2g 装药量模型试块，说明当空孔直径相同时增加装药量可以提高岩石的破碎度及槽腔体积。关键词：爆破掘进；掏槽爆破；空孔直径；分形维数；破碎效果；数值模拟中图分类号：TD235文献标志码：A文章编号：02532336（2023）10005510Influence of diameter of empty hole on the fragmentation effect of parallel cutblasting under confining pressureZHANGXiantang1,2,MALi1,2,YUHui1,2,ZHANGJingshuang3,ZHUBaohe4,ZHOUHongmin1,2,WANGXiangyang1,2（1.Shandong Key Laboratory of Disaster Prevention and Mitigation in Civil Engineering,Shandong University of Science and Technology,Qingdao 266590,China;2.College of Civil Engineering and Architecture,Shandong University of Science and Technology,Qingdao 266590,China;3.Engineering ResearchCenter of Underground Mine Construction,Ministry of Education,Anhui University of Science and Technology,Huainan 232001,China;4.China Construc-tion Fifth Engineering Division Corp.,Ltd.,Changsha 410004,China）Abstract:Cutblastingisthekeytechnologyforroadwayblastingexcavation.Highconfiningpressurehasobviousinhibitioneffectoncutblastingeffect,whileemptyholeisofgreatsignificancetoimprovecutblastingefficiency,buttherelevantresearchisnotsufficientsofar.收稿日期：20220830责任编辑：常琛DOI：10.13199/ki.cst.2022-1357基金项目：国家自然科学基金资助项目（51874189）；山东省自然科学基金资助项目（ZR2023ME106）；矿山地下工程教育部工程研究中心开放基金资助项目(JYBGCZX2021102)作者简介：张宪堂（1973）,男，河北井陉人，教授，博士生导师。E-mail：第51卷第10期煤炭科学技术Vol.51No.102023年10月CoalScienceandTechnologyOct.2023张宪堂，马力，余辉，等.空孔直径对有围压直眼掏槽爆破破碎效果的影响J.煤炭科学技术，2023，51（10）：5564.ZHANGXiantang，MALi，YUHui，et al.Influenceofdiameterofemptyholeonthefragmentationeffectofparal-lelcutblastingunderconfiningpressureJ.CoalScienceandTechnology，2023，51（10）：5564.55Tostudythecrushingeffectofthecutblastingwithdifferentemptyholediametersundertheconfiningpressure,cubiccementmortartestblockswhosesidelengthwas600mmweremadewithsandcementratioof2:1.Thetestblocksweresubjectedtoequalbi-axialconfiningpressurethroughanelectro-hydraulicservoloadingsystem,andtheparallelcutblastingcementmortarphysicalmodeltestwascarriedout.Aftertheblastingfragmentationwasscreenedandstatisticallyanalyzed,thefragmentationdistributionandaverageparticlesizeofmodelfragmentsunderdifferentemptyholediameterswereanalyzedandinvestigated.Withfractalfeaturesoffragmentscharacterizedbyfractaldimension,thefragmentationdegreewasquantitativelyanalyzed.ANSYS/LSDYNAsoftwarewasusedtoestablishnumericalmodelsbasedontheparametersofthetestmodel,andthestresscontouratdifferenttimeswithdifferentholediameterswasanalyzed.There-searchresultshowsthat,whentwoemptyholesarechargedwithdifferentchargeof2gand3g,respectively,themaximumfractaldimen-sionofthemodelfragmentcanbeattainedwithemptyholesindiameterof24mmand16mm,respectively,withtheparticlesizeofthefragmentmainlyconcentratedinthemiddlesegmentbetween19to75mm,andthelargesttotalmassofthefragmentis3.128kgand3.615kg,respectively;Whenthesizesofemptyholesaredifferent,thestressconcentrationeffectandthereflectedtensilewaveproducedneartheemptyholearealsodifferent,whichbringsthedifferenceofcutblastingfragmentationeffectandcaneffectivelybenefitfromtheemptyholeeffect.Thepeakvalueandactionrangeofstresswaveneartheemptyholearebigger,withalongactiontime,suggestingabet-tercrushingeffect;thediameteroftheemptyholeistoosmall,thestressaroundtheemptyholeissmall.Anoversizedemptyholewillleadtotheoverflowofthestresswavefromthehole.Thecompressivestressissmallafterthesuperpositionofstressconcentration,andtheactionareaoftensilestressisalsosmall.Theproportionoflargerockpieceswillbehighwithextremelylargeorsmallemptyhole.Whenchargedwith3g,thefractaldimensiondifferencebetweenemptyholemodelfragmentsof32mmand16mmisgreaterthanthedif-ferenceofmodelfragmentsbetween16mmandwithoutemptyhole,indicatingthattheinfluenceonthefragmentationofrockisgreaterwhentheemptyholediameteristoolarge.Whentheemptyholediameteristhesameas16mm,theaverageparticlesizeofthefragmentsinthemodelchargedwith3gislessthanthatinthemodelchargedwith2g,withthegreaterfractaldimensionandtotalmassofthefrag-mentsthanthatinthetestmodelchargedwith2g.Thisshowsthatanincreasedchargecancontributetorockfragmentationandthevolumeofthecutcavitywiththesamediameterofemptyholes.Key words:blastingexcavation；cutblasting；diameterofemptyhole；fractaldimension；fragmentationeffect；numericalsimulation0引言我国地下资源开采深度逐渐增加，钻孔爆破依然是岩巷掘进的主要手段1，围压条件下的爆破开采技术有待跟进，掏槽爆破可为后续巷道爆破提供新的自由面，有利于为后续爆破克服围压条件下岩石的夹制作用，是深部巷道爆破掘进的关键技术2。许多学者对于掏槽爆破进行了大量研究，也取得了诸多进展。王远来等3通过理论计算和LSDYNA软件数值模拟方法，进行了直眼掏槽爆破参数的优化研究，结果表明空孔具有导向和应力集中作用，孔径为 0.051m，装药孔与空孔间距为 0.30m 时，爆破效果较好。李兰等4通过理论分析和现场试验研究，改善装药系数和炮孔堵塞，增设大直径空孔，提高了掏槽爆破循环进尺。单仁亮等5分析了掏槽爆破作用机理，发现岩石爆破掏槽是由冲击波、应力波和爆生气体 3 者共同作用的结果。黄丹等6提出了中心双空孔超深滞后抛掷直眼掏槽，该爆破方法在降低炸药单耗的同时提高了炮孔利用率和有效进尺，爆破块度均匀。GAO 等7利用 ANSYS/LSDYNA 软件模拟了不同侧压力系数下复杂五孔掏槽爆破，分析了爆破裂缝侧压力系数的演化规律。左进京等8研究了掏槽爆破不同中空孔直径对槽腔体积、爆破块度的影响，试验结果表明空孔直径越大，槽腔体体积越小，岩石块度越小。ZHANG 等9通过数值模拟分析了不同空孔个数对掏槽爆破效果的影响，发现空孔具有导向和应力集中作用，多空孔布置时，相邻空孔导向方向不一致。郝英剑10在巷道爆破掘进中增设了大直径空孔，提高了掏槽爆破效率，缩短了工期。皇新宇等11基于 RHT 岩石本构模型，进行了不同地应力下四孔掏槽爆破数值模拟，分析了岩石损伤裂纹的演化规律和炮孔周围应力场变化过程。笔者等12给出了考虑碎块间挤压碰撞作用的节理裂隙岩体爆破块度分布的理论模型，并进行了数值模拟和实验验证。杨洋等13开展霍普金森杆冲击试验，发现随着冲击荷载的增大，爆破块度的平均粒径越来越小。于冰冰等14分析了顾北煤矿岩石巷道断面掘进原始方案，提出炮孔超深优化方案，并通过了现场试验验证。上述研究大多是针对无围压作用下的岩体爆破效果分析或围压作用下的数值模拟分析，关于围压作用下空孔直径对岩石破碎效果的试验分析还缺乏深入研究。采用自主设计的大型双榀组合式加载框架及电液伺服加载系统，以水泥砂浆为模型材料建立了双轴围压条件下掏槽爆破模型，分析了初始地应力作用下不同空孔直径时，岩石破碎块度分布特2023年第10期煤炭科学技术第51卷56征和分形维数。1掏槽爆破模型试验1.1试验方案试验掏槽眼布置形式为菱形布孔，模型中心设空孔。水泥砂浆爆破模型试块的尺寸为 600mm600mm600mm，根据几何相似理论确定模型几何参数，炮孔间距为 50mm，炮孔直径为 16mm，炮孔到空孔中心距离为 35mm，炮孔与中空孔深均为 160mm，装药深度为 80mm，炮孔堵塞长度为 80mm，实际最小抵抗线为 160mm。以含有 DDNP 起爆药的数码电子雷管代替炸药，每发电子雷管含有 1g 炸药，炮孔底部近似耦合装药，四炮孔同时起爆。炮孔填塞选用细砂和土作为填塞物，填塞时将填塞物捣实，避免在填塞时造成炮孔卡堵、填塞不实的现象。为能近似满足爆炸动力学相似，遵循“炸药爆炸能量相似”理论15设计装药量，确定炸药单耗为0.588kg/m3，单炮孔炸药用量为 2.3g。选择每炮孔装填 2 发、3 发电子雷管，即单炮孔装药量为 2、3g。通过 3g 装药量模型试验确定适宜的空孔直径大致范围，后续减小装药量为 2g，并在前期试验基础上细化空孔直径范围继续试验研究。炮孔布置示意如图 1 所示，具体试验方案见表 1。爆破模型炮孔空孔1111600600Hv图1炮孔布置示意Fig.1Layoutdiagramofblastholes采用自主设计的大型双榀组合式加载框架及电液伺服加载系统，对试验模型施加双向围压。试验加载区域面积为 400mm340mm，考虑到水泥砂浆为低强度脆性材料，应力取模型抗压强度的 10%，v=H=1.0MPa。表 1 爆破模型试验方案Table 1 Blasting model test scheme模型编号炮孔直径/mm炮孔间距/mm炮孔深度/mm空孔直径/mm单炮孔药量/gH01650160无空孔3H161650160163H321650160323E161650160162E241650160242E281650160282E3216501603221.2模型制作选用 M10 水泥砂浆作为模型试验的相似材料，选用 P.O325 普通硅酸盐水泥，砂子选用中砂，按照砂胶比为 21，浇筑水泥砂浆作为模型，浇筑过程中边浇筑边振捣，以减少试块中的气泡，常温下养护 28d（表 2）。表 2 水泥砂浆力学参数Table 2 Mechanical parameters of cement mortar密度/(kgm3)泊松比抗压强度/MPa抗拉强度/MPa弹性模量/GPa21750.2310.01.3224.172中空孔直眼掏槽破岩机理1)耦合应力场炮孔周边应力演化规律。围压作用下的爆破，岩石受到初始地应力场和爆炸动荷载的耦合作用。地应力能够抑制掏槽爆破损伤裂纹的扩展，随着地应力的增大，掏槽损伤区域减小，切向拉应力峰值降低11。耦合应力场下，炮孔周围某一点径向应力 r和切向应力 见式(1)：r=12(H+v)(1r21r2)+12(Hv)(1r21r2)(13r21r2)cos2+pc(r1r)=12(H+v)(1+r21r2)12(Hv)(1+3r21r2)cos2dpc(r1r)（1）式中，H、v分别为该点的水平应力与垂直应力，MPa；为计算点与炮孔中心的连线和水平方向的夹角；r1为炮孔半径，m；r 为该点至炮孔中心距离，m；为应力波传播过程的衰减系数；d为侧向应力系数。2)空孔应力集中效应。掏槽孔中炸药爆炸后会产生冲击波，随着在岩石中传播距离的增加，冲击波张宪堂等：空孔直径对有围压直眼掏槽爆破破碎效果的影响2023年第10期57逐渐衰减为应力波在岩石中继续传播，当到达中空孔孔壁时，部分应力波会发生反射，空孔附近发生应力重分布，此时空孔附近岩石中的应力值大于无空孔时的应力值，从而更有利于岩石的破碎，此为空孔的应力集中效应16。空孔附近的叠加应力场17表示为rr=12(1k2)(r)+(14k2+3k4)(+r)cos2=12(1+k2)(r)+(1+3k2)(+r)cos2r=12(1+2k23k4)(+r)cos2k=r2rb（2）式中，rr为应力重分布后计算点的径向应力，MPa；为应力重分布后计算点的切向应力，MPa；为应力重分布后计算点的剪应力，MPa；r为计算点与炮孔和空孔连心线夹角，()；r2为空孔半径，m；rb为空孔附近某点距离空孔中心距离，m（图 2）。3)自由面效应。空孔可为掏槽爆破提供一个新的自由面，当爆炸应力波传播到自由面时，应力波部分反射形成拉伸波，拉伸波作用于自由面处的岩石，由于岩石的抗拉强度仅为抗压强度的 1/81/15，如果反射的拉伸应力波的强度大于岩石抗拉强度，则自由面处岩石被拉伸破坏，此即为空孔的自由面效应18。3爆破块度筛分统计3.1碎块块度筛分析在掏槽爆破中，岩石破碎效果是反应爆破效果的一个重要指标，爆破块度过大会对碎岩的清理、搬运过程产生阻碍，爆破块度过小则说明炸药有剩余能量未被完全利用，使得炸药利用效率低，造成资源的浪费（图 3）。HHvvrbr1r2Lr炮孔空孔图2空孔应力集中效应示意Fig.2Schematicofstressconcentrationeffectofemptyhole表 3 不同空孔直径模型试验碎块块度筛分统计Table 3 Lumpiness screening statistics of model tests with different diameters of empty holes碎岩块度分布/mm90H0质量/kg0.0660.0690.3210.3580.5940.4770.579累计质量/kg0.0660.1350.4560.8141.4081.8852.464质量占比/%2.72.813.014.524.119.423.5累计质量占比/%2.75.518.533.057.176.5100.0H16质量/kg0.2410.5450.7660.6380.5540.5120.359累计质量/kg0.2410.7861.5522.192.7443.2563.615质量占比/%6.715.121.217.615.314.29.9累计质量占比/%6.721.843.060.675.990.1100.0H32质量/kg0.0280.0720.1360.3210.540.4760.866累计质量/kg0.0280.1000.2360.5571.0971.5732.439质量占比/%1.13.05.613.222.119.535.5累计质量占比/%1.14.19.722.945.064.5100.0E16质量/kg0.1110.1640.4020.4170.2000.3160.339累计质量/kg0.1110.2750.6771.0941.2941.6101.949质量占比/%5.78.420.621.410.316.217.4累计质量占比/%5.714.134.756.166.482.6100.0E24质量/kg0.2760.3820.5570.5780.5380.3650.432累计质量/kg0.2760.6581.2151.7932.3312.6963.128质量占比/%8.812.217.818.517.211.713.8累计质量占比/%8.821.038.857.374.586.2100.0E28质量/kg0.0390.0740.1080.1680.4460.5470.736累计质量/kg0.0390.1130.2210.3890.8351.3822.118质量占比/%1.83.55.17.921.125.834.8累计质量占比/%1.85.310.418.339.465.2100.02023年第10期煤炭科学技术第51卷58爆破试验完成后，根据文献 1921 选取破碎区半径为 0.08m 范围内的水泥砂浆碎块进行筛分。采用国家标准的等级石子筛将碎块块度分为 7 个级配，对每个级配下的碎块进行称重，对不同空孔直径模型的碎块质量、累计质量及质量占比等进行统计，得到碎块块度筛分见表 3。为更直观表征各试验方案碎块块度分布情况，根据块度筛分析统计结果绘制碎岩块度分布折线图如图 4 所示。图3部分爆破碎块Fig.3PartialblastingfragmentsE16E24E289005101520253035质量百分比/%碎块块度分布/mmH0H16H32(a)3 g 装药量模型试验9005101520253035质量百分比/%碎块块度分布/mm(b)2 g 装药量模型试验图4碎岩块度分布Fig.4Distributionoffragmentationdegree由图 4 中可以看出，6 种试验方案在各个筛分等级内均有碎块分布，且大块度碎块质量占比差异较大。由于围压的存在，掏槽爆破需要更多的能量克服岩石的夹制作用形成槽腔，6 种试验方案的小块度碎块均较少。炮孔装药量为 3g 时，H0 和 H32 模型试块碎块块度分布趋势大致相同，小块度碎块较少而大块度碎块较多，碎块块度分布主要集中在粒径大于 53mm 区间，此区间碎块质量占比分别为 67.0%、77.1%。这是由于 H0 试块未布设空孔，相比其他试块需要更多的能量克服岩石的夹制作用，导致没有足够的能量使岩石进一步破碎。空孔直径为 16mm时，碎块块度分布较为均匀。炮孔装药量为 2g 时，空孔直径为 28mm 的模型试块小块度碎块最少，E16 和 E24 碎块块度整体分布趋势相近，相比于 E16 试验方案，E24 小块度碎块较多而大块度碎块较少。但碎块粒径小于 19mm时，E24 方案碎块各筛分等级质量占比均大于 E16方案，碎块粒径大于 75mm 时，E24 方案碎块各筛分等级质量占比均小于 E16 方案。E24 试验方案碎块块度分布主要集中在 1975mm，质量占比为 53.5%，超过半数，其中 3 种筛分等级质量占比接近。E28 试验方案粒径小于 53mm 的碎块所占质量百分比较少，大部分为粒径大于 53mm 的大块度碎块。3.2碎岩平均粒径分析为进一步分析碎块粒径分布情况，利用平均破碎粒径特征值 ds来表示掏槽爆破后的岩石破碎程度，计算式如式(3)：ds=ridiri（3）式中，ds为围岩岩石破碎平均粒径特征值；di为不同孔径下岩石粒径尺寸；ri为粒径为 di时，对应粒径质量占比。由式(3)计算得各试验方案模型碎块平均粒径见表 4。表 4 碎块平均粒径Table 4 Average particle size of fragments模型编号H0H16H32E16E24E28平均粒径/mm63.746.271.452.248.273.0由表 4 中可以看出，E16 和 E24 试验方案碎块平均粒径较小。H16 平均粒径小于 E16，当空孔直径相同时增加装药量可以提高岩石的破碎度。E28 试验方案碎块平均粒径大于 E16，H32 试验方案碎块平均粒径大于 H0，说明当空孔直径过大时对碎块块度的影响大于空孔直径较小或没有空孔时。张宪堂等：空孔直径对有围压直眼掏槽爆破破碎效果的影响2023年第10期594分形维数分析已有研究表明2224，岩石破碎后的块度具有分形特征，分形维数是衡量岩石破碎程度的标尺，可以清晰地反映破碎过程与破碎产物的联系，是一种很好的研究岩石破碎的方法。GGS 分布函数是碎岩块度分形分析中常用的计算方法25。GGS 分布函数的计算公式为Y=(rrm)（4）式中，Y 为筛分等级 r 以下的碎岩累计质量占比；rm为碎岩最大尺寸；为碎岩分布指数。Y=M(r)MT=(rrm)（5）式中，M(r)为尺寸小于 r 的碎岩累积质量；MT为碎岩总质量；对上式求导得：dM(r)r1dr（6）设 M 为粒径大于等于 r 的碎岩质量，则 M=MTM(r)：dM=dM(r)r1dr（7）代入碎岩块度个数与质量关系：dM r3dN（8）碎岩分形维数 D 与特征尺寸 r 及等效粒径大于等于 r 的碎块个数 N 存在以下关系：N rD（9）对式(9)求导得：dN rD1dr（10）将式(10)代入式(8)得：dM(r)r3D1dr（11）结合式(7)，式(11)得：D=3（12）对式(5)两边求取对数，可得近似一元一次方程表达形式，y=ax+b。因此可将各筛分等级的碎岩质量累积百分数M(r)/MT与碎岩特征尺寸 r 在双对数坐标系下进行线性拟合，所得直线方程的斜率 a 即为碎岩分布指数。分维数D值越大，表明爆后生成的小块度岩石占比越大26。当 D=2时，各尺度区间的碎块质量比例相等；当 0D2时，大尺度区间碎块所占质量比例较大，破碎程度低；当 2D3时，小尺度区间碎块所占质量比例较大，破碎程度高。将各试验方案分形参数拟合得到图 5，为便于观察分形维数与空孔直径变化关系，拟合结果的相关系数 R2均大于 0.96，拟合效果良好。由拟合参数表 5 和分形维数与空孔直径关系图 6 可以看出，两组试验方案分形维数随着空孔直径的增大，均呈现先增大后减小的趋势，分形维数大于 2 的试验方案为 E16、E24、H16，H32 试验方案分形维数最小，大块度碎块比例最高。0.500.751.001.251.501.752.0000.51.01.52.0H0H16H32lg(M(r)/MT)lg(M(r)/MT)lg r0.51.01.52.0E16E24E28(a)3 g 装药量模型试验0.500.751.001.251.501.752.00lg r(b)2 g 装药量模型试验图5不同试验方案分形维数拟合曲线Fig.5Fractaldimensionfittingcurvesofdifferenttestschemes两组试验中 E24、H16 分形维数最大，碎块破碎程度最好，且产生的碎块总质量最大，分别为 3.128、3.615kg，即形成的槽腔体积最大。炮孔装药量分别为 2g 和 3g 时，掏槽爆破破碎效果较好的试块中空孔直径为 24mm 和 16mm。E32 未能产生槽腔及 H32 分形维数最小的原因，分析是由于空孔直径过大导致过多的爆炸应力波和爆生气体从空孔溢出，剩余能量不能使岩石进一步破碎，使得大块度碎块占比过大。E16 分形维数小于 E24 的原因，分析是由于空孔直径较小时，爆炸应力波到达空孔壁反射后形成的拉伸应力波较少，爆炸产生的能量不足以使岩石更好破碎。因此，适宜的空孔直径能够更好地发挥空孔效应，有效利用爆破能量，促进岩石的破碎，提高碎岩小块率，降低大块率，使得碎岩块度更加均匀。而空孔直径过大或2023年第10期煤炭科学技术第51卷60过小均不能有效利用爆炸能量，导致碎岩块度较大，不利于实际生产中的排矸。40481216202428321.51.61.71.81.92.02.12.2分形维数空孔直径/mmEH图6空孔直径与分形维数关系Fig.6Relationshipbetweenemptyholediameterandfractaldimension5掏槽爆破数值模拟分析各试验模型破碎程度的不同是由于中空孔直径不同导致空孔附近应力不同，从而影响了爆破破碎效果，而模型试验无法直观分析空孔附近应力分布情况，因此有必要通过数值模拟手段分析应力分布，并与试验结果相互验证。5.1掏槽爆破数值模型V=H=11)模型建立。为验证试验结果及分析的正确性，根据试验所用的模型尺寸，利用 ANSYS/LSDYNA有限元软件建立了装药量 2g，4 种不同空孔直径的掏槽爆破数值模型，分别记为 N16、N24、N28、N32。模型尺寸为 600mm600mm10mm，取双轴围压大小MPa（图 7）。同时建立空孔直径 24mm的无围压模型，分析围压对爆炸应力的影响。2)模型材料与状态方程。岩石单元类型设置为SOLID164号单元，岩石材料使用 RHT 材料模型。模型关键字为:*MAT_RHT。数值模拟岩石力学参数与相似模型试验材料保持一致，并通过参考相关文献确定27-28。炸药模型关键字为*MAT_HIGH_EXPLOS-IVE_BURN，采用 JWL 状态方程，JWL 状态方程用于反映爆炸后爆轰结果的体积和压力变化，JWL 状态方程如式(13)所示。=A(1R1V)eR1V+B(1R2V)eR2V+E0V（13）式中，为爆轰产物的压力；V 为爆轰产物的相对体积；E0为初始内能密度；A、B、R1、R2、为状态方程参数（表 6）。图7数值模型Fig.7Numericalmodel表 6 炸药及其状态方程参数Table 6 Explosive and its state equation parameter/(kgm3)D/(ms1)A/GPaB/GPaR1R2E0/GPa12003202290.1824.20.90.154.1925.2模拟结果分析由图 8 分析可知，炸药起爆后，各炮孔爆炸应力波以圆形向外传播，空孔附近出现应力叠加，应力波在 40s 时到达空孔孔壁，并在空孔壁处发生反射，形成反射拉伸波，之后爆炸应力波在空孔附近形成多次反射应力波，应力波作用范围主要集中在炮孔及空孔附近的掏槽区域，220s 后应力波逐渐消散，表 5 碎块 lg rlg(M(r)/MT)拟合参数Table 5 Fitting parameters of fragment lg rlg(M(r)/M)参数试验编号拟合曲线a(斜率)b(截距)分形维数D调整后R2lgrlg(M(r)/MT)拟合参数H0y=1.254x0.521.2540.5201.7460.973H16y=0.888x+0.2910.8880.2912.1120.977H32y=1.497x1.0481.4971.0481.5030.991E16y=1.016x+0.0650.9630.1102.0370.981E24y=0.811x+0.4080.8110.4082.1890.994E28y=1.302x0.7241.3090.7361.6910.968张宪堂等：空孔直径对有围压直眼掏槽爆破破碎效果的影响2023年第10期61试块破碎完成形成槽腔。N24 模型应力峰值最大，且应力峰值作用时间最长，作用面积最大，因此破碎效果较好，N16 次之，N32 模型由于空孔较大，应力波由空孔溢出，应力集中叠加后的压应力较小，且拉伸应力作用区域较小，因此未能破碎试块。220s 后应力波逐渐消散，试块破碎完成形成槽腔。模型应力状态与试验破碎效果相吻合，较好地验证了试验结果。由有无围压时的应力云图 9 对比分析可知，无围压时应力传播速度较快，相同时刻下无围压模型的应力云图面积较大，且空孔和炮孔附近的压应力及拉应力相较于有围压时均更大，围压对爆炸应力波的传播有抑制作用。6结论1)空孔直径不同时，空孔附近产生的应力集中效应与反射拉伸波效果不同，使得掏槽爆破破碎效果不同，适宜的空孔直径能够有效利用空孔效应，空孔附近应力波峰值和作用范围较大，应力波作用时间较长，因此破碎效果较好；空孔直径过小，空孔附近应力较小，空孔直径过大时，应力波由空孔溢出，t=40 st=100 st=160 st=220 s(a)模型 N16 不同时刻应力云图t=40 st=100 st=160 st=220 s(b)模型 N24 不同时刻应力云图t=40 st=100 st=160 st=220 s(c)模型 N28 不同时刻应力云图t=40 st=100 st=160 st=220 s(d)模型 N32不同时刻应力云图应力/MPa321.0288.7256.4224.1191.9159.6127.394.9962.7130.421.867应力/MPa348.7213.7278.6243.5208.5173.4138.4103.368.2333.171.886应力/MPa278.0250.0222.0194.1166.1138.1110.182.1254.1426.151.829应力/MPa240.1215.9191.7167.5143.2119.094.8170.5946.3722.152.065图8各模型不同时刻应力云图Fig.8Stressclouddiagramatdifferenttimes2023年第10期煤炭科学技术第51卷62应力集中叠加后的压应力较小，且拉伸应力作用区域较小，空孔直径过大或过小都会使大块度碎块占比较高。2)空孔直径为 32mm 与 16mm 模型碎块的分形维数差值大于空孔直径 16mm 与无空孔模型碎块分形维数的差值，说明当空孔直径过大时对碎岩块度的影响大于空孔直径较小或没有空孔时。3)空孔直径同为 16mm 时，3g 装药量模型试块的碎块平均粒径小于 2g 装药量模型试块，分形维数及碎块总质量大于 2g 装药量模型试块，说明当空孔直径相同时增加装药量可以提高岩石的破碎度及槽腔体积。4)相同时刻下无围压模型的应力云图面积较大，且空孔和炮孔附近的压应力及拉应力相较于有围压时均更大，围压对爆炸应力波的传播和峰值有抑制作用。参考文献(References)：杨仁树，张召冉，安晨，等.煤矿岩巷掘进爆破掏槽孔超深问题探讨J.煤炭科学技术，2020，48（1）：1023.YANGRenshu，ZHANGZhaoran，ANChen，et al.Discussiononultra-deepdepthproblemofslotholeinblastingexcavationofrockroadway in coal mineJ.Coal Science and Technology，2020，48（1）：1023.1李刚，牛磊，李文龙.煤矿硬岩巷道快速掘进技术研究J.煤炭科学技术，2018，46（11）：1320.LIGang，NIULei，