CN111845775A - 一种分布式驱动电动汽车行驶状态与惯性参数联合估计方法 - Google Patents
一种分布式驱动电动汽车行驶状态与惯性参数联合估计方法 Download PDFInfo
- Publication number
- CN111845775A CN111845775A CN202010701509.4A CN202010701509A CN111845775A CN 111845775 A CN111845775 A CN 111845775A CN 202010701509 A CN202010701509 A CN 202010701509A CN 111845775 A CN111845775 A CN 111845775A
- Authority
- CN
- China
- Prior art keywords
- vehicle
- mass
- inertia
- calculating
- tire
- 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
Links
Images
Classifications
-
- 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
- B60W50/00—Details of control systems for road vehicle drive control not related to the control of a particular sub-unit, e.g. process diagnostic or vehicle driver interfaces
-
- 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/10—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 vehicle motion
-
- 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
- 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
- B60W2040/1315—Location of the centre of gravity
-
- 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
- B60W50/00—Details of control systems for road vehicle drive control not related to the control of a particular sub-unit, e.g. process diagnostic or vehicle driver interfaces
- B60W2050/0001—Details of the control system
- B60W2050/0019—Control system elements or transfer 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
- B60W50/00—Details of control systems for road vehicle drive control not related to the control of a particular sub-unit, e.g. process diagnostic or vehicle driver interfaces
- B60W2050/0001—Details of the control system
- B60W2050/0043—Signal treatments, identification of variables or parameters, parameter estimation or state estimation
-
- 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
-
- 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
- B60W2530/00—Input parameters relating to vehicle conditions or values, not covered by groups B60W2510/00 or B60W2520/00
- B60W2530/10—Weight
Landscapes
- Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Transportation (AREA)
- Mechanical Engineering (AREA)
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Human Computer Interaction (AREA)
- Control Of Driving Devices And Active Controlling Of Vehicle (AREA)
Abstract
本发明涉及一种分布式驱动电动汽车行驶状态与惯性参数联合估计方法,在建立分布式驱动电动汽车动力学***模型过程中,考虑车辆质量、车辆横摆转动惯量等车辆惯性参数及几何参数变化车对车辆动力学估计***的影响,采用时滞的非线性动态Dugoff轮胎模型来描述来瞬态侧偏力。利用分布式驱动电动汽车四轮驱/制动转矩、四轮转速和四轮转角能够准确获得的多信息源优势,以及纵/侧向加速度、横摆角速度、方向盘转角等低成本传感器,利用了容积积分数值计算能高效计算多维函数积分效率的特点,对高维非线性车辆动力学估计***有更强的适应性,具有低成本、强实时性、精度高、可靠性强等优点。
Description
技术领域
本发明属于分布式驱动电动汽车主动安全控制领域,具体涉及一种分布式驱动电动汽车行驶状态与惯性参数联合估计方法。
背景技术
与传统内燃机汽车或集中式驱动式电动汽车相比,分布式驱动电动汽车使用轮毂电机直接驱动四个车轮,以线控***取代传统的机械传动链接,在大幅简化底盘结构、提高传动效率的同时,可实现各车轮转矩的独立控制与快速响应,为车辆动力学主动安全***例如电子稳定程序(ESP)、侧倾稳定性控制***(RSC)等的控制提供了独特的优势,将极大地改善车辆操纵舒适性和紧急工况下车辆平顺性。然而,准确而实时地获得电动汽车行驶过程中的状态参数信息是实现车辆底盘动力学主动安全控制的前提。由于测量车辆状态参数信息的传感器价格比较昂贵,且信号的可靠性等技术问题尚未完全解决,车辆运行的关键状态(如,车辆质心侧偏角等)只能被观测或估计。从成本节约和工程实际应用来看,如何利用现有的车载传感器测量信息,准确获得不易直接测量的状态参数,是电动汽车成本控制所面临的急迫问题。
目前国内外学者对车辆动力学***状态估计研究多侧重于车辆速度、车辆质心侧偏角等车辆状态估计,但对车辆惯性参数估计关注相对较少,实际上,分布式驱动电动汽车舍弃了传统内燃机、传统的变速箱等复杂的传动机构以及轮毂电机与电池等的安装,导致分布式驱动电动汽车的簧载质量发生了较大的变化,其惯性参数(车辆质量、车辆横摆转动惯量)的变化对电动汽车左右车轮上垂直载荷的重新分配影响着轮胎与路面之间相互作用的轮胎侧向力等侧偏特性及侧倾状态,进而决定着电动汽车的操纵稳定性、平顺性等。因此,对电动汽车行驶过程中车辆惯性参数及车辆侧倾、瞬态侧偏力等状态的准确联合估计也成为电动汽车动力学安全控制亟待解决的问题。
另外,车辆动力学***本质上是一个高维的非线性***,尤其当车辆运行于极限操纵工况下,车辆动力学呈现出高维强非线性特征,虽然扩展式卡尔曼滤波(EKF)在车辆工程估计领域使用较多,但EKF运用Tayfor级数对车辆非线性***进行一阶泰勒近似带来的线性化高阶截断估计误差明显,如何适应车辆高维非线性动力学的状态参数实时估计也是应该考虑的重要问题。
发明内容
本发明的目的是针对已有技术的不足,提供一种分布式驱动电动汽车行驶状态与惯性参数联合估计方法,基于建立的受惯性参数影响的车辆非线性侧倾动力学估计模型,利用分布式驱动电动汽车提供的车载信息测量量,设计并联双容积卡尔曼滤波算法,对电动汽车行驶过程中车辆侧倾、瞬态侧偏力等状态以及整车质量、横摆转动惯量等车辆惯性参数进行联合估计,具有强实时性、精度高、可靠性强等优点。
本发明采用以下技术方案:
一种分布式驱动电动汽车行驶状态与惯性参数联合估计方法,采用双容积卡尔曼滤波算法,具体包括如下步骤:
(1)构建考虑车辆惯性参数变化的侧倾动力学估计模型:
构建考虑车辆质量、车辆横摆转动惯量等惯性参数变化的八自由度车辆非线性侧倾动力学模型,建立能反映车辆质心侧倾角、轮胎侧向力以及整车质量、横摆转动惯量变化的包括车辆纵向、侧向、横摆运动、侧倾运动及四轮转动在内整车动力学估计模型;
(2)建立时滞的非线性动态Dugoff轮胎模型:
利用分布式驱动电动汽车驱动力矩可直接测量,可计算车辆在行驶过程中的纵向力;而轮胎侧偏力则采用时滞的非线性动态Dugoff轮胎模型来描述,其具有动态反应瞬态侧偏力的特点;
(3)设计双容积卡尔曼滤波器的状态和惯性参数估计***:
根据构建的受惯性参数影响的八自由度车辆动力学模型和轮胎模型,设计双容积卡尔曼滤波器的状态和惯性参数估计***;
(4)确定双容积卡尔曼滤波观测器具体步骤:
进一步,所述步骤(1)的所建立的考虑惯性参数变化的包括车辆纵向、侧向、横摆运动、侧倾运动及四轮转动在内的分布式驱动电动汽车整车动力学方程如下:
上述方程中:
上述式中,Vx、Vy分别为车辆质心的纵向和侧向速度;ωz为车辆质心的横摆角速度;β、分别为车辆质心侧偏角和侧倾角;Ms、M分别表示车辆簧载质量和总质量;Fxij、Fyij分别是车辆第i、j轮胎的纵向、侧向力,其中i=f、r,j=l、r,且fl、fr、rl、rr分别表示左前、右前、左后、右后轮胎;Fw、Ff分别是车辆空气阻力与地面轮胎滚动阻力;C为空气阻力系数;ρ为空气密度;A为汽车正面迎风面积;ax、ay分别为车辆纵向与侧向加速度;μ为已知路面附着系数;δfl、δfr分别为前轮左右轮的转向角;Lf、Lr分别为质心到车辆前后轴的水平距离;Bl、Br分别为质心到左右车轮中心的水平距离;hs是簧载质量质心到侧倾轴垂直距离,Kφf,Kφr是前后侧倾刚度,Kφf,Kφr是前后侧倾阻尼系数。
当电动汽车货物或乘客被加载,车辆的质心发生改变,假加载的质量质心位置相对于原始坐标系坐标矢量为:
当载荷质量Mp被加载,车辆的侧倾转动惯量改变为
式中,Ixxo为车辆未加载时的侧倾转动惯量;
车辆的横摆与侧倾运动的协转动惯量改变为
Ixz=xp(zp+hs)Mp
车辆的横摆转动惯量改变为
式中,Izzo为车辆未加载时的横摆转动惯量。
同时,在考虑载荷参数变化的情况下,车辆动力学模型的质心也相应发生改变,位于原来坐标系下的新的质心位置坐标为:
式中xp、yp、zp分别为载荷在原车辆坐标系下的纵向、横向、垂向坐标值;xn、yn、zn为车辆加载时的质心纵向、横向、垂向坐标。
当车辆载荷参数发生变化后,其新的车辆几何坐标位置参数计算如下式:
式中,hc为侧倾轴到地面的高度;L、B为车辆前后轴的水平距离和车辆左右车轮的水平距离;Lf、Lr分别为未加载时车辆前后轴到质心的水平距离;Bf、Br分别为车辆未加载时左右轮到质心的水平距离。
分布式驱动电动汽车每个车轮上的电机给各车轮提供驱动力矩,当车辆行驶时,车轮旋转动力学方程为:
其中,ωij为各车轮旋转角速度,Iωij为各车轮转动惯量,Tdij为各车轮驱动力矩,Rij为车轮滚动半径。
进一步,所述步骤(2)的轮胎力模型建立,利用分布式驱动电动汽车驱动力矩可直接测量,根据车轮旋转动力学可计算车辆在行驶过程中的纵向力;而轮胎侧向力则采用非线性瞬态Dugoff轮胎模型来计算。
当车辆载荷改变时,电动汽车轮胎与路面之间相互作用的轮胎侧向力瞬态侧偏特征可采用非线性瞬态Dugoff轮胎描述为:
非线性Dugoff轮胎的准静态侧向力为
其中,
上式中,Cxij,Cyij分别为轮胎纵向、侧向刚度,λij、αij分别为轮胎滑移率侧偏角;
当车辆载荷发生变化,车辆的车轮垂直载荷Fzij变化表示为:
上述式,musf、musr为车辆前后非簧载质量;huf、hur为前、后轴侧倾中心离地高度;
车轮轮胎侧偏角αij计算为:
进一步,所述步骤(3)双容积卡尔曼滤波器的状态估计***的状态方程和观测方程,离散化后可表示为如下形式:
上式中
其中,
上式是中Ts采样时间。
上述式中,xp(k)=[M,Izz]T分别为车辆非线性动力学观测器***的状态矢量和参数矢量,u(k)=[δf,wij,Tij]T和y(k)=[ωz,ax]T分别为车辆非线性动力学观测器***的输入矢量和量测矢量,qs、rs分别为***的过程噪音和量测噪音,两者为***互不相关且均值为零的高斯白噪声。
进一步,惯性参数估计***可以进一步被构造:
在上述参数估计***中,qp、rp分别为***的过程噪音和量测噪音,d(k)=[ωz,ax]T为量测矢量。
其中:
同时
所述步骤(4)的双容积卡尔曼滤波观测器具体包括以下步骤:
2)参数变量的时间更新:
2-1)误差协方差矩阵Pp(t-1)进行奇异值分解,并得到xpj(t-1);
2-4)计算协方差预测值Pp(t|t-1);
3)状态变量的时间更新:
3-1)误差协方差矩阵Pp(t-1)进行奇异值分解,并得到xsj(t-1);
3-4)计算协方差预测值Ps(t|t-1);
4)参数变量的量测更新:
4-1)误差协方差矩阵Pp(t|t-1)进行奇异值分解,并得到xpj(t|t-1);
4-2)计算容积点dpj(t|t-1);
4-4)计算新息协方差预测值Ppdd(t|t-1);
4-5)计算互协方差Ppxd(t|t-1);
4-6)计算增益矩阵Kpt;
4-8)计算误差协方差矩阵Pp(t);
5)状态变量的量测更新:
5-1)误差协方差矩阵Ps(t|t-1)进行奇异值分解,并得到xsj(t|t-1);
5-2)计算容积点ysj(t|t-1);
5-4)计算新息协方差预测值Psyy(t|t-1);
5-5)计算互协方差Psxy(t|t-1);
5-6)计算增益矩阵Kst;
5-8)计算误差协方差矩阵Ps(t)。
所述步骤(4)的Carsim与Matlab/Simulink联合仿真是由Carsim提供实车运行过程中的仿真模型数据,传送给Simulink,其中在Simulink中搭建分布式驱动***,以给电动汽车提供动力源。
与现有技术相比,本发明具有如下显而易见的突出性实质性特点和显著的进步:
1.本发明在建立分布式驱动电动汽车动力学联合估计***模型过程中,考虑车辆质量、车辆横摆转动惯量等车辆惯性参数及几何参数变化车对车辆动力学模型及轮胎模型的影响,发展了受惯性参数与几何参数变化影响的车辆非线性侧倾动力学估计模型;
2.在对分布式驱动电动汽车行使状态与惯性参数联合估计过程中,本发明引入贝叶斯非线性滤波框架下三阶球面-相径容积准则求解贝叶斯滤波中的“非线性函数×高斯密度”数值积分的容积卡尔曼滤波算法,采用三阶球面-相径容积准则对高斯加权积分进行近似,充分利用了容积积分数值计算能高效计算多维函数积分效率的特点,与其他滤波相比,对高维非线性车辆动力学估计***有更强的适应性;
3.本发明利用分布式驱动电动汽车提供的车载信息测量量,设计并联双容积卡尔曼滤波算法,对电动汽车行驶过程中车辆侧倾、瞬态侧偏力等状态以及整车质量、横摆转动惯量等车辆惯性参数进行联合估计,在车辆行驶状态估计的过程中考虑自身到惯性参数的变化使两者相互校正,能够用于分布式驱动电动汽车主动安全控制对车辆行驶过程中状态与惯性参数的测量与估计需要,具有低成本、强实时性、精度高、可靠性强等优点。
附图说明
图1为本发明能分布式驱动电动汽车行驶状态与惯性参数联合估计方法流程图。
图2为本发明分布式驱动电动汽车行驶状态与惯性参数联合估计方法的原理图。
图3为本发明的受惯性参数影响的非线性车辆动力学侧倾模型。
图4为本发明的车轮动力学模型。
图5为本发明的双容积卡尔曼滤波算法流程图。
具体实施方式
下面将结合附图对本发明的优选实施例作详细说明。
实施例一:
在本实施例中,参加图1,一种分布式驱动电动汽车行驶状态与惯性参数联合估计方法,所述方法包括以下过程:
步骤S1、构建考虑车辆惯性参数变化的侧倾动力学估计模型:
构建考虑车辆质量、车辆横摆转动惯量的惯性参数变化的八自由度车辆非线性侧倾动力学模型,建立能反映车辆质心侧倾角、轮胎侧向力以及整车质量、横摆转动惯量变化的包括车辆纵向、侧向、横摆运动、侧倾运动及四轮转动在内整车动力学估计模型;
步骤S2、建立时滞的非线性动态Dugoff轮胎模型:
利用分布式驱动电动汽车驱动力矩进行直接测量,计算车辆在行驶过程中的纵向力;而轮胎侧偏力则采用时滞的非线性动态Dugoff轮胎模型来描述,其具有动态反应瞬态侧偏力的特点;
步骤S3、设计双容积卡尔曼滤波器的状态和惯性参数估计***:
根据构建的受惯性参数影响的八自由度车辆动力学模型和轮胎模型,设计双容积卡尔曼滤波器的状态和惯性参数估计***;
步骤S4、确定双容积卡尔曼滤波观测器具体步骤:
实施例二:
本实施例与实施例一基本相同,特别之处在于:
在本实施例中,参加图1和图2,一种分布式驱动电动汽车行驶状态与惯性参数联合估计方法,步骤如下:
1.建立八自由度受惯性参数影响的非线性车辆动力学侧倾模型:
分布式驱动电动汽车动力学模型及车轮动力学模型如图3、图4所示,定义车辆坐标系的原点位于整车质心(CG)处,建立的考虑惯性参数变化的包括车辆纵向、侧向、横摆运动、侧倾运动及四轮转动在内的八自由度整车动力学方程如下:
纵向运动:
侧向运动:
横摆运动:
侧倾运动:
上述方程中:
Fw=CρAVx 2/2,Ff=μMg
上述式中,Vx、Vy分别为车辆质心的纵向和侧向速度;ωz为车辆质心的横摆角速度;β、分别为车辆质心侧偏角和侧倾角;Ms、M分别表示车辆簧载质量和总质量;Fxij、Fyij分别是车辆第i、j轮胎的纵向、侧向力,其中i=f、r,j=l、r,且fl、fr、rl、rr分别表示左前、右前、左后、右后轮胎;Izz、Ixx为分别为车辆的横摆转动惯量和侧倾转动惯量;Ixz为车辆的横摆与侧倾运动的协转动惯量;Fw、Ff分别是车辆空气阻力与地面轮胎滚动阻力;C为空气阻力系数;ρ为空气密度;A为汽车正面迎风面积;ax、ay分别为车辆纵向与侧向加速度;μ为已知路面附着系数;δfl、δfr分别为前轮左右轮的转向角;Lf、Lr分别为质心到车辆前后轴的水平距离;Bl、Br分别为质心到左右车轮中心的水平距离;hs是簧载质量质心到侧倾轴垂直距离,Kφf,Kφr是前后侧倾刚度,Kφf,Kφr是前后侧倾阻尼系数。
当电动汽车货物或乘客被加载,车辆的质心发生改变,假加载的质量质心位置相对于原始坐标系坐标矢量为:
当载荷质量Mp被加载,利用并行轴原理,车辆的侧倾转动惯量改变为:
式中,Ixxo为车辆未加载时的侧倾转动惯量;
同理,车辆的横摆与侧倾运动的协转动惯量改变为:
Ixz=xp(zp+hs)Mp
同理,车辆的横摆转动惯量改变为:
式中,Izzo为车辆未加载时的横摆转动惯量。同时,在考虑载荷参数变化的情况下,车辆动力学模型的质心也相应发生改变,利用杠杆原理,位于原来坐标系下的新的质心位置坐标为:
式中xp、yp、zp分别为载荷在原车辆坐标系下的纵向、横向、垂向坐标值;xn、yn、zn为车辆加载时的质心纵向、横向、垂向坐标。
当得车辆载荷参数发生变化后,其新的车辆几何坐标位置参数计算如下式:
式中,hc为侧倾轴到地面的高度;L、B为车辆前后轴的水平距离和车辆左右车轮的水平距离;Lf、Lr分别为未加载时车辆前后轴到质心的水平距离;Bf、Br分别为车辆未加载时左右轮到质心的水平距离。
分布式驱动电动汽车每个车轮上的电机给各车轮提供驱动力矩,当车辆行驶时,根据力矩平衡原理可得车轮旋转动力学方程为:
其中,ωij为各车轮旋转角速度,Iωij为各车轮转动惯量,Tdij为各车轮驱动力矩,Rij为车轮滚动半径。
2.构建时滞的非线性动态Dugoff轮胎模型:
利用分布式驱动电动汽车驱动力矩可直接测量,根据车轮旋转动力学可计算车辆在行驶过程中的纵向力;而轮胎侧向力则采用Dugoff轮胎模型来计算。
当车辆载荷改变时,为刻画电动汽车左右车轮上垂直载荷的重新分配影响轮胎与路面之间相互作用的轮胎侧向力瞬态侧偏特征,利用半物理半经验方法引入时滞的“松弛长度(Relaxation Length)”对Dugoff稳态轮胎模型进行修正来反映轮胎的动态特性,非线性瞬态Dugoff轮胎侧向力描述为:
非线性Dugoff轮胎的准静态侧向力如下
其中,
上式中,Cxij,Cyij分别为轮胎纵向、侧向刚度,λij、αij分别为轮胎滑移率侧偏角;
当车辆载荷发生变化,汽车非簧载质量侧倾也会产生载荷转移,因此车辆的车轮垂直载荷变化表示为:
上述式,musf、musr为车辆前后非簧载质量;huf、hur为前、后轴侧倾中心离地高度;
车轮轮胎侧偏角αij计算如下:
3.设计双容积卡尔曼滤波器的状态和惯性参数估计***:
1)对建立的受惯性参数影响到车辆动力学估计模型,得离散化后的状态估计***状态方程和观测方程:
具体为:
其中,
上式是中Ts采样时间。
离散化后的参数估计***为:
具体为:
4.设计双容积卡尔曼滤波观测器,算法过程如图5所示,具体步骤如下:
1)初始化:
2)参数变量的时间更新
2-1)误差协方差矩阵Pp(t-1)进行奇异值分解
2-4)计算协方差预测值Pp(t|t-1)
其中,Qp为过程噪声协方差矩阵。
3)状态变量的时间更新
3-1)误差协方差矩阵Ps(t-1)进行奇异值分解
3-4)计算协方差预测值Ps(t|t-1)
其中,Qs为过程噪声协方差矩阵。
4)计算参数变量的量测更新
4-1)误差协方差矩阵Pp(t|t-1)进行奇异值分解
4-2)计算容积点dpj(t|t-1)
4-4)计算新息协方差预测值Ppdd(t|t-1)
其中,Rp为测量噪声协方差矩阵。
4-5)计算互协方差Ppxd(t|t-1)
4-6)计算增益矩阵Kpt
4-8)误差协方差矩阵Pp(t)
5)状态变量的量测更新
5-1)误差协方差矩阵Ps(t|t-1)进行奇异值分解
5-2)计算容积点ysj(t|t-1)
5-4)计算新息协方差预测值Psyy(t|t-1)
其中,Rs为测量噪声协方差矩阵。
5-5)计算互协方差Psxy(t|t-1)
5-6)计算增益矩阵Kst
5-8)误差协方差矩阵Ps(t)
5.为验证双容积卡尔曼滤波的状态估计算法的有效性,采用Carsim与Matlab/Simulink联合仿真的平台进行试验。由Carsim提供实车运行过程中的仿真模型数据,传送给Simulink,其中在Simulink中搭建分布式驱动***,以给电动汽车提供动力源。
综合上述实施例可知,本发明涉及一种分布式驱动电动汽车行驶状态与惯性参数联合估计方法,在建立分布式驱动电动汽车动力学***模型过程中,考虑车辆质量、车辆横摆转动惯量等车辆惯性参数及几何参数变化车对车辆动力学估计***的影响,发展了受惯性参数等变化影响的包括车辆纵向、侧向、横摆运动、侧倾运动及四轮转动在内的车辆非线性侧倾动力学估计模型,采用时滞的非线性动态Dugoff轮胎模型来描述来瞬态侧偏力。利用分布式驱动电动汽车四轮驱/制动转矩、四轮转速和四轮转角能够准确获得的多信息源优势,以及纵/侧向加速度、横摆角速度、方向盘转角等低成本传感器,应用动力学理论和信息融合技术,基于贝叶斯非线性滤波框架下三阶球面-相径容积准则求解贝叶斯滤波中的“非线性函数×高斯密度”数值积分的双容积卡尔曼滤波理论设计车辆行驶状态与惯性参数联合估计算法,充分利用了容积积分数值计算能高效计算多维函数积分效率的特点,对高维非线性车辆动力学估计***有更强的适应性,具有低成本、强实时性、精度高、可靠性强等优点。
上面对本发明实施例结合附图进行了说明,但本发明不限于上述实施例,还可以根据本发明的发明创造的目的做出多种变化,凡依据本发明技术方案的精神实质和原理下做的改变、修饰、替代、组合或简化,均应为等效的置换方式,只要符合本发明的发明目的,只要不背离本发明的技术原理和发明构思,都属于本发明的保护范围。
Claims (6)
1.一种分布式驱动电动汽车行驶状态与惯性参数联合估计方法,其特征在于,所述方法包括以下过程:
步骤S1、构建考虑车辆惯性参数变化的侧倾动力学估计模型:
构建考虑车辆质量、车辆横摆转动惯量的惯性参数变化的八自由度车辆非线性侧倾动力学模型,建立能反映车辆质心侧倾角、轮胎侧向力以及整车质量、横摆转动惯量变化的包括车辆纵向、侧向、横摆运动、侧倾运动及四轮转动在内整车动力学估计模型;
步骤S2、建立时滞的非线性动态Dugoff轮胎模型:
利用分布式驱动电动汽车驱动力矩进行直接测量,计算车辆在行驶过程中的纵向力;而轮胎侧偏力则采用时滞的非线性动态Dugoff轮胎模型来描述,其具有动态反应瞬态侧偏力的特点;
步骤S3、设计双容积卡尔曼滤波器的状态和惯性参数估计***:
根据构建的受惯性参数影响的八自由度车辆动力学模型和轮胎模型,设计双容积卡尔曼滤波器的状态和惯性参数估计***;
步骤S4、确定双容积卡尔曼滤波观测器具体步骤:
2.根据权利要求1所述分布式驱动电动汽车行驶状态与惯性参数联合估计方法,其特征在于,在所述步骤S1中,非线性车辆动力学模型的方程为:
上述方程中:
∑Mzi=(Fyflsinδfl-Fxflcosδfl)Bln+(Fxflsinδfl+Fyflcosδfl)Lfn+(Fxfrcosδfr-Fyfrsinδfr)Brn+(Fxfrsinδfr+Fyfrcosδfr)Lfn+(FxrrBrn-FxrlBln)-(Fyrr+Fyrl)Lrn
上述式中,Vx、Vy分别为车辆质心的纵向和侧向速度;ωz为车辆质心的横摆角速度;β、分别为车辆质心侧偏角和侧倾角;Ms、M分别表示车辆簧载质量和总质量;Fxij、Fyij分别是车辆第i、j轮胎的纵向、侧向力,其中i=f、r,j=l、r,且fl、fr、rl、rr分别表示左前、右前、左后、右后轮胎;Fw、Ff分别是车辆空气阻力与地面轮胎滚动阻力;C为空气阻力系数;ρ为空气密度;A为汽车正面迎风面积;ax、ay分别为车辆纵向与侧向加速度;μ为已知路面附着系数;δfl、δfr分别为前轮左右轮的转向角;Lf、Lr分别为质心到车辆前后轴的水平距离;Bl、Br分别为质心到左右车轮中心的水平距离;hs是簧载质量质心到侧倾轴垂直距离,Kφf,Kφr是前后侧倾刚度,Kφf,Kφr是前后侧倾阻尼系数。
3.根据权利要求1所述分布式驱动电动汽车行驶状态与惯性参数联合估计方法,其特征在于,在所述步骤S1中,在非线性车辆动力学模型的方程中,当电动汽车货物或乘客被加载,车辆的质心发生改变,假加载的质量质心位置相对于原始坐标系坐标矢量为:
当载荷质量Mp被加载,车辆的侧倾转动惯量改变为
式中,Ixxo为车辆未加载时的侧倾转动惯量;
车辆的横摆与侧倾运动的协转动惯量改变为
Ixz=xp(zp+hs)Mp
车辆的横摆转动惯量改变为
式中,Izzo为车辆未加载时的横摆转动惯量;
同时,在考虑载荷参数变化的情况下,车辆动力学模型的质心也相应发生改变,位于原来坐标系下的新的质心位置坐标为:
式中xp、yp、zp分别为载荷在原车辆坐标系下的纵向、横向、垂向坐标值;xn、yn、zn为车辆加载时的质心纵向、横向、垂向坐标;
当车辆载荷参数发生变化后,其新的车辆几何坐标位置参数计算如下式:
式中,hc为侧倾轴到地面的高度;L、B为车辆前后轴的水平距离和车辆左右车轮的水平距离;Lf、Lr分别为未加载时车辆前后轴到质心的水平距离;Bf、Br分别为车辆未加载时左右轮到质心的水平距离;
分布式驱动电动汽车每个车轮上的电机给各车轮提供驱动力矩,当车辆行驶时,车轮旋转动力学方程为:
其中,ωij为各车轮旋转角速度,Iωij为各车轮转动惯量,Tdij为各车轮驱动力矩,Rij为车轮滚动半径。
4.根据权利要求1所述分布式驱动电动汽车行驶状态与惯性参数联合估计方法,其特征在于,在所述步骤S2中,进行轮胎力模型建立时,利用分布式驱动电动汽车驱动力矩进行直接测量,根据车轮旋转动力学计算车辆在行驶过程中的纵向力;而轮胎侧向力则采用非线性瞬态Dugoff轮胎模型来计算;
当车辆载荷改变时,对于电动汽车轮胎与路面之间相互作用的轮胎侧向力瞬态侧偏特征,采用非线性瞬态Dugoff轮胎描述为:
非线性Dugoff轮胎的准静态侧向力为
其中,
上式中,Cxij,Cyij分别为轮胎纵向、侧向刚度,λij、αij分别为轮胎滑移率侧偏角;
当车辆载荷发生变化,车辆的车轮垂直载荷Fzij变化表示为:
上述式,musf、musr为车辆前后非簧载质量;huf、hur为前、后轴侧倾中心离地高度;
车轮轮胎侧偏角αij计算为:
5.根据权利要求1所述分布式驱动电动汽车行驶状态与惯性参数联合估计方法,其特征在于,在所述步骤S3中,双容积卡尔曼滤波器的状态估计***的状态方程和观测方程,离散化后表示为如下:
上式中
其中,
∑Fyij(t-1)=Fyfl(t-1)cos(δfl(t-1))+Fyfr(t-1)cos(δfr(t-1))+(Fyrl(t-1)+Fyrr(t-1))+Fxfl(t-1)sin(δfl(t-1))+Fxfr(t-1)sin(δfr(t-1))
∑Mzi(t-1)=[Fyfl(t-1)sinδfl-Fxfl(t-1)cosδfl]Bln+[Fxfl(t-1)sinδfl+Fyfl(t-1)cosδfl]Lfn+[Fxfr(t-1)cosδfr-Fyfr(t-1)sinδfr]Brn+[Fxfr(t-1)sinδfr+Fyfr(t-1)cosδfr]Lfn+Fxrr(t-1)Brn-Fxrl(t-1)Bln-[Fyrr(t-1)+Fyrl(t-1)]Lrn
上式是中Ts采样时间;
上述式中,xp(k)=[M,Izz]T分别为车辆非线性动力学观测器***的状态矢量和参数矢量,u(k)=[δf,wij,Tij]T和y(k)=[ωz,ax]T分别为车辆非线性动力学观测器***的输入矢量和量测矢量,qs、rs分别为***的过程噪音和量测噪音,两者为***互不相关且均值为零的高斯白噪声;
惯性参数估计***被构造如下:
在上述参数估计***中,qp、rp分别为***的过程噪音和量测噪音,d(k)=[ωz,ax]T为量测矢量;
其中:
同时
∑Fxij(t-1)=Fxfl(t-1)cosδfl-Fyfl(t-1)sinδfl+Fxfr(t-1)cosδfr-Fyfr(t-1)sinδfr+Fxrl(t-1)+Fxrr(t-1)
6.根据权利要求1所述的一种分布式驱动电动汽车行驶状态与惯性参数联合估计方法,其特征在于,在所述步骤S4中,所述双容积卡尔曼滤波观测器具体包括以下步骤:
2)参数变量的时间更新:
2-1)误差协方差矩阵Pp(t-1)进行奇异值分解,并得到xpj(t-1);
2-4)计算协方差预测值Pp(t|t-1);
3)状态变量的时间更新:
3-1)误差协方差矩阵Pp(t-1)进行奇异值分解,并得到xsj(t-1);
3-4)计算协方差预测值Ps(t|t-1);
4)参数变量的量测更新:
4-1)误差协方差矩阵Pp(t|t-1)进行奇异值分解,并得到xpj(t|t-1);
4-2)计算容积点dpj(t|t-1);
4-4)计算新息协方差预测值Ppdd(t|t-1);
4-5)计算互协方差Ppxd(t|t-1);
4-6)计算增益矩阵Kpt;
4-8)计算误差协方差矩阵Pp(t);
5)状态变量的量测更新:
5-1)误差协方差矩阵Ps(t|t-1)进行奇异值分解,并得到xsj(t|t-1);
5-2)计算容积点ysj(t|t-1);
5-4)计算新息协方差预测值Psyy(t|t-1);
5-5)计算互协方差Psxy(t|t-1);
5-6)计算增益矩阵Kst;
5-8)计算误差协方差矩阵Ps(t)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010701509.4A CN111845775B (zh) | 2020-07-20 | 2020-07-20 | 一种分布式驱动电动汽车行驶状态与惯性参数联合估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010701509.4A CN111845775B (zh) | 2020-07-20 | 2020-07-20 | 一种分布式驱动电动汽车行驶状态与惯性参数联合估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111845775A true CN111845775A (zh) | 2020-10-30 |
CN111845775B CN111845775B (zh) | 2022-01-07 |
Family
ID=73001678
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010701509.4A Active CN111845775B (zh) | 2020-07-20 | 2020-07-20 | 一种分布式驱动电动汽车行驶状态与惯性参数联合估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111845775B (zh) |
Cited By (14)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112622876A (zh) * | 2021-03-09 | 2021-04-09 | 天津所托瑞安汽车科技有限公司 | 车辆质心位置确定方法、装置、设备和存储介质 |
CN112622536A (zh) * | 2020-12-25 | 2021-04-09 | 中国农业大学 | 车辆轮胎工作状态监测传感装置及方法 |
CN113008229A (zh) * | 2021-02-26 | 2021-06-22 | 南京理工大学 | 一种基于低成本车载传感器的分布式自主组合导航方法 |
CN113104040A (zh) * | 2021-04-21 | 2021-07-13 | 江苏大学 | 一种考虑观测信息时域衰减的轮胎-路面附着系数获取方法 |
CN113203422A (zh) * | 2021-04-14 | 2021-08-03 | 武汉理工大学 | 一种基于尺寸测量装置的货车状态惯性参数联合估计方法 |
CN113335263A (zh) * | 2021-07-19 | 2021-09-03 | 中国第一汽车股份有限公司 | 一种分布式四驱扭矩控制方法 |
CN113341997A (zh) * | 2021-06-28 | 2021-09-03 | 湖南大学无锡智能控制研究院 | 一种基于多状态参数协同估计的横向控制方法及*** |
CN113408050A (zh) * | 2021-06-07 | 2021-09-17 | 美通重机有限公司 | 轮式装载机行驶稳定性分析方法 |
CN113442906A (zh) * | 2021-06-25 | 2021-09-28 | 江苏大学 | 一种分布式驱动电动汽车横向稳定性分层控制***及方法 |
CN113602279A (zh) * | 2021-08-27 | 2021-11-05 | 吾驾之宝汽车技术(江苏)有限公司 | 分布式驱动电动汽车质心侧偏角与轮胎侧向力估计方法 |
CN113911124A (zh) * | 2021-09-28 | 2022-01-11 | 武汉理工大学 | 基于Bayes的路面附着系数估计方法及电子设备 |
CN114093218A (zh) * | 2021-11-08 | 2022-02-25 | 中国农业大学 | 一种多功能拖拉机驾驶模拟方法及*** |
CN116279523A (zh) * | 2022-12-09 | 2023-06-23 | 聊城大学 | 一种四轮转向驱动电动汽车状态估计方法 |
CN116992697A (zh) * | 2023-09-26 | 2023-11-03 | 华东交通大学 | 一种智能电动汽车行驶状态信息估计方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106250591A (zh) * | 2016-07-21 | 2016-12-21 | 辽宁工业大学 | 一种考虑侧倾影响的汽车行驶状态估计方法 |
CN106515740A (zh) * | 2016-11-14 | 2017-03-22 | 江苏大学 | 基于icdkf的分布式电驱动汽车行驶状态参数估计算法 |
CN108128308A (zh) * | 2017-12-27 | 2018-06-08 | 长沙理工大学 | 一种分布式驱动电动汽车的车辆状态估计***及方法 |
CN108284841A (zh) * | 2017-12-11 | 2018-07-17 | 江苏大学 | 一种分布式驱动电动汽车行驶状态自适应迭代估计方法 |
CN108394413A (zh) * | 2018-01-26 | 2018-08-14 | 辽宁工业大学 | 一种四轮独立驱动与转向的电动汽车状态与参数校正方法 |
CN109849899A (zh) * | 2019-02-01 | 2019-06-07 | 江苏大学 | 一种电动轮汽车电液复合车身稳定控制***及方法 |
-
2020
- 2020-07-20 CN CN202010701509.4A patent/CN111845775B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106250591A (zh) * | 2016-07-21 | 2016-12-21 | 辽宁工业大学 | 一种考虑侧倾影响的汽车行驶状态估计方法 |
CN106515740A (zh) * | 2016-11-14 | 2017-03-22 | 江苏大学 | 基于icdkf的分布式电驱动汽车行驶状态参数估计算法 |
CN108284841A (zh) * | 2017-12-11 | 2018-07-17 | 江苏大学 | 一种分布式驱动电动汽车行驶状态自适应迭代估计方法 |
CN108128308A (zh) * | 2017-12-27 | 2018-06-08 | 长沙理工大学 | 一种分布式驱动电动汽车的车辆状态估计***及方法 |
CN108394413A (zh) * | 2018-01-26 | 2018-08-14 | 辽宁工业大学 | 一种四轮独立驱动与转向的电动汽车状态与参数校正方法 |
CN109849899A (zh) * | 2019-02-01 | 2019-06-07 | 江苏大学 | 一种电动轮汽车电液复合车身稳定控制***及方法 |
Non-Patent Citations (2)
Title |
---|
金贤建: "分布式驱动电动汽车状态参数估计与侧向稳定性鲁棒控制研究", 《中国博士学位论文全文数据库 工程科技Ⅱ辑》 * |
金贤建等: "分布式驱动电动汽车双无迹卡尔曼滤波状态参数联合观测", 《机械工程学报》 * |
Cited By (23)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112622536A (zh) * | 2020-12-25 | 2021-04-09 | 中国农业大学 | 车辆轮胎工作状态监测传感装置及方法 |
CN112622536B (zh) * | 2020-12-25 | 2023-10-31 | 中国农业大学 | 车辆轮胎工作状态监测传感装置及方法 |
CN113008229A (zh) * | 2021-02-26 | 2021-06-22 | 南京理工大学 | 一种基于低成本车载传感器的分布式自主组合导航方法 |
CN113008229B (zh) * | 2021-02-26 | 2024-04-05 | 南京理工大学 | 一种基于低成本车载传感器的分布式自主组合导航方法 |
CN112622876A (zh) * | 2021-03-09 | 2021-04-09 | 天津所托瑞安汽车科技有限公司 | 车辆质心位置确定方法、装置、设备和存储介质 |
CN112622876B (zh) * | 2021-03-09 | 2021-05-14 | 天津所托瑞安汽车科技有限公司 | 车辆质心位置确定方法、装置、设备和存储介质 |
CN113203422A (zh) * | 2021-04-14 | 2021-08-03 | 武汉理工大学 | 一种基于尺寸测量装置的货车状态惯性参数联合估计方法 |
CN113104040A (zh) * | 2021-04-21 | 2021-07-13 | 江苏大学 | 一种考虑观测信息时域衰减的轮胎-路面附着系数获取方法 |
CN113408050A (zh) * | 2021-06-07 | 2021-09-17 | 美通重机有限公司 | 轮式装载机行驶稳定性分析方法 |
CN113442906A (zh) * | 2021-06-25 | 2021-09-28 | 江苏大学 | 一种分布式驱动电动汽车横向稳定性分层控制***及方法 |
CN113341997A (zh) * | 2021-06-28 | 2021-09-03 | 湖南大学无锡智能控制研究院 | 一种基于多状态参数协同估计的横向控制方法及*** |
CN113341997B (zh) * | 2021-06-28 | 2022-04-05 | 湖南大学无锡智能控制研究院 | 一种基于多状态参数协同估计的横向控制方法及*** |
CN113335263B (zh) * | 2021-07-19 | 2022-04-12 | 中国第一汽车股份有限公司 | 一种分布式四驱扭矩控制方法 |
CN113335263A (zh) * | 2021-07-19 | 2021-09-03 | 中国第一汽车股份有限公司 | 一种分布式四驱扭矩控制方法 |
CN113602279A (zh) * | 2021-08-27 | 2021-11-05 | 吾驾之宝汽车技术(江苏)有限公司 | 分布式驱动电动汽车质心侧偏角与轮胎侧向力估计方法 |
CN113911124A (zh) * | 2021-09-28 | 2022-01-11 | 武汉理工大学 | 基于Bayes的路面附着系数估计方法及电子设备 |
CN113911124B (zh) * | 2021-09-28 | 2024-04-30 | 武汉理工大学 | 基于Bayes的路面附着系数估计方法及电子设备 |
CN114093218B (zh) * | 2021-11-08 | 2022-07-26 | 中国农业大学 | 一种多功能拖拉机驾驶模拟方法及*** |
CN114093218A (zh) * | 2021-11-08 | 2022-02-25 | 中国农业大学 | 一种多功能拖拉机驾驶模拟方法及*** |
CN116279523A (zh) * | 2022-12-09 | 2023-06-23 | 聊城大学 | 一种四轮转向驱动电动汽车状态估计方法 |
CN116279523B (zh) * | 2022-12-09 | 2024-03-08 | 金华市壹维科技有限公司 | 一种四轮转向驱动电动汽车状态估计方法 |
CN116992697A (zh) * | 2023-09-26 | 2023-11-03 | 华东交通大学 | 一种智能电动汽车行驶状态信息估计方法 |
CN116992697B (zh) * | 2023-09-26 | 2023-12-15 | 华东交通大学 | 一种智能电动汽车行驶状态信息估计方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111845775B (zh) | 2022-01-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111845775B (zh) | 一种分布式驱动电动汽车行驶状态与惯性参数联合估计方法 | |
CN102556075B (zh) | 一种基于改进扩展卡尔曼滤波的车辆运行状态估计方法 | |
Acosta et al. | Virtual tyre force sensors: An overview of tyre model-based and tyre model-less state estimation techniques | |
Li et al. | Side-slip angle estimation based lateral dynamics control for omni-directional vehicles with optimal steering angle and traction/brake torque distribution | |
CN111547059A (zh) | 一种分布式驱动电动汽车惯性参数估计方法 | |
CN113359457B (zh) | 智能车底盘域控制器的高维动力学模型解算装置及方法 | |
CN111796522B (zh) | 一种车辆状态估计方法 | |
WO2022266824A1 (zh) | 一种转向控制方法及装置 | |
CN113602279B (zh) | 分布式驱动电动汽车质心侧偏角与轮胎侧向力估计方法 | |
CN111959516B (zh) | 一种车辆状态与路面附着系数联合估计的方法 | |
CN114103967A (zh) | 四轮独立驱动电动汽车质心侧偏角与轮胎侧向力估计方法 | |
Michałek et al. | Modular kinematic modelling of articulated buses | |
Wenzel et al. | Kalman filter as a virtual sensor: applied to automotive stability systems | |
Rill | A modified implicit Euler algorithm for solving vehicle dynamic equations | |
Li et al. | Vehicle velocity estimation for real-time dynamic stability control | |
Cheng et al. | A hierarchical estimation scheme of tire-force based on random-walk SCKF for vehicle dynamics control | |
CN115946707A (zh) | 四轮毂电机驱动全线控电动汽车轮胎力估计方法及*** | |
CN115406669A (zh) | 一种多轴特种车辆侧翻指标优化方法 | |
CN116049608A (zh) | 一种非结构性路面下的路面峰值附着系数估计方法及*** | |
Azman et al. | Influence of anti-dive and anti-squat geometry in combined vehicle bounce and pitch dynamics | |
Du et al. | Performance analysis for the model-aided inertial navigation system with different vehicle models | |
Tan et al. | Vehicle state estimation of steer by wire system based on multi sensor fusion | |
CN111469855A (zh) | 一种车辆运动参数的计算方法 | |
Cadiou et al. | Non-linear tyre forces estimation based on vehicle dynamics observation in a finite time | |
Shin et al. | New vehicle dynamics model for yaw rate estimation |
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 |