王宗偉，任志彬，常志鵬，劉建設(shè)，崔鵬飛

?

基于ANSYS和AMESim的螺管電磁鐵吸力仿真研究

王宗偉，任志彬，常志鵬，劉建設(shè)，崔鵬飛

（空間物理重點實驗室，北京，100076）

對常開式螺管電磁鐵，分別采用有限元軟件ANSYS、系統(tǒng)參數(shù)仿真軟件AMESim對螺管電磁鐵通電工作時的電磁場分布、磁感應(yīng)強度分布以及電磁鐵的吸力特性進行仿真計算，得到電磁吸力隨源電壓的變化曲線及電磁吸力隨工作氣隙的變化曲線，同時將仿真分析結(jié)果與吸力試驗實測結(jié)果進行了對比分析。分析結(jié)果表明：ANSYS和AMESim的仿真結(jié)果與試驗實測值基本一致，均可實現(xiàn)對電磁鐵吸力特性的可靠分析，而AMESim的優(yōu)勢在于對電磁鐵性能的快速評估，進而縮短設(shè)計周期。

螺管電磁鐵；電磁吸力；仿真分析

0 引 言

在運載火箭/導(dǎo)彈/航天器動力系統(tǒng)中，電磁閥的主要作用是通過通電、斷電的工作模式，實現(xiàn)流體通路的開啟和關(guān)閉。其典型應(yīng)用包括增壓輸送系統(tǒng)中控制高壓氣體對貯箱增壓，姿控發(fā)動機中控制推進劑流動以實現(xiàn)發(fā)動機啟動和關(guān)機，或控制其它氣動閥實現(xiàn)各種程序動作[1,2]。應(yīng)用在液體軌姿控推進系統(tǒng)中的電磁閥一般為螺管式電磁閥，主要由閥體、閥芯、線圈、圓柱彈簧以及蝶形彈簧等組成[3]。螺管電磁鐵作為電磁閥的驅(qū)動器，其作用是在擋鐵和閥芯之間產(chǎn)生電磁吸力，進而控制電磁閥的打開和關(guān)閉，其性能對整個電磁閥的特性有重要的影響。電磁吸力是電磁閥的重要指標(biāo)之一，無論是穩(wěn)態(tài)控制精度和動態(tài)響應(yīng)性能，還是抗干擾能力或工作可靠性都在很大程度上取決于電磁鐵的吸力特性的綜合性能[4]。相關(guān)行業(yè)內(nèi)對電磁閥的電磁吸力一般采用工程算法進行計算分析，采用有限元軟件或AMESim系統(tǒng)仿真軟件對電磁吸力特性進行研究并不十分廣泛。龐末紅等利用電磁有限元軟件Ansoft Maxwell對電磁鐵的靜特性和動特性進行仿真計算，得出磁感應(yīng)強度、磁力線分布圖和靜磁場下的吸力特性曲線[5]；婁路亮對計算直流螺線管式電磁鐵電磁吸力的經(jīng)驗公式法、磁路分割法和有限元方法進行了對比，分析了3種方法的計算難易程度、計算精度等[6]。

本文針對一種常開式螺管電磁鐵，利用ANSYS有限元分析軟件和AMESim系統(tǒng)參數(shù)仿真軟件對其靜特性進行仿真分析，得到電磁鐵的磁感應(yīng)強度、磁力線分布和吸力特性曲線，并將仿真結(jié)果與試驗實測值進行了對比分析。

1 電磁鐵結(jié)構(gòu)及方案設(shè)計

螺管電磁鐵主要由殼體、線圈、擋鐵、銜鐵和套筒組成，其中殼體為電工純鐵DT4，擋鐵、銜鐵材料為軟磁合金1J50，套筒材料為不銹鋼，其結(jié)構(gòu)示意圖如圖1所示。本文采用常開式螺管電磁鐵方案，其工作原理為：電磁鐵通電時，線圈組件與銜鐵組成磁路產(chǎn)生的吸力使銜鐵向上運動，電磁閥進氣通道打開；電磁鐵斷電時，線圈組件與銜鐵組成磁路產(chǎn)生的吸力消失，銜鐵在副彈簧力的作用下向下運動，電磁閥進氣通道關(guān)閉。

