CN109933918B - 旋转轴垂直度误差的误差矢量建模方法 - Google Patents

旋转轴垂直度误差的误差矢量建模方法 Download PDF

Info

Publication number
CN109933918B
CN109933918B CN201910203988.4A CN201910203988A CN109933918B CN 109933918 B CN109933918 B CN 109933918B CN 201910203988 A CN201910203988 A CN 201910203988A CN 109933918 B CN109933918 B CN 109933918B
Authority
CN
China
Prior art keywords
axis
error
motion
perpendicularity
matrix
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
CN201910203988.4A
Other languages
English (en)
Other versions
CN109933918A (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.)
Southwest Jiaotong University
Original Assignee
Southwest Jiaotong 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 Southwest Jiaotong University filed Critical Southwest Jiaotong University
Priority to CN201910203988.4A priority Critical patent/CN109933918B/zh
Publication of CN109933918A publication Critical patent/CN109933918A/zh
Application granted granted Critical
Publication of CN109933918B publication Critical patent/CN109933918B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Numerical Control (AREA)

Abstract

本发明公开了一种旋转轴垂直度误差的误差矢量建模方法,包括S1、根据垂直度误差定义结合旋转轴运动性质分析垂直度误差对旋转轴运动的影响;S2、得到A轴、B轴和C轴垂直度误差;S3、根据不同旋转轴运动性质构建A轴、B轴和C轴的理想指数运动矩阵;S4、根据不同旋转轴垂直度误差分布,采用指数积理论建立A轴、B轴和C轴垂直度误差影响下的实际指数运动矩阵;S5、根据运动轴误差矩阵、理想运动矩阵和实际运动矩阵之间关系构建不同旋转轴的垂直度误差变换矩阵;S6、将A轴、B轴和C轴的垂直度误差变换矩阵与运动轴基本几何误差项的误差矩阵进行比较,得到不同旋转轴的垂直度误差的影响;S7、建立不同旋转轴垂直度误差的误差矢量。

Description

