CN111185915A - 一种基于刚度性能的机器人钻削加工系统布局方法 - Google Patents

一种基于刚度性能的机器人钻削加工系统布局方法 Download PDF

Info

Publication number
CN111185915A
CN111185915A CN202010024900.5A CN202010024900A CN111185915A CN 111185915 A CN111185915 A CN 111185915A CN 202010024900 A CN202010024900 A CN 202010024900A CN 111185915 A CN111185915 A CN 111185915A
Authority
CN
China
Prior art keywords
robot
matrix
rigidity
stiffness
layout
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
Application number
CN202010024900.5A
Other languages
English (en)
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.)
Beijing Transpacific Technology Development Ltd
University of Shanghai for Science and Technology
Original Assignee
Beijing Transpacific Technology Development 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 Beijing Transpacific Technology Development Ltd filed Critical Beijing Transpacific Technology Development Ltd
Priority to CN202010024900.5A priority Critical patent/CN111185915A/zh
Publication of CN111185915A publication Critical patent/CN111185915A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • BPERFORMING OPERATIONS; TRANSPORTING
    • B25HAND TOOLS; PORTABLE POWER-DRIVEN TOOLS; MANIPULATORS
    • B25JMANIPULATORS; CHAMBERS PROVIDED WITH MANIPULATION DEVICES
    • B25J11/00Manipulators not otherwise provided for
    • B25J11/005Manipulators for mechanical processing tasks
    • BPERFORMING OPERATIONS; TRANSPORTING
    • B25HAND TOOLS; PORTABLE POWER-DRIVEN TOOLS; MANIPULATORS
    • B25JMANIPULATORS; CHAMBERS PROVIDED WITH MANIPULATION DEVICES
    • B25J9/00Programme-controlled manipulators
    • B25J9/16Programme controls
    • B25J9/1602Programme controls characterised by the control system, structure, architecture
    • B25J9/1605Simulation of manipulator lay-out, design, modelling of manipulator

Landscapes

  • Engineering & Computer Science (AREA)
  • Robotics (AREA)
  • Mechanical Engineering (AREA)
  • Automation & Control Theory (AREA)
  • Manipulator (AREA)
  • Numerical Control (AREA)

Abstract

本发明提供一种基于刚度性能的机器人钻削加工系统布局方法,首先基于机器人雅可比矩阵及关节刚度,建立了机器人末端操作刚度矩阵。通过建立机器人刚度性能指标—瑞利商,计算出每个布局点机器人的瑞利商的值,给机器人刚度矩阵提供了一个标量测度。通过对计算得出的瑞利商值的分析,得到机器人刚度性能的变化规律,并得到机器人钻削系统的布局分布,结合工厂的实际要求,最终得到机器人钻削系统的合理布局,在此布局下,机器人刚度性能最优,将此布局方法应用到工业化加工系统中,显著提升了机器人的加工精度。

Description

