+高级检索
 引 en
×

分享给微信好友或者朋友圈

使用微信“扫一扫”功能。

基于麦克风阵列小波分析的中低雷诺数尾缘噪声研究

  • 宋章辰

    机 构：

    北京航空航天大学陆士嘉实验室(航空气动声学工业和信息化部重点实验室),北京 100191

    ×
  • ，郭昊

    机 构：

    北京航空航天大学陆士嘉实验室(航空气动声学工业和信息化部重点实验室),北京 100191

    ×
  • ，刘沛清

    机 构：

    北京航空航天大学陆士嘉实验室(航空气动声学工业和信息化部重点实验室),北京 100191

    ×

北京航空航天大学陆士嘉实验室(航空气动声学工业和信息化部重点实验室),北京 100191；

An investigation into trailing edge noise of middle and low Reynolds number based on wavelet analysis of microphone array

  • SONG Zhangchen

    Affiliation：

    Beihang University (Key Laboratory of Aeroacoustics of Ministry of Industry and Information Technology), Beijing 100191 , China

    ×
  • ，GUO Hao

    Affiliation：

    Beihang University (Key Laboratory of Aeroacoustics of Ministry of Industry and Information Technology), Beijing 100191 , China

    ×
  • ，LIU Peiqing

    Affiliation：

    Beihang University (Key Laboratory of Aeroacoustics of Ministry of Industry and Information Technology), Beijing 100191 , China

    ×

Beihang University (Key Laboratory of Aeroacoustics of Ministry of Industry and Information Technology), Beijing 100191 , China；

作者简介:

宋章辰,男,博士。主要研究方向:翼型尾缘噪声、阵列测试技术。E-mail:songdbs_work@buaa.edu.cn;

郭昊,男,副教授,博士生导师。主要研究方向:湍流理论及飞机减阻降噪流动控制技术研究等领域。E-mail:guohao@buaa.edu.cn;

刘沛清,男,教授,博士生导师。主要研究方向:旋涡分离流与流动控制、飞行器大迎角空气动力学、鸭式布局涡系干扰与控制、高速层流控制技术、现代轻质高效螺旋桨设计与优化、飞行器低Re数流动机理与控制、大型飞机起飞着落气动性能、水上迫降性能、气动噪声研究等领域。E-mail:lpq@buaa.edu.cn

通讯作者:

刘沛清,E-mail:lpq@buaa.edu.cn

中图分类号:V211;TK263

文献标识码:A

DOI:10.19416/j.cnki.1674-9804.2022.01.003

