CN112298193B - 快速实时的后轮主动转向预测控制方法 - Google Patents

快速实时的后轮主动转向预测控制方法 Download PDF

Info

Publication number
CN112298193B
CN112298193B CN202011114070.1A CN202011114070A CN112298193B CN 112298193 B CN112298193 B CN 112298193B CN 202011114070 A CN202011114070 A CN 202011114070A CN 112298193 B CN112298193 B CN 112298193B
Authority
CN
China
Prior art keywords
vehicle
value
model
rear wheel
point
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
Application number
CN202011114070.1A
Other languages
English (en)
Other versions
CN112298193A (zh
Inventor
王萍
刘行行
许娟
陈虹
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Jilin University
Original Assignee
Jilin University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Jilin University filed Critical Jilin University
Priority to CN202011114070.1A priority Critical patent/CN112298193B/zh
Publication of CN112298193A publication Critical patent/CN112298193A/zh
Application granted granted Critical
Publication of CN112298193B publication Critical patent/CN112298193B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT 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/00Estimation 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/10Estimation 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
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT 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/00Details 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
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60WCONJOINT 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/00Details 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/0001Details of the control system
    • B60W2050/0019Control system elements or transfer functions
    • B60W2050/0028Mathematical models, e.g. for simulation
    • B60W2050/0031Mathematical model of the vehicle
    • B60W2050/0034Multiple-track, 2D vehicle model, e.g. four-wheel model

Landscapes

  • Engineering & Computer Science (AREA)
  • Automation & Control Theory (AREA)
  • Transportation (AREA)
  • Mechanical Engineering (AREA)
  • Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Human Computer Interaction (AREA)
  • Steering Control In Accordance With Driving Conditions (AREA)

Abstract

一种快速实时的后轮主动转向预测控制方法,属于车辆控制技术领域。本发明的目的是使得车辆在低附着路面行驶工况下,能够在所设计的控制器的作用下跟踪期望的横摆角速度同时能够抑制质心侧偏角从而保证车辆稳定性的快速实时的后轮主动转向预测控制方法。本发明步骤是:后轮主动转向车辆的模型搭建与仿真工况的构建、模型预测控制器设计、优化问题快速求解算法设计。本发明进一步提高了车辆在低附着行驶工况下的稳定性,有效降低求解时间,提高算法实时性,降低了控制器的实现难度。

Description

快速实时的后轮主动转向预测控制方法
技术领域
本发明属于车辆控制技术领域。
背景技术
后轮主动转向技术一种新兴的技术形式,使用这种技术的车辆在行驶过程中后轮也可转动一定的角度因而具有更高的机动性和灵活性。例如,在低速行驶时可以通过后轮转向减小转弯半径,提高车辆灵活性,高速行驶时可以通过后轮转角提高稳定性,因此有必要对后轮主动转向车辆的控制问题进行研究。与传统车辆相比,后轮主动转向技术使得车辆的结构更加复杂,控制量也更多,后轮主动转向也进一步增加了车辆动力学的复杂程度,如何良好地发挥后轮主动转向的作用改善车辆在高速、低附着路面等行驶工况下的稳定性能成为了目前的研究热点。由于车辆动力学具有高度非线性的特点,且车辆是一个典型的快变系统,因此在车辆的稳定性控制中要求控制器有较高的控制精度来满足控制准确性的要求,同时控制器还要有较高的实时性来满足车辆对控制器的快速响应的要求,除此之外,在车辆行驶过程中还要考虑相应的安全性约束以及在未来一段时间内车辆状态的变化规律,进行及时有效的控制。由于预测控制能够在考虑系统的约束和非线性特性的条件下,对系统未来一段时间内的状态进行预测,并处理多目标优化问题,因此在同时满足上述要求的条件下,预测控制是一个较好地解决方案。模型预测控制基于系统的模型对系统的未来状态进行预测,同时考虑所设计的目标函数以及相应的状态约束和控制约束,将控制需求转化为数学优化问题,并通过在线求解优化问题来得到所需的控制量。一方面由于车辆是一个典型的快变系统,其失稳时间往往都是毫秒级的,一旦车辆失去稳定性则无法再通过控制作用使其回到稳定状态,因此车辆对车载控制器的实时性与快速响应能力有较高的要求;另一方面求解优化问题需要的计算量较大,同时车载控制器的计算能力却相当有限,在求解优化问题时需要较长的求解时间。因此,当前的面向后轮主动转向车辆的模型预测控制存在着以下问题:
1.相比于传统车辆,后轮主动转向车辆的控制机构更多,复杂程度也更高。因此,当行驶在在高速和低附着路面等复杂工况下时,后轮主动转向车辆具有更高的灵活性和稳定性。而在此类复杂工况下,车辆往往都进入了高度非线性区域,此时所需要的模型精度也更高,因此需要采用非线性模型来提高模型精度,这会造成计算负担大幅增加。
2.随着计算负担的增加,若采用传统的优化方法(例如已被广泛应用的内点法-IPOPT),可能无法满足车载实时性的要求,从而导致控制器性能不能满足要求。
3.虽然车载控制器计算能力不足的问题可以通过安装更高性能的硬件设备来解决,但此举势必会造成车辆成本及能耗的大幅增加。若能够改进求解算法,使得能够在现有硬件设备的基础上提高求解速度,满足车辆对控制器的实时性的要求,便可节省大量成本同时还不增加整车的能耗。
发明内容
本发明的目的是使得车辆在低附着路面行驶工况下,能够在所设计的控制器的作用下跟踪期望的横摆角速度同时能够抑制质心侧偏角从而保证车辆稳定性的快速实时的后轮主动转向预测控制方法。
本发明步骤是:
步骤一、后轮主动转向车辆的模型搭建与仿真工况的构建:在CarSim软件中选择车辆模型,并将车辆的运动状态参数读取到Simulink中,基于所选的车辆模型构建低附着路面下行驶的仿真工况,模拟实际车辆的横摆运动和侧向运动特征;
步骤二、模型预测控制器设计:
1)建立考虑后轮主动转向的车辆动力学模型;
2)建立能够产生期望状态参考值的参考模型;
3)根据控制器设计中跟踪期望横摆角速度、抑制质心侧偏角以及限制控制作动能量的需求,并结合相应的约束和车辆动力学模型构建目标函数;
步骤三、优化问题快速求解算法设计
1)引入松弛函数,将初始的带约束的非线性优化问题转化为无约束非线性优化问题;
2)根据目标函数和车辆动力学模型构建哈密顿函数;
3)根据PMP原理推导出沿最优轨线的递推关系式,以及最优终端必要性条件;
4)根据Nelder-Mead单纯形算法原理设计最优协态变量初值搜索算法;
上述步骤详细过程如下:
1)车辆运动模型
根据车辆动力学的理论,简化后的车辆二自由度模型可由如下方程描述:
Figure GDA0003153790570000021
其中,
Figure GDA0003153790570000022
Figure GDA0003153790570000023
分别表示车辆的质心侧偏角的导数和车辆的横摆角速度的导数,V表示车辆的纵向速度,Fyf和Fyr则分别表示前后轮胎的轮胎侧向力,Lf和Lr分别表示前后轴到车辆质心的距离,m是车辆的质量,Iz是车辆绕质心旋转的转动惯量,ΔMz为附加的横摆力矩;
2)前轮轮胎模型建立
当前轮侧偏角αf很小时,有tan(αf)≈αf,之后该前轮模型可近似为:
Figure GDA0003153790570000024
其中μ为路面附着系数,Fz为垂直载荷,Cf为前轮侧偏刚度;
式中的前轮侧偏角αf可由下式进行计算:
Figure GDA0003153790570000025
其中δf车辆的前轮转角;
3)后轮轮胎模型建立
考虑路面的附着系数的影响,该模型可描述为:
Figure GDA0003153790570000026
其中的μ为路面附着系数,Fz为垂直载荷,Cr为后轮侧偏刚度;
式中的后轮侧偏角αr可由下式进行计算:
Figure GDA0003153790570000031
其中的δr为后轮转角;
4)车辆参考模型建立
前轮转角δf到车辆的期望横摆角速度γref以及期望质心侧偏角βref的传递函数为:
Figure GDA0003153790570000032
定义车辆的稳定因子为
Figure GDA0003153790570000033
其中L=Lf+Lr为车辆的前后轴距;
系统的振荡频率为
Figure GDA0003153790570000034
系统的阻尼系数为
Figure GDA0003153790570000035
横摆角速度稳态增益为
Figure GDA0003153790570000036
质心侧偏角稳态增益为
Figure GDA0003153790570000037
其中的横摆角速度微分系数定义为
Figure GDA0003153790570000038
质心侧偏角微分系数定义为
Figure GDA0003153790570000039
定义横摆角速度的上限值为
Figure GDA00031537905700000310
横摆角速度的参考值γref应被约束在|γref|≤γup,定义质心侧偏角的上限值为
Figure GDA00031537905700000311
横摆角速度的参考值βref应被约束在|βref|≤βup
5)控制器设计
定义系统的状态向量为x=[x1,x2]T=[β,γ]T,控制量定义为u=[u1,u2]=[δr,ΔMz],定义ΔMmax和δmax分别是附加横摆力矩和后轮主动转角的最大约束,将(1)式中的状态空间方程进行欧拉离散化得在每个采样时刻kTs离散后的状态空间方程为:
Figure GDA0003153790570000041
其中的前轮侧向力Fyf由式(2)中的非线性轮胎模型计算得到,后轮侧向力Fyr则由式(4)中的线性轮胎模型计算得到,βup和γup分别为质心侧偏角和横摆角速度的上限值,ΔMmax为附加横摆力矩的最大值;
在每个时刻k+1≤ki≤k+N+1定义代价函数为:
Figure GDA0003153790570000042
其中L1(ki)、L'2(ki)和L3(ki)分别用于跟踪横摆角速度参考值、跟踪质心侧偏角参考值和抑制作动能量;
状态约束为|x1(ki)|≤βup,控制量限幅为|u1(ki-1)|≤δmax以及|u2(ki-1)|≤ΔMmax
得到非线性模型预测控制的目标函数为:
Figure GDA0003153790570000043
满足的约束为|u1(ki-1)|≤δmax,|u2(ki-1)|≤ΔMmax以及|x1(ki)|≤βup,其中Γβ和Γu=[ΓrM]分别为质心侧偏角和控制量的权重系数,N为预测时域;
6)基于PMP原理的求解
由PMP原理可知,对于此类的问题得到:
Figure GDA0003153790570000044
若存在控制量u*(t)和以及终端时间
Figure GDA0003153790570000045
使得性能指标中的目标函数在该控制量的作用下取到最小的最优解,且在该最优解的作用下最优状态轨线x*(t)所形成的最优轨迹也存在,则可导出如下必要条件:
(1)x(t)和λ(t)满足正则方程:
Figure GDA0003153790570000051
其中λ=[λ12]T为拉格朗日乘子的向量函数,L(x,u)为目标函数中的各项,
H(x,λ,u)=L(x,u)+λT(t)f(x,u)为哈密顿函数;
(2)状态轨线x(t)和协态变量λ(t)满足边界条件:
Figure GDA0003153790570000052
(3)哈密顿函数H(x,λ,u)在最优控制量u*(t)的作用下取到绝对极小值:
Figure GDA0003153790570000053
(4)同时哈密顿函数在最优状态轨线的末端满足:
H[x*(t),λ(t),u*(t)]=H[x*(tf),λ(tf),u*(tf)]=constant (14)
由上述的状态约束|x1(ki)|≤1的基础上,引入终端松弛因子对状态约束进行处理:
Figure GDA0003153790570000054
其中的κ代表了函数的松弛程度,而ν则为一个较大的数,
最终得到目标函数的子项L2(ki)为:
L2(ki)=L′2(ki)+ζ(ki) (16)
此时式(10)中的目标函数则转化为:
Figure GDA0003153790570000055
该目标函数可满足状态约束|x1(ki)|≤1;
定义在时刻k+1≤ki≤k+N+1的哈密顿函数为:
Figure GDA0003153790570000056
其中F1(x(ki))和F2(x(ki))分别定义为:
Figure GDA0003153790570000061
在上述方程中,λ1(k)和λ2(k)分别表示拉格朗日乘子,根据庞特里亚金最小值原理可得到最优必要性条件为:
Figure GDA0003153790570000062
终端条件为:
Figure GDA0003153790570000063
在各个时刻都有最优控制律u*(ki)使得哈密顿函数最小化:
Figure GDA0003153790570000064
将哈密顿函数重新整理为关于二维控制量u=[u1,u2]T的二元二次函数:
Figure GDA0003153790570000065
其中g[x(ki),δf(k)]为与当前控制量无关的余项,其定义如下:
Figure GDA0003153790570000066
哈密顿函分别对两个控制量求偏导可得:
Figure GDA0003153790570000067
Figure GDA0003153790570000071
分别令上述两个偏导等于零便可求得哈密顿函数的驻点为
Figure GDA0003153790570000072
以及
Figure GDA0003153790570000073
对哈密顿函数相对于控制量求二阶偏导,其二阶偏导如下所示:
Figure GDA0003153790570000074
Figure GDA0003153790570000075
Figure GDA0003153790570000076
Figure GDA0003153790570000077
得AC-B2>0,且A>0,C>0可知哈密顿函数在驻点处取得极小值;
根据在驻点处的结果,便可给出使哈密顿函数取得极小值的最优控制律:
Figure GDA0003153790570000078
Figure GDA0003153790570000079
其中的第一个控制量u1为后轮主动转角,该控制量直接施加给后轮转向执行机构,第二个控制量u2为附加横摆力矩,该控制量首先需要通过如下转换:
Figure GDA00031537905700000710
其中ΔTfl(k),ΔTrl(k),ΔTfr(k),ΔTrr(k)分别表示左前轮、左后轮、右前轮以及右后轮的附加转矩,Re表示轮胎的滚动半径,d表示车辆车身的宽度,之后该控制量便转换为四个轮胎的附加转矩施加给车辆;
最优控制问题此时转化为了两点边值问题:
Figure GDA0003153790570000081
该两点边值问题通过Nelder-Mead算法进行求解;
7)基于Nelder-Mead的最优初值搜索
(1)检查给定初始测试点的收敛情况
三个初始测试点p1,p2,p3,每个测试点都是二维向量;各个测试点的终端收敛值记作
Figure GDA0003153790570000082
之后按照终端收敛值从小到大的顺序进行排序,现假设排序情况为:
Figure GDA0003153790570000083
(2)计算质心点
Figure GDA0003153790570000084
(3)计算反射点
pr=p0+α(p0-p3) (36)
其中α=0.5为反射系数,之后按照图三所示的迭代过程计算反射点pr的终端收敛值
Figure GDA0003153790570000085
之后进行如下判断;
(3.1)反射
如果反射点的收敛值
Figure GDA0003153790570000086
小于第二个测试点的收敛值
Figure GDA0003153790570000087
且大于第一个测试点的收敛值
Figure GDA0003153790570000088
此时将pr的值赋予p3并回到第(1)步重新计算排序;
(3.2)扩张
如果反射点的收敛值
Figure GDA0003153790570000089
小于第一个测试点的收敛值
Figure GDA00031537905700000810
即表示找到了一个更好的收敛方向,此时需要继续沿着该方向搜索,之后计算扩展点:
pe=p0+η(pr-p0) (37)
其中η=2为扩展系数,之后计算扩展点的收敛值
Figure GDA00031537905700000811
如果
Figure GDA00031537905700000812
则将pe的值赋予p3并回到第(1)步重新计算排序,否则将pr的值赋予p3并回到第(1)步重新计算排序;(3.3)紧缩
如果反射点收敛值满足
Figure GDA00031537905700000813
则计算紧缩点:
pc=p0+ρ(p3-p0) (38)
其中ρ=0.5为紧缩系数,如果紧缩点的收敛值满足
Figure GDA0003153790570000091
则将pc的值赋予p3并回到第(1)步重新计算排序;
(3.4)收缩
如果以上条件均不满足,则证明此时的测试点已经很接近最优点了,此时需要沿着当前最优点的方向进行收缩,即将第二和第三个测试点按照如下方式进行替换:
pi=p1+σ(pi-p1),i=2,3 (39)
其中σ=0.5为收缩系数,替换完成后回到第(1)步重新计算排序。
本发明的有益效果是:
1、本发明针对低附着路面行驶工况,面向后轮主动转向车辆设计了非线性模型预测控器,与传统车辆相比,在车辆能够跟踪期望的横摆角速度的条件下,同时能够大幅抑制质心侧偏角,进一步提高了车辆在低附着行驶工况下的稳定性;
2、本发明根据PMP原理结合后轮主动转向车辆的动力学模型,设计了面向后轮主动转向车辆的模型预测控制中的非线性优化问题的间接求解算法,与传统的IPOPT方法相比,本发明所提出的求解算法可以有效降低求解时间,提高算法实时性;
3、本发明根据PMP原理将原最优控制问题进行转化,并使用Nelder-Mead单纯形算法进行搜索求解,该方法在搜索过程中无需数学求导,也不需要复杂的数值优化运算,降低了控制器的实现难度。
附图说明
图1是发明中所设计的控制器框图;
图2是车辆模型示意图;
图3是求解过程中的迭代关系示意图;
图4是本发明所设计的快速求解算法的控制效果;
图5是使用传统IPOPT求解算法的控制效果;
图6是在本发明中所设计的仿真工况中不施加控制的效果图;
图7是在仿真过程中车辆的质心侧偏角对比图;
图8是在仿真过程中不同求解算法的平均求解时间;
图9是在仿真过程中在不同预测时域下基于传统IPOPT求解算法的求解时间;
图10是在仿真过程中在不同预测时域下基于本发明所设计的快速求解算法的求解时间。
具体实施方式
本发明涉及一种快速实时的后轮主动转向预测控制器设计方法,具体涉及一种模型预测控制框架下的后轮主动转向控制方法,同时设计了面向车载应用的控制器快速实时求解方法。更具体地说,由于在车辆的稳定性控制中,不仅需要同时考虑车辆动力学的非线性特性和车辆的安全性约束,还需要考虑车辆在未来一段时间内的运动状态从而防止车辆失稳。基于以上考虑,本发明基于后轮主动转向技术设计了后轮主动转向预测控制器。但由于模型预测控制需要在线求解优化问题,计算负担较大,不能够满足车载应用的实时性要求。因此本发明同时设计了一种预测控制快速实时求解的方法,使得所设计的模型预测控制器能够加快计算速度,满足车载实时性要求。
由于在复杂行驶工况下,后轮主动转向车辆具有更好的灵活性和更高的稳定性,因此本发明在考虑车辆动力学非线性特性和安全性约束的条件下,针对在低附着路面下的后轮主动转向车辆的稳定性控制问题设计了一种模型预测控制器,使得车辆在低附着路面行驶工况下,能够在所设计的控制器的作用下跟踪期望的横摆角速度同时能够抑制质心侧偏角从而保证车辆的稳定性。但在通常情况下,出于控制成本和整车能耗等方面的考虑,车载控制器的计算能力往往较为有限,在求解较为复杂的优化问题时不能够满足车辆的实时性要求。因此在求解本发明中所设计的模型预测控制中的优化问题时,若采用传统的优化算法,例如序列二次规划(SQP)和内点法(IPOPT),计算时间会随着预测时域呈指数级增加,不能满足车载实时性要求。为了解决这一问题,本发明同时设计了一种预测控制的快速实时求解方法。具体为基于庞特里亚金极值原理(PMP)将模型预测控制中的非线性优化问题转化为显式的迭代计算问题进行求解,同时为了寻找满足最优性终端必要条件的最优协态变量初值,本发明采用了Nelder-Mead单纯形算法来搜索最优协态变量的初值。通过PMP和Nelder-Mead算法的结合,可以大大加快模型预测中优化问题的求解速度。为了证明本发明所设计的控制器的有效性以及优化算法的快速性,需要设计相应的仿真实验来进行验证与对比。本发明设计了一组车辆在低附着路面上的仿真实验来验证控制器的有效性。而在求解算法的对比中,采用了应用较为广泛的内点法(IPOPT)来和本发明所设计的快速求解算法进行对比。
本发明包括以下步骤:
步骤一、后轮主动转向车辆的模型搭建与仿真工况的构建:在CarSim软件中选择车辆模型,并将车辆的运动状态参数读取到Simulink中,基于所选的车辆模型构建低附着路面下行驶的仿真工况,模拟实际车辆的横摆运动和侧向运动特征;
步骤二、模型预测控制器设计:
1)建立考虑后轮主动转向的车辆动力学模型;
2)建立能够产生期望状态参考值的参考模型;
3)根据控制器设计中跟踪期望横摆角速度、抑制质心侧偏角以及限制控制作动能量的需求,并结合相应的约束和车辆动力学模型构建目标函数;
步骤三、优化问题快速求解算法设计
1)引入松弛函数,将初始的带约束的非线性优化问题转化为无约束非线性优化问题;
2)根据目标函数和车辆动力学模型构建哈密顿函数;
3)根据PMP原理推导出沿最优轨线的递推关系式,以及最优终端必要性条件;
4)根据Nelder-Mead单纯形算法原理设计最优协态变量初值搜索算法。
为详细说明本发明的技术内容、构造特点、实现目的等下面结合附图对本发明进行全面解释:
基于本发明所述的控制器框图如图1所示,其中驾驶员为CarSim中自带的驾驶员模型,在仿真验证中负责操纵车辆完成设定的工况。驾驶员在驾驶车辆时会打出相应的方向盘转角,参考模型根据驾驶的方向盘转角计算出期望的横摆角速度和质心侧偏角,并将期望值送入控制器。控制器根据期望值和车辆当前的状态计算出相应的控制量,在本发明中的控制量为附加横摆力矩和后轮主动转角。其中附加横摆力矩会送入转矩分配模块,将附加横摆力矩转化为四个轮胎的附加转矩并施加在车辆上。由控制器计算所得的后轮转角则直接施加给车辆的后轮转向机构。在仿真实验中,车辆模型和仿真工况均是在CarSim中构建的,控制器则是在Simulink中构建的。
本发明所述的控制器是通过软件系统的联合仿真实现并进行验证的1、软件选择
该快速求解算法及其相应的控制器和由控制器控制的被控对象的仿真模型分别通过软件Matlab/Simulink和高保真车辆动力学仿真软件CarSim进行搭建,软件版本分别为MatlabR2019b和CarSim2016.1,求解器选择为ODE1。仿真步长为0.001s。其中CarSim用于提供高保真的车辆动力学模型以及相应的仿真工况;MATLAB/Simulink软件则是用于控制器的搭建,即通过Simulink编程来完成该方法中控制器的运算。
2、联合仿真设置
要实现两者的联合仿真,首先需要在Matlab的路径设置中添加CarSim的路径;其次在CarSim界面中添加输出接口模块;然后将CarSim中的模型信息经过系统编译之后以CarSimS-function的形式保留在Simulink中,最后再进行Simulink中CarSim模块的参数设置。
为了对本发明所述的面向后轮主动转向车辆的模型预测控制快速求解算法进行验证,选择了低附着路面下的车辆稳定性控制的工况来进行验证。首先推导出后轮主动转向车辆的数学模型;其次,在CarSim中选择合适的车辆模型并获取相应参数;之后在低附着路面下构建仿真工况;然后设计面向后轮主动转向车辆的模型预测控制器并使用本发明所述的快速求解算法进行控制量的求解;最后,在联合仿真实验中对本发明所述的方法进行验证,同时与传统的IPOPT求解算法进行对比,以说明本发明的有益效果。
本发明具体包括以下步骤:
步骤一、后轮主动转向车辆模型搭建:车辆模型模拟真实的被控对象,主要作用是能够精确的模拟实际的后轮主动转向车辆的横摆运动和侧向运动特征。
在本发明中,由于使用了联合仿真,因此在CarSim中,主要用到的是车辆模型选择以及仿真工况的构建。
首先选择典型的乘用车模型,之后对模型的相关参数进行修改并获取,将车辆模型参数添加到Simulink仿真模型中。车辆的主要模型参数有车辆质量、前后轴距轮胎侧偏刚度等。在选择相应的车辆模型和参数之后,需要构建相应的仿真工况,仿真工况中可以选择车辆的行驶路线,行驶环境以及驾驶员模型等。由于在本发明中只是对通过附加横摆力矩来对车辆的稳定性进行控制,因此选用CarSim中自带的驾驶员模型,并将车辆的运动状态参数读取到Simulink中,基于所选的车辆模型构建低附着路面下的仿真工况。
步骤二、基于快速求解算法的控制器设计:基于模型预测控制原理对所需要的控制问题进行描述,描述为一个非线性优化问题,并基于本发明所述的快速求解算法进行控制器的设计。
本发明的被控对象是行驶在低附着路面上的后轮主动转向车辆,因此控制目标就是根据当前车辆的状态信息以及当前的参考值计算出所需要的后轮转角以及附加横摆力矩,提高车辆的稳定性。主要设计过程描述如下。首先,针对后轮主动转向车辆的横摆运动和侧向运动进行建模,建立合适的数学模型。
1)车辆运动模型
在本发明中采用了已被广泛应用的二自由度模型来对车辆的运动进行描述,在该模型中只考虑了车辆的侧向运动和横摆运动。如图2所示,车辆被描述为一个只有单个前轮和单个后轮的车辆模型。在本发明中,驾驶员通过方向盘操纵前轮,而后轮的转向则由相应的后轮转向机构来操纵。此时车辆模型便可以简化为车辆二自由度模型。同时根据车辆动力学的理论,简化后的车辆二自由度模型可由如下方程描述:
Figure GDA0003153790570000111
其中,
Figure GDA0003153790570000112
Figure GDA0003153790570000113
分别表示车辆的质心侧偏角的导数和车辆的横摆角速度的导数,V表示车辆的纵向速度,Fyf和Fyr则分别表示前后轮胎的轮胎侧向力,Lf和Lr分别表示前后轴到车辆质心的距离,m是车辆的质量,Iz是车辆绕质心旋转的转动惯量,ΔMz为附加的横摆力矩。
2)前轮轮胎模型建立
在本发明中,为了提高模型精度,前轴轮胎的侧向力是由一个非线性模型来描述的,采用了Fiala轮胎模型进行描述。在该模型中,使用了前轴轮胎的侧偏角作为内部变量。当前轮侧偏角αf很小时,有tan(αf)≈αf,之后该前轮模型可近似为:
Figure GDA0003153790570000114
其中μ为路面附着系数,Fz为垂直载荷,Cf为前轮侧偏刚度。
式中的前轮侧偏角αf可由下式进行计算:
Figure GDA0003153790570000121
其中δf车辆的前轮转角。
3)后轮轮胎模型建立
在本发明中,为了保证优化问题的哈密顿函数能够在驻点处取到极小值,同时也为了在保证精度的情况下降低模型的复杂度,后轮的轮胎模型选择为一个线性的轮胎模型。与前轮轮胎模型类似,在后轮轮胎模型中同样选取轮胎的侧偏角作为内部变量。同时考虑路面的附着系数的影响,该模型可描述为:
Figure GDA0003153790570000122
其中的μ为路面附着系数,Fz为垂直载荷,Cr为后轮侧偏刚度。
式中的后轮侧偏角αr可由下式进行计算:
Figure GDA0003153790570000123
其中的δr为后轮转角,在本发明中后轮转角有控制器计算得到并直接施加给后轮转向的执行机构。
4)车辆参考模型建立
在本发明中,车辆在行驶过程中需要计算出当前的期望状态,其中包括期望的横摆角速度和期望的质心侧偏角。在本发明中,采用了已被广泛应用的二阶参考模型来计算车辆的期望状态。在参考模型中,根据当前的前轮转角以及相应的传递函数来计算期望状态。
其中前轮转角δf到车辆的期望横摆角速度γref以及期望质心侧偏角βref的传递函数为:
Figure GDA0003153790570000124
此时我们定义车辆的稳定因子为
Figure GDA0003153790570000125
其中L=Lf+Lr为车辆的前后轴距。系统的振荡频率为
Figure GDA0003153790570000126
系统的阻尼系数为
Figure GDA0003153790570000131
横摆角速度稳态增益为
Figure GDA0003153790570000132
质心侧偏角稳态增益为
Figure GDA0003153790570000133
其中的横摆角速度微分系数定义为
Figure GDA0003153790570000134
质心侧偏角微分系数定义为
Figure GDA0003153790570000135
由于车辆的二阶参考模型没有考虑车辆所行驶的道路环境,因此计算所得到的期望值还需要进一步处理。当路面的附着系数较低时,轮胎所能产生的轮胎力的最大值不足以支撑所需要的较大的横摆角速度。此时我们需要对横摆角速度的参考值进行适当的限制,来适应路面的摩擦系数。为了达到满意的性能,首先定义横摆角速度的上限值为
Figure GDA0003153790570000136
横摆角速度的参考值γref应被约束在|γref|≤γup。类似地,我们定义质心侧偏角的上限值为
Figure GDA0003153790570000137
横摆角速度的参考值βref应被约束在|βref|≤βup
5)控制器设计
针对后轮主动转向车辆,为了提高车辆的稳定性,可以通过控制作用使车辆能够跟踪期望的横摆角速度和质心侧偏角。因此在本发明所设计的控制器中有两个控制量,分别是作用于车辆附加横摆力矩和作用于后轮转向机构的后轮主动转角。控制器的框图如图1所示,其中计算所得的附加横摆力矩通过转矩分配模块分别施加到四个轮胎上,后轮转角则直接施加到后轮转向机构上。
首先定义系统的状态向量为x=[x1,x2]T=[β,γ]T。控制量定义为u=[u1,u2]=[δr,ΔMz]。同时定义ΔMmax和δmax分别是附加横摆力矩和后轮主动转角的最大约束。由于硬件设备需要进行离散计算,首先需要将状态空间方程进行离散化,将(1)式中的状态空间方程进行欧拉离散化可得在每个采样时刻kTs离散后的状态空间方程为:
Figure GDA0003153790570000138
其中的前轮侧向力Fyf由式(2)中的非线性轮胎模型计算得到,后轮侧向力Fyr则由式(4)中的线性轮胎模型计算得到,βup和γup分别为质心侧偏角和横摆角速度的上限值,ΔMmax为附加横摆力矩的最大值。
在本发明中,目标函数由三项组成,分别为跟踪期望横摆角速度,跟踪质心侧偏角参考值以及考虑作动能量的最小化。因此在每个时刻k+1≤ki≤k+N+1定义代价函数为:
Figure GDA0003153790570000141
其中L1(ki)、L'2(ki)和L3(ki)分别用于跟踪横摆角速度参考值、跟踪质心侧偏角参考值和抑制作动能量。处于车辆的行驶安全性考虑,需要对车辆的状态进行约束,可得状态约束为|x1(ki)|≤βup,控制量限幅为|u1(ki-1)|≤δmax以及|u2(ki-1)|≤ΔMmax
此时便可以得到非线性模型预测控制的目标函数为:
Figure GDA0003153790570000142
满足的约束为|u1(ki-1)|≤δmax,|u2(ki-1)|≤ΔMmax以及|x1(ki)|≤βup,其中Γβ和Γu=[ΓrM]分别为质心侧偏角和控制量的权重系数,N为预测时域。
6)基于PMP原理的求解算法
通过上述过程可知,本发明中所设计的目标函数实质上是一种典型的积分型性能指标目标函数,且控制受约束的最优控制问题。由PMP原理可知,对于此类的问题可以得到:
Figure GDA0003153790570000143
若存在控制量u*(t)和以及终端时间
Figure GDA0003153790570000144
使得性能指标中的目标函数在该控制量的作用下取到最小的最优解。
且在该最优解的作用下最优状态轨线x*(t)所形成的最优轨迹也存在,则可导出如下必要条件:
(1)x(t)和λ(t)满足正则方程:
Figure GDA0003153790570000145
其中λ=[λ12]T为拉格朗日乘子的向量函数,L(x,u)为目标函数中的各项,H(x,λ,u)=L(x,u)+λT(t)f(x,u)为哈密顿函数。
(2)状态轨线x(t)和协态变量λ(t)(也可称为拉格朗日乘子)满足边界条件:
Figure GDA0003153790570000151
(3)哈密顿函数H(x,λ,u)在最优控制量u*(t)的作用下取到绝对极小值:
Figure GDA0003153790570000152
(4)同时哈密顿函数在最优状态轨线的末端满足:
H[x*(t),λ(t),u*(t)]=H[x*(tf),λ(tf),u*(tf)]=constant (53)
由上述基于PMP原理导出的必要性条件可知,基于PMP原理的最优控制是无法处理系统状态约束的,而出于安全性的考虑,系统状态约束不得不考虑,因此需要对原状态约束进行转化处理。由上述的状态约束|x1(ki)|≤1的基础上,引入终端松弛因子对状态约束进行处理:
Figure GDA0003153790570000153
其中的κ代表了函数的松弛程度,而ν则为一个较大的数,用来保证最优的状态轨迹是在约束范围之内的。通过该松弛因子的定义可知,在求解过程中若状态接近约束值,则松弛因子的值将会急剧增大,而当状态离边界约束较远时,松弛因子的值约等于零,此时对目标函数整体的影响可忽略不计。
最终可得到目标函数的子项L2(ki)为:
L2(ki)=L′2(ki)+ζ(ki) (55)
此时式(10)中的目标函数则转化为:
Figure GDA0003153790570000154
该目标函数可满足状态约束|x1(ki)|≤1。
根据PMP原理可知,与系统模型的离散化类似,需要将上述的最优控制必要条件离散化处理。再根据离散化的系统状态空间模型,现定义在时刻k+1≤ki≤k+N+1的哈密顿函数为:
Figure GDA0003153790570000155
其中F1(x(ki))和F2(x(ki))分别定义为:
Figure GDA0003153790570000161
在上述方程中,λ1(k)和λ2(k)分别表示拉格朗日乘子,根据庞特里亚金最小值原理可得到最优必要性条件为:
Figure GDA0003153790570000162
终端条件为:
Figure GDA0003153790570000163
在各个时刻都有最优控制律u*(ki)使得哈密顿函数最小化:
Figure GDA0003153790570000164
基于必要性条件,我们便可以给出从初始状态到终端状态的映射关系,如图3所示。
基于上述讨论,在初始条件λ(ki)和x(ki)已知的条件下,便可基于PMP原理,给出优化控制问题的最优解。为了简便起见,现将哈密顿函数重新整理为关于二维控制量u=[u1,u2]T的二元二次函数:
Figure GDA0003153790570000165
其中g[x(ki),δf(k)]为与当前控制量无关的余项,其定义如下:
Figure GDA0003153790570000166
哈密顿函分别对两个控制量求偏导可得:
Figure GDA0003153790570000171
Figure GDA0003153790570000172
分别令上述两个偏导等于零便可求得哈密顿函数的驻点为
Figure GDA0003153790570000173
以及
Figure GDA0003153790570000174
在求得哈密顿函数的驻点之后,还需要进一步判断在驻点处哈密顿函数是否可以取到极小值。
由上述哈密顿函数的形式可知,该哈密顿函数为一个二元二次函数,由二元函数极值判定条件可知,为了判断二元二次函数在驻点处是否取得极值,以及取到的极值为极大值或极小值,需要对哈密顿函数相对于控制量求二阶偏导,其二阶偏导如下所示:
Figure GDA0003153790570000175
Figure GDA0003153790570000176
Figure GDA0003153790570000177
Figure GDA0003153790570000178
可得AC-B2>0,且A>0,C>0可知哈密顿函数在驻点处取得极小值。根据在驻点处的结果,便可给出使哈密顿函数取得极小值的最优控制律:
Figure GDA0003153790570000179
Figure GDA00031537905700001710
其中的第一个控制量u1为后轮主动转角,该控制量直接施加给后轮转向执行机构,第二个控制量u2为附加横摆力矩。
该控制量首先需要通过如下转换:
Figure GDA0003153790570000181
其中ΔTfl(k),ΔTrl(k),ΔTfr(k),ΔTrr(k)分别表示左前轮、左后轮、右前轮以及右后轮的附加转矩,Re表示轮胎的滚动半径,d表示车辆车身的宽度,之后该控制量便转换为四个轮胎的附加转矩施加给车辆。
通过上述过程可知,原最优控制问题在经过PMP原理转化后,变成了求解最优拉格朗日乘子初值λ*(k)的问题,若可以找到最优初值,使得在该初值条件下哈密顿函数取最小值且终端协态变量可以收敛到零,便可得到最优控制序列。换句话说,原最优控制问题此时转化为了两点边值问题:
Figure GDA0003153790570000182
该两点边值问题在本发明中通过Nelder-Mead算法进行求解。
7)基于Nelder-Mead的最优初值搜索算法
为了找到满足最优必要性条件的最优初值λ*(k),本发明采用了Nelder-Mead算法进行搜索最优初值。在该搜索算法中,可以根据给定的测试点进行搜索,经过事先定义的基本搜索环节后便可找到满足最优必要性条件的最优初值。
在本发明中,搜索的详细过程如下:
(1)检查给定初始测试点的收敛情况
由上述过程可知,在本发明中后轮主动转向车辆是一个二阶系统,因此需要找到一个二阶的最优初值使其满足最优必要性条件。因此需要给定三个初始测试点p1,p2,p3,每个测试点都是二维向量。按照图3所示的迭代过程计算各个测试点的终端收敛值记作
Figure GDA0003153790570000183
之后按照终端收敛值从小到大的顺序进行排序,现假设排序情况为:
Figure GDA0003153790570000184
(2)计算质心点
在搜索过程中需要使用到测试点的质心点,质心点的计算方式如下:
Figure GDA0003153790570000185
(3)计算反射点
基于上一步得到的质心点计算反射点,其计算过程如下:
pr=p0+α(p0-p3) (75)
其中α=0.5为反射系数,之后按照图三所示的迭代过程计算反射点pr的终端收敛值
Figure GDA0003153790570000186
之后进行如下判断。
(3.1)反射
如果反射点的收敛值
Figure GDA0003153790570000191
小于第二个测试点的收敛值
Figure GDA0003153790570000192
且大于第一个测试点的收敛值
Figure GDA0003153790570000193
此时将pr的值赋予p3并回到第(1)步重新计算排序。
(3.2)扩张
如果反射点的收敛值
Figure GDA0003153790570000194
小于第一个测试点的收敛值
Figure GDA0003153790570000195
即表示找到了一个更好的收敛方向,此时需要继续沿着该方向搜索,之后计算扩展点:
pe=p0+η(pr-p0) (76)
其中η=2为扩展系数,之后计算扩展点的收敛值
Figure GDA0003153790570000196
如果
Figure GDA0003153790570000197
则将pe的值赋予p3并回到第(1)步重新计算排序,否则将pr的值赋予p3并回到第(1)步重新计算排序。
(3.3)紧缩
如果反射点收敛值满足
Figure GDA0003153790570000198
则计算紧缩点:
pc=p0+ρ(p3-p0) (77)
其中ρ=0.5为紧缩系数,如果紧缩点的收敛值满足
Figure GDA0003153790570000199
则将pc的值赋予p3并回到第(1)步重新计算排序。
(3.4)收缩
如果以上条件均不满足,则证明此时的测试点已经很接近最优点了,此时需要沿着当前最优点的方向进行收缩,即将第二和第三个测试点按照如下方式进行替换:
pi=p1+σ(pi-p1),i=2,3 (78)
其中σ=0.5为收缩系数,替换完成后回到第(1)步重新计算排序。
通过上述讨论可知,基于Nelder-Mead的搜索过程是一个给定初始值后反复迭代的搜索过程,因此需要设定迭代过程的相关参数。在本发明中的迭代终止条件设定为迭代次数达到上限值或三个测试点的标准差小于设定阈值。其中,迭代次数上限为100次,标准差阈值设定为0.05,迭代的初始测试点给定为[0,0;1,0;0,1]T,之后通过上述搜索过程便可找出最优初始值λ*(k),进而求得最优控制量u*(k)。
仿真实验验证与对比
为了验证本发明所设计的面向后轮主动转向车辆的快速实时预测控制算法的有效性,对行驶在低附着路面上后轮主动转向车辆进行稳定性控制。为了进行对比,同样的控制问题还是用了传统的IPOPT方法进行求解。在仿真验证中的各个参数设定为质量m=1430kg,车辆前半轴距Lf=1.05m,车辆后半轴距Lr=1.61m,车辆绕质心的转动惯量Iz=2059.2kg·m2,前轮等效侧偏刚度Cf=90700N/rad,后轮等效侧偏刚度Cr=109000N/rad,车轮半径Re=0.32m5,车辆车体宽度d=1.5m5,路面附着系数μ=0.5。状态约束分别为γup=0.205ra8d s和βup=0.0376rad,车辆纵向速度V=60km/h,采样时间Ts=0.01s,预测时域N=5,最大后轮转角约束δmax=0.05rad,最大附加横摆力矩ΔMmax=2000Nm。目标函数中的权重系数分别为Γβ=0.2和Γu=[ΓrM]=[100,1/8000]。
双移线工况实验
在仿真实验中,所选择的是双移线工况,车辆速度为60km/h且在整个工况中保持不变。图4和图5分别给出了在不同求解算法下的控制效果图以及相应的控制量,从图中可以看到,在不同的求解算法下均可求得合理的控制量,且在控制器作用下,车辆可以准确地跟踪期望的横摆角速度,提高了车辆的稳定性。图6给出了车辆在不施加控制作用时的效果图,可以看到车辆不能够跟踪上期望的横摆角速度,更进一步证明了本发明所设计的控制器的有效性。图7则给出了质心侧偏角的对比结果,可以发现在控制器的作用下,质心侧偏角可以被抑制在一个更小的范围内,车辆的稳定性有较大的提升。通过以上对比可以证明,本发明所提出的面向后轮主动转向车辆的快速实时预测控制可以满足控制器的控制需求,且在求解精度上与传统的求解算法类似,能够求得满足要求的控制量。
图8给出了在不同预测时域下本发明所设计的快速实时求解算法和传统IPOPT求解算法的平均计算时间,可以看到随着预测时域的增加,本发明所设计的求解算法的平均计算时间和预测时域基本呈线性关系,增长较为缓慢;而传统的IPOPT求解算法的平均计算时间和预测时域呈指数关系,增长较快,且其平均计算时间也大于本发明所设计的求解算法的计算时间。图9和图10分别给出了在不同预测时域下两种求解算法的瞬态求解时间。从图中可以看出,不论是基于传统IPOPT算法还是基于本发明所设计的快速实时求解算法,他们的计算时间都会随着预测时域的增加而增加。但在同样的预测时域下,本发明所设计的快速实时求解算法的求解时间要大大小于传统IPOPT算法的求解时间。当优化问题较为复杂时,例如在第4秒到第11秒期间,本发明所设计的快速实时求解算法的求解时间可比IPOPT算法的求解时间快约十倍。以上两点充分体现出了本发明所述的快速求解算法的有益性。
通过仿真示例可以看出,本发明所设计的面向后轮主动转向车辆的快速实时预测控制算法能够在大幅提高计算速度的情况下还能保证较高的求解精度。从而满足了车载的实时性要求,同时也降低了车载控制器预测控制的实现难度。

