CN105717869B - 工业机器人操作空间路径复合限制求解方法 - Google Patents

工业机器人操作空间路径复合限制求解方法 Download PDF

Info

Publication number
CN105717869B
CN105717869B CN201610146940.0A CN201610146940A CN105717869B CN 105717869 B CN105717869 B CN 105717869B CN 201610146940 A CN201610146940 A CN 201610146940A CN 105717869 B CN105717869 B CN 105717869B
Authority
CN
China
Prior art keywords
parameter
path
operating space
restrictive condition
industrial robot
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
Application number
CN201610146940.0A
Other languages
English (en)
Other versions
CN105717869A (zh
Inventor
韩建欢
韩峰涛
张雷
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Rokae (beijing) Technology Co Ltd
Original Assignee
Rokae (beijing) Technology Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Rokae (beijing) Technology Co Ltd filed Critical Rokae (beijing) Technology Co Ltd
Priority to CN201610146940.0A priority Critical patent/CN105717869B/zh
Publication of CN105717869A publication Critical patent/CN105717869A/zh
Application granted granted Critical
Publication of CN105717869B publication Critical patent/CN105717869B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B19/00Programme-control systems
    • G05B19/02Programme-control systems electric
    • G05B19/18Numerical 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/41Numerical 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 interpolation, e.g. the computation of intermediate points between programmed end points to define the path to be followed and the rate of travel along that path
    • G05B19/4103Digital interpolation
    • GPHYSICS
    • G05CONTROLLING; REGULATING
    • G05BCONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B2219/00Program-control systems
    • G05B2219/30Nc systems
    • G05B2219/34Director, elements to supervisory
    • G05B2219/34169Coarse interpolator, path calculator delivers position, speed, acceleration blocks

Landscapes

  • Engineering & Computer Science (AREA)
  • Computing Systems (AREA)
  • Theoretical Computer Science (AREA)
  • Human Computer Interaction (AREA)
  • Manufacturing & Machinery (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Automation & Control Theory (AREA)
  • Numerical Control (AREA)
  • Manipulator (AREA)

Abstract

本发明提出了一种工业机器人操作空间路径复合限制求解方法,包括:设置工业机器人的操作空间路径;计算各轴的速度限制条件和加速度限制条件,根据各轴的速度限制条件和加速度限制条件计算参数的二阶导数的限制条件;根据刚体动力学模型和旋转部分路径,计算各轴力矩限制条件,并根据力矩限制条件计算参数的一阶导数的限制条件;根据工业机器人的操作空间的六维速度旋量计算操作工件的速度限制条件和加速度限制条件,根据速度限制条件和加速度限制条件计算对参数的限制条件;将对工业机器人的复合限制条件统一制对路径参数的限制条件的计算上。本发明有效地求解关节空间和操作空间的多种限制条件,将复合限制统一的表示出来。

Description

工业机器人操作空间路径复合限制求解方法
技术领域
本发明涉及工业机器人技术领域,特别涉及一种工业机器人操作空间路径复合限制求解方法。
背景技术
工业机器人末端执行器的轨迹规划在工业机器人控制系统中占重要地位。轨迹规划生成伺服驱动器的指令值,从而控制机器人各轴电机产生相应的运动,各轴的运动组成末端执行器的操作空间的运动。若使得末端执行器产生用户期望的运动,如跟踪用户指定的曲线,轨迹规划就必须满足伺服驱动器和电机的输入要求,主要体现为:轨迹规划的输出不能使得各轴电机超速;不能使得各轴电机超出容许的最大力矩;不能使得各轴电机超过伺服驱动器容许的加速度和加加速度。同时生成的操作空间运动也不能超过用户指定的最大速度和最大加速度。因此,机器人的轨迹规划,作为规划问题而言,其约束条件是复合的。
现有的轨迹规划方法通常将操作空间的运动映射到机器人的关节空间,映射方法通常是先采样然后在关节空间应用样条曲线进行插补,用这样的方法将电机和伺服驱动器的约束条件直接体现在关节空间,再根据这些条件进行规划。
然而,现有的这种方法有如下缺点:将操作空间路径在关节空间用样条曲线表示,不可避免会与原来的路径有所差别,特别是当采样点较大时,差别尤为明显。并且,这种方法不能体现出操作空间的约束条件。
专利CN 104742127A公开了一种机器人的控制方法及机器人,有效控制了机械臂末端的运行轨迹,提高了机械臂末端的控制精度。但是,该技术方案仍是基于关节空间进行路径规划,采用插补的形式约束机械臂末端的运行轨迹,这种方式不能体现出操作空间的约束条件。
发明内容
本发明的目的旨在至少解决所述技术缺陷之一。
为此,本发明的目的在于提出一种工业机器人操作空间路径复合限制求解方法,能够避免样条曲线插补造成的路径跟踪误差,有效地求解关节空间和操作空间的多种限制条件,将这些限制条件统一地变换为路径参数s的限制,即将复合限制统一的表示出来。
为了实现上述目的,本发明一方面的实施例提供一种工业机器人操作空间路径复合限制求解方法,包括如下步骤:
步骤S1,设置工业机器人的操作空间路径,所述操作空间路径包括:平动部分路径和旋转部分路径q(s);
步骤S2,根据平动部分路径和旋转部分路径q(s),计算所述工业机器人关节空间的各轴的速度限制条件和加速度限制条件,根据所述各轴的速度限制条件和加速度限制条件计算参数s的二阶导数的限制条件;
步骤S3,根据所述工业机器人的刚体动力学模型和旋转部分路径q(s),计算所述工业机器人的各轴力矩限制条件,并根据所述力矩限制条件计算参数s的一阶导数的限制条件;
步骤S4,根据所述工业机器人的操作空间的六维速度旋量计算操作工件的速度限制条件和加速度限制条件,根据所述速度限制条件和加速度限制条件计算对参数s的限制条件;
步骤S5,根据步骤S2、S3和S4中计算得到对参数s、参数s的一阶导数、参数s的二阶导数的限制条件,将对所述工业机器人的复合限制条件统一至对路径参数s的限制条件的计算上。
进一步,在所述步骤S1中,
所述平动部分路径为三维欧式空间的一条曲线,用参数化几何曲线的形式来表示如下:
所述旋转部分路径q(s)采用参数化曲线的形式,以四元数表示工业机器人的末端执行器的姿态
q(s)=slerp(q0,q1,s),
其中,s为路径参数,为平动部分的起点坐标,为平动部分的终点坐标,q0为旋转部分的起点坐标,q1为旋转部分的终点坐标。
进一步,在所述步骤S2中,
设置所述工业机器人的关节空间和操作空间速度有雅克比矩阵映射关系如下:
其中,v为操作空间的六维速度旋量,可以表示为:v=vd·sd,sd表示参数s的一阶导数,
计算得到各轴的速度限制可表示为:
设置关节空间和操作空间加速度的映射关系:
其中,操作空间的六维加速度旋量,参数化形式表示为:sdd表示路径参数s的二阶导数,
计算关节空间各轴的加速度限制为:
当机器人沿某路径运动时,运动状态为(s,sd),此时机器人可用的沿路径的加速度
sdd,需要满足下述不等式,求解该不等式,得到对参数s的二阶导数的限制条件。
其中,
进一步,在所述步骤S3中,
设置所述工业机器人的刚体动力学模型为:
将参数s带入上式,可改写为参数s的形式,
(M(s)α(s))sdd+(M(s)β(s)+C(s))sd2+G(s)=τ
由于力矩约束条件为:τ≤τmax,则α(s)sdd+β(s)sd2+C(s)≤0;
上述不等式在(sd2,sdd)平面构成了一个可行域,通过一般的线性规划方法求得各轴力矩限制对参数s的一阶导数的约束条件。
进一步,在所述步骤S4中,
设置操作空间的六维速度旋量为:v=vd·sd,每个速度分量的限制可以表示为v≤vmax,计算上述不等式得到操作空间速度分量的约束条件对参数s的一阶导数的限制条件;
计算平动速度大小的约束对参数s的限制条件为:|vdp|·sd≤|v|pmax,旋转速度大小的约束对参数s的限制条件为:|vdr|·sd≤|v|rmax,求解得到操作空间速度约束对路径参数s的限制条件;
设置操作空间的六维加速度旋量的参数化形式如下:每个加速度分量的约束条件可以表示为α≤αmax,求解该不等式以得到操作空间每个方向的加速度约束对参数s的限制条件。
根据本发明实施例的工业机器人操作空间路径复合限制求解方法,结合关节空间和操作空间约束条件对路径参数限制,将机器人轨迹规划所涉及的各种约束条件统一地体现到操作空间的路径上,作为沿路径速度规划的输入,能够避免样条曲线插补造成的路径跟踪误差,有效地求解关节空间和操作空间的多种限制条件,将这些限制条件统一地用变换为路径参数s的限制,即将复合限制统一的表示出来,从而为精确的轨迹跟踪等应用提供了先决条件。
本发明附加的方面和优点将在下面的描述中部分给出,部分将从下面的描述中变得明显,或通过本发明的实践了解到。
附图说明
本发明的上述和/或附加的方面和优点从结合下面附图对实施例的描述中将变得明显和容易理解,其中:
图1为根据本发明实施例的工业机器人操作空间路径复合限制求解方法的流程图;
图2为根据本发明实施例的工业机器人操作空间路径复合限制求解方法的示意图;
图3为根据本发明实施例的生成的操作空间的一条圆弧轨迹的对于参数s的一阶导数的复合限制的示意图;
图4为根据本发明实施例的操作空间的路径示意图。
具体实施方式
下面详细描述本发明的实施例,实施例的示例在附图中示出,其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的,旨在用于解释本发明,而不能理解为对本发明的限制。
本发明提出一种工业机器人操作空间路径复合限制求解方法,该方法使用统一参数对操作空间路径进行参数化,并将机器人复合限制条件统一到此参数上。
如图1所示,本发明实施例的工业机器人操作空间路径复合限制求解方法,包括如下步骤:
步骤S1,设置工业机器人的操作空间路径,操作空间路径包括:平动部分路径和旋转部分路径q(s)。
具体地,机器人的末端执行器可视为刚体,操作空间路径可视为刚体的运动路径。其中,操作空间路径包括:平动部分路径和旋转部分路径q(s),两部分和在一起组成SE(3)。
在本步骤中,平动部分路径为三维欧式空间的一条曲线,用参数化几何曲线的形式来表示。操作空间直线路径的平动部分可以表示为,起点坐标加上沿起点到终点的方向增量的形式,如下所示:
其中,s为从0到1的参数,当s从0变到1时直线从起点变到终点。
旋转部分路径q(s)采用参数化曲线的形式,以四元数表示工业机器人的末端执行器的姿态
q(s)=slerp(q0,q1,s),
其中,平动部分的起点坐标,为平动部分的终点坐标,q0为旋转部分的起点坐标,q1为旋转部分的终点坐标,s为路径参数,从0到1的参数,slerp为标准的四元数线性插补函数,当s从0变到1时末端执行器的姿态从初始姿态变化到终点姿态。所以,可以将工业机器人操作空间路径表示为统一的参数化的形式P(s)。这样,工业机器人的多种限制条件就可以转化为对s的约束条件。
步骤S2,根据平动部分路径和旋转部分路径q(s),计算工业机器人关节空间的各轴的速度限制条件和加速度限制条件,根据各轴的速度限制条件和加速度限制条件计算参数s的二阶导数的限制条件。(算法1)
具体地,参考图2,设置工业机器人的关节空间和操作空间速度有雅克比矩阵映射关系如下:
其中,v为操作空间的六维速度旋量,可以表示为:v=vd·sd,即速度方向(包括平动和旋转)乘以速度大小sd表示参数s的一阶导数,
计算得到各轴的速度限制可表示为:
得到2×dof(自由度)个不等式,通过解这些方程就可以得到sd的取值范围。
设置关节空间和操作空间加速度的映射关系:
其中,操作空间的六维加速度旋量,参数化形式表示为:sdd表示路径参数s的二阶导数,即切向加速度(切向方向乘sdd,sdd表示路径参数s的二阶导数)加上法向加速度(速度方向的导数乘速度平方)。带入得到,
其中,
从而得到,
可以写成,
计算关节空间各轴的加速度限制为:
带入得到2×dof个不等式,这些不等式在(sd2,sdd)平面内构成可行域。在此可行域中可以求得轴的加速度限制造成的操作路径上的对参数s的限制。
当机器人沿某路径运动时,运动状态为(s,sd),此时机器人可用的沿路径的加速度sdd,需要满足下述不等式,求解该不等式,得到对参数s的二阶导数的限制条件。
其中,
步骤S3,根据工业机器人的刚体动力学模型和旋转部分路径q(s),计算工业机器人的各轴力矩限制条件,并根据力矩限制条件计算参数s的一阶导数的限制条件。(算法2)
具体地,参考图2,电机提供的力矩是有限的,轨迹规划必须考虑电机的最大输出力矩,以使得伺服驱动器较好的控制电机,避免出现力矩饱和的情况。
机器人的刚体动力学模型为:
轨迹所需各轴力矩由三个部分组成:分别为惯性项、柯氏力项和重力项。将参数s带入,可改写为参数s的形式,
M(s)(α(s)sdd+β(s)sd2)+C(s)sd2+G(s)=τ
将上式改写成,
(M(s)α(s))sdd+(M(s)β(s)+C(s))sd2+G(s)=τ
力矩约束条件为:τ≤τmax
基于上述力矩约束条件,可得到下式:α(s)sdd+β(s)sd2+C(s)≤0
同样的,此式在(sd2,sdd)平面,构成了一个可行域,通过一般的线性规划方法可以求得各轴力矩限制对参数s的一阶导数的约束条件。
步骤S4,根据工业机器人的操作空间的六维速度旋量计算操作工件的速度限制条件和加速度限制条件,根据速度限制条件和加速度限制条件计算对参数s的限制条件。(算法3)
具体地,参考图2,操作空间的六维速度旋量可以表示为速度方向乘以速度大小,即
v=vd·sd
每个速度分量的限制可以表示为:v≤vmax
通过求解上式代表的不等式组即可得到操作空间速度分量的约束条件对参数s的一阶导数的约束。
在实际中,往往要求限制平动速度或者旋转速度的大小,将上式改写为:
vp=vdp·sd
vr=vdr·sd
分别代表了参数s与平动速度和旋转速度之间的关系。
在本发明的一个实施例中,平动速度大小的约束对参数s的限制为:
|vdp|·sd≤|v|pmax
旋转速度大小的约束对参数s的限制为:
|vdr|·sd≤|v|rmax
解上述不等式即可得到操作空间速度约束对参数s的限制。
设置操作空间的六维加速度旋量可以表示为如下的参数化形式,
每个加速度分量的约束条件可以表示为:α≤αmax
采用与求解力矩约束相同的方法解上式代表的不等式组,可得到操作空间每个方向的加速度约束对参数s的限制。
与速度类似,实际中往往对加速度大小的约束较为关心。
具体地,平动加速度大小的约束可以表示为:
因为加速度为切向分量和法向分量的矢量加和,上式可以写成,
旋转加速度与此类似,表示如下:由此可以计算得到加速度大小的约束对参数s的限制。
步骤S5,根据步骤S2、S3和S4中计算得到对参数s、参数s的一阶导数、参数s的二阶导数的限制条件,将对工业机器人的复合限制条件统一制对路径参数s的限制条件的计算上,从而避免了样条曲线插值带来的误差,由此生成的轨迹具有路径不变的性质。
图3为根据本发明实施例的生成的操作空间的一条圆弧轨迹的对于参数s的一阶导数的复合限制的示意图。其中,曲线A、B和C分别对应参数s的不同的限制条件。其中,A是关节空间各轴加速度对参数s一阶导数的限制,B是关节空间电机力矩对参数s一阶导数的限制,C是关节空间各轴转速对参数s一阶导数的限制。
图4为根据本发明实施例的操作空间的路径示意图。其中,该路径示意图包括位置和姿态的变化。
根据本发明实施例的工业机器人操作空间路径复合限制求解方法,结合关节空间和操作空间约束条件对路径参数限制,将机器人轨迹规划所涉及的各种约束条件统一地体现到操作空间的路径上,作为沿路径速度规划的输入,能够避免样条曲线插补造成的路径跟踪误差,有效地求解关节空间和操作空间的多种限制条件,将这些限制条件统一地用变换为路径参数s的限制,即将复合限制统一的表示出来,从而为精确的轨迹跟踪等应用提供了先决条件。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中,对上述术语的示意性表述不一定指的是相同的实施例或示例。而且,描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例或示例中以合适的方式结合。
尽管上面已经示出和描述了本发明的实施例,可以理解的是,上述实施例是示例性的,不能理解为对本发明的限制,本领域的普通技术人员在不脱离本发明的原理和宗旨的情况下在本发明的范围内可以对上述实施例进行变化、修改、替换和变型。本发明的范围由所附权利要求极其等同限定。

Claims (5)

1.一种工业机器人操作空间路径复合限制求解方法,其特征在于,包括如下步骤:
步骤S1,设置工业机器人的操作空间路径,所述操作空间路径包括:平动部分路径和旋转部分路径q(s);
步骤S2,根据平动部分路径和旋转部分路径q(s),计算所述工业机器人关节空间的各轴的速度限制条件和加速度限制条件,根据所述各轴的速度限制条件和加速度限制条件计算参数s的二阶导数的限制条件;
步骤S3,根据所述工业机器人的刚体动力学模型和旋转部分路径q(s),计算所述工业机器人的各轴力矩限制条件,并根据所述力矩限制条件计算参数s的一阶导数的限制条件;
步骤S4,根据所述工业机器人的操作空间的六维速度旋量计算操作工件的速度限制条件和加速度限制条件,根据所述速度限制条件和加速度限制条件计算对参数s的限制条件;
步骤S5,根据步骤S2、S3和S4中计算得到对参数s、参数s的一阶导数、参数s的二阶导数的限制条件,将对所述工业机器人的复合限制条件统一至对路径参数s的限制条件的计算上。
2.如权利要求1所述的工业机器人操作空间路径复合限制求解方法,其特征在于,在所述步骤S1中,
所述平动部分路径为三维欧式空间的一条曲线,用参数化几何曲线的形式来表示如下:
所述旋转部分路径q(s)采用参数化曲线的形式,以四元数表示工业机器人的末端执行器的姿态
q(s)=slerp(q0,q1,s),
其中,s为路径参数,为平动部分的起点坐标,为平动部分的终点坐标,q0为旋转部分的起点坐标,q1为旋转部分的终点坐标。
3.如权利要求1所述的工业机器人操作空间路径复合限制求解方法,其特征在于,在所述步骤S2中,
设置所述工业机器人的关节空间和操作空间速度有雅克比矩阵映射关系如下:
<mrow> <mi>v</mi> <mo>=</mo> <mi>J</mi> <mover> <mi>q</mi> <mo>&amp;CenterDot;</mo> </mover> </mrow>
其中,v为操作空间的六维速度旋量,可以表示为:v=vd·sd,sd表示参数s的一阶导数,
计算得到各轴的速度限制可表示为:
<mrow> <msubsup> <mover> <mi>q</mi> <mo>&amp;CenterDot;</mo> </mover> <mi>k</mi> <mi>min</mi> </msubsup> <mo>&amp;le;</mo> <msub> <mi>&amp;alpha;</mi> <mi>k</mi> </msub> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> <mi>s</mi> <mi>d</mi> <mo>&amp;le;</mo> <msubsup> <mover> <mi>q</mi> <mo>&amp;CenterDot;</mo> </mover> <mi>k</mi> <mi>max</mi> </msubsup> <mo>,</mo> <mi>k</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2...</mn> <mi>d</mi> <mi>o</mi> <mi>f</mi> </mrow>
设置关节空间和操作空间加速度的映射关系:
其中,操作空间的六维加速度旋量,参数化形式表示为:sdd表示路径参数s的二阶导数,
计算关节空间各轴的加速度限制为:
<mrow> <msubsup> <mover> <mi>q</mi> <mo>&amp;CenterDot;&amp;CenterDot;</mo> </mover> <mi>k</mi> <mi>min</mi> </msubsup> <mo>&amp;le;</mo> <msub> <mover> <mi>q</mi> <mo>&amp;CenterDot;&amp;CenterDot;</mo> </mover> <mi>k</mi> </msub> <mo>&amp;le;</mo> <msubsup> <mover> <mi>q</mi> <mo>&amp;CenterDot;&amp;CenterDot;</mo> </mover> <mi>k</mi> <mi>max</mi> </msubsup> <mo>,</mo> <mi>k</mi> <mo>=</mo> <mn>1</mn> <mo>,</mo> <mn>2...</mn> <mi>d</mi> <mi>o</mi> <mi>f</mi> </mrow>
当机器人沿某路径运动时,运动状态为(s,sd),此时机器人可用的沿路径的加速度sdd,需要满足下述不等式,求解该不等式,得到对参数s的二阶导数的限制条件;
<mrow> <msubsup> <mover> <mi>q</mi> <mo>&amp;CenterDot;&amp;CenterDot;</mo> </mover> <mi>k</mi> <mi>min</mi> </msubsup> <mo>-</mo> <msub> <mi>&amp;beta;</mi> <mi>k</mi> </msub> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> <msup> <mi>sd</mi> <mn>2</mn> </msup> <mo>&amp;le;</mo> <msub> <mi>&amp;alpha;</mi> <mi>k</mi> </msub> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> <mi>s</mi> <mi>d</mi> <mi>d</mi> <mo>&amp;le;</mo> <msubsup> <mover> <mi>q</mi> <mo>&amp;CenterDot;&amp;CenterDot;</mo> </mover> <mi>k</mi> <mi>max</mi> </msubsup> <mo>-</mo> <msub> <mi>&amp;beta;</mi> <mi>k</mi> </msub> <mrow> <mo>(</mo> <mi>s</mi> <mo>)</mo> </mrow> <msup> <mi>sd</mi> <mn>2</mn> </msup> <mo>,</mo> </mrow>
其中,
4.如权利要求1所述的工业机器人操作空间路径复合限制求解方法,其特征在于,在所述步骤S3中,
设置所述工业机器人的刚体动力学模型为:
将参数s带入上式,可改写为参数s的形式,
(M(s)α(s))sdd+(M(s)β(s)+C(s))sd2+G(s)=τ
由于力矩约束条件为:τ≤τmax,则α(s)sdd+β(s)sd2+C(s)≤0;
上述不等式在(sd2,sdd)平面构成了一个可行域,通过一般的线性规划方法求得各轴力矩限制对参数s的一阶导数的约束条件。
5.如权利要求1所述的工业机器人操作空间路径复合限制求解方法,其特征在于,在所述步骤S4中,
设置操作空间的六维速度旋量为:v=vd·sd,每个速度分量的限制可以表示为v≤vmax,计算上述不等式得到操作空间速度分量的约束条件对参数s的一阶导数的限制条件;
计算平动速度大小的约束对参数s的限制条件为:|vdp|·sd≤|v|pmax,旋转速度大小的约束对参数s的限制条件为:|vdr|·sd≤|v|rmax,求解得到操作空间速度约束对路径参数s的限制条件;
设置操作空间的六维加速度旋量的参数化形式如下:每个加速度分量的约束条件可以表示为α≤αmax,求解该不等式以得到操作空间每个方向的加速度约束对参数s的限制条件。
CN201610146940.0A 2016-03-15 2016-03-15 工业机器人操作空间路径复合限制求解方法 Active CN105717869B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610146940.0A CN105717869B (zh) 2016-03-15 2016-03-15 工业机器人操作空间路径复合限制求解方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610146940.0A CN105717869B (zh) 2016-03-15 2016-03-15 工业机器人操作空间路径复合限制求解方法

Publications (2)

Publication Number Publication Date
CN105717869A CN105717869A (zh) 2016-06-29
CN105717869B true CN105717869B (zh) 2018-05-29

Family

ID=56157778

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610146940.0A Active CN105717869B (zh) 2016-03-15 2016-03-15 工业机器人操作空间路径复合限制求解方法

Country Status (1)

Country Link
CN (1) CN105717869B (zh)

Families Citing this family (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108544493B (zh) * 2018-04-11 2019-06-11 广州启帆工业机器人有限公司 一种针对机器人操作空间的时间最优轨迹规划方法和装置
EP3569367B1 (de) * 2018-05-17 2022-08-03 Siemens Aktiengesellschaft Rechnergestütztes ermitteln einer bewegung einer vorrichtung
CN109634100B (zh) * 2018-12-30 2021-11-02 深圳市优必选科技有限公司 仿人机器人行走加速度补偿方法、装置及仿人机器人
CN110757455B (zh) * 2019-10-25 2022-11-22 上海新时达机器人有限公司 一种速度规划方法
CN112338918B (zh) * 2020-10-29 2022-05-27 珠海格力智能装备有限公司 六轴机器人工作范围的确定方法及装置
CN112847335A (zh) * 2020-12-24 2021-05-28 珞石(山东)智能科技有限公司 机器人伺服的滑膜变结构控制方法及装置

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1780679A1 (en) * 2002-02-12 2007-05-02 The University of Tokyo Method for generating a motion of a human type link system
CN102510793A (zh) * 2010-06-04 2012-06-20 中国科学院自动化研究所 加速度连续的机器人轨迹生成系统和方法
JP2013184245A (ja) * 2012-03-07 2013-09-19 Canon Inc ロボット制御装置、ロボット装置、ロボット制御方法、プログラム及び記録媒体
JP2013184233A (ja) * 2012-03-06 2013-09-19 Toyota Motor Corp ロボットおよびその軌道生成方法
CN103970139A (zh) * 2014-05-09 2014-08-06 上海交通大学 一种机器人连续点位运动规划方法及其运动控制器
CN104742127A (zh) * 2015-04-08 2015-07-01 深圳市山龙科技有限公司 一种机器人的控制方法及机器人
CN105159096A (zh) * 2015-10-10 2015-12-16 北京邮电大学 一种基于粒子群算法的冗余度空间机械臂关节力矩优化方法

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1780679A1 (en) * 2002-02-12 2007-05-02 The University of Tokyo Method for generating a motion of a human type link system
CN102510793A (zh) * 2010-06-04 2012-06-20 中国科学院自动化研究所 加速度连续的机器人轨迹生成系统和方法
JP2013184233A (ja) * 2012-03-06 2013-09-19 Toyota Motor Corp ロボットおよびその軌道生成方法
JP2013184245A (ja) * 2012-03-07 2013-09-19 Canon Inc ロボット制御装置、ロボット装置、ロボット制御方法、プログラム及び記録媒体
CN103970139A (zh) * 2014-05-09 2014-08-06 上海交通大学 一种机器人连续点位运动规划方法及其运动控制器
CN104742127A (zh) * 2015-04-08 2015-07-01 深圳市山龙科技有限公司 一种机器人的控制方法及机器人
CN105159096A (zh) * 2015-10-10 2015-12-16 北京邮电大学 一种基于粒子群算法的冗余度空间机械臂关节力矩优化方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
《A Unified Approach for Motion and Force Control of Robot Manipulators:The Operational Space Formulation》;OUSSAMA KHATIB;《IEEE JOURNAL OF ROBOTICS AND AUTOMATION》;19870228(第1期);全文 *
《工业机器人姿态规划及轨迹优化研究》;季晨;《中国优秀硕士学位论文全文数据库 信息科技辑》;20140415(第4期);全文 *

Also Published As

Publication number Publication date
CN105717869A (zh) 2016-06-29

Similar Documents

Publication Publication Date Title
CN105717869B (zh) 工业机器人操作空间路径复合限制求解方法
JP5327722B2 (ja) ロボットの負荷推定装置及び負荷推定方法
CN104589304A (zh) 机器人控制装置以及机器人
WO2001038048A1 (fr) Unite de commande de robot
Lee et al. Development, analysis, and control of series elastic actuator-driven robot leg
US20100198404A1 (en) Control device and method for a manipulator
Lasnier et al. Workspace based force sensorless bilateral control with multi-degree-of-freedom motion systems
CN106426163B (zh) 电动机的控制装置
SE516720C2 (sv) Utrustning för styrning av en industrirobot och förfarande för programmering och/eller justering av robotens rörelse
JP3296585B2 (ja) 位置および力を制御するハイブリッドシステム
Tumari et al. The direct strain feedback with PID control approach for a flexible manipulator: Experimental results
Devie et al. Impedance control using a cascaded loop force control
JP4269150B2 (ja) ロボット制御装置
Miro et al. Modelling an industrial manipulator a case study
JPS61201304A (ja) ロボツトの位置制御法
Kim et al. Joint compliance error compensation for robot manipulator using body frame
Jun et al. Continuous terminal sliding mode control of a 6-DOF wire-driven parallel robot
Filaretov et al. Adaptive force/position control of robot manipulators
Winkler et al. Force-guided motions of a 6-dof industrial robot with a joint space approach
Sharan et al. Simulating the arm movements of a stepper motor controlled pick-and-place robot using the stepper motor model
JPH0760667A (ja) ロボットのウィービング制御装置
Lee et al. Interactive force control of an elastically actuated bi-articular two-link manipulator
Luong et al. Continuous terminal sliding mode control with perturbation estimation for a stewart platform
Imani et al. High-speed contouring enhanced with PH curves
Zhang Research and analysis on the robot trajectory interpolation methods

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant