CN113408214A - 油耗与排放一体化优化的商用车队列合并方法 - Google Patents

油耗与排放一体化优化的商用车队列合并方法 Download PDF

Info

Publication number
CN113408214A
CN113408214A CN202110664644.0A CN202110664644A CN113408214A CN 113408214 A CN113408214 A CN 113408214A CN 202110664644 A CN202110664644 A CN 202110664644A CN 113408214 A CN113408214 A CN 113408214A
Authority
CN
China
Prior art keywords
vehicle
vehicles
ith
queue
speed
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
Application number
CN202110664644.0A
Other languages
English (en)
Other versions
CN113408214B (zh
Inventor
刘迪
胡云峰
张辉
宫洵
孙耀
李勇
陈虹
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Jilin University
Original Assignee
Jilin University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Jilin University filed Critical Jilin University
Priority to CN202110664644.0A priority Critical patent/CN113408214B/zh
Publication of CN113408214A publication Critical patent/CN113408214A/zh
Application granted granted Critical
Publication of CN113408214B publication Critical patent/CN113408214B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/28Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • GPHYSICS
    • G08SIGNALLING
    • G08GTRAFFIC CONTROL SYSTEMS
    • G08G1/00Traffic control systems for road vehicles
    • G08G1/22Platooning, i.e. convoy of communicating vehicles
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04LTRANSMISSION OF DIGITAL INFORMATION, e.g. TELEGRAPHIC COMMUNICATION
    • H04L67/00Network arrangements or protocols for supporting network services or applications
    • H04L67/01Protocols
    • H04L67/12Protocols specially adapted for proprietary or special-purpose networking environments, e.g. medical networks, sensor networks, networks in vehicles or remote metering networks
    • H04L67/125Protocols specially adapted for proprietary or special-purpose networking environments, e.g. medical networks, sensor networks, networks in vehicles or remote metering networks involving control of end-device applications over a network
    • HELECTRICITY
    • H04ELECTRIC COMMUNICATION TECHNIQUE
    • H04WWIRELESS COMMUNICATION NETWORKS
    • H04W4/00Services specially adapted for wireless communication networks; Facilities therefor
    • H04W4/30Services specially adapted for particular environments, situations or purposes
    • H04W4/40Services specially adapted for particular environments, situations or purposes for vehicles, e.g. vehicle-to-pedestrians [V2P]
    • H04W4/46Services specially adapted for particular environments, situations or purposes for vehicles, e.g. vehicle-to-pedestrians [V2P] for vehicle-to-vehicle communication [V2V]

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Signal Processing (AREA)
  • Computing Systems (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Mathematical Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Health & Medical Sciences (AREA)
  • General Health & Medical Sciences (AREA)
  • Medical Informatics (AREA)
  • Fluid Mechanics (AREA)
  • Algebra (AREA)
  • Control Of Vehicle Engines Or Engines For Specific Uses (AREA)

Abstract

一种油耗与排放一体化优化的商用车队列合并方法,属于智能交通技术领域。本发明的目的是通过车联网技术,获得同路车辆目的地、车速、加速度等信息,构建综合油耗与排放的优化问题,优化车辆瞬时编队行为的油耗与排放一体化优化的商用车队列合并方法。本发明的步骤是:车速模型搭建、油耗与排放模型、控制目标的确定、问题约束的确定。本发明对参与队列合并的车辆的车速进行了规划,获得了可以使队列整体油耗与排放综合最优的每辆车的车速以及对应该车速的控制输入。

Description

油耗与排放一体化优化的商用车队列合并方法
技术领域
本发明属于智能交通技术领域。
背景技术
车辆智能网联(车联网)是由车辆位置、速度和路线等信息构成的巨大交互网络,车辆可以通过车联网完成自身环境和状态信息的采集以及自身信息的上传。研究表明重型货运车辆以很小间距列队行驶时能有效减少跟随车的空气阻力,从而降低风阻油耗。车联网技术的发展使得车辆在编队行驶时能够通过车辆与车辆信息交互(V2V),车辆与网络信息交互(V2N),车辆与交通设施信息交互(V2I)获得道路信息,同路车辆的加速度、速度、位置等等智能信息,从而使用这些广泛超前的信息优化车辆编队行为,获得车辆在编队时整体油耗与排放综合最优的效果。现有研究在对车辆编队控制时,具有以下缺点:(1)为了简化问题,将时域问题转化为距离域问题,只得到车辆的平均最优速度,没有优化车辆的瞬时速度。(2)排放只是作为降低油耗带来的附加效益,没有考虑车辆行驶时排放和油耗的综合控制。(3)为了最大限度地减小风阻,采用了风阻最小的距离,但没有考虑车辆行驶的安全性。(4)在车辆层面优化中没有添加发动机系统信息,可能导致发动机系统不能满足更高层次优化的要求或增加发动机控制系统的负担。
发明内容
本发明的目的是通过车联网技术,获得同路车辆目的地、车速、加速度等信息,构建综合油耗与排放的优化问题,优化车辆瞬时编队行为的油耗与排放一体化优化的商用车队列合并方法。
本发明的步骤是:
S1、车速模型搭建
以柴油商用车为研究对象,当道路上有N辆车准备进行编队行驶时,选取行驶位移和车速为状态变量,构建参与编队的第i辆车的纵向动力学系统方程如下
Figure BDA0003116346470000011
Figure BDA0003116346470000012
其中i代表参与编队的第i辆车,各类符号分别为第i辆车的:si为第i辆车的车辆位移,vi为第i辆车的车速,Ft,i为第i辆车的车辆控制输入力,Mi为第i辆车的车辆质量,
Figure BDA0003116346470000013
为第i辆车的风阻力带来的减速度,αr,i=gCr,icos(θi)为第i辆车的滚动阻力带来的减速度,αg,i=gsin(θi)为第i辆车的坡度阻力带来的减速度,ρ为空气密度,Cd,i为第i辆车的风阻系数,Af,i为第i辆车的迎风面积,g为重力加速度,Cr,i为第i辆车的阻力摩擦系数,θi为第i辆车的道路坡度;其中Cd,i与车辆间距有关:
Figure BDA0003116346470000014
其中δ(i,i+1)为第i辆车与第i+1辆车之间的车辆间距即δ(i,i+1)=si+1-si,si+1为第i+1辆车的位移,
Figure BDA0003116346470000021
为没有前车影响时车辆的风阻系数,
Figure BDA0003116346470000022
Figure BDA0003116346470000023
为拟合参数;
S2、油耗与排放模型:
动能在时间域上的体现即为车辆的驱动功率,第i辆车的驱动功率公式为:
Pac,i=Ft,i·vi (3)
其中Pac,i为第i辆车的驱动功率;
因此发动机处于工作状态时,第i辆车的车辆的瞬时油耗建模为:
Fac,i=k1,iPac,i+k2,i (4)
其中,Fac,i为第i辆车的瞬时油耗,k1,i与k2,i为拟合参数;
温度与剩余氧气均与车辆的动能有关,因此第i辆车的排放模型建立为:
Figure BDA0003116346470000024
其中Nac,i为第i辆车的排放,c1,i,c2,i,c3,i为拟合参数;
S3、控制目标的确定
1)为简化设计,给出几点假设:
①合并为车队的过程中,仅考虑空气阻力对跟随车辆的影响,不考虑其他因素的影响;
②参与编队的车辆物理参数完全相同;
③车辆合并过程中不存在车辆超车行为;
2)基于公式(3)(4),给出参与N辆车准备进行编队行驶时的目标函数:
Figure BDA0003116346470000025
其中,J表示目标函数的值,因为Σ表示加和,所以
Figure BDA0003116346470000026
为从第i辆车的油耗一直累加到第N辆车的油耗,即参与编队的N辆车的总油耗项,
Figure BDA0003116346470000027
为从第i辆车的NOx排放一直累加到第N辆车的NOx排放,即参与编队的N辆车的NOx排放项,
Figure BDA0003116346470000028
为从第i辆车的舒适度一直累加到第N辆车的舒适度,即参与编队的N辆车的总舒适度项,ai为第i辆车的加速度,ωf待调的油耗权重参数,ωN为待调的排放权重参数,ωa为待调的舒适度的权重参数,t0为队列合并开始的时间,tf为队列合并结束的时间;
S4、问题约束的确定
对车辆队列形成的位置进行约束:
dm∈[dmin,dmax] (7)
其中dm这N辆车合并为队列时候的位置,dmin为这N辆合并为队列时候位置的最小限值,dmax为这N辆车合并为队列时候位置的最大限值,本专利中dmin=0,基于工程经验
Figure BDA0003116346470000031
df为基于智能网联信息获得的参与编队的车辆从当前时刻开始,将要共同行使路段的长度;车辆单独行驶时或者以队列行驶时,速度必须满足法规约束:
vi,vp∈[vr,min,vr,max] (8)
其中vp为车辆合并为队列后的队列车速,vr,min为法规限制的车速最小值,vr,max为法规限制的车速最大值;
每辆车与前车的车距均由前车以及本车的参数以及状态决定:
Figure BDA0003116346470000032
其中vi+1为第i+1辆车的速度,sbr,i为第i辆车的紧急刹车位移,sbr,i+1为第i辆车前方第i+1辆车的紧急刹车位移,Treact,i为第i辆车的驾驶员应急反应时间,abmax,i第i辆车的最大减速度,abmax,i+1第i+1辆车的最大减速度;
避免发生碰撞的最小安全车距为:
ssafe,(i,i+1)=max(viTreact,i,sbr,i-sbr,i+1) (10)
其中,ssafe,(i,i+1)为第i辆车与该车前方最近的第i+1辆车之间的安全车距;
为了最大程度减小风阻,因此需要车辆之间保持最小车距行驶:
δ(i,i+1)=ssafe,(i,i+1) (11)
根据以上公式,结合道路法规车速算出车辆安全车距范围为[ssafe,min ssafe,max],为了最大能力减小风阻油耗,选择约束下限作为最终队内车距ssafe,(i,i+1)=ssafe,min
给出车辆的加减速度约束:
amin,i≤ai≤amax,i (12)
其中,ai为第i辆车的加速度
Figure BDA0003116346470000033
amin,i为第i辆车的最小减速度,由制动器性能直接给出,amax,i为第i辆车的最大加速度,基于车辆发动机性能给出:
Figure BDA0003116346470000041
其中,ηt,i为第i辆车的传动系统的传动效率,rt,i为第i辆车的轮胎半径,Te,max,i为第i辆车的发动机能输出的最大有效力矩;
控制器量的确定
选择车辆的控制输入力为控制输入量;
优化问题求解终端时间自由的优化问题,即公式(6)中的tf∈[tmin,tmax],结合前面内容给出队列合并的时间范围:
Figure BDA0003116346470000042
Figure BDA0003116346470000043
其中tmin为队列合并的最终时间的最小值,tmax为队列合并的最终时间的最大值;
将可能的队列合并的最终时间区域网格化,每个时间网格节点即表示一个可能的队列合并结束时间:
Figure BDA0003116346470000044
其中,j表示时间区间[tmin,tmax]网格化的第j个时间节点,Δ表示将时间区域网格化的间隔,这样在每个时间节点tf,j上,目标函数(6)中的tf都取值为tf,j,这样目标函数(6)就转化为一个终端时间固定的目标函数,从而可以进行计算,计算这个终端时间固定的目标函数,获得每个时间节点tf,j处的目标函数的值J;并根据对比选出其中最小的J及其对应的车辆的状态以及控制输入力,即为通过瞬时速度优化获得的这些车辆编队过程中的状态以及这些车辆编队过程中所需要的控制输入力。
本发明首先建立起基于参与队列合并时单车车速的排放与油耗模型,其次根据智能交通信息中获得的本车与前车的相对车距,相对车速与共同行使的路程等信息以及车辆底层发动机系统的性能信息给出车辆队列合并时的加速度、速度与安全距离约束,对参与队列合并的车辆的车速进行了规划,获得了可以使队列整体油耗与排放综合最优的每辆车的车速以及对应该车速的控制输入(车辆输入控制力)。最后为了验证本方法的有效性,在Matlab中进行了仿真验证。
附图说明
图1是风阻系数模型验模曲线图;
图2是油耗模型相关性分析图;
图3是NOx模型相关性分析图;
图4是参与队列合并的车辆轨迹、速度、风阻系数、车间距曲线图;
图5是单独车辆油耗对比柱状图;
图6是单独车辆NOx排放对比柱状图;
图7是车队油耗对比柱状图;
图8是车队NOx排放对比柱状图。
具体实施方式
本发明使用智能交通信息,对参与队列合并的商用车车速与总驱动力进行规划,获得队列整体排放与油耗综合最优的队列合并行为。本专利首先建立起基于参与队列合并时单车车速的排放与油耗模型,其次根据智能交通信息中获得的本车与前车的相对车距,相对车速与共同行使的路程等信息以及车辆底层发动机系统的性能信息给出车辆队列合并时的加速度、速度与安全距离约束,对参与队列合并的车辆的车速进行了规划,获得了可以使队列整体油耗与排放综合最优的每辆车的车速以及对应该车速的控制输入(车辆输入控制力)。最后为了验证本方法的有效性,在Matlab中进行了仿真验证。
研究表明重型货运车辆以很小间距列队行驶时能有效减少跟随车的空气阻力,从而降低风阻油耗。智能交通系统的发展使得车辆在编队行驶时能够获得道路信息,同路车辆信息等智能信息,从而优化车辆编队行为,获得车辆编队时整体油耗与排放综合最优的效果。本专利通过智能网联信息,获得同路车辆目的地、车速、加速度等信息,在车辆编队的瞬态行为优化问题中,基于降低风阻以及安全需求给出约束条件,构建综合油耗与排放的优化问题,优化车辆瞬时速度的编队行为。最后通过仿真实验,与现有的优化车辆平均速度的编队行为对比,验证了本发明提出方法的有效性。
本发明的具体方法如下:
1.车速模型搭建
以柴油商用车为研究对象,当道路上有N辆车准备进行编队行驶时,选取行驶位移和车速为状态变量,构建参与编队的第i辆车的纵向动力学系统方程如下
Figure BDA0003116346470000051
Figure BDA0003116346470000052
其中i代表参与编队的第i辆车,各类符号分别为第i辆车的:si为第i辆车的车辆位移,vi为第i辆车的车速,Ft,i为第i辆车的车辆控制输入力,Mi为第i辆车的车辆质量,
Figure BDA0003116346470000053
为第i辆车的风阻力带来的减速度,αr,i=gCr,icos(θi)为第i辆车的滚动阻力带来的减速度,αg,i=gsin(θi)为第i辆车的坡度阻力带来的减速度。ρ为空气密度,Cd,i为第i辆车的风阻系数,Af,i为第i辆车的迎风面积,g为重力加速度,Cr,i为第i辆车的阻力摩擦系数,θi为第i辆车的道路坡度。
其中Cd,i与车辆间距有关:
Figure BDA0003116346470000054
其中δ(i,i+1)为第i辆车与第i+1辆车之间的车辆间距即δ(i,i+1)=si+1-si,si+1为第i+1辆车的位移,
Figure BDA0003116346470000055
为没有前车影响时车辆的风阻系数,
Figure BDA0003116346470000056
Figure BDA0003116346470000057
为拟合参数。
2.油耗与排放模型:
根据热能与机械原理,发动机的燃油燃烧产生的能量通过传动系统的传递驱动车辆运行,因此车辆的动能与发动机的油耗存在正比例关系,而动能在时间域上的体现即为车辆的驱动功率,第i辆车的驱动功率计算公式为:
Pac,i=Ft,i·vi (3)
其中Pac,i为第i辆车的驱动功率。
因此发动机处于工作状态时,第i辆车的车辆的瞬时油耗可以建模为:
Fac,i=k1,iPac,i+k2,i (4)
其中,Fac,i为第i辆车的瞬时油耗,k1,i与k2,i为拟合参数。
发动机只要运行就会产生排放,NOx在发动机内部处于高温富养的环境下会加速生成,温度与剩余氧气均与车辆的动能有关,因此第i辆车的排放模型可以建立为:
Figure BDA0003116346470000061
其中Nac,i为第i辆车的排放,c1,i,c2,i,c3,i为拟合参数。
3.控制问题描述
3.1控制目标的确定
在N辆车进行编队时,前方车辆可以通过适当的减速,同时后方车辆适当加速来加快车队的合并过程,这在样车队形成后,队内车辆将能够以更小的风阻行驶更长的距离,以此降低了更多的风阻油耗。因此这种方式比前方车速度不变或者加速,仅靠后方车辆加速追击形成车队具有更大的节油潜能。在这种形成队列方式下,前方车辆减速行驶时会带来动能转化为热能消散的能量损失,后方车辆加速行驶会产生额外的油耗与排放,这是因为后车加速行驶,发动机处于高温富养的工作环境下,所以排放的增加不可避免,因此本专利将排放优化放在控制车辆合并为队列过程中,控制参与编队的车辆以油耗和排放综合最优为目标形成车队。为简化设计,给出几点假设:
①合并为车队的过程中,仅考虑空气阻力对跟随车辆的影响,不考虑其他因素的影响;
②参与编队的车辆物理参数完全相同;
③车辆合并过程中不存在车辆超车行为。
基于公式(3)(4),给出参与N辆车准备进行编队行驶时的目标函数:
Figure BDA0003116346470000062
其中,J表示目标函数的值。因为Σ表示加和,所以
Figure BDA0003116346470000063
为从第i辆车的油耗一直累加到第N辆车的油耗,即参与编队的N辆车的总油耗项。
Figure BDA0003116346470000064
为从第i辆车的NOx排放一直累加到第N辆车的NOx排放,即参与编队的N辆车的NOx排放项。
Figure BDA0003116346470000065
为从第i辆车的舒适度一直累加到第N辆车的舒适度,即参与编队的N辆车的总舒适度项。ai为第i辆车的加速度,ωf待调的油耗权重参数,ωN为待调的排放权重参数,ωa为待调的舒适度的权重参数,t0为队列合并开始的时间,tf为队列合并结束的时间。
3.2问题约束的确定
参与编队的车辆确定后,这些车辆只有在一定范围内形成队列才有节油的潜力,否则车辆进行编队过程中增加的油耗会超过以队列行驶节省的风阻油耗,因此需要对车辆队列形成的位置进行约束:
dm∈[dmin,dmax] (7)
其中dm这N辆车合并为队列时候的位置,dmin为这N辆合并为队列时候位置的最小限值,dmax为这N辆车合并为队列时候位置的最大限值,本专利中dmin=0,基于工程经验
Figure BDA0003116346470000071
df为基于智能网联信息获得的参与编队的车辆从当前时刻开始,将要共同行使路段的长度。
车辆单独行驶时或者以队列行驶时,速度必须满足法规约束:
vi,vp∈[vr,min,vr,max] (8)
其中vp为车辆合并为队列后的队列车速,vr,min为法规限制的车速最小值,vr,max为法规限制的车速最大值。
由于队列内部的车间距并不相等,每辆车与前车的车距均由前车以及本车的参数以及状态决定:
Figure BDA0003116346470000072
其中vi+1为第i+1辆车的速度,sbr,i为第i辆车的紧急刹车位移,sbr,i+1为第i辆车前方第i+1辆车的紧急刹车位移,Treact,i为第i辆车的驾驶员应急反应时间,abmax,i第i辆车的最大减速度,abmax,i+1第i+1辆车的最大减速度。
所以避免发生碰撞的最小安全车距为:
ssafe,(i,i+1)=max(viTreact,i,sbr,i-sbr,i+1) (10)
其中,ssafe,(i,i+1)为第i辆车与该车前方最近的第i+1辆车之间的安全车距。
为了最大程度减小风阻,因此需要车辆之间保持最小车距行驶:
δ(i,i+1)=ssafe,(i,i+1) (11)
根据以上公式,结合道路法规车速算出车辆安全车距范围为[ssafe,min ssafe,max],为了最大能力减小风阻油耗,本专利选择约束下限作为最终队内车距ssafe,(i,i+1)=ssafe,min
由于车辆加速度取决于底层发动机的动力性能,车辆的减速度取决于制动器的性能,所以根据底层发动机动力性能以及制动器性能,给出车辆的加减速度约束:
amin,i≤ai≤amax,i (12)
其中,ai为第i辆车的加速度
Figure BDA0003116346470000081
amin,i为第i辆车的最小减速度,由制动器性能直接给出,amax,i为第i辆车的最大加速度。
基于车辆发动机性能给出:
Figure BDA0003116346470000082
其中,ηt,i为第i辆车的传动系统的传动效率,rt,i为第i辆车的轮胎半径,Te,max,i为第i辆车的发动机能输出的最大有效力矩。
控制器量的确定由于车辆在路上形式过程中,直接表现的状态为位移,速度,加速度。而位移为速度的积分,加速度为速度的微分,所以选取车辆的状态为速度,根据公式(1b)可知,对车速的影响为车辆的控制输入力,这个控制输入力为正时表示发动机工作驱动车辆前进,控制输入力为负时表示发动机不工作,制动器工作,所以选择车辆的控制输入力为控制输入量。
优化问题求解由于队列合并时,只对最终队列形成位置有约束要求,对队列合并的最终时间并没有要求,所以这个终端时间自由的优化问题,即公式(6)中的tf∈[tmin,tmax]。
结合前面内容给出队列合并的时间范围:
Figure BDA0003116346470000083
Figure BDA0003116346470000084
其中tmin为队列合并的最终时间的最小值,tmax为队列合并的最终时间的最大值。
但是并不知道tf的具体值,因此目标函数(6)无法计算。为了计算这个终端时间自由的目标函数,将可能的队列合并的最终时间区域网格化,这样每个时间网格节点即表示一个可能的队列合并结束时间:
Figure BDA0003116346470000085
其中,j表示时间区间[tmin,tmax]网格化的第j个时间节点,Δ表示将时间区域网格化的间隔,这样在每个时间节点tf,j上,目标函数(6)中的tf都取值为tf,j,这样目标函数(6)就转化为一个终端时间固定的目标函数,从而可以进行计算。计算这个终端时间固定的目标函数,我们可以获得每个时间节点tf,j处的目标函数的值J。并根据对比选出其中最小的J及其对应的车辆的状态以及控制输入力,即为通过瞬时速度优化获得的这些车辆编队过程中的状态以及这些车辆编队过程中所需要的控制输入力。
4.仿真验证及分析
从GT-Power中采集初来车辆实际运行时的数据,根据第2章建模方法,搭建了车辆排放面向控制的模型,根据GT的数据辨识得到模型参数,获得最终的在智能网联环境下的车辆油耗与排放控制模型,在相同的车辆工况下,对GT数据与建立的模型的输出数据进行相关分析,如图1、2、3所示,可以看到风阻系数模型、油耗模型与NOx排放模型表征相关性的R2均大于0.7,表明模型可以较好地描述线管的车辆动态,可以作为车辆面向控制的模型进行下一步的优化规划。
求解得到使得油耗与排放综合最优的车辆速度,与对应的车辆输入控制力,其中输入控制力为负值时,表明此时输入控制力为制动力,车辆的制动器工作,此时制动力可以直接提供给车辆模型的制动器,当求解得到的输入控制力为正值时,表明此时输入控制力为驱动力,则通过传动系统的速比转化为发动机控制器的目标有效力矩,由发动机控制器控制发动机提供此需求力矩,并经过传动系统的转化驱动车辆前进。
最后的求解结果与现有的基于平均速度优化的队列合并方法进行了对比,以验证本专利提出方法的有效性。实验图如附图中所示,实验数据与简单分析如附表中所示。
从图4的前两幅图可以看出,本专利提出的方法可以有效地将车辆编队,并且从图4的后两幅图看出,车辆编队行驶后风阻确实会降低,因此风阻带来的油耗也会降低。结合图4、5、7以及表3可以看出,相比于平均速度优化的车辆队列合并方法,本专利提出的瞬时车速优化的车辆队列合并方法,可以动态调整参与编队的每辆车的速度,前面车辆会减速等待后方车辆,后方车辆则会加速追及前方车辆,所以虽然参与编队的车辆的油耗有的高有的低,但是会达到队列整体油耗最低的效果。同样结合图4、6、8以及表4可以看出,相比于平均速度优化的车辆队列合并方法,本专利提出的瞬时车速优化的车辆队列合并方法,可以动态调整参与编队的每辆车的速度,前面车辆会减速等待后方车辆,后方车辆则会加速追及前方车辆,所以虽然参与编队的车辆的排放有的高有的低,但是会达到队列整体油耗最低的效果。单车油耗与排放变化的百分比与车队整体油耗与排放变化的百分比均整理在表3、4中。需要特别解释的是,车辆2与车辆3之所以会油耗与排放最低,是因为在瞬时速度优化时,车辆2与车辆3处于减速状态,其发动机则维持在怠速转速,所以车辆2、3的油耗与排放最低。
表1模型参数
根据假设②,第i辆车的数值即为i=1,2,…,N这些参与编队的每辆车的数值。
本申请中选取5辆车作为研究对象,所以N=5
Figure BDA0003116346470000091
Figure BDA0003116346470000101
表2参与编队车辆初始状态
车辆序号 初始速度(m/s) 初始位置(m)
车辆1(i=1) 28 250
车辆2(i=2) 26 210
车辆3(i=3) 24 190
车辆4(i=4) 22 160
车辆5(i=5) 20 130
表3油耗实验结果
表中“数据分析”中百分比均为瞬时速度优化相对于平均速度优化的计算结果,即(瞬时速度优化油耗结果-平均速度优化油耗结果)/平均速度优化油耗结果
Figure BDA0003116346470000102
Figure BDA0003116346470000111
表4NOx排放实验结果
表中“数据分析”中百分比均为瞬时速度优化相对于平均速度优化的计算结果,即(瞬时速度优化NOx结果-平均速度优化NOx结果)/平均速度优化NOx结果
Figure BDA0003116346470000112

