CN108297858A - 一种行星混联混合动力汽车发动机需求功率计算方法 - Google Patents

Abstract

本发明公开了一种行星混联混合动力汽车发动机需求功率计算方法，涉及汽车技术领域。发动机需求功率按用途分为超级电容需求充电功率、动力传输过程中的损失功率以及主减速器输入端整车驱动需求功率。该方法重点包括提出基于超级电容SOC实际值低于控制目标程度及SOC变化速率双反馈计算超级电容充电损失功率，以保证SOC在维持均衡、变化平缓的目的；将损失功率细致划分为电路径损失功率和机械路径损失功率并在两电机不同工作状态下精细化计算。基于精细化考虑损失和合理估算超级电容需求充电功率的计算结果能够准确反映动力系统的需求，应用在控制策略中能够提高其可靠性和响应速度，保证汽车性能。

Description

一种行星混联混合动力汽车发动机需求功率计算方法

技术领域

本发明属于混合动力汽车技术领域，特别涉及一种行星混联混合动力汽车发动机需求功率计算方法。

背景技术

随着能源与环境问题日益突出，节能与新能源汽车成为汽车行业发展的方向。纯电动汽车电池较低的能量密度，燃料电池汽车昂贵的制造、使用和维护成本以及复杂的配套设施限制了其各自发展，混合动力汽车作为向新能源汽车过渡的方案，既可保证动力性，又能改善发动机工作状况、减少有害物排放、节约燃油，经过多年的发展技术已经成熟可靠。油电混合动力汽车将电机和发动机结合在一起，常用的形式有串联式、并联式和混联式，其中混联式包括开关混联和行星混联。行星混联式混合动力汽车的发动机和两电机通过行星齿轮机构连接并向外输出动力，该形式可以通过电机调节发动机的转速和转矩，实现发动机转速转矩的解耦，使发动机在不同的工况需求下都能工作在最节油的状态。因此，行星混联油电混合动力汽车是目前各大汽车公司发展的优选车型。

现有技术中的行星混联式混合动力汽车的动力耦合系统可分为单行星排式和多行星排式，其中单行星排式耦合系统结构简单、成本较低、控制容易，因而获得了较为广泛的应用。针对行星混联式混合动力汽车的发动机、电动机与发电机三个动力源分别固定连接在行星排的行星架、齿圈与太阳轮上，齿圈同时用于向主减速器输出动力。针对行星混联构型的特点，在制定动力源能量分配控制策略时需要考虑复杂的耦合和解耦关系。现有控制策略以发动机最优为基本思想，非插电式(使用超级电容)则主要依靠发动机驱动车辆，超级电容、两电机起到调节发动机工作点的作用，即尽量让发动机承担整车驱动需求功率，同时负责给超级电容充电，因此对发动机需求功率的正确计算对于系统的快速响应和动力性的发挥以及超级电容电量的平衡起到至关重要的作用。

中国专利公开号为CN106894901A，公开日为2017-06-27，公开了一种发动机功率控制方法及系统，该方法仅适用传统内燃机汽车，直接根据加速踏板信号和车速解析发动机需求功率，较为简单；中国专利公开号为CN104343552A，公开日为2015-02-11，公开了一种发动机控制系统及方法，该方法基于对车辆载荷的识别，目的在于通过控制发动机功率在保证动力性的前提下实现较低的燃油消耗，同样只适用于传统内燃机汽车；而在多动力源且耦合关系较复杂的行星混联混合动力汽车领域则尚无相关发明。

发明内容

为克服现有技术存在的不足，本发明提供一种行星混联混合动力汽车发动机需求功率计算方法，充分考虑发动机发出功率在电路径和机械路径的损失，并提出基于超级电容SOC当前值低于控制目标程度及SOC变化速率双反馈的超级电容需求充电功率的计算，使得所求发动机需求功率接近系统真实需求，在提高控制系统响应速度和动力性发挥程度的同时保证超级电容电量的平衡。

为实现上述目的，根据本发明实施例的一种行星混联混合动力汽车发动机需求功率计算方法，将发动机需求功率分为超级电容需求充电功率、动力传输过程中的损失功率以及主减速器输入端整车驱动需求功率，包括以下内容：

