一种四轮转向联合差动制动的智能车辆分层控制方法

文档序号:1636413 发布日期:2020-01-17 浏览:35次 >En<

阅读说明:本技术 一种四轮转向联合差动制动的智能车辆分层控制方法 (Intelligent vehicle layered control method combining four-wheel steering and differential braking ) 是由 朱梁鹏 刘成晔 赵景波 朱燕 于 2019-09-19 设计创作,主要内容包括:本发明公开一种四轮转向联合差动制动的无人驾驶车辆分层控制方法,以非线性模型预测(MPC)算法设计联合控制器,对车辆进行四轮转向和差动制动控制。该方法具有三层控制结构,上层为轨迹规划层,实现对车辆纵向、横向以及横摆角跟踪,输出理想横摆角速度和质心侧偏角;中间层为MPC集中控制层,实现车辆后轮转角控制和附加横摆力矩控制;下层为制动力分配器,通过附加横摆力矩的计算,转化为对单个车轮的制动控制;最终,MPC联合控制器输出的四轮转角和差动制动力作用于整车,实现目标路径的稳定性跟踪。本方法利用MPC控制方法滚动优化的特性,在提高跟踪精度的基础上,实现了四轮转向和差动制动的联合控制,提高了高速车况下智能汽车轨迹跟踪的稳定性。(The invention discloses a four-wheel steering and differential braking combined unmanned vehicle hierarchical control method, which designs a combined controller by a nonlinear Model Prediction (MPC) algorithm to control four-wheel steering and differential braking of a vehicle. The method has a three-layer control structure, the upper layer is a track planning layer, so that the longitudinal, transverse and yaw angles of the vehicle are tracked, and an ideal yaw angular speed and a centroid side slip angle are output; the middle layer is an MPC centralized control layer to realize the control of the rear wheel steering angle and the control of the additional yaw moment of the vehicle; the lower layer is a braking force distributor which is converted into braking control on a single wheel through the calculation of an additional yaw moment; finally, the MPC is combined with the four-wheel corner and the differential braking force output by the controller to act on the whole vehicle, so as to realize the stability tracking of the target path. The method utilizes the rolling optimization characteristic of the MPC control method, realizes the combined control of four-wheel steering and differential braking on the basis of improving the tracking precision, and improves the stability of intelligent automobile track tracking under the condition of high-speed automobiles.)

一种四轮转向联合差动制动的智能车辆分层控制方法

技术领域

本发明涉及到智能车辆控制领域,特别是一种四轮转向联合差动制动的无人驾驶车辆轨迹跟踪分层控制方法。

背景技术

目前无人驾驶车辆的路径跟踪控制大多通过前轮转向来进行横向的跟踪控制,由于车辆模型的强非线性导致车辆在高速工况时仅仅依靠前轮转向控制是无法保证车辆的稳定性,而且过于集成的控制算法设计及应用难度太大,因此综合考虑四轮转向联合差动制动分层控制是解决无人驾驶车辆路径跟踪稳定性的有效途径。

发明内容

为解决无人驾驶车辆路径跟踪控制精度以及车辆稳定性问题,本发明提供一种四轮转向联合差动制动的分层控制方法,其目的在于:以非线性模型预测(MPC)算法设计联合控制器,对车辆进行四轮转向和差动制动控制,实现高速车况下对期望路径的稳定跟踪。

为实现上述目的,本发明采用的技术方案是:

一种四轮转向联合差动制动的分层控制方法,具有三层控制结构,具体为:

上层:轨迹规划层,实现对车辆纵向、横向以及横摆角跟踪,输出理想横摆角速度和质心侧偏角:;

中间层:MPC集中控制层,实现车辆后轮转角控制和附加横摆力矩控制;

下层:制动力分配器,通过附加横摆力矩的计算,转化为对单个车轮的制动控制;

最终,MPC联合控制器输出的四轮转角和差动制动力作用于整车,实现目标路径的稳定性跟踪。

具体地:

第一层:基于运动学MPC控制

1建立上层轨迹规划层MPC控制器,实现前轮转向控制

1.1根据车辆单轨模型,建立车辆运动学模型,并将其写成状态空间表达形式为:

Figure BDA0002207863810000011

其中状态量为:

Figure BDA0002207863810000021

式中:X为车辆纵向位移,Y为车辆横向位移,

Figure BDA0002207863810000022

为车辆横摆角;

控制量选取为:

u=[v,δf]T

式中:v为车辆后轴中心车速,δf为前轮偏转角

1.2将公式所示的车辆运动学非线性模型进行线性处理,得到线性时变系统为:

Figure BDA0002207863810000023

1.3将公式采用一阶差商的方法进行离散化处理,得到离散的状态空间表达式:

Figure BDA0002207863810000024

其中:Ak=I+A*T;Bk=I+B*T;

1.4设车辆轨迹跟踪过程某一时刻为t时刻,为求t时刻的前轮转角,建立如下目标函数:

Figure BDA0002207863810000025

式中:NP为预测时域;Nc为控制时域;ρ为权重系数;ε为松弛因子

1.5公式采用控制增量取代控制量进行系统的控制,因此将公式改写为控制增量式的状态空间方程:

Figure BDA0002207863810000026

Figure BDA0002207863810000027

其中:

Figure BDA0002207863810000028

1.6对控制过程中的控制量和控制增量进行约束:

umin(k+1)<u<umax(k+1)

Δumin(k+1)<Δu<Δumax(k+1)

1.7对目标函数求解后得到控制时域内的一系列控制输入增量

Figure BDA0002207863810000029

将t时刻该控制序列中的第一个元素作为控制输入增量作用到系统,即:

u(t)=u(t-1)+Δu(t)

1.8在t+1时刻重复步骤,完成对期望轨迹的跟踪,得到一系列的控制量v,δf第二层:理想跟随模型

2.建立理想车辆跟随模型,计算出轨迹跟踪下的理想质心侧偏角和横摆角速度

2.1车辆理想质心侧偏角模型

车辆稳态时的理想质心侧偏角为:

Figure BDA0002207863810000031

式中:L为轴距;lf为质心到前轮距离;lr为质心到后轮距离;kr为后轮侧偏刚度轮胎最大附着极限下的质心侧偏角为:

式中:μ为附着率;g为重力加速度

理想的车辆质心侧偏角为:

βd=min(|βideal|,|βmax|)·sgn(βideal)

2.2车辆理想横摆角速度模型

Figure BDA0002207863810000033

轮胎附着极限下的最大横摆加速度为:

Figure BDA0002207863810000034

理想的车辆横摆角速度为:

rd=min(|rideal|,|rmax|)·sgn(rideal)

第三层:中间控制层

3.建立中间层MPC控制器,计算制动力产生的附加横摆力矩,实现后轮转角控制

3.1建立车辆四轮转向和差动制动联合动力学模型,并改写为状态空间表达形式:

Figure BDA0002207863810000036

其中状态量为:

x=[β,r]T

式中:β为车辆质心侧偏角,r为车辆横摆角速度,

Figure BDA0002207863810000035

为前轮偏转角;

控制量选取为:

u=[δr,Ms]T

式中:Ms为制动产生的附加横摆力矩,δr为后轮偏转角

3.2将公式所示的车辆运动学非线性模型进行线性处理,得到线性时变系统为:

Figure BDA0002207863810000041

3.3将公式采用一阶差商的方法进行离散化处理,得到离散的状态空间表达式:

其中:Ak=I+A*T;Bk=I+B*T;

3.4设车辆轨迹跟踪过程某一时刻为t时刻,为求t时刻的前轮转角,建立如下目标函数:

Figure BDA0002207863810000043

式中:NP为预测时域;Nc为控制时域;ρ为权重系数;ε为松弛因子;

3.5公式采用控制增量取代控制量进行系统的控制,因此将公式改写为控制增量式的状态空间方程:

Figure BDA0002207863810000044

其中:

3.6对控制过程中的控制量和控制增量进行约束:

umin(k+1)<u<umax(k+1)

Δumin(k+1)<Δu<Δumax(k+1)

3.7对目标函数求解后得到控制时域内的一系列控制输入增量

Figure BDA0002207863810000047

