顧 林，張 合，石秀麗

（南京理工大學 智能彈藥技術(shù)國防重點實驗室，江蘇 南京，210094）

電磁脈沖彈、磁爆加載彈藥等磁爆彈藥的能量和動力都是由彈載磁通量壓縮發(fā)生器提供的。磁通量壓縮發(fā)生器往往采用高速斷路開關(guān)將其電感儲存的能量加載到負載上[1]。高速斷路開關(guān)和磁通量壓縮發(fā)生器的聯(lián)合運行，可使得能量儲存階段和能量提取階段分割開來，并保障了阻抗匹配和脈沖成形。電爆炸絲斷路開關(guān)（Electro Explosive Breaker，EEB）的通斷時間短、斷路電流大、功率大、成本低，故在磁通量壓縮器脈沖形成網(wǎng)絡(luò)中應(yīng)用最為廣泛[2]。

電爆炸斷路開關(guān)是一次性、破壞性開關(guān)，數(shù)值計算、模擬仿真對其設(shè)計有著重要作用。解放軍理工大學工程兵工程學院朱翼超等人采用電阻率-比作用量模型，對含電爆炸絲元件的電路仿真模擬了整個電爆炸過程[3]。中國工程物理研究院蔣吉昊等人對電爆炸絲進行了一維磁流體模型數(shù)值模擬[4]。國防科技大學楊漢武等人則建立了PSpice模型進行了斷流分析[2]。英國 Loughborough大學開展了脈沖變壓器結(jié)合電爆炸絲斷路開關(guān)技術(shù)的研究[5]。這些建立的模型盡管能很好地解釋和符合實驗數(shù)據(jù)，但模型建立復雜，調(diào)試參數(shù)不方便，不利于開關(guān)設(shè)計。SIMULINK仿真模型相比較于數(shù)值模擬直觀，無需書寫大量程序公式函數(shù)，同時也避免了數(shù)值模擬計算的復雜性和不通用性，故本文利用SIMULINK軟件基于時變電阻模型，對并聯(lián)電爆炸絲斷路開關(guān)斷流特性進行仿真分析，研究電爆炸絲的長度、根數(shù)和直徑等參數(shù)變化對斷路時間以及斷路電流的影響。

1 電爆炸絲斷路開關(guān)物理模型

電爆炸絲斷路開關(guān)由許多根金屬細絲并聯(lián)組成，并聯(lián)的電爆炸絲快速通過大電流時，金屬絲先后經(jīng)歷固態(tài)、熔化、液態(tài)加熱、汽化、電弧成長及爆炸等幾個過程，見圖1。爆炸斷路開關(guān)斷路初期，電阻率隨著注入能量的增大而增大。當注入能量足夠大時，金屬導體材料解體，氣化波陣面向?qū)w內(nèi)部運動，并發(fā)生爆炸，從而開關(guān)斷開[6]。

圖1 典型電爆炸絲爆炸物理過程[6]Fig.1 Typical physical process of electrical explosion wires [6]

式（1）中：d為電爆炸絲直徑， (0)R為電阻初始值，vv為氣化波陣面的運動速度。式（1）的計算模型表明氣化波在銅爆炸絲中傳播速度vv和銅絲橫截面積決定著電爆炸斷路開關(guān)的電流通斷時間。但是vv由于計算比較復雜，不適合做數(shù)值模擬。本文從Tucker模型來推導電爆炸絲的電阻率具體計算模型。

忽略熱傳導，由焦耳定律[7]：

式（2）中：γ為爆炸金屬絲的能量密度；ρ為爆炸金屬絲的電阻率；β為輸入的電流密度。

爆炸金屬絲的能量密度和爆炸金屬絲的電阻率兩者呈一定比例的關(guān)系：

式（3）中：下標τ=1,2，分別對應(yīng)固態(tài)與液態(tài)；ρτ為最初電阻率；ετ為加熱的效率。

由式（2）～（3）推導得出：

在固態(tài)加熱結(jié)束時刻與液態(tài)加熱結(jié)束時刻，式（4）可寫為：

式（5）中：gmax為固態(tài)加熱結(jié)束時刻的比作用量；ρmax為固態(tài)加熱結(jié)束時刻的電阻率。又有：

將式（6）代入式（5）可得：

