CN110949344A - 一种分布式驱动电动汽车制动防抱死控制方法及*** - Google Patents

一种分布式驱动电动汽车制动防抱死控制方法及*** Download PDF

Info

Publication number
CN110949344A
CN110949344A CN201911308436.6A CN201911308436A CN110949344A CN 110949344 A CN110949344 A CN 110949344A CN 201911308436 A CN201911308436 A CN 201911308436A CN 110949344 A CN110949344 A CN 110949344A
Authority
CN
China
Prior art keywords
wheel
braking
current
automobile
adhesion coefficient
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
CN201911308436.6A
Other languages
English (en)
Other versions
CN110949344B (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.)
Bitnei Co ltd
Beijing Institute of Technology BIT
Original Assignee
Bitnei Co ltd
Beijing Institute of Technology BIT
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 Bitnei Co ltd, Beijing Institute of Technology BIT filed Critical Bitnei Co ltd
Priority to CN201911308436.6A priority Critical patent/CN110949344B/zh
Publication of CN110949344A publication Critical patent/CN110949344A/zh
Application granted granted Critical
Publication of CN110949344B publication Critical patent/CN110949344B/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/172Determining control parameters used in the regulation, e.g. by calculations involving measured or detected parameters
    • 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
    • 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/32Arrangements for adjusting wheel-braking force to meet varying vehicular or ground-surface conditions, e.g. limiting or varying distribution of braking force responsive to a speed condition, e.g. acceleration or deceleration

Landscapes

  • Engineering & Computer Science (AREA)
  • Transportation (AREA)
  • Mechanical Engineering (AREA)
  • Regulating Braking Force (AREA)

Abstract

本发明公开了一种分布式驱动电动汽车制动防抱死控制方法及***。该方法包括:采用无迹卡尔曼滤波对汽车的当前车速进行滤波;根据汽车的当前车速确定各车轮的当前滑移率;计算汽车各车轮的当前利用附着系数;将当前滑移率和当前利用附着系数输入自适应神经网络模型,得到当前路面的最优滑移率和最大利用附着系数;根据Thb=(μmaxFz‑Freg)r计算各车轮需要液压制动***提供的制动力矩Thb的大小,记为第一制动力矩;根据Teb=TSMC‑Thb计算各车轮需要其轮毂电机提供的制动力矩Teb的大小,记为第二制动力矩;在汽车的制动防抱死过程中,根据第一制动力矩和第二制动力矩,控制液压制动***输出的制动力矩和各车轮轮毂电机输出的制动力矩。本发明提高了紧急制动的效果。

Description

一种分布式驱动电动汽车制动防抱死控制方法及***
技术领域
本发明涉及制动防抱死控制技术领域,特别是涉及一种分布式驱动电动汽车制动防抱死控制方法及***。
背景技术
目前大多数电动汽车在防抱死制动时依然沿用传统汽车的ABS(Anti-lockBraking System)***,即液压防抱死***。其主要包括制动踏板、真空助力器、制动主缸、ABS压力调节模块、制动管路及制动卡钳。除此之外防抱死装置还有气压防抱死***,其与液压防抱死***的区别是工作介质为高压空气。以上两种防抱死***均为成熟的商用产品。
在液压制动***中,由于制动管路建立油压需要一定时间,且由于液体的可压缩性和制动管路内壁非绝对光滑等因素,油压建立过程并非线型增长且会出现波动,对于紧急制动的稳定性会造成影响。液压制动***建立油压的过程如图1所示。除此之外,气压防抱死***体积、重量均较大,多用于载重量较大的商用货车,并不可以用在电动乘用汽车。
电动汽车的电机可以提供一定的再生制动力,且响应迅速,制动力平稳,但是其再生制动能力受很多条件的影响,诸如电机最大制动功率、储能装置最大供电电流等,因此电动汽车在制动时产生的再生制动力常常很有限,不足以达到整车的总制动需求。
发明内容
本发明的目的是提供一种分布式驱动电动汽车制动防抱死控制方法及***,以提高紧急制动的效果。
为实现上述目的,本发明提供了如下方案:
一种分布式驱动电动汽车制动防抱死控制方法,包括:
采用无迹卡尔曼滤波对汽车的当前车速进行滤波;
根据所述汽车的当前车速确定各车轮的当前滑移率;
计算所述汽车各车轮的当前利用附着系数;
将所述当前滑移率和所述当前利用附着系数输入自适应神经网络模型,得到当前路面的最优滑移率和最大利用附着系数;
根据Thb=(μmaxFz-Freg)r计算各车轮需要液压制动***提供的制动力矩Thb的大小,记为第一制动力矩,其中,Freg=αFregmax,Fregmax为当前电机可产生的最大回馈制动力,α为回馈系数,Freg为当前电机可回馈的制动力,μmax为所述最大利用附着系数;
根据Teb=TSMC-Thb计算各车轮需要其轮毂电机提供的制动力矩Teb的大小,记为第二制动力矩,其中,TSMC为车轮的滑模制动总力矩;
在所述汽车的制动防抱死过程中,根据所述第一制动力矩和所述第二制动力矩,控制液压制动***输出的制动力矩和各车轮轮毂电机输出的制动力矩。
可选的,所述方法还包括:根据
Figure BDA0002323830830000021
计算各车轮的滑模制动总力矩,其中,
Figure BDA0002323830830000022
ksat(s)为引入的边界层,φ为边界层的厚度,s=λ-λd,λ为当前滑移率,λd为最优滑移率,Fz为车轮垂向力,μ为当前利用附着系数,r为车轮滚动半径,J为车轮转动惯量,
Figure BDA0002323830830000025
为车轮加速度,k为系数。
可选的,所述采用无迹卡尔曼滤波对汽车的当前车速进行滤波,具体包括:
Figure BDA0002323830830000023
为状态方程,以
Figure BDA0002323830830000024
为观测方程,采用无迹卡尔曼滤波对汽车的当前车速进行滤波,状态方程中的ax、ay、γ、Mz根据动力学模型计算得到,其中,ax为汽车纵向加速度,ay为汽车侧向加速度,γ为汽车横摆角度,ax测量为传感器测量得到的汽车纵向加速度,ay测量为传感器测量得到的汽车侧向加速度,γ测量为传感器测量得到的汽车横摆角度,vx为汽车纵向行驶速度,vy为汽车侧向行驶速度,Mz为绕z轴的旋转力矩,T为传感器信号的采样频率,k为迭代次数。
可选的,所述方法还包括:
获取标准路面的滑移率与利用附着系数关系曲线,记为标准曲线,所述标准路面包括干沥青路面、湿沥青路面、湿鹅卵石路面和冰路面;
获取样本数据,所述样本数据包括滑移率和与其相对应的利用附着系数;
确定所述样本数据与各标准曲线的相似度;
根据所述相似度确定所述样本数据对各标准曲线的权重系数;
根据各所述标准曲线以及各所述权重系数,确定所述样本数据对应的滑移率与利用附着系数关系曲线,记为样本曲线;
确定所述样本曲线上利用附着系数最大的点以及该点所对应的滑移率,分别记为最大利用附着系数和最优滑移率;并将所述最大利用附着系数和所述最优滑移率作为所述样本数据的标签;
根据所述样本数据以及其标签对神经网络进行训练,得到自适应神经网络模型。
本发明还提供了一种分布式驱动电动汽车制动防抱死控制***,包括:
车速滤波模块,用于采用无迹卡尔曼滤波对汽车的当前车速进行滤波;
滑移率确定模块,用于根据所述汽车的当前车速确定各车轮的当前滑移率;
利用附着系数计算模块,用于计算所述汽车各车轮的当前利用附着系数;
最优参数确定模块,用于将所述当前滑移率和所述当前利用附着系数输入自适应神经网络模型,得到当前路面的最优滑移率和最大利用附着系数;
第一制动力矩计算模块,用于根据Thb=(μmaxFz-Freg)r计算各车轮需要液压制动***提供的制动力矩Thb的大小,记为第一制动力矩,其中,Freg=αFregmax,Fregmax为当前电机可产生的最大回馈制动力,α为回馈系数,Freg为当前电机可回馈的制动力,μmax为所述最大利用附着系数;
第二制动力矩计算模块,用于根据Teb=TSMC-Thb计算各车轮需要其轮毂电机提供的制动力矩Teb的大小,记为第二制动力矩,其中,TSMC为车轮的滑模制动总力矩;
制动控制模块,用于在所述汽车的制动防抱死过程中,根据所述第一制动力矩和所述第二制动力矩,控制液压制动***输出的制动力矩和各车轮轮毂电机输出的制动力矩。
可选的,所述***还包括:
滑模制动总力矩计算模块,用于根据
Figure BDA0002323830830000041
计算各车轮的滑模制动总力矩,其中,
Figure BDA0002323830830000042
ksat(s)为引入的边界层,φ为边界层的厚度,s=λ-λd,λ为当前滑移率,λd为最优滑移率,Fz为车轮垂向力,μ为当前利用附着系数,r为车轮滚动半径,J为车轮转动惯量,
Figure BDA0002323830830000045
为车轮加速度,k为系数。
可选的,所述车速滤波模块,具体包括:
车速滤波单元,用于以
Figure BDA0002323830830000043
为状态方程,以
Figure BDA0002323830830000044
为观测方程,采用无迹卡尔曼滤波对汽车的当前车速进行滤波,状态方程中的ax、ay、γ、Mz根据动力学模型计算得到,其中,ax为汽车纵向加速度,ay为汽车侧向加速度,γ为汽车横摆角度,ax测量为传感器测量得到的汽车纵向加速度,ay测量为传感器测量得到的汽车侧向加速度,γ测量为传感器测量得到的汽车横摆角度,vx为汽车纵向行驶速度,vy为汽车侧向行驶速度,Mz为绕z轴的旋转力矩,T为传感器信号的采样频率,k为迭代次数。
可选的,所述***还包括:
标准曲线获取模块,用于获取标准路面的滑移率与利用附着系数关系曲线,记为标准曲线,所述标准路面包括干沥青路面、湿沥青路面、湿鹅卵石路面和冰路面;
样本数据获取模块,用于获取样本数据,所述样本数据包括滑移率和与其相对应的利用附着系数;
相似度确定模块,用于确定所述样本数据与各标准曲线的相似度;
权重系数确定模块,用于根据所述相似度确定所述样本数据对各标准曲线的权重系数;
样本曲线确定模块,用于根据各所述标准曲线以及各所述权重系数,确定所述样本数据对应的滑移率与利用附着系数关系曲线,记为样本曲线;
标签确定模块,用于确定所述样本曲线上利用附着系数最大的点以及该点所对应的滑移率,分别记为最大利用附着系数和最优滑移率;并将所述最大利用附着系数和所述最优滑移率作为所述样本数据的标签;
神经网络训练模块,用于根据所述样本数据以及其标签对神经网络进行训练,得到自适应神经网络模型。
根据本发明提供的具体实施例,本发明公开了以下技术效果:本发明提供的分布式驱动电动汽车制动防抱死控制方法及***,在汽车的制动防抱死过程中,根据汽车各车轮的滑移率和利用附着系数,采用自适应神经网络模型确定当前路面的最优滑移率和最大利用附着系数,各车轮采用根据其求取的最优滑移率和最大利用附着系数确定其需要液压制动***提供的制动力矩以及需要其自身轮毂电机提供的制动力矩,最后,根据各车轮的需求,控制液压制动***以及各车轮的轮毂电机分别提供相应的制动力矩,实现对汽车制动防抱死过程的复合控制。本发明采用液压制动***和轮毂电机的复合制动控制,既具有液压制动防抱死的可靠性,又具有电机响应迅速的优势,提高了电动汽车的紧急制动效果。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为现有技术中液压制动管路油压建立过程图;
图2为本发明实施例复合制动防抱死控制结构简图;
图3为本发明实施例复合ABS控制算法简图;
图4为本发明实施例复合ABS控制算法流程图;
图5为本发明实施例分布式驱动电动汽车制动防抱死控制方法流程图;
图6为7自由度车辆动力学模型结构示意图;
图7为本发明实施例标准路面利用附着系数与滑移率关系曲线图;
图8为本发明实施例车轮制动模型图;
图9为本发明实施例分布式驱动电动汽车制动防抱死控制***结构图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
本发明的第一方面提供了一种分布式驱动电动汽车制动防抱死控制方法,如图1所示,该方法基于一种具有复合制动***的电动汽车,该复合制动***既包括液压制动***还包括电机制动,如图2所示,该复合制动***具体包括:制动踏板、制动油压传感器、制动管路、ABS压力调节模块、液压控制器、制动卡钳、轮速传感器、加速度传感器、整车控制器、复合制动控制器、电机控制器、驱动电机。
本发明的控制算法简图如图3所示,包括车速估计模块、路面识别模块、复合ABS控制策略模块。
本发明的电液复合ABS控制算法的流程如图4所示。整车控制器采集整车的轮速信号、加速度信号、制动油压信号以及电机的实时转矩并发送给复合制动器共同享用,然后通过自适应无迹卡尔曼滤波(AUKF)估算实时车速并进行整车动力学计算。将估算得到的各个车轮的滑移率λij和计算得到的当前各个车轮的利用附着系数μij传递给路面识别模块,通过自适应神经网络算法识别出当前路面的最大利用附着系数μmax和最优滑移率λopt。以最优滑移率为目标,通过算法独立进行电机力矩与液压制动力矩的控制,将控制结果返回给整车,不断迭代计算,直至完成整个电液复合制动防抱死过程。
具体包括以下步骤,如图5所示:
步骤501:采用无迹卡尔曼滤波对汽车的当前车速进行滤波;
步骤502:根据所述汽车的当前车速确定各车轮的当前滑移率;
步骤503:计算所述汽车各车轮的当前利用附着系数;
步骤504:将所述当前滑移率和所述当前利用附着系数输入自适应神经网络模型,得到当前路面的最优滑移率和最大利用附着系数;
步骤505:根据Thb=(μmaxFz-Freg)r计算各车轮需要液压制动***提供的制动力矩Thb的大小,记为第一制动力矩,其中,Freg=αFregmax,Fregmax为当前电机可产生的最大回馈制动力,α为回馈系数,Freg为当前电机可回馈的制动力,μmax为所述最大利用附着系数;
步骤506:根据Teb=TSMC-Thb计算各车轮需要其轮毂电机提供的制动力矩Teb的大小,记为第二制动力矩,其中,TSMC为车轮的滑模制动总力矩;
步骤507:在所述汽车的制动防抱死过程中,根据所述第一制动力矩和所述第二制动力矩,控制液压制动***输出的制动力矩和各车轮轮毂电机输出的制动力矩。
在上述实施例的基础上,本实施例的步骤501具体包括:
(1)根据动力学模型计算相关参数:ax、ay、γ、Mz
根据图6所示的7自由度动力学模型结构,对车体进行整车的纵向、侧向和横摆动力学分析,根据达朗贝尔原理,各动力学方程如下:
纵向运动(车辆沿x轴线方向的运动):
Figure BDA0002323830830000081
max=(Fx_fl+Fx_fr)cosδ-(Fy_fl+Fy_fr)sinδ+Fx_rl+Fx_rr (2)
侧向运动(车辆沿y轴线方向的运动):
Figure BDA0002323830830000082
may=(Fx_fl+Fx_fr)sinδ+(Fy_fl+Fy_fr)cosδ+Fy_rl+Fy_rr (4)
横摆运动(车辆绕通过质心的z轴线的转动):
Figure BDA0002323830830000083
Figure BDA0002323830830000084
Figure BDA0002323830830000085
式中,ij=fl,fr,rl,rr,分别表示车的前左、前右、后左、后右轮;Fxij为纵向力;Fyij为侧向力;Mz为旋转力矩;vx为汽车纵向行驶速度;vy为汽车侧向行驶速度;Jij为各车轮的转动惯量;Bf、Br分别为前、后轴轴距;a、b分别为质心到前、后轴的距离;ax为车体纵向加速度;ay为侧向加速度;Tij为各车轮制动转矩;ωij为各车轮旋转角速度;δ为前轮转角;γ为车体横摆角速度。
Figure BDA0002323830830000086
Figure BDA0002323830830000091
式中,αij为各个车轮的轮胎侧偏角;Cx、Cy分别为轮胎的纵向刚度和侧偏刚度;λij为各车轮的滑移率;μij_max为当前路面的峰值利用附着系数。
(2)车速估计
Figure BDA0002323830830000092
为状态方程,以
Figure BDA0002323830830000093
为观测方程,采用无迹卡尔曼滤波对汽车的当前车速进行滤波,状态方程中的ax、ay、γ、Mz根据动力学模型计算得到,其中,ax为汽车纵向加速度,ay为汽车侧向加速度,γ为汽车横摆角度,ax测量为传感器测量得到的汽车纵向加速度,ay测量为传感器测量得到的汽车侧向加速度,γ测量为传感器测量得到的汽车横摆角度,vx为汽车纵向行驶速度,vy为汽车侧向行驶速度,Mz为绕z轴的旋转力矩,T为传感器信号的采样频率,k为迭代次数。
1)初始化
设待估状态变量为Xs=[vx vy γ ax ay Mz]T,状态变量的初始估计值为
Figure BDA0002323830830000097
控制输入u=[Tij ωij δ]T;测量输出y=[ax ay γ];初始状态估计误差协方差为Ps,0,即
Figure BDA0002323830830000094
Figure BDA0002323830830000095
2)状态变量对称采样
假设状态变量个数为L,采用对称采样方法行程2L+1个Sigma点集。
3)状态预测
Figure BDA0002323830830000096
Figure BDA0002323830830000101
Figure BDA0002323830830000102
式中,f为离散化后的***状态方程;qs,k-1为状态估计过程噪声均值;Qs,k-1为状态估计过程噪声协方差。与传统无迹卡尔曼滤波算法所不同的是,此处qs,k-1和Qs,k-1是两个变量,利用Sage-Husa噪声估计器对过程噪声的统计信息进行估计和更新,从而使算法无需提前获知噪声的信息,普适性大大提高。
4)状态更新
Figure BDA0002323830830000103
Figure BDA0002323830830000104
Figure BDA0002323830830000105
Figure BDA0002323830830000106
Figure BDA0002323830830000107
Figure BDA0002323830830000108
Figure BDA0002323830830000109
式中:Rs为测量误差协方差矩阵,ys,k为输出信号,此处为ys,k=[ax,k ay,k γk]T。车速估计值v便可获取。
在上述实施例的基础上,本实施例中的步骤502具体包括:
车速估计值v确定之后,便可获得各个轮胎的实时滑移率λij
Figure BDA00023238308300001010
其中,ωij为车轮轮速,r为车轮半径,v为估算车速,ij=fl,fr,rl,rr,分别表示左前轮、右前轮、左后轮、右后轮。
在上述实施例的基础上,本实施例中的步骤503具体包括:
利用附着系数μij满足:
Figure BDA0002323830830000111
其中,Fx_ij、Fz_ij分别为各个车轮的纵向力和垂向力。
在上述实施例的基础上,本实施例提供的方法还包括:自适应神经网络模型的训练,具体包括:
获取标准路面的滑移率与利用附着系数关系曲线,记为标准曲线,所述标准路面包括干沥青路面、湿沥青路面、湿鹅卵石路面和冰路面;
获取样本数据,所述样本数据包括滑移率和与其相对应的利用附着系数;
确定所述样本数据与各标准曲线的相似度;
根据所述相似度确定所述样本数据对各标准曲线的权重系数;
根据各所述标准曲线以及各所述权重系数,确定所述样本数据对应的滑移率与利用附着系数关系曲线,记为样本曲线;
确定所述样本曲线上利用附着系数最大的点以及该点所对应的滑移率,分别记为最大利用附着系数和最优滑移率;并将所述最大利用附着系数和所述最优滑移率作为所述样本数据的标签;
根据所述样本数据以及其标签对神经网络进行训练,得到自适应神经网络模型。
将Burckhardt提出的曲线表达式作为路面识别算法的基础并以此获得4种标准路面的最大利用附着系数和最优滑移率。
Figure BDA0002323830830000112
式中,C1,C2,C3为拟合系数,不同的路面采用不同的拟合系数。
本发明选择干沥青、湿沥青、湿鹅卵石、冰4种路面作为标准路面,具体的参数见下表所示:
表一4种标准路面的拟合参数
Figure BDA0002323830830000113
Figure BDA0002323830830000121
上述4种标准路面的μ-λ曲线如图7所示。
根据计算得到的实时滑移率λij和利用附着系数μij,通过通过对比4条标准路面,获取当前路面的最大利用附着系数μmax和最优滑移率λopt
Figure BDA0002323830830000122
Figure BDA0002323830830000123
其中,μmax为当前路面的最大利用附着系数;λopt为当前路面的最优滑移率;μ1,μ2,μ3,μ4和λ1,λ2,λ3,λ4分别为干沥青路面、湿沥青路面、湿鹅卵石路面和冰路面的最大利用附着系数和最优滑移率;x1,x2,x3,x4为当前路面与4种标准路面的相似程度(取值范围为[0,1])。
在上述实施例的基础上,本实施例提供的方法还包括:
根据
Figure BDA0002323830830000124
计算各车轮的滑模制动总力矩,其中,
Figure BDA0002323830830000125
ksat(s)为引入的边界层,φ为边界层的厚度,s=λ-λd,λ为当前滑移率,λd为最优滑移率,Fz为车轮垂向力,μ为当前利用附着系数,r为车轮滚动半径,J为车轮转动惯量,
Figure BDA0002323830830000126
为车轮加速度,k为系数。
图8为制动时的车轮模型,其动力学方程为:
Figure BDA0002323830830000127
Figure BDA0002323830830000131
TSMC=Thb+Teb (29)
式中,ω为车轮角速度,Thb为液压制动力矩,Teb为电机制动力矩,Fx为车轮纵向力,r为车轮滚动半径,J为车轮转动惯量,Fz为车轮垂向力,μ(λ)为滑移率为λ时的路面利用附着系数,v为车速,TSMC为制动总力矩。
选取滑模控制的切换函数如下:
s=λ-λopt (30)
式中,λopt为目标滑移率。
对式(28)滑移率λ求导,得
Figure BDA0002323830830000132
当***在滑模面上滑动时,满足
Figure BDA0002323830830000133
则有
Figure BDA0002323830830000134
将式(27)、(28)、(29)、(30)带入上式,得
Figure BDA0002323830830000135
则有滑模制动总力矩
Figure BDA0002323830830000136
式中,ksat(s)为引入的边界层,用于消除抖振。其中
Figure BDA0002323830830000137
式中,φ为边界层的厚度。
为保证***能够达到稳定,按等速趋近律η设计的切换函数应满足
Figure BDA0002323830830000141
式中,η>0。将式(18)、(19)带入式(22)可得
Figure BDA0002323830830000142
在上述实施例的基础上,本实施例的步骤505具体包括:
确定每个车轮的液压制动力矩Thb
在电液复合制动***中,由液压制动提供基础制动力。液压制动力矩Thb如下式:
Thb=(μFz-Freg)r (38)
式中,Freg为当前电机可回馈的制动力。Freg满足
Freg=αFregmax (39)
式中,Fregmax为当前电机可产生的最大回馈制动力,α为回馈系数,由动力电池SOC、电机转速、电机温度等共同决定。
在控制过程中,为了减少液压制动的频繁控制,设置了液压控制死区,即在目标液压力的附近设置不进行压力控制的门限值。若测得的轮缸压力位于目标轮缸压力两侧的控制死区内,则不进行液压调节,目标偏差由回馈制动补充,这样可以减少液压制动的控制频率,改善防抱过程中的舒适性。
在上述实施例的基础上,本实施例的步骤506具体包括:
确定每个车轮的电机制动力矩Teb
根据式(17),可得电机制动力矩
Teb=TSMC-Thb (40)
目前电动汽车电机回馈制动力矩虽然不能满足总制动转矩需求,但是可以满足滑模控制策略中转矩波动的需求。经过步骤501至步骤507的不断迭代,可完成电液复合ABS,从而缩短制动距离,提高安全性。由于该算法可减少液压制动的控制频率,因此可以减少对整车的冲击,提高舒适性。
本发明的第二方面提供了一种分布式驱动电动汽车制动防抱死控制***,如图9所示,该***包括:
车速滤波模块901,用于采用无迹卡尔曼滤波对汽车的当前车速进行滤波;
滑移率确定模块902,用于根据所述汽车的当前车速确定各车轮的当前滑移率;
利用附着系数计算模块903,用于计算所述汽车各车轮的当前利用附着系数;
最优参数确定模块904,用于将所述当前滑移率和所述当前利用附着系数输入自适应神经网络模型,得到当前路面的最优滑移率和最大利用附着系数;
第一制动力矩计算模块905,用于根据Thb=(μmaxFz-Freg)r计算各车轮需要液压制动***提供的制动力矩Thb的大小,记为第一制动力矩,其中,Freg=αFregmax,Fregmax为当前电机可产生的最大回馈制动力,α为回馈系数,Freg为当前电机可回馈的制动力,μmax为所述最大利用附着系数;
第二制动力矩计算模块906,用于根据Teb=TSMC-Thb计算各车轮需要其轮毂电机提供的制动力矩Teb的大小,记为第二制动力矩,其中,TSMC为车轮的滑模制动总力矩;
制动控制模块907,用于在所述汽车的制动防抱死过程中,根据所述第一制动力矩和所述第二制动力矩,控制液压制动***输出的制动力矩和各车轮轮毂电机输出的制动力矩。
在上述实施例的基础上,本实施例提供的***还包括:
滑模制动总力矩计算模块,用于根据
Figure BDA0002323830830000151
计算各车轮的滑模制动总力矩,其中,
Figure BDA0002323830830000152
ksat(s)为引入的边界层,φ为边界层的厚度,s=λ-λd,λ为当前滑移率,λd为最优滑移率,Fz为车轮垂向力,μ为当前利用附着系数,r为车轮滚动半径,J为车轮转动惯量,
Figure BDA0002323830830000153
为车轮加速度,k为系数。
在上述实施例的基础上,本实施例的车速滤波模块901,具体包括:
车速滤波单元,用于以
Figure BDA0002323830830000161
为状态方程,以
Figure BDA0002323830830000162
为观测方程,采用无迹卡尔曼滤波对汽车的当前车速进行滤波,状态方程中的ax、ay、γ、Mz根据动力学模型计算得到,其中,ax为汽车纵向加速度,ay为汽车侧向加速度,γ为汽车横摆角度,ax测量为传感器测量得到的汽车纵向加速度,ay测量为传感器测量得到的汽车侧向加速度,γ测量为传感器测量得到的汽车横摆角度,vx为汽车纵向行驶速度,vy为汽车侧向行驶速度,Mz为绕z轴的旋转力矩,T为传感器信号的采样频率,k为迭代次数。
在上述实施例的基础上,本实施例提供的***还包括:
标准曲线获取模块,用于获取标准路面的滑移率与利用附着系数关系曲线,记为标准曲线,所述标准路面包括干沥青路面、湿沥青路面、湿鹅卵石路面和冰路面;
样本数据获取模块,用于获取样本数据,所述样本数据包括滑移率和与其相对应的利用附着系数;
相似度确定模块,用于确定所述样本数据与各标准曲线的相似度;
权重系数确定模块,用于根据所述相似度确定所述样本数据对各标准曲线的权重系数;
样本曲线确定模块,用于根据各所述标准曲线以及各所述权重系数,确定所述样本数据对应的滑移率与利用附着系数关系曲线,记为样本曲线;
标签确定模块,用于确定所述样本曲线上利用附着系数最大的点以及该点所对应的滑移率,分别记为最大利用附着系数和最优滑移率;并将所述最大利用附着系数和所述最优滑移率作为所述样本数据的标签;
神经网络训练模块,用于根据所述样本数据以及其标签对神经网络进行训练,得到自适应神经网络模型。
本发明既保留了传统液压制动防抱死的可靠性,又体现出电机力矩响应迅速、大小精确可调的优势,通过对二者联合控制,可以极大地提升制动防抱死的性能,提高紧急制动的稳定性和缩短刹车距离。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的***而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。

Claims (8)

1.一种分布式驱动电动汽车制动防抱死控制方法,其特征在于,包括:
采用无迹卡尔曼滤波对汽车的当前车速进行滤波;
根据所述汽车的当前车速确定各车轮的当前滑移率;
计算所述汽车各车轮的当前利用附着系数;
将所述当前滑移率和所述当前利用附着系数输入自适应神经网络模型,得到当前路面的最优滑移率和最大利用附着系数;
根据Thb=(μmaxFz-Freg)r计算各车轮需要液压制动***提供的制动力矩Thb的大小,记为第一制动力矩,其中,Freg=αFregmax,Fregmax为当前电机可产生的最大回馈制动力,α为回馈系数,Freg为当前电机可回馈的制动力,μmax为所述最大利用附着系数;
根据Teb=TSMC-Thb计算各车轮需要其轮毂电机提供的制动力矩Teb的大小,记为第二制动力矩,其中,TSMC为车轮的滑模制动总力矩;
在所述汽车的制动防抱死过程中,根据所述第一制动力矩和所述第二制动力矩,控制液压制动***输出的制动力矩和各车轮轮毂电机输出的制动力矩。
2.根据权利要求1所述的分布式驱动电动汽车制动防抱死控制方法,其特征在于,所述方法还包括:根据
Figure FDA0002323830820000011
计算各车轮的滑模制动总力矩,其中,
Figure FDA0002323830820000012
ksat(s)为引入的边界层,φ为边界层的厚度,s=λ-λd,λ为当前滑移率,λd为最优滑移率,Fz为车轮垂向力,μ为当前利用附着系数,r为车轮滚动半径,J为车轮转动惯量,
Figure FDA0002323830820000013
为车轮加速度,k为系数。
3.根据权利要求1所述的分布式驱动电动汽车制动防抱死控制方法,其特征在于,所述采用无迹卡尔曼滤波对汽车的当前车速进行滤波,具体包括:
Figure FDA0002323830820000021
为状态方程,以
Figure FDA0002323830820000022
为观测方程,采用无迹卡尔曼滤波对汽车的当前车速进行滤波,状态方程中的ax、ay、γ、Mz根据动力学模型计算得到,其中,ax为汽车纵向加速度,ay为汽车侧向加速度,γ为汽车横摆角度,ax测量为传感器测量得到的汽车纵向加速度,ay测量为传感器测量得到的汽车侧向加速度,γ测量为传感器测量得到的汽车横摆角度,vx为汽车纵向行驶速度,vy为汽车侧向行驶速度,Mz为绕z轴的旋转力矩,T为传感器信号的采样频率,k为迭代次数。
4.根据权利要求1所述的分布式驱动电动汽车制动防抱死控制方法,其特征在于,所述方法还包括:
获取标准路面的滑移率与利用附着系数关系曲线,记为标准曲线,所述标准路面包括干沥青路面、湿沥青路面、湿鹅卵石路面和冰路面;
获取样本数据,所述样本数据包括滑移率和与其相对应的利用附着系数;
确定所述样本数据与各标准曲线的相似度;
根据所述相似度确定所述样本数据对各标准曲线的权重系数;
根据各所述标准曲线以及各所述权重系数,确定所述样本数据对应的滑移率与利用附着系数关系曲线,记为样本曲线;
确定所述样本曲线上利用附着系数最大的点以及该点所对应的滑移率,分别记为最大利用附着系数和最优滑移率;并将所述最大利用附着系数和所述最优滑移率作为所述样本数据的标签;
根据所述样本数据以及其标签对神经网络进行训练,得到自适应神经网络模型。
5.一种分布式驱动电动汽车制动防抱死控制***,其特征在于,包括:
车速滤波模块,用于采用无迹卡尔曼滤波对汽车的当前车速进行滤波;
滑移率确定模块,用于根据所述汽车的当前车速确定各车轮的当前滑移率;
利用附着系数计算模块,用于计算所述汽车各车轮的当前利用附着系数;
最优参数确定模块,用于将所述当前滑移率和所述当前利用附着系数输入自适应神经网络模型,得到当前路面的最优滑移率和最大利用附着系数;
第一制动力矩计算模块,用于根据Thb=(μmaxFz-Freg)r计算各车轮需要液压制动***提供的制动力矩Thb的大小,记为第一制动力矩,其中,Freg=αFregmax,Fregmax为当前电机可产生的最大回馈制动力,α为回馈系数,Freg为当前电机可回馈的制动力,μmax为所述最大利用附着系数;
第二制动力矩计算模块,用于根据Teb=TSMC-Thb计算各车轮需要其轮毂电机提供的制动力矩Teb的大小,记为第二制动力矩,其中,TSMC为车轮的滑模制动总力矩;
制动控制模块,用于在所述汽车的制动防抱死过程中,根据所述第一制动力矩和所述第二制动力矩,控制液压制动***输出的制动力矩和各车轮轮毂电机输出的制动力矩。
6.根据权利要求5所述的分布式驱动电动汽车制动防抱死控制***,其特征在于,所述***还包括:
滑模制动总力矩计算模块,用于根据
Figure FDA0002323830820000031
计算各车轮的滑模制动总力矩,其中,
Figure FDA0002323830820000032
ksat(s)为引入的边界层,φ为边界层的厚度,s=λ-λd,λ为当前滑移率,λd为最优滑移率,Fz为车轮垂向力,μ为当前利用附着系数,r为车轮滚动半径,J为车轮转动惯量,
Figure FDA0002323830820000033
为车轮加速度,k为系数。
7.根据权利要求5所述的分布式驱动电动汽车制动防抱死控制***,其特征在于,所述车速滤波模块,具体包括:
车速滤波单元,用于以
Figure FDA0002323830820000041
为状态方程,以
Figure FDA0002323830820000042
为观测方程,采用无迹卡尔曼滤波对汽车的当前车速进行滤波,状态方程中的ax、ay、γ、Mz根据动力学模型计算得到,其中,ax为汽车纵向加速度,ay为汽车侧向加速度,γ为汽车横摆角度,ax测量为传感器测量得到的汽车纵向加速度,ay测量为传感器测量得到的汽车侧向加速度,γ测量为传感器测量得到的汽车横摆角度,vx为汽车纵向行驶速度,vy为汽车侧向行驶速度,Mz为绕z轴的旋转力矩,T为传感器信号的采样频率,k为迭代次数。
8.根据权利要求5所述的分布式驱动电动汽车制动防抱死控制***,其特征在于,所述***还包括:
标准曲线获取模块,用于获取标准路面的滑移率与利用附着系数关系曲线,记为标准曲线,所述标准路面包括干沥青路面、湿沥青路面、湿鹅卵石路面和冰路面;
样本数据获取模块,用于获取样本数据,所述样本数据包括滑移率和与其相对应的利用附着系数;
相似度确定模块,用于确定所述样本数据与各标准曲线的相似度;
权重系数确定模块,用于根据所述相似度确定所述样本数据对各标准曲线的权重系数;
样本曲线确定模块,用于根据各所述标准曲线以及各所述权重系数,确定所述样本数据对应的滑移率与利用附着系数关系曲线,记为样本曲线;
标签确定模块,用于确定所述样本曲线上利用附着系数最大的点以及该点所对应的滑移率,分别记为最大利用附着系数和最优滑移率;并将所述最大利用附着系数和所述最优滑移率作为所述样本数据的标签;
神经网络训练模块,用于根据所述样本数据以及其标签对神经网络进行训练,得到自适应神经网络模型。
CN201911308436.6A 2019-12-18 2019-12-18 一种分布式驱动电动汽车制动防抱死控制方法及*** Active CN110949344B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201911308436.6A CN110949344B (zh) 2019-12-18 2019-12-18 一种分布式驱动电动汽车制动防抱死控制方法及***

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201911308436.6A CN110949344B (zh) 2019-12-18 2019-12-18 一种分布式驱动电动汽车制动防抱死控制方法及***

