不同含氮量NC對CMDB推進劑力學性能的影響

張亞俊，李吉禎，唐秋凡，李 偉，王 可，樊學忠，屈 蓓，張正中，鄢海濤，鮑遠鵬

(1.西安近代化學研究所，陜西 西安 710065；2.宜賓北方川安化工有限公司，四川 宜賓 644219)

引 言

改性雙基(CMDB)推進劑裝藥具有特征信號(煙霧、尾焰)低的顯著優點，能夠有效降低對激光、紅外制導信號的干擾，是現役戰術武器的重要推進劑[1-3]。目前，CMDB推進劑中唯一成熟使用的黏合劑是半剛性的大分子硝化棉(NC)，但其玻璃化溫度較高，低溫容易脆化[4]。當推進劑中加入了大量高能炸藥RDX、HMX及燃燒劑Al粉等固體填料后，其力學性能會下降，“高溫變軟、低溫變脆”的問題更加突出[5]，極大地限制了CMDB推進劑特別是高能CMDB推進劑的發展和應用[6]。

目前，解決上述問題的主要途徑有NC改性和包覆[7-9]以及對高能填料進行適當表面改性和包覆[10]。20世紀70～80年代，前蘇聯科研人員為了改善纖維素大分子的柔順性，先制備具有分子內增塑支鏈的多羥基纖維素醚，再硝化合成熱塑性好、能量高的新型纖維素基含能黏合劑，并曾多次報道纖維素醚的硝化物能改善推進劑的力學性能[11-12]。近年來，張有德、吳艷光等[13-14]也對多種纖維素醚硝酸酯、GAP型NC進行了大量研究，發現許多改性后的NC能夠改善推進劑力學性能。丁海琴等[15]曾用溶液共混法獲得GAPE增韌NC共混體系，其拉伸強度可達43.4MPa，延伸率達33.5%。可見，NC的改性一直是CMDB推進劑力學性能研究的熱點[16]，但由于缺乏相匹配的理論研究，針對NC表面的修飾及改性都是通過大量實驗來實現，費時費力。

因此，本實驗從最基礎的NC氮含量展開研究，采用分子動力學模擬不同含氮量的NC對CMDB推進劑力學性能的影響，通過計算彈性常數、回轉半徑和徑向分布函數等參數，從分子水平上認識不同含氮量的NC與推進劑其他組分相互作用的實質，為探索結構更復雜的NC對推進劑力學性能的影響奠定基礎。同時結合拉伸試驗，測試含有不同含氮量NC的CMDB推進劑的力學性能，為理論研究提供數據支撐，并為改善CMDB推進劑的力學性能奠定數據基礎。

1 分子動力學模擬

1.1 物理建模

依據NG、NC和HMX的化學結構式，采用美國Accelrys公司開發的Materials Studio 7.0 軟件的Visualizer 模塊和Amorphous Cell 模塊，建立相應的分子物理模型[17-18]，各組分分子模型如圖1所示，其中，碳原子為灰色，氫原子為白色，氮原子為藍色，氧原子為紅色。

1.2 MD模擬細節

用Smart Minimizer方法對所構建的模型進行能量優化，然后選用COMPASS力場[19-21]，利用Smart Minimizer方法對其進行幾何構型優化[22]。應用周期性邊界條件，即以立方元胞為中心，周圍有26個相鄰的鏡像立方元胞，以達到利用較少分子模擬宏觀性質的目的。各分子起始速度由Maxwell-Boltzmann隨機分布給定，用Velocity Verlet算法進行求解。對分子間的范德華(vdW)和靜電(coulomb)作用力計算分別采用Atom-based方法和Ewald方法，非鍵截取半徑0.95nm，樣條寬度(Spline width)取0.1nm，緩沖寬度(Buffer width)取0.05nm。

構建等質量比的NC/NG/HMX共混體系模型，其中NC分子鏈(每條分子鏈有20個聚合單元組成)，數目為3條(分別標記為紅色、藍色和黃色)，考慮實際推進劑配方比，NC、NG和HMX質量比為26∶33∶33，混合模型的結構如圖2所示。采用Smart Minimizer方法對上述混合體系模型進行5000步能量最小化優化。接著，采用Forcite模塊，在系綜為NPT系綜(101.325kPa)，分別在293K和323K條件下，進行1 000ps的分子動力學模擬以獲得平衡密度，時間步長為1fs。

2 實 驗

2.1 原料與儀器

