徐志剛，白鑫林，王軍義，王　昊

(1．中國科學(xué)院沈陽自動化研究所，遼寧沈陽　110016；2．中國科學(xué)院大學(xué)，北京　100039)

0　引言

目前，航天器平面運動測量多采用慣性導(dǎo)航陀螺及視覺等測量方式來實現(xiàn)[1]，但是陀螺以及視覺測量裝置測量頻率較低，無法用于閉環(huán)反饋控制?，F(xiàn)開發(fā)的測量頻率較高的平面運動測量系統(tǒng)多為利用角度測量裝置以及距離測量裝置通過機構(gòu)組合而構(gòu)成的接觸式測量系統(tǒng)，在一些運動精度要求較高，且不允許有運動干擾的系統(tǒng)中，測量系統(tǒng)自身機構(gòu)的運動會給被測系統(tǒng)運動帶來影響[2]，因此迫切要求一種新型的非接觸式平面運動測量裝置。

針對這一情況，文中提出一種新型非接觸式并聯(lián)平面激光測距傳感器，通過組合測量與計算求解，得到模樣件的位姿信息。激光測距傳感器具有較高的測量精度[7]，但是影響總體精度的因素也有很多[3]，對主要的8個影響因素，用蒙特卡洛法進行了相同概率分布下的隨機抽樣，并根據(jù)所建立的誤差數(shù)學(xué)模型，得到被測量誤差的概率分布情況。最后以某一實際待測情況為例，用蒙特卡洛法進行了相同概率分布下的隨機抽樣，并根據(jù)所建立的誤差數(shù)學(xué)模型，得到被測量誤差的概率分布情況。

1　傳感器測量原理與安裝標定

1．1基本結(jié)構(gòu)與測量原理

非接觸式平面激光測距傳感器由3個相同的激光測距傳感器、安裝基座和被測對象組成，如圖1所示。

其測量原理為3個相同的非接觸式激光測距傳感器測量相對于安裝基準的距離值，通過對所獲測量數(shù)據(jù)的組合算法，從而實現(xiàn)對物體平面運動兩個方向相對位移的測量及繞Z軸轉(zhuǎn)動角度的測量。非接觸式激光測距傳感器測量頻率較高，能滿足系統(tǒng)閉環(huán)控制的要求。

1．2傳感器安裝與標定

安裝基座上安裝3個相同的非接觸式激光測距傳感器，即在X軸方向安裝2個，在Y軸方向安裝一個，而安裝基座位于在模樣件相互垂直的相鄰側(cè)面上。

圖1　非接觸式并聯(lián)平面激光測距傳感器結(jié)構(gòu)

標定方法：模樣件的安裝面經(jīng)過精密加工，其垂直度誤差很小，可以忽略不計。設(shè)在X、Y軸方向上的2個激光測距傳感器的初始測量值分別為：y01、y02、x01，Y向傳感器測量理論值為x02。在該傳感器測量前，需要進行調(diào)整，直到滿足下式：

y01=y02，x01=x02

2　測量對象數(shù)學(xué)模型建立

設(shè)被測量為G，則其與輸入量之間的函數(shù)關(guān)系為：

G=f(x1，x2，l1，l2，l3，y1，|o1o2|，θ，t)

式中：l1、l2、l3為激光測距傳感器測量值；x1、x2、y1為與激光測距傳感器安裝位置相關(guān)的參數(shù)；|o1o2|、θ1為與模樣件質(zhì)心相關(guān)的參數(shù)；t為環(huán)境溫度。

由于環(huán)境溫度變化不大，且對被測量G的影響較小，忽略其影響[4]，因此，G=f(x1，x2，l1，l2，l3，y1，|o1o2|，θ1)。

測量X、Y軸方向傳感器安裝點位置坐標分別為：a1(x1，0)、a2(x2，0)、b(0，y1)，則被測物體上傳感器的坐標分別為：

u(x1，l1)，v(x2，l2)，w(l3，y1)

根據(jù)圖2中三點坐標及幾何關(guān)系o2o3⊥o2o4，得到直線o2o3、o2o4的方程分別為：

解上面的方程組得交點o2的x、y坐標為：

由交點o2的坐標、被測對象上任一點o1到o2的距離|o1o2|、直線o1o2與直線o2o3間的夾角θ1，可求得被測物體相對于安裝基座及激光測距傳感器X向、Y向位移和繞Z軸的轉(zhuǎn)動運動參數(shù)如下。

圖2　測量原理圖

3　誤差分析

在不考慮溫度變化、模樣件安裝面及測量基準板直線度等因素對被測量的影響下，用全微分法建立被測量誤差的數(shù)學(xué)模型[5]，并用蒙特卡洛法對測量誤差進行了定量分析。

3．1誤差數(shù)學(xué)模型建立

由于被測物體是三自由度的平面運動，故激光測距帶來的測量誤差由三部分組成，分別是：X、Y向位移誤差ΔX、ΔY、繞Z向的轉(zhuǎn)角誤差Δθ，誤差矩陣為：

3．2測量誤差蒙特卡洛法模擬

誤差的蒙特卡洛模擬，就是在計算機對具有不同分布的誤差隨機量進行抽樣，再按測量誤差的數(shù)學(xué)模型算出大量的抽樣數(shù)據(jù)，最終借助統(tǒng)計分析定量的獲取測量誤差的分布規(guī)律的一種方式[6]。

3．2．1偽隨機數(shù)生成

生成(0，1)之間均勻分布的偽隨機數(shù)的算法主要有：線性同余(MLCG)法、移位寄存器(或Tausworthe)法、Fibonacci法和混合法，這里使用線性同余法。

線性同余法由D．H．Lehmer提出，其基本形式是對任意整數(shù)X0，由以下公式確定：

X1=(aXi-1+c)modM

ξi=Xi/M，i=1，2，…

式中：乘子a為小于模M的正整數(shù)；c為非負整數(shù)，常選為零；ξ為偽隨機數(shù)。

當(dāng)c=0時，稱為乘同余法，其滿周期為：T=M-1。有研究表明，乘同余法優(yōu)于加同余法，而混合同余法與簡單的乘同余法相比并無顯著優(yōu)點[6]。目前，具有參數(shù)為M=231-1，以及a=16 807或630 360 016的乘同余方法得到廣泛應(yīng)用。

3．2．2偽隨機數(shù)的隨機性檢驗

偽隨機數(shù)序列的隨機性檢驗，就是按照從(0，1)上均勻分布隨機所應(yīng)具有的性質(zhì)和規(guī)律，來研究所產(chǎn)生的偽隨機數(shù)序列是否與這些性質(zhì)和規(guī)律相符合。

偽隨機數(shù)通過的檢驗越多，說明其越可靠。這里對所產(chǎn)生的偽隨機數(shù)用矩檢驗法進行參數(shù)檢驗，用χ2檢驗法檢驗偽隨機數(shù)的均勻性，用游程檢驗法檢驗偽隨機數(shù)的獨立性。

3．2．3隨機變量抽樣

隨機變量的抽樣是為了生成與原始誤差有相同概率分布和統(tǒng)計特征的隨機變量的抽樣值。由于G=f(x1，x2，l1，l2，l3，y1，|o1o2|，θ1)的輸入?yún)?shù)之間是獨立的，由于各隨機變量在有限區(qū)間都有上界，因此可以用簡單分布的舍選抽樣方法。具有概率密度f0(x)的隨機變量x∈(m，n)的抽樣公式為：

式中：η1、η2為(0，1)上的偽隨機數(shù)；fmax為f0(x)的最大值。

對f0(x)的舍選抽樣方法的流程圖如圖3所示。

圖3　抽樣流程圖

3．2．4算例誤差的數(shù)值計算

若某激光測距傳感器的最小采樣時間為：20 μs，位移分辨率為：1 μm，重復(fù)精度為：+2 μm，其安裝精度為：0．02 mm，則對于如下輸入?yún)?shù)：

u(520.000±0.01，342.019±0.002)(mm)

v(680.000±0.01，416.629±0.002)(mm)

w(379.324±0.002，600.000±0.01)(mm)

|o1o2|=(283.285±0.005)(mm)，θ1=45.00°±0.01°

用乘同余法生成(0，1)區(qū)間的偽隨機數(shù)，并對偽隨機數(shù)進行檢驗，偽隨機數(shù)的檢驗結(jié)果見表1，從表中可以看出所產(chǎn)生的10 000個偽隨機數(shù)滿足隨機性、均勻性、獨立性要求。

表1　偽隨機數(shù)檢驗結(jié)果

圖4　X向位移測量誤差分布

圖5　Y向位移測量誤差分布

圖6　Δθ的誤差分布

4　結(jié)束語

文中提出了一種基于并聯(lián)的新型非接觸式平面激光測距傳感器，通過組合算法，能夠測量出平面運動物體平移和轉(zhuǎn)動量，實現(xiàn)了物體平面運動的非接觸式測量。采用蒙特卡洛方法模擬出這種新型傳感器平移和轉(zhuǎn)動誤差的概率分布，得到X、Y向測量精度分別高于±0.056 mm、±0.029 5 mm，轉(zhuǎn)動精度高于5.6×10-3°，為這種傳感器的實際應(yīng)用提供了依據(jù)。

參考文獻：

[1]姜春英，房立金，徐志剛，等．一種平面運動位姿組合測量方法研究．機械工程學(xué)報，2006，17(4)：362-366．

[2]王志乾，謝幕君，喬彥峰．一種非接觸式測量方法的研究．光學(xué)技術(shù)，2003，7(4)：431-433．

[3]羅相杰，段純剛，胡雙啟，等．某軍用三坐標激光測量系統(tǒng)的安全性及精度控制方法研究．激光雜志，2006，27(3)：18-20．

[4]BAROUDI，DJEBAR，THIBERT．An instrumented structure to measure avalanche impact pressure：Error analysis from Monte Carlo simulations．EI．Cold Regions Science and Technology，2009，59(2)，242-250．

[5]費業(yè)泰．誤差理論與數(shù)據(jù)處理．北京：中國機械工業(yè)出版社，2010．

[6]徐鐘濟．蒙特卡羅方法．上海：上?？茖W(xué)技術(shù)出版社，1985-06-ISBN．

[7]Haima YANG，Ruiting LI，Qingzhong WEI，et al．The study of high accuracy capacitive displacement sensor used in non-contact precision displacement measurement．The Ninth International Conference on Electronic Measurement & Instruments，Bei jing，2009．