999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

4,5-二(1H-5-四唑基)-1,2,3-三唑及其含能鹽的合成與表征

2017-05-18 07:45:24張艷芳熊華林林秋漢程廣斌楊紅偉
火炸藥學報 2017年2期

張艷芳,熊華林,林秋漢,程廣斌,楊紅偉

(南京理工大學化工學院,江蘇 南京 210094)

4,5-二(1H-5-四唑基)-1,2,3-三唑及其含能鹽的合成與表征

張艷芳,熊華林,林秋漢,程廣斌,楊紅偉

(南京理工大學化工學院,江蘇 南京 210094)

以4,5-二氰基-1,2,3-三唑為原料,經過重氮化偶聯反應和[2+3]偶氮環加成反應合成了4,5-二(1H-5-四唑基)-1,2,3-三唑,并經過復分解反應合成了7種高氮含能鹽。采用IR、1H NMR、13C NMR對化合物的結構進行了表征;采用DSC法研究了它們的熱行為;基于B3LYP/6-311G**方法計算了各化合物的理論密度、標準生成焓,并通過EXPLO5預估了各化合物的爆速和爆壓。結果表明,8種化合物均具有較好的熱穩定性;生成焓為647.4~1579.4kJ/mol;爆速為7652~8389m/s;爆壓除化合物2外均高于20 GPa;除化合物1、7、8外,其他化合物標準狀況下的氣體產生量均大于800 L/kg,是潛在的氣體發生劑。

有機化學;高氮化合物;四唑;三唑;含能離子鹽

引 言

高氮含能化合物作為一種新型含能材料受到世界各國的關注[1-4],具有良好的應用前景[5-6]。高氮含能化合物具有高的生成焓,爆炸產生大量的N2和H2O,不會對環境造成污染。唑類化合物由于氮含量高、熱穩定性好等優點已在高氮含能化合物合成領域得到廣泛研究[7]。唑類高氮含能離子鹽是一類獨特的高氮含量的含能材料,含能離子鹽與非離子分子相比具有較低的蒸汽壓而更具有優勢[8]。Thomas M. Klap?tke等[9]首先合成1,5-聯四唑,由于氮含量高被認為是推進劑的潛在組分;翟連杰、樊學忠等[10]報道了多種合成4,5-二(1H-5-四唑基)-1,2,3-三唑的方法;隨后Thomas M. Klap?tke等[11]報道了4,5-二(1H-5-四唑基)-1,2,3-三唑及其銨鹽、羥胺鹽、氨基胍鹽的合成,并預估了各化合物的爆轟性能。其爆速介于7785~9022m/s,各化合物還具有較高的生成焓(591~846kJ/mol)和熱分解溫度(160~301°C)。任曉雪[12]研究表明高氮化合物是具有良好應用前景的煙火型氣體發生劑原料。由于四唑類衍生物氮含量和能量均高、毒性低,因而成為高氮化合物的典型代表。

為了進一步擴展唑類化合物作為氣體發生劑的應用潛力,本研究在前人工作的基礎上合成了4,5-二(1H-5-四唑基)-1,2,3-三唑的一系列高氮含能鹽。該系列化合物氮含量高、生成焓高、氣體產生量大且不對環境造成污染,符合當代氣體發生劑、低特征信號推進劑、低煙或無煙煙火劑等所追求的高性能、綠色環保的特點。

1 實 驗

1.1 試劑與儀器

鹽酸,上海凌峰化學試劑有限公司;氯化鋅,天津市瑞金特化學品有限公司;亞硝酸鈉、乙醚,國藥集團化學試劑有限公司;硝酸銀,上海麥克林生化科技有限公司;水合肼(質量分數80%),上海展云化工有限公司;3-氨基-1,2,4-三氮唑,上海晶純生化科技股份有限公司;疊氮化鈉,山東西亞化學股份有限公司;乙醇、氨水,上海泰坦科技股份有限公司;2,3-二氨基-2-丁烯二腈、羥胺(質量分數50%)、碳酸二酰肼、1,2,4-三氮唑,薩恩化學技術(上海)有限公司。以上試劑均為分析純。

