Research on microphone phase array design based on surrogate model
-
摘要: 气动声学风洞试验过程中,针对目标声源特性进行麦克风相位阵列改进时,为保证风洞试验效率,必须在较短的时间内完成阵列优化设计工作。为满足这一试验需求,最大限度提升麦克风相位阵列设计效率,引入了基于Kriging代理模型的优化设计方法。采用基于点扩散函数的计算程序进行阵列性能分析,获取阵列最大旁瓣水平和分辨率。通过对样本点响应值进行计算,建立Kriging代理模型,进而以计算速度极高的Kriging代理模型作为阵列性能分析方法开展优化搜索,避免了大量调用阵列性能计算程序导致计算耗时过高的问题,显著提升了麦克风相位阵列设计效率。该阵列设计方法能够有效改善麦克风阵列的测量性能,满足声学风洞试验的特殊应用需求。Abstract: In aero-acoustic wind tunnel tests, the improvement of the microphone phase array according to the noise source features must be conducted as fast as possible to ensure the required efficiency of the tests. An optimization method based on the Kriging surrogate model is proposed in this study. The maximum sidelobe level and the resolution of the microphone phase array are obtained by analyzing the results computed from the point spread function. A Kringing surrogate model is constructed based on the response values computed on sample points. Then the surrogate model was utilized to fast assess the maximum sidelobe level and the resolution of the microphone phase array, which avoided the high computational cost caused by computing the point spread function. This method can improve the performance of the microphone phase array effectively for aero-acoustic wind tunnel tests.
-
Key words:
- microphone phase array /
- optimization /
- surrogate model /
- resolution /
- maximum sidelobe level
-
表 1 Kriging模型预测误差
Table 1. Prediction error of Kriging model
最大旁瓣水平/dB 分辨率/m 阵列性能分析结果 -14.696 0.1467 代理模型预测结果 -14.680 0.1460 相对误差 0.07% 0.45% 表 2 分辨率优化阵列设计变量
Table 2. Design variables of resolution optimized array
设计变量 n k a q 变量范围 8~10 0.60~1.40 0.050~0.200 0.600~1.400 初始阵列变量值 9 1.00 0.150 1.000 优化阵列变量值 9 0.89 0.147 0.716 表 3 初始阵列与分辨率优化麦克风阵列性能
Table 3. Performance of initial array and resolution optimized microphone array
初始阵列 优化阵列 频率/
kHz最大旁瓣水平/dB 分辨率
/m最大旁瓣水平/dB 分辨率
/m1 -17.412 0.388 -14.237 0.341 2 -16.882 0.187 -14.330 0.164 3 -15.535 0.123 -11.868 0.108 4 -15.252 0.089 -11.934 0.080 5 -14.771 0.070 -11.992 0.061 6 -14.188 0.057 -10.235 0.048 表 4 初始阵列与最大旁瓣水平优化麦克风阵列性能
Table 4. Performance of initial array and MSL optimized microphone array
初始阵列 优化阵列 频率
/kHz最大旁瓣水平/dB 分辨率
/m最大旁瓣水平/dB 分辨率
/m1 -17.412 0.388 -20.537 0.412 2 -16.882 0.187 -21.922 0.198 3 -15.535 0.123 -19.346 0.129 4 -15.252 0.089 -16.105 0.094 5 -14.771 0.070 -15.698 0.074 6 -14.188 0.057 -14.009 0.061 表 5 初始阵列与优化矩形麦克风阵列性能
Table 5. Performance of initial array and optimized rectangular microphone array
初始阵列 优化阵列 频率
/kHz最大旁瓣水平/dB 分辨率
/m最大旁瓣水平/dB 分辨率
/m1 -12.545 0.302 -16.101 0.320 2 -11.786 0.146 -16.012 0.155 3 -11.798 0.094 -15.450 0.101 4 -11.800 0.068 -15.252 0.073 5 -11.806 0.053 -15.147 0.056 6 -6.523 0.042 -11.271 0.046 -
[1] Brooks T F, Marcolini M A, Pope D S. A directional array approach for the measurement of rotor noise source distributions with controlled spatial resolution[J]. Journal of Sound and Vibration, 1987, 112(1):192-197. doi: 10.1016/S0022-460X(87)80105-0 [2] Ahmad D V, Radhakrishnan S. Review of acoustic measurement techniques in wind tunnels[R]. AIAA-2000-2204, 2000. [3] Wei Q K. Compressive sensing based beamforming and its application in aeroacoustic experiment[R]. AIAA-2014-2918, 2014. [4] 周家检, 郝璇, 张卫民, 等.相阵列技术在民机机体气动噪声研究中的应用[J].空气动力学学报, 2016, 34(1):91-97. http://d.old.wanfangdata.com.cn/Periodical/kqdlxxb201601016Zhou J J, Hao X, Zhang W M, et al. Application of phased array technique in the research of civil airplane airframe noise[J]. Acta Aerodynamic Sinica, 2016, 34(1):91-97. http://d.old.wanfangdata.com.cn/Periodical/kqdlxxb201601016 [5] 陈正武, 王勋年, 李征初, 等.基于声学风洞的麦克风阵列测试技术应用研究[J].实验流体力学, 2012, 26(3):84-90. doi: 10.3969/j.issn.1672-9897.2012.03.016Chen Z W, Wang X N, Li Z C, et al. Application investigation of microphone array measuring and testing technique in anechoic wind tunnel[J]. Journal of Experiments in Fluid Mechanics, 2012, 26(3):84-90. doi: 10.3969/j.issn.1672-9897.2012.03.016 [6] Gade S, Hald J. Array designs optimized for both low-frequency NAH and high-frequency Beamforming[R]. SAE Technical Paper 2005-01-4014, 2015. [7] Underbrink J R. Practical considerations in focused array design for passive broadband source mapping applications[D]. Pennsylvania: Pennsylvania State University, 1995. [8] 黄奔.气动噪声源的麦克风阵列识别定位技术研究[D].绵阳: 中国空气动力研究与发展中心, 2014.Huang B. Investigation of aerodynamic noise sources identification technique based on microphone arrays[D]. Mianyang: China Aerodynamics Research and Development Center, 2014. [9] William M, Humphreys J R. Design and use of microphone direction arrays for aeroacoustic measurement[R]. AIAA-98-0471, 1998. [10] Thomas F B, William M, Humphreys J R. Effect of directional array size on the measurement of airframe noise components[R]. AIAA-99-1958, 1999. [11] Takeshi Ito. Aeroacoustic noise measurements in aerodynamic low-speed wind tunnels[C]//Proc of the 26th International Congress of the Aeronautical Sciences. 2008. [12] Sijtsma P. Phased array beamforming in wind tunnels[R]. NLR-TP-2006-732, 2006. [13] 陈宝, 李周复, 谭啸, 等.声衬试验段环境下航空声学定位试验技术研究[J].实验流体力学, 2015, 22(5):78-83. http://www.syltlx.com/CN/abstract/abstract10879.shtmlChen B, Li Z F, Tan X, et al. Investigation of aeroacoustic localization technique in lining test section[J]. Journal of Experiments in Fluid Mechanics, 2015, 29(5):78-83. http://www.syltlx.com/CN/abstract/abstract10879.shtml [14] 郭旺柳, 宋文萍, 许建华, 等.旋翼桨尖气动/降噪综合优化设计研究[J].西北工业大学报, 2012, 30(1):73-79.Guo W L, Song W P, Xu J H, et al. An effective aerodynamics/acoustic optimization of blade tip planform for helicopter rotors[J]. Journal of Northwestern Polytechnical University, 2012, 30(1):73-79. [15] 聂雪媛, 刘中玉, 杨国伟.基于Kriging代理模型的飞行器结构刚度气动优化设计[J].气体物理, 2017, 2(2):8-16. http://d.old.wanfangdata.com.cn/Periodical/qtwl201702002Nie X Y, Liu Z Y, Yang G W. Aircraft structure stiffness and aerodynamics optimization design based on Kriging surrogate model[J]. Physics of Gases, 2017, 2(2):8-16. http://d.old.wanfangdata.com.cn/Periodical/qtwl201702002 [16] Kanazaki M, Tanaka K, Jeong S, et al. Multi-objective aerodynamic optimization of elements' setting for high-lift airfoil using Kriging model[R]. AIAA-2006-1471, 2006. [17] Han Z H. Improving adjoint-based aerodynamic optimization via gradient-enhanced Kriging[R]. AIAA-2012-0670, 2012. [18] 刘俊.基于代理模型的高效气动优化设计方法及应用[D].西安: 西北工业大学, 2015. http://cdmd.cnki.com.cn/Article/CDMD-10699-1016909342.htmLiu J. Efficient surrogate-based optimization method and its application in aerodynamic design[D]. Xi'an: Northwestern Polytechnical University, 2015. http://cdmd.cnki.com.cn/Article/CDMD-10699-1016909342.htm [19] 韩忠华. Kriging模型及代理优化算法研究进展[J].航空学报, 2016, 37(11):3197-3225. http://d.old.wanfangdata.com.cn/Periodical/hkxb201611001Han Z H. Kriging surrogate model and its application to design optimization:a review of recent progress[J]. Acta Aeronautica et Astronautica Sinica, 2016, 37(11):3197-3225. http://d.old.wanfangdata.com.cn/Periodical/hkxb201611001 [20] 夏露, 王丹.基于Kriging自适应代理模型的气动优化方法[J].航空计算技术, 2013, 43(1):13-17. doi: 10.3969/j.issn.1671-654X.2013.01.004Xia L, Wang D. Aerodynamic optimization method based on Kriging adaptive surrogate model[J]. Aeronautical Computing Technique, 2013, 43(1):13-17. doi: 10.3969/j.issn.1671-654X.2013.01.004 [21] Zhang Y, Han Z H, Shi L X, et al. Multi-round surrogate-based optimization for benchmark aerodynamic design problems[R]. AIAA-2016-1545, 2016. [22] Thomas J M. Aeroacoustic measurements[M]. Berlin:Springer, 2002.