旋转轴垂直度误差的误差矢量建模方法
技术领域
本发明属于数控机床的技术领域,具体涉及一种旋转轴垂直度误差的误差矢量建模方法。
背景技术
五轴数控机床已经广泛应用于航天、航空、航海、汽车、国防等各种领域,五轴数控机床高精度制造是衡量一个国家先进制造水平的重要标志。高精度是数控机床发展的趋势之一,旋转轴作为五轴数控机床重要功能部件,旋转轴垂直度误差成为机床几何误差项中重要组成部分之一,是影响机床精度的关键因素之一。在机床误差建模过程中,旋转轴垂直度误差的几何性质使得其难以描述,很大程度上影响了机床综合几何误差模型精度。另外,垂直度误差模型关系到误差测量辨识的精度,同时也影响到后续误差补偿效果。
常用的数控机床几何误差建模方法是基于多体系统理论的齐次变换矩阵(HTM)建模方法。第一种常用方法是将垂直度误差转化为相应旋转轴的线性误差来建立齐次变换矩阵,由于旋转轴运动量为旋转角度,该方法不符合旋转轴运动性质。第二种方法是将垂直度误差模型作为旋转轴转角误差的一部分进行建模,第三种方法是将垂直度误差作为运动轴与位置无关的误差进行处理,建立位置误差矩阵,这两种方法都将垂直度误差作为角度误差处理,不符合垂直度误差的几何定义,也没有反映垂直度误差对机床精度的影响。
发明内容
本发明的目的在于针对现有技术中的上述不足,提供一种旋转轴垂直度误差的误差矢量建模方法,以解决或改善上述问题。
为达到上述目的,本发明采取的技术方案是:
一种旋转轴垂直度误差的误差矢量建模方法,其包括:
S1、根据垂直度误差定义结合旋转轴运动性质分析垂直度误差对旋转轴运动的影响;
S2、分析旋转轴垂直度误差定义,得到A轴、B轴和C轴垂直度误差;
S3、根据不同旋转轴运动性质构建A轴、B轴和C轴的理想指数运动矩阵;
S4、根据不同旋转轴垂直度误差分布,采用指数积理论建立A轴、B轴和C轴垂直度误差影响下的实际指数运动矩阵;
S5、根据运动轴误差矩阵、理想运动矩阵和实际运动矩阵之间关系构建不同旋转轴的垂直度误差变换矩阵;
S6、将A轴、B轴和C轴的垂直度误差变换矩阵与运动轴基本几何误差项的误差矩阵进行比较,得到不同旋转轴的垂直度误差的影响;
S7、建立不同旋转轴垂直度误差的误差矢量。
优选地,步骤S1中旋转轴运动性质是驱动旋转轴围绕轴本身坐标系中某条坐标轴进行转动,本身坐标系原点位置不发生变化;
垂直度误差对旋转轴运动的影响是偏转旋转轴的转动轴线方向,使得旋转轴围绕另一条由垂直度误差决定的直线旋转。
优选地,步骤S2中分析机床不同旋转轴垂直度误差定义,得到A轴、B轴和C轴的垂直度误差的分布为:
A轴垂直度误差为与Y轴之间的垂直度误差Say和与Z轴之间的垂直度误差Saz;B轴垂直度误差为与X轴之间的垂直度误差Sbx和与Z轴之间的垂直度误差Sbz;C轴垂直度误差为与X轴之间的垂直度误差Scx和与Y轴之间的垂直度误差Scy
优选地,步骤S3中构建A轴、B轴和C轴的理想指数运动矩阵的方法为:
S3.1、根据A轴理想运动为围绕理想坐标系x轴旋转且旋转轴线是理想坐标系x轴的运动性质,得到A轴运动方向向量ωai=[1,0,0]T,建立A轴运动旋量
Figure GDA0002024902050000031
为:
Figure GDA0002024902050000032
根据指数积理论得到A轴理想指数运动矩阵Tai为:
Figure GDA0002024902050000033
其中,e表示指数积理论中旋量的指数矩阵函数,α表示A轴运动量;
S3.2、结合B轴理想运动为围绕理想坐标系y轴旋转且旋转轴线是理想坐标系y轴的运动性质,得到B轴运动方向向量ωbi=[0,1,0]T,建立B轴运动旋量
Figure GDA0002024902050000034
为:
Figure GDA0002024902050000035
根据指数积理论得到B轴理想指数运动矩阵Tbi为:
Figure GDA0002024902050000036
其中,β表示B轴运动量;
S3.3、结合C轴理想运动为围绕理想坐标系z轴旋转且旋转轴线是理想坐标系z轴的运动性质,得到C轴运动方向向量ωci=[0,0,1]T,建立C轴运动旋量
Figure GDA0002024902050000037
为:
Figure GDA0002024902050000038
根据指数积理论得到C轴理想指数运动矩阵Tci为:
Figure GDA0002024902050000039
其中,γ表示C轴运动量。
优选地,步骤S4中建立A轴、B轴和C轴垂直度误差影响下的实际指数运动矩阵的方法为:
S4.1、得到A轴垂直度误差影响下的运动性质为旋转轴线方向发生偏转但是轴线原点位置不变,得到A轴垂直度误差影响下运动方向向量ωar为:
ωar=[cos(Saz)cos(Say),-cos(Saz)sin(Say),-sin(Saz)]T
建立A轴垂直度误差影响下的运动旋量
Figure GDA0002024902050000041
为:
Figure GDA0002024902050000042
根据指数积理论得到A轴垂直度误差影响下的实际指数运动矩阵Tar为:
Figure GDA0002024902050000043
S4.2、得到B轴垂直度误差影响下的运动性质为为旋转轴线方向发生偏转但是轴线原点位置不变,得到B轴运动方向向量ωbr为:
ωbr=[-sin(Sbx),cos(Sbx)cos(Sbz),-cos(Sbx)sin(Sbz)]T
建立B轴垂直度误差影响下的运动旋量
Figure GDA0002024902050000044
为:
Figure GDA0002024902050000045
根据指数积理论得到B轴垂直度误差影响下的实际指数运动矩阵Tbr为:
Figure GDA0002024902050000046
S4.3、得到C轴垂直度误差影响下的运动性质为为旋转轴线方向发生偏转但是轴线原点位置不变,得到C轴运动方向向量ωcr为:
ωcr=[-sin(Scx),-cos(Scx)sin(Scy),cos(Scx)cos(Scy)]T
建立C轴垂直度误差影响下的运动旋量
Figure GDA0002024902050000047
为:
Figure GDA0002024902050000048
根据指数积理论得到C轴垂直度误差影响下的实际指数运动矩阵Tcr为:
Figure GDA0002024902050000049
优选地,步骤S5中构建不同旋转轴的垂直度误差变换矩阵的方法为:
S5.1、根据运动轴误差矩阵、理想运动矩阵和实际运动矩阵之间关系建立运动轴误差矩阵Tme计算公式为:
Figure GDA00020249020500000410
其中,Tmi表示运动轴m的理想运动矩阵,Tmr表示运动轴m的实际运动矩阵;
S5.2、代入A轴理想指数运动矩阵和A轴垂直度误差影响下的实际指数运动矩阵,忽略高阶误差项建立A轴垂直度误差变换矩阵Tae为:
Figure GDA0002024902050000051
S5.3、代入B轴理想指数运动矩阵和B轴垂直度误差影响下的实际指数运动矩阵,忽略高阶误差项建立B轴垂直度误差变换矩阵Tbe为:
Figure GDA0002024902050000052
S5.4、代入C轴理想指数运动矩阵和C轴垂直度误差影响下的实际指数运动矩阵,忽略高阶误差项建立C轴垂直度误差变换矩阵Tce为:
Figure GDA0002024902050000053
优选地,步骤S6中得到不同旋转轴的垂直度误差的影响的方法为:
S6.1、建立运动轴基本几何误差项的误差矩阵为:
Figure GDA0002024902050000054
其中,Tmse表示运动轴m的基本几何误差项的误差矩阵,δkm表示运动轴m在k方向上的线性几何误差,εkm表示运动轴m在k方向上的角度几何误差,k=x,y,z;
S6.2、比较A轴垂直度误差变换矩阵与运动轴基本几何误差项的误差矩阵,得到A轴垂直度误差的影响为:引起y方向上的角度几何误差为–Saz+Saz cosα+Saysinα,引起z方向上的角度几何误差为Say–Say cosα–Saz sinα;
S6.3、比较B轴垂直度误差变换矩阵与运动轴基本几何误差项的误差矩阵,得到B轴垂直度误差的影响为:引起x方向上的角度几何误差为Sbz–Sbz cosβ–Sbxsinβ,引起z方向上的角度几何误差为–Sbx+Sbx cosβ–Sbz sinβ;
S6.4、比较C轴垂直度误差变换矩阵与运动轴基本几何误差项的误差矩阵,得到B轴垂直度误差的影响为:引起x方向上的角度几何误差为–Scy+Scy cosγ–Scxsinγ,引起y方向上的角度几何误差为Scx–Scxcosγ–Scysinγ。
优选地,步骤S7中建立不同旋转轴垂直度误差的误差矢量的方法为:
S7.1、建立运动轴基本几何误差项的误差矢量ΔEm表达式:
ΔEm=[δxmymzmxmymzm]T
S7.2、代入A轴垂直度误差的影响,建立A轴垂直度误差的误差矢量ΔEAe为:
ΔEAe=[0,0,0,0,-Saz+Sazcosα+Saysinα,Say-Saycosα-Sazsinα]T
S7.3、代入B轴垂直度误差的影响,建立B轴垂直度误差的误差矢量ΔEBe为:
ΔEBe=[0,0,0,Sbz-Sbzcosβ-Sbxsinβ,0,-Sbx+Sbxcosβ-Sbzsinβ]T
S7.4、代入C轴垂直度误差的影响,建立C轴垂直度误差的误差矢量ΔECe为:
ΔECe=[0,0,0,-Scy+Scycosγ-Scxsinγ,Scx-Scxcosγ-Scysinγ,0]T
本发明提供的旋转轴垂直度误差的误差矢量建模方法,具有以下有益效果:
本发明得到的垂直度误差模型考虑了垂直度误差对机床精度的影响,反映了垂直度误差的几何定义,可进一步得到更加精确的旋转轴包含垂直度误差的误差矩阵,可进一步提高机床综合几何误差模型的精度。
附图说明
图1为旋转轴垂直度误差的误差矢量建模方法的流程图。
图2为旋转轴垂直度误差的误差矢量建模方法C轴垂直度误差对C轴运动的影响示意图。
图3为旋转轴垂直度误差的误差矢量建模方法C轴垂直度误差分布图。
具体实施方式
下面对本发明的具体实施方式进行描述,以便于本技术领域的技术人员理解本发明,但应该清楚,本发明不限于具体实施方式的范围,对本技术领域的普通技术人员来讲,只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内,这些变化是显而易见的,一切利用本发明构思的发明创造均在保护之列。
根据本申请的一个实施例,参考图1,本方案的旋转轴垂直度误差的误差矢量建模方法,包括:
S1、根据垂直度误差定义结合旋转轴运动性质分析垂直度误差对旋转轴运动的影响;
S2、分析旋转轴垂直度误差定义,得到A轴、B轴和C轴垂直度误差;
S3、根据不同旋转轴运动性质构建A轴、B轴和C轴的理想指数运动矩阵;
S4、根据不同旋转轴垂直度误差分布,采用指数积理论建立A轴、B轴和C轴垂直度误差影响下的实际指数运动矩阵;
S5、根据运动轴误差矩阵、理想运动矩阵和实际运动矩阵之间关系构建不同旋转轴的垂直度误差变换矩阵;
S6、将A轴、B轴和C轴的垂直度误差变换矩阵与运动轴基本几何误差项的误差矩阵进行比较,得到不同旋转轴的垂直度误差的影响;
S7、建立不同旋转轴垂直度误差的误差矢量。
以下对上述各个步骤进行详细描述
S1、根据垂直度误差定义结合旋转轴运动性质分析垂直度误差对旋转轴运动的影响;
其中,旋转轴运动性质是驱动旋转轴围绕轴本身坐标系中某条坐标轴进行转动,本身坐标系原点位置不会发生变化。
垂直度误差对旋转轴运动的影响是偏转旋转轴的转动轴线方向,使得旋转轴围绕另一条由垂直度误差决定的直线旋转,图2为C轴垂直度误差对C轴运动的影响示意图。
S2、分析旋转轴垂直度误差定义,得到A轴、B轴和C轴垂直度误差的分布为:
A轴垂直度误差为与Y轴之间的垂直度误差Say和与Z轴之间的垂直度误差Saz;B轴垂直度误差为与X轴之间的垂直度误差Sbx和与Z轴之间的垂直度误差Sbz;C轴垂直度误差为与X轴之间的垂直度误差Scx和与Y轴之间的垂直度误差Scy;附图3为C轴垂直度误差分布图。
S3、根据不同旋转轴运动性质构建A轴、B轴和C轴的理想指数运动矩阵,其具体步骤包括:
S3.1、根据A轴理想运动为围绕理想坐标系x轴旋转且旋转轴线是理想坐标系x轴的运动性质,得到A轴运动方向向量ωai=[1,0,0]T,建立A轴运动旋量
Figure GDA0002024902050000081
为:
Figure GDA0002024902050000082
根据指数积理论得到A轴理想指数运动矩阵Tai为:
Figure GDA0002024902050000083
其中,e表示指数积理论中旋量的指数矩阵函数,α表示A轴运动量;
S3.2、结合B轴理想运动为围绕理想坐标系y轴旋转且旋转轴线是理想坐标系y轴的运动性质,得到B轴运动方向向量ωbi=[0,1,0]T,建立B轴运动旋量
Figure GDA0002024902050000084
为:
Figure GDA0002024902050000091
根据指数积理论得到B轴理想指数运动矩阵Tbi为:
Figure GDA0002024902050000092
其中,β表示B轴运动量;
S3.3、结合C轴理想运动为围绕理想坐标系z轴旋转且旋转轴线是理想坐标系z轴的运动性质,得到C轴运动方向向量ωci=[0,0,1]T,建立C轴运动旋量
Figure GDA0002024902050000093
为:
Figure GDA0002024902050000094
根据指数积理论得到C轴理想指数运动矩阵Tci为:
Figure GDA0002024902050000095
其中,γ表示C轴运动量。
S4、根据不同旋转轴垂直度误差分布,采用指数积理论建立A轴、B轴和C轴垂直度误差影响下的实际指数运动矩阵,其具体步骤包括:
S4.1、得到A轴垂直度误差影响下的运动性质为旋转轴线方向发生偏转但是轴线原点位置不变,得到A轴垂直度误差影响下运动方向向量ωar为:
ωar=[cos(Saz)cos(Say),-cos(Saz)sin(Say),-sin(Saz)]T
建立A轴垂直度误差影响下的运动旋量
Figure GDA0002024902050000096
为:
Figure GDA0002024902050000097
根据指数积理论得到A轴垂直度误差影响下的实际指数运动矩阵Tar为:
Figure GDA0002024902050000098
S4.2、得到B轴垂直度误差影响下的运动性质为为旋转轴线方向发生偏转但是轴线原点位置不变,得到B轴运动方向向量ωbr为:
ωbr=[-sin(Sbx),cos(Sbx)cos(Sbz),-cos(Sbx)sin(Sbz)]T
建立B轴垂直度误差影响下的运动旋量
Figure GDA0002024902050000099
为:
Figure GDA0002024902050000101
根据指数积理论得到B轴垂直度误差影响下的实际指数运动矩阵Tbr为:
Figure GDA0002024902050000102
S4.3、得到C轴垂直度误差影响下的运动性质为为旋转轴线方向发生偏转但是轴线原点位置不变,得到C轴运动方向向量ωcr为:
ωcr=[-sin(Scx),-cos(Scx)sin(Scy),cos(Scx)cos(Scy)]T
建立C轴垂直度误差影响下的运动旋量
Figure GDA0002024902050000103
为:
Figure GDA0002024902050000104
根据指数积理论得到C轴垂直度误差影响下的实际指数运动矩阵Tcr为:
Figure GDA0002024902050000105
S5、根据运动轴误差矩阵、理想运动矩阵和实际运动矩阵之间关系构建不同旋转轴的垂直度误差变换矩阵,其具体步骤包括:
S6、将A轴、B轴和C轴的垂直度误差变换矩阵与运动轴基本几何误差项的误差矩阵进行比较,得到不同旋转轴的垂直度误差的影响,其具体步骤包括:
S6.1、建立运动轴基本几何误差项的误差矩阵为:
Figure GDA0002024902050000106
其中,Tmse表示运动轴m的基本几何误差项的误差矩阵,δkm表示运动轴m在k方向上的线性几何误差,εkm表示运动轴m在k方向上的角度几何误差,k=x,y,z;
S6.2、比较A轴垂直度误差变换矩阵与运动轴基本几何误差项的误差矩阵,得到A轴垂直度误差的影响为:引起y方向上的角度几何误差为–Saz+Saz cosα+Say sinα,引起z方向上的角度几何误差为Say–Say cosα–Saz sinα;
S6.3、比较B轴垂直度误差变换矩阵与运动轴基本几何误差项的误差矩阵,得到B轴垂直度误差的影响为:引起x方向上的角度几何误差为Sbz–Sbz cosβ–Sbx sinβ,引起z方向上的角度几何误差为–Sbx+Sbx cosβ–Sbz sinβ;
S6.4、比较C轴垂直度误差变换矩阵与运动轴基本几何误差项的误差矩阵,得到B轴垂直度误差的影响为:引起x方向上的角度几何误差为–Scy+Scy cosγ–Scx sinγ,引起y方向上的角度几何误差为Scx–Scxcosγ–Scysinγ。
S7、建立不同旋转轴垂直度误差的误差矢量,其具体步骤包括:
S7.1、建立运动轴基本几何误差项的误差矢量ΔEm表达式:
ΔEm=[δxmymzmxmymzm]T
S7.2、代入A轴垂直度误差的影响,建立A轴垂直度误差的误差矢量ΔEAe为:
ΔEAe=[0,0,0,0,-Saz+Sazcosα+Saysinα,Say-Saycosα-Sazsinα]T
S7.3、代入B轴垂直度误差的影响,建立B轴垂直度误差的误差矢量ΔEBe为:
ΔEBe=[0,0,0,Sbz-Sbzcosβ-Sbxsinβ,0,-Sbx+Sbxcosβ-Sbzsinβ]T
S7.4、代入C轴垂直度误差的影响,建立C轴垂直度误差的误差矢量ΔECe为:
ΔECe=[0,0,0,-Scy+Scycosγ-Scxsinγ,Scx-Scxcosγ-Scysinγ,0]T
本发明在得到不同旋转轴垂直度误差的影响和垂直度误差的误差矢量后,可以进一步结合旋转轴基本误差项的误差矩阵建立旋转轴包含垂直度误差的误差矩阵,这样可以进一步提高误差模型的精度,对于建立高精度机床误差模型,提出高精度误差补偿方法具有重要意义。
虽然结合附图对发明的具体实施方式进行了详细地描述,但不应理解为对本专利的保护范围的限定。在权利要求书所描述的范围内,本领域技术人员不经创造性劳动即可做出的各种修改和变形仍属本专利的保护范围。

