汽车驱动桥NVH 性能分析与优化
虽然上述研究对汽车驱动桥NVH性能分析与优化做了很多工作,取得许多成果,但仍然存在一些不足:(1)系统研究不足,仅从CAD/CAE技术和NVH实车测试等分别展开研究,没有根据以上研究方法特性,建立系统的分析方法;(2)CAD/CAE技术分析法所使用的软件各异,具体分析过程各异,分析结果的可信度各异,没有对不同软件分析过程进行优劣分析;(3)理论分析与实际效果验证未充分关联,以致理论分析结果可信度待定。
本文中主要通过建立驱动桥NVH性能系统分析的流程及方法,对模态分析、噪声辐射分析等CAE分析所使用软件及其分析过程方法进行研究,以选择恰当分析软件及方法,并通过实车驱动桥NVH最终测试来验证以上两项研究目标是否达成。驱动桥NVH性能系统分析流程如图1所示,本文中以某客车在60~65 km/h加速行驶工况出现噪声大的问题为例进行分析研究。
1 模态分析
对后驱动桥进行模态分析,目的是得到各阶模态频率,来确认其是否与其他激励源产生共振。前期研究结果表明,后桥噪声主要是主减速器齿轮啮合冲击通过轴承传至后桥壳产生振动引起的辐射噪声,差速器在普通工况下一般不起作用[24],本文中主要是对后桥壳进行模态分析。
1.1 有限元建模
采用UG软件系统建模,网格划分过程中,主减速器壳选取四面体单元划分,单元质量主要控制参数如表1所示,最后给各个部件赋相应的厚度和材料属性,如表2所示。
表1 单元质量控制参数
1.2 有限元分析结果
前期研究得出驱动桥噪声频率主要分布在2 000 Hz以内[24],将建好的驱动桥壳模型导入ansys workbench,设置频率1~2 000 Hz,计算固有频率1~2 000 Hz的模态。计算得到各阶模态振型及其对应的频率,根据经验,4阶模态后对共振的影响很小,在此仅列出1~4阶模态,如图3~图6所示。
表3列出了1~4阶模态频率以及振型描述。来自路面的随机激励一般在几赫兹到十几赫兹,而该后桥最低固有频率为367 Hz,所以后桥桥壳模态不会被路面激励激起产生共振,而该车型无发动机,也不会激起后桥模态。
主减速器齿轮啮合频率由式(1)计算得到。
在62.3 km/h加速行驶工况时,主减速器齿轮啮合频率为367 Hz。根据图3分析结果可知驱动桥壳1阶模态为367 Hz,主减速器齿轮1阶啮合频率与桥壳的第1阶模态重合,发生了共振,是该车型在60~65 km/h左右加速直线工况后桥噪声较大的根本原因。在后续的计算中,关注驱动桥壳第1阶模态频率处的振动噪声。
表3 2 000 Hz以内各阶模态频率及振型
2 振动响应分析
振动响应分析的目的是确定响应较大部位,以实现后续精准优化。频率响应分析是指结构对某载荷(可以是冲击载荷,也可能是一频率在一定范围内的载荷)的响应。根据驱动桥噪声机理,以及驱动桥NVH性能分析需要,在进行频率响应分析前,需要先计算其轴承的载荷。
2.1 差速器轴承载荷计算
为确定后桥内轴承的载荷,须先对准双曲面齿轮进行啮合受力分析,图7为主动齿轮齿面受力图。该主动小齿轮旋向为左旋,沿X轴方向看去,汽车前进时旋转方向为顺时针方向,被动齿轮作用在主动齿轮节锥面上齿面宽中点M处齿面法向力为F T。在法平面内,F T可分解为位于节锥切面内的力F N和与OM垂直的力F f,F N和F f相互垂直,F f和F T之间的夹角为法向压力角α。在节锥的切面内,F f可分解为沿节锥母线即OM方向的力F S和切线方向的圆周力F,F S和F也相互垂直,F和F f之间的夹角则为螺旋角β。
根据图7主动齿轮受力分析可得
表4 主被动齿轮所受轴向力和径向力
主减速器轴承布置如图8所示,其中包括A、B、C和D 4个圆锥滚子轴承。尺寸a为A轴承受力中心至主动齿轮齿面节圆距离,尺寸b为A与B轴承受力中心距离,尺寸c为C轴承受力中心与被动齿轮齿面节圆距离,尺寸d为D轴承受力中心与被动齿轮齿面节圆距离,D m1为主动齿轮齿面节圆直径,D m2为被动齿轮齿面节圆直径。
在确定了齿面上的圆周力、轴向力和径向力以后,根据齿轮和轴承布置的相关尺寸和参数,在整车坐标系中对其进行受力分析,图9~图12分别为主动齿轮轴、差速器轴受力分析简化图。
F A X、F A Y分别为轴承A所受的轴向力、径向力,F B X、F B Y分别为轴承B所受的轴向力、径向力,F A Z、F B Z分别为A、B轴承在Z向承受的力,F C X、F C Y分别为轴承C所受的轴向力、径向力,F D X、F D Y分别为轴承D所受的轴向力、径向力,F C Z、F D Z分别为C、D轴承在Z向承受的力。
通过受力分析计算,可得A、B、C、D轴承所承受的径向力F RA、F RB、F RC、F RD:
当该车前进行驶时主动齿轮从X轴方向看为顺时针,旋转方向为左旋,根据表4中的公式代入相应的数据可得F ap=56.22 N,F Rp=21.94 N,F aG=50.11 N,F RG=42.06 N。再将得到的齿轮径向力和轴向力代入到轴承径向力公式中,F A y=-70.34 N,F A z=-56.94 N,得F RA=90.498 N,F B y=14.11 N,F B z=8.086 N,F RB=16.263 N,F C x=21.94 N,F C z=22.05 N,F RC=31.106 N,F D x=10.82 N,F D z=-26.81 N,F RD=28.911 N。
F A x=10.63 N,同理,查得B轴承计算系数Y=2,B轴承派生轴向力F B x=21.75 N,C、D轴承的计算系数Y=1.6,C、D轴承的派生轴向力分别为F C y=35.81 N,F D y=20.4 N。对于主减速器端轴承对,F ap+F B x>F A x,靠近主减速器齿轮A轴承被压紧,F A x=F ap+F B x=77.97 N,同理,F aG+F Cy>F D y,D轴承被压紧,修正F D y=F aG+F C y=85.92 N。
使用模态分析结果,计算桥壳振动响应,求解已知1~2 000 Hz频段的所有结果。将频率范围设成1~2 000 Hz。前面己经算出4个轴承的受力分配情况,在做频响分析时须在4个轴承处施加动态激励,采用谐响应来进行分析。
选择模态叠加法来进行分析,ANSYS workbench求解,得到结果文件。
图13为367 Hz激励下振动位移分布云图,该频率点处恰有第1阶模态,振型表现为后油盖局部扭转,振动响应也相对较大,最大位移达到了0.641 53 mm。另外此处频率相对较高,对声压有一定的影响。振动较大区域还是分布在后油盖处,再一次印证后油盖为主要的声辐射区域,后续的优化部位应该为此处。
3 声学仿真分析
由于驱动桥壳及减速器壳局部结构不同,噪声辐射结果体现在表面及外部的也是不一样的,在桥壳频率响应分析完成后,再对其进行声学仿真分析,可以确定噪声辐射较大部位,以实现后续更精准的结构优化。
3.1 声学仿真分析流程与方法
根据噪声源物理特性,可将噪声分为气体动力噪声、机械噪声和电磁性噪声;根据噪声的频率成分分类,可将噪声分为高频噪声、中频噪声和低频噪声。研究方法主要有边界元法、有限元法和统计能量法。主流声学软件Virtual.Lab Acoustic以边界元为主,具备振动噪声的有限元解算器;Actran以有限元为基础,擅长气动/流动噪声、声振耦合、隔声降噪与声疲劳等问题分析;VA-One以统计能量法为基础,适用于高频振动噪声问题。驱动桥噪声属中低频振动噪声,宜采用边界元法分析,因此选Virtual.Lab Acoustic软 件 的Harmonic/Transient Acoustic FEM/BEM模块。具体分析流程如图14所示。
3.2 声学边界元网格模型的建立
Virtual.Lab Acoustic在进行声学计算时,须保证最小波长内有6个单元,最大的单元边长要小于最高计算频率点处波长的1/6,对于二次单元这个数值就是1/3。假设声音在某流体介质中的传播速度为c,模型计算最高频率点为f max,边界元网格单元长度为L,则单元的长度应满足:
3.3 桥壳辐射噪声结果
图16是频率为367 Hz时桥壳表面声场分布云图。从图中可以看出,后油盖区域声辐射较大,是主要的辐射区域,桥壳的振动响应较大处也是这个区域,再对比桥壳在367 Hz时的模态振型,发现这个区域表现出了主要的振型,说明在367 Hz激励下引起了共振,使桥壳产生了较大的振动响应,然后辐射出较大的噪声,最大声压级达到了72.3 dB。图17是频率为367 Hz时场点声场分布云图,桥弓附近区域声压级较大,另外后油盖区域声压级也较大,与桥壳表面声场分布是一致的,最大声压级为62.5 dB。
4 后桥壳的拓扑优化及分析
结构优化设计是用系统的、目标定向的过程与方法代替传统设计,其目的在于寻求既经济又适用的结构形式,以最少的材料、最低的造价实现结构的最佳性能。
拓扑优化设计是在给定材料品质和设计域内,通过优化设计方法得到既满足约束条件又使目标函数最优的结构布局形式或构件尺寸,是结构优化设计的常用方法。在复杂结构的选型及轻量化设计工作中,拓扑优化比形状、尺寸优化更有价值。拓扑优化的研究领域主要分为连续体拓扑优化和离散结构拓扑优化。
4.1 拓扑分析的流程与方法
拓扑优化设计流程如图18所示,目前连续体拓扑优化方法主要有均质化方法、相对密度法、进化结构优化等主要的拓扑优化方法,后桥壳应采用连续体拓扑优化,并采取相对密度法进行计算。美国Altair公司Hyperworks软件中的Optistruct、德国Fe-design公司的Tosca和ansys workbench 3种软件理论上均具备完成该类分析的功能,但在汽车承载系统拓扑优化领域内采用optisruct较多,此次拓扑优化分析采用Optisruct软件。
4.2 后桥壳基于模态的拓扑优化与改进
前面己经得出,在加速直线60~65 km/h工况下该后桥噪声较大的原因为在367 Hz处啮合冲击与第1阶模态产生共振。在不影响整车模态匹配的前提下,通过适当减质量的方式来降低桥壳模态从而降低振动噪声。驱动桥壳以及后盖是主要的振动辐射区域,这些区域为设计变量空间,以模态为约束响应,以体积最小为目标,利用Optistruct对驱动桥壳及后油盖进行拓扑分析,利用分析结果指导结构优化。
为保证拓扑优化结果的实用性,将该桥壳的第1阶模态上下限分别设为300和350 Hz,既降低了模态,又不至于模态过低引起其他匹配问题。拓扑优化结果如图19所示。从分析结果可知,驱动桥壳可减厚30%,后油盖可减厚15%。
图22为改进后的结构在367 Hz处桥壳表面声场分布图,由于在该频率处没有模态振型,桥壳表面声场最大声压级为60.1 dB,相比原结构减少了12.2 dB。图23为改进后的结构在367 Hz处场点声场分布图,与表面声场分布类似,最大声压级有了明显的降低,降到了54.3 dB。因此,在60~65 km/h加速直线工况下桥壳辐射噪声有了较大的降低。
5 优化后实车NVH测试结果
采用声压级分析测试的方法对车辆加速工况测试,可更全面地分析,且同时测量滑行工况。每个工况在驾驶员右耳、后排乘客右耳两个位置采集数据。在下面4个声压级图中,1#线条代表整车总声压级(overall level),2#线条代表改进前驱动桥声压级,3#、4#线条代表改进后驱动桥声压级,下面是声压级测试情况。
(1)车辆滑行驾驶员右耳处声压级
由图24可知,该工况下,改进前后声压级水平基本一致,未发现差异。这是由于整体噪声较小,非顾客抱怨点,因此没有对其开展优化工作。
(2)车辆滑行后排乘客右耳处声压级
由图25可知,该工况下,改进前后声压级水平基本一致,未发现差异。但是,可以发现在1 823 r/min时,优化前后驱动桥声压级均存在明显峰值,应为后桥与其他系统匹配问题,较大可能是与传动轴匹配问题,虽然滑行工况使用不多,一般消费者难以发现,但是产品性能仍需提高,是后续改善点。
(3)车辆加速驾驶员右耳处声压级
(4)车辆加速后排乘客右耳处声压级
从图27可知,该工况下,改进前驱动桥在2 000和3 000 r/min时,均存在明显峰值,最大峰值接近70 dB,消费者抱怨是正常的,优化后的两台驱动桥,在2 000和3 000 r/min时,均不存在峰值,曲线比较平缓,噪声水平降低大于10 dB,说明此次驱动桥结构优化工作是成功的。
6 结论
本文中建立了系统的、理论与实际相结合的汽车驱动桥NVH性能分析与优化的流程与方法,并对其有效性进行了验证,研究工作取得以下研究成果:(1)所建立的汽车驱动桥NVH性能分析与优化的流程与方法是准确、可行的;(2)通过有限元、振动响应、声学仿真分析,并恰当选取相应的算法及软件,可以精准找到驱动桥具体需要优化的部位;(3)采用Optisruct软件,以模态为约束,以体积最小为目标的方法局部优化驱动桥,可实现以减质量的方式来降低模态从而降低噪声。此次研究,仅局限在驱动桥自身的激励,对外部激励源对驱动桥NVH影响的研究是后续研究的方向。
-
汽车测试网V课堂
-
微信公众号
-
汽车测试网手机站
最新资讯
-
NVIDIA 发布 2025 财年第三季度财务报告
2024-11-21 13:30
-
Mack卡车为买家推出创新的虚拟现场探索体验
2024-11-21 13:29
-
氢燃料电池卡车从1到100要多长时间?戴姆勒
2024-11-21 13:28
-
聚焦消费者用车极限环境,2024中国汽研汽车
2024-11-21 13:21
-
新能源汽车高寒环境可靠性行驶试验研究
2024-11-21 13:19