CN103661395B - 动态道路坡度估计 - Google Patents
动态道路坡度估计 Download PDFInfo
- Publication number
- CN103661395B CN103661395B CN201310390790.4A CN201310390790A CN103661395B CN 103661395 B CN103661395 B CN 103661395B CN 201310390790 A CN201310390790 A CN 201310390790A CN 103661395 B CN103661395 B CN 103661395B
- Authority
- CN
- China
- Prior art keywords
- vehicle
- vehicular system
- instruction
- rge
- estimated
- 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.)
- Active
Links
- 230000001133 acceleration Effects 0.000 claims abstract description 166
- 238000000034 method Methods 0.000 claims description 169
- 230000004044 response Effects 0.000 claims description 17
- 238000009825 accumulation Methods 0.000 claims description 9
- 238000006073 displacement reaction Methods 0.000 claims description 9
- 238000002360 preparation method Methods 0.000 claims description 6
- 238000004422 calculation algorithm Methods 0.000 description 100
- 230000003068 static effect Effects 0.000 description 90
- 230000008859 change Effects 0.000 description 32
- 238000010586 diagram Methods 0.000 description 29
- 238000011156 evaluation Methods 0.000 description 25
- 230000006870 function Effects 0.000 description 25
- 230000033001 locomotion Effects 0.000 description 17
- 238000001914 filtration Methods 0.000 description 14
- 230000006866 deterioration Effects 0.000 description 13
- 230000010354 integration Effects 0.000 description 13
- 230000000694 effects Effects 0.000 description 12
- 230000008569 process Effects 0.000 description 11
- 230000001141 propulsive effect Effects 0.000 description 10
- 230000005540 biological transmission Effects 0.000 description 9
- 230000005484 gravity Effects 0.000 description 7
- 238000013441 quality evaluation Methods 0.000 description 7
- 238000005096 rolling process Methods 0.000 description 7
- 230000035945 sensitivity Effects 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 5
- 239000000446 fuel Substances 0.000 description 5
- 238000005259 measurement Methods 0.000 description 5
- 238000012544 monitoring process Methods 0.000 description 5
- 230000007704 transition Effects 0.000 description 5
- 240000006313 Alternanthera bettzickiana Species 0.000 description 4
- 235000014138 Caesalpinia decapetala Nutrition 0.000 description 4
- 230000008901 benefit Effects 0.000 description 4
- 238000013461 design Methods 0.000 description 4
- 239000011159 matrix material Substances 0.000 description 4
- 230000009467 reduction Effects 0.000 description 4
- 230000001105 regulatory effect Effects 0.000 description 4
- 230000006641 stabilisation Effects 0.000 description 4
- 238000011105 stabilization Methods 0.000 description 4
- 239000000725 suspension Substances 0.000 description 4
- 230000001808 coupling effect Effects 0.000 description 3
- 230000002123 temporal effect Effects 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 101150044627 ECM14 gene Proteins 0.000 description 2
- 238000003745 diagnosis Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000002360 explosive Substances 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000010355 oscillation Effects 0.000 description 2
- 230000001052 transient effect Effects 0.000 description 2
- 101000802640 Homo sapiens Lactosylceramide 4-alpha-galactosyltransferase Proteins 0.000 description 1
- 102100035838 Lactosylceramide 4-alpha-galactosyltransferase Human genes 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 230000001143 conditioned effect Effects 0.000 description 1
- 239000013256 coordination polymer Substances 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 230000008030 elimination Effects 0.000 description 1
- 238000003379 elimination reaction Methods 0.000 description 1
- 238000005562 fading Methods 0.000 description 1
- 238000010304 firing Methods 0.000 description 1
- 238000002955 isolation Methods 0.000 description 1
- 238000011068 loading method Methods 0.000 description 1
- 238000004806 packaging method and process Methods 0.000 description 1
- 230000002572 peristaltic effect Effects 0.000 description 1
- 239000000047 product Substances 0.000 description 1
- 230000001172 regenerating effect Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
Classifications
-
- 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
-
- 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/02—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 ambient conditions
- B60W40/06—Road conditions
- B60W40/076—Slope angle of the road
-
- 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/12—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 parameters of the vehicle itself, e.g. tyre models
- B60W40/13—Load or weight
-
- 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
- B60W2520/00—Input parameters relating to overall vehicle dynamics
- B60W2520/10—Longitudinal speed
-
- 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
- B60W2520/00—Input parameters relating to overall vehicle dynamics
- B60W2520/10—Longitudinal speed
- B60W2520/105—Longitudinal acceleration
-
- 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
- B60W2520/00—Input parameters relating to overall vehicle dynamics
- B60W2520/14—Yaw
-
- 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
- B60W2520/00—Input parameters relating to overall vehicle dynamics
- B60W2520/30—Wheel torque
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Mechanical Engineering (AREA)
- Transportation (AREA)
- Automation & Control Theory (AREA)
- Theoretical Computer Science (AREA)
- Databases & Information Systems (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Software Systems (AREA)
- Data Mining & Analysis (AREA)
- Control Of Driving Devices And Active Controlling Of Vehicle (AREA)
- Electric Propulsion And Braking For Vehicles (AREA)
Abstract
本发明涉及动态道路坡度估计。一种车辆和车辆系统设置有控制器,控制器被配置成:产生指示车辆质量估计的输出。当纵向加速度、车辆速度和横摆角中的至少一个指示发生合格的事件时,车辆质量估计基于纵向加速度和车轮扭矩。控制器还被配置成:基于车辆速度、车轮扭矩和车辆质量估计而产生指示动态道路坡度估计的输出。
Description
技术领域
一个或多个实施例涉及一种车辆系统以及用于估计道路坡度和车辆质量的方法。
背景技术
道路坡度和车辆质量是影响多个车辆控制功能的两个参数。这样的功能包括动力传动系控制、能量管理控制、变速器控制、制动控制及稳定性控制。例如,传统的车辆包括提供用于推进车辆的扭矩的内燃发动机。即使当发动机怠速运转时,发动机也通过变速器将扭矩提供给驱动车轮。这样的扭矩被称为蠕动扭矩。当车辆在斜面上停止时,该蠕动扭矩对由制动系统施加的扭矩(制动扭矩)进行补充,以使车辆保持在静止位置。如果在车辆在斜面上停止时关闭发动机,则发动机的蠕动扭矩可消失,且会需要附加制动扭矩以保持车辆的位置并防止车辆向后滚动。需要的附加制动扭矩的量基于道路坡度和车辆质量。
微混合动力车辆可启用停止/起动策略,以用于在驱动事件期间起动和停止车辆的发动机。如果不需要功率(例如,在交通信号灯处等待时),则关闭发动机。只要需要功率,就自动重新起动发动机。通过避免不必要的发动机怠速运转,将提高车辆的燃料经济性。具有起动/停止功能的微混合动力车辆可调节制动压力,以当在车辆在斜面上停止的同时发动机关闭时保持车辆的位置。这样的微混合动力车辆估计道路坡度和车辆质量,并基于这些估计来调节制动压力。
发明内容
在一个实施例中,一种车辆设置有控制器,控制器被配置成:产生指示车辆质量估计的输出。当纵向加速度、车辆速度和横摆角中的至少一个指示发生合格的事件时,车辆质量估计基于纵向加速度和车轮扭矩。车辆质量估计独立于道路坡度估计。
一种车辆包括控制器,控制器被配置成:当纵向加速度、车辆速度和横摆角中的至少一个指示发生合格的事件时,基于纵向加速度和车轮扭矩产生指示车辆质量估计的输出,其中,车辆质量估计独立于道路坡度估计。
合格的事件包括加速事件和制动事件中的至少一个。
加速事件还包括:指示车辆速度大于推进入口速度阈值的输入;指示纵向加速度大于正加速度阈值的输入;指示横摆角的幅度小于横摆角阈值的输入。
制动事件还包括:指示车辆速度大于制动入口速度阈值的输入;指示纵向加速度小于负加速度阈值的输入;指示横摆角的幅度小于横摆角阈值的输入。
所述车辆还包括:惯性传感器,用于提供指示纵向加速度的信号。
控制器还被配置成:将指示车辆质量估计的输出设置为默认值,直到发生预定数量的合格的事件为止。
控制器还被配置成:使用动态上边界限制车辆质量估计,所述动态上边界的初始值基于最大车辆总重量。
控制器还被配置成:使用动态下边界限制车辆质量估计,所述动态下边界的初始值基于车辆整备重量。
所述动态上边界和动态下边界中的至少一个被配置成响应于合格的事件的累积而收敛。
控制器还被配置成:基于纵向加速度、车轮扭矩和加速度偏移量使用递归最小二乘估计来估计车辆质量估计,其中,加速度偏移量包括横向速度、横向偏移量和车辆俯仰角中的至少一个。
控制器还被配置成:基于车辆质量估计、车辆速度和车轮扭矩而产生指示动态道路坡度估计的输出。
在另一实施例中,一种车辆系统设置有控制器,控制器被配置成:响应于合格的事件,基于纵向加速度和车轮扭矩而产生指示车辆质量估计的输出。控制器还被配置成:基于车辆速度、车轮扭矩和车辆质量估计而产生指示动态道路坡度估计的输出。
在另一实施例中,提供一种用于估计动态道路坡度的方法。响应于合格的事件,基于纵向加速度和车轮扭矩而产生指示车辆质量估计的输出。还基于车辆速度、车轮扭矩和车辆质量估计而产生指示动态道路坡度估计的输出。
一种车辆系统包括控制器,控制器被配置成:响应于合格的事件,基于纵向加速度和车轮扭矩而产生指示车辆质量估计的输出;基于车辆速度、车轮扭矩和车辆质量估计而产生指示动态道路坡度估计的输出。
控制器还被配置成:将指示车辆质量估计的输出设置为默认值,直到发生预定数量的合格的事件为止。
控制器还被配置成:响应于检测到合格的事件而使事件计数加1;当事件计数大于合格的事件阈值时,产生指示车辆质量估计的输出。
控制器还被配置成:使用动态上边界和动态下边界限制车辆质量估计,所述动态上边界的初始值基于最大车辆总重量,所述动态下边界的初始值基于车辆整备重量。
所述动态上边界和动态下边界中的至少一个被配置成响应于合格的事件的累积而收敛。
控制器还被配置成:提供指示动态道路坡度估计的可用性的动态品质因子。
车辆质量估计包括车辆有效负荷。
一种用于估计动态道路坡度的方法,所述方法包括:响应于合格的事件,基于纵向加速度和车轮扭矩而产生指示车辆质量估计的输出;基于车辆速度、车轮扭矩和车辆质量估计而产生指示动态道路坡度估计的输出。
所述方法还包括:将指示车辆质量估计的输出设置为默认值,直到发生预定数量的合格的事件为止。
这样,通过使用事件查找策略独立于道路坡度估计而估计车辆质量,车辆系统提供优点。在特定输入的信噪比高的车辆状况或合格的事件下估计质量。然后,当发生查找事件时,车辆系统使用递归最小二乘(RLS)算法估计车辆质量。在车辆状况对应于质量可能改变时(例如,当车辆停止时),车辆系统重新初始化质量估计算法中的特定参数,以增加算法的灵敏度。车辆系统基于估计的车辆质量来估计动态道路坡度(RGEdyn)。一般来说,与静态RGE和运动学RGE相比,动态RGE的品质状况更低,其原因是动态RGE在进行估计时涉及更多的信号,且每个信号中的噪声和误差将传送到最终RGE输出中。然而,动态RGE不依赖纵向加速度信号,因此,动态RGE的信号覆盖范围不同于静态估计和运动学估计的信号覆盖范围。另外,独立于道路坡度估计,车辆系统估计车辆质量,从而将RGEdyn中的噪声与车辆质量估计隔离。
附图说明
图1是根据一个或多个实施例的具有用于估计道路坡度和车辆质量的车辆系统的车辆的示意图;
图2是图1的车辆被示出为位于以一定坡度倾斜的道路上的侧视图;
图3是图1的车辆的仰视图;
图4是示出根据一个或多个实施例的如由图1的车辆系统执行的用于仲裁道路坡度估计的方法的示意性框图;
图5A至图5E是示出由图4的方法控制的道路坡度估计的时间图;
图6是示出根据一个或多个实施例的用于仲裁道路坡度估计的方法的流程图;
图7是示出车辆速度和加速度的时间图;
图8是示出根据一个或多个实施例的用于评价加速度计收敛状态的方法的流程图;
图9是示出根据一个或多个实施例的用于锁定静态道路坡度输出的方法的流程图;
图10A至图10C是示出如由图9的方法控制的车辆速度、加速度和道路坡度估计的时间图;
图11是示出根据一个或多个实施例的用于运动学道路坡度估计的方法的示意性框图;
图12是示出根据一个或多个实施例的用于估计运动学道路坡度的方法的流程图;
图13A至图13D是示出如由图12的方法控制的车辆速度、加速度和道路坡度估计的时间图;
图14是示出根据一个或多个实施例的用于确定车辆驱动状况的方法的流程图;
图15是示出根据一个或多个实施例的用于估计车辆横向速度的方法的流程图;
图16是示出根据一个或多个实施例的用于估计车身侧倾角的方法的流程图;
图17是根据一个或多个实施例的用于选择悬挂参数的方法的视图;
图18是示出根据一个或多个实施例的用于车辆质量估计的方法的示意性框图;
图19是示出根据一个或多个实施例的用于估计车辆质量的方法的流程图;
图20是示出根据一个或多个实施例的用于使车辆质量估计参数复位的方法的流程图;
图21是示出根据一个或多个实施例的用于更新车辆质量估计参数的方法的流程图;
图22是示出车辆质量估计的时间图;
图23是示出车辆质量估计的另一时间图;
图24是示出车辆质量估计的又一时间图;
图25是示出根据一个或多个实施例的用于估计估计器品质的方法的流程图;
图26A至图26C是示出如由图25的方法控制的车辆速度、加速度和道路坡度估计的时间图;
图27A至图27G是示出如由图25的方法控制的车辆速度和品质因子的时间图。
具体实施方式
根据需要,在此公开本发明的具体实施例;然而,应该理解到,公开的实施例仅仅是本发明的示例,本发明可以以多种和可选的形式实施。附图并不一定按照比例绘制;可夸大或最小化一些特征以示出特定部件的细节。因此,在此公开的具体结构和功能性细节不被解释为限制,而仅仅作为用于教导本领域的技术人员以多种方式使用本发明的代表性基础。
参照图1,示出了根据一个或多个实施例的用于估计道路坡度和车辆质量的车辆系统,该车辆系统总体上由标号10指示。车辆系统10被描述为在车辆12内。车辆系统10包括控制器,例如彼此通信的发动机控制模块(ECM)14、内燃发动机(ICE)16及车辆系统控制器(VSC)18。VSC18接收对应于车辆速度(Vx)、加速度(a)、横摆角(r)及扭矩(Tpwt,Tbrk)的输入,然后响应于所述输入而提供对应于道路坡度估计(R6E)和估计的车辆质量(Mv)的输出。
示出的实施例将车辆12描述为微混合动力车辆,微混合动力车辆是一种由发动机16推进的车辆,且发动机16反复地起动和停止以节省燃料。增强的起动电动机20结合到发动机的曲轴。起动电动机20接收电功率并将输出扭矩提供给曲轴以用于起动发动机16。虽然在上下文中示出和描述了微混合动力车辆12,但是应该理解到,车辆系统10的实施例可在其他类型的车辆(例如,传统的动力传动系车辆、由电机驱动而不存在发动机辅助的电池电动车辆(BEV)以及由电动机和发动机驱动的混合动力电动车辆(HEV))上实施。
车辆12包括变速器22,以用于调节发动机16的输出扭矩。来自发动机16的扭矩通过变速器22经变速器输出轴26传递到差速器24。半轴28从差速器24延伸到一对驱动轮30,以提供驱动扭矩来推进车辆12。
车辆12包括换档杆32,以用于手动选择变速器的档位。换档杆32包括传感器(未示出),以用于提供对应于选择的变速器档位(例如,PRNDL)的输出信号。变速器控制模块(TCM)34与换档杆32和变速器22通信,以用于基于换档杆的选择来调节变速器的速比。可选地,换档杆32可机械地连接到变速器22,以用于调节变速器的速比。
车辆12包括制动系统,制动系统包括制动踏板36、增压器和主缸,在图1中,增压器和主缸总体上由制动激活模块38指示。制动系统还包括ABS制动控制模块40,ABS制动控制模块40通过一系列液压管路44连接到车轮制动组件42和制动激活模块38,以实施摩擦制动。车轮制动组件42布置在每个车轮30处,且可被构造为卡钳或鼓式制动组件。每个车轮制动组件42包括传感器(未示出),以用于提供车轮速度信号(Ws)。在一个或多个实施例中,ABS制动控制模块40基于车轮速度信号将计算的车辆速度信号(Vx)提供给VSC18。
制动系统还包括传感器,以用于提供对应于当前制动特性的信息。制动系统包括位置开关,以用于提供对应于制动踏板位置(例如,踩下或松开)的制动踏板状态(Sbp)信号。在其他实施例中,制动系统包括位置传感器(未示出),以用于测量踏板位置。制动系统还包括一个或多个压力传感器,以用于提供对应于制动系统内实际制动压力值(例如,制动管路压力、主缸压力或车轮组件处的压力)的制动压力(Pbrk)信号。
车辆12包括加速踏板48,加速踏板48带有位置传感器,以用于提供对应于驾驶员请求推进的加速踏板位置(APP)信号。ECM14基于APP信号控制发动机16的节气门。
车辆12包括能量储存装置,例如电池50。电池50将电能供应到车辆控制器和起动电动机20以及用于变速器22的换档的电动辅助泵46,如总体上由图1中的虚线指示的。车辆12可包括单个电池50(例如,传统的低电压电池)或包括包含高电压电池的多个电池。另外,车辆12可包括其他类型的能量储存装置,例如电容器或燃料电池。
车辆12还包括一个或多个惯性传感器52,惯性传感器52提供对应于车辆加速度(a)的输出。在一个实施例中,惯性传感器52提供对应于纵向加速度的输出。在另一实施例中,惯性传感器52提供对应于纵向加速度、横向加速度和垂向加速度的输出,在图1中,所有这样的加速度总体上由变量a指示。在一个或多个实施例中,惯性传感器52还提供对应于横摆角(r)的输出。在其他实施例中,传感器52是倾角仪(未示出),该倾角仪提供对应于与车辆俯仰角相关的道路坡度的输出。
车辆12包括导航系统54,以用于从驾驶员接收目的地信息。导航系统54还与在车辆外部的多个系统/网络通信。例如,导航系统54可包括用于(通过蜂窝网络、硬线连接或射频(RF)波)(未示出)与卫星(例如,GPS)通信的收发器和计算机。这样的网络在由P1eet等人的第12/744,092号美国申请中详细描述,该申请通过引用被包含于此。导航系统54确定车辆12的当前位置和由驾驶员选择的目的地之间的距离,并将该距离和与车辆12相关的速度、加速度以及海拔信息、车辆12的目标目的地或其他相关的GPS路点一起提供给VSC18。
VSC18与其他车辆系统、传感器和控制器通信,以用于协调它们的功能。如在示出的实施例中所示,VSC18从各种车辆传感器和控制器接收多个输入信号(例如,Sbp、Pbrk、发动机速度(Ne)、车辆速度(Vx)、动力传动系扭矩(Tpt)、制动扭矩(Tbr)、a和r)。虽然VSC18被示出为单个控制器,但是VSC18可包括多个控制器,所述多个控制器可用于根据总的车辆控制逻辑或软件控制多个车辆系统。包括VSC18的车辆控制器通常包括任意数量的微处理器、ASIC、IC、存储器(例如,FLASH、ROM、RAM、EPROM和/或EEPROM)及软件代码,它们彼此共同作用以执行一系列操作。控制器还包括预定数据或“查找表”,所述“查找表”基于计算和测试数据并存储在存储器内。VSC18利用通用总线协议(例如,CAN和LIN)通过一个或多个有线车辆连接或无线车辆连接与其他车辆系统和控制器(例如,ECM14、TCM34等)通信。
参照图2和图3,车辆12被示出为多个力、加速度和力矩作用在车辆12上。车辆12被示出为参考以车辆的重心(CG)为中心的坐标系。坐标系包括纵轴(x)、横轴(y)和竖轴(z)。参照图2,车辆12位于斜面上。斜面的坡度(道路坡度)由符号αr表示。另外,车辆相对于道路坡度的俯仰角由符号0表示。车辆的重量由Mg表示,Mg被示出为作用在车辆12的CG上,其中,M是估计的车辆质量,g表示由于重力而导致的加速度(9.8m/s2)。车辆12的纵向速度由符号Vx表示,且根据一个或多个实施例,Vx基于车轮速度信号(在图1中示出的Ws)确定。车辆12的纵向加速度和垂向加速度分别由变量和表示,其中,上标字母(s)指示由传感器(例如,在图1中示出的惯性传感器52)提供的信号。参照图3,横摆角围绕竖轴(z)作用且由字母r表示。惯性传感器52安装在车辆12上从CG偏移一定距离的位置,车辆系统10基于该偏移量计算CG处的加速度值。在图3中,传感器52和CG之间的纵向偏移量由ld表示,传感器52和CG之间的横向偏移量由lc表示。
图4描述了根据一个或多个实施例的用于道路坡度估计仲裁的方法410的简化框图。根据一个或多个实施例,使用包含在VSC内的软件代码实现该方法。在其他实施例中,在其他车辆控制器中实现方法410,或者方法410分布在多个车辆控制器中。车辆系统包括用于估计道路坡度的多个算法,其中,每个算法适合于不同的车辆状况。方法410基于当前车辆状况选择道路坡度估计算法。
在信号处理框412处,VSC接收输入。所述输入包括车辆状态和传感器信号(例如,r)、获得的或计算的值(例如,Vx、Tpwt、Tbrk、Fpwt、Fbrk)及品质因子,所有这样的输入总体上由输入414表示。在一些实施例中,车辆系统从导航系统接收对应于车辆速度和加速度的输入,所述输入总体上由输入418表示。
在框420处,VSC使用分别由子框422、424和426表示的静态估计算法、运动学估计算法和动态估计算法来估计道路坡度(RGE)。RGE值被提供给RGE仲裁和诊断框428,RGE仲裁和诊断框428基于多个车辆状况选择一个估计(RGEst、RGEkin或RGEdyn)作为输出(RGEout)提供给其他车辆系统。
静态RGE算法422基于纵向加速度输入估计道路坡度。静态RGE适合于车辆静止状况,例如,当微混合动力车辆在斜面上停止时的发动机起动/停止功能。当车辆在斜面上停止(如图2所示)时,由于车辆悬挂系统内的挠度(该挠度由车辆俯仰角(0)表示),导致纵向加速度值存在偏差。静态RGE算法对车辆俯仰角进行补偿,如参照图7至图10C详细描述的。一般来说,静态RGE是三种RGE估计算法中最精确的算法,其原因是静态RGE算法的估计涉及最少估计输入且静态RGE算法最直接地获得估计结果。然而,静态RGE算法仅在车辆静止状况下可用。
运动学RGE算法424基于对应于车辆的运动学特性(包括车辆速度、加速度和横摆角之间的关系)的输入来估计RGE。该算法使用卡尔曼滤波器处理所述输入,如参照图11至图13D详细描述的。运动学RGE算法424适合于正常车辆运动状况(例如,Vx>5kph),其原因是在车辆速度低时一些基本计算的精度低。另外,卡尔曼滤波器是动态滤波器,通常不可能获得良好的估计,直到估计状态收敛为止。瞬态时间取决于初始状态、真实道路坡度状态以及估计状况(例如:信号频率分量丰富度和信噪比等)之间的差异。
动态RGE算法426基于对应于车辆速度、车轮扭矩(制动扭矩和动力传动系输出扭矩)以及作用在车辆上的其他拖曳力(例如,气动阻力、滚动阻力和道路坡度负荷)的输入来估计RGE。车辆系统包括用于估计车辆质量的算法,该算法包括在框426中且参照图18至图24详细描述。一般来说,与静态RGE和运动学RGE相比,动态RGE的品质状况更低,其原因是动态RGE在进行估计时涉及更多的信号,且每个信号中的噪声和误差将传送到最终RGE输出中。然而,动态RGE不依赖纵向加速度信号因此,动态RGE的信号覆盖范围不同于静态估计和运动学估计的信号覆盖范围。例如,在传感器出故障的情况下,可能不能获得纵向加速度信号然后,动态RGE是唯一可用的估计。
在框430处,车辆系统执行驱动状态监测和估计算法调节。估计算法调节包括对框424的RGEkin算法进行前馈输入补偿,这参照图14至图17详细描述。
在框432处,车辆系统评价车辆操作状况(例如,车辆稳定性、振荡和速度),并将对应的判断提供给框428。在框434处,车辆系统针对每种RGE算法确定品质因子(QFst、QFkin、QFdyn)以及总的RGE品质因子(QFRGE)。在一个或多个实施例中,QFRGE与RGEout一起提供给其他车辆系统。参照图25至图27G详细描述这样的品质评价。一般来说,品质因子评价估计的信息的可用性和精度。品质因子从0至3的范围内选择,其中,品质值3指示全品质,品质值2指示品质存在一些降低,品质值1指示品质差,品质值0指示估计在当前不可用。
图5A至图5E示出了图4的用于仲裁道路坡度估计的方法410的效果。图5A至图5E包括在实际道路坡度已知的道路上在共同的时间段内获得的测量数据的5个图形。图5A示出了车辆速度(Vx)和静态道路坡度估计(RGEst)。图5B示出了车辆速度(Vx)和运动学道路坡度估计(RGEkin)。图5C示出了车辆速度(Vx)和道路坡度估计输出(RGEout)。图5D示出了道路坡度估计输出(RGEout)和实际道路坡度(αr)。图5E示出了RGEout和αr之间的百分比误差。
图5A至图5C示出了在车辆多次停止的驱动循环期间,方法410如何在静态道路坡度估计和运动学道路坡度估计之间进行仲裁。在图5A至图5C中的每个图中示出了车辆速度(Vx)。在示出的实施例中,车辆停止三次,这三次停止由标号512、514和516指示。
如图5A描述的,当车辆停止或者静止时,静态道路坡度估计(RGEst)可用。在第一次停止512,RGEst指示道路坡度约为20%,如标号518表示的。在第二次停止514,RGEst指示道路坡度约为0%,如标号520表示的,以及在第三次停止516,RGEst指示道路坡度约为-12%,如标号522表示的。根据一个或多个实施例,当车辆运动时,静态道路坡度估计不可用,RGEst波形等于零。
如图5B描述的,当车辆运动时或者在正常车辆运动期间,运动学道路坡度估计(RGEkin)可用。在第一次停止512之前,RGEkin指示道路坡度从约0%增加到20%,如标号524表示的。在第一次停止512,RGEkin不可用,RGEkin锁止/锁定到最后那个RGEkin值,如标号526表示的。在第一次停止512和第二次停止514之间,RGEkin从约20%减小到0%,如标号528表示的。在第二次停止514和第三次停止516之间,RGEkin从约0%减小到-12%,如标号530表示的。在第三次停止516之后,RGEkin从-12%增加到0%,如标号532表示的。
如图5C描述的,如参照图4描述的道路坡度估计输出(RGEout)通常在静止状况下选择RGEst以及在车辆运动时选择RGEkin。图5D示出了RGEout和实际道路坡度(αr)。图5E描述了RGEout和αr之间的百分比误差(误差)。对于本示例,误差约为+/-4%。
表A描述了与用于估计道路坡度的算法相关的输入和车辆状况的比较,如下所示:
表A的第1列示出了三种道路坡度估计算法:静态、运动学和动态。第2列至第4列描述了每种算法所使用以进行估计的原始输入。例如,车辆系统基于车辆加速度估计静态道路坡度。车辆系统基于车辆速度和加速度估计运动学道路坡度,车辆系统基于车辆速度和车轮扭矩/力估计动态道路坡度。如参照图1描述的,这些输入可以是测量值或者获得的值。另外,在一些实施例中,对于对应的道路坡度估计,车辆系统除了使用在图5中描述的输入之外还使用其他输入。
第5列至第8列描述了每种道路坡度估计算法可用的车辆状况。例如,当车辆静止时,静态道路坡度估计算法可用。当车辆正常运动(例如,车辆速度大于5mph)时,运动学道路坡度估计算法可用。当车辆静止时或者在正常车辆运动期间,动态道路坡度估计算法可用。虽然表1的第6列指示在车辆速度极低(例如,小于5mph)时没有估计算法可用,但是车辆系统仍然在这样的状况期间提供精确的道路坡度估计,如下面参照图25至图26C描述的。
第7列描述了每种估计算法的收敛阶段。收敛阶段涉及在车辆状况之间过渡(例如,从静止过渡到速度极低,过渡到正常运动)期间车辆系统使用每种算法估计道路坡度所花费的时间的量。
参照图6,示出了根据一个或多个实施例的用于仲裁道路坡度估计的方法,该方法总体上由标号610指示。根据一个或多个实施例,使用包含在VSC内的软件代码实现该方法。在其他实施例中,在其他车辆控制器中实现方法610,或者方法610分布在多个车辆控制器中。在一个或多个实施例中,方法610的操作对应于参照图4描述的方法410的操作。车辆系统包括用于估计道路坡度的多种算法,其中,每种算法适合于不同的车辆状况。方法610基于估计状态、估计准备和估计品质因子选择道路坡度估计算法。
用于每种RGE算法(静态、运动学和动态)的估计状态因子被限定为布尔变量(“真”/“假”),且由变量(STATUSst,STATUSkin,STATUSdyn)表示。估计状态因子指示是否满足估计条件以及估计功能是否在对应的RGE通道(RGEst、RGEkin和RGEdyn)更新其状态。状态值“真”指示估计方法主动监测车辆动态状态(例如,和r),且正在更新道路坡度信息及其估计通道。
用于每种RGE算法(静态、运动学和动态)的估计准备因子也被限定为布尔变量(“真”/“假”),且由变量(READYst、READYkin,READYdyn)表示。估计准备通常指示RGE功能的当前有用性。更具体地说,静态估计准备(READYst)指示RGEst是否已经针对于当前车辆位置提供道路坡度估计。运动学估计准备值(READYkin)和动态估计准备值(READYdyn)分别指示在品质因子大于2(降低的品质或全品质)的情况下,它们的估计状态在超过阈值时间段的时间内持续有效。在存在不合格的估计状况时,在这样的状况下累积的时间未超过另一时间阈值,使得估计状态未偏离降低的边界。
用于RGE输出和每种RGE算法(静态、运动学和动态)的估计品质因子对应于0至3之间的值,且由变量(QFRGE、QFst、QFkin、QFdyn)表示。估计品质因子评价估计的道路坡度信息的可用性和精度。QFRGE是普通RGE输出信号品质评价,其基于选择的RGE算法的品质因子,且在更新事件之后根据车辆的动态状况被进一步处理。参照图25至图26C详细描述这样的品质评价。对于来自各个通道的评价的品质因子QFst、QFkin和QFdyn,值3指示全品质,值2指示估计仍然可用,但是当前值可能在指定的精度范围之外。
在操作612处,车辆系统评价静态RGE功能的当前估计状态,以确定STATUSst是否为等于“真”。如果在操作612处确定是肯定的,则车辆系统进行到操作614。在操作614处,车辆系统评价静态RGE功能的当前估计准备,以确定READYst是否等于“真”。如果在操作614处确定是肯定的,则车辆系统进行到操作616。在操作616处,车辆系统评价静态RGE功能的当前估计品质,以确定QFst是否等于3(全品质)。如果在操作616处确定是肯定的,则车辆系统进行到操作618,并通过将RGE输出设置为等于RGEst来选择静态道路坡度估计。如果在操作612处确定是否定的,则车辆系统进行到操作620。
在操作620处,车辆系统评价运动学RGE功能的当前估计状态,以确定STATUSkin是否等于“真”。如果在操作620处确定是肯定的,则车辆系统进行到操作622。在操作622处,车辆系统评价运动学RGE功能的当前估计准备,以确定READYkin是否等于“真”。如果在操作622处确定是肯定的,则车辆系统进行到操作624。在操作624处,车辆系统评价运动学RGE功能的当前估计品质,以确定QFkin是否等于3(全品质)。如果在操作624处确定是肯定的,则车辆系统进行到操作626,并通过将RGE输出设置为等于RGEkin来选择运动学道路坡度估计。
静态道路坡度估计和运动学道路坡度估计优于RGEdyn,其原因是:由于在RGEdyn进行估计时需要附加信号导致RGEdyn的品质状况低,且每个信号中的噪声将传送到最终RGE输出中。因此,即使静态品质因子和运动学品质因子小于3,车辆系统仍然可选择RGEst或RGEkin,而非选择RGEdyn。
如果在操作616处确定是否定的,则QFst小于3,车辆系统进行到操作628。在操作628处,车辆系统将临时品质因子值(QFtmp)设置为等于QFst,以及将临时道路坡度估计(RGEtmp)设置为等于RGEst。类似地,如果在操作624处确定是否定的,则QFkin小于3,车辆系统进行到操作630。在操作630处,车辆系统将临时品质因子值(QFtmp)设置为等于QFkin,以及将临时道路坡度估计(RGEtmp)设置为等于RGEkin。在操作628或630之后,车辆系统进行到操作632,以评价总的RGE品质因子(QFRGE)。
在操作632处,车辆系统确定QFRGE是否小于3。如果在操作632处确定是否定的,则车辆系统进行到操作634,不更新RGE。如果总的RGE具有全品质,则车辆系统不去选择品质更低的不同的RGE。如果在操作632处确定是肯定的,则车辆系统进行到操作636。
在操作636处,车辆系统评价动态RGE功能的当前估计状态和当前估计准备,以确定STATUSdyn和READYdyn是否等于“真”。如果在操作636处确定是否定的,则车辆系统进行到操作638。在操作638处,车辆系统评价临时品质因子QFtmp,以确定QFtmp是否等于或大于2。如果在操作638处确定是否定的,则车辆系统进行到操作634,不更新RGE。如果在操作638处确定是肯定的,则车辆系统进行到操作640,并通过将RGE输出(RGEout)设置为等于RGEtmp来选择临时道路坡度(静态或运动学)。如果在操作636处确定是肯定的,则车辆系统进行到操作642。
在操作642处,车辆系统评价临时品质因子(QFtmp)和动态品质因子(QFdyn),以确定QFtmp和QFdyn中的最小值是否等于或大于2。如果在操作642处确定是肯定的,则车辆系统进行到操作644。在操作644处,车辆系统确定QFtmp是否大于QFdyn。如果在操作644处确定是肯定的,则车辆系统进行到操作640,并通过将RGE输出(RGEout)设置为等于RGEtmp来选择临时道路坡度(静态或运动学)。如果在操作644处确定是否定的(QFtmp不大于QFdyn),则车辆系统进行到操作646,并通过将RGE输出(RGEout)设置为等于RGEdyn来选择动态道路坡度。
如果在操作614、620或622处确定是否定的,则车辆系统进行到操作648。在操作648处,车辆系统评价动态RGE功能的当前估计状态和当前估计准备,以确定STATUSdyn和READYdyn是否均等于“真”。
如果在操作648处确定是肯定的,则车辆系统进行到操作650。在操作650处,车辆系统评价动态RGE功能的当前估计品质,以确定QFdyn是否等于或大于2。如果在操作650处确定是肯定的,则车辆系统进行到操作646,并通过将RGE输出设置为等于RGEdyn来选择动态道路坡度估计。如果在操作648或650处确定是否定的,则车辆系统进行到操作652,并通过将RGE输出值锁止到之前的RGE输出值而不更新RGE输出。在操作618、626、640、646或652之后,车辆系统进行到操作654。在操作654处,车辆系统提供RGE输出和QF,然后返回到操作612,以进行方法610的另一次迭代。
这样,通过从多个不同的估计(静态、运动学和动态)中仲裁或选择道路坡度估计,使车辆系统优于现有的系统。每个道路坡度估计适合于不同的车辆状况。通过在这些不同的估计之间进行仲裁,车辆系统在宽的车辆状况范围内提供道路坡度估计。
回顾图2,静态道路坡度估计和运动学道路坡度估计基于车辆运动学方程。运动学通常涉及物体在空间和时间中如何运动的描述。车辆运动学方程基于由惯性传感器52提供的车辆加速度值坐标系(x,y,z)相对于道路且围绕车辆的重心(CG)定位。然而,惯性传感器52被安装为从坐标系(x,y,z)偏移一定距离,其中,lc是横向偏移量以及ld是纵向偏移量。另外,惯性传感器52提供由于车辆悬挂导致的从坐标系偏移的输出信号。车辆俯仰角(θ)表示从纵轴(x)偏移的角度。车辆侧倾角(φ)(未示出)表示从横轴(y)偏移的角度,车辆横摆角(r)表示从竖轴(z)偏移的角度。车辆系统估计道路坡度(αr),因此从道路坡度(αr)分离车辆偏移量(例如,车辆俯仰角(θ))。
起主导作用的车架运动学方程如下所示:
其中,和是车辆纵向加速度、横向加速度和垂向加速度。lc是在转弯机动时传感器从车辆旋转中心(CG)横向偏移的量,ld是传感器纵向偏移的量。Vx和Vy是纵向速度和横向速度,是计算的纵向加速度、横向加速度和垂向加速度;θ是车辆俯仰角,r是横摆角。φ是总的车辆侧倾角,其对应于车身侧倾角φv和道路坡度角φb之和。另外,naz表示由车辆悬挂振荡和路面起伏导致的所有的垂向加速度扰动。
车辆系统基于车辆纵向加速度和等式1估计静态道路坡度。由于车辆是静态的或者车辆静止,所以纵向速度(Vx)、纵向加速度横摆角、横摆角的倒数被设置为零,且简化满足运动学方程。等式4示出了用于在静止状况下计算纵向加速度的等式,如下所示:
由于悬挂的挠度导致的车辆俯仰角(θ)与道路坡度角成比例θ=Kp(gsinαr),其中,Kp表示在稳定状态下车身的刚度。Kp是基于的幅度的预定值。在一个或多个实施例中,VSC包括与预校准数据的查找表对应的存储数据,其中,在的正操作范围和负操作范围内,针对于的矢量值对Kp值进行插值。等式5示出了中间推导步骤,等式6提供用于在静止状况下计算道路坡度(αr)的等式,如下所示:
一般来说,静态RGE是三种RGE估计算法中最精确的算法,其原因是静态RGE涉及最少估计输入且静态RGE最直接地获得估计结果。然而,静态RGE仅在车辆静止状况下可用。
图7示出了根据一个或多个实施例的当车辆在坡上停止时获得的测量数据。图7示出了输入信号:车速(Vx)和纵向加速度以及根据等式6计算的道路坡度数据(αr)。在时间t1,车辆在减速,如由Vx的下降斜坡指示的。在时间t2,车辆停止。在车辆停止之后,由于车身振荡导致在纵向加速度信号上存在噪声,如由标号712指示的。通过道路坡度计算放大该噪声,如由道路坡度αr波形示出的且由标号714指示的。在时间t3,纵向加速度信号已经足够收敛,而提供道路坡度估计的品质信息,如由标号716指示的。在示出的实施例中,t2和t3之间的时间差约为3s。因此,对于RGE,依赖于道路坡度(αr)计算的车辆系统在使用该值之前可等待大约3s。这样的延时导致使用发动机起动/停止功能的微混合动力车辆的燃料效率降低。
参照图8,示出了根据一个或多个实施例的用于评价加速度计收敛状态的方法,该方法总体上由标号810指示。根据一个或多个实施例,使用包含在VSC内的软件代码实现该方法。虽然车辆系统估计静态道路坡度,但是可连续监测其估计状态和品质。当加速度计输出收敛且保持这样的稳定状态持续了最小量的时间时,发生高品质的估计。加速度计输出被分配收敛状态因子(CONst),CONst在三级标签代码中指示原始信号是否处于稳定状态。其中,CONst等于零指示收敛不稳定,CONst等于1指示收敛在某种程度上稳定,CONst等于2指示收敛极其稳定。
另外,收敛状态信号(CONst)还用于评价静态RGE输出品质。例如,在一个或多个实施例中,CONst为2(极其稳定)对应于最高的或全品质估计(QFst=3)。CONst为1(在某种程度上稳定)对应于品质劣化的估计(QFst=2)。CONst为0(不稳定)对应于不合格的估计(QFst=1)或者无数据(QFst=0)。收敛状态信号(CONst)基于纵向加速度的导数。
在操作812处,车辆系统评价车辆速度(Vx)和收敛状态因子(CONst),以确定Vx是否小于低速阈值以及加速度计的信号品质是否良好(例如,是否是CONst=2)。如果在操作812处确定是否定的,则车辆系统进行到操作814。在操作814处,车辆系统评价收敛状态因子(CONst),以确定CONst是否指示信号稳定(例如,是否是CONst=2)。如果在操作814处确定是否定的,则车辆系统进行到操作816,通过将稳定计时器(Stable_Timer)、收敛计时器(Converging_Timer)、不稳定计时器(Nonstable_Timer)以及收敛状态因子(CONst)均设置为零(CONst=不稳定),而使计时器和标志复位。当存在可用的RGE值时,车辆系统还将标志Grade_Found设置为等于“假”。Grade_Found指示对于当前车辆停止事件是否已经提供有效的估计(RGEst)。Grade_Found等于“真”,指示已经提供有效的RGEst,Grade_Found等于“假”,指示没有提供有效的RGEst或者RGEst不再有效。如果在操作812处确定是肯定的,则车辆系统进行到操作818。
Grade_Found意思是在当前车辆停止事件中在违背估计条件之前,从静态估计算法获得有效的坡度估计。“违背估计条件”意思是对于估计来说状况不再有效。对于静态RGE,这样的一个事件是车辆速度增加到速度阈值(例如,5kph)之E。
在操作818处,车辆系统计算纵向加速度的时间导数或者使用低通滤波器计算纵向加速度的一阶导数。在操作820处,车辆系统将该加速度的时间导数与收敛阈值进行比较,以确定是否小于收敛阈值。如果在操作820处确定是否定的,则车辆系统进行到操作822,使不稳定计时器累加循环时间Dt。如果在操作814处确定是肯定的,则车辆系统也进行到操作822。在操作824处,车辆系统评价不稳定计时器,以确定不稳定计时器是否大于不稳定计时器阈值。如果在操作824处确定是肯定的,则车辆系统进行到操作826,将收敛计时器、稳定计时器和收敛状态因子(CONst)均设置为零(CONst=不稳定)。如果在操作820处确定是肯定的,则车辆系统进行到操作828。
在操作828处,车辆系统评价纵向加速度的时间导数以确定是否小于稳定阈值。如果在操作828处确定是否定的,则车辆系统进行到操作830,使收敛计时器累加循环时间Dt。在操作832处,车辆系统评价收敛计时器,以确定收敛计时器是否大于收敛计时器阈值。如果在操作832处确定是肯定的,则车辆系统进行到操作834,通过将不稳定计时器和稳定计时器设置为等于零而使不稳定计时器和稳定计时器复位;将收敛状态因子设置为在某种程度上稳定(CONst=1);将Grade_Found设置为“真”。如果在操作828处确定是肯定的,则车辆系统进行到操作836。
在操作836处,使稳定计时器累加循环时间Dt。在操作838处,车辆系统评价稳定计时器,以确定稳定计时器是否大于稳定计时器阈值。如果在操作838处确定是肯定的,则车辆系统进行到操作840,通过将不稳定计时器和收敛计时器设置为等于零而使不稳定计时器和收敛计时器复位;将收敛状态因子设置为极其稳定(CONst=2);将Grade_Found设置为“真”。
参照图9,示出了根据一个或多个实施例的用于锁定静态道路坡度输出的方法,该方法总体上由标号910指示。根据一个或多个实施例,使用包含在VSC内的软件代码实现该方法。一般来说,一旦车辆系统确定收敛状态极其稳定(例如,CONst=2),则将RGEst锁定以避免从未来车辆扰动导致的任何冲击。当车辆系统检测到车辆在运动时,锁定被释放。STATIClock信号指示RGEst是否被锁定、解除锁定或者隔离。被隔离的锁定状态指示锁定和解除锁定之间的中间状态,此时,车辆系统等待附加信息。
在操作912处,车辆系统进入初始化状态,将STATIClock信号设置为等于解除锁定。在操作914处,车辆系统评价STATIClock信号,以确定静态道路坡度(RGEst)是否被锁定。如果在操作914处确定是否定的,则车辆系统进行到操作916。在操作916处,车辆系统确定STATIClock信号是否被隔离。如果在操作916处确定是肯定的,则车辆系统进行到操作918。在操作918处,车辆系统确定Grade_Found是否为“真”。Grade_Found信号基于当前RGEout信号。如果在操作918处确定是肯定的,则车辆系统进行到操作920。在操作920处,车辆系统评价收敛状态信号,以确定CONst是否“极其”稳定(例如,是否是CONst=2)。如果在操作920处确定是肯定的,则车辆系统进行到操作922,将STATIClock信号设置为锁定。然后,在操作924处,静态道路坡度估计(RGEst)被锁定。如果在操作914处确定是肯定的(例如,STATIClock=锁定),则车辆系统进行到操作926。
在操作926处,车辆系统评价获得的坡度(GRADE_Found),以确定是否是GRADE_Found=“真”。如果在操作926处确定是肯定的,则车辆系统进行到操作928,评价收敛状态信号以确定CONst是否不稳定(例如,是否是CONst=0)。如果在操作928处确定是肯定的,则车辆系统进行到操作930,将STATIClock设置为隔离。在操作930之后,并且如果在操作928或934处确定是否定的,则车辆系统进行到操作924,并锁定静态道路坡度估计(RGEst)。如果在操作916处确定是肯定的,则车辆系统进行到操作932。
在操作932处,车辆系统评价获得的坡度信号,以确定是否是GRADE_Found=“真”。如果在操作932处确定是肯定的,则车辆系统进行到操作934,评价收敛状态信号,以确定CONst是否极其稳定(例如,是否是CONst=2)。如果在操作934处确定是肯定的,则车辆系统进行到操作922,将STATIClock设置为锁定。如果在操作932、918、920或926处确定是否定的,则车辆系统进行到操作936,将STATIClock设置为解除锁定。在操作936或924之后,车辆系统进行到操作938,然后返回到操作914。
参照图10A至图10C,车辆系统基于图8的方法810和图9的方法910估计静态道路坡度(RGEst)。图10A至图10C包括在共同的时间段内获得的测量数据的3个图形,这3个图形总体上由标号1010指示。另外,根据一个或多个实施例,车辆系统基于运动学道路坡度估计或动态道路坡度估计来估计RGEst的初始值,然后使用带宽可变的低通滤波器对纵向加速度信号进行滤波,以确定RGEst。
图10A示出了车辆速度(Vx)和表示低通滤波器的可变带宽的信号(flt),其中,当flt等于1时,滤波器具有窄的带宽,当flt等于零时,滤波器具有宽的带宽。图10B示出了测量的纵向加速度值图10C示出了静态道路坡度估计(RGEst)和运动学道路坡度估计(RGEkin)。在一个或多个实施例中,沿着在0.04Hz和2.5Hz之间的频率范围分配带宽;其中,宽的带宽对应于频率范围的高端,窄的带宽对应于频率范围的低端。在其他实施例中,宽的带宽对应于大约1Hz,窄的带宽对应于大约0.05Hz。
车辆系统基于运动学道路坡度估计或动态道路坡度估计来估计RGEst的初始值。如果没有估计可用,则车辆系统将使用默认值。由于静态RGE功能仅在车辆速度小于预定速度阈值(例如,1kph)之后开始启用,所以估计算法的初始值是收敛时间或延时的因子,直到RGEst收敛到实际道路坡度为止。幸运的是,车辆系统包括多种道路坡度算法。其他运动学估计算法和动态估计算法主要在该状况下同时车辆在运动时起作用。因此,在静态RGE算法开始启用之前,关于在当前时间之前的几秒钟内的道路坡度的信息或者关于在车辆(假设车辆向前行驶)之后的几米内的道路坡度的信息应该可从其他RGE算法获得,这些信息是静态RGE算法开始启用的有效参考。在此之后,静态RGE算法接管了主要的估计操作,以继续监测道路坡度变化,同时在不存在合格的估计条件时暂停其他两种算法。更具体地说,当满足估计条件时,根据运动学估计算法的输出和动态估计算法的输出中的哪一个具有最高品质评价,静态RGE算法将运动学估计算法的输出或动态估计算法的输出作为其初始值。当运动学估计算法的输出和动态估计算法的输出均不满足品质时,将使用默认值替代。合格的静态估计条件是车辆速度将小于速度阈值,这样的速度阈值将小于用于确定其他两种RGE算法暂停的速度阈值。
在时间t1,车辆开始减速,如由Vx的下降斜坡指示的。在时间t2,车辆停止,RGEst启用。与从0%开始并收敛到实际道路坡度不同的是,RGEst从最后那个RGEkin值开始,如在点1012处所示,在点1012处,RGEst等于RGEkin,二者均约等于-22%。然后,随着RGEst累积数据,在时间t3,RGEst收敛到约-20%的实际道路坡度,如由标号1014表示的。
在车辆停止之后,由于车身振荡导致在纵向加速度信号上存在噪声,如由标号1016指示的。车辆系统使用带宽可变的低通滤波器对纵向加速度信号进行滤波,如由1018处的flt波形指示的。车辆系统通过使用基于的一阶导数和更高阶导数选择的时间常数对数据进行动态滤波,而计算滤波后的纵向加速度值。在一个实施例中,车辆系统将纵向加速度的时间导数或纵向加速度的一阶导数与阈值时间导数进行比较。如果该加速度的时间导数大于阈值时间导数,则为低通滤波器选择对应于更低带宽的时间常数。相反,如果该加速度的时间导数小于阈值时间导数,则为低通滤波器选择对应于更高带宽的时间常数。例如,阈值时间导数可以是在上面参照图8描述的收敛阈值820或稳定阈值828。在示出的实施例中,flt波形对应于带宽,其中,当应用低带宽时,flt等于1,当选择高带宽时,flt等于0。
在另一实施例中,基于纵向加速度的二阶导数从预定数据(例如,查找表)选择带宽。通过对进行滤波,道路坡度估计(RGEst)相对快速地收敛。例如,在时间t3,RGEst收敛,如由标号1014指示的,在t2和t3之间的延时约为1s。
回顾图7,在进行道路坡度估计之前,用于评价道路坡度的现有的系统可等待纵向加速度信号收敛。因此,对于在图7中示出的车辆状况,在估计道路坡度之前,这样的系统将等待大约3秒。然而,如图10A至图10C所示,通过将基于其他估计的初始值用于RGEst,并通过对纵向加速度进行动态滤波,车辆系统减少延时。
这样,通过使用运动学方程和纵向加速度的动态滤波来估计静态道路坡度,使车辆系统优于现有的系统。车辆系统还将静态道路坡度估计的初始值设置为等于由运动学道路坡度估计或动态道路坡度估计提供的道路坡度。在可提供合格的静态道路坡度估计值之前,通过对纵向加速度进行动态滤波,并通过针对初始值使用另一估计,车辆系统减少延时,这进而允许在起动/停止应用中发动机更早地关闭,由此使得燃料经济性提高。
参照图11,示出了根据一个或多个实施例的说明用于运动学道路坡度估计的方法的示意性框图,该框图总体上由标号1110指示。方法1110对应于在图4中示出的总的RGE仲裁方法410的运动学道路坡度估计部分。框1112对应于图4的框412。在框1112处,车辆系统接收输入,所述输入包括加速度横摆角(r)和车辆速度(Vx)。在一个实施例中,车辆系统仅接收纵向加速度和车辆速度。在其他实施例中,车辆系统还接收GPS信息(例如,车辆速度、加速度和海拔)。可从参照图1描述的传感器直接接收输入,或者可在CAN总线上间接接收输入。
框1114和1116对应于图4的RGEkin框424。在框1114处,车辆系统确定分别表示系统输出和输入的扩展卡尔曼滤波器(EKF)参数Yk和uk。输入参数基于车辆速度(Vx)、纵向加速度和加速度偏移量对于使用GPS信息的实施例,参数还可基于可选状态信息。在框1116处,车辆系统使用EKF估计RGEkin并确定运动学品质因子(QFkin)。
输入(和r)提供给对应于图4的框430的RGEkin输入补偿框1118。在框1118处,车辆系统确定加速度偏移量提供加速度偏移量作为框1114和1116的运动学道路坡度估计的前馈输入补偿值。参照图14至图17详细描述确定加速度偏移量。
车辆系统基于车架运动学方程1至4来估计运动学道路坡度,并使用积分式滤波策略。通过忽略更高阶次的角度项,车架运动学方程可简化为下面示出的等式7至9:
用于纵向加速度的等式可通过将小角度逼近应用于等式8而获得,如下所示:
其中,表示加速度偏移量。
根据一个或多个实施例,车辆系统使用卡尔曼滤波器来估计运动学道路坡度。卡尔曼滤波器对噪声输入数据流进行递归操作,以提供基本系统状态的统计上的优化估计。基本上,卡尔曼滤波器是这样一种算法,该算法使用在一定时间段内观察的且包含随机变化(噪声)的一系列测量值,提供对于未知变量的估计,这比只进行测量更精确。卡尔曼滤波器可被认为是基于可测量的输入和输出来估计“黑盒子”的状态的方法。
卡尔曼滤波器算法使用包括时间更新预测步骤和测量值更新步骤的两步处理来估计状态。在预测步骤中,卡尔曼滤波器提供当前状态变量及其不确定性的估计。测量值更新步骤是反馈控制的形式。一旦卡尔曼滤波器观察到下一个测量值,则使用加权平均更新估计,其中,更大的权重给予确定性更高的估计。
从协方差计算权重,协方差是对于系统状态的预测的估计不确定性的测量值。每个时间步骤重复卡尔曼滤波器处理,其中,新的估计及其协方差通知用于后续迭代的预测。卡尔曼滤波器是递归估计器。其意思是只有来自先前时间步骤的估计状态和当前测量值用于计算当前状态的估计。
对于在过渡系统和测量系统中存在附加的独立白噪声的线性系统模型来说,卡尔曼滤波器是优化的估计。然而,大多数系统是非线性的。扩展卡尔曼滤波器(EKF)用于提供非线性系统的估计。EKF采用包括多元泰勒级数展开的技术,以线性化关于工作点的计算。
车辆系统限定扩展状态由于重力加速度是恒定值(9.8m/s2),所以扩展状态的一阶动态变化等于零扩展状态使由道路坡度的变化引起的车辆加速度动态变化模型化,道路坡度在低频时相对大些。
用于估计的总体系统模型如下所示,且由等式11指示:
Xk=f(Xk-1,uk)+ωk
Yk=[10]Xk=Vx
ωk=[ωv(k)ωa(k)]T(11)
其中,Xk是在时间步骤(k)处的状态,Xk-1是在先前的时间步骤处的状态。uk是系统输入,Yk是系统输出。ωk是过程噪声,ωv(k)是与速度相关的过程噪声,ωa(k)是与加速度相关的过程噪声。Vx是纵向速度。是由于重力导致的扩展状态加速度,是测量的纵向加速度,是加速度偏移量。τs是用于上述离散时间实现模型的采样时间,T指示矩阵变换。
基于在等式11中示出的系统模型,观测仪被设计为估计扩展状态并且根据如下所示的等式12估计运动学道路坡度:
使用标准扩展卡尔曼滤波器(EKF)系统模型设计观测仪,该系统模型由在上面示出的等式11中的函数和等式13表示,这两个等式如下所示:
Xk=f(Xk-1,uk)+ωk
Yk=CXk+ek(13)
其中,C是矩阵,该矩阵在等式11中是相对简单的矩阵,ek是未知的噪声。
EKF用于非线性系统。非线性系统模型在每个时间步骤围绕当前状态被线性化,并提供转移矩阵Fk。然后,获得的转移矩阵Fk用于完成标准EKF递归的步骤。这些递归由两个更新步骤(预测步骤和测量值更新步骤)来描述。在预测步骤中,系统输入用于预测系统的未来状态如下面的等式14和15所示:
其中,Pk|k-1是估计的误差协方差,Qk是过程噪声协方差,其等于Fk是雅可比行列式,由下面示出的等式16提供:
测量值更新步骤由下面示出的等式17至19提供:
Kk=Pk|k-1CT(CPk|k-1CT+Rk)-1(17)
Pk|k=(I-KkC)Pk|k-1(19)
其中,Kk是EKF增益,Rk是测量值噪声协方差,其等于
一般来说,对于EKF系统,过程噪声协方差(Q)和测量值噪声协方差(R)不是已知的先验值或者与经验无关。然而,状况良好的信号(例如,Kk和Pk|k)将收敛到与Q和R相关的值。因此,Kk的收敛值可从预定数据获得,且该收敛值将直接用于实时估计。
在一个或多个实施例中,车辆系统在其估计运动学道路坡度时包括垂向位移。垂向位移(Z)可由图1的导航系统基于GPS信息来提供。在球形惯性系中起主导作用的车辆垂向运动学方程根据如下所示的等式20计算:
等式(11)的EKF系统模型被修改以包括垂向位移Z(k),如下面的等式21所示:
其中,表示基于GPS的加速度偏移量,等于在使用GPS信息的实施例中,提供作为前馈输入补偿。等式21与等式11的不同之处在于:等式21包括两种扩展状态:和其中,是由于车辆俯仰角导致的加速度。
参照图12,示出了根据一个或多个实施例的用于估计运动学道路坡度的方法,该方法总体上由标号1210指示。根据一个或多个实施例,使用包含在VSC内的软件代码实现该方法。在其他实施例中,在其他车辆控制器中实现方法1210,或者方法1210分布在多个车辆控制器中。
在操作1212处,车辆系统初始化估计器参数值。估计器参数值包括Xk、uk和Yk。通过将估计器参数值设置为等于默认值而初始化估计器参数值。在操作1214处,车辆系统确定运动学估计状态(STATUSkin)是否是开启。估计状态因子被限定为布尔变量(“真”/“假”),其中,“真”指示状态是开启,“假”指示状态是未开启。状态值“真”指示估计方法主动监测车辆动态状态(例如,和r),道路坡度信息及其估计通道在更新。如果在操作1214处确定是肯定的,则车辆系统进行到操作1216。
在操作1216处,车辆系统评价当前车速(Vx),以确定Vx是否等于或大于阈值车辆入口速度(vspd_entry)。在一个实施例中,阈值车辆入口速度约为7kph(4.3mph)。如果在操作1216处确定是否定的,则车辆系统进行到操作1218。在操作1218处,车辆系统评价Vx,以确定Vx是否小于阈值车辆出口速度(vspd_exit)。在一个实施例中,阈值车辆出口速度约为3kph(1.9mph)。如果在操作1218处确定是肯定的,则车辆系统进行到操作1220。
在操作1220处,车辆系统将运动学估计状态(STATUSkin)设置为暂停。通过将STATUSkin设置为暂停,车辆系统停止估计,并保留当前RGEkin值,而非使该值复位到默认值。在操作1222处,车辆系统评价当前速度(Vx),以确定Vx是否小于阈值车辆静止速度(vspd_stand)。阈值车辆静止速度对应于由表A的第6列指示的极低的车辆速度。在一个实施例中,阈值车辆静止速度约为1kph(0.6mph)。如果在操作1222处确定是肯定的,则车辆系统进行到操作1224,将运动学估计状态(STATUSkin)设置为关闭(“假”)。如果在操作1216处确定是肯定的,或者在操作1218处确定是否定的,则车辆系统进行到操作1226。
操作1226至1236涉及扩展卡尔曼滤波器操作,该操作在上面参照等式11至21描述过。首先,在操作1226处,车辆系统计算补偿的估计器输入(uk)。在操作1228处,车辆系统获得新的系统输出(Yk)。在操作1230处,车辆系统基于输出预测误差更新EKF状态(Xk),如在上面参照等式14至16描述的。在操作1232处,车辆系统更新EKF参数值(Xk、uk和Yk)。在操作1234处,车辆系统预测下一个系统动态输出Yk。然后,在操作1236处,车辆系统提供运动学道路坡度估计。
如果在操作1214处确定是否定的,则车辆系统进行到操作1238。在操作1238处,车辆系统评价当前车辆速度(Vx),以确定Vx是否等于或大于阈值车辆入口速度(vspd_entry),这类似于在操作1216处的确定。在一个实施例中,阈值车辆入口速度约为7kph(4.3mph)。如果在操作1238处确定是肯定的,则车辆系统进行到操作1240,将运动学估计状态(STATUSkin)设置为开启(“真”)。在操作1240之后,车辆系统进行到操作1226,以开始EKF操作。
在操作1224或1236之后,以及在操作1222或1238处确定是否定之后,车辆系统进行到操作1242。在操作1242之后,车辆系统返回到操作1214,以进行另一次迭代。
图13A至图13D是示出当车辆在斜面上停止时用于估计运动学道路坡度的方法1210的效果的图形。图13A至图13D包括在共同的时间段内获得的测量数据的4个图形,这4个图形总体上由标号1310指示。图13A示出了根据一个实施例的在主缸内测量的制动压力(Pbrk)。在其他实施例中,在制动管路内或者在车轮组件处测量Pbrk。图13B示出了车辆速度(Vx)和测量的纵向加速度值。图13C示出了档位选择(PRNDL),其中,值3对应于驱动(D)档位选择。图13C还示出了制动踏板的状态(Sbp),其中,正值指示施加制动,零值指示释放制动。图13D示出了静态道路坡度估计(RGEst)和运动学道路坡度估计(RGEkin)。
在时间t1,车辆在加速,如由Vx的上升斜坡指示的。在时间t2,车辆开始减速,如由Vx的下降斜坡指示的。在时间t3,车辆停止。在时间t4,车辆再次开始加速,在时间t5,车辆再次开始减速。运动学道路坡度估计(RGEkin)相对快速地收敛,如由标号1312和1314指示的。
如参照操作1222和1224描述的,如果车辆速度减小到车辆静止阈值速度(例如,1kph)之下,则运动学估计状态(STATUSkin)被设置为关闭。这样的条件由RGEkin示出且由标号1320指示。在t3之后,估计状态被设置为关闭,如在1320处在RGEkin上不存在任何可察觉的噪声所描述的。另外,一旦车辆停止,则静态道路坡度估计的状态(STATUSst)和RGEst的状态收敛,如由标号1322指示的。如在1320处所示,即使运动学状态关闭,RGEkin仍然在精确的道路坡度估计处锁定,RGEst和RGEkin均提供大致相等的坡度估计(17.56%)。
使用运动学方法估计道路坡度的现有的方法使用导数滤波策略。例如,一些现有的方法基于下面的车辆纵向运动学方程来估计道路坡度:
其中,是车辆纵向速度的导数,ax是车辆纵向加速度,r是横摆角,v是车辆横向速度。然而,由于导数计算导致这样的方法放大了信号上的任何噪声。
这样,通过使用EKF估计运动学道路坡度,提供积分式滤波策略,使车辆系统优于现有的方法。这样的积分策略不会放大信号上的噪声,因此,这样的积分策略提供改进的估计精度同时相对快速地收敛。
车辆系统使用对应于加速度偏移量的前馈输入补偿来估计运动学道路坡度。根据如下所示的等式23计算加速度偏移量:
其中,lc是横向偏移量,r是横摆角,Vy是横向速度,0是车辆俯仰角。在一个或多个实施例中,车辆系统不接收这些值(lc、Vy和0)作为输入。相反,车辆系统使用多种补偿策略从其他输入估计这些值。在一个或多个实施例中,车辆系统使用科里奥利效应(CoriolisEffect)补偿策略估计等式23的横向速度(Vy)项,使用传感器横向位置偏移效应补偿策略估计横向偏移量(lc)项,以及使用车身相对俯仰效应补偿策略估计车辆俯仰角(θ)。
车辆系统使用科里奥利效应补偿策略估计等式23的rVy项,该rVy项涉及车辆的转弯机动。图1的惯性传感器52提供横摆角,然而,通过车辆系统估计横向速度(Vy)。提供两种补偿方法来估计横向速度(Vy):第一种方法是基于估计的方法;第二种方法是基于导数的方法。分别基于横向动力等式8和车辆侧偏刚度参数提供这两种方法。此外,在连续变化的车辆驱动状况下对这两种方法进行积分,以提供估计的横向速度(Vy)。
用于估计横向速度(Vy)的第一种方法(基于估计的方法)通过将小角度逼近应用于等式8而获得,且第一种方法基于测量的横向加速度根据如下所示的等式24计算横向速度(Vy):
道路坡度角(φr)未知。因此,道路坡度角被处理为扰动且从计算中去除。车身侧倾角φv可由测量的横向加速度进行近似,如等式25所示:
其中,Kφ表示车身侧倾刚度。
等式24可简化为如下面的等式26所示那样:
通过对等式26的两侧进行积分,获得用于计算估计的横向速度的等式,如等式27所示:
测量的横向加速度通常受到加速度扰动和坡度角的影响。因此,由于偏向低频的加速度分量,使得在等式24至27中描述的Vy估计方法适合于短时间间隔的车辆动态机动,其中,存在幅度大的横向加速度和横摆角(r)。
与第一方法不同,用于估计横向速度(Vy)的第二方法(基于导数的方法)不依赖于积分。从车辆横摆角(r)到车辆侧滑角(β)的稳定状态增益(Gr2β)可从线性车辆模型获得,如等式28所示:
其中,Mr是分配到后轴的车辆质量,Cαr是集总的后轴侧偏刚度。
由于横向速度(Vy)大致等于纵向速度(Vx)与车辆侧滑角(β)的乘积,所以可获得用于计算横向速度的等式,如下所示:
其中,kU是侧偏刚度,等于Mr除以Cαr。
因此,横向速度可通过如下所示的等式30估计:
下面提供科里奥利效应补偿方程:
参照等式28至31描述的第二方法适合于侧滑角(β)小的稳定状态车辆转弯机动。
车辆系统实现综合的RGEkin架构,该架构包含(基于估计的)第一方法和(基于导数的)第二方法。然后,车辆系统基于由横摆角(r)体现的车辆动态状况而选择这两种方法中的哪一种方法来应用。通常的设计构思是:在主要是向前驱动时以及在稳定状态转动情况下使用来自等式30的基于导数的Vy,但是在存在大的车辆横向动力(速度和加速度)时切换到来自等式27的基于估计的Vy。由于不期望执行估计方法持续长的时间段,所以在车辆连续机动操作期间,每当可观察到稳定状态状况时,估计的Vy必须复位到起源状态。估计的Vy还可以预定的上限和下限为边界。通过低通滤波器提供输出Vy。
参照图14,示出了根据一个或多个实施例的用于确定车辆驱动状况的方法,该方法总体上由标号1410指示。根据一个或多个实施例,使用包含在VSC内的软件代码实现该方法。在其他实施例中,在其他车辆控制器中实现方法1410,或者方法1410分布在多个车辆控制器中。
方法1410包括评价对应于当前车辆驱动状况的横向动力指数(ILD)。ILD限定为:
其中,变量(a)是换算参数。车辆系统将ILD与两个阈值参数进行比较,这两个阈值参数是:对应于导数方法的稳定状态阈值(SST);对应于估计方法的动态处理阈值(DHT),其中,SST<DHT。还与两个计时器一起进行确定,这两个计时器是对车辆在每种状况下停留的持续时间进行计数的稳定状态计时器(Tss)和动态驱动计时器(Tdd)。函数LPF(x)指示普通低通滤波器函数。滑动平均滤波器是低通滤波器的示例,其可用于平滑不连续的输出过渡效应。
在操作1412中,车辆系统确定车辆是否处于稳定状态驱动状况。车辆系统接收与驱动状况相关的标志(SSD_Flag)。SSD_Flag等于1指示当前是稳定状态驱动状况,SSD_Flag等于2指示当前是动态处理状况。因此,在操作1412处,系统评价SSD_Flag,以确定是否是SSD_Flag=2。如果在操作1412处确定是肯定的,则车辆系统进行到操作1414。
在操作1414处,车辆系统评价横向动力指数(ILD),以确定ILD是否小于稳定状态阈值(SST)。SST对应于在上面参照等式28至31描述的用于确定横向速度(Vy)的导数方法。如果在操作1414处确定是否定的,则车辆系统进行到操作1416,并通过将稳定状态计时器(Tss)设置为等于0而使Tss复位。如果在操作1414处确定是肯定的,则车辆系统进行到操作1418。
在操作1418处,车辆系统通过使稳定状态计时器(Tss)加上一个执行循环时间(ts)而使稳定状态计时器(Tss)增加。在操作1420处,系统评价稳定状态计时器(Tss),以确定Tss是否等于或大于稳定时间阈值(tss_settle)。如果在操作1420处确定是肯定的,则车辆系统进行到操作1422,使动态驱动计时器(Tdd)复位。然后,在操作1424处,车辆系统通过将SSD_Flag设置为等于1而将车辆驱动状况设置为稳定状态。
如果在操作1412处确定是否定的(未在稳定状态情况下驱动车辆),则车辆系统进行到操作1426。在操作1426处,车辆系统评价横向动力指数(ILD),以确定ILD是否大于动态处理阈值(DHT)。DHT对应于在上面参照等式24至27描述的用于确定横向速度(Vy)的估计方法。如果在操作1426处确定是否定的,则车辆系统进行到操作1428,通过将动态驱动计时器(Tdd)设置为等于0而使Tdd复位。如果在操作1426处确定是肯定的,则车辆系统进行到操作1430。
在操作1430处,车辆系统通过使动态驱动计时器(Tdd)加上一个执行循环时间(ts)而使动态驱动计时器(Tdd)增加。在操作1432处,车辆系统评价动态驱动计时器(Tdd),以确定Tdd是否等于或大于激发时间阈值(tdd_excite)。如果在操作1432处确定是肯定的,则车辆系统进行到操作1434,使稳定状态计时器(Tss)复位。然后,在操作1436处,车辆系统通过将SSD_Flag设置为等于2而将车辆驱动状况设置为动态驱动。在操作1416、1424、1428或1436之后,以及在操作1420或1432处确定是否定之后,车辆系统进行到操作1438,然后返回到操作1412,以进行另一次迭代。
参照图15,示出了根据一个或多个实施例的用于基于当前车辆驱动状况来估计车辆横向速度(Vy)的方法,该方法总体上由标号1510指示。根据一个或多个实施例,使用包含在VSC内的软件代码实现该方法。在其他实施例中,在其他车辆控制器中实现方法1510,或者方法1510分布在多个车辆控制器中。方法1510通过选择用于估计横向速度(Vy)的方法,来仲裁上面参照等式24至31描述的两种方法中的一种(基于积分估计的方法或基于导数估计的方法)。
在操作1512中,车辆系统确定车辆是否处于动态驱动状况。车辆系统评价SSD_Flag以确定是否是SSD_Flag=2。如果在操作1512处确定是肯定的,则车辆系统进行到操作1514,使用在等式27中提供的积分估计来估计横向速度在操作1514之后,车辆系统进行到操作1516。在操作1516处,车辆系统使用低通滤波器(LPF)对估计的横向速度进行滤波,并更新横向速度输出(Vy)。如果在操作1512处确定是否定的,则车辆系统进行到操作1518。在操作1518处,车辆系统使用在等式30中提供的导数估计来估计横向速度在操作1518之后,车辆系统进行到操作1516,使用低通滤波器(LPF)对估计的横向速度进行滤波,并基于导数估计来更新横向速度输出(Vy)。在操作1516之后,车辆系统进行到操作1520,提供更新了的横向速度输出(Vy),然后返回到操作1512。
在一个或多个实施例中,车辆系统使用传感器横向位置偏移效应补偿策略来估计横向偏移量(lc)项。横向偏移量(lc)由下面的等式限定:
其中,是静态传感器安装偏移量。变量htφv是车身动态侧倾分量,其中,ht是从簧上质量C.G.到侧倾中心的距离。车辆侧倾角φv可从侧倾动态模型获得或者(在具有悬挂挠度传感器的实施例中)从悬挂的挠度获得。在一个或多个实施例中,从处于稳定状态的横摆角(r)获得处于稳定状态的φv,如下所示:
φv=Gr2φr(34)
将等式33和34合并,并使每侧乘以横摆角的导数提供下面的等式:
在动态驱动状况下,根据如下所示的等式36从横向加速器近似车身侧倾角(φv)。
将等式35和36合并,提供下面的等式:
其中,Gay2φ等于Kphi。
参照图16,示出了根据一个或多个实施例的用于基于当前车辆驱动状况来估计车身侧倾角(φv)的方法,该方法总体上由标号1610指示。根据一个或多个实施例,使用包含在VSC内的软件代码实现该方法。在其他实施例中,在其他车辆控制器中实现方法1610,或者方法1610分布在多个车辆控制器中。方法1610通过选择用于估计车身侧倾角(φv)的方法,来仲裁上面参照等式33至37描述的两种方法中的一种(使用横向加速度或横摆角)。
在操作1612中,车辆系统确定车辆是否处于动态驱动状况。车辆系统评价SSD_Flag,以确定是否是SSD_Flag=2。如果在操作1612处确定是肯定的,则车辆系统进行到操作1614,基于如等式36所示的横向加速度来估计车身侧倾角在操作1614之后,车辆系统进行到操作1616。在操作1616处,车辆系统使用低通滤波器(LPF)对估计的车身侧倾角进行滤波,并更新车身侧倾角输出(φv)。如果在操作1612处确定是否定的,则车辆系统进行到操作1618。在操作1618处,车辆系统基于在等式34中提供的横摆角来估计车身侧倾角在操作1618之后,车辆系统进行到操作1616,使用低通滤波器(LPF)对估计的车身侧倾角进行滤波,并更新车身侧倾角输出(φv)。在操作1616之后,车辆系统进行到操作1620,提供更新了的车身侧倾角输出(φv),然后返回到操作1612。
在一个或多个实施例中,车辆系统使用三种补偿策略中的一种来估计车身俯仰角(θ)。不同的策略对应于传感器的可用性和精度。例如,第一策略基于车辆纵向速度(Vx)和车辆纵向加速度来估计俯仰角(θ)。第二策略基于Vx、和垂向加速度来估计俯仰角(θ)。第三策略基于Vx、 以及由于由扩展卡尔曼滤波器提供的重力分量而导致的估计加速度来估计俯仰角(θ)。一般来说,使用的输入信号越多以及源信息越精确,在消除估计的耦合效应时补偿策略越好,以及总体道路坡度估计越好。
用于估计车身俯仰角(θ)的第一补偿策略被称为动态俯仰补偿。对于该策略,车辆系统基于车辆纵向速度(Vx)和车辆纵向加速度估计俯仰角(θ)。θ的低频分量可由一阶或二阶线性动态系统进行近似。造成车辆动态俯仰的推力通常根据下面的等式由作用在车身上的纵向力导致:
其中,θ与成比例,Ms是车辆簧上质量,其等于车辆质量和悬挂上的有效负荷,hcg是沿着竖直方向车辆CG在地面之上的高度。Kf和Kr分别是在前轴和后轴上的预定悬挂刚度。是等效反作用加速度,且等于ax+gsinαr,。输入到动态俯仰模型的初始推力基于如下所示:
一阶动态相对俯仰由如下所示的等式40提供:
其传递函数为:
其中,Dpitch。表示悬挂的普通阻尼效应,Kpitch表示从纵向力到俯仰推力的悬挂刚度。
图17描述了用于使用相平面增益安排规则选择Dpitch和Kpitch的图,该图总体上由标号1710指示。图l710示出了8个区域,这8个区域以笛卡尔坐标系的原点为中心。x轴对应于y轴对应于俯仰角(0)。如所示出的,俯仰角θ以θmax和θmin为边界。一般来说,在普通区域中存在线性关系,然而在更刚性的区域中存在非线性关系。在向前俯仰和向后俯仰之间,值稍微不同。
车辆系统基于Vx、和垂向加速度使用第二补偿策略估计俯仰角(θ)。根据一个或多个实施例,图l的惯性传感器52提供垂向加速度
为了方便起见,变量(ψ)可如下所示那样限定:
可通过将等式7、9和42合并并替入(αr+θ)而获得下面的等式:
可通过忽略含θ4的更高阶次的项而简化等式43,如下所示:
通过求解等式44的方根而获得用于根据第二策略计算相对俯仰角(θ)的等式,如等式45所示:
车辆系统基于Vx、以及由扩展卡尔曼滤波器提供的重力分量使用第三补偿策略估计俯仰角(θ)。下面提供近似等式:
可基于下面的等式通过回归来迭代地识别俯仰角(θ)与估计的道路坡度(αr):
其中,来自卡尔曼滤波器的输出的基于的初始补偿值。该值可等于零(0)、任何其他容许值或者来自先前循环的回归结果。
用于根据第三补偿策略估计道路坡度和相对俯仰角的等式由如下所示的等式48和49提供:
其中,和(∞)指示和随着k→∞通过回归过程收敛的值。通常在5至10步迭代内发生收敛。
这样,通过使用输入补偿的加速度偏移量来估计运动学道路坡度,以提高精度从而消除由车辆横向动力导致的偏移量和不确定性及其耦合效应,使车辆系统优于现有的方法。
车辆系统基于对应于估计的车辆质量、车辆速度、车轮扭矩(制动扭矩和动力传动系输出扭矩)以及作用在车辆上的其他拖曳力(例如,气动阻力、滚动阻力和道路坡度负荷)的输入,估计动态道路坡度(RGEdyn)。一般来说,与静态RGE和运动学RGE相比,动态RGE的品质状况更低,其原因是动态RGE在进行估计时涉及更多的信号,且每个信号中的噪声和误差将传送到最终RGE输出中。然而,动态RGE不依赖纵向加速度信号因此,动态RGE的信号覆盖范围不同于静态估计和运动学估计的信号覆盖范围。虽然动态道路坡度估计是次要的,但是对于很多车辆系统来说,车辆质量估计部分是有用的参数。
在特定车辆状况下,独立于道路坡度估计,车辆系统使用递归最小二乘(RLs)算法估计车辆质量。一般来说,当物体(例如,乘客)增加到车辆或者从车辆去除时,只有车辆质量发生改变。通常,仅在车辆停止时发生这样的改变。独立于道路坡度估计,在车辆使用车辆动态方程运动时,车辆系统估计车辆质量。相对于现有的方法,通过独立地或者并行地估计车辆质量和道路坡度,车辆系统提高每种估计的精度。车辆系统使用事件查找策略估计车辆质量,由此在特定输入的信噪比高的车辆状况下估计质量。然后,当查找事件发生时,车辆系统使用递归最小二乘(RLS)算法估计车辆质量。在车辆状况对应于质量可能改变时(例如,当车辆停止时),车辆系统重新初始化质量估计算法中的特定参数,以增加算法的灵敏度。
车辆系统使用如下提供的基于牛顿第二定律(∑F=ma)的车架动态方程来估计车辆质量:
其中,Me是在不存在有效负荷(乘客或货物)时估计的车辆质量,是车辆速度的导数,Fpwt是动力传动系的输出处的推进力,Fbrk是制动力,Fdrag是拖曳力,Fg是由于重力导致的力。
可重写等式50以求解Me,等式50还可基于车辆是否处于正常状态而分成多个等式,如下所示:
其中,Mv是估计的车辆质量,其包括在不存在有效负荷时的车辆质量(Me)和车辆有效负荷。术语“处于正常状态”通常指的是动力传动系是否连接到车轮。Jw和Jp分别是车轮和动力传动系处的转动惯量。Rw是车轮半径。rp是集总的动力传动系扭矩传递比,ηp是集总的扭矩传递效率。Jp、rp和ηp根据变速器的速比改变。对于具有多种推进装置(例如,电动机和发动机)的混合动力电动车辆来说,这些值还改变。
根据如下所示的等式52计算取决于惯性的车辆质量(Mi):
可根据如下所示的等式53基于总的动力传动系推进扭矩(Tpwt_whl)或总的再生制动扭矩(Tbrk_whl)计算作用在车轮上的等效纵向轮胎力(Fpwt):
其中,Teng和Tmot是发动机输出扭矩和电动机输出扭矩,Teng_whl和Tmot_whl是在车轮处对应的发动机输出扭矩和电动机输出扭矩。Rpe、rpm、ηpe和ηpm分别是集总的发动机扭矩传递比、集总的电动机扭矩传递比、集总的发动机的效率参数以及集总的电动机的效率参数。
可根据如下所示的等式54基于总的摩擦制动扭矩计算作用在车轮上的等效纵向制动力(Fbrk):
其中,Pbrk是在车轮组件处在车轮制动腔室中有效制动管路压力。这是在所有车轮上评价的等效项,其原因是独立的车轮制动腔室可被施加有不同的制动管路压力。Kbrk是摩擦制动效率因子,其确定在给定了特定制动压力的情况下在制动垫上可实现多大的摩擦制动扭矩。针对制动系统状况(例如,温度、湿度等),Kbrk的值改变。在一个实施例中,Kbrk对应于预定的常数。
拖曳力(Fdrag)可根据如下所示的等式55计算:
其中,ρd和Cd是恒定的气动参数,Ad是车辆横截面积,fr是滚动阻力系数。
由于重力而作用在车辆上的负荷(Fg)可根据如下所示的等式56计算:
Fg=Mvgsinαr(56)
其中,αr是道路坡度。
等式56可转换到扭矩域,如下面等式57所示:
Tg=MvgRwsinαr(57)
其中,Rw是车轮半径。
提供第一方法,用于基于下面的参数化模型使用运动学方程使用最小二乘类型的估计器来估计车辆质量同时估计道路坡度:
基于等式58的模型获得车辆质量(Mv)和道路坡度(αr)的估计结果,如下所示:
然而,同时求解两个未知数(αr和Mv)给估计带来额外的误差。
用于使用运动学方程估计车辆质量的第二方法基于下面的数学模型使车辆质量估计与道路坡度估计分开:
等式可如下简化,如在上面参照等式10描述的且在下面再现的:
其中,表示加速度偏移量,并且等于lc是横向偏移量,r是横摆角,Vy是横向速度,θ是车辆俯仰角。在一个或多个实施例中,车辆系统不接收这些值(lc、Vy和θ)作为输入。相反,车辆系统使用多种补偿策略从其他输入估计这些值。
下面的等式通过将等式61和等式10合并而获得:
变量(aδg)表示从道路坡度变量引起的加速度扰动,可根据如下所示的等式63计算:
其中,Mv是车辆质量Me和车辆有效负荷之和。
加速度扰动(aδg)造成估计误差。当车辆有效负荷小时,可忽略aδg。还可在估计的车辆质量处于稳定状态且收敛之后回归补偿aδg。加速度偏移量可被确定和用作前馈输入补偿,如在上面参照图14至图17描述的。
可重新整理等式62,以求解Me和Mv,如下所示:
补偿拖曳力(Fdrag)和道路坡度(αr),以获得车辆质量估计。然而,由于信号噪声导致难以精确地补偿这些项。与由车辆纵向动力表现出的变化相比,这些项表现出缓慢的动态变化。因此,在短时间段估计事件中,这些项可被处理为未知的常数。可获得基于差异的估计等式,如下所示:
其中,在时间指数(i)处的纵向加速度等于在车轮处的力(Fwh1)等于(Tpwt_whl-KbrkPbrk)/Rw。变量(t)是在合格的估计事件内的时间实例。Δt是预定时间间隔参数。
对于Fdrag(t)≈Fdrag(t+Δt)和增量式车辆纵向动态模型可通过在等式66和67之间相减而获得,如下所示:
其中,ΔFwhl和是在合格的估计事件内在时间间隔Δt中估计的Fwhl的值之间的差以及的值之间的差。
使用这样的方法估计车辆质量(Mv),使车辆质量估计与噪声因子(例如,气动噪声,滚动阻力和道路坡度)隔开。随着Δt→0,等式68可根据如下所示的等式69简化成导数形式:
参照图18,示出了根据一个或多个实施例的示意性框图,该框图示出了用于车辆质量估计的方法,该方法总体上由标号1810指示。根据一个或多个实施例,方法1810包括在图4的动态道路坡度估计(RGEdyn)框426内。在框1812处,车辆系统接收输入,所述输入包括加速度横摆角(r)、车辆速度(Vx)、制动压力(Pbrk)、动力传动系扭矩(Tpwt)和制动扭矩(Tbrk)。可从参照图1描述的传感器直接接收输入,或者可在CAN总线上间接接收输入。
在框1814处,车辆系统确定动态质量输入参数Yk和φk,Yk和φk分别表示系统输出和系统输入。输入参数基于计算的车辆速度(Vx)、纵向加速度动力传动系力(Fpwt)、制动力(Fbrk)和加速度偏移量基于Tpwt和Tbrk计算Fpwt和Fbrk。在框1816处,车辆系统估计动态质量(Mv)。在框1818处,车辆系统确定加速度偏移量提供加速度偏移量作为前馈RGE输入补偿。在框1820处,车辆系统估计道理坡度αr。
根据一个或多个实施例,在框1816处,车辆系统根据等式65、68或69使用递归最小二乘(RLS)估计算法来估计动态车辆质量(Mv)。RLS估计算法基于下面的参数化估计模型:
y=φ*θ(70)
递归最小二乘(RLS)估计算法如下所示那样演变:
P(k)=[1-L(k)φ(k)]P(k-1)(73)
其中,是在第k次执行循环时估计的参数θ变形。RLS估计器是滤波策略,该策略使基于测量值和预测模型输出之间的误差而更新的估计数据平均化。当待估计的参数本质上恒定时,可使用这样的估计策略。为了开始估计,将P(0)设置为等于1。
RLS等式71至73的参数根据使用哪一个车辆质量估计等式(例如,等式65、68或69)而改变。
当基于车辆质量估计等式65计算RLS等式(71至73)时,使用下面的参数:
θ=Me
当基于以微分为基础的车辆质量估计等式68计算RLS等式(71至73)时,使用下面的参数:
θ=Me
y=ΔFwhl|Δt
可选地,当基于以微分为基础的算法等式69计算RLS等式(71至73)时,使用下面的参数:
在特定合格的车辆事件下,车辆系统使用RLS来估计车辆质量。基于车辆动态变化主要在纵向上的状况而开发质量估计器。因此,车辆系统查找这样的状况,在该状况下,车辆的纵向惯性效应显著影响相关的噪声因子,或者信噪比高。在车辆制动/减速期间或者在车辆加速期间,通常发生这些合格的事件。
在一个或多个实施例中,当满足下面的条件中的每个条件时,车辆系统检测合格的车辆制动/减速事件:
1、Vx≥C_DynMassBrkEntry;
2、
3、|r|≤C_small_yaw;
4、(Fwhl<C_FBrk且且或者(Fwhl<C_FBrk2且
在一个或多个实施例中,当满足下面的条件中的每个条件时,车辆系统检测合格的车辆加速事件:
1、Vx≥C_DynMassPropEntry;
2、
3、|r|≤C_small_yaw;
4、(Fwhl>C_FProp且 且 且 且 或者(Fwhl>C_FProp2且 且
在上面的条件中,以“C”开头的变量指示校准参数。校准参数是被设计为以特定估计事件为目标的预定值。这些条件给合格的起主导作用的车辆纵向运动提供可调的构造。参数越是精确地朝着动力学意义调节,越是可实现更高的估计精度。然而,这样的精确参数将减小可用的估计时间间隔,其原因是合格的事件将更不频繁地出现。因此,合格的估计事件取决于特定程度的激进驾驶。车辆加速或制动越是有力或快速,惯性效应越是在将引起车辆动态变化时起主导作用。在估计精度和事件发生之间权衡合格的估计参数的调节,其原因是在正常驾驶过程期间,激进驾驶行为更不常见。
参照图19,示出了根据一个或多个实施例的用于估计车辆质量的方法,该方法总体上由标号1910指示。根据一个或多个实施例,使用包含在VSC内的软件代码实现该方法。在其他实施例中,在其他车辆控制器中实现方法1910,或者方法1910分布在多个车辆控制器中。方法1910可应用于加速事件和制动/减速事件,且对应于在上面列举的事件条件。
在操作1912中,车辆系统评价车辆速度(Vx),以确定Vx是否大于估计的入口速度阈值。操作1912对应于第一事件条件。对于制动事件,估计的入口速度阈值对应于校准的制动入口速度阈值(C_DynMassBrkEntry),对于加速事件,估计的入口速度阈值对应于校准的推进入口速度阈值(C_DynMassPropEntry)。如果在操作1912处确定是肯定的,则车辆系统进行到操作1914。
在操作1914处,车辆系统基于选择的车辆质量估计算法更新估计模型输入(φ)和输出(Y),如参照等式74至76描述的。在操作1916处,车辆系统计算估计条件信号(和),如参照等式75和76描述的。
在操作1918处,车辆系统评价纵向加速度和横摆角(r),以确定在横摆角足够小时车辆动态变化是否主要在纵向上。操作1918对应于第二事件条件和第三事件条件。对于制动事件,如果纵向加速度小于或等于校准的负加速度阈值(C_aRGEAxNeg),则车辆动态变化主要在纵向上。对于加速事件,如果纵向加速度等于或大于校准的正加速度阈值(C_aRGEAxPos),则车辆动态变化主要在纵向上。对于制动事件和加速事件来说,如果横摆角的幅度(|r|)小于或等于校准的小横摆角阈值(C_small_yaw),则横摆角(r)足够小。如果在操作1918处确定是肯定的,则车辆系统进行到操作1920。
在操作1920处,车辆系统评价车轮处的力(Fwhl)以及车轮处的力的更高阶导数(和),以确定是否存在合格的估计事件。操作1920对应于第四事件条件。
对于制动事件,如果:车轮处的力(Fwhl)小于校准的制动力(C_Fbrk);车轮处的力的导数小于校准的制动力的导数(C_DFbrk);车轮处的力的二阶导数小于校准的制动力的二阶导数(C_D2Fbrk),则满足第四制动条件。可选地,如果:车轮处的力(Fwhl)小于第二校准的制动力(C_Fbrk2);车轮处的力的导数小于第二校准的制动力的导数(C_DFbrk2),则也可满足第四制动条件。
对于加速事件,如果:车轮处的力(Fwhl)大于校准的推进力(C_FProp);车轮处的力的导数小于校准的上坡推进力的导数(C_DFPropUp);车轮处的力的导数大于校准的下坡推进力的导数(C_DFPropDn);车轮处的力的二阶导数小于校准的上坡推进力的二阶导数(C_D2FPropUp);车轮处的力的二阶导数小于校准的下坡推进力的二阶导数(C_D2FPropDn),则满足第四加速条件。可选地,如果:车轮处的力(Fwhl)大于第二校准的推进力(C_FProp2);车轮处的力的导数大于第二校准的推进力的导数(C_DFProp2);车轮处的力的二阶导数大于第二校准的推进力的二阶导数(C_D2FProp2),则也可满足第四加速条件。如果在操作1920处确定是肯定的,则车辆系统进行到操作1922。
在操作1922处,车辆系统执行RLS车辆质量估计算法(例如,使用等式71至73)。在操作1924处,车辆系统估计车辆质量(例如,使用等式65、68或69)。在操作1924之后,或者响应于在操作1912、1918或1920处确定是否定的,则车辆系统进行到操作1926,然后返回到操作1912。
当车辆在运动时,车辆质量不可能改变。在存在合格的估计事件时,车辆系统继续监测车辆状态并执行估计更新。假设车辆质量恒定且估计噪声是高斯零均值,则估计的车辆质量在整个估计实例上平均。然而,当车辆在缓慢运动时,以及当车辆停止持续特定的时间段时,车辆质量可改变。在这样的车辆状况下,车辆系统将RLS估计器的灵敏度调节到更灵敏的状态,从而容易检测车辆质量的改变。
参照图20,示出了根据一个或多个实施例的用于使车辆质量估计参数复位的方法,该方法总体上由标号2010指示。根据一个或多个实施例,使用包含在VSC内的软件代码实现该方法。在其他实施例中,在其他车辆控制器中实现方法2010,或者方法2010分布在多个车辆控制器中。在用于调节RLS估计器的灵敏度的特定车辆状况下,方法2010使车辆质量估计参数复位。
在操作2012处,车辆系统评价车辆速度(Vx),以确定Vx是否小于车辆停止速度阈值。如果在操作2012处确定是肯定的,则车辆系统进行到操作2014,使车辆停止/静止计时器(Tmr_stand)增加或向上计数一个循环时间,停止/静止计时器将开始记录车辆停止持续时间。在操作2016处,车辆系统通过将驾车离开计时器(Tmr_drive)设置为等于零而使Tmr_drive复位。
在操作2018处,车辆系统评价合格的估计事件的总计数(Cntr_EstTot),以确定Cntr_EstTot是否大于校准的估计的事件阈值(C_QuaEstEvnt)。如果在操作2018处确定是肯定的,则车辆系统进行到操作2020,通过将RLS估计存储器(Vmass_memo)设置为等于当前车辆质量估计(Vmass_est)而更新Vmass_memo。根据一个或多个实施例,Vmass_est对应于在图19的操作1924处估计的车辆质量(Vmass_rls)以及等式65、68和69处估计的车辆质量(Mv)。在操作2020之后,或者响应于在操作2018处确定是否定的,车辆系统进行到操作2022。
在操作2022处,车辆系统评价停止/静止计时器(Tmr_stand),以确定Tmr_stand是否大于校准的静止计时器(C_StandstilTm)。如果在操作2022处确定是肯定的,则车辆系统进行到操作2024,使RLS参数和质量估计边界复位到它们的原始值。例如,车辆系统可将RLS估计参数P(k)(来自等式73)设置为其初始值P(0)。在操作2026处,车辆系统通过将总的合格估计实例计数器(Enter_EstTot)设置为等于零而使Enter_EstTot复位。如果在操作2012处确定是否定的,则车辆系统进行到操作2028。
在操作2028处,车辆系统评价车辆速速(Vx),以确定Vx是否大于驾车离开速度阈值。如果在操作2028处确定是肯定的,则车辆系统进行到操作2030,使车辆驾车离开计时器(Tmr_drive)增加或向上计数一个循环时间,驾车离开计时器将开始记录车辆运动持续时间。在操作2032处,车辆系统评价驾车离开计时器(Tmr_drive)和停止/静止计时器(Tmr_stand),以确定Tmr_drive是否大于校准的驾车离开计时器(C_StandstilTm)以及Tmr_stand是否大于零。如果在操作2032处确定是肯定的,则车辆系统进行到操作2034。在操作2034处,车辆系统通过将停止/静止计时器(Tmr_stand)设置为等于零而使Tmr_stand复位。在操作2026或2034之后,或者响应于在操作2022、2028或2032处确定是否定的,车辆系统进行到操作2036,然后返回到操作2012。
短暂的或瞬时的合格的估计事件不足以进行车辆质量估计。一般来说,瞬时事件不满足与纵向惯性相关的第二条件。然而,即使这样的瞬时事件满足第二条件,车辆系统仍然从估计结果去除瞬时事件,其原因是在估计的车辆质量开始收敛之前车辆系统将花费一定的时间进行估计。当已经观察到足够数量的估计实例时,更新估计的车辆质量。
参照图21,示出了根据一个或多个实施例的用于更新车辆质量估计参数的方法,该方法总体上由标号2110指示。根据一个或多个实施例,使用包含在VSC内的软件代码实现该方法。在其他实施例中,在其他车辆控制器中实现方法2110,或者方法2110分布在多个车辆控制器中。当合格的估计事件超过预定的时间间隔时,方法2110更新车辆质量估计状态。
在操作2112处,车辆系统确定是否存在合格的估计事件。操作2112对应于图19的操作1912至1920。如果在操作2112处确定是肯定的,则车辆系统进行到操作2114,计算Vmass_rls,如在操作1924中描述的。在操作2116处,车辆系统使估计实例计数器(Cntr_est)增加一个循环时间。在操作2118处,车辆系统评价估计实例计数器(Cntr_est),以确定Cntr_est是否等于或大于校准的时间间隔阈值(C_EstEvnt4Update),C_EstEvnt4Update对应于估计事件的足够的估计实例。如果在操作2118处确定是肯定的,则车辆系统进行到操作2120,通过将估计更新状态(Vmass_update)设置为等于当前车辆质量估计状态(Vmass_rls)而更新Vmass_update。
在操作2122处,车辆系统评价估计实例计数器(Cntr_est),以确定Cntr_est是否等于校准的时间间隔阈值(C_EstEvnt4Update)。如果在操作2122处确定是否定的,则车辆系统进行到操作2124,使估计实例计数器(Cntr_est)增加一个循环时间。如果在操作2122处确定是肯定的,则车辆系统进行到操作2126,使累积计数器(Cntr_EstTot)增加Cntr_est,在合格的估计事件中在每个实例处使Cntr_EstTot向上计数。
在操作2128处,车辆系统评价合格的估计事件的累积计数(Cntr_EstTot),以确定Cntr_EstTot是否大于校准的合格的估计事件阈值(C_QuaEstEvnt)。如果在操作2128处确定是肯定的,则车辆系统进行到操作2130,将动态质量估计输出(Vmass_est)设置为等于RLS估计器状态(Vmass_rls)。如果在操作2128处确定是肯定的,则车辆系统进行到操作2132,将动态质量估计输出(Vmass_est)设置为等于存储的估计存储器值(Vmass_memo)。
如果在操作2112处确定是否定的,则不存在合格的估计事件,然后车辆系统进行到操作2134。在操作2134处,车辆系统使估计实例计数器(Cntr_est)减小一个循环时间。在操作2136处,车辆系统评价估计实例计数器(Cntr_est),以确定Cntr_est是否小于或等于零。如果在操作2136处确定是肯定的,则车辆系统进行到操作2138,将当前车辆质量估计状态(Vmass_rls)设置为等于估计更新状态(Vmass_update)。在操作2130、2132或2138之后,以及响应于在操作2136处确定是否定的,车辆系统进行到操作2140,然后返回到操作2112。
参照图22,根据一个或多个实施例,车辆系统使用动态边界限制车辆质量估计。单个不精确的车辆质量估计实例的存在可能劣化估计状态。因此,车辆系统建立边界,以在估计过程中避免由特定的噪声因子导致的偏离的和不合理的估计结果。如在示出的实施例中所示,车辆系统设置上边界(Ubd)和下边界(Lbd)。根据一个或多个实施例,上边界(Ubd)对应于包括车辆、乘客和货物的最大车辆总重量(GVW);下边界(Lbd)对应于小于仅包括车辆的车辆“整备”重量的值。如图22所示,在初始状况下,Ubd对应于大约2500Kg的6VW,Lbd对应于大约1700Kg的整备重量,在(时间=0s)二者之间相差大约800Kg。
由于边界朝着估计的车辆质量收敛,所以边界是动态变化的。由于在车辆运动时车辆质量通常不改变,所以边界朝着估计的车辆质量的移动平均值收敛或变窄,使得随着更多的估计信息变得可用,这两个边界之间的差减小,以滤除不合理的车辆质量估计。例如,如在示出的实施例中所示,在时间150s,Ubd减小到大约2200Kg,Lbd增加到大约1750Kg二者之间相差大约450Kg。
如在上面参照图20描述的,当满足复位条件时,估计边界将复位到原始的上限和下限。另外,在车辆累积合格的估计更新时,边界将收敛。如在上面参照图21描述的,车辆系统包括计数器(Cntr_est),每当合格的加速或减速事件发生时,使Cntr_est增加。质量估计事件标志(Mevent)在图22中示出,与其他信号按比例绘制,使得当未检测到合格的事件时Mevent等于“NQ”,当(例如,在点2212处)检测到合格的加速事件时Mevent等于“QA”,当(例如,在点2214处)检测到合格的制动事件时Mevent等于“QB”。一旦检测到合格的事件,则计数器开始每一个执行循环就向上加1计数。一旦计数器的总数超过预定值,则车辆系统将初始估计的车辆质量设置为等于RLS估计的质量,如由标号2216指示的。
参照图23,随着车辆系统累积更多的车辆质量信息,估计的质量朝着实际车辆质量收敛。图23描述了在质量已知的车辆上获得的测试结果,所述质量由实际质量波形(Mass_actual)表示。图23还包括车辆速度(Vx),合格的事件标志(Mevent)波形。当未检测到合格的事件时Mevent等于“NQ”,当检测到合格的加速事件时Mevent等于“QA”,当检测到合格的制动事件时Mevent等于“QB”。合格的制动事件对应于车辆减速,如由标号2312指示的。初始估计的车辆质量(Vmass_est)由标号2314指示。随着累积附加信息,估计的车辆质量(Vmass_est)朝着实际质量(Mass_actual)收敛,如由标号2316指示的。
图24是示出车辆质量估计的另一图形。虽然由于检测到附加信息而使车辆质量估计(Vmass_est)随着时间朝着实际质量(Mass_actual)收敛,但是当在输入信号中存在噪声时,可将车辆质量估计调节为远离实际质量,如由标号2412指示的。例如,在一个或多个实施例中,基于如在上面参照等式65描述的拖曳力(Fdrag)估计车辆质量。在特定车辆状况(例如,大风)下,拖曳力估计的精度减小。参照图22描述的边界通过限制估计的车辆质量而减小这样的噪声状况的影响。
车辆系统可基于估计的车辆质量例如通过使用等式60估计动态道路坡度。另外,估计的车辆质量可用于其他车辆系统,例如,用于具有起动/停止功能的微混合动力电动车辆的制动施加和释放检测(BARD)。
这样,通过独立于道路坡度估计使用事件查找策略估计车辆质量,车辆系统提供优点。在特定输入的信噪比高的车辆状况或者合格的事件下估计质量。然后,当发生查找事件时,车辆系统使用递归最小二乘(RLS)算法估计车辆质量。在车辆状况对应于质量可能改变时(例如,当车辆停止时),车辆系统重新初始化质量估计算法中的特定参数,以增加算法的灵敏度。车辆系统基于估计的车辆质量来估计动态道路坡度(RGEdyn)。一般来说,与静态RGE和运动学RGE相比,动态RGE的品质状况更低,其原因是动态RGE在进行估计时涉及更多的信号,且每个信号中的噪声和误差将传送到最终RGE输出中。然而,动态RGE不依赖纵向加速度信号,因此,动态RGE的信号覆盖范围不同于静态估计和运动学估计的信号覆盖范围。另外,独立于道路坡度估计,车辆系统估计车辆质量,从而将RGEdyn中的噪声与车辆质量估计隔离。
如上所述,参照图4,车辆系统确定针对静态RGE算法、运动学RGE算法和动态RGE算法中的每种的品质因子,静态RGE算法、运动学RGE算法和动态RGE算法中的每种的品质因子分别由变量QFst、QFkin、QFdym表示。车辆系统还确定总的RGE品质因子(QFRGE)。由于每种RGE算法依赖于不同的输入信号,所以每种RGE算法具有不同的范围和可用性。
针对RGE输出以及每种RGE算法(静态、运动学和动态)的估计品质因子对应于0至3之间的值,且由变量(QFRGE、QFst、QFkin和QFdyn)表示。品质因子评价估计的信息的可用性和精度。品质因子(QF)从0至3的范围内选择,其中,品质因子是3指示全品质,在全品质的情况下满足估计条件。品质因子是2指示品质劣化,此时不完全满足估计条件但是估计状态未过远地偏离其设定点。品质因子是1指示在瞬时点附近未实现合格的估计,且估计状态未实现足够的收敛,品质因子是0指示估计不可用。
车辆系统评价RGE输入信号品质和处理,如由图4中的标号412和414所指示的。车辆纵向加速度是输入信号414,且从惯性传感器52(在图1中示出)测量。如果信号的品质因子(QF)等于3且信号在特定的操作范围(例如,+/-2g)内,则该信号被认为是良好的信号。当输入信号的品质劣化时,输入信号被锁定到最后那个高品质值,当品质劣化持续超过阈值时间时,输入信号将复位到默认值。均以相同的方式处理其他输入信号,例如,制动扭矩、动力传动系扭矩、车辆速度、横摆角等。
静态RGE算法基于纵向加速度输入估计道路坡度。静态RGE适合于车辆静止状况,例如当微混合动力车辆在斜面上停止时发动机起动/停止功能。车辆系统基于加速度传感器测量信号的稳定性评价静态RGE品质。传感器稳定性由信号与阈值相比的振荡和变化水平指示,由信号的滤波的一阶导数指示。当车辆静止且车身稳定时,获得最高品质,使得加速度计的输出接近收敛的稳定状态。在当前车辆停止事件中,更高的振荡水平将降低估计品质,且当先前的运动学估计未能设置静态估计的初始值时,更高的振荡水平变得不可用。
运动学RGE算法基于对应于车辆的运动学特性(包括车辆速度、加速度和横摆角之间的关系)的输入来估计RGE。该算法使用卡尔曼滤波器处理所述输入。运动学RGE算法适合于正常车辆运动状况(例如,Vx>5kph),其原因是在车辆速度低时,一些基本计算的精度低。另外,卡尔曼滤波器是动态滤波器,通常直到估计状态收敛才能能获得良好的估计。瞬态时间取决于初始状态、真实道路坡度状态以及估计状况(例如:信号频率分量丰富度和信噪比)之间的差异。车辆系统基于评价其输入信号的品质、其估计状态、车辆状态以及估计稳定时间段来评价运动学RGE的质量。
动态RGE算法426基于对应于车辆速度、车轮扭矩(制动扭矩和动力传动系输出扭矩)以及作用在车辆上的其他拖曳力(例如,气动阻力、滚动阻力和道路坡度负荷)的输入来估计RGE。一般来说,与静态RGE和运动学RGE相比,动态RGE的品质状况更低,其原因是动态RGE在进行估计时涉及更多的信号,且每个信号中的噪声和误差将传送到最终RGE输出中。与运动学RGE算法类似,动态RGE估计的品质也以其输入信号、估计状态、车辆状态、估计品质状况以及估计收敛/稳定时间段来评价。
提供给基于卡尔曼滤波器的估计器(例如,RGEkin)的信号被分成两组:主要信号和支持信号。主要信号用于计算估计器系统等式的输入和输出。支持信号用于调节估计器模式或参数,以适应于不同的驱动状况。
当任意的主信号不是良好(例如,QF<3)时,估计结果将开始偏离真实值。车辆系统通过对劣化的主信号存在的时间段(由变量Tdgd指示)进行计数来评价劣化的信号对估计状态的影响。车辆系统还针对Tdgd对劣化影响因子Cdgi进行插值。这样的插值可基于存储在查找表中的预定数据。一般来说,确定Cdgi的原理在于:劣化的主信号存在的时间越长,对估计结果的影响越坏。估计品质指数变量根据如下所示的等式77计算:
可选地,可使用qest的预先校准的值,在存在劣化的主信号输入到估计器时,将停止估计过程。
在所有主信号均恢复良好品质之后,估计不会立即往回收敛到真实估计目标状态。在收敛阶段发生之后,可放心使用估计结果。为了表示收敛阶段,设计渐消因子fdgi,fdgi取0至1之间的值,以减小从先前的劣化信号事件累积的qest的值,使得:
在每个控制实施循环中重复该计算,以逐渐将qest的值向下减小到零。
在每个时刻,可与来自先前的事件的影响减小同时发生的是,当前发生的影响在增加。因此,最终的估计品质指数由和中的最大值确定,即,
参照图25,示出了根据一个或多个实施例的用于估计卡尔曼估计器或RLS动态质量估计器的质量的方法,该方法总体上由标号2510指示。
在操作2512中,车辆系统评价道路坡度估计(RGE)信号,以确定估计是否被锁定。当RGE不在更新时,RGE被锁定。当RGE的品质劣化时,RGE可被锁定到其最后那个高品质值。如果在操作2512处确定是肯定的,则车辆系统进行到操作2514,将评价的估计品质因子(QF_est)锁定到其当前值。变量QF_est是对应于用于运动学道路坡度估计或动态道路坡度估计的品质因子的普通值。当估计算法不满足合格的估计条件且估计被锁定时,其相关的品质评价也被锁定,直到将发生下一次估计事件为止。如果在操作2512处确定是否定的,则车辆系统进行到操作2516。
在操作2516处,车辆系统评价道路坡度估计(RGE)信号,以确定估计是否起作用。如果RGE在更新,则RGE起作用。如果在操作2516处确定是否定的,则车辆系统进行到操作2518,将品质因子设置为对应于“无数据”的值(例如,QF=0)。如果在操作2516处确定是肯定的,则车辆系统进行到操作2520。
在操作2520处,车辆系统将评价的品质指数(qest)(如在上面参照等式77和78描述的)与良好阈值参数(Q_est_ok)进行比较,以确定qest是否小于Q_est_ok,且支持信号没有劣化。如果在操作2520处确定是肯定的,则车辆系统进行到操作2522,将品质因子设置为对应于“良好”的值(例如,QF=3)。如果在操作2520处确定是否定的,则车辆系统进行到操作2524。
在操作2524处,车辆系统将评价的品质指数(qest)与故障阈值参数(Q_est_degrade)进行比较,以确定qest是否等于或大于Q_est_degrade。如果在操作2524处确定是肯定的,则车辆系统进行到操作2526,将品质因子设置为对应于“故障”的值。如果在操作2524处确定是否定的,则车辆系统进行到操作2528,将品质因子设置为对应于“劣化”的值(例如,QF=2)。在操作2514、2518、2522、2526或2528之后,车辆系统进行到操作2530,然后返回到操作2512。
参照图26A至图26C,车辆系统监测运动学道路坡度估计(RGEkin)中的变化。图26A至图26C包括在共同的时间段内获得的测量数据的3个图形,这3个图形总体上由标号2610指示。图26A示出了车辆速度(Vx),图26B示出了实际道路坡度(RGEref)和估计的道路坡度(RGEEST)。根据一个或多个实施例,RGEEST对应于对RGEst、RGEkin和RGEdyn仲裁的输出(RGEout)。图26C示出了在RGEEST和RGEref之间计算的百分比误差。
基于卡尔曼滤波器的估计(例如,用于RGEkin的估计)假设在短的时间段或距离内道路坡度状态稳定且接近于恒定状态。该假设通常是真实存在的,其原因是大多数道路的道路坡度不会快速地改变。然而,如图26B和图26C所示,在不同的道路坡度之间过渡时,可发生大的估计误差。例如,道路坡度从大约0%改变到大约7%,如由标号2612指示的。在这种过渡期间存在通常大的估计误差,如由标号2614指示的。
车辆系统包括一种方法(未示出),该方法用于监测估计变化,以基于道路坡度变化的幅度确定估计质量(精度)。在该方法中,计算估计的道路坡度的导数当导数大于阈值估计且这样的状况在比预定时间段长的时间内真实存在时,识别过渡时间段。接下来,如果同时不存在其他劣化因子,则品质因子(QF)将减小到劣化水平(QF=2)。在另一稳定时间段内,当导数信号稳定到低值时,将消除这样的劣化。一般来说,当道路坡度在短的时间段内显著改变时,该方法通知使用RGE输出的其他车辆系统:在过渡时间段内估计不完全精确。
在车辆速度极低时,道路坡度估计不可用,如表A的第6列所示。虽然在车辆速度极低(例如,小于5mph)时,没有估计可用,但是在这样的状况期间,车辆系统仍然提供精确的道路坡度估计。在车辆速度低时,由传感器(例如,车轮速度传感器或GPS)信号提供的或者从传感器信号获得的车辆速度不精确。在这样的速度,车辆系统停止估计,以避免将这样的不精确性传送到道路坡度估计(例如,RGEkin和RGEdyn)中,并将估计器状态锁定到其最后那个合格的估计结果。
车辆可以以这样的极低的车辆速度行驶延长的距离。因此,在RGE停止之后,实际道路坡度可改变。在这样的状况下,锁定的RGE输出的精度将减小。因此,车辆系统相应地降低相关的估计品质。车辆系统限定两个车辆速度阈值(Vstp和Vexit),Vstp和Vexit分别表示车辆停止速度阈值和车辆出口速度阈值。其中,Vstp大于Vexit。当车辆速度(Vx)小于Vstp时,车辆系统(使用基于卡尔曼滤波器(或者最小二乘类型的滤波器)的估计算法)通过锁定所有的估计状态而暂停或停止估计。然而,在估计停止时刻记录数据之后,在车辆速度(Vx)进一步向下减小到Vexit之前或者在Vx增加到Vstp之上之前,估计算法继续运行。当在Vexit<Vx<Vstp之间的速度范围内操作时,新的估计输出(RGEout)将继续与锁定的估计输出比较。如果当前RGE状态和锁定的RGE之间的差大于第一阈值差,则估计品质将减小到更低的水平,而指示在这样的新状况下估计输出可能不够精确。然而,如果当前RGE状态和锁定的RGE之间的差大于第二阈值差,则估计品质将减小到“无数据”水平(例如,QF=0),而指示估计结果不能使用。在车辆速度(Vx)进一步减小到Vexit之下之后,估计将完全停止,且将不再执行估计操作。一旦车辆速度(Vx)增加到阈值Vstp之上,则锁定的估计状态和记录的估计状态将重新载入到估计算法,以恢复正常估计功能。
在一个或多个实施例中,一旦车辆速度(Vx)小于Vexit,则车辆系统停止道路坡度估计器,其原因是在这样的状况下估计不精确。Vexit可被称为“蠕动”速度。然而,车辆可以以这样的蠕动速度继续行驶。
参照图27A至图27G,在车辆速度增加到更高的速度级之上之后,运动学道路坡度估计器或动态道路坡度估计器将变得可用,以及在车辆完全停止之后,静态RGE估计器可接管过来。图27A至图27G包括在共同的时间段内获得的测量数据的7个图形,这7个图形总体上由标号2710指示。图27B示出了车辆速度(Vx),图27E示出了静态品质因子(QFst),图27F示出了运动学品质因子(QFkin),图27G示出了仲裁的道路坡度估计的品质因子(QFRGE)。当车辆停止时,静态RGE变得可用。车辆速度为零由标号2712指示,对应的高品质RGEst由标号2714指示。另外,当车辆在正常运动(如由标号2716指示的)时,运动学RGE变得可用,对应的高品质RGEkin由标号2718指示。仲裁的RGE结合RGEst和RGEkin,以在大多数车辆状况下提供高品质输出,如由标号2720指示的。然而,在特定的极低车辆速度(如由标号2722指示的)时,可能不存在可用的高品质RGE,如总体上由标号2724指示的。
然而,在0和Vexit之间的低速度范围内,可能不存在可用的RGE算法。因此,在车辆以这样的低速度运动时,可能不会检测和估计道路坡度变化。为了在报告道路坡度信息时确保RGE的稳健性,车辆系统包括保守的对策,以避免在不提供关于RGE品质的信息的情况下使用道路坡度估计。
根据一个或多个实施例,车辆系统包括一种方法(未示出),该方法用于在蠕动事件中监测车辆速度(Vx)。当Vx小于Vexit(例如,小于5kph)且当前估计器品质评价良好(例如,QF=3)时,发生蠕动事件。车辆系统通过对车辆速度进行积分来计算蠕动距离(Dc)。当Dc大于第一距离阈值时,车辆系统将品质因子降低到劣化水平(QF=2)。当Dc大于第二距离阈值时,车辆系统将品质因子降低到无数据(QF=0)。当任意的RGE估计器开始提供合格的输出时,车辆系统使蠕动距离(Dc)复位到零。同时,RGE的品质将被仲裁为当前起作用的估计器的品质评价结果,以替代在蠕动事件中评价的品质。当在速度范围(Vexit<Vx<Vstp)内品质因子(QF)劣化时,车辆系统将重新开始蠕动距离计算。另外,即使车辆速度大于Vexit,当计算的估计结果与锁定的估计结果偏离超过特定阈值时,车辆系统也重新开始Dc计算。
这样,通过评价输入信号的品质、道路坡度估计的品质以及质量估计的品质并提供与每个输入和估计相关的品质因子,车辆系统优于现有的系统。车辆系统部分基于品质因子仲裁或选择估计。然后,车辆系统基于选择的估计通过对应的品质因子提供道路坡度估计和车辆质量估计。根据相关的品质因子,其他车辆系统可不同地使用道路坡度估计和车辆质量估计。
如由在上面描述的实施例验证的,通过从多个不同的估计(静态、运动学和动态)仲裁或选择道路坡度估计,车辆、车辆系统和方法优于现有技术。通过在这些不同的估计之间进行仲裁,车辆系统在宽范围的车辆状况下提供道路坡度估计。通过使用纵向加速度的动态滤波估计静态道路坡度,以及通过使用针对初始值的另一估计(其减小了在可提供合格的静态道路坡度估计值之前的延时),车辆系统还优于现有的系统。通过使用EKF估计运动学道路坡度(其提供积分式滤波策略),车辆系统还优于现有的方法。这样的积分策略不会放大信号上的噪声,因此,这样的积分策略提供改进的估计精度同时相对快速地收敛。通过使用输入补偿的加速度偏移量估计运动学道路坡度以提高精度,来消除由车辆横向动力导致的偏移量和不确定性及其耦合效应,车辆系统还优于现有的方法。通过独立于道路坡度估计使用事件查找策略估计车辆质量,车辆系统提供优点。在车辆状况对应于质量可能改变时(例如,当车辆停止时),车辆系统重新初始化质量估计算法中的特定参数,以增加算法的灵敏度。车辆系统基于估计的车辆质量来估计动态道路坡度(RGEdyn)。通过评价输入信号的品质、道路坡度估计的品质以及质量估计的品质并提供与每个输入和估计相关的品质因子,车辆系统还优于现有的系统。车辆系统部分基于品质因子仲裁或选择估计。然后,车辆系统基于选择的估计提供具有对应的品质因子的道路坡度估计和车辆质量估计。根据相关的品质因子,其他车辆系统可不同地使用道路坡度估计和车辆质量估计。
虽然已经详细描述了最佳方式,但是熟悉本领域的技术人员将认识到在权利要求的范围内的各种可选的设计和实施例。另外,各个实施的实施例的特征可被结合,以形成本发明的进一步的实施例。虽然各个实施例可能已被描述为提供优点或者在一个或多个期望的特性方面优于其他实施例或现有技术的实施方式,但是本领域的普通技术人员将认识到,一个或多个特点或特性可被折衷,以实现期望的系统属性,期望的系统属性取决于具体的应用和实施方式。这些属性可包括但不限于:成本、强度、耐久性、生命周期成本、可销售性、外观、包装、尺寸、维护保养方便性、重量、可制造性、装配容易性等。因此,在此被描述为在一个或多个特性方面不如其他实施例或现有技术的实施方式的实施例并不在本公开的范围之外,并且可被期望用于特定的应用。另外,各个实施的实施例的特征可被结合,以形成本发明的进一步的实施例。
Claims (16)
1.一种车辆,包括:
控制器,被配置成:当纵向加速度、车辆速度和横摆角中的至少一个指示发生合格的事件时,通过使用事件查找策略,基于纵向加速度、加速度偏移量和车轮扭矩使用递归最小二乘估计来产生指示车辆质量估计的输出,
其中,车辆质量估计独立于道路坡度估计,
其中,加速度偏移量包括横向速度、横向偏移量和车辆俯仰角中的至少一个,
其中,合格的事件包括加速事件和制动事件中的至少一个,
其中,加速事件还包括:指示车辆速度大于推进入口速度阈值的输入;指示纵向加速度大于正加速度阈值的输入;指示横摆角的幅度小于横摆角阈值的输入,
其中,制动事件还包括:指示车辆速度大于制动入口速度阈值的输入;指示纵向加速度小于负加速度阈值的输入;指示横摆角的幅度小于横摆角阈值的输入。
2.根据权利要求1所述的车辆,所述车辆还包括:
惯性传感器,用于提供指示纵向加速度的信号。
3.根据权利要求1所述的车辆,其中,控制器还被配置成:
将指示车辆质量估计的输出设置为默认值,直到发生预定数量的合格的事件为止。
4.根据权利要求1所述的车辆,其中,控制器还被配置成:
使用动态上边界限制车辆质量估计,所述动态上边界的初始值基于最大车辆总重量。
5.根据权利要求4所述的车辆,其中,控制器还被配置成:
使用动态下边界限制车辆质量估计,所述动态下边界的初始值基于车辆整备重量。
6.根据权利要求5所述的车辆,其中,所述动态上边界和动态下边界中的至少一个被配置成响应于合格的事件的累积而收敛。
7.根据权利要求1所述的车辆,其中,控制器还被配置成:
基于车辆质量估计、车辆速度和车轮扭矩而产生指示动态道路坡度估计的输出。
8.一种车辆系统,包括:
控制器,被配置成:
响应于合格的事件,通过使用事件查找策略,基于纵向加速度和车轮扭矩使用递归最小二乘估计而产生指示车辆质量估计的输出;
基于车辆速度、车轮扭矩和车辆质量估计而产生指示动态道路坡度估计的输出,
其中,合格的事件包括加速事件和制动事件中的至少一个,
其中,加速事件还包括:指示车辆速度大于推进入口速度阈值的输入;指示纵向加速度大于正加速度阈值的输入;指示横摆角的幅度小于横摆角阈值的输入,
其中,制动事件还包括:指示车辆速度大于制动入口速度阈值的输入;指示纵向加速度小于负加速度阈值的输入;指示横摆角的幅度小于横摆角阈值的输入。
9.根据权利要求8所述的车辆系统,其中,控制器还被配置成:
将指示车辆质量估计的输出设置为默认值,直到发生预定数量的合格的事件为止。
10.根据权利要求8所述的车辆系统,其中,控制器还被配置成:
响应于检测到合格的事件而使事件计数加1;
当事件计数大于合格的事件阈值时,产生指示车辆质量估计的输出。
11.根据权利要求8所述的车辆系统,其中,控制器还被配置成:
使用动态上边界和动态下边界限制车辆质量估计,所述动态上边界的初始值基于最大车辆总重量,所述动态下边界的初始值基于车辆整备重量。
12.根据权利要求11所述的车辆系统,其中,所述动态上边界和动态下边界中的至少一个被配置成响应于合格的事件的累积而收敛。
13.根据权利要求8所述的车辆系统,其中,控制器还被配置成:
提供指示动态道路坡度估计的可用性的动态品质因子。
14.根据权利要求8所述的车辆系统,其中,车辆质量估计包括车辆有效负荷。
15.一种用于估计动态道路坡度的方法,所述方法包括:
响应于合格的事件,通过使用事件查找策略,基于纵向加速度和车轮扭矩使用递归最小二乘估计而产生指示车辆质量估计的输出;
基于车辆速度、车轮扭矩和车辆质量估计而产生指示动态道路坡度估计的输出,
其中,合格的事件包括加速事件和制动事件中的至少一个,
其中,加速事件还包括:指示车辆速度大于推进入口速度阈值的输入;指示纵向加速度大于正加速度阈值的输入;指示横摆角的幅度小于横摆角阈值的输入,
其中,制动事件还包括:指示车辆速度大于制动入口速度阈值的输入;指示纵向加速度小于负加速度阈值的输入;指示横摆角的幅度小于横摆角阈值的输入。
16.根据权利要求15所述的方法,所述方法还包括:
将指示车辆质量估计的输出设置为默认值,直到发生预定数量的合格的事件为止。
Applications Claiming Priority (4)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
US201261695886P | 2012-08-31 | 2012-08-31 | |
US61/695,886 | 2012-08-31 | ||
US13/735,367 | 2013-01-07 | ||
US13/735,367 US8793035B2 (en) | 2012-08-31 | 2013-01-07 | Dynamic road gradient estimation |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103661395A CN103661395A (zh) | 2014-03-26 |
CN103661395B true CN103661395B (zh) | 2016-01-13 |
Family
ID=50098662
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310390790.4A Active CN103661395B (zh) | 2012-08-31 | 2013-08-30 | 动态道路坡度估计 |
Country Status (3)
Country | Link |
---|---|
US (1) | US8793035B2 (zh) |
CN (1) | CN103661395B (zh) |
DE (1) | DE102013216667A1 (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107247824A (zh) * | 2017-05-23 | 2017-10-13 | 重庆大学 | 考虑刹车和转弯影响的汽车质量‑道路坡度联合估计方法 |
Families Citing this family (94)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB2447672B (en) | 2007-03-21 | 2011-12-14 | Ford Global Tech Llc | Vehicle manoeuvring aids |
WO2012138286A1 (en) * | 2011-04-04 | 2012-10-11 | Scania Cv Ab | Estimation of road inclination |
US9708000B2 (en) | 2011-04-19 | 2017-07-18 | Ford Global Technologies, Llc | Trajectory planner for a trailer backup assist system |
US9683848B2 (en) | 2011-04-19 | 2017-06-20 | Ford Global Technologies, Llc | System for determining hitch angle |
US10196088B2 (en) | 2011-04-19 | 2019-02-05 | Ford Global Technologies, Llc | Target monitoring system and method |
US9555832B2 (en) | 2011-04-19 | 2017-01-31 | Ford Global Technologies, Llc | Display system utilizing vehicle and trailer dynamics |
US8909426B2 (en) | 2011-04-19 | 2014-12-09 | Ford Global Technologies | Trailer path curvature control for trailer backup assist |
US9499200B2 (en) | 2011-04-19 | 2016-11-22 | Ford Global Technologies, Llc | Trailer backup assist system with object detection |
US9434414B2 (en) | 2011-04-19 | 2016-09-06 | Ford Global Technologies, Llc | System and method for determining a hitch angle offset |
US9783230B2 (en) | 2011-04-19 | 2017-10-10 | Ford Global Technologies, Llc | Trailer backup assist system with off-shoot correction |
US9513103B2 (en) | 2011-04-19 | 2016-12-06 | Ford Global Technologies, Llc | Hitch angle sensor assembly |
US9493187B2 (en) | 2011-04-19 | 2016-11-15 | Ford Global Technologies, Llc | Control for trailer backup assist system |
US9937953B2 (en) | 2011-04-19 | 2018-04-10 | Ford Global Technologies, Llc | Trailer backup offset determination |
US9969428B2 (en) | 2011-04-19 | 2018-05-15 | Ford Global Technologies, Llc | Trailer backup assist system with waypoint selection |
FR3005296B1 (fr) * | 2013-05-03 | 2016-10-07 | Renault Sa | Procede d'optimisation de la consommation energetique d'un vehicule hybride |
JP6236884B2 (ja) * | 2013-06-04 | 2017-11-29 | アイシン精機株式会社 | 車両制御装置 |
ITTO20130584A1 (it) * | 2013-07-11 | 2015-01-12 | Fiat Ricerche | Stima della massa di un veicolo e della pendenza della strada |
US10023103B2 (en) * | 2013-09-13 | 2018-07-17 | J.W. Speaker, Corporation | Systems and methods for illumination control and distribution during a vehicle bank |
WO2015045502A1 (ja) * | 2013-09-30 | 2015-04-02 | 日立オートモティブシステムズ株式会社 | 車両の走行制御装置 |
US9616898B2 (en) * | 2013-12-12 | 2017-04-11 | Ford Global Technologies, Llc | System and method for determining effective road grade characteristic |
US9327709B2 (en) * | 2014-02-12 | 2016-05-03 | Ford Global Technologies, Llc | Cancelling creep torque in a hybrid vehicle |
US9340212B2 (en) * | 2014-04-01 | 2016-05-17 | GM Global Technology Operations LLC | System and method for estimating road grade based on an output of a longitudinal acceleration sensor in a vehicle |
US9623904B2 (en) | 2014-06-03 | 2017-04-18 | Ford Global Technologies, Llc | Trailer curvature control with adaptive trailer length estimation |
US9744970B2 (en) * | 2014-07-14 | 2017-08-29 | Ford Global Technologies, Llc | Estimating a trailer road grade |
US9963004B2 (en) | 2014-07-28 | 2018-05-08 | Ford Global Technologies, Llc | Trailer sway warning system and method |
US9517668B2 (en) | 2014-07-28 | 2016-12-13 | Ford Global Technologies, Llc | Hitch angle warning system and method |
US9540043B2 (en) | 2014-07-30 | 2017-01-10 | Ford Global Technologies, Llc | Trailer backup assist system with active trailer braking for curvature control |
SE538252C2 (sv) * | 2014-08-13 | 2016-04-19 | Scania Cv Ab | Förfarande och system för säkerhetsförbättring vid framförande av ett motorfordon |
US9725093B2 (en) * | 2014-09-23 | 2017-08-08 | Cummins Inc. | Vehicle controls including dynamic vehicle mass and road grade estimation during vehicle operation |
US10099575B2 (en) | 2014-09-29 | 2018-10-16 | Ford Global Technlogies, Llc | Creep compensation using road grade and mass estimation |
US9579991B2 (en) | 2014-10-02 | 2017-02-28 | Ford Global Technologies, Llc | Vehicle system and method for controlling torque delivery during transmission engagements with road grade and mass estimation |
US9533683B2 (en) | 2014-12-05 | 2017-01-03 | Ford Global Technologies, Llc | Sensor failure mitigation system and mode management |
US9607242B2 (en) | 2015-01-16 | 2017-03-28 | Ford Global Technologies, Llc | Target monitoring system with lens cleaning device |
US9522699B2 (en) | 2015-02-05 | 2016-12-20 | Ford Global Technologies, Llc | Trailer backup assist system with adaptive steering angle limits |
US10286950B2 (en) | 2015-02-10 | 2019-05-14 | Ford Global Technologies, Llc | Speed optimized trajectory control for motor vehicles |
GB2535773B (en) * | 2015-02-27 | 2019-10-30 | Jaguar Land Rover Ltd | Vehicle mass estimation |
GB2535775B (en) * | 2015-02-27 | 2019-10-30 | Jaguar Land Rover Ltd | Vehicle mass estimation |
GB2535774B (en) * | 2015-02-27 | 2019-10-30 | Jaguar Land Rover Ltd | Vehicle mass estimation |
US9616923B2 (en) * | 2015-03-03 | 2017-04-11 | Ford Global Technologies, Llc | Topographical integration for trailer backup assist system |
US9804022B2 (en) | 2015-03-24 | 2017-10-31 | Ford Global Technologies, Llc | System and method for hitch angle detection |
US9623859B2 (en) | 2015-04-03 | 2017-04-18 | Ford Global Technologies, Llc | Trailer curvature control and mode management with powertrain and brake support |
US9840240B2 (en) | 2015-04-09 | 2017-12-12 | Ford Global Technologies, Llc | Trailer backup aid speed limiting via braking |
US9744972B2 (en) | 2015-04-09 | 2017-08-29 | Ford Global Technologies, Llc | Trailer backup aid speed limiting via braking |
GB2538806B (en) | 2015-05-29 | 2021-04-07 | Sevcon Ltd | Method and apparatus |
US9676377B2 (en) | 2015-06-17 | 2017-06-13 | Ford Global Technologies, Llc | Speed limiting comfort enhancement |
US9896126B2 (en) | 2015-07-08 | 2018-02-20 | Ford Global Technologies, Llc | Jackknife detection for vehicle reversing a trailer |
DE102016116297A1 (de) | 2015-09-11 | 2017-03-16 | Ford Global Technologies, Llc | Anhängerrücksetzhilfe durch Geschwindigkeitsbegrenzung mittels Bremsen |
US10435007B2 (en) * | 2015-09-23 | 2019-10-08 | Cummins, Inc. | Systems and methods of engine stop/start control of an electrified powertrain |
US9981662B2 (en) | 2015-10-15 | 2018-05-29 | Ford Global Technologies, Llc | Speed limiting comfort enhancement |
US10611407B2 (en) | 2015-10-19 | 2020-04-07 | Ford Global Technologies, Llc | Speed control for motor vehicles |
US10384607B2 (en) | 2015-10-19 | 2019-08-20 | Ford Global Technologies, Llc | Trailer backup assist system with hitch angle offset estimation |
US9836060B2 (en) | 2015-10-28 | 2017-12-05 | Ford Global Technologies, Llc | Trailer backup assist system with target management |
US10094704B2 (en) * | 2015-10-30 | 2018-10-09 | Caterpillar Inc. | System for estimating a mass of a payload in a hauling machine |
CN106677913B (zh) * | 2015-11-10 | 2019-11-01 | 中联重科股份有限公司 | 动力与负载匹配控制方法、装置、系统及运输车辆 |
US10017115B2 (en) | 2015-11-11 | 2018-07-10 | Ford Global Technologies, Llc | Trailer monitoring system and method |
US9827818B2 (en) | 2015-12-17 | 2017-11-28 | Ford Global Technologies, Llc | Multi-stage solution for trailer hitch angle initialization |
US9610975B1 (en) | 2015-12-17 | 2017-04-04 | Ford Global Technologies, Llc | Hitch angle detection for trailer backup assist system |
US10155478B2 (en) | 2015-12-17 | 2018-12-18 | Ford Global Technologies, Llc | Centerline method for trailer hitch angle detection |
US9798953B2 (en) | 2015-12-17 | 2017-10-24 | Ford Global Technologies, Llc | Template matching solution for locating trailer hitch point |
US9934572B2 (en) | 2015-12-17 | 2018-04-03 | Ford Global Technologies, Llc | Drawbar scan solution for locating trailer hitch point |
US9796228B2 (en) | 2015-12-17 | 2017-10-24 | Ford Global Technologies, Llc | Hitch angle detection for trailer backup assist system |
US10011228B2 (en) | 2015-12-17 | 2018-07-03 | Ford Global Technologies, Llc | Hitch angle detection for trailer backup assist system using multiple imaging devices |
US9850839B2 (en) | 2016-01-29 | 2017-12-26 | Ford Global Technologies, Llc | System and method for hill ascent speed assistance |
US9919739B2 (en) | 2016-02-04 | 2018-03-20 | Ford Global Technologies, Llc | Indication of vehicle direction of travel |
US10005492B2 (en) | 2016-02-18 | 2018-06-26 | Ford Global Technologies, Llc | Trailer length and hitch angle bias estimation |
WO2017175840A1 (ja) * | 2016-04-06 | 2017-10-12 | ヤマハ発動機株式会社 | 姿勢推定装置及び輸送機器 |
US10106193B2 (en) | 2016-07-01 | 2018-10-23 | Ford Global Technologies, Llc | Enhanced yaw rate trailer angle detection initialization |
US10046800B2 (en) | 2016-08-10 | 2018-08-14 | Ford Global Technologies, Llc | Trailer wheel targetless trailer angle detection |
US10773721B2 (en) | 2016-10-21 | 2020-09-15 | Ford Global Technologies, Llc | Control method using trailer yaw rate measurements for trailer backup assist |
US10222804B2 (en) | 2016-10-21 | 2019-03-05 | Ford Global Technologies, Llc | Inertial reference for TBA speed limiting |
US10543846B2 (en) * | 2016-10-28 | 2020-01-28 | Hall Labs Llc | Determining vehicle power |
US10678237B2 (en) * | 2016-12-22 | 2020-06-09 | Panasonic Intellectual Property Corporation Of America | Information processing apparatus, operated vehicle, information processing method, and recording medium storing program |
US10604184B2 (en) | 2017-08-31 | 2020-03-31 | Ford Global Technologies, Llc | Adaptive steering control for robustness to errors in estimated or user-supplied trailer parameters |
US10710585B2 (en) | 2017-09-01 | 2020-07-14 | Ford Global Technologies, Llc | Trailer backup assist system with predictive hitch angle functionality |
US10730553B2 (en) | 2017-09-27 | 2020-08-04 | Ford Global Technologies, Llc | Adaptive steering control for robustness to errors in estimated or user-supplied trailer parameters |
CN108437998B (zh) * | 2018-01-09 | 2019-09-13 | 南京理工大学 | 基于纵向动力学的纯电动汽车坡度识别方法 |
GB2571322B (en) * | 2018-02-26 | 2020-05-27 | Jaguar Land Rover Ltd | A controller for a vehicle and a method |
US10784841B2 (en) | 2018-03-08 | 2020-09-22 | Cnh Industrial America Llc | Kalman filter for an autonomous work vehicle system |
JP6939660B2 (ja) * | 2018-03-13 | 2021-09-22 | トヨタ自動車株式会社 | 車両走行制御システム |
WO2020055677A1 (en) | 2018-09-13 | 2020-03-19 | Crown Equipment Corporation | System and method for controlling a maximum vehicle speed for an industrial vehicle based on a calculated load |
CN110161513B (zh) * | 2018-09-28 | 2023-11-14 | 腾讯科技(北京)有限公司 | 估计道路坡度的方法、装置、存储介质和计算机设备 |
CN111114551B (zh) * | 2018-10-31 | 2020-12-29 | 广州汽车集团股份有限公司 | 一种车辆坡道坡度识别方法和装置 |
US11077795B2 (en) | 2018-11-26 | 2021-08-03 | Ford Global Technologies, Llc | Trailer angle detection using end-to-end learning |
US10814912B2 (en) | 2018-11-28 | 2020-10-27 | Ford Global Technologies, Llc | Trailer backup assist system having advanced user mode with selectable hitch angle limits |
US10829046B2 (en) | 2019-03-06 | 2020-11-10 | Ford Global Technologies, Llc | Trailer angle detection using end-to-end learning |
WO2021022473A1 (zh) * | 2019-08-06 | 2021-02-11 | 深圳市大疆创新科技有限公司 | 控制方法、控制设备、可移动平台以及存储介质 |
US11420634B2 (en) * | 2020-03-25 | 2022-08-23 | Fca Us Llc | Kalman filter based road grade estimation method using accelerometer, gyroscope, and vehicle velocity |
CN111469674A (zh) * | 2020-05-06 | 2020-07-31 | 浙江吉利新能源商用车集团有限公司 | 一种用于新能源载货车辆的制动能量回收控制方法 |
CN111731308B (zh) * | 2020-06-17 | 2021-09-28 | 东风汽车有限公司 | 坡度估算方法、存储介质及电子设备 |
CN111731309B (zh) * | 2020-06-28 | 2021-08-10 | 中国第一汽车股份有限公司 | 坡度估计方法、装置、设备及车辆 |
CN112896119A (zh) * | 2021-03-11 | 2021-06-04 | 蔚来汽车科技(安徽)有限公司 | 坡道刹车压力确定方法以及确定系统、车辆 |
CN114684159A (zh) * | 2022-03-21 | 2022-07-01 | 潍柴动力股份有限公司 | 一种车辆质量估算方法、装置、电子设备和存储介质 |
US20230382405A1 (en) * | 2022-05-27 | 2023-11-30 | Curtis Instruments, Inc. | Grade severity detection |
CN116039648B (zh) * | 2023-04-03 | 2023-06-27 | 成都赛力斯科技有限公司 | 一种基于权重的坡度计算方法、装置及车辆 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102627108A (zh) * | 2012-04-11 | 2012-08-08 | 清华大学 | 一种基于高频信息提取的整车质量估算方法 |
Family Cites Families (20)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6167357A (en) | 1998-04-23 | 2000-12-26 | Cummins Engine Company, Inc. | Recursive vehicle mass estimation |
US6714851B2 (en) | 2002-01-07 | 2004-03-30 | Ford Global Technologies, Llc | Method for road grade/vehicle pitch estimation |
US6556908B1 (en) | 2002-03-04 | 2003-04-29 | Ford Global Technologies, Inc. | Attitude sensing system for an automotive vehicle relative to the road |
US6684140B2 (en) | 2002-06-19 | 2004-01-27 | Ford Global Technologies, Llc | System for sensing vehicle global and relative attitudes using suspension height sensors |
US7668645B2 (en) | 2004-10-15 | 2010-02-23 | Ford Global Technologies | System and method for dynamically determining vehicle loading and vertical loading distance for use in a vehicle dynamic control system |
US7715965B2 (en) * | 2004-10-15 | 2010-05-11 | Ford Global Technologies | System and method for qualitatively determining vehicle loading conditions |
JP4631549B2 (ja) * | 2005-06-01 | 2011-02-16 | 株式会社アドヴィックス | 車両運動安定化制御装置 |
JP4929944B2 (ja) * | 2006-09-20 | 2012-05-09 | 日産自動車株式会社 | 車両の駆動力配分制御装置 |
US7499784B2 (en) | 2007-04-09 | 2009-03-03 | General Motors Corporation | Method of selecting a transmission shift schedule |
JP2009040308A (ja) | 2007-08-10 | 2009-02-26 | Denso Corp | 路面勾配推定装置、車両用制御装置、及び車両用制御システム |
US8010252B2 (en) * | 2007-10-05 | 2011-08-30 | Ford Global Technologies | Trailer oscillation detection and compensation method for a vehicle and trailer combination |
US20090187324A1 (en) | 2008-01-23 | 2009-07-23 | Jianbo Lu | Vehicle Stability Control System and Method |
JP4582170B2 (ja) | 2008-03-27 | 2010-11-17 | トヨタ自動車株式会社 | 勾配情報演算装置、車両走行制御装置、ナビゲーションシステム |
TWI337585B (en) | 2008-10-16 | 2011-02-21 | Univ Nat Chiao Tung | Road angle estimation system and its method |
WO2010062671A1 (en) * | 2008-10-30 | 2010-06-03 | Ford Global Technologies, Llc | Vehicle and method for advising driver of same |
JP4691174B2 (ja) * | 2009-04-21 | 2011-06-01 | 本田技研工業株式会社 | 車両用勾配推定装置および車両用走行制御装置 |
GB2473436B (en) | 2009-09-09 | 2016-02-17 | Gm Global Tech Operations Inc | Method and apparatus for road surface friction estimation based on the self aligning torque |
US9403415B2 (en) | 2009-10-12 | 2016-08-02 | Ford Global Technologies | GPS based pitch sensing for an integrated stability control system |
US8579765B2 (en) | 2009-12-09 | 2013-11-12 | Ford Global Technologies, Llc | Hybrid electric vehicle hold function |
KR101230903B1 (ko) | 2010-12-08 | 2013-02-07 | 현대자동차주식회사 | 전기자동차의 등판밀림방지를 위한 크립토크 제어 방법 |
-
2013
- 2013-01-07 US US13/735,367 patent/US8793035B2/en active Active
- 2013-08-22 DE DE102013216667.2A patent/DE102013216667A1/de active Pending
- 2013-08-30 CN CN201310390790.4A patent/CN103661395B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102627108A (zh) * | 2012-04-11 | 2012-08-08 | 清华大学 | 一种基于高频信息提取的整车质量估算方法 |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107247824A (zh) * | 2017-05-23 | 2017-10-13 | 重庆大学 | 考虑刹车和转弯影响的汽车质量‑道路坡度联合估计方法 |
Also Published As
Publication number | Publication date |
---|---|
US8793035B2 (en) | 2014-07-29 |
US20140067155A1 (en) | 2014-03-06 |
DE102013216667A1 (de) | 2014-03-06 |
CN103661395A (zh) | 2014-03-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103661395B (zh) | 动态道路坡度估计 | |
CN103661393B (zh) | 运动学道路坡度估计 | |
CN103661394B (zh) | 道路坡度估计仲裁 | |
CN103661352B (zh) | 静态道路坡度估计 | |
US20220016941A1 (en) | System and method for predicting wear progression for vehicle tires | |
CN103129560B (zh) | 用于估计车辆的质量的系统和方法 | |
US9403415B2 (en) | GPS based pitch sensing for an integrated stability control system | |
EP2610605B1 (en) | Device for estimating vehicle body vibration and controller for suppressing vehicle body vibration using same | |
EP2380795B1 (en) | Vehicle condition estimating device | |
Kidambi et al. | Methods in vehicle mass and road grade estimation | |
CN105073526B (zh) | 车辆参考速度确定方法和利用这种方法的车辆控制装置 | |
GB2535773B (en) | Vehicle mass estimation | |
CN106515739B (zh) | 一种电动汽车坡道识别装置及方法 | |
CN106232446B (zh) | 用于确定惯性传感器的误差的方法 | |
EP1994388A1 (en) | Method for determining the centre of gravity for an automotive vehicle | |
WO2008050782A1 (fr) | Dispositif de commande d'amortissement des vibrations pour véhicule | |
CN109606378A (zh) | 面向非高斯噪声环境的车辆行驶状态估计方法 | |
CN104354700A (zh) | 一种基于无迹卡尔曼滤波的车辆参数在线估计方法 | |
EP2858870B1 (en) | Vehicle and method of control thereof | |
De Bruyne et al. | Online estimation of vehicle inertial parameters for improving chassis control systems | |
Lundquist et al. | Recursive identification of cornering stiffness parameters for an enhanced single track model | |
JP5337090B2 (ja) | 車両特性情報推定装置及びそれを用いた警告装置 | |
GB2535775B (en) | Vehicle mass estimation | |
GB2535774B (en) | Vehicle mass estimation | |
US20230077673A1 (en) | Systems and methods for predicting and detecting vehicle instability |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |