CN108594815B - 一种分阶段的轮式机器人移动路径规划方法 - Google Patents
一种分阶段的轮式机器人移动路径规划方法 Download PDFInfo
- Publication number
- CN108594815B CN108594815B CN201810361638.6A CN201810361638A CN108594815B CN 108594815 B CN108594815 B CN 108594815B CN 201810361638 A CN201810361638 A CN 201810361638A CN 108594815 B CN108594815 B CN 108594815B
- Authority
- CN
- China
- Prior art keywords
- speed
- wheeled robot
- target
- current
- tangent point
- 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
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
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Aviation & Aerospace Engineering (AREA)
- Radar, Positioning & Navigation (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Manipulator (AREA)
- Control Of Position, Course, Altitude, Or Attitude Of Moving Bodies (AREA)
Abstract
本发明公开一种分阶段的轮式机器人移动路径规划方法,将轮式机器人从一个已知的初始全状态转变到一个已知的目标全状态的这个过程划分为三个阶段,进行规划控制,速度压制阶段在轮式机器人当前速度方向与当前目标位置方向不一致时,压制速度并调整方向,以尽快朝着目标方向移动;直线段阶段中,直线路径与小圆弧路径相切,以尽可能直线向着目标切点接近,同时控制到达的速度与姿态;小圆弧段阶段中利用最大离心加速度限制圆弧路径的半径,使得轮式机器人沿着此圆弧形路径接近目标并且以期望的角度切入目标。本发明不仅能够实现轮式机器人全状态的改变,并且对周遭环境没有空间大小方面的要求,具有较强的普适性和稳定性。
Description
技术领域
本发明属于机器人运动控制技术领域,具体涉及一种分阶段的轮式机器人移动路径规划方法。
背景技术
随着工业自动化的不断发展,诸如智能小车、移动机器人等轮式机器人在日常生活中的应用也越来越广泛,引导运动主体的移动,已经不仅仅是位置的移动,更包括姿态和速度的改变。高效且稳定的移动方式是后续上层工作的基础,如何能够便捷并且流畅的完成移动过程一直以来都是争相探讨的话题。
为了解决初、末全状态改变中,末速度方向与初、末位置矢量方向不一致的问题,通常情况下是通过一段从起始点到目标点,并且与始末速度方向相切的大圆弧路径完成的。这种方式只有一个过程,相对简单,但是,沿着大圆弧段路径行进,所需的空间太大,对于环境空间有一定的要求,条形空间或较窄的环境下无法使用,正是由于这一弊端,迫切的需要一种更具有普适性的移动方法,本文就是基于这一点,提出了新的移动路径规划方法,可以很好的解决以上问题。
另一方面,本方法能够同步改变姿态、速度、位置,与目前广为采用的将姿态转变、位置转变、速度转变独立依次完成的方法相区别,本方法更加流畅快速。
其次,移动系统的部件控制延迟会导致行进路径产生偏离。如,当轮式机器人本身就具有较大的速度,设置轮式机器人进行转向并且从初始位置到达设定值的这一过程是需要时间的,在这个过程中路径将产生一个较大的偏离,如果没有得到很好的处理,当偏差较大时,可能导致无法到达预设目标状态,甚至轮式机器人失控。本文就这一问题,提出了速度压制方案,可以有效的削弱转向渐变产生的偏离。
发明内容
为了解决现有技术存在的问题,在已知独立控制车轮方向和速度大小的轮式机器人可以实现位置、速度、姿态的同步改变,并且车轮不会侧向横移的前提下,本发明提供一种分阶段的轮式机器人移动路径规划方法,核心技术就是创新性的采用直线+小圆弧的移动路径,从而能够快速平稳的引导移动主体朝着预定的目的进行位置、速度、姿态改变,同时综合考虑了系统部件延迟,提高了各种情景下的适应性、稳定性。
本发明所采用的技术方案提供一种分阶段的轮式机器人移动路径规划方法,将轮式机器人从一个已知的初始全状态转变到一个已知的目标全状态的这个过程划分为以下三个阶段,进行规划控制,所述轮式机器人的全状态包括轮式机器人的位置、速度和姿态,
1)速度压制阶段,该阶段在轮式机器人当前速度方向与当前目标位置方向不一致时,压制速度并调整方向,以尽快朝着目标方向移动;
2)直线段阶段,该阶段直线路径与小圆弧路径相切,以尽可能直线向着目标切点接近,同时控制到达的速度与姿态;
3)小圆弧段阶段,该阶段利用最大离心加速度限制圆弧路径的半径,使得轮式机器人沿着此圆弧形路径接近目标并且以期望的角度切入目标。
而且,不断修正进行移动路径规划,包括每历元根据当前时刻位置、速度和姿态,判断执行对应阶段的处理,直到最终到达目标状态。
而且,不断修正进行移动路径规划,采用以下流程实现,
步骤S1,输入主体的起始状态以及目标状态,分别包括位置、速度、姿态三类信息;
步骤S2,根据目标速度以及轮式机器人最大离心加速度的限制得到小圆弧半径的最小值,并根据轮式机器人当前的状态信息得到切点目标信息;
步骤S3,若此时上一历元的状态为空,直接进入S3-1,若不为空,则通过相邻两个历元的状态信息对当前时刻轮式机器人所处的阶段进行判断,并做出相应处理;
步骤S4,记录轮式机器人当前历元的状态信息,与上一历元状态信息联合,判断目的地处是否处于相邻两历元中某位置,若是,则已经到达终点;否则,将当前历元状态信息传递到后续过程,并从步骤S2开始重新进行下一历元的处理。
而且,判断阶段的实现如下,
1)当相邻两历元位置变化方向与切点目标位置方向相差大于预设的许可角度偏差阈值,且与目的地的距离大于小圆弧段最小半径时,标识第一阶段;
2)当相邻两历元位置变化方向与切点目标位置方向相差小于等于预设的许可角度偏差阈值,且与目的地的距离大于小圆弧段最小半径时,标识第二阶段;
3)切点目标位置位于相邻两历元某中间位置时,标识第三阶段。
本发明的优点为:
(1)此方法适用于所有轮式机器人,是一种通用性较强的方法。
(2)本发明中直线+小圆弧的方式,不仅能够很好的实现轮式机器人全状态的改变,并且对周遭环境不存在空间大小方面的要求,解决了常规几何关系解算中大圆弧模型无法适应较小空间的情况,该模型适用于各种情景,具有较强的普适性。
(3)考虑到了部件控制延迟,在应用中,当移动主体给定初始速度,由于设置转向角度是逐渐变化的过程,因此无法达到预期的直线效果,而会形成一段弧线路径,采用对速度进行压制的方式,削弱转角渐变产生的弧线效应,使系统更加可控。
(4)考虑到判断到达目的与否条件的合理性,采用了相邻时间节点对比的方式,相较于阈值判断,具有更好的稳定性。
(5)综合考虑了移动主体的离心加速度,最大加速度以及最大速度等物理特性,最大限度的保证了轮式机器人快速稳定的运动,更加高效平稳,在实际应用中具有很大的优势。
(6)每历元都对状态进行重新计算、判断,增加了移动系统的容错性、可靠性。
附图说明
图1为本发明实施例的全局坐标系,局部坐标系下轮式机器人姿态、速度示意图。
图2为本发明实施例的分阶段移动路径规划路径示意图。
图3为本发明实施例的计算全局坐标系下速度方向几何关系示意图。
图4为本发明实施例的确定圆弧圆心的几何关系示意图。
图5为本发明实施例的确定切点位置的几何关系示意图。
图6为本发明实施例的移动路径规划方案流程示意图。
具体实施方式
以下结合附图和实施例详细说明本发明的技术方案。
本发明将轮式机器人的位置、速度、姿态三者合称为轮式机器人的全状态,此轮式机器人可以同时独立调控速度大小和方向,本发明将轮式机器人从一个已知的初始全状态转变到一个已知的目标全状态的这个过程划分为三个阶段。为了解决初、末全状态改变中,末速度方向与初、末位置矢量方向不一致的问题,区别于传统做法中需要较宽的空间进行速度方向的切入,采用本方法可以实现在较为狭长的空间环境下任意全状态改变;另一方面,本方法能够同步改变姿态、速度、位置,与目前广为采用的将姿态转变、位置转变、速度转变独立依次完成的方法相区别,本方法更加流畅快速。
三个阶段如下:
1)速度压制阶段:考虑到行进速度与目标方向不一致的情况,压制速度并调整方向,以尽快朝着目标方向移动
2)直线段:直线路径与小圆弧路径相切,以尽可能直线向着目标切点接近,同时控制到达的速度与姿态
3)小圆弧段:利用最大离心加速度限制圆弧路径的半径,使得轮式机器人沿着此圆弧形路径接近目标并且以期望的角度切入目标
第一阶段速度压制阶段,通过计算出轮式机器人当前速度方向与当前目标位置方向之间的夹角,记为θ,当θ较大时,以最大加速度进行减速,从而削弱当速度方向与目标方向不一致时,转向渐变引发的垂直于路径方向的偏移。
第二阶段直线路径,为了保证移动过程的平稳性,直线路径与第三阶段中的小圆弧路径尽量相切,直线路径的目标点为切点,可以通过一定的几何关系解算得到。另一方面,为了最大程度上优化轮式机器人的功能,在满足最大加速度以及最大速度约束的前提下,先以最大加速度加速到最大速度再以负向的最大加速度减速到目标速度。
第三阶段小圆弧路径,此阶段的目标即为最终目的地,轮式机器人的速度大小保持目标速度大小恒定不变,在满足最大离心加速度的条件下,找到最小半径,尽量减小整个路径占用的空间。
整体上看,轮式机器人的移动路径规划是将整个过程分为了三个阶段,但实际上移动路径规划是在不断修正进行的,每历元根据当前时刻位置,速度,姿态,判断具体执行第几个阶段的命令,直到最终到达目标状态。
因此本发明提出,对移动路径规划的不断修正以及每历元进行阶段判断,实现方式如下:1)当相邻两历元位置变化方向与切点目标位置方向相差较大,且与目的地的距离大于小圆弧段最小半径时,标识第一阶段。
2)当相邻两历元位置变化方向与切点目标位置方向相差较小,且与目的地的距离大于小圆弧段最小半径时,标识第二阶段。
3)切点目标位置位于相邻两历元某中间位置时,标识第三阶段。
如图1所示,实施例定义了此运动模型涉及的两个坐标系:
1)局部坐标系
以轮式机器人的中心P为坐标原点,轮式机器人的前方为Y轴正方向,右方为X轴正方向,在轮式机器人的运动过程中,该坐标系随着轮式机器人的移动而移动。用(X,Y)Local标识
2)全局坐标系
以空间中一点Q为坐标原点,水平方向为X轴,垂直方向为Y轴,在轮式机器人的运动过程中,该坐标系不发生平移和旋转。用(X,Y)Global标识
同时,定义了轮式机器人的速度和姿态
1)轮式机器人的速度
在局部坐标系下,表示为(v,α)。α∈(-π,π],为轮式机器人速度矢量与局部坐标系X轴正方向的夹角,一、二象限为正,三、四象限为负,v是速度大小。
2)轮式机器人的姿态
在全局坐标系下,表示为β。β∈(-π,π],移动主体正方向与全局坐标系X轴正方向的夹角,一、二象限为正,三、四象限为负。
如图2所示,实施例中整体移动路径规划的路径如下:
1)A点为轮式机器人的瞬时位置,其质心坐标表示为:A:(XA,YA)Global,速度表示为(VA,αA)Local,姿态表示为(βA)Global
2)B点为目的地位置,其坐标表示为:B:(XB,YB)Global,速度表示为(VB,αB)Local,姿态表示为(βB)Global
3)C点为切点目标位置,其坐标表示为:C:(XC,YC)Global,速度表示为(VC,αC)Local,姿态表示为(βC)Global
在速度压制部分,通过最大加速度进行减速只能使得垂直方向的偏移削弱,无法消除。
在直线路径过程中,包含着一个加减速过程,加减速阶段的判断通过当前速度,目标速度,当前轮式机器人与目标位置之间的距离,以及最大加速度来进行判断。
对于轮式机器人来说,当接近于切点C时,轮式机器人的姿态、速度大小都已经达到目的地的要求,小圆弧阶段保持固定的姿态、固定的速度大小,只调整速度的方向。
实施例中,移动路径规划的具体步骤如下:
1)输入初始位置、速度、姿态,以及目标位置、速度、姿态
本步骤确定运动主体的起始状态以及目标状态,均包括位置、速度、姿态三类信息。
为更加周全细致的考虑各种情况的适用性,对运动过程进行初步预估,并对特殊情况进行处理。
a)根据始、末位置连线的垂向空间大小来限制小圆弧段最大半径RMax,并且由最大离心加速度ac的约束,通过
计算得到当前环境下允许的小圆弧段最大速度VB_Max,若VB>VB_Max,则令VB=VB_Max。此时由于空间受限,小车无法达到目标速度。
b)当目标速度为0时,将其设置为一个接近于0的较小正值,避免接近目的地的过程中运动的太慢,长时间无法抵达
2)计算小圆弧半径以及对于小圆弧路径圆心的确定,进而唯一确定切点
本步骤根据目标速度以及轮式机器人最大离心加速度的限制得到小圆弧半径的最小值,并根据轮式机器人当前的状态信息在一定的几何关系基础上计算出切点目标信息。
a)小圆弧半径R
VB是目标速度大小
ac是轮式机器人最大离心加速度,取决于轮式机器人的物理特性
b)确定小圆弧圆心坐标O:(XO,YO)Global
如图3中已知目的地B点的状态信息:B:(XB,YB)Global、(VB,αB)Local、(βB)Local。则全局坐标系下B点速度方向αB_Glabal的计算如下:
继而根据圆弧半径以及目标速度方向可以得到如图4中所示两个大小相同且相切的圆。圆心O1、O2的位置信息如下:
R为a)中计算得到的小圆弧段半径
在二者中就近选择一个圆心:
AngleABX=arctan(YA-YB,XA-XB)
将AngleO1BX归算到(-π,π):
if AngleO1BX>2π,AngleO1BX=AngleO1BX-2π;
if AngleO1BX<-2π,AngleO1BX=AngleO1BX+2π;
AngleABO1=AngleABX-AngleO1BX
将AngleABO1归算到(-π,π):
if AngleABO1>2π,AngleABO1=AngleABO1-2π;
if AngleABO1<-2π,AngleABO1=AngleABO1+2π;
其本质上是判断O1、O2中处于起始点和目标点中间的圆心位置,将其记作
O:(XO,YO)Global
c)计算切点的全状态信息
如图5中所示切点几何关系示意图
切点C的速度大小Vc为:
Vc=VB
切点C的姿态βc为:
βc=βB
下面求得切点C的位置C:(XC,YC)Global
为了简化解算上下切点的坐标,不考虑姿态,建立新坐标系。如图4中所示。新坐标系以轮式机器人质心为原点,轮式机器人质心到圆弧圆心方向为X轴正方向,X轴逆时针旋转90度到达Y轴正方向。用(X,Y)Temp标识
记全局坐标系下圆心O到轮式机器人瞬时质心A的距离为L:
在新坐标系下,圆心点坐标O:(L,0)Temp,圆弧半径为R,上方、下方切点坐标分别为:
新坐标系下坐标逆时针旋转θ后,平移(XA,YA)Global,即可转换到全局坐标系下:
由以下公式可将新坐标系下的坐标转换到全局坐标系中:
新坐标系下目标速度方向为αB_Temp
αB_Temp=αB_Glabal+θ
为了使轮式机器人能够以小半圆弧行进,且达到目标要求的速度方向,切点C的位置判断依据如下:
ifαB_Temp>0,则切点在圆心下方,否则,切点在圆心上方
由此得到唯一的切点C在新坐标系下的坐标C:(XC_Temp,YC_Temp)Temp
因此,全局坐标系下,唯一的切点C的坐标为:
XC=XC_Tempcosθ+YC_Tempsinθ+XA
YC=YC_Tempcosθ-XC_Tempsinθ+YA
3)判断轮式机器人此时所处的阶段并执行相应的指令
若此时上一历元的状态为空,直接进入S3-1,若不为空,则通过相邻两个历元的状态信息对当前时刻轮式机器人所处的阶段进行判断,并做出相应指令。判断条件及执行指令如下:
S3-1当相邻两历元位置变化方向与切点目标位置方向相差较大,且与目的地的距离大与小圆弧段最小半径时,标识第一阶段。
执行动作:通过上下两个历元位置变化计算出轮式机器人当前速度方向与切点目标位置方向之间的夹角,记为θ,当角度不为零时,将cosθ作为权值与轮式机器人运动的最大速度相乘,计算得到当前情况下的理想速度,以最大加速度进行加减速,使轮式机器人达到此理想速度。
S3-2当相邻两历元位置变化方向与切点目标位置方向相差较小,且与目的地的距离大于小圆弧段最小半径时,标识第二阶段。
执行动作:通过最终目标速度和当前时刻轮式机器人速度大小以及当前时刻与切点目标之间的距离对轮式机器人速度大小进行调控更新,实现加减速过程。
S3-3当目的地之间的距离小于等于小圆弧段最小半径时,标识第三阶段。
执行动作:调整轮式机器人转向角度,并保持速度大小不变,走小圆弧段路径。
实施例中,将轮式机器人上一历元位置标记为A1:(XA1,YA1)Global,当前历元位置标记为A2:(XA2,YA2)Global,切点位置为C:(XC,YC)Global,目的地为B:(XB,YB)Global。
全局坐标系下:
轮式机器人当前时刻的运动方向γA12为
切点位置方向γAC为
轮式机器人与切点目标之间的距离DAC为:
相邻两历元轮式机器人到切点目标方向的夹角为:
由于轮式机器人的实际运动状态受到各种复杂情况的影响,难免产生微小的偏差,在判断当前速度方向是否与切点目标位置方向相同时,不能够简单将两者完全相等作为判断条件,需要设立一个较小的许可角度偏差阈值,使得两角度偏差在此误差许可范围内时,则视为两角度相等。
本实施例假设δ为许可角度偏差阈值,为一个较小且接近于0的正数,具体实施时,可以预先设置为经验值。基于许可角度偏差阈值判断得到三个阶段实现如下:
a)若
sin|γA12-γAC|>δ且DAC>R
则表示轮式机器人当前运动方向与切点目标位置方向不同,且与目的地的距离大于小圆弧段最小半径,标识为第一阶段
b)若
sin|γA12-γAC|≤δ且DAC>R
则表示轮式机器人当前运动方向与切点目标位置方向相同,且与目的地的距离大于小圆弧段最小半径,标识第二阶段。
c)若满足
cosγA1_C_A2<0
则表示A1、A2与切点C形成的张角为钝角,也就是A1、A2位于切点C异侧,此时标识第三阶段。
根据轮式机器人此时所处的阶段,执行相应的任务实现如下:
a)若为第一阶段,则进行速度压制,如图6中所示。
当前时刻轮式机器人的运动速度大小为VA,轮式机器人目标速度大小为VB。根据3)中得到的轮式机器人与切点目标之间的距离DAC。计算AC段的最大速度VAC_Max,其中,a为轮式机器人运动的最大加速度,由轮式机器人物理特性决定。由公式:
得到:
若上一历元情况未知,则通过当前历元轮式机器人的姿态βA以及速度方向αA计算轮式机器人当前时刻的运动方向γA12:
之前已得到轮式机器人当前时刻的运动方向γA12和切点位置方向γAC,计算当前时刻轮式机器人运动速度方向与切点目标位置方向的角度差异,得到理想速度的权值K为:
K=cos|γA1-γAC|
此时理想速度VA′为:
VA′=VAC_Max×K
对当前速度VA进行加减速,使其到达理想速度,其中t为相邻两历元时间间隔,a为轮式机器人运动的最大加速度:
ifVA+at<VA′,VA=VA+at
ifVA<VA′<VA+at,VA=VA′
ifVA-at<VA′<VA,VA=VA′
ifVA′<VA-at,VA=VA-at
b)直线路径下的加减速详细判断过程,其中,a为轮式机器人运动的最大加速度,VAC_Max为a)中计算得到的AC段轮式机器人运动的最大速度。
已知当前轮式机器人的运动速度VA,以及目标速度VC的情况下,减速所需要的最小距离S为:
通过判断此时轮式机器人至切点的距离DAC与S的大小以及VA与VAC_Max的大小关系,来确定速度大小的变化:
若DAC<S则以-a减速
若DAC>S且VA<VAC_Max,则以a加速
若DAC>S且VA≥VAC_Max,则将速度大小设置为VAC_Max
同时,边走直线加减速边调整姿态,将轮式机器人质心速度大小和方向分解到各转向控制部件,保持速度方向不改变的前提下,各转向控制部分在原有速度大小基础上进行微小的改变,来实现姿态的调整。
c)小圆弧路径部分保持速度大小不变,调整速度方向
轮式机器人质心速度、姿态保持不变,通过平台物理结构,将质心速度投射到各个转向控制部分,计算得到能够保持轮式机器人质心沿固定弧线行进的轮子速度以及转向角度数据。优选建议的轮子参数计算步骤如下:
1、计算当前时刻小车质心运动的速度方向。已知起始点A,末位置点B,末速度方向,若末速度方向与AB方向平行,则A点时刻小车的速度方向即为末速度方向;否则计算经过A、B点,且与末速度方向相切的圆形路径,A点时刻小车的速度方向与此圆相切。
2、计算当前时刻小车质心角速度大小。根据当前速度以及距离,近似得到运动到目的地所需要的时间。初、末姿态的角度变化大小除上此时长得到当前时刻小车质心运动角度。
3、此时,将A点时刻小车质心的速度方向,以及质心速度大小和质心角速度大小,通过小车轮子与质心之间的位置关系对车轮侧向受力情况进行约束,使车轮满足侧向不受力的条件,确定各个车轮当前的速度大小和转向角度。
4)阶段任务执行完毕后,判断是否抵达最终目的地
本步骤记录轮式机器人当前历元的状态信息,与上一历元状态信息联合,判断目的地是否处于相邻两历元中某位置,若是,则已经到达终点;否则,将当前历元状态信息传递到后续过程,并从步骤2开始,进行下一历元的处理。
相邻两历元轮式机器人到目的地方向的夹角为
若满足:
cosγA1_B_A2<0
则A1、A2位于目的地B异侧,表明轮式机器人到达了目的地;否则从步骤1开始继续进行,直到最终抵达目的地。
具体实施时,可采用计算机软件实现以上流程的自动运行。
应当理解的是,上述针对较佳实施例的描述较为详细,并不能因此而认为是对本发明专利保护范围的限制,本领域的普通技术人员在本发明的启示下,在不脱离本发明权利要求所保护的情况下,还可以做出替换或变形,均落入本发明的保护范围之内,本发明的请求保护范围应以所附权利要求为准。
Claims (5)
1.一种分阶段的轮式机器人移动路径规划方法,其特征在于:用于实现在狭长的空间环境下任意全状态改变,将轮式机器人从一个已知的初始全状态转变到一个已知的目标全状态的这个过程划分为以下三个阶段,进行规划控制,所述轮式机器人的全状态包括轮式机器人的位置、速度和姿态,
1)速度压制阶段,该阶段在轮式机器人当前速度方向与当前目标位置方向不一致时,压制速度并调整方向,以尽快朝着目标方向移动;
2)直线段阶段,该阶段直线路径与小圆弧路径相切,以尽可能直线向着目标切点接近,同时控制到达的速度与姿态;
3)小圆弧段阶段,该阶段利用最大离心加速度限制圆弧路径的半径,使得轮式机器人沿着此圆弧形路径接近目标并且以期望的角度切入目标;
轮式机器人移动时不断修正进行移动路径规划,包括每历元根据当前时刻位置、速度和姿态,判断执行对应阶段的处理,直到最终到达目标状态;
不断修正进行移动路径规划,包括定义坐标系如下,
1)局部坐标系
以轮式机器人的中心P为坐标原点,轮式机器人的前方为Y轴正方向,右方为X轴正方向,在轮式机器人的运动过程中,该坐标系随着轮式机器人的移动而移动,用(X,Y)Local标识;
2)全局坐标系
以空间中一点Q为坐标原点,水平方向为X轴,垂直方向为Y轴,在轮式机器人的运动过程中,该坐标系不发生平移和旋转,用(X,Y)Global标识;
定义整体移动路径规划的路径如下,
1)A点为轮式机器人的瞬时位置,其质心坐标表示为:A:(XA,YA)Global,速度表示为(VA,αA)Local,姿态表示为(βA)Global;
2)B点为目的地位置,其坐标表示为:B:(XB,YB)Global,速度表示为(VB,αB)Local,姿态表示为(βB)Global;
3)C点为切点目标位置,其坐标表示为:C:(XC,YC)Global,速度表示为(VC,αC)Local,姿态表示为(βC)Global;
并采用以下流程实现,
步骤S1,输入主体的起始状态以及目标状态,分别包括位置、速度、姿态三类信息;
对运动过程进行初步预估,并对特殊情况进行处理如下,
a)根据始、末位置连线的垂向空间大小限制小圆弧段最大半径RMax,并且由最大离心加速度ac约束,通过
计算得到当前环境下允许的小圆弧段最大速度VB_Max,若VB>VB_Max,则令VB=VB_Max,此时由于空间受限,小车无法达到目标速度;
b)当目标速度为0时,将其设置为一个接近于0的较小正值,避免接近目的地的过程中运动的太慢,长时间无法抵达;
步骤S2,根据目标速度以及轮式机器人最大离心加速度的限制得到小圆弧半径的最小值,并根据轮式机器人当前的状态信息得到切点目标信息;实现如下,
a)求取小圆弧半径R,
VB是目标速度大小;
ac是轮式机器人最大离心加速度;
b)根据目的地的状态信息确定小圆弧圆心坐标,继而根据圆弧半径以及目标速度方向得到两个大小相同且相切的圆,圆心记为O1、O2,判断O1、O2中处于起始点和目标点中间的圆心位置;
c)计算切点的全状态信息如下,
切点C的速度大小Vc为:
Vc=VB
切点C的姿态βc为:
βc=βB
设切点C的位置为:C:(XC,YC)Global
为了简化解算上下切点的坐标,不考虑姿态,建立新坐标系;新坐标系以轮式机器人质心为原点,轮式机器人质心到圆弧圆心方向为X轴正方向,X轴逆时针旋转90度到达Y轴正方向,用(X,Y)Temp标识;
记全局坐标系下圆心O:(XO,YO)Global,O到轮式机器人瞬时质心A的距离为L:
在新坐标系下,圆心点坐标O:(L,0)Temp,圆弧半径为R,上方、下方切点坐标分别为:
新坐标系下坐标逆时针旋转θ后,平移(XA,YA)Global,转换到全局坐标系下:
由以下公式将新坐标系下的坐标C:(XC_Temp,YC_Temp)Temp转换到全局坐标系中:
新坐标系下目标速度方向为αB_Temp
αB_Temp=αB_Global+θ
为了使轮式机器人能够以小半圆弧行进,且达到目标要求的速度方向,切点C的位置判断依据如下:
if αB_Temp>0,则切点在圆心下方,否则,切点在圆心上方
由此得到唯一的切点C在新坐标系下的坐标C:(XC_Temp,YC_Temp)Temp
全局坐标系下,唯一的切点C的坐标为:
XC=XC_Tempcosθ+YC_Tempsinθ+XA
YC=YC_Tempcosθ-XC_Tempsinθ+YA
步骤S3,若此时上一历元的状态为空,直接进入S3-1,若不为空,则通过相邻两个历元的状态信息对当前时刻轮式机器人所处的阶段进行判断,并做出相应处理;
当处于速度压制阶段,执行动作:通过上下两个历元位置变化计算出轮式机器人当前速度方向与切点目标位置方向之间的夹角,记为θ’,当角度不为零时,将cosθ’作为权值与轮式机器人运动的最大速度相乘,计算得到当前情况下的理想速度,以最大加速度进行加减速,使轮式机器人达到此理想速度;
当处于直线段阶段,执行动作:通过最终目标速度和当前时刻轮式机器人速度大小以及当前时刻与切点目标之间的距离对轮式机器人速度大小进行调控更新,实现加减速过程;
当处于小圆弧段阶段,执行动作:调整轮式机器人转向角度,并保持速度大小不变,走小圆弧段路径;
将轮式机器人上一历元位置标记为A1:(XA1,YA1)Global,当前历元位置标记为A2:(XA2,YA2)Global,切点位置为C:(XC,YC)Global,目的地为B:(XB,YB)Global;
全局坐标系下:
轮式机器人当前时刻的运动方向γA12为
切点位置方向γAC为
轮式机器人与切点目标之间的距离DAC为:
相邻两历元轮式机器人到切点目标方向的夹角为:
步骤S4,记录轮式机器人当前历元的状态信息,与上一历元状态信息联合,判断目的地处是否处于相邻两历元中某位置,若是,则已经到达终点;否则,将当前历元状态信息传递到后续过程,并从步骤S2开始进行下一历元的处理。
2.根据权利要求1所述分阶段的轮式机器人移动路径规划方法,其特征在于:判断阶段的实现如下,
1)当相邻两历元位置变化方向与切点目标位置方向相差大于预设的许可角度偏差阈值,且与目的地的距离大于小圆弧段最小半径时,标识第一阶段,表示速度压制阶段;
2)当相邻两历元位置变化方向与切点目标位置方向相差小于等于预设的许可角度偏差阈值,且与目的地的距离大于小圆弧段最小半径时,标识第二阶段,表示直线段阶段;
3)切点目标位置位于相邻两历元某中间位置时,标识第三阶段,表示小圆弧段阶段。
3.根据权利要求2所述分阶段的轮式机器人移动路径规划方法,其特征在于:设δ为许可角度偏差阈值,为一个较小且接近于0的正数,基于许可角度偏差阈值判断得到三个阶段实现如下:
若
sin|γA12-γAC|>δ且DAC>R
则表示轮式机器人当前运动方向与切点目标位置方向不同,且与目的地的距离大于小圆弧段最小半径,标识为第一阶段;
若
sin|γA12-γAC|≤δ且DAC>R
则表示轮式机器人当前运动方向与切点目标位置方向相同,且与目的地的距离大于小圆弧段最小半径,标识第二阶段;
若满足
cosγA1_C_A2<0
则表示A1、A2与切点C形成的张角为钝角,也就是A1、A2位于切点C异侧,此时标识第三阶段;
根据轮式机器人此时所处的阶段,执行相应的任务实现如下:
a)若为第一阶段,则进行速度压制,当前时刻轮式机器人的运动速度大小为VA,轮式机器人目标速度大小为VB,根据轮式机器人与切点目标之间的距离DAC;计算AC段的最大速度VAC_Max,其中,a为轮式机器人运动的最大加速度,由公式:
得到:
若上一历元情况未知,则通过当前历元轮式机器人的姿态βA以及速度方向αA计算轮式机器人当前时刻的运动方向γA12:
根据轮式机器人当前时刻的运动方向γA12和切点位置方向γAC,计算当前时刻轮式机器人运动速度方向与切点目标位置方向的角度差异,得到理想速度的权值K为:
K=cos|γA12-γAC|
此时理想速度V′A为:
V′A=VAC_Max×K
对当前速度VA进行加减速,使其到达理想速度,其中t为相邻两历元时间间隔,a为轮式机器人运动的最大加速度:
if VA+at<V′A,VA=VA+at
if VA<V′A<VA+at,VA=V′A
if VA-at<V′A<VA,VA=V′A
if V′A<VA-at,VA=VA-at
b)直线路径下的加减速详细判断过程,其中,a为轮式机器人运动的最大加速度,VAC_Max为a)中计算得到的AC段轮式机器人运动的最大速度;
已知当前轮式机器人的运动速度VA,以及目标速度VC的情况下,减速所需要的最小距离S为:
通过判断此时轮式机器人至切点的距离DAC与S的大小以及VA与VAC_Max的大小关系,来确定速度大小的变化:
若DAC<S则以-a减速
若DAC>S且VA<VAC_Max,则以a加速
若DAC>S且VA≥VAC_Max,则将速度大小设置为VAC_Max
同时,边走直线加减速边调整姿态,将轮式机器人质心速度大小和方向分解到各转向控制部件,保持速度方向不改变的前提下,各转向控制部分在原有速度大小基础上进行微小的改变,来实现姿态的调整;
c)小圆弧路径部分保持速度大小不变,调整速度方向。
4.根据权利要求3所述分阶段的轮式机器人移动路径规划方法,其特征在于:
轮式机器人质心速度、姿态保持不变,通过平台物理结构,将质心速度投射到各个转向控制部分,计算得到能够保持轮式机器人质心沿固定弧线行进的轮子速度以及转向角度数据;轮子参数计算步骤如下:
1)计算当前时刻小车质心运动的速度方向,已知起始点A,末位置点B,末速度方向,若末速度方向与AB方向平行,则A点时刻小车的速度方向即为末速度方向;否则计算经过A、B点,且与末速度方向相切的圆形路径,A点时刻小车的速度方向与此圆相切;
2)计算当前时刻小车质心角速度大小,根据当前速度以及距离,近似得到运动到目的地所需要的时间,初、末姿态的角度变化大小除以时长得到当前时刻小车质心运动角度;
3)此时,将A点时刻小车质心的速度方向,以及质心速度大小和质心角速度大小,通过小车轮子与质心之间的位置关系对车轮侧向受力情况进行约束,使车轮满足侧向不受力的条件,确定各个车轮当前的速度大小和转向角度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810361638.6A CN108594815B (zh) | 2018-04-20 | 2018-04-20 | 一种分阶段的轮式机器人移动路径规划方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810361638.6A CN108594815B (zh) | 2018-04-20 | 2018-04-20 | 一种分阶段的轮式机器人移动路径规划方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108594815A CN108594815A (zh) | 2018-09-28 |
CN108594815B true CN108594815B (zh) | 2021-02-02 |
Family
ID=63613523
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810361638.6A Active CN108594815B (zh) | 2018-04-20 | 2018-04-20 | 一种分阶段的轮式机器人移动路径规划方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108594815B (zh) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109375624B (zh) * | 2018-11-12 | 2022-02-15 | 楚天智能机器人(长沙)有限公司 | 一种双舵轮agv圆弧路径生成方法、装置及介质 |
CN109828599B (zh) * | 2019-01-08 | 2020-12-15 | 苏州极目机器人科技有限公司 | 飞行器作业路径规划方法以及控制装置和控制设备 |
CN111714031B (zh) * | 2019-03-20 | 2023-10-03 | 北京奇虎科技有限公司 | 扫地设备执行清扫的方法、装置及电子设备 |
CN110244718B (zh) * | 2019-06-04 | 2022-01-28 | 淮阴工学院 | 一种可自动避障的巡视智能小车 |
CN113126607B (zh) * | 2019-12-31 | 2024-03-29 | 深圳市优必选科技股份有限公司 | 一种机器人及其运动控制方法和装置 |
CN111329399B (zh) * | 2020-04-09 | 2021-09-10 | 湖南格兰博智能科技有限责任公司 | 一种基于有限状态机的扫地机目标点导航方法 |
CN112650234B (zh) * | 2020-12-16 | 2022-05-17 | 浙江大学 | 一种双足机器人的路径规划方法 |
CN112729305B (zh) * | 2020-12-21 | 2022-03-08 | 中国人民解放军军事科学院国防科技创新研究院 | 一种基于单飞行器导引头图像信息的多目标定位方法 |
CN113252040B (zh) * | 2021-05-08 | 2022-10-18 | 云南财经大学 | 一种改进的agv小车二维码弧线导航方法 |
CN113377104A (zh) * | 2021-06-02 | 2021-09-10 | 北京布科思科技有限公司 | 基于差速模型的机器人位置控制方法、装置 |
CN114859941B (zh) * | 2022-07-05 | 2022-09-23 | 上海仙工智能科技有限公司 | 一种在弯曲光滑路径下的速度规划方法、存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105082156A (zh) * | 2015-08-12 | 2015-11-25 | 珞石(北京)科技有限公司 | 一种基于速度最优控制的空间轨迹平滑方法 |
CN105737838A (zh) * | 2016-02-22 | 2016-07-06 | 广东嘉腾机器人自动化有限公司 | 一种agv路径跟踪方法 |
CN107030697A (zh) * | 2017-04-28 | 2017-08-11 | 广州大学 | 一种机器人笛卡尔空间平滑轨迹的规划方法 |
CN107092266A (zh) * | 2017-06-26 | 2017-08-25 | 东北大学 | 一种移动车轨迹跟踪控制方法 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103454978B (zh) * | 2013-08-30 | 2016-03-02 | 成都乐创自动化技术股份有限公司 | 点胶、涂胶设备圆弧过渡轨迹速度规划方法 |
CN105500354B (zh) * | 2016-02-02 | 2017-05-17 | 南京埃斯顿机器人工程有限公司 | 一种工业机器人应用的过渡轨迹规划方法 |
-
2018
- 2018-04-20 CN CN201810361638.6A patent/CN108594815B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105082156A (zh) * | 2015-08-12 | 2015-11-25 | 珞石(北京)科技有限公司 | 一种基于速度最优控制的空间轨迹平滑方法 |
CN105737838A (zh) * | 2016-02-22 | 2016-07-06 | 广东嘉腾机器人自动化有限公司 | 一种agv路径跟踪方法 |
CN107030697A (zh) * | 2017-04-28 | 2017-08-11 | 广州大学 | 一种机器人笛卡尔空间平滑轨迹的规划方法 |
CN107092266A (zh) * | 2017-06-26 | 2017-08-25 | 东北大学 | 一种移动车轨迹跟踪控制方法 |
Non-Patent Citations (3)
Title |
---|
A Method of Path Planning for Tractor-Trailer Mobile Robot Based on the Concept of Global-Width;Qiang Han; Yalou Huang;Jing Yuan;Yewei Kang;《Fifth World Congress on Intelligent Control and Automation (IEEE Cat. No.04EX788)》;20040630;全文 * |
基于圆弧_直线_圆弧理论的移动机器人路径规划;王毅 等;《科学技术与工程》;20180131;第18卷(第2期);第111-117页 * |
搬运机器人笛卡尔空间轨迹规划研究;林仕高;《中国优秀硕士学位论文全文数据库信息科技辑》;20140115(第01期);正文第21-22,44,-45页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108594815A (zh) | 2018-09-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108594815B (zh) | 一种分阶段的轮式机器人移动路径规划方法 | |
CN106647282B (zh) | 一种考虑末端运动误差的六自由度机器人轨迹规划方法 | |
CN107030697B (zh) | 一种机器人笛卡尔空间平滑轨迹的规划方法 | |
CN107615194B (zh) | 加工轨迹平滑转接的方法及加工装置 | |
WO2018126355A1 (zh) | 机器人运动轨迹规划方法及相关装置 | |
WO2017113219A1 (zh) | 折弯跟随轨迹规划方法、装置及系统 | |
CN105082156B (zh) | 一种基于速度最优控制的空间轨迹平滑方法 | |
CN103760908B (zh) | 一种巡视器闭环跟踪控制方法 | |
WO2018126354A1 (zh) | 机器人运动轨迹规划方法及相关装置 | |
CN108940759B (zh) | 连续加工路径的恒定胶量的控制方法及系统 | |
CN109623825B (zh) | 一种移动轨迹规划方法、装置、设备和存储介质 | |
CN112975992B (zh) | 一种误差可控的机器人轨迹同步优化方法 | |
Yang et al. | Smooth trajectory planning along Bezier curve for mobile robots with velocity constraints | |
CN112578671B (zh) | 一种基于u模型优化smc的agv轨迹跟踪控制方法 | |
Oliveira et al. | Trajectory generation using sharpness continuous dubins-like paths with applications in control of heavy-duty vehicles | |
CN111580453A (zh) | 基于plc的提高二维圆弧插补轨迹精度的实现方法 | |
Qinpeng et al. | Path tracking control of wheeled mobile robot based on improved pure pursuit algorithm | |
CN106094737A (zh) | 一种指定加工误差条件下的数控加工速度优化控制方法 | |
Fang et al. | Study on high-speed and smooth transfer of robot motion trajectory based on modified S-shaped acceleration/deceleration algorithm | |
JP4577247B2 (ja) | 移動体及びその制御方法 | |
CN103901815A (zh) | 数值控制装置和数值控制方法 | |
CN115167429A (zh) | 一种移动机器人位置和姿态同时规划方法 | |
CN111267093A (zh) | 一种scara机器人连续轨迹控制实现方法 | |
US7110853B2 (en) | Processes and devices for computer-aided adaptation of an application program for a machine tool | |
CN114253241B (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 |