圖1 螺管電磁鐵結(jié)構(gòu)示意

根據(jù)電磁場原理，計算室溫下電磁吸力的Maxwell公式如下[7~9]：

考慮到電磁鐵的實際結(jié)構(gòu)，引入相關(guān)設(shè)計系數(shù)，推導(dǎo)得到如下電磁吸力計算公式：

式中M為側(cè)向吸力摩擦力系數(shù)，無量綱；為電源電壓，V；CT為磁勢降系數(shù)，無量綱；20為20 ℃時線圈電阻，Ω；為線圈匝數(shù)，匝；0為真空磁導(dǎo)率，0=4π×10-7H/m；1為銜鐵半徑，mm；0為銜鐵中心孔半徑，mm；K為工作氣隙，mm。

電磁鐵電磁吸力的設(shè)計值為121N，根據(jù)式（2）對電磁鐵結(jié)構(gòu)參數(shù)進行計算并開展結(jié)構(gòu)設(shè)計。

2 電磁鐵的仿真計算

為了研究螺管電磁鐵通電工作時的電磁場分布、磁感應(yīng)強度分布以及電磁鐵的吸力特性，分別采用有限元軟件ANSYS和系統(tǒng)仿真軟件AMESim對電磁鐵進行仿真計算。

2.1 有限元軟件仿真

利用ANSYS軟件對螺管電磁鐵靜態(tài)電磁力進行仿真計算[10~12]，并按照下列步驟建模仿真：

a）定義單元類型，包括兩種不同自由度的PLANE53單元，均為軸對稱。

b）為自由度為AZ CURR的PLANE53單元定義實常數(shù)，包括線圈截面積、線圈匝數(shù)、線圈充填系數(shù)等。

c）定義材料屬性，其中空氣、線圈與套筒為不導(dǎo)磁材料，殼體、擋鐵和銜鐵為導(dǎo)磁材料，電工純鐵DT4和軟磁合金1J50的-曲線如圖2所示。

a）電工純鐵DT4

b）軟磁合金1J50

圖2 材料的曲線

d）按照設(shè)計方案的結(jié)構(gòu)參數(shù)建立幾何模型，該電磁鐵為軸對稱結(jié)構(gòu)，故建立電磁鐵的二維軸對稱模型，具體結(jié)構(gòu)如圖3所示。

e）為各面域賦予不同材料并指定單元類型。

f）劃分四邊形網(wǎng)格，網(wǎng)格結(jié)構(gòu)如圖4所示。

圖3 幾何模型

圖4 二維網(wǎng)格模型

g）將銜鐵單元定義為組件，并施加磁力邊界條件。

h）將線圈單元定義為組件，對其節(jié)點耦合電流自由度，并對單元施加電壓激勵。

i）對模型的外部節(jié)點施加磁力線平行邊界條件。

j）運行穩(wěn)態(tài)仿真。

根據(jù)仿真計算，當(dāng)工作氣隙為0.75 mm、源電壓為24 V時，計算得到螺管電磁鐵的磁力線和磁感應(yīng)強度分布如圖5所示。

a）磁力線

b）磁感應(yīng)強度

圖5 磁力線和磁感應(yīng)強度分布

采用相同的計算方法，對工作氣隙為0.75 mm時不同源電壓下的電磁吸力和對源電壓為24 V時不同工作氣隙下的電磁吸力分別計算，得到電磁鐵吸力隨源電壓的變化曲線如圖6所示，電磁鐵吸力隨工作氣隙的變化曲線如圖7所示。

圖6 吸力隨源電壓的變化

圖7 吸力隨工作氣隙的變化

2.2 系統(tǒng)參數(shù)仿真

