程雷 秦東晨 王耀凱 王迎佳

（鄭州大學(xué)，鄭州 450001）

基于Modelica的純電動(dòng)客車建模仿真研究＊

程雷 秦東晨 王耀凱 王迎佳

（鄭州大學(xué)，鄭州 450001）

為實(shí)現(xiàn)純電動(dòng)客車多領(lǐng)域統(tǒng)一建模與仿真分析，采用模塊化劃分方法對(duì)純電動(dòng)客車進(jìn)行了結(jié)構(gòu)劃分，利用面向?qū)ο蟮亩囝I(lǐng)域物理系統(tǒng)一建模語(yǔ)言Modelica構(gòu)建了純電動(dòng)客車動(dòng)力系統(tǒng)主要部件模型、底盤(pán)總成核心機(jī)械部件模型及整車控制模型。參照某城市公交車設(shè)置仿真模型參數(shù)，對(duì)純電動(dòng)客車多領(lǐng)域模型的綜合性能進(jìn)行了標(biāo)準(zhǔn)工況下的跟隨性仿真與分析。結(jié)果表明，所建立的純電動(dòng)客車模型具有良好的跟隨性，且動(dòng)力性和經(jīng)濟(jì)性與參照公交車基本一致。

1 前言

純電動(dòng)汽車是一個(gè)涉及機(jī)械、電子、液壓、控制等多個(gè)領(lǐng)域的復(fù)雜物理系統(tǒng)[1]，在其研究與開(kāi)發(fā)過(guò)程中，廣泛采用Advisor、Cruise、PAST和HEVC等軟件進(jìn)行建模與仿真，但這些軟件只側(cè)重于單個(gè)部件或單一系統(tǒng)的建模與仿真，很難實(shí)現(xiàn)多領(lǐng)域統(tǒng)一建模與仿真，因此在對(duì)整車進(jìn)行建模與仿真時(shí)，大多采取多個(gè)軟件建模后再進(jìn)行聯(lián)合仿真的方法，這將導(dǎo)致仿真結(jié)果與實(shí)際結(jié)果偏差較大[2]。目前，多領(lǐng)域統(tǒng)一建模與仿真技術(shù)得到了快速發(fā)展，如逐漸發(fā)展成熟的多領(lǐng)域統(tǒng)一建模語(yǔ)言Modelica已應(yīng)用于航空、航天、船舶、車輛等各行各業(yè)[3]，基于Mod?elica建立的模型庫(kù)已經(jīng)囊括了電學(xué)、車輛動(dòng)力學(xué)、熱動(dòng)力學(xué)以及模糊控制等諸多領(lǐng)域[4]，為實(shí)現(xiàn)純電動(dòng)汽車的統(tǒng)一建模與仿真提供了基礎(chǔ)。

本文以純電動(dòng)客車為研究對(duì)象，采用Modelica建模語(yǔ)言分別建立純電動(dòng)客車的動(dòng)力電池、永磁同步電機(jī)、機(jī)械系統(tǒng)部件、控制器等各系統(tǒng)模型及整車仿真模型，并參照某城市純電動(dòng)公交車的實(shí)際運(yùn)營(yíng)情況設(shè)置仿真參數(shù)，在MWorks仿真平臺(tái)中對(duì)所建整車模型進(jìn)行了標(biāo)準(zhǔn)試驗(yàn)工況及特定城市工況下的跟隨性仿真分析，通過(guò)仿真結(jié)果驗(yàn)證了所建模型的正確性。

2 純電動(dòng)客車結(jié)構(gòu)及模塊劃分

2.1 結(jié)構(gòu)分析

純電動(dòng)客車以電池為動(dòng)力系統(tǒng)的能量來(lái)源，采用電動(dòng)機(jī)進(jìn)行動(dòng)力驅(qū)動(dòng)，通過(guò)電池管理系統(tǒng)（（Battery Man?agement System，BMS）和電機(jī)控制器實(shí)現(xiàn)電能與機(jī)械能的能量轉(zhuǎn)化和信息傳遞[5]，以滿足整車的動(dòng)力需求。純電動(dòng)客車的懸架、制動(dòng)系統(tǒng)等底盤(pán)機(jī)械部件與傳統(tǒng)客車基本相同；對(duì)于純電動(dòng)客車的整車控制，是通過(guò)CAN控制總線及各通訊串口來(lái)實(shí)現(xiàn)各子系統(tǒng)信息的交換與控制。純電動(dòng)客車結(jié)構(gòu)（簡(jiǎn)化）如圖1所示。

