提高三点法测量精度的研究

HEATS

高等学校博士学科点专项科研基金资助课题为提高三点法^［1］的测量精度，前人倾注了巨大精力，但其研究大都限于仅从数学上考虑了如何避免三点法的谐波加权函数ω(k)=1+a₂e^-jkφ₂+a₃e^-jkφ₃造成的谐波抑制^［2～4］，而忽视了测量误差对测量精度的影响.文献［5］虽然注意到上述问题，但未能透彻分析测头读数及定位误差与圆度测量误差的关系.作者重点分析测头读数及定位误差对圆度测量精度的影响，研究系统的优化方法及谐波测量误差的传播规律，以进一步提高三点法测量精度.

1　三点法测量原理及精度分析

1.1　三点法测量原理

9 z2 K1 I9 c# P/ o5 [

+ z9 h" S- N0 m; s6 M4 P

图1　三点法原理

三点法测量原理如图1所示.在工作圆周上布置3个测头P_i，它们与y轴的夹角为φ_i(φ₁=0)，工作圆度为r(θ_n)，轴系径向运动误差的x，y向分量为x(θ_n)，y(θ_n)，各测头的输出为

) ]$ _* X3 b! o; z* R

其中　R_i为测头P_i的零起读数； θ_n为采样点，θ_n=2nπ/N，n=0,1,…，N-1.构造线性方程

9 z0 u* ~5 _' \" {3 P$ Q

记　

设s₀(θ_n)的Fourier系数为F_k和G_k，为分离x(θ_n)，y(θ_n)与r(θ_n)，使α₁=β₁=0，由此确定a₂，a₃.经适当变换，得圆度各次谐波的Fourier系数为

/ L' {. K3 o$ `+ O* \# Z

W3 o& Y$ \- s+ w

由此得到工件圆度

/ X' G' Q* x* r6 H+ B o+ L0 q; ?

其中　K_m为谐波最高级次.由方程(1)可求得轴系径向运动误差.
1.2　三点法测量系统精度分析
　　影响测圆精度的主要误差有： ① 三个测头读数误差δ_1n,δ_2n,δ_3n，它影响F_k,G_k； ② 测头角位置误差δ(φ₂)，δ(φ₃)，它影响α_k,β_k.由式(4)得

% A) l: s3 `/ F% u. F

1.2.1　测头读数误差对测量精度的影响
　　此时忽略测头角位置误差的影响，则

　　设3个测头读数误差的方差均为σ²，由随机变量均值与方差的关系^［6］，A_k，B_k的方差为

( R1 r9 v) ]# A# i7 H6 S; X* }! X/ V

　　由测头读数误差引入的圆周第p点的圆度测量误差为

$ c& ?+ \, `- g1 r% Y2 I; v. j

其中　.
则

1.2.2　测头角位置误差对测量精度的影响
　　此时，若忽略读数误差的影响，则

7 |6 J, K6 Y! D3 h" |

　　由得第k次谐波的幅值相对误差为

　　由tanφ_k=A_k/B_k得第k次谐波的相位误差

