CN114435369A - 一种混合动力汽车能耗优化方法 - Google Patents
一种混合动力汽车能耗优化方法 Download PDFInfo
- Publication number
- CN114435369A CN114435369A CN202210129337.7A CN202210129337A CN114435369A CN 114435369 A CN114435369 A CN 114435369A CN 202210129337 A CN202210129337 A CN 202210129337A CN 114435369 A CN114435369 A CN 114435369A
- Authority
- CN
- China
- Prior art keywords
- individual
- torque
- battery
- motor
- population
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 61
- 238000005457 optimization Methods 0.000 title claims abstract description 53
- 238000005265 energy consumption Methods 0.000 title claims abstract description 42
- 238000004364 calculation method Methods 0.000 claims abstract description 21
- 230000002068 genetic effect Effects 0.000 claims abstract description 21
- 230000035772 mutation Effects 0.000 claims abstract description 10
- 239000000446 fuel Substances 0.000 claims description 26
- 238000010248 power generation Methods 0.000 claims description 14
- 230000008569 process Effects 0.000 claims description 10
- 238000012217 deletion Methods 0.000 claims description 7
- 230000037430 deletion Effects 0.000 claims description 7
- 230000009467 reduction Effects 0.000 claims description 7
- 239000003638 chemical reducing agent Substances 0.000 claims description 6
- 230000001186 cumulative effect Effects 0.000 claims description 6
- 238000007599 discharging Methods 0.000 claims description 6
- OIGNJSKKLXVSLS-VWUMJDOOSA-N prednisolone Chemical compound O=C1C=C[C@]2(C)[C@H]3[C@@H](O)C[C@](C)([C@@](CC4)(O)C(=O)CO)[C@@H]4[C@@H]3CCC2=C1 OIGNJSKKLXVSLS-VWUMJDOOSA-N 0.000 claims description 5
- 230000005611 electricity Effects 0.000 claims description 4
- 230000005540 biological transmission Effects 0.000 claims description 3
- 239000000295 fuel oil Substances 0.000 claims description 3
- 239000011159 matrix material Substances 0.000 description 3
- 230000001174 ascending effect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000003912 environmental pollution Methods 0.000 description 2
- 238000007726 management method Methods 0.000 description 2
- 239000003921 oil Substances 0.000 description 2
- 241001441754 Heteralocha acutirostris Species 0.000 description 1
- 230000001133 acceleration Effects 0.000 description 1
- 230000003044 adaptive effect Effects 0.000 description 1
- 238000007792 addition Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000011084 recovery Methods 0.000 description 1
Images
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
- B60W40/00—Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models
-
- 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
- B60W40/00—Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models
- B60W40/10—Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models related to vehicle motion
-
- 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
-
- 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/0004—In digital systems, e.g. discrete-time systems involving sampling
- B60W2050/0005—Processor details or data handling, e.g. memory registers or chip architecture
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/60—Other road transportation technologies with climate change mitigation effect
- Y02T10/62—Hybrid vehicles
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/80—Technologies aiming to reduce greenhouse gasses emissions common to all road transportation technologies
- Y02T10/84—Data processing systems or methods, management, administration
Landscapes
- Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Transportation (AREA)
- Mechanical Engineering (AREA)
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Human Computer Interaction (AREA)
- Electric Propulsion And Braking For Vehicles (AREA)
- Hybrid Electric Vehicles (AREA)
Abstract
本发明提供了一种混合动力汽车能耗优化方法,属于汽车技术领域。它解决了现有的策略很难保证全局的能耗最低的问题。本混合动力汽车能耗优化方法包括:利用遗传算法将车辆运行参数进行初始化、编码,进而随机产生N个个体,组成初始种群;利用适应度函数计算初始种群中每个个体的适应度值;判断优化代数是否达到预设的最大优化代数,若为是,停止计算,取初始种群中适应度值最小的个体为最优个体,并解码输出;否则,根据选择概率从初始种群中随机选择一个个体进行交叉、变异操作,获得新一代种群并进行适应度值计算,将优化代数加1并返回优化代数判断。本发明能够实现对车辆参数的全局优化,使车辆运行于最低能耗点。
Description
技术领域
本发明属于汽车技术领域,涉及一种混合动力汽车能耗优化方法。
背景技术
面对能源危机和环境污染,开发新能源和相关节能技术以减少对不可再生能源的依赖和环境污染已成为全球发展共识。而混合动力汽车与传统汽车相比,其节油潜力巨大。首先混合动力汽车可依据驾驶员行车需求保证发动机工作在最佳工况范围,发动机有效功率高,油耗控制在最低的范围内;其次驱动电机可利用车辆制动、下坡、发动机剩余功率实现能量回收,减少了能源损耗;在起步时可利用电机单独驱动,避免了发动机怠速或低负荷工况下油耗过高的问题。
而目前的三电机混合动力系统是由发动机以及电机P1、P2、P4这四个动力源组成,由发动机+3DHT混动变速箱(内置P1+P2电机)组成前驱系统,由P4电机+固定速比减速驱动桥组成后驱系统,该混合动力系统具有更好的四驱驾驶体验,车辆的动力性与通过性也更好,但由于存在四个动力源,结构复杂,驾驶模式较多,自由度大,而现有混合动力汽车的控制都是基于既定规则的能量管理策略,其主要是依据工程经验制定的控制规则,但是城市道路往往由于城市路况不同所处的工况也千差万别,实际运营过程中,每趟线路的载客量、城市道路维修情况、早晚上下班高峰期等影响因素都会对车辆的行驶产生一定的影响。因此,基于既定规则的能量管理策略会存在不合理情况,且难以达到理想的效果,很难保证全局的能耗最低。
发明内容
本发明的目的是针对现有技术存在的上述问题,提出了一种混合动力汽车能耗优化方法,其所要解决的技术问题是:如何实现对车辆参数的全局优化,使车辆运行于最低能耗点。
本发明的目的可通过下列技术方案来实现:一种混合动力汽车能耗优化方法,包括如下步骤:
A、根据车辆运行工况获得轮端需求扭矩;
B、基于轮端需求扭矩,利用遗传算法实时生成车辆运行参数,所述车辆运行参数包括工作模式、挡位、P1电机扭矩、P2电机扭矩以及发动机转速和扭矩,将车辆运行参数进行初始化、编码,进而随机产生N个个体,组成初始种群;
C、利用适应度函数计算初始种群中每个个体的适应度值,并记录为1次优化代数;
D、判断优化代数是否达到预设的最大优化代数,若为是,停止计算,取初始种群中适应度值最小的个体,即总燃油消耗量最小的个体为最优个体,并解码输出该最优个体所对应的车辆运行参数;否则,进入步骤E;
E、根据选择概率,从初始种群中随机选择一个个体进行交叉、变异操作,获得新的个体,新的个体组成新一代种群,利用适应度函数计算新一代种群中每个个体的适应度值,并将优化代数加1,之后返回步骤D。
本混合动力汽车能耗优化方法的工作原理为:首先实时监测车辆运行工况与电池剩余电量,车辆运行工况包括起步、加速、等速、减速、转弯、上下坡、停车等,根据车辆运行工况与电池剩余电量,获得轮端需求扭矩,之后基于轮端需求扭矩,利用遗传算法实时生成车辆运行参数,车辆运行参数如工作模式、挡位、P1电机扭矩、P2电机扭矩、P4电机扭矩、发动机转速以及扭矩确定后,将车辆运行参数进行初始化、编码,进而随机产生N个个体,该N个个体组成初始种群,之后利用适应度函数对初始种群中的每个个体进行求解获得适应度值,此时记录为1次优化代数,将该优化代数与最大优化代数进行比较,在优化代数小于最大优化代数时,根据选择概率,从初始种群中随机选择一个个体进行交叉、变异操作,获得新的个体,新的个体组成新一代种群,将新一代种群返回适应度值计算,并将优化代数加1,重新判断优化代数是否达到最大优化代数;在优化代数达到最大优化代数时,停止计算,取初始种群中适应度值最小的个体,即总燃油消耗量最小的个体为最优个体,并解码输出该最优个体所对应的车辆运行参数,包括工作模式、挡位、P1电机扭矩、P2电机扭矩、P4电机扭矩以及发动机转速和扭矩。本方法基于自适应遗传算法对车辆的运行模式、档位、扭矩分配等参数进行优化,实现了车辆参数的全局优化,使车辆运行于最低能耗点。
在上述的混合动力汽车能耗优化方法中,在所述步骤B中,还包括基于过往时刻的挡位与工作模式进行预规划的操作,包括:
标定时间限值T_eng和时间限值T_gear;
在换挡时间T_gearact小于时间限值T_gear时,不允许更换,保持挡位为换挡开始时刻挡位值,不进入遗传算法的寻优操作;
在发动机运行时间T_engrun小于时间限值T_eng时,发动机继续处于RUN阶段,不更换工作模式,不进入遗传算法的寻优操作。其中,时间限值T_eng和时间限值T_gear均可设置为5s。这一步的操作,避免频繁切换挡位或者发动机频繁起机等带来的车辆NVH问题,NVH指Noise(噪声)、Vibration(振动)和Harshness(声振粗糙度),也可以通俗地理解为不平顺性。
在上述的混合动力汽车能耗优化方法中,在所述步骤B中,利用遗传算法将车辆运行参数进行初始化、编码的操作包括:
初始化车辆运行参数,设置种群大小为N,变异率为Pm,交叉率为Pc,最大优化代数为Tmax,寻优区间[Lb,Ub];
采用自适应二进制编码方法对车辆运行参数进行编码,在利用遗传算法优化过程中动态缩小寻优区间,操作为:
将第k次优化代数所获得的最优个体适应度值与当前种群的适应度平均值进行差值计算,获得差值ΔF,在差值ΔF小于预设值Ck时,以当前种群的适应度最小的值X_best作为中点,然后将寻优区间缩小i倍,即[X_best-(Ub-Lb)/i,X_best+(Ub-Lb)/i]。寻优区间[Lb,Ub]是一个矩阵,每一个参数的寻优区间都在里面,是个多阶的矩阵。将寻优区间缩小,结果将逐渐收敛于全局最优点。其中,在本方法中,将N设置为100;Pm设置为80%;Pc设置为80%;Tmax设置为1O0;k设置为50,Ck设置为0.5;i设置为10。
在上述的混合动力汽车能耗优化方法中,在所述步骤C中,利用适应度函数计算初始种群中每个个体的适应度值的操作包括:
其中,QElc为电池电能消耗;ηbat,charge为电池发电效率;ηgen_p1为P1电机发电效率;ηgen_P1_inverter为P1电机发电时inverter效率;ηeng为发动机热效率;VHL为燃油热值。
在上述的混合动力汽车能耗优化方法中,在所述步骤C中,所述电池电能消耗QElc由电源电流Ibat来确定;所述电源电流Ibat由以下公式获得:
Ibat=IP1+IP2+IP4+ILS
其中:Ibat为电源电流;IP1为P1电机电流;IP2为P2电机电流;IP4为P4电机电流;ILS是控制器低压负载电流,为常数;
在电源电流Ibat大于0时,电池处于放电状态,则电池电能消耗QElc为:
在电源电流Ibat小于0时,电池处于充电状态,则电池电能消耗QElc为:
QElc=Ibat×U×ηbat,charge
其中,U为回路电压;ηbat,charge为电池发电效率。
在上述的混合动力汽车能耗优化方法中,在所述步骤C中,还包括利用惩罚函数对种群中的个体进行优化删选的操作,包括:
在车辆处于纯电工作模式下时,判断种群个体中的发动机扭矩Tqeng是否大于零,在发动机扭矩Tqeng大于零时,惩罚去除该个体;
在车辆处于串联工作模式下时,判断种群个体中的发动机扭矩Tqeng是否等于P1电机扭矩TqP1,在发动机扭矩Tqeng不等于P1电机扭矩TqP1时,惩罚去除该个体;
在车辆处于并联工作模式下时,判断种群个体中的发动机转速Neng是否等于P2电机转速NP2,在发动机转速Neng不等于P2电机转速NP2时,惩罚去除该个体;
反之,在发动机扭矩Tqeng等于零、发动机扭矩Tqeng等于P1电机扭矩TqP1以及发动机转速Neng等于P2电机转速NP2时,输出该个体的适应度值。惩罚函数主要是为了去除一些不合理的情况,比如车辆如果处于纯电动模式下,发动机存在正扭矩,那么这一组数据就是不合理的,就要把它给舍弃,不让结果算到里面,这样在优化的过程中就会排除这些情况,提高优化结果的准确性。
在上述的混合动力汽车能耗优化方法中,在所述步骤C中,利用惩罚函数对种群中的个体进行优化删选的操作,还包括:
利用传动模型计算P4电机扭矩TqP4与转速NP4以及P2电机转速NP2;
预设扭矩最大限值Tp_P4_max和扭矩最小限值Tp_P4_min;
判断P4电机扭矩TqP4是否小于扭矩最小限值Tp_P4_min,或者P4电机扭矩TqP4是否大于扭矩最大限值Tp_P4_max,或者P4电机扭矩TqP4与P2电机扭矩TqP2的乘积是否小于零;
在上述任一条件为是时,启动惩罚函数进行惩罚去除;反之,输出适应度值;
所述P4电机扭矩TqP4由以下公式计算获得:
其中:iGear为DHT减速比;iFinal为后轴主减速器减速比;Tqwhl为轮端驱动扭矩;Tqeng为发动机扭矩;TqP1为P1电机扭矩;TqP2为P2电机扭矩;Tq_transloss为扭矩损失,公式为:
Tqtransloss=(Tqeng+TqP1+TqP2)×ηtranns×iGear+TqP4×ηfinal×iFinal
其中:η_trans为DHT效率;ηfinal为后轴主减速器效率。
在上述的混合动力汽车能耗优化方法中,在所述步骤C中,利用惩罚函数对种群中的个体进行优化删选的操作,还包括:
预设最小电池电量限值SOCmin和最大电池电量限值SOCmax;
判断电池剩余电量SOC是否在最小电池电量限值SOCmin和最大电池电量限值SOCmax的区间内,若是,输出该个体的适应度值;反之,利用惩罚函数去除该个体;
所述电池剩余电量SOC的计算包括:
在电源电流Ibat小于0时,电池处于充电状态,电池在第i时刻的剩余电量SOCi为:
在电源电流Ibat大于0时,电池处于放电状态,电池在i时刻的剩余电量SOCi为:
其中:SOCi-1为电池在i的上一个时刻的剩余电量;ηbat,discharqe为电池放电效率;Ebat为电池最大容量;U为回路电压;ηbat,charge为电池发电效率。其中,最小电池电量限值SOCmin可设置为总电量的15%,最大电池电量限值SOCmax可设置为总电量的95%。判断电池剩余电量SOC是否在最小电池电量限值SOCmin和最大电池电量限值SOCmax的区间内的操作,能够防止电池SOC过低以及电池过充的问题。
在上述的混合动力汽车能耗优化方法中,在所述步骤D中,判断是否停止计算的操作还包括:
预设适应度函数偏差阈值为TolFun;
将每个个体的适应度值进行平均值计算获得适应度平均值;
将适应度平均值与最优个体的适应度值进行差值计算获得适应度函数偏差;
在判断适应度函数偏差小于TolFun,则停止计算,取初始种群中适应度值最小的个体作为最优个体,解码输出该最优个体所对应的车辆运行参数;否则,进入步骤E。
在上述的混合动力汽车能耗优化方法中,在所述步骤E中,从初始种群中随机选择一个个体进行交叉、变异的操作包括:
E01、计算每一个个体的选择概率为:
E02、基于选择概率pk计算每个个体的累计概率:
其中,qk为第k个个体的累计概率;
E03、在[0,1]内产生一个均匀分布的随机数r;
E04、若qk-1<r<qk,则选择个体k,重复M次步骤E03和E04;
E05、根据交叉率Pc进行两两交叉处理,获得交叉种群,两两交叉处理的操作为:
akj′=akj(1-b)+aljb
alj′=alj(1-b)+akjb其中:alj为第l个个体的第j位;akj为第k个个体的第j位;akj′为交叉处理后第k个个体的第j位;alj′为交叉处理后第1个个体的第j位;b∈[0,1];
E06、从交叉种群中随机选择一个个体,按变异率Pm变异得到新的个体,新的个体组成新一代种群,再将其返回步骤C进行适应度计算。
与现有技术相比,本混合动力汽车能耗优化方法具有以下优点:
1、本发明将运行模式也作为优化目标,而非按照既定规则选择运行模式,而且本发明在通过遗传算法生成车辆的运行模式、变速箱档位、扭矩分配、转速分配等之后,还通过添加惩罚函数,对遗传算法生成的状态进行筛选与淘汰,可以避免一些不合理的工况的出现,如:发动机的突然起机,变速箱档位的突变等,实现了对车辆参数的全局优化,优化效果更好,优化结果更合理。
2、本发明在遗传算法寻优时加入自适应编码代替传统的二进制编码方式,计算过程中动态缩小优化区间,减少迭代次数,提高了寻优时的收敛次数,提高了遗传算法的搜索效率。
附图说明
图1是本发明实施例一的控制流程图。
图2是本发明实施例二的控制流程图。
具体实施方式
以下是本发明的具体实施例并结合附图,对本发明的技术方案作进一步的描述,但本发明并不限于这些实施例。
实施例一:
如图1所示,车辆在行驶过程中,本混合动力汽车能耗优化方法启动,实时监测车辆运行工况,如起步、加速、等速、减速、转弯、上下坡、停车等参数,同时实时监测电池剩余电量,根据车辆运行工况和电池剩余电量获得轮端需求扭矩;具体为:在车辆行驶时,当存在加速需求时,根据加速踏板的开度与车辆时速,结合标定的加速踏板MAP获得轮端需求扭矩。当存在减速需求时,根据制动踏板的开度与车辆时速,结合标定的制动踏板MAP获得轮端需求扭矩。
之后根据电池剩余电量,预规划车辆的运行模式,若电池剩余电量低于最小电池电量限值SOCmin,则退出纯电模式,通过P1电机给电池充电。
之后进入基于过往时刻的挡位与工作模式进行预规划的操作,包括:
标定时间限值T_eng和时间限值T_gear;
在换挡时间T_gearact小于时间限值T_gear时,不允许更换,保持挡位为换挡开始时刻挡位值,该参数不进入遗传算法的寻优操作;
在发动机运行时间T_engrun小于时间限值T_eng时,发动机继续处于RUN阶段,不更换工作模式,该参数不进入遗传算法的寻优操作。
通过上面两步的预规划,避免频繁切换挡位或者发动机频繁起机等带来的车辆NVH问题,也能防止电量过放的问题,还能提高寻优效率。
之后基于轮端需求扭矩,利用遗传算法实时生成车辆运行参数,车辆运行参数包括工作模式、挡位、P1电机扭矩、P2电机扭矩以及发动机转速和扭矩,将车辆运行参数进行初始化、编码,若上一步骤中挡位不允许更换,则该参数不进行初始化和编码操作;
利用传动模型计算P4电机扭矩TqP4与转速NP4,P2电机转速NP2;
本车型动力总成系统由发动机与P1电机、P2电机、P4电机,加上DHT与主减速器组成,可得出轮端驱动扭矩Tqwnl方程:
Tqwnl=(Tqeng+TqP1+TqP2)×iGear+TqP4×iFinal+Tq_transloss
其中:iGear为DHT减速比;iFinal为后轴主减速器减速比;Tq_transloss为扭矩损失;Tqeng为发动机扭矩;TqP1为P1电机扭矩;TqP2为P2电机扭矩;
则P4电机扭矩TqP4为:
轮端转速Nwhl可由车速与轮胎半径计算获得:
Nwhl=vi/3.6*60/2/π/R
可知,P4电机转速NP4,P2电机转速NP2为:
NP4=Nwhl×iFinal
NP2=Nwhl×iGear
初始化和编码的操作为:
初始化车辆运行参数,包括工作模式、挡位、P1电机扭矩、P2电机扭矩、P4电机扭矩以及发动机转速和扭矩,设置种群大小为N,变异率为Pm,交叉率为Pc,最大优化代数为Tmax,寻优区间[Lb,Ub];
采用自适应二进制编码方法对车辆运行参数进行编码,在利用遗传算法优化过程中动态缩小寻优区间,操作为:
将优化代数为k时获得的最优个体的适应度值与当前种群的适应度平均值进行差值计算,获得差值ΔF,在差值ΔF小于预设值Ck时,以当前种群的适应度最小的值X_best作为中点,然后将寻优区间缩小i倍,即[X_best-(Ub-Lb)/i,X_best+(Ub-Lb)/i]。
经过初始化和编码操作后,在寻优区间[Lb,Ub]内随机产生N个个体,组成初始种群;比如令初始种群为V,则A={a1,a2,…ak,a1…aN},其中第k个个体为ak=(ak1,ak2,…akj),ak1,ak2,…aki分别表示第k个个体的工作模式、挡位、P1电机扭矩、P2电机扭矩、P4电机扭矩以及发动机转速和扭矩,将当前优化代数设置为T=1:
之后利用适应度函数计算初始种群中每个个体的适应度值,并记录优化代数为1次;具体为:
在电机处于Q1驱动阶段时,Q1驱动阶段是指电机加速前进状态,此时电机的电流I为:
I=Tq*N/9550/ηmot/ηmot_inverter/U
在电机处于Q4发电阶段时,Q4发电阶段为前进制动状态,电机的电流为:
I=Tq*N/9550/ηgen/ηgen_inverter/U
其中:Tq为电机扭矩;N为电机转速;ηmot为电机驱动效率;ηmot_inverter为电机驱动时inverter效率;ηgen为电机发电效率;ηgen_inverter为电机发电时inverter效率;U为回路电压。
将P2电机扭矩和P2电机转速代入上式即可获得P2电机电流IP2;将P4电机扭矩和P4电机转速代入上式即可获得P4电机电流IP4;
P1电机主要应用在发电状态,因此P1电机扭矩和转速与发动机扭矩和转速一致,由此也可获得P1电机电流IP1。
由此可得出,电源电流Ibat为:
Ibat=IP1+IP2+IP4+ILS
其中:Ibat为电源电流;IP1为P1电机电流;IP2为P2电机电流;IP4为P4电机电流;ILS是控制器低压负载电流,为常数;
电池电能消耗QElc的计算操作为:
在电源电流Ibat大于0时,电池处于放电状态,则电池电能消耗QElc为:
在电源电流Ibat小于0时,电池处于充电状态,则电池电能消耗QElc为:
QElc=Ibat×U×ηbat,charge
其中,U为回路电压;ηbat,charge为电池发电效率。
其中,QElc为电池电能消耗;ηbat,charge为电池发电效率;ηgen_p1为P1电机发电效率;ηgen_P1_inverter为P1电机发电时inverter效率;ηeng为发动机热效率;VHL为燃油热值。
在计算每个个体的适应度值时,还包括利用惩罚函数对种群中的个体进行优化删选的操作,包括:
在车辆处于纯电工作模式下时,判断种群个体中的发动机扭矩Tqeng是否大于零,在发动机扭矩Tqeng大于零时,惩罚去除该个体;
在车辆处于串联工作模式下时,判断种群个体中的发动机扭矩Tqeng是否等于P1电机扭矩TqP1,在发动机扭矩Tqeng不等于P1电机扭矩TqP1时,惩罚去除该个体;
在车辆处于并联工作模式下时,判断种群个体中的发动机转速Neng是否等于P2电机转速NP2,在发动机转速Neng不等于P2电机转速NP2时,惩罚去除该个体;
反之,在发动机扭矩Tqeng等于零、发动机扭矩Tqeng等于P1电机扭矩TqP1以及发动机转速Neng等于P2电机转速NP2时,输出该个体的适应度值。
还包括:预设扭矩最大限值Tp_P4_max和扭矩最小限值Tp_P4_min;
判断P4电机扭矩TqP4是否小于扭矩最小限值Tp_P4_min,或者P4电机扭矩TqP4是否大于扭矩最大限值Tp_P4_max,或者P4电机扭矩TqP4与P2电机扭矩TqP2的乘积是否小于零;
在上述任一条件为是时,启动惩罚函数进行惩罚去除;反之,输出适应度值;
还包括:预设最小电池电量限值SOCmin和最大电池电量限值SOCmax;
判断电池剩余电量SOC是否在最小电池电量限值SOCmin和最大电池电量限值SOCmax的区间内,若是,输出该个体的适应度值;反之,利用惩罚函数去除该个体;
电池剩余电量SOC的计算包括:
在电源电流Ibat小于0时,电池处于充电状态,电池在第i时刻的剩余电量SOCi为:
在电源电流Ibat大于0时,电池处于放电状态,电池在i时刻的剩余电量SOCi为:
其中:SOCi-1为电池在i的上一个时刻的剩余电量;ηbat,discharge为电池放电效率;Ebat为电池最大容量;U为回路电压;ηbat,charge为电池发电效率。
之后判断是否停止计算;首先判断优化代数是否达到预设的最大优化代数Tmax,若为是,停止计算,取初始种群中适应度值最小的个体,即总燃油消耗量最小的个体为最优个体,并解码输出该最优个体所对应的车辆运行参数;
否则,根据选择概率,从初始种群中随机选择一个个体进行交叉、变异操作,获得新的个体,新的个体组成新一代种群,对新一代种群进行适应度值计算,并将优化代数加1,之后重新判断优化代数是否达到最大优化代数Tmax。新一代种群组成的具体操作为:
首先计算每一个个体的选择概率为:
再基于选择概率pk计算每个个体的累计概率:
其中,qk为第k个个体的累计概率;
之后在[0,1]内产生一个均匀分布的随机数r;
若qk-1<r<qk,则选择个体k,重复M次该步骤和上一步;
之后根据交叉率Pc进行两两交叉处理,获得交叉种群,两两交叉处理的操作为:
akj′=akj(1-b)+aljb
alj′=alj(1-b)+akjb
其中:alj为第1个个体的第j位;akj为第k个个体的第j位;akj′为交叉处理后第k个个体的第j位;alj′为交叉处理后第1个个体的第j位;b∈[0,1];b表示0-1之间随机的交叉位置。
这里涉及到两个二进制数的部分位的交换,比如:
个体A:1001↑111→1001000新个体
个体B:0011↑000→0011111新个体
从交叉种群中随机选择一个个体,按变异率Pm变异得到新的个体,新的个体组成新一代种群,再返回对每个个体进行适应度计算的步骤,由此获取车辆运行参数的最优结果。
本方法优化了四个动力源的扭矩分配,在优化P1、P2、P4电机扭矩分配的同时,可以兼顾发动机的运行工况,将发动机尽可能地工作在最佳运行线附近,降低发动机的燃油消耗。
实施例二:
如图2所示,本实施例中的技术方案与实施例一中的技术方案基本相同,不同之处在于,判断停止计算的操作为:
预设适应度函数偏差阈值为TolFun;
将每个个体的适应度值进行平均值计算获得适应度平均值;
将适应度平均值与最优个体的适应度值进行差值计算获得适应度函数偏差;
在判断适应度函数偏差小于TolFun,则停止计算,取初始种群中适应度值最小的个体作为最优个体,解码输出该最优个体所对应的车辆运行参数;否则根据选择概率,从初始种群中随机选择一个个体进行交叉、变异操作,获得新的个体,新的个体组成新一代种群,对新一代种群进行适应度值计算,并将优化代数加1,之后重新判断优化代数是否达到最大优化代数Tmax。
本文中所描述的具体实施例仅仅是对本发明精神作举例说明。本发明所属技术领域的技术人员可以对所描述的具体实施例做各种各样的修改或补充或采用类似的方式替代,但并不会偏离本发明的精神或者超越所附权利要求书所定义的范围。
Claims (10)
1.一种混合动力汽车能耗优化方法,其特征在于,包括如下步骤:
A、根据车辆运行工况获得轮端需求扭矩;
B、基于轮端需求扭矩,利用遗传算法实时生成车辆运行参数,所述车辆运行参数包括工作模式、挡位、P1电机扭矩、P2电机扭矩以及发动机转速和扭矩,将车辆运行参数进行初始化、编码,进而随机产生N个个体,组成初始种群;
C、利用适应度函数计算初始种群中每个个体的适应度值,并记录为1次优化代数;
D、判断优化代数是否达到预设的最大优化代数,若为是,停止计算,取初始种群中适应度值最小的个体,即总燃油消耗量最小的个体为最优个体,并解码输出该最优个体所对应的车辆运行参数;否则,进入步骤E;
E、根据选择概率,从初始种群中随机选择一个个体进行交叉、变异操作,获得新的个体,新的个体组成新一代种群,利用适应度函数计算新一代种群中每个个体的适应度值,并将优化代数加1,之后返回步骤D。
2.根据权利要求1所述的混合动力汽车能耗优化方法,其特征在于,在所述步骤B中,还包括基于过往时刻的挡位与工作模式进行预规划的操作,包括:
标定时间限值T_eng和时间限值T_gear;
在换挡时间T_gearact小于时间限值T_gear时,不允许更换,保持挡位为换挡开始时刻挡位值,该参数不进入遗传算法的寻优操作;
在发动机运行时间T_engrun小于时间限值T_eng时,发动机继续处于RUN阶段,不更换工作模式,该参数不进入遗传算法的寻优操作。
3.根据权利要求1或2所述的混合动力汽车能耗优化方法,其特征在于,在所述步骤B中,利用遗传算法将车辆运行参数进行初始化、编码的操作包括:
初始化车辆运行参数,设置种群大小为N,变异率为Pm,交叉率为Pc,最大优化代数为Tmax,寻优区间[Lb,Ub];
采用自适应二进制编码方法对车辆运行参数进行编码,在利用遗传算法优化过程中动态缩小寻优区间,操作为:
将第k次优化代数所获得的最优个体适应度值与当前种群的适应度平均值进行差值计算,获得差值ΔF,在差值ΔF小于预设值Ck时,以当前种群的适应度最小的值X_best作为中点,然后将寻优区间缩小i倍,即[X_best一(Ub-Lb)/i,X_best+(Ub-Lb)/i]。
6.根据权利要求1所述的混合动力汽车能耗优化方法,其特征在于,在所述步骤C中,还包括利用惩罚函数对种群中的个体进行优化删选的操作,包括:
在车辆处于纯电工作模式下时,判断种群个体中的发动机扭矩Tqeng是否大于零,在发动机扭矩Tqeng大于零时,惩罚去除该个体;
在车辆处于串联工作模式下时,判断种群个体中的发动机扭矩Tqeng是否等于P1电机扭矩TqP1,在发动机扭矩Tqeng不等于P1电机扭矩TqP1时,惩罚去除该个体;
在车辆处于并联工作模式下时,判断种群个体中的发动机转速Neng是否等于P2电机转速NP2,在发动机转速Neng不等于P2电机转速NP2时,惩罚去除该个体;
反之,在发动机扭矩Tqeng等于零、发动机扭矩Tqeng等于P1电机扭矩TqP1以及发动机转速Neng等于P2电机转速NP2时,输出该个体的适应度值。
7.根据权利要求6所述的混合动力汽车能耗优化方法,其特征在于,在所述步骤C中,利用惩罚函数对种群中的个体进行优化删选的操作,还包括:
利用传动模型计算P4电机扭矩TqP4与转速NP4以及P2电机转速NP2;
预设扭矩最大限值Tp_P4_max和扭矩最小限值Tp_P4_min;
判断P4电机扭矩TqP4是否小于扭矩最小限值Tp_P4_min,或者P4电机扭矩TqP4是否大于扭矩最大限值Tp_P4_max,或者P4电机扭矩TqP4与P2电机扭矩TqP2的乘积是否小于零;
在上述任一条件为是时,启动惩罚函数进行惩罚去除;反之,输出适应度值;
所述P4电机扭矩TqP4由以下公式计算获得:
其中:iGear为DHT减速比;iFinal为后轴主减速器减速比;Tqwhl为轮端驱动扭矩;Tqeng为发动机扭矩;TqP1为P1电机扭矩;TqP2为P2电机扭矩;Tq_transloss为扭矩损失,公式为:
Tqtransloss=(Tqeng+TqP1+TqP2)×ηtrans×iGear+TqP4×ηfinal×iFinal
其中:η_trans为DHT效率;ηfinal为后轴主减速器效率。
8.根据权利要求6或7所述的混合动力汽车能耗优化方法,其特征在于,在所述步骤C中,利用惩罚函数对种群中的个体进行优化删选的操作,还包括:
预设最小电池电量限值SOCmin和最大电池电量限值SOCmax;
判断电池剩余电量SOC是否在最小电池电量限值SOCmin和最大电池电量限值SOCmax的区间内,若是,输出该个体的适应度值;反之,利用惩罚函数去除该个体;
所述电池剩余电量SOC的计算包括:
在电源电流Ibat小于0时,电池处于充电状态,电池在第i时刻的剩余电量SOCi为:
在电源电流Ibat大于0时,电池处于放电状态,电池在i时刻的剩余电量SOCi为:
其中:SOCi-1为电池在i的上一个时刻的剩余电量;ηbat,discharge为电池放电效率;Ebat为电池最大容量;U为回路电压;ηbat,charge为电池发电效率。
9.根据权利要求1或2所述的混合动力汽车能耗优化方法,其特征在于,在所述步骤D中,判断是否停止计算的操作还包括:
预设适应度函数偏差阈值为TolFun;
将每个个体的适应度值进行平均值计算获得适应度平均值;
将适应度平均值与最优个体的适应度值进行差值计算获得适应度函数偏差;
在判断适应度函数偏差小于TolFun,则停止计算,取初始种群中适应度值最小的个体作为最优个体,解码输出该最优个体所对应的车辆运行参数;否则,进入步骤E。
10.根据权利要求1或2所述的混合动力汽车能耗优化方法,其特征在于,在所述步骤E中,从初始种群中随机选择一个个体进行交叉、变异的操作包括:
E01、计算每一个个体的选择概率为:
E02、基于选择概率pk计算每个个体的累计概率:
其中,qk为第k个个体的累计概率;
E03、在[0,1]内产生一个均匀分布的随机数r;
E04、若qk-1<r<qk,则选择个体k,重复M次步骤E03和E04;
E05、根据交叉率Pc进行两两交叉处理,获得交叉种群,两两交叉处理的操作为:
akj′=akj(1-b)+aljb
alj′=alj(1-b)+akjb
其中:alj为第1个个体的第j位;akj为第k个个体的第j位;akj′为交叉处理后第k个个体的第j位;alj′为交叉处理后第1个个体的第j位;b∈[0,1];
E06、从交叉种群中随机选择一个个体,按变异率Pm变异得到新的个体,新的个体组成新一代种群,再将其返回步骤C进行适应度计算。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210129337.7A CN114435369B (zh) | 2022-02-11 | 2022-02-11 | 一种混合动力汽车能耗优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210129337.7A CN114435369B (zh) | 2022-02-11 | 2022-02-11 | 一种混合动力汽车能耗优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114435369A true CN114435369A (zh) | 2022-05-06 |
CN114435369B CN114435369B (zh) | 2023-12-22 |
Family
ID=81370986
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210129337.7A Active CN114435369B (zh) | 2022-02-11 | 2022-02-11 | 一种混合动力汽车能耗优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114435369B (zh) |
Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016101536A1 (zh) * | 2014-12-22 | 2016-06-30 | 中兴通讯股份有限公司 | 一种并联交错buck变换器及控制方法 |
US20160304080A1 (en) * | 2015-04-17 | 2016-10-20 | Mitsubishi Electric Corporation | Vehicle energy management device |
CN108000501A (zh) * | 2017-11-22 | 2018-05-08 | 湖北工业大学 | 一种用于串联机器人的新型轨迹规划方法 |
CN109346787A (zh) * | 2018-09-21 | 2019-02-15 | 北京机械设备研究所 | 一种电动汽车动力电池自适应优化充电方法 |
CN109849897A (zh) * | 2019-02-28 | 2019-06-07 | 江苏大学 | 一种考虑耦合传动系统动态效率的混合动力能量管理方法 |
CN109981103A (zh) * | 2019-04-19 | 2019-07-05 | 华北水利水电大学 | 一种双二阶广义积分器锁频环的参数优化方法及系统 |
CN110673489A (zh) * | 2019-10-25 | 2020-01-10 | 国网山东省电力公司电力科学研究院 | 一种商用建筑房间的热负荷辨识方法 |
CN110834623A (zh) * | 2019-11-11 | 2020-02-25 | 常熟理工学院 | 一种基于顺序选择遗传算法的全混动汽车能效优化方法 |
CN110834624A (zh) * | 2019-11-11 | 2020-02-25 | 常熟理工学院 | 一种自适应遗传算法的全混动汽车能效优化控制方法 |
CN112590760A (zh) * | 2020-12-22 | 2021-04-02 | 江苏大学 | 一种考虑模式切换频繁度的双电机混合动力汽车能量管理系统 |
CN113007012A (zh) * | 2019-12-19 | 2021-06-22 | 新疆金风科技股份有限公司 | 扭矩控制系数的寻优方法和装置、风力发电机组 |
CN113911101A (zh) * | 2021-10-14 | 2022-01-11 | 燕山大学 | 一种基于同轴并联式结构的在线能量分配方法 |
-
2022
- 2022-02-11 CN CN202210129337.7A patent/CN114435369B/zh active Active
Patent Citations (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2016101536A1 (zh) * | 2014-12-22 | 2016-06-30 | 中兴通讯股份有限公司 | 一种并联交错buck变换器及控制方法 |
US20160304080A1 (en) * | 2015-04-17 | 2016-10-20 | Mitsubishi Electric Corporation | Vehicle energy management device |
CN108000501A (zh) * | 2017-11-22 | 2018-05-08 | 湖北工业大学 | 一种用于串联机器人的新型轨迹规划方法 |
CN109346787A (zh) * | 2018-09-21 | 2019-02-15 | 北京机械设备研究所 | 一种电动汽车动力电池自适应优化充电方法 |
CN109849897A (zh) * | 2019-02-28 | 2019-06-07 | 江苏大学 | 一种考虑耦合传动系统动态效率的混合动力能量管理方法 |
CN109981103A (zh) * | 2019-04-19 | 2019-07-05 | 华北水利水电大学 | 一种双二阶广义积分器锁频环的参数优化方法及系统 |
CN110673489A (zh) * | 2019-10-25 | 2020-01-10 | 国网山东省电力公司电力科学研究院 | 一种商用建筑房间的热负荷辨识方法 |
CN110834623A (zh) * | 2019-11-11 | 2020-02-25 | 常熟理工学院 | 一种基于顺序选择遗传算法的全混动汽车能效优化方法 |
CN110834624A (zh) * | 2019-11-11 | 2020-02-25 | 常熟理工学院 | 一种自适应遗传算法的全混动汽车能效优化控制方法 |
CN113007012A (zh) * | 2019-12-19 | 2021-06-22 | 新疆金风科技股份有限公司 | 扭矩控制系数的寻优方法和装置、风力发电机组 |
CN112590760A (zh) * | 2020-12-22 | 2021-04-02 | 江苏大学 | 一种考虑模式切换频繁度的双电机混合动力汽车能量管理系统 |
CN113911101A (zh) * | 2021-10-14 | 2022-01-11 | 燕山大学 | 一种基于同轴并联式结构的在线能量分配方法 |
Also Published As
Publication number | Publication date |
---|---|
CN114435369B (zh) | 2023-12-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109606348B (zh) | 一种插电式行星混联汽车能量管理控制方法 | |
CN106080585B (zh) | 一种双行星排式混合动力汽车非线性模型预测控制方法 | |
CN111301398B (zh) | 一种cvt插电式混合动力汽车能量管理方法及系统 | |
CN104760591B (zh) | 混合动力综合控制系统 | |
CN104590248B (zh) | 一种基于串并联混合动力系统的控制方法 | |
CN102556056A (zh) | 一种混合动力汽车的双模糊能量控制管理系统 | |
CN101708722A (zh) | 基于模糊逻辑的串联式混合动力电动汽车整车控制方法 | |
CN107458369A (zh) | 一种同轴并联式混合动力电动汽车能量管理方法 | |
CN113911101B (zh) | 一种基于同轴并联式结构的在线能量分配方法 | |
CN107264293A (zh) | 一种基于模糊控制的并联式混合动力车的能量回收方法 | |
CN105128855A (zh) | 一种双轴并联混合动力城市公交车的控制方法 | |
Zhang et al. | Powertrain design and energy management of a novel coaxial series-parallel plug-in hybrid electric vehicle | |
CN110834624B (zh) | 一种自适应遗传算法的全混动汽车能效优化控制方法 | |
CN202413783U (zh) | 一种混合动力汽车的双模糊能量控制管理系统 | |
CN113147412A (zh) | 一种后驱纯电动商用车制动能量回收策略 | |
WO2023098465A1 (zh) | 一种电动汽车制动换挡控制方法及系统 | |
CN112498332B (zh) | 一种并联式混合动力汽车模糊自适应能量管理控制方法 | |
AU2020102967A4 (en) | A parallel – series electric hybrid vehicles based on geometry formulation | |
CN112677771A (zh) | 一种基于模糊控制的前驱电动汽车的再生制动控制方法 | |
CN114435369A (zh) | 一种混合动力汽车能耗优化方法 | |
Liu et al. | Robust energy management strategy for a range extender electric vehicle via genetic algorithm | |
CN105922987B (zh) | Hev/ev再生制动过程中的降挡策略及机电系统 | |
CN111891109B (zh) | 基于非合作博弈论的混合动力汽车能量优化分配控制方法 | |
Zhong et al. | An optimal torque distribution strategy for an integrated starter—generator parallel hybrid electric vehicle based on fuzzy logic control | |
CN113997926A (zh) | 基于分层强化学习的并联式混合动力汽车能量管理方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |