結(jié)構(gòu)參數(shù)約束下的磁場補償式雙層時柵角位移傳感器研究

楊繼森 周 潤 張?zhí)旌?盧 渝 吳 灼 張 迪 張 靜

重慶理工大學(xué)機械檢測技術(shù)與裝備教育部工程研究中心,重慶,400054

0 引言

近年來,我國制造業(yè)正迅速向高質(zhì)量制造轉(zhuǎn)型升級,高質(zhì)量制造對位移傳感器的測量精度提出了更高要求,研究高精度位移傳感器顯得愈發(fā)重要與緊迫[1-5]。時柵位移傳感器作為一種新型位移傳感器,以勻速運動的場為參考系,將空間位置測量轉(zhuǎn)化為相對時間測量,提供了一種新型位移傳感思路,提高了傳感器測量精度。傳感器以抗干擾能力強、精度高、抗油污、抗抖動等特點廣泛用于數(shù)控機床、工業(yè)機器人以及軍用領(lǐng)域等[6-8]。

場作為時柵位移傳感器測量過程中的重要部分,場的勻速性與均勻性直接影響了傳感測量精度,在研究進程中一直作為重點研究對象。如今印刷電路板(printed circuit board,PCB)工藝成熟,布線精度高,靈活度強,制造周期短,有助于優(yōu)化傳感器布線結(jié)構(gòu),增強場的勻速性,提高傳感器的測量精度[9-11]。時柵位移傳感器也從以前“開槽繞線式”結(jié)構(gòu)轉(zhuǎn)變?yōu)橛蒔CB工藝制造的平面線圈陣列結(jié)構(gòu)。文獻[12]中提出了反向串聯(lián)“回”字形陣列結(jié)構(gòu),對傳感器結(jié)構(gòu)上進行了一定的優(yōu)化,消除了齒槽效應(yīng),提高了行波磁場的運動勻速性,減小了測量結(jié)果的四次諧波誤差;文獻[13]在文獻[12]的結(jié)構(gòu)基礎(chǔ)上采用非線性規(guī)劃遺傳算法對其結(jié)構(gòu)參數(shù)進行優(yōu)化,使得感應(yīng)信號諧波畸變率減小了90%,測量精度得到了提高;文獻[14]通過分析測量模型得到感應(yīng)信號與平板間隔無關(guān)的規(guī)律,通過一系列結(jié)構(gòu)變換將雙排差動平板結(jié)構(gòu)優(yōu)化為單排差動平板結(jié)構(gòu),消除串?dāng)_與安裝所帶來的影響,提高傳感器的信號穩(wěn)定性、抗干擾能力,在200 mm的測量范圍內(nèi)測量誤差值從300 nm減小至150 nm。筆者所在課題組在優(yōu)化時柵位移傳感器工作中,通過分析磁場式時柵的時變磁場產(chǎn)生原理,提出采用雙層互補式線圈的時柵傳感器設(shè)計方案,有效抑制了時變磁場“端部效應(yīng)”,提高了測量精度[15]。后續(xù)的研究中發(fā)現(xiàn),仿真過程均在理想的情況下進行,忽略了上下層線圈間距,進一步研究發(fā)現(xiàn),線圈間距會引起上下層線圈在感應(yīng)線圈平面產(chǎn)生的磁場不對稱,使得理論上均勻?qū)ΨQ的磁場在感應(yīng)平面內(nèi)分布不對稱,該現(xiàn)象將導(dǎo)致傳感器感應(yīng)電動勢在互補位置的幅值不等,最終影響測量精度。針對上述問題,本文提出了一種新優(yōu)化思路,采用磁場補償驅(qū)動方案提高磁場對稱性和均勻性,彌補結(jié)構(gòu)上的不足,以達到提高傳感器測量精度的目的。

1 測量原理及數(shù)學(xué)模型搭建

1.1 雙層時柵角位移傳感器測量原理

雙層時柵角位移傳感器總體結(jié)構(gòu)如圖1所示,傳感器由激勵部分與感應(yīng)部分組成,兩部分均采用雙層互補式結(jié)構(gòu)。激勵部分由節(jié)距為w的矩形線圈組成,共有6個扇區(qū),相鄰扇區(qū)相差四分之一個節(jié)距,以實現(xiàn)傳感器激勵線圈在空間上的正交,感應(yīng)部分由節(jié)距為w的正弦型結(jié)構(gòu)組成。

