薄片炸藥與固體靶沖量耦合的計算模型＊

盧 強(qiáng)，王占江，劉曉新，郭志昀，吳玉蛟

（西北核技術(shù)研究所，陜西西安710024）

薄片炸藥與固體靶沖量耦合的計算模型＊

盧 強(qiáng)，王占江，劉曉新，郭志昀，吳玉蛟

（西北核技術(shù)研究所，陜西西安710024）

薄片炸藥可用于模擬強(qiáng)脈沖X射線輻照下空間結(jié)構(gòu)的響應(yīng)，其與固體靶的沖量耦合作用是進(jìn)行模擬載荷設(shè)計的基礎(chǔ)。本文中基于有限元模擬技術(shù)，研究不同直徑薄片炸藥的邊界稀疏作用和安裝支座的反射作用對薄片炸藥與固體靶沖量耦合的影響，結(jié)果表明：邊界稀疏作用對薄片炸藥的比沖量起到減弱的作用，由邊界稀疏作用引起的薄片炸藥比沖量的變化量和薄片炸藥直徑的倒數(shù)近似成線性關(guān)系；薄片炸藥安裝支座的反射作用對薄片炸藥的比沖量起到增強(qiáng)的作用，由支座反射作用引起的薄片炸藥比沖量的變化量和薄片炸藥直徑的平方近似成反比；隨著薄片炸藥直徑的增加，薄片炸藥比沖量趨于定值，邊界稀疏作用和支座反射作用可以忽略；存在一個臨界直徑，當(dāng)薄片炸藥的直徑為臨界直徑時，薄片炸藥的邊界稀疏作用和支座反射作用達(dá)到平衡，其比沖量大小和薄片炸藥的直徑為無窮大時的比沖量值相等。

爆炸力學(xué)；沖量耦合；比沖量；薄片炸藥

高空核爆炸約有70%～80%的能量以X射線的形式釋放，當(dāng)強(qiáng)脈沖X射線輻照到空間結(jié)構(gòu)體表面時，在空間結(jié)構(gòu)體迎光面表層產(chǎn)生能量沉積，并引起表層材料熔化、汽化，從而引起物質(zhì)的反沖噴射，形成作用在空間結(jié)構(gòu)體上的汽化反沖沖量，繼而引起空間結(jié)構(gòu)的動態(tài)響應(yīng)、破壞或失效，這屬于脈沖X射線對空間結(jié)構(gòu)的力學(xué)效應(yīng)［1－2］。在實(shí)驗(yàn)室內(nèi)可以采用脈沖電子束輻照、激光輻照等模擬脈沖X射線的輻照效應(yīng)［2－3］，采用光敏炸藥、薄片炸藥或柔爆索等模擬脈沖X射線引起空間飛行器的結(jié)構(gòu)響應(yīng)［3］。

R.A.Benham［4］和W.G.Rivera等［5－6］發(fā)展了一種通過閃光起爆光敏炸藥的加載技術(shù)，可以實(shí)現(xiàn)大面陣同時加載，能較好地模擬X射線輻照空間結(jié)構(gòu)殼體的脈沖載荷特征，其技術(shù)難點(diǎn)在于光敏炸藥的安全控制以及噴涂工藝。W.E.Alzheimer等［7］、B.W.Franklin［8］、H.E.Lindberg等［9］和S.W.Kirkpatrick等［10］通過對按照一定規(guī)律離散的薄片炸藥條一端起爆的方式實(shí)現(xiàn)對結(jié)構(gòu)的滑移爆轟加載，從而模擬脈沖X射線輻照下的載荷。毛勇建等［11－13］對薄片炸藥加載技術(shù)進(jìn)行了較為系統(tǒng)的研究，從數(shù)值模擬的角度分析了薄片炸藥加載的載荷特征以及對X射線力學(xué)效應(yīng)的等效性問題。相比于光敏炸藥加載技術(shù)，薄片炸藥加載技術(shù)的優(yōu)點(diǎn)是安全性高、技術(shù)成本較低，但由于炸藥的極限爆轟尺寸的限制，較難實(shí)現(xiàn)對空間結(jié)構(gòu)載荷均勻性高、低比沖量的加載模擬。趙國民等［14］采用按照一定規(guī)律離散分布的鉛皮柔爆索一端起爆的方式來模擬脈沖X射線輻照下的載荷。此技術(shù)主要是利用柔爆索爆炸產(chǎn)生的飛散鉛粒子與空間結(jié)構(gòu)碰撞產(chǎn)生的相互作用，雖然載荷的可設(shè)計性強(qiáng)，但同樣較難實(shí)現(xiàn)對空間結(jié)構(gòu)載荷均勻性高、低比沖量的加載模擬。

