首页/文章/ 详情

【整车风噪开发解析】——有限元法风噪仿真关键技术解密(3)

4月前浏览8970

本文摘要:(由ai生成)

本文探讨了使用可压缩CFD结果进行车辆风噪仿真的方法,重点介绍了波数分解和薄膜模态分解技术,以区分和提取声学贡献和湍流压力贡献。文章还讨论了这些技术在Actran中的应用,及如何将分解后的压力加载到车内模型进行内声场计算。车内噪声以声学贡献为主。最后总结了压力分解的目的与特点。懿朵科技为整车风噪开发提供全面解决方案,包括指标定义、仿真优化、测试整改等。

车辆高速行驶在车窗上产生的压力可分解为湍流壁面脉动压力和声壁面脉动压力,湍流部分与湍流漩涡和边界层结构相关的小尺度结构有关,声学部分由汽车附近的气动声源产生,具有较大的波长特征。从量级上来看,尽管声学振幅与湍流贡献相比较小,但车窗对声学部分贡献更加敏感。  
 
图 1.湍流压(左)和声压(右)  
 
有限元法进行风噪仿真 主要有两种策略。第一种方法,使用不可压缩CFD来确定湍流压力脉动,用气动声学类比(Lighthill,Mohring)以确定声压脉动;另一种方法,采用可压缩的CFD方法解决流体动力学和声学问题, 该方法需要更多的计算资源,能单独进行气动噪声分析,但不能区分湍流压和声压两者的贡献。  
前两篇文章对第一种方法进行了详细介绍,本文将主要介绍利用可压缩CFD结果进行风噪仿真的方法,以及从压力场中识别声学贡献的两种技术:波数分解(WND)和薄膜模态分解(PMD)。  

01

波数分解(WND)

波数分解数学上是对平面表面压力场进行空间傅里叶变换,得到含有对应于不同波长的波数成分的信息谱结果,从中识别“声波”和“湍流波”。波数分解基于以下几个假设:压力场在空间上是周期性的谐波,粘度忽略不计。  
①空间周期性:  
 
②谐波求解:  
 
因此,可以将压力场投影为矩形模态  
 
 
图 2 压力场投影矩形模态  
假设压力场是完全对流的  
 
可得到以下色散关系:  
 
得到湍流压力贡献在波数谱中位于为垂直于流体速度的直线区域(图2绿色直线)。  
 
对于声学信号,根据亥姆霍兹方程:  
 
可得到以下色散关系:  
 
得到声学贡献集中在以(kx,ky)=(0,0)为中心的圆上(图2蓝色圆),半径为 2πf/c 。  
 
 
图 3 湍流贡献和声学贡献在波数谱中的位置(左湍流贡献,右声学贡献)  
 
因此,可以根据波数谱(kx,ky)中的位置来区分湍流压和声压。声学贡献为理论轨迹圆内的压力之和,湍流贡献为WND域的剩余部分。  
 
图 4压力场波数分解为湍流压和声压  
 
        在Actran中通过PRESSURE DECOMPOSITION命令可对所有关注频率执行波数分解操作,湍流压和声压在波数谱中的位置可通过瀑布图得到,同时能分别得到声压激励与湍流激励对玻璃的贡献量,可以看到对车窗压力的主要贡献是由于湍流激励引起(图5)。  
 
图 5车窗压力与声学贡献  

02

薄膜模态分解(PMD)

车窗表面压力场(声压或湍流压力场)可以表示为:  
 
其中,P车窗压力, ai为参与因子,xp为薄膜模态。  
薄膜模态分解(PMD)通过提取玻璃表面的薄层空气的声学模态进行求解,将压力场投影到声学薄膜模态中并检索参与因子,根据模态参与因子的不同来区分声学贡献和湍流贡献。薄膜模态参与因子可以在模态分解图中进行查看(如图6,颜色深浅代表了模态参与因子的大小)。  
 
图 6 薄膜模态分解图  
 
在薄膜模态分解图中可以看到两个区域,如图7:  
声压区域:根据定义,表面压力场中的声压分量应该具有基于声波波长的模态,所以声压风量应该为接近所考虑声波频率的固有频率的模态提供更高的参与因子。因此,在薄膜模态分解图中,声学分量应该集中在斜率为1的轨迹附近。  
湍流压力区域:表面压力场的湍流分量应具有近似于频率除以马赫数(图中Ma≈0.11)的波长,故湍流分量应该对于具有相似波长的薄膜模态提供更高的参与因子,位置集中在斜率等于马赫数的轨迹附近。  
因此可以根据模态分解图中两个区域的位置对湍流压和声压进行分解。  
 
图 7 薄膜模态分解轨迹识别  
 
波数分解只能应用于矩形表面,薄膜模态分解可以应用任何形状表面,更适应实际的车窗。在完整车窗上取矩形子区域(图8),可以看到在相同表面形状上执行WND和PMD计算,分解结果非常相似。但与全尺寸车窗的PMD结果存在差异,主要是因为全尺寸PMD使用较大平面进行分解计算,表面涉及其他湍流产生的附加声学贡献。  
 
图 8 两种分解方法声学贡献量对比  
 
需要注意的是,薄膜模态的数量对PMD结果影响很大,实际上,声学贡献用几种模态就可以表示(如图9,73和652模态数声学贡献结果一致),但小波长的湍流压力贡献需要大量模态进行评估,尽管模态提取计算不需要大量CPU和RAM资源,但压力场投影积分计算量很大。可以将薄膜模态提取限制为分析频率的1.5倍减少误差,也可以使用少量的薄膜模态来识别声学贡献,滤除声学贡献后,通过压力场推导出湍流压力贡献,如图10,少量模态数(73modes)获得的湍流压力不足以准确表示湍流贡献,增加模态数后(652modes)湍流贡献与采用少量模态数(73modes)滤除声学贡献后的湍流压力结果相似。  
 
图 9 薄膜模态数对结果的影响  
 
 
图 10薄膜模态数对湍流压力影响  
 
 

03

内场声学计算

 
利用可压CFD结果进行风噪仿真的技术路线如下:  
①可压CFD瞬态计算,获得车窗表面压力分布;  
②通过薄膜模态分解/波数分解将车窗表面压力分为湍流压力(TWPF)和声压(AWPF);  
③将分解后的两个压力载荷通过耦合面加载到车内振动声学模型上,计算车内声传播。  
 
图 11 可压缩风噪仿真技术路线  
 
    对比内场结果可以发现,630Hz以上的中高频,噪声贡献主要以声学贡献为主,630Hz以下则是湍流压力贡献为主,而前文压力分解结论为车窗表面压力贡主要由湍流压力引起,究竟是什么原因导致车内两者贡献产生了差异?  
 
图 12 内场仿真与试验结果对比与声学贡献  
 
这里需要提到弯曲波长、声学波长和声学波长 三个概念,  
弯曲波长:   
 
声学波长:     
 
湍流波长:     
 
外部的激励可以激发车窗玻璃中的弯曲波,当外部激励波长等于玻璃弯曲波长时,将会发生吻合效应,引起玻璃振动,向外辐射噪声。声波在全频范围内波长与车窗玻璃弯曲波长接近,更容易激发玻璃中的弯曲波,通过玻璃向外辐射声波,故声压脉动激励通过玻璃的传递性较高,车窗对声学贡献也更敏感;湍流波长除在较低频率的耦合频率处与玻璃弯曲波长接近外,其他频率范围远大于玻璃弯曲波长,所以湍流脉动激励通过玻璃的传递性较差。因此,虽然车外湍流压力幅值远大于声压,但车内噪声的贡献主要以声学贡献为主,具体表现为低频时湍流脉动压力占主导,高频时声压脉动为主。  

04

总结

本文主要介绍了波数分解、薄膜模态分解两种压力分解方法,对比了两种方法的差异,以及采用可压缩CFD结果进行风噪仿真的方法,总结压力分解目的与特点:  
1)可压缩CFD结果的后处理,从整体压力信号中提取声学信号;  
2)作为NVH实验部门的输入数据 (避免声学波长与侧窗结构波长重合);  
3)作为CAE仿真部门的输入数据。  

懿朵科技为客户提供国际领先水平整车风噪开发解决方案:(指标定义,仿真分析与优化,测试与整改)

如需了解专题其他文章,点击以下链接查看:

【整车风噪开发解析】——专题一  
【整车风噪开发解析】——SNGR在整车风噪开发中的高效应用

【整车风噪开发解析】——有限元法风噪仿真关键技术解密(1)

【整车风噪开发解析】——有限元法风噪仿真关键技术解密(2)

敬请期待:《动态密封开发与控制策略》

                  《LBM在整车风噪开发中的应用》