将t时刻该控制序列中的第一个元素作为控制输入增量作用到系统,即:

u(t)=u(t-1)+Δu(t)

3.8在t+1时刻重复步骤,完成对期望轨迹的跟踪,得到一系列的控制量Msr第四层:制动力分配层

4.1车辆转向情况判别

Figure BDA0002207863810000048

式中:

Figure BDA0002207863810000049

当K1=0,

Figure BDA0002207863810000051

中性转向:K=0,

Figure BDA0002207863810000052

为中值;

过度转向:K<0,

Figure BDA0002207863810000053

增大;

不足转向:K>0,

Figure BDA0002207863810000054

减小;

4.1单轮制动压力计算

Figure BDA0002207863810000055

4.2制动力施加决策

Figure BDA0002207863810000056

相对于现有技术,本发明具有如下有益效果:

1.本发明设计了一种考虑车辆稳定性的四轮转向联合差动制动的智能车辆轨迹跟踪分层控制方法,通过上层轨迹规划层对期望路径进行跟踪,中间层利用上层规划出的前轮转角、期望横摆角速度和质心侧偏角优化计算出后轮转角和附加横摆力矩控制量,下层控制器对制动力进行分配,实现无人驾驶车辆四轮转角和差动制动的主动调节,改善车辆轨迹跟踪的稳定性。

2.本发明将四轮转向控制纳入控制算法之中,改善了车辆高速工况下的操纵性和稳定性。

3.本发明实现了轨迹跟踪的差动制动控制,能够主动调节汽车后轮转向角和制动力矩,从而优化汽车的质心侧偏角、横摆角速度和侧向位移,提高了智能车辆在高速工况下的轨迹跟踪稳定性。

4.本发明采用模型预测控制器作为控制算法,对控制目标进行实时预测和滚动优化,处理多约束优化问题,极大提高了轨迹跟踪的稳定性。

附图说明

通过阅读参照以下附图对非限制性实施例所作的详细描述,本发明的其它特征、目的和优点将会变得更明显:

图1为本发明车辆运动学模型图;

图2为本发明一种四轮转向联合差动制动的智能车辆分层控制算法框架图;

图3为本发明控制策略仿真图;

图4为本发明轨迹跟踪仿真结果图;

图5为本发明四轮转角仿真结果图;

图6为本发明四轮纵向力仿真结果图。

具体实施方式

为了使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例,对本发明进行进一步详细说明。应当理解,此处所描述的具体实施例仅仅用以解释本发明,并不用于限定本发明。

参阅图1-5,本发明一种四轮转向联合差动制动的无人驾驶车辆分层控制方法,该控制方法具有三层控制结构。上层为轨迹规划层,在实现对车辆纵向、横向以及横摆角跟踪,输出理想横摆角速度和质心侧偏角。中间层为MPC集中控制层,实现车辆后轮转角控制和附加横摆力矩控制。下层为制动力分配器,通过附加横摆力矩的计算,转化为对单个车轮的制动控制。最终,MPC联合控制器输出的四轮转角和差动制动力作用于整车,实现目标路径的稳定性跟踪。

第一层:上层轨迹跟踪控制层

下面结合图1对上层模型预测轨迹跟踪控制器设计进行详细说明:

1.车辆运动学建模

Figure BDA0002207863810000061

Figure BDA0002207863810000062

Figure BDA0002207863810000063

式中:l为车辆轴距;v为车辆后轴中心速度;δf为前轮偏转角

2.根据公式,可以将其推导成状态空间表达形式为:

Figure BDA0002207863810000064

其中状态量为:

式中:X为车辆纵向位移,Y为车辆横向位移,

Figure BDA0002207863810000072

为车辆横摆角;

控制量选取为:

u=[v,δf]T

式中:v为车辆后轴中心车速,δf为前轮偏转角

3.建立线性误差方程

3.1由于车辆模型具有非线性特点,需要对模型进行线性化处理,由于期望输入无法给出所有状态点的信息,因此本发明采用近似线性化方法对模型进行线性化,对系统施加某个工作点的常数控制量,根据该工作点的轨迹和实际状态的偏差设计控制算法,考虑系统的某个工作点为[x0,u0],x0为始终施加控制量为u0后得到的系统状态量,则