Claims (8)

1.一种旋转轴垂直度误差的误差矢量建模方法,其特征在于,包括:
S1、根据垂直度误差定义结合旋转轴运动性质分析垂直度误差对旋转轴运动的影响;
S2、分析旋转轴垂直度误差定义,得到A轴、B轴和C轴垂直度误差;
S3、根据不同旋转轴运动性质构建A轴、B轴和C轴的理想指数运动矩阵;
S4、根据不同旋转轴垂直度误差分布,采用指数积理论建立A轴、B轴和C轴垂直度误差影响下的实际指数运动矩阵;
S5、根据运动轴误差矩阵、理想运动矩阵和实际运动矩阵之间关系构建不同旋转轴的垂直度误差变换矩阵;
S6、将A轴、B轴和C轴的垂直度误差变换矩阵与运动轴基本几何误差项的误差矩阵进行比较,得到不同旋转轴的垂直度误差的影响;
S7、建立不同旋转轴垂直度误差的误差矢量。
2.根据权利要求1所述的旋转轴垂直度误差的误差矢量建模方法,其特征在于,所述步骤S1中旋转轴运动性质是驱动旋转轴围绕轴本身坐标系中某条坐标轴进行转动,本身坐标系原点位置不发生变化;
垂直度误差对旋转轴运动的影响是偏转旋转轴的转动轴线方向,使得旋转轴围绕另一条由垂直度误差决定的直线旋转。
3.根据权利要求1所述的旋转轴垂直度误差的误差矢量建模方法,其特征在于,所述步骤S2中分析旋转轴垂直度误差定义,得到A轴、B轴和C轴垂直度误差的分布为:
A轴垂直度误差为与Y轴之间的垂直度误差Say和与Z轴之间的垂直度误差Saz;B轴垂直度误差为与X轴之间的垂直度误差Sbx和与Z轴之间的垂直度误差Sbz;C轴垂直度误差为与X轴之间的垂直度误差Scx和与Y轴之间的垂直度误差Scy
4.根据权利要求1所述的旋转轴垂直度误差的误差矢量建模方法,其特征在于,所述步骤S3中构建A轴、B轴和C轴的理想指数运动矩阵的方法为:
S3.1、根据A轴理想运动为围绕理想坐标系x轴旋转且旋转轴线是理想坐标系x轴的运动性质,得到A轴运动方向向量ωai=[1,0,0]T,建立A轴运动旋量
Figure FDA0002919786850000021
为:
Figure FDA0002919786850000022
根据指数积理论得到A轴理想指数运动矩阵Tai为:
Figure FDA0002919786850000023
其中,e表示指数积理论中旋量的指数矩阵函数,α表示A轴运动量;
S3.2、结合B轴理想运动为围绕理想坐标系y轴旋转且旋转轴线是理想坐标系y轴的运动性质,得到B轴运动方向向量ωbi=[0,1,0]T,建立B轴运动旋量
Figure FDA0002919786850000024
为:
Figure FDA0002919786850000025
根据指数积理论得到B轴理想指数运动矩阵Tbi为:
Figure FDA0002919786850000026
其中,β表示B轴运动量;
S3.3、结合C轴理想运动为围绕理想坐标系z轴旋转且旋转轴线是理想坐标系z轴的运动性质,得到C轴运动方向向量ωci=[0,0,1]T,建立C轴运动旋量
Figure FDA0002919786850000027
为:
Figure FDA0002919786850000028
根据指数积理论得到C轴理想指数运动矩阵Tci为:
Figure FDA0002919786850000029
其中,γ表示C轴运动量。
5.根据权利要求1所述的旋转轴垂直度误差的误差矢量建模方法,其特征在于,所述步骤S4中建立A轴、B轴和C轴垂直度误差影响下的实际指数运动矩阵的方法为:
S4.1、得到A轴垂直度误差影响下的运动性质为旋转轴线方向发生偏转但是轴线原点位置不变,得到A轴垂直度误差影响下运动方向向量ωar为:
ωar=[cos(Saz)cos(Say),-cos(Saz)sin(Say),-sin(Saz)]T
建立A轴垂直度误差影响下的运动旋量
Figure FDA0002919786850000031
为:
Figure FDA0002919786850000032
根据指数积理论得到A轴垂直度误差影响下的实际指数运动矩阵Tar为:
Figure FDA0002919786850000033
S4.2、得到B轴垂直度误差影响下的运动性质为旋转轴线方向发生偏转但是轴线原点位置不变,得到B轴运动方向向量ωbr为:
ωbr=[-sin(Sbx),cos(Sbx)cos(Sbz),-cos(Sbx)sin(Sbz)]T
建立B轴垂直度误差影响下的运动旋量
Figure FDA0002919786850000034
为:
Figure FDA0002919786850000035
根据指数积理论得到B轴垂直度误差影响下的实际指数运动矩阵Tbr为:
Figure FDA0002919786850000036
S4.3、得到C轴垂直度误差影响下的运动性质为旋转轴线方向发生偏转但是轴线原点位置不变,得到C轴运动方向向量ωcr为:
ωcr=[-sin(Scx),-cos(Scx)sin(Scy),cos(Scx)cos(Scy)]T
建立C轴垂直度误差影响下的运动旋量
Figure FDA0002919786850000037
为:
Figure FDA0002919786850000038
根据指数积理论得到C轴垂直度误差影响下的实际指数运动矩阵Tcr为:
Figure FDA0002919786850000039
6.根据权利要求1所述的旋转轴垂直度误差的误差矢量建模方法,其特征在于,所述步骤S5中构建不同旋转轴的垂直度误差变换矩阵的方法为:
S5.1、根据运动轴误差矩阵、理想运动矩阵和实际运动矩阵之间关系建立运动轴误差矩阵Tme计算公式为:
Figure FDA0002919786850000041
其中,Tmi表示运动轴m的理想运动矩阵,Tmr表示运动轴m的实际运动矩阵;
S5.2、代入A轴理想指数运动矩阵和A轴垂直度误差影响下的实际指数运动矩阵,忽略高阶误差项建立A轴垂直度误差变换矩阵Tae为:
Figure FDA0002919786850000042
S5.3、代入B轴理想指数运动矩阵和B轴垂直度误差影响下的实际指数运动矩阵,忽略高阶误差项建立B轴垂直度误差变换矩阵Tbe为:
Figure FDA0002919786850000043
S5.4、代入C轴理想指数运动矩阵和C轴垂直度误差影响下的实际指数运动矩阵,忽略高阶误差项建立C轴垂直度误差变换矩阵Tce为:
Figure FDA0002919786850000044
7.根据权利要求1所述的旋转轴垂直度误差的误差矢量建模方法,其特征在于,所述步骤S6中得到不同旋转轴的垂直度误差的影响的方法为:
S6.1、建立运动轴基本几何误差项的误差矩阵为:
Figure FDA0002919786850000051
其中,Tmse表示运动轴m的基本几何误差项的误差矩阵,δkm表示运动轴m在k方向上的线性几何误差,εkm表示运动轴m在k方向上的角度几何误差,k=x,y,z;
S6.2、比较A轴垂直度误差变换矩阵与运动轴基本几何误差项的误差矩阵,得到A轴垂直度误差的影响为:引起y方向上的角度几何误差为–Saz+Sazcosα+Saysinα,引起z方向上的角度几何误差为Say–Saycosα–Sazsinα;
S6.3、比较B轴垂直度误差变换矩阵与运动轴基本几何误差项的误差矩阵,得到B轴垂直度误差的影响为:引起x方向上的角度几何误差为Sbz–Sbzcosβ–Sbxsinβ,引起z方向上的角度几何误差为–Sbx+Sbxcosβ–Sbzsinβ;
S6.4、比较C轴垂直度误差变换矩阵与运动轴基本几何误差项的误差矩阵,得到B轴垂直度误差的影响为:引起x方向上的角度几何误差为–Scy+Scycosγ–Scxsinγ,引起y方向上的角度几何误差为Scx–Scxcosγ–Scysinγ。
8.根据权利要求1所述的旋转轴垂直度误差的误差矢量建模方法,其特征在于,所述步骤S7中建立不同旋转轴垂直度误差的误差矢量的方法为:
S7.1、建立运动轴基本几何误差项的误差矢量△Em表达式:
ΔEm=[δxmymzmxmymzm]T
S7.2、代入A轴垂直度误差的影响,建立A轴垂直度误差的误差矢量△EAe为:
ΔEAe=[0,0,0,0,-Saz+Sazcosα+Saysinα,Say-Saycosα-Sazsinα]T
S7.3、代入B轴垂直度误差的影响,建立B轴垂直度误差的误差矢量△EBe为:
ΔEBe=[0,0,0,Sbz-Sbzcosβ-Sbxsinβ,0,-Sbx+Sbxcosβ-Sbzsinβ]T
S7.4、代入C轴垂直度误差的影响,建立C轴垂直度误差的误差矢量△ECe为:
ΔECe=[0,0,0,-Scy+Scycosγ-Scxsinγ,Scx-Scxcosγ-Scysinγ,0]T
CN201910203988.4A 2019-03-18 2019-03-18 旋转轴垂直度误差的误差矢量建模方法 Active CN109933918B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910203988.4A CN109933918B (zh) 2019-03-18 2019-03-18 旋转轴垂直度误差的误差矢量建模方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910203988.4A CN109933918B (zh) 2019-03-18 2019-03-18 旋转轴垂直度误差的误差矢量建模方法

