Action time sequence of underwater explosion shock waves and shaped charge projectiles
-
摘要: 聚能装药水下爆炸过程中会产生高速聚能侵彻体和强间断冲击波等多种毁伤元。由于聚能侵彻体和冲击波的作用时间接近，且聚能装药水下爆炸作用时序的理论并不完善，因而认识两者的作用时序对聚能型战斗部作用下舰船结构的毁伤研究具有重要意义。首先，基于接触爆炸理论和牛顿第二定律，推导药型罩压垮后加速度和速度公式的基本形式。随后，基于欧拉控制方程，建立聚能装药空中和水下爆炸数值模型，得到装药和药型罩交界面处压力时程曲线，定量地确定药型罩压垮的加速度和速度公式，通过理论公式可解决不同炸高下聚能侵彻体和直达冲击波先后到达目标的问题。为了验证理论公式的可靠性，讨论了空气域长度为5倍装药半径时的复杂工况，数值模拟结果和理论推导结果基本一致：当空气域长度为5倍装药半径时，炸高在3倍装药半径之外，冲击波先于侵彻体。提出了药型罩压垮的加速度和速度理论公式的形式和求解聚能侵彻体和冲击波作用时序问题的思路，为分析聚能装药水下爆炸聚能侵彻体和冲击波的作用时序提供了理论依据。Abstract: Damage elements such as high-speed explosively-formed projectiles and strong discontinuous shock waves are generated during the process of the shaped charge associated with the underwater explosion. The theory on the action time sequence of the explosively-formed projectile and the shock wave should be refined because their action time is close. Therefore, it is of great significance to investigate the action time sequence of different loads and their damage on ship structures. First, the formulations of acceleration and velocity equations are deduced in the forming process of the explosively-formed projectile, based on the contact explosion theory and Newton’s second law. Subsequently, based on the Eulerian governing equations, numerical models of air and water explosions of shaped charges are established. The evolution of pressure at the interaction of the charge and the liner is obtained. The acceleration and velocity equations of the explosively-formed projectile are presented quantitatively in this paper as a result. Besides, the obtained theoretical formulation can be utilized to solve the problem of the action time sequence of the explosively-formed projectile and direct shock wave. In order to verify the reliability of this theoretical formulation, the case in which the air cavity length is five times of the charge radius is studied. The numerical results are in general agreement with those of the theoretical derivation. The results show that when the length of the air cavity is five times of the charge radius and the stand-off distance is greater than three times of the charge radius, the shock wave precedes the explosively-formed projectile. The basic form of theoretical formulas is presented for the acceleration and velocity of the explosively-formed projectile. Moreover, the idea of solving the action time sequence problem of these two loads provides a theoretical basis for analyzing the action time sequence of underwater explosions.
-
中国西南地区大型水电工程高边坡工程规模巨大、地质环境复杂、工程作用强烈，高边坡的变形与稳定控制是高坝工程建设与运行安全的重要保障[1]。钻孔爆破是水电工程岩石高边坡开挖的主要手段，节理岩质边坡在爆破开挖过程中常发生位移突变，影响边坡的稳定[2-4]。预应力锚索作为一种主动支护技术，因其施工方便、见效快和造价低等优点在边坡支护中得到了广泛的应用[5-6]，有效地控制了边坡岩体的位移突变。
对于节理岩质边坡的位移突变问题，国内外学者从卸荷松动和爆炸荷载扰动的角度进行了广泛的研究。Yoshida等[7]、Sheng等[8]通过数值模拟与现场实测数据的对比分析，认为开挖引起边坡岩体内应力释放及重分布是造成边坡表面岩体位移增加、节理张开及滑移的主要原因。卢文波等[9]的研究表明，爆破开挖过程中岩体初始地应力的释放具有瞬态特性，地应力的瞬态释放会导致岩体节理张开；罗忆等[10]采用室内试验验证了地应力瞬态卸载条件下的节理岩体位移突变问题，并认为突变位移与初始地应力的平方成正比。Gran等[11]通过对节理岩体模型进行爆破试验，发现爆炸应力波传播至节理处发生反射会将节理拉开；金李等[12]利用一维应力波理论对节理岩体的爆破松动机理进行了分析；代金豪等[3]则从能量的角度研究了节理岩体的爆破松动问题，认为爆炸荷载作用期间岩体积累的应变能及其释放速率决定了岩体的爆破松动位移。以往的研究多将地应力瞬态卸荷和爆炸荷载扰动效应分开讨论，较少涉及二者耦合作用下节理岩质边坡的位移突变特征。
在工程现场，边坡爆破开挖施工与锚固支护常穿插进行。众多学者在预应力锚索对节理岩质边坡变形控制的机理方面开展了研究，如丁秀丽等[5]提出预应力锚索锚固效应使得岩体力学性能劣化在边坡开挖过程中得到缓解，进而抑制了边坡变形；Yang等[13]通过数值模拟研究发现，锚索预应力的施加可有效改善岩体的应力状态并提高软弱面的抗滑力，从而达到控制边坡变形的目的；李剑等[6]则认为对于高陡岩质边坡，预应力锚索主要在于限制潜在可动滑体在开挖过程中的位移发展。上述有关预应力锚索控制边坡变形的机理研究多针对静态或拟静态力学条件进行。对于深切河谷底部地应力集中区的边坡岩体爆破开挖，存在爆炸荷载与地应力瞬态卸荷2种动力扰动。有关爆炸荷载与地应力瞬态卸荷耦合作用下的锚固节理岩质边坡位移突变机理，目前研究还较为少见。
本文中，首先，对白鹤滩水电站左岸坝基节理岩质边坡爆破开挖过程中的岩体位移及锚索轴力进行现场监测资料分析；然后，建立三维节理岩质边坡数值模型，分析爆炸荷载与地应力瞬态卸荷耦合作用下节理岩质边坡的位移突变特性及其能量机理，并从能量的角度研究预应力锚索对节理岩质边坡位移突变的控制机理。
1. 节理岩质边坡位移突变与预应力锚索轴力监测分析
1.1 工程概况
白鹤滩水电站位于金沙江下游四川省宁南县和云南省巧家县境内，左岸坝基边坡岩体由峨眉山玄武岩构成，其中834～820 m高程的岩体主要为弱风化下段、弱卸荷的Ⅲ2类岩体，820～600 m高程的岩体主要为弱风化下段～微新无卸荷的Ⅲ1类岩体，600 m高程以下的岩体主要是微新无卸荷的Ⅱ类岩体，如图1所示。680 m高程出露的NE向断层F17斜切坝基边坡，在F17上盘岩体中发育有NNW向张扭性断层f108和f110、层间错动带C3和C3-1以及层内错动带LS342、LS411和LS423等。在F17下盘岩体中存在有顺坡向的层内错动带LS3319、LS331和LS3301等，还发育有纵横交错且易松弛的柱状节理和微裂隙。
图 1 白鹤滩水电站左岸坝基边坡地质剖面图Figure 1. Geological section of the left-bank dam foundation slope at the Baihetan hydropower station1.2 监测点布置
随着左岸坝基边坡爆破开挖施工向河谷底部推进，LS3301、LS3302、LS331等顺坡向层内错动带和柱状节理岩体在620～580 m高程逐渐出露，这些软弱结构面交错切割岩体形成的潜在可滑动块体与临空面构成了不利组合，严重威胁边坡的稳定性。为此，现场垂直于建基面布置长短不等的预应力锚索对坝基边坡潜在不稳定块体进行加固，锚索预应力等级为3000 kN。
为掌握开挖施工过程中坝基边坡的位移变化情况及预应力锚索的工作性态，现场开展了边坡变形与预应力锚索轴力的联合监测，如图2所示。在615 m高程处布置了一个表面位移计，编号为D1；考虑到层内错动带LS3301在开挖过程中可能会产生过大的滑动而引起边坡失稳，在1号帷幕洞与1号排水洞的侧壁上分别布置了一套裂缝计监测裂缝的开合度和错动位移，编号分别为F1和F2；在607和600 m高程处3000 kN级的预应力锚索上分别安装了一套锚索测力计，编号分别为C1和C2。
图 2 边坡位移与预应力锚索轴力监测点布置Figure 2. Arrangement of the observation points for monitoringthe slope displacement and anchor cable axial force1.3 边坡位移及锚索轴力监测数据分析
2016年6月1日至2016年8月23日，白鹤滩水电站左岸坝基边坡从605 m高程开挖至580 m高程，607和600 m高程处3000 kN级预应力锚索分别在2016年6月1日和7月3日前安装完成。该段时间内边坡表观测点D1处的水平位移、内部裂缝监测点F1处的错动位移和F2处的张开位移随时间变化如图3所示。从图3可以看到，爆破开挖施工期间，各监测点的位移增长较快，特别是在595～590 m高程和590～585 m高程梯段爆破时，由于爆区距F1和F2监测点较近，F1监测点的错动位移和F2监测点的张开位移均有在爆破后大幅突然增加的现象，即位移突变。与此同时，边坡表面监测点D1也发生了多次较明显的位移突变。尽管位移监测点所在高程的预应力锚索已经安装完成，但位移突变现象仍有发生。2016年8月23日梯段爆破施工停止后，各监测点的位移增长较为缓慢，且逐渐趋于稳定。截至2016年10月29日，边坡表面测点D1处的水平位移为29.9 mm，内部监测点F1处的错动位移为7.0 mm、F2处的张开位移为6.8 mm。
图 3 各监测点位移随时间变化曲线Figure 3. Slope displacement histories at the different observation points锚索监测点C1和C2处的轴力变化如图4所示，在595～590 m高程梯段爆破开挖期间，607 m高程处C1监测点因距爆区较近且爆破规模较大，该监测点的锚索轴力出现了明显的突增现象，增量为54.6 kN，在随后的2个台阶爆破开挖过程中，C1监测点的锚索轴力也发生了不同程度的突增。位于600 m高程的C2监测点在590～580 m高程岩体爆破开挖过程中发生了较大幅度的锚索轴力突增，最大增量为53 kN。待开挖至580 m高程，爆破施工停止后，预应力锚索轴力增大变缓，轴力逐渐趋于稳定。C1和C2监测点锚索的轴力最终分别稳定在3383.9和3321.2 kN，较初始值分别增大16.2%和8.9%。
图 4 各监测点预应力锚索轴力随时间变化曲线Figure 4. Axial force histories of the prestressed anchor cables at the different observation points综上，节理岩质边坡位移的变化与锚索轴力的变化具有一致性，二者在附近有大规模岩体爆破时均出现了不同程度的突变，预应力锚索通过自身轴力的增大抑制边坡出现过大的位移突变。毫无疑问，岩体中结构面交错分布形成的具有后缘拉裂面、底滑面和侧滑面的潜在滑动块体是边坡发生位移突变的内在条件，而爆破开挖扰动是诱发节理岩质边坡位移突变的外在条件。众所周知，深切河谷具有应力降低、升高和原始三带的“驼峰形”及谷底存在“高应力包”的地应力分布基本规律，“高应力包”深度可达谷底以下150～200 m[14]。刘国锋等[15]的研究表明，白鹤滩水电站河床谷底应力在20 MPa以上。对河床谷底高地应力岩体的爆破开挖，伴随着炸药爆轰、岩体破碎及新开挖面的形成，开挖边界上的岩体地应力在岩体破碎瞬间也随之突然释放；该过程为一瞬态卸荷力学过程，在边坡保留岩体中激发瞬态卸载应力波，引起近邻开挖面的岩体应力动态调整[16-18]。因此，对于谷底高地应力岩体爆破开挖，诱发节理岩质边坡位移突变的动力扰动除爆炸荷载外，还包括地应力瞬态卸荷过程，二者几乎同步发生。然而，由于监测设备的局限性，上述实测数据仅能反映每日的边坡位移及锚索轴力变化值，无法反映爆破开挖动力扰动下边坡位移与锚索轴力的瞬态调整过程。下文将采用数值模拟方法，研究爆炸荷载与地应力瞬态卸荷联合作用下边坡位移与锚索轴力的动态变化过程及相关能量机理。
2. 数值模拟
2.1 数值模型与材料参数
采用有限差分软件FLAC3D，研究爆破开挖动力扰动下锚固节理岩质边坡的位移突变特性及能量机制。白鹤滩左岸坝基边坡地质条件复杂，岩体结构面众多，影响边坡位移突变的因素多而复杂，若将众多因素均考虑在内，很难揭示爆破开挖扰动下锚固节理岩质边坡位移突变的力学本质。为此，采用概化的节理岩质边坡模型进行数值模拟研究，以揭示爆破开挖扰动下锚固节理岩质边坡位移突变的能量机理。为便于进行机理分析，考虑爆炸应力波和瞬态卸载应力波垂直入射节理面的情况，建立含竖直和水平节理的概化边坡模型，如图5所示。模型整体尺寸为40 m×6 m×30 m，包括3个高5 m的直立台阶，上2层台阶已开挖完成，第3层台阶正在进行爆破开挖施工。第2层和第3层台阶岩体中存在一个由一组后缘竖直节理面、一组底部水平节理面和两组侧部竖直节理面切割而成的潜在滑动体，其尺寸为8 m×2 m×10 m。在第2层台阶上施加2根垂直于边坡坡面的预应力锚索，锚索总长设置为20 m，自由段与锚固段的长度分别为12和8 m，2根锚索到第3层台阶顶部的距离分别为2和4 m。卢文波等[9]、罗忆等[10]和金李等[12]也采用类似的概化模型进行了相关的研究，并取得了有益的结论；此外，卢文波等[9]的研究还表明，如底滑面为顺层的倾斜滑面时，节理岩体的位移突变特征与底滑面为水平时的情况类似，只是最终的位移值会更大一些。篇幅所限，本文中仅考虑底滑面为水平节理面的情况。
白鹤滩左岸坝肩及坝基边坡主要出露玄武岩中部P2 β3及P2 β4层，总厚度约350 m，岩性主要有斜斑玄武岩、隐晶质玄武岩、柱状节理玄武岩、杏仁状玄武岩、角砾熔岩、凝灰岩6种。高程660～555 m坝基部位出露
P2β33 层第1类柱状节理玄武岩，该层厚度约为55 m。本文中研究610～580 m高程坝基边坡岩体开挖，该高程范围内坝基边坡岩体只含有第1类柱状节理玄武岩这一种岩层，根据现场地质勘察和室内试验结果[19]，第1类柱状节理玄武岩岩石及节理物理力学参数如表1所示。数值模拟中将岩石假定为线弹性材料；对于岩体中的节理面，采用FLAC3D提供的基于库伦剪切本构模型的无厚度interface单元进行模拟，interface单元可以模拟荷载作用下2个接触面之间的错动滑移、张开、闭合等变形[20]，模拟中不考虑节理面的抗拉强度和黏聚力；预应力锚索则采用cable单元进行模拟[20]，锚索的物理力学参数如表2所示。根据Blair的研究[21]，为满足动力分析精度，数值模型的网格尺寸不应超过应力波波长的1/12～1/6，本模型的最小波长(Rayleigh波的波长)约为2.2 m，因此动力加载边界附近的单元尺寸划分为0.25 m，模型总共划分为1 793 136个单元。表 1 岩石与节理的物理力学参数Table 1. Physical and mechanical parameters of the rock and joints岩石 节理 弹性模量/GPa 泊松比 密度/(kg·m−3) 内摩擦角/(°) 法向刚度/GPa 剪切刚度/GPa 40 0.25 2700 25 20 10 表 2 锚索的物理力学参数Table 2. Physical and mechanical parameters of the anchor cable锚索索体 锚固剂 弹性模量/GPa 密度/(kg·m−3) 横截面积/mm2 黏聚力/(N·m−1) 内摩擦角/(°) 剪切刚度/GPa 锚固外圈周长/mm 锚固剂厚度/mm 200 7 890 3 465.8 2×107 38 10 518.1 49.2 2.2 动力荷载
由于预裂孔距边坡开挖轮廓最近，且预裂孔爆破后形成的预裂缝阻碍了主爆孔爆破产生的应力波向边坡岩体内传播，预裂孔爆破对边坡保留岩体的动力扰动往往是最强烈的[16]。因此，本文中仅研究预裂孔爆破动力扰动作用下锚固节理岩质边坡的位移突变。预裂爆破时每段4个炮孔同时起爆，由于炮孔直径远小于数值模型的尺寸，若将所有的炮孔均建立在模型中，为保证计算精度，数值模型的单元数将达到上千万，从而导致计算效率低甚至无法计算。为提高计算效率，本文中将炮孔壁受到的爆炸荷载压力等效施加在图5(a)所示的爆破开挖面上。根据圣维南原理，这种等效处理方法在计算炮孔附近的岩体动力响应时会产生一定的偏差，而在计算炮孔中远区的岩体动力响应时则与爆炸荷载压力施加在炮孔壁上的计算结果基本一致，该方法对于模拟中远区岩体动力响应问题的有效性已被证实[12,16,22-23]。本文中研究的竖直节理面距炮孔较远，属于中远区问题，采用该等效方法是可行的。对于预裂孔爆破，炮孔中心连线与炮孔轴线所组成的平面即为爆破开挖面。根据圣维南原理，爆破开挖面上的等效爆炸荷载压力pe(t)与炮孔壁上的爆炸荷载压力pw(t)满足如下关系：
pe(t)=(db/dbaa)pw(t) (1) 式中：db为炮孔直径，a为炮孔间距。
根据凝聚炸药爆轰波的CJ理论，炸药爆炸后作用在炮孔壁上的爆炸荷载峰值压力pw0为：
pw0=ρev2d2(γ+1)(dcdb)2υ (2) 式中：ρe为炸药密度；vd为炸药爆轰速度；γ为炸药的等熵指数，一般取3；υ为爆生气体的绝热膨胀常数，近似取1.5；dc为装药直径。
根据白鹤滩水电站左岸坝基边坡现场爆破采用的钻孔及炸药设计，取ρe=1300 kg/m3, vd=4500 m/s, db=76 mm, dc=25 mm, a=0.6 m，炮孔装药长度lc=4.5 m。由式(1)和(2)计算可得爆破开挖面上的等效爆炸荷载峰值压力pe0=15 MPa。目前获取爆炸荷载压力时程曲线的方法主要有3种：(1)半经验半理论的计算公式；(2)通过简化模型分析爆生气体的运动过程；(3)采用数值模拟软件（如LS-DYNA）模拟炸药爆轰过程。由于炸药爆炸产生的高温高压环境，现有的测试手段还很难从现场准确获取炮孔壁上或爆破开挖面上的爆炸荷载压力时程曲线。在实际工程爆破的数值模拟中，目前大多采用半经验半理论的爆炸荷载压力时程曲线，其中以双指数函数和进一步简化的三角形函数压力时程曲线应用最为广泛。本文中采用三角形荷载曲线，如图6所示，取爆炸荷载压力上升时间tri=1 ms、持续时间td=9 ms，其计算过程详见文献[16]和[23]，爆炸荷载压力时程曲线pe(t)的表达式为：
图 6 爆破开挖面上等效爆炸荷载与地应力瞬态卸荷历程曲线Figure 6. Time-histories of the equivalent blasting pressure and the transient unloading of in-situ stress on the blasting excavation boundarype(t)={pe0ttri0≤t＜tripe0td−ttd−tritri≤t≤td0t＞td (3) 同一排炮孔同时起爆时，爆生裂纹优先在炮孔连线方向上扩展，当裂纹完全贯通、开挖面上的爆炸荷载压力衰减至与开挖面上的地应力大小相等时，开挖面岩体才出现宏观上的卸荷效应[9]。根据应力连续条件，爆破开挖面上的地应力瞬态卸荷历程在岩体出现宏观卸荷后与爆炸荷载作用曲线重合，如图6所示，地应力瞬态卸荷时程曲线σ(t)可表示为：
σ(t)={σ00≤t＜tbσ0td−ttd−tbtb≤t≤td0t＞td (4) 式中：σ0为爆破开挖面上的地应力；tb为卸荷开始时间。
第3层最后一个梯段预裂爆破时，取爆破开挖面上的地应力σ0=5 MPa，根据图6，地应力瞬态卸荷持续时间tu=td−tb=2.7 ms。需要说明的是，虽然白鹤滩水电站河床谷底的初始地应力较高（20 MPa以上）[15]，但随着开挖工作面逐步向坝基边坡轮廓面推进，边坡轮廓面上的初始地应力也逐步释放，因此在最后一个梯段爆破时，坝基边坡轮廓面上的地应力取一个较小的值。
数值模拟过程中，首先进行地应力和预应力锚索作用下的静力求解，模型四周设置为法向位移约束，待计算达到平衡后打开动力求解模块，将上述爆炸荷载和地应力瞬态卸荷作用历程用FISH语言编写后施加在爆破开挖面上进行动力求解。为防止动力荷载激发的应力波传播至模型四周边界发生反射而影响计算结果，动力计算中将模型四周设置为可吸收入射波的黏性边界。岩体阻尼采用FLAC3D提供的局部阻尼，局部阻尼系数设置为0.157。
3. 节理岩质边坡位移突变特征与能量机理
为对比分析预应力锚索对节理岩质边坡位移突变的影响，本文数值模拟分2类工况进行研究，一类未施加预应力锚索，另一类为施加预应力锚索工况。由于岩体出现朝向河谷的位移时对边坡稳定最为不利，因此本文在潜在滑动块体的前缘和后缘各选取一个监测点A和B（见图5）分析边坡的水平位移及应变能密度特性。
3.1 未施加预应力锚索工况
未施加预应力锚索时，爆破开挖动力扰动下两监测点的水平位移和应变能密度时程曲线如图7所示。在初始地应力σ0的作用下，节理岩体本身储存有一定量的应变能；爆炸荷载作用后，节理岩体被压缩，岩体中进一步积聚应变能。随着爆炸荷载压力降低以及地应力释放，岩体开始回弹，岩体中的应变能得以释放。在t=9 ms左右时，爆炸荷载压力降为零，被压缩的岩体回弹恢复到爆炸荷载作用前的状态，此时监测点的位移为零。但由于爆炸荷载压力快速降低和地应力瞬态释放所产生的惯性力作用，潜在滑动块体继续向河谷方向运动，导致后缘节理张开，产生了节理张开位移。由于底部及两侧滑面的摩擦作用，在t=30 ms左右时，岩块停止了运动，A和B两监测点的最终位移分别为6.5和5.6 mm，该位移即为爆破开挖动力扰动过程中节理岩质边坡的突变位移。爆破开挖完成后，A和B两监测点的最终水平位移之差为0.9 mm；而爆破前，在初始地应力σ0的作用下，A、B两点的水平位移也正好相差0.9 mm（即该岩块的压缩变形量为0.9 mm），这表明A和B两点最终位移之差是由于地应力释放、岩体松弛所产生的回弹位移。对于河谷底部高地应力边坡岩体爆破开挖，边坡表面岩体的突变位移包括节理张开位移和回弹位移两部分。对于本算例中边坡表面的A监测点，总突变位移为6.5 mm，节理张开位移为5.6 mm，回弹位移为0.9 mm，节理张开位移占到了总突变位移的86.2%，是节理岩质边坡位移突变的主要组成部分。
图 7 无预应力锚索时各监测点的岩体位移和应变能密度时程曲线Figure 7. Rock mass displacement and strain energy density histories at the observation points without prestressed anchor cables从能量的观点来看，岩体由于地应力作用本身所存储的应变能以及爆炸荷载挤压作用所集聚的应变能为节理边坡岩体位移突变提供了能量来源，上述2种应变能的快速释放是节理岩体位移突的直接诱因。代金豪等[3]的相关研究表明，岩体存储和集聚的应变能越高、应变能释放持续时间越短，产生的位移突变就越大，即位移突变与应变能释放率正相关。应变能释放率
rSER 为：rSER=ESER, maxtSER (5) 式中：
ESER,max 为岩体积聚的体积应变能峰值，tSER 为应变能释放持续时间。因此，在保证开挖岩体正常破碎的前提下，通过优化爆破方案，降低岩体应变能的释放率有助于减小节理岩质边坡的位移突变。如采用低密度、低爆轰波速的炸药降低爆炸荷载峰值压力，从而降低岩体集聚的应变能；采用小抵抗线降低开挖面上的地应力，从而降低岩体存储的应变能；加强炮孔堵塞减缓爆生气体从孔口逸出，从而延长爆炸荷载压力衰减时间，即延长应变能释放持续时间。
3.2 施加预应力等级锚索工况
除爆破方案优化外，工程中还主要采用锚杆、锚索等支护方式来控制节理岩质边坡的位移突变。为探讨预应力锚索对节理岩质边坡位移突变的控制机理，设置锚索预应力等级分别为1 000、1 500、 2 000、2 500和3 000 kN。图8给出了不同等级预应力锚索作用下边坡表面A点的位移和应变能密度变化时程曲线。可以看到，不同等级预应力锚索作用下，节理岩体的位移总体特征基本相同，均具有压缩→回弹→突变→稳定这4个阶段。不同预应力等级下节理岩体的压缩与回弹路径完全重合，这表明预应力锚索并不影响节理张开前岩体的位移与能量转化，不同等级预应力锚索作用下A点的应变能密度时程曲线完全重合也证明了这一点。在t=9 ms节理张开发生位移突变后，锚索预应力等级越高，节理岩体的突变位移越小。当预应力超过1 000 kN时，节理岩体在达到最大突变位移后会出现明显的“被拉回”现象，最终使突变位移稳定在较小的值。当预应力等级超过1 500 kN时，随着锚索预应力等级的提高，节理岩体的最终突变位移虽有所减小，但减小幅度不大。图9给出了不同等级预应力锚索作用下A监测点的节理最终张开位移和岩体最终回弹位移，随着预应力等级的提高，节理张开位移同样呈现先显著后缓慢降低的趋势，而岩体回弹位移基本保持不变，这表明预应力锚索主要是通过抑制节理张开来降低边坡岩体的位移突变。结合前述实测资料分析结果，预应力锚索的施加虽不能完全避免节理岩质边坡在爆破开挖过程中产生位移突变，但采用适当预应力等级的锚索进行加固后可有效减少岩体内节理的张开位移，从而防止边坡岩体位移突变过大。
图 8 不同等级预应力锚索作用下监测点A的位移和应变能密度的时程曲线Figure 8. Rock mass displacement and strain energy density histories at the monitoring point A under different levels of the prestressed anchor cables
图 9 不同预应力等级下的节理最终张开位移与岩体最终回弹位移Figure 9. Joint opening displacement and rebound displacement under different levels of the prestressed anchor cables针对预应力锚索在控制节理岩质边坡位移突变所发挥的作用，从能量的角度对其作用机理开展讨论。对于预应力锚索，其总的应变能W由锚索索体的拉压弹性应变能Wca和锚固剂的剪切弹性应变能Wg组成[20]，分别由下式计算：
W=Wca+Wg (6) Wca=M∑i=1F2iLi2EcaS (7) Wg=N∑j=1τ2jVj2G (8) 式中：Fi为第i个索体单元的轴力；Li为第i个索体单元的长度；Eca为锚索索体的弹性模量；S为锚索索体横截面面积；τj为第j个锚固剂单元的剪切应力；Vj第j个锚固剂单元的体积；G为锚固剂剪切模量；M和N分别为索体单元和锚固剂单元的数量。
通过FLAC3D中的FISH语言对锚索总应变能W涉及的变量进行实时监测和计算，得到不同等级预应力锚索总应变能的时程变化曲线，如图10所示（篇幅所限，仅以1000和3000 kN等级为例）。锚索总应变能在开始突增前经历了降低和升高至初始值这2个阶段，这与节理岩体在位移突变前发生压缩和回弹相对应；在锚索应变能突增开始后存在吸能和释能2个阶段，定义锚索的平均吸能速率或释能速率K为：
图 10 不同等级预应力锚索的总应变能时程变化曲线Figure 10. Time-histories of the total strain energy of the anchor cables under different prestress levelsK=1nn∑k=1|ΔWkΔtk| (9) 式中：Δtk为第k个时间段(k=1, 2, 3,···, n)；ΔWk为第k个时间段开始与结束时锚索的总应变能之差。锚索平均吸能和释能速率分别表示为Ka和Kre。
结合图8(a)和图10可以看到，节理岩体发生位移突变后，节理岩体与锚索相互作用，节理岩体的动能转化为锚索的应变能，锚索应变能急剧增加，表现为吸能，此时锚索被拉伸且轴力增大，这与实测资料反映的规律相对应，即节理岩体位移突变与锚索轴力突增具有同步性；锚索的预应力等级越高，吸能速率越高，节理岩体以更高的速率达到较小的突变位移峰值；在节理岩体的位移达到最大值后，锚索的应变能开始释放，释放的应变能作用于节理岩体，部分转化为节理岩体的动能，将滑动的岩块逐渐“拉回”，随着底部及两侧滑面的摩擦耗能，岩块最终停止运动，位移趋于稳定；锚索的预应力等级越高，释能速率也越高，节理岩体“被拉回”的效应也就越明显，即节理岩体以更高的速率稳定至较小的最终位移值。综上，预应力锚索通过自身具有的吸能与释能特性来控制边坡节理岩体的位移突变。
受构造应力、边坡结构、开挖方案等多方面影响，爆破开挖面上的初始地应力大小多有不同，此外，不同工程的岩体变形参数也各不相同。不同地应力水平和岩体弹性模量条件下监测点A的最终突变位移如图11所示。从图11(a)可以看到，随着地应力水平的提高，节理岩质边坡的突变位移增大，这是由于较高的地应力会使岩体内储存更多的应变能，相同释放时间情况下具有更高的应变能释放率；当爆破开挖面上的初始地应力由2.5 MPa升高12.5 MPa时，与未施加预应力锚索相比，施加3 000 kN等级的预应力锚索后边坡突变位移分别减小了64.2%和58.7%，控制效果显著。从图11(b)可以看出，随着岩体弹性模量的提高，节理岩质边坡的突变位移逐渐减小，这是由于在相同的应力作用下，较高弹性模量的岩体应变较小，岩体所存储的应变能较低，从而使应变能释放率降低。对于弹性模量为20 GPa的岩体，未施加预应力锚索和施加3 000 kN等级的预应力锚索情况下，边坡突变位移由11.2 mm减小为5.9 mm，减小了47.3%；而对于弹性模量为60 GPa的岩体，施加3 000 kN等级的预应力锚索后位移突变降低了59.2%，预应力锚索对岩体弹性模量不同的岩质边坡均体现出了较好的控制效果。由此可见，对于地应力水平较高或岩体弹性模量较低的节理岩质边坡爆破开挖，尤其要注意节理岩体的位移突变问题，可考虑采用适当等级的预应力锚索进行加固。
图 11 岩体地应力和弹性模量对不同等级预应力锚索作用下监测点A的最终突变位移的影响Figure 11. Effects of in-situ stress level and elastic modulus of the rock mass on final mutation displacement of the monitoring point A under different levels of the prestressed anchor cables4. 结 论
（1）白鹤滩水电站左岸坝基锚固节理岩质边坡在爆破开挖过程中出现了明显的岩体位移突变和锚索轴力突增的现象，二者具有同步性。
（2）对于深切河谷底部高地应力边坡岩体爆破开挖，爆炸荷载对岩体挤压作用所集聚的应变能以及地应力作用下岩体所存储的应变能为节理岩质边坡位移突变提供了能量来源，爆破过程中这2种应变能的快速释放是节理岩体位移突变的直接诱因。
（3）深切河谷底部高地应力边坡岩体在爆破开挖过程中产生的突变位移包括节理张开位移和岩体回弹位移，爆破开挖面上初始地应力水平越高、岩体弹性模量越低，产生的突变位移越大；预应力锚索主要控制节理张开位移，锚索预应力等级越高，其吸能和释能速率越高，位移突变控制效果越明显，但当锚索的预应力等级高到一定程度后，节理岩体的突变位移不再明显减小。
为揭示爆破开挖扰动下锚固节理岩质边坡位移突变的能量机理，本文的数值模拟研究中做了一些简化和假定，岩体结构面实际形态、爆破开挖扰动下岩体及锚固结构的损伤对边坡位移突变的影响还需要开展进一步的研究。
-
图 1 药型罩顶部微元和整体受力分析
Figure 1. The top micro-element and overall force analysis of the liner
图 2 聚能侵彻体处于不同阶段示意图
Figure 2. Schematic diagram of an explosively-formed projectile in different stages
图 3 聚能侵彻体/直达冲击波速度随时间的变化
Figure 3. Velocity versus time curves of explosively-formed projectile and direct shock wave
图 6 侵彻体和直达冲击波速度随时间的变化曲线
Figure 6. Velocity versus time curves of explosively-formed projectile and direct shock wave
图 7 侵彻体和直达冲击波传播距离随时间的变化曲线
Figure 7. Propagation distance versus time curves of explosively-formed projectile and direct shock wave
图 8 侵彻体和直达冲击波的传播
Figure 8. Propagation of explosively-shaped projectile and direct shock wave
A1/GPa B1/GPa R1 R2 ω1 ρe/(kg·m−3) DCJ/(m·s−1) Ee/(GJ·m−3) pCJ/GPa 778.28 7.07 4.20 1.00 0.30 1891 9110 10.5 42.0 
下载: 导出CSV
表 2 不同数值模型对应的工况
Table 2. Different cases corresponding to numerical models
工况 介质 空气域长度 1 空气 无限 2 水 20倍装药半径 3 水 5倍装药半径
下载: 导出CSV
-
[1] LI J T, HUANG C, MA T, et al. Numerical investigation of composite laminate subjected to combined loadings with blast and fragments [J]. Composite Structures, 2019, 214: 335–347. DOI: 10.1016/j.compstruct.2019.02.019. [2] ZHANG L, WANG X, WANG Y T, et al. High-hardness polyurea coated steel plates subjected to combined loadings of shock wave and fragments [J]. Latin American Journal of Solids and Structures, 2022, 19(2): e433. DOI: 10.1590/1679-78256882. [3] WANG C L, ZHOU G, MA K, et al. Research on damage of typical water-containing composite structure under different types of charge [J]. Journal of Physics: Conference Series, 2020, 1507(3): 032066. DOI: 10.1088/1742-6596/1507/3/032066. [4] 侯海量, 张成亮, 李茂, 等. 冲击波和高速破片联合作用下夹芯复合舱壁结构的毁伤特性 [J]. 爆炸与冲击, 2015, 35(1): 116–123. DOI: 10.11883/1001-1455(2015)01-0116-08.HOU H L, ZHANG C L, LI M, et al. Damage characteristics of sandwich bulkhead under the impact of shock and high-velocity fragments [J]. Explosion and Shock Waves, 2015, 35(1): 116–123. DOI: 10.11883/1001-1455(2015)01-0116-08. [5] 张成亮, 朱锡, 侯海量, 等. 爆炸冲击波与高速破片对夹层结构的联合毁伤效应试验研究 [J]. 振动与冲击, 2014, 35(15): 184–188. DOI: 10.13465/j.cnki.jvs.2014.15.032.ZHANG C L, ZHU X, HOU H L, et al. Tests for combined damage effect of blast waves and high-velocity fragments on composite sandwich plates [J]. Journal of Vibration and Shock, 2014, 35(15): 184–188. DOI: 10.13465/j.cnki.jvs.2014.15.032. [6] 李茂, 高圣智, 侯海量, 等. 空爆冲击波与破片群联合作用下聚脲涂覆陶瓷复合装甲结构毁伤特性 [J]. 爆炸与冲击, 2020, 40(11): 111403. DOI: 10.11883/bzycj-2019-0119.LI M, GAO S Z, HOU H L, et al. Damage characteristics of polyurea coated ceramic/steel composite armor structures subjected to combined loadings of blast and high-velocity fragments [J]. Explosion and Shock Waves, 2020, 40(11): 111403. DOI: 10.11883/bzycj-2019-0119. [7] NYSTRÖM U, GYLLTOFT K. Numerical studies of the combined effects of blast and fragment loading [J]. International Journal of Impact Engineering, 2009, 36(8): 995–1005. DOI: 10.1016/j.ijimpeng.2009.02.008. [8] 吴华鑫, 丁剑锋, 程瑞琪, 等. 舰船防护结构抗近爆毁伤作用数值模拟 [J]. 兵器装备工程学报, 2021, 42(4): 91–96. DOI: 10.11809/bqzbgcxb2021.04.017.WU H X, DING J F, CHENG R Q, et al. Numerical simulation of near explosion damage effect on ship protective structure [J]. Journal of Ordnance Equipment Engineering, 2021, 42(4): 91–96. DOI: 10.11809/bqzbgcxb2021.04.017. [9] 夏冰寒, 王金相, 周楠, 等. 柱状装药预制破片缩比战斗部爆炸冲击波和破片的作用时序 [J]. 高压物理学报, 2020, 34(1): 015101. DOI: 10.11858/gywlxb.20190780.XIA B H, WANG J X, ZHOU N, et al. Blast wave and time sequence of prefabricated fragments for scaled warhead with cylindrical charge [J]. Chinese Journal of High Pressure Physics, 2020, 34(1): 015101. DOI: 10.11858/gywlxb.20190780. [10] 陈长海, 侯海量, 朱锡, 等. 破片式战斗部空中爆炸下冲击波与破片的耦合作用 [J]. 高压物理学报, 2018, 32(1): 015104. DOI: 10.11858/gywlxb.20170849.CHEN C H, HOU H L, ZHU X, et al. Coupling action spans for air-blast waves and fragments by fragmentation warheads exploding in air [J]. Chinese Journal of High Pressure Physics, 2018, 32(1): 015104. DOI: 10.11858/gywlxb.20170849. [11] 郑红伟, 陈长海, 侯海量, 等. 爆炸冲击波和高速破片载荷的复合作用特性及判据研究 [J]. 振动与冲击, 2019, 38(3): 24–31, 38. DOI: 10.13465/j.cnki.jvs.2019.03.004.ZHENG H W, CHEN C H, HOU H L, et al. Multiple impact features of blast shock waves and high-velocity fragments on clamped square plates and a criterion to judge if multiple impact happens [J]. Journal of Vibration and Shock, 2019, 38(3): 24–31, 38. DOI: 10.13465/j.cnki.jvs.2019.03.004. [12] 吕勇, 石全, 钱芳, 等. 不同时序破片和冲击波在对天线的复合毁伤分析 [J]. 火力与指挥控制, 2014, 39(3): 128–132,139. DOI: 10.3969/j.issn.1002-0640.2014.03.032.LÜ Y, SHI Q, QIAN F, et al. Analysis of multiplex damages for fragment and shock wave in different time sequence to antenna [J]. Fire Control & Command Control, 2014, 39(3): 128–132,139. DOI: 10.3969/j.issn.1002-0640.2014.03.032. [13] 赵德辉, 许金余, 张燕. 武器爆炸条件下冲击波与弹片共同作用数值模拟分析 [J]. 西北工业大学学报, 2006, 24(3): 334–337. DOI: 10.3969/j.issn.1000-2758.2006.03.014.ZHAO D H, XU J Y, ZHANG Y. Velocities and displacements of shrapnel and shock wave during blast [J]. Journal of Northwestern Polytechnical University, 2006, 24(3): 334–337. DOI: 10.3969/j.issn.1000-2758.2006.03.014. [14] 侯俊亮, 蒋建伟, 李应波, 等. 杀爆战斗部破片和冲击波对目标的耦合作用 [J]. 火炸药学报, 2020, 43(3): 335–340. DOI: 10.14077/j.issn.1007-7812.201909009.HOU J L, JIANG J W, LI Y B, et al. Study on the coupling effect of fragments and shock wave of the blast-fragmentation warhead on typical targets [J]. Chinese Journal of Explosives and Propellants, 2020, 43(3): 335–340. DOI: 10.14077/j.issn.1007-7812.201909009. [15] 李峰, 石全, 张芳, 等. 破片和冲击波复合作用下装甲板毁伤效应预测 [J]. 火力与指挥控制, 2020, 45(11): 101–105. DOI: 10.3969/j.issn.1002-0640.2020.11.019.LI F, SHI Q, ZHANG F, et al. Prediction of damage effect of armor plate under combined action of fragments and shock wave [J]. Fire Control and Command Control, 2020, 45(11): 101–105. DOI: 10.3969/j.issn.1002-0640.2020.11.019. [16] 王远飞. 接触爆炸作用金属材料中冲击波的衰减 [D]. 南京: 南京理工大学, 2007: 5–6. [17] 王作山, 刘玉存, 郑敏, 等. 爆轰冲击波在有机隔板中衰减模型的研究 [J]. 应用基础与工程科学学报, 2001, 9(4): 316–319. DOI: 10.3969/j.issn.1005-0930.2001.04.004.WANG Z S, LIU Y C, ZHENG M, et al. Study on the attenuating model of detonation shock wave in the PMMA gap [J]. Journal of Basic Science and Engineering, 2001, 9(4): 316–319. DOI: 10.3969/j.issn.1005-0930.2001.04.004. [18] 李金河, 尚海林, 傅华. 平面冲击波在有机玻璃中的衰减规律 [J]. 高压物理学报, 2018, 32(4): 044201. DOI: 10.11858/gywlxb.20170643.LI J H, SHANG H L, FU H. Attenuation of plane shock wave in PMMA [J]. Chinese Journal of High Pressure Physics, 2018, 32(4): 044201. DOI: 10.11858/gywlxb.20170643. [19] 刘飞. 爆炸成型弹丸(EFP)研制及其工程破坏效应研究 [D]. 安徽合肥: 中国科学技术大学, 2006: 26–30.LIU F. The explosively formed penetrator and its engineering damage effects research [D]. Hefei, Anhui, China: University of Science and Technology of China, 2006: 26–30. [20] 杨莉, 张庆明, 巨圆圆. 爆炸成型弹丸对含水复合装甲侵彻的实验研究 [J]. 北京理工大学学报, 2009, 29(3): 197–200. DOI: 10.15918/j.tbit1001-0645.2009.03.007.YANG L, ZHANG Q M, JU Y Y. Experimental study on the penetration of explosively formed projectile against water-partitioned armor [J]. Transactions of Beijing Institute of Technology, 2009, 29(3): 197–200. DOI: 10.15918/j.tbit1001-0645.2009.03.007. [21] 黄正祥. 聚能装药理论与实践 [M]. 北京: 北京理工大学出版社, 2014: 19–59.HUANG Z X. Theory and practice of shaped charge [M]. Beijing, China: Beijing Institute of Technology Press, 2014: 19–59. [22] Century Dynamics Incorporated. Interactive non-linear dynamic analysis software AUTODYNTM user manual [M]. Concord, USA: Century Dynamics Incorporated, 2003. [23] 张晓伟, 段卓平, 张庆明. 钛合金药型罩聚能装药射流成型与侵彻实验研究 [J]. 北京理工大学学报, 2014, 34(12): 1229–1233. DOI: 10.15918/j.tbit1001-0645.2014.12.004.ZHANG X W, DUAN Z P, ZHANG Q M. Experimental study on the jet formation and penetration of conical shaped charges with titanium alloy liner [J]. Transactions of Beijing Institute of Technology, 2014, 34(12): 1229–1233. DOI: 10.15918/j.tbit1001-0645.2014.12.004. 期刊类型引用(12)
1. 孙鹏昌，杨广栋，卢文波，范勇，孟海利，薛里. 考虑岩体破坏分区的岩石爆破爆炸荷载历程研究. 爆炸与冲击. 2024(03): 171-186 . 
本站查看2. 孙鹏昌，杨招伟，卢文波，孟海利，薛里. 白鹤滩水电站坝肩边坡开挖爆破振动概率分布研究. 工程科学与技术. 2024(03): 179-188 . 
百度学术3. 李聪然，马明星，潘懿，宋军涛，杜晓飞，刘文进，岳祖润. 季冻区露天矿山边坡稳定性分析及加固措施研究. 河北地质大学学报. 2024(03): 73-79 . 百度学术4. 和铁柱，邱海燕，王基禹，孟凡豹，刘康琦. 节理倾角对岩质边坡稳定性影响的数值模拟研究. 工程建设与设计. 2024(11): 33-36 . 百度学术5. 段宝福，刘凯强，柴明星，孙宗军，侯亚洪，陈佳华. 隧道爆破对燃气管道的振动响应分析. 工程爆破. 2024(04): 142-149 . 百度学术6. 冷振东，贾永胜，卢文波，孙金山，谢先启. 高孔隙软岩爆破能量传输规律与爆破大块形成机制. 应用基础与工程科学学报. 2024(05): 1420-1433 . 百度学术7. 邰成群，杨春喜，张祝九. 三孔爆破作用下岩石损伤特点研究. 山西冶金. 2023(08): 67-70 . 百度学术8. 朱颂波，徐振洋，张祚富，李政，郇宝乾. 考虑炮孔密集系数及高程影响下质点峰值速度的预测模型研究. 有色金属工程. 2023(09): 139-149 . 百度学术9. 刘洋，杨天鸿，李华，杨意德，赵永，李金多，邓文学. 研山铁矿西帮反倾岩质边坡变形监测及数值模拟分析. 金属矿山. 2023(11): 261-267 . 百度学术10. 何少云，赵修龙，谭支超，张帅. 环境敏感区边坡无爆破切割开挖及支护稳定性评价. 水利科学与寒区工程. 2022(06): 129-132 . 百度学术11. 段卫娜. 水利工程中高边坡开挖支护技术研究. 中国新技术新产品. 2022(13): 118-120 . 百度学术12. 张霄，于昊，李壮，刘衍顺，张紫东，籍鑫雨，李相辉. 循环加载强化作用对花岗岩细观破坏影响的离散元研究. 浙江大学学报(工学版). 2022(11): 2303-2312 . 百度学术其他类型引用(5)
-
下载:
百度学术
计量
- 文章访问数: 547
- HTML全文浏览量: 139
- PDF下载量: 309
- 被引次数: 17