司驥躍，龐兆君，付 杰，杜忠華

(1. 南京理工大學機械工程學院, 南京 210094; 2. 北京機械設備研究所, 北京, 100854; 3. 上海航天控制技術研究所, 上海 201109)

0 引 言

自1957年第一顆人造衛(wèi)星發(fā)射以來，日趨頻繁的空間活動使得空間碎片的數(shù)量顯著增加，嚴重威脅著國際空間站(ISS)等關鍵的空間設施的運行安全。研究表明，每年須主動清除5到10個空間碎片才能維持空間環(huán)境的長期穩(wěn)定[1]?？臻g繩網(wǎng)系統(tǒng)是一種具有良好應用前景的空間碎片主動清除技術[2]，相比于傳統(tǒng)的剛性捕獲技術[3]，其具有遠距離抓捕、輕質(zhì)靈活和費效比高等優(yōu)勢[4]。

搭載空間繩網(wǎng)系統(tǒng)的母平臺通過發(fā)射柔性繩網(wǎng)對空間碎片進行可靠包裹，并通過平臺與繩網(wǎng)之間的拖曳繩將空間碎片拖至墳墓軌道或大氣層內(nèi)燒毀。從早期“ROGER”項目對于空間繩網(wǎng)系統(tǒng)的方案設計及任務規(guī)劃[5]，到2012-2017年期間的拋物線飛行試驗[6]、真空展開試驗[7]、地面展開試驗[8]和繩網(wǎng)收口試驗[9]，再到2018年9月英國薩里大學首次在軌完成空間繩網(wǎng)抓捕試驗[10]，空間繩網(wǎng)捕獲技術的可行性得到了較好的驗證。

空間繩網(wǎng)的理論研究主要集中在繩網(wǎng)動力學建模方法、繩網(wǎng)展開動力學和繩網(wǎng)碰撞動力學等方面。彈簧質(zhì)點法[11]和絕對節(jié)點坐標法[12]是主要的繩網(wǎng)動力學建模方法，兩者各有優(yōu)劣[13]。繩網(wǎng)展開動力學方面主要關注繩網(wǎng)發(fā)射參數(shù)，如發(fā)射速度、發(fā)射角度等對于繩網(wǎng)展開性能的影響[14-16]。同時一些學者也從發(fā)射模式、繩網(wǎng)構(gòu)型等方面入手，開展了空間繩網(wǎng)高效展開方法方面的研究[17-18]。繩網(wǎng)的碰撞動力學建模方法主要有罰函數(shù)法[19]和沖量法[20-21]，其中罰函數(shù)法應用較為廣泛。捕獲階段的研究主要集中在對具有規(guī)則外形、不規(guī)則外形目標的抓捕過程仿真，以及目標靜止或旋轉(zhuǎn)時的抓捕過程[22-23]。空間繩網(wǎng)的自碰撞算法也被提出，并用于抓捕后繩網(wǎng)纏繞現(xiàn)象的研究[24]。

現(xiàn)有文獻中，關于繩網(wǎng)捕獲性能及其影響因素的研究較少。Endo等[25]研究了發(fā)射平臺與抓捕目標之間存在位置偏移時的抓捕過程，給出了特定工況下繩網(wǎng)抓捕可接受的位置偏移量，但是沒有考慮抓捕過程中繩網(wǎng)與目標之間的摩擦力。Salvi等[22]研究了發(fā)射速度、發(fā)射角度等發(fā)射參數(shù)不對稱對繩網(wǎng)展開性能的影響，但是研究僅關注了繩網(wǎng)的展開階段，沒有涉及繩網(wǎng)與目標接觸碰撞過程。另一方面，現(xiàn)階段研究一般認為繩網(wǎng)展開率超過80%時能抓捕成功[14-15,26]。雖然這一標準在衡量繩網(wǎng)展開、抓捕性能方面提供了統(tǒng)一的衡量標準，但是缺少理論或仿真支持。

