宋大朝，繆聰雄，周文峰，姜武華

（合肥工業(yè)大學(xué) 機(jī)械與汽車工程學(xué)院，安徽 合肥 230009）

0 引 言

中國大學(xué)生純電動(dòng)方程式汽車大賽（Formula Student Electric China，簡稱FSEC）是一項(xiàng)由高等院校的汽車工程或汽車相關(guān)專業(yè)的在校學(xué)生組隊(duì)參加的汽車設(shè)計(jì)與制造比賽。這項(xiàng)比賽起源于中國大學(xué)生方程式汽車大賽，至今已舉辦過 2屆。大學(xué)生方程式賽車在國際上被視為“學(xué)界的F1方程式賽車”。該賽事是一項(xiàng)整車設(shè)計(jì)制造方面的比賽，也是一場(chǎng)“工程教育式的體驗(yàn)”。

懸架作為汽車上連接車輪與車身的彈性連接系統(tǒng)，主要由彈簧、減震器、導(dǎo)向機(jī)構(gòu)和緩沖塊組成[1]。懸架設(shè)計(jì)的好壞，直接影響整車的操縱穩(wěn)定性和行駛平順性，對(duì)于FSEC賽車而言，懸架設(shè)計(jì)更需要考慮整車的操縱穩(wěn)定性。汽車懸架有著非常復(fù)雜的幾何運(yùn)動(dòng)學(xué)和彈性運(yùn)動(dòng)學(xué)特性，由于FSEC賽車都通過鉸鏈來連接桿件，沒有襯套作為連接體，它的懸架幾何運(yùn)動(dòng)學(xué)特性就顯得更加重要[2]。文中對(duì)某高校一款FSEC賽車雙橫臂后懸架進(jìn)行了詳細(xì)的幾何運(yùn)動(dòng)學(xué)分析，并對(duì)懸架的前束、外傾角、輪心縱向和側(cè)向位移變化率進(jìn)行了優(yōu)化。在此基礎(chǔ)上，對(duì)懸架搖臂設(shè)計(jì)進(jìn)行了分析，以懸架搖臂各桿件的長度為參數(shù)，得到相應(yīng)的杠桿比變化曲線，確定其對(duì)賽車的杠桿比、偏頻的影響，驗(yàn)證了采用可調(diào)式搖臂改變杠桿比，對(duì)懸架剛度進(jìn)行調(diào)節(jié)的可行性。

1 模型建立

以某高校車隊(duì)中推桿形式的雙橫臂懸架數(shù)據(jù)為基礎(chǔ)，建立懸架、推桿、搖臂以及避震器模型。硬點(diǎn)坐標(biāo)參數(shù)見表1。

表1 雙橫臂后懸架硬點(diǎn)坐標(biāo)mm

續(xù)表1

根據(jù)硬點(diǎn)信息和結(jié)構(gòu)尺寸參數(shù)，建立懸架模型[3]。其中懸架桿件都通過約束副連接，整個(gè)懸架都以剛體形式建模，忽略桿件柔性對(duì)懸架運(yùn)動(dòng)學(xué)特性的影響，也不考慮實(shí)車鉸鏈的間隙、剛度、阻尼等的影響，所建立的懸架系統(tǒng)如圖1所示。

2 懸架幾何運(yùn)動(dòng)學(xué)仿真

懸架運(yùn)動(dòng)學(xué)分析及優(yōu)化包括幾何運(yùn)動(dòng)學(xué)和彈性運(yùn)動(dòng)學(xué) 2個(gè)方面，主要通過懸架的 K&C（Kinematica&Compliance）特性進(jìn)行定量分析[4]。懸架K&C特性決定車輪跳動(dòng)時(shí)的姿態(tài)、位置和車輛行駛性能。K&C特性通常分為2個(gè)方面：

1）懸架幾何運(yùn)動(dòng)學(xué)特性：車輪發(fā)生垂直位移時(shí)懸架幾何學(xué)位置引起的其他 5個(gè)自由度定位變化，簡稱K特性，主要與懸架桿件運(yùn)動(dòng)的幾何參數(shù)有關(guān)；

2）懸架彈性運(yùn)動(dòng)學(xué)特性：車輪受力時(shí)襯套撓曲變形以及導(dǎo)向桿件彎曲引起的車輪 6自由度定位變化，也稱順從特性，簡稱C特性，主要與懸架彈性元件的剛度特性有關(guān)。

因該懸架不存在橡膠襯套之類的彈性元件，且不考慮懸架桿件在受力時(shí)發(fā)生的彈性變形，所以在 Adams/Car中只對(duì)懸架的幾何運(yùn)動(dòng)學(xué)特性進(jìn)行分析和優(yōu)化[5]。在四輪定位參數(shù)中，前束和外傾角隨輪跳的變化率對(duì)整車的操縱穩(wěn)定性影響很大，應(yīng)盡量減小車輪上下跳動(dòng)時(shí)這 2個(gè)參數(shù)的變化，保證整車的操縱穩(wěn)定性；其次，輪跳時(shí)軸距和輪距的變化也會(huì)影響整車的操縱穩(wěn)定性[6]。