BOMEM MB系列154S系列傅里葉變換紅外光譜儀,加拿大Abb Bomem公司;Bruker Avance 300M數字化核磁共振譜儀,瑞士Bruker公司;TGA/ADTA851e型熱重分析儀、DSC823e型差示掃描量熱儀,Mettler-Toledo國際貿易有限公司;ZF-1型三用紫外分析儀,上海市安亭電子儀器廠。

1.2 實驗過程

1.2.1 8種高氮含能化合物的合成路線如圖1所示。

圖1 化合物1~8的合成路線Fig.1 Synthetic route of compounds 1-8

1.2.2 4,5-二氰基-1,2,3-三唑的合成[13]

將10.81 g(0.10mol)2,3-二氨基-2-丁烯二腈加入到100mL濃度為1mol/L的鹽酸溶液中。將6.89g(0.10mol)亞硝酸鈉溶于125mL水中,充分攪拌至完全溶解。在低溫條件下,將上述亞硝酸鈉的水溶液緩慢滴入反應液中,滴加過程中保持反應液溫度介于0~5℃。加料完畢后在該溫度下反應0.5h,常溫反應1h。經抽濾、乙醚萃取、干燥并旋干乙醚溶液得到10.20g橘色粉末狀固體,收率85.61%。

1.2.3 4,5-二(1H-5-四唑基)-1,2,3-三唑(化合物1)的合成[14]

向250mL單口瓶中加入2.38g(37mmol)疊氮化鈉、5.00g(37mmol)氯化鋅、100mL乙醇,充分攪拌。緩慢加入2.00g(17mmol) 4,5-二氰基-1,2,3-三唑回流24h,TLC監測反應終點。反應結束后,將反應液倒入等體積的冰水中,濃鹽酸調節pH值為1。攪拌2h ,經過濾、洗滌、干燥,得到3.01g白色塊狀固體,收率86.49%。

1H NMR ( DMSO-d6,300 M ),δ:4.27 (s);13C NMR (DMSO-d6,75 M ),δ: 148.37,132.01;IR(KBr), ν(cm-1):3343,3088,2469,1909,1620,1027,980, 886;Tdonset=296.14℃。元素分析(C4H3N11,%):計算值,C 23.42,H 1.47,N 75.11;實測值, C 23.37,H 1.50,N 75.13。

1.2.4 4,5-二(1H-5-四唑基)-1,2,3-三唑鹽的合成

(1) 4,5-二(1H-5-四唑基)-1,2,3-三唑銨鹽(化合物2)

向25mL 的單口瓶中加入0.20g(0.99mmol) 4,5-二(1H-5-四唑基)-1,2,3-三唑,少量水,升溫至50 ℃,緩慢加入0.27mL(1.96mmol)氨水,攪拌2h,過濾、洗滌、干燥得到0.21g白色粉末狀固體,收率90.05%。

13C NMR (DMSO-d6,75M),δ:152.92,133.15;IR(KBr),ν(cm-1):3175,3055,2832,1882,1445,1412,994,576;Td,onset=297.75℃。元素分析(C4H9N13,%):計算值,C 20.08,H 3.77,N 76.15;實測值,C 20.12,H 3.81,N 76.07。

(2) 4,5-二(1H-5-四唑基)-1,2,3-三唑肼鹽(化合物3)

向25mL 的單口瓶中加入0.20g(0.99mmol) 4,5-二(1H-5-四唑基)-1,2,3-三唑和少量乙醇。緩慢滴入0.13mL(2.14mmol) 質量分數80%的水合肼溶液,常溫攪拌2h,抽濾、洗滌、干燥得到0.23g白色粉末狀固體,收率88.47%。

1H NMR ( DMSO-d6,300M ),δ:5.62 (s);13C NMR (DMSO-d6,75M),δ:153.35,132.36;IR(KBr),ν(cm-1):3343,2952,2644,1533,1135,927,765,442;Td,onset=301.14℃。元素分析(C4H11N15,%):計算值,C 17.84,H 4.09,N 70.07;實測值,C 17.90,H 4.16,N 77.94。