AMESim是基于鍵合圖的液壓/機械系統(tǒng)建模、仿真及動力學(xué)分析軟件。本文使用AMESim軟件對電磁鐵的性能進行仿真計算[13~15]。利用軟件提供的相關(guān)庫中的元件搭建電磁鐵的仿真模型，具體仿真模型如圖8所示。根據(jù)設(shè)計方案為各元件賦予不同導(dǎo)磁材料屬性，并按照螺管電磁鐵的結(jié)構(gòu)尺寸進行參數(shù)設(shè)置，主要參數(shù)設(shè)置如表1所示。

圖8 AMESim仿真模型

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

將源電壓設(shè)置為階躍輸入進行系統(tǒng)仿真計算，得到線圈電流隨時間變化的曲線，計算結(jié)果如圖9所示。工作氣隙為0.75 mm時不同源電壓下的電磁吸力如圖6所示，源電壓為24 V時不同工作氣隙下的電磁吸力如圖7所示。

圖9 電流隨時間的變化

3 吸力試驗

將6件電磁鐵產(chǎn)品與試驗工裝按照如圖10所示測力裝置進行組裝，用吊重法進行吸力試驗。調(diào)節(jié)銜鐵與擋鐵的氣隙，依次增加砝碼，給線圈通電24 VDC，銜鐵運動時的砝碼質(zhì)量即為電磁力值。試驗數(shù)據(jù)與仿真結(jié)果對比如圖7所示。

圖10 測力裝置

4 計算結(jié)果分析與對比

分別采用有限元軟件ANSYS、系統(tǒng)仿真軟件AMESim和吸力試驗3種方法對電磁鐵的吸力特性進行分析和對比研究。

根據(jù)ANSYS計算結(jié)果可知，當(dāng)電磁鐵工作氣隙為0.75 mm、源電壓為24 V時，電磁吸力為117.5 N，氣隙磁感應(yīng)強度約為1.1 T。由AMESim計算結(jié)果可知，當(dāng)電磁鐵工作氣隙為0.75 mm、源電壓為24 V時，電磁吸力為112 N；當(dāng)源電壓為階躍輸入時，線圈電流由0 A升至峰值需要約0.2 s，其中吸合過程（上升段）和斷開過程（下降段）出現(xiàn)局部波峰，其原因為銜鐵運動產(chǎn)生與線圈電流方向相反的電動勢，從而導(dǎo)致電流短暫下降（吸合）或上升（斷開），此現(xiàn)象為螺管電磁鐵的典型特征，如圖9所示。

針對同一工作氣隙不同源電壓下的電磁吸力，采用ANSYS計算和AMESim計算方法得到的電磁吸力隨源電壓的變化趨勢相同（見圖6）。其中，源電壓小于24 V時的電磁吸力ANSYS的計算結(jié)果略大于AMESim的計算結(jié)果，當(dāng)源電壓大于24 V時的電磁吸力ANSYS的計算值較AMESim的計算值偏低。其原因為ANSYS計算可以根據(jù)電磁鐵結(jié)構(gòu)模型將幾何分布考慮在內(nèi)，對于氣隙處漏磁情況的計算相對準確，而AMESim仿真模型中將漏磁系數(shù)設(shè)置為常數(shù)，在較高的磁勢下AMESim仿真模型的計算準確性低于ANSYS模型。將仿真結(jié)果與吸力試驗結(jié)果對比可以看出，在額定磁勢條件下，仿真計算的吸力值結(jié)果與實測的吸力值基本一致。在氣隙較小的情況下，仿真計算值略大于試驗值；在氣隙較大的情況下，仿真計算值略小于試驗值。

