固体中的应力波

固体中的应力波李清中国矿业大学（北京）精品文档参考书：1 王礼立 . 应力波基础 第 2 版(2005 年 8 月 1 日) ，国防工业出版社2 李玉龙 . 应力波基础简明教程第 1 版 (2007 年 4 月 1 日) ，西北工业大学3 丁启财 ( 美国 ). 固体中的非线性波，中国友谊出版公司4 宋守志 . 固体中的应力波，煤炭工业出版社5 杨善元 . 岩石爆破动力学基础 ，煤炭工业出版社6 莱茵哈特 ( 杨善元译 ). 固体中的应力瞬变，煤炭工业出版社7 徐小荷 . 冲击凿岩的理论基础与电算方法 ，东工出版社8 郭自强 . 固体中的波，地震出版社。1 欢迎下载精品文档目录第 0 章绪论 .11波动现象.12应力波的概念 .13应力波分类 .34应力波理论与其它力学理论的关系.35应力波理论的发展 .36应力波理论在岩土工程中的应用.3第 1 章一维应力波基础 .41.1波动方程及其解 .41.1.1一维纵波的波动方程 .41.1.2波的传播速度 .41.1.3波动方程的解 .51.1.4解的物理意义 .61.2应力波的几个基本参量 .71.3应力波的能量 .71.4波的衰减 .81.4.1原因 .81.4.2度量 .81.4.3衰减率的测定 .91.5考虑杆的横向效应的波动方程.101.6杆中的扭转波与弯曲波 .121.6.1扭转波 .121.6.2弯曲波 .13第 2 章二维和三维弹性波理论基础.142.1弹性体的运动微分方程 .142.2弹性体的无旋波与等容波.152.2.1无旋波 ( 纵波、 P 波 ) .152.2.2等容波 ( 横波、 S 波 ) .162.3平面波的传播 .172.3.1平面纵波 (V/c) .172.3.2平面横波 (V c) .182.4薄板中的应力波 .192.4.1控制方程 .192.4.2纵波 .202.4.3横波 .212.4.4各种波速关系 .212.5球面波 .222.5.1波动方程及其解 .222.6柱面波 .23第 3 章应力波的相互作用 .243.1一维应力波在界面的反射和透射.243.1.1应力波在不同介质界面的反射和透射.253.1.2应力波在变截面杆中的反射和透射.26。2 欢迎下载精品文档3.2两杆相撞的入射波 .273.3传播图与状态图 .293.3.1传播图 . .293.3.2状态图 .303.4弹性杆中波的传播（图解法举例）.323.4.1冲锤撞击杆件应力波的传播.323.4.2双圆柱活塞撞击钎杆应力波传播.333.5平面波的边界效应 .363.5.1平面波在界面上的垂直入射.363.5.2平面波在界面上的倾斜入射.373.6应力波引起的破裂 .393.6.1金属丝冲击波拉伸断裂.393.6.2Hopkinson压杆与飞片 .413.6.3断裂准则.413.6.4简单反射拉伸波引起的层裂或剥裂.423.6.5物体形状对应力波引起破裂的影响.453.7冲击波基本问题.45第 4 章固体中的非线性波基础 .484.1弹塑性加载波及其相互作用.484.1.1强间断弹塑性波的迎面加载.484.1.2弱间断弹塑性波的迎面加载 .504.2卸载波的控制方程和特征线.51第 5 章 岩石动态力学性质与应力波的相互作用.535.1岩石动态本构关系与动态强度.535.2岩石动态力学参数测试 .545.3本构关系对应力波传播的影响.545.4应变率相关的应力波理论.555.4.1Voigt体. .555.4.2Maxwell体 .55第 6 章 应力波在岩土工程中的应用.556.1应力波在冲击凿岩中的应用.556.1.1冲击凿岩的应力波的传递.556.1.2凿岩机的凿入机理 .556.1.3入射波形对凿入效果的影响.566.1.4冲击凿岩的破坏原理.566.2应力波在爆破工程中的应用.566.3应力波在土动力学中的应用.566.3.1绪论 .566.3.2土的动应力 - 应变关系及其描述.586.4应力波在地震工程学的应用.58第 7 章 应力波测试分析技术简介了解 .607.1膨胀环测试技术.607.2 Hopkinson杆测试技术 .607.3 Taylor 圆柱测试技术 .617.4高速冲击载荷的实验技术.61。3 欢迎下载精品文档第0章绪论1 波动现象波动现象：水波、声波、电磁波、光波等。波是一种扰动或状态在介质中的传播，波动是非常普遍存在的一种运动形式，一般可分为两大类：机械波和电磁波。这里所述的应力波属于机械波，是机械扰动在连续介质中的传播过程。机械波产生于可变形介质的强迫运动，通过质点在平衡位置附近的振动来传递能量。2 应力波的概念介质的某部分受力发生了一种状态的扰动，离开初始平衡位置，与相邻介质质点发生相对运动 ( 变形 ) ，并和周围介质产生压力差，这种压力差将导致周围介质质点投入运动，但由于介质质点具有惯性，而使某相邻质点运动滞后，外载荷在表面上的扰动就这样在介质中由近及远地传播出去而形成应力波。应力波理论主要研究力、位移、速度等物理量在固体中传播的规律以及它们对固体的作用效应。理论力学中，物理被认为是刚体( 不变形 ) ，遵循牛顿惯性定律F=ma材料力学、弹性力学，研究物理变形，但不考虑变形而产生的物理运动，不考虑物理的惯性，遵循虎克定律现实的物体，惯性和弹性兼而有之，当它受力时，既改变它的速度又改变它的形状 。物理受力部位的质点， 克服惯性，发生速度的变化， 这种变化遵循惯性定律 ( 牛顿定律 ) ，速度的变化必然导致变形，变形阻碍速度变化；反过来说，物理受力部位，由于弹性的作用，必定会有变形， 这种变形符合虎克定律， 但在实现变形时， 质点会出现变速运动，变速运又阻碍变形的发展。由此可见，物理内部同时存在着弹性和惯性，相互作用，导致物理中形变和速度的转移，这就是应力波。应力波得以在连续介质中传播的基本条件是介质的可变形性和惯性 。对于不可变形的刚体，局部的扰动（力或位移）可立即传到整个物体的每一部分。若介质没有惯性，则扰动的传递也是瞬时完成的，一切实际材料都具备这两个条件，所以一切实际材料都能传播应力波。固体中的应力波的研究主要用于地震、爆作、高速撞击、爆破、超生波等应力波的发生和传播过程。应力波波阵面 介质中扰动的区域和扰动未波及的区域的界面。1 欢迎下载精品文档分析波阵面的前后状态参量的变化关系，有两种类型。间断波波阵面 前后质点微团的状态参量有一个有限的差值。状态参量发生跃变，数学上叫强间断。连续波波阵面 前后质点微团的状态参量的差值为无限小。 状态参量的分布是连续的，数学上叫弱间断。强调一点间断波和连续波是相互转化。弥散波：如介质的性质使得高应力水平增量波具有较低传播速度，波形在传播过程中会逐渐拉长、散开的连续波。汇聚波：如介质的性质使得高应力水平增量波具有较高传播速度，那么处于后面的高波速的增量波不断追赶前面的较低波速的增量波，使得连续波波形逐渐缩短。冲击波：一定条件下， 后面具有高波幅的增量波赶上前面波幅的较低的增量波形成以统一波速传播的强间断波波阵面，连续波转化为冲击波。间断波中除了冲击波之外，还有一种等熵的间断波，这就是弹性间断波，因为弹性变形是可逆的过程， 弹性间断波只是在波形上与连续波不相同，二者在本质上没有区别。最后介绍关于 加载波与卸载波的概念。固体介质不但能承受压力， 而且能承受拉力。对介质加压，使介质压密就是加载；对已经受压后的介质减压，使介质稀疏就是卸载。当波阵面通过一个介质微团时，其效果是使微团压密的就是加载波（压缩波） ；其效果是使微团稀疏的就是卸载波（拉伸波） 。加载波和卸载波的波形如图示。2 欢迎下载精品文档3 应力波分类(1) 按力的特征分 拉伸、压缩波（稀疏波或纵波） ；弯曲波、剪切波（横波）(2) 按波阵面的形状分 平面波、柱面波、球面波(3) 按变形特征分 无旋波（膨胀波） 、等容波（畸变波）(4) 按介质的物理特征分 弹性波、塑性波、粘弹波、粘塑波(5) 按介质的几何特性分 一维波 （杆波）、二维波（平面波）、三维波（空间波）4 应力波理论与其它力学理论的关系应力波理论是固体动力学的分支。 但目前的固体动力学往往集中研究材料在高应变率下的动态力学性能，而把已知材料的动态力学性能、介质受到外部动载作用的规律研究让位于应力波理论。但二者是相互依赖而发展，一方面应力波理论的发展必须建立在对材料动态力学性能的了解之上；另一方面，材料的动态力学性能往往必须通过应力波的测试与分析才能得到。应力波理论与其它力学理论的区别(1) 动力学研究载荷的 早期效应 （瞬时效应，着重研究质点的运动和变形等物理量随时间的变化过程以及在物理中的传递） ，静力学研究的是载荷的 后期效应 （只研究在力的作用下达到平衡之后的状态） ；(2) 动力学研究载荷对介质的 局部效应 ，静力学研究载荷对介质的 整体效应 ；(3) 动力学研究的动载有明显的 耦合效应 ，静力学研究的静载作用于固体的 应力分布不随介质而变。5 应力波理论的发展线弹性波传播的数学理论早在上个世纪中叶由柯西(Cauchy) 、泊松 (Poisson) 、斯托克斯 (Stokes) 等解决，可直到本世纪四十年代，由于电子技术的发展，人们才直观地在固体中“看到”波，应力波理论才开始在一些工程领域得到应用。与此同时，Donnell 、Taylor 等人在理论上又发展了塑性波理论。五十年代前后， 考虑应变率效应的粘塑性波理论又得到了发展。应力波理论特别在地球物理勘探中的“实时采集与数据处理”技术得到了迅速发展。6 应力波理论在岩土工程中的应用爆破工程、凿岩工程、桩基工程、岩石动态力学、土动力学、地震工程与抗震工程、地球化学勘探。3 欢迎下载精品文档第 1 章一维应力波基础1.1 波动方程及其解一维纵波的波动方程如图 1-1 所示，在一等截面的一维杆中取一微段dx，截面面积为 A。基本假设为杆的横截面在变形过程中保持平面，不考虑横向扩展效应，杆上只分布沿截面均匀分布的轴向应力，因而 位移 u、工程应变、质点速度 v 和应力都只是 x 和 t 的函数。其左截面 mn与右截面 m n 上的作用力 F、 F分别为F AAEAEu ， vu1uxxtxx ttFFF dxAEuAE2 u2 dx1xxx微元体 dx 所受的惯性力为： maA dx2 ut 2由牛顿定律 ，得微元体平衡方程FFA2udx2tAE2u dxAdx2ux2t 2即2 u E2uc 22 u , cE（1-1 ）t 2x 2x方程 1-1 是一维纵波的波动方程，c 为纵波的波速。波的传播速度波动方程中的 c 为称为波速， cE （纵波波速）。根据能量守衡定律 和冲量定理可以推导。4 欢迎下载精品文档取一单位面积细长的杆，一端受到撞击，撞击后杆端的初速度为v0，受力为 0。经过时间 t 后，扰动扩展到长为 l 的区域，在此范围质点的运动速度均为v0 ，内力为 0，则时间 t 内有外力作功w0l0 v0t扰动区域的动能Ed1 lv022扰动区域的 位能El102lE2根据能量守衡定律有0v0t1 lv02102l E（ 1）22根据冲量定理有0tlv（ 2）所以有v01 (v002 E v0 )2v00（ 3）E可见细长杆受撞击后，力与质点速度v0 并非无关，而是成正比。将（ 3）式代入（ 2）式，得到纵波的传播速度l0Ecv0t同样可以弹性横波的波速cTGE2(1 )以上波速推导中，应用了线弹性本构关系E ，事实上只要( ) ，既应力是应变的单值函数，而与应变率无关，波动方程（1-1 ）就成立 。1.1.3波动方程的解方程（ 1-1 ）的解法有分离变量法（驻波法） 、积分变换法及行波法等，其中行波法对求解波动方程最为有效。令x ct 、x ct ，则 x、 t，故22c。5 欢迎下载精品文档u(x,t )u(2,)w(,)2cuwwwwxxx2u(ww )(ww)x2xx2 w2 w2 w2 w（ 1）22uwtwc(ww )tt2uc( ww )c( ww )t 2tt2 w2 w2 w2 w（ 2）c2(22 )将（ 1）和（ 2）式代入波动方程，得2 w（ 3）0对（ 3）式先对积分，得wf ( )对（ 3）式先对积分，得wf ( ) dg( )f ( ) g ( )于是有u=f(x+ct)+g(x-ct)（1-2 ）公式（ 1-2 ）是一维波动方程的通解。解的物理意义在某范围内（ - l ，l ）的扰动，在 t 时间之后，分为右行波和左行波。右行波的前后沿分别 沿着 直线 x kt l 和 x kt l 前进，而 左行 波的 前后 沿分别 沿着 直线 x kt l 和 x kt l 前进，经过 t 时间之后，分解右行波和左行波的形状不随时间和dx位置的变化而变化，前进的速度均为|kc 。6 欢迎下载精品文档1.2 应力波的几个基本参量右行波 u=f(x-ct)：（1）应变uf ( xct )x（2）质点速度tc2vuf ( x ct ),v（1-3 ）（3）应力cv（1-4 ）（4）波阻率c ：表示单位质点速度所需要的应力如v，v101210 32000，=10m/sc5000500E2.1 1062 1034200 kg / cm2 ， FA21顿 （B25钎杆）对于钢材7800N / m35100m/s39.78106Pa/(m/ ) 39.78(N/mm2) /(/ )ccsm sg10m / s2部分岩石的应力波参量如下表所示。部分岩石的应力波参量岩石纵波速度 m/s岩石比重32T/m波阻率 (N/mm)/(m/s)砂岩140040002.55310页岩140030002.3037大理岩350060002.65916石英岩500055002.651316.5板岩350055002.65914.5花岗岩300050002.65813玄武岩450060002.501318岩石的波阻率只有几到十几，远比钢材要低，表明岩石中产生单位质点速度所需要的应力比钢材中要低得多。1.3 应力波的能量一维杆中位能（弹性能）11A21Ac 2123dwe2Adx2Ecdt2Edt2AEcdt一维杆中动能。7 欢迎下载精品文档dw d1Adxv21Acdt (c ) 21 A c 32 dt1 AEc2 dt1Ac 2 dt22222EwdwAEc2dtAc2dt(15)0E 0dx 微段所具有的总能量为Ac 2 dtdwdwedwdAEc 2 dtE应力波的总能量为wdw AEc2 dtAc2 dt（1-5 ）0E0结论（ 1）波动过程中，对任一微段在任一时刻具有的动能和位能都相等 ，这是波动和振动的本质区别，一个作纯振动的系统，尽管总能量不变，但总存在着动能和位能的转换。结论（ 2）任一微段所具有的总能量是应力或应变的函数，应力和应变是一种波动过程，因此 能量也是一个波动过程 ，对任一微段来说，其能量是不守衡，沿着波传播的方向，该段波源获得能量，使其能量逐渐增大，又逐渐把自身能量传递给后面的介质，能量随着波动过程而有规律地传播， 波动是能量传递的一种方式 。1.4 波的衰减原因真实物质很少是理想的弹性体，而常常是弹塑性或粘弹性等。当波在粘弹性介质中传播时，因存在内摩擦，将产生能量的损耗；当波在热弹性体中传播时，在应力波通过时，固体一部分受压，另一部分发生膨胀，压缩部分温度升高，膨胀部分温度降低，这种温度梯度的出现，将在固体中引起热的传递，并伴随着不可逆过程的发生，使应力波因热耗散而发生衰减。总之波的衰减来源于内摩擦和外摩擦的作用：内摩擦由材料的粘弹和热弹性决定；外摩擦决定于材料所处的工作环境。度量度量波的衰减程度通常采用衰减系数、损耗因子-1。Q（ 1）衰减率振幅衰减率、应力衰减率d /（1-6 ）dx表示应力波在单位长度上的振幅衰减。ddx ，0 e x（1-7 ）。8 欢迎下载精品文档其中， 0 为在x=0 处的应力波幅值； 为在任一位置x 处的应力波幅值， 相当于前面的 max。（ 2）损耗因子 Q-1Q品质因素Q 11w（1-8 ）2w式中w一次应力循环所损失的能量；w应变达到最大时所贮存的弹性能。（ 3）能量衰减率2 d ) /(2dw / w()2dEEdxdxdx2（1-9 ）由此可见能量衰减率是振幅衰减率的两倍。（ 4）损耗因子 Q-1 与衰减率的关系Q 11w 11c Tcc（ 1-10 ）2w 222 ff衰减率的测定0e x ，e x ， lnx ， ln0x ，1 ln01 ln000xn1ln0（ 1-11 ）2nl但在岩土工程中， 品质因素 Q或损耗因子 Q-1的应用更为普及， 大量科研表明损耗因子 Q-1 对岩石物理性质变化的反映比声速更明显和可靠。由于衰减系数在很宽的频率范围内是线性函数，因而Q-1 与频率无关，这时 Q 1r c 。损耗因子的测试大多采用频谱振幅比法，该方法涉及到频谱分析技术，比较复杂。参考应用声学 1987（ 4）或地震研究 1983（ 4）等有关文献。9 欢迎下载精品文档1.5 考虑杆的横向效应的波动方程前面讨论的一维应力纵波理论都假设杆的平截面在变形后仍保持平截面，并在平截面上只作用着均匀的轴向应力x。这时实际忽略了杆中质点横向运动的惯性作用，即忽略了杆的横向收缩或膨胀，因而是一种近似理论，通常称为初等理论或工程理论。由材料力学可知，广义虎克定律如下x1 x(yz ) ， y1 y( zx ) ， zEE当一维杆受x 作用时，即y= z=0 时1x ， zxx ，yxxEEE即1 Ez(xy )xu yx yy u x ， uzx zz uxxxu yyv x ， vzuzzxzv xvytxttx单位体积的平均横向动能为11dydz(vy2vz2 )12 rg2 (tx ) 2A A 22式中r g截面对 x 轴的转动半径， r g21( y 2z2 )dAA A上面在推导一维波动方程时，是从分析杆中微段的受力着手，在运动着的微元体上作用着一对静力平衡的力A和一非静力平衡的力Adx 。非静力平衡的力与微元体x的纵向惯性有关，其所作的功转化为微元体的纵向动能，单位时间所作的功等于纵向动能的增加率，即Adx v(1Adxv2 )xt2整理后，微元杆端的运动方程为v（ 1-12 ）tx一对静力平衡的力 A所作的功，在初等理论中全部转化为应变能。在计及横向运动的情况下，下可看作由两部分组成：一部分使微元体应变能增加 ，另一部分转变成了杆的横向动能 。这样，就单位时间、单位体积而言，有(1E 2)t( 12 rg2 ( )2 )tt22t。10 欢迎下载精品文档2 rg22Et 2（ 1-13 ）当横向动能相关的第二项可忽略时，上式可以作这样的理解： 在考虑了横向惯性后，Hook 定律应被上式所表示的新的应力应变关系所代替。2既然横向修正项与t2 成正比，显然只有在载荷随时间有十分显著变化的情况下，这一修正才是必要。将（ 1-13 ）式代入（ 1-12 ）式，并将u 代入，得x2uE2 u2 24 u（ 1-14 ）t 2x 2gtx22与（ 1-1 ）对照可知，等号右边的第二项代表横向惯性效应，有了这一项，杆中的弹性纵波将不再如初等理论中那样恒速c0E 传播，而是对不同频率 f （或波长）的谐波将以不同的波速（相速）c 传播。取 uAe ik ( x ct ) ，式中 k 为波数， k2 ； c 为相速度，代入式（ 1-14 ），得2222k 2rgk 4 c2k c0c02c02化简得(1 k 22r g2 )( c ) 21c0令 kk rg ， cc （相对波速），则c0(1k 2 )c 211c(1k 2 )1 / 2。11 欢迎下载精品文档当 k 2k 22 rg21时，近似有c1(k g )21 222(r g212)c0对于半径为 R 的圆柱杆， rgR，得2c12 2(R)2（ 1-15 ）c0式（ 1-15 ）是考虑横向惯性修正的Rayleigh 近似解。当 k0.3 时，可不考虑横向效应，对圆截面杆，有Rd2r g2， =0.29 ， k22d20.290.322即2.1527d由式（ 1-15 ）可知，高频波（短波）的传播速度低，而低频波（长波）的传播速度较高。对于线弹性波来说，既然任意波形的波总可看作由不同频率的谐波分量迭加组成，而不同频率的谐波分量现在将各自按自己的相速传播， 因此波形不能再保持原形而必定在传播过程中分散开来，即发生所谓波的弥散 ，又叫几何弥散，不同于非线性本构弥散和粘性弥散。1.6 杆中的扭转波与弯曲波杆中的横波包括扭转波和弯曲波。扭转波杆中的扭转波波动方程的推导与纵波波动方程（1-1 ）完全类似，形式也相同。2G22 ct22G（ 1-16 ）t2xx2 , ctct 为弹性扭转波波速，与无限介质中的剪切波波速相一致。12 欢迎下载精品文档因为 ctE， =0 0.5 ，所以 c2(1)2 3 ，即纵波传播速度是2(1)ct扭转波或剪切波的1.4 1.7倍。钢材中的纵波速度为5100m/s，而剪切波为 3220m/s。需要指出的是： 在圆截面杆中，不同频率的扭转波都以相同的相速ctG ，不发0生弥散现象。弯曲波杆中的弯曲波波动方程为2 wEI4 w（ 1-17）t 2Ax 4弯曲波的波速是变化的，波形也是变化，是一种比较复杂的波。不作要求。13 欢迎下载精品文档第 2 章 二维和三维弹性波理论基础 2.1 弹性体的运动微分方程