微型冷氣推力器電磁驅(qū)動方案優(yōu)化設(shè)計

張 曦 郭加利 張鑫彬 王 玥

1.上海航天控制技術(shù)研究所,上海,2011092.上海伺服系統(tǒng)工程技術(shù)研究中心，上海，201109

微型冷氣推力器電磁驅(qū)動方案優(yōu)化設(shè)計

張 曦1,2郭加利1張鑫彬1,2王 玥1

1.上海航天控制技術(shù)研究所,上海,2011092.上海伺服系統(tǒng)工程技術(shù)研究中心，上海，201109

對螺線管驅(qū)動式微型冷氣推力器電磁驅(qū)動方案相關(guān)設(shè)計指標進行分析和梳理，確定電磁驅(qū)動線圈的評價指標。將CAE分析與綜合評價方法相結(jié)合，利用經(jīng)驗證有效的微型冷氣推力器電磁驅(qū)動仿真模型，完成符合指標要求的驅(qū)動方案集設(shè)計，進而采用組合賦權(quán)法對方案集中各方案進行綜合評價，并根據(jù)評價結(jié)果指導(dǎo)樣機改進，樣機測試結(jié)果與設(shè)計結(jié)果一致性良好。

微型冷氣推力器；電磁驅(qū)動；動態(tài)響應(yīng)；方案優(yōu)化

0 引言

現(xiàn)代小衛(wèi)星技術(shù)的快速發(fā)展,特別是小衛(wèi)星編隊飛行和小衛(wèi)星軌道機動的要求,使得微推進系統(tǒng)和微推進器的技術(shù)發(fā)展越來越受到關(guān)注[1]。微型冷氣推進系統(tǒng)具有結(jié)構(gòu)簡單、功耗低、成本低、安全性和可靠性高等優(yōu)點，不僅適用于高性能微小衛(wèi)星，而且是其他高要求空間任務(wù)的備選方案[2]。英國薩利衛(wèi)星技術(shù)有限公司開發(fā)的微衛(wèi)星和小衛(wèi)星平臺都采用了冷氣微推進系統(tǒng)[3]；日本PROCYON探測器采用冷氣微推進技術(shù)進行方向控制和軌跡修正[4]；意大利Thales Alenia Space公司的高精度冷氣微推進系統(tǒng)被應(yīng)用于GAIA衛(wèi)星[5]；2008年10月，閃蒸射流式冷氣推進技術(shù)被成功應(yīng)用于國內(nèi)某伴飛衛(wèi)星[6]。

微推力器是微型冷氣推進系統(tǒng)的核心元件，其性能的優(yōu)化是提升系統(tǒng)整體性能的基礎(chǔ)。螺線管驅(qū)動式微推力器技術(shù)成熟度高，是微小衛(wèi)星推進系統(tǒng)重要備選產(chǎn)品之一[7]，其電磁驅(qū)動方案設(shè)計結(jié)果直接影響產(chǎn)品的功耗和響應(yīng)時間，對提升工作條件(溫度、壓力等)變化時的工作可靠性，減小產(chǎn)品外包絡(luò)尺寸和質(zhì)量亦具有一定作用。本文對螺線管驅(qū)動式微型冷氣推力器電磁驅(qū)動方案相關(guān)設(shè)計指標進行分析和梳理，確定其評價指標，建立微推力器電磁驅(qū)動部分的動態(tài)仿真模型，對驅(qū)動方案進行設(shè)計，進而采用組合賦權(quán)法完成方案評價和優(yōu)化設(shè)計。

1 微型冷氣推力器結(jié)構(gòu)及工作原理

微型冷氣推力器主要由殼體(骨架)、線圈、彈簧、外殼、銜鐵組件、噴管等組成。與較大規(guī)格冷氣推力器不同，為減小產(chǎn)品體積和質(zhì)量，微型冷氣推力器的噴管與電磁閥通常為一體式，其噴管入口即為電磁閥閥口，通過驅(qū)動帶密封塊的銜鐵組件控制推力器工作，其典型結(jié)構(gòu)如圖1[8]所示。

圖1 微型冷氣推力器典型結(jié)構(gòu)[8]Fig.1 Typical structure of a cold gas micro-thruster[8]

未通電時，銜鐵組件在彈簧和氣壓的作用下，關(guān)閉噴管的進氣口，使噴管進氣截止；線圈通電后，產(chǎn)生的電磁力克服彈簧和介質(zhì)壓力，驅(qū)動銜鐵組件運動，閥門打開，氣體由噴管噴出產(chǎn)生推力。

2 電磁驅(qū)動方案評價指標

在產(chǎn)品結(jié)構(gòu)和材料設(shè)計方案基本確定時，電磁驅(qū)動方案設(shè)計主要影響微推力器的響應(yīng)時間、功耗、最低工作電壓和線圈厚度等指標。

微推力器響應(yīng)時間包括開啟響應(yīng)時間和關(guān)閉響應(yīng)時間，開啟響應(yīng)時間為給推力器通電至推力增大到90%額定推力的時間，關(guān)閉響應(yīng)時間為給推力器斷電至推力減小到10%額定推力的時間[9]。微推力器推力水平較低，測試難度較大、費用高，而微推力器相對于一般電磁閥僅在出口處存在較大差異，故在推力噴管及閥口結(jié)構(gòu)一定時，空載響應(yīng)時間主要由電磁驅(qū)動設(shè)計方案決定，并直接影響推力器響應(yīng)時間指標的優(yōu)劣。由于螺線管電磁閥關(guān)閉響應(yīng)時間受放電回路影響較大[10-11]，故本文以空載開啟響應(yīng)時間作為微推力器電磁驅(qū)動方案評價的指標之一。

對于微推力器，在結(jié)構(gòu)與材料相同的情況下，一般來說，功耗越高響應(yīng)性能越好；然而，功耗亦為產(chǎn)品的重要指標和設(shè)計約束，功耗與響應(yīng)性能的權(quán)衡和取舍是電磁驅(qū)動方案評價的重要組成部分。

最低工作電壓是各類電磁閥的重要使用指標之一。螺線管線圈阻值會隨工作溫度的變化而改變，從而導(dǎo)致閥的啟動性能變化；溫度變化導(dǎo)致密封材料的熱脹冷縮和變形，亦會在一定程度上影響推力器的啟動性能；推力器在使用過程中，銜鐵組件中的密封材料發(fā)生的塑性變形導(dǎo)致的工作氣隙變化也會使推力器開啟性能劣化，因此，在產(chǎn)品設(shè)計時必須留有足夠的裕量，以保證推力器在規(guī)定的環(huán)境條件下和使用壽命內(nèi)能正常工作。

對于相同殼體(骨架)，不同的線圈設(shè)計方案下，線圈厚度不同。在外殼尺寸(外包絡(luò))一定的情況下，線圈薄便于實施灌封等后續(xù)工序；在灌封工藝需求一定的情況下，線圈薄可使產(chǎn)品外殼尺寸相應(yīng)減小，對于螺線管結(jié)構(gòu)，外包絡(luò)徑向尺寸的減小能有效優(yōu)化產(chǎn)品的質(zhì)量指標。

綜上，本文以空載條件下開啟響應(yīng)時間T0、功耗P、最低工作電壓U0和線圈厚度H作為評價指標，建立電磁驅(qū)動方案的評價體系。

3 電磁驅(qū)動方案設(shè)計

3.1 線圈計算

對于給定的殼體結(jié)構(gòu)，驅(qū)動線圈的厚度

(1)

式中，n為層數(shù)；k為線圈層厚系數(shù)；d為漆包線外徑。

線長可由下式進行近似計算：

(2)

式中，l為殼體有效繞線長度；D為殼體內(nèi)徑。

線圈的電阻可根據(jù)漆包線每米電阻值[12]與線長相乘得到，進而得到額定工作電壓時推力器的功耗值。

3.2 動態(tài)性能建模與驗證

采用仿真技術(shù)對電磁閥進行動態(tài)建模分析，與以往設(shè)計電磁閥利用經(jīng)驗公式與解析計算確定參數(shù),再通過樣機試驗進行驗證和修正的傳統(tǒng)方式相比，可有效縮短產(chǎn)品的研制周期、降低成本和提升效率[13-15]。本文利用AMESim仿真平臺建立微推力器電磁驅(qū)動部分的動態(tài)仿真模型。

動態(tài)性能仿真模型由勵磁電路、磁路及機械運動部分組成，如圖2所示。

圖2 電磁驅(qū)動性能動態(tài)仿真模型圖Fig.2 Dynamic simulation model for the electromagnetic drive mechanism

對于勵磁電路，有

(3)

式中，U為輸入電壓；i為線圈電流；R為線圈回路總電阻；ψ為線圈磁鏈。

對于驅(qū)動磁路，有

Ni=∑λi

(4)

式中，λi為磁路各部分磁動勢；N為線圈匝數(shù)。

根據(jù)微推力器磁路結(jié)構(gòu)，將磁路模型劃分為工作氣隙磁阻Rg0、非工作氣隙磁阻Rg1、磁路磁阻Rm以及漏磁。其中，磁路磁阻由磁路沿程軟磁材料導(dǎo)磁結(jié)構(gòu)構(gòu)成，包括徑向磁阻(R2、R4～R7)和軸向磁阻(R1、R3)兩類。

對于磁阻，有

R=lm/(μA)

(5)

式中，R為磁阻；lm為磁阻特征長度；μ為磁導(dǎo)率；A為磁阻特征面積。

通過磁阻的磁通

Φ=λ/R

(6)

式中，Φ為磁通；λ為作用在磁阻上的磁動勢。

工作氣隙處電磁力用于驅(qū)動銜鐵組件動作，根據(jù)麥克斯韋電磁力公式，有

(7)

式中，F(xiàn)m為銜鐵受到的電磁力；Φg0為工作氣隙處的磁通；μ0為真空磁導(dǎo)率；Ag0為銜鐵底面面積。

對銜鐵組件，有

(8)

式中，m為銜鐵組件的質(zhì)量；x為氣隙長度；Fs為彈簧力；Ff為摩擦力。

銜鐵組件的動作較難直接檢測，可通過線圈電流信號進行間接檢測，其檢測原理如下：由式(6)可知，作用于工作氣隙的磁動勢

λg0=Φg0Rg0

(9)

對于活動的工作氣隙，由式(5)和式(9)有

(10)

圖3所示為樣機在額定工作電壓12V時開啟響應(yīng)時間測試結(jié)果與模型仿真結(jié)果對比。開啟響應(yīng)時間的仿真結(jié)果為3.68ms，樣機測試結(jié)果為3.78ms。通過調(diào)整線圈輸入電壓對推力器最低開啟電壓進行仿真計算，推力器空載開啟電壓為8.74V，實測值為9V(樣機額定工作壓力條件下推力器開啟電壓較空載開啟電壓高約0.5V，本文以空載條件下開啟電壓作為評價指標)。仿真模型與樣機實測結(jié)果一致性良好，能較真實地反映樣機性能。

圖3 樣機響應(yīng)性能測試結(jié)果與仿真結(jié)果Fig.3 Test data and simulation result of the prototype

3.3 方案設(shè)計

利用試制樣機的測試數(shù)據(jù)校驗?zāi)Ｐ蜏蚀_性后，結(jié)合驅(qū)動方案基本參數(shù)對待定方案進行設(shè)計。取有效方案各指標閾值如下：開啟響應(yīng)時間T0≤3.5 ms，功耗P≤3 W，最低開啟電壓U0≤8.5 V，線圈厚度H≤1.8 mm，通過計算和仿真得到符合要求的9套待定方案，見表1。

表1 驅(qū)動方案Tab.1 Driving schemes

對表1中不同方案各指標進行初步計算可知，在不改變推力器機械結(jié)構(gòu)的條件下，僅改進驅(qū)動線圈的方案，符合設(shè)計指標要求的不同驅(qū)動方案的各個指標之間差異明顯，其中，推力器空載開啟響應(yīng)時間差異最大可達25.8%，功耗差異可達40.4%，最低啟動電壓差異可達32.1%，線圈厚度差異可達44%。由方案設(shè)計結(jié)果可以明顯看出，功耗、開啟響應(yīng)、最低啟動電壓、線圈厚度等指標間存在制約關(guān)系，故有必要對不同方案進行合理評價。

以表1中方案及對應(yīng)的各個指標值為基礎(chǔ)建立分析矩陣，其中，矩陣元素xij對應(yīng)第i個方案中第j個指標，并記xjmax和xjmin分別為各指標列中最大與最小的指標值。不同指標的量綱不同，根據(jù)指標評價特征，對矩陣進行量綱一化，令

可得決策矩陣

決策矩陣Y中第1至第9行分別對應(yīng)表1中方案A1至A9，Y中第1至第4列分別對應(yīng)T0、P、U0、H四個評價指標。Y中各元素值對應(yīng)指標的優(yōu)劣，1為最優(yōu)，0為最差。

4 驅(qū)動方案評價

微推力器驅(qū)動方案的評價是已知方案進行多指標決策的問題，常用賦權(quán)方法有主觀賦權(quán)法、客觀賦權(quán)法和組合賦權(quán)法等。其中，主觀賦權(quán)法客觀性較差，客觀賦權(quán)法確定的權(quán)值有時與指標的實際重要程度相悖，而組合賦權(quán)法將主觀和客觀賦權(quán)法相結(jié)合，避免了單一方法在賦權(quán)中存在的局限性[16-17]。本文采用線性加權(quán)法對方案進行組合賦權(quán)與評價。組合權(quán)重可表示為

wj=αaj+(1-α)bj

(11)

式中，wj為第j個指標的最終權(quán)重；α為權(quán)系數(shù)；aj為第j個指標的主觀權(quán)重；bj為第j個指標的客觀權(quán)重。

各方案的綜合評價值

di=wj∑yij

(12)

根據(jù)各指標設(shè)計重要性，取指標主觀權(quán)重矩陣W′=[0.4 0.3 0.15 0.15]?？陀^權(quán)重采用均方差法計算，有W″=[0.246 0.242 0.260 0.252 ]。取α=0.5，可得組合權(quán)重矩陣W=[0.323 0.271 0.205 0.201 ]。

根據(jù)式(12)對決策矩陣進行計算，可得各方案綜合評價矩陣D=[0.376 0.455 0.525 0.585 0.528 0.563 0.578 0.504 0.415]T，即方案綜合評價排序結(jié)果由優(yōu)到劣依次為：A4，A7，A6，A5，A3，A8，A9，A1。

采用方案A4進行樣機制作，推力器線圈厚度約1.35 mm，功耗為2.38 W，響應(yīng)時間為2.85 ms，最低開啟電壓為7.3 V，測試結(jié)果與設(shè)計結(jié)果一致性良好。

5 結(jié)論

(1)根據(jù)微推力器電磁驅(qū)動相關(guān)指標，結(jié)合測試便利性，對驅(qū)動方案評價指標進行分析，確定通過空載開啟響應(yīng)時間、功耗、最低啟動電壓和線圈厚度等指標對電磁驅(qū)動方案進行評價。

