張廣華，馮順山，葛 超，王 超，方國(guó)燚

（北京理工大學(xué) 爆炸科學(xué)與技術(shù)國(guó)家重點(diǎn)實(shí)驗(yàn)室，北京100081）

導(dǎo)電纖維絲團(tuán)是導(dǎo)電纖維彈中的有效裝填物，用于毀傷目標(biāo)。絲團(tuán)一般以子母彈的形式裝填在母彈當(dāng)中，根據(jù)母彈彈體結(jié)構(gòu)的不同，裝填絲團(tuán)的個(gè)數(shù)一般為幾十個(gè)到上百個(gè)不等。每個(gè)絲團(tuán)上都纏繞有幾十m長(zhǎng)的導(dǎo)電絲束，用于搭掛于目標(biāo)上并引發(fā)短路效應(yīng)。

絲團(tuán)被拋出后，在運(yùn)動(dòng)慣性力、重力及氣動(dòng)阻力的聯(lián)合作用下，一邊向下運(yùn)動(dòng)一邊展開(kāi)；展開(kāi)后的絲束在空中飄移下落一段距離后以一定的概率搭掛于變電站或高壓線等供電設(shè)施上，通過(guò)其短路引弧效應(yīng)導(dǎo)致目標(biāo)乃至整個(gè)區(qū)域內(nèi)的電力系統(tǒng)癱瘓。

絲團(tuán)空中展開(kāi)及飄落的運(yùn)動(dòng)規(guī)律直接影響到展開(kāi)后的絲束能否準(zhǔn)確地搭落到被攻擊的目標(biāo)之上，并對(duì)其進(jìn)行有效的毀傷，并且可以對(duì)導(dǎo)電纖維彈的終端彈道設(shè)計(jì)（拋絲團(tuán)高度、速度等）提供依據(jù)。目前，國(guó)內(nèi)外關(guān)于絲團(tuán)及導(dǎo)電絲束空中運(yùn)動(dòng)規(guī)律的研究不多［1－3］，因此，本文建立了一種用于描述絲團(tuán)空中運(yùn)動(dòng)的彈道模型，闡述其空中運(yùn)動(dòng)特性。

1 絲團(tuán)運(yùn)動(dòng)模型

1.1 假設(shè)條件

導(dǎo)電纖維絲團(tuán)的結(jié)構(gòu)如圖1所示。

絲團(tuán)從母彈中被拋出時(shí)的運(yùn)動(dòng)姿態(tài)及空中運(yùn)動(dòng)時(shí)受到的風(fēng)力影響是不確定的，如果考慮所有因素，將會(huì)給研究帶來(lái)較大的困難。因此，為了便于建立空中運(yùn)動(dòng)模型，特做如下假設(shè)：①絲團(tuán)被拋出母彈時(shí)的運(yùn)動(dòng)姿態(tài)為水平運(yùn)動(dòng)；②絲束展開(kāi)過(guò)程中只考慮絲軸繞絲束的旋轉(zhuǎn)運(yùn)動(dòng)而忽略其它運(yùn)動(dòng)；③絲軸及絲束受到的風(fēng)力及風(fēng)向固定不變，風(fēng)向?yàn)樗椒较?；④不考慮絲束在展開(kāi)過(guò)程中的彎曲；⑤忽略阻力葉片繞絲束的轉(zhuǎn)動(dòng)。

圖1 導(dǎo)電纖維絲團(tuán)結(jié)構(gòu)圖

1.2 絲團(tuán)展開(kāi)模型

以大地坐標(biāo)系為基準(zhǔn)坐標(biāo)系，絲團(tuán)展開(kāi)過(guò)程示意圖如圖2所示。

圖2 絲團(tuán)展開(kāi)過(guò)程示意圖

絲束展開(kāi)長(zhǎng)度隨時(shí)間的變化規(guī)律：

式中：l為絲束展開(kāi)長(zhǎng)度為絲束展開(kāi)速度，R為絲團(tuán)的半徑為絲軸旋轉(zhuǎn)角速度。

絲軸旋轉(zhuǎn)時(shí)的角速度及角加速度的變化規(guī)律如下：

式中：y為絲軸（絲軸連同未展開(kāi)的絲束）在豎直方向的運(yùn)動(dòng)距離，yc為已展開(kāi)絲束在豎直方向的運(yùn)動(dòng)距離，θ為絲束與y軸間的夾角。

絲團(tuán)運(yùn)動(dòng)過(guò)程中，隨著絲束的不斷展開(kāi)，纏繞在絲軸上的絲團(tuán)半徑R也在不斷變化，R的變化規(guī)律為

式中：a為絲軸的長(zhǎng)度，一般取為1.5～2cm，d為單根絲束的直徑，一般取為0.2～0.5mm。

1.3 未展開(kāi)絲團(tuán)的運(yùn)動(dòng)模型

未展開(kāi)絲團(tuán)包括絲軸及纏繞在絲軸上的未展開(kāi)的絲束，其受力情況如圖3所示。

圖3 未展開(kāi)絲團(tuán)受力圖

圖3 中，F(xiàn)1x，F(xiàn)1y分別為絲團(tuán)受到的水平方向和豎直方向的空氣阻力；G1，G2分別為絲軸及未展開(kāi)絲束所受到的重力；FT為已展開(kāi)絲束對(duì)絲團(tuán)的拉力；Ff1為絲團(tuán)受到的風(fēng)動(dòng)力。根據(jù)圖3所示未展開(kāi)絲團(tuán)的受力情況，建立其運(yùn)動(dòng)方程：

式中：x，y分別為絲團(tuán)在基準(zhǔn)坐標(biāo)系下的水平運(yùn)動(dòng)距離和豎直運(yùn)動(dòng)距離；J為絲軸旋轉(zhuǎn)時(shí)的轉(zhuǎn)動(dòng)慣量；m為絲軸及未展開(kāi)絲束的質(zhì)量和，m＝mz－lρl，其中，mz為絲團(tuán)初始質(zhì)量，ρl為絲束線密度。

1.4 已展開(kāi)絲束的運(yùn)動(dòng)模型

已展開(kāi)絲束的受力情況如圖4所示。

圖4 已展開(kāi)絲束受力圖

圖4 中，F(xiàn)2x和F2y分別為已展開(kāi)絲束受到的水平方向和豎直方向的空氣阻力；F3x，F(xiàn)3y分別為阻力葉片受到的水平方向和豎直方向的空氣阻力；Ff2和Ff3分別為已展開(kāi)絲束及阻力葉片受到的風(fēng)的作用力；G3為已展開(kāi)絲束受到的重力。A點(diǎn)為已展開(kāi)絲束與絲軸的交點(diǎn)，C點(diǎn)為已展開(kāi)絲束的中心點(diǎn)，B點(diǎn)為阻力葉片與絲束的交點(diǎn)。

根據(jù)圖4所示已展開(kāi)絲束的受力分析，建立其運(yùn)動(dòng)方程：

式中：xc為已展開(kāi)絲束的水平運(yùn)動(dòng)距離；mc為已展開(kāi)絲束的質(zhì)量，mc＝lρl；J1為已展開(kāi)絲束繞絲束與絲軸的接觸點(diǎn)A轉(zhuǎn)動(dòng)時(shí)的轉(zhuǎn)動(dòng)慣量。

1.5 絲團(tuán)展開(kāi)及飄落運(yùn)動(dòng)方程組

聯(lián)立式（1）～式（10），輔以絲團(tuán)運(yùn)動(dòng)過(guò)程中運(yùn)動(dòng)距離、速度之間的微分關(guān)系式，得到了基于假設(shè)條件下的絲團(tuán)空中展開(kāi)與飄落外彈道模型：