(3) 4,5-二(1H-5-四唑基)-1,2,3-三唑羥胺鹽(化合物4)

向25mL的單口瓶中加入0.20g(0.99mmol) 4,5-二(1H-5-四唑基)-1,2,3-三唑和少量乙醇。緩慢滴入0.14mL(2.19mmol) 質量分數50%的羥胺水溶液,常溫攪拌2h,抽濾、洗滌、干燥得到0.23g白色粉末狀固體,收率87.71%。

1H NMR ( DMSO-d6,300M ),δ:9.42 (s),3.14 (s);13C NMR (DMSO-d6,75M),δ:150.43, 130.57;IR(KBr),ν(cm-1):3074,2832,2664,1405,994,887,576,489;Td,onset=287.85℃。元素分析(C4H9N13O2,%):計算值,C 1.48,H 3.32,N 67.16,;實測值, C 1.52,H 3.36,N 67.25。

(4) 4,5-二(1H-5-四唑基)-1,2,3-三唑碳酸二酰肼鹽(化合物5)

向25mL的單口瓶中加入0.20g(0.99mmol)4,5-二(1H-5-四唑基)-1,2,3-三唑和少量乙醇。將0.20g(2.22mmol)碳酸二酰肼溶于少量水,緩慢滴入到反應液中,常溫攪拌2h,過濾、洗滌、干燥得到0.32g白色粉末狀固體,收率85.01%。

1H NMR ( DMSO-d6,300 M ),δ:5.01 (s),3.41 (s);13C NMR ( DMSO-d6,75 M ),δ:159.77,150.03,130.60;IR(KBr),ν(cm-1) :3337,3000,2832,2677,1587,1439,987,576;Td,onset=278.22℃。元素分析(C6H15N19O2,%):計算值,C 18.70,H 3.90,N 69.09;實測值,C 18.72,H 3.95,N 68.98。

(5) 4,5-二(1H-5-四唑基)-1,2,3-三唑氨基胍鹽(化合物6)

向25mL的單口瓶中加入0.20g(0.84mmol) 4,5-二(1H-5-四唑基)-1,2,3-三唑銨鹽和少量水,加熱至50 ℃,銨鹽逐漸溶解。將0.40g(2.35mmol)硝酸銀溶于水,緩慢滴入反應液中攪拌2h,抽濾、洗滌、干燥得到0.34g白色固體,收率84.23%。

向25mL 的單口瓶中加入0.20g(0.48mmol) 4,5-二(1H-5-四唑基)-1,2,3-三唑銀鹽和少量水,加熱至回流。0.10g(1.05mmol)鹽酸胍溶于適量水中,緩慢滴入反應液,回流、抽濾、洗滌、干燥得到0.15g白色固體,收率88.93%。

1H NMR ( DMSO-d6,300 M ),δ:8.23 (br),4.76 (s),3.19 (s);13C NMR ( DMSO-d6,75 M ),δ:159.48,156.44,133.79;IR(KBr),ν(cm-1) :3357,3182,3013,1680,1621,1128,980,698;Td,onset=274.41℃。元素分析(C6H15N19,%):計算值,C 20.40,H 4.28,N 75.32;實測值, C 20.42,H 4.30,N 75.28。

(6) 4,5-二(1H-5-四唑基)-1,2,3-三唑1,2,4-三氮唑鹽(化合物7)

向25mL 的單口瓶中加入0.20g(0.99mmol) 4,5-二(1H-5-四唑基)-1,2,3-三唑和10mL甲醇。緩慢加入0.15g(2.17mmol) 1,2,4-三氮唑回流,過濾、洗滌、干燥得到 0.29g白色粉末狀固體,收率86.65%。

1H NMR ( DMSO-d6,300 M ),δ:9.41 (s),8.86 (s),8.16 (s);13C NMR ( DMSO-d6,75 M ),δ:148.33,143.10,132.07;IR(KBr),ν(cm-1):3512,3243,3082,2462,1903,1028,960,900;Td,onset=300.05℃。元素分析(C8H9N17,%):計算值,C 28.07,H 2.62,N 69.39;實測值, C 28.07,H 2.55,N 69.38。

(7) 4,5-二(1H-5-四唑基)-1,2,3-三唑3-氨基-1,2,4-三氮唑鹽(化合物8)

向25mL的單口瓶中加入0.20g(0.99mmol) 4,5-二(1H-5-四唑基)-1,2,3-三唑和10mL甲醇。緩慢加入0.17g(2.02mmol) 3-氨基-1,2,4-三氮唑回流,過濾、洗滌、干燥得到0.27g白色粉末狀固體,收率74.19%。

1H NMR ( DMSO-d6,300 M ),δ:9.64 (tr),7.91 (s),3.14 (s);13C NMR ( DMSO-d6, 75 M ),δ:154.52, 149.87, 143.61, 130.69; IR(KBr),ν(cm-1) :3472,3108,2759,1680,1250,953,725,583;Td,onset=309.06℃。元素分析(C8H11N19,%):計算值,C 25.74,H 2.97,N 71.29;實測值, C 25.80,H 2.89,N 71.31。

2 結果與討論

2.1 熱力學性能分析

在氮氣流量30mL/min、升溫速率10℃/min的條件下,對4,5-二(1H-5-四唑基)-1,2,3-三唑及其含能鹽進行DSC分析,結果見圖2。

圖2 化合物1~8的DSC曲線Fig.2 DSC curves of compounds 1-8

由圖2可知,化合物1~8的熱分解溫度介于274.41~309.06℃,均高于RDX(210℃),表明化合物具有良好的熱穩定性。其中化合物3、化合物7、化合物8的熱分解溫度均高于300℃,優于TNT(295℃)。該系列化合物熱穩定性較高是由于這類高氮化合物所含的氮原子電負性高,氮雜環體系形成類似苯環結構的大π鍵體系,而大π鍵體系的存在通常使熱穩定性較高。化合物1在109.48℃時開始吸熱,到126.44℃出現一個峰值,為其吸熱峰。從296.14℃開始,經歷一個緩慢的放熱過程,在301.83℃時放熱速率達到最大。化合物2從297.75℃開始放熱,在316.20℃放熱速率達到最大。化合物3沒有經過熔化吸熱的過程,而是直接發生劇烈的分解反應釋放大量的熱,分解峰呈尖銳狀。從301.14℃開始放熱,在304.60℃放熱速率達到最大。化合物4~7均有一個很小的吸收峰。分別從287.85、278.22、274.41、300.05℃開始放熱,分別在304.03、288.94、317.97、306.97℃放熱速率達到最大。化合物8在257.63℃時開始吸熱,到271.57℃出現吸熱峰的峰值。從309.06℃開始放熱,在319.95℃時放熱速率達到最大。

2.2 生成熱和爆轟性能的理論計算

運用Gaussian09程序,在B3LYP/6-311+G**水平下對化合物結構進行優化,得到其在勢能面上的穩定結構,經振動分析無虛頻,并對各化合物的性質進行計算。所有化合物的生成熱均基于設計等鍵反應和運用B3LYP/6-311+G**方法計算,所設計的等鍵反應如圖3所示。各化合物總能量(E0)、零點能(ZPE)、溫度校正系數(ΔHT)以及生成熱(ΔfHm)見表1。

(a) 化合物1

(b) 化合物2~8陰離子圖3 化合物1和化合物2~8陰離子的等鍵方程Fig.3 Isodesmic reaction for compound 1 and anions of compounds 2-8

表1 各化合物的總能量、零點能、溫度校正系數以及生成熱

Table 1E0,ZPE,HTand ΔfHmof different compound

化合物E0/(a.u.)ZPE/(kJ·mol-1)ΔHT/(kJ·mol-1)ΔfHm/(kJ·mol-1)1-755.43219.1129.10774.862-56.92124.719.98636.743-112.24166.0212.04789.464-132.09137.0611.20694.985-336.31269.1919.25790.496-261.14255.6117.03745.977-242.59179.3912.861020.998-298.01223.8916.17937.07