圖1 雙層時柵位移傳感器結(jié)構(gòu)示意圖Fig.1 Schematic diagram of double layer time gate displacement sensor

相鄰扇區(qū)激勵線圈分別通入相位差為π/4的兩路激勵信號,以實現(xiàn)驅(qū)動信號在時間上的正交。傳感器的測量原理如圖2所示。在感應(yīng)線圈所在平面上,根據(jù)畢奧-薩法爾定律對通電導(dǎo)線產(chǎn)生的磁場進行疊加形成時變磁場,圖2a為感應(yīng)線圈平面時變磁場分布與感應(yīng)線圈內(nèi)部磁通量變化示意圖,其中圓點表示磁場方向為垂直平面向外,叉號表示磁場方向為垂直平面向內(nèi),Hc為感應(yīng)線圈高度。由圖2a可知,當(dāng)感應(yīng)線圈沿著X軸正方向移動θ時,此時感應(yīng)線圈內(nèi)部有正負磁場強度同時作用,負磁場所穿過的面積抵消了部分正磁場所穿過的面積,正好對應(yīng)正弦函數(shù)陰影部分半個周期的積分形式。根據(jù)法拉第電磁感應(yīng)定律,感應(yīng)線圈上將感應(yīng)生成一路感應(yīng)電動勢,感應(yīng)電動勢與磁通量和位移量的關(guān)系為

(1)

(a)時變磁場平面分布與感應(yīng)線圈內(nèi)部磁通量示意圖

(b)傳感器感應(yīng)信號計時原理圖2 傳感器測量原理示意圖Fig.2 Schematic diagram of sensor measurement principle

式中,ΔΦ為穿過感應(yīng)線圈的磁通量的變化量;Δtx為位移前后時間差;K1為感應(yīng)電動勢的系數(shù);ω為角頻率;t為時間;w為感應(yīng)線圈節(jié)距;θ為感應(yīng)線圈與激勵線圈的相對角位移量。

傳感器激勵信號在時間上正交,激勵線圈結(jié)構(gòu)在空間上正交,由此得出另外一路感應(yīng)電動勢表達式為

(2)

由于傳感器感應(yīng)線圈相互串聯(lián),因此傳感器輸出的行波信號是兩路感應(yīng)電動勢之和,則行波信號表達式為

(3)

由式(3)可知,位移量的表達式為行波信號的相位值,因此測量行波信號與參考信號的相位差能實現(xiàn)位移量的解算。行波信號經(jīng)過放大、濾波、整形處理之后,形成一路帶有位置信息的方波信號,方波信號在處理器中與參考信號作數(shù)字信號處理,采用高頻時鐘插補技術(shù)測量該方波信號周期T以及與參考信號的時間差Δt,Δt與T可得相位差,從而實現(xiàn)對位移量θ的解算。計數(shù)原理如圖2b所示。θ與Δt和T的關(guān)系可表示為

(4)

式中,M為移動節(jié)距數(shù)。

1.2 雙層時柵角位移傳感器數(shù)學(xué)模型

在1.1節(jié)測量原理介紹中,時變磁場均在理想的條件下進行計算,忽略了上下層激勵線圈間距。本文在考慮線圈間距的基礎(chǔ)上,重新建立雙層正交型結(jié)構(gòu)的數(shù)學(xué)模型,對時變磁場作進一步分析。由于傳感器結(jié)構(gòu)精細,在較大對極數(shù)的情況下,外徑線圈弧度與內(nèi)徑線圈弧度近似于平行線,因此根據(jù)圖1中激勵線圈建立磁場仿真模型,如圖3所示。線圈沿Y方向長度為2b,沿X軸方向長度為2a,兩層激勵線圈相距為c,在感應(yīng)線圈平面內(nèi)將任意場點設(shè)為P(m,n,k),場點P到導(dǎo)線B′C′的距離為H,P′為場點投影到OXY平面的坐標點。原有驅(qū)動方案是將上下層激勵線圈進行串聯(lián)并通入幅值為I的電流,P點的磁場強度由6根通電導(dǎo)線BC、CD、DA、D′A′、A′B′、B′C′共同疊加而成。Q為通電導(dǎo)線B′C′中的一點,Q點電流元在P點處產(chǎn)生的元磁場(元磁場方向為垂直于平面QPC′)根據(jù)畢奧-薩法爾定律可矢量表達為

