机床空间加工误差的非确定性描述及预测方法
【专利摘要】机床空间加工误差的非确定性描述及预测方法，是属于机床精度设计领域。首先，根据多体系统理论建立起机床的误差模型，在误差模型的基础上，对误差项进行合理的削减至三个方向的“当量误差”。当量误差中也同样存在着非确定性的波动，在本项发明中，加工平面时所拥有随机波动可以根据随机过程理论进行描述和预测。波动的范围也应该被限制在一定范围内；此外对加工误差波动有较大影响的关键误差项会被甄别出来，根据得到的结论，提出一些对于机床零件进行改善的地方。改善后的试验机，可以从测量结果中清晰的看到，精度的提升及波动范围的减小。这对于精密及超精密加工有着至关重要的指导意义。
【专利说明】机床空间加工误差的非确定性描述及预测方法
【技术领域】
[0001]本发明是一种关于机床空间加工误差的非确定性描述及预测方法，属于机床精度设计领域。
【背景技术】
[0002]随着科学技术和社会经济的快速发展，数控机床在现代加工制造和高性能装备制造的重要组成部分。如何更好的提高数控机床的加工精度成为国内外学术界的交流热点。影响机床加工精度的因素有很多，例如说几何误差、压力变形误差、热误差和动态误差等等。机床几何误差是影响加工精度的最重要部分，几乎占到所有误差的30% -40%，特别在精密及超精密的加工情况下。又因为几何误差很少受外界环境的影响，从而建立几何误差丰吴型并进行分析对于提闻加工精度有着很深的意义。
[0003]机床的几何误差最主要来源于其导轨的制造精度还有安装精度及本身的直线度等误差。由于在安装过程中几何误差存在着一定的随机性，所以在不同的位置也会存在着一定的波动。几何误差都可以被分为两部分，一个确定量部分和符合某种概率特征分布的围绕着确定值波动的随机量部分。定义确定部分可以被补偿，随机部分也可以被控制在更小的范围内，这对加工精度的提高有着至关重要的意义。为了更好的提高数控机床的精度，误差模型的建立也是十分重要的，稳健精确的误差模型也是误差纠正和补偿的第一步。然而，随机波动的误差过大，在生产一些需要控制在特定范围内的零件，就会有面临超差报废的浪费现象。因此，如何更好的表达和分析几何误差的随机部分对于提高加工精度也是十分重要的。误差模型的建立好比绘制了一幅关于几何误差源的地图，这也是进行精度设计还有补偿的最初始最重要的部分。数十年前，科研人员主要解决的问题是机床的误差模型；在近十年，大多数的研究内容主要针对于几何误差模型建立的方法。以更稳健、简洁和精准的描述误差模型，是成为实现误差补偿最基础的要求。这些工作也大多在三轴机床上开展的。国内外专家学者一直在建立数控机床空间误差模型领域进行不懈的探索和研究，开展了多方面的工作。例如三角关系建模法、误差矩阵法、二次关系模型法、机构学建模法、刚体运动学法等。多体系统运动特征分析方法采用齐次列阵表示点的位置和矢量的姿态，在多体系统中建立广义坐标系，将机床抽象为多体系统，将在理想条件下和实际条件下的静态和动态过程中的体间的相对位置和姿态变化以及误差情况作了统一的、完整的描述，使多体系统误差的分析变得简单、迅速、明了和普遍适用，从而为实现计算机快速建模提供基础。然而，在加工过程中，由于机床加工环境具有着大量的不确定性因素，加之机床本身，零件材料等，因此一个零件拥有十分精确的尺寸是不现实的。所以，对于大多数生产过程而言，随机部分误差都被定义成可重复性的即在总体样本中取众数，而随机部分在围绕重复测量的均数以2-3倍的方差来控制随机部分的偏差。
[0004]本项发明以经典的三轴高精度数控机床为例，对于机床的非确定性波动开展预测分析。首先，根据多体系统理论建立起机床的误差模型，在误差模型的基础上，对误差项进行合理的削减至三个方向的“当量误差”。当量误差中也同样存在着非确定性的波动，加工平面时的随机波动在本项发明中，可以根据随机过程理论进行描述和预测。波动的范围也应该被限制在一定范围内；此外对加工误差波动有较大影响的关键误差项会被甄别出来，根据得到的结论，提出一些对于机床零件进行改善的地方。改善后的试验机，可以从测量结果中清晰的看到，精度的提升及波动范围的减小。这对于精密及超精密加工有着至关重要的指导意义。
【发明内容】
[0005]本发明的目的在于提供一种机床空间加工误差的非确定性描述及预测方法，现有研究技术方法中，机床的空间误差被分为两部分:确定性误差和围绕确定性误差的符合某种概率分布的随机波动。虽然随机波动时不能被补偿的，但是尽可能减少波动范围这对于精密和超精密加工过程中，也是至关重要的。
[0006]为实现上述目的，本发明采用的技术方案为机床空间加工误差的非确定性描述及预测方法，首先，根据多体系统理论建立起机床的误差模型，在误差模型的基础上，对误差项进行合理的削减至三个方向的“当量误差”。当量误差中也同样存在着非确定性的波动，在本发明中，加工平面时所拥有随机波动可以根据随机过程理论进行描述和预测。波动的范围也应该被限制在一定范围内；此外对加工误差波动有较大影响的关键误差项会被甄别出来，根据得到的结论，提出一些对于机床零件进行改善的地方。改善后的试验机，可以从测量结果中清晰的看到，精度的提升及波动范围的减小。这对于精密及超精密加工有着至关重要的指导意义。
[0007]如图1所示，本方法的具体实施步骤如下，
[0008]步骤一为三轴机床设置广义坐标系，并建立机床的空间误差模型。
[0009]基于多体系统运动学理论，采用低序体阵列描述抽象机床系统的拓扑结构，在多体系统中建立广义坐标系，用矢量及其列向量表达位置关系，用齐次变换矩阵表示多体系统间的相互关系；
[0010]步骤1.1建立三轴机床的拓扑结构
[0011]分析机床的结构，定义三轴机床的各个组成部件，以及刀具和工件为“典型体”，用“Bj”表示,其中j = 0,1,2...n, j表示各典型体的序号,n-ι表示机床所包含典型体的个数。
[0012]典型体的编号规则如下:
[0013]I)选定床身为典型体“B。”
[0014]2)将三轴机床分为刀具分支和工件分支，共两个分支。首先对刀具分支沿远离床身的方向，按照自然增长数列，对各典型体进行编号。再对工件分支沿远离床身的方向，按照自然增长数列，对各典型体进行编号，如图2，其中m表示刀具分支中典型体的个数，n+1表示机床总共包含的典型体的个数。
[0015]3)任选系统中典型体B」，体Bj的R阶低序体的序号定义为:
[0016]Lr(j) = i (I)
[0017]当Bj体为Bi体的r阶高序体(或Bj体为Bi体的相邻高序体)时，会满足:
[0018]Lr(j) =L(Lr-1U)) (2)
[0019]式中L——低序体算子；
[0020]r、j——自然数
[0021]且补充定义:
[0022]L0(J) = j, Lr(O) = 0 (3)，(4)
[0023]步骤1.2建立三轴机床的特征矩阵。
[0024]该方法所研究的三轴数控机床几何误差项(如图3所示)的几何意义及其表达式如表1所示
[0025]表1:几何误差释义表
[0026]
【权利要求】
1.机床空间加工误差的非确定性描述及预测方法，其特征在于:首先，根据多体系统理论建立起机床的误差模型，在误差模型的基础上，对误差项进行合理的削减至三个方向的“当量误差”；当量误差中也同样存在着非确定性的波动，在本发明中，加工平面时所拥有随机波动可以根据随机过程理论进行描述和预测；波动的范围也应该被限制在一定范围内；此外对加工误差波动有较大影响的关键误差项会被甄别出来，根据得到的结论，提出一些对于机床零件进行改善的地方； 本方法的具体实施步骤如下， 步骤一为三轴机床设置广义坐标系，并建立机床的空间误差模型； 基于多体系统运动学理论，采用低序体阵列描述抽象机床系统的拓扑结构，在多体系统中建立广义坐标系，用矢量及其列向量表达位置关系，用齐次变换矩阵表示多体系统间的相互关系； 步骤1.1建立三轴机床的拓扑结构 分析机床的结构，定义三轴机床的各个组成部件，以及刀具和工件为“典型体”，用“B/’表示，其中j = O, I, 2...n, j表示各典型体的序号，η-1表示机床所包含典型体的个数； 典型体的编号规则如下: 1)选定床身为典型体“BJ 2)将三轴机床分为刀具分支和工件分支，共两个分支；首先对刀具分支沿远离床身的方向，按照自然增长数列，对各典型体进行编号；再对工件分支沿远离床身的方向，按照自然增长数列，对各典型体进行编号，其中m表示刀具分支中典型体的个数，n+1表示机床总共包含的典型体的个数； 3)任选系统中典型体体B」的R阶低序体的序号定义为: Lr(j) = i (I) 当B」体为Bi体的r阶高序体(或B」体为Bi体的相邻高序体)时，会满足: Lr(j) =L(LrtU)) (2) 式中L——低序体算子； r> j——自然数 且补充定义: L0(J) = j，Lr(O) =0 (3)，(4) 步骤1.2建立三轴机床的特征矩阵； 该方法所研究的三轴数控机床几何误差项的几何意义及其表达式如表I所示 表1:几何误差释义表
在床身Btl和所有部件Bj上均建立起与其固定联接的右手直角笛卡尔三维坐标系O0-X0Y0Z0和这些坐标系的集合称为广义坐标系，各体坐标系称为子坐标系，每个坐标系的三个正交基按右手定则分别取名为X，Y, Z轴；各个子坐标系的相对应的坐标轴分别对应平行；坐标轴的正方向与其所对应的运动轴的正方向相同； 将各体之间的运动和静止情况，看作坐标系之间的运动和静止情况；根据两相邻典型体之间的静止和运动情况，在理想运动特征矩阵和误差特征矩阵表中选择相应的运动特征矩阵，如表2 ； 表2:理想运动特征矩阵和运动误差特征矩阵表
其中:TijS表示典型体Bj相对于典型体Bi运动的理想运动特征矩阵； Δ Tijs表示典型体B」相对于典型体Bi运动的运动误差特征矩阵； Xs表示沿X轴平移的距离； Is表示沿Y轴平移的距离； Zs表示沿Z轴平移的距离； 其余参数均已在表1中列出； 若相邻的典型体Bi与典型体B」之间不存在相对运动，则理想运动特征矩阵Tus = I4x4,运动误差特征矩阵ATus = I4x4, I4x4表示4X4的单位矩阵；由于本发明是有关于机床空间加工误差的非确定性描述及预测方法，故使用过程中忽略除几何误差之外的所有误差因素，因此典型体间的体间静止特征矩阵均为Tup = I4x4 ； 根据相邻典型体在静止状态下的实际位置关系，确定典型体间的体间静止误差特征矩阵 ATijp 步骤1.3建立机床的空间误差模型 刀具成型点实际运动位置与理想运动位置的偏差即为机床的空间误差； 设刀具加工点在刀具坐标系中的坐标为:
Pt = [xt, yt, zt, 0]τ (5) 其中Xt表示刀具加工点在刀具坐标系中X轴方向的坐标值； yt表示刀具加工点在刀具坐标系中Y轴方向的坐标值； Zt表示刀具加工点在刀具坐标系中Z轴方向的坐标值； 下标t表示刀具 机床在理想状态时成型点的运动位置:
表示典型体与典型体Bi之间的体间静止特征矩阵； Tijs表示典型体B」与典型体Bi之间的理想运动特征矩阵； Pt表示刀具加工点在刀具坐标系中的坐标； Pwideal表示理想条件下成型点在工件坐标系中的坐标， m+1表示刀具分支中典型体的个数； n+1表示三轴机床所包含的典型体的总个数； 机床在实际状态时成型点的运动位置:
Tijp表示典型体B」与典型体Bi之间的体间静止特征矩阵； Δ Tijs表示典型体B」与典型体Bi之间的体间静止误差特征矩阵； Tijs表示典型体B」与典型体Bi之间的理想运动特征矩阵； Δ Tijs表示典型体B」与典型体Bi之间的运动误差特征矩阵； Pt表示刀具加工点在刀具坐标系中的坐标； 则机床的空间误差模型表示为:
E 一 Pwideal-Pw ⑶ 步骤1.4误差项合理削减及当量误差方程的建立 本发明此步将以空间误差模型为基础，进一步将机床的所有误差项进行合理削减；机床的误差均值模型可以表示为:
F = F (E，G, Pff, U, Uff, Ut, Gv) (9)其中: F = Lf1, f2,..., fr]T:r个独立方程组成的向量； E = [Ex, Ey, Ez, 0]τ:机床的空间误差向量； G = [g1； g2,......，gn]T:n个机床各零部件几何误差组成的向量；
Gv = [ Δ y xy, Δ β χζ, Δ a yz, 1]τ:二主轴间姿态形式误差； Pw = [Pwx, Pwy, Pwz, Ut:工件上成形点在工件坐标系中的坐标向量； U = [X，y, ζ, Β]τ:机床各运动轴的位置向量； Uw = [xw, yw, zw, 1]τ:工件位置坐标向量； Ut = [xt, yt, zt, 1]τ:刀具位置坐标向量； 本发明中定义Pw，U，Uw，Ut是不存在误差；因此，可以进一步写为: F = F(E，G，GV) (10) 其中G的表达式可书写为:
若存在空间误差项，可采用的方法利用激光干涉仪、球杆仪和五坐标测量仪工具来得出；其中对于机床测量方法而言，最常用的方法就是激光干涉仪；优点是可以通过一个轴的测量测出该方向上的6个误差项，总类可以分为直线度误差和线性误差，如果定义有一个与该轴运动形式趋势一致的激光干涉仪测量，此时产生的某些线性误差和直线度误差会有一定的相关性，因此，本发明中定义一个相关系数P表示其中的关系； 在激光干涉仪测量X向的六项基本误差时，并与此同时，在Y的方向，再加之一个激光干涉仪，运动趋势与X向运动一致，此时产生的Y项的6项基本误差就会产生一定的重叠项；x轴沿Y项的线性误差Ayx与Y轴的定位误差Ayy从空间来看二者是存在一定的关系，定义P = Cov ( Ayx, Ayy)便为二者的相关系数；一般情况下，设P = Cov (AIj, AJi)为误差与误差之间的相关系数，其中任意两项定位误差的相关系数为零；同理则可定义出其他误差项之间的相关性，矩阵:
当量误差，由于机床几何误差最终体现在定位精度上，本发明中定义一种新的误差含义:即将空间误差量，投影到在各个轴线上的误差分量；
其中: Λ Xx:Χ项上的当量误差； AYy:Υ项上的当量误差； Λ Zz:Ζ项上的当量误差； 最后得到当量误差方程；步骤二:数控机床各几何误差的测量及其测量数据的整理 激光干涉仪被频繁的用于机床误差检测上，本发明通过定点多测的方法在x，Y，z三个方向上进行测量；分别在各轴50-600mm的行程上，以每20mm为一节点，进行测量重复9次并计算均值；只保留误差值:tr = Tr-D (14) D:目标点； Tr:激光干涉仪测量值； tr:误差值； 使用垂直度测量仪测量机床的三项垂直度误差； 定义各项几何误差均符合&~Ν(μ，。2)均符合高斯分布的独立同分布；
μ:为误差均值； σ2:为误差的方差； 步骤三:计算当量误差并利用随机过程对加工轴及面的随机性波动进行描述及预测 步骤3.1计算当量误差并进行线条拟合 本发明中，认为ΛΧχ，AYy,八乙被设为独立同分布的；根据实验数据的均值，可以计算出三向的当量误差；利用B-spline曲线进行数据在位置点的拟合；拟合原理如下:
其中: U:表示当量误差； P:表示阶数； 步骤3.2轴向随机性描述及预测原理 针对其中一项误差的随机过程，可以将其称之为“高斯白噪声序列”，由白噪声过程定义可知，其中任意两点过程I^n2两点的相关函数爲n|,ll2)与其协方差函数&^相同均为σ2δ Oi1, η2)，并且在移动过程中的任意时刻，均为不相关的，并且任意时刻均为Ν(0，σ 2),于是在这个过程中得到任一点的概率密度函数为:
其中:AXx1:为某一方向上的当量误差； nn:为某一方向上的位置点； 步骤3.3在平面上的随机性描述及预测原理 任意两个当量误差((ΛΧχ，AYy), (AYy, Λ Zz)和(ΛΧχ，Δ Zz))都是独立的随机变量，并且均符合Ν(0，σ 2)分布；定义在X-Y平面上加工一平面，根据随机过程理论；可将平面上的任意点的误差点的误差预测为:
{ΧΥ (η) = Δ XxCos ω η+Δ Yysin ω η, n e (-οο，+ οο )} (19) Λ Xx:Χ向当量误差； AYy:Υ向当量误差； ω:相对加工平面坐标系任一点与远点的矢量角； Exy(η)属于联合高斯过程，从而也可以求出:
Exy (t) = E Δ Xx X cos ω t+E Δ YyX sin ω t = 0 (20) 在机床运作过程中，任意两个过程点I^n2时，可得到他们的相关函数与其协方差函数^?,.4是相等的且均为:
由于、各个点为独立同分布的，因此对于可得到相关系数:
而 ΛΧχ，AYy 是服从 Ν(0, σ 2 ；0, σ 2 ;cos ω (ηι_η2))，其二维密度函数为:
依据本方法，同样可以到得到在Y-z、x-z面的联合概率密度函数； 步骤四:关键误差识别与修改意见 当量误差及其波动作为空间误差项的反应结果，如何将对空间误差项影响较大的误差甄别出来，并减少波动范围就成为此步骤的重点；控制波动范围，最直观的方法是控制影响该项的方差，根据步骤1.4提出的均值误差模型则有:
由于本发明只针对机床的几何误差项则有:
其中偏微分M=αF/aΔ是用来具体识别出具体对于加工影响较大的误差项的，可将其就 某一方向上展开归一化处理:
mni的总量为I ;在某一方向上的mni表示了该项误差对于结果影响的大小； 并可以根据此原理来进行可削减波动范围的关键误差项识别工作； 本发明中，为了验证预测及比较随机性效果，在各轴50-600mm的行程上，以每3mm为一节点随机记录下一组数据。
2.根据权利要求1所述的机床空间加工误差的非确定性描述及预测方法，其特征在于:本发明以三轴精密立式加工中心为例，对上述的机床空间加工误差的非确定性描述及预测方法进行验证； 步骤一:为三轴机床设置广义坐标系，并建立机床的空间误差模型； 基于多体系统运动学理论，采用低序体阵列描述抽象机床系统的拓扑结构，在多体系统中建立广义坐标系，用矢量及其列向量表达位置关系，用齐次变换矩阵表示多体系统间的相互关系； 步骤1.1建立三轴机床的拓扑结构 该机床包括X轴、刀具、工件、Y轴、Z轴、床身； 该三轴数控机床的成型系统由X轴平动单元、Y轴平动单元、Z轴平动单元组成；在数控机床成型运动中，本发明考虑机床的几何误差；本机床共有21项几何误差，包括X，Y，Z轴的六项几何误差(ΔχχΔYχΔζχΔ αχΔ βχΔ Y χ Δ xy Δ yy Δ zy Δ α yA YyAxzAyzAzzAαζΔ βζΔ Yζ)和三项垂直度误差(Λ ΥχYΛ βχζΛ αYζ); 根据多体理论的基本原理将该机床抽象对多体系统，该机床主要由6个典型体组成，定义三轴机床的各个组成部件，以及刀具和工件为“典型体”，用“Β/’表示，其中j =O, I, 2，3，4，5，j表不各典型体的序号，n+1表不机床所包含典型体的个数； 根据编号规则选定床身为典型体“BJ，将三轴机床分为刀具分支和工件分支，共两个分支；首先对刀具分支沿远离床身的方向，按照自然增长数列，对各典型体进行编号；再对工件分支沿远离床身的方向，按照自然增长数列，对各典型体进行编号；步骤1.2建立三轴机床的特征矩阵； 在床身Btl和所有部件Bj上均建立起与其固定联接的右手直角笛卡尔三维坐标系O0-X0Y0Z0和这些坐标系的集合称为广义坐标系，各体坐标系称为子坐标系，每个坐标系的三个正交基按右手定则分别取名为X，Y, Z轴；各个子坐标系的相对应的坐标轴分别对应平行；坐标轴的正方向与其所对应的运动轴的正方向相同； 将各体之间的运动和静止情况，看作坐标系之间的运动和静止情况；根据两相邻典型体之间的静止和运动情况，在理想运动特征矩阵和运动误差特征矩阵表中选择相应的运动特征矩阵；选择结果如表4 ； 表4:该三轴机床的运动特征矩阵和运动误差特征矩阵表
由于B5相对于B。无相对运动,则T5tis = I4x4AT5tis = I4X4 ; B4相对于B3无相对运动,则T34s = I4x4 Δ T34s = I4x4 ； 由于本发明是一种关于机床空间加工误差的非确定性描述及预测方法，在使用过程中忽略除几何误差之外的所有误差因素；根据相邻典型体在静止状态下的位置关系，确定典型体间静止特征矩阵和静止误差特征矩阵；结果如表5 ; 表5:该三轴机床的静止特征矩阵和静止误差特征矩阵表
步骤1.3建立机床的空间误差模型刀具成型点实际运动位置与理想运动位置的偏差即为机床的空间误差；设刀具加工点在刀具坐标系中的坐标为:Pt = [xt, yt, zt, 0]τ (27)其中Xt表示刀具加工点在刀具坐标系中X轴方向的坐标值；yt表示刀具加工点在刀具坐标系中Y轴方向的坐标值；Zt表示刀具加工点在刀具坐标系中Z轴方向的坐标值；下标t表示刀具机床在理想状态时成型点的运动位置:
式中Tup表示典型体B」与典型体Bi之间的体间静止特征矩阵；Tijs表示典型体B」与典型体Bi之间的理想运动特征矩阵；Pt表示刀具加工点在刀具坐标系中的坐标；Pwideal表示理想条件下成型点在工件坐标系中的坐标，机床在实际状态时成型点的运动位置:Pw= [T05]-1 [T01X T12 X T23 X T34] Pt (29)
其中 Tij = Tijp.ATijp.Tijs.ATijs Tijp表示典型体B」与典型体Bi之间的体间静止特征矩阵； Δ Tijs表示典型体B」与典型体Bi之间的体间静止误差特征矩阵； Tijs表示典型体B」与典型体Bi之间的理想运动特征矩阵； Δ Tijs表示典型体B」与典型体Bi之间的运动误差特征矩阵； Pt表示刀具加工点在刀具坐标系中的坐标； 则机床的空间误差模型表示为:
E 一 Pwideal-Pw (30) 步骤1.4误差项合理削减及当量误差方程的建立 本发明此步将以空间误差模型为基础，进一步将机床的所有误差项进行削减；机床的误差均值模型可以表示为:
F = F(E，G，Pw，U，Uw，Ut，Gv) (31) 其中: F = Lf1, f2,..., fr]T:r个独立方程组成的向量； E = [Ex, Ey, Ez, 0]τ:机床的空间误差向量； G = [g1； g2,......，gn]T:n个机床各零部件几何误差组成的向量；
Gv = [ Δ y xy, Δ β χζ, Δ a yz, 1]τ:二主轴间姿态形式误差； Pw = [Pwx, Pwy, Pwz, Ut:工件上成形点在工件坐标系中的坐标向量； U = [X，y, ζ, Β]τ:机床各运动轴的位置向量； Uw = [xw, yw, zw, 1]τ:工件位置坐标向量； Ut = [xt, yt, zt, 1]τ:刀具位置坐标向量； 由于在实际加工过程中，装夹误差及刀具装夹误差是必然要存在误差项的，故本发明中定义Pw，U是没有误差的；因此，可以进一步写为:
F = F(E，G，Gv，Uw，Ut) (32) 其中G的表达式可书写为:
空间误差项，利用激光干涉仪、球杆仪和五坐标测量仪来得出；其中对于机床测量方法而言，最常用的方法就是激光干涉仪；优点是可以通过一个轴的测量测出该方向上的6个误差项，总类可以分为直线度误差和线性误差，如果有一个与该轴运动形式趋势一致的激光干涉仪测量，此时产生的某些线性误差和直线度误差会有一定的相关性，因此，本发明中定义一个相关系数P表示其中的关系； 在激光干涉仪测量X向的六项基本误差时，并与此同时，在Y的方向，再加之一个激光干涉仪，运动趋势与X向运动一致，此时产生的Y项的6项基本误差就会产生一定的重叠项；X轴沿Y项的线性误差Ayx与Y轴的定位误差Ayy从空间来看二者是存在一定的关系，定义P = Cov ( Ayx, Ayy)便为二者的相关系数；一般情况下，设P = Cov (AIj, AJi)为误差与误差之间的相关系数，其中任意两项定位误差的相关系数为零；同理则可定义出其他误差项之间的相关性，矩阵:
当量误差，由于机床几何误差最终体现在定位精度上，本发明中定义一种新的误差含义:即将空间误差量，投影到在各个轴线上的误差分量；
其中: Λ Xx:Χ项上的当量误差； AYy:Υ项上的当量误差； ΛΖζ:Ζ项上的当量误差； 最后得到当量误差方程:
ΔΧχ = Δ Xz- Δ Xx- Δ Xy- Δ Xwd+Z δ β χ-ζ Δ β wd+y Δ y wd_z Δ β y (36)
Δ Yy = ζ [ ( Δ α χ+Δ α y) - ( Δ yx+ Δ yy) ] -χ ( Δ Y wd+Δ Y y+Δ Y xy) - Δ ywd+z Aawd (37)
Δ Zz = χ ( Δ β wd+ Δ β y) + Δ ζζ+ Δ zt+ Δ y+ Δ yt~z Δ a z+y Δ a wd_ Δ Zwd (38) 步骤二:数控机床各几何误差的测量及其测量数据的整理 激光干涉仪被频繁的用于机床误差检测上，本发明中，通过定点多测的方法在X，Y，Z三个方向上进行测量；分别在各轴50-600mm的行程上，以每20mm为一节点，进行测量重复9次并计算均值；只保留误差值:tr = Tr-D (39) D:目标点； Tr:激光干涉仪测量值； tr:误差值； 使用垂直度测量仪测量机床的三项垂直度误差；定义各项几何误差均符合k~Ν(μ，。2)均符合高斯分布的独立同分布；
μ:为误差均值； σ2:为误差的方差； 本发明中，为了验证预测及比较随机性效果，在各轴50-600mm的行程上，以每3mm为一节点随机记录下一组数据；表6~9为50-600mm的行程上，以每20mm为一节点，测量9次并取均值；由于篇幅限制，只列举数据的一部分表6 X轴几何误差测量值均值(mm)
表7 Y轴几何误差测量值(mm)
表8 Z轴几何误差测量值(mm)
表9单元间误差测量值(mm)
步骤三:计算当量误差并利用随机过程对加工轴及面的随机性波动进行描述及预测 步骤3.1计算当量误差并进行线条拟合 本发明中，定义ΛΧχ，AYy,八4被设为独立同分布的；根据实验数据的均值，可以计算出三向的当量误差；利用B-spline曲线进行数据在位置点的拟合；拟合原理如下:
其中: U:表示当量误差； P:表示阶数； 步骤3.2轴向随机性描述及预测原理 针对其中一项误差的随机过程，可以将其称之为“高斯白噪声序列”，由白噪声过程定义可知，其中任意两点过程I^n2两点的相关函数萬,与其协方差函数<^~相同均为σ2δ Oi1, η2)，并且在移动过程中的任意时刻，均为不相关的，并且任意时刻均为Ν(0，σ 2)于是在这个过程中得到任一点的概率密度函数为:
其中:AXx1:为某一方向上的当量误差； ηη:为某一方向上的位置点； 本发明就ΛΧχ，AYy, ΛΖζ三向当量误差添加高斯白噪声序列，并以此来描述和预测机床的几何误差非确定性波动，其波动范围在±3σ之间； 步骤3.3在平面上的随机性描述及预测原理任意两个当量误差((ΛΧχ，AYy), (AYy, Λ Zz)和(ΛΧχ，Δ Zz))都是独立的随机变量，并且均符合Ν(0，σ 2)分布；定义在X-Y平面上加工一平面，根据随机过程理论；可将平面上的任意点的误差点的误差预测为:
{ΧΥ (η) = Δ XxCos ω η+Δ Yysin ω η, n e (-οο，+ οο )} (44) Λ Xx:Χ向当量误差； AYy:Υ向当量误差； ω:相对加工平面坐标系任一点与远点的矢量角； Exy(η)属于联合高斯过程，从而也可以求出:
Exy (t) = E Δ Xx X cos ω t+E Δ YyX sin ω t = 0 (45) 在机床运作过程中，任意两个过程点I^n2时，可得到他们的相关函数夂~_与其协方差函数^一是相等的且均为:
由于、各个点为独立同分布的，因此对于可得到相关系数:
而 ΛΧχ，AYy 是服从 Ν(0, σ 2 ；0, σ 2 ;cosco Cn1-1i2))，其二维密度函数为:
依据本专利的方法，同样可以到得到在Y-z、X-Z面的联合概率密度函数，其波动范围也应在±3 O之间； 步骤四:关键误差识别与修改意见 在本项发明的前面步骤中，已经提及过当量误差及波动预测的求解方法；当量误差及其波动作为空间误差项的反应结果，如何将对空间误差项影响较大的误差甄别出来，并减少波动范围就成为此步骤的重点；控制波动范围，最直观的方法是控制影响该项的方差，根据步骤1.4提出的均值误差模型则有:
由于本项发明只针对机床的几何误差项则有:
其中偏微分
是用来具体识别出具体对于加工影响较大的误差项的，可将其就
mni的总量为I ;在某一方向上的mni表示了该项误差对于结果影响的大小；图15~17分别表示了，在各个方向上对误差结果即非确定性波动范围影响较大的误差项；分别是在X 向上，Δ χζ, Δ χχ, Δ xy, Δ β χ, Δ β y ;Y 方向上，Δ yx, Δ yy, Δ α χ, Δ a y, Δ xy ;Ζ 方向上，Δζζ, ΔΥζ, Δ αζ, Λ i3y对加工结果有着较大影响，本发明为了更为直观的看到的其影响，并可以根据此原理来进行可削减波动范围的关键误差项识别工作。
【文档编号】G06F19/00GK104200063SQ201410384172
【公开日】2014年12月10日 申请日期:2014年8月6日 优先权日:2014年8月6日
【发明者】程强, 冯秋男, 蔡力钢, 章子玲, 赵永胜, 刘志峰 申请人:北京工业大学