本文以繩網(wǎng)邊長為度量確定了抓捕目標的極限包絡，并定義了理想和非理想工況下的繩網(wǎng)成功抓捕的判據(jù)。針對理想抓捕工況，確定了繩網(wǎng)可以成功抓捕的展開率范圍；針對非理想工況，分別研究了存在發(fā)射速度誤差、發(fā)射不同步和抽出阻力差異時的容許誤差范圍，并在抓捕過程中考慮了繩網(wǎng)與目標間的摩擦力。

1 繩網(wǎng)動力學模型

1.1 繩網(wǎng)建模方法

使用彈簧質(zhì)點模型對網(wǎng)格為正方形的四邊形繩網(wǎng)進行建模，繩網(wǎng)網(wǎng)角依次連接四個牽引體。將繩網(wǎng)物理節(jié)點視為集中質(zhì)量，節(jié)點間的繩段視為彈簧和阻尼器并聯(lián)的無質(zhì)量單元，繩網(wǎng)的構(gòu)型及離散方法如圖1所示，繩網(wǎng)中心位于坐標原點處。

圖1 繩網(wǎng)構(gòu)型Fig.1 Configuration of the tether-net

質(zhì)點i的動力學方程為：

(1)

式中:mi為質(zhì)點i的等效質(zhì)量，根據(jù)質(zhì)點在繩網(wǎng)內(nèi)的位置分布有如下定義：

(2)

Fext,i=Fn·en+Ft·et

(3)

式中:Fn是法向接觸力;Ft是切向接觸力;en和et分別是是法向接觸力和切向接觸力的方向向量。Fint,i是相鄰繩段傳遞給質(zhì)點的內(nèi)力，表示如下：

(4)

當繩段端點距離大于繩段原始長度時，繩段中產(chǎn)生拉力;當繩段端點距離小于繩段原始長度時繩段內(nèi)張力為零。因此在以i、j為端點的繩段中，i點所受張力如下：

(5)

式中:lij是繩段的真實長度;l0是繩段的原始長度;eij是從i指向j的單位向量;vij是質(zhì)點i與質(zhì)點j的相對速度;剛度系數(shù)kij和阻尼系數(shù)dij表示如下：

(6)

E和ζ分別是繩網(wǎng)選用材料的楊氏模量和阻尼比;A是網(wǎng)繩的橫截面積。

1.2 碰撞動力學模型

繩網(wǎng)碰撞動力學建模的核心問題是碰撞檢測及碰撞響應。碰撞檢測分析在仿真中的每一步都會進行，用來判斷檢測目標之間是否有碰撞發(fā)生，也為碰撞響應的計算提供了依據(jù)。本文算例中目標包絡為圓球形，繩網(wǎng)與目標之間的碰撞檢測可以通過判斷繩網(wǎng)節(jié)點中心與圓球中心的距離實現(xiàn)。如圖2所示，當目標與節(jié)點之間出現(xiàn)穿透時，則節(jié)點與目標發(fā)生碰撞。記δ為碰撞穿透量，則判定碰撞發(fā)生的依據(jù)如下：

δ=(Ri+Robj)-|xi-xobj|>0

(7)

其中:Ri和Robj分別是節(jié)點i和目標的半徑;xi和xobj分別是當前時刻繩網(wǎng)節(jié)點i與目標的位置矢量。

圖2 碰撞檢測及穿透量定義Fig.2 Collision detection and definition of penetration

碰撞響應主要有沖量法和罰函數(shù)法兩類。沖量法計算效率較高但是無法計算碰撞過程中的接觸力，碰撞速度也不連續(xù)。而罰函數(shù)法可以得到碰撞過程中碰撞力的時間變化歷程，同時這種方法適用于涉及多體接觸的復雜接觸場景，但是為了防止出現(xiàn)過大的接觸力，數(shù)值仿真中步長必須設置的足夠小。本文采用罰函數(shù)法中的非線性阻尼模型[27]建立繩網(wǎng)與目標的碰撞模型。發(fā)生碰撞時，法向接觸力取為：

(8)

(9)

式中:r為與幾何構(gòu)型有關的等效半徑;ν和E分別是碰撞體的泊松比和楊氏模量;α是與恢復系數(shù)有關的經(jīng)驗參數(shù)。

繩網(wǎng)與目標之間還存在切向接觸力，即摩擦力。根據(jù)庫侖定律，動摩擦與法向接觸力成正比，引入庫倫摩擦系數(shù)μd。并引入切向速度閾值vs作為動靜摩擦的判定依據(jù)，當繩網(wǎng)與目標的切向速度模量vt大于等于vs時認為兩者之間為動摩擦，反之為靜摩擦。因此切向接觸力可以表示為[19,29]：

(10)

實驗證明當vs取值在10-4到10-6m/s之間時，能夠在計算效率和求解精度之間達到較好的平衡[30]。

2 理想工況下繩網(wǎng)捕獲性能分析

2.1 繩網(wǎng)可抓捕目標的極限尺寸

本文研究邊長為L的具有正方形網(wǎng)格的四邊形繩網(wǎng)，首先根據(jù)繩網(wǎng)尺寸確定理想工況下可以抓捕的目標尺寸范圍。繩網(wǎng)的鏤空結(jié)構(gòu)使得其具有輕便的質(zhì)量，而其較大的外形尺寸使繩網(wǎng)相比剛性抓捕具備更大的抓捕容錯率。抓捕目標的最小尺寸應不小于繩網(wǎng)的網(wǎng)目大小，以防止捕獲后目標從繩網(wǎng)網(wǎng)孔內(nèi)脫離。因此定義繩網(wǎng)網(wǎng)目邊長不大于目標包絡圓的內(nèi)接正方體邊長。繩網(wǎng)還應確保將目標完全包裹，故繩網(wǎng)的邊長應不小于目標包絡圓的周長。綜上，繩網(wǎng)可抓捕的最小目標半徑為：

(11)

可抓捕的最大目標半徑為：

(12)

其中:L為繩網(wǎng)邊長;l為網(wǎng)目邊長。繩網(wǎng)及其可抓捕目標的對比圖如圖3所示。

圖3 空間繩網(wǎng)可抓捕目標的極限尺寸Fig.3 Limit size of the target

2.2 成功捕獲判據(jù)

相比繩網(wǎng)展開階段，繩網(wǎng)的抓捕過程更為復雜，因此首先對抓捕條件進行合理假設和簡化：抓捕過程中不考慮牽引體的外形，將其作為質(zhì)點處理；繩網(wǎng)關于網(wǎng)體中心完全對稱；目標的長徑比接近1。對長徑比進行限定主要有兩點原因：一方面，長徑比較大時目標呈現(xiàn)細長形狀，其包絡圓的尺寸可能較大,但是仍有可能在捕獲過程中從網(wǎng)目中逃脫；另一方面若目標外形為細桿狀或太陽帆等平板狀結(jié)構(gòu)，其與繩網(wǎng)的接觸作用過程可能異于本文的研究目標。

理想工況是指繩網(wǎng)及其發(fā)射系統(tǒng)具有較好的對稱性及一致性，且繩網(wǎng)中心和目標中心的連線與繩網(wǎng)飛行方向重合。理想工況下的抓捕過程為：繩網(wǎng)網(wǎng)體首先與目標接觸，牽引體繼續(xù)飛行直到達到極限位置后回彈并帶動網(wǎng)角對目標形成包裹的趨勢。隨后依靠繩網(wǎng)的自纏繞或通過收口機構(gòu)對目標形成可靠包裹。根據(jù)這一過程將理想工況下繩網(wǎng)成功抓捕的判據(jù)定義為：每個牽引體到達的極限位置均超過目標的包絡上限，且這一過程中牽引體與目標不應發(fā)生碰撞，以防止產(chǎn)生新的碎片。包絡上限定義為目標包絡遠離發(fā)射點的一側(cè)，包絡下限相反，如圖4所示。