在 Adams/Car中，設(shè)置雙輪同向跳動(dòng)，范圍為-30～+30 mm，前束初始設(shè)置為0，外傾角設(shè)置為-1°。查看前束、外傾角、輪心縱向和側(cè)向位移隨輪跳的變化率。

2.1 前束的變化

前束角變化過大會(huì)影響賽車直線行駛的穩(wěn)定性并使輪胎磨損加劇，所以在FSEC賽車設(shè)計(jì)時(shí)希望前束角的變化量越小越好。前束隨輪胎垂向跳動(dòng)的變化率如圖2所示。

2.2 外傾角的變化

車輪跳動(dòng)時(shí)外傾角的變化對(duì)賽車的穩(wěn)態(tài)響應(yīng)特性有很大影響。為了最大限度地發(fā)揮輪胎性能，使賽車在轉(zhuǎn)彎中能夠獲得最大側(cè)向力，設(shè)計(jì)時(shí)常把外傾角設(shè)為負(fù)值，并且希望其隨車輪跳動(dòng)的變化盡可能小[7]。外傾角隨輪胎垂向跳動(dòng)變化率如圖3所示。

2.3 輪心縱向位移的變化

車輪跳動(dòng)時(shí)，輪心縱向位移即賽車軸距也會(huì)隨之改變，變化情況如圖4所示。

2.4 輪心側(cè)向位移的變化

車輪跳動(dòng)時(shí)，輪心側(cè)向位移即賽車輪距也會(huì)隨之改變，變化情況如圖5所示。

3 幾何運(yùn)動(dòng)學(xué)優(yōu)化

對(duì)于設(shè)計(jì)初期的懸架，在運(yùn)動(dòng)學(xué)仿真時(shí)，會(huì)存在某些參數(shù)變化過大的情況，對(duì)此，需運(yùn)用Adams/Insight或手動(dòng)調(diào)試來優(yōu)化參數(shù)[8]。在上述仿真中，分別查看了前束、外傾角、輪心縱向和側(cè)向位移隨輪胎垂向跳動(dòng)的變化率，據(jù)此將優(yōu)化變量設(shè)置為懸架上下橫臂內(nèi)側(cè)與車架連接點(diǎn)的硬點(diǎn)坐標(biāo)，即：hpl_UCA_inner_front、hpl_UCA_inner_rear、hpl_LCA_inner_front、hpl_LCA_inner_rear這4個(gè)硬點(diǎn)的X、Y、Z三向坐標(biāo)，共12個(gè)變量，設(shè)置在±5mm范圍內(nèi)變化，以各輸出曲線的絕對(duì)最大值為目標(biāo)函數(shù)，分析硬點(diǎn)對(duì)這些參數(shù)的靈敏度。通過Adams/Insight分析后，各硬點(diǎn)參數(shù)對(duì)輸出曲線的絕對(duì)最大值影響最大的前 5個(gè)參數(shù)如圖6所示（以前束靈敏度分析為例）。

對(duì)12個(gè)變量進(jìn)行靈敏度分析后，發(fā)現(xiàn)這些變量對(duì)前束和輪心縱向位移的影響較大，對(duì)外傾角和輪心側(cè)向位移的影響較小，其中，hpl_LCA_inner_rear_z增加可對(duì)上述目標(biāo)函數(shù)進(jìn)行改善，特別是對(duì)前束的改善非常明顯；hpl_UCA_inner_rear_z增加可對(duì)外傾角、輪心側(cè)向和縱向位移進(jìn)行優(yōu)化，對(duì)前束變化起了反作用，但影響不大；hpl_UCA_inner_front_z增加僅對(duì)輪心縱向位移影響較大；hpl_LCA_ inner_front_z減小可改善前束和輪心縱向位移，但對(duì)外傾角和輪心側(cè)向位移的影響較小。綜合考慮上述條件和整車布置，對(duì)3個(gè)硬點(diǎn)坐標(biāo)的Z向進(jìn)行調(diào)整，見表2。

表2 優(yōu)化坐標(biāo)值mm

其中 hpl_LCA_inner_rear_z的增加改善了全部目標(biāo)函數(shù)，hpl_UCA_inner_rear_z的增加改善了外傾角、輪心側(cè)向和縱向位移，但造成了輪心縱向位移調(diào)整過度，所以對(duì) hpl_UCA_inner_front_z進(jìn)行反向調(diào)整，抵消前 2個(gè)變量調(diào)整造成的輪心縱向位移調(diào)整過度的情況。優(yōu)化前后對(duì)比如圖7～圖10所示。

4 懸架搖臂分析

首先以不可調(diào)搖臂為例，在軟件中設(shè)置輪胎以滿載時(shí)輪胎接地為基準(zhǔn)，上下跳動(dòng)值分別為 30 mm，測(cè)量彈簧壓縮值，并將輪胎上跳值與彈簧壓縮值的比值作為不可調(diào)搖臂的杠桿比曲線輸出。搖臂模型如圖11所示。

在圖11中，搖臂三角形的A點(diǎn)為搖臂與推桿連接點(diǎn)，B點(diǎn)為搖臂與車架連接點(diǎn)，C點(diǎn)為搖臂與避震器連接點(diǎn)。AB長40 mm，BC長65 mm，AC長 70 mm，3點(diǎn)在 Adams/Car中的坐標(biāo)值為：A（1500.0，289.286，321.965），B（1500.0，269.843，287.008），C（1500.0，219.512，327.887）。

通過仿真，可以得出初始搖臂的杠桿比曲線如圖12所示。

接著對(duì)搖臂參數(shù)進(jìn)行修改，將AB加長5 mm至45 mm，在仿真模型中，點(diǎn)B、點(diǎn)C坐標(biāo)不動(dòng)，對(duì)點(diǎn)A坐標(biāo)進(jìn)行修改，修改后的搖臂杠桿比曲線如圖13所示。

從圖13可以看出，加長后的杠桿比明顯增大，繼續(xù)對(duì)AB加長5 mm至50 mm，在仿真模型中，點(diǎn)B、點(diǎn)C坐標(biāo)不動(dòng)，對(duì)點(diǎn)A坐標(biāo)進(jìn)行修改，修改后的搖臂杠桿比曲線如圖14所示。

在相同的懸架模型基礎(chǔ)上，通過對(duì) 3組搖臂三角形的結(jié)構(gòu)設(shè)計(jì)進(jìn)行比較，可以看到，只對(duì)搖臂三角形的參數(shù)進(jìn)行細(xì)微改變，就可使杠桿比在較大的范圍內(nèi)發(fā)生改變，如圖15所示，AB邊在初始長度時(shí)杠桿比范圍為1.78～2.05，加長5 mm后的杠桿比范圍為2.08～2.34，加長10 mm后的杠桿比范圍為2.43～2.68。這在實(shí)際運(yùn)用中是非常不錯(cuò)的變化范圍，具有很強(qiáng)的現(xiàn)實(shí)意義。

5 結(jié) 論

在基于 Adams/Car的 FSEC懸架運(yùn)動(dòng)學(xué)分析及優(yōu)化和懸架搖臂設(shè)計(jì)的仿真分析中，通過對(duì)懸架硬點(diǎn)進(jìn)行調(diào)整，改善輪胎垂向跳動(dòng)時(shí)前束、外傾角、輪心縱向與側(cè)向位移的變化率；并通過調(diào)整懸架搖臂的尺寸參數(shù)，對(duì)FSEC懸架的杠桿比變化范圍進(jìn)行分析，驗(yàn)證了通過調(diào)整搖臂尺寸來改變懸架杠桿比的可行性，對(duì)FSEC懸架搖臂設(shè)計(jì)起到了預(yù)判和指導(dǎo)作用，具有較大的實(shí)用價(jià)值。

[1]米奇克，瓦倫托維茲. 汽車動(dòng)力學(xué)[M]. 陳蔭三，余強(qiáng)，譯. 北京：清華大學(xué)出版社，2009.

[2]迪克遜. 減振器手冊(cè)[M]. 李惠彬，孫振蓮，金婷，譯. 北京：機(jī)械工業(yè)出版社，2011.

[3]范成建，熊光明，周明飛.虛擬樣機(jī)軟件MSC．Adams應(yīng)用與提高[M]. 北京：機(jī)械工業(yè)出版社，2006.

[4]Milliken W F．Race Car Vehicle Dynamics [M]．United States of America：Society of Automotive Engineers，1995：61l-620.

[5]劉虹，王其東.基于 Adams雙橫臂獨(dú)立懸架的運(yùn)動(dòng)仿真分析［J］．合肥工業(yè)大學(xué)學(xué)報(bào)，2007，30（1）：57-59.

[6]劉美燕.FSAE賽車懸架仿真分析及操縱穩(wěn)定性虛擬試驗(yàn)[D].長沙：湖南大學(xué)，2008：39-47.

[7]貝爾恩德·海森英，漢斯·于爾根·布蘭德耳. 汽車行駛動(dòng)力學(xué)性能的主觀評(píng)價(jià)[M]. 石曉明，陳禎福，譯. 北京：人民交通出版社，2010：45-86.

[8]倪俊，徐彬. FSAE賽車雙橫臂前懸架運(yùn)動(dòng)學(xué)仿真及優(yōu)化[J]. 車輛與動(dòng)力技術(shù)，2011（4）：51-54.