Publications (2)

Publication Number Publication Date
CN110949344A true CN110949344A (zh) 2020-04-03
CN110949344B CN110949344B (zh) 2021-01-29

Family

ID=69982554

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201911308436.6A Active CN110949344B (zh) 2019-12-18 2019-12-18 一种分布式驱动电动汽车制动防抱死控制方法及***

Country Status (1)

Country Link
CN (1) CN110949344B (zh)

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111824095A (zh) * 2020-06-14 2020-10-27 长春理工大学 四轮轮毂电动汽车电液复合制动防抱死协调优化控制方法
CN111976677A (zh) * 2020-09-02 2020-11-24 哈尔滨理工大学 一种纯电动汽车复合制动防抱死控制***及控制方法
CN112298135A (zh) * 2020-10-30 2021-02-02 重庆长安汽车股份有限公司 一种车辆制动距离优化方法
CN112406824A (zh) * 2020-11-27 2021-02-26 中车唐山机车车辆有限公司 列车制动方法、装置、设备及存储介质
CN112810588A (zh) * 2021-01-26 2021-05-18 北京理工大学 一种分布式驱动电动汽车电液复合制动防抱死方法及***
CN114228721A (zh) * 2022-01-29 2022-03-25 苏州挚途科技有限公司 路面附着系数的计算方法、装置和***
CN114537346A (zh) * 2021-12-11 2022-05-27 清华大学 电动轮滑移率确定方法、装置、计算机设备和存储介质
CN114701463A (zh) * 2022-04-11 2022-07-05 浙江理工大学 区间二型模糊神经网络复合制动防抱死控制***及方法
CN114715093A (zh) * 2022-04-24 2022-07-08 郑州轻工业大学 一种基于神经网络自适应估计的汽车防抱死制动方法
CN114771476A (zh) * 2022-05-27 2022-07-22 烟台大学 一种基于非线性自适应滑模的汽车防抱死制动方法
CN116118687A (zh) * 2023-03-31 2023-05-16 吉林大学 8×8电驱动车辆紧急制动工况下滑移率控制***及方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH07156782A (ja) * 1993-12-07 1995-06-20 Honda Motor Co Ltd 路面状態検出装置およびこの装置を利用したアンチロックブレーキシステム
CN109421552A (zh) * 2017-09-04 2019-03-05 郑州宇通客车股份有限公司 分布式驱动电动汽车路面自适应驱动防滑控制方法及***
CN109664774A (zh) * 2018-03-28 2019-04-23 北京理工大学 一种用于分布式驱动电动汽车的驱动防滑控制方法及***

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH07156782A (ja) * 1993-12-07 1995-06-20 Honda Motor Co Ltd 路面状態検出装置およびこの装置を利用したアンチロックブレーキシステム
CN109421552A (zh) * 2017-09-04 2019-03-05 郑州宇通客车股份有限公司 分布式驱动电动汽车路面自适应驱动防滑控制方法及***
CN109664774A (zh) * 2018-03-28 2019-04-23 北京理工大学 一种用于分布式驱动电动汽车的驱动防滑控制方法及***

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
张厚忠,等: "电动轮汽车电液复合ABS控制的联合仿真", 《广西大学学报(自然科学版)》 *
张磊: "电动轿车再生制动与防抱死协调控制算法研究", 《工程科技II辑》 *
彭晓燕,等: "基于RBF神经网络的最佳滑移率在线计算方法", 《机械工程学报》 *
潘劲,等: "分布式驱动电动汽车电液复合ABS控制研究", 《汽车技术》 *

