CN112373477B - 自动驾驶系统的冗余控制方法、自动驾驶系统、汽车、控制器及计算机可读存储介质 - Google Patents

自动驾驶系统的冗余控制方法、自动驾驶系统、汽车、控制器及计算机可读存储介质 Download PDF

Info

Publication number
CN112373477B
CN112373477B CN202011317326.9A CN202011317326A CN112373477B CN 112373477 B CN112373477 B CN 112373477B CN 202011317326 A CN202011317326 A CN 202011317326A CN 112373477 B CN112373477 B CN 112373477B
Authority
CN
China
Prior art keywords
vehicle
lane
function
formula
safe
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
CN202011317326.9A
Other languages
English (en)
Other versions
CN112373477A (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.)
Chongqing Changan Automobile Co Ltd
Original Assignee
Chongqing Changan Automobile Co Ltd
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 Chongqing Changan Automobile Co Ltd filed Critical Chongqing Changan Automobile Co Ltd
Priority to CN202011317326.9A priority Critical patent/CN112373477B/zh
Priority to RU2021104308A priority patent/RU2759453C1/ru
Priority to PCT/CN2020/133607 priority patent/WO2022104925A1/zh
Publication of CN112373477A publication Critical patent/CN112373477A/zh
Application granted granted Critical
Publication of CN112373477B publication Critical patent/CN112373477B/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
    • B60W40/00Estimation or calculation of non-directly measurable driving parameters for road vehicle drive control systems not related to the control of a particular sub unit, e.g. by using mathematical models
    • 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
    • B60W60/00Drive control systems specially adapted for autonomous road vehicles
    • 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/0002Automatic control, details of type of controller or control system architecture

Abstract

本方案涉及一种自动驾驶系统的冗余控制方法、自动驾驶系统、汽车、控制器及计算机可读存储介质,用于在自动驾驶系统故障且故障未超过容限时,仍然能够保持沿着车道中心线行驶。该方法包括:在自动驾驶功能激活时,根据前视摄像头和雷达采集到的信息拟合出使车辆沿车道中心线行驶的安全路径方程以及确定车辆当前所在车道的车道宽度;根据安全路径方程的系数和确定出的车道宽度,确定车辆当前具体所处的偏置区域;根据车辆所处的偏置区域和车辆在上一控制周期是否激活第一冗余安全控制模块功能,确定主控制子系统当前的故障等级;根据主控制子系统当前的故障等级,激活对应等级的冗余控制功能,并按照对应等级的冗余控制功能进行车辆控制。

Description

自动驾驶系统的冗余控制方法、自动驾驶系统、汽车、控制器 及计算机可读存储介质
技术领域
本发明用于汽车自动驾驶控制和安全控制领域,具体是一种自动驾驶系统在传感器/主控制子系统/执行器存在故障且故障有界情况下的使车辆紧急车道保持和减速的冗余控制方法、自动驾驶系统、汽车、控制器及计算机可读存储介质。
背景技术
中/高等级(2.5级以上)自动驾驶系统安全性依赖于各感知传感器的正确输出。作为车道保持的主要传感器的前视摄像头,对于车道线的识别一般多基于深度学习系统,或深度学习系统和经典图像识别技术的结合,容易因为光照,明暗交替,路面积水反射,道路坡度,车道标识等因素影响,而对车道线尤其是车道中心线曲率出现错误的识别。并且,多传感器感知融合系统复杂度高,出现硬件或软件故障和计算错误(比如部分寄存器状态的随机翻转)可能性不低。这导致了中/高等级自动驾驶汽车在车道保持过程中,容易因系统中的故障和错误引起安全问题;而在中/高等级自动驾驶车辆中,驾驶员/安全员容易因过于依赖自动驾驶系统而分散注意力,因此这种安全问题更加致命。因此需要实现一种在传感器对车道中心线曲率检测不准确,或在系统部分故障且故障未超过容限时,仍然能够保持车道,将车辆控制在安全区域内的方法。
发明内容
本发明实施例提供了一种自动驾驶系统的冗余控制方法、自动驾驶系统、汽车、控制器及计算机可读存储介质,用于在自动驾驶系统故障且故障未超过容限时,仍然能够保持沿着车道中心线行驶。
本发明的技术方案为:
本发明实施例提供了一种自动驾驶系统的冗余控制方法,所述自动驾驶系统包括主控制子系统和冗余控制子系统,所述方法应用于冗余控制子系统,所述方法包括:
在自动驾驶系统的自动驾驶功能激活时,确定车辆在当前所在车道内具体所处的偏置区域;偏置区域是指车道中心线左侧方向和车道中心线右侧方向各自与车道中心线之间的距离位于预定距离范围内时所组成的区域;
根据车辆当前所处的偏置区域和冗余控制子系统在上一控制周期输出的控制结果,确定主控制子系统当前的故障等级;
根据主控制子系统当前的故障等级,激活对应等级的冗余控制功能,并按照对应等级的冗余控制功能进行车辆控制。
优选地,确定车辆在当前所在车道内具体所处的偏置区域的具体步骤为:
根据前视摄像头和雷达输出到的多源信息,拟合使车辆保持车道中心线行驶的安全路径方程C0(x),并确定车辆当前所在车道的车道宽度Lw;
根据所规划的安全路径方程C0(x)中的系数C0和车道宽度Lw,确定车辆在当前所在车道内具体所处的偏置区域。
优选地,摄像头和雷达输出到的多源信息包括:相对于车辆坐标系原点规划的车道线曲线方程L0(x)、路沿曲线方程K0(x)、护栏曲线方程G0(x)和/或车流曲线方程T0(x);根据前视摄像头和雷达输出到的多源信息,拟合使车辆保持车道中心线行驶的安全路径方程C0(x),并确定车辆当前所在车道的车道宽度Lw的具体步骤为:
根据车辆的车速向前等间隔时长的选取N个纵向位置点,通过公式:xi=v×i×tLA(i=0,...,N-1)计算第i个纵向位置点与车辆坐标系原点之间的纵向距离xi,v为车辆的车速,tLA为设定的时间间隔时长;
通过公式:
E0Left(xi)=-min[|L0Left(xi)|,|K0Left(xi)|,|G0Left(xi)|,|T0Left(xi)|](i=0,...,N-1)确定与第i个纵向位置点对应的左侧横向位置点E0Left(xi);
通过公式:
E0Right(xi)=min[|L0Right(xi)|,|K0Right(xi)|,|G0Right(xi)|,|T0Right(xi)|](i=0,...,N-1)确定与第i个纵向位置点对应的右侧横向位置点E0Left(xi);
使用最小二乘拟合法将取得的车辆所在车道的左侧边缘的N个位置点(xi,E0Left(xi))拟合为车辆所在车道的左侧边缘曲线E0Left(x);
使用最小二乘拟合法将取得的车辆所在车道的右侧边缘的N个位置点(xi,E0Right(xi))拟合为车辆所在车道的右侧边缘曲线E0Right(x);
对左侧边缘曲线和右侧边缘曲线的各系数进行均值处理,通过公式:
Figure GDA0004001832220000021
获得使车辆沿着车道中心线行驶的安全路径方程C0(x);C0为车辆所在车道的车道中心线相对于车辆坐标系原点的横向偏置距离,C1为车辆所在车道的车道中心线相对车辆坐标系原点的夹角,C2为车道中心线曲率的1/2,C3为车道中心线曲率变化率的1/6;
通过公式:
Lw=min(Lw(x0),...,Lw(xN-1),Lwmin(xLB));
计算车辆当前所在车道的车道宽度Lw;Lwmin(xLB)表示在前tLB时间内的历史最小车道宽度,Lw(xi)=|E0Right(xi)-E0Left(xi)|(i=0,...,N-1)。
优选地,根据所规划的安全路径C0(x)中的系数C0和车道宽度Lw,确定车辆在当前所在车道内具体所处的偏置区域的具体步骤为:
通过公式:
Figure GDA0004001832220000031
计算以车道宽度Lw和横向偏置距离C0为基础的故障函数F(Δ);Lw为车辆当前所在车道的车道宽度,系数C0为车道中心线相对于车辆坐标系原点的横向偏置距离;
若故障函数F(Δ)小于第一设定比值D1;则确定车辆在当前所在车道内具体所处的偏置区域为第一偏置区域Z1;
若故障函数F(Δ)大于或等于第一设定比值D1且小于第二设定比值D2;则确定车辆在当前所在车道内具体所处的偏置区域为第二偏置区域Z2;
若故障函数F(Δ)大于或等于第二设定比值D2;则确定车辆在当前所在车道内具体所处的偏置区域为第三偏置区域Z3。
优选地,根据车辆当前所处的偏置区域和冗余控制子系统在上一控制周期输出的控制结果,确定主控制子系统当前的故障等级的具体步骤为:
若冗余控制子系统在上一控制周期输出的控制结果为激活第一等级冗余控制功能,或者,若冗余控制子系统在上一控制周期输出的控制结果为激活第二等级冗余控制功能或未激活冗余控制功能且车辆在当前所在车道内具体所处的偏置区域为第三偏置区域Z3,则确定主控制子系统当前的故障等级为高故障等级;
若冗余控制子系统在上一控制周期输出的控制结果为激活第二等级冗余控制功能或未激活冗余控制功能且车辆在当前所在车道内具体所处的偏置区域为第二偏置区域Z2,则进一步比对车辆姿态的第一判别函数G(S)与第二设定比值D2的大小;在车辆姿态的第一判别函数G(S)大于或等于第二设定比值D2时,确定主控制子系统当前的故障等级为高故障等级;在车辆姿态的第一判别函数G(S)小于第二设定比值D2时,确定主控制子系统当前的故障等级为中故障等级;
若冗余控制子系统在上一控制周期输出的控制结果为激活第二等级冗余控制功能或未激活冗余控制功能且车辆在当前所在车道内具体所处的偏置区域为第一偏置区域Z1,则确定主控制子系统当前的故障等级为低故障等级。
优选地,车辆姿态的第一判别函数G(S)通过公式:
G(S)=λS
计算获得;λ为识别信度,λ为常数或者基于车辆的车速v、车辆的偏航率Yr、方向盘转角θSteer、车道中心线曲率2C2以及前视摄像头识别到的车辆当前所在车道的路面类型所对应的路面附着系数
Figure GDA0004001832220000032
进行RBF神经网络拟合得到;当λ为常数时λ=0.7;当λ通过RBF神经网络拟合得到时,λ的取值范围为:λ∈[0.35,1];
S为滑模变量,通过公式:
Figure GDA0004001832220000041
计算获得;滑模超平面方程满足:
Figure GDA0004001832220000042
Lw为车辆当前所在车道的车道宽度;C0为车道中心线相对于车辆坐标系原点的横向偏置距离;C1为车辆的航向角;δ为对车道宽度的一半
Figure GDA0004001832220000043
所取的横向偏置距离安全系数,δ∈[0.5,0.9];θ为车辆偏离车道中心线行驶时,使车辆渐进的回到第一偏置区域Z1内行驶的安全航向角;安全航向角θ通过公式:/>
Figure GDA0004001832220000044
计算获得;θ∈[θminmax],θmax为保证车辆渐进的趋近平衡点的安全上限,θmin为保证车辆渐进的趋近平衡点的安全下限;aymax为车辆最大安全侧向加速度;车辆最大安全侧向加速度aymax通过公式:
Figure GDA0004001832220000045
计算获得;T为车辆的轮距,h为车辆的质心高度,g为重力加速度,
Figure GDA0004001832220000046
为前视摄像头识别到的车辆当前所在车道的路面类型对应的路面附着系数。
优选地,根据主控制子系统当前的故障等级,激活对应等级的冗余控制功能的具体步骤为:
若主控制子系统当前的故障等级为高故障等级,则激活第一等级冗余控制功能;若主控制子系统当前的故障等级为中故障等级,则激活第二等级冗余控制功能;若主控制子系统当前的故障等级为低故障等级,则不激活第一等级冗余控制功能或第二等级冗余控制功能。
优选地,按照对应等级的冗余控制功能进行车辆控制的具体步骤为:
若激活第一等级冗余控制功能,则按照计算出的综合减速度a控制车身电子稳定系统ESP减速至第一目标减速度,并按照计算出的目标方向盘转角θsteer和最大安全力矩Tmax控制电动助力转向系统EPS在柔性控制时间tSIL内进行方向盘转角控制;
若激活第二等级冗余控制功能,则根据车辆姿态的第二判别函数H(S)计算出的安全减速度a1和主控制子系统输出的减速度a2控制车身电子稳定系统ESP减速至第二目标减速度;
综合减速度a为根据车辆姿态的第二判别函数H(S)计算出的安全减速度a1和主控制子系统输出的减速度a2中的较大值;第一目标速度与当前行驶环境工况相关;第二目标速度与主控制子系统输出的目标速度相关,且第二目标速度小于主控制子系统输出的目标速度。
优选地,通过公式:
Figure GDA0004001832220000051
计算目标方向盘转角θsteer;2C2为车道中心线曲率;C1为车辆的航向角;Lw为车辆当前所在车道的车道宽度;Kθδ为车辆的转向系统传动比;
Figure GDA0004001832220000052
为自适应参数;l为车辆的轴距;v为车辆的车速;ε为自适应参数初始值,ε为常数;δ为对车道宽度的一半/>
Figure GDA0004001832220000053
所取的横向偏置距离安全系数,δ∈[0.5,0.9];S为滑模变量,通过公式/>
Figure GDA0004001832220000054
计算获得,,滑模超平面方程满足:/>
Figure GDA0004001832220000055
θ为车辆偏离车道中心线行驶时,使车辆渐进的回到第一偏置区域Z1内行驶的安全航向角;安全航向角θ通过公式:
Figure GDA0004001832220000056
计算获得;θ∈[θminmax],θmax为保证车辆渐进的趋近平衡点的安全上限,θmin为保证车辆渐进的趋近平衡点的安全下限;aymax为车辆最大安全侧向加速度;车辆最大安全侧向加速度aymax通过公式:
Figure GDA0004001832220000061
计算获得;T为车辆的轮距,h为车辆的质心高度,g为重力加速度,
Figure GDA0004001832220000062
为车辆传感系统识别到的车辆当前所在车道的路面类型对应的路面附着系数;
自适应参数
Figure GDA0004001832220000063
通过公式:
Figure GDA0004001832220000064
计算获得;
C0为车道中心线相对于车辆坐标系原点的横向偏置距离;
Figure GDA0004001832220000065
为死区为[0,η]的死区函数,[0,η]函数为零;Yr为车辆的偏航率;δ为对车道宽度的一半
Figure GDA0004001832220000066
所取的横向偏置距离安全系数,δ∈[0.5,0.9];μ为防止车辆的相对偏航率过大所设置的安全偏航率参数,K为自适应速率参数,γmax为自适应上界。
优选地,通过公式:
Figure GDA0004001832220000067
计算目标方向盘转角θsteer,2C2为车道中心线曲率;C1为车辆的航向角;Lw为车辆当前所在车道的车道宽度;Kθδ为车辆的转向系统传动比;
Figure GDA0004001832220000068
为自适应参数;l为车辆的轴距;v为车辆的车速;ε为自适应参数初始值,ε为常数;δ对车道宽度的一半/>
Figure GDA0004001832220000069
所取的横向偏置距离安全系数,δ∈[0.5,0.9];S为滑模变量,通过公式/>
Figure GDA00040018322200000610
计算获得,,滑模超平面方程满足:/>
Figure GDA0004001832220000071
θ为车辆偏离车道中心线行驶时,使车辆渐进的回到第一偏置区域Z1内行驶的安全航向角;安全航向角θ通过公式:
Figure GDA0004001832220000072
计算获得;θ∈[θminmax],θmax为保证车辆渐进的趋近平衡点的安全上限,θmin为保证车辆渐进的趋近平衡点的安全下限;aymax为车辆最大安全侧向加速度;车辆最大安全侧向加速度aymax通过公式:
Figure GDA0004001832220000073
计算获得;T为车辆的轮距,h为车辆的质心高度,g为重力加速度,
Figure GDA0004001832220000074
为前视摄像头识别到的车辆当前所在车道的路面类型对应的路面附着系数;
自适应参数
Figure GDA0004001832220000075
通过公式:
Figure GDA0004001832220000076
计算获得;
C0为车道中心线相对于车辆坐标系原点的横向偏置距离;
Figure GDA0004001832220000077
为死区为[0,η]的死区函数,[0,η]函数为零;Yr为车辆的偏航率;δ为对车道宽度的一半/>
Figure GDA0004001832220000078
所取的横向偏置距离安全系数,δ∈[0.5,0.9];μ为防止车辆的相对偏航率过大所设置的安全偏航率参数,K为自适应速率参数,γmax为自适应上界;
饱和函数为
Figure GDA0004001832220000081
类饱和函数为/>
Figure GDA0004001832220000082
Figure GDA0004001832220000083
是任意连续函数,满足f(0)=0,f(1)=1,f(-1)=-1;φ为准滑模边界层宽度,φ为常数。
优选地,根据车辆姿态的第二判别函数H(S)通过公式:
Figure GDA0004001832220000084
计算出安全减速度a1
g为重力加速度;
Figure GDA0004001832220000085
为前视摄像头识别到的车辆当前所在车道的路面类型估计出的对应的路面附着系数;a1max为车辆的最大安全减速度,/>
Figure GDA0004001832220000086
a1min为车辆的最小安全减速度,a1min为常数;tD1为从车辆从第一偏置区域Z1进入第二偏置区域Z2的时间点;tNow为当前时间点;Yr为车辆的偏航率;v为车辆的车速;
σ为根据车辆的转向能力和底盘稳定能力所确定的减速度系数,通过公式:
Figure GDA0004001832220000087
计算获得;LTRCrit为安全临界荷载偏移率,LTRCrit为常数;LTRavg为平均荷载偏移率,通过公式:/>
Figure GDA0004001832220000088
计算获得,tnow为当前时间点,tCrit为判定时间,T为车辆的轮距,h为车辆的质心高度;
α为故障估计量化积分参数,通过公式:
Figure GDA0004001832220000089
计算获得,ΔθCrit为安全临界航向角增量,ΔθCrit常数;
S为滑模变量,通过公式:
Figure GDA00040018322200000810
计算;C0表示车道中心线相对于车辆坐标系原点的横向偏置距离;C1为车辆的航向角;2C2为车道中心线曲率;Lw为车辆当前所在车道的车道宽度;δ为对车道宽度的一半/>
Figure GDA0004001832220000091
所取的横向偏置距离安全系数,δ∈[0.5,0.9];θ为车辆偏离车道中心线行驶时,使车辆渐进的回到第一偏置区域Z1内行驶的安全航向角;安全航向角θ通过公式:
Figure GDA0004001832220000092
计算获得;θ∈[θminmax],θmax为保证车辆渐进的趋近平衡点的安全上限,θmin为保证车辆渐进的趋近平衡点的安全下限;aymax为车辆最大安全侧向加速度;车辆最大安全侧向加速度aymax通过公式:
Figure GDA0004001832220000093
计算获得;T为车辆的轮距,h为车辆的质心高度,g为重力加速度,
Figure GDA0004001832220000094
为车辆传感系统识别到的车辆当前所在车道的路面类型对应的路面附着系数。
优选地,通过公式:
Tmax=min(TeRequire,Te-Limit)计算最大安全力矩Tmax
第一力矩TeRequire通过根据车辆的车速v从车速与第一力矩TeRequire的预设关系表查表获得,第二力矩Te-Limit由电动助力转向系统EPS输出;
通过公式:
Figure GDA0004001832220000095
计算得到柔性控制时间tSIL,且所述柔性控制时间tSIL满足:2s≥tSIL≥100ms。
优选地,按照计算出的方向盘转角θsteer和最大安全力矩Tmax控制电动助力转向系统EPS在柔性控制时间tSIL内进行方向盘转角控制的具体步骤为:
通过sigmoid函数插值,确定不同插值时刻t对应的插值方向盘转角指令值θSteerInterp;通过公式:
Figure GDA0004001832220000101
计算插值方向盘转角指令值θSteerInterp,插值时间t∈[0,tSIL];θSteerDifference为第一等级冗余控制功能激活的时间点对应的车辆的方向盘转角实际反馈值θSteerFB与目标方向盘转角θsteer之差;
将在柔性控制时间tSIL内的每一插值时刻t对应的插值方向盘转角θSteerInterp输出给电动助力转向系统EPS,使电动助力转向系统EPS在所述最大安全力矩Tmax条件下将车辆的方向盘转角由方向盘转角实际反馈值θSteerFB调节为所述插值方向盘转角指令值θSteerInterp
优选地,所述方法还包括:
针对自动驾驶系统为L4或L5高等级自动驾驶系统,在按照第一等级冗余控制功能或第二等级冗余控制功能进行车辆控制时,判断是否满足:车辆当前所处的偏置区域稳定为第一偏置区域Z1、车辆的车速是否小于车辆当前所在车道的最低限制车速且车辆当前所在车道的路况为良好路况;
若满足,输出提示安全员进行车辆接管的提示信息;
判断安全员是否接管车辆;
若接管,则关闭第一等级冗余控制功能或第二等级冗余控制功能,并停止向电动助力转向系统EPS和车身电子稳定系统ESP输出控制指令;
针对自动驾驶系统为L2.5或L3中等级自动驾驶系统,在按照第一等级冗余控制功能或第二等级冗余控制功能进行车辆控制时,输出提示驾驶员进行车辆接管的提示信息;
判断驾驶员是否接管车辆;
若接管,则关闭第一等级冗余控制功能或第二等级冗余控制功能,并停止向电动助力转向系统EPS和车身电子稳定系统ESP输出控制指令。
本发明实施例还提供了一种自动驾驶系统,所述自动驾驶系统包括主控制子系统和冗余控制子系统,所述冗余控制子系统包括:
第一确定模块,用于在自动驾驶系统的自动驾驶功能激活时,确定车辆在当前所在车道内具体所处偏置区域;偏置区域是指车道中心线左侧方向和车道中心线右侧方向各自与车道中心线之间的距离位于预定距离范围内时所组成的区域;
第二确定模块,用于根据车辆当前所处的偏置区域和冗余控制子系统在上一控制周期输出的控制结果,确定主控制子系统当前的故障等级;
仲裁模块,用于根据主控制子系统当前的故障等级,激活对应等级的冗余控制功能,并按照对应等级的冗余控制功能进行车辆控制。
优选地,第一确定模块包括:
拟合单元,用于根据前视摄像头和雷达输出到的多源信息,拟合使车辆保持车道中心线行驶的安全路径方程C0(x),并确定车辆当前所在车道的车道宽度Lw;
第一确定单元,根据所规划的安全路径方程C0(x)中的系数C0和车道宽度Lw,确定车辆在当前所在车道内具体所处的偏置区域。
优选地,摄像头和雷达输出到的多源信息包括:相对于车辆坐标系原点规划的车道线曲线方程L0(x)、路沿曲线方程K0(x)、护栏曲线方程G0(x)和/或车流曲线方程T0(x);拟合单元具体用于:
根据车辆的车速向前等间隔时长的选取N个纵向位置点,通过公式:xi=v×i×tLA(i=0,...,N-1)计算第i个纵向位置点与车辆坐标系原点之间的纵向距离xi,v为车辆的车速,tLA为设定的时间间隔时长;
通过公式:
E0Left(xi)=-min[|L0Left(xi)|,|K0Left(xi)|,|G0Left(xi)|,|T0Left(xi)|](i=0,...,N-1)确定与第i个纵向位置点对应的左侧横向位置点E0Left(xi);
通过公式:
E0Right(xi)=min[|L0Right(xi)|,|K0Right(xi)|,|G0Right(xi)|,|T0Right(xi)|](i=0,...,N-1)确定与第i个纵向位置点对应的右侧横向位置点E0Left(xi);
使用最小二乘拟合法将取得的车辆所在车道的左侧边缘的N个位置点(xi,E0Left(xi))拟合为车辆所在车道的左侧边缘曲线E0Left(x);
使用最小二乘拟合法将取得的车辆所在车道的右侧边缘的N个位置点(xi,E0Right(xi))拟合为车辆所在车道的右侧边缘曲线E0Right(x);
对左侧边缘曲线和右侧边缘曲线的各系数进行均值处理,通过公式:
Figure GDA0004001832220000111
获得使车辆沿着车道中心线行驶的安全路径方程C0(x);
通过公式:
Lw=min(Lw(x0),...,Lw(xN-1),Lwmin(xLB));
计算车辆当前所在车道的车道宽度Lw;Lwmin(xLB)表示在前tLB时间内的历史最小车道宽度,Lw(xi)=|E0Right(xi)-E0Left(xi)|(i=0,...,N-1)。
优选地,第一确定单元包括:
第一计算子单元,用于通过公式:
Figure GDA0004001832220000121
计算以车道宽度Lw和横向偏置距离C0为基础的故障函数;Lw为车辆当前所在车道的车道宽度,系数C0为车道中心线相对于车辆坐标系原点的横向偏置距离;
第一确定子单元,用于若故障函数F(Δ)小于第一设定比值D1;则确定车辆在当前所在车道内具体所处的偏置区域为第一偏置区域Z1;
第二确定子单元,用于若故障函数F(Δ)大于或等于第一设定比值D1且小于第二设定比值D2;则确定车辆在当前所在车道内具体所处的偏置区域为第二偏置区域Z2;
第三确定子单元,用于若故障函数F(Δ)大于或等于第二设定比值D2;则确定车辆在当前所在车道内具体所处的偏置区域为第三偏置区域Z3。
优选地,第二确定模块包括:
第二确定单元,用于若冗余控制子系统在上一控制周期输出的控制结果为激活第一等级冗余控制功能,或者,若冗余控制子系统在上一控制周期输出的控制结果为激活第二等级冗余控制功能或未激活冗余控制功能且车辆在当前所在车道内具体所处的偏置区域为第三偏置区域Z3,则确定主控制子系统当前的故障等级为高故障等级;
第三确定单元,包括:比对子单元,用于若冗余控制子系统在上一控制周期输出的控制结果为激活第二等级冗余控制功能或未激活冗余控制功能且车辆在当前所在车道内具体所处的偏置区域为第二偏置区域Z2,则进一步比对车辆姿态的第一判别函数G(S)与第二设定比值D2的大小;第四确定子单元,用于在车辆姿态的第一判别函数G(S)大于或等于第二设定比值D2时,确定主控制子系统当前的故障等级为高故障等级;第五确定子单元,用于在车辆姿态的第一判别函数G(S)小于第二设定比值D2时,确定主控制子系统当前的故障等级为中故障等级;
第四确定单元,用于若冗余控制子系统在上一控制周期输出的控制结果为激活第二等级冗余控制功能或未激活冗余控制功能且车辆在当前所在车道内具体所处的偏置区域为第一偏置区域Z1,则确定主控制子系统当前的故障等级为低故障等级。
优选地,比对子单元具体用于:
通过公式:
G(S)=λS
计算获得车辆姿态的判别函数G(S);λ为识别信度,λ为常数或者基于车辆的车速v、车辆的偏航率Yr、方向盘转角θSteer、车道中心线曲率2C2以及前视摄像头识别到的车辆当前所在车道的路面类型所对应的路面附着系数
Figure GDA0004001832220000131
进行RBF神经网络拟合得到;当λ为常数时λ=0.7;当λ通过RBF神经网络拟合得到时,λ的取值范围为:λ∈[0.35,1];
S为滑模变量,通过公式:
Figure GDA0004001832220000132
计算获得;滑模超平面方程满足:
Figure GDA0004001832220000133
Lw为车辆当前所在车道的车道宽度;C0为车道中心线相对于车辆坐标系原点的横向偏置距离;C1为车辆的航向角;δ为对车道宽度的一半
Figure GDA0004001832220000134
所取的横向偏置距离安全系数,δ∈[0.5,0.9];θ为车辆偏离车道中心线行驶时,使车辆渐进的回到第一偏置区域Z1内行驶的安全航向角;安全航向角θ通过公式:
Figure GDA0004001832220000135
计算获得;θ∈[θminmax],θmax为保证车辆渐进的趋近平衡点的安全上限,θmin为保证车辆渐进的趋近平衡点的安全下限;aymax为车辆最大安全侧向加速度;车辆最大安全侧向加速度aymax通过公式:
Figure GDA0004001832220000136
计算获得;T为车辆的轮距,h为车辆的质心高度,g为重力加速度,
Figure GDA0004001832220000137
为前视摄像头识别到的车辆当前所在车道的路面类型对应的路面附着系数。
优选地,仲裁模块包括:
第一激活单元,用于若主控制子系统当前的故障等级为高故障等级,则激活第一等级冗余控制功能;
第二激活单元,用于若主控制子系统当前的故障等级为中故障等级,则激活第二等级冗余控制功能;
关闭单元,用于若主控制子系统当前的故障等级为低故障等级,则不激活第一等级冗余控制功能或第二等级冗余控制功能。
优选地,仲裁模块还包括:
第一控制单元,用于若激活第一等级冗余控制功能,则按照计算出的综合减速度a控制车身电子稳定系统ESP减速至第一目标减速度,并按照计算出的目标方向盘转角θsteer和最大安全力矩Tmax控制电动助力转向系统EPS在柔性控制时间tSIL内进行方向盘转角控制;
第二控制单元,用于若激活第二等级冗余控制功能,则根据车辆姿态的第二判别函数H(S)计算出的安全减速度a1和主控制子系统输出的减速度a2控制车身电子稳定系统ESP减速至第二目标减速度;
综合减速度a为根据车辆姿态的第二判别函数H(S)计算出的安全减速度a1和主控制子系统输出的减速度a2中的较大值;第一目标速度与当前行驶环境工况相关;第二目标速度与主控制子系统输出的目标速度相关,且第二目标速度小于主控制子系统输出的目标速度。
优选地,第一控制单元具体用于:
通过公式:
Figure GDA0004001832220000141
计算目标方向盘转角θsteer;2C2为车道中心线曲率;C1为车辆的航向角;Lw为车辆当前所在车道的车道宽度;Kθδ为车辆的转向系统传动比;
Figure GDA0004001832220000142
为自适应参数;l为车辆的轴距;v为车辆的车速;ε为自适应参数初始值,ε为常数;δ为对车道宽度的一半/>
Figure GDA0004001832220000143
所取的横向偏置距离安全系数,δ∈[0.5,0.9];S为滑模变量,通过公式/>
Figure GDA0004001832220000144
计算获得,滑模超平面方程满足:/>
Figure GDA0004001832220000145
θ为车辆偏离车道中心线行驶时,使车辆渐进的回到第一偏置区域Z1内行驶的安全航向角;安全航向角θ通过公式:
Figure GDA0004001832220000151
计算获得;θ∈[θminmax],θmax为保证车辆渐进的趋近平衡点的安全上限,θmin为保证车辆渐进的趋近平衡点的安全下限;aymax为车辆最大安全侧向加速度;车辆最大安全侧向加速度aymax通过公式:
Figure GDA0004001832220000152
计算获得;T为车辆的轮距,h为车辆的质心高度,g为重力加速度,
Figure GDA0004001832220000153
为车辆传感系统识别到的车辆当前所在车道的路面类型对应的路面附着系数;
自适应参数
Figure GDA0004001832220000154
通过公式:
Figure GDA0004001832220000155
计算获得;
C0为车道中心线相对于车辆坐标系原点的横向偏置距离;
Figure GDA0004001832220000156
为死区为[0,η]的死区函数,[0,η]函数为零;Yr为车辆的偏航率;δ为对车道宽度的一半/>
Figure GDA0004001832220000157
所取的横向偏置距离安全系数,δ∈[0.5,0.9];μ为防止车辆的相对偏航率过大所设置的安全偏航率参数,K为自适应速率参数,γmax为自适应上界。
优选地,第一控制单元具体用于:
通过公式:
Figure GDA0004001832220000158
计算目标方向盘转角θsteer,2C2为车道中心线曲率;C1为车辆的航向角;Lw为车辆当前所在车道的车道宽度;Kθδ为车辆的转向系统传动比;
Figure GDA0004001832220000161
为自适应参数;l为车辆的轴距;v为车辆的车速;ε为自适应参数初始值,ε为常数;δ对车道宽度的一半/>
Figure GDA0004001832220000162
所取的横向偏置距离安全系数,δ∈[0.5,0.9];S为滑模变量,通过公式/>
Figure GDA0004001832220000163
计算获得,滑模超平面方程满足:/>
Figure GDA0004001832220000164
θ为车辆偏离车道中心线行驶时,使车辆渐进的回到第一偏置区域Z1内行驶的安全航向角;安全航向角θ通过公式:
Figure GDA0004001832220000165
计算获得;θ∈[θminmax],θmax为保证车辆渐进的趋近平衡点的安全上限,θmin为保证车辆渐进的趋近平衡点的安全下限;aymax为车辆最大安全侧向加速度;车辆最大安全侧向加速度aymax通过公式:
Figure GDA0004001832220000166
计算获得;T为车辆的轮距,h为车辆的质心高度,g为重力加速度,
Figure GDA0004001832220000167
为前视摄像头识别到的车辆当前所在车道的路面类型对应的路面附着系数;
自适应参数
Figure GDA0004001832220000168
通过公式:
Figure GDA0004001832220000169
计算获得;
C0为车道中心线相对于车辆坐标系原点的横向偏置距离;
Figure GDA0004001832220000171
为死区为[0,η]的死区函数,[0,η]函数为零;Yr为车辆的偏航率;δ为对车道宽度的一半/>
Figure GDA0004001832220000172
所取的横向偏置距离安全系数,δ∈[0.5,0.9];μ为防止车辆的相对偏航率过大所设置的安全偏航率参数,K为自适应速率参数,γmax为自适应上界;
饱和函数为
Figure GDA0004001832220000173
类饱和函数为/>
Figure GDA0004001832220000174
Figure GDA0004001832220000175
是任意连续函数,满足f(0)=0,f(1)=1,f(-1)=-1;φ为准滑模边界层宽度,φ为常数。
优选地,第一控制单元和第二控制单元具体用于:
通过公式:
Figure GDA0004001832220000176
计算出安全减速度a1
g为重力加速度;
Figure GDA0004001832220000177
为前视摄像头识别到的车辆当前所在车道的路面类型估计出的对应的路面附着系数;a1max为车辆的最大安全减速度,/>
Figure GDA0004001832220000178
a1min为车辆的最小安全减速度,a1min为常数;tD1为从车辆从第一偏置区域Z1进入第二偏置区域Z2的时间点;tNow为当前时间点;Yr为车辆的偏航率;v为车辆的车速;
σ为根据车辆的转向能力和底盘稳定能力所确定的减速度系数,通过公式:
Figure GDA0004001832220000179
计算获得;LTRCrit为安全临界荷载偏移率,LTRCrit为常数;LTRavg为平均荷载偏移率,通过公式:/>
Figure GDA00040018322200001710
计算获得,tnow为当前时间点,tCrit为判定时间,T为车辆的轮距,h为车辆的质心高度;
α为故障估计量化积分参数,通过公式:
Figure GDA0004001832220000181
计算获得,ΔθCrit为安全临界航向角增量,ΔθCrit常数;
S为滑模变量,通过公式:
Figure GDA0004001832220000182
计算;C0表示车道中心线相对于车辆坐标系原点的横向偏置距离;C1为车辆的航向角;2C2为车道中心线曲率;Lw为车辆当前所在车道的车道宽度;δ为对车道宽度的一半/>
Figure GDA0004001832220000183
所取的横向偏置距离安全系数,δ∈[0.5,0.9];θ为车辆偏离车道中心线行驶时,使车辆渐进的回到第一偏置区域Z1内行驶的安全航向角;安全航向角θ通过公式:
Figure GDA0004001832220000184
计算获得;θ∈[θminmax],θmax为保证车辆渐进的趋近平衡点的安全上限,θmin为保证车辆渐进的趋近平衡点的安全下限;aymax为车辆最大安全侧向加速度;车辆最大安全侧向加速度aymax通过公式:
Figure GDA0004001832220000185
计算获得;T为车辆的轮距,h为车辆的质心高度,g为重力加速度,
Figure GDA0004001832220000186
为车辆传感系统识别到的车辆当前所在车道的路面类型对应的路面附着系数。
优选地,第一控制单元具体用于:
通过公式:
Tmax=min(TeRequire,Te-Limit)计算最大安全力矩Tmax
第一力矩TeRequire通过根据车辆的车速v从车速与第一力矩TeRequire的预设关系表查表获得,第二力矩Te-Limit由电动助力转向系统EPS输出;
通过公式:
Figure GDA0004001832220000191
计算得到柔性控制时间tSIL,且所述柔性控制时间tSIL满足:2s≥tSIL≥100ms。
优选地,第一控制单元具体用于:
通过sigmoid函数插值,确定不同插值时刻t对应的插值方向盘转角指令值θSteerInterp;通过公式:
Figure GDA0004001832220000192
计算插值方向盘转角指令值θSteerInterp,插值时间t∈[0,tSIL];θSteerDifference为第一等级冗余控制功能激活的时间点对应的车辆的方向盘转角实际反馈值θSteerFB与目标方向盘转角θsteer之差;
将在柔性控制时间tSIL内的每一插值时刻t对应的插值方向盘转角θSteerInterp输出给电动助力转向系统EPS,使电动助力转向系统EPS在所述最大安全力矩Tmax条件下将车辆的方向盘转角由方向盘转角实际反馈值θSteerFB调节为所述插值方向盘转角指令值θSteerInterp
优选地,所述仲裁模块还包括:
第一判断单元,用于针对自动驾驶系统为L4或L5高等级自动驾驶系统,在按照第一等级冗余控制功能或第二等级冗余控制功能进行车辆控制时,判断是否满足:车辆当前所处的偏置区域稳定为第一偏置区域Z1、车辆的车速是否小于车辆当前所在车道的最低限制车速且车辆当前所在车道的路况为良好路况;
第一输出单元,用于若满足,输出提示驾驶员进行车辆接管的提示信息;
第二判断单元,用于判断驾驶员是否接管车辆;
第一关闭单元,用于若接管,则关闭第一等级冗余控制功能或第二等级冗余控制功能,并停止向电动助力转向系统EPS和车身电子稳定系统ESP输出控制指令;
第二输出单元,用于针对自动驾驶系统为L2.5或L3中等级自动驾驶系统,在按照第一等级冗余控制功能或第二等级冗余控制功能进行车辆控制时,输出提示驾驶员进行车辆接管的提示信息;
第三判断单元,用于判断驾驶员是否接管车辆;
第二关闭单元,用于若接管,则关闭第一等级冗余控制功能或第二等级冗余控制功能,并停止向电动助力转向系统EPS和车身电子稳定系统ESP输出控制指令。
本发明实施例还提供了一种自动驾驶汽车,包括上述的自动驾驶系统。
本发明实施例还提供了一种控制器,包括存储器和处理器,所述存储器中存储有至少一条指令、至少一段程序、代码集或指令集,所述至少一条指令、至少一段程序、代码集或指令集由所述处理器加载并执行,以实现如上述的冗余控制方法。
本发明实施例还提供了一种计算机可读存储介质,所述存储介质上存储有至少一条指令、至少一段程序、代码集或指令集,所述至少一条指令、至少一段程序、代码集或指令集由所述处理器加载并执行,以实现如上述的冗余控制方法。
本发明的有益效果为:
当车身外侧和车道的边缘线接近重合(处于Z2区域),轮胎离车道的边缘线仅剩下一个安全距离裕量;或者车身处于(处于Z1区域),车身姿态不佳时(G(S)≥D2);通过激活冗余控制功能,执行冗余控制。本方案中,执行冗余控制的具体方式为:基于根据自适应率和滑模趋近率,渐进的增加方向盘转角,补偿方向为车道内侧,使车辆朝向靠近车道中心线一侧运动;同时协同纵向控制来降低车速,提高车辆横向控制过程中的稳定性。
对于中等级自动驾驶车辆,当故障、错误或扰动未超过设计容错上界,本方案中的冗余控制功能可以使将车辆保持在安全区域内行驶;当故障、错误或扰动超过设计容错上界,车身可能会偏离出车道边缘线,此时,本实施例中的冗余控制功能可以延长车辆和侧方护栏、路沿或车辆的碰撞时间,为驾驶员接管自动驾驶车辆延长时间和创造条件。
对于等级自动驾驶车辆,安全员不太可能全程注意力在环或完全无安全员,当故障、错误或扰动未超过设计容错上界,冗余控制功能可以使车辆进行紧急车道保持,控制车辆保持在安全区域内行驶。当故障,错误或扰动超过设计容错上界,通过本方案中的冗余控制功能使得车辆保持在安全区域内,或从法律和法规上处于有利位置(即车辆未随意变道,不会与其他车道内车辆碰撞,同时,车辆已执行减速可避免碰撞前方目标)。
附图说明
图1为本发明实施例的自动驾驶系统的构成图;
图2为本发明实施例中用于说明本实施例的冗余控制功能效果的图;
图3是用于说明偏置区域和故障等级的图;
图4是实施例中应用于中等级自动驾驶系统的方法详细流程图;
图5是实施例中应用于高等级自动驾驶系统的方法详细流程图;
图6为本发明中等级自动驾驶系统的冗余控制子系统的模块框图;
图7为本发明高等级自动驾驶系统的冗余控制子系统的模块框图;
图8为本发明自适应滑模控制方法的相平面、滑模超平面和相轨迹示意图。
具体实施方式
首先应当说明的是,本发明实施例中的各计算公式均采用车辆重心为坐标原点,方向按照SAE定义。在车辆坐标系中,横轴为y轴(右侧为正),纵轴为x轴(前侧为正)。偏航率顺时针方向为正,方向盘转角为左正右负。进行车辆坐标系的变更不影响本方案的计算结果的物理含义和方案的实施。
本发明实施例中提供了一种应用于高等级自动驾驶系统(4级或5级)和中等级自动驾驶系统(2.5级或3级)的冗余控制方法。
对于高等级自动驾驶系统和中等级自动驾驶系统来说,都包括有如图1中的主控制子系统和冗余控制子系统。当自动驾驶系统的自动驾驶功能被激活时,主控制子系统即按照既定设计策略运行,主控制子系统未故障时,仅依靠主控制子系统就能控制车辆进行车道保持路径/轨迹跟踪。当主控制子系统发生异常时,则由冗余控制子系统继续执行自动控制。具体来说,在冗余控制子系统检测到因为传感器识别、主控制子系统的控制器计算引起的有界的故障时而引起的车辆偏移时,冗余控制子系统能够在主控制子系统故障或错误的情况下,补偿有界的故障并输出正确的目标方向盘转角和减速度,使得车辆能够紧急保持车道并执行减速。
对于本实施例中的主控制子系统来说,其执行的自动控制逻辑同现有技术一致,具体来说:主控制子系统基于作为车道保持的主要传感器(前视摄像头和雷达等)识别到的信息拟合出基于车道中心线的三阶曲线方程作为车辆行驶路径,再基于拟合出的车辆行驶路径进行车辆横向控制;以及根据车载雷达输出的目标车辆信息,识别威胁车辆行驶的目标车辆,并基于识别结果进行车辆纵向控制。主控制子系统进行自动驾驶控制的原理与现有技术中的技术原理一致。
本发明实施例提供了一种针对中等级或高等级自动驾驶系统的冗余控制方法,该方法通过自动驾驶系统中新增的冗余控制子系统来实现,本方法设计目的在于在传感器对车道中心线曲率检测不准确、或在自动驾驶系统的主控制子系统部分出现故障且故障未超过容限时,仍然能够将车辆控制在安全区域(本实施例中的第一偏置区域)内行驶,进而实现尽可能地保护驾驶员和车辆;以及在不可避免的事故中,使得车辆仍然能够保持在车辆当前所在车道内,并避免追尾前方车辆,从法律和法规上处于有利位置(原因在于:车辆并未随意变道,因而不会与其他车道内的车辆碰撞;同时由于车辆已执行纵向减速尽可能避免碰撞前方目标)。
本实施例中提供的该方法仅适用于自动驾驶系统的车道保持功能,不适用于自动驾驶系统中的主控制子系统因用户动作触发自动变道和主控制子系统控制自主自动变道功能。在车辆在自动变道过程中,本方法所属的冗余控制功能处于关闭状态。
在图4的实施例中,该方法应用于搭载中等级自动驾驶(2.5级/3级)系统的车辆中。中等级自动驾驶系统采用人机共驾模式,驾驶员作为自动驾驶车辆的安全责任主体,因此在冗余控制方面需要优先考虑驾驶员的接管(即在执行冗余控制的过程中,需要直接提醒驾驶员进行车辆接管)。
结合参照图1和图4来说,本实施例中的应用于中等级自动驾驶系统的该方法具体执行下面的步骤。
S0步:上电初始化,进入S1步。
S1步:冗余控制子系统接收自动驾驶系统的状态信号,确定自动驾驶系统的自动驾驶功能是否激活。具体来说,S1步包括:
S1.1步,冗余控制子系统接收自动驾驶系统的工作状态信号;
S1.2步,冗余控制子系统确定自动驾驶系统的自动驾驶功能是否激活,如激活则进入S2步。如没有激活,则进入S1.3步。
S1.3步,判定整车是否下电,如下电结束程序,否则循环执行S1步。
S2步,冗余控制子系统接收自动驾驶系统的各输入传感器(前视摄像头,雷达等)所输出的多源信息(主要包括车道线曲线方程,路沿曲线方程、护栏曲线和车流曲线),再进入S3步。各传感器所采集的信息进行一些初步的处理等步骤作为现有技术常见的方式,本实施例1中同样也需要这些处理方式,不对此详细描述。
S3步:执行S3.1和S3.2,完成后进入S4步。
S3.1步,通过前视摄像头和雷达输出的多源信息,即车道线曲线方程L0(x),路沿曲线方程K0(x)、护栏曲线G0(x)和车流曲线T0(x),计算位于车辆所在车道的左侧边缘曲线方程和右侧边缘曲线方程。该步骤具体包括:
根据车速向前等间隔时长的选取N(N≥5)个纵向位置点,通过公式:
xi=v×i×tLA(i=0,...,N-1)
计算第i个纵向位置点与车辆坐标系原点之间的纵向距离xi;其中v为车辆的车速,tLA为设定的时间间隔时长(tLA为常数)。
根据前视摄像头和雷达识别到的车道线曲线方程L0(x),路沿曲线方程K0(x)、护栏曲线K0(x)和车流曲线方程T0(x),通过公式:
E0Left(xi)=-min[|L0Left(xi)|,|K0Left(xi)|,|G0Left(xi)|,|T0Left(xi)|](i=0,...,N-1)确定与第i个纵向位置点对应的左侧横向位置点E0Left(xi);
通过公式:
E0Right(xi)=min[|L0Right(xi)|,|K0Right(xi)|,|G0Right(xi)|,|T0Right(xi)|](i=0,...,N-1)确定与第i个纵向位置点对应的右侧横向位置点E0Left(xi);
使用最小二乘拟合法将取得的车辆所在车道的左侧边缘的N个位置点(xi,E0Left(xi))拟合为车辆所在车道的左侧边缘曲线E0Left(x),其中最小二乘拟合法为现有数值计算技术;
使用最小二乘拟合法将取得的车辆所在车道的右侧边缘的N个位置点(xi,E0Right(xi))拟合为车辆所在车道的右侧边缘曲线E0Right(x),其中最小二乘拟合法为现有数值计算技术;
对左侧边缘曲线和右侧边缘曲线的各系数进行均值处理,通过公式:
Figure GDA0004001832220000221
获得使车辆沿着车道中心线行驶的安全路径方程C0(x);
C0为车辆所在车道的车道中心线相对于车辆坐标系原点的横向偏置距离,C1为车辆所在车道的车道中心线相对车辆坐标系原点的夹角(即车辆的航向角),2C2为车道中心线曲率,6C3为车道中心线曲率变化率,x为以车辆重心为坐标原点的车辆直角坐标系的纵向坐标。
S3.2步,通过左右道路边缘曲线和历史车道宽度,计算车辆当前所在车道的车道宽度Lw,具体通过前N(N≥5)个边沿曲线点,计算N个的道路宽度,
Lw(xi)=E0Right(xi)-E0Left(xi)(i=0,...,N-1)
通过前N个位置点处的道路宽度和前tLB(tLB为设定时间常数)时间内的历史最小车道宽度Lwmin(xLB)共N+1个车道宽度取最小值,得到道路宽度Lw:
Lw=min(Lw(x0),...,Lw(xN-1),Lwmin(xLB))
S4步,根据规划出的安全路径方程C0(x)的系数C0和车道宽度Lw,识别车辆当前所处的目标偏置区域。该步骤具体执行S4.1步和S4.2步。
S4.1步,根据车道宽度Lw,以及预先设定的各偏置区域与车道宽度的设定百分比,划分各偏置距离。如图2和3,本实施例中:第一偏置区域Z1是指在左侧方向和右侧方向各自与车道中心线之间的距离小于或等于M1所围成的区域,第二偏置区域Z2是指在左侧方向和右侧方向各自与车道中心线之间的距离小于或等于M2且大于M1所围成的区域;第三偏置区域Z3是指在左侧方向和右侧方向各自与车道中心线之间的距离大于M2且小于或等于
Figure GDA0004001832220000231
所围成的区域。
S4.2步,计算以车道宽度Lw和横向偏置距离C0为基础的故障函数F(Δ)如下:
Figure GDA0004001832220000232
其中Lw为车道宽度,C0为车道中心线相对于车辆坐标系原点的横向偏置距离。
S4.2步,根据判别函数F(Δ)识别车辆具体所处的目标偏置区域,
1)、车辆位于第一偏置区域Z1(见图3)内需满足:
F(Δ)=|Δ|<D1
D1为第一设定比值,具体为百分比值,为(预先设置的使车辆保持车道中心线行驶的最大横向偏差M1/0.5倍车道宽度)×100%。
2)、车辆位于第二偏置区域Z2(见图3)内需要满足:
D1≤F(Δ)=|Δ|<D2
D2为第二设定比值,具体为百分比值,为(预先设定的车辆保持与侧方车道/护栏/路沿横向安全距离的最大横向偏差M2/0.5倍车道宽度)×100%。
3)、车辆位于第三偏置区域Z3(见图3)内需要满足:
F(Δ)=|Δ|≥D2
S5步,根据车辆当前具体所处的偏置区域Z1,Z2,Z3,以及冗余控制子系统在上一控制周期输出的控制结果(具体为冗余控制子系统在上一控制周期是否激活第一等级冗余控制功能),来判定主控制子系统的故障等级。该步骤S5具体包括分支S5B1步-S5B4步。
S5B1步,冗余控制子系统确定在上一控制周期是否激活第一等级冗余控制功能,如激活,则确定主控制子系统的故障等级为高故障等级,进一步进入S6B1。
S5B2步,如冗余控制子系统确定在上一控制周期未激活第一等级冗余控制功能但车辆当前处于第三偏置区域Z3,同样确定主控制子系统的故障等级为高故障等级,进入S6B2步。
S5B3步,如冗余控制子系统确定在上一控制周期未激活第一等级冗余控制功能但车辆当前处于第二偏置区域Z2,则需要进一步计算车辆姿态的第一判别函数G(S),通过公式:
Figure GDA0004001832220000241
计算该第一判别函数,C0为车道中心线相对于车辆坐标系原点的横向偏置距离;C1为车道中心线与车辆的夹角(即航向角)。δ为车辆的横向偏置距离安全系数,δ具体为对车道宽度一半的
Figure GDA0004001832220000242
所取横向偏置距离安全系数,一般取δ∈(0.5,0.9),当整车和转向系统的制造一致性较差时,δ应适当减小。θ为车辆偏离车道中心线行驶时,使车辆渐进的回到第一偏置区域Z1内行驶的安全航向角;安全航向角θ通过公式:
Figure GDA0004001832220000243
计算获得;θmax为保证车辆渐进的趋近平衡点的安全上限,θmin为保证车辆渐进的趋近平衡点的安全下限;aymax为车辆最大安全侧向加速度;车辆最大安全侧向加速度aymax通过公式:
Figure GDA0004001832220000244
计算获得;T为车辆的轮距,h为车辆的质心高度,g为重力加速度,
Figure GDA0004001832220000245
为车辆传感系统识别到的车辆当前所在车道的路面类型对应的路面附着系数。
λ为识别信度,在无自动驾驶大数据支撑的正常情况下,λ为常数,取λ=0.7;在有大数据支撑的情况下,可通过运营和用户数据提高识别准确率,具体来说:将车辆的车速v(由车辆的车身电子稳定系统ESP输出),车辆偏航率(由车辆的车身电子稳定系统ESP)输出)Yr,方向盘转角θSteer,车道中心线曲率2C2,根据前视摄像头识别到的路面类型确定出的对应的路面附着系数
Figure GDA0004001832220000251
(系统中预设有不同路面类型和路面附着系数的对应关系),通过RBF神经网络拟合得到关于参数λ的函数:
Figure GDA0004001832220000252
为保证拟合的安全性,本实施例1中,限定λ的范围λ∈[0.35,1]。
如车辆姿态的第一判别函数G(S)≥D2,则判定主控制子系统的故障等级为高故障等级,进一步进入S6B2;否则判定主控制子系统的故障等级为中故障等级,进一步进入S6B3步。
S5B4步,如冗余控制子系统确定在上一控制周期未激活冗余控制功能但车辆当前处于第一偏置区域Z1,则确定主控制子系统的故障等级为低故障等级,进一步进入S6B4步。
S6步:根据S5步选择执行分支S5B1-S5B4,分别执行对应S6B1-S6B4,完成后再分别进入对应S7B1-S7B4步。
S6B1步,如进入分支S6B1步,即确定主控制子系统的故障等级为高故障等级,则冗余控制子系统激活第一等级冗余控制功能,激活第一等级冗余控制功能表示需要计算:1)、使车身电子稳定系统ESP进行纵向控制的所需综合减速度a;2),使电动助力转向系统EPS进行横向控制的最大安全力矩Tmax、目标方向盘转角θsteer和柔性控制时间tSIL;在计算出这些参数信息后再进入S7B1步。
本实施例中,对于使车身电子稳定系统ESP进行纵向控制的综合减速度a的具体计算方式如下:
1)、根据车辆姿态的第一判别函数H(S),通过下列公式计算出安全减速度a1
Figure GDA0004001832220000253
式中,g为重力加速度。
Figure GDA0004001832220000254
为前视摄像头识别到的路面类型估计出的对应的路面附着系数。/>
Figure GDA0004001832220000255
为车辆的最大安全减速度。a1min=0.5m/s,表示车辆的最小安全减速度。
σ为根据车辆的转向能力和底盘稳定能力所确定的减速度系数,计算公式:
Figure GDA0004001832220000261
LTRCrit为安全临界荷载偏移率(LTRCrit为常数)。LTRavg为平均荷载偏移率:/>
Figure GDA0004001832220000262
tnow为当前时间点,tCrit为判定时间,T为轮距,h为车辆质心高度。
α为故障估计量化积分参数,计算公式:
Figure GDA0004001832220000263
ΔθCrit为安全临界航向角增量(ΔθCrit为常数)。tD1为从车辆从第一偏置区域Z1进入第二偏置区域Z2的时间点;tNow为当前时间点;Yr为车辆的偏航率;v为车辆的车速。
S为滑模变量,通过公式:
Figure GDA0004001832220000264
计算,C0表示车道中心线相对于车辆坐标系原点的横向偏置距离;C1为车辆的航向角;2C2为车道中心线曲率;Lw为车辆当前所在车道的车道宽度;δ为对车道宽度的一半/>
Figure GDA0004001832220000265
所取的横向偏置距离安全系数,δ∈[0.5,0.9];θ为车辆偏离车道中心线行驶时,使车辆渐进的回到第一偏置区域Z1内行驶的安全航向角;安全航向角θ通过公式:/>
Figure GDA0004001832220000266
计算获得;θ∈[θminmax],θmax为保证车辆渐进的趋近平衡点的安全上限,θmin为保证车辆渐进的趋近平衡点的安全下限;aymax为车辆最大安全侧向加速度;车辆最大安全侧向加速度aymax通过公式:
Figure GDA0004001832220000267
计算获得;T为车辆的轮距,h为车辆的质心高度,g为重力加速度,
Figure GDA0004001832220000271
为前视摄像头识别到的车辆当前所在车道的路面类型对应的路面附着系数。
2)、获得主控制子系统在当前控制周期输出的减速度a2,该减速度a2用于保证不和前方车辆道目标碰撞。主控制子系统对减速度a2的计算方式为现有技术。
3)、计算综合减速度a,在上述减速度a1和a2中取其中绝对值大的,综合减速度a不超过系统最大能力。综合减速度a满足:
a=-max(|a1|,|a2|)
此外,本实施例1中,对于使电动助力转向系统EPS执行横向控制所需的目标方向盘转角θsteer的具体计算方式为:
Figure GDA0004001832220000272
Figure GDA0004001832220000273
Lw为车辆当前所在车道的车道宽度,v为车辆的车速,δ为车辆的横向偏置距离安全系数,θ为车辆偏离车道中心线时,使车辆渐进的回到车道中心线附近Z1区行驶的安全航向角;Kθδ为车辆的转向系统传动比,/>
Figure GDA0004001832220000274
为自适应参数,l为车辆的轴距;ε为滑模自适应参数初始值,ε为常数。
本实施例中,采用自适应滑模控制原理来计算所需的目标方向盘转角θsteer。首先,定义使得电动助力转向系统EPS进行横向控制所需的滑模函数S(C0,C1)为:
Figure GDA0004001832220000275
如图8,滑模超平面满足:
Figure GDA0004001832220000276
/>
上式中,Lw为车辆当前所在车道的车道宽度,δ为车辆的横向偏置距离安全系数,θ为车辆的安全航向角。
δ具体为对车道宽度一半的
Figure GDA0004001832220000277
所取的横向偏置距离安全系数,一般取δ∈(0.5,0.9),当整车和转向系统的制造一致性较差时,δ应适当减小。当δ设计确定后,控制车辆稳定在滑模超平面上时应满足:/>
Figure GDA0004001832220000281
因此,θ(安全航向角)决定了当横向偏置距离C0一定时的航向角C1。例如,如果车道中心线相对于车辆坐标系原点的横向偏置距离为/>
Figure GDA0004001832220000282
由/>
Figure GDA0004001832220000283
求解可得到航向角C1=-θ。
进一步地,定义横向控制李雅普诺夫函数J:
Figure GDA0004001832220000284
根据滑模控制稳定性理论,横向控制稳定,李雅普诺夫函数J则满足:
Figure GDA0004001832220000285
进一步地,采用自适应滑模趋近律
Figure GDA0004001832220000286
Figure GDA0004001832220000287
根据滑模控制理论,
Figure GDA0004001832220000288
越大,自适应滑模趋近律/>
Figure GDA0004001832220000289
的容错能力和抗扰动抗摄动能力越强,因此用于估计自适应参数/>
Figure GDA00040018322200002810
的一种自适应率方法为:
Figure GDA00040018322200002811
上面自适应率中的自适应参数
Figure GDA00040018322200002812
的初值可取大于零的实数,如γ0=1,式中
Figure GDA00040018322200002813
为死区为[0,η]的死区函数,[0,η]函数为零,Yr为车辆的偏航率,δ为车辆的横向偏置距离安全系数,θ为车辆的安全航向角;μ为防止车辆的相对偏航率过大所设置的安全偏航率参数(μ为常数),K为自适应速率参数(K为常数);γmax为自适应上界(γmax为常数),由传感器/控制器/执行器的故障误差总体设计容错上界确定。
想要满足
Figure GDA00040018322200002814
Figure GDA0004001832220000291
/>
得出滑模控制率(即为本实施例1中的目标方向盘转角θsteer)为:
Figure GDA0004001832220000292
此外,在本实施例中,对发送给电动助力转向系统EPS进行纵向控制所需要的最大安全力矩Tmax的计算公式如下:
Tmax=min(TeRequire,Te-Limit)
Te-Limit为电动助力转向系统EPS根据直流电压、最大电流、电机温度、变换器温度求出的第二力矩(其计算方式为现有技术)。TeRequire为通过车速查表计算得到的第一力矩。
TeRequire=map(v)
其中查表函数map(v)通过车型实车试验数据拟合得到,如下表所示:
v(km/h) 0 20 40 60 80 100 120 140 160
T<sub>eRequire</sub>(Nm) 20 17.5 12.5 9.5 7.25 6.5 5.75 4.65 3.75
此外,本实施例1中,柔性控制时间tSIL由冗余控制子系统在当前控制周期激活冗余控制功能时刻所对应的横向安全距离(当v>1m/s,C1>0时,该横向安全距离为
Figure GDA0004001832220000293
当v>1m/s,C1<0时,该横向安全距离为/>
Figure GDA0004001832220000294
)除以横向速度|vsinC1|得到,并限制2s≥tSIL≥100ms,具体通过公式:
Figure GDA0004001832220000295
计算获得。进一步地,确定电动助力转向系统EPS在进行转向转角柔性控制(横向控制)而使车辆从当前的方向盘转角实际反馈值θSteerFB变化至上述的目标方向盘转角θsteer。这一过程中不同插值时刻t对应的插值方向盘转角θSteerInterp。本实施例1中,通过Sigmoid函数插值,来使电动助力转向系统EPS柔性启动横向控制。
其中,插值期间,不同插值时刻t对应的插值方向盘转角指令值θSteerInterp的计算公式为:
Figure GDA0004001832220000301
其中参数θsteer的计算方式为S7B1.2中的目标方向盘转角的计算方式;t∈[0,tSIL];θSteerDifference为冗余控制子系统在当前控制周期激活第一等级冗余控制功能时刻对应的方向盘转角实际反馈值θSteerFB与目标方向盘转角θsteer之差。
S6B2步,如进入分支S6B2步,即确认主控制子系统的故障等级为高故障等级,冗余控制子系统激活第一等级冗余控制功能。激活冗余控制功能表示需要计算:1)、使车身电子稳定系统ESP进行纵向控制的所需综合减速度a;2),使电动助力转向系统EPS进行横向控制的最大安全力矩Tmax、目标方向盘转角θsteer和柔性控制时间tSIL;在计算出这些参数信息后再进入S7B2步。其中综合减速度a、最大安全力矩Tmax、目标方向盘转角θsteer和柔性控制时间tSIL的计算方式与S6B1步中的方式相同。
S6B3步,如进入分支S6B3步,即确认主控制子系统的故障等级为中故障等级,则冗余控制子系统激活第二等级冗余控制功能,激活第二等级冗余控制功能需要计算安全减速度a1,然后进入S7B3步。其中安全减速度a1的计算的方式与S6B1步中的方式相同。
S6B4步,如进入分支S6B4步,确认主控制子系统的故障等级为低故障等级,冗余控制子系统不激活冗余控制功能,关闭安全减速度a1计算,然后进入S7B4步。
S7步:根据S6选择执行分支S6B1-S6B4,分别执行对应S7B1-S7B4,完成后进入S8步。
本实施例中,S7B1步,包括执行S7B1.1步和S7B1.2步。
S7B1.1步,1)、冗余控制子系统向车身电子稳定系统ESP发出纵向控制指令,具体为将综合减速度a作为最终输出给车身电子稳定系统ESP的最终加速度控制指令,使车身电子稳定系统ESP执行纵向控制。ESP车身电子稳定系统执行纵向控制的具体步骤为:执行综合减速度a至第一目标车速。其中,不同行驶环境工况条件下对应的第一目标车速不同,本实施例中,第一目标车速的确定方式包括:在非堵车工况下,若前车车速大于道路法规允许的车辆当前所在车道的最低安全车速,则执行综合减速度a至车辆当前所在车道允许的最低安全车速(此时,该最低安全车速即为第一目标车速);在车行缓慢工况下,执行综合减速度a至前车车速(该前车车速即作为第一目标车速),使本车与前车之间保持特定的纵向安全距离;在堵车工况下,若本车的当前车速在设定车速(如30Km/h)以下,则执行综合减速度a至车辆停止(此时,第一目标车速为0)。2)、冗余控制子系统向电子助力转向系统EPS发出横向控制指令,该横向控制指令中包含需要输出给电动助力转向系统EPS的基于自适应滑模控制原理计算出的目标方向盘转角θsteer,以及使电动助力转向系统EPS进行转向转角柔性控制的柔性控制时间内tSIL的不同插值时刻t对应的插值方向盘转角θSteerInterp
S7B1.2,冗余控制子系统发出提醒信号,该提醒信号作为(针对中等级自动驾驶系统)提醒驾驶员第一等级冗余控制功能激活的信号。
S7B2步,执行S7B2.1和S7B2.2。
S7B2.1,包括:
1)、向车身电子稳定系统ESP发出纵向控制指令,使得车身电子稳定系统ESP执行纵向控制,即按照步骤S7B1.1中记录的方式进行减速控制。
2)、向EPS(电子助力转向系统)发出横向控制指令,使得电动助力转向系统EPS执行横向控制。具体来说,该横向控制指令中包含需要输出给电动助力转向系统EPS的基于自适应滑模控制原理计算出的目标方向盘转角θsteer,以及使电动助力转向系统EPS进行转向转角柔性控制的柔性控制时间内tSIL的不同插值时刻t对应的插值方向盘转角θSteerInterp
S7B2.2冗余控制子系统发出提醒信号,该提醒信号作为(针对中等级自动驾驶系统)提醒驾驶员第一等级冗余控制功能激活的信号。
S7B3步,执行S7B3.1和S7B3.2。
S7B3.1步,冗余控制子系统(中等级自动驾驶)将根据计算出的安全加速度a1和主控制子系统输出的减速度a2确定输出给车身电子稳定系统ESP的最终加速度。
具体方法为,在安全减速度a1和主控制子系统输出的减速度a2中取最小值,将此作为最终输出给车身电子稳定系统ESP的最终减速度,仲裁模块进行减速控制的第二目标速度为Dv=70%(Dv为设定百分比)乘以当前主控制子系统输出的目标速度。
S7B3.2步,冗余控制子系统输出信号提醒驾驶员集中注意力,观察道路和车辆状态。
S7B4步,执行S7B4.1和S7B4.2。
S7B4.1冗余控制子系统关闭冗余控制功能。
S7B4.2冗余控制子系统(中等级自动驾驶)向主控制子系统输出由主控制子系统进行控制的主功能控制信号,使主控制子系统按照既定策略进行车辆自动驾驶控制。
S8步(中等级自动驾驶),判定驾驶员是否接管车辆。
S9步(中等级自动驾驶),驾驶员如未接管车辆,冗余控制子系统继续执行冗余控制功能直至驾驶员接管。驾驶员如已接管,则关闭冗余控制功能,停止向执行机构,电动助力转向系统EPS或车身电子稳定系统ESP输出任何控制指令。
本发明的该实施例中提供的方法,冗余控制子系统中的一部分功能模块(拟合模块和冗余模块)持续性的与主控制子系统一起工作,这部分功能模块主要用于来识别车辆当前所在的车辆偏置区域。冗余控制子系统中的另一部分功能模块(仲裁模块)根据前述模块的识别结果以及其自身在上一控制周期是否激活第一等级冗余控制功能来确定主控制子系统当前的故障等级,进而选择性的激活或不激活本实施例中的冗余控制功能。在激活第一等级冗余控制功能时,通过控制电动助力转向系统EPS进行横向控制和控制车身电子稳定系统ESP进行纵向控制;在激活第一等级冗余控制功能时,通过控制车身电子稳定系统ESP进行纵向控制;来使得车辆能够在主控制子系统出现错误或故障时仍能够保持在车道内的安全行驶区域内行驶。
本方法的安全路径规划部分,规划安全路径曲线方程时,经过了护栏,路沿,车流的多源信息校验,即使传感器和规划出的安全路径曲线方程出现错误,也极大概率使得错误的影响有界。
本方法的横向控制部分,作为输入的信息至少包括车辆的横向偏置距离、车辆航向角、车道中心线曲率、道路宽度四个输入,再基于滑模或准滑模补偿控制和自适应的控制,滑模超平面至少由横向偏置距离、车辆航向角、道路宽度三个变量得到。只要车辆横向偏移和车辆航向角的误差在一定范围内,而传感器曲率检测不准确、或者控制器、执行器或通信的故障或错误引起的方向盘转向角的错误偏角没有超过设计容错上界,那么该方法一定能干全部补偿这个故障或错误的影响;当故障,错误或扰动超过设计容错上界,该方法也能部分补偿这个故障或错误的影响。
本方法的纵向控制部分,为根据航向角、横向偏置距离等状态变量计算安全减速度,并结合前方威胁目标执行减速协同控制,利用车辆减速降低车辆侧向动能,降低发生侧滑的概率,提高车辆横向控制稳定性。在非堵车工况,前方无威胁目标时,尽可能的不采用刹停的策略,而是降低车速至高速公路车辆道最低速度限制,减小在高速公路上执行安全控制时被追尾概率,极大可能提高安全性。
如图5,本实施例中的方法执行逻辑与上述图4中的方法执行逻辑一致。区别仅在于,图5中方法应用在高等级自动驾驶系统(4级/5级)中,此时的自动驾驶系统采用全自动驾驶模式,自动驾驶感知和主控制子系统是自动驾驶车辆的安全责任主体,安全员可能注意力完全转移,而全程非在环。因此在冗余控制方面以冗余控制优先,只有在冗余控制稳定之后,安全员才能接管,其中S0至S7与图4中的方法执行步骤相同,区别在于步骤S7B1.2、S7B2.2,S7B3.2,S8和S9与实施例1中的对应步骤不同。
应用于高等级自动驾驶系统的冗余控制方法的S7B1.2和S7B2.2步为:冗余控制子系统仅发出提醒安全员第一等级冗余控制功能激活的信号,而在车辆稳定、车速降低到安全车速且路况良好前不提醒安全员进行接管。
S7B3.2步为:冗余控制子系统仅输出安全减速度a1;但不输出提醒信号去提醒安全员集中注意力。
S8步,在冗余控制子系统激活第一冗余控制功能后,在冗余控制子系统执行步骤S7B1和步骤S7B2后,先执行S8d1。
S8d1步,冗余控制子系统需要先判定自动驾驶车辆是否处于稳定安全状态,判断车辆处于稳定安全状态的具体步骤为:1)、车辆坐标系原点稳定位于第一偏置区域Z1内;2)、当前道路路况良好;3)、车辆的车速低于车辆当前所在车道允许的最低安全车速。在满足上述这些条件后,再执行步骤S8,确定是否重新启动主控制子系统进行主功能控制或请求安全员接管。
S9步,先执行S9d1。
S9d1步,安全员未接管车辆或主控制子系统未按照冗余控制子系统发送的主控制功能信号去进行车辆控制,则尝试重新向主控制子系统发送主功能控制信号请求主控制子系统激活主功能控制。执行步骤S9d1步后,在执行步骤S9。
S9步,安全员或远程控制如已接管,则冗余控制子系统关闭冗余控制功能,停止向电动助力转向系统EPS和车身电子稳定系统ESP输出控制指令。
本实施例中,目标方向盘转角θsteer还可以采用另外一种方式进行求解获得,此处滑模控制率(即目标方向盘转角θsteer的求解公式)为:
Figure GDA0004001832220000331
饱和函数为
Figure GDA0004001832220000332
类饱和函数为/>
Figure GDA0004001832220000333
Figure GDA0004001832220000334
是任意连续函数,满足f(0)=0,f(1)=1,f(-1)=-1;φ为准滑模边界层宽度(其为常数)。
参照图6和图7,本发明实施例还提供了一种自动驾驶系统,所述自动驾驶系统包括主控制子系统和冗余控制子系统,所述冗余控制子系统包括:
第一确定模块,用于在自动驾驶系统的自动驾驶功能激活时,确定车辆在当前所在车道内具体所处偏置区域;偏置区域是指车道中心线左侧方向和车道中心线右侧方向各自与车道中心线之间的距离位于预定距离范围内时所组成的区域;
第二确定模块,用于根据车辆当前所处的偏置区域和冗余控制子系统在上一控制周期输出的控制结果,确定主控制子系统当前的故障等级;
仲裁模块,用于根据主控制子系统当前的故障等级,激活对应等级的冗余控制功能,并按照对应等级的冗余控制功能进行车辆控制。
优选地,第一确定模块包括:
拟合单元,用于根据前视摄像头和雷达输出到的多源信息,拟合使车辆保持车道中心线行驶的安全路径方程C0(x),并确定车辆当前所在车道的车道宽度Lw;
第一确定单元,根据所规划的安全路径方程C0(x)中的系数C0和车道宽度Lw,确定车辆在当前所在车道内具体所处的偏置区域。
优选地,摄像头和雷达输出到的多源信息包括:相对于车辆坐标系原点规划的车道线曲线方程L0(x)、路沿曲线方程K0(x)、护栏曲线方程G0(x)和/或车流曲线方程T0(x);拟合单元具体用于:
根据车辆的车速向前等间隔时长的选取N个纵向位置点,通过公式:xi=v×i×tLA(i=0,...,N-1)计算第i个纵向位置点与车辆坐标系原点之间的纵向距离xi,v为车辆的车速,tLA为设定的时间间隔时长;
通过公式:
E0Left(xi)=-min[|L0Left(xi)|,|K0Left(xi)|,|G0Left(xi)|,|T0Left(xi)|](i=0,...,N-1)确定与第i个纵向位置点对应的左侧横向位置点E0Left(xi);
通过公式:
E0Right(xi)=min[|L0Right(xi)|,|K0Right(xi)|,|G0Right(xi)|,|T0Right(xi)|](i=0,...,N-1)确定与第i个纵向位置点对应的右侧横向位置点E0Left(xi);
使用最小二乘拟合法将取得的车辆所在车道的左侧边缘的N个位置点(xi,E0Left(xi))拟合为车辆所在车道的左侧边缘曲线E0Left(x);
使用最小二乘拟合法将取得的车辆所在车道的右侧边缘的N个位置点(xi,E0Right(xi))拟合为车辆所在车道的右侧边缘曲线E0Right(x);
对左侧边缘曲线和右侧边缘曲线的各系数进行均值处理,通过公式:
Figure GDA0004001832220000341
获得使车辆沿着车道中心线行驶的安全路径方程C0(x);
通过公式:
Lw=min(Lw(x0),...,Lw(xN-1),Lwmin(xLB));
计算车辆当前所在车道的车道宽度Lw;Lwmin(xLB)表示在前tLB时间内的历史最小车道宽度,Lw(xi)=|E0Right(xi)-E0Left(xi)|(i=0,...,N-1)。
优选地,第一确定单元包括:
第一计算子单元,用于通过公式:
Figure GDA0004001832220000342
计算以车道宽度Lw和横向偏置距离C0为基础的故障函数;Lw为车辆当前所在车道的车道宽度,系数C0为车道中心线相对于车辆坐标系原点的横向偏置距离;
第一确定子单元,用于若故障函数F(Δ)小于第一设定比值D1;则确定车辆在当前所在车道内具体所处的偏置区域为第一偏置区域Z1;
第二确定子单元,用于若故障函数F(Δ)大于或等于第一设定比值D1且小于第二设定比值D2;则确定车辆在当前所在车道内具体所处的偏置区域为第二偏置区域Z2;
第三确定子单元,用于若故障函数F(Δ)大于或等于第二设定比值D2;则确定车辆在当前所在车道内具体所处的偏置区域为第三偏置区域Z3。
优选地,第二确定模块包括:
第二确定单元,用于若冗余控制子系统在上一控制周期输出的控制结果为激活第一等级冗余控制功能,或者,若冗余控制子系统在上一控制周期输出的控制结果为激活第二等级冗余控制功能或未激活冗余控制功能且车辆在当前所在车道内具体所处的偏置区域为第三偏置区域Z3,则确定主控制子系统当前的故障等级为高故障等级;
第三确定单元,包括:比对子单元,用于若冗余控制子系统在上一控制周期输出的控制结果为激活第二等级冗余控制功能或未激活冗余控制功能且车辆在当前所在车道内具体所处的偏置区域为第二偏置区域Z2,则进一步比对车辆姿态的第一判别函数G(S)与第二设定比值D2的大小;第四确定子单元,用于在车辆姿态的第一判别函数G(S)大于或等于第二设定比值D2时,确定主控制子系统当前的故障等级为高故障等级;第五确定子单元,用于在车辆姿态的第一判别函数G(S)小于第二设定比值D2时,确定主控制子系统当前的故障等级为中故障等级;
第四确定单元,用于若冗余控制子系统在上一控制周期输出的控制结果为激活第二等级冗余控制功能或未激活冗余控制功能且车辆在当前所在车道内具体所处的偏置区域为第一偏置区域Z1,则确定主控制子系统当前的故障等级为低故障等级。
优选地,比对子单元具体用于:
通过公式:
G(S)=λ|S|
计算获得车辆姿态的判别函数G(S);λ为识别信度,λ为常数或者基于车辆的车速v、车辆的偏航率Yr、方向盘转角θSteer、车道中心线曲率2C2以及前视摄像头识别到的车辆当前所在车道的路面类型所对应的路面附着系数
Figure GDA0004001832220000351
进行RBF神经网络拟合得到;当λ为常数时λ=0.7;当λ通过RBF神经网络拟合得到时,λ的取值范围为:λ∈[0.35,1];
S为滑模变量,通过公式:
Figure GDA0004001832220000352
/>
计算获得;滑模超平面方程满足:
Figure GDA0004001832220000353
Lw为车辆当前所在车道的车道宽度;C0为车道中心线相对于车辆坐标系原点的横向偏置距离;C1为车辆的航向角;δ为对车道宽度的一半
Figure GDA0004001832220000354
所取的横向偏置距离安全系数,δ∈[0.5,0.9];θ为车辆偏离车道中心线行驶时,使车辆渐进的回到第一偏置区域Z1内行驶的安全航向角;安全航向角θ通过公式:
Figure GDA0004001832220000361
计算获得;θ∈[θminmax],θmax为保证车辆渐进的趋近平衡点的安全上限,θmin为保证车辆渐进的趋近平衡点的安全下限;aymax为车辆最大安全侧向加速度;车辆最大安全侧向加速度aymax通过公式:
Figure GDA0004001832220000362
计算获得;T为车辆的轮距,h为车辆的质心高度,g为重力加速度,
Figure GDA0004001832220000363
为前视摄像头识别到的车辆当前所在车道的路面类型对应的路面附着系数。
优选地,仲裁模块包括:
第一激活单元,用于若主控制子系统当前的故障等级为高故障等级,则激活第一等级冗余控制功能;
第二激活单元,用于若主控制子系统当前的故障等级为中故障等级,则激活第二等级冗余控制功能;
关闭单元,用于若主控制子系统当前的故障等级为低故障等级,则不激活第一等级冗余控制功能或第二等级冗余控制功能。
优选地,仲裁模块还包括:
第一控制单元,用于若激活第一等级冗余控制功能,则按照计算出的综合减速度a控制车身电子稳定系统ESP减速至第一目标减速度,并按照计算出的目标方向盘转角θsteer和最大安全力矩Tmax控制电动助力转向系统EPS在柔性控制时间tSIL内进行方向盘转角控制;
第二控制单元,用于若激活第二等级冗余控制功能,则根据车辆姿态的第二判别函数H(S)计算出的安全减速度a1和主控制子系统输出的减速度a2控制车身电子稳定系统ESP减速至第二目标减速度;
综合减速度a为根据车辆姿态的第二判别函数H(S)计算出的安全减速度a1和主控制子系统输出的减速度a2中的较大值;第一目标速度与当前行驶环境工况相关;第二目标速度与主控制子系统输出的目标速度相关,且第二目标速度小于主控制子系统输出的目标速度。
优选地,第一控制单元具体用于:
通过公式:
Figure GDA0004001832220000371
计算目标方向盘转角θsteer;2C2为车道中心线曲率;C1为车辆的航向角;Lw为车辆当前所在车道的车道宽度;Kθδ为车辆的转向系统传动比;
Figure GDA0004001832220000372
为自适应参数;l为车辆的轴距;v为车辆的车速;ε为自适应参数初始值,ε为常数;δ为对车道宽度的一半/>
Figure GDA0004001832220000373
所取的横向偏置距离安全系数,δ∈[0.5,0.9];S为滑模变量,通过公式/>
Figure GDA0004001832220000374
计算获得,滑模超平面方程满足:/>
Figure GDA0004001832220000375
θ为车辆偏离车道中心线行驶时,使车辆渐进的回到第一偏置区域Z1内行驶的安全航向角;安全航向角θ通过公式:
Figure GDA0004001832220000376
计算获得;θ∈[θminmax],θmax为保证车辆渐进的趋近平衡点的安全上限,θmin为保证车辆渐进的趋近平衡点的安全下限;aymax为车辆最大安全侧向加速度;车辆最大安全侧向加速度aymax通过公式:
Figure GDA0004001832220000377
计算获得;T为车辆的轮距,h为车辆的质心高度,g为重力加速度,
Figure GDA0004001832220000378
为车辆传感系统识别到的车辆当前所在车道的路面类型对应的路面附着系数;
自适应参数
Figure GDA0004001832220000379
通过公式:
Figure GDA0004001832220000381
计算获得;
C0为车道中心线相对于车辆坐标系原点的横向偏置距离;
Figure GDA0004001832220000382
为死区为[0,η]的死区函数,[0,η]函数为零;Yr为车辆的偏航率;δ为对车道宽度的一半/>
Figure GDA0004001832220000383
所取的横向偏置距离安全系数,δ∈[0.5,0.9];μ为防止车辆的相对偏航率过大所设置的安全偏航率参数,K为自适应速率参数,γmax为自适应上界。
优选地,第一控制单元具体用于:
通过公式:
Figure GDA0004001832220000384
计算目标方向盘转角θsteer,2C2为车道中心线曲率;C1为车辆的航向角;Lw为车辆当前所在车道的车道宽度;Kθδ为车辆的转向系统传动比;
Figure GDA0004001832220000385
为自适应参数;l为车辆的轴距;v为车辆的车速;ε为自适应参数初始值,ε为常数;δ对车道宽度的一半/>
Figure GDA0004001832220000386
所取的横向偏置距离安全系数,δ∈[0.5,0.9];S为滑模变量,通过公式/>
Figure GDA0004001832220000387
计算获得,滑模超平面方程满足:/>
Figure GDA0004001832220000388
θ为车辆偏离车道中心线行驶时,使车辆渐进的回到第一偏置区域Z1内行驶的安全航向角;安全航向角θ通过公式:
Figure GDA0004001832220000389
计算获得;θ∈[θminmax],θmax为保证车辆渐进的趋近平衡点的安全上限,θmin为保证车辆渐进的趋近平衡点的安全下限;aymax为车辆最大安全侧向加速度;车辆最大安全侧向加速度aymax通过公式:
Figure GDA0004001832220000391
计算获得;T为车辆的轮距,h为车辆的质心高度,g为重力加速度,
Figure GDA0004001832220000392
为前视摄像头识别到的车辆当前所在车道的路面类型对应的路面附着系数;
自适应参数
Figure GDA0004001832220000393
通过公式:/>
Figure GDA0004001832220000394
计算获得;
C0为车道中心线相对于车辆坐标系原点的横向偏置距离;
Figure GDA0004001832220000395
为死区为[0,η]的死区函数,[0,η]函数为零;Yr为车辆的偏航率;δ为对车道宽度的一半/>
Figure GDA0004001832220000396
所取的横向偏置距离安全系数,δ∈[0.5,0.9];μ为防止车辆的相对偏航率过大所设置的安全偏航率参数,K为自适应速率参数,γmax为自适应上界;
饱和函数为
Figure GDA0004001832220000397
类饱和函数为/>
Figure GDA0004001832220000398
Figure GDA0004001832220000399
是任意连续函数,满足f(0)=0,f(1)=1,f(-1)=-1;φ为准滑模边界层宽度,φ为常数。
优选地,第一控制单元和第二控制单元具体用于:
通过公式:
Figure GDA00040018322200003910
计算出安全减速度a1
g为重力加速度;
Figure GDA0004001832220000401
为前视摄像头识别到的车辆当前所在车道的路面类型估计出的对应的路面附着系数;a1max为车辆的最大安全减速度,/>
Figure GDA0004001832220000402
a1min为车辆的最小安全减速度,a1min为常数;tD1为从车辆从第一偏置区域Z1进入第二偏置区域Z2的时间点;tNow为当前时间点;Yr为车辆的偏航率;v为车辆的车速;
σ为根据车辆的转向能力和底盘稳定能力所确定的减速度系数,通过公式:
Figure GDA0004001832220000403
计算获得;LTRCrit为安全临界荷载偏移率,LTRCrit为常数;LTRavg为平均荷载偏移率,通过公式:/>
Figure GDA0004001832220000404
计算获得,tnow为当前时间点,tCrit为判定时间,T为车辆的轮距,h为车辆的质心高度;
α为故障估计量化积分参数,通过公式:
Figure GDA0004001832220000405
计算获得,ΔθCrit为安全临界航向角增量,ΔθCrit常数;
S为滑模变量,通过公式:
Figure GDA0004001832220000406
计算;C0表示车道中心线相对于车辆坐标系原点的横向偏置距离;C1为车辆的航向角;2C2为车道中心线曲率;Lw为车辆当前所在车道的车道宽度;δ为对车道宽度的一半/>
Figure GDA0004001832220000407
所取的横向偏置距离安全系数,δ∈[0.5,0.9];θ为车辆偏离车道中心线行驶时,使车辆渐进的回到第一偏置区域Z1内行驶的安全航向角;安全航向角θ通过公式:
Figure GDA0004001832220000411
计算获得;θ∈[θminmax],θmax为保证车辆渐进的趋近平衡点的安全上限,θmin为保证车辆渐进的趋近平衡点的安全下限;aymax为车辆最大安全侧向加速度;车辆最大安全侧向加速度aymax通过公式:
Figure GDA0004001832220000412
计算获得;T为车辆的轮距,h为车辆的质心高度,g为重力加速度,
Figure GDA0004001832220000413
为车辆传感系统识别到的车辆当前所在车道的路面类型对应的路面附着系数。
优选地,第一控制单元具体用于:
通过公式:
Tmax=min(TeRequire,Te-Limit)计算最大安全力矩Tmax
第一力矩TeRequire通过根据车辆的车速v从车速与第一力矩TeRequire的预设关系表查表获得,第二力矩Te-Limit由电动助力转向系统EPS输出;
通过公式:
Figure GDA0004001832220000414
计算得到柔性控制时间tSIL,且所述柔性控制时间tSIL满足:2s≥tSIL≥100ms。
优选地,第一控制单元具体用于:
通过sigmoid函数插值,确定不同插值时刻t对应的插值方向盘转角指令值θSteerInterp;通过公式:
Figure GDA0004001832220000421
计算插值方向盘转角指令值θSteerInterp,插值时间t∈[0,tSIL];θSteerDifference为第一等级冗余控制功能激活的时间点对应的车辆的方向盘转角实际反馈值θSteerFB与目标方向盘转角θsteer之差;
将在柔性控制时间tSIL内的每一插值时刻t对应的插值方向盘转角θSteerInterp输出给电动助力转向系统EPS,使电动助力转向系统EPS在所述最大安全力矩Tmax条件下将车辆的方向盘转角由方向盘转角实际反馈值θSteerFB调节为所述插值方向盘转角指令值θSteerInterp
优选地,所述仲裁模块还包括:
第一判断单元,用于针对自动驾驶系统为L4或L5高等级自动驾驶系统,在按照第一等级冗余控制功能或第二等级冗余控制功能进行车辆控制时,判断是否满足:车辆当前所处的偏置区域稳定为第一偏置区域Z1、车辆的车速是否小于车辆当前所在车道的最低限制车速且车辆当前所在车道的路况为良好路况;
第一输出单元,用于若满足,输出提示驾驶员进行车辆接管的提示信息;
第二判断单元,用于判断驾驶员是否接管车辆;
第一关闭单元,用于若接管,则关闭第一等级冗余控制功能或第二等级冗余控制功能,并停止向电动助力转向系统EPS和车身电子稳定系统ESP输出控制指令;
第二输出单元,用于针对自动驾驶系统为L2.5或L3中等级自动驾驶系统,在按照第一等级冗余控制功能或第二等级冗余控制功能进行车辆控制时,输出提示驾驶员进行车辆接管的提示信息;
第三判断单元,用于判断驾驶员是否接管车辆;
第二关闭单元,用于若接管,则关闭第一等级冗余控制功能或第二等级冗余控制功能,并停止向电动助力转向系统EPS和车身电子稳定系统ESP输出控制指令。
参照图1,本发明实施例还提供了一种自动驾驶汽车,包括上述的自动驾驶系统。
本发明实施例还提供了一种控制器,包括存储器和处理器,所述存储器中存储有至少一条指令、至少一段程序、代码集或指令集,所述至少一条指令、至少一段程序、代码集或指令集由所述处理器加载并执行,以实现如上述的冗余控制方法。
本发明实施例还提供了一种计算机可读存储介质,所述存储介质上存储有至少一条指令、至少一段程序、代码集或指令集,所述至少一条指令、至少一段程序、代码集或指令集由所述处理器加载并执行,以实现如上述的冗余控制方法。