9 r; h( u# C. ^- l2 I$ z6 U

　　设,q_k是表示测头读数及角位置误差对圆度及其各次谐波测量精度影响大小的重要指标.显然，为降低测头读数误差及测头角位置误差对测量结果的影响，应使q_k尽量小.计算机仿真表明：圆度测量误差主要由测头读数误差引起，角位置误差影响较小，当测头角位置的测量精度高于1′时，测头位置误差所引入的测量误差可以忽略.

. i. s9 r2 m" [" z$ W* j

2　三点法测量系统参数的优化

　　优化的目的是选择一组最佳的测头角位置，使测头读数及角位置误差对圆度各次谐波测量精度的影响为最小.优化过程为一个多目标函数的优化决策问题.为降低测头读数误差及测头角位置误差对测圆精度的影响，对于给定的级次K_m，求φ₂，φ₃，使所有各次谐波的q_k中的最大值为最小.
　　以q_k构成的目标函数特性复杂，有许多局部极小值，常用的单峰函数寻优方法无法得到全局最优值，而采用枚举法^［5，7］搜索耗时很长.为此，作者基于Monte Carlo随机试验思想，并参考单纯形模式寻优方法^［7］，编制了随机模式搜索方法的算法程序.其方法如下：
　　① 在由φ_imin＜φ_i＜φ_imax所决定的平面内产生一个随机试验点φ⁽⁰⁾；
　　② 以初始点φ⁽⁰⁾为中心，以A₁和A₂为半径作一空心圆域，在这空心圆域中按均匀分布产生随机变量搜索5～10次，若不成功，则反向搜索一次，求出新点φ⁽¹⁾，使f(φ⁽¹⁾)＜f(φ⁽⁰⁾)；
　　③ 进行一次模式寻优，令φ=2φ⁽¹⁾-φ⁽⁰⁾，若f(φ)＜f(φ⁽¹⁾),表示模式寻优成功，转向步骤②，以新点为中心，继续寻找更好的点；
　　④ 若f(φ)＞f(φ⁽¹⁾)，表示模式寻优失败，将φ向φ⁽¹⁾收缩，直到f(φ)＜f(φ⁽¹⁾)为止，重复步骤②；
　　⑤ 返回①，重复上述试验2000～5000次，比较各次试验结果，选取最佳结果，以确保得到全局最优点.
　　上述方法的可靠性及算法误差可通过两种途径考察： ① 取寻优结果的相对偏差小于10^-31为截止判据； ② 优于文献［5］中的结果.图2给出了K_m为20和160时的最佳角度及最佳角度下的q_k，由此可看出三点法各次谐波测量误差的传播规律.作者计算出的最佳角度下的q_k值，不仅明显优于按传统方法确定的最佳角度下的q_k(见图2b与图2d)，而且明显优于文献［5］(见图2a与图2c).试验表明，随着所取K_m的增大，对应的各次谐波的q_k变大.使用某K_m下确定的系统测量大于此谐波上限的谐波时，由于误差传递系数很大，较小的读数误差就会使圆度测量结果产生较大失真.因此正确估计待测件谐波上限并以此确定测头的最佳角位置至关重要.

0 {# f( w- \& M/ }

(a)优化方法结果

. f) m/ j* P4 {# V' d# G9 P. s9 m

% }; z; f" h9 S5 H; n( I g* ^

(b)优化方法结果

* w5 p# g5 p/ m0 I8 i, x- T

(c)优化方法结果

4 ]8 x! h; g& l9 B) l3 k0 x

(d)传统方法结果

图2　三点法测头最佳角位置及圆度各次谐波测量误差传递系统

3　测头角位置测量误差及测头灵敏度不一致的校正

3 K% _7 x% i) A7 u5 W3 x7 ?

　　实际测量中，考虑到： ① 测头角位置存在测量误差Δφ₂，Δφ₃； ② 3个测头的灵敏度不一致，设它们的比值为1∶p₂∶p₃.则圆度与轴系误差分离的基本条件应为

2 m( t) H) ?+ o5 Q2 j1 @

: r3 q {6 ?" W. z! X7 I

显然，由α₁=β₁=0得出的a₂，a₃无法保证上式成立.这导致： ① 圆度与轴系误差不能彻底分离，残留的轴系运动误差计入到圆度中； ② 计算得出的圆度各次Fourier系数被虚假地放大，造成较大的测量误差.
　　当满足误差分离的基本条件时，α₁=β₁=0，此时，S₀(θ_n)的一次谐波幅值C²₁=F²₁+G²₁=0.当φ₂，φ₃有误差时，｜C₁｜＞0.据此可以对φ₂，φ₃的直接测量值进行优化以确定其真值.方法是：当φ₂，φ₃变化时，a₂，a₃随之变化，C₁的值也随之变化，如果｜C₁｜变小，则说明φ₂，φ₃趋近于其真值.即：求φ₂，φ₃，φ₂min≤φ₂≤φ₂max,φ3min≤φ₃≤φ3max，使｜C₁(φ₂，φ₃)｜最小.其中，φ_imin，φ_imax由直接测量误差限确定.采用直接搜索的方法，即可求得φ₂，φ₃的真值.
　　实验表明：通过优化，可以使｜C₁｜的数值明显减小，有效地提高了φ₂，φ₃的测量精度.表1给出了5次不同测量列S₀(θ_n)的优化结果，显然，优化结果具有较好的重复性精度.

3 H& }2 \, b( i

表1　测头角间距φ₂，φ₃的优化结果

/ J) e2 h+ |! f* q& C5 e; y4 c- G# H$ B9 @- ?2 {4 o8 Q( L n& h' F3 X% o2 ~9 T! b1 V0 G# d; n+ e/ k! O0 C2 `9 d/ o9 V) G; D4 B l! i2 I I: y% _ |! j3 v+ ?+ r9 e+ }! W/ ~* W5 Z6 P+ j3 A$ m0 @4 X1 N% C2 B. R7 F1 [& A8 _, v( Y9 ` W* N2 p2 l+ _ L% P5 f3 i8 s/ `& A! R& a9 P4 r; S* y( v$ `% T# u* ]& H9 q! ^+ K E6 V. | u9 a: ] \/ ~+ `+ c7 d- |4 P0 S, g) ~! u3 w) C" H' B+ H- s' Z8 Q4 C2 t3 {- W; H' b& o, e- n" A3 @$ h3 a& a, a% V" M* i5 i5 y% D5 H7 |: a, u( U A9 \" @' P( `9 m. z, U* Q0 J! ?: v# }9 Q