Publications (2)

Publication Number Publication Date
CN109933918A CN109933918A (zh) 2019-06-25
CN109933918B true CN109933918B (zh) 2021-04-27

Family

ID=66987513

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910203988.4A Active CN109933918B (zh) 2019-03-18 2019-03-18 旋转轴垂直度误差的误差矢量建模方法

Country Status (1)

Country Link
CN (1) CN109933918B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112526927B (zh) * 2021-02-18 2021-06-08 成都飞机工业(集团)有限责任公司 面向五轴数控机床旋转轴空间定位误差快速优化补偿方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106959667A (zh) * 2017-04-11 2017-07-18 西南交通大学 一种机床平动轴垂直度误差建模方法
CN107066726A (zh) * 2017-04-11 2017-08-18 西南交通大学 一种数控机床旋转轴垂直度误差建模方法
CN108563186A (zh) * 2018-04-02 2018-09-21 西南交通大学 一种五轴球头铣削几何误差补偿方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106959667A (zh) * 2017-04-11 2017-07-18 西南交通大学 一种机床平动轴垂直度误差建模方法
CN107066726A (zh) * 2017-04-11 2017-08-18 西南交通大学 一种数控机床旋转轴垂直度误差建模方法
CN108563186A (zh) * 2018-04-02 2018-09-21 西南交通大学 一种五轴球头铣削几何误差补偿方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
Accuracy enhancement of five-axis machine tool based on differential motion matrix:Geometric error modeling,identification and compensation;Guoqiang Fu,et al.;《International Journal of Machine Tools and Manufacture》;20141128;第89卷;170-181 *
基于指数积理论和坐标系微分运动关系的数控机床几何误差建模与补偿方法研究;付国强;《中国博士学位论文全文数据库 工程科技I辑》;20190115;B022-86 *

