CN107878453A - 一种躲避动障碍物的汽车紧急避撞一体式控制方法 - Google Patents
一种躲避动障碍物的汽车紧急避撞一体式控制方法 Download PDFInfo
- Publication number
- CN107878453A CN107878453A CN201711081210.8A CN201711081210A CN107878453A CN 107878453 A CN107878453 A CN 107878453A CN 201711081210 A CN201711081210 A CN 201711081210A CN 107878453 A CN107878453 A CN 107878453A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- mtd
- mtr
- automobile
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 230000004888 barrier function Effects 0.000 title claims abstract description 56
- 238000000034 method Methods 0.000 title claims abstract description 37
- 238000005457 optimization Methods 0.000 claims abstract description 9
- 238000013461 design Methods 0.000 claims description 17
- 230000033001 locomotion Effects 0.000 claims description 16
- 230000035772 mutation Effects 0.000 claims description 13
- 230000001133 acceleration Effects 0.000 claims description 9
- 239000005364 simax Substances 0.000 claims description 7
- 238000004458 analytical method Methods 0.000 claims description 6
- 238000002474 experimental method Methods 0.000 claims description 5
- 238000012938 design process Methods 0.000 claims description 4
- 238000013459 approach Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 claims description 3
- 238000002513 implantation Methods 0.000 claims description 3
- 230000009286 beneficial effect Effects 0.000 claims description 2
- 238000005259 measurement Methods 0.000 claims description 2
- XXXSILNSXNPGKG-ZHACJKMWSA-N Crotoxyphos Chemical compound COP(=O)(OC)O\C(C)=C\C(=O)OC(C)C1=CC=CC=C1 XXXSILNSXNPGKG-ZHACJKMWSA-N 0.000 description 3
- 208000012661 Dyskinesia Diseases 0.000 description 3
- 238000012360 testing method Methods 0.000 description 3
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 230000001953 sensory effect Effects 0.000 description 2
- 238000000926 separation method Methods 0.000 description 2
- 102100034112 Alkyldihydroxyacetonephosphate synthase, peroxisomal Human genes 0.000 description 1
- 101000799143 Homo sapiens Alkyldihydroxyacetonephosphate synthase, peroxisomal Proteins 0.000 description 1
- 238000000848 angular dependent Auger electron spectroscopy Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005096 rolling process Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
Classifications
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60W—CONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
- B60W30/00—Purposes of road vehicle drive control systems not related to the control of a particular sub-unit, e.g. of systems using conjoint control of vehicle sub-units
- B60W30/08—Active safety systems predicting or avoiding probable or impending collision or attempting to minimise its consequences
- B60W30/09—Taking automatic action to avoid collision, e.g. braking and steering
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60W—CONJOINT 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
- B60W10/00—Conjoint control of vehicle sub-units of different type or different function
- B60W10/18—Conjoint control of vehicle sub-units of different type or different function including control of braking systems
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60W—CONJOINT 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
- B60W10/00—Conjoint control of vehicle sub-units of different type or different function
- B60W10/20—Conjoint control of vehicle sub-units of different type or different function including control of steering systems
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60W—CONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
- B60W30/00—Purposes of road vehicle drive control systems not related to the control of a particular sub-unit, e.g. of systems using conjoint control of vehicle sub-units
- B60W30/08—Active safety systems predicting or avoiding probable or impending collision or attempting to minimise its consequences
- B60W30/095—Predicting travel path or likelihood of collision
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60W—CONJOINT CONTROL OF VEHICLE SUB-UNITS OF DIFFERENT TYPE OR DIFFERENT FUNCTION; CONTROL SYSTEMS SPECIALLY ADAPTED FOR HYBRID VEHICLES; ROAD VEHICLE DRIVE CONTROL SYSTEMS FOR PURPOSES NOT RELATED TO THE CONTROL OF A PARTICULAR SUB-UNIT
- B60W30/00—Purposes of road vehicle drive control systems not related to the control of a particular sub-unit, e.g. of systems using conjoint control of vehicle sub-units
- B60W30/08—Active safety systems predicting or avoiding probable or impending collision or attempting to minimise its consequences
- B60W30/095—Predicting travel path or likelihood of collision
- B60W30/0953—Predicting travel path or likelihood of collision the prediction being responsive to vehicle dynamic parameters
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60W—CONJOINT 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
- B60W2520/00—Input parameters relating to overall vehicle dynamics
- B60W2520/10—Longitudinal speed
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60W—CONJOINT 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
- B60W2520/00—Input parameters relating to overall vehicle dynamics
- B60W2520/12—Lateral speed
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60W—CONJOINT 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
- B60W2554/00—Input parameters relating to objects
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60W—CONJOINT 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
- B60W2710/00—Output or target parameters relating to a particular sub-units
- B60W2710/18—Braking system
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B60—VEHICLES IN GENERAL
- B60W—CONJOINT 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
- B60W2710/00—Output or target parameters relating to a particular sub-units
- B60W2710/20—Steering systems
Landscapes
- Engineering & Computer Science (AREA)
- Transportation (AREA)
- Mechanical Engineering (AREA)
- Automation & Control Theory (AREA)
- Chemical & Material Sciences (AREA)
- Combustion & Propulsion (AREA)
- Steering Control In Accordance With Driving Conditions (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
本发明涉及一种躲避动障碍物的汽车紧急避撞一体式控制方法,其特征在于,该方法是利用路径动态规划与实时跟踪控制模块,根据实时采集的障碍物信息、目标点坐标、汽车行驶状态信息,实时优化得出汽车的前轮转角和四个车轮滑移率,进而控制汽车实现避撞;其中,障碍物信息包括由雷达传感器测量获得的障碍物外形轮廓的离散点坐标,汽车行驶状态信息包括由车速传感器测量获得的汽车纵向速度和侧向速度以及由陀螺仪测量获得的横摆角速度;在控制避撞过程中,通过电动助力转向(Electric Power Steering,EPS)力矩补偿模块根据车速、前轮附加转角,确定力矩补偿控制增益,将转向盘突变力矩控制在理想范围内,实现人机和谐的汽车紧急避撞。
Description
技术领域
本发明涉及汽车先进驾驶辅助技术领域,具体涉及一种躲避动障碍物的汽车紧急避撞一体式控制方法。
背景技术
汽车能给人们带来方便与快捷,其行驶安全性已成为全球性的社会问题。为了进一步提高道路交通安全性,帮助驾驶员减少错误操作,近年来以先进驾驶辅助系统(Advanced Driver Assistance Systems,ADAS)为代表的智能汽车安全技术逐渐得到重视和发展。汽车紧急避撞系统通过执行器的主动干预,辅助驾驶员调整汽车的运动轨迹,实现避撞。它能够在危急时刻拯救驾乘者的生命,有着良好的市场前景。
汽车紧急避撞控制方面已有很多研究成果,能够较好的解决避撞控制问题,但这些研究成果主要针对静止障碍物。在考虑运动障碍物的汽车紧急避撞控制方面,文献[Ackermann C,Isermann R,Min S,etal.Collision avoidance with automatic brakingand swerving[J].IFAC Proceedings Volumes,2014,47(3):10694-10699.]考虑障碍物纵向运动情况,检测汽车与运动障碍物的速度差,决策出避撞的转向时机,即是否可以进行转向避撞,但在避撞过程中没有考虑运动障碍物位置的动态变化,并且没有考虑障碍物侧向运动情况。公开号为CN105539586A的中国专利公开了一种用于自主驾驶的汽车躲避移动障碍物的统一的运动规划方法,该方法考虑障碍物的纵向及侧向运动情况,但只用来决策出避撞的转向时机及避撞路径,同样没有在避撞过程中考虑运动障碍物位置的动态变化。
汽车紧急避撞控制离不开转向系统的主动干预。欧洲现有法规规定转向盘与转向车轮之间必须有机械连接,所以主动前轮转向系统(Active Front Steering,AFS)作为今后线控转向系统(Steering-by-wire,SBW)的过渡产品应运而生。文献[Sumio Sugita,Masayoshi Tomizuka.Cancellation ofUnnatural Reaction Torque in Variable-Gear-Ratio[J].Journal of Dynamic Systems Measurement&Control,2012,134(2):021019.]和文献[AtsushiOshima,XuChen,Sumio Sugita,Masayoshi Tomizuka.Control designfor cancellation of unnatural reaction torque and vibrations in variable-gear-ratio steering system[C].ASME 2013Dynamic Systems and ControlConference.American Society of Mechanical Engineers,2013-3797,V001T11A003:10pages.]论证了AFS在改变系统位移传递特性的同时,也会影响转向系统的力传递特性,引起转向盘力矩的突变。过大的转向盘突变力矩会加剧驾驶员的紧张心理,容易使驾驶员产生误操作,不利于驾驶安全。适当的转向盘突变力矩却有利于驾驶员感知汽车的姿态变化,并起到警示作用。但驾驶员对转向盘突变力矩的可接受程度却因人而异。
发明内容
为了解决现有紧急避撞方法存在的在避撞过程中没有考虑运动障碍物位置的动态变化而导致避撞过程不安全的技术问题,以及现有紧急避撞方法存在的转向盘突变力矩不可控容易导致驾驶员误操作的技术问题,本发明提供一种躲避动障碍物的汽车紧急避撞一体式控制方法,能够辅助驾驶员实现安全可靠的避撞,在紧急关头拯救驾驶员生命。
本发明解决技术问题所采取的技术方案如下:
1、一种躲避动障碍物的汽车紧急避撞一体式控制方法,其特征在于,该方法是利用路径动态规划与实时跟踪控制模块,根据实时采集的障碍物信息、目标点坐标、汽车行驶状态信息,实时优化得出汽车的前轮转角和四个车轮滑移率,进而控制汽车实现避撞;其中,障碍物信息包括由雷达传感器测量获得的障碍物外形轮廓的离散点坐标,汽车行驶状态信息包括由车速传感器测量获得的汽车纵向速度和侧向速度以及由陀螺仪测量获得的横摆角速度;在控制避撞过程中,通过电动助力转向(Electric Power Steering,EPS)力矩补偿模块根据车速、前轮附加转角,确定力矩补偿控制增益,将转向盘突变力矩控制在理想范围内,实现人机和谐的汽车紧急避撞;
该方法包括以下步骤:
步骤1、汽车紧急避撞控制的性能指标设计过程包括如下子步骤:
步骤1.1、用预测时域内预测轨迹的终点坐标与目标点坐标误差的二范数作为跟踪性能指标,体现汽车的轨迹跟踪特性,其表达式如下:
其中,Hp为预测时域,(Xt+Hp,Yt+Hp)为预测时域内预测轨迹的终点坐标,由汽车模型迭代获得,避撞时汽车要达到的目标点坐标(Xg,Yg);
所述汽车动力学模型为:
Fxi=fxicos(δi)-fyisin(δi),i∈{1,2,3,4} (15)
Fyi=fxisin(δi)+fyicos(δi),i∈{1,2,3,4} (16)
其中,Fxi、Fyi分别为四个车轮沿着车身坐标方向的纵向分力与侧向分力;fxi、fyi分别为四个车轮沿车轮坐标方向的分力,其中fxi为四个车轮滑移率和车轮垂直载荷的函数,fyi为前轮转角和车轮垂直载荷的函数,具体数值可由魔术公式确定;分别为汽车纵向速度和纵向加速度;分别为汽车侧向速度和侧向加速度;分别为汽车横摆角、横摆角速度和横摆角加速度;lf、lr分别为汽车质心到前、后轴的距离,ls为轮距大小的一半;Jz为绕过汽车质心的铅垂轴的横摆转动惯量;M为汽车质量;X、Y分别为大地坐标系中汽车质心位置的横纵坐标;δi为四个车轮转角,这里汽车为前轮转向,故δ3=δ4=0;
所述魔术公式的参数由试验拟合得出,具体表达式如下:
其中,V为当前汽车纵向速度;αf、αr分别为前轮侧偏角和后轮侧偏角;Fz,f、Fz,r分别为汽车前、后轴载荷;si为汽车四个车轮滑移率;Axi、Bxi、Cxi、Dxi、Exi和Ayi、Byi、Cyi、Dyi、Eyi为试验拟合参数,具体参数由下表所示:
表3魔术公式参数
a0 | a1 | a2 | a3 | a4 | a5 | a6 | ||
1.75 | 0 | 1000 | 1289 | 7.11 | 0.0053 | 0.1925 | ||
b0 | b1 | b2 | b3 | b4 | b5 | b6 | b7 | b8 |
1.57 | 35 | 1200 | 60 | 300 | 0.17 | 0 | 0 | 0.2 |
步骤1.2、用控制量变化率的二范数作为转向制动平滑指标,体现避撞过程中的执行器的转向制动平滑特性,控制量u为汽车前轮转角δ和汽车四个车轮滑移率si i∈{1,2,3,4},建立离散二次型转向制动平滑指标为:
其中,Hc为控制时域,t表示当前时刻,Δu为控制量变化率,w为Δu的权重系数;
步骤2、考虑运动障碍物的汽车紧急避撞控制的约束设计过程包括如下子步骤:
步骤2.1、设置执行器物理约束,满足执行器要求;
利用线性不等式限制前轮转角和四个车轮滑移率的上下限,分别得到转向、制动执行器的物理约束,其数学表达式为:
δmin<δk,t<δmax k=t,t+1……t+Hc-1 (3)
simin<sik,t<simax i∈{1,2,3,4}k=t,t+1……t+Hc-1 (4)
其中,δmin为前轮转角下限,δmax为前轮转角上限,simin为四个车轮滑移率下限,simax为四个车轮滑移率上限;
步骤2.2、设置位置约束,保证避撞过程中不会与障碍物碰撞;
t时刻障碍物的位置信息可表征为N个离散点的集合,这些信息可由雷达传感器测量获得,其中第j个离散点的坐标表示为(Xj,t,Yj,t),t时刻的汽车质心坐标记为(Xk,t,Yk,t),可由步骤1.1所述的汽车模型计算得出,位置约束定为
其中,a为汽车质心到车头的距离;b为汽车质心到车尾的距离;c为汽车车宽的一半;为以t时刻为起点预测时域内k时刻汽车的横摆角;Dx,j,t为障碍物第j个离散点在汽车坐标系中到汽车质心的纵向距离,Dy,j,t为障碍物第j个离散点在汽车坐标系中到汽车质心的横向距离;
假定在预测时域内障碍物沿Y方向以恒定速度运动,公式(5)表征了汽车与障碍物N个离散点的接近程度,值越大,说明汽车与障碍物对应离散点的距离越接近,也就越危险;定义t时刻值最大的障碍物离散点j为当前采样周期内的危险点,记为(Xj,t,Yj,t),在预测时域内基于此危险点对障碍物运动进行预测,迭代关系表示为:
其中,(Xj,t-1,Yj,t-1)为危险点在t-1时刻的坐标;(Xj,k,Yj,k)为预测时域内k时刻危险点的坐标;
通过迭代的方式更新公式(5)中的障碍物离散点坐标,在预测时域内将障碍物的位置变化整合为模型预测控制算法的位置约束;
步骤3、构建汽车紧急避撞多目标优化控制问题,求解多目标优化控制问题,以动态约束形式制定汽车行驶的不碰撞路径,实现考虑运动障碍物的汽车紧急避撞控制,其包括如下子步骤:
步骤3.1、通过雷达传感器获取障碍物信息,通过车速传感器和陀螺仪获取汽车行驶状态信息,并将获取的障碍物信息和汽车行驶状态信息输入避撞控制器;
步骤3.2、利用线性加权法将步骤1.1所述跟踪性能指标和步骤1.2所述转向制动平滑指标转化为单一指标,构建汽车紧急避撞多目标优化控制问题,该问题要同时满足转向、制动执行器的物理约束和位置约束,且保证紧急避撞系统输入输出符合步骤1.1所述的汽车动力学模型特性:
服从于
i)汽车动力学模型
ii)约束条件为公式(3)~(9)
步骤3.3、在紧急避撞控制器中,调用SQP算法,求解多目标优化控制问题(10),得到最优开环控制u*为:
服从于
i)汽车动力学模型
ii)约束条件为公式(3)~(9)
步骤3.4、利用当前时刻最优开环控制u*(0)进行反馈,实现闭环控制,实现了考虑运动障碍物的汽车紧急避撞控制。
步骤4、设计植有转向盘突变力矩人性化调节算法的EPS力矩补偿模块,EPS力矩补偿模块根据车速、前轮附加转角,确定力矩补偿控制增益,将转向盘突变力矩控制在理想范围;其中,前轮附加转角是路径动态规划与实时跟踪控制模块优化出的前轮转角与驾驶员转向输入产生的前轮转角的差值,由AFS系统来实现;设计过程包括如下子步骤:
步骤4.1、EPS力矩补偿模块的设计方法为:选取多名驾驶员进行实车调试,首先通过调试定车速、定前轮附加转角下的力矩补偿控制增益,实验员根据驾驶员的主观感受进行反复调试,保证转向盘突变力矩能够被驾驶员接受;
步骤4.2、改变前轮附加转角,实验员调试力矩补偿控制增益使不同前轮附加转角干预下的转向盘突变力矩均可被驾驶员接受,进而确定出该车速下的力矩补偿控制增益;
步骤4.3、采用相同的方法确定出不同车速、不同前轮附加转角干预下的力矩补偿控制增益,完成车速、前轮附加转角、力矩补偿控制增益三维MAP图的确定,使用力矩补偿控制增益三维数表进行力矩补偿控制,将转向盘突变力矩控制在理想范围内,实现转向盘突变力矩人性化调节的汽车紧急避撞。
步骤4.4、将EPS力矩补偿控制增益三维MAP图植入EPS控制器中,EPS控制器控制EPS助力电机达到力矩补偿的控制效果。
本发明的有益效果是:该方法通过构建多目标优化问题,解决了汽车紧急避撞时路径动态规划与实时跟踪问题,并同时考虑了动障碍物的情况,实现了安全最优避撞。该方法基于模型预测控制构建多目标优化问题,再以动态约束的形式制定了不碰撞路径,与传统分层式方法相比实时性高,路径满足动力学约束,避撞过程更为可靠。该方法通过障碍物坐标变化的方式,将障碍物运动情况转化为避撞控制优化求解的动态约束,解决了避障控制中的运动障碍物问题;该方法通过EPS力矩补偿控制器,将转向盘突变力矩控制在驾驶员可接受的范围,本方法使用主观评测的方式对EPS力矩补偿控制增益进行反复调试,实现了人性化突变力矩调节。
附图说明
图1是本发明躲避运动障碍物的汽车紧急避撞控制一体式方法的原理示意图。
图2是汽车与障碍物位置关系示意图。
图3是汽车与障碍物运动关系示意图。
图4是本发明汽车模型图。
图5是本发明的EPS力矩补偿控制器实验流程。
图6是本发明EPS力矩补偿控制增益三维MAP图。
具体实施方式
下面结合附图和实例对本发明做进一步详细说明。
如图1所示,本发明一种躲避运动障碍物的汽车紧急避撞控制一体式方法是:路径动态规划与实时跟踪控制模块1根据实时采集的障碍物信息、目标点坐标、汽车行驶状态信息,实时优化得出汽车2的前轮转角和四个车轮滑移率,控制汽车2实现避撞;其中,障碍物信息包括障碍物外形轮廓的离散点坐标,由雷达传感器测量获得;汽车行驶状态信息包括汽车纵向速度、侧向速度、横摆角速度,汽车纵向速度和侧向速度由车速传感器测量获得,汽车横摆角速度由陀螺仪测量获得。在控制避撞过程中,通过EPS力矩补偿模块3根据车速、前轮附加转角,确定力矩补偿控制增益,将转向盘突变力矩控制在驾驶员4可接受理想范围内,实现人机和谐的汽车紧急避撞。
本发明中的路径动态规划与实时跟踪控制模块1包括三部分内容:1)汽车紧急避撞控制的性能指标设计;2)躲避运动障碍物的汽车紧急避撞控制的约束设计;3)控制律滚动时域求解。
下面以某轿车为平台,具体说明本发明的方法,试验轿车的主要参数如表1所示:
表1试验轿车的主要参数
在1)部分内容中,汽车紧急避撞控制的性能指标设计包括以下两部分:1.1、用预测时域内预测轨迹的终点坐标与目标点坐标误差的二范数作为跟踪性能指标,体现汽车的轨迹跟踪特性;1.2、利用控制量变化率的二范数作为转向制动平滑指标,体现执行器的转向制动平滑特性。
在1.1部分中,跟踪性能指标以预测时域内预测轨迹的终点坐标与目标点坐标误差的二范数为评价标准,表达式如下:
其中,Hp为预测时域,(Xt+Hp,Yt+Hp)为预测时域内预测轨迹的终点坐标,由汽车模型迭代获得,避撞时汽车要达到的目标点坐标(Xg,Yg),即障碍物后方一安全点。
在1.2部分中,利用控制量变化率的二范数描述避撞过程中的执行器的转向制动平滑特性,其中控制量u为汽车前轮转角δ和四个车轮滑移率si i∈{1,2,3,4},建立离散二次型转向制动平滑指标为:
其中,Hc为控制时域,t表示当前时刻,Δu为控制量变化率,w为Δu的权重系数。紧急避撞控制器设计参数如表2所示,其中Ts为采样周期。
表2紧急避撞控制器设计参数
控制器参数 | 参数值 | 控制器参数 | 参数值 |
Hp | 4 | δmin | -6deg |
w | 0.5 | δmax | 6deg |
Ts | 0.5s | simin | 0 |
Hc | 3 | simax | 0.25 |
在2)部分内容中,躲避运动障碍物的汽车紧急避撞控制的约束设计包括两部分:2.1、设置执行器物理约束,满足执行器要求;2.2、设置位置约束,保证避撞过程中不会与障碍物碰撞。
在2.1部分中,利用线性不等式限制前轮转角和四个车轮滑移率的上下限,分别得到转向、制动执行器的物理约束,其数学表达式为:
δmin<δk,t<δmax k=t,t+1……t+Hc-1 (3)
simin<sik,t<simax i∈{1,2,3,4}k=t,t+1……t+Hc-1 (4)
其中δmin为前轮转角下限,δmax为前轮转角上限,simin为四个车轮滑移率下限,simax为四个车轮滑移率上限。
在2.2部分中,如图2所示,t时刻障碍物的位置信息可表征为N个离散点的集合,这些信息可由雷达测量获得,其中第j个离散点的坐标表示为(Xj,t,Yj,t),t时刻的汽车质心坐标记为(Xk,t,Yk,t),可由汽车模型计算得出,位置约束定为
其中,a为汽车质心到车头的距离;b为汽车质心到车尾的距离;c为汽车车宽的一半;为以t时刻为起点预测时域内k时刻汽车的横摆角,;Dx,j,t为障碍物第j个离散点在汽车坐标系中到汽车质心的纵向距离;Dy,j,t为障碍物第j个离散点在汽车坐标系中到汽车质心的横向距离。
如图3所示,汽车行驶过程中,障碍物可能以运动的方式突然出现;考虑障碍物沿Y方向运动情况,假定在预测时域内障碍物以恒定速度运动。
公式(5)表征了汽车与障碍物N个离散点的接近程度,值越大,说明汽车与障碍物对应离散点的距离越接近,也就越危险。为了保证算法实时性,定义t时刻值最大的障碍物离散点j为当前采样周期内的危险点,记为(Xj,t,Yj,t),在预测时域内基于此危险点对障碍物运动进行预测,迭代关系表示为:
其中,(Xj,t-1,Yj,t-1)为危险点在t-1时刻的坐标;(Xj,k,Yj,k)为预测时域内k时刻危险点的坐标。
通过迭代的方式更新公式(5)中的障碍物离散点坐标,在预测时域内将障碍物的位置变化整合为模型预测控制算法的位置约束,优化求解运动障碍物下的紧急避撞问题。
在3)部分内容中,控制律滚动时域求解包括以下步骤:
3.1、从雷达和车载传感器获取障碍物信息和汽车行驶状态信息,并将信息输入避撞控制器;
3.2、利用线性加权法将跟踪性能指标和转向制动平滑指标转化为单一指标,构建紧急避撞多目标优化控制问题,该问题要同时满足转向、制动执行器的物理约束和位置约束,且保证紧急避撞系统输入输出符合汽车动力学模型特性:
服从于
i)汽车动力学模型
ii)约束条件为公式(3)~(9)
3.3、在紧急避撞控制器中,调用SQP算法,求解多目标优化控制问题(10),得到最优开环控制u*为:
服从于
i)汽车动力学模型
ii)约束条件为公式(3)~(9)
3.4、利用当前时刻最优开环控制u*(0)进行反馈,实现闭环控制;
如图4所示,本发明上述的汽车动力学模型为:
Fxi=fxicos(δi)-fyisin(δi),i∈{1,2,3,4} (15)
Fyi=fxisin(δi)+fyicos(δi),i∈{1,2,3,4} (16)
其中,Fxi、Fyi分别为四个车轮沿着车身坐标方向的纵向分力与侧向分力;fxi、fyi分别是四个车轮沿车轮坐标方向的分力,其中fxi为四个车轮滑移率和车轮垂直载荷的函数,fyi为前轮转角和车轮垂直载荷的函数,具体数值可由魔术公式确定;分别为汽车纵向速度和纵向加速度;分别为汽车侧向速度和侧向加速度;分别为汽车横摆角、横摆角速度和横摆角加速度;lf、lr分别为汽车质心到前、后轴的距离,ls为轮距大小的一半;Jz为绕过汽车质心的铅垂轴的横摆转动惯量;M为汽车质量;X、Y分别为大地坐标系中汽车质心位置的横纵坐标;δi为四个车轮转角,这里汽车为前轮转向,故δ3=δ4=0。
上述魔术公式的参数由试验拟合得出,具体表达式如下:
其中,V为当前汽车纵向速度;αf、αr分别为前轮侧偏角和后轮侧偏角;Fz,f、Fz,r分别为汽车前、后轴载荷;si为汽车四个车轮滑移率;Axi、Bxi、Cxi、Dxi、Exi和Ayi、Byi、Cyi、Dyi、Eyi为试验拟合参数,具体参数由下表所示:
表3魔术公式参数
a0 | a1 | a2 | a3 | a4 | a5 | a6 | ||
1.75 | 0 | 1000 | 1289 | 7.11 | 0.0053 | 0.1925 | ||
b0 | b1 | b2 | b3 | b4 | b5 | b6 | b7 | b8 |
1.57 | 35 | 1200 | 60 | 300 | 0.17 | 0 | 0 | 0.2 |
本发明中的EPS力矩补偿模块3的设计方法为:选取30名驾驶员,按照性别、熟练程度分为以下四类:熟练男驾驶员、熟练女驾驶员、不熟练男驾驶员、不熟练女驾驶员。驾驶员按照预先分类分别进行实车调试,调试流程如图5所示,首先将车速定为60km/h,前轮附加转角定为3deg,实验员根据驾驶员的对转向盘突变力矩接受程度的反馈信息,反复调试力矩补偿控制增益,当驾驶员感觉突变力矩过大时,实验员将力矩补偿控制增益降低,当驾驶员感觉突变力矩过小时,实验员则将力矩补偿控制增益调大,最终保证转向盘突变力矩能够被驾驶员所接受,并记录此时的力矩补偿控制增益数值;其次,车速仍定为60km/h,前轮附加转角范围为-6deg到6deg,间隔为2deg,由于汽车转向时左右两侧是对称的,前轮附加转角相同幅值情况下左右两侧产生的转向盘突变力矩是相同的,因此只需调整前轮附加转角范围为0deg到6deg下的力矩补偿控制增益即可。试验时实验员根据驾驶员对转向盘突变力矩的接受程度调试0deg到6deg范围内各个转角干预下的力矩补偿控制增益,使各前轮附加转角干预下的转向盘突变力矩均被驾驶员接受,进而确定出车速60km/h不同转角干预下的力矩补偿控制增益,并记录力矩补偿控制增益的具体数值;最后,采用相同的方法调试出不同车速不同转角干预下的力矩补偿控制增益,车速范围为10km/h到100km/h,车速间隔为20km/h,最终确定出车速、前轮附加转角、力矩补偿控制增益的三维数表,图6是本发明EPS力矩补偿控制增益三维MAP图。最后将EPS力矩补偿控制增益三维MAP图植入EPS控制器中,EPS控制器控制EPS助力电机达到力矩补偿的控制效果。
Claims (1)
1.一种躲避动障碍物的汽车紧急避撞一体式控制方法,其特征在于,该方法是利用路径动态规划与实时跟踪控制模块,根据实时采集的障碍物信息、目标点坐标、汽车行驶状态信息,实时优化得出汽车的前轮转角和四个车轮滑移率,进而控制汽车实现避撞;其中,障碍物信息包括由雷达传感器测量获得的障碍物外形轮廓的离散点坐标,汽车行驶状态信息包括由车速传感器测量获得的汽车纵向速度和侧向速度以及由陀螺仪测量获得的横摆角速度;在控制避撞过程中,通过EPS力矩补偿模块根据车速、前轮附加转角,确定力矩补偿控制增益,将转向盘突变力矩控制在理想范围内,实现人机和谐的汽车紧急避撞;
该方法包括以下步骤:
步骤1、汽车紧急避撞控制的性能指标设计过程包括如下子步骤:
步骤1.1、用预测时域内预测轨迹的终点坐标与目标点坐标误差的二范数作为跟踪性能指标,体现汽车的轨迹跟踪特性,其表达式如下:
<mrow>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>X</mi>
<mi>g</mi>
</msub>
<mo>-</mo>
<msub>
<mi>X</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<msub>
<mi>H</mi>
<mi>p</mi>
</msub>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>Y</mi>
<mi>g</mi>
</msub>
<mo>-</mo>
<msub>
<mi>Y</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<msub>
<mi>H</mi>
<mi>p</mi>
</msub>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,Hp为预测时域,(Xt+Hp,Yt+Hp)为预测时域内预测轨迹的终点坐标,由汽车模型迭代获得,避撞时汽车要达到的目标点坐标(Xg,Yg);
所述汽车动力学模型为:
Fxi=fxicos(δi)-fyisin(δi),i∈{1,2,3,4} (15)
Fyi=fxisin(δi)+fyicos(δi),i∈{1,2,3,4} (16)
其中,Fxi、Fyi分别为四个车轮沿着车身坐标方向的纵向分力与侧向分力;fxi、fyi分别为四个车轮沿车轮坐标方向的分力,其中fxi为四个车轮滑移率和车轮垂直载荷的函数,fyi为前轮转角和车轮垂直载荷的函数,具体数值可由魔术公式确定;分别为汽车纵向速度和纵向加速度;分别为汽车侧向速度和侧向加速度;分别为汽车横摆角、横摆角速度和横摆角加速度;lf、lr分别为汽车质心到前、后轴的距离,ls为轮距大小的一半;Jz为绕过汽车质心的铅垂轴的横摆转动惯量;M为汽车质量;X、Y分别为大地坐标系中汽车质心位置的横纵坐标;δi为四个车轮转角,这里汽车为前轮转向,故δ3=δ4=0;
所述魔术公式的参数由实验拟合得出,具体表达式如下:
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>z</mi>
<mo>,</mo>
<mi>f</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>Mgl</mi>
<mi>r</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>l</mi>
<mi>f</mi>
</msub>
<mo>+</mo>
<msub>
<mi>l</mi>
<mi>r</mi>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>21</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>z</mi>
<mo>,</mo>
<mi>r</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>Mgl</mi>
<mi>f</mi>
</msub>
</mrow>
<mrow>
<msub>
<mi>l</mi>
<mi>f</mi>
</msub>
<mo>+</mo>
<msub>
<mi>l</mi>
<mi>r</mi>
</msub>
</mrow>
</mfrac>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>22</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>f</mi>
<mrow>
<mi>x</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>D</mi>
<mrow>
<mi>x</mi>
<mi>i</mi>
</mrow>
</msub>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>x</mi>
<mi>i</mi>
</mrow>
</msub>
<mi>arctan</mi>
<mo>(</mo>
<mrow>
<msub>
<mi>A</mi>
<mrow>
<mi>x</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>E</mi>
<mrow>
<mi>x</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>A</mi>
<mrow>
<mi>x</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<mi>arctan</mi>
<mi> </mi>
<msub>
<mi>A</mi>
<mrow>
<mi>x</mi>
<mi>i</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>A</mi>
<mrow>
<mi>x</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>B</mi>
<mrow>
<mi>x</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>s</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
<mi>i</mi>
<mo>&Element;</mo>
<mo>{</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>}</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>23</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>C</mi>
<mrow>
<mi>x</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>b</mi>
<mn>0</mn>
</msub>
<mo>;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>D</mi>
<mrow>
<mi>x</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>b</mi>
<mn>1</mn>
</msub>
<msup>
<msub>
<mi>F</mi>
<mrow>
<mi>z</mi>
<mo>,</mo>
<mi>r</mi>
</mrow>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>2</mn>
</msub>
<msub>
<mi>F</mi>
<mrow>
<mi>z</mi>
<mo>,</mo>
<mi>r</mi>
</mrow>
</msub>
<mo>;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>B</mi>
<mrow>
<mi>x</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>b</mi>
<mn>3</mn>
</msub>
<msup>
<msub>
<mi>F</mi>
<mrow>
<mi>z</mi>
<mo>,</mo>
<mi>r</mi>
</mrow>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>4</mn>
</msub>
<msub>
<mi>F</mi>
<mrow>
<mi>z</mi>
<mo>,</mo>
<mi>r</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>C</mi>
<mrow>
<mi>x</mi>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>D</mi>
<mrow>
<mi>x</mi>
<mi>i</mi>
</mrow>
</msub>
<msup>
<mi>e</mi>
<mrow>
<msub>
<mi>b</mi>
<mn>5</mn>
</msub>
<msub>
<mi>F</mi>
<mrow>
<mi>z</mi>
<mo>,</mo>
<mi>r</mi>
</mrow>
</msub>
</mrow>
</msup>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>E</mi>
<mrow>
<mi>x</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>b</mi>
<mn>6</mn>
</msub>
<msub>
<mi>F</mi>
<mrow>
<mi>z</mi>
<mo>,</mo>
<mi>r</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>7</mn>
</msub>
<msub>
<mi>F</mi>
<mrow>
<mi>z</mi>
<mo>,</mo>
<mi>r</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>8</mn>
</msub>
<mo>;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mi>i</mi>
<mo>&Element;</mo>
<mo>{</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>}</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>24</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>f</mi>
<mrow>
<mi>y</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>D</mi>
<mrow>
<mi>y</mi>
<mi>i</mi>
</mrow>
</msub>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>y</mi>
<mi>i</mi>
</mrow>
</msub>
<mi>arctan</mi>
<mo>(</mo>
<mrow>
<msub>
<mi>A</mi>
<mrow>
<mi>y</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>E</mi>
<mrow>
<mi>y</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>A</mi>
<mrow>
<mi>y</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<mi>arctan</mi>
<mi> </mi>
<msub>
<mi>A</mi>
<mrow>
<mi>y</mi>
<mi>i</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>A</mi>
<mrow>
<mi>y</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>B</mi>
<mrow>
<mi>y</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>&alpha;</mi>
<mi>f</mi>
</msub>
<mo>,</mo>
<mi>i</mi>
<mo>&Element;</mo>
<mo>{</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>}</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>25</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>C</mi>
<mrow>
<mi>y</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>a</mi>
<mn>0</mn>
</msub>
<mo>;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>D</mi>
<mrow>
<mi>y</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<msup>
<msub>
<mi>F</mi>
<mrow>
<mi>z</mi>
<mo>,</mo>
<mi>f</mi>
</mrow>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<msub>
<mi>F</mi>
<mrow>
<mi>z</mi>
<mo>,</mo>
<mi>f</mi>
</mrow>
</msub>
<mo>;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>B</mi>
<mrow>
<mi>y</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>a</mi>
<mn>3</mn>
</msub>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>arctan</mi>
<mo>(</mo>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>z</mi>
<mo>,</mo>
<mi>f</mi>
</mrow>
</msub>
<mo>/</mo>
<msub>
<mi>a</mi>
<mn>4</mn>
</msub>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>C</mi>
<mrow>
<mi>y</mi>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>D</mi>
<mrow>
<mi>y</mi>
<mi>i</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>E</mi>
<mrow>
<mi>y</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>a</mi>
<mn>5</mn>
</msub>
<msub>
<mi>F</mi>
<mrow>
<mi>z</mi>
<mo>,</mo>
<mi>f</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>6</mn>
</msub>
<mo>;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mi>i</mi>
<mo>&Element;</mo>
<mo>{</mo>
<mn>1</mn>
<mo>,</mo>
<mn>2</mn>
<mo>}</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>26</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>f</mi>
<mrow>
<mi>x</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>D</mi>
<mrow>
<mi>x</mi>
<mi>i</mi>
</mrow>
</msub>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>x</mi>
<mi>i</mi>
</mrow>
</msub>
<mi>arctan</mi>
<mo>(</mo>
<mrow>
<msub>
<mi>A</mi>
<mrow>
<mi>x</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>E</mi>
<mrow>
<mi>x</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>A</mi>
<mrow>
<mi>x</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<mi>arctan</mi>
<mi> </mi>
<msub>
<mi>A</mi>
<mrow>
<mi>x</mi>
<mi>i</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>A</mi>
<mrow>
<mi>x</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>B</mi>
<mrow>
<mi>x</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>s</mi>
<mi>i</mi>
</msub>
<mo>,</mo>
<mi>i</mi>
<mo>&Element;</mo>
<mo>{</mo>
<mn>3</mn>
<mo>,</mo>
<mn>4</mn>
<mo>}</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>27</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>C</mi>
<mrow>
<mi>x</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>b</mi>
<mn>0</mn>
</msub>
<mo>;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>D</mi>
<mrow>
<mi>x</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>b</mi>
<mn>1</mn>
</msub>
<msup>
<msub>
<mi>F</mi>
<mrow>
<mi>z</mi>
<mo>,</mo>
<mi>r</mi>
</mrow>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>2</mn>
</msub>
<msub>
<mi>F</mi>
<mrow>
<mi>z</mi>
<mo>,</mo>
<mi>r</mi>
</mrow>
</msub>
<mo>;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>B</mi>
<mrow>
<mi>x</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>b</mi>
<mn>3</mn>
</msub>
<msup>
<msub>
<mi>F</mi>
<mrow>
<mi>z</mi>
<mo>,</mo>
<mi>r</mi>
</mrow>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>4</mn>
</msub>
<msub>
<mi>F</mi>
<mrow>
<mi>z</mi>
<mo>,</mo>
<mi>r</mi>
</mrow>
</msub>
</mrow>
<mrow>
<msub>
<mi>C</mi>
<mrow>
<mi>x</mi>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>D</mi>
<mrow>
<mi>x</mi>
<mi>i</mi>
</mrow>
</msub>
<msup>
<mi>e</mi>
<mrow>
<msub>
<mi>b</mi>
<mn>5</mn>
</msub>
<msub>
<mi>F</mi>
<mrow>
<mi>z</mi>
<mo>,</mo>
<mi>r</mi>
</mrow>
</msub>
</mrow>
</msup>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>E</mi>
<mrow>
<mi>x</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>b</mi>
<mn>6</mn>
</msub>
<msub>
<mi>F</mi>
<mrow>
<mi>z</mi>
<mo>,</mo>
<mi>r</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>7</mn>
</msub>
<msub>
<mi>F</mi>
<mrow>
<mi>z</mi>
<mo>,</mo>
<mi>r</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>b</mi>
<mn>8</mn>
</msub>
<mo>;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mi>i</mi>
<mo>&Element;</mo>
<mo>{</mo>
<mn>3</mn>
<mo>,</mo>
<mn>4</mn>
<mo>}</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>28</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>f</mi>
<mrow>
<mi>y</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mo>-</mo>
<msub>
<mi>D</mi>
<mrow>
<mi>y</mi>
<mi>i</mi>
</mrow>
</msub>
<mi>s</mi>
<mi>i</mi>
<mi>n</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>C</mi>
<mrow>
<mi>y</mi>
<mi>i</mi>
</mrow>
</msub>
<mi>arctan</mi>
<mo>(</mo>
<mrow>
<msub>
<mi>A</mi>
<mrow>
<mi>y</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>E</mi>
<mrow>
<mi>y</mi>
<mi>i</mi>
</mrow>
</msub>
<mrow>
<mo>(</mo>
<mrow>
<msub>
<mi>A</mi>
<mrow>
<mi>y</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>-</mo>
<mi>arctan</mi>
<mi> </mi>
<msub>
<mi>A</mi>
<mrow>
<mi>y</mi>
<mi>i</mi>
</mrow>
</msub>
</mrow>
<mo>)</mo>
</mrow>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>A</mi>
<mrow>
<mi>y</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>B</mi>
<mrow>
<mi>y</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>&CenterDot;</mo>
<msub>
<mi>&alpha;</mi>
<mi>r</mi>
</msub>
<mo>,</mo>
<mi>i</mi>
<mo>&Element;</mo>
<mo>{</mo>
<mn>3</mn>
<mo>,</mo>
<mn>4</mn>
<mo>}</mo>
</mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>29</mn>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>C</mi>
<mrow>
<mi>y</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>a</mi>
<mn>0</mn>
</msub>
<mo>;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>D</mi>
<mrow>
<mi>y</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>a</mi>
<mn>1</mn>
</msub>
<msup>
<msub>
<mi>F</mi>
<mrow>
<mi>z</mi>
<mo>,</mo>
<mi>f</mi>
</mrow>
</msub>
<mn>2</mn>
</msup>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>2</mn>
</msub>
<msub>
<mi>F</mi>
<mrow>
<mi>z</mi>
<mo>,</mo>
<mi>f</mi>
</mrow>
</msub>
<mo>;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>B</mi>
<mrow>
<mi>y</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<mfrac>
<mrow>
<msub>
<mi>a</mi>
<mn>3</mn>
</msub>
<mi>sin</mi>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mi>arctan</mi>
<mo>(</mo>
<mrow>
<msub>
<mi>F</mi>
<mrow>
<mi>z</mi>
<mo>,</mo>
<mi>f</mi>
</mrow>
</msub>
<mo>/</mo>
<msub>
<mi>a</mi>
<mn>4</mn>
</msub>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<msub>
<mi>C</mi>
<mrow>
<mi>y</mi>
<mi>i</mi>
</mrow>
</msub>
<msub>
<mi>D</mi>
<mrow>
<mi>y</mi>
<mi>i</mi>
</mrow>
</msub>
</mrow>
</mfrac>
<mo>;</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>E</mi>
<mrow>
<mi>y</mi>
<mi>i</mi>
</mrow>
</msub>
<mo>=</mo>
<msub>
<mi>a</mi>
<mn>5</mn>
</msub>
<msub>
<mi>F</mi>
<mrow>
<mi>z</mi>
<mo>,</mo>
<mi>f</mi>
</mrow>
</msub>
<mo>+</mo>
<msub>
<mi>a</mi>
<mn>6</mn>
</msub>
<mo>;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>,</mo>
<mi>i</mi>
<mo>&Element;</mo>
<mo>{</mo>
<mn>3</mn>
<mo>,</mo>
<mn>4</mn>
<mo>}</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>30</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,V为当前汽车纵向速度;αf、αr分别为前轮侧偏角和后轮侧偏角;Fz,f、Fz,r分别为汽车前、后轴载荷;si为汽车四个车轮滑移率;Axi、Bxi、Cxi、Dxi、Exi和Ayi、Byi、Cyi、Dyi、Eyi为实验拟合参数,具体参数由下表所示:
表3 魔术公式参数
步骤1.2、用控制量变化率的二范数作为转向制动平滑指标,体现避撞过程中的执行器的转向制动平滑特性,控制量u为汽车前轮转角δ和汽车四个车轮滑移率si i∈{1,2,3,4},建立离散二次型转向制动平滑指标为:
<mrow>
<mi>w</mi>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>t</mi>
<mo>+</mo>
<msub>
<mi>H</mi>
<mi>c</mi>
</msub>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>&Delta;u</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>|</mo>
<msub>
<mo>|</mo>
<mn>2</mn>
</msub>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>2</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,Hc为控制时域,t表示当前时刻,Δu为控制量变化率,w为Δu的权重系数;
步骤2、考虑运动障碍物的汽车紧急避撞控制的约束设计过程包括如下子步骤:
步骤2.1、设置执行器物理约束,满足执行器要求;
利用线性不等式限制前轮转角和四个车轮滑移率的上下限,分别得到转向、制动执行器的物理约束,其数学表达式为:
δmin<δk,t<δmax k=t,t+1……t+Hc-1 (3)
simin<sik,t<simax i∈{1,2,3,4}k=t,t+1……t+Hc-1 (4)
其中,δmin为前轮转角下限,δmax为前轮转角上限,simin为四个车轮滑移率下限,simax为四个车轮滑移率上限;
步骤2.2、设置位置约束,保证避撞过程中不会与障碍物碰撞;
t时刻障碍物的位置信息可表征为N个离散点的集合,这些信息可由雷达传感器测量获得,其中第j个离散点的坐标表示为(Xj,t,Yj,t),t时刻的汽车质心坐标记为(Xk,t,Yk,t),可由步骤1.1所述的汽车模型计算得出,位置约束定为
<mrow>
<mi>l</mi>
<mo>=</mo>
<mfenced open = "{" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<mi>a</mi>
<mo>-</mo>
<mo>|</mo>
<msub>
<mi>D</mi>
<mrow>
<mi>x</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>|</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mi>f</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<mo>(</mo>
<msub>
<mi>D</mi>
<mrow>
<mi>x</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>></mo>
<mi>a</mi>
<mo>)</mo>
<mo>&cup;</mo>
<mo>(</mo>
<msub>
<mi>D</mi>
<mrow>
<mi>x</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo><</mo>
<mo>-</mo>
<mi>b</mi>
<mo>)</mo>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mi>c</mi>
<mo>-</mo>
<mo>|</mo>
<msub>
<mi>D</mi>
<mrow>
<mi>y</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>|</mo>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>i</mi>
<mi>f</mi>
</mrow>
</mtd>
<mtd>
<mrow>
<msub>
<mi>D</mi>
<mrow>
<mi>x</mi>
<mo>,</mo>
<mi>j</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>&Element;</mo>
<mo>&lsqb;</mo>
<mo>-</mo>
<mi>b</mi>
<mo>,</mo>
<mi>a</mi>
<mo>&rsqb;</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>6</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,a为汽车质心到车头的距离;b为汽车质心到车尾的距离;c为汽车车宽的一半;为以t时刻为起点预测时域内k时刻汽车的横摆角;Dx,j,t为障碍物第j个离散点在汽车坐标系中到汽车质心的纵向距离,Dy,j,t为障碍物第j个离散点在汽车坐标系中到汽车质心的横向距离;
假定在预测时域内障碍物沿Y方向以恒定速度运动,公式(5)表征了汽车与障碍物N个离散点的接近程度,值越大,说明汽车与障碍物对应离散点的距离越接近,也就越危险;定义t时刻值最大的障碍物离散点j为当前采样周期内的危险点,记为(Xj,t,Yj,t),在预测时域内基于此危险点对障碍物运动进行预测,迭代关系表示为:
<mrow>
<mtable>
<mtr>
<mtd>
<mrow>
<mo>(</mo>
<msub>
<mi>X</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>Y</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
</mrow>
</msub>
<mo>)</mo>
<mo>=</mo>
<mo>(</mo>
<msub>
<mi>X</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>Y</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>k</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
<mo>+</mo>
<msub>
<mi>&Delta;d</mi>
<mi>t</mi>
</msub>
</mrow>
</mtd>
<mtd>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mi>t</mi>
<mo>,</mo>
<mi>t</mi>
<mo>+</mo>
<mn>1</mn>
<mo>,</mo>
<mo>...</mo>
<mo>,</mo>
<mi>t</mi>
<mo>+</mo>
<msub>
<mi>H</mi>
<mi>p</mi>
</msub>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<msub>
<mi>&Delta;d</mi>
<mi>t</mi>
</msub>
<mo>=</mo>
<mrow>
<mo>(</mo>
<msub>
<mi>X</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>X</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>t</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>,</mo>
<msub>
<mi>Y</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>-</mo>
<msub>
<mi>X</mi>
<mrow>
<mi>j</mi>
<mo>,</mo>
<mi>t</mi>
<mo>-</mo>
<mn>1</mn>
</mrow>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
<mtd>
<mrow></mrow>
</mtd>
</mtr>
</mtable>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>9</mn>
<mo>)</mo>
</mrow>
</mrow>
其中,(Xj,t-1,Yj,t-1)为危险点在t-1时刻的坐标;(Xj,k,Yj,k)为预测时域内k时刻危险点的坐标;
通过迭代的方式更新公式(5)中的障碍物离散点坐标,在预测时域内将障碍物的位置变化整合为模型预测控制算法的位置约束;
步骤3、构建汽车紧急避撞多目标优化控制问题,求解多目标优化控制问题,以动态约束形式制定汽车行驶的不碰撞路径,实现考虑运动障碍物的汽车紧急避撞控制,其包括如下子步骤:
步骤3.1、通过雷达传感器获取障碍物信息,通过车速传感器和陀螺仪获取汽车行驶状态信息,并将获取的障碍物信息和汽车行驶状态信息输入避撞控制器;
步骤3.2、利用线性加权法将步骤1.1所述跟踪性能指标和步骤1.2所述转向制动平滑指标转化为单一指标,构建汽车紧急避撞多目标优化控制问题,该问题要同时满足转向、制动执行器的物理约束和位置约束,且保证紧急避撞系统输入输出符合步骤1.1所述的汽车动力学模型特性:
<mrow>
<munder>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
<mi>u</mi>
</munder>
<mo>{</mo>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>X</mi>
<mi>g</mi>
</msub>
<mo>-</mo>
<msub>
<mi>X</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<msub>
<mi>H</mi>
<mi>p</mi>
</msub>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>Y</mi>
<mi>g</mi>
</msub>
<mo>-</mo>
<msub>
<mi>Y</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<msub>
<mi>H</mi>
<mi>p</mi>
</msub>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>+</mo>
<mi>w</mi>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>t</mi>
<mo>+</mo>
<msub>
<mi>H</mi>
<mi>c</mi>
</msub>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>&Delta;u</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>|</mo>
<msub>
<mo>|</mo>
<mn>2</mn>
</msub>
<mo>}</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>10</mn>
<mo>)</mo>
</mrow>
</mrow>
服从于
i)汽车动力学模型
ii)约束条件为公式(3)~(9)
步骤3.3、在紧急避撞控制器中,调用SQP算法,求解多目标优化控制问题(10),得到最优开环控制u*为:
<mrow>
<msup>
<mi>u</mi>
<mo>*</mo>
</msup>
<mo>=</mo>
<mi>arg</mi>
<munder>
<mrow>
<mi>m</mi>
<mi>i</mi>
<mi>n</mi>
</mrow>
<mi>u</mi>
</munder>
<mo>{</mo>
<msqrt>
<mrow>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>X</mi>
<mi>g</mi>
</msub>
<mo>-</mo>
<msub>
<mi>X</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<msub>
<mi>H</mi>
<mi>p</mi>
</msub>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mo>+</mo>
<msup>
<mrow>
<mo>(</mo>
<msub>
<mi>Y</mi>
<mi>g</mi>
</msub>
<mo>-</mo>
<msub>
<mi>Y</mi>
<mrow>
<mi>t</mi>
<mo>+</mo>
<msub>
<mi>H</mi>
<mi>p</mi>
</msub>
</mrow>
</msub>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
</mrow>
</msqrt>
<mo>+</mo>
<mi>w</mi>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>k</mi>
<mo>=</mo>
<mi>t</mi>
</mrow>
<mrow>
<mi>t</mi>
<mo>+</mo>
<msub>
<mi>H</mi>
<mi>c</mi>
</msub>
<mo>-</mo>
<mn>1</mn>
</mrow>
</munderover>
<mo>|</mo>
<mo>|</mo>
<msub>
<mi>&Delta;u</mi>
<mrow>
<mi>k</mi>
<mo>,</mo>
<mi>t</mi>
</mrow>
</msub>
<mo>|</mo>
<msub>
<mo>|</mo>
<mn>2</mn>
</msub>
<mo>}</mo>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>11</mn>
<mo>)</mo>
</mrow>
</mrow>
服从于
i)汽车动力学模型
ii)约束条件为公式(3)~(9)
步骤3.4、利用当前时刻最优开环控制u*(0)进行反馈,实现闭环控制,实现了躲避运动障碍物的汽车紧急避撞控制。
步骤4、设计植有转向盘突变力矩人性化调节算法的EPS力矩补偿模块,EPS力矩补偿模块根据车速、前轮附加转角,确定力矩补偿控制增益,将转向盘突变力矩控制在理想范围;设计过程包括如下子步骤:
步骤4.1、EPS力矩补偿模块的设计方法为:选取多名驾驶员进行实车调试,首先通过调试定车速、定前轮附加转角下的力矩补偿控制增益,实验员根据驾驶员的主观感受进行反复调试,保证转向盘突变力矩能够被驾驶员接受;
步骤4.2、改变前轮附加转角,实验员调试力矩补偿控制增益使不同前轮附加转角干预下的转向盘突变力矩均可被驾驶员接受,进而确定出该车速下的力矩补偿控制增益;
步骤4.3、采用相同的方法确定出不同车速、不同前轮附加转角干预下的力矩补偿控制增益,完成车速、前轮附加转角、力矩补偿控制增益三维MAP图的确定,使用力矩补偿控制增益三维数表进行力矩补偿控制,将转向盘突变力矩控制在理想范围内,实现转向盘突变力矩人性化调节的汽车紧急避撞。
步骤4.4、将EPS力矩补偿控制增益三维MAP图植入EPS控制器中,EPS控制器控制EPS助力电机达到力矩补偿的控制效果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711081210.8A CN107878453B (zh) | 2017-11-07 | 2017-11-07 | 一种躲避动障碍物的汽车紧急避撞一体式控制方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711081210.8A CN107878453B (zh) | 2017-11-07 | 2017-11-07 | 一种躲避动障碍物的汽车紧急避撞一体式控制方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107878453A true CN107878453A (zh) | 2018-04-06 |
CN107878453B CN107878453B (zh) | 2019-07-30 |
Family
ID=61779047
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711081210.8A Expired - Fee Related CN107878453B (zh) | 2017-11-07 | 2017-11-07 | 一种躲避动障碍物的汽车紧急避撞一体式控制方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107878453B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109808687A (zh) * | 2019-01-17 | 2019-05-28 | 福瑞泰克智能系统有限公司 | 一种驾驶辅助方法、装置及设备 |
CN112068445A (zh) * | 2020-09-23 | 2020-12-11 | 北京理工大学 | 自动驾驶车辆路径规划与路径跟踪集成控制方法及系统 |
CN114071013A (zh) * | 2021-10-13 | 2022-02-18 | 浙江大华技术股份有限公司 | 一种用于车载摄像机的目标抓拍与跟踪方法及装置 |
WO2024152489A1 (zh) * | 2023-01-18 | 2024-07-25 | 清华大学 | 智能车辆紧急避撞方法与系统 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010023775A (ja) * | 2008-07-23 | 2010-02-04 | Toyota Motor Corp | 車両用動力伝達装置の制御装置 |
CN104977933A (zh) * | 2015-07-01 | 2015-10-14 | 吉林大学 | 一种自主驾驶车辆的区域型路径跟踪控制方法 |
CN105313891A (zh) * | 2015-10-22 | 2016-02-10 | 清华大学 | 一种多车协同避撞方法及装置 |
CN105539586A (zh) * | 2014-08-29 | 2016-05-04 | 通用汽车环球科技运作有限责任公司 | 用于自主驾驶的车辆躲避移动障碍物的统一的运动规划 |
CN105691388A (zh) * | 2016-01-14 | 2016-06-22 | 南京航空航天大学 | 一种汽车主动避撞系统及其轨迹规划方法 |
CN105930611A (zh) * | 2016-05-10 | 2016-09-07 | 长春工业大学 | 一种汽车电动助力转向特征车速区间划分方法 |
US20170003686A1 (en) * | 2015-07-03 | 2017-01-05 | Commissariat A L'energie Atomique Et Aux Energies Alternatives | Automatic control method for the insertion and the extraction of a vehicle into and from a receiving station, and control device implementing a method of this kind |
CN106564495A (zh) * | 2016-10-19 | 2017-04-19 | 江苏大学 | 融合空间和动力学特性的智能车辆安全驾驶包络重构方法 |
-
2017
- 2017-11-07 CN CN201711081210.8A patent/CN107878453B/zh not_active Expired - Fee Related
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2010023775A (ja) * | 2008-07-23 | 2010-02-04 | Toyota Motor Corp | 車両用動力伝達装置の制御装置 |
CN105539586A (zh) * | 2014-08-29 | 2016-05-04 | 通用汽车环球科技运作有限责任公司 | 用于自主驾驶的车辆躲避移动障碍物的统一的运动规划 |
CN104977933A (zh) * | 2015-07-01 | 2015-10-14 | 吉林大学 | 一种自主驾驶车辆的区域型路径跟踪控制方法 |
US20170003686A1 (en) * | 2015-07-03 | 2017-01-05 | Commissariat A L'energie Atomique Et Aux Energies Alternatives | Automatic control method for the insertion and the extraction of a vehicle into and from a receiving station, and control device implementing a method of this kind |
CN105313891A (zh) * | 2015-10-22 | 2016-02-10 | 清华大学 | 一种多车协同避撞方法及装置 |
CN105691388A (zh) * | 2016-01-14 | 2016-06-22 | 南京航空航天大学 | 一种汽车主动避撞系统及其轨迹规划方法 |
CN105930611A (zh) * | 2016-05-10 | 2016-09-07 | 长春工业大学 | 一种汽车电动助力转向特征车速区间划分方法 |
CN106564495A (zh) * | 2016-10-19 | 2017-04-19 | 江苏大学 | 融合空间和动力学特性的智能车辆安全驾驶包络重构方法 |
Non-Patent Citations (1)
Title |
---|
杨莉; 吴晓东; 许敏; 叶昌: "基于位置-力矩混合方法的线控转向系统双向控制", 《汽车技术》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109808687A (zh) * | 2019-01-17 | 2019-05-28 | 福瑞泰克智能系统有限公司 | 一种驾驶辅助方法、装置及设备 |
CN109808687B (zh) * | 2019-01-17 | 2021-09-21 | 福瑞泰克智能系统有限公司 | 一种驾驶辅助方法、装置及设备 |
CN112068445A (zh) * | 2020-09-23 | 2020-12-11 | 北京理工大学 | 自动驾驶车辆路径规划与路径跟踪集成控制方法及系统 |
US11634146B2 (en) | 2020-09-23 | 2023-04-25 | Beijing Institute Of Technology | Method and system for integrated path planning and path tracking control of autonomous vehicle |
CN114071013A (zh) * | 2021-10-13 | 2022-02-18 | 浙江大华技术股份有限公司 | 一种用于车载摄像机的目标抓拍与跟踪方法及装置 |
WO2024152489A1 (zh) * | 2023-01-18 | 2024-07-25 | 清华大学 | 智能车辆紧急避撞方法与系统 |
Also Published As
Publication number | Publication date |
---|---|
CN107878453B (zh) | 2019-07-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113386795B (zh) | 一种自动驾驶车辆智能决策及局部轨迹规划方法及其决策系统 | |
CN110286681B (zh) | 一种变曲率弯道的动态自动驾驶换道轨迹规划方法 | |
CN107839683B (zh) | 一种考虑运动障碍物的汽车紧急避撞控制方法 | |
CN106004870B (zh) | 一种基于变权重模型预测算法的车辆稳定性集成控制方法 | |
Attia et al. | Coupled longitudinal and lateral control strategy improving lateral stability for autonomous vehicle | |
US20210213935A1 (en) | Safety and Stability Control Method against Vehicle Tire Burst | |
CN107867290B (zh) | 一种考虑运动障碍物的汽车紧急避撞分层式控制方法 | |
CN107885932B (zh) | 一种考虑人机和谐的汽车紧急避撞分层式控制方法 | |
Cai et al. | Implementation and development of a trajectory tracking control system for intelligent vehicle | |
Wnag et al. | Path planning and stability control of collision avoidance system based on active front steering | |
CN104176054B (zh) | 一种汽车主动防碰撞自动换道控制系统及其工作方法 | |
CN107878453B (zh) | 一种躲避动障碍物的汽车紧急避撞一体式控制方法 | |
CN107856737A (zh) | 一种基于危险程度变权重的人机协同转向控制方法 | |
CN107380162B (zh) | 基于功能分配与多目标模糊决策的协同避撞方法 | |
Hongbo et al. | Lane‐keeping control based on an improved artificial potential method and coordination of steering/braking systems | |
Smith et al. | Effects of model complexity on the performance of automated vehicle steering controllers: controller development and evaluation | |
Dandiwala et al. | Vehicle dynamics and active rollover stability control of an electric narrow three-wheeled vehicle: A review and concern towards improvement | |
Mohajer et al. | Effects of road path profiles on autonomous vehicles’ handling behaviour | |
CN107856733B (zh) | 一种面向人机和谐的汽车躲避动障碍物控制方法 | |
Yue et al. | Automated hazard escaping trajectory planning/tracking control framework for vehicles subject to tire blowout on expressway | |
Xu et al. | Local path planning and tracking control of vehicle collision avoidance system | |
CN107885931A (zh) | 一种转向盘突变力矩人性化调节的汽车紧急避撞控制方法 | |
Zhu et al. | Trajectory planning algorithm using Gauss pseudo-spectral method based on vehicle-infrastructure cooperative system | |
CN110516353B (zh) | 一种山区高速公路弯道设计缺陷快速识别方法 | |
Lowe et al. | FMVSS 126 sine with dwell ESC regulation test for autonomous 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 | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20190730 Termination date: 20201107 |
|
CF01 | Termination of patent right due to non-payment of annual fee |