(5)

式中,B為磁感應(yīng)強度;μ0為真空磁導(dǎo)率;Idl為導(dǎo)線上Q點的線電流元;dl為電流的微小線元素;QP為Q點指向場點P的向量。

圖3 激勵線圈磁場仿真模型Fig.3 Magnetic field simulation model of excitation coil

根據(jù)圖3設(shè)點Q的坐標為(a,y,0),y為Y軸上的位移,則dB可用向量的方式表示為

dB=(dBX,dBY,dBZ) =

(6)

式中,dBX、dBY、dBZ分別為沿坐標軸三個方向的分量。

由于感應(yīng)線圈平面與激勵線圈平面相平行,因此只需要考慮感應(yīng)線圈平面的法線方向(Z軸)上的磁場dBZ。

根據(jù)三角公式可得

(7)

(8)

式中,θ1為線段B′P與線段B′C′的夾角;θ2為線段C′P與線段B′C′的夾角。

采用積分法聯(lián)立式(7)和式(8)可求得通電導(dǎo)線B′C′在場點P處沿Z軸方向上產(chǎn)生的磁場強度BB′C′Z為

(9)

其中,x、z分別為X軸和Z軸上的位移。同理可得其余5條導(dǎo)線在場點P處沿Z軸方向上產(chǎn)生的磁場強度BA′B′Z、BD′A′Z、BBCZ、BCDZ、BDAZ。

由磁場疊加定理可得場點P沿Z軸方向的磁場強度為

BZ=BB′C′Z+BA′B′Z+BD′A′Z+BBCZ+BCDZ+BDAZ

(10)

2 磁場補償驅(qū)動方案設(shè)計

2.1 磁場仿真

根據(jù)式(10)所建立的數(shù)學(xué)模型不能直觀地反映時變磁場的不對稱性,因此在數(shù)學(xué)模型的基礎(chǔ)上結(jié)合MATLAB軟件進行磁場仿真,仿真參數(shù)如表1所示,仿真結(jié)果如圖4所示。

表1 仿真參數(shù)設(shè)置

(a)磁場強度三維圖

(b)沿Y軸方向平面磁場強度分布圖圖4 感應(yīng)線圈平面磁場分布圖Fig.4 Planar magnetic field distribution diagram of induction coil

由圖4可知,該平面中心區(qū)域的磁場強度為5.3×10-8T,在Y方向上越靠近線圈端部磁場分布愈發(fā)不對稱,若以中心區(qū)域的磁場強度為平均值時,Y方向上位于兩邊端部位置的磁場強度相對誤差最大為10%,因此考慮線圈間距時激勵線圈所產(chǎn)生的時變磁場分布不對稱。

(a)磁場強度三維圖(中部四根通電導(dǎo)線)

(b)沿Y軸方向平面磁場強度分布圖(中部四根通電導(dǎo)線)

(c)磁場強度三維圖(端部兩根通電導(dǎo)線)

(d)沿Y軸方向平面磁場強度分布圖(端部兩根通電導(dǎo)線)圖5 分層分析時感應(yīng)線圈平面磁場分布圖Fig.5 Planar magnetic field distribution diagram of induction coil during stratified analysis

