CN111055837B - 一种基于启发式算法的车辆横向稳定性控制方法 - Google Patents

一种基于启发式算法的车辆横向稳定性控制方法 Download PDF

Info

Publication number
CN111055837B
CN111055837B CN201911190692.XA CN201911190692A CN111055837B CN 111055837 B CN111055837 B CN 111055837B CN 201911190692 A CN201911190692 A CN 201911190692A CN 111055837 B CN111055837 B CN 111055837B
Authority
CN
China
Prior art keywords
automobile
parameters
vehicle
model
moment
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
CN201911190692.XA
Other languages
English (en)
Other versions
CN111055837A (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.)
Guangdong University of Technology
Original Assignee
Guangdong University of Technology
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 Guangdong University of Technology filed Critical Guangdong University of Technology
Priority to CN201911190692.XA priority Critical patent/CN111055837B/zh
Publication of CN111055837A publication Critical patent/CN111055837A/zh
Priority to AU2020102015A priority patent/AU2020102015A4/en
Application granted granted Critical
Publication of CN111055837B publication Critical patent/CN111055837B/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
    • B60W30/00Purposes 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/02Control of vehicle driving stability
    • 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
    • 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/12Estimation 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

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)
  • Control Of Driving Devices And Active Controlling Of Vehicle (AREA)

Abstract

本发明公开了一种基于启发式算法的车辆横向稳定性控制方法,包括下述步骤:采集汽车的参数;需要采集的参数包括固有参数和实时参数,其中:固有参数:包括汽车总质量m,转动惯量Iz,前、后轴到汽车重心的距离lf、lr,以及前、后轮的侧偏刚度Cf、Cr;这些参数是汽车固有的,可提前采集并存储;本发明采用启发式算法来求解输出反馈控制器,避免了传统求解方法中的复杂迭代;考虑汽车纵向速度和轮胎侧偏刚度为不确定参数,本发明所设计的控制器鲁棒性更好,能够计算出合适的横摆力矩,控制质心侧偏角和横摆角速度在一个合适的范围内,提高汽车的操纵性能,改善汽车的横向稳定性能,降低汽车运行过程中发生危险情况的可能。

Description

