一种数控机床位置相关几何误差项的建模方法技术

本发明专利技术公开了一种数控机床位置相关几何误差项的建模方法，属于数控机床加工技术领域。包括：根据数控机床位置相关几何误差项性质，确定位置相关几何误差项n次多项式形式；根据数控机床位置相关几何误差项离散数据，采用最小二乘法拟合原理得到n次多项式系数计算矩阵M，计算n次多项式系数；结合数控机床位置相关几何误差项离散数据计算从1次到m次的一系列多项式，最高次数m设定为几何误差项离散数据数目的一半；根据回归分析F检验方法计算1次到m次的一系列多项式的F值；比较找出最大F值，选择最大F值对应的多项式为数控机床位置相关几何误差项的优化多项式模型。主要用于实现机床综合误差建模的自动化和程序化。

A method for modeling the position dependent geometric error of NC machine tools

The invention discloses a method for modeling position related geometric error items of numerically controlled machine tools, which belongs to the technical field of numerical control machine tools processing. Including: according to the position of NC machine tool geometric error properties, determine the location of geometric error n polynomial; according to the position of CNC machine tool geometric error of the discrete data fitting principle n polynomial coefficient matrix M by least square method, calculating n polynomial coefficient; calculation from 1 to a series of M polynomials the combination of CNC machine tool geometric error location of discrete data, the maximum number of M is set to half the number of errors of geometric discrete data; regression analysis based on the F test method 1 times to a series of polynomial m of F-measure; compared to find the largest F-measure, select the maximum F-measure corresponding polynomial for the position of NC machine tool geometric errors optimal polynomial model. It is mainly used to automate and program the integrated error modeling of machine tools.

