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

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

?

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

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

（空間物理重點(diǎn)實(shí)驗(yàn)室，北京，100076）

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

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

0 引 言

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

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

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

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

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

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

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

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

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

2 電磁鐵的仿真計(jì)算

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

2.1 有限元軟件仿真

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

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

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

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

a）電工純鐵DT4

b）軟磁合金1J50

圖2 材料的曲線

d）按照設(shè)計(jì)方案的結(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é)點(diǎn)耦合電流自由度，并對單元施加電壓激勵。

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

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

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

a）磁力線

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

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

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

圖6 吸力隨源電壓的變化

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

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

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

圖8 AMESim仿真模型

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

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

圖9 電流隨時間的變化

3 吸力試驗(yàn)

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

圖10 測力裝置

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

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

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

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

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

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

5 結(jié) 論

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

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

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

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

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

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

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

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

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

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

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

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

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

[9] 孫雨施, 王素菊等. 直流磁系統(tǒng)的計(jì)算與分析, 模型?算法?程序[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)彈與航運(yùn)載技術(shù), 2015(5): 83-88.

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

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

[16] 蔡偉, 張碩, 戴民強(qiáng), 等. 基于AMEsim的電磁閥仿真與試驗(yàn)驗(yàn)證[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-），男，工程師，主要研究方向?yàn)轱w行器動力系統(tǒng)設(shè)計(jì)