一種實(shí)現(xiàn)彈體感應(yīng)磁場補(bǔ)償?shù)姆抡娣椒ㄑ芯?/h1>
2014-02-28 05:10:30王麗黎尹文超
網(wǎng)絡(luò)安全與數(shù)據(jù)管理 2014年2期
王麗黎，尹文超
(西安理工大學(xué) 自動(dòng)化與信息工程學(xué)院，陜西 西安 710048)
在地磁匹配導(dǎo)航系統(tǒng)中，載體自身干擾磁場的分析與補(bǔ)償是一個(gè)關(guān)鍵問題。根據(jù)飛行器上的磁干擾的類型可以建立磁干擾的數(shù)學(xué)模型，模型中含有一系列的系數(shù)，在飛行器上的儀器都固定的條件下，這些系數(shù)是固定的。測得需要的數(shù)據(jù)就可以推導(dǎo)出磁補(bǔ)償系數(shù)，利用這些系數(shù)可以求出飛行器周圍引起的干擾磁場大小，從而在地磁場未知的情況下得出真實(shí)的地磁場值。
目前國內(nèi)在磁場補(bǔ)償方面的研究主要集中在概略分析[1]和算法的研究[2-4]，而針對基于電磁場仿真模型和數(shù)值計(jì)算方面的研究卻很少，少數(shù)單位僅僅研究了縮比模型周圍的磁場分布[5],還沒有基于真實(shí)模型的仿真補(bǔ)償研究。針對目前在磁補(bǔ)償方面還沒有電磁仿真研究的情況，本文使用Ansoft Maxwell有限元軟件，仿真出物體在不同磁導(dǎo)率的情況下以及不同磁力計(jì)放置位置時(shí)的測量總場值，依據(jù)這些測量值可以求出彈體外表面任意點(diǎn)及材料在任意磁導(dǎo)率下的補(bǔ)償系數(shù)，并利用這些系數(shù)估計(jì)出背景地磁場。
1 導(dǎo)彈磁場補(bǔ)償模型分析與軟件建模
1.1 彈體磁場模型求解分析
根據(jù)飛行器磁干擾的Tolles-Lawson數(shù)學(xué)模型及表達(dá)式[6]作為分析模型，本模型以仿真感應(yīng)磁場為主，暫不考慮恒定場和渦流場，建立彈體坐標(biāo)系如圖1所示。L軸與彈體縱軸平行，首部為正；T軸與L軸垂直，V軸與TOL平面正交，向下為正。α、β、γ分別為 He與3個(gè)坐標(biāo)軸之間的夾角， 其方向余弦分別為 cosα、cosβ、cosγ。He為地磁場矢量，感應(yīng)出來的干擾場投影到地磁場方向大小為 Hid，化簡后為表達(dá)式(1)：
其中感應(yīng)場有 6個(gè)補(bǔ)償系數(shù) (TT-LL、VV-LL、LL、LT+TL、VT+TV、VL+LV),磁力計(jì)測量總場 Hz減去地磁場 He可得干擾場Hid，再利用地磁場矢量與彈體坐標(biāo)系3軸夾角和干擾場大小來求出補(bǔ)償系數(shù)。
圖1 彈體坐標(biāo)系示意圖
1.2 模型的建立
在Maxwell中建立Tolles-Lawson模型，彈體總長度為7 m，前半部分彈頭為圓錐型總長度3 m，后半部分彈體包括尾翼總長度4 m，半徑為30 cm，大小和形狀接近真實(shí)模型，使用永磁體加上適當(dāng)?shù)某C頑力來形成一定方向的均勻磁場，如圖2所示。彈體在均勻磁場會(huì)受到磁化并在周圍產(chǎn)生感應(yīng)磁場。由于軟件所仿真出的磁場值結(jié)果輸出都以Global坐標(biāo)系為準(zhǔn)，在軟件仿真中通過改變地磁場對物體的作用方向來模擬彈體在空間內(nèi)與地磁場形成任意角度，將軟件輸出值作為模擬磁力計(jì)的測量值。模型的材料及屬性參數(shù)如表1所示。
圖2 仿真模型
表1 模型的材料屬性
2 磁場補(bǔ)償技術(shù)的仿真與計(jì)算
2.1 仿真原理和計(jì)算步驟
磁干擾表達(dá)式(1)中有6個(gè)未知補(bǔ)償系數(shù)，求出方程的解最少需要6個(gè)方程，將式(1)寫成矩陣形式(2),左列為彈體在地磁場n種姿態(tài)時(shí)的干擾場大小，右側(cè)為對應(yīng)的姿態(tài)矩陣和感應(yīng)系數(shù)的乘積，這里使用最小二乘法求解。
仿真時(shí)令初始磁場方向垂直XOY面向下,再令磁場方向相對Y軸逆時(shí)針旋轉(zhuǎn)形成參數(shù)為α的角度，之后相對Z軸逆時(shí)針旋轉(zhuǎn)形成參數(shù)為β角度,如圖3所示。
圖3 模型的側(cè)面和縱截面
本文研究磁力計(jì)安裝在彈體表面的情況，通過軟件來實(shí)現(xiàn)多點(diǎn)采樣。這樣可以找出磁力計(jì)安裝在某個(gè)位置使得補(bǔ)償效果達(dá)到最好，本例取距離彈體表面5 cm沿彈體縱軸線的一條路經(jīng)作為磁力計(jì)放置一系列的測量位置，目的是要算出這條路徑上所有采樣點(diǎn)的補(bǔ)償系數(shù)，最后再求出所有位置的補(bǔ)償誤差，確定補(bǔ)償效果最佳的磁力計(jì)放置位置。
在軟件求解前需要進(jìn)行參數(shù)設(shè)置,如表2所示。
表2 求解參數(shù)設(shè)置
感應(yīng)磁場是在靜態(tài)的恒定磁場作用下產(chǎn)生的，因此選用靜磁場求解器。零邊界條件保證了外部沒有磁場并且內(nèi)部磁場方向單一且垂直于邊界面。參數(shù)化α1和β1的目的是形成對物體不同作用方向的磁場，每個(gè)組合對應(yīng)一個(gè)彈體在地磁場中的一個(gè)姿態(tài)，由于彈體模型分別關(guān)于XOY面和XOZ面對稱并且地磁場反向作用時(shí)僅僅改變最后磁場值的符號，所以這里參數(shù)化0～90代表地磁場在一個(gè)象限內(nèi)作用物體，通過maxwell的參數(shù)設(shè)置step可以包含任意磁場對彈體的作用方向，以此來模擬飛行器在空間中的姿態(tài)。參數(shù)化 α2和β2形成的磁場作為驗(yàn)證，利用前部分求出的系數(shù)和后部分求出的總場來估計(jì)地磁場大小，最后將估計(jì)出的磁場值與形成的背景地磁場相比較來評價(jià)補(bǔ)償結(jié)果。將上述參數(shù)設(shè)置好后即可利用軟件進(jìn)行數(shù)值計(jì)算。
2.2 數(shù)據(jù)的后處理
導(dǎo)出測磁補(bǔ)償路徑上所有點(diǎn)在相對磁導(dǎo)率為1～5時(shí)彈體處于地磁場中所有姿態(tài)的地磁場三分量BXBYBZ和模值MagB，在這條路徑上采樣700個(gè)點(diǎn)作為磁探頭的測量點(diǎn)，磁場大小為50 266 nt，同時(shí)計(jì)算出式(2)中的cosX、cosY、cosZ,表達(dá)式左邊的干擾場大小等于該磁探頭位置的測量值MagB減去地磁場大小得出。利用角度 α1、β1、step=10時(shí)參數(shù)化的結(jié)果求出所有位置的 6個(gè)補(bǔ)償系數(shù)，再利用角度 α2、β2、step=36時(shí)參數(shù)化的結(jié)果求出所有位置的補(bǔ)償系數(shù)來估計(jì)背景地磁場，并與真值He對比得出每個(gè)位置補(bǔ)償后的剩余標(biāo)準(zhǔn)差，最后對補(bǔ)償結(jié)果進(jìn)行分析。
3 對仿真方法可靠性的驗(yàn)證
下面建立一個(gè)簡單的靜態(tài)磁偶極子模型對軟件的計(jì)算精度進(jìn)行驗(yàn)證和分析。之后,以一個(gè)簡單實(shí)驗(yàn)來驗(yàn)證利用本文的補(bǔ)償方法還原真實(shí)地磁場。
3.1 驗(yàn)證精度
Ansoft公司的Maxwell 3D是一個(gè)功能強(qiáng)大、結(jié)果精確、易于使用的三維電磁場有限元分析軟件。
根據(jù)電磁場理論，磁偶極子模型[7]在側(cè)向空間中形成的磁感應(yīng)強(qiáng)度大小為式(3)，軸線方向形成的磁感應(yīng)強(qiáng)度大小為式(4)：
根據(jù)這一理論在磁場中建立磁偶極子模型，其磁場分布如圖4所示。將仿真出來的有限元解與理論值比較，如表3所示。有限元解與理論值的誤差百分比很低，說明數(shù)值計(jì)算精度較高，有限元解較接近實(shí)際理論值。
3.2 實(shí)驗(yàn)分析
實(shí)驗(yàn)在一個(gè)可以形成單一方向均勻磁場的磁屏蔽箱內(nèi)進(jìn)行。
實(shí)驗(yàn)時(shí)保證磁力計(jì)與物體相對位置不變，通過增加亥姆霍茲線圈的電流控制最終形成了51 807 nt的電磁場，磁場沿磁屏蔽箱的軸線方向，通過不斷改變物體在磁場內(nèi)的位置狀態(tài)測得每個(gè)狀態(tài)對應(yīng)的磁場三分量和總場強(qiáng)度。這個(gè)實(shí)驗(yàn)在一個(gè)平面內(nèi)進(jìn)行一些角度的旋轉(zhuǎn)，測量值、補(bǔ)償值和真實(shí)值的結(jié)果如圖5所示。
圖5 測量和補(bǔ)償結(jié)果
圖4 磁偶極子模型在空間中產(chǎn)生的磁感應(yīng)強(qiáng)度
表3 有限元解與理論解比較
圖5顯示了物體在磁場中36個(gè)狀態(tài)的補(bǔ)償結(jié)果，平均誤差比較小，直觀地反映了本文的補(bǔ)償方法是可行的。
該實(shí)驗(yàn)說明，通過不斷改變被補(bǔ)償物體在磁場中的不同狀態(tài)并利用最小二乘法可以近似地估計(jì)外部地磁場的大小。
4 仿真結(jié)果及其分析
在Maxwell軟件中仿真生成的背景地磁場He的大小為50 266 nt，在補(bǔ)償路徑上采樣700個(gè)測量點(diǎn)作為測磁點(diǎn)，同時(shí)計(jì)算了物體在磁導(dǎo)率從2～5每個(gè)點(diǎn)的補(bǔ)償結(jié)果，對結(jié)果的分析評價(jià)如圖6所示，圖表的縱坐標(biāo)為對彈體表面點(diǎn)在參數(shù)化(α2和β2)姿態(tài)時(shí)的補(bǔ)償剩余標(biāo)差。
圖6 某路徑在同一地磁場作用下不同磁導(dǎo)率下的補(bǔ)償結(jié)果
通過對仿真數(shù)據(jù)的處理，得出了沿著彈體表面某路徑在同一地磁場作用下不同磁導(dǎo)率的補(bǔ)償結(jié)果，同時(shí)得出了材料在相對磁導(dǎo)率在2～5時(shí)對多種飛行姿態(tài)補(bǔ)償后的平均誤差。其一，磁導(dǎo)率越大，補(bǔ)償后的誤差相對偏大，但并不是所有的補(bǔ)償位置的補(bǔ)償誤差都隨著磁導(dǎo)率的增大而增加；其二，基于本模型，彈體首部和彈體表面的補(bǔ)償誤差偏小，基本小于100 nT，而尾部的補(bǔ)償誤差較大。
本文通過電磁場仿真軟件Ansoft maxwell仿真飛行載體空間中的感應(yīng)磁場分布，并使用最小二乘法(LSE)進(jìn)行參數(shù)估計(jì)，求出彈體表面空間某曲線上所有采樣點(diǎn)的補(bǔ)償系數(shù)并進(jìn)行分析。本文提供了一種利用電磁場仿真軟件建立Tolles-Lawson模型并補(bǔ)償?shù)姆椒?，可以在不進(jìn)行試驗(yàn)的情況下仿真研究感應(yīng)磁場分布，實(shí)現(xiàn)了磁場補(bǔ)償，并對多種姿態(tài)的補(bǔ)償結(jié)果進(jìn)行了綜合評價(jià)。通過軟件參數(shù)化設(shè)置，對仿真數(shù)據(jù)進(jìn)行處理之后可以找到不同磁導(dǎo)率下模型的最佳補(bǔ)償點(diǎn)。本仿真方法可以對任意形狀模型在任意磁導(dǎo)率的情況下找出最佳補(bǔ)償位置，以及磁力計(jì)在同一安裝位置下找出補(bǔ)償效果最好的材料相對磁導(dǎo)率，在工程中可為材料磁導(dǎo)率的選取和磁力計(jì)的安裝位置提供依據(jù)。
[1]趙雙求.UAV航磁系統(tǒng)及補(bǔ)償研究[D].長沙：中南大學(xué)，2011.
[2]張堅(jiān),林春生,羅青,等.基于預(yù)測殘差平方和的飛機(jī)干擾磁場模型求解[J].探測與控制學(xué)報(bào)，2010,32(5):74-78.
[3]龐學(xué)亮,林春生,張寧.飛機(jī)磁場模型系數(shù)的截?cái)嗥娈愔捣纸夥ü烙?jì)[J].探測與控制學(xué)報(bào),2009,31(5):48-51.
[4]伍東凌,陳正想,王秀.基于遺傳算法的磁干擾補(bǔ)償方法[J].探測與控制學(xué)報(bào)，2012,34(6):16-20.
[5]李永慧，張曉明，劉俊,等.某常規(guī)彈載體感應(yīng)磁場分布的仿真技術(shù)研究[J].彈箭與制導(dǎo)學(xué)報(bào)，2012,32(4):181-184.
[6]LELIAK P.Identification and evaluation of magnetic field sources of magnetic airborne detector equipped aircraft[J].IRE Transactions on Airspace and Navigational Electronics，1961(8)：95-105.
[7]任來平,趙俊生,侯世喜.磁偶極子磁場空間分布模式[J].海洋測繪,2002,22(2):18-21.
[8]趙博，張洪亮.Ansoft 12在工程電磁場中的應(yīng)用[M].北京：中國水利水電出版社，2010.
猜你喜歡
尾錐角對彈體斜侵徹過程中姿態(tài)的影響研究兵器裝備工程學(xué)報(bào)(2022年7期)2022-08-10 03:37:18 寬頻高磁導(dǎo)率R10k軟磁材料的開發(fā)山東冶金(2022年2期)2022-08-08 01:50:52 橢圓截面彈體斜侵徹金屬靶體彈道研究*爆炸與沖擊(2022年2期)2022-03-17 07:28:44 基于FEMM的永磁電機(jī)動(dòng)態(tài)凍結(jié)磁導(dǎo)率并行仿真及程序防爆電機(jī)(2021年3期)2021-07-21 08:12:50 基于遞推最小二乘法加速度計(jì)信息輔助的磁力計(jì)標(biāo)定方法?傳感技術(shù)學(xué)報(bào)(2021年12期)2021-03-13 07:17:52 六軸IMU補(bǔ)償?shù)拇帕τ?jì)動(dòng)態(tài)穩(wěn)定校準(zhǔn)儀表技術(shù)與傳感器(2021年1期)2021-02-25 11:52:36 基于低功耗藍(lán)牙的磁力計(jì)在線校準(zhǔn)系統(tǒng)設(shè)計(jì)中國科技縱橫(2019年8期)2019-06-19 01:07:24 STOPAQ粘彈體技術(shù)在管道施工中的應(yīng)用上海煤氣(2018年6期)2018-03-07 01:03:22 Advanced Frozen Permeability Technique and Applications in Developing High Performance Electrical Machines電工技術(shù)學(xué)報(bào)(2016年20期)2016-11-17 07:33:45 鋼板磁導(dǎo)率變化對船舶感應(yīng)磁場的影響中國艦船研究(2014年6期)2014-05-14 06:45:29
王麗黎，尹文超
(西安理工大學(xué) 自動(dòng)化與信息工程學(xué)院，陜西 西安 710048)
在地磁匹配導(dǎo)航系統(tǒng)中，載體自身干擾磁場的分析與補(bǔ)償是一個(gè)關(guān)鍵問題。根據(jù)飛行器上的磁干擾的類型可以建立磁干擾的數(shù)學(xué)模型，模型中含有一系列的系數(shù)，在飛行器上的儀器都固定的條件下，這些系數(shù)是固定的。測得需要的數(shù)據(jù)就可以推導(dǎo)出磁補(bǔ)償系數(shù)，利用這些系數(shù)可以求出飛行器周圍引起的干擾磁場大小，從而在地磁場未知的情況下得出真實(shí)的地磁場值。
目前國內(nèi)在磁場補(bǔ)償方面的研究主要集中在概略分析[1]和算法的研究[2-4]，而針對基于電磁場仿真模型和數(shù)值計(jì)算方面的研究卻很少，少數(shù)單位僅僅研究了縮比模型周圍的磁場分布[5],還沒有基于真實(shí)模型的仿真補(bǔ)償研究。針對目前在磁補(bǔ)償方面還沒有電磁仿真研究的情況，本文使用Ansoft Maxwell有限元軟件，仿真出物體在不同磁導(dǎo)率的情況下以及不同磁力計(jì)放置位置時(shí)的測量總場值，依據(jù)這些測量值可以求出彈體外表面任意點(diǎn)及材料在任意磁導(dǎo)率下的補(bǔ)償系數(shù)，并利用這些系數(shù)估計(jì)出背景地磁場。
1 導(dǎo)彈磁場補(bǔ)償模型分析與軟件建模
1.1 彈體磁場模型求解分析
根據(jù)飛行器磁干擾的Tolles-Lawson數(shù)學(xué)模型及表達(dá)式[6]作為分析模型，本模型以仿真感應(yīng)磁場為主，暫不考慮恒定場和渦流場，建立彈體坐標(biāo)系如圖1所示。L軸與彈體縱軸平行，首部為正；T軸與L軸垂直，V軸與TOL平面正交，向下為正。α、β、γ分別為 He與3個(gè)坐標(biāo)軸之間的夾角， 其方向余弦分別為 cosα、cosβ、cosγ。He為地磁場矢量，感應(yīng)出來的干擾場投影到地磁場方向大小為 Hid，化簡后為表達(dá)式(1)：
其中感應(yīng)場有 6個(gè)補(bǔ)償系數(shù) (TT-LL、VV-LL、LL、LT+TL、VT+TV、VL+LV),磁力計(jì)測量總場 Hz減去地磁場 He可得干擾場Hid，再利用地磁場矢量與彈體坐標(biāo)系3軸夾角和干擾場大小來求出補(bǔ)償系數(shù)。
圖1 彈體坐標(biāo)系示意圖
1.2 模型的建立
在Maxwell中建立Tolles-Lawson模型，彈體總長度為7 m，前半部分彈頭為圓錐型總長度3 m，后半部分彈體包括尾翼總長度4 m，半徑為30 cm，大小和形狀接近真實(shí)模型，使用永磁體加上適當(dāng)?shù)某C頑力來形成一定方向的均勻磁場，如圖2所示。彈體在均勻磁場會(huì)受到磁化并在周圍產(chǎn)生感應(yīng)磁場。由于軟件所仿真出的磁場值結(jié)果輸出都以Global坐標(biāo)系為準(zhǔn)，在軟件仿真中通過改變地磁場對物體的作用方向來模擬彈體在空間內(nèi)與地磁場形成任意角度，將軟件輸出值作為模擬磁力計(jì)的測量值。模型的材料及屬性參數(shù)如表1所示。
圖2 仿真模型
表1 模型的材料屬性
2 磁場補(bǔ)償技術(shù)的仿真與計(jì)算
2.1 仿真原理和計(jì)算步驟
磁干擾表達(dá)式(1)中有6個(gè)未知補(bǔ)償系數(shù)，求出方程的解最少需要6個(gè)方程，將式(1)寫成矩陣形式(2),左列為彈體在地磁場n種姿態(tài)時(shí)的干擾場大小，右側(cè)為對應(yīng)的姿態(tài)矩陣和感應(yīng)系數(shù)的乘積，這里使用最小二乘法求解。
仿真時(shí)令初始磁場方向垂直XOY面向下,再令磁場方向相對Y軸逆時(shí)針旋轉(zhuǎn)形成參數(shù)為α的角度，之后相對Z軸逆時(shí)針旋轉(zhuǎn)形成參數(shù)為β角度,如圖3所示。
圖3 模型的側(cè)面和縱截面
本文研究磁力計(jì)安裝在彈體表面的情況，通過軟件來實(shí)現(xiàn)多點(diǎn)采樣。這樣可以找出磁力計(jì)安裝在某個(gè)位置使得補(bǔ)償效果達(dá)到最好，本例取距離彈體表面5 cm沿彈體縱軸線的一條路經(jīng)作為磁力計(jì)放置一系列的測量位置，目的是要算出這條路徑上所有采樣點(diǎn)的補(bǔ)償系數(shù)，最后再求出所有位置的補(bǔ)償誤差，確定補(bǔ)償效果最佳的磁力計(jì)放置位置。
在軟件求解前需要進(jìn)行參數(shù)設(shè)置,如表2所示。
表2 求解參數(shù)設(shè)置
感應(yīng)磁場是在靜態(tài)的恒定磁場作用下產(chǎn)生的，因此選用靜磁場求解器。零邊界條件保證了外部沒有磁場并且內(nèi)部磁場方向單一且垂直于邊界面。參數(shù)化α1和β1的目的是形成對物體不同作用方向的磁場，每個(gè)組合對應(yīng)一個(gè)彈體在地磁場中的一個(gè)姿態(tài)，由于彈體模型分別關(guān)于XOY面和XOZ面對稱并且地磁場反向作用時(shí)僅僅改變最后磁場值的符號，所以這里參數(shù)化0～90代表地磁場在一個(gè)象限內(nèi)作用物體，通過maxwell的參數(shù)設(shè)置step可以包含任意磁場對彈體的作用方向，以此來模擬飛行器在空間中的姿態(tài)。參數(shù)化 α2和β2形成的磁場作為驗(yàn)證，利用前部分求出的系數(shù)和后部分求出的總場來估計(jì)地磁場大小，最后將估計(jì)出的磁場值與形成的背景地磁場相比較來評價(jià)補(bǔ)償結(jié)果。將上述參數(shù)設(shè)置好后即可利用軟件進(jìn)行數(shù)值計(jì)算。
2.2 數(shù)據(jù)的后處理
導(dǎo)出測磁補(bǔ)償路徑上所有點(diǎn)在相對磁導(dǎo)率為1～5時(shí)彈體處于地磁場中所有姿態(tài)的地磁場三分量BXBYBZ和模值MagB，在這條路徑上采樣700個(gè)點(diǎn)作為磁探頭的測量點(diǎn)，磁場大小為50 266 nt，同時(shí)計(jì)算出式(2)中的cosX、cosY、cosZ,表達(dá)式左邊的干擾場大小等于該磁探頭位置的測量值MagB減去地磁場大小得出。利用角度 α1、β1、step=10時(shí)參數(shù)化的結(jié)果求出所有位置的 6個(gè)補(bǔ)償系數(shù)，再利用角度 α2、β2、step=36時(shí)參數(shù)化的結(jié)果求出所有位置的補(bǔ)償系數(shù)來估計(jì)背景地磁場，并與真值He對比得出每個(gè)位置補(bǔ)償后的剩余標(biāo)準(zhǔn)差，最后對補(bǔ)償結(jié)果進(jìn)行分析。
3 對仿真方法可靠性的驗(yàn)證
下面建立一個(gè)簡單的靜態(tài)磁偶極子模型對軟件的計(jì)算精度進(jìn)行驗(yàn)證和分析。之后,以一個(gè)簡單實(shí)驗(yàn)來驗(yàn)證利用本文的補(bǔ)償方法還原真實(shí)地磁場。
3.1 驗(yàn)證精度
Ansoft公司的Maxwell 3D是一個(gè)功能強(qiáng)大、結(jié)果精確、易于使用的三維電磁場有限元分析軟件。
根據(jù)電磁場理論，磁偶極子模型[7]在側(cè)向空間中形成的磁感應(yīng)強(qiáng)度大小為式(3)，軸線方向形成的磁感應(yīng)強(qiáng)度大小為式(4)：
根據(jù)這一理論在磁場中建立磁偶極子模型，其磁場分布如圖4所示。將仿真出來的有限元解與理論值比較，如表3所示。有限元解與理論值的誤差百分比很低，說明數(shù)值計(jì)算精度較高，有限元解較接近實(shí)際理論值。
3.2 實(shí)驗(yàn)分析
實(shí)驗(yàn)在一個(gè)可以形成單一方向均勻磁場的磁屏蔽箱內(nèi)進(jìn)行。
實(shí)驗(yàn)時(shí)保證磁力計(jì)與物體相對位置不變，通過增加亥姆霍茲線圈的電流控制最終形成了51 807 nt的電磁場，磁場沿磁屏蔽箱的軸線方向，通過不斷改變物體在磁場內(nèi)的位置狀態(tài)測得每個(gè)狀態(tài)對應(yīng)的磁場三分量和總場強(qiáng)度。這個(gè)實(shí)驗(yàn)在一個(gè)平面內(nèi)進(jìn)行一些角度的旋轉(zhuǎn)，測量值、補(bǔ)償值和真實(shí)值的結(jié)果如圖5所示。
圖5 測量和補(bǔ)償結(jié)果
圖4 磁偶極子模型在空間中產(chǎn)生的磁感應(yīng)強(qiáng)度
表3 有限元解與理論解比較
圖5顯示了物體在磁場中36個(gè)狀態(tài)的補(bǔ)償結(jié)果，平均誤差比較小，直觀地反映了本文的補(bǔ)償方法是可行的。
該實(shí)驗(yàn)說明，通過不斷改變被補(bǔ)償物體在磁場中的不同狀態(tài)并利用最小二乘法可以近似地估計(jì)外部地磁場的大小。
4 仿真結(jié)果及其分析
在Maxwell軟件中仿真生成的背景地磁場He的大小為50 266 nt，在補(bǔ)償路徑上采樣700個(gè)測量點(diǎn)作為測磁點(diǎn)，同時(shí)計(jì)算了物體在磁導(dǎo)率從2～5每個(gè)點(diǎn)的補(bǔ)償結(jié)果，對結(jié)果的分析評價(jià)如圖6所示，圖表的縱坐標(biāo)為對彈體表面點(diǎn)在參數(shù)化(α2和β2)姿態(tài)時(shí)的補(bǔ)償剩余標(biāo)差。
圖6 某路徑在同一地磁場作用下不同磁導(dǎo)率下的補(bǔ)償結(jié)果
通過對仿真數(shù)據(jù)的處理，得出了沿著彈體表面某路徑在同一地磁場作用下不同磁導(dǎo)率的補(bǔ)償結(jié)果，同時(shí)得出了材料在相對磁導(dǎo)率在2～5時(shí)對多種飛行姿態(tài)補(bǔ)償后的平均誤差。其一，磁導(dǎo)率越大，補(bǔ)償后的誤差相對偏大，但并不是所有的補(bǔ)償位置的補(bǔ)償誤差都隨著磁導(dǎo)率的增大而增加；其二，基于本模型，彈體首部和彈體表面的補(bǔ)償誤差偏小，基本小于100 nT，而尾部的補(bǔ)償誤差較大。
本文通過電磁場仿真軟件Ansoft maxwell仿真飛行載體空間中的感應(yīng)磁場分布，并使用最小二乘法(LSE)進(jìn)行參數(shù)估計(jì)，求出彈體表面空間某曲線上所有采樣點(diǎn)的補(bǔ)償系數(shù)并進(jìn)行分析。本文提供了一種利用電磁場仿真軟件建立Tolles-Lawson模型并補(bǔ)償?shù)姆椒?，可以在不進(jìn)行試驗(yàn)的情況下仿真研究感應(yīng)磁場分布，實(shí)現(xiàn)了磁場補(bǔ)償，并對多種姿態(tài)的補(bǔ)償結(jié)果進(jìn)行了綜合評價(jià)。通過軟件參數(shù)化設(shè)置，對仿真數(shù)據(jù)進(jìn)行處理之后可以找到不同磁導(dǎo)率下模型的最佳補(bǔ)償點(diǎn)。本仿真方法可以對任意形狀模型在任意磁導(dǎo)率的情況下找出最佳補(bǔ)償位置，以及磁力計(jì)在同一安裝位置下找出補(bǔ)償效果最好的材料相對磁導(dǎo)率，在工程中可為材料磁導(dǎo)率的選取和磁力計(jì)的安裝位置提供依據(jù)。
[1]趙雙求.UAV航磁系統(tǒng)及補(bǔ)償研究[D].長沙：中南大學(xué)，2011.
[2]張堅(jiān),林春生,羅青,等.基于預(yù)測殘差平方和的飛機(jī)干擾磁場模型求解[J].探測與控制學(xué)報(bào)，2010,32(5):74-78.
[3]龐學(xué)亮,林春生,張寧.飛機(jī)磁場模型系數(shù)的截?cái)嗥娈愔捣纸夥ü烙?jì)[J].探測與控制學(xué)報(bào),2009,31(5):48-51.
[4]伍東凌,陳正想,王秀.基于遺傳算法的磁干擾補(bǔ)償方法[J].探測與控制學(xué)報(bào)，2012,34(6):16-20.
[5]李永慧，張曉明，劉俊,等.某常規(guī)彈載體感應(yīng)磁場分布的仿真技術(shù)研究[J].彈箭與制導(dǎo)學(xué)報(bào)，2012,32(4):181-184.
[6]LELIAK P.Identification and evaluation of magnetic field sources of magnetic airborne detector equipped aircraft[J].IRE Transactions on Airspace and Navigational Electronics，1961(8)：95-105.
[7]任來平,趙俊生,侯世喜.磁偶極子磁場空間分布模式[J].海洋測繪,2002,22(2):18-21.
[8]趙博，張洪亮.Ansoft 12在工程電磁場中的應(yīng)用[M].北京：中國水利水電出版社，2010.
猜你喜歡
兵器裝備工程學(xué)報(bào)(2022年7期)2022-08-10 03:37:18
山東冶金(2022年2期)2022-08-08 01:50:52
爆炸與沖擊(2022年2期)2022-03-17 07:28:44
防爆電機(jī)(2021年3期)2021-07-21 08:12:50
傳感技術(shù)學(xué)報(bào)(2021年12期)2021-03-13 07:17:52
儀表技術(shù)與傳感器(2021年1期)2021-02-25 11:52:36
中國科技縱橫(2019年8期)2019-06-19 01:07:24
上海煤氣(2018年6期)2018-03-07 01:03:22
電工技術(shù)學(xué)報(bào)(2016年20期)2016-11-17 07:33:45
中國艦船研究(2014年6期)2014-05-14 06:45:29
