插电式混合动力汽车动力传动系参数优化
王庆年1, 段本明1, 王鹏宇1, 拱印生2, 朱庆林2
1.吉林大学 汽车仿真与控制国家重点实验室,长春130022
2.启明信息技术股份有限公司 电子服务中心,长春130122
通讯作者:王鹏宇(1979-),男,副教授,博士.研究方向:节能与新能源汽车.E-mail:wangpy@jlu.edu.cn

作者简介:王庆年(1952-),男,教授,博士生导师.研究方向:节能与新能源汽车.E-mail:wqn@jlu.edu.cn

摘要

在既定的动力总成功率参数及整车控制策略前提下,为充分挖掘插电式混合动力汽车(PHEV)节能潜力,以满足整车性能为约束条件,以百公里行驶成本为响应,利用最优拉丁超立方设计方法探索了其传动系统所有因子的空间响应特性,辨别了系统关键设计因子,给出了最优局部区域。基于此,建立了系统径向基(RBF)神经网络模型,并充分利用非线性二次规划算法较强的局部优化能力,在上述局部区域内得到了传动系参数全局最优组合以及对应的百公里行驶成本。结果显示:基于近似模型的优化方法精度较高,误差为1.06%;百公里行驶成本降低了9.72%。

关键词: 车辆工程; 插电式混合动力汽车; 百公里行驶成本; 最优拉丁超立方设计; 径向基神经网络模型; 非线性二次规划算法
中图分类号:U469.7 文献标志码:A 文章编号:1671-5497(2017)01-0001-07
Optimization of powertrain transmission parameters of plug-in hybrid electric vehicle
WANG Qing-nian1, DUAN Ben-ming1, WANG Peng-yu1, GONG Yin-sheng2, ZHU Qing-lin2
1.State Key Laboratory of Automotive Simulation and Control,Jilin University,Changchun 130022,China
2.Electric Service Center,Qiming Information Technology Co.,Ltd.,Changchun 130122,China
Abstract

In order to fully explore the energy saving potential of plug-in hybrid electric vehicle with determined power parameters of powertrain and control strategy, and in meeting the vehicle performance constraints, the running cost of 100 km was proposed as the response. The optimal Latin hypercube design was used to investigate the spatial response characteristics of all factors of the transmission system, based on which the key design factors of the system were identified and the optimal local area was confirmed. Combined with making full use of the advantages of strong local optimization ability of nonlinear quadratic programming algorithm, the Radial Basis (RBF) neural network model was established to obtain global optimal solution and the best response of the optimization problem. The model error is 1.06% which shows that the optimization method based on approximation model is of high accuracy. The running cost of 100 km is reduced by 9.72%.

Keyword: vehicle engineering; plug-in hybrid electric vehicle(PHEV); running costs of 100 km; optimal latin hypercube design; RBF neural network; nonlinear quadratic programming algorithm
0 引 言

目前, 针对插电式混合动力汽车(Plug-in hybrid electric vehicle, PHEV)[1, 2]的研究多集中在开发较先进的优化算法, 例如利用动态规划算法、模糊逻辑控制算法、模型预测控制算法等对动力系统能量管理策略进行优化[3]; 利用遗传算法、自适应模拟退火算法等全局探索方法对动力传动系统参数自动离线寻优等[4]。上述研究方法仍存在一些问题:并没有考虑动力传动系统设计参数与目标函数之间的空间响应特性, 缺少对其准确的把握, 使之对工程实际应用的参考价值得不到充分的体现; 常见优化算法中目标函数确定方法是将车辆的动力性指标、经济性指标进行加权转化为单目标优化问题或者直接利用多目标优化算法进行优化[5], 其优化效率较低且评价指标具有一定的人为主观性。