初始值	测量值	优化值					平均值	方差
		1	2	3	4	5
φ2/(°)	112.45	112.49	112.45	112.47	112.49	112.47	112.48	0.014
φ3/(°)	230.23	229.91	229.89	229.90	229.86	229.88	229.89	0.018
｜C1｜/10-6	3521	891	957	1264	1018	984

' I; ]9 h3 k8 @! |/ Q, q: O* [

& v" {. Q: b% o+ s4 _9 K9 s" ~& h) j5 Q7 s1 m m

2 P/ ~# h. `$ X* @6 Y/ H2 z: K1 p 4　实验结果及分析 　　分别采用某文献推荐的φ2=84.24°，φ3=168.78°及本文中取Km=160时优化出的角度安置测头，对一个直径为110mm的标准件进行测量，并与Talyrond圆度仪测得结果相比较，见表2. 表2　测量结果比较

: g) X" Z7 ?1 I - ~5 y' W5 |8 j5 k% ] a* |( i9 I7 P! [' k- J+ @; F6 ]4 }$ u) `5 c" Q2 w$ p( a6 I$ r2 w/ E( I: Y$ ?8 b+ C" @. A6 R/ J! h8 _! {) C6 E% d' ]! C- @: a. J- R9 b- c# A1 U6 V, b4 I. @. ^! }& t/ O) n6 ?! J6 j# _5 Q; M% t9 m5 |4 M- J0 r/ v" r, T: `9 i% e; b( J/ o, k) I6 o: i4 `) t% h) C0 M; h, V9 g2 ~ E7 L0 [0 o3 i2 `5 @4 F+ m* Q2 I3 H% y* A" |7 p/ \+ P# y- E# W/ I3 N4 n& c' W8 k; M8 d+ L9 H0 x& G3 ~6 g; C$ n7 d; ?- J2 X8 j+ o6 X- b. K7 w# J+ k9 ?& d$ h8 \( @! a- ?: d4 ~7 k* u4 [, l7 A2 ^( N; F

方　法	圆　　　度				谐波误差		圆周各点测量误差
	均值/μm	改善程度/%	σ/μm	改善程度/%	最大值/μm	改善程度/%	最大值/μm	改善程度/%
圆度仪	1.253		0.050
传统三点法	1.292		0.035		0.072		0.250
优化三点法	1.282	26	0.026	26	0.024	67	0.160	36

5 C2 N) F% K9 k W/ f$ H7 e

1 C3 G S$ j6 c# T1 l3 g/ Q/ j, d4 g* S5 \- _- _( Z% ~" t" l$ O; R- A0 c

结果表明： 　　① 传统三点法因某些级次的误差传递系数qk较大，读数误差被大大放大，从而使测量结果中有明显的高次谐波误差(见图2d与图3a). 　　② 传统三点法测得圆度以及各次谐波值的离散程度较大，单圈测量值与圆度仪测得值相差较大，必须取多圈平均才能求得较精确的结果. 　　③ 优化三点法圆度以及各次谐波值的测量不确定性明显降低，尤其高次谐波测量误差明显降低(见图3b).

0 Y( p# {7 ]/ x K, s* G

(a)传统三点法结果

. `* X0 h: g% q8 V+ s3 u) L

(b)优化三点法测量结果

; L7 _, `4 n% W8 h. S% ^ X

图3　圆度和各次谐波幅值测量误差

! Q& {4 z; H. a1 g& M4 `. p* U