嚴(yán)啟龍，宋振偉，安 亭，張曉宏，趙鳳起

(西安近代化學(xué)研究所燃燒與爆炸技術(shù)重點(diǎn)實(shí)驗(yàn)室，陜西 西安 710065)

?

含能材料物理化學(xué)性能理論預(yù)估研究進(jìn)展

嚴(yán)啟龍，宋振偉，安 亭，張曉宏，趙鳳起

(西安近代化學(xué)研究所燃燒與爆炸技術(shù)重點(diǎn)實(shí)驗(yàn)室，陜西 西安 710065)

從含能材料領(lǐng)域的最近發(fā)展成果出發(fā)，討論了該領(lǐng)域的主要研究方向,重點(diǎn)論述了當(dāng)前含能材料物理化學(xué)性能理論預(yù)估的最新成果，主要包括量子化學(xué)、分子動(dòng)力學(xué)或者半經(jīng)驗(yàn)QSPR建模的方法預(yù)估含能材料的感度、燃燒爆轟性能、反應(yīng)活性、固化機(jī)制與力學(xué)性能的研究進(jìn)展?？偨Y(jié)了目前存在的主要技術(shù)壁壘，包括缺乏完備統(tǒng)一的含能材料性能標(biāo)準(zhǔn)實(shí)驗(yàn)數(shù)據(jù)庫，沒有自主知識(shí)產(chǎn)權(quán)的商業(yè)化含能材料性能計(jì)算軟件，且國(guó)際上商業(yè)軟件對(duì)含能材料的物理化學(xué)性能的可靠預(yù)測(cè)僅局限于爆轟性能和燃燒性能。文獻(xiàn)調(diào)研表明，我國(guó)需要進(jìn)一步加強(qiáng)該領(lǐng)域研究，最終建立一個(gè)能評(píng)價(jià)含能材料性能與安全的綜合軟件平臺(tái)。附參考文獻(xiàn)90篇。

含能材料；物理化學(xué)性能；理論計(jì)算；量子化學(xué)；分子模擬

引 言

鑒于國(guó)防的重要性，世界各國(guó)對(duì)研發(fā)新型含能材料進(jìn)行了長(zhǎng)期持續(xù)的投入。至20世紀(jì)70年代，已基本形成了獨(dú)立的含能材料學(xué)科體系(國(guó)內(nèi)專業(yè)名稱為軍事化學(xué)與煙火技術(shù))，并成為軍械武器系統(tǒng)技術(shù)的重要支撐。尤其是自20世紀(jì)90年代以來，隨著世界局勢(shì)的緊張變化及軍事工業(yè)的發(fā)展，作為武器系統(tǒng)最重要的動(dòng)力和毀傷能量源的含能材料得到了迅猛發(fā)展。調(diào)查表明[1]，近年比較活躍的含能材料研究方向包括：含能化合物的分子設(shè)計(jì)、含能材料納米化、共晶化與表面改性、含能材料熱分解、點(diǎn)火燃燒和爆轟、極限條件下含能材料的響應(yīng)、含能材料配方設(shè)計(jì)及工程化應(yīng)用、含能材料安全性能與綠色利用及含能材料物理化學(xué)性能理論預(yù)估(包括量化計(jì)算、多尺度分子模擬)等。本文簡(jiǎn)要總結(jié)了含能材料物理化學(xué)性能理論預(yù)估方面的最新研究進(jìn)展。

1 含能材料感度預(yù)估

含能材料在受到撞擊、摩擦、靜電火花或沖擊波等刺激時(shí)，容易發(fā)生燃燒或爆炸[2]。撞擊感度是評(píng)判炸藥安全性的主要指標(biāo)，在理論和工程應(yīng)用中具有重要意義[3]。解決炸藥感度理論判據(jù)問題，僅靠量子化學(xué)的方法是不夠的[4]。炸藥在外界刺激作用下發(fā)生爆炸是一個(gè)十分復(fù)雜的過程，涉及機(jī)械、物理和化學(xué)等諸多因素[5]。在同類刺激源作用下，由于炸藥的摩擦系數(shù)、彈塑性、硬度和模量等性質(zhì)不同，炸藥所吸收的機(jī)械功也不同。即使在炸藥吸收的機(jī)械功相同的情況下，由于炸藥的熔點(diǎn)、熔化焓、比熱、導(dǎo)熱系數(shù)等參數(shù)不同，炸藥內(nèi)所產(chǎn)生的熱點(diǎn)溫度也不同。只有在熱點(diǎn)臨界尺寸、溫度，以及持續(xù)時(shí)間相同的情況下，炸藥發(fā)生爆炸的難易程度才取決于該分子的反應(yīng)活性，即可通過量子化學(xué)計(jì)算的結(jié)構(gòu)參數(shù)和熱力學(xué)數(shù)據(jù)表征。然而，科研人員最感興趣的是決定安全性能的分子反應(yīng)活性[6-8]。肖鶴鳴課題組[9]成功運(yùn)用量子化學(xué)方法，完成諸多開拓性的工作。近來他們又對(duì)含能混合體系分子間的相互作用展開了探索，以期對(duì)高聚物黏結(jié)炸藥或固體推進(jìn)劑的配方設(shè)計(jì)提供一定的理論指導(dǎo)[10]。

撞擊感度通常采用落錘實(shí)驗(yàn)法獲得，它以在常壓室溫(20℃)下、特定質(zhì)量的落錘作用下含能材料的爆炸概率(百分比)或50%爆炸概率下的特性落高H50來表征(勢(shì)能值)[11]。然而此類實(shí)驗(yàn)存在一定局限性：(1)實(shí)驗(yàn)具有危險(xiǎn)性；(2)實(shí)驗(yàn)結(jié)果受外界條件和人為因素影響，重復(fù)性差；(3)無法通過實(shí)驗(yàn)獲得新含能材料的撞擊感度。因此，完全靠實(shí)驗(yàn)來確定含能材料的H50已不能滿足日益劇增的新型含能分子設(shè)計(jì)的需要[12]，有必要根據(jù)理論方法對(duì)含能材料的撞擊感度進(jìn)行預(yù)估。首先，物質(zhì)結(jié)構(gòu)可決定其性質(zhì)，同時(shí)物質(zhì)性質(zhì)可反映其分子結(jié)構(gòu)。尋求撞擊感度與炸藥分子結(jié)構(gòu)參數(shù)之間的關(guān)聯(lián)已成為當(dāng)前炸藥撞擊感度理論研究的一個(gè)重要方向[13]。國(guó)內(nèi)最初的研究主要以硝基含能材料等簡(jiǎn)單分子的撞擊感度預(yù)測(cè)為主，提出相關(guān)的計(jì)算原理和方法[14]。國(guó)際上有關(guān)撞擊感度與結(jié)構(gòu)的研究最早始于20世紀(jì)中葉，Bowden等[15]對(duì)含能材料撞擊感度與其晶體結(jié)構(gòu)之間的關(guān)系作了探討。Delpuech等[16]首先發(fā)現(xiàn)了仲硝基類含能材料的沖擊波感度和熱穩(wěn)定性與其分子電子結(jié)構(gòu)、以及C-NO2或N-NO2鍵能之間存在重要的關(guān)聯(lián)性。Kamlet等[17]則提出了一種基于氧平衡估算撞擊感度的方法。通過多年的發(fā)展，形成了以下幾種主要的撞擊感度理論預(yù)測(cè)方法[18]。

1.1 量子化學(xué)方法(QM)

量子力學(xué)理論的不斷完善，尤其是自洽場(chǎng)方法與密度泛函理論(DFT)的建立與完善，及高速計(jì)算技術(shù)的發(fā)展，使得人們能夠借助高水平量子化學(xué)方法在微觀結(jié)構(gòu)層面上研究物質(zhì)結(jié)構(gòu)與性質(zhì)的內(nèi)在聯(lián)系。撞擊感度的量子化學(xué)理論研究也隨之得到了發(fā)展。Politzer等[19]通過對(duì)CaHbNcOd炸藥量子化學(xué)計(jì)算，發(fā)現(xiàn)靜電勢(shì)等分子結(jié)構(gòu)參數(shù)與其撞擊感度存在一定的相關(guān)性。近期，他們又發(fā)現(xiàn)撞擊感度與其晶格的可壓縮性或晶體內(nèi)自由體積存在重大關(guān)聯(lián)，同時(shí)也證實(shí)了關(guān)于含能材料晶體空穴受沖擊壓縮時(shí)產(chǎn)生熱點(diǎn)的起爆理論[20]。Keshavarz等[21]在此基礎(chǔ)上也進(jìn)行了相關(guān)研究，發(fā)現(xiàn)C-NO2鍵區(qū)域的靜電勢(shì)在一定程度上反映了其不穩(wěn)定性，從而用于標(biāo)識(shí)其敏感度[22-23]。Liu等[24]則證實(shí)了硝基化合物中硝基的電荷值決定了高能材料的機(jī)械感度。Mullay等[25]也發(fā)現(xiàn)不含羥基的18種硝基芳香化合物的撞擊感度與C-N鍵的近似靜電勢(shì)值存在較大相關(guān)性。Ren[26]進(jìn)一步發(fā)現(xiàn)可由分子的靜電勢(shì)獲得環(huán)狀結(jié)構(gòu)炸藥分子，如硝基環(huán)丙烷、環(huán)丁烷硝基、硝基環(huán)戊烷和硝基環(huán)己烷等的撞擊感度。Rice和Hare[27]則選取了化學(xué)鍵中點(diǎn)處的靜電勢(shì)近似值作為關(guān)聯(lián)值，用以計(jì)算CaHbNcOd含能分子的撞擊感度和爆熱。而Politzer等[28]則認(rèn)為撞擊感度與含能材料的理論最大爆熱存在必然聯(lián)系，而與爆速、爆壓的關(guān)聯(lián)性小。而根據(jù)熱點(diǎn)起爆理論，所有的失控化學(xué)反應(yīng)都始于熱點(diǎn)引發(fā)的分解反應(yīng)。Zohari等[29]的研究表明，CaHbNcOd系列含能分子的撞擊感度不僅與H/O摩爾比有關(guān)，還與熱分解活化能存在明確的關(guān)系。依據(jù)這一觀點(diǎn)，Mathieu等[30]通過分解反應(yīng)速率常數(shù)估算了硝基化合物的H50值，所得結(jié)果與實(shí)驗(yàn)值非常接近(相關(guān)系數(shù)約為0.8)。該結(jié)果表明含能材料的感度取決于在熱點(diǎn)分散前分解反應(yīng)的自蔓延能力。此外，Tan等[31]研究表明，相對(duì)上述決定性因素，取決于含能材料分子的化學(xué)鍵與非鍵耦合(應(yīng)變能)的分子剛度對(duì)其感度的影響最大。結(jié)合以上多種因素，Keshavarz[32]開發(fā)了一套可以計(jì)算含能材料機(jī)械感度的Visual Basic程序，對(duì)硝基吡啶、硝基咪唑、硝基吡唑、硝基呋咱、硝基三唑、硝基嘧啶、多硝基芳烴、苯并呋咱、硝胺，硝酸脂、含其他官能團(tuán)硝基脂肪族和硝酸高能化合物的撞擊感度計(jì)算精度較高。他們的預(yù)測(cè)結(jié)果對(duì)Mullay等[33]關(guān)于CaHbNcOd炸藥撞擊感度與分子內(nèi)部電荷不平衡程度的相關(guān)性模型給予了支持。Mullay的表面靜電勢(shì)參量模型有5個(gè): 模型1采用每個(gè)鍵中點(diǎn)的靜電勢(shì)的近似值來計(jì)算感度；模型2應(yīng)用等靜電勢(shì)面上正電荷與負(fù)電荷平均值的差值來計(jì)算感度；模型3應(yīng)用與靜電勢(shì)相關(guān)的統(tǒng)計(jì)參量(平衡參數(shù)ν)來計(jì)算感度；模型4運(yùn)用單分子量子化學(xué)信息來估算其爆熱Qdet，然后通過爆熱來計(jì)算撞擊感度；模型5用平衡參數(shù)ν與爆熱結(jié)合來計(jì)算撞擊感度；對(duì)于硝胺化合物，在熱源、沖擊波和機(jī)械撞擊所引發(fā)的分解過程中，雖然在一些情況下會(huì)存在其他起主導(dǎo)作用的反應(yīng)路徑，但N-NO2的斷裂仍然看做初始反應(yīng)步驟。Edwards等[34]同樣采納了模型4，并輔以PM3和DFT兩種級(jí)別量子化學(xué)方法來計(jì)算幾種硝胺炸藥的爆熱。他們發(fā)現(xiàn)，在DFT理論水平，感度隨著最高占有軌道(HOMO)和最低空軌道(LOMO)能量的增加呈指數(shù)遞減。張朝陽等[35]也在DFT理論計(jì)算的基礎(chǔ)上，發(fā)現(xiàn)撞擊感度與硝基所帶電荷之間存在較大的相關(guān)性。他們采用廣域梯度近似(GGA)的方法，基于Beck 混合泛函計(jì)算了硝基上的Mulliken電荷，并與硝基化合物的撞擊感度進(jìn)行關(guān)聯(lián)。結(jié)果表明，當(dāng)硝基上的負(fù)電荷小于0.23 時(shí)，該化合物較為敏感，即H50≤40 cm (2.5 kg)。因?yàn)橄趸鶐щ姾芍悼捎脕砉浪沔I能、氧平衡和分子靜電勢(shì)等一系列結(jié)構(gòu)參數(shù)，且硝基上的Mulliken凈電荷越多，該分子就越鈍感。需要指出的是，他們的方法僅適用于含有弱鍵C-NO2、N-NO2或O-NO2的硝基化合物。

1.2 定量結(jié)構(gòu)-性質(zhì)相關(guān)性法(QSPR)

QSPR通過選用合適的分子結(jié)構(gòu)來描述分子的結(jié)構(gòu)特征，結(jié)合各種統(tǒng)計(jì)建模工具，研究有機(jī)物的結(jié)構(gòu)與其各種物理化學(xué)性質(zhì)之間的定量關(guān)系。分子結(jié)構(gòu)可用反映其特征的各種參數(shù)來描述，即有機(jī)物的各類性質(zhì)都可以用化學(xué)結(jié)構(gòu)的某個(gè)函數(shù)來表示。通過對(duì)分子結(jié)構(gòu)參數(shù)和所研究性質(zhì)的實(shí)驗(yàn)數(shù)據(jù)之間的內(nèi)在定量關(guān)系進(jìn)行關(guān)聯(lián)，建立分子結(jié)構(gòu)參數(shù)和性質(zhì)之間的關(guān)系模型?？煽康亩拷Y(jié)構(gòu)-性質(zhì)相關(guān)模型可用來預(yù)測(cè)尚未合成的化合物的各種性質(zhì)。Nefati等[36]首次嘗試了用人工神經(jīng)網(wǎng)絡(luò)方法來預(yù)測(cè)含能材料的撞擊感度。他們選取204個(gè)含能材料分子作為樣本集，同時(shí)設(shè)置了3類共39個(gè)描述參量來表征其分子結(jié)構(gòu)(包括拓?fù)鋮⒘俊缀螛?gòu)型參數(shù)和電子參數(shù))。通過計(jì)算這204個(gè)含能化合物的39個(gè)結(jié)構(gòu)參數(shù)，并對(duì)它們進(jìn)行自由組合。然后經(jīng)過多元線性回歸(MLR)、偏最小二乘法(PLS)和BP神經(jīng)網(wǎng)絡(luò)法(BP-ANN)等算法，最終確立了撞擊感度的預(yù)測(cè)模型。研究結(jié)果表明，相對(duì)傳統(tǒng)的線性方法，使用非線性神經(jīng)網(wǎng)絡(luò)方法(MLR和PLS)可獲得更優(yōu)化的模型。最佳的神經(jīng)網(wǎng)絡(luò)模型共采用13個(gè)參數(shù)作為輸入神經(jīng)元，包括隱含層的2個(gè)神經(jīng)元。在Nefati等人的研究基礎(chǔ)上，Cho等[37]做了進(jìn)一步的優(yōu)化和改進(jìn)，并預(yù)測(cè)了234種含能化合物的撞擊感度。他們選取了與文獻(xiàn)[36]不同的參量來描述含能化合物的結(jié)構(gòu)，并根據(jù)參數(shù)種類和不同組合將它們分成7個(gè)子集。通過構(gòu)建3層BP神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)，對(duì)每個(gè)描述參量子集進(jìn)行建模。結(jié)果發(fā)現(xiàn)，最好的網(wǎng)絡(luò)結(jié)構(gòu)為17-2-1，即采用含分子組成及拓?fù)漕愋偷?7種分子參量作為輸入層變量，神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)隱含層需包含2個(gè)神經(jīng)元。他們指出，包含分子組成及拓?fù)涿枋龇淖蛹群须娮訁?shù)如LUMO、HOMO和偶極矩的子集能獲得更精確的預(yù)測(cè)結(jié)果。Wang等[38]則進(jìn)一步利用遺傳算法和基于電拓?fù)鋺B(tài)指數(shù)的人工神經(jīng)網(wǎng)絡(luò)方法的QSPR模型預(yù)測(cè)了非雜環(huán)硝基化合物的撞擊感度，獲得的最佳BP神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)為16-12-1，預(yù)測(cè)結(jié)果與實(shí)測(cè)值最接近。

隨后，Keshavarz等[39]僅選取了10個(gè)比較重要的分子結(jié)構(gòu)特征參數(shù)，并利用神經(jīng)網(wǎng)絡(luò)算法通過MATLAB編程，預(yù)測(cè)了大量CaHbNcOd炸藥分子的感度。結(jié)果表明，該模型通過275個(gè)實(shí)驗(yàn)樣本訓(xùn)練后得到了最優(yōu)化網(wǎng)絡(luò)結(jié)構(gòu)，大幅提高了其計(jì)算精度，明顯優(yōu)于Rice等[27]采用5個(gè)量子化學(xué)模型所預(yù)測(cè)的結(jié)果。此外，王睿等[40]選用原子型電性拓?fù)錉顟B(tài)指數(shù)，表征了20種均三硝基苯類含能化合物。他們采用MLR方法進(jìn)行擬合，所創(chuàng)建的4參數(shù)線性模型預(yù)測(cè)效果較好。隨后，他們?cè)诖嘶A(chǔ)上采用原子電性拓?fù)錉顟B(tài)指數(shù)和基團(tuán)電性拓?fù)渲笖?shù)共同來表征包括硝基芳香化合物、硝酸酯和硝胺在內(nèi)的41種含能硝基化合物的分子結(jié)構(gòu)[41]。并采用逐步回歸MLR法成功建立了5變量線性預(yù)測(cè)模型。初步研究表明，電性拓?fù)錉顟B(tài)指數(shù)不僅可以反映硝基含能化合物的拓?fù)浣Y(jié)構(gòu)，還包含了其分子中電子狀態(tài)。為了擴(kuò)大樣本數(shù)，他們還采用156個(gè)硝基非雜環(huán)含能化合物進(jìn)行訓(xùn)練標(biāo)定，形成了基于MLR、PLS和BP神經(jīng)網(wǎng)絡(luò)等3種建模方法的預(yù)測(cè)模型。結(jié)果表明，非線性的BP神經(jīng)網(wǎng)絡(luò)方法構(gòu)建的預(yù)測(cè)模型無論在穩(wěn)定性、內(nèi)部及外部預(yù)測(cè)能力以及泛化性能方面都優(yōu)于線性方法(MLR與PLS)。近來，Morrill等[42]則利用該軟件在AM1半經(jīng)驗(yàn)水平計(jì)算了227個(gè)化合物的結(jié)構(gòu)參數(shù)，然后結(jié)合軟件集成的最優(yōu)多元線性回歸(BMLR)算法，從大量算符中篩選出8個(gè)并建立線性模型，取得了較好的結(jié)果。肖鶴鳴等[43]采用HMO、CNDO/2、MINDO/3和MNDO等分子軌道算法，對(duì)苯、甲苯、苯胺及苯酚等4類分子的硝基衍生物進(jìn)行了系統(tǒng)研究。結(jié)果表明，對(duì)結(jié)構(gòu)相似物，其分子中最弱鍵的鍵級(jí)(如π鍵、Mullikan鍵或Wiberg鍵)或雙原子作用能與其撞擊感度或熱安定性之間存在漸變關(guān)系。因而，根據(jù)炸藥熱分解和起爆機(jī)理，他們提議以基態(tài)分子最弱鍵的鍵級(jí)或該鍵所連接的兩原子之間的相互作用能作為判據(jù)來判別同系物炸藥的熱安定性或撞擊感度的相對(duì)大小。判斷方法有兩種:(1)由II級(jí)鍵估算化學(xué)鍵的離解能來判別；(2)根據(jù)鍵級(jí)和雙原子作用能的線性相關(guān)性來判斷。這些研究對(duì)炸藥的撞擊感度影響因素有更深層次的認(rèn)識(shí)，對(duì)炸藥其他爆炸性能的預(yù)測(cè)也有重大的指導(dǎo)意義。

Kim等[44]根據(jù)范德華分子表面靜電勢(shì)(MSEP)的定量結(jié)構(gòu)-性質(zhì)關(guān)系(QSPR)更精確地預(yù)測(cè)含能材料的撞擊感度。他們從MSEP衍生的各種三維描述，利用總和為正MSEP的變化，并結(jié)合其他3個(gè)參數(shù)，確立了新的QSPR方程。在此基礎(chǔ)上建立了如下6種不同精度的模型：

(1)

(2)

模型3:h50=a1+a2exp(a3υ)

(3)

模型4:h50=a1+a2exp[-a3(Q-a4)]

(4)

模型5:h50=a1exp[a2υ-a3(Q-a4)]

(5)

(6)

(7)

(8)

(9)

(10)

從以上表達(dá)式可以看出，含能化合物分子表面電荷分布與其撞擊感度大體呈指數(shù)關(guān)系，而大多數(shù)QSPR研究所得的結(jié)論是簡(jiǎn)單的線性關(guān)系。綜上所述，量子化學(xué)方法可提供精確的結(jié)構(gòu)數(shù)據(jù)，但是需要耗費(fèi)大量的CPU時(shí)間，對(duì)計(jì)算機(jī)硬件要求較高。而QSPR方法則可以系統(tǒng)全面地描述含能材料分子結(jié)構(gòu)參數(shù)與其撞擊感度之間的內(nèi)在聯(lián)系，并建立相應(yīng)的預(yù)測(cè)模型。但是，一般采用的描述參量精度都集中在經(jīng)驗(yàn)、半經(jīng)驗(yàn)水平，精確度稍差。同時(shí)，QSPR方法所需要的實(shí)測(cè)感度訓(xùn)練數(shù)據(jù)源差別較大、可靠性不能得到保證，給研究帶來了一定的不確定性。在確定新含能化合物的撞擊感度后，可以進(jìn)一步確定其靜電火花感度，因?yàn)楦鶕?jù)熱點(diǎn)起爆理論，這兩者本質(zhì)上存在一定的相關(guān)性[45]。盡管如此，含能材料的靜電火花感度產(chǎn)生機(jī)制還有待于進(jìn)一步驗(yàn)證。

2 含能材料燃燒性能預(yù)估

含能材料的燃燒性能預(yù)估主要采用綜合性預(yù)估軟件ALE3D[46]。該軟件是利用任意拉格朗日-歐拉(ALE)技術(shù)進(jìn)行多物理場(chǎng)數(shù)值模擬計(jì)算的重要工具。它通過采用混合有限元和有限體積非結(jié)構(gòu)化網(wǎng)格來解決二維或三維流體的彈塑性響應(yīng)問題[47]。相比局限于拉格朗日或歐拉(平流)的工具，ALE和網(wǎng)絡(luò)松弛方法的引入拓寬了其應(yīng)用范圍。同時(shí)，在多物理場(chǎng)和復(fù)雜幾何形狀下能保持高精度模擬和高效率運(yùn)算。除了用作流體力學(xué)和結(jié)構(gòu)化學(xué)計(jì)算，ALE3D還可通過用戶分工合作方式整合各種嵌套軟件，因而具備多物理場(chǎng)模擬的能力。ALE3D附加的功能包括熱傳導(dǎo)、化學(xué)動(dòng)力學(xué)和擴(kuò)散傳質(zhì)、不可壓縮流體、廣泛的材料模型、化學(xué)反應(yīng)模擬、多相流、以及磁流體力學(xué)(隱式)、短(顯式)時(shí)間尺度方面的應(yīng)用[48-49]。可模擬爆轟、爆燃、對(duì)流燃燒，主要針對(duì)高能炸藥和推進(jìn)劑的熱、機(jī)械刺激響應(yīng)[50-51]。它用有限元法增強(qiáng)動(dòng)量守恒，使很多實(shí)際問題可由“拉格朗日+重映射”法來解決，還允許用戶在全拉格朗日法或全歐拉法模式下運(yùn)行。它可以通過任意拉格朗日-歐拉計(jì)算，模擬含能材料斷裂和戰(zhàn)斗部破片過程，從而準(zhǔn)確模擬殼體材料的表面形貌及其變化。

對(duì)于傳熱過程，ALE3D主要用于鋁合金制造行業(yè)。它通過繪制溫度曲線來描述從鋁合金板坯傳熱到滾筒過程中的溫度分布，這可以分析鋁板內(nèi)應(yīng)力松弛率[52]。除了集成材料復(fù)雜損傷演化模型，ALE3D代碼還可模擬表面摩擦滑動(dòng)和接觸熱阻抗。同時(shí)，可以模擬炸藥爆炸沖擊波和鋼筋混凝土柱接觸時(shí)的復(fù)雜相互作用，ALE框架內(nèi)允許完全耦合的流固耦合法。其中磁流體動(dòng)力學(xué)(MHD)模塊能解瞬態(tài)磁場(chǎng)對(duì)流擴(kuò)散方程、磁力與流體力學(xué)耦合方程、焦耳熱與傳熱耦合方程等。它采用單胞為中心的戈東諾夫型有限體積法來模擬多相或多個(gè)分散物質(zhì)，每個(gè)階段都擁有其特定的速度和狀態(tài)參數(shù)，其中每個(gè)階段都被視為連續(xù)流體[53]。

對(duì)于推進(jìn)劑燃燒性能(主要是能量性能和燃速)的計(jì)算還有幾種軟件：如美國(guó)的NASA-CEA、葡萄牙的THOR、加拿大的CERV、德國(guó)的ICT熱力學(xué)代碼、俄羅斯的REAL等，它們可對(duì)雙組分及以上的推進(jìn)劑進(jìn)行循環(huán)計(jì)算，采用標(biāo)準(zhǔn)試樣驗(yàn)證其精確度[54-55]。這些代碼都是以熱力學(xué)計(jì)算為基礎(chǔ)，在既定燃燒產(chǎn)物的前提下進(jìn)行的計(jì)算，其精度取決于材料的組分[56-57]，對(duì)含金屬燃料、尤其是含有機(jī)金屬化合物的配方計(jì)算誤差較大。

3 含能材料爆轟性能預(yù)估

含能材料的爆轟性能與爆熱密切相關(guān)，而爆轟過程的反應(yīng)熱可由爆轟反應(yīng)產(chǎn)物的生成熱(HOFs)直接求得。通過選用適當(dāng)?shù)姆椒ǎ瑴?zhǔn)確計(jì)算含能材料分子的生成熱是量子化學(xué)的主要優(yōu)勢(shì)。最早使用半經(jīng)驗(yàn)分子軌道方法(如MNDO、AM1、PM3等)可以快捷地計(jì)算生成熱[58]。但這些半經(jīng)驗(yàn)方法主要依據(jù)一些代表性小分子和烴類的熱力學(xué)和光譜數(shù)據(jù)進(jìn)行參數(shù)化，對(duì)含多種取代基或特殊結(jié)構(gòu)的高能量密度材料，生成熱的計(jì)算往往存在較大誤差。一些采用改進(jìn)的半經(jīng)驗(yàn)方法如成對(duì)距離定向高斯修正法(PDDG)可減少誤差[59]。對(duì)622個(gè)含CHNO化合物的生成熱計(jì)算表明，理論值與實(shí)驗(yàn)值的絕對(duì)誤差由PM3法的21.3 kJ/mol下降至PDDG/PM3方法的13.4 kJ/mol[60]。

隨著計(jì)算技術(shù)的發(fā)展，能夠?qū)崿F(xiàn)對(duì)絕大多數(shù)高能量密度材料的第一性原理計(jì)算。因此半經(jīng)驗(yàn)方法已逐漸被第一性原理方法所代替。但第一性原理只能求分子的總能量，而無法直接計(jì)算生成熱。這就需要設(shè)計(jì)等鍵反應(yīng)，利用參考物的實(shí)驗(yàn)生成熱，借助Hess定律，求得目標(biāo)分子的生成熱。在等鍵反應(yīng)中，反應(yīng)物和產(chǎn)物的電子環(huán)境相近，電子相關(guān)能造成的誤差可以大部分抵消，可大幅降低計(jì)算的生成熱誤差。在設(shè)計(jì)等鍵反應(yīng)時(shí)，通常根據(jù)鍵分離規(guī)則(BSR)把分子分解成一系列與所求物質(zhì)具有相同化學(xué)鍵類型的小分子(已知生成熱的參考物)。為進(jìn)一步減少誤差，應(yīng)盡量保持母體骨架或原有分子的化學(xué)鍵。由下列公式：

ΔH298=∑ΔfHP-∑ΔfHR

(11)

ΔH298=ΔE298+Δ(PV)=ΔE0+ΔZPE+

ΔHT+ΔnRT

(12)

式中：∑ΔfHP和∑ΔfHR分別為298K下等鍵反應(yīng)產(chǎn)物和反應(yīng)物的生成熱之和；ΔE0和ΔZPE分別為0 K時(shí)產(chǎn)物與反應(yīng)物的總能量和零點(diǎn)能(ZPE)之差；ΔHT為從0 K到298 K的焓值溫度校正項(xiàng)；Δ(PV)在理想狀態(tài)條件下為ΔnRT；對(duì)于等鍵反應(yīng)，Δn=0，故Δ(PV)=0。

聯(lián)立上述兩式，可由參考物的生成熱求得目標(biāo)分子的氣相生成熱。第一性原理計(jì)算方法很多，但對(duì)高能量密度材料計(jì)算方法很有限。雖然閉合層的限制性Hartree-Fock(HF)方法可得到準(zhǔn)確的分子幾何構(gòu)型，但所得分子能量與實(shí)驗(yàn)誤差通常在200kJ/mol以內(nèi)，且個(gè)別誤差高達(dá)700kJ/mol。因HF方法忽略了電子相關(guān)效應(yīng)，導(dǎo)致能量值出現(xiàn)系統(tǒng)正誤差。為此，可用微擾法(MP2和MP4等)、多組態(tài)法(如CISD、CCSD和CASSCF)及DFT法(如B3LYP、B3PW91)校正電子相關(guān)效應(yīng)。微擾法和DFT法能量誤差通常在40kJ/mol以內(nèi)，而多組態(tài)法誤差一般小于5～8kJ/mol[61]。然而，多組態(tài)法對(duì)中等體系的計(jì)算量也非常大，因此通常用DFT方法處理電子相關(guān)效應(yīng)。雖然DFT方法的能量絕對(duì)誤差較大(40kJ/mol)，但是由于生成熱的計(jì)算是通過設(shè)計(jì)等鍵反應(yīng)來實(shí)現(xiàn)的，反應(yīng)物與產(chǎn)物的能量誤差大部分相互抵消，即生成熱的真實(shí)誤差通常比40kJ/mol小得多。關(guān)于DFT具體泛函的選用，通常要以同類型化合物的準(zhǔn)確實(shí)驗(yàn)生成熱為基準(zhǔn)，比較并檢驗(yàn)各種泛函計(jì)算結(jié)果對(duì)該類化合物的準(zhǔn)確性，以確定最佳方案。

量子化學(xué)計(jì)算得到的是氣相生成熱，結(jié)合遺傳算法(Genetic Algorithm, GA)、靜電勢(shì)法和神經(jīng)網(wǎng)絡(luò)法(Neutral Network, NN, 例如Back-propagation Neural Networks, BPNN)等手段，計(jì)算固體升華熱后即可求得固相生成熱。舒遠(yuǎn)杰等[62]采用上述方法對(duì)72種籠狀和橋環(huán)類分子的固態(tài)標(biāo)準(zhǔn)生成熱進(jìn)行計(jì)算，結(jié)果表明，其標(biāo)準(zhǔn)誤差約為20.9kJ/mol。相對(duì)于第一性原理，靜電勢(shì)法因其物理意義明確、結(jié)果誤差小，在高能量密度材料升華熱計(jì)算中得到了廣泛應(yīng)用[63-64]。對(duì)于有機(jī)分子所組成的固體而言，其升華熱取決于分子間相互作用能，相互作用能越大則升華熱越大。因此，升華熱與分子表面靜電勢(shì)有關(guān)。由靜電勢(shì)求得升華熱的原理及相關(guān)公式可參見文獻(xiàn)[65]。

目前，基于1384個(gè)化合物(含172個(gè)化學(xué)基團(tuán))的升華熱，建立了3層前饋神經(jīng)網(wǎng)絡(luò)(feed forward neural network)計(jì)算模型[66]。采用該模型計(jì)算的升華熱值與實(shí)驗(yàn)結(jié)果的相關(guān)系數(shù)平方、平均誤差和均方根誤差分別為0.9854、3.54%和4.21kJ/mol[67]。后者明顯小于靜電勢(shì)法所得誤差(11.7J/mol)，具有一定的普適性和精確性。但由于樣本極少涉及高能量密度材料分子，且高能量密度材料的特殊性(強(qiáng)分子內(nèi)基團(tuán)的相互作用、較大環(huán)張力等)，該方法是否能準(zhǔn)確預(yù)測(cè)其升華熱還有待檢驗(yàn)。在獲得生成熱后，如果已知高能量密度材料分子的密度，即可由Kamlet經(jīng)驗(yàn)公式預(yù)測(cè)其爆速和爆壓[68]。此外，基于化合物優(yōu)化構(gòu)型，用Monte-Carlo方法求得分子周圍0.001e/bohr3等電子密度面所包圍的體積，即可求得摩爾體積(V)。因該法求得的體積值波動(dòng)范圍較大，通常要取重復(fù)100次以上的平均值。由分子質(zhì)量與平均摩爾體積之比可計(jì)算該分子的理論密度(ρ)。研究表明，在B3LYP/6-31G**水平下計(jì)算的理論密度與實(shí)測(cè)值吻合度較好[69]。將所求得的生成熱及密度值代入Kamlet公式，即可得到單質(zhì)含能材料的爆速和爆壓。第一性原理計(jì)算結(jié)合等鍵反應(yīng)不僅可求得生成熱，還可計(jì)算基團(tuán)相互作用能和環(huán)張力，從分子水平闡明結(jié)構(gòu)與性能的關(guān)系。當(dāng)然，計(jì)算目的不同，所設(shè)計(jì)的等鍵反應(yīng)也不同。

為了有效地預(yù)估新含能材料及其配方的性能，科研人員一直在探索精確計(jì)算爆轟參數(shù)的算法和軟件。除了最早Rothstein(1979和1981年)和Kamlet(1968年)分別創(chuàng)立的簡(jiǎn)便方法之外，目前已經(jīng)發(fā)展建立了一系列更精確的狀態(tài)方程。通常，熱化學(xué)計(jì)算比量子化學(xué)計(jì)算更加便捷、成本較低、且準(zhǔn)確性和可靠性也越來越高。比較常用的熱力學(xué)代碼有BKW-Fortran、ARPEGE、Ruby、TIGER[70]、CHEETAH、EXPLO-5、MWEQ、BARUT-X和ZMWNI[71]。目前，美國(guó)的TIGER和CHEETAH及克羅地亞的EXPLO-5是被廣泛使用的3種預(yù)測(cè)推進(jìn)劑和炸藥性能的可靠熱化學(xué)代碼。它們都可以由研究人員來任意增添組分或改變配方，以優(yōu)化所設(shè)計(jì)產(chǎn)品的性能。最新版本的CHEETAH 7.0是通過美國(guó)能源部和國(guó)防部與澳大利亞防御局(DSTO)武器發(fā)展與防御數(shù)據(jù)相互交換協(xié)議下共同開發(fā)的[72]。然而，CHEETAH不僅用于高能炸藥的爆轟參數(shù)計(jì)算，還可以預(yù)測(cè)許多復(fù)雜材料的熱力學(xué)行為，包括塑料、有機(jī)液體混合物、固-液混合物等在沖擊波、激光或地球內(nèi)部等極端條件下的性能。CHEETAH通過不斷拓展其功能，最新版本可以模擬化學(xué)和物理動(dòng)力學(xué)，與先進(jìn)的流體力學(xué)聯(lián)合使用可解決復(fù)雜高能材料體系很多應(yīng)用基礎(chǔ)問題。TIGER和CHEETAH都是源于美國(guó)的Ruby和BKW。目前，美國(guó)已將這兩個(gè)軟件內(nèi)部合并為CTH-TIGGER。新代碼包含了H2O的最新狀態(tài)方程，并新增了200多種離子化合物相關(guān)數(shù)據(jù)。

EXPLO-5是基于化學(xué)方程式、生成熱和密度，預(yù)測(cè)高能炸藥、推進(jìn)劑和煙火藥爆炸參數(shù)的一款熱化學(xué)計(jì)算程序。EXPLO-5在含能材料的合成、配方優(yōu)化和數(shù)學(xué)建模中是一個(gè)非常重要的工具。EXPLO-5運(yùn)用自由能最小化方法在指定的溫度和壓力下計(jì)算平衡組成和熱力性質(zhì)，數(shù)據(jù)結(jié)合Chapman-Jouguet爆炸理論，能夠計(jì)算爆炸參數(shù)，如爆速、爆壓、爆熱和爆溫等。從平衡組成到熱力學(xué)狀態(tài)參數(shù)和等熵膨脹，軟件通過內(nèi)置擬合程序計(jì)算Jones-Wilkins-Lee (JWL)狀態(tài)方程中的系數(shù)。通過解樣品的熱力學(xué)方程和恒壓燃燒條件下的守恒方程來預(yù)測(cè)固體推進(jìn)劑的燃燒性能，如比沖、推力系數(shù)和定容推力等。該程序采用氣態(tài)爆轟產(chǎn)物的Becker-Kistiakowsky- Wilson(BKW)狀態(tài)方程和Jacobs-Cowperthwaite- Zwisler (JCZ3)狀態(tài)方程、理想氣體方程、維也里狀態(tài)方程，以及Murnaghan狀態(tài)方程等。EXPLO-5的現(xiàn)行版本由于應(yīng)用JCZ3 EOS方程，爆炸參數(shù)的預(yù)測(cè)準(zhǔn)確性得到提高。EXPLO-5的數(shù)據(jù)庫目前包含260種反應(yīng)物，超過330種產(chǎn)物(包括同一產(chǎn)物的不同相態(tài))，并包含了32種元素：C, H, N, O, Al, Cl, Si, F, B, Ba, Ca, Na, P, Li, K, S, Mg, Mn, Zr, Mo, Cu, Fe, Ni, Pb, Sb, Hg, Be, Ti, I, Xe, U,和W。

除EXPLO-5外，波蘭華沙技術(shù)大學(xué)最近也開發(fā)了一款名為ZMWNI的熱力學(xué)軟件。它可以計(jì)算含能材料的燃燒、爆炸和熱力學(xué)參數(shù)，并確定爆轟產(chǎn)物的JWL等熵膨脹曲線[73]和爆轟能[74]等。此外，ZMWNI代碼能夠確定不同溫度下，配方在非平衡狀態(tài)的爆轟參數(shù)。程序基于最小化學(xué)勢(shì)法計(jì)算平衡或非平衡組分的反應(yīng)性體系。最終數(shù)據(jù)采集是通過求解線性方程組和最陡降法獲得的[75]。氣體的物理性質(zhì)則由BKW方程狀態(tài)描述。上述的TIGER和CHEETAH都采用凝聚相組分OLD狀態(tài)方程，用于計(jì)算給定條件下HEDM的燃燒或爆炸的平衡狀態(tài)和定容爆轟參數(shù)。ZMWNI熱力學(xué)軟件則不需要給定條件，可以計(jì)算非理想狀態(tài)下的任何參數(shù)?？傊?，在熱力學(xué)計(jì)算領(lǐng)域，我國(guó)的發(fā)展水平相對(duì)落后，這在一定程度上制約我國(guó)新型含能材料的理論設(shè)計(jì)、開發(fā)和應(yīng)用。

4 含能材料反應(yīng)活性的分子動(dòng)力學(xué)模擬

原子級(jí)別的模擬方法有：分子動(dòng)力學(xué)(MD)和蒙特卡諾(MC)模擬。它們考慮單個(gè)原子或官能團(tuán)，可以預(yù)測(cè)含能材料晶體、晶體與黏合劑之間界面的動(dòng)態(tài)性能。分子動(dòng)力學(xué)模擬中原子間作用力可以通過經(jīng)驗(yàn)力場(chǎng)參數(shù)來計(jì)算，而力場(chǎng)又是基于電子結(jié)構(gòu)和體系瞬態(tài)電子的波函數(shù)。其中分子動(dòng)力學(xué)模擬能量場(chǎng)中的擴(kuò)散比較有效。目前主要商業(yè)軟件包括美國(guó)的Material Studio和荷蘭的ADF軟件(后面提到的反應(yīng)力場(chǎng)ReaxFF模塊也已經(jīng)嵌入該軟件)。分子間相互作用遵循一定的近似物理規(guī)律，在原子水平上有一定的復(fù)雜性。因此，“力場(chǎng)”被用來模擬結(jié)構(gòu)、振動(dòng)、構(gòu)象和熱物理特性。含能材料研究中常用的力場(chǎng)是“簡(jiǎn)明相優(yōu)化分子勢(shì)原子力場(chǎng)”(COMPASS)。 COMPASS力場(chǎng)可準(zhǔn)確地模擬孤立分子和固相材料分子的性能。分子模型COMPASS力場(chǎng)可以克服現(xiàn)實(shí)實(shí)驗(yàn)中遇到的困難，已成功用于模擬增塑劑擴(kuò)散性以及它們與HTPB推進(jìn)劑黏結(jié)體系的相容性[76]。COMPASS力場(chǎng)還被用來模擬HTPB對(duì)Al、氧化鋁吸附能和它們之間的界面力學(xué)性能(如彈性系數(shù)、模量和泊松比等)。通過分析吸附能量和相關(guān)函數(shù)，可以模擬界面相互作用的性質(zhì)。在涉及建模和仿真熱力學(xué)領(lǐng)域，國(guó)內(nèi)已開展了一定的工作。蕭忠良等[77]完成了一些發(fā)射藥氣體體積流量、氧平衡和能量釋放量的建模。

含能材料分子動(dòng)力學(xué)研究的熱點(diǎn)還包括采用ReaxFF來模擬含能材料的熱分解和燃燒過程[78]。ReaxFF基于從頭量子力學(xué)(QM)計(jì)算，可準(zhǔn)確描述烴和各種CHNO體系的反應(yīng)性，如RDX同系物的分解途徑及后續(xù)氣相化學(xué)反應(yīng)[79-80]。在模擬高能材料分解的化學(xué)途徑方面，“HE”力場(chǎng)是專門為硝銨化合物構(gòu)建的[81-82]。在“HE”力場(chǎng)基礎(chǔ)上發(fā)展了另一個(gè)普適于含能化合物的力場(chǎng)“CHONSSi-lg”，該力場(chǎng)考慮了倫敦色散力的影響[83]，但它不包含氟原子的數(shù)據(jù)。為了更好地比較含氟配方(如含氟聚物)，可采用“TiOCHNCl”力場(chǎng)[84]。除ReaxFF可以模擬計(jì)算含能材料化學(xué)反應(yīng)機(jī)理之外，還可采用CFD軟件進(jìn)行建模計(jì)算。如北京化工大學(xué)開發(fā)了基于渦級(jí)聯(lián)和分形理論的現(xiàn)象學(xué)概念來模擬湍流化學(xué)反應(yīng)的分形模型(FM)，采用湍流燃燒化學(xué)反應(yīng)和渦耗散概念(EDC)模型，可以模擬含能材料燃燒反應(yīng)機(jī)理[85]。王國(guó)青等[86]使用DFT和高斯03軟件，研究了TNT與硫酸根的反應(yīng)機(jī)理。在他們的研究中，所有反應(yīng)物分子、過渡態(tài)結(jié)構(gòu)和產(chǎn)物均在B3LYP/6-31G(d)水平下進(jìn)行了幾何優(yōu)化。除了研究反應(yīng)機(jī)理，數(shù)值建模和仿真也被用于反應(yīng)動(dòng)力學(xué)研究。如機(jī)械損傷對(duì)PBX起爆和爆轟行為的影響。該反應(yīng)模型可以分析含缺陷的炸藥爆轟機(jī)制，其孔隙率和粒徑隨沖擊載荷而改變[87]。建模和仿真還可以用于模擬推進(jìn)劑的燃燒。如封鋒等[88]開發(fā)的SPRS軟件可以計(jì)算推進(jìn)劑的化學(xué)成分及其在指定壓力下的燃速和壓力指數(shù)。固體推進(jìn)劑的組分可以根據(jù)所需的燃燒速率和壓力指數(shù)來確定。此數(shù)值模擬縮短固體推進(jìn)劑的研發(fā)周期，進(jìn)而節(jié)約成本。發(fā)射藥方面，張江波等[89]開發(fā)了基于經(jīng)典內(nèi)彈道理論產(chǎn)生的多層推進(jìn)劑裝藥內(nèi)彈道數(shù)學(xué)模型模擬程序。發(fā)射藥和推進(jìn)劑數(shù)值模擬研究包括以下幾個(gè)方面：(1)由于對(duì)發(fā)射藥的發(fā)射負(fù)荷引起的內(nèi)孔變形；(2)變?nèi)妓侔l(fā)射藥的逐層可燃性；(3)環(huán)境溫度、組分、粒徑和硼粉氧化膜厚度對(duì)點(diǎn)火性能的影響；(4)硼粉團(tuán)聚對(duì)富燃料固體推進(jìn)劑能量釋放速率和燃速性能的影響。此外，蕭忠良等[90]發(fā)展了推進(jìn)劑定容燃燒行為的數(shù)學(xué)模型，并在此基礎(chǔ)上研究了腹板厚度和非同步點(diǎn)火對(duì)推進(jìn)劑燃速系數(shù)的影響。數(shù)值建模和仿真還可用于某些特殊研究領(lǐng)域。如使用熱機(jī)械耦合有限元法研究在固體推進(jìn)劑整形加工過程的危險(xiǎn)性。采用基于黏彈性積分本構(gòu)關(guān)系的方法，確定了由整形過程中壓力和摩擦產(chǎn)生的熱量，并且對(duì)危險(xiǎn)點(diǎn)進(jìn)行計(jì)算。數(shù)值模擬的另一種應(yīng)用是基于固化時(shí)間理論的修正蠕變模型，該模型用以模擬塑料黏結(jié)炸藥(PBX)組分短時(shí)間內(nèi)的蠕變行為。該模型還被用來模擬HMX-PBX炸藥在受真空吸力過程中的緩慢蠕變行為。數(shù)值模擬和理論分析表明，改進(jìn)的固化時(shí)間蠕變模型適用于模擬PBX瞬時(shí)蠕變行為。

