石 錦 杜登惠 蔣 晨 岳亮亮
(泛亞汽車技術(shù)有限公司)
汽車保險杠系統(tǒng)的整體模態(tài)是其前期開發(fā)階段需要關(guān)注的重要性能指標。若模態(tài)指標不滿足要求,將在整車疲勞耐久試驗或客戶使用過程中出現(xiàn)保險杠系統(tǒng)結(jié)構(gòu)損壞等問題。
在汽車保險杠產(chǎn)品設(shè)計研發(fā)階段,因沒有物理樣件,所以通常采用CAE虛擬仿真來評估保險杠系統(tǒng)模態(tài),預測其結(jié)構(gòu)的動態(tài)特性,因而較高CAE仿真結(jié)果的置信度成為保險杠開發(fā)能力的體現(xiàn)。本文通過對虛擬分析結(jié)果與物理試驗結(jié)果進行對比分析,合理改進CAE分析方法,進一步提高保險杠系統(tǒng)模態(tài)仿真結(jié)果的置信度。
保險杠系統(tǒng)的模態(tài)分析方法主要有理論模態(tài)分析方法和試驗?zāi)B(tài)分析方法。
在實際工程中,通常采用數(shù)值方法進行理論模態(tài)分析,由于有限元法具有很強的實用性和通用性,區(qū)域離散的有限元法在眾多數(shù)值方法中最受歡迎[1],因此,將采用有限元法對汽車前保險杠進行理論模態(tài)分析,求解基本原理如下。
對于具有多個自由度的系統(tǒng)的強迫振動方程為[2]:
式中,[M]為質(zhì)量矩陣;[C]為阻尼矩陣;[K]為剛度矩陣;{x}為節(jié)點位移列陣;{f}為結(jié)構(gòu)載荷列陣。
由于需要分析結(jié)構(gòu)的固有特征,所以取{f}=0;同時因結(jié)構(gòu)的阻尼很小,對固有頻率和振型的影響可忽略不計,則無阻尼自由振動方程為:
式(2)為常系數(shù)齊次微分方程,設(shè)其解為:
式中,xM為振幅列矢量;j2=-1;t為時間;w為振動固有頻率,令w2=λ。
將式(3)代入式(2),并令其有非零解的條件系數(shù)行列式必須等于零,則
由式(4)可解得特征值所對應(yīng)的第i階特征矢量 xMi,即對應(yīng)于各階固有頻率 w1,w2,,…,wn,則有各階固有振型 xM1,xM2,…,xMn。
模態(tài)試驗測試系統(tǒng)包括激勵設(shè)備、傳感系統(tǒng)和分析設(shè)備三部分。通常采用力錘作為激振,以單點激振、多點拾振來采集數(shù)據(jù)。采集的數(shù)據(jù)經(jīng)過專業(yè)軟件處理即可得到模態(tài)振型和頻率。前保險杠系統(tǒng)的模態(tài)測試在整車上進行,其試驗流程見圖 1。
采用hypermesh軟件截取半車模型進行有限元建模,包括前保險杠系統(tǒng)、部分白車身、副車架、翼子板和大燈等,如圖2所示。對于前保險杠系統(tǒng)的各零件,用網(wǎng)格準確表現(xiàn)其幾何特征和安裝孔位;由于塑料零件可能存在厚度變化的情況,根據(jù)不同幾何數(shù)模分組并賦予不同的厚度。
前保險杠系統(tǒng)在車身側(cè)的安裝點同樣是建模的關(guān)鍵。對于螺栓、螺釘、推釘?shù)冗B接方式,各安裝孔采用了兩圈四邊形單元模擬,以確保內(nèi)圈的外徑與墊片尺寸相同,并采用MPC單元連接內(nèi)圈節(jié)點;對于塑料卡扣,根據(jù)搭接位置和面積選取節(jié)點并采用MPC單元進行連接;采用ACM模擬點焊,輸入實際焊點的直徑;采用adhesives模擬涂膠;采用共節(jié)點的方式模擬T型部位的超聲波焊接。
利用有限元軟件Nastran計算前保險杠系統(tǒng)模態(tài),得到第1階頻率為22.1 Hz,振型為上下方向的彎曲振動,如圖3所示。由于前保險杠近似懸臂系統(tǒng)安裝在車身前端,加之集成了許多功能件,具有一定質(zhì)量,因此上下方向的彎曲振動是其主要振型之一。
根據(jù)有限元分析結(jié)果,將三向加速度傳感器布置于非低階模態(tài)振型節(jié)點位置,并確保傳感器的位置和數(shù)量可描述保險杠系統(tǒng)的大致外形,以測量和記錄測點在整車坐標系中的位置。本文所研究的前保險杠系統(tǒng)中布置了30個傳感器,包括車身側(cè)4個測量點、副車架上3個測量點、蒙皮上16個測量點、下腹板上6個測量點和小腿保護桿上1個測量點,見圖4和圖5。
考慮到保險杠系統(tǒng)總成尺寸及特性,利用Maxwell互異性原理,采用錘擊法輸入能量,同時通過調(diào)整合適的輸入能量確定系統(tǒng)線性范圍,在系統(tǒng)線性范圍內(nèi)測得試驗結(jié)果。
首次測量得到的原始狀態(tài)的前保險杠系統(tǒng)模態(tài)試驗頻響結(jié)果如圖6所示,可知試驗的第1階模態(tài)為28.9 Hz,為彎曲振型。雖然第1階模態(tài)振型與有限元分析結(jié)果一致,但由于頻率與有限元分析結(jié)果相差達6.8 Hz,因此需要進行試驗與有限元模型對比研究。
通過對試驗車狀態(tài)、整車CAD數(shù)模及有限元模型進行檢測發(fā)現(xiàn),在實車中,吸能塊在垂直和水平方向被緊密地擠壓在前保險杠蒙皮與防撞梁之間,小腿保護桿前端也與前保險杠緊密接觸;在CAD數(shù)模中,吸能塊與防撞梁之間及小腿保護桿與前保險杠蒙皮之間存在設(shè)計間隙,見圖7。
在設(shè)計前保險杠系統(tǒng)結(jié)構(gòu)時,為保證部件之間的配合會留有一定間隙,以便吸收制造偏差,避免在裝配過程中發(fā)生干涉,影響功能實現(xiàn)或造成安全隱患[3],通常吸能塊與防撞梁之間的設(shè)計間隙約為5 mm。而試驗車上的吸能塊由于制造偏差的緣故吸收了間隙,與防撞梁密實壓緊。小腿保護桿與吸能塊的情況相同,在制造偏差的影響下與蒙皮發(fā)生了接觸。然而由于設(shè)計間隙的存在,在初始的有限元模型中并未考慮上述接觸關(guān)系。
為了在試驗中考查吸能塊和小腿保護桿對前保險杠系統(tǒng)模態(tài)的貢獻量,設(shè)計了2組試驗:第1組試驗去除吸能塊,第2組試驗去除吸能塊和小腿保護桿。試驗結(jié)果見表1、圖8和圖9。
表1 前保險杠系統(tǒng)模態(tài)試驗結(jié)果 Hz
由試驗結(jié)果可知,去除吸能塊后,第1階模態(tài)頻率由28.9 Hz下降至24.6 Hz,其在該車型中對前保險杠系統(tǒng)模態(tài)的貢獻量為4.3 Hz,約占原始模型第1階頻率的15%;而去除吸能塊和小腿保護桿后,第1階頻率為24.9 Hz,比僅拆除吸能塊升高0.3 Hz,說明小腿保護桿對模態(tài)的貢獻量為負。
為在有限元分析中考查吸能塊和小腿保護桿對前保險杠系統(tǒng)模態(tài)的貢獻量,采取去除吸能塊與同時去除吸能塊和小腿保護桿,并分別在考慮接觸和不考慮接觸2種情況下進行分析計算,計算結(jié)果見表2。
表2 前保險杠系統(tǒng)模態(tài)分析計算結(jié)果 Hz
由表2可知,試驗結(jié)果與有限元分析結(jié)果的第1階頻率的振型均一致,為彎曲模態(tài)。當考慮接觸時,頻率的計算結(jié)果與試驗結(jié)果吻合度較高;不考慮接觸時,計算結(jié)果與試驗結(jié)果相差較大。
當原始模型中未考慮吸能塊和小腿保護桿與周邊零件的接觸時,模態(tài)計算值為22.1 Hz;當考慮接觸后,模態(tài)計算值為28.3 Hz,與試驗結(jié)果僅相差0.6 Hz。進一步分析吸能塊和小腿保護桿對模態(tài)的貢獻量,先拆除吸能塊,并考慮小腿保護桿接觸蒙皮,則計算結(jié)果為24.8 Hz,與試驗結(jié)果僅相差0.2 Hz;然后同時拆除吸能塊和小腿保護桿,計算結(jié)果為24.1 Hz,與試驗結(jié)果相差0.8 Hz。
由上述可知,吸能塊的接觸導致第1階模態(tài)頻率上升3.5 Hz,貢獻量達12.4%;小腿保護桿的接觸使得第1階模態(tài)頻率上升2.4 Hz,貢獻量達8.5%。由此表明吸能塊和小腿保護桿在前保險杠系統(tǒng)中的接觸關(guān)系對提高模態(tài)頻率起了很大作用。
通過對某車型前保險杠系統(tǒng)模態(tài)分析結(jié)果和試驗結(jié)果進行相關(guān)性分析,研究了接觸在前保險杠系統(tǒng)模態(tài)分析中的貢獻量,結(jié)果表明:
a.汽車前保險杠中部件之間的接觸關(guān)系對模態(tài)結(jié)果有很大影響。在有限元模型中適當考慮吸能塊和小腿保護桿與防撞梁和蒙皮的接觸,可顯著提高計算結(jié)果的精準度以及有限元模型的可靠性。
b.在進行模態(tài)計算結(jié)果與試驗結(jié)果相關(guān)性分析時,除要檢查主要零件清單及主要連接部位外,還應(yīng)考慮產(chǎn)品實際狀態(tài)與設(shè)計狀態(tài)之間的差異,并據(jù)此修正有限元模型,保證分析結(jié)果與試驗結(jié)果的一致性。
1 黃金陵主編.汽車車身設(shè)計.北京:機械工業(yè)出版社,2007.
2 G.R.Liu,S.S.Quek,著;龍述堯,等譯.有限元法實用教程.長沙:湖南大學出版社,2004.
3 曹渡主編.汽車內(nèi)外飾設(shè)計與實戰(zhàn).北京:機械工業(yè)出版社,2011.