![Actran声学仿真分析标准教程](https://wfqqreader-1252317822.image.myqcloud.com/cover/625/52717625/b_52717625.jpg)
1.2.2 声学波动方程及一般解
N-S方程(Navier-Stokes Equation)阐释了质量守恒定律的连续方程以及牛顿第二定律的动量方程。
为了推导出声学波动方程,将对连续方程和动量方程施加下面5个假设条件。
(1)无黏性流体:假设声波传播的流体介质是无黏性的。
将无黏性流体条件公式(1-7)代入动量方程(1-11),可以得到欧拉方程:
![](https://epubservercos.yuewen.com/E1E8DE/31310243904416906/epubprivate/OEBPS/Images/tx00715.jpg?sign=1739567582-2s27ccJonO7m3MR6I8qD5o792dcgDsEk-0-36b00328ee341ecf4994f8f0a634210f)
(1-14)
(2)小扰动量:将压力场和密度场分解为两部分,即在时间和空间上恒定的参考值()和比参考值小很多的波动值(
)。
(3)线性本构关系:利用状态方程,把压力与密度联系起来,如图1-5所示。
在理想气体的某压力与密度组合状态附近,假设线性化的关系为:
![](https://epubservercos.yuewen.com/E1E8DE/31310243904416906/epubprivate/OEBPS/Images/tx00739.jpg?sign=1739567582-tEm9Z1y4w3mW6wrR8KfoCHy5rp38H63z-0-d357af9bd78c55b2427cfa55cc7114fa)
(1-15)
其中,数值c2被定义为:
![](https://epubservercos.yuewen.com/E1E8DE/31310243904416906/epubprivate/OEBPS/Images/tx00750.jpg?sign=1739567582-xf4JSIkEAGJYX33ORSuZoyIVVZOCxtF1-0-34d0f38e95aede97e48cf207aa24fd73)
(1-16)
![](https://epubservercos.yuewen.com/E1E8DE/31310243904416906/epubprivate/OEBPS/Images/tx00752.jpg?sign=1739567582-DSamOcXI76CJE3k4pyGGAQcAva7unTqk-0-1db21cb2a221b4d3fcc81937912a203c)
图1-5 理想气体的压力-密度关系以及在参考状态附近的线性化
下文中会得出,数值c对应压力波在流体中传播的速度。
(4)流体宏观静止条件:流体的运动仅限于由压力波动引起的微小扰动。换句话说,流体不存在平均流速,而且声学速度(流体粒子振动速度)比声波传播速度小很多。
(5)对方程的线性化:将小扰动量假设分别代入连续方程和欧拉方程中,对方程进行线性化,并忽略质量源。方程组可以写为如下形式:
![](https://epubservercos.yuewen.com/E1E8DE/31310243904416906/epubprivate/OEBPS/Images/tx00753.jpg?sign=1739567582-uRuhbByoQQpzcXA5KihJf4XpPCj9OLw0-0-d5184b14bc0d5c4da959669d0c61aa47)
(1-17)
![](https://epubservercos.yuewen.com/E1E8DE/31310243904416906/epubprivate/OEBPS/Images/tx00754.jpg?sign=1739567582-Utnp9Mj7JzUmXRobYoBf26LsYvCiSxdG-0-5a7daae4456e9d5e0f716661d6b2a31b)
(1-18)
上述第一个方程称为波动方程或达朗贝尔方程。第二个方程显示了在不存在分布源项的情况下,流体的加速度与声压的梯度成正比。
通过对N-S方程的线性化,可得到用来描述部分声学波动的达朗贝尔方程。事实上,大部分的声学现象是遵循线性规律的。线性规律意味着物理现象的效果与产生此现象的原因是成一定比例的,而且若干现象之间遵循线性叠加的原则。
波动方程的一维形式为:
![](https://epubservercos.yuewen.com/E1E8DE/31310243904416906/epubprivate/OEBPS/Images/tx00755.jpg?sign=1739567582-D2Kt3R5hI7eF6j0d2G8qiqFp1g9p38Ff-0-146e5acbbe512ff4598e4b4a02029dcb)
(1-19)
可以证明所有写成以下形式的函数均是上述方程的解:
![](https://epubservercos.yuewen.com/E1E8DE/31310243904416906/epubprivate/OEBPS/Images/tx00764.jpg?sign=1739567582-8Su9YcHauCZp3ey264DH0QjIMahKxfcf-0-c00075878f8b12d48abf02dd33429c30)
(1-20)
一维波动方程的一般解为两项之和(见图1-6)。其中,第一项表示从左向右传播的声压分布,其在传播过程中信号没有任何变形。由此可以看出这一项在x点或在x+Δ点取值相同,但两点之间存在Δ/c的时间差:
![](https://epubservercos.yuewen.com/E1E8DE/31310243904416906/epubprivate/OEBPS/Images/tx00777.jpg?sign=1739567582-tjb9JefFaAUzg3YzuRAss15kAI9rmn1L-0-f16943cb84aee04c381586d3a24b233e)
(1-21)
一般解的第二项表示压力分布从右向左地传播:
![](https://epubservercos.yuewen.com/E1E8DE/31310243904416906/epubprivate/OEBPS/Images/tx00792.jpg?sign=1739567582-tZl1gghefJ8MaTxV9Zkw6Ozjc1RHs8an-0-7ecdfad4ebc55c82765c5e78668df61f)
(1-22)
![](https://epubservercos.yuewen.com/E1E8DE/31310243904416906/epubprivate/OEBPS/Images/tx00799.jpg?sign=1739567582-arH5iB5mV0ctbtGWxdKUbr8cg7VFH6Xi-0-bc883b18a1e36c35c44722bfd6f4e1bb)
图1-6 一维波动方程的一般解是两项之和:以速度c向相反方向传播的两个波动量