1、汽车气动噪声外辐射声场的数值仿真概要2013年(第35卷第1期汽车工程 Automotive Engineering汽车气动噪声外辐射声场的数值仿真牵郑拯宇1,李人宪1(1.西南交通大学机械工程学院,成都610031;2.重庆理工大学重庆汽车学院,重庆4000542013018摘要从LighthiH声类比理论出发,将流体动力学技术与边界元法结合起来,在某轿车边界元模型中,导入 流场边界脉动压力数据,并经转换和计算获得汽车表面附近的气动偶极子声源边界条件;采用直接边界元算法进行 汽车气动噪声外辐射声场的数值仿真。结果表明:轿车表面的偶极子声源强度随频率增大而降低;在120km/h车 速和2000
2、Hz频率时后视镜附近声场的气动噪声声压级可达78dB左右;在同一频率下,轿车在纵向对称面上的气动 声源辐射强度要大于地平面上的气动声源辐射强度。关键词:汽车气动噪声;偶极子声源;辐射声场;边界元法Numerical Simulation on the External Radiation SoundField of Vehicle Aerodynamic NoiseZheng Zhengyul2&Li Renxianl1.School ofMechanical Engineering,Southwest Jiaotong Univecsity,Chengdu 610031;2.Chongqln
3、g Automobile College,Chongqing Unb口ersity ofTechnology,Chongqing 400054AbstractBased on Lighthill acoustic analogy theory and combining boundary element method(BEMwith computational fluid dynamics(CFD,the pulsating pressure data of flow field boundary are input into the BEMmodel for apassenger cal,a
4、ndthe boundary conditions of aerodynamic dipole sound source near car surface are ob-tained through conversion and calculation.A numerical simulation on the external radiation sound field of vehicle aerodynamic noise is conducted by using direct BEM.The results show that the intensity of dipole soun
5、d source on car surface decreases with the increase of frequency,the sound pressure level of aerodynamic noise near rearview mirrors is up to some 78dB at 2,000Hz with a vehicle speed of 120km/h,and at the same frequency,the radiation intensity of aerodynamic sound source on the longitudinal symmetr
6、ic plane of vehicle is larger than that 013horizontal plane.Keywords:vehicle aerodynamic noise;dipole sound source;radiation sound field;BEM日!j吞气动噪声又称气流噪声。当汽车高速行驶时, 周围的空气在汽车表面附近形成流动状态十分复杂 的湍流场,而流场中的涡运动将在车身表面附近形 成一种压力脉动,这种压力脉动是诱发汽车产生气 动噪声的主要原因。研究表明:气动噪声是与车速 的6次方成正比五1。随着车辆工程技术的发展, 汽车其它噪声(诸如发动机噪声、传动系噪声
7、等得 到有效控制,而气动噪声则随着车速的不断提高而 增强,并逐步成为汽车高速行驶时的主要噪声源之 一。这种气动噪声对于轿车乘坐舒适性以及对高速 公路周边的环境噪声污染等指标产生重要影 响5J。国内外学者已在相关领域做了大量的研究 工作,并取得一定进展怕4J。但由于气动噪声是在+高等学校博士学科点专项科研基金(20100184110002、重庆市科委自然科学基金项目(CSTC,2008BB7341和重庆市教 委科学技术研究项目(KJl00S11资助。原稿收到日期为2011年4月6口,修改稿收到日期为2011年9月22日。万方数据郑拯宇,等:汽车气动噪声外辐射声场的数值仿真 里!二一个宽频带范围呈
8、连续分布,故从实车实道现场采 样数据中提取并分析气动噪声信号存在较大技术困 难;静音风洞数量稀少、成本高;而目前相关的数值 仿真技术则主要致力于车辆内部封闭空间的声场分 析,或针对车辆外部流场中有限个离散点的噪声信 号进行预测,或局限于选择某种基于稳态流场建立 起来的宽带声源模型理论对车辆表面气动声源强弱 的分布进行大致预判,而这些宽带声源模型在其理 论推导过程中做了一系列必要的简化假设,忽略了 湍流中涡旋在时间历程上对整体流场声能的影响, 不仅精度不高,而且不能获得基于车辆表面气动声 源的外辐射声场。上述研究均不能对车辆外部空间 气动噪声场进行较为全面、直观的认识和理解。 本文中从Light
9、hill声类比理论0。11o出发,利用 气动偶极子声源原理,将流体动力学数值分析技术 与边界元声场分析方法结合起来,在某轿车边界元 模型中,导人流场脉动压力数据并转换成气动偶极 子声源边界条件,采用直接边界元算法进行基于气 动偶极子声源的轿车气动噪声外辐射声场的数值仿 真分析,为进一步研究轿车车身优化设计和高速公 路沿线噪声环境控制的可行性方案提供数值依据。 文献中指出:气动噪声中四极子声源与偶极子声源 强度之比正比于马赫数的平方j,而地面行驶的汽 车运动仍属于低速运动(远小于马赫数,其四极子 源噪声强度远小于偶极子声源,故四极子源项(第1项亦可忽略不计。因此适合车辆流场的FwH方 程可简化为
10、等哪:等一砉【p,6差】 c2, 式(2表明:汽车高速行驶时,车身表面附近流 场脉动压力p。i所形成的车身表面偶极子声源决定了 汽车外部的气动噪声辐射声场。1.2边界元基本原理边界元方法(boundary element method,BEM是 分析声场问题最常用的方法,基本思想是基于格林 公式的应用,即将一个区域上的积分转化为该区域 边界上的积分,将问题的维数降低一维,从而使求解 的方程数目大大降低。1卜14J。边界元法分为直接边 界元法与间接边界元法两种,其中直接边界元法用 于分析封闭结构的内声场或外声场,以解决声波的 辐射和散射问题,其起点是外部Helmholz边界积分 方程:1Kght
11、hill声源理论与边界元基本原理 P(笕=巧p(y望学一G(戈,笔ds(, (31.1车辆气动噪声场中的声源理论针对流场中运动存在运动固壁的情况,Ffowcs Williams和Hawkings扩展了Lighthill方程的解,得 到如下Ffowcs Williams-Hawkings方程(简称FwH 方程12|:ic32p r 矗 警=堕 O x i O x 一 珈艿期+珈 舢鞠ji川。面。一面【p户瓦J+面【助%w面J (1 式中:p为流场压力脉动量;a。为声传播速度;瓦为 Lighthill张量pii为作用在运动固壁表面附近的压 力;6为狄拉克-Delta函数;风是未受扰动的流体 密度或
12、均值。式(1方程右边的第1项是四极子声源项,只存 在于运动固体表面之外;第2项表示由表面脉动压 力引起的声源,是偶极子声源项;第3项表示由表面 加速度引起的声源,是单极子声源项。第2、3项仅 在固体表面上产生。对于在空气流场中高速行驶的 车辆,车身表面可视为刚性,体积脉动量几乎为零, 所以单极子声源项(第3项积分可不必考虑;另外, 式中:p(茁和p(y分别为观测点x处的辐射压力和 声源点Y处的脉动声压;刀为辐射源表面.s上Y点外 法线向量。在其推导过程中,Helmholz积分方程严 格自动满足无界区域的Sommerfeld辐射条件: 1imf,f竺+i纫(,11:o (4 7L、or ,J式中
13、:r为声源指向观测点的矢量;p(,为helmholtz 方程的解;k=2wf/Co为波数/为声源自振频率。 由此,外部辐射声场中任意一点戈处的声压与 辐射源表面|s上的声压和声压梯度值之间就可以利 用式(3联系起来。本文中首先将所计算的车身表面附近流场的脉 动压力数据导入声振耦合软件SYSNOISE中,利用 式(5转换成声辐射压力;然后利用快速傅里叶变 换(FFTr将流场中时域脉动信号变换成频域信号, 再根据式(6求出压脉动频率5|,从而获得气动噪 声辐射场中的边界偶极子声源条件:%=p(y掣ds, (5万方数据汽车工程 2013年(第35卷第1期fanalysis=k=q忐 (6 式中:为辐
14、射声压;p(Y为流场中运动固壁边界上 脉动压力;为采样点数;矗为第k个采样点。最后利用直接边界元法,采用式(7计算辐射 声场中任一点P处的声压:PP=81P+b1h(7 式中12和b是网格数据导人SYSNOISE后所生 成的影响矩阵A和B中相对应的向量。2表面气动偶极子声源边界条件建立计算以Toyota Yaris 1.6RS轿车为实例,其车身 基本数据为:长3750ram、宽1695mm、高1545mm、 轴距2460mm、前轮距1460mm、后轮距1460mm,其 几何模型如图1所示。首先建立轿车外部空气流场模型,利用流体动力学数值分析工具,在非稳态下采用大涡模拟(LES湍流模型,对车速为
15、120km/h的轿车外部图1轿车几何模型流场进行仿真分析,算出轿车车身表面附近流场的 脉动压力数据。建立用以生成声源边界条件的轿车 表面边界元声学网格,将所获得的脉动压力数据导 入并在声学网格上转换成偶极子声源边界条件?其 中声学网格尺寸的大小由所确定的分析频率上限决 定。同时,由于耳道频率共振的原因,人耳最敏感的声波频率在10003000Hz,因此本文中轿车气动噪声分析频率范围设定在此区域内,采用最大单元长度为18ram的四边形声学网格,以满足声场中上限频率图2轿车车头声学网格(3000Hz的计算精度要求,轿车车头部分声学网格 如图2所示。在120km/h车速下频率设定为1000Hz时,由
16、轿车外流场边界脉动压力在轿车表面声学网格上所 转换生成的偶极子声源边界条件如图3所示。由图可见:由于轿车周围涡系和尾部涡系的相 互作用,在车体侧后部形成了压力脉动较强的复杂 紊流区,是轿车车身表面主要的分布噪声源,其声压偶极子声压级/dB592E+001一。e。叫刚圉 194E+001 995E+001I|,。e e +。 。,l l 597E+001圈 398E。001I;j图3偶极子声源边界条件(频率1000Hz,速度120km/h级(SPL在7196dB范围内,最高声源声压级为 95.92dB;前后轮附近流场紊乱,所产生的气动声源 对于整车噪声能贡献量较大,最高达95.1dB。在车 速相同的条
copyright@ 2008-2022 冰豆网网站版权所有
经营许可证编号:鄂ICP备2022015515号-1