Claims (1)

1.一种快速实时的后轮主动转向预测控制方法,其特征在于:其步骤是:
步骤一、后轮主动转向车辆的模型搭建与仿真工况的构建:在CarSim软件中选择车辆模型,并将车辆的运动状态参数读取到Simulink中,基于所选的车辆模型构建低附着路面下行驶的仿真工况,模拟实际车辆的横摆运动和侧向运动特征;
步骤二、模型预测控制器设计:
1)建立考虑后轮主动转向的车辆动力学模型;
2)建立能够产生期望状态参考值的参考模型;
3)根据控制器设计中跟踪期望横摆角速度、抑制质心侧偏角以及限制控制作动能量的需求,并结合相应的约束和车辆动力学模型构建目标函数;
步骤三、优化问题快速求解算法设计
1)引入松弛函数,将初始的带约束的非线性优化问题转化为无约束非线性优化问题;
2)根据目标函数和车辆动力学模型构建哈密顿函数;
3)根据PMP原理推导出沿最优轨线的递推关系式,以及最优终端必要性条件;
4)根据Nelder-Mead单纯形算法原理设计最优协态变量初值搜索算法;
上述步骤详细过程如下:
1)车辆动力学模型
根据车辆动力学的理论,简化后的车辆二自由度模型可由如下方程描述:
Figure FDA0003153790560000011
其中,
Figure FDA0003153790560000012
Figure FDA0003153790560000013
分别表示车辆的质心侧偏角的导数和车辆的横摆角速度的导数,V表示车辆的纵向速度,Fyf和Fyr则分别表示前后轮胎的轮胎侧向力,Lf和Lr分别表示前后轴到车辆质心的距离,m是车辆的质量,Iz是车辆绕质心旋转的转动惯量,ΔMz为附加的横摆力矩;
2)前轮轮胎模型建立
当前轮侧偏角αf很小时,有tan(αf)≈αf,之后该前轮轮胎模型可近似为:
Figure FDA0003153790560000014
其中μ为路面附着系数,Fz为垂直载荷,Cf为前轮侧偏刚度;
式中的前轮侧偏角αf可由下式进行计算:
Figure FDA0003153790560000015
其中δf车辆的前轮转角;
3)后轮轮胎模型建立
考虑路面的附着系数的影响,该模型可描述为:
Figure FDA0003153790560000021
其中的μ为路面附着系数,Fz为垂直载荷,Cr为后轮侧偏刚度;
式中的后轮侧偏角αr可由下式进行计算:
Figure FDA0003153790560000022
其中的δr为后轮转角;
4)车辆参考模型建立
前轮转角δf到车辆的期望横摆角速度γref以及期望质心侧偏角βref的传递函数为:
Figure FDA0003153790560000023
定义车辆的稳定因子为
Figure FDA0003153790560000024
其中L=Lf+Lr为车辆的前后轴距;
系统的振荡频率为
Figure FDA0003153790560000025
系统的阻尼系数为
Figure FDA0003153790560000026
横摆角速度稳态增益为
Figure FDA0003153790560000027
质心侧偏角稳态增益为
Figure FDA0003153790560000028
其中的横摆角速度微分系数定义为
Figure FDA0003153790560000029
质心侧偏角微分系数定义为
Figure FDA00031537905600000210
定义横摆角速度的上限值为
Figure FDA00031537905600000211
横摆角速度的参考值γref应被约束在|γref|≤γup,定义质心侧偏角的上限值为
Figure FDA00031537905600000212
横摆角速度的参考值βref应被约束在|βref|≤βup
5)控制器设计
定义系统的状态向量为x=[x1,x2]T=[β,γ]T,控制量定义为u=[u1,u2]=[δr,ΔMz],定义ΔMmax和δmax分别是附加横摆力矩和后轮主动转角的最大约束,将(1)式中的状态空间方程进行欧拉离散化得在每个采样时刻kTs离散后的状态空间方程为:
Figure FDA0003153790560000031
其中的前轮侧向力Fyf由式(2)中的前轮轮胎模型计算得到,后轮侧向力Fyr则由式(4)中的后轮轮胎模型计算得到,βup和γup分别为质心侧偏角和横摆角速度的上限值,ΔMmax为附加横摆力矩的最大值;
在每个时刻k+1≤ki≤k+N+1定义代价函数为:
Figure FDA0003153790560000032
其中L1(ki)、L'2(ki)和L3(ki)分别用于跟踪横摆角速度参考值、跟踪质心侧偏角参考值和抑制作动能量;
状态约束为|x1(ki)|≤βup,控制量限幅为|u1(ki-1)|≤δmax以及|u2(ki-1)|≤ΔMmax
得到非线性模型预测控制的目标函数为:
Figure FDA0003153790560000033
分别为质心侧偏角和控制量的权重系数,N为预测时域;
6)基于PMP原理的求解
由PMP原理可知,对于此类的问题得到:
Figure FDA0003153790560000034
若存在控制量u*(t)和以及终端时间
Figure FDA0003153790560000035
使得性能指标中的目标函数在该控制量的作用下取到最小的最优解,且在该最优解的作用下最优状态轨线x*(t)所形成的最优轨迹也存在,则可导出如下必要条件:
(1)x(t)和λ(t)满足正则方程:
Figure FDA0003153790560000041
其中λ=[λ12]T为拉格朗日乘子的向量函数,L(x,u)为目标函数中的各项,H(x,λ,u)=L(x,u)+λT(t)f(x,u)为哈密顿函数;
(2)状态轨线x(t)和协态变量λ(t)满足边界条件:
Figure FDA0003153790560000042
(3)哈密顿函数H(x,λ,u)在最优控制量u*(t)的作用下取到绝对极小值:
Figure FDA0003153790560000043
(4)同时哈密顿函数在最优状态轨线的末端满足:
H[x*(t),λ(t),u*(t)]=H[x*(tf),λ(tf),u*(tf)]=constant (14)
由上述的状态约束|x1(ki)|≤1的基础上,引入终端松弛因子对状态约束进行处理:
Figure FDA0003153790560000044
其中的κ代表了函数的松弛程度,而ν则为一个较大的数,
最终得到目标函数的子项L2(ki)为:
L2(ki)=L′2(ki)+ζ(ki) (16)
此时式(10)中的目标函数则转化为:
Figure FDA0003153790560000045
该目标函数可满足状态约束|x1(ki)|≤1;
定义在时刻k+1≤ki≤k+N+1的哈密顿函数为:
Figure FDA0003153790560000046
其中F1(x(ki))和F2(x(ki))分别定义为:
Figure FDA0003153790560000051
在上述方程中,λ1(k)和λ2(k)分别表示拉格朗日乘子,根据庞特里亚金最小值原理可得到最优必要性条件为:
Figure FDA0003153790560000052
终端条件为:
Figure FDA0003153790560000053
在各个时刻都有最优控制律u*(ki)使得哈密顿函数最小化:
Figure FDA0003153790560000054
将哈密顿函数重新整理为关于二维控制量u=[u1,u2]T的二元二次函数:
Figure FDA0003153790560000055
其中g[x(ki),δf(k)]为与当前控制量无关的余项,其定义如下:
Figure FDA0003153790560000056
哈密顿函分别对两个控制量求偏导可得:
Figure FDA0003153790560000061
Figure FDA0003153790560000062
分别令上述两个偏导等于零便可求得哈密顿函数的驻点为
Figure FDA0003153790560000063
以及
Figure FDA0003153790560000064
对哈密顿函数相对于控制量求二阶偏导,其二阶偏导如下所示:
Figure FDA0003153790560000065
Figure FDA0003153790560000066
Figure FDA0003153790560000067
Figure FDA0003153790560000068
得AC-B2>0,且A>0,C>0可知哈密顿函数在驻点处取得极小值;
根据在驻点处的结果,便可给出使哈密顿函数取得极小值的最优控制律:
Figure FDA0003153790560000069
Figure FDA00031537905600000610
其中的第一个控制量u1为后轮主动转角,该控制量直接施加给后轮转向执行机构,第二个控制量u2为附加横摆力矩,该控制量首先需要通过如下转换:
Figure FDA0003153790560000071
其中ΔTfl(k),ΔTrl(k),ΔTfr(k),ΔTrr(k)分别表示左前轮、左后轮、右前轮以及右后轮的附加转矩,Re表示轮胎的滚动半径,d表示车辆车身的宽度,之后该控制量便转换为四个轮胎的附加转矩施加给车辆;
最优控制问题此时转化为了两点边值问题:
Figure FDA0003153790560000072
该两点边值问题通过Nelder-Mead算法进行求解;
7)基于Nelder-Mead的最优初值搜索
(1)检查给定初始测试点的收敛情况
三个初始测试点p1,p2,p3,每个测试点都是二维向量;各个测试点的终端收敛值记作
Figure FDA0003153790560000073
之后按照终端收敛值从小到大的顺序进行排序,现假设排序情况为:
Figure FDA0003153790560000074
(2)计算质心点
Figure FDA0003153790560000075
(3)计算反射点
pr=p0+α(p0-p3) (36)
其中α=0.5为反射系数,之后按照图三所示的迭代过程计算反射点pr的终端收敛值
Figure FDA0003153790560000076
之后进行如下判断;
(3.1)反射
如果反射点的收敛值
Figure FDA0003153790560000077
小于第二个测试点的收敛值
Figure FDA0003153790560000078
且大于第一个测试点的收敛值
Figure FDA0003153790560000079
此时将pr的值赋予p3并回到第(1)步重新计算排序;
(3.2)扩张
如果反射点的收敛值
Figure FDA00031537905600000710
小于第一个测试点的收敛值
Figure FDA00031537905600000711
即表示找到了一个更好的收敛方向,此时需要继续沿着该方向搜索,之后计算扩展点:
pe=p0+η(pr-p0) (37)
其中η=2为扩展系数,之后计算扩展点的收敛值
Figure FDA00031537905600000712
如果
Figure FDA00031537905600000713
则将pe的值赋予p3并回到第(1)步重新计算排序,否则将pr的值赋予p3并回到第(1)步重新计算排序;
(3.3)紧缩
如果反射点收敛值满足
Figure FDA0003153790560000081
则计算紧缩点:
pc=p0+ρ(p3-p0) (38)
其中ρ=0.5为紧缩系数,如果紧缩点的收敛值满足
Figure FDA0003153790560000082
则将pc的值赋予p3并回到第(1)步重新计算排序;
(3.4)收缩
如果以上条件均不满足,则证明此时的测试点已经很接近最优点了,此时需要沿着当前最优点的方向进行收缩,即将第二和第三个测试点按照如下方式进行替换:
pi=p1+σ(pi-p1),i=2,3 (39)
其中σ=0.5为收缩系数,替换完成后回到第(1)步重新计算排序。
CN202011114070.1A 2020-10-18 2020-10-18 快速实时的后轮主动转向预测控制方法 Active CN112298193B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011114070.1A CN112298193B (zh) 2020-10-18 2020-10-18 快速实时的后轮主动转向预测控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011114070.1A CN112298193B (zh) 2020-10-18 2020-10-18 快速实时的后轮主动转向预测控制方法