圖1 純電動(dòng)客車結(jié)構(gòu)（簡(jiǎn)化）示意

2.2 模塊劃分

依據(jù)純電動(dòng)客車各系統(tǒng)所屬學(xué)科領(lǐng)域的不同，采用模塊化劃分的方法將整車劃分為動(dòng)力系統(tǒng)模塊、機(jī)械部件模塊、控制模塊及其它附件模塊，如圖2所示。其中，動(dòng)力系統(tǒng)模塊包括動(dòng)力電池、驅(qū)動(dòng)電動(dòng)機(jī)；機(jī)械部件模塊包括輪胎、懸架、制動(dòng)系等底盤(pán)總成部件；控制模塊主要包括BMS、電機(jī)控制系統(tǒng)等；附件模塊包括道路、環(huán)境等。

圖2 電動(dòng)客車模塊劃分結(jié)果

3 動(dòng)力系統(tǒng)模型

3.1 電池模型

目前，電池模型主要包括等效電路模型、電化學(xué)模型、神經(jīng)網(wǎng)絡(luò)模型等，其中等效電路模型以電池充、放電的工作原理為基礎(chǔ)，能夠準(zhǔn)確地模擬出電池工作時(shí)的基本特性?，F(xiàn)有的等效電路模型如Thevenin電池模型、PNGV電池模型、GNL電池模型中，PNGV電池模型更貼近實(shí)際[6]，其電路方程推導(dǎo)如下：

荷電狀態(tài)方程為：

式中，SOC0為電池SOC的初始值；C0為電池額定電量。

考慮到電池放電過(guò)程中的極化現(xiàn)象、熱效應(yīng)等對(duì)電池電壓、電流及SOC值的影響，基于電池的PNGV電路方程和Modelica的電學(xué)基礎(chǔ)庫(kù)Electrical中的電器組件，建立如圖3所示的電池單體模型。

圖3 電池單體的等效電路模型

在該電池單體模型中，各可變組件的參數(shù)值通過(guò)參數(shù)表進(jìn)行參數(shù)標(biāo)識(shí)，電池的SOC值通過(guò)對(duì)電路模型中電流數(shù)值積分計(jì)算獲得。在建立了電池單體的PNGV等效電路模型后，通過(guò)Modelica標(biāo)準(zhǔn)庫(kù)中Gain組件進(jìn)行電池單體的串/并聯(lián)構(gòu)成電池組，并建立控制總線（Control?Bus）接口的子控制總線（BatteryControlBus），進(jìn)行封裝后得到具有基本管理系統(tǒng)的完整電池箱。

3.2 電動(dòng)機(jī)模型

目前，純電動(dòng)客車所使用的電動(dòng)機(jī)主要有異步電機(jī)、永磁同步電機(jī)和開(kāi)關(guān)磁阻電機(jī)等3種類型[7]，其中永磁同步電機(jī)因具有較高的能量密度和效率以及較寬的調(diào)速范圍而被廣泛應(yīng)用[8]，所以本文采用Modelica標(biāo)準(zhǔn)庫(kù)中的永磁同步電機(jī)作為驅(qū)動(dòng)電機(jī)模型本體，并根據(jù)永磁同步電動(dòng)機(jī)的逆變器、控制器等部件的物理結(jié)構(gòu)，進(jìn)一步構(gòu)建出電動(dòng)機(jī)的DC/AC逆變器模型、矢量控制（FOC）模型等，得到永磁同步電機(jī)及控制器模型，如圖4所示。

4 機(jī)械部件與控制模型

4.1 輪胎模型

輪胎模型采用被廣泛認(rèn)可的“魔術(shù)公式（MagicFor?mula）”輪胎特性模型，也就是荷蘭Delft理工大學(xué)H B Pacejka教授所描述的HBPacejka輪胎模型[9]，該模型利用一套形式相同的函數(shù)方程描述輪胎在不同行駛狀態(tài)下的縱向力、側(cè)向力以及回正力矩，其方程如下：