為進一步分析導(dǎo)致該現(xiàn)象的因素,對激勵線圈進行了分層討論,上層激勵線圈與下層激勵線圈各自單獨作用時的磁場分布如圖5所示。上層線圈lD′A′和lB′C′與下層線圈lBC和lDA產(chǎn)生的磁場分布分別如圖5a和圖5b所示,從圖中可知中部四根通電導(dǎo)線在感應(yīng)線圈平面產(chǎn)生的磁場強度約為5.3×10-8T,四根導(dǎo)線所產(chǎn)生的磁場關(guān)于坐標軸對稱,由此可知這四根導(dǎo)線并不是導(dǎo)致磁場不對稱的原因。圖5c、圖5d所示分別為導(dǎo)線lA′B′和lCD通入相同激勵電流時在感應(yīng)平面所產(chǎn)生的磁場強度,從圖中可以看出,兩根通電導(dǎo)線在感應(yīng)平面產(chǎn)生的磁場強度左右不對稱,這是由于lA′B′和lCD位于端部遠離中心區(qū)域,在線圈中部影響不大,只對接近端部的磁場影響較大,且上層端部線段lA′B′產(chǎn)生的磁場強度較大,約為1.7×10-8T,因此可以得出影響磁場不對稱的因素為不在同一平面的兩根通電導(dǎo)線lA′B′和lCD。

2.2 磁場補償式驅(qū)動方案與結(jié)構(gòu)參數(shù)的關(guān)系

根據(jù)磁場強度表達式,感應(yīng)線圈所在平面的磁場強度與激勵線圈高、激勵線圈寬、空氣間隙、上下層線圈間距以及驅(qū)動電流有關(guān)。傳感器結(jié)構(gòu)與各個參數(shù)相對應(yīng),其中線圈高、線圈寬和上下層線圈間距與傳感器線圈結(jié)構(gòu)有關(guān)。當(dāng)傳感器結(jié)構(gòu)確定時,感應(yīng)平面上的磁場強度僅與電流相關(guān)。

通過對6條通電導(dǎo)線的磁場分析可知,遠離感應(yīng)線圈平面的通電導(dǎo)線lCD產(chǎn)生的磁場強度弱于lA′B′產(chǎn)生的磁場強度,最終導(dǎo)致時變磁場不對稱。為了解決時變磁場不對稱的問題,通過增大下層激勵線圈電流使得下層感應(yīng)線圈磁場強度增大,以抵消通電導(dǎo)線與lCD所產(chǎn)生的不對稱磁場,并以此作為磁場補償式驅(qū)動方案的核心思想。當(dāng)傳感器的結(jié)構(gòu)確定時各個參數(shù)均已確定,則感應(yīng)平面場點磁場強度表達式為關(guān)于x與y的二元二次表達式,可以根據(jù)關(guān)于二元二次的極值求解方法求得中部磁場強度的最大值及其對應(yīng)坐標A(x1,y1),同理可得整個結(jié)構(gòu)在端部所產(chǎn)生的磁場最大值及其對應(yīng)坐標B(x2,y2),將A和B點的磁場強度分別作為中間部分磁場強度和端部的磁場強度的參考點,增大下層線圈電流使得A、B兩點磁場強度相等,以達到抵消端部磁場的不對稱性的目的。

將下層激勵線圈電流設(shè)為I2,則生成的新的磁場強度BZ1為

BZ1=BB′C′Z+BA′B′Z+BD′A′Z+

(BBCZ+BCDZ+BDAZ)I2/I

(11)

根據(jù)磁場補償驅(qū)動方案的核心思想可得

BZ1(x1,y1)=BZ1(x2,y2)

(12)

結(jié)合式(11)和式(12)即可得電流I與I2的關(guān)系表達式為

I2=I(BB′C′Z(x1,y1)+BA′B′Z(x1,y1)+

BD′A′Z(x1,y1)-BB′C′Z(x2,y2)+BA′B′Z(x2,y2)+

BD′A′Z(x2,y2))/(BBCZ(x2,y2)+BCDZ(x2,y2)+

BDAZ(x2,y2)-BBCZ(x1,y1)+BCDZ(x1,y1)+

BDAZ(x1,y1))=If(a,b,c,z)

(13)

其中,f(a,b,c,z)表示關(guān)于傳感器結(jié)構(gòu)參數(shù)的函數(shù)。將表1中的結(jié)構(gòu)參數(shù)代入式(13),可求得I2為118 mA。

為了驗證磁場補償思路的可行性,將得出的電流幅值結(jié)合仿真模型再次進行磁場仿真,設(shè)置的參數(shù)如表2所示,磁場分布如圖6所示,由圖可知,當(dāng)下層激勵線圈電流為118 mA時磁場強度增大到5.8×10-8T,磁場強度相對誤差減小為0.3%。對比圖4可知,與原有的驅(qū)動方案相比,利用這種磁場補償驅(qū)動方案得到的磁場分布更加對稱,可以抵消端部導(dǎo)線所帶來的磁場不對稱性。