Publications (2)

Publication Number Publication Date
CN112298193A CN112298193A (zh) 2021-02-02
CN112298193B true CN112298193B (zh) 2021-08-24

Family

ID=74327658

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011114070.1A Active CN112298193B (zh) 2020-10-18 2020-10-18 快速实时的后轮主动转向预测控制方法

Country Status (1)

Country Link
CN (1) CN112298193B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113867330B (zh) * 2021-05-11 2023-05-02 吉林大学 一种基于多自由度预测模型在任意路径下实现车辆漂移的控制方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1920981A1 (de) * 2006-11-08 2008-05-14 Bayerische Motoren Werke Aktiengesellschaft Verfahren zum Ermitteln der von einem Fahrzeugreifen übertragenen Seitenkräfte oder Längskräfte mit einem sogenannten Reifenmodell
CN102707282A (zh) * 2012-06-08 2012-10-03 北京航空航天大学 一种对非协同目标高精度的isar运动补偿方法
CN107117170A (zh) * 2017-04-28 2017-09-01 吉林大学 一种基于经济性驾驶的实时预测巡航控制系统
CN111158264A (zh) * 2020-01-09 2020-05-15 吉林大学 面向车载应用的模型预测控制快速求解方法
CN111267834A (zh) * 2020-02-17 2020-06-12 北京理工大学 一种车辆横摆稳定预测控制方法及系统

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1920981A1 (de) * 2006-11-08 2008-05-14 Bayerische Motoren Werke Aktiengesellschaft Verfahren zum Ermitteln der von einem Fahrzeugreifen übertragenen Seitenkräfte oder Längskräfte mit einem sogenannten Reifenmodell
CN102707282A (zh) * 2012-06-08 2012-10-03 北京航空航天大学 一种对非协同目标高精度的isar运动补偿方法
CN107117170A (zh) * 2017-04-28 2017-09-01 吉林大学 一种基于经济性驾驶的实时预测巡航控制系统
CN111158264A (zh) * 2020-01-09 2020-05-15 吉林大学 面向车载应用的模型预测控制快速求解方法
CN111267834A (zh) * 2020-02-17 2020-06-12 北京理工大学 一种车辆横摆稳定预测控制方法及系统

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
Design and Experimental Verification of Real-Time Nonlinear Predictive;Ping Wang el;《IEEE TRANSACTIONS ON CONTROL SYSTEMS TECHNOLOGY》;20200916;全文 *

