張 烽，申 麟，吳勝寶，王小錠，王書廷

（中國(guó)運(yùn)載火箭技術(shù)研究院研究發(fā)展中心，北京，100076）

電動(dòng)力繩系離軌技術(shù)性能與任務(wù)適應(yīng)性分析

張 烽，申 麟，吳勝寶，王小錠，王書廷

（中國(guó)運(yùn)載火箭技術(shù)研究院研究發(fā)展中心，北京，100076）

電動(dòng)力繩系離軌技術(shù)屬于被動(dòng)離軌方式，適用于火箭末級(jí)、失效衛(wèi)星等廢棄飛行器快速離軌。針對(duì)這項(xiàng)新穎的離軌技術(shù)，通過建立多場(chǎng)耦合攝動(dòng)下的離軌動(dòng)力學(xué)模型，研究并分析了該技術(shù)的離軌性能和適用的離軌任務(wù)類型，用于指導(dǎo)工程實(shí)踐及控制方案設(shè)計(jì)。首先建立多場(chǎng)耦合攝動(dòng)下的電動(dòng)力繩系離軌動(dòng)力學(xué)模型，基于數(shù)值仿真，驗(yàn)證該項(xiàng)技術(shù)的有效性及優(yōu)越性，并重點(diǎn)研究了不同影響因素下的離軌性能，據(jù)此提出了系統(tǒng)參數(shù)的設(shè)置原則。再通過設(shè)置合理的系統(tǒng)參數(shù)，分析電動(dòng)力繩系離軌技術(shù)的可適用任務(wù)類型。分析結(jié)果表明，電動(dòng)力繩系離軌技術(shù)適用于非極軌道的廢棄飛行器離軌任務(wù)。

電動(dòng)力繩系離軌技術(shù)；動(dòng)力學(xué)模型；性能分析；任務(wù)適應(yīng)性

0 引 言

隨著人類空間活動(dòng)日益頻繁，在宇宙中遺留大量的太空垃圾不可避免，空間碎片清理已成為世界航天界關(guān)注的重大課題[1～3]。為實(shí)現(xiàn)空間碎片清除，目前主要有推力離軌方法、網(wǎng)狀機(jī)構(gòu)捕獲方法[1]、空間機(jī)器人方法[2]等，但是這些方法都存在耗能大、成本高、結(jié)構(gòu)復(fù)雜等缺點(diǎn)。

近年來，隨著空間繩系技術(shù)的發(fā)展，基于電動(dòng)力繩系的飛行器離軌技術(shù)成為了一種新穎的被動(dòng)離軌方式，特別適用于火箭末級(jí)、廢棄衛(wèi)星等飛行器[4]。其基本原理是：導(dǎo)電繩系在地球磁場(chǎng)中運(yùn)動(dòng)時(shí)，會(huì)由于切割磁力線而在繩系中產(chǎn)生電流，電流和地球磁場(chǎng)相互作用會(huì)在繩系上產(chǎn)生洛侖茲力，形成離軌阻力，實(shí)現(xiàn)快速離軌[5～8]。

與傳統(tǒng)的使用推進(jìn)劑進(jìn)行離軌的方式相比，電動(dòng)力繩系離軌技術(shù)不需要攜帶推進(jìn)劑，可以極大減少離軌系統(tǒng)的質(zhì)量，顯著節(jié)約發(fā)射成本，具有高效率、低質(zhì)量的優(yōu)點(diǎn)。然而，電動(dòng)力繩系離軌系統(tǒng)進(jìn)行離軌時(shí)，系統(tǒng)處于一個(gè)復(fù)雜的耦合多場(chǎng)狀態(tài)，因此，離軌性能也受到多種因素制約。目前的研究大多集中于繩系的控制方案研究[6,7]、硬件方案分析與設(shè)計(jì)[8,9]、繩系動(dòng)力學(xué)建模[10]等問題，鮮有考慮空間多場(chǎng)耦合因素影響，從任務(wù)規(guī)劃級(jí)層面開展電動(dòng)力繩系離軌技術(shù)的離軌性能分析及任務(wù)適應(yīng)性分析，以指導(dǎo)工程設(shè)計(jì)。

本文綜合考慮電磁場(chǎng)、大氣阻力及地球引力場(chǎng)等多種因素，建立了多場(chǎng)耦合攝動(dòng)下的電動(dòng)力繩系離軌動(dòng)力學(xué)模型。隨后，基于數(shù)值仿真，重點(diǎn)研究并分析了不同影響因素下的離軌性能，據(jù)此提出了系統(tǒng)參數(shù)的設(shè)置原則。再通過設(shè)置合理的系統(tǒng)參數(shù)，分析了電動(dòng)力繩系離軌技術(shù)的可適用任務(wù)類型。

1 電動(dòng)力繩系離軌系統(tǒng)動(dòng)力學(xué)

1.1 離軌動(dòng)力學(xué)

電動(dòng)力繩系離軌系統(tǒng)由待離軌飛行器M，末端載荷m和電動(dòng)力繩構(gòu)成，如圖1所示。為便于動(dòng)力學(xué)建模，本文采用下述3組坐標(biāo)系：

a）地心慣性坐標(biāo)系O-XEYEZE：坐標(biāo)系原點(diǎn)O位于地心，OXE軸位于赤道面內(nèi)并指向春分點(diǎn)，OZE軸與地球自轉(zhuǎn)軸重合，OYE軸由右手定則確定。

b）軌道坐標(biāo)系OS-xyz：坐標(biāo)系原點(diǎn)OS位于系統(tǒng)質(zhì)心，x軸位于軌道面內(nèi)，與地心徑重合，并指向天頂方向，z軸與軌道法向重合，y軸指向系統(tǒng)飛行方向，并由右手定則確定。

考慮到系統(tǒng)離軌過程中，繩系擺動(dòng)過程與軌道變化過程時(shí)間尺度差異較大，本文不考慮繩系擺動(dòng)影響，將飛行器與末端載荷視為質(zhì)點(diǎn)，電動(dòng)繩系視為無質(zhì)量的剛性桿，并始終沿軌道徑向保持穩(wěn)定。

基于慣性坐標(biāo)系O-XEYEZE，為保證軌道推算時(shí)不發(fā)生奇異現(xiàn)象，引入非奇異軌道要素{p,ξ,η, h, k, L}，對(duì)經(jīng)典軌道六要素{a, e, i, Ω, ω, v}做如下變換：

式中 a為軌道半長(zhǎng)軸；e為軌道偏心率；i為軌道傾角；Ω為軌道升交點(diǎn)赤經(jīng)；ω為近地點(diǎn)幅角；v為真近點(diǎn)角。

基于上述非奇異變換，系統(tǒng)的非奇異軌道動(dòng)力學(xué)可描述如下：

式中 S，T，W分別是沿軌道坐標(biāo)系x軸，y軸和z軸的空間多場(chǎng)耦合攝動(dòng)加速度分量，具體分析參見1.2節(jié)，w和s2滿足：

1.2 多場(chǎng)耦合動(dòng)力學(xué)分析

電動(dòng)力繩離軌系統(tǒng)在運(yùn)行過程中的主要軌道攝動(dòng)力來源包括：大氣阻力、地球的不均勻性和扁平率、電動(dòng)力繩切割地磁場(chǎng)產(chǎn)生的洛倫茲力、日月等天體引力、太陽輻射光壓。

由于天體引力及太陽輻射光壓所致攝動(dòng)力在近地軌道上較前3種攝動(dòng)力的影響小。因此，本文僅考慮前3種攝動(dòng)力作用。

a）大氣阻力。

系統(tǒng)運(yùn)行于近地軌道時(shí)，飛行器、末端載荷和電動(dòng)力繩系的表面受到大氣阻力的作用。因此，整個(gè)系統(tǒng)受到的大氣阻力dF可以表示為

式中 Fd,M，F(xiàn)d,m，F(xiàn)d,t分別為飛行器、末端載荷和電動(dòng)力繩受到的大氣阻力。

為簡(jiǎn)化分析，將飛行器、末端載荷和電動(dòng)力繩視為圓柱型，并假設(shè)三者在飛行過程中成線形排列，則，

式中dρ為大氣密度；dC為阻力系數(shù)；Md，md和td分別為飛行器、末端載荷及繩系的直徑；ML，mL和l分別為飛行器、末端載荷及繩系的長(zhǎng)度；α為系統(tǒng)軸向矢量tu與相對(duì)速度rv之間的夾角；rv滿足：

式中 v為離軌系統(tǒng)的質(zhì)心速度；Eω為地球自轉(zhuǎn)角速度；r為系統(tǒng)地心徑矢量。

b）地球扁率。

由于地球扁率的存在，地心引力在系統(tǒng)的飛行軌道上不一致，從而會(huì)引起軌道參數(shù)的變化。假設(shè)地球?yàn)橐粋€(gè)均質(zhì)橢球體，考慮地球扁率引起的J2項(xiàng)攝動(dòng)，那么地球引力場(chǎng)的勢(shì)函數(shù)近似表示為

式中 R為地球在赤道處的半徑；μE為地球引力常數(shù)；J2=1.08 263× 10-3。

進(jìn)而，能夠?qū)С龅厍虮饴室鸬臄z動(dòng)力Fg在軌道坐標(biāo)系下的分量：

式中 M0為系統(tǒng)總質(zhì)量，M0=M+m+mt。

c）電動(dòng)力（洛倫茲力）。

系統(tǒng)在離軌過程中，電動(dòng)力繩切割地磁場(chǎng)產(chǎn)生洛倫茲力進(jìn)而形成離軌阻力，具體描述如下[8～10]：

式中 ut為沿繩系徑向的單位矢量；s為繩系微元ds相對(duì)于系統(tǒng)質(zhì)心的距離；B為地磁場(chǎng)產(chǎn)生的磁場(chǎng)強(qiáng)度矢量；()I s為微元ds的電流強(qiáng)度。

根據(jù)繩系材料不同，電動(dòng)力繩存在2種類型：

a）絕緣繩系。

針對(duì)絕緣繩系情形，繩系中的電流保持不變，則洛倫茲力可以表示為

這種情形下，電流I可根據(jù)硬件性能進(jìn)行設(shè)計(jì)，

b）裸繩系。

裸繩系是通過裸露于空間環(huán)境的繩系自身收集電子，再通過繩系末端的電荷接觸裝置將電子射入空間形成回路，從而可以提高電荷交換效率。但是這種類型的電動(dòng)力繩不具備調(diào)節(jié)電流的能力，其電流大小很大程度上取決于空間中的電子密度及地磁場(chǎng)分布。

裸繩系的電流強(qiáng)度可表示為[8～10]

式中 Vp為動(dòng)生電動(dòng)勢(shì)；Vt為繩系電動(dòng)勢(shì)；me為負(fù)電子質(zhì)量；im為正電子質(zhì)量；n∞為無擾離子密度；電壓差ΔV=Vp-Vt滿足下述方程：

式中 σ為繩系電導(dǎo)率；mE為電場(chǎng)強(qiáng)度；A為橫截面積。

2 離軌性能分析

通過第1節(jié)的動(dòng)力學(xué)建模過程得出，當(dāng)待離軌飛行器參數(shù)確定的情況下，影響電動(dòng)力繩系離軌性能的主要因素有：

a） 任務(wù)參數(shù)：初始軌道高度0H，初始軌道傾角0i；

b）系統(tǒng)參數(shù)：末端載荷質(zhì)量m，繩系長(zhǎng)度l，繩系類型（絕緣繩系或裸繩系）。

本節(jié)基于第1節(jié)建立的多場(chǎng)耦合因素下的系統(tǒng)離軌動(dòng)力學(xué)模型，借助數(shù)值仿真，研究并分析不同系統(tǒng)參數(shù)下的系統(tǒng)離軌性能，并據(jù)此給出合理的系統(tǒng)參數(shù)選擇范圍，為研究電動(dòng)力繩系離軌技術(shù)的任務(wù)適應(yīng)性分析提供依據(jù)。

2.1 基本參數(shù)設(shè)置

設(shè)定待離軌飛行器質(zhì)量為4 000 kg，長(zhǎng)度為11 m，直徑為3.35 m。末端載荷長(zhǎng)度為0.5 m，直徑為0.3 m，繩系直徑為0.001 m。初始軌道參數(shù)設(shè)置為：初始軌道高度H0=700 km，軌道傾角i0=0°，偏心率e0=0，升交點(diǎn)赤經(jīng)Ω0=0°，近地點(diǎn)幅角ω0=0°，真近點(diǎn)角v0=0°。

