苗飛超, 張向榮, 李東偉, 江濤, 周霖

(1.安徽理工大學(xué) 化學(xué)工程學(xué)院, 安徽 淮南 232001; 2.北京理工大學(xué) 機電學(xué)院, 北京 100081;3.重慶紅宇精密工業(yè)集團有限公司, 重慶 402760)

0 引言

霍普金森桿廣泛應(yīng)用于測試材料在高應(yīng)變率條件下的力學(xué)性能。當(dāng)采用霍普金森桿對試樣進行測試時,試樣應(yīng)處于應(yīng)力平衡和常應(yīng)變率狀態(tài)。然而,在采用傳統(tǒng)的霍普金森桿對炸藥材料進行測試時,試樣兩端很難達到應(yīng)力平衡[1]。這主要是因為炸藥屬于準脆性材料,其失效應(yīng)變很小,可能在材料失效時還未達到應(yīng)力平衡狀態(tài)。此外,常應(yīng)變率狀態(tài)要比動態(tài)應(yīng)力平衡更難實現(xiàn)。采用霍普金森桿對炸藥材料進行測試時,斜坡型的入射波是合適的,因此需要對傳統(tǒng)的霍普金森桿產(chǎn)生的矩形波進行整形。三桿技術(shù)[2]、異型子彈[3]和整形器[4]都可以對入射波波形進行調(diào)整。在入射桿和子彈之間添加整形器是最簡單、實用的入射波整形方法。整形器的塑性變形能夠增加入射波的應(yīng)力上升時間,有利于實現(xiàn)應(yīng)力平衡和常應(yīng)變率。

采用整形器調(diào)整入射波波形的試驗方法已經(jīng)成功應(yīng)用于測試金屬、橡膠和混凝土等多種類型材料[5-7]。整形器不僅應(yīng)用于傳統(tǒng)的霍普金森壓桿(SHPB),還應(yīng)用于霍普金森拉桿[8]裝置,用來測試材料的動態(tài)拉伸性能以及控制入射波卸載段的波形[9]。整形器的材料要根據(jù)待測試樣的力學(xué)性能選擇,銅具有優(yōu)良的延展性,是最常用的整形器材料[10]。在某些特殊情況下,還會采用兩種不同材料整形器組合形成的復(fù)合整形器[11]。整形器增加了入射波的上升時間,使得試樣兩端應(yīng)力平衡變得簡單,然而實現(xiàn)常應(yīng)變率要相對困難。為實現(xiàn)常應(yīng)變率,需要合理設(shè)計試驗條件如整形器材料和尺寸以及子彈的速度和長度。LS-DYNA和Abaqus等商業(yè)軟件也被應(yīng)用于研究整形器對入射波波形的影響[12-14]。然而,采用數(shù)值模擬的方法進行試驗設(shè)計時,需要知道待測試樣的力學(xué)性能,這正是需要試驗測量的。因此,依然需要采用試錯法對待測試樣進行試驗,以確定合適的試驗條件。

Nemat-Nasser等[15]基于一維應(yīng)力波理論首次建立了整形器理論模型,模型假設(shè)整形器是塑性不可壓縮的,在變形過程中處于應(yīng)力平衡狀態(tài),并采用冪律型響應(yīng)函數(shù)描述整形器的應(yīng)力-應(yīng)變關(guān)系,預(yù)測了整形器變形對入射波波形的影響,計算值與試驗值吻合較好。Frew等[16]進一步考慮了整形器的大變形和彈性卸載,使得模型可以預(yù)測更加真實試驗條件下的入射波波形。但是,采用模型預(yù)測需要給定整形器卸載發(fā)生的階段。整形器在不同階段卸載,控制方程也有所不同。然而,在進行計算之前并不知道整形器的卸載時間。因此,Nemat-Nasser模型和Frew模型擴展到雙整形器時,不能預(yù)測整形器卸載段的變形情況[17],也就得不到完整的入射波波形。在入射桿撞擊端粘貼整形器后,相同長度的子彈產(chǎn)生的入射波長度增加,粘貼在桿中部的應(yīng)變片記錄的入射波和反射波信號可能發(fā)生疊加。因此,預(yù)測完整的入射波波形可以避免信號疊加的發(fā)生,對試驗數(shù)據(jù)的有效性至關(guān)重要。

本文基于應(yīng)力波理論分析一維應(yīng)力波在子彈-整形器-入射桿中的傳播過程,充分利用波的疊加原理構(gòu)建整形器變形控制方程,據(jù)此建立能夠預(yù)測預(yù)測完整入射波波形的整形器模型。采用T2純銅作為整形器材料,開展SHPB試驗和整形器試驗,標定T2純銅在一維應(yīng)力狀態(tài)下的本構(gòu)模型。通過對比模型預(yù)測值和試驗值驗證模型的可行性和預(yù)測結(jié)果的準確性。

