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

    The skeleton of 5,7-fused bicyclic imidazole-diazepine for heatresistant energetic materials

    2023-10-09 04:29:50XiaoxiaoZhengYubingXueChanghaoDaiHongweiYangGuangbinCheng
    Defence Technology 2023年9期

    Xiaoxiao Zheng,Yubing Xue,Changhao Dai,Hongwei Yang ,Guangbin Cheng

    School of Chemistry and Chemical Engineering,Nanjing University of Science and Technology,Nanjing,210094,China

    Keywords:Imidazole-diazepine Heat-resistant material 5,7-Fused skeleton Energetic materials

    ABSTRACT In light of the low yields and complex reaction routes of some well-known 5,5-fused and 5,6-fused bicyclic compounds,a series of 5,7-fused bicyclic imidazole-diazepine compounds were developed with high yields by only two efficient steps.Significantly,the seven-membered heterocyclic ring has a stable energetic skeleton with multiple modifiable sites.However,the 5,7-fused bicyclic energetic compounds were rarely reported in the area of energetic materials.Three neutral compounds 1,2 and 4 were synthesized in this work.To improve the detonation performances of the 5,7-fused neutral compounds,corresponding perchlorate 1a and 2a were further developed.The physicochemical and energetic performances of all newly developed compounds were experimentally determined.All newly prepared energetic compounds exhibit high decomposition temperatures (Td: 243.8-336 ℃) and low mechanical sensitivities (IS: >15 J, FS: >280 N).Among them,the velocities performances of 1a(Dv=7651 m/s)and 4(Dv=7600 m/s)are comparable to that of typical heat-resistant energetic material HNS(Dv=7612 m/s).Meanwhile,the high decomposition temperature and low mechanical sensitivities(Td=336 ℃; IS=32 J; FS >353 N) of 4 are superior to that of HNS (Td=318 ℃; IS=5 J; FS=250 N).Hence,the 5,7-fused bicyclic compounds with high thermostability,low sensitivities and adjustable detonation performance have a clear tendency to open up a new space for the development of heatresistant energetic materials.

    1.Introduction

    Energetic materials (EMs) are a class of compounds or formulations with oxidants and combustible components,which can release tremendous energy independently during the chemical reactions [1,2].Generally,to meet the demand for rapid and massive energy release,energetic groups such as nitro (-NO2),nitramine (-NHNO2) and azide (-N3) are installed to decorate the skeleton of Ems[3].However,the preoccupation with energy might have undesirable consequences,giving the material unsatisfactory stabilities.In other words,the two objectives of high detonation performance and excellent stabilities for EMs are mutually incompatible.For a long time,researchers pay more attention to improving detonation performance at the expense of stability.As a result,the heat-resistant materials (Td>300 ℃) account for less than 10% of the reported energetic materials (the samples of EMs>600) [4].In the context of the increasing demand for the exploitation of scarce resources,heat-resistant materials have been concerned gradually in recent years.However,the poor thermostability of EMs can restrict their application in several task-specific such as hypersonic weapons,deep mining and aerospace engineering [5].In this way,the innovation of the design and development for novel heat-resistant energetic materials is of great significance.

    Fused-ring energetic skeletons are generally considered to be the polycyclic system composed of two or more heterocycles sharing at least one ring edge [6].In the past few decades,the fused-ring derivatives had been chiefly applied in pharmaceutical chemistry,bioluminescence and photoelectric [7-9].In recent years,fused-ring energetic compounds with high thermostability have been designed gradually [10].The fused-ring energetic compounds could lead to a high decomposition temperature due to their large π-conjugated system,which has a tendency of application in the area of heat-resistant energetic materials [11-13].Additionally,the amino group as an insensitive substituent group plays a vital role in maintaining molecular stability,which can provide hydrogen donors to form intra-and intermolecular hydrogen bond interactions [14-16].For example,DNPP(Td=336 ℃),PTX (Td=246 ℃) and DNPPDA (Td=325 ℃) all demonstrate that the fused energetic skeleton with amino group is an outstanding frame structure to beget the thermostable energetic materials (Fig.1) [17-19].Regrettably,the development of fusedring energetic materials is not plain sailing due to the obstacles of practical factors.The prohibitively high price of the reaction substrates,low yields,time-consuming synthesis steps or complicated reaction routes are some principal causes that hinder the practical application of fused-ring energetic compounds [20,21].

    Fig.1.Comparison of previous work (5,5-fused and 5,6-fused bicyclic energetic compounds) and current work (5,7-fused bicyclic energetic compound 4).

    At present,the bicyclic fused-ring framework used in the field of EMs is dominated by 5,5-fused and 5,6-fused bicyclic skeletons[22].The 5,7-fused bicyclic energetic skeleton was rarely reported.Considering that the seven-membered heterocyclic ring is a stable energetic skeleton with multiple modifiable sites[23],we focus our attention on the design and development of 5,7-fused bicyclic compounds.In this study,a facile synthetic method of 5,7-fused bicyclic compounds in high yields via only a two-step reaction was reported.Furthermore,in view of the costly price of reaction substrates for a certain number of reported 5,5-fused and 5,6-fused compounds,we selected the commercially available and economical starting materials for the design of 5,7-fused compounds.Here,compounds1and2were developed based on the skeleton of 5,7-fused bicyclic imidazole-diazepine.To obtain the compound with high energy properties,the 5,7-fused compound4with a nitro group was further synthesized [24].Significantly,the unusual skeleton of 5,7-fused bicyclic imidazole-diazepine has unique advantages in heat-resistant energetic materials due to its powerful π-delocalization and symmetric modifiable sites.All the compounds in this work gather the performances of high thermostability and low mechanical sensitivity.Among them,the decomposition temperatures of 2 and 4 are higher than 300 ℃.Overall,the high decomposition temperatures and facile synthesis of 5.7-fused bicyclic compounds are competitive and attractive in the field of heat-resistant materials.

    2.Experimental part

    2.1.Safety precaution

    In the work,all newly prepared energetic compounds are a class of explosive materials with high energy.So,necessary safety precautions should be noted.The experiment ought to be handled in the fume hood during the whole experiments.Additionally,a protective helmet,gloves,and explosion-proof glass ought to be equipped all the time.The small scale trials are necessary.Last,the mechanical force such as impact,friction and spark should be avoided.

    2.2.Synthesis

    2.2.1.Synthesis of imidazo[4,5-e][1,3]diazepine-4,6,8-triamine (1)

    1was synthesized according to the previous Ref.[25].

    2.2.2.Synthesis of 4,6,8-triaminoimidazo[4,5-e][1,3]diazepin-1-ium(1a)

    The mixture of1(0.354 g,2.0 mmol) in 10 mL distilled water was heated in an oil bath at 60 ℃.The strong acid HClO4was added to the mixture dropwise under constant stirring until the emergence of a clear solution,and then cooled down the reaction.The white crystal was precipitated after standing for half an hour.The filtration was performed to obtain1a(0.487 g,90%).1H NMR chemical shifts (δ,ppm;500 MHz,DMSO-d6):8.50 (s,2H),8.32(s,2H),8.21 (s,2H),7.89 (s,2H);13C NMR chemical shifts (δ,ppm;125 MHz,DMSO-d6):163.80,155.90,140.50,131.20.IR(KBr):3240,2900,1896,1769,1552,1494,1448,1323,1272,1209,1006,962,929,830,775,696,583,551 cm-1.Element analysis (%): C6H9N7O4Cl(MW: 278.63),calculated result: C 25.86,H 3.26,N 35.19;found: C 25.91,H 3.31,N 35.04.

    2.2.3.Synthesis of 5-(N′-nitrocarbamimidoyl)-1H-imidazole-4-carboxamide (1b)

    5 mL 100% HNO3was added into the 25 mL flask under the lowtemperature bath of -10 ℃.Compound1(0.5 g,2.5 mmol) was patted into the cold flask in small amounts under constant stirring.After completing the feeding,the mixture was stirred for another 3 h at 0 ℃.Finally,the reaction system was quenched by 8 g crushed ice to precipitate a white solid.The product was collected after filtration,washed with water (3 × 3 mL),and dried in air to afford powder (0.25 g,44.7%).1H NMR chemical shifts (δ,ppm;500 MHz,DMSO-d6):8.29(s,2H),8.05(s,H),6.91(s,2H);13C NMR chemical shifts (δ,ppm;125 MHz,DMSO-d6): 163.38,162.73,157.42,144.21,133.03.IR (KBr): 3437,3307,1671,1499,1351,1160,848,657,483 cm-1.Element analysis (%): C5H6N6O3(MW: 198.14),calculated result:C 30.31,H 3.05,N 42.42;found:C 30.20,H 3.31,N

    42.44.

    2.2.4.Synthesis of imidazo[4,5-e][1,3]diazepine-2,4,6,8-tetraamine(2)

    Guanidine:The solid of sodium methoxide(1.080 g,20.0 mmol)was added to 25 mL MeOH and stirred until a clear solution appeared,then guanidine hydrochloride (1.910 g,20.0 mmol) was mixed into the solution.The precipitates emerged and were filtered off after the neutralization reaction for 1 h.The filtrate was collected for the next step.

    2-amino-1H-imidazole-4,5-dicarbonitrile (1.330 g,10.0 mmol)was suspended in the above fresh filtrate,then the reaction mixture was warmed slowly in an oil bath until reflux for 20 h.A plentiful of white deposits appeared at the termination of the reaction.The solids were filtrated,washed by MeOH(4×5 mL),and dried in the air to yield a pure yellow product2(1.576 g,82%).1H NMR chemical shifts(δ,ppm;500 MHz,DMSO-d6):7.31(s,2H),7.01(s,2H),6.55(s,2H),5.06 (s,2H);13C NMR chemical shifts (δ,ppm;125 MHz,DMSO-d6):164.21,159.11,156.95,135.96.IR(KBr):3465,3019,1609,1489,1386,1076,730,571,482 cm-1.Element analysis(%):C6H8N8(MW: 192.19),calculated result: C 37.50,H 4.20,N 58.31;found: C 37.31,H 4.31,N 58.45.

    2.2.5.Synthesis of 2,6,8-triamino-4-oxo-4,5-dihydroimidazo[4,5-e][1,3]diazepine-1,7-diium (2a)

    The related synthetic method of2ais similar to that of compound1a.The white filter cake was confirmed to be2a(0.600 g,76%).1H NMR chemical shifts(δ,ppm;500 MHz,DMSO-d6):9.20(s,2H),8.80(s,2H),8.20(s,2H),7.06(s,2H);13C NMR chemical shifts(δ,ppm;125 MHz,DMSO-d6):163.38,162.73,157.42,144.21,133.03.IR(KBr):3323,3148,2803,1598,1328,738,581,464 cm-1.Element analysis (%):C6H9N7O9Cl2(MW:394.08),calculated result:C 18.29,H 2.03,N 24.88;found: C 18.03,H 2.10,N 25.04.

    2.2.6.Synthesis of 2-nitro-1H-imidazole-4,5-dicarbonitrile (3)

    3was synthesized according to the previous Ref.[26].

    2.2.7.Synthesis of 2-nitroimidazo[4,5-e][1,3]diazepine-4,6,8-triamine (4)

    2-nitro-1H-imidazole-4,5-dicarbonitrile (1.631 g,10.0 mmol)was suspended in a methanol solution of guanidine.The reaction was kept at 70 ℃ in an oil bath for 20 h.A great many yellow precipitates were formed and filtered out in the end.The filter cakes were collected,washed by MeOH(4×5 mL),and air dried to obtain a yellow product4(1.667 g,80%).1H NMR chemical shifts(δ,ppm;500 MHz,DMSO-d6): 8.04 (s,4H),7.54 (s,2H);13C NMR chemical shifts (δ,ppm;125 MHz,DMSO-d6): 164.04,158.74,156.96,135.80.IR(KBr):3438,2617,1713,1584,857,802,537 cm-1.Element analysis (%): C6H6N8O2(MW: 222.17),calculated result: C 32.44,H 2.72,N 50.44;found: C 32.47,H 2.80,N 50.60.

    3.Results and discussion

    3.1.Synthesis

    In this study,the synthetic routes of all prepared compounds were presented in Fig.2.Compound1was prepared according to the reported literature.Then,various nitration conditions such as fuming HNO3,100% HNO3,H2SO4/98% HNO3and KNO3/H2SO4were tried for the nitration of1,but none of the nitramino groups were introduced successfully.Unexpectedly,a ring-opening product1bwas obtained accidently by treating1with 100% HNO3.To improve the density and detonation performances of1,perchlorate1awas synthesized by using cheap and oxygen-rich perchloric acid [27].Since high nitrogen content can generally lead to an increase in enthalpy of formation [28],the introduction of an amino group at the unmodified position ofp1was tried to gain an optimized substratep2.Hence,compound2was synthesized by the cyclization of two cyano-groups attached top2in the presence of guanidine and methanol.Considering that the nitro group is a favorable energetic substituent in realizing high energy density materials[29],we tried multiple oxidation systems (H2SO4/NaNO2,Oxone and H2O2/H2SO4) to oxidize the amino groups of2,but the long trials ended in failure.Therefore,based on the successful case of1a,perchlorate2awas obtained by using perchloric acid as the neutralization reagent to improve its detonation performance.Interestingly,one of the amino groups in the seven-membered ring of2was converted to a carbonyl group,which might be due to the strong oxidability of perchloric acid.Significantly,another strategy was proposed to introduce the nitro group in the imidazole ring based on optimization of the reaction substrate[30].First,2-nitro-1H-imidazole-4,5-dicarbonitrile (3) was synthesized by oxidating the amino group in substratep2,and then the optimized substrate with a nitro group was participated in the heterocyclization reaction to obtain compound4.All newly prepared compounds were characterized by NMR spectra (Nuclear Magnetic Resonance:1H and13C),EA (Element Analysis) and IR (Infrared Spectroscopy).

    Fig.2.Synthesis of 5,7-fused energetic compounds 1-4 and 1a-2a.

    3.2.Single crystal X-ray diffraction

    Single crystals of compounds1a,2aand4were cultivated and characterized by single crystal X-ray diffraction to confirm their special structures and further explore their structure-property relationships(see Fig.3).Crystal1a?H2O and2a?2H2O were afforded by slow cooling of the reaction system,while crystal4?DMSO grew from its saturated solution after slow evaporation at atmospheric temperature.

    Fig.3.(a)-(c) The X-ray crystal structures for 1a?H2O,2a?2H2O,4?DMSO,respectively;(d)-(f) The packing diagrams and intermolecular hydrogen bonding (H-bonding,green dotted line) of 1a?H2O,2a?2H2O,4?DMSO,respectively (the solvents were removed for clarity);(g)-(h) The enlarged view form packing diagrams of 1a?H2O,2a?2H2O,4?DMSO,respectively.

    Crystal1a?H2O belongs to the triclinic crystal system,the P-1 space group,and its crystal density was calculated to be 1.754 g/cm3at 170 K.In Fig.3(d),the packing diagram of the cationic components in1acan be described as face-to-face stacking,and the attached anions are embedded in those layers at regular intervals.To obtain more details about crystal packing arrangements,the enlarged view of packing diagrams is displayed in Fig.3(g).The two adjacent layers are arranged in parallel with a tight distance of 3.14 ? (d1,2).Actually,each horizontal layer consists of a discontinuous tier with a tiny distance of 0.19 ?(d1,3).All atoms of cationic components in crystal1a?H2O lie in the same plane,which can be confirmed reasonably by the torsion angles of C1-N1-C3-C4 179.1(2)o,C1-N2-C2-C3 0.1(3)o,C6-N4-C5-C3 0.9(5)o(ESI).The crystal system of2a?2H2O is belong to the monoclinic crystal system,P21/n space group with a calculated density of 1.873 g/cm3at 170 K.The seven-membered ring is slightly distorted due to the difference between the carbonyl and amino groups in polarity(Fig.3(b)).The packing diagram of the cationic components in2acan be regarded as wave-like stacking.Among them,the neighboring fragments are linked head to tail by related anions(Fig.4(e)).In the zoomed-in view of the packing diagram of crystal2a?2H2O(Fig.3(h)),the distance between the adjacent parallel plane is measured to be 3.23 ?(d1,2=d3,4)and the dihedral angle between two intersecting planes is measured to be 60°.The crystal system of4?DMSO is the same as that of1a?H2O,which belongs to the triclinic crystal system,P-1 space group,and its crystal density was calculated to be 1.646 g/cm3at 170 K (Fig.3(c)).The low density may be caused by the presence of DMSO.As shown in the packing diagram(Fig.3(i)),crystal4can be viewed as flat“Z”stacking.The interlayer distance between two horizontally staggered planes is measured to be 3.29 ?(d1,2).Besides,there is a slight slide between two contiguous layers connected by DMSO with a distance of 1.68 ?(d2,3).

    Fig.4.DSC plots of newly developed compounds.

    To further explore the structure-property relationship of these compounds,the intermolecular interaction and stacking distances of crystals1a,2aand4were analyzed to elucidate their effects on the energetic properties [31].The hydrogen-bonding interaction(green dotted line) between N and O atoms in molecules can be observed in Fig.3(d)-(f).Among them,2ahas the largest number of hydrogen bonds,resulting in its highest density.Furthermore,the interlayer distances of these packing molecules are all within the scope of the interaction distance of π-π stacking(Dd<4 ?)[32].The strong hydrogen bond interaction and π-π stacking interaction in molecules contribute to the thermostabilities and insensitivities of compounds,thus the desirable thermostability and the low sensitivities of compounds 1a,2a and 4 can be explained.

    3.3.Detonation property and mechanical sensitivity

    A series of physicochemical properties of the developed compounds were displayed in Table 1.All synthetic products were fully dried at the standard condition to remove solvents before performance tests.

    Table 1 The energetic and sensitive properties of 1-4 and 1a-2a.

    Thermal stability is one of the important parameters that should be concerned in practical applications [33].Generally,the decomposition temperature of contemporary EMs is required more than 180 ℃ [34-36],because the low thermostability may limit their practical application in transportation,storage and use safety[37].In this work,decomposition temperatures of all developed compounds were measured depending on DSC (Differential Scanning Calorimetry) device in a nitrogen atmosphere at a heating rate of 5oC/min.All compounds showed good thermostability and their onset decomposition temperatures were in the range of 243.8-336 ℃.Among them,neutral compound4displays the highest onset decomposition temperature of 336 ℃,which is superior to that of the typical heat-resistant energetic material HNS(Td=318 ℃).In term of energetic salts,the decomposition temperature of perchlorate2a(Td=293.0 ℃) shows an increase of 49.2 ℃ compared with perchlorate1a(Td=243.8 ℃)(Fig.4).In addition to thermal stability,mechanical sensitivities are another factor that needs to be considered at the security level,including impact sensitivity (IS) and friction sensitivity (FS) [38].The mechanical sensitivities of these developed compounds were put to the test through the standard BAM method.All compounds showed low mechanical sensitivities (IS: >15 J,FS: >280 N),which are all superior to that of the classical heat-resistant energetic material HNS (IS: 5 J,FS: >250 N).

    To get a more profound understanding of the high heat resistance and the low mechanical sensitivity of these compounds,the analysis of Hirshfeld surfaces and 2D fingerprint plots for compounds1a,2aand4were employed respectively (the solvent molecules of the three compounds were not selected in the formation process of Hirshfeld surfaces and fingerprint plots)[39].The weak intermolecular interactions are visualized graphically in Fig.5(a)-Fig.5(c),the prominent red spots can be found on the side edges of the plates,representing that there are hydrogen bonding interactions in the corresponding regions [40],and more detailed information about the weak intermolecular interaction is presented in Fig.5(d)-Fig.5(g).As can be seen in the 2D fingerprint plots,several needle-like shapes could be observed in the bottom lefthand corner of the images[41],which are formed by the N…H and O…H interactions between molecules.Further combined with the data statistics from Fig.5(g),the hydrogen bonding interactions of the three compounds are dominant in the weak interactions and the percentage of 1a and 2a account for well over half of the total.Overall,the extensive hydrogen bonding interactions can not only stabilize the stacking system to improve the heat resistance performance,but also short the interlayer distance of molecular packing in crystal to decrease the mechanical sensitivity [42].Therefore,the high heat resistant and the low mechanical sensitivities of compounds 1a,2a,and 4 can be reasonably explained.

    Fig.5.(a),(b)Hirshfeld surfaces for 1a,2a and 4(view from the front,vertical,and side,respectively);(c),(d)2D fingerprint plots for 1a,2a and 4,respectively;(e)The bar graph of individual atomic contacts percentage for 1a,2a and 4.

    Detonation velocity (Dv) and detonation pressure (P) are two essential indicators of EMs,because they can give us the basis for evaluating the detonation power of the developed explosives.First of all,density (ρ) and enthalpy of formation (ΔHf) should be available before calculating the detonation performance [43].In this work,the densities of all compounds ranged from 1.63 to 1.87 g/cm,which were tested by a gas pycnometer in a helium atmosphere.Furthermore,all compounds except carbonyl ionic compounds 2a present positive enthalpies of formation.After the two main variables (ρ and ΔHf) were obtained,the detonation properties of all newly synthesized energetic compounds were calculated by the EXPL05 v6.01 program.Among them,the detonation velocities of compounds 1a (Dv=7651 m/s) and 4 (Dv=7600 m/s) are comparable to that of HNS (Dv=7612 m/s).

    4.Conclusions

    Considering the shortcoming of some reported 5,5-fused and 5,6-fused bicyclic compounds with intricate synthesis and low yield,we designed and developed a class of 5,7-fused bicyclic compounds with high output and convenient synthesis (two-step synthesis).As yet,the energetic skeleton of 5,7-fused bicyclic imidazole-diazepine is rarely reported in the field of EMs.In this study,three neutral 5,7-fused imidazole-diazepine compounds1,2,and4were synthesized by the cyclization of two cyano-groups attached to related imidazole substrates,respectively.Since several nitrification and oxidation reactions had been attempted but failed,perchlorate1aand2awere further synthesized to improve their detonation performances.Furthermore,the physicochemical and energetic properties of the developed compounds were experimentally and theoretically ascertained.Among them,the detonation performances of compounds1a(Dv=7651 m/s)and4 (Dv=7600 m/s) are comparable to that of the classical heatresistant energetic material HNS (Dv=7612 m/s).Moreover,4possesses the high decomposition temperature and low mechanical sensitivities (Td=336 ℃;IS=32 J;FS>353 N),which are superior to that of HNS (Td=318 ℃;IS=5 J;FS=250 N).On the one hand,the large conjugation of the 5,7-fused bicyclic skeleton can enhance the thermal stability and facilitate the intermolecular π-π stacking to reduce their sensitivity.On the other hand,the properties of 5,7-fused bicyclic compounds can be reasonably adjusted due to the multiple modifiable sites of the diazepine ring.In brief,the application of the 5,7-fused bicyclic energetic skeleton can provide a new pioneering guide for the design of novel energetic materials and effectively alleviate the growing demand for heat-resistant energetic materials.

    Declaration of competing interest

    The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.

    Acknowledgments

    This work was completed with support from the National Natural Science Foundation of China (Grant No.22075143,21875110),as well as the Science Challenge Project(Grant No.TZ2018004).H.Yang thanks the Qing Lan Project for the grant.

    Appendix B.Supplementary data

    Supplementary data to this article can be found online at https://doi.org/10.1016/j.dt.2022.09.003.

    中文字幕最新亚洲高清| 黄片无遮挡物在线观看| 一区二区三区免费毛片| 97在线视频观看| 18+在线观看网站| 免费不卡的大黄色大毛片视频在线观看| 国产亚洲精品第一综合不卡 | 欧美日韩av久久| 看免费成人av毛片| 在线观看国产h片| 大又大粗又爽又黄少妇毛片口| 中文欧美无线码| 18禁动态无遮挡网站| 久久99精品国语久久久| 丰满少妇做爰视频| 国产一区亚洲一区在线观看| 两个人的视频大全免费| 欧美激情 高清一区二区三区| 18+在线观看网站| 少妇被粗大的猛进出69影院 | 国产精品麻豆人妻色哟哟久久| 亚洲第一av免费看| 热99久久久久精品小说推荐| 伦精品一区二区三区| av在线app专区| 91成人精品电影| 精品久久久噜噜| 涩涩av久久男人的天堂| 午夜免费男女啪啪视频观看| 国产精品国产av在线观看| 亚洲av综合色区一区| 午夜精品国产一区二区电影| 在线观看免费视频网站a站| 欧美日韩视频高清一区二区三区二| 在线观看国产h片| 婷婷成人精品国产| 99九九在线精品视频| av不卡在线播放| 大香蕉97超碰在线| 一级片'在线观看视频| 亚洲综合色网址| 久久精品国产亚洲av涩爱| 亚洲,欧美,日韩| 999精品在线视频| 又大又黄又爽视频免费| 十分钟在线观看高清视频www| 十分钟在线观看高清视频www| 免费人妻精品一区二区三区视频| 少妇被粗大猛烈的视频| 欧美3d第一页| 国产亚洲欧美精品永久| 精品少妇久久久久久888优播| 各种免费的搞黄视频| 七月丁香在线播放| 国产精品三级大全| 日日爽夜夜爽网站| 哪个播放器可以免费观看大片| 成人黄色视频免费在线看| 丝袜喷水一区| 中国美白少妇内射xxxbb| 美女国产视频在线观看| 色5月婷婷丁香| 国产熟女午夜一区二区三区 | 能在线免费看毛片的网站| 在线播放无遮挡| 黄色欧美视频在线观看| xxxhd国产人妻xxx| 亚洲综合色惰| 精品国产露脸久久av麻豆| 精品少妇内射三级| 久热这里只有精品99| a级毛片在线看网站| 日韩 亚洲 欧美在线| 午夜久久久在线观看| 国产免费一区二区三区四区乱码| 少妇被粗大的猛进出69影院 | freevideosex欧美| 久久精品夜色国产| 国语对白做爰xxxⅹ性视频网站| 日本wwww免费看| 一区二区av电影网| 九草在线视频观看| 久久精品人人爽人人爽视色| 夜夜看夜夜爽夜夜摸| 夫妻午夜视频| 男女边吃奶边做爰视频| 久久婷婷青草| 久久精品国产鲁丝片午夜精品| 亚洲欧美一区二区三区国产| 夫妻午夜视频| 最黄视频免费看| 国产一区二区在线观看日韩| av福利片在线| 99热全是精品| 国产一区二区三区av在线| 亚洲美女视频黄频| av又黄又爽大尺度在线免费看| 国产成人91sexporn| 亚洲精品美女久久av网站| 久热这里只有精品99| 成年女人在线观看亚洲视频| 久久精品熟女亚洲av麻豆精品| 亚洲第一av免费看| 毛片一级片免费看久久久久| 在线观看国产h片| 国产免费现黄频在线看| 赤兔流量卡办理| 欧美精品高潮呻吟av久久| 免费人成在线观看视频色| 亚洲av在线观看美女高潮| 亚洲综合色网址| 男女无遮挡免费网站观看| 高清不卡的av网站| 成年人免费黄色播放视频| 亚洲怡红院男人天堂| 欧美日韩视频精品一区| 天堂中文最新版在线下载| 日韩av在线免费看完整版不卡| 午夜视频国产福利| 国产欧美亚洲国产| 亚洲国产精品成人久久小说| 国产欧美日韩一区二区三区在线 | 狠狠婷婷综合久久久久久88av| 在线精品无人区一区二区三| 国产精品免费大片| 亚洲经典国产精华液单| 最新中文字幕久久久久| 国产黄频视频在线观看| 永久网站在线| 18在线观看网站| 最近最新中文字幕免费大全7| 亚洲精品久久久久久婷婷小说| 国产女主播在线喷水免费视频网站| 亚洲内射少妇av| 午夜免费男女啪啪视频观看| 大陆偷拍与自拍| 十分钟在线观看高清视频www| 精品人妻熟女毛片av久久网站| 免费大片黄手机在线观看| 国产成人精品久久久久久| 国产亚洲av片在线观看秒播厂| 精品人妻一区二区三区麻豆| 边亲边吃奶的免费视频| 精品久久久久久久久av| 永久免费av网站大全| 亚洲精品美女久久av网站| 在线 av 中文字幕| 亚洲国产色片| 中国国产av一级| 成年美女黄网站色视频大全免费 | 亚洲精品aⅴ在线观看| 国产精品久久久久久久电影| 制服丝袜香蕉在线| 日韩中字成人| 欧美日韩视频高清一区二区三区二| 男女国产视频网站| 中文天堂在线官网| 精品少妇黑人巨大在线播放| 午夜激情福利司机影院| 另类亚洲欧美激情| 一级爰片在线观看| 一级爰片在线观看| 国产免费现黄频在线看| 熟女电影av网| 插逼视频在线观看| 国产精品偷伦视频观看了| 成人国产麻豆网| 中文精品一卡2卡3卡4更新| 久久人妻熟女aⅴ| 日韩一区二区视频免费看| 99国产精品免费福利视频| 久久韩国三级中文字幕| 亚洲第一av免费看| 亚洲av在线观看美女高潮| 国产片内射在线| 一本色道久久久久久精品综合| 亚洲国产成人一精品久久久| 日韩成人av中文字幕在线观看| 日产精品乱码卡一卡2卡三| 日韩成人伦理影院| a级片在线免费高清观看视频| 欧美日韩国产mv在线观看视频| 国产一区二区在线观看日韩| 亚洲欧洲国产日韩| 纯流量卡能插随身wifi吗| 女人精品久久久久毛片| 国产有黄有色有爽视频| 久久青草综合色| 亚洲情色 制服丝袜| 一级毛片aaaaaa免费看小| 在线观看国产h片| 少妇丰满av| 久久国产精品大桥未久av| videossex国产| 亚洲精品久久久久久婷婷小说| 午夜福利,免费看| 18禁在线播放成人免费| av在线播放精品| 久久精品久久久久久噜噜老黄| 日韩在线高清观看一区二区三区| 欧美人与性动交α欧美精品济南到 | 亚洲精品一二三| 午夜福利,免费看| 大码成人一级视频| 成人国产麻豆网| .国产精品久久| 欧美xxⅹ黑人| 免费日韩欧美在线观看| 国产日韩欧美视频二区| 免费高清在线观看视频在线观看| 国产日韩欧美在线精品| 丝袜在线中文字幕| 男女免费视频国产| av在线app专区| 亚洲精品456在线播放app| kizo精华| 人人妻人人澡人人爽人人夜夜| 黄色欧美视频在线观看| 国产精品久久久久久精品电影小说| av网站免费在线观看视频| 久久人人爽av亚洲精品天堂| 亚洲色图 男人天堂 中文字幕 | 久久久久久久精品精品| 亚洲精品日本国产第一区| 丁香六月天网| 丝袜在线中文字幕| 国产亚洲午夜精品一区二区久久| 国产一级毛片在线| 一本久久精品| 亚洲精品久久久久久婷婷小说| 亚洲欧美成人精品一区二区| 少妇高潮的动态图| 黄片无遮挡物在线观看| 日韩精品有码人妻一区| 久久久精品94久久精品| 热re99久久精品国产66热6| 欧美变态另类bdsm刘玥| 婷婷色麻豆天堂久久| 日本黄大片高清| 看非洲黑人一级黄片| 国产成人av激情在线播放 | 日韩av不卡免费在线播放| 欧美精品国产亚洲| 成年av动漫网址| 久久99一区二区三区| 国产男人的电影天堂91| 国产亚洲精品久久久com| 涩涩av久久男人的天堂| 在线天堂最新版资源| 久久人人爽人人爽人人片va| 99精国产麻豆久久婷婷| 狂野欧美激情性xxxx在线观看| xxx大片免费视频| 尾随美女入室| 91成人精品电影| 国产片内射在线| 亚洲怡红院男人天堂| 一级毛片电影观看| 亚洲婷婷狠狠爱综合网| 国产精品国产三级国产av玫瑰| 国产精品偷伦视频观看了| 久久久亚洲精品成人影院| 亚洲国产色片| 51国产日韩欧美| 丰满迷人的少妇在线观看| 熟女电影av网| 亚洲精华国产精华液的使用体验| 精品视频人人做人人爽| 十八禁网站网址无遮挡| 岛国毛片在线播放| 国产精品久久久久久精品电影小说| 3wmmmm亚洲av在线观看| 看十八女毛片水多多多| 久久免费观看电影| 街头女战士在线观看网站| 少妇丰满av| 一区二区三区精品91| 日本黄大片高清| 亚洲av日韩在线播放| 国产精品 国内视频| 国产成人免费观看mmmm| 亚洲综合色惰| av又黄又爽大尺度在线免费看| 国产精品欧美亚洲77777| 一级a做视频免费观看| 搡老乐熟女国产| 久久午夜综合久久蜜桃| 丰满少妇做爰视频| 在线免费观看不下载黄p国产| 三上悠亚av全集在线观看| tube8黄色片| 久久精品人人爽人人爽视色| 亚洲av国产av综合av卡| 久久久久久久久久久久大奶| 国产精品无大码| 久久久久久伊人网av| 精品久久蜜臀av无| 91久久精品国产一区二区成人| 天美传媒精品一区二区| 狂野欧美激情性xxxx在线观看| 观看美女的网站| 久久影院123| 国语对白做爰xxxⅹ性视频网站| 色哟哟·www| 日韩强制内射视频| 国产片内射在线| 国产淫语在线视频| 少妇精品久久久久久久| 国产又色又爽无遮挡免| 久久久久久久久久久免费av| 久久热精品热| 一区二区av电影网| 日本wwww免费看| 色哟哟·www| 美女cb高潮喷水在线观看| 亚洲成人一二三区av| 老女人水多毛片| 热99久久久久精品小说推荐| 久久精品熟女亚洲av麻豆精品| 欧美最新免费一区二区三区| av播播在线观看一区| 久久 成人 亚洲| 乱码一卡2卡4卡精品| 国产一区二区三区综合在线观看 | 成人综合一区亚洲| kizo精华| 久久女婷五月综合色啪小说| 精品视频人人做人人爽| 18+在线观看网站| 久久精品国产鲁丝片午夜精品| 久久久久久久久久久丰满| 国产精品偷伦视频观看了| 在线观看三级黄色| 亚洲美女视频黄频| 一边亲一边摸免费视频| 久久国产亚洲av麻豆专区| 亚洲精品成人av观看孕妇| 国产黄色视频一区二区在线观看| 最近中文字幕高清免费大全6| 日日撸夜夜添| 国产成人免费无遮挡视频| 国产成人aa在线观看| 十八禁高潮呻吟视频| 免费看不卡的av| 亚洲中文av在线| 婷婷成人精品国产| 国产日韩欧美亚洲二区| 自线自在国产av| 一级毛片aaaaaa免费看小| 中文字幕人妻丝袜制服| 亚洲天堂av无毛| 全区人妻精品视频| 91精品伊人久久大香线蕉| 天堂8中文在线网| 一级毛片我不卡| 一二三四中文在线观看免费高清| a 毛片基地| 欧美亚洲日本最大视频资源| 成人二区视频| 在线观看一区二区三区激情| 精品人妻一区二区三区麻豆| 另类亚洲欧美激情| 亚洲不卡免费看| 国产伦理片在线播放av一区| 亚洲精品aⅴ在线观看| 久久久久精品久久久久真实原创| 久久午夜综合久久蜜桃| av电影中文网址| 高清av免费在线| 欧美少妇被猛烈插入视频| 亚洲精品国产av成人精品| 热re99久久精品国产66热6| 久久久久久人妻| 高清视频免费观看一区二区| 日日摸夜夜添夜夜添av毛片| 亚洲精品久久成人aⅴ小说 | 校园人妻丝袜中文字幕| 国产黄色视频一区二区在线观看| 三级国产精品片| 黄色怎么调成土黄色| 国产探花极品一区二区| 国产在线视频一区二区| 亚洲精品自拍成人| 亚洲精品色激情综合| 桃花免费在线播放| 久久午夜综合久久蜜桃| av女优亚洲男人天堂| 夜夜骑夜夜射夜夜干| 人妻系列 视频| 少妇人妻精品综合一区二区| 亚洲美女搞黄在线观看| 欧美日本中文国产一区发布| 精品少妇黑人巨大在线播放| 我的老师免费观看完整版| 国产一级毛片在线| 女性被躁到高潮视频| 极品少妇高潮喷水抽搐| 亚洲少妇的诱惑av| 精品亚洲乱码少妇综合久久| 欧美激情 高清一区二区三区| 人妻制服诱惑在线中文字幕| 天天躁夜夜躁狠狠久久av| 日本黄大片高清| 国产成人91sexporn| 日本爱情动作片www.在线观看| 国产av国产精品国产| 少妇 在线观看| 精品卡一卡二卡四卡免费| 中文字幕av电影在线播放| 视频中文字幕在线观看| 成人无遮挡网站| 国产免费视频播放在线视频| 日韩亚洲欧美综合| 国产无遮挡羞羞视频在线观看| 亚洲不卡免费看| 在线精品无人区一区二区三| 男女免费视频国产| 成人手机av| 亚洲av.av天堂| 99九九在线精品视频| 日日摸夜夜添夜夜爱| 99热全是精品| 高清毛片免费看| 中国美白少妇内射xxxbb| 久久精品人人爽人人爽视色| 久久精品国产亚洲av涩爱| 欧美日韩在线观看h| 七月丁香在线播放| 制服诱惑二区| 国产精品欧美亚洲77777| 日本欧美国产在线视频| 国产一区有黄有色的免费视频| 高清视频免费观看一区二区| 日韩av在线免费看完整版不卡| 国产一区二区在线观看日韩| 亚洲欧美中文字幕日韩二区| 成人午夜精彩视频在线观看| 久久人妻熟女aⅴ| 3wmmmm亚洲av在线观看| 制服丝袜香蕉在线| av播播在线观看一区| 亚洲综合色网址| 精品久久久久久久久av| 熟女人妻精品中文字幕| 两个人免费观看高清视频| 国产精品成人在线| 欧美bdsm另类| 丰满少妇做爰视频| 精品99又大又爽又粗少妇毛片| 全区人妻精品视频| 97超视频在线观看视频| 国产精品三级大全| 国产精品久久久久久精品电影小说| 蜜桃久久精品国产亚洲av| 中国国产av一级| 午夜福利影视在线免费观看| 简卡轻食公司| 国产精品一区二区在线不卡| 欧美亚洲 丝袜 人妻 在线| 日韩制服骚丝袜av| 一本大道久久a久久精品| 久久久国产欧美日韩av| av专区在线播放| 欧美精品一区二区免费开放| 91在线精品国自产拍蜜月| 老女人水多毛片| 久久久亚洲精品成人影院| 人人澡人人妻人| 亚洲国产精品国产精品| 久久人妻熟女aⅴ| 男女高潮啪啪啪动态图| 亚洲精品中文字幕在线视频| 亚洲国产日韩一区二区| 国产一区二区三区av在线| 王馨瑶露胸无遮挡在线观看| 夫妻午夜视频| 七月丁香在线播放| 日日啪夜夜爽| 最新中文字幕久久久久| 亚洲av综合色区一区| 婷婷色麻豆天堂久久| 我要看黄色一级片免费的| 美女脱内裤让男人舔精品视频| 久久久久人妻精品一区果冻| 久久婷婷青草| 国产亚洲欧美精品永久| 午夜福利视频在线观看免费| 欧美bdsm另类| 国产精品久久久久久精品古装| 亚洲第一区二区三区不卡| 一区二区三区免费毛片| 欧美人与善性xxx| 中国国产av一级| 久久久欧美国产精品| 高清av免费在线| 久久久久久久亚洲中文字幕| 色婷婷av一区二区三区视频| 18禁观看日本| 日韩伦理黄色片| 日本黄色片子视频| 久久久久视频综合| 国产视频内射| 欧美精品人与动牲交sv欧美| 一区二区三区免费毛片| 亚洲精品美女久久av网站| 99久久中文字幕三级久久日本| 久热久热在线精品观看| 在线观看免费日韩欧美大片 | 免费人成在线观看视频色| 成人国语在线视频| 国产亚洲精品第一综合不卡 | 在线观看三级黄色| 国产欧美日韩一区二区三区在线 | 最近的中文字幕免费完整| 91久久精品电影网| 亚洲成人手机| 国产免费现黄频在线看| 久久韩国三级中文字幕| 国产精品一区二区三区四区免费观看| 日本av手机在线免费观看| 国产探花极品一区二区| 亚洲色图综合在线观看| 伦精品一区二区三区| 亚洲欧洲精品一区二区精品久久久 | 日韩欧美一区视频在线观看| 色婷婷av一区二区三区视频| 热re99久久精品国产66热6| 如日韩欧美国产精品一区二区三区 | 国产精品国产三级国产专区5o| 亚洲精品乱码久久久久久按摩| 在线免费观看不下载黄p国产| 亚洲激情五月婷婷啪啪| 久久人人爽人人片av| 只有这里有精品99| 啦啦啦啦在线视频资源| 久久99蜜桃精品久久| 在线 av 中文字幕| 母亲3免费完整高清在线观看 | 熟女电影av网| 2021少妇久久久久久久久久久| 国产精品久久久久久久电影| 精品久久久精品久久久| 九九在线视频观看精品| 国产成人av激情在线播放 | 91aial.com中文字幕在线观看| 在线观看www视频免费| 亚洲美女黄色视频免费看| 啦啦啦啦在线视频资源| 国产一区二区三区av在线| 国产精品久久久久久久久免| 99热国产这里只有精品6| 乱人伦中国视频| 亚洲怡红院男人天堂| 亚洲av在线观看美女高潮| 色吧在线观看| 国产色婷婷99| 久久影院123| 日本黄色日本黄色录像| 国产精品久久久久久av不卡| 美女脱内裤让男人舔精品视频| 国内精品宾馆在线| 国产免费一级a男人的天堂| 国产精品一区www在线观看| av网站免费在线观看视频| 人妻系列 视频| 69精品国产乱码久久久| 亚洲精品国产av蜜桃| 草草在线视频免费看| 国产高清三级在线| 亚洲国产欧美在线一区| 日韩一区二区视频免费看| 免费av中文字幕在线| 国产精品国产av在线观看| 免费大片18禁| 大香蕉久久网| 亚洲精品国产av成人精品| 欧美精品一区二区大全| 在线免费观看不下载黄p国产| 只有这里有精品99| 另类精品久久| 日韩中字成人| 午夜福利视频在线观看免费| 亚洲av福利一区| 亚洲精品一二三| 国产成人免费无遮挡视频| 欧美激情极品国产一区二区三区 | xxxhd国产人妻xxx| 九草在线视频观看| 一区二区三区免费毛片| 大陆偷拍与自拍| 中文乱码字字幕精品一区二区三区| 男的添女的下面高潮视频| 国产精品一二三区在线看| 午夜av观看不卡| 如何舔出高潮| 成人国语在线视频| 久热这里只有精品99| 精品一区在线观看国产| 亚洲在久久综合| 男男h啪啪无遮挡| 久久国产精品大桥未久av| 中文乱码字字幕精品一区二区三区| av网站免费在线观看视频| 久久久精品区二区三区| 日韩中字成人| 成年av动漫网址| 欧美精品国产亚洲| 熟女av电影| 一区二区三区免费毛片| 亚洲精品aⅴ在线观看| 热re99久久精品国产66热6| 菩萨蛮人人尽说江南好唐韦庄| 午夜影院在线不卡| 一本—道久久a久久精品蜜桃钙片| 赤兔流量卡办理|