基于实际道路工况的混合动力客车能量管理策略多目标优化
………………………………… 高建平 孙家辉 郗建国 任德轩 （
【摘要】为解决以标准工况为基础优化的混合动力客车（ HEB）能量管理策略在实际道路工况下燃油经济性、排放性能未能充分发挥的问题，采集某市S122路公交车的实际行驶数据，通过主成分分析与模糊聚类相结合的工况构建法解析该公交线路的代表性典型工况，并在该工况下，采用多目标粒子群优化算法对整车控制参数进行优化，对优化后的能量管理策略进行驾驶员在环半实物仿真验证，结果表明，优化后整车燃油消耗量下降8.19%，综合污染物排放量下降8.63%。主题词：混合动力客车 行驶工况 多目标粒子群优化 半实物仿真验证U469.7 A 10.19620/j.cnki.1000-3703.20181273中图分类号： 文献标识码： DOI:
1 前言
Hybrid Electric Bus，HEB）混合动力电动客车（ 含有发动机和电机两种动力源，能量管理策略是分配并调节各动力源功率输出的关键技术，直接影响整车的动力性、经济性和排放性。能量管理策略的多目标优化可在不牺牲动力性的前提下，提高整车的燃油经济性，降低污染物排放量，并有效延长电池的使用寿命。近年来，国内外学者提出了很多针对控制策略的多目标优化方法，主要分为两类，一类是通过加权油耗、排放等将多目标归一化为单目标的优化法[1- 3]，另一类是同时对各目
标并行优化的多目标优化方法[4- 7]。然而，前者的优化结果难以客观反映各子目标的真实情况，且每次优化只1能获得组最优解。目前，针对控制策略的多目标并行优化大多是基于标准典型循环工况的，与车辆实际道路行驶工况存在较大差异，导致车辆在实际路况行驶中的燃油经济性、排放性能不能达到最佳状态。基于上述分析，本文以行驶线路固定的混合动力城MATLAB市公交车为研究对象，在中编写相应程序对采集到的有效行驶数据进行运动学片段划分、主成分分
C-析、模糊均值聚类，从而构建出符合该公交线路特征的实际典型循环工况。同时以整车的燃油经济性和污
染物排放量作为优化目标，基于构建的典型循环工况对控制策略进行多目标并行优化，最后通过半实物仿真试验验证优化后控制策略的可行性。
2 车辆典型行驶工况构建及验证2.1 车辆行驶数据采集及主成分分析
综合考虑该城市道路的实际特点和公交车的工作1特征及运行区域类型等因素，最终选择如图所示的公交车固定行驶线路，将车载终端数据采集设备安装到混3合动力公交车上，以自主行驶法进行了为期个月的数据采集。
为便于数据处理，定义运动学片段为车辆从一个怠速状态开始到下一个怠速状态开始之间的运动过程[8]。对采集的有效行驶数据进行运动学片段划分，共得到9 643
条运动学片段。为了充分表达各运动学片段的特征，本文选取运行距离L、平均速度vm、最大速度vmax、运行时间t、怠速时间ti、加速时间ta、减速时间td、匀速时间tc、最大加速度amax、最小加速度amin、加速段平均加速度aa、减速段平均加速度ad等12
个具有代表性的特征参数，构
1建出了短行程特征参数矩阵M9 643×12，其中的元素如表
12
所示。在选取的 个特征参数中，有些变量间存在关联，并不相互独立。为消除选取参数之间的冗余性，减少特征参数数量，提高后期聚类准确度，对M9 进行主成643×12
2分分析，得到其特征值的贡献率和累积贡献率如表 所
2 4
示。从表 可以看出，前个主成分的累积贡献率为84.97%， 80%，
已超过 能够较好地反映原有信息。在进行主成分分析的同时也得到了主成分得分矩阵S9 643×4， 3
其元素如表 所示。
2.2 模糊C-均值聚类算法
C- C-
模糊 均值算法由普通 均值算法改进形成，是一种柔性的模糊划分，能准确地反映出实际样本的分C-
布情况[9]，因此本文选择模糊 均值聚类对得分矩阵C-
S9 进行聚类。模糊 均值聚类算法是一种基于目643×4标函数的聚类算法，需要求取样本点到各聚类中心的距离，以及样本点的隶属度，具体步骤为： a.
确定类别个数V，初始化隶属度矩阵，选取模糊加权参数m，迭代停止误差Error以及最大迭代步数Mis。b. 初始化聚类中心C，计算初始的距离矩阵。c. 更新聚类中心和隶属度：式中， Cj为第j类聚类中心； Xi为第i个样本； uij为样本Xi dij(j= 1,2,…, V) dik(k= 1,2,…, V)隶属于第j类的隶属度； 、 分别为样本Xi到第j类、第k类中心的距离。d. 更新距离矩阵，并计算目标函数值J： e. 2若达到最大迭代次数 或者连续 次目标函数差值的绝对值小于Error则停止，否则转到步骤c
。
C- 4
采用模糊 均值聚类将前 个主成分的得分矩阵
3 Ⅰ
聚成 类，分别为 类低速工况（平均速度相对较低，怠Ⅱ
速比例较高）、 类中速工况（平均速度适中，怠速比例Ⅲ
适中）、 类高速工况（平均速度较高，怠速比例小）。
2.3 行驶工况合成
1 200 s，通常，构建的典型城市循环工况时长约为可利用各类总运行时间在总体数据中所占的时间比来确定各类运动学片段在最终拟合工况中所占的时间[10]，然后从每类中选取与此类综合特征值相关系数最大的代表性运动学片段，来构建车辆行驶工况。将各类综合特征参数值及运动学片段的特征参数值视为变量GP、的分布，二者的相关系数δ为： δ越大，说明该运动片段与该类越相关，越具有代表性。
Ⅰ 5 Ⅱ完成上述计算后，从 类中取 条低速片段，从
9 Ⅲ 1
类中取 条中速片段，从 类中取 条高速片段，构建出1270 s 46.8 km/h了持续时间为 、最高车速为 的行驶工况DC_S122， 2
如图 所示。
2.4 车辆典型行驶工况验证
选取某公司开发的混合动力公交车进行工况验3证，该车辆的动力系统结构如图 所示，关键部件参数4
如表 所示。该动力系统各工作模式间相互切换的逻4
辑条件如图 所示，其中， Treq为整车需求转矩， Te- off为发动机转矩下限， Te- max为发动机转矩上限， SOC、SOClow分别为电池荷电状态及其下限。
DC_S122
为验证 典型工况的合理性，随机选取该路线上某日公交车的实际行驶数据，并选择基于中国典型DC_S122城市公交工况优化后的控制策略，与基于 工况Interface的控制策略进行对比。通过 接口将在AVL_Cruise MATLAB/中建立的整车动力学模型集成于
Simulink Cruise中，将上述工况加载到整车模型中进行
Simulink CO NOx HC
与 的联合仿真，燃油消耗量和 、 和 排5
放量的仿真结果如表 所示。
5
由表 可知，中国典型城市公交工况下的燃油消耗量、污染物排放量与实际工况相比差异显著，而DC_S122
行驶工况下的燃油消耗量和污染物排放量与实际工况结果相比误差较小，能够代表该路线公交车的真实行驶状况。
3 基于Pareto解集的多目标优化设计3.1 优化参数及约束条件选择
HEB
的优化涉及很多参数，本文在车辆各部件参数已经确定的条件下，以对整车动力性、燃油经济性及排放性能影响较大的控制参数为变量进行优化，优化变
6 0~18 km/h ≤3 s；量如表 所示。约束条件为： 加速时间
0~50 km/h ≤20 s；15 km/h ≥12%；加速时间 最大爬坡度
≤3 km/h实际车速与期望车速误差的绝对值 。
3.2 多目标优化数学模型
HEB
控制策略的优化是典型的多目标数学问题，求解的目的是使相互冲突的多个目标函数，在可行域内Pareto Pareto得到最优解集，用 原理表示，称为 解集。通常，多目标数学模型可表述为：式中， 为多目标函数； 、 为互相冲突的子目标g(x)
函数； 为约束条件；为约束条件个数；为设计变量i j个数； xj、xj, min、xj, max分别设计变量及其下限和上限。依据多目标数学模型，建立以发动机燃油消耗量与污染物排放量为多目标优化的数学模型，定义Qfuel为车f1(x)辆的燃油消耗量，则车辆燃油消耗量函数 的数学模型为： CO NOx HC QCO(x) QNOx(x)定义 、 、 的排放量分别为 、 、QHC(x)， CO NOx HC
由于 的排放量约为 、 两种污染物排放10 f2(x)
量的 倍，因此污染物排放量的综合评价指标 的数学模型为：
3.3 多目标粒子群优化算法
Particle Swarm Optimization，PSO）粒子群优化（ 算法与其他优化算法相比，具有规则简单、容易实现、精度PSO高、收敛快等优点，因此本文选取 对控制参数进行PSO
优化。 算法通过社会认知和自我认知，进行个体速度和位置的更新，并逐代搜索寻优，最终获得全局最优解。速度和位置更新公式分别为：式中， V k 、X k分别为粒子i的第k代速度和位置； V k +1 、i i i ( k+ 1)
X k +1分别为粒子i的第 代速度和位置； w为惯性因i
子； C1、C2为学习因子，通常C1=C2；Pbesti为粒子i的历史最rand( ) 0~1
佳位置； 为随机数生成函数，可生成 之间的随机数； Gbest为种群中的历史最佳位置。多目标粒子群优化方法与求解单目标的优化算法相Pareto比，引入了精英集策略，即选取非劣解“粒子”构成精英集（每次迭代过程中的非劣解集合）。利用多目标粒5子群优化算法对控制参数优化的具体流程如图 所示。
4 优化仿真结果分析及半实物仿真验证4.1 优化结果分析
DC_S122
在 工况下进行粒子群多目标离线优化，仿100 km
真得到 燃油消耗量和综合排放量，仿真结果的Pareto 6 Pareto解集及其收敛曲线如图 所示。将得到的
9
解集与优化前的结果进行比较，得到组最优支配解7
集，如表 所示，其中最后一组为控制策略优化前的燃油消耗量、污染物排放量综合评价指标，以该组数据作为基准。
9 Pareto
与优化前相比，组 最优支配解集下的目标6.8%，值均得到了优化，其中燃油消耗量最多下降了 污6.43%染物排放量综合指标平均降低了 。多目标优化后可得到多组最优解，可以根据关注目标的偏重个性化定制控制策略，本文最终选取燃油经济2
性最优组（第 组）的控制策略进行半实物仿真验证。
4.2半实物仿真验证
为使仿真结果更加接近于实际车辆运行情况，对优化前、后的控制策略进行基于驾驶员在环的半实物仿真dSPACE CANoe验证。该仿真平台主要包括 、 、驾驶员模MATLAB AVL Cruise ControlDesk拟器等硬件和 、 、 等软7
件[11]，仿真平台如图 所示。仿真试验工况是从车辆实28 min际行驶中随机选取的一段时长为 的工况，将Cruise 0.65，整车的电池模型SOC初值设为 仿真试验结8
果如图 所示。
8a
从图 中可以看出，在优化后的能量管理策略下，驾驶员很好地完成了车速跟踪，实际车速与期望车速基
3 km/h本相同，车速跟随误差的最大绝对值控制在 以内，满足整车动力性需求，验证了控制策略的可靠性。8b
由图 可知，在优化后的控制策略下动力电池的SOC比优化前的波动更小，满足混合动力汽车对电量均衡的8c 8d要求，有利于减少电池的损耗。由图 和图 可知，控制策略优化后发动机起动参与系统工作的时间、频次比优化前均有所下降，改善了整车的燃油经济性。
28.94 L
优化前、后百公里燃油消耗量分别为 和26.57 L， 127.12 g/km污染物排放量综合评价指标分别为
116.15 g/km
和 。优化后燃油消耗量和污染物排放量综8.19% 8.63%，合评价指标分别降低 和 表明多目标优化后HEB
的整车控制策略在保证 最佳燃油经济性的基础上，最大程度地降低了污染物的排放量，从而验证了多目标优化后能量管理策略的可行性。
5 结束语
C-本文运用主成分分析和模糊 均值聚类相结合的DC_S122，方法拟合出了某公交线路的代表性循环工况
与中国典型城市公交工况进行仿真对比，结果表明， DC_S122
典型循环工况能够更准确地表征该路线公交DC_S122
车的实际行驶状况。在 工况下，采用多目标粒子群优化算法优化了对整车能量管理策略有重要影响的关键控制参数，燃油经济性和污染物排放性综合评价
Pareto
指标均明显提高，得到的 解集具有良好的收敛性。选取燃油经济性最优组的控制策略进行半实物仿真试验，从而验证了优化后控制策略的可行性，解决了因控制策略与实际行驶工况不匹配导致混合动力客车在实际行驶中燃油经济性、排放性能无法达到最佳状态的问题。
参考文献
[1] , , .王永宽 钱立军 牛礼民 基于免疫算法的四驱插电式混[J]. ,合动力汽车控制策略多目标优化 中国机械工程2017, 28(14): 1683-1689. [2] , , , .李军 朱亚洲 徐杨蛟 等 混合动力汽车动力传动系统控[J]. , 2014(3): 138-141.制策略优化研究 机械设计与制造[3] , .吴光强 陈慧勇 基于遗传算法的混合动力汽车参数多目[J]. , 2009, 31(1): 61-64.标优化 汽车工程[4] .周云山，贾杰锋 基于正交试验设计和多目标遗传算法的HEV [J]. , 2014, 5(4): 324参数优化 汽车安全与节能学报330. [5] , , , . NSGA-II邓涛 林椿松 李亚南 等 采用 算法的混合动力[J]. ,能量管理控制多目标优化方法 西安交通大学学报2015, 49(10): 144-150. [6] . HEV赵亚光 多目标粒子群优化算法研究及在 参数优化[D]. : , 2016.中应用 南京 南京信息工程大学[7] KITAYAMA S, SAIKYO M, NISHIO Y, et al. Torque Control Strategy and Optimization for Fuel Consumption and Emission Reduction in Parallel Hybrid Electric Vehicles[J]. Structural & Multidisciplinary Optimization, 2015, 52(3): 595-611. [8] , , , .胡志远 秦艳 谭丕强 等 基于大样本的上海市乘用车行[J]. ( ), 2015, 43(10):驶工况构建 同济大学学报 自然科学版1523-1527. [9] , . w- [J].张瑞丽 张继福 基于 距离均值的模糊聚类算法 计, 2012, 32(7): 1978-1982.算机应用[10] , , , . k秦大同 詹森 漆正刚 等 基于 均值聚类算法的行驶[J]. ( ), 2016, 46(2):工况构建方法 吉林大学学报 工学版383-389. [11] , , , .高建平 赵金宝 杨松 等 基于驾驶意图的插电式混合[J]. , 2016, 52(24):动力公交车控制策略 机械工程学报107-114. （责任编辑 斛畔） 2018 11 14修改稿收到日期为 年 月 日。