CN105128855B - 一种双轴并联混合动力城市公交车的控制方法 - Google Patents
一种双轴并联混合动力城市公交车的控制方法 Download PDFInfo
- Publication number
- CN105128855B CN105128855B CN201510604187.0A CN201510604187A CN105128855B CN 105128855 B CN105128855 B CN 105128855B CN 201510604187 A CN201510604187 A CN 201510604187A CN 105128855 B CN105128855 B CN 105128855B
- Authority
- CN
- China
- Prior art keywords
- torque
- motor
- engine
- tin
- battery
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 36
- 238000009826 distribution Methods 0.000 claims abstract description 12
- 239000000446 fuel Substances 0.000 claims description 33
- ATJFFYVFTNAWJD-UHFFFAOYSA-N Tin Chemical compound [Sn] ATJFFYVFTNAWJD-UHFFFAOYSA-N 0.000 claims description 25
- 230000008859 change Effects 0.000 claims description 23
- 230000005540 biological transmission Effects 0.000 claims description 15
- 239000011159 matrix material Substances 0.000 claims description 12
- 238000003860 storage Methods 0.000 claims description 12
- 238000010586 diagram Methods 0.000 claims description 10
- 230000008878 coupling Effects 0.000 claims description 8
- 238000010168 coupling process Methods 0.000 claims description 8
- 238000005859 coupling reaction Methods 0.000 claims description 8
- 230000005611 electricity Effects 0.000 claims description 8
- 239000000295 fuel oil Substances 0.000 claims description 6
- 238000012886 linear function Methods 0.000 claims description 6
- 230000001133 acceleration Effects 0.000 claims description 5
- 239000003921 oil Substances 0.000 claims description 4
- 230000001172 regenerating effect Effects 0.000 claims description 4
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 238000013459 approach Methods 0.000 claims description 3
- 238000011156 evaluation Methods 0.000 claims description 3
- 230000033001 locomotion Effects 0.000 claims description 3
- 238000005096 rolling process Methods 0.000 claims description 3
- 230000007704 transition Effects 0.000 claims description 3
- 238000005457 optimization Methods 0.000 abstract description 12
- 230000007547 defect Effects 0.000 abstract description 2
- 238000002474 experimental method Methods 0.000 abstract description 2
- 230000008901 benefit Effects 0.000 description 4
- 238000011217 control strategy Methods 0.000 description 4
- 230000008450 motivation Effects 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 3
- 210000004556 brain Anatomy 0.000 description 1
- 238000002485 combustion reaction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 239000002283 diesel fuel Substances 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000003344 environmental pollutant Substances 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 231100000719 pollutant Toxicity 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60W—CONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
- B60W20/00—Control systems specially adapted for hybrid vehicles
- B60W20/10—Controlling the power contribution of each of the prime movers to meet required power demand
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60W—CONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
- B60W50/00—Details of control systems for road vehicle drive control not related to the control of a particular sub-unit, e.g. process diagnostic or vehicle driver interfaces
- B60W50/0097—Predicting future conditions
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60W—CONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
- B60W50/00—Details of control systems for road vehicle drive control not related to the control of a particular sub-unit, e.g. process diagnostic or vehicle driver interfaces
- B60W2050/0001—Details of the control system
- B60W2050/0002—Automatic control, details of type of controller or control system architecture
- B60W2050/0013—Optimal controllers
Landscapes
- Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Transportation (AREA)
- Mechanical Engineering (AREA)
- Human Computer Interaction (AREA)
- Electric Propulsion And Braking For Vehicles (AREA)
- Hybrid Electric Vehicles (AREA)
Abstract
本发明公开了一种双轴并联混合动力城市公交车的控制方法,包括以下步骤:建立动力传动***模型;进行线性拟合;预确定荷电状态SOC可达区域;建立预测控制模型;求解混合整数线性规划问题。本发明通过工作模式的划分和混合逻辑动态的特点,使混合动力公交车可以根据行驶工况和驾驶员行为进行工作模式的自动切换,从而用最优化的方法进行发动机和电动机的转矩分配。有效地解决了基于规则的控制方法所带来的无法根据行驶工况自动优化转矩分配的问题。本发明通过预测模型,有效地解决了瞬时优化中只能考虑当前时刻的缺陷;并且预测视距可根据实验找到一个使预测结果和理想值相对误差较小的时长,有效地解决了全局优化中对整个行驶工况的苛刻要求。
Description
技术领域
本发明涉及一种双轴并联混合动力公交车的能量控制策略,特别涉及到一种双轴并联混合动力城市公交车的控制方法。
背景技术
当前,发展以节能、环保为目标的新能源汽车成为世界各国的共识。混合动力汽车是第一个实现量产的新能源汽车。现有的混合动力汽车一般是以传统的内燃机为主要动力源,以电池电动机***为辅助动力源,通过优化控制***使发动机工作在高效率区,并且实现再生制动能量回收,与传统液体燃料汽车相比,大大提高了燃油经济性,降低了污染物排放。因此,在目前的条件下,混合动力汽车成为世界各国新能源汽车发展的首要选择。控制策略如同混合动力汽车的大脑,指挥发动机、电动机两个动力源及其它部件的协调工作,是结合传统燃油汽车和纯电动汽车优势的纽带,决定着混合动力汽车综合性能优势的充分发挥,是决定其整车燃油经济性、排放和效率的关键,对于提高整车性能、降低成本具有重要意义。按何种规律和方法分配发动机和电动机功率才能使混合动力***综合性能实现最佳是混合动力汽车控制策略研究的主要问题。借助车载导航***提供的汽车未来运行状态信息的混合动力汽车预测控制策略研究在提高汽车的燃油经济性方面具有重要的意义。
目前投入使用的双轴并联混合动力城市公交车通过发动机和电动机的共同驱动,混合动力车通过控制***,可以选择不同的工作模式。目前混合动力汽车控制方法都是对电动机和发动机转矩分配的控制,大致可以分为两类:基于规则的和基于优化方法的。总的来说,基于规则的控制方法由于功率的分配没有经过优化,不能充分发挥混合动力公交车的性能优点。其在规则中所设计的各种控制参数,比如启动发动机的速度门限值,严重依赖特定的混合动力***以及行驶工况。而且,在一种工况下表现良好的参数并不能很好的工作于另一种工况之下。因此,目前混合动力公交车控制方法的研究更多的集中在基于优化的控制方法上。基于优化的控制方法主要分为全局优化和瞬时优化两大类:全局优化控制方法因为需要准确知道整个行驶工况,要求过于苛刻;基于瞬时等效燃油最小的能量控制方法又只是对当前时刻的优化,对车辆未来行驶工况信息考虑不够,并且该种控制方法比较容易造成蓄电池SOC的过放电现象,造成蓄电池SOC处于一个较低值,反而无法发挥混合动力公交车燃油经济性的最大优势。总的来说,现有双轴并联混合动力城市公交车存在的问题主要是对行驶工况的预知要求过于苛刻,并且蓄电池的过放电现象导致无法发挥混合动力公交车燃油经济性的最大优势。
发明内容
为了克服现有混合动力公交车上所存在的一些问题,本发明要提出一种既无需准确知道整个行驶工况、又可以使电动机和发动机高效配合从而提升混合动力公交车燃油经济性的双轴并联混合动力公交车控制方法。
为实现上述目的,本发明的技术方案如下:一种双轴并联混合动力城市公交车的控制方法,所述的双轴并联混合动力城市公交车的动力传动***包括发动机和电动机两套单独驱动***,具体部件包括发动机、电动机、电池、转矩耦合器、变速器、主减速器和车轮。所述的电动机由电池供电能驱动车辆,也具有发电机的功能并向电池充电。所述的发动机和电动机分别安装在两根轴上,并用专门的转矩耦合器将其转矩耦合,再传递给驱动轴,转矩耦合比为2。
所述的控制方法,包括以下步骤:
A、建立动力传动***模型:
建立双轴并联混合动力城市公交车动力结构数学模型,模型如下:
A1、计算车轮的转速和转矩:
A2、计算变速器的转速转矩:
A3、建立双轴并联转矩耦合器力学模型:
A4、计算发动机燃油消耗率:
A5、计算电动机等效燃油消耗率:
A6、计算电池输出电流:
绘出混合动力公交车电池电动机***的等效电路图,并按下式计算输出电流:
上述各式中:ωW为车轮转速,单位为r/min;v为车速,单位为m/s;TW为车轮半径,单位为m;TW为车轮转矩,单位为N·m;Af为迎风面积,单位为平方米;Cd为风阻系数,fr为轮胎滚动阻力系数;m为汽车质量,单位kg;J为汽车总的转动惯量,单位kg·m2;α为车辆行驶偏移角。ωin为变速箱输入端的需求转速,单位为r/min;i0为变速箱处于某档时车辆总的传动比;Tin为变速箱输入端的需求转矩,单位为N·m;η为变速箱处于某档时从变速箱输入端到车轮总的传动效率。为发动机转矩,单位为N·m;为电动机转矩,单位为N·m;为发动机转速,单位为r/min;为电动机转速,单位为r/min。为发动机的燃油消耗率,它是转矩和转速的函数。为电动机的用电量的消耗等效的燃油消耗率,和wch分别表示电动机作为电动机消耗电池能量和电动机作为发电机给电池充电时的等效燃油因子,R为柴油质量热值常数,取33000kJ/kg。I为混合动力公交车电池电动机***的等效电路图中的电池输出电流,单位为A;U0为电池的开路电压,单位为V;Ri为电池等效内阻,单位为欧姆;ɑ为汽车的加速度,单位为m/s2。
B、进行线性拟合:
用最小二乘法将发动机不同转速下的燃油消耗和发动机扭矩之间的关系进行线性拟合;同理将单位时间内电池的荷电变化和电动机转矩变化之间的关系进行线性拟合。
B1、建立发动机燃油消耗分段线性模型:
根据发动机燃油效率数据,将发动机转速的范围以100rpm为区间长度划分、发动机转矩范围以10Nm区间长度划分,得到不同转速下发动机燃油消耗率随转矩的变化曲线图,用最小二乘法近似拟合为一次函数,如式(7),其中p1和p0分别为一次项系数和常数项。并令P0=a0,P1=a1。
B2、建立电池荷电变化率分段线性模型:
将步骤A6中的电流输出转换成电池荷电变化率,再将电动机转速的范围以1000rpm为区间长度划分,用矩阵实验室软件MATLAB绘出一定转速下电池荷电变化率随转矩的变化曲线图,用最小二乘法近似拟合为一次函数,如式(8),其中q1和q0分别为一次项系数和常数项。当Tm>0时,q1=b1,q0=b0;Tm<0时,q1=c1,q0=c0。
C、预确定荷电状态SOC可达区域:
将蓄电池荷电状态SOC作为混合动力公交车的状态变量,需要计算出预测区域内荷电状态SOC的可达区域。当汽车行驶在k时刻时,当前时刻的蓄电池荷电状态SOC值及车速已知,用以下方法计算出电动机在当前时刻能提供的转矩最值:
C1、确定车辆处于驱动状态时电动机转矩最值:
C11、电动机最大正转矩取下面情况的最小值:
C111、电动机提供车轮需求的所有驱动功率时电动机的转矩;
C112、以纯电动驱动,或当前蓄电池荷电状态SOC下电动机本身能提供的最大正转矩;
C12、发电机的最小转矩取下面情况的最大值:
C121、满足汽车行驶驱动所需要发动机提供的转矩与在当前转速下发动机能提供的最大转矩的差值;
C122、在当前荷电状态下发电机本身可提供的最小转矩;
C2、确定车辆处于制动状态时电动机转矩最值:
C21、电动机的最大正转矩为0;发电机的最小转矩取下面情况的最大值:
C211、减速或下坡制动产生的全部再生制动能量通过发电机向蓄电池充电;
C212、前荷电状态SOC下发电机本身可提供的最小转矩。
进而用式(6)和式(9)计算出在k+1时刻的荷电状态SOC最小、最大值。以此类推,反复计算得到在预测区域p内所有时刻的荷电状态SOC的最大、最小值,这样就得到预测区域p内***状态变量——蓄电池荷电状态SOC的可达区域。
D、建立预测控制模型:
基于混合逻辑动态标准模型,结合模型预测控制思想建立混合动力城市公交车的转矩控制模型,包括状态转移方程和评价指标方程,如式(10)所示;六种工作模式的约束方程如下:
式中:
B2=[b0(k )0b0(k) c0(k) c0(k)0];
B3=[b1(k) 0b1(k) c1(k) c1(k)0];
D2=[0a0(k)+a1(k)*Tin(k)a0(k)+a1(k)*Tin(k)a0(k)+a1(k)*Tin(k)00];
其中,x(k)为k时刻电池的荷电状态SOC的值,y(k)为k时刻等效燃油消耗率,δ(k)为6*1的矩阵逻辑变量,为k时刻的工作模式矩阵,每个分量只能取0或1且分量和为1,z(k)为辅助变量,z(k)=δ(k)·u(k)。
建立城市公交车工况的六种工作模式的约束方程,包括:
纯电动模式:2Minδ1≤2Min+Tin-ε;
2Minδ1≤2Min-Tin+2Tm_max;
-Tinδ1+2z1≤0;
Tinδ1-2z1≤0;
油动模式:Minδ2≤Min+Tin-ε;
Minδ2≤Min-Tin+Te_max;
z2≤0;
-z2≤0;
油电混动模式:Minδ3≤-2u+Min+Tin-ε;
Minδ3≤2u+Min-Tin+Te_max;
Mmδ3≤u+Mm-ε;
Mmδ3≤-u+Mm-Tm_max;
行车充电模式:Minδ4≤-2u+Min+Tin-ε;
Minδ4≤2u+Min-Tin+Te_max;
Mmδ4≤u+Mm+Tm_max;
Mmδ4≤-u+Mm-ε;
制动回收模式:minδ5≤min-Tin-ε;
停车模式:Tinδ6≤0;
-Tinδ6≤0;
z6≤0;
-z6≤0。
其中,Mm为电动机转矩最大值,Min为最大需求转矩,Te_max为当前时刻发动机最大转矩,Tm_max为当前时刻电动机最大转矩,ε为机器精度,选取ε=0.0001;Tin为当前时刻的需求转矩,δi为选取第i个工作模式,即向量δt的第i个分量为1,其余分量为0;zi为选取第i个工作模式。
将以上不等式约束引入时间序列k,统一为矩阵不等式:E2δ(k)+E3z(k)≤Elu(k)+E5,其中,
E、求解混合整数线性规划问题:
为减少控制算法的搜索空间,引入启发式算法如式(11)-(13),将步骤D中建立的模型转化为混合整数线性规划问题。
当Tin=0时,δ6=1,u=0; (11)
当Tin<0时,δ5=1,u=Tin; (12)
当Tin>Te_max时,δ3=1,u=Tin-Te max。 (13)
当车辆停车时,发动机和电动机转矩均为0,采用停车模式;
当车辆制动时,需求转矩为负数,直接采用制动回收模式;
当车辆处于加速的高动力行驶工况时,直接采用油电混动模式,同时对电动机和发动机的扭矩分配也有限制:发动机提供其最大转矩,剩余需求转矩由电动机提供。
在城市交通中,车辆处在不断的加速、制动运行工况时,即高功率需求和制动的运行工况下。那么,主要集中在车辆中、小功率需求时,车辆在纯电动模式、油动模式、行车充电模式这三种工作模式间的选择和转矩在发动机和电动机间的分配。
求解该混合整数线性规划问题,得出随控制时间变化的工作模式序列和电动机转矩控制序列。根据式(3)发动机与电动机的转矩关系,最终确定发动机与电动机的转矩分配从而对双轴并联混合动力公交车进行控制。
与现有技术相比,本发明具有以下有益效果:
1、本发明通过工作模式的划分和混合逻辑动态的特点,使混合动力公交车可以根据行驶工况和驾驶员行为进行工作模式的自动切换,从而用最优化的方法进行发动机和电动机的转矩分配。有效地解决了基于规则的控制方法所带来的无法根据行驶工况自动优化转矩分配的问题。
2、本发明通过预测模型,基于当前的工作状况和未来一段时间的需求转矩,对下一时刻的转矩和工作模式进行预测的思想,使得车辆全面参考了未来一定时间段的工况信息量,有效地解决了瞬时优化中只能考虑当前时刻的缺陷;并且预测视距可根据实验找到一个使预测结果和理想值相对误差较小的时长,有效地解决了全局优化中对整个行驶工况的苛刻要求。
3、本发明通过预测蓄电池SOC的可达区域,使得蓄电池SOC在整个工况中始终保持相对平稳,有效地解决了蓄电池SOC的过放电现象。
附图说明
本发明共有附图3张,其中:
图1是双轴并联混合动力公交车整车模型模块示意图。
图2是本发明双轴并联混合动力城市公交车动力传动***示意图。
图3是混合动力城市公交车电池电动机***的等效电路图。
具体实施方式
下面结合附图对本发明进行进一步地描述。
图1所示为双轴并联混合动力公交车整车模型模块示意图,其动力传动***如图2所示,所述的车辆为双轴并联混合动力城市公交车,由发动机和电动机两套单独驱动***共同驱动,电动机为交流异步电动机,由于是双轴并联公交车,可以最大程度地发挥交流异步电动机结构简单,有较高的转速极限的特点。传动***中只有一台电动机,可以当作电动机由电池供电驱动车辆,也可以当作发电机向电池充电。发动机为柴油发动机,功率大,最大输出转矩大,可以为城市公交车提供足够的转矩。其动力传动***还包括电池、转矩耦合器、变速器、主减速器和车轮等。发动机通过离合器与转矩耦合器连接,离合器控制发动机的转矩是否传递到转矩耦合器上,电动机与电池相连,并将转矩通过齿轮啮合传递到转矩耦合器上与发动机转矩耦合。车辆的发动机和电动机分别在两根轴上,并用专门的转矩耦合器将其转矩耦合,转矩耦合比为2,再通过5档变速箱传递给驱动轴,驱动车轮转动。
以下详细叙述该发明基于预测控制的双轴并联混合动力公交车的控制方法。具体包括以下步骤:
根据对公交车行驶工况分析,综合运用汽车理论知识和汽车仿真软件ADVISOR,对双轴并联混合动力公交车动力总成进行合理的参数匹配,选取数据配置合适的发动机、电动机、电池、转矩耦合器等部件使之满足公交车的动力及经济性能指标,为进一步研究能量管理策略打好基础。
A、建立动力传动***模型:
建立双轴并联混合动力城市公交车动力结构数学模型,模型如下:
A1、计算车轮的转速和转矩:
A2、计算变速器的转速转矩:
A3、建立双轴并联转矩耦合器力学模型:
A4、计算发动机燃油消耗率:
A5、计算电动机等效燃油消耗率:
A6、计算电池输出电流:
绘出图3所示的混合动力公交车电池电动机***的等效电路图,并按下式计算输出电流:
上述各式中:ωW为车轮转速,单位为r/min;v为车速,单位为m/s;TW为车轮半径,单位为m;TW为车轮转矩,单位为N·m;Af为迎风面积,单位为平方米;Cd为风阻系数,fr为轮胎滚动阻力系数;m为汽车质量,单位kg;J为汽车总的转动惯量,单位kg-m2;α为车辆行驶偏移角。ωin为变速箱输入端的需求转速,单位为r/min;i0为变速箱处于某档时车辆总的传动比;Tin为变速箱输入端的需求转矩,单位为N·m;η为变速箱处于某档时从变速箱输入端到车轮总的传动效率。T0为发动机转矩,单位为N·m;为电动机转矩,单位为N·m;ω0为发动机转速,单位为r/min;ωm为电动机转速,单位为r/min。为发动机的燃油消耗率,它是转矩和转速的函数。为电动机的用电量的消耗等效的燃油消耗率,wdis和wch分别表示电动机作为电动机消耗电池能量和电动机作为发电机给电池充电时的等效燃油因子,R为柴油质量热值常数,取33000kJ/kg。I为混合动力公交车电池电动机***的等效电路图中的电池输出电流(如附图3),单位为A;U0为电池的开路电压,单位为V;Ri为电池等效内阻,单位为欧姆。
B、进行线性拟合:
用最小二乘法将发动机不同转速下的燃油消耗和发动机扭矩之间的关系进行线性拟合;同理将单位时间内电池的荷电变化和电动机转矩变化之间的关系进行线性拟合。
B1、建立发动机燃油消耗分段线性模型:
根据发动机燃油效率数据,将发动机转速的范围以100rpm为区间长度划分、发动机转矩范围以10Nm区间长度划分,得到不同转速下发动机燃油消耗率随转矩的变化曲线图,用最小二乘法近似拟合为一次函数,如式(7),其中p1和p0分别为一次项系数和常数项。并令p0=a0,p1=a1。
B2、建立电池荷电变化率分段线性模型:
将步骤A6中的电流输出转换成电池荷电变化率,再将电动机转速的范围以1000rpm为区间长度划分,用矩阵实验室软件MATLAB绘出一定转速下电池荷电变化率随转矩的变化曲线图,用最小二乘法近似拟合为一次函数,如式(8),其中q1和q0分别为一次项系数和常数项。当Tm>0时,q1=b1,q0=b0;Tm<0时,q1=c1,q0=c0。
C、预确定荷电状态SOC可达区域:
将蓄电池荷电状态SOC作为混合动力公交车的状态变量,需要计算出预测区域内荷电状态SOC的可达区域。当汽车行驶在k时刻时,当前时刻的蓄电池荷电状态SOC值及车速已知,用以下方法计算出电动机在当前时刻能提供的转矩最值:
C1、确定车辆处于驱动状态时电动机转矩最值:
C11、电动机最大正转矩取下面情况的最小值:
C111、电动机提供车轮需求的所有驱动功率时电动机的转矩;
C112、以纯电动驱动,或当前蓄电池荷电状态SOC下电动机本身能提供的最大正转矩;
C12、发电机的最小转矩取下面情况的最大值:
C121、满足汽车行驶驱动所需要发动机提供的转矩与在当前转速下发动机能提供的最大转矩的差值;
C122、在当前荷电状态下发电机本身可提供的最小转矩;
C2、确定车辆处于制动状态时电动机转矩最值:
C21、电动机的最大正转矩为0;发电机的最小转矩取下面情况的最大值:
C211、减速或下坡制动产生的全部再生制动能量通过发电机向蓄电池充电;
C212、前荷电状态SOC下发电机本身可提供的最小转矩。
进而用式(6)和式(9)计算出在k+1时刻的荷电状态SOC最小、最大值。以此类推,反复计算得到在预测区域p内所有时刻的荷电状态SOC的最大、最小值,这样就得到预测区域p内***状态变量——蓄电池荷电状态SOC的可达区域。
D、建立预测控制模型:
基于混合逻辑动态标准模型,结合模型预测控制思想建立混合动力城市公交车的转矩控制模型,包括状态转移方程和评价指标方程,如式(10)所示;六种工作模式的约束方程如下:
式中:
B2=[b0(k) 0b0(k) c0(k) c0(k)0];
B3=[b1(k) 0b1(k )c1(k) c1(k)0];
D2=[0a0(k)+a1(k)*Tin(k)a0(k)+a1(k)*Tin(k)a0(k)+a1(k)*Tin(k)00];
其中,x(k)为k时刻电池的荷电状态SOC的值,y(k)为k时刻等效燃油消耗率,δ(k)为6*1的矩阵逻辑变量,为k时刻的工作模式矩阵,每个分量只能取0或1且分量和为1,z(k)为辅助变量,z(k)=δ(k)·u(k)。
建立城市公交车工况的六种工作模式的约束方程,包括:
纯电动模式:2Minδ1≤2Min+Tin-ε;
2Minδ1≤2Min-Tin+2Tm_max;
-Tinδ1+2z1≤0;
Tinδ1-2z1≤0;
油动模式:Minδ2≤Min+Tin-ε;
Minδ2≤Min-Tin+Te_max;
z2≤0;
-z2≤0;
油电混动模式:minδ3≤-2u+Min+Tin-ε;
Minδ3≤2u+Min-Tin+Te_max;
Mmδ3≤u+Mm-ε;
Mmδ3≤-u+Mm-Tm_max;
行车充电模式:Minδ4≤-2u+Min+Tin-ε;
Minδ4≤2u+Min-Tin+Te_max;
Mmδ4≤u+Mm+Tm_max;
Mmδ4≤-u+Mm-ε;
制动回收模式:Minδ5≤Min-Tin-ε;
停车模式:Tinδ6≤0;
-Tinδ6≤0;
z6≤0;
-z6≤0。
其中,Mm为电动机转矩最大值,Min为最大需求转矩,Te_max为当前时刻发动机最大转矩,Tm_max为当前时刻电动机最大转矩,ε为机器精度,选取ε=0.0001;Tin为当前时刻的需求转矩,δi为选取第i个工作模式,即向量δt的第i个分量为1,其余分量为0;zi为选取第i个工作模式。
将以上不等式约束引入时间序列k,统一为矩阵不等式:E2δ(k)+E3z(k)≤E1u(k)+E5,其中,
E、求解混合整数线性规划问题:
为减少控制算法的搜索空间,引入启发式算法如式(11)-(13),将步骤D中建立的模型转化为混合整数线性规划问题。
当Tin=0时,δ6=1,u=0; (11)
当Tin<0时,δ5=1,u=Tin; (12)
当Tin>Te_max时,δ3=1,u=Tin-Te_max。 (13)
当车辆停车时,发动机和电动机转矩均为0,采用停车模式;
当车辆制动时,需求转矩为负数,直接采用制动回收模式;
当车辆处于加速的高动力行驶工况时,直接采用油电混动模式,同时对电动机和发动机的扭矩分配也有限制:发动机提供其最大转矩,剩余需求转矩由电动机提供。
在城市交通中,车辆处在不断的加速、制动运行工况时,即高功率需求和制动的运行工况下。那么,主要集中在车辆中、小功率需求时,车辆在纯电动模式、油动模式、行车充电模式这三种工作模式间的选择和转矩在发动机和电动机间的分配。
求解该混合整数线性规划问题,得出随控制时间变化的工作模式序列和电动机转矩控制序列。根据式(3)发动机与电动机的转矩关系,最终确定发动机与电动机的转矩分配从而对双轴并联混合动力公交车进行控制。
Claims (1)
1.一种双轴并联混合动力城市公交车的控制方法,所述的双轴并联混合动力城市公交车的动力传动***包括发动机和电动机两套单独驱动***,具体部件包括发动机、电动机、电池、转矩耦合器、变速器、主减速器和车轮;所述的电动机由电池供电能驱动车辆,也具有发电机的功能并向电池充电;所述的发动机和电动机分别安装在两根轴上,并用专门的转矩耦合器将其转矩耦合,再传递给驱动轴,转矩耦合比为2;
其特征在于:所述的控制方法,包括以下步骤:
A、建立动力传动***模型:
建立双轴并联混合动力城市公交车动力结构数学模型,模型如下:
A1、计算车轮的转速和转矩:
A2、计算变速器的转速转矩:
A3、建立双轴并联转矩耦合器力学模型:
A4、计算发动机燃油消耗率:
A5、计算电动机等效燃油消耗率:
A6、计算电池输出电流:
绘出混合动力公交车电池电动机***的等效电路图,并按下式计算输出电流:
上述各式中:ωW为车轮转速,单位为r/min;v为车速,单位为m/s;rw为车轮半径,单位为m;Tw为车轮转矩,单位为N·m;Af为迎风面积,单位为平方米;Cd为风阻系数,fr为轮胎滚动阻力系数;m为汽车质量,单位kg;J为汽车总的转动惯量,单位kg·m2;α为车辆行驶偏移角;ωin为变速箱输入端的需求转速,单位为r/min;i0为变速箱处于某档时车辆总的传动比;Tin为变速箱输入端的需求转矩,单位为N·m;η为变速箱处于某档时从变速箱输入端到车轮总的传动效率;Te为发动机转矩,单位为N·m;Tm为电动机转矩,单位为N·m;ωe为发动机转速,单位为r/min;ωm为电动机转速,单位为r/min;为发动机的燃油消耗率,它是转矩和转速的函数;为电动机的用电量的消耗等效的燃油消耗率,wdis和wch分别表示电动机作为电动机消耗电池能量和电动机作为发电机给电池充电时的等效燃油因子,R为柴油质量热值常数,取33000kJ/kg;I为混合动力公交车电池电动机***的等效电路图中的电池输出电流,单位为A;U0为电池的开路电压,单位为V;Ri为电池等效内阻,单位为欧姆;ɑ为汽车的加速度,单位为m/s2;
B、进行线性拟合:
用最小二乘法将发动机不同转速下的燃油消耗和发动机扭矩之间的关系进行线性拟合;同理将单位时间内电池的荷电变化和电动机转矩变化之间的关系进行线性拟合;
B1、建立发动机燃油消耗分段线性模型:
根据发动机燃油效率数据,将发动机转速的范围以100rpm为区间长度划分、发动机转矩范围以10Nm区间长度划分,得到不同转速下发动机燃油消耗率随转矩的变化曲线图,用最小二乘法近似拟合为一次函数,如式(7),其中p1和P0分别为一次项系数和常数项;并令P0=a0,p1=a1;
B2、建立电池荷电变化率分段线性模型:
将步骤A6中的电流输出转换成电池荷电变化率,再将电动机转速的范围以1000rpm为区间长度划分,用矩阵实验室软件MATLAB绘出一定转速下电池荷电变化率随转矩的变化曲线图,用最小二乘法近似拟合为一次函数,如式(8),其中q1和q0分别为一次项系数和常数项;当Tm>0时,q1=b1,q0=b0;Tm<0时,q1=c1,q0=c0;
C、预确定荷电状态SOC可达区域:
将蓄电池荷电状态SOC作为混合动力公交车的状态变量,需要计算出预测区域内荷电状态SOC的可达区域;当汽车行驶在k时刻时,当前时刻的蓄电池荷电状态SOC值及车速已知,用以下方法计算出电动机在当前时刻能提供的转矩最值:
C1、确定车辆处于驱动状态时电动机转矩最值:
C11、电动机最大正转矩取下面情况的最小值:
C111、电动机提供车轮需求的所有驱动功率时电动机的转矩;
C112、以纯电动驱动,或当前蓄电池荷电状态SOC下电动机本身能提供的最大正转矩;
C12、发电机的最小转矩取下面情况的最大值:
C121、满足汽车行驶驱动所需要发动机提供的转矩与在当前转速下发动机能提供的最大转矩的差值;
C122、在当前荷电状态下发电机本身可提供的最小转矩;
C2、确定车辆处于制动状态时电动机转矩最值:
C21、电动机的最大正转矩为0;发电机的最小转矩取下面情况的最大值:
C211、减速或下坡制动产生的全部再生制动能量通过发电机向蓄电池充电;
C212、前荷电状态SOC下发电机本身可提供的最小转矩;
进而用式(6)和式(9)计算出在k+1时刻的荷电状态SOC最小、最大值;以此类推,反复计算得到在预测区域p内所有时刻的荷电状态SOC的最大、最小值,这样就得到预测区域p内***状态变量——蓄电池荷电状态SOC的可达 区域;
D、建立预测控制模型:
基于混合逻辑动态标准模型,结合模型预测控制思想建立混合动力城市公交车的转矩控制模型,包括状态转移方程和评价指标方程,如式(10)所示;六种工作模式的约束方程如下:
式中:
B2=[b0(k) 0 b0(k) c0(k) c0(k) 0];
B3=[b1(k) 0 b1(k) c1(k) c1(k) 0];
D2=[0 a0(k)+a1(k)*Tin(k) a0(k)+a1(k)*Tin(k) a0(k)+a1(k)*Tin(k) 0 0];
D3=[2ω′disωin(k)-2a1(k)2ω′disωin(k)-2a1(k)2w′chωin(k)-2a1(k)2w′chωin(k)0];
其中,x(k)为k时刻电池的荷电状态SOC的值,y(k)为k时刻等效燃油消耗率,δ(k)为6*1的矩阵逻辑变量,为k时刻的工作模式矩阵,每个分量只能取0或1且分量和为1,z(k)为辅助变量,z(k)=δ(k)·u(k);
建立城市公交车工况的六种工作模式的约束方程,包括:
纯电动模式:2Min δ1≤2Min+Tin-ε;
2Min δ1≤2Min-Tin+2Tm_max;
-Tin δ1+2z1≤0;
Tin δ1-2z1≤0;
油动模式:Min δ2≤Min+Tin-ε;
Min δ2≤Min-Tin+Te_max;
z2≤0;
-z2≤0;
油电混动模式:Min δ3≤-2u+Min+Tin-ε;
Min δ3≤2u+Min-Tin+Te_max;
Mm δ3≤u+Mm-ε;
Mm δ3≤-u+Mm-Tm_max;
行车充电模式:Min δ4≤-2u+Min+Tin-ε;
Min δ4≤2u+Min-Tin+Te_max;
Mm δ4≤u+Mm+Tm_max;
Mm δ4≤-u+Mm_ε;
制动回收模式:Min δ5≤Min_Tin_ε;
停车模式:Tin δ6≤0;
-Tin δ6≤0;
z6≤0;
-z6≤0;
其中,Mm为电动机转矩最大值,Min为最大需求转矩,Te_max为当前时刻发动机最大转矩,Tm_max为当前时刻电动机最大转矩,ε为机器精度,选取ε=0.0001;Tin为当前时刻的需求转矩,δi为选取第i个工作模式,即向量δi的第i个分量为1,其余分量为0;zi为选取第i个工作模式;
将以上不等式约束引入时间序列k,统一为矩阵不等式:E2δ(k)+E3z(k)≤E1u(k)+E5,其中,
E、求解混合整数线性规划问题:
为减少控制算法的搜索空间,引入启发式算法如式(11)-(13),将步骤D中建立的模型转化为混合整数线性规划问题;
当Tin=0时,δ6=1,u=0; (11)
当Tin<0时,δ5=1,u=Tin; (12)
当Tin>Te_max时,δ3=1,u=Tin-Te_max; (13)
当车辆停车时,发动机和电动机转矩均为0,采用停车模式;
当车辆制动时,需求转矩为负数,直接采用制动回收模式;
当车辆处于加速的高动力行驶工况时,直接采用油电混动模式,同时对电动机和发动机的扭矩分配也有限制:发动机提供其最大转矩,剩余需求转矩由 电动机提供;
在城市交通中,车辆处在不断的加速、制动运行工况时,即高功率需求和制动的运行工况下;那么,主要集中在车辆中、小功率需求时,车辆在纯电动模式、油动模式、行车充电模式这三种工作模式间的选择和转矩在发动机和电动机间的分配;
求解该混合整数线性规划问题,得出随控制时间变化的工作模式序列和电动机转矩控制序列;根据式(3)发动机与电动机的转矩关系,最终确定发动机与电动机的转矩分配从而对双轴并联混合动力公交车进行控制。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510604187.0A CN105128855B (zh) | 2015-09-21 | 2015-09-21 | 一种双轴并联混合动力城市公交车的控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510604187.0A CN105128855B (zh) | 2015-09-21 | 2015-09-21 | 一种双轴并联混合动力城市公交车的控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105128855A CN105128855A (zh) | 2015-12-09 |
CN105128855B true CN105128855B (zh) | 2017-07-21 |
Family
ID=54714554
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510604187.0A Expired - Fee Related CN105128855B (zh) | 2015-09-21 | 2015-09-21 | 一种双轴并联混合动力城市公交车的控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105128855B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105691182B (zh) * | 2016-01-21 | 2017-12-29 | 厦门市福工动力技术有限公司 | 基于amt的混合动力***及其控制方法 |
CN106527515B (zh) * | 2016-12-12 | 2019-12-24 | 中国燃气涡轮研究院 | 一种对轮盘转速进行精确控制的方法 |
CN106926841B (zh) * | 2017-03-10 | 2018-04-24 | 江苏大学 | 一种双行星排式混合动力汽车能量管理控制方法 |
CN107891864B (zh) * | 2017-10-12 | 2020-05-15 | 北京动力源新能源科技有限责任公司 | 并联混合动力***的等效油电折算系数获取方法及装置 |
CN109334654A (zh) * | 2018-09-21 | 2019-02-15 | 江苏大学 | 一种带变速箱档位控制的并联式混合动力汽车能量管理方法 |
CN110509914B (zh) * | 2019-09-16 | 2020-08-04 | 重庆邮电大学 | 一种并联式混合动力汽车的能耗优化方法 |
CN111259493B (zh) * | 2020-02-09 | 2022-08-02 | 吉林大学 | 适用于智能网联车排放控制的车辆排放模型建模方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPH0998515A (ja) * | 1995-07-25 | 1997-04-08 | Nippon Soken Inc | ハイブリッド車のエンジン制御装置 |
KR100494919B1 (ko) * | 2003-03-05 | 2005-06-13 | 현대자동차주식회사 | 병렬형 하이브리드 전기자동차의 모터 제어방법 및 그제어장치 |
US8612078B2 (en) * | 2011-08-08 | 2013-12-17 | Bae Systems Controls Inc. | Parallel hybrid electric vehicle power management system and adaptive power management method and program therefor |
CN102431550B (zh) * | 2011-10-17 | 2013-12-25 | 大连理工大学 | 一种混合动力汽车控制方法 |
CN102416950B (zh) * | 2011-10-31 | 2013-06-19 | 大连理工大学 | 一种基于最小等效燃油消耗的混合动力汽车控制方法 |
GB2517469A (en) * | 2013-08-21 | 2015-02-25 | Jaguar Land Rover Ltd | Hybrid electric vehicle controller and method |
CN103606271B (zh) * | 2013-11-27 | 2015-10-28 | 大连理工大学 | 一种混合动力城市公交车控制方法 |
-
2015
- 2015-09-21 CN CN201510604187.0A patent/CN105128855B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN105128855A (zh) | 2015-12-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105128855B (zh) | 一种双轴并联混合动力城市公交车的控制方法 | |
CN107839468B (zh) | 高充放电效率的混合动力***及其控制方法 | |
CN106080585B (zh) | 一种双行星排式混合动力汽车非线性模型预测控制方法 | |
CN101898557B (zh) | 电动四驱混合动力车辆的控制方法 | |
Geng et al. | A cascaded energy management optimization method of multimode power-split hybrid electric vehicles | |
Nguyen et al. | Optimal sizing and energy management of an electric vehicle powertrain equipped with two motors and multi-gear ratios | |
CN102069804B (zh) | 一种混合动力汽车行驶状态预测控制方法 | |
CN106240336B (zh) | 一种插电式双电机四驱混合动力***及控制方法 | |
CN106394548B (zh) | 一种插电式四驱混合动力汽车分层协调能效控制方法 | |
CN112959992B (zh) | 基于能效分析与效率最优的混合动力汽车能量管理方法 | |
CN104742898A (zh) | 一种分速汇矩式混合动力功率流控制方法 | |
CN102556056A (zh) | 一种混合动力汽车的双模糊能量控制管理*** | |
CN104477051B (zh) | 一种双驱轴双电机纯电动车驱动电机功率差异化匹配方法 | |
CN107512261B (zh) | 基于双动力源协同的并联phev换挡控制方法 | |
CN104175891A (zh) | 纯电动汽车能量回收再生制动控制方法 | |
Wang et al. | Hybrid electric vehicle modeling accuracy verification and global optimal control algorithm research | |
CN110348063A (zh) | 一种混合动力汽车控制策略的联合仿真方法 | |
Zhang et al. | Powertrain design and energy management of a novel coaxial series-parallel plug-in hybrid electric vehicle | |
CN110667564B (zh) | 并联式混合动力汽车自主队列行驶能量智能管理方法 | |
CN103568814B (zh) | 一种混合动力汽车的驱动*** | |
Esfahanian et al. | Hybrid electric haulage trucks for open pit mining | |
Dong et al. | Performance comparison of series–parallel hybrid transmissions with multiple gears and modes based on efficiency model | |
CN107499303A (zh) | 一种基于双动力源协同的并联phev三参数换挡控制方法 | |
Vajedi et al. | Traction-motor power ratio and speed trajectory optimization for power split PHEVs using route information | |
Sciarretta et al. | Fundamentals of vehicle modeling |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170721 Termination date: 20200921 |
|
CF01 | Termination of patent right due to non-payment of annual fee |