CN102279563A - 非标数字装备操作机参数化运动学模型库建立方法 - Google Patents
非标数字装备操作机参数化运动学模型库建立方法 Download PDFInfo
- Publication number
- CN102279563A CN102279563A CN2010101998337A CN201010199833A CN102279563A CN 102279563 A CN102279563 A CN 102279563A CN 2010101998337 A CN2010101998337 A CN 2010101998337A CN 201010199833 A CN201010199833 A CN 201010199833A CN 102279563 A CN102279563 A CN 102279563A
- Authority
- CN
- China
- Prior art keywords
- alpha
- digital equipment
- kinematics
- kinematics model
- nonstandard
- 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.)
- Pending
Links
Images
Landscapes
- Numerical Control (AREA)
Abstract
一种非标数字装备运动学参数化通用模型库建立方法。目前的非标数字装备控制器多是专用控制器,只能是单一控制器对于单一运动学模型的控制器,一旦控制结构改变就要重新设计控制器,不具通用性、互换性和可重构能力。本发明提出一种非标数字装备运动学模型库建立方法,参数化运动学模型的建立,不是针对某一台非标数字装备操作机而是针对某一类非标数字装备操作机,建立的一种抽象操作机运动学模型。参数化输入包括自由度的参数化,连杆的参数化,关节运动副的参数化等。以Paul等人提出的代数法求解操作机运动学反解,并将运动正反解编写成函数库。用户对模型的选择与输入的D-H参数作为参数传递给函数,使运动学模型与特定的操作机分离。本发明能够减少重复建模工作,提高编程效率,使非标数字装备的控制器具有通用性、互换性和可重构能力。
Description
【技术领域】:本发明属于先进制造技术和运动控制技术领域,具体涉及一种非标数字装备运动学参数化通用模型库建立方法。
【背景技术】:非标数字装备操作机是由若干连杆通过转动或移动关节串联而成,可看成一个开式运动链。开链的一端固定在基座上,另一端为末端执行器,用以操作物体,完成各种作业。操作机运动学研究的是各连杆之间运动的几何关系。操作机运动学主要研究两类问题:①正向运动学,即根据各关节变量的值和连杆的几何参数计算末端执行器相对工作站的位姿;②反向运动学,即根据末端执行器的位姿计算相应的关节变量。非标数字装备控制器是非标数字装备***中最关键、最重要的部分,非标数字装备运动学模型是指非标数字装备操作机运动学正解和逆解,是非标数字装备控制器的基础,直接决定了非标数字装备控制器的应用对象及其互换性。由于非标数字装备操作机结构形式的不同使得其运动学模型呈现多样化。现有的非标数字装备操作机控制器多是专用控制器,为封闭式结构,只能是适用于单一运动学模型的单一控制器,一旦控制结构改变就要重新设计控制器,不具通用性、互换性和可重构能力。
【发明内容】:本发明的目的是克服现有非标数字装备控制器的不足,提供一种非标数字装备参数化运动学模型库建立方法,实现非标数字装备控制器的通用化,使得如自动化专机、生产线、立体仓库、包装机、涂胶机、自动焊机、检测设备、2-3自由度搬运机械手等非标数字装备的控制器具有通用性、互换性和可重构能力。
本发明建立了非标数字装备操作机参数化运动学模型。该模型具备了参数化的特征,满足了非标数字装备操作机控制器具备通用性的要求。该模型库的建立方法是:
(1)构建非标数字装备操作机运动学模型;
(2)将运动学模型建立成一一对应的运动学模型库函数;
(3)自由度、连杆、关节运动副的结构参数可由用户输入并存储在计算机上,作为运动学模型库函数的参数由控制器调用,进而实现非标数字装备操作机的轨迹规划,完成非标数字装备操作机的运动控制。
本发明提供的参数化运动学模型的建立,不是针对某一台非标数字装备操作机,而是针对某一类非标数字装备操作机建立的一种抽象操作机运动学模型。参数输入包括自由度、连杆、关节运动副的参数。在该抽象运动学模型中,自由度数既可能为1个,也可能是2个或3个,任何一个运动关节既可能是移动关节,也可能是旋转关节。用户可以选择操作机运动学模型,并进行自由度、连杆、关节运动副的结构参数输入。如3自由度运动学模型可以选择的结构类型有8种,分别为RRR、RRP、RPR、RPP、PPR、PPP、PRR、PRP(P-表示移动关节,R-表示旋转关节)。
本发明以Denavit-Hartenberg法建立了非标数字装备操作机参数化运动学模型,以Paul等人提出的代数法求解操作机运动学反解,并将运动正反解编写成运动学模型库函数。用户对模型的选择与输入的D-H参数作为参数传递给模型库函数,使运动学模型与特定的操作机分离,能够减少重复建模工作,提高编程效率。
本发明的优点和积极效果:
本发明所提非标数字装备参数化运动学模型库建模方法,实现了运动学模型的通用化。本发明具有以下技术效果:
(1)本发明提出的方法使得非标数字装备控制器具有通用性,根据非标数字装备操作机形式的不同,只须在人机交互界面选择对应的运动学模型并输入D-H参数即可。
(2)当非标数字装备操作机某关节损坏时,控制器可自动识别该关节,并自动选择运动学模型完成任务,具有一定的容错能力。
(3)当非标数字装备操作机结构改变时,只须在选择对应的运动学模型并输入D-H参数即可,应用本方法的控制器拥有可重构能力。
(4)本发明可用于非标数字装备的控制,也可用于其它运动控制***。
【附图说明】:
图1是连杆的长度及扭角;
图2是关节变量及偏置量;
图3是连杆坐标系的设定;
图4是非标数字装备操作机控制器结构示意图。
【具体实施方式】:
实施例1:
下面结合具体实施例对本发明作进一步的说明:
为了研究操作机各连杆之间的几何关系,可在每个连杆上固结一个坐标系,然后描述这些坐标系之间的关系。用一4×4的齐次变换矩阵描述相邻两连杆的空间关系,从而推导出“手爪坐标系”相对于“参考系”的等价齐次变换矩阵,建立操作机的运动学方程。
1、连杆参数
1).连杆参数
如图1所示,与操作机运动学有关的几何量只有两个:
(1)连杆的长度ai——两关节转轴之间的最短距离,即两轴线之间公共垂线之长度。当两轴线相交于一点时,ai=0;当两轴线平行时,有无穷多相等的公共垂线长度;两轴线交错时,只有唯一的公共垂线。所以这样定义连杆长度是唯一的。但是机座及末端连杆,由于它们只有一个关节,特别的规定它们的长度为0,即a0=an=0。对于一端为旋转关节另一端为移动关节的连杆,其长度也规定为零。
(2)连杆的扭角αi——将同一连杆的任一轴线向另一轴线移动,使之相交,则此二直线决定一个与连杆的长度ai垂直的平面,此二直线的平面交角就是该连杆的扭角αi。显然机座杆件和末端杆件的扭角为0。
ai及αi的下标i=0,1,2,…,n,是连杆的编号。编号顺序如下:机座i=0,与机座相连的连杆i=1,与第1号连杆相连接的连杆i=2,……。比i号连杆编号大或小的称i号连杆的上或下编号连杆。
除机座及末端连杆外,其余每个连杆均有两个关节,i号连杆与i-1号连杆相连接的关节称为i号连杆的下关节,另一关节称为i号连杆的上关节。下关节编号为i,上关节编号为i+1。
2).关节变量及偏置量
关节变量是指两相邻连杆相对位置的变化量,一般用广义坐标qi表示。当两连杆以旋转关节相连时,qi为转角θi,当两连杆以移动关节相连时,qi是两者沿关节轴线的相对线位移di。
如图2所示,连杆i与连杆i-1在关节i处相连。若关节i是旋转关节,则关节变量为θi。由于连杆i及连杆i-1的杆长ai及ai-1一般是不相交的,故两者沿i号关节轴线存在一距离,称为连杆的偏置量。若关节i是移动关节,则关节变量就是di。di和θi都带正负号。di表示ai-1与轴线i的交点到ai与轴线i的交点问的距离,沿轴线i测量;θi表示ai-1与ai之间的夹角,绕轴线i由ai-1到ai测量。
d1和θ1的确定方法如下:
(1)若关节1是转动关节,则θ1是可变的,称为关节变量,规定θ1=0为连杆1的零位。习惯约定d1=0;
(2)若关节1是移动关节,则d1是可变的,称为关节变量,规定d1=0为连杆1的零位。习惯约定θ1=0。
上面的约定方法对于dn和θn同样适用。
2、连杆坐标系
为了确定各连杆之间的相对运动和位姿关系,在每一连杆的上关节固结一个坐标系。与基座(连杆0)固结的坐标系称基坐标系,与连杆1固结的坐标系称坐标系{1},与连杆i固结的坐标系称为坐标系{i}。
1)、用连杆坐标系规定连杆参数
利用连杆坐标系,则连杆参数可以明确地定义为:
ai为从Zi到Zi+1沿Xi测量的距离;
αi为从Zi到Zi+1绕Xi旋转的角度;
di为从Xi-1到Xi沿Zi测量的距离;
θi为从Xi-1到Xi绕Zi旋转的角度。
2)、连杆坐标系建立的步骤
对于给定的操作机,它的各个连杆坐标系建立的步骤如下:
(1)找出并画出各个关节轴线。
(2)找出并画出相邻两轴线i和i+1的公垂线ai,或两轴线的交点。求出公垂线ai与轴线i+1的交点令这交点为坐标系{i}的原点Oi。
(3)规定Zi轴与关节i轴重合。
(4)规定Xi轴与公垂线ai重合,若Zi与Zi+1相交,则规定Xi是Zi-1和Zi所成平面的法线。
(5)按右手法则决定Yi=Zi×Xi。
(6)当最后一个关节变量为零时,规定{n-1}与{n}重合。对于首端坐标系{0},原点和X0的方向可任意选取。但是,所选择的坐标系{0}使连杆参数尽可能为零。
3、A矩阵
A矩阵(Denavit-Hartenberg Matrix)是二相邻连杆坐标系的齐次坐标变换矩阵,它将上编号连杆坐标系向下编号连杆坐标系变换。连杆变换Ai可以看成坐标系{i}经以下四个子变换得到的:
(1)绕Zi-1轴转θi角;
(2)沿Zi-1轴移动di;
(3)沿Xi轴移动ai;
(4)绕Xi轴转αi角。
所以
Ai=Rot(zi-1,θi)Trans(0,0,di)Trans(ai,0,0)Rot(xi,αi) (3-1)
即
若关节i是转动关节,则式中di是常量;若关节i是移动关节,则式中di是关节变量,且θi=0,因坐标系{i}与坐标系{i-1}无相对转动。
以非标数字装备操作机中三自由度RRR类操作机为例说明参数化运动学模型建立的过程,如下:
连杆坐标系的设定如图3所示,相应的连杆参数如表1所示。
表1 RRR类操作机的连杆参数
(1)RRR类操作机运动学正解:
根据表1所列连杆参数,即可得到各个连杆变换矩阵:
其中,参数θ1、θ2、θ3为各关节变量,参数α1、α2、a1、a2、d2为用户输入常量。
RRR类操作机运动学方程为
式中
式中简写符号为ci=cosθi,si=sinθi,i=1,2,3。
已知操作机各连杆参数和关节变量θi的基础上就可求出末端执行器的位姿T3,这就是RRR类操作机运动学正解。
(2)RRR类操作机运动学逆解:
A1×T3=A2×A3 (5)
即
利用三角代换:
式中正负号对应于θ1的两个可能解。
当时,由矩阵方程(6)两端的元素(1,4)和(2,4)分别对应相等,得
c1px+s1py-a1=a2c2 (10)
-s1px+c1py=a2s2 (11)
式(10)与(11)的平方和为
c1px+s1py=(px 2+py 2-a1 2+a2 2)/2a1 (12)
同理,令(px 2+py 2-a1 2+a2 2)/2a1=c得θ1的解为
θ1有两个可能解,选定其中之一后,由式(10)与(11),a2≥=0得
由式(15)与式(16),得
c1nx+s1ny=c23 (18)
cxox+s1oy=-s23 (19)
由式(18)与式(19),得
θ3=θ23-θ2=A tan2(-c1ox-s1oy,c1nx+s1ny]-θ2 (20)
根据给定末端执行器位姿T3的情况,就求解出了其对应的各关节变量θ1、θ2、θ3,即RRR类操作机运动学逆解。
建立非标数字装备操作机的运动学模型后,将其编写为非标数字装备操作机参数化运动学模型库函数,其中包括运动学正解、反解等函数,部分代码示例如下:
int RRDkinematic(double DH[2][4],double RRT[4][4])
{…}; //2R操作机运动学正解
int RR Ikinematic(double RRT[4][4],double RR I[2])
{…}; //2R操作机运动学反解
如图4所示,为非标数字装备操作机控制器结构图。将建立的参数化运动学模型库嵌入到运动控制***中,即可实现非标数字装备控制器的通用化。
Claims (1)
1.一种非标数字装备操作机运动学参数化模型库建立方法,其特征在于具体实现方法是:
(1)构建非标数字装备操作机运动学模型;
(2)将运动学模型建立成一一对应的运动学模型库函数;
(3)自由度、连杆、关节运动副的结构参数可由用户输入并存储在计算机上,作为运动学模型库函数的参数由控制器调用,进而实现非标数字装备操作机的轨迹规划,完成非标数字装备操作机的运动控制。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010101998337A CN102279563A (zh) | 2010-06-13 | 2010-06-13 | 非标数字装备操作机参数化运动学模型库建立方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN2010101998337A CN102279563A (zh) | 2010-06-13 | 2010-06-13 | 非标数字装备操作机参数化运动学模型库建立方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN102279563A true CN102279563A (zh) | 2011-12-14 |
Family
ID=45105062
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN2010101998337A Pending CN102279563A (zh) | 2010-06-13 | 2010-06-13 | 非标数字装备操作机参数化运动学模型库建立方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102279563A (zh) |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103792885A (zh) * | 2014-01-20 | 2014-05-14 | 北京理工大学 | 一种数控弯管加工仿真方法及装置 |
CN105788443A (zh) * | 2016-05-20 | 2016-07-20 | 南通大学 | 基于d-h参数法的机器人运动学教具 |
CN108161943A (zh) * | 2017-12-14 | 2018-06-15 | 西安航天精密机电研究所 | 弹丸取放特种机器人的建模与控制方法 |
CN109238188A (zh) * | 2018-08-29 | 2019-01-18 | 苏州伟信奥图智能科技有限公司 | 一种三维运动机构的参数标定方法 |
CN109902069A (zh) * | 2019-03-04 | 2019-06-18 | 重庆科技学院 | 一种智能数学模型储存***及方法 |
CN110587596A (zh) * | 2019-07-30 | 2019-12-20 | 武汉恒新动力科技有限公司 | 多轴构型装置远程控制方法、装置、终端设备及存储介质 |
CN113836607A (zh) * | 2021-09-24 | 2021-12-24 | 广东工业大学 | 一种立体仓库的参数化仿真设计方法及其仿真设计*** |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1161753A (zh) * | 1994-10-24 | 1997-10-08 | 奥林公司 | 模型预测控制装置与方法 |
CN1224515A (zh) * | 1997-02-21 | 1999-07-28 | 贝克·休斯公司 | 自适应面向对象的优化软件*** |
CN1618044A (zh) * | 2001-12-18 | 2005-05-18 | Mts***公司 | 确定控制***控制参数的方法 |
EP2107439A1 (en) * | 2008-04-04 | 2009-10-07 | Honeywell International Inc. | Method and system for the design and implementation of optimal multivariable model predictive controllers for fast-sampling constrained dynamic systems |
-
2010
- 2010-06-13 CN CN2010101998337A patent/CN102279563A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1161753A (zh) * | 1994-10-24 | 1997-10-08 | 奥林公司 | 模型预测控制装置与方法 |
CN1224515A (zh) * | 1997-02-21 | 1999-07-28 | 贝克·休斯公司 | 自适应面向对象的优化软件*** |
CN1618044A (zh) * | 2001-12-18 | 2005-05-18 | Mts***公司 | 确定控制***控制参数的方法 |
EP2107439A1 (en) * | 2008-04-04 | 2009-10-07 | Honeywell International Inc. | Method and system for the design and implementation of optimal multivariable model predictive controllers for fast-sampling constrained dynamic systems |
Non-Patent Citations (1)
Title |
---|
赵连玉等: "非标数字装备操作机运动学分析与应用", 《自合机床与自动化加工技术》, 30 November 2007 (2007-11-30) * |
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103792885A (zh) * | 2014-01-20 | 2014-05-14 | 北京理工大学 | 一种数控弯管加工仿真方法及装置 |
CN103792885B (zh) * | 2014-01-20 | 2016-07-13 | 北京理工大学 | 一种数控弯管加工仿真方法及装置 |
CN105788443A (zh) * | 2016-05-20 | 2016-07-20 | 南通大学 | 基于d-h参数法的机器人运动学教具 |
CN105788443B (zh) * | 2016-05-20 | 2018-05-08 | 南通大学 | 基于d-h参数法的机器人运动学教具 |
CN108161943A (zh) * | 2017-12-14 | 2018-06-15 | 西安航天精密机电研究所 | 弹丸取放特种机器人的建模与控制方法 |
CN108161943B (zh) * | 2017-12-14 | 2020-05-05 | 西安航天精密机电研究所 | 弹丸取放特种机器人的建模与控制方法 |
CN109238188A (zh) * | 2018-08-29 | 2019-01-18 | 苏州伟信奥图智能科技有限公司 | 一种三维运动机构的参数标定方法 |
CN109902069A (zh) * | 2019-03-04 | 2019-06-18 | 重庆科技学院 | 一种智能数学模型储存***及方法 |
CN110587596A (zh) * | 2019-07-30 | 2019-12-20 | 武汉恒新动力科技有限公司 | 多轴构型装置远程控制方法、装置、终端设备及存储介质 |
CN113836607A (zh) * | 2021-09-24 | 2021-12-24 | 广东工业大学 | 一种立体仓库的参数化仿真设计方法及其仿真设计*** |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102279563A (zh) | 非标数字装备操作机参数化运动学模型库建立方法 | |
CN107589934A (zh) | 一种关节型机械臂逆运动学解析解的求取方法 | |
CN104493816B (zh) | 一种变胞并联码垛机器人机构 | |
CN104965517A (zh) | 一种机器人笛卡尔空间轨迹的规划方法 | |
US8733143B2 (en) | Method of incremental forming with successive wrap surfaces | |
CN201936524U (zh) | 自动加工与装配生产线实训考核设备 | |
CN102581849B (zh) | 一种基于nc代码的工业机器人轨迹规划方法 | |
CN105066831A (zh) | 一种单一或多机器人***协同工作坐标系的标定方法 | |
Pieskä et al. | Simulation and programming experiences of collaborative robots for small-scale manufacturing | |
CN101927495A (zh) | 一种冗余度机械臂重复运动规划方法 | |
CN101804627A (zh) | 一种冗余度机械臂运动规划方法 | |
Irani et al. | Analysis of manipulators using SDRE: A closed loop nonlinear optimal control approach | |
CN105171734A (zh) | 高速四轴码垛机器人及码垛方法 | |
CN111791234A (zh) | 一种狭窄空间内多机器人空间工作位置防撞控制算法 | |
Sotnik et al. | Modern industrial robotics industry | |
CN104589340A (zh) | 一种新型六自由度机械手 | |
CN105127988B (zh) | 一种三轴混联机构工业机器人 | |
CN209579568U (zh) | 一种移动式七自由度智能协作机器人 | |
KR20160003774A (ko) | 다관절 링크 기구의 역운동학 해법, 및 이 역운동학 해법을 이용한 교시 데이터 작성 장치 | |
CN207264177U (zh) | 一种基于伺服***的高频焊管成型设备自动控制装置 | |
CN202825825U (zh) | 一种机器人 | |
Brinker et al. | Optimization of the reconfiguration planning of handling systems based on parallel manipulators with Delta-like architecture | |
Li et al. | Analysis of workspace and kinematics for a tomato harvesting robot | |
CN207172076U (zh) | 一种双机械臂*** | |
CN204800636U (zh) | 一种新型六自由度机械手 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C02 | Deemed withdrawal of patent application after publication (patent law 2001) | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20111214 |