基于非奇异终端滑模控制的汽车横向稳定性控制方法

技术领域

本发明涉及一种汽车稳定性的控制方法，特别涉及一种基于非奇异终端滑模控制的汽车横向稳定性控制方法，属于车辆控制技术领域。

背景技术

随着社会的飞速发展，电动汽车在汽车和交通运输行业中占据越来越重要的作用。其中横向稳定性控制是分布式驱动电动汽车安全行驶的关键技术之一。通常横向稳定性控制是指汽车在极限工况下，例如在冰雪路面行驶面对突发情况急转弯时可能会出现失稳的情况下回稳的控制。

目前，用户对车辆行驶稳定性和安全性能的要求不断提高，要求车辆在较为恶劣的行驶工况下依然保持良好的稳定性和可靠性。然而，由于车辆的失稳现象与多个参数相关，在以往的研究中未考虑的部分会导致控制性能受限，横摆角速度和质心侧偏角作为关键参数在不同的工况下参考值并不相同，因此从系统综合控制的角度出发，实时计算考虑多个关键因素的最优直接横摆力矩对车辆的横向动力学控制是极为重要的。在下层转矩分配中，多数研究采用的高精度算法大多需要复杂的迭代，然后在实车中受限于车载ECU的计算能力，复杂的高精度算法的使用会受限。

发明内容

发明目的：针对现有技术中极限工况下分布式驱动电动汽车的横向稳定性控制精度和安全性问题，本发明提出一种更高精度的控制模型和更可靠的基于非奇异终端滑模控制的汽车横向稳定性控制方法。

技术方案：一种基于非奇异终端滑模控制的汽车横向稳定性控制方法，包括以下步骤：

步骤1、获取影响车辆稳定性的关键因素，建立车辆动力学模型，通过动态仿真得到影响车辆稳定性的关键因素为横摆角速度和质心侧偏角；

步骤2、获取横摆角速度和质心侧偏角的阈值，建立车辆横向运动和横摆运动的车辆二自由度模型，通过动态仿真得到横摆角速度和质心侧偏角的阈值；

步骤3、获取最优附加横摆力矩，通过采集车辆运行过程中的横摆角速度和质心侧偏角实测值与步骤2的阈值比较，根据横摆角速度和质心侧偏角实测值与阈值的偏差值采用非奇异终端滑模控制求解得到横摆角速度和质心侧偏角的最优附加横摆力矩；

步骤4、获取最终附加横摆力矩，根据步骤3中横摆角速度和质心侧偏角的最优附加横摆力矩，通过赋予权重的方法合并获得最终附加横摆力矩；

步骤5、设计最优四轮转矩分配算法，以四轮的负荷率之和最低为目标函数，结合四轮权重系数、质心侧偏角速度与质心侧偏角相图的稳定性判别表达式，通过等式约束带入目标函数的方法对其求偏导，得到最优四轮转矩的分配算法；

步骤6、装车应用，根据步骤1-5中的方法设计控制系统刷写至分布式驱动电动汽车车载控制单元上，通过车上现有的传感器得到控制器所需要的参数判定车辆的稳定性状态，控制系统计算出实时所需横摆力矩，再计算出四个车轮转矩分配给各车轮。

本发明针对极限工况下分布式驱动电动汽车的横向稳定性控制精度和安全性问题，提供更高精度的控制模型和更可靠的基于非奇异终端滑模控制的汽车横向稳定性控制方法，提高在极限行驶工况下分布式驱动电动汽车的横向稳定性能。

优选项，为了获取影响车辆稳定性的关键因素，建立步骤1中所述的车辆动力学模型包括车辆的纵向运动、横向运动、横摆运动以及四个车轮旋转自由度的七自由度车辆模型和轮胎非线性力学模型，

所述七自由度车辆模型的数学表达式如下：

其中F

所述轮胎非线性力学模型数学表达式如下：

其中，F

根据上述方程式，组合条件下的轮胎纵向力F

在纵向加速度和横向加速度的影响下，每个车轮的垂直载荷的表达式如下：

其中，a

基于车轮动力学分析，各个轮胎侧滑角可描述为：

其中α

各车轮转动的动力学方程为：

式中，J为车轮转动惯量，ω

通过以上对车辆各状态参数的计算，可以得到各车轮的中心速度：

v

通过求解上面各车轮的中心速度，得到纵向滑移率的表达式：

λ

优选项，为了获取横摆角速度和质心侧偏角的阈值，所述步骤2中车辆横向运动和横摆运动的车辆二自由度模型的状态方程如下：

其中，k

横摆角速度阈值的中间值r

其中，μ为路面附着系数，L为车辆轴距，k

K为稳定性因数表达式为：

质心侧偏角阈值的取值范围用两条平行线表示的稳定性区域，其表达式如下：

其中，B

通过车辆横向运动和横摆运动的车辆二自由度模型的仿真获取不同车速下的横摆角速度阈值的取值范围为θ。

优选项，为了获取最优附加横摆力矩，所述步骤3中获取最优附加横摆力矩的方法如下：

(1)根据横摆角速度实测值与阈值的偏差值，设计非奇异终端滑模控制方法，输入值为横摆角偏差，输出值为基于横摆角速度的最优附加横摆力矩；

设计横摆角偏差和非奇异终端滑模面为：

其中，e

s

根据步骤1车辆动力学建模所述，车辆横摆角速度对时间的导数

其中，F

通过上式，滑模面对时间的导数

对于车辆的横向稳定性控制，在非奇异终端滑模控制方法中选择了指数趋近律，指数趋近律被表示为：

其中，K

基于横摆角速度得到的最优附加横摆力矩表达式为：

(2)根据质心侧偏角实测值与阈值的偏差值，设计非奇异终端滑模控制方法，输入值为质心侧偏角偏差，输出值为基于质心侧偏角的最优附加横摆力矩；

质心侧偏角偏差和非奇异终端滑模面设计为：

e

其中，e

在质心侧偏角非奇异终端滑模控制方法中也选择了指数趋近律，滑模面s

其中K

结合上式，基于质心侧偏角误差的最优附加横摆力矩定义为：

优选项，为了获取最终附加横摆力矩，所述步骤4中获取最终附加横摆力矩的方法如下：

基于横摆角速度和质心侧偏角的最优附加横摆力矩，通过赋予权重的方法合并获得最终附加横摆力矩，权重系数由粒子群算法求解得到，性能指标包括横摆角速度和质心侧偏角的横向稳定性判定方法和车辆的侧向加速度；

(1)基于最佳权重系数的最终附加横摆力矩表示为

ΔM＝PΔM

其中，P是两个计算值的权重系数，当P较大时，基于横摆角速度跟踪误差求解出的最优横摆力矩M

(2)建立以下指标函数用于确定分配系数P

J＝w

式中，W

式中，r为横摆角速度，r

J

优选项，为了设计最优四轮转矩分配算法，所述步骤5中设计最优四轮转矩分配算法如下：

以四轮的负荷率之和最低为目标函数，结合四轮权重系数、质心侧偏角速度与质心侧偏角相图的稳定性判别表达式，通过等式约束带入目标函数的方法对其求偏导，得到最优四轮转矩的分配算法；

(1)根据实际速度和目标速度之间的差异，使用比例积分(PI)控制器计算总扭矩，控制过程可以表示如下

其中，v

(2)定义各个轮胎的负荷率表示为：

由于电机只能控制纵向力，同时将F

以各个轮胎的负荷率并赋于不同权重之和为目标函数：

minJ＝minξ

其中，minJ为各个车轮的轮胎利用率权重总和，minξ

定义各个车轮的权重系数分别为：

其中，B

总的驱动力矩与横摆力矩和各个车轮力矩的关系如下：

两个前轮转矩可以被求解为：

把上述公式带入目标函数得到：

把目标函数J对T

令

有益效果：本发明通过完成车辆动力学的建模，在此基础上完成了整车横向动力学系统控制模型的构建，基于横摆角速度和质心侧偏角实际值与参考值的偏差，通过非奇异终端滑模控制算法求解得到两个横摆力矩，并利用粒子群算法对两者进行自适应分配得到最终横摆力矩；最终，建立了一种简洁可靠的最优转矩分配算法计算出实时的四轮最优转矩，提高在特殊行驶工况下智能汽车的横向稳定性能。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据提供的附图获得其他的附图；

图1为本发明横向稳定性控制流程图；

图2为本发明粒子群算法的流程图；

图3为本发明整车控制的流程图。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

在本发明的描述中，需要理解的是，术语“上”、“下”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明的限制。

