CN110789504A - 基于电磁制动迟滞补偿的非线性复合abs控制方法 - Google Patents

基于电磁制动迟滞补偿的非线性复合abs控制方法 Download PDF

Info

Publication number
CN110789504A
CN110789504A CN201910993608.1A CN201910993608A CN110789504A CN 110789504 A CN110789504 A CN 110789504A CN 201910993608 A CN201910993608 A CN 201910993608A CN 110789504 A CN110789504 A CN 110789504A
Authority
CN
China
Prior art keywords
braking torque
mathematical model
brake
braking
tire
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
Application number
CN201910993608.1A
Other languages
English (en)
Other versions
CN110789504B (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.)
Dragon Totem Technology Hefei Co ltd
Harbin Meiqirui Technology Information Consulting Service Co ltd
Original Assignee
Jiangsu University
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 Jiangsu University filed Critical Jiangsu University
Priority to CN201910993608.1A priority Critical patent/CN110789504B/zh
Publication of CN110789504A publication Critical patent/CN110789504A/zh
Application granted granted Critical
Publication of CN110789504B publication Critical patent/CN110789504B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B60VEHICLES IN GENERAL
    • B60TVEHICLE BRAKE CONTROL SYSTEMS OR PARTS THEREOF; BRAKE CONTROL SYSTEMS OR PARTS THEREOF, IN GENERAL; ARRANGEMENT OF BRAKING ELEMENTS ON VEHICLES IN GENERAL; PORTABLE DEVICES FOR PREVENTING UNWANTED MOVEMENT OF VEHICLES; VEHICLE MODIFICATIONS TO FACILITATE COOLING OF BRAKES
    • B60T8/00Arrangements for adjusting wheel-braking force to meet varying vehicular or ground-surface conditions, e.g. limiting or varying distribution of braking force
    • B60T8/17Using electrical or electronic regulation means to control braking
    • B60T8/176Brake regulation specially adapted to prevent excessive wheel slip during vehicle deceleration, e.g. ABS
    • B60T8/1761Brake regulation specially adapted to prevent excessive wheel slip during vehicle deceleration, e.g. ABS responsive to wheel or brake dynamics, e.g. wheel slip, wheel acceleration or rate of change of brake fluid pressure
    • 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
    • B60W2050/0028Mathematical models, e.g. for simulation
    • B60W2050/0031Mathematical model of the vehicle
    • 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
    • B60W2050/0028Mathematical models, e.g. for simulation
    • B60W2050/0037Mathematical models of vehicle sub-units

Landscapes

  • Engineering & Computer Science (AREA)
  • Transportation (AREA)
  • Mechanical Engineering (AREA)
  • Physics & Mathematics (AREA)
  • Fluid Mechanics (AREA)
  • Automation & Control Theory (AREA)
  • Human Computer Interaction (AREA)
  • Regulating Braking Force (AREA)

Abstract

本发明公开了一种基于电磁制动迟滞补偿的非线性复合ABS控制方法,包括建立车辆动力学数学模型、轮胎数学模型、电磁制动器数学模型、电子液压制动器数学模型和参考滑移率数学模型;根据电子液压制动器数学模型和电磁制动器数学模型设计执行模块1和执行模块2,根据车辆动力学数学模型、轮胎数学模型和参考滑移率数学模型计算模块1计算出所需的制动力矩,执行模块1提供所需制动力矩;计算模块2根据电子液压制动器计算出所需要的迟滞补偿制动力矩,执行模块2提供所需的迟滞补偿制动力矩。有益效果:不仅减少了力矩传感器的成本,能够有效提高汽车ABS的控制精度和响应速度,使滑移率迅速收敛到目标值,保证汽车的制动安全。

Description

基于电磁制动迟滞补偿的非线性复合ABS控制方法
技术领域
本发明涉及汽车防抱死制动***的控制方法,特别涉及一种基于电磁制动迟滞补偿的非线性复合ABS控制方法,属于汽车制动***领域。
背景技术
防抱死制动***(ABS)是车辆最基本的安全控制,通过机电液***自动控制制动器摩擦力的大小,使车轮不被抱死并保持在最佳滑移率附近。但是传统制动***受限于其结构和工作原理难以满足汽车电动化和智能化等对制动***提出的许多新的功能需求,难以实现人机制动解耦和快速精确的制动等。电子液压制动能够满足上述要求,在汽车未来制动***具有广阔的应用前景,纵观制动***技术发展历程,从ABS到ESC的推出,一般都是渐进功能式的叠加,这将导致功能冗余和处理器的浪费。因此,利用电子液压制动***实现ABS的控制具有十分重要的意义。
但是,电子液压制动***根据现有的最优控制技术进行ABS的控制,效果并不能实现较好的效果,甚至低于传统的ABS控制,最重要的一个原因之一就是电子液压制动***自身存在的迟滞特性将会影响到最优控制。由于执行器特性的原因,最优控制技术在ABS控制领域受到制约,因此,如何改善制动***执行器的响应特性,是各种先进的最优控制技术应用到 ABS的关键。
中国专利CN201410225506.2一种电动车辆的复合制动***及其复合制动方法,提出的复合制动方法,选择优先进行再生制动力矩,当再生制动力矩不足,用摩擦制动力矩进行补偿,尽管能够有效的回收能量,但忽略了两种制动***存在明显的响应差异,没有考虑两者之间的迟滞特性将会导致制动***的震荡。
中国专利CN201811349994.2一种电液复合制动***及其控制方法,是将实际液压制动力与理想液压制动力的差值叠加到电制动力矩的目标输入中,来补偿液压制动力由于滞后产生的差值。上述技术属于反馈控制,需要较高的传感器精度,且测量得到实际制动力矩需要一定时间,再经过差值反馈给电机控制器进行补偿。经测量计算带来的时滞影响,制约了制动***的响应特性的改善。
发明内容
发明目的:针对现有技术中存在不足,本发明提供了一种基于电磁制动迟滞补偿的非线性复合ABS控制方法。
本发明是通过利用电磁制动力矩补偿电子液压制动***制动力矩的迟滞特性,改善ABS 的响应特性,提高制动性能。通过积分反馈非线性最优控制器求解所需制动力矩,利用电子液压制动***进行制动力矩跟踪控制,利用电磁制动进行制动力矩的迟滞补偿。
技术方案:一种基于电磁制动迟滞补偿的非线性复合ABS控制方法,主要包含如下步骤:
步骤一、建立数学模型:分别建立车辆动力学数学模型、轮胎数学模型、电磁制动器数学模型、电子液压制动器数学模型和参考滑移率数学模型;
步骤二、根据电子液压制动器数学模型和电磁制动器数学模型设计执行模块1和执行模块2,根据车辆动力学数学模型、轮胎数学模型和参考滑移率数学模型设计计算模块1,计算模块1计算出所需的制动力矩,执行模块1提供所需制动力矩;计算模块2根据电子液压制动器一阶惯性常数和动态特性计算出所需要的迟滞补偿制动力矩,执行模块2提供所需的迟滞补偿制动力矩。
所述车辆动力学数学模型为:
Figure RE-GDA0002349571570000021
其中R为车轮半径,为车轮总惯性矩,为车辆纵向加速度,
Figure RE-GDA0002349571570000024
为车轮角加速度,Tb为制动力矩,Fx为轮胎纵向力,It为轮胎等效转动惯量;
mt为车辆的总质量的四分之一:
Figure RE-GDA0002349571570000025
其中,mvs为车辆簧载质量,mw为车轮质量;
1/4车模型下的轮胎垂直载荷为:
Figure RE-GDA0002349571570000026
其中l为轴距,hcg为簧载质量质心的高度,FL为制动传递的动荷载,公式(1)代入公式(4)求解非线性代数方程来计算Fz
定义轮胎的滑移率λ:
Figure RE-GDA0002349571570000027
其中,λ表示轮胎滑移率,v表示车辆纵向速度;
对公式(5)关于时间求导:
Figure RE-GDA0002349571570000031
其中,表示
Figure RE-GDA0002349571570000032
表示滑移率的导数,
将公式(1)和(2)代入公式(6)得
所述轮胎数学模型为:
采用非线性Dugoff轮胎模型:
Figure RE-GDA0002349571570000034
其中
Figure RE-GDA0002349571570000035
其中,Ci为轮胎纵向刚度,Cα为轮胎侧向刚度,u为路面摩擦系数,εr为路面附着系数,α为轮胎侧偏角。
所述电磁制动器数学模型为:
电磁制动器施加在制动盘上的制动力矩为 Te=2BIebL (11)
其中,
Figure RE-GDA0002349571570000037
S=ab;
因此,
Figure RE-GDA0002349571570000038
式中,Te为电磁制动力矩,B为穿过制动盘的磁感应强度,a为磁极截面宽度,b为磁极面长度,L为制动盘中心到磁极中心的距离,Ie为在制动盘上产生的电涡流有效值,i为在制动盘上产生的电涡流瞬时值,ρ′为制动盘电阻率,Δh为涡流在制动盘上的集肤深度,μr为相对磁导率,ω为制动盘转动角速度,μ0为真空磁导率,N为铁芯上所扰线圈的匝数,I为电磁制动器通电电流强度,lg为气隙间距,ke为折算系数;
一阶惯性环节来表征电磁制动力矩迟滞特性,即
Figure RE-GDA0002349571570000041
式中:Kbm为电磁制动力矩常数,Tm(s)和Twm(s)分别为电磁制动的实际制动力矩和参考制动力矩,时间常数τ2反映电子液压制动******惯性,s表示频域坐标轴的频率值。
所述电子液压制动器数学模型为:
同样采用一阶惯性环节来表征电子液压制动力矩的变化产生的时间延迟,即
Figure RE-GDA0002349571570000042
式中:Kbh为电子液压制动力矩常数,Th(s)和Thw(s)分别为电子液压制动***实际制动力矩和参考制动力矩,时间常数τ2反映电子液压制动******惯性;
所述参考滑移率数学模型:
其中λopt=0.15表示最优滑移率,a=20为时间常数。对方程两边求拉普拉斯逆变换,解一阶零初始条件微分方程,得
λd(t)=λoptopte-at (16)。
所述步骤二中计算模块1和执行模块1的设计包括以下步骤:
步骤2.1、用泰勒级数展开法预测下一个时间区间λ(t+h)的滑移率非线性响应:
用公式(1)、(7)构建制动工况下的非线性车辆***动力学状态方程,考虑轮胎滑移率作为***输出,写成状态空间形式
Figure RE-GDA0002349571570000044
Figure RE-GDA0002349571570000045
选取输出向量y1为λ,即
y1=x2 (19)
x=[v λ]T为***状态向量,y1为***输出向量,Tb表示控制输入,非线性轮胎模型(8) 已合并到函数f1和f2中;
因此,定义新的状态变量x3为:
Figure RE-GDA0002349571570000051
控制***的目的为控制车轮滑移率x2=λ和其积分x3=∫λdt收敛于参考响应;将状态变量 y=[x2 x3]T作为***的输出,构建优化下一时刻跟踪错误的性能指标:
Figure RE-GDA0002349571570000052
即:
Figure RE-GDA0002349571570000053
w1和w2分别为轮胎滑移率和其积分的加权系数,为实现最高效的跟踪,控制输入Tb的加权项未包含在性能指标中;
在t时刻处的k阶泰勒级数近似为
Figure RE-GDA0002349571570000054
一阶泰勒级数对x2和x3进行泰勒级数展开:
Figure RE-GDA0002349571570000056
对参考滑移率的状态向量进行增广
Figure RE-GDA0002349571570000057
Figure RE-GDA0002349571570000058
步骤2.2以预测跟踪误差最小化为优化目标,计算当前控制Tb(t):
通过将公式(24)-(27)带入到公式(22)得到以控制输入为变量的性能指标函数,根据最优理论,性能指标函数最优化的必要条件为
Figure RE-GDA0002349571570000061
求得
Figure RE-GDA0002349571570000062
其中,e2和e3为当前输出向量的跟踪误差:e2=x2(t)-x2d(t);e3=x3(t)-x3d(t);β为权重比:
Figure RE-GDA0002349571570000063
式(29)是积分反馈的非线性控制***的输出控制。
所述步骤二中基计算模块2和执行模块2的设计方法为:
制动***将滑移率信号和车速信号发送给作为上层控制器的计算模块1,计算模块1求解得到制动力矩Tb,并将制动力矩Tb信号发送到执行模块1,执行模块1实现制动力矩Tb的跟踪控制;作为下层控制器的计算模块2通过采样制动力矩Tb信号,求解得到ABS制动***中液压制动力的变化率,采用一阶泰勒级数前馈补偿控制方法求解电子液压制动***迟滞补偿制动力Tc
Figure RE-GDA0002349571570000064
迟滞补偿控制器将补偿制动力Tc信号发送到执行模块2信号,底层执行模块2实现制动力矩Tc的跟踪控制;最后,将执行模块1产生的制动力矩Th和执行模块2产生的制动力矩Tc相加得到迟滞补偿后的复合制动力矩Tbc作用于***:
Figure RE-GDA0002349571570000065
有益效果:不仅减少了力矩传感器的成本,基于电磁制动迟滞补偿的积分反馈非线性最优复合ABS控制能够有效提高汽车ABS的控制精度和响应速度,使滑移率迅速收敛到目标值,保证汽车的制动安全。
附图说明
图1为本发明1/4车制动模型示意图;
图2为本发明控制原理图;
图3为滑移率和电子液压制动控制的滑移率对比图;
图4为电磁制动迟滞补偿的制动力矩。
具体实施方式
下面结合附图以及具体实施例对本发明作进一步的说明,但本发明的保护范围并不限于此。
1、建立数学模型
1.1、建立车辆模型
本发明将1/4车模型作为研究对象,车辆制动模型见图1。此模型有纵向车速和轮速两个自由度。
车辆的动力学数学模型为:
Figure RE-GDA0002349571570000071
Figure RE-GDA0002349571570000072
其中R为车轮半径,为车轮总惯性矩,
Figure RE-GDA0002349571570000073
为车辆纵向加速度,
Figure RE-GDA0002349571570000074
为车轮角加速度,Tb为制动力矩,Fx为轮胎纵向力。
mt为给出的四分之一车辆的总质量:
其中mvs为车辆簧载质量,mw为车轮质量。
作用在轮胎上的纵向力取决于轮胎的垂直载荷,而垂直载荷由两部分组成:一是由于车辆质量分布而产生的静态载荷,二是由于制动过程中产生的轮胎动载荷。因此,1/4车模型下的轮胎垂直载荷为:
Figure RE-GDA0002349571570000076
其中l为轴距,hcg为簧载质量质心的高度。FL为制动传递的动荷载。在实际应用中,轮胎的垂直载荷是通过测量轮胎的纵向加速度a来计算的。但在仿真研究中,将(1)代入代入 (4)求解非线性代数方程来计算Fz
定义轮胎的滑移率:
Figure RE-GDA0002349571570000077
对公式(5)关于时间求导,得到滑移率关于时间的导数为:
Figure RE-GDA0002349571570000078
将公式(1)和(5)代入公式(6)得
Figure RE-GDA0002349571570000081
方程(1)和(7)构成状态空间运动形式的控制方程。车速v和车轮纵向滑移率λ为状态向量,制动力矩Tb为控制向量。在推导上述方程时,忽略了制动对车身俯仰和横摆的影响,仅考虑无转向时的直线制动工况。
1.2、轮胎模型建立
车轮纵向力Fx被描述为轮胎纵向滑移率的函数。当纵向滑移量较小时,纵向力与滑移量呈线性关系,但随滑移率增加,在某处轮胎纵向力达到最大饱和值,滑移率再增加轮胎纵向力反而有所下降。轮胎力的饱和特性是引起车辆运动非线性特性的主要原因,同样也是造成车辆不安全的主要原因。
为了考虑轮胎力的饱和特性,本发明采用了基于摩擦椭圆思想的非线性Dugoff轮胎模型。在该模型中,轮胎纵向力为:
其中
Figure RE-GDA0002349571570000083
Figure RE-GDA0002349571570000084
其中,Ci为轮胎纵向刚度,Cα为轮胎侧向刚度,u为路面摩擦系数,εr为路面附着系数,α为轮胎侧偏角。
Dugoff轮胎模型描述了轮胎在制动和转向时的轮胎特性。因此,该模型考虑了轮胎纵向力和横向力的相互作用。根据模型中使用的摩擦椭圆概念,当施加制动力时,由于接触面积引起的附加滑移,侧向力逐渐减小。本发明仅考虑没有转向的直线制动工况,对轮胎施加纯纵向力。
1.3、电涡流制动器模型建立
电涡流制动器制动力矩Te计算:电磁制动器施加在制动盘上的制动力矩为
Te=2BIebL (11)
其中,
Figure RE-GDA0002349571570000085
S=ab,S=ab。
因此,
Figure RE-GDA0002349571570000091
式中,Te为电涡流制动力矩,B为穿过制动盘的磁感应强度,a为磁极截面宽度,b为磁极面长度,L为制动盘中心到磁极中心的距离,Ie为在制动盘上产生的电涡流有效值,i为在制动盘上产生的电涡流瞬时值,ρ′为制动盘电阻率,Δh为涡流在制动盘上的集肤深度,μr为相对磁导率,ω为制动盘转动角速度,μ0为真空磁导率,N为铁芯上所扰线圈的匝数,I 为电磁制动器通电电流强度,lg为气隙间距,ke为折算系数,通常取1.5。
并用一阶惯性环节来表征电磁制动力矩的变化产生的时间延迟,即
Figure RE-GDA0002349571570000092
式中:Kbm为电磁制动力矩常数,Tm(s)和Twm(s)分别为电磁制动的实际制动力矩和参考制动力矩,时间常数τ1反映电磁制动***惯性,即一阶***的惯性越小,其响应过程越快;反之惯性越大响应越慢。
1.4、电子液压制动***模型建立
电子液压制动***由于保留原有成熟可靠的液压制动***,作为线控制动***,实现与再生制动协调控制,既能够回收制动能量,又保证制动***的稳定性。同样采用一阶惯性环节来表征电子液压制动力矩的变化产生的时间延迟,即
Figure RE-GDA0002349571570000093
式中:Kbh为电子液压制动力矩常数,Th(s)和Thw(s)分别为电子液压制动***实际制动力矩和参考制动力矩,时间常数τ2反映电子液压制动******惯性。
1.5、参考滑移率模型的建立
基于Dugoff轮胎模型,为在参考模型中包含车轮滑移率的瞬态响应,避免较大跟踪误差和制动初期制动力矩突增,建立车轮滑移率的参考模型
Figure RE-GDA0002349571570000094
其中λopt=0.15表示最优滑移率,a=20为时间常数。其实,式(15)假设期望滑移的阶跃响应为一阶***。对方程两边求拉普拉斯逆变换,解一阶零初始条件微分方程,得到
λd(t)=λoptopte-at (16)
式(16)描述了时域内的车轮滑移率参考模型,下节将设计积分反馈非线性控制器对参考滑移率进行跟踪。
2、电磁制动迟滞补偿的积分反馈非线性复合ABS控制
步骤2.1设计计算模块1
用(1)、(7)构建制动工况下的非线性车辆***动力学状态方程,考虑轮胎滑移率作为***输出,写成状态空间形式
Figure RE-GDA0002349571570000101
Figure RE-GDA0002349571570000102
y1=x2 (19)
x=[v λ]T为***状态向量,y1为***输出向量,Tb表示控制输入。非线性轮胎模(8) 已合并到函数f1和f2中。
本发明计算模块1的设计思路是基于一种优化的非线性控制律。首先,用泰勒级数展开法预测下一个时间区间λ(t+h)的滑移率非线性响应,预测周期h类似于预测控制中的预测时域;然后,在连续区间内以预测跟踪误差最小化为优化目标,找出当前控制Tb(t);最后,为进一步提高滑移率控制器的鲁棒性,采用积分反馈技术,将滑移率的积分附加到状态方程中。
因此,定义新的状态变量x3为:
Figure RE-GDA0002349571570000103
控制***的目的为控制车轮滑移率x2=λ和其积分x3=∫λdt收敛于参考响应。将状态变量 y=[x2 x3]T作为***的输出,构建优化下一时刻跟踪错误的性能指标:
Figure RE-GDA0002349571570000104
即:
w1和w2分别为轮胎滑移率和其积分的加权系数,为实现最高效的跟踪,控制输入Tb的加权项未包含在性能指标中。
在t时刻处的k阶泰勒级数近似为
Figure RE-GDA0002349571570000111
泰勒级数预测的关键问题在于阶数的选择,阶数越高近似程度越高,但控制***能量将会增加,而较低的阶数导致预测的误差增大。因此,控制阶数被认为是一个设计参数,在性能和输入能量需求之间做出妥协。泰勒级数预测的充分条件为不低于预测向量的阶数,一阶泰勒级数对于x2的展开是充分的,而对于x3的展开至少需要二阶泰勒级数。
上式用到了公式(19)中
Figure RE-GDA0002349571570000114
的关系。
同理,可对参考滑移率的状态向量进行增广
Figure RE-GDA0002349571570000115
Figure RE-GDA0002349571570000116
因此,通过将公式(24)-(27)带入到(22)得到以控制输入为变量的性能指标函数,根据最优理论,性能指标函数最优化的必要条件为
求得
其中,e2和e3为当前输出向量的跟踪误差:e2=x2(t)-x2d(t);e3=x3(t)-x3d(t);β为权重比:
Figure RE-GDA0002349571570000119
式(29)是积分反馈的非线性控制***的输出控制。
步骤2.2设计计算模块2
迟滞特性将会影响控制***的性能,所以必须对其进行补偿,才能实现稳定精密的控制。迟滞补偿控制器的思路为通过积分反馈非线性控制器求解所需制动力矩,利用电子液压制动***进行制动力矩跟踪控制,利用电磁制动进行制动力矩的迟滞补偿。
迟滞补偿控制器的原理图如图2所示,制动***将滑移率信号和车速信号发送给计算模块1,该上层控制器计算模块1求解得到目标制动力矩Tb,并将目标制动力矩Tb信号发送到执行模块1,执行模块1实现目标制动力矩Tb的跟踪控制。但是由于执行模块1存在明显的迟滞特性,无法实现跟踪高频信号的跟踪控制。计算模块2通过采样执行模块1的信号,求解得到ABS制动***中液压制动力的变化率,采用一阶泰勒级数前馈补偿控制方法求解电子液压制动***迟滞补偿制动力Tc。
计算模块1将补偿制动力Tc信号发送到执行模块2,底层执行模块2实现制动力矩Tc的跟踪控制。最后,将执行模块1产生的制动力矩Th和执行模块2产生的制动力矩Tc相加得到迟滞补偿后的复合制动力矩Tbc,并作用于制动***。
其中,计算模块2迟滞补偿控制器利用微分概念,假设时间常数τ为无穷小,电子液压制动力矩下一时刻的一阶惯性迟滞环节补偿量等于迟滞惯性力的斜率与时间常数τ之积,因此τ时刻之后经补偿的电子液压制动力矩为
液压制动力矩的一阶惯性环节为自身固有特性,在力的跟踪控制时只能实现Th部分,自身无法实现
Figure RE-GDA0002349571570000122
的迟滞补偿。但电磁制动具有响应速度快的特点,响应速度达到毫秒级别,能够对液压制动力的迟滞特性进行有效补偿。因此液压制动力矩补偿量通过电磁制动Tc来完成,即
实施例:
为验证上述控制方法的有效性,本节分别对高附着系数路面的ABS制动工况进行仿真,对采用基于电磁制动迟滞补偿的积分反馈非线性复合制动与电子液压制动ABS控制策略进行对比分析。实际应用时所需的参数:u为路面摩擦系数选为0.8,R=0.326m,v=30m/s, Ci=30000,Cα=50000N/rad,mw=40kg,It=1.7kgm2,mvs=415kg,a=20,h=0.01,β=20,τ1=150ms,τ2=5ms,εr=0.015s/m,hcg=0.5m。
图3为基于电磁制动迟滞补偿的积分反馈非线性复合ABS控制的滑移率和电子液压制动控制的滑移率对比图。通过对比发现,基于电磁制动迟滞补偿的积分反馈非线性复合ABS控制的滑移率能够更快的收敛到目标值,能够有效提高电动汽车ABS的控制精度和响应速度。
图4为电磁制动迟滞补偿的制动力矩。充分利用电磁制动力矩的响应特性来改善制动***的响应特性。
所述实施例为本发明的优选的实施方式,但本发明并不限于上述实施方式,在不背离本发明的实质内容的情况下,本领域技术人员能够做出的任何显而易见的改进、替换或变型均属于本发明的保护范围。