Also Published As

Publication number Publication date
CN109933918A (zh) 2019-06-25

Similar Documents

Publication Publication Date Title
CN108372428B (zh) 五轴机床结构误差自动测量补偿的方法及校正装置
Fu et al. Product of exponential model for geometric error integration of multi-axis machine tools
CN109765848B (zh) 五轴机床运动轴几何误差贡献值计算及影响评估方法
Fu et al. Product-of-exponential formulas for precision enhancement of five-axis machine tools via geometric error modeling and compensation
Peng et al. Total differential methods based universal post processing algorithm considering geometric error for multi-axis NC machine tool
CN109709892B (zh) 一种多轴联动数控机床空间误差在线补偿方法
CN112558547B (zh) 一种五轴数控机床平动轴几何误差补偿数据快速优化方法
CN107421442A (zh) 一种外部测量辅助的机器人定位误差在线补偿方法
Huang et al. Identification of two different geometric error definitions for the rotary axis of the 5-axis machine tools
CN109933920B (zh) 一种旋转轴位置偏差的误差矢量建模方法
CN107066726B (zh) 一种数控机床旋转轴垂直度误差建模方法
Wang et al. A reconfigurable mechanism model for error identification in the double ball bar tests of machine tools
CN110837246A (zh) 五轴数控机床双旋转轴几何误差敏感度分析方法
Ding et al. Novel method for position-independent geometric error compensation of five-axis orthogonal machine tool based on error motion
CN109839920B (zh) 一种五轴机床运动轴灵敏度分析方法
Li et al. Sensitivity analysis of relationship between error motions and machined shape errors in five-axis machining center-Peripheral milling using square-end mill as test case
CN114216456A (zh) 一种基于imu与机器人本体参数融合的姿态测量方法
CN109933918B (zh) 旋转轴垂直度误差的误差矢量建模方法
CN110703689A (zh) 一种数控机床空间误差补偿方法及系统
CN113359609B (zh) 五轴数控机床关键几何误差优化配比补偿方法
Yin et al. Vision-based measurement for decoupling identification of geometric errors of rotating axes for five-axis platform
Lin et al. Modeling and measurement of active parameters and workpiece home position of a multi-axis machine tool
Fan et al. Research on geometric error modeling and compensation method of CNC precision cylindrical grinding machine based on differential motion theory and Jacobian matrix
Zou et al. Error distribution of a 5-axis measuring machine based on sensitivity analysis of geometric errors
Tang et al. A geometric errors identification method for the rotating axis of five-axis welding equipment

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