I.超级电容需求充电功率基于超级电容SOC当前值低于超级电容SOC目标控制值程度及超级电容SOC变化速率双反馈计算，具体包括以下步骤：

①参数标定：

基于超级电容SOC平衡原则和寿命要求设定基础充电功率P_{chg_bas}。、超级电容SOC目标控制值SOC_control及SOC维持下限SOC_low；

超级电容SOC当前值SOC_real低于超级电容SOC目标控制值SOC_control时即开始充电，根据不同的SOC_real低于SOC_control的程度ΔSOC和不同的超级电容SOC变化速率δSOC，使用不同的充电倍率C；

当超级电容SOC处于下降段时，制定充电倍率C关于ΔSOC和δSOC的二维插值矩阵M1，该矩阵满足充电倍率C随着ΔSOC增大而增大，随着δSOC增大而增大；当超级电容SOC处于上升段时，超级电容需求充电功率的计算不考虑δSOC，即充电倍率C只和ΔSOC有关，同样制定C关于ΔSOC的一维插值矩阵M2，该矩阵满足充电倍率C随着ΔSOC增大而增大；

②数据获取：

超级电容状态检测，包括超级电容SOC当前值SOC_real、超级电容SOC上一步长值SOC_pre、当前允许的最大充电功率P_{chg_max}；SOC_real、SOC_pre及P_{chg_max}可从超级电容控制器处直接获取；由此计算ΔSOC与δSOC，公式如下：

ΔSOC＝SOC_real-SOC_control (1)

δSOC＝(SOC_real-SOC_pre)/Δt (2)

式中Δt为计算步长；

③超级电容需求充电功率计算：

先利用ΔSOC和δSOC根据M1和M2插值得到充电倍率C，再结合基础充电功率P_{chg_base}得到需求充电功率P_{chg_req}，公式如下：

P_{chg_req}＝C*P_{chg_base} (3)

④对需求充电功率P_{chg_req}进行修正：

SOC_real＞SOC_control时，P_{chg_req}=0

SOC_real＜SOC_low时，P_{chg_req}＝P_{chg_max}

II.发动机通过曲轴发出的功率在传输过程中的损失细致划分为电路径上的损失功率和机械路径上的损失功率，其中电路径损失功率P_{loss_e}分为电机MG1损失功率P_{loss_MG1}和电机MG2损失功率P_{loss_MG2}，计算公式如下：

P_{loss_e}＝P_{loss_MG1}+P_{loss_MG2} (4)

P_{loss_MG1}及P_{loss_MG2}均根据电机MG1功率P_MG1及其能量转换效率η_MG1、电机MG2功率P_MG2及其能量转换效率η_MG2计算得到；P_MG1、η_MG1、P_MG2、η_MG2可分别从两电机各自控制器处直接获取；P_MG1、P_MG2均为电功率且规定电机电动状态下功率为正，发电状态下功率为负；两电机损失功率P_{loss_MG2}、P_{loss_MG1}的计算分为以下四种情况：

A.电机MG1发电、电机MG2电动：

电机MG1所发出的电功率全部来自发动机，其损失功率计算公式如下：

P_{loss_MG1}＝|P_MG1|/η_MG1-|P_MG1| (5)

电机MG2损失功率用两电机功率绝对值的较小值结合电机MG2的效率得到，计算公式如下：

P_{loss_MG2}＝min{|P_MG1|，|P_MG2|}*(1-η_MG2) (6)

B.电机MG1电动、电机MG2发电：

电机MG2所发出的电功率全部来自发动机，其损失功率计算公式如下：

P_{loss_MG2}＝|PM_G2|/η_MG2-|P_MG2| (7)

电机MG1损失功率用两电机功率绝对值的较小值结合电机MG1的效率得到，计算公式如下：

P_{loss_MG1}＝min{|P_MG1|，|P_MG2|}*(1-η_MG1) (8)

C.电机MG1电动、电机MG2电动：

此时发动机在电路径上没有功率传输，故电路径损失为零：

P_{loss_MG1}＝P_{oss_MG2}＝0 (9)

D.电机MG1发电、电机MG2发电：

两电机发出的电功率均全部来自发动机，其损失功率计算公式如下：

P_{loss_MG1}＝|P_MG1|/η_MG1-|P_MG1| (10)