圖4 繩網(wǎng)抓捕閾值定義Fig.4 Definition of capture threshold value of tether-net

繩網(wǎng)的展開面積Sdep一般定義為繩網(wǎng)4個網(wǎng)角所圍成四邊形的面積，繩網(wǎng)的展開面積與設計面積Sori的比值為展開率。在以往的繩網(wǎng)展開性能研究中，一般認為繩網(wǎng)展開率超過80%時，繩網(wǎng)可以抓捕成功。這一抓捕閾值為不同的展開性能研究提供了相同的衡量標準，但是實際抓捕過程中，80%的設定并沒有明確的依據(jù)。因此使用上文定義的抓捕成功判據(jù)，對不同展開率下的繩網(wǎng)的抓捕性能進行研究。

繩網(wǎng)達到最大展開率時會發(fā)生回彈，其展開率先增大后減小，變化過程是一個單峰曲線。因此以最大展開率為界，將繩網(wǎng)展開率增大階段可以成功抓捕的最小抓捕展開率稱為第一閾值，并記為D1；將展開率減小階段可以成功抓捕的最小抓捕展開率稱為第二閾值，記為D2，如圖4所示。理想工況下，只要繩網(wǎng)抓捕目標時的展開率處于D1及D2之間，繩網(wǎng)便能成功抓捕目標。且D1和D2越小可抓捕區(qū)間越大，抓捕性能越高。需要說明的是，在繩網(wǎng)展開率增大階段，抓捕展開率為網(wǎng)角到達目標包絡下限時的展開率；在繩網(wǎng)展開率減小階段，抓捕展開率為繩網(wǎng)網(wǎng)體與目標接觸時刻的展開率。下文將對抓捕極限閾值的確定進行討論。

2.3 抓捕極限閾值的確定

首先討論D1的確定，此過程繩網(wǎng)展開率由小變大。若繩網(wǎng)在此過程對目標進行抓捕，必然是牽引體首先越過目標包絡下限，之后繩網(wǎng)與目標接觸并對目標進行包裹?？紤]到牽引體不能與目標發(fā)生碰撞，因此理想工況下牽引繩延伸線與目標包絡相切，且繩網(wǎng)網(wǎng)角恰好到達目標下包絡線時繩網(wǎng)的展開率即為D1。這一階段繩網(wǎng)網(wǎng)角被牽引體逐漸抽出且呈現(xiàn)為束狀，由于發(fā)射初期牽引體拖動的繩網(wǎng)質(zhì)量較低，因此束狀網(wǎng)角與繩網(wǎng)發(fā)射方向之間的角度可以認為與發(fā)射角相同。達到D1時繩網(wǎng)與目標的位置關系如圖5所示。

圖5 第一閾值確定示意Fig.5 Determination of D1

將目標半徑記為Robj，發(fā)射角記為θ，處于對角位置的網(wǎng)角連線長度記為s。則此時繩網(wǎng)的展開率為：

(13)

根據(jù)圖5可知：

(14)

因此，發(fā)射點O到目標中心Oobj的距離可以表示為：

(15)

可知

(16)

可得

(17)

取最惡劣的情況，即目標半徑為最大，可將式(12)、式(13)和式(17)聯(lián)立得：

(18)

通過式(18)可以發(fā)現(xiàn)繩網(wǎng)的第一閾值僅與發(fā)射角度有關，可以繪制出第一閾值與發(fā)射角度的關系，如圖6所示。在給定的發(fā)射角變化范圍內(nèi)，隨著發(fā)射角度的增加D1逐漸減小，D1的最大值不超過4.5%。說明理想工況下，繩網(wǎng)在回彈前的抓捕過程中具有較強的容錯能力。

圖6 D1與發(fā)射角度的關系Fig.6 Relationship betweenD1 and shooting angle

繩網(wǎng)達到最大展開率后中心回彈，展開率減小。此過程繩網(wǎng)形態(tài)與目標包絡沒有明確的幾何關系，因此采用控制變量法研究不同工況下第二閾值D2與繩網(wǎng)發(fā)射參數(shù)之間的關系。影響繩網(wǎng)展開性能的發(fā)射參數(shù)主要有發(fā)射速度vb、發(fā)射角度θ和牽引體總質(zhì)量與繩網(wǎng)質(zhì)量的比值η，將三個參數(shù)設定一系列水平進行研究。其中發(fā)射速度和牽引體質(zhì)量對繩網(wǎng)與目標接觸時的動能和形態(tài)影響較大，記為一組進行控制變量研究。發(fā)射角度對繩網(wǎng)的展開時間和展開位移有影響，單獨研究。仿真參數(shù)見表1及表2。

表1 仿真參數(shù)Table 1 Simulation parameters

表2 繩網(wǎng)的輸入?yún)?shù)Table 2 Input parameters of tether-net

在給定發(fā)射參數(shù)后，根據(jù)2.2節(jié)中的定義“在繩網(wǎng)展開率減小階段，抓捕展開率為繩網(wǎng)網(wǎng)體與目標接觸時刻的展開率”，通過調(diào)整發(fā)射點與目標之間的距離調(diào)整抓捕展開率，通過一系列仿真得到的可以成功抓捕的最小抓捕展開率即為D2。目標包絡半徑按照式(12)取最大為1.59 m。首先保持發(fā)射角為40°繪制不同發(fā)射速度和質(zhì)量系數(shù)下的D2，如圖7所示?？梢园l(fā)現(xiàn)，發(fā)射速度對于D2的大小影響不大，而隨著牽引體質(zhì)量的增加D2逐漸減小。盡管提高牽引體質(zhì)量可以減低D2，但是這種方法存在著較大的局限性。一方面提高牽引體質(zhì)量會提升系統(tǒng)質(zhì)量及發(fā)射難度，另一方面當質(zhì)量系數(shù)超過7時，D2穩(wěn)定在78%附近，抓捕性能提升效果不佳。

圖7 D2的影響因素Fig.7 The influencing factors ofD2

保持發(fā)射速度為10 m/s，質(zhì)量系數(shù)為3繪制發(fā)射角與D2的關系，如圖8所示。可以看出隨著發(fā)射角的增加D2逐漸增加，抓捕性能逐漸下降。在發(fā)射角超過50°后，D2穩(wěn)定在95%附近。相比增加牽引體質(zhì)量，降低發(fā)射角對于提升繩網(wǎng)回彈之后的抓捕性能效果更好。當發(fā)射角為20°時，D2的值僅為38.2%。

圖8 D2與發(fā)射角的關系Fig.8 Relationship betweenD2and shooting angle

通過推導和仿真得到了以展開率D1和D2為界的繩網(wǎng)可抓捕區(qū)間。D1和D2的值越小，則繩網(wǎng)的可抓捕區(qū)間越大，繩網(wǎng)的抓捕性能越好。D1僅與發(fā)射角θ負相關，較大的發(fā)射角可以獲得較小的D1值。但是D1在各發(fā)射角下變化幅度很小。D2對于牽引體質(zhì)量和發(fā)射角度較為敏感，較大的牽引體質(zhì)量和較小的發(fā)射角可以獲得較好的抓捕性能，同時提高牽引體質(zhì)量弊端較多。綜上，建議使用較小的發(fā)射角以獲得最佳的抓捕區(qū)間。

3 非理想工況下繩網(wǎng)捕獲性能分析

3.1 非理想工況設定