無論采用何種方式來模擬X射線力學(xué)效應(yīng)，載荷的比沖量特征都是針對相應(yīng)空間結(jié)構(gòu)進(jìn)行加載設(shè)計的關(guān)鍵技術(shù)參數(shù)。目前國內(nèi)對X射線力學(xué)效應(yīng)的模擬仍主要采用薄片炸藥加載來實(shí)現(xiàn)［11－13，15－17］。采用薄片炸藥模擬脈沖X射線引起空間飛行器結(jié)構(gòu)響應(yīng)時，需通過一定方式對薄片炸藥進(jìn)行離散以滿足模擬脈沖X射線輻照空間結(jié)構(gòu)產(chǎn)生的反沖沖量的需要。對薄片炸藥的離散加載設(shè)計的前提是要對薄片炸藥與固體靶的沖量耦合作用進(jìn)行研究。岳曉紅等［15］采用薄片炸藥與固體靶的沖量耦合作用測量的擺錘結(jié)構(gòu)及薄片炸藥比沖量標(biāo)定方法，但起爆炸藥條、薄片炸藥邊界的稀疏作用對比沖量的影響沒有提及。林鵬等［16］采用直線運(yùn)動式?jīng)_量探頭測得了0.4mm厚薄片炸藥的比沖量，實(shí)驗(yàn)中薄片炸藥由柔爆索引爆，并未考慮起爆方式對薄片炸藥比沖量的影響，在數(shù)值模擬研究中沒有考慮薄片炸藥邊界稀疏作用。郭志昀等［17］對0.4和0.5mm厚薄片炸藥與固體靶的沖量耦合特性進(jìn)行研究，指出薄片炸藥邊界稀疏作用和安裝支座反射沖量對比沖量測試具有影響，但未進(jìn)一步討論。

本文中基于有限元模擬技術(shù)，考慮不同直徑薄片炸藥邊界稀疏以及安裝支座反射沖量對超薄片狀炸藥與固體靶的沖量耦合作用的影響，建立了薄片炸藥與固體靶沖量耦合的計算分析模型。

1 薄片炸藥比沖量測試技術(shù)

圖1 薄片炸藥沖量測量裝置Fig.1 Impulse measurement apparatus for sheet explosive

2 薄片炸藥與固體靶沖量耦合模型

薄片炸藥與固體靶沖量耦合作用測試中不同因素對沖量耦合作用的影響主要有：（1）薄片炸藥的尺寸效應(yīng)。薄片炸藥沖量測試的理想實(shí)驗(yàn)條件是應(yīng)滿足一維應(yīng)變式的加載，實(shí)際情況下難以做到，一般采用具有一定大小的薄片炸藥比沖量來代替理想的薄片炸藥比沖量，但有限尺寸薄片炸藥邊界產(chǎn)生的稀疏波會減弱薄片炸藥的加載能力，使得測試得到的薄片炸藥比沖量偏?。唬?）由于柔爆索和薄片炸藥都是柔性加載元件，一般需采用1個連接支座過渡，如圖2所示。使柔爆索頂端可以垂直頂在薄片炸藥表面，以提高薄片炸藥起爆的可靠性。薄片炸藥起爆后會在連接支座和薄片炸藥之間產(chǎn)生1個較高的反射壓，同時柔爆索爆炸也會對固體靶產(chǎn)生1個附加沖量，這會使得測試得到的薄片炸藥比沖量偏大。基于上述考慮，采用有限元方法，以不同尺寸的連接支座及不同厚度、不同直徑的薄片炸藥進(jìn)行數(shù)值模擬，分析薄片炸藥邊界稀疏作用和連接支座沖量反射作用對薄片炸藥與固體靶沖量耦合的影響，由此建立薄片炸藥比沖量特性的計算模型。