在計(jì)算大氣阻力引起的攝動(dòng)加速度時(shí)，采用NRLMSISE-00大氣模型計(jì)算大氣密度，阻力系數(shù)Cd=2.2。

在計(jì)算洛倫茲力時(shí)，采用第11代國(guó)際地磁參考場(chǎng)模型（IGRF11）計(jì)算13階地磁場(chǎng)磁場(chǎng)強(qiáng)度，同時(shí)，采用國(guó)際參考電離層模型IRI2007計(jì)算裸繩系模型所需的電子密度參數(shù)。

2.2 不同系統(tǒng)參數(shù)下的離軌性能分析

基于2.1節(jié)的基本參數(shù)設(shè)置，借助數(shù)值仿真，研究不同系統(tǒng)參數(shù)下的系統(tǒng)離軌性能，其中，離軌結(jié)束認(rèn)為軌道遠(yuǎn)地點(diǎn)高度降落至120 km以內(nèi)。

a）末端載荷質(zhì)量。

考慮末端載荷質(zhì)量分別為m =40 kg，100 kg，200 kg，400 kg，其它系統(tǒng)參數(shù)設(shè)置為：l =5 km，I =0.5 A（絕緣繩系）。

圖2為不同末端載荷質(zhì)量下的遠(yuǎn)地點(diǎn)高度變化曲線。由圖2可以看出，末端載荷質(zhì)量從40 kg增大至400 kg，離軌時(shí)間并沒有顯著減少。因此，末端載荷質(zhì)量對(duì)離軌性能的影響較小。需要說明的是，當(dāng)遠(yuǎn)地點(diǎn)高度降低至250 km時(shí)，軌道高度出現(xiàn)陡降。這是由于在軌道高度降低至接近地球時(shí)，大氣阻力的作用將更為明顯。

b）繩系長(zhǎng)度。

考慮繩系長(zhǎng)度分別為l =1 km，5 km，10 km，15 km，20 km，其它系統(tǒng)參數(shù)設(shè)置為：m =40 kg，I =0.5 A（絕緣繩系）。

圖3為不同繩系長(zhǎng)度下的遠(yuǎn)地點(diǎn)高度變化曲線。

由圖3可以看出，繩系長(zhǎng)度為1 km時(shí)，離軌時(shí)間最慢，繩系越長(zhǎng)，離軌時(shí)間越短，但是當(dāng)長(zhǎng)度大于10 km后，離軌時(shí)間的降幅并不大。

c）繩系類型。

繩系類型分別為裸繩系，I =0.01 A（絕緣繩系），I =0.1 A（絕緣繩系），I =0.5 A（絕緣繩系），I =1 A（絕緣繩系），其它系統(tǒng)參數(shù)設(shè)置為：m =40 kg，l=5 km。

圖4為不同繩系類型的遠(yuǎn)地點(diǎn)變化曲線。由圖4可以看出，對(duì)于絕緣繩系，電流越大，離軌時(shí)間越短。相比而言，裸繩系的離軌性能與設(shè)定電流為0.5 A的絕緣繩系的離軌性能相當(dāng)。

根據(jù)圖4曲線的變化趨勢(shì)可推測(cè)，當(dāng)電流設(shè)定為0時(shí)，即在常規(guī)的僅依靠大氣阻力及地球引力的離軌方式下，離軌時(shí)間要遠(yuǎn)遠(yuǎn)大于電動(dòng)力繩系離軌系統(tǒng)的離軌時(shí)間，從而間接驗(yàn)證了電動(dòng)力繩系離軌技術(shù)的有效性和優(yōu)越性。

根據(jù)上述仿真結(jié)果，圖5匯總了不同系統(tǒng)參數(shù)下的離軌時(shí)間變化趨勢(shì)。

由圖5可得以下分析結(jié)論：

a）系統(tǒng)離軌時(shí)間受末端載荷質(zhì)量影響較小。因此，為降低運(yùn)載器的有效載荷質(zhì)量，應(yīng)盡可能小地選取末端載荷質(zhì)量。

b）系統(tǒng)離軌時(shí)間受繩系長(zhǎng)度影響較大。但是，繩系增長(zhǎng)至一定程度時(shí)，離軌時(shí)間降幅變緩，同時(shí)，繩系越長(zhǎng)，在空間中受到破壞的概率也就越大。因此，需要根據(jù)任務(wù)的離軌時(shí)間要求，選取合適的繩系長(zhǎng)度。