將抓捕時機設定為繩網(wǎng)展開率最大的時刻，抓捕目標選為最大可抓捕尺寸，對非理想工況下繩網(wǎng)捕獲性能進行研究。繩網(wǎng)發(fā)射過程中可能出現(xiàn)不可控誤差的情景主要有：1)由于動力源的推力誤差所引起的發(fā)射速度不一致；2)由控制系統(tǒng)或擊發(fā)系統(tǒng)誤差引起的發(fā)射時間不同步；3)由于繩網(wǎng)折疊工藝引起的繩網(wǎng)抽出阻力差異。將這三種難以徹底排除的非理想工況作為研究對象，對繩網(wǎng)抓捕過程進行研究，得出保證捕獲成功的誤差邊界。另外，還以存在發(fā)射速度誤差為例，討論了抓捕距離誤差對抓捕性能的影響。

理想工況下成功抓捕的判據(jù)是：每個牽引體到達的極限位置均超過目標的包絡上限，且這一過程中牽引體與目標不發(fā)生碰撞。而非理想工況下繩網(wǎng)中心與目標中心會出現(xiàn)偏移，繩網(wǎng)可能從目標表面滑過造成抓捕失效。因此在理想判據(jù)的前提下增加一項約束條件作為成功抓捕判據(jù)：出現(xiàn)誤差的牽引體回落到目標包絡上限內(nèi)之前，繩網(wǎng)中心點在飛行方向(Z軸正向)不超過目標中心。將目標中心在Z軸的坐標記為zobj，將抓捕過程中繩網(wǎng)中心點在Z軸所到達到的距離最大值記為zcen。當zobj-zcen>0時繩網(wǎng)中心點在Z軸正向不超過目標中心，抓捕有效；當zobj-zcen<0時抓捕失敗。

3.2 數(shù)值仿真

以發(fā)射速度為10 m/s，發(fā)射角度為40°，質(zhì)量比為4作為基本輸入，目標包絡尺寸按照式(12)取最大。目標中心位于Z軸正方向12.28 m處，此距離下目標與繩網(wǎng)接觸時繩網(wǎng)處于最大展開率。其他仿真條件見表1。根據(jù)不同的非理想工況，在基礎工況上疊加誤差并進行動力學仿真，最終得到繩網(wǎng)系統(tǒng)可接受的最大誤差值。在本研究中針對不同的非理想工況設定了兩個子工況，分別是單個牽引體出現(xiàn)誤差、兩個相鄰牽引體出現(xiàn)誤差。兩個相鄰牽引體出現(xiàn)誤差時認為其誤差值相同。下文中單個牽引體的誤差施加在牽引體1上，相鄰牽引體的誤差施加在牽引體1和牽引體2上，牽引體編號方法如圖1所示。另外，當處于對角位置的牽引體出現(xiàn)相同誤差時繩網(wǎng)中心與目標中心不會出現(xiàn)位置偏移，因此這種工況在本文不作考慮。本節(jié)考慮摩擦，摩擦系數(shù)為0.2，切向速度閾值vs為1×10-4m/s。

首先對發(fā)射速度誤差進行研究。根據(jù)仿真結(jié)果，以誤差值為橫坐標，以目標中心與繩網(wǎng)中心在Z軸的距離差值為縱軸繪圖，如圖9所示。可以發(fā)現(xiàn)當單個牽引體的發(fā)射速度低于其他牽引體時，繩網(wǎng)系統(tǒng)可接受的速度最大偏差為-3.15 m/s，此時牽引體1的發(fā)射速度為6.85 m/s。此誤差下，各牽引體在發(fā)射方向的飛行距離如圖10(a)所示，由于對稱性牽引體2與牽引體4的飛行軌跡基本重合。4個牽引體在Z軸的飛行距離隨時間而增加，4個牽引體具有相似的飛行軌跡，但是牽引體1由于發(fā)射速度較小，其飛行距離小于其他牽引體。繩網(wǎng)在約1.86 s時與目標發(fā)生碰撞并逐漸包裹目標，之后牽引體越過目標包絡上限繼續(xù)飛行，直到約3.82 s時達到極限位置開始回彈。此過程中每個牽引體到達的極限位置均超過目標的包絡上限，且繩網(wǎng)中心點在Z軸的距離總是小于目標中心，因此抓捕有效。單個牽引體的發(fā)射速度高于其他牽引體時，繩網(wǎng)系統(tǒng)可接受的速度最大偏差為3.55 m/s，此時牽引體1的發(fā)射速度為13.55 m/s，各牽引體飛行軌跡如圖10(b)所示。