P_{loss_MG2}＝|P_MG2|/η_MG2-|P_MG2| (11)

III.主减速器输入端需求整车驱动功率计算公式如下：

式中，6为旋转质量换算系数，单位1；m为整车质量，单位kg；a为车辆加速度，单位m/s²；f为滚动阻力系数，α为坡度，单位rad；C_D为车辆空气阻力系数，单位1；A为车辆迎风面积，单位m²；ρ为空气密度，单位kg/m³；v为车辆速度，单位m/s；η_final为主减速器传动效率，单位1；

IV.发动机需求功率计算公式如下：

P_req＝(P_{chg_req}+P_{loss_e}+P_drv)/η_pg (13)

式中η_pg为行星排传动效率，单位1。

本发明与现有技术相比，更加适合行星混联混合动力汽车多动力源的特点；提出基于超级电容SOC当前值低于控制目标程度及SOC变化速率双反馈的超级电容需求充电功率的计算，能够保证超级电容SOC在各工况下的平衡及快速恢复；充分考虑发动机通过曲轴发出功率在电路径和机械路径的损失，使得所求发动机需求功率接近系统真实需求，在提高控制系统响应速度的同时保证动力性充分发挥。

附图说明

本发明的上述和/或附加的方面和优点结合下面附图对实施例的描述中将变得明显和容易理解，其中：

图1为根据本发明实施例的超级电容SOC处于下降阶段时，充电倍率C的二维插值矩阵M1示意图；

图2为根据本发明实施例的超级电容SOC处于上升阶段时，充电倍率C的一维插值矩阵M2示意图；

图3为根据本发明实施例的超级电容SOC处于下降阶段时，充电倍率C关于超级电容SOC低于控制目标程度ΔSOC及超级电容SOC变化速率δSOC的变化曲面示意图；

图4为根据本发明实施例的超级电容SOC处于上升阶段时，充电倍率C关于超级电容SOC低于控制目标程度ΔSOC的变化曲线示意图；

图5为根据本发明实施例的电机MG1发电、电机MG2电动时能量流动示意图；

图6为根据本发明实施例的电机MG1电动、电机MG2发电时能量流动示意图；

图7为根据本发明实施例的电机MG1电动、电机MG2电动时能量流动示意图；

图8为根据本发明实施例的电机MG1发电、电机MG2发电时能量流动示意图；

具体实施方式

下面通过参考附图描述的实施例是示例性的，仅用于解释本发明，而不能理解为对本发明的限制。在本发明的描述中，需要理解的是，术语“行”、“列”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的数据必须具有特定的方位，因此不能理解为对本发明的限制。

由于行星混联混合动力汽车和行星齿轮的基本结构为本领域技术人员所悉知的，因此在此不再赘述。

本发明所述的一种行星混联混合动力汽车发动机需求功率计算方法，基于单行星排构型，发动机与行星架连接，电机MG1与太阳轮连接，电机MG2与齿圈连接，同时齿圈与主减速器输入端连接；将发动机需求功率分为主减速器输入端整车驱动需求功率、超级电容需求充电功率及动力传输过程中的损失功率，包括以下内容：

I.超级电容需求充电功率基于SOC低于控制目标程度ΔSOC及SOC变化速率δSOC双反馈计算，总体思路为先标定基础充电功率P_{chg_base}，然后在基础充电功率P_{chg_base}的基础上乘以充电倍率C得到P_{chg_req}，因此重点在于充电倍率C的计算。根据超级电容SOC平衡要求，应该遵循超级电容SOC低于SOC_control程度ΔSOC越大，充电倍率C越大；然而，在ΔSOC较小但超级电容SOC下降速率很快时，若只根据ΔSOC计算C，SOC会迅速降低，致使恢复缓慢，因此增加超级电容SOC下降段变化速率δSOC这一反馈。具体包括以下步骤：

①参数标定：

基于超级电容SOC平衡原则和寿命要求设定基础充电功率P_{chg_base}、超级电容SOC目标控制值SOC_control及SOC维持下限SOC_low；P_{chg_base}的大小应适应超级电容充电平稳、恢复迅速的要求；SOC_control的大小应能使超级电容始终保持较大的输出能力，且能够给再生制动充电预留足够的空间，即SOC_control应在50％到80％之间；SOC_low为超级电容SOC的下限，标定值应起到保护超级电容不过度放电、不损害寿命的作用。