一种基于启发式算法的车辆横向稳定性控制方法
技术领域
本发明涉及车辆稳定性控制技术领域,具体涉及一种基于启发式算法的车辆横向稳定性控制方法。
背景技术
随着人民生活水平的不断提高,中国的汽车保有量迅速上升,而随着汽车数量的増多,交通安全问题己经成为了一个社会性问题。中国是世界上交通事故率最高的国家之一,这引起了人们对汽车行驶安全性能的高度重视。汽车安全性防护可根据其作用时间的差异分为被动安全和主动安全两方面:被动安全防护设有安全气囊、高强度车身等防护措施,通过减少事故发生时对乘客的冲击伤害来达到防护效果;主动安全技术包括以提高车辆的稳定性为主的底盘动力学控制,最为典型的是车身电子稳定性控制系统和四轮转向系统,其主要根据汽车的动力学原理,对汽车的内部构造和控制方案进行设计,实现对汽车运行姿态进行检测、修正、报警等方式来主动预防事故的产生,其中ABS、ASR等技术是主动安全技术典型的代表。早期,在汽车安全方面的研究主要集中在被动安全性方面,但随着社会发展,汽车越来越多,交通压力也日益变大,被动安全性方面的技术已不能很好的满足汽车安全性的要求,而以预防为核心的主动安全性己经成为如今交通的迫切要求,且已经成为了本世纪研究汽车安全性能的重点。汽车直接横摆力矩控制(简称DYC,Direct Yaw-momentControl)的研究就是旨在淡化驾驶人员的操作技能对车辆运动安全性的影响,在车辆的各种行驶状态下通过对每个车轮的受力进行调节,产生横摆力矩,克服过多转向或不足转向,从而主动地对车辆进行动力学控制,提高汽车在高速和恶劣道路等极限条件下行驶时的操纵稳定性。
现有技术中,中国专利CN107215329A公开了一种基于ATSM的分布式驱动电动汽车横向稳定性控制方法,该方法也是针对汽车横向动力学稳定性的控制方法,所用方法为二阶滑模变结构算法中的螺旋滑模控制算法,但是,没有考虑汽车参数不确定性因素,因此所设计的控制器在汽车速度发生较大变化时并不能保证稳定性。
发明内容
本发明的目的在于克服现有技术的缺点与不足,提供一种基于启发式算法的车辆横向稳定性控制方法,该方法能够提高汽车的操纵性能,改善汽车的横向稳定性能,降低汽车运行过程中发生危险情况的可能。
本发明的目的通过下述技术方案实现:
一种基于启发式算法的车辆横向稳定性控制方法,包括下述步骤:
步骤一,采集汽车的参数;
需要采集的参数包括固有参数和实时参数,其中:
固有参数:包括汽车总质量m,转动惯量Iz,前、后轴到汽车重心的距离lf、lr,以及前、后轮的侧偏刚度Cf、Cr;这些参数是汽车固有的,可提前采集并存储;
实时参数:包括前轮转向角δ,汽车纵向速度vx,横摆角速度r,质心侧偏角β;这些参数中,前轮转向角可通过传感器采集方向盘转角,然后通过转向传动机构的参数来计算获得;横摆角速度可通过陀螺仪测算获得;车身纵向速度和质心侧偏角可通过状态参数估计器结合方向盘转角、横摆角速度、车身各方向加速度计算得到,车身各方向加速度可由陀螺仪得到;
步骤二,利用汽车的参数,建立包含参数不确定性的汽车状态反馈控制模型;
在该步骤中,首先建立汽车二自由度模型,根据牛顿力学定律,可得到以下方程:
Figure BDA0002293493230000031
在上述方程中:
Fyf=Cfαf,Fyr=Crαr
Figure BDA0002293493230000032
其中,Cf、Cr分别为前后轮胎的侧偏刚度,αfr分别为前、后轮侧偏角;
选取质心侧偏角β和横摆角速度r作为控制模型的状态变量,可得到如下汽车横向运动模型:
Figure BDA0002293493230000033
其中:
x=[β r]T,u(t)=ΔMz,w(t)=δ,
Figure BDA0002293493230000034
考虑上述模型中的车辆纵向速度vx为不确定参数,令q1=1/vx,q2=1/vx 2,则系统可以描述为:
Figure BDA0002293493230000041
其中:
Figure BDA0002293493230000042
假设汽车纵向速度在一定在范围内变化,则qi在范围
Figure BDA0002293493230000043
内变化,i=1,2;则变量q的值就会在一个多胞体范围内变化,可以记为:
Figure BDA0002293493230000044
其中:
Figure BDA0002293493230000045
Figure BDA0002293493230000046
上式中,参数上部的横线表示可取的最大值,下部的横线表示可取的最小值;因为考虑纵向速度为不确定因素,因此假设纵向速度在一定范围内取值,如q 1表示q1可取的最小值即纵向速度最大时,
Figure BDA0002293493230000047
表示q1可取的最大值即纵向速度最小时;
然后,考虑上述模型中的轮胎侧偏刚度为不确定参数,即有:
Cf=Cf0+ΔCf,Cr=Cr0+ΔCr
其中Cf0,Cr0为基础值,ΔCf,ΔCr为不确定部分,且有
Figure BDA0002293493230000048
ρfr为时变的系数,且取值有|ρf|≤1,|ρr|≤1,
Figure BDA0002293493230000049
为不确定量的最大值;
基于以上分析,可得系统方程为:
Figure BDA0002293493230000051
其中:
Figure BDA0002293493230000052
Figure BDA0002293493230000053
ΔA,ΔB1可以表示为以下形式:
ΔA=HΓEa(q),ΔB1=HΓEb(q),H=[I I],Γ=diag{ρfr}.
Ea(q)=[Eaf(q) Ear(q)]T,Eb(q)=[Ebf(q) 01×2]T
Figure BDA0002293493230000054
Figure BDA0002293493230000055
也可写为如下形式:
Figure BDA0002293493230000056
为了获得更好的操纵性和稳定性,质心侧偏角应趋近于0,横摆角速度应趋近于以下参考值:
Figure BDA0002293493230000057
其中Ku为车身相关的一个固有参数;
选取系统输出变量为
Figure BDA0002293493230000061
其中cβ,cr分别为两个被控变量的权重系数;C1=diag{cβ,cr},D1(q)=[0,-crq1/(lf+lr)Ku]T
设系统的状态反馈控制器为usf=Ksf·x,则闭环系统可表示为:
Figure BDA0002293493230000062
其中Asf=A+B2Ksf
Figure BDA0002293493230000063
步骤三,求解状态反馈控制器增益矩阵Ksf
所述步骤二中建立了汽车控制系统模型,针对该系统模型,需要确定状态反馈增益矩阵Ksf的取值;选取H参数来表征输出z(t),设
Figure BDA0002293493230000064
||T1||表示系统的H增益,||z||2,||δ||2分别表示z和δ的二范数,由于实际情况中前轮转向角δ是有一定限度的,因此δ∈L2,L2表示前轮转向角的取值范围;此外,考虑控制力执行机构的饱和,控制力应满足条件||u||2≤umax
基于以上分析,通过求解以下两个线性矩阵不等式,即可解出Ksf
Figure BDA0002293493230000065
Figure BDA0002293493230000066
其中,umax和γ为给定的正标量,ε,η为需要求解的正标量,W和X为需要求解的矩阵变量,且X为正定对称矩阵,i=1,2,3,4.,通过求解出W和X,可得Ksf=WX-1
步骤四,建立汽车输出反馈的控制模型,并求解输出反馈控制器;
类似步骤二,考虑静态输出反馈控制器u=Ksof·y=KsofCyx(t),C_y=diag{0,1},则闭环系统可表达为:
Figure BDA0002293493230000071
其中Asof=A(q)+B2KsofCy.,其余参数与步骤二中相同;
然后求解以下线性矩阵不等式:
Figure BDA0002293493230000072
其中:
Figure BDA0002293493230000073
ζi=[K 0 -I]T[LiCy 0 -F],
Figure BDA0002293493230000074
γ为给定的正标量,K=Ksf,Ksf为步骤三中求解出的状态反馈控制增益矩阵,λ为需求解的正标量,Li,F,P为需求解的适当维度的矩阵变量,且P为正定对称矩阵,i=1,2,3,4.最终可得输出反馈控制增益为:
Figure BDA0002293493230000075
然后通过下式计算汽车的横摆力矩ΔMz
u=ΔMz=KsofCyx(t);
通过力矩分配系统计算出需要分配到每个车轮的力矩,再通过动力系统动作来产生所需力矩,以实现汽车横向稳定性控制。
本发明与现有技术相比具有以下的有益效果:
本发明采用启发式算法来求解输出反馈控制器,避免了传统求解方法中的复杂迭代;考虑汽车纵向速度和轮胎侧偏刚度为不确定参数,本发明所设计的控制器鲁棒性更好,能够计算出合适的横摆力矩,控制质心侧偏角和横摆角速度在一个合适的范围内,提高汽车的操纵性能,改善汽车的横向稳定性能,降低汽车运行过程中发生危险情况的可能;考虑存在难以在线测量的状态变量,本发明采用静态输出反馈的方法来设计控制器。
附图说明
图1为本发明线性二自由度汽车模型示意图;
图2为本发明控制系统示意图。
具体实施方式
下面结合实施例及附图对本发明作进一步详细的描述,但本发明的实施方式不限于此。
汽车行驶过程中,为了保证安全性和舒适性,需要控制汽车的质心侧偏角和横摆角速度在一定的范围内。当汽车高速运行的时候,要求汽车的响应速度要快、稳定性要好、安全性要高,当汽车低速运行的时候,要求汽车在保持横向稳定性的同时还能提高电动汽车的机动性能,因此要求汽车质心侧偏角的响应值应该尽量趋近于零。而横摆角速度是指汽车绕垂直轴的偏转,该偏转的大小代表汽车的稳定程度。如果偏转角速度达到一个阈值,说明汽车发生测滑或甩尾等危险工况的可能性较大。因此汽车质心侧偏角和横摆角速度这两个参数在汽车运行过程中是需要加以关注的。考虑到车辆运行中一些参数难以实时测量,本发明提出了一种基于输出反馈的横摆力矩控制方法,并利用一种新型的启发式算法来求解输出反馈控制器。考虑汽车速度和轮胎侧偏刚度为不确定变量,通过设计的控制器计算出需要提供的合适的横摆力矩,控制上述两个参数(质心侧偏角和横摆角速度)在一个合适的范围内,提高汽车的操纵性能,改善汽车的横向稳定性能,降低汽车运行过程中发生危险情况的可能。在本发明中,需要先求解系统的状态反馈控制增益,然后用状态反馈控制增益来求解输出反馈控制增益。
具体来说,如图1~2所示,一种基于启发式算法的车辆横向稳定性控制方法,包括下述步骤:
步骤一,采集汽车的参数;
需要采集的参数包括固有参数和实时参数,其中:
固有参数:包括汽车总质量m,转动惯量Iz,前、后轴到汽车重心的距离lf、lr,以及前、后轮的侧偏刚度Cf、Cr;这些参数是汽车固有的,可提前采集并存储;
实时参数:包括前轮转向角δ,汽车纵向速度vx,横摆角速度r,质心侧偏角β;这些参数中,前轮转向角可通过传感器采集方向盘转角,然后通过转向传动机构的参数来计算获得;横摆角速度可通过陀螺仪测算获得;车身纵向速度和质心侧偏角可通过状态参数估计器结合方向盘转角、横摆角速度、车身各方向加速度计算得到,车身各方向加速度可由陀螺仪得到;
步骤二,利用汽车的参数,建立包含参数不确定性的汽车状态反馈控制模型;
在该步骤中,首先建立汽车二自由度模型,该模型的示意图如图1所示,CG表示汽车的重心,m为汽车总质量,β为汽车质心侧偏角,r为横摆角速度,
Figure BDA0002293493230000101
为横摆角加速度(横摆角速度的一阶导数),Fyf和Fyr分别为前、后轮的侧向力,vy和vx分别是汽车横向和纵向速度,
Figure BDA0002293493230000102
为汽车横向加速度,Iz为车身转动惯量,ΔMz是外部提供的横摆力矩,可以通过转向系统和制动系统提供,lf、lr分别是前后轴到汽车重心的距离。根据牛顿力学定律,可得到以下方程:
Figure BDA0002293493230000103
Figure BDA0002293493230000104
在上述方程中:
Fyf=Cfαf,Fyr=Crαr
Figure BDA0002293493230000105
其中,Cf、Cr分别为前后轮胎的侧偏刚度,αfr分别为前、后轮侧偏角;
选取质心侧偏角β和横摆角速度r作为控制模型的状态变量,可得到如下汽车横向运动模型:
Figure BDA0002293493230000106
其中:
x=[βr]T,u(t)=ΔMz,w(t)=δ,
Figure BDA0002293493230000107
考虑上述模型中的车辆纵向速度vx为不确定参数,令q1=1/vx,q2=1/vx 2,则系统可以描述为:
Figure BDA0002293493230000111
其中:
Figure BDA0002293493230000112
假设汽车纵向速度在一定在范围内变化,则qi在范围
Figure BDA0002293493230000118
内变化,i=1,2;则变量q的值就会在一个多胞体范围内变化,可以记为:
Figure BDA0002293493230000113
其中:
Figure BDA0002293493230000114
Figure BDA0002293493230000115
上式中,参数上部的横线表示可取的最大值,下部的横线表示可取的最小值;因为考虑纵向速度为不确定因素,因此假设纵向速度在一定范围内取值,如q 1表示q1可取的最小值即纵向速度最大时,
Figure BDA0002293493230000119
表示q1可取的最大值即纵向速度最小时;
然后,考虑上述模型中的轮胎侧偏刚度为不确定参数,即有:
Cf=Cf0+ΔCf,Cr=Cr0+ΔCr
其中Cf0,Cr0为基础值,ΔCf,ΔCr为不确定部分,且有
Figure BDA0002293493230000116
ρfr为时变的系数,且取值有|ρf|≤1,|ρr|≤1,
Figure BDA0002293493230000117
为不确定量的最大值;
基于以上分析,可得系统方程为:
Figure BDA0002293493230000121
其中:
Figure BDA0002293493230000122
Figure BDA0002293493230000123
ΔA,ΔB1可以表示为以下形式:
ΔA=HΓEa(q),ΔB1=HΓEb(q),H=[I I],Γ=diag{ρfr}.
Ea(q)=[Eaf(q) Ear(q)]T,Eb(q)=[Ebf(q) 01×2]T
Figure BDA0002293493230000124
Figure BDA0002293493230000125
也可写为如下形式:
Figure BDA0002293493230000126
为了获得更好的操纵性和稳定性,质心侧偏角应趋近于0,横摆角速度应趋近于以下参考值:
Figure BDA0002293493230000127
其中Ku为车身相关的一个固有参数;
选取系统输出变量为
Figure BDA0002293493230000131
其中cβ,cr分别为两个被控变量的权重系数;C1=diag{cβ,cr},D1(q)=[0,-crq1/(lf+lr)Ku]T
设系统的状态反馈控制器为usf=Ksf·x,则闭环系统可表示为:
Figure BDA0002293493230000132
其中Asf=A+B2Ksf
Figure BDA0002293493230000133
步骤三,求解状态反馈控制器增益矩阵Ksf
所述步骤二中建立了汽车控制系统模型,针对该系统模型,需要确定状态反馈增益矩阵Ksf的取值;选取H参数来表征输出z(t),设
Figure BDA0002293493230000134
||T1||表示系统的H增益,||z||2,||δ||2分别表示z和δ的二范数,由于实际情况中前轮转向角δ是有一定限度的,因此δ∈L2,L2表示前轮转向角的取值范围;此外,考虑控制力执行机构的饱和,控制力应满足条件||u||2≤umax
基于以上分析,通过求解以下两个线性矩阵不等式,即可解出Ksf
Figure BDA0002293493230000135
Figure BDA0002293493230000136
其中,umax和γ为给定的正标量,ε,η为需要求解的正标量,W和X为需要求解的矩阵变量,且X为正定对称矩阵,i=1,2,3,4.,通过求解出W和X,可得Ksf=WX-1
步骤四,建立汽车输出反馈的控制模型,并求解输出反馈控制器;
类似步骤二,考虑静态输出反馈控制器u=Ksof·y=KsofCyx(t),C_y=diag{0,1},则闭环系统可表达为:
Figure BDA0002293493230000141
其中Asof=A(q)+B2KsofCy.,其余参数与步骤二中相同;
然后求解以下线性矩阵不等式:
Figure BDA0002293493230000142
其中:
Figure BDA0002293493230000143
ζi=[K 0 -I]T[LiCy 0 -F],
Figure BDA0002293493230000144
γ为给定的正标量,K=Ksf,Ksf为步骤三中求解出的状态反馈控制增益矩阵,λ为需求解的正标量,Li,F,P为需求解的适当维度的矩阵变量,且P为正定对称矩阵,i=1,2,3,4.最终可得输出反馈控制增益为:
Figure BDA0002293493230000145
然后通过下式计算汽车的横摆力矩ΔMz
u=ΔMz=KsofCyx(t);
通过力矩分配系统计算出需要分配到每个车轮的力矩,再通过动力系统动作来产生所需力矩,以实现汽车横向稳定性控制,如图2所示。
需要注意的是,在实际应用中,控制器增益Ksof需要离线计算完成,而需要实时计算的参数有前轮转向角δ,纵向速度vx和横摆角速度r。
本发明采用启发式算法来求解输出反馈控制器,避免了传统求解方法中的复杂迭代;考虑汽车纵向速度和轮胎侧偏刚度为不确定参数,本发明所设计的控制器鲁棒性更好,能够计算出合适的横摆力矩,控制质心侧偏角和横摆角速度在一个合适的范围内,提高汽车的操纵性能,改善汽车的横向稳定性能,降低汽车运行过程中发生危险情况的可能;考虑存在难以在线测量的状态变量,本发明采用静态输出反馈的方法来设计控制器。
上述为本发明较佳的实施方式,但本发明的实施方式并不受上述内容的限制,其他的任何未背离本发明的精神实质与原理下所作的改变、修饰、替代、组合、简化,均应为等效的置换方式,都包含在本发明的保护范围之内。