圖9 繩網(wǎng)可接受的發(fā)射速度誤差范圍Fig.9 The acceptable error of shooting velocity of tether-net

圖10 單牽引體速度偏差下牽引體在Z軸的飛行距離Fig.10 The distance of single bullet in Z-axis under shooting velocity error

從圖9中還可知，誤差存在于相鄰牽引體時繩網(wǎng)系統(tǒng)可接受的發(fā)射速度偏差范圍為-2.31～2.94 m/s，偏差量明顯低于單牽引體工況。牽引體1、2的發(fā)射速度誤差為-2.31 m/s時牽引體的飛行軌跡如圖11(a)，牽引體1、2的速度誤差為2.94 m/s時軌跡如圖11(b)所示。由于對稱性，牽引體1與牽引體2的軌跡基本一致，牽引體3與牽引體4的軌跡基本一致。牽引體1、2回落到目標包絡上限內(nèi)之前，繩網(wǎng)中心點在Z軸正向不超過目標中心，抓捕有效。

圖11 雙牽引體速度偏差下牽引體在Z軸的飛行距離Fig.11 The distance of two bullets in Z-axis under shooting velocity error

接著對發(fā)射同步性誤差進行研究，繪制了繩網(wǎng)可接受的發(fā)射同步性誤差范圍，如圖12所示。單個牽引體存在發(fā)射同步性誤差時，可接受的發(fā)射時間誤差范圍為-0.62～0.78 s。相鄰牽引體存在發(fā)射同步性誤差時，可接受的時間誤差范圍為-0.46～0.46 s。所研究繩網(wǎng)系統(tǒng)的發(fā)射同步性容錯范圍為百毫秒級別，具有較好的抓捕容錯能力。

圖12 繩網(wǎng)可接受的發(fā)射同步性誤差范圍Fig.12 The acceptable error of shooting synchronization of tether-net

發(fā)射同步性存在誤差時，各牽引體飛行軌跡與圖10、圖11類似，因此不再贅述。圖13繪制出了牽引體1延時發(fā)射0.62 s時繩網(wǎng)的抓捕歷程。由于牽引體1延遲發(fā)射，2 s前其在發(fā)射方向的位移明顯小于其他三個牽引體，但是得益于繩網(wǎng)良好的抓捕容錯能力，2 s之后繩網(wǎng)逐漸包裹目標，并最終完成抓捕。

圖13 發(fā)射同步性誤差下繩網(wǎng)成功抓捕歷程Fig.13 Capture sequence with shooting synchronization error

最后研究繩網(wǎng)抽出阻力不一致對于抓捕過程的影響。為了便于對比，認為理想情況下繩網(wǎng)對各牽引體造成的抽出阻力均為0 N。根據(jù)圖14可知，若牽引體1所受抽出阻力大于其他牽引體時，此抽出阻力差值最大為1.41 N。相鄰牽引體的抽出阻力存在誤差時，最大的阻力差值為1.3 N。高慶玉等[31]對繩網(wǎng)抽出阻力的研究顯示，繩網(wǎng)從網(wǎng)包抽出時可以將總的抽出阻力控制在5 N之內(nèi)。若按照4個牽引體計算，則每個牽引體的平均抽出阻力為1.25 N，小于本研究得出的最大誤差值，說明了空間繩網(wǎng)在抽出阻力差異方面具有較好的容錯能力。表3總結(jié)了上述各工況下可接受的誤差范圍。

