弹性波波动方程及其解课件.ppt

《弹性波波动方程及其解课件.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、第六章弹性波波动方程及其解§6.1线性弹性动力学的基本方程基本方程运动微分方程几何方程本构方程分析已知方程数——15个；未知数——15个；边界条件和初始条件边界条件给定了弹性体在其边界面上所满足的条件。边界条件分类位移边界条件：当S=SU时应力边界条件：当S=St时混合边界条件：当S=SU+St时在St上在SU上初始条件初始条件给定了弹性体在时刻t=0时的位移和速度，称为初始条件。在V+S的弹性体上有：定解条件边界条件+初始条件弹性波动力学的求解路线弹性波动力学问题的表述：弹性体的形状、大小以及其物理性质(即密度和弹性系数)；弹性体所受外来作用的体力及表

2、面力；弹性体所受的约束性；弹性体各点的初始位移及初始速度，求:弹性体内的位移场、应变场及应力场。按应力求解按位移求解2三维三分量波动方程各向同性介质三维三分量波动方程纳维方程是线性弹性假设条件下得到的各向同性弹性体中的弹性波最基本方程。指标表示的纳维方程向量表示的纳维方程线性弹性体的整个理论就是在定解条件下解纳维方程。三维三分量的含义位移场ui是空间坐标x1、x2、x3的函数，因此，是三维立体空间的；位移场ui有u1、u2、u3三个分量，因此，是三分量的；体力为零时的纳维方程通常在研究地震波的传播时，认为地壳所受的体力为零，即：f=0。此时，纳维方程可表

3、示为：3.弹性流体介质中的波动方程弹性流体介质中的基本方程几何方程运动微分方程(不计体力)本构方程？弹性流体中的本构方程黏滞力：在实际流体中两层流体间的相互滑动时，流体间有相互作用的阻力，称其为黏滞力或内摩擦力。理想流体介质：可以将黏滞力忽略的流体称为理想流体。在理想流体中只存在胀缩力，而不存在剪切力。右图为理想弹性体内的某个体元，其只受外法线方向相反的正应力，而无剪应力。即：上式中P为压力，当体元为单位体元时，P可视为压强。PP弹性流体中的波动方程将上式代入运动微分方程，得上式两边取散度将压力与应变关系代入上式§6.2无旋波和无散波斯托克斯-亥姆霍兹矢

4、量定理任何一个足够平滑的矢量场都可以分解成无旋的部分和无散的部分，这称为斯托克斯-亥姆霍兹矢量定理。结论无旋位移场的散度对应弹性休的胀缩应变场；无散位移场的旋度对应弹性体的切应变情形；在非稳定条件下，这两种场分别以波的形式运动着，故分别叫做无旋波和无散波，也称之为胀缩波与等体积波。无旋位移场波动方程试证对于任意一阶张量都成立！结论：在均匀各向同性弹性体内，膨胀扰动以速度VP向外传播。这种膨胀波称为纵波或P波。其传播速度为VP。如果只研究纵波的传播问题，可以不考虑振源的影响，此时设外力为”0”，则方程可简化为：无散位移场的波动方程结论：在均匀各向同性弹性体

5、内，切变扰动以速度VS向外传播。这种切变波称为横波或S波。其传播速度为VS。如果只研究横波的传播问题，可以不考虑振源的影响，此时设外力为”0”，则方程可简化为：体应变表示的纵波方程此时只有胀缩波，波有旋转波！为什么？如果只研究纵波的传播问题，可以不考虑振源的影响，此时设外力为”0”，则方程可简化为：上式表示波场是以速度VP向外传播的无旋场。转动矢量表示的横波方程弹性体发生剪切形变时，由于转动很小，由矢量分析可知，定义转动矢量：如果只研究横波的传播问题，可以不考虑振源的影响，此时设外力为”0”，则方程可简化为：§6.3标量势与矢量势拉梅势(Lame＇)根据

6、斯托克斯—亥姆霍兹矢量分解定理，位移矢量场可以分解成无旋场与无散场两个部分。如果位移矢量表示的纳维方程二次可微，则存在一个标量势函数f和一个矢量势函数y，此时位移场可表示为：同理，对于体力也存在一个标量势函数F和一个矢量势函数A，此时体力可表示为：表示无旋场？表示无散场？以势函数表示的波动方程将势函数代入上式，则：因为在空间区域V及任意时间区域中有上式成立，则有：如果对标量势表示的波动方程两边取梯度如果对矢量势表示的波动方程两边取旋度如果F=A=0，则因此，用标量势函数和矢量势函数表示纵波和横波的传播是合理的！弹性波波动方程的一般表达式由以上的讨论可知，

7、位移矢量u、位移势函数f和y、体积应变q及转动角矢量W都具有相同的方程形式，为了方便起见，可以用统一的方程表示。当有体力时，可表示为达朗贝尔方程当只考虑波的传播，不考虑体力时此时的纵波波动方程为横波波动方程§6.4三维三分量波动方程的退化处理以上讲到的波动方程都是三维三分量的，与当前地震勘探的实际情况有此差别。由于条件的限制，目前在进行地震勘探时所有的检波器都放在地面上，而且也主要是接收垂直地面的垂向分量的地震波。因此，在处理很多问题时我们将三维三分量波动方程进行退化处理。所谓退化处理就是人为地降低波动方程的维数或分量。二维单垂向分量波动方程(2D1VC

8、)二维单垂向分量是目前常规二维地震勘探数据采集的观测方式，即在地表直测线上采集地