Claims (1)

1.一种油耗与排放一体化优化的商用车队列合并方法,其特征在于:
S1、车速模型搭建
以柴油商用车为研究对象,当道路上有N辆车准备进行编队行驶时,选取行驶位移和车速为状态变量,构建参与编队的第i辆车的纵向动力学系统方程如下
Figure FDA0003116346460000011
Figure FDA0003116346460000012
其中i代表参与编队的第i辆车,各类符号分别为第i辆车的:si为第i辆车的车辆位移,vi为第i辆车的车速,Ft,i为第i辆车的车辆控制输入力,Mi为第i辆车的车辆质量,
Figure FDA0003116346460000013
为第i辆车的风阻力带来的减速度,αr,i=gCr,icos(θi)为第i辆车的滚动阻力带来的减速度,αg,i=gsin(θi)为第i辆车的坡度阻力带来的减速度,ρ为空气密度,Cd,i为第i辆车的风阻系数,Af,i为第i辆车的迎风面积,g为重力加速度,Cr,i为第i辆车的阻力摩擦系数,θi为第i辆车的道路坡度;其中Cd,i与车辆间距有关:
Figure FDA0003116346460000014
其中δ(i,i+1)为第i辆车与第i+1辆车之间的车辆间距即δ(i,i+1)=si+1-si,si+1为第i+1辆车的位移,
Figure FDA0003116346460000015
为没有前车影响时车辆的风阻系数,
Figure FDA0003116346460000016
Figure FDA0003116346460000017
为拟合参数;
S2、油耗与排放模型:
动能在时间域上的体现即为车辆的驱动功率,第i辆车的驱动功率公式为:
Pac,i=Ft,i·vi (3)
其中Pac,i为第i辆车的驱动功率;
因此发动机处于工作状态时,第i辆车的车辆的瞬时油耗建模为:
Fac,i=k1,iPac,i+k2,i (4)
其中,Fac,i为第i辆车的瞬时油耗,k1,i与k2,i为拟合参数;
温度与剩余氧气均与车辆的动能有关,因此第i辆车的排放模型建立为:
Figure FDA0003116346460000018
其中Nac,i为第i辆车的排放,c1,i,c2,i,c3,i为拟合参数;
S3、控制目标的确定
1)为简化设计,给出几点假设:
①合并为车队的过程中,仅考虑空气阻力对跟随车辆的影响,不考虑其他因素的影响;
②参与编队的车辆物理参数完全相同;
③车辆合并过程中不存在车辆超车行为;
2)基于公式(3)(4),给出参与N辆车准备进行编队行驶时的目标函数:
Figure FDA0003116346460000021
其中,J表示目标函数的值,因为Σ表示加和,所以
Figure FDA0003116346460000022
为从第i辆车的油耗一直累加到第N辆车的油耗,即参与编队的N辆车的总油耗项,
Figure FDA0003116346460000023
为从第i辆车的NOx排放一直累加到第N辆车的NOx排放,即参与编队的N辆车的NOx排放项,
Figure FDA0003116346460000024
为从第i辆车的舒适度一直累加到第N辆车的舒适度,即参与编队的N辆车的总舒适度项,ai为第i辆车的加速度,ωf待调的油耗权重参数,ωN为待调的排放权重参数,ωa为待调的舒适度的权重参数,t0为队列合并开始的时间,tf为队列合并结束的时间;
S4、问题约束的确定
对车辆队列形成的位置进行约束:
dm∈[dmin,dmax] (7)
其中dm这N辆车合并为队列时候的位置,dmin为这N辆合并为队列时候位置的最小限值,dmax为这N辆车合并为队列时候位置的最大限值,本专利中dmin=0,基于工程经验
Figure FDA0003116346460000025
df为基于智能网联信息获得的参与编队的车辆从当前时刻开始,将要共同行使路段的长度;
车辆单独行驶时或者以队列行驶时,速度必须满足法规约束:
vi,vp∈[vr,min,vr,max] (8)
其中vp为车辆合并为队列后的队列车速,vr,min为法规限制的车速最小值,vr,max为法规限制的车速最大值;
每辆车与前车的车距均由前车以及本车的参数以及状态决定:
Figure FDA0003116346460000026
其中vi+1为第i+1辆车的速度,sbr,i为第i辆车的紧急刹车位移,sbr,i+1为第i辆车前方第i+1辆车的紧急刹车位移,Treact,i为第i辆车的驾驶员应急反应时间,abmax,i第i辆车的最大减速度,abmax,i+1第i+1辆车的最大减速度;
避免发生碰撞的最小安全车距为:
ssafe,(i,i+1)=max(viTreact,i,sbr,i-sbr,i+1) (10)
其中,ssafe,(i,i+1)为第i辆车与该车前方最近的第i+1辆车之间的安全车距;
为了最大程度减小风阻,因此需要车辆之间保持最小车距行驶:
δ(i,i+1)=ssafe,(i,i+1) (11)
根据以上公式,结合道路法规车速算出车辆安全车距范围为[ssafe,min ssafe,max],为了最大能力减小风阻油耗,选择约束下限作为最终队内车距ssafe,(i,i+1)=ssafe,min
给出车辆的加减速度约束:
amin,i≤ai≤amax,i (12)
其中,ai为第i辆车的加速度
Figure FDA0003116346460000031
amin,i为第i辆车的最小减速度,由制动器性能直接给出,amax,i为第i辆车的最大加速度,基于车辆发动机性能给出:
Figure FDA0003116346460000032
其中,ηt,i为第i辆车的传动系统的传动效率,rt,i为第i辆车的轮胎半径,Te,max,i为第i辆车的发动机能输出的最大有效力矩;
控制器量的确定
选择车辆的控制输入力为控制输入量;
优化问题求解
终端时间自由的优化问题,即公式(6)中的tf∈[tmin,tmax],结合前面内容给出队列合并的时间范围:
Figure FDA0003116346460000033
Figure FDA0003116346460000034
其中tmin为队列合并的最终时间的最小值,tmax为队列合并的最终时间的最大值;
将可能的队列合并的最终时间区域网格化,每个时间网格节点即表示一个可能的队列合并结束时间:
Figure FDA0003116346460000041
其中,j表示时间区间[tmin,tmax]网格化的第j个时间节点,Δ表示将时间区域网格化的间隔,这样在每个时间节点tf,j上,目标函数(6)中的tf都取值为tf,j,这样目标函数(6)就转化为一个终端时间固定的目标函数,从而可以进行计算,计算这个终端时间固定的目标函数,获得每个时间节点tf,j处的目标函数的值J;并根据对比选出其中最小的J及其对应的车辆的状态以及控制输入力,即为通过瞬时速度优化获得的这些车辆编队过程中的状态以及这些车辆编队过程中所需要的控制输入力。
CN202110664644.0A 2021-06-16 2021-06-16 油耗与排放一体化优化的商用车队列合并方法 Active CN113408214B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110664644.0A CN113408214B (zh) 2021-06-16 2021-06-16 油耗与排放一体化优化的商用车队列合并方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110664644.0A CN113408214B (zh) 2021-06-16 2021-06-16 油耗与排放一体化优化的商用车队列合并方法

