基于MATLAB的純電動汽車動力性仿真分析

李 凱，冶存良，劉 彬

（寧夏理工學院 機械工程學院，寧夏 石嘴山 753000）

隨著世界科技革命和產業變革進程加快，能源短缺、環境污染等問題受到越來越多的重視。為應對環境問題、促進綠色發展，發展新能源汽車是一項行之有效的戰略[1]。目前由于純電動汽車具有無污染、結構簡單和能量轉換率高的優勢，其已成為新能源汽車的重點發展方向。動力性作為純電動汽車重要性能之一，進行純電動汽車動力性仿真，對提高純電動汽車設計水平、促進新能源汽車的發展有重要意義。

本文以比亞迪某純電動汽車作為研究對象，選取最高車速、加速能力和爬坡能力作為動力性評價指標，根據新能源汽車理論和依靠MATLAB軟件強大的計算及數據分析能力，選取車型相關實際參數（如表1所示），采用MATLAB編寫車輛動力模擬仿真程序，對純電動汽車動力性進行分析，獲得動力性仿真曲線，進而分析影響電動汽車動力的主要因素。

表1 整車主要技術參數

1 純電動汽車動力模型建立

1.1 電動汽車行駛方程

純電動汽車的動力性主要取決于電動汽車驅動力和行駛阻力，電動汽車的驅動力和傳統燃油車來源不同，電動汽車的驅動力主要來自電機的輸出轉矩。電動汽車驅動力計算公式為

式中，Ft為電動汽車驅動力；Te為電機輸出轉矩，it為總傳動比；ηt為傳動效率；r為滾動半徑。

電機輸出轉矩和轉速之間有一定的關系，這也是電動汽車動力性計算的主要依據。電機輸出轉矩的計算公式為

式中，Pe為電機高速功率；Tc為電機低速功率；n為電機轉速；nb為電機基速。

電動汽車的行駛速度計算公式為

式中，u為電動汽車的行駛速度；r為滾動半徑。

純電動汽車在行駛過程中受到的阻力主要有滾動阻力、空氣阻力、坡道阻力和加速阻力[2]。電動汽車滾動阻力與輪胎行駛單位距離能量損失有關，是輪胎和路面的變形引起的。滾動阻力計算公式為

式中，Ff為滾動阻力；m為電動汽車整車質量；f為滾動阻力系數；αG為坡度角；g為重力加速度，取g=9.8 m/s2。

電動汽車的空氣阻力和行駛速度、迎風面積和空阻系數有關，計算公式如下：

式中，Fw為電動汽車空氣阻力；CD為電動汽車空阻系數；A為電動汽車迎風面積；u為電動汽車的行駛速度。

電動汽車行駛過程中遇到上坡，需要克服重力沿坡道方向的坡道阻力，計算公式為

式中，Fi為電動汽車受到的坡道阻力。

加速阻力是指電動汽車在正常行駛中加速行駛時要克服的汽車本身質量加速運動時的慣性力，計算公式為

式中，Fj為電動汽車加速阻力；δ為電動汽車的旋轉質量換算系數。

由式（1）—式（7）可知，電動汽車克服行駛阻力行駛的動力學方程為

通過式（8）進而可以求得汽車動力性能的三大評價指標，分別是最高車速、加速能力和最大爬坡度[3]。

1.2 電動汽車最高車速

電動汽車最高車速和傳統燃油汽車的最高車速定義類似，在不考慮風速的條件下，電動汽車在平坦路面上行駛所能達到的最大車速為最高車速umax，此時電動汽車電機處于恒功率區間運行，驅動力和滾動阻力保持平衡[4]。電動汽車能達到的最高車速越高，則表明其動力越強。根據式（5）可得到電動汽車最高車速為

上述情況下驅動力曲線與滾動阻力-空氣阻力之和曲線有交點，交點所對應的車速就是電動汽車的最高車速。

當電動汽車驅動力大于滾動阻力和驅動力之和時，電動汽車的最高車速由電機的決定，這種情況下驅動力曲線與滾動阻力-空氣阻力之和曲線沒有交點[5]。電動汽車最高車速的計算公式為

式中，nmax為電機最大轉速。

1.3 電動汽車加速能力

電動汽車在平坦路面上行駛的加速度aj根據式（7）可求得

電動汽車的加速能力是指電動汽車的速度從u1加速到u2所耗費的時間，時間越短說明電動汽車的加速能力越強。計算公式為

式中，t為電動汽車加速時間。

1.4 電動汽車爬坡能力

電動汽車的爬坡能力用其滿負載情況下在良好不顛簸路面勻速行駛的最大爬坡度來表征，爬坡角度越大，說明電動汽車的爬坡能力越強[6]。

2 程序設計和仿真分析

2.1 純電動汽車動力性仿真程序設計

根據建立的式（9）—式（12）純電動汽車動力模型，使用MATLAB軟件編寫純電動汽車動力性仿真程序并運行，求得比亞迪某電動車型的最高速度、加速時間和最大爬坡度，并輸出仿真圖形。

2.2 純電動汽車仿真分析

在MATLAB R2022a編譯器中輸入程序并運行，計算得到該型純電動汽車的最高車速為189.10 km/h，加速時間為 7.90 s，最大爬坡度為38.46%。輸出的該型電動汽車驅動力-行駛阻力平衡圖、加速時間曲線和爬坡度曲線分別如圖1—圖3所示。

圖1 驅動力-行駛阻力平衡圖

圖2 加速時間曲線

圖3 最大爬坡度曲線

3 純電動汽車動力性影響因素分析

3.1 最高車速的影響因素分析

最高車速的影響因素主要有電機峰值功率、整車質量和空阻系數等。當車輛滿載且運行在平坦路面時，最高車速umax所對應的電機峰值功率Pm為

由式（9）可得空阻系數CD和最高車速umax之間的關系為

根據式（13）—式（14），編寫電機峰值功率和空阻系數對該型電動汽車最高車速影響的MATLAB程序，運行得到該型電動汽車最高車速與電機峰值功率和空阻系數這二者關系的仿真曲線，如圖4和圖5所示。

圖4 最高車速和峰值功率關系曲線

圖5 最高車速和空阻系數關系曲線

由圖4可以看出電機峰值功率越大，電動汽車的最高車速越高。從圖5可以看出空阻系數越大，電動汽車的最高車速越小。

3.2 加速能力的影響因素分析

電動汽車的加速能力由其百公里加速時間來表征，編寫電機峰值功率和空阻系數對電動汽車加速能力影響的MATLAB程序，運行得到該型電動汽車最高車速與電機峰值功率和空阻系數這二者關系的仿真曲線，如圖6和圖7所示。

圖6 加速時間和峰值功率關系曲線

圖7 加速時間和空阻系數關系曲線

由圖6可以看出電機峰值功率越大，電動汽車的加速時間越短，加速能力越強。從圖7可以看出空阻系數越大，電動汽車的加速時間越長。

3.3 爬坡能力的影響因素分析

電動汽車的爬坡能力由其最大爬坡度進行表征，編寫在三種不同的行車速度下，電機峰值功率和空阻系數對電動汽車爬坡能力影響的MATLAB程序，運行得到該型電動汽車最高車速與電機峰值功率和空阻系數這二者關系仿真曲線，如圖8和圖9所示。

圖8 最大爬坡度和峰值功率關系曲線

圖9 最大爬坡度和空阻系數關系曲線

由圖8可以看出在一定車速下電機峰值功率越大，電動汽車的最大爬坡度越高，其爬坡能力越強；在一定的峰值功率條件下，電動汽車的速度越高，其對應的最大爬坡度越高。由圖9可以看出在一定車速下空阻越大，電動汽車的最大爬坡度越低，其爬坡能力越弱。

4 結論

根據純電動汽車動力學理論，針對具體車型建立電動汽車動力性模型，并利用MATLAB計算分析程序，得到在日常工況下純電動汽車的最高車速、加速時間和最大爬坡度，并繪制電動汽車動力性仿真曲線。最后具體分析電車峰值功率和空阻系數對汽車動力性的影響。通過該模型能夠為純電動汽車動力性設計提供理論依據，同時該模型具有廣泛的適用性，更改相關基礎參數就可以實現不同車型、不同條件下的動力性分析，為今后研究純電動汽車動力性能提供了有價值的參考。