声学中波动方程的建立

1、田佳星海洋技术12020041049今天我介绍一下声学中波动方程的建立。我们首先介绍一下声学的基本概念。声波是机械振动状态在介质中的传播。存在声波的空间称为声场。理论上描述声场需要引入一些物理量：声压、位移、振速、密度压缩量和相位等。通常采用上述各物理量的时空分布函数描述声场。下面对这些物理量作简要介绍。基本概念1）声压（标量）声波为压缩波。描述“压缩”过程的一个物理量是压强。然而，声波是声扰动（如振动源）引起介质中的压强发生变化的部分。因此，我们引入声压的概念：声压p为介质压强的变化量：p=PP（2-1）0其中，P是压强，P是介质中的静态压强。0声压是描述波动的物理量。为使用方便，还由声压引

2、入了瞬时声压p、峰值声压p和有效声压p。0e声场中某瞬时的声压称为瞬时声压。一定时间间隔内的最大瞬时声压称为峰值声压。瞬时声压在一定时间间隔内的均方根值称为有效声压，即对简谐声波，p、p和p相互之间的关系和电压可作相同类比，即0ep二pexpjt0p=p2。e0一般仪器仪表测得是有效声压。2）位移和振速（矢量）质点位移是指介质质点离开其平衡位置的距离。质点振速是介质质点瞬时振lmrbsmrbsmrbsmrb.(2-3)动的速度。两者均是有大小和方向的量，即矢量，相互关系为u=ddt对简谐振动，位移和振速都满足如下关系：(2-4a)(2-4b)=expjot，0u=expjt，其中，和u分别为位

3、移幅值和振速幅值。00需要注意的是区分质点振速和声传播速度。声传播速度是指振动状态在介质中传播的速度，而质点振速是指在给定时间和给定空间位置的某一质点的振动速度。密度和压缩量(2-5)密度的变化也是描述声波的一个物理量。这里引入压缩量的概念：其中，,密度，,为静态密度，,=pp为密度改变量。010压缩量s的含义为介质密度的相对变化量。相位为描写简谐振动而引入的物理量。它描述质点简谐振动的状态。质点振动的一个周期对应着相位0-2n。相位和质点振动状态有对应的关系。声波是振动状态在介质中的传播，而相位描述的是质点简谐振动的状态。由此可见相位在声场描述中的重要性。以上物理量并不是独立的，如根据位移由

4、(2-3)式可以求出振速。实际应用时可根据需要选择使用哪些物理量来描述，如对简谐声波，只需要位移幅值和相位就可导出振速、加速度等基本物理量；更进一步，如果已知介质条件，只要知道位移幅值和相位的初值，就可计算声场的时空分布函数了。理想流体介质中的小振幅波本节先建立描述声波的基本方程-波动方程，并讨论波动方程的线性特性；然后分别介绍波动方程在几种简单介质条件下的解-行波解、平面波解、球面波解和柱面波解，并对各种解中相关的物理量，如声场中的能量、介质特性阻抗和声阻抗率、相速度和群速度等概念，进行讨论，并重点分析在水声物理中应用较多的平面波在两种不同均匀介质界面上的反射和折射现象。一、波动方程2.1建

5、立波动方程为更清楚地了解声波的物理本质，我们先对介质条件和声波做出一定的限制，而得到形式简洁的波动方程，并通过它认识声波的物理本质。在后续的学习和研究过程中，将不断引入更为复杂的介质条件和放宽对声波的限制，再进行研究。这也是物理中研究常用的方法之一。假设条件：介质静止、均匀、连续的；介质是理想流体介质，即忽略粘滞性和热传导；声波是小振幅波。连续性方程0理论推导见教材。连续性方程即质量守恒定律：介质流入体元的净质量等于密度变化引起的体元内质量的增加。-V(pu)at(2-6)根据假设条件有：-V(pu)ato(2-7)事实上，当介质本身有流动时，u中含有介质流动速度的影响，相关理论可参阅朗道著连

6、续介质力学。考虑到假设介质是静止的，(2-6)式和(2-7)式中没有考虑介质流动速度的影响。状态方程_在理想流体介质声传播过程中，还没有来得及进行热交换，声波传播(介质的压缩和膨胀)的力学过程已经完成，这一过程近似为绝热过程，即无热传导。绝热过程中其中，定义dPI-d,c2dp丿s(2-8)0PP0(2-10)(2-9)P为压强，p为密度，下标s表示绝热过程。本节后面讨论波动方程的解时，可知c为声波在介质中的传播速度。声速是介质固有的特性，是由介质的物理参数所确定的。下面由理想流体介质的绝热状态方程导出声速c和介质参数的关系。一定质量理想流体的绝热状态方程为00其中，P和p分别为平衡态下流体的

7、压强和密度，Y为流体定压比热和定容比00热之比。(2-10)式表明，在绝热状态下，流体压强只是密度的函数。对(2-10)式微分得(2-11)dPpr-1,Yd,P0由(2-8)、(2-10)和(2-11)式可推得Pc，一利用小振幅近似可以给出声速的近似表式。由(2-12a)式可知，c只是密度P的函数。将c2在p附近展开得0(2-12a)更进步，dPdP丿s,0(d2P、(dP2丿s,0忽略小量(p-p)之后的项近似得0s,0(P-P。)+式中c是在平衡态时的声速。因此，0(2-12b)c，P引入平衡态下的绝热压缩系数(单位压强变化引起的体积相对变化)卩，-必，sVdP利用质量不变表式pV，PV

8、0和小振幅近似，1s0c二c二0V以上几种声速表式，可根据使用方便选用。另外，以后除特殊说明,用c表示声速，而不用c0(2-12c)以后只3)运动方程运动方程实质上连续流体介质中的牛顿第二定律。理论推导见教材。duP忑，P，-p(2-13)式中，理本地加速度、(u)迁移加速度。对于小振幅波，迁移加速度项()为二阶小量，略去后再根据假设条柿:P遁，-pSt(2-14)4)波动方程0对(2-14)求散度(两端同时左乘),将(2-7)式两端对时间求偏导数，注意到时、空导数次序对调不变，消去振速项，再利用(2-8)式可得波动方程(2-15)其中，在直角坐标系中：2,竺+竺+竺dx2dy2dz2在球坐标

9、系中：2=r2dr丫2d、+1dIsin4!kdrJr2sindkd丿d1d2d2r1+。drJr2d2dz2在柱坐标系中：2,-rdr1d2+r2sin2d200引入速度势物理量(2-16)它表示单位介质质量具有的声扰动冲量。由(2-16)式可知，声压、质点振速与速度势的关系是：(2-17)dp=P0dt(2-18)将(2-17)和(2-18)式分别代入连续性方程、状态方程和运动方程，作和(2-15)式类似的推导，可得关于速度势的波动方程：2P，竺。(2-19)c2dt2(2-15)和(2-19)两式具有相同的形式，但却是不同物理量满足的声纳方程。实际使用时，两者边界条件和初始条件的表述不同，可根据解算的方便选择不同的物理量来描述声场，当然也就需要选择相应的波动方程了。速度式和声压、振速等物理量间的转换关系由(2-16)、(2-