地震波动方程.docxVIP

第三章地震波动方程 现在，我们用前一章提出的应力和应变理论来建立和解在均匀全空间里弹性 波传播的地震波动方程。这章涉及欠量运算和复数，附录2对一些数学问题进行 了复习。 3.1 运动方程（Equation of Motion ） 前一章考虑了在静力平衡和不随时间变化情况下的应力、 应变和位移场。然 而，因为地震波动是速度和加速度随时间变化的现象， 因此，我们必须考虑动力 学效应，为此，我们把牛顿定律（F ma）用丁连续介质。 3.1.1 一维空间之振动方程式 质点面上由于应力差的存在而使质点产生振动。 如图1-3所示，考虑一薄棒向x轴延 伸，其位移量为u： 则其作用力为“应力” X “其所在的质点面积”，所以其两边的作用力差为ds x dx x惯量(inertia )为 则其作用力为“应力” X “其所在的质点面积”，所以其两边的作用力差为 ds x dx x 惯量(inertia )为 2 udxds —2 t2 所以得出 2 _u t2 x — dxds x (3-1) 其中p为密度（density ） , r为应力（ stress) = E 3-1式表示，物体因介质中的应力梯度（ stress gradient）而得到加速度。如果p与 E 为常数，则3-1式可写为 2 2u 1 2 2 u 1 u ~2 ~2 2 x c t (3-2) 其中c E 运用分离变量法求解(3-2)式，设u=F(x)T(t),(3-2)式可以变为 运用分离变量法求解( 1 X T —XT c TOC \o 1-5 \h \z 方 c2X T 2 设 X T ― 「t 则可得：T e ,X e c 考虑欧拉公式：e1 t cos( t) i sin( t), e 1 t cos( t) i sin( t) i — x ct i_ x ct i_ x ct i — x ct u Ae c Be c Ce c De c (3-3) 其中A,B,C,D为根据初始条件和边界条件确定的常数。 考虑到 可正可负，方程式的解具有 u f x ct g x ct的形式，其中f及g 为波的函数，以c的波行速度向+x与-x方向传递。 我们可以采用如下程序模拟地震波的传播。剪切波的齐次微分方程可表达为：2 我们可以采用如下程序模拟地震波的传播。 剪切波的齐次微分方程可表达为： 2u下 这里u是位移。对100公里的波长和假定 差分法解这方程的计算机程序。用长度间距 定在u (50公里)震源时间函数的形式为: u50 t sin2 t 5 用u (0公里)的应力自由边界条件和u 差分图解来近似二次导数： 2 u ui 1 2 x 以4秒的间隔画出1-33秒的图。 M = moviein(101); dx=1;dt=0.1;tlen=3;beta=4; % 初始化变量,tlen u1=zeros(101,1);u2=u1;u3=u1; 平面波在均匀介质里沿x方向传播， 2 2 u 2 x 4公里/秒的情况，我们写出用有限 dx 1公里，时间间距dt 0.1秒。假 0 t 5 秒 (100公里)的固定边界条件。用有限 2ui ui 1 dx2 为震源持续时间，beta为波传播的速度 %u伪前一个时刻的各点的位移， u2为当前时刻的位移，u3为下一个时刻的位移值，开始均假定 为零 t=0; jj=0； while (t=33) % 模拟的最长时间为33秒 for ii=2:100 rhs=betaA2*(u2(ii+1)-2*u2(ii)+u2(ii-1))/dxA2; % 方程的解 u3(ii)=dtA2*rhs+2*u2(ii)-u1(ii); % 对时间求导数 end %左边为自由边界条件，右边为固定边界条件 u3(1)=u3(2); % 左边为自由边界条件 u3(101)=0.0; % 右边为固定边界条件 %左右两边为自由边界条件 % u3(1)=u3⑵；% 左边为自由边界条件 % u3(101)=u3(100); % 右边为自由边界条件 %左右两边为固定边界条件 % u3(1)=0.0; % 左边为固定边界条件 % u3(101)=0.0; % 右边为固定边界条件 if(t=tlen) u3(51)=(sin(pi*t/tlen)).A2; % 地震震源时间函数 end for ii=1:101 u1(ii)=u2(ii);u2(ii)=u3(ii); % 时刻的更新 end plot(u2); % 绘制目前的波形图 ylim([-1.2 1.2]); M(:,jj+1) = getframe; % 获得当前的图像 t=t+dt; % 时间延长 end movie(M) %演示波形传播 3.1.2三维空间之振动方程式 推导三维空间之振动方程式