孫鐘阜,汪文強(qiáng),鄭 健,陳 雄,許進(jìn)升
(1.海軍駐上海水聲導(dǎo)航代表室,上海201108;2.武漢濱湖電子有限公司,武漢430000;3.南京理工大學(xué) 機(jī)械工程學(xué)院,南京210094)
武器工程應(yīng)用中,為改善含能材料的能量和力學(xué)特性,往往通過在其基體中添加含能晶粒組分(如高氯酸銨AP、金屬燃燒劑)的方式進(jìn)行改性。其中改性雙基推進(jìn)劑、塑性粘結(jié)炸藥等含能材料多以混合體系的形式復(fù)合而成,含能晶體必然是混合體系中的核心組元,因而對材料的力學(xué)性能、燃燒特性有顯著的影響[1]。含能材料在服役過程中,在沖擊載荷帶來的高應(yīng)變率作用下,含能晶體之間因相互擠壓極易發(fā)生變形甚至破裂,顆粒聚集體的微觀破壞行為對含能材料的力學(xué)性能有顯著的影響。
沖擊載荷作用下含能材料的強(qiáng)度與斷裂行為是武器物理研究和工程應(yīng)用中的基礎(chǔ)問題。近年來,針對固體推進(jìn)劑和含能炸藥的動態(tài)力學(xué)性能,國內(nèi)外相關(guān)領(lǐng)域?qū)W者開展了一系列的研究工作。趙玉剛等[2]結(jié)合平臺巴西實(shí)驗和SHPB實(shí)驗技術(shù)研究了3種PBX炸藥的動態(tài)拉伸力學(xué)性能。盧芳云[3]對2種不同的炸藥進(jìn)行了動態(tài)壓縮實(shí)驗,并利用電鏡掃描技術(shù)分析了試件破壞后的微觀結(jié)構(gòu),同時建立了相應(yīng)的本構(gòu)模型以描述材料在高應(yīng)變下的力學(xué)特性。李俊玲[4]利用SHPB實(shí)驗技術(shù)建立了PBX炸藥含細(xì)觀損傷的本構(gòu)關(guān)系。文獻(xiàn)[5]研究了HTPB推進(jìn)劑在點(diǎn)火沖壓載荷下的動態(tài)力學(xué)特性。文獻(xiàn)[6]開展了HTPB推進(jìn)劑在高應(yīng)變率條件下的非線性累積損傷研究,并建立了推進(jìn)劑在應(yīng)變率為103~104s-1范圍內(nèi)的損傷本構(gòu)模型。
通常情況下復(fù)合材料的拉伸破壞強(qiáng)度明顯小于其壓縮強(qiáng)度。往往固體推進(jìn)劑在服役過程中最先在拉伸應(yīng)力作用下失效,因而研究固體推進(jìn)劑的動態(tài)拉伸力學(xué)特性對評估推進(jìn)劑的使用安全性具有一定的研究意義。
本文利用SHPB實(shí)驗技術(shù)對改性雙基推進(jìn)劑進(jìn)行了動態(tài)沖擊實(shí)驗,結(jié)合掃描電鏡技術(shù),基于微觀結(jié)構(gòu)分析討論了CMDB推進(jìn)劑內(nèi)部微觀結(jié)構(gòu)的破壞機(jī)理對其力學(xué)性能的影響。
實(shí)驗中所研究的CMDB推進(jìn)劑基礎(chǔ)配方為:硝化棉(NC)20.5%,硝化甘油(NG)21%,黑索今(RDX)32.5%,高氯酸銨(AP)13.7%,鋁粉(Al)8.4%,催化劑2.2%,炭黑(CB)0.4%,及其他添加成分。實(shí)驗樣品均采用螺壓工藝制備。試件結(jié)構(gòu)圖及相關(guān)尺寸見圖1。
圖1 試件結(jié)構(gòu)示意圖
動態(tài)沖擊實(shí)驗在某研究室SHPB實(shí)驗系統(tǒng)上進(jìn)行。圖2為SHPB實(shí)驗裝置圖。沖擊條件下,由于固體推進(jìn)劑的破壞應(yīng)變較小,實(shí)驗中所采用的試件結(jié)構(gòu)不再是一維單軸壓縮試件,其拉伸破壞應(yīng)變不能直接通過壓桿上的測量信號計算得到。因此,實(shí)驗中采用在試件表面直接粘貼應(yīng)變片的方式進(jìn)行測量,如圖2所示。
圖2 SHPB實(shí)驗裝置圖及應(yīng)變片粘貼方式實(shí)物圖
SHPB實(shí)驗要求試件在整個沖擊加載過程中滿足動態(tài)應(yīng)力/應(yīng)變平衡條件,即εi(t)+εr(t)=εt(t),式中:εi(t)、εr(t)、εt(t)分別為入射波、反射波、透射波在壓桿內(nèi)傳播時所引起的應(yīng)變。依據(jù)這一基本假設(shè)對實(shí)驗數(shù)據(jù)的合理有效性進(jìn)行校驗。
為使試件在承載過程中滿足應(yīng)力/應(yīng)變均勻假設(shè),采用脈沖整形技術(shù)[7]獲取良好的加載脈沖。圖3為典型實(shí)驗脈沖波形圖。顯然,從圖3中可以看出,經(jīng)過整形后獲取的加載脈沖變寬,其上升沿變緩,有利于材料在脈沖的上升過程中盡快達(dá)到應(yīng)力平衡狀態(tài)[8]。為校驗數(shù)據(jù)的合理性,對每個實(shí)驗數(shù)據(jù)進(jìn)行波形分離檢測,剔除不合理的數(shù)據(jù)。圖4給出了實(shí)驗數(shù)據(jù)波形分離圖,從圖4中可以看出,計算透射信號與實(shí)驗透射吻合良好;結(jié)合圖4還可以看出,反射信號出現(xiàn)典型的很長一段平臺區(qū),這意味著試件在整個沖擊加載時域內(nèi)基本實(shí)現(xiàn)了恒應(yīng)變率加載要求,并滿足應(yīng)力平衡狀態(tài)。
圖3 典型實(shí)驗脈沖波形圖
圖4 實(shí)驗數(shù)據(jù)波形分離校驗圖
圖5 給出了相同沖擊速度條件下試件中心位置承受拉應(yīng)力σt的重復(fù)性實(shí)驗數(shù)據(jù),可以看出,同一實(shí)驗條件下所獲取的實(shí)驗數(shù)據(jù)在誤差允許的范圍內(nèi)具有良好的重復(fù)性,尤其在上升沿時段內(nèi)幾乎重合,這說明通過采用合理的整形技術(shù)和嚴(yán)格控制沖擊桿的速度,能夠保證實(shí)驗數(shù)據(jù)的可靠性。
圖5 同組實(shí)驗數(shù)據(jù)重復(fù)性校驗圖
圖6 為子彈沖擊速度達(dá)到10.5m/s時試件表面應(yīng)變片測得的結(jié)果,可以看出,試件中心位置的拉伸應(yīng)變在A點(diǎn)(31.6μs)達(dá)到最大值,對應(yīng)的測量應(yīng)變?yōu)?.015 3,在119μs出現(xiàn)第二峰值點(diǎn)B,對應(yīng)的測量應(yīng)變?yōu)?.015 1,此后測量信號急劇下降。這說明在A點(diǎn)之前,推進(jìn)劑的動態(tài)拉伸應(yīng)變呈線性遞增趨勢;AB段測量結(jié)果呈現(xiàn)先下降后急劇上升的趨勢,由于這一時域內(nèi)應(yīng)力波在試件中來回傳播,試件中心部位在壓縮波和拉伸波的作用下,應(yīng)變出現(xiàn)波動回復(fù)現(xiàn)象;點(diǎn)B處出現(xiàn)第二峰值后,測量信號急劇下降并趨近于0,顯然,試件中心位置在B點(diǎn)達(dá)到最大應(yīng)變后已經(jīng)發(fā)生破壞或嚴(yán)重?fù)p傷,導(dǎo)致應(yīng)變片破壞失效。值得注意的是,應(yīng)變片從監(jiān)測到應(yīng)力波的作用開始到出現(xiàn)極值的時間為31.6μs,顯然滯后于試件開始滿足恒應(yīng)變率狀態(tài)所需時間(見圖4),這從一定程度上說明試件是在滿足應(yīng)力平衡狀態(tài)后才開始出現(xiàn)失效現(xiàn)象。需要說明的是,試件表面上所粘貼的應(yīng)變片在不超過其有效量程的前提下,實(shí)時監(jiān)測推進(jìn)劑材料的屈服應(yīng)變。AB段測量應(yīng)變急劇下降后并出現(xiàn)第二峰值,這一時域內(nèi)應(yīng)力波在試件內(nèi)部反復(fù)傳播,試件內(nèi)部應(yīng)力和應(yīng)變分布十分復(fù)雜,本文不作深入分析。
圖6 動態(tài)拉伸應(yīng)變測量結(jié)果
圖7 為改性雙基推進(jìn)劑試件承受沖擊載荷后的典型破壞形貌??梢钥闯?,裂紋直接從試件中心位置沿直徑方向向兩加載端擴(kuò)展,這說明推進(jìn)劑材料是在拉應(yīng)力的作用下破壞失效的。
圖7 推進(jìn)劑典型破壞形貌
由于實(shí)驗中所采用的試件結(jié)構(gòu)形狀不是一維單軸圓柱試件,不能通過壓桿上測量到的沖擊載荷直接獲取試件的拉伸應(yīng)力,因而需要對這一幾何構(gòu)形的試件進(jìn)行受力分析,計算得到試件中心位置拉伸應(yīng)力。依據(jù)彈塑性力學(xué)理論推導(dǎo)出試件中心位置的拉應(yīng)力[9]:
式中:p為施加在試件承載端面的載荷,D、δ分別為直徑和厚度;k為尺寸結(jié)構(gòu)因子,其值大小由試件承載端面所對應(yīng)的圓心角θ決定。此次研究中,試件加載平臺對應(yīng)的圓心角為20°。
依據(jù)一維應(yīng)力波理論,試件承載端面的載荷[10]為
式中:E、A分別為壓桿的彈性模量和橫截面積;εi(t)、εr(t)、εt(t)分別為入射、反射、透射測量信號。
依據(jù)式(1)~式(3),圖8給出了推進(jìn)劑動態(tài)拉伸載荷-時間圖。從圖中可以看出,其拉伸載荷-時間曲線可分為3個階段:Ⅰ初始線性承載階段;Ⅱ動態(tài)載荷平臺區(qū),即載荷-時間曲線偏離初始線性關(guān)系,開始出現(xiàn)屈服現(xiàn)象,此后出現(xiàn)在一段較長時間范圍內(nèi)的近似載荷平臺區(qū);Ⅲ破壞失效階段,隨應(yīng)力波的連續(xù)作用,直至試件出現(xiàn)宏觀的裂紋,導(dǎo)致應(yīng)力驟降。
圖8 典型動態(tài)拉伸載荷圖
從圖8中可觀察到,材料開始出現(xiàn)第一峰值對應(yīng)的時間為32μs左右,應(yīng)變片監(jiān)測到試件中心位置應(yīng)變開始突變的時間為31.6μs(圖6),這是由于SHPB系統(tǒng)中的各應(yīng)變片粘貼位置的空間位置差異,導(dǎo)致測量信號不能滿足嚴(yán)格的同步性,存在測量時差是無法避免的。實(shí)驗所測應(yīng)變峰值對應(yīng)的時刻與應(yīng)力峰值對應(yīng)的時刻基本吻合。說明圖8中第一峰值點(diǎn)A可以判定為推進(jìn)劑的初始屈服點(diǎn)。
此外,圖8中還觀察到一段近似動態(tài)載荷平臺區(qū)(Ⅱ區(qū)),因試件自身結(jié)構(gòu)在承載時為非一維方向承載,初始屈服后并未立即失效,在此期間內(nèi),應(yīng)力波在推進(jìn)劑內(nèi)部來回反復(fù)傳播,不斷對其微觀結(jié)構(gòu)構(gòu)成損傷效應(yīng),當(dāng)損傷達(dá)到一定程度時,材料才完全失效。由圖8觀察到推進(jìn)劑材料發(fā)生失效行為對應(yīng)的時間約為116μs,這與應(yīng)變片測定的第二峰值點(diǎn)B對應(yīng)的時間相差并不大。這表明試件表面上的應(yīng)變片能夠?qū)崟r監(jiān)測到推進(jìn)劑材料的屈服和失效時刻。需要說明的是,由于應(yīng)力波的來回反復(fù)傳播效應(yīng)的影響,圖6中所監(jiān)測到的實(shí)時應(yīng)變信號中AB段時域內(nèi),試件內(nèi)部動態(tài)應(yīng)變過程十分復(fù)雜,因此應(yīng)變片在這一時域內(nèi)所測信號只能與壓桿上應(yīng)變片所測信號進(jìn)行定性對比,對確定材料的初始屈服和失效點(diǎn)起到一定的參考作用。
對實(shí)驗數(shù)據(jù)可靠性進(jìn)行校驗后,同組實(shí)驗數(shù)據(jù)中挑選有效數(shù)據(jù),并結(jié)合試件上應(yīng)變片所測信號確定不同加載率下材料的屈服點(diǎn),依據(jù)式(1)~式(3)計算得到CMDB推進(jìn)劑的動態(tài)拉伸屈服強(qiáng)度值與加載率對應(yīng)的關(guān)系,如圖9所示。
圖9 動態(tài)拉伸屈服強(qiáng)度與加載率對應(yīng)關(guān)系
從圖9中可以觀察到CMDB推進(jìn)劑的動態(tài)拉伸強(qiáng)度隨加載率的增加而增大,表現(xiàn)出明顯的加載率敏感性;此外,從圖9中還可以看出,隨加載的逐漸增加,相同加載率下獲得的屈服強(qiáng)度值散差明顯增大。由于實(shí)驗時加載桿速度由氣罐控制,撞擊速度增大,難免存在初始誤差,這一現(xiàn)象也有可能與推進(jìn)劑內(nèi)部含能晶體在高加載率下的變形和破壞方式有關(guān)。
高應(yīng)變率沖擊載荷下,含能晶體顆粒的破壞行為十分復(fù)雜。Delhaye V[11]研究發(fā)現(xiàn)復(fù)合材料的動態(tài)力學(xué)特性與其內(nèi)部微觀結(jié)構(gòu)的破壞機(jī)理有直接的聯(lián)系。為從微觀角度分析推進(jìn)劑在高應(yīng)變率條件下所表現(xiàn)出的率敏感性,采用SEM掃描技術(shù)對破壞后的材料內(nèi)部的含能晶粒組分進(jìn)行微觀表征。
圖10給出了CMDB推進(jìn)劑在4種不同加載率下的微觀失效形貌電鏡圖。
圖10 不同加載率下試件斷面微觀電鏡圖
圖10 (a)為加載率達(dá)到400GPa/s時的斷面電鏡圖,從圖中可以看到少數(shù)AP顆粒出現(xiàn)損傷和破壞現(xiàn)象,大部分AP顆粒均已脫離基體,幾乎裸露在基體表面,這與準(zhǔn)靜態(tài)單軸拉伸條件下有明顯的區(qū)別,文獻(xiàn)[12]研究發(fā)現(xiàn),準(zhǔn)靜態(tài)拉伸條件下推進(jìn)劑發(fā)生斷裂行為時,AP顆粒依然鑲嵌在基體中。這說明動態(tài)條件下,試件主要以基體材料(硝化棉和硝化甘油)承載,其分子鏈承受高應(yīng)力,因熱激活而斷裂,導(dǎo)致周圍分子鏈上的應(yīng)力重新分布,應(yīng)力集中區(qū)域的分子鏈?zhǔn)紫劝l(fā)生斷裂行為,并積累成微孔洞。實(shí)驗中推進(jìn)劑在滿足應(yīng)力平衡條件下,微孔洞數(shù)量積累需要一定的時間,材料屈服后并不會立即失效。
從圖8可以看出,材料在屈服點(diǎn)后出現(xiàn)一近似應(yīng)力平臺區(qū),隨時間變化,應(yīng)力呈逐漸減小趨勢,這一時域內(nèi)是推進(jìn)劑內(nèi)部微孔洞不斷積累的過程區(qū)。當(dāng)微孔洞積累到一定數(shù)量后,就會在應(yīng)力集中區(qū)擴(kuò)展為如圖9中的宏觀裂紋,導(dǎo)致材料最終失效。
圖10(b)為加載率達(dá)到450GPa/s時的失效形面電鏡圖,顯然隨加載率增大,該視場內(nèi)可觀察到大部分大粒徑AP顆粒從基體中被拔出而留下清晰的凹坑,且部分AP顆粒已出現(xiàn)明顯的破損現(xiàn)象,說明這一加載率下,應(yīng)力波直接對推進(jìn)劑內(nèi)部微觀結(jié)構(gòu)構(gòu)成損傷效應(yīng),此時推進(jìn)劑材料不僅僅是限于基體材料承載,其內(nèi)部大粒徑尺寸的微觀顆粒結(jié)構(gòu)在拉伸應(yīng)力作用下發(fā)生變形和損傷破壞,因此無疑會消耗更多的能量。
圖10(c)為加載率達(dá)到490GPa/s時的失效斷面電鏡圖,從圖中可以明顯觀察到,鑲嵌在基體內(nèi)部相當(dāng)數(shù)量的AP顆粒發(fā)生穿晶斷裂行為,這說明隨加載率的增強(qiáng),應(yīng)力波引起推進(jìn)劑材料內(nèi)部的微觀顆粒結(jié)構(gòu)的損傷效應(yīng)明顯加劇。
圖10(d)為加載率達(dá)到660GPa/s時的斷面電鏡圖,從圖中可以直觀看到,大粒徑AP顆粒結(jié)構(gòu)破壞程度進(jìn)一步惡化,斷面內(nèi)幾乎觀察不到完整的顆粒結(jié)構(gòu)。這說明應(yīng)力波在推進(jìn)劑內(nèi)部的微觀顆粒結(jié)構(gòu)和基體界面之間來回傳播,導(dǎo)致顆粒結(jié)構(gòu)在拉伸和壓縮應(yīng)力波的作用下發(fā)生脆性破壞現(xiàn)象。
顯然,根據(jù)所觀察到的不同加載率下推進(jìn)劑內(nèi)部的微觀結(jié)構(gòu)進(jìn)行分析,推進(jìn)劑在高應(yīng)變條件下表現(xiàn)出的明顯的加載率敏感性與應(yīng)力波對其內(nèi)部微觀結(jié)構(gòu)的破壞機(jī)理有直接的聯(lián)系。推進(jìn)劑材料在加載率持續(xù)增強(qiáng)時,孔洞和微裂紋的孕育、成長所需時間減短,進(jìn)而直接導(dǎo)致推進(jìn)劑的整個破壞過程弛豫時間縮短;基體內(nèi)部顆粒微觀結(jié)構(gòu)在應(yīng)變率持續(xù)增大的情況下,變形時間縮短,而裂紋擴(kuò)展區(qū)是典型的應(yīng)力集中部位,因而大粒徑微觀顆粒在來不及充分變形的情況下直接發(fā)生斷裂甚至破碎行為,而不同形式的斷裂行為無疑會消耗額外的斷裂能,所以推進(jìn)劑的宏觀動態(tài)屈服失效強(qiáng)度相應(yīng)增強(qiáng)。因此,從能量損耗角度分析,推進(jìn)劑在高應(yīng)變率下拉伸強(qiáng)度存在率敏感性,與AP顆粒破壞數(shù)量和程度有直接的聯(lián)系。
CMDB推進(jìn)劑在高應(yīng)變率動態(tài)拉伸條件下,其動態(tài)載荷力學(xué)特征曲線由3個階段表征:線彈性上升階段、損傷屈服階段、破壞失效階段。
動態(tài)沖擊條件下,隨著加載率增大,CMDB推進(jìn)劑材料的屈服強(qiáng)度相應(yīng)增大,表現(xiàn)出明顯的率敏感特性。
微觀結(jié)構(gòu)分析表明,推進(jìn)劑在高應(yīng)變率下拉伸強(qiáng)度存在率敏感性,與其基體內(nèi)部微觀顆粒結(jié)構(gòu)破壞數(shù)量和程度直接相關(guān)。
[1]李吉禎,樊學(xué)忠,鐘 雷,等.NC/NG/AP/Al復(fù)合改性雙基推進(jìn)劑力學(xué)性能研究[J].含能材料,2007,15(4):345-348.LI Ji-zhen,F(xiàn)AN Xue-zhong,ZHONG Lei,et al.Mechanical properties of NC/NG/AP/AL composite modified double base propellant[J].Chinese Journal of Energetic Materials,2007,15(4):345-348.(in Chinese)
[2]趙玉剛,傅 華,李俊玲,等.三種PBX炸藥的動態(tài)拉伸力學(xué)性能[J].含能材料,2011,19(2):194-199.ZHAO Yu-gang,F(xiàn)U Hua,LI Jun-ling et al.Dynamic tensile mechanical properties of three types of PBX[J].Chinese Journal of Energetic Materials,2011,19 (2):194 - 199.(in Chinese)
[3]盧芳云,吳會民,王曉燕,等.兩種炸藥材料本構(gòu)行為的應(yīng)變率效應(yīng)分析[J].含能材料,2004,12(1):280-285.LU Fang-yun,WU Hui-min,WANG Xiao-yan,et al.Analysis of strain-rate effects on constitutive behaviors of two explosives[J].Chinese Journal of Energetic Materials,2004,12(1):280-285.
[4]李俊玲.PBX炸藥裝藥的力學(xué)性能及損傷破壞研究[D].長沙:國防科學(xué)技術(shù)大學(xué),2012.LI Jun-ling.Study on PBX’s mechanical behavior and damage feature[D].Changsha:National University of Defense Technology,2012.(in Chinese)
[5]CHYUAN S W.Dynamic analysis of solid propellant grains subjected to ignition pressurization loading[J].Journal of Sound and Vibration,2003,268(3):465-483.
[6]HO S Y.High strain-rate constitutive models for solid rocket propellants[J].Journal of Propulsion and Power,2002,18(5):1 106-1 111.
[7]FREW D J,F(xiàn)ORRESTAL M J,CHEN W.Pulse shaping techniques for testing brittle materials with a split Hopkinson pressure bar[J].Experimental Mechanics,2002,42(1):93-106.
[8]SONG B,CHEN W.Loading and unloading split Hopkinson pressure bar pulse-shaping techniques for dynamic hysteretic loops[J].Experimental Mechanics,2004,44(6):622-627.
[9]WANG Q Z,JIA X M,KOU S Q,et al.The flattened Brazilian disc specimen used for testing elastic modulus,tensile strength and fracture toughness of brittle rocks:analytical and numerical results[J].International Journal of Rock Mechanics and Mining Sciences,2004,41(2):245-253.
[10]CHEN X,WU S,ZHOU J.Quantification of dynamic tensile behavior of cement-based materials[J].Construction and Building Materials,2014,51:15-23.
[11]DELHAYE V,CLAUSEN A H,MOUSSY F,et al.Influence of stress state and strain rate on the behavior of a rubber-particle reinforced polypropylene[J].International Journal of Impact Engineering,2011,38(4):208-218.
[12]IDE K M,HO S Y,WILLIAMS D R G.Fracture behavior of accelerated aged solid rocket propellants[J].Journal of Materials Science,1999,34(17):4 209-4 218.