Claims (18)

1.一种自动驾驶系统的冗余控制方法,所述自动驾驶系统包括主控制子系统和冗余控制子系统,所述方法应用于冗余控制子系统,其特征在于,所述方法包括:
在自动驾驶系统的自动驾驶功能激活时,确定车辆在当前所在车道内具体所处的偏置区域;偏置区域是指车道中心线左侧方向和车道中心线右侧方向各自与车道中心线之间的距离位于预定距离范围内时所组成的区域;
根据车辆当前所处的偏置区域和冗余控制子系统在上一控制周期输出的控制结果,确定主控制子系统当前的故障等级;
根据主控制子系统当前的故障等级,激活对应等级的冗余控制功能,并按照对应等级的冗余控制功能进行车辆控制。
2.根据权利要求1所述的方法,其特征在于,确定车辆在当前所在车道内具体所处的偏置区域的具体步骤为:
根据前视摄像头和雷达输出到的多源信息,拟合使车辆保持车道中心线行驶的安全路径方程C0(x),并确定车辆当前所在车道的车道宽度Lw;
根据所规划的安全路径方程C0(x)中的系数C0和车道宽度Lw,确定车辆在当前所在车道内具体所处的偏置区域。
3.根据权利要求2所述的方法,其特征在于,摄像头和雷达输出到的多源信息包括:相对于车辆坐标系原点规划的车道线曲线方程L0(x)、路沿曲线方程K0(x)、护栏曲线方程G0(x)和/或车流曲线方程T0(x);根据前视摄像头和雷达输出到的多源信息,拟合使车辆保持车道中心线行驶的安全路径方程C0(x),并确定车辆当前所在车道的车道宽度Lw的具体步骤为:
根据车辆的车速向前等间隔时长的选取N个纵向位置点,通过公式:xi=v×i×tLA(i=0,...,N-1)计算第i个纵向位置点与车辆坐标系原点之间的纵向距离xi,v为车辆的车速,tLA为设定的时间间隔时长;
通过公式:
E0Left(xi)=-min[|L0Left(xi)|,|K0Left(xi)|,|G0Left(xi)|,|T0Left(xi)|](i=0,...,N-1)确定与第i个纵向位置点对应的左侧横向位置点E0Left(xi);
通过公式:
E0Right(xi)=min[|L0Right(xi)|,|K0Right(xi)|,|G0Right(xi)|,|T0Right(xi)|](i=0,...,N-1)确定与第i个纵向位置点对应的右侧横向位置点E0Left(xi);
使用最小二乘拟合法将取得的车辆所在车道的左侧边缘的N个位置点(xi,E0Left(xi))拟合为车辆所在车道的左侧边缘曲线E0Left(x);
使用最小二乘拟合法将取得的车辆所在车道的右侧边缘的N个位置点(xi,E0Right(xi))拟合为车辆所在车道的右侧边缘曲线E0Right(x);
对左侧边缘曲线和右侧边缘曲线的各系数进行均值处理,通过公式:
Figure FDA0004001832210000021
获得使车辆沿着车道中心线行驶的安全路径方程C0(x);C0为车辆所在车道的车道中心线相对于车辆坐标系原点的横向偏置距离,C1为车辆所在车道的车道中心线相对车辆坐标系原点的夹角,C2为车道中心线曲率的1/2,C3为车道中心线曲率变化率的1/6;
通过公式:
Lw=min(Lw(x0),...,Lw(xN-1),Lwmin(xLB));
计算车辆当前所在车道的车道宽度Lw;Lwmin(xLB)表示在前tLB时间内的历史最小车道宽度,Lw(xi)=|E0Right(xi)-E0Left(xi)|(i=0,...,N-1)。
4.根据权利要求2或3所述的方法,其特征在于,根据所规划的安全路径方程C0(x)中的系数C0和车道宽度Lw,确定车辆在当前所在车道内具体所处的偏置区域的具体步骤为:
通过公式:
Figure FDA0004001832210000022
计算以车道宽度Lw和横向偏置距离C0为基础的故障函数F(Δ);Lw为车辆当前所在车道的车道宽度,系数C0为车道中心线相对于车辆坐标系原点的横向偏置距离;
若故障函数F(Δ)小于第一设定比值D1;则确定车辆在当前所在车道内具体所处的偏置区域为第一偏置区域Z1;
若故障函数F(Δ)大于或等于第一设定比值D1且小于第二设定比值D2;则确定车辆在当前所在车道内具体所处的偏置区域为第二偏置区域Z2;
若故障函数F(Δ)大于或等于第二设定比值D2;则确定车辆在当前所在车道内具体所处的偏置区域为第三偏置区域Z3。
5.根据权利要求4所述的方法,其特征在于,根据车辆当前所处的偏置区域和冗余控制子系统在上一控制周期输出的控制结果,确定主控制子系统当前的故障等级的具体步骤为:
若冗余控制子系统在上一控制周期输出的控制结果为激活第一等级冗余控制功能,或者,若冗余控制子系统在上一控制周期输出的控制结果为激活第二等级冗余控制功能或未激活冗余控制功能且车辆在当前所在车道内具体所处的偏置区域为第三偏置区域Z3,则确定主控制子系统当前的故障等级为高故障等级;
若冗余控制子系统在上一控制周期输出的控制结果为激活第二等级冗余控制功能或未激活冗余控制功能且车辆在当前所在车道内具体所处的偏置区域为第二偏置区域Z2,则进一步比对车辆姿态的第一判别函数G(S)与第二设定比值D2的大小;在车辆姿态的第一判别函数G(S)大于或等于第二设定比值D2时,确定主控制子系统当前的故障等级为高故障等级;在车辆姿态的第一判别函数G(S)小于第二设定比值D2时,确定主控制子系统当前的故障等级为中故障等级;
若冗余控制子系统在上一控制周期输出的控制结果为激活第二等级冗余控制功能或未激活冗余控制功能且车辆在当前所在车道内具体所处的偏置区域为第一偏置区域Z1,则确定主控制子系统当前的故障等级为低故障等级。
6.根据权利要求5所述的方法,其特征在于,车辆姿态的第一判别函数G(S)通过公式:
G(S)=λ|S|
计算获得;λ为识别信度,λ为常数或者基于车辆的车速v、车辆的偏航率Yr、方向盘转角θSteer、车道中心线曲率2C2以及前视摄像头识别到的车辆当前所在车道的路面类型所对应的路面附着系数
Figure FDA0004001832210000031
进行RBF神经网络拟合得到;当λ为常数时λ=0.7;当λ通过RBF神经网络拟合得到时,λ的取值范围为:λ∈[0.35,1];
S为滑模变量,通过公式:
Figure FDA0004001832210000032
计算获得;滑模超平面方程满足:
Figure FDA0004001832210000033
Lw为车辆当前所在车道的车道宽度;C0为车道中心线相对于车辆坐标系原点的横向偏置距离;C1为车辆的航向角;δ为对车道宽度的一半
Figure FDA0004001832210000034
所取的横向偏置距离安全系数,δ∈[0.5,0.9];θ为车辆偏离车道中心线行驶时,使车辆渐进的回到第一偏置区域Z1内行驶的安全航向角;安全航向角θ通过公式:
Figure FDA0004001832210000041
计算获得;θ∈[θminmax],θmax为保证车辆渐进的趋近平衡点的安全上限,θmin为保证车辆渐进的趋近平衡点的安全下限;aymax为车辆最大安全侧向加速度;车辆最大安全侧向加速度aymax通过公式:
Figure FDA0004001832210000042
计算获得;T为车辆的轮距,h为车辆的质心高度,g为重力加速度,
Figure FDA0004001832210000043
为前视摄像头识别到的车辆当前所在车道的路面类型对应的路面附着系数。
7.根据权利要求6所述的方法,其特征在于,根据主控制子系统当前的故障等级,激活对应等级的冗余控制功能的具体步骤为:
若主控制子系统当前的故障等级为高故障等级,则激活第一等级冗余控制功能;
若主控制子系统当前的故障等级为中故障等级,则激活第二等级冗余控制功能;
若主控制子系统当前的故障等级为低故障等级,则不激活第一等级冗余控制功能或第二等级冗余控制功能。
8.根据权利要求7所述的方法,其特征在于,按照对应等级的冗余控制功能进行车辆控制的具体步骤为:
若激活第一等级冗余控制功能,则按照计算出的综合减速度a控制车身电子稳定系统ESP减速至第一目标减速度,并按照计算出的目标方向盘转角θsteer和最大安全力矩Tmax控制电动助力转向系统EPS在柔性控制时间tSIL内进行方向盘转角控制;
若激活第二等级冗余控制功能,则根据车辆姿态的第二判别函数H(S)计算出的安全减速度a1和主控制子系统输出的减速度a2控制车身电子稳定系统ESP减速至第二目标减速度;
综合减速度a为根据车辆姿态的第二判别函数H(S)计算出的安全减速度a1和主控制子系统输出的减速度a2中的较大值;第一目标速度与当前行驶环境工况相关;第二目标速度与主控制子系统输出的目标速度相关,且第二目标速度小于主控制子系统输出的目标速度。
9.根据权利要求8所述的方法,其特征在于,
通过公式:
Figure FDA0004001832210000051
计算目标方向盘转角θsteer;2C2为车道中心线曲率;C1为车辆的航向角;Lw为车辆当前所在车道的车道宽度;Kθδ为车辆的转向系统传动比;
Figure FDA0004001832210000052
为自适应参数;l为车辆的轴距;v为车辆的车速;ε为自适应参数初始值,ε为常数;δ为对车道宽度的一半/>
Figure FDA0004001832210000053
所取的横向偏置距离安全系数,δ∈[0.5,0.9];S为滑模变量,通过公式/>
Figure FDA0004001832210000054
计算获得,滑模超平面方程满足:/>
Figure FDA0004001832210000055
θ为车辆偏离车道中心线行驶时,使车辆渐进的回到第一偏置区域Z1内行驶的安全航向角;安全航向角θ通过公式:
Figure FDA0004001832210000056
计算获得;θ∈[θminmax],θmax为保证车辆渐进的趋近平衡点的安全上限,θmin为保证车辆渐进的趋近平衡点的安全下限;aymax为车辆最大安全侧向加速度;车辆最大安全侧向加速度aymax通过公式:
Figure FDA0004001832210000057
计算获得;T为车辆的轮距,h为车辆的质心高度,g为重力加速度,
Figure FDA0004001832210000058
为车辆传感系统识别到的车辆当前所在车道的路面类型对应的路面附着系数;
自适应参数
Figure FDA0004001832210000059
通过公式:/>
Figure FDA0004001832210000061
计算获得;C0为车道中心线相对于车辆坐标系原点的横向偏置距离;/>
Figure FDA0004001832210000062
为死区为[0,η]的死区函数,[0,η]函数为零;Yr为车辆的偏航率;δ为对车道宽度的一半
Figure FDA0004001832210000063
所取的横向偏置距离安全系数,δ∈[0.5,0.9];μ为防止车辆的相对偏航率过大所设置的安全偏航率参数,K为自适应速率参数,γmax为自适应上界。
10.根据权利要求8所述的方法,其特征在于,通过公式:
Figure FDA0004001832210000064
计算目标方向盘转角θsteer,2C2为车道中心线曲率;C1为车辆的航向角;Lw为车辆当前所在车道的车道宽度;Kθδ为车辆的转向系统传动比;
Figure FDA0004001832210000065
为自适应参数;l为车辆的轴距;v为车辆的车速;ε为自适应参数初始值,ε为常数;δ对车道宽度的一半/>
Figure FDA0004001832210000066
所取的横向偏置距离安全系数,δ∈[0.5,0.9];S为滑模变量,通过公式/>
Figure FDA0004001832210000067
计算获得,滑模超平面方程满足:/>
Figure FDA0004001832210000068
θ为车辆偏离车道中心线行驶时,使车辆渐进的回到第一偏置区域Z1内行驶的安全航向角;安全航向角θ通过公式:
Figure FDA0004001832210000069
计算获得;θ∈[θminmax],θmax为保证车辆渐进的趋近平衡点的安全上限,θmin为保证车辆渐进的趋近平衡点的安全下限;aymax为车辆最大安全侧向加速度;车辆最大安全侧向加速度aymax通过公式:
Figure FDA0004001832210000071
计算获得;T为车辆的轮距,h为车辆的质心高度,g为重力加速度,
Figure FDA0004001832210000072
为前视摄像头识别到的车辆当前所在车道的路面类型对应的路面附着系数;/>
自适应参数
Figure FDA0004001832210000073
通过公式:
Figure FDA0004001832210000074
C0为车道中心线相对于车辆坐标系原点的横向偏置距离;
Figure FDA0004001832210000075
为死区为[0,η]的死区函数,[0,η]函数为零;Yr为车辆的偏航率;δ为对车道宽度的一半
Figure FDA0004001832210000076
所取的横向偏置距离安全系数,δ∈[0.5,0.9];μ为防止车辆的相对偏航率过大所设置的安全偏航率参数,K为自适应速率参数,γmax为自适应上界;
饱和函数为
Figure FDA0004001832210000077
类饱和函数为/>
Figure FDA0004001832210000078
是任意连续函数,满足f(0)=0,f(1)=1,f(-1)=-1;φ为准滑模边界层宽度,φ为常数。
11.根据权利要求10所述的方法,其特征在于,根据车辆姿态的第二判别函数H(S)通过公式:
Figure FDA0004001832210000079
计算出安全减速度a1
g为重力加速度;
Figure FDA0004001832210000081
为前视摄像头识别到的车辆当前所在车道的路面类型估计出的对应的路面附着系数;a1max为车辆的最大安全减速度,/>
Figure FDA0004001832210000082
a1min为车辆的最小安全减速度,a1min为常数;tD1为从车辆从第一偏置区域Z1进入第二偏置区域Z2的时间点;tNow为当前时间点;Yr为车辆的偏航率;v为车辆的车速;
σ为根据车辆的转向能力和底盘稳定能力所确定的减速度系数,通过公式:
Figure FDA0004001832210000083
计算获得;LTRCrit为安全临界荷载偏移率,LTRCrit为常数;LTRavg为平均荷载偏移率,通过公式:/>
Figure FDA0004001832210000084
计算获得,tnow为当前时间点,tCrit为判定时间,T为车辆的轮距,h为车辆的质心高度;
α为故障估计量化积分参数,通过公式:
Figure FDA0004001832210000085
计算获得,ΔθCrit为安全临界航向角增量,ΔθCrit常数;/>
S为滑模变量,通过公式:
Figure FDA0004001832210000086
计算;C0表示车道中心线相对于车辆坐标系原点的横向偏置距离;C1为车辆的航向角;2C2为车道中心线曲率;Lw为车辆当前所在车道的车道宽度;δ为对车道宽度的一半/>
Figure FDA0004001832210000087
所取的横向偏置距离安全系数,δ∈[0.5,0.9];θ为车辆偏离车道中心线行驶时,使车辆渐进的回到第一偏置区域Z1内行驶的安全航向角;安全航向角θ通过公式:
Figure FDA0004001832210000091
计算获得;θ∈[θminmax],θmax为保证车辆渐进的趋近平衡点的安全上限,θmin为保证车辆渐进的趋近平衡点的安全下限;aymax为车辆最大安全侧向加速度;车辆最大安全侧向加速度aymax通过公式:
Figure FDA0004001832210000092
计算获得;T为车辆的轮距,h为车辆的质心高度,g为重力加速度,
Figure FDA0004001832210000093
为车辆传感系统识别到的车辆当前所在车道的路面类型对应的路面附着系数。
12.根据权利要求11所述的方法,其特征在于,通过公式:
Tmax=min(TeRequire,Te-Limit)计算最大安全力矩Tmax
第一力矩TeRequire通过根据车辆的车速v从车速与第一力矩TeRequire的预设关系表查表获得,第二力矩Te-Limit由电动助力转向系统EPS输出;通过公式:
Figure FDA0004001832210000094
计算得到柔性控制时间tSIL,且所述柔性控制时间tSIL满足:2s≥tSIL≥100ms。
13.根据权利要求10所述的方法,其特征在于,按照计算出的方向盘转角θsteer和最大安全力矩Tmax控制电动助力转向系统EPS在柔性控制时间tSIL内进行方向盘转角控制的具体步骤为:
通过sigmoid函数插值,确定不同插值时刻t对应的插值方向盘转角指令值θSteerInterp;通过公式:
Figure FDA0004001832210000101
计算插值方向盘转角指令值θSteerInterp,插值时间t∈[0,tSIL];θSteerDifference为第一等级冗余控制功能激活的时间点对应的车辆的方向盘转角实际反馈值θSteerFB与目标方向盘转角θsteer之差;
将在柔性控制时间tSIL内的每一插值时刻t对应的插值方向盘转角θSteerInterp输出给电动助力转向系统EPS,使电动助力转向系统EPS在所述最大安全力矩Tmax条件下将车辆的方向盘转角由方向盘转角实际反馈值θSteerFB调节为所述插值方向盘转角指令值θSteerInterp
14.根据权利要求1所述的方法,其特征在于,所述方法还包括:
针对自动驾驶系统为L4或L5高等级自动驾驶系统,在按照第一等级冗余控制功能或第二等级冗余控制功能进行车辆控制时,判断是否满足:车辆当前所处的偏置区域稳定为第一偏置区域Z1、车辆的车速是否小于车辆当前所在车道的最低限制车速且车辆当前所在车道的路况为良好路况;
若满足,输出提示安全员进行车辆接管的提示信息;
判断安全员是否接管车辆;
若接管,则关闭第一等级冗余控制功能或第二等级冗余控制功能,并停止向电动助力转向系统EPS和车身电子稳定系统ESP输出控制指令;
针对自动驾驶系统为L2.5或L3中等级自动驾驶系统,在按照第一等级冗余控制功能或第二等级冗余控制功能进行车辆控制时,输出提示驾驶员进行车辆接管的提示信息;
判断驾驶员是否接管车辆;
若接管,则关闭第一等级冗余控制功能或第二等级冗余控制功能,并停止向电动助力转向系统EPS和车身电子稳定系统ESP输出控制指令。
15.一种自动驾驶系统,所述自动驾驶系统包括主控制子系统和冗余控制子系统,其特征在于,所述冗余控制子系统包括:
第一确定模块,用于在自动驾驶系统的自动驾驶功能激活时,确定车辆在当前所在车道内具体所处偏置区域;偏置区域是指车道中心线左侧方向和车道中心线右侧方向各自与车道中心线之间的距离位于预定距离范围内时所组成的区域;
第二确定模块,用于根据车辆当前所处的偏置区域和冗余控制子系统在上一控制周期输出的控制结果,确定主控制子系统当前的故障等级;
仲裁模块,用于根据主控制子系统当前的故障等级,激活对应等级的冗余控制功能,并按照对应等级的冗余控制功能进行车辆控制。
16.一种自动驾驶汽车,其特征在于,包括权利要求15所述的自动驾驶系统。
17.一种控制器,其特征在于,包括存储器和处理器,所述存储器中存储有至少一条指令、至少一段程序、代码集或指令集,所述至少一条指令、至少一段程序、代码集或指令集由所述处理器加载并执行,以实现如权利要求1至14任一所述的冗余控制方法。
18.一种计算机可读存储介质,其特征在于,所述存储介质上存储有至少一条指令、至少一段程序、代码集或指令集,所述至少一条指令、至少一段程序、代码集或指令集由处理器加载并执行,以实现如权利要求1至14任一所述的冗余控制方法。
CN202011317326.9A 2020-11-23 2020-11-23 自动驾驶系统的冗余控制方法、自动驾驶系统、汽车、控制器及计算机可读存储介质 Active CN112373477B (zh)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN202011317326.9A CN112373477B (zh) 2020-11-23 2020-11-23 自动驾驶系统的冗余控制方法、自动驾驶系统、汽车、控制器及计算机可读存储介质
RU2021104308A RU2759453C1 (ru) 2020-11-23 2020-12-03 Способ управления избыточностью для системы автоматического вождения, автоматическая система вождения, транспортное средство, контроллер и машиночитаемый носитель информации
PCT/CN2020/133607 WO2022104925A1 (zh) 2020-11-23 2020-12-03 自动驾驶系统的冗余控制方法、自动驾驶系统、汽车、控制器及计算机可读存储介质

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011317326.9A CN112373477B (zh) 2020-11-23 2020-11-23 自动驾驶系统的冗余控制方法、自动驾驶系统、汽车、控制器及计算机可读存储介质