以上研究中目標被固定在繩網(wǎng)展開率最大處，目標的中心位置距離發(fā)射平臺12.28 m。如果發(fā)射平臺與目標之間的距離存在誤差，則會放大發(fā)射誤差對于捕獲性能的影響。針對這一情況，以雙牽引體存在發(fā)射速度誤差為例，研究繩網(wǎng)可接受的抓捕距離誤差范圍。牽引體的速度誤差參照表3取為2.5 m/s，則牽引體1、2的發(fā)射速度為12.5 m/s，牽引體3、4的發(fā)射速度為10 m/s。發(fā)射距離以12.28 m為基數(shù)，當發(fā)射距離小于12.28時記誤差為負，大于時記誤差為正。繪制了繩網(wǎng)可接受的發(fā)射距離誤差范圍，如圖15所示。當發(fā)射距離誤差小于-34.1%時，雖然由速度誤差所造成的繩網(wǎng)中心偏移量還不大，但是此時繩網(wǎng)的展開率較小，對于發(fā)射誤差較為敏感，并未滿足捕獲成功判據(jù)。當距離誤差介于28.3%和85%之間時，出現(xiàn)誤差的牽引體回落到目標包絡上限內(nèi)之前，繩網(wǎng)中心點在飛行方向超過了目標中心，判定為抓捕失敗。當距離誤差大于85%時，由于速度誤差的累積，繩網(wǎng)中心與目標中心在飛行方向上已經(jīng)存在了較大的偏移量，繩網(wǎng)與目標接觸后直接從目標一側(cè)滑過，造成脫靶。此工況下，繩網(wǎng)可接受的抓捕距離誤差范圍為-34.1%～60.7%，對于抓捕距離具有較好的容錯范圍。

圖14 繩網(wǎng)可接受的抽出阻力誤差范圍Fig.14 The acceptable error of extraction force of tether-net

表3 繩網(wǎng)可接受的誤差范圍Table 3 Acceptable error range of tether-net

圖15 繩網(wǎng)可接受的抓捕距離誤差范圍Fig.15 The acceptable error of capture distance of tether-net

4 結(jié) 論

本文重點關注了空間繩網(wǎng)的抓捕動力學模型，并對繩網(wǎng)的捕獲性能進行了研究。使用彈簧質(zhì)點法建立了繩網(wǎng)的動力學模型，基于非線性阻尼模型建立了繩網(wǎng)與目標的碰撞動力學模型，同時根據(jù)庫倫定律在碰撞模型中記入了摩擦力。根據(jù)繩網(wǎng)邊長確定了可抓捕目標的最大、最小包絡。提出了理想工況和非理想工況下界定繩網(wǎng)抓捕成功的判據(jù)，并研究了發(fā)射速度誤差、發(fā)射同步性、抽出阻力差異和抓捕距離誤差對繩網(wǎng)抓捕性能的影響，給出了繩網(wǎng)系統(tǒng)可接受的誤差邊界。通過研究得出以下結(jié)論：理想工況下以展開率80%作為抓捕成功的界定標準過于嚴格，繩網(wǎng)的抓捕性能在回彈前與發(fā)射角有關，在回彈后與發(fā)射角和牽引體質(zhì)量關系較大；相比增加牽引體質(zhì)量，減小發(fā)射角對于提升繩網(wǎng)回彈之后的抓捕性能效果更好；非理想工況下發(fā)射速度誤差、發(fā)射同步性和抽出阻力差異均會影響繩網(wǎng)的抓捕性能，兩個牽引體出現(xiàn)誤差時可接受的誤差范圍劣于單個牽引體；過大的抓捕距離誤差也會造成捕獲失敗，特別是抓捕距離較大時，繩網(wǎng)有可能直接脫靶。