5 結(jié)束語

含能材料的發(fā)展日臻完善，主要研究領(lǐng)域包括含能化合物的分子設(shè)計(jì)、含能材料納米化、共晶化與表面改性、含能材料分解點(diǎn)火燃燒與爆轟、極限條件下含能材料的響應(yīng)、含能材料配方設(shè)計(jì)及工程化應(yīng)用、含能材料安全性能與綠色利用，及含能材料的量化計(jì)算與多尺度分子模擬等。對(duì)于含能材料的物理化學(xué)性能預(yù)估，可以采用量子化學(xué)、分子動(dòng)力學(xué)或滿足半經(jīng)驗(yàn)QSPR模型的方法來預(yù)估含能材料的感度、燃燒性能、爆轟性能、反應(yīng)活性、固化機(jī)制與力學(xué)性能等。目前存在的主要發(fā)展瓶頸有：

(1)沒有建立統(tǒng)一標(biāo)準(zhǔn)的含能材料性能實(shí)驗(yàn)數(shù)據(jù)庫，無法滿足匹配大量模擬和模型訓(xùn)練的需求。且諸多數(shù)據(jù)的測(cè)試標(biāo)準(zhǔn)不統(tǒng)一，人為因素影響了結(jié)果的可靠性；

(2)建模和參數(shù)有待進(jìn)一步系統(tǒng)化和標(biāo)準(zhǔn)化，國(guó)內(nèi)未能創(chuàng)立商業(yè)化計(jì)算軟件，且對(duì)國(guó)外相關(guān)商業(yè)化軟件的二次開發(fā)較少。有待于進(jìn)一步加大模擬仿真方面的投入力度，確保建立更加精確、可靠，普適性強(qiáng)的預(yù)測(cè)模型；

(3)目前商業(yè)軟件僅局限于預(yù)測(cè)含能材料的爆轟性能和燃燒性能，對(duì)于安全性能的理論計(jì)算還處于探索階段。有待以建模仿真和量子化學(xué)計(jì)算為基礎(chǔ)創(chuàng)立一個(gè)能評(píng)價(jià)含能材料理化與安全性能的綜合體系。

[1] Short J M,Kavetsky R A, Davinder K.Energetics Science and Technology in Central Europe [M].Maryland:CALCEEPSC Press,2013.

[2] 王睿, 蔣軍成, 潘勇, 硝基含能材料撞擊感度的預(yù)測(cè)研究進(jìn)展[J], 工業(yè)安全與環(huán)保, 2010, 36(7):19-22.

WANG Rui, JIANG Jun-cheng, PAN Yang. Research on the prediction of impact sensitivity of nitro energetic materials [J]. Industrial Safety and Environment Protection, 2010, 36(7):19-22.

[3] 王澤山.含能材料概論 [M]. 哈爾濱: 哈爾濱工業(yè)大學(xué)出版社, 2006.

[4] 董海山. 評(píng)介《四唑化學(xué)的現(xiàn)代理論》[J]. 含能材料, 2002(2):95-96.

[5] Zeman S, Jungová M. Sensitivity and performance of energetic materials[J]. Propellants, Explosives, Pyrotechnics,2016,41:426-451.

[6] 房偉, 王建華, 劉玉存,等. 基于分子基團(tuán)預(yù)測(cè)硝基含能材料撞擊感度[J].火工品, 2014(5):34-37.

FANG Wei, WANG Jian-hua, LIU Yu-cun, et al.Prediction on impact sensitivity of nitro-energetic materials by molecular groups [J]. Initiators & Pyrotechnics,2014(5):34-37.

[7] 劉歡, 姜峰, 于國(guó)強(qiáng),等.遺傳-神經(jīng)網(wǎng)絡(luò)方法在炸藥撞擊感度預(yù)測(cè)中的應(yīng)用研究[J], 火工品, 2010(6):42-45.

LIU Huan, JIANG Feng, YU Guo-qiang,et al. Application research of genetic-neural network method in the prediction of explosives′ impact sensitivity [J].Initiators & Pyrotechnics, 2010(6):42-45.

[8] 王睿, 蔣軍成, 潘勇,等. 均三硝基苯類化合物撞擊感度與電性拓?fù)渲笖?shù)的QSPR研究[J].含能材料, 2008, 16(1):90-93.

WANG Rui, JIANG Jun-cheng, PAN Yong, et al. QSPR study of correlation between impact sensitivity of m-nitroaromatics and electrotopological state indices [J]. Chinese Journal of Energetic Materials, 2008, 16(1):90-93.

[9] 肖鶴鳴.四唑化學(xué)的現(xiàn)代理論 [M].北京：科學(xué)出版社, 2000.

[10] 肖鶴鳴, 朱衛(wèi)華, 肖繼軍,等.含能材料感度判別理論研究——從分子、晶體到復(fù)合材料[J].含能材料, 2012, 20(5):514-527.

XIAO He-ming, ZHU Wei-hua, XIAO Ji-jun. et al. Theoretical studies on sensitivity criterion of energetic materials-from molecules, crystals, to composite materials[J]. Chinese Journal of Energetic Materials, 2012, 20(5):514-527.

[11] Nefati H,Cense J M, Legendre J J. Prediction of the impact sensitivity by neural networks [J].Journal of Chemistry Information Computer Science, 1996, 36: 804.

[12] Smirnov A, Voronko O, Korsunsky B, et al. Impact and friction sensitivity of energetic materials:methodical evaluation of technological safety features [J]. Chinese Journal of Explosives and Propellants, 2015, 38 (3): 1-8.

[13] 金韶華, 王偉, 松全才. 含能材料機(jī)械撞擊感度判據(jù)的認(rèn)識(shí)和發(fā)展 [J]. 爆破器材, 2006, 35(6):11-14.

JIN Shao-hua, WANG Wei, SONG Quan-cai, Understanding and development of criteria of impact sensitivity of energetic materials [J]. Explosive Materials, 2006, 35(6):11-14.

[14] 李金山, 曹剛, 肖鶴鳴, 等. 多硝基芳香化合物撞擊感度的量子化學(xué)研究 [J]. 火炸藥學(xué)報(bào), 1997,20(2):57-58.

LI Jin-shan, CAO Gang, XIAO He-ming,et al. Quantum chemical study on the impact sensitivity of polynitroaromatics[J].Chinese Journal of Explosives & Propellants(Huozhayao Xuebao), 1997,20(2):57-58.

[15] Bowden F, Yoffe Y. Initiation and Growth of explosion in Liquids and Solids [M].London:Cambridge University Press, 1952.

[16] Delpuech A, Cherville J, Michaud C. Molecular electronic structure and initiation of secondary explosives[C]∥7th Symposium (International) on Detonation. Annapolis:[s.n.],1981.

[17] Kamlet M J, Adolph H G. The relationship of impact sensitivity with structure of organic high explosices [J]. Propellants, Explosives, Pyrotechnics, 1979, 4: 30.

[18] Chen Z X, Xiao H M. Quantum chemistry derived criteria for impact sensitivity [J].Propellants, Explosives, Pyrotechnics, 2014, 39 (4):487-495.

[19] Murray J S, Lane P, Politzer P.Relationships between impact sensitivities and molecular surface electrostatic potentials of nitroaromatic and nitroheterocyclic molecules [J]. Molecular Physics, 1995, 85 (1): 1-8.

[20] Politzer P, Murray J S. Impact sensitivity and crystal lattice compressibility/free space [J].Journal of Molecular Modeling, 2014,20(5):2223.

[21] Keshavarz M H. Calculations of bond dissociation energies and dipole moments in energetic materials using density-functional methods [J]. Journal of Hazardous Material, 2007, 148: 658-662.

[22] Keshavarz M H, Motamedoshariati H,Pouretedal H R,et al. Prediction of shock sensitivity of explosives based on small-scale gap test [J]. Journal of Hazardous Material, 2007, 145: 109-112.

[23] Keshavarz M H, Shokrolahi A, Esmailpoor K,et al, Prediction of heats of sublimation of energetic compounds using their molecular structures [J].Energetic Materials, 2008, 177(1-3):648-659.