x0(k+1)=f(x0(k),u0)

x0(0)=x0

经过推导可以得到:

Figure BDA0002207863810000073

式中:

Figure BDA0002207863810000074

3.2将公式采用一阶差商的方法进行离散化处理,得到离散的状态空间表达式:

Figure BDA0002207863810000075

y(k)=Ckx(k)

其中:

Figure BDA0002207863810000077

Figure BDA0002207863810000078

3.3公式采用控制增量取代控制量进行系统的控制,因此将公式改写为控制增量式的状态空间方程:

Figure BDA0002207863810000079

Figure BDA0002207863810000081

其中:

Figure BDA0002207863810000082

3.4经过推导,可以得到系统的预测输出表达式:

Figure BDA0002207863810000089

Figure BDA0002207863810000083

Figure BDA0002207863810000085

Figure BDA0002207863810000087

4.建立滚动时域优化控制器

4.1针对轨迹跟踪控制问题,建立如下控制目标函数

Figure BDA0002207863810000088

式中:NP为预测时域;Nc为控制时域;ρ为权重系数;ε为松弛因子;Q,R为权重矩阵。

4.2为了增加轨迹跟踪的平稳性,需要对控制量、控制增量进行约束:

umin(t+k)≤u(t+k)≤umax(t+k)

Δumin(t+k)≤Δu(t+k)≤Δumax(t+k)

式中:Δu为控制增量,umin为控制增量极小值,umax为控制增量极大值本发明将控制量极限设置为umin=[-0.2-0.6]T,umax=[0.2 0.6]T,Δumin=[-0.34-0.02]T,Δumax=[0.02]TYmax=[-3-0.3]T,Ymin=[3 0.3]T

4.3将目标函数转化为标准二次型形式并结合约束条件,解决以下优化问题:

J(ξ(t),u(t-1),ΔU(t))=[ΔU(t)T,ε]Ht[ΔU(t)T,ε]+Gt[ΔU(t)T,ε]

Figure BDA0002207863810000091

Figure BDA0002207863810000092

Figure BDA0002207863810000093

式中:

Figure BDA0002207863810000094

et为预测时域内的跟踪误差

4.3对目标函数求解后得到控制时域内的一系列控制输入增量

Figure BDA0002207863810000095

将t时刻该控制序列中的第一个元素作为控制输入增量作用到系统,即:

u(t)=u(t-1)+Δu(t)

在t+1时刻重复步骤,完成对期望轨迹的跟踪,得到一系列的控制量v,δf

第二层:理想跟随模型

2.建立理想车辆跟随模型,计算出轨迹跟踪下的理想质心侧偏角和横摆角速度

2.1车辆理想质心侧偏角模型

车辆稳态时的理想质心侧偏角为:

Figure BDA0002207863810000096

式中:L为轴距;lf为质心到前轮距离;lr为质心到后轮距离;kr为后轮侧偏刚度轮胎最大附着极限下的质心侧偏角为:

Figure BDA0002207863810000097

式中:μ为附着率;g为重力加速度

理想的车辆质心侧偏角为:

βd=min(|βideal|,|βmax|)·sgn(βideal)

2.2车辆理想横摆角速度模型

Figure BDA0002207863810000101

轮胎附着极限下的最大横摆加速度为:

Figure BDA0002207863810000102

理想的车辆横摆角速度为:

rd=min(|rideal|,|rmax|)·sgn(rideal)

第三层:中间控制层

1.车辆运动学建模

Figure BDA0002207863810000103

式中:I为转动惯量;lf,lr为前后轴距;kf,kr为前后轮侧偏刚度

2.根据公式,可以将其推导成状态空间表达形式为:

Figure BDA0002207863810000105

其中状态量为:

x=[β,r]T

式中:β为车辆质心侧偏角,r为车辆横摆角速度,为前轮偏转角;

控制量选取为:

u=[δr,Ms]T

式中:Ms为制动产生的附加横摆力矩,δr为后轮偏转角

3.建立线性误差方程