1#NC(含氮量12.5%～12.7%)、2#NC(含氮量11.88%～12.40%)、3#NC(含氮量11.75%～12.10%)，四川北方硝化棉股份有限公司；NG，西安近代化學研究所，阿貝爾安定性試驗合格；HMX，粒徑(d50)為12～25μm，甘肅白銀銀光化學材料廠。

INSTRON萬能材料試驗機，得米特(天津)科技有限公司；行星式捏合機，湖北航鵬化學動力科技有限公司。

2.2 樣品制備

以含質量分數26%的3#NC球形藥的CMDB推進劑配方為基礎，其他組分含量不變，分別以1#NC球形藥和2#NC球形藥部分取代3#NC球形藥(分別取代質量分數2%、4%、6%、8%和10%的3#NC球形藥)，研究不同含氮量的球形藥含量對CMDB推進劑力學性能的影響。推進劑樣品均采用淤漿澆鑄工藝制備。將質量分數33%NG、26%NC、33%HMX及8%其他組分在2立升行星式捏合機中混合1h左右，出料后經70℃固化3d，退模。

2.3 力學性能測試

按照GJB770B-2005 方法431.1 單向拉伸法進行推進劑拉伸試驗。將厚度為10mm的片狀推進劑制成啞鈴形狀，在20℃和50℃條件下，在制備好的試樣縱軸方向施加靜態單向拉伸載荷，拉伸速率100mm/min，測試其所能承受的最大拉伸強度及相應的延伸率。

3 結果與討論

3.1 不同含氮量NC與NG/HMX共混體系彈性常數模擬

根據彈性體的統計力學原理[23-24]，胡克定律通常可寫為以下公式:

(1)

在各向同性體系中，具有兩個獨立的彈性常數，C11和C12。使C12=λ，C11-C12=μ。Cij記為如下形式：

(2)

式中：λ和μ為常數。

拉伸模量(E)、剪切模量(G)、體積模量(K)和泊松比(υ)見式(3)～(6)：

(3)

G=μ

(4)

(5)

(6)

不同含氮量NC與NG/HMX共混體系彈性常數和各向同性的力學性能參數計算值見表1。

表1 不同含氮量NC與NG/HMX共混體系彈性常數和力學性能參數Table 1 Elastic coefficient and mechanical property parameters of NC/NG/HMX blends with different content of nitrogen

由表1可以看出，不同含氮量NC與NG/HMX共混體系的拉伸模量、剪切模量和體積模量與純NC相比均有所降低，這表明增塑劑可降低NC的剛性并增加其韌性。同時，與純NC相比，彈性常數的變化趨勢大小為NC(w(N)=11.8%)/NG>NC(w(N)=12.6%)/NG>NC(w(N)=12.2%)/NG。柯西壓值隨著NG的加入而減小，表明NG的加入能夠降低NC的脆性，且體系延展性變好。含氮量較低的NC/NG體系，柯西壓值降低更加明顯，NC(w(N)=11.8%)體系的柯西壓值從純NC的2.06降至0.560。說明在相同條件下含氮量11.8%的NC更能改善混合體系的彈性常數。

3.2 不同含氮量NC在共混體系的回轉半徑計算

由于NC分子含有不同的含氮量，體系中其他組分與其之間的相互作用必然會受到影響，因此通過回轉半徑來研究不同體系中NC分子鏈構象受不同含氮量NC分子影響的變化。不同NC/NG/HMX共混體系模型中，NC的回轉半徑(Rg)見表2。

表2 不同NC/NG/HMX共混體系模型中NC的回轉半徑Table 2 Radius of gyration in different NC/NG/HMX blend models

由表2中數據可知，含不同含氮量的NC分子鏈的回轉半徑隨著體系中NG和HMX的加入而逐漸增大，這是由于增塑劑的加入促進NC鏈的伸展，導致高分子鏈尺寸增大，回轉半徑增加。同時對比不同NC/NG/HMX混合體系中NC的回轉半徑，發現隨著NC分子中含氮量的增加，回轉半徑也逐漸增大。NC是由葡萄糖酐環狀殘基組成的線性聚合物，具有一定的剛性，但其環間的醚鏈又使鏈節的內旋轉比較容易，所以NC分子鏈具有一定的內旋自由度，但這種內旋自由度會受到NC分子鏈中各原子間作用力的抑制。因此，推斷上述NC分子鏈回轉半徑增加的原因，可能是由于NC分子鏈中硝酸酯基團的增加(即NC體系中含氮量的增加)，使其內旋自由度增大，內旋轉越容易，NC分子的柔性增加，NC分子鏈尺寸增大所致。已有研究表明[25]，NG分子對NC分子鏈內部的作用力有弱化作用，且隨著NG分子數量的增加，弱化作用增強。