【技术实现步骤摘要】
一种数控机床位置相关几何误差项的建模方法
本专利技术涉及数控机床加工控制
，尤其涉及数控机床位置相关几何误差项优化多项式建模技术。
技术介绍
高精度加工是数控机床发展的必然趋势，几何误差是数控机床在零件加工过程中的主要误差源之一，几何误差重复性高、系统性好、易测量的特点使得几何误差补偿成为提高机床精度的重要措施之一。几何误差建模是进行误差补偿的基础，其中建立几何误差项的数学模型是误差建模和补偿的重要环节，直接关系到误差补偿精度。根据数控机床几何误差项性质定义可将其分为两类：位置不相关几何误差和位置相关几何误差。位置相关几何误差反映了运动轴的运动精度，误差值会随着运动轴进给量的变化而变化。数控机床位置相关几何误差是通过间接测量方法或者直接测量方法来测量运动轴轴线上或者机床工作空间中均匀分布的测量点处的几何误差值并进行辨识得到的，所以几何误差项数据是离散的。一般来说，位置相关几何误差项数值呈现非线性变化且无规律。多项式是常用的几何误差项表示形式，多项式次数不仅关系到模型精度，同时对后续计算效率有很大影响，需要进行进一步确定。Lee等人根据旋转轴基本几何误差项的性质，采用C1类连续函数进行建模，然后通过误差测量辨识得到模型中的参数(参见LeeKI,LeeDM,YangSH(2012)Parametricmodelingandestimationofgeometricerrorsforarotaryaxisusingdoubleball-bar.TheInternationalJournalofAdvancedManufacturingTechnology,62(5-8):741-750)。但是该方法的拟合精度有待商榷，且多项式的次数等需要进一步研究。Fan等人采用正交多项式对基本几何误差项进行建模，将多项式拟合转化为多线性回归问题(参见FanK,YangJ,YangL(2013)Orthogonalpolynomials-basedthermallyinducedspindleandgeometricerrormodelingandcompensation.TheInternationalJournalofAdvancedManufacturingTechnology65(9-12):1791-1800)。但是需要查正交多项式F表和方差分析表来确定系数。这样几何误差数据数目的不确定性增加了查表的复杂性，无法实现几何误差元素建模的程序化。
技术实现思路
本专利技术的目的是提供一种数控机床位置相关几何误差项的建模方法，它能有效地解决建立几何误差项的数学表达式的问题。本专利技术的目的是通过以下技术方案来实现的：一种数控机床位置相关几何误差项的建模方法，包括如下步骤：步骤1、根据数控机床位置相关几何误差项的性质，确定位置相关几何误差项n次多项式形式为；其中est表示机床运动轴t在s方向上的位置相关几何误差项；t表示运动轴t运动量，t＝x、y、z、α、β、γ；s＝x、y、z；e＝δ、ε，δ表示线性误差，ε表示角度误差；n表示多项式次数；fn表示n次多项式；bj表示多项式j次项系数。步骤2、根据某项数控机床位置相关几何误差项离散数据，采用最小二乘法拟合原理得到n次多项式系数计算矩阵M，n次多项式系数计算公式表示为：B＝M-1·Y其中M表示n次多项式系数计算矩阵，B表示n次多项式系数矩阵，Y表示位置相关几何误差项离散数据矩阵。步骤3、根据步骤2，结合数控机床位置相关几何误差项的某项离散数据计算得到从1次到m次的一系列多项式；其中m表示一系列多项式中最高次数；步骤4、根据回归分析F检验方法，结合数控机床位置相关几何误差项离散数据，计算步骤3中得到的1次到m次的一系列多项式的F值；步骤5、比较1次到m次的一系列多项式的F值，找出最大F值，选择最大F值对应的多项式为该数控机床位置相关几何误差项的优化多项式模型。作为优选，所述步骤1中数控机床位置相关几何误差项的性质为误差数值随着运动轴位置变化而变化，在运动轴零位置处位置相关几何误差项数值为零，从而确定位置相关几何误差项多项式形式中常数项为零。作为优选，所述步骤2中采用最小二乘法拟合原理计算n次多项式系数的具体方法为：步骤2.1、最小二乘法采用偏差平方和最小的原则，偏差平方和表示为：其中Δ2表示偏差平方和；fi表示机床位置相关几何误差项第i个离散数据，ti表示第i个离散数据对应的自变量；fn(ti)表示ti对应的拟合多项式计算值，表示第i个离散数据的偏差平方；p表示数控机床位置相关几何误差项离散数据数目；步骤2.2、根据偏差平方和表达式依次对bj求偏导数得到等式：步骤2.3、对等式进行简化后得到n次多项式系数与数控机床位置相关几何误差项离散数据之间关系，表示为：步骤2.4、将n次多项式系数与数控机床位置相关几何误差项离散数据之间关系转换为矩阵形式：步骤2.5、建立n次多项式系数计算矩阵及位置相关几何误差项离散数据矩阵，得到n次多项式系数计算公式：B＝M-1·Y其中作为优选，所述步骤3中结合某项数控机床位置相关几何误差项离散数据计算得到从1次到m次的一系列多项式时，最高次数m的值小于等于8。设定为几何误差项离散数据数目的一半，如果离散数据数目的一半大于8，m值设定为8，即多项式次数不高于8。拟合多项式的表示形式为：式中，m的值小于等于8。作为优选，所述步骤4中根据回归分析F检验方法计算1次到m次的一系列多项式的F值的方法为：其中Fn表示n次多项式的F值，表示机床位置相关几何误差项离散数据的平均值；作为优选，所述数控机床包括三轴数控机床、五轴数控机床等多轴数控机床。作为优选，本专利技术适合于五轴数控机床的30项位置相关几何误差项。作为优选，本专利技术适合于三轴数控机床的18项位置相关几何误差项。作为优选，所述数控机床位置相关几何误差项数据由激光干涉仪测量得到。作为优选，所述数控机床位置相关几何误差项数据由球杆仪测量得到。作为优选，所述数控机床位置相关几何误差项数据由激光跟踪仪测量得到。本专利技术与现有技术相比，具体的有益效果是：本专利技术得到的优化多项式精度高，符合位置相关几何误差性质，且鲁棒性强，适应于不同位置相关几何误差项，可实现建模的自动化和程序化，提高建模的通用性，可实现误差建模的自动化和程序化。附图说明图1为本专利技术数控机床位置相关几何误差项优化多项式建模方法流程图；图2为某五轴数控机床的结构示意图；图3为本专利技术在五轴数控机床X轴x方向线性误差δxx的离散数据；图4为本专利技术在五轴数控机床X轴x方向线性误差δxx一系列多项式F值分布图；图5为本专利技术在五轴数控机床X轴x方向线性误差δxx的优化多项式模型；图6a为本专利技术的数控机床Z轴的位置相关线性几何误差离散数据与优化多项式模型；图6b为本专利技术的数控机床Z轴的位置相关角度几何误差离散数据与优化多项式模型；具体实施方式下面结合附图和具体实施例对本专利技术做进一步说明。附图1所示为本专利技术数控机床位置相关几何误差项的建模方法流程图。附图2所示为某五轴数控机床的结构示意图，附图3所示为X轴x方向线性误差δxx的离散数据，以附图2中机床X运动轴x方向的线性几何误差δxx为例来阐述本专利技术数控机床位置相关几何误差项建模方法。步骤1、根据数控机床位置相关几何误差本文档来自技高网...