在本发明中，除非另有明确的规定和限定，第一特征在第二特征之“上”或之“下”可以包括第一和第二特征直接接触，也可以包括第一和第二特征不是直接接触而是通过它们之间的另外的特征接触。而且，第一特征在第二特征“之上”、“上方”和“上面”包括第一特征在第二特征正上方和斜上方，或仅仅表示第一特征水平高度高于第二特征。第一特征在第二特征“之下”、“下方”和“下面”包括第一特征在第二特征正下方和斜下方，或仅仅表示第一特征水平高度小于第二特征。

如图所示，一种基于非奇异终端滑模控制的汽车横向稳定性控制方法，包括以下步骤：

步骤1、获取影响车辆稳定性的关键因素，建立车辆动力学模型，通过动态仿真得到影响车辆稳定性的关键因素为横摆角速度和质心侧偏角；

所述的车辆动力学模型包括车辆的纵向运动、横向运动、横摆运动以及四个车轮旋转自由度的七自由度车辆模型和轮胎非线性力学模型，

所述七自由度车辆模型的数学表达式如下：

其中F

考虑到当DYC系统开始工作时，轮胎的力学特性总是表现出明显的非线性行为，因此必须准确地反映轮胎纵向力和轮胎横向力之间的耦合关系，即应反映轮胎滑移率对轮胎侧向力的影响以及轮胎侧滑角对轮胎纵向力的影响。在这项工作中，轮胎非线性力学采用以下魔术公式建模：

其中，F

根据上述方程式，组合条件下的轮胎纵向力F

纵向和横向自由度被考虑在内，在纵向加速度和横向加速度的影响下，每个车轮的垂直载荷是不同，每个车轮的垂直载荷的表达式如下：

其中，a

基于车轮动力学分析，各个轮胎侧滑角可描述为：

其中α

各车轮转动的动力学方程为：

式中，J为车轮转动惯量，ω

通过以上对车辆各状态参数的计算，可以得到各车轮的中心速度：

v

轮胎的纵向滑移率对轮胎的纵向受力也有很大影响；通过求解上面各车轮的中心速度，得到纵向滑移率的表达式：

λ

步骤2、获取横摆角速度和质心侧偏角的阈值，建立车辆横向运动和横摆运动的车辆二自由度模型，通过动态仿真得到横摆角速度和质心侧偏角的阈值；

所述步骤2中车辆横向运动和横摆运动的车辆二自由度模型的状态方程如下：

其中，k

同时，从车辆主动安全性和稳定性角度出发，尽量保证车辆轨迹的稳定跟随，横摆角速度阈值的中间值r

其中，μ为路面附着系数，L为车辆轴距，k

K为稳定性因数表达式为：

基于横摆角速度，质心侧偏角实际值与理想值的差进行综合分析设计了失稳判定方法，

DYC系统的主要功能是确保车辆横向运动的相位轨迹处于稳定区域。在湿滑路面的稳定性分析中dβ-β相平面对稳定性分析的精度较好。以道路附着系数μ＝0.2为例，通过给出不同的初始值(β(0),r(0))得到了质心侧偏角的相轨迹。在相图中，存在一个可以用两条平行线表示的稳定性区域，其表达式如下：

其中，B

此外，如果实际横摆角速度与其参考值之间的误差过大，车辆轨迹将与计划轨迹有很大偏差，因此车辆横向稳定性也无法得到保证。为了提高所提出的DYC系统的稳定性控制性能，通过车辆横向运动和横摆运动的车辆二自由度模型的仿真获取不同车速下的横摆角速度阈值的取值范围为θ。

步骤3、获取最优附加横摆力矩，通过采集车辆运行过程中的横摆角速度和质心侧偏角实测值与步骤2的阈值比较，根据横摆角速度和质心侧偏角实测值与阈值的偏差值采用非奇异终端滑模控制求解得到横摆角速度和质心侧偏角的最优附加横摆力矩；

所述步骤3中获取最优附加横摆力矩的方法如下：

(1)根据横摆角速度实测值与阈值的偏差值，设计非奇异终端滑模控制方法，输入值为横摆角偏差，输出值为基于横摆角速度的最优附加横摆力矩；

设计横摆角偏差和非奇异终端滑模面为：

其中，e

s

根据步骤1车辆动力学建模所述，车辆横摆角速度对时间的导数

其中，F

通过上式，滑模面对时间的导数

对于车辆的横向稳定性控制，在非奇异终端滑模控制方法中选择了指数趋近律，这能更快的减小横摆角误差，指数趋近律被表示为：

其中，K

基于横摆角速度得到的最优附加横摆力矩表达式为：

(2)根据质心侧偏角实测值与阈值的偏差值，设计非奇异终端滑模控制方法，输入值为质心侧偏角偏差，输出值为基于质心侧偏角的最优附加横摆力矩；

质心侧偏角偏差和非奇异终端滑模面设计为：

e

其中，e

在质心侧偏角非奇异终端滑模控制方法中也选择了指数趋近律，这能更快的减小质心侧偏角误差，滑模面s

其中K

结合上式，基于质心侧偏角误差的最优附加横摆力矩定义为：

步骤4、获取最终附加横摆力矩，根据步骤3中横摆角速度和质心侧偏角的最优附加横摆力矩，通过赋予权重的方法合并获得最终附加横摆力矩；

所述步骤4中获取最终附加横摆力矩的方法如下：

基于横摆角速度和质心侧偏角的最优附加横摆力矩，通过赋予权重的方法合并获得最终附加横摆力矩，权重系数由粒子群算法求解得到，性能指标包括横摆角速度和质心侧偏角的横向稳定性判定方法和车辆的侧向加速度；

(1)基于最佳权重系数的最终附加横摆力矩表示为

ΔM＝PΔM

其中，P是两个计算值的权重系数，权重系数P在确定最终横摆力矩中起着重要作用，当P较大时，基于横摆角速度跟踪误差求解出的最优横摆力矩M

(2)建立以下指标函数用于确定分配系数P

J＝w

式中，W

式中，r为横摆角速度，r

J

步骤5、设计最优四轮转矩分配算法，以四轮的负荷率之和最低为目标函数，结合四轮权重系数、质心侧偏角速度与质心侧偏角相图的稳定性判别表达式，通过等式约束带入目标函数的方法对其求偏导，得到最优四轮转矩的分配算法；

所述步骤5中设计最优四轮转矩分配算法如下：

以四轮的负荷率之和最低为目标函数，结合四轮权重系数、质心侧偏角速度与质心侧偏角相图的稳定性判别表达式，通过等式约束带入目标函数的方法对其求偏导，得到最优四轮转矩的分配算法；

(1)车辆的纵向速度应实时跟踪该路径下速度的目标值，因此设计了速度跟踪控制器以确定作用在四个车轮上的总扭矩。根据实际速度和目标速度之间的差异，使用比例积分(PI)控制器计算总扭矩，控制过程可以表示如下

其中，v

(2)以各个轮胎利用率之和为目标函数，轮胎利用率可以用来表征车辆的稳定裕度。当J越大，稳定裕度越小，车辆发生失稳的可能性越大。相反，车辆则有更多的稳定裕度，以克服可能出现的失稳状况。定义各个轮胎的负荷率表示为：

由于电机只能控制纵向力，同时将F

以各个轮胎的负荷率并赋于不同权重之和为目标函数：

min J＝minξ

其中，min J为各个车轮的轮胎利用率权重总和，minξ

定义各个车轮的权重系数分别为：

其中，B

总的驱动力矩与横摆力矩和各个车轮力矩的关系如下：

两个前轮转矩可以被求解为：

把上述公式带入目标函数得到：

把目标函数J对T

令

步骤6、装车应用，根据步骤1-5中的方法设计控制系统刷写至分布式驱动电动汽车车载控制单元上，通过车上现有的传感器得到控制器所需要的参数判定车辆的稳定性状态，控制系统计算出实时所需横摆力矩，再计算出四个车轮转矩分配给各车轮。

本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其他实施例的不同之处，各个实施例之间相同相似部分互相参见即可。对于实施例公开的装置而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。

对所公开的实施例的上述说明，使本领域专业技术人员能够实现或使用本发明。对这些实施例的多种修改对本领域的专业技术人员来说将是显而易见的，本文中所定义的一般原理可以在不脱离本发明的精神或范围的情况下，在其它实施例中实现。因此，本发明将不会被限制于本文所示的这些实施例，而是要符合与本文所公开的原理和新颖特点相一致的最宽的范围。