各化合物的理論密度是通過分子摩爾質量(M)除以平均摩爾體積(V)的方法計算。對各分子的0.001e/bohr3等電子密度面所包圍的體積空間,用Monte-Carlo方法求V,進而求得各分子的密度(ρ)。根據各化合物理論計算所得的生成焓和密度,利用EXPLO5程序[15-16]對關鍵爆轟參數爆速(D)和爆壓(p)進行計算,結果見表2。由表2可知,化合物1~8的氮含量高于60 %。除化合物2和化合物6外,其他化合物的爆熱均高于TNT。化合物1~8的生成焓均為正,介于647.4~1597.4kJ/mol。這是因為其結構中含有大量的N-N鍵及N=N鍵,含能化合物骨架中N-N、C=N、N=N鍵的存在對化合物的生成焓均有較大的貢獻,有利于提高分子的爆熱。此外,除化合物1和化合物7外的其他化合物在標準狀況下的氣體產生量均高于RDX(739L/kg),其中化合物5的氣體產生量為854L/kg。各化合物的爆壓介于19.4~25.7GPa,爆速為7652~8389m/s,優于TNT的爆轟性能(爆速6881m/s,爆壓19.5GPa)。除化合物2、化合物7外的其他化合物的爆速均高于8000m/s,優于TNT。

表2 化合物1~8物理性質及爆轟性能

注:M為分子質量;w(N)為氮含量;Ω為氧平衡;Tdec為熱分解溫度;ρ為密度;ΔfHm為標準摩爾生成焓;D為爆速;p為爆壓;Q為放熱量;V0為標況下的氣體產生量。

3 結 論

(1)以4,5-二氰基-1,2,3-三唑為原料,經過重氮化偶聯反應、[2+3]偶氮環加成反應合成了4,5-二(1H-5-四唑基)-1,2,3-三唑,并經過復分解反應合成了一系列高氮含能鹽。DSC分析表明,化合物1~8的熱分解溫度介于274.41~309.06℃,表明所合成的化合物均具有良好的熱穩定性。

(2)由于大量N-N鍵及N=N鍵的存在,化合物1~8均具有較高的正生成焓及較好的爆轟性能。各化合物的爆速均大于TNT,除化合物2爆壓與TNT相當外,其他化合物的爆壓均大于TNT,綜合比較各化合物,化合物5的性能較佳。

(3)化合物1~8氮含量高、生成焓高、氣體產生量大且不對環境造成污染,是潛在的氣體發生劑。

[1]HervéG,RousselC,GraindorgeH.Selectivepreparationof3,4,5-trinitro-1H-pyraz-ole:astableall-carbon-nitratedarene[J].AngewandteChemieInternationalEdition, 2010, 49(18):3177 -3181.

[2]LiS,WangY,QiC,etal. 3Denergeticmetal-organicframeworks:synthesisandpropertiesofhighenergymaterials[J].AngewandteChemieInternationalEdition, 2013, 52(52): 14031-14035.

[3]ZhangQ,ShreeveJM.Growingcatenatednitrogenatomchains[J].AngewandteChemieInternationalEdition, 2013, 52(34): 8792-8794.

[4]ZhangM,EatonPE,GilardiR.Hepta-andoctanitrocubanes[J].AngewandteChemieInternationalEdition, 2000, 39(2): 401-404.

[5]MandalAK,PantCS,KasarSM,etal.ProcessoptimizationforsynthesisofCL-20 [J].JournalofEnergeticMaterials, 2009, 27(4):231-246.

[6] 畢衛宇, 劉愆, 張志忠, 等.TNAZ的硝基甲烷法合成研究 [J]. 火炸藥學報, 2002, 25(4):31-32.BIWei-yu,LIUQian,ZHANGZhi-zhong,etal.SynthesisofTNAZvianitromethane[J].ChineseJournalExplosives&Propellants(HuozhayaoXuebao), 2002, 25(4):31-32.

