围压卸载速率对砂岩力学性质影响数值试验研究.pdf

第5 4卷 第4期2 0 2 3年7月 太原理工大学学报J OUR NA L O F TA I YUAN UN I V E R S I T Y O F T E CHNO L OG Y V o l.5 4 N o.4 J u l.2 0 2 3 引文格式:张健,郭亮亮.围压卸载速率对砂岩力学性质影响数值试验研究J.太原理工大学学报,2 0 2 3,5 4(4):7 0 0-7 0 9.Z HAN G J i a n,GUO L i a n g l i a n g.N u m e r i c a l e x p e r i m e n t a l i n v e s t i g a t i o n s o f t h e e f f e c t o f c o n f i n i n g p r e s s u r e u n l o a-d i n g r a t e o n s a n d s t o n e m e c h a n i c a l p r o p e r t i e sJ.J o u r n a l o f T a i y u a n U n i v e r s i t y o f T e c h n o l o g y,2 0 2 3,5 4(4):7 0 0-7 0 9.收稿日期:2 0 2 2-0 9-1 6;修回日期:2 0 2 3-0 3-1 0 基金项目:中国博士后科学基金资助项目(2 0 2 0 T 1 3 0 3 9 0,2 0 1 9 M 6 6 1 0 5 3);山西省基础研究计划资助项目(2 0 2 1 0 3 0 2 1 2 4 4 8 5)第一作者:张健(1 9 8 6-),博士,讲师,主要从事矿井水文地质与工程地质方面的研究,(E-m a i l)5 0 c e n t 0 8 1 01 6 3.c o m 通信作者:郭亮亮(1 9 8 8-),博士,副教授,主要从事水文地质与工程地质方面的研究,(E-m a i l)7 9 2 0 8 7 6 8 5q q.c o m围压卸载速率对砂岩力学性质影响数值试验研究张 健1,2,郭亮亮2(1.山西省地质调查院有限公司,太原 0 3 0 0 0 6;2.太原理工大学 水利科学与工程学院,太原 0 3 0 0 2 4)摘 要:【目的】为了排除天然岩石试件离散性与试验过程多种应力路径相互转换对试验结果的叠加影响,进一步探讨围压卸载速率对岩石破坏失稳的影响机理。【方法】利用离散元法开展均质数值试件在恒定轴压与两种固定围压卸载速率条件的数值试验研究,对比与讨论了两种围压卸载速率下的试件宏观与细观力学特性与规律。【结果】结果表明:两种条件下,宏观和细观破裂分布具有一定相似性和规律性,慢速围压卸荷条件较快速条件的试件承载强度更高,快速卸荷条件局部破裂分布密度相对更大。两种卸载速率试验过程出现轴压和围压的“跌落”和“上升”,反映试件经历了多阶段突发式渐进破坏过程,因局部破坏不断累积导致试件向整体失稳方向发展。两种卸载速率条件下的张拉破裂数量显著大于剪切破裂数量,张拉破裂增长速率大于剪切破裂,试件表现为张拉破坏为主,快速卸载条件破裂总数量相对更大。慢速卸载条件相对充分地调动了试件局部承载能力,使试件极限承载能力有所提高。【结论】两种条件下,2 0次以上声发射事件构成的裂纹分布程度最不均一,随着组成每条裂纹的声发射事件数增加,震级整体逐渐增大,快速卸荷形成裂纹过程总体释放能量相对更高,而慢速卸荷以相对较慢的速率传递应力,使裂纹形成过程释放能量相对更低、平稳且均一。关键词:围压卸载速率;数值试验;力学特性;散元法;砂岩中图分类号:TU 4 5 文献标识码:AD O I:1 0.1 6 3 5 5/j.c n k i.i s s n 1 0 0 7-9 4 3 2 t y u t.2 0 2 3.0 4.0 1 4 文章编号:1 0 0 7-9 4 3 2(2 0 2 3)0 4-0 7 0 0-1 0N u m e r i c a l E x p e r i m e n t a l I n v e s t i g a t i o n s o f t h e E f f e c t o f C o n f i n i n g P r e s s u r e U n l o a d i n g R a t e o n S a n d s t o n e M e c h a n i c a l P r o p e r t i e sZ H A N G J i a n1,2,G U O L i a n g l i a n g2(1.S h a n x i I n s t i t u t e o f G e o l o g i c a l S u r v e y C o.,L t d.,T a i y u a n 0 3 0 0 0 6,C h i n a;2.C o l l e g e o f W a t e r R e s o u r c e s S c i e n c e a n d E n g i n e e r i n g,T a i y u a n U n i v e r s i t y o f T e c h n o l o g y,T a i y u a n 0 3 0 0 2 4,C h i n a)A b s t r a c t:【P u r p o s e s】T h i s w o r k i s c o n d u c t e d t o e x c l u d e t h e s u p e r i m p o s e d e f f e c t o f n a t u r a l r o c k s p e c i m e n d i s c r e t e n e s s a n d m u l t i p l e s t r e s s p a t h c o n v e r s i o n i n t e s t p r o c e s s o n r e s u l t s,a n d t o f u r t h e r i n v e s t i g a t e t h e i n f l u e n c e m e c h a n i s m o f t h e c o n f i n i n g p r e s s u r e u n l o a d i n g r a t e o n r o c k f a i l-u r e a n d i n s t a b i l i t y.【M e t h o d s】T h e d i s c r e t e e l e m e n t m e t h o d w a s u s e d t o c a r r y o u t h o m o g e n e o u s n u m e r i c a l s p e c i m e n s u n d e r c o n s t a n t a x i a l p r e s s u r e a n d t w o f i x e d c o n f i n i n g p r e s s u r e u n l o a d i n g r a t e s.T h e m a c r o s c o p i c a n d m e s o s c o p i c m e c h a n i c a l p r o p e r t i e s a n d l a w s o f t h e s p e c i m e n s u n d e r t w o c o n f i n i n g p r e s s u r e u n l o a d i n g r a t e s a r e c o m p a r e d a n d d i s c u s s e d.【F i n d i n g s】T h e r e s u l t s s h o w t h a t t h e m a c r o s c o p i c a n d m e s o s c o p i c f r a c t u r e d i s t r i b u t i o n s h a v e c e r t a i n s i m i l a r i t y a n d r e g u l a r i t y u n d e r t h e t w o c o n d i t i o n s,t h e l o a d-b e a r i n g s t r e n g t h o f t h e s p e c i m e n i n s l o w c o n f i n i n g p r e s s u r e u n l o a d i n g r a t e i s h i g h e r t h a n t h a t i n f a s t c o n d i t i o n,a n d t h e l o c a l f r a c t u r e d i s t r i b u t i o n d e n s i t y a t t h e f a s t u n l o a d i n g c o n d i t i o n i s r e l a t i v e l y l a r g e r.T h e a x i a l a n d c o n f i n i n g p r e s s u r e“f a l l s”a n d“r i-s e s”d u r i n g t h e t w o u n l o a d i n g r a t e t e s t s,w h i c h r e f l e c t s t h a t t h e s p e c i m e n s h a v e u n d e r g o n e a m u l t i-s t a g e s u d d e n a n d p r o g r e s s i v e f a i l u r e p r o c e s s.T h e c o n t i n u o u s a c c u m u l a t i o n o f l o c a l d a m a g e l e a d s t o t h e o v e r a l l i n s t a b i l i t y o f t h e s p e c i m e n s.T h e n u m b e r o f t e n s i l e f r a c t u r e s u n d e r t h e t w o u n l o a d i n g r a t e s i s s i g n i f i c a n t l y l a r g e r t h a n t h a t o f s h e a r f r a c t u r e s,a n d t h e g r o w t h r a t e o f t e n s i l e f r a c t u r e s i s g r e a t e r t h a n t h a t o f s h e a r f r a c t u r e s.T h e s l o w u n l o a d i n g c o n d i t i o n r e l a t i v e l y f u l l y m o-b i l i z e d t h e l o c a l b e a r i n g c a p a c i t y o f t h e s p e c i m e n a n d i m p r o v e d t h e u l t i m a t e b e a r i n g c a p a c i t y o f t h e s p e c i m e n.【C o n c l u s i o n s】U n d e r t h e t w o c o n d i t i o n s,t h e d i s t r i b u t i o n o f f r a c t u r e s c o m p o s e d o f m o r e t h a n 2 0 A E e v e n t s i s t h e m o s t u n e v e n.W i t h t h e i n c r e a s e o f A E e v e n t n u m b e r s c r e a t i n g e a c h f r a c t u r e,t h e o v e r a l l s e i s m i c m a g n i t u d e i n c r e a s e s g r a d u a l l y,a n d t h e o v e r a l l e n e r g y r e l e a s e d u r i n g t h e f a s t u n l o a d i n g p r o c e s s t o f o r m f r a c t u r e s b e c o m e s r e l a t i v e l y h i g h e r,w h i l e s l o w u n l o a-d i n g t r a n s f e r s s t r e s s a t a r e l a t i v e l y s l o w r a t e,r e s u l t i n g i n a r e l a t i v e l y l o w e r,s m o o t h s,a n d u n i-f o r m e n e r g y r e l e a s e d u r i n g f r a c t u r e f o r m a t i o n.K e y w o r d s:c o n f i n i n g p r e s s u r e u n l o a d i n g r a t e;n u m e r i c a l e x p e r i m e n t;m e c h a n i c a l p r o p e r t i e s;d i s c r e t e e l e m e n t m e t h o d;s a n d s t o n e 常见的深部地下工程灾害(如岩爆、煤与瓦斯突出等动力灾害)与开挖诱发的围压卸荷密切相关,岩体卸荷破坏所引发的灾变越来越受到科研和工程技术人员的关注。在煤体或岩石开挖过程中,掌子面应力状态可简化为轴压相对稳定而围压卸载。如何设置合理的工作面开挖速率以控制灾变程度,是工程实践中面临的现实问题,其背后的科学问题是围压卸载速率对岩石损伤破坏的影响机理。由于岩石材料自身的特性,加荷与卸荷条件力学响应与机理具有差异。国外学者曾利用三轴仪研究围压卸荷对岩石强度的影响,并提出卸荷应力路径对岩石强度不具影响的结论1-2。随着实验仪器的发展,研究人员通过现场和室内试验发现,应力路径显著影响着岩石强度、变形和破裂模式,卸荷条件相对加荷条件下,岩石破坏表现更加强烈,沿卸荷方向扩容现象也更加明显3-5。还有学者提出,节理岩体卸荷损伤破坏是由于裂隙非稳定扩展导致的脆性断裂,但相对加荷条件,其更具有突发性6。但上述成果并未对比 卸荷速率 相对大 小 对 试 验 结 果 的影响。不同开挖速率引起卸荷过程中岩石力学特性的表现不同7-8,卸载速率快慢显著影响着岩石物理力学特性9-1 3。人们注意到,通过调整开挖速率并辅助加固措施,可有效控制由开挖卸荷引起的工程灾害威胁1 4-1 5,如高地应力环境下通过调整开挖速率在一定程度上控制岩爆发生的频率与烈度。黄润秋等1 6利用大理石试件建立了卸荷破坏损伤变量与卸荷速率的拟合关系,提出在一定围压卸载速率范围内,卸荷破坏损伤变量存在最小值,围压卸荷过程岩石破裂以侧向张拉破坏为主,岩石破坏是强烈扩容所致。借助声发射监测仪器,学者们从能量转换的角度开展了不同卸载速率下的真三轴试验,发现声发 射 能 量 释 放、岩 爆 程 度 与 卸 载 速 率 关 系 密切1 7,卸荷应力路径显著影响着岩石破裂模式、能量存储与转换过程1 8-1 9。学者们关于岩体卸荷对其自身强度的影响,提出了不同看法。由于当时实验条件的限制,曾有学者提出卸荷应力路径对岩石强度的影响并不明显2 0-2 1。而吴刚等2 2-2 3认为岩体强度随卸荷速率的提高而不断增加。文献9-1 0 利用大理岩试件开展了不同围压卸载速率下的室内试验研究,对多组不同初始围压条件下的测试结果进行了拟合,同样认为,随着围压卸载速率的提高,试件极限承载强度也随之不断增加。然而容易被忽略的是,上述成果虽考虑到天然试件的离散性对试验结果的影响,并尝试将该影响控制在一定范围内,但仍然无法完全排除这种干扰。学者们结合工程现场实际,考虑了包含卸荷在内的多类型应力路径条件对岩体力学特性的影响。S I e t a l2 4使用钻进装置模拟了天然应力状态下隧道施工过程,将其与预制钻孔试件经过加载的条件107 第4期 张 健,等:围压卸载速率对砂岩力学性质影响数值试验研究进行了对比,认为预制钻孔后种试验条件出现的试件强度强化现象,比天然条件引起的试件强度软化现象更加明显。WANG e t a l2 5通过对经疲劳加载后的试件进行围压卸载试验,提出疲劳加载引起的初始损伤积累影响着试件破裂模式。雷华阳等2 6详细探讨了开挖面失稳过程中的土拱效应,提出基于土体三向应力时空变化的开挖面失稳破坏机理。但上述成果综合考虑了试件多阶段受载过程的叠加效果,故无法辨别其中卸载阶段对试件的单独影响机理。数值试验方面,YANG e t a l2 7通过有限元模型,分析了预制裂隙试件突发卸荷条件的破裂过程和类型,发现预制裂隙倾角与长度对试验结果影响非常明显。L I U e t a l2 8借助有限离散元法,建立了实验室和场尺度两种数值模型,分析了加载和卸载速率对于裂纹扩展的影响,并分别给出了两种模型网格单元尺寸选择建议。李竹等2 9通过精细化模拟巷道掘进过程围岩加载和卸载应力路径变化,归纳了细砂岩损伤和破裂规律。然而,现有数值试验中对单纯围压卸荷引发试件破坏失稳的细观机理尚未探明。综上,现有关于围压卸荷速率对岩石物理力学特性和破坏失稳影响机理研究程度不足,仍有待继续深入探讨。本文借助离散元法开展了恒定轴压与两种固定速率围压卸载条件下的数值试验研究,研究目的是对比与分析卸荷速率相对快慢对数值试件强度、变形和破裂特征的影响机理,旨在对卸围压条件下岩体破坏失稳机制分析提供理论和工程实际参考。1 试验设置1.1 模型构建与参数标定若模型颗粒越小,则计算效果越佳,但笔者前期对球体颗粒直径范围为1.01.5 mm的三维数值圆柱型标准试件进行了试算,模型平衡过程与后续三维加载过程消耗的计算时间成本较高。由于建立的是三维数值模型,为了平衡计算精度和时间,经过反复试算调整,并参考已有文献3 0-3 2,使颗粒最大粒径比真实试件尺寸小一个数量级。确定的模型粒径区间为1.54.0 mm,颗粒数量为1 3 7 8 8,内尺度比数量级为1 0-2,符合上述文献的建议范围内。由于研究区现场取样工作一直未能达成,但研究区与文献3 3-3 4 处于同一区域,且生产条件较为相似。文献中记录了砂岩不同围压条件下的三轴压缩实验,以及巴西劈裂实验结果。参考研究区深部区域地应力水平以及围压3为1 5 MP a的条件,试件峰 值 轴 向 应 力7 4.4 2 MP a,对 应 轴 向 应 变 为0.0 0 5 1.利用上述结果对数值模型细观参数进行标定,标定后的模型与文献中试验结果吻合,细观参数合理,限于篇幅,文中不再赘述文献结果。模型接触类型选取平行黏结模型。标定后的细观参数如表1所示。标定后的数值模型在三轴压缩条件下应力-应变曲线特征,及侧表面破裂分布情况如图1所示。表1 数值模型细观参数T a b l e 1 M e s o s c o p i c p a r a m e t e r s o f t h e n u m e r i c a l m o d e l数值模型细观参数赋值密度/(k gm-3)2 6 0 0孔隙度0.0 8粒径范围/mm1.54.0平行黏结模量/G P a2 2.6平行黏结刚度比3法向黏结强度/MP a4 0.3 8切向黏结强度/MP a3 5.1 1摩擦系数0.6 8 7内摩擦角/()3 4.5 4剪切破裂张拉破裂80706050403020100（13）/?MPa0.00100.002 0.003 0.004 0.005 0.006 0.007 0.0081图1 校正后的数值模型与破裂分布F i g.1 C a l i b r a t e d n u m e r i c a l m o d e l a n d f r a c t u r e d i s t r i b u t i o n1.2 试验方案本文旨在对比围压卸荷速率差异对工程岩体损伤破坏机理的影响,分析同一数量级围压卸载速率条件下,相对快速和相对慢速两种结果的差异,并非针对某一特定实际工况。文中设置两种围压卸荷速率的依据主要是计算时间成本和工程施工实际。首先,为了保证计算精度,本文设置的数值计算时间步长为1 0-7秒/步(s/s t e p),当程序计算完成1 07步(s t e p)仅相当于现实中的1 s,耗时巨大。若要实现数值试验与现实中的围压卸载速率完全对等,需要超大型计算机的硬件支撑,而研究团队目前暂不具备该条件,尽管如此,团队现有计算能力仍然可以对表象内部的规律与机理进行探究。通过试算,如果设置卸荷速率数量级为1 01 P a/s t e p,计算时间太长难以承受,如遇计算机断电或死机等异常情况,则前功尽弃;如果设置的207太 原 理 工 大 学 学 报 第5 4卷 数量级为1 03 P a/s t e p,通过5组相对快慢两种速率的试算结果显示,差异并不明显,因此选择的数量级为1 02 P a/s t e p,便于后文数据分析。其次,本文研究目的是在高地应力条件下,施工速率在一个数量级变化范围内,影响工作面岩体损伤破坏机理的差异,通过设置相对快速和相对慢速两种条件,前者是后者卸荷速率的约1.5倍。在深部地下工程开挖中,如煤矿巷道掘进施工,在地质条件较为稳定的区块内,掘进施工速率亦保持相对稳定,速率快慢相差1.5倍是相对合理的比值,能够在一定程度上反映工程实际当中高、低两种围压卸载速率。因此本文考虑工程实际,兼顾平衡计算时间和精度,经过反复试算确定了两种卸荷速率:1 7 8.5 7 P a/s t e p和1 1 3.6 4 P a/s t e p,分别代表快速和慢速围压卸荷速率,计算时间步长数量级设置为1 0-7以保证计算精度。两种速率条件均采用相同细观参数的均质模型计算,如前文所述,如果试验过程轴压显著变化,则无法客观反映卸围压速率快慢对试件变形与破坏的影响,故试验过程保持恒定轴压,以两种固定速度卸载围压。初始轴压的设置需超过试件弹性极限值,同时考虑试件卸荷过程变形破坏的可观测性,初始轴压设置过高或过低都难以达到试验目的,经过试算,设置初始轴压为围压1 5 MP a三轴压缩条件下峰值轴压的8 0%,初始围压1 5 MP a,分别以两种固定速率卸载围压。2 结果与讨论2.1 初始损伤分布特征如图2所示,试件达到预设轴压与围压后的初始破裂分布情况,主要分布在试件的顶面和底面,侧表面相对较少,共计6 7 4个破裂。两种围压卸荷速率条件都需加载至相同轴压与围压,形成了相同的试件初始损伤,但从围压卸荷开始后,其破裂分布与类型各有特征。剪切破裂张拉破裂图2 数值试件初始损伤分布F i g.2 I n i t i a l d a m a g e d i s t r i b u t i o n o f t h e n u m e r i c a l s p e c i m e n2.2 黏结破裂分布特征如图3所示,第一行和第二行分别代表快速和慢速条件破坏情况,每一行从左到右依次展示试件顶面、侧表面和底面,侧表面的前一幅视角沿试件轴向顺时针旋转9 0 即为后一幅视角,后一幅视角可观察到前一幅视角右半侧图像,第4幅视角顺时针旋转9 0 即 得 到 第1视 角。颗 粒 设 置 透 明 度 为9 0%,便于观察试件内部一定范围破裂分布情况。剪切破裂张拉破裂图3 两种围压卸荷速率下破裂分布特征F i g.3 B o n d-b r e a k i n g d i s t r i b u t i o n c h a r a c t e r i s t i c s o f t w o c o n f i n i n g p r e s s u r e u n l o a d i n g r a t e s两种速率条件的试件破裂位置分布具有一定差异性和相似性。对于顶面(如红色虚线圈所示),快速较慢速条件的局部破裂密度更大,而黑色虚线圈内前者比后者破裂密度要小。对于底面(如黄色虚线圈所示),两种速率的局部破裂类型、位置和数量几乎相同,该现象是由于试件加载至预设条件后,形成的相同初始损伤所致。在侧表面上,快速卸荷条件形成集中破裂分布面积显著大于慢速条件,且出现了一条横贯试件的明显破坏带;而慢速卸荷条件在相同范围内形成的破坏带内的破裂分布密度明显减小且较为分散。快速和慢速条件的总黏结破裂数量分别为5 9 1 4和5 3 2 0,可以看出,前者较后者的破裂分布相对更加集中,是由大量局部破裂连通导致试件破坏,而后者是由于试件内部不同范围相对分散的局部破裂彼此连通造成整体失稳破坏。2.3 强度与变形特征2.3.1 应力-应变曲线图4(a)和(b)分别展示了试件加载至预设轴压和围压,然后保持轴压稳定,围压按照两种固定速率卸载至峰后阶段过程的应力-应变曲线。两条曲线在围压卸载前阶段完全一致,区别在于黑色箭头后307 第4期 张 健,等:围压卸载速率对砂岩力学性质影响数值试验研究的卸荷阶段。6050403020100（13）/?（MPa178.57?Pa/step-0.00100.0010.0020.0030.00413（a）快速卸荷应力-应变曲线6050403020100（13）/?MPa113.64?Pa/step-0.00100.0010.0020.0030.00413（b）慢速卸荷应力-应变曲线图4 两种速率卸荷应力-应变曲线F i g.4 S t r e s s-s t r a i n c u r v e s o f t w o c o n f i n i n g p r e s s u r e u n l o a d i n g r a t e s图5(a)和(b)分别展示了两种速率围压卸载过程的应力-应变曲线,二者峰值强度具有区别。但在围压开始卸荷后,二者曲线具有显著区别。快速卸荷条件峰值偏应力约为5 7.5 0 MP a,而慢速条件下约为5 7.6 5 MP a,峰值偏应力有所增加,试件极限承载强度得以提高,结果与文献9,2 2-2 3 不同。两种卸载速率下接触模型类型和细观参数完全相同,但其峰值强度却有一定差异,其原因在于不同围压卸载速率对于试件内部接触元件驱动机制具有影响,随围压卸载速率降低,试件颗粒之间的接触应力,在相对更长的时间步中相对更加充分地传递、调整,由局部承载逐渐过渡为整体承载,这一过程能够使颗粒有愈来愈长的时间参与应力调整与重新分布,从而提高了试件的极限承载强度。另外,真实岩石试件中天然存在结构面,不同部位材料强度具有差异,导致真实试件承载力在内部并不均一,但数值试件能够做到承载力完全均一,且数值试验卸荷速率高于真实物理试验,造成了试验结果与文献9,2 2-2 3 不同。2.3.2 轴压与围压变化特征在数值试验过程中,预加载后轴压恒定而围压5857565554（13）/?MPa178.57?Pa/step-0.00100.0010.0020.0030.00413（a）快速卸荷开始后应力-应变曲线5857565554（13）/?MPa113.64?Pa/step-0.00100.0010.0020.0030.00413（b）慢速卸荷开始后应力-应变曲线图5 两种速率卸荷开始后应力-应变曲线特征F i g.5 S t r e s s-s t r a i n c u r v e s o f c o n f i n i n g p r e s s u r e u n l o a d i n g p r o c e s s a t t w o r a t e s卸载,则偏应力不断增加,那么试件将发生破裂,从而引起试件轴向承载力的变化,这样就必然造成轴压波动,但可通过轴向应力控制以满足相对恒定轴压条件,如设置轴压波动范围小于固定值的1%,而围压则按照设计固定速率卸载,这一过程是数值程序控制实现的,数值试验的伺服原理与真实室内试验的伺服程序在本质上并无差别,不同之处是数值试件与真实试件的材料结构无法完全相同。图6展示了围压卸载开始至达到试件峰值强度9 5%停止试验的过程中,轴压与围压随时间的波动特征。停止试验的依据在下文中进行解释。两种速率下的轴压值均保持稳定,而围压波动有所差异,从围压波动次数看,试验结束前快速条件为4次,而慢速条件为5次;从围压波动发生时间看,快速条件所发生的4次波动均早于慢速条件;从围压波动的平均幅度看,快速条件较慢速条件相对更大。试验结束时,快速较慢速条件的围压值大。可以看出,慢速条件在以一种相对平稳的方式进行卸荷,慢速条件的围压卸载程度相对更高,反映出慢速条件下试件颗粒间作用相对更加充分。两种卸载速率下的共同特点是,轴压和围压的波动方向相反,轴压出现的小幅跌落而同时使得围压反弹上升,这407太 原 理 工 大 学 学 报 第5 4卷 是由于试件突发破裂轴向承载力瞬时减小,而为了保持轴压恒定,围压瞬时增大以维持轴压相对稳定,轴压在围压卸荷过程始终处于动态调整状态,这也验证了伺服程序的正确性。慢速条件1 1 3.6 4 P a/s t e p的试验时间只有约0.1 3 s,而在笔者的数值计算过程中却消耗了接近2 4 h,其时间成本可见一斑。正如前文所述,由于时间成本高昂,数值试验的围压卸载速率无法与真实试验完全对等。如果将数值计算的卸载速率设置的非常小,那么在足够长的试验时间内围压反弹的频率和幅度将会被平均“抹平”,这样就使得围压与时间呈现相对意义上的线性关系。（a）慢速卸荷过程轴压与围压随时间变化曲线7060504030201001?/?MPa0.020.040.060.080.10时间/?s7060504030201003?/?MPa13（a）快速卸荷过程轴压与围压随时间变化曲线7060504030201001?/?MPa0.020.040.060.080.10时间/?s7060504030201003?/?MPa130.120.14图6 两种速率下轴压与围压随时间变化曲线特征F i g.6 A x i a l p r e s s u r e a n d c o n f i n i n g p r e s s u r e c h a n g e s w i t h t i m e o f u n l o a d i n g p r o c e s s a t t w o r a t e s两种条件试验过程中轴压和围压的波动,表明试件经历了多阶段突发性渐进破裂过程,试件脆性特点表现明显。如果延长峰后阶段试验时间,例如持续至峰值应力7 5%停止试验,那么在数值程序的控制下,随着试验的进行二者波动幅度也会越来越大,这样难以满足恒轴压和固定速率卸围压试验条件,故试验达到峰值强度9 5%即停止,这充分反映两种卸载速率条件下对应特征的变化规律。2.4 细观损伤特征2.4.1 微裂纹数量变化特征图7为两种速率卸荷过程黏结破裂类型和数量随试验时间变化规律。6?0005?0004?0003?0002?0001?0000破裂数量破裂总数张拉破裂剪切破裂0.020.040.060.080.100.120.14时间/?s178.57?Pa/step（a）快速卸荷过程破裂类型和数量变化曲线6?0005?0004?0003?0002?0001?0000破裂数量破裂总数张拉破裂剪切破裂0.020.040.060.080.100.120.14时间/?s113.64?Pa/step（b）慢速卸荷过程破裂类型和数量变化曲线图7 两种速率下破裂类型和数量变化曲线特征F i g.7 F r a c t u r e n u m b e r s a n d t y p e s o f c o n f i n i n g p r e s s u r e u n l o a d i n g p r o c e s s a t d i f f e r e n t r a t e s两种卸荷速率条件下,围压不断卸载而轴压稳定,使试件径向扩容引起张拉破裂的数量约为剪切破裂数量的2倍多,试件均以张拉破坏为主。快速卸荷条件产生的破裂数量大于慢速条件,破裂增长速率也相对更大。黏结破裂的增长均表现为多阶段“平台”式,在每个“平台”阶段破裂数量基本稳定,当产生突发性变形时,破裂数量瞬时增加。快速条件较慢速条件形成的“平台”数量要小,且平均“平台”高度更大,反映出试件的短时破坏失稳的特点,而慢速条件则为多阶段破裂累积引发失稳。2.4.2 声发射事件分布特征如图8所示,快速卸荷条件较慢速条件的平均单位体积声发射事件分布密度要高,声发射事件密集部位为破裂部位,慢速条件的声发射事件分布相对分散。快速卸荷条件下,试件颗粒依靠整体均一的承载能力迅速调整其内部应力状态并产生局部密集破裂,声发射事件分布较为不均匀;而慢速卸荷条件下,颗粒有相对较长的时间由两端部向内部传递应力,使试件整体相对缓慢地承受荷载,局部拥有相对更长的时间参与应力调整与分配,相对充分地调动了局部承载能力,从而使试件极限承载强度有所提高。2.4.3 单裂纹与震级分布特征根据地震学中的矩张量理论可对震源信息进行507 第4期 张 健,等:围压卸载速率对砂岩力学性质影响数值试验研究（b）慢速卸荷过程声发射事件分布0.040.050.020.010-0.01-0.02-0.03-0.040.040.030.020.010-0.01-0.02-0.03-0.040.040.020-0.02-0.04178.57?Pa/step张拉破裂剪切破裂（a）快速卸荷过程声发射事件分布0.040.050.020.010-0.01-0.02-0.03-0.040.040.030.020.010-0.01-0.02-0.03-0.040.040.020-0.02-0.04113.64?Pa/step张拉破裂剪切破裂长/?m宽/?m高/?m长/?m宽/?m高/?m图8 两种速率卸荷过程声发射事件分布F i g.8 A c o u s t i c e m i s s i o n e v e n t d i s t r i b u t i o n s o f c o n f i n i n g p r e s s u r e u n l o a d i n g p r o c e s s a t d i f f e r e n t r a t e s分析。将两个颗粒间的黏结破坏视为一次微破裂事件,即声发射事件。由于黏结破坏后颗粒位移使接触变形,引起接触力变化。微破裂事件的影响区域以微破裂中心为圆心,影响半径为最大颗粒直径,将颗粒周围所有接触的接触力变化量乘以对应的力臂,然后求和可得到矩张量Mi j3 5:Mi j=S(FiRj).(1)式中:Fi为接触力变化量的第i个分量;Rj为颗粒接触点与微裂纹中心距离的第j个分量;S为颗粒周围接触数量。但在实际中为了提高计算效率,采用矩张量的最大标量力矩值代替单次声发射事件的