本文以MATLAB軟件為基礎(chǔ)，以四階龍格－庫(kù)塔算法為求解方法，編制了相應(yīng)的程序?qū)Ψ匠探M（11）進(jìn)行數(shù)值求解［4］。結(jié)合相關(guān)試驗(yàn)（地面拋撒試驗(yàn)、靶場(chǎng)飛行試驗(yàn)）結(jié)果及方程組（11）反推出了方程組（11）中的相關(guān)氣動(dòng)參數(shù)，并根據(jù)文獻(xiàn)［5－6］提供的方法對(duì)絲團(tuán)及絲束所受氣動(dòng)力進(jìn)行計(jì)算；轉(zhuǎn)動(dòng)慣量根據(jù)文獻(xiàn)［7］提供的方法計(jì)算。

2 試驗(yàn)研究

為了對(duì)本文所建立模型的合理性進(jìn)行驗(yàn)證，進(jìn)行了單個(gè)絲團(tuán)的定向拋撒試驗(yàn)及靶場(chǎng)飛行試驗(yàn)。其中，拋撒試驗(yàn)主要用于研究絲團(tuán)的空中展開(kāi)規(guī)律，靶場(chǎng)飛行試驗(yàn)主要用于研究絲團(tuán)的空中飄落規(guī)律。

2.1 試驗(yàn)絲團(tuán)的技術(shù)狀態(tài)

試驗(yàn)絲團(tuán)的技術(shù)狀態(tài)如表1所示，表中，la和lb分別為阻力葉片的長(zhǎng)度和寬度，lss為纏繞在絲軸上的絲束長(zhǎng)度。

表1 試驗(yàn)絲團(tuán)的技術(shù)狀態(tài)

2.2 地面拋撒試驗(yàn)

地面拋撒試驗(yàn)的現(xiàn)場(chǎng)布局如圖5所示。

圖5 地面拋撒試驗(yàn)布局示意圖

地面發(fā)射裝置將絲團(tuán)發(fā)射出去后，通過(guò)高速攝像設(shè)備記錄其空中運(yùn)動(dòng)狀態(tài)；將高速攝像設(shè)備拍攝的結(jié)果導(dǎo)入到與之配套的后臺(tái)軟件上，并對(duì)其進(jìn)行相應(yīng)的后處理，即可得到絲團(tuán)空中展開(kāi)時(shí)的運(yùn)動(dòng)參數(shù)隨時(shí)間的變化情況。

2.3 靶場(chǎng)試驗(yàn)

2012年6月28日于某靶場(chǎng)進(jìn)行了導(dǎo)電纖維絲團(tuán)空中拋撒的飛行試驗(yàn)。

試驗(yàn)絲團(tuán)裝填于特定的拋撒裝置中進(jìn)行發(fā)射；拋撒裝置運(yùn)動(dòng)到一定高度后將絲團(tuán)從空中進(jìn)行拋撒，并通過(guò)三坐標(biāo)跟蹤雷達(dá)測(cè)量得到了絲團(tuán)空中運(yùn)動(dòng)時(shí)的各項(xiàng)參數(shù)。

3 結(jié)果分析

根據(jù)拋撒試驗(yàn)及靶場(chǎng)飛行試驗(yàn)的試驗(yàn)條件對(duì)方程組（11）中的各個(gè)變量賦予如表2所示的初始值。

表2 方程組（11）變量初始值

3.1 絲束展開(kāi)規(guī)律

選取典型時(shí)刻的絲束展開(kāi)長(zhǎng)度值為分析對(duì)象，將數(shù)值計(jì)算結(jié)果與拋撒試驗(yàn)結(jié)果進(jìn)行對(duì)比，對(duì)模型進(jìn)行驗(yàn)證。因?yàn)榈孛鎾伻鲈囼?yàn)的發(fā)射裝置距離地面較低，試驗(yàn)時(shí)絲團(tuán)在空中運(yùn)動(dòng)0.48s后便已著地，所以只能對(duì)0.48s之前的相關(guān)數(shù)據(jù)進(jìn)行對(duì)比分析。試驗(yàn)及數(shù)值計(jì)算的結(jié)果如表3所示，表中，le為試驗(yàn)值，ls為計(jì)算值，Δl為絕對(duì)誤差，Δl／le為相對(duì)誤差。

