周廣宇，胡時(shí)勝

（1.山東水利職業(yè)學(xué)院，山東 日照 276826；2.中國(guó)科學(xué)技術(shù)大學(xué)中國(guó)科學(xué)院材料力學(xué)行為和設(shè)計(jì)重點(diǎn)實(shí)驗(yàn)室，安徽 合肥 230026）

現(xiàn)代高科技戰(zhàn)爭(zhēng)彈藥技術(shù)的發(fā)展，對(duì)火工品的安全性提出了更高的要求，而耐過(guò)載實(shí)驗(yàn)成為檢驗(yàn)火工品在沖擊環(huán)境下或經(jīng)沖擊后使用安全性的重要手段，構(gòu)成其加載環(huán)境的主要因素［1］包括沖擊的時(shí)間歷程（沖擊波形）、沖擊脈沖的峰值加速度和沖擊脈沖的持續(xù)時(shí)間（脈沖寬度）。目前，Hopkinson壓桿技術(shù)［2］已被廣泛用來(lái)評(píng)估火工品的抗過(guò)載能力、標(biāo)定高g 值加速度傳感器［3－5］。T.C.Togami等［4－5］用改進(jìn)后的Hopkinson壓桿裝置得到加速度脈沖的峰值200 000g（g＝9.81m／s2），并對(duì)比了不同子彈速度下的加速度曲線。R.D.Sill［6］、T.C.Togami等［4］、李玉龍等［7］及盛黨紅等［8］采用 Hopkinson壓桿技術(shù)標(biāo)定了高g值加速度傳感器。沈瑞琪等［9］則將Hopkinson壓桿技術(shù)應(yīng)用于火工品安全性的評(píng)估。雖然采用Hopkinson壓桿技術(shù)可以得到高于200 000g的沖擊加速度峰值，但沒(méi)有考慮波形彌散、加速度脈寬對(duì)火工品抗過(guò)載性能的影響，且加速度脈寬也只在10～25μs之間。鄧強(qiáng)等［10］將Hopkinson壓桿中的波形整形技術(shù)應(yīng)用到火工品安全性檢測(cè)中，試圖得到不同的沖擊脈沖。進(jìn)一步表明火工品的破壞臨界加速度值隨著加速度脈寬的增大而減小，研究火工品的抗過(guò)載性能時(shí)，外界環(huán)境的沖擊加速度值和加速度脈寬同時(shí)影響著火工品的安全性能。因此，研究波形整形技術(shù)對(duì)加速度值和加速度脈寬的影響規(guī)律，來(lái)實(shí)現(xiàn)更高的加速度值和更長(zhǎng)的加速度脈寬，是火工品安全測(cè)試中面臨的問(wèn)題，更是迫切需要解決的問(wèn)題。本文中，基于Hopkinson壓桿（HPB）中的波形整形技術(shù)，針對(duì)高g值加速度發(fā)生器火工品安全性檢測(cè)的要求，通過(guò)數(shù)值模擬，系統(tǒng)研究子彈及波形整形器等對(duì)加載脈沖的影響。

1 實(shí)驗(yàn)裝置

本文中，采用的實(shí)驗(yàn)裝置實(shí)際為一套Hopkinson壓桿裝置，如圖1所示。

圖1 高g值沖擊加速度發(fā)生器實(shí)驗(yàn)裝置示意圖Fig.1 Schematic diagram of high－g－value impact acceleration generator

高壓氣槍發(fā)射的子彈以一定的速度撞擊Hopkinson壓桿，在該壓桿中產(chǎn)生彈性應(yīng)力波。此應(yīng)力波沿壓桿并通過(guò)桿端傳遞到試件上，對(duì)試件進(jìn)行高速加載。由于Hopkinson壓桿長(zhǎng)徑比大，可將桿中的應(yīng)力波看作一維應(yīng)力波。根據(jù)一維應(yīng)力波理論及牛頓第二定律，得到壓桿端面試件的加速度［11］。

2 數(shù)值模擬

2.1 幾何模型

為了系統(tǒng)研究子彈、波形整形器等相關(guān)因素對(duì)火工品沖擊實(shí)驗(yàn)中沖擊脈沖的影響規(guī)律，本文中采用ABAQUS／Explicit有限元軟件［12］對(duì)高g值沖擊加速度發(fā)生器實(shí)驗(yàn)裝置產(chǎn)生高g值加速度沖擊脈沖的工作過(guò)程進(jìn)行了數(shù)值模擬。

