本说明书涉及用于设计影响车辆动力学的至少一个变量的虚拟传感器的技术,其能够使用从电子稳定性控制系统或在生产汽车和重型车辆时通常可用的其他装置测量的数据实时地给出所关心的变量的可靠估计。
这些技术特别地意在估计车辆的侧滑角和/或横向和纵向速度和/或轮胎-道路摩擦系数。在生产车辆时实时地获得这些变量的可靠估计可以允许实现了超越当前设计的显著的性能提升的对稳定性、牵引和制动控制系统的设计。
背景技术:
在过去一些年中已经提出了许多关于本发明的技术问题目的(即,影响车辆动力学的变量(例如侧滑角)的估计)的方法,所述方法的直接测量需要使用对它们在生产车辆中使用而言太复杂和/或昂贵的传感器。所有这些方法基于两步过程:首先,车辆的合适的模型(该模型包括例如具有给定的自由度的动力学或运动学公式)被推导;随后,基于所推导的模型设计合适的估计算法(其可以除了别的以外基于卡尔曼滤波器,或滑动模式观测器,或水平移动估计器,或粒子滤波器),并且所述估计算法被实现在用于使用可从esc系统或其他车辆装置获得的数据(例如,关于转向角、车轮速度、偏航率、横向加速度等等的信号)来实时地估计所关心的变量(例如侧滑角)的电子板上。
这些两步方法具有严重的缺点,其被综合地描述如下。所述估计算法在识别的模型上操作,所述模型仅仅是真实车辆动力学行为的近似描述。即使可以获得非常精确的模型,当所识别的模型是非线性时,找到最优估计(例如最小方差)在计算上是难处理的,并且随后使用在计算上易处理但是必然是近似的方法。归因于所述建模和估计步骤中的所述近似,不存在用于评估所述两步方法可以如何精确的方法。甚至对于复杂系统不容易实现估计误差的界限。甚至更相关地,所述车辆模型取决于可以根据不同的操作条件(干路或湿路、轮胎磨损状态、车辆负载等等)而改变的参数,所述不同的操作条件的实时值在正常生产车辆中没有被检测到。归因于这些问题(模型和估计近似、可变的操作条件),上面所讨论的方法看来似乎均没有已经达到设计实现可接受的估计精度的针对侧滑角、纵向和横向速度、轮胎-路面摩擦系数的估计算法的能力。这也通过如下事实被证明:使用这些方法设计的这些变量的估计器没有已经在商品化的车辆上被使得是可用的。
技术实现要素:
一个或多个实施例的目的是克服根据现有技术可实现的解决方案中的固有的限制。
根据一个或多个实施例,该目的由于具有权利要求1中所指定的特征的用于估计影响车辆动力学的变量的方法而被实现。一个或多个实施例可以涉及相应的系统、涉及虚拟传感器模块、涉及配备有这样的系统或虚拟传感器模块的车辆、以及涉及能够被加载到至少一个计算机的存储器中并且包括软件代码的一些部分的计算机程序产品,所述软件代码在所述产品在至少一个计算机上被运行时能够执行所述方法的步骤。如此处所使用的,提及这样的计算机程序产品被理解为等同于提及包含用于控制所述处理系统以便协调根据所述实施例的方法的实现的指令的计算机可读装置。提及“至少一个计算机”显然意在突出以模块化和/或分布式形式实现所提出的实施例的可能性。所述至少一个计算机可以例如在车辆的电子控制板或包括用于实现上述方法的软件代码的一些部分的所谓的电子控制单元(ecu)的级别处。
权利要求形成此处关于各种实施例提供的技术教导的组成部分。
根据此处所描述的解决方案,所述方法包括:
通过估计过程来计算所述至少一个变量的估计,所述估计过程包括:
考虑在所述车辆在各个时间间隔之上的运动期间测量的一组动力学变量,以及
在所述一组测量的动力学变量上应用关于所述变量计算的至少一个最优非线性回归函数,以估计来获得所述变量的估计,
所述最优非线性回归函数通过离线最优计算过程而被获得,所述离线最优计算过程包括:
基于获取的一组参考数据并且基于在所述车辆的运动期间测量的所述一组动力学变量,针对期望的精度等级在给定的一组操作条件中找到回归函数,该回归函数提供低于或等于所述期望的精度等级的估计误差,
通过在所述给定的一组操作条件中获取变量的一组参考数据,所述获取的一组参考数据被获得,所述变量包括对应于所测量的所述车辆的动力学变量以及所述车辆的横向和纵向速度的变量。
在各种实施例中,从测试车辆和/或从车辆模拟器获取的所述测量的变量包括转向角、横向加速度、四轮速度、偏航率、纵向加速度,并且获取一组参考数据的操作包括获取与所述测量的变量以及与横向和纵向速度相关的数据,从在测试车辆上运行的特定传感器和/或从车辆模拟器获得后者。
在各种实施例中,通过在测试车辆上的测试和/或通过所述车辆的模拟器来获取所述一组参考数据。
在各种实施例中,获得非线性回归函数的操作包括:在给定属于具有给定的衰退记忆的一类函数的函数的情况下,在这样的类中找到如下这样的函数:其利用横向速度与纵向速度的比率的反正切使在先前的时刻处测量的变量的数据向量之上计算的所述函数的差的模量的测量时间间隔之上的最大值最小化,将所找到的所述函数分配为非线性回归函数以获得所述估计。
在各种实施例中,所述侧滑角被计算为所述横向速度和所述纵向速度的估计(尤其是它们的比率的反正切)的函数、,所述纵向速度的估计被获得,作为基于与所述纵向速度相关的所述一组测量的变量的子集计算的最优非线性回归函数,所述子集包括与转向角、四轮速度和纵向加速度相关的数据,所述横向速度的估计被获得,作为基于包括横向加速度和偏航率的所述一组测量的变量的另一子集以及先前获得的所述纵向速度的估计计算的最优非线性回归函数。
在各种实施例中,用于估计影响车辆动力学的至少一个变量的传感器被实现在处理模块上,并且被配置为计算描述车辆的运动动力学的所述至少一个变量的估计,考虑在所述车辆的运动期间测量的所述一组动力学变量,并且向所述一组测量的动力学变量应用至少最优非线性回归函数以获得所述估计,所述最优非线性回归函数根据前述实施例中的任一个的方法而被获得。
在各种实施例中,所述传感器具有被包括在所述电子控制单元中或所述车辆的电子控制板中的处理模块。
在各种实施例中,所述传感器从被配置为测量包括esc(电子稳定性控制)系统的车辆的动力学变量的模块接收在所述车辆的运动期间测量的所述一组动力学变量。
在各种实施例中,所述传感器被集成在用于估计影响车辆动力学的至少一个变量的系统中,包括被配置为在其运动期间测量所述车辆的动力学变量的所述模块,并且所述处理模块被配置为基于这样的测量的动力学变量实时地计算影响车辆动力学的所述至少一个变量的估计。
附图说明
现在将仅通过参考附图的非限制性实例的方式来描述实施例,在其中:
图1显示了虚拟传感器系统的框图,所述虚拟传感器系统在车辆上在线地实现根据一个实施例的方法,以便实时地估计侧滑角;
图2表示由所述方法离线地执行以便设计所述虚拟传感器的操作的框图。
具体实施方式
随后的描述示出了旨在深入理解实施例的各种特定的细节。所述实施例可以在没有所述特定的细节中的一个或多个的情况下被实现,或者使用其他方法、部件、材料等等被实现。在其他情况下,已知的结构、材料或操作没有被详细地示出或描述,以致所述实施例的各个方面将不会是模糊不清的。
在本说明书的框架中提及“实施例”或“一个实施例”意在指示关于该实施例描述的特定配置、结构或特性被包括在至少一个实施例中。同样地,可能存在于本说明书的各个点中的诸如“在实施例中”或“在一个实施例中”的短语不一定指的是同一个实施例。此外,特定的构造、结构或特性可以在一个或多个实施例中被适当地组合。
此处所使用的参考仅仅意在方便,并且因此不限定保护的范围或所述实施例的范围。
在图1中,显示了示意性地表示用于实时地估计车辆的侧滑角的虚拟传感器的实施例的框图。
用参考标号10指示的框代表要被测量的侧滑角β与其相关联的车辆。esc系统12或所述车辆的另一系统基于所述车辆10的在其运动期间的变量的值的测量来控制所述车辆10的稳定性,其以本身已知的方式测量所述车辆10在其运动期间的下列测量的变量mq:转向角α、横向加速度ay、四轮速度ws、偏航率
所述虚拟传感器11实现用于估计影响车辆的运动的变量(尤其是在如下这样的侧滑角的情况下:所述侧滑角通过在图2中简要勾画的方法而被设计)的过程dvsβ。该方法包括如下步骤210:通过来自测试车辆10’的测量值或通过经由车辆模拟器10”的模拟,或者通过测量所述数据的一部分并且通过模拟所述数据的其他部分,获取参考数据集dd。关于通过测量参考数据进行的获取,所述测试车辆10’必须是与在其上实施所述传感器11的车辆10相同类型的车辆。
所述参考数据集dd随后被传递至所述虚拟传感器dvsβ的设计的操作220,所述虚拟传感器dvsβ随后在步骤230中被实现在软件模块11中。该模块(其被嵌入在车辆10上可用的ecu中)实时地仅接收关于由所述车辆电子稳定性控制12测量的所测量的变量mq的信号,即
在此处所描述的实施例中,所述一组参考采样的数据dd如下:
其中t表示所述获取的时刻,从δt变化至n*δt,其中δt是采样时间,并且n是获取的样本的数量。为了符号简化的目的,在下面考虑采样时间δt=1。该参考数据集dd可以包括典型地在测试车辆(尤其是汽车或重型车辆)上测量的实验数据,用于评估所考虑的车辆的动力学性能,尤其是由用于构建和测试所述车辆模型的已知的两步骤方法使用的变量的数据。然而,这样的参考数据或该参考数据的一部分也可以不是通过直接测量而是通过模拟被获取:如果所述车辆的可靠的模拟器是可用的,则可以通过这样的模拟器生成等式(1)的参考数据集dd的数据。在两种情况下,所述参考数据集dd的这样的数据必须与在针对速度、制动、道路-轮胎摩擦条件、车辆负载、驾驶风格和其他参数的给定的范围内执行的(被实施的或被模拟的)驾驶测试相关,其限定所述车辆10的不同的操作条件oc(尤其是所关心的操作条件)。
根据等式(1),所述参考数据集dd包括对应于在所述车辆运动期间在所述车辆10上测量的由所述esc系统12提供给所述虚拟传感器11的变量mq的所获取的参考数据
此处一般地描述了用于影响车辆动力学的变量的估计的解决方案。尽管下面描述的实施例主要意在获得侧滑角的估计,在这些实施例中的一些中,此处所描述的方法可以被用于估计在描述车辆的运动时所关心的其它变量,而不必也获得所述侧滑角。
此处所描述的解决方案的实施例一般基于虚拟传感器dvsβ的实现,作为离散时间非线性回归方程,其形式为:
即,给定在其处被测量的变量mq的时间t,通过在变元向量rβ(t)处评估的非线性回归函数fβ的值获得在后续的时刻t+1处的所估计的侧滑角
通过找到是下面的ε-鲁棒设计问题的解的函数fβ,设计所述虚拟传感器dvsβ的操作220被获得:利用被包含在所述参考数据集dd中以及关于它们的测量精度的信息的数据(即每个测量的值用其测量精度或误差而被考虑),针对期望的精度水平ε,找到回归函数fβ,所述回归函数fβ给出被定界为
此处呈现三个实施例实例。
在导致被指示为dvsβ的虚拟传感器的第一实施例中,所述离线操作220与找到下面的最优化问题的回归函数解
其中被包含在所述参考数据集dd中的数据被用于计算(3)的右端侧。
f是具有给定的衰退记忆的一类函数,其被定义为:
即,f类的函数f与针对每个时间t的lipschitz条件有关。所述lipschitz条件的常数是参数γ、ρ与第k次方的乘积,其中0≤γ<∞,0≤ρ<∞,并且k从1变化至m,其中m=max[ny,nw,nψ,nx,nα],即所述变量mq的不同测量间隔之中的最大间隔。
回顾前面:按照定义,所述侧滑角是β(t)=arctan(vy(t)/vx(t)和
等式(3)和(4)中的ρ和m的值是估计过程dvsβ的设计参数:越大的值被选择,则能够被实现的估计误差ε越低,但是所述估计
根据c.novara、f.ruiz、m.milanese的文章“directfiltering:anewapproachtooptimalfilterdesignfornonlinearsystems(直接滤波:针对非线性系统的最优滤波器设计的新方法)”(ieee论文集,关于自动化控制,58,89-99页,2013)中的章节d中所描述的过程的步骤5b,参数γ的值可以被选择,在所述参考数据dd上运算。
在同一文章的章节ii.d中,用于求解所述最优化问题(3)的方法被描述。
最后,让ε*被计算为:
其中,被包含在所述参考数据集dd中的数据被用于计算(5)的右端侧。
针对被包括在所述数据集dd中的所有时间t和操作条件oc的整个范围,由所述函数
使用被包含在所述参考数据集dd中的数据,通过离线操作220执行从(3)至(5)的所有计算。
在所述车辆上实时地实施此实施例的虚拟传感器(由图1中的11表示)在当前时间t处如下地计算所述侧滑角β(t)的估计β*(t):
其中,从所述esc模块12获取所需的时间处的ay,ws,ax,
因此,所述第一实施例涉及用于估计至少一个变量的方法,该变量在所述实施例中是影响所述车辆10的动力学的侧滑角β,包括:在所述车辆10的运动期间测量所述车辆10的变量mq,基于这样的测量的变量mq实时地计算由所述侧滑角β表示的所述变量的估计β*,执行通过估计过程dvsβ计算所述侧滑角β的估计的步骤230,包括:考虑一组变量mq,即从所述esc系统12获取的变量(其在所述车辆10的所述运动期间在相应的时间间隔ny、nw、nψ、nx、nα之上被测量)的向量rβ(t),并且在这样的一组测量的变量上应用关于这样的侧滑角β计算的最优非线性回归函数
导致所述车辆侧滑角的虚拟传感器的第二实施例(被指示为dvsβv)基于分别被指示为dvsvx和dvsvy的纵向速度νx和横向速度νy的两个虚拟传感器的设计。
这两个虚拟传感器由利用参考数据dd的离线操作220以如下方式被设计:
-获得如下形式的虚拟传感器dvsvx:
rx(t-1)=ws(t-1),.,ws(t-nw-1),ax(t-1),.,ax(t-nx-1),α(t-1),.,α(t-nα-1)
rx是由所述测量的变量mq的子集(即,所述四轮速度w、纵向加速度ax和转向角α)组成的向量。在各种实施例中,针对所述组成向量rx的子集的其他选择是可能的。
其中,被包含在所述参考数据集dd中的数据被用于计算(8)的右端侧。
上述最优化问题的
-获得如下形式的虚拟传感器dvsvy:
ry是由所述测量的变量mq的子集(即横向加速度ay、偏航率
其中,被包含在所述参考数据集dd中的数据被用于计算(10)的右端侧。
上述最优化问题的
-由根据此第二实施例设计的所述虚拟传感器dvsβv提供的侧滑角β的估计
即,被定义为所述横向速度vy的最优估计
-所述量值
其中,被包含在所述参考数据集dd中的数据被用于计算(11)的右端侧。
该量值提供了对根据此第二实施例设计的所述虚拟传感器dvsβv的估计误差的约束,即,针对所有时间t和被包括在所述数据集dd中的操作条件oc的整个范围,
从(7)至(11)的所有计算由离线操作220使用被包含在如上面所描述的参考数据集dd中的数据来执行。
在所述车辆上实时地实施此实施例的所述虚拟传感器dvsβv(在图1中由11表示)在当前时间t处计算执行如下操作的侧滑角β(t)的估计β*(t):
-当前时间t处的纵向速度vx(t)的估计
其中,在所需的时间处的ws,ax,α的值被从所述esc模块12在线地获取,并且
-当前时间t处的横向速度vx(t)的估计
其中,在所需的时间处的ay,
-由所述虚拟传感器dvsβv提供的当前时间t处的侧滑角β(t)的估计
其中,
该第二实施例除了给出可能比由所述第一实施例虚拟传感器dvsβ给出的更精确的侧滑角β的虚拟传感器dvsβv之外,也提供所述虚拟传感器dvsvx和dvsvy,分别给出车辆纵向和横向速度的估计。注意的是:获得这些变量的可靠估计表示相关的技术方面本身,其是针对优化牵引和制动控制系统、碰撞避免系统等等具有巨大价值的νx(t)和νy(t)的知识。
由于给出横向速度
导致被指示为dvsβμ的侧滑角的虚拟传感器的设计的方法的第三实施例设想在步骤210处操作所述参考数据集dd在l个参考子集dd1、dd2、...ddl中的划分。由所述给定的一组操作条件oc中的确定的所关心的操作条件识别的所述子集中的每一个包含针对所关心的这样的操作条件的相同值(或值范围)获取的数据。在此处所描述的实施例中,如下被考虑:根据可以被认为是影响车辆动力学的最相关的操作条件oc的道路-轮胎摩擦系数μ1,μ2,..,μl执行所述划分。为了阐述简化的目的,考虑l=2个子集中的划分,但是该方法可以被容易地扩展到l的更大值,如下面所显示的。假设在所述参考数据集dd中,利用摩擦系数
在这样的划分步骤之后,在所述离线操作220中,通过找到是ε-鲁棒设计问题的解的函数来评估下面的估计器,即:
-所述侧滑角的第一估计器以下面的形式被评估:
其中,rβ(t)在等式(2)中被给出并且
其中,被包含在所述第一参考数据集dd1中的数据被用于计算(16)的右端侧。
-所述侧滑角的第二估计器以下面的形式被评估:
其中
其中,被包含在所述第二参考数据集dd2中的数据被用于计算(17)的右端侧。
-通过所述操作220评估下面形式的偏航率的第一估计器:
其中,
其中,被包含在所述第一参考数据集dd1中的数据被用于计算(18)的右端侧。
-所述偏航率的第二估计器以下面的形式被评估:
其中,
其中,被包含在所述第二参考数据集dd2中的数据被用于计算(19)的右端侧。
从(16)至(19)的所有计算由离线操作220使用被包含在如上面所描述的参考数据集dd1和dd2中的数据来执行。所述最优化问题(16)、(17)、(18)、(19)的f*解可以通过使用先前所引用的novara-ruiz-milanese的文章中描述的算法而被获得。
在所述车辆上实时地实施此实施例的虚拟传感器(由图1中的11表示)在当前时间t处执行下面的操作。
-当前时间t处的偏航率
其中,所需的时间处的ws,ay,α的值被从所述esc模块12获取,并且
-下面的最优化问题的解λ*被计算:
即,在给定所述esc模块12处实时地测量的当前时间t处的偏航率
所述当前时间t处的侧滑角β(t)的两个估计
其中,所需的时间处的ay,ws,ax,
-最后,当前时间t处的侧滑角β(t)的估计
即,被计算为具有最优参数λ*的函数作为系数的侧滑角的所述第一和第二估计的线性组合,即所述最优参数λ*和其对一的补数的负值。
因此,刚刚描述的第三实施例通过所述偏航率的估计获得所述侧滑角β的估计
等式(20)和(21)定义了实时地给出所述偏航率的估计的两个虚拟传感器。
然而,所述偏航率不需要被估计,因为它实际上由所述esc系统12实时地测量。所述偏航率估计实际上被利用以实时地检测所述轮胎-道路摩擦系数μ的值,利用通过将所述两个偏航率估计与由所述esc系统12实际测量的值相比较而获得的作为所述最优化问题(22)的解的所述最优参数λ*。
如果所述车辆10例如正在具有实际摩擦系数μ1的道路上运行,则所述第一偏航率估计
除了上面所讨论的在虚拟传感器dvsβ和dvsβv之上的估计精度提高之外,该第三虚拟传感器dvsβμ的另外的有趣的特征是:该第三实施例允许获得被指示为dvsμ的虚拟传感器,实时地提供了当前时间t处的轮胎-道路摩擦系数μ(t)的估计μ*(t),其被计算如下:
μ*(t)=λ*(t)μ1(t)-(1-λ*(t))μ2(t)(26)
所述轮胎-道路摩擦系数的所述实时估计代表相关的技术方面本身。该信息确实对于许多车辆动力学问题具有相关价值,例如牵引和制动控制、车辆动力学控制、碰撞避免等等,但是目前其在正常生产车辆上是不可用的。用于仅估计所述轮胎-道路摩擦系数的虚拟传感器dvsμ的设计需要仅仅最优化问题(18)和(19)的离线解以及仅仅操作(20)、(21)、(22)和(26)的在线执行。
如所述的,刚刚针对l=2个子集中的划分描述的操作可以被容易地扩展到参考子集的数目l的更大值。在所述离线操作220中,从相似于等式(16)、(17)和(18)、(19)的l个等式计算l个函数
已经参考作为变量
然而,清楚的是:参考所述第三实施例描述的所述操作可以更一般地被用于基于所述系数向量λ*来检测所关心的操作条件oc的实时值:
根据所述给定的操作条件集oc中的所关心的操作条件的值,在多个(l个)参考子集dd1、dd2、...、ddl中划分所述参考数据集dd;
通过获得最优非线性回归函数的所述操作220获得影响车辆动力学的变量的对应于所述参考子集(dd1、dd2、...、ddl)的多个估计,其值在所述车辆10上被实时地测量并且取决于所述所关心的条件oc;
获得使所述测量的变量和影响所述车辆动力学的变量的这样的估计的线性组合之间的差最小化的最优化问题的系数的向量λ*的解,该变量具有由所述系数的向量确定的系数。
根据此处所描述的各种实施例的解决方案允许获得下面的优点。
根据此处所描述的各种实施例的方法和虚拟传感器允许获得所关心的变量的实时估计,其与现有技术的方法不同,允许在正常生产车辆中实现。特别地,该方法是计算易处理的,并且其允许保证在生产车辆必须运行的、其实时值在生产车辆中不被检测的不同的操作条件(例如,道路-轮胎摩擦系数、负载、轮胎状态)中可以被实现的估计精度。
当然,在不损害所述实施例的原理的情况下,结构和所述实施例的细节可以相对于纯粹作为实例已在此处描述和示出的内容而广泛地变化,而不由此脱离如随附的权利要求限定的所呈现的实施例的范围。
尽管参考由esc系统测量的变量描述了所述解决方案,该解决方案也可以使用在所述车辆上实时可用的其他测量而被实现。
在特别关注了车辆侧滑角、纵向和横向速度、轮胎-路面摩擦系数的估计的情况下,此处所描述和要求权利的解决方案被揭露,这些变量是与车辆稳定性、转向、牵引和制动的监测和控制相关的针对汽车安全系统的相关变量。然而,参考这些特定的变量不是意在被理解为以任何方式限制所述实施例,这些实施例本身也适用于其它车辆动力学变量(例如与悬架控制系统相关的垂直和侧倾变量)、与汽车领域不同的领域中相关的变量(例如航空航天和海洋交通工具中的姿态估计)。