[7]HeP,ZhangJ,YinX,etal.EnergeticsaltsbasedontetrazoleN-oxide[J].Chemistry, 2016, 22(23):7670-7685.

[8] 劉曉建, 林秋漢, 龐思平, 等. 唑類含能離子化合物的合成研究進展 [J]. 火炸藥學報, 2010, 33(1):6-10.LIUXiao-jian,LINQiu-han,PANGSi-ping,etal.Progressofstudyonthesynthesisofazoleenergeticioniccompounds[J].ChineseJournalExplosives&Propellants(HuozhayaoXuebao), 2010, 33(1):6-10.

[9] 黃曉川, 郭濤, 葛忠學, 等. 聯唑類含能化合物及其含能離子鹽研究進展 [J]. 含能材料. 2015, 23(3):291-301.HUANGXiao-chuan,GUOTao,GEZhong-xue.Reviewonbis-azolesanditsenergeticionderivatives[J].ChineseJournalofEnergeticMaterials, 2015, 23(3):291-301.

[10] 翟連杰, 樊學忠, 畢福強, 等. 4,5-二(1H-5-四唑基)-1,2,3-三唑的合成方法:CN,105037330[P]. 2015.ZHAILian-jie,FANXue-zhong,BIFu-qiang,etal.Thesynthesisof4,5-bis(1H-tetrazol-5-yl)-2H-1,2,3-trizole:CN, 105037330[P]. 2015.

[11]DippoldAA,IzsákD,Klap?tkeTM,etal.Combiningtheadvantagesoftetrazolesand1,2,3-triazoles: 4,5-bis(tetrazol-5-yl)-1,2,3-triazole, 4,5-bis(1-hydroxytetrazol-5-yl)-1,2,3-triazole,andtheirenergeticderivatives[J].EuropeanJournalofChemistry, 2016, 22(5):1768 -1778.

[12] 任曉雪. 國外富氮化合物技術研究進展 [J]. 飛航導彈,2016(7):87-90.RENXiao-xue.Researchprogressofnitrogen-richcompoundsabroad[J].AerodynamicMissileJournal, 2016(7):87-90.

[13]CrawfordMJ,KaraghiosoffK,Klap?tkeTM,etal.Synthesisandcharacterizationof4,5-dicyano-2H-1,2,3-tri-azoleanditssodium,ammonium,andguanidiniumsalts[J].InorganicChemistry, 2009, 48(4):1731-1743.

[14]Klap?tkeTM,MayrN,StierstorferJ,etal.Maximumcompactionofionicorganicexplosives:bis(hydroxylammonium)5.5′-dinitromethyl-3,3′-bis(1,2,4-oxadiazolate)anditsderivatives[J].Chemistry-AEuropeanJournal, 2014, 20(5):1410-1417.

[15]SuceskaM.EvaluationofdetonationenergyfromEXPLO5computer[J].Propellants,Explosives,Pyrotechnics, 1999, 24(5):280-285.

[16]SuceskaM.CalculationofdetonationparametersbyEXPLO5computerprogram[J].MaterialsScienceForum, 2004, 465/466:325-330.

[17]TangY,YangH,ShenJ,etal. 4-(1-Amino-5-aminotetrazolyl)-methyleneimino-3-methylfuroxananditsderivatives:synthesis,characterization,andenergeticproperties[J].EuropeanJournalofInorganicChemistry, 2014, 2014(7):1231-1238.

Synthesis and Characterization of 4,5-Bis(tetrazol-5-yl)-1,2,3-triazole and Its Energetic Salts

ZHANG Yan-fang, XIONG Hua-lin, LIN Qiu-han, CHENG Guang-bin, YANG Hong-wei

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

