CN105416276B - 基于高阶滑模的电动汽车稳定性直接横摆力矩控制方法 - Google Patents
基于高阶滑模的电动汽车稳定性直接横摆力矩控制方法 Download PDFInfo
- Publication number
- CN105416276B CN105416276B CN201510922607.XA CN201510922607A CN105416276B CN 105416276 B CN105416276 B CN 105416276B CN 201510922607 A CN201510922607 A CN 201510922607A CN 105416276 B CN105416276 B CN 105416276B
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- mtd
- mtr
- mfrac
- 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
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000012706 support-vector machine Methods 0.000 claims abstract description 11
- 238000011217 control strategy Methods 0.000 claims abstract description 10
- 230000006641 stabilisation Effects 0.000 claims abstract description 8
- 238000011105 stabilization Methods 0.000 claims abstract description 7
- 230000001133 acceleration Effects 0.000 claims abstract description 3
- 238000013461 design Methods 0.000 claims description 20
- 230000008859 change Effects 0.000 claims description 7
- 239000011159 matrix material Substances 0.000 claims description 6
- 230000008569 process Effects 0.000 claims description 6
- 238000012549 training Methods 0.000 claims description 4
- 238000013507 mapping Methods 0.000 claims description 3
- 238000010276 construction Methods 0.000 claims description 2
- 230000005484 gravity Effects 0.000 claims description 2
- 238000005457 optimization Methods 0.000 claims description 2
- 230000008901 benefit Effects 0.000 description 5
- 238000009415 formwork Methods 0.000 description 4
- 230000004044 response Effects 0.000 description 4
- 235000021170 buffet Nutrition 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 206010044565 Tremor Diseases 0.000 description 1
- 238000005299 abrasion Methods 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000003213 activating effect Effects 0.000 description 1
- 230000032683 aging Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 239000004035 construction material Substances 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000010612 desalination reaction Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012417 linear regression Methods 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000035935 pregnancy Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 230000001737 promoting effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000035882 stress Effects 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 230000001960 triggered effect Effects 0.000 description 1
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
- B60W30/00—Purposes of road vehicle drive control systems not related to the control of a particular sub-unit, e.g. of systems using conjoint control of vehicle sub-units
- B60W30/02—Control of vehicle driving stability
- B60W30/045—Improving turning performance
-
- 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
- B60W40/114—Yaw movement
-
- 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
- B60W2540/00—Input parameters relating to occupants
- B60W2540/18—Steering angle
-
- 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
- B60W2720/00—Output or target parameters relating to overall vehicle dynamics
- B60W2720/14—Yaw
Landscapes
- Engineering & Computer Science (AREA)
- Automation & Control Theory (AREA)
- Transportation (AREA)
- Mechanical Engineering (AREA)
- Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Steering Control In Accordance With Driving Conditions (AREA)
Abstract
基于高阶滑模的电动汽车稳定性直接横摆力矩控制方法涉及电动汽车稳定性的控制领域,根据信号采集与调理电路检测到方向盘转角和纵向车速,求得理想横摆角速度值;根据检测到的车辆当前时刻的横摆加速度和实际横摆角速度经基于主动控制与自适应估计的鲁棒滑模观测器获得质心侧偏角估计值;然后以横摆角速度与理想横摆角速度的偏差、实际车辆的质心侧偏角两个参数作为输入变量,采用高阶滑模控制策略,获得满足汽车稳定性的直接横摆力矩;最后以车辆稳定裕度为目标函数和约束条件,利用支持向量机算法进行驱动力或制动力的分配。本发明实现了汽车稳定性直接横摆力矩控制系统的有限时间收敛,提高了汽车高速和恶劣道路等极限条件下的行驶稳定性。
Description
技术领域
本发明涉及电动汽车稳定性的控制领域,特别涉及一种基于高阶滑模的电动汽车稳定性直接横摆力矩控制方法。
背景技术
通常情况下,电动汽车转向操纵都能实现较平稳的转向控制,但遇到极限情况时极易造成车辆进入动力学不稳定状态,出现严重的过多转向或不足转向问题,甚至甩尾。汽车直接横摆力矩控制(以下简称DYC Direct yaw moment control)的开发就是旨在淡化驾驶人员的操作技能对车辆运动安全性的影响,在车辆的各种行驶状态下通过对每个车轮的受力进行调节,克服过多转向或不足转向,从而主动地对车辆进行动力学控制,提高电动汽车在高速和恶劣道路等极限条件下行驶时的操纵稳定性。
目前,横摆力矩控制策略常采用的方法有滑模控制、预测控制和模糊控制等,并取得了不错的效果。但是,由于电动汽车系统中的燃烧及能量转化过程的复杂性、机械部件的磨损和变形特性,轮胎结构材料力学性能的复杂性和非线性、行驶工况的多样性,使得整车稳定性控制系统具有高阶、非线性和强耦合的特征,很难找到精确的机理模型来准确的描述其动力学特性,目前已有的模型大多是忽略很多因素简化得到的,在此基础上设计的控制策略是难以保证整车全工况条件下性能指标最优的;随着设备的老化或环境、负载的作用,设备的性能会随之退化,当累积到一定程度时,导致退化失效,使得原有的基于模型的控制器控制品质下降,随着使用年限和行驶里程的增加,排放浓度和油耗呈现不断增加趋势。因此,需要选用具有高鲁棒性的控制策略,模糊控制是基于专家经验的控制方法,控制品质与专家经验密切相关,很难保证全工况下整车性能最优。如何根据受控系统的动态特性,去选择控制规律,才能使得系统按照一定的技术要求运行,并使得描述系统性能或品质的某个指标在一定的意义下达到最优值,解决最优控制问题的主要方法有古典变分法、极大值原理和动态规划等。
考虑在实际的电动汽车稳定性控制系统设计中很难得到精确、完整的运动模型。因此在建立数学模型时,常做合理的近似处理,忽略一些不确定性的因素,诸如参数误差、未建模动态、观测噪声以及不确定性的外界干扰等等,然而这些不确定性的存在可能会引起控制系统品质恶化,甚至成为系统不稳定的原因。
滑模变结构控制具有快速响应、对参数变化及扰动不灵敏、无需系统在线辨识、实现简单等优点受到各国学者的重视。传统的滑模变结构控制采用线性滑模,系统状态与给定轨迹之间的偏差渐近收敛。与线性滑模相比,终端滑模通过在滑模中有目的的引入非线性项,改善了系统的收敛特性,使得系统状态在有限时间内收敛到给定轨迹。因此终端滑模具有动态响应速度快,有限时间收敛,稳态跟踪精度高等优点,特别适用于高精度的控制。但是这两种滑模变结构控制都存在一个严重的缺点,即抖振。由于抖振很容易激发系统的未建模特性,从而影响了系统的控制性能,给滑模变结构控制的实际应用带来了困难。应运而生的高阶滑模保持了传统滑模的优点,抑制了抖振,消除了相对阶的限制并提高了控制精度,因此高阶滑模理论特别适用于高精度的控制,并且在实际工程中逐渐得到了推广和应用。
发明内容
为了解决电动汽车稳定性直接横摆力矩控制系统在高速和恶劣道路等极限条件下出现的过多转向或转向不足,参数摄动、外界环境条件变化导致所设计的控制策略失效而引发的电动汽车行驶稳定性问题,本发明提出了一种基于高阶滑模的电动汽车稳定性直接横摆力矩控制方法,有效提高了系统动态响应的快速性和鲁棒性,提高了电动汽车高速和恶劣道路等极限条件下的行驶稳定性。
本发明解决技术问题所采取的技术方案如下:
基于高阶滑模的电动汽车稳定性直接横摆力矩控制方法,其方法包括如下步骤:
步骤一、理想横摆角速度计算器根据信号采集与调理电路检测到实际车辆的方向盘转角和纵向车速,求得理想横摆角速度值;
步骤二、根据检测到的实际车辆当前时刻的横摆角速度和步骤一求得的理想横摆角速度值,经基于主动控制与自适应估计的鲁棒观测器获得质心侧偏角估计值;
步骤三、高阶滑模控制器以横摆角速度与理想横摆角速度的偏差、由步骤二获得的实际车辆的质心侧偏角估计值两个参数作为输入变量,采用汽车稳定性直接横摆力矩的高阶滑模控制策略,计算求得满足汽车稳定性的直接横摆力矩;
步骤四、横摆力矩分配器以车辆稳定裕度为目标函数,以步骤三求得的直接横摆力矩限制、车辆电机的输出限制、路面附着条件等为约束条件,利用支持向量机算法计算获得实际车辆四个车轮的驱动力矩或制动力矩,使汽车实际运行路径与期望路径保持一致,实现对汽车稳定性直接横摆力矩的控制。
本发明的有益效果如下:
1)本发明基于主动控制设计的鲁棒观测器,其未知参数的变化率采用自适应估计的方法获得,通过修改设计参数可调整观测器跟踪系统状态的收敛速度,使状态估计达到预期的指标,解决了在极限工况下质心侧偏角较难测量及精确估计的问题,具有鲁棒性强、可靠性高、动态响应速度快、稳态跟踪精度高等优点。
2)本发明基于高阶滑模的电动汽车稳定性直接横摆力矩控制器,以横摆角速度与理想横摆角速度的偏差作为输入变量,基于三阶滑模理论获得维持电动汽车稳定运行所需的直接横摆力矩,解决了系统因建模不精确、参数摄动以及外界环境条件变化时所设计的控制策略失效引发的电动汽车行驶稳定性问题,提高系统鲁棒性的同时降低了滑模控制中固有的抖阵问题,通过调整控制器中的设计参数可实现电动汽车稳定性直接横摆力矩控制系统的有限时间收敛,提高了直接横摆力矩控制系统的快速性和鲁棒性。
3)本发明以车辆稳定裕度为目标函数,以高阶滑模控制器的输出限制、车辆电机输出限制、路面附着条件等为约束条件,利用支持向量机算法实现横摆力矩分配,从而获得四个车轮的驱动力矩或制动力矩,以保证车辆的稳定运行。
4)本发明基于质心侧偏角观测的电动汽车稳定系统直接横摆力矩控制方法提高了电动汽车在高速和恶劣道路等极限条件下的行驶稳定性,克服行驶过程中出现的过多转向或转向不足的问题。
5)本发明的方法简单易于实现,适宜广泛推广应用。
附图说明
图1是本发明基于高阶滑模的电动汽车稳定性直接横摆力矩控制方法的原理图。
具体实施方式
下面结合附图和具体实施例对本发明作进一步详细说明。
如图1所示,本发明基于高阶滑模的电动汽车稳定性直接横摆力矩控制方法由信号采集与调理电路2、鲁棒观测器3、理想横摆角速度计算器4、高阶滑模控制器5以及横摆力矩分配器6实现。基于主动控制的鲁棒观测器3的输入为信号采集与调理电路2检测到的实际车辆1方向盘转角和纵向车速,设计主动控制策略,未知参数的变化率采用自适应估计的方法获得,通过修改设计参数可调整观测器跟踪系统状态的收敛速度,使状态估计达到预期的指标,以获取准确的质心侧偏角估计值。通过理想横摆角速度计算器4确定理想横摆角速度,信号采集与调理电路2检测到驾驶员输入的方向盘转角δ和车速vx,经过理想横摆角速度计算器4计算出理想横摆角速度rd。基于高阶滑模的电动汽车稳定性直接横摆力矩控制器即高阶滑模控制器5以横摆角速度与理想横摆角速度的偏差作为输入变量,基于三阶滑模理论获得维持电动汽车稳定运行所需的直接横摆力矩。横摆力矩分配器6采取不同工况下汽车稳定性的直接横摆力矩分配策略,汽车直线行驶时采用平均分配方式,即在车辆行驶过程中将调整力矩和驱动力平均分配给各个车轮;在车辆转向工况运行过程中,以车辆稳定裕度为目标函数,以高阶滑模控制器5的输出限制、车辆电机输出限制、路面附着条件等为约束条件,利用支持向量机算法实现横摆力矩分配,使四个车轮获得最优的驱动力矩或制动力矩。
本发明基于高阶滑模的电动汽车稳定性直接横摆力矩控制方法具体实施步骤如下:
1)鲁棒观测器3
信号采集与调理电路2检测到实际车辆1中的侧向加速度和实际横摆角速度,经基于主动控制的鲁棒观测器获得质心侧偏角估计值。
算法如下:
考虑实际应用中,总是存在或多或少的诸如建模误差、噪声干扰等不确定因素,因此根据牛顿定律,建立包含参数不确定以及干扰、噪声影响的线性二自由度车辆动力学方程
y=Cx+Dδ (1)式中,
其中,β为车辆质心侧偏角,r为横摆角速度,m和Iz分别为车辆质量和横摆转动惯量,a和b为前轴和后轴到质心的距离,vx为车辆纵向速度,δ为前轮转角,Ci(i=F/R前/后)为轮胎侧偏刚度。未知函数f=ΔA+ΔB1δ+ξ,满足有界。其中,ΔA和ΔB1参数不确定项,ξ为考虑干扰、噪声影响的不确定函数项。
基于主动控制的思想设计如下形式的鲁棒观测器:
式中,为观测器估计的系统状态,G为设计参数矩阵,A0=A-GC,其中设计参数矩阵G满足A0为Hurwitz矩阵,v为观测器的控制输入。
式中,设计参数矩阵η0和τ0的选取满足其特征根为负。m=3,n=5,式(1)中f的估计值的变化率采用自适应估计的方法获得,即ζ>0为设计参数。
则鲁棒观测器3可以渐近估计出式(1)系统的车辆质心侧偏角。
2)理想横摆角速度计算器4
信号采集与调理电路2检测到实际车辆1中驾驶员输入的方向盘转角δ和车速vx,由理想横摆角速度计算器4确定理想横摆角速度rd。算法如下:
式中,K为稳定性系数,K=[m/(a+b)2][a/CR-b/CF];μ为路面附着系数;g为重力加速度,sgn为符号函数。其他参数的含义与式(1)中的定义相同。
3)高阶滑模控制器5
基于高阶滑模的直接横摆力矩控制器设计,在二自由度模型上增加一个附加的横摆力矩Tz,采用电动汽车稳定性直接横摆力矩控制系统的3阶滑模策略,在提高系统鲁棒性的同时降低甚至消除横摆力矩Tz的抖振。
算法如下:
考虑到建模误差以及参数摄动,在二自由度模型上增加一个附加的横摆力矩Tz,式(2)可修改为如下形式:
式中,B2=[1,0]T,Tz为横摆力矩。
设计直接横摆力矩Tz的高阶滑模控制策略:
Tz=Tzeq+Tzn (6)
式中,
Tzn由滤波器获得:
式中,e=rd-r,k、η2、ε、ρ、η为设计参数,满足k>0,η2>0,ε>0,ρ>0,η>0,其中η、ε满足[0,1;-η/ε,-1/ε]的特征值为负值,p,q为奇数,且1<p/q<2。
4)横摆力矩分配器6
高阶滑模控制器5的输出为电动汽车稳定性的直接横摆力矩,以车辆稳定裕度为指标,经基于支持向量机的横摆力矩分配器6计算获得四个车轮的驱动力矩或制动力矩。
(1)当直线工况下,电动汽车直线行驶时采用平均分配方式,即在车辆行驶过程中将调整力矩和驱动力平均分配给各个车轮,分配公式如式(9)所示。
式中Ti(i=1,2,3,4)分别表示左前轮、左后轮、右前轮、右后轮的驱动力力矩;Tr表示总的驱动力矩,由踏板开度决定,Tz表示高阶滑模控制器计算求得的横摆力矩。
(2)在转向极限工况时,考虑路面附着条件、轴荷转移、电机输出限制的情况下纵向力的变化对车辆转向稳定性的影响,对高阶滑模控制器5的输出即电动汽车稳定性的直接横摆力矩,以车辆稳定裕度为指标,采用支持向量机的方法实现转矩的分配,具体实现方法如下:
当前轮转角较小时,近似有
w=Bu
式中,w=[Fx,Tz],Fx表示车辆总的横摆力;u=[Fx1,Fx2,Fx3,Fx4]T,Fx1、Fx2、Fx3和Fx4分别表示各轮胎的横摆力;lt表示车轴的轮间距,R表示车轮半径。
支持向量机在解决小样本、非线性及高维模式识别中表现出许多特有的优势,它通过一个非线性映射将数据映射到一个高维特征空间,并在此空间进行线性回归。本发明采用支持向量机算法对各轮胎的横摆力进行分配:
式中,Fxi表示各轮胎的横摆力,i=1,2,3,4,是从输入空间到高维特征空间的非线性映射向量,本发明选择为径向基核函数,向量ωs和偏置b是待定参数,输入量x为四个车轮的实际角速度ω1、ω2、ω3、ω4与理想角速度的差值及变化率即
通过对电动汽车的动力学原理分析,轮胎利用率越低则车辆的稳定性裕度越高。因此本发明利用表征车辆的稳定性裕度的轮胎利用率来确定目标函数和约束条件。将4个轮胎利用率平方和的最小值即车辆的稳定性裕度最大值作为车辆稳定性的目标函数,来对各轮胎应受力矩进行分配,以保证车辆的稳定运行。在只考虑优化横摆力分配、电机输出限制、路面附着条件限制,利用结构风险最小化准则(SRM)构造最小化目标函数:
式中,em为训练数据精度,l为训练次数,m为迭代次数,C为容量因子,ε权重系数,用于协调目标函数中各项的作用,Vu=diag(1/(μFzi)2),其中μ为路面附着系数,Fzi为四个车轮各轮胎垂直方向的力。
定义Lagrange函数为:
式中,αm为Lagrange乘子。
根据KTT条件:
对于m=1,…,l,消去ω和em,得到如下方程:
式中,el=[e1...el]T,m,j=1,...l,FxiL=[Fxi1,...,Fxim,...,Fxil]T,α=[α1...αl]T。
根据最小二乘法求出αm和b,由此得到各轮胎的横摆力为:
Claims (5)
1.基于高阶滑模的电动汽车稳定性直接横摆力矩控制方法,其特征在于,该方法包括如下步骤:
步骤一、理想横摆角速度计算器(4)根据信号采集与调理电路(2)检测到实际车辆(1)的方向盘转角和纵向车速,求得理想横摆角速度值;
步骤二、根据检测到的实际车辆(1)当前时刻的横摆角速度和步骤一求得的理想横摆角速度值,经基于主动控制与自适应估计的鲁棒观测器(3)获得质心侧偏角估计值;
步骤三、高阶滑模控制器(5)以横摆角速度与理想横摆角速度的偏差、由步骤二获得的实际车辆(1)的质心侧偏角估计值两个参数作为输入变量,采用汽车稳定性直接横摆力矩的高阶滑模控制策略,计算求得满足汽车稳定性的直接横摆力矩;
步骤四、横摆力矩分配器(6)以车辆稳定裕度为目标函数,以步骤三求得的直接横摆力矩限制、车辆电机的输出限制、路面附着条件为约束条件,利用支持向量机算法计算获得实际车辆(1)四个车轮的驱动力矩或制动力矩,使汽车实际运行路径与期望路径保持一致,实现对汽车稳定性直接横摆力矩的控制。
2.如权利要求1所述的基于高阶滑模的电动汽车稳定性直接横摆力矩控制方法,其特征在于,步骤二所述的鲁棒观测器(3)获得质心侧偏角估计值的过程如下:
根据牛顿定律,建立包含参数不确定以及干扰、噪声影响的线性二自由度车辆动力学方程如下:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<mi>A</mi>
<mi>x</mi>
<mo>+</mo>
<msub>
<mi>B</mi>
<mn>1</mn>
</msub>
<mi>&delta;</mi>
<mo>+</mo>
<mi>f</mi>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>y</mi>
<mo>=</mo>
<mi>C</mi>
<mi>x</mi>
<mo>+</mo>
<mi>D</mi>
<mi>&delta;</mi>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,
<mrow>
<mi>x</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>&beta;</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>r</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mi>A</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>A</mi>
<mn>11</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>A</mi>
<mn>12</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>A</mi>
<mn>21</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>A</mi>
<mn>22</mn>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mn>2</mn>
<mfrac>
<mrow>
<msub>
<mi>C</mi>
<mi>F</mi>
</msub>
<mo>+</mo>
<msub>
<mi>C</mi>
<mi>R</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>v</mi>
<mi>x</mi>
</msub>
<mi>m</mi>
</mrow>
</mfrac>
</mrow>
</mtd>
<mtd>
<mrow>
<mn>2</mn>
<mfrac>
<mrow>
<mo>-</mo>
<msub>
<mi>aC</mi>
<mi>F</mi>
</msub>
<mo>+</mo>
<msub>
<mi>bC</mi>
<mi>R</mi>
</msub>
</mrow>
<mrow>
<msup>
<msub>
<mi>v</mi>
<mi>x</mi>
</msub>
<mn>2</mn>
</msup>
<mi>m</mi>
</mrow>
</mfrac>
<mo>-</mo>
<mn>1</mn>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mn>2</mn>
<mfrac>
<mrow>
<mo>-</mo>
<msub>
<mi>aC</mi>
<mi>F</mi>
</msub>
<mo>+</mo>
<msub>
<mi>bC</mi>
<mi>R</mi>
</msub>
</mrow>
<msub>
<mi>I</mi>
<mi>Z</mi>
</msub>
</mfrac>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mn>2</mn>
<mfrac>
<mrow>
<msup>
<mi>a</mi>
<mn>2</mn>
</msup>
<msub>
<mi>C</mi>
<mi>F</mi>
</msub>
<mo>+</mo>
<msup>
<mi>b</mi>
<mn>2</mn>
</msup>
<msub>
<mi>C</mi>
<mi>R</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>v</mi>
<mi>x</mi>
</msub>
<msub>
<mi>I</mi>
<mi>Z</mi>
</msub>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<msub>
<mi>B</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>B</mi>
<mn>11</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>B</mi>
<mn>12</mn>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mfrac>
<mrow>
<mn>2</mn>
<msub>
<mi>C</mi>
<mi>F</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>v</mi>
<mi>x</mi>
</msub>
<mi>m</mi>
</mrow>
</mfrac>
</mtd>
</mtr>
<mtr>
<mtd>
<mfrac>
<mrow>
<mn>2</mn>
<msub>
<mi>C</mi>
<mi>F</mi>
</msub>
<mi>a</mi>
</mrow>
<msub>
<mi>I</mi>
<mi>Z</mi>
</msub>
</mfrac>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
<mrow>
<mi>y</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>&alpha;</mi>
<mi>y</mi>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>r</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mi>C</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>C</mi>
<mn>11</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>C</mi>
<mn>12</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>C</mi>
<mn>21</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>C</mi>
<mn>22</mn>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mn>2</mn>
<mfrac>
<mrow>
<msub>
<mi>C</mi>
<mi>F</mi>
</msub>
<mo>+</mo>
<msub>
<mi>C</mi>
<mi>R</mi>
</msub>
</mrow>
<mi>m</mi>
</mfrac>
</mrow>
</mtd>
<mtd>
<mrow>
<mn>2</mn>
<mfrac>
<mrow>
<mo>-</mo>
<msub>
<mi>aC</mi>
<mi>F</mi>
</msub>
<mo>+</mo>
<msub>
<mi>bC</mi>
<mi>R</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>v</mi>
<mi>x</mi>
</msub>
<mi>m</mi>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mi>D</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>D</mi>
<mn>11</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>D</mi>
<mn>12</mn>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mfrac>
<mrow>
<mn>2</mn>
<msub>
<mi>C</mi>
<mi>F</mi>
</msub>
</mrow>
<mi>m</mi>
</mfrac>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
</mrow>
其中,β为车辆质心侧偏角,r为横摆角速度,m为车辆质量,Iz为横摆转动惯量,a为前轴到质心的距离,b为后轴到质心的距离,vx为车辆纵向速度,δ为前轮转角,Ci为轮胎侧偏刚度,i=F或R,即前或后;未知函数f=ΔA+ΔB1δ+ξ,满足有界;其中,ΔA和ΔB1为参数不确定项,ξ为考虑干扰、噪声影响的不确定函数项;
基于主动控制的思想设计如下形式的鲁棒观测器:
<mrow>
<mover>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<mi>A</mi>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mo>+</mo>
<msub>
<mi>B</mi>
<mn>1</mn>
</msub>
<mi>&delta;</mi>
<mo>-</mo>
<mi>G</mi>
<mrow>
<mo>(</mo>
<mi>C</mi>
<mover>
<mi>x</mi>
<mo>^</mo>
</mover>
<mo>+</mo>
<mi>D</mi>
<mi>&delta;</mi>
<mo>-</mo>
<mi>y</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mi>v</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,为观测器估计的系统状态,G为设计参数矩阵,A0=A-GC,其中设计参数矩阵G满足A0为Hurwitz矩阵,v为观测器的控制输入,如下:
<mrow>
<mi>v</mi>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>A</mi>
<mn>0</mn>
</msub>
<mi>e</mi>
<mo>-</mo>
<msub>
<mi>&eta;</mi>
<mn>0</mn>
</msub>
<mi>e</mi>
<mo>+</mo>
<mover>
<mi>f</mi>
<mo>^</mo>
</mover>
<mo>-</mo>
<msub>
<mi>&tau;</mi>
<mn>0</mn>
</msub>
<msup>
<mi>e</mi>
<mrow>
<mi>m</mi>
<mo>/</mo>
<mi>n</mi>
</mrow>
</msup>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>3</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,设计参数矩阵η0和τ0的选取满足其特征根为负;m=3,n=5,式(1)中f的估计值的变化率采用自适应估计的方法获得,即ζ>0为设计参数,则鲁棒观测器(3)可以渐近估计出式(1)系统的车辆质心侧偏角。
3.如权利要求2所述的基于高阶滑模的电动汽车稳定性直接横摆力矩控制方法,其特征在于,步骤一所述的理想横摆角速度计算器(4)求得理想横摆角速度值的过程如下:
信号采集与调理电路(2)检测到实际车辆(1)中驾驶员输入的方向盘转角δ和车速vx,由理想横摆角速度计算模器(4)确定理想横摆角速度rd,算法如下:
<mrow>
<msub>
<mi>r</mi>
<mi>d</mi>
</msub>
<mo>=</mo>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
<mfenced open = "{" close = "}">
<mtable>
<mtr>
<mtd>
<mo>|</mo>
<mfrac>
<mrow>
<msub>
<mi>v</mi>
<mi>x</mi>
</msub>
<mo>/</mo>
<mrow>
<mo>(</mo>
<mi>a</mi>
<mo>+</mo>
<mi>b</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mn>1</mn>
<mo>+</mo>
<msup>
<msub>
<mi>Kv</mi>
<mi>x</mi>
</msub>
<mn>2</mn>
</msup>
</mrow>
</mfrac>
<mo>|</mo>
<mi>&delta;</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mo>|</mo>
<mn>0.85</mn>
<mfrac>
<mrow>
<mi>&mu;</mi>
<mo>&CenterDot;</mo>
<mi>g</mi>
</mrow>
<msub>
<mi>v</mi>
<mi>x</mi>
</msub>
</mfrac>
<mo>|</mo>
<mi>sgn</mi>
<mo>(</mo>
<mi>&delta;</mi>
<mo>)</mo>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>4</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,K为稳定性系数,K=[m/(a+b)2][a/CR-b/CF];μ为路面附着系数;g为重力加速度,sgn为符号函数,a为前轴到质心的距离,b为后轴到质心的距离。
4.如权利要求2所述的基于高阶滑模的电动汽车稳定性直接横摆力矩控制方法,其特征在于,步骤三所述的高阶滑模控制器(5)计算求得满足汽车稳定性的直接横摆力矩的过程如下:
考虑到建模误差以及参数摄动,在二自由度模型上增加一个附加的横摆力矩Tz,式(2)可修改为如下形式:
<mrow>
<mover>
<mi>x</mi>
<mo>&CenterDot;</mo>
</mover>
<mo>=</mo>
<mi>A</mi>
<mi>x</mi>
<mo>+</mo>
<msub>
<mi>B</mi>
<mn>1</mn>
</msub>
<mi>&delta;</mi>
<mo>+</mo>
<msub>
<mi>B</mi>
<mn>2</mn>
</msub>
<msub>
<mi>T</mi>
<mi>Z</mi>
</msub>
<mo>+</mo>
<mi>f</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>5</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,B2=[1,0]T,Tz为横摆力矩;
设计直接横摆力矩Tz的高阶滑模控制策略,如下:
Tz=Tzeq+Tzn (6)
式中,
<mrow>
<msub>
<mi>T</mi>
<mrow>
<mi>z</mi>
<mi>e</mi>
<mi>q</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mover>
<mi>r</mi>
<mo>&CenterDot;</mo>
</mover>
<mi>d</mi>
</msub>
<mo>-</mo>
<msub>
<mi>A</mi>
<mn>21</mn>
</msub>
<mover>
<mi>&beta;</mi>
<mo>^</mo>
</mover>
<mo>-</mo>
<msub>
<mi>A</mi>
<mn>22</mn>
</msub>
<mi>r</mi>
<mo>-</mo>
<msub>
<mi>B</mi>
<mn>12</mn>
</msub>
<mi>&delta;</mi>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>7</mn>
<mo>)</mo>
</mrow>
</mrow>
Tzn由滤波器获得:
<mrow>
<msub>
<mi>&tau;T</mi>
<mrow>
<mi>z</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mover>
<mi>T</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mi>z</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>+</mo>
<mi>&epsiv;</mi>
<msub>
<mover>
<mi>T</mi>
<mo>&CenterDot;&CenterDot;</mo>
</mover>
<mrow>
<mi>z</mi>
<mi>n</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mi>&rho;</mi>
</mfrac>
<mfrac>
<mi>q</mi>
<mi>p</mi>
</mfrac>
<msup>
<mover>
<mi>s</mi>
<mo>&CenterDot;</mo>
</mover>
<mrow>
<mn>2</mn>
<mo>-</mo>
<mi>p</mi>
<mo>/</mo>
<mi>q</mi>
</mrow>
</msup>
<mo>+</mo>
<mrow>
<mo>(</mo>
<mi>k</mi>
<mo>+</mo>
<msub>
<mi>&eta;</mi>
<mn>2</mn>
</msub>
<mo>)</mo>
</mrow>
<mi>sgn</mi>
<mrow>
<mo>(</mo>
<mi>l</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>8</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,e=rd-r,k、η2、ε、ρ、η为设计参数,满足k>0,η2>0,ε>0,ρ>0,η>0,其中η、ε满足[0,1;-η/ε,-1/ε]的特征值为负值,p,q为奇数,且1<p/q<2。
5.如权利要求1所述的基于高阶滑模的电动汽车稳定性直接横摆力矩控制方法,其特征在于,步骤四所述的横摆力矩分配器(6)计算获得实际车辆(1)四个车轮的驱动力矩或制动力矩的过程如下:
当直线工况下,电动汽车直线行驶时采用平均分配方式,即在车辆行驶过程中将调整力矩和驱动力平均分配给各个车轮,分配公式如式(9)所示:
<mrow>
<msub>
<mi>T</mi>
<mn>1</mn>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>T</mi>
<mi>r</mi>
</msub>
<mn>4</mn>
</mfrac>
<mo>-</mo>
<mfrac>
<msub>
<mi>T</mi>
<mi>z</mi>
</msub>
<mn>4</mn>
</mfrac>
<mo>,</mo>
<msub>
<mi>T</mi>
<mn>2</mn>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>T</mi>
<mi>r</mi>
</msub>
<mn>4</mn>
</mfrac>
<mo>+</mo>
<mfrac>
<msub>
<mi>T</mi>
<mi>z</mi>
</msub>
<mn>4</mn>
</mfrac>
<mo>,</mo>
<msub>
<mi>T</mi>
<mn>3</mn>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>T</mi>
<mi>r</mi>
</msub>
<mn>4</mn>
</mfrac>
<mo>-</mo>
<mfrac>
<msub>
<mi>T</mi>
<mi>z</mi>
</msub>
<mn>4</mn>
</mfrac>
<mo>,</mo>
<msub>
<mi>T</mi>
<mn>4</mn>
</msub>
<mo>=</mo>
<mfrac>
<msub>
<mi>T</mi>
<mi>r</mi>
</msub>
<mn>4</mn>
</mfrac>
<mo>+</mo>
<mfrac>
<msub>
<mi>T</mi>
<mi>z</mi>
</msub>
<mn>4</mn>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
式中:T1表示左前轮的驱动力力矩,T2表示左后轮的驱动力力矩,T3表示右前轮的驱动力力矩,T4表示右后轮的驱动力力矩;Tr表示总的驱动力矩,由踏板开度决定,Tz表示高阶滑模控制器计算求得的横摆力矩;
在转向极限工况时,对高阶滑模控制器(5)计算求得的电动汽车稳定性的直接横摆力矩,以车辆稳定裕度为指标,采用支持向量机的方法实现转矩的分配,具体实现方法如下:
当前轮转角较小时,近似有
w=Bu
<mrow>
<mi>B</mi>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mfrac>
<mn>1</mn>
<mi>R</mi>
</mfrac>
</mtd>
<mtd>
<mfrac>
<mn>1</mn>
<mi>R</mi>
</mfrac>
</mtd>
<mtd>
<mfrac>
<mn>1</mn>
<mi>R</mi>
</mfrac>
</mtd>
<mtd>
<mfrac>
<mn>1</mn>
<mi>R</mi>
</mfrac>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>-</mo>
<mfrac>
<msub>
<mi>l</mi>
<mi>t</mi>
</msub>
<mrow>
<mn>2</mn>
<mi>R</mi>
</mrow>
</mfrac>
</mrow>
</mtd>
<mtd>
<mfrac>
<msub>
<mi>l</mi>
<mi>t</mi>
</msub>
<mrow>
<mn>2</mn>
<mi>R</mi>
</mrow>
</mfrac>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mfrac>
<msub>
<mi>l</mi>
<mi>t</mi>
</msub>
<mrow>
<mn>2</mn>
<mi>R</mi>
</mrow>
</mfrac>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>-</mo>
<mfrac>
<msub>
<mi>l</mi>
<mi>t</mi>
</msub>
<mrow>
<mn>2</mn>
<mi>R</mi>
</mrow>
</mfrac>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,w=[Fx,Tz],Fx表示车辆总的横摆力;u=[Fx1,Fx2,Fx3,Fx4]T,Fx1、Fx2、Fx3和Fx4分别表示各轮胎的横摆力;lt表示车轴的轮间距,R表示车轮半径;
采用支持向量机算法对各轮胎的横摆力进行分配,如下:
式中,Fxi表示各轮胎的横摆力,i=1,2,3,4,是从输入空间到高维特征空间的非线性映射向量,选择为径向基核函数,向量ωs和偏置b是待定参数,输入量x为四个车轮的实际角速度ω1、ω2、ω3、ω4与理想角速度的差值及变化率即
将4个轮胎利用率平方和的最小值即车辆的稳定性裕度最大值作为车辆稳定性的目标函数,来对各轮胎应受力矩进行分配;在只考虑优化横摆力分配、电机输出限制、路面附着条件限制,利用结构风险最小化准则构造最小化目标函数:
<mrow>
<mtable>
<mtr>
<mtd>
<munder>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
<mrow>
<mi>&omega;</mi>
<mo>,</mo>
<mi>b</mi>
<mo>,</mo>
<mi>e</mi>
</mrow>
</munder>
</mtd>
<mtd>
<mrow>
<msup>
<mi>J</mi>
<mo>*</mo>
</msup>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mn>2</mn>
</mfrac>
<mrow>
<mo>(</mo>
<msup>
<msub>
<mi>&omega;</mi>
<mi>s</mi>
</msub>
<mi>T</mi>
</msup>
<msub>
<mi>&omega;</mi>
<mi>s</mi>
</msub>
<mo>+</mo>
<msup>
<mi>u</mi>
<mi>T</mi>
</msup>
<msub>
<mi>V</mi>
<mi>u</mi>
</msub>
<mi>u</mi>
<mo>+</mo>
<mi>&epsiv;</mi>
<msup>
<mrow>
<mo>(</mo>
<mrow>
<mi>B</mi>
<mi>u</mi>
<mo>-</mo>
<mi>w</mi>
</mrow>
<mo>)</mo>
</mrow>
<mi>T</mi>
</msup>
<msub>
<mi>V</mi>
<mi>u</mi>
</msub>
<mo>(</mo>
<mrow>
<mi>B</mi>
<mi>u</mi>
<mo>-</mo>
<mi>w</mi>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<mi>C</mi>
<mn>2</mn>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>m</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>l</mi>
</munderover>
<msubsup>
<mi>e</mi>
<mi>m</mi>
<mn>2</mn>
</msubsup>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>13</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,em为训练数据精度,l为训练次数,m为迭代次数,C为容量因子,ε为权重系数,Vu=diag(1/(μFzi)2),其中μ为路面附着系数,Fzi为四个车轮各轮胎垂直方向的力;
定义Lagrange函数为:
式中,αm为Lagrange乘子;
根据KTT条件:
对于m=1,…,l,消去ω和em,得到如下方程:
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mrow>
<msup>
<msub>
<mi>e</mi>
<mi>l</mi>
</msub>
<mi>T</mi>
</msup>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>e</mi>
<mi>l</mi>
</msub>
</mtd>
<mtd>
<mrow>
<mi>Q</mi>
<mo>+</mo>
<msup>
<mi>C</mi>
<mrow>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msup>
<mi>I</mi>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>b</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>&alpha;</mi>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>F</mi>
<mrow>
<mi>x</mi>
<mi>i</mi>
<mi>l</mi>
</mrow>
</msub>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>16</mn>
<mo>)</mo>
</mrow>
</mrow>
式中,el=[e1...el]T,m,j=1,...l,FxiL=[Fxi1,...,Fxim,...,Fxil]T,α=[α1...αl]T;
根据最小二乘法求出αm和b,由此得到各轮胎的横摆力为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510922607.XA CN105416276B (zh) | 2015-12-14 | 2015-12-14 | 基于高阶滑模的电动汽车稳定性直接横摆力矩控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510922607.XA CN105416276B (zh) | 2015-12-14 | 2015-12-14 | 基于高阶滑模的电动汽车稳定性直接横摆力矩控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105416276A CN105416276A (zh) | 2016-03-23 |
CN105416276B true CN105416276B (zh) | 2017-09-22 |
Family
ID=55494957
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510922607.XA Active CN105416276B (zh) | 2015-12-14 | 2015-12-14 | 基于高阶滑模的电动汽车稳定性直接横摆力矩控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105416276B (zh) |
Families Citing this family (28)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105835889B (zh) * | 2016-04-12 | 2018-02-27 | 江苏大学 | 一种基于二阶滑模观测器的车辆质心侧偏角的估计方法 |
CN107662595B (zh) * | 2016-07-30 | 2020-07-10 | 比亚迪股份有限公司 | 一种车辆制动稳定控制方法及系统 |
CN106347361B (zh) * | 2016-10-19 | 2018-12-11 | 长春工业大学 | 一种冗余驱动车辆动力学控制分配方法 |
CN106502091B (zh) * | 2016-10-19 | 2019-06-07 | 长春工业大学 | 一种车辆动力学控制的优化分配方法 |
CN106527139B (zh) * | 2016-11-29 | 2019-07-26 | 桂林航天工业学院 | 通信受限情况下的车辆横摆运动鲁棒控制器设计方法 |
CN107415939B (zh) * | 2017-03-17 | 2020-01-24 | 江苏大学 | 一种分布式驱动电动汽车转向稳定性控制方法 |
CN107253453B (zh) * | 2017-07-05 | 2019-08-27 | 厦门大学 | 一种分布式电动汽车横向稳定性自适应控制系统及方法 |
CN108099900B (zh) * | 2017-12-18 | 2019-09-03 | 长春工业大学 | 一种极限工况下保持汽车横向稳定的四轮转向控制方法 |
CN108099902B (zh) * | 2017-12-18 | 2019-08-30 | 长春工业大学 | 一种体现汽车非线性特性的横摆稳定性控制方法 |
CN108107732B (zh) * | 2017-12-18 | 2020-02-11 | 长春工业大学 | 主动前轮转向和直接横摆力矩联合的汽车稳定性控制方法 |
CN108227491B (zh) * | 2017-12-28 | 2021-11-16 | 重庆邮电大学 | 一种基于滑模神经网络的智能车轨迹跟踪控制方法 |
CN108482363B (zh) * | 2018-04-09 | 2020-03-17 | 吉林大学 | 车辆横摆稳定性预测模型控制方法 |
CN108762083B (zh) * | 2018-06-13 | 2021-04-02 | 长春萨米特光电科技有限公司 | 一种基于加速度观测器的自动控制系统 |
CN109017804B (zh) * | 2018-07-31 | 2022-01-18 | 大连民族大学 | 力矩分配控制器为车辆各轮毂电机分配驱动力矩的方法 |
CN109017759B (zh) * | 2018-07-31 | 2022-04-08 | 大连民族大学 | 期望路径车辆横摆控制方法 |
CN109291932B (zh) * | 2018-10-16 | 2021-03-12 | 浙江东方机电有限公司 | 基于反馈的电动汽车横摆稳定性实时控制装置及方法 |
CN109941246B (zh) * | 2018-11-13 | 2020-07-03 | 清华大学 | 一种集成式线控液压制动系统及其车辆稳定性控制方法 |
CN109522666B (zh) * | 2018-11-27 | 2023-07-14 | 上海埃维汽车技术股份有限公司 | 一种分布式电动汽车稳定性控制方法 |
CN110096077B (zh) * | 2019-05-31 | 2022-01-25 | 西南石油大学 | 开关磁阻电机非奇异快速终端滑模转速控制方法及系统 |
CN110271534B (zh) * | 2019-06-14 | 2021-05-28 | 百度在线网络技术(北京)有限公司 | 自动驾驶车辆的控制方法、装置、计算机设备和存储介质 |
CN110341714B (zh) * | 2019-06-26 | 2021-02-12 | 江苏大学 | 一种同时估计车辆质心侧偏角和扰动的方法 |
CN111055921B (zh) * | 2019-12-31 | 2021-06-29 | 吉林大学 | 一种基于数据驱动的四轮转向模型预测控制方法 |
CN111703417B (zh) * | 2020-06-24 | 2023-09-05 | 湖北汽车工业学院 | 一种高低速统一预瞄滑膜驾驶控制方法及控制系统 |
CN112590737B (zh) * | 2020-12-22 | 2023-12-22 | 安徽卡思普智能科技有限公司 | 一种直接横摆力矩控制器的设计方法 |
CN113978263A (zh) * | 2021-11-02 | 2022-01-28 | 长春工业大学 | 一种驱动轮防滑与转矩优化融合的电动汽车稳定性控制方法 |
CN113830075B (zh) * | 2021-11-30 | 2022-03-11 | 天津所托瑞安汽车科技有限公司 | 车辆稳定控制方法、装置、电子设备及介质 |
CN114684109A (zh) * | 2022-02-07 | 2022-07-01 | 北京新能源汽车股份有限公司 | 一种电动汽车横摆稳定性控制方法及装置 |
CN117681858B (zh) * | 2024-02-04 | 2024-05-17 | 中国第一汽车股份有限公司 | 车辆侧向的安全控制方法、存储介质及车辆 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101783563A (zh) * | 2010-03-10 | 2010-07-21 | 江苏大学 | 电动汽车用多相容错多电机及其解耦控制方法 |
CN102267460A (zh) * | 2011-05-26 | 2011-12-07 | 上海理工大学 | 一种基于轮胎垂向载荷分配的车辆稳定性控制方法 |
CN102407846A (zh) * | 2011-06-14 | 2012-04-11 | 上海理工大学 | 基于轮胎垂向载荷分配的车辆稳定性控制方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR101309508B1 (ko) * | 2011-10-10 | 2013-09-24 | 주식회사 만도 | 차량의 차선 변경 보조 시스템 및 그 방법 |
TWI471239B (zh) * | 2011-10-21 | 2015-02-01 | Automotive Res & Testing Ct | Vehicle Yaw Stability Control Method and Its System |
-
2015
- 2015-12-14 CN CN201510922607.XA patent/CN105416276B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101783563A (zh) * | 2010-03-10 | 2010-07-21 | 江苏大学 | 电动汽车用多相容错多电机及其解耦控制方法 |
CN102267460A (zh) * | 2011-05-26 | 2011-12-07 | 上海理工大学 | 一种基于轮胎垂向载荷分配的车辆稳定性控制方法 |
CN102407846A (zh) * | 2011-06-14 | 2012-04-11 | 上海理工大学 | 基于轮胎垂向载荷分配的车辆稳定性控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN105416276A (zh) | 2016-03-23 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105416276B (zh) | 基于高阶滑模的电动汽车稳定性直接横摆力矩控制方法 | |
Xiong et al. | Vehicle dynamics control of four in-wheel motor drive electric vehicle using gain scheduling based on tyre cornering stiffness estimation | |
Ni et al. | Dynamics control of autonomous vehicle at driving limits and experiment on an autonomous formula racing car | |
Ren et al. | MPC-based yaw stability control in in-wheel-motored EV via active front steering and motor torque distribution | |
Wang et al. | Integrated optimal dynamics control of 4WD4WS electric ground vehicle with tire-road frictional coefficient estimation | |
Zhu et al. | Braking/steering coordination control for in-wheel motor drive electric vehicles based on nonlinear model predictive control | |
Song et al. | Direct‐yaw‐moment control of four‐wheel‐drive electrical vehicle based on lateral tyre–road forces and sideslip angle observer | |
Novi et al. | Real-time control for at-limit handling driving on a predefined path | |
CN107415939B (zh) | 一种分布式驱动电动汽车转向稳定性控制方法 | |
Cai et al. | Implementation and development of a trajectory tracking control system for intelligent vehicle | |
CN103895704B (zh) | 基于后轮主动转向的变传动比控制方法 | |
CN107015477A (zh) | 基于状态反馈的车辆路径跟踪h∞控制方法 | |
Katsuyama et al. | A state-of-the-art review: toward a novel vehicle dynamics control concept taking the driveline of electric vehicles into account as promising control actuators | |
Chen et al. | Design of vehicle running states-fused estimation strategy using Kalman filters and tire force compensation method | |
Huang et al. | Lateral Stability Control of Four‐Wheel Independent Drive Electric Vehicles Based on Model Predictive Control | |
CN101867339A (zh) | 一种电子机械制动系统的电机控制方法 | |
Goodarzi et al. | Stability enhancement and fuel economy of the 4-wheel-drive hybrid electric vehicles by optimal tyre force distribution | |
Chen et al. | Implementation of MPC-based trajectory tracking considering different fidelity vehicle models | |
Shi et al. | DLMPCS‐based improved yaw stability control strategy for DDEV | |
Anderson et al. | Modelling minimum-time manoeuvering with global optimisation of local receding horizon control | |
Fu et al. | Overview of longitudinal and lateral control for intelligent vehicle path tracking | |
Hu et al. | MTTE-based motion stabilization control for in-wheel motor electric vehicles | |
Li et al. | Design and application of a tire slip energy control scheme based on optimal guaranteed cost theory for four‐in‐wheel drive electric vehicles | |
Yue-Lin et al. | A research on adaptive neural network control strategy of vehicle yaw stability | |
Zhao et al. | Integrated longitudinal and lateral control system design and case study on an electric vehicle |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20221028 Address after: 134300 No. 1601, Changbaishan street, Badaojiang District, Baishan City, Jilin Province (Development Zone) Patentee after: STATE GRID JILINSHENG ELECTRIC POWER SUPPLY COMPANY BAISHAN POWER SUPPLY Co. Address before: No.2055 Yan'an street, Chaoyang District, Changchun City, Jilin Province Patentee before: Changchun University of Technology |