波动声学基础 - 简正波声场计算及波形预报仿真程序说明文档.pdf

简正波声场计算及波形预报 简介 本算例根据简正波理论计算水下声场（同时仿真了绝对软海底和绝对硬海底 两种情况），并对时域波形进行预报，本算例对应水声学原理第三章 3.2 节的部 分内容。 1.1 基本理论 根据柱坐标系下的波动方程： 2 1   p   p 2 r  k p 4A r r   2 0  s  r r  r  z 利用分离变数法，并代入边界条件，求得理想边界条件下声场声压的解析解 为： 2 p (r ,z ) j  sin(k z )sin(k z )H (2) ( r ) zn zn 0 0 n H n 1  其中绝对硬海底条件下简正波的本征值为：kzn (n  ) n 1,2,3, 2 H 绝对软海底条件下本征值为：kzn n  n 1,2,3 H 其中 k 2 k 2 ， 为波导深度。 H n 0 zn 上述得到的声场实际上是稳态场，也即声场声压的频域解，为此将信道考虑 成一个二端网络，计算某一空间位置不同频率的声场结果，得到的结果即为信道 的传输函数，再利用卷积定理与傅立叶变换（逆变换）即可预报声压的时域波形。 1.2 数值仿真 仿真参数：发射信号频率25Hz ；水中声速1500m/s；波导深度200m ；声源 深度100m；声场预报水平最大距离10000m。 1.2.1 简正波声场计算 仿真结果： 图1 10km 范围内绝对硬海底情况下声场的空间分布 绝对硬海底 30m深度处传播损失曲线 30 40 50 60 m / B d 70 80 90 100 0 1000 2000 3000 4000 5000 6000 7000 8000 9000 10000 Range/m 图2 绝对硬海底情况下深度30m 处的传播衰减 图3 10km 范围内绝对软海底情况下声场的空间分布 绝对软海底 30m深度处传播损失曲线 30 40 50