圖2 柔爆索起爆薄片炸藥的支撐結(jié)構(gòu)Fig.2 Structural support for sheet explosive detonated by mild detonating fuse

圖3 爆索柔起爆薄片炸藥的支撐結(jié)構(gòu)的簡化模型Fig.3 Simplified model of the structural support for sheet explosive detonated by mild detonating fuse

薄片炸藥加載沖擊擺的擺錘過程約為10～100μs量級，之后擺錘即可達(dá)到恒定速度，而沖擊擺的周期為1s量級，可以認(rèn)為薄片炸藥加載擺錘質(zhì)量塊的過程瞬間完成，因此可對沖擊擺模型進(jìn)行簡化，只考慮薄片炸藥加載質(zhì)量塊的過程，忽略了接觸爆炸時質(zhì)量塊產(chǎn)生的塑性變形對沖量的影響。另外柔爆索附加沖量的影響在數(shù)值模擬中先不予考慮，在后面分析中再談及這個問題。由此圖2所示的柔爆索起爆薄片炸藥的支撐結(jié)構(gòu)可以簡化為圖3所示的模型。由于結(jié)構(gòu)的對稱性，建立1／4有限元模型進(jìn)行計算，柔爆索和薄片炸藥之間的連接座采用有機(jī)玻璃材料，質(zhì)量塊采用鋼材料，如圖4所示。所有模擬中均采用相同的網(wǎng)格尺寸，其中薄片炸藥沿厚度方向的網(wǎng)格尺寸為0.1mm。

2.1 邊界稀疏波對沖量耦合作用的影響

為考慮薄片炸藥與固體靶沖量耦合作用的影響，在計算模型中未對連接支座建模，薄片炸藥直接貼敷在質(zhì)量塊表面，起爆點(diǎn)設(shè)置在圓形薄片炸藥上表面中心位置。選取了厚度δ為0.2、0.3和0.5mm的薄片炸藥進(jìn)行模擬，其直徑D為10～60mm。設(shè)質(zhì)量塊的質(zhì)量為m0，薄片炸藥爆炸后質(zhì)量塊的穩(wěn)定運(yùn)動速度為v，則薄片炸藥的比沖量i為：

圖4 薄片炸藥比沖量測試的有限元模型Fig.4 Finite model for measurement of specific impulse for sheet explosive

通過計算，圖5給出了僅考慮邊界稀疏作用下的薄片炸藥比沖量iRA的模擬結(jié)果。從圖5中可以看出，對應(yīng)同一直徑的薄片炸藥，其厚度越大，比沖量也越大。對應(yīng)同一厚度的薄片炸藥，隨著直徑的增加，比沖量大小趨于恒定值。以1／D 為橫坐標(biāo)、比沖量iRA為縱坐標(biāo)對圖5進(jìn)行了重新繪制，如圖6所示。圖6中不同形狀的點(diǎn)代表不同厚度薄片炸藥比沖量數(shù)值模擬結(jié)果，直線為對應(yīng)的比沖量擬合線?？梢钥闯?，通過模擬計算給出的薄片炸藥比沖量iRA和1／D近似成線性關(guān)系（3種厚度薄片炸藥比沖量的擬合線性度均在0.996以上），不同厚度薄片炸藥比沖量擬合直線的截距代表了對應(yīng)厚度薄片炸藥的極限比沖量值（D 趨于無窮大時的比沖量）。設(shè)極限比沖量值為i0（δ），則圖6中的擬合直線可以統(tǒng)一寫為：