本文以实际项目中某PHEV为研究对象, 以百公里行驶成本作为优化目标[6, 7], 利用最优拉丁超立方设计方法(Optimal Latin hypercube design, OptLHD)对优化参数响应空间特性进行探索, 从全局把握了最优设计区域。建立了传动系参数与百公里行驶成本之间的径向基神经网络(RBF nueral network)近似模型。最后充分利用非线性二次规划算法(NLPQL)较强的局部优化能力在最优设计区域对近似模型寻优, 得到了最优参数设计组合。通过以上研究, 既把握了因子与响应之间的空间特性, 又提高了优化效率。

1 传动系参数优化对象
1.1 PHEV动力系统构型

某插电式混合动力汽车结构如图1所示, 其中, i0i1i2i3分别为主减速系传动比、驱动电机与主减速器之间的传动比、发动机与主减速器之间的传动比、发动机与ISG电机之间的传动比。可以看出是一种典型的混联式构型。驱动电机、发动机可以分别单独驱动或者联合驱动整车; 启动/发电一体机(ISG)可以启动发动机并与发动机组成发动机-发电机组为电池充电; 在制动能量回收模式下, 驱动电机转变为发电机将可回收能量存储到电池中。

图1 某插电式混合动力汽车系统结构示意图Fig.1 Configuration of plug-in hybrid electric vehicle system

1.2 PHEV控制策略的设计

整车控制策略基于以下原则:首先, 将该PHEV工作过程分为电量消耗(Charge depleting, CD)和电量保持(Charge sustaining, CS)两个工作阶段; 其次, 结合该插电式混合动力汽车构型特点, 采用低速纯电动驱动、高速发动机优先参与工作的控制逻辑。具体控制思想如下:

(1)当车速 vv0且SOC≥ SOC0时, 驱动电机单独驱动。

(2)当车速 vv0且SOC< SOC0时, 驱动电机单独驱动, 发动机-发电机组发电。

(3)当车速 v> v00< Treq< Ten_low且SOC≥ SOC0时, 驱动电机单独驱动。

(4)当车速 v> v0Ten_lowTreqTen_h且SOC≥ SOC0时, 发动机单独驱动。

(5)当车速 v> v0TreqTen_h且SOC< SOC0时, 发动机单独驱动并发电。

(6)当车速 v> v0Treq> Ten_h时, 发动机与驱动电机联合驱动。

(7)当制动踏板踩下且 v> vr时, 驱动电机转换为发电机进行制动能量回收。

其中, v0为车速门限值, km/h; SOC0为电池目标荷电状态; Treq为整车需求发动机提供转矩, N· m; Ten_low为发动机最优工作区间下限, N· m; Ten_h为发动机外特性曲线扭矩, N· m; vr为进入制动能量回收模式车速最低值, km/h。

1.3 整车系统的简化模型

(1)发动机以及电机模型

在已知发动机和电机准确的实验数据前提下, 通常将发动机和电机简化为参数插值查表模型以获得其输出特性, 并结合其自身特性对模型进行修正。在一定的发动机转速 we(r/min)(电机转速 wm(r/min))下, 可以通过查表得到发动机转矩 Te(N· m)(电机转矩 Tm(N· m))、发动机功率 Pe(kW)(电机功率 Pm(kW))、发动机燃油消耗率be(g/(kW· h))之间的序列。故发动机的稳态油耗、稳态功率、电机的功率可以表示为:

be=ke·f(Te, we)(1)

Pe=kp·f(Te, we)(2)

Pm=kpm·f(Tm, wm)(3)

式中: kekpkpm为插值修正系数, 由插值方法确定。

(2)电池模型

若已知电池在一定温度 T(℃)下电池的荷电状态SOC、充放电电阻R(Ω )、电压 U(V)电流 I(A), 根据电池的特性, 电池功率可以表示为:

Pb=UIηb-I2Rηb2(4)

式中: ηb为库伦效率; 电池SOC可以表示为:

SOC=SOC0+13600Ct0tIdt(5)

式中: C为电池容量, Ah。

(3)整车动力学模型

根据PHEV的构型, 结合汽车方程可知:汽车行驶阻力转矩Treq 应等于动力总成驱动转矩Tt, 分别表示为:

Treq=mgf+CDAua221.15+σmduadtr(6)

Tt=(Tmi1η1+Tei2η2)i0η0(7)

式中: m为汽车整备质量, kg; g为重力加速度, m/s2; f为轮胎滚动阻力系数; CD为空气阻力系数; A为迎风面积, s2; ua为汽车的行驶速度, km/h; σ为车辆旋转质量的换算系数; r为车轮滚动半径, m; η1η2η0分别为传动比 i1i2i0对应的传递效率。

2 PHEV传动系参数优化设计
2.1 问题描述

(1)优化目标确定及参数选择

PHEV运行时分为电量消耗和电量维持两个阶段, 其燃油经济性应充分考虑以上两个阶段。传统评价方法主要有两种:一种是分别给出两个阶段的燃油经济性, 该种方法并不能直观体现PHEV总的燃油经济性; 另一种是等效燃油消耗量方法, 但等效转化建立在一定的估算效率基础上, 而效率容易受各种环境因素的影响, 故该种方法亦不能准确地反映出整车燃油经济性。提出以百公里行驶成本为评价指标, 将油耗、电耗统一转化为总消耗费用:

f=fe+ffuel=Qe·pe+Qa·pfuel(8)

式中: feffuel分别为百公里电耗和油耗费用, 元; Qe为百公里电能消耗量, kW· h/100km; pe为当前每度电的价格, 元/(kW· h); Qa为百公里燃油消耗量, L/100 km; pfuel为当前每升燃油的价格, 元/L。

基于PHEV构型, 选择 i1i2i3以及 i0作为传动系优化参数。

(2)约束条件确定

整车最大总转矩需要满足最大爬坡度需求:

Tα=rii·ηTFt(9)

Ft=mgfcosα+mgsinα+CDAui221.15(10)

式中: ηT=η1·η0; ii=i0·i1; Ft为轮胎驱动力, N; r为车轮滚动半径, m; ii为爬坡时传动比; ui为爬坡速度, km/h; ηT为爬坡时传动效率; α为最大爬坡度对应的角度, rad。

最高车速要求传动系匹配需满足以下条件:

i0·i10.377·nm_max·rua_maxi0·i20.377·nen_max·rua_max(11)

式中: ua_max为整车最大要求速度, km/h; nm_max为驱动电机最高转速, r/min; nen_max为发动机最高转速, r/min。

由以上分析可知, i1i2i0需要满足式(9)(10)(11), i3需要考虑两者的高效区范围。

图2和图3为已选定发动机的万有特性图和ISG电机效率特性图。由图2~图3可以看出, 发动机和ISG电机高效率区间对应的转速区间大致分别为1300~3500 r/min和6000~13 000 r/min。综合考虑两者最优工作转速区间以及转矩特性, 并结合式(9)~(11)得到各传动比取值范围为:

5.6406i0·i17.8919i0·i24.27480.433i30.55(12)

图2 发动机万有特性曲线Fig.2 Curves of engine fuel consumption characteristics

图3 ISG电机效率特性曲线Fig.3 Efficiency map of ISG

(3)整车参数以及整车性能要求

整备质量m=1600 kg; 空气阻力系数CD=0.35; 迎风面积A=2.1 m2; 车轮滚动阻力系数 f=0.013; 车轮滚动半径r=0.314 m; 最大车速umax=180 km/h; 最大爬坡度imax≥ 30%; 爬坡车速ui=30 km/h; 传动系参数组合初选值为 i0=3.837, 1.87, 0.991, 0.435T, 对应的百公里行驶成本为 p0=14.385元。

2.2 传动系设计空间的初步探索

试验设计方法(Design of experiments, DOE)是数理统计学的一个分支。因其能够辨别系统关键设计因子、分析因子与响应之间的关系和趋势, 在工程和科研中应用非常广泛[8]