Claims (8)

1.一种基于电磁制动迟滞补偿的非线性复合ABS控制方法,其特征在于,主要包含如下步骤:
步骤一、建立数学模型:分别建立车辆数学模型、轮胎数学模型、电磁制动器数学模型、电子液压制动器数学模型和参考滑移率数学模型;
步骤二、根据电子液压制动器数学模型和电磁制动器数学模型设计执行模块1和执行模块2,根据车辆数学模型、轮胎数学模型和参考滑移率数学模型设计计算模块1,计算模块1计算出所需的制动力矩,执行模块1提供所需制动力矩;计算模块2根据电子液压制动器一阶惯性常数和动态特性计算出所需要的迟滞补偿制动力矩,执行模块2提供所需的迟滞补偿制动力矩。
2.根据权利要求1所述的基于电磁制动迟滞补偿的非线性复合ABS控制方法,其特征在于:所述车辆数学模型为:
Figure FDA0002239067690000012
其中R为车轮半径,为车轮总惯性矩,
Figure FDA0002239067690000013
为车辆纵向加速度,为车轮角加速度,Tb为制动力矩,Fx为轮胎纵向力,It为轮胎等效转动惯量;
mt为车辆的总质量的四分之一:
Figure FDA0002239067690000015
其中,mvs为车辆簧载质量,mw为车轮质量;
1/4车模型下的轮胎垂直载荷为:
其中l为轴距,hcg为簧载质量质心的高度,FL为制动传递的动荷载,公式(1)代入公式(4)求解非线性代数方程来计算Fz
定义轮胎的滑移率λ:
Figure FDA0002239067690000017
其中,λ表示轮胎滑移率,v表示车辆纵向速度;
对公式(5)关于时间求导:
其中,表示
Figure FDA0002239067690000022
表示滑移率的导数,
将公式(1)和(2)代入公式(6)得
Figure FDA0002239067690000023
3.根据权利要求2所述的基于电磁制动迟滞补偿的非线性复合ABS控制方法,其特征在于,所述轮胎数学模型为:
采用非线性Dugoff轮胎模型:
Figure FDA0002239067690000024
其中
Figure FDA0002239067690000025
Figure FDA0002239067690000026
其中,Ci为轮胎纵向刚度,Cα为轮胎侧向刚度,u为路面摩擦系数,εr为路面附着系数,α为轮胎侧偏角。
4.根据权利要求3所述的基于电磁制动迟滞补偿的非线性复合ABS控制方法,其特征在于,所述电磁制动器数学模型为:
电磁制动器施加在制动盘上的制动力矩为
Te=2BIebL (11)
其中,
Figure FDA0002239067690000027
S=ab;
因此,
Figure FDA0002239067690000028
式中,Te为电磁制动力矩,B为穿过制动盘的磁感应强度,a为磁极截面宽度,b为磁极面长度,L为制动盘中心到磁极中心的距离,Ie为在制动盘上产生的电涡流有效值,i为在制动盘上产生的电涡流瞬时值,ρ′为制动盘电阻率,Δh为涡流在制动盘上的集肤深度,μr为相对磁导率,ω为制动盘转动角速度,μ0为真空磁导率,N为铁芯上所扰线圈的匝数,I为电磁制动器通电电流强度,lg为气隙间距,ke为折算系数;
一阶惯性环节来表征电磁制动力矩迟滞特性,即
式中:Kbm为电磁制动力矩常数,Tm(s)和Twm(s)分别为电磁制动的实际制动力矩和参考制动力矩,时间常数τ2反映电子液压制动******惯性,s表示频域坐标轴的频率值。
5.根据权利要求4所述的基于电磁制动迟滞补偿的非线性复合ABS控制方法,其特征在于,所述电子液压制动器数学模型为:
同样采用一阶惯性环节来表征电子液压制动力矩的变化产生的时间延迟,即
Figure FDA0002239067690000032
式中:Kbh为电子液压制动力矩常数,Th(s)和Thw(s)分别为电子液压制动***实际制动力矩和参考制动力矩,时间常数τ2反映电子液压制动******惯性。
6.根据权利要求5所述的基于电磁制动迟滞补偿的非线性复合ABS控制方法,其特征在于,所述参考滑移率数学模型:
Figure FDA0002239067690000033
其中λopt=0.15表示最优滑移率,a=20为时间常数。对方程两边求拉普拉斯逆变换,解一阶零初始条件微分方程,得
λd(t)=λoptopte-at(16)。
7.根据权利要求6所述的基于电磁制动迟滞补偿的非线性复合ABS控制方法,其特征在于,所述步骤二中计算模块1和执行模块1的设计包括以下步骤:
用泰勒级数展开法预测下一个时间区间λ(t+h)的滑移率非线性响应:
用公式(1)、(7)构建制动工况下的非线性车辆***动力学状态方程,考虑轮胎滑移率作为***输出,写成状态空间形式
Figure FDA0002239067690000041
Figure FDA0002239067690000042
选取输出向量y1为λ,即
y1=x2(19)
x=[v λ]T为***状态向量,y1为***输出向量,Tb表示控制输入,非线性轮胎模型(8)已合并到函数f1和f2中;
因此,定义新的状态变量x3为:
Figure FDA0002239067690000043
控制***的目的为控制车轮滑移率x2=λ和其积分x3=∫λdt收敛于参考响应;将状态变量y=[x2 x3]T作为***的输出,构建优化下一时刻跟踪错误的性能指标:
Figure FDA0002239067690000044
即:
w1和w2分别为轮胎滑移率和其积分的加权系数,为实现最高效的跟踪,控制输入Tb的加权项未包含在性能指标中;
在t时刻处的k阶泰勒级数近似为
Figure FDA0002239067690000046
一阶泰勒级数对x2和x3进行泰勒级数展开:
Figure FDA0002239067690000047
Figure FDA0002239067690000048
对参考滑移率的状态向量进行增广
Figure FDA0002239067690000051
Figure FDA0002239067690000052
以预测跟踪误差最小化为优化目标,计算当前控制Tb(t):
通过将公式(24)-(27)带入到公式(22)得到以控制输入为变量的性能指标函数,根据最优理论,性能指标函数最优化的必要条件为
求得
Figure FDA0002239067690000054
其中,e2和e3为当前输出向量的跟踪误差:e2=x2(t)-x2d(t);e3=x3(t)-x3d(t);β为权重比:
Figure FDA0002239067690000055
式(29)是积分反馈的非线性控制***的输出控制。
8.根据权利要求7所述的基于电磁制动迟滞补偿的非线性复合ABS控制方法,其特征在于,所述步骤二中计算模块2和执行模块2的设计方法为:
制动***将滑移率信号和车速信号发送给作为上层控制器的计算模块1,计算模块1求解得到制动力矩Tb,并将制动力矩Tb信号发送到执行模块1,执行模块1实现制动力矩Tb的跟踪控制;作为下层控制器的计算模块2通过采样制动力矩Tb信号,求解得到ABS制动***中液压制动力的变化率,采用一阶泰勒级数前馈补偿控制方法求解电子液压制动***迟滞补偿制动力Tc
Figure FDA0002239067690000056
迟滞补偿控制器将补偿制动力Tc信号发送到执行模块2信号,底层执行模块2实现制动力矩Tc的跟踪控制;最后,将执行模块1产生的制动力矩Th和执行模块2产生的制动力矩Tc相加得到迟滞补偿后的复合制动力矩Tbc作用于***:
Figure FDA0002239067690000057
CN201910993608.1A 2019-10-18 2019-10-18 基于电磁制动迟滞补偿的非线性复合abs控制方法 Active CN110789504B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910993608.1A CN110789504B (zh) 2019-10-18 2019-10-18 基于电磁制动迟滞补偿的非线性复合abs控制方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910993608.1A CN110789504B (zh) 2019-10-18 2019-10-18 基于电磁制动迟滞补偿的非线性复合abs控制方法