Also Published As

Publication number Publication date
CN112298193A (zh) 2021-02-02

Similar Documents

Publication Publication Date Title
CN111258323B (zh) 一种智能车辆轨迹规划与跟踪的联合控制方法
Min et al. Trajectory optimization of an electric vehicle with minimum energy consumption using inverse dynamics model and servo constraints
CN111497826B (zh) 一种电动汽车横摆稳定性协调控制方法及系统
CN111158264B (zh) 面向车载应用的模型预测控制快速求解方法
CN111332277B (zh) 极限工况下基于稳定区域的车辆侧向稳定性控制方法
CN112622903B (zh) 一种车辆跟随驾驶环境下自主车辆的纵向和横向控制方法
CN107561942A (zh) 基于模型补偿的智能车辆轨迹跟踪模型预测控制方法
CN111391822B (zh) 一种极限工况下汽车横纵向稳定性协同控制方法
Singh et al. Trajectory tracking and integrated chassis control for obstacle avoidance with minimum jerk
CN112298193B (zh) 快速实时的后轮主动转向预测控制方法
CN116184835A (zh) 考虑输入时延与速度时变的无人驾驶车辆跟踪控制方法
Li et al. Integrated coordination control for distributed drive electric vehicle trajectory tracking
Gu et al. Genetic algorithm based LQR control for AGV path tracking problem
Júnior et al. Vehicle stability upper-level-controller based on parameterized model predictive control
Chu et al. Autonomous high-speed overtaking of intelligent chassis using fast iterative model predictive control
WO2024114037A1 (zh) 智能电动汽车路径跟踪控制方法及装置
Kanchwala et al. Improving handling performance of an electric vehicle using model predictive control
Wohner et al. Hierarchical nonlinear model predictive control for an autonomous racecar
Jan et al. Decoupling of vehicle lateral dynamics using four-wheel steering system
Li et al. Path tracking control based on the prediction of tire state stiffness using the optimized steering sequence
CN116279409A (zh) 一种四轮独立驱动与转向电动汽车的协同控制方法
Li et al. Time-Optimal Trajectory Planning and Tracking for Autonomous Vehicles
Zhang et al. Low-speed vehicle path-tracking algorithm based on model predictive control using QPKWIK solver
Wang et al. Trajectory Tracking Control of Four-Wheel Independent Drive Electric Vehicles
Sill et al. Cascaded predictive control of tire force saturation levels for vehicle stability

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