1 整形器模型

1.1 單整形器

采用Nemat-Nasser等[15]和Frew等[16-17]的假設(shè),認為整形器是塑性不可壓縮的,并且在變形過程中處于應(yīng)力平衡狀態(tài),則有

(1)

(2)

式中:a0和a分別為整形器初始時刻和當(dāng)前時刻的橫截面積;eps(t)為工程應(yīng)變;σin、σst、σps分別為入射桿、子彈和整形器的應(yīng)力(壓正拉負);A為霍普金森桿的橫截面積。整形器的應(yīng)力σps是工程應(yīng)變eps的函數(shù),因此式(1)和式(2)建立了σin和eps的關(guān)系。整形器的工程應(yīng)變率為

(3)

式中:h0為整形器初始厚度;v3和v4分別為子彈/整形器界面和整形器/入射桿界面質(zhì)點3和質(zhì)點4的速度(見圖1),可以通過對應(yīng)力波的傳播過程進行分析得到,從而進一步得到入射波σin。圖1中,1、2、3和4分別為子彈左側(cè)界面、入射桿右側(cè)界面、子彈/整形器界面和整形器/入射桿界面。

圖1 單整形器示意圖

1.2 波傳播過程分析

整形器的厚度遠小于子彈和入射桿的長度,而且銅的聲速較高,波在其中的傳播時間較短,可以認為整形器兩端應(yīng)力平衡,因此應(yīng)力波在整形器中的具體傳播過程可以不用考慮[16]。子彈以速度v0撞擊粘貼在入射桿撞擊端的整形器上,產(chǎn)生幅值相同、方向相反的壓縮波,傳入子彈和入射桿中。傳入子彈中的壓縮波在子彈自由面處反射一個拉伸波,向子彈撞擊端傳播。由于整形器為彈塑性材料,而且在變形過程中其橫截面積逐漸變大,傳入子彈的壓縮波并非矩形波,進而子彈撞擊端反射的拉伸波并不能將子彈中的應(yīng)力卸載至0 MPa。

圖2 應(yīng)力波在子彈-整形器-入射桿中的傳播過程示意圖

相鄰兩次反射的時間間隔τ=2Lst/cb,其中Lst為子彈長度,cb為霍普金森桿的聲速。子彈和入射桿均為線彈性的,則在[nτ,(n+1)τ]時間段,對應(yīng)力和質(zhì)點速度增量應(yīng)用疊加原理,可得

(4)

(5)

(6)

(7)

動態(tài)連續(xù)性條件[18]給出Δσ=±ρbcbΔv,ρb為霍普金森桿的密度,負號表示左行波,Δv為質(zhì)點速度增量。利用動態(tài)連續(xù)性條件以及式(5),可以將式(6)和式(7)改寫為

(8)

(9)

至此已經(jīng)得到了v3和v4的表達式(式(8)和式(9)),將其代入式(3)中,就可以得到整形器應(yīng)變率

(10)

(11)

等號右側(cè)第2項可以按照式(11)格式展開如下:

(12)

然后,將式(12)代入式(11)中,得

(13)

重復(fù)以上兩個步驟,直至式(11)等號右側(cè)第2項消去,最終得到

(14)

最后,將式(14)、式(1)和式(2)代入式(10)中,得到整形器工程應(yīng)變的表達式:

(15)

式中:整形器應(yīng)力σps是整形器工程應(yīng)變eps函數(shù),則求解此常微分方程就可以得到eps。進而可以根據(jù)式(1)和式(2)得到入射波波形σin(t)。

在以上入射波波形表達式的推導(dǎo)過程中只用到了應(yīng)力波基本原理。這些基本原理在一維應(yīng)力條件下普遍成立,不會因為整形器處于加載階段還是卸載階段而有所改變。因此,由式(15)、式(1)和式(2)可以計算得到包含加載段和卸載段的完整入射波,而不必像Nemat-Nasser等[15]和Frew等[16]的模型一樣對卸載段單獨處理,而且在不同階段卸載對應(yīng)不同的公式。在雙整形器和多整形器情況下,這種優(yōu)點體現(xiàn)得更加明顯。

1.3 雙整形器和多整形器

對脆性材料而言,合理地調(diào)整整形器尺寸和其他試驗條件,可以使試樣處于應(yīng)力平衡和常應(yīng)變率狀態(tài)。但是,當(dāng)待測試樣為金屬材料時需要兩個整形器[19]。兩個整形器一般采用兩種不同的材料,其中一種材料強度較高如鋼,另一種材料強度相對較低如銅。兩個整形器用一個可以近似為剛體的高強度、高硬度的金屬圓片隔開。兩個整形器和金屬片的相對位置如圖3所示。圖3中,3、4、5和6分別為子彈/整形器1界面、整形器1/金屬圓片界面、金屬圓片/整形器2界面和整形器2/入射桿界面。

圖3 雙整形器示意圖

兩個整形器的工程應(yīng)變率表達式為

(16)

(17)

式中:上標1和2表示整形器編號;v5(t)和v6(t)分別為金屬圓片/整形器2界面和整形器2/入射桿界面的質(zhì)點速度。隔開兩個整形器的金屬圓片硬度高、強度高,可以將其作為剛體處理,則金屬圓片兩端質(zhì)點速度相同,即v4(t)=v5(t)。將此式代入式(16)和式(17)中,得

(18)

式中等號右側(cè)的質(zhì)點速度項可以按照1.1節(jié)的方法消去,消去此項后的表達式為

(19)

(20)

以上雙整形器的推導(dǎo)過程很容易擴展到多個整形器的情形。類似于式(19)和式(20),多整形器的控制方程為

(21)

式中:M為整形器個數(shù)。為方便求解式(21),開發(fā)了多整形器(MPS)代碼。MPS采用輸入文件的形式,將桿件信息和整形器的材料特性輸入到文件中。

2 SHPB與整形器試驗

為采用整形器模型對入射波波形進行預(yù)測,首先需要知道整形器的應(yīng)力-應(yīng)變關(guān)系。在一維應(yīng)力狀態(tài)下,可以采用簡單的冪函數(shù)表達式

(22)

描述整形器的應(yīng)力-應(yīng)變關(guān)系。式中:εen為工程應(yīng)變;σ0、n和m均為待標定參數(shù)。

為得到T2純銅的真應(yīng)力-工程應(yīng)變關(guān)系,對T2純銅樣品進行標準SHPB試驗,并對不同尺寸的T2純銅進行整形器試驗。試驗在重慶紅宇精密工業(yè)集團有限公司動態(tài)加載實驗室完成。試驗采用鋼桿對T2純銅進行測試,使得T2純銅應(yīng)變盡可能大,進而擴展整形器的應(yīng)力-應(yīng)變關(guān)系的適用范圍,提高整形器的預(yù)測能力。SHPB的桿件材料參數(shù)如表1所示。SHPB試驗采用的樣品尺寸分別為φ5.0 mm×2.5 mm和φ5.0 mm×5.0 mm,相應(yīng)的應(yīng)變率分別為5 000 s-1和3 500 s-1,試驗結(jié)果如圖4所示。

表1 桿件材料參數(shù)

圖4 T2純銅的真應(yīng)力-工程應(yīng)變關(guān)系

整形器試驗給出每次試驗整形器的峰值應(yīng)力和峰值應(yīng)變。由式(1)和式(2)可知,整形器中的峰值應(yīng)力可以由應(yīng)變片記錄的入射桿的峰值應(yīng)力計算。

整形器的峰值應(yīng)變則不能直接測量。但是整形器由峰值應(yīng)力減小到0 MPa的過程中,整形器彈性卸載,相應(yīng)的變形可以忽略。計算時取T2純銅的彈性模量為100 GPa,在試驗的應(yīng)力范圍內(nèi)忽略彈性卸載變形導(dǎo)致的誤差小于2%。因此,在試驗后回收整形器測量其尺寸就可以計算峰值應(yīng)變。每發(fā)試驗可以得到1個試驗點,針對T2純銅進行了多發(fā)試驗,試驗結(jié)果如圖4所示。

對圖4中T2純銅的SHPB和整形器試驗結(jié)果進行標定,標定的參數(shù)為σ0=435.6 MPa,n=0.113 8,m=3.447。參數(shù)標定完畢后,可以采用模型對子彈撞擊整形器所產(chǎn)生的入射波進行預(yù)測。圖5 計算了兩種不同工況下的入射波波形并與相應(yīng)試驗波形進行了對比。試驗采用鋼桿,子彈長度為200 mm。兩種工況下子彈速度和整形器尺寸分別為:工況1,子彈速度8.03 m/s,整形器尺寸φ3.0 mm×1.6 mm;工況2,子彈速度9.48 m/s,整形器尺寸φ5.0 mm×2.5 mm。由圖5可以看出,兩種工況條件下,模型預(yù)測的入射波峰值點和卸載點等特征點均與測量值的偏差較小,入射波波形與試驗測量波形整體吻合較好,驗證了整形器模型的可行性和整形器應(yīng)力-應(yīng)變關(guān)系參數(shù)標定的準確性。

圖5 兩種不同工況下MPS計算值與試驗值對比

