吳 穹，王瑞林，向紅軍

(軍械工程學院 彈藥工程系，河北 石家莊 050003)

基于正交試驗法的單級磁阻線圈發(fā)射器驅動線圈參數優(yōu)化

吳 穹，王瑞林，向紅軍

(軍械工程學院 彈藥工程系，河北 石家莊050003)

為了減少單級磁阻線圈發(fā)射器仿真過程中仿真試驗的任務量，將正交試驗法和仿真分析有機結合。正交試驗法能夠用少數組合搭配均衡表示全部試驗，利用正交表設計驅動線圈參數優(yōu)化仿真試驗，從全部27組試驗中選出具有代表性的9組進行有限元分析，根據仿真結果選出最優(yōu)的參數進行試驗驗證。結果表明：正交試驗法可以應用于單級磁阻線圈發(fā)射器的仿真分析，仿真分析結果為今后的單級磁阻線圈發(fā)射器提供參考依據,也為多級磁阻型線圈發(fā)射器的仿真分析提供了新思路。

電磁學；單級磁阻線圈發(fā)射器；有限元法；正交試驗法；參數優(yōu)化

電磁發(fā)射器利用電磁力推進射彈，主要分為導軌型、線圈型和重接型，其中線圈型又可以分為感應型和磁阻型。磁阻型線圈發(fā)射器是利用線圈中鐵磁磁路的磁阻變化吸引鐵芯向磁阻最低的位置運動來加速鐵芯彈丸的[1]。磁阻型線圈發(fā)射器本質是一種直線電機，它具有速度可控性好，發(fā)射過程隱蔽性和安全性高等優(yōu)點,并且可以通過調節(jié)電壓、電容值、觸發(fā)級數、觸發(fā)位置等參數調節(jié)彈丸出口速度[2]。上述參數與彈丸出口速度之間的關系無法用數學模型表示，需要進行大量仿真與試驗來實現(xiàn)參數的優(yōu)化。目前，遺傳算法[3-4]、蟻群算法[5]、退火算法[6]等都已經應用于線圈發(fā)射器的優(yōu)化之中。但隨著參數和級數增加，其優(yōu)化時間變長，且優(yōu)化算法無法反映單個參數對系統(tǒng)影響。采用單因素分析法對電參數、驅動線圈結構、彈丸結構以及觸發(fā)位置的研究，往往是在特定的發(fā)射器結構情況下，研究單一參數對系統(tǒng)性能的影響，沒有考慮到各個參數間的相互作用。將最優(yōu)單因素簡單組合難以得到最優(yōu)方案。而正交試驗法不但能夠得到較優(yōu)的設計方案，還能分析各個因素對整個系統(tǒng)的影響[7]。

1 磁阻型線圈發(fā)射器原理

單級磁阻線圈發(fā)射器由驅動線圈、鐵磁性彈丸、發(fā)射管、脈沖儲能電容器、二極管和一個手動觸發(fā)開關組成。依據磁阻最小原理，磁通總是趨向于經過磁阻最小的路徑。驅動線圈、鐵磁性彈丸和空氣組成的磁路中磁阻的變化恰好滿足磁阻最小原理，可以吸引彈丸加速運動。由于鐵磁性彈丸的磁導率遠大于空氣的磁導率，因此在線圈、彈丸和空氣組成的磁路里，彈丸會向著磁阻最小的方向運動。亦可理解為被磁化的鐵磁性彈丸中的磁化電流與驅動線圈中的電流具有相同方向，從而使彈丸受吸力運動。

當電流流經線圈，磁能存儲在線圈周圍的磁場中。磁通由磁動勢的表達式定義

Fm=Ni=ΦR

(1)

式中：Fm為磁動勢;N為驅動線圈安匝數;i為電流大小;Ф為磁通;R為線圈周圍磁路的磁阻。

在該模型中，磁阻R是常量，磁通Ф與Ni的關系是線性的，系統(tǒng)儲存的磁能Em[8]是磁通Ф與Ni關系曲線圍成面積，定義為

(2)

磁阻R的一般方程為

(3)

式中：l為磁路長度；μ0為真空磁導率；μr為相對磁導率；A為磁通穿過的面積(磁路截面)。

這里描述的模型，幾乎所有的磁阻都出現(xiàn)在驅動線圈與鐵磁性彈丸之間的空氣間隙中。當彈丸中心與線圈中心重合時，考慮到空氣的相對磁導率μr為1，則磁阻R的表達式可以表示為

(4)

式中：g為驅動線圈和鐵磁性彈丸之間的氣隙；d是發(fā)射管直徑；p為鐵磁性彈丸長度；μ0為真空磁導率。

此時鐵磁性彈丸獲得的加速力可以表示為

(5)

由于式(5)在推導過程中使用線性化方法解決非線性問題，并且沒有考慮漏磁的影響，所以只能用于非精確計算彈丸受力情況，若想精確分析還需使用有限元法。

2 正交試驗法描述

正交試驗法[9]是由日本田口玄一等根據試驗優(yōu)化規(guī)律發(fā)明的一種基于正交表安排與分析多因素試驗的一種設計方法。該方法選用具有代表性的水平組合進行試驗，通過對結果分析了解全面試驗情況，得出最優(yōu)組合。正交試驗法具有正交性、代表性和綜合可比性，用正交表設計的試驗具有均衡分散和整齊可比的特點[10-13]。磁阻型電磁發(fā)射器的驅動線圈是多個參數相互作用的系統(tǒng)，應用正交試驗法進行參數優(yōu)化，對結果的極差和方差進行分析，可在考慮線圈各參數相互作用的基礎上得到線圈參數對彈丸出口速度的影響，較單因素分析更為準確。正交試驗法的設計流程圖如圖1所示。

3 有限元模型建立

驅動線圈自感、彈丸自感以及二者互感等參數相互影響的非線性計算較為復雜，利用Ansoft軟件可以解決非線性計算問題[14]。驅動線圈和彈丸均為軸對稱結構，假設彈丸與驅動線圈的軸線始終重合，即在加速過程中不考慮彈丸偏離驅動線圈軸線的情況，采用Ansoft軟件中的軸對稱坐標系來分析問題，可以將發(fā)射器的三維結構簡化為二維軸對稱模型，這種方法在減少了計算量的同時又確保了仿真結果的可信度[15]，發(fā)射器仿真模型圖如圖2所示。

模型由驅動線圈、彈丸、運動域和求解域4部分構成。彈丸直徑為9.2mm，長度為34mm，材料為steel_1010，驅動線圈根據試驗選取不同的直徑，設置不同的長度和層數，進而確定匝數，材料設置為銅；運動域與求解域材料設置為空氣。為確保仿真精度，在進行網格劃分時，彈丸網格最為密集，其次為驅動線圈和運動域，求解域網格劃分最稀疏。仿真過程中，外部電路如圖3所示，為線圈繞組提供激勵源。由于不考慮電參數的影響，所以電容的電容值和充電電壓值保持不變，驅動線圈的電阻值根據測量設定。仿真時間始于0時刻，止于6ms，步長為0.02ms。

在進行激勵源設置時，選擇由外部電路為線圈繞組提供激勵源，其中外電路的電路電阻仿真軟件無法自行計算，采用TE2810C型LCR數字電橋對電路電阻進行測量以提高仿真精度。

4 正交試驗結果及分析

試驗的變化參數有銅線直徑A、繞制層數B和線圈長度C，每個因素有3個水平1、2、3，根據試驗因素和水平選取L9(34)正交表[10]構成本試驗的因素水平表，如表1所示。

表1 因素水平表

利用正交設計助手軟件，進行仿真試驗，仿真結果的直觀分析表如表2所示。

表2 直觀分析表

依據表2設計的9組試驗進行仿真試驗，仿真過程中彈丸速度隨時間變化規(guī)律如圖4所示。

正交試驗的分析方法通常有極差分析法和方差分析法。前者具有簡單直觀，計算量小等優(yōu)點，但是極差分析法不能估計誤差的大小，不能精確估計各因素對試驗結果影響的重要程度，特別是對水平數大于等于3且考慮交互作用的試驗，不便使用。如果再對結果進行方差分析，就能彌補極差分析法的不足[11]。根據表2中極差與方差的計算和分析結果可知：