表3 絲束展開(kāi)長(zhǎng)度隨時(shí)間變化值

通過(guò)表3可以看出，采用方程組（11）求解出的結(jié)果與實(shí)驗(yàn)值之間的絕對(duì)誤差在0.03～0.09m之間，相對(duì)誤差在1.96%～4.66%之間，二者吻合情況較好。

當(dāng)拋撒速度為38m／s，風(fēng)力條件為無(wú)風(fēng)時(shí)，對(duì)方程組（11）進(jìn)行求解，求解時(shí)間設(shè)定為20s，得出絲束展開(kāi)長(zhǎng)度l隨時(shí)間t的變化情況如圖6所示。

圖6 絲束展開(kāi)長(zhǎng)度－時(shí)間關(guān)系

通過(guò)圖6可以看出，運(yùn)動(dòng)初期，絲束以接近均勻速度展開(kāi)，展開(kāi)速度為4.5～5m／s，運(yùn)動(dòng)到4.5s左右時(shí)完全展開(kāi)。

3.2 絲團(tuán)飄落規(guī)律

由于靶場(chǎng)試驗(yàn)時(shí)測(cè)試?yán)走_(dá)距離拋撒裝置較遠(yuǎn)，加之絲團(tuán)體積很小，想要對(duì)絲團(tuán)空中運(yùn)動(dòng)的全過(guò)程進(jìn)行跟蹤測(cè)試具有比較大的難度，所以只對(duì)其空中拋出位置及落地后的位置進(jìn)行了測(cè)量，由此可以確定絲團(tuán)在給定拋撒高度下的水平飄移距離。試驗(yàn)及數(shù)值計(jì)算的相關(guān)數(shù)據(jù)如表4所示，表中，yp為拋絲團(tuán)高度；xe，xs分別為絲團(tuán)飄移距離試驗(yàn)值和計(jì)算值，Δx為絕對(duì)誤差，Δx／xe為相對(duì)誤差。

表4 絲團(tuán)飄移距離

通過(guò)表4可以看出，絲團(tuán)飄移距離的計(jì)算值與試驗(yàn)值之間的絕對(duì)誤差在2.7～5.4m／s之間，相對(duì)誤差在4.48%～6%之間，二者吻合程度較好。

當(dāng)拋撒速度為139.9m／s，風(fēng)力條件為無(wú)風(fēng)時(shí)，對(duì)方程組（11）進(jìn)行求解，求解時(shí)間設(shè)定為20s，得出絲團(tuán)的飄落規(guī)律如圖7、圖8所示。絲團(tuán)空中全展后的絲束長(zhǎng)度有幾十m，其運(yùn)動(dòng)狀態(tài)為絲軸連帶著已展開(kāi)的絲束共同飄落。因此，絲軸的飄落狀態(tài)可以反映出絲束的飄落狀態(tài)，為了便于對(duì)絲團(tuán)的飄落規(guī)律進(jìn)行分析，下文中以絲軸的水平及豎直運(yùn)動(dòng)距離為研究對(duì)象，對(duì)其飄落規(guī)律進(jìn)行研究。

圖7 絲軸水平運(yùn)動(dòng)距離－時(shí)間關(guān)系

圖8 絲軸豎直運(yùn)動(dòng)距離－時(shí)間關(guān)系

通過(guò)圖7可以看出，絲軸在水平方向是以減速度狀態(tài)運(yùn)動(dòng)的，這是由空氣阻力及已展開(kāi)絲束對(duì)絲軸的拉力作用造成的；運(yùn)動(dòng)到10s左右時(shí)絲軸的水平運(yùn)動(dòng)距離變化越來(lái)越緩慢，說(shuō)明在10s左右時(shí)其水平方向的運(yùn)動(dòng)速度已經(jīng)衰減到了接近于0的狀態(tài)。在20s的時(shí)間內(nèi)，絲軸在水平方向共運(yùn)動(dòng)了92.3m。

