陳達(dá)亮李洪亮,顧燦松王海洋任　斌

（1.中國(guó)汽車(chē)技術(shù)研究中心，天津300300；2.河北工業(yè)大學(xué)，天津300130）

?

副車(chē)架系統(tǒng)邊界約束對(duì)計(jì)算模態(tài)分析精度的影響研究

陳達(dá)亮1李洪亮1,2顧燦松1王海洋1任斌2

（1.中國(guó)汽車(chē)技術(shù)研究中心，天津300300；2.河北工業(yè)大學(xué)，天津300130）

【摘要】針對(duì)某車(chē)型“剛性”副車(chē)架系統(tǒng)建立有限元模型，以約束處理方式為設(shè)計(jì)變量，以0～200 Hz范圍內(nèi)的約束模態(tài)頻率計(jì)算值與試驗(yàn)值的誤差作為評(píng)價(jià)指標(biāo)進(jìn)行正交試驗(yàn)設(shè)計(jì)，分析得到不同約束處理方式對(duì)副車(chē)架約束模態(tài)分析精度的影響規(guī)律。結(jié)果表明：副車(chē)架與車(chē)身連接點(diǎn)處車(chē)身局部剛度對(duì)模態(tài)分布的影響較大，懸架擺臂次之，襯套動(dòng)剛度對(duì)模態(tài)分布的影響相對(duì)較弱；基于正交試驗(yàn)遴選出的最佳邊界約束處理方式，可得約束模態(tài)頻率計(jì)算值與試驗(yàn)結(jié)果的相對(duì)誤差為0.45％，達(dá)到了理想精度。

1　副車(chē)架系統(tǒng)及其邊界約束

副車(chē)架連接車(chē)身和懸架，在一定程度上起到隔振降噪的作用。圖1為某型汽車(chē)的“剛性”副車(chē)架系統(tǒng)，其特點(diǎn)是副車(chē)架結(jié)構(gòu)通過(guò)螺栓“剛性”地固定于車(chē)身。副車(chē)架系統(tǒng)主要受2種激勵(lì)，即來(lái)自不平路面的激勵(lì)和來(lái)自發(fā)動(dòng)機(jī)的激勵(lì)。

一般情況下2種激勵(lì)都會(huì)通過(guò)副車(chē)架系統(tǒng)傳至車(chē)身，嚴(yán)重時(shí)會(huì)削弱車(chē)輛的乘坐舒適性。因此，副車(chē)架系統(tǒng)動(dòng)態(tài)特性直接影響汽車(chē)NVH性能。工程上常用的研究方法是通過(guò)模態(tài)分析來(lái)探究副車(chē)架系統(tǒng)的動(dòng)態(tài)特性。

圖1　某型汽車(chē)“剛性”副車(chē)架系統(tǒng)

計(jì)算機(jī)輔助工程（Computer Aided Engineering，CAE）具有高效率、低成本等優(yōu)點(diǎn)，因而成為研究結(jié)構(gòu)動(dòng)態(tài)特性的常用方法，在車(chē)架及副車(chē)架方面已有相關(guān)研究成果［1～5］。有限元建模過(guò)程中對(duì)約束邊界的不同處理方式會(huì)得到不同精度的仿真結(jié)果，文獻(xiàn)［6］通過(guò)有限元法計(jì)算某車(chē)型前副車(chē)架在自由狀態(tài)、忽略襯套連接的剛性約束狀態(tài)以及安裝狀態(tài)下（包含襯套連接）的模態(tài)并將計(jì)算值與試驗(yàn)值進(jìn)行對(duì)比，結(jié)果表明前副車(chē)架與車(chē)身搭接件安裝處的剛度是決定前副車(chē)架振動(dòng)頻率高低的主要原因。但單一副車(chē)架模型沒(méi)有考慮到下擺臂對(duì)副車(chē)架系統(tǒng)的影響，這會(huì)給計(jì)算結(jié)果帶來(lái)誤差，并且計(jì)算副車(chē)架系統(tǒng)安裝狀態(tài)下的約束模態(tài)時(shí)需要依托整車(chē)模型建模，這不僅增大了工作量而且降低了計(jì)算效率。文獻(xiàn)［7］利用AD?AMS軟件分析了副車(chē)架-橡膠襯套-車(chē)身12自由度剛?cè)狁詈险駝?dòng)系統(tǒng)的特性，結(jié)果表明副車(chē)架與車(chē)身連接處的襯套剛度應(yīng)相互匹配才能達(dá)到較好的效果，這說(shuō)明當(dāng)副車(chē)架的約束邊界涉及彈性部件時(shí)，彈性部件對(duì)激勵(lì)響應(yīng)的滯后程度也會(huì)影響系統(tǒng)的振動(dòng)傳遞特性。因此，對(duì)副車(chē)架約束邊界條件進(jìn)行細(xì)致分析具有重要的意義。

本文針對(duì)圖1所示的副車(chē)架系統(tǒng)采用有限元法建立模型，探究不同的邊界約束處理方式對(duì)副車(chē)架系統(tǒng)有限元模態(tài)分析精度的影響。即以邊界約束的處理方式為設(shè)計(jì)變量，以0～200 Hz（低頻）內(nèi)副車(chē)架系統(tǒng)計(jì)算約束模態(tài)頻率與試驗(yàn)值的誤差作為評(píng)價(jià)指標(biāo)進(jìn)行試驗(yàn)設(shè)計(jì)，分析不同的邊界約束處理方式與模態(tài)分析精度的關(guān)系，以期遴選出最佳的邊界約束處理方式。

2　邊界約束處理與有限元建模

2.1副車(chē)架系統(tǒng)的約束邊界

如圖1所示，副車(chē)架結(jié)構(gòu)通過(guò)螺栓固定于車(chē)身，通過(guò)橡膠襯套與橫向穩(wěn)定桿、轉(zhuǎn)向機(jī)殼體以及懸架擺臂連接，擺臂另一端通過(guò)球鉸與轉(zhuǎn)向節(jié)連接。副車(chē)架結(jié)構(gòu)及其約束邊界構(gòu)成了模態(tài)分析的對(duì)象系統(tǒng)。

懸架下擺臂一端與副車(chē)架通過(guò)橡膠襯套連接，其剛度較弱，使得被連接部件處于相對(duì)浮動(dòng)狀態(tài)；下擺臂另一端則通過(guò)球鉸與轉(zhuǎn)向節(jié)連接，具有完全的轉(zhuǎn)動(dòng)自由度。模型中，對(duì)懸架下擺臂的處理方式有如下2種。

a.由于懸架擺臂與副車(chē)架之間約束較弱，因此將擺臂舍去，使副車(chē)架在與擺臂連接點(diǎn)位置處于自由狀態(tài)；

b.將懸架擺臂視為副車(chē)架系統(tǒng)的一部分，擺臂與轉(zhuǎn)向節(jié)的球鉸連接因而視為副車(chē)架系統(tǒng)的約束邊界。

副車(chē)架通過(guò)橡膠襯套與橫向穩(wěn)定桿、轉(zhuǎn)向機(jī)殼體連接，而橫向穩(wěn)定桿和轉(zhuǎn)向機(jī)殼體的剛度遠(yuǎn)大于襯套剛度，并且它們與車(chē)身之間為剛性連接，因此，模型中將橫向穩(wěn)定桿和轉(zhuǎn)向機(jī)殼體視為車(chē)身的一部分。

所謂“剛性”副車(chē)架，是其與車(chē)身連接處不設(shè)置彈性減振與緩沖元件，因而具有相對(duì)較大的局部剛度。針對(duì)這一特點(diǎn)，模型中對(duì)副車(chē)架與車(chē)身連接處約束的處理有2種方式：

a.將副車(chē)架與車(chē)身視為剛性連接；

b.用襯套模擬副車(chē)架與車(chē)身連接，襯套剛度賦值為副車(chē)架與車(chē)身連接處的車(chē)身局部剛度。

橡膠襯套剛度特性主要分為靜剛度和動(dòng)剛度，對(duì)襯套剛度的處理，此處提供2種方式：

a.只考慮襯套靜剛度，這種處理方式的優(yōu)點(diǎn)是數(shù)據(jù)測(cè)量簡(jiǎn)單且在軟件中賦值方便，不足之處是沒(méi)有考慮襯套的動(dòng)態(tài)特性；

b.考慮襯套動(dòng)剛度，這種處理方式的優(yōu)點(diǎn)是可使模型更符合實(shí)際情況，缺點(diǎn)是數(shù)據(jù)測(cè)量比較困難且部分有限元分析軟件不支持橡膠襯套的動(dòng)剛度曲線輸入。