表2 磁場仿真參數(shù)設(shè)置

(a)補償后磁場強度三維圖

(b)補償后沿Y軸方向平面磁場強度分布圖圖6 補償后感應(yīng)線圈平面磁場分布圖Fig.6 Magnetic field distribution diagram of induction coil after compensation

3 有限元仿真與誤差分析

3.1 有限元對比仿真

根據(jù)第2節(jié)的磁場仿真結(jié)果,實際情況下傳感器節(jié)距內(nèi)時變磁場存在不對稱,但在測量原理中磁場強度均是采用理想狀態(tài)下的均勻磁場進行計算的。根據(jù)式(1)所示的感應(yīng)電動勢的表達式,當(dāng)傳感器感應(yīng)線圈處于激勵線圈節(jié)距內(nèi)的互補位置時,面積變化量不變,但傳感器激勵線圈節(jié)距內(nèi)時變磁場并不是呈理想的對稱分布,從而造成感應(yīng)信號在互補位置時幅值不相等。實際情況下兩路感應(yīng)電動勢為

(14)

(15)

式中,kc為互補位置的幅值差值。

由于行波信號是由兩路感應(yīng)電動勢合成的,則合成后的行波信號可表示為

(16)

式中,K為行波信號的幅值;φ1、φ2為行波信號的相位值。

根據(jù)三角函數(shù)變化可得

(17)

根據(jù)式(16)、式(17)得出,不對稱時變磁場影響了相位差,進一步導(dǎo)致解算出來的位移量出現(xiàn)偏差,因此為了驗證磁場補償驅(qū)動方案的可行性,采用分析軟件Maxwell對傳感器的結(jié)構(gòu)模型進行有限元仿真,分析該方案對感應(yīng)電動勢幅值以及測量誤差的影響,仿真結(jié)構(gòu)模型如圖1所示。本文的仿真采用兩種不同驅(qū)動方案進行對比,仿真參數(shù)如表3所示,仿真結(jié)果如圖7所示。

表3 有限元仿真參數(shù)

(a)正弦信號激勵感應(yīng)電動勢

(b)余弦信號激勵感應(yīng)電動勢

(c)行波信號圖圖7 感應(yīng)電動勢與行波信號圖Fig.7 Induced electromotive force and traveling wave signal diagram

3.2 誤差分析

仿真產(chǎn)生兩路感應(yīng)電動勢,每一路有50個采樣點。根據(jù)感應(yīng)電動勢的表達式,感應(yīng)電動勢的幅值與空間位移量和節(jié)距有關(guān),一個節(jié)距內(nèi)轉(zhuǎn)子在互補位置時幅值相等符號相反。取25對互補位置的感應(yīng)電動勢幅值進行取模作差,兩種不同驅(qū)動方式下互補位置的感應(yīng)電動勢幅值相對誤差如圖8所示。

圖8 互補位置感應(yīng)電動勢幅值誤差Fig.8 Amplitude error of complementary position induced electromotive force

由圖8可以看出,不同驅(qū)動方案下感應(yīng)電動勢幅值差值曲線的變化趨勢相似,原有驅(qū)動方案激勵下,幅值的相對誤差整體偏大,最大值出現(xiàn)在1.3°,磁場補償驅(qū)動方案激勵時,互補位置的感應(yīng)電動勢幅值相對誤差明顯減小,說明磁場補償驅(qū)動方案可以減小轉(zhuǎn)子在互補位置上感應(yīng)電動勢的幅值相對誤差。

通過將兩路感應(yīng)電動勢相加可得帶有位置信息的行波信號,如圖7c所示。根據(jù)圖2中的求解原理可得行波信號的相位差與傳感器的位移量。故可將0°位置處行波信號的相位作為初始相位值,以此計算每個角度行波信號的相位與初始相位的相位差,進而求解出各個角度相對于0°位置的相對位移。將測量位移量與真實位移量作比較可得每個位置的角度誤差。兩種不同驅(qū)動方案激勵時傳感器模型的角度誤差曲線見圖9。