Publications (2)

Publication Number Publication Date
CN112373477A CN112373477A (zh) 2021-02-19
CN112373477B true CN112373477B (zh) 2023-03-28

Family

ID=74587301

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011317326.9A Active CN112373477B (zh) 2020-11-23 2020-11-23 自动驾驶系统的冗余控制方法、自动驾驶系统、汽车、控制器及计算机可读存储介质

Country Status (2)

Country Link
CN (1) CN112373477B (zh)
WO (1) WO2022104925A1 (zh)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113044063B (zh) * 2021-03-31 2022-09-06 重庆长安汽车股份有限公司 用于高级自动驾驶的功能冗余控制系统
CN113511221B (zh) * 2021-05-20 2022-10-11 重庆长安汽车股份有限公司 一种监控横向控制能力的方法、系统、车辆及存储介质
CN113341940B (zh) * 2021-06-21 2022-09-16 安徽江淮汽车集团股份有限公司 并联式自动驾驶系统
CN114348016B (zh) * 2021-11-23 2023-09-19 重庆长安汽车股份有限公司 一种自动驾驶纵向减速度合理性校验方法及系统
CN114162068B (zh) * 2021-12-31 2023-12-15 阿维塔科技(重庆)有限公司 车辆智能驾驶功能的管理方法、装置以及车辆
CN114885103A (zh) * 2022-06-14 2022-08-09 智道网联科技(北京)有限公司 多路车载摄像头视频数据的处理方法及装置
CN115158445B (zh) * 2022-07-11 2023-06-23 岚图汽车科技有限公司 一种车辆转向冗余控制方法、系统、电子设备及存储介质
CN116150771B (zh) * 2023-04-19 2023-07-04 北京集度科技有限公司 一种冗余处理的装置、方法、系统及智能车辆
CN116229723B (zh) * 2023-05-06 2023-08-11 山东纵云信息技术有限公司 一种智能化交通信息数据分析管理系统
CN116339120A (zh) * 2023-05-30 2023-06-27 天津德科智控股份有限公司 Eps双冗余系统及故障处理方法
CN116383330B (zh) * 2023-06-06 2023-08-11 中航信移动科技有限公司 一种轨迹拟合方法、存储介质及电子设备

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108422998A (zh) * 2018-03-14 2018-08-21 合肥工业大学 一种车道偏离辅助控制系统及方法
JP2018177181A (ja) * 2017-04-12 2018-11-15 トヨタ自動車株式会社 車線変更支援装置