因?yàn)檎麄€(gè)實(shí)驗(yàn)裝置以及載荷等都關(guān)于桿的中心線軸對(duì)稱(chēng)，所以可以使用軸對(duì)稱(chēng)單元進(jìn)行二維分析。如圖2所示，單元類(lèi)型均設(shè)為CAX4R軸對(duì)稱(chēng)單元，子彈與整形器、整形器與入射桿之間采用面－面的有限滑移接觸，接觸屬性設(shè)定為無(wú)摩擦硬接觸。為避免單元形狀發(fā)生畸變、計(jì)算無(wú)法收斂并控制脈沖波形質(zhì)量，對(duì)整形器采用ALE自適應(yīng)網(wǎng)格。

圖2 高g值沖擊加速度發(fā)生器有限元模型Fig.2 The finite element model for high－g－value impact acceleration generator

2.2 材料模型

子彈、壓桿在工作過(guò)程中均處于彈性狀態(tài)且使用相同材料，材料本構(gòu)關(guān)系選擇為各向同性彈性材料模型，材料模型參數(shù)為：彈性模量E＝210GPa，密度ρ＝7.8t／m3，泊松比μ＝0.3。軟鋁［13］、紫銅［14］、LY12鋁［14］及45鋼［15］為波形整形器選用材料，采用Johnson－Cook（J－C）本構(gòu)模型。

3 結(jié)果與分析

3.1 子彈對(duì)加速度脈沖的影響

3.1.1 梭形子彈與圓柱形子彈

基于?14.5mm的HPB裝置，將子彈設(shè)計(jì)為如圖3所示的梭形子彈，由削去頂端的圓錐體與圓柱體組合而成，其中D1＝5.0mm，D2＝14.5mm，D3＝2.0mm，L1＝64mm，L2＝40mm，L3＝112mm。

圖3 梭形子彈平面圖Fig.3 Plan of spindle－shaped bullet

該梭形子彈以v＝15m／s的速度撞擊壓桿，得到的加載波（應(yīng)變波）ε波形及加速度a脈沖如圖4所示。

圖4 用梭形子彈和圓柱形子彈得到的加載波波形和加速度脈沖Fig.4 Loading wave shapes and acceleration pusles obtained by spindle－shaped and cylindrical bullets

從圖4可知：采用梭形子彈能夠明顯改善加載波波形上升沿，上升時(shí)間比采用圓柱形子彈時(shí)長(zhǎng)得多，峰值平臺(tái)變窄甚至不出現(xiàn)，形狀也由矩形波變成近似半正弦波，所測(cè)波形基本上消除了由彌散而帶來(lái)的波形振蕩；另外，雖然加速度脈沖的脈寬增加了，但加速度脈沖的峰值降低了很多。

3.1.2 梭形子彈撞擊端面積的影響

改變梭形子彈撞擊端的端面直徑D1，用撞擊端端面直徑不同的梭形子彈以v＝15m／s的速度分別直接撞擊壓桿，得到的加載波波形及加速度脈沖如圖5所示。

圖5 用撞擊端端面直徑不同的梭形子彈得到的加載波波形和加速度脈沖Fig.5 Loading wave shapes and acceleration pusles obtained by spindle－shaped bullets with different－diameter impact ends

從圖5可知：隨著梭形子彈撞擊端面面積的增加，加載波波形被逐漸壓扁，即峰值加大，寬度減小，上升沿變得陡峭；另一方面，加速度脈沖的脈寬減小，峰值提高。

3.1.3 梭形子彈撞擊端錐度的影響

改變梭形子彈端部的長(zhǎng)度L1，獲得不同撞擊端錐度。用撞擊端錐度不同的梭形子彈以v＝15m／s的速度分別直接撞擊壓桿，得到的加載波波形及加速度脈沖如圖6所示。

從圖6可知：隨著梭形子彈撞擊端部錐體錐度的增加，加載波的上升沿略微變陡，應(yīng)變幅值稍微提高；另一方面，對(duì)應(yīng)的加速度峰值也有提升，而加速度脈寬也略微變窄。

圖6 用撞擊端錐度不同的梭形子彈得到的加載波波形和加速度脈沖Fig.6 Loading wave shapes and acceleration pusles obtained by spindle－shaped bullets with different－cone－angle impact ends

3.1.4 梭形子彈尾部的影響

分別改變梭形子彈的尾部直徑D3和尾部長(zhǎng)度L3，用不同尾部直徑和不同尾部長(zhǎng)度的梭形子彈以v＝15m／s的速度分別直接撞擊入射桿，得到的加載波波形及加速度脈沖分別如圖7、8所示。

圖7 用尾部直徑不同的梭形子彈得到的加載波波形和加速度脈沖Fig.7 Loading wave shapes and acceleration pusles obtained by spindle－shaped bullets with different－diameter afterbodies

圖8 用尾部長(zhǎng)度不同的梭形子彈得到的加載波波形和加速度脈沖Fig.8 Loading wave shapes and acceleration pusles obtained by spindle－shaped bullets with different－length afterbodies

