CN106934108B - 一种并联构型主轴头的动态误差传递特性评价方法 - Google Patents

一种并联构型主轴头的动态误差传递特性评价方法 Download PDF

Info

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
parallel configuration
main shaft
error
deviation
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
CN201710089785.8A
Other languages
English (en)
Other versions
CN106934108A (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.)
Tsinghua University
Original Assignee
Tsinghua University
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 Tsinghua University filed Critical Tsinghua University
Priority to CN201710089785.8A priority Critical patent/CN106934108B/zh
Publication of CN106934108A publication Critical patent/CN106934108A/zh
Application granted granted Critical
Publication of CN106934108B publication Critical patent/CN106934108B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design

Abstract

本发明公开了一种并联构型主轴头的动态误差传递特性评价方法,属于机床技术领域。该方法首先推导并联构型主轴头各驱动轴跟随误差与终端动平台运动误差间的传递关系;进一步基于并联构型主轴头刀尖点位置偏差,建立各驱动轴跟随误差与刀尖点位置偏差间的传递关系,得到表征并联构型主轴头动态误差传递特性的无量纲方阵;基于该无量纲方阵,提出无量纲误差传递指标,该误差传递指标的大小表征各驱动轴跟随误差对刀尖点位置偏差影响的强弱。本发明克服了传统动态误差传递特性评价方法物理意义不明确的问题,通过提出无量纲误差传递指标,能够对并联构型主轴头的动态误差传递特性进行更为有效的评价。

Description

一种并联构型主轴头的动态误差传递特性评价方法
技术领域
本发明属于机床技术领域,特别涉及并联构型主轴头的动态误差传递特性评价问题。
背景技术
运动精度是面向工业加工的并联构型主轴头最重要的指标。相对于传统串联构型主轴头,并联构型主轴头因没有误差累积,在理论上具有更高的运动精度,因此有很多并联构型主轴头方案被相继提出,但大多数并联构型主轴头在实际中动态精度过低,不能满足加工需求。因此需要对并联构型主轴头的精度问题开展更深入的研究,这其中,并联构型主轴头的动态误差传递特性评价是一项核心内容。
在并联构型主轴头的运动过程中,各驱动轴的跟随误差是造成刀尖点位置偏差的主要因素,刀尖点位置偏差进而会对加工精度造成影响,因此需要对并联构型主轴头的动态误差传递特性进行有效评价。传统雅克比矩阵代表了并联构型主轴头各驱动轴跟随误差与终端动平台运动误差之间的映射关系,通常使用雅克比矩阵对并联构型主轴头的动态误差传递特性进行描述和评价,但是对于面向工业加工的并联构型主轴头,由于其终端动平台同时耦合了转动及平动自由度并存在伴生运动,导致雅克比矩阵存在以下两个显著问题:1)量纲问题,雅克比矩阵中元素的物理单位不一致;2)维数问题,雅克比矩阵不是一个方阵。由于量纲及维数问题,传统雅克比矩阵不再具有明确的物理意义,同时基于雅克比矩阵提出的评价指标不能够很好的评价并联构型主轴头的动态误差传递特性,如雅克比矩阵的最大奇异值、最小奇异值和条件数。
目前针对并联构型主轴头动态误差传递特性的评价方法较少,已有的评价方法及指标不能精确评价并联构型主轴头各驱动轴跟随误差对刀尖点位置偏差影响的强弱。因此为更好地分析、评价并联构型主轴头的动态误差传递特性,需要基于刀尖点位置偏差,建立新的误差传递矩阵,并进一步提出新的评价指标。
目前还没有普遍适用于并联构型主轴头的动态误差传递特性评价方法,针对并联构型主轴头的特点,提出一种广泛有效的动态误差传递特性评价方法对改善并联构型主轴头的动态精度具有重要意义。
发明内容
本发明的目的是为克服已有评价方法的不足之处,提出了一种并联构型主轴头的动态误差传递特性评价方法,用来解决并联构型主轴头的动态误差传递特性评价问题,从而更精确地评价并联构型主轴头各驱动轴跟随误差对刀尖点位置偏差的影响。
本发明的技术方案如下:
一种并联构型主轴头的动态误差传递特性评价方法,包括如下步骤:
1)建立并联构型主轴头各驱动轴跟随误差与终端动平台运动误差间的传递关系:
dX=Gadq
式中,dX=[dx dy dz dθxyz]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是一个无量纲方阵,具体表示为:
Figure GDA0002180961560000021
该无量纲方阵用于表征并联构型主轴头的动态误差传递特性;
3)基于步骤2)中的无量纲方阵T,提出无量纲误差传递指标ETIj
Figure GDA0002180961560000022
式中,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)建立并联构型主轴头刀尖点位置方程:
Figure GDA0002180961560000031
式中,Px、Py、Pz分别为刀尖点沿X方向、Y方向和Z方向的位置;x、y、z分别为并联构型主轴头动平台沿X方向、Y方向和Z方向的位置;L为刀具长度;欧拉角α与β用于表示并联构型主轴头刀尖点的位姿;sα=sin(α),sβ=sin(β),cα=cos(α),cβ=cos(β);
(2)对步骤(1)中所建立的并联构型主轴头刀尖点位置方程进行一阶摄动,建立刀尖点位置偏差方程:
Figure GDA0002180961560000032
式中,
Figure GDA0002180961560000033
dα与dβ分别为并联构型主轴头动平台沿欧拉角α与β方向的转动偏差;
(3)建立转角偏差映射关系:
Figure GDA0002180961560000034
式中,
Figure GDA0002180961560000035
dγ为并联构型主轴头动平台沿欧拉角γ方向的转动偏差;进一步有:
Figure GDA0002180961560000036
式中,
Figure GDA0002180961560000041
(4)联立以下三个方程:
dX=Gadq
Figure GDA0002180961560000042
Figure GDA0002180961560000043
最终建立并联构型主轴头各驱动轴跟随误差与刀尖点位置偏差间的传递关系:
Figure GDA0002180961560000044
本发明具有以下优点及突出性的技术效果:本发明通过建立各驱动轴跟随误差与刀尖点位置偏差间的传递关系,得到表征并联构型主轴头动态误差传递特性的无量纲方阵;基于该无量纲方阵,提出无量纲误差传递指标,该误差传递指标的大小表征各驱动轴跟随误差对刀尖点位置偏差影响的强弱;该方法克服了传统动态误差传递特性评价方法物理意义不明确的问题,通过提出无量纲误差传递指标,能够对并联构型主轴头的动态误差传递特性进行更为有效的评价;从数值上有效评价并联构型主头各驱动轴跟随误差对刀尖点偏差造成的影响,进而达到动态误差传递特性评价的目的。
附图说明
图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θxyz]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)的两端,建立并联构型主轴头的驱动方程:
Figure GDA0002180961560000051
式中,
Figure GDA0002180961560000052
为驱动雅克比矩阵,
Figure GDA0002180961560000053
3)用向量u1、u2、u3依次点乘方程(1)的两端,u1、u2、u3分别为第1被动杆件、第2被动杆件和第3被动杆件的转动轴轴向单位方向向量,建立并联构型主轴头的约束方程:
Figure GDA0002180961560000061
式中,
Figure GDA0002180961560000062
为约束雅克比矩阵;
4)联立式(2)、(3)得:
Figure GDA0002180961560000063
式中,
Figure GDA0002180961560000064
为广义雅克比矩阵;
5)由式(4)进一步得到:
dX=Gadq (5)
式中,矩阵Ga具有6行3列元素,表征并联构型主轴头各驱动轴跟随误差与终端动平台运动误差间的传递关系;
6)建立并联构型主轴头刀尖点位置方程:
Figure GDA0002180961560000065
式中,Px、Py、Pz分别为刀尖点沿X方向、Y方向和Z方向的位置;x、y、z分别为并联构型主轴头动平台沿X方向、Y方向和Z方向的位置;L为刀具长度;欧拉角α与β用于表示并联构型主轴头刀尖点的位姿,sα=sin(α),sβ=sin(β),cα=cos(α),cβ=cos(β);
7)对式(6)两端进行一阶摄动,建立刀尖点位置偏差方程:
Figure GDA0002180961560000066
式中,
Figure GDA0002180961560000071
dα与dβ分别为并联构型主轴头动平台沿欧拉角α与β方向的转动偏差;
8)建立转角偏差映射关系:
Figure GDA0002180961560000072
式中,
Figure GDA0002180961560000073
dγ为并联构型主轴头动平台沿欧拉角γ方向的转动偏差;进一步有:
Figure GDA0002180961560000074
式中,
Figure GDA0002180961560000075
9)联立式(5)、(7)、(9),最终建立并联构型主轴头各驱动轴跟随误差与刀尖点位置偏差间的传递关系:
Figure GDA0002180961560000076
式中,T是一个无量纲方阵,表征并联构型主轴头的动态误差传递特性,具体表示为:
Figure GDA0002180961560000077
10)基于式(11),提出无量纲误差传递指标ETIj
Figure GDA0002180961560000078
式中,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θxyz]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是一个无量纲方阵,具体表示为:
Figure FDA0002180961550000011
该无量纲方阵用于表征并联构型主轴头的动态误差传递特性;
3)基于步骤2)中的无量纲方阵T,提出无量纲误差传递指标ETIj
Figure FDA0002180961550000012
式中,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)建立并联构型主轴头刀尖点位置方程:
Figure FDA0002180961550000021
式中,Px、Py、Pz分别为刀尖点沿X方向、Y方向和Z方向的位置;x、y、z分别为并联构型主轴头动平台沿X方向、Y方向和Z方向的位置;L为刀具长度;欧拉角α与β用于表示并联构型主轴头刀尖点的位姿;sα=sin(α),sβ=sin(β),cα=cos(α),cβ=cos(β);
(2)对步骤(1)中所建立的并联构型主轴头刀尖点位置方程进行一阶摄动,建立刀尖点位置偏差方程:
Figure FDA0002180961550000022
式中,
Figure FDA0002180961550000023
dα与dβ分别为并联构型主轴头动平台沿欧拉角α与β方向的转动偏差;
(3)建立转角偏差映射关系:
Figure FDA0002180961550000024
式中,
Figure FDA0002180961550000025
dγ为并联构型主轴头动平台沿欧拉角γ方向的转动偏差;进一步有:
Figure FDA0002180961550000031
式中,
Figure FDA0002180961550000032
(4)联立以下三个方程:
dX=Gadq
Figure FDA0002180961550000033
Figure FDA0002180961550000034
最终建立并联构型主轴头各驱动轴跟随误差与刀尖点位置偏差间的传递关系:
Figure FDA0002180961550000035
CN201710089785.8A 2017-02-20 2017-02-20 一种并联构型主轴头的动态误差传递特性评价方法 Active CN106934108B (zh)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

