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

    1H-1,2,4-三唑衍生物爆轟性能的理論計算

    2014-08-22 08:21:24馬忠亮王建龍李景滿
    火炸藥學(xué)報 2014年4期
    關(guān)鍵詞:三唑熱穩(wěn)定性基團

    馬忠亮,周 華,王建龍,李景滿

    (1.中北大學(xué)化工與環(huán)境學(xué)院,山西太原030051;2.四川通達(dá)化工有限責(zé)任公司,四川達(dá)州635000)

    1H-1,2,4-三唑衍生物爆轟性能的理論計算

    馬忠亮1,周 華1,王建龍1,李景滿2

    (1.中北大學(xué)化工與環(huán)境學(xué)院,山西太原030051;2.四川通達(dá)化工有限責(zé)任公司,四川達(dá)州635000)

    采用密度泛函理論(DFT)B3LYP方法,在6-31G(d,p)基組水平下對1H-1,2,4-三唑衍生物進(jìn)行幾何構(gòu)型全優(yōu)化,設(shè)計了合適的等鍵反應(yīng)計算1 H-1,2,4-三唑衍生物的生成焓?;?.001e·bohr-3等電子密度面所包圍的體積空間求得分子平均摩爾體積(V)和理論密度(ρ),運用Kamlet-Jacobs方程估算1H-1,2,4-三唑衍生物的爆速(D)和爆壓(p),結(jié)果表明,1H-1,2,4-三唑發(fā)生N-取代或引入-N3基團都能顯著提高衍生物的生成焓,-NF2或-ONO2基團能夠增強衍生物的爆轟性能。通過計算弱鍵的化學(xué)解離能(BDE),得出部分1H-1,2,4三唑衍生物具有較好的熱穩(wěn)定性。

    有機化學(xué);量子化學(xué);密度泛函理論;1H-1,2,4三唑;生成焓;爆轟性能;熱穩(wěn)定性

    引 言

    低感度、高能量密度化合物是含能材料研究領(lǐng)域的重要方向。三唑是含有3個氮原子的五元環(huán)化合物,分子結(jié)構(gòu)中的高氮(50.9%)低碳?xì)浜坎粌H使其具有較高的密度,而且更容易達(dá)到氧平衡,其衍生物有可能作為高能量密度材料[1]。鈍感炸藥5-氨基-3-硝基-1,2,4-三唑(ANTA)和3-硝基-1,2,4-三唑-5-酮(NTO)是1,2,4-三唑重要的衍生物,ANTA可用于火箭推進(jìn)劑,NTO及其鹽類是具有廣泛應(yīng)用前景的鈍感含能材料。

    近年來,三唑衍生物作為含能材料或含能材料中間體備受關(guān)注。Krishna kumar等[2]研究了3,5-二硝基-1,2,4-三唑分子間的相互作用;XU L N等[3]研究了NTO分子間的相互作用;Siker A K等[4]研究了ANTA及其衍生物的熱穩(wěn)定性;朱佳平等[5]計算了DNMT的爆轟性能;魯亞琳[6]研究了1,2,3-三唑及1,2,4-三唑的前線軌道能及熱力學(xué)性質(zhì)。為了深入研究取代基團對1 H-1,2,4-三唑性能的影響,本研究對1H-1,2,4-三唑衍生物的爆轟性能進(jìn)行了系統(tǒng)的計算,以期為1 H-1,2,4-三唑的實驗研究提供參考。

    1 計算方法

    根據(jù)取代基種類、位置及數(shù)量的不同,將1H-1,2,4-三唑衍生物分為A、B、C、D和E五類(如圖1和表1所示)。運用Gaussian03程序[7],采用密度泛函理論(DFT)B3LYP[8-9]方法結(jié)合6-31G(d,p)基組,對1 H-1,2,4-三唑衍生物進(jìn)行分子幾何構(gòu)型優(yōu)化計算。經(jīng)振動頻率分析,確定優(yōu)化構(gòu)型具有勢能面上極小點(無虛頻),得到分子的總能量。

    圖1 1H-1,2,4三唑衍生物的分子結(jié)構(gòu)Fig.1 Molecular structure of 1H-1,2,4-triazole derivatives

    表1 1H-1,2,4-三唑衍生物的取代基Table 1 Substituent of 1H-1,2,4-triazole derivatives

    等鍵反應(yīng)廣泛用于計算化合物的生成焓[10-11]。本研究設(shè)計如下等鍵反應(yīng):

    根據(jù)式(1)所確定的等鍵反應(yīng)計算化合物的氣相生成焓,由于含能材料大部分以固相形式存在,在估算爆轟性能時應(yīng)先計算其固相生成焓。根據(jù)蓋斯定律:ΔHf,solid=ΔHf,gas-ΔHsub,因此,只需計算升華焓,即可得到目標(biāo)化合物的固相生成焓。Politzer[12-14]計算升華焓的經(jīng)驗公式為:

    式中:a=11.166×10-4kJ·mol-1A-4;b=6.900 kJ·mol-1;c=12.404 kJ·mol-1[15]

    運用Monte-Carlo方法由0.001e/Bohr3等電子密度面所包圍的體積求得理論密度,采用公式(3)修正求得其晶體密度:

    式中:α=0.9183;β=0.0028;γ=0.0443[16]

    在求得固相生成焓和晶體密度的基礎(chǔ)上,采用Kamlet-Jacobs公式[17-18]估算了1 H-1,2,4-三唑衍生物的爆速和爆壓,篩選出爆轟性能較好的化合物。

    為了考察化合物的熱穩(wěn)定性,以鍵解離能(BDE)作為衡量化合物熱穩(wěn)定性的一個指標(biāo),采用公式(6)計算1 H-1,2,4-三唑衍生物中較弱化學(xué)鍵的鍵解離能:

    2 結(jié)果與討論

    2.1 生成焓

    表1列出7種相關(guān)化合物的氣相生成焓,其中1 H-1,2,4-三唑(C2H3N3)、CH4、CH3NH2、CH3NO2和CH3ONO2生成焓的實驗值來自文獻(xiàn)[19-21]。由表1可看出,1H-1,2,4-三唑、CH4、CH3NN2、CH3NO2和CH3ONO2的氣相生成焓計算值與實驗值的相對誤差分別為0.39%、3.77%、0.71%、4.32%、8.89%,表明在G2水平下相關(guān)化合物的計算值可信。

    根據(jù)等鍵反應(yīng)式(1)計算1 H-1,2,4-三唑衍生物的氣相生成焓、升華焓以及固相生成焓,結(jié)果見表3。從表3可以看出,1 H-1,2,4-三唑衍生物的生成焓受取代基種類、位置以及數(shù)量的影響。觀察A、B、C系列可見,發(fā)生C-取代時引入-N3基團能大幅度提高衍生物的生成焓,引入-NH2、-NF2或-ONO2基團則導(dǎo)致衍生物的生成焓降低。發(fā)生N-取代時(D系列)引入-N3、-NO2、-NH2、-NF2或-ONO2基團都能提高衍生物生成焓,引入基團能夠提高衍生物生成焓的順序為:-N3>-NO2>-NH2>-NF2>-ONO2。

    表2 相關(guān)化合物的氣相生成焓及相對誤差Table 2 Heats of formation and relative error for the related compounds

    表3 1H-1,2,4-三唑衍生物的總能量、零點能、熱校正和生成焓的計算值Table 2 Calculated values of total energies,zero-point energies,thermal correction,heats of formation for the 1H-1,2,4-triazole derivatives

    從表3看出,發(fā)生C-取代的衍生物,取代基位置對衍生物生成焓的影響并不明顯,而發(fā)生N-取代的衍生物生成焓則大幅提高。

    1 H-1,2,4-三唑衍生物引入基團數(shù)與生成焓的關(guān)系見表4。從表4可以看出,隨著C上的H逐漸被-NH2、-NF2或-ONO2取代,衍生物的生成焓越來越??;而被-N3取代時,衍生物的生成焓隨著-N3基團數(shù)的增多而增大;當(dāng)取代基為-NO2時,隨著-NO2基團數(shù)的增加,衍生物生成焓呈先減小后增大的規(guī)律。對-NH2、-N3和-ONO2三種基團衍生物,比較P-Mono(C3)+P-Mono(C5)和P-Di(C,C)的值以及P-Mono(C3)+P-Mono(C5)+PMono(N)和P-Tri(C,C,N)的值,發(fā)現(xiàn)1 H-1,2,4-三唑引入-NH2,-N3和-ONO2基團后其衍生物的生成焓有較好的基團加和性。而當(dāng)取代基為-NF2基團時由于位阻效應(yīng),導(dǎo)致衍生物生成焓不符合基團加和性。當(dāng)引入-NO2基團時,既有位阻效應(yīng)又能生成弱氫鍵,同時影響衍生物的生成焓。

    表4 1H-1,2,4-三唑引入取代基數(shù)目與生成焓的關(guān)系Table 4 Relationship of heats of formation with the substituents number join to 1H-1,2,4-triazole

    2.2 能量性能

    采用Kamlet-Jacobs公式估算爆速(D)和爆壓(p)。從Kamlet-Jacobs公式可以看出,爆速、爆壓與密度和爆熱有關(guān)。運用Monte-Carlo方法由0. 001e/Bohr3等電子密度面所包圍的體積求得理論密度,通過修正公式(3)求得其晶體密度,根據(jù)得到的固相生成焓計算其爆熱。衍生物的密度、爆熱、氧平衡、爆速和爆壓計算值見表5。

    表5可以顯示,取代基的種類、位置以及基團數(shù)量對衍生物的密度有一定的影響。除A1和D1外,衍生物的密度都大于1 H-1,2,4-三唑的密度,基團對密度增大的貢獻(xiàn)順序為:-NF2>-ONO2>-NO2>-N3>-NH2。同一個基團發(fā)生C-取代比發(fā)生N-取代得到的衍生物的密度大,而在發(fā)生C-取代時取代位置對密度的影響卻不大。衍生物密度隨著取代基數(shù)量的增加而增大,這種增大的趨勢隨著取代基數(shù)量的不斷增大而變緩。

    從表5可以看出,-NF2、-ONO2、-NO2和-N3基團發(fā)生C-取代或N-取代都能提高1H-1,2,4-三唑衍生物的爆熱。而-NH2基團發(fā)生N-取代時1H-1,2,4-三唑衍生物的爆熱增加,但發(fā)生C-取代爆熱降低。

    表5 1H-1,2,4-三唑衍生物的密度、爆熱、氧平衡、爆速和爆壓計算值Table 5 Predicted densities,heats of detonation,oxygen balance,detonation velocities and detonation pressures of 1H-1,2,4-triazole derivatives

    續(xù)表5

    從表5可以看出,密度和爆熱受取代基團類型、取代位置和取代數(shù)量的影響而變化,根據(jù)Kamlet-Jacobs公式,密度和爆熱的變化決定了爆速和爆壓的變化。因此,爆速和爆壓也與取代基的類型、取代位置和取代數(shù)量有關(guān)。-NF2、-NO2、-NO2和-N3基團不但能提高衍生物的密度,同時還能提高其爆熱,因此這些基團的引入可大幅度提高1 H-1,2,-4-三唑衍生物的爆轟性能。

    由表5可以看出,C3,C5,E3和E5的爆轟性能高于RDX,但是只有C3和E3的爆轟性能優(yōu)于HMX。

    2.3 熱穩(wěn)定性

    在計算過程中根據(jù)Mulliken鍵級選擇相對較弱的鍵,由式(6)計算其化學(xué)鍵解離能(BDE)。表6列出1 H-1,2,4-三唑衍生物中弱鍵的解離能和Mulliken鍵級。對于發(fā)生C-取代時(A、B和C系列),含-NH2取代基時環(huán)內(nèi)NI-N2鍵的BDE小于C-R鍵,說明衍生物開環(huán)反應(yīng)是誘導(dǎo)反應(yīng);相反,含-NO2和-N3基團的衍生物環(huán)內(nèi)N1-N2鍵的BDE大于C3-R鍵,說明在衍生物中取代基的脫落是誘導(dǎo)反應(yīng);含-NF2和-ONO2基團的衍生物,N-F和O-NO2鍵的斷裂是解離過程中最容易發(fā)生的。而發(fā)生N-取代(D和E系列),并且只有-ONO2基團引入到N1位置時,O-NO2鍵最弱;當(dāng)-NH2、-N3、-NF2或-NO2引入到N1位置時,N1-R鍵的BDE遠(yuǎn)小于其他鍵,說明1H-1,2,4-三唑衍生物中只要發(fā)生N-取代,化合物中最可能發(fā)生斷裂的鍵是N1-R鍵,而-ONO2基團引入到N1位置時,最容易發(fā)生斷裂的鍵則是ONO2鍵。

    Chung[23]和Su[24]等提出,作為新的高能材料候選物的首要條件是其鍵解離能應(yīng)超過80k J/mol。從表5數(shù)據(jù)可知,所有目標(biāo)化合物的鍵解離能高于80 kJ/mol,但在實際應(yīng)用中有待進(jìn)一步實驗驗證。

    表6 1H-1,2,4-三唑衍生物中弱鍵的解離能和Mulliken鍵級Table 6 Bond dissociation energies and Mulliken bond orders of the relatively weak bonds in the 1H-1,2,4-triazole derivatives

    續(xù)表6

    3 結(jié) 論

    (1)運用DFT-B3LYP方法研究了1H-1,2,4-三唑衍生物的生成焓、爆轟性能和熱穩(wěn)定性。取代基的種類和位置影響衍生物的生成焓,引入-N3基團能顯著增加衍生物的生成焓。對于同類取代基發(fā)生N-取代衍生物生成焓值高于發(fā)生C-取代衍生物。

    (2)-NF2、-ONO2、-NO2和-N3基團不僅能提高衍生物的密度,還能提高其爆熱,有效提高1H-1,2,4-三唑衍生物的爆轟性能。所有衍生物的熱穩(wěn)定性較好。綜合考慮爆轟性能和熱穩(wěn)定性,C3、C5、E3和E5有可能成為潛在的高能量密度材料。

    [1] 陽世清,岳守體.國外四嗪四唑類高氮含能材料研究進(jìn)展[J].含能材料,2003,11(4):231-235.

    YANG Shi-qing,YUE Shou-ti.Progress in high-nitrogen energetic materials derails derived from tetrazine and tetrazole[J].Chinese Journal of Energetic Materials,2003,11(4):231-235.

    [2] Krishna K V,Keresztury G,Sundius T,et al.Hydrogen bonding and molecular vibrations of 3,5-diamino-1,2,4-triazole[J].Spectrochim Acta A,2005,61:261-267.

    [3] Xu L N,Xiao H M,F(xiàn)ang G Y,et al.Theoretical study on intermolecular interactions of 3-nitro-1,2, 4triazole-5-one-dimers[J].Acta Chimica Sinica,2005,63(12):1062-1068.

    [4] Sikder A K,Geetha M,Sarwade D B,et al.Studies on characterization and thermal behaviour of 3-amino-5-nitro-1,2,4-triazole and its derivatives[J].Hazard,2001,82(1):1-12.

    [5] 朱佳平,劉勝楠,李永祥.幾種唑類炸藥分子的結(jié)構(gòu)和爆轟性能的密度泛函理論研究[J].中北大學(xué)學(xué)報(自然科學(xué)版),2011,2:200-206.

    ZHU Jia-ping,LIU Sheng-nan,LI Yong-xiang. Structures and detonation performance for several explosive molecule of thiiazoles by density functional theory[J].Chinese Journal of North University of China(Natural Science Edition),2011,2:200-206.

    [6] Lu Y L,Gong X D,Ju X H.structures and properties of 1,2,3-triazoles and 1,2,4-triazoles[J].Chinese J Struct Chem,2006,25(5):582-588.

    [7] Frisch M J,Trucks G W,Schlegel H B.et al.Gaussian 03,Revision B.03.[M].Pittsburgh:Gaussian Inc,2003.

    [8] Lee C,Yang W,Parr R G.Development of the collsalvetti correlation-energy formula into a functional of the electron density[J].Phys Rev B,1988,37:785-789.

    [9] Becke A D.Density-functional thermochemistry.Ⅲ. the role of exact exchange[J].Chem Phys,1993,98:5648-5652.

    [10]Zhang X W,Zhu W H,Xiao H M.Comparative theoretical studies of energetic substituted carbon-and ni-trogen-bridged difurazans[J].Phys Chem A,2010,114:603-612.

    [11]Ju X H,Wang X,Bei F L.Substituent effects on heats of formation,group interactions,and detonation properties of polyazidocubanes[J].Comput Chem,2005,26(12):1263-1269.

    [12]Politzer P,Lane P,Murray J S.Computational characterization of a potential energetic compound:1,3,5,7-Tetranitro-2,4,6,8-tetraazacubane[J].Cent Eur J Energ Mater,2011,8:39-52.

    [13]Politzer P,Murray J S.Some perspectives on estimating detonation properties of C,H,N,O compounds[J].Cent Eur J Energ Mater,2011,8:209-220.

    [14]Politzer P,Murray J S,Grice M E,et al.Calculation of heats of sublimation and solid phase heats of formation[J].Mol Phys,1997,91:923-928.

    [15]Byrd E F C,Rice B M.Improved prediction of heats of formation of energetic materials using quantum chemical calculations[J].Phys Chem A,2006,110:1005-1013.

    [16]Politzer P,Martinez J,Murray JS.et al.An electrostatic interaction correction for improved crystal density predictions[J].Mol Phys,2009,107:2095-2101.

    [17]Kamlet M J,Jacobs S J.Chemistry of detonation.I.A simple method for calculating detonation properties of C,H,N,O explosives[J].Chem Phys,1968,48:23-35.

    [18]Kamlet M J,Ablard J E.Chemistry of detonations. II.buffered equilibria[J].Chem Phys,1968,48:36-42.

    [19]Dean J A.Lange′s Handbook of Chemistry[M].15th Ed.New York:McGraw-Hill,1999.

    [20]David R L.Handbook of chemistry and physics[M]. 84th edn.Boca Raton:CRC Press,2003-2004.

    [21]Afeefy H Y,Liebman J F,Stein S E.Neutral Thermochemical Data[M].Gaithersburg:National Institute of Standards and Technology,2000.

    [22]Talawar M B,Sivabalan R,Mukundan T,et al.Environmentally compatible next generation green energetic materials(GEMs)[J].Hazard Mater,2009,161:589-607.

    [23]Chung G S,Schimidt M W,Gordon M S.An ab initio study of potential energy surfaces for N8 isomers[J]. Phys Chem A,2000,104:5647-5650.

    [24]Su X F,Cheng X L,Meng C M.et al.Quantum chemical study on nitroimidazole,polynitroimidazole and their methyl derivatives[J].Hazard Mater,2009,161(1):551-558.

    Theoretical Calculation on Detonation Performance of 1H-1,2,4-Triazole Derivatives

    MA Zhong-liang1,ZHOU Hua1,WANG Jian-long1,LI Jing-man2
    (1.School of Chemical Engineering and Environment,North University of China,Taiyuan 030051,China;2.Sichuan Tongda Chemical Industry Co.,Ltd,Dazhou Sichuan 635000,China)

    The geometries of 1H-1,2,4-triazole derivatives were obtained by density functional theory at DFTB3LYP/6-31G(d,p)level.The heats of formation(HOF)were predicted by designed isodesmic reactions.The average molar volume(V)and theoretical density(ρ)were estimated based on 0.001 e·bohr-3density space.The detonation velocity(D)and pressure(p)of 1H-1,2,4-triazole derivatives were estimated by using the Kamlet-Jacbos equation.It was found that N-substituent in 1H-1,2,4-triazole can improve the HOF values of the derivatives,and-N3is an effective structural unit for improving the HOF values of the derivatives.The-NF2or-ONO2group is very useful for enhancing the detonation performance of these derivatives.The analysis of the bond dissociation energies for several relatively weak bonds suggest that most of the derivatives have good thermal stability.

    organic chemistry;quantum chemistry;density functional theory;1H-1,2,4-triazole derivatives;heat of formation;detonation performance;thermal stability

    TJ55;O<641 文獻(xiàn)標(biāo)志碼:A class="emphasis_bold">641 文獻(xiàn)標(biāo)志碼:A 文章編號:1007-7812(2014)04-0038-07641 文獻(xiàn)標(biāo)志碼:A

    1007-7812(2014)04-0038-07

    A 文章編號:1007-7812(2014)04-0038-07

    2014-03-10;

    2014-05-30

    馬忠亮(1967-),男,教授,博導(dǎo),從事新型發(fā)射藥及裝藥研究。

    猜你喜歡
    三唑熱穩(wěn)定性基團
    R基團篩選技術(shù)用于HDACIs的分子設(shè)計
    芳烴ArCOR的構(gòu)象分析和基團對親電取代反應(yīng)的定位作用
    PVC用酪氨酸鑭的合成、復(fù)配及熱穩(wěn)定性能研究
    中國塑料(2016年7期)2016-04-16 05:25:52
    不同濃度三唑錫懸浮劑防治效果研究
    中國果菜(2016年9期)2016-03-01 01:28:41
    提高有機過氧化物熱穩(wěn)定性的方法
    可聚合松香衍生物的合成、表征和熱穩(wěn)定性?
    內(nèi)含雙二氯均三嗪基團的真絲織物抗皺劑的合成
    三組分反應(yīng)高效合成1,2,4-三唑烷類化合物
    對羥基安息香醛苯甲酰腙的合成、表征及熱穩(wěn)定性
    1,1′-二(硝氧甲基)-3,3′-二硝基-5,5′-聯(lián)-1,2,4-三唑的合成及性能
    久久这里只有精品19| 国产亚洲最大av| 欧美成人精品欧美一级黄| 老司机影院成人| 一级毛片我不卡| 婷婷色综合www| av电影中文网址| 国产老妇伦熟女老妇高清| 成年人午夜在线观看视频| 一本久久精品| 青春草国产在线视频| 亚洲国产欧美日韩在线播放| 韩国精品一区二区三区| 女人被躁到高潮嗷嗷叫费观| 黄色视频在线播放观看不卡| 国产精品免费视频内射| 国产爽快片一区二区三区| 国产又爽黄色视频| 国产伦理片在线播放av一区| 成年女人在线观看亚洲视频| 欧美97在线视频| 国产深夜福利视频在线观看| 午夜影院在线不卡| 人妻 亚洲 视频| 巨乳人妻的诱惑在线观看| 亚洲精品第二区| 精品少妇一区二区三区视频日本电影 | 久久精品国产亚洲av涩爱| 亚洲国产看品久久| 999久久久国产精品视频| 一区二区三区激情视频| 亚洲国产精品成人久久小说| 久久 成人 亚洲| 国产av国产精品国产| 9191精品国产免费久久| 在现免费观看毛片| 99热全是精品| 老熟女久久久| 黄色 视频免费看| 亚洲综合精品二区| videossex国产| 咕卡用的链子| 99久久精品国产国产毛片| 久久ye,这里只有精品| 久久综合国产亚洲精品| 黄片播放在线免费| av国产精品久久久久影院| 日本午夜av视频| 午夜av观看不卡| 亚洲成色77777| 亚洲精华国产精华液的使用体验| 美女中出高潮动态图| 精品少妇一区二区三区视频日本电影 | 美女午夜性视频免费| 老司机影院成人| 99久久精品国产国产毛片| 夜夜骑夜夜射夜夜干| 亚洲情色 制服丝袜| 啦啦啦在线观看免费高清www| 午夜av观看不卡| 亚洲精品久久久久久婷婷小说| 在线天堂最新版资源| 久久亚洲国产成人精品v| 日韩精品免费视频一区二区三区| 老熟女久久久| 国产日韩欧美视频二区| 欧美日韩精品网址| 国产精品99久久99久久久不卡 | 青春草亚洲视频在线观看| av网站在线播放免费| 亚洲av成人精品一二三区| 亚洲精品日韩在线中文字幕| 在线观看免费日韩欧美大片| 69精品国产乱码久久久| 在线观看三级黄色| 成人影院久久| 色播在线永久视频| 国产 精品1| 久久久久久人妻| 亚洲精品久久成人aⅴ小说| 黄色一级大片看看| 香蕉丝袜av| 少妇熟女欧美另类| 美国免费a级毛片| 欧美另类一区| 黑人猛操日本美女一级片| 啦啦啦在线免费观看视频4| 精品少妇内射三级| 夫妻午夜视频| 久久久国产一区二区| 国产片内射在线| av国产久精品久网站免费入址| 精品少妇一区二区三区视频日本电影 | 在线天堂最新版资源| 国产女主播在线喷水免费视频网站| 亚洲熟女精品中文字幕| 午夜久久久在线观看| 丝袜美足系列| 性色av一级| 国产日韩欧美亚洲二区| 亚洲在久久综合| 免费观看av网站的网址| 永久免费av网站大全| 大香蕉久久成人网| 亚洲伊人久久精品综合| 在线观看国产h片| 性色avwww在线观看| 热99国产精品久久久久久7| 99九九在线精品视频| 久久午夜福利片| 国产成人a∨麻豆精品| 一本久久精品| 精品少妇一区二区三区视频日本电影 | 久久女婷五月综合色啪小说| 亚洲一码二码三码区别大吗| 免费观看a级毛片全部| www.精华液| 一本色道久久久久久精品综合| 久久人人爽av亚洲精品天堂| 九色亚洲精品在线播放| 亚洲第一av免费看| 国产在线免费精品| 免费观看av网站的网址| 国产精品无大码| 一区二区av电影网| 爱豆传媒免费全集在线观看| 2018国产大陆天天弄谢| 亚洲欧洲精品一区二区精品久久久 | 国产在视频线精品| 色视频在线一区二区三区| 一二三四中文在线观看免费高清| 亚洲,一卡二卡三卡| 高清不卡的av网站| 国产在视频线精品| 麻豆精品久久久久久蜜桃| 伊人久久大香线蕉亚洲五| 成人午夜精彩视频在线观看| 波野结衣二区三区在线| 欧美精品人与动牲交sv欧美| 一本—道久久a久久精品蜜桃钙片| xxxhd国产人妻xxx| 欧美成人午夜免费资源| 一边摸一边做爽爽视频免费| 日韩av不卡免费在线播放| 精品人妻在线不人妻| 精品少妇一区二区三区视频日本电影 | 九草在线视频观看| 久久精品aⅴ一区二区三区四区 | 国精品久久久久久国模美| 亚洲国产欧美在线一区| av免费在线看不卡| 久久精品亚洲av国产电影网| 天天躁夜夜躁狠狠久久av| 免费观看在线日韩| 成年美女黄网站色视频大全免费| av在线app专区| 亚洲三级黄色毛片| 蜜桃国产av成人99| 免费观看无遮挡的男女| 制服人妻中文乱码| √禁漫天堂资源中文www| 中文字幕色久视频| 少妇人妻久久综合中文| 国产黄频视频在线观看| 日韩人妻精品一区2区三区| 久久久久久久久久人人人人人人| 性高湖久久久久久久久免费观看| 尾随美女入室| 日韩一区二区视频免费看| 久久精品国产a三级三级三级| 妹子高潮喷水视频| 亚洲五月色婷婷综合| 我要看黄色一级片免费的| 亚洲欧美一区二区三区久久| 啦啦啦在线免费观看视频4| 青青草视频在线视频观看| 男女边吃奶边做爰视频| 天天躁夜夜躁狠狠久久av| 好男人视频免费观看在线| 日韩av免费高清视频| 欧美人与性动交α欧美精品济南到 | 尾随美女入室| 久久精品国产亚洲av涩爱| 男女午夜视频在线观看| 狠狠精品人妻久久久久久综合| 久久久久国产精品人妻一区二区| 久久久精品区二区三区| 国产野战对白在线观看| 一区二区日韩欧美中文字幕| 国产野战对白在线观看| 精品人妻一区二区三区麻豆| 精品一区二区三卡| 久久亚洲国产成人精品v| www.av在线官网国产| 午夜免费观看性视频| 三级国产精品片| 我要看黄色一级片免费的| 亚洲四区av| 日韩一区二区三区影片| 黑人巨大精品欧美一区二区蜜桃| 国产免费又黄又爽又色| 国产色婷婷99| 如何舔出高潮| 精品国产一区二区三区四区第35| 亚洲精品一区蜜桃| 久久精品国产a三级三级三级| 成年女人在线观看亚洲视频| 18+在线观看网站| 人人妻人人澡人人爽人人夜夜| videosex国产| 精品视频人人做人人爽| 婷婷色av中文字幕| 久久久久国产精品人妻一区二区| 欧美日韩亚洲高清精品| 日韩视频在线欧美| 欧美日韩精品网址| 午夜福利视频精品| 91精品国产国语对白视频| 99九九在线精品视频| 春色校园在线视频观看| 国产精品嫩草影院av在线观看| 王馨瑶露胸无遮挡在线观看| 天堂8中文在线网| 国语对白做爰xxxⅹ性视频网站| av在线播放精品| 久久人人爽av亚洲精品天堂| 亚洲国产精品一区三区| 亚洲欧洲精品一区二区精品久久久 | 国产一区二区 视频在线| 精品国产国语对白av| 三级国产精品片| tube8黄色片| 国产成人午夜福利电影在线观看| 国产精品免费视频内射| 亚洲欧美一区二区三区久久| videossex国产| 一区二区三区激情视频| 波多野结衣av一区二区av| 不卡视频在线观看欧美| 波野结衣二区三区在线| 国产男女超爽视频在线观看| 母亲3免费完整高清在线观看 | 免费av中文字幕在线| 啦啦啦视频在线资源免费观看| av国产精品久久久久影院| 日韩制服骚丝袜av| 又大又黄又爽视频免费| 91aial.com中文字幕在线观看| 日本vs欧美在线观看视频| 欧美精品av麻豆av| av不卡在线播放| 少妇熟女欧美另类| 建设人人有责人人尽责人人享有的| 伊人久久大香线蕉亚洲五| 宅男免费午夜| 精品酒店卫生间| 天堂中文最新版在线下载| 香蕉国产在线看| 国产男女内射视频| 午夜老司机福利剧场| 在线观看www视频免费| 国产亚洲一区二区精品| 99久久综合免费| 欧美日韩av久久| 丁香六月天网| av有码第一页| 一区在线观看完整版| 久久久久久伊人网av| 久久这里只有精品19| av网站在线播放免费| 黑人猛操日本美女一级片| 综合色丁香网| av国产久精品久网站免费入址| 高清欧美精品videossex| 免费在线观看视频国产中文字幕亚洲 | 大陆偷拍与自拍| 久久99一区二区三区| 叶爱在线成人免费视频播放| 久久精品国产亚洲av高清一级| 亚洲精品自拍成人| 成人国产麻豆网| 黑人猛操日本美女一级片| 免费观看无遮挡的男女| 五月天丁香电影| 少妇 在线观看| 久久精品国产亚洲av天美| 久久精品aⅴ一区二区三区四区 | 久久久久久久久久久久大奶| 亚洲美女黄色视频免费看| 菩萨蛮人人尽说江南好唐韦庄| av一本久久久久| 久久99热这里只频精品6学生| av国产久精品久网站免费入址| 亚洲精品第二区| av在线老鸭窝| 人人妻人人澡人人爽人人夜夜| 两性夫妻黄色片| 亚洲av男天堂| 在线观看国产h片| 国产成人免费无遮挡视频| 美女福利国产在线| 亚洲少妇的诱惑av| 看十八女毛片水多多多| 午夜av观看不卡| 两个人看的免费小视频| 亚洲精品国产av成人精品| 欧美成人精品欧美一级黄| 国产伦理片在线播放av一区| 久久99一区二区三区| 国产精品亚洲av一区麻豆 | 日本色播在线视频| 看非洲黑人一级黄片| 亚洲精品乱久久久久久| 制服人妻中文乱码| 激情视频va一区二区三区| 欧美激情 高清一区二区三区| 久久免费观看电影| 欧美另类一区| 最近手机中文字幕大全| 久热久热在线精品观看| 久久久欧美国产精品| 视频在线观看一区二区三区| 日本色播在线视频| 亚洲国产欧美在线一区| 在现免费观看毛片| 97在线人人人人妻| 人体艺术视频欧美日本| 各种免费的搞黄视频| 欧美人与性动交α欧美软件| 自线自在国产av| 香蕉精品网在线| 久久久久久人人人人人| 亚洲视频免费观看视频| 日韩av在线免费看完整版不卡| √禁漫天堂资源中文www| 一级a爱视频在线免费观看| 在线观看人妻少妇| 久久鲁丝午夜福利片| 视频在线观看一区二区三区| 国产男人的电影天堂91| 高清欧美精品videossex| 久久精品国产亚洲av天美| 久久精品国产亚洲av涩爱| 久久久亚洲精品成人影院| 桃花免费在线播放| 免费av中文字幕在线| 国产成人免费无遮挡视频| 久久精品国产a三级三级三级| 日韩熟女老妇一区二区性免费视频| 高清黄色对白视频在线免费看| 91国产中文字幕| 大片电影免费在线观看免费| 日韩免费高清中文字幕av| 久久久久久久国产电影| 少妇的逼水好多| 亚洲国产精品一区二区三区在线| 桃花免费在线播放| 黄色毛片三级朝国网站| 蜜桃国产av成人99| 免费在线观看完整版高清| 交换朋友夫妻互换小说| 亚洲av.av天堂| 99久久综合免费| 天天躁夜夜躁狠狠躁躁| 老汉色av国产亚洲站长工具| 国产精品偷伦视频观看了| 免费久久久久久久精品成人欧美视频| 夜夜骑夜夜射夜夜干| 韩国精品一区二区三区| 大香蕉久久网| 亚洲天堂av无毛| 丰满迷人的少妇在线观看| 国产黄色视频一区二区在线观看| 在线观看一区二区三区激情| 99精国产麻豆久久婷婷| 日韩中字成人| 久久韩国三级中文字幕| 久久久久久久国产电影| 亚洲国产日韩一区二区| 丝袜人妻中文字幕| 亚洲综合色网址| 欧美人与性动交α欧美软件| 日韩精品免费视频一区二区三区| 国产av码专区亚洲av| 亚洲精品乱久久久久久| 国产片特级美女逼逼视频| 晚上一个人看的免费电影| 欧美黄色片欧美黄色片| 国产高清不卡午夜福利| 一区在线观看完整版| 免费观看性生交大片5| 美女视频免费永久观看网站| 久久国产亚洲av麻豆专区| 欧美日韩精品网址| 在线观看人妻少妇| a级毛片在线看网站| 91在线精品国自产拍蜜月| 久久久国产欧美日韩av| 成年美女黄网站色视频大全免费| 亚洲内射少妇av| 一二三四在线观看免费中文在| 国产伦理片在线播放av一区| 国产极品粉嫩免费观看在线| 中文字幕精品免费在线观看视频| 亚洲第一青青草原| 亚洲欧洲国产日韩| 五月天丁香电影| 80岁老熟妇乱子伦牲交| 久久久久精品人妻al黑| 汤姆久久久久久久影院中文字幕| 搡老乐熟女国产| 9热在线视频观看99| 黄片小视频在线播放| 国产 一区精品| 咕卡用的链子| 两个人免费观看高清视频| 久久精品国产自在天天线| av卡一久久| 亚洲图色成人| 肉色欧美久久久久久久蜜桃| 国产精品免费大片| 亚洲三区欧美一区| 99国产精品免费福利视频| 桃花免费在线播放| 韩国av在线不卡| 99国产综合亚洲精品| 成人影院久久| 高清在线视频一区二区三区| 纵有疾风起免费观看全集完整版| 久久影院123| 性色av一级| 美女中出高潮动态图| 国产黄频视频在线观看| 看免费av毛片| 国产淫语在线视频| 日本av免费视频播放| 春色校园在线视频观看| kizo精华| 男的添女的下面高潮视频| 美女xxoo啪啪120秒动态图| 欧美日韩亚洲高清精品| 久久国产精品大桥未久av| 亚洲av福利一区| 人妻系列 视频| 黄色一级大片看看| 久久久久久伊人网av| 国产成人aa在线观看| 亚洲精品美女久久av网站| 色吧在线观看| 在线免费观看不下载黄p国产| www.av在线官网国产| 一区二区av电影网| 亚洲av免费高清在线观看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 欧美激情高清一区二区三区 | 大香蕉久久成人网| 我的亚洲天堂| 日日啪夜夜爽| 边亲边吃奶的免费视频| 国产日韩一区二区三区精品不卡| 久久久久久人妻| 中文天堂在线官网| 成人影院久久| 久久精品国产鲁丝片午夜精品| 熟女电影av网| 一本—道久久a久久精品蜜桃钙片| 2021少妇久久久久久久久久久| 成年动漫av网址| 免费高清在线观看视频在线观看| 国产淫语在线视频| 老汉色av国产亚洲站长工具| 一级毛片 在线播放| 亚洲三区欧美一区| 一级毛片黄色毛片免费观看视频| 最近中文字幕高清免费大全6| 欧美在线黄色| 超色免费av| 久久午夜福利片| 成人午夜精彩视频在线观看| 精品一区二区三区四区五区乱码 | 欧美精品高潮呻吟av久久| 汤姆久久久久久久影院中文字幕| 亚洲欧美一区二区三区久久| 热99久久久久精品小说推荐| 亚洲欧洲国产日韩| 免费黄色在线免费观看| 中文字幕人妻丝袜制服| av在线播放精品| 欧美变态另类bdsm刘玥| 日韩一卡2卡3卡4卡2021年| 午夜激情av网站| 成年人免费黄色播放视频| 国产男人的电影天堂91| 久久久久久人人人人人| 国产成人精品久久久久久| 久久青草综合色| 久久这里有精品视频免费| 亚洲欧洲国产日韩| 在线免费观看不下载黄p国产| 国产男女超爽视频在线观看| av一本久久久久| 老女人水多毛片| 国产精品偷伦视频观看了| 中文字幕精品免费在线观看视频| 亚洲精品美女久久久久99蜜臀 | 欧美成人午夜精品| 90打野战视频偷拍视频| 精品一区在线观看国产| 亚洲国产成人一精品久久久| 在线看a的网站| 亚洲精品久久午夜乱码| 亚洲三级黄色毛片| 看非洲黑人一级黄片| 国产精品久久久久久精品古装| 最近最新中文字幕大全免费视频 | 久久精品熟女亚洲av麻豆精品| 观看av在线不卡| 久久久久久久精品精品| 欧美日韩亚洲高清精品| 一级片免费观看大全| 精品国产国语对白av| 在线精品无人区一区二区三| 亚洲av福利一区| 女的被弄到高潮叫床怎么办| 一级片免费观看大全| 婷婷色综合www| 女人久久www免费人成看片| 青春草亚洲视频在线观看| 80岁老熟妇乱子伦牲交| 日韩中文字幕视频在线看片| 午夜精品国产一区二区电影| 久久免费观看电影| 日韩免费高清中文字幕av| 看免费成人av毛片| 国产亚洲午夜精品一区二区久久| 老熟女久久久| 日本wwww免费看| 纯流量卡能插随身wifi吗| 久久99蜜桃精品久久| 亚洲五月色婷婷综合| 母亲3免费完整高清在线观看 | 久热久热在线精品观看| 91午夜精品亚洲一区二区三区| 久久精品国产综合久久久| 国产精品.久久久| 日日摸夜夜添夜夜爱| 捣出白浆h1v1| 丰满迷人的少妇在线观看| 免费看不卡的av| 91久久精品国产一区二区三区| 午夜福利乱码中文字幕| 爱豆传媒免费全集在线观看| 国产在线一区二区三区精| 国产 一区精品| 亚洲欧美日韩另类电影网站| 免费女性裸体啪啪无遮挡网站| 中文字幕av电影在线播放| 人人妻人人爽人人添夜夜欢视频| 国产男女内射视频| 少妇的逼水好多| 在线观看美女被高潮喷水网站| 国产精品二区激情视频| 久久ye,这里只有精品| 80岁老熟妇乱子伦牲交| 成人影院久久| 国产不卡av网站在线观看| 亚洲精品国产一区二区精华液| tube8黄色片| 又粗又硬又长又爽又黄的视频| 在线观看一区二区三区激情| 美女国产视频在线观看| 亚洲图色成人| 哪个播放器可以免费观看大片| 性色av一级| 青草久久国产| 亚洲精品第二区| av网站免费在线观看视频| 超碰成人久久| 免费观看av网站的网址| 夜夜骑夜夜射夜夜干| 永久网站在线| 欧美国产精品va在线观看不卡| 国产成人精品久久久久久| www日本在线高清视频| 成人毛片60女人毛片免费| 美女国产高潮福利片在线看| 满18在线观看网站| 又黄又粗又硬又大视频| 午夜福利在线免费观看网站| 亚洲精品av麻豆狂野| 亚洲,欧美,日韩| 青春草亚洲视频在线观看| 免费在线观看完整版高清| 国产综合精华液| 精品人妻熟女毛片av久久网站| 免费观看性生交大片5| 亚洲国产日韩一区二区| 精品少妇久久久久久888优播| 久久av网站| 久久午夜综合久久蜜桃| 蜜桃在线观看..| 精品第一国产精品| 亚洲成av片中文字幕在线观看 | 大片免费播放器 马上看| videosex国产| 欧美日韩国产mv在线观看视频| 人人妻人人澡人人爽人人夜夜| 精品卡一卡二卡四卡免费| 看十八女毛片水多多多| 亚洲婷婷狠狠爱综合网| 亚洲成人av在线免费| 亚洲欧美清纯卡通| 免费观看在线日韩| 在现免费观看毛片| 最近最新中文字幕大全免费视频 |