• <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-三唑的合成及性能
    中文字幕高清在线视频| 99热网站在线观看| 99久久久亚洲精品蜜臀av| 91狼人影院| 在线播放无遮挡| 99久久成人亚洲精品观看| 国产久久久一区二区三区| 亚洲欧美精品综合久久99| 特大巨黑吊av在线直播| 美女 人体艺术 gogo| 免费观看的影片在线观看| 免费观看的影片在线观看| a在线观看视频网站| 俺也久久电影网| 桃红色精品国产亚洲av| 日本免费a在线| 国产精品三级大全| 亚洲内射少妇av| 国产黄a三级三级三级人| 精品乱码久久久久久99久播| 亚洲av成人av| 欧美精品国产亚洲| 亚洲自偷自拍三级| 亚洲综合色惰| 丰满的人妻完整版| 亚洲美女视频黄频| 国产一区二区在线观看日韩| 国产乱人视频| 大型黄色视频在线免费观看| 久久热精品热| 国产爱豆传媒在线观看| 国产高清不卡午夜福利| 女同久久另类99精品国产91| 亚洲av二区三区四区| 亚洲国产精品合色在线| 老熟妇乱子伦视频在线观看| 九九在线视频观看精品| 网址你懂的国产日韩在线| 亚洲最大成人中文| 尾随美女入室| 亚洲国产精品合色在线| 精品福利观看| 婷婷精品国产亚洲av| 一个人看的www免费观看视频| 国语自产精品视频在线第100页| 欧美日韩瑟瑟在线播放| 人妻制服诱惑在线中文字幕| 天天一区二区日本电影三级| 天天一区二区日本电影三级| 亚洲色图av天堂| 成年版毛片免费区| 日韩中文字幕欧美一区二区| 男女啪啪激烈高潮av片| 男人舔女人下体高潮全视频| 淫秽高清视频在线观看| 久久久久国内视频| 性色avwww在线观看| 最近在线观看免费完整版| 美女cb高潮喷水在线观看| 人妻少妇偷人精品九色| 国产精品一区www在线观看 | 乱系列少妇在线播放| 国产精品一区www在线观看 | 精品一区二区三区人妻视频| 三级国产精品欧美在线观看| 一区二区三区高清视频在线| 亚洲专区中文字幕在线| 国产精品综合久久久久久久免费| 亚洲一区二区三区色噜噜| av在线蜜桃| 国产精品无大码| 国产高清视频在线观看网站| 国产一区二区在线av高清观看| 狂野欧美激情性xxxx在线观看| 观看美女的网站| 最后的刺客免费高清国语| 日本一二三区视频观看| 一本一本综合久久| 久久久色成人| 午夜a级毛片| 欧美日韩综合久久久久久 | 窝窝影院91人妻| 午夜亚洲福利在线播放| 啦啦啦观看免费观看视频高清| 99国产精品一区二区蜜桃av| 国产黄色小视频在线观看| 国产精华一区二区三区| 桃色一区二区三区在线观看| 高清在线国产一区| 国产精品电影一区二区三区| 人人妻人人看人人澡| 亚洲精华国产精华液的使用体验 | 国内毛片毛片毛片毛片毛片| 最近在线观看免费完整版| 色尼玛亚洲综合影院| 国产精品99久久久久久久久| 亚洲国产精品合色在线| 国产高清不卡午夜福利| 三级国产精品欧美在线观看| 一区二区三区高清视频在线| 22中文网久久字幕| 国产精品久久久久久亚洲av鲁大| 久久亚洲精品不卡| 久久久久久久久久黄片| 丰满乱子伦码专区| 国产精品自产拍在线观看55亚洲| 国产精品女同一区二区软件 | 欧美最黄视频在线播放免费| 欧美日本视频| 国产白丝娇喘喷水9色精品| a级毛片a级免费在线| 国产一区二区三区av在线 | 欧美极品一区二区三区四区| 乱码一卡2卡4卡精品| 亚洲在线观看片| 日日摸夜夜添夜夜添小说| 中国美白少妇内射xxxbb| 亚洲精品一卡2卡三卡4卡5卡| 国产极品精品免费视频能看的| 国产三级中文精品| 午夜福利在线观看吧| 色视频www国产| 日日夜夜操网爽| 成人国产综合亚洲| 日韩中文字幕欧美一区二区| 1000部很黄的大片| 最好的美女福利视频网| 久久久成人免费电影| 亚洲内射少妇av| 久久久久久久精品吃奶| 国内精品久久久久久久电影| 亚洲人成网站在线播| 大又大粗又爽又黄少妇毛片口| 少妇的逼水好多| 国产欧美日韩一区二区精品| 美女 人体艺术 gogo| 国产精品一区二区免费欧美| 在线国产一区二区在线| 熟女电影av网| 国产主播在线观看一区二区| 精品午夜福利在线看| 亚洲av一区综合| 美女cb高潮喷水在线观看| 麻豆av噜噜一区二区三区| 亚洲人成网站在线播| 成人美女网站在线观看视频| 精品福利观看| 91av网一区二区| 少妇的逼好多水| 久久久久九九精品影院| 国产v大片淫在线免费观看| 国产探花极品一区二区| 国产高清视频在线播放一区| 亚洲成人精品中文字幕电影| 久久婷婷人人爽人人干人人爱| 国产一区二区在线av高清观看| 日本 欧美在线| 97超级碰碰碰精品色视频在线观看| 天堂网av新在线| 一级黄色大片毛片| 国产av在哪里看| 免费高清视频大片| 日本a在线网址| 成人鲁丝片一二三区免费| 真实男女啪啪啪动态图| 成人高潮视频无遮挡免费网站| 神马国产精品三级电影在线观看| 久久中文看片网| 露出奶头的视频| 欧美色视频一区免费| 在线看三级毛片| 人妻丰满熟妇av一区二区三区| 亚洲av五月六月丁香网| 免费看a级黄色片| 欧美精品国产亚洲| 亚洲国产精品成人综合色| 亚洲欧美日韩高清专用| 亚洲,欧美,日韩| 全区人妻精品视频| 男女下面进入的视频免费午夜| 欧美日韩中文字幕国产精品一区二区三区| 成人无遮挡网站| 国产精品伦人一区二区| 麻豆一二三区av精品| 日本黄色视频三级网站网址| 男人舔女人下体高潮全视频| 黄色视频,在线免费观看| 精品一区二区三区视频在线| 日本在线视频免费播放| 成人永久免费在线观看视频| ponron亚洲| 免费不卡的大黄色大毛片视频在线观看 | 国产一区二区在线观看日韩| 舔av片在线| 欧美色视频一区免费| 精品久久久久久久久久久久久| 不卡视频在线观看欧美| 又粗又爽又猛毛片免费看| 波多野结衣高清作品| 最近最新免费中文字幕在线| 97超级碰碰碰精品色视频在线观看| 真人做人爱边吃奶动态| 国内精品美女久久久久久| 亚洲美女视频黄频| 69av精品久久久久久| 麻豆精品久久久久久蜜桃| 精品一区二区三区av网在线观看| 黄色日韩在线| 黄片wwwwww| 天堂动漫精品| 欧美zozozo另类| 桃红色精品国产亚洲av| 禁无遮挡网站| 婷婷丁香在线五月| 色5月婷婷丁香| h日本视频在线播放| 999久久久精品免费观看国产| 男女啪啪激烈高潮av片| 亚洲精品日韩av片在线观看| 琪琪午夜伦伦电影理论片6080| 男插女下体视频免费在线播放| 精品久久久久久久久久久久久| 色综合色国产| 91午夜精品亚洲一区二区三区 | 欧美日韩亚洲国产一区二区在线观看| 波多野结衣巨乳人妻| 亚洲一区二区三区色噜噜| 日韩欧美三级三区| 亚洲五月天丁香| 男人狂女人下面高潮的视频| 美女被艹到高潮喷水动态| 中文亚洲av片在线观看爽| 18禁黄网站禁片免费观看直播| 乱人视频在线观看| 人妻少妇偷人精品九色| 久久久久久久亚洲中文字幕| 精品不卡国产一区二区三区| 国产高清激情床上av| 岛国在线免费视频观看| 国产亚洲欧美98| 搡老熟女国产l中国老女人| 亚洲av中文av极速乱 | av黄色大香蕉| 亚洲国产精品合色在线| 亚洲成人精品中文字幕电影| 国国产精品蜜臀av免费| 亚洲欧美日韩高清专用| 国产aⅴ精品一区二区三区波| 国产精品国产高清国产av| 亚洲成人精品中文字幕电影| a级毛片a级免费在线| 国产精品一区www在线观看 | 欧美黑人欧美精品刺激| 亚洲精品久久国产高清桃花| 日本三级黄在线观看| 亚洲va日本ⅴa欧美va伊人久久| 中文字幕av成人在线电影| 波多野结衣高清作品| 亚洲欧美清纯卡通| 亚洲第一电影网av| h日本视频在线播放| 亚洲美女视频黄频| 久久久久久久精品吃奶| 久久精品国产99精品国产亚洲性色| 国产午夜精品论理片| 22中文网久久字幕| 国产av不卡久久| 好男人在线观看高清免费视频| 天堂影院成人在线观看| av福利片在线观看| 色综合亚洲欧美另类图片| 一区二区三区四区激情视频 | 2021天堂中文幕一二区在线观| 亚洲狠狠婷婷综合久久图片| 国产成人aa在线观看| a级一级毛片免费在线观看| 动漫黄色视频在线观看| 他把我摸到了高潮在线观看| 久久天躁狠狠躁夜夜2o2o| 久久久久精品国产欧美久久久| 波野结衣二区三区在线| 日本 av在线| 久久6这里有精品| 亚洲av不卡在线观看| 99久久九九国产精品国产免费| 美女免费视频网站| 欧美zozozo另类| 直男gayav资源| 国产精品1区2区在线观看.| 香蕉av资源在线| 国产亚洲精品久久久com| 欧美成人一区二区免费高清观看| 欧美高清性xxxxhd video| 国产白丝娇喘喷水9色精品| 他把我摸到了高潮在线观看| 国产成人影院久久av| 久久精品国产鲁丝片午夜精品 | 狠狠狠狠99中文字幕| 亚洲精品乱码久久久v下载方式| 国产黄片美女视频| 99久久精品热视频| 在线观看午夜福利视频| 日本一本二区三区精品| 日本黄大片高清| 九色成人免费人妻av| 国产美女午夜福利| 免费观看精品视频网站| 成人国产一区最新在线观看| 日韩强制内射视频| 夜夜夜夜夜久久久久| 天天一区二区日本电影三级| 男女做爰动态图高潮gif福利片| videossex国产| 老师上课跳d突然被开到最大视频| 在线观看免费视频日本深夜| 美女黄网站色视频| 好男人在线观看高清免费视频| 日韩一区二区视频免费看| 国产男靠女视频免费网站| 嫩草影院新地址| 日本三级黄在线观看| 在线a可以看的网站| 亚洲18禁久久av| 简卡轻食公司| 尾随美女入室| 久久国产乱子免费精品| 嫩草影院入口| 看免费成人av毛片| 成人高潮视频无遮挡免费网站| 深夜a级毛片| 最近最新免费中文字幕在线| 国产高清三级在线| 国产一级毛片七仙女欲春2| 欧美性猛交╳xxx乱大交人| 成人鲁丝片一二三区免费| 国产毛片a区久久久久| 88av欧美| 亚洲人与动物交配视频| 日韩,欧美,国产一区二区三区 | 在线免费观看的www视频| 亚洲欧美日韩高清专用| 中国美白少妇内射xxxbb| 国产精品av视频在线免费观看| 亚洲欧美清纯卡通| 亚洲av中文字字幕乱码综合| av在线观看视频网站免费| 丝袜美腿在线中文| 男女那种视频在线观看| 日韩精品青青久久久久久| 国产一区二区三区在线臀色熟女| 女人被狂操c到高潮| 国产乱人伦免费视频| 久久精品国产清高在天天线| 好男人在线观看高清免费视频| 欧美日本视频| 男女下面进入的视频免费午夜| 禁无遮挡网站| 色综合婷婷激情| 午夜福利欧美成人| 成年人黄色毛片网站| 成人欧美大片| 亚洲精品成人久久久久久| 最近最新免费中文字幕在线| 成人一区二区视频在线观看| 国产真实伦视频高清在线观看 | 又黄又爽又免费观看的视频| 我的女老师完整版在线观看| 麻豆国产av国片精品| 国产精品人妻久久久久久| 男女下面进入的视频免费午夜| 日韩亚洲欧美综合| 在线播放无遮挡| 久久精品夜夜夜夜夜久久蜜豆| 99在线视频只有这里精品首页| 国产精品国产高清国产av| 久久人人爽人人爽人人片va| 久久久国产成人精品二区| 波野结衣二区三区在线| 变态另类丝袜制服| 99久国产av精品| 日本黄大片高清| 美女xxoo啪啪120秒动态图| 欧洲精品卡2卡3卡4卡5卡区| 老熟妇乱子伦视频在线观看| 国产精品人妻久久久影院| 欧美中文日本在线观看视频| h日本视频在线播放| 女生性感内裤真人,穿戴方法视频| 成人二区视频| 免费观看在线日韩| 免费观看人在逋| 真实男女啪啪啪动态图| 国产精品无大码| 美女cb高潮喷水在线观看| 尾随美女入室| 国产伦精品一区二区三区视频9| 日韩亚洲欧美综合| 欧美精品啪啪一区二区三区| 精品午夜福利在线看| 久久久国产成人精品二区| 永久网站在线| 久久久精品欧美日韩精品| 老司机福利观看| av.在线天堂| 日韩强制内射视频| 一本精品99久久精品77| 日韩欧美国产一区二区入口| 欧美激情在线99| 日韩一区二区视频免费看| 午夜a级毛片| 久久精品国产亚洲网站| 国产伦在线观看视频一区| 国产蜜桃级精品一区二区三区| 精品欧美国产一区二区三| 欧美绝顶高潮抽搐喷水| 哪里可以看免费的av片| 国产v大片淫在线免费观看| av在线蜜桃| 69av精品久久久久久| 小蜜桃在线观看免费完整版高清| 日韩中字成人| 国产成人aa在线观看| 在线a可以看的网站| 我要搜黄色片| 男女那种视频在线观看| 亚洲精品一卡2卡三卡4卡5卡| 国国产精品蜜臀av免费| 亚洲不卡免费看| 嫩草影视91久久| 亚洲精品久久国产高清桃花| 日韩欧美在线乱码| 国产精品,欧美在线| 干丝袜人妻中文字幕| 我的老师免费观看完整版| 丝袜美腿在线中文| 亚洲av中文字字幕乱码综合| 亚洲av.av天堂| 极品教师在线视频| ponron亚洲| 国产精品久久视频播放| 精品人妻一区二区三区麻豆 | 国产精品电影一区二区三区| 国产主播在线观看一区二区| 免费在线观看成人毛片| 国产视频一区二区在线看| а√天堂www在线а√下载| 免费观看精品视频网站| 精品一区二区免费观看| 国产精品伦人一区二区| 老熟妇仑乱视频hdxx| 两个人视频免费观看高清| 免费一级毛片在线播放高清视频| 色综合亚洲欧美另类图片| 性色avwww在线观看| 麻豆成人午夜福利视频| 日韩精品中文字幕看吧| 淫秽高清视频在线观看| 大又大粗又爽又黄少妇毛片口| 国产伦精品一区二区三区视频9| 久久99热6这里只有精品| 两个人的视频大全免费| 亚洲aⅴ乱码一区二区在线播放| 波多野结衣巨乳人妻| 美女高潮喷水抽搐中文字幕| 一卡2卡三卡四卡精品乱码亚洲| 男人舔奶头视频| 性色avwww在线观看| 欧美人与善性xxx| 天堂动漫精品| 在线播放国产精品三级| 久久国产乱子免费精品| 成人性生交大片免费视频hd| 麻豆成人av在线观看| 最近中文字幕高清免费大全6 | 亚洲专区中文字幕在线| 高清毛片免费观看视频网站| 内射极品少妇av片p| 亚洲av日韩精品久久久久久密| 亚洲第一电影网av| 精华霜和精华液先用哪个| 美女 人体艺术 gogo| 亚洲第一区二区三区不卡| 日日干狠狠操夜夜爽| 中国美白少妇内射xxxbb| 亚洲熟妇熟女久久| 久久久久精品国产欧美久久久| 国产 一区精品| 日本撒尿小便嘘嘘汇集6| av在线蜜桃| 啪啪无遮挡十八禁网站| 91在线观看av| 国产一区二区亚洲精品在线观看| 国产在线精品亚洲第一网站| 精品一区二区免费观看| 色av中文字幕| 最好的美女福利视频网| 午夜爱爱视频在线播放| 一区二区三区激情视频| 欧美成人a在线观看| 日韩,欧美,国产一区二区三区 | 午夜免费激情av| 两性午夜刺激爽爽歪歪视频在线观看| 国产一区二区亚洲精品在线观看| 久久久久久久亚洲中文字幕| 丰满的人妻完整版| 日韩欧美在线二视频| 日韩中字成人| 51国产日韩欧美| 国内少妇人妻偷人精品xxx网站| 日本一二三区视频观看| 国语自产精品视频在线第100页| 午夜亚洲福利在线播放| 色5月婷婷丁香| 国产一区二区三区在线臀色熟女| 人妻丰满熟妇av一区二区三区| 久久久久久久久大av| 别揉我奶头 嗯啊视频| 久久精品国产自在天天线| 成年女人毛片免费观看观看9| av在线观看视频网站免费| 午夜精品在线福利| av天堂中文字幕网| 毛片一级片免费看久久久久 | 亚洲精品一区av在线观看| 午夜激情福利司机影院| 大型黄色视频在线免费观看| 精品欧美国产一区二区三| 精品一区二区三区视频在线观看免费| 性色avwww在线观看| 欧美不卡视频在线免费观看| 欧美zozozo另类| 久久久久国产精品人妻aⅴ院| 校园人妻丝袜中文字幕| 久久久久久久久久成人| av在线天堂中文字幕| 看黄色毛片网站| 久久久久久大精品| 波多野结衣巨乳人妻| 九九热线精品视视频播放| 中文字幕av在线有码专区| 久久久久久国产a免费观看| 午夜视频国产福利| 精品久久久久久久久av| 国产激情偷乱视频一区二区| 美女高潮喷水抽搐中文字幕| 淫妇啪啪啪对白视频| 欧美精品国产亚洲| 啦啦啦观看免费观看视频高清| 亚州av有码| 九色国产91popny在线| 日韩欧美三级三区| 真实男女啪啪啪动态图| 国产成年人精品一区二区| 国产 一区 欧美 日韩| 亚洲精品色激情综合| 久久国产精品人妻蜜桃| 久久久久久久久大av| 别揉我奶头 嗯啊视频| 久久99热这里只有精品18| 99久久成人亚洲精品观看| 国产精品国产三级国产av玫瑰| 在线播放国产精品三级| 91麻豆精品激情在线观看国产| 国产精品永久免费网站| 久久久久久久久大av| 乱码一卡2卡4卡精品| 国产精品,欧美在线| 亚洲欧美精品综合久久99| www日本黄色视频网| 蜜桃亚洲精品一区二区三区| 国产 一区精品| 亚洲久久久久久中文字幕| 免费电影在线观看免费观看| 日本在线视频免费播放| 91麻豆精品激情在线观看国产| 观看免费一级毛片| 日韩欧美三级三区| 日韩大尺度精品在线看网址| 国产成人福利小说| 免费一级毛片在线播放高清视频| or卡值多少钱| 国产精品99久久久久久久久| 日日干狠狠操夜夜爽| 亚洲av二区三区四区| 午夜免费激情av| 我的老师免费观看完整版| 热99re8久久精品国产| 午夜福利在线在线| 久久精品国产99精品国产亚洲性色| 在线观看一区二区三区| 99久久成人亚洲精品观看| 国产精品一及| 日韩大尺度精品在线看网址| 精品久久国产蜜桃| 国产精品不卡视频一区二区| а√天堂www在线а√下载| 黄色配什么色好看| 久久午夜福利片| 国产一级毛片七仙女欲春2| 国产美女午夜福利| 亚洲美女黄片视频| 亚洲av电影不卡..在线观看| 国产蜜桃级精品一区二区三区| 午夜久久久久精精品| 久久人人爽人人爽人人片va| 中文字幕av在线有码专区| 校园人妻丝袜中文字幕| 伦理电影大哥的女人| 欧美xxxx黑人xx丫x性爽| 国产一区二区三区视频了| 午夜精品久久久久久毛片777| 日日啪夜夜撸| 欧美成人免费av一区二区三区| 国产精品av视频在线免费观看| 亚洲人成网站高清观看|