[24] Liu X, Su Z, Ji W,et al. Structure, physicochemical properties, and density functional theory calculation of high-energy-density materials constructed with intermolecular interaction: nitro group charge determines sensitivity [J]. Journal of Physical Chemistry C, 2014, 118 (41):23487-23498.

[25] Mullay J. A relationship between impact sensitivity and molecular electronegativity [J]. Propellants, Explosives, Pyrotechnics, 1987, 12 (2): 60-63.

[26] Ren F D, Cao D L, Shi W J,et al. A theoretical prediction of the relationships between the impact sensitivity and electrostatic potential in strained cyclic explosive and application to H-bonded complex of nitrocyclohydrocarbon [J].Journal of Molecular Modeling, 2016, 22(4).

[27] Rice B M, Hare J J, Phys J.A quantum mechanical investigation of the relation between impact sensitivity and the charge distribution in energetic molecules[J].Journal of Physical Chemistry Chem A,2002, 106, 1770-1783.

[28] Politzer P, Murray J S. Impact sensitivity and the maximum heat of detonation[J]. Journal of Molecular Modelling, 2015, 21 (10):262:11.

[29] Zohari N, Keshavarz M H, Seyedsadjadi S A A link between impact sensitivity of energetic compounds and their activation energies of thermal decomposition [J]. J Therm Anal Calorim, 117(1):423-432.

[30] Mathieu D, Alaime T. Predicting impact sensitivities of nitro compounds on the basis of a semi-empirical rate constant [J]. Journal of Physical Chemistry A, 2014, 118 (41):9720-9726.

[31] Tan B S, Huang M, Li J S,et al. A new sensitivity criterion of explosives: bonding & nonbonding coupling related molecular rigidity and flexibility [J]. Chinese Journal of Energetic Materials, 2016,24 (1):10-18.

[32] Keshavarz M H, Motamedoshariati H, Moghayadnia R,et al. Prediction of sensitivity of energetic compounds with a new computer code [J]. Propellants, Explosives, Pyrotechnics, 2014, 39(1): 95-101.

[33] Mullay J. Relationships between impact sensitivity and molecular electronic structure [J]. Propellants, Explosives, Pyrotechnics, 1987, 12 (4): 121-124.

[34] Edwards J, Eybl C, Johnson B. Correlation between sensitivity and approximated heats of detonation of several nitroamines using quantum mechanical methods [J]. International Journal of Quantum Chemistry, 2004,100 (5):713-719.

[35] 曹霞, 向斌, 張朝陽. 炸藥分子和晶體結(jié)構(gòu)與其感度的關(guān)系[J]. 含能材料, 2012(5): 643-649.

CAO Xia, XIANG Bin, ZHANG Chao-yang. Review on relationships between the molcular and crystal structure of explosives and their sensitivities [J]. Chinese Journal of Energetic Materials, 2012(5): 643-649.

[36] Nefati H, Cense J M, Legendre J J. Prediction of the impact sensitivity by neural networks[J]. Journal of Chemical Information and Computer Sciences, 1996, 36 (4): 804-810.

[37] Cho S G., No K T , Goh E M, et al. Optimization of neural networks architecture for impact sensitivity of energetic molecules[J]. Bulletin of the Korean Chemical Society, 2005, 26 (3):399-408.

[38] Wang R, Jiang J, Pan Y. Prediction of impact sensitivity of nonheterocyclic nitroenergetic compounds using genetic algorithm and artificial neural network [J]. Journal of Energetic Materials, 2012,30 (2):135-155.

[39] Keshavarz M H, Jaafari M. Investigation of the various structure parameters for predicting impact sensitivity of energetic molecules via artificial neural network[J]. Propellants, Explosives, Pyrotechnics, 2006, 31 (3):216-225.

[40] 王睿, 蔣軍成, 潘勇, 等. 均三硝基苯類化合物撞擊感度與電性拓?fù)鋺B(tài)指數(shù)的QSPR的研究[J]. 含能材料, 2008, 16(1): 90-93.

WANG Rui, JIANG Jun-cheng,PAN Yong, et al. QSPR study of correlation between impact sensitivity ofm-nitroaromatics and electrotopological state indices[J].Chinese Journal of Energetic Materials,2008,16(1):90-93.

[41] 王睿, 蔣軍成, 潘勇, 等. 電性拓?fù)鋺B(tài)指數(shù)預(yù)測(cè)硝基類含能材料撞擊感度[J].固體火箭技術(shù), 2008, 31(6): 647-662.

WANG Rui, JIANG Jun-cheng,PAN Yong, et al. Predicton on impact sensitivity of nitro energetic materials by means of electrotopological state indices[J]. Journal of Solid Rocket Technology,2008,31(6):647-662.

[42] Jason A Morrill,Edward F C Byrd.Development of quantitative structure-property relationships for predictive modeling and design of energetic materials[J].Journal of Molecular Graphics and Modeling,2008,27(3):349-355.

[43] 肖鶴鳴, 王遵堯, 姚劍敏. 芳香族硝基炸藥感度和安定性的量子化學(xué)研究-I. 苯胺類硝基衍生物[J].化學(xué)學(xué)報(bào), 1985, 43 (1): 14-18.

XIAO He-ming,WANG Zun-yao,YAO Jian-min. Quantum chemical study on sensitivity and stability of aromatic nitro explosives i. nitro derivatives of aminobenzenes[J]. Acta Chimica Sinica,1985,43(1):14-18.

[44] Kim C K, Cho S G, Li J, et al. QSPR studies on impact sensitivities of high energy density molecules [J]. Bulletin Korean Chemical Society, 2011, 32(12):4341-46.

[45] Keshavarz M H, Keshavarz Z. Relation between electric spark sensitivity and impact sensitivity of nitroaromatic energetic compounds [J]. Zeitschrift fur Anorganische and Allgemeine Chemie, 2016, 642 (4):335-342.

[46] Cazamias J U, Schraml S J. Modeling of non-eroding penetration using ALE3D and Zapotec [C]∥ AIP Conference Proceedings. [S.l.]:Inverse Problems International Association, 2007, 955:1281-1284.

[47] Pierce T, Rodrigue G. A parallel two-sided contact algorithm in ALE3D [J]. Computer Methods in Applied Mechanics and Engineering, 2005,194 (27-29):3127-3146.

[48] McClelland M A, Maienschein J L, Howard W M, et al. ALE3D simulation of heating and violence in a fast cookoff experiment with LX-10 [C]∥ Proceedings of the 13th International Detonation Symposium.[S.l.]:Office of Naval Research,2006:606-616.

[49] Kuo I W, Bastea S, Fried L E. Reactive flow modeling of liquid explosives via ale3d/cheetah simulations [C]∥ Proceedings of 14th International Detonation Symposium. [S.l.]:Office of Naval Research, 2010:333-337.

[50] O′Daniel J L. Modeling contact detonations with ALE3D and PARADYN [J]. Computational Fluid and Solid Mechanics ,2003:537-539.

[51] Howard W M, McClelland M A, Knap J, et al. Ale3d simulations of gap closure and surface ignition for cookoff modeling [C]∥ Proceedings of the 13th International Detonation Symposium.[S.l.]:Office of Naval Research,2006: 1247-1255.

[52] Knap J, McClelland M A, Maienschein J L, et al. Measurement and ale3d simulation of violence in a deflagration experiment with lx-10 and aermet-100 alloy [C]∥Proceedings of the 13th International Detonation Symposium.[S.l.]:Office of Naval Research,2006:662-671.

[53] Stewart J J P. MOPAC: A semiempirical molecular orbital program [J]. Journal of Computer-Aided Molecular Design, 1990, 4(1): 1-103.

[54] Gicquel L Y M, Staffelbach G, Poinsot T. Large Eddy Simulations of gaseous flames in gas turbine combustion chambers[J]. Progress in Energy and Combustion Science,2012, 38 (6): 782-817.

[55] Beckstead M W, Liang Y, Pudduppakkam K V. Numerical simulation of single aluminum particle combustion [J]. Combustion, Explosion and Shock Waves, 2005, 41 (6): 622-638.

[56] Ermolin N E, Zarko V E. Simulation of combustion of cyclic nitramines [J].Fizika Goreniya i Vzryva, 1998, 34 (5): 3-22.

[57] Wang J, Wight C A. Use of kinetic models for solid state reactions in combustion simulations [J]. Theoretical Computational Chemistry, 2003,13: 351-371.

[58] Stewart J J P. MOPAC: A semiempirical molecular orbital program [J]. J Computer-Aided Mol Des, 1990, 4(1): 1-103.

[59] Bredow T, Jug K. Theory and range of modern semiempirical molecular orbital methods[J]. Theoretical Chemistry Accounts, 2005, 113(1):1-14.

[60] Repasky M P, Chandrasekhar J, Jorgensen W L. PDDG/PM3 and PDDG/MNDO: improved semiempirical methods [J]. Journal of Computational Chemistry, 2002, 23: 1601-1622.

[61] Dorsett H, White A. Overview of molecular modelling and Ab initio molecular orbital methods suitable for use with energetic materials[R].Salisbury: DSTO Aeronautical and Maritime Research Laboratory， 2000.

[62] 舒遠(yuǎn)杰. 含能材料的理論研究與數(shù)值模擬進(jìn)展 [M]. 綿陽：中國(guó)工程物理研究院科技年報(bào)(2009年版) , 2010： 115-117.

[63] Byrd E F C，Rice B M. Improved prediction of heats of formation of energetic materials using quantum mechanical calculations [J]. Journal of Physical Chemistry A,2006, 110(3): 1005-1013.

[64] Mohammad H K. Improved prediction of heats of sublimation of energetic compounds using their molecular structure[J]. Journal of Hazardous Materials, 2010, 177(1-3): 648-659.

[65] Politzer P, Murray J S. The fundamental nature and role of the electrostatic potential in atoms and molecules[J]. THEO CHEM ACCS, 2002, 108(3):134-142.

[66] Mathieu D. Simple alternative to neural networks for predicting sublimation enthalpies from fragment contributions[J]. Ind Eng Chem Res, 2012, 51(6):2814-2819.

[67] Farhad Gharagheizi, Mehdi Sattari, Behnam Tirandazi. Prediction of crystal lattice energy using enthalpy of sublimation: a group contribution-based model[J]. Ind Eng Chem Res, 2011, 50:2482-2486.

[68] Kamlet M J, Jacobs S J. Chemistry of detonations. I. A simple method for calculating detonation properties of C-H-N-O explosives[J]. Journal of Chemical Physics, 1968, 48:23-35.

[69] 邱玲. 氮雜環(huán)硝胺類高能量密度材料(HEDM)的分子設(shè)計(jì)[D]. 南京：南京理工大學(xué), 2007.

[70] Persson P A. TIGER WIN—a window PC code for computing explosive performance and thermodynamic properties [C]∥ Proceedings of 2000 High-tech Seminar, State-of-the Art Blasting Technology and Explosive Applications. Orlando:[s.n.],2000:541.

[71] Grys S, Trzciński W A.Characteristics of Energetic Materials[J].Central European Journal of Energetic Materials,2010,7(2):97-113.

[72] Lu J P. Evaluation of the Thermochemical Code - CHEETAH 2.0 for Modelling Explosives Performance, DSTO-TR-1199[R]. Australia:DSTO,2001:1-34.

[73] Lee E L, Horning H C, Kury J W. Adiabatic expansion of high explosive detonation products,UCRL-50422[R]. Livermore : Lawrence Livermore National Laboratory, 1968.

[74] Jacobs S J. Energy of detonation, NAVORD-4366[R].[S.l.]: Naval Ordnance Laboratory, 1956.

[75] White W B, Johnson S M, Danzig G B. Chemical equilibrium in complex mixtures [J].J Chem Phys,1958, 28: 751-755.

[76] Li H X, Qiang H F, Wu W M. Molecular simulation on plasticizer migration in the bond system of HTPB propellant [J]. Chinese Journal of Explosives & Propellants(Huozhayao Xuebao), 2008, 31(5):74-78.

[77] 代淑蘭,許厚謙,肖忠良,等. 帶制退器的膛口燃燒流場(chǎng)并行數(shù)值模擬 [J]. 彈道學(xué)報(bào), 2009,21(4):84-87.

DAI Shu-lan, XU Hou-qian, XIAO Zhong-liang. Numerical simulation of muzzle combustion flow field with brake by parallel computation [J]. Journal of Ballistics, 2009,21(4):84-87.

[78] Duin van A C T, Dasgupta S, Lorant F, et al. A reactive force field for hydrocarbons [J]. Journal of Physical Chemistry A, 2001, 105 (41):9396-9409.

[79] Strachan A, van Duin A C T, Chakraborty D, et al. Shock waves in high-energy materials: the initial chemical events in nitramine RDX [J]. Physical Reveiw Letter, 2003, 91: 098301.

[80] Strachan A, Kober E.M, van Duin Adri C T, et al. Reactive molecular dynamics [J]. Journal of Physical Chemistry, 2005, 122:054502.

[81] Li Z, Lang C, Chen W, et al. Molecular dynamics study of the effect of H2O on the thermal decomposition of α-phase CL-20 [J]. Acta Physial Chimica Sinica, 2013, 29(6):1145-1153.

[82] Zhang L, Zybin S, van Duin A C T, et al. Carbon cluster formation during thermal decomposition of octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine and 1,3,5-Triamino-2,4,6- trinitrobenzene high explosives from ReaxFF reactive molecular dynamics simulations [J]. Journal of Physical Chemistry A, 2009, 113: 10619-10640.

[83] Liu L, Liu Y, Zybin S V, et al. Correction of the ReaxFF reactive force field for London Dispersion, with applications to the equations of state for energetic materials [J]. Journal of Physical Chemistry A, 2011, 115 (40): 11016-11022.

[84] Chenoweth K, Cheung S, van Duin A C T, et al. Simulations on the thermal decomposition of a poly(dimethylsiloxane) polymer using the ReaxFF reactive force field [J]. J Am Chem Soc, 2005, 127 (19): 7192-7202.

[85] 張建文，王艷飛, 湍流化學(xué)反應(yīng)的分形數(shù)值模擬[J]. 火炸藥學(xué)報(bào), 2007, 30(3): 5-8, 15.

ZHANG Jian-wen., WANG Yan-fei. Numerical simulation for turbulent chemical reaction [J]. Chinese Journal Explosives & Propellants(Huozhayao Xuebao), 2007, 30(3):5-8, 15.

[86] 王國(guó)青,吳玉凱,侯慶偉,等.硫酸自由基與TNT反應(yīng)的密度泛函理論[J]. 火炸藥學(xué)報(bào), 2010, 33(2):10-12, 18.

WANG Guo-qing, WU Yu-kai, HOU Qing-wei, et al. A density functional theory of the reaction of TNT with sulfate radical [J]. Chinese Journal of Explosives & Propellants(Huozhayao Xuebao), 2010, 33(2):10-12, 18.

[87] 梁增友,黃風(fēng)雷,張振宇,等. PBX炸藥二維沖擊起爆機(jī)理的數(shù)值模擬[J]. 火炸藥學(xué)報(bào), 2008, 31(5):15-18.

LIAGN Zeng-you, HUANG Feng-lei, ZHANG Zhen-yu, et al. Numerical simulation of two dimension shock initiation and detonation behavior of PBX [J]. Chinese Journal of Explosives & Propellants(Huozhayao Xuebao), 2008, 31(5):15-18.

[88] 封鋒, 陳軍, 鄭亞,等.基于一維氣相穩(wěn)態(tài)反應(yīng)流的燃速預(yù)估軟件研究[J]. 火炸藥學(xué)報(bào), 2009, 32(3):58-61, 65.

FENG Feng, CHEN Jun, ZHENG Ya, et al. Study on software of solid propellant burning rate prediction based on one-dimensional steady-state reaction gas flow [J]. Chinese Journal of Explosives & Propellants(Huozhayao Xuebao),2009, 32(3):58-61, 65.

[89] 張江波,張玉成,蔣樹君,等.多層發(fā)射藥內(nèi)彈道模型及數(shù)值求解[J]. 火炸藥學(xué)報(bào), 2009, 32(3):83-86.

ZHAGN Jiang-bo, ZHANG Yu-cheng, JIANG Shu-jun,et al. The interior ballistic model and numerical simulation of multilayer gun propellant charge [J]. Chinese Journal of Explosives & Propellants(Huozhayao Xuebao), 2009, 32(3):83-86.

[90] 呂秉峰，劉幼平，董鳳云,等.定容條件下火藥實(shí)際燃燒規(guī)律的數(shù)值模擬[J]. 火炸藥學(xué)報(bào), 2007, 30(6):72-74.

Lü Bing-feng, LIU You-ping, DONG Feng-yun, et al. The mathematical model of propellant actual combustion law on the constant volume condition[J]. Chinese Journal of Explosives & Propellants(Huozhayao Xuebao), 2007, 30(6):72-74.

Research Progress in Theoretical Prediction of Physicochemical Properties for Energetic Materials

YAN Qi-long, SONG Zhen-wei, AN Ting, ZHANG Xiao-hong, ZHAO Feng-qi

(Science and Technology on Combustion and Explosion Laboratory, Xi′an Modern Chemistry Research Institute,Xi′an 710065, China)

Starting from the related achievements of recent development in the field of energetic materials, the main research direction in this field was discussed, and the latest achievements of theoretical prediction of the physicochemical properties of energetic materials were described in particular, mainly including quantum chemistry, molecular dynamics or semi empirical QSPR modeling approaches to predict the research progress in sensitivity, combustion and detonation performances, reaction activity, curing mechanism and mechanical properties of energetic materials.The main technical barriers existed at present were summarized, including the lack of complete and unified standard experimental database on performances of energetic materials, no commercial software with independent intellectual property rights to calculate the energetic material properties, and international commercial software with reliable predictors of physical and chemical properties of energetic materials is limited to the detonation performance and combustion performance. Literature research shows that China needs to further strengthen the research in this field, and finally build a comprehensive software platform that can evaluate the performance and security of energetic materials.With 90

.

10.14077/j.issn.1007-7812.2016.05.001

國(guó)家自然科學(xué)基金(No.21303133;No.21473130)

嚴(yán)啟龍(1983-)，男，博士，高級(jí)工程師，從事含能材料性能及表征研究。E-mail:terry.well@163.com

趙鳳起(1963-)，男，博士，博士生導(dǎo)師，研究員，從事固體推進(jìn)劑配方與性能研究。E-mail:zhaofqi@163.com

TJ55；TQ560

A

1007-7812(2016)05-0001-12