圖4 永磁同步電機(jī)及控制器模型

式（4）中，對(duì)應(yīng)于輪胎在行駛狀態(tài)下受到的縱向力、側(cè)向力以及回正力矩等不同情況，Y()X可分別代表側(cè)向力、縱向力與回正力矩，X則可以代表輪胎的縱向滑移率或側(cè)偏角；B為剛度因子；C為y曲線形狀因子；D為巔因子。表示y曲線的最大值；E為y曲線的曲率因子；SV、Sh分別為車輪垂直方向和水平方向的偏移量。

根據(jù)“魔術(shù)公式”所描述的輪胎特性方程，將組件Revolute（轉(zhuǎn)動(dòng)副）、BodyShape（‘體’形狀）、WorldForce（相對(duì)力）、Torque（絕對(duì)力矩）、WorldTorque（相對(duì)力矩）等通過(guò)方程建立關(guān)系，構(gòu)建輪胎模型如圖5所示。

圖5 輪胎模型

4.2 懸架模型

在建立懸架模型時(shí)，依據(jù)懸架系統(tǒng)結(jié)構(gòu)的拓?fù)潢P(guān)系，采用Modelica的機(jī)械基礎(chǔ)庫(kù)Mechanics中的基本組件，按照懸架結(jié)構(gòu)的層次關(guān)系[10]先建立雙橫臂獨(dú)立懸架1/2的雙橫臂獨(dú)立懸架模型，然后根據(jù)懸架在縱向平面上的對(duì)稱性，通過(guò)坐標(biāo)對(duì)稱的方法得到整個(gè)雙橫臂懸架模型并進(jìn)行編譯和求解，所建懸架模型如圖6所示。

圖6 懸架多體模型

4.3 制動(dòng)器模型

客車制動(dòng)方式包括盤(pán)式制動(dòng)和鼓式制動(dòng)兩種，由于鼓式制動(dòng)器抗熱衰退性能較差，會(huì)影響純電動(dòng)客車的再生制動(dòng)功能，所以純電動(dòng)客車中多使用盤(pán)式制動(dòng)。根據(jù)盤(pán)式制動(dòng)器的物理結(jié)構(gòu)和制動(dòng)原理，利用Modelica標(biāo)準(zhǔn)庫(kù)中的基本組件Brake、FlangeWithBear?ing、ConstantTorque、Mounting1D等構(gòu)建盤(pán)式制動(dòng)器模型，如圖7所示。

圖7 盤(pán)式制動(dòng)器模型

4.4 底盤(pán)模型

在建立底盤(pán)的Modelica模型時(shí)，首先需要將所建立的懸架模型、制動(dòng)模型和輪胎模型之間的接口定義為通訊接口[11]，然后再基于純電動(dòng)客車的動(dòng)力學(xué)特性，跟據(jù)定義的通訊接口進(jìn)行懸架模型、制動(dòng)模型和輪胎模型的集成和封裝。所建立的純電動(dòng)客車底盤(pán)模型如圖8所示。

圖8 底盤(pán)模型

4.5 控制模型

整車控制模型主要包括仿真工況曲線數(shù)據(jù)讀取、電池的充/放電控制、電動(dòng)機(jī)轉(zhuǎn)速/轉(zhuǎn)矩控制、行駛與制動(dòng)狀態(tài)控制等。其中制動(dòng)踏板和油門(mén)踏板的控制是通過(guò)駕駛員對(duì)制動(dòng)力矩和電機(jī)轉(zhuǎn)速/轉(zhuǎn)矩信號(hào)的控制實(shí)現(xiàn)的，由此基于Modelica語(yǔ)言建立的整車控制模型如圖9所示。

圖9 整車控制模型

根據(jù)純電動(dòng)客車的結(jié)構(gòu)特點(diǎn)及各物理系統(tǒng)的耦合特性，將建立的動(dòng)力系統(tǒng)模型、機(jī)械部件模型、控制模型以及其它附件模型通過(guò)建模時(shí)定義好的接口進(jìn)行耦合連接，得到完整的純電動(dòng)客車多領(lǐng)域仿真模型，如圖10所示。

圖10 整車仿真模型

5 仿真與分析