式中：斜率k（δ）為與薄片炸藥厚度相關(guān)的量，其表示薄片炸藥比沖量隨著1／D變化的快慢，k（δ）／D代表了稀疏作用引入到薄片炸藥比沖量中的大小。從圖6可以看出，隨著薄片炸藥厚度 增加而增加，這說明薄片炸藥厚度越大，邊界稀疏作用的影響也越大。

圖5 考慮邊界稀疏作用下的薄片炸藥比沖量iRA隨D的變化Fig.5 Specific impulse iRAfor sheet explosive varying with Din consideration of the rarefaction by the border

圖6 考慮邊界稀疏作用下的薄片炸藥比沖量iRA隨1／D的變化Fig.6 Specific impulse iRAfor sheet explosive varying with 1／Din consideration of the rarefaction by the border

2.2 連接支座對薄片炸藥與固體靶沖量耦合作用的影響

要想防止和解決物流企業(yè)財務(wù)風(fēng)險，就必須做好建立、健全企業(yè)內(nèi)部財務(wù)控制機(jī)制體系的有效措施。物流企業(yè)財務(wù)內(nèi)部監(jiān)控工作建立在企業(yè)各項資金預(yù)算的基礎(chǔ)上，并且符合企業(yè)運(yùn)作的整體利益。企業(yè)在充分運(yùn)用實(shí)時財務(wù)信息的條件下，應(yīng)該及時監(jiān)督和控制公司整體經(jīng)營目標(biāo)進(jìn)展情況。同時，企業(yè)需要加強(qiáng)對資金流量加以監(jiān)控，其目的為確保企業(yè)經(jīng)營目標(biāo)實(shí)現(xiàn)的情況下提升資金使用的效率。

在實(shí)際應(yīng)用中，薄片炸藥和固體靶之間需用支座進(jìn)行連接，安裝示意圖如圖2所示。柔爆索穿過有機(jī)玻璃支座中心孔直接頂?shù)奖∑ㄋ幧?，柔爆索直徑?mm，裝藥（RDX）的線密度為0.2g／m，柔爆索引爆薄片炸藥可近似為點(diǎn)起爆。在數(shù)值模擬中忽略了柔爆索附加沖量的作用，把有機(jī)玻璃支座和薄片炸藥接觸面的中心直接設(shè)置為起爆點(diǎn)，如圖3所示。為對比說明有機(jī)玻璃支座對薄片炸藥比沖量的影響，采用表1所示的幾種工況進(jìn)行數(shù)值模擬。

表1 考慮連接支座反射影響時的模擬工況Table 1 Working condition for simulation in consideration of the reflection by the structural support

圖7給出了考慮支座反射和邊界稀疏聯(lián)合作用時和僅考慮邊界稀疏作用時的薄片炸藥比沖量模擬結(jié)果，分別用虛線和實(shí)線表示。用iRE表示考慮了有機(jī)玻璃支座反射后薄片炸藥的比沖量。當(dāng)考慮支座反射沖量的影響后，iRE大于iRA，并且iRE隨D的增加而變小并趨于恒定值，這和僅考慮邊界稀疏作用時給出的iRA的變化規(guī)律是相反的。值得注意的是隨著D的增加，iRE和iRA越來越接近，這說明隨著D的增加，無論是薄片炸藥的邊界稀疏作用還是有機(jī)玻璃支座的反射作用都對薄片炸藥比沖量測試的影響在減小。定義Δi為：

式中：Δi反映了扣除稀疏作用后有機(jī)玻璃支座反射引起比沖量的增加量，d和h分別為有機(jī)玻璃支座的直徑和厚度。如圖7中的灰色區(qū)域所示。有機(jī)玻璃支座引入薄片炸藥與固體靶沖量耦合作用中的總沖量增量設(shè)為ΔI，則ΔI可以寫為：