Patent Citations (4)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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
CN108015808B (zh) 一种混联机器人的运动学标定方法
CN102566577B (zh) 一种工业机器人简易标定方法
CN110161850A (zh) 一种工业机器人变参数刚度辨识与建模方法
CN109434829B (zh) 一种立体石雕机器人加工系统的变形预测及补偿方法
CN109877840B (zh) 一种基于相机光轴约束的双机械臂标定方法
CN107598919B (zh) 一种基于五点标定法的两轴变位机标定方法
CN109176494A (zh) 绳驱多关节柔性机器人自标定方法及系统、存储介质
CN108789404A (zh) 一种基于视觉的串联机器人运动学参数标定方法
CN107471257B (zh) 基于单拉线编码器的机器人几何标定方法
CN105014677A (zh) 基于Camshift视觉跟踪和D-H建模算法的视觉机械臂控制装置及方法
CN104890013A (zh) 一种基于拉线编码器的工业机器人标定算法
CN113580148B (zh) 基于等效运动链的并联机器人运动学标定方法
CN107791248B (zh) 基于不满足Pieper准则的六自由度串联机器人的控制方法
CN106934108B (zh) 一种并联构型主轴头的动态误差传递特性评价方法
CN113400088B (zh) Ac双转台五轴机床的位置无关几何误差建模和辨识方法
CN113715016B (zh) 一种基于3d视觉的机器人抓取方法、系统、装置及介质
CN105856231A (zh) 一种特定构型六轴工业机器人的运动控制方法
CN111679629A (zh) 一种多主轴头加工的空行程无干涉轨迹规划方法
CN108339995B (zh) 一种可调角度的机器人末端执行器安装机构及其调整方法
Wang et al. Research on the error transfer characteristics of a 3-DOF parallel tool head
CN109866224A (zh) 一种机器人雅可比矩阵计算方法、装置及存储介质
CN110428436B (zh) 一种钢卷贴标喷码六轴机器人选型与安装位置的确定方法
CN106735581B (zh) Ac五轴水切割误差补偿的方法及系统
CN112276940A (zh) 一种六自由度非球型手腕机器人逆运动学求解方法
CN116276990A (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