一种摇篮式五轴机床双旋转轴位置无关误差的辨识方法技术

一种摇篮式五轴机床双旋转轴位置无关误差的辨识方法。本发明专利技术公开了一种关于任意转轴位置摇篮式五轴机床双转台位置无关误差的检测方法，其特征在于，针对任意转台轴线位置结构的摇篮式五轴机床利用球杆仪设计了通用的双转台联动轨迹，根据单位对偶四元数的全局表征形式与误差定义建立了工件端和刀具端间的运动学误差模型，通过对运动轨迹的均分化处理解决球杆仪运动与其采集速率不同步的问题，基于最小二乘算法并结合球杆仪实验数据对简化后的误差运动学模型进行解耦以辨识双转台的位置无关结合误差。本发明专利技术方法通过单条实验轨迹进行位置无关几何误差辨识，适用于各类任意旋转轴位置结构的摇篮式五轴机床，方法简单快捷，实用性好。实用性好。实用性好。

【技术实现步骤摘要】
一种摇篮式五轴机床双旋转轴位置无关误差的辨识方法


[0001]本专利技术属于数控机床检测
，特别涉及一种关于任意转轴位置结构的摇篮式五轴机床双旋转轴位置无关几何误差的辨识方法。
技术背景
[0002]相较于三轴机床，五轴机床通过转动旋转轴可以调整刀具与加工平面间的相对方向，使其具有更好的灵活性和更高的效率。摇篮式五轴机床是小型复杂曲面零件加工中应用最为广泛的机床类型，其根据双旋转轴线位置的不同可分为相交结构和非相交结构。相交共面结构作为非相交结构的特殊情况，是众多研究和实验方法的研究对象，然而这些方法大多并不适用于双旋转轴轴线非相交结构。
[0003]针对机床位置无关几何误差的辨识主要分为机床建模和误差测量两部分。现有的建模方法理论有齐次变换矩阵、旋量理论和指数积公式等，但以上方法中在运动学模型建立过程中仍局限于误差矩阵间的相乘，一定程度上导致计算的效率较低。误差测量方面使用的仪器包括激光干涉仪、激光跟踪仪、球杆仪和R
‑
test等仪器。其中球杆仪由于便携性、结构简单、使用方便等优点，被广泛应用于误差测量。但是由于球杆仪的测量原理限制，在多轴同步运动过程中会存在球杆仪合成运动与球杆仪采样频率不同步的问题，很少有一种通用且简便的测量方法适用于任意旋转轴线位置的摇篮式五轴机床双旋转轴位置无关几何误差的测量。

技术实现思路