3.1由于车辆模型具有非线性特点,需要对模型进行线性化处理,由于期望输入无法给出所有状态点的信息,因此本发明采用近似线性化方法对模型进行线性化,对系统施加某个工作点的常数控制量,根据该工作点的轨迹和实际状态的偏差设计控制算法,考虑系统的某个工作点为[x0,u0],x0为始终施加控制量为u0后得到的系统状态量,则

x0(k+1)=f(x0(k),u0)

x0(0)=x0

经过推导可以得到:

Figure BDA0002207863810000111

式中:

Figure BDA0002207863810000112

3.2将公式采用一阶差商的方法进行离散化处理,得到离散的状态空间表达式:

Figure BDA0002207863810000113

其中:

Figure BDA0002207863810000114

Figure BDA0002207863810000115

3.3公式采用控制增量取代控制量进行系统的控制,因此将公式改写为控制增量式的状态空间方程:

Figure BDA0002207863810000116

Figure BDA0002207863810000117

其中:

Figure BDA0002207863810000118

公式采用控制增量取代控制量进行系统的控制,因此将公式改写为控制增量式的状态空间方程:

Figure BDA0002207863810000119

Figure BDA00022078638100001110

其中:

Figure BDA00022078638100001111

3.4经过推导,可以得到系统的预测输出表达式:

Figure BDA00022078638100001112

Figure BDA00022078638100001113

Figure BDA0002207863810000121

Figure BDA0002207863810000122

Figure BDA0002207863810000124

4.建立滚动时域优化控制器

4.1针对轨迹跟踪控制问题,建立如下控制目标函数

Figure BDA0002207863810000125

式中:NP为预测时域;Nc为控制时域;ρ为权重系数;ε为松弛因子;Q,R为权重矩阵。

4.2为了增加轨迹跟踪的平稳性,需要对控制量、控制增量进行约束:

umin(t+k)≤u(t+k)≤umax(t+k)

Δumin(t+k)≤Δu(t+k)≤Δumax(t+k)

式中:Δu为控制增量,umin为控制增量极小值,umax为控制增量极大值本发明将控制量极限设置为:

umin=[-0.6 -100]T,umax=[0.2 100]T,Δumin=[-0.34 -5]T,Δumin=[0.34 5]T,Ymax=[-1.5 -10]T,Ymin=[1.5 10]T

4.3将目标函数转化为标准二次型形式并结合约束条件,解决以下优化问题:

J(ξ(t),u(t-1),ΔU(t))=[ΔU(t)T,ε]Ht[ΔU(t)T,ε]+Gt[ΔU(t)T,ε]

Figure BDA0002207863810000132

Figure BDA0002207863810000133

式中:et为预测时域内的跟踪误差

4.3对目标函数求解后得到控制时域内的一系列控制输入增量

Figure BDA0002207863810000135

将t时刻该控制序列中的第一个元素作为控制输入增量作用到系统,即:

u(t)=u(t-1)+Δu(t)

在t+1时刻重复步骤,完成对期望轨迹的跟踪,得到一系列的控制量δr,Ms

第四层:制动力分配层

4.1车辆转向情况判别

Figure BDA0002207863810000136

式中:

Figure BDA0002207863810000137

当K1=0,

Figure BDA0002207863810000138

中性转向:K=0,

Figure BDA0002207863810000139

为中值;

过度转向:K<0,

Figure BDA00022078638100001310

增大;

不足转向:K>0,

Figure BDA00022078638100001311

减小;

4.1单轮制动压力计算

Figure BDA00022078638100001312

4.2制动力施加决策

Figure BDA00022078638100001313

Figure BDA0002207863810000141

本发明方案所公开的技术手段不仅限于上述实施方式所公开的技术手段,还包括由以上技术特征任意组合所组成的技术方案。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,均属于本发明保护的范围。

20页详细技术资料下载
上一篇:一种医用注射器针头装配设备
下一篇:用于产生影响汽车车轮滑转调节的指令信号的设备和方法

网友询问留言

已有0条留言

还没有人留言评论。精彩留言会获得点赞!

精彩留言,会给你点赞!

技术分类