爆炸金屬導體在經(jīng)歷固態(tài)轉(zhuǎn)變?yōu)橐簯B(tài)、液態(tài)轉(zhuǎn)變?yōu)槠麘B(tài)的過程中能量的表達式為：

式（8）～（9）中：K為相變潛熱；S為金屬導體的截面積；m為金屬導體原初的相質(zhì)量；R為相1與相2的總電阻；R1、R2為相1與相2電阻。結(jié)合式（7）得出電阻率與比作用量的函數(shù)關(guān)系[3,8,9]：

圖2為文獻[6]的比作用量和電阻率的實驗數(shù)據(jù)，反映了ρ與g的大致函數(shù)關(guān)系，可以作為數(shù)值模擬仿真模型建立的參考。

利用電容器放電電路對電爆炸絲的斷流特性進行分析，電容放電回路電爆炸絲起爆時間et的經(jīng)驗公式[10]：

式（11）中：L為回路電感；C為電容值；有關(guān)kj、kr的計算文獻[3]中有詳細方法，這里不再贅述。

圖2 ρ與g實驗曲線(銅)Fig.2 Experimental curve of resistivity and specific action(Cu)

2 模型建立

采用Matlab軟件的SIMULINK模塊，建立電爆炸絲斷流仿真模型，部分如圖3所示。

圖3 SIMULINK電爆炸絲斷流仿真模型（部分）Fig.3 SIMULINK simulation model of EEB (parts)

根據(jù)上節(jié)的理論分析與推導的公式，利用積分模塊對電流平方進行積分，經(jīng)過增益得到相應(yīng)的比作用量。銅電爆炸絲電阻故再次用增益模塊得到相應(yīng)的電阻。根據(jù)UIR= ，電壓源由其兩端的電流來控制。導通時，電容(SIMULINK中電容充電電壓值可直接預置，無需充電)通過儲能電感、回路電阻和電爆炸絲進行放電，簡化等效為RLC回路，研究電爆炸過程的電學特性。

模型可以通過切斷電壓和切斷電流波形直觀地觀察電爆炸絲的工作過程，從而可以分析電爆炸絲的切斷電流、切斷電壓、切斷時間以及切斷效果之間的關(guān)系。模型取電容10μF、放電預置電壓50kV、負載電感2.5μH。

3 斷路特性仿真與分析

利用模型分析長度l、根數(shù)n以及直徑d等參數(shù)單一變量對電爆炸絲的斷路時間和斷路電流等斷流特性的影響。首先固定電爆炸絲根數(shù)n=30、直徑d=0.05mm，改變長度l為100～500mm，研究電爆炸絲的斷流特性，仿真結(jié)果見圖4和表1。

圖4 長度對電爆炸絲斷流及電感兩端電壓影響Fig.4 Influence of length on EEB and inductor’s voltage

表1 長度l對斷流特性影響Tab.1 Influence of length on current breaking characteristics

從表1可知長度l的變化對切斷電流的時間以及切斷電流峰值影響并不是太顯著。圖4顯示，長度l影響了電爆炸絲斷路開關(guān)切斷電流的過程。長度l小于一定值時，內(nèi)部產(chǎn)生的高電壓會使導體擊穿，將不會切斷，出現(xiàn)后期不能徹底切斷電流，無法恢復到零值的現(xiàn)象，仿真斷流曲線準確地反映了這一點。儲能電感兩端的電壓峰值隨著電爆炸絲長度的增加，出現(xiàn)了變大再變小的變化曲線，從電壓角度看長度l存在一個最優(yōu)值。

固定直徑d=0.05mm和長度l=500mm，改變根數(shù)n為30～90根，研究電爆炸絲的斷流特性，仿真結(jié)果見表2和圖5。

表2 根數(shù)對斷流特性的影響Tab.2 Influence of amount on current breaking characteristics

圖5 根數(shù)n對電爆炸絲斷流及電感兩端電壓的影響Fig.5 Influence of amount on EEB and inductor’s voltage

