CN114707366A - 车辆行驶阻力计算方法及系统 - Google Patents

车辆行驶阻力计算方法及系统 Download PDF

Info

Publication number
CN114707366A
CN114707366A CN202210627373.6A CN202210627373A CN114707366A CN 114707366 A CN114707366 A CN 114707366A CN 202210627373 A CN202210627373 A CN 202210627373A CN 114707366 A CN114707366 A CN 114707366A
Authority
CN
China
Prior art keywords
vehicle
resistance
test
vehicle speed
temperature
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
CN202210627373.6A
Other languages
English (en)
Other versions
CN114707366B (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.)
Jiangling Motors Corp Ltd
Original Assignee
Jiangling Motors Corp Ltd
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 Jiangling Motors Corp Ltd filed Critical Jiangling Motors Corp Ltd
Priority to CN202210627373.6A priority Critical patent/CN114707366B/zh
Publication of CN114707366A publication Critical patent/CN114707366A/zh
Application granted granted Critical
Publication of CN114707366B publication Critical patent/CN114707366B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08—Thermal analysis or thermal optimisation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
    • 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
    • Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Algebra (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Operations Research (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Biology (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Probability & Statistics with Applications (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • Control Of Transmission Device (AREA)
  • Control Of Driving Devices And Active Controlling Of Vehicle (AREA)

Abstract

本发明公开了一种车辆行驶阻力计算方法及系统。该车辆行驶阻力计算方法，包括：获取车辆的车速v、车辆行驶的环境温度T、实时车重m、车辆的行驶坡度η；根据车速v及实时车重m，计算基础阻力F；根据环境温度T和车速v，在预先标定的三轴MAP中，通过查表的方式，得出温度补偿系数，并根据温度补偿系数对基础阻力F进行修正，以得到修正后的基础阻力F1；根据行驶坡度η计算坡度阻力Fp；根据修正后的基础阻力F1和坡度阻力Fp计算出最终的整车阻力Ff。本发明能够解决现有技术计算出的行驶阻力与实际行驶阻力存在较大误差的技术问题。

Description

车辆行驶阻力计算方法及系统
技术领域
本发明涉及数据处理技术领域，特别是涉及一种车辆行驶阻力计算方法及系统。
背景技术
在车辆行驶过程中，会受到空气阻力、滚动阻力及坡度阻力等多方面的阻力影响，获取相对精准且实时的车辆阻力可把控车辆的行驶状态，并为车速计算、加速度计算、所需驱动力计算等一系列后续计算过程提供数据输入。
当前获取车辆实时阻力的方法主要是通过输入的车速，并根据滑行阻力曲线查表得到，其中GB/T 12534、GB 18352、ETA-HTP01、SAE J1263、SAE J2263等标准系统的介绍了道路滑行阻力的测试方法。通常汽车道路滑行阻力测试可以使用带GPS定位的数据采集系统（如VBOX、Dewetron等）在平直、无风（或者风速不足以干扰试验结果）的试验道路进行，记录在高速状态断开动力系统输出使得车辆由目标车速空档滑行至低车速甚至0车速的过程，采集滑行全程数据数据。
在测试数据处理方面，测试采集空档滑行过程中的时间和车速，通过微分处理进一步获取车辆在滑行过程中的阻力和滑行车速之间的二次函数关系，即有y=Cx2+Bx+A的关系，其中y是滑行阻力，x为车速，C、B、A分别是阻力系数。通过测试获取C、B、A系数后只需要输入车速即可得到对应的阻力。
但该种方法存在较大局限性，首先随着车的总质量增加，车轮对路面的压力会增加，在路面滚动摩擦系数不变的情况下车辆所受的滚动阻力会出现增加，其次车辆上坡时需克服额外的坡度阻力，环境温度的影响也会影响路面的滚动阻力，随着环境温度的下降，滚动阻力也会出现相应的增加，以上三个变量在当前行驶阻力的计算过程中均未纳入考虑，这会导致计算出的行驶阻力与实际行驶阻力存在较大误差。
发明内容
为此，本发明的一个实施例提出一种车辆行驶阻力计算方法，以解决现有技术计算出的行驶阻力与实际行驶阻力存在较大误差的技术问题。
根据本发明一实施例的车辆行驶阻力计算方法，包括：
获取车辆的车速v、车辆行驶的环境温度T、实时车重m、车辆的行驶坡度η；
根据车速v及实时车重m，计算基础阻力F；
根据环境温度T和车速v，在预先标定的三轴MAP中，通过查表的方式，得出温度补偿系数，并根据温度补偿系数对基础阻力F进行修正，以得到修正后的基础阻力F1；
根据行驶坡度η计算坡度阻力Fp；
根据修正后的基础阻力F1和坡度阻力Fp计算出最终的整车阻力Ff。
根据本发明实施例提供的车辆行驶阻力计算方法，首先获取车辆的车速v、车辆行驶的环境温度T、实时车重m、车辆的行驶坡度η，然后根据车速v及实时车重m，计算基础阻力F，再在预先标定的三轴MAP中，通过查表的方式，得出环境温度T和车速v对应的温度补偿系数，以对基础阻力F进行修正，最后根据修正后的基础阻力F1和坡度阻力Fp计算出最终的整车阻力Ff，该整车阻力Ff综合考虑了车速、环境温度、实时车重以及行驶坡度的影响，即通过车速、环境温度、实时车重及行驶坡度对行驶阻力进行相关校准，能够有效提高行驶阻力计算的准确性。
此外，根据本发明实施例提供的车辆行驶阻力计算方法，还具有以下技术特征：
进一步地，所述三轴MAP通过以下方式标定得到：
S1，选取平直路面标定环境，测试车辆半载；
S2，将测试车辆加速至测试车辆的最高车速，随后将档位挂入N档，使测试车辆在断开动力系统输出的条件下由最高车速滑行至完全静止，滑行全程记录车速随时间变化的曲线；
S3，调转测试车辆车头至与上一次测试车辆行驶方向相反，在相同的测试路面上重复步骤S1~S2的测试过程并记录车速随时间变化的曲线，两次测试结束后取两次平均值作为最终的车速随时间变化的曲线；
S4，通过对最终的车速随时间变化的曲线中的车速进行微分计算，得出滑行过程中加速度与车速的关系曲线，加速度与车速的关系曲线中，加速度与测试车辆的车重的乘积即为测试车辆的基础阻力，将测试车辆的基础阻力与车速相关联，即得到测试车辆的基础阻力与车速之间的二轴MAP；
S5，在测试车辆允许工作的环境温度下从最低温度至最高温度，每间隔预设温度重复步骤S1~S4的测试过程，以得出测试车辆的基础阻力与车速及环境温度三者之间的关系曲线；
S6，根据标定车速V标及半载工况下测试车辆的车重m半，计算测试车辆的标准阻力F标；
S7，将标准阻力F标和标定车速V标代入测试车辆的基础阻力与车速及环境温度三者之间的关系曲线中，得出基准温度；
S8，根据测试车辆在基准温度以外的环境温度下的阻力与标准阻力F标计算温度补偿系数；
S9，将环境温度、车速以及温度补偿系数相关联，以得到三轴MAP。
进一步地，基础阻力F的计算公式如下：
Figure 263343DEST_PATH_IMAGE001
其中，CD为车辆的风阻系数，A为车辆前向行驶时的迎风面积，a、b均为滚动阻力对应的相关系数，g表示重力加速度；
修正后的基础阻力F1的计算公式如下：
F1=β*F；
其中，β表示在预先标定的三轴MAP中，环境温度T和车速v对应的温度补偿系数。
进一步地，实时车重m通过以下方式得到：
获取轴荷传感器采集到的前轴轴荷m1及后轴轴荷m2；
通过预先标定的整车簧下质量mh，结合前轴轴荷m1、后轴轴荷m2，计算实时车重m；
其中，m=m1+m2+mh；
坡度阻力Fp的计算公式如下：
Fp=m×g×sinη
其中，g表示重力加速度。
进一步地，整车阻力Ff的计算公式如下：
Ff= F1+Fp。
本发明的另一个实施例提出一种车辆行驶阻力计算系统，以解决现有技术计算出的行驶阻力与实际行驶阻力存在较大误差的技术问题。
根据本发明实施例的车辆行驶阻力计算系统，包括：
获取模块，用于获取车辆的车速v、车辆行驶的环境温度T、实时车重m、车辆的行驶坡度η；
第一计算模块，用于根据车速v及实时车重m，计算基础阻力F；
查找修正模块，用于根据环境温度T和车速v，在预先标定的三轴MAP中，通过查表的方式，得出温度补偿系数，并根据温度补偿系数对基础阻力F进行修正，以得到修正后的基础阻力F1；
第二计算模块，用于根据行驶坡度η计算坡度阻力Fp；
第三计算模块，用于根据修正后的基础阻力F1和坡度阻力Fp计算出最终的整车阻力Ff。
根据本发明实施例提供的车辆行驶阻力计算系统，首先获取车辆的车速v、车辆行驶的环境温度T、实时车重m、车辆的行驶坡度η，然后根据车速v及实时车重m，计算基础阻力F，再在预先标定的三轴MAP中，通过查表的方式，得出环境温度T和车速v对应的温度补偿系数，以对基础阻力F进行修正，最后根据修正后的基础阻力F1和坡度阻力Fp计算出最终的整车阻力Ff，该整车阻力Ff综合考虑了车速、环境温度、实时车重以及行驶坡度的影响，即通过车速、环境温度、实时车重及行驶坡度对行驶阻力进行相关校准，能够有效提高行驶阻力计算的准确性。
此外，根据本发明实施例提供的车辆行驶阻力计算系统，还具有以下技术特征：
进一步地，所述三轴MAP通过以下方式标定得到：
S1，选取平直路面标定环境，测试车辆半载；
S2，将测试车辆加速至测试车辆的最高车速，随后将档位挂入N档，使测试车辆在断开动力系统输出的条件下由最高车速滑行至完全静止，滑行全程记录车速随时间变化的曲线；
S3，调转测试车辆车头至与上一次测试车辆行驶方向相反，在相同的测试路面上重复步骤S1~S2的测试过程并记录车速随时间变化的曲线，两次测试结束后取两次平均值作为最终的车速随时间变化的曲线；
S4，通过对最终的车速随时间变化的曲线中的车速进行微分计算，得出滑行过程中加速度与车速的关系曲线，加速度与车速的关系曲线中，加速度与测试车辆的车重的乘积即为测试车辆的基础阻力，将测试车辆的基础阻力与车速相关联，即得到测试车辆的基础阻力与车速之间的二轴MAP；
S5，在测试车辆允许工作的环境温度下从最低温度至最高温度，每间隔预设温度重复步骤S1~S4的测试过程，以得出测试车辆的基础阻力与车速及环境温度三者之间的关系曲线；
S6，根据标定车速V标及半载工况下测试车辆的车重m半，计算测试车辆的标准阻力F标；
S7，将标准阻力F标和标定车速V标代入测试车辆的基础阻力与车速及环境温度三者之间的关系曲线中，得出基准温度；
S8，根据测试车辆在基准温度以外的环境温度下的阻力与标准阻力F标计算温度补偿系数；
S9，将环境温度、车速以及温度补偿系数相关联，以得到三轴MAP。
进一步地，基础阻力F的计算公式如下：
Figure 725548DEST_PATH_IMAGE002
其中，CD为车辆的风阻系数，A为车辆前向行驶时的迎风面积，a、b均为滚动阻力对应的相关系数，g表示重力加速度；
修正后的基础阻力F1的计算公式如下：
F1=β*F；
其中，β表示在预先标定的三轴MAP中，环境温度T和车速v对应的温度补偿系数。
进一步地，实时车重m通过以下方式得到：
获取轴荷传感器采集到的前轴轴荷m1及后轴轴荷m2；
通过预先标定的整车簧下质量mh，结合前轴轴荷m1、后轴轴荷m2，计算实时车重m；
其中，m=m1+m2+mh；
坡度阻力Fp的计算公式如下：
Fp=m×g×sinη
其中，g表示重力加速度。
进一步地，整车阻力Ff的计算公式如下：
Ff= F1+Fp。
本发明的附加方面和优点将在下面的描述中部分给出，部分将从下面的描述中变得明显，或通过本发明的实施例了解到。
附图说明
本发明实施例的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解，其中：
图1是根据本发明一实施例的车辆行驶阻力计算方法的流程图；
图2是根据本发明一实施例的车辆行驶阻力计算系统的结构框图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
请参阅图1，本发明一实施例提出的车辆行驶阻力计算方法，包括步骤S101~S105：
S101，获取车辆的车速v、车辆行驶的环境温度T、实时车重m、车辆的行驶坡度η。
其中，可以通过车速传感器采集车速v，通过环境温度传感器采集环境温度T。
通过坡度传感器采集车辆行驶路面与水平面的夹角η即为行驶坡度，若车辆上坡行驶则坡度η为正值，若车辆下坡行驶则坡度η为负值。
其中，实时车重m通过以下方式得到：
获取轴荷传感器采集到的前轴轴荷m1及后轴轴荷m2；
通过预先标定的整车簧下质量mh，结合前轴轴荷m1、后轴轴荷m2，计算实时车重m；
其中，m=m1+m2+mh。整车簧下质量mh是指悬吊系统包括车轮、车轴等的总质量，可根据车辆物理属性进行标定。
S102，根据车速v及实时车重m，计算基础阻力F。
其中，基础阻力F的计算公式如下：
Figure 322620DEST_PATH_IMAGE002
其中，CD为车辆的风阻系数，根据实车物理属性进行标定，A为车辆前向行驶时的迎风面积，根据车辆物理属性进行标定，a、b均为滚动阻力对应的相关系数，可根据实车车型及物理属性对a和b进行标定，例如实车车型为小型乘用车时，a取0.0116，b取0.000142，g表示重力加速度。
S103，根据环境温度T和车速v，在预先标定的三轴MAP中，通过查表的方式，得出温度补偿系数，并根据温度补偿系数对基础阻力F进行修正，以得到修正后的基础阻力F1。
其中，所述三轴MAP通过以下方式标定得到：
S1，选取平直路面标定环境，标定场地要求无风或微风（风速不足以干扰试验结果），测试过程中维持环境温度无过大变化，测试车辆半载；
S2，将测试车辆加速至测试车辆的最高车速，随后将档位挂入N档，使测试车辆在断开动力系统输出的条件下由最高车速滑行至完全静止，滑行全程记录车速随时间变化的曲线；
S3，调转测试车辆车头至与上一次测试车辆行驶方向相反，在相同的测试路面上重复步骤S1~S2的测试过程并记录车速随时间变化的曲线，两次测试结束后取两次平均值作为最终的车速随时间变化的曲线；
S4，通过对最终的车速随时间变化的曲线中的车速进行微分计算，得出滑行过程中加速度与车速的关系曲线，加速度与车速的关系曲线中，加速度与测试车辆的车重的乘积即为测试车辆的基础阻力，将测试车辆的基础阻力与车速相关联，即得到测试车辆的基础阻力与车速之间的二轴MAP；
S5，在测试车辆允许工作的环境温度下从最低温度至最高温度，每间隔预设温度（例如每间隔10℃）重复步骤S1~S4的测试过程，以得出测试车辆的基础阻力与车速及环境温度三者之间的关系曲线；
S6，根据标定车速V标及半载工况下测试车辆的车重m半，计算测试车辆的标准阻力F标；
S7，将标准阻力F标和标定车速V标代入测试车辆的基础阻力与车速及环境温度三者之间的关系曲线中，得出基准温度；
S8，根据测试车辆在基准温度以外的环境温度下的阻力与标准阻力F标计算温度补偿系数。具体用测试车辆在基准温度以外的环境温度下的阻力除以标准阻力F标就得到了相应温度的温度补偿系数，例如基准温度是25℃，其对应标准阻力F标，测试车辆在26℃的环境温度下的阻力为F26，则F26/F标就是26℃的环境温度对应的温度补偿系数；
S9，将环境温度、车速以及温度补偿系数相关联，以得到三轴MAP。
根据环境温度T和车速v，在预先标定的三轴MAP中，通过查表的方式，能够得出环境温度T和车速v对应的温度补偿系数β，则修正后的基础阻力F1的计算公式如下：
F1=β*F。
S104，根据行驶坡度η计算坡度阻力Fp。
其中，坡度阻力Fp的计算公式如下：
Fp=m×g×sinη
其中，g表示重力加速度。
S105，根据修正后的基础阻力F1和坡度阻力Fp计算出最终的整车阻力Ff。
其中，整车阻力Ff的计算公式如下：
Ff= F1+Fp。
根据本发明提供的车辆行驶阻力计算方法，首先获取车辆的车速v、车辆行驶的环境温度T、实时车重m、车辆的行驶坡度η，然后根据车速v及实时车重m，计算基础阻力F，再在预先标定的三轴MAP中，通过查表的方式，得出环境温度T和车速v对应的温度补偿系数，以对基础阻力F进行修正，最后根据修正后的基础阻力F1和坡度阻力Fp计算出最终的整车阻力Ff，该整车阻力Ff综合考虑了车速、环境温度、实时车重以及行驶坡度的影响，即通过车速、环境温度、实时车重及行驶坡度对行驶阻力进行相关校准，能够有效提高行驶阻力计算的准确性。
请参阅图2，本发明一实施例提出的车辆行驶阻力计算系统，包括：
获取模块，用于获取车辆的车速v、车辆行驶的环境温度T、实时车重m、车辆的行驶坡度η；
第一计算模块，用于根据车速v及实时车重m，计算基础阻力F；
查找修正模块，用于根据环境温度T和车速v，在预先标定的三轴MAP中，通过查表的方式，得出温度补偿系数，并根据温度补偿系数对基础阻力F进行修正，以得到修正后的基础阻力F1；
第二计算模块，用于根据行驶坡度η计算坡度阻力Fp；
第三计算模块，用于根据修正后的基础阻力F1和坡度阻力Fp计算出最终的整车阻力Ff。
本实施例中，所述三轴MAP通过以下方式标定得到：
S1，选取平直路面标定环境，测试车辆半载；
S2，将测试车辆加速至测试车辆的最高车速，随后将档位挂入N档，使测试车辆在断开动力系统输出的条件下由最高车速滑行至完全静止，滑行全程记录车速随时间变化的曲线；
S3，调转测试车辆车头至与上一次测试车辆行驶方向相反，在相同的测试路面上重复步骤S1~S2的测试过程并记录车速随时间变化的曲线，两次测试结束后取两次平均值作为最终的车速随时间变化的曲线；
S4，通过对最终的车速随时间变化的曲线中的车速进行微分计算，得出滑行过程中加速度与车速的关系曲线，加速度与车速的关系曲线中，加速度与测试车辆的车重的乘积即为测试车辆的基础阻力，将测试车辆的基础阻力与车速相关联，即得到测试车辆的基础阻力与车速之间的二轴MAP；
S5，在测试车辆允许工作的环境温度下从最低温度至最高温度，每间隔预设温度重复步骤S1~S4的测试过程，以得出测试车辆的基础阻力与车速及环境温度三者之间的关系曲线；
S6，根据标定车速V标及半载工况下测试车辆的车重m半，计算测试车辆的标准阻力F标；
S7，将标准阻力F标和标定车速V标代入测试车辆的基础阻力与车速及环境温度三者之间的关系曲线中，得出基准温度；
S8，根据测试车辆在基准温度以外的环境温度下的阻力与标准阻力F标计算温度补偿系数；
S9，将环境温度、车速以及温度补偿系数相关联，以得到三轴MAP。
本实施例中，基础阻力F的计算公式如下：
Figure 325211DEST_PATH_IMAGE002
其中，CD为车辆的风阻系数，A为车辆前向行驶时的迎风面积，a、b均为滚动阻力对应的相关系数，g表示重力加速度；
修正后的基础阻力F1的计算公式如下：
F1=β*F；
其中，β表示在预先标定的三轴MAP中，环境温度T和车速v对应的温度补偿系数。
本实施例中，实时车重m通过以下方式得到：
获取轴荷传感器采集到的前轴轴荷m1及后轴轴荷m2；
通过预先标定的整车簧下质量mh，结合前轴轴荷m1、后轴轴荷m2，计算实时车重m；
其中，m=m1+m2+mh；
坡度阻力Fp的计算公式如下：
Fp=m×g×sinη
其中，g表示重力加速度。
本实施例中，整车阻力Ff的计算公式如下：
Ff= F1+Fp。
根据本发明提供的车辆行驶阻力计算系统，首先获取车辆的车速v、车辆行驶的环境温度T、实时车重m、车辆的行驶坡度η，然后根据车速v及实时车重m，计算基础阻力F，再在预先标定的三轴MAP中，通过查表的方式，得出环境温度T和车速v对应的温度补偿系数，以对基础阻力F进行修正，最后根据修正后的基础阻力F1和坡度阻力Fp计算出最终的整车阻力Ff，该整车阻力Ff综合考虑了车速、环境温度、实时车重以及行驶坡度的影响，即通过车速、环境温度、实时车重及行驶坡度对行驶阻力进行相关校准，能够有效提高行驶阻力计算的准确性。
在流程图中表示或在此以其他方式描述的逻辑和/或步骤，例如，可以被认为是用于实现逻辑功能的可执行指令的定序列表，可以具体实现在任何计算机可读介质中，以供指令执行系统、装置或设备（如基于计算机的系统、包括处理器的系统或其他可以从指令执行系统、装置或设备取指令并执行指令的系统）使用，或结合这些指令执行系统、装置或设备而使用。就本说明书而言，“计算机可读介质”可以是任何可以包含、存储、通讯、传播或传输程序以供指令执行系统、装置或设备或结合这些指令执行系统、装置或设备而使用的装置。
计算机可读介质的更具体的示例（非穷尽性列表）包括以下：具有一个或多个布线的电连接部（电子装置），便携式计算机盘盒（磁装置），随机存取存储器（RAM），只读存储器（ROM），可擦除可编辑只读存储器（EPROM或闪速存储器），光纤装置，以及便携式光盘只读存储器（CDROM）。另外，计算机可读介质甚至可以是可在其上打印所述程序的纸或其他合适的介质，因为可以例如通过对纸或其他介质进行光学扫描，接着进行编辑、解译或必要时以其他合适方式进行处理来以电子方式获得所述程序，然后将其存储在计算机存储器中。
应当理解，本发明的各部分可以用硬件、软件、固件或它们的组合来实现。在上述实施方式中，多个步骤或方法可以用存储在存储器中且由合适的指令执行系统执行的软件或固件来实现。例如，如果用硬件来实现，和在另一实施方式中一样，可用本领域公知的下列技术中的任一项或他们的组合来实现：具有用于对数据信号实现逻辑功能的逻辑门电路的离散逻辑电路，具有合适的组合逻辑门电路的专用集成电路，可编程门阵列（PGA），现场可编程门阵列（FPGA）等。
在本说明书的描述中，参考术语“一个实施例”、“一些实施例”、 “示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
尽管已经示出和描述了本发明的实施例，本领域的普通技术人员可以理解：在不脱离本发明的原理和宗旨的情况下可以对这些实施例进行多种变化、修改、替换和变型，本发明的范围由权利要求及其等同物限定。

Claims (10)

1.一种车辆行驶阻力计算方法，其特征在于，包括：
获取车辆的车速v、车辆行驶的环境温度T、实时车重m、车辆的行驶坡度η；
根据车速v及实时车重m，计算基础阻力F；
根据环境温度T和车速v，在预先标定的三轴MAP中，通过查表的方式，得出温度补偿系数，并根据温度补偿系数对基础阻力F进行修正，以得到修正后的基础阻力F1；
根据行驶坡度η计算坡度阻力Fp；
根据修正后的基础阻力F1和坡度阻力Fp计算出最终的整车阻力Ff。
2.根据权利要求1所述的车辆行驶阻力计算方法，其特征在于，所述三轴MAP通过以下方式标定得到：
S1，选取平直路面标定环境，测试车辆半载；
S2，将测试车辆加速至测试车辆的最高车速，随后将档位挂入N档，使测试车辆在断开动力系统输出的条件下由最高车速滑行至完全静止，滑行全程记录车速随时间变化的曲线；
S3，调转测试车辆车头至与上一次测试车辆行驶方向相反，在相同的测试路面上重复步骤S1~S2的测试过程并记录车速随时间变化的曲线，两次测试结束后取两次平均值作为最终的车速随时间变化的曲线；
S4，通过对最终的车速随时间变化的曲线中的车速进行微分计算，得出滑行过程中加速度与车速的关系曲线，加速度与车速的关系曲线中，加速度与测试车辆的车重的乘积即为测试车辆的基础阻力，将测试车辆的基础阻力与车速相关联，即得到测试车辆的基础阻力与车速之间的二轴MAP；
S5，在测试车辆允许工作的环境温度下从最低温度至最高温度，每间隔预设温度重复步骤S1~S4的测试过程，以得出测试车辆的基础阻力与车速及环境温度三者之间的关系曲线；
S6，根据标定车速V标及半载工况下测试车辆的车重m半，计算测试车辆的标准阻力F标；
S7，将标准阻力F标和标定车速V标代入测试车辆的基础阻力与车速及环境温度三者之间的关系曲线中，得出基准温度；
S8，根据测试车辆在基准温度以外的环境温度下的阻力与标准阻力F标计算温度补偿系数；
S9，将环境温度、车速以及温度补偿系数相关联，以得到三轴MAP。
3.根据权利要求1所述的车辆行驶阻力计算方法，其特征在于，基础阻力F的计算公式如下：
Figure 252238DEST_PATH_IMAGE001
其中，CD为车辆的风阻系数，A为车辆前向行驶时的迎风面积，a、b均为滚动阻力对应的相关系数，g表示重力加速度；
修正后的基础阻力F1的计算公式如下：
F1=β*F；
其中，β表示在预先标定的三轴MAP中，环境温度T和车速v对应的温度补偿系数。
4.根据权利要求1所述的车辆行驶阻力计算方法，其特征在于，实时车重m通过以下方式得到：
获取轴荷传感器采集到的前轴轴荷m1及后轴轴荷m2；
通过预先标定的整车簧下质量mh，结合前轴轴荷m1、后轴轴荷m2，计算实时车重m；
其中，m=m1+m2+mh；
坡度阻力Fp的计算公式如下：
Fp=m×g×sinη
其中，g表示重力加速度。
5.根据权利要求4所述的车辆行驶阻力计算方法，其特征在于，整车阻力Ff的计算公式如下：
Ff= F1+Fp。
6.一种车辆行驶阻力计算系统，其特征在于，包括：
获取模块，用于获取车辆的车速v、车辆行驶的环境温度T、实时车重m、车辆的行驶坡度η；
第一计算模块，用于根据车速v及实时车重m，计算基础阻力F；
查找修正模块，用于根据环境温度T和车速v，在预先标定的三轴MAP中，通过查表的方式，得出温度补偿系数，并根据温度补偿系数对基础阻力F进行修正，以得到修正后的基础阻力F1；
第二计算模块，用于根据行驶坡度η计算坡度阻力Fp；
第三计算模块，用于根据修正后的基础阻力F1和坡度阻力Fp计算出最终的整车阻力Ff。
7.根据权利要求6所述的车辆行驶阻力计算系统，其特征在于，所述三轴MAP通过以下方式标定得到：
S1，选取平直路面标定环境，测试车辆半载；
S2，将测试车辆加速至测试车辆的最高车速，随后将档位挂入N档，使测试车辆在断开动力系统输出的条件下由最高车速滑行至完全静止，滑行全程记录车速随时间变化的曲线；
S3，调转测试车辆车头至与上一次测试车辆行驶方向相反，在相同的测试路面上重复步骤S1~S2的测试过程并记录车速随时间变化的曲线，两次测试结束后取两次平均值作为最终的车速随时间变化的曲线；
S4，通过对最终的车速随时间变化的曲线中的车速进行微分计算，得出滑行过程中加速度与车速的关系曲线，加速度与车速的关系曲线中，加速度与测试车辆的车重的乘积即为测试车辆的基础阻力，将测试车辆的基础阻力与车速相关联，即得到测试车辆的基础阻力与车速之间的二轴MAP；
S5，在测试车辆允许工作的环境温度下从最低温度至最高温度，每间隔预设温度重复步骤S1~S4的测试过程，以得出测试车辆的基础阻力与车速及环境温度三者之间的关系曲线；
S6，根据标定车速V标及半载工况下测试车辆的车重m半，计算测试车辆的标准阻力F标；
S7，将标准阻力F标和标定车速V标代入测试车辆的基础阻力与车速及环境温度三者之间的关系曲线中，得出基准温度；
S8，根据测试车辆在基准温度以外的环境温度下的阻力与标准阻力F标计算温度补偿系数；
S9，将环境温度、车速以及温度补偿系数相关联，以得到三轴MAP。
8.根据权利要求6所述的车辆行驶阻力计算系统，其特征在于，基础阻力F的计算公式如下：
Figure 44613DEST_PATH_IMAGE001
其中，CD为车辆的风阻系数，A为车辆前向行驶时的迎风面积，a、b均为滚动阻力对应的相关系数，g表示重力加速度；
修正后的基础阻力F1的计算公式如下：
F1=β*F；
其中，β表示在预先标定的三轴MAP中，环境温度T和车速v对应的温度补偿系数。
9.根据权利要求6所述的车辆行驶阻力计算系统，其特征在于，实时车重m通过以下方式得到：
获取轴荷传感器采集到的前轴轴荷m1及后轴轴荷m2；
通过预先标定的整车簧下质量mh，结合前轴轴荷m1、后轴轴荷m2，计算实时车重m；
其中，m=m1+m2+mh；
坡度阻力Fp的计算公式如下：
Fp=m×g×sinη
其中，g表示重力加速度。
10.根据权利要求9所述的车辆行驶阻力计算系统，其特征在于，整车阻力Ff的计算公式如下：
Ff= F1+Fp。
CN202210627373.6A 2022-06-06 2022-06-06 车辆行驶阻力计算方法及系统 Active CN114707366B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210627373.6A CN114707366B (zh) 2022-06-06 2022-06-06 车辆行驶阻力计算方法及系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210627373.6A CN114707366B (zh) 2022-06-06 2022-06-06 车辆行驶阻力计算方法及系统

Publications (2)

Publication Number Publication Date
CN114707366A true CN114707366A (zh) 2022-07-05
CN114707366B CN114707366B (zh) 2022-08-26

Family

ID=82177741

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210627373.6A Active CN114707366B (zh) 2022-06-06 2022-06-06 车辆行驶阻力计算方法及系统

Country Status (1)

Country Link
CN (1) CN114707366B (zh)

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000027674A (ja) * 1998-07-15 2000-01-25 Nissan Motor Co Ltd 車両駆動力制御装置
EP2591967A2 (de) * 2011-11-11 2013-05-15 Volkswagen Aktiengesellschaft Verfahren zum Betreiben eines Fahrzeugs, Steuereinrichtung und Fahrzeug
CN103439040A (zh) * 2013-08-21 2013-12-11 吴明 基准温度汽车传动系台架阻力检测方法
CN109882587A (zh) * 2019-01-31 2019-06-14 中国北方车辆研究所 基于等效阻力系数的特种车辆自动换挡控制系统及方法
CN109990766A (zh) * 2019-04-24 2019-07-09 重庆理工大学 基于路面综合阻力系数的道路坡度检测方法与系统
CN110441068A (zh) * 2019-07-31 2019-11-12 辽宁龙马安泰网络技术有限责任公司 汽车行驶带挡滑行系统阻力检测方法
CN112440999A (zh) * 2019-08-16 2021-03-05 郑州宇通客车股份有限公司 一种新能源车辆及其驱动力修正方法和装置
CN112697429A (zh) * 2020-12-09 2021-04-23 朱卫萍 乘用车空挡传动系阻力检测方法
CN113335290A (zh) * 2021-07-22 2021-09-03 中国第一汽车股份有限公司 一种车辆滚动阻力获取方法、获取模块及存储介质
CN113531114A (zh) * 2021-09-09 2021-10-22 北汽福田汽车股份有限公司 自动变速箱的换挡方法、装置、存储介质、电子设备及车辆
CN114347793A (zh) * 2022-01-24 2022-04-15 摩登汽车有限公司 汽车续航的估算方法及系统

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000027674A (ja) * 1998-07-15 2000-01-25 Nissan Motor Co Ltd 車両駆動力制御装置
EP2591967A2 (de) * 2011-11-11 2013-05-15 Volkswagen Aktiengesellschaft Verfahren zum Betreiben eines Fahrzeugs, Steuereinrichtung und Fahrzeug
CN103439040A (zh) * 2013-08-21 2013-12-11 吴明 基准温度汽车传动系台架阻力检测方法
CN109882587A (zh) * 2019-01-31 2019-06-14 中国北方车辆研究所 基于等效阻力系数的特种车辆自动换挡控制系统及方法
CN109990766A (zh) * 2019-04-24 2019-07-09 重庆理工大学 基于路面综合阻力系数的道路坡度检测方法与系统
CN110441068A (zh) * 2019-07-31 2019-11-12 辽宁龙马安泰网络技术有限责任公司 汽车行驶带挡滑行系统阻力检测方法
CN110702427A (zh) * 2019-07-31 2020-01-17 辽宁龙马安泰网络技术有限责任公司 汽车行驶带挡滑行系统阻力检测方法
CN112440999A (zh) * 2019-08-16 2021-03-05 郑州宇通客车股份有限公司 一种新能源车辆及其驱动力修正方法和装置
CN112697429A (zh) * 2020-12-09 2021-04-23 朱卫萍 乘用车空挡传动系阻力检测方法
CN113335290A (zh) * 2021-07-22 2021-09-03 中国第一汽车股份有限公司 一种车辆滚动阻力获取方法、获取模块及存储介质
CN113531114A (zh) * 2021-09-09 2021-10-22 北汽福田汽车股份有限公司 自动变速箱的换挡方法、装置、存储介质、电子设备及车辆
CN114347793A (zh) * 2022-01-24 2022-04-15 摩登汽车有限公司 汽车续航的估算方法及系统

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
HUI LIU: "Adaptive equivalent consumption minimisation strategy and dynamic control allocation-based optimal power management strategy for four-wheel drive hybrid electric vehicles", 《JOURNAL OF AUTOMOBILE ENGINEERING》 *
刘耀华: "轻型汽车道路行驶阻力的测量及影响因素探讨", 《科技与企业》 *
王功博: "某纯电动轻型物流车能量回收策略优化", 《客车技术与研究》 *
陈寅: "基于转鼓试验的微型汽车传动系统阻力测试与研究", 《中国优秀硕士学位论文全文数据库 工程科技II辑》 *

Also Published As

Publication number Publication date
CN114707366B (zh) 2022-08-26

Similar Documents

Publication Publication Date Title
CN111806449A (zh) 一种纯电动车的整车质量和路面坡度的估算方法
CN102649433B (zh) 用于提高燃料经济性指标计算的道路坡度估计的方法
US5548536A (en) Method for determining quantities which characterize the driving behavior
Stephant et al. Virtual sensor: Application to vehicle sideslip angle and transversal forces
CN113002549B (zh) 一种车辆状态估算方法、装置、设备及存储介质
CN108871788B (zh) 一种自动变速器换挡品质测试台架的校验方法
CN107406076A (zh) 行驶控制装置及行驶控制方法
CN105890914A (zh) 电动汽车续航里程测试方法及装置
CN110053431B (zh) 一种测量车辆轮胎压力的误差补偿方法与装置
CN108819950A (zh) 汽车稳定性控制系统的车速估计方法及系统
CN113799783B (zh) 应用于车辆的道路横向坡度测量方法及系统
CN106198044A (zh) 电动汽车续航里程测试方法及装置
US8948959B2 (en) Method of determining the stress that should be applied to a tyre during an indoor endurance bench test
CN110426011A (zh) 车辆转向角测量系统及方法
CN114707366B (zh) 车辆行驶阻力计算方法及系统
CN105579824B (zh) 用于模拟转弯行驶的方法
CN114620052A (zh) 一种车辆信号处理方法、装置和驾驶设备
CN113063414A (zh) 一种用于视觉惯性slam的车辆动力学预积分构建方法
CN109342078B (zh) 电动四驱车型的经济性测试方法、装置及测试装置
CN115406669A (zh) 一种多轴特种车辆侧翻指标优化方法
CN111695197B (zh) 一种罐车侧翻阈值高可靠动态估计方法
CN114932909A (zh) 一种基于互补滤波实现加速度校正的坡度估算方法
Golsch et al. Improved CFD Methodology for Class 8 Tractor-Trailer Coastdown Correlation
Kitahama et al. Measurement method of normalized cornering stiffness
Stasewitsch et al. „Scope of application of a kinematic model for path tracking control “

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