超级电容SOC当前值SOC_real低于超级电容SOC目标控制值SOC_control时即开始充电，根据不同的SOC_real低于SOC_control的程度ΔSOC和不同的超级电容SOC变化速率δSOC，使用不同的充电倍率C。

当超级电容SOC处于下降段时，制定充电倍率C关于ΔSOC和δSOC的二维插值矩阵M1，该矩阵满足充电倍率C随着ΔSOC增大而增大，随着δSOC增大而增大；当超级电容SOC处于上升段时，超级电容需求充电功率的计算不考虑δSOC，即充电倍率C只和ΔSOC有关，同样制定C关于ΔSOC的一维插值矩阵M2，该矩阵满足充电倍率C随着ΔSOC增大而增大。

参阅附图1-2，二维插值矩阵M1的不同行对应ΔSOC，不同列对应δSOC，输出数据为C，行列互换成转置矩阵亦可。二维插值矩阵中ΔSOC的元素覆盖每个计算步长中SOC可能变化的大小；δSOC的元素覆盖每个计算步长中SOC可能变化的速率大小；一维插值矩阵M2中的每一个ΔSOC对应一个充电倍率C，ΔSOC的元素同样覆盖每个计算步长中SOC可能变化的大小。

矩阵输出数据C的取值范围下限为0，上限在标定过程中根据电池控制器反馈的超级电容SOC变化调整，应保证合理的充电功率，使超级电容SOC能较快恢复且不出现过充；根据不同的ΔSOC和δSOC，可根据该矩阵插值得到对应的充电倍率C。

参阅附图3，可看出当超级电容SOC处于下降段时充电倍率C关于ΔSOC和δSOC变化的趋势示意，C应随ΔSOC的增大或δSOC的增大呈线性或非线性增大趋势，具体变化关系根据实际情况人为标定调整。

参阅附图4，可看出当超级电容SOC处于上升段时充电倍率C关于ΔSOC变化的趋势示意，C应随ΔSOC的增大呈线性或非线性增大趋势，具体变化关系根据实际情况人为标定调整。

②数据获取：

超级电容状态检测，包括超级电容SOC当前值SOC_real、超级电容SOC上一步长值SOC_pre、当前允许的最大充电功率P_{chg_max}；SOC_real、SOC_pre及P_{chg_max}可从超级电容控制器处直接获取；由此计算ΔSOC与δSOC，公式如下：

ΔSOC＝SOC_real-SOC_control (1)

δSOC＝(SOC_real-SOC_pre)/Δt (2)

式中Δt为计算步长。

③超级电容需求充电功率计算：

先利用ΔSOC和δSOC插值得到充电倍率C，再结合基础充电功率P_{chg_base}得到需求充电功率P_{chg_req}，公式如下：

P_{chg_req}＝C*P_{chg_base} (3)

④对需求充电功率P_{chg_req}进行修正：

SOC_real＞SOC_control时，P_{chg_req}=0

SOC_real＜SOC_low时，P_{chg_req}＝P_{chg_max}

II.发动机通过曲轴发出的功率在传输过程中的损失细致划分为电路径上的损失功率和机械路径上的损失功率，其中电路径损失功率P_{loss_e}分为电机MG1损失功率P_{loss_MG1}和电机MG2损失功率P_{loss_MG2}，计算公式如下：

P_{loss_e}＝P_{loss_MG1}+P_{loss_MG2} (4)

P_{loss_MG1}根据电机MG1功率P_MG1及其能量转换效率η_MG1计算得到，P_{loss_MG2}根据电机MG2功率P_MG2及其能量转换效率η_MG2计算得到，P_MG1、η_MG1、P_MG2、η_MG2可从两电机各自控制器处直接获取；规定电机电动状态下功率为正，发电状态下功率为负，计算所得损失功率始终为正，功率单位均为W。

发动机需求功率为当前计算步长中系统状态对下一计算步长的需求，电路径损失功率应为下一步长中的实际值，但由于下一计算步长两电机功率及其能量转换效率无法预测，而一个计算步长内电路径损失功率变化微小，可认为不变，所以用当前计算步长的电路径损失功率代替下一步长的电路径损失功率。