綜上所述，處理副車(chē)架系統(tǒng)約束邊界主要從3個(gè)角度出發(fā)：是否將麥弗遜懸架下擺臂視為副車(chē)架系統(tǒng)的一部分；是否引入車(chē)身局部剛度以模擬副車(chē)架系統(tǒng)與車(chē)身的連接；是否考慮橡膠襯套動(dòng)剛度。

2.2基礎(chǔ)參數(shù)獲取與有限元建模

根據(jù)上述對(duì)約束邊界的討論，確定需要通過(guò)測(cè)量手段獲取的建?；A(chǔ)參數(shù)主要包括橡膠襯套靜剛度、動(dòng)剛度和副車(chē)架與車(chē)身連接點(diǎn)處的局部剛度。分別采用MTS831.50彈性體試驗(yàn)臺(tái)和MTS832襯套試驗(yàn)臺(tái)測(cè)量橡膠襯套靜剛度和動(dòng)剛度，并以動(dòng)剛度曲線在0～200 Hz范圍內(nèi)的RMS值來(lái)近似代替動(dòng)剛度。相應(yīng)的測(cè)量和計(jì)算結(jié)果如表1和表2所示，其中各坐標(biāo)皆處于整車(chē)坐標(biāo)系下。另一方面，通過(guò)試驗(yàn)測(cè)試副車(chē)架與車(chē)身連接處的局部剛度，結(jié)果如表3所示。限于測(cè)試條件，連接部位的角剛度無(wú)法測(cè)量，因此在模型中帶入經(jīng)驗(yàn)值。

表1　橡膠襯套靜剛度

表2　橡膠襯套動(dòng)剛度RMS值

表3　車(chē)身安裝點(diǎn)處車(chē)身鋼板局部動(dòng)剛度　N/mm

在劃分有限元網(wǎng)格前，對(duì)幾何模型進(jìn)行特征清理并刪除對(duì)模態(tài)影響較小的部分圓角和小孔特征。對(duì)幾何模型抽取中面并采用四邊形殼單元?jiǎng)澐志W(wǎng)格，經(jīng)多次調(diào)試對(duì)比最終確定網(wǎng)格大小為8 mm。該副車(chē)架材料為45#鋼。所建立有限元模型的單元數(shù)為12 891，節(jié)點(diǎn)數(shù)為16 575，并可靈活設(shè)置不同的邊界約束，如圖2所示。

圖2　副車(chē)架有限元模型

3　邊界約束處理方式的試驗(yàn)設(shè)計(jì)

3.1基本思想

正交試驗(yàn)設(shè)計(jì)方法是20世紀(jì)50年代至70年代之間由日本統(tǒng)計(jì)學(xué)家田口玄一提出的一種試驗(yàn)設(shè)計(jì)方法。應(yīng)用這一方法研究變量對(duì)試驗(yàn)結(jié)果的影響時(shí)，可以在不影響分析精度的前提下顯著減少試驗(yàn)次數(shù)［8～10］。這里，以副車(chē)架系統(tǒng)的邊界約束處理方式作為試驗(yàn)的設(shè)計(jì)變量，以0～200 Hz低頻范圍內(nèi)的有限元計(jì)算模態(tài)頻率誤差作為評(píng)價(jià)指標(biāo)，進(jìn)行邊界約束處理方式的試驗(yàn)設(shè)計(jì)。此處，計(jì)算模態(tài)頻率的誤差定義為其與相應(yīng)的試驗(yàn)?zāi)B(tài)分析結(jié)果之間的差異。因此，首先要對(duì)副車(chē)架系統(tǒng)進(jìn)行相應(yīng)的模態(tài)測(cè)試。

3.2副車(chē)架約束模態(tài)試驗(yàn)

在實(shí)車(chē)安裝狀態(tài)下，通過(guò)測(cè)試獲取副車(chē)架系統(tǒng)試驗(yàn)?zāi)B(tài)分析結(jié)果。傳感器布置如圖3所示，測(cè)試采用單輸入多輸出的方式進(jìn)行。測(cè)試結(jié)果表明，在0～200 Hz頻率范圍內(nèi)，副車(chē)架系統(tǒng)只有1階模態(tài)，模態(tài)頻率點(diǎn)為175.9 Hz，模態(tài)振型如圖4所示。

圖3　副車(chē)架系統(tǒng)約束模態(tài)試驗(yàn)示意

圖4　副車(chē)架系統(tǒng)試驗(yàn)約束模態(tài)振型

3.3試驗(yàn)設(shè)計(jì)結(jié)果分析

