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

IEC安全火花試驗裝置點燃因素研究

2023-04-12 00:00:00劉樹林于昌隆
西安科技大學學報 2023年1期

摘 要:為研究IEC火花試驗裝置的影響因素,建立IEC火花試驗裝置放電間隙電壓、電流數學模型,得到有效點燃能量表達式。在此基礎上,利用Fluent軟件建立IEC火花試驗裝置中鎢絲與鎘盤放電間隙模型,采用UDF編程模塊模擬鎢絲向鎘盤運動過程,并設置計算出的放電間隙有效點燃能量及其相關參數,仿真模擬有效點燃能量、點火半徑、電極材質3個因素對電火花引燃甲烷/空氣混合氣體的影響。結果表明:點火半徑不變的情況下,鎢電極在運動至電火花之前,溫度有短暫上升并保持不變,當與火花接觸后會瞬間吸收大量熱量,通過仿真計算電火花能量在鎢電極上的損耗占有很大比重;得到不同電容值對應的最小點火半徑,最小點火半徑隨著電容值的增大逐漸減小,并擬合出最小點火半徑與有效點燃能量之間的函數關系式。

關鍵詞:IEC安全火花試驗裝置;電容電路;短路放電;爆炸特性;數值模擬;點火

Study on igniting factors of IEC safety spark test apparatus

LIU Shulin,YU Changlong

(College of Electrical and Control Engineering,Xi’an University of Science and Technology,Xi’an 710054,China)

Abstract:In order to study the influence factors of the initiation of methane/air mixture in IEC safety spark test apparatus,the mathematical model of spark discharge gap voltage and current was established and the effective ignition energy expression was obtained.On this basis,the Fluent software was used to establish the discharge gap model between tungsten wire and cadmium disk in IEC safety spark test apparatus.Then UDF programming module was used to simulate the motion of tungsten wire to cadmium disk,and the calculated effective ignition energy and the related parameters were achieved.Finally,the effects of effective ignition energy,ignition radius and electrode material on the ignition of methane/air mixture by electric spark were simulated.The results indicate:when the ignition radius is constant,the temperature of tungsten electrode rises briefly before moving to the spark and remains unchanged;when the tungsten electrode is in contact with spark,it will instantly absorb a large amount of heat.Through simulation calculation,the loss of spark energy on tungsten electrode occupies a large proportion.The minimum ignition radius corresponding to different capacitance values is obtained.The minimum ignition radius decreases with the increase of capacitance values,and the function relation between the minimum ignition radius and effective ignition energy is fitted.

Key words:IEC safety spark test apparatus;capacitance circuit;short-circuit discharge;explosion characteristic;numerical simulation;ignition

0 引 言

隨著煤礦信息化、智能化、集成化的發展,現有本安電源越來越難以滿足煤礦井下監測、報警等設備的供電需求,為提高本安電源的輸出功率同時保證本安性能不變,需要對電路進行檢測。應用IEC(International Electrotechnical Commission)安全火花試驗裝置進行爆炸性試驗是檢驗電路本質安全性能的有效方法,但試驗結果受各種因素的影響。因此,利用IEC安全火花試驗裝置得到電容電路短路放電特性,并建立相應的數學模型,計算出電火花點燃有效能量,除了與電火花能量有關,還與點火半徑、試驗電極材質密切相關。為探究不同因素對電容短路放電點燃甲烷/空氣混合氣體的影響,迫切需要對電火花能量、點火半徑和電極材質進行研究。

胡家龍等利用Fluent軟件建立二維軸對稱計算域,設置非預混燃燒、輻射、NOx生成模型,研究甲烷在不同氧濃度、壓力下的火焰溫度、形態和生成NOx特性[1]。LIU等研究電容放電過程中電火花、導線、開關、等離子體的能量消耗,得出電容放電過程中的能量分布和電火花的動態過程及特征參數[2]。劉慶明等比較電容不同放電過程中3種放電能量計算方法,研究電容殘余能量和電路損耗能量的損耗率及其變化規律[3]。王珂利用Fluent軟件建立小火焰伴流燃燒器模型,設置24步甲烷燃燒簡化反應機理,研究甲烷在不同過量空氣系數下的燃燒溫度、燃燒速率、主要污染物生成情況[4]。孫昊等采用針-板反應器,研究甲烷氣體在納秒脈沖火花放電下的放電特性、產物特性、發射光譜特性、主要粒子光譜發展過程和轉化效率[5]。陳明仙等研究在不同點火位置情況下,甲烷泄爆產生的內部超壓、火焰演化和固有頻率等特征對甲烷泄爆容器結構響應的影響機制[6]。李暢等建立甲烷/空氣混合氣體電火花點火的數學模型,研究環境溫度、最小點火能、著火延遲時間、火焰傳播速度、點火能量對電火花點燃混合氣體的影響[7]

