![水工建筑物抗震计算基础](https://wfqqreader-1252317822.image.myqcloud.com/cover/484/40937484/b_40937484.jpg)
1.6 地震波的衰减
1.5节介绍了线弹性介质中传播的地震体波和面波,没有考虑使波衰减(Attenuation)的任何因素,波幅因此不会发生改变。实际材料中波的传播不会发生这种情况,地震波是随着传播距离的增加而衰减的。这种衰减归因于两方面:一方面是岩土介质本身存在材料阻尼而要耗散一部分波动能量;另一方面则是存在几何扩散效应,使得单位体积的弹性能(称为比能)随传播面向远处的扩散而减小。
![img](https://epubservercos.yuewen.com/0C0364/21277078208321506/epubprivate/OEBPS/Images/txt001_164.jpg?sign=1739542046-Q8JgfYJRAzYidXtkiQaVul9R5pq2Xh7e-0-c9cbb9eb20703a0770e7dfd24f22a148)
图1.6.1 模拟土体的开尔文-伏格特(Kelvin-Voigt)模型
1.6.1 材料阻尼
在实际地质材料中,传播的行波所携带的弹性能总是有一部分转换成热能,使得波幅逐渐衰减。黏滞阻尼(Viscous damping)因其数学上的方便,常用来表示弹性能的这种耗散(即转换成热能)性质。因此,在黏弹性波动分析中,岩土体通常被模拟成所谓的开尔文-伏格特(Kelvin-Voigt)型固体材料,这种本构模型是最常用的流变模型(Rheological model)之一。如图1.6.1所示,该模型是由弹性元件(弹簧)和黏性元件(阻尼器)并联而成,土在动力作用下的应力由弹性恢复力和黏性阻尼力组成。
由于地震波的主要能量由剪切波所携带,因此通常考虑剪切波的作用效应。为方便计算,考虑SH波。考虑图1.6.1中的薄土层单元,土体剪切应力τ与剪应变γ、剪应变率(剪应变随时间的变化)的关系为
![img](https://epubservercos.yuewen.com/0C0364/21277078208321506/epubprivate/OEBPS/Images/txt001_166.jpg?sign=1739542046-1NLaQdQ24zlws8s4T7U29e11TyKrBeTg-0-ffe9209fe32e37d1219efa5ba2c2c2f3)
式(1.6.1)的物理意义是,总剪应力τ可归为弹性部分(正比于剪应变γ)与黏性部分(正比于剪应变率之和,η为介质黏性阻尼系数。若η=0,则τ=Gγ,满足胡克定律。这里假定土的剪切模量G和阻尼系数η是常数(实际上,剪切模量与阻尼系数随土的动应变幅值而变。一般地,剪切模量随剪应变幅值增加而降低,阻尼随剪应变幅值的增加而增加)。
若剪应变具有简谐形式:
![img](https://epubservercos.yuewen.com/0C0364/21277078208321506/epubprivate/OEBPS/Images/txt001_168.jpg?sign=1739542046-dLPU7nefsgAt35ZVRFRsUaIZBGOWsEOw-0-0b23320ffaa5f86cca899de930bc9dfb)
则剪应力为
![img](https://epubservercos.yuewen.com/0C0364/21277078208321506/epubprivate/OEBPS/Images/txt001_169.jpg?sign=1739542046-GxGJZShabR7G0QDGbAEFVHKI9t6w5Z6p-0-7bef223b82b89a1e283f67467206db2a)
联立式(1.6.2)、式(1.6.3),可得
![img](https://epubservercos.yuewen.com/0C0364/21277078208321506/epubprivate/OEBPS/Images/txt001_170.jpg?sign=1739542046-UBEvXrbjmy9qHpOfjdKfpYFODd6EiSNW-0-b5c4fe86e20d07ee95e0585de6f23d9c)
式(1.6.4)表明,在一个振动循环中,应力τ与应变γ所围成的图形为如图1.6.2所示的椭圆。
![img](https://epubservercos.yuewen.com/0C0364/21277078208321506/epubprivate/OEBPS/Images/txt001_171.jpg?sign=1739542046-PoQae0cJzQ7In836bxX1ynqNKAD3XLdC-0-22f56cadba6a6c0fbbc9f305062bbd56)
图1.6.2 应力-应变滞回圈与阻尼比ζ的关系
上图中的椭圆就是通常所说的滞回圈。在一个振动循环中所耗散的能量可由椭圆面积给出,即
![img](https://epubservercos.yuewen.com/0C0364/21277078208321506/epubprivate/OEBPS/Images/txt001_172.jpg?sign=1739542046-jrmflNmS8ls63SyWcFYTRoWmKZztzSVN-0-7b94f20c819501bc8e9fe38e260a52a7)
式(1.6.5)表明,耗散的能量ΔW与加载频率ω成正比。但实际土体是通过晶粒间的相对滑移来实现能量的耗散,并不对加载频率ω有敏感性。
当剪应变达到最大值γ0时,剪应变率为零(对应速度为零,则动能为零),则在一个循环中所储存的最大弹性应变能(势能)为
![img](https://epubservercos.yuewen.com/0C0364/21277078208321506/epubprivate/OEBPS/Images/txt001_173.jpg?sign=1739542046-eAkeI4lHvkJ2dm5BBy7XYiK45CHs26q2-0-c7fccb4612452ec02836ad4ffe1ac77b)
引入按式(1.6.7)定义的阻尼比ζ:
![img](https://epubservercos.yuewen.com/0C0364/21277078208321506/epubprivate/OEBPS/Images/txt001_174.jpg?sign=1739542046-JOdBrvfdHhGKPIJ7y15okeQPm5Gf3tJR-0-7c76f786761456fc639a4e2a951895b3)
因此,将式(1.6.4)与式(1.6.5)代入式(1.6.7)得
![img](https://epubservercos.yuewen.com/0C0364/21277078208321506/epubprivate/OEBPS/Images/txt001_175.jpg?sign=1739542046-hwdttFwK8pKFzYdYHyl38B4HBdswzDbd-0-22e9c5edbee7dc1211e1c05d4de00fa5)
从而有
![img](https://epubservercos.yuewen.com/0C0364/21277078208321506/epubprivate/OEBPS/Images/txt001_176.jpg?sign=1739542046-9GMgoB0qjOqNmuvSxOovzHEnJIKZCWlY-0-f9dce952d45521b2b8cfadf35560ea97)
代入式(1.6.5),得
![img](https://epubservercos.yuewen.com/0C0364/21277078208321506/epubprivate/OEBPS/Images/txt001_177.jpg?sign=1739542046-T4MySt1tkeggamwvjfMJkQE30jUbNgt5-0-a668b9dc22e5259577f09baaab410524)
式(1.6.10)表明,若黏性阻尼系数η取式(1.6.9)表示的值,则耗散的能量ΔW就与加载频率ω无关了。
如图1.6.3所示,考虑沿竖向向上(+z向)传播的剪切波,假定仅引起水平x方向的位移u(z,t)。
根据运动时土体单元受力的平衡,有
![img](https://epubservercos.yuewen.com/0C0364/21277078208321506/epubprivate/OEBPS/Images/txt001_178.jpg?sign=1739542046-J3BF7782Nw2VIY5tHSvGZ7bDQGAj9KJN-0-fb7e7e80019945794102f437299b823f)
将式(1.6.1)代入,并注意到剪应变γ=u/z,则
![img](https://epubservercos.yuewen.com/0C0364/21277078208321506/epubprivate/OEBPS/Images/txt001_179.jpg?sign=1739542046-eOpiY1GBavG66yuLg6nndE9bhrcsZhKJ-0-48514487f99de7129b3657e714af842f)
这是介质满足黏弹性应力-应变关系[如式(1.6.1)]的波动方程。若η=0,则波动方程变为式(1.5.17)的标准形式。
![img](https://epubservercos.yuewen.com/0C0364/21277078208321506/epubprivate/OEBPS/Images/txt001_180.jpg?sign=1739542046-tX2ST94faBea9p7sTCN5lkWW94ZX8kyC-0-c20a2c29d080b286d0ef67d9bcfebe8d)
图1.6.3 剪切波作用下土体单元受力示意图
假定式(1.6.12)的位移运动为简谐形式,同1.5.2.2节里的分析类似,这里也采用复数表示法(仅为了推导方便,以避免采用三角函数表示法所带来的表达式过于冗长),利用欧拉公式有
![img](https://epubservercos.yuewen.com/0C0364/21277078208321506/epubprivate/OEBPS/Images/txt001_181.jpg?sign=1739542046-E9dEH5fHfrlzojCa8WyIrVZLK0dfIegU-0-4f19a42cbb9ff8770ac93e029f9b44a6)
式中i=为虚数单位,U(z)为沿深度变化的幅值。将式(1.6.13)代入式(1.6.12),得
![img](https://epubservercos.yuewen.com/0C0364/21277078208321506/epubprivate/OEBPS/Images/txt001_183.jpg?sign=1739542046-mhq0tCZtpgHhekigkYjAWMCPKBm2QUGZ-0-de8c2d252a68d37cd3a74f5568f15285)
式中G*=G+iωη,一般称为复剪切模量。解上面的标准常微分方程,可得U(z),再代入式(1.6.13),可得位移为
![img](https://epubservercos.yuewen.com/0C0364/21277078208321506/epubprivate/OEBPS/Images/txt001_184.jpg?sign=1739542046-pmkpZe2q9XMPq2CXwOHoaNcvS7UvOJtR-0-10a63e8386c867196e9002fcba8187ad)
系数A,B由边界条件确定,而k*=一般称为复波数,可写成
![img](https://epubservercos.yuewen.com/0C0364/21277078208321506/epubprivate/OEBPS/Images/txt001_186.jpg?sign=1739542046-oQrmHoJEqNJmld4xegb1scBBHNWV9LhC-0-b232f686d8c105a6a68ed8d3d40478c5)
其中
![img](https://epubservercos.yuewen.com/0C0364/21277078208321506/epubprivate/OEBPS/Images/txt001_187.jpg?sign=1739542046-2yTKI4gKKFKaDMDErUkR3XZ4sgay8tMk-0-16b66bcae447449536726dc6e304670d)
仅正的k1和负的k2具有物理意义。实际位移可只取式(1.6.15)的实部。注意,对于无黏性情况,η=ζ=0,k2=0,k1=k。
对于一个沿+z向传播的波,位移式(1.6.15)的第一项可写成
![img](https://epubservercos.yuewen.com/0C0364/21277078208321506/epubprivate/OEBPS/Images/txt001_188.jpg?sign=1739542046-5owSNXS62SOUdeBuXiIpa4wOSVzKyJ3E-0-8c231e36225d288063b6c6cf10925ecf)
由于k2<0,式(1.6.18)意味着材料阻尼ζ的存在使得波幅A随距离的增加而以指数ek2z规律衰减。
1.6.2 辐射阻尼
地震波从震源向外传播,即使介质没有材料阻尼起衰减作用,单位体积的能量也随着传播范围的扩大而逐渐衰减,这就是所谓的介质辐射阻尼(Radiation damping,也称几何扩散阻尼)在起作用。从震源向外传播得越远,能量就衰减得越多。下面通过一个纵波衰减的例子对此进行说明。
如果发震断层破裂范围很有限,断层区尺寸仅约数千米,那么震源可用一个点源来表示。从震源发出的波将沿所有方向以波速c向外传播,波动的前沿(称为波前)是一系列不断扩大的球面(图1.5.5)。在球半径r足够远处,波前可认为是平面。在球内部取一微小夹角为α的无限长楔形体,在r处取一长为dr的微元ABCD,如图1.6.4所示。由于α足够小,可认为在AB和CD面上作用的正应力(此处考虑的是纵波)为均匀分布。假定质点的径向位移为u(r,t),以向外扩散传播为正。这里不考虑材料阻尼,可列出微元ABCD在半径方向的运动平衡方程:
![img](https://epubservercos.yuewen.com/0C0364/21277078208321506/epubprivate/OEBPS/Images/txt001_189.jpg?sign=1739542046-F2m1tdMARdmtgfEzTL8kWvgM7nJSe5tG-0-69488a95ff2489248d6871f2e1a98ac0)
![img](https://epubservercos.yuewen.com/0C0364/21277078208321506/epubprivate/OEBPS/Images/txt001_190.jpg?sign=1739542046-y4zzZYODFNqcRRF8lm30jsBViB6DCcqJ-0-486554f960202515863fa816acb5396c)
图1.6.4 辐射阻尼的存在(以点源发出的纵波向外扩散传播为例)
化简式(1.6.19),并由应力-应变关系(图1.5.4):
![img](https://epubservercos.yuewen.com/0C0364/21277078208321506/epubprivate/OEBPS/Images/txt001_191.jpg?sign=1739542046-UFs9pxMda99FYfL7qEaTcZXrQC3BS0Kj-0-617618e69dfad0584ad89062ac03a2a9)
得
![img](https://epubservercos.yuewen.com/0C0364/21277078208321506/epubprivate/OEBPS/Images/txt001_192.jpg?sign=1739542046-IoKHUHcQ9NnPKQd7kAAqbuUz0zTDVEow-0-ba93792897682e5f271792b61b9ec51c)
式(1.6.21)正是标准波动方程,其解具有形式
![img](https://epubservercos.yuewen.com/0C0364/21277078208321506/epubprivate/OEBPS/Images/txt001_193.jpg?sign=1739542046-mz4mY2Amx5B7FDeWAjS7tcPM3ixa0sHY-0-d5de4c75ab662e54381c73531674c718)
由于考虑的是从震源向外传播的波,因此,舍去表示向内传播的波g(r+ct),位移为
![img](https://epubservercos.yuewen.com/0C0364/21277078208321506/epubprivate/OEBPS/Images/txt001_194.jpg?sign=1739542046-7gMMmd3yIjib1bt0IC8OJ4b5Wp9S6ZTw-0-e7918ff267cb3f7d4ec853476de6fffd)
式(1.6.24)表明,位移随传播距离r的增加而减小。这就相当于尽管没有考虑介质的材料阻尼作用,介质中仍存在某种阻尼而使运动衰减,这种阻尼称为辐射阻尼或几何扩散阻尼。
对于面波的类似研究表明,位移随传播距离以的规律衰减,比体波衰减要慢得多。这也就是在震中距远的地点所观测到的面波运动分量要比体波大的原因。
在很多问题当中,比如地震波能量从有限大尺度的断层释放,以及地基基础受到动力机器激励而振动的情况,介质中的辐射阻尼对波动的衰减常比材料阻尼要大。