c）系統(tǒng)離軌時(shí)間受繩系類型影響較大。一方面，對(duì)于絕緣繩系而言，設(shè)置繩系電流越大，離軌時(shí)間越短，但電流的設(shè)定值受到系統(tǒng)硬件能力的制約。而且，絕緣繩系需要添加額外的電荷收集裝置，硬件結(jié)構(gòu)較為復(fù)雜。另一方面，對(duì)于裸繩系而言，電荷交換效率較高，且無需額外的電荷收集裝置，但是，電流無法調(diào)節(jié)，且大小依賴于空間中的電子密度。因此，需要綜合考慮硬件能力及任務(wù)背景，選取合適的繩系類型。

基于上述分析，電動(dòng)力繩系離軌系統(tǒng)的末端載荷可選取飛行器質(zhì)量的1%，繩系長(zhǎng)度選取為1～10 km，繩系類型可選取為裸繩系或絕緣繩系（I =0.5 A）。

3 任務(wù)適應(yīng)性分析

根據(jù)第2節(jié)的分析結(jié)論，合理選取下述2組系統(tǒng)參數(shù)作為標(biāo)準(zhǔn)工況，進(jìn)一步開展電動(dòng)力繩系離軌技術(shù)的任務(wù)適應(yīng)性分析：

a）工況1：末端載荷質(zhì)量40 kg，繩系長(zhǎng)度5 km，繩系類型為裸繩系；

b）工況2：末端載荷質(zhì)量40 kg，繩系長(zhǎng)度5 km，繩系類型為絕緣繩系（I =0.5 A）。

經(jīng)過數(shù)值仿真，表1統(tǒng)計(jì)了2種工況不同初始軌道參數(shù)下的系統(tǒng)離軌時(shí)間。

表1 不同任務(wù)參數(shù)下的系統(tǒng)離軌時(shí)間

由表1可以看出：a）離軌任務(wù)的軌道高度影響較大，初始軌道高度越高，離軌時(shí)間越長(zhǎng)。b）離軌任務(wù)的初始軌道傾角影響較大，當(dāng)初始軌道近似為極軌道（i ≈90°）時(shí)，離軌時(shí)間最長(zhǎng)，初始軌道傾角在0°和170°情形下的離軌性能基本一致，這主要是由于地磁場(chǎng)圍繞地磁軸近似均勻分布，但地磁軸與地軸并不重合，存在約11.5°的夾角。

依據(jù)仿真結(jié)果，可以分析得到如下結(jié)論：

a）對(duì)于500 km以下，非極軌道的電動(dòng)力繩系離軌任務(wù)，裸繩系與絕緣繩系系統(tǒng)離軌性能相當(dāng)。相比而言，優(yōu)先選用裸繩系，其無需額外的電荷收集裝置，電荷交換效率較高，其離軌時(shí)間不超過1年。若離軌時(shí)間不苛刻，可以縮短繩系長(zhǎng)度，以降低硬件成本。

b）對(duì)于500 km以上，非極軌道的電動(dòng)力繩系離軌任務(wù)，由于軌道變高，空間中的電子濃度稀薄，裸繩系系統(tǒng)離軌性能較差，因此，優(yōu)先選用絕緣繩系，其對(duì)于1 000 km以下的離軌任務(wù)耗時(shí)不超過2年。

c）對(duì)于極軌道附近的電動(dòng)力繩系離軌任務(wù)，由于地磁場(chǎng)的磁軸與軌道面法線接近垂直，電磁效應(yīng)較弱，產(chǎn)生的洛倫茲力較小，離軌性能受到較大制約。相比而言，優(yōu)先選用絕緣繩系，其對(duì)于1 000 km的離軌任務(wù)，將會(huì)耗時(shí)5～6年。若離軌時(shí)間要求較高，可以合理增加繩系長(zhǎng)度，縮短離軌時(shí)間。

4 結(jié) 論