【技术保护点】
一种数控机床位置相关几何误差项的建模方法，包括如下步骤：步骤1、根据数控机床位置相关几何误差项的性质，确定位置相关几何误差项n次多项式形式为；

【技术特征摘要】
2016.10.08 CN 20161087698001.一种数控机床位置相关几何误差项的建模方法，包括如下步骤：步骤1、根据数控机床位置相关几何误差项的性质，确定位置相关几何误差项n次多项式形式为；其中，est表示机床运动轴t在s方向上的位置相关几何误差项；t表示运动轴t运动量，t＝x、y、z、α、β、γ；s＝x、y、z；e＝δ、ε，δ表示线性误差，ε表示角度误差；n表示多项式次数；fn表示n次多项式；bj表示多项式j次项系数；步骤2、根据数控机床位置相关几何误差项离散数据，采用最小二乘法拟合原理得到n次多项式系数计算矩阵M，n次多项式系数计算公式表示为：B＝M-1·Y其中，M表示n次多项式系数计算矩阵，B表示n次多项式系数矩阵，Y表示位置相关几何误差项离散数据矩阵；步骤3、根据步骤2，结合数控机床位置相关几何误差项离散数据计算得到从1次到m次的一系列多项式，其中m表示一系列多项式中最高次数；步骤4、根据回归分析F检验方法，结合数控机床位置相关几何误差项离散数据，计算步骤3中得到的1次到m次的一系列多项式的F值；步骤5、比较1次到m次的一系列多项式的F值，找出最大F值，选择最大F值对应的多项式为该数控机床位置相关几何误差项的优化多项式模型。2.根据权利要求1所述的一种数控机床位置相关几何误差项的建模方法，其特征在于：所述步骤1中数控机床位置相关几何误差项的性质，在运动轴零位置处的位置相关几何误差项数值为零，从而确定位置相关几何误差项多项式形式中常数项为零。3.根据权利要求1所述的一种数控机床位置相关几何误差项的建模方法，其特征在于：所述步骤2中采用最小二乘法拟合原理计算n次多项式系数的具体方法为：步骤2.1、最小二乘法采用偏差平方和最小的原则，偏差平方和表示为：其中，Δ2表示偏差平方和；fi表示数控机床位置相关几何误差项第i个离散数据，ti表示第i个离散数据对应的自变量；fn(ti)表示ti对应的拟合多项式计算值，表示第i个离散数据的偏差平方；p表示数控机床位置相关几何误差项离散数据数目；步骤2.2、根据偏差平方和表达式依次对bj求偏导数得到等式：步骤2.3、对等式进行简化后得到n次多项式系数与数控机床位置相关几何误差项离散数据之间关系，表示为：

【专利技术属性】
技术研发人员：付国强，高宏力，宋兴国，曹中清，
申请(专利权)人：西南交通大学，
类型：发明
国别省市：四川,51