一种基于刚度性能的机器人钻削加工系统布局方法
技术领域
本发明属于机器人加工应用技术领域,涉及一种基于刚度性能的机器人钻削加工系统布局方法。
背景技术
随着机器人的位置精度、负载能力的不断提高,机器人技术受到了广大航空企业的关注。自动钻削是航空领域应用最广泛的机器人技术,机器人具有灵活性好,智能化程度高,成本低等优点。在飞机装配领域,钻削工作量巨大,其位置精度与质量十分重要。而用于钻削的机器人多为串联机器人,串联机器人具有空间悬臂梁结构,所以刚度性能低,其末端的误差是各个关节和连杆误差的累积和放大,末端定位精度低。针对刚度的研究现状,通过对机器人加工姿态进行优化能有效解决因加工姿态导致的刚度低问题,机器人加工姿态取决于机器人加工系统的布局以及机器人末端的姿态,然而,布局对刚度性能的影响往往都被忽略了,目前机器人钻削加工系统的布局通常是基于设计者的经验,没有定量的选择标准。因此,提出一种合理的机器人钻削系统布局方法,对于提高机器人的刚度性能,对机器人的加工精度提升,提升工业机器人在加工中的通用性具有重要意义。
发明内容
本发明的目的在于解决机器人钻削系统加工中的不足,提出一种基于刚度性能的机器人钻削加工系统布局方法以提升机器人加工精度。基于机器人雅可比矩阵及关节刚度,建立了机器人末端操作刚度矩阵。通过建立机器人刚度性能指标—瑞利商,计算出每个布局点机器人的瑞利商的值,给机器人刚度矩阵提供了一个标量测度。通过对计算得出的瑞利商值的分析,得到机器人刚度性能的变化规律,并得到机器人钻削系统的布局分布,结合工厂的实际要求,最终得到机器人钻削系统的合理布局,在此布局下,机器人刚度性能最优。
为达到上述目的,本发明采用如下述技术方案:
一种基于刚度性能的机器人钻削加工系统布局方法,具体步骤如下:
步骤1:基于机器人运动学矩阵建立机器人刚度模型:
根据等效转换原理,用等效扭转弹簧代替柔性接头;在该情况下,关节刚度由弹簧常数表示,六自由度机器人关节刚度由对角矩阵表示为:
Kq=diag(k1,k2,k3,k4,k5,k6) (1)
其中,Kq是机器人关节刚度矩阵,k1至k6为机器人的关节刚度;
通过机器人雅可比矩阵J,推出关节刚度矩阵到操作刚度矩阵的映射关系为:
K=J(q)-TKqJ(q)-1 (2)
六自由度机器人末端操作矩阵一般形式如下:
Figure BDA0002362118410000021
其中,K是机器人末端操作刚度6×6矩阵,
Figure BDA0002362118410000022
是机器人工具坐标系到机器人基坐标系的转换矩阵,NP是由机器人末端指向工具坐标系的位置矢量,
Figure BDA0002362118410000023
为矩阵
Figure BDA0002362118410000024
的斜对称矩阵;
步骤2:机器人刚度性能评价指标:
建立代表产生单位变形所需的力的瑞利商作为机器人刚度性能的评价指标,将机器人末端操作刚度矩阵根据力和变形的关系分为四个子矩阵:
Figure BDA0002362118410000025
其中,Kfd是关于力与线位移的刚度矩阵;K是关于力与角位移的刚度矩阵;Knd是力矩与线位移的刚度矩阵;K是力矩与角位移的刚度矩阵;
选择表示力与线位移的刚度矩阵Kfd,写出如下关系式:
Figure BDA0002362118410000026
用矢量的长度来表示刚度矩阵瑞丽商,推导机器人刚度矩阵子矩阵Kfd的瑞丽商表达式,它是机器人末端广义力矢量f的长度平方与末端变形矢量d的长度平方之比:
Figure BDA0002362118410000027
Figure BDA0002362118410000028
则称QK(d)为机器人刚度矩阵的瑞利商,并把瑞利商作为机器人刚度矩阵的评价指标,则:
|f|=Qk(d)×|d| (7)
式(7)表示产生单位变形所需的外力值,瑞利商越大,产生单位变形所需的外力就越大,即机器人抵抗变形的能力更强;QK(d)是d的函数,随d的方向变化而变化,通过
Figure BDA0002362118410000029
的特征值来研究其性质,因此,选择矩阵的最小特征值作为机器人刚度性能指标的评价指标I,根据I的值看出在当前布局下机器人钻削加工系统的刚度性能优劣;
步骤3:机器人钻削加工系统的布局:
通过改变机器人与钻削加工点的位置关系,计算得出每个布局点的瑞利商大小,进而得到刚度的变化规律;瑞利商值增大,则刚度性能增强;瑞利商值减小,则刚度性能减弱,机器人更容易受到外力而产生变形;根据工厂的场地、机器人及钻具的具体要求,选择合理的布局方式。
所述步骤1中,假设所有关节都是刚体,整个传动系统的刚度都集中在关节上,通过机器人关节刚度矩阵和雅可比矩阵,推导出机器人末端操作刚度矩阵。
所述步骤2中,机器人刚度矩阵是刚度性能的表达式,但它属于一种不能直观评价刚度性能好坏的张量,因此,应选择一种刚度性能评价指标,通过该指标定量反映结构的刚度性能。
本发明与现有技术相比较,具有以下实质性特点和优点:
本发明对机器人钻削加工系统的布局确定了定量的标准,将此布局方法应用到工业化加工系统中,显著提升了机器人的加工精度。
附图说明
图1是基于刚度性能的机器人钻削加工系统布局方法流程图。
图2是机器人关节刚度示意图。
图3是机器人钻削加工系统布局结构图。
具体实施方式
以下结合附图及实施例,对本发明进行进一步详细说明。
如图1所示,一种基于刚度性能的机器人钻削加工系统布局方法,具体步骤如下:
步骤1:基于机器人运动学矩阵建立机器人刚度模型:
根据等效转换原理,用等效扭转弹簧代替柔性接头;在该情况下,关节刚度由弹簧常数表示,六自由度机器人关节刚度由对角矩阵表示为:
Kq=diag(k1,k2,k3,k4,k5,k6) (1)
其中,Kq是机器人关节刚度矩阵,k1至k6为机器人的关节刚度;
通过机器人雅可比矩阵J,推出关节刚度矩阵到操作刚度矩阵的映射关系为:
K=J(q)-TKqJ(q)-1 (2)
六自由度机器人末端操作矩阵一般形式如下:
Figure BDA0002362118410000031
其中,K是机器人末端操作刚度6×6矩阵,
Figure BDA0002362118410000032
是机器人工具坐标系到机器人基坐标系的转换矩阵,NP是由机器人末端指向工具坐标系的位置矢量,
Figure BDA0002362118410000033
为矩阵
Figure BDA0002362118410000034
的斜对称矩阵。
步骤2:机器人刚度性能评价指标:
机器人刚度矩阵是刚度性能的表达式,但它属于一种不能直观评价刚度性能好坏的张量,因此,应选择一种刚度性能评价指标,通过该指标可以定量反映结构的刚度性能。
建立代表产生单位变形所需的力的瑞利商作为机器人刚度性能的评价指标;将机器人末端操作刚度矩阵根据力和变形的关系分为四个子矩阵:
Figure BDA0002362118410000041
其中,Kfd是关于力与线位移的刚度矩阵;K是关于力与角位移的刚度矩阵;Knd是力矩与线位移的刚度矩阵;K是力矩与角位移的刚度矩阵;
选择表示力与线位移的刚度矩阵Kfd,写出如下关系式:
Figure BDA0002362118410000042
用矢量的长度来表示刚度矩阵瑞丽商,推导机器人刚度矩阵子矩阵Kfd的瑞丽商表达式,它是机器人末端广义力矢量f的长度平方与末端变形矢量d的长度平方之比:
Figure BDA0002362118410000043
Figure BDA0002362118410000044
则称QK(d)为机器人刚度矩阵的瑞利商,并把瑞利商作为机器人刚度矩阵的评价指标,则:
|f|=Qk(d)×|d| (7)
式(7)表示产生单位变形所需的外力值,瑞利商越大,产生单位变形所需的外力就越大,即机器人抵抗变形的能力更强;QK(d)是d的函数,随d的方向变化而变化,通过
Figure BDA0002362118410000045
的特征值来研究其性质,因此,选择矩阵的最小特征值作为机器人刚度性能指标的评价指标I,根据I的值看出在当前布局下机器人钻削加工系统的刚度性能优劣。
步骤3:机器人钻削加工系统的布局:
在连杆刚性的假定下,机器人末端刚度有关节刚度、机器人姿态和作用力共同决定。其中,机器人的姿态对末端刚度具有决定性影响,一个典型的机器人钻削加工系统需要机器人的布局和姿态的配合,这在整个制孔系统的性能中起着至关重要的作用。因此,本发明将机器人布局与姿态之间的耦合问题进行解耦。在整个机器人的工作空间中找到一个最好的姿态是有一定难度的,但是对于一个固定的位置,机器人的姿态因为约束条件是有限的。因此,必须要先进行合理的布局分布,使机器人保持最优的刚度性能。
通过改变机器人与钻削加工点的位置关系,可以通过计算得出每个布局点的瑞利商大小,进而得到刚度的变化规律。瑞利商值增大,则刚度性能增强;瑞利商值减小,则刚度性能减弱,机器人更容易受到外力而产生变形。根据以上结论,在实际加工中,可根据工厂的场地、机器人及钻具的具体要求,选择合理的布局。
如图2所示,为实施例机器人关节刚度示意图。以六自由度机器人为例,假设所有关节都是刚体,整个传动系统的刚度都集中在关节上。根据等效转换原理,用等效扭转弹簧代替柔性接头。在该情况下,关节刚度由弹簧常数表示,机器人关节刚度由对角矩阵表示。
如图3所示,为实施例机器人钻削加工系统布局结构图。机器人与钻削加工点之间的距离定义为L和相对角θ。通过不断改变距离和相对角度,可以得到每个位置的瑞利商值。在整个工作空间中,最小瑞利商的位置是刚度最弱,更容易受外力而产生变形的位置。