(2)建立了微推力器電磁驅(qū)動部分的動態(tài)仿真模型，仿真結(jié)果與試制樣機測試結(jié)果一致性良好。

(3)對驅(qū)動方案進行設(shè)計，結(jié)果表明，對于特定的推力器結(jié)構(gòu)，符合設(shè)計要求的不同驅(qū)動線圈方案，各主要指標差異明顯,且各指標相互之間存在一定的制約。

(4)采用組合賦權(quán)法對各方案進行評價，并根據(jù)評價結(jié)果進行樣機制作，樣機測試結(jié)果符合設(shè)計預(yù)期要求。

(5)微型冷氣推力器電磁驅(qū)動部分采用的直動式電磁閥結(jié)構(gòu)在航天及工業(yè)領(lǐng)域有廣泛的應(yīng)用，本文設(shè)計與評價方法對類似產(chǎn)品的設(shè)計具有一定參考意義。

[1] 林來興.現(xiàn)代小衛(wèi)星的微推進系統(tǒng)[J].航天器工程，2010，19(6)：13-20. LIN Laixing. Micro-propulsion System for Modern Small Satellites[J]. Spacecraft Engineering，2010，19(6)：13-20.

[2] NGUYEN H，K?HLER J, STENMARK L.The Merits of Cold Gas Micropropulsion in State-of-the-art Space Mission[R]. Houston:American Institute of Aeronautics and Astronautics，2002.

[3] 吳漢基，蔣遠大，張志遠，等.微小衛(wèi)星的在軌推進技術(shù)[J].火箭推進，2006，32(3)：40-43. WU Hanji，JIANG Yuanda，ZHANG Zhiyuan，et al.On-board Propulsion Technologies for Micro/Minisatellites[J]. Journal of Rocket Propulsion，2006，32(3)：40-43.

[4] KAWAHARA H.Ground Experiment for the Small Unified Propulsion System：I-COUPS Installed on the small space probe：PROCYON [C]//Joint Conference of 30th ISTS，34th IEPC and 6th NSAT.Kobe-Hyogo，2015：1-13.

[5] MATTICARI G，NOCI G，SICILIANO P，et al.Cold Gas Micro Propulsion Prototype for Very Fine Spacecraft Attitude/Position Control[R]. Sacramento:American Institute of Aeronautics and Astronautics，2006.

[6] 魏青，郭尚群.閃蒸射流推進的應(yīng)用[J].火箭推進，2010，36(3)：19-23. WEI Qing，GUO Shangqun. Application of the Flashing Jet Propulsion[J]. Journal of Rocket Propulsion，2010，36(3)：19-23.

[7] MUELLER J, ZIEMER J, HOFER R，et al.A Survey of Micro Thrust Propulsion Options for Microspacecraft and Formation Flying Missions[C]//5th Annual Cube Sat Developers Workshop. San Luis Obispo，2008：1-19.

[8] SMITH P，EDWARDS S, SOLWAY N.Cryosat Cold Gas System and Component Development[R]. Fort Lauderdale: American Institute of Aeronautics and Astronautics，2004.

[9] BZIBZIAK S R.Miniature Cold Gas Thrusters[R]. Nashville:American Institute of Aeronautics and Astronautics，1992.

[10] 沈公槐.改進電磁閥動態(tài)特性的一種方法[J].航天控制，2001，19(3)：75-81. SHEN Gonghuai.A Method Improving the Dynamic Characteristic of Solenoid Valve[J]. Aerospace Control，2001，19 (3)：75-81.

[11] 姜華，陳健，汪旭東.拍合式電磁閥動態(tài)響應(yīng)仿真分析與試驗驗證[J].空間控制技術(shù)與應(yīng)用，2013，39(4)：23-27. JIANG Hua，CHEN Jian，WANG Xudong.Simulation Analysis of Clapper Type Solenoid Valve’s Response Characteristics and Experimental Verification[J]. Aerospace Control and Application，2013，39(4)：23-27.