電動(dòng)力繩系離軌技術(shù)是一種新穎的離軌方式。本文綜合考慮電磁場(chǎng)、大氣阻力及地球引力場(chǎng)等多種因素，建立了多場(chǎng)耦合攝動(dòng)下的電動(dòng)力繩系離軌動(dòng)力學(xué)模型。隨后，基于數(shù)值仿真，驗(yàn)證了電動(dòng)力繩系離軌技術(shù)的有效性及優(yōu)越性，并分析了不同任務(wù)參數(shù)和系統(tǒng)參數(shù)對(duì)離軌性能的影響，同時(shí)提出了系統(tǒng)參數(shù)設(shè)置原則。最后，通過任務(wù)適應(yīng)性分析表明，一方面，電動(dòng)力繩系離軌技術(shù)更適用于非極軌道的離軌任務(wù)，另一方面，對(duì)于極軌道離軌任務(wù)，地磁場(chǎng)的特殊分布及空間電荷的分布情況制約了電動(dòng)力繩系離軌技術(shù)的性能，但這一情況可以通過改變系統(tǒng)參數(shù)（如繩系長(zhǎng)度）使之得以改善。

[1] Zhai G, Zhang J R. Space tether net system for debris capture and removal[C]. Washington: International Conference on Intelligent Human-Machine Systems and Cybernetics, 2012.

[2] Rouleau G, Verma S, Sharf I, Martin E. Vision-based tracking and trajectory generation for robotic capture of objects in space[C]. California: AIAA Guidance, Navigation and Control Conference, 2005.

[3] Aslanov V, Yudintsnev V. Dynamics of large space debris removal using tethered space tug[J]. Acta Astronautica, 2013, 91(10): 149-156.

[4] 蔡洪, 楊育偉, 郭才發(fā). 電動(dòng)力繩系研究進(jìn)展[J]. 宇航學(xué)報(bào), 2014(11): 1223-1232.

[5] Forward R, Hoyt R. Application of the terminator tether electrodynamic drag technology to the deorbit of constellation spacecraft[C]. Cleveland: AIAA/ASME/SAE/ASEE Joint Propulsion Conference, 1998.

[6] Stevens R, Baker W. Optimal control of a librating electrodynamic tether performing a multirevolution orbit change[J]. Guidance, Control and Dynamics, 2009, 32(5): 1497-1507.

[7] Wen H, Jin D, Hu H. Advances in dynamics and control of tethered satellite systems[J]. Acta Berneticsasdfasdf Mechanica Sinica, 2008, 24(3): 229-241.

[8] Sanmartin J, Lorenzini E, Martinez-Sanchez M. Electrodynamictether applications and constraints[J]. Journal of Spacecraft and Rockets, 2010, 47(3): 442-456.

[9] Sanmartin J, Martinez-Sanchez M, Ahedo E. Bare wire anodes of electrodynamic tethers[J]. Journal of Propulsion and Power, 1993, 9(3): 353-360.

[10] Zhong R, Zhu Z H. Dynamics of nanosatellite deorbit by bare electrodynamic tether in low earth orbit[J]. Journal of Spacecraft and Rocket, 2013, 50(3): 691-700.

The Deorbit Performance and Applicability Analysis of the Electrodynamics Te-thered Deorbit System

Zhang Feng, Shen Lin, Wu Sheng-bao, Wang Xiao-ding, Wang Shu-ting
(R&D Center, China Academy of Launch Vehicle Technology, Beijing, 100076)

The electrodynamics tethered deorbit technology is a novel way to remove abandoned spacecraft like upper stages or unusable satellites. The deorbit performance and applicability of the electrodynamics tethered deorbit system are analyzed. The electrodynamics tethered deorbit dynamics with multi-perturbation is formulated. Furthermore, the key impact factors, including mission parameters and system parameters, are analyzed. The following numerical simulations are undertaken to show the influence of the impact factors to the deorbit performance, based on which the applicability of the electrodynamics tethered deorbit technology are further analyzed. The analysis concludes that this technology shows a better deorbit performance for missions with non-polar orbits.

The electrodynamics tethered deorbit technology; Dynamics modeling; Deorbit performance analysis; Applicability analysis

V412.4

A

1004-7182(2017)03-0006-06

10.7654/j.issn.1004-7182.20170302

2015-10-08；

2016-02-02

張 烽（1985-），男，博士，高級(jí)工程師，主要研究方向?yàn)轱w行器制導(dǎo)與控制