地震波动方程 (一).pdfVIP

第三章地震波动方程

现在，我们用前一章提出的应力和应变理论来建立和解在匀称全

空间里弹性波传的地震波动方程。这章涉与矢量运算和复数，附录

2对一些数学问题进行了复习。

3.1运动方程(EquationofMotion)

前一章考虑了在静力平衡和不随时间改变状况下的应力、应变和

位移场。然而，因为地震波动是速度和加速度随时间改变的现象，因

此，我们必需考虑动力学效应，为此，我们把牛顿定律(”=〃也)用

于连续介质。

3.1.1一维空间之振动方程式

质点面上由于应力差的存在而使质点产生振动。如图1-3所示,

考虑一薄棒向x轴延长，其位移量为u：

dx、ds

u

Fig3-1

则其作用力为“应力”X”其所在的质点面积”，所以其两边的作用

力差为

ds(cr(x+dx)-cr(x))=—dxds

dx

惯量(inertia)为

,.d2u

pdxds旅

所以得出

P—Z-=——

ordx

(3-1)

其中P为密度(density),。为应力(stress)o

dx

3-1式表示，物体因介质中的应力梯度(stressgradient)而

得到加速度。假如P与E为常数，则3-1式可写为

22

du1du

222

dxcdt

(3-2)

运用分别变量法求解(3-2)式,设u=F(x)T(t),(3-2)式可以变为

v*r_VT*

则可得：Toc/rxocj

考虑欧拉公式：*=cos(cot)+isin(rtX),=cos(cw)-isin(由)

/—(x+rr)i-(x-ci)(x-r/)/、

c

u=Ae+Be，+Ce，+De，(3-3)

其中A,B,C,D为依据初始条件和边界条件确定的常数。

考虑到①可正可负，方程式的解具有u=f(x-cf)+g(x+cf)的形

式，其中f与g为波的函数，以c的波行速度向+x与-x方向传递。

我们可以采纳如下程序模拟地震波的传。平面波在匀称介质里

沿工方向传，剪切波的齐次微分方程可表达为：

这里〃是位移。对100公里的波长和假定夕=4公里/秒的状况，我们写

出用有限差分法解这方程的计算机程序。用长度间距公=1公里，时间

间距力=0.1秒。假定在〃(50公里)震源时间函数的形式为：

〃5o(/)=sin2(/r〃5)0V/5秒

用“(0公里)的应力自由边界条件和〃(100公里)的固定边界条件。

用有限差分图解来近似二次导数：

22

dx~dx

以4秒的间隔画出『33秒的图。

M=moviein(101);

dx=l;dt=O.1;tlen=3;beta=4;%初始化变量，tlen为震源持续时间，

beta为波传