Family Cites Families (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3760827B2 (ja) * 2001-09-28 2006-03-29 日産自動車株式会社 車線逸脱防止装置
JP4389567B2 (ja) * 2003-12-03 2009-12-24 日産自動車株式会社 車線逸脱防止装置
JP5672966B2 (ja) * 2010-10-29 2015-02-18 株式会社デンソー 車両運動制御システム
WO2012068331A1 (en) * 2010-11-19 2012-05-24 Magna Electronics Inc. Lane keeping system and lane centering system
DE102013100577B4 (de) * 2012-01-25 2024-02-29 Denso Corporation Spurhalte-Steuersystem
JP5543501B2 (ja) * 2012-01-27 2014-07-09 株式会社日本自動車部品総合研究所 車両制御装置
CN103832433B (zh) * 2012-11-21 2016-08-10 中国科学院沈阳计算技术研究所有限公司 车道偏离及前车防碰撞报警系统及其实现方法
JP5967224B2 (ja) * 2013-01-10 2016-08-10 日産自動車株式会社 車線内走行支援装置
AT514754B1 (de) * 2013-09-05 2018-06-15 Avl List Gmbh Verfahren und Vorrichtung zur Optimierung von Fahrassistenzsystemen
CN104442814B (zh) * 2014-10-31 2016-07-06 重庆长安汽车股份有限公司 基于偏差预测算法的车道保持及自动对中系统和方法
DE102015201569A1 (de) * 2015-01-29 2016-08-04 Continental Teves Ag & Co. Ohg Fahrzeugsteuervorrichtung und verfahren
CN105711588B (zh) * 2016-01-20 2018-05-11 奇瑞汽车股份有限公司 一种车道保持辅助系统和车道保持辅助方法
CN110191832A (zh) * 2017-01-24 2019-08-30 本田技研工业株式会社 车辆控制装置、车辆控制方法及车辆控制程序
JP6801591B2 (ja) * 2017-06-06 2020-12-16 トヨタ自動車株式会社 操舵支援装置
KR102371610B1 (ko) * 2017-08-10 2022-03-07 현대자동차주식회사 차로 유지 제어 방법 및 시스템
WO2019066949A1 (en) * 2017-09-29 2019-04-04 Intel Corporation RANDOMIZATION OF TRACK MOVEMENT OF AUTOMATED VEHICLES
CN108647572B (zh) * 2018-04-04 2022-03-29 华南理工大学 一种基于霍夫变换的车道偏离预警方法
JP7205768B2 (ja) * 2019-03-08 2023-01-17 スズキ株式会社 車両の走行制御装置
CN109760677B (zh) * 2019-03-13 2020-09-11 广州小鹏汽车科技有限公司 一种车道保持辅助方法及系统
CN110979305B (zh) * 2019-03-18 2021-06-22 毫末智行科技有限公司 车辆异常换道控制方法、装置及系统
CN109917779A (zh) * 2019-03-26 2019-06-21 中国第一汽车股份有限公司 面向l3自动驾驶的冗余控制系统
JP7274327B2 (ja) * 2019-03-28 2023-05-16 株式会社Subaru 自動運転支援装置
CN111487971B (zh) * 2020-04-23 2023-04-28 重庆长安汽车股份有限公司 一种车辆自动驾驶横向控制方法及系统

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2018177181A (ja) * 2017-04-12 2018-11-15 トヨタ自動車株式会社 車線変更支援装置
CN108422998A (zh) * 2018-03-14 2018-08-21 合肥工业大学 一种车道偏离辅助控制系统及方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
郭烈 ; 葛平淑 ; 夏文旭 ; 秦增科 ; .基于人机共驾的车道保持辅助控制系统研究.中国公路学报.2019,(第12期),50-61. *

Also Published As

Publication number Publication date
CN112373477A (zh) 2021-02-19
WO2022104925A1 (zh) 2022-05-27

Similar Documents

Publication Publication Date Title
CN112373477B (zh) 自动驾驶系统的冗余控制方法、自动驾驶系统、汽车、控制器及计算机可读存储介质
US11059462B2 (en) Method and device for controlling a movement of a vehicle, and vehicle movement control system
US9180862B2 (en) Vehicle dynamic control apparatus and vehicle dynamic control system using the same
US8170739B2 (en) Path generation algorithm for automated lane centering and lane changing control system
US8190330B2 (en) Model based predictive control for automated lane centering/changing control systems
US10814906B2 (en) Method for operating a power steering system
US8296033B2 (en) Curve-related accident mitigation
EP3674156B1 (en) Method and device for generating target path for autonomous vehicle
EP3674154B1 (en) Method and apparatus for controlling travel of driving-assisted vehicle
CN107640151A (zh) 用于确定车轮对准改变的装置和方法
JP2009061878A (ja) 走行制御装置
CN115023380A (zh) 非对称性的故障安全的系统结构
KR20200089275A (ko) 적어도 부분적으로 자동으로 주행하는 자동차를 위한 운전자 보조 시스템, 자동차 및 자동차 동역학을 제어하기 위한 방법
US10988166B2 (en) System and method for adapting steering torque based on a motor vehicle present driving situation
Ali et al. Predictive prevention of loss of vehicle control for roadway departure avoidance
US11897452B2 (en) Automotive electronic lateral dynamics control system for a self-driving motor vehicle
US20220324466A1 (en) Device for predictively controlling the movement of a motor vehicle
LeBlanc et al. A warning and intervention system to prevent road-departure accidents
RU2759453C1 (ru) Способ управления избыточностью для системы автоматического вождения, автоматическая система вождения, транспортное средство, контроллер и машиночитаемый носитель информации
US20230391399A1 (en) Method for controlling the lateral position of a vehicle on a traffic lane
Lee Automatic speed control system for vehicle approaching and driving on a curve
Liang et al. Realization of Emergency Lane Keeping System by Adaptive Control based on the Finite State Machine
US20230294717A1 (en) Method for Determining a Trajectory for Controlling a Vehicle
US11794777B1 (en) Systems and methods for estimating heading and yaw rate for automated driving
WO2023037758A1 (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