CN110509915B - 一种基于时变速度的四轮驱动汽车横向稳定性控制方法 - Google Patents

一种基于时变速度的四轮驱动汽车横向稳定性控制方法 Download PDF

Info

Publication number
CN110509915B
CN110509915B CN201910137939.5A CN201910137939A CN110509915B CN 110509915 B CN110509915 B CN 110509915B CN 201910137939 A CN201910137939 A CN 201910137939A CN 110509915 B CN110509915 B CN 110509915B
Authority
CN
China
Prior art keywords
vehicle
automobile
time
wheel
wheel drive
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
CN201910137939.5A
Other languages
English (en)
Other versions
CN110509915A (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 CN201910137939.5A priority Critical patent/CN110509915B/zh
Publication of CN110509915A publication Critical patent/CN110509915A/zh
Application granted granted Critical
Publication of CN110509915B publication Critical patent/CN110509915B/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, or advanced driver assistance systems for ensuring comfort, stability and safety or drive control systems for propelling or retarding the vehicle
    • 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
    • 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

Abstract

本发明涉及汽车稳定性控制技术领域,具体涉及一种基于时变速度的四轮驱动汽车横向稳定性控制方法,包括以下步骤:T1采集汽车的固有参数和实时参数;T2利用T1的参数,建立汽车控制模型系统;T3对T2系统考虑时变速度设计LPV控制器;T4根据T3结果对四轮驱动汽车的转矩分配;本发明考虑汽车速度为时变的,通过设计的控制器计算出需要提供的合适的横摆力矩,控制上述两个参数在一个合适的范围内;然后通过下层控制器将计算出的横摆力矩分配到四个车轮上,提高汽车的操纵性能,改善汽车的横向稳定性能,降低汽车运行过程中发生危险情况的可能,提高乘坐舒适度,具有很强的创造性。

Description

一种基于时变速度的四轮驱动汽车横向稳定性控制方法
技术领域
本发明涉及汽车稳定性控制技术领域,具体涉及一种基于时变速度的四轮驱动汽车横向稳定性控制方法。
背景技术
能源短缺和生态恶化是人类面临的两大难题。石油资源作为地球上有限的自然资源,正在逐渐消耗殆尽。而采用汽油作为燃料的车辆,其尾气排放也是不可避免的问题。随着人民生活的富裕,传统汽车的保有量迅速上升,随之而来的是汽车尾气排放的剧增,这已成为雾霾的主要成因之一。而雾霾对人和生态环境都是一种巨大的威胁。因此,开发新能源汽车,减少排放已经成为全球范围内的一种大趋势。而电能作为一种可靠、来源广泛、对环境友好的能源,将是未来的汽车的主要发展方向。且采用电能的汽车可以做到零排放,对生态环境也是非常友好的。电动汽车的发展方向中,四轮驱动的轮毂电机驱动汽车是一个重要的分支。轮毂电机驱动电动汽车具有很多优点,主要包括:(1)底盘主动安全集成,轮毂电机驱动电动汽车可将动力装量灵活地布置在四个轮、前轮或者后轮上,通过对轮毂电机的单独控制或者协调控制,可实现更加有效地控制,提高了整车动力性和稳定性;(2)整车空间布置灵活,省略了传统汽车上的换挡装置、离合器、变速器、传动轴和差速器等器件,降低了整车质量,使得利用空间更大; (3)节能环保,具有制动能量回收,采用的轮毂电机具有零排放且低噪声等优点。
而汽车的安全性技术经过了被动安全性技术的发展,现在主要是针对主动安全性技术进行研究。其中ABS,ASR是主动安全技术典型的代表。汽车横向稳定性是保证汽车高速安全行驶的一项非常重要的性能,以提高汽车主动安全性为目的的汽车横向稳定性控制系统的研究己成为汽车技术方面的研究的一个重点。
ABS/ASR可在车辆紧急制动或加速行驶状态时通过控制车轮滑移率来提高车辆的安全性,但对于紧急转向或高速换线行驶的车辆则无法提供安全保障。汽车直接横摆力矩控制(简称DYC,Direct Yaw-moment Control)的研究就是旨在淡化驾驶人员的操作技能对车辆运动安全性的影响,在车辆的各种行驶状态下通过对每个车轮的受力进行调节,产生横摆力矩,克服过多转向或不足转向,从而主动地对车辆进行动力学控制,提高汽车在高速和恶劣道路等极限条件下行驶时的操纵稳定性。
与本发明最近似的实现方案:CN105691381B-一种四轮独立驱动电动汽车稳定性控制方法及系统。此已有的发明专利也是针对汽车横向动力学稳定性的控制方法,但所用方法为滑模控制。其用滑模控制的方法求得理想的横摆角速度,但是没有考虑汽车行驶中速度时变的情况。
汽车行驶过程中,为了保证安全性和舒适性,需要控制汽车的质心侧偏角和横摆角速度在一定的范围内。当汽车高速运行的时候,要求汽车的响应速度要快、稳定性要好、安全性要高,当汽车低速运行的时候,要求汽车在保持横向稳定性的同时还能提高电动汽车的机动性能,因此要求汽车质心侧偏角的响应值应该尽量趋近于零。而横摆角速度是指汽车绕垂直轴的偏转,该偏转的大小代表汽车的稳定程度。如果偏转角速度达到一个阈值,说明汽车发生测滑或甩尾等危险工况的可能性较大。因此汽车质心侧偏角和横摆角速度这两个参数在汽车运行过程中是需要加以关注的。
发明内容
针对现有技术的不足,本发明公开了一种基于时变速度的四轮驱动汽车横向稳定性控制方法,通过设计的控制器计算出需要提供的合适的横摆力矩,控制上述两个参数在一个合适的范围内;然后通过下层控制器将计算出的横摆力矩分配到四个车轮上,提高汽车的操纵性能,改善汽车的横向稳定性能,降低汽车运行过程中发生危险情况的可能,提高乘坐舒适度。
本发明通过以下技术方案予以实现:
一种基于时变速度的四轮驱动汽车横向稳定性控制方法,其特征在于,所述控制方法包括以下步骤:
T1采集汽车的固有参数和实时参数;
T2利用T1的参数,建立汽车控制模型系统;
T3对T2系统考虑时变速度设计LPV控制器;
T4根据T3结果对四轮驱动汽车的转矩分配。
优选的,所述T1中,
所述固有参数包括汽车总质量m,转动惯量Iz,前、后轴到汽车重心的距离lf、 lr,以及前、后轮的侧偏刚度Cf、Cr;这些参数是汽车固有的,提前采集并存储;
所述实时参数包括前轮转向角δ,汽车纵向速度vx,横摆角速度r,质心侧偏角β;
所述实时参数中,前轮转向角通过传感器采集方向盘转角,然后通过转向传动机构的参数来计算获得;横摆角速度通过陀螺仪测算获得;车身纵向速度和质心侧偏角通过状态参数估计器结合方向盘转角、横摆角速度、车身各方向加速度计算得到,车身各方向加速度由陀螺仪得到。
优选的,所述T2中,首先建立汽车二自由度模型,CG表示汽车的重心,m为汽车总质量,β为汽车质心侧偏角,r为横摆角速度,
Figure RE-GDA0002229102490000034
为横摆角加速度,Fyf和Fyr分别为前、后轮的侧向力,vy和vx分别是汽车横向和纵向速度,
Figure RE-GDA0002229102490000035
为汽车横向加速度, Iz为车身转动惯量,ΔMz是外部提供的横摆力矩,通过转向系统和制动系统提供, lf、lr分别是前后轴到汽车重心的距离;
根据牛顿力学定律,得到以下方程:
Figure RE-GDA0002229102490000031
Figure RE-GDA0002229102490000032
在上面的方程中:
Fyf=Cfαf,Fyr=Crαr
Figure RE-GDA0002229102490000033
其中,Cf、Cr分别为前后轮胎的侧偏刚度,αfr分别为前、后轮侧偏角。
优选的,选取质心侧偏角β和横摆角速度r作为控制模型的状态变量,得到如下汽车横向运动模型:
Figure RE-GDA0002229102490000045
其中:
Figure RE-GDA0002229102490000041
Figure RE-GDA0002229102490000042
Figure RE-GDA0002229102490000043
为了获得更好的操纵性和稳定性,选取质心侧偏角和横摆角速度为控制输出变量,则有:
z1(t)=β(t)
z2(t)=r(t)
其中,β(t)表示质心侧偏角β随时间t的变化,r(t)表示横摆角速度r随时间t的变化;
令q1=1/vx,q2=1/vx 2,则系统描述为:
Figure RE-GDA0002229102490000044
z1(t)=C1x(t)
z2(t)=C2x(t)
其中:
C1=[1 0],C2=[0 1]
Figure RE-GDA0002229102490000051
优选的,所述T3中,假设汽车纵向速度vx是实时测的,且在范围[vx min vx max] 内变化,则qi(i=1,2)在范围[qi min qi max]内变化,设变量
Figure RE-GDA0002229102490000052
则q的值就会在一个多胞体范围内变化,记为:
Figure RE-GDA0002229102490000053
其中:
α1=x(1-y),α2=(1-x)y,α3=xy,α4=(1-x)(1-y),
Figure RE-GDA0002229102490000054
Figure RE-GDA0002229102490000055
考虑控制器结构为u=Kx(t),则系统描述为
Figure RE-GDA0002229102490000056
z1(t)=C1x(t)
z2(t)=C2x(t)
其中
Figure RE-GDA0002229102490000057
K为状态反馈增益矩阵;
Figure RE-GDA0002229102490000058
为x(t)的一阶导数,外部横摆力矩ΔMz=u。
优选的,假设汽车纵向速度是时变的,建立的系统模型就变为连续时间的多胞体LPV系统;
分别选取H和GH2范数来表征输出z1和z2,由于实际情况中前轮转向角δ是有一定限度的,因此考虑δ∈L2,L2表示前轮转向角的取值范围;||z1||2,||δ||2分别表示z1和δ的二范数,||z2||表示z2的无穷范数;||T1||表示系统的H增益,
Figure RE-GDA0002229102490000059
表示系统的GH2增益(广义H2增益);
Figure RE-GDA0002229102490000061
状态反馈增益矩阵K的取值由如下方法确定:针对上述LPV系统,若存在矩阵 X>0和Wi>0,i=1,2,3,4.使得以下LMI成立,则存在LPV控制器使得系统二次稳定,且满足
Figure RE-GDA0002229102490000062
Figure RE-GDA0002229102490000063
Figure RE-GDA0002229102490000064
Figure RE-GDA0002229102490000065
通过求解以上矩阵不等式,求得X和Wi,然后求得
Figure RE-GDA0002229102490000066
其中 Ki=WiX-1;然后通过下式计算汽车的横摆力矩ΔMz
ΔMz=Kx(t)。
优选的,所述T4中,四轮驱动汽车的转矩分配采用基于目标优化的方法进行转矩分配;
对车辆运动模型进行分析,得到:
ΔMz=Fx1(-cosδ·ls+sinδ·lf)+Fx2(cosδ·ls+sinδ·lf)-Fx3·ls+Fx4·ls
其中Fxi,i=1,2,3,4为四个车轮的纵向轮胎力,δ为前轮转向角,ls为车辆横轴长度的一半
转矩分配策略还需满足纵向驱动合力Fx和四个车轮的纵向轮胎力Fxi之间的关系如下:
Fx=Fx1·cosδ+Fx2·cosδ+Fx3+Fx4
对上面的式子进行整理,得y=Gξ,其中:
y=[ΔMz Fx]T,ξ=[Fx1 Fx2 Fx3 Fx4]T
Figure RE-GDA0002229102490000071
假设要使车辆稳定所需要的外部横摆力矩为ΔMdes,Fx=max,那么y的期望值为yd=[ΔMdes Fx]T
优选的,转矩分配优化目标函数由两部分组成:
(1)考虑转矩的分配误差,得到目标函数的第一部分
J1=(Gξ-yd)TQ1(Gξ-yd),其中Q1=diag{q11 q12}为权重因子;
(2)考虑轮胎的轮胎利用率为
Figure RE-GDA0002229102490000072
其中i=1,2,3,4.μi为第i个轮胎与路面之间的摩擦系数,Fzi为每个轮胎的垂向载荷,分别为:
Figure RE-GDA0002229102490000073
Figure RE-GDA0002229102490000074
其中g为重力加速度,ax和ay分别为车辆的纵向加速度和横向加速度,h 为车辆的重心高度,L=lf+lr;然后得目标函数的第二部分为:
Figure RE-GDA0002229102490000075
其中
Figure RE-GDA0002229102490000076
最终的目标函数为J=J1+J2;通过最优化此目标函数得到理想的转矩。
优选的,考虑车轮的纵向力与转矩之间的关系为:
Ti=Fxi·R,i=1,2,3,4.R为车轮转动半径,Ti为每个轮内电机产生的转矩,在实际运行中,假设每个轮内电机产生的最大转矩为Tu,最小为Tl,即Tl≤Ti≤Tu,那么得到轮胎力的约束为:
Figure RE-GDA0002229102490000081
综上将转矩分配问题描述为优化问题如下:
Figure RE-GDA0002229102490000082
上式利用matlab软件进行优化求解,得到ξ,也即得到了四个轮内机应输出的力的大小,通过电子系统控制电机按需求输出。
本发明的有益效果为:
与现有技术相比,本发明主要考虑了汽车的速度作为一个时变参数,同时现有技术大多采用滑模控制方法,而本发明采用了LPV控制器,并且考虑了转矩在四轮上的分配。
1.考虑汽车纵向速度为时变参数,设计了多目标LPV控制器,鲁棒性更好;
2.用最优化方法来设计下层横摆力矩分配系统,减小分配误差并最大化轮胎利用率,保证稳定;
3.本发明方法有效改善了汽车的横向稳定性能,降低汽车运行过程中发生危险情况的可能,提高乘坐舒适度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还根据这些附图获得其他的附图。
图1是线性二自由度汽车模型;
图2是四轮车辆动力模型图;
图3是汽车横摆力矩稳定性控制组成结构图;
图4是控制系统结构框图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例;基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本实施例提出了一种基于时变速度的横摆力矩控制方法,考虑汽车速度为时变的,通过设计的控制器计算出需要提供的合适的横摆力矩,控制上述两个参数在一个合适的范围内;然后通过下层控制器将计算出的横摆力矩分配到四个车轮上,提高汽车的操纵性能,改善汽车的横向稳定性能,降低汽车运行过程中发生危险情况的可能,提高乘坐舒适度。
步骤1:采集汽车的参数
需要采集的参数包括固有参数和实时参数,其中:
固有参数:包括汽车总质量m,转动惯量Iz,前、后轴到汽车重心的距离lf、lr,以及前、后轮的侧偏刚度Cf、Cr;这些参数是汽车固有的,可提前采集并存储。
实时参数:包括前轮转向角δ,汽车纵向速度vx,横摆角速度r,质心侧偏角β;这些参数中,前轮转向角可通过传感器采集方向盘转角,然后通过转向传动机构的参数来计算获得;横摆角速度可通过陀螺仪测算获得;车身纵向速度和质心侧偏角可通过状态参数估计器结合方向盘转角、横摆角速度、车身各方向加速度计算得到,车身各方向加速度可由陀螺仪得到。
步骤2:利用汽车的参数,建立汽车控制模型系统
在该步骤中,首先建立汽车二自由度模型,模型的示意图如图2所示;图中, CG表示汽车的重心,m为汽车总质量,β为汽车质心侧偏角,r为横摆角速度,r 为横摆角加速度(横摆角速度的一阶导数),Fyf和Fyr分别为前、后轮的侧向力,vy和vx分别是汽车横向和纵向速度,
Figure RE-GDA0002229102490000108
为汽车横向加速度,Iz为车身转动惯量,ΔMz是外部提供的横摆力矩,可以通过转向系统和制动系统提供,lf、lr分别是前后轴到汽车重心的距离。
根据牛顿力学定律,可得到以下方程:
Figure RE-GDA0002229102490000101
Figure RE-GDA0002229102490000102
在上面的方程中:
Fyf=Cfαf,Fyr=Crαr
Figure RE-GDA0002229102490000103
其中,Cf、Cr分别为前后轮胎的侧偏刚度,αfr分别为前、后轮侧偏角;
选取质心侧偏角β和横摆角速度r作为控制模型的状态变量,可以得到如下汽车横向运动模型:
Figure RE-GDA0002229102490000104
其中:
Figure RE-GDA0002229102490000105
Figure RE-GDA0002229102490000106
Figure RE-GDA0002229102490000107
为了获得更好的操纵性和稳定性,选取质心侧偏角和横摆角速度为控制输出变量,则有:
z1(t)=β(t)
z2(t)=r(t)
其中,β(t)表示质心侧偏角β随时间t的变化,r(t)表示横摆角速度r随时间t的变化;
令q1=1/vx,q2=1/vx 2,则系统可以描述为:
Figure RE-GDA0002229102490000111
z1(t)=C1x(t)
z2(t)=C2x(t)
其中:
C1=[1 0],C2=[0 1]
Figure RE-GDA0002229102490000112
步骤3:考虑时变速度的LPV控制器设计
假设汽车纵向速度vx是实时可测的,且在范围[vx min vx max]内变化,则 qi(i=1,2)在范围[qi min qi max]内变化,设变量
Figure RE-GDA0002229102490000113
则q的值就会在一个多胞体范围内变化,可以记为:
Figure RE-GDA0002229102490000114
其中:
α1=x(1-y),α2=(1-x)y,α3=xy,α4=(1-x)(1-y),
Figure RE-GDA0002229102490000115
Figure RE-GDA0002229102490000116
考虑控制器结构为u=Kx(t),则系统可以描述为
Figure RE-GDA0002229102490000121
z1(t)=C1x(t)
z2(t)=C2x(t)
其中
Figure RE-GDA0002229102490000122
K为状态反馈增益矩阵。
Figure RE-GDA0002229102490000123
为x(t)的一阶导数,外部横摆力矩ΔMz=u。
我们假设汽车纵向速度是时变的,上面所建立的系统模型就变为连续时间的多胞体LPV系统。分别选取H和GH2范数来表征输出z1和z2,由于实际情况中前轮转向角δ是有一定限度的,因此考虑δ∈L2,L2表示前轮转向角的取值范围。 ||z1||2,||δ||2分别表示z1和δ的二范数,||z2||表示z2的无穷范数。||T1||表示系统的 H增益,
Figure RE-GDA0002229102490000124
表示系统的GH2增益(广义H2增益)。
Figure RE-GDA0002229102490000125
状态反馈增益矩阵K的取值可由如下方法确定:针对上述LPV系统,若存在矩阵 X>0和Wi>0,i=1,2,3,4.使得以下LMI成立,则存在LPV控制器可使得系统二次稳定,且满足
Figure RE-GDA0002229102490000126
Figure RE-GDA0002229102490000127
Figure RE-GDA0002229102490000128
Figure RE-GDA0002229102490000129
通过求解以上矩阵不等式,可以求得X和Wi,然后可求得
Figure RE-GDA00022291024900001210
其中Ki=WiX-1。然后通过下式计算汽车的横摆力矩ΔMz
ΔMz=Kx(t)
步骤4:转矩分配
四轮驱动汽车的转矩分配主要有两种方法,第一种是按设置好的既定规则分配,第二种是基于目标优化的方法分配。本发明采用基于目标优化的方法进行转矩分配。
如图3所示,对车辆运动模型进行分析,可以得到:
ΔMz=Fx1(-cosδ·ls+sinδ·lf)+Fx2(cosδ·ls+sinδ·lf)-Fx3·ls+Fx4·ls
其中Fxi,i=1,2,3,4为四个车轮的纵向轮胎力,δ为前轮转向角,ls为车辆横轴长度的一半,其余参数与图2中所示含义相同。
转矩分配策略还需满足纵向驱动合力Fx和四个车轮的纵向轮胎力Fxi之间的关系如下:
Fx=Fx1·cosδ+Fx2·cosδ+Fx3+Fx4
对上面的式子进行整理,可得y=Gξ,其中:
y=[ΔMz Fx]T,ξ=[Fx1 Fx2 Fx3 Fx4]T
Figure RE-GDA0002229102490000131
我们假设要使车辆稳定所需要的外部横摆力矩为ΔMdes,Fx=max,那么y的期望值为yd=[ΔMdes Fx]T
转矩分配优化目标函数由两部分组成:
(3)考虑转矩的分配误差,我们得到目标函数的第一部分
J1=(Gξ-yd)TQ1(Gξ-yd),其中Q1=diag{q11 q12}为权重因子。
(4)考虑轮胎的轮胎利用率为
Figure RE-GDA0002229102490000132
其中i=1,2,3,4.μi为第i个轮胎与路面之间的摩擦系数,Fzi为每个轮胎的垂向载荷,分别为:
Figure RE-GDA0002229102490000141
Figure RE-GDA0002229102490000142
其中g为重力加速度,ax和ay分别为车辆的纵向加速度和横向加速度,h 为车辆的重心高度,L=lf+lr。然后可得目标函数的第二部分为:
Figure RE-GDA0002229102490000143
其中
Figure RE-GDA0002229102490000144
最终的目标函数为J=J1+J2
通过最优化此目标函数可以得到理想的转矩。
然后考虑车轮的纵向力与转矩之间的关系为:
Ti=Fxi·R,i=1,2,3,4.R为车轮转动半径,Ti为每个轮内电机产生的转矩,在实际运行中,假设每个轮内电机产生的最大转矩为Tu,最小为Tl,即Tl≤Ti≤Tu,那么可以得到轮胎力的约束为:
Figure RE-GDA0002229102490000145
综上可将转矩分配问题描述为优化问题如下:
Figure RE-GDA0002229102490000146
上式可利用matlab软件进行优化求解,得到ξ,也即得到了四个轮内电机应输出的力的大小,通过电子系统可控制电机按需求输出。
注:对于目标函数J,对其求导可得:
Figure RE-GDA0002229102490000151
由于GTw1G+w2>0,因此在给定ξ的情况下,目标函数存在最小值,即上述优化问题有解,这表明我们的力矩分配策略是可行的。
以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然对前述各实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (8)

1.一种基于时变速度的四轮驱动汽车横向稳定性控制方法,其特征在于,所述控制方法包括以下步骤:
T1采集汽车的固有参数和实时参数;
T2利用T1的参数,建立汽车控制模型系统;
T3对T2系统考虑时变速度设计LPV控制器;
T4根据T3结果对四轮驱动汽车的转矩分配;
所述T4中,四轮驱动汽车的转矩分配采用基于目标优化的方法进行转矩分配;
对车辆运动模型进行分析,得到:
ΔMz=Fx1(-cosδ·ls+sinδ·lf)+Fx2(cosδ·ls+sinδ·lf)-Fx3·ls+Fx4·ls
其中Fxi,i=1,2,3,4为四个车轮的纵向轮胎力,δ为前轮转向角,ls为车辆横轴长度的一半,ΔMz是外部提供的横摆力矩,lf是前轴到汽车重心的距离;
转矩分配策略还需满足纵向驱动合力Fx和四个车轮的纵向轮胎力Fxi之间的关系如下:
Fx=Fx1·cosδ+Fx2·cosδ+Fx3+Fx4
对上面的式子进行整理,得y=Gξ,其中:
y=[ΔMz Fx]T,ξ=[Fx1 Fx2 Fx3 Fx4]T
Figure FDA0002723489470000011
假设要使车辆稳定所需要的外部横摆力矩为ΔMdes,Fx=max,那么y的期望值为yd=[ΔMdes Fx]T
2.根据权利要求1所述的基于时变速度的四轮驱动汽车横向稳定性控制方法,其特征在于,所述T1中,
所述固有参数包括汽车总质量m,转动惯量Iz,前、后轴到汽车重心的距离lf、lr,以及前、后轮的侧偏刚度Cf、Cr;这些参数是汽车固有的,提前采集并存储;
所述实时参数包括前轮转向角δ,汽车纵向速度vx,横摆角速度r,质心侧偏角β;
所述实时参数中,前轮转向角通过传感器采集方向盘转角,然后通过转向传动机构的参数来计算获得;横摆角速度通过陀螺仪测算获得;车身纵向速度和质心侧偏角通过状态参数估计器结合方向盘转角、横摆角速度、车身各方向加速度计算得到,车身各方向加速度由陀螺仪得到。
3.根据权利要求1所述的基于时变速度的四轮驱动汽车横向稳定性控制方法,其特征在于,所述T2中,首先建立汽车二自由度模型,CG表示汽车的重心,m为汽车总质量,β为汽车质心侧偏角,r为横摆角速度,
Figure FDA0002723489470000025
为横摆角加速度,Fyf和Fyr分别为前、后轮的侧向力,vy和vx分别是汽车横向和纵向速度,
Figure FDA0002723489470000026
为汽车横向加速度,Iz为车身转动惯量,ΔMz是外部提供的横摆力矩,通过转向系统和制动系统提供,lf、lr分别是前后轴到汽车重心的距离;
根据牛顿力学定律,得到以下方程:
Figure FDA0002723489470000021
Figure FDA0002723489470000022
在上面的方程中:
Fyf=Cfαf,Fyr=Crαr
Figure FDA0002723489470000023
其中,Cf、Cr分别为前后轮胎的侧偏刚度,αfr分别为前、后轮侧偏角。
4.根据权利要求3所述的基于时变速度的四轮驱动汽车横向稳定性控制方法,其特征在于,选取质心侧偏角β和横摆角速度r作为控制模型的状态变量,得到如下汽车横向运动模型:
Figure FDA0002723489470000024
其中:
Figure FDA0002723489470000031
u(t)=ΔMz,w(t)=δ
Figure FDA0002723489470000032
Figure FDA0002723489470000033
为了获得更好的操纵性和稳定性,选取质心侧偏角和横摆角速度为控制输出变量,则有:
z1(t)=β(t)
z2(t)=r(t)
其中,β(t)表示质心侧偏角β随时间t的变化,r(t)表示横摆角速度r随时间t的变化;
令q1=1/vx,q2=1/vx 2,则系统描述为:
Figure FDA0002723489470000034
z1(t)=C1x(t)
z2(t)=C2x(t)
其中:
C1=[1 0],C2=[0 1]
Figure FDA0002723489470000035
5.根据权利要求1所述的基于时变速度的四轮驱动汽车横向稳定性控制方法,其特征在于,所述T3中,假设汽车纵向速度vx是实时测的,且在范围[vxmin vxmax]内变化,则qi(i=1,2)在范围[qimin qimax]内变化,设变量
Figure FDA0002723489470000041
则q的值就会在一个多胞体范围内变化,记为:
Figure FDA0002723489470000042
其中:
α1=x(1-y),α2=(1-x)y,α3=xy,α4=(1-x)(1-y),
Figure FDA0002723489470000043
Figure FDA0002723489470000044
考虑控制器结构为u=Kx(t),则系统描述为
Figure FDA0002723489470000045
z1(t)=C1x(t)
z2(t)=C2x(t)
其中
Figure FDA0002723489470000046
K为状态反馈增益矩阵;
Figure FDA00027234894700000410
为x(t)的一阶导数,外部横摆力矩ΔMz=u。
6.根据权利要求5所述的基于时变速度的四轮驱动汽车横向稳定性控制方法,其特征在于,假设汽车纵向速度是时变的,建立的系统模型就变为连续时间的多胞体LPV系统;
分别选取H和GH2范数来表征输出z1和z2,由于实际情况中前轮转向角δ是有一定限度的,因此考虑δ∈L2,L2表示前轮转向角的取值范围;||z1||2,||δ||2分别表示z1和δ的二范数,||z2||表示z2的无穷范数;||T1||表示系统的H增益,
Figure FDA0002723489470000047
表示系统的GH2增益(广义H2增益);
Figure FDA0002723489470000048
状态反馈增益矩阵K的取值由如下方法确定:针对上述LPV系统,若存在矩阵X>0和Wi>0,i=1,2,3,4,使得以下LMI成立,则存在LPV控制器使得系统二次稳定,且满足
Figure FDA0002723489470000049
Figure FDA0002723489470000051
Figure FDA0002723489470000052
Figure FDA0002723489470000053
通过求解以上矩阵不等式,求得X和Wi,然后求得
Figure FDA0002723489470000054
其中Ki=WiX-1;然后通过下式计算汽车的横摆力矩ΔMz
ΔMz=Kx(t)。
7.根据权利要求1所述的基于时变速度的四轮驱动汽车横向稳定性控制方法,其特征在于,转矩分配优化目标函数由两部分组成:
(1)考虑转矩的分配误差,得到目标函数的第一部分
J1=(Gξ-yd)TQ1(Gξ-yd),其中Q1=diag{q11 q12}为权重因子;
(2)考虑轮胎的轮胎利用率为
Figure FDA0002723489470000055
其中i=1,2,3,4.μi为第i个轮胎与路面之间的摩擦系数,Fzi为每个轮胎的垂向载荷,分别为:
Figure FDA0002723489470000056
Figure FDA0002723489470000057
其中g为重力加速度,ax和ay分别为车辆的纵向加速度和横向加速度,h为车辆的重心高度,L=lf+lr;然后得目标函数的第二部分为:
Figure FDA0002723489470000058
其中
Figure FDA0002723489470000059
最终的目标函数为J=J1+J2;通过最优化此目标函数得到理想的转矩。
8.根据权利要求1所述的基于时变速度的四轮驱动汽车横向稳定性控制方法,其特征在于,考虑车轮的纵向力与转矩之间的关系为:
Ti=Fxi·R,i=1,2,3,4.R为车轮转动半径,Ti为每个轮内电机产生的转矩,在实际运行中,假设每个轮内电机产生的最大转矩为Tu,最小为Tl,即Tl≤Ti≤Tu,那么得到轮胎力的约束为:
Figure FDA0002723489470000061
综上将转矩分配问题描述为优化问题如下:
Figure FDA0002723489470000062
上式利用matlab软件进行优化求解,得到ξ,也即得到了四个轮内机应输出的力的大小,通过电子系统控制电机按需求输出。
CN201910137939.5A 2019-02-25 2019-02-25 一种基于时变速度的四轮驱动汽车横向稳定性控制方法 Active CN110509915B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910137939.5A CN110509915B (zh) 2019-02-25 2019-02-25 一种基于时变速度的四轮驱动汽车横向稳定性控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910137939.5A CN110509915B (zh) 2019-02-25 2019-02-25 一种基于时变速度的四轮驱动汽车横向稳定性控制方法

Publications (2)

Publication Number Publication Date
CN110509915A CN110509915A (zh) 2019-11-29
CN110509915B true CN110509915B (zh) 2021-03-09

Family

ID=68622433

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910137939.5A Active CN110509915B (zh) 2019-02-25 2019-02-25 一种基于时变速度的四轮驱动汽车横向稳定性控制方法

Country Status (1)

Country Link
CN (1) CN110509915B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110979303B (zh) * 2019-12-24 2021-06-11 大连理工大学 行驶车辆决策附加横摆力矩的方法
CN111267835B (zh) * 2020-03-26 2021-04-27 桂林电子科技大学 基于模型预测算法的四轮独立驱动汽车稳定性控制方法
CN112590483B (zh) * 2021-01-05 2023-06-16 广东工业大学 一种基于观测器的汽车横向稳定性与主动悬架多目标集成控制方法
CN113176733B (zh) * 2021-04-27 2023-06-16 广东工业大学 一种基于切换控制的自主车辆路径跟踪与稳定性控制方法
CN113635780B (zh) * 2021-08-17 2023-09-08 长春工业大学 四轮驱动电动汽车驱动节能与操稳底盘一体化控制方法

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
DE10132440A1 (de) * 2001-07-04 2003-01-23 Bosch Gmbh Robert System und Verfahren zum Überwachen des Fahrverhaltens eines Fahrzeugs
FR2922846B1 (fr) * 2007-10-30 2009-11-20 Renault Sas Systeme de commande du braquage et du freinage d'un vehicule automobile.
CN103303157B (zh) * 2013-06-19 2015-04-08 电子科技大学 一种四轮驱动电动汽车转矩分配方法
DE102014223646A1 (de) * 2014-11-19 2016-05-19 Robert Bosch Gmbh Steuergerät für eine Momentenberechnung für einen Verbrennungsmotor
CN108099900B (zh) * 2017-12-18 2019-09-03 长春工业大学 一种极限工况下保持汽车横向稳定的四轮转向控制方法
CN109017760B (zh) * 2018-07-31 2022-04-12 大连民族大学 车辆期望轨迹跟踪方法、装置及滚动时域优化算法
CN109367532B (zh) * 2018-08-20 2020-04-17 广东工业大学 一种基于速度依赖的汽车横向稳定性控制方法

Also Published As

Publication number Publication date
CN110509915A (zh) 2019-11-29

Similar Documents

Publication Publication Date Title
CN110509915B (zh) 一种基于时变速度的四轮驱动汽车横向稳定性控制方法
CN110466359B (zh) 轮毂四驱纯电动汽车扭矩矢量控制系统及控制方法
Zhai et al. Electronic stability control based on motor driving and braking torque distribution for a four in-wheel motor drive electric vehicle
CN110254405B (zh) 一种面向自动驾驶与智能辅助驾驶的汽车线控制动控制系统及其控制方法
CN104786804B (zh) 一种车辆及其轮边驱动系统和轮边驱动扭矩分配方法
CN110606079B (zh) 一种分层控制的车辆防侧翻方法及多轴分布式驱动车辆
CN109606133A (zh) 基于双层控制的分布式驱动电动汽车转矩矢量控制方法
CN111267835B (zh) 基于模型预测算法的四轮独立驱动汽车稳定性控制方法
Li et al. Design of a vehicle lateral stability control system via a fuzzy logic control approach
US20110307129A1 (en) Vehicle steerability and stability control via independent wheel torque control
Huh et al. Handling and driving characteristics for six-wheeled vehicles
CN113353081B (zh) 一种四驱车的前后轴扭矩分配系统和方法
JP4911111B2 (ja) 車両の制御装置
CN105857304A (zh) 基于四轮驱动汽车力矩分配控制系统
US20230339456A1 (en) Method and apparatus for front and rear driving torque distribution of vehicle, and vehicle
CN111152661A (zh) 一种四轮分布式驱动客车电驱动系统失效控制方法
WO2014054148A1 (ja) 車輪制御装置、車両、車輪制御方法
Shen et al. Multi-axle vehicle dynamics stability control algorithm with all independent drive wheel
WO2024012089A1 (zh) 分布式三电机车辆的控制方法、装置、电动车和介质
JP2023543189A (ja) 電気自動車の四駆トルク配分方法、システム及び車両
CN116512934A (zh) 实现三电机四驱电动汽车能耗优化的扭矩分配控制方法
Li et al. Adaptive sliding mode control of lateral stability of four wheel hub electric vehicles
Shuang et al. Skid steering in 4-wheel-drive electric vehicle
CN110239499A (zh) 车辆的控制装置及车辆的控制方法
JPH11151956A (ja) 車両運動制御装置

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