在拱腰爆炸荷载作用下锚杆动态响应数值分析.pdf
第3 2 卷第2 期 爆破 V o l 3 2N 0 2 2 0 1 5 年6 月B L A S T I N G J u n .2 0 1 5 d o i 1 0 .3 9 6 3 /j .i s s n .1 0 0 1 4 8 7 X .2 0 1 5 .0 2 .0 0 6 在拱腰爆炸荷载作用下锚杆动态响应数值分析 王光勇1 , 2 ,郭晓燕3 ,王超1 ,余永强1 1 .河南理工大学土木工程学院,焦作4 5 4 0 0 0 ;2 .中国科学技术大学近代力学系,合肥2 3 0 0 2 6 ; 3 .安阳师范学院审计处,安阳4 5 5 0 0 0 摘要为研究拱腰集中装药下锚杆的动载响应特性,利用有限元分析程序K S D Y N A3 D ,对锚杆的轴向 应力分布规律进行了显示动力分析。对比爆心径向同一距离测点的实测压应力时程曲线,模拟结果曲线峰 值、形态和趋势与之较为一致。迎爆侧受压强度明显大于背爆侧,最强的受压部位在迎爆侧拱腰;迎爆侧拱 顶和拱腰受拉强度比较大,最强受拉部位在拱顶。受压峰值和受拉峰值从锚头到锚端都是先增加后减小,但 最大峰值位置不同,并利用应力波原理和最大拉应力瞬时断裂准则分析其原因。 关键词地下工程;锚杆;拱腰;爆炸荷载;动态响应;数值分析 中图分类号T D 3 5 3文献标识码A文章编号1 0 0 l 一4 8 7 X 2 0 1 5 0 2 0 0 3 3 一0 6 N u m e r i c a lA n a l y s i so fD y n a m i cR e s p o n s eo fR o c k B o l t su n d e rE x p l o s i o nL o a df r o mH a n c e W A N GG u a n g .y o n 9 1 ”,G U OX i a o y a n 3 ,W A N GC h a 0 1 ,Y UY o n g q i a n 9 1 1 .S c h o o lo fC i v i lE n g i n e e r i n g ,H e n a nP o l y t e c h n i cU n i v e r s i t y ,J i a o z u o4 5 4 0 0 3 ,C h i n a ; 2 .D e p a r t m e n to fM o d e mM e c h a n i c s ,U n i v e r s i t yo fS c i e n c ea n dT e c h n o l o g yo fC h i n a ,H e f e i 2 3 0 0 2 6 ,C h i n a ;3 .A u d i t i n gD e p a r t m e n t ,A n y a n gN o r m a lU n i v e r s i t y ,A n y a n g4 5 5 0 0 0 ,C h i n a A b s t r a c t I no r d e rt os t u d yd y n a m i cr e s p o n s eo fr o c kb o l t su n d e re x p l o s i o nt o a df r o mh a n c e ,t h ed i s t r i b u t i o no f a x i a lp r e s s u r eo fr o c kb o l t sw a sa n a l y z e db yL s D Y N A3D .C o m p a r i n gp r e s s u r ec u r v e sf r o mt h es i m u l a t i o na n dt h e t e s t ,t h ep e a kv a l u e s ,s h a p ea n dt r e n d so ft h er e s u l t sw e r ei nu n i f o r m .T h ep e a kp r e s s u r eo nf a c eo fb l a s tw a v ew a s s t r o n g e rt h a nt h a to nb a c ko fb l a s tw a v e ,a n dt h el o c a t i o no ft h em a x i m u mp e a kp r e s s u r el i e di nt h eh a n c e .C o m p a r e d w i t ho t h e rp a r t so ft h et u n n e l ,t h el a r g e rt e n s i o nw a si nt h er a n g eo fv a u l ta n dh a n c e ,a n dt h em a x i m u mt e n s i o nl i e di n t h ev a u l t .F r o ma n c h o rh e a dt oa n c h o r ,t h ep e a kp r e s s u r ea n dp e a kt e n s i o ni n c r e a s e df i r s t l y ,a n dt h e nd e c r e a s e d ,b u t t h el o c a t i o n sw e r ed i f f e r e n t .T h ep r i n c i p l e so fs t r e s sw a v e sa n di n s t a n t a n e o u sf r a c t u r eo ft h em a x i m u mt e n s i o nw e r e a p p l i e dt od i s c l o s et h ei n n e rm e c h a n i s m . K e yw o r d s u n d e r g r o u n de n g l n e e r i n g ;r o c kb o l t s ;h a n c e ;e x p l o s i o nl o a d ;d y n a m i cr e s p o n s e ;n u m e r i c a la n a l y s i s 随着锚杆加固技术在地下工程中的广泛应 用⋯,关于锚杆加固机理及其自身力学性能方面的 收稿日期2 0 1 5 一O l 一2 4 作者简介王光勇 1 9 7 7 一 ,男,副教授、博士后,2 0 0 0 年毕业于河南 理工大学土木工程学院矿井建设专业,主要从事岩土工程 动载试验方面的教学与研究工作, E m a i l w g y 2 0 0 3 m a i l .u s t c .e d u .c n 。 基金项目国家自然科学基金 4 1 0 7 2 2 2 5 ;河南省科技攻关 1 4 2 1 0 2 2 1 0 43 3 ;河南省教育厅科学技术研究重点项目 科技攻关计划 1 4 A 4 4 0 0 0 3 研究一直是科研热点。对于静载下的情况,科研人 员通过大量研究已有了较成熟的认识。2 一。动载情 况下锚杆的加固机理及响应特性研究方面,国内外 学者也做了一些有益研究∞‘1 川,但多是针对吸能锚 杆设计和锚杆质量的无损检测,并没有考虑与围岩 相互作用和动载作用方向,即使考虑了与围岩相互 作用也多是在顶爆情况下进行,关于拱腰侧爆下锚 杆动载响应的研究资料却相对很少。拱腰侧爆是一 万方数据 3 4爆破 2 0 1 5 年6 月 种比顶爆更普遍的情况,爆炸源位置的变化必然导 致锚杆动载响应特性有所差异。由于爆炸动力学过 程非常复杂,很难进行精确解析分析,目前最常用的 方法有数值分析与模型试验。模型试验周期长、成 本高,而且对锚杆动载响应参数的实测难度较大。 相比之下,有限单元法和有限差分法等数值方法却 能有效处理结构与介质的相互作用问题,给出大量 的结构动载响应信息‘18 I 。因此,采用有限单元方法 模拟拱腰侧爆集中装药下,全长粘结式锚杆与应力 波的相互作用,并对该类型锚杆的动载响应特性进 行研究。 1 数值模型概况 1 .1 数值计算模型 计算模型选自文献[ 1 5 ] 中室内试验模型的一 部分,取宽高厚 2 4 0c m 2 3 0c m 4c m ,作近 似二维数值模拟,模型硐室跨度为6 0c m ,加固硐室 的锚杆长度为2 4c m ,间距4c m ,具体见图1 。为计 算简便,建立单排锚杆来模拟加固洞室。为了模拟 出半无限体中应力波的传播规律,模型左右边界、前 后边界、下边界设置为透射边界,左右边界和前后边 界水平位移设为零,底边界固定,下边界设置三个方 向都固定。上表面为自由边界,忽略岩体自重对应 力波传播造成的影响。炸药、围岩单元类型选用3 D S o l i d l 6 4 单元,加固锚杆单元类型为B e a m l 6 1 单元, 炸药、围岩和锚杆用节点法划分网格。 图1数值计算模型 单位e n l F i g .1 M o d e lo fn u m e r i c a la n a l y s i s u n i t C n l 1 .2 材料模型及参数 围岩材料模型采用压力增强的弹塑性本构模型 M A T S O I L C O N C R E T E ;锚杆材料选用弹性本构模 型M A T P L A S T I C K I N E M A T I C ;炸药模型使用L S D Y N A 3 0 程序提供的高性能炸药材料模型M A T H I G H E X P L O S I V E B U R N ,高性能炸药爆轰产物压 力一体积关系采用J w L 状态方程9 | ,其P y 关系 如式 1 P A 一南 e 幽r 日 一南 e 吨V 等 1 式中P 、E 。、V 分别表示爆轰产物单元压力、初始内 能密度和相对体积;相关模型材料计算参数见表1 , 其中A 、B 、尺1 、R 2 、∞为J W L 方程常数;p 1 、P 、P 3 分别 为炸药、围岩介质、锚杆的材料密度;P 。.为炸药爆轰 压力;D 。为炸药的爆速;E 。为爆轰初始内能;G 、K 为炸药材料的剪切模量、体积弹性模量;E 、肛为锚杆 材料的弹性模量、泊松比。 表1T N T 炸药、J w L 方程、模型及锚杆材料参数 T a b l e1P a r a m e t e r so fT N T 。J W L ,m o d e la n db o l tm a t e r i a l G /G P aK /G P aP 3 / g c m 。 E /G P a “ 0 .8 5 60 .9 5 82 .7 27 60 .3 4 2 数值模拟结果分析 2 .1 模拟、实测压应力时程曲线对比 爆压是造成洞室破坏的重要因素之一,根据爆 压时程曲线可看出应力波的作用规律。在爆源与迎 爆侧拱腰连心线上,提取距爆心4 0c m 点P 处 见 图1 围岩单元的压应力时程曲线,并进行炸药当量 折算,然后与实测爆心径向同一距离上测点的压应 力时程曲线进行对比,其结果如图2 所示。从峰值 上看,模拟与实测压应力峰值比较接近,模拟峰值稍 大些;从形态上看,模拟与实测较为相似,均由上升 与下降阶段组成,压应力变化趋势都是先增大后减 小,最后在零点附近达到平衡。其中,模拟时程曲线 上升、下降比实测的要快,应力波的作用时间比实测 的要短,残余压应力较实测的要小。这主要与围岩 密实程度有关,模拟中围岩为均匀整体,而试验中围 岩是用水泥砂浆制成,难免有孑L 隙存在,应力波在孔 隙中发生反射与衍射,既延长了其在模型中的作用 时间,也有一定的强度衰减。因而,实测时程曲线峰 值要小,作用时间要长。 万方数据 第3 2 卷第2 期王光勇,郭晓燕,王超,等在拱腰爆炸荷载作用下锚杆动态响应数值分析 3 5 r 间/m s 图2 模拟、实测压应力时程曲线对比 F i g .2C o m p a r i s o np r e s s u r eC u r V e So fa n a l y s i sa n dt e s t 2 .2 锚杆轴向应力分布规律 要想提高锚固围岩抗爆加固效果,就必须知道 在动载作用下锚杆的受力情况,数值分析提取每根 锚杆距锚头0m m 、6 0m m 、2 0m m 、1 8 0m m 、2 4 0m m 的数据,图3 是6 号锚杆5 个选取点轴应力时程曲 线,每个选取点的曲线都是先受压后受拉,最后趋于 稳定受拉。锚杆是先受压后受拉。图4 和图5 分别 是锚杆轴向受压和受拉峰值应力分布图,从图中分 析可知在迎爆侧拱腰锚杆受压作用最强,超过了 1 2 0M P a ,比其它锚杆受压峰值高6 倍以上;随着与 锚头距离增大,受压峰值先增加后减小,锚杆6 离锚 头1 8 0m m ,比其它位置的锚杆更偏向锚端;受拉峰 值也是先增加后减小,锚杆6 离锚头6 0m m ,比其它 位置的锚杆更偏向锚头。这主要是因为拱部集中装 药爆炸初始时刻,迎爆侧拱腰围岩最先受应力波作 用,围岩受压发生位移,近爆源围岩位移较大,而近 洞室锚杆加固区位移则相对较小,从而使得锚杆先 受压,并且锚杆受压越来越强。随着应力波向洞室 方向推进,距洞室表面越近的地方受到压应力作用 而向洞室内发生的位移越大,从而使前后位移差越 来越小,导致锚杆受压减弱。当应力波传至洞室自 由面边界时会发生反射,形成的反射拉伸波使近洞 室表层围岩受拉产生层裂现象,如图6 所示。从层 裂产生机理可以知道离洞室表面一定的厚度地方会 产生裂纹,此处的受拉强度最大,根据最大拉应力瞬 时断裂准则公式 2 ,层裂的厚度与应力波的波长 成正比旧⋯,与应力波幅值成反比,而锚杆6 应力波 作用最强,波长最短,故拱腰锚杆6 比其它位置的锚 杆更偏向锚头。 6 - 虿A 。瓦O “ c 2 式中6 ,是第一层裂的厚度;A 是波长;盯。为动态断 裂强度;O r 。应力波幅值。 背爆侧的锚杆轴向应力分布规律如图7 和图8 所示,从图轴力分布规律可以发现背爆侧几乎不受 压,只受拉;受压和受拉的最大峰值都比迎爆侧有大 幅度的减小,这是波与洞室作用和衰减的结果。从 锚头到锚端,受拉峰值先增加后减小,并且与迎爆侧 受拉峰值位置相比,更远离锚头。这是由于随着应 力波绕过洞室,波的绕射和反射会使应力波的波长 增加,幅值减弱,根据应力反射拉伸原理,背爆侧的 受拉峰值位置更远离锚头。 一0 ;趣~蜘一三磐 ∥爵龇l 4 ⋯ 图36 号锚杆轴应力时程曲线 F i g .3 A x i a lp r e s s u r ec u r v e so ft h eN o .6r o c kb o l t 皇一2 R 一4 望一6 擤一8 嘉芝 图4 迎爆侧锚杆轴向压应力分布 F i g .4 D i s t r i b u t i o no fa x i a lp r e s s u r eo fr o c kb o l t so nf a c e 图5 迎爆侧锚杆轴向拉应力分布 F i g .5 D i s t r i b u t i o no fa x i a lt e n s i l es t r e s s o fr o c kb o l t sO Nf a c e 2 .3 锚杆时程曲线分析 为了全面了解锚杆的动态响应情况,分别对迎、 背爆相应位置的锚杆进行分析,见图l 。从锚杆的 时程曲线可以看出整个形态与文献[ 1 4 ] 的试验所 得锚杆时程曲线相似,最后残余强度都是拉应力,说 明锚杆起到加固的效果。从前面的叙述可以知道受 拉和受压的最大值在锚杆中间位置附近,为了比较 0 O O 0 O 0 O O 0 O 0 0 o o 4 之4 “矗m 他H 孟窖、R倒厘霹士掌擐 万方数据 爆破2 0 1 5 年6 月 不同锚杆最强部位的受力特点,故取其中间位置的 锚杆作比较。本次模拟共计算了1 0m s ,结果显示 集中装药产生的球面应力波强度在4m s 后已衰减 趋于稳定,故只取计算中前4m s 的时程曲线进行分 析。图9 和图1 0 分别是迎、背爆侧不同位置锚杆中 间单元轴向应力时程曲线对比图,锚杆位置见图1 。 从图中可见迎爆侧除直墙脚锚杆 9 号 只受拉外, 从拱顶到直墙腰处的锚杆 5 、6 、7 、8 号 都是先受压 后受拉;背爆侧所有监测锚杆 1 、2 、3 、4 号 主要表 现为受拉应力作用。受压强度方面,距离爆源最近 的迎爆拱腰锚杆中间单元压应力最大,达到 1 0 0M P a 左右,其余锚杆的受压作用要小很多,最大 强度在1 3M P a 左右,相差约7 倍,这是由于球形爆 炸应力波在传播过程中,其强度会随爆心作用距离 的增大而迅速衰减,同时与应力波作用的角度发生 改变,在洞室表面发生波的反射和绕射,造成了锚杆 5 、7 、8 、9 相对于6 受压小得多。在受拉强度方面, 不仅迎爆侧拱部锚杆受拉作用明显,背爆侧拱部锚 杆受拉作用也很强。在迎爆侧,从拱顶到直墙脚 5 9 ,锚杆受拉强度整体上呈现减小趋势,其中锚 杆6 和7 以及8 和9 相差比较小,最大受拉峰值在 拱顶,为9 2M P a ;背爆侧受拉情况从拱腰到直墙脚 变化逐渐减小,但在侧墙的锚杆 1 ,2 变化不大,都 是8M P a 左右,这是因为侧墙波的传播距离长,同时 反射和绕射比较彻底。因此,在拱腰爆炸荷载作用 下,在洞室加固设计、施工时,应该对拱顶和拱腰围 岩重点加强。 图6 层裂图 F i g .6S p a l l i n gg r a p h 3 结论 利用数值分析软件,分析了锚固洞室在拱腰爆 炸荷载作用锚杆动态响应,得出以下几点结论 1 对比爆心径向同一距离上测点的压应力时 程曲线,模拟与实测时程曲线的峰值较为接近,形态 与趋势较为一致,表明该数值模拟结果可靠。 日 凸_ 塞 、 R 倒 匠 暴 妊 擐 锚杆长度/m m 图7 背爆侧锚杆轴向压应力分布 F i g .7 D i s t r i b u t i o no fa x i a lp r e s s u r eo fr o c kb o l t so nb a c k 日 山 茎 \ R 创 匠 暴 蜒 掘 锚杆长度/m m 图8 背爆侧锚杆轴向拉应力分布 F i g .8 D i s t r i b u t i o no fa x i a lt e n s i l e s t r e s so fr o c kb o l t so nb a c k 图9 迎爆侧锚杆中间单元轴向应力对比 F i g .9 A x i a lp r e s s u r ec o m p a r i s o n0 1 1 m i d d l ee l e m e n to ff a c er o c kb o l t s 日 山 岂 \ R 僖 厘 霹 套 据 图1 0 背爆侧锚杆中间单元轴向应力对比 F i g .1 0 A x i a lp r e s s u r ec o m p a r i s o no n m i d d l ee l e m e n to fb a c kr o c kb o l t s 万方数据 第3 2 卷第2 期王光勇,郭晓燕,王超,等在拱腰爆炸荷载作用下锚杆动态响应数值分析 3 7 2 在迎爆侧拱腰锚杆整体受压作用最强,随 着与锚头距离增大,受压峰值先增加后减小,锚杆6 离锚头1 8 0m m ,比其它位置的锚杆更偏向锚端;受 拉峰值也是先增加后减小,锚杆6 离锚头6 0m m ,比 其它位置的锚杆更偏向锚头。 3 背爆侧几乎不受压,只受拉;受压和受拉的 最大峰值都比迎爆侧有大幅度的减小;从锚头到锚 端,受拉峰值先增加后减小,并且与迎爆侧受拉峰值 位置相比,更远离锚头。 4 迎爆侧除直墙脚锚杆只受拉外,从拱顶到 直墙腰处的锚杆都是先受压后受拉;背爆侧锚杆主 要表现为受拉应力作用。在迎爆侧,从拱顶到直墙 脚,锚杆受拉强度整体上呈现减小趋势,最大受拉峰 值在拱顶;背爆侧受拉情况从拱腰到直墙脚变化逐 渐减小。 [ 2 ] [ 2 ] [ 3 ] [ 4 ] [ 5 ] [ 6 ] [ 7 ] [ 8 ] [ 9 ] [ 1 0 ] [ 1 2 ] 参考文献 R e f e r e n c e s [ 1 2 ] 漆泰岳.锚杆与围岩相互作用的数值模拟[ M ] .徐州 中国矿业大学出版社,2 0 1 2 . 方从严,卓家寿.锚杆加固机理的试验研究现状[ J ] . 河海大学学报 自然科学版 ,2 0 0 5 ,2 3 6 6 9 6 7 0 0 . F A N GC o n g y a n ,Z H U OJ i a s h o u .E x p e r i m e n t a ls t u d ya c t u a l i t yo fr o c kb o l t sr e i n f o r c e m e n tm e c h a n i s m 『J ] .J o u r n a l o fH o h a iU n i v e r s i t y N a t u r a lS c i e n c e ,2 0 0 5 ,3 3 6 6 9 6 7 0 0 . i nC h i n e s e T A N N A N T D D ,B R U M M E RRK ,Y IX .R o c k b o hb e h a v i o r u n d e rd y n a m i cl o a d i n gf i e l dt e s t sa n dm o d e l i n g [ J ] .I n t e r - n a t i o n a lJ o u r n a lo fR o c kM e c h a n i c sa n dM i n i n gS c i e n c e s a n dG e o m e e h a n i c sA b s t r a c t s ,1 9 9 5 ,3 2 6 5 3 7 5 5 0 . A N D E R SA .D y n a m i ct e s t i n go fs t e e lf o ran e wt y p eo fe n e r g ya b s o r b i n gr o c kb o l t [ J ] .J o u r n a lo f C o n s t r u c t i o n a l S t e e lR e s e a r c h ,2 0 0 6 ,6 2 5 5 0 1 - 5 1 2 . A N D E R SA .L a b o r a t o r yt e s t i n go fan e wt y p eo fe n e r g y a b s o r b i n gr o c kb o l t [ J ] .T u n n e l i n ga n dU n d e r g r o u n d S p a c eT e c h n o l o g y ,2 0 0 5 ,2 0 4 2 9 1 3 0 0 . S T P I E R R EL .H A S S A N IFP ,R A D Z I S Z E W S K IPH ,e t a 1 .D e v e l o p m e n to fad y n a m i cm o d e lf o rac o n eb o l t [ J ] . I n t e r n a t i o n a lJ o u r n a lo fR o c kM e c h a n i c sa n dM i n i n gS c i e n c e s ,2 0 0 9 ,4 6 1 1 0 7 11 4 . C H U N L I NL iC .An e we n e r g y a b s o r b i n gb o l t f o rr o c k s u p p o r ti nh i g hs t r e s sr o c km a s s e s [ J ] .I n t e r n a t i o n a lJ o u r n a l o fR o c kM e c h a n i c s a n dM i n i n gS c i e n c e s ,2 0 1 0 , 4 7 3 3 9 6 - 4 0 4 . L ICC ,D O U C E TC .P e r f o r m a n c eo fd - b o l t su n d e rd y n a m i cl o a d i n g [ J ] .R o c kM e c h a n i c sa n dR o c kE n g i n e e r i n g , 2 0 1 2 ,4 5 2 1 9 3 2 0 4 . W A N GC ,H EW ,N I N GJ ,e ta 1 .P r o p a g a t i o np r o p e r t i e so f [ 1 3 ] [ 1 4 ] [ 1 4 ] [ 1 5 ] [ 1 5 ] [ 1 6 ] [ 1 6 ] [ 1 7 ] g u i d e dw a v ei nt h ea n c h o r a g es t r u c t u r eo fr o c kb o h s [ J ] . J o u r n a lo fA p p l i e dG e o p h y s i c s ,2 0 0 9 ,6 9 3 1 3 1 - 1 3 9 . I V A N O V IA ,N E I L S O NRD .I n f l u e n c eo fg e o m e t r ya n d m a t e r i a lp r o p e r t i e so nt h ea x i a lv i b r a t i o no far o c kb o l t JI .I n t e m a t i o n a lJ o u r n a lo fR o c kM e c h a n i c sa n dM i n - i n gS c i e n c e s ,2 0 0 8 ,4 5 6 9 4 1 - 9 5 1 . 薛亚东,张世平,康天合.回采巷道锚杆动载响应的 数值分析[ J ] .岩石力学与工程学报,2 0 0 3 ,2 2 1 1 1 9 0 3 .1 9 0 6 . X U EY a d o n g ,Z H A N GS h i p i n g ,K A N GT i a n h e .N u m e r i e a la n a l y s i so fd y n a m i cr e s p o n s eo fr o c kb o l t si n m i n i n gr o a d w a y s [ J ] .C h i n e s eJ o u m a lo fR o c kM e c h a n i c sa n dE n 舀n e e f i n g ,2 0 0 3 ,2 2 11 1 9 0 3 1 9 0 6 . i nC h i n e s e 荣耀,许锡宾,赵明阶,等.锚杆对应力波传播影响 的有限元分析[ J ] .地下空间与工程学报,2 0 0 6 , 2 1 1 1 5 1 1 9 . R O N GY a o ,X UX i b i n g ,Z H A OM i n g - j i e ,e ta 1 .N u m e r i c a ls i m u l a t i n gp r o p a g a t i o no fs t e s sw a v e si nt h ec o u n t r y r o c k w i t hb o l t i n g [ J ] .C h i n e s eJ o u r n a lo fU n d e r g r o u n d S p a c ea n dE n g i n e e r i n g ,2 0 0 6 ,2 1 1 1 5 1 1 9 . i nC h i n e s e Z H A N GCS ,Z O UDH ,M A D E N G AV .N u m e r i c a ls i m u l a t i o no fw a v ep r o p a g a t i o ni ng r o u t e dr o c kb e l t sa n dt h e e f f e c t so fm e s hd e n s i t ya n dw a v ef r e q u e n e y [ J ] .I n t e r n a t i o n a lJ o u r n a lo fR o c kM e c h a n i c sa n dM i n i n gS c i e n c e s , 2 0 0 6 ,4 3 4 6 3 4 6 3 9 . 余永强,顾金才,杨小林,等.动载下洞室加固锚杆受 力的实验研究[ J ] .兵工学报,2 0 0 9 ,3 0 S 2 2 6 3 2 6 6 . Y UY o n g q i a n g ,G UJ i n c a i ,Y A N GX i a o l i n ,e ta 1 .E x p e r i m e n t a ls t u d yo na p p l i e df o r c eo fr e i n f o r c e m e n tc a v i t y r o c kb o l t su n d e rd y n a m i cl o a d i n gc o n d i t i o n s [ J ] .A c t a A r m a m e n t a r i i ,2 0 0 9 ,3 0 S 2 2 6 3 - 2 6 6 . i nC h i n e s e 王光勇,顾金才,陈安敏,等.拱顶端部加密锚杆支护 洞室抗爆加固效果模型试验研究[ J ] .岩土工程学 报,2 0 0 9 ,3 1 3 3 7 8 3 8 3 . W A N GG u a n g y o n g ,G UJ i n c a i ,C H E NA n m i n ,e ta 1 . M o d e lt e s t so na n t i e x p l o s i o na n c h o r i n ge f f e c to ft u n n e l s r e i n f o r c e db yd e n s eb o h sa ta r c ht o p [ J ] .C h i n e s eJ o u r n a lo fG e o t e c h n i e a lE n g i n e e r i n g ,2 0 0 9 ,3 1 3 3 7 8 3 8 3 . i nC h i n e s e 杨自友,顾金才,杨本水,等.锚杆对围岩的加固效果 和动载响应的数值分析[ J ] .岩土力学,2 0 0 9 ,3 0 9 2 8 0 5 - 2 8 0 9 . Y A N GZ i y o u ,G UJ i n c a i ,Y A N GB e n s h u i ,e ta 1 .N u m e r i c a la n a l y s i so fr e i n f o r c e m e n te f f e c t sa n dr e s p o n s et o d y n a m i cl o a d sc h a r a c t e r i s t i c so fr o c kb o l t s 『J ] .R o c ka n d S o i lM e c h a n i c s ,2 0 0 9 ,3 0 9 2 8 0 5 - 2 8 0 9 . i nC h i n e s e 王光勇,顾金才,陈安敏,等.顶爆作用下锚杆破坏形 万方数据 3 8 爆破2 0 1 5 年6 月 上接第1 6 页 2 理论计算得出水炮泥的填塞长度为 7 .7 5c m ,试验值与理论计算值相近,说明水炮泥填 塞试验模型具有一定的合理性。 3 由于模型试验时,有边界效应对试验结果 的影响以及操作误差的存在,使得试验结果存在一 定的误差。 [ 1 ] [ 2 I [ 2 ] [ 3 ] [ 3 ] [ 4 ] [ 4 ] [ 5 ] [ 5 ] 参考文献 R e f e r e n c e s 李夕兵.凿岩爆破工程[ M1 .长沙中南大学出版社, 2 0 1 1 . 曾新枝.矿岩爆破效果综合评价体系研究与实现[ D ] . 武汉武汉理工大学,2 0 1 2 . Z E N GX i n z h i .T h es t u d ya n dr e a l i z a t i o no ft h eo r ea n d r o c kb l a s t i n ge f f e c tc o m p r e h e n s i v ee v a l u a t i o ns y s t e m [ D ] . W u h a n W u h a nU n i v e r s i t yo fT e c h n o l o g y ,2 0 1 2 . i nC h i n e s e 孙英翔.露天深孔爆破炮孔填塞长度及填塞物运动规 律的研究[ D ] .沈阳东北大学,2 0 1 1 . S U NY i n g x i a n g .S t u d yo fb l a s t - h o l es t e m m i n gl e n g t ha n d s t e m m i n gm o v e m e n tr u l ei nd e e p - h o l eb l a s t i n g [ D ] .S h e n y a n g N o r t h e a s t e r nU n i v e r s i t y ,2 0 1 1 . i nC h i n e s e 李向东,孙萌苑.新型水炮泥降低爆破烟尘的试验 [ J ] .煤炭科学技术,2 0 1 1 ,3 9 1 5 3 - 5 6 . L IX i a n g d o n g ,S U NY i n y u a n .E x p e r i m e n to nn e ww a t e r s t e m m i n gt or e d u c eb l a s t i n gs m o k ea n dd u s t [ J ] .C o a lS c i e n c ea n dT e c h n o l o g y ,2 0 1 1 ,3 9 1 5