CN109050659A - 一种基于时变动力学模型的四轮转向汽车稳定性控制方法 - Google Patents
一种基于时变动力学模型的四轮转向汽车稳定性控制方法 Download PDFInfo
- Publication number
- CN109050659A CN109050659A CN201810735047.0A CN201810735047A CN109050659A CN 109050659 A CN109050659 A CN 109050659A CN 201810735047 A CN201810735047 A CN 201810735047A CN 109050659 A CN109050659 A CN 109050659A
- Authority
- CN
- China
- Prior art keywords
- tire
- automobile
- model
- slip angle
- yaw velocity
- 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.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 49
- 230000033001 locomotion Effects 0.000 claims abstract description 14
- 238000005457 optimization Methods 0.000 claims abstract description 10
- 238000005096 rolling process Methods 0.000 claims description 8
- 238000004088 simulation Methods 0.000 claims description 6
- 230000005540 biological transmission Effects 0.000 claims description 4
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 241000252147 Atractosteus spatula Species 0.000 claims description 3
- 238000013016 damping Methods 0.000 claims description 3
- 238000010586 diagram Methods 0.000 description 5
- 230000000694 effects Effects 0.000 description 5
- 238000004458 analytical method Methods 0.000 description 2
- 238000005312 nonlinear dynamic Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B62—LAND VEHICLES FOR TRAVELLING OTHERWISE THAN ON RAILS
- B62D—MOTOR VEHICLES; TRAILERS
- B62D6/00—Arrangements for automatically controlling steering depending on driving conditions sensed and responded to, e.g. control circuits
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/10—Internal combustion engine [ICE] based vehicles
- Y02T10/40—Engine management systems
Landscapes
- Engineering & Computer Science (AREA)
- Transportation (AREA)
- Mechanical Engineering (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Automation & Control Theory (AREA)
- Steering Control In Accordance With Driving Conditions (AREA)
- Control Of Driving Devices And Active Controlling Of Vehicle (AREA)
Abstract
一种基于时变动力学模型的四轮转向汽车稳定性控制方法,其特征在于,该方法包括参考模型、轮胎数据处理器、MPC控制器、CarSim汽车模型;参考模型用于确定期望的汽车横摆角速度和质心侧偏角;轮胎数据处理器用于确定轮胎的侧偏角、侧向力和侧向力梯度;CarSim汽车模型用于输出汽车的实际运动状态信息,包括汽车纵向速度、横摆角速度、质心侧偏角和路面附着系数;MPC控制器根据期望的汽车横摆角速度、质心侧偏角和汽车的实际运动状态信息,优化求解出汽车的后轮转角,输出给CarSim汽车模型,控制汽车实现稳定性控制。
Description
技术领域:
本发明涉及汽车稳定性控制领域,特别是关于一种基于时变动力学模型的四轮转向汽车稳定性控制方法。
背景技术:
随着人们对汽车行驶安全性越来越重视,汽车主动安全***得到快速发展,其中四轮转向(Four Wheel Steering,4WS)作为一种提高汽车稳定性的有效手段已得到了广泛的认可,运用4WS技术,可以有效地减小低速行驶时汽车的转弯半径,提高高速行驶时汽车的稳定性。
目前,4WS所采用的控制方法主要有PID控制、鲁棒性控制、神经网络控制和模型预测控制(Model Predictive Control,MPC)等,其中模型预测控制因能较好地处理多目标任务以及***约束,在汽车稳定性控制领域得到了广泛的应用。
根据其采用的预测模型以及优化方法的不同,MPC可分为非线性MPC和线性MPC。非线性MPC的计算负担重,实时性差,很难应用于实际,而线性MPC计算负担小,计算速度快,因而得到了广泛使用。论文[Hang P,Luo F,Fang S,et al.Path Tracking Control of aFour-Wheel-Independent-Steering Electric Vehicle based on Model PredictiveControl[C]//中国控制会议.2017:9360-9366.]假设轮胎侧偏角处于线性区域,将非线性的轮胎模型近似成线性轮胎模型,通过线性MPC实现了四轮转向控制,然而线性MPC不能表征轮胎的非线性特性侧偏特性,因此在极限工况下控制器的控制效果差。论文[陈杰,李亮,宋健.基于LTV-MPC的汽车稳定性控制研究[J].汽车工程,2016,38(3):308-316.]基于线性化的轮胎模型采用线性时变的MPC方法实现了汽车稳定性控制,同时兼顾了***的非线性特性和计算负担。但是,论文中对轮胎侧向力线性化处理的方法过于简单,无法表征轮胎侧向力的实际变化趋势,在极限工况下控制器的控制效果不理想;此外,该方法采用的预测模型在预测时域里保持不变,在滚动预测过程中不能体现汽车实际的变化趋势。
发明内容:
为了解决现有的线性时变MPC预测模型不能体现汽车的非线性动力学特性,导致极限工况下4WS***控制效果差的问题。本发明提供一种基于时变动力学模型的四轮转向汽车稳定性控制方法,采用线性时变的方法将非线性预测控制问题转换成线性预测控制问题,并且在预测时域内根据轮胎侧向力的变化趋势自动调节预测模型,在减小***的计算负担的同时准确表征汽车的非线性动力学特性,保证极限工况下4WS控制器的稳定性,实现汽车的稳定性控制。
本发明解决技术问题所采取的技术方案如下:
一种基于时变动力学模型的四轮转向汽车稳定性控制方法,其特征在于,该方法包括参考模型、轮胎数据处理器、MPC控制器、CarSim汽车模型;参考模型用于确定期望的汽车横摆角速度和质心侧偏角;轮胎数据处理器用于确定轮胎的侧偏角、侧向力和侧向力梯度;CarSim汽车模型用于输出汽车的实际运动状态信息,包括汽车纵向速度、横摆角速度、质心侧偏角和路面附着系数;MPC控制器根据期望的汽车横摆角速度、质心侧偏角和汽车的实际运动状态信息,优化求解出汽车的后轮转角,输出给CarSim汽车模型,控制汽车实现稳定性控制;
该方法包括以下步骤:
步骤1、建立参考模型,确定期望的汽车横摆角速度和质心侧偏角,其过程包括如下子步骤:
步骤1.1、采用线性二自由度汽车模型作为参考模型,其运动微分方程表达式如下:
其中:β是汽车质心侧偏角;γ是汽车横摆角速度;Iz是绕汽车质心铅垂轴的横摆转动惯量;Ux是汽车纵向速度;lf和lr分别是汽车质心至前、后轴的距离;Cf和Cr分别是汽车前、后轮轮胎的侧偏刚度;δf是驾驶员操纵方向盘产生的汽车前轮转角;
步骤1.2、将公式(1a)转换成传递函数,形式如下式:
基于公式(2)得到期望的汽车横摆角速度:
其中:γref是期望的汽车横摆角速度;wn是***的固有频率;ξ是***阻尼;Gω(s)是传递函数增益;wd=k1wn,ξd=k2ξ,Gkω(s)=k3Gω(s);k1、k2、k3是改善***相位延迟和响应速度的参数;
步骤1.3、设定期望的质心侧偏角为:
βref=0 (4)
步骤2、设计轮胎数据处理器,其过程包括如下子步骤:
步骤2.1、设计轮胎侧偏角计算模块,前、后轮轮胎侧偏角通过下式计算获得:
其中:αf和αr分别是汽车前、后轮轮胎的侧偏角;δf是汽车前轮转角,δr是汽车后轮转角;
步骤2.2、设计轮胎侧向力和轮胎侧向力梯度计算模块,为了获得轮胎的非线性特性,基于Pacejka轮胎模型,获取不同路面附着系数下的轮胎侧向力与轮胎侧偏角的关系曲线,得到轮胎侧偏特性三维图;获取不同路面附着系数下的轮胎侧向力对轮胎侧偏角导数的关系曲线,得到轮胎侧向力梯度三维图;轮胎数据处理器将当前时刻实际的轮胎侧偏角和路面附着系数分别输入到轮胎侧偏特性三维图和轮胎侧向力梯度三维图中,通过线性插值法分别获得当前时刻的轮胎侧向力和轮胎侧向力梯度,并输出给MPC控制器;在每个控制周期轮胎数据处理器更新一次轮胎侧向力和轮胎侧向力梯度值;
其中:Pacejka轮胎模型如下:
Fy,j=μDsin(Catan(Bαj-E(Bαj-αjtan(Bαj))))
其中:j=f,r,表示前轮和后轮;Fy,j是轮胎侧向力,αj是轮胎侧偏角;B,C,D和E取决于车轮垂直载荷Fz;a0=1.75;a1=0;a2=1000;a3=1289;a4=7.11;a5=0.0053;a6=0.1925;
步骤3、设计MPC控制器,其过程包括如下子步骤:
步骤3.1、建立预测模型,其过程包括如下子步骤:
步骤3.1.1、线性化轮胎模型,其表达式如下:
其中:是在当前侧偏角的轮胎侧向力梯度值;是轮胎的残余侧向力,通过如下公式计算:
其中:是基于轮胎侧偏特性三维图,通过线性插值法获得的轮胎侧向力;是基于轮胎侧偏刚度特性三维图,通过线性插值法获得的轮胎侧向力梯度;是当前时刻实际的轮胎侧偏角;
基于公式(6),在滚动预测过程中,设计轮胎侧向力表达式如下:
其中:
其中:P是预测时域;上标“k+i|k”表示在当前时刻k预测的将来第i时刻;ρk+i|k和ξk+i|k是调节和变化的权重因子;
步骤3.1.2、建立预测模型,其运动微分方程表达式为:
将公式(9)带入公式(10),得到在滚动预测过程中的预测模型为:
步骤3.1.3、建立预测方程,用于预测***未来输出,将公式(11)写成状态空间方程,用于设计预测方程,具体如下:
其中:
为了实现汽车横摆角速度和质心侧偏角的跟踪控制,将连续时间***的预测模型转换成离散时间***的增量式模型:
其中:取样时间k=int(t/Ts),t是仿真时间,Ts是仿真步长;
步骤3.2、设计优化目标及约束条件,其过程包括如下子步骤:
步骤3.2.1、用期望的汽车横摆角速度、质心侧偏角和实际的汽车横摆角速度、质心侧偏角误差的二范数作为横摆角速度、质心侧偏角跟踪性能指标,体现汽车的轨迹跟踪特性,其表达式如下:
其中:γref是期望的汽车横摆角速度;γ是实际的汽车横摆角速度;βref是期望的汽车质心侧偏角;β是实际的汽车质心侧偏角;P是预测时域;k表示当前时刻;Q1、Q2是加权因子;
步骤3.2.2、用控制量变化率的二范数作为转向平滑指标,体现横摆角速度和质心侧偏角跟踪过程中的转向平滑特性,控制量u是汽车后轮转角,建立离散二次型转向平滑指标为:
其中:M是控制时域;△u是控制量的变化量;k表示当前时刻;S是加权因子;
步骤3.2.3、设置执行器物理约束,满足执行器要求:
利用线性不等式限制后轮转角及其变化量的上下限,得到转向执行器的物理约束,其数学表达式为:
其中:δrmin是后轮转角下限,δrmax是后轮转角上限;△δrmin是后轮转角变化量的下限;△δrmax是后轮转角变化量的上限;
步骤3.3、求解***预测输出,其过程包括如下子步骤:
步骤3.3.1、利用线性加权法将步骤3.2.1所述跟踪性能指标和步骤3.2.2所述转向平滑指标转化为单一指标,构建汽车稳定性多目标优化控制问题,该问题要满足转向执行器的物理约束,且输入输出符合预测模型:
服从于
i)预测模型
ii)约束条件为公式(16)
步骤3.3.2、在控制器中,采用二次规划算法,求解多目标优化控制问题(17),得到最优开环控制序列△δr为:
选取当前时刻最优开环控制序列中的第一个元素△δr(0)进行反馈,与前一时刻进行线性叠加,输出给CarSim汽车模型,实现汽车的稳定性控制。
本发明的有益效果是:本方法使用线性时变的方法将非线性预测控制问题转换成线性预测控制问题,可以减小***的计算负担;本方法根据轮胎侧向力的变化趋势在滚动预测过程中自适应调节***的预测模型,能够达到非线性MPC的控制效果,提升极限工况下4WS的控制效果。
附图说明
图1是本发明的控制***结构示意图。
图2是线性二自由度汽车模型示意图。
图3是轮胎侧偏特性三维图。
图4是轮胎侧向力梯度三维图。
图5是轮胎模型线性化示意图。
图6是在滚动预测过程中轮胎模型线性化示意图。
图7是四轮转向汽车模型示意图。
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
图1是本发明一种基于时变动力学模型的四轮转向汽车稳定性控制方法的***结构示意图,该***主要包括参考模型1、轮胎数据处理器2、MPC控制器3、CarSim汽车模型4;参考模型1用于确定期望的汽车横摆角速度和质心侧偏角;轮胎数据处理器2用于确定轮胎的侧偏角、侧向力和侧向力梯度;CarSim汽车模型4用于输出汽车的实际运动状态信息,包括汽车纵向速度、横摆角速度、质心侧偏角和路面附着系数;MPC控制器3根据期望的汽车横摆角速度、质心侧偏角和汽车的实际运动状态信息,优化求解出汽车的后轮转角,输出给CarSim汽车模型4,控制汽车实现稳定性控制。
下面以CarSim汽车仿真软件某车型为平台,具体说明本发明的方法,其主要参数如表1所示:
表1 CarSim汽车的主要参数
参考模型1的建立包括三部分:1.1建立线性二自由度汽车模型;1.2确定期望的汽车横摆角速度;1.3确定期望的汽车质心侧偏角;
在1.1部分中,线性二自由度汽车模型如图2所示,其运动微分方程表达式如下:
其中:β是汽车质心侧偏角;γ是汽车横摆角速度;Iz是绕汽车质心铅垂轴的横摆转动惯量;Ux是汽车纵向速度;lf和lr分别是汽车质心至前、后轴的距离;Cf和Cr分别是汽车前、后轮轮胎的侧偏刚度;δf是驾驶员操纵方向盘产生的汽车前轮转角。
在1.2部分中,将公式(1a)转换成传递函数,形式如下式:
基于公式(2)得到期望的汽车横摆角速度:
其中:γref是期望的横摆角速度;wn是***的固有频率;ξ是***阻尼;Gω(s)是传递函数增益;wd=k1wn,ξd=k2ξ,Gkω(s)=k3Gω(s);k1、k2、k3是改善***相位延迟和响应速度的参数;wn、ξ、Gω(s)、Kω的计算过程如下:
在1.3部分中,设定期望的质心侧偏角为:
βref=0 (4)
轮胎数据处理器2的设计包括两部分:2.1设计轮胎侧偏角计算模块;2.2设计轮胎侧向力和轮胎侧向力梯度计算模块;
在2.1部分中,前、后轮轮胎侧偏角通过下式计算获得:
其中:αf和αr分别是汽车前、后轮轮胎的侧偏角;δf是驾驶员操纵方向盘产生的汽车前轮转角,δr是控制器优化求解出的汽车后轮转角。
在2.2部分中,为了获得轮胎的非线性特性,基于Pacejka轮胎模型,获取不同路面附着系数下的轮胎侧向力与轮胎侧偏角的关系曲线,得到轮胎侧偏特性三维图,如图3;获取不同路面附着系数下的轮胎侧向力对轮胎侧偏角导数的关系曲线,得到轮胎侧向力梯度三维图,如图4。轮胎数据处理器2将当前时刻实际的轮胎侧偏角和路面附着系数分别输入到轮胎侧偏特性三维图和轮胎侧向力梯度三维图中,通过线性插值法分别获得当前时刻的轮胎侧向力和轮胎侧向力梯度,并输出给MPC控制器3。在每个控制周期轮胎数据处理器2更新一次轮胎侧向力和轮胎侧向力梯度值。
其中:Pacejka轮胎模型如下:
Fy,j=μDsin(Catan(Bαj-E(Bαj-αjtan(Bαj))))
其中:j=f,r,表示前轮和后轮;Fy,j是轮胎侧向力,αj是轮胎侧偏角;B,C,D和E取决于车轮垂直载荷Fz;a0=1.75;a1=0;a2=1000;a3=1289;a4=7.11;a5=0.0053;a6=0.1925。
MPC控制器3的设计包括三部分:3.1建立预测模型及预测方程;3.2设计优化目标及约束条件;3.3求解***预测输出;
在3.1部分中,预测模型及预测方程的建立包括三部分:3.1.1线性化轮胎模型;3.1.2建立预测模型;3.1.3建立预测方程;
在3.1.1部分中,在当前侧偏角处,如图5所示,对轮胎模型进行线性化,其表达式如下:
其中:是在当前侧偏角的轮胎侧向力梯度值;是轮胎的残余侧向力,通过如下公式计算:
其中:是基于轮胎侧偏特性三维图(图3),通过线性插值法获得的轮胎侧向力;是基于轮胎侧偏刚度特性三维图(图4),通过线性插值法获得的轮胎侧向力梯度;是当前时刻实际的轮胎侧偏角。
基于公式(5),在滚动预测过程中,如图6,设计轮胎侧向力表达式如下:
其中:
其中:P是预测时域;上标“k+i|k”表示在当前时刻k预测的将来第i时刻;ρk+i|k和ξk+i|k是调节和变化的权重因子。
在3.1.2部分中,预测模型采用图7所示的四轮转向汽车模型,其运动微分方程表达式为:
将公式(9)带入公式(10),得到在滚动预测过程中的预测模型为:
在3.1.3部分中,将公式(11)写成状态空间方程,用于设计预测方程,具体如下:
其中:
为了实现汽车横摆角速度和质心侧偏角的跟踪控制,将连续时间***的预测模型转换成离散时间***的增量式模型:
其中:取样时间k=int(t/Ts),t是仿真时间,Ts是仿真步长;
在3.2部分中优化目标及约束条件的设计包括三部分:3.2.1设计横摆角速度和质心侧偏角跟踪性能指标;3.2.2设计转向平滑指标;3.2.3设置执行器物理约束;
在3.2.1部分中,用期望的汽车横摆角速度、质心侧偏角和实际的汽车横摆角速度、质心侧偏角误差的二范数作为横摆角速度、质心侧偏角跟踪性能指标,体现汽车的轨迹跟踪特性,其表达式如下:
其中:γref是期望的汽车横摆角速度;γ是实际的汽车横摆角速度;βref是期望的汽车质心侧偏角;β是实际的汽车质心侧偏角;P是预测时域;k表示当前时刻;Q1、Q2是加权因子。
在3.2.2部分中,用控制量变化率的二范数作为转向平滑指标,体现横摆角速度和质心侧偏角跟踪过程中的转向平滑特性,控制量u是汽车后轮转角,建立离散二次型转向平滑指标为:
其中:M是控制时域;△u是控制量的变化量;k表示当前时刻;S是加权因子;
在3.2.3部分中,利用线性不等式限制后轮转角及其变化量的上下限,得到转向执行器的物理约束,其数学表达式为:
其中:δrmin是后轮转角下限,δrmax是后轮转角上限;△δrmin是后轮转角变化量的下限;△δrmax是后轮转角变化量的上限
在3.3部分中,***预测输出的求解包括两部分:3.3.1构建汽车稳定性多目标优化控制问题;3.3.2求解多目标优化控制问题;
在3.3.1部分中,利用线性加权法将公式(14)的横摆角速度和质心侧偏角跟踪性能指标和公式(15)的转向平滑指标转化为单一指标,构建汽车稳定性多目标优化控制问题,该问题要满足转向执行器的物理约束,且输入输出符合预测模型:
服从于
i)预测模型
ii)约束条件为公式(16)
在3.3.2部分中,在控制器中,采用二次规划算法,求解多目标优化控制问题(17),得到最优开环控制序列△δr为:
选取当前时刻最优开环控制序列中的第一个元素△δr(0)进行反馈,与前一时刻进行线性叠加,输出给CarSim汽车模型4,实现汽车的稳定性控制。
Claims (1)
1.一种基于时变动力学模型的四轮转向汽车稳定性控制方法,其特征在于,该方法包括参考模型、轮胎数据处理器、MPC控制器、CarSim汽车模型;参考模型用于确定期望的汽车横摆角速度和质心侧偏角;轮胎数据处理器用于确定轮胎的侧偏角、侧向力和侧向力梯度;CarSim汽车模型用于输出汽车的实际运动状态信息,包括汽车纵向速度、横摆角速度、质心侧偏角和路面附着系数;MPC控制器根据期望的汽车横摆角速度、质心侧偏角和汽车的实际运动状态信息,优化求解出汽车的后轮转角,输出给CarSim汽车模型,控制汽车实现稳定性控制;
该方法包括以下步骤:
步骤1、建立参考模型,确定期望的汽车横摆角速度和质心侧偏角,其过程包括如下子步骤:
步骤1.1、采用线性二自由度汽车模型作为参考模型,其运动微分方程表达式如下:
其中:β是汽车质心侧偏角;γ是汽车横摆角速度;Iz是绕汽车质心铅垂轴的横摆转动惯量;Ux是汽车纵向速度;lf和lr分别是汽车质心至前、后轴的距离;Cf和Cr分别是汽车前、后轮轮胎的侧偏刚度;δf是驾驶员操纵方向盘产生的汽车前轮转角;
步骤1.2、将公式(1a)转换成传递函数,形式如下式:
基于公式(2)得到期望的汽车横摆角速度:
其中:γref是期望的汽车横摆角速度;wn是***的固有频率;ξ是***阻尼;Gω(s)是传递函数增益;wd=k1wn,ξd=k2ξ,Gkω(s)=k3Gω(s);k1、k2、k3是改善***相位延迟和响应速度的参数;
步骤1.3、设定期望的质心侧偏角为:
βref=0 (4)
步骤2、设计轮胎数据处理器,其过程包括如下子步骤:
步骤2.1、设计轮胎侧偏角计算模块,前、后轮轮胎侧偏角通过下式计算获得:
其中:αf和αr分别是汽车前、后轮轮胎的侧偏角;δf是汽车前轮转角,δr是汽车后轮转角;
步骤2.2、设计轮胎侧向力和轮胎侧向力梯度计算模块,为了获得轮胎的非线性特性,基于Pacejka轮胎模型,获取不同路面附着系数下的轮胎侧向力与轮胎侧偏角的关系曲线,得到轮胎侧偏特性三维图;获取不同路面附着系数下的轮胎侧向力对轮胎侧偏角导数的关系曲线,得到轮胎侧向力梯度三维图;轮胎数据处理器将当前时刻实际的轮胎侧偏角和路面附着系数分别输入到轮胎侧偏特性三维图和轮胎侧向力梯度三维图中,通过线性插值法分别获得当前时刻的轮胎侧向力和轮胎侧向力梯度,并输出给MPC控制器;在每个控制周期轮胎数据处理器更新一次轮胎侧向力和轮胎侧向力梯度值;
其中:Pacejka轮胎模型如下:
Fy,j=μDsin(Catan(Bαj-E(Bαj-αjtan(Bαj))))
其中:j=f,r,表示前轮和后轮;Fy,j是轮胎侧向力,αj是轮胎侧偏角;B,C,D和E取决于车轮垂直载荷Fz;a0=1.75;a1=0;a2=1000;a3=1289;a4=7.11;a5=0.0053;a6=0.1925;
步骤3、设计MPC控制器,其过程包括如下子步骤:
步骤3.1、建立预测模型,其过程包括如下子步骤:
步骤3.1.1、线性化轮胎模型,其表达式如下:
其中:是在当前侧偏角的轮胎侧向力梯度值;是轮胎的残余侧向力,通过如下公式计算:
其中:是基于轮胎侧偏特性三维图,通过线性插值法获得的轮胎侧向力;是基于轮胎侧偏刚度特性三维图,通过线性插值法获得的轮胎侧向力梯度;是当前时刻实际的轮胎侧偏角;
基于公式(6),在滚动预测过程中,设计轮胎侧向力表达式如下:
其中:
其中:P是预测时域;上标“k+i|k”表示在当前时刻k预测的将来第i时刻;ρk+i|k和ξk+i|k是调节和变化的权重因子;
步骤3.1.2、建立预测模型,其运动微分方程表达式为:
将公式(9)带入公式(10),得到在滚动预测过程中的预测模型为:
步骤3.1.3、建立预测方程,用于预测***未来输出,将公式(11)写成状态空间方程,用于设计预测方程,具体如下:
y(k)=Cx(t) (12)
其中:
为了实现汽车横摆角速度和质心侧偏角的跟踪控制,将连续时间***的预测模型转换成离散时间***的增量式模型:
其中:取样时间k=int(t/Ts),t是仿真时间,Ts是仿真步长;
步骤3.2、设计优化目标及约束条件,其过程包括如下子步骤:
步骤3.2.1、用期望的汽车横摆角速度、质心侧偏角和实际的汽车横摆角速度、质心侧偏角误差的二范数作为横摆角速度、质心侧偏角跟踪性能指标,体现汽车的轨迹跟踪特性,其表达式如下:
其中:γref是期望的汽车横摆角速度;γ是实际的汽车横摆角速度;βref是期望的汽车质心侧偏角;β是实际的汽车质心侧偏角;P是预测时域;k表示当前时刻;Q1、Q2是加权因子;
步骤3.2.2、用控制量变化率的二范数作为转向平滑指标,体现横摆角速度和质心侧偏角跟踪过程中的转向平滑特性,控制量u是汽车后轮转角,建立离散二次型转向平滑指标为:
其中:M是控制时域;△u是控制量的变化量;k表示当前时刻;S是加权因子;
步骤3.2.3、设置执行器物理约束,满足执行器要求:
利用线性不等式限制后轮转角及其变化量的上下限,得到转向执行器的物理约束,其数学表达式为:
其中:δrmin是后轮转角下限,δrmax是后轮转角上限;△δrmin是后轮转角变化量的下限;△δrmax是后轮转角变化量的上限;
步骤3.3、求解***预测输出,其过程包括如下子步骤:
步骤3.3.1、利用线性加权法将步骤3.2.1所述跟踪性能指标和步骤3.2.2所述转向平滑指标转化为单一指标,构建汽车稳定性多目标优化控制问题,该问题要满足转向执行器的物理约束,且输入输出符合预测模型:
服从于
i)预测模型
ii)约束条件为公式(16)
步骤3.3.2、在控制器中,采用二次规划算法,求解多目标优化控制问题(17),得到最优开环控制序列△δr为:
选取当前时刻最优开环控制序列中的第一个元素△δr(0)进行反馈,与前一时刻进行线性叠加,输出给CarSim汽车模型,实现汽车的稳定性控制。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810735047.0A CN109050659B (zh) | 2018-07-06 | 2018-07-06 | 一种基于时变动力学模型的四轮转向汽车稳定性控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810735047.0A CN109050659B (zh) | 2018-07-06 | 2018-07-06 | 一种基于时变动力学模型的四轮转向汽车稳定性控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109050659A true CN109050659A (zh) | 2018-12-21 |
CN109050659B CN109050659B (zh) | 2020-09-18 |
Family
ID=64818891
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810735047.0A Expired - Fee Related CN109050659B (zh) | 2018-07-06 | 2018-07-06 | 一种基于时变动力学模型的四轮转向汽车稳定性控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109050659B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110217229A (zh) * | 2019-06-25 | 2019-09-10 | 长春工业大学 | 一种适用于高速极限工况的路径跟踪控制方法 |
CN110539752A (zh) * | 2019-06-26 | 2019-12-06 | 江苏大学 | 一种智能汽车多预测范围模型预测轨迹跟踪控制方法及*** |
CN111055921A (zh) * | 2019-12-31 | 2020-04-24 | 吉林大学 | 一种基于数据驱动的四轮转向模型预测控制方法 |
CN111452801A (zh) * | 2019-01-21 | 2020-07-28 | 上海汽车集团股份有限公司 | 一种四轮转向汽车的鲁棒自适应控制方法和装置 |
CN117125079A (zh) * | 2023-08-28 | 2023-11-28 | 广州汽车集团股份有限公司 | 车辆的控制方法及装置、电子设备、存储介质 |
Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN2435298Y (zh) * | 1999-01-07 | 2001-06-20 | 杨应捷 | 四轮汽车的四轮转向机构 |
US20010008986A1 (en) * | 1999-12-21 | 2001-07-19 | Brown Todd Allen | Roll over stability control for an automotive vehicle |
JP2007237933A (ja) * | 2006-03-08 | 2007-09-20 | Toyota Motor Corp | 車体姿勢制御装置 |
CN101400558A (zh) * | 2006-03-31 | 2009-04-01 | 丰田自动车株式会社 | 车辆运动控制装置及控制方法 |
CN101565043A (zh) * | 2008-04-25 | 2009-10-28 | 福特全球技术公司 | 横摆稳定性控制*** |
JP4648283B2 (ja) * | 2006-10-12 | 2011-03-09 | 本田技研工業株式会社 | 車両制御装置 |
WO2012023305A1 (ja) * | 2010-08-20 | 2012-02-23 | 株式会社ユニバンス | 自動車 |
JP4929944B2 (ja) * | 2006-09-20 | 2012-05-09 | 日産自動車株式会社 | 車両の駆動力配分制御装置 |
CN106427957A (zh) * | 2015-08-11 | 2017-02-22 | 比亚迪股份有限公司 | 基于四轮驱动的电动汽车稳定控制***及方法及电动汽车 |
CN107089226A (zh) * | 2017-03-28 | 2017-08-25 | 江苏大学 | 一种用于车辆防抱死***的lqg滑移率控制器及其设计方法 |
CN108099901A (zh) * | 2017-12-18 | 2018-06-01 | 长春工业大学 | 一种提高极限工况下主动前轮转向控制精度的方法 |
-
2018
- 2018-07-06 CN CN201810735047.0A patent/CN109050659B/zh not_active Expired - Fee Related
Patent Citations (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN2435298Y (zh) * | 1999-01-07 | 2001-06-20 | 杨应捷 | 四轮汽车的四轮转向机构 |
US20010008986A1 (en) * | 1999-12-21 | 2001-07-19 | Brown Todd Allen | Roll over stability control for an automotive vehicle |
JP2007237933A (ja) * | 2006-03-08 | 2007-09-20 | Toyota Motor Corp | 車体姿勢制御装置 |
CN101400558A (zh) * | 2006-03-31 | 2009-04-01 | 丰田自动车株式会社 | 车辆运动控制装置及控制方法 |
JP4929944B2 (ja) * | 2006-09-20 | 2012-05-09 | 日産自動車株式会社 | 車両の駆動力配分制御装置 |
JP4648283B2 (ja) * | 2006-10-12 | 2011-03-09 | 本田技研工業株式会社 | 車両制御装置 |
CN101565043A (zh) * | 2008-04-25 | 2009-10-28 | 福特全球技术公司 | 横摆稳定性控制*** |
WO2012023305A1 (ja) * | 2010-08-20 | 2012-02-23 | 株式会社ユニバンス | 自動車 |
CN106427957A (zh) * | 2015-08-11 | 2017-02-22 | 比亚迪股份有限公司 | 基于四轮驱动的电动汽车稳定控制***及方法及电动汽车 |
CN107089226A (zh) * | 2017-03-28 | 2017-08-25 | 江苏大学 | 一种用于车辆防抱死***的lqg滑移率控制器及其设计方法 |
CN108099901A (zh) * | 2017-12-18 | 2018-06-01 | 长春工业大学 | 一种提高极限工况下主动前轮转向控制精度的方法 |
Non-Patent Citations (1)
Title |
---|
秦绪柏: "控制器参数对四轮转向和汽车稳定性控制的影响", 《机械制造与自动化》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111452801A (zh) * | 2019-01-21 | 2020-07-28 | 上海汽车集团股份有限公司 | 一种四轮转向汽车的鲁棒自适应控制方法和装置 |
CN111452801B (zh) * | 2019-01-21 | 2021-05-18 | 上海汽车集团股份有限公司 | 一种四轮转向汽车的鲁棒自适应控制方法和装置 |
CN110217229A (zh) * | 2019-06-25 | 2019-09-10 | 长春工业大学 | 一种适用于高速极限工况的路径跟踪控制方法 |
CN110539752A (zh) * | 2019-06-26 | 2019-12-06 | 江苏大学 | 一种智能汽车多预测范围模型预测轨迹跟踪控制方法及*** |
CN111055921A (zh) * | 2019-12-31 | 2020-04-24 | 吉林大学 | 一种基于数据驱动的四轮转向模型预测控制方法 |
CN117125079A (zh) * | 2023-08-28 | 2023-11-28 | 广州汽车集团股份有限公司 | 车辆的控制方法及装置、电子设备、存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN109050659B (zh) | 2020-09-18 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108099900B (zh) | 一种极限工况下保持汽车横向稳定的四轮转向控制方法 | |
CN108099902B (zh) | 一种体现汽车非线性特性的横摆稳定性控制方法 | |
CN108944866A (zh) | 一种改善转向与制动协同控制的自适应模型预测控制算法 | |
CN109050659A (zh) | 一种基于时变动力学模型的四轮转向汽车稳定性控制方法 | |
CN109050658A (zh) | 基于模型预测控制的汽车主动前轮转向自适应调节方法 | |
CN108082182B (zh) | 一种拓宽汽车稳定域的主动后轮转向控制方法 | |
CN108099901B (zh) | 一种提高极限工况下主动前轮转向控制精度的方法 | |
CN108107731A (zh) | 一种基于轮胎非线性特性的汽车稳定性控制方法 | |
Sharp et al. | A mathematical model for driver steering control, with design, tuning and performance results | |
You et al. | Multi-objective control synthesis: an application to 4WS passenger vehicles | |
CN108107732A (zh) | 主动前轮转向和直接横摆力矩联合的汽车稳定性控制方法 | |
Chen et al. | Robust control of reaction wheel bicycle robot via adaptive integral terminal sliding mode | |
CN110481343A (zh) | 四轮轮毂电机驱动汽车力矩补偿的组合二阶滑模控制方法 | |
CN114967475B (zh) | 一种无人驾驶车辆轨迹跟踪与稳定性鲁棒控制方法及*** | |
CN108099877A (zh) | 一种紧急制动工况下滑移率跟踪控制方法 | |
CN109094644A (zh) | 极限工况下的主动后轮转向与直接横摆力矩控制方法 | |
CN110239519A (zh) | 一种车辆爆胎过程的协调控制方法 | |
Cui et al. | Design optimization of a steering and suspension integrated system based on dynamic constraint analytical target cascading method | |
CN111679575A (zh) | 一种基于鲁棒模型预测控制的智能汽车轨迹跟踪控制器及其构造方法 | |
Ye et al. | The stability and accuracy analysis of automatic steering system with time delay | |
CN110126816A (zh) | 一种汽车主动前轮转向***的容错控制方法 | |
CN117369287A (zh) | 一种用于多轴车辆的鲁棒超螺旋滑模控制方法 | |
Fu et al. | Overview of longitudinal and lateral control for intelligent vehicle path tracking | |
CN116373846A (zh) | 一种基于bp神经网络优化的后轮转向车辆稳定性控制方法 | |
Zhai et al. | Handling stability control strategy for four‐wheel hub motor‐driven vehicle based on adaptive control for road adhesion |
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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20200918 Termination date: 20210706 |
|
CF01 | Termination of patent right due to non-payment of annual fee |