Cited By (17)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111824095B (zh) * 2020-06-14 2022-07-05 长春理工大学 四轮轮毂电动汽车电液复合制动防抱死协调优化控制方法
CN111824095A (zh) * 2020-06-14 2020-10-27 长春理工大学 四轮轮毂电动汽车电液复合制动防抱死协调优化控制方法
CN111976677A (zh) * 2020-09-02 2020-11-24 哈尔滨理工大学 一种纯电动汽车复合制动防抱死控制***及控制方法
CN112298135A (zh) * 2020-10-30 2021-02-02 重庆长安汽车股份有限公司 一种车辆制动距离优化方法
CN112298135B (zh) * 2020-10-30 2022-07-08 重庆长安汽车股份有限公司 一种车辆制动距离优化方法
CN112406824A (zh) * 2020-11-27 2021-02-26 中车唐山机车车辆有限公司 列车制动方法、装置、设备及存储介质
CN112810588A (zh) * 2021-01-26 2021-05-18 北京理工大学 一种分布式驱动电动汽车电液复合制动防抱死方法及***
CN114537346A (zh) * 2021-12-11 2022-05-27 清华大学 电动轮滑移率确定方法、装置、计算机设备和存储介质
CN114537346B (zh) * 2021-12-11 2022-12-09 清华大学 电动轮滑移率确定方法、装置、计算机设备和存储介质
CN114228721A (zh) * 2022-01-29 2022-03-25 苏州挚途科技有限公司 路面附着系数的计算方法、装置和***
CN114701463A (zh) * 2022-04-11 2022-07-05 浙江理工大学 区间二型模糊神经网络复合制动防抱死控制***及方法
CN114715093A (zh) * 2022-04-24 2022-07-08 郑州轻工业大学 一种基于神经网络自适应估计的汽车防抱死制动方法
CN114715093B (zh) * 2022-04-24 2023-02-21 郑州轻工业大学 一种基于神经网络自适应估计的汽车防抱死制动方法
CN114771476A (zh) * 2022-05-27 2022-07-22 烟台大学 一种基于非线性自适应滑模的汽车防抱死制动方法
CN114771476B (zh) * 2022-05-27 2023-10-03 烟台大学 一种基于非线性自适应滑模的汽车防抱死制动方法
CN116118687A (zh) * 2023-03-31 2023-05-16 吉林大学 8×8电驱动车辆紧急制动工况下滑移率控制***及方法
CN116118687B (zh) * 2023-03-31 2024-05-07 吉林大学 8×8电驱动车辆紧急制动工况下滑移率控制***及方法