考慮有機(jī)玻璃支座反射影響時，測試得到的總沖量可寫為：

根據(jù)式（4）～（5），圖8給出了總沖量I和有機(jī)玻璃支座引入薄片炸藥與固體靶沖量耦合作用中的總沖量增量ΔI隨薄片炸藥直徑的變化規(guī)律?？梢钥闯鲭S著D的增加，ΔI變化不大，可近似視為ΔI不隨D變化，即ΔI＝ΔI（δ，d，h）。結(jié)合式（4），能夠得到Δi隨著D的變化滿足冪指數(shù)衰減，且衰減指數(shù)為2。由于薄片炸藥的邊界稀疏作用和有機(jī)玻璃支座的反射作用是相互耦合的，同時有機(jī)玻璃支座自身的邊界也引入了一定的稀疏作用，通過式（3）計算的Δi并不能完全扣除邊界稀疏作用造成的影響。若考慮附加柔爆索沖量的影響，可認(rèn)為近似為點(diǎn)起爆的柔爆索爆炸對質(zhì)量塊的附加沖量為定值，可以附加到ΔI中去考慮，其引起的薄片炸藥比沖量的變化符合式（4）的描述。

圖7 考慮支座反射作用下的比沖量i隨D的變化關(guān)系Fig.7 Relationship between the specific impulse i and D in consideration of the reflection by the structural support

圖8 考慮支座反射作用下的沖量I隨D的變化關(guān)系Fig.8 Relationship between the total impulse I and D in consideration of the reflection by the structural support

2.3 薄片炸藥與固體靶沖量耦合的計算模型及分析

從上述數(shù)值模擬結(jié)果的分析可以知道，隨著D的增加，邊界稀疏及有機(jī)玻璃支座反射作用對薄片炸藥與固體靶沖量耦合的影響會減弱，假定邊界稀疏和支座反射作用對沖量耦合作用的影響相互獨(dú)立，結(jié)合式（2）～（4）能夠給出考慮有機(jī)玻璃支座時的比沖量i計算公式：

圖9 比沖量計算模型的擬合曲線和數(shù)值模擬對比Fig.9 Comparison of specific impulse between model fitting results and simulation

現(xiàn)在討論比沖量i計算模型的適用性，具體結(jié)果可見如表2所示。

表2 模型擬合結(jié)果和數(shù)值模擬結(jié)果對比Table 2 Comparison between model fitting results and numerical simulation results

對于δ＝0.2、0.3、0.5mm這3種情況，采用計算模型擬合給出的薄片炸藥比沖量近似值i0和i0（δ）十分接近，兩者的最大相對誤差為0.71%。通過計算模型的擬合曲線給出的臨界直徑近似值和由數(shù)值模擬結(jié)果插值給出i0（δ）對應(yīng)的臨界直徑D0均隨著薄片炸藥厚度的增加而變小，但在δ＝0.2，0.3時計算模型給出的誤差較大，在δ＝0.5mm時兩者給出的結(jié)果較為接近。這主要是因?yàn)棣模?0.5mm時，式（6）表示的關(guān)于1／D的拋物線對稱軸在2D0＝29.2mm處，當(dāng)D的模擬區(qū)間為10～60mm時，計算結(jié)果充分反映了i（δ，D0，d，h）隨D先減小而后增大的趨勢。當(dāng)δ＝0.2，0.3mm時，2D0＝103.2，65.0mm，當(dāng)D的模擬區(qū)間為10～60mm時，計算結(jié)果只能反映i（δ，D0，d，h）隨D減小的趨勢，這使得采用計算模型擬合的比沖量曲線誤差較大。可見，若能確定D0，則在進(jìn)行比沖量測試實(shí)驗(yàn)時僅用1發(fā)實(shí)驗(yàn)即可確定薄片炸藥的比沖量i0（δ）。但在實(shí)際應(yīng)用中，對一種新配方或新工藝生產(chǎn)的薄片炸藥而言，是無法提前確定一個合適的D0進(jìn)行實(shí)驗(yàn)的，那么為了測得某一厚度為δ的薄片炸藥的比沖量i0（δ），則需要進(jìn)行一系列不同直徑薄片炸藥與固體靶沖量耦合作用的測試實(shí)驗(yàn)，以得到反映薄片炸藥比沖量隨D先減小而后增大現(xiàn)象的測試曲線。