由表2及圖5可以看出，電爆炸絲根數(shù)n與切斷時間近似存在線性關(guān)系，即電爆炸絲總橫截面積與切斷時間存在線性關(guān)系。n變多，電爆炸時間變長，且切斷的電流值也變大。根數(shù)n變多，總橫截面積增大，電流密度減小，達到爆炸所需的比作用量g的時間變長。當n大于一定數(shù)量時，電流密度太小，所供給的能量將不足以使電爆炸絲被切斷。同時，跟長度l一樣，隨著電爆炸絲根數(shù)n的增多，儲能電感兩端的電壓峰值出現(xiàn)了變大再變小的變化曲線。

固定根數(shù)n=30和長度l=500mm，改變直徑d為0.05～0.10mm，研究電爆炸絲的斷流特性，仿真結(jié)果見表3和圖6。

表3 直徑d對斷流特性影響Tab.3 Influence of diameter on current breaking characteristics

圖6 直徑d對電爆炸絲斷流及電感兩端電壓影響Fig.6 Influence of diameter on EEB and inductor’s voltage

由表3和圖6可見，隨著電爆炸絲直徑d的增大，切斷電流的時間變長，切斷的電流值也變大。儲能電感兩端的電壓峰值隨著電爆炸絲直徑d的增大，先變大后變小。增大電爆炸絲直徑d與電爆炸絲根數(shù)n的效果很接近，都等效增大了電爆炸絲的總橫截面積。但實驗表明[8]，在總橫截面積一定時，選擇n大d小的效果較好，因為這樣的設(shè)計具有較低的電感[2,8]。因此可以對模型進行進一步的優(yōu)化，增加電爆炸絲電感分布的不均勻性與隨機性。

4 結(jié)束語

本文為了研究銅導體做并聯(lián)電爆炸絲斷路開關(guān)的斷流特性，首先理論分析了電爆炸絲斷路開關(guān)物理模型，在這基礎(chǔ)上采用 SIMULINK建立了基于時變電阻模塊的仿真模型。模型采用電容器脈沖放電對銅電爆炸絲斷路特性進行分析，得出電爆炸絲根數(shù)、長度以及直徑參數(shù)對斷路時間和斷路電流的影響規(guī)律。與其他模型相比較，基于 SIMULINK的電爆炸絲斷流仿真模型更為直觀便捷。仿真結(jié)果與文獻[2]基于PSpice的仿真模型結(jié)果對應(yīng)，根數(shù)、長度以及直徑等參數(shù)對斷流特性的影響規(guī)律也一致；并且也能很好地吻合文獻[11]的實驗結(jié)果，故此仿真模擬能夠準確反映電爆炸絲的一般規(guī)律。本文的理論分析推導和仿真研究為設(shè)計并聯(lián)電爆炸絲斷路開關(guān)的參數(shù)選擇提供了參考依據(jù)。

[1]孫承緯,周之奎,譯.磁通量壓縮發(fā)生器[M].北京：國防工業(yè)出版社,2008.

[2]楊漢武,鐘輝. PSpice模型用于電爆炸絲的數(shù)值模擬[J].國防科技大學學報,2000(22):38-42.

[3]朱翼超,羅根新,方向,陳宇.電爆炸絲的電氣特性[J].爆炸與沖擊,2011,31(6):664-668.

[4]蔣吉昊,楊宇,胡熙靜,王鋼華.電爆炸絲 1維磁流體模型數(shù)值模擬[J].強激光與粒子束,2006,18(3):463-466.

[5]B.M.Novac, I.R.Smith, K.Gregory. High-power high-voltage generator with a multi-pulse twin output[C]//13th IEEE Internatinal Pulse Power Conference. Las Vegas :IEEE Press,2001.

[6]龔興根.電爆炸斷路開關(guān)[J].強激光與粒子束,2002,14(4):577-582.

[7]Tucker T J，Toth R P. A . Computer code for the prediction of the behavior of electrical circuits containing exploding wire element[R].Albuquerque: Sandia Laboratories,1975.

[8]石秀麗.磁通量壓縮發(fā)生器關(guān)鍵技術(shù)研究[D].南京:南京理工大學,2011.

[9]楊漢武. 爆炸磁壓縮發(fā)生器及其脈沖功率調(diào)制研究[D].長沙:國防科學技術(shù)大學, 2002.

[10]曾正中. 實用脈沖功率技術(shù)引論[M].西安：陜西科學技術(shù)出版社,2003.

[11]Howard S. Leoplod. Effect of bridgewire parameters on explosive initiation[J]. Exploding Wires, 1964(3):125-152.