CN106934108B - 一种并联构型主轴头的动态误差传递特性评价方法 - Google Patents
一种并联构型主轴头的动态误差传递特性评价方法 Download PDFInfo
- Publication number
- CN106934108B CN106934108B CN201710089785.8A CN201710089785A CN106934108B CN 106934108 B CN106934108 B CN 106934108B CN 201710089785 A CN201710089785 A CN 201710089785A CN 106934108 B CN106934108 B CN 106934108B
- Authority
- CN
- China
- Prior art keywords
- parallel configuration
- spindle head
- error
- tool nose
- configuration spindle
- 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
Links
- 238000012546 transfer Methods 0.000 title claims abstract description 44
- 238000000034 method Methods 0.000 title claims abstract description 17
- 239000011159 matrix material Substances 0.000 claims abstract description 33
- 230000005540 biological transmission Effects 0.000 claims abstract description 22
- 230000033001 locomotion Effects 0.000 claims abstract description 12
- 238000013519 translation Methods 0.000 claims description 6
- 238000013507 mapping Methods 0.000 claims description 4
- 238000011156 evaluation Methods 0.000 abstract description 16
- 150000001875 compounds Chemical class 0.000 description 9
- 239000013598 vector Substances 0.000 description 5
- 238000012994 industrial processing Methods 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000003754 machining Methods 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research 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
Landscapes
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Numerical Control (AREA)
Abstract
本发明公开了一种并联构型主轴头的动态误差传递特性评价方法,属于机床技术领域。该方法首先推导并联构型主轴头各驱动轴跟随误差与终端动平台运动误差间的传递关系;进一步基于并联构型主轴头刀尖点位置偏差,建立各驱动轴跟随误差与刀尖点位置偏差间的传递关系,得到表征并联构型主轴头动态误差传递特性的无量纲方阵;基于该无量纲方阵,提出无量纲误差传递指标,该误差传递指标的大小表征各驱动轴跟随误差对刀尖点位置偏差影响的强弱。本发明克服了传统动态误差传递特性评价方法物理意义不明确的问题,通过提出无量纲误差传递指标,能够对并联构型主轴头的动态误差传递特性进行更为有效的评价。
Description
技术领域
本发明属于机床技术领域,特别涉及并联构型主轴头的动态误差传递特性评价问题。
背景技术
运动精度是面向工业加工的并联构型主轴头最重要的指标。相对于传统串联构型主轴头,并联构型主轴头因没有误差累积,在理论上具有更高的运动精度,因此有很多并联构型主轴头方案被相继提出,但大多数并联构型主轴头在实际中动态精度过低,不能满足加工需求。因此需要对并联构型主轴头的精度问题开展更深入的研究,这其中,并联构型主轴头的动态误差传递特性评价是一项核心内容。
在并联构型主轴头的运动过程中,各驱动轴的跟随误差是造成刀尖点位置偏差的主要因素,刀尖点位置偏差进而会对加工精度造成影响,因此需要对并联构型主轴头的动态误差传递特性进行有效评价。传统雅克比矩阵代表了并联构型主轴头各驱动轴跟随误差与终端动平台运动误差之间的映射关系,通常使用雅克比矩阵对并联构型主轴头的动态误差传递特性进行描述和评价,但是对于面向工业加工的并联构型主轴头,由于其终端动平台同时耦合了转动及平动自由度并存在伴生运动,导致雅克比矩阵存在以下两个显著问题:1)量纲问题,雅克比矩阵中元素的物理单位不一致;2)维数问题,雅克比矩阵不是一个方阵。由于量纲及维数问题,传统雅克比矩阵不再具有明确的物理意义,同时基于雅克比矩阵提出的评价指标不能够很好的评价并联构型主轴头的动态误差传递特性,如雅克比矩阵的最大奇异值、最小奇异值和条件数。
目前针对并联构型主轴头动态误差传递特性的评价方法较少,已有的评价方法及指标不能精确评价并联构型主轴头各驱动轴跟随误差对刀尖点位置偏差影响的强弱。因此为更好地分析、评价并联构型主轴头的动态误差传递特性,需要基于刀尖点位置偏差,建立新的误差传递矩阵,并进一步提出新的评价指标。
目前还没有普遍适用于并联构型主轴头的动态误差传递特性评价方法,针对并联构型主轴头的特点,提出一种广泛有效的动态误差传递特性评价方法对改善并联构型主轴头的动态精度具有重要意义。
发明内容
本发明的目的是为克服已有评价方法的不足之处,提出了一种并联构型主轴头的动态误差传递特性评价方法,用来解决并联构型主轴头的动态误差传递特性评价问题,从而更精确地评价并联构型主轴头各驱动轴跟随误差对刀尖点位置偏差的影响。
本发明的技术方案如下:
一种并联构型主轴头的动态误差传递特性评价方法,包括如下步骤:
1)建立并联构型主轴头各驱动轴跟随误差与终端动平台运动误差间的传递关系:
dX=Gadq
式中,dX=[dx dy dz dθx dθy dθz]T,dx、dy、dz分别为并联构型主轴头终端动平台沿X方向、Y方向和Z方向的平动误差,这里X、Y、Z为传统笛卡尔坐标系的三个方向;dθx、dθy、dθz分别为并联构型主轴头终端动平台绕X方向、Y方向和Z方向的转动误差;dq=[dq1dq2 dq3]T,dq1、dq2、dq3分别为并联构型主轴头第1驱动轴、第2驱动轴和第3驱动轴的跟随误差;矩阵Ga具有6行3列元素,表征并联构型主轴头各驱动轴跟随误差与终端动平台运动误差间的传递关系;
2)基于并联构型主轴头刀尖点位置偏差,建立并联构型主轴头各驱动轴跟随误差与刀尖点位置偏差间的传递关系:
dP=Tdq
式中,dP=[dPx dPy dPz]T,dPx、dPy、dPz分别为刀尖点沿X方向、Y方向和Z方向的平动误差;T是一个无量纲方阵,具体表示为:
该无量纲方阵用于表征并联构型主轴头的动态误差传递特性;
3)基于步骤2)中的无量纲方阵T,提出无量纲误差传递指标ETIj:
式中,ETIj为并联构型主轴头第j驱动轴的误差传递指标,j=1,2,3;tij为无量纲方阵T相应位置上的元素,i=1,2,3,j=1,2,3;S为无量纲方阵T的最大奇异值;∑||代表对元素取绝对值后进行求和;
4)利用3)中提出的无量纲误差传递指标ETIj对并联构型主轴头动态误差传递特性进行评价,ETIj的大小表征并联构型主轴头第j驱动轴的跟随误差对刀尖点位置偏差影响的强弱,ETIj的值越大,代表并联构型主轴头第j驱动轴的跟随误差对刀尖点位置偏差的影响越大。
本发明的上述方法中,建立并联构型主轴头各驱动轴跟随误差与刀尖点位置偏差间的传递关系的方法包括如下步骤:
(1)建立并联构型主轴头刀尖点位置方程:
式中,Px、Py、Pz分别为刀尖点沿X方向、Y方向和Z方向的位置;x、y、z分别为并联构型主轴头动平台沿X方向、Y方向和Z方向的位置;L为刀具长度;欧拉角α与β用于表示并联构型主轴头刀尖点的位姿;sα=sin(α),sβ=sin(β),cα=cos(α),cβ=cos(β);
(2)对步骤(1)中所建立的并联构型主轴头刀尖点位置方程进行一阶摄动,建立刀尖点位置偏差方程:
(3)建立转角偏差映射关系:
(4)联立以下三个方程:
dX=Gadq
最终建立并联构型主轴头各驱动轴跟随误差与刀尖点位置偏差间的传递关系:
本发明具有以下优点及突出性的技术效果:本发明通过建立各驱动轴跟随误差与刀尖点位置偏差间的传递关系,得到表征并联构型主轴头动态误差传递特性的无量纲方阵;基于该无量纲方阵,提出无量纲误差传递指标,该误差传递指标的大小表征各驱动轴跟随误差对刀尖点位置偏差影响的强弱;该方法克服了传统动态误差传递特性评价方法物理意义不明确的问题,通过提出无量纲误差传递指标,能够对并联构型主轴头的动态误差传递特性进行更为有效的评价;从数值上有效评价并联构型主头各驱动轴跟随误差对刀尖点偏差造成的影响,进而达到动态误差传递特性评价的目的。
附图说明
图1为本发明一种并联构型主轴头的动态误差传递特性评价方法的流程图。
图2为一种典型的并联构型主轴头。
图3为使用本发明得到的并联构型主轴头动态误差传递特性评价结果。
图2中:1-第一滑块;2-第二滑块;3-第三滑块;4-动平台;5-第一杆件;6-第二杆件;7-第三杆件。
具体实施方式
以下结合附图及具体实施例对本发明作进一步详细说明。
图1所示为本发明一种并联构型主轴头的动态误差传递特性评价方法的流程图。图2所示为一钟典型三自由度并联构型主轴头,该主轴头通过第一滑块1、第二滑块2与第三滑块3的运动带动终端动平台4的运动,第一滑块1、第二滑块2与第三滑块3由相应电机驱动,动平台4与第一滑块1间通过第一杆件5进行连接,动平台4与第二滑块2间通过第二杆件6进行连接,动平台4与第三滑块3间通过第三杆件7进行连接;滑块1所在轴为第1驱动轴,滑块2所在轴为第2驱动轴,滑块3所在轴为第3驱动轴。将所提出的一种并联构型主轴头的动态误差传递特性评价方法应用于该并联构型主轴头,具体方法步骤如下:
1)建立并联构型主轴头的误差传递方程:
dt+dθ×ai=dqie3+dθi×li (1)
式中,dt=[dx dy dz]T,dx、dy、dz分别为并联构型主轴头终端动平台沿X方向、Y方向和Z方向的平动误差,这里X、Y、Z为传统笛卡尔坐标系的三个方向;dθ=[dθx dθy dθz]T,dθx、dθy、dθz分别为并联构型主轴头终端动平台绕X方向、Y方向和Z方向的转动误差;dqi为并联构型主轴头第i驱动轴的跟随误差,i=1,2,3;e3=[0 0 1]T;ai为动平台第i端点的位置向量,i=1,2,3;li为第i被动杆件的位置向量,i=1,2,3;
2)用向量li点乘方程(1)的两端,建立并联构型主轴头的驱动方程:
3)用向量u1、u2、u3依次点乘方程(1)的两端,u1、u2、u3分别为第1被动杆件、第2被动杆件和第3被动杆件的转动轴轴向单位方向向量,建立并联构型主轴头的约束方程:
4)联立式(2)、(3)得:
5)由式(4)进一步得到:
dX=Gadq (5)
式中,矩阵Ga具有6行3列元素,表征并联构型主轴头各驱动轴跟随误差与终端动平台运动误差间的传递关系;
6)建立并联构型主轴头刀尖点位置方程:
式中,Px、Py、Pz分别为刀尖点沿X方向、Y方向和Z方向的位置;x、y、z分别为并联构型主轴头动平台沿X方向、Y方向和Z方向的位置;L为刀具长度;欧拉角α与β用于表示并联构型主轴头刀尖点的位姿,sα=sin(α),sβ=sin(β),cα=cos(α),cβ=cos(β);
7)对式(6)两端进行一阶摄动,建立刀尖点位置偏差方程:
8)建立转角偏差映射关系:
9)联立式(5)、(7)、(9),最终建立并联构型主轴头各驱动轴跟随误差与刀尖点位置偏差间的传递关系:
式中,T是一个无量纲方阵,表征并联构型主轴头的动态误差传递特性,具体表示为:
10)基于式(11),提出无量纲误差传递指标ETIj:
式中,ETIj为并联构型主轴头第j驱动轴的误差传递指标,j=1,2,3;tij为无量纲方阵T相应位置上的元素,i,j=1,2,3;S为无量纲方阵T的最大奇异值;∑||代表对元素取绝对值后进行求和;
11)利用无量纲误差传递指标ETIj对并联构型主轴头动态误差传递特性进行评价,ETIj的大小表征并联构型主轴头第j驱动轴的跟随误差对刀尖点位置偏差影响的强弱,ETIj的值越大,代表并联构型主轴头第j驱动轴的跟随误差对刀尖点位置偏差的影响越大。如果ETIj的值大于或接近1,说明此时并联构型主轴头第j驱动轴的跟随误差对刀尖点位置偏差的影响能力与最大奇异值相当。
为进一步说明所提出评价方法,图3所示为该并联构型主轴头沿β方向的ETIj变化情况,图像的横坐标代表摆角β变化范围,这里摆角β从-π/4变化至π/4,纵坐标为ETIj的数值。从图3中可以清楚看到,在摆角β从-π/4变化至π/4的过程中,ETI1与ETI2的数值始终大于ETI3的数值,因此,相较于第3驱动轴,第1驱动轴与第2驱动轴的跟随误差会对刀尖点位置偏差造成更大影响。此外,ETI1与ETI2的数值在一定范围内均接近或超过1,此时相应驱动轴的跟随误差对终端刀尖点位置偏差的影响能力与最大奇异值相当。从结果中可以看到,采用本发明提出的一种并联构型主轴头的动态误差传递方法能从数值上有效评价并联构型主轴头各驱动轴跟随误差对刀尖点偏差造成的影响,进而达到动态误差传递特性评价的目的。
Claims (2)
1.一种并联构型主轴头的动态误差传递特性评价方法,其特征在于该方法包括如下步骤:
1)建立并联构型主轴头各驱动轴跟随误差与终端动平台运动误差间的传递关系:
dX=Gadq
式中,dX=[dx dy dz dθx dθy dθz]T,dx、dy、dz分别为并联构型主轴头终端动平台沿X方向、Y方向和Z方向的平动误差,这里X、Y、Z为传统笛卡尔坐标系的三个方向;dθx、dθy、dθz分别为并联构型主轴头终端动平台绕X方向、Y方向和Z方向的转动误差;dq=[dq1 dq2 dq3]T,dq1、dq2、dq3分别为并联构型主轴头第1驱动轴、第2驱动轴和第3驱动轴的跟随误差;矩阵Ga具有6行3列元素,表征并联构型主轴头各驱动轴跟随误差与终端动平台运动误差间的传递关系;
2)基于并联构型主轴头刀尖点位置偏差,建立并联构型主轴头各驱动轴跟随误差与刀尖点位置偏差间的传递关系:
dP=Tdq
式中,dP=[dPx dPy dPz]T,dPx、dPy、dPz分别为刀尖点沿X方向、Y方向和Z方向的平动误差;T是一个无量纲方阵,具体表示为:
该无量纲方阵用于表征并联构型主轴头的动态误差传递特性;
3)基于步骤2)中的无量纲方阵T,提出无量纲误差传递指标ETIj:
式中,ETIj为并联构型主轴头第j驱动轴的误差传递指标,j=1,2,3;tij为无量纲方阵T相应位置上的元素,i=1,2,3,j=1,2,3;S为无量纲方阵T的最大奇异值;∑||代表对元素取绝对值后进行求和;
4)利用3)中提出的无量纲误差传递指标ETIj对并联构型主轴头动态误差传递特性进行评价,ETIj的大小表征并联构型主轴头第j驱动轴的跟随误差对刀尖点位置偏差影响的强弱,ETIj的值越大,代表并联构型主轴头第j驱动轴的跟随误差对刀尖点位置偏差的影响越大。
2.根据权利要求1所述的一种并联构型主轴头的动态误差传递特性评价方法,其特征在于,步骤2)中所述的建立并联构型主轴头各驱动轴跟随误差与刀尖点位置偏差间的传递关系的方法包括如下步骤:
(1)建立并联构型主轴头刀尖点位置方程:
式中,Px、Py、Pz分别为刀尖点沿X方向、Y方向和Z方向的位置;x、y、z分别为并联构型主轴头动平台沿X方向、Y方向和Z方向的位置;L为刀具长度;欧拉角α与β用于表示并联构型主轴头刀尖点的位姿;sα=sin(α),sβ=sin(β),cα=cos(α),cβ=cos(β);
(2)对步骤(1)中所建立的并联构型主轴头刀尖点位置方程进行一阶摄动,建立刀尖点位置偏差方程:
(3)建立转角偏差映射关系:
(4)联立以下三个方程:
dX=Gadq
最终建立并联构型主轴头各驱动轴跟随误差与刀尖点位置偏差间的传递关系:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710089785.8A CN106934108B (zh) | 2017-02-20 | 2017-02-20 | 一种并联构型主轴头的动态误差传递特性评价方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710089785.8A CN106934108B (zh) | 2017-02-20 | 2017-02-20 | 一种并联构型主轴头的动态误差传递特性评价方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106934108A CN106934108A (zh) | 2017-07-07 |
CN106934108B true CN106934108B (zh) | 2020-05-15 |
Family
ID=59423271
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710089785.8A Active CN106934108B (zh) | 2017-02-20 | 2017-02-20 | 一种并联构型主轴头的动态误差传递特性评价方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106934108B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108038337B (zh) * | 2017-12-29 | 2020-03-31 | 中国科学院沈阳自动化研究所 | 一种自动锁紧释放机构传动链精度分析方法 |
CN109100987B (zh) * | 2018-08-23 | 2020-03-31 | 西安交通大学 | 一种应用于数控机床的自主调整伺服参数的智能方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104772628A (zh) * | 2015-03-04 | 2015-07-15 | 清华大学 | 一种三自由度并联主轴头机构 |
CN105127840A (zh) * | 2015-09-08 | 2015-12-09 | 清华大学 | 一种五轴联动机床主轴头姿态角测量装置及测量方法 |
CN105302068A (zh) * | 2015-11-12 | 2016-02-03 | 西安交通大学 | 一种提高机床加工精度的设计方法 |
CN105892294A (zh) * | 2016-04-14 | 2016-08-24 | 清华大学 | 一种并联构型主轴头的伺服系统控制参数设计方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6775013B2 (en) * | 2001-03-07 | 2004-08-10 | Optodyne, Inc. | Method and apparatus for measuring displacement or motion error |
-
2017
- 2017-02-20 CN CN201710089785.8A patent/CN106934108B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104772628A (zh) * | 2015-03-04 | 2015-07-15 | 清华大学 | 一种三自由度并联主轴头机构 |
CN105127840A (zh) * | 2015-09-08 | 2015-12-09 | 清华大学 | 一种五轴联动机床主轴头姿态角测量装置及测量方法 |
CN105302068A (zh) * | 2015-11-12 | 2016-02-03 | 西安交通大学 | 一种提高机床加工精度的设计方法 |
CN105892294A (zh) * | 2016-04-14 | 2016-08-24 | 清华大学 | 一种并联构型主轴头的伺服系统控制参数设计方法 |
Non-Patent Citations (3)
Title |
---|
Dong Wang 等.A method for designing control parameters of a 3-DOF parallel tool head.《Mechatronics》.2016, * |
岳素平 等.7自由度虚拟轴串并联机床概念设计.《华北工学院学报》.2004, * |
王冬 等.3-PRS并联机器人惯量耦合特性研究.《力学学报》.2016, * |
Also Published As
Publication number | Publication date |
---|---|
CN106934108A (zh) | 2017-07-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109176494B (zh) | 绳驱多关节柔性机器人自标定方法及系统、存储介质 | |
CN107589934B (zh) | 一种关节型机械臂逆运动学解析解的求取方法 | |
Liao et al. | Region-based toolpath generation for robotic milling of freeform surfaces with stiffness optimization | |
CN107995885B (zh) | 一种坐标系标定方法 、系统及装置 | |
CN103495977B (zh) | 一种6r型工业机器人负载识别方法 | |
CN108015808B (zh) | 一种混联机器人的运动学标定方法 | |
CN113400088B (zh) | Ac双转台五轴机床的位置无关几何误差建模和辨识方法 | |
CN105014677A (zh) | 基于Camshift视觉跟踪和D-H建模算法的视觉机械臂控制装置及方法 | |
CN105014667A (zh) | 一种基于像素空间优化的相机与机器人相对位姿标定方法 | |
CN107471257B (zh) | 基于单拉线编码器的机器人几何标定方法 | |
CN106272416A (zh) | 基于力觉和视觉的机器人细长轴精密装配系统及方法 | |
CN106335057A (zh) | 一种基于实时力控的装配机器人全空间柔顺插孔控制方法 | |
CN110549338B (zh) | 一种圆-长方形复合孔类零件的机器人自动装配方法 | |
CN106934108B (zh) | 一种并联构型主轴头的动态误差传递特性评价方法 | |
CN105619389B (zh) | 一种无奇异位形的对称并联机构 | |
CN114131611A (zh) | 机器人重力位姿分解的关节误差离线补偿方法、系统及终端 | |
CN105598984A (zh) | 一种冗余度机械臂加速度层运动规划的初始化方法 | |
Wang et al. | Research on the error transfer characteristics of a 3-DOF parallel tool head | |
CN108717262A (zh) | 一种基于矩特征学习神经网络的异形曲面跟踪方法及系统 | |
CN108582047B (zh) | 一种六自由度混联抛光机器人位姿精度校准装置及方法 | |
CN206216693U (zh) | 一种针对解决形变问题的机器人磨削系统 | |
CN114147720B (zh) | 一种多自由度机械臂的逆运动学通用求解方法及装置 | |
Ghaffari et al. | Kinematic design of a novel 4-DOF parallel mechanism for turbine blade machining | |
CN110900608B (zh) | 基于最优测量构型选择的机器人运动学标定方法 | |
CN109866224A (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 |