通過(guò)圖8可以看出，在運(yùn)動(dòng)初期3～4s的時(shí)間內(nèi)，絲軸在豎直方向是以加速狀態(tài)運(yùn)動(dòng)的，這是由其自身重力大于豎直方向的空氣阻力所造成的；隨著豎直方向運(yùn)動(dòng)速度的不斷增大，絲軸所受空氣阻力也逐漸增大，最后，其自身重力與阻力達(dá)到平衡的狀態(tài)。所以，在4s過(guò)后，絲軸以接近勻速的狀態(tài)下落，下落速度為8.5～9m／s。在20s的時(shí)間內(nèi)，絲軸在豎直方向共下落了195.5m。

根據(jù)圖7、圖8所示內(nèi)容，繪制絲軸的運(yùn)動(dòng)軌跡，如圖9所示，圖9所示內(nèi)容與試驗(yàn)時(shí)觀測(cè)到的絲軸運(yùn)動(dòng)軌跡吻合。

圖9 絲軸空中運(yùn)動(dòng)軌跡

4 結(jié)論

①通過(guò)對(duì)導(dǎo)電纖維絲團(tuán)空中展開(kāi)與飄落過(guò)程中的受力情況進(jìn)行分析，建立了基于一定假設(shè)條件下的絲團(tuán)空中運(yùn)動(dòng)外彈道模型。此模型可以對(duì)絲團(tuán)在不同拋撒條件（風(fēng)速、拋撒絲團(tuán)高度、拋撒速度等）下的展開(kāi)及飄落情況進(jìn)行預(yù)測(cè)。

②為了驗(yàn)證模型的合理性，進(jìn)行了單個(gè)絲團(tuán)的定向拋撒試驗(yàn)及靶場(chǎng)飛行試驗(yàn)，計(jì)算結(jié)果與試驗(yàn)結(jié)果吻合較好，表明本文所建立的絲團(tuán)彈道模型合理。

［1］COSTELLO M F，F(xiàn)ROST G W.Simulation of projectiles connected by a flexible tether，ARL－CR－456［R］.2000.

［2］張國(guó)偉，馮順山，孫學(xué)清，等.導(dǎo)電絲束空中展開(kāi)規(guī)律［J］.彈道學(xué)報(bào)，2006，18（1）：48－50.ZHANG Guo－wei，F(xiàn)ENG Shun－shan，SUN Xue－qing，et al.The spreading law of conductive fiber in the air［J］.Journal of Ballistics，2006，18（1）：48－50.（in Chinese）

［3］劉顯江，孫旭，趙非玉，等.美軍CBU－94型斷電彈藥投放落點(diǎn)分析［J］.光電技術(shù)應(yīng)用，2007，22（3）：21－23.LIU Xian－jiang，SUN Xu，ZHAO Fei－yu，et al.Analysis of the point fall of CBU－94blackout bomb［J］.Electro－Optic Technology Application，2007，22（3）：21－23.（in Chinese）

［4］張德豐.MATLAB數(shù)值計(jì)算方法［M］.北京：機(jī)械工業(yè)出版社，2010：225－229.ZHANG De－feng.Numerical computation method of MATLAB［M］.Beijing：China Machine Press，2010：225－229.（in Chinese）

［5］浦發(fā).外彈道學(xué)［M］.北京：國(guó)防工業(yè)出版社，1980：21－22.PU Fa.External ballistics［M］.Beijing：National Defense Industry Press，1980：21－22.（in Chinese）

［6］楊啟仁.子母彈飛行動(dòng)力學(xué)［M］.北京：國(guó)防工業(yè)出版社，1999：20－24.YANG Qi－ren.Flight dynamics of cluster bomb［M］.Beijing：National Defense Industry Press，1999：20－24.（in Chinese）

［7］梅鳳祥，周際平，水小平.工程力學(xué)（下冊(cè)）［M］.北京：高等教育出版社，2003：412－416.MEI Feng－xiang，ZHOU Ji－ping，SHUI Xiao－ping.Engineering mechanics（Volume Two）［M］.Beijing：High Education Press，2003：412－416.（in Chinese）