本站小編為你精心準備了NEPE推進劑固化劑的分子模擬參考范文，愿這些范文能點燃您思維的火花，激發您的寫作靈感。歡迎深入閱讀并收藏。

《固體火箭技術雜志》2015年第五期

目前，在nepe推進劑體系中廣泛使用的是聚氨酯粘合劑，雖然只占總質量的6%～8%，但卻起著重要的作用。這類粘合劑主要成分有環氧乙烷-四氫呋喃共聚醚(PET)，由于這種共聚醚主鏈結構的不規整性以及醚鍵的存在，使其具有不易結晶、力學和工藝性能好的特點，是NEPE推進劑的一種較為理想的聚醚多元醇。工程應用中，向粘合劑中加入一定量的固化劑，使低分子聚合物或單體經一定的化學反應生成高分子化合物，或使線型高分子化合物交聯成體型高分子化合物［1］。聚氨酯粘合劑中采用的固化劑主要為多官能度異氰酸酯(N-100)、甲苯二異氰酸(TDI)、異佛爾酮二異氰酸(IPDI)等，固化劑的選擇不僅要注重交聯固化的效果，還要與其他組分具有很好的相容性，這也是固體推進劑成型過程得以實現的前提。相容性直接影響著推進劑的力學性能，精確地表征相容性是研究推進劑共混體系的基礎［3］，實驗方法過程復雜，且耗時耗力。分子模擬技術克服了實驗研究的缺點，成為研究固體推進劑等含能材料微觀性質的一種非常有效的手段。焦東明［4］采用分子動力學(MD)方法計算了丁羥推進劑粘合劑、增塑劑的溶度參數，來考察組分間的相容性，模擬數值與實驗基本一致。蘭艷花等應用MD方法模擬計算了HTPB固體推進劑粘合劑與增塑劑組分DOS、DOA、DBP和DOP所組成共混物的密度、結合能和徑向分布函數。以及純物質的溶度參數，并根據模擬計算結果，評價了HTPB與不同增塑劑相容性的優劣，次序為HTPB/DOS＞HTPB/DOA＞HTPB/DOP＞HTPB/DBP;黃銳、姚維尚等對纖維素基、聚酯類等不同粘合劑進行了模擬計算，為制得力學性能較好的高能粘合劑提供了依據;楊月誠［9］對HTPB粘合劑及其常用的增塑劑、固化劑組分的溶度參數進行了模擬計算，發現相容性對固化效果有一定影響。但以上研究所建模型采用的原子數量較少，計算結果精度較低，且分析方法較單一，不能直觀展示組分相容性。

本文采用MD方法對PET和不同固化劑的純物質及其組成的共混體系進行了模擬計算，針對以上不足，建模時適當增加模型原子數，對溶度參數模擬進行了驗證，提高了計算結果的準確性，在分析徑向分布函數的基礎上，增加了Blends分析方法，不僅計算出混合能，而且通過結合能分布圖更加直觀地展示組分相容性。運用MD方法很容易得到組分溶度參數，而徑向分布函數分析法和Blends方法的結果更加直觀，本文使用這3種分析方法，可從不同的角度去分析組分相容性，其模擬結果能夠為NEPE推進劑及其他含能材料的配方設計提供參考依據。

1基本理論與方法

1．1理論依據從熱力學觀點來看，2種物質混合時，有以下熱力學關系。兩組分能否相容并得到均相共混物，取決于ΔGM，這里討論的相容性是熱力學相容性，即分子水平上的均相體系。只有滿足ΔGM＜0時，兩組分互容，且溶解過程能夠自發進行;溶解中熵總是增加的，即ΔSM＞0，對于推進劑組分內弱極性或非極性高分子體系的混合溶解過程一般是吸熱過程，即ΔHM＞0;為了滿足ΔGM＜0，ΔHM值越小或越趨近于零較好。從式(3)可看出，δ1和δ2的差越小，則ΔHM值越小，越有利于溶解，這就是溶度參數相近原則。實驗表明，對非晶態高分子材料來說，若分子間沒有強極性基團或氫鍵基團，2種組分只要滿足|δ1－δ2|＜1．7～2．0(J/cm3)1/2［12］，則2種組分相容。

1．2模擬方法分子動力學(MD)方法通過計算分子的汽化熱預測材料的溶度特性，首先使用MaterialsVisualizer構建分子的空間構型，采用Discover模塊中的SmartMini-mizer進行幾何優化，而后構建其無定形體，并進行優化，優化方法采用智能最小優化方法，在正則系綜(NVT)條件下，進行250ps的分子動力學模擬，前200ps用于平衡，后50ps用于計算。動力學運算后，采取Analysis來獲取溶度參數，分析徑向分布函數，判斷相容性。其中，力場采用COMPASS力場，能量最優化采用最陡下降法(SteepestDescent)和共軛梯度法(ConjugateGradient)［13］，溫度控制用Anderson方法，非鍵截斷(cutoff)取0．95nm，樣條寬度(splinewidth)取0．1nm，緩沖寬度(bufferwidth)取0．05nm，時間步長1fs，靜電和范德華力計算采取Eward和Atom-based［16］長程加和方法。共混(Blends)方法是通過MonteCarlo模擬后加權分子空間取向排布計算結合能來判斷組分相容性，將優化后的不同組分的分子模型設置為screen和base作為共混模塊輸入條件，在COMPASS力場條件下，提交到mixing任務進行共混模擬，計算各個組分的結合能，通過分析共混結合能分布曲線，判斷組分相容性。