針對(duì)2.1節(jié)提出的約束邊界的3種處理方式，在試驗(yàn)設(shè)計(jì)中將其視為影響建模精度的3種因素并對(duì)每種因素設(shè)置2個(gè)水平：對(duì)于襯套剛度，數(shù)字1表示考慮靜剛度，數(shù)字2表示考慮動(dòng)剛度，即模型中帶入橡膠襯套動(dòng)剛度在0～200 Hz內(nèi)的RMS值；對(duì)于擺臂，數(shù)字1表示在模型中不考慮副車(chē)架與懸架擺臂的連接關(guān)系，即副車(chē)架與擺臂連接點(diǎn)位置處于自由狀態(tài)，數(shù)字2表示將擺臂視為副車(chē)架系統(tǒng)的一部分，并將擺臂與轉(zhuǎn)向節(jié)的球鉸連接視為副車(chē)架系統(tǒng)的約束邊界；對(duì)于車(chē)身安裝點(diǎn)剛度，數(shù)字1表示將車(chē)身視為純剛體，即副車(chē)架與車(chē)身連接處鋼板的剛度為無(wú)窮大，數(shù)字2表示考慮車(chē)身安裝點(diǎn)處車(chē)身鋼板的局部剛度，并在模型中用橡膠襯套模擬副車(chē)架與車(chē)身的連接關(guān)系。由此進(jìn)行試驗(yàn)設(shè)計(jì)。

在所有約束條件下均可將車(chē)身視為與大地固定的剛體，在有限元模型中分別添加第1、2、3、4組約束條件，計(jì)算約束模態(tài)并與試驗(yàn)值對(duì)比，以模態(tài)頻率誤差作為評(píng)價(jià)指標(biāo)，計(jì)算分析結(jié)果如表4和表5所示。

表4　正交試驗(yàn)計(jì)算結(jié)果

表5　正交試驗(yàn)分析結(jié)果

表5中，yk為某因素k（k=1/2）水平對(duì)應(yīng)的試驗(yàn)指標(biāo)（即模態(tài)頻率誤差）之和；為yk的平均值；R為某因素的極差，即該因素不同水平對(duì)應(yīng)平均值中最大值與最小值的差值。應(yīng)用極差分析法可以判斷每個(gè)因素對(duì)指標(biāo)影響程度的大小。從表5中的極差R可以看出，副車(chē)架與車(chē)身連接關(guān)系的處理對(duì)建模精度的影響相對(duì)最大，橡膠襯套是否考慮動(dòng)剛度特性對(duì)建模精度的影響相對(duì)最小，而擺臂的處理方式對(duì)建模精度的影響介于兩者之間。

3.4試驗(yàn)設(shè)計(jì)結(jié)果驗(yàn)證

根據(jù)正交試驗(yàn)結(jié)果選出約束邊界處理方法的最優(yōu)組合如表5中“優(yōu)水平”一欄所示：考慮襯套動(dòng)剛度特性、將擺臂引入副車(chē)架系統(tǒng)以及將副車(chē)架與車(chē)身連接視為柔性。依據(jù)這一組合修改有限元模型約束邊界后重新計(jì)算約束模態(tài)結(jié)果如圖5所示。

圖5　計(jì)算模態(tài)振型

依據(jù)正交試驗(yàn)結(jié)果對(duì)副車(chē)架約束邊界進(jìn)行修改后，有限元模型計(jì)算約束模態(tài)頻率為176.7 Hz，試驗(yàn)值為175.9 Hz，誤差為0.45％。相比于正交試驗(yàn)中其它約束邊界處理方式組合，經(jīng)正交試驗(yàn)優(yōu)化遴選出的約束邊界處理方式在模型中的計(jì)算約束模態(tài)頻率與試驗(yàn)值誤差最小。

另一方面，對(duì)比圖5和圖4中的副車(chē)架系統(tǒng)模態(tài)振型，可以發(fā)現(xiàn)兩者具有很高的相似性：它們均為副車(chē)架繞Y軸（整車(chē)坐標(biāo)系）的俯仰振動(dòng)，沒(méi)有沿Y軸的平動(dòng)。

綜上所述，經(jīng)正交試驗(yàn)選出的邊界約束處理方式應(yīng)用到有限元模型中后，計(jì)算約束模態(tài)與試驗(yàn)約束模態(tài)存在較好的吻合。

4　結(jié)束語(yǔ)