2.4 連接支座直徑大小對臨界直徑的影響

從上述分析可以看出，臨界直徑D0越小越有利于獲得薄片炸藥的極限比沖量i0（δ），因此需要通過一定的設(shè)計來減小臨界直徑D0。由于，則減小D0有2種方式：（1）增加k（δ）；（2）降低A（δ，d，h）。對于特定形狀的薄片炸藥而言，k（δ）只和薄片炸藥的厚度有關(guān)，其表征的是特定厚度薄片炸藥的邊界稀疏效應(yīng)，薄片炸藥厚度一定時，k（δ）可視為常數(shù)，因此只能通過降低表征反射沖量大小的A（δ，d，h）來減少臨界直徑k（δ）的大小。在采用7mm×0.5mm尺寸的有機(jī)玻璃支座進(jìn)行薄片炸藥與固體靶沖量耦合作用模擬的基礎(chǔ)上，增加了5mm×0.5mm和10mm×0.5mm這2種尺寸有機(jī)玻璃支座的模擬，薄片炸藥厚度取0.5mm。

圖10給出了采用不同尺寸有機(jī)玻璃支座作用時薄片炸藥與固體靶沖量耦合作用的模擬結(jié)果?？梢钥闯?，隨著有機(jī)玻璃支座直徑d的減少，薄片炸藥比沖量值隨之減少，這說明支座反射沖量對薄片炸藥比沖量的影響降低，反映在式（6）中即A（δ，d，h）也隨之降低。當(dāng)d＝7，10mm時，對應(yīng)的臨界直徑D0為14.6和35.7mm。當(dāng)d＝5mm時，在D為10～60mm的模擬區(qū)間內(nèi)得到的薄片炸藥比沖量均小于i0（δ）＝1 770.39Pa·s，則可知D0必小于10mm。利用此組數(shù)據(jù)結(jié)合式（6）可以給出D≈7mm，而通過數(shù)值模擬得到的D≈5.7 mm，造成此種現(xiàn)象的主要原因是當(dāng)薄片炸藥直徑和有機(jī)玻璃支座的直徑接近時，薄片炸藥邊界稀疏作用和支座反射作用耦合作用增大，邊界稀疏作用和支座反射作用對沖量耦合作用的影響相互獨(dú)立的假定不再成立，這造成了計算模型的偏差。上述分析結(jié)果證實(shí)了通過降低有機(jī)玻璃支座的直徑d的大小可以有效降低薄片炸藥比沖量測試中臨界直徑D0的大小。

圖10 不同尺寸有機(jī)玻璃支座作用下比沖量曲線Fig.10 Specific impulses for sheet explosive with different PMMA support sizes

3 結(jié) 論

本文中針對薄片炸藥與固體靶沖量耦合特性進(jìn)行了研究，建立薄片炸藥比沖量特性的計算模型，該模型對于指導(dǎo)薄片炸藥比沖量的測定以及對模擬強(qiáng)脈沖X射線施加到空間結(jié)構(gòu)上的載荷進(jìn)行理論設(shè)計具有一定實(shí)用價值。基于分析可以得到以下結(jié)論：

（1）薄片炸藥邊界稀疏對沖量耦合作用起到減弱的效果，其對比沖量測試的影響隨著薄片炸藥直徑的增加而減少，邊界稀疏引起的薄片炸藥比沖量的減少量和1／D近似成線性關(guān)系；

