一种基于固定角向标定及补偿的加工精度提升方法与流程

1.本发明涉及航空航天数控加工技术领域，特别涉及一种基于固定角向标定及补偿的加工精度提升方法。


背景技术：

2.随着航空制造产业的飞速发展，机匣、整体叶盘、叶片等航空零部件结构更加复杂的同时，加工精度要求越来越高，国内早期引入的数控设备大部分已经难以满足新产品的加工精度要求。然而，航空制造企业在有百年发展历程的同时，也面临着设备老旧的问题，航空航天制造企业80％的数控设备在1970～2000年之间采购，并且受限于企业经济发展需求和厂房空间面积，短期内不会进行大规模设备引进或批量设备更换。面对日益激烈的市场竞争压力，如何提高机床加工精度，保证零部件加工质量，降低相关生产成本，加大经济效益，是企业的发展目标。
3.加工中心由于零件制造、安装工艺、使用磨损等原因，每台设备都不可避免一些误差现象存在。由于误差对于设备的各项精度有着严重的影响，对数控立式加工中心的误差分析也在逐步的深入。国内开展机床精度补偿研究的主要有两个方向：一部分学者通过镭射仪检测数控立式加工中心线性定位误差，采用数控系统对数控加工中心几何误差进行补偿；另一部分学者，从减少发热、控制温升、热误差补偿 3个层次对卧式加工中心的热变形与热误差进行补偿。国内的研究对于机床误差的检测及精度的提升具有促进作用，但是对于一般精度或是精度不稳定的五坐标加工中心加工航空发动机复杂零部件时，引发的各向异性误差、多角度误差积累等问题，尚未得到解决。到目前为止，尚没有公开的基于定角向标定及补偿的加工精度提升方法。


技术实现要素：

