第39卷第8期 天津大学学报 Vo1.39 No.8 Aug.2006 2006年8月 Journal of Tianjin University 提高激光干涉测量系统精度的方法与途径 齐永岳,赵美蓉,林玉池 (天津大学精密测试技术及仪器国家重点实验室,天津300072) 摘要:激光干涉系统的测量精度会受到系统本身机械结构及光学元件布局的影响,因此采用了耦合差动干涉的 方法,使影响测量精度的各种因素尽量由干涉系统自身予以消除,提高了干涉仪的分辨率和稳定性.利用以Heyde— mann修正模型为基础的误差补偿方法,对干涉信号中存在的非正交误差、不等幅误差及直流电平漂移误差进行了 修正.根据多点拟合的原理,提出了减小运算量的新算法.设计了以C8051FO05单片机为核心的电路补偿模块,实 现了以上3种误差的实时补偿.实验结果表明:通过采取以上措施,该干涉测量系统在10 mm的测量范围内取得了 10~12 nm的测量精度. 关键词:激光干涉;测量精度;细分误差;误差补偿 中图分类号:TB96 文献标志码:A 文章编号:0493—2137(2006)08—0989—05 Methods of Improving Accuracy of Laser Interferometry System QI Yong—yue,ZHAO Mei—rong,LIN Yu—chi (State Key Laboratory of Precision Me ̄emem Technology and Instrument,Tianjin University,Tianjin 300072,China) Abstract:The accuracy of laser interferometry system are influenced by its mechanical construction and the distribution of the optical elements.A coupled differential interferometry method was adopted,which eliminated the affecting factors by the system itself.The resolution and the stability of the interferometer were improved.In addition,the three e ̄ors,namely,quadrature phase shift eⅡor,unequal amplitude of two signals,and zero offsets,were compensated by the Heydemann mode1.A fast algorithm Was proposed to speed up the calculation of these en’ors.The real time—error compensation was accomplished wih the C805 t1 F005 single chip microcom— puter.Then the measuring accuracy of lO—l2 nm was achieved within the range of 10 nun. Keywords:lser interferometaer;measuring accuracy;division enDfs;error compensation 单频激光干涉系统被广泛应用于长度测量领域, 其优点是测量范围广且测量精度高,但同时对测量环 1 干涉测量系统结构设计的优化 激光干涉测量系统的结构原理如图1所示.组成 境条件要求较高,其稳定性、分辨率和测量精度容易受 到环境温度、湿度、气压、振动以及仪器结构设计等因 素的影响 j.因此,提高激光干涉测量系统的精度要 从以下两方面人手:①合理设计机械结构及光学元件 布局,尽量使影响测量精度的各种因素由干涉仪自身 予以消除 ;②在干涉信号中还存在非正交误差、不 等幅误差及直流电平漂移误差 j,严重影响着细分精 度,需要对其进行补偿修正,以保证整套系统的测量 精度. 该系统的光学元件,除角锥棱镜2、5固定在导轨滑架 两端以外 ,其余均固定在横跨导轨的主支架上.从 图中可以看出,激光器发出的光束经分光镜1分成两 束光:一路反射光经过角锥棱镜2、3多次反射,平移一 段距离后返回分光镜1,形成一路测量光束a;另一路 透射光,经反射镜4反射,到达角锥棱镜5、6,经过多 次反射平移一段距离后回到反射镜4,再经反射回到 分光镜1,形成另一路测量光束b.两路光在分光镜1 收稿日期:2005—07—20;修回日期:2006—02—20. 基金项目:国家自然科学基金资助项目(50105016). 作者简介:齐永岳(1978一 ),男,博士研究生,qyy1978qyy@163.COIII. 维普资讯 http://www.cqvip.com
・990・ 天 津 大 学 学 报 第39卷第8期 重新汇合,产生干涉.干涉条纹信号由反射镜7反射, 被分光镜8分成两部分,由光电接收器D 、D 来接收 转换,获得相位差为90。的两路按正弦规律变化的电 信号. D l一分光镜;2,3,5,6一角锥棱镜;4,7一反射镜;8一分光镜 图1干涉光路 Fig.1 lnterferometric optical path 当角锥棱镜2和5所在的导轨滑架向左移动△L 时,测量光束a的光程增加了4AL,而测量光束b的光 程减少了4AL,这样光程差 D叩=4AL一(一4AL)=8AL (1) 因此可以产生A。/(8n)的光学细分,即导轨滑架每移 动A。/(8n),光电接收器输出信号就会变化一个周期, 于是便得到了这套干涉系统测长的基本公式,即 AL=N・A。/(8n) (2) 式中:n为测量时的空气折射率;A。为真空中的激光波 长;Ⅳ为干涉条纹明暗变化次数. 如果用计算机对干涉信号再进行200细分,则该 干涉测量系统的理论分辨率 SRE=ALl(N・200)=A。/(1 600 n) (3) 与普通的迈克尔逊干涉仪相比,该系统的结构布 局具有特点:①将长度信息直接转化为多倍光程差,即 通过光学细分的方法,使干涉系统的分辨力提高了4 倍,这将有助于提高干涉系统的抗干扰能力和测量精 度;②干涉系统受到的扰动主要来自于环境温度、湿度 以及气压等引起的空气折射率的变化 J.在图1所示 的系统中,空气折射率变化 引起的光程变化为 6。=6 ・(P。一P6) (4) 式中:P。为测量光束a的光程;P 为测量光束b的光 程.因光路布局对称,即P。一P ,这一影响将趋近于 零.由此可见,这套干涉系统具有良好的耐环境变化的 特性. 2干涉信号的误差补偿 2.1干涉信号的处理 由于干涉条纹对比度不为1并受到外界干扰,从 光电接收器输出的信号中除含有交流成分外,还含有 直流成分,所以必须首先调整信号中的直流分量.为此 采用如图2所示的直流电平调整电路.通过改变放大 器同相输人端的电压,达到消去直流电平的目的. 图2直流电平调整 Fig.2 Zero offset adjustment 为了消除高频噪声等干扰信号,对输出信号还必 须进行滤波处理.为此,采用如图3所示的二阶压控电 压源切比雪夫型低通滤波器.取截至频率 =10 kHz, 用查表归一快速设计方法确定滤波器的参数为R = 9.1 kQ,R2=25.4 kQ,Cl=330 PF,C=1 000 PF. C 图3低通滤波 Fig.3 Low pass f'dter 2.2信号三差对细分精度的影响 在增量式计量仪器中,一般以相位差为90。的两 路信号进行辨向计数和细分 J.而在这两路信号中普 遍存在着非正交误差、不等幅误差和直流电平漂移误 差(以下简称为三差),这些误差的存在严重地影响了 测量结果的精度. 在理想情况下,一对等幅正交的正余弦信号,其李 沙茹图形为理想的圆.其数学模型可表示为 l【uH2=R sin 0 1 c0s 0 (5) 式中0为信号的相位角.而实际信号中总是不可避免 地存在上述3种误差,同时考虑这3种误差的影响,信 号的数学模型则应该为 f【H/t 2d=Hd (/ +t2 COSp Ot—Hl sin )/G+q (6) 式中iot为非正交误差;G为两路信号幅值之比;P和q 分别为两路信号的直流电平.当对/t 、/Z2d进行细分 时,相位角的计算值为 0=arctan(Uld/u2d) (7) 由于H 和H 偏离了理想信号,必然使0的计算值引 维普资讯 http://www.cqvip.com
2006年8月 齐永岳等:提高激光干涉测量系统精度的方法与途径 ・99l・ 入误差.由微分定理可知0的计算误差为 68=(M16u2一u28 1)/R (8) =arcsin(一C 、 G= ) 由式(5)和(6)得 r6u1=p R—sin 0(cos a-G)-R cos 0 sin a p (2BD—EC)/(C +4B) q=(一2E—DC)/(C +4 ) +g (15) :R:[ 差的影响. , ]÷ (9) 将式(9)代人式(8),得 0。一 。 这样,对信号实测后就可以通过误差补偿消除三 r 1 ld—P sin 0 COS 0(COS Ot—G)一COS 0 sin Ot. G (COS 0・q—sin 0・P)/R 6 (1一G)sin 20/2一Ot COS + (10) {, M 1 sin +G(M2d—g) l//,,=一 (16) 通常G一1, 一0,忽略二次以上误差,上式简化为 L ‘ COS OL 多元线性回归方程根据随机变量Y及玑个自变 量 , 一, 一1的n组观测数据( ∽…, _1 .) sin arctan詈) (11) (i=0,1,…,n一1),用式(17)所示的线性表达式对观 测数据进行回归分析. Y=Ⅱ0 0+Ⅱ1 1+…+Ⅱm一1 m一1+am (17) 从式(11)中可以看出,不等幅引入的误差最大值 为(1一G)/2 rad,最大值出现的位置为0=,n'/4,3,n'/4, …;非正交误差引入的误差最大值为Ot(rad),最大值 出现的位置为0=0,1T,…;直流电平漂移引入的误差 最大值办 +g /R,出现的位置随p和g的相对大小 式中 ,…,a 为回归系数.为确定回归系数,通常采 用最小二乘原理.为使 而定.由于3种误差的峰值彼此错开,因此三差的综合 影响小于各自的最大误差之和. 由此可见,信号三差的存在严重影响着细分精度。 通过改进硬件也难以彻底消除,因此采用了Heyde— mann修正模型对非正交误差、不等幅误差和直流电平 漂移误差进行修正补偿 ]. 2.3信号三差的补偿原理 (18) 达到最小,回归系数应满足方程组 (CCT)(00…am) =C(y0…Y ) 一 (19) 其中 X00 10 两路理想的正余弦信号的矢量合成图,即李沙茹 图形是一个圆.当存在上述3种误差时,李沙茹图形是 个椭圆.椭圆方程与三差的关系可描述为 一C= : ● (20) m一1.0 1 ( p)z+『 d二 1 :Rz (12) 把式(17)与式(13)比较可知: Y =M21d2(i) (i) Ⅱ0=B 0i:M2dⅡ1=C 式中M 和 分别为带有三差的实际信号.可将上式 改写为 21 :u1d(i)M2d(i) Ⅱ2=D 2 :M1d(i) Ⅱ =E (21) M1d=A+ M +cM1dM2d+DM1d+EM2d A=R COS d—P 一G q 一2Gpq sin B:一G (13) 3 :M2d(i) Ⅱ4=A 由此生成矩阵C,然后通过矩阵运算生成实际的 方程组(19);通过求解这个方程组求出系数A、 、C、 (14) C=一2G sin Ot D=2p 4-2Cq sin D、E.由式(15)从系数A、 、C、D、E,求得 、G.p、q. 2.4实时补偿的实现 E=2G q+2cp sin Ot 由于在测量过程中采用的是逐点显示方式,即在 此式符合多元线性回归的形式,通过运算可以求 出系数A、 、C、D、E的最佳值,根据计算结果反算出 0、G,p、g的大小. 导轨滑架运动的同时把它的位移值显示出来,因此采 用如下算法:当导轨滑架开始运动时,连续采集17,组数 据,用于首次三差补偿.其后每次只采1组数据,用最 维普资讯 http://www.cqvip.com
・992・ 天 津 大 学 学 报 第39卷第8期 新采集的这组数据代替前面It组数据中最早采集的1 组数据,将新组成的 组数据进行三差拟合,然后显示 最新数据对应的位移值.如此周而复始.因为显示时问 问隔短,在短时问内要完成三差拟合、补偿和显示,所 以要求改进算法,减小运算量.设式(19)中A:CC , D=C(Y。…Y ) ,则A为(m+1)×(m+1)阶矩 阵,D为(m+1)行矩阵.矩阵内各元素分别为 =∑ √ t=l aji ai,,j (22) i=1,2,…,m+1 J=i,i+1,…,,n+1 =∑xt:1 t,jY =1,2,…,m+1 (23) 从中可以看出线性回归计算量主要由两部分组 成:①计算系数矩阵A和D;②解方程组(19).前者需 要It(m+1) /2次乘除法运算,后者解方程,若采用高 斯消去法约需要(m+1) /3次运算.由于It>>m,因此 计算系数矩阵A、D是线性回归运算的最主要部分.只 有减小求系数矩阵的运算量,才能提高线性回归的速 度.当导轨滑架运动时,计算机不断采集新数据,使用 于评定三差的数据依次被刷新.当n组数据中的1组 发生变化时,根据式(22)和(23),系数矩阵A、D的每 一个元素都将发生变化.由于其他(n一1)组数据并没 有发生变化,因此不必按照式(22)、(23)完全重新计 算,而只需根据其变化量计算新值.设n组数据中的第 P组被替代,则系数矩阵新旧值有如下关系: at tj=at i一.Xpi.XPj+ P Pj ,D i=Oj -Xp|,j P ) , , 、 0 Ji a i,,j i=1,2,…,,n+1 =i,i+1,…,,n+1 式中上角标 表示第P组数据被替代后对应的新值(未 注角标的表示被替代前的数据).式(24)称为系数矩 阵的迭代运算法. 比较式(22)、(23)与式(24)发现,迭代法运算量 为(m+1) /2,比直接算法减小运算量n倍.加上解方 程组(19)的运算量,每一次回归只需约80次乘除运 算,基本上可以满足系统逐点显示的需要. 2.5 补偿电路设计及测试结果 采用美国Cygnal公司生产的C8051F005单片机 完成信号的采集和数据处理工作.电路原理框图如图 4所示.两路带有三差的实际信号u u 经 C8051F005单片机的片内ADC转换成数字量,按照 Heydemann算法进行数据处理,得出消除三差以后的 信号u ,、 ,:,再经片内DAC转换成模拟量输出 图4电路原理 Fig.4 Circuit diagram 测试结果如图5所示.从图中可以明显看出,信号 三差得到了有效补偿,从而保证了细分精度. (a)输入信号Mld和M2d (b)输出信号H l和H 2 图5补偿效果 Fig.5 Compensation result 3 系统实验 采用与精度水平相近的仪器作比对的方法对该干 涉测量系统的精度进行了检测.实验在常温下进行,选 用天津大学研制开发的JDC.2000型高精度电容测微 仪和本干涉系统进行比对.JDC精密电容式测微仪是 一种非接触式精密测量仪器,其分辨率为4 nm. 测量结果如表1所示.其中, 表示电容测微仪的 输出值, 表示本干涉系统的测量值. 维普资讯 http://www.cqvip.com
2006年8月 齐永岳等:提高激光干涉测量系统精度的方法与途径 ・993・ 表1单次比对结果 Tab.1 Contrast result ,、 上/ m ,、 上/ m 0.552 0.605 5.316 5.561 0.908 0.973 5.804 6.O66 1.364 1.436 6.288 6.576 1.936 2.029 6.794 7.113 2.412 2.5l9 7.260 7.598 2.832 2.976 7.696 8.048 3.476 3.633 8.188 8.549 3.922 4.120 8.720 9.O96 4.548 4.754 9.112 9.500 4.956 5.185 9.576 9.981 注:线性回归方程为L=1.033 282 U+0.002 154;标准误差 = 0.010. 从表1中可以看出, 与 呈近似线性关系,于 是可以通过一元线性回归求出 与 之间的线性回 归方程.首先假设 是没有误差或误差可以忽略,将 所有误差都归结在L方向,此时 对 的回归方程为 L=a・U+b (25) 接下来考察另一种极端情况,即假设 是没有误 差或误差可以忽略,而将所有误差都归结在 方向, 则 对 的回归方程为 U=a0・L+b0 (26) 将上式改写为 , L : . 一 f27) n0 口0 由式(25)和(27)发现,通过一批点能够作2条回 归直线,在这种情况下,则可以采用2条回归线的平均 线作为测量的回归方程.在确定了 与 的线性回归 方程之后,便可求出测量数据的标准不确定度 等精 度参数. 4 结语 通过与电容测微仪在不同时间内进行的5次比对 测试,得出5次测量的标准不确定度(1 )为10~12 由此可以看出:优化干涉仪机械结构及光学元件 布局,可使干涉系统具有良好的耐环境变化挣l生;对信 号三差进行动态补偿则有效地保证了细分精度.通过 采用以上措施切实提高了干涉系统的测量精度. 参考文献: [1] Brand U,Herrmann K.A laser measurement system for the high precision calibration of displacement transducers[J]. Measurement Science and Technology,1996,7(6):911m 917. 『2]Lee C K,Wu T W.Diferentila laser interferometer for nanometer displacement measurements[J].AIAA Journal, 1995,33(9):1675--1680. [3] Wang C H,Augousti A T,Mason J.Real time evaluation and correction of nonlinear errors in single ̄equeney inter. ferometers[J].Transactions ofthe Isntitute ofMesaureemnt and Control,2000,22(5):405—1412. [4] Zhao Meirong,Qi Yongyue,Lin Mingchun.et a1.Practical system for long—range nano—positioning and nano—measure— ment[C]//2nd Internatoinal Symposium on Isntrumentation Science and Technology.Jinan,China,2002:105—109. [5] 黎永前,朱名铨.光干涉技术在纳米精度测量中的应用 及发展[J].航空精密制造技术,2002,38(2):32—36. Li Yongqian,Zhu Mingquan.Optical interferometry techno1. ogy and its development in nanometric accuracy measurement [J].Aviatoin Precision Manufacturing Technology.2002. 38(2):32--36(in Chinese). [6]韩旭东,艾 华.共光路移相单频激光干涉测长系统 [J].光学技术,2004,30(2):195—198. Han Xudong,Ai Hua.Commonpath and phaseshifting single frequency laser interferometer for length measurement[J]. Optical Technique,2004,30(2):195—198(in Chinese). [7]Heydemann Peter L M.Determination and correction of quadrature firnge measurement errors in interferometers[J]. Appleid Optcis,1981,20(19):3382—3384. [8] 陈卫兵.从系统结构看Cygnal C8051FXXX系列单片机 [J].微电子技术,2003,31(6):40 1. Chen Weibing.Seeing Cygnal C805 1 FXXX single chip processor from system construction[J].Microelectronic Tech— onlogy,2003,31(6):4O 1(in Chinese).
因篇幅问题不能全部显示,请点此查看更多更全内容