（2）用于連接薄片炸藥和柔爆索的支座反射作用對沖量耦合起到增強(qiáng)的效果，其對比沖量測試的影響隨著薄片炸藥直徑的增加而減少，支座反射引起的薄片炸藥比沖量的增加量隨著薄片炸藥直徑的增加近似滿足冪指數(shù)衰減的變化規(guī)律；

（3）根據(jù)式（6）建立的薄片炸藥比沖量計算模型可以知道，薄片炸藥比沖量和1／D成拋物線關(guān)系，且存在一個薄片炸藥的臨界直徑，當(dāng)D＝D0時，有i（δ，D，d，h）＝i0（δ）；

（4）臨界直徑 較小時有利于獲得薄片炸藥的極限比沖量i0（δ），對于特定厚度和邊界形狀的薄片炸藥，降低連接支座反射沖量的大小可以有效降低臨界直徑D0。

［1］Longley R W.Analytical relationships for estimating the effects of X－rays on materials［R］.AD－786926，1974.

［2］周南，喬登江.脈沖束輻照材料動力學(xué)［M］.北京：國防工業(yè)出版社，2002.

［3］毛勇建，鄧宏見，何榮建.強(qiáng)脈沖軟X光噴射沖量的幾種模擬加載技術(shù)［J］.強(qiáng)度與環(huán)境，2003，30（2）：55－64.Mao Yongjian，Deng Hongjian，He Rongjian.Several simulation techniques of blow－off impulses by intense pulsed cold X－ray［J］.Structure and Environment Engineering，2003，30（2）：55－64.

［4］Benham R A.Preliminary experiments using light－initiated high explosive for driving thin flyer plates［R］.SAND79－1847，1980.

［5］Rivera W G，Benham R A，Duggins B D，et al.Explosive technique for impulse loading of space structures［C］∥The 70th Shock and Vibration Symposium.1999.

［6］Rivera W G.Cosine distributed impulsive loading of a thin aluminum ring by a light initiated high explosive driven flyer plate［C］∥The 77th SAVIAC Shock and Vibration Symposium.2006.

［7］Alzheimer W E，F(xiàn)orrestal M J.Response of a circular elastic shell to moving and simultaneous loads［J］.AIAA Journal，1970，8（5）：970－971.

［8］Franklin B W.Development of a high－energy flexible sheet explosive［R］.AD－786510，1974.

［9］Lindberg H E，Colton J D.Sheet explosive simulation for combined shock and structural response［R］.AFWL－TR－69－124，1970.

［10］Kirkpatrick S W，Hohmes B S.Structural response of thin cylindrical shells subjected to impulsive external loads［J］.AIAA Journal，1988，26（1）：96－103.

［11］毛勇建，李玉龍，陳穎，等.炸藥條加載圓柱殼的數(shù)值模擬（Ⅰ）：流固耦合模擬［J］.高壓物理學(xué)報，2012，26（2）：155－162.Mao Yongjian，Li Yulong，Chen Ying，et al.Numerical simulation of cylindrical shell loaded by explosive rods（Ⅰ）：Fluid－structure interaction simulation［J］.Chinese Journal of High Pressure Physics，2012，26（02）：155－162.

［12］毛勇建，李玉龍，陳穎，等.炸藥條加載圓柱殼的數(shù)值模擬（Ⅱ）：解耦分析與實(shí)驗(yàn)驗(yàn)證［J］.高壓物理學(xué)報，2013，27（1）：76－82.Mao Yongjian，Li Yulong，Chen Ying，et al.Numerical simulation of cylindrical shell loaded by explosive rods（Ⅱ）：Decoupling analysis and experimental validation［J］.Chinese Journal of High Pressure Physics，2013，27（1）：76－82.