Also Published As

Publication number Publication date
CN110949344B (zh) 2021-01-29

Similar Documents

Publication Publication Date Title
CN110949344B (zh) 一种分布式驱动电动汽车制动防抱死控制方法及***
CN110382326B (zh) 高速正常行驶状况下的轮胎路面摩擦系数估算方法和轮胎路面摩擦系数估算装置
CN109606133B (zh) 基于双层控制的分布式驱动电动汽车转矩矢量控制方法
CN107719372B (zh) 基于动态控制分配的四驱电动汽车动力学多目标控制***
CN108128308B (zh) 一种分布式驱动电动汽车的车辆状态估计***及方法
CN111391822B (zh) 一种极限工况下汽车横纵向稳定性协同控制方法
Xia et al. Estimation of maximum road friction coefficient based on Lyapunov method
CN103612634B (zh) 分散式轮毂电机驱动电动汽车路面附着系数的估算方法
CN104773170A (zh) 一种车辆稳定性集成控制方法
CN112537307B (zh) 一种四轮轮毂电机自寻优驱动防滑控制方法和***
US11648933B2 (en) Method for controlling wheel slip of vehicle
CN111469670A (zh) 一种基于路面识别的电动汽车再生制动控制策略
CN105606530A (zh) 一种路面峰值附着系数测试装置及方法
Savitski et al. Influence of the tire inflation pressure variation on braking efficiency and driving comfort of full electric vehicle with continuous anti-lock braking system
CN112810588B (zh) 一种分布式驱动电动汽车电液复合制动防抱死方法及***
CN109733400A (zh) 用于在车辆中分配驱动扭矩的方法、装置和设备
KR20110125280A (ko) 경사로의 적응 순항 제어 방법
Pusca et al. Fuzzy-logic-based control applied to a hybrid electric vehicle with four separate wheel drives
Le Solliec et al. Regenerative braking optimization and wheel slip control for a vehicle with in-wheel motors
CN112660107B (zh) 一种改善分布式电驱动底盘功率循环的转矩分配方法
CN113682152B (zh) 一种分布式驱动汽车牵引力控制方法
CN111762172A (zh) 一种基于电-液参数的路面附着估计方法
CN115891961A (zh) 用于四驱汽车电动助力转向与主动动力分配集成控制方法
Chen et al. Traction control for electric vehicles: A novel control scheme
JP2008049996A (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