圖9 兩種不同驅(qū)動方案下傳感器誤差曲線Fig.9 Sensor error curves under two different driving schemes

由圖9可知,兩種驅(qū)動方案驅(qū)動時誤差曲線呈相同的周期性變化,當(dāng)磁場補償驅(qū)動方案作用時,轉(zhuǎn)子在1.3°～2.4°與2.9°～4°時測量誤差明顯減小。從仿真數(shù)據(jù)分析中可以得到,沒有補償電流時,測量誤差范圍為(-31.6″,26.9″),當(dāng)有補償電流作用時測量誤差范圍為(-22.4″,14.2″)。由此可以得出,當(dāng)磁場補償驅(qū)動方案激勵時測量誤差峰峰值從之前的58.5″變?yōu)?6.6″,誤差峰峰值減小了37.4%。

4 實驗驗證

4.1 實驗臺的搭建

根據(jù)圖1中的傳感器結(jié)構(gòu)模型,采用PCB工藝加工傳感器樣機。安裝時柵角位移傳感器并搭建對比實驗平臺,為測量時柵角位移傳感器精度,采用更高精度等級的海德漢光柵進行誤差標定,光柵型號為RON-866,傳感器整周刻有36 000個柵線,系統(tǒng)精度為±1″;采用電子細分盒處理光柵輸出信號,可以將光柵分辨率提高100倍。安裝時將角位移編碼器與時柵角位移傳感器同軸安裝,從而使兩者同步運動,信號處理電路同步采集兩個傳感器數(shù)據(jù),通過串口上傳到上位機進行誤差標定,實驗現(xiàn)場如圖10所示。

圖10 雙層時柵角位移傳感器對比實驗臺Fig.10 Double layer time grating angular displacement sensor contrast experimental platform

4.2 對比實驗

雙層時柵角位移傳感器驅(qū)動頻率為10 kHz,根據(jù)磁場補償驅(qū)動方案可得兩層驅(qū)動電壓分別為4.25 V和5.02 V,本次實驗將兩層驅(qū)動電壓設(shè)置為4.3 V和5.0 V進行驅(qū)動,與原有驅(qū)動方案中上下層驅(qū)動電壓均為4.3 V展開對比實驗。設(shè)置一個節(jié)距內(nèi)采樣點個數(shù)為50,得到不同驅(qū)動方案下的誤差曲線見圖11。

圖11 節(jié)距內(nèi)實驗誤差曲線對比Fig.11 Comparison of error curves in polar experiment

實驗誤差曲線表明,節(jié)距內(nèi)測量出的誤差從(-25.5″,18.3″)減小到(-20.8″,12.6″),磁場補償驅(qū)動方案提高了傳感器測量精度,實驗結(jié)果與仿真結(jié)果基本一致。

5 結(jié)論

(1)本次研究通過對雙層時柵角位移傳感器結(jié)構(gòu)建模與時變磁場的仿真,發(fā)現(xiàn)上層激勵線圈與下層激勵線圈的間隙會對時變磁場產(chǎn)生影響,從而導(dǎo)致在互補位置時感應(yīng)電動勢幅值不等。由此引出一種磁場補償驅(qū)動方案,該方案將上下兩層激勵線圈通入不同幅值電壓以消除間隙對時變磁場的影響,通過磁場補償思路與傳感器的結(jié)構(gòu)參數(shù)確定了與之對應(yīng)的上下層激勵線圈驅(qū)動電壓關(guān)系。同時經(jīng)過數(shù)學(xué)分析、建模仿真、實驗對比,驗證了磁場補償驅(qū)動方案的有效性。該驅(qū)動方案提高了行波磁場的對稱性,經(jīng)對比研究后,節(jié)距內(nèi)測量誤差從(-25.5″,18.3″)減小到(-20.8″,12.6″)。

(2)本次研究更貼近實際情況下的傳感器研究,通過優(yōu)化驅(qū)動方案來彌補結(jié)構(gòu)上的不足,并且適用于不同結(jié)構(gòu)參數(shù)的雙層時柵位移傳感器,為高精度傳感器優(yōu)化設(shè)計提供了新思路。