参阅附图5-8，实线为机械功率传递路线，虚线为电功率传递路线。两电机损失功率P_{loss_MG2}、P_{loss_MG1}的计算分为下面四种情况。需要说明的是，这里的两电机损失功率均是在电机传递由发动机发出的功率时产生的，不讨论电机再生制动时情况。

A.电机MG1发电、电机MG2电动：

参阅附图5，电机MG1所发出的电功率全部来自发动机，其损失功率计算公式如下：

P_{loss_MG1}＝|P_MG1|/η_MG1-|P_MG1| (5)

如果|P_MG1|＞|P_MG2|，那么只有大小为|P_MG2|的来自发动机的功率通过路径A传递给电机MG2，电机MG1发出的功率剩余部分通过路径B给电池充电，路径c上功率为零；如果|P_MG1|＜|P_MG2|，那么只有大小为|P_MG1|的来自发动机的功率通过路径A传递给电机MG2，不足部分由超级电容通过路径c补充，路径B上的功率为零。综上，MG2损失功率计算公式如下：

P_{loss_MG2}＝min{|P_MG1|，|P_MG2|}*(1-η_MG2) (6)

B.电机MG1电动、电机MG2发电：

参阅附6，两电机损失功率的分析与A中类似，不再赘述，可得如下结论：

电机MG2所发出的电功率全部来自发动机，其损失功率计算公式如下：

P_{loss_MG2}＝|P_MG2|/η_MG2-|P_MG2| (7)

用两电机功率绝对值的较小值结合电机MG1的效率得到，计算公式如下：

P_{loss_MG1}＝min{|P_MG1|，|P_MG2|}*(1-η_MG1) (8)

c.电机MG1电动、电机MG2电动：

参阅附图7，此时发动机在电路径上没有功率传输，故电路径损失为零；

P_{loss_MG1}＝P_{loss_MG2}＝0 (9)

D.电机MG1发电、电机MG2发电：

参阅附图8，电机MG1、电机MG2发出的电功率均全部来自发动机，其损失功率计算公式如下：

P_{loss_MG1}＝|P_MG1|/η_MG1-|P_MG1| (10)

P_{loss_MG2}＝|P_MG2|/η_MG2-|P_MG2| (11)

III.主减速器输入端需求整车驱动功率计算公式如下：

式中，δ为旋转质量换算系数，单位1；m为整车质量，单位kg；a为车辆加速度，单位m/s²；f为滚动阻力系数，α为坡度，单位rad；C_D为车辆空气阻力系数，单位1；A车辆迎风面积，单位m²；ρ为空气密度，单位kg/m³；v为车辆速度，单位m/s；η_final为主减速器效率，单位1；

IV.发动机需求功率计算公式如下：

P_req＝(P_{chg_req}+P_{loss_MG1}+P_{loss_MG2}+P_drv)/η_pg (13)

式中η_pg为行星排传动效率，单位1。

机械路径的损失功率只存在于主减速器处和行星排处，由式(12)计算得到的主减速器输入端需求整车驱动功率P_drv中已包含主减速器处机械损失功率；由式(13)计算得到的发动机需求功率P_req中已包含行星排处机械损失功率。

Claims (1)

1.一种行星混联混合动力汽车发动机需求功率计算方法，将发动机需求功率分为超级电容需求充电功率、动力传输过程中的损失功率以及主减速器输入端整车驱动需求功率，其特征在于：包括以下内容：

I.超级电容需求充电功率基于超级电容SOC当前值低于超级电容SOC目标控制值程度及超级电容SOC变化速率双反馈计算，具体包括以下步骤：

①参数标定：

基于超级电容SOC平衡原则和寿命要求设定基础充电功率P_{chg_base}、超级电容SOC目标控制值SOC_control及SOC维持下限SOC_low；

超级电容SOC当前值SOC_real低于超级电容SOC目标控制值SOC_control时即开始充电，根据不同的SOC_real低于SOC_control的程度ΔSOC和不同的超级电容SOC变化速率δSOC，使用不同的充电倍率C；