[0004]为解决上述问题，本专利技术的目的在于提出了一种关于任意旋转轴线位置结构的摇篮式五轴机床双旋转轴位置无关几何误差的辨识方法，通过基于单位对偶四元数建立全局坐标系下的五轴机床运动学误差模型，并结合球杆仪提出适用于任意旋转轴线位置结构的摇篮式五轴机床的双旋转轴同步协调运动轨迹，以更为简洁、快捷准确地辨识A轴和C轴的位置无关几何误差。具体方法步骤如下：
[0005]步骤1、基于单位对偶四元数变换原理对旋转轴位置无关几何误差进行重新定义表征，并建立五轴机床理想和实际的运动学模型。
[0006]步骤2、根据摇篮式五轴机床双旋转轴间的位置偏移并结合理想的运动学模型设计旋转轴A轴和C轴的同步协调运动轨迹。
[0007]步骤3、建立轨迹均分算法解决双旋转轴A轴和C轴同步协调运动过程中合运动与球杆仪采样频率间的不同步问题。
[0008]步骤4、简化实际的运动学误差模型，结合最小二乘算法进行双旋转轴位置无关几何误差的解耦。
[0009]步骤1中根据单位对偶四元数变换原理重新定义了摇篮式五轴机床双旋转轴的位置无关几何误差，并基于此建立了五轴机床理想和实际的运动学模型，包括步骤：
[0010]步骤1.1、基于单位对偶四元数，给出其旋量运动的指数表示形式：
[0011][0012]其中对偶角度和对偶向量分别表示旋量运动的运动过程和所绕螺旋轴的位姿，θ表示绕螺旋轴的旋转角度，d表示沿螺旋轴的位移距离，n表示螺旋轴的单位方向向量，m＝p
×
n表示螺旋轴的矩，p是螺旋轴上任意一点。
[0013]步骤1.2、根据泰勒展开式将等式(1)展开为
[0014][0015]步骤1.3、根据五轴机床的运动轴类型，由式(2)分别给出线性轴和旋转轴运动的单位对偶四元数表征形式：
[0016][0017]步骤1.4、基于式(3)，将梧州机床旋转轴的名义的和实际的运动轴线间的位置无关几何误差变换如图2所示。以旋转轴C轴为例，其理想轴线在基坐标系b下通过式(1)中的对偶向量(也被称为Pl
ü
cker line)表示为
[0018][0019]步骤1.5、对应得到轴线与基坐标系的XOY平面上的交点坐标为
[0020][0021]步骤1.6、将名义轴线和实际轴线间的位置无关几何误差被转化一个为绕位于 XOY平面Pl
ü
cker line的微量角度的旋转变换
[0022][0023]以及一个XOY平面上的微量位移变换
[0024][0025]其中表示微量位移向量。
[0026]步骤1.7、基于最小角理论，由于角度参数非常小，式(6)简化为
[0027][0028]其中表示微量旋转向量。
[0029]步骤1.8、名义和实际轴线间的变换过程总结表示为
[0030][0031]其中形如表示单位对偶四元数的乘法运算，是其伴随表示形式。
[0032]以旋转C轴为例的位置无关几何误差参数可以表征为
[0033][0034]步骤1.9、基于如图3所示的五轴机床结构式意图和机床运动链示意图，将机床的基坐标系b设置与C轴旋转轴中心重合，各个运动轴在基坐标系下初始的Pl
ü
cker line位姿分别为
[0035][0035][0035][0035][0035][0036]其中y
ab
和z
ab
分别表示A轴在基坐标系下Y和Z方向的位置偏移。
[0037]步骤1.10、基于单位对偶数的乘法运算及其误差表征，建立工件在基坐标系下的名义运动学模型为
[0038][0039]其中其中分别表示刀具链中第i个运动轴在基础坐标系b中运动驱动量和Pl
ü
cker线位姿表征，表示刀具在基坐标系下的位置单位对偶四元数。
[0040]步骤1.11、建立刀具在基坐标系下的名义运动学模型为
[0041][0042]其中其中分别表示工件链中第j个运动轴在基础坐标系b中运动驱动
量和Pl
ü
cker线位姿表征，表示工件在基坐标系下的位置单位对偶四元数。
[0043]步骤1.12、建立工件在刀具坐标系下名义下的运动学模型表示为
[0044][0045]步骤1.13、结合旋转轴位置无关几何误差定义，建立五轴机床实际的运动学误差模型为
[0046][0047]步骤2、针对双旋转轴线位置不相交的摇篮式五轴机床，结合理想的运动学模型设计了通用的A、C旋转轴的球杆仪同步协调运动轨迹。球杆仪两端分别吸附在主轴和工作台上，具体步骤如下：
[0048]步骤2.1、根据式(14)的名义运动学模型，在不考虑线性轴的运动情况下，刀具工具杯保持静止，得到名义的工件在刀具坐标系下的位置变化为
[0049][0050]进而得到相对坐标表示为
[0051][0052]其中s θ
c
＝sin θ
c
，c θ
c
＝cos θ
c
，s θ
a
＝sin θ
a
，c θ
a
＝cos θ
a
；t
wb
(0)＝[x
wb y
wb z
wb
]T
和 t
tb
(0)＝[x
tb y
tb z
tb
]T
分别表示工件工具杯和刀具工具杯在基坐标系下的位置向量。为了保证球杆仪在运动过程中不掉落，本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种摇篮式五轴机床双旋转轴位置无关误差的辨识方法，其特征在于，包括如下步骤：步骤1、基于单位对偶四元数变换原理重新表征旋转轴的位置无关几何误差，并建立五轴机床理想和实际的运动学模型。步骤2、根据摇篮式五轴机床双旋转轴线间的位置偏移，并结合理想的运动学模型设计A轴和C轴的同步协调运动轨迹。步骤3、建立轨迹均分算法，解决双旋转轴A轴和C轴同步协调运动过程中的合运动与球杆仪采样频率间的不同步问题。步骤4、简化实际的运动学误差模型，结合最小二乘算法进行双旋转轴位置无关几何误差的解耦。2.根据权利要求1所述的关于摇篮式五轴机床双旋转轴位置无关误差的辨识方法，其特征在于，在步骤1中，基于单位对偶四元数变换理论，将五轴机床每个旋转轴的位置无关几何误差重新定义表征为两项微量旋转参量和两项微量位移参量。误差定义过程步骤如下：步骤1.1、基于单位对偶四元数，给出其旋量运动的指数表示形式：其中对偶角度和对偶向量分别表示旋量运动的运动过程和所绕螺旋轴的位姿，θ表示绕螺旋轴的旋转角度，d表示沿螺旋轴的位移距离，n表示螺旋轴的单位方向向量，m＝p
×
n表示螺旋轴的矩，p是螺旋轴上任意一点。步骤1.2、根据泰勒展开式将等式(1)展开为步骤1.3、根据五轴机床的运动轴类型，由式(2)分别给出线性轴和旋转轴运动的单位对偶四元数表征形式：步骤1.4、基于式(3)，如图2所示，以旋转轴C轴为例，建立名义和实际运动轴线位置无关几何误差间的变换关系。其名义轴线在基坐标系b下表示为式(1)中对偶向量(也被称为Pl
ü
cker line)的形式步骤1.5、对应得到轴线与基坐标系的XOY平面上的交点坐标为
步骤1.6、将名义轴线和实际轴线间的位置无关几何误差被转化一个为绕位于XOY平面Pl
ü
cker的微量角度的旋转变换以及一个XOY平面上的微量位移变换其中表示微量位移向量。步骤1.7、基于最小角理论，由于角度参数非常小，将式(6)简化为其中表示微量旋转向量。步骤1.8、名义和实际轴线间的变换过程总结表示为其中形如表示单位对偶四元数的乘法运算，是其伴随表示形式。以旋转C轴为例的位置无关几何误差参数可以表征为对应旋转C轴的位置无关几何误差参数表征，给出A轴和B轴的误差定义如表1所示。表1 旋转轴位置无关几何误差表示。3.根据权利要求1所述的关于摇篮式五轴机床双旋转轴位置无关误差的辨识方法，其特征在于，在步骤1当中，基于单位对偶数的乘法运算及其误差表征建立五轴机床名义和实际的运动学模型，包括步骤：步骤1.10、基于单位对偶数的乘法运算及其误差表征，建立工件在基坐标系下的名义运动学模型为
其中分别表示刀具链中第i个运动轴在基础坐标系b中运动驱动量和Pl
ü
cker线位姿表征，表示刀具在基坐标系下的位置单位对偶四元数。步骤1.11、建立刀具在基坐标系下的名义运动学模型为其中分别表示工件链中第j个运动轴在基础坐标系b中运动驱动量和Pl
ü
cker线位姿表征，表示工件在基坐标系下的位置单位对偶四元数。步骤1.12、建立工件在刀具坐标系下名义下的运动学模型表示为步骤1.13、结合旋转轴位置无关几何误差定义，建立五轴机床实际的运动学误差模型为。4.根据权利要求1所述的关于摇篮式五轴机床双旋转轴位置无关误差的辨识方法，其特征在于，在步骤2当中，结合理想的运动学模型设计了通用的A、C旋转轴的球杆仪同步协调运动轨迹，包括步骤：步骤2.1、根据式(14)的名义运动学模型，在不考虑线性轴的运动情况下，刀具工具杯保持静止，得到名义的工件在刀具坐标系下的位置变化为进而得到相对坐标表示为其中sθ
c
＝sinθ
c
，cθ
c
＝cosθ
c
，sθ
a
＝sinθ
a
，cθ
a
＝cosθ
a
；t
wb
(0)＝[x
wb y
wb z
wb
]
T
和t
tb
(0)＝[x
tb y
tb z
tb
]
T
分别表示工件工具杯和刀具工具杯在基坐标系下的位置向量。为了保证球杆仪在运动过程中不掉落，式(17)表示的相对坐标各方向分量需要满足
步骤2.2、考虑使用的DMU 85摇篮式五轴机床双旋转轴线间的位置偏移，为了简化实验安装和误差模型，使用长度为L
dbb
＝100mm的球杆仪进行实验，分别对球杆仪刀具杯和工具杯的位置进行以下设置：设置1：当A轴和C轴转角都为0
°
时，工件工具杯和刀具工具杯均位于基坐标系的z＝0平面上，即z
wb
＝z
tb
＝0；设置2：将工件工具杯和刀具工具杯分别设置在基坐标系x轴和y轴上，即t
wb

【专利技术属性】
技术研发人员：蒋晓耕，王浩，郭瑞东，
申请(专利权)人：天津工业大学，
类型：发明
国别省市：