從圖7～8可知：改變梭形子彈的尾部直徑或尾部長(zhǎng)度，得到的加速度脈沖峰值及波形基本一致，只是脈寬稍長(zhǎng)。由此說(shuō)明，梭形子彈尾部幾何參數(shù)相對(duì)于撞擊端而言對(duì)加速度脈沖的影響很小。

3.1.5 子彈彈長(zhǎng)的影響

改變梭形子彈中間段的長(zhǎng)度L2，用中間段長(zhǎng)度不同的梭形子彈以v＝15m／s的速度分別直接撞擊壓桿，得到的加載波波形及加速度脈沖如圖9所示。

圖9 用中間段長(zhǎng)度不同的梭形子彈得到的加載波波形和加速度脈沖Fig.9 Loading wave shapes and acceleration pusles obtained by spindle－shaped bullets with different－length middle parts

從圖9可知：用不同彈長(zhǎng)的子彈所得到的波形上升沿走勢(shì)基本一致，脈沖寬度隨著子彈長(zhǎng)度的增大而增加。由此說(shuō)明，子彈長(zhǎng)度對(duì)加速度峰值的影響很小，僅對(duì)加速度脈沖的脈寬有一定影響，隨子彈長(zhǎng)度的增大而略有增加。

3.2 波形整形器對(duì)加速度脈沖的影響

3.2.1 整形器材料的影響

在保持彈形不變（圖3所示）的情況下，在壓桿前端面（撞擊端）加設(shè)一個(gè)金屬墊塊作為整形器。在梭形子彈以相同的速度撞擊下，得到的加載波波形及加速度脈沖如圖10所示。

從圖10可知：加入整形器后，加載波上升沿均明顯被拉寬、峰值減小，其對(duì)應(yīng)的加速度脈沖峰值下降、脈寬增加；整形器材料越軟，加載波形的上升沿越寬，而加速度脈沖峰值越低，脈寬越大。

圖10 不同材料的整形器對(duì)應(yīng)的加載波波形和加速度脈沖Fig.10 Loading wave shapes and acceleration pusles corresponding to the shapers of different materials

3.2.2 整形器尺寸的影響

以圓柱形紫銅作為整形器，分別考慮整形器的直徑d和厚度l對(duì)加載波波形及加速度脈沖的影響，如圖11～12所示。

圖11 不同直徑的整形器對(duì)應(yīng)的加載波波形和加速度脈沖Fig.11 Loading wave shapes and acceleration pusles corresponding to the shapers of different diameters

圖12 不同厚度的整形器對(duì)應(yīng)的加載波波形和加速度脈沖Fig.12 Loading wave shapes and acceleration pusles corresponding to the shapers of different thicknesses

從圖11～12可知：隨著紫銅墊塊直徑的增加，加載波上升時(shí)間減小，而加速度脈沖的峰值增大；隨著紫銅墊塊厚度的增加，加載波上升時(shí)間增大，而加速度脈沖的峰值減小。

4 火工品沖擊實(shí)驗(yàn)加載環(huán)境的數(shù)值模擬

實(shí)驗(yàn)表明：火工品的破壞臨界加速度值隨著加速度脈寬的增大而減小［16］，從而證明研究火工品的抗過(guò)載性能時(shí)，必須考慮加速度脈寬的影響，外界環(huán)境的沖擊加速度峰值和加速度脈寬同時(shí)影響火工品的安全性能。

為此，在上述分析的基礎(chǔ)上，設(shè)計(jì)合理的子彈形狀，選擇適當(dāng)?shù)恼纹鳎烧呦嚓P(guān)因素有效地組合在一起，計(jì)算得到的加載波波形及加速度脈沖如圖13所示。數(shù)值模擬中，子彈形狀為梭形，改變子彈長(zhǎng)度和子彈的撞擊速度；整形器材料為45鋼，直徑均為6mm，改變整形器的厚度。

由圖13可知：通過(guò)控制子彈和整形器的相關(guān)因素，加載波上升沿被拉寬，其形狀近似為半正弦波，加速度脈沖脈寬約為100μs，同時(shí)得到100 000g的加速度峰值。

圖13 不同組合對(duì)應(yīng)的加載波波形和加速度脈沖Fig.13 Loading wave shapes and acceleration pusles corresponding to different combinations

5 結(jié) 論

（1）在一定的撞擊速度范圍內(nèi)，梭形子彈所產(chǎn)生的波形近似于半正弦波，子彈撞擊端幾何參數(shù)對(duì)加載波波形的影響遠(yuǎn)大于子彈尾部。加載波上升時(shí)間與子彈彈長(zhǎng)有關(guān)，并隨彈長(zhǎng)增加而增加。

（2）改變波形整形器的材料、直徑和厚度等任何一個(gè)條件，都會(huì)獲得不同的入射波波形和加速度脈沖，并根據(jù)模擬結(jié)果總結(jié)出其各自的影響規(guī)律。