Claims (1)

1.一种基于启发式算法的车辆横向稳定性控制方法,其特征在于,包括下述步骤:
步骤一,采集汽车的参数;
需要采集的参数包括固有参数和实时参数,其中:
固有参数:包括汽车总质量m,转动惯量Iz,前、后轴到汽车重心的距离lf、lr,以及前、后轮的侧偏刚度Cf、Cr;这些参数是汽车固有的,可提前采集并存储;
实时参数:包括前轮转向角δ,汽车纵向速度vx,横摆角速度r,质心侧偏角β;这些参数中,前轮转向角可通过传感器采集方向盘转角,然后通过转向传动机构的参数来计算获得;横摆角速度可通过陀螺仪测算获得;车身纵向速度和质心侧偏角可通过状态参数估计器结合方向盘转角、横摆角速度、车身各方向加速度计算得到,车身各方向加速度可由陀螺仪得到;
步骤二,利用汽车的参数,建立包含参数不确定性的汽车状态反馈控制模型;
在该步骤中,首先建立汽车二自由度模型,根据牛顿力学定律,可得到以下方程:
Figure FDA0002988946430000011
Figure FDA0002988946430000012
在上述方程中:
Fyf=Cfαf,Fyr=Crαr
Figure FDA0002988946430000021
其中,Cf、Cr分别为前后轮胎的侧偏刚度,αfr分别为前、后轮侧偏角;
选取质心侧偏角β和横摆角速度r作为控制模型的状态变量,可得到如下汽车横向运动模型:
Figure FDA0002988946430000022
其中:
x(t)=[βr]T,u(t)=ΔMz,w(t)=δ,
Figure FDA0002988946430000023
考虑上述模型中的车辆纵向速度vx为不确定参数,令q1=1/vx,q2=1/vx 2,则系统可以描述为:
Figure FDA0002988946430000024
其中:
Figure FDA0002988946430000025
假设汽车纵向速度在一定在范围内变化,则qi在范围
Figure FDA0002988946430000037
内变化,i=1,2;则变量q的值就会在一个多胞体范围内变化,可以记为:
Figure FDA0002988946430000031
其中:
Figure FDA0002988946430000032
Figure FDA0002988946430000033
上式中,参数上部的横线表示可取的最大值,下部的横线表示可取的最小值;因为考虑纵向速度为不确定因素,因此假设纵向速度在一定范围内取值,如q 1表示q1可取的最小值即纵向速度最大时,
Figure FDA0002988946430000036
表示q1可取的最大值即纵向速度最小时;
然后,考虑上述模型中的轮胎侧偏刚度为不确定参数,即有:
Cf=Cf0+ΔCf,Cr=Cr0+ΔCr
其中Cf0,Cr0为基础值,ΔCf,ΔCr为不确定部分,且有
Figure FDA0002988946430000034
ρfr为时变的系数,且取值有|ρf|≤1,|ρr|≤1,
Figure FDA0002988946430000035
为不确定量的最大值;
基于以上分析,可得系统方程为:
Figure FDA0002988946430000041
其中:
Figure FDA0002988946430000042
Figure FDA0002988946430000043
ΔA,ΔB1可以表示为以下形式:
ΔA=HΓEa(q),ΔB1=HΓEb(q),H=[I I],Γ=diag{ρfr}.
Ea(q)=[Eaf(q) Ear(q)]T,Eb(q)=[Ebf(q) 01×2]T
Figure FDA0002988946430000044
Figure FDA0002988946430000045
也可写为如下形式:
Figure FDA0002988946430000046
为了获得更好的操纵性和稳定性,质心侧偏角应趋近于0,横摆角速度应趋近于以下参考值:
Figure FDA0002988946430000051
其中Ku为车身相关的一个固有参数;
选取系统输出变量为
Figure FDA0002988946430000052
其中cβ,cr分别为两个被控变量的权重系数;
C1=diag{cβ,cr},D1(q)=[0,-crq1/(lf+lr)Ku]T
设系统的状态反馈控制器为usf=Ksf·x,则闭环系统可表示为:
Figure FDA0002988946430000053
其中Asf=A+B2Ksf
Figure FDA0002988946430000054
步骤三,求解状态反馈控制器增益矩阵Ksf
所述步骤二中建立了汽车控制系统模型,针对该系统模型,需要确定状态反馈增益矩阵Ksf的取值;选取H参数来表征输出z(t),设
Figure FDA0002988946430000055
||T1||表示系统的H增益,||z||2,||δ||2分别表示z和δ的二范数,由于实际情况中前轮转向角δ是有一定限度的,因此δ∈L2,L2表示前轮转向角的取值范围;此外,考虑控制力执行机构的饱和,控制力应满足条件||u||2≤umax
基于以上分析,通过求解以下两个线性矩阵不等式,即可解出Ksf
Figure FDA0002988946430000061
Figure FDA0002988946430000062
其中,umax和γ为给定的正标量,ε,η为需要求解的正标量,W和X为需要求解的矩阵变量,且X为正定对称矩阵,i=1,2,3,4.,通过求解出W和X,可得Ksf=WX-1
步骤四,建立汽车输出反馈的控制模型,并求解输出反馈控制器;
类似步骤二,考虑静态输出反馈控制器u=Ksof·y=KsofCyx(t),C_y=diag{0,1},则闭环系统可表达为:
Figure FDA0002988946430000063
其中Asof=A(q)+B2KsofCy.,其余参数与步骤二中相同;
然后求解以下线性矩阵不等式:
Figure FDA0002988946430000064
其中:
Figure FDA0002988946430000071
ζi=[K 0 -I]T[LiCy 0 -F],
Figure FDA0002988946430000072
γ为给定的正标量,K=Ksf,Ksf为步骤三中求解出的状态反馈控制增益矩阵,λ为需求解的正标量,Li,F,P为需求解的适当维度的矩阵变量,且P为正定对称矩阵,i=1,2,3,4.最终可得输出反馈控制增益为:
Figure FDA0002988946430000073
然后通过下式计算汽车的横摆力矩ΔMz
u=ΔMz=KsofCyx(t);
通过力矩分配系统计算出需要分配到每个车轮的力矩,再通过动力系统动作来产生所需力矩,以实现汽车横向稳定性控制。
CN201911190692.XA 2019-11-28 2019-11-28 一种基于启发式算法的车辆横向稳定性控制方法 Active CN111055837B (zh)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN201911190692.XA CN111055837B (zh) 2019-11-28 2019-11-28 一种基于启发式算法的车辆横向稳定性控制方法
AU2020102015A AU2020102015A4 (en) 2019-11-28 2020-08-27 A vehicle lateral stability control method based on the heuristic algorithm

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911190692.XA CN111055837B (zh) 2019-11-28 2019-11-28 一种基于启发式算法的车辆横向稳定性控制方法