a.“剛性”副車(chē)架系統(tǒng)與車(chē)身連接點(diǎn)處局部剛度對(duì)有限元模態(tài)分析精度的影響相對(duì)較大，擺臂次之，橡膠襯套剛度對(duì)模型精度的影響相對(duì)前兩者要弱。因此，在模型中需要細(xì)致處理副車(chē)架與車(chē)身連接點(diǎn)處的約束關(guān)系，使連接點(diǎn)處的局部剛度特性能夠在模型中正確體現(xiàn)。

b.經(jīng)試驗(yàn)設(shè)計(jì)分析遴選出的邊界約束處理方式顯著提高了模型分析計(jì)算的精度，表明了該方法的有效性，并可為類(lèi)似問(wèn)題的研究提供參考。

參考文獻(xiàn)

1鄭灝.汽車(chē)前副車(chē)架有限元分析及優(yōu)化：［學(xué)位論文］.武漢：武漢理工大學(xué)，2012.

2劉寶生.基于有限元法的摩托車(chē)車(chē)架振動(dòng)響應(yīng)與可靠性分析：［學(xué)位論文］.天津：天津大學(xué)，2010.

3鄭松林，趙德彪，馮金芝，等.基于強(qiáng)度與模態(tài)靈敏度分析的轎車(chē)前副車(chē)架輕量化設(shè)計(jì).機(jī)械設(shè)計(jì)，2012，29（2）：92～96.

4王銳，蘇小平.汽車(chē)副車(chē)架強(qiáng)度模態(tài)分析及結(jié)構(gòu)優(yōu)化.機(jī)械設(shè)計(jì)與制造，2015（4）：152～154.

5王曉明，孟令超，沈勇，等.裝載機(jī)后車(chē)架動(dòng)態(tài)特性分析及結(jié)構(gòu)改進(jìn).建筑機(jī)械，2014（8）：63～70.

6孫風(fēng)蔚，徐昊，陳杰龍，等.某車(chē)前副車(chē)架模態(tài)分析與改進(jìn).汽車(chē)工程學(xué)報(bào)，2012，2（4）：303～306.

7趙帥.轎車(chē)副車(chē)架隔振性能的動(dòng)態(tài)建模研究：［學(xué)位論文］.長(zhǎng)春：吉林大學(xué)，2014.

8姜立標(biāo)，侯文超，谷方德.基于正交試驗(yàn)的汽車(chē)穩(wěn)態(tài)回轉(zhuǎn)性能仿真.北京航空航天大學(xué)學(xué)報(bào)，2011，37（8）：968～972.

9張晗，過(guò)學(xué)迅，楊波，等.基于正交試驗(yàn)的載貨汽車(chē)轉(zhuǎn)向回正性仿真試驗(yàn)研究.機(jī)械設(shè)計(jì)與制造，2013（5）：137～139.

10廖鴻胡，趙云飛，成艾國(guó)，等.基于正交試驗(yàn)的汽車(chē)滑移門(mén)平順性優(yōu)化研究.機(jī)電工程，2010，27（12）：49～52.

（責(zé)任編輯斛畔）

修改稿收到日期為2015年8月19日。

主題詞:副車(chē)架系統(tǒng)有限元法邊界約束模態(tài)分布正交試驗(yàn)設(shè)計(jì)

The Effect of Subframe System Boundary Constraints on the Accuracy of Computational Modal Analysis

Chen Daliang1, Li Hongliang1,2, Gu Cansong1, Wang Haiyang1, Renbin2
（1.China Automotive Technology and Research Center, Tianjin 300300; 2. Hebei University of Technology, Tianjin 300130）

【Abstract】FE model of a vehicle "rigid" subframe system is established, the influence law of different constraint processing modes on the accuracy of subframe constraint modal analysis is obtained through analysis with the constraint processing mode as design variable, and the error between constraint computational modal frequency and experimental value in the range of 0～200 Hz as evaluation criteria for the orthogonal experimental design. The results show that the local stiffness of the joining part between subframe and body is the most important factor to affect modal distribution, followed by the suspension swing arm, whereas the sleeve dynamic stiffness affects modal distribution slightly; the optimal boundary constraint processing mode is selected based on orthogonal experiment, and it is concluded that the relative error of the computational modal frequency compared with the experimental one is 0.45％, thus a relative ideal accuracy is obtained.

Key words:Subframe system, Finite element method, Boundary constraint, Modal distrubution, Orthogonal experimental design

中圖分類(lèi)號(hào):U463.83

文獻(xiàn)標(biāo)識(shí)碼:A

文章編號(hào):1000-3703（2016）04-0027-04