CN111965977A - 一种基于轮胎均等后备能力的汽车稳定性控制方法 - Google Patents
一种基于轮胎均等后备能力的汽车稳定性控制方法 Download PDFInfo
- Publication number
- CN111965977A CN111965977A CN202010786034.3A CN202010786034A CN111965977A CN 111965977 A CN111965977 A CN 111965977A CN 202010786034 A CN202010786034 A CN 202010786034A CN 111965977 A CN111965977 A CN 111965977A
- Authority
- CN
- China
- Prior art keywords
- tire
- force
- automobile
- vehicle
- follows
- 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 51
- 238000005457 optimization Methods 0.000 claims description 12
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000004422 calculation algorithm Methods 0.000 claims description 3
- 238000005096 rolling process Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 2
- 238000010586 diagram Methods 0.000 description 4
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B13/00—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion
- G05B13/02—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric
- G05B13/04—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators
- G05B13/042—Adaptive control systems, i.e. systems automatically adjusting themselves to have a performance which is optimum according to some preassigned criterion electric involving the use of models or simulators in which a parameter or coefficient is automatically adjusted to optimise the performance
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Artificial Intelligence (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Regulating Braking Force (AREA)
- Control Of Driving Devices And Active Controlling Of Vehicle (AREA)
Abstract
一种基于轮胎均等后备能力的汽车稳定性控制方法,其特征在于,该方法包括参考模型模块、轮胎数据处理模块、MPC控制器、CarSim汽车模型;参考模型模块用于计算参考横摆角速度和质心侧偏角;轮胎数据处理模块用于计算轮胎侧向力和轮胎状态刚度;CarSim汽车模型用于输出汽车的实际状态量;MPC控制器优化出前轮转角以及四个车轮的制动力矩,输入给CarSim汽车模型,实现转向和制动控制权的合理分配。
Description
技术领域:
本发明涉及汽车稳定性控制领域,更具体地说,涉及一种基于轮胎均等后备能力的汽车稳定性控制方法。
背景技术:
随着汽车行驶速度越来越高,汽车的行驶稳定性正受到人们越来越多的关注。其中,基于主动前轮转向(Active front steering,AFS)和直接横摆力矩(Direct yaw-moment control,DYC)集成的稳定性控制方法被广泛研究,也是今后汽车稳定性控制的发展方向。
AFS***能够通过前轮转角的主动干预纠正车辆的不足转向或过多转向,避免汽车发生侧向失稳,但轮胎侧向力接近附着极限时,AFS的控制性能也将接近极限。此时,可以利用DYC***产生的横摆力矩纠正车身姿态。因此,AFS与DYC的集成控制可以充分利用两者的优势,进一步提高车辆的侧向稳定性。但AFS与DYC对车辆的运动控制存在相互干涉和耦合,且产生转向的侧向力和产生制动的纵向力之间也存在相互影响,使得AFS与DYC集成控制时转向和制动的控制权分配问题一直缺乏有效的解决方法。
现有的AFS与DYC集成控制主要分为分层式集成控制和一体式集成控制两种。论文[Guo H,Feng L,Fang X,et al.Nonlinear model predictive lateral stabilitycontrol of active chassis for intelligent vehicles and its FPGAimplementation[J].IEEE Transactions on Systems Man&Cybernetics Systems,2017,49(1):2-13.]提出一种基于非线性模型预测控制的分层式AFS与DYC集成控制方法,上层根据参考横摆角速度计算出前轮转角和广义横摆力矩,下层根据制动策略将广义横摆力矩分配给相应的车轮进行制动。但前轮转角和横摆力矩的分配权重是恒定值,使得转向和制动的分配为固定比值。论文[Jalali M,Khosravani S,Khajepour A,et al.Model predictive control ofvehicle stability using coordinated active steering and differential brakes[J],2017,48(1):30-41.]提出了一种一体式车辆稳定性控制框架,直接将前轮转角和四个车轮的制动力矩作为控制输入,并通过软约束来避免车辆发生侧滑,为控制输入的全局最优性提供了可能。但上述方法没有考虑轮胎纵向力和侧向力的最优利用和分配,因此若能在现有一体式集成控制框架基础上,充分考虑轮胎力的最优利用和分配,将会使转向和制动的控制权分配更加合理,进一步提高汽车的侧向稳定性。
发明内容:
为解决现有汽车一体式稳定性集成控制时,因未考虑轮胎纵向力和侧向力的最优利用和分配而导致的汽车失稳问题,本发明提供一种基于轮胎均等后备能力的稳定性控制方法。首先,设计了一种基于轮胎均等后备能力的轮胎纵向力和侧向力分配方法,进行转向和制动的控制权分配。其次,采用线性时变的模型预测控制(记为LTV-MPC)设计基于轮胎均等后备能力的AFS和DYC一体式集成稳定性控制器,提高车辆的行驶稳定性。
本发明解决技术问题所采取的技术方案如下:
一种基于轮胎均等后备能力(记为ERCT)的稳定性控制方法,其特征在于,该方法包括参考模型模块、轮胎数据处理模块、MPC控制器、CarSim汽车模型;参考模型模块根据驾驶员输入δf,dri计算出参考横摆角速度γref和质心侧偏角βref;轮胎数据处理模块用于计算轮胎的侧向力和轮胎状态刚度;CarSim汽车模型用于输出汽车的实际状态量,包括汽车纵向速度、侧向速度、横摆角速度和质心侧偏角;MPC控制器根据参考横摆角速度γref、质心侧偏角βref和汽车的实际状态量,计算出前轮转角δf以及四个车轮的制动力矩Tb,ij,输入给CarSim汽车模型,从而实现转向和制动控制权的合理分配。
该方法包括以下步骤:
步骤1、设计参考模型模块,确定期望横摆角速度和质心侧偏角,其过程包括如下子步骤:
步骤1.1、采用线性二自由度汽车模型作为参考模型,其运动微分方程表达式如下:
其中:β是汽车质心侧偏角;γ是汽车横摆角速度;Iz是绕汽车质心的铅垂轴的横摆转动惯量;Ux是汽车纵向速度;lf和lr分别是汽车质心至前、后轴的距离;Cf和Cr分别是汽车前、后轮轮胎的侧偏刚度;δf是汽车前轮转角。
步骤1.2、根据步骤1.1中的运动微分方程(1),设计期望的汽车横摆角速度,其表达式如下:
步骤1.3、设计期望的汽车质心侧偏角,其表达式如下所示:
βref=0 (3)
其中,βref是汽车的期望质心侧偏角;质心侧偏角响应体现了车辆的稳定性,其值越小表示车辆侧向滑移越小,故将其参考值设为零;
步骤2、设计轮胎数据处理模块,其过程主要包括如下子步骤:
步骤2.1、设计非线性UniTire轮胎模型,如下:
其中,Fy为轮胎侧向力;为无量纲总切向力;φx为相对纵向滑移率;φn为修正后的相对综合滑移率;μy为侧向摩擦系数;Fz为轮胎垂直载荷;E为综合曲率因子;φ为相对综合滑移率;Ex为纵向力曲率因子;Ey为侧向力曲率因子;φy为相对侧向滑移率;λ为总切向力方向因子;Kx为纵滑刚度;Ky为侧偏刚度;Sx为纵向滑移率;Sy为侧向滑移率;μx为纵向摩擦系数;Fzn为轮胎无量纲垂直载荷;ω为车轮角速度;Vx为轮胎接地印迹中心纵向滑移速度;Vy为轮胎接地印迹中心侧向滑移速度;Re为轮胎有效滚动半径;κ为ISO轮胎坐标系中的纵向滑移率;α为ISO轮胎坐标系中的轮胎侧偏角;Fz0为轮胎标称载荷;模型参数:η=4,φc=1,pu1=0.89,pu2=-0.289,pu3=-0.289,su1=1.106,su2=-0.36,su3=0.051,pe1=-2.98,pe2=9.37,se1=-2.08,se2=9.37,pk1=0.019,pk2=-0.00000001,pk3=-0.00000001,sk1=0.046,sk2=0.008,sk3=0.006;
步骤2.2轮胎状态刚度预测,其过程包括如下子步骤:
步骤2.2.1、计算轮胎侧偏角,前、后轮轮胎侧偏角通过下式计算获得:
其中:αf和αr分别是汽车前、后轮轮胎的侧偏角;
步骤2.2.2为了获得轮胎的非线性特性,基于UniTire模型,获取不同轮胎载荷下轮胎状态刚度与轮胎侧偏角的关系曲线,得到轮胎状态刚度特性三维图。轮胎数据处理模块将轮胎的载荷和轮胎侧偏角分别输入到轮胎状态刚度特性三维图,通过线性插值法获得轮胎的轮胎状态刚度并输出给MPC控制器。在每个控制周期,轮胎数据处理模块更新一次轮胎的状态刚度。
步骤2.2.3轮胎模型线性化采用基于状态刚度的线性化方法,在预测时域内对轮胎模型线性化设计如下:
式中,
n=0,1,2……P
步骤3、设计MPC控制器,其过程包括如下子步骤:
步骤3.1.1、建立汽车动力学模型,其微分方程表达式如下:
其中,lf和lr分别为汽车质心到前轴和后轴的距离,w为轮距,Fy和Fx分别是纵向和侧向轮胎力,下标fl,fr,rl和rr分别指左前、右前、左后和右后车轮。
步骤3.1.2、建立车轮动力学模型,其表达式如下:
式中,Fb=-Fx;其中,Re为车轮的半径,Tb为制动力矩,ω为车轮的角速度。
将步骤3.1.1中的式(18)和步骤3.1.2中的(17)代入步骤2.2.1中的式(16),可得到基于线性时变的一体式AFS与DYC集成稳定性控制器模型如下:
式中,Cf=Cfl+Cfr。
将上式整理成标准的状态空间形式,如下:
对步骤3.1.2中的式(20)进行离散化,可以得到增量式的离散预测模型如下:
步骤3.2、设计基于轮胎均等后备能力的轮胎纵向力和侧向力分配方法,包括如下子步骤:
步骤3.2.1计算轮胎纵向力和侧向力的后备力,具体表达式如下:
步骤3.2.2轮胎后备力的分配,对于前轮转向汽车,轮胎纵向力和侧向力分配仅存在于车辆的前轮,转向和制动的相互干扰主要发生在车辆的横摆运动中,且由于力臂w/2和lf不等导致纵向力和侧向力对车辆的横摆运动贡献不同,因此需调整前轮转向和制动的控制权分配问题,使车辆前轮的纵向力和侧向力的后备力对保持车辆横摆稳定的能力相同,其公式如下:
将步骤3.2.1中的式(23)和步骤3.2.2中的式(24)联立,得到基于ERCT轮胎力分配方法中轮胎纵向力和侧向力之间的关系,如下所示:
基于附着椭圆曲线,得到基于ERCT的轮胎力分配平衡线。
步骤3.2.3根据车辆的运动状态确定是否需要轮胎力的分配;当车辆只有转向而无制动时,不需要进行轮胎力的分配。当前轮轮胎力达到基于ERCT轮胎力平衡线的初始平衡点时,轮胎力开始沿平衡线变化。
步骤3.2.4对于前轮转向的汽车,轮胎纵向力和侧向力分配仅存在于车辆的前轮。因此,左前轮和右前轮的轮胎力分配指标如下所示:
该指标的值越接近零,表示轮胎力的变化越接***衡线,即轮胎力分配越接近设计目标。
步骤3.2.5由于后轮轮胎不需要考虑轮胎力的分配,只需要满足轮胎力的利用率尽可能小,因此定义以下指标来描述后轮轮胎动力的利用情况,其表达式如下:
使用雅克比矩阵将前后轮的轮胎力分配指标进行线性化,具体如下:
式中,JF=[Jfl,Jfr,Jrl,Jrr]T,
将步骤3.1.2中的式(21)和步骤3.2.3中的式(30)联立,最终得到基于线性时变的一体式稳定性控制器***模型:
式中,
步骤3.3、计算预测输出,根据模型预测控制理论,取预测时域为P,控制时域为M,可以得到在当前k时刻的预测输出为:
Y(k+1|k)=Sdv,ξ·Δx(k)+I·y(k)+Sdv,u·ΔU(k)+Sdv,d·d(k) (32)
Sdv,u=[Su1,Su2,…,SuM]P×M
其中,
步骤3.4、设计优化目标及约束条件,其过程包括如下子步骤:
步骤3.4.1、用期望的横摆角速度和质心侧偏角与汽车实际的横摆角速度和质心侧偏角偏差的二范数,作为汽车稳定性能指标,其表达式如下:
其中:τγ和Γβ分别是对汽车横摆角速度和质心侧偏角跟踪性能的加权因子;
步骤3.4.2、使轮胎力变化尽可能的接***衡线,即轮胎力分配指标Jfl、Jfr、Jrl和Jrr值尽可能小,其表达式如下:
当汽车左前轮和右前轮需要制动时,则需要进行轮胎力的分配,对于暂时不需要轮胎力分配的车轮,可以将其在目标函数中对应分配指标的权重系数设为零。因此,目标函数中分配指标Jfl和Jfr的权重系数根据车辆状态设计为:
其中,υ是为避免不必要的频繁控制而设置的阈值,取υ=0.15·γref,max;ρk-1是上一步的ρ值。
步骤3.4.3、用控制量变化量二范数作为转向、制动平滑指标,其表达式如下:
步骤3.4.4、设置执行器物理约束,满足执行器要求:
利用线性不等式限制前轮转角及制动力变化量的上下限,得到转向、制动执行器的物理约束,其数学表达式为:
其中:-δfmax是前轮转角下限,δfmax是前轮转角上限;-Δδfmax是前轮转角变化量的下限;Δδfmax是前轮转角变化量的上限;-Tb,ijmax是制动力矩下限,Tb,ijmax是制动力矩上限;-ΔTb,ijmax是制动力矩变化量的下限;ΔTb,ijmax是制动力矩变化量的上限;
步骤3.4.5、设置控制输出约束,满足车辆的稳定性约束要求:
其中,-βmax是汽车质心侧偏角下限,βmax是汽车质心侧偏角上限;-γmax是汽车横摆角速度下限;γmax是汽车横摆角速度上限;
步骤3.5、求解***控制输入,其过程包括如下子步骤:
步骤3.5.1、利用线性加权法转化为单一指标,构建多目标优化控制问题:
服从于
i)预测模型式(31)
ii)约束条件式(35)~(36)
步骤3.5.2、采用二次规划算法,求解多目标优化控制问题(39),得到最优开环控制序制序列[Δδf,Tb,ij]T为:
本发明的有益效果是:
本发明提出了一种基于轮胎均等后备能力的新型轮胎纵向力和侧向力分配方法,使轮胎纵向力和侧向力的后备力对保持车辆横摆稳定的能力相同,实现了轮胎纵向力和侧向力的最优利用和分配,并设计了基于线性时变的一体式AFS与DYC集成控制器,能够有效的解决现有集成控制方法中转向和制动的控制权分配不合理的问题,明显提高车辆的侧向稳定性。
附图说明
图1是本发明的控制***结构示意图。
图2是汽车动力学模型示意图。
图3是车轮动力学模型示意图。
图4是轮胎状态刚度特性三维图。
图5是基于轮胎均等后备能力的轮胎力分配示意图。
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
一种基于轮胎均等后备能力的汽车稳定性控制方法,其特征在于,该方法包括参考模型模块1、轮胎数据处理模块2、MPC控制器3、CarSim汽车模型4;参考模型模块1用于根据驾驶员输入δf,dri计算出参考横摆角速度γref和质心侧偏角βref;轮胎数据处理模块2用于计算轮胎的侧向力和轮胎状态刚度;CarSim汽车模型4用于输出汽车的实际状态量,包括汽车纵向速度、侧向速度、横摆角速度和质心侧偏角;MPC控制器3根据参考横摆角速度γref、质心侧偏角βref和汽车的实际状态量,计算出前轮转角δf以及四个车轮的制动力矩Tb,ij,输入给CarSim汽车模型4,从而实现转向和制动控制权的合理分配。
下面以CarSim汽车仿真软件某车型为平台,具体说明本发明的方法,其主要参数如表1所示:
表1 CarSim汽车的主要参数
参考模型模块1的设计包括三部分:1.1建立线性二自由度汽车模型作为参考模型;1.2设计期望的汽车横摆角速度;1.3设计期望的汽车质心侧偏角;
在1.1部分中,采用线性二自由度汽车模型作为参考模型,其运动微分方程表达式如下:
其中:β是汽车质心侧偏角;γ是汽车横摆角速度;Iz是绕汽车质心的铅垂轴的横摆转动惯量;Ux是汽车纵向速度;lf和lr分别是汽车质心至前、后轴的距离;Cf和Cr分别是汽车前、后轮轮胎的侧偏刚度;δf是汽车前轮转角。
在1.2部分中,根据式(1)设计期望的汽车横摆角速度,其表达式如下:
在1.3部分中,设计期望的汽车质心侧偏角,其表达式如下所示:
βref=0 (3)
其中,βref是汽车的期望质心侧偏角;质心侧偏角响应体现了车辆的稳定性,其值越小表示车辆侧向滑移越小,故将其参考值设为零;
轮胎数据处理器模块2的设计包括两部分:2.1设计非线性轮胎模型;2.2轮胎状态刚度的预测;
在2.1部分中,设计非线性UniTire轮胎模型,如下:
其中,Fy为轮胎侧向力;为无量纲总切向力;φx为相对纵向滑移率;φn为修正后的相对综合滑移率;μy为侧向摩擦系数;Fz为轮胎垂直载荷;E为综合曲率因子;φ为相对综合滑移率;Ex为纵向力曲率因子;Ey为侧向力曲率因子;φy为相对侧向滑移率;λ为总切向力方向因子;Kx为纵滑刚度;Ky为侧偏刚度;Sx为纵向滑移率;Sy为侧向滑移率;μx为纵向摩擦系数;Fzn为轮胎无量纲垂直载荷;ω为车轮角速度;Vx为轮胎接地印迹中心纵向滑移速度;Vy为轮胎接地印迹中心侧向滑移速度;Re为轮胎有效滚动半径;κ为ISO轮胎坐标系中的纵向滑移率;α为ISO轮胎坐标系中的轮胎侧偏角;Fz0为轮胎标称载荷;模型参数η=4,φc=1,pu1=0.89,pu2=-0.289,pu3=-0.289,su1=1.106,su2=-0.36,su3=0.051,pe1=-2.98,pe2=9.37,se1=-2.08,se2=9.37,pk1=0.019,pk2=-0.00000001,pk3=-0.00000001,sk1=0.046,sk2=0.008,sk3=0.006;
在2.2部分中,轮胎状态刚度预测主要包括三部分:2.2.1计算轮胎前后轮的侧偏角;2.2.2轮胎状态刚度三维特性图的设计;2.2.3基于状态刚度法对轮胎模型进行线性化;
在2.2.1部分中,计算轮胎侧偏角,前、后轮轮胎侧偏角通过下式计算获得:
其中:αf和αr分别是汽车前、后轮轮胎的侧偏角。
在2.2.2部分中,为了获得轮胎的非线性特性,基于Unitire轮胎模型,获取不同轮胎载荷下的轮胎侧向刚度相对轮胎侧偏角的关系曲线,得到轮胎状态刚度特性三维图,如图4所示;轮胎数据处理器模块将轮胎的载荷和轮胎侧偏角分别输入到轮胎状态刚度特性三维图,通过线性插值法获得轮胎的轮胎状态刚度并输出给MPC控制器。在每个控制周期轮胎数据处理器模块更新一次轮胎的状态刚度。
在2.2.3部分中,轮胎模型线性化采用基于状态刚度的线性化方法,在预测时域内对轮胎模型线性化设计如下:
式中,
n=0,1,2……P
MPC控制器3的设计包括五部分:3.1建立预测模型;3.2基于轮胎均等后备能力的轮胎力分配方法的设计;3.3计算预测输出;3.4设计优化目标及约束条件;3.5求解***控制输入;
在3.1.1部分中,首先建立车辆动力学模型,如图2所示,其微分方程表达式如下:
其中,lf和lr分别为汽车质心到前轴和后轴的距离,w为轮距,Fy和Fx分别是纵向和侧向轮胎力,下标fl,fr,rl和rr分别指左前、右前、左后和右后车轮。
在3.1.2部分中,建立车轮动力学模型,如图3所示,其表达式如下:
式中,Fb=-Fx;其中,Re为车轮的半径,Tb为制动力矩,ω为车轮的角速度。
将式(18)和(17)代入式(16),可得到基于线性时变的一体式AFS与DYC集成稳定性控制器模型如下:
式中,Cf=Cfl+Cfr。
将式(20)整理成标准的状态空间形式,如下:
将式(21)进行离散化,可以得到增量式的离散预测模型如下:
在3.2部分中,设计基于轮胎均等后备能力的轮胎纵向力和侧向力分配方法,主要包括如下步骤:3.2.1计算轮胎纵向力和侧向力的后备力;3.2.2轮胎均等后备力的分配;3.2.3根据车辆的运动状态确定是否需要轮胎力的分配;3.2.4设计左前轮和右前轮的轮胎力分配指标;3.2.5设计后轮轮胎动力的利用指标;
在3.2.1部分中,计算轮胎纵向力和侧向力的后备力,具体表达式如下:
在3.2.2部分中,轮胎均等后备力的分配,对于前轮转向汽车,轮胎纵向力和侧向力分配仅存在于车辆的前轮,转向和制动的相互干扰主要发生在车辆的横摆运动中,且由于力臂w/2和lf不等导致纵向力和侧向力对车辆的横摆运动贡献不同,因此需调整前轮转向和制动的控制权分配问题,使车辆前轮的纵向力和侧向力的后备力对保持车辆横摆稳定的能力相同,其公式如下:
将式(23)和式(24)联立,得到基于ERCT轮胎力分配方法中轮胎纵向力和侧向力之间的关系,如下所示:
基于附着椭圆曲线,得到基于ERCT的轮胎力分配平衡线。
在3.2.3部分中,根据车辆的运动状态确定是否需要轮胎力的分配;当车辆只有转向而无制动时,不需要进行轮胎力的分配。当前轮轮胎力达到基于ERCT的轮胎力分平衡线线的初始平衡点时,轮胎力开始沿平衡线变化。
在3.2.4部分中,对于前轮转向的汽车,轮胎纵向力和侧向力分配仅存在于车辆的前轮。因此,左前轮和右前轮的轮胎力分配指标如下所示:
该指标的值越接近零,表示轮胎力的变化越接***衡线,即轮胎力分配越接近设计目标。
在3.2.5部分中,由于后轮轮胎不需要考虑轮胎力的分配,只需要满足轮胎力的利用率尽可能小,因此定义以下指标来描述后轮轮胎动力的利用情况,其表达式如下:
采用雅克比矩阵将前后轮的轮胎力分配指标进行线性化,具体如下:
式中,JF=[Jfl,Jfr,Jrl,Jrr]T,
将式(21)和式(30)联立,最终得到基于LTV-MPC的一体式稳定性控制器***模型:
式中,
在3.3部分中,计算预测输出,根据模型预测控制理论,取预测时域为P,控制时域为M,可以得到在当前k时刻的预测输出为:
Y(k+1|k)=Sdv,ξ·Δx(k)+I·y(k)+Sdv,u·ΔU(k)+Sdv,d·d(k) (32)
Sdv,u=[Su1,Su2,…,SuM]P×M
其中,
在3.4部分中,优化目标及约束条件设计包括三个部分:3.4.1设计跟踪性能指标;3.4.2设计轮胎力分配指标;3.4.3设计转向和制动平滑指标;3.4.4设置执行器物理约束;3.4.5设置控制输出约束;
在3.4.1部分中,用期望的横摆角速度和质心侧偏角与汽车实际的横摆角速度和质心侧偏角偏差的二范数作为汽车稳定性指标,其表达式如下:
其中:τγ和Γβ分别是对汽车横摆角速度和质心侧偏角跟踪性能的加权因子;
在3.4.2部分中,使轮胎力变化尽可能的接***衡线,其值应尽可能的小,其表达式如下:
当汽车左前轮和右前轮需要制动时,则需要进行轮胎力的分配,对于暂时不需要轮胎力分配的车轮,可以将其在目标函数中对应分配指标的权重系数设为零。因此,目标函数中分配指标Jfl和Jfr的权重系数根据车辆状态设计为:
其中,υ是为避免不必要的频繁控制而设置的阈值,取υ=0.15·γref,max;ρk-1是上一步的ρ值。
在3.4.3部分中,用控制量变化量二范数作为转向和制动平滑指标,其表达式如下:
在3.4.4部分中,设置执行器物理约束,满足执行器要求;利用线性不等式限制前轮转角及其变化量的上下限,得到转向、制动执行器的物理约束,其数学表达式为:
其中:-δfmax是前轮转角下限,δfmax是前轮转角上限;-Δδfmax是前轮转角变化量的下限;Δδfmax是前轮转角变化量的上限;-Tb,ijmax是制动力矩下限,Tb,ijmax是制动力矩上限;-ΔTb,ijmax是制动力矩变化量的下限;ΔTb,ijmax是制动力矩变化量的上限;
在3.4.5部分中,设置控制输出约束,满足车辆的稳定性约束要求:
其中,-βmax是汽车质心侧偏角下限,βmax是汽车质心侧偏角上限;-γmax是汽车横摆角速度下限;γmax是汽车横摆角速度上限;
在3.5部分中,求解***控制输入包括两部分:3.5.1构建多目标优化控制问题;3.5.2求解多目标优化控制问题;
在3.5.1部分中,利用线性加权法转化为单一指标,构建多目标优化控制问题:
服从于
i)预测模型式(31)
ii)约束条件式(35)~(36)
在3.5.2部分中,采用二次规划算法,求解多目标优化控制问题(37),得到最优开环控制序列[Δδf,Tb,ij]T为:
Claims (1)
1.一种基于轮胎均等后备能力的稳定性控制方法,其特征在于,该方法包括参考模型模块、轮胎数据处理模块、MPC控制器、CarSim汽车模型;参考模型模块根据驾驶员输入δf,dri计算出参考横摆角速度γref和质心侧偏角βref;轮胎数据处理模块用于计算轮胎的侧向力和轮胎状态刚度;CarSim汽车模型用于输出汽车的实际状态量,包括汽车纵向速度、侧向速度、横摆角速度和质心侧偏角;MPC控制器根据参考横摆角速度γref、质心侧偏角βref和汽车的实际状态量,计算出前轮转角δf以及四个车轮的制动力矩Tb,ij,输入给CarSim汽车模型,从而实现转向和制动控制权的合理分配;
该方法包括以下步骤:
步骤1、设计参考模型模块,确定期望横摆角速度和质心侧偏角,其过程包括如下子步骤:
步骤1.1、采用线性二自由度汽车模型作为参考模型,其运动微分方程表达式如下:
其中:β是汽车质心侧偏角;γ是汽车横摆角速度;Iz是绕汽车质心的铅垂轴的横摆转动惯量;Ux是汽车纵向速度;lf和lr分别是汽车质心至前、后轴的距离;Cf和Cr分别是汽车前、后轮轮胎的侧偏刚度;δf是汽车前轮转角;
步骤1.2、根据步骤1.1中的运动微分方程(1),设计期望的汽车横摆角速度,其表达式如下:
步骤1.3、设计期望的汽车质心侧偏角,其表达式如下所示:
βref=0 (3)
其中,βref是汽车的期望质心侧偏角;质心侧偏角响应体现了车辆的稳定性,其值越小表示车辆侧向滑移越小,故将其参考值设为零;
步骤2、设计轮胎数据处理模块,其过程主要包括如下子步骤:
步骤2.1、设计非线性UniTire轮胎模型,如下:
其中,Fy为轮胎侧向力;为无量纲总切向力;φx为相对纵向滑移率;φn为修正后的相对综合滑移率;μy为侧向摩擦系数;Fz为轮胎垂直载荷;E为综合曲率因子;φ为相对综合滑移率;Ex为纵向力曲率因子;Ey为侧向力曲率因子;φy为相对侧向滑移率;λ为总切向力方向因子;Kx为纵滑刚度;Ky为侧偏刚度;Sx为纵向滑移率;Sy为侧向滑移率;μx为纵向摩擦系数;Fzn为轮胎无量纲垂直载荷;ω为车轮角速度;Vx为轮胎接地印迹中心纵向滑移速度;Vy为轮胎接地印迹中心侧向滑移速度;Re为轮胎有效滚动半径;κ为ISO轮胎坐标系中的纵向滑移率;α为ISO轮胎坐标系中的轮胎侧偏角;Fz0为轮胎标称载荷;模型参数:η=4,φc=1,pu1=0.89,pu2=-0.289,pu3=-0.289,su1=1.106,su2=-0.36,su3=0.051,pe1=-2.98,pe2=9.37,se1=-2.08,se2=9.37,pk1=0.019,pk2=-0.00000001,pk3=-0.00000001,sk1=0.046,sk2=0.008,sk3=0.006;
步骤2.2轮胎状态刚度预测,其过程包括如下子步骤:
步骤2.2.1、计算轮胎侧偏角,前、后轮轮胎侧偏角通过下式计算获得:
其中:αf和αr分别是汽车前、后轮轮胎的侧偏角;
步骤2.2.2为了获得轮胎的非线性特性,基于UniTire模型,获取不同轮胎载荷下轮胎状态刚度与轮胎侧偏角的关系曲线,得到轮胎状态刚度特性三维图;轮胎数据处理模块将轮胎的载荷和轮胎侧偏角分别输入到轮胎状态刚度特性三维图,通过线性插值法获得轮胎的轮胎状态刚度并输出给MPC控制器;在每个控制周期,轮胎数据处理模块更新一次轮胎的状态刚度;
步骤2.2.3轮胎模型线性化采用基于状态刚度的线性化方法,在预测时域内对轮胎模型线性化设计如下:
式中,
n=0,1,2…P
步骤3、设计MPC控制器,其过程包括如下子步骤:
步骤3.1.1、建立汽车动力学模型,其微分方程表达式如下:
其中,lf和lr分别为汽车质心到前轴和后轴的距离,w为轮距,Fy和Fx分别是纵向和侧向轮胎力,下标fl,fr,rl和rr分别指左前、右前、左后和右后车轮;
步骤3.1.2、建立车轮动力学模型,其表达式如下:
式中,Fb=-Fx;其中,Re为车轮的半径,Tb为制动力矩,ω为车轮的角速度;
将步骤3.1.1中的式(18)和步骤3.1.2中的(17)代入步骤2.2.1中的式(16),可得到基于线性时变的一体式AFS与DYC集成稳定性控制器模型如下:
式中,Cf=Cfl+Cfr;
将上式整理成标准的状态空间形式,如下:
对步骤3.1.2中的式(20)进行离散化,可以得到增量式的离散预测模型如下:
步骤3.2、设计基于轮胎均等后备能力(equal reserve capacity of tire,ERCT)的轮胎纵向力和侧向力分配方法,包括如下子步骤:
步骤3.2.1计算轮胎纵向力和侧向力的后备力,具体表达式如下:
步骤3.2.2轮胎后备力的分配,对于前轮转向汽车,轮胎纵向力和侧向力分配仅存在于车辆的前轮,转向和制动的相互干扰主要发生在车辆的横摆运动中,且由于力臂w/2和lf不等导致纵向力和侧向力对车辆的横摆运动贡献不同,因此需调整前轮转向和制动的控制权分配问题,使车辆前轮的纵向力和侧向力的后备力对保持车辆横摆稳定的能力相同,其公式如下:
将步骤3.2.1中的式(23)和步骤3.2.2中的式(24)联立,得到基于ERCT的轮胎力分配方法中轮胎纵向力和侧向力之间的关系,如下所示:
基于附着椭圆曲线,得到基于ERCT的轮胎力分配平衡线;
步骤3.2.3根据车辆的运动状态确定是否需要轮胎力的分配;当车辆只有转向而无制动时,不需要进行轮胎力的分配;当前轮轮胎力达到基于ERCT的轮胎力平衡线的初始平衡点时,轮胎力开始沿平衡线变化;
步骤3.2.4对于前轮转向的汽车,轮胎纵向力和侧向力分配仅存在于车辆的前轮;因此,左前轮和右前轮的轮胎力分配指标如下所示:
该指标的值越接近零,表示轮胎力的变化越接***衡线,即轮胎力分配越接近设计目标;
步骤3.2.5由于后轮轮胎不需要考虑轮胎力的分配,只需要满足轮胎力的利用率尽可能小,因此定义以下指标来描述后轮轮胎动力的利用情况,其表达式如下:
使用雅克比矩阵将前后轮的轮胎力分配指标进行线性化,具体如下:
式中,JF=[Jfl,Jfr,Jrl,Jrr]T,
将步骤3.1.2中的式(21)和步骤3.2.3中的式(30)联立,最终得到基于线性时变的一体式稳定性控制器***模型:
式中,
步骤3.3、计算预测输出,根据模型预测控制理论,取预测时域为P,控制时域为M,可以得到在当前k时刻的预测输出为:
Y(k+1|k)=Sdv,ξ·Δx(k)+I·y(k)+Sdv,u·ΔU(k)+Sdv,d·d(k) (32)
Sdv,u=[Su1,Su2,…,SuM]P×M
其中,
步骤3.4、设计优化目标及约束条件,其过程包括如下子步骤:
步骤3.4.1、用期望的横摆角速度和质心侧偏角与汽车实际的横摆角速度和质心侧偏角偏差的二范数,作为汽车稳定性能指标,其表达式如下:
其中:τγ和Γβ分别是对汽车横摆角速度和质心侧偏角跟踪性能的加权因子;
步骤3.4.2、使轮胎力变化尽可能的接***衡线,即轮胎力分配指标Jfl、Jfr、Jrl和Jrr值尽可能小,其表达式如下:
当汽车左前轮和右前轮需要制动时,则需要进行轮胎力的分配,对于暂时不需要轮胎力分配的车轮,可以将其在目标函数中对应分配指标的权重系数设为零;因此,目标函数中分配指标Jfl和Jfr的权重系数根据车辆状态设计为:
其中,υ是为避免不必要的频繁控制而设置的阈值,取υ=0.15·γref,max;ρk-1是上一步的ρ值;
步骤3.4.3、用控制量变化量二范数作为转向、制动平滑指标,其表达式如下:
步骤3.4.4、设置执行器物理约束,满足执行器要求:
利用线性不等式限制前轮转角及制动力变化量的上下限,得到转向、制动执行器的物理约束,其数学表达式为:
其中:-δfmax是前轮转角下限,δfmax是前轮转角上限;-Δδfmax是前轮转角变化量的下限;Δδfmax是前轮转角变化量的上限;-Tb,ijmax是制动力矩下限,Tb,ijmax是制动力矩上限;-ΔTb,ijmax是制动力矩变化量的下限;ΔTb,ijmax是制动力矩变化量的上限;
步骤3.4.5、设置控制输出约束,满足车辆的稳定性约束要求:
其中,-βmax是汽车质心侧偏角下限,βmax是汽车质心侧偏角上限;-γmax是汽车横摆角速度下限;γmax是汽车横摆角速度上限;
步骤3.5、求解***控制输入,其过程包括如下子步骤:
步骤3.5.1、利用线性加权法转化为单一指标,构建多目标优化控制问题:
服从于
i)预测模型式(31)
ii)约束条件式(35)~(36)
步骤3.5.2、采用二次规划算法,求解多目标优化控制问题(39),得到最优开环控制序制序列[Δδf,Tb,ij]T为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010786034.3A CN111965977B (zh) | 2020-08-06 | 2020-08-06 | 一种基于轮胎均等后备能力的汽车稳定性控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010786034.3A CN111965977B (zh) | 2020-08-06 | 2020-08-06 | 一种基于轮胎均等后备能力的汽车稳定性控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111965977A true CN111965977A (zh) | 2020-11-20 |
CN111965977B CN111965977B (zh) | 2023-01-10 |
Family
ID=73365084
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010786034.3A Active CN111965977B (zh) | 2020-08-06 | 2020-08-06 | 一种基于轮胎均等后备能力的汽车稳定性控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111965977B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113954873A (zh) * | 2021-11-01 | 2022-01-21 | 北京科技大学 | 一种基于降维模型的转向和制动一体化路径跟踪控制方法 |
CN113954821A (zh) * | 2021-11-01 | 2022-01-21 | 北京科技大学 | 一种转向和扭矩矢量一体化车辆稳定性控制方法 |
CN114368369A (zh) * | 2022-03-22 | 2022-04-19 | 北京航空航天大学 | 一种适应路面摩擦系数的底盘制动***集成控制方法 |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020109402A1 (en) * | 2000-11-20 | 2002-08-15 | Toyota Jidosha Kabushiki Kaisha | Vehicle motion control device and method |
US20070106446A1 (en) * | 2005-02-11 | 2007-05-10 | Honeywell International Inc. | Method and system using tire stretch data to control braking |
JP2010081720A (ja) * | 2008-09-25 | 2010-04-08 | Mazda Motor Corp | 車両用駆動力制御装置 |
KR20100093672A (ko) * | 2009-02-17 | 2010-08-26 | 주식회사 만도 | 차량 자세제어장치 및 그 제어방법 |
CN106458210A (zh) * | 2014-09-29 | 2017-02-22 | 日立汽车***株式会社 | 车辆的运动控制方法、装置以及搭载该装置的车辆 |
CN106467111A (zh) * | 2015-08-20 | 2017-03-01 | 比亚迪股份有限公司 | 车身稳定控制方法、***及汽车 |
CN108107731A (zh) * | 2017-12-18 | 2018-06-01 | 长春工业大学 | 一种基于轮胎非线性特性的汽车稳定性控制方法 |
CN108107732A (zh) * | 2017-12-18 | 2018-06-01 | 长春工业大学 | 主动前轮转向和直接横摆力矩联合的汽车稳定性控制方法 |
CN108181811A (zh) * | 2017-12-18 | 2018-06-19 | 长春工业大学 | 一种基于线性时变的滑移率跟踪控制方法 |
CN108944866A (zh) * | 2018-07-06 | 2018-12-07 | 长春工业大学 | 一种改善转向与制动协同控制的自适应模型预测控制算法 |
US20190210591A1 (en) * | 2018-01-08 | 2019-07-11 | Sri International | System and Method to Enhance the Driving Performance of a Leanable Vehicle |
CN111267835A (zh) * | 2020-03-26 | 2020-06-12 | 桂林电子科技大学 | 基于模型预测算法的四轮独立驱动汽车稳定性控制方法 |
CN111391822A (zh) * | 2020-03-27 | 2020-07-10 | 吉林大学 | 一种极限工况下汽车横纵向稳定性协同控制方法 |
-
2020
- 2020-08-06 CN CN202010786034.3A patent/CN111965977B/zh active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020109402A1 (en) * | 2000-11-20 | 2002-08-15 | Toyota Jidosha Kabushiki Kaisha | Vehicle motion control device and method |
US20070106446A1 (en) * | 2005-02-11 | 2007-05-10 | Honeywell International Inc. | Method and system using tire stretch data to control braking |
JP2010081720A (ja) * | 2008-09-25 | 2010-04-08 | Mazda Motor Corp | 車両用駆動力制御装置 |
KR20100093672A (ko) * | 2009-02-17 | 2010-08-26 | 주식회사 만도 | 차량 자세제어장치 및 그 제어방법 |
CN106458210A (zh) * | 2014-09-29 | 2017-02-22 | 日立汽车***株式会社 | 车辆的运动控制方法、装置以及搭载该装置的车辆 |
CN106467111A (zh) * | 2015-08-20 | 2017-03-01 | 比亚迪股份有限公司 | 车身稳定控制方法、***及汽车 |
CN108107731A (zh) * | 2017-12-18 | 2018-06-01 | 长春工业大学 | 一种基于轮胎非线性特性的汽车稳定性控制方法 |
CN108107732A (zh) * | 2017-12-18 | 2018-06-01 | 长春工业大学 | 主动前轮转向和直接横摆力矩联合的汽车稳定性控制方法 |
CN108181811A (zh) * | 2017-12-18 | 2018-06-19 | 长春工业大学 | 一种基于线性时变的滑移率跟踪控制方法 |
US20190210591A1 (en) * | 2018-01-08 | 2019-07-11 | Sri International | System and Method to Enhance the Driving Performance of a Leanable Vehicle |
CN108944866A (zh) * | 2018-07-06 | 2018-12-07 | 长春工业大学 | 一种改善转向与制动协同控制的自适应模型预测控制算法 |
CN111267835A (zh) * | 2020-03-26 | 2020-06-12 | 桂林电子科技大学 | 基于模型预测算法的四轮独立驱动汽车稳定性控制方法 |
CN111391822A (zh) * | 2020-03-27 | 2020-07-10 | 吉林大学 | 一种极限工况下汽车横纵向稳定性协同控制方法 |
Non-Patent Citations (9)
Title |
---|
BINGTAOREN等: "MPC-based yaw stability control in in-wheel-motored EV via active front steering and motor torque distribution", 《MECHATRONICS》 * |
LU DENG 等: "A Multi-Point Tire Model for Studying Bridge–Vehicle Coupled Vibration", 《INTERNATIONAL JOURNAL OF STRUCTURAL STABILITY AND DYNAMICS》 * |
张亮修等: "自主车辆线性时变模型预测路径跟踪控制", 《同济大学学报(自然科学版)》 * |
张细政等: "基于转矩协调分配的分布式驱动电动汽车稳定性控制", 《中国机械工程》 * |
李绍松 等: "极限工况下主动前轮转向汽车稳定性控制", 《汽车工程》 * |
池成: "基于线控全轮转向驱动协调的轮毂电动汽车操控稳定性控制研究", 《中国优秀硕士学位论文全文数据库 (工程科技Ⅱ辑)》 * |
郑顺航: "汽车高速紧急避让转向制动协调控制研究", 《中国优秀硕士学位论文全文数据库 (工程科技Ⅱ辑)》 * |
郭孔辉等: "UniTire统一轮胎模型", 《机械工程学报》 * |
陈杰等: "基于LTV-MPC的车辆稳定性控制研究", 《汽车工程》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113954873A (zh) * | 2021-11-01 | 2022-01-21 | 北京科技大学 | 一种基于降维模型的转向和制动一体化路径跟踪控制方法 |
CN113954821A (zh) * | 2021-11-01 | 2022-01-21 | 北京科技大学 | 一种转向和扭矩矢量一体化车辆稳定性控制方法 |
CN114368369A (zh) * | 2022-03-22 | 2022-04-19 | 北京航空航天大学 | 一种适应路面摩擦系数的底盘制动***集成控制方法 |
CN114368369B (zh) * | 2022-03-22 | 2022-06-07 | 北京航空航天大学 | 一种适应路面摩擦系数的底盘制动***集成控制方法 |
Also Published As
Publication number | Publication date |
---|---|
CN111965977B (zh) | 2023-01-10 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111965977B (zh) | 一种基于轮胎均等后备能力的汽车稳定性控制方法 | |
CN108944866B (zh) | 一种改善转向与制动协同控制的自适应模型预测控制算法 | |
CN111806427B (zh) | 一种四轮毂电机驱动车辆的综合控制方法 | |
CN111959500B (zh) | 一种基于轮胎力分配的汽车路径跟踪性能提升方法 | |
CN110696793B (zh) | 一种四轮转向联合差动制动的智能车辆分层控制方法 | |
CN108107732B (zh) | 主动前轮转向和直接横摆力矩联合的汽车稳定性控制方法 | |
CN108177692B (zh) | 一种电动轮驱动汽车差动助力转向与稳定性协调控制方法 | |
CN108107731B (zh) | 一种基于轮胎非线性特性的汽车稳定性控制方法 | |
CN112572411B (zh) | 一种考虑轮胎侧偏特性的车辆底盘协调控制方法及*** | |
CN109291932B (zh) | 基于反馈的电动汽车横摆稳定性实时控制装置及方法 | |
CN112026533B (zh) | 一种极限工况下的四轮独立驱动电动汽车牵引力控制方法 | |
CN112829766B (zh) | 一种基于分布式驱动电动车辆的自适应路径跟踪方法 | |
JP2007269295A (ja) | 車両運動制御装置及び制御方法 | |
CN112373459B (zh) | 一种四轮毂电机驱动车辆上层运动状态控制方法 | |
WO2022266824A1 (zh) | 一种转向控制方法及装置 | |
CN113978263A (zh) | 一种驱动轮防滑与转矩优化融合的电动汽车稳定性控制方法 | |
CN113147422A (zh) | 一种分布式驱动电动汽车直接横摆力矩控制***及方法 | |
CN116279525A (zh) | 基于模型预测控制的分布式电动汽车集中式姿态控制方法 | |
JP2009502621A (ja) | 車両の操舵方向の制御方法 | |
JP4796480B2 (ja) | 車両運動制御装置及び制御方法 | |
CN113830075B (zh) | 车辆稳定控制方法、装置、电子设备及介质 | |
CN114834263A (zh) | 电动汽车主动前轮转向与扭矩矢量的协调控制方法及装置 | |
CN110758404B (zh) | 一种针对车辆稳定***执行器故障的容错控制方法 | |
JP3212134B2 (ja) | 車両の統合制御装置 | |
CN113954821B (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 |