CN106503373B - 一种基于b样条曲线的双机器人协调装配轨迹规划方法 - Google Patents
一种基于b样条曲线的双机器人协调装配轨迹规划方法 Download PDFInfo
- Publication number
- CN106503373B CN106503373B CN201610965825.6A CN201610965825A CN106503373B CN 106503373 B CN106503373 B CN 106503373B CN 201610965825 A CN201610965825 A CN 201610965825A CN 106503373 B CN106503373 B CN 106503373B
- Authority
- CN
- China
- Prior art keywords
- assembly
- robot
- master
- obstacle avoidance
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 45
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 28
- 238000005457 optimization Methods 0.000 claims abstract description 12
- 239000002245 particle Substances 0.000 claims description 34
- 238000001514 detection method Methods 0.000 claims description 15
- 101100129590 Schizosaccharomyces pombe (strain 972 / ATCC 24843) mcp5 gene Proteins 0.000 claims description 12
- 230000009466 transformation Effects 0.000 claims description 8
- 239000011159 matrix material Substances 0.000 claims description 6
- 230000001133 acceleration Effects 0.000 claims description 5
- 238000011156 evaluation Methods 0.000 claims description 5
- 230000036461 convulsion Effects 0.000 claims description 4
- 238000006073 displacement reaction Methods 0.000 claims description 4
- 230000010349 pulsation Effects 0.000 claims description 3
- 241000009298 Trigla lyra Species 0.000 claims description 2
- 238000006243 chemical reaction Methods 0.000 claims 2
- 230000006870 function Effects 0.000 description 8
- 239000013598 vector Substances 0.000 description 8
- 230000009977 dual effect Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 238000013178 mathematical model Methods 0.000 description 3
- 230000004888 barrier function Effects 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 2
- 238000011426 transformation method Methods 0.000 description 2
- 230000005428 wave function Effects 0.000 description 2
- 238000000342 Monte Carlo simulation Methods 0.000 description 1
- 230000000712 assembly Effects 0.000 description 1
- 238000000429 assembly Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 238000005315 distribution function Methods 0.000 description 1
- 230000009191 jumping Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
- G06Q10/047—Optimisation of routes or paths, e.g. travelling salesman problem
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Pure & Applied Mathematics (AREA)
- Entrepreneurship & Innovation (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Mathematical Optimization (AREA)
- Development Economics (AREA)
- Mathematical Analysis (AREA)
- Game Theory and Decision Science (AREA)
- Computational Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Marketing (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Manipulator (AREA)
- Numerical Control (AREA)
Abstract
本发明提供了一种基于B样条曲线的双机器人协调装配轨迹规划方法,属于轨迹规划领域;其包括有任务模块、主‑从机器人关节轨迹规划模块、避障策略模块、多目标优化模块、工业机器人驱动器。所述主‑从机器人关节轨迹规划模块发明了一种基于B样条曲线的双机器人协调装配关节轨迹规划方法和协调装配误差补偿的方法;所述避障策略模块发明了一种双机器人深度搜索避障策略;多目标优化模块发明了一种基于群体智能算法求解双机器人协调装配时间较少,装配轨迹平滑的运动学约束条件下的多目标规划方法;采用本发明的技术方案,可以得到装配时间较短,装配段关节轨迹平滑,成功避开障碍物且无干涉的双机器人协调装配轨迹,降低了成本,避免了装配件的损坏。
Description
技术领域
本发明涉及轨迹规划领域。更具体地,涉及一种基于B样条曲线的双机器人协调装配的轨迹规划方法。
背景技术
由于单机器人装配存在工作空间有限,容易存在装配盲区,灵活性不够,难以实现全自动化装配作业等不足。因此,为了适应装配任务的复杂性和智能型,弥补单机器人装配的不足,提出双机器人协调装配***。双机器人协调装配***虽然具有灵活性好、适应性高的优点,但其复杂性也随之增加,双机器人***的协调运动方案、机器人之间的避障算法、双机器人的轨迹规划及其优化算法都是亟待解决的问题。
发明内容
为了解决上述问题,本发明旨在公开一种双机器人协调装配轨迹规划方法,尤指一种基于B样条曲线获得运行时间最少、装配轨迹平滑、避开障碍物的多目标双机器人装配轨迹规划的方法。
为实现上述目的,本发明所采用的技术方案主要包括以下步骤:
步骤1、建立双工业机器人的运动学模型,求得主-从工业机器人的正反解;
所述主机器人是指在协调装配阶段,关节位移、速度和加速度根据用运动学约束来确定的夹持一个装配体的工业机器人;从机器人是指在协调装配阶段,关节运动随主机器人的关节运动而变化的夹持另一个装配体的工业机器人。
步骤2、由所述的任务模块获取主-从机器人的任务点位姿信息;
所述任务模块是指给定主-从机器人具体的离散示教路径点,利用步骤1的结果可以得到主-从机器人的任务点位姿信息。
步骤3、建立主-从机器人协调***的坐标系,通过坐标转换的方法得到主-从机器人的末端坐标系相对于自身基座标系的变换矩阵。
步骤4、采用B样条曲线插值的方法对主-从机器人协调装配关节位置点进行插值,求得双机器人协调装配关节轨迹,并进行协调装配误差的补偿,即关节轨迹规划模块;
所述B样条曲线插值方法的数学模型为:
公式(2)中:Q1(u)∈RM×1和Q2(u)∈RM×1分别为归一化时刻u处主-从机器人关节的位移,其中M为机器人的关节数目,dj,1(u)∈RM×1和dj,2(u)∈RM×1分别为主-从机器人B样条关节轨迹曲线的控制点;为主从机器人B样条关节轨迹曲线的归一化时间向量,n1+1和n2+1分别为主从机器人的任务插值点个数,Nj,k,1(u)和Nj,k,2(u)分别为主-从机器人B样条关节轨迹曲线的k次规范B样条基函数,按照递推公式(7)定义:
公式(3)中,R=1表示主机器人,R=2表示从机器人,i为样条序号,B样条基函数Ni,k,R(u)的支承区间为[ui,ui+k+1],i=0,1,...n;若运行时间为ttotal,根据主-从机器人关节轨迹的B样条插值曲线的控制点和节点矢量,将的德布尔递推公式变换可以得到曲线上时刻t的r阶导矢,即:
当公式(4)中R=1,2、r=0,1,2,3时可以分别得到主-从机器人B样条关节轨迹随时间变化的关节位移、速度、加速度和加加速度曲线;
由于装配体之间的位姿约束和相对运动约束难以保证,会产生装配误差,因此基于B样条曲线插值后需要进行装配误差补偿;
所述协调装配误差补偿方法为:首先定义装配的误差建立装配体之间的配合坐标系,以装配体A的轴线为Z轴,装配体A在轴线方向装配路径段的起始点O1为原点建立装配体A的坐标系;以无装配误差时装配体B的配合面轴线为Z轴,装配体B在轴线方向装配路径段的起始点O2为原点建立装配体B的坐标系,建立轴孔简化模型,如图2所示,无误差装配情况下装配体B坐标系的Z轴始终与装配体A坐标系的Z轴重合,装配体A和装配体B之间相对运动约束为:装配体A沿着装配体B坐标系的Z轴做直线运动;根据装配任务要求在笛卡尔空间给定一系列相对世界坐标系的装配任务点序列:开始进入协调装配时O1与O2重合;装配完成后O1与Pn重合,O2与P1重合;图2装配体A的原点O1在装配体B的Z轴上的垂足;F为装配体B的原点O2在WP1的Z轴上的垂足;为了在双机器人协调装配策略中消除装配误差,提高装配成功率,定义装配位置误差和装配姿态误差;在装配路径段,选取两垂线距离(O1,E)至(O2,F)距离的最大值定义为装配位置误差PositionError,将当前装配体A和装配体B之间姿态变换的欧拉角与理想欧拉角的差值取绝对值求和定义为装配体定义装配姿态误差AttitudeError,具体为:
公式(5)中PositionError为位置误差,AttitudeError为姿态误差,分别为两垂线距离(参考图2),α、β、γ、α*、β*、γ*为欧拉角。理想装配情况下,装配体B坐标系的Z轴始终与装配体A坐标系的Z轴重合,所以装配位置误差PositionError为0;当前装配体A和装配体B之间装配姿态误差AttitudeError也为0。
步骤6、执行双机器人深度搜索避障策略模块;
所述双机器人深度搜索避障策略中机器人与障碍物之间的避障策略如下:在机器人的连杆轴线上每隔一定间距设置一个干涉检测点,并将干涉检测点各个区间段连杆的最大包络半径叠加到障碍物的厚度上,从而将机器人与障碍物之间的碰撞转化为由一系列干涉检测点与扩充后的障碍物间的位置关系判定。
所述一系列干涉检测点与扩充后的障碍物间的位置关系判定方法如下:将机器人简化为一系列干涉检测点后,判断整个工作时间内干涉检测点是否在扩充后的障碍物区域内;若是,则连杆和障碍物有干涉情况,判定机器人与障碍物发生碰撞;否则继续判断其他干涉检测点在整个工作时间内是否在扩充后的障碍物区域内;若所有干涉检测点在整个工作时间内都不在扩充后的障碍物区域内,则判定机器人与障碍物没有发生碰撞,否则机器人与障碍物发生碰撞(参考图3)。
所述双机器人深度搜索避障策略中机器人与机器人之间的避障策略如下:将机器人简化为由连杆轴组成的直线段,从而将机器人与机器人之间的碰撞简化为判断整个工作空间内空间中两条直线段的最短距离与相应两条连杆的最大包络半径之和的大小判定关系。
所述简化判断整个工作空间内空间中两条直线段的最短距离与相应两条连杆的最大包络半径之和的大小判定关系方法如下:计算主机器人连杆i的简化直线段与从机器人连杆j的简化直线段之间的最短距离(计算方法参考空间中两条线段之间的最短距离的计算方法)。判断在整个工作时间内最短距离是否小于连杆i和连杆j的最大包络半径之和,若是,则连杆i与连杆j发生碰撞,判定双机器人之间发生碰撞;否则连杆i与连杆j没有发生碰撞,继续判断主机器人其他连杆i与从机器人的连杆是否发生碰撞;当主机器人所有连杆与从机器人所有连杆都没有发生碰撞,则判定双机器人之间没有发生碰撞(参考图4)。
所述双机器人深度搜索避障策略如下步骤:
S1、设置主-从机器人各自的起始点与目标点,令num1=num2=0;
S2、根据任务空间中机器人与障碍物存在的位置关系人为确定大致的避障中间点位置,num1=num1+1;
S3、在避障中间点附近指定的搜索范围内进行启发式搜索,确定避障中间点,num2=num2+1;
S4、分别对主-从机器人的起始点、中间点和目标点进行B样条插值,生成主-从机器人的关节轨迹;
S5、通过避障检测方法判断整个工作时间内机器人和障碍物、双机器人之间是否发生碰撞,若否,则输出主-从机器人的无碰关节轨迹,fob=0终止算法;若是,则执行步骤S6;
S6、判断num2是否小于num2max,若是,则执行步骤S3,若否,则执行步骤S7;
S7、判断num1是否小于num1max,若是,则执行步骤S2,若否,则没有成功搜索到主-从机器人的无碰撞关节轨迹,fob=+∞,算法结束;
其中num1max为人为确定避障中间点位置的尝试次数,num2max为在避障中间点附近范围内随机确定避障中间点的最大允许次数,具体的算法执行流程参考图5。
步骤5、确定协调装配轨迹的优化目标和约束条件,结合量子粒子群算法和避障策略进行求解,即多目标规划模块;
所述多目标数学模型为:
目标函数:FG=ω1.ft+ω2.fj+ω3.fob (6)
约束条件:
公式(6~10)中:ft为运行时间评价指标,衡量装配效率;fj为关节平均脉动,作为轨迹的平滑度评价指标;fob为碰撞评价指标;ω1、ω2、ω3分别为ft、fj、fob三个评价指标的权值系数。
所述量子粒子群算法数学模型为:
假设N是种群粒子数,粒子i的当前位置记为Xi=(xi1,xi2,…,xiD),D是粒子的维数。在每一次迭代中,粒子i搜索的最好位置记为pbesti=(pbesti1,pbesti2,…,pbestiD),整个种群搜索的最好位置记为gbesti=(gbesti1,gbesti2,…,gbestiD)。其中:
公式(11~12)中t是算法的当前迭代次数,r1d(t)和r2d(t)是(0,1)间的随机数;c1和c2为PSO算法的学习因子,一般取c1=c2,此时是均匀分布在(0,1)上的随机数;在粒子i的当前最优位置pbesti和群体的全局最优位置gbesti为顶点构成的超矩形中按比例随机选取一个位置作为粒子i的下一代的最好位置;此方法称为粒子群算法。
假设粒子群中粒子具有量子动力学的行为,由波函数ψ来确定每个粒子的位置,|ψ|2是粒子位置的概率密度函数。假定在第t次迭代,粒子i在D维空间运动,该粒子在第d维的势阱为pbestid(t),则可以得到粒子i在第(t+1)次迭代的波函数为:
则粒子的概率密度函数Q为:
概率分布函数T为:
采用蒙特卡罗方法对粒子位置进行随机采样,得到在第(t+1)次迭代时,第i个粒子第d维的位置分量:
公式(16)中uid(t)是均匀分布在(0,1)上的随机数。Lid(t)的值由公式(17)确定:
Lid(t)=2α(t)|mbestd(t)-xid(t)| (17)
其中mbest称为平均最优位置,它是所有粒子自身最优位置的中心点,可由式(18)计算得到:
则粒子的位置更新公式为:
xid(t+1)=pbestid(t)±α(t)|mbestid(t)-xid(t)|ln[1/uid(t)] (19)
公式(19)中,参数α称为压缩-扩张因子,直接影响粒子的收敛速度。该参数的值可以是固定的,也可以是动态变化的,通常由下式确定:
公式(20)中a<b,b和a分别是控制参数α的初始值和终止值。
粒子的当前最优位置pbesti的更新方式为:
全局最优位置gbest的更新方式为:
公式(22)中f为目标函数。将公式(19)作为粒子位置更新公式的粒子群算法即为具有量子粒子群优化算法。
所述结合量子粒子群算法和避障策略进行求解方法流程参考图6。
步骤7、将多目标规划模块得到的关节空间信息发送给主-从机器人的驱动器。
采用上述双机器人协调装配轨迹规划方法,可以得到一条速度、加速度、加加速度连续,装配时间较短,装配段关节轨迹平滑,成功避开障碍物且无干涉的双机器人协调装配轨迹,降低了成本,避免了装配件的损坏。
附图说明
图1是本发明提出的基于量子粒子群算法和B样条曲线的双机器人协调装配的轨迹规划方法的流程图;
图2是本发明提出的基于B样条曲线的双机器人协调装配的轨迹规划方法的协调装配误差定义示意图;
图3是本发明提出的基于B样条曲线的双机器人协调装配的轨迹规划方法的机器人与障碍物之间避障简化示意图;
图4是本发明提出的基于B样条曲线的双机器人协调装配的轨迹规划方法的机器人与机器人之间避障简化示意图;
图5是本发明提出的基于B样条曲线的双机器人协调装配的轨迹规划方法的双机器人深度搜索避障策略流程图;
图6是本发明提出的基于B样条曲线的双机器人协调装配的轨迹规划方法的量子粒子群算法执行的流程图;
具体实施方式
为了使本发明的目的、技术方案和优点更加清楚,下面结合附图详细说明本发明的具体实施方式:
参考图1,本发明公开了一种基于B样条曲线的双工业机器人协调装配的轨迹规划方法,其包括有任务模块、主-从机器人关节轨迹规划模块、避障策略模块、多目标优化模块、工业机器人驱动器。以下就以任务为实现典型的双工业机器人(六自由度)轴孔协调装配任务为例进行描述。其他符合要求的任务均可按照本方案实施,具体的实施主要步骤如下:
步骤1、建立双工业机器人的运动学模型,求得主-从工业机器人的正反解;
所述运动学模型是指可以通过D-H方法分别建立主-从机器人各个连杆坐标系,通过D-H参数表和坐标系变换矩阵推导出运动学正解,然后通过PIPER法推导出运动学逆解。
步骤2、由所述的任务模块获取主-从机器人的任务点位姿信息。
双机器人协调装配任务一般分为两个阶段:第一阶段为夹取工件;第二阶段为协调装配阶段;第一阶段只需要双机器人各自绕过障碍物夹取装配体并且运动到指定的装配位置;第二阶段需要双机器人在保持既定的相对装配约束下必须严格经过离散的装配路径点。通过示教获取离散装配路径点,并通过步骤1得到的双机器人运动学逆解公式去求任务点姿态。
步骤3、建立主-从机器人协调***的坐标系,通过坐标转换的方法得到主-从机器人的末端坐标系相对于自身基座标系的变换矩阵。
所述主-从机器人协调***包括双机器人的基座标系、末端坐标系、工具坐标系、装配任务点坐标系;主-从机器人组成两条运动链所形成的封闭环链:一条是“主机器人的基座-主机器人末端—装配体A-装配任务点”的运动链,另外一条是“主机器人基座-从机器人基座-从机器人末端-装配体B-装配任务点”的运动链,根据两条运动链,可以得到主-从机器人末端坐标系相对于自身基座标系的变换矩阵,进而根据主-从机器人的逆运动学求解出相应的关节角,从而驱使双机器人进行协调装配运动。
步骤4、采用B样条曲线插值的方法对主-从机器人协调装配关节位置点进行插值,求得双机器人协调装配关节轨迹,并进行误差补偿,即关节轨迹规划模块;
根据步骤1~3可以得到主-从机器人装配任务路径点相对应的关节角度,为了实现主-从机器人关节轨迹的平滑性,采用B样条曲线插值的方法进行插值。具体插值步骤为:
S3、结合主机器人切矢边界条件,确定主机器人的B样条关节轨迹;
S4、同步骤S1,采用累积弦长法,得到从机器人k次B样条轨迹曲线的节点向量;
S6、同步骤S3,结合从机器人切矢边界条件,确定从机器人的B样条关节轨迹;
S7、将主-从机器人运动学约束转化为B样条控制点的约束;
S8、遍历从机器人B样条轨迹的时间区间内的最大装配位置误差,若小于预设的精度,则执行从机器人的B样条曲线插值功能;若大于预设精度,则从机器人根据运动学约束随从主机器人运动。
步骤5、执行双机器人深度搜索避障策略模块;
S1、设置主-从机器人各自的起始点与目标点,令num1=num2=0;
S2、根据任务空间中机器人与障碍物存在的位置关系人为确定大致的避障中间点位置,num1=num1+1;
S3、在避障中间点附近指定的搜索范围内进行启发式搜索,确定避障中间点,num2=num2+1;
S4、分别对主-从机器人的起始点、中间点和目标点进行B样条插值,生成主-从机器人的关节轨迹;
S5、通过避障检测方法判断整个工作时间内机器人和障碍物、双机器人之间是否发生碰撞,若否,则输出主-从机器人的无碰关节轨迹,fob=0终止算法;若是,则执行步骤S6;
S6、判断num2是否小于num2max,若是,则执行步骤S3,若否,则执行步骤S7;
S7、判断num1是否小于num1max,若是,则执行步骤S2,若否,则没有成功搜索到主-从机器人的无碰撞关节轨迹,fob=+∞,算法结束;
其中num1max为人为确定避障中间点位置的尝试次数,num2max为在避障中间点附近范围内随机确定避障中间点的最大允许次数,具体的实现算法流程参考图5。
步骤6、确定协调装配轨迹的优化目标和约束条件,结合量子粒子群算法和避障策略进行求解,即多目标规划模块;
所述结合量子粒子群算法和避障策略进行求解方法涉及到算法,可参考图6,具体的算法执行步骤如下:
S1、初始化设置:在给定范围内对所有粒子的位置进行随机赋值;
S2、判断是否达到最大迭代次数,若是,则跳出算法;若否,则至S3;
S3、执行双机器人协调装配轨迹规划方法所构造的主-从B样条插值轨迹,即步骤4;
S4、执行双机器人深度搜索策略算法,即步骤5;
S5、判断碰撞检测指标fob是否为0,;若是,则令粒子群算法适应度值为无穷大,至S9;若否,至S6;
S8、用公式FG=ω1.ft+ω2.fj+ω3.fob计算适应度函数值;
S9、用公式xid(t+1)=pbestid(t)±α(t)|mbestid(t)-xid(t)|ln[1/uid(t)]更新粒子位置,至S2;
通过算法的循环迭代,将会得到一条速度、加速度、加加速度连续,执行时间较短,关节轨迹平滑的关节轨迹。
步骤7、将多目标规划模块得到的关节空间信息发送给主-从机器人的驱动器。
Claims (5)
1.一种基于B样条曲线的双机器人协调装配的轨迹规划方法,其特征在于以下几个步骤:
步骤1、通过D-H方法分别建立主-从机器人各个连杆坐标系,通过D-H参数表和坐标系变换矩阵推导出运动学正解,然后通过PIPER法推导出运动学逆解;
步骤2、通过示教获取离散装配路径点,并通过步骤1的双机器人运动学逆解公式求任务点姿态;
步骤3、建立主-从机器人协调***的坐标系,通过坐标转换的方法得到主-从机器人的末端坐标系相对于自身基座标系的变换矩阵;
步骤4、以装配体A的轴线为Z1轴,装配体A在轴线方向装配路径段的起始点O1为原点建立装配体A的坐标系,以无装配误差时装配体B的配合面轴线为Z2轴,装配体B在轴线方向装配路径段的起始点O2为原点建立装配体B的坐标系,E为装配体A的原点O1在装配体B的Z2轴上的垂足;F为装配体B的原点O2在装配体A的Z1轴上的垂足,在装配路径段,选取两垂线距离(O1,E)至(O2,F)距离的最大值定义为装配位置误差,将装配体A和装配体B之间姿态变换的欧拉角与理想欧拉角的差值取绝对值求和定义为装配体姿态欧拉角误差;
步骤5、采用B样条曲线插值方法分别对主-从机器人协调装配关节位置点进行插值,其步骤为首先利用德布尔k次非均匀有理B样条递推公式及r阶导矢变换公式构造主-从机器人B样条关节轨迹随时间变化的关节位移、速度、加速度和加加速度曲线,然后获得双机器人协调装配关节轨迹,并对装配误差进行补偿;
步骤6、执行双机器人深度搜索避障策略;
步骤7、以运行时间评价指标、关节平均脉动和碰撞评价指标为协调装配轨迹的优化目标,利用量子粒子群算法和避障策略进行求解;
步骤8、将多目标优化得到的关节空间信息发送给主-从机器人的驱动器。
2.根据权利要求1中所述的基于B样条曲线的双机器人协调装配的轨迹规划方法,其特征在于:所述步骤6中,机器人与障碍物之间的避障策略如下:在机器人的连杆轴线上每隔一定间距设置一个干涉检测点,并将干涉检测点各个区间段连杆的最大包络半径叠加到障碍物的厚度上,从而将机器人与障碍物之间的碰撞转化为由一系列干涉检测点与扩充后的障碍物间的位置判定关系。
3.根据权利要求1中所述的基于B样条曲线的双机器人协调装配的轨迹规划方法,其特征在于:所述步骤6中,机器人与机器人之间的避障策略如下:将机器人简化为由连杆轴组成的直线段,从而将机器人与机器人之间的碰撞简化为判断整个工作空间内空间中两条直线段的最短距离与相应两条连杆的最大包络半径之和的大小判定关系。
4.根据权利要求1中所述的基于B样条曲线的双机器人协调装配的轨迹规划方法,其特征在于:所述步骤6中,双机器人深度搜索避障策略如下步骤:
S1、设置主-从机器人各自的起始点与目标点,令num1=num2=0;
S2、根据任务空间中机器人与障碍物存在的位置关系人为确定大致的避障中间点位置,num1=num1+1;
S3、在避障中间点附近指定的搜索范围内进行启发式搜索,确定避障中间点,num2=num2+1;
S4、分别对主-从机器人的起始点、中间点和目标点进行B样条插值,生成主-从机器人的关节轨迹;
S5、通过避障检测方法判断整个工作时间内机器人和障碍物、双机器人之间是否发生碰撞,若否,则输出主-从机器人的无碰关节轨迹,终止算法;若是,则执行步骤S6;
S6、判断num2是否小于num2max,若是,则执行步骤S3,若否,则执行步骤S7;
S7、判断num1是否小于num1max,若是,则执行步骤S2,若否,则没有成功搜索到主-从机器人的无碰撞关节轨迹,算法结束;
其中num1max为人为确定避障中间点位置的尝试次数,num2max为在避障中间点附近范围内随机确定避障中间点的最大允许次数。
5.根据权利要求1中所述的基于B样条曲线的双机器人协调装配的轨迹规划方法,其特征在于:所述步骤5中,对装配误差进行补偿的方法为:在协调装配路径段,主机器人始终按照B样条插值关节轨迹运动;遍历从机器人最大装配误差,若小于预设精度,从机器人按照从机器人的B样条插值关节轨迹运动;若大于预设精度,则从机器人按照运动学约束随从主机器人运动,按照位置误差和姿态误差对轨迹进行补偿,以免损坏装配体,并更新从机器人装配误差最大值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610965825.6A CN106503373B (zh) | 2016-11-04 | 2016-11-04 | 一种基于b样条曲线的双机器人协调装配轨迹规划方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610965825.6A CN106503373B (zh) | 2016-11-04 | 2016-11-04 | 一种基于b样条曲线的双机器人协调装配轨迹规划方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106503373A CN106503373A (zh) | 2017-03-15 |
CN106503373B true CN106503373B (zh) | 2020-07-28 |
Family
ID=58321896
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610965825.6A Expired - Fee Related CN106503373B (zh) | 2016-11-04 | 2016-11-04 | 一种基于b样条曲线的双机器人协调装配轨迹规划方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106503373B (zh) |
Families Citing this family (27)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2019009350A1 (ja) * | 2017-07-05 | 2019-01-10 | オムロン株式会社 | 経路出力方法、経路出力システム及び経路出力プログラム |
DE102017213651A1 (de) * | 2017-08-07 | 2019-02-07 | Robert Bosch Gmbh | Handhabungseinrichtung mit einem Roboter sowie Verfahren und Computerprogramm |
CN108009680A (zh) * | 2017-11-30 | 2018-05-08 | 航天科工智能机器人有限责任公司 | 基于多目标粒子群算法的仿人机器人步态规划方法 |
CN108009012B (zh) * | 2017-12-14 | 2021-12-14 | 中南大学 | 一种基于任务模型的多智能体动态任务分配方法 |
US10261491B1 (en) * | 2018-03-05 | 2019-04-16 | Rockwell Automation Technologies, Inc. | Method and apparatus for coordinating external device motion with motion of independent movers in a linear drive system |
CN108549423B (zh) * | 2018-04-12 | 2020-11-27 | 江南大学 | 一种加速度上限可变的差分驱动移动机器人速度插值方法 |
CN108436915A (zh) * | 2018-04-17 | 2018-08-24 | 上海达野智能科技有限公司 | 双机器人运动控制方法 |
CN111208807A (zh) * | 2018-11-06 | 2020-05-29 | 苏州艾吉威机器人有限公司 | 一种基于b样条曲线的agv运动控制方法 |
CN109298682A (zh) * | 2018-11-16 | 2019-02-01 | 重庆邮电大学 | 一种基于抛物线插值算法的双机器人协同轨迹优化方法 |
CN109955250A (zh) * | 2019-01-21 | 2019-07-02 | 中国船舶重工集团公司第七一六研究所 | 应用于工业机器人的追踪与实时避障反应规划算法 |
CN109676610B (zh) * | 2019-01-25 | 2021-07-06 | 温州大学 | 一种断路器装配机器人及其实现工作轨迹优化的方法 |
CN110103220B (zh) * | 2019-05-20 | 2021-11-26 | 华南理工大学 | 机器人高速高精度运动轨迹规划方法、装置、设备及介质 |
CN110297490B (zh) * | 2019-06-17 | 2022-06-07 | 西北工业大学 | 基于强化学习算法的异构模块化机器人自重构规划方法 |
CN110455290B (zh) * | 2019-07-17 | 2020-12-08 | 爱克斯维智能科技(苏州)有限公司 | 一种智能液压挖掘机的最优轨迹规划方法 |
CN110757455B (zh) * | 2019-10-25 | 2022-11-22 | 上海新时达机器人有限公司 | 一种速度规划方法 |
CN111230267B (zh) * | 2020-02-11 | 2021-06-22 | 长春华翔轿车***有限责任公司 | 一种双工位协作混焊焊接生产线 |
CN112549025B (zh) * | 2020-11-27 | 2022-06-21 | 北京工业大学 | 一种基于融合人体运动学约束的协调图双臂协作控制方法 |
CN113232021B (zh) * | 2021-05-19 | 2022-09-06 | 中国科学院自动化研究所苏州研究院 | 一种机械臂抓取路径碰撞检测方法 |
CN113246131B (zh) * | 2021-05-27 | 2022-10-28 | 广东智源机器人科技有限公司 | 运动捕捉方法、装置、电子设备以及机械臂控制*** |
CN113561994B (zh) * | 2021-08-13 | 2022-04-05 | 北京三快在线科技有限公司 | 轨迹规划方法、装置、存储介质及电子设备 |
CN114019912B (zh) * | 2021-10-15 | 2024-02-27 | 上海电机学院 | 一种群体机器人运动规划控制方法及*** |
CN114030222B (zh) * | 2021-11-04 | 2023-12-01 | 温州大学 | 可视化冲压线多机协调轨迹优化方法及*** |
CN114407019B (zh) * | 2022-02-21 | 2024-04-05 | 中国科学院宁波材料技术与工程研究所 | 一种工业机器人关节空间的无碰撞轨迹规划方法 |
CN114625148B (zh) * | 2022-03-23 | 2024-06-28 | 中国人民解放军空军军医大学 | 包含模型约束和实时误差***的消毒机器人路径规划方法 |
CN114851207B (zh) * | 2022-06-09 | 2023-08-25 | 珠海格力电器股份有限公司 | 一种双机器人的防碰撞控制方法、***和多机器人*** |
CN116728437B (zh) * | 2023-08-08 | 2023-10-13 | 江苏集萃智能制造技术研究所有限公司 | 基于欧氏空间距离的康养机器人轨迹采样滤波方法及*** |
CN117681213B (zh) * | 2024-02-04 | 2024-05-31 | 之江实验室 | 一种工业机器人运动轨迹规划及评价方法和装置 |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101966617A (zh) * | 2010-08-30 | 2011-02-09 | 东南大学 | 一种用于焊接机器人连续运动的视图模型简约描述方法 |
CN102646148A (zh) * | 2012-04-25 | 2012-08-22 | 浙江大学 | 仿人机器人机械臂避碰运动轨迹规划方法 |
CN102662350A (zh) * | 2012-05-31 | 2012-09-12 | 东南大学 | 主从式多机器人协作***的轨迹示教与规划方法 |
CN103235862A (zh) * | 2013-05-10 | 2013-08-07 | 北京理工大学 | 选择拆卸序列规划方法及装置 |
CN103235513A (zh) * | 2013-04-24 | 2013-08-07 | 武汉科技大学 | 一种基于遗传算法的移动机械臂轨迹规划优化方法 |
CN104215206A (zh) * | 2014-09-28 | 2014-12-17 | 东南大学 | 一种双机器人协作***的基坐标标定方法 |
CN105479490A (zh) * | 2015-12-24 | 2016-04-13 | 华中科技大学 | 一种双机器人实时动态避障装置及其避障方法 |
CN105751196A (zh) * | 2016-04-12 | 2016-07-13 | 华南理工大学 | 一种基于主从式的工业机器人协作作业方法 |
CN105751245A (zh) * | 2016-03-30 | 2016-07-13 | 广东工业大学 | 一种用于标定多机器人***基坐标系的方法及其设备 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP5774361B2 (ja) * | 2011-04-28 | 2015-09-09 | 本田技研工業株式会社 | 軌道計画方法、軌道計画システム及び軌道計画・制御システム |
US8768622B2 (en) * | 2012-09-14 | 2014-07-01 | The United States Of America, As Represented By The Secretary Of The Navy | System and method for maneuver plan for satellites flying in proximity using apocentral coordinate system |
-
2016
- 2016-11-04 CN CN201610965825.6A patent/CN106503373B/zh not_active Expired - Fee Related
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101966617A (zh) * | 2010-08-30 | 2011-02-09 | 东南大学 | 一种用于焊接机器人连续运动的视图模型简约描述方法 |
CN102646148A (zh) * | 2012-04-25 | 2012-08-22 | 浙江大学 | 仿人机器人机械臂避碰运动轨迹规划方法 |
CN102662350A (zh) * | 2012-05-31 | 2012-09-12 | 东南大学 | 主从式多机器人协作***的轨迹示教与规划方法 |
CN103235513A (zh) * | 2013-04-24 | 2013-08-07 | 武汉科技大学 | 一种基于遗传算法的移动机械臂轨迹规划优化方法 |
CN103235862A (zh) * | 2013-05-10 | 2013-08-07 | 北京理工大学 | 选择拆卸序列规划方法及装置 |
CN104215206A (zh) * | 2014-09-28 | 2014-12-17 | 东南大学 | 一种双机器人协作***的基坐标标定方法 |
CN105479490A (zh) * | 2015-12-24 | 2016-04-13 | 华中科技大学 | 一种双机器人实时动态避障装置及其避障方法 |
CN105751245A (zh) * | 2016-03-30 | 2016-07-13 | 广东工业大学 | 一种用于标定多机器人***基坐标系的方法及其设备 |
CN105751196A (zh) * | 2016-04-12 | 2016-07-13 | 华南理工大学 | 一种基于主从式的工业机器人协作作业方法 |
Non-Patent Citations (3)
Title |
---|
Cooperative task planning of multi-robot systems with temporal constraints;Feng-Li Lian 等;《2003 IEEE International Conference on Robotics and Automation》;20030919;第2504-2509页 * |
双机器人协调操作运动规划及控制研究;林风云;《万方数据库-学位论文库》;20060918;全文 * |
基于遗传算法的B样条曲线优化在机器人轨迹规划中应用;陈丹;《万方数据库-学位论文库》;20071226;全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN106503373A (zh) | 2017-03-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106503373B (zh) | 一种基于b样条曲线的双机器人协调装配轨迹规划方法 | |
KR101732902B1 (ko) | 로봇의 경로 계획 장치 및 그 방법 | |
US11813753B2 (en) | Collision avoidance motion planning method for industrial robot | |
CN111300425B (zh) | 一种超冗余机械臂末端轨迹运动规划方法 | |
US8600554B2 (en) | System and method for robot trajectory generation with continuous accelerations | |
KR101667031B1 (ko) | 로봇의 경로 계획 장치 및 그 방법 | |
CN105183009B (zh) | 一种冗余机械臂轨迹控制方法 | |
CN107966907B (zh) | 一种应用于冗余度机械臂的障碍物躲避解决方法 | |
CN109397292B (zh) | 一种基于解析解的7自由度机械臂控制方法与*** | |
CN110609547A (zh) | 一种基于可视图引导的移动机器人规划方法 | |
CN110561419B (zh) | 臂型线约束柔性机器人轨迹规划方法及装置 | |
CN105511266A (zh) | 基于引力搜索粒子群算法的Delta机器人轨迹规划方法 | |
CN112809665B (zh) | 一种基于改进rrt算法的机械臂运动规划方法 | |
CN107538493A (zh) | 一种避免线缆干涉约束的工业机器人轨迹规划方法 | |
CN110231821B (zh) | 多机器人编队的改进自适应零空间行为融合方法 | |
CN112497216B (zh) | 一种基于深度学习的工业机器人位姿精度补偿方法 | |
He et al. | Immune optimization based multi-objective six-DOF trajectory planning for industrial robot manipulators | |
CN113325799A (zh) | 一种曲面工件的点焊机器人操作空间平滑路径规划方法 | |
Hou et al. | Kinematics analysis and self-collision detection of Truss type multi-robot cooperative welding platform | |
Shao et al. | Rrt-goalbias and path smoothing based motion planning of mobile manipulators with obstacle avoidance | |
Furukawa | Time-subminimal trajectory planning for discrete non-linear systems | |
US20210323154A1 (en) | Disassembly based assembly planning | |
Riboli et al. | Collision-free and smooth motion planning of dual-arm Cartesian robot based on B-spline representation | |
CN111515954B (zh) | 一种机械臂高质量运动路径生成方法 | |
Li et al. | Stiffness-maximum trajectory planning of a hybrid kinematic-redundant robot machine |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | 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: 20200728 |