提高三点法测量精度的研究
<FONT face=宋体>高等学校博士学科点专项科研基金资助课题为提高三点法<SUP>[1]</SUP>的测量精度,前人倾注了巨大精力,但其研究大都限于仅从数学上考虑了如何避免三点法的谐波加权函数ω(k)=1+a<SUB>2</SUB>e<SUP>-j</SUP>kφ<SUB>2</SUB>+a<SUB>3</SUB>e<SUP>-j</SUP>kφ<SUB>3</SUB>造成的谐波抑制<SUP>[2~4]</SUP>,而忽视了测量误差对测量精度的影响.文献[5]虽然注意到上述问题,但未能透彻分析测头读数及定位误差与圆度测量误差的关系.作者重点分析测头读数及定位误差对圆度测量精度的影响,研究系统的优化方法及谐波测量误差的传播规律,以进一步提高三点法测量精度.</FONT><P align=left><STRONG><FONT face=宋体>1 三点法测量原理及精度分析</FONT></STRONG></P>
<P align=left><STRONG><FONT face=宋体>1.1 三点法测量原理</FONT></STRONG></P>
<P align=center><IMG src="http://www.chmcw.com/upload/news/RCL/13220_2lghvn2008324135258.gif"></P>
<P align=center><FONT face=宋体><STRONG>图1 三点法原理</STRONG></FONT></P>
<P align=left><FONT face=宋体>三点法测量原理如图1所示.在工作圆周上布置3个测头P<SUB>i</SUB>,它们与y轴的夹角为φ<SUB>i</SUB>(φ<SUB>1</SUB>=0),工作圆度为r(θ<SUB>n</SUB>),轴系径向运动误差的x,y向分量为x(θ<SUB>n</SUB>),y(θ<SUB>n</SUB>),各测头的输出为</FONT></P>
<P><IMG src="http://www.chmcw.com/upload/news/RCL/13220_9znyrj2008324135315.gif"></P>
<P align=left><FONT face=宋体>其中 R<SUB>i</SUB>为测头P<SUB>i</SUB>的零起读数; θ<SUB>n</SUB>为采样点,θ<SUB>n</SUB>=2nπ/N,n=0,1,…,N-1.构造线性方程</FONT></P>
<P><IMG src="http://www.chmcw.com/upload/news/RCL/13220_5ify262008324135330.gif"></P>
<P><FONT face=宋体>记 </FONT><IMG src="http://www.chmcw.com/upload/news/RCL/13220_paymmp2008324135344.gif"></P>
<P align=left><FONT face=宋体>设s<SUB>0</SUB>(θ<SUB>n</SUB>)的Fourier系数为F<SUB>k</SUB>和G<SUB>k</SUB>,为分离x(θ<SUB>n</SUB>),y(θ<SUB>n</SUB>)与r(θ<SUB>n</SUB>),使α<SUB>1</SUB>=β<SUB>1</SUB>=0,由此确定a<SUB>2</SUB>,a<SUB>3</SUB>.经适当变换,得圆度各次谐波的Fourier系数为</FONT></P>
<P><IMG src="http://www.chmcw.com/upload/news/RCL/13220_ddt1g82008324135359.gif"></P>
<P align=left><FONT face=宋体>由此得到工件圆度</FONT></P>
<P><IMG src="http://www.chmcw.com/upload/news/RCL/13220_mykjr52008324135412.gif"></P>
<P align=left><FONT face=宋体>其中 K<SUB>m</SUB>为谐波最高级次.由方程(1)可求得轴系径向运动误差.<BR></FONT><FONT face=宋体><STRONG><STRONG>1.2 三点法测量系统精度分析</STRONG></STRONG></FONT><FONT face=宋体><BR> 影响测圆精度的主要误差有: ① 三个测头读数误差δ<SUB>1n</SUB>,δ<SUB>2n</SUB>,δ<SUB>3n</SUB>,它影响F<SUB>k</SUB>,G<SUB>k</SUB>; ② 测头角位置误差δ(φ<SUB>2</SUB>),δ(φ<SUB>3</SUB>),它影响α<SUB>k</SUB>,β<SUB>k</SUB>.由式(4)得</FONT></P>
<P><IMG src="http://www.chmcw.com/upload/news/RCL/13220_huw1t12008324135426.gif"></P>
<P align=left><FONT face=宋体><STRONG>1.2.1 </STRONG>测头读数误差对测量精度的影响<BR> 此时忽略测头角位置误差的影响,则</FONT></P>
<P><IMG src="http://www.chmcw.com/upload/news/RCL/13220_x0cysg2008324135438.gif"></P>
<P align=left><FONT face=宋体> 设3个测头读数误差的方差均为σ<SUP>2</SUP>,由随机变量均值与方差的关系<SUP>[6]</SUP>,A<SUB>k</SUB>,B<SUB>k</SUB>的方差为</FONT></P>
<P><IMG src="http://www.chmcw.com/upload/news/RCL/13220_tgacay2008324135453.gif"></P>
<P align=left><FONT face=宋体> 由测头读数误差引入的圆周第p点的圆度测量误差为</FONT></P>
<P><IMG src="http://www.chmcw.com/upload/news/RCL/13220_r6ndsn2008324135526.gif"></P>
<P align=left><FONT face=宋体>其中 <IMG src="http://www.chmcw.com/upload/news/RCL/13220_hjjszo2008324135536.gif">.<BR>则</FONT></P>
<P><IMG src="http://www.chmcw.com/upload/news/RCL/13220_kw2sim2008324135554.gif"></P>
<P align=left><FONT face=宋体><STRONG>1.2.2 </STRONG>测头角位置误差对测量精度的影响<BR> 此时,若忽略读数误差的影响,则</FONT></P>
<P><IMG src="http://www.chmcw.com/upload/news/RCL/13220_k4lswf2008324135559.gif"></P>
<P align=left><FONT face=宋体> 由<IMG src="http://www.chmcw.com/upload/news/RCL/13220_dgmoxu2008324135620.gif">得第k次谐波的幅值相对误差为</FONT></P>
<P><IMG src="http://www.chmcw.com/upload/news/RCL/13220_gsr7yh2008324135613.gif"> </P>
<P align=left><FONT face=宋体> 由tanφ<SUB>k</SUB>=A<SUB>k</SUB>/B<SUB>k</SUB>得第k次谐波的相位误差</FONT></P>
<P><IMG src="http://www.chmcw.com/upload/news/RCL/13220_reu1nx2008324135637.gif"></P>
<P align=left><FONT face=宋体> 设<IMG src="http://www.chmcw.com/upload/news/RCL/13220_htms2j2008324135656.gif">,q<SUB>k</SUB>是表示测头读数及角位置误差对圆度及其各次谐波测量精度影响大小的重要指标.显然,为降低测头读数误差及测头角位置误差对测量结果的影响,应使q<SUB>k</SUB>尽量小.计算机仿真表明:圆度测量误差主要由测头读数误差引起,角位置误差影响较小,当测头角位置的测量精度高于1′时,测头位置误差所引入的测量误差可以忽略.</FONT></P>
<P align=left><FONT face=宋体><STRONG>2 三点法测量系统参数的优化</STRONG></FONT></P>
<P align=left><FONT face=宋体> 优化的目的是选择一组最佳的测头角位置,使测头读数及角位置误差对圆度各次谐波测量精度的影响为最小.优化过程为一个多目标函数的优化决策问题.为降低测头读数误差及测头角位置误差对测圆精度的影响,对于给定的级次K<SUB>m</SUB>,求φ<SUB>2</SUB>,φ<SUB>3</SUB>,使所有各次谐波的q<SUB>k</SUB>中的最大值为最小.<BR> 以q<SUB>k</SUB>构成的目标函数特性复杂,有许多局部极小值,常用的单峰函数寻优方法无法得到全局最优值,而采用枚举法<SUP>[5,7]</SUP>搜索耗时很长.为此,作者基于Monte Carlo随机试验思想,并参考单纯形模式寻优方法<SUP>[7]</SUP>,编制了随机模式搜索方法的算法程序.其方法如下:<BR> ① 在由φ<SUB>i</SUB>min<φ<SUB>i</SUB><φ<SUB>i</SUB>max所决定的平面内产生一个随机试验点φ<SUP>(0)</SUP>;<BR> ② 以初始点φ<SUP>(0)</SUP>为中心,以A<SUB>1</SUB>和A<SUB>2</SUB>为半径作一空心圆域,在这空心圆域中按均匀分布产生随机变量搜索5~10次,若不成功,则反向搜索一次,求出新点φ<SUP>(1)</SUP>,使f(φ<SUP>(1)</SUP>)<f(φ<SUP>(0)</SUP>);<BR> ③ 进行一次模式寻优,令φ=2φ<SUP>(1)</SUP>-φ<SUP>(0)</SUP>,若f(φ)<f(φ<SUP>(1)</SUP>),表示模式寻优成功,转向步骤②,以新点为中心,继续寻找更好的点;<BR> ④ 若f(φ)>f(φ<SUP>(1)</SUP>),表示模式寻优失败,将φ向φ<SUP>(1)</SUP>收缩,直到f(φ)<f(φ<SUP>(1)</SUP>)为止,重复步骤②;<BR> ⑤ 返回①,重复上述试验2000~5000次,比较各次试验结果,选取最佳结果,以确保得到全局最优点.<BR> 上述方法的可靠性及算法误差可通过两种途径考察: ① 取寻优结果的相对偏差小于10<SUP>-31</SUP>为截止判据; ② 优于文献[5]中的结果.图2给出了K<SUB>m</SUB>为20和160时的最佳角度及最佳角度下的q<SUB>k</SUB>,由此可看出三点法各次谐波测量误差的传播规律.作者计算出的最佳角度下的q<SUB>k</SUB>值,不仅明显优于按传统方法确定的最佳角度下的q<SUB>k</SUB>(见图2b与图2d),而且明显优于文献[5](见图2a与图2c).试验表明,随着所取K<SUB>m</SUB>的增大,对应的各次谐波的q<SUB>k</SUB>变大.使用某K<SUB>m</SUB>下确定的系统测量大于此谐波上限的谐波时,由于误差传递系数很大,较小的读数误差就会使圆度测量结果产生较大失真.因此正确估计待测件谐波上限并以此确定测头的最佳角位置至关重要.</FONT></P>
<P align=center><IMG src="http://www.chmcw.com/upload/news/RCL/13220_dx14t32008324135715.gif"></P>
<P align=center><FONT face=宋体>(a)优化方法结果</FONT></P>
<P align=center><IMG src="http://www.chmcw.com/upload/news/RCL/13220_85xyfe2008324135721.gif"></P>
<P align=center><FONT face=宋体>(b)优化方法结果</FONT></P>
<P align=center><IMG src="http://www.chmcw.com/upload/news/RCL/13220_i5o7tb2008324135835.gif"></P>
<P align=center><FONT face=宋体>(c)优化方法结果</FONT></P>
<P align=center><IMG src="http://www.chmcw.com/upload/news/RCL/13220_dmmue42008324135849.gif"></P>
<P align=center><FONT face=宋体>(d)传统方法结果</FONT></P>
<P align=center><FONT face=宋体><STRONG>图2 三点法测头最佳角位置及圆度各次谐波测量误差传递系统</STRONG></FONT></P>
<P align=left><STRONG><FONT face=宋体>3 测头角位置测量误差及测头灵敏度不一致的校正</FONT></STRONG></P>
<P align=left><STRONG><FONT face=宋体> </FONT></STRONG><FONT face=宋体>实际测量中,考虑到: ① 测头角位置存在测量误差Δφ<SUB>2</SUB>,Δφ<SUB>3</SUB>; ② 3个测头的灵敏度不一致,设它们的比值为1∶p<SUB>2</SUB>∶p<SUB>3</SUB>.则圆度与轴系误差分离的基本条件应为</FONT><STRONG></P>
<P></STRONG><IMG src="http://www.chmcw.com/upload/news/RCL/13220_t5akoc2008324135933.gif"></P>
<P align=left><FONT face=宋体>显然,由α<SUB>1</SUB>=β<SUB>1</SUB>=0得出的a<SUB>2</SUB>,a<SUB>3</SUB>无法保证上式成立.这导致: ① 圆度与轴系误差不能彻底分离,残留的轴系运动误差计入到圆度中; ② 计算得出的圆度各次Fourier系数被虚假地放大,造成较大的测量误差.<BR> 当满足误差分离的基本条件时,α<SUB>1</SUB>=β<SUB>1</SUB>=0,此时,S<SUB>0</SUB>(θ<SUB>n</SUB>)的一次谐波幅值C<SUP>2</SUP><SUB>1</SUB>=F<SUP>2</SUP><SUB>1</SUB>+G<SUP>2</SUP><SUB>1</SUB>=0.当φ<SUB>2</SUB>,φ<SUB>3</SUB>有误差时,|C<SUB>1</SUB>|>0.据此可以对φ<SUB>2</SUB>,φ<SUB>3</SUB>的直接测量值进行优化以确定其真值.方法是:当φ<SUB>2</SUB>,φ<SUB>3</SUB>变化时,a<SUB>2</SUB>,a<SUB>3</SUB>随之变化,C<SUB>1</SUB>的值也随之变化,如果|C<SUB>1</SUB>|变小,则说明φ<SUB>2</SUB>,φ<SUB>3</SUB>趋近于其真值.即:求φ<SUB>2</SUB>,φ<SUB>3</SUB>,φ<SUB>2</SUB>min≤φ<SUB>2</SUB>≤φ<SUB>2</SUB>max,φ3min≤φ<SUB>3</SUB>≤φ3max,使|C<SUB>1</SUB>(φ<SUB>2</SUB>,φ<SUB>3</SUB>)|最小.其中,φ<SUB>i</SUB>min,φ<SUB>i</SUB>max由直接测量误差限确定.采用直接搜索的方法,即可求得φ<SUB>2</SUB>,φ<SUB>3</SUB>的真值.<BR> 实验表明:通过优化,可以使|C<SUB>1</SUB>|的数值明显减小,有效地提高了φ<SUB>2</SUB>,φ<SUB>3</SUB>的测量精度.表1给出了5次不同测量列S<SUB>0</SUB>(θ<SUB>n</SUB>)的优化结果,显然,优化结果具有较好的重复性精度.</FONT></P>
<P align=center><FONT face=宋体><STRONG>表1 测头角间距φ<SUB>2</SUB>,φ<SUB>3</SUB>的优化结果</STRONG></FONT></P>
<CENTER>
<TABLE cellSpacing=0 border=1>
<TBODY>
<TR>
<TD align=middle rowSpan=2><FONT face=宋体 size=2>初始值</FONT></TD>
<TD align=middle rowSpan=2><FONT face=宋体 size=2>测量值</FONT></TD>
<TD align=middle colSpan=5><FONT face=宋体 size=2>优化值</FONT></TD>
<TD align=middle rowSpan=2><FONT face=宋体 size=2>平均值</FONT></TD>
<TD align=middle rowSpan=2><FONT face=宋体 size=2>方差</FONT></TD></TR>
<TR>
<TD align=middle><FONT face=宋体 size=2>1</FONT></TD>
<TD align=middle><FONT face=宋体 size=2>2</FONT></TD>
<TD align=middle><FONT face=宋体 size=2>3</FONT></TD>
<TD align=middle><FONT face=宋体 size=2>4</FONT></TD>
<TD align=middle><FONT face=宋体 size=2>5</FONT></TD></TR>
<TR>
<TD align=middle><FONT face=宋体 size=2>φ<SUB>2</SUB>/(°)</FONT></TD>
<TD align=middle><FONT face=宋体 size=2>112.45</FONT></TD>
<TD align=middle><FONT face=宋体 size=2>112.49</FONT></TD>
<TD align=middle><FONT face=宋体 size=2>112.45</FONT></TD>
<TD align=middle><FONT face=宋体 size=2>112.47</FONT></TD>
<TD align=middle><FONT face=宋体 size=2>112.49</FONT></TD>
<TD align=middle><FONT face=宋体 size=2>112.47</FONT></TD>
<TD align=middle><FONT face=宋体 size=2>112.48</FONT></TD>
<TD align=middle><FONT face=宋体 size=2>0.014</FONT></TD></TR>
<TR>
<TD align=middle><FONT face=宋体 size=2>φ<SUB>3</SUB>/(°)</FONT></TD>
<TD align=middle><FONT face=宋体 size=2>230.23</FONT></TD>
<TD align=middle><FONT face=宋体 size=2>229.91</FONT></TD>
<TD align=middle><FONT face=宋体 size=2>229.89</FONT></TD>
<TD align=middle><FONT face=宋体 size=2>229.90</FONT></TD>
<TD align=middle><FONT face=宋体 size=2>229.86</FONT></TD>
<TD align=middle><FONT face=宋体 size=2>229.88</FONT></TD>
<TD align=middle><FONT face=宋体 size=2>229.89</FONT></TD>
<TD align=middle><FONT face=宋体 size=2>0.018</FONT></TD></TR>
<TR>
<TD align=middle><FONT face=宋体><FONT size=2>|C<SUB>1</SUB>|/10<SUP>-6</SUP></FONT></FONT></TD>
<TD align=middle><FONT face=宋体 size=2>3521</FONT></TD>
<TD align=middle><FONT face=宋体 size=2>891</FONT></TD>
<TD align=middle><FONT face=宋体 size=2>957</FONT></TD>
<TD align=middle><FONT face=宋体 size=2>1264</FONT></TD>
<TD align=middle><FONT face=宋体 size=2>1018</FONT></TD>
<TD align=middle><FONT face=宋体 size=2>984</FONT></TD></TR></TBODY></TABLE></CENTER>
<P>
<TABLE width="90%" border=0>
<TBODY>
<TR>
<TD>
<P><FONT face=宋体 size=2><STRONG>4 实验结果及分析</STRONG></FONT></P>
<P><FONT face=宋体><FONT size=2> 分别采用某文献推荐的φ<SUB>2</SUB>=84.24°,φ<SUB>3</SUB>=168.78°及本文中取K<SUB>m</SUB>=160时优化出的角度安置测头,对一个直径为110mm的标准件进行测量,并与Talyrond圆度仪测得结果相比较,见表2.</FONT></FONT></P>
<P align=center><FONT face=宋体 size=2><STRONG>表2 测量结果比较</STRONG></FONT></P></TD></TR></TBODY></TABLE></P>
<CENTER>
<TABLE cellSpacing=0 border=1>
<TBODY>
<TR>
<TD align=middle rowSpan=2><FONT face=宋体 size=2>方 法</FONT></TD>
<TD align=middle colSpan=4><FONT face=宋体 size=2>圆 度</FONT></TD>
<TD align=middle colSpan=2><FONT face=宋体 size=2>谐波误差</FONT></TD>
<TD align=middle colSpan=2><FONT face=宋体 size=2>圆周各点测量误差</FONT></TD></TR>
<TR>
<TD align=middle><FONT face=宋体 size=2>均值/μm</FONT></TD>
<TD align=middle><FONT face=宋体 size=2>改善程度/%</FONT></TD>
<TD align=middle><FONT face=宋体 size=2>σ/μm</FONT></TD>
<TD align=middle><FONT face=宋体 size=2>改善程度/%</FONT></TD>
<TD align=middle><FONT face=宋体 size=2>最大值/μm</FONT></TD>
<TD align=middle><FONT face=宋体 size=2>改善程度/%</FONT></TD>
<TD align=middle><FONT face=宋体 size=2>最大值/μm</FONT></TD>
<TD align=middle><FONT face=宋体 size=2>改善程度/%</FONT></TD></TR>
<TR>
<TD align=middle><FONT face=宋体 size=2>圆度仪</FONT></TD>
<TD align=middle><FONT face=宋体 size=2>1.253</FONT></TD>
<TD align=middle><FONT size=2> </FONT></TD>
<TD align=middle><FONT face=宋体 size=2>0.050</FONT></TD></TR>
<TR>
<TD align=middle><FONT face=宋体 size=2>传统三点法</FONT></TD>
<TD align=middle><FONT face=宋体 size=2>1.292</FONT></TD>
<TD align=middle><FONT size=2> </FONT></TD>
<TD align=middle><FONT face=宋体 size=2>0.035</FONT></TD>
<TD align=middle><FONT size=2> </FONT></TD>
<TD align=middle><FONT face=宋体 size=2>0.072</FONT></TD>
<TD align=middle><FONT size=2> </FONT></TD>
<TD align=middle><FONT face=宋体 size=2>0.250</FONT></TD></TR>
<TR>
<TD align=middle><FONT face=宋体 size=2>优化三点法</FONT></TD>
<TD align=middle><FONT face=宋体 size=2>1.282</FONT></TD>
<TD align=middle><FONT face=宋体 size=2>26</FONT></TD>
<TD align=middle><FONT face=宋体 size=2>0.026</FONT></TD>
<TD align=middle><FONT face=宋体 size=2>26</FONT></TD>
<TD align=middle><FONT face=宋体 size=2>0.024</FONT></TD>
<TD align=middle><FONT face=宋体 size=2>67</FONT></TD>
<TD align=middle><FONT face=宋体 size=2>0.160</FONT></TD>
<TD align=middle><FONT face=宋体 size=2>36</FONT></TD></TR></TBODY></TABLE></CENTER>
<P>
<TABLE width="90%" border=0>
<TBODY>
<TR>
<TD><FONT face=宋体 size=2> 结果表明:<BR> ① 传统三点法因某些级次的误差传递系数q<SUB>k</SUB>较大,读数误差被大大放大,从而使测量结果中有明显的高次谐波误差(见图2d与图3a).<BR> ② 传统三点法测得圆度以及各次谐波值的离散程度较大,单圈测量值与圆度仪测得值相差较大,必须取多圈平均才能求得较精确的结果.<BR> ③ 优化三点法圆度以及各次谐波值的测量不确定性明显降低,尤其高次谐波测量误差明显降低(见图3b).</FONT></TD></TR></TBODY></TABLE></P>
<P align=center><IMG src="http://www.chmcw.com/upload/news/RCL/13220_bcdmv1200832414033.gif"></P>
<P align=center><FONT face=宋体>(a)传统三点法结果</FONT></P>
<P align=center><IMG src="http://www.chmcw.com/upload/news/RCL/13220_8fzimy200832414050.gif"></P>
<P align=center><FONT face=宋体>(b)优化三点法测量结果</FONT></P>
<P align=center><FONT face=宋体><STRONG>图3 圆度和各次谐波幅值测量误差</STRONG></FONT></P>
页:
[1]