RAJIV等基于IEC火花試驗裝置建立電火花放電數學模型并研究電火花的熱效應[8]。李洪波等基于簡化的化學反應動力學機理,研究甲烷/空氣混合氣體中電火花溫度、半徑和點火源位置對爆炸壓力和溫度的影響,并計算出臨界點火溫度和臨界點火半徑[9]。馬秋菊等建立甲烷/空氣混合氣體電火花點火的數學模型,通過數值模擬預測電火花點火的臨界溫度及最小點火能量[10]。高志崇根據燃燒產生的發光現象,研究燃燒反應的火焰溫度,并提出烴燃燒反應的機理[11]。蔣利橋等在試驗條件下研究微尺度圓管內甲烷/空氣火焰熄滅極限特點,并擬合出淬熄速度與圓管直徑的關系式[12]。趙黛青等以數值模擬的方法研究微尺度圓管內甲烷擴散火焰結構和熄滅尺度[13]。李艷霞等以數值模擬的方法,研究微尺度下矩形截面通道的幾何尺寸對甲烷/空氣混合氣體燃燒和傳熱[14]。上述文獻主要研究不同燃燒模型下的甲烷燃燒特性、生成產物,但并未研究電容短路放電時,有效點燃能量、點火半徑、電極材質對點燃甲烷空氣混合氣體的影響。

為研究不同因素對電容電路短路電火花點燃甲烷/空氣混合氣體的影響,依據試驗電路圖計算釋放到火花間隙中的電火花能量;利用Fluent軟件數值模擬有效點燃能量、點火半徑、電極材料3個因素對電火花引爆甲烷/空氣混合氣體的影響。為提高本安電源輸出功率奠定理論基礎。

1 IEC安全火花試驗裝置模型

1.1 模型建立

簡單電容電路短路放電原理和安全火花試驗裝置如圖1所示,其中,U為電源電壓,V;R1為充電電阻,Ω;C為被測試的電容,μF;uC為電容短路放電時兩端電壓,V;iC為電容短路放電時放電電流,A;R為放電回路中的電阻,Ω;G為安全火花試驗裝置,電容兩端分別連接鎢絲和鎘盤電極;uG為安全火花試驗裝置放電電壓,V;iG為安全火花試驗裝置放電電流,A。IEC安全火花試驗裝置放電間隙結構如圖2所示[8]。電容短路放電過程可分為火花產生階段、火花維持階段、電極閉合火花熄滅3個階段[15]。電容短路放電的3個階段中,火花產生階段是點燃甲烷/空氣混合氣體的主要能量來源[16-17],因此將主要研究火花產生階段的能量變化。

如圖2所示,鎢絲下伸到鎘盤2 mm處,當鎢絲運動到與鎘盤間距小于1 mm時,電容短路放電產生電火花引爆腔室內甲烷/空氣混合氣體,為簡化模型的復雜程度,采用二維坐標系建立鎢絲-鎘盤放電間隙模型[18],如圖3所示。

通過圖1可得電容短路放電回路方程

uC=uR+uG "(1)

式中 uG為安全火花試驗裝置兩電極間電壓;uC為電容兩端電壓,uR為回路電阻R兩端電壓。

式中 iC為電容短路放電時放電電流,A;iG為安全火花試驗裝置放電電流,A。

根據對大量電容性電路閉合火花放電波形的分析,發現火花放電電壓在擊穿后呈指數規律下降[19-24],因此,可以作如下假設

式中 u0為放電維持電壓值,V;a0,a1,a2,a3為未知系數;t為電容短路放電持續時間,μs。

聯立式(2)(3)(4)求解,可得

則式(3)(4)可寫為