1)最大彈丸初速：由直觀分析表可以看出，實驗3最大彈丸初速為22.45 m/s，發(fā)射器的最優(yōu)水平組合初選表如表3所示。

表3 最優(yōu)水平組合初選表

2)極差方差分析：3個因素的極差和方差計算結果如表2所示，各因素的效應曲線圖如圖5所示。各個因素對發(fā)射效果即彈丸初速的影響重要性主次順序依次為：繞制層數、銅線直徑、線圈長度。

由圖5分析可知，隨著繞制層數增加，彈丸初速明顯增大；在相同線圈尺寸的情況下銅線直徑的減小可以增加線圈安匝數，從而使彈丸初速增大，但增加效果不如繞制層數明顯；線圈長度對發(fā)射效果影響不明顯。其中銅線直徑和線圈長度的效應曲線圖變化不明顯，這可能是由于水平數量過少或二者之間存在交互作用導致，可以設計水平數較多的單因素試驗和用考慮交互作用的正交試驗進行研究分析。

3)預計最優(yōu)參數驗證：根據結果分析表和效應曲線圖分析另兩組水平組合也可能達到最優(yōu)，分別是a組：銅線直徑0.6 mm，繞制9層，線圈長度55 mm；b組：銅線直徑1.0 mm，繞制9層，線圈長度60 mm。分別對這兩組水平組合進行計算，得到的速度分別為22.87 m/s和21.55 m/s，a組的速度大于試驗3的速度，其原因是，線圈放電時間過長，會有短暫的反向力作用，表現(xiàn)為彈丸達到最大速度后略有減小，最終達到穩(wěn)定。將最優(yōu)水平組合整理成表格如表4所示。

表4 最優(yōu)水平組合表

以上結果表明，在進行單級磁阻線圈發(fā)射器驅動線圈參數選擇時，應首先考慮繞制層數，因為該因素相較于其他兩個因素對發(fā)射效果影響最大。在實際工程應用過程中線圈的繞制層數和線圈長度受發(fā)射器外形以及集成時難易影響，減小銅線直徑可以減小驅動線圈外形,但要考慮安全電流和散熱問題。因此，在設計過程中，確定繞制層數和銅線直徑這兩個因素是后續(xù)設計的基礎。

5 結論

筆者基于正交試驗法，利用Ansoft仿真軟件對單級磁阻發(fā)射器的驅動線圈的參數進行優(yōu)化和分析，得出以下結論：

1)正交試驗法可以應用于單級磁阻發(fā)射器驅動線圈的參數優(yōu)化，得到較為理想的設計參數，成功減少了試驗次數。這對今后其他參量的優(yōu)化和多級磁阻線圈發(fā)射器的設計提供了新的思路。

2)繞制層數對彈丸速度影響最大，隨著繞制層數的增加，彈丸初速明顯增加；銅線直徑次之，隨著銅線直徑的減小，彈丸初速有所增加；線圈長度對彈丸初速影響最小，無明顯的遞變規(guī)律。

References)

[1] 王瑩，肖峰.電炮原理[M].北京：國防工業(yè)出版社，1995:5-15. WANG Ying, XIAO Feng. Principle of electric gun[M].Beijing: National Defense Industry Press,1995:5-15.(in Chinese)

[2] 支彬安，雷彬，李治源，等.非線性等效電感對磁阻驅動器驅動性能的影響[J].彈箭與制導學報，2011，31(1)：173-176. ZHI Bin’an, LEI Bin, LI Zhiyuan, et al. Equivalent inductance’s influence on the performance of reluctance driver[J]. Journal of Projectiles，Rockets，Missiles and Guidance,2011,31(1):173-176. (in Chinese)

[3] GUO L M ,GUO N N,WANG S H. Optimization for capacitordriven coilgun based on equivalent circuit model and genetic algorithm[C]∥2009IEEE Conference on Energy Conversion Congress and Exposition. San Jose， CA：IEEE,2009:234-239.

[4] KIM S W,JUNG H K,HAHN S Y. An optimal design of capacitor-driven coilgun[J].IEEE Transactions on Magnetics,1994,30(2):207-211.