最优拉丁超立方设计(Optimal Latin hypercube design, OptLHD)可以使所有的试验点尽量均匀的分布在设计空间, 具有非常好的空间填充性和均衡性[8]。具体做法是在 n维空间中, 将每一维坐标区间 xkmin, xkmax, k1, n分别均匀等分为 m个区间, 每一个小区间记为 xki-1, xki, i1, m在其中随机选取 m个点, 选择过程保证一个因子的每个水平只被研究一次, 并要求均匀分布, 即构成 n维空间, 可以记为 m×nLHD

为了更全面地把握参数设计空间特性, 在试验设计中设置跨度较大的变量上、下限, 选择较大的抽样空间。在Matlab中编写最优拉丁超立方设计算法, 设定抽样空间为500, 生成样本初始空间。样本点分布情况如图4所示。

图4 OptLHD方法初始样本点分布Fig.4 Distribution of initial sample points generated by OptLHD

将以上初始样本点限制在约束范围之内, 得到最终样本点分布如图5所示, 共有90组样本点满足条件。图4和图5中红色代表 i0, 蓝色代表 i1, 黑色代表 i2, 绿色代表 i3

图5 考虑约束的OptLHD方法样本点分布Fig.5 Constraints-considered distribution of sample points generated by OptLHD

由图5可以看出, 最优拉丁超立方设计矩阵能够均匀地覆盖整个设计空间。以90组样本点为自变量输入整车模型, 结合整车控制策略在NEDC工况下进行仿真分析。设计过程中响应变化曲线如图6所示。

图6 OptLHD响应变化过程Fig.6 Changes of response by OptLHD

通过对以上所有满足条件的试验设计点进行试验研究, 可以初步得到上述设计参数组最优响应值为 p'min=13.342元, 对应第62组的试验样本, 即试验设计得到的PHEV最优传动系参数组合为 i'min=3.859, 2.2715, 0.8056, 0.54531T, 与初始样本点相比节约成本7.25%。

为进一步确定传动比优化空间, 需要了解因子的主效应, 其绘制方法如下:当某因子从较低水平增加到最大水平时, 用每个水平和其他因子的所有可能组合得到的响应的平均值分别作为横坐标和纵坐标可以得到该因子的主效应图。利用此方法得到所有因子主效应图, 如图7所示。

图7 因子对响应的主效应Fig.7 Main effects graph of factors to response

由图7可以看出因子 i0i1i2对响应的影响较为明显, 为主要影响因子, i3对响应的影响最小。

为进一步寻优得到全局最佳参数组合, 需要缩小可优化空间。缩小因子优化空间的原则是:在该空间内, 响应值应保持在较低的水平之内。在这里以主效应响应不超过最优样本点响应值 p'min为约束, 因子全局较优水平空间最终确定如下:

1.5i03.512.02i14.250.5i21.810.433i30.55(13)

2.3 基于RBF网络模型的传动系参数优化

为进一步优化得到该PHEV动力系统的最佳传动比设计组合, 需先建立各传动比与百公里行驶成本之间的数学模型。插电式混合动力汽车的动力系统集电器、机械、化学和电控技术于一体, 是一种多输入多输出的高度非线性系统。各系统之间存在耦合关系, 难以严格按照解析法要求建立理论模型。近似模型能够起到减少仿真程序调用的时间、提高优化效率的作用, 有利于更快地收敛到全局最优点[8]。径向基神经网络模型较其他近似模型有较强的非线性逼近能力, 无需数学假设, 对函数的逼近能力是最优的, 可以以任意精度逼近任意连续函数[9, 10]。RBF网络模型如图8所示。

图8 径向基神经网络模型示意图Fig.8 Model of RBF nueral network

RBF网络中除输入输出层之外仅有一个隐层, 该层中的转换函数一般为是局部响应的高斯函数, 隐层中的神经元越多, 逼近越精确[11]

RBF网络输入为:

P=(p0, p1, p2, p3)T=(i0, i1, i2, i3)T(14)

在本问题中, 输出层只有一个输出。故RBF网络输出层节点函数为:

hj(p)=φp-cjrj, j=1, 2, , m(16)

p-cj=(p-cj)T(p-cj)(17)

式中: wj为隐层与输出层之间第 j个隐层节点的权值; b为输出层阀值; m为隐层节点数; hj(p)为输入为 p时第 j个神经元的输出函数; cj为第 j个隐层节点中心; rj为第 j个隐层节点的宽度; p-cj为输入 p与输出节点 cj中心之间的欧几里德距离; φ为隐层径向基函数, 在这里利用高斯函数作为径向基函数, 即:

φ(r)=exp-r22σ2(18)

网络输出 T'(p)与期望输出 T(p)之间的误差评价指标为:

e=12i=1N(Ti(p)-T'i(p))2(19)

式中: N为网络训练输入组数。

本文采用Matlab编程, 以 e=0为目标, 以DOE试验设计样本点训练RBF神经网络, 建立了动力传动系统中各传动比与百公里行驶成本之间的神经网络近似模型。其训练过程的误差变化过程如图9所示。

图9 径向基神经网络训练过程误差变画曲线Fig.9 MSE curve during training process of RBF nueral network

从误差变化曲线可以看出, 最终拟合均方误差为1.5× 10-24, 与目标均方误差0非常接近, 说明该神经网络模型精度非常高, 可以作为传动系参数优化模型。图10为RBF网络模型部分三维特性图。

图10 RBF网络模型对输入(i0-i1)的响应Fig.10 Response to (i0-i1) of RBF nueral network

图10表示RBF网络模型中传动系参数 i2i3保持不变的前提下, 模型响应随输入 i0i1的变化而变化的规律。仔细观察可以看出当 i0i1分别在[3, 4]和[2, 3]区间时响应较小, 这与之前的DOE试验设计结果一致。

2.4 近似模型的数值优化

数值优化算法能够利用梯度、导数等数学特性非常高效地解决优化问题, 但是非常依赖初始设计点, 如果初始点选择不当会造成最终解落入局部解[8]。参考前文提供的非常准确的全局较优区域和可优化初始点, 利用数值优化方法中的非线性二次规划算法(NLPQL), 在Matlab中编程实现了基于RBF网络模型的优化方法。迭代过程如图11所示。

图11 非线性二次规划算法迭代过程Fig.11 Iterative process during NLPQL

从数值优化结果可以看出最优传动比组合为: imin=3.0295, 2.7518, 0.95155, 0.50282T, 对应的响应值为 pmin=13.126元。

将以上最优传动系参数组合作为原整车动力学模型的输入, 计算得到百公里行驶成本为12.9874元, 与RBF网络预测结果有1.06%的误差, 在允许范围之内。故通过基于RBF网络预测模型的数值优化方法得到的最优传动系参数组合, 可以节约成本。

p0-pminp0×100%=9.72%(20)

将以上最优速比工程化得到最终传动比组合为 imin=3.03, 2.75, 0.95, 0.5029T经计算验证, 工程化前后细小差异不会严重影响整车的性能优化目标, 可以忽略不计。

3 结束语

本文围绕某插电式混合动力汽车传动系参数匹配优化环节, 为避免人为主观评价因素, 提出以百公里行驶成本为系统响应。通过最优拉丁超立方设计探索了因子的空间响应特性、辨别了系统关键设计因子, 并定位了全局最优区域。为提高优化效率, 建立了误差较小的系统径向基神经网络近似模型。基于此, 利用局部优化能力较强的非线性二次规划算法在已定位的全局最优参数设计区域寻优, 得到了全局最优传动系参数组合。结果显示:除发动机与ISG电机之间的传动比 i3之外的所有传动系参数对系统响应分布影响均较大, 为系统关键设计因子; 全局最优传动系参数组合下的百公里行驶成本较初始值降低9.72%; 经验证近似模型与原系统模型相比误差为1.06%, 证明了本文方法的可行性。本文将试验设计方法、近似模型方法、数值优化算法相结合, 不仅提高了混合动力汽车传动系参数设计效率, 而且建立了因子与响应之间的定量关系, 为工程实际传动系参数设计提供了理论指导。

The authors have declared that no competing interests exist.

参考文献
[1] Golbuff S. Optimization of a plug-in hybrid electric vehicle[D]. Atlanta: Georgia Institute of Technology, 2006. [本文引用:1]
[2] Li L, Zhang Y H, Yang C, et al. Hybrid genetic algorithm-based optimization of powertrain and control parameters of plug-in hybrid electric bus[J]. Journal of the Fanklin Institute-Engineering and Applied Mathematics, 2014, 352(3): 776-801. [本文引用:1]
[3] Sun Y, Chen Z, Yan B J, et al. A learning method for energy optimization of the plug-in hybrid electric bus[J]. Science China Technological Sciences, 2015, 58(7): 1242-1249. [本文引用:1]
[4] 张抗抗, 徐梁飞, 华剑锋, . 基于多目标优化的纯电动车动力系统参数匹配方法[J]. 汽车工程, 2015, 37(7): 757-765.
Zhang Kang-kang, Xu Liang-fei, Hua Jian-feng, et al. A parameter matching method for the power train of battery electric vehicle based on multi-objective optimization[J]. Automotive Engineering, 2015, 37(7): 757-765. [本文引用:1]
[5] 尹安东, 赵韩, 杨亚娟, . 多目标遗传算法的混合动力传动系参数优化[J]. 中国机械工程. 2013, 24(4): 552-556.
Yin An-dong, Zhao Han, Yang Ya-juan, et al. Parametric optimization of hybrid powertrain based on multi-objective genetic algorithm[J]. China Mechanical Engineering, 2013, 24(4): 552-556. [本文引用:1]
[6] 王庆年, 韩彪, 王鹏宇, . 基于成本优化的插电式混合动力参数匹配[J]. 吉林大学学报: 工学版, 2016, 46(2): 340-347.
Wang Qing-nian, Han Biao, Wang Peng-yu, et al. Parameter matching for plug-in hybrid electric vehicle based on the cost optimization[J]. Journal of Jilin University Engineering and Technology Edition, 2016, 46(2): 340-347. [本文引用:1]
[7] 吴为理, 兰凤崇, 陈吉清, . Plug-in混合动力汽车动力系统参数设计成本效益分析[J]. 机械科学与技术, 2013, 32(3) : 383-387.
Wu Wei-li, Lan Feng-chong, Chen Ji-qing, et al. Cost-benefit analysis of parameter design of powertrain for plug-in hybrid electric vehicle[J]. Mechanical Science and Technology for Aerospace Engineering, 2013, 32(3): 383-387. [本文引用:1]
[8] 赖宇阳. Isight参数优化理论与实例详解[M]. 北京: 北京航空航天大学出版社, 2012. [本文引用:4]
[9] Chang W Y. Estimation of the state of charge for a LFP battery using a hybrid method that combines a RBF neural network, an OLS algorithm and AGA[J]. International Journal of Electrical Power & Energy Systems, 2013, 53: 603-611. [本文引用:1]
[10] 张庆洪. 基于径向基神经网络的汽车油耗软测量及监测系统研究[D]. 长沙: 湖南大学机械与运载工程学院, 2014.
Zhang Qing-hong. Soft measurement and monitoring system of vehicle fuel consumption based on the RBF neural network[D]. Changsha: College of Mechanical and Vehicle Engineering, Hunan University, 2014. [本文引用:1]
[11] 张万兴. 电动汽车动力电池剩余容量和续驶里程预测研究[D]. 合肥: 合肥工业大学汽车与交通工程学院, 2012.
Zhang Wan-xing. Forecasting for battery remaining capacity and driving ranges of electric vehicle[D]. Hefei: College of Automotive and Traffic Engineering, Hefei University of Technology, 2012. [本文引用:1]