1、 空气中声音传播方程的推导数学与应用数学 2011 级作者： 王锴丰，王保山，王冬羽1.问题提出：由物理性质可知声音的本质是介质的机械震动，由波动方程与震动现象的联系，物理属性的时空分布可以由波动方程导出，弦的波动方程也是声音发生的一种情况之一，现讨论声音在空气中的传播，由声学基础可知声场的的物理特征可由密度 ，压强 ，速度 刻画，根据模型所tzyx,tzyxP,tzyxv,满足的基本假设和遵循的物理规律，从而导出声波方程，并求解出普通声源产生声波在均匀空间中的传播的状态方程。2.理想假设：1.空气在压强为101.325kPa、温度为20的条件下，空气动力粘度和运动粘度为：空气 ， sPa6-
2、smv26108.4可以近似为理想流体，所谓理想流体介质，就是介质在运动过程中没有能量的损耗，即介质是无粘的不考虑切向力。2.还必须假设介质是连续、静态和均匀的流体因此媒质中质点速度 为零，静态压强 ，静态密度 v0P都是常数。3.由理想气体的状态方程 温度变化会影响压强，未排除这种干扰，将声音在空气中传播nRTPV的过程视为等熵绝热过程，介质之间不会产生热交换。4.讨论的都是小振幅声波即各类声学变量都是一级微量， （ 4 ）媒质中传播的是小振幅声波，各声学参量都是一级微量，声压 甚小于媒质中静态压强 ；质点速度 甚小于声速 ；质点位移 甚p0Pvc小于声波波长 ；媒质密度增量甚小于静态密度；
3、或密度的相对增量 小于 1。 0tPS5.声源相对于传播空间来说很小，可视为一个半径非常小的球体，故由对称性声音在空间中的各点的特征（声压，质点速度，以及密度）相对于声源来说球对称，.只与半径 r 有关。3.模型建立即方程推导理想介质中声波传播的基本规律可以通过三个方程表示，即：连续性方程、运动方程和状态方程。下面给出所用到的物理参数， ct vPpP， 声 速 ：时 间 ： 质 点 速 度 ：气 体 密 度 ：总 压 强 ：变 化 声 压 ：静 态 压 强 ： ,0连续性方程：声场中任意一点 ，以 点为中心选取一边长分别为 的体积元，则体积元的zyxB,Bdzyx、体积为 。假设某一瞬时，介
4、质质点流过 点的速度向量为 ， 点的密dVBtv,B度为 ，则单位时间内从 各方向流入体积元的质量为，tzyx,zyx、(1.1)dxyzzvyxvdm式中 为 点的速度向量 在三个坐标轴上的投影。流入体积元的质量必然引zyxv、 Btzv,起体积元内密度的增加，单位时间内体积元内介质密度的增量为 ，则，t(1.2)dxyztm根据质量守恒定律，得到，(1.3)zvyxvt 即为：(1)0tv状态方程：在声波作用下介质的状态发生了变化。由假设2，根据热力学关系，一定质量的气体介质的压强是密度和熵的函数，即 ，这里 表示熵。由于压力和密度变化很小，因此由泰勒级数展开得到，sfP,(2.1)高 阶
5、 小 量dfPps00式中下标 代表等熵绝热过程。系数 ， 就是气体介质中小振幅声波的传播速度，故0s 20cfs理想介质的 状态方程为，(2)dcp2运动方程:假设微小质团中心坐标为 ，体积为 ，介质原处于静止状态 ，当声zyxB,dxyzV0v波通过时，由假设一，不考虑切向力，只考虑压力，质团在各个方向上的受力都不均衡，假设压力分布为 ，则作用在 和 面上的总压力分别为，tzyxP,2dx(3.1)dyzxPydzPzyxdx 2,2(3.2)yzxyzzyxdx ,2因而沿 轴正方向的合力为，xdxyzPdFzx,(3.3）同样可以得到质团沿 轴和 轴正方向的合力，介质受到的总的合力为，
6、yz(3.4)dVzyxPF,由于静压强 为常数，因此压强的微小变化也就是声压的微小变化，即 ，所以上式可以改0P pP写为，(3.5)dVpF根据牛顿运动定律，得到，(3.6)dtvVF由以上两式可以得到欧拉方程，(3.7)pdtv式中 是质点 的加速度，它包括本地加速度和迁移加速度两部分，dtv，0B(3.8)vtdv由假设4振动速度远小于声传播的速度，所以 项可以忽略，结合以上两式，得到小振幅声场中的 运动方程为，(3)ptv0最后，将（1）式对 t 求偏导，再用（2） （3）式中的变量替换就可以获得理想流体介质中小振幅波传播的 声波方程如下，（4 ）221tpcp通过以上推导终于得出声
7、压满足波动方程，这是该问题的核心。4.定解条件提出由假设5相对于较大的空间， 声波以球面波形式传播时，可看作半径非常小的球体，半径为 ，声r0源，声压 只与球面坐标 有关，而与角度无关，因此球面波的波动方程可以简化为，pr(4.1)221tpcrp初始条件为： ， ),(0,0, zyxzyxzzyxp）（）（其中 时022),(),( r当其中在即问题为柯西问题求解5.问题求解由三维问题的柯西问题的解法用到球平均法，这里不多赘述，用泊松公式可以得到原问题的解：（5）),(),(),( tzyxtMczyxttzyxp其中 MctMctSSdctzyxt2241),(),(故压强场为：（6）0
8、),(PtzyxpP由公式（2）可知密度场为：（7 ）2cP由牛顿运动定理可知质点振动速度场为：（8）tPdtzyxv01),(6.结论：本小组通过流体力学知识，推导证明出了声波在空气中传播满足波动方程，并且对声波在空间中一点扩散模型的解析解进行了求解，在求解过程中由数学物理方程可知，三维问题的柯西问题，用到泊松积分公式，求出了空间中的压强场，密度场，与速度场。参考文献：数学物理方程 ，谷超豪，李大潜，陈恕行，郑宋穆，谭永基（高等教育出版社)2002.7大学物理学 （上） ，赵近芳，颜晓红 （北京邮电大学出版社）2008.6声学基础第三版，杜功焕，朱哲民 ，龚秀芬（南京大学出版社）2012.5
