CN108549324A - 用于高速分拣系统的工件跟随抓取轨迹规划方法及系统 - Google Patents
用于高速分拣系统的工件跟随抓取轨迹规划方法及系统 Download PDFInfo
- Publication number
- CN108549324A CN108549324A CN201810465164.XA CN201810465164A CN108549324A CN 108549324 A CN108549324 A CN 108549324A CN 201810465164 A CN201810465164 A CN 201810465164A CN 108549324 A CN108549324 A CN 108549324A
- Authority
- CN
- China
- Prior art keywords
- workpiece
- crawl
- section
- solution
- locus
- 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
Classifications
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B19/00—Programme-control systems
- G05B19/02—Programme-control systems electric
- G05B19/18—Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form
- G05B19/408—Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form characterised by data handling or data format, e.g. reading, buffering or conversion of data
- G05B19/4086—Coordinate conversions; Other special calculations
-
- G—PHYSICS
- G05—CONTROLLING; REGULATING
- G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
- G05B2219/00—Program-control systems
- G05B2219/30—Nc systems
- G05B2219/35—Nc in input of data, input till input file format
- G05B2219/35356—Data handling
Landscapes
- Engineering & Computer Science (AREA)
- Human Computer Interaction (AREA)
- Manufacturing & Machinery (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Automation & Control Theory (AREA)
- Manipulator (AREA)
- Numerical Control (AREA)
Abstract
本发明公开了一种用于高速分拣系统的工件跟随抓取轨迹规划方法及系统,包括:建立跟随抓取轨迹模型,分别获得工件坐标,机械手末端坐标以及机械手运动参数;根据工件坐标确定工件轨迹所处情况,确定当前求解区间;根据当前求解区间,确定该区间位移与时间公式,确定该区间内单一变量的取值范围;根据工件所处工作范围的上下半区确定抓取轨迹方程;将位移与时间公式代入轨迹方程;使用二分法在取值范围内求解轨迹方程的解,若有解,则根据方程解规划抓取轨迹,执行抓取;若无解,转到下一求解区间重复上述求解过程。本发明有益效果:充分利用机械手工作范围,充分利用S曲线加减速算法特性,提高抓取效率,改善抓取精度。
Description
技术领域
本发明涉及计算机数控系统运动控制领域,尤其涉及一种用于高速分拣系统的工件跟随抓取轨迹规划方法及系统。
背景技术
工件分拣是现代工厂流水线生产线上的一个重要环节。机器人分拣系统一般包含以下设备:机械手,传送带,工件托盘,控制器,工业相机。工业相机通过采集传送带上工件图像获取工件位置信息并将工件信息传递到机器人控制器,控制器进行抓取路径规划与速度规划,然后控制机械手周期插补规划好的抓取轨迹。传统的抓取路径一般采用“门”字型轨迹,在获取抓取指令时,机械手末端首先沿Z轴方向上升固定距离,然后水平运动到传送带上方的预测抓取位置,最后下降至工件上方,再控制吸盘执行抓取动作,吸盘吸气时抓取工件,吸盘吹气则放下工件。此方法的好处是,上升一定高度后再进行水平运动,可以有效避开运动过程中的障碍物,比如传送带边界。但是,吸盘吸气与吹气一般是由控制器控制相应IO完成,IO的控制延迟要远高于机械手插补的1ms,IO控制的延迟导致在机械手到达工件上方时不能及时执行抓取动作,工件继续沿着传送带方向运动一段距离,在IO执行抓取时,抓取位置已经偏离原来的计算位置,导致抓取位置不在工件中心,降低了抓取位置精度,同时放置精度也随之降低;另一方面,由于传送带表面凹凸不平,工件在机器人坐标系内的高度信息存在波动,且难以判断,波动的高度信息给抓取高度的设定带来了极大的困难,为了保证吸盘能够成功吸取工件,吸盘的高度必须尽量靠近工件上表面,从而保证抓取位置精度不受影响。
现有技术提出了根据工件在传送带上的分布密度来调整传送带速度的控制思想,以保证机器人总是处在最快抓取速度状态,并进行了简单的理论分析,但该方法实现较为困难,并且不符合对生产节拍有要求的生产过程。
现有技术提出使用修正梯形算法作为机器人加减速算法,虽然该算法计算简单,但由于修正梯形算法自身时间分配的不合理性,导致它在短时间内的运动具有跳跃性,长距离运动表现疲软的特点,动态特性不足,难以保证高效率的抓取。
综上所述,针对现有技术如何进一步有效动态跟踪运动目标以及提高抓取效率的问题,尚缺乏有效的解决方案。
发明内容
为克服现有技术的不足,本发明提出一种用于高速分拣系统的工件跟随抓取轨迹规划方法及系统,该方法一方面在抓取轨迹的规划过程中,将机器人工作范围和S曲线加减速算法的特点相结合,能够充分利用机器人工作范围;另一方面在保证运动速度的同时,采用圆弧轨迹改变运动方向,实现对于工件的匀速跟随。
为实现上述目的,本发明的具体方案如下:
本发明的第一目的是公开一种用于高速分拣系统的工件跟随抓取轨迹规划方法,包括:
建立跟随抓取轨迹模型,分别获得工件坐标,机械手末端坐标以及机械手运动参数;
根据工件坐标确定工件轨迹所处情况,确定当前求解区间;
根据当前求解区间,确定该区间位移与时间公式,确定该区间内单一变量的取值范围;
根据工件所处工作范围的上下半区确定抓取轨迹方程;
将位移与时间公式代入轨迹方程;
使用二分法在取值范围内求解轨迹方程的解,若有解,则根据方程解规划抓取轨迹,执行抓取;若无解,转到下一求解区间重复上述求解过程。
进一步地,所述建立跟随抓取轨迹模型,具体为:
将机器人的工件抓取轨迹分为四段:竖直上升段,水平直线运动段,水平圆弧过渡段和匀速跟随段;
其中,水平直线运动段采用S曲线加减速算法进行速度规划,通过水平圆弧过渡段实现在保持匀速的前提下将机械手末端运动方向过渡到与工件运动方向保持一致。
进一步地,所述机械手运动参数具体为:起点速度、指令速度、终点速度、最大速度以及待插补的位移。
进一步地,所述的根据工件坐标确定工件轨迹所处情况,确定当前求解区间,具体为:
基于S曲线加减速算法,将机械手的抓取工作范围划分为四个部分;
根据工件的坐标确定工件的运动轨迹;
根据工件的运动轨迹与机械手的抓取工作范围的相交点,将工作范围划分为若干个求解区间。
进一步地,所述基于S曲线加减速算法,将机械手的抓取工作范围划分为四个部分,具体为:
设定S型速度曲线的起点速度vs=0,终点速度ve>0,指令速度vcom,加速度限制amax,加加速度限制Jmax;并且直线插补终点速度ve、圆弧插补速度和传送带速度保持一致;
假设S曲线的最大速度为加速度限制与加加速度限制的比值,判断最大速度是否小于指令速度,如果是,分别计算S曲线加速段位移和减速段位移,得到第一部分的工作半径;
假设S曲线的最大速度为加速度限制与加加速度限制的比值与终点速度之和,判断最大速度是否小于指令速度,如果是,分别计算S曲线加速段位移和减速段位移,得到第二部分的工作半径;
假设S曲线的最大速度为指令速度,分别计算S曲线加速段位移和减速段位移,得到第三部分的工作半径。
进一步地,所述根据当前求解区间,确定该区间位移与时间公式,具体为:
如果当前求解区间属于第一部分或者第二部分或者第三部分抓取工作范围,则该求解区间直线段运动时间为加速段时间和减速段时间之和;直线段位移为最大速度与终点速度的均值与减速段时间的乘积,加上最大速度与起点速度的均值与加速段时间的乘积;上述求解区间内单一变量为最大速度;
如果当前求解区间属于第四部分抓取工作范围,则该求解区间直线段运动时间为加速段时间、匀速段时间和减速段时间之和;此时,直线段位移的最大速度是指令速度,直线段位移为最大速度与终点速度的均值与减速段时间的乘积,加上最大速度与起点速度的均值与加速段时间的乘积,再加上最大速度与匀速段时间的乘积;上述求解区间内单一变量为匀速段时间。
进一步地,根据工件所处工作范围的上半区确定抓取轨迹方程,具体为:
如果跟随运动起点的X坐标小于机械手末端X坐标,则设为工作范围的上半区;
该种情况下的抓取轨迹方程为:
其中,E1(xE1,yE1)为上半区的机械手末端坐标,D1(xD1,yD1)为上半区的工件起始点,θ3为上半区的圆弧插补的圆心角,TE1A1为上半区的S曲线加减速算法规划的直线段时间,LE1A1为上半区的为S曲线加减速算法规划的直线段位移;Trise为抓取轨迹的上升段时间,R为圆弧插补半径,vconveyor为传送带速度。
进一步地,根据工件所处工作范围的下半区确定抓取轨迹方程,具体为:
如果跟随运动起点的X坐标大于机械手末端X坐标,则设为工作范围的下半区;
该种情况下的抓取轨迹方程为:
其中,E2(xE2,yE2)为下半区的机械手末端坐标,D2(xD2,yD2)为下半区的工件起始点,θ6为下半区的圆弧插补的圆心角,TE2A2为下半区的S曲线加减速算法规划的直线段时间,LE2A2为下半区的为S曲线加减速算法规划的直线段位移;Trise为抓取轨迹的上升段时间,R为圆弧插补半径,vconveyor为传送带速度。
本发明的第二目的是公开一种用于高速分拣系统的工件跟随抓取轨迹规划系统,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现以下步骤:
建立跟随抓取轨迹模型,分别获得工件坐标,机械手末端坐标以及机械手运动参数;
根据工件坐标确定工件轨迹所处情况,确定当前求解区间;
根据当前求解区间,确定该区间位移与时间公式,确定该区间内单一变量的取值范围;
根据工件所处工作范围的上下半区确定抓取轨迹方程;
将位移与时间公式代入轨迹方程;
使用二分法在取值范围内求解轨迹方程的解,若有解,则根据方程解规划抓取轨迹,执行抓取;若无解,转到下一求解区间重复上述求解过程。
本发明的第三目的是公开一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时执行以下步骤:
建立跟随抓取轨迹模型,分别获得工件坐标,机械手末端坐标以及机械手运动参数;
根据工件坐标确定工件轨迹所处情况,确定当前求解区间;
根据当前求解区间,确定该区间位移与时间公式,确定该区间内单一变量的取值范围;
根据工件所处工作范围的上下半区确定抓取轨迹方程;
将位移与时间公式代入轨迹方程;
使用二分法在取值范围内求解轨迹方程的解,若有解,则根据方程解规划抓取轨迹,执行抓取;若无解,转到下一求解区间重复上述求解过程。
本发明的有益效果:
本发明将基于S曲线加减速算法的位移分段方法应用于抓取轨迹的规划过程中,将机器人工作范围和S曲线加减速算法的特点相结合,能够充分利用机器人工作范围;
本发明提出基于圆弧过渡的跟随抓取轨迹规划方案,在保证运动速度的同时,采用圆弧轨迹改变运动方向,实现对于工件的匀速跟随。
本发明抓取轨迹方程的公式推导能够证明方程有且仅有一个未知变量,并且采用二分法的数值分析方法验证了该发明提出的抓取轨迹规划方法的可行性。
本发明相比于现有技术,能够充分利用机械手工作范围,充分利用S曲线加减速算法特性,提高抓取效率,改善抓取精度。
附图说明
图1为本发明抓取轨迹规划方法工作流程图;
图2为本发明上半区抓取运动示意图;
图3为本发明下半区抓取运动示意图;
图4为七段S曲线加减速示意图;
图5为位移分段流程图;
图6为S曲线加减速位移分段图;
图7为工件轨迹与机械手工作范围相交情况示意图;
图8为轨迹方程分段求解过程示意图;
图9为本发明机器人分拣系统结构示意图;
图10为本发明机器人分拣系统俯视图。
具体实施方式:
下面结合附图对本发明进行详细说明:
需要注意的是,这里所使用的术语仅是为了描述具体实施方式,而非意图限制根据本申请的示例性实施方式。此外,还应当理解的是,当在本说明书中使用术语“包含”和/或“包括”时,其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
实施例一
为了解决背景技术提出的问题,本发明公开了一种用于高速分拣系统的工件跟随抓取轨迹规划方法,具体过程如图1所示,包括以下步骤:
Step1:工业相机拍照获得工件坐标,控制器获得机械手末端坐标,输入机械手运动参数;
Step2:根据工件坐标确定工件轨迹所处情况,确定当前求解区间;
Step3:根据当前求解区间,确定该区间位移与时间公式,确定该区间内单一变量的取值范围;
Step4:根据工件所处工作范围的上下半区确定抓取轨迹方程;
Step5:将位移与时间公式代入轨迹方程;
Step6:使用二分法在取值范围内求解轨迹方程的解,若有解,则输出方程解转到Step7,若无解,转到下一区间执行Step3-Step6;
Step7:根据方程解规划抓取轨迹,执行抓取。
本发明方法所基于的机器人分拣系统如图9和图10所示,包括:工业相机,四轴高速机械手,吸盘,工件传送带,托盘,托盘传送带以及机器人控制器;工件放置在工件传送带上,工业相机通过采集传送带上工件图像获取工件位置信息并将工件信息传递到机器人控制器,控制器进行抓取路径规划与速度规划,然后控制四轴高速机械手周期插补规划好的抓取轨迹。
下面对上述过程的各个步骤进行详细的介绍。
本实施中,机器人的工件抓取轨迹分为四段:竖直上升段,水平直线运动段,水平圆弧过渡段,匀速跟随段。其中,水平直线运动采用S曲线加减速算法进行速度规划,水平圆弧过渡段是为了在保持匀速的前提下将机械手末端运动方向过渡到与工件运动方向保持一致。
图2和图3代表两种可能出现的抓取轨迹平面示意图。以图1为例,D1点是工件在传送带上运动的起点,E1点时机械手末端运动的起点。D1F1是工件运动轨迹,E1A1是机器人末端直线运动部分轨迹,A1F1是机器人末端圆弧运动部分轨迹,A1是圆弧起点,O1是圆弧圆心点,F1是圆弧终点,也是机器人末端和工件相遇点,F1F2是机器人末端匀速跟随工件运动部分轨迹,机器人末端在F1点时开始下降,直至F2下降至合适抓取高度。
抓取过程:在开始时刻,机械手末端位于点E1,工件位于点D1,并沿着D1F2方向匀速运动,机械手末端首先在E1点竖直上升固定高度,然后沿E1A1段进行直线运动,该段运动基于S曲线加减速算法规划速度,然后沿A1F1段进行圆弧运动,为了实现匀速跟随工件,圆弧插补速度、水平直线段终点速度和传送带速度三者相同。在机械手末端结束圆弧插补到达F1点时,工件也在该时刻运行至圆弧终点F1点,机械手末端位于工件正上方,并跟随工件朝向F2点运动,并在跟随过程中下降高度,运行至F2点时,机械手末端执行抓取动作,然后减速为0。至此便完成该工件抓取任务。
图4所示为S曲线加减速全过程中的速度、加速度以及加加速度曲线。加速过程包含加加速段、匀加速段及减加速段;减速过程包含加减速段、匀减速段及减减速段;加速段和减速段之间通过匀速运动衔接,整个过程分为7段。
本实施例中,定义S型速度曲线的起点速度vs=0,终点速度ve>0,指令速度vcom,且ve<vcom,加速度限制amax,加加速度限制Jmax,S型速度曲线七段运行时间分别为t1,t2,t3,t4,t5,t6,t7。事实上,S曲线加减速可能不完全包含上述7段。根据给定的起点速度、指令速度、终点速度以及待插补的位移S,S曲线可能不包含匀速段、匀加速段或匀减速段。
本实施例中的S曲线加减速算法在使用上有所不同,传统的用法是根据已知的位移,起点速度,终点速度计算各段时间。而本发明假设最大速度,根据已知起点速度和终点速度,计算位移。假设的最大速度决定位移的大小,本发明提出的位移分段原则如图5所示,位移被分成四种情况,分段的关键在于根据假设的最大速度vmax计算相应位移R1、R2、R3。
在本实施例中,S曲线加减速算法应用于抓取轨迹的水平直线段,在执行水平直线插补之前,机械手竖直上升至固定高度,故水平直线段起点速度vs=0,直线段插补完成后,转入圆弧插补和匀速直线跟随,故本发明设定直线插补终点速度ve、圆弧插补速度和传送带速度保持一致。
假设S曲线的最大速度为判断最大速度是否小于指令速度,如果是,分别计算S曲线加速段位移和减速段位移,得到第一部分的工作半径;
假设S曲线的最大速度为判断最大速度是否小于指令速度,如果是,分别计算S曲线加速段位移和减速段位移,得到第二部分的工作半径;
假设S曲线的最大速度为指令速度,分别计算S曲线加速段位移和减速段位移,得到第三部分的工作半径。
因此,如图5所示,R1、R2、R3均存在,S曲线位移被分成①、②、③、④四部分:
第①部分:0<S≤R1、vmax∈(ve,vmax1];
第②部分:R1<S≤R2、vmax∈(vmax1,vmax2];
第③部分:R2<S≤R3、vmax∈(vmax2,vcom];
第④部分:S>R3、vmax=vcom。
根据位移分段的四种情况,绘制S曲线加减速位移分段图如图6所示,结合机械手末端的工作范围,以机械手末端点E为圆心,以R1、R2、R3为半径,设机器人工作范围为半径R4,因为机械手工作范围不是典型的圆形,故使用虚线表示位移范围。
(1)第①部分:0<S≤R1、vmax∈(ve,vmax1];
在该部分,运动过程不包含匀速段、匀加速段和匀减速段,加速段时间tacc和减速段tdec时间计算如下:
直线段运动时间T表示为:
T=tacc+tdec (2)
直线段位移L便可以按照公式计算:
(2)第②部分:R1<S≤R2、vmax∈(vmax1,vmax2];
在该部分,运动过程不包含匀速段,匀减速段,但是包含匀加速段,加速段时间和减速段时间计算如下:
直线段时间可以按照公式(2)计算,位移L按照公式(3)计算。
(3)第③部分:R2<S≤R3、vmax∈(vmax2,vcom];
在该部分,运动过程不包含匀速段,但是包含匀加速段和匀减速段,加速段和减速段计算时间如下:
直线段时间可以按照公式(2)计算,位移L按照公式(3)计算。
(4)第④部分:S>R3、vmax=vcom。
在该部分,运动过程包含匀速段,匀加速段和匀减速段,加速段和减速段计算时间可以按照公式(6)计算,假设tcon为匀速段时间,
直线段时间可以计算为:
T=tacc+tdec+tcon (7)
直线段位移L按照以下公式计算:
综上所述,本发明中由S曲线算法规划的直线段位移和时间均可以由单变量vmax或tcon表示:在第①、②、③部分,根据公式(1)—(6)可得,直线段位移L和时间T公式中有且仅有vmax一个未知变量;在第④部分,根据公式(7)和(8)可得,直线段位移L和时间T公式中有且仅有tcon一个未知变量。
本实施例中,根据跟随运动起点的X坐标大于机械手末端X坐标的大小,对工作范围进行分区,进而建立抓取轨迹计算模型。
如图2所示,本发明设定工件起始点坐标D1(xD1,yD1),机械手末端坐标E1(xE1,yE1),A1F1段起点坐标A1(xA1,yA1),终点坐标F1(xF1,yF1),也是机械手与工件在传送带上的相遇点,跟随段的起点。抓取点坐标为F2(xF2,yF2),A1F1圆心点坐标为O1(xO1,yO1),辅助坐标B1(xB1,yB1)。
同样,在图3中,设定工件起始点坐标D2(xD2,yD2),机械手末端坐标E2(xE2,yE2),A2F3段起点坐标A2(xA2,yA2),终点坐标F3(xF3,yF3),也是机械手与工件在传送带上的相遇点,跟随段的起点。抓取点坐标为F4(xF4,yF4),A2F3圆心点坐标为O2(xO2,yO2),辅助坐标B2(xB2,yB2)。设抓取轨迹的上升段时间为Trise,Trise设为定值。
一般来说,存在两种抓取轨迹,第一种,xF1<xE1,即跟随运动起点的X坐标小于机械手末端X坐标,此时抓取轨迹如图2所示,此种情况该发明设为工作范围的上半区;
第二种,xF3>xE2,即跟随运动起点的X坐标大于机械手末端X坐标,此时抓取轨迹如图3所示,此种情况设为工作范围的下半区。下面分两种情况分别讨论。
(1)工作范围的上半区
首先计算圆弧插补半径,根据向心加速度公式且an=amax,可得圆弧半径公式:
在图2中,设∠θ1=∠A1O1E1,∠θ2=∠B1O1E1,∠θ3=∠A1O1F1,E1A1段插补位移为LE1A1,E1A1段S曲线加减速方法规划的运动时间为TE1A1,A1F1段插补圆弧长度为LA1F1,插补时间为TA1F1.在本问题中,机械手末端点E1(xE1,yE1),工件起始点D1(xD1,yD1),传送带速度vconveyor及运动限制参数vcom,amax,Jmax已知,且对于水平直线段E1A1,vs=0,ve=vconveyor。所以该问题中关键在于求解点F1的坐标,因为yF1=yD1,最终需要求解xF1。
针对机械手末端,机械手末端需要依次插补直线段E1A1和圆弧段A1F1后运动到点F1,点F1的X坐标推导过程如下:
根据角度关系,表示出∠θ1的正切值和∠θ2的余弦值:
根据勾股定理,可得:
其中,
LE1B1=yE1-yD1-R (13)
根据公式(10)、(11)、(12)、(13)可用LE1A1表示cos(∠θ1)和cos(∠θ2)如下
根据图1的X坐标关系,表示出xF1:
然后将公式(10)带入公式(11)可得:
针对工件的运动轨迹,要满足的条件是在机械手末端到达点F1时,工件也必须到达点F1,所以点F1的X坐标推导过程如下:
xF1=xD1+vconveyor(TE1A1+TA1F1+Trise) (18)
其中,TE1A1通过S曲线加减速规划获得,且
其中,
LA1F1=R∠θ3 (20)
并且
综合公式(17)和公式(18)可得方程:
将公式(14)、(15)、(19)带入公式(22)可得:
分析公式(23)不难看出,方程中未知数共有3个:∠θ3、TE1A1和LE1A1,其中∠θ3根据公式(14)、(15)和(21)可用LE1A1表示,方程(23)中存在的未知数只剩下TE1A1和LE1A1,即S曲线算法规划的直线段位移与时间。
(2)工作范围的下半区
圆弧半径计算方法与上半区一致,如图3所示,设∠θ4=∠A2O2E2,∠θ5=∠B2O2E2,∠θ6=∠A2O2F3,E2A2段插补位移为LE2A2,E2A2段S曲线加减速方法规划的运动时间为TE2A2,A2F3段插补圆弧长度为LA2F3,插补时间为TA2F3.在本问题中,机械手末端点E2(xE2,yE2),工件起始点D2(xD2,yD2)。
针对机械手末端,机械手末端需要依次插补直线段E2A2和圆弧段A2F3后运动到点F3,点F3F1的X坐标推导过程如下:
根据角度关系,表示出∠θ4的正切值和∠θ5的余弦值:
其中
LO2B2=yE2-yD2-R (25)
根据公式(24)、(25)、(26)可用LE2A2表示cos(∠θ4)和sin(∠θ5)如下
根据图3的X坐标关系,表示出xF3:
然后可得:
针对工件的运动轨迹,要满足的条件是在机械手末端到达点F3时,工件也必须到达点F3,所以点F3F1的X坐标推导过程如下:
xF3=xD2+vconveyor(TE2A2+TA2F3+Trise)(31)
其中,TE2A2通过S曲线加减速规划获得,且
其中,
LA2F3=R∠θ6 (33)
并且
∠θ6=π-(∠θ4+∠θ5) (34)
综合公式(30)和公式(31)可得方程:
将公式(27)、(28)、(32)带入公式(35)可得:
分析公式(36)不难看出,方程中未知数共有3个:∠θ6、TE2A2和LE2A2,其中∠θ6根据公式(27)、(28)和(34)可用LE2A2表示,方程(36)中存在的未知数只剩下TE2A2和LE2A2,即S曲线算法规划的直线段位移与时间。
综合上半区和下半区推导过程可得,抓取轨迹计算方程求解的关键在于求解S曲线加减速算法规划的直线段时间T和位移L。
方程(23)或(36)给出了计算抓取轨迹的计算模型,模型中未知变量为抓取轨迹水平运动直线段的运动时间T和位移L,同时,由于直线段的运动时间T和位移L可以由单一变量vmax或tcon表示,因此,抓取轨迹方程可以通过数值分析vmax或tcon的方法求解。本发明通过采用二分法分析变量vmax或tcon来求解轨迹方程,且本发明提出的S曲线加减速算法的位移分段方法提供了一种讨论依据,下面根据工件轨迹的分布详细讨论求解方法。
如图6所示,基于S曲线加减速算法,本发明将机械手的抓取工作范围划分为四个部分。另一方面,工件在传送带上沿X轴正方向以匀速度运动,工件的Y坐标决定工件在传送带上的轨迹,也能确定与图6所示工作范围的交点个数。其中,ai(i=1,2,3,4)是四种轨迹的相交情况的上限,也是抓取工作范围的上限;hi(i=1,2,3,4)是四种轨迹的相交情况的下限,也是抓取工作范围的下限。
所以,根据工件的Y坐标可以将工件轨迹分为图7所示的四种情况:
(1)工件轨迹属于情况α
在此情况下,工件轨迹与机械手工作空间的四个部分(①、②、③、④)都相交,a1、b1、c1、d1、f1、g1、h1是两者的相交点。
(2)工件轨迹属于情况β
在此情况下,工件轨迹与机械手工作范围的三个部分(②、③、④)相交,a2、b2、c2、f2、g2、h2是两者的相交点。
(3)工件轨迹属于情况γ
在此情况下,工件轨迹与机械手工作范围的两个部分(③、④)相交,a3、b3、g3、h3是两者的相交点。
(4)工件轨迹属于情况δ
在此情况下,工件轨迹与机械手的工作范围只在第④部分相交,a4、h4是两者的相交点。
本实施例中,在实际抓取过程中,为了提高抓取效率,当工件进入机械手工作范围内,机械手便开始规划抓取轨迹。本发明根据工件轨迹与工作范围的相交情况,将工作范围划分为若干个求解区间,每个区间对应不同的求解公式。当工件进入工作范围后,根据工件X坐标由小到大的方向依次根据工件所在区间求解轨迹方程,若当前区间内无解,则转向下一区间,直至求出轨迹方程的解。求解过程如图8所示。
实施例二
本实施例中,公开了一种用于高速分拣系统的工件跟随抓取轨迹规划系统,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现以下步骤:
建立跟随抓取轨迹模型,分别获得工件坐标,机械手末端坐标以及机械手运动参数;
根据工件坐标确定工件轨迹所处情况,确定当前求解区间;
根据当前求解区间,确定该区间位移与时间公式,确定该区间内单一变量的取值范围;
根据工件所处工作范围的上下半区确定抓取轨迹方程;
将位移与时间公式代入轨迹方程;
使用二分法在取值范围内求解轨迹方程的解,若有解,则根据方程解规划抓取轨迹,执行抓取;若无解,转到下一求解区间重复上述求解过程。
实施例三
本实施例中,公开了一种计算机可读存储介质,其上存储有计算机程序,该程序被处理器执行时执行以下步骤:
建立跟随抓取轨迹模型,分别获得工件坐标,机械手末端坐标以及机械手运动参数;
根据工件坐标确定工件轨迹所处情况,确定当前求解区间;
根据当前求解区间,确定该区间位移与时间公式,确定该区间内单一变量的取值范围;
根据工件所处工作范围的上下半区确定抓取轨迹方程;
将位移与时间公式代入轨迹方程;
使用二分法在取值范围内求解轨迹方程的解,若有解,则根据方程解规划抓取轨迹,执行抓取;若无解,转到下一求解区间重复上述求解过程。
上述虽然结合附图对本发明的具体实施方式进行了描述,但并非对本发明保护范围的限制,所属领域技术人员应该明白,在本发明的技术方案的基础上,本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本发明的保护范围以内。
Claims (10)
1.用于高速分拣系统的工件跟随抓取轨迹规划方法,其特征在于,包括:
建立跟随抓取轨迹模型,分别获得工件坐标,机械手末端坐标以及机械手运动参数;
根据工件坐标确定工件轨迹所处情况,确定当前求解区间;
根据当前求解区间,确定该区间位移与时间公式,确定该区间内单一变量的取值范围;
根据工件所处工作范围的上下半区确定抓取轨迹方程;
将位移与时间公式代入轨迹方程;
使用二分法在取值范围内求解轨迹方程的解,若有解,则根据方程解规划抓取轨迹,执行抓取;若无解,转到下一求解区间重复上述求解过程。
2.如权利要求1所述的用于高速分拣系统的工件跟随抓取轨迹规划方法,其特征在于,所述建立跟随抓取轨迹模型,具体为:
将机器人的工件抓取轨迹分为四段:竖直上升段,水平直线运动段,水平圆弧过渡段和匀速跟随段;
其中,水平直线运动段采用S曲线加减速算法进行速度规划,通过水平圆弧过渡段实现在保持匀速的前提下将机械手末端运动方向过渡到与工件运动方向保持一致。
3.如权利要求1所述的用于高速分拣系统的工件跟随抓取轨迹规划方法,其特征在于,所述机械手运动参数具体为:起点速度、指令速度、终点速度、最大速度以及待插补的位移。
4.如权利要求1所述的用于高速分拣系统的工件跟随抓取轨迹规划方法,其特征在于,所述的根据工件坐标确定工件轨迹所处情况,确定当前求解区间,具体为:
基于S曲线加减速算法,将机械手的抓取工作范围划分为四个部分;
根据工件的坐标确定工件的运动轨迹;
根据工件的运动轨迹与机械手的抓取工作范围的相交点,将工作范围划分为若干个求解区间。
5.如权利要求4所述的用于高速分拣系统的工件跟随抓取轨迹规划方法,其特征在于,所述基于S曲线加减速算法,将机械手的抓取工作范围划分为四个部分,具体为:
设定S型速度曲线的起点速度vs=0,终点速度ve>0,指令速度vcom,加速度限制amax,加加速度限制Jmax;并且直线插补终点速度ve、圆弧插补速度和传送带速度保持一致;
假设S曲线的最大速度为加速度限制与加加速度限制的比值,判断最大速度是否小于指令速度,如果是,分别计算S曲线加速段位移和减速段位移,得到第一部分的工作半径;
假设S曲线的最大速度为加速度限制与加加速度限制的比值与终点速度之和,判断最大速度是否小于指令速度,如果是,分别计算S曲线加速段位移和减速段位移,得到第二部分的工作半径;
假设S曲线的最大速度为指令速度,分别计算S曲线加速段位移和减速段位移,得到第三部分的工作半径。
6.如权利要求1所述的用于高速分拣系统的工件跟随抓取轨迹规划方法,其特征在于,所述根据当前求解区间,确定该区间位移与时间公式,具体为:
如果当前求解区间属于第一部分或者第二部分或者第三部分抓取工作范围,则该求解区间直线段运动时间为加速段时间和减速段时间之和;直线段位移为最大速度与终点速度的均值与减速段时间的乘积,加上最大速度与起点速度的均值与加速段时间的乘积;上述求解区间内单一变量为最大速度;
如果当前求解区间属于第四部分抓取工作范围,则该求解区间直线段运动时间为加速段时间、匀速段时间和减速段时间之和;直线段位移为最大速度与终点速度的均值与减速段时间的乘积,加上最大速度与起点速度的均值与加速段时间的乘积,再加上最大速度与匀速段时间的乘积;上述求解区间内单一变量为匀速段时间。
7.如权利要求1所述的用于高速分拣系统的工件跟随抓取轨迹规划方法,其特征在于,根据工件所处工作范围的上半区确定抓取轨迹方程,具体为:
如果跟随运动起点的X坐标小于机械手末端X坐标,则设为工作范围的上半区;
该种情况下的抓取轨迹方程为:
其中,E1(xE1,yE1)为上半区的机械手末端坐标,D1(xD1,yD1)为上半区的工件起始点,θ3为上半区的圆弧插补的圆心角,TE1A1为上半区的S曲线加减速算法规划的直线段时间,LE1A1为上半区的为S曲线加减速算法规划的直线段位移;Trise为抓取轨迹的上升段时间,R为圆弧插补半径,vconveyor为传送带速度。
8.如权利要求1所述的用于高速分拣系统的工件跟随抓取轨迹规划方法,其特征在于,根据工件所处工作范围的下半区确定抓取轨迹方程,具体为:
如果跟随运动起点的X坐标大于机械手末端X坐标,则设为工作范围的下半区;
该种情况下的抓取轨迹方程为:
其中,E2(xE2,yE2)为下半区的机械手末端坐标,D2(xD2,yD2)为下半区的工件起始点,θ6为下半区的圆弧插补的圆心角,TE2A2为下半区的S曲线加减速算法规划的直线段时间,LE2A2为下半区的为S曲线加减速算法规划的直线段位移;Trise为抓取轨迹的上升段时间,R为圆弧插补半径,vconveyor为传送带速度。
9.用于高速分拣系统的工件跟随抓取轨迹规划系统,其特征在于,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述程序时实现以下步骤:
建立跟随抓取轨迹模型,分别获得工件坐标,机械手末端坐标以及机械手运动参数;
根据工件坐标确定工件轨迹所处情况,确定当前求解区间;
根据当前求解区间,确定该区间位移与时间公式,确定该区间内单一变量的取值范围;
根据工件所处工作范围的上下半区确定抓取轨迹方程;
将位移与时间公式代入轨迹方程;
使用二分法在取值范围内求解轨迹方程的解,若有解,则根据方程解规划抓取轨迹,执行抓取;若无解,转到下一求解区间重复上述求解过程。
10.一种计算机可读存储介质,其上存储有计算机程序,其特征在于,该程序被处理器执行时执行以下步骤:
建立跟随抓取轨迹模型,分别获得工件坐标,机械手末端坐标以及机械手运动参数;
根据工件坐标确定工件轨迹所处情况,确定当前求解区间;
根据当前求解区间,确定该区间位移与时间公式,确定该区间内单一变量的取值范围;
根据工件所处工作范围的上下半区确定抓取轨迹方程;
将位移与时间公式代入轨迹方程;
使用二分法在取值范围内求解轨迹方程的解,若有解,则根据方程解规划抓取轨迹,执行抓取;若无解,转到下一求解区间重复上述求解过程。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810465164.XA CN108549324B (zh) | 2018-05-16 | 2018-05-16 | 用于高速分拣系统的工件跟随抓取轨迹规划方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810465164.XA CN108549324B (zh) | 2018-05-16 | 2018-05-16 | 用于高速分拣系统的工件跟随抓取轨迹规划方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108549324A true CN108549324A (zh) | 2018-09-18 |
CN108549324B CN108549324B (zh) | 2019-07-05 |
Family
ID=63494950
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810465164.XA Active CN108549324B (zh) | 2018-05-16 | 2018-05-16 | 用于高速分拣系统的工件跟随抓取轨迹规划方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108549324B (zh) |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109623821A (zh) * | 2018-12-26 | 2019-04-16 | 深圳市越疆科技有限公司 | 机器手抓取物品的视觉引导方法 |
CN110480630A (zh) * | 2019-07-16 | 2019-11-22 | 绍兴文理学院 | 一种工业机器人的自适应目标定位系统及方法 |
CN112068586A (zh) * | 2020-08-04 | 2020-12-11 | 西安交通大学 | 一种空间时间联合优化的四旋翼无人飞行器轨迹规划方法 |
CN112476437A (zh) * | 2020-12-08 | 2021-03-12 | 北京航天新立科技有限公司 | 一种柔性的自动抓取设备控制方法 |
CN113253678A (zh) * | 2021-04-29 | 2021-08-13 | 中国科学院微电子研究所 | 精密运动平台的运动轨迹规划方法 |
CN113941783A (zh) * | 2021-11-23 | 2022-01-18 | 深圳软动智能控制有限公司 | 激光轴蛙跳实现方法、装置、激光设备和存储介质 |
CN114078162A (zh) * | 2022-01-19 | 2022-02-22 | 湖南视比特机器人有限公司 | 钢板切割后工件的桁架分拣方法及系统 |
CN114227655A (zh) * | 2021-12-30 | 2022-03-25 | 深圳市英威腾电气股份有限公司 | 确定规划路径的方法、装置、scara机器人及介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104240263A (zh) * | 2014-07-30 | 2014-12-24 | 华南理工大学 | 一种用于Delta并联机械手的运动特征分析系统 |
CN105014666A (zh) * | 2015-07-13 | 2015-11-04 | 广州霞光技研有限公司 | 一种多自由度机械手自主抓取逆解工程算法 |
CN106926233A (zh) * | 2015-12-29 | 2017-07-07 | 楚天科技股份有限公司 | 一种平面机械手运动路径的规划方法 |
CN107671859A (zh) * | 2017-10-20 | 2018-02-09 | 山东大学 | 基于s曲线加减速的运动目标动态跟踪方法、系统及装置 |
-
2018
- 2018-05-16 CN CN201810465164.XA patent/CN108549324B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104240263A (zh) * | 2014-07-30 | 2014-12-24 | 华南理工大学 | 一种用于Delta并联机械手的运动特征分析系统 |
CN105014666A (zh) * | 2015-07-13 | 2015-11-04 | 广州霞光技研有限公司 | 一种多自由度机械手自主抓取逆解工程算法 |
CN106926233A (zh) * | 2015-12-29 | 2017-07-07 | 楚天科技股份有限公司 | 一种平面机械手运动路径的规划方法 |
CN107671859A (zh) * | 2017-10-20 | 2018-02-09 | 山东大学 | 基于s曲线加减速的运动目标动态跟踪方法、系统及装置 |
Non-Patent Citations (1)
Title |
---|
倪鹤鹏 等: "基于机器视觉的Delta机器人分拣系统算法", 《机器人》 * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109623821A (zh) * | 2018-12-26 | 2019-04-16 | 深圳市越疆科技有限公司 | 机器手抓取物品的视觉引导方法 |
CN109623821B (zh) * | 2018-12-26 | 2022-04-01 | 日照市越疆智能科技有限公司 | 机械手抓取物品的视觉引导方法 |
CN110480630A (zh) * | 2019-07-16 | 2019-11-22 | 绍兴文理学院 | 一种工业机器人的自适应目标定位系统及方法 |
CN112068586A (zh) * | 2020-08-04 | 2020-12-11 | 西安交通大学 | 一种空间时间联合优化的四旋翼无人飞行器轨迹规划方法 |
CN112476437A (zh) * | 2020-12-08 | 2021-03-12 | 北京航天新立科技有限公司 | 一种柔性的自动抓取设备控制方法 |
CN113253678A (zh) * | 2021-04-29 | 2021-08-13 | 中国科学院微电子研究所 | 精密运动平台的运动轨迹规划方法 |
CN113941783A (zh) * | 2021-11-23 | 2022-01-18 | 深圳软动智能控制有限公司 | 激光轴蛙跳实现方法、装置、激光设备和存储介质 |
CN113941783B (zh) * | 2021-11-23 | 2023-11-21 | 深圳软动智能控制有限公司 | 激光轴蛙跳实现方法、装置、激光设备和存储介质 |
CN114227655A (zh) * | 2021-12-30 | 2022-03-25 | 深圳市英威腾电气股份有限公司 | 确定规划路径的方法、装置、scara机器人及介质 |
CN114078162A (zh) * | 2022-01-19 | 2022-02-22 | 湖南视比特机器人有限公司 | 钢板切割后工件的桁架分拣方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN108549324B (zh) | 2019-07-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108549324B (zh) | 用于高速分拣系统的工件跟随抓取轨迹规划方法及系统 | |
CN107671859B (zh) | 基于s曲线加减速的运动目标动态跟踪方法、系统及装置 | |
CN106313047B (zh) | 一种基于Bezier样条的机器人实时拐角过渡方法 | |
CN110653826B (zh) | 一种面向传送带目标跟随的机器人实时轨迹规划方法 | |
CN105883616B (zh) | 桥式吊车最短时间防摆轨迹实时生成方法 | |
CN111158318B (zh) | 一种非对称性四次曲线柔性加减速规划方法 | |
CN107160389B (zh) | 一种工业机器人的力矩控制方法 | |
CN106945042A (zh) | 一种机械手的离散s型曲线速度控制方法 | |
CN104589357B (zh) | 基于视觉跟踪的delta机器人控制系统和方法 | |
CN107479581A (zh) | 货箱码放控制装置、系统、方法和存储介质 | |
CN107505918A (zh) | 一种切割机的速度规划方法 | |
CN102172774B (zh) | 一种选择性激光烧结扫描方法 | |
CN106227154B (zh) | 一种两轴联动机械臂运动控制的同步误差补偿方法 | |
CN110147077B (zh) | 一种工业机器人在空间直线运行下的余弦插补方法 | |
CN104570735A (zh) | 一种码垛机器人弧形运动轨迹算法 | |
CN106444645A (zh) | 一种基于嵌入式裁床控制器的多轴联动运动控制方法 | |
CN113031621B (zh) | 一种桥式起重机安全避障路径规划方法及系统 | |
CN106743340A (zh) | 物流上包控制方法及控制系统 | |
CN110377065A (zh) | 速度控制方法、装置及计算机可读存储介质 | |
CN107807681A (zh) | 一种基于编码器的传动设备自动定位方法及装置 | |
CN109271689A (zh) | 圆轨迹下空间包络成形包络模轨迹精确计算方法 | |
CN111630461B (zh) | 数控加工方法及系统、具有存储功能的装置 | |
CN110398993A (zh) | 速度控制方法、装置及计算机可读存储介质 | |
CN117466145A (zh) | 基于深度强化学习的桥门式起重机防摇控制方法及装置 | |
CN109842333B (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 |