原创力文档- 1、本文档共61页，可阅读全部内容。
- 2、原创力文档（book118）网站文档一经付费（服务费），不意味着购买了该文档的版权，仅供个人/单位学习、研究之用，不得用于商业用途，未经授权，严禁复制、发行、汇编、翻译或者网络传播等，侵权必究。
- 3、本站所有内容均由合作方或网友上传，本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺！文档内容仅供研究参考，付费前请自行鉴别。如您付费，意味着您自己接受本站规则且自行承担风险，本站不退款、不进行额外附加服务；查看《如何避免下载的几个坑》。如果您已付费下载过本站文档，您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等，请点击“版权申诉”（推荐），也可以打举报电话：400-050-0827(电话支持时间：9:00-18:30)。
查看更多
第八章 震源理论
前几章我们叙述了模拟地震波传播的方法，但主要忽略了这些波从那里来和与震源物理性质有关的地震能量是怎样辐射的问题。如果唯一的兴趣在于研究震源区以外区域地球结构的细节，例如速度结构的走时研究，那么就可以忽视这些问题。然而，在许多情况下，分辨地震结构需要有震源特征的某些知识。当然，分辨震源性质是任何真正了解地震的基础。因为震源理论是很复杂的，在这章我们将不拘谨于方程的推导，而是对地震学中实际应用的一些重要结果进行归纳。关于详细的理论论述，请读者查看Aki和Richards的书。
8.1地震的静态位移场
1、集中力所产生的静态位移场
我们想确定各向同性均匀无限弹性介质中由于体力在O点产生的在P点的位移，介质的密度为ρ，拉梅系数为λ，剪切模量为μ。对于离震源较大的距离处，位移u=0，我们定义点源力F为：
(1)
这里f为单位质量上的力，ρf为单位体积的体力，为小体积元。这里介绍三维δ函数:
(2)
运用高斯定理，可以给出：
(3)
其证明如下：
运用此式我们可以用径向坐标(r-1)的空间导数来表示δ函数。
我们知道，在弹性介质内，基本的运动方程为：
(4)
在上述基本的方程中，我们考虑静态问题，即：上式左边的项为零，有：
(5)
考虑在原点有一大小为F的点力：
(6)
这里a为力的方向矢量，我们运用了场论中的公式：。代入运动方程可得：
(7)
我们假定位移解为：
(8)
这里我们基于任何位移场都可以表示一个无散场和无旋场的叠加，我们有： (9)
代入（7）式中可得：
(10)
这样，就得到：
(11)
如果我们令Ap=Apa, As=Asa，我们可以得到如下泊松方程：
(12)
因为,我们可以对上式积分得：
(13)
这些解是方程（11）的势，我们将上式代入（8）式来计算位移，运爱因斯坦求和符号，并设在j方向上的单位力（F=1）在i方向产生的位移为,可得：
(14)
表示成以下简洁形式
(15)
这里，对泊松固体有。按照力的方向建立迪卡尔坐标系（即，力的方向沿着坐标轴的方向），就得到了单位力的位移。上式被称为Somigliana张量。注意，张量是对称的，即，,对于力F，六个独立的分量为：
(16)
我们考虑位于x1方向上的集中力的情况，在极坐标系中有：
（17）
运用雅可比坐标变换，有：
(18)
将位移表示为球坐标系中为：
(19)
将集中力产生的静态位移各个分量辐射花样绘成图1。可见，集中力在无限弹性空间产生的静态位移的径向分量在力的作用方向及其相反方向上达到最大。在与其垂直的方向上为零，这是很容易理解的。切向分量在x1,x2所在平面上为零，在x3轴上达到最大，其切向分量的幅度比径向分量的幅度大1.33倍。切向分量只与角有关。在x1轴上达到最大，其最大幅度与径向分量相同。
%force.m
n=50;
fai = (-n:2:n)/n*pi; %-pi~pi
theta = (0:2:2*n)/n*pi/2; %0~pi
cosfai = cos(fai);
sinfai=sin(fai);
costheta=cos(theta);sintheta=sin(theta);
sintheta = sin
文档评论（0）