[12] 全國電線電纜標準化技術(shù)委員會.GB/T 6109.1-2008 漆包圓繞組線第1部分：一般規(guī)定[S].北京：中國標準出版社，2008. SAC/TC 213.GB/T 6109.1-2008 Enamelled Round Winding Wire-Part 1:General Requirements[S].Beiing：Standards Press of China，2008.

[13] 張榛.電磁閥動態(tài)響應(yīng)特性的有限元仿真與優(yōu)化設(shè)計[J].空間控制技術(shù)與應(yīng)用，2008，34(5)：53-56. ZHANG Zhen.FEA Simulation of Dynamic Response of Solenoid Valve and Its Optimal Design[J]. Aerospace Control and Application，2008，34(5)：53-56.

[14] 袁洪濱，張民慶，孫彥堂.基于AMESim的直動式電磁閥動態(tài)仿真研究[J].火箭推進，2011,37(5)：30-35. YUAN Hongbin，ZHANG Minqing，SUN Yantang. AMESim-based Simulation Analysis of Dynamic Characteristics of Direct-acting Solenoid Valve[J]. Journal of Rocket Propulsion，2011，37(5)：30-35.

[15] 龐末紅，楊倫奎，陳成峰.基于 AMESim 的電磁閥動態(tài)響應(yīng)特性仿真研究[J].導(dǎo)彈與航天運載技術(shù)，2015(5)：83-88. PANG Mohong, YANG Lunkui, CHEN Chengfeng. AMESim-based Simulation on the Dynamic Response Characteristics of Solenoid Valve[J]. Missiles and Space Vehicles，2015(5)：83-88.

[16] 王明濤.多指標綜合評價中權(quán)系數(shù)確定的一種綜合分析方法[J].系統(tǒng)工程，1999，17(2)：56-61. WANG Mingtao. A Comprehensive Analysis Method on Determining the Coefficients in Multi-index Evaluation[J]. Systems Engineering，1999，17(2)：56-61.

[17] 宋光興，楊德禮.基于決策者偏好及賦權(quán)法一致性的組合賦權(quán)法[J]. 系統(tǒng)工程與電子技術(shù)，2004，26(9)：1226-1230. SONG Guangxing，YANG Deli. Combination Weighting Approach Based on the Decision Maker’s Preference and Consistency of Weighting Methods[J]. Systems Engineering and Electronics，2004，26(9)：1226-1230.

(編輯 陳 勇)

Optimization Design of Electromagnetic-driving Schemes of a Cold Gas Micro-thruster

ZHANG Xi1,2GUO Jiali1ZHANG Xinbin1,2WANG Yue1

1.Shanghai Institute of Spaceflight Control Technology, Shanghai,201109 2.Shanghai Engineering Research Center of Servo Systems, Shanghai,201109

Evaluation indicators of electromagnetic-driving project for solenoid-based cold gas micro-thruster were discussed. CAE analysis and comprehensive evaluation method were used to optimize the solenoid design for a cold gas micro-thruster. At first, a dynamic simulation model which could match the test data well was built to establish the solenoid design matrix. Then, the optimal design was acquired by using the combination assigning method. The test results of the prototype by using the excellent scheme are proved to be accord with the expectations.

cold gas micro-thruster; electromagnetic-driving; dynamic response; scheme optimization

2016-10-17

TH138.52

10.3969/j.issn.1004-132X.2017.09.005

張 曦，男，1982年生。上海航天控制技術(shù)研究所工程師。主要研究方向為電液伺服系統(tǒng)、微推進系統(tǒng)組件技術(shù)。E-mail：messial@163.com。郭加利，男，1987年生。上海航天控制技術(shù)研究所工程師。張鑫彬，男，1980年生。上海航天控制技術(shù)研究所高級工程師。王 玥，女，1983年生。上海航天控制技術(shù)研究所工程師。