2021年第2期
【摘要】为完成中国重型商用车辆行驶工况(CHTC )的开发,采集并处理多个城市的车辆实际运行数据,建立了短片段
数据库,分析对比各类型车辆数据库特征,确定工况的体系构成和开发路线,通过聚类和加权确定各车型工况曲线的速度区间结构,通过以速度-加速度联立分布为基础的卡方检验选择最优片段进行组合构建曲线,最终完成CHTC 的开发。对比实测工况数据及其他典型工况,证明了CHTC 的可靠性。
主题词:重型商用车车辆运行特征速度-加速度分布工况开发
中图分类号:U467.12
文献标识码:A DOI:10.19620/jki.1000-3703.20191355
Development of the China Heavy-Duty Commercial Vehicle Test Cycle
Hu Xi,Liu Yu,Li Mengliang
(China Automotive Technology &Research Center Co.,Ltd.,Tianjin 300300)采埃孚
【Abstract 】To complete development of China Heavy Truck Cycle (CHTC),vehicle real running data of several cities are collected and processed,short fragment database is established,characteristics of vehicle database of different varieties of vehicles are analyzed and compared,thus the driving cycle system construction and development route are defined.Speed range structure of cycle curves of different vehicle models is defined by clustering and weighting.The optimal fragment is selected by chi-square test based on speed-acceleration distribution,thus the CHTC cycle is constructed.Measured cycle data and other typical driving cycle are compared,which proves reliability of CHTC.
Key words:Heavy-duty commercial vehicle,Vehicle driving characteristic,Speed-acceleration
distribution,Driving cycle development
胡熙
刘昱
李孟良
(中国汽车技术研究中心有限公司,天津300300)
*基金项目:国家重点研发计划项目(2018YFB0106402)。
中国重型商用车行驶工况开发*
汽车技术·Automobile Technology
1前言
汽车行驶工况是汽车行业的一项重要基础标准,重型商用车工况是该类型车性能优化时的主要基准,也是进行排放和油耗认证的基础。现阶段我国进行重型车油耗认证时采用的是以世界重型商用车辆瞬态循环(World Transient Vehicle Cycle ,WTVC )工况为基础调整生成的C-WTVC 工况[1]。
近年来,国内对重型车行驶工况开发进行了一系列研究。武汉理工大学的吴其伟等用多元统计理论构建了武汉公交车行驶工况[2];中国环境科学研究院的王军方等根据北京市乘用车、公交车和重型载货汽车的运行数据建立了北京市典型行驶工况[3];秦大同等提出一种基于C 均值聚类算法的城市循环工况构建方法[4]。
研究表明,上述重型商用车认证和研究存在一定的不足[5-6]:C-WTVC 工况更多地采用了欧美国家的数据,其他工况研究采用的数据局限在特定城市或区域,均不能很好地代表中国综合实际道路行驶状况;重型商用车类型繁多,体系复杂,不同类型车辆的运行范围、使用方式和车辆性能等差异显著,
现有工况无法充分反映这些差异。
工信部于2015年委托中国汽车技术研究中心牵头实施“中国新能源汽车产品检测工况研究和开发”(简称“中国工况”)项目,研究我国车辆实际行驶工况特征并开发各类车辆的全国标准行驶循环工况(China Automotive Test Cycle ,CATC )。中国重型商用车辆行驶工况(China Heavy-Duty Commercial Vehicle Test Cycle ,
CHTC )是该项目的重要组成部分。
本文通过短片段处理分析对比不同类型重型商用
【引用格式】胡熙,刘昱,李孟良.中国重型商用车行驶工况开发[J].汽车技术,2021(2):51-57.
HU X,LIU Y,LI M L.Development of the China Heavy-Duty Commercial Vehicle Test Cycle[J].Automobile Technology,2021(2):51-57.
-
-51
车的运行数据库特征,进而采用聚类、加权及卡方检验的方法构建工况,最终完成包含多条曲线的CHTC工况体系的开发,令其能够全面覆盖、真实反映各类型重型车行驶情景。
2工况开发
以采集到的多个城市、多种类型的重型车队运行数据为基础,开发重型商用车辆行驶工况,总体流程如图1所示。
图1重型商用车辆行驶工况开发流程
2.1数据采集和处理
中国工况项目组建的车队包含1216辆重型商用车辆,在我国41个代表城市采用自主驾驶法采集数据,即不规划道路,车辆自由行驶,以保证数据真实性。
采集的数据包括车辆定位、动力总成、环境信息等多种类参数,采样频率为1Hz。数据在城市、气候、道路和车型等方面具有广泛代表性和覆盖度。
经过1年多的稳定行驶,采集了累积里程接近23×106km的各类型重型商用车辆运行数据,车型构成和相应里程如表1所示。
表1重型商用车车队构成和相应里程
CHTC以短片段为基础构建:将车辆从一次起步到下一次停车的运动定义为运动片段;一次停车到下一次起步的运动定义为怠速片段。一个短片段包含一个运动片段和与之相邻的怠速片段,如图2所示。
在数据采集过程中,会出现一些不合理或无效数据,根据车辆具体运行情况制定包含运行片段持续时长、怠速时长、速度范围、加速度范围和数据缺失率在内的筛选规则,选出有效短片段建立片段库,用于后续工况开发。
图2短片段示意
对于每一个短片段,计算其运行时间、加速比例、减速比例、匀速比例、怠速比例、运行距离、最大速度、平均速度、运行平均速度、最大加速度、最大减速度、平均正加速度、平均负加速度、速度标准差、加速度标准差共15个特征参数,并归一化作为片段特征。不同特征的短片段反映车辆在不同道路和交通状况下的行驶情景。
2.2工况体系构成规划
重型商用车辆的运行方式受用途影响较大,首先按照用途将其初步划分为客车和载货汽车。分别提取
这2类车型各细分车辆类型的短片段库,计算并对比各片段库的总体特征参数,判断是否为特定车型单独构建工况。
2.2.1客车运行特征分析和对比
将客车划分为城市客车和普通客车。城市客车即公交车,工况构建采用的车队数据覆盖不同等级的22个城市。普通客车以在城市间运行的长途客车和在城郊运行的班车为主,数据覆盖14个城市。
计算两者短片段库的关键特征参数进行对比分析,结果如图3所示。
(a)速度特征
(b )时间比例
图3客车片段库特征参数对比
特定类型重型车的行驶数据采集短行程切割、筛选
特征值求取
主成分分析、聚类分析
各类特征计算、分析、比较
工况曲线各速度
区间长度确定n类:运行/怠速
片段库n类:速度-加速
度分布
确定片段数量及长度片段组合
抽取
卡方检验
工况构建
包头市违章查询
CHTC
重型车载货特征调查出行特征调查
参数
数量/辆里程/×106km 城市客车
568
14.61
普通客车
216
3.40
货车
265
2.38
自卸汽车
86
0.43
半挂牵引车
81
1.79
时间
怠速片段1运动片段1怠速片段2运动片段2加速
减速
普通客车城市客车
40
30
20
10
/
k
m
宝马4s
·
h
-
1
/
m
·
s
-
2
0.6
0.4
0.2
-0.2
-0.4
-
0.6
-0.8
平均速度正加速度负加速度
020*********
时间比例/%
城市客车
普通客车
加速减速匀速怠速
17.80%
20.86%
-
-
52
汽车技术
2021年
第2期
从图3中可以看出,城市客车平均速度约为16.8km/h ,显著低于普通客车,怠速比例为34%,显著偏高。二者差异的主要原因在于运行范围不同:城市客车主要在交通拥堵的市区行驶,普通客车运行范围相对广泛。鉴于2类车运行特征之间的较大差异,有必要分别构建行驶工况。
此外,比较城市客车车队在各城市的速度分布情况,如图4所示。
(a )城市客车的速度分布(第1组)
(b )城市客车的速度分布(第2组)
图4
城市客车的速度分布
从图4中可以看出,各城市间速度分布存在较大差异,反映了不同城市的城市客车具有独特的行驶工况特征,构建城市客车工况时需考虑上述差异,获取各城市权重,对采集数据进行加权处理。
各城市普通客车的运行特征差异不大,构建工况时不需要进行城市加权。
2.2.2
载货汽车运行特征分析和对比
各城市载货汽车的运行特征差异不大,构建工况时不需要进行城市加权。
不同类型载货汽车的运行范围差异较大:货车运行区域较广泛,自卸汽车主要在市区和城郊行驶,半挂牵引车在市内运行受限,主要在城郊和高速公路上行驶。与客车类似,计算几种类型载货汽车的短片段库综合特
征进行对比,如图5所示。
(a )速度特征
(b )时间比例
图5
载货汽车片段库特征参数对比
从图5中可以看出,3种车型的平均速度差异显著,半挂牵引车怠速比例显著偏低,匀速比例显著偏高,有
必要分别构建工况曲线。
根据相关法规,总质量(Gross Vehicle Weight ,GVW )较小的货车可以较自由地进入市区,运行区域与其他货车存在较大差异,需要判断是否分别构建工况循环。
分别对轻型货车(总质量≤5500kg )与其他质量等级的货车(称为重型货车)片段库进行简单的速度区间划分,对比各区间的片段库特征参数,如图6所示。
(a )运行平均速度和运动片段持续时间
(b )时间比例
图6
各质量级别货车分区间片段库特征参数
0~<
1
010~<2020~<3030~<4040~<5050~<6060~<7070~<8080~<90
车速区间/km·h -1
35302520151050
分布频率/%深圳
福州芜湖
天津合肥
呼和浩特
昆明株洲海口
青岛镇江
0~
<1
10
~<
20
20
~<
30
30
~<
2013捷达
40
40
~<
50
50
~<
60
60
~<
70
70
~
<
80
80
~<
90
车速区间/km·h -1
35302520151050
分布频率/%北京成都太原
广州武汉南宁
上海西安兰州
重庆
厦门
货车
自卸汽车
50403020100
平均速度/k m ·h -1
平均正、负加速度/m ·s -2
0.6
0.40.20-0.2-0.4-0.6
平均速度
正加速度
负加速度
半挂牵引车
20
4060
80
100时间比例/%
半挂牵引车自卸汽车
货车
加速
减速
匀速
怠速
6.90%18.13%14.25%
低速
中速
7060
50403020100
运行平均速度/k m ·h -1
运动片段持续时间度/s
500
4003002001000
轻型货车运行平均速度轻型货车运行时间
高速
重型货车运行平均速度重型货车运行时间
20
406080100
时间比例/%
轻型货车重型货车轻型货车重型货车轻型货车重型货车
加速
减速
匀速
怠速
5.10%低速中速高速
5.61%11.33%8.13%
29.37%
40.06%
-
-53
从图6中可以看出:轻型货车在低速区间速度偏高,怠速比例偏低,说明它在该区间行驶时受交通流限制较小;在中速区间速度略低,怠速比例显著偏高;在高速区间速度偏高。中、高速区间轻型货车运动片段持续时间显著偏低,与其相对频繁的装、卸货有关。3个区间轻型货车的匀速比例都显著偏低。此外,两者在各区间的样本比例有一定差异。
综上所述,有必要为轻型货车和重型货车分别构建工况曲线。
2.2.3重型商用车辆行驶工况体系
最终确定的中国重型商用车辆行驶工况包含6条工况曲线,分别为:中国城市客车行驶工况(China Heavy-Duty Commercial Vehicle Test Cycle for Bus,CHTC-B)、中国普通客车行驶工况(China Heavy-Duty Commercial Vehicle Test Cycle for Coach,CHTC-C)、中国轻型货车行驶工况(China
Heavy-Duty Commercial Vehicle Test Cycle for Light Truck,CHTC-LT)、中国重型货车行驶工况(China Heavy-Duty Commercial Vehicle Test Cycle for Heavy Truck,CHTC-HT)、中国自卸汽车行驶工况(China Heavy-Duty Commercial Vehicle Test Cycle for Dumper,CHTC-D)、中国半挂牵引车行驶工况(China Heavy-Duty Commercial Vehicle Test Cycle for Tractor-Trailer,CHTC-TT)。
2.3曲线结构确定
以城市客车为例,对CHTC中曲线的具体构建过程进行说明。
2.3.1速度区间划分
为确定CHTC-B的曲线结构,根据片段特征对城市客车的全国短片段库进行聚类,划分成2个区间。区间的主要工况特征指标差异显著:低速区间平均速度为8km/h,高速区间为20km/h,随着速度提高,运动片段持续时间和里程增加,怠速比例显著降低,加速度参数略有升高,可以认为2个区间分别反映城市客车在拥堵和顺畅交通状况下的运行情况。根据聚类特征将城市客车在每个城市采集的短片段划分进2个区间,进而计算各城市各区间的综合特征参数,如图7所示。
从图7中可以看出,同一城市在2个区间的特征参数差异显著,不同城市在同一区间基本在同一范围内,验证了区间划分的合理性。
最终确定CHTC-B曲线由低速和高速2个工况区间构成。
2.3.2速度区间时间比例确定
城市客车在各城市运行特征差异显著,为了反映全国综合行驶工况,需引入代表城市客车交通流量的城市权重,结合各城市数据库信息来确定曲线的区间时长比例。
(a)平均速度
(b)怠速比例
图7城市客车各城市、区间特征参数对比
城市客车每条路线每天的运行量相近,调研各城市所有线路数量作为评估城市客车交通流量的参考指标,其与所有城市的路线数量总和的比值即为该城市权重,如表2所示。
表2各城市的城市客车线路数及权重
提取各城市在各区间的城市客车短片段库,统计行驶总时长,与该城市的行驶总时长的比值即为城市的区间权重,如表3所示。
表3各城市的速度区间权重%对于2个速度区间,各城市权重与该城市该区间权
特北
广
西
口25
20
15
10
5
/
k
m
·
h
-
1
低速区高速区总体均值
广
西
口50
40
30
20
10
/
%
低速区高速区总体均值
城市
北京
广州
南宁
线路数量/条
2004
1632
201
线路数量权重/%
14.78
12.04
1.48
城市
北京
广州
南宁
速度区间权重
低速区间
20.13
29.14
38.39
高速区间
79.78
70.86
61.61
-
-54
汽车技术
重的乘积按城市累加即得到该区间在工况曲线中所占时长比例。低速、高速区间时长比例分别为30.5%和69.5%。
2.4工况曲线构建
CHTC-B持续时间设置为1310s,既能满足统计学的代表性,又具备在试验室内进行排放和油耗测试的可行性。
2.4.1区间短片段时长和数量确定
工况曲线按各区间权重进行分配,得到低速区间时长t l=399s,高速区间时长t h=911s。
统计各城市在低速区间数据库中运行片段持续时长的平均值与该城市权重的乘积,按城市累加即可得到低速区间运动片段平均时长t st,l=36.39s;同理计算高速区间运动片段平均时长t st,h=69s;
同理,计算得到低、高速区间怠速片段平均时长分别为t id,l=7.31s、t id,h=13.32s。
结合区间时长,计算低速区间的运动片段数量n st,l 及怠速片段数量n id,l:
n st,l=(t l-t st,l)/(t st,l+t id,l)≈7(1)
n id,l=n st,l+1=8(2)计算低速片段库中运动片段时长的累计分布频率,以此为基础将片段库均分为n st,l个区间,每个区间中样本占比为1/n st,l。对第i(i=0,1,2,…,n st,l-1)个频率区间(i/n st,l,(i+1)/n st,l],定位其第50百分位对应的片段时长,作为构成低速工况第i个运动片段的选择基准。最终确定区间的7个运动片段的时长依次为15s、21s、27s、33s、40s、49s和64s。
同理,确定高速区间由8个运动片段和9个怠速片段构成,运动片段时长依次为30s、37s、43s、49s、55s、61s、68s、76s、85s、96s和117s。
类似地,根据各区间怠速片段时长的累计频率确定怠速片段的时长集合。
2.4.2区间短片段组合选择
对于特定速度区间,获取各城市在该区间运行数据库的速度-加速度的联合分布(v-a分布),用相应城市权重以及城市在该区间的权重对其进行加权,可以生成能够反映该区间全国工况的统一v-a分布,作为工况片段选择的基础。
根据确定的片段时长基准,从库中挑选出候选短片段集合,按照笛卡尔积进行自由组合,将组合后片段的v-a分布与速度区间的统一v-a分布进行卡方检验:
p=∑(V v-a-W v-a)/W v-a(3)式中,p为卡方检验值;V v-a为候选片段组合的v-a分布值;W v-a为区间的统一v-a分布值。
选择卡方检验结果最优的,即卡方检验值较小的片段组合。
2.4.3工况输出
从短到长交替排列最优片段组合中的运动和怠速片段作为区间工况曲线。将低速、高速区间曲线排序组合,生成城市客车工况,如图8所示。
图8城市客车工况曲线
采用类似方法为其他重型商用车构建工况曲线。最终完成包含6条工况曲线的中国重型商用车行驶工况[7]。
3与实际采集数据及其他典型工况的对比计算CHTC体系中各条曲线的关键参数对其特征进行描述。为了说明工况的可靠性,与曲线对应车型的实际采集数据库进行对比;为了说明工况开发的必要性,与其他已有典型工况进行对比。
3.1客车工况特征对比
图9所示为CHTC-B曲线与该车型实际数据库和典型工况的关键特征对比结果。典型工况包括按车型特征里程比例加权后的C-WTVC工况及中国城市客车循环(China City Bus Cycle,CCBC)工况[8]。
(a)速度参数
(b)时间比例
图9CHTC-B、实际数据、CCBC、C-WTVC特征对比
200400600800100012001400
时间/s
80
60
40
20
/
k
m
·
h
-
1
1310s
1部(399s)2部(911s )
总体分布
30
20
10
/
k
m
·
h
-
1
平均速度
平均正加速度
沈阳交通违章查询
CHTC-B CCBC C-WTVC城区
0.8
0.6
0.4
0.2
/
m
·
蒸发器原理
s
-
2
运行平均速度
平均负加速度
0102030405060708090100
时间比例/%
C-WTVC城区
CCBC
CHTC-B
总体分布
加速减速匀速怠速
16.67%
28.40%
23.99%
24.21%
-
-55
2021年第2期