5.1 仿真參數(shù)

參照在某城市運(yùn)營(yíng)的純電動(dòng)公交車的實(shí)車參數(shù)，設(shè)置整車仿真模型參數(shù)，如表1所列。

表1 仿真參數(shù)設(shè)置

5.2 仿真結(jié)果分析

采用NEDC（New European Drive Cycle）標(biāo)準(zhǔn)循環(huán)工況和特定城市工況對(duì)整車模型進(jìn)行跟隨性仿真。NEDC循環(huán)工況包括4個(gè)市區(qū)循環(huán)工況（1部）和1個(gè)市郊運(yùn)行工況（2部），圖11為NEDC循環(huán)工況仿真結(jié)果。由圖11可看出，市區(qū)循環(huán)中，平均車速為19 km/h，最高車速為50 km/h，每個(gè)市區(qū)循環(huán)的有效行駛時(shí)間是195 s，行程為1.013 km；市郊循環(huán)中，有效行駛時(shí)間為400 s，平均車速為62.6 km/h，行程為6.955 km。整個(gè)NEDC循環(huán)工況仿真時(shí)間為1 180 s，行程為9.3 km，最高車速為92.41 km/h，仿真結(jié)果表明，整車模型的跟隨性較好，其中市郊循環(huán)時(shí)最高車速偏差是由于設(shè)計(jì)時(shí)最高車速限制導(dǎo)致的。

圖11 NEDC工況跟蹤性仿真結(jié)果

圖12為特定城市工況的仿真結(jié)果。由圖12可看出，在269.20 s、300.40 s和368.08 s三個(gè)時(shí)刻，仿真車速依次為42.36 km/h、21.46 km/h和30.52 km/h，對(duì)應(yīng)的特定城市工況下的車速分別為61.53 km/h、44.14 km/h和47.69 km/h，車速偏差較大，其中300.40 s時(shí)刻時(shí)車速偏差最大，為22.68 km/h，這是因?yàn)榭蛙囐|(zhì)量較大且在較短時(shí)間內(nèi)啟停過(guò)于頻繁導(dǎo)致的，但在隨后的304.40 s時(shí)刻，仿真速度曲線與城市工況曲線接近重合；在其余的仿真時(shí)間（269.20 s、368.08 s時(shí)刻）車速偏差均小于3 km/h，偏差比率基本在5%以內(nèi)，由此表明所建客車模型的整體跟隨性較好。

圖12 特定城市工況仿真結(jié)果

圖13為特定城市工況下純電動(dòng)客車電池電流及電壓仿真曲線。由圖13可看出，隨著車速及加速度的增大，放電電流逐漸增大，電池電壓降低；在制動(dòng)時(shí)，電池的放電電流變?yōu)樨?fù)值（圖13a），即電池處于充電狀態(tài)，此時(shí)電機(jī)轉(zhuǎn)為發(fā)電狀態(tài)，電壓曲線有升高趨勢(shì)（圖13b），表明電池能量消耗速率降低，說(shuō)明所建立客車模型的經(jīng)濟(jì)性良好。

圖13 電池電流、電壓仿真曲線

圖14為特定城市工況下行駛里程仿真曲線。由圖14可看出，仿真所行駛路程為8.3 km，最高車速為82.3 km/h（城市公交車行駛的最高車速一般不超過(guò)80 km/h），說(shuō)明所建立客車模型動(dòng)力性較好。

圖14 行駛里程仿真曲線

圖15為電池SOC仿真曲線。由圖15可看出，仿真時(shí)電池的SOC值從理想狀態(tài)1.0下降到0.871，耗能為3.504 kW·h，可計(jì)算得能耗為465.48 Wh/km，與所參照的純電動(dòng)公交車實(shí)測(cè)的能耗平均值480 Wh/km相比，誤差為3%，經(jīng)濟(jì)性與實(shí)車相近，說(shuō)明所構(gòu)建的多領(lǐng)域模型是正確的。

圖15 電池SOC仿真曲線

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