電容在火花產生階段釋放出的能量即為點燃甲烷/空氣混合氣體的能量,但電容短路放電釋放出的能量除了用于點燃甲烷/空氣混合氣體,還有一部分消耗于電路和電極上[2,3,25]。因而電容短路放電用于點燃的有效能量可用下式表示

Ereal=Etotal-Eloss "(8)

式中 Ereal為電容短路放電用于點燃的有效能量,μJ;Etotal為電容短路放電總能量,μJ;Eloss為電路中消耗的能量,μJ。

其中,Etotal,Eloss表達式分別為

式中 t1為電容在火花產生階段放電時間[17],μs。

由式(8)(9)(10)(11)可得電容有效點燃能量。

1.2 模型簡化

為簡化甲烷燃燒模型作出如下假設。

1)忽略間隙中電場和磁場對模型的影響。

2)間隙中的混合氣體是理想氣體。

3)反應過程中導熱系數、定壓比熱容、活化能及指前因子等參數不隨反應的進行而改變。

4)忽略等離子體對燃燒的影響。

5)電火花的放電功率恒定。

6)在能量守恒方程中忽略黏性耗散。

1.3 控制方程

狀態方程

式中 ρ為密度,kg/m3;R′為普適氣體常數,J/(mol·K);T為溫度,K;Yi為i組分的質量分數;mi為i組分的相對分子質量,kg/mol。

質量守恒方程

式中 t為時間,s;r為計算節點到火花中心的距離,m;u為速度,m/s。

動量守恒方程

式中 Ρ為壓力,Pa;μ為黏度,Pa/s。

能量守恒方程

式中 CP為定壓比熱容,J/(kg·K);λ為熱傳導系數,J/(m·K·s);k0為指前因子,s-1;Ea為表觀活化能,J/mol;Hv為燃燒反應放熱,J/kg;rs為電火花半徑,m;

Ereal為電火花能量,J;t1為電火花放電時間,s;Yf為燃料質量分數。

組分守恒方程

式中 Di為有效擴散系數,m2/s;ni為i組分燃燒1 kg燃料消耗/生成對應組分的質量,kg。

1.4 初始條件和邊界條件

Fluent軟件內置電火花設置模塊,可設置電火花在計算區域內的位置、火花半徑、持續時間、放電能量。假設電火花產生的初始恒定高溫區域由火花半徑確定,而大于火花半徑的區域認為氣體溫度為室溫。

式中 T1為點火半徑內的初始溫度,K;Ta為初始時刻氣體的溫度,即室溫,K;C1,C2為邊界兩側氣體的定容比熱容,J/(kg·K)。

在圖3所示的二維計算域中,3個邊界認為是外邊界,其中鎢絲、鎘盤為固體邊界、其余部分為氣體邊界,邊界初始溫度設置為300 K,甲烷濃度為8.3%,初始溫度為300 K[26],壓力為常壓。整個計算域的網格尺寸設置為0.05 mm,由于電容短路放電時間是微秒量級,設置時間步長為1×10-6 s,總的時間步由式(11)計算得到。在Fluent軟件中,求解模式為組分輸運方式,為了較準確模擬放電過程中鎢絲運動至鎘盤的動態過程,采用UDF編程定義鎢絲運動軌跡、動網格實現模型的動態變化。采用PISO算法只需求解一次動量方程,即可對壓力方程及速度修正方程進行迭代循環從而達到收斂,減少迭代次數。

2 結果和討論

為研究最危險的情況,放電回路電阻取值為零歐姆,則電路總電阻等于電路內阻和火花試驗裝置內阻之和,IEC安全火花裝置試驗電路參數:選取電容值C=47 μF,依據Ⅰ類電容短路最小點燃電壓曲線,電容初始充電電壓U=19 V,根據多次試驗結果電路內阻取平均值為R=0.3 Ω,維持電壓取平均值為u0=10 V。代入式(11),可得一次放電時間為t1=50.29 μs,同理,根據式(8)、式(9)、式(10)可得Ereal=3 449 μJ,試驗結果與仿真計算結果基本一致,如圖4所示。

用Fluent軟件求解電火花點燃甲烷空氣混合氣體的過程,得到計算域內溫度隨時間的變化情況,混合氣體發生爆炸時溫度會急劇上升,因此用溫度變化判斷混合氣體是否發生爆炸。如圖5的點火云圖和圖6所示的溫度-時間曲線,在第16 μs時,溫度突然急劇上升,則認為發生爆炸;如圖7的點火云圖和圖8所示的溫度-時間曲線,整個放電過程中溫度雖然一直在升高,但沒有產生溫度的突然躍升,則認為沒有發生爆炸。