3.3 徑向分布函數模擬

通過分析25℃下H-N以及H-O之間的徑向分布函數(RDF)，研究不同含氮量NC和NG之間的相互作用。NC分子中的H用H(1)表示，NG中的O和N原子分別用O(2)和N(2)表示。RDF曲線如圖3所示。

一般來說，氫鍵的作用距離是0.20～0.31nm，強的范德華力作用距離在0.31～0.50nm之間，弱的范德華力作用距離大于0.50nm。在圖3中，所有曲線的第一個峰出現在0.20～0.25nm，說明NC分子中的H和NG分子中的N以及O之間都存在氫鍵作用。且圖中混合體系NC(w(N)=11.8%)/NG/HMX的H(NC)和N(NG)的峰值明顯高于NC(w(N)=12.2%)/NG/HMX和NC(w(N)=12.6%)/NG/HMX體系；混合體系NC(w(N)=11.8%)/NG/HMX的H(NC)和O(NG)的峰值略高于NC(w(N)=12.2%)/NG/HMX和NC(w(N)=12.6%)/NG/HMX體系。這表明在理論模擬的含氮量范圍內，無論是NG分子中的N還是NG分子中的O都與NC分子中H之間的相互作用，隨著NC氮含量的降低而增強。即在計算范圍內，NC氮含量的降低對NC分子與NG分子間的作用力有強化作用，且含氮量低的NC與NG的相互作用更強。

3.4 不同含氮量的NC對推進劑力學性能的影響

由拉伸試驗測得1#NC和2#NC含量對CMDB推進劑力學性能的影響結果見表3和表4。

表3 1#NC含量對推進劑拉伸強度和延伸率的影響Table 3 Effect of 1#NC content on the tensile strength and elongation of propellants

表4 2#NC含量對推進劑的拉伸強度和延伸率的影響Table 4 Effect of 2#NC content on the tensile strength and elongation of propellants

由表3可見，隨著1#NC球形藥取代3#NC球形藥的含量增加，CMDB推進劑常溫和高溫(20℃和50℃)拉伸強度(σm)增大，而延伸率(εm)減小；當1#NC球形藥取代3#NC球形藥的質量分數超過8%時，推進劑常溫和高溫的拉伸強度減小，延伸率繼續減小。由表4可見，隨著2#NC球形藥取代3#NC球形藥的含量的增加，CMDB推進劑常溫和高溫(20℃和50℃)的拉伸強度(σm)和延伸率(εm)增大；當2#NC球形藥取代3#NC球形藥的質量分數超過4%時，推進劑的拉伸強度(σm)和延伸率(εm)反而減小。

上述結果表明，在澆鑄CMDB推進劑中添加適量的1#NC球形藥和2#NC球形藥均有助于提高澆鑄CMDB推進劑拉伸強度。結合理論研究結果，在計算模擬的含氮量范圍內，含氮量越低，NC與NG相互作用越強，同時對NC/NG/HMX共混體系的彈性系數改善效果越好。但是，含氮量高的NC分子，可能對NC分子鏈內部的作用力均有弱化作用，使其內旋自由度增大，回轉半徑增加，可促進NC分子的延展。因此用1#NC球形藥和2#NC球形藥取代3#NC球形藥后，隨著1#NC球形藥和2#NC球形藥的含量增加，拉伸強度顯示先增加后減少的現象。這種現象是NC分子鏈內部作用與NC與NG分子間作用力共同影響的結果。

4 結 論

(1)理論研究表明，NC體系中含氮量的增加，對NC分子鏈內部的作用力有弱化作用，使其內旋自由度增大，NC分子鏈尺寸增大。表明含氮量的增加可以改善NC的延展性。

(2)徑向分布函數計算結果顯示，NC含氮量越低，NC與NG的相互作用越強，共混體系NC(w(N)=11.8%)/NG/HMX中NC與NG的相互作用最強。

(3)試驗結果表明,在澆鑄CMDB推進劑中添加適量的1#NC球形藥和2#NC球形藥均有助于提高澆鑄CMDB推進劑的拉伸強度。