謝華林,龔運息,堯瀟雪,黃大偉
(1.廣西科技大學 機械與交通工程學院,廣西 柳州 545616;2.佛吉亞(柳州)排氣控制技術(shù)有限公司,廣西 柳州 545000)
排氣系統(tǒng)是汽車振動分析、噪聲分析的重要組成部分之一。它的熱端(三元催化氧化器)與發(fā)動機相連,冷端(波紋管、消聲器、排氣尾管等)通過吊鉤與車身連接。當發(fā)動機工作時,所產(chǎn)生的振動會通過排氣系統(tǒng)傳遞到吊耳掛鉤,然后傳遞到汽車底盤,從而影響汽車的乘坐舒適性。因此,吊鉤位置布置是否合理,是對于排氣系統(tǒng)能否有效地降低車身的振動,降低車內(nèi)振動和噪聲水平,以及確保進行后續(xù)整車NVH的研究的重要前提。目前,國內(nèi)外學者主要采用平均驅(qū)動自由度位移法[1-4],也就是(ADDOFD)法對排氣系統(tǒng)吊鉤位置進行設計,但是這種方法在數(shù)據(jù)處理過程中較為繁瑣,其需要關(guān)注各潛在節(jié)點在不同階振型下的位移輸出與每一階系統(tǒng)振型。而均方根值(RMS)法相對與ADDOFD法在數(shù)據(jù)處理上更加簡單快速,只需將設計吊點位置函數(shù)進行積分求均方根值,找到振動能量最小值點即可。
由此,基于上述兩種方法對某車型排氣系統(tǒng)懸掛位置進行了設計與對比分析。通過UG和Hyperworks軟件建立有限元模型并進行驗證,基于ADDOFD法與均方根值法得到了兩種方案的吊鉤懸掛位置,結(jié)合排氣系統(tǒng)與整車底盤的安裝情況和仿真驗證(模態(tài)分析、靜力分析以及動力響應分析),確定了吊鉤的最終位置,且較為合理地預測了此排氣系統(tǒng)的性能,從而為之后整車NVH性能的改善提供了一定的理論依據(jù)。
通過UG建立排氣系統(tǒng)三維實體模型,然后導入Hypermesh軟件中。為了不影響其模型求解精度,減少計算時間,作出了以下結(jié)構(gòu)簡化處理[5]。其中,對結(jié)構(gòu)厚度與長度相差太大的組件采取抽中面(midsurface)的方法,將排氣系統(tǒng)連接管、消聲器殼體、消聲器內(nèi)部進、排氣管、隔熱板等簡化為三角形或四邊形的殼體單元;消聲器內(nèi)部穿孔管、隔板上的小孔不予考慮;波紋管使用一個彈簧單元和兩個質(zhì)量點來代替,兩端用RBE2剛性單元與連接管道相連。材料為SH409L,彈性模量2.02×105MPa,泊松比0.3,密度為7.72×10-9t/mm3。
設定單元尺寸為5 mm,對該排氣系統(tǒng)進行網(wǎng)格劃分,法蘭之間的連接采用螺栓緊固,可以看成是剛性連接,用RBE2進行模擬連接[6],剩下的其余部分組件結(jié)構(gòu)之間也使用RBE2進行連接,得到有限元模型如圖1所示。
圖1 排氣系統(tǒng)有限元模型
假設單點激勵,根據(jù)多自由度系統(tǒng)模態(tài)分析理論[7],響應點l與激勵點p之間的頻響函數(shù)為:
(1)
式中:φlr為第l個測點、第r個模態(tài)振型系統(tǒng);Mr和ξr分別為模態(tài)質(zhì)量和模態(tài)阻尼比。
如果激勵力的頻率為ωr,則近似地有:
(2)
對于線性系統(tǒng), 位移響應的幅值和頻響函數(shù)的幅值成正比, 即
(3)
進一步假設振型以質(zhì)量矩陣歸一化,各階模態(tài)阻尼近似相等,則:
(4)
為預測某個自由度在一般激勵情況下 (在某個頻率范圍所有的N個模態(tài)均被激發(fā))位移響應的相對大小, 定義第j個自由度的ADDOFD為:
(5)
在波紋管末端,沿X軸方向每隔50mm選取一個潛在吊掛點位置進行編號,然后利用平均驅(qū)動自由度位移法計算,得到排氣系統(tǒng)0~200 Hz以內(nèi)的ADDOFD值,如圖2所示。
圖2 排氣系統(tǒng)潛在懸掛點的平均驅(qū)動自由度位移結(jié)果
結(jié)合排氣系統(tǒng)的結(jié)構(gòu)和平均驅(qū)動自由度位移理論,得到排氣系統(tǒng)潛在懸掛點編號,如圖2中位置1~5所示的點[8]。由于位置1處于彎管中段,考慮車身結(jié)構(gòu)以及安裝情況,所以位移1不作為潛在選取點。因此,最后確定的潛在位置點有4個,分別為圖中的位置2、3、4、5。由此可以得到掛鉤的設計方案如圖3所示。
圖3 吊鉤位置設計方案1
均方根值,又被稱為方根均值或有效值。在一段時間T內(nèi),對于隨機信號f(t),將其所有值平方求和然后再開平方[9],即:
(6)
自相關(guān)函數(shù)Rf反映與隨機變量在時域范圍內(nèi)特性相關(guān)的信息:
(7)
式中:Sf為功率譜密度函數(shù),反映時域范圍內(nèi)的相關(guān)信息;ω為信號頻率;τ為時間差。
當τ=0時,式(7)可改寫為:
(8)
式(6)、(7)和(8)表明,時域范圍內(nèi)的均方根值(RMS)與頻域范圍內(nèi)頻響曲線下方區(qū)域的面積相等,對頻響曲線在頻域范圍內(nèi)進行積分即可得到時域范圍內(nèi)的均方根值(RMS)。對于振動信號,此均方根值(RMS)可表征系統(tǒng)振動能量的大小[10]。
沿波紋管末端,在X軸方向每隔50 mm選取一個潛在吊掛點位置進行編號,約束排氣系統(tǒng)前端法蘭除Z向外的所有自由度,并在前端法蘭處施加范圍為0~200 Hz的Z向激勵[11],得到69個點的Z向頻率響應函數(shù)。然后通過調(diào)用RMS函數(shù)對69條曲線在頻域范圍內(nèi)進行積分,從而得到各點的振動能量大小如圖4所示。
圖4 潛在吊鉤位置點的均方根值
由方均根值理論我們可以知道,選擇振動能量較小的位置放置吊鉤比較合適,即圖4中位于波谷及波谷附近的點,如17號、28號、40號、56號以及69號點。由于17號點與28號點位置靠的比較近,且它們都屬于前消聲器總成部分,而17號點位置相對于28號點比較靠前,再結(jié)合整體受力情況,很明顯可以看出17號點和40號點比28號點和40號點的組合更為合適。最終擬定的設計方案如圖5所示。
圖5 吊鉤位置設計方案2
(1) 約束模態(tài)分析
對上述兩種方案排氣系統(tǒng)進行約束模態(tài)分析,約束邊界條件:約束三元催化氧化裝置法蘭與發(fā)動機懸置點,約束所有車身側(cè)吊鉤。吊鉤使用六面體實體單元進行計算模型的建立,排氣系統(tǒng)懸掛膠采用彈簧連接[12],排氣系統(tǒng)懸掛的剛度在x、y和z方向分別為3 N/mm、3 N/mm和12 N/mm。求解頻率范圍設置為20~200 Hz, 利用Hypermesh軟件中的OptiStruct求解器進行約束模態(tài)求解計算,得到兩種方案的約束模態(tài)固有頻率如表1所列。
表1 兩種方案約束模態(tài)分析結(jié)果對比
(2) 靜力分析
約束三元催化氧化裝置法蘭和排氣系統(tǒng)吊鉤位置橡膠吊耳處與車身相連的吊鉤端,然后加載重力加速度進行靜力分析。從而得到如圖6兩種方案下的各個被動側(cè)吊鉤支反力與圖7兩種方案下的排氣系統(tǒng)整體位移圖,以及圖8橡膠吊耳的最大位移值。
圖6 兩種方案下的各個被動側(cè)吊鉤支反力
圖7 兩種方案下的排氣系統(tǒng)整體位移圖
通過對比兩種方案的約束模態(tài)以及重力工況下靜強度計算情況,得出以下結(jié)論。
① 約束模態(tài):在頻率20~200 Hz范圍下,兩種方案的各階模態(tài)值與其相對應的振型圖相差不大。
② 靜力分析計算:如圖7、8可以看出,方案二排氣系統(tǒng)的整體位移量比方案一小,且方案2吊鉤位置的支反力分布較均勻,最大靜應力較小。
③ 兩種方案下的四組橡膠吊耳最大位移均小于3 mm,排氣系統(tǒng)整體載荷分布均勻,兩種設計方案均滿足重力載荷下的相關(guān)要求。
圖8 兩種方案下的橡膠吊耳位移
動力響應分析的目的是分析傳遞到車體上的力,簡而言之就是分析橡膠吊耳上的傳遞力,對于普通的轎車而言,目標一般不大于10 N[13]。由于受到掛鉤受力大小和空間位置的限制,并不是所有掛鉤都能處于位移最小值點,所以我們選擇兩種方案中位置一致的掛鉤1、掛鉤2和掛鉤4作為研究對象。參考企業(yè)標準,在發(fā)動機質(zhì)心加載1×105 N·mm的扭矩,激勵頻率范圍在0~120 Hz,得到兩種方案下橡膠吊耳的傳遞力如下圖10所示。結(jié)果表明,兩種方案下的吊鉤傳遞力都小于10 N,滿足設計要求。
圖9 兩種方案下的吊耳傳遞力
基于ADDOFD法與均方根值法理論基礎(chǔ),通過CAE仿真得到了兩種不同的吊鉤位置設計方案,并進行了一系列仿真驗證,得出以下結(jié)論:
(1) 從兩種方法確定的吊鉤位置計算驗證得到:約束模態(tài)分析結(jié)果相差不大;靜力分析方案2的系統(tǒng)最大位移較小且掛鉤的支反力較均勻,橡膠吊耳的最大位移量都小于3mm;動力響應分析吊耳傳遞力都小于10mm;綜合來看方案2比方案1更好。
(2) 對于排氣系統(tǒng)吊鉤位置的設計,均方根值法相較于傳統(tǒng)的平均驅(qū)動自由度法在數(shù)據(jù)處理上速度更快,位置點確定較容易,且設計驗證都滿足排氣系統(tǒng)相關(guān)設計要求。