来源:懿朵科技
ACT振动气动噪声湍流汽车声学理论NVH控制
著作权归作者所有,欢迎分享,未经许可,不得转载
首次发布时间:2024-04-20
最近编辑:4月前
懿朵科技
签名征集中
获赞 22粉丝 14文章 81课程 0
点赞
收藏
作者推荐

【整车风噪开发解析】——SNGR在整车风噪开发中的高效应用

本文摘要:(由ai生成)SNGR方法是一种高效气动噪声计算技术,它利用稳态流场结果进行噪声分析,可快速识别噪声源位置和评估不同设计的噪声水平,适合中高频段预测,对CFD网格尺度不敏感,适用于整车前期风噪水平快速评估。懿朵科技利用该方法进行整车风噪项目开发,实现了模块化和自动化,提高了计算效率,并为整车风噪开发提供了国际领先解决方案,包括指标定义、仿真分析优化、测试整改等,但在低频噪声预测上存在局限。 某车企的私人讨论组又发出了“叮咚”声! 小甲:“听说了吗?今天CFD和造型部门又干起来了!”小乙:“啊….见怪不怪咯!一天到晚提各种噪声问题,造型设计部门能不火么?”小丙:“其实也不是,主要是CFD那边效率太低了啊!开发周期那么紧张,造型那边等一版结果黄花菜都凉了,能不火么!”小丁:“哎…难办哦!”“哈哈哈…我们好好吃瓜就行啦!”“……”“好奇怪啊,最近CFD和造型咋这么和谐啊?”“你不知道么?SNGR来了啊!” SNGR难道有通天的能力?SNGR(Stochastic Noise Generation and Radiation)方法采用稳态流场结果进行气动噪声计算, 一天内完成一轮优化方案的分析,高效完成外造型优化方案的快速评估;极大地缩短了计算时间增加效率,用极少的资源就可对外造型优化方案进行风噪水平定性评估,非常适合整车前期外造型阶段风噪评估和优化。01 SNGR的原理和流程 (1)SNGR的原理SNGR方法采用雷诺平均Navier-Stokes(RANS)计算得到的时均流场的速度、湍流动能分布特性,通过添加随机扰动的方法重新合成含有时间项的流场数据,从而实现从RANS结果识别Lighthill体积源。SNGR方法对湍流脉动速度 采用Bailly随机模型,不可压缩空间湍流速度场可由傅里叶空间内N个傅里叶模态求解得到,具体如下: ACTRAN的声类比方法计算声源和噪声传播,湍流动能谱可以写成以下形式: 由于RANS为稳态结果,为了保证可以输出频域的噪声结果,需要添加与时间相关的项(图1中的随机函数),从而合成脉动速度。 图1 SNGR基本原理(2)SNGR的流程基于稳态CFD和Actran的风噪仿真流程如图2所示,为了计算声源项,需要从RANS计算中输出以下物理量:速度矢量、湍动能和湍流耗散率。在此基础上,湍动能谱可以使用von Karman-Pao 谱或用户自动义谱处理。 图2 SNGR仿真流程02 SNGR参数的影响 (1)湍流阀值湍流阀值的设定影响声源项考虑区域与湍流速度合成的时间,在阀值之上的源项将被考虑到声源计算中,通过设置湍流阀值可只考虑湍流贡献最大的区域,提高CPU和内存使用率,减小计算时间。 a 原始湍动能 b 设定湍流阀值后湍动能 图3 湍动能对比 湍流阀值增加,考虑的源项区域减小,计算时间减小。对比如下5种不同湍流阀值计算结果和计算资源,最大差值为2dB。实际应用中可根据计算资源和精度要求合理设置湍流阀值,减少仿真计算时间。表1 不同湍流阀值计算资源消耗 图4 不同湍流阀值结果(2)合成速度场(声源)的数量随着合成速度场数量的增加,仿真计算时间将会线性增加,计算结果也将逐渐趋于一致,实际应用中可以根据计算资源需求和精度选择合成速度场数量。如图5,合成速度场的数量为20和60在500Hz最大误差为3dB,合成速度场数量为40误差小于1dB。 图5 不同合成速度场数量结果(3) CFD网格SNGR方法预测风噪时对CFD网格的尺度不敏感,不需要采用特别密集的网格来进行稳态CFD分析,仅需CFD稳态流场的收敛和湍动能的准确统计即可。当湍流场轮廓 (等值面) 是平顺的,可以认为计算得到的湍流场是正确的。如果云图并不平顺,说明CFD网格还需要优化。 图6 不同CFD网格湍动能和耗散率 图7 不同CFD网格结果在瞬态CFD模拟中,当CFD的网格尺度过大,数值耗散较大时,CFD的数据就很难反映出高频的声波信息。如图8。网格尺度与截止频率关系为: 其中,F是声学截止频率,ε是湍动能耗散率,∆是网格尺度。SNGR方法则不需要考虑CFD网格带来的截止频率的影响,4mm的CFD网格可以很好的预测400Hz~3000Hz频段汽车风噪声大小。因此,中高频段利用Actran中 SNGR方法快速、高效且准确的预估风噪声大小是一个不错选择。需要注意的是SNGR方法在低频有些信息缺失,无法准确捕捉低频结果。 图8 仿真实验结果对比03 SNGR的二次开发 懿朵科技利用SNGR方法进行过多个整车风噪项目的开发,结合自身风噪开发和SNGR方法的应用经验,将稳态CFD+SNGR方法模块化,进行联合二次开发。整个流程只需要在计算前将流场面网格和声场网格进行导入,便可自动实现网格加密、体网格生成、计算以及后处理,从流场到声场的计算也不需要人为的干预,流场计算完成后便可自动进入声场计算,界面会实时监控和显示整个计算流程进度(绿色表示进度已完成、黄色表示正在运行),实现了一键自动计算和后处理,减少人为干预易出错的同时,减少1/3的工作量,再次提高计算效率。目前已经应用于懿朵科技风噪开发中。 图9 流场和声场计算模块04 SNGR的实际应用 1、三菱汽车。基于SNGR + ICFD的 Lighthil 声源识别 进行噪声源定位并同风洞中麦克风声阵列测试对比,识别声源位置与麦克风声阵列结果基本一致。 图10 噪声源识别2、Volvo汽车。采用CFD软件与Actran混合计算方法研究V70车型风噪声,声学计算分别采用了稳态CFD结果(SNGR)和瞬态LES结果作对比,仿真与试验有良好的一致性。 图11 Volvo案例3、Magna。采用稳态CFD+Actran混合方法(SNGR),研究外造型对侧窗表面风噪贡献,对比不同后视镜对内部风噪风噪影响,并通过稳态SNGR方法和非稳态CFD结果+声学计算两者结合预测乘员舱内部风噪。 图12 Magna 仿真模型 图13 声学计算结果4、Volvo。采用SNGR方法对重型卡车进行噪声源识别,对比不同后视镜方案对风噪的影响,利用SNGR得到的结果与非稳态DES混合方法趋势一致。 图14 Volvo卡车案例05 SNGR总结 SNGR方法用时间较短的RANS结果进行气动噪声计算,可极大地节省计算资源和时间,可用于汽车风噪、起落架噪音、HVAC管道噪声等各个方面。在汽车风噪方面可以通过SNGR方法快速识别噪声源位置,同时可以快速评估不同外造型设计的噪声相对水平,相比利用非稳态CFD结果进行气动声学计算,SNGR方法优势主要有以下几点:(1)对CFD网格尺寸不敏感,降低了气动噪声预测对CFD的要求,大幅缩短优化周期; (2)采用稳态CFD计算,对计算资源要求低,适用不同造型设计的噪声相对水平快速评估; (3)稳态CFD+SNGR混合方法能快速识别噪声源位置; (4)SNGR方法不受湍流截止频率影响,可以在较少的计算资源中实现中高频段的精确模拟。 SNGR方法虽然在减小计算时间、缩短开发周期上有极大地优势,但由于SNGR方法采用冯卡门谱合成湍流速度场,用其预测噪声的绝对值,需要经过多次修正拟合,且SNGR在低频有部分信息缺失,难以准确捕捉低频结果。因此,SNGR方法更适用于多方案之间噪声相对水平的快速评估。 懿朵科技为客户提供国际领先水平整车风噪开发解决方案:(指标定义,仿真分析与优化,测试与整改)如需了解专题其他文章,点击以下链接查看: 【整车风噪开发解析】——综述篇 【整车风噪开发解析】——噪声类型与开发策略 敬请期待: 《有限元法整车风噪开发关键技术解密》 《LBM在整车风噪开发中的应用》 来源:懿朵科技

未登录
还没有评论
课程
培训
服务
行家
VIP会员 学习 福利任务 兑换礼品
下载APP
联系我们
帮助与反馈