Claims (3)

1.一种基于刚度性能的机器人钻削加工系统布局方法,其特征在于,具体步骤如下:
步骤1:基于机器人运动学矩阵建立机器人刚度模型:
根据等效转换原理,用等效扭转弹簧代替柔性接头;在该情况下,关节刚度由弹簧常数表示,六自由度机器人关节刚度由对角矩阵表示为:
Kq=diag(k1,k2,k3,k4,k5,k6) (1)
其中,Kq是机器人关节刚度矩阵,k1至k6为机器人的关节刚度;
通过机器人雅可比矩阵J,推出关节刚度矩阵到操作刚度矩阵的映射关系为:
K=J(q)-TKqJ(q)-1 (2)
六自由度机器人末端操作矩阵一般形式如下:
Figure FDA0002362118400000011
其中,K是机器人末端操作刚度6×6矩阵,
Figure FDA0002362118400000012
是机器人工具坐标系到机器人基坐标系的转换矩阵,NP是由机器人末端指向工具坐标系的位置矢量,
Figure FDA0002362118400000013
为矩阵
Figure FDA0002362118400000014
的斜对称矩阵;
步骤2:机器人刚度性能评价指标:
建立代表产生单位变形所需的力的瑞利商作为机器人刚度性能的评价指标;将机器人末端操作刚度矩阵根据力和变形的关系分为四个子矩阵:
Figure FDA0002362118400000015
其中,Kfd是关于力与线位移的刚度矩阵;K是关于力与角位移的刚度矩阵;Knd是力矩与线位移的刚度矩阵;K是力矩与角位移的刚度矩阵;
选择表示力与线位移的刚度矩阵Kfd,写出如下关系式:
Figure FDA0002362118400000016
用矢量的长度来表示刚度矩阵瑞丽商,推导机器人刚度矩阵子矩阵Kfd的瑞丽商表达式,它是机器人末端广义力矢量f的长度平方与末端变形矢量d的长度平方之比:
Figure FDA0002362118400000017
Figure FDA0002362118400000018
则称QK(d)为机器人刚度矩阵的瑞利商,并把瑞利商作为机器人刚度矩阵的评价指标,则:
|f|=Qk(d)×|d| (7)
式(7)表示产生单位变形所需的外力值,瑞利商越大,产生单位变形所需的外力就越大,即机器人抵抗变形的能力更强;QK(d)是d的函数,随d的方向变化而变化,通过
Figure FDA0002362118400000021
的特征值来研究其性质,因此,选择矩阵的最小特征值作为机器人刚度性能指标的评价指标I,根据I的值看出在当前布局下机器人钻削加工系统的刚度性能优劣;
步骤3:机器人钻削加工系统的布局:
通过改变机器人与钻削加工点的位置关系,计算得出每个布局点的瑞利商大小,进而得到刚度的变化规律;瑞利商值增大,则刚度性能增强;瑞利商值减小,则刚度性能减弱,机器人更容易受到外力而产生变形;根据工厂的场地、机器人及钻具的具体要求,选择合理的布局方式。
2.根据权利要求1所述的基于刚度性能的机器人钻削加工系统布局方法,其特征在于,所述步骤1中,假设所有关节都是刚体,整个传动系统的刚度都集中在关节上,通过机器人关节刚度矩阵和雅可比矩阵,推导出机器人末端操作刚度矩阵。
3.根据权利要求1所述的基于刚度性能的机器人钻削加工系统布局方法,其特征在于,所述步骤2中,机器人刚度矩阵是刚度性能的表达式,但它属于一种不能直观评价刚度性能好坏的张量,因此,应选择一种刚度性能评价指标,通过该指标定量反映结构的刚度性能。
CN202010024900.5A 2020-01-10 2020-01-10 一种基于刚度性能的机器人钻削加工系统布局方法 Pending CN111185915A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010024900.5A CN111185915A (zh) 2020-01-10 2020-01-10 一种基于刚度性能的机器人钻削加工系统布局方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010024900.5A CN111185915A (zh) 2020-01-10 2020-01-10 一种基于刚度性能的机器人钻削加工系统布局方法

