基于matlab的汽車動力學(xué)仿真計算

鄭戰(zhàn)光，汪兆亮，王佳祥，朱 帥，黃世鵬

（廣西制造系統(tǒng)與先進制造技術(shù)重點實驗室廣西大學(xué)機械工程學(xué)院，廣西南寧530004）

汽車是日常生活中應(yīng)用最廣泛的交通工具，認(rèn)識和了解汽車的各種性能對于汽車性能的綜合評價有著十分重要的意義。汽車的一些常用性能包括：汽車的動力性、汽車的燃油經(jīng)濟性、汽車的制動性、汽車的操縱穩(wěn)定性、汽車的平順性、汽車的通過性。汽車性能試驗周期長、耗時耗力，需要特定的場地給汽車性能測試帶來了諸多不便，同時，為避免人為操作不當(dāng)?shù)纫虻挠绊?，目前常用matlab編程或者結(jié)合simulink仿真等建立相關(guān)數(shù)學(xué)模型的方法共同實現(xiàn)。本文著眼于汽車各種性能中最基本也是最重要的性能——汽車的動力性研究做了部分工作，為相關(guān)工作人員提供了一些參考。

1 汽車的動力性相關(guān)分析

汽車動力性評價指標(biāo)有最高車速、加速時間、最大爬坡度等，與之對應(yīng)的試驗內(nèi)容有最高車速的測試、汽車起步連續(xù)換擋加速時間與超車加速時間的測試和汽車最大爬坡度的測試[1]。另外，按照我國標(biāo)準(zhǔn)，動力性評價試驗均在滿載情況下進行。轎車的最高車速大，加速時間短，經(jīng)常在較好的道路上行駛，一般不強調(diào)它的爬坡能力；然而，它的Ⅰ擋加速能力大，故爬坡能力也強[2]。貨車在各種地區(qū)的各種道路上行駛，所以必須具有足夠的爬坡能力，一般最大爬坡度在30%即16.7°左右。需要進一步說明的是：最大爬坡度代表了汽車的極限爬坡能力，它應(yīng)該比實際行駛中遇到的道路最大坡度超出很多，這是因為應(yīng)考慮到實際坡道行駛時，在坡道上停車順利起步加速、克服松軟坡道路面的阻力、克服坡道上崎嶇不平路面的局部最大阻力等要求的緣故。本文利用matlab軟件編程實現(xiàn)汽車動力學(xué)模型的建立與求解，為用戶依據(jù)汽車相關(guān)參數(shù)來計算汽車動力性指標(biāo)提供了方法參考，同時也為開發(fā)車輛動力學(xué)軟件仿真奠定了一定的基礎(chǔ)。

2 汽車的最高車速、加速時間、最大爬坡度求解

動力性數(shù)學(xué)模型

汽車行駛方程：Ft=Ff+Fw+Fi+Fj

行駛方程具體化：

式中：Ft為驅(qū)動力（N）；Ff為滾動阻力（N）；Fw為空氣阻力（N）；Fi為坡度阻力（N）；Fj為加速阻力（N）；Ttq為發(fā)動機輸出轉(zhuǎn)矩（N·m）；ig、i0為變速器傳動比、主減速器傳動比；ηt為傳動系的機械效率；r為車輪滾動半徑（m）；G 為汽車重量（kg）；f為滾動阻力系數(shù)；i為道路坡度；CD為空氣阻力系數(shù)；A為迎風(fēng)面積（m2）；ua為車速（km/h）；δ 為汽車旋轉(zhuǎn)質(zhì)量換算系數(shù)；m為汽車質(zhì)量（kg）為加速度（m/s2）[4].

所研究汽車動力模型的參數(shù)如表1所列，模型變速器傳動比如表2所列。

表1 汽車動力模型參數(shù)[3]

表2 變速器傳動比（ig）[3]

程序設(shè)計流程：依據(jù)汽油發(fā)動機的使用外特性Tq-n曲線的擬合公式[3]計算出最低轉(zhuǎn)速與最高轉(zhuǎn)速范圍內(nèi)不同轉(zhuǎn)速下的轉(zhuǎn)矩大小，根據(jù)計算出來的不同轉(zhuǎn)速對應(yīng)下的轉(zhuǎn)矩、表1、表2所提供的汽車相關(guān)參數(shù)結(jié)合驅(qū)動力公式計算出變速器不同檔位下、不同車速所對應(yīng)的驅(qū)動力，結(jié)合計算出來的滾動阻力與空氣阻力大小，便可得到汽車驅(qū)動力-行駛阻力平衡圖，進而得到最高車速（在驅(qū)動力-行駛阻力平衡圖上Ft曲線與Ff+Fw曲線交點對應(yīng)的車速，即最高車速U－max）；依據(jù)計算出的不同檔位、不同車速下的驅(qū)動力、空氣阻力、汽車重量可求得汽車的動力因數(shù)，進而得到汽車在各檔位下的動力因數(shù)與車速的關(guān)系曲線（動力特性圖），汽車1檔下的最大動力因數(shù)求出最大爬坡度；依據(jù)求得的不同檔位及車速下的驅(qū)動力、滾動阻力、空氣阻力、質(zhì)量換算系數(shù)，可求出汽車的加速度大小，便得到加速度關(guān)于車速曲線、加速度導(dǎo)數(shù)關(guān)于車速曲線，進而可以求得汽車加速時間。

根據(jù)以上發(fā)動機外特性的參數(shù)可以計算出多組不同轉(zhuǎn)速下的轉(zhuǎn)矩，如表3所示。

表3 不同轉(zhuǎn)速下的轉(zhuǎn)矩

再根據(jù)發(fā)動機的外特性Tq-n曲線方程用matlab繪制出該汽車發(fā)動機輸出轉(zhuǎn)矩特性曲線，如圖1所示。

圖1 汽車發(fā)動機輸出轉(zhuǎn)矩特性曲線

汽車動力裝置參數(shù)一旦確定了，汽車的動力性基本上可以得到初步的認(rèn)識。該輕型貨車的相關(guān)參數(shù)出自文獻[3]，如表4所示。

表4 整車配置參數(shù)[3]

汽車在水平路面上勻速行駛時阻力包括滾動阻力和空氣阻力，運用式（1）可求出不同擋位下汽車的行駛速度。

繼而利用式（2）和式（3）分別求出滾動阻力和空氣阻力。

Ff+Fw與車速的關(guān)系為行駛阻力曲線，再結(jié)合式（4）所描繪的各檔驅(qū)動力曲線就可得到汽車驅(qū)動力—行駛阻力平衡關(guān)系，再用matlab繪制出汽車驅(qū)動力與阻力平衡圖，如圖2所示。

圖2 驅(qū)動力與阻力平衡圖

在汽車驅(qū)動力與阻力平衡圖中將第Ft5曲線與Ft+Fw曲線的交點便求出最高車速uamax，此時汽車最高車速=99.7353.

同樣，在汽車驅(qū)動力與阻力平衡圖中采用matlab進行仿真分析可以得到汽車最大爬坡度：1檔的最大爬坡度為I1max=0.428，1擋的最大爬坡角為alphamax=23.15°.

現(xiàn)利用圖3來求出汽車從0km/h全力加速至70km/h所用的加速時間。圖3為基于MATLAB軟件繪制的汽車動力特性圖。

圖3 汽車動力特性圖

由行駛方程式（5）并結(jié)合圖3所示的各檔節(jié)氣門全開時的加速度曲線，再由運動學(xué)知識可得式（6），再對式（6）積分就可得式（7）。

此時，加速時間可用計算機進行數(shù)值積分法或圖解積分法求出，曲線下兩個速度區(qū)間的面積就是通過此速度區(qū)的加速時間，采用的方法也就是數(shù)值積分法。經(jīng)過數(shù)值求解可得汽車從0 km/h全力加速至70 km/h所用的加速時間經(jīng)計算約為25.8223 s.

3 結(jié)束語

利用動力學(xué)公式在matlab中建立汽車動力性試驗的仿真模型，可得到汽車動力性的各項參數(shù)。設(shè)計合理的數(shù)學(xué)模型及高效的仿真程序，能得出接近真實試驗的結(jié)果，為工作人員提供了重要參考，有很強的實用性[5]。

[1]何耀華.汽車試驗技術(shù)[M].北京：機械工業(yè)出版社，2010.

[2]楊萬福.汽車?yán)碚揫M].廣州：華南理工大學(xué)出版社，2010.

[3]余志生.汽車?yán)碚揫M].北京：機械工業(yè)出版社，2009.

[4]肖啟瑞，樊明明.車輛工程仿真與分析—基于MATLAB的實現(xiàn)[M].北京：機械工業(yè)出版社，2010.

[5]王 薇.MATLAB從基礎(chǔ)到精通[M].北京：電子工業(yè)出版社，2012.