Publications (2)

Publication Number Publication Date
CN111055837A CN111055837A (zh) 2020-04-24
CN111055837B true CN111055837B (zh) 2021-05-11

Family

ID=70299408

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911190692.XA Active CN111055837B (zh) 2019-11-28 2019-11-28 一种基于启发式算法的车辆横向稳定性控制方法

Country Status (2)

Country Link
CN (1) CN111055837B (zh)
AU (1) AU2020102015A4 (zh)

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111824146A (zh) * 2020-06-19 2020-10-27 武汉理工大学 一种路径跟随模型预测控制方法、系统、装置及存储介质
CN112590483B (zh) * 2021-01-05 2023-06-16 广东工业大学 一种基于观测器的汽车横向稳定性与主动悬架多目标集成控制方法
CN113002527B (zh) * 2021-03-01 2022-09-30 东北大学 一种自主电动汽车横向稳定性鲁棒容错控制方法
CN113176733B (zh) * 2021-04-27 2023-06-16 广东工业大学 一种基于切换控制的自主车辆路径跟踪与稳定性控制方法
CN113815602B (zh) * 2021-07-19 2023-10-24 北京理工大学 一种车辆横向自适应控制方法
CN113753080B (zh) * 2021-08-31 2023-09-26 的卢技术有限公司 一种自动驾驶汽车横向运动自适应参数控制方法
CN114253128B (zh) * 2021-12-22 2023-06-30 合肥工业大学 包含反馈时滞的汽车自动转向系统的控制参数优化方法
CN114368380B (zh) * 2022-01-06 2023-02-17 上海宏景智驾信息科技有限公司 一种适应不同载重的自动驾驶半挂卡车横向控制方法
CN117360486B (zh) * 2023-12-08 2024-03-08 福州大学 一种用于多轴线控底盘的抗干扰直接横摆力矩控制方法
CN117369287B (zh) * 2023-12-08 2024-03-12 福州大学 一种用于多轴车辆的鲁棒超螺旋滑模控制方法

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE19650691C2 (de) * 1996-12-07 1998-10-29 Deutsch Zentr Luft & Raumfahrt Verfahren zur Lenkunterstützung eines Fahrers eines Straßenfahrzeugs
US7835846B2 (en) * 2007-02-28 2010-11-16 Gm Global Technology Operations, Inc. Nonlinear vehicle yaw/roll/sideslip command interpreter
US8725310B2 (en) * 2009-04-29 2014-05-13 Koninklijke Philips N.V. Laser diode based self-mixing sensor for a vehicle electronic stability program
CN105172790B (zh) * 2015-10-30 2017-09-12 吉林大学 一种基于三步法的车辆横摆稳定性控制方法
CN109367532B (zh) * 2018-08-20 2020-04-17 广东工业大学 一种基于速度依赖的汽车横向稳定性控制方法
CN109849898B (zh) * 2018-12-27 2020-06-26 合肥工业大学 基于遗传算法混合优化gpc的车辆横摆稳定性控制方法

Also Published As

Publication number Publication date
CN111055837A (zh) 2020-04-24
AU2020102015A4 (en) 2020-10-01

Similar Documents

Publication Publication Date Title
CN111055837B (zh) 一种基于启发式算法的车辆横向稳定性控制方法
CN109522666B (zh) 一种分布式电动汽车稳定性控制方法
CN104773170B (zh) 一种车辆稳定性集成控制方法
CN106080553B (zh) 一种融合车速变化的四轮转向汽车防侧翻控制系统及方法
CN109367532B (zh) 一种基于速度依赖的汽车横向稳定性控制方法
CN107215329B (zh) 一种基于atsm的分布式驱动电动汽车横向稳定性控制方法
CN103303367B (zh) 一种四轮驱动电动车车身稳定控制方法
CN110606079B (zh) 一种分层控制的车辆防侧翻方法及多轴分布式驱动车辆
CN109606133A (zh) 基于双层控制的分布式驱动电动汽车转矩矢量控制方法
CN106394561A (zh) 一种车辆的纵向车速的估计方法和装置
CN110509915B (zh) 一种基于时变速度的四轮驱动汽车横向稳定性控制方法
EP4253181A1 (en) Vehicle front and rear drive torque distribution method and apparatus, and vehicle
CN106515716A (zh) 轮式驱动电动汽车底盘集成控制系统协调控制装置及方法
CN105936273A (zh) 车用主动转矩轮间、轴间分配方法
CN109291932A (zh) 基于反馈的电动汽车横摆稳定性实时控制装置及方法
CN110161854B (zh) 一种高速公路重型卡车编队纵向行驶控制方法
IT202100015182A1 (it) Metodo di controllo degli ammortizzatori attivi di un veicolo stradale che prevede la regolazione dell'angolo di rollio e dell'angolo di beccheggio
CN111891118A (zh) 基于模型预测控制算法的电动汽车防侧翻控制方法
IT202100015170A1 (it) Metodo di controllo degli ammortizzatori attivi di un veicolo stradale che prevede l'abbassamento del baricentro
CN111731267B (zh) 一种装备非充气弹性车轮的分布式电动汽车稳定性控制系统及方法
KR20150062779A (ko) 인휠 전기자동차의 휠 토크 제어 시스템 및 방법
CN111002976A (zh) 一种基于模糊自适应pid控制的智能车辆抗侧风控制方法
CN109850015B (zh) 一种控制参数可自动调节的电动车主动前轮转向控制方法
CN112590483B (zh) 一种基于观测器的汽车横向稳定性与主动悬架多目标集成控制方法
CN113928306B (zh) 汽车集成增稳控制方法及系统

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