劉 雙,張錫銘,張 京,鄭夢澤,劉文皓,羅運軍
(1.北京理工大學(xué) 材料學(xué)院,北京 100081;2.高能量密度材料教育部重點實驗室,北京 100081)
固體推進劑在運輸、裝載及長期貯存過程中,由于自重及受到的外部載荷的作用,會產(chǎn)生軸向和徑向的蠕變形變,導(dǎo)致固體推進劑藥型和尺寸發(fā)生變化,使推進劑的通氣面積變小,工作時壓力峰劇增,嚴(yán)重影響其燃燒與安全系數(shù)[1],影響固體推進劑的使用壽命,因此,研究固體推進劑蠕變性能具有重要意義。目前,國內(nèi)外對多種類型的固體推進劑蠕變性能進行了研究。胡義文等[2]研究了加載應(yīng)力、溫度、黏合劑組成對熱固性PBT基固體推進劑蠕變性能的影響,并通過Burgers模型研究其蠕變機制;Zhang J B等[3]研究了常溫和高溫下應(yīng)力對雙基推進劑蠕變性能的影響,并比較了5種模型對推進劑蠕變曲線的擬合程度;李東等[4]研究了應(yīng)力、溫度對雙基推進劑蠕變性能的影響,并構(gòu)建了在20℃以及5MPa加載應(yīng)力下的蠕變?nèi)崃恐髑€。Feng S X等[5]采用冪律本構(gòu)模型、線性黏彈性模型和Schapery非線性本構(gòu)模型研究了HTPB推進劑的蠕變過程。
GAP基含能熱塑性彈性體(GAP-ETPE)推進劑具有能量高、力學(xué)性能優(yōu)異、工藝簡單等諸多優(yōu)點,同時作為一種綠色推進劑,GAP-ETPE基固體推進劑實現(xiàn)了重回收、重循環(huán)和重利用[6-7],解決了傳統(tǒng)熱固性固體推進劑的廢棄問題[8],具有良好的應(yīng)用前景,是復(fù)合固體推進劑領(lǐng)域研究的前沿與熱點。然而針對其蠕變性能的研究較少,尤其是未見針對其高溫蠕變性能以及長期蠕變行為的相關(guān)報道。
GAP-ETPE基固體推進劑是一種高填充的熱塑性聚合物復(fù)合材料,其蠕變性能同樣對溫度與加載應(yīng)力具有很強的依賴性。由于熱塑性材料宏觀力學(xué)性能主要依賴于微觀結(jié)構(gòu)間的物理作用形成的物理交聯(lián)點[9],其本身具有在高溫下物理交聯(lián)網(wǎng)絡(luò)被應(yīng)力破壞而獲得可加工性能的特點。同時,目前研究發(fā)現(xiàn)以疊氮聚合物為黏合劑的推進劑在高溫下軟化傾向程度增加[2],這將對推進劑抗蠕變性能產(chǎn)生不利影響?;谏鲜龇治觯狙芯酷槍AP-ETPE基推進劑進行了不同溫度和加載應(yīng)力水平的拉伸蠕變實驗,對GAP-ETPE基推進劑的蠕變過程和機制進行分析,開展了對GAP-ETPE基推進劑長期蠕變行為的預(yù)測,為GAP-ETPE基高能固體推進劑的應(yīng)用提供參考。
GAP基含能熱塑性彈性體(GAP-ETPE),自制,數(shù)均分子質(zhì)量為 26500g/mol;Ⅲ類AP,平均粒徑125μm,大連氯酸鉀廠;8類RDX,平均粒徑 40μm,甘肅銀光化學(xué)工業(yè)集團有限公司;Ⅳ類鋁粉(Al),平均粒徑4μm,湖南金天鋁業(yè)有限公司;以上固體填料使用前均在60℃恒溫水浴烘箱中干燥7天。四氫呋喃,分析純,北京通廣精細(xì)化工公司,直接使用。GAPA,數(shù)均分子質(zhì)量為767g/mol,黎明化工研究設(shè)計院有限責(zé)任公司。
壓延混合采用廣東省東莞市錫華精密儀器有限公司的XH-401C油加熱開煉機。雙輥壓延參數(shù)為:雙輥溫度60℃;前后輥的剪切速率分別為5.8、5.1r/s,輥速比為1.14;開煉壓延時間為20min。
蠕變性能測試采用法國Metrovib 1000+動態(tài)黏彈譜儀。推進劑試樣制成矩形薄片,有效尺寸為15mm×10mm×2mm;采用拉伸模式分別測試樣品在不同溫度和應(yīng)力水平下的蠕變及其恢復(fù)。實驗前預(yù)平衡20min,蠕變及恢復(fù)時間均為40min,加載時間為0.4s。
制備固含量為80%的GAP-ETPE基固體推進劑,根據(jù)能量優(yōu)化的結(jié)果,包含以下組分(質(zhì)量分?jǐn)?shù)):19.5%GAP-ETPE,0.5%GAPA,7%AP,25%RDX,18%Al。其理論比沖達(dá)到269.2s,氧系數(shù)為0.45。按照配方制備的固體推進劑實測爆熱6920kJ/kg,實測密度為1.72g/cm3。
制備工藝為:將計量好的GAP-ETPE以及增塑劑GAPA加入四氟杯中,加入四氫呋喃,在60℃下攪拌溶解;待完全溶解后依次加入固體填料Al、RDX和AP;60℃下攪拌揮發(fā)溶劑。待體系不再流動后放入60℃真空烘箱中保持負(fù)壓24h繼續(xù)去除溶劑;將混合體系放入雙棍開煉機中壓延20min,將得到的片狀固體推進劑負(fù)重壓平冷卻。按照尺寸要求制備測試試樣,如圖1所示。靜置24h去除加工過程中的熱應(yīng)力。
圖1 GAP-ETPE基推進劑實物圖Fig.1 Pictures of GAP-ETPE-based propellant
不同加載應(yīng)力對GAP-ETPE基固體推進劑的蠕變性能的影響如圖2所示。由圖2可以看出,GAP-ETPE基推進劑表現(xiàn)出3個蠕變階段[10]:在第一階段,蠕變應(yīng)變迅速增加,形變瞬時產(chǎn)生,這一階段對應(yīng)的是鍵長鍵角的瞬時響應(yīng);在第二階段,蠕變應(yīng)變繼續(xù)增加,而蠕變應(yīng)變速率持續(xù)降低,這一階段對應(yīng)的是一系列具有不同松弛時間的分子鏈鏈段的運動;在第三階段,蠕變應(yīng)變速率達(dá)到恒定,蠕變應(yīng)變隨時間線性增加,這一階段對應(yīng)的主要是具有極長松弛時間的鏈段,其構(gòu)象隨時間進行緩慢調(diào)整,以及產(chǎn)生了部分分子鏈的滑移。圖2中小圖給出了從室溫單軸拉伸實驗中獲得的推進劑應(yīng)力—應(yīng)變曲線,由圖可見,本研究所進行的蠕變實驗均在GAP-ETPE基推進劑線性黏彈性區(qū)域[11]。
圖2 不同加載應(yīng)力下GAP-ETPE基推進劑蠕變曲線Fig.2 The creep curves of GAP-ETPE-based propellant at different loading stresses
由圖2還可以看出,隨著加載應(yīng)力水平的增加,推進劑蠕變應(yīng)變顯著增加,在相同的應(yīng)力增長幅度下,蠕變應(yīng)變的增長幅度越來越大,推進劑更遲進入恒定蠕變速率階段,同時形變恢復(fù)能力顯著降低。具體的蠕變實驗結(jié)果見表1。
表1 不同加載應(yīng)力下GAP-ETPE基推進劑蠕變實驗參數(shù)Table 1 The creep experimental parameters of GAP-ETPE-based propellant at different loading stresses
由表1可以看出,隨著加載應(yīng)力水平的增加,蠕變應(yīng)變由0.4%提高到4.9%,并且恒定蠕變速率顯著提高;其主要的原因為:一方面,外力的作用效果相當(dāng)于對高分子鏈段的運動起到活化作用,降低鏈段沿外力作用方向運動的活化能,從而沿外力作用方向鏈段運動能力增強,蠕變應(yīng)變與恒定蠕變速率均隨著應(yīng)力提高而增大;另一方面,較大的應(yīng)力使得黏合劑基體變形嚴(yán)重,推進劑脫濕嚴(yán)重,填料與基體的相互作用減弱,物理交聯(lián)網(wǎng)絡(luò)部分被破壞,蠕變應(yīng)變增加幅度越來越大,殘余應(yīng)變增大,恒定蠕變速率上升。
文獻(xiàn)報道[3]恒定蠕變階段的速率將決定著最終蠕變的結(jié)果。為了描述并預(yù)測不同的環(huán)境應(yīng)力對推進劑蠕變性能的影響,根據(jù)表1建立了加載應(yīng)力對最終恒定蠕變速率的定量關(guān)系,如圖3所示,擬合優(yōu)度(R2)為0.998,得到的關(guān)系式為:
Rs=-8.23×10-7+8.87×10-7exp(σ/0.77)
(1)
式中:Rs為恒定蠕變速率;σ為加載的應(yīng)力, MPa。
結(jié)合拉伸曲線,由式(1)可以看出,延長蠕變時間或者增大加載應(yīng)力,推進劑總會出現(xiàn)蠕變破壞階段。
圖3 不同加載應(yīng)力下GAP-ETPE基推進劑恒定蠕變速率曲線Fig.3 The constant creep rate curve of GAP-ETPE-based propellant at different loading stresses
GAP-ETPE基固體推進劑在不同溫度下的蠕變測試結(jié)果如圖4所示。由圖4可知,隨著溫度提高,蠕變應(yīng)變增加,形變恢復(fù)能力明顯下降,推進劑更遲進入恒定蠕變階段。具體的實驗參數(shù)見表2。
圖4 不同溫度下GAP-ETPE基推進劑蠕變曲線Fig.4 The creep curves of GAP-ETPE-based propellant at different temperatures
表2 不同溫度下GAP-ETPE基推進劑蠕變實驗參數(shù)Table 2 The creep experimental parameters of GAP-ETPE-based propellant at different temperatures
由表2可知,蠕變應(yīng)變從1.3%提高到6.1%,恢復(fù)相同時間后殘余應(yīng)變從0.2%提高到了2.3%,同時可以發(fā)現(xiàn),在加工溫度(60℃)以上,GAP-ETPE基推進劑的蠕變性能變化幅度增大。在加工溫度以下時,相同恢復(fù)時間下殘余應(yīng)變(26~46℃)提高了0.3%,而高溫(56~76℃)殘余應(yīng)變提高了1.3%,相差4倍以上。當(dāng)溫度接近加工溫度時,GAP-ETPE基推進劑的抗蠕變性能迅速下降,產(chǎn)生了更多不可逆形變,其主要原因為:一方面溫度的提高使得鏈段自由體積增大,鏈段熱運動能力增強,鏈段松弛時間減小,在更短時間內(nèi)達(dá)到更大的形變;另一方面,溫度提高使得體系內(nèi)各種物理相互作用減弱,相同應(yīng)力下更多物理交聯(lián)點被破壞,鏈段熱運動能力增強,并且更容易發(fā)生分子鏈間的滑移。
同時由表2可以看到,恒定蠕變速率隨溫度先增大后減小,在加工溫度附近恒定蠕變速率出現(xiàn)了極值點,這是因為:一方面溫度提高導(dǎo)致鏈段運動能力提高,松弛時間變短,蠕變速率提高;另一方面,在微觀運動中,分子鏈整鏈的運動也是由于鏈段協(xié)同躍遷引起的[12],因此,在恒定蠕變速率階段,其發(fā)生的不僅僅是分子鏈的滑移,也存在著較長松弛時間鏈段的運動。由圖4還可以看出,隨著溫度的升高,蠕變過程推遲進入恒定蠕變速率階段,第二蠕變階段延長,即隨著溫度升高,原本發(fā)生在恒定蠕變速率階段的運動形式,由于松弛時間變短,更多地出現(xiàn)在了第二階段。并且在加工溫度以上時,結(jié)構(gòu)中大量原本凍結(jié)的運動單元開始緩慢運動,二者綜合導(dǎo)致在恒定蠕變速率階段難以運動的分子鏈占比明顯增大,因此最終達(dá)到的恒定蠕變速率反而減小。
為了更好地描述GAP-ETPE基推進劑的蠕變過程,聯(lián)系宏觀過程與微觀結(jié)構(gòu)變化的規(guī)律,預(yù)測在貯存、運輸及裝載使用過程中推進劑尺寸的形變,選擇合適的黏彈性本構(gòu)模型對GAP-ETPE基推進劑蠕變過程進行分析十分重要。由于實際貯存期間推進劑受到的應(yīng)力較小,約為0.013MPa[5],本研究所進行的一系列蠕變實驗,蠕變應(yīng)變均低于6%,由圖2中小圖可以看出,此時材料仍然處于線性黏彈區(qū)[11],因此,本研究采用線黏彈性本構(gòu)模型分析推進劑在不同環(huán)境溫度與外力載荷下的蠕變過程。其中最常見的是四元件模型,也稱Burgers模型,見圖5(a)。該模型包括3部分:(1)彈性模量為E1的硬彈簧,用來表示鍵長鍵角的變化;(2)模量為E2的軟彈簧與黏度為η2的黏壺并聯(lián)的Kelvin模型用來描述松弛時間為τ2的鏈段運動;(3)黏度為η1的黏壺模擬分子鏈滑移引起的黏性流動。
圖5 蠕變模型示意圖Fig.5 Schematic diagram of creep model
結(jié)合Zhang J B等[3]使用Burgers模型對雙基推進劑的分析結(jié)果,F(xiàn)eng S X等[5]對HTPB推進劑的模擬結(jié)果,以及朱國翠[13]對HTPE推進劑的模擬分析可以發(fā)現(xiàn),該模型對于蠕變初期擬合較差。這是因為實際推進劑體系是非常復(fù)雜的,大量的固體填料以及各類添加劑使得各運動單元實際處在不同的環(huán)境中,尤其是對于熱塑性推進劑,其力學(xué)性能主要依靠的是種類復(fù)雜、強度不同的大量物理交聯(lián)點,因而其力學(xué)松弛過程不止一個松弛時間,而是分布較寬的松弛時間譜[12]。研究發(fā)現(xiàn),在Burgers模型基礎(chǔ)上繼續(xù)串聯(lián)一個Kelvin模型就能使得擬合程度達(dá)到0.995以上,見圖5(b)。其中Kelvin模型Ⅰ模擬受到束縛較小,松弛時間較短的鏈段運動;Kelvin模型Ⅱ模擬受到束縛較大,松弛時間較長的鏈段運動。由胡克定律與牛頓黏性定律可以得出總?cè)渥兞繛椋?/p>
ε(t)=[1/E1+t/η1+(1-e-E2t/η2)/E2+
(1-e-E3t/η3)/E3]σ0
(2)
使用Matlab軟件對GAP-ETPE基推進劑在室溫、不同加載應(yīng)力下以及0.5MPa加載應(yīng)力、不同環(huán)境溫度下的蠕變數(shù)據(jù)進行擬合,擬合結(jié)果如圖6所示,擬合程度均達(dá)到0.995以上,說明廣義的Kelvin模型能夠很好地模擬GAP-ETPE基推進劑的蠕變過程。
圖6 推進劑蠕變數(shù)據(jù)及廣義Kelvin模型擬合曲線Fig.6 Creep data of propellant and fitting curves of generalized Kelvin model
在模型擬合精度非常高且基本一致的前提下,模型中的參數(shù)變化反映了材料宏觀性能的規(guī)律變化。
Ⅰ階段本身表示松弛時間較小的鏈段運動,因而其參數(shù)相比于Ⅱ階段變化不大。不同加載應(yīng)力下廣義Kelvin模型擬合參數(shù)見表3。由表3可以看出,隨著加載應(yīng)力的增大,η1、η3、E1、E3明顯降低,這主要因為基體變形嚴(yán)重,從而導(dǎo)致推進劑發(fā)生脫濕的現(xiàn)象,沿加載方向填料與基體的界面產(chǎn)生微裂紋與微孔洞[11],部分物理交聯(lián)點被破壞,隨著應(yīng)力增大,破壞更加嚴(yán)重,從而模量與黏度降低。同時,隨著應(yīng)力水平的增大,分子鏈沿外力作用方向伸展程度增大,分子鏈間內(nèi)摩擦減弱,黏度下降。并且隨著加載應(yīng)力越大,下降的幅度也更大,這也反映了大應(yīng)力對物理交聯(lián)網(wǎng)絡(luò)的破壞越來越嚴(yán)重。
不同溫度下廣義Kelvin模型擬合參數(shù)見表4。由表4可以看出,隨著溫度升高,模量參數(shù)E1、E2、E3均下降。這主要是因為溫度升高,分子鏈鏈段運動能力增強,使得鏈段更容易沿外力作用方向進行構(gòu)象調(diào)整,導(dǎo)致模量參數(shù)降低。在加工溫度以下(26~56℃),黏度η1、η3降低,這主要是因為溫度提高使得鏈段自由體積增大,同時各種物理相互作用減弱,導(dǎo)致分子間摩擦力降低,黏度參數(shù)降低。而在加工溫度以上 出現(xiàn)了本體黏度參數(shù)η1上升的現(xiàn)象,與前文恒定蠕變速率減小的現(xiàn)象相對應(yīng),原因是在恒定蠕變速率階段分子間滑移比例增加明顯,以及大量原本凍結(jié)的運動單元開始運動,表現(xiàn)為在該階段的運動單元運動能力降低。由此可見,溫度升高,推進劑抗蠕變性能下降嚴(yán)重。
表3 不同加載應(yīng)力下廣義Kelvin模型擬合參數(shù)Table 3 Fitting parameters of generalized Kelvin model at different loading stresses
表4 不同溫度下廣義Kelvin模型擬合參數(shù)Table 4 Fitting parameters of generalized Kelvin model at different temperatures
推進劑的貯存是一個長期的過程,在此期間,由于受到自重或者環(huán)境應(yīng)力的影響,推進劑將發(fā)生長期的力學(xué)松弛行為,表現(xiàn)為其尺寸緩慢發(fā)生變化,長期蠕變甚至?xí)顾幮桶l(fā)生明顯變化,影響燃燒和安全系數(shù)[1],因此,對于推進劑長期蠕變行為的預(yù)測有利于很好地描述其尺寸穩(wěn)定性。但是實際長期蠕變實驗需要耗費大量的人力物力,對實驗條件要求高,本研究借助TTSP原理,將不同溫度下的短期蠕變數(shù)據(jù)進行平移,反映GAP-ETPE基推進劑長期蠕變性能的預(yù)測結(jié)果。
在蠕變實驗中,蠕變?nèi)崃縅是一維情況下線性黏彈材料的蠕變函數(shù),是衡量材料黏彈性的重要參數(shù),也是固體發(fā)動機裝藥結(jié)構(gòu)完整性分析的重要參數(shù)[14],其表達(dá)式如下:
(3)
式中:J(t)為t時刻的蠕變?nèi)崃?;?t)為推進劑在t時刻的應(yīng)變;σ(t)為推進劑在t時刻所受應(yīng)力。
選取室溫26℃為參考溫度,將各柔量曲線平移,得到室溫及加載應(yīng)力為0.5MPa條件下,GAP-ETPE基推進劑的蠕變?nèi)崃恐髑€如圖7所示,其變化符合七次多項式的形式,R2為0.997,其擬合公式為:
J(t)=-8.59×10-4+2.98×10-2lgt-2.51×10-2(lgt)2+1.12×10-2(lgt)3-2.45×10-3(lgt)4+2.73×10-4(lgt)5-1.24×10-5(lgt)6+6.33×10-8(lgt)7
(4)
圖7 GAP-ETPE基推進劑蠕變?nèi)崃恐髑€Fig.7 The creep compliance master curve of GAP-ETPE-based propellant
蠕變?nèi)崃恐髑€將GAP-ETPE基推進劑蠕變的考察時間擴大到了108.1s,相當(dāng)于4年。由圖7可知,在此時間尺度內(nèi),蠕變?nèi)崃咳匀惠^小,GAP-ETPE基推進劑未發(fā)生蠕變破壞。蠕變?nèi)崃啃再|(zhì)與推進劑結(jié)構(gòu)密切相關(guān),通過柔量主曲線能夠預(yù)測推進劑長期蠕變行為以及尺寸穩(wěn)定性,為推進劑結(jié)構(gòu)完整性分析以及裝藥設(shè)計提供支撐。
同樣,在不同溫度下的蠕變?nèi)崃科揭七^程中得到了不同溫度下的位移因子,如圖8所示。在26~76℃范圍內(nèi),GAP-ETPE基推進劑的位移因子與溫度呈現(xiàn)線性關(guān)系,擬合優(yōu)度R2為0.994,擬合后其函數(shù)式關(guān)系為:
lga(T)=28.828-0.097T
(5)
圖8 GAP-ETPE基推進劑位移因子Fig.8 The displacement factor of GAP-ETPE-based propellant
在針對其他類型的固體推進劑拉伸蠕變研究方面,胡義文等[2]制備了75%固含量、增塑比為1、R值為1.06的熱固性PBT基復(fù)合固體推進劑;朱國翠[12]制備了80%固含量、增塑比為1的HTPE基固體推進劑;王鴻麗制備了改性雙基推進劑[15],NC/NG/RDX/AP/Al/添加劑體系質(zhì)量分?jǐn)?shù)分別為20.5%、20%、32.5%、13.7%、8.4%、3.9%。各自對應(yīng)的蠕變性能如表5所示。
表5 不同種類固體推進劑蠕變性能對比Table 5 Comparison of creep properties of different types of solid propellants
由表5可以看出,GAP-ETPE基推進劑的抗蠕變性能優(yōu)于常見的熱固性復(fù)合固體推進劑,略微低于改性雙基推進劑。這是因為熱固性復(fù)合固體推進劑制備過程中往往需要加入較多的增塑劑以提高工藝性能,否則藥漿黏度較高,推進劑容易產(chǎn)生孔洞缺陷,而大量的增塑劑使推進劑的抗蠕變性能嚴(yán)重下降;改性雙基推進劑中的硝化棉具有較高的模量,這使得推進劑具有較好的抗蠕變性能。同時也可以看出,即使在高溫條件下,GAP-ETPE基推進劑依然具有良好的抗蠕變性能,優(yōu)于熱固性的PBT基復(fù)合固體推進劑。
(1) 較高的加載應(yīng)力與溫度使推進劑抗蠕變性能下降嚴(yán)重;室溫下恒定蠕變速率與加載應(yīng)力具有良好的指數(shù)函數(shù)依賴形式;同時對比發(fā)現(xiàn),GAP-ETPE基推進劑室溫及高溫的蠕變性能優(yōu)于常見的熱固性復(fù)合固體推進劑。
(2) 廣義的Kelvin模型能夠很好地對GAP-ETPE基推進劑的蠕變過程進行描述,擬合結(jié)果均達(dá)到0.995以上,模型中的參數(shù)在一定程度上反映了推進劑蠕變過程中內(nèi)部結(jié)構(gòu)的變化。
(3) 推進劑在4年內(nèi)未發(fā)生蠕變破壞,蠕變?nèi)崃咳韵鄬^小,同時位移因子與溫度呈現(xiàn)線性關(guān)系。這為GAP-ETPE基推進劑的實際應(yīng)用奠定了基礎(chǔ)。