張大軍,劉合朋,劉 軍,閔躍軍,鐘正虎
(北京航天控制儀器研究所,北京100039)
基于雙頻激光干涉儀的角速率檢測系統(tǒng)精度分析
張大軍,劉合朋,劉 軍,閔躍軍,鐘正虎
(北京航天控制儀器研究所,北京100039)
為了能夠精確檢測到測試轉(zhuǎn)臺在低速運行下的角速率,設(shè)計了一種基于雙頻激光干涉儀的角速率精度檢測系統(tǒng)。首先,對本檢測系統(tǒng)各項誤差進行機理分析。然后,綜合各項誤差項建立總體相對誤差模型。最終,通過仿真分析得到各個誤差因子對相對誤差項的影響,為保證檢測系統(tǒng)相對精度提供理論依據(jù),能夠滿足對高精度慣性器件測試轉(zhuǎn)臺低速段0.0001(°)/s~1(°)/s速率檢測要求。
測試轉(zhuǎn)臺;低速;角速率測量;雙頻激光干涉;精度分析
Abstract:In order to accurately measure angular rate of the test turn?table running at low speeds,an angular rate measurement system based on dual?frequency laser interferometer is designed.Therefore,themechanism of the errors of the detection system is analyzed.Then,itswhole errormodel integrating relative error terms is established.Finally,calculated by simulation,the influence of error factors to relative error terms is obtained.The theory basis for ensuringmeasurement systematic accuracy is provided.At last,when high precision inertial components test turntable is running at low speed of 0.0001(°)/s~1(°)/s,angular ratemeasurement is achieved.
Key w ords:test turn?table;low speed;angular ratemeasurement;dual?frequency laser interferometer;precision a?nalysis
轉(zhuǎn)臺在慣性器件測量領(lǐng)域發(fā)揮著重要作用,其中角速率精度又是重要指標(biāo)。速率精度可分為低速0.0001(°)/s~1(°)/s、中速1(°)/s~10(°)/s、高速10(°)/s~200(°)/s這3段。高精度慣性器件測試轉(zhuǎn)臺精度要求達到5×10-5(0.0001(°)/s~1(°)/s,1°平均)、5×10-6(1(°)/s~10(°)/s,10°平均)、5× 10-7(10(°)/s~200(°)/s,360°平均)的國際領(lǐng)先水平。采用傳統(tǒng)的測量方法測量轉(zhuǎn)臺低速階段角速率時,比如光電軸角編碼器法,要想滿足其測量相對精度,光電編碼器的絕對精度要到0.06″,目前技術(shù)達不到。因此,需要研制新的速率精度檢測方法。本文利用雙頻激光干涉儀與光學(xué)差動結(jié)構(gòu)組合測量角速度的特點[1?2],設(shè)計了一種基于雙頻激光干涉儀的轉(zhuǎn)臺低速速率精度檢測系統(tǒng),并通過反正弦解算的方法,低速速率檢測相對精度可達1.7×10-5,解決了無法對轉(zhuǎn)臺高精度低速率進行有效檢測的技術(shù)難題。因此,本文要分析基于雙頻激光干涉儀的角速率檢測系統(tǒng)精度,為保證檢測系統(tǒng)相對精度提供理論依據(jù)。
轉(zhuǎn)臺角速率測量基本原理如式(1)所示:
式中,ω為角速率,Δα為角位移,ΔT為轉(zhuǎn)動角位移Δα所經(jīng)歷的時間。時間間隔ΔT可采用高精度時鐘進行采樣滿足精度要求,角位移Δα的高精度測量可根據(jù)雙頻激光干涉儀和光學(xué)差動結(jié)構(gòu)共同實現(xiàn)。
雙頻激光干涉角速率檢測系統(tǒng)工作原理如圖1所示。雙頻激光發(fā)生器1產(chǎn)生2個相互垂直的具有不同頻率的偏振光,頻率分別為f1和f2,一部分偏振光在內(nèi)部實現(xiàn)干涉,被內(nèi)部光電接收器檢測到基準信號(f1-f2)。同時,一部分光通過偏振分光鏡4,按照偏振方向不同,兩偏振光分離,f1光透過偏振分光鏡4,通過反射鏡5兩次反射,投射到正弦臂8上的一個角錐棱鏡7;f2光從偏振分光鏡4反射到另一反射鏡5,再反射投射到正弦臂8上的另一個角錐棱鏡7。
式中,N1為在t時間內(nèi)由偏振光f1引起的周期數(shù),N2為在t時間內(nèi)由偏振光f2引起的周期數(shù)。
角錐棱鏡沿光路移動L表達式為:
式中,λ1為偏振光f1的波長,λ2為偏振光f2的波長,c為光速。
根據(jù)式(2)和式(3)可推導(dǎo)得出沿光路移動L與周期數(shù)N之間的關(guān)系:
根據(jù)正弦定理得:
式中,α為轉(zhuǎn)動角位移,R為角錐棱鏡光學(xué)中心到旋轉(zhuǎn)中心的距離。
根據(jù)式(4)和式(5),λ1與λ2兩波長近似相等,記作λ。因此,可得到轉(zhuǎn)動角位移α與周期數(shù)N的關(guān)系為:
式(6)為基于雙頻激光干涉儀的轉(zhuǎn)臺角速率檢測系統(tǒng)角位置數(shù)學(xué)模型,可以根據(jù)該數(shù)學(xué)模型精確反正弦解算出轉(zhuǎn)臺的角位移。同時,使用具有高精度晶振時基的多功能數(shù)字相位卡進行高速數(shù)據(jù)實時采樣,數(shù)據(jù)處理后可得到角速率信息。
2.1 正弦零位誤差引起的相對誤差
由于正弦零位位置是個理想位置,在實際工程中,是存在正弦零位誤差的,如圖2所示。
由圖2可知,實際角位移α~與周期數(shù)N的關(guān)系為:
式中,ε為正弦零位誤差。
根據(jù)式(6)和式(7),得到實際角位移α~與測量角位移α、正弦零位誤差ε之間的關(guān)系:
可得出正弦零位誤差引起的相對誤差項ΔJ:
2.2 正弦臂差動結(jié)構(gòu)偏心安裝誤差引起的相對誤差
正弦臂安裝誤差存在以下幾種情況,分別如圖3~圖5所示,其誤差對檢測系統(tǒng)相對精度的影響程度是不一樣的,針對這些情況,進行精度分析。
如果正弦臂差動結(jié)構(gòu)偏心安裝如圖3所示,當(dāng)正弦臂中心與轉(zhuǎn)臺的旋轉(zhuǎn)中心沿著徑向方向有偏置時,實際角位移α~可表示為:
式中,Δ為偏心距離,L1、L2分別為兩角錐棱鏡沿光路方向移動的距離。
根據(jù)式(2)和式(10)整理得到:
通過式(11),采用平均真空波長為633nm的雙頻激光頭,拍頻為1MHz~4MHz,也就是說λ1+λ2?λ1-λ2,已經(jīng)滿足系統(tǒng)的相對精度的要求。因此,2R(λ1+λ2)?2(λ1-λ2)Δ,使得該方式偏心安裝引起的系統(tǒng)測量誤差可以忽略不計。
如果正弦臂差動結(jié)構(gòu)偏心安裝如圖4所示,當(dāng)正弦臂中心與轉(zhuǎn)臺的旋轉(zhuǎn)中心沿著周方向有偏置時,實際角位移α~可表示為:
式中,θ為由于偏置Δ所產(chǎn)生的夾角。
根據(jù)式(12)和式(13)整理得到:
通過式(14),采用平均真空波長為633nm的雙頻激光頭,拍頻為1MHz~4MHz,也就是說λ1≈λ2。因此,可以得出偏心安裝誤差引起的系統(tǒng)測量誤差可忽略不計。
如果正弦臂差動結(jié)構(gòu)與轉(zhuǎn)臺平面有安裝夾角,如圖5所示。實際角位移可通過式(15)得到:
式中,β為正弦臂差動結(jié)構(gòu)與轉(zhuǎn)臺平面有安裝夾角。
根據(jù)式(6)和式(15),得到實際角位移α~與測量角位移α、平行度β之間的關(guān)系:
可得出正弦臂差動結(jié)構(gòu)偏心安裝誤差引起的相對誤差項ΔA:
在安裝正弦臂時,盡量做到正弦臂差動結(jié)構(gòu)與轉(zhuǎn)臺平面平行。本文采用現(xiàn)場臂長標(biāo)定方法,在保證檢測系統(tǒng)精度的同時,對正弦臂的安裝精度要求降低。
2.3 正弦臂臂長相對誤差
環(huán)境溫度的變化使正弦臂長發(fā)生變化,進而影響測角精度。同時,安裝時正弦臂不平行于臺面可以看作影響臂長來影響測角精度。因此,如何保證現(xiàn)場的臂長精度是十分必要的。本文采用現(xiàn)場臂長標(biāo)定方法,即在檢測轉(zhuǎn)臺低速角速率之前,先使用轉(zhuǎn)臺位置控制模式,對光學(xué)正弦臂臂長進行標(biāo)定。利用轉(zhuǎn)臺的絕對位置精度高,標(biāo)定臂長時增加轉(zhuǎn)臺轉(zhuǎn)過角度。
根據(jù)正弦定理得:
式中,?為標(biāo)定時轉(zhuǎn)臺的角位移。
根據(jù)式(8),?對R求微分,整理得到:
式中,?為標(biāo)定時轉(zhuǎn)臺的角位移,Δ?為轉(zhuǎn)臺角位移絕對精度,ΔR是標(biāo)定的正弦臂長誤差。
因此,得到現(xiàn)場標(biāo)定的臂長相對誤差項ΔR與標(biāo)定時選擇的轉(zhuǎn)臺轉(zhuǎn)角和轉(zhuǎn)臺絕對精度的關(guān)系。
2.4 角錐棱鏡的反射光學(xué)特性及其引起的相對誤差
角錐棱鏡又稱角體合作目標(biāo)[3],是由3個互相垂直的平面組成3個直角。頂點是3個互相垂直棱的交點,底部切割成一等邊三角形。直角表面為反射面,角體合作目標(biāo)的表面不平度、直角誤差均影響反射光束的發(fā)散角和反射光束的遠場衍射花樣[4?6]。下面對角體合作目標(biāo)的有關(guān)性質(zhì)進行初步的精度分析。
角體合作目標(biāo)的內(nèi)部光線軌跡:光線從角體合作目標(biāo)的底面入射后,入射線將在三角面上依次反射,出射線將平行于入射線射出,幾何相的對稱中心是三角面頂點。
根據(jù)光學(xué)反射定律,入射光經(jīng)過角體合作目標(biāo)的3個直角反射面全反射恰好可以等效為光線從入射點關(guān)于頂點的對稱點入射,通過角體合作目標(biāo)內(nèi)部射出后方向正好與入射光方向相反[7?9]。因此,角錐棱鏡等效光路圖可以看成如圖6所示,將檢測系統(tǒng)中正弦臂兩端的兩角錐棱鏡展開成平板玻璃,并分別建立坐標(biāo)系,令正弦臂臂長O1O2為2R。
當(dāng)正弦臂在水平面內(nèi)繞正弦臂回轉(zhuǎn)中心M旋轉(zhuǎn)一個微小角度α?xí)r,根據(jù)正弦原理,光程差H的表達式可表示為:
正弦原理沒有考慮正弦臂兩端角錐棱鏡的安裝和制造不對稱引入的光程差,當(dāng)考慮其不對稱性時,實際測試中的光程差ΔI可根據(jù)式(21)得到:
式中,ΔH為正弦臂兩端角錐鏡不對稱性引入的光程差。
經(jīng)分析,正弦臂兩端角錐鏡的不對稱性可以歸納成2種情況:
1)兩角錐棱鏡均為標(biāo)準的相同尺寸的角錐鏡,但兩者安裝在鏡組中時前表面不平行,存在微小夾角ε;
2)兩角錐棱鏡均為標(biāo)準角錐棱鏡,但兩者厚度不等,分別記為Ha和Hb,表現(xiàn)為將角錐棱鏡展開成平板玻璃以后的厚度不相等。
圖6是初位置角錐棱鏡等效光路,激光器發(fā)出的光束垂直于正弦臂中角錐棱鏡前表面入射。但是由于兩角錐棱鏡安裝時前表面不平行,當(dāng)入射光束垂直于角錐鏡a的前表面時,必然與角錐鏡b前表面存在一個夾角ε。
兩角錐棱鏡的光路長度分別為:
式中,n為角錐棱鏡折射率,S為O1O2距基準面的初始距離。推導(dǎo)出兩光路的光路差I(lǐng)0為:
同理,當(dāng)正弦臂繞回轉(zhuǎn)中心M轉(zhuǎn)過微小角度α后,角錐棱鏡與激光器入射光束的相對位置發(fā)生變化。圖7為正弦臂轉(zhuǎn)過α角后兩角錐棱鏡的等效光路圖。經(jīng)分析,正弦臂轉(zhuǎn)過α角后,光束通過兩角錐棱鏡的光路長度分別為:
式中,S1為轉(zhuǎn)動α角后O1距基準面的距離,S2為轉(zhuǎn)動α角時O2距基準面的距離。
推導(dǎo)出兩光路的光路差I(lǐng)′為:
可看出2(S2-S1)=4R sinα,則正弦臂轉(zhuǎn)過α角后,由雙頻激光干涉測量系統(tǒng)測得的光程差ΔI為:
因此,根據(jù)式(20)、式(21)和式(28)可計算得到角錐棱鏡不對稱性引起的相對誤差項ΔZ:
根據(jù)式(1)可以得出,角速率檢測系統(tǒng)誤差由測角誤差和采樣時鐘誤差共同決定??捎檬剑?0)表示:
式中,ΔT為采樣周期的相對誤差,Δα為角位移α相對誤差。
將式(6)全微分,略去高次項可得:
式中,Δα為轉(zhuǎn)角變化量,Δλ為波長不確定量,ΔN為周期數(shù)不確定量,ΔR為臂長不確定量。
其中,時鐘基準誤差項ΔT、電氣相對誤差項ΔN和光波波長相對誤差項Δλ都可以做得很小,滿足檢測系統(tǒng)相對精度的要求[10]。因此,本文重點對其他重要誤差項進行了機理研究,并進行了仿真計算。
使用Matlab仿真分析各個誤差因子對檢測系統(tǒng)各個誤差項的相對精度影響。
圖8是正弦零位誤差項ΔJ跟隨正弦零位絕對誤差和轉(zhuǎn)臺轉(zhuǎn)角的變化曲面。由圖8可知,正弦零位誤差控制在40″范圍內(nèi),測量范圍在±5°以內(nèi),正弦零位引起的相對誤差項ΔJ可控制在1×10-5范圍內(nèi)。
通過分析可知,采用光學(xué)差動結(jié)構(gòu)能夠減小安裝誤差對測量系統(tǒng)精度的影響。各類安裝誤差影響最大的是正弦臂差動結(jié)構(gòu)與轉(zhuǎn)臺平面平行度,其他的安裝誤差都能夠滿足1×10-5相對精度。因此,針對該誤差因子進行仿真分析。圖9是相對誤差項ΔA跟隨轉(zhuǎn)臺轉(zhuǎn)角和平行度β的變化曲面。由圖9可知,正弦臂差動結(jié)構(gòu)與轉(zhuǎn)臺平面平行度控制在0.3°范圍內(nèi),測量范圍在±5°以內(nèi),正弦臂安裝誤差引起的相對誤差項ΔA可控制在1×10-5范圍內(nèi)。
本文提供現(xiàn)場臂長標(biāo)定方案,利用轉(zhuǎn)臺角位置的絕對精度高的優(yōu)勢來提高正弦臂臂長的相對精度,同時對正弦臂的安裝精度要求降低。通過Matlab計算繪制現(xiàn)場標(biāo)定的臂長相對誤差跟隨標(biāo)定時所選的轉(zhuǎn)臺轉(zhuǎn)角和轉(zhuǎn)臺位置精度關(guān)系變化曲面,如圖10所示。通過分析可知轉(zhuǎn)臺位置精度在0.3″內(nèi),現(xiàn)場標(biāo)定的臂長相對精度 ΔR可控制在1×10-5范圍內(nèi)。
由于正弦臂兩端角錐鏡的不對稱性造成的光程差對檢測系統(tǒng)相對精度具有一定影響,配置各類參數(shù)后(R=32mm,n=1.4,Ha=10mm,α=1°),用Matlab仿真繪制相對誤差項ΔZ跟隨厚度差和兩錐面的夾角變化曲面,如圖11所示。分析可知,角錐棱鏡的加工工藝誤差以及在正弦臂中的安裝偏差都會對雙頻激光干涉正弦測角帶來測量誤差。由圖11可知,安裝偏差(兩角錐棱鏡前表面不平行度)在20″范圍內(nèi),制造工藝誤差(兩角錐棱鏡的厚度差)在0.1mm范圍內(nèi),相對誤差項ΔZ可控制在1×10-5內(nèi)。
本文建立了基于雙頻激光干涉儀的角速率檢測系統(tǒng)總體誤差模型,分析各個誤差因子與相對誤差項的關(guān)系,并利用Matlab對其進行了誤差仿真分析,可以得到以下結(jié)論:檢測系統(tǒng)測量范圍在±5°內(nèi),正弦零位誤差控制在40″范圍內(nèi),正弦零位引起的相對誤差項 ΔJ可控制在1×10-5范圍內(nèi);轉(zhuǎn)臺位置精度在0.3″內(nèi),添加標(biāo)定時轉(zhuǎn)臺轉(zhuǎn)角選擇5°,現(xiàn)場標(biāo)定的臂長相對精度ΔR可控制在1×10-5范圍內(nèi);安裝偏差(兩角錐鏡前表面不平行度)在20″范圍內(nèi),制造工藝誤差(兩角錐鏡的厚度差)在0.1mm范圍內(nèi),相對誤差項ΔZ可控制在1×10-5內(nèi)。綜上所述,該角速率檢測系統(tǒng)相對精度可以達到1.7×10-5,滿足了對高精度慣性器件測試轉(zhuǎn)臺低速段0.0001(°)/s~1(°)/s速率檢測要求。同時,本文的誤差分析也為進一步提高該檢測系統(tǒng)精度提供了理論依據(jù)。
[1]羨一民.雙頻激光干涉儀的原理與應(yīng)用(一)[J].工具技術(shù),1996,30(4):44?46.XIAN Yi?min.Principle and application of dual frequency laser interferometer(1)[J].Tool Technology,1996,30(4):44?46.
[2]羨一民,俞成棟,于光華,等.雙頻激光干涉儀在機床檢定中的應(yīng)用[J].制造技術(shù)與機床,1995(5):8?11+3.XIAN Yi?min,YU Cheng?dong,YU Guang?hua,et al.Ap?plication of dual frequency laser interferometer in machine checking[J].Manufacturing Technology& Machine Tool,1995(5):8?11+3.
[3]鐘聲遠.關(guān)于角體合作目標(biāo)性質(zhì)的理論計算[J].激光與紅外,1973,3(12):1?15.ZHONG Sheng?yuan.Theoretical calculation of the nature of the cooperative target[J].Laser&Infrared,1973,3(12):1?15.
[4]Zurasky JL.Cube corner retroreflector test and analysis[J].Applied Optics,1976,15(15):445?452.
[5]何川,林家明,鄒桂蘭,等.雙頻激光干涉測角中角錐鏡的不對稱性對測角精度影響的分析[J].光學(xué)技術(shù),2007,33(S1):289?290+292.HE Chuan,LIN Jia?ming,ZOU Gui?lan,et al.Effect of optical corner cubes'asymmetry on angle measurement[J].Optical Technique,2007,33(S1):289?290+292.
[6]JJG 998?2005,激光小角度測量儀[S].北京:中國標(biāo)準出版社,2005.JJG 998?2005,Laser small angle measuring instrument[S].Beijing:China Standard Press,2005.
[7]Huang L H,Li S,Zhou H.A method for measuring the active reflecting area of the cube corner prism and its error analysis[J].Journal of Geomatics,2005,30(1):45?47.
[8]Chen H,Tan J B.Effect of beam incident angles on tri?hedral corner reflector laser ranging[J].Journal of Optoe?lectronics Laser,2006,17(8):986?988.
[9]Ullrich D,Lubcck E,Wetzig V.Static and dynamic angular measurements with a ring laser gyro[C].Symposium Gyro Technology,Sutttgart,Germany,1988:150?159.
[10]何京微,許素安,謝敏,等.基于LabVIEW的激光雙頻外差干涉納米位移測量系統(tǒng)[J].測控技術(shù),2014,33(5):31?34.HE Jing?wei,XU Su?an,XIE Min,et al.A heterodyne interferometric nanometric displacement measurement system based on LabVIEW [J].Measurement&Control Technology,2014,33(5):31?34.
Precision Analysis on Angu lar Rate M easurem en t System Based on Dual?frequency Laser In terferom eter
ZHANG Da?jun,LIU He?peng,LIU Jun,MIN Yue?jun,ZHONG Zheng?hu
(Beijing Institute of Aerospace Control Devices,Beijing 100039)
U666.1
A
1674?5558(2017)03?01340
10.3969/j.issn.1674?5558.2017.04.012
張大軍,男,碩士,精密儀器及機械專業(yè),研究方向為精密測量技術(shù)。
2016?11?25
國家重大科學(xué)儀器設(shè)備開發(fā)專項(編號:2013YQ310737)