Publications (1)

Publication Number Publication Date
CN111185915A true CN111185915A (zh) 2020-05-22

Family

ID=70684586

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010024900.5A Pending CN111185915A (zh) 2020-01-10 2020-01-10 一种基于刚度性能的机器人钻削加工系统布局方法

Country Status (1)

Country Link
CN (1) CN111185915A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114193436A (zh) * 2021-12-07 2022-03-18 珠海格力智能装备有限公司 一种机器人工作空间优化方法、装置、存储介质和设备

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2917322A1 (fr) * 2007-06-15 2008-12-19 Commissariat Energie Atomique Procede de calibration de la position d'un systeme poly-articule, notamment d'un robot
US20140067117A1 (en) * 2012-08-31 2014-03-06 Honda Motor Co., Ltd. Actuating apparatus
CN107414834A (zh) * 2017-08-10 2017-12-01 上海大学 一种多机器人协同系统静刚度实时性能评价方法
CN107703748A (zh) * 2017-10-09 2018-02-16 东南大学 一种基于偏置板设计的重载机器人静态刚度辨识方法
CN109434829A (zh) * 2018-11-07 2019-03-08 华侨大学 一种立体石雕机器人加工系统的变形预测及补偿方法
CN110161850A (zh) * 2019-04-24 2019-08-23 南京航空航天大学 一种工业机器人变参数刚度辨识与建模方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2917322A1 (fr) * 2007-06-15 2008-12-19 Commissariat Energie Atomique Procede de calibration de la position d'un systeme poly-articule, notamment d'un robot
US20140067117A1 (en) * 2012-08-31 2014-03-06 Honda Motor Co., Ltd. Actuating apparatus
CN107414834A (zh) * 2017-08-10 2017-12-01 上海大学 一种多机器人协同系统静刚度实时性能评价方法
CN107703748A (zh) * 2017-10-09 2018-02-16 东南大学 一种基于偏置板设计的重载机器人静态刚度辨识方法
CN109434829A (zh) * 2018-11-07 2019-03-08 华侨大学 一种立体石雕机器人加工系统的变形预测及补偿方法
CN110161850A (zh) * 2019-04-24 2019-08-23 南京航空航天大学 一种工业机器人变参数刚度辨识与建模方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
曲巍崴等: "机器人加工系统刚度性能优化研究", 《航空学报》 *
赵延治等: "弹性铰平面并联三自由度机器人连续刚度映射研究", 《燕山大学学报》 *
赵铁石等: "空间并联机构连续刚度非线性映射", 《机械工程学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114193436A (zh) * 2021-12-07 2022-03-18 珠海格力智能装备有限公司 一种机器人工作空间优化方法、装置、存储介质和设备
CN114193436B (zh) * 2021-12-07 2023-12-26 珠海格力智能装备有限公司 一种机器人工作空间优化方法、装置、存储介质和设备