為了分析子彈撞擊過程中整形器的變形特性,圖6給出了工況1條件下200 mm鋼桿撞擊整形器時,整形器兩個端面的質(zhì)點速度和工程應(yīng)變隨時間的變化情況。由圖6可見:子彈撞擊整形器后,v3>v4,故整形器開始發(fā)生壓縮變形,應(yīng)變增加的同時應(yīng)力也在增加,波在子彈中傳播一個來回的時間τ=80 μs;當(dāng)t≈3τ時,v3>v4整形器應(yīng)變開始減小,應(yīng)力開始逐漸卸載到0 MPa;從應(yīng)變最高點到整形器完全卸載,應(yīng)變變化為1.28%,表明以上采用回收整形器的應(yīng)變代替峰值應(yīng)變是可行的。

圖6 整形器變形程度隨時間的變化

此外,圖6也體現(xiàn)了建立整形器模型的必要性。若采用常規(guī)的SHPB,即沒有添加整形器,則由應(yīng)力波理論可知入射波的波長應(yīng)為τ。添加整形器后,入射波波長明顯增加。由τ=2Lst/cb可知,若子彈長度Lst較大,則添加整形器后,粘貼在桿中部的應(yīng)變片記錄的入射波和反射波信號可能發(fā)生疊加。因此,預(yù)測完整的入射波波形可以避免信號疊加的發(fā)生,為整形器設(shè)計提供指導(dǎo)。

為進一步驗證整形器模型的適用性,圖7采用MPS對Frew等[17]的銅/鋼雙整形器試驗進行了預(yù)測,詳細的試驗條件參見文獻[17]。圖7表明,MPS和Frew模型的計算值均與試驗值吻合較好,但是Frew模型只能預(yù)測加載段,而MPS可以預(yù)測包含卸載段在內(nèi)的完整的入射波波形。由于Frew模型不能預(yù)測卸載段,文獻[17]中并未給出整形器材料的彈性模量。圖7計算時銅和鋼在卸載時的彈性模量分別為100 GPa和200 GPa。銅/鋼雙整形器在試驗過程中的變形特性與單個整形器相比更加復(fù)雜。

圖7 雙整形器時MPS計算值與試驗值[17]對比

圖8給出了在圖7(a)工況下兩個整形器端面的質(zhì)點速度隨時間的變化情況。質(zhì)點v4=v5速度(詳見1.2節(jié)),故圖8中v4和v5位于同一條曲線。結(jié)合式(16)和式(17)可知,整形器的應(yīng)變與兩端質(zhì)點速度所圍成的面積呈正比關(guān)系。銅的強度小于鋼,因此子彈撞擊整形器后,銅整形器開始產(chǎn)生塑性變形,而鋼整形器兩端的質(zhì)點速度相差極小,處于彈性變形階段。隨著應(yīng)力的不斷增加,在t=83 μs時鋼整形器開始發(fā)生塑性變形,此時銅整形器的工程應(yīng)變已經(jīng)達到了89%。此后,應(yīng)力增加的幅度開始減小(見圖7(a)),鋼整形器的應(yīng)變開始大幅增加,而銅整形器的應(yīng)變盡管也在增加,但是增加的幅值不超過2%。波在子彈中傳播一個來回后,即t>τ=123.2 μs時整形器仍處于加載階段,直至t=198.8 μs兩個整形器均開始卸載。當(dāng)應(yīng)力完全卸載后,銅和鋼整形器的工程應(yīng)變分別為89.5%和4.8%。通過以上分析可知,采用兩種不同材料作為整形器時,入射波波形的初始階段主要由強度低的整形器控制,當(dāng)強度高的整形器開始發(fā)生塑性變形后,則由其控制入射波波形。

圖8 銅/鋼雙整形器端面質(zhì)點速度隨時間的變化

3 結(jié)論

1) 通過分析一維應(yīng)力波在子彈-整形器-入射桿中的傳播過程,發(fā)現(xiàn)充分利用波的疊加原理可以得到整形器加載-卸載的統(tǒng)一控制方程,據(jù)此建立了整形器模型,實現(xiàn)了整形器彈性加載-塑性加載-彈性卸載的自動轉(zhuǎn)換,使模型具備預(yù)測包含卸載段在內(nèi)的完整的入射波波形的能力。

3) 采用整形器模型預(yù)測了單整形器和雙整形器試驗的入射波波形,結(jié)果表明模型預(yù)測的入射波峰值點和卸載點等特征點均與試驗值的偏差較小,入射波波形與試驗測量波形整體吻合較好,驗證了整形器模型的可行性和整形器應(yīng)力-應(yīng)變關(guān)系參數(shù)標定的準確性。