[5] LIU W B,CAO Y J,ZHANG Y, et al. Parameters optimization of synchronous induction coilgun based on ant colony algorithm[J].IEEE Transactions on Plasma Scie-nce,2011,39(1):100-104.

[6] HAINSWORTH G,RODGER D. Design optimization of coilgun[J].IEEE Transactions on Magnetics,1995,31(1):473-477.

[7] 張亞東，肖剛，王青子，等.基于正交試驗法的同步感應線圈發(fā)射器關鍵設計參數選擇[J].高電壓技術，2016，42(9):2843-2849. ZHANG Yadong, XIAO Gang, WANG Qingzi,et al. Key design parameters selection of synchronous induction coil launcher based on orthogonal experiment[J].High Voltage Engineering,2016,42(9):2843-2849. (in Chinese)

[8] BRESIE D A,ANDREWS J A.Design of reluctance accele-rator[J].IEEE Transaction on Magnetics ,1991,27(1):623-627.

[9] 邱軼兵.試驗設計與數據處理[M].合肥：中國科學技術大學出版社，2008. QIU Yibing. Experiment design and data processing[M]. Hefei:Press of Univessity of Science and Technology of China,2008. (in Chinese)

[10] 中國科學院.正交實驗方法[M].北京：人民教育出版社，1975. Chinese Academy of Science. Orthogonal experiment method[M].Beijing: People’s Education Press,1975. (in Chinese)

[11] 徐天爽.基于面向對象有限元程序設計的正交試驗及其在結構優(yōu)化中的研究[D].長春：吉林大學，2009. XU Tianshuang. Study of orthogonal experiments based on object-oriented finite element programing and its application in structure optimization[D].Changchun:Jilin University,2009. (in Chinese)

[12] LI Z H,MALIK O P. An orthogonal test approach based control parameter optimization and its application to a hydro-turbine governor[J].IEEE Transaction on Energy Conversion,1997,12(4):388-393.

[13] BANERJI S. Orthogonal array approach for test case optimization[J].Intenational Journal of Advanced Research in Computer and Communication Engineering,2012,1(9):613-621.

[14] 趙科義，李治源，程樹康.單級感應線圈炮工作過程的動態(tài)仿真[J].高壓電技術，2008,34(8):1667-1671. ZHAO Keyi, LI Zhiyuan, CHENG Shukang. Dynamic simulation of working process of the single-stage induction coil-gun[J]. High Voltage Engineering,2008,34(8):1667-1671. (in Chinese)

[15] 趙博，張洪亮.Ansoft12在工程電磁場中的應用[M].北京：中國水利水電出版社，2010:130-153. ZHAO Bo, ZHANG Hongliang. Application of Ansoft12on engineering electromagnetic field[M]. Beijing: China Water & Power Press,2010:130-153. (in Chinese)

ParameterOptimizationofSingle-stageReluctanceCoilLauncherBasedonOrthogonalExperiment

WU Qiong , WANG Ruilin , XIANG Hongjun

(Department of Ammunition Engineering, Ordnance Engineering College, Shijiazhuang050003, Hebei，China)

Orthogonal test method and simulation analysis are combined organically for the purposes of reducing the work of simulation experiment in single-stage reluctance coil launcher simulation. A small number of combinations with a balanced representation of all experiments can be used in orthogonal experiment. Orthogonal experiments are used to design driving coil parameter optimization simulation experiment. A representative group of9groups was selected from all27groups to perform finite element analysis. According to the simulation results, experimenting the optimal parameters. The results show that orthogonal experiments can be applied to the simulation of single-stage reluctance coil launcher. The simulation results provide references for the future single-stage reluctance coil launcher, and also provide a new train of thought for the simulation of multi-stage reluctance coil launcher.

electromagnetics; single-stage reluctance coil launcher; finite element method; orthogonal experiment; parameter optimization

TJ27+1

： A

：1673-6524(2017)03-0007-05

10.19323/j.issn.1673-6524.2017.03.002

2016-12-30

吳穹(1993—)，男，碩士研究生，主要從事輕武器和電磁技術研究。E-mail：wuqiongpanda@163.com