• <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.

    亚洲精品aⅴ在线观看| 乱码一卡2卡4卡精品| 91精品伊人久久大香线蕉| 少妇丰满av| .国产精品久久| 欧美丝袜亚洲另类| 色哟哟·www| 最近中文字幕2019免费版| av专区在线播放| 久久97久久精品| 亚洲电影在线观看av| 国产片特级美女逼逼视频| 国产成年人精品一区二区| 九九在线视频观看精品| 亚洲av在线观看美女高潮| 美女cb高潮喷水在线观看| 国产成年人精品一区二区| 欧美xxⅹ黑人| 国产亚洲精品av在线| 18禁动态无遮挡网站| 免费看美女性在线毛片视频| 国产熟女欧美一区二区| 伊人久久国产一区二区| 久久午夜福利片| 婷婷色综合大香蕉| 国产精品一区二区三区四区免费观看| 韩国av在线不卡| 舔av片在线| 深夜a级毛片| 国产女主播在线喷水免费视频网站 | 一级二级三级毛片免费看| 免费电影在线观看免费观看| 能在线免费观看的黄片| 听说在线观看完整版免费高清| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产精品嫩草影院av在线观看| 97超视频在线观看视频| 我的老师免费观看完整版| 日韩欧美一区视频在线观看 | 亚洲成人久久爱视频| 69人妻影院| 爱豆传媒免费全集在线观看| 汤姆久久久久久久影院中文字幕 | 在线观看av片永久免费下载| 久久这里有精品视频免费| 日本黄色片子视频| 久久久久九九精品影院| 亚洲欧洲日产国产| 日韩欧美精品v在线| 亚洲成色77777| 国产又色又爽无遮挡免| 亚洲精品国产av蜜桃| 久久久色成人| 啦啦啦中文免费视频观看日本| 精品午夜福利在线看| 亚洲精品自拍成人| 欧美日韩综合久久久久久| 日本熟妇午夜| 乱码一卡2卡4卡精品| 麻豆精品久久久久久蜜桃| 99久国产av精品国产电影| 我要看日韩黄色一级片| 五月伊人婷婷丁香| 中文字幕av成人在线电影| h日本视频在线播放| 国产色婷婷99| 三级国产精品片| 午夜老司机福利剧场| 男女国产视频网站| 日韩欧美三级三区| 五月伊人婷婷丁香| 亚洲av中文av极速乱| 国产有黄有色有爽视频| 成人亚洲精品av一区二区| 午夜福利视频精品| 午夜激情久久久久久久| 亚洲自拍偷在线| 国产午夜福利久久久久久| 小蜜桃在线观看免费完整版高清| 18禁在线无遮挡免费观看视频| 黄片wwwwww| 你懂的网址亚洲精品在线观看| 午夜免费观看性视频| 韩国高清视频一区二区三区| av国产久精品久网站免费入址| 亚洲精品成人久久久久久| 女人久久www免费人成看片| 亚洲国产精品国产精品| 日本欧美国产在线视频| 国国产精品蜜臀av免费| 亚洲激情五月婷婷啪啪| 免费av不卡在线播放| 真实男女啪啪啪动态图| 免费在线观看成人毛片| 亚洲四区av| 日韩成人av中文字幕在线观看| 亚洲av免费在线观看| 男人舔女人下体高潮全视频| 精品国产露脸久久av麻豆 | 成人午夜精彩视频在线观看| 欧美成人午夜免费资源| 麻豆成人午夜福利视频| videossex国产| 免费电影在线观看免费观看| 亚洲欧美精品专区久久| 国产白丝娇喘喷水9色精品| 久久鲁丝午夜福利片| 一级爰片在线观看| 晚上一个人看的免费电影| 秋霞在线观看毛片| 干丝袜人妻中文字幕| 国产69精品久久久久777片| 极品教师在线视频| 免费在线观看成人毛片| 熟女人妻精品中文字幕| 国产黄色视频一区二区在线观看| 久久国产乱子免费精品| 国产高清不卡午夜福利| 欧美激情久久久久久爽电影| 91精品一卡2卡3卡4卡| 夜夜爽夜夜爽视频| 精品人妻一区二区三区麻豆| 99久久人妻综合| 久久精品国产亚洲网站| 日韩亚洲欧美综合| 美女内射精品一级片tv| av在线天堂中文字幕| 国产精品av视频在线免费观看| 精品不卡国产一区二区三区| 色尼玛亚洲综合影院| 日韩 亚洲 欧美在线| 免费看美女性在线毛片视频| 嫩草影院精品99| 狂野欧美激情性xxxx在线观看| 午夜激情久久久久久久| 亚洲国产精品专区欧美| 亚洲欧洲国产日韩| 亚洲成人中文字幕在线播放| 少妇人妻一区二区三区视频| 成人毛片a级毛片在线播放| 欧美日韩综合久久久久久| 看十八女毛片水多多多| 国产精品一区二区在线观看99 | 寂寞人妻少妇视频99o| 黄色配什么色好看| 久久精品国产亚洲av涩爱| 成人欧美大片| 在线观看人妻少妇| 99热这里只有精品一区| 午夜激情欧美在线| 精品人妻熟女av久视频| 在线 av 中文字幕| 成年人午夜在线观看视频 | 高清视频免费观看一区二区 | 高清日韩中文字幕在线| 亚洲av成人精品一二三区| 欧美日韩在线观看h| 少妇人妻一区二区三区视频| ponron亚洲| 1000部很黄的大片| 免费观看在线日韩| 人妻系列 视频| 国产精品一区www在线观看| 毛片女人毛片| 99久久精品国产国产毛片| 2021少妇久久久久久久久久久| 最近手机中文字幕大全| 欧美成人一区二区免费高清观看| 麻豆成人午夜福利视频| 在线 av 中文字幕| 在线观看免费高清a一片| 国国产精品蜜臀av免费| 国产黄a三级三级三级人| 国产精品一区二区三区四区久久| 日本wwww免费看| 国产黄频视频在线观看| 国产有黄有色有爽视频| 精品酒店卫生间| 久久精品熟女亚洲av麻豆精品 | 真实男女啪啪啪动态图| 麻豆乱淫一区二区| 亚洲四区av| 日日摸夜夜添夜夜爱| 嫩草影院精品99| 免费观看性生交大片5| 国产视频首页在线观看| 亚洲最大成人av| 69人妻影院| 丝袜美腿在线中文| 嫩草影院新地址| 久久久精品欧美日韩精品| 亚洲欧洲日产国产| 国产欧美另类精品又又久久亚洲欧美| xxx大片免费视频| 国产欧美另类精品又又久久亚洲欧美| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 久久久久久久久久久丰满| 欧美高清性xxxxhd video| 久久人人爽人人片av| a级毛色黄片| 天堂俺去俺来也www色官网 | 欧美日韩综合久久久久久| 久久久久久久久久人人人人人人| 成人亚洲欧美一区二区av| 国产精品不卡视频一区二区| 久久鲁丝午夜福利片| 伦理电影大哥的女人| 成人高潮视频无遮挡免费网站| 亚洲无线观看免费| 看十八女毛片水多多多| 日本-黄色视频高清免费观看| 欧美精品一区二区大全| 亚洲精品亚洲一区二区| 国产亚洲5aaaaa淫片| 久久精品国产自在天天线| 亚洲熟妇中文字幕五十中出| 深爱激情五月婷婷| 色视频www国产| 国产高清三级在线| 欧美高清成人免费视频www| 特级一级黄色大片| 国产高潮美女av| 久久久久久久久大av| 一级毛片久久久久久久久女| 99热全是精品| 日韩一区二区视频免费看| 亚洲电影在线观看av| 精品一区二区免费观看| 超碰av人人做人人爽久久| 乱系列少妇在线播放| 真实男女啪啪啪动态图| 九草在线视频观看| 国产免费福利视频在线观看| 午夜日本视频在线| 亚洲欧洲国产日韩| 一区二区三区四区激情视频| 一本一本综合久久| 麻豆av噜噜一区二区三区| 精品久久久久久成人av| 男女那种视频在线观看| 久久久久久久久久成人| 精品熟女少妇av免费看| 2021天堂中文幕一二区在线观| 99热网站在线观看| 91久久精品国产一区二区三区| 91久久精品国产一区二区成人| 老司机影院成人| 男人狂女人下面高潮的视频| 一级爰片在线观看| 日本熟妇午夜| 国产精品综合久久久久久久免费| 国产不卡一卡二| 色哟哟·www| 大片免费播放器 马上看| 国产黄色视频一区二区在线观看| 狠狠精品人妻久久久久久综合| 亚洲丝袜综合中文字幕| 亚洲av二区三区四区| 午夜日本视频在线| 一级毛片 在线播放| 亚洲av电影在线观看一区二区三区 | 免费大片黄手机在线观看| av在线蜜桃| eeuss影院久久| 国产成人精品婷婷| 丝袜美腿在线中文| 国产黄色视频一区二区在线观看| 午夜精品国产一区二区电影 | 男女视频在线观看网站免费| 亚洲国产精品成人久久小说| 国产精品久久久久久久久免| 一级毛片aaaaaa免费看小| 国内精品一区二区在线观看| 国产永久视频网站| 日韩伦理黄色片| 搡老乐熟女国产| 亚洲av免费高清在线观看| 中文字幕亚洲精品专区| 亚洲国产日韩欧美精品在线观看| 你懂的网址亚洲精品在线观看| 欧美日韩精品成人综合77777| 久久久久久国产a免费观看| 丰满少妇做爰视频| 蜜桃亚洲精品一区二区三区| 日韩电影二区| 欧美日韩一区二区视频在线观看视频在线 | 成人无遮挡网站| 国产精品久久久久久久久免| 欧美高清成人免费视频www| 精品一区二区三区人妻视频| 热99在线观看视频| 18禁在线播放成人免费| 亚洲人与动物交配视频| 日本一本二区三区精品| 欧美最新免费一区二区三区| 精品人妻视频免费看| 欧美另类一区| 亚洲最大成人手机在线| 欧美高清性xxxxhd video| 麻豆成人av视频| 男女啪啪激烈高潮av片| 国产精品一区二区三区四区久久| 精品久久久精品久久久| 日本三级黄在线观看| 免费看光身美女| 哪个播放器可以免费观看大片| 最近视频中文字幕2019在线8| av女优亚洲男人天堂| 丝袜喷水一区| 亚洲成人av在线免费| 美女主播在线视频| 有码 亚洲区| 亚洲精品日本国产第一区| 欧美97在线视频| 成年女人看的毛片在线观看| 99热网站在线观看| 毛片一级片免费看久久久久| 精品久久久久久成人av| 久久精品国产亚洲av涩爱| 国产黄频视频在线观看| 国产精品日韩av在线免费观看| 国产老妇伦熟女老妇高清| 十八禁国产超污无遮挡网站| 大陆偷拍与自拍| 久久久久久伊人网av| 国模一区二区三区四区视频| 免费黄网站久久成人精品| 国产白丝娇喘喷水9色精品| 黄片wwwwww| 国产成人a∨麻豆精品| 蜜桃久久精品国产亚洲av| 久久鲁丝午夜福利片| 精华霜和精华液先用哪个| 视频中文字幕在线观看| 在线免费观看不下载黄p国产| 久久久精品免费免费高清| 91久久精品电影网| 我要看日韩黄色一级片| 少妇人妻精品综合一区二区| 日韩欧美一区视频在线观看 | 女人久久www免费人成看片| 2021天堂中文幕一二区在线观| 国产视频首页在线观看| 熟女人妻精品中文字幕| 可以在线观看毛片的网站| av在线蜜桃| 色哟哟·www| 最近的中文字幕免费完整| 日韩欧美一区视频在线观看 | 国产亚洲5aaaaa淫片| 亚洲国产欧美在线一区| av国产久精品久网站免费入址| 精品久久久久久成人av| 欧美不卡视频在线免费观看| 日本wwww免费看| 丰满少妇做爰视频| 熟妇人妻不卡中文字幕| 在线播放无遮挡| 国产永久视频网站| 99九九线精品视频在线观看视频| 久久人人爽人人片av| 人人妻人人看人人澡| ponron亚洲| 久久久久久久午夜电影| 在线 av 中文字幕| 成人美女网站在线观看视频| 国产免费福利视频在线观看| 哪个播放器可以免费观看大片| 搡老乐熟女国产| 色视频www国产| 日韩在线高清观看一区二区三区| 亚洲av日韩在线播放| 看非洲黑人一级黄片| 色播亚洲综合网| 亚洲久久久久久中文字幕| 精品少妇黑人巨大在线播放| 亚洲av国产av综合av卡| 国产亚洲午夜精品一区二区久久 | 最近中文字幕高清免费大全6| 最近的中文字幕免费完整| 国产伦理片在线播放av一区| 国产真实伦视频高清在线观看| 久久99热这里只频精品6学生| 一级av片app| 在线a可以看的网站| 免费电影在线观看免费观看| 亚洲精品国产av成人精品| 天堂网av新在线| 国产精品综合久久久久久久免费| 有码 亚洲区| 午夜爱爱视频在线播放| 午夜免费观看性视频| 久久热精品热| 亚洲欧洲日产国产| 亚洲精品第二区| 天堂av国产一区二区熟女人妻| 男人舔奶头视频| 天美传媒精品一区二区| 国产91av在线免费观看| 免费av不卡在线播放| 我的女老师完整版在线观看| 国产老妇伦熟女老妇高清| 在线播放无遮挡| 3wmmmm亚洲av在线观看| 欧美xxxx黑人xx丫x性爽| 亚洲人成网站高清观看| 2022亚洲国产成人精品| 一二三四中文在线观看免费高清| 又黄又爽又刺激的免费视频.| 日韩av在线免费看完整版不卡| av在线播放精品| 在线观看人妻少妇| 中文天堂在线官网| 最近视频中文字幕2019在线8| 欧美 日韩 精品 国产| 黄色一级大片看看| 免费av毛片视频| 成人av在线播放网站| 80岁老熟妇乱子伦牲交| 黄片无遮挡物在线观看| 大片免费播放器 马上看| av在线老鸭窝| 亚洲av国产av综合av卡| 成人午夜高清在线视频| 午夜免费观看性视频| 亚洲精品成人久久久久久| 能在线免费看毛片的网站| 国产欧美日韩精品一区二区| 亚洲欧洲国产日韩| 男人爽女人下面视频在线观看| 国产在线男女| 精品人妻视频免费看| 欧美成人a在线观看| 欧美97在线视频| av国产免费在线观看| 男人爽女人下面视频在线观看| 久久久久久久久久成人| 免费电影在线观看免费观看| 最近2019中文字幕mv第一页| 亚洲最大成人中文| 欧美精品国产亚洲| 女人被狂操c到高潮| 国产精品日韩av在线免费观看| 精品一区二区三卡| 99久久中文字幕三级久久日本| 国产单亲对白刺激| 91久久精品电影网| 在现免费观看毛片| 夫妻性生交免费视频一级片| 久久久久久国产a免费观看| 99热这里只有精品一区| 久久精品人妻少妇| 亚洲av中文字字幕乱码综合| 欧美三级亚洲精品| 亚洲伊人久久精品综合| 天天躁日日操中文字幕| 黄色欧美视频在线观看| 天堂中文最新版在线下载 | 身体一侧抽搐| 乱人视频在线观看| 国产高清三级在线| 亚洲国产精品sss在线观看| 美女cb高潮喷水在线观看| 91狼人影院| 日韩大片免费观看网站| 99热这里只有是精品在线观看| 亚洲成人av在线免费| 欧美性猛交╳xxx乱大交人| 国产 一区 欧美 日韩| 国产av在哪里看| 中文乱码字字幕精品一区二区三区 | 精品熟女少妇av免费看| av在线天堂中文字幕| 国产欧美日韩精品一区二区| 国产 亚洲一区二区三区 | 三级毛片av免费| 亚洲av二区三区四区| 久久精品综合一区二区三区| 高清在线视频一区二区三区| 久久99蜜桃精品久久| 狂野欧美白嫩少妇大欣赏| 亚洲色图av天堂| 免费大片18禁| 久久久精品免费免费高清| xxx大片免费视频| av在线老鸭窝| 美女国产视频在线观看| 精品99又大又爽又粗少妇毛片| 亚洲精品成人久久久久久| 女人被狂操c到高潮| 久久99精品国语久久久| 中文乱码字字幕精品一区二区三区 | 亚洲精品视频女| 高清欧美精品videossex| 欧美不卡视频在线免费观看| 91久久精品国产一区二区三区| 国产 一区 欧美 日韩| 男人舔女人下体高潮全视频| 97精品久久久久久久久久精品| 你懂的网址亚洲精品在线观看| 日日干狠狠操夜夜爽| 国产伦精品一区二区三区四那| 久久99热这里只有精品18| 麻豆av噜噜一区二区三区| 91久久精品电影网| 一级毛片我不卡| 九草在线视频观看| 国产高清国产精品国产三级 | 国产成人91sexporn| 免费看日本二区| 亚洲精品亚洲一区二区| 久久国内精品自在自线图片| 欧美不卡视频在线免费观看| 国产精品精品国产色婷婷| 精品国内亚洲2022精品成人| 亚洲伊人久久精品综合| 女的被弄到高潮叫床怎么办| 亚洲精品乱码久久久久久按摩| 18禁裸乳无遮挡免费网站照片| 女人久久www免费人成看片| 欧美性猛交╳xxx乱大交人| 精品一区二区三卡| 久久久久久久久大av| 亚洲国产成人一精品久久久| 免费看a级黄色片| 久久久色成人| 欧美日韩综合久久久久久| 嫩草影院新地址| 国产亚洲精品av在线| 日本猛色少妇xxxxx猛交久久| av黄色大香蕉| 亚洲av成人精品一二三区| 久久草成人影院| 亚洲精品久久久久久婷婷小说| 国产精品嫩草影院av在线观看| 国产精品一区二区在线观看99 | 一级黄片播放器| 三级国产精品欧美在线观看| 欧美xxxx黑人xx丫x性爽| 免费看光身美女| 成人漫画全彩无遮挡| 国产在线男女| 亚洲欧洲日产国产| 精品久久久噜噜| 美女国产视频在线观看| 中文字幕制服av| 哪个播放器可以免费观看大片| 亚洲精品第二区| 久久久久久九九精品二区国产| 春色校园在线视频观看| 欧美最新免费一区二区三区| av在线播放精品| 麻豆成人av视频| 国产高清国产精品国产三级 | 国产黄a三级三级三级人| 中文在线观看免费www的网站| 亚洲美女视频黄频| 国产久久久一区二区三区| 亚洲精品成人久久久久久| 麻豆av噜噜一区二区三区| 亚洲欧美精品专区久久| 麻豆乱淫一区二区| 国产乱人视频| 久99久视频精品免费| 午夜福利网站1000一区二区三区| 18禁动态无遮挡网站| 免费看a级黄色片| 国产精品国产三级国产专区5o| 久久久亚洲精品成人影院| 久99久视频精品免费| 一区二区三区高清视频在线| 在现免费观看毛片| 国产免费福利视频在线观看| 毛片女人毛片| 国产免费一级a男人的天堂| av国产久精品久网站免费入址| 伦理电影大哥的女人| 亚洲欧美日韩卡通动漫| 黄色配什么色好看| 日日干狠狠操夜夜爽| 看免费成人av毛片| 日日干狠狠操夜夜爽| www.av在线官网国产| 一夜夜www| 黄色配什么色好看| 91精品伊人久久大香线蕉| 国产成人精品福利久久| 亚洲国产精品成人综合色| 午夜老司机福利剧场| 欧美区成人在线视频| 最近手机中文字幕大全| 精品国产露脸久久av麻豆 | 亚洲精品成人av观看孕妇| 男女啪啪激烈高潮av片| 狂野欧美白嫩少妇大欣赏| 久久久亚洲精品成人影院| 在线免费观看不下载黄p国产| 亚洲精品自拍成人| 99热6这里只有精品| 国产极品天堂在线| 最近的中文字幕免费完整| 国产极品天堂在线| 干丝袜人妻中文字幕| 日本午夜av视频| 欧美极品一区二区三区四区| 日韩精品有码人妻一区| 伊人久久国产一区二区| 国产精品国产三级国产av玫瑰| 天堂网av新在线| 一级毛片电影观看| 国产成人一区二区在线| 亚洲欧美日韩卡通动漫| 肉色欧美久久久久久久蜜桃 | 午夜老司机福利剧场| 国产免费一级a男人的天堂|