4. 将记录的数值输入系统之内，系统会自动生成X轴Z向的直线度误差（图2）
# ]9 D2 N" L/ K2 C* U. y: x/ t9 a+ ?4 p) ]. J
图1 X轴Z向直线度检测示意图
& t3 B9 n m3 z2 @# U- D" E0 ^; p u
" |& d5 w+ O8 i+ l5 T
6 r+ u( [4 j) X+ Z! a8 ?图2 检验结果样例
5 T( Z( h0 D# j( u+ b5 p
# N8 k% x, ?3 R- U5 |: E' G
问题的提出：
# e/ O s P: s3 E4 i公司规定的该处直线度误差应为11μm之内，据车间及工艺反映：，自从该机型量产后的几十台检验结果，该处直线度误差一直在超差边缘，机械调整非常困难，公司内部不允许在工厂阶段进行软件优化补偿，因为现在我主管该产品了，恰好之前按照总部的要求将导轨及丝杠都做了优化改进，铸件也都做了更改，超差的两台机器正好是前两台样机，于是质量部门召开了一次会议，要求从机械上想办法对该问题做出一些分析，并提出改进的办法。
" o1 z4 F8 x0 N& B+ o: e! Z. [& E, u解决问题的思路：
& M& e( }6 o$ c) V6 L( U联想到该机床的基本结构，立柱中间为了放置转台开了一个很大的口子，这必然影响到整个机床在此处的刚性，反映到上面提到的x轴z方向直线度便是两侧高而中间低的情况，这在软件上也非常清楚的显示了出来（图3），解决的办法就是在立柱的导轨上加工出一个中凸的曲线来抵消重力对直线度的影响。
% U3 n& a( t$ L1 R, \9 s3 m
) t( c/ m. p6 k% `$ F( n图3 直线度误差曲线图
3 ?5 m1 [4 j/ D/ d- N5 N
问题的一般化：
6 B# V3 H g; u2 n7 a) s% V0 N我的雄心绝不仅仅是随便给出一个结果这么简单，我想把这个问题更加一般化:
& N* ]8 }2 K1 L# `# m
1.如何根据输入的几个点求整条轴的直线度?
/ g0 N9 Q& M2 I: B9 _7 x2.如果该直线度的走向有规律可循，如何找到该规律？换言之，该直线度的走向能否被一条曲线拟合，如果能，那曲线的方程应该是什么？
. Q: i- D: t! K. f# f
3.如果我手上有几十台该直线度的检测数据，我如何来求得最适合的一条拟合曲线？
! N* l* E; p2 u1 r# M# \
4.如何求得3中提出来的拟合曲线的补偿曲线？
: d* h7 }# C e6 r2 ~5 X8 X8 ~- ~5.补偿曲线求好后，如何反推其正确性（通过已有的机床直线度数据反求）
e F5 }' X, L具体的解决思路：
) @% M. i9 f% C2 W为了便于理解，我对机床轴和坐标轴的对应关系做一个描述：令机床轴X轴的中点作为坐标的原点，坐标轴X轴的正向为机床轴X轴行程的正向，这样，这个问题中700mm的行程可以在坐标轴上对应为为X轴的-350，-320，…，0，…，320，350，单位为毫米，坐标轴Y轴和机床Z轴重合，方向一直，单位为微米。
7 {' `2 J* I# k# y1.如何根据输入的几个点求整条轴的直线度?
0 |! K* a, \. \. O6 E
这个问题比较简单，只要在坐标轴上画出每个点对应的X轴Z方向的直线度测量数据，然后通过最小二乘法拟合出一条直线，每个点与该直线的竖直距离（注意是Z方向差值的绝对值，而不是点到直线的垂直距离，这个具体可以参考最小二乘法的定义及求法）的差值即为该点的偏差值，这个行程上点与点之间偏差值的差的最大值就是X轴在Z方向的直线度。下面就以图2的检测数据为例编程求出该机床的X轴Z方向直线度（图4）。
$ ^/ k7 J _/ t7 f9 j, T1 |$ [5 F图4 X轴Z方向的直线度求法模拟
& c3 Z* v% V- F
由图中可以看出，用该方法求出的偏差值及直线度误差与系统给出的完全一致，这里给出的精度甚至更高。（具体的编程代码我就不给出了，熟悉编程的朋友可以自己去编个看看）
2.如果该直线度的走向有规律可循，如何找到该规律？换言之，该直线度的走向能否被一条曲线拟合，如果能，那曲线的方程应该是什么？
2 X# N) l9 C7 t5 u上面既然求出了每个点的实际偏差值，问题就简单了：把实际偏差值输入到坐标轴中，用合适的曲线模拟出偏差的趋势或走向，用的当然还是最小二乘法。至于拟合曲线的选取，我这里选择了二次曲线，当然读者也可以选更高次的曲线模拟。（图5）
: A/ [7 b7 r' G# y0 a& r8 H
: E# R- U0 f6 a% @9 T. s( o1 L* r
图5 用二次曲线拟合直线度偏差点
* \5 J: T0 K8 m4 D U: [) c: B' R/ c. s- [1 e
从图中可见，模拟的效果还是不错的，点的位置基本上符合了二次曲线的规律。
) H' a9 X2 G9 W( a/ F, O" n
8 w8 Z; Z E6 W: T' {3.如果我手上有几十台该直线度的检测数据，我如何来求得最适合的一条拟合曲线？
( C/ @7 l6 f T2 j% |从问题2中我们可以求得每台机床的直线度情况，但是供应商加工铸件的精度在一定范围内是随机的，这就需要我们从大数据中求得最适合的一条曲线，通过该条曲线的加工补偿，可以使得每台机床直线度都不至于偏离太远。这个最适合应该怎么定义呢？我还是使用最小二乘法的原理，在轴的每个特定位置求出刚刚拟合出的几十条曲线的直线度偏差值，然后找到该位置上的一个偏差值使得该值到其他机床在该位置偏差的距离的平方和最小，如此再做最小二乘拟合，所得的二次曲线便是几十台机床直线度误差的平均值了。（图6）
3 [+ R s. P: l9 K4 r# e( U3 f
9 C8 r6 ]' j( z1 R3 |) z+ R% N
图6 多台机床直线度数据拟合
1 ~, L3 m$ G1 g8 p( ~, f/ E
- B O/ o% R2 ]% G: c
本次我一共拟合了10台的数据，得到了上述的表达式和函数图像
3 e' W& i9 |3 y3 ^
4.如何求得3中提出来的拟合曲线的补偿曲线？
! k( ]1 H! v% i2 d |9 a
上面的曲线方程出来以后，这一步也就可以顺利解答了，只要把上面的方程按X轴对称即可得到补偿函数，当然这个函数还不是我们要的，为了方便加工，给出的函数需要假设两头为0，中间高，从而让加工厂不必费力再去计算各点的差值，直接从函数中就可以计算出每个点的补偿值。（图7）当然考虑的实际的精度，中间一项可以删除。从而得到更加简洁的函数表达式。
0 x- Q6 M1 R+ |- }/ Y3 M图7 最终得到的可用于加工的补偿函数及图像
( i$ k0 G0 K; a
, c/ U: N$ C( [
5.补偿曲线求好后，如何反推其正确性（通过已有的机床直线度数据反求）
把上面求出的补偿曲线与每台机床的直线度误差曲线合并后，再使用步骤1的方法可以得到补偿后的结果。（图8）
! L. w0 O7 E0 j0 O! K) z7 k6 v- {3 b2 e7 K6 O/ w
图8 补偿前后直线度误差对比
* T+ Z; S. o1 P/ }/ }# ?! ?2 ?
+ i1 {9 I, o) S5 D: k4 R
由上图可以看出，补偿后的直线度误差没有超过6微米的，除了第5台机床，其他的都要明显好于补偿前的结果。
! [" D0 }% {7 L! o' E! }. {" N n
结论：. K+ |$ t" d3 C c$ e7 p- k/ g
通过机械的补偿可以达到预期的效果，但是值得注意的是，补偿后的误差仍然可以通过系统插值法来进一步补偿。而软件补偿是针对单台机床来说的，从下图可以看到，如果需要，二次补偿后直线度在1微米之内甚至更高精度要求内完全没有问题。
# `0 ~( g, D( M& j. }0 J$ U8 V8 O
% K6 {( @6 }' D: \5 X G5 O" O图9 补偿后的偏差仍然可以通过二次曲线在软件中补偿
( A8 X. Z: A& J5 q& \3 K
* c z9 j0 N' b8 G6 s) d1 t4 a0 U4 K: U
5 [" W( ]) I6 e0 j0 q5 t, Y. n" f
( Z7 t& {! {/ O8 v9 b
发表于 2016-3-5 23:24:41