4.为解决上述技术问题，提出了一种基于固定角向标定及补偿的加工精度提升方法，具体技术方案如下：
5.一种基于固定角向标定及补偿的加工精度提升方法，包括如下步骤：
6.步骤1，机床主轴回到机床零点；
7.步骤2，完成测头的标定及补偿；
8.测头标定指机床通过在线测量系统测量标准特征，计算测量偏差；所述的补偿指将根据测量标准特征计算的测量误差补偿到机床的坐标基准中；
9.步骤3，读入标准的数控程序；
10.步骤4，对程序进行检索，确认程序的最大行程角度范围；
11.所述的最大行程角度范围，指在数控程序中检索最大a角∠a
max
、最小a角∠a
min
，最大b角∠b
max
、最小b角∠b
min
；
12.步骤5，计算程序角向约束；
13.首先，计算出程序的回转中心∠a
cen
、∠b
cen
，∠a
cen
＝(∠a
max
∠a
min
)/2，∠b
cen
＝
(∠b
max
∠b
min
)/2，则程序的角向约束范围∠a＝[∠a
max-∠a
cen
]，∠b＝[∠b
max-∠b
cen
]；
[0014]
步骤6，以标准球的球心为原点，完成基于角向约束的标准球25 个均布点测量；
[0015]
将标准球的球心定为(0，0，0)，球面上∠a、∠b上的点坐标为：
[0016][0017]
其中，0≤∠a《90、0≤∠b《90；
[0018]
步骤7，完成定角向精确配准计算；
[0019]
所述的精确配准，指将实际测量的点坐标模型向理论的点坐标模型进行变换、重合，求解变换过程中移动、旋转的数值；
[0020]
步骤8，将求解的旋转arot和平移trans数值输入到机床的零件加工基准中，完成误差补偿；
[0021]
步骤9，在补偿后的坐标基准下，完成零件数控加工；
[0022]
步骤10，将零件加工基准补偿清零，用于下一个叶片定角向分析及误差补偿。
[0023]
所述的一种基于固定角向标定及补偿的加工精度提升方法，其优选方案为，步骤2所述的标准特征，包括标准孔特征、标准球特征；标准球特征的测量，指对半球表面进行25点均布测量；标准孔特征的测量，指对工作台孔内径进行 x向、-x向、 y向、-y向，定向测量；在线测量时，同一点在测头0度、测头180度条件下测量2次，计算每一点的平均值，保证每一点测量的准确性。
[0024]
所述的一种基于固定角向标定及补偿的加工精度提升方法，其优选方案为，步骤6所述的25个均布点测量，指主轴在∠a
cen
、∠b
cen
的角向下，在限定角向的约束范围内均匀分布25点，计算出点坐标，并且完成在线测量；其中，25点的点坐标及点位矢量为：
[0025]
原点：
[0026]
p
(0，0)
(0，0，r，0，0)
[0027]
x轴不变，绕y轴旋转，测量x轴轴线上的6个点坐标：
[0028]
p
(-3，0)
(-rsin∠b，0，rcos∠b，0，-∠b)
[0029][0030][0031][0032][0033][0034]
y轴不变，绕x轴旋转，测量y轴轴线上的6个点坐标：
[0035]
p
(0，-3)
(0，-rsin∠a，rcos∠a，-∠a，0)
[0036]
[0037][0038][0039][0040][0041]
第一象限内，x、y轴联动，测量3个点的点坐标：
[0042][0043][0044][0045]
第二象限内，x、y轴联动，测量3个点的点坐标：
[0046]
[0047][0048][0049]
第三象限内，x、y轴联动，测量3个点的点坐标：
[0050][0051][0052][0053]
第四象限内，x、y轴联动，测量3个点的点坐标：
[0054][0055][0056][0057]
所述的一种基于固定角向标定及补偿的加工精度提升方法，其优选方案为，步骤7所述的定角向精确配准计算，指以25点的质心作为相对坐标系原点，将实际测量的25点质心和主方向行旋转arot和平移 trans，与理论25点质心和主方向重合，完成定角向配准；
[0058]
其中，arot为3
×
3的旋转矩阵，trans为3
×
1的平移矩阵；act为 25点实测值，质心定义为p
act
，主方向为(v1，v2，v3)；org为25点理论值，质心定义为p
org
，主方向为(u1，u2，u3)；p
acti
为点act中为 i的点，p
orgi
为点云org中为i的点；则理论点的质心及协方差计算公式为：
[0059]
p
orgi
＝(u
ix
,u
iy
,u
iz
)，则
[0060]
计算点云协方差矩阵
[0061]
实测点的质心及协方差计算公式为：
[0062]
p
acti
＝(v
ix
,v
iy
,v
iz
)，则
[0063]
计算点云协方差矩阵
[0064]
根据协方差得到的每个点的特征值，从大到小排序，选择三个特征向量，作为点云的主方向，构成两个点云的协方差矩阵t
act
'，t
org
'：
[0065][0066]
旋转矩阵arot：arot＝t
org
'
×
t
act
'-1
[0067]
平移矩阵trans：trans＝p
org-arot
×
p
act
。
[0068]
本发明的有益效果：
[0069]
本发明首次提出基于固定角向标定及补偿的加工精度提升方法，成功将相关技术应用于大型风扇钛合金整体叶盘的自适应数控加工中，在原有设备0.06mm的在线测量精度条件下，实现线性摩擦焊整体叶盘焊后叶片0.03mm接刀痕控制，完成大型风扇钛合金整体叶盘 20个叶片凸台部位精加工，填补了基于定角向标定及补偿的加工精度提升的技术空白。该方法的实现，解决了机床在加工大型发动机零部件各向特征的精度偏差问题、周向特征积累问题，有较强的通用型和实用性，在为企业提升核心创新能力和研发效率的同时创造巨大的经济效益。
[0070]
本方法经实际加工验证，按照本发明对提取加工程序五轴联动范围、基于角向约束计算标准球校准程序、定角向误差分析、叶片加工基准补偿，在原有设备0.06mm的在线测量精度条件下，实现线性摩擦焊整体叶盘焊后叶片0.03mm接刀痕控制。按照单台设备精度提升后，投入大型零件精密切削产生的经济效益计算，每当1台设备精度提升后，每年创造的效益按照“效益＝单台设备全年生产时间*该状态设备现有数量*高精度设备每小时占用费用”计算，即(285天/年*10h/ 天)/台*20台*20元/h＝101.4万/年。
附图说明
[0071]
图1为本发明的一种基于固定角向标定及补偿的加工精度提升方法流程图；
[0072]
图2为本发明中数控程序在最大角向联动环境下，机床上标准球 p点测量点示意图；
[0073]
图中，球面上∠a为机床绕x轴旋转的角度，∠b为机床绕y轴旋转的角度；
[0074]
图3为本发明标准球25个均布点测量示意图；
[0075]
图4为本发明定角向移动及旋转精确配准示意图；
[0076]
图中，act为25点实测值，质心定义为p
act
，主方向为(v1，v2， v3)；org为25点理论值，质心定义为p
org
，主方向为(u1，u2，u3)， arot为3
×
3的旋转矩阵，trans为3
×
1的平移矩阵；
[0077]
图5为本发明测头的标定及补偿示意图；
[0078]
图中，tool_number为调用刀具库中的刀具号，9999为机床中的测头；r8为测头基准x向误差补偿变量、r9为测头基准y向误差补偿变量、r7为测头半径补偿、r19为测头基准z向误差补偿；
[0079]
图6为本发明数控程序最大行程角度范围检索的逻辑关系示意图；
[0080]
图中，n为数控程序代码的总行数，i代表读取到数控代码的当前行数，a代表数控程序中a轴的角度，b代表数控程序中b轴的角度，∠a
max
为数控程序中最大a角，∠a
min
∠a
max
为数控程序中最小a 角，∠b
max
∠a
max
为数控程序中最大b角，∠b
min
∠a
max
为数控程序中最小b角；
[0081]
图7为本发明基于固定角向最大范围约束的精确配准计算结果；
[0082]
图8为本发明的机床基准误差补偿示意图。
具体实施方式
[0083]
本发明采用基于固定角向标定及补偿的加工精度提升方法，成功原有设备0.06mm的在线测量精度条件下，实现线性摩擦焊整体叶盘焊后叶片0.03mm接刀痕控制，完成大型风扇钛合金整体叶盘20个叶片凸台部位精加工；本专利以公司五轴机床为例，结合附图1-8和实施过程对本发明作进一步说明。
[0084]
一种基于固定角向标定及补偿的加工精度提升方法，包括如下步骤：
[0085]
1)机床主轴回到机床原点；
[0086]
2)完成测头的标定及补偿；
[0087]
针对d25标准球，编制出25点覆盖半球面的3轴测量程序，从机床刀具库中调出测头及探针，在机床基准下完成标准球的在机测量，分析测量偏差：
[0088]
x向偏差-0.0003mm，y向偏差-0.0034mm，z向偏差0.2685mm
[0089]
由于z向偏差过大，需要进行补偿，首先将计算出的x向、y 向、z向的偏差反向输入到机床基准中进行复测，分析标准球各项偏差：
[0090]
x向偏差-0.0014mm，y向偏差-0.0036mm，z向偏差0.0015mm
[0091]
标准球各向精度良好，满足测量条件，可以用于在机测量；
[0092]
3)读入标准的数控程序；
[0093]
读入数控程序，包括apt文件、g-code文件；
[0094]
4)完成程序检索；
[0095]
按照图6，通过压栈的方式，完成程序的最大行程角度范围的检索。以机床加工某型号航空发动机整体叶盘第n片叶片的数控程序为例，通过程序检索，该段程序的最大行程角度范围为：
[0096]
最大a角∠a
max
＝29.9981
°
；
[0097]
最小a角∠a
min
＝-29.998
°
；
[0098]
最大b角∠b
max
＝55.7293
°
；
[0099]
最小b角∠b
min
＝-30.6237
°
；
[0100]
5)计算程序角向约束范围；
[0101]
按照权利要求书中步骤5所述的方法，首先计算出程序的回转中心：
[0102]
∠a
cen
＝(∠a
max
∠a
min
)/2＝0
°
；
[0103]
∠b
cen
＝(∠b
max
∠b
min
)/2＝15.6028
°
；
[0104]
则程序的角向约束范围：
[0105]
∠a＝[∠a
max-∠a
cen
]＝30
°
；
[0106]
∠b＝[∠b
max-∠b
cen
]＝45
°
；
[0107]
即机床主轴在a角、b角度的前提下，以a0
°
(30
°
)的角度范围联动；
[0108]
6)计算标准球上的检测点坐标；
[0109]
以标准球的球心为原点，完成基于角向约束的标准球25个均布点测量；将标准球的球心定为(0，0，0)，∠a＝30
°
、∠b＝45
°
时，球面 25点的点坐标为：
[0110]
原点：
[0111]
p
(0，0)
(0，0，12.5，0，0)
[0112]
x轴不变，绕y轴旋转，测量x轴轴线上的6个点坐标：
[0113]
p
(-3，0)
(-8.84，0，8.84，0，-45)
[0114]
p
(-2，0)
(-6.25，0，10.825，0，-30)
[0115]
p
(-1，0)
(-3.235，0，12.074，0，-15)
[0116]
p
(1，0)
(3.235，0，12.074，0，15)
[0117]
p
(2，0)
(6.25，0，10.825，0，30)
[0118]
p
(3，0)
(8.84，0，8.84，0，45)
[0119]
y轴不变，绕x轴旋转，测量y轴轴线上的6个点坐标：
[0120]
p
(0，-3)
(0，-6.25，10.825，-30，0)
[0121]
p
(0，-2)
(0，-4.275，11.746，-20，0)
[0122]
p
(0，-1)
(0，-2.171，12.31，-10，0)
[0123]
p
(0，1)
(0，2.171，12.31，10，0)
[0124]
p
(0，2)
(0，4.275，11.746，20，0)
[0125]
p
(0，3)
(0，6.25，10.825，30，0)
[0126]
第一象限内，x、y轴联动，测量3个点的点坐标：
[0127]
p
(1，1)
(3.189，2.099，11.903，10，15)
[0128]
p
(2，2)
(5.961，3.758，10.325，20，30)
[0129]
p
(3，3)
(8.183，6.25，8.183，30，45)
[0130]
第二象限内，x、y轴联动，测量3个点的点坐标：
[0131]
p
(-1，1)
(-3.189，2.099，11.903，10，-15)
[0132]
p
(-2，2)
(-5.961，3.758，10.325，20，-30)
[0133]
p
(-3，3)
(-8.183，6.25，8.183，30，-45)
[0134]
第三象限内，x、y轴联动，测量3个点的点坐标：
[0135]
p
(-1，-1)
(-3.189，-2.099，11.903，-10，-15)
[0136]
p
(-2，-2)
(-5.961，-3.758，10.325，-20，-30)
[0137]
p
(-3，-3)
(-8.183，-6.25，8.183，-30，-45)
[0138]
第四象限内，x、y轴联动，测量3个点的点坐标：
[0139]
p
(1，-1)
(3.189，-2.099，11.903，-10，15)
[0140]
p
(2，-2)
(5.961，-3.758，10.325，-20，30)
[0141]
p
(3，-3)
(8.183，-6.25，8.183，-30，45)
[0142]
7)完成定角向精确配准计算；
[0143]
按照权利要求书中步骤7所述的方法，完成基于定角向约束的精确配准计算；
[0144]
所述的精确配准，指将实际测量的点坐标模型向理论的点坐标模型进行变换、重合，求解变换过程中移动、旋转的数值；
[0145]
通过配准，计算出旋转arot和平移trans数值，如图7所示：
[0146][0147]
8)完成误差补偿；
[0148]
将求解的旋转arot和平移trans数值输入到机床的零件加工基准中，完成误差补偿，如图8所示；
[0149]
9)完成数控加工；
[0150]
在补偿后的坐标基准下，完成零件数控加工；
[0151]
10)零件加工基准补偿清零；
[0152]
将零件加工基准补偿清零，用于下一个叶片定角向分析及误差补偿。