CN114510063B - 一种无人履带车辆及其轨迹跟踪控制方法和系统 - Google Patents
一种无人履带车辆及其轨迹跟踪控制方法和系统 Download PDFInfo
- Publication number
- CN114510063B CN114510063B CN202210352604.7A CN202210352604A CN114510063B CN 114510063 B CN114510063 B CN 114510063B CN 202210352604 A CN202210352604 A CN 202210352604A CN 114510063 B CN114510063 B CN 114510063B
- Authority
- CN
- China
- Prior art keywords
- tracked vehicle
- vehicle
- unmanned
- unmanned tracked
- track
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 45
- 230000004044 response Effects 0.000 claims abstract description 18
- 238000013528 artificial neural network Methods 0.000 claims abstract description 12
- 238000004804 winding Methods 0.000 claims description 21
- 230000001133 acceleration Effects 0.000 claims description 18
- 230000005484 gravity Effects 0.000 claims description 15
- 239000000126 substance Substances 0.000 claims description 15
- 239000011159 matrix material Substances 0.000 claims description 8
- 238000004458 analytical method Methods 0.000 claims description 6
- 238000002474 experimental method Methods 0.000 claims description 6
- 238000006243 chemical reaction Methods 0.000 claims description 5
- 239000000463 material Substances 0.000 claims 1
- 238000004364 calculation method Methods 0.000 abstract description 6
- 238000012545 processing Methods 0.000 abstract description 2
- 230000006870 function Effects 0.000 description 44
- 239000002245 particle Substances 0.000 description 9
- 238000004422 calculation algorithm Methods 0.000 description 8
- 230000008569 process Effects 0.000 description 7
- 238000005457 optimization Methods 0.000 description 6
- 230000008859 change Effects 0.000 description 5
- 230000003044 adaptive effect Effects 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000000605 extraction Methods 0.000 description 3
- 238000012549 training Methods 0.000 description 3
- ORILYTVJVMAKLC-UHFFFAOYSA-N Adamantane Natural products C1C(C2)CC3CC1CC2C3 ORILYTVJVMAKLC-UHFFFAOYSA-N 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 230000007246 mechanism Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000013468 resource allocation Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000012795 verification Methods 0.000 description 2
- 230000004913 activation Effects 0.000 description 1
- 230000006399 behavior Effects 0.000 description 1
- 238000005452 bending Methods 0.000 description 1
- 230000002146 bilateral effect Effects 0.000 description 1
- 230000001149 cognitive effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 230000009347 mechanical transmission Effects 0.000 description 1
- 238000003062 neural network model Methods 0.000 description 1
- -1 oxy Chemical group 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000007306 turnover Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course or altitude of land, water, air, or space vehicles, e.g. automatic pilot
- G05D1/02—Control of position or course in two dimensions
- G05D1/021—Control of position or course in two dimensions specially adapted to land vehicles
- G05D1/0212—Control of position or course in two dimensions specially adapted to land vehicles with means for defining a desired trajectory
- G05D1/0221—Control of position or course in two dimensions specially adapted to land vehicles with means for defining a desired trajectory involving a learning process
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05D—SYSTEMS FOR CONTROLLING OR REGULATING NON-ELECTRIC VARIABLES
- G05D1/00—Control of position, course or altitude of land, water, air, or space vehicles, e.g. automatic pilot
- G05D1/02—Control of position or course in two dimensions
- G05D1/021—Control of position or course in two dimensions specially adapted to land vehicles
- G05D1/0212—Control of position or course in two dimensions specially adapted to land vehicles with means for defining a desired trajectory
- G05D1/0223—Control of position or course in two dimensions specially adapted to land vehicles with means for defining a desired trajectory involving speed control of the vehicle
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/10—Internal combustion engine [ICE] based vehicles
- Y02T10/40—Engine management systems
Abstract
本发明涉及一种无人履带车辆及其轨迹跟踪控制方法和系统,属于车辆轨迹跟踪控制技术领域。本发明先基于无人履带车辆的运动学模型生成的轨迹跟踪控制器的目标函数,以及基于这一函数和无人履带车辆的动力学模型生成的轨迹跟踪控制器的约束函数形成MPC模型,然后,根据无人履带车辆的位姿特征、道路的曲率特征和实时工况状态信息生成参考轨迹后,将参考轨迹和车辆轨迹跟踪响应状态信息输入MPC模型得到参数组合,接着,将参数组合输入MLP神经网络得到控制参数,最后,根据控制参数完成无人履带车辆轨迹的跟踪控制,进而有效的提升无人履带车辆轨迹跟踪的工况适应性,实现轨迹跟踪精度、车辆行驶稳定性以及计算时间成本的较好均衡。
Description
技术领域
本发明涉及车辆轨迹跟踪控制技术领域,特别是涉及一种无人履带车辆及其轨迹跟踪控制方法和系统。
背景技术
近年来,地面无人平台在港口运输、疫情防控、野外搜救以及战场环境中得到了越来越多的应用,轨迹跟踪控制作为无人驾驶的关键技术,其发展也得到了越来越多学者的关注,这其中以模型预测的控制方法应用较多。然而,目前大多数都是基于定参数的轨迹跟踪控制方法,很难保证车辆在复杂多变工况下的跟踪精度和稳定性。基于此,有学者提出了自适应参数调整的方法,包括更精确的车辆模型参数在线估计、更合理的控制系统结构设计以及更高效的优化算法架构设计等。由于越野工况下地面-履带的相互作用关系难以精确表征,因此基于车辆模型参数在线估计的方法很难为系统提供稳定、可靠的输入,且在线估计对系统的实时资源配置要求较高。部分研究成果以MPC为基本框架,通过增加前馈、数据学习等方式提升轨迹跟踪性能,但是缺乏在实际高速履带平台的验证。
发明内容
为解决现有技术存在的上述问题,本发明提供了一种无人履带车辆及其轨迹跟踪控制方法和系统。
为实现上述目的,本发明提供了如下方案:
一种无人履带车辆的轨迹跟踪控制方法,包括:
获取无人履带车辆的运动学模型和无人履带车辆的动力学模型;
基于所述无人履带车辆的运动学模型生成轨迹跟踪控制器的目标函数;
基于所述无人履带车辆的动力学模型和所述轨迹跟踪控制器的目标函数生成轨迹跟踪控制器的约束函数;
基于所述轨迹跟踪控制器的目标函数和所述轨迹跟踪控制器的约束函数生成MPC模型;
根据无人履带车辆的位姿特征、道路的曲率特征和实时工况状态信息生成参考轨迹;
将所述参考轨迹和车辆轨迹跟踪响应状态信息输入所述MPC模型得到参数组合;所述车辆轨迹跟踪响应状态信息通过实验获得;所述参数组合包括:轨迹曲率特征向量、无人履带车辆姿态参考向量和无人履带车辆速度信息向量;
将所述参数组合输入MLP神经网络得到控制参数;
根据所述控制参数完成无人履带车辆轨迹的跟踪控制。
优选地,基于传统的无人履带车辆运动学模型,通过对无人履带车辆接地段的运动学分析,建立了基于瞬时转向中心模型的无人履带车辆的运动学模型。
优选地,所述基于瞬时转向中心模型的无人履带车辆的运动学模型为:
其中,ξ为无人履带车辆的状态向量,为无人履带车辆的横向速度,为无人履带车辆的纵向速度,为无人履带车辆的横摆角速度,v j1为一侧履带的卷绕速度,v j2为另一侧履带的卷绕速度,φ为履带车辆的几何中心在固连于地面的笛卡尔坐标系下的航向角,为一侧履带的横向方向,为另一侧履带的横向方向。
优选地,所述无人履带车辆的动力学模型为:
其中,x c 为无人履带车辆的重心与几何中心的横向距离,y c 为无人履带车辆的重心与几何中心的纵向距离,F cx 为离心力的横向分量、F cy 为离心力的纵向分量,F r1为一侧履带接地段所受的地面变形阻力,F r2为另一侧履带接地段所受的地面变形阻力,F y1为一侧履带制动力,F y2为另一侧履带驱动力,F x1为一侧履带所受的横向阻力的合力、F x2为另一侧履带所受的横向阻力的合力,m为无人履带车辆的整车质量,δ为无人履带车辆旋转质量换算系数,M μ 为转向阻力矩,M μ =M μ1+M μ2,M μ1为一侧履带所受的转向阻力矩、M μ2为另一侧履带所受的转向阻力矩,J为无人履带车辆绕z’轴的转动惯量,B为履带中心距,为车辆的横向加速度,为车辆的纵向加速度,为车辆的横摆角加速度。
优选地,所述轨迹跟踪控制器的目标函数为:
其中,J(ξ(k),U k )为轨迹跟踪控制器的目标函数,ξ(k)为k时刻的车辆状态向量,U k 为控制序列,N p 为预测时域的个数、N c 为控制时域的个数,且N c ≤N p ,i为个数变量,l e (k+i)为k+i时刻的轨迹跟踪的横向偏差,h e (k+i)为k+i时刻的轨迹跟踪的航向偏差,v(k+i)为k+i时刻的实际车速,l desired (k+i)为k+i时刻的期望横向偏差,h desired (k+i)为k+i时刻的期望航向偏差,v desired (k+i)为k+i时刻的期望车速,Q、R、P、M和N均为权重系数矩阵,v j1(*)为*时刻的一侧履带的卷绕速度,v j2(*)为*时刻的另一侧履带的卷绕速度。
优选地,所述轨迹跟踪控制器的约束函数包括:车辆纵向速度约束、转向中心纵向偏移量约束、电机转速约束和电机功率约束。
对应于上述提供的无人履带车辆的轨迹跟踪控制方法,本发明还提供了一种无人履带车辆的轨迹跟踪控制系统,该轨迹跟踪控制系统包括:
模型获取模块,用于获取无人履带车辆的运动学模型和无人履带车辆的动力学模型;
目标函数生成模块,用于基于所述无人履带车辆的运动学模型生成轨迹跟踪控制器的目标函数;
约束函数生成模块,用于基于所述无人履带车辆的动力学模型和所述轨迹跟踪控制器的目标函数生成轨迹跟踪控制器的约束函数;
MPC模型生成模块,用于基于所述轨迹跟踪控制器的目标函数和所述轨迹跟踪控制器的约束函数生成MPC模型;
参考轨迹生成模块,用于根据无人履带车辆的位姿特征、道路的曲率特征和实时工况状态信息生成参考轨迹;
参数组合确定模块,用于将所述参考轨迹和车辆轨迹跟踪响应状态信息输入所述MPC模型得到参数组合;所述车辆轨迹跟踪响应状态信息通过实验获得;所述参数组合包括:轨迹曲率特征向量、无人履带车辆姿态参考向量和无人履带车辆速度信息向量;
控制参数确定模块,用于将所述参数组合输入MLP神经网络得到控制参数;
跟踪控制模块,用于根据所述控制参数完成无人履带车辆轨迹的跟踪控制。
此外,还提供了一种无人履带车辆,该无人履带车辆包括有上述轨迹跟踪控制系统。
根据本发明提供的具体实施例,本发明公开了以下技术效果:
本发明提供的无人履带车辆的轨迹跟踪控制方法,以无人履带车辆在越野工况下的大量行驶数据为依托,以轨迹跟踪的高精度、高稳定性、低时间成本为目标,完成最优参数组合的离线提取和最优控制参数的在线预测,为基于模型预测的轨迹跟踪控制器提供时变双参数输入(控制时域长度、时间步长),实现无人履带车辆越野工况下的自适应轨迹跟踪控制,进而有效的提升无人履带车辆轨迹跟踪的工况适应性,实现轨迹跟踪精度、车辆行驶稳定性以及计算时间成本的较好均衡。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明提供的无人履带车辆的轨迹跟踪控制方法的流程图;
图2为本发明实施例提供的实施无人履带车辆的轨迹跟踪控制方法的数据框图;
图3为本发明实施例提供的MLP算法原理图;
图4为本发明提供的无人履带车辆的轨迹跟踪控制系统的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的目的是提供一种无人履带车辆及其轨迹跟踪控制方法和系统,以为无人履带车辆轨迹的跟踪控制提供稳定、可靠的输入,同时降低无人履带车辆轨迹跟踪控制的资源配置,进而能够提升无人履带车辆轨迹跟踪控制性能。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
如图1所示,本发明提供的无人履带车辆的轨迹跟踪控制方法,包括:
步骤100:获取无人履带车辆的运动学模型和无人履带车辆的动力学模型。其中,基于传统的无人履带车辆运动学模型,通过对无人履带车辆接地段的运动学分析,建立了基于瞬时转向中心模型的无人履带车辆的运动学模型。建立得到的基于瞬时转向中心模型的无人履带车辆的运动学模型为:
无人履带车辆的动力学模型为:
其中,x c 为无人履带车辆的重心与几何中心的横向距离,y c 为无人履带车辆的重心与几何中心的纵向距离,F cx 为离心力的横向分量、F cy 为离心力的纵向分量,F r1为一侧履带接地段所受的地面变形阻力,F r2为另一侧履带接地段所受的地面变形阻力,F y1为一侧履带制动力,F y2为另一侧履带驱动力,F x1为一侧履带所受的横向阻力的合力、F x2为另一侧履带所受的横向阻力的合力,m为无人履带车辆的整车质量,δ为无人履带车辆旋转质量换算系数,M μ 为转向阻力矩,M μ =M μ1+M μ2,M μ1为一侧履带所受的转向阻力矩、M μ2为另一侧履带所受的转向阻力矩,J为无人履带车辆绕z’轴的转动惯量,B为履带中心距,为车辆的横向加速度,为车辆的纵向加速度,为车辆的横摆角加速度。
步骤101:基于无人履带车辆的运动学模型生成轨迹跟踪控制器的目标函数。其中,轨迹跟踪控制器的目标函数为:
其中,J(ξ(k),U k )为轨迹跟踪控制器的目标函数,ξ(k)为k时刻的车辆状态向量,U k 为控制序列,N p 为预测时域的个数、N c 为控制时域的个数,且N c ≤N p ,i为个数变量,l e (k+i)为k+i时刻的轨迹跟踪的横向偏差,h e (k+i)为k+i时刻的轨迹跟踪的航向偏差,v(k+i)为k+i时刻的实际车速,l desired (k+i)为k+i时刻的期望横向偏差,h desired (k+i)为k+i时刻的期望航向偏差,v desired (k+i)为k+i时刻的期望车速,Q、R、P、M和N均为权重系数矩阵,v j1(*)为*时刻的一侧履带的卷绕速度,v j2(*)为*时刻的另一侧履带的卷绕速度。
步骤102:基于无人履带车辆的动力学模型和轨迹跟踪控制器的目标函数生成轨迹跟踪控制器的约束函数。轨迹跟踪控制器的约束函数包括:车辆纵向速度约束、转向中心纵向偏移量约束、电机转速约束和电机功率约束。其中各约束条件的具体设置过程在下述实施例中予以说明。
步骤103:基于轨迹跟踪控制器的目标函数和轨迹跟踪控制器的约束函数生成MPC模型。
步骤104:根据无人履带车辆的位姿特征、道路的曲率特征和实时工况状态信息生成参考轨迹。
步骤105:将参考轨迹和车辆轨迹跟踪响应状态信息输入MPC模型得到参数组合。车辆轨迹跟踪响应状态信息通过实验获得。参数组合包括:轨迹曲率特征向量、无人履带车辆姿态参考向量和无人履带车辆速度信息向量。
步骤106:将参数组合输入MLP神经网络得到控制参数。
步骤107:根据控制参数完成无人履带车辆轨迹的跟踪控制。
下面基于PSO-MLP算法,以控制时域长度t l 和控制时间步长t s 为优化目标,以离线和在线相结合的方式,对上述提供的无人履带车辆的轨迹跟踪控制方法进行详细说明。离线部分,首先对不同控制参数下的履带车辆轨迹跟踪响应状态信息进行提取,得到包含运动基元类型、参考轨迹S r 、轨迹跟踪响应信息、车辆状态ξ(k)以及控制参数的车辆历史数据集合DH。然后,基于改进PSO粒子群算法得到不同运动基元下的最优控制参数组合,并构建训练数据集,利用Adam学习率自适应优化算法完成对MLP神经网络模型的训练。在线部分,首先对运动基元类型和车辆状态进行特征匹配,基于MLP神经网络对控制时间步长t s 和控制时域长度t l 进行预测,为轨迹跟踪控制模块提供自适应双参数输入。然后,结合自适应参数输入、车辆模型、约束条件c以及代价函数等建立轨迹跟踪控制器,得到最优控制量u,实现了双侧独立电驱动履带平台的自适应轨迹跟踪控制。基于如图2所示的数据处理框架,该方法的具体实施过程如下:
步骤1:建立履带车辆运动学模型。
首先对履带车辆的相关坐标系进行定义:oxy为固连于地面的空间笛卡尔坐标系,坐标轴x、y所在平面与路面重合;为固连于履带车辆并随履带车辆共同运动的空间笛卡尔坐标系,位于履带车辆的几何中心,为车辆的横向方向,为车辆的纵向方向。为履带车辆的瞬时转向中心,分别为低速侧履带、高速侧履带的瞬时转向中心。
为了给步骤3中模型预测控制器目标函数提供可靠的模型支撑、准确的履带车辆状态空间表达,需要首先建立履带车辆的运动学模型。传统的履带车辆运动学模型可以表示为:
传统的履带车辆运动学模型并没有考虑到履带车辆的滑移滑转特性,而滑移滑转特性是履带车辆在行进过程中不可避免的,因此传统的履带车辆运动学模型精度较差。基于传统的履带车辆运动学模型,通过对履带车辆以及履带接地段的运动学分析,建立了基于瞬时转向中心模型进行履带车辆运动学模型,可以真实的反应履带车辆的运动学特性,为轨迹跟踪控制方法提供可靠、稳定的模型支撑,基于瞬时转向中心的履带车辆运动学微分方程可以表示为:
其中,ξ为无人履带车辆的状态向量,为无人履带车辆的横向速度,为无人履带车辆的纵向速度,为无人履带车辆的横摆角速度,v j1为一侧履带的卷绕速度,v j2为另一侧履带的卷绕速度,φ为履带车辆的几何中心在固连于地面的笛卡尔坐标系下的航向角,为一侧履带的横向方向,为另一侧履带的横向方向。
上述建立的履带车辆运动学模型是后续生成MPC控制器状态空间的基础,该模型的精确性也决定了模型在进行仿真计算时的结果与真实情况的准确性,反映了仿真结果的可靠性。
步骤2:建立履带车辆动力学模型
为了提升履带车辆模型的精度,保证履带车辆在行进过程中的动力性、安全性,需要结合动力学模型对步骤4中轨迹跟踪控制器约束设置施加安全约束。为了简化建模,进行如下所示的假设:
(1)车辆在转向过程中,不考虑履带的沉陷作用和推土作用;
(2)履带载荷沿着履带长度方向均匀分布;
(3)在进行阻力分析时,忽略空气阻力的作用;
(4)在转向过程中的行驶阻力系数与直驶过程相同;
基于达朗贝尔原理,对车辆几何中心建立履带车辆动力学方程组,如下所示:
其中,x c 为无人履带车辆的重心与几何中心的横向距离,y c 为无人履带车辆的重心与几何中心的纵向距离,F cx 为离心力的横向分量、F cy 为离心力的纵向分量,F r1为一侧履带接地段所受的地面变形阻力,F r2为另一侧履带接地段所受的地面变形阻力,F y1为一侧履带制动力,F y2为另一侧履带驱动力,F x1为一侧履带所受的横向阻力的合力、F x2为另一侧履带所受的横向阻力的合力,m为无人履带车辆的整车质量,δ为无人履带车辆旋转质量换算系数,M μ 为转向阻力矩,M μ =M μ1+M μ2,M μ1为一侧履带所受的转向阻力矩、M μ2为另一侧履带所受的转向阻力矩,J为无人履带车辆绕z’轴的转动惯量,B为履带中心距,为车辆的横向加速度,为车辆的纵向加速度,为车辆的横摆角加速度。
步骤3:建立模型预测控制器目标函数
对于步骤1中建立的履带车辆运动学微分方程(公式(2)),使用矩阵代数中有微分方程到状态空间的方法,可以建立该履带车辆的状态空间表达式。另外,在进行代价函数的构建时,需要综合考虑轨迹跟踪精度、稳定性。
为求得在k时刻对于车辆的最优控制序列,需要求使得代价函数取最小值的控制量,从而为步骤6中的最优参数组合数据集的建立提供轨迹跟踪响应状态信息,因此基于离散非线性模型的模型预测控制(MPC)优化问题可以描述为:
优化问题当中的目标函数是进行控制量求解的必须条件,目标函数中各个项目的设置也会直接影响轨迹跟踪的效果以及步骤6中轨迹跟踪效应状态信息的构建。
满足时域约束条件:
其中目标函数为:
其中,ξ(k)为k时刻的车辆状态向量,ξ(k)∈R n ,u(k)为k时刻的控制输入向量,u(k)∈R m ,U k 为控制序列,y c (k)为k时刻的控制输出向量,y b (k)为约束输出向量,y b (k)∈R b ,u min 为控制量下限值,u max 为控制量上限值,Δu min 为控制量变化率的下限值,Δu max 为控制量变化率的上限值,N p 为预测时域的个数,N c 为控制时域的个数,且N c ≤N p ,y min (k+i)为k+i时刻约束输出的下限值,y max (k+i)为k+i时刻约束输出的上限值,l e (k+i)为k+i时刻的轨迹跟踪的横向偏差,h e (k+i)为k+i时刻的轨迹跟踪的航向偏差,v(k+i)为k+i时刻的实际车速,l desired (k+i)为k+i时刻的期望横向偏差,h desired (k+i)为k+i时刻的期望航向偏差,v desired (k+i)为k+i时刻的期望车速,Q、R、P、M和N均为权重系数矩阵,v j1(*)为*时刻的一侧履带的卷绕速度,v j2(*)为*时刻的另一侧履带的卷绕速度。
步骤4:轨迹跟踪控制器约束设置
基于步骤1中的履带车辆运动学模型,进行步骤3中轨迹跟踪控制器的设计,由于并未考虑到履带车辆的动力学特性以及驱动电机的外特性,履带车辆可能会出现侧翻、侧滑、后备功率不足等情况。为了保证车辆的动力性和安全性,需要借助履带车辆的动力学模型以及电机外特性,为轨迹跟踪控制器提供可靠的约束设置。
基于步骤2中建立的履带车辆动力学模型及其约束条件,为保证履带车辆不发生侧翻,车速需要满足:
其中,g为重力加速度,R为履带车辆转向半径,h g 为重心高度,β为车辆的侧倾角度,k r 为侧翻阈值系数。
为了保证车辆在转弯时,不发生失控的情况,车辆的转向中心纵向偏移量需要满足:
其中,L为履带接地段长度,x c 为无人履带车辆的重心与几何中心的横向距离,y c 为无人履带车辆的重心与几何中心的纵向距离,λ c 为转向中心纵向偏移量。
整理得到使履带车辆安全行进的转向中心纵向速度的约束条件:
在规划轨迹已知的条件下,由步骤2中的履带车辆动力学模型,在入弯之前,根据当前弯道对于车辆转向中心纵向速度的约束,对车辆转向的转矩需求进行预测:
其中,T 1、T 2分别为低速侧、高速侧履带转矩需求,η为机械传动效率,i为综合传动比,r为驱动轮半径,T(v y ,ρ)为基于规划轨迹的电机最大转矩需求,ρ为道路的曲率。这些参数表征了履带车辆转弯时需要的驱动力,结合实际电机的外特性曲线,能够对于临界极限转矩工况进行预测,防止极限工况下出现轨迹不可控等危险行为。
将电机的转矩和功率作为约束输出,并设置转矩保护系数k T 和后备功率保护系数k p ,可得:
步骤5:参考轨迹的提取
轨迹规划模块根据车辆位姿特征、道路的曲率特征,基于最优控制理论进行履带车辆运动基元的生成,并结合感知模块提供的实时工况状态信息进行运动基元的拼接,得到具备车辆位姿信息以及运动语义信息的参考轨迹,对车辆的期望路径进行预测:
其中,为参考轨迹点在车体坐标下的坐标,为参考轨迹的车辆三维姿态信息,为参考轨迹的曲率;为车辆的参考速度,为履带车辆运动基元的类型。研究中结合算法验证的实际需求,基于车辆位姿及道路特征,将车辆运动基元分为类直线运动基元、类小曲率转向运动基元、类单移线运动基元、类U型弯运动基元四类。
步骤6:基于改进PSO算法的最优参数组合提取
基于步骤5中提取出来的参考轨迹以及基于实车试验得到的车辆轨迹跟踪响应状态信息将输入到MPC模型(步骤3、步骤4)的控制时域长度、时间步长作为粒子。
适应度函数即目标函数的表达形式如式下所示:
粒子的最新位置和速度需要利用状态更新方程进行迭代求解,其状态更新方程如下式所示:
其中,t、t+1为迭代次数,、分别为第i个粒子在第t、t+1轮迭代的速度,分别为第i个粒子在第t、t+1轮迭代的位置;为第i个粒子累计到第t轮迭代的局部最优值,为粒子群内所有粒子累计到第t轮迭代的全局最优值,c 1、c 2为学习因子,分别代表认知系数和社会系数,通常c 1=c 2∈(0,4),r 1、r 2为(0,1)之间的随机数。为第t轮迭代的惯性因子,w>0,为了解决算法早熟以及迭代后期易在最优解附近震荡的情况,这里采用线性可变权重系数的方法进行权重系数的设置:
其中,w max 、w min 分别为权重系数的最大值和最小值,t max 为最大迭代次数。
步骤7:基于MLP的最优控制参数预测
基于步骤6中提取的最优参数组合,进行MLP神经网络的离线训练和最优控制参数在线预测。本模块的工作原理如图3所示。
MLP神经网络具体算法流程如下所示:
(3)反向传播计算。
未达到迭代终止条件时,需要根据误差的梯度方向对不同网络深度的权重系数矩阵和偏置参数进行更新,然后再重复步骤(2),直至达到迭代收敛条件。
本发明采用Adam学习率自适应优化算法进行MLP神经网络的训练,具体流程如表1所示。
基于上述描述,相对于现有技术,本发明提供的无人履带车辆的轨迹跟踪控制方法具有以下优点:
1、本发明可以保证越野工况下履带车辆轨迹跟踪的行驶安全性,通过建立的运动学模型和安全约束,考虑了车辆侧翻、侧滑失控、电机外特性边界等执行器约束,可以保证履带车辆工作在安全区间。
2、本发明可以提升越野工况下履带车辆轨迹跟踪的精度、稳定性、计算时间成本等综合性能。其中,通过将轨迹跟踪的精度、控制量变化率作为惩罚项进行控制量的求解,保证较高的跟踪精度和较为平滑的控制量变化。通过以跟踪精度、行驶稳定性和计算时间成本为目标,进行最优控制参数的求解,保证车辆轨迹跟踪性能的均衡,提升越野工况下履带车辆轨迹跟踪控制的工况适应性。
对应于上述提供的无人履带车辆的轨迹跟踪控制方法,本发明还提供了一种无人履带车辆的轨迹跟踪控制系统,如图3所示,该轨迹跟踪控制系统包括:
模型获取模块400,用于获取无人履带车辆的运动学模型和无人履带车辆的动力学模型。
目标函数生成模块401,用于基于无人履带车辆的运动学模型生成轨迹跟踪控制器的目标函数。
约束函数生成模块402,用于基于无人履带车辆的动力学模型和轨迹跟踪控制器的目标函数生成轨迹跟踪控制器的约束函数。
MPC模型生成模块403,用于基于轨迹跟踪控制器的目标函数和轨迹跟踪控制器的约束函数生成MPC模型。
参考轨迹生成模块404,用于根据无人履带车辆的位姿特征、道路的曲率特征和实时工况状态信息生成参考轨迹。
参数组合确定模块405,用于将参考轨迹和车辆轨迹跟踪响应状态信息输入MPC模型得到参数组合。车辆轨迹跟踪响应状态信息通过实验获得。参数组合包括:轨迹曲率特征向量、无人履带车辆姿态参考向量和无人履带车辆速度信息向量。
控制参数确定模块406,用于将参数组合输入MLP神经网络得到控制参数。
跟踪控制模块407,用于根据控制参数完成无人履带车辆轨迹的跟踪控制。
此外,还提供了一种无人履带车辆,该无人履带车辆包括有上述轨迹跟踪控制系统。
本说明书中各个实施例采用递进的方式描述,每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
本发明中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。
Claims (3)
1.一种无人履带车辆的轨迹跟踪控制方法,其特征在于,包括:
步骤100:获取无人履带车辆的运动学模型和无人履带车辆的动力学模型;其中,基于传统的无人履带车辆运动学模型,通过对无人履带车辆接地段的运动学分析,建立了基于瞬时转向中心模型的无人履带车辆的运动学模型;
所述基于瞬时转向中心模型的无人履带车辆的运动学模型为:
其中,ξ为无人履带车辆的状态向量,为无人履带车辆的横向速度,为无人履带车辆的纵向速度,为无人履带车辆的横摆角速度,v j1为一侧履带的卷绕速度,v j2为另一侧履带的卷绕速度,φ为履带车辆的几何中心在固连于地面的笛卡尔坐标系下的航向角,为一侧履带的横向方向,为另一侧履带的横向方向;
所述无人履带车辆的动力学模型为:
其中,x c 为无人履带车辆的重心与几何中心的横向距离,y c 为无人履带车辆的重心与几何中心的纵向距离,F cx 为离心力的横向分量、F cy 为离心力的纵向分量,F r1为一侧履带接地段所受的地面变形阻力,F r2为另一侧履带接地段所受的地面变形阻力,F y1为一侧履带制动力,F y2为另一侧履带驱动力,F x1为一侧履带所受的横向阻力的合力、F x2为另一侧履带所受的横向阻力的合力,m为无人履带车辆的整车质量,δ为无人履带车辆旋转质量换算系数,M μ 为转向阻力矩,M μ =M μ1+M μ2,M μ1为一侧履带所受的转向阻力矩、M μ2为另一侧履带所受的转向阻力矩,J为无人履带车辆绕z’轴的转动惯量,B为履带中心距,为车辆的横向加速度,为车辆的纵向加速度,为车辆的横摆角加速度;
步骤101:基于所述无人履带车辆的运动学模型生成轨迹跟踪控制器的目标函数;所述轨迹跟踪控制器的目标函数为:
其中,J(ξ(k),U k )为轨迹跟踪控制器的目标函数,ξ(k)为k时刻的车辆状态向量,U k 为控制序列,N p 为预测时域的个数、N c 为控制时域的个数,且N c ≤N p ,i为个数变量,l e (k+i)为k+i时刻的轨迹跟踪的横向偏差,h e (k+i)为k+i时刻的轨迹跟踪的航向偏差,v(k+i)为k+i时刻的实际车速,l desired (k+i)为k+i时刻的期望横向偏差,h desired (k+i)为k+i时刻的期望航向偏差,v desired (k+i)为k+i时刻的期望车速,Q、R、P、M和N均为权重系数矩阵,v j1(*)为*时刻的一侧履带的卷绕速度,v j2(*)为*时刻的另一侧履带的卷绕速度;
步骤102:基于所述无人履带车辆的动力学模型和所述轨迹跟踪控制器的目标函数生成轨迹跟踪控制器的约束函数;所述轨迹跟踪控制器的约束函数包括:车辆纵向速度约束、转向中心纵向偏移量约束、电机转速约束和电机功率约束;
步骤103:基于所述轨迹跟踪控制器的目标函数和所述轨迹跟踪控制器的约束函数生成MPC模型;
步骤104:根据无人履带车辆的位姿特征、道路的曲率特征和实时工况状态信息生成参考轨迹;
步骤105:将所述参考轨迹和车辆轨迹跟踪响应状态信息输入所述MPC模型得到参数组合;所述车辆轨迹跟踪响应状态信息通过实验获得;所述参数组合包括:轨迹曲率特征向量、无人履带车辆姿态参考向量和无人履带车辆速度信息向量;
步骤106:将所述参数组合输入MLP神经网络得到控制参数;
步骤107:根据所述控制参数完成无人履带车辆轨迹的跟踪控制。
2.一种无人履带车辆的轨迹跟踪控制系统,其特征在于,包括:
模型获取模块,用于获取无人履带车辆的运动学模型和无人履带车辆的动力学模型;其中,基于传统的无人履带车辆运动学模型,通过对无人履带车辆接地段的运动学分析,建立了基于瞬时转向中心模型的无人履带车辆的运动学模型;
所述基于瞬时转向中心模型的无人履带车辆的运动学模型为:
其中,ξ为无人履带车辆的状态向量,为无人履带车辆的横向速度,为无人履带车辆的纵向速度,为无人履带车辆的横摆角速度,v j1为一侧履带的卷绕速度,v j2为另一侧履带的卷绕速度,φ为履带车辆的几何中心在固连于地面的笛卡尔坐标系下的航向角,为一侧履带的横向方向,为另一侧履带的横向方向;
所述无人履带车辆的动力学模型为:
其中,x c 为无人履带车辆的重心与几何中心的横向距离,y c 为无人履带车辆的重心与几何中心的纵向距离,F cx 为离心力的横向分量、F cy 为离心力的纵向分量,F r1为一侧履带接地段所受的地面变形阻力,F r2为另一侧履带接地段所受的地面变形阻力,F y1为一侧履带制动力,F y2为另一侧履带驱动力,F x1为一侧履带所受的横向阻力的合力、F x2为另一侧履带所受的横向阻力的合力,m为无人履带车辆的整车质量,δ为无人履带车辆旋转质量换算系数,M μ 为转向阻力矩,M μ =M μ1+M μ2,M μ1为一侧履带所受的转向阻力矩、M μ2为另一侧履带所受的转向阻力矩,J为无人履带车辆绕z’轴的转动惯量,B为履带中心距,为车辆的横向加速度,为车辆的纵向加速度,为车辆的横摆角加速度;
目标函数生成模块,用于基于所述无人履带车辆的运动学模型生成轨迹跟踪控制器的目标函数;所述轨迹跟踪控制器的目标函数为:
其中,J(ξ(k),U k )为轨迹跟踪控制器的目标函数,ξ(k)为k时刻的车辆状态向量,U k 为控制序列,N p 为预测时域的个数、N c 为控制时域的个数,且N c ≤N p ,i为个数变量,l e (k+i)为k+i时刻的轨迹跟踪的横向偏差,h e (k+i)为k+i时刻的轨迹跟踪的航向偏差,v(k+i)为k+i时刻的实际车速,l desired (k+i)为k+i时刻的期望横向偏差,h desired (k+i)为k+i时刻的期望航向偏差,v desired (k+i)为k+i时刻的期望车速,Q、R、P、M和N均为权重系数矩阵,v j1(*)为*时刻的一侧履带的卷绕速度,v j2(*)为*时刻的另一侧履带的卷绕速度;
约束函数生成模块,用于基于所述无人履带车辆的动力学模型和所述轨迹跟踪控制器的目标函数生成轨迹跟踪控制器的约束函数;所述轨迹跟踪控制器的约束函数包括:车辆纵向速度约束、转向中心纵向偏移量约束、电机转速约束和电机功率约束;
MPC模型生成模块,用于基于所述轨迹跟踪控制器的目标函数和所述轨迹跟踪控制器的约束函数生成MPC模型;
参考轨迹生成模块,用于根据无人履带车辆的位姿特征、道路的曲率特征和实时工况状态信息生成参考轨迹;
参数组合确定模块,用于将所述参考轨迹和车辆轨迹跟踪响应状态信息输入所述MPC模型得到参数组合;所述车辆轨迹跟踪响应状态信息通过实验获得;所述参数组合包括:轨迹曲率特征向量、无人履带车辆姿态参考向量和无人履带车辆速度信息向量;
控制参数确定模块,用于将所述参数组合输入MLP神经网络得到控制参数;
跟踪控制模块,用于根据所述控制参数完成无人履带车辆轨迹的跟踪控制。
3.一种无人履带车辆,其特征在于,包括如权利要求2所述的轨迹跟踪控制系统。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210352604.7A CN114510063B (zh) | 2022-04-06 | 2022-04-06 | 一种无人履带车辆及其轨迹跟踪控制方法和系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210352604.7A CN114510063B (zh) | 2022-04-06 | 2022-04-06 | 一种无人履带车辆及其轨迹跟踪控制方法和系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114510063A CN114510063A (zh) | 2022-05-17 |
CN114510063B true CN114510063B (zh) | 2022-07-15 |
Family
ID=81555032
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210352604.7A Active CN114510063B (zh) | 2022-04-06 | 2022-04-06 | 一种无人履带车辆及其轨迹跟踪控制方法和系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114510063B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115342817B (zh) * | 2022-10-20 | 2023-02-03 | 北京理工大学 | 无人履带车辆轨迹跟踪状态监测方法、系统、设备、介质 |
CN115933643A (zh) * | 2022-11-21 | 2023-04-07 | 深圳市智绘科技有限公司 | 机器人的控制方法、装置、电子设备及存储介质 |
CN117068185B (zh) * | 2023-10-18 | 2024-01-02 | 中汽研(天津)汽车工程研究院有限公司 | 履带车辆轨迹预测方法、设备和介质 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105739305A (zh) * | 2016-01-29 | 2016-07-06 | 北京理工大学 | 一种履带车辆的控制方法 |
CN109799814A (zh) * | 2018-12-29 | 2019-05-24 | 中国北方车辆研究所 | 基于运动学模型的履带无人车辆轨迹跟踪控制器设计方法 |
CN111338354A (zh) * | 2020-03-27 | 2020-06-26 | 淮安信息职业技术学院 | 履带式车辆轨迹跟踪控制方法、装置及系统 |
CN111703432A (zh) * | 2020-06-28 | 2020-09-25 | 湖南大学 | 一种智能履带车辆滑动参数实时估计方法 |
CN111812974A (zh) * | 2020-05-28 | 2020-10-23 | 北京理工大学 | 一种双侧电机驱动履带车辆的综合控制方法 |
CN113467470A (zh) * | 2021-07-23 | 2021-10-01 | 西安理工大学 | 一种无人自主小车的轨迹跟踪控制方法 |
-
2022
- 2022-04-06 CN CN202210352604.7A patent/CN114510063B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105739305A (zh) * | 2016-01-29 | 2016-07-06 | 北京理工大学 | 一种履带车辆的控制方法 |
CN109799814A (zh) * | 2018-12-29 | 2019-05-24 | 中国北方车辆研究所 | 基于运动学模型的履带无人车辆轨迹跟踪控制器设计方法 |
CN111338354A (zh) * | 2020-03-27 | 2020-06-26 | 淮安信息职业技术学院 | 履带式车辆轨迹跟踪控制方法、装置及系统 |
CN111812974A (zh) * | 2020-05-28 | 2020-10-23 | 北京理工大学 | 一种双侧电机驱动履带车辆的综合控制方法 |
CN111703432A (zh) * | 2020-06-28 | 2020-09-25 | 湖南大学 | 一种智能履带车辆滑动参数实时估计方法 |
CN113467470A (zh) * | 2021-07-23 | 2021-10-01 | 西安理工大学 | 一种无人自主小车的轨迹跟踪控制方法 |
Non-Patent Citations (3)
Title |
---|
Research on Lateral-Longitudinal Coupling Trajectory Tracking Control Method for Bilateral Electric Drive Tracked Vehicle;Jiaxing Lu等;《2021 China Automation Congress (CAC)》;20220314;第3235-3240页 * |
井下掘进机行进纠偏调度规划与控制研究;瞿圆媛等;《矿业科学学报》;20200430;第5卷(第02期);第194-202页 * |
基于模型预测控制的履带式无人平台轨迹跟踪;杨福威等;《舰船电子工程》;20180331;第38卷(第3期);第44-50页 * |
Also Published As
Publication number | Publication date |
---|---|
CN114510063A (zh) | 2022-05-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN114510063B (zh) | 一种无人履带车辆及其轨迹跟踪控制方法和系统 | |
CN111413966B (zh) | 一种递进式模型预测无人驾驶规划跟踪协同控制方法 | |
CN110377039B (zh) | 一种车辆避障轨迹规划与跟踪控制方法 | |
CN107561942A (zh) | 基于模型补偿的智能车辆轨迹跟踪模型预测控制方法 | |
Su et al. | Attitude control of underwater glider combined reinforcement learning with active disturbance rejection control | |
CN113703319B (zh) | 基于强化学习的关节模组不等式约束最优鲁棒控制方法 | |
CN113433827B (zh) | 一种质心变化喷砂除锈并联机器人移动平台轨迹跟踪控制方法 | |
CN111812974B (zh) | 一种双侧电机驱动履带车辆的综合控制方法 | |
Kebbati et al. | Lateral control for autonomous wheeled vehicles: A technical review | |
CN115303289A (zh) | 一种基于深度高斯车辆动力学模型、训练方法、智能汽车轨迹跟踪控制方法及终端设备 | |
CN114684199A (zh) | 一种机理分析-数据驱动的车辆动力学串联混合模型、智能汽车轨迹跟踪控制方法及控制器 | |
Yue et al. | Path tracking control of skid-steered mobile robot on the slope based on fuzzy system and model predictive control | |
Yu et al. | Nonlinear model predictive planning and control for high-speed autonomous vehicles on 3D terrains | |
Xiu et al. | Path-following based on fuzzy line-of-sight guidance for a bionic snake robot with unknowns | |
CN114228690A (zh) | 一种基于ddpg和迭代控制的自动驾驶车辆侧倾控制方法 | |
Kabtoul et al. | Proactive and smooth maneuvering for navigation around pedestrians | |
Jiang et al. | Learning based predictive error estimation and compensator design for autonomous vehicle path tracking | |
Fukushima et al. | An optimal control method based on the energy flow equation | |
Gao et al. | Receding horizon tracking control for wheeled mobile robots with time-delay | |
Wohner et al. | Hierarchical nonlinear model predictive control for an autonomous racecar | |
CN114179818A (zh) | 基于自适应预瞄时间和滑模控制的智能汽车横向控制方法 | |
Regolin et al. | A path tracking approach for autonomous driving on slippery surfaces | |
CN111857112A (zh) | 一种汽车局部路径规划方法及电子设备 | |
CN112882460A (zh) | 基于模型预测控制的无人驾驶车辆动态避障轨迹规划方法 | |
CN113050666B (zh) | 水下自主航行器深度和纵向解耦控制方法及系统 |
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 |