4,5-Bis(tetrazol-5-yl)-1,2,3-triazole was synthesized by diazotization coupling reaction and [2+3] azo cycloaddition reaction and seven types of energetic nitrogen-rich salts were synthesized via. metathesis reactions using 4,5-dicyano-1,2,3-triazole as raw material. The structures of compounds were characterized by IR,1H NMR and13C NMR.Their thermal behaviors were studied by DSC. Theoretical densities and standard enthalpies of formation of each compound were calculated based on B3LYP/6-311G**method, and the detonation velocities and detonation pressures of compounds were estimated by EXPLO5. The results show that eight types of compounds have better thermal stability. Their enthalpies of formation are between 647.4 kJ/mol and 1579.4kJ/mol.Their detonation velocities are between 7652m/s and 8389m/s.The detonation pressures are higher than 20 GPa except for compound 2. The amounts of gas production of other compounds are greater than 800L/kg under the standard condition, except for compounds 1,7 and 8, which are potential gas generators.

organic chemistry; high nitrogen compound; tetrazole; triazole;energetic ion salt

10.14077/j.issn.1007-7812.2017.02.007

2016-10-10;

2016-11-04

國家自然科學基金(No.21376121;No.21676147)

張艷芳(1992-),女,碩士研究生,從事含能材料的合成研究。E-mail: zhangyanfang000@163.com

楊紅偉(1979-),女,副教授,從事含能材料的合成與應用研究。E-mail: hyang@mail.njust.edu.cn

TJ55;O621

A

1007-7812(2017)02-0041-06

主站蜘蛛池模板: 亚洲欧洲日产国码无码av喷潮| 72种姿势欧美久久久大黄蕉| 欧美h在线观看| 亚洲成在人线av品善网好看| 国产色偷丝袜婷婷无码麻豆制服| 欧洲av毛片| 欧美一区二区啪啪| 国产精品不卡永久免费| 亚洲高清无码久久久| 992Tv视频国产精品| 在线观看av永久| 国产极品粉嫩小泬免费看| 国产精品欧美激情| 久热中文字幕在线| 亚洲精品动漫| 免费网站成人亚洲| 伊在人亚洲香蕉精品播放| 91成人在线免费观看| 亚洲V日韩V无码一区二区| 看你懂的巨臀中文字幕一区二区 | 国产本道久久一区二区三区| 97在线观看视频免费| 麻豆精品视频在线原创| 国产高潮流白浆视频| 国产流白浆视频| 国产美女无遮挡免费视频网站| av在线手机播放| 亚洲欧美日韩另类在线一| 亚洲香蕉在线| 亚洲精品无码高潮喷水A| 香蕉伊思人视频| 国产主播一区二区三区| 久久青青草原亚洲av无码| 国产亚洲精品资源在线26u| 中美日韩在线网免费毛片视频| 亚洲区一区| 亚洲中文无码av永久伊人| 不卡网亚洲无码| 亚洲资源站av无码网址| 日韩毛片免费| 国产九九精品视频| 婷婷成人综合| 中文无码精品a∨在线观看| 日本高清有码人妻| 亚洲性日韩精品一区二区| 直接黄91麻豆网站| 一区二区午夜| 国产69囗曝护士吞精在线视频| 亚洲二区视频| 在线看片免费人成视久网下载| 久久无码免费束人妻| 国产成人精品综合| 精品久久国产综合精麻豆| 国产精品免费p区| 亚洲欧美精品日韩欧美| 国产另类视频| 精品欧美一区二区三区久久久| 久草视频精品| 日韩精品资源| 视频一本大道香蕉久在线播放| 亚洲资源在线视频| 亚洲第一成网站| 成年看免费观看视频拍拍| 97狠狠操| 欧美国产精品不卡在线观看 | 国产欧美综合在线观看第七页| 国产黄色爱视频| 欧美视频在线播放观看免费福利资源 | 国产95在线 | 亚洲成在线观看| 国内丰满少妇猛烈精品播| 青青草91视频| 国产高清精品在线91| 中文无码精品A∨在线观看不卡 | 91福利片| 污污网站在线观看| 天天躁夜夜躁狠狠躁躁88| 欧美人与动牲交a欧美精品| 国产成人AV大片大片在线播放 | 国产精鲁鲁网在线视频| 亚洲天堂网站在线| 午夜高清国产拍精品|