［13］毛勇建，李玉龍，陳穎，等.炸藥條加載圓柱殼的數(shù)值模擬（Ⅲ）：對X射線力學(xué)效應(yīng)的模擬等效性分析［J］.高壓物理學(xué)報，2013，27（5）：711－718.Mao Yongjian，Li Yulong，Chen Ying，et al.Numerical simulation of cylindrical shell loaded by explosive rods（Ⅲ）：Fidelity for simulating X－ray mechanical effects［J］.Chinese Journal of High Pressure Physics，2013，27（5）：711－718.

［14］趙國民，王占江，張若棋.用柔爆索構(gòu)成沿圓柱殼體周向呈余弦分布的沖量載荷［J］.爆炸與沖擊，2008，28（6）：557－560.Zhao Guomin，Wang Zhanjiang，Zhang Ruoqi.To form a cosine distribution impulse around cylindrical shell circumference by using mild detonating fuse［J］.Explosion and Shock Waves，2008，28（6）：557－560.

［15］岳曉紅，毛勇建，何榮建，等.片狀炸藥的比沖量標(biāo)定技術(shù)：火炸藥技術(shù)及鈍感彈藥學(xué)術(shù)研討會［C］∥火炸藥技術(shù)及鈍感彈藥學(xué)術(shù)研討會.2002：396－398.

［16］林鵬，王長利，王等旺.撓性炸藥比沖量的數(shù)值模擬與實(shí)驗(yàn)研究［J］.火炸藥學(xué)報，2011，34（4）：30－33.Lin Peng，Wang Changli，Wang Dengwang.Numerical simulation and experimental studies on impulse of flexible explosive［J］.Chinese Journal of Explosives and Propellants，2011，34（4）：30－33.

［17］郭志昀，王占江，盧強(qiáng)，等.柔爆索引爆薄片炸藥的比沖量研究［J］.兵工學(xué)報，2013，34（增刊1）：182－186.Guo Zhiyun，Wang Zhanjiang，Lu Qiang，et al.Research on specific impulse of the sheet explosive ignited by mild detonating fuses［J］.Acta Armamentarii，2013，34（suppl 1）：182－186.

A computational model for impulse coupling between sheet explosive and target

Lu Qiang，Wang Zhanjiang，Liu Xiaoxin，Guo Zhiyun，Wu Yujiao
（Northwest Institute of Nuclear Technology，Xi’an 710024，Shaanxi，China）

The sheet explosive can be used to simulate the space structure response caused by the highpower pulse X－rays，and the impulse coupling of the sheet explosive and the target is the basis for designing the simulated load.Based on the numerical results by the nonlinear finite element code，the computational modeling for the impulse coupling of the sheet explosive and the target is formulated in consideration of the effect on the specific impulse of the sheet explosive caused by the rarefaction of the border of the sheet explosive and the reflection of the support for connecting the mild detonating fuse and the sheet explosive.The analytic results reveal that the increment of the specific impulse caused by the rarefaction of the border is approximately in inverse relation to the diameter of the sheet explosive，and that caused by the reflection of the support is approximately in inverse relation to the square of the diameter of the sheet explosive.The specific impulse for the sheet explosive converges to the limit specific impulse related to the depth of the sheet explosive as the diameter of the sheet explosive reaches infinity.Both the rarefaction for the border of the sheet explosive and the reflection of the support are negligible with the diameter of the sheet explosive increasing.There exists a critical diameter and when it is equaled by the diameter of the sheet explosive，the specific impulse for the sheet explosive equals the limit specific impulse as the diameter of the sheet explosive reaches infinity.

mechanics of explosion；impulse coupling；specific impulse；sheet explosive

O383國標(biāo)學(xué)科代碼：13035

A

10.11883／1001－1455（2017）01－0084－08

（責(zé)任編輯 王易難）

2015－09－10；

2016－01－18

盧 強(qiáng)（1984— ），男，博士研究生；通信作者：王占江，wangzhanjiang＠nint.ac.cn。