一种用于直升机旋翼地面声场预测的方法
【技术领域】
[0001] 本发明属于直升机技术领域，特别涉及一种直升机旋翼地面声场高效预测方法， 用于模拟直升机飞行过程中旋翼噪声对地面声场的影响。
【背景技术】
[0002] 随着直升机的广泛应用，直升机噪声大的缺点越来越引起人们的重视，并且成为 其在军民领域进一步应用的主要限制因素。因此建立能够快速评估直升机飞行过程中噪声 对地面影响的预测方法，对于直升机的定型、适航以及使用过程中的降噪都具有重要的支 撑作用。
[0003] 目前地面声场预测只能采用"直接计算法"，即首先进行旋翼的流场计算，得到旋 翼的载荷信息；其次采用相应的噪声计算方法计算桨叶旋转一周内，桨叶表面每个声源点 对地面观测点的噪声影响。对于一般性地评估，桨叶上的声源点和地面观测点的数量均会 达到十万以上，且随着需要评估范围的扩大，地面观测点数量会明显增加。以100000 X 100000(声源点数X观测点数)为例，仅一次评估时间就达3个小时以上。在实际的直升机噪 声评估中飞行状态很多，因而计算时间就非常耗时。
[0004] 此外，现有的地面声场预测方法中并没有考虑直升机真实飞行过程中大气、地面 等环境因素对噪声传播特性的影响。因此，噪声的评估结果会与实际结果相差较大。
【发明内容】
[0005] 针对上述现有技术的不足，为了能够准确、快速地预测直升机真实飞行环境中旋 翼噪声对地面的影响，本发明提出了一种用于直升机旋翼地面声场预测的方法。
[0006] 为实现上述目的，本发明采用以下技术方案：
[0007] -种用于直升机旋翼地面声场预测的方法，包括以下步骤：
[0008] 步骤一，旋翼气动载荷计算，获得桨叶表面气动载荷分布;本发明中对于旋翼气动 载荷及声源的计算预测不限定方法，可以采用简单的叶素理论或者复杂的计算流体力学方 法；
[0009] 步骤二，噪声声辐射球的计算，在获得旋翼气动载荷分布的基础上，建立以位于桨 毂处的声辐射球代替实际桨叶声源，声辐射球上等角度间隔的观测点噪声代表不同方位旋 翼的噪声辐射;选择以旋翼中心为球心、半径为25R的球面，通过基于声学类比法的FW-H方 程计算球心上离散点的声压时间历程以及声压级；
[0010] 步骤三，地面观测域点的离散，根据观测域尺寸及评估精度的需要，对地面观测域 进行离散划分，采用平均划分的方法，或者在需要密切关注的区域进行观测点加密；
[0011] 步骤四，噪声传播模型建立，建立基于频域噪声传播模型，按地面观测点循环，连 线旋翼的桨毂与地面观测点，并计算观测点对应的二级声辐射点，采用二维线性插值方法 插值计算二级声辐射点的噪声;运用时-频变换方法将时域噪声转换为频域噪声，并将频域 信号输入噪声传播模型中，进行噪声传播计算;为了提高计算精度，本发明建立了基于频域 噪声传播模型，其一方面可以考虑由于传播距离增加造成的几何损失，另一方面可以考虑 大气、地面等环境因素对噪声传播特性的影响；
[0012] 步骤五，地面声场的计算，利用步骤二建立的噪声声辐射球以及步骤四建立的噪 声传播模型，直接计算辐射球对地面观测点的噪声;地面声场输出包括:地面观测面噪声的 声压级等高图、各观测点声压-时间历程图及噪声频域图。
[0013] 步骤一中，通过计算获得旋翼桨叶表面的气动载荷，如果计算得到的是展向集中 载荷，则通过弦向平均化获得沿弦向的分布。
[0014] 步骤二中，根据球声源的辐射特征，用二级声辐射球来替代真实的旋翼桨叶声源。
[0015] 步骤四中，运用傅里叶变换地方法将时域噪声转换为频域噪声。
[0016] 步骤四中，噪声传播模型中噪声几何衰减、大气吸声、地面吸声及反射的相关计算 公式为：
[0017]
[0018]
[0019]
[0020] 兵屮，Pf f u，t;，pgr u，t;，pra u，t;力、別衣不μ栄戸按叹总按収到的直达声波、反射 声波和总声波;ri，r2分别表不直达声波和反射声波的传播距离;x，t分别表不观测点接受位 置和接收时刻;α， Π ( Φ 1，ω )表不大气吸声系数和地面反射系数;frl,C)，frl,N分别表不氧气 和氮气的弛豫频率;f表示噪声频率;T表示环境温度；Φ 1表示直达声波入射角；表示地面 的孔隙率。
[0021] 本发明的有益效果是：
[0022] (1)本发明的用于直升机旋翼地面声场预测的方法一方面可在现有直接计算法基 础上，极大降低旋翼噪声对地面影响的计算时间；
[0023] (2)本发明中考虑了大气、地面等环境因素对噪声传播特性影响的噪声传播模型。 相比现有方法，可显著提高预测结果的准确性。
【附图说明】
[0024] 图1是本发明的技术解决方案示意图；
[0025] 图2是二级声辐射模型；
[0026]图3是大气吸声效应对噪声传播特性影响的计算仿真结果；
[0027] 图4是地面吸声效应对噪声传播特性影响的计算仿真结果；
[0028] 图5是地面反射效应对噪声传播特性影响的计算仿真结果。
【具体实施方式】
[0029]下面结合附图和具体实施例对本发明作更进一步的说明。
[0030] 实施例
[0031] 本实施采用AH-1/0LS试验旋翼。
[0032] 本发明的用于直升机旋翼地面声场预测的方法，包括以下几个步骤：
[0033] 第一步，计算状态参数输入，本部分采用的计算状态参数为:拉力系数0.0054，前 进比0.164，总距6.2°，纵横向周期变距分别为-2.4°、1.0°。
[0034]第二步，计算该状态下AH-1旋翼距桨毂中心25R半径处的二级声辐射球，声辐射球 上噪声观测点沿桨叶旋转方向和仰角方向等角度间隔（间隔5°)分布。观测点噪声采用自由 尾迹/FW-H方程进行求解;上述两个过程完成了地面声场预测的前期计算。
[0035] 生成相关的地面观测点分布模型，观测点沿X、Y轴间隔2.5R分布，观测面X、Y轴分 别为120m和60m，如图2所示。
[0036] 第三步，按地面观测点循环，连线桨毂与地面观测