Similar Documents

Publication Publication Date Title
Briot et al. Are parallel robots more accurate than serial robots?
CN110193829B (zh) 一种耦合运动学与刚度参数辨识的机器人精度控制方法
CN109304730B (zh) 一种基于激光测距仪的机器人运动学参数标定方法
CN107703748B (zh) 一种基于偏置板设计的重载机器人静态刚度辨识方法
US6690999B2 (en) Manipulator control method
CN112571452B (zh) 一种基于刚度模型的机器人刚度设计方法与终端
CN111055273A (zh) 一种用于机器人的两步误差补偿方法
CN104890013A (zh) 一种基于拉线编码器的工业机器人标定算法
CN115229796B (zh) 一种面向航空大构件铣削的机器人末端位姿跟踪补偿方法
CN104537172B (zh) 一种通过优化夹具布局控制孔组复合位置度误差的方法
CN113927599B (zh) 绝对精度补偿方法及系统、设备和计算机可读存储介质
CN113580148B (zh) 基于等效运动链的并联机器人运动学标定方法
CN113500599B (zh) 一种基于神经网络的绳驱机械臂的轨迹跟踪方法及系统
CN108527368B (zh) 柔性支撑串联工业机器人作业最优初始位姿确定方法
CN111185915A (zh) 一种基于刚度性能的机器人钻削加工系统布局方法
CN114474003A (zh) 一种基于参数辨识的车载建筑机器人误差补偿方法
Lou et al. Improved and modified geometric formulation of POE based kinematic calibration of serial robots
CN111390914A (zh) 一种机器人零位和工具坐标标定方法
Yang et al. Self-calibration of three-legged modular reconfigurable parallel robots based on leg-end distance errors
CN114833834A (zh) 一种基于多源误差建模的工业机器人精度补偿方法
CN109352649B (zh) 一种基于深度学习的机械手控制方法及系统
Xu et al. Design, analysis and optimization of Hex4, a new 2R1T overconstrained parallel manipulator with actuation redundancy
CN114474070A (zh) 一种基于刚柔耦合的建筑机器人误差分析方法及验证方法
CN112091972B (zh) 一种基于刚度性能的多机器人系统姿态与布局方法
Mehta et al. A general approach for optimal kinematic design of 6-DOF parallel manipulators

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20200522