工作氣隙為0.75 mm時電磁鐵電磁吸力的設(shè)計值為121 N，ANSYS仿真結(jié)果為117.5 N，比設(shè)計值低2.9%；AMESim仿真結(jié)果為112 N，比設(shè)計值低7.4%；6件電磁鐵產(chǎn)品電磁吸力實測值也均小于設(shè)計值。根據(jù)圖5中的磁力線分布可以看出，電磁鐵內(nèi)存在漏磁現(xiàn)象，漏磁將導(dǎo)致磁感應(yīng)強度降低，從而導(dǎo)致電磁吸力減小。電磁鐵結(jié)構(gòu)形式及結(jié)構(gòu)參數(shù)對電磁鐵的漏磁現(xiàn)象產(chǎn)生影響，而工程計算中對漏磁現(xiàn)象的計算并不十分精確，因而仿真計算得到的電磁吸力與設(shè)計值存在偏差。另外，工程計算以及數(shù)值模擬仿真對于氣體阻尼、鍍膜、銜鐵側(cè)向吸力等因素產(chǎn)生的摩擦力將無法準確計算，一般根據(jù)經(jīng)驗對相關(guān)因素的影響選取一定值的系數(shù)予以考慮，這些因素的影響將導(dǎo)致實際產(chǎn)品電磁吸力實測值與設(shè)計值、仿真值存在偏差。因此，電磁鐵進行工程設(shè)計中應(yīng)對磁感應(yīng)強度留有適當(dāng)?shù)挠嗔?，從而保證電磁吸力滿足使用要求。

通過上述對比分析，采用ANSYS和AMESim軟件均可以實現(xiàn)對螺管電磁鐵吸力特性的仿真。AMESim的仿真結(jié)果與設(shè)計值的偏差略大于ANSYS仿真結(jié)果與設(shè)計值的偏差，但采用AMESim對電磁鐵吸力特性進行仿真不需要對電磁鐵進行詳細的三維實體建模、網(wǎng)格劃分等工作，只需要利用AMESim元件庫中已有元件對電磁鐵的仿真模型進行構(gòu)建，并根據(jù)電磁鐵的結(jié)構(gòu)參數(shù)對子模型的參數(shù)進行設(shè)置，可快速完成仿真計算。因此，采用AMESim軟件對電磁鐵進行仿真分析的優(yōu)勢在于無需精確三維建模即可實現(xiàn)對電磁鐵性能的快速評估，在方案設(shè)計初始階段對電磁鐵結(jié)構(gòu)參數(shù)的確定和優(yōu)化可以實現(xiàn)快速迭代，大幅縮短設(shè)計周期。

5 結(jié) 論

通過對ANSYS、AMESim仿真結(jié)果與吸力試驗結(jié)果進行對比與分析，可以得到以下結(jié)論：

a）在額定磁勢條件下，ANSYS、AMESim仿真計算結(jié)果與試驗實測的吸力值基本一致；氣隙較小時，仿真計算結(jié)果略大于試驗實測值；氣隙較大時，仿真計算結(jié)果略小于試驗實測值。

b）針對同一氣隙不同源電壓下的電磁吸力，ANSYS和AMESim仿真得到的電磁吸力隨源電壓的變化趨勢相同。源電壓小于24 V時，ANSYS仿真結(jié)果略大于AMESim仿真結(jié)果；源電壓大于24 V時，ANSYS仿真結(jié)果略小于AMESim仿真結(jié)果。

c）ANSYS、AMESim仿真結(jié)果與吸力試驗實測值均小于電磁吸力的設(shè)計值，其原因為電磁鐵內(nèi)存在漏磁現(xiàn)象導(dǎo)致磁感應(yīng)強度有所降低，從而導(dǎo)致電磁吸力偏小。

d）采用ANSYS和AMESim軟件均可以實現(xiàn)對電磁鐵吸力特性進行仿真，但AMESim軟件的優(yōu)勢在于無需精確三維建模即可實現(xiàn)對電磁鐵性能的快速評估，在方案設(shè)計初始階段對電磁鐵結(jié)構(gòu)參數(shù)的確定和優(yōu)化可以實現(xiàn)快速迭代，大幅縮短設(shè)計周期。

[1] 任志彬, 王宗偉, 常志鵬, 等. 電磁閥用自鎖電磁鐵吸力與電壓關(guān)系研究[J]. 推進技術(shù), 2016, 37(7): 1372-1379.

[2] 朱寧昌, 劉國球, 等. 液體火箭發(fā)動機設(shè)計(下)[M]. 北京: 宇航出版社, 1994.

[3] 楊亦婷, 李進賢, 黃啟龍. 軌∕姿控發(fā)動機用電磁閥動態(tài)特性仿真[J]. 計算機仿真, 2015, 32(9): 48-52.

[4] 張宏, 廉自生, 熊曉燕, 等. 先導(dǎo)閥電磁鐵的吸力特性數(shù)值模擬與試驗[J]. 煤炭學(xué)報, 2008, 33(3): 339-342.

[5] 龐末紅, 楊倫奎, 陳成峰. 基于Maxwell 方程的電磁閥特性仿真研究[J]. 導(dǎo)彈與航天運載技術(shù), 2015(6): 78-83.

[6] 婁路亮, 王海洲. 電磁閥設(shè)計中電磁力的工程計算方法[J]. 導(dǎo)彈與航運載技術(shù), 2007(1): 40-45.

[7] 張冠生, 陸儉國. 電磁鐵與自動電磁元件[M]. 北京: 機械工業(yè)出版社, 1982.

[8] 葛偉亮. 電磁元件及其控制[M]. 北京: 機械工業(yè)出版社, 1998.

[9] 孫雨施, 王素菊等. 直流磁系統(tǒng)的計算與分析, 模型?算法?程序[M]. 北京: 國防工業(yè)出版社, 1987.

[10] Wang S M, Miyano T, Hubbard M. Electromagnetic field analysis and dynamic simulation of a two-valve solenoid actuator[J]. IEEE Transactions on Magnetics, 1993, 29(2): 1741-1746.

[11] Lee G S, Sung H J, Kim H C. Multiphysics analysis of a linear control solenoid valve[J]. Journal of Fluids Engineering, 2013, 135: 011104-1-011104-10.

[12] 申衛(wèi)兵, 周靜, 李永振. 基于ANSYS的電磁閥仿真與分析[J]. 流體傳動與控制, 2010(5): 4-5.

[13] 龐末紅, 楊倫奎, 陳成峰. 基于AMESim的電磁閥動態(tài)響應(yīng)特性仿真研究[J]. 導(dǎo)彈與航運載技術(shù), 2015(5): 83-88.

[14] 袁洪濱, 張民慶, 孫彥堂. 基于AMESim的直動式電磁閥動態(tài)仿真研究[J]. 火箭推進, 2011(5): 30-35.

[15] 陳亮, 梁國柱. 基于AMESim的電磁閥工作過程動態(tài)特性建模與仿真[J]. 導(dǎo)彈與航天運載技術(shù), 2014(3): 49-54.

[16] 蔡偉, 張碩, 戴民強, 等. 基于AMEsim的電磁閥仿真與試驗驗證[J]. 科學(xué)技術(shù)與工程, 2016, 16(28): 80-84.

Simulation of the Solenoid Force Based on ANSYS and AMESim

Wang Zong-wei, Ren Zhi-bin, Chang Zhi-peng, Liu Jian-she, Cui Peng-fei

(Science and Technology on Space Physics Laboratory, Beijing, 100076)

The electromagnetic field distribution, the magnetic induction distribution and the solenoid force characteristic of a normal opened solenoid magnet is simulated by ANSYS and AMESim. The change curve of solenoid force with voltage and the change curve of solenoid force with working air gap is concluded. And then the results between simulation and solenoid force test is compared and analyzed. The results show that the simulation results of ANSYS and AMESim are basically consistent with the experiment. The simulation results of the solenoid force can be realized. But AMESim can complete the simulation quickly and then the design cycle can be shortened.

Solenoid electromagnet; Solenoid force; Simulation analysis

1004-7182(2017)06-0093-05

10.7654/j.issn.1004-7182.20170620

V43

A

2017-08-23；

2017-10-27

王宗偉（1985-），男，工程師，主要研究方向為飛行器動力系統(tǒng)設(shè)計