Publications (2)

Publication Number Publication Date
CN113408214A true CN113408214A (zh) 2021-09-17
CN113408214B CN113408214B (zh) 2022-06-10

Family

ID=77684176

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110664644.0A Active CN113408214B (zh) 2021-06-16 2021-06-16 油耗与排放一体化优化的商用车队列合并方法

Country Status (1)

Country Link
CN (1) CN113408214B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114282804A (zh) * 2021-12-23 2022-04-05 同济大学 一种基于燃油节约率的卡车编队方案评估方法及装置
CN116820080A (zh) * 2023-01-05 2023-09-29 北京交通大学 高速公路自动驾驶货车编队运行调度方法、系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20190171227A1 (en) * 2017-12-01 2019-06-06 Cummins Inc. Optimization of mission efficiency through platoon opportunity assessment
US20190228649A1 (en) * 2018-01-23 2019-07-25 International Business Machines Corporation Detection of vehicle queueing events on a road
CN110329258A (zh) * 2019-07-23 2019-10-15 吉林大学 智能驾驶汽车节能减排协调控制方法
CN111422192A (zh) * 2020-01-21 2020-07-17 吉林大学 智能网联商用车跟车行为下的排放与油耗协调控制系统
CN111439260A (zh) * 2020-04-27 2020-07-24 吉林大学 面向个性化需求的网联商用柴油车巡航行驶优化控制系统

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20190171227A1 (en) * 2017-12-01 2019-06-06 Cummins Inc. Optimization of mission efficiency through platoon opportunity assessment
US20190228649A1 (en) * 2018-01-23 2019-07-25 International Business Machines Corporation Detection of vehicle queueing events on a road
CN110329258A (zh) * 2019-07-23 2019-10-15 吉林大学 智能驾驶汽车节能减排协调控制方法
CN111422192A (zh) * 2020-01-21 2020-07-17 吉林大学 智能网联商用车跟车行为下的排放与油耗协调控制系统
CN111439260A (zh) * 2020-04-27 2020-07-24 吉林大学 面向个性化需求的网联商用柴油车巡航行驶优化控制系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
刘迪: "高速公路多车协同驾驶控制策略研究", 《CNKI优秀硕士论文全文库》 *
陈虹 郭露露 宫洵 高炳钊 张琳: "智能时代的汽车控制", 《自动化学报》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114282804A (zh) * 2021-12-23 2022-04-05 同济大学 一种基于燃油节约率的卡车编队方案评估方法及装置
CN114282804B (zh) * 2021-12-23 2022-11-15 同济大学 一种基于燃油节约率的卡车编队方案评估方法及装置
CN116820080A (zh) * 2023-01-05 2023-09-29 北京交通大学 高速公路自动驾驶货车编队运行调度方法、系统
CN116820080B (zh) * 2023-01-05 2024-04-30 北京交通大学 高速公路自动驾驶货车编队运行调度方法、系统