為實(shí)現(xiàn)純電動(dòng)客車多領(lǐng)域物理系統(tǒng)的統(tǒng)一建模與仿真，本文分析了純電動(dòng)客車與傳統(tǒng)熱動(dòng)力客車的結(jié)構(gòu)差異，并采用模塊化劃分的方法將純電動(dòng)客車劃分為動(dòng)力系統(tǒng)模塊、機(jī)械部件模塊、控制模塊以及其它附件模塊，同時(shí)基于多領(lǐng)域統(tǒng)一建模語(yǔ)言Modelica分別建立了純電動(dòng)客車各系統(tǒng)模型，搭建了純電動(dòng)客車的仿真模型。參照某城市純電動(dòng)公交車實(shí)際運(yùn)行參數(shù)，對(duì)所建模型進(jìn)行了NEDC工況和特定城市公交工況下的跟隨性仿真分析，仿真結(jié)果表明，所建模型的動(dòng)力性和經(jīng)濟(jì)性與實(shí)際公交車相符，驗(yàn)證了所建模型的正確性。

1 王俊蘭,吳義忠,熊會(huì)元.純電動(dòng)汽車整車建模仿真研究.計(jì)算機(jī)仿真,2015（10）:190～195.

2 王淑慧.汽車饋能式懸架性能仿真分析：[學(xué)位論文].武漢：武漢理工大學(xué),2013.

3 王祺.面向Modelica的多體系統(tǒng)笛卡爾方法建模：[學(xué)位論文].武漢：華中科技大學(xué),2011.

4 田麗娟.純電動(dòng)轎車多領(lǐng)域統(tǒng)一建模與操控穩(wěn)定性仿真分析：[學(xué)位論文].廣州：廣東工業(yè)大學(xué),2014.

5 Clemens Satzger,Jonathan Brembeck.Framework for the Evaluation of Wheel Torque Blending Algorithms.IFAC Proceedings Volumes,2013,4621.

6 李思.電動(dòng)汽車鋰離子電池等效電路模型的參數(shù)辨識(shí)研究：[學(xué)位論文].北京：北京理工大學(xué),2015.

7 戴彥.電動(dòng)汽車永磁同步電機(jī)直接轉(zhuǎn)矩控制策略研究.機(jī)械設(shè)計(jì)與制造,2016（4）:102～106.

8 裴素榮.面向電動(dòng)汽車的永磁同步電機(jī)驅(qū)動(dòng)控制方法研究：[學(xué)位論文].東北大學(xué),2011.

9 王先云,劉艷華,馬斌,劉志敏,劉立國(guó).基于魔術(shù)公式的輪胎特征函數(shù)計(jì)算方法.科學(xué)技術(shù)與工程,2016（7）:265～269.

10 劉煒,吳義忠,陳立平,熊會(huì)元.基于Modelica的電動(dòng)汽車懸架系統(tǒng)建模與仿真分析.汽車技術(shù),2014（9）:43～47.

11 楊俊杰.基于Modelica的純電動(dòng)客車動(dòng)力系統(tǒng)建模與仿真研究：[學(xué)位論文].鄭州大學(xué),2016.、

（責(zé)任編輯 文 楫）

修改稿收到日期為2017年6月6日。

Modeling and Simulation of Pure Electric Bus Based on Modelica

Cheng Lei,Qin Dongchen,Wang Yaokai,Wang Yingjia
（Zhengzhou University,Zhengzhou 450001）

In order to realize multi-domain unified modeling and simulation analysis of pure electric buses,structure of the pure electric bus was divided by the modularization method.Modelica,the unified modeling language of multidomain physics for object oriented,was used to build the model of main power system parts,model of other core mechanical parts of the chassis assembly and the control model of the vehicle.Following performance simulation and analysis of the synthetic performance of the multi-domain model of pure electric buses were carried out with reference of simulation model parameters of a city bus.The results show that the model of pure electric bus has good following performance,and the power and economy are basically identical with the reference bus.

Pure electric bus,Module division,Multi-domain modeling,Modelica

純電動(dòng)客車 模塊劃分 多領(lǐng)域建模 Modelica

U469.7 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：1000-3703（2017）08-0043-06

國(guó)家重大科技成果轉(zhuǎn)化項(xiàng)目（財(cái)建[2012]258）、2014年河南省重點(diǎn)科技攻關(guān)項(xiàng)目（142102210103）和鄭州市科技領(lǐng)軍人才項(xiàng)目（10LJRC188）。