第二章 岩体内的应力波

1、第二章 应力波理论与岩体内的爆炸应力 波 2.1 热力学的基本概念 2.2 液体和气体中的应力波（自习） 2.3 应力波 2.4 岩体内的爆炸应力波 2.5 弹性波的波动方程（自习） 2.6 平面波、球面波、柱面波的波动方程（自习） 2.7 纵波在两种岩石界面的折射和反射 2.8 应力波的干涉 2.1 热力学的基本概念 2.1.1基本概念 理想液体和气体 理想液体或气体是其粒子之间既无剪力又无摩擦力作用的 一种连续介质。 热力学基本参量 流体的状态参数常用压力P，比容V或密度,温度T,比内 能e，比熵s等热力学状态参量来描述。其中压力、密度和温度 是可以直接测量的量，称为热力学基本参量。 热力

2、学零定律 流体热力学体系化学成分不因混合或扩散而改变时，流体 的状态可由任何两个独立的热力学参数来确定，该状态原理有 时又称热力学第零定律。 状态方程 热力学参量P、V、T、e、s中任何两个都可以选作为独立变 量，其他量都可用它们表示，这些热力学量之间的关系称为状 态方程。 状态方程决定三个参数之间的关系，通常具有以下形式： , , (22) 根据气体和液体的运动学理论： (23) 其中： ， 都是函数 。 0),(TPF),(SVPP ),(TPP )( )( V V TVf P )(Vf)(V 熵 熵指熵指的是混乱的程度热力学中表征物质状态的参量之一，通常 用符号S表示。在经典热力学中，

3、熵的定义： 或 (21) 式中：比热； 绝对温度； 介质的两种状态。 T dQ ds B A AB T dQ SS dQ TBA, 体系的热力学关系式 物质的基本粒子（原子、分子）的运动能量和相互 的效应称为物质的内能。单位质量的内能称为比内能， 用E表示，并由状态的参数决定： ， ， (24) 热力学第一定律表示为能量守恒： (25) ),(TPEE ),(VEE ),(TVE dAdEdQ 根据方程（25）可知，体系内的介质所吸收的热量 dQ分别消耗于完成内能与功的变化量dE和dA。量Q、E 和A都是就单位质量而言，功常用dA=PdV表示，于是 有： (26) 将式： (27) 代入（26

4、）式得： (28) PdVdEdQ dV V E dT T E dE TV PdVdA dVP V E T dT T E dQ V 体积恒定时，则（28）式为： (29) 压力恒定时的比热为： C CV V是在体积恒定情况下，当单位质量的物质温度升高一 个单位时的内能增量; C CP P是在压力恒定情况下，当单位质量的物质温度升高一 度时的热量增量，在此等压过程中，体积是变化的。对 于所有物质，CPCV ，这些定义说明，差（CP-CV）等于 单位质量的物质温度升高一度时膨胀所作的外功。 dT T E dQ V VV V T E dV dQ C P P dT dQ C (210) 热量的总变化可

5、分为内部的变化（粒子之间的摩 擦）和外部的变化（加热和放热）即： 所以 对于绝热过程（ ）： 因为粒子之间的摩擦的功恒为正（ ），因 此，对于绝热过程始终有 。 对于理想物质： 绝热过程称为理想绝热或等熵绝热过程。 ACC vp ext dQdQdQ int TdQdQTdQdS ext / )(/ int 0 ext dQ TdQdS/ int 0 int dQ 0dS 0dS 绝热系统中熵无变化的过程称为可逆过程。在不 可逆过程中，熵是增大的。使介质的熵增加的不可逆 过程发生在冲击波的波阵面上；波阵面上的压力和其 它量的变化为阶跃式的，压力跃变值越大，冲击波波 阵面上的热损耗和介质的熵的增

6、加越多。 在等熵运动情形中，S=S0=常数，方程（22）式 有形式： 即压力仅与体积或密度有关，等熵并不意味着温 度不变。 )()(),( 0 PVPSVPP v2.1.2理想气体与正压介质 1、理想介质 理想气体是指在其粒子之间既无摩擦力又无结合 力作用的气体，因此，它的内势能为零，而整个内能 等于运动粒子的动能，所以内能与绝对温度有关。 若在特殊情况下，内能与温度成正比，这种理想 气体称为多方气体。 或 (212) 热力学第一定律（25）式为： (213) )(TEE dTCdE V TCE V dAdEdQ 在可逆绝热过程中，不与外部介质发生热交换， 于是对多方气体，我们有 ， 或 (2

7、14) 克拉伯龙门捷列夫理想气体的状态方程为： (215) 式中0R是气体常数，当P=常数时，对式（215）进行 微分，得到 ，热力学第一定律取以下形式： (216) 式（216）两边同除以 得迈耶公式： (217) 0dQ 0 PdVdTCV RTPV 0 RdTPdV 0 RdTdTCdQ v 0 dT RCC Vp 0 对于多方气体，由于式（21）、（215）、（2 16）得到： (218) (219) (220) v多方气体的熵是两个变量：温度和体积或压力和体 积的函数。解（220）得： v （221） VRTC V dV R T dT CdS VV lnln 00 0V RC0 V0

8、 VTlnRlnVlnTCS-S P V 0 V C C 0 R C 0 0 V R P lnV R PV lnS-S k* (S)AP V P C C k V 00* C S-S Rexp(S)A 多方气体的状态方程也可以表示为方程（2 21）形式，在等熵过程中，熵是常数。A*(S)=A*=常 数，多方气体状态方程变为： 或 （222） 这里A*是常数，k是等熵指数。 声波的传播速度为： （223） 若 与 分别表示压力为时介质的声速与密度，由 （222）、（223）得： （224） k AP* k VAP * 1 -k* S z kA P C 1-k 0 2 z0* k C A k C P

9、 0 2 z0 0 0Z C 0 v利用上述关系，理想气体的内能表示为： （225） 1-kk C PV 1-k 1 PV C-C C R PVC TCE 2 z VP V 0 V V v2、正压介质 定义：状态方程中不包含熵的介质称为正压介质。因 此，这种介质的状态方程形式是： （2-26） 在这种介质中，熵的变化对压力的影响可以忽略 不计，水、饱和土以及诸如油水气、石油水、 石油气、水气等多相介质都属于这种介质。 v平均压缩率： v （2-27） 当压力为P0时，比容为，压力为P1时，比容为 。 )()(PVPP 01 10 P-P V-V 1 V v瞬时压缩率： (228) v式中： (

10、229) vCz 称为介质的特征阻抗或波阻抗。 v对于不可压缩物体，0，A 故声速Cz v压缩率（228）表示P=P（V）曲线上一点的某性质， 若利用平压缩率，这个性质就定义在一个压力区间。 v根据热力学第一定律，并考虑到（21）式，对正 压介质的内能便得到关系： v (230) 2 z 2 C 1 A 1 P V - 2 2 z SS z V 1 -C- VV P - V P -CA PdVTdSdE v内能和压力可以表示成为取作基本量的两个参数 体积和熵的函数E=E(V,S) v (231) v比较（30）式与（31）式得： ， v根据状态方程，压力只与体积有关，因此： v (232) v

11、由此可知正压介质是内能分成两个函数的介质，其 一只与熵，另一只与体积有关。 v有些正压介质的状态比较复杂，水和其它均质液体 的状态方程一般表示为： v (233) dV V E dS S E dE T S E P V E )(),(S S E V V E )()( 21 SEVEE * BAP k v系数A*与B*是试验确定的常数。 v根据公式（223）、（224）得： k1时，Cz随密度的增加而增大，k1时，Cz与密度无 关，而Cz1时，随密度的增大而减少。 v P=P0时, 因而 v (234) v式（233）可表示为： v v (235) 1*2 k z kA d dP C 1 0 *2

12、 k z kAC 1 0 2 0* k z k C A * 1 0 0 2 0 B k C P k k z O 0 0 2 0* P k C B Z 11 0 0 2 00 0 0 2 00 0 k z k z V V kV C P k C PP 2.3应力波 v2.3.1基本概念 v在介质中传播的扰动称为波.由于任何有界或无界介质 的质点是相互联系着的,其中任何一处的质点受到外界 作用而产生变形和扰动时,就要向其他部分传播,这种 在应力状态下介质质点的运动或扰动的传播,称为应力 波. v固体爆炸应力波的定义：装药在岩体或其他固体介质 中所激起的应力扰动的传播，称为爆炸应力波，即介 质处于某一

13、应力状态下的部分和介质其余部分的界面 称为波阵面。 v按应力波传播途径： v应力波可分为体积波和表面波 v体积波：是指在固体介质中传播的应力波，体积波 又可分为纵波（P波）和横波（S波） v纵波是由波源向外传播的一种压缩波，其前进方向 与质点的振动方向一致图 2.2(a)，介质在其作用下被 压缩和膨胀，其特点是周期短、振幅小、波速快，引 起地面竖向振动，破坏性较弱； v横波是由波源向外传播的一种剪切波，其前进方向 与质点振动方向垂直图2.2(b)，其特点是周期长、振 幅相对纵波大、波速慢引起地面水平振动，破坏性较 强 。 图2-7 体积波传播示意图 v在地层表面或介质体表面传播的地震波称为面波

14、， 它通常被认为是体波经地层界面的多次反射后形成的 次生波，主要包括瑞利波(R 波)和洛夫波(L 波)。 v在瑞利波的传播过程中，介质质点在地表面内及波 的传播方向上作椭圆运动，但在与该平面垂直的水平 方向没有振动，因而瑞利波质点运动轨迹在地面上呈 滚动形态图 2.3(a) v在勒夫波的传播过程中，质点仅在与传播方向相垂 直的水平方向运动，在地面上呈蛇形运动形式图 2.3(b)。其特点是振幅大、周期长，在地面传播时衰 减慢，但其携带的能量比较大，影响范围较广，也是 造成建筑物发生强烈振动甚至破坏的主要因素。 在整个爆破过程中，体波是导致岩石产生破裂的主要 因素，而面波则是引起爆破地震效应的主要

15、因素。 图2-8 面波传播示意图 v按波阵面形状分类 v球面波、柱面波和平面波 v按传播介质变形性质不同分类 v由于固体介质变形性质不同，在固体中传播的应力 波可以分为以下几种： v弹性波：在弹性介质传播的波，服从虎克定律 v黏弹性波：在非线性弹性体中传播的波。这种波 除了弹性变形产生的弹性应力外，还产生摩擦应力或 黏滞应力 v塑性波。应力超过弹性极限的波。在能够传播塑性 波的介质中，应力在未超过弹性极限前仍然是弹性的。 应力超过弹性极限后，出现屈服应力，起传播速度比 弹性应力传播速度小得多。 v冲击波。如果介质的变形性质能使大扰动的传播速 度远远大于小扰动的传播速度，在介质中就会形成波 阵面

16、陡峭的，已超声波传播的冲击波。 2.3.2爆炸应力波的传播和特性 v1、应力波的传播和作用 v爆炸应力波传播过程和范围如图所示 图2-8爆炸应力波的传播和作用范围 R0R2 R1 图3-2 爆破的内部作用 R0-药包半径；R1-粉碎区半径；R2-破裂区半径 v冲击波的作用： 冲击波具有陡峭，以超音速进行传播，波阵面上的 状态参数突跃变化，波峰压力大大超过介质的抗压强 度，此时介质产生塑性变形或粉碎，形成压碎区，冲 击波的作用范围一般为37倍装药半径。 v应力波的作用 冲击波衰减后变成压缩波，波阵面上介质参数变化 也明显平稳，传播速度等于音速，在此范围内，应力 波峰压力小于介质的抗压强度，只能造

17、成介质的变性， 裂隙或拉断，该区域因而也成为裂隙区。作用范围为 7120150R。是爆破的主要作用区。 v地震波的作用 压缩波对介质连续做功，能量进一步损耗，直至 衰减到不对介质起任何破坏作用，只能引起质点的弹 性振动，形成地震波。地震波衰减慢，作用范围大， 对介质没有直接的破坏作用，该区域也称弹性振动区。 v2、爆炸应力波的特性 应力波特性，亦即它们在介质中的性状和变化， 首先取决于介质本身的性质。为了透彻了解应力波的 机制，必须知道介质的一般基本特性，介质的基本特 性可用应力-应变关系表示。 设理想情况下，波的传播方向与质点运动方向一致， 在固体取一微元体，扰动前处于静止状态，初始应力为

18、，质点运动速度为零u0，扰动后密度和应力变化分 别为，质点获得运动速度du，扰动传播速度为D，设t时 刻为与截面AA的波头，在t+dt时刻传播至截面CC，传 播的距离为Ddt，在同一时刻dt时间内，截面AA的质点 移至截面BB，产生位移为du、dt。 由动量守恒得： 化简得： (2-91) 质量守恒： 忽略方程中高阶微量得： （2-92） 根据比容V，密度和体积变形的关系有： dtdDdtdu dDdu dtdDdt duDd （293） 由（2-91）和（292）式得： （2-95） v v（295）式表明，固体中微幅应力波的传播速度决定 于应力随应变的变化率即变形模量。 ddv 2 1 1

19、 V d V d V dV d 2 （294） d du dV d V d d d d D 1 图2-10 介质中应力扰动的传播 A B A B CBA C c 0 图2-9 固体在冲击荷载作用下的 变形曲线 v因有限幅应力波可看作是一系列微幅应力波的叠加，要了解 固体中爆炸应力波的特性，须知道固体材料在冲击载荷下的 应力和应变关系。 v图2-9为固体在冲击和在作用下的典型变形曲线。图中0A为 弹性区，A为屈服点，AB 为弹塑性区，B点以后材料进入 类似流体状态。 v若应力值不超过弹性极限或屈服点。在固体中传播的是弹性 应力波，应力和应变为线性关系，变形模量为常数，与扰动 强度无关。弹性应力波

20、的传播速度等于未扰动固体中的声速， 按下式计算，有 v上式只适用于无侧限应力的一维平面波。 E c v当应力超过屈服点而处于弹塑性变形区内时，紧跟 弹性应力波后面传播的是弹塑性波。该区内变形模量 不是常数，而是随应力值增大而减小，因此，高应力 处的微幅应力扰动要比低应力处的微幅应力扰动传播 慢，波头在传播过程中逐渐变缓。通常弹塑性变形区 可近似用直线表示，并以该直线的变形模量计算弹塑 性的平均波速，其值小于未扰动固体中的声速。 应力值超过B点，固体的变形特性类似于流体，其变形模 量随应力增大或压缩率减小而增加。与弹塑性波相反，高应力 处的微幅应力扰动要比低应力处的微幅应力扰动传播快，其结 果将

21、形成陡峭波头（见2.4.3冲击波的形成过程）。但在BC变 形段内，变形模量仍小于弹性区变形模量，故在该区内，虽能 形成陡峭波头，但波头传播速度不是超声速的（与固体介质中 的声速比较），不能把它看作是冲击波。为与冲击波区别起 见，把此波称为不稳定冲击波（非稳态冲击波）。 v应力值超过C点后，不仅能形成陡峭波头，而且波头传播速 度是超声速的，这时将可以形成冲击波。 v装药在岩体中爆炸时，若作用在岩体上的冲击载荷超过C点 应力（临界应力），首先形成的是冲击波，尔后衰减为非稳 态冲击波、弹塑性波、弹性应力波、地震波（图2-11）。 图2-11 不同应力条件下介质中传播的各种爆炸应力波 2.4岩体内的爆

22、炸应力波 v2.4.1爆炸应力波的初始压力 炸药爆炸后作用于周围介质中的爆炸应力波初 始压力，因装药的耦合条件不同，直接或间接地作 用于周围介质，激起不同强度的冲击应力波，按声 学近似的方法或按爆轰产物等熵膨胀后碰撞炮孔壁 方法，在耦合装药情况下，初始径向应力峰值为： (297) 式中： 冲击应力波初始峰值压力，Pa 被作用介质密度，kg/m3 e p p p DC C p 00 0 2 0 p 炸药装药密度，kg/m3 介质中纵波传播速度，m/s 炸药的爆速，m/s 炸药的爆轰压，Pa 不耦合装药： 孔壁压力 式中： - 碰撞次数=8-10； - 装药半径 -炮孔半径 - 待定系数，与药包形

23、状有关。 0 D e P 0 p C 2 00 4 1 DPe nrr P P bc e )( 2 1 b r c r n v2.4.2岩体内的冲击波 冲击波通过时，岩石类似于流体，波头上状态参 数将发生突跃变化。设岩石初始参数为 ， ， ， =0 ，冲击波波速为D，波头上岩石状态参数突变为P、V、 E、u，（P、V、E、u分别为压力、比容、内能、质点 速度）这些参数满足下列冲击波基本方程： （298） （2-99） （2-100） 0 P 0 V 0 E 0 u V V uD D 0 )()( 00 uDuD 0 0 PP V Du VVPPEE 000 2 1 )( 0000 uuuDPP

24、 （298）式至（2100）式三个基本方程，四个未知 数，解上述方程必须知道状态函数E=(P,V)，这很困 难，因此一般通过实验给出岩石的雨果尼奥曲线， 即，最常采用的RH方程为： (2101) 式中：B试验确定的系数，或按下列公式计算： （2102） n试验指数，大多数岩石n4 或： a,b试验确定的常数，某些岩石值见表2-1。 1 0 n BP 2 4 1 p CB buaD v2.4.3应力波衰减规律 随着应力波的传播，由于几何衰减和能量损失造 成的衰减，随应力波传播距离增大，应力峰值将不断 减少，径向应力峰值与距离的关系可表示为： (2103) 式中： 距离装药中心距离与炮孔半径的比值

25、。称为 对比距离 应力波衰减指数 对冲击波取正，对应 力波取负。 泊松比 r P rm 1 r 1 切向拉应力峰值（绝对值）可通过径向压应力峰 值来表示： （2104） 系数B与岩石的泊松比有关，可用下式计算： 因此 （2-105） maxmaxr B 1 B maxmax 1 r 2.7 纵波在两种岩石界面的反射和折射 v 2.7.1入射波、发射波、折射波的性质 入射波：从波源向外发射，有一种介质向另一种介质 传播的波称为入射波，炸药爆炸产生的冲击波、应力 波、地震波，均以压缩波的形式从爆源向外发射传播。 向自由边界垂直和斜射入射的波，称为垂直入射波和 倾斜入射波。 反射波：入射波遇到两种不

26、同介质界面则被反射，这 部分被反射回来的波称为反射波。 折射波：入射波遇到两种不同介质的界面时，除了一 部分反射为反射波外，另一部分波透过界面进入另一 种介质的波，称为折射波（透射波）。 v2.7.2 折射定律 在OAO中，有： 在OBO中，有： 图2-17 入射波在不同介质中的折射 两边相除，即得的波的折射定律： 上式称为斯涅耳定律 1 sin OO tV OO AO i 2 sin OO tV OO BO i 2 1 sin sin V V i i 2.7.3 纵波在两种岩石界面的反射和折射 因装药在岩石中爆炸后直接激起的应力波主要是纵 波，我们主要讨论纵波在两种岩石界面上的反射和折射。

27、图2-18 纵波在两种岩石界面上的反射和折射 设两种岩石界面为YX平面，波源在第一种岩内， 其波阻抗 ，第二种岩石的波阻抗为 ，若界 面上两变岩石的粘结力不能使岩石沿界面产生相对滑 移，则应力波入射到界面时所应满足的条件是界面两 边法向位移、切向位移、法向应力、切向应力均连续 相等。为满足边界条件，在界面上将同时产生四种不 同的新波：反射纵波和反射横波、折射纵波和折射横 波。 11p C 22p C v入射角、反射角、折射角与波速关系遵循斯涅耳定 律 式中： 第一种岩石纵波和横波的波速 第二种岩石纵波和横波的波速 v各波位移中幅值满足下列方程（规定在波传播方向 产生的位移为正） v法向位移：

28、（2148） v切向位移： （2149） 2 3 2 3 1 2 1 2 1 1 sinsinsinsinsin spspp CCCCC 11,sp CC 22 , sp CC 0sincossincos)( 353423121 AAAAA 0cossincossin 353423121 AAAAA v 法向应力： （2150） v 切向应力 （2151） 由以上四个方程可求得入射波位移幅值表示的反射波和折射波 位移。 正入射时，入射波和所有其它角度均等于零，在这种情况 下 ，即没有横波产生，只产生反射纵波和折射纵波。 解方程得： （2152） 02sin2cos2sin2cos 3 1 2

29、253 1 2 242132121 spsp CACACACAA 02cos2sin2cos2sin 2 1 53 2 1 4 2 222 2 1 3121 2 11 P p p p s s p s C C A C C AC C C AAAC 1122 1122 1 2 pp pp CC CC A A 0 53 AA （2153） 反射应力R，折射应力T与入射应力I的比值为： (2-154) (2-155) 1122 11 1 4 2 pp p CC C A A 0 421 AAA 0 2241121 pp CACAA VC p 1122 1122 1 2 pp pp I R CC CC A

30、A 1122 22 11 22 1 4 2 pp p p p I T CC C C C A A 11p I I C V 11p R R C V 22p T T C V 1122 1122 1 2 pp pp I R CC CC A A 1122 22 11 22 1 4 2 pp p p p I T CC C C C A A 1122 2211 pp pp I R CC CC V V 1122 11 2 pp p I t CC C V V 如果将炸药看作是第种介质，岩石看作是第种介 质，炸药内传播的爆轰波为弹性波，利用(2-155)求得 的折射压力： （2156） 式中： P1为炸药的爆轰压力

31、： 为炸药冲击阻抗， 为岩石的波阻抗。 pm e C D PP 0 12 1 2 2 01 4 1 e DP e D 0 pmC 或：应力波垂直入射时，情形简单。满足连续和平衡 条件：在边界的两侧，其应力必须相等；垂直于边 界面方向的质点速度必须相等。在单位面积下有： RIT VVV RIT 1122 1122 1 2 pp pp I R CC CC A A 1122 22 11 22 1 4 2 pp p p p I T CC C C C A A v讨论： v 、 、 说明反射波为拉伸波， 质点运动与入射波质点运动相同，而与反射波传播方 向相反，透射波的应力幅值小于入射波。若第二种介 质为空

32、气，或界面为自由面时，空气的波阻抗可忽略 不计（ 0），则有： 、 、 、 即反射波的应力和质点速度值与入射波的相同，但应 力极性倒转；自由面上的应力为零，质点速度为入射 波速度的2倍。 1122pp CC 0 R 0 R V 22p C IR IR VV 0 T IT VV2 v(2)若 ，即在阻抗匹配情况下，则有： 、 、 、 说明没有反射波产生，入射波 就像不存在有界面一样直接通过。 根据能量守恒定律，反射波与折射波的能量总和 应等于入射波的能量。因此，在阻抗匹配的情况下， 入射波能量也将全部传给折射波。 v(3)若 ,反射波应力为正,而质点速度为负， 说明反射波为压缩波,质点运动方向与

33、入射波质点运动 方向相反,而与反射波传播方向相同。若第 种介质为 刚体或界面为刚性界面， 与 相比可忽略不记， 或认为 ，则有： 21 21PP CC 0 R 0 R V IT IT VV 21 21PP CC 11p C 22p C 2 2P C 、 、 、 即反射波的应力和质点速度值与入射波的相同，但质 点速度极性倒转；界面上的质点速度为零，应力为入 射波应力的2倍。 IR IR VV IT 2 0 T V 2.8应力波的干涉 两个或多个应力波相遇发生相互作用的现象称为应 力波的干涉。 v2.8.1应力波干涉的类型 v应力波干涉的情况是多种多样的，其基本类型有以下 几种： v按两个干涉波的种类区分有：两个纵波、两个横波 或一个纵波与一个横波之间的干涉。 v按两个干涉波的来源区分有：有个独立波源产生应 力波的干涉，例如相邻两装药在岩体内爆炸激起应力 波的干涉