参考文献 1
张青,闫国华,武耀罡,等.噪声适航标准严格化趋势探究[J].噪声与振动控制,2013,33(3):235-240.
查找原文
参考文献 2
中国民用航空总局适航审定司.中国民用航空规章第36部航空器型号和适航合格审定噪声规定[S].北京:中国民用航空总局,2005.
查找原文
参考文献 3
MIKHAIL Y Z,VICTOR K.Fly-over noise source localization during acoustic flight tests of advanced passenger aircraft:25th AIAA/CEAS Aeroacoustics Conference[C].Delft:AIAA,2019.
查找原文
参考文献 4
PARK C.Radiation enhancement by nonequilibrium in earth’s atmosphere [J].Journal of Spacecraft and Rockets,1985,22:27-36.
查找原文
参考文献 5
WERNER D.Almost 40 years of airframe noise research:what did we achieve[J].Journal of Aircraft,2012,47(2):353-367.
查找原文
参考文献 6
ALLEN C S,BLAKE W K,DOUGHERTY R P,et al.Aeroacoustic Measurements[M].[S.l.]:Springer,2002.
查找原文
参考文献 7
TAM C K W.Discrete tones of isolated airfoils[J].Journal of the Acoustical Society of America,1974,55(6):1173.
查找原文
参考文献 8
FINK M R.Prediction of airfoil tone frequencies[J].Journal of Aircraft 1975,12:118-120.
查找原文
参考文献 9
ARBEY H,BATAILLE J.Noise generated by airfoil profiles placed in a uniform laminar flow[J].Journal of Fluid Mechanics,1983,134:33-47.
查找原文
参考文献 10
Pröbsting S,SERPIERI J,SCARANO F.Experimental investigation of aerofoil tonal noise generation[J].Journal of Fluid Mechanics,2014,747(1):656-687.
查找原文
参考文献 11
RICCIARDI T R,ARIAS-RAMIREZ W,WOLF W R.On secondary tones arising in trailing-edge noise at moderate Reynolds numbers[J].European Journal of Mechanics-B/Fluids,2020,79:54-66.
查找原文
参考文献 12
徐亮,胡鹏,张永斌,等.可用于相干声源的快速反卷积声源成像算法[J].机械工程学报,2018,(23):82-92.
查找原文
参考文献 13
徐康乐,陈迎春,江渊,等.基于相控麦克风阵列的民机主起落架气动噪声源识别技术研究[J].空气动力学学报,2015,33(4):523-529.
查找原文
参考文献 14
周家检,郝璇,张卫民,等.相阵列技术在民机机体气动噪声研究中的应用[J].空气动力学学报,2016,34(1):91-97.
查找原文
参考文献 15
李征初,李勇,陈正武,等.CLEAN-SC算法在风洞声源定位与识别中的应用研究[J].实验流体力学,2016,30(3):104-109.
查找原文
参考文献 16
WILLIAM H,KEVIN J,THOMAS A,et al.Measurements of 26%-scale 777airframe noise in the NASA ames 40-by 80 foot wind tunnel:AIAA 2005-2810[R].USA;AIAA,2005.
查找原文
参考文献 17
PASCIONI K A,CATTAFESTA L N.Aeroacoustic measurements of leading-edge slat noise:22th AIAA/CEAS aeroacoustics Conference[C].Lyon:AIAA,2016.
查找原文
参考文献 18
SIJTSMA P.CLEAN based on spatial source coherence:13th AIAA/CEAS Aeroacoustics Conference[C].USA,AIAA,2007.
查找原文
参考文献 19
MALLAT S G.A theory for multiresolution signal decomposition:the wavelet representation[J].IEEE Transactions on Pattern Analysis & Machine Intelligence,1989,11(4):674-693.
查找原文
参考文献 20
SIJTSMA P,MERINO-MARTINEZ R,MALGOEZAR A M,et al.High-resolution CLEAN-SC:Theory and experimental validation[J].International Journal of Aeroacoustics,2017,16(4-5):274-298.
查找原文
目录contents

    摘要

    随着机体噪声逐渐成为民机降噪的关键,作为典型气动噪声源的翼型噪声是未来更低噪声民机必须关注的部分。中低雷诺数下尾缘离散噪声存在阶梯的频率速度关系、层流要求以及偶极子噪声等特征,其噪声机理涉及稳定性理论和反馈回路。尽管尾缘噪声的产生出现和强度都得到比较广泛的研究,但尾缘噪声的空间分布研究尚不充分。基于二维平面麦克风阵列,小波系数和CLEAN-SC算法对于尾缘噪声进行时频域分析。实验在北航D5气动声学风洞进行,使用Kevlar布构成的闭口实验段对于300 mm的NACA0012翼型进行噪声研究。通过实验中发现的固定频率范围内存在强度维持现象,基于小波分析拆解出对应的空间时域上声源强度的位置间歇性移动,提供可能存在的反馈回路内部更为复杂噪声机理的依据。

    Abstract

    Since airframe noise gradually plays a key role on civil aircraft noise reduction, airfoil noise, which is a typical source of aerodynamic noise, must be paid attention to in further noise reduction in civil aircraft industry in the future. The trailing edge discrete noise at low and medium Reynolds numbers has characteristics such as ‘ladder’ frequency speed relationship, laminar flow requirements, and dipole noise. Its noise mechanism is concerned with stability theory and feedback loop. Although the appearance and intensity of trailing edge noise have been extensively studied, the research on the spatial distribution of trailing edge noise is still insufficient. Based on wavelet coefficients and CLEAN-SC algorithm, a two-dimensional planar microphone array was applied to analyze the trailing edge noise in the time-frequency domain. The experiment was carried out in the Beihang D5 aeroacoustics wind tunnel. The closed experimental section composed of Kevlar cloth was used to conduct noise research on the NACA0012 airfoil with a chord length of 300 mm. Through the phenomenon of intensity maintenance in the fixed frequency range found in the experiment, the position of the sound source intensity in the space and time domain corresponding to the wavelet analysis disassembled intermittently moves, which further confirms the more complicated mechanism that may exist in the feedback loop.

  • 0 引言

  • 飞机噪声标准[1-2]越发严格使得飞机降噪成为研究热点,而随着发动机噪声降低,机体噪声愈加凸显[3]。翼型尾缘离散噪声是一个经典的噪声现象[4],其一般在中低雷诺数下受限于尾缘厚度是中低频噪声,且当采取人工转捩后压力面处于湍流边界层不会出现离散噪声[5-6]。噪声主频频率存在阶梯现象,其背后的机理一般认为是反馈回路导致[7-9],另一个现象即离散噪声存在等间距的频率分布,推测其可能是由于振幅调制[10]或者第二反馈回路[11]导致的。离散噪声机理存在时空变化,对于振幅调制存在间歇性,对于反馈回路存在反馈回路长度变化引起的空间位置变化,因此有必要研究声源强度的空间分布。

  • 麦克风阵列测试技术[12-14]作为一种确定声源强度信息在空间上的分布的噪声测量技术手段,被国内外学者广泛应用到气动声源噪声实验中[15-17]。基于麦克风阵列研究尾缘离散噪声的空间分布有助于进一步研究尾缘离散噪声机理。为了获得高分辨率的频域下空间分布解以及时间上的瞬时变化结果,本文引入了CLEAN-SC算法[18]以及小波分析[19],前者以较高的计算效率提高了对于声源强度空间分布精度,后者可以获得瞬时麦克风在不同频率上的幅值和相位信息。

  • 本文基于上述方法进行翼型自噪声的阵列识别研究,确定尾缘离散噪声特性。通过采用小波分析获得瞬时的各个麦克风频域的幅值和相位信息的复数解,最后在CLEAN-SC算法的基础上获得瞬时声源面上强度分布,进而分析不同离散噪声频率随时间的变化规律。

  • 1 实验设置

  • 1.1 风洞介绍

  • 本次实验在北京航空航天大学沙河校区陆士嘉实验室的D5气动声学风洞进行。北航D5气动声学风洞是一座低速、低湍流度、低噪声回流气动声学风洞,如图1所示。总体长25.58m,宽9.2m,高3.0m。D5风洞开口试验段风速为0m/s~80m/s,闭口试验段的风速为0m/s~100m/s,风洞来流湍流度小于0.08%。D5风洞消声室内部尺寸为长7m,宽6m,高6m,截止频率为200Hz,内壁由低频吸声隔声板制成,用于模拟自由声场环境,可以吸收99%以上的反射声,满足声传播的平方反比规律。

  • 图1 D5气动声学风洞

  • 1.2 闭口段模型布置

  • 本次实验采用NACA0012翼型,根据BANC#3的缩比模型,采用模型弦长为300mm,翼展为1 000mm。实验采用Kevlar布作为侧壁面,在兼顾流场信息稳定的同时保证了流场内气动噪声可以几乎不受阻碍地穿过布面到麦克风处。该方面已经得到比较广泛地应用,且已经在30P30N翼型等多组实验研究[31-32]中应用该方法。

  • 本次实验中选择G.R.A.S 40PH阵列麦克风构成的麦克风阵列作为主要的噪声采集和声源定位测量研究工具,翼型和麦克风阵列面的布置如图2所示。阵列正对着翼型旋转中心位置,距离为1 000mm。翼型旋转中心距离尾缘为180mm。实验工况设置如表1所示。AOA为来流攻角,U为来流速度。

  • 图2 实验布置示意图

  • 表1 工况设置

  • 1.3 阵列性能以及算法

  • 本次实验中使用了64通道麦克风阵列测试平台,其空间分布和基本性能如图3所示。

  • 图3 阵列空间布局以及基本性能

  • 本实验中采用的传统波束成形(Conventional Beamforming,简称CB)算法以及CLEAN-SC算法的原理简单介绍如下,具体可以参考SIJTSMA的文章[20]。对于CB算法,根据各个麦克风的频谱结果可以得到扫描面声源强度分布:

  • A=w*Rw
    (1)

  • 式中,A为扫描面上各个位置点的强度,w为权重向量,R为麦克风频域结果的互谱密度矩阵(CSM)。在本文中采用了去对角线(Diagonal Removal)的方法降低了麦克风自噪声的干扰,本文所有的CB都指代DR Beamforming。权重向量是通过声传播方程的格林函数解对应的指向向量归一化得到的,一般指向向量e是一个指向某一个扫描点的M×1维向量,M为麦克风数目,具体公式为:

  • em=Rm/Rcexp⁡-ikRm-Rc
    (2)

  • 式中,m为麦克风编号,Rm为扫描点到麦克风m的距离,Rc为扫描点到参考位置的距离,一般设置为阵列中心,k为波数。

  • 对于CLEAN-SC算法,其将互谱密度矩阵分为若干个迭代中声源强度最大位置处声源的求和,同时考虑到每次迭代中,声源强度最大位置点处于全部扫描面上的点存在的空间相干关系,将确保下一步的迭代中不再包含上一步位置的强度关系。具体的情况举例如下。

  • 在扫描面上基于CB算法用互谱密度矩阵R求得最大声源强度为Qmax,位置为Xmax,对应的指向向量为emax,对于扫描面上任意的指向向量e满足如下公式:

  • eRemax=eDemax
    (3)

  • 式中,m为麦克风编号,D为需要从R中单独取出对应位置Xmax的所有相干声源对于麦克风阵列互谱密度矩阵的影响程度,对于单点声源模型,可以定义Xmax处声源对于CSM贡献D满足公式:

  • D=Qmaxhh∘
    (4)

  • 式中,h为设定的一维向量。通过将公式(4)代入公式(3),求解迭代可以获得每一次R中单独分解的D,从而Ri+1=R-φD,最终求得所有最强声源点的位置和强度信息。φ为迭代系数。

  • 对于小波分析算法,采用的是复数Morlet小波基,计算结果是通过小波系数作为阵列麦克风瞬时接受到的信号幅值和相位,通过采样平均将得到的结果进行上述阵列算法求解瞬时声源面强度分布。

  • 2 实验结果与讨论

  • 2.1 离散噪声阶梯现象

  • 相比较单一麦克风由于空间位置关系导致各个不同离散噪声主频感受性不同,这里选择将阵列上各个麦克风自谱的平均值作为尾缘噪声测量结果,获得的平均噪声强度更为可信,如图4所示,在来流攻角为0°、1°和2°的情况下,可以发现其存在离散噪声频率随来流速度变化的“阶梯”现象。

  • 图4 阵列麦克风平均噪声强度频谱图

  • 为了显示对比方便,进一步提取各个角度下主频随来流速度变化进行作图,如图5所示,可以观察到目前存在两段“台阶”,疑似存在第三级但是由于此时雷诺数足够大,压力面层流边界层已经转捩,观察不到明显的离散噪声。出现阶梯“跳跃”现象的速度在13m/s~15m/s之间。

  • 图5 主频随来流速度变化

  • 2.2 阵列识别结果分析

  • 针对离散噪声存在的多峰现象,现有的通过单一麦克风研究结果中可以通过小波分析将主频时间上间歇性与对应的旋涡流动频率联系起来[10]。本实验结果中选取来流速度在18m/s下的结果进行详细分析。图6是麦克风测量离散噪声长时间平均结果以及存在间歇性的归一化C-Morlet小波系数结果。

  • 图6 离散主频幅值和小波系数结果

  • 针对452Hz和516Hz这两个频率的噪声进行小波系数求解,麦克风频域上随时间变化的复数小波系数作为幅值和相位的出处,获得流向上平均强度随时间的变化,如图7所示。

  • 选择最强声源点的流向位置为研究对象作图,如图8所示,可以发现,对于452Hz和516Hz的离散噪声,其在空间位置上存在相对较大的位置移动,且516Hz的移动要比452Hz的要剧烈。

  • 基于其在流向上位置移动的近似速度进行多普勒频率修正,可以确认,对于516Hz的噪声,如果设置对应的速度为2m/s~6m/s且远离麦克风,那么在18m/s来流风速下,对应产生的频率为464.4Hz~384Hz,存在较小强度的408Hz和较大强度的452Hz,对于反馈回路内的离散噪声的多个频率情况,很有可能是这种声源移动所导致的。

  • 图7 CLEAN-SC算法下流向上展向平均强度随时间变化

  • 图8 声源流向上相对于弦长位置随时间变化

  • 3 结论

  • 本文针对NACA0012翼型尾缘离散噪声进行研究,实验中基于多个阵列麦克风结果即在来流迎角2°下离散噪声频率改变而幅值基本不变的现象,基于小波分析得到麦克风时域上相位和幅值结果,通过CLEAN-SC算法捕捉声源在空间位置上强度分布随时间的变化。通过分析特定频率的噪声随时间变化,可以确认在不同频率存在声源位置的变化,而这一变化可能是对于高频高强度的噪声产生低频低强度的噪声,最终反映到麦克风接受的噪声上。

  • 参考文献

    • [1] 张青,闫国华,武耀罡,等.噪声适航标准严格化趋势探究[J].噪声与振动控制,2013,33(3):235-240.

    • [2] 中国民用航空总局适航审定司.中国民用航空规章第36部航空器型号和适航合格审定噪声规定[S].北京:中国民用航空总局,2005.

    • [3] MIKHAIL Y Z,VICTOR K.Fly-over noise source localization during acoustic flight tests of advanced passenger aircraft:25th AIAA/CEAS Aeroacoustics Conference[C].Delft:AIAA,2019.

    • [4] PARK C.Radiation enhancement by nonequilibrium in earth’s atmosphere [J].Journal of Spacecraft and Rockets,1985,22:27-36.

    • [5] WERNER D.Almost 40 years of airframe noise research:what did we achieve[J].Journal of Aircraft,2012,47(2):353-367.

    • [6] ALLEN C S,BLAKE W K,DOUGHERTY R P,et al.Aeroacoustic Measurements[M].[S.l.]:Springer,2002.

    • [7] TAM C K W.Discrete tones of isolated airfoils[J].Journal of the Acoustical Society of America,1974,55(6):1173.

    • [8] FINK M R.Prediction of airfoil tone frequencies[J].Journal of Aircraft 1975,12:118-120.

    • [9] ARBEY H,BATAILLE J.Noise generated by airfoil profiles placed in a uniform laminar flow[J].Journal of Fluid Mechanics,1983,134:33-47.

    • [10] Pröbsting S,SERPIERI J,SCARANO F.Experimental investigation of aerofoil tonal noise generation[J].Journal of Fluid Mechanics,2014,747(1):656-687.

    • [11] RICCIARDI T R,ARIAS-RAMIREZ W,WOLF W R.On secondary tones arising in trailing-edge noise at moderate Reynolds numbers[J].European Journal of Mechanics-B/Fluids,2020,79:54-66.

    • [12] 徐亮,胡鹏,张永斌,等.可用于相干声源的快速反卷积声源成像算法[J].机械工程学报,2018,(23):82-92.

    • [13] 徐康乐,陈迎春,江渊,等.基于相控麦克风阵列的民机主起落架气动噪声源识别技术研究[J].空气动力学学报,2015,33(4):523-529.

    • [14] 周家检,郝璇,张卫民,等.相阵列技术在民机机体气动噪声研究中的应用[J].空气动力学学报,2016,34(1):91-97.

    • [15] 李征初,李勇,陈正武,等.CLEAN-SC算法在风洞声源定位与识别中的应用研究[J].实验流体力学,2016,30(3):104-109.

    • [16] WILLIAM H,KEVIN J,THOMAS A,et al.Measurements of 26%-scale 777airframe noise in the NASA ames 40-by 80 foot wind tunnel:AIAA 2005-2810[R].USA;AIAA,2005.

    • [17] PASCIONI K A,CATTAFESTA L N.Aeroacoustic measurements of leading-edge slat noise:22th AIAA/CEAS aeroacoustics Conference[C].Lyon:AIAA,2016.

    • [18] SIJTSMA P.CLEAN based on spatial source coherence:13th AIAA/CEAS Aeroacoustics Conference[C].USA,AIAA,2007.

    • [19] MALLAT S G.A theory for multiresolution signal decomposition:the wavelet representation[J].IEEE Transactions on Pattern Analysis & Machine Intelligence,1989,11(4):674-693.

    • [20] SIJTSMA P,MERINO-MARTINEZ R,MALGOEZAR A M,et al.High-resolution CLEAN-SC:Theory and experimental validation[J].International Journal of Aeroacoustics,2017,16(4-5):274-298.

  • 基本信息

    中图分类号: V211;TK263
    文献标识码: A
    DOI: 10.19416/j.cnki.1674-9804.2022.01.003

    基金信息

    国家自然科学基金(12072016, 11772033, 1217021666 和 11721202)；

    引用信息

    宋章辰,郭昊,刘沛清.基于麦克风阵列小波分析的中低雷诺数尾缘噪声研究[J].民用飞机设计与研究,2022(1):38-44.

    SONG Z C, GUO H, LIU P Q. An investigation into trailing edge noise of middle and low Reynolds number based on wavelet analysis of microphone array[J].Civil Aircraft Design and Research,2022(1):38-44(in Chinese).

    稿件历史

  • 参考文献

    • [1] 张青,闫国华,武耀罡,等.噪声适航标准严格化趋势探究[J].噪声与振动控制,2013,33(3):235-240.

    • [2] 中国民用航空总局适航审定司.中国民用航空规章第36部航空器型号和适航合格审定噪声规定[S].北京:中国民用航空总局,2005.

    • [3] MIKHAIL Y Z,VICTOR K.Fly-over noise source localization during acoustic flight tests of advanced passenger aircraft:25th AIAA/CEAS Aeroacoustics Conference[C].Delft:AIAA,2019.

    • [4] PARK C.Radiation enhancement by nonequilibrium in earth’s atmosphere [J].Journal of Spacecraft and Rockets,1985,22:27-36.

    • [5] WERNER D.Almost 40 years of airframe noise research:what did we achieve[J].Journal of Aircraft,2012,47(2):353-367.

    • [6] ALLEN C S,BLAKE W K,DOUGHERTY R P,et al.Aeroacoustic Measurements[M].[S.l.]:Springer,2002.

    • [7] TAM C K W.Discrete tones of isolated airfoils[J].Journal of the Acoustical Society of America,1974,55(6):1173.

    • [8] FINK M R.Prediction of airfoil tone frequencies[J].Journal of Aircraft 1975,12:118-120.

    • [9] ARBEY H,BATAILLE J.Noise generated by airfoil profiles placed in a uniform laminar flow[J].Journal of Fluid Mechanics,1983,134:33-47.

    • [10] Pröbsting S,SERPIERI J,SCARANO F.Experimental investigation of aerofoil tonal noise generation[J].Journal of Fluid Mechanics,2014,747(1):656-687.

    • [11] RICCIARDI T R,ARIAS-RAMIREZ W,WOLF W R.On secondary tones arising in trailing-edge noise at moderate Reynolds numbers[J].European Journal of Mechanics-B/Fluids,2020,79:54-66.

    • [12] 徐亮,胡鹏,张永斌,等.可用于相干声源的快速反卷积声源成像算法[J].机械工程学报,2018,(23):82-92.

    • [13] 徐康乐,陈迎春,江渊,等.基于相控麦克风阵列的民机主起落架气动噪声源识别技术研究[J].空气动力学学报,2015,33(4):523-529.

    • [14] 周家检,郝璇,张卫民,等.相阵列技术在民机机体气动噪声研究中的应用[J].空气动力学学报,2016,34(1):91-97.

    • [15] 李征初,李勇,陈正武,等.CLEAN-SC算法在风洞声源定位与识别中的应用研究[J].实验流体力学,2016,30(3):104-109.

    • [16] WILLIAM H,KEVIN J,THOMAS A,et al.Measurements of 26%-scale 777airframe noise in the NASA ames 40-by 80 foot wind tunnel:AIAA 2005-2810[R].USA;AIAA,2005.

    • [17] PASCIONI K A,CATTAFESTA L N.Aeroacoustic measurements of leading-edge slat noise:22th AIAA/CEAS aeroacoustics Conference[C].Lyon:AIAA,2016.

    • [18] SIJTSMA P.CLEAN based on spatial source coherence:13th AIAA/CEAS Aeroacoustics Conference[C].USA,AIAA,2007.

    • [19] MALLAT S G.A theory for multiresolution signal decomposition:the wavelet representation[J].IEEE Transactions on Pattern Analysis & Machine Intelligence,1989,11(4):674-693.

    • [20] SIJTSMA P,MERINO-MARTINEZ R,MALGOEZAR A M,et al.High-resolution CLEAN-SC:Theory and experimental validation[J].International Journal of Aeroacoustics,2017,16(4-5):274-298.

    • 出版刊物

    微信公众号二维码

    手机版网站二维码

    本日访问次数 您是第位访问者

    民用飞机设计与研究 ® 2024 版权所有
    技术支持：北京勤云科技发展有限公司
    我要投稿 投稿指南 联系我们 二维码
    TOP