当超级电容SOC处于下降段时，制定充电倍率C关于ΔSOC和δSOC的二维插值矩阵M1，该矩阵满足充电倍率C随着ΔSOC增大而增大，随着δSOC增大而增大；当超级电容SOC处于上升段时，超级电容需求充电功率的计算不考虑δSOC，即充电倍率C只和ΔSOC有关，同样制定C关于ΔSOC的一维插值矩阵M2，该矩阵满足充电倍率C随着ΔSOC增大而增大；

②数据获取：

超级电容状态检测，包括超级电容SOC当前值SOC_real、超级电容SOC上一步长值SOC_pre、当前允许的最大充电功率P_{chg_max}；SOC_real、SOC_pre及P_{chg_max}可从超级电容控制器处直接获取；由此计算ΔSOC与δSOC，公式如下：

ΔSOC＝SOC_real-SOC_control (1)

δSOC＝(SOC_real-SOC_pre)/Δt (2)

式中Δt为计算步长；

③超级电容需求充电功率计算：

先利用ΔSOC和δSOC根据M1和M2插值得到充电倍率C，再结合基础充电功率P_{chg_base}得到需求充电功率P_{chg_req}，公式如下：

P_{chg_req}＝C*P_{chg_base} (3)

④对需求充电功率P_{chg_req}进行修正:

SOC_real>SOC_control时，P_{chg_req}＝0

SOC_real<SOC_low时，P_{chg_req}＝P_{chg_max}

II.发动机通过曲轴发出的功率在传输过程中的损失细致划分为电路径上的损失功率和机械路径上的损失功率，其中电路径损失功率P_{loss_e}分为电机MG1损失功率P_{loss_MG1}和电机MG2损失功率P_{loss_MG2}，计算公式如下：

P_{loss_e}＝P_{loss_MG1}+P_{loss_MG2} (4)

P_{loss_MG1}及P_{loss_MG2}均根据电机MG1功率P_MG1及其能量转换效率η_MG1、电机MG2功率P_MG2及其能量转换效率η_MG2计算得到；P_MG1、η_MG1、P_MG2、η_2G2可分别从两电机各自控制器处直接获取；P_MG1、P_MG2均为电功率且规定电机电动状态下功率为正，发电状态下功率为负；两电机损失功率P_{loss_MG2}、P_{loss_MG1}的计算分为以下四种情况：

A.电机MG1发电、电机MG2电动：

电机MG1所发出的电功率全部来自发动机，其损失功率计算公式如下：

P_{loss_MG1}＝|P_MG1|/η_MG1-|P_MG1| (5)

电机MG2损失功率用两电机功率绝对值的较小值结合电机MG2的效率得到，计算公式如下：

P_{loss_MG2}＝min{|P_MG1|，|P_MG2|}*(1-η_MG2) (6)

B.电机MG1电动、电机MG2发电：

电机MG2所发出的电功率全部来自发动机，其损失功率计算公式如下：

P_{loss_MG2}＝|P_MG2|/η_MG2-|P_MG2| (7)

电机MG1损失功率用两电机功率绝对值的较小值结合电机MG1的效率得到，计算公式如下：

P_{loss_MG1}＝min{|P_MG1|，|P_MG2|}*(1-η_MG1) (8)

C.电机MG1电动、电机MG2电动：

此时发动机在电路径上没有功率传输，故电路径损失为零：

P_{loss_MG1}＝P_{loss_MG2}＝0 (9)

D.电机MG1发电、电机MG2发电：

两电机发出的电功率均全部来自发动机，其损失功率计算公式如下：

P_{loss_MG1}＝|P_MG1|/η_MG1-|P_MG1| (10)

P_{loss_MG2}＝|P_MG2|/η_MG2-|P_MG2| (11)

III.主减速器输入端需求整车驱动功率计算公式如下：

式中，δ为旋转质量换算系数，单位1；m为整车质量，单位kg；a为车辆加速度，单位m/s²；f为滚动阻力系数，α为坡度，单位rad；C_D为车辆空气阻力系数，单位1；A为车辆迎风面积，单位m²；ρ为空气密度，单位kg/m³；v为车辆速度，单位m/s；η_final为主减速器传动效率，单位1；

IV.发动机需求功率计算公式如下：

P_req＝(P_{chg_req}+P_{loss_e}+P_drv)/η_pg (13)

式中η_pg为行星排传动效率，单位1。