Also Published As

Publication number Publication date
CN113408214B (zh) 2022-06-10

Similar Documents

Publication Publication Date Title
CN108475465B (zh) 用于移动、特别是控制或者调整车队的方法
CN111439260B (zh) 面向个性化需求的网联商用柴油车巡航行驶优化控制系统
Liang et al. Heavy-duty vehicle platoon formation for fuel efficiency
CN113408214B (zh) 油耗与排放一体化优化的商用车队列合并方法
CN112437412B (zh) 一种基于车路协同的混行驾驶车辆编队控制方法
Huang et al. Speed trajectory planning at signalized intersections using sequential convex optimization
CN107187447B (zh) 一种基于车联网的车辆自适应巡航控制系统及其控制方法
EP2867560B1 (en) Transmission control system
CN110893853B (zh) 一种基于前方坡度信息的车辆控制方法以及系统
CN107628029A (zh) 一种网联汽车队列的节能型稳定性运动控制方法
CN111422192B (zh) 智能网联商用车跟车行为下的排放与油耗协调控制方法
Kotiev et al. Energy efficient motion control of the electric bus on route
CN110188401B (zh) 一种基于改进pso的有轨电车运行能耗优化方法
CN112731806B (zh) 智能网联汽车随机模型预测控制实时性优化方法
CN113821966A (zh) 高速磁浮列车运行节能优化方法、系统及存储介质
Kaku et al. Model predictive control for ecological vehicle synchronized driving considering varying aerodynamic drag and road shape information
D'Amato et al. Development of a cruise controller based on current road load information with integrated control of variable velocity set-point and gear shifting
Bentaleb et al. Energy-optimal control for eco-driving on curved roads
Held et al. Optimal freewheeling control of a heavy-duty vehicle using mixed integer quadratic programming
Ganji et al. Look-ahead intelligent energy management of a parallel hybrid electric vehicle
CN113110022A (zh) 一种基于非线性pid的多车列队纵向跟随控制方法及装置
Gong et al. Benchmark study on real-time energy optimization of hevs under connected environment
CN116118705B (zh) 一种跟驰场景下插电式混合动力公交车能量管理控制方法
Pan et al. Energy-optimized adaptive cruise control strategy design at intersection for electric vehicles based on speed planning
Junkai et al. Speed planning and energy optimal control of hybrid electric vehicles based on internet of vehicles

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