為探究鎢絲電極的運動對電火花點燃甲烷/空氣混合氣體的影響程度,監測鎢絲電極運動過程中溫度的變化,由此推斷鎢絲電極在整個放電過程中吸收的能量。如圖9和圖10所示,根據圖11的監測結果可計算出鎢絲電極消耗的能量占放電總能量的43.9%。

為探究電火花半徑對點燃甲烷/空氣混合氣體的影響,對不同電容值設置不同半徑的電火花,并用Fluent軟件進行仿真,可得到最小點火半徑,當火花半徑大于最小點火半徑時,混合氣體被引爆,當火花半徑小于最小點火半徑時,混合氣體無法被引爆。由此得到不同電容值對應的最小點火半徑,如圖12所示。隨著電容值的增大最小點火半徑逐漸變小,這是因為電容值增大后,其短路放電火花釋放出的能量也相應增大,點火過程中能量損失的影響相對越來越小,導致最小點火半徑不斷變小,如圖13所示。利用Matlab軟件可得擬合表達式。

r=-8.315×10-8×E1.654real+0.179 2 "(20)

式中 r為最小點火半徑,mm;Ereal為有效點燃能量,μJ。

3 結 論

1)建立電容短路時放電間隙的電壓、電流數學模型,得到有效點燃能量的表達式。

2)點火半徑不變的情況下,鎢電極運動至電火花之前,溫度有短暫上升之后保持不變,在與火花接觸后會瞬間吸收大量的能量,又通過仿真計算可得出鎢電極上損耗的能量,在電火花能量中占有很大比重。

3)得到不同電容值對應的最小點火半徑,且點火半徑隨著電容值的增大逐漸減小,并擬合出最小點火半徑與有效點燃能量之間的函數關系式。

參考文獻(References):

[1] 胡家龍,任航,聶曉康,等.氧濃度、壓力對甲烷/空氣層流擴散燃燒特性影響的數值模擬研究[J].東北電力大學學報,2021,41(3):10-18.HU Jialong,REN Hang,NIE Xiaokang,et al.Numerical simulation of effects of oxygen concentration and pressure on methane/air laminar diffusion combustion characteristics[J].Journal of Northeast Electric Power University,2021,41(3):10-18.

[2]LIU Q M,HUANG J X,SHAO H G,et al.Electric ignition energy evaluation and the energy distribution structure of energy released in electrostatic discharge process[J].Chinese Physics B,2017,26(10):1-6.

[3]劉慶明,汪建平,李磊,等.電火花放電能量及其損耗的計算[J].高電壓技術,2014,40(4):1255-1260.LIU Qingming,WANG Jianping,LI Lei,et al.Calculation of electric spark discharge energy and its energy loss[J].High Voltage Engineering,2014,40(4):1255-1260.

[4]王珂.過量空氣系數對甲烷燃燒特性的影響[J].工業加熱,2020,49(11):40-43.WANG Ke.Effect of excess air coefficient on methane combustion characteristics[J].Industrial Heating,2020,49(11):40-43.

[5]孫昊,張帥,韓偉,等.納秒脈沖火花放電高效轉化甲烷的實驗研究[J].電工技術學報,2019,34(4):880-888.SUN Hao,ZHANG Shuai,HAN Wei,et al.An experimental investigation of nanosecond pulsed spark discharge for high-efficient methane conversion[J].Transactions of China Electrotechnical Society,2019,34(4):880-888.

[6]陳明仙,郭進,羅飛云,等.點火位置對甲烷/空氣預混泄爆容器結構響應影響研究[J].西安科技大學學報,2021,41(5):800-807.CHEN Mingxian,GUO Jin,LUO Feiyun,et al.Influence of ignition position on structural response of methane/air premixed explosion venting vessel[J].Journal of Xi’an University of Science and Technology,2021,41(5):800-807.

[7]李暢,馬澤鵬,苑藝笑,等.甲烷/空氣混合物的火花著火敏感性分析[J].安全與環境學報,2022,22(4):1913-1918.LI Chang,MA Zepeng,YUAN Yixiao,et al.Theoretical study on electric spark ignition sensitivity of methane/air mixture[J].Journal of Safety and Environment,2022,22(4):1913-1918.

[8]RAJIV S,CARSTEN U.Modelling of sparking contacts for hazardous area applications[C]//2015 IEEE 61st Holm" Conference on Electrical Contacts(Holm),2015:347-352.

[9]李洪波,崔淦,李自力,等.基于簡化機理的甲烷/空氣混合氣體爆炸特性研究[J].消防科學與技術,2017,36(9):1197-1202.LI Hongbo,CUI Gan,LI Zili,et al.Study on the explosion characteristics of methane/air mixture using simplified chemical kinetics[J].Fire Science and Technology,2017,36(9):1197-1202.

[10]馬秋菊,張奇,龐磊.甲烷/空氣最小點火能量預測理論模型[J].高壓物理學報,2012,26(3):301-305.MA Qiuju,ZHANG Qi,PANG Lei.Theoretical model of minimum ignition energy prediction for methane/air mixture[J].Chinese Journal of High Pressure Physics,2012,26(3):301-305.

[11]高志崇.甲烷燃燒反應的火焰溫度[J].山西大學學報,2004,27(1):32-34.GAO Zhichong.Flame temperature from CH4 combustion reaction[J].Journal of Shanxi University,2004,27(1):32-34.

[12]蔣利橋,趙黛青,汪小憨.微尺度甲烷擴散火焰及其熄滅特性[J].燃燒科學與技術,2007(2):183-186.JIANG Liqiao,ZHAO Daiqing,WANG Xiaohan.Structure and extinction characteristics of methane

micro-diff-usion flames[J].Journal of Combustion Science and Technology,2007(2):183-186.

[13]趙黛青,黃顯峰,山下博史,等.微尺度甲烷擴散火焰特性的數值解析[J].燃燒科學與技術,2005(6):530-534.ZHAO Daiqing,HUANG Xianfeng,YAMASHITA Hiroshi,et al.Numerical analysis of methane micro diffusion flames[J].Journal of Combustion Science and Techno-logy,2005(6):530-534.

[14]李艷霞,劉中良,桑麗霞.微小型通道幾何參數對甲烷燃燒的影響[J].工程熱物理學報,2012,33(12):2155-2158.LI Yanxia,LIU Zhongliang,SANG Lixia.Effects of geometry of micro-channel on the methane/air combustion[J].Journal of Engineering Thermophysics,2012,33(12):2155-2158.

[15]劉樹林,劉健.本質安全開關變換器[M].北京:科學出版社,2008.

[16]鐘久明,劉樹林,崔強.IEC火花試驗裝置的電容短路放電特性數學仿真分析[J].電工電能新技術,2014,33(2):29-34.ZHONG Jiuming,LIU Shulin,CUI Qiang.Short circuit discharge behavior of capacitive circuit and its mathematical simulation analysis[J].Advanced Technology of Electrical Engineering and Energy,2014,33(2):29-34.

[17]劉樹林,鐘久明,樊文斌,等.電容電路短路火花放電特性及其建模研究[J].煤炭學報,2012,37(12):2123-2128.LIU Shulin,ZHONG Jiuming,FAN Wenbin,et al.Short circuit discharge characteristics of the capacitive circuit and its mathematical model[J].Journal of China Coal Society,2012,37(12):2123-2128.

[18]WEI R N,LI C Y,ZHOU R F,et al.Experimental study on the effects of the ignition parameters on the spark characteristics and the flame propagation of premixed methane air mixtures[J].

Proceedings of the" Institution of Mechanical Engineers,Part D:Journal of Automobile Engineering,2016,230(7):914-927.

[19]劉輝.隔爆兼本安直流穩壓開關電源的研究[D].西安:西安科技大學,2005.LIU Hui.Research on flame-proof amp; intrinsically safe DC switching power supply[D].Xi’an:Xi’an University of Science and Technology,2005.

[20]JITHIN S,RAUT A,JOSHI S S,et al.Finite element model for topography prediction of electrical discharge textured surfaces considering multi-discharge phenomenon[J].International Journal of Mechanical Sciences,2020,177(6):1-16.

[21]JITHIN S,BHANDARKAR U V,JOSHI S S.Multi-spark model for predicting surface roughness of electrical discharge textured surfaces[J].The International Journal of Advanced Manufacturing Technology,2020,106(9):3741-3758.

[22]WANG G Q,MEI B W,LIU X C,et al.Investigation on spherically expanding flame temperature of nbutane/air mixtures with tunable diode laser absorption spectroscopy[J].Proceedings of the Combustion Institute,2019,37(2):1589-1596.

[23]ZHEN X D,LI X Y,WANG Y,et al.Effects of the initial flame kernel radius and EGR rate on the performance,combustion and emission of high-compression spark ignition methanol engine[J].Fuel,2020,262:1-12.

[24]BREQUIGNY P,ENDOUARD C,FOUCHER F.Impr-ovement of turbulent burning velocity measurements by schlieren technique,for high pressure isooctane-air premixed flames[J].Combustion Science and Technology,2020,192(3):416-432.

[25]張云明,劉慶明,宇燦,等.大能量電點火系統設計與火花放電特性實驗研究[J].高電壓技術,2014,40(4):1267-1274.ZHANG Yunming,LIU Qingming,YU Can,et al.Design of a big energy electric ignition system and experimental study on its spark discharge characteristics[J].High Voltage Engineering,2014,40(4):1267-1274.

[26]陳國華,董浩宇,張強,等.狹長受限空間甲烷/空氣爆炸事故研究評述[J].安全與環境學報,2020,20(3):946-959.CHEN Guohua,DONG Haoyu,ZHANG Qiang,et al.Review on the methane/air explosion accidents in the narrow confined space[J].Journal of Safety and Environment,2020,20(3):946-959.

(責任編輯:高 佳)

主站蜘蛛池模板: 高清亚洲欧美在线看| 国产欧美又粗又猛又爽老| 欧美高清国产| 国产精品成人第一区| 五月天综合婷婷| 九九热视频在线免费观看| 久久永久精品免费视频| 国产成+人+综合+亚洲欧美 | 日本AⅤ精品一区二区三区日| 综合亚洲网| 久久亚洲欧美综合| 美女被狂躁www在线观看| 国产精品久久自在自线观看| 亚洲欧美极品| 精品亚洲麻豆1区2区3区 | 亚洲天堂免费观看| 国产亚洲美日韩AV中文字幕无码成人| 国产高清无码第一十页在线观看| 欧美视频免费一区二区三区| 亚洲swag精品自拍一区| 中文字幕亚洲精品2页| 最新精品国偷自产在线| 久久无码av三级| 亚洲开心婷婷中文字幕| 凹凸国产熟女精品视频| 日韩欧美成人高清在线观看| 精品偷拍一区二区| 国产精品v欧美| 国产麻豆精品久久一二三| 国产在线观看高清不卡| 久久久久久国产精品mv| 一级毛片无毒不卡直接观看| 三区在线视频| 成人午夜天| 国产亚洲精品自在线| 人妻丰满熟妇啪啪| 91丨九色丨首页在线播放| 亚洲妓女综合网995久久| 乱码国产乱码精品精在线播放| 日韩国产精品无码一区二区三区| 国产va欧美va在线观看| 亚洲h视频在线| 人禽伦免费交视频网页播放| 香蕉久久国产超碰青草| 亚洲永久色| 国产区成人精品视频| 婷婷丁香在线观看| 午夜国产在线观看| 99久视频| 欧美乱妇高清无乱码免费| 视频二区中文无码| 黄色网站在线观看无码| 国产一国产一有一级毛片视频| 狠狠亚洲五月天| 尤物视频一区| 青青草原国产| 97se综合| 国产成人亚洲日韩欧美电影| 天堂成人在线| 亚洲精品视频网| 中文字幕在线看视频一区二区三区| 88av在线播放| 欧美a√在线| 欧美精品高清| 日韩高清在线观看不卡一区二区| 日本欧美精品| 毛片卡一卡二| 亚洲综合激情另类专区| 亚洲欧洲日产国码无码av喷潮| 最新午夜男女福利片视频| 天天视频在线91频| 91免费国产高清观看| 久久精品国产一区二区小说| 欧美色图久久| 亚洲欧美另类色图| 少妇极品熟妇人妻专区视频| 久久综合九色综合97婷婷| 2021国产v亚洲v天堂无码| 青青青国产视频| 国产福利免费视频| 国内精品伊人久久久久7777人| 亚洲欧美另类视频|