譚繼錦,戚東元
(合肥工業(yè)大學 機械與汽車工程學院,合肥 230009)
轎車帶內飾車身模態(tài)試驗分析
譚繼錦,戚東元
(合肥工業(yè)大學 機械與汽車工程學院,合肥 230009)
針對帶內飾車身模態(tài)參數(shù)難以識別的問題,在結構線性度檢驗的基礎上,分別采用頻域識別和時域識別兩種方法對某轎車帶內飾車身進行模態(tài)試驗。介紹了試驗方案與信號采集技術,重點對識別方法進行比較研究,其中時域隨機子空間法的較好應用效果為車身模態(tài)試驗提供新思路。最后通過有限元分析與試驗結果對比,驗證了內飾車身建模分析的準確性,所確定的內飾車身模態(tài)參數(shù)為車身NVH性能分析奠定基礎。
振動與波;帶內飾車身;模態(tài)試驗;頻域識別;時域識別;有限元分析
內飾車身(Trimmed Body)與白車身(Body In White)一樣,在車身性能開發(fā)中都有具體的模態(tài)頻率目標,其中白車身模態(tài)分析應用較為成熟,不論是試驗還是仿真均可識別出較為準確的模態(tài)參數(shù)。相對于白車身,內飾車身的結構性能更接近于整車,也是整車模態(tài)頻率的基礎,因而研究內飾車身的動態(tài)性能十分必要。
閉合件、座椅、轉向系統(tǒng)和內飾通過鉸鏈、密封膠、結構膠和卡扣等與白車身連接,使得車身非線性因素增加,模態(tài)參數(shù)密集,阻尼增大,局部與整體振型耦合較多,大大增加了識別難度。為解決此問題,文獻[1]提出取車身左右兩側若干對稱點,經過VTF分析,通過合成加速度和差值加速度均值曲線來識別內飾車身模態(tài)參數(shù)。文獻[2]分別采用最小二乘復指數(shù)法(LSCE)和最小二乘復頻域法(LSFD)對內飾車身模態(tài)試驗進行參數(shù)識別,并比較了兩種方法的優(yōu)劣。文獻[3]提出一種新的整體平順算法(GST),用于分析多輸入多輸出系統(tǒng)所得的頻響函數(shù)曲線,使得對于大阻尼高耦合的結構模態(tài)參數(shù)得以識別。
文中將從頻域和時域兩種識別方法對帶內飾車身模態(tài)試驗進行研究和對比,提出內飾車身的試驗方案,闡述解決模態(tài)參數(shù)識別困難的辦法,并將有限元分析與試驗進行對比。
1.1 結構線性度檢驗
模態(tài)識別以系統(tǒng)線性為前提,任何結構都必須在滿足一定線性度的前提下才可以進行模態(tài)分析。內飾車身非線性度加大,試驗之前必須對其進行結構線性度檢驗,以確保此結構模態(tài)測試的可靠性。設置不同等級的激振力,以同一測點頻響函數(shù)的重合度進行檢驗。圖1為車身某測點在4 V、8 V和10V激勵下的頻響曲線,可以看出100 Hz內,整體重合度良好,車身保持一定的線性度。僅在80 Hz之后出現(xiàn)微量偏差,后續(xù)研究表明此頻段由于內外飾件振動復雜密集,使得車身非線性增大,但其仍在可識別范圍內。
圖1 帶內飾車身線性度檢驗
1.2 試驗方案確定
試驗方案的確定要結合內飾車身的實際,具體包括以下幾個方面:
(1)測點選擇。測點布置應能反映車身整體輪廓,同時要注意避開局部剛度較弱的區(qū)域。內飾車身由于有門蓋系統(tǒng)和內飾附件等,測點布置較為困難。采用兩種測點方案,對于頻域法采用密集型布置,車身共布置249個測點,密集測點布置可以在識別車身整體振型的同時,識別出局部振型。為使測點布置均勻,先對車身表面進行100 mm×100 mm網格劃分,再選取測點,測點方案如圖2所示。
圖2 頻域法測點方案
(2)支撐方式。模態(tài)試驗結構安置需盡量接近自由狀態(tài),白車身多采用橡皮繩懸吊方式,內飾車身質量較大,采用空氣彈簧支撐方式,便于安置車身與安放傳感器。
(3)激勵方式。模態(tài)試驗激勵有單點激勵和多點激勵兩種,采用何種方法取決于試件相對于整體的激振難度,內飾車身單點激勵即可獲得任意點的響應,且響應幅度足夠大,故采用單點激勵多點識振的方式,激振點設于前保險杠右側(如圖2所示)。
對于激勵信號的選擇,內飾車身力錘敲擊能量不足,正弦激勵效果較差,經比較猝發(fā)隨機為最佳激勵方式。激勵力大小應具有足夠的能量激發(fā)內飾車身結構,但是不能超出結構可接受的線性區(qū)域內,經調試定為8 V。
(4)試驗頻段。相對于白車身,內飾的添加會使車身整體一彎一扭頻率有所下降,綜合考慮,試驗頻段選取0~100 Hz。試驗前對加速度傳感器在頻段內進行標定。
1.3 測試系統(tǒng)的組建與數(shù)據的采集
試驗測試系統(tǒng)采用DHDAS動態(tài)信號采集分析系統(tǒng),系統(tǒng)框圖見圖3。
圖3 試驗測試系統(tǒng)框圖
試驗系統(tǒng)建立后,調試測試參數(shù),選擇最優(yōu)參數(shù)設置。然后需對其進行可靠性驗證,常用的驗證方法為頻響函數(shù)的互異性、重復性和相干性檢驗等?;ギ愋赃x擇前保險杠左右兩側對稱處激勵,檢查頻響函數(shù)重合度。每組測點逐一檢查相干性,圖4為某組4號通道的相干函數(shù)曲線,從圖中可以看出盡管在15 Hz之前有部分泄露,但在全頻段內基本保持為1,表明頻響函數(shù)具有較高性噪比,系統(tǒng)相干性在合理范圍內。
圖4 某組4號通道相干函數(shù)
249個測點分批進行測量,每組10個測點。測試分析點數(shù)為2 048,線性平均次數(shù)為30,重疊率為30%,輸入通道加力窗,輸出通道加指數(shù)窗。測試時實時觀察記錄頻響函數(shù)曲線和相干函數(shù)曲線質量,確保數(shù)據可行性,如果出現(xiàn)某一個通道或是某組通道異常,需停止試驗進行通道調試檢查。每組試驗均需要檢查傳感器接線狀況是否良好,在放置傳感器時,為減小壓電傳感器的接地回路噪聲,做好系統(tǒng)接地與絕緣工作[4]。所有通道測試前清零平衡。
1.4 車身模態(tài)參數(shù)識別
基于頻域的模態(tài)識別方法最大的優(yōu)點在于使用平均技術,很大程度抑制了噪聲,使得定階問題易于解決?,F(xiàn)有方法主要有最小二乘復頻域法(LSCF)、峰值拾取法、多項式擬合法、頻域分解法等,其中最小二乘復頻域法在應對密集型模態(tài)時,即便頻響函數(shù)受到部分噪聲干擾,依然可以建立清晰的穩(wěn)態(tài)圖,識別效果較優(yōu)。經比較采用LSCF方法作為內飾車身頻域模態(tài)參數(shù)識別方法,多項式擬合法作為參考,圖5為LSCF法穩(wěn)態(tài)圖。
圖5 基于LSCF法的穩(wěn)態(tài)圖極點識別
以穩(wěn)態(tài)圖為依據,選擇模態(tài)階數(shù),穩(wěn)態(tài)圖中S表示假定極點增加時在給定精度內頻響函數(shù)在頻率、阻尼和模態(tài)參與因子處(V是在頻率和模態(tài)參與因子處)保持穩(wěn)定狀態(tài)[5]。以S縱向較為一致為定階條件。表1給出了提取出的19階模態(tài),圖5-圖6為車身垂向1階彎曲和1階扭轉振型。
圖6 車身整體1階扭轉
表1 LSCF法模態(tài)參數(shù)識別結果
對于門蓋系統(tǒng)等組件局部振型難以辨識的情況,為避免漏階,可將組件的測試數(shù)據進行單獨處理,同時將組件模型單獨創(chuàng)建。圖7為發(fā)動機艙蓋的1階扭轉和1階彎曲振型,單獨處理后振型清晰。
圖7 車身垂向1階彎曲
圖8 發(fā)動機艙蓋1階扭轉與1階彎曲
為確保識別結果的可信度,進行振型相關性驗證,圖8為MAC值驗證結果,可以看出提取的19階模態(tài)具有較好的正交性。
從頻域法振型識別結果來看,車身整體振型與局部振型相互耦合,并會出現(xiàn)局部異常測點,對振型識別產生很大干擾,甚至造成誤判。針對此問題,以下對內飾車身進一步展開時域法識別模態(tài)試驗研究。
基于頻域的模態(tài)識別方法是機械類結構模態(tài)分析中的主流方法,對于材料單一的簡單結構,頻域法具有較大優(yōu)勢。而內飾車身具有一定非線性度,耦合振型較多,振型的識別較為困難,為此借用基于環(huán)境激勵的工作模態(tài)[6]的試驗方法對內飾車身進行時域法參數(shù)識別,并將試驗結果與頻域法結果進行對比。
2.1 試驗方案
基于時域法的模態(tài)試驗測試系統(tǒng)與前述相同,鑒于篇幅,僅敘述試驗方案的差異性。區(qū)別于環(huán)境激勵,車身位于試驗臺架上,對其進行激勵使得車身振動,激勵信號采用隨機白噪聲。需注意的是,采集系統(tǒng)不對輸入信號進行采集,關閉力采集通道。試驗重新布置測點方案,共37個測點,如圖9所示。
僅在車身骨架上布置測點,目的在于減少局部振型的干擾,著重識別車身整體模態(tài)參數(shù)。測試以1號測點為參考點,采用移動測試法,保持1號測點不動,分4組完成。
圖9 LSCF法振型有效性驗證
圖10 時域法測點方案
2.2 基于隨機子空間法的車身模態(tài)參數(shù)識別
時域識別方法也有多種,目前主要有隨機子空間法(SSI)、最小二乘復指數(shù)法(LSCE)、特征系統(tǒng)實現(xiàn)法(ERA)、Ibrahim時域法(ITD)等。在工程應用中隨機子空間法效果較優(yōu),其基本思路如下:首先,在白噪聲激勵的前提下,建立線性振動系統(tǒng)離散狀態(tài)下的空間模型[7]。
式中xk∈Rn為離散時間狀態(tài)向量,yk∈Rl×1為輸出向量,A∈Rn×n為離散狀態(tài)矩陣,C∈Rl×n為輸出矩陣,E為數(shù)學期望符號,δpq為kronecker矩陣,wk∈Rn為過程噪聲,vk∈Rl為測量矩陣,兩種噪聲均假定為白噪聲,故
將測量數(shù)據構造為Hankel矩陣H,按矩陣行數(shù)的不同將H分為“過去”與“將來”兩個子矩陣,即QR分解:
然后對加權投影矩陣做奇異值分解,通過系統(tǒng)矩陣識別模態(tài)參數(shù)[8]。應用此方法識別結果如表2所示,對應振型見圖10-圖15。
表2 基于隨機子空間法的模態(tài)參數(shù)識別結果
圖11 車身1階扭轉
圖12 車身底板局部扭轉
圖13 車身垂向1階彎曲
圖14 頂蓋前部局部振型
圖15 頂蓋后部局部振型
2.3 結果對比與分析
從表2可以看出,隨機子空間法識別與頻域法識別頻率誤差在5%以內,振型一致,表明基于時域辨識方法可行有效。比較而言,試驗方案更為簡單,無需采集激勵信號,響應信號也無需進行傅里葉變換,直接作用于加速度信號,避免了變換誤差。并且由于測點布置于車身骨架,車身整體振型更加清晰。
利用Hypermesh搭建帶內飾車身有限元模型,見圖16。其中車身骨架主要由板單元與焊點組成,網格大小選擇8 mm,忽略尺寸較小的工藝孔,重要區(qū)域適當加密,采用鋼質材料,焊點采用ACM焊接模擬,螺栓采用梁單元,焊縫采用剛性連接。
圖16 帶內飾車身有限元模型
門蓋系統(tǒng)主要有板件、玻璃和膠粘連接組成。采用實體單元模擬膠體,再采用CWELD單元連接膠體與板件。轉向系、座椅主要抽取其鋼質材料進行建模,其余內外飾組件多采取附加質量方式模擬。
模型搭建完成后共有1 729 565個單元、1 864 775個節(jié)點、5 703個焊點,三角形單元有75 498個,占總數(shù)的4.3%。模型總重725 kg,試驗車身實際重741 kg,考慮底漆與小件質量的存在,誤差在合理范圍內。
將內飾車身計算結果與頻域法識別結果進行對比,表3列出了主要階次的模態(tài)對比。
表3 計算結果與試驗結果對比
從表3可知,有限元仿真結果中大部分模態(tài)頻率與試驗結果對應,且誤差在許可范圍之內,表明內飾車身整體模型達到要求。需要注意的是車身某些局部振型,如頂蓋在115 Hz之后才出現(xiàn)較大振動,與試驗結果相差較大,需分析原因調整模型。
帶內飾車身屬于模態(tài)密集的復雜振動系統(tǒng),內飾的添加、材料的豐富、以及連接方式的多樣使得車身的非線性度增加。為避開相互干擾,可將車身整體與各部件分開識別,僅以車身骨架布置測點可較容易識別車身整體模態(tài)。最小二乘復頻域法在內飾車身模態(tài)識別時穩(wěn)態(tài)圖較為清晰,但是振型較難甄別。工作模態(tài)試驗方法用于內飾車身模態(tài)參數(shù)識別,提供了新思路與相互驗證手段。有限元計算結果與試驗結果對比,驗證了內飾車身建模分析的準確性。內飾車身模態(tài)數(shù)據為其動力學評價與整車NVH性能控制奠定了基礎。
鑒于內飾車身模態(tài)測試難度較大,識別振型較難,且仍然與實際整車狀態(tài)有所差別,建議可以從整車工作模態(tài)角度識別車身振型與模態(tài)參數(shù),通過路試,直接面對車身NVH問題,也不失為一種解決辦法。研究整車-內飾車身-白車身模態(tài)關聯(lián)技術,這方面的對比研究正在進行。
[1]藍浩倫,李飛,喬鑫,等.汽車Trimmed Body有限元模態(tài)分析中彎曲和扭轉模態(tài)的識別[J].汽車技術,2015(4):26-29.
[2]SCZIBOR V,ALVES P S L,BRINGHENTI I.Modal parameter estimation on automotive development[J].SAE Technical Paper,2012.
[3]RUOTOLO R,STORER D M,GATTI S.Experimental modal analysis to represent structural components in vehicle models[J].StrategicManagementJournal, 2001,31(13):1498-1526.
[4]劉煥光,陳朝陽,譚繼錦,等.商務車白車身試驗模態(tài)分析研究[J].合肥工業(yè)大學學報,2007,30(12):1615-1619.
[5]謝小平,韓旭.基于PolyMAX方法的某轎車白車身實驗模態(tài)分析[J].汽車工程,2009,31(5):440-447.
[6]M?LLER N,GADE S.Application of operational modal analysis on cars[J].SAE Transactions,2003,112(6): 1970-1977.
[7]肖祥,任偉新.實時工作模態(tài)參數(shù)數(shù)據驅動隨機子空間識別[J].振動與沖擊,2009,28(8):148-153.
[8]辛峻峰,王樹青,劉福順.數(shù)據驅動與協(xié)方差驅動隨機子空間法差異性分析[J].振動與沖擊,2013,32(9):1-4.
Analysis and Research of Modal Test for a Trimmed Car Body
TAN Ji-jin,QI Dong-yuan
(Collage of Mechanical andAutomotive Engineering,Hefei University of Technology, Hefei 230009,China)
In order to solve the problem that modals of the trimmed car body are difficult to identify,both frequencydomain and time-domain identifications are used for the modal test of the trimmed body of a car on the basis of the structure’s linearity inspection.The test scheme and signal acquisition technique are introduced and the identification methods are comparatively studied.Of these methods,the successful application of the stochastic subspace identification provides a new way for the body modal test.Finally,the comparison of results from the test and FEA verifies the accuracy of the trimmed body modeling analysis.The determined modals of the trimmed body lay a foundation for the analysis of NVH performance.
vibration and wave;trimmed body;modal test;frequency domain identification;time domain identification;FEA
U467;O241.82
:A
:10.3969/j.issn.1006-1335.2016.06.039
1006-1355(2016)06-0197-05+205
2016-07-18
戚東元(1991-),男,江蘇省連云港市人,碩士研究生,研究方向為振動噪聲控制與整車性能分析。
譚繼錦,男,碩士生導師。E-mail:jsdh_qdy@126.com