（3）通過(guò)有效控制子彈和波形整形器等相關(guān)因素，可以獲得所需的加速度脈沖，為檢驗(yàn)火工品在沖擊環(huán)境下或經(jīng)沖擊后性能可靠度的實(shí)驗(yàn)設(shè)計(jì)、測(cè)試等提供依據(jù)。

［1］柳繼昌，席蘭霞.火工品沖擊試驗(yàn)［J］.火工品，1994（4）：32－37.Liu Ji－chang，Xi Lan－xia.Shock test for initiating explosive device［J］.Initiators and Pyrotechnics，1994（4）：32－37.

［2］胡時(shí)勝.霍普金森壓桿技術(shù)［J］.兵器材料科學(xué)與工程，1991（11）：40－47.Hu Shi－sheng.Hopkinson pressure bar technology［J］.Ordnance Material Science and Engineering，1991（11）：40－47.

［3］Bateman V I，Brown F A，Davie N T.Use of a beryllium Hopkinson bar to characterize a piezoresistive accelerometer in shock environments［J］.Journal of the Institute of Environmental Sciences，1996，39（6）：33－39.

［4］Togami T C，Baker W E，F(xiàn)orrestal M J.A split Hopkinson bar technique to evaluate the performance of accelerometers［J］.Journal of Applied Mechanics，1996，63（2）：353－356.

［5］Togami T C，Brown F A，F(xiàn)orrestal M J，et al.Performance evaluation of accelerometers to 200 000g［J］.Journal of Applied Mechanics，1998，65（1）：266－268.

［6］Sill R D.Testing techniques involved with the development of high shock acceleration sensors［R］.ADP002693，1983.

［7］李玉龍，郭偉國(guó)，賈德新，等.高g 值加速度傳感器校準(zhǔn)系統(tǒng)的研究［J］.爆炸與沖擊，1997，17（1）：90－96.Li Yu－long，Guo Wei－guo，Jia De－xin，et al.An equipment for calibrating high shock acceleration sensors［J］.Explosion and Shock Waves，1997，17（1）：90－96.

［8］盛黨紅.高a 值加速度傳感器標(biāo)定方法的理論探討與分析［J］.無(wú)損檢測(cè)，1999，21（12）：544－545.Sheng Dang－h(huán)ong.Theoretical study and analysis on the method for calibrating high a value accelerometers［J］.Nondestructive Testing，1999，21（12）：544－545.

［9］張學(xué)舜，沈瑞琪.火工品動(dòng)態(tài)著靶模擬仿真技術(shù)研究［J］.火工品，2003（4）：1－4.Zhang Xue－shun，Shen Rui－qi.Study on dynamic touch－target analog simulation technique for initiating explosive devices［J］.Initiators and Pyrotechnics，2003（4）：1－4.

［10］鄧強(qiáng)，吳艷霞，沈瑞琪，等.波形整形器在火工品高過(guò)載試驗(yàn)中的應(yīng)用［J］.火工品，2005（5）：12－15.Deng Qiang，Wu Yan－xia，Shen Rui－qi，et al.Application of pulse shaper in over－loading test for initiator［J］.Initiators and Pyrotechnics，2005（5）：12－15.

［11］Forrestal M J，Togami T C，Baker W E，et al.Performance evaluation of accelerometers used for penetration experiments［J］.Experimental Mechanics，2003，43（1）：90－96.

［12］莊茁.ABAQUS非線性有限元分析與應(yīng)用［M］.3版.北京：清華大學(xué)出版社，2005：49－57.

［13］盧劍鋒.沖擊載荷作用下材料和結(jié)構(gòu)力學(xué)行為有限元模擬［D］.北京：清華大學(xué)，2003：51－56.

［14］Johnson G R，Cook W H.A constitutive model and data for metals subjected to large strains，high strain rates and high temperatures［C］∥Proceedings of the Seventh International Symposium on Ballistics.Netherlands：Hague，1983：541－547.

［15］陳剛，陳忠富，陶俊林，等.45鋼動(dòng)態(tài)塑性本構(gòu)參量與驗(yàn)證［J］.爆炸與沖擊，2005，25（5）：451－456.Chen Gang，Chen Zhong－fu，Tao Jun－lin，et al.Investigation and validation on plastic constitutive parameters of 45steel［J］.Explosion and Shock Waves，2005，25（5）：451－456.

［16］鄧瓊，李玉龍，索濤，等.火工品高過(guò)載動(dòng)態(tài)力學(xué)性能測(cè)試方法研究［J］.火工品，2007（1）：28－31.Deng Qiong，Li Yu－long，Suo Tao，et al.Test method on dynamic mechanical behavior of initiating explosive device under high acceleration［J］.Initiators and Pyrotechnics，2007（1）：28－31.