Publications (2)

Publication Number Publication Date
CN110789504A true CN110789504A (zh) 2020-02-14
CN110789504B CN110789504B (zh) 2021-06-22

Family

ID=69439442

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910993608.1A Active CN110789504B (zh) 2019-10-18 2019-10-18 基于电磁制动迟滞补偿的非线性复合abs控制方法

Country Status (1)

Country Link
CN (1) CN110789504B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112329143A (zh) * 2020-11-04 2021-02-05 浙江天行健智能科技有限公司 一种制动能量回收研发***及其控制方法
CN113581152A (zh) * 2021-07-20 2021-11-02 浙江万安科技股份有限公司 一种线控制动***及一种主缸液压力时滞补偿控制方法
CN114771477A (zh) * 2022-05-27 2022-07-22 烟台大学 一种基于摩檫力随动补偿的汽车防抱死制动方法

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2927065A1 (en) * 2014-04-03 2015-10-07 The Goodyear Tire & Rubber Company Road surface friction and surface type estimation system and method
CN105682980A (zh) * 2013-11-18 2016-06-15 Ntn株式会社 防抱死制动器控制装置
DE102016219727A1 (de) * 2016-10-11 2018-04-12 Robert Bosch Gmbh Hydraulisches Bremssystem für ein Fahrzeug und bistabiles Magnetventil
CN108058601A (zh) * 2017-12-18 2018-05-22 长春工业大学 一种基于线性时变的电动汽车防抱死控制方法
KR20180094646A (ko) * 2017-02-16 2018-08-24 주식회사 만도 전자식 브레이크 시스템 및 그 제어 방법
CN108583554A (zh) * 2018-05-23 2018-09-28 江苏大学 一种基于混杂理论的电磁与摩擦集成制动方法
JP2019016086A (ja) * 2017-07-05 2019-01-31 日立オートモティブシステムズ株式会社 自動車用電子制御装置
CN109501790A (zh) * 2018-10-24 2019-03-22 江苏大学 挂车制动压力延迟补偿与车轮滑移率联合控制***及方法

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105682980A (zh) * 2013-11-18 2016-06-15 Ntn株式会社 防抱死制动器控制装置
EP2927065A1 (en) * 2014-04-03 2015-10-07 The Goodyear Tire & Rubber Company Road surface friction and surface type estimation system and method
DE102016219727A1 (de) * 2016-10-11 2018-04-12 Robert Bosch Gmbh Hydraulisches Bremssystem für ein Fahrzeug und bistabiles Magnetventil
KR20180094646A (ko) * 2017-02-16 2018-08-24 주식회사 만도 전자식 브레이크 시스템 및 그 제어 방법
JP2019016086A (ja) * 2017-07-05 2019-01-31 日立オートモティブシステムズ株式会社 自動車用電子制御装置
CN108058601A (zh) * 2017-12-18 2018-05-22 长春工业大学 一种基于线性时变的电动汽车防抱死控制方法
CN108583554A (zh) * 2018-05-23 2018-09-28 江苏大学 一种基于混杂理论的电磁与摩擦集成制动方法
CN109501790A (zh) * 2018-10-24 2019-03-22 江苏大学 挂车制动压力延迟补偿与车轮滑移率联合控制***及方法

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112329143A (zh) * 2020-11-04 2021-02-05 浙江天行健智能科技有限公司 一种制动能量回收研发***及其控制方法
CN112329143B (zh) * 2020-11-04 2023-09-26 浙江天行健智能科技有限公司 一种制动能量回收研发***及其控制方法
CN113581152A (zh) * 2021-07-20 2021-11-02 浙江万安科技股份有限公司 一种线控制动***及一种主缸液压力时滞补偿控制方法
CN114771477A (zh) * 2022-05-27 2022-07-22 烟台大学 一种基于摩檫力随动补偿的汽车防抱死制动方法
CN114771477B (zh) * 2022-05-27 2024-02-09 烟台大学 一种基于摩檫力随动补偿的汽车防抱死制动方法

