• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    不同組分比例的HMX/RDX混合物的熱穩(wěn)定性和力學(xué)性能研究

    2024-01-01 00:00:00李沙常雙君陳星妍魯娜
    爆破器材 2024年6期
    關(guān)鍵詞:模擬熱穩(wěn)定性力學(xué)性能

    [摘 要]為了研究不同組分比例的HMX/RDX混合物的熱穩(wěn)定性和力學(xué)性能,采用分子動(dòng)力學(xué)(MD)方法對(duì)不同HMX/RDX混合物在不同溫度下的引發(fā)鍵最大鍵長Lmax、內(nèi)聚能密度EC和力學(xué)性能參數(shù)進(jìn)行了計(jì)算;利用差示掃描量熱法(DSC)測(cè)試了不同HMX/RDX混合物在升溫速率為2、5、10、20 K/min下的熱分解性能,并采用EXPLO 5對(duì)不同HMX/RDX混合物的爆轟參數(shù)進(jìn)行了理論計(jì)算。結(jié)果表明:HMX/RDX中HMX與RDX的質(zhì)量比分別為90∶10、80∶20、70∶30時(shí),隨著溫度升高,HMX/RDX混合物的Lmax逐漸增大,EC逐漸減??;其中,HMX90/RDX10的Lmax較小,EC較大,熱感度較低。由DSC測(cè)試結(jié)果計(jì)算得到不同比例混合物的分解峰溫Tp0和熱爆炸臨界溫度Tb;其中,HMX90/RDX10的Tp0和Tb相對(duì)較高,分別為268.12℃和277.66℃,表明熱穩(wěn)定性較好。使用Forcite計(jì)算力學(xué)性能參數(shù),HMX90/RDX10的拉伸模量E、體積模量K、切變模量G相對(duì)較大,表明硬度和斷裂強(qiáng)度較大,力學(xué)性能相對(duì)較好。HMX90/RDX10的爆速、爆壓和爆轟總能量分別為9 192.19 m/s、38.43 GPa和10.98 kJ/cm3,能量相對(duì)較高。

    [關(guān)鍵詞]HMX/RDX混合物;熱穩(wěn)定性;力學(xué)性能;分子動(dòng)力學(xué)(MD)模擬

    [分類號(hào)]TQ560.7; TJ510

    Thermal Stability and Mechanical Properties of HMX/RDX Mixtures with

    Different Component Proportions

    LI Sha, CHANG Shuangjun, CHEN Xingyan, LU Na

    School of Environmental and Safety Engineering, North University of China (Shanxi Taiyuan, 030051)

    [ABSTRACT]The thermal stability and mechanical properties of HMX/RDX mixtures with different component proportions were studied. The maximum bond length Lmax, cohesive energy density EC , and mechanical properties of different HMX/RDX mixtures at different temperatures were calculated by molecular dynamics (MD) method. The thermal decomposition performance of different HMX/RDX mixtures was tested by differential scanning calorimetry (DSC) at heating rates of 2, 5, 10, and 20 K/min. And the detonation parameters of different HMX/RDX mixtures were theoretically calculated by EXPLO 5. The results show that when m(HMX)∶m(RDX) are 90∶10, 80∶20, and 70∶30, Lmax of the HMX/RDX mixtures gradually increases and EC gradually decreases with increasing temperature. Lmax and EC of HMX90/RDX10 are relatively small, with low thermal sensitivity and good thermal stability. The decomposition peak temperature Tp0 and thermal explosion critical temperature Tb of different mixtures were calculated based on DSC test results. Tp0 and Tb of HMX90/RDX10 are relatively high, which are 268.12℃ and 277.66℃ respectively, indicating good thermal stability. The mechanical property parameters were calculated by Forcite. The tensile modulus E, bulk modulus K, and shear modulus G of HMX90/RDX10 are relatively high, with high hardness and fracture strength, and relatively good mechanical pro-perties. The detonation velocity, detonation pressure, and total detonation energy of HMX90/RDX10 are 9 192.19 m/s, 38.43 GPa, and 10.98 kJ/cm3, respectively, indicating relatively high energy.

    [KEYWORDS]HMX/RDX mixture; thermal stability; mechanical property; molecular dynamics (MD) simulation

    0 引言

    奧克托今(HMX,C4H8N8O8)和黑索今(RDX,C3H6N6O6)是常見的高能單質(zhì)炸藥,具有高爆熱、高爆速和高爆壓等特點(diǎn),可作為混合炸藥、傳爆藥、固體推進(jìn)劑和發(fā)射藥配方中的高能組分[1,在國防、軍事和航天工業(yè)中發(fā)揮著重要作用。與RDX相比,HMX成本較高,將2種單質(zhì)炸藥混合使用,可降低成本。因此,需要對(duì)HMX和RDX混合物的相關(guān)性能進(jìn)行研究,得到一種熱穩(wěn)定性和力學(xué)性能較好、能量較高的HMX/RDX混合物。

    隨著計(jì)算機(jī)技術(shù)的發(fā)展,分子模擬方法應(yīng)用于含能材料分子間相互作用的研究取得了一定進(jìn)展。肖繼軍等[2提出的理論可以用來判斷HMX及其混合體系的熱穩(wěn)定性。Han等[3分析了摩爾比對(duì)HMX/MDNI(1-甲基-4,5-二硝基咪唑)體系的結(jié)合能和內(nèi)聚能密度的影響;結(jié)果表明,MDNI降低了HMX的感度。劉冬梅等[4通過分子動(dòng)力學(xué)(MD)方法研究了引發(fā)鍵最大鍵長Lmax、內(nèi)聚能密度EC與熱感度之間的關(guān)系,使用HMX、RDX和CL-20/HMX等晶體進(jìn)行理論計(jì)算,得出的結(jié)果和實(shí)驗(yàn)測(cè)試結(jié)果一致。悅征等[5通過差示掃描量熱法(DSC)測(cè)試了HMX/RDX復(fù)合物細(xì)化前、后的熱分解特性,得出細(xì)化后的熱穩(wěn)定性較好。Ye等[6采用MD模擬方法研究了HMX/RDX高能共晶和混合物模型的熱穩(wěn)定性和力學(xué)性能,得出共晶的熱穩(wěn)定性和力學(xué)性能優(yōu)于混合物模型。

    采用MD方法對(duì)不同組分比例的HMX/RDX混合物進(jìn)行理論計(jì)算,研究HMX/RDX混合物的熱穩(wěn)定性和力學(xué)性能;并通過DSC測(cè)試熱分解性能;采用EXPLO 5軟件對(duì)不同HMX/RDX混合物的爆轟參數(shù)進(jìn)行理論計(jì)算。通過以上3種方法對(duì)HMX/RDX混合物進(jìn)行研究,為HMX/RDX混合炸藥的配方設(shè)計(jì)提供理論參考。

    1 數(shù)值模擬

    1.1 HMX/RDX混合模型搭建

    依據(jù)CCDC(劍橋晶體數(shù)據(jù)中心)數(shù)據(jù)庫,獲得HMX和RDX的晶體結(jié)構(gòu),用做本研究的初始結(jié)構(gòu)。HMX的晶胞參數(shù)a=6.526 ?,b=11.037 ?,c=7.364 ?,α=γ=90°,β=102.67°,空間群14P21/N;RDX的晶胞參數(shù)a=13.182 ?,b=11.574 ?,c=10.709 ?,α=γ=β=90°,空間群61PBCA。HMX和RDX晶胞模型如圖1(a)所示。選用Compass[7力場(chǎng),在298 K采用Smart方法對(duì)HMX和RDX晶胞模型進(jìn)行結(jié)構(gòu)優(yōu)化。按照HMX、RDX以及HMX與RDX的質(zhì)量比分別為90∶10、80∶20、70∶30的要求構(gòu)建HMX和RDX超級(jí)晶胞,使用Build Layers工具建立RDX和HMX混合模型。HMX/RDX混合物中,HMX與RDX的質(zhì)量比為90∶10、80∶20和70∶30的樣品分別命名為HMX90/RDX10、HMX80/RDX20和HMX70/RDX30。構(gòu)建好的HMX90/RDX10混合物模型如圖1(b)所示。

    1.2 計(jì)算方法

    采用MD中的Forcite模塊對(duì)HMX、RDX、HMX90/RDX10、HMX80/RDX20和HMX70/RDX30共5種模型進(jìn)行幾何優(yōu)化,并對(duì)優(yōu)化后的結(jié)構(gòu)進(jìn)行理論計(jì)算。選用恒原子數(shù)、恒體積、恒溫(NVT)系綜,采用Andersen控溫方法,范德華和靜電作用分別用Atom-based和Ewald加和方法,截?cái)喟霃饺?2.5×10-10 m,并進(jìn)行截?cái)辔膊啃U?。溫度依次選取295、345、395、445、495 K。初始的分子運(yùn)動(dòng)速度按Maxwell-Boltzmann分布確定,積分采用Verlet方法。時(shí)間步長為1 fs,總模擬時(shí)間400 ps。其中,前200 ps用于體系平衡,后200 ps用于統(tǒng)計(jì)分析。

    2 熱分解實(shí)驗(yàn)

    2.1 材料與儀器

    HMX和RDX,甘肅銀光化學(xué)工業(yè)公司;HCT-1型微機(jī)質(zhì)量動(dòng)態(tài)分析儀,北京恒久實(shí)驗(yàn)設(shè)備有限公司?;旌衔镏?,HMX和RDX分別以90∶10、80∶20和70∶30的質(zhì)量比進(jìn)行混合。

    2.2 實(shí)驗(yàn)過程

    使用DSC對(duì)HMX、RDX、HMX90/RDX10、HMX80/RDX20和HMX70/RDX30的熱分解進(jìn)行測(cè)試。試樣質(zhì)量為3~5 mg,升溫區(qū)間為25~400℃,升溫速率為2、5、10、20 K/min。

    3 結(jié)果與討論

    3.1 體系平衡的判別

    HMX、RDX和HMX/RDX混合物的模擬在NVT系綜下進(jìn)行,當(dāng)溫度和能量在5%~10%波動(dòng)時(shí),可判定達(dá)到平衡[8。HMX90/RDX10混合物在298 K的NVT-MD模擬的能量和溫度曲線見圖2。

    圖2(a)中,混合體系的勢(shì)能和總能量分別在-29 662.24 kJ/mol和-20 193.87 kJ/mol上下波動(dòng),波動(dòng)幅度均小于10%,則該體系已達(dá)能量平衡;圖2(b)中,混合物的溫度波動(dòng)約±20 K左右,波動(dòng)幅度為6%,低于10%,表明溫度已達(dá)到平衡。其他4種模型采用同樣的方法進(jìn)行平衡判定。

    3.2 引發(fā)鍵鍵長

    引發(fā)鍵是高能化合物分子所有化學(xué)鍵中最弱的鍵。在外界作用刺激下,該鍵將最先斷裂,最易引發(fā)高能化合物的分解和起爆。N—NO2是硝銨類高能化合物熱分解或起爆的引發(fā)鍵[9。用MD方法計(jì)算引發(fā)鍵N—NO2鍵長的統(tǒng)計(jì)分布,對(duì)于討論HMX、RDX和HMX/RDX混合物的起爆和熱安全性具有重要意義。HMX、RDX和HMX/RDX混合物中N—NO2鍵的最大鍵長Lmax隨溫度的變化如圖3所示。

    從圖3中看出,在295~495 K溫度下,HMX、RDX和HMX/RDX混合物的引發(fā)鍵N—NO2鍵的Lmax隨著溫度的升高而增大。由此可見,引發(fā)鍵的Lmax可作為評(píng)判體系熱感度大小的理論依據(jù)之一。在相同溫度下,HMX/RDX中HMX與RDX的質(zhì)量比依次為90∶10、80∶20、70∶30時(shí),Lmax逐漸增大,更易斷裂;且3種HMX/RDX混合物的Lmax中,HMX90/RDX10的平均Lmax為1.620 ?,相對(duì)較小,表明HMX90/RDX10的熱感度較低,熱穩(wěn)定性較好,不易發(fā)生引發(fā)鍵的斷裂。

    3.3 內(nèi)聚能密度

    內(nèi)聚能密度EC是單位體積內(nèi)1 mol凝聚體克服分子間作用力變?yōu)闅鈶B(tài)時(shí)所需的能量[10。

    式中:ΔH為汽化熱;R為摩爾氣體常數(shù),8.314 J/(mol·K);T為溫度;Vm為摩爾體積。

    MD模擬中,EC是范德華力EvdW和靜電力Ee之和。不同溫度下HMX、RDX和HMX/RDX混合物的EC及分量見表1。

    由表1可知,隨著溫度的升高,HMX、RDX和HMX/RDX混合物的EC、EvdW和Ee均減小,表明由晶態(tài)變?yōu)闅鈶B(tài)所需能量變小,這與溫度升高、體系的熱感度增大的實(shí)驗(yàn)事實(shí)相符[4。所以,不同溫度下,HMX、RDX和HMX/RDX混合物中的EC在一定條件下可用于衡量體系熱感度相對(duì)大小的理論判據(jù)。HMX和RDX的平均EC分別為0.865 kJ/cm3和0.795 kJ/cm3,HMX的EC大于RDX,表明EC與熱感度之間存在負(fù)相關(guān)關(guān)系,所以HMX的熱感度比RDX低。在相同溫度下,HMX/RDX中HMX和RDX的質(zhì)量比依次為90∶10、80∶20、70∶30時(shí),EC逐漸減小。比較3種HMX/RDX混合物,HMX90/RDX10的熱感度相對(duì)較低,熱穩(wěn)定性較好。

    3.4 熱分解過程分析

    3.4.1 非等溫動(dòng)力學(xué)分析

    熱分解性能是含能材料的重要指標(biāo),采用DSC測(cè)試了HMX、RDX和HMX/RDX混合物在2、5、10 K/min和20 K/min升溫速率下的熱分解特性,并計(jì)算了活化能Ea。Ea可以用來表示化學(xué)反應(yīng)的難易程度,Ea越大,化學(xué)反應(yīng)越難進(jìn)行。對(duì)HMX與RDX

    混合前、后的熱穩(wěn)定性進(jìn)行研究,分別通過式(2)的Kissinger公式和式(3)的Ozawa公式計(jì)算了5種樣品在2、5、10 K/min和20 K/min升溫速率下的熱分解反應(yīng)動(dòng)力學(xué)參數(shù)。相應(yīng)數(shù)據(jù)如表2所示。

    式中:βi為升溫速率;Tpi為熱分解反應(yīng)的峰溫;A為指前因子;R為摩爾氣體常數(shù),8.314 J/(mol·K);EK為Kissinger法計(jì)算的活化能;EO為Ozawa法計(jì)算的活化能;G(α)為機(jī)理函數(shù)的積分式;α為熱分解反應(yīng)的程度或熱分解反應(yīng)的轉(zhuǎn)化率;Ea為EK和EO之和的平均值。

    從表2中看出,利用Kissinger法和Ozawa法計(jì)算出的5種樣品的熱分解表觀活化能一致性較好,線性相關(guān)系數(shù)RK或RO都接近0.99,說明計(jì)算結(jié)果的可信。當(dāng)HMX與RDX的質(zhì)量比分別為90∶10、80∶20、70∶30時(shí),Ea分別為289.78、360.09、436.63 kJ/mol,其中,HMX70/RDX30的Ea較大。HMX的熱分解則由RDX分解時(shí)放出的熱量所影響,分解機(jī)理可能發(fā)生了變化,從而影響了不同升溫速率下分解峰溫出現(xiàn)的先后順序。所以,需要進(jìn)一步計(jì)算出HMX/RDX混合物的初始分解峰溫和熱爆炸臨界溫度,分析熱穩(wěn)定性。

    3.4.2 熱力學(xué)參數(shù)計(jì)算

    HMX/RDX混合物在加熱速率趨于0時(shí)的初始分解峰溫Tp0的計(jì)算公式:

    Tpi=Tp0+bβi+cβ2i,i=1,2,3,4。(5)

    式中:b、c為常數(shù);βi為升溫速率;Tpi為不同升溫速率下的熱分解峰溫Tp。將數(shù)據(jù)βi和Tpi代入式(5),求得Tp0。

    采用Hu-Zhao-Gao法[11計(jì)算HMX/RDX混合物的熱爆炸臨界溫度Tb

    式中:A0、b為Berthelot方程的2個(gè)系數(shù);G(α)為積分形式的動(dòng)力學(xué)函數(shù)。

    以式(6)中的lnβi對(duì)Tpi作圖,通過線性擬合,由斜率求得b,代入式(7),求得Tb。

    根據(jù)非等溫得到的動(dòng)力學(xué)參數(shù)和熱力學(xué)關(guān)系式計(jì)算活化焓ΔH、活化熵ΔS及活化自由能ΔG等熱力學(xué)參數(shù)12

    式中:KB為Boltzmann常數(shù),1.380 7×10-23J/K;h為Plank常數(shù),6.626×10-34 J·s;ΔG為活化自由能;ΔH為活化焓;ΔS為活化熵;Ea為Ozaw法與Kissinger法得出的表觀活化能的平均值;A為指前因子。根據(jù)式(5)~式(10),分別計(jì)算出Tp0、Tb、ΔG、ΔH和ΔS,見表3。

    由表3可知,HMX的Tp0和Tb分別為268.97℃和277.98℃,比RDX大,所以HMX的熱穩(wěn)定性較好。HMX/RDX中HMX與RDX的質(zhì)量比分別為90∶10、80∶20、70∶30時(shí),對(duì)3種HMX/RDX混合物進(jìn)行比較。其中,HMX90/RDX10的Tp0和Tb相對(duì)較大,分別為268.12℃和277.66℃,熱穩(wěn)定性較好,常溫下不易發(fā)生分解。

    3.5 力學(xué)性能

    力學(xué)性能是高能材料最重要的性能之一。依據(jù)Voigt理論與Reuss理論方法,根據(jù)拉伸模量E、體積模量K、切變模量G、泊松比ν之間的關(guān)系,得到E與ν的計(jì)算公式:

    模量可以用來衡量剛性的強(qiáng)弱,材料的硬度與E和G密切相關(guān);E和G越大,材料硬度越高。材料的斷裂強(qiáng)度的大小由K度量。K越大,表明材料的斷裂強(qiáng)度越大[13-14。表4中,列出了不同HMX/RDX混合物在NVT系綜下MD模擬得到的力學(xué)性能參數(shù)。

    表4中,ν在0.330~0.379之間,屬于0.200~0.400的范圍內(nèi),說明HMX/RDX混合物具有一定的塑性。溫度的變化將會(huì)對(duì)材料硬度和斷裂強(qiáng)度的強(qiáng)弱產(chǎn)生影響。HMX/RDX混合物的E、K和G隨著溫度的升高而減小,則HMX/RDX混合物的硬度和斷裂強(qiáng)度減弱;在相同溫度下,當(dāng)HMX/RDX中HMX與RDX質(zhì)量比為90∶10時(shí),E、K和G較高,硬度和斷裂強(qiáng)度較大,力學(xué)性能相對(duì)較好。

    3.6 爆轟參數(shù)

    由非理想氣體的維里(Virial)理論方程推導(dǎo)出VLW狀態(tài)方程的表達(dá)式[15

    式中:B*為無量綱第二維里系數(shù);T*為無量綱溫度,T*=KT/ε;b0=0.67πNAσ3;ε和σ為Lennard-Jones勢(shì)參數(shù)。

    在計(jì)算HMX/RDX混合物的爆熱、爆壓、爆溫、爆速和爆轟總能量等爆轟參數(shù)時(shí),除VLW狀態(tài)方程外,還需要下列方程組:

    式中:D為C-J爆轟狀態(tài)下的爆速;u為C-J爆轟狀態(tài)下的質(zhì)點(diǎn)速度;E為C-J爆轟狀態(tài)下的內(nèi)能;ρ為C-J爆轟狀態(tài)下的密度;c為C-J面音速;下標(biāo)“0”為初始狀態(tài)下的對(duì)應(yīng)值。

    基于VLW狀態(tài)方程,采用EXPLO 5軟件對(duì)不同HMX/RDX混合物的爆轟參數(shù)進(jìn)行理論計(jì)算,結(jié)果見表5。

    從表5計(jì)算結(jié)果可知,3種HMX/RDX混合物中,HMX90/RDX10的爆速為9 192.19 m/s,爆壓為38.43 GPa,爆轟總能量為10.98 kJ/cm3,爆轟能量相對(duì)較高。

    4 結(jié)論

    采用MD方法計(jì)算不同HMX/RDX混合物的Lmax、EC和力學(xué)性能參數(shù),通過DSC測(cè)試并計(jì)算了Tp0和Tb,并采用EXPLO 5軟件對(duì)爆轟參數(shù)進(jìn)行了理論計(jì)算,得出的結(jié)論如下:

    1)隨著溫度的升高,HMX/RDX混合物的Lmax逐漸增大,EC逐漸減小,表明熱感度逐漸增大。相同溫度下,HMX90/RDX10的Lmax較小,EC較大,熱感度較小。由DSC測(cè)試結(jié)果可知,3種HMX/RDX混合物中,HMX90/RDX10的Tp0和Tb相對(duì)較高,分別為268.12℃和277.66℃,表明熱穩(wěn)定性較好。

    2)溫度升高,HMX/RDX混合物的E、K和G逐漸減小,表明材料的剛性減弱,硬度和斷裂強(qiáng)度降低。在相同溫度下,HMX90/RDX10的E、K和G相對(duì)較大,硬度和斷裂強(qiáng)度較大,力學(xué)性能相對(duì)較好。

    3)EXPLO 5計(jì)算結(jié)果表明,3種不同HMX/RDX混合物中,HMX90/RDX10的爆速、爆壓和爆轟總能量分別為 9 192.19 m/s、38.43 GPa和10.98 kJ/cm3,能量相對(duì)較高。

    參考文獻(xiàn)

    [1]GE N N, WEI Y K, JI G F, et al. Initial decomposition of the condensed β-HMX under shock waves: molecular dynamics simulations [J]. Journal of Physical Chemistry B, 2012, 116(46): 13696-13704.

    [2]朱偉, 劉冬梅, 肖繼軍, 等. 多組分高能復(fù)合體系的感度判據(jù)、熱膨脹和力學(xué)性能的MD研究[J]. 含能材料, 2014, 22(5): 582-587.

    ZHU W, LIU D M, XIAO J J, et al. Molecular dynamics study on sensitivity criterion,thermal expansion and mechanical properties of multi-component high energy systems [J]. Chinese Journal of Energetic Materials, 2014, 22(5): 582-587.

    [3]HAN G, GOU R J, REN F D, et al. Theoretical investigation into the influence of molar ratio on binding energy, mechanical property and detonation performance of 1, 3, 5, 7-tetranitro-1, 3, 5, 7-tetrazacyclo octane (HMX)/1-methyl-4, 5-dinitroimidazole (MDNI) cocrystal explosive[J]. Computational and Theoretical Chemistry, 2017, 1109: 27-35.

    [4]劉冬梅, 趙麗, 肖繼軍, 等. 不同溫度下HMX和RDX晶體的感度判別和力學(xué)性能預(yù)估: 分子動(dòng)力學(xué)比較研究[J]. 高等學(xué)?;瘜W(xué)學(xué)報(bào), 2013, 34(11): 2558-2565.

    LIU D M, ZHAO L, XIAO J J, et al. Sensitivity crite-rion and mechanical properties prediction of HMX and RDX crystals at different temperatures:comparative study with molecular dymamics simulation [J]. Chemistry Journal of Chinese Universities, 2013, 34(11): 2558-2565.

    [5]悅征, 李志華, 魏成龍, 等. HMX/RDX復(fù)合物細(xì)化與熱分解特性研究[J]. 火工品, 2021(1): 39-42.

    YUE Z, LI Z H, WEI C L, et al. Refining and thermal decomposition of HMX/RDX complex [J]. Initiators amp; Pyrotechnics, 2021(1): 39-42.

    [6]YE B S, GONG J, HU X Y, et al. Comparative investigation on the thermostability, sensitivity, and mechanical performance of RDX/HMX energetic cocrystal and its mixture[J]. Journal of Molecular Modeling, 2020, 26(7): 176.

    [7]陶俊, 王曉峰, 趙省向, 等. CL-20/HMX共晶與共混物的分子動(dòng)力學(xué)模擬[J].含能材料, 2016, 24(4): 324-330.

    TAO J, WANG X F, ZHAO S X, et al. Molecular dynamics simulation of CL-20/HMX cocrystal and blends[J]. Chinese Journal of Energetic Materials, 2016, 24(4): 324-330.

    [8]ZHANG L, YU Y, XIANG M Z. A study of the shock sensitivity of energetic single crystals by large-scale ab initio molecular dynamics simulations[J]. Nanomate-rials, 2019, 9(9): 1251.

    [9]高杰, 王紅星, 金大勇, 等. DNAN/DNIF二元共熔物在熱和機(jī)械刺激下的安全性研究[J]. 爆破器材, 2021, 50(3): 35-39.

    GAO J, WANG H X, JIN D Y, et al. Safety of DNAN/DNIF binary eutectic under thermal stimulation and mechanical stimulation [J]. Explosive Materials, 2021, 50(3): 35-39.

    [10]HANG G Y, YU W L, WANG T, et al. Theoretical investigations into effects of adulteration crystal defects on properties of CL-20/TNT cocrystal explosive [J]. Computational Materials Science, 2019, 156: 77-83.

    [11]趙宏安, 胡榮祖, 王喜軍, 等. 1,3,3-三硝基氮雜環(huán)丁烷的熱安全性[J]. 化學(xué)學(xué)報(bào), 2009, 67(22): 2536-2540.

    ZHAO H A, HU R Z, WANG X J, et al. Thermal safety of 1,3,3-trinitroazetidine [J]. Acta Chimica Sinica, 2009, 67(22): 2536-2540.

    [12]胡榮祖, 趙鳳起, 高紅旭, 等.量熱學(xué)基礎(chǔ)與應(yīng)用[M]. 北京: 科學(xué)出版社, 2011.

    [13]楊犁, 余庚澤, 余晨, 等. 基于TKX-50的PBXs含能材料力學(xué)性能計(jì)算模擬[J]. 實(shí)驗(yàn)技術(shù)與管理, 2021, 38(3): 61-68.

    YANG L, YU G Z, YU C, et al. Computational simulation of mechanical properties of PBXs energetic materials based on TKX-50 [J]. Experiment Technology and Management, 2021, 38(3): 61-68.

    [14]孫婷, 肖繼軍, 趙鋒, 等. CL-20/DNB共晶基PBXs兼容性、界面作用和力學(xué)性能的MD模擬[J]. 含能材料, 2015, 23(4): 309-314.

    SUN T, XIAO J J, ZHAO F, et al. Molecular dyna-mics simulation of compatibility, interface interactions and mechanical properties of CL-20/DNB cocrystal based PBXs[J]. Chinese Journal of Energetic Mate-rials, 2015, 23(4): 309-314.

    [15]吳雄, 龍新平,何碧, 等. VLW爆轟產(chǎn)物狀態(tài)方程[J]. 中國科學(xué)B輯:化學(xué), 2008, 38(12): 1129-1132.

    WU X, LONG X P, HE B, et al. VLW equation of state for detonation products [J]. Science in China Series B: Chemistry, 2008, 38(12): 1129-1132.

    猜你喜歡
    模擬熱穩(wěn)定性力學(xué)性能
    Pr對(duì)20MnSi力學(xué)性能的影響
    云南化工(2021年11期)2022-01-12 06:06:14
    Mn-Si對(duì)ZG1Cr11Ni2WMoV鋼力學(xué)性能的影響
    山東冶金(2019年3期)2019-07-10 00:54:00
    一個(gè)高分子模擬計(jì)算網(wǎng)格的作業(yè)管理
    工業(yè)機(jī)器人模擬仿真技術(shù)在職業(yè)教育中的應(yīng)用淺析
    淺析柔道運(yùn)動(dòng)員的模擬實(shí)戰(zhàn)訓(xùn)練
    虛擬機(jī)局域網(wǎng)組建技術(shù)應(yīng)用初探
    PVC用酪氨酸鑭的合成、復(fù)配及熱穩(wěn)定性能研究
    中國塑料(2016年7期)2016-04-16 05:25:52
    提高有機(jī)過氧化物熱穩(wěn)定性的方法
    可聚合松香衍生物的合成、表征和熱穩(wěn)定性?
    INCONEL625+X65復(fù)合管的焊接組織與力學(xué)性能
    焊接(2015年9期)2015-07-18 11:03:53
    国产成人a∨麻豆精品| 久久久久久久国产电影| kizo精华| 亚洲av欧美aⅴ国产| 成人国语在线视频| 男男h啪啪无遮挡| 国产一级毛片在线| 久久久精品免费免费高清| 成人亚洲欧美一区二区av| 后天国语完整版免费观看| av又黄又爽大尺度在线免费看| 大片免费播放器 马上看| 亚洲欧美日韩另类电影网站| 精品视频人人做人人爽| 国产av精品麻豆| 操出白浆在线播放| 久久人人97超碰香蕉20202| 999久久久国产精品视频| 深夜精品福利| 免费观看a级毛片全部| 精品一区在线观看国产| 高清欧美精品videossex| 亚洲精品国产av蜜桃| 在线精品无人区一区二区三| 中文字幕人妻熟女乱码| 日日爽夜夜爽网站| av线在线观看网站| 夫妻性生交免费视频一级片| 日韩伦理黄色片| 国产成人av激情在线播放| 又紧又爽又黄一区二区| 高清视频免费观看一区二区| 一区在线观看完整版| 国产日韩一区二区三区精品不卡| 狂野欧美激情性bbbbbb| a级片在线免费高清观看视频| 天天操日日干夜夜撸| 91麻豆av在线| 亚洲精品一卡2卡三卡4卡5卡 | 天天躁狠狠躁夜夜躁狠狠躁| 满18在线观看网站| 男女边吃奶边做爰视频| 国产黄色视频一区二区在线观看| 熟女av电影| 婷婷色麻豆天堂久久| 国产精品 欧美亚洲| 成年人午夜在线观看视频| 一个人免费看片子| 国产精品久久久久久人妻精品电影 | 女人被躁到高潮嗷嗷叫费观| 国产三级黄色录像| 少妇人妻久久综合中文| 99九九在线精品视频| 亚洲图色成人| 国产精品一区二区在线观看99| 国产国语露脸激情在线看| 亚洲精品av麻豆狂野| 丝袜在线中文字幕| 国产高清不卡午夜福利| 中文精品一卡2卡3卡4更新| 欧美性长视频在线观看| 精品一品国产午夜福利视频| 人成视频在线观看免费观看| 精品第一国产精品| 中文字幕亚洲精品专区| 99国产综合亚洲精品| 国产麻豆69| 9191精品国产免费久久| 亚洲国产毛片av蜜桃av| 99国产精品一区二区蜜桃av | 在线精品无人区一区二区三| 日本色播在线视频| xxxhd国产人妻xxx| 精品人妻一区二区三区麻豆| 国产女主播在线喷水免费视频网站| 啦啦啦在线观看免费高清www| 超碰成人久久| 午夜福利一区二区在线看| 麻豆av在线久日| 欧美精品av麻豆av| 精品国产超薄肉色丝袜足j| 午夜福利,免费看| 制服诱惑二区| 91九色精品人成在线观看| 国产亚洲精品第一综合不卡| 男的添女的下面高潮视频| 黑人猛操日本美女一级片| 侵犯人妻中文字幕一二三四区| 热99久久久久精品小说推荐| av福利片在线| 免费在线观看视频国产中文字幕亚洲 | 赤兔流量卡办理| 丝袜人妻中文字幕| 国产日韩一区二区三区精品不卡| 成人国产一区最新在线观看 | 777久久人妻少妇嫩草av网站| 人妻 亚洲 视频| 久久精品成人免费网站| 日日夜夜操网爽| 亚洲av成人不卡在线观看播放网 | 国产精品亚洲av一区麻豆| 国产在视频线精品| 男女无遮挡免费网站观看| 一级,二级,三级黄色视频| 久久午夜综合久久蜜桃| 99久久人妻综合| 肉色欧美久久久久久久蜜桃| 男女国产视频网站| 欧美日韩亚洲高清精品| 免费观看av网站的网址| √禁漫天堂资源中文www| 97精品久久久久久久久久精品| 尾随美女入室| 亚洲伊人久久精品综合| 手机成人av网站| 国产欧美日韩一区二区三 | 亚洲伊人久久精品综合| 国产精品一区二区精品视频观看| 日韩,欧美,国产一区二区三区| 黄色片一级片一级黄色片| 精品国产国语对白av| 日本91视频免费播放| 纵有疾风起免费观看全集完整版| 校园人妻丝袜中文字幕| 午夜免费成人在线视频| 麻豆av在线久日| 欧美日韩视频高清一区二区三区二| 欧美 日韩 精品 国产| 90打野战视频偷拍视频| 黑人猛操日本美女一级片| 丰满饥渴人妻一区二区三| 日韩,欧美,国产一区二区三区| 亚洲国产av新网站| 国产精品99久久99久久久不卡| 国产精品国产av在线观看| 日本av免费视频播放| 在线观看免费高清a一片| 欧美av亚洲av综合av国产av| 国产精品成人在线| 一本大道久久a久久精品| 在线观看免费视频网站a站| 国产精品成人在线| 2021少妇久久久久久久久久久| 在线观看免费视频网站a站| 大香蕉久久成人网| 午夜福利影视在线免费观看| 精品人妻在线不人妻| 免费人妻精品一区二区三区视频| 欧美乱码精品一区二区三区| 赤兔流量卡办理| 国产熟女午夜一区二区三区| 热99国产精品久久久久久7| 亚洲人成网站在线观看播放| 99精品久久久久人妻精品| 丝袜美足系列| 国产午夜精品一二区理论片| 久9热在线精品视频| 午夜日韩欧美国产| 1024香蕉在线观看| av片东京热男人的天堂| 国产日韩欧美视频二区| 欧美激情高清一区二区三区| 中文字幕最新亚洲高清| 免费在线观看影片大全网站 | av一本久久久久| 色播在线永久视频| 男人添女人高潮全过程视频| 交换朋友夫妻互换小说| 欧美另类一区| 成年av动漫网址| 满18在线观看网站| 91精品伊人久久大香线蕉| 午夜影院在线不卡| 操出白浆在线播放| 国产亚洲av高清不卡| 国产亚洲欧美在线一区二区| 亚洲精品国产色婷婷电影| 午夜视频精品福利| 久久免费观看电影| 精品少妇内射三级| 欧美激情 高清一区二区三区| 一区二区av电影网| 久久精品久久久久久久性| 久久久国产精品麻豆| 777米奇影视久久| 久久精品久久久久久久性| 美女脱内裤让男人舔精品视频| 日本猛色少妇xxxxx猛交久久| 免费看av在线观看网站| 欧美精品人与动牲交sv欧美| 日日爽夜夜爽网站| 少妇被粗大的猛进出69影院| 最新在线观看一区二区三区 | 交换朋友夫妻互换小说| netflix在线观看网站| 在线观看免费视频网站a站| 国产熟女午夜一区二区三区| 久久国产精品影院| 你懂的网址亚洲精品在线观看| 国产高清国产精品国产三级| av线在线观看网站| 亚洲色图 男人天堂 中文字幕| 最近手机中文字幕大全| 欧美激情极品国产一区二区三区| 丰满人妻熟妇乱又伦精品不卡| 飞空精品影院首页| 免费看十八禁软件| www日本在线高清视频| 91精品三级在线观看| 少妇的丰满在线观看| 啦啦啦啦在线视频资源| 亚洲欧美精品自产自拍| 亚洲 欧美一区二区三区| 国产在线观看jvid| 日韩制服骚丝袜av| 亚洲精品中文字幕在线视频| 日韩人妻精品一区2区三区| 久久 成人 亚洲| av欧美777| 免费在线观看黄色视频的| 丝袜脚勾引网站| 五月开心婷婷网| 午夜久久久在线观看| 成年av动漫网址| 色精品久久人妻99蜜桃| 夜夜骑夜夜射夜夜干| 亚洲男人天堂网一区| 首页视频小说图片口味搜索 | 成人国产av品久久久| 一级毛片女人18水好多 | 天堂8中文在线网| 亚洲国产毛片av蜜桃av| 一区二区av电影网| 欧美成人精品欧美一级黄| 欧美精品一区二区大全| 天天躁日日躁夜夜躁夜夜| 亚洲av美国av| 精品一区在线观看国产| 国产野战对白在线观看| 日韩伦理黄色片| 97精品久久久久久久久久精品| 成人亚洲欧美一区二区av| 久久国产精品人妻蜜桃| 涩涩av久久男人的天堂| 香蕉国产在线看| 亚洲精品日本国产第一区| 大香蕉久久网| 青青草视频在线视频观看| 日韩视频在线欧美| 一本色道久久久久久精品综合| 国产亚洲一区二区精品| 一级片免费观看大全| 精品卡一卡二卡四卡免费| 久久天堂一区二区三区四区| 2018国产大陆天天弄谢| 波多野结衣av一区二区av| 婷婷色综合大香蕉| 日韩 亚洲 欧美在线| 成人亚洲欧美一区二区av| 黄色a级毛片大全视频| 男女无遮挡免费网站观看| 国产成人一区二区三区免费视频网站 | 国产又色又爽无遮挡免| 熟女少妇亚洲综合色aaa.| 国产深夜福利视频在线观看| 精品国产乱码久久久久久小说| 水蜜桃什么品种好| 高清黄色对白视频在线免费看| xxxhd国产人妻xxx| 欧美人与善性xxx| 欧美中文综合在线视频| 91成人精品电影| av福利片在线| 老司机靠b影院| 亚洲av国产av综合av卡| 国产精品99久久99久久久不卡| 国产精品av久久久久免费| 免费在线观看日本一区| 国产精品久久久人人做人人爽| 日韩大片免费观看网站| 午夜免费成人在线视频| 国产麻豆69| 热re99久久精品国产66热6| 男的添女的下面高潮视频| 女人爽到高潮嗷嗷叫在线视频| 99国产精品免费福利视频| www.精华液| 丁香六月天网| 操美女的视频在线观看| 午夜免费成人在线视频| 菩萨蛮人人尽说江南好唐韦庄| 另类精品久久| 大香蕉久久成人网| 韩国高清视频一区二区三区| 天天躁夜夜躁狠狠久久av| 亚洲精品国产av成人精品| av有码第一页| 久久青草综合色| av视频免费观看在线观看| 中文字幕色久视频| av在线老鸭窝| 国产成人91sexporn| 欧美 日韩 精品 国产| 国产成人影院久久av| 国产精品久久久久久精品电影小说| 美女脱内裤让男人舔精品视频| 黄色毛片三级朝国网站| 在线天堂中文资源库| 久久影院123| 色网站视频免费| 亚洲欧美日韩高清在线视频 | 国产淫语在线视频| 日本vs欧美在线观看视频| 激情视频va一区二区三区| 亚洲黑人精品在线| 热99国产精品久久久久久7| 色网站视频免费| 美女福利国产在线| 国产成人精品在线电影| 丝袜美足系列| 国产男女超爽视频在线观看| 国产成人av激情在线播放| 国产精品国产av在线观看| 一级片'在线观看视频| 日本91视频免费播放| 亚洲国产av新网站| 女人高潮潮喷娇喘18禁视频| 日本vs欧美在线观看视频| 婷婷色av中文字幕| 亚洲av国产av综合av卡| 国产精品国产三级国产专区5o| 男的添女的下面高潮视频| 在现免费观看毛片| 免费女性裸体啪啪无遮挡网站| 亚洲欧美一区二区三区久久| 欧美日韩综合久久久久久| 男男h啪啪无遮挡| 岛国毛片在线播放| 91精品伊人久久大香线蕉| 19禁男女啪啪无遮挡网站| 校园人妻丝袜中文字幕| 19禁男女啪啪无遮挡网站| 人人妻人人添人人爽欧美一区卜| 中文字幕人妻丝袜一区二区| 久久久久久久精品精品| 欧美日本中文国产一区发布| 中文字幕av电影在线播放| 91精品三级在线观看| 91老司机精品| 黄色视频不卡| 91精品国产国语对白视频| 99国产精品免费福利视频| 久久综合国产亚洲精品| 捣出白浆h1v1| 美女福利国产在线| 最近中文字幕2019免费版| av国产久精品久网站免费入址| 亚洲中文av在线| 看免费av毛片| 19禁男女啪啪无遮挡网站| 脱女人内裤的视频| 久久精品人人爽人人爽视色| 美女福利国产在线| 国产精品av久久久久免费| 天天躁日日躁夜夜躁夜夜| 国产高清不卡午夜福利| 亚洲自偷自拍图片 自拍| 国产伦人伦偷精品视频| 免费日韩欧美在线观看| 久久久久久久久久久久大奶| 亚洲欧美日韩另类电影网站| 欧美精品一区二区免费开放| 欧美成人精品欧美一级黄| 日本午夜av视频| 999精品在线视频| 午夜av观看不卡| 999精品在线视频| 久9热在线精品视频| 欧美成人精品欧美一级黄| 别揉我奶头~嗯~啊~动态视频 | 午夜福利免费观看在线| netflix在线观看网站| 成人亚洲欧美一区二区av| 一区二区三区乱码不卡18| 国产男女超爽视频在线观看| av不卡在线播放| 满18在线观看网站| 又紧又爽又黄一区二区| 亚洲色图 男人天堂 中文字幕| 老汉色∧v一级毛片| 久9热在线精品视频| 午夜视频精品福利| 少妇人妻 视频| 精品亚洲乱码少妇综合久久| 在线av久久热| avwww免费| svipshipincom国产片| 亚洲视频免费观看视频| 黄色a级毛片大全视频| 国产一区二区三区综合在线观看| 亚洲中文日韩欧美视频| 久久人妻熟女aⅴ| 亚洲一区中文字幕在线| 交换朋友夫妻互换小说| 黄网站色视频无遮挡免费观看| 亚洲美女黄色视频免费看| 婷婷成人精品国产| 久久久精品国产亚洲av高清涩受| 天天躁狠狠躁夜夜躁狠狠躁| 夫妻午夜视频| 欧美黄色片欧美黄色片| 脱女人内裤的视频| 欧美日韩综合久久久久久| 高清欧美精品videossex| 久久99一区二区三区| 国产精品99久久99久久久不卡| 丝袜人妻中文字幕| 欧美日韩亚洲高清精品| 一级毛片 在线播放| 搡老岳熟女国产| av网站在线播放免费| 两个人免费观看高清视频| 在线观看国产h片| 国产精品二区激情视频| 亚洲中文日韩欧美视频| 大话2 男鬼变身卡| 一边摸一边做爽爽视频免费| 99香蕉大伊视频| 人人妻,人人澡人人爽秒播 | 91成人精品电影| 又粗又硬又长又爽又黄的视频| 久久精品aⅴ一区二区三区四区| 久久ye,这里只有精品| 欧美人与善性xxx| 99国产精品99久久久久| 亚洲国产精品一区二区三区在线| 国产亚洲av高清不卡| 日本色播在线视频| 涩涩av久久男人的天堂| 多毛熟女@视频| 亚洲,欧美,日韩| 人人妻人人添人人爽欧美一区卜| 天天添夜夜摸| 精品福利观看| 成年人黄色毛片网站| 男人舔女人的私密视频| 亚洲国产欧美在线一区| 宅男免费午夜| 亚洲中文日韩欧美视频| 99国产精品一区二区蜜桃av | 日本vs欧美在线观看视频| 免费少妇av软件| 欧美精品av麻豆av| 国产精品.久久久| 欧美人与性动交α欧美软件| 这个男人来自地球电影免费观看| 国产欧美日韩一区二区三区在线| 成年动漫av网址| 精品熟女少妇八av免费久了| 丝袜喷水一区| 久久久久精品国产欧美久久久 | 一区二区日韩欧美中文字幕| 久久久久久久精品精品| 久久人妻福利社区极品人妻图片 | 婷婷色综合大香蕉| 国产精品麻豆人妻色哟哟久久| 午夜久久久在线观看| 人妻一区二区av| 精品高清国产在线一区| 天天操日日干夜夜撸| 在线观看www视频免费| 精品国产一区二区久久| 日韩精品免费视频一区二区三区| 欧美日韩视频高清一区二区三区二| 国产成人精品在线电影| a 毛片基地| 欧美成狂野欧美在线观看| 人人妻,人人澡人人爽秒播 | 亚洲第一av免费看| 国产精品一区二区免费欧美 | 在线av久久热| 18禁国产床啪视频网站| 99国产精品99久久久久| 国产精品久久久久久人妻精品电影 | 欧美另类一区| 女人久久www免费人成看片| 中文字幕精品免费在线观看视频| 视频在线观看一区二区三区| 成人午夜精彩视频在线观看| av在线播放精品| 你懂的网址亚洲精品在线观看| 亚洲专区国产一区二区| 黄色片一级片一级黄色片| 国产xxxxx性猛交| 中国国产av一级| 午夜福利在线免费观看网站| 97人妻天天添夜夜摸| 婷婷色综合大香蕉| 国产真人三级小视频在线观看| 亚洲情色 制服丝袜| 色网站视频免费| 久久久久国产精品人妻一区二区| 黄色片一级片一级黄色片| 美女午夜性视频免费| 天堂8中文在线网| 国产熟女欧美一区二区| 亚洲中文字幕日韩| 七月丁香在线播放| 精品国产一区二区三区四区第35| 我的亚洲天堂| 欧美xxⅹ黑人| 亚洲精品中文字幕在线视频| 久久久久视频综合| 中文字幕人妻丝袜一区二区| 我要看黄色一级片免费的| 最新的欧美精品一区二区| 成年av动漫网址| 国产激情久久老熟女| 啦啦啦 在线观看视频| av不卡在线播放| 欧美成狂野欧美在线观看| 亚洲成人免费av在线播放| 久久午夜综合久久蜜桃| 叶爱在线成人免费视频播放| 亚洲欧美激情在线| 免费不卡黄色视频| 蜜桃在线观看..| 亚洲中文日韩欧美视频| 日韩av不卡免费在线播放| 久久人人97超碰香蕉20202| www.精华液| 女性生殖器流出的白浆| 亚洲欧美一区二区三区国产| 久久这里只有精品19| 亚洲人成网站在线观看播放| netflix在线观看网站| 日韩大码丰满熟妇| 亚洲国产精品国产精品| 男女国产视频网站| 大陆偷拍与自拍| 一边摸一边做爽爽视频免费| e午夜精品久久久久久久| 久久毛片免费看一区二区三区| 亚洲国产欧美网| 国产黄色免费在线视频| 新久久久久国产一级毛片| 久久影院123| 一级毛片黄色毛片免费观看视频| 欧美日韩黄片免| 亚洲国产欧美在线一区| 啦啦啦 在线观看视频| 亚洲国产精品成人久久小说| 黄色片一级片一级黄色片| av视频免费观看在线观看| netflix在线观看网站| 欧美亚洲 丝袜 人妻 在线| 男女午夜视频在线观看| 人人澡人人妻人| 日本a在线网址| 晚上一个人看的免费电影| 高清av免费在线| 成人免费观看视频高清| 一二三四在线观看免费中文在| 国产精品秋霞免费鲁丝片| 亚洲人成网站在线观看播放| 美国免费a级毛片| 日本av免费视频播放| 精品国产一区二区三区四区第35| 国产1区2区3区精品| 亚洲一区二区三区欧美精品| av天堂在线播放| 久久国产精品人妻蜜桃| 亚洲七黄色美女视频| 三上悠亚av全集在线观看| 国产亚洲午夜精品一区二区久久| 成人免费观看视频高清| 亚洲欧美激情在线| 精品久久蜜臀av无| av有码第一页| 国产日韩欧美视频二区| 可以免费在线观看a视频的电影网站| 捣出白浆h1v1| 国产黄色视频一区二区在线观看| 两人在一起打扑克的视频| 国产成人一区二区在线| 亚洲精品国产av成人精品| 亚洲天堂av无毛| 一区二区日韩欧美中文字幕| 女人爽到高潮嗷嗷叫在线视频| 国产一区二区三区av在线| 一区福利在线观看| 精品国产国语对白av| 两性夫妻黄色片| 精品一品国产午夜福利视频| 亚洲国产最新在线播放| 波多野结衣一区麻豆| 一区福利在线观看| 精品久久久精品久久久| 日本午夜av视频| 少妇被粗大的猛进出69影院| 最近手机中文字幕大全| 中文字幕人妻熟女乱码| 免费看不卡的av| 欧美国产精品va在线观看不卡| 免费在线观看黄色视频的| 宅男免费午夜| 亚洲一区中文字幕在线| 少妇人妻久久综合中文| 精品人妻熟女毛片av久久网站| 国产精品免费视频内射| 狂野欧美激情性xxxx| 久久精品人人爽人人爽视色| 啦啦啦啦在线视频资源| 国产亚洲午夜精品一区二区久久| 欧美日韩亚洲综合一区二区三区_|