1．3方法驗證使用MD方法模擬計算常見聚合物的溶度參數，以聚乙烯(PE)、聚丙烯(PP)和端羥基聚丁二烯(HT-PB)為例，聚合物的物理性質(298K)和計算出的聚合物溶度參數見表1。將MD模擬計算的結果δMD與實驗數據δexp及焦東明采用的Synthia方法模擬計算得到的結果δSynthia進行對比。可看出，MD方法得到的數值較好吻合實驗值和采用Synthia方法得到的模擬值，且誤差在允許范圍內。因此，可用來模擬計算組分的溶度參數，進而研究組分相容性。

2模擬計算與結果分析

2．1模型構建運用MS軟件中的Visualizer模塊，根據PET、N-100、TDI、IPDI的分子簡式建立相應的分子模型，如圖1所示。利用Amorphouscell模塊，構建上述4種組分的純物質以及PET與固化劑混合體系的無定形分子模型。模型中原子總數約為1000個。

2．2溶度參數分析方法對4種組分純物質的無定形模型進行優化后，進行分子動力學運算并進行分析，組分的物理性質及模擬計算，結果見表2。從表2可看出，對PET、N-100、TDI、IPDI溶度的模擬，能夠定量地得到溶度參數。4種組分的溶度參數值均較接近，經過計算得到溶度參數差值見表3。從表3可看出，PET與3種固化劑的溶度參數差值的絕對值均小于1．7(J/cm3)1/2，根據溶度參數相近原則，說明PET與這3種固化劑能很好的相容，且PET與N-100的相容性好于PET與IPDI的相容性，PET與IPDI的相容性好于PET與TDI的相容性。

2．3徑向分布函數分析方法對PET與3種固化劑組成的共混物進行優化，并進行MD運算后，分析其徑向分布函數g(r)。g(r)是反映材料微觀結構的特征物理量，它表示在一個分子周圍距離為r的地方出現另一個分子的概率密度相對于隨機分布概率密度的比值。分子間徑向分布函數可揭示非鍵原子間相互作用的方式和本質，氫鍵作用范圍為0．26～0．31nm，范德華力(vanderWaals，vdw)作用范圍為0．31～0．50nm［17］，從圖2可看出，4種純物質分子間主要作用方式為vdw作用，PET與3種固化劑的分子間徑向分布函數曲線在vdw作用范圍內距離越近，說明PET與之相容性好，在vdw作用范圍內，距離PET最近的是N-100，其次是IPDI，最遠的是TDI，說明固化劑與PET相容性優劣順序為N-100＞IPDI＞TDI。對于混合體系而言，分子間徑向分布函數也可用來判斷混合物的相容性。混合物的分子間徑向分布函數g(r)越是高于單一物質的徑向分布函數，相容性越好;反之，則會發生相分離［3］。

2．4Blends分析方法在Blends模塊中，將PET分子模型設置為base，將N-100、TDI和IPDI設置為screen，分別組成共混體系，提交Mixing任務進行計算，得到混合能如表4所示，共混結合能分布曲線如圖4所示。從圖3可看出，混合體系分子間主要作用方式也是vdw作用，在PET/N-100混合體系中，PET-(N-100)的分子間g(r)值高于(N-100)-(N-100)和PET-PET的分子間g(r)值，而在PET/TDI混合體系和PET/IPDI混合體系中同組分的分子間g(r)值均高于兩組分分子間g(r)值，但在PET/TDI混合體系中，g(r)差值較大，故可判斷相容性結論為PET/N-100＞PET/IPDI＞PET/TDI。從圖4可看出，PET/N-100共混時，共混結合能Ebb、Ess和Ebs這3個值分布較接近，其次是PET/IPDI共混體系、PET/TDI共混體系，共混結合能分布越接近，體系的混合能就越小。所以，可得出PET與3種固化劑的相容性優劣順序為PET/N-100＞PET/IPDI＞PET/TDI。

3結論

(1)PET與固化劑N-100、TDI及IPDI的相容性優劣順序為PET/N-100＞PET/IPDI＞PET/TDI。從相容性角度考慮，NEPE推進劑粘合劑中，選用N-100作為固化劑，效果較好。(2)溶度參數法簡單易實現，徑向分布函數分析與Blends方法顯示結果直觀明了，這3種方法能夠從不同角度去分析組分相容性，綜合運用這3種方法，可準確直觀地判斷組分間相容性。(3)本文模擬方法可作為預測NEPE推進劑其他組分相容性的有利工具，也可為固體推進劑和炸藥的配方設計提供參考依據。

作者：王廣 侯世海 武文明 李紅霞 單位：第二炮兵工程大學 601 室