Also Published As

Publication number Publication date
CN110789504B (zh) 2021-06-22

Similar Documents

Publication Publication Date Title
CN110789504B (zh) 基于电磁制动迟滞补偿的非线性复合abs控制方法
Wang et al. Coordination control of differential drive assist steering and vehicle stability control for four-wheel-independent-drive EV
CN107719372B (zh) 基于动态控制分配的四驱电动汽车动力学多目标控制***
WO2021077727A1 (zh) 一种电动卡车线控转向***及其网络不确定控制方法
CN108107731B (zh) 一种基于轮胎非线性特性的汽车稳定性控制方法
Yu et al. A global approach to vehicle control: Coordination of four wheel steering and wheel torques
CN109050658A (zh) 基于模型预测控制的汽车主动前轮转向自适应调节方法
CN108107732A (zh) 主动前轮转向和直接横摆力矩联合的汽车稳定性控制方法
Li et al. Integrated longitudinal and lateral vehicle stability control for extreme conditions with safety dynamic requirements analysis
CN111391822A (zh) 一种极限工况下汽车横纵向稳定性协同控制方法
CN108099877A (zh) 一种紧急制动工况下滑移率跟踪控制方法
CN108181811A (zh) 一种基于线性时变的滑移率跟踪控制方法
CN110217227A (zh) 一种适用于冰雪道路工况的转向制动联合避撞控制方法
Lu et al. Global integrated control of vehicle suspension and chassis key subsystems
Hang et al. Path-tracking controller design for a 4WIS and 4WID electric vehicle with steer-by-wire system
Saikia et al. Vehicle stability enhancement using sliding mode based active front steering and direct yaw moment control
Gong et al. Research on electronic brake force distribution and anti-lock brake of vehicle based on direct drive electro hydraulic actuator
CN111027136B (zh) 基于分数阶极值搜索的非线性abs控制方法
Liu et al. Redundancy control of anti-lock braking system based on electro-hydraulic braking system
Yun et al. Brake performance evaluation of ABS with sliding mode controller on a split road with driver model
Bardawil et al. Integrated vehicle lateral stability control with different coordination strategies between active steering and differential braking
Zhang et al. Research on coordinated control of braking mode switching for enhancing switching quality considering communication delay and time-varying model parameters
Kejun et al. Application study of fuzzy PID control with S-function on automotive ABS
Lei et al. Vehicle direct yaw moment control based on tire cornering stiffness estimation
Shi et al. A novel integral terminal sliding mode control of yaw stability for steer-by-wire vehicles

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
TR01 Transfer of patent right

Effective date of registration: 20231219

Address after: No. 3, 3rd Floor, Unit 2, Building 4-7, Hanshui Road Science and Technology Apartment, Nangang Concentration District, Harbin Economic Development Zone, Heilongjiang Province, 150000 RMB

Patentee after: Harbin Meiqirui Technology Information Consulting Service Co.,Ltd.

Address before: 230000 floor 1, building 2, phase I, e-commerce Park, Jinggang Road, Shushan Economic Development Zone, Hefei City, Anhui Province

Patentee before: Dragon totem Technology (Hefei) Co.,Ltd.

Effective date of registration: 20231219

Address after: 230000 floor 1, building 2, phase I, e-commerce Park, Jinggang Road, Shushan Economic Development Zone, Hefei City, Anhui Province

Patentee after: Dragon totem Technology (Hefei) Co.,Ltd.

Address before: Zhenjiang City, Jiangsu Province, 212013 Jingkou District Road No. 301

Patentee before: JIANGSU University

TR01 Transfer of patent right