弹性波方程和Christoffel方程.docVIP

弹性波方程和Christoffel方程 Sdhizhj 介质在直角坐标系中的运动方程 设想连续弹性介质中小立方体Δv=Δx1Δx2Δx3，考察x1=0面和x1= Δx1面上的沿x1方向的应力，则x1=0面上的应力为T11，x1= Δx1面上的应力为； 作用在与x2垂直的两个面上的沿x1方向的应力分别为T12和； 作用在与x3垂直的两个面上的沿x1方向的应力分别为T13和； 则作用在小立方体6个面上沿x1方向的应力和为： 根据牛顿第二定律，有： 式中，ρ为介质的体密度，u1为质点沿x1方向上的位移，化简得： (1a) 同理有： (1b) (1c) 式中，ui为质点沿xi方向上的位移，上面几式表示介质在直角坐标系中的运动方程，可以用下式概括表示： （i,j=1,2,3） (2) 引用爱因斯坦求和表示，为： （i,j=1,2,3） (3) 上式中约定，当物理量脚标出现重复时就自动求和。 非压电弹性介质中的波动方程 根据胡克定律，弹性介质中，应力与应变有如下关系： T=cS Tij = cijkl . Skl (4) 对于各向异性介质，应力T与应变S为张量，有6个独立分量，而弹性刚度系数则是四阶张量，有36个独立分量。 而应变Skl与位移ui之间有关系 （k,l=1,2,3） (5) 代入胡克定律，考虑对k,l求和的时候k,l将会分别遍历所有坐标，因此 所以有 （i,j,k,l=1,2,3） (6) 上式就是非压电弹性介质的波动方程，求解这一方程就能得到非压电弹性介质中弹性波的基本特征。 Christoffel方程 考虑沿任意方向传播的均匀平面波。设传播方向的方向余弦为(l1,l2,l3)，ω＝2πf，为圆频率；为波数矢量，，li为在I方向的方向余弦，q=2π/λ；速度。则载波介质中的质点位移矢量ui可以表示为 （i=1,2,3） (7) 其中，Ai为质点位移得振幅。 代入波动方程，得 （i,j,k,l=1,2,3） 上式可以变换为 （i,j,k,l=1,2,3） (8a) 上式实际上就是Christoffel方程。 其中，， Christoffel方程的矩阵形式为： （i,j,k,l=1,2,3） (8b) (9a) 化为缩减坐标形式为： (9b)