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

    氣體泄爆壓力分步計算模型及其湍流修正*

    2019-06-21 02:55:06高康華邱艷宇王明洋
    爆炸與沖擊 2019年5期
    關(guān)鍵詞:腔體湍流燃氣

    孫 松,高康華,邱艷宇,王明洋

    (陸軍工程大學爆炸沖擊防災減災國家重點實驗室,江蘇 南京 210007)

    可燃性氣體因其來源廣、熱能高、污染小等特點已經(jīng)成為人們生產(chǎn)、生活中最常用的能源之一;與此同時,頻繁發(fā)生的氣體爆炸事故也給人們的生命財產(chǎn)安全帶來了巨大威脅。泄爆是工程實際中最常用的防護措施,研究者們對于泄爆條件下空間內(nèi)的爆炸壓力進行了大量的實驗研究[1-5],對氣體爆炸的物理過程及各影響因素有了較全面的了解。

    同時,出于安全性、經(jīng)濟性和實用性的考慮,研究者們以大量實驗為基礎(chǔ)結(jié)合氣體爆炸物理過程進行了理論模型的研究,以期能夠預測空間內(nèi)氣體爆炸壓力的變化情況。Bradley 等[6]根據(jù)守恒方程、理想氣體狀態(tài)方程和等熵壓縮方程推導出關(guān)于升壓速率的偏微分方程組,成為其后理論模型發(fā)展的基礎(chǔ)。Han 等[7]通過對比模型計算值與實驗結(jié)果發(fā)現(xiàn),壓力計算模型可以較好地描述建筑空間內(nèi)的壓力變化,隨后通過模型計算探討了泄放面積、泄放壓力、屋室尺寸、氣體濃度和點火位置對爆燃壓力的影響。Ugarte 等[8]根據(jù)不同點火位置、容器形狀在一維等熵泄放模型的基礎(chǔ)上對火焰面積、火焰?zhèn)鞑ニ俣群屯獠勘級毫M行修正,對比發(fā)現(xiàn)修正后的模型與實驗結(jié)果吻合更好。美國防火規(guī)范NFPA68-2013[9]中總結(jié)大量實驗數(shù)據(jù)得到不同工況下泄爆壓力及泄爆面積的簡化確定方法,為實際工程防護提供了重要指導。Sustek 等[10]使用常用經(jīng)驗與半經(jīng)驗公式對不同實驗工況進行計算,對比發(fā)現(xiàn),由于各經(jīng)驗公式中的參數(shù)很多都是基于實驗數(shù)據(jù)擬合得到,因此各公式適用范圍較窄,使用經(jīng)驗公式計算不同工況可能會產(chǎn)生較大誤差。

    現(xiàn)有的氣體爆燃壓力計算方法多是以大量實驗數(shù)據(jù)為基礎(chǔ)評估空間內(nèi)的最大峰值超壓,無法準確描述空間內(nèi)的壓力變化過程;數(shù)值模擬方法計算成本高,耗時長,不利于工程快速評估應用;而常用的壓力計算模型[11-13]多為關(guān)于升壓速率的偏微分方程組,其求解過程較復雜,同時這些模型大多假設泄放構(gòu)件開啟后先泄放已燃氣體,再泄放未反應氣體,無法計算同時泄放已燃氣體和未反應氣體的工況。本文中,通過將爆燃過程進行微分,假設每個微小時間段內(nèi)爆燃泄放過程均按照先燃燒、再泄放、最后壓力平衡的過程獨立分步進行,得到了泄爆壓力分步計算模型;同時在尺寸為2 m×1.2 m×0.6 m 的爆炸腔體一端安裝不同面積的泄壓結(jié)構(gòu)進行泄爆實驗,對分步壓力計算模型進行驗證,并根據(jù)實驗結(jié)果對模型進行湍流修正。該模型清晰地反映了氣體爆燃泄放機理,通過對微分時段內(nèi)的分步計算過程進行迭代,可以較好地描述泄爆條件下空間內(nèi)的壓力及各狀態(tài)參數(shù)的變化過程,計算過程簡單,可為工程實際中快速評估內(nèi)爆炸壓力提供參考;經(jīng)過湍流加速修正的壓力計算模型可以較好地描述因泄放截面改變引起的湍流區(qū)域內(nèi)的壓力變化情況。

    1 理論模型

    1.1 基本假設

    由于爆燃過程非常復雜,涉及到燃燒學、熱動力學、流體動力學等多方面內(nèi)容,因此在模型計算時假設可燃氣體密度均勻充滿密閉空間且符合理想氣體狀態(tài)方程;同時假設密閉容器為絕熱壁面,忽略熱交換及內(nèi)部反應的能量耗散過程。

    泄放條件下氣體爆燃的超壓變化過程受燃燒產(chǎn)生的壓力和泄放降低的壓力共同影響,變化過程較復雜,為簡化計算將燃燒過程和泄放過程單獨考慮,將氣體爆燃泄放過程分成3 個獨立階段:(1)只燃燒不泄放(confined combustion stage)。將此階段簡化為密閉空間內(nèi)的爆燃過程,不考慮泄放對于爆燃壓力的影響。(2)泄放過程(venting stage)。第1 階段燃燒結(jié)束后室內(nèi)壓力上升,高于外部環(huán)境壓力,此時燃燒產(chǎn)物和未燃氣體將分別發(fā)生泄放。(3)泄放后壓力平衡階段(pressure equilibrium stage)。該階段泄放過程結(jié)束,此時空間內(nèi)部燃燒產(chǎn)物與未燃氣體壓力不同,空間內(nèi)部會發(fā)生氣體流動最終達到壓力平衡。密閉空間內(nèi)氣體爆燃泄放過程如圖1 所示,圖1(a)~(c)為燃燒反應階段,圖1(d)~(e)為泄放階段,圖1(f)為壓力平衡階段。

    圖 1 爆燃泄放分步計算模型示意圖Fig.1 Diagrams of the sub-steps model for combustion venting

    1.2 控制方程

    1.2.1 燃燒階段

    通常假設預混氣體以規(guī)則球形燃燒并向四周擴散,則在 Δt時間內(nèi)所燃燒的體積為:

    式中: VB,j為第j 個時間段內(nèi)燃燒的氣體體積;S 為火焰?zhèn)鞑ニ俣龋?Vh,j-1為j-1 時刻已燃氣體體積。

    現(xiàn)定義燃盡率 Δμj為 Δt時間內(nèi)燃燒的氣體質(zhì)量與可燃氣體總質(zhì)量之比,則:

    式中: ρc,j-1為 j-1 時刻未燃氣體密度, V0為 爆炸空間總體積, ρ0為初始氣體密度。

    又因為根據(jù)等熵壓縮方程,在絕熱過程中有:

    式中: γu為未燃氣體絕熱指數(shù), pj-1為j-1 時刻空間內(nèi)的壓力,p0為空間內(nèi)的初始壓力。

    則 Δμj可以寫成:

    由于空間內(nèi)的壓力隨氣體燃燒分數(shù)的增大而升高,假設室內(nèi)壓力的升高與燃燒產(chǎn)物生成量成線性關(guān)系[14],則圖1 中(c)階段室內(nèi)壓力為:

    式中: pb為 密閉空間內(nèi)氣體爆燃產(chǎn)生的最大壓力,該值與氣體種類和濃度有關(guān); pn,j為j 時刻密閉燃燒階段結(jié)束時內(nèi)部壓力值。

    由于整個反應過程為絕熱過程且壓力 pn,j已知,則根據(jù)等熵壓縮方程:

    可得圖1 中j 時刻(c)階段未燃氣體體積為:

    已燃氣體體積為:

    式中: Vc,j-1為 j-1 時刻未燃氣體體積; Vh,j-1為j-1 時刻已燃氣體的體積。

    由質(zhì)量守恒可得圖1 中(c)階段未燃氣體密度為:

    式中: ρn,c,j為j 時刻燃燒階段結(jié)束時未燃氣體密度, ρc,j-1為j-1 時刻未燃氣體密度。

    已燃氣體密度為:

    式中: ρn,h,j為j 時刻燃燒階段結(jié)束時燃燒產(chǎn)物密度, ρh,j-1為j-1 時刻燃燒產(chǎn)物密度。

    此時第j 個時間段內(nèi)燃燒階段結(jié)束,開始進行泄放階段的壓力計算。

    1.2.2 泄放階段

    考慮實際情況,假設該階段未燃氣體與燃燒產(chǎn)物分別發(fā)生泄放且外界氣壓恒為標準大氣壓,則以次臨界速度( p0/ pj-1≥βk)泄放時[6,14],單位時間內(nèi)泄放的燃燒產(chǎn)物和未燃氣體體積分別為:

    以臨界速度( p0/ pj-1≤βk)泄放時,單位時間內(nèi)泄放的燃燒產(chǎn)物和未燃氣體體積分別為:

    式中: Ah,j和 Ac,j分別為 Δt 時間內(nèi)燃燒產(chǎn)物和未燃產(chǎn)物的泄放流動面積; γb為 燃燒產(chǎn)物的絕熱指數(shù); βk為臨界比值,其計算公式為 βk=[2/(γu+1)]γu/(γu-1)。

    由于發(fā)生氣體泄放,未燃氣體和燃燒產(chǎn)物的密度及壓力均發(fā)生變化,根據(jù)等熵壓縮方程可得圖1 中(e)階段j 時刻泄放階段結(jié)束時空間內(nèi)未燃氣體和燃燒產(chǎn)物的壓力分別為:

    j 時刻泄放階段結(jié)束時空間內(nèi)未燃氣體和燃燒產(chǎn)物的密度分別為:

    在第2 階段中假設只發(fā)生氣體的泄放流動,但空間內(nèi)未燃氣體與已燃氣體所占體積不發(fā)生改變,此時第2 階段結(jié)束,開始進行空間內(nèi)壓力平衡階段的計算。

    1.2.3 壓力平衡階段

    由等熵壓縮方程可知:

    將式(19)代入式(20)并變形后得完整反應過程結(jié)束時空間內(nèi)的壓力:

    式中: β= pr,j/pj。

    求得圖1 中(f)階段的壓力pj后,根據(jù)狀態(tài)方程與等熵壓縮方程可以求得(f)階段的氣體物理參數(shù):

    式中: ρh,j和 ρc,j分別為j 時刻整個爆燃泄放過程結(jié)束后空間內(nèi)未燃氣體和燃燒產(chǎn)物的密度; Vh,j和 Vc,j分別為此時空間內(nèi)未燃氣體和燃燒產(chǎn)物的體積。

    此時第3 階段結(jié)束,一個微分時間段內(nèi)腔體中氣體壓力、密度、體積的完整變化過程結(jié)束,開始迭代進行下一個時間步的計算。

    反應后期若空間內(nèi)混合氣體已全部燃盡,則燃燒過程結(jié)束,此時氣體繼續(xù)由內(nèi)部流向外部空間,直至內(nèi)外壓力平衡。該階段在迭代計算過程中省略圖1 中(b)、(c)過程,即該情況下完整反應流程圖變?yōu)閳D1(a)、(d)、(e)、(f)。

    2 實驗驗證

    2.1 實驗裝置

    實驗裝置由配氣系統(tǒng)、長方體爆炸容器、點火系統(tǒng)、泄爆系統(tǒng)與數(shù)據(jù)采集系統(tǒng)構(gòu)成,如圖2 所示。

    圖 2 實驗系統(tǒng)示意圖Fig.2 Schematic of the experimental system

    爆炸腔體為自行設計的長方體容器,長、寬、高分別為2、1.2、0.6 m,最大可承受1 MPa 的壓力。爆炸容器一端為開敞的用于安裝泄爆結(jié)構(gòu)的卡箍式法蘭盤,其尺寸為1.2 m×0.6 m,泄爆構(gòu)件示意圖如圖3所示;另一端壁面中部通過法蘭結(jié)構(gòu)與電阻絲點火器連接,構(gòu)成點火系統(tǒng),其最大點火溫度可達500 ℃。實驗中選用乙烯作為可燃氣體,通過在爆炸腔體側(cè)面等距離設置進氣孔和抽氣孔運用壓力分配法配氣,氣體濃度由流量計和壓力表控制。進氣孔內(nèi)配置有細長進氣管道,管道上均勻分布有小孔,這種均布直管式進氣系統(tǒng)可以保證腔體內(nèi)預混氣體混合較均勻。在爆炸腔體頂部裝有PCB113B26 系列壓電式高頻壓力傳感器,分別距泄放端口1.2、0.1 m。數(shù)據(jù)采集系統(tǒng)將采集到的電信號轉(zhuǎn)化為壓力信號并輸出,其數(shù)據(jù)采樣頻率為200 kHz。高速攝影儀置于距泄爆口側(cè)面8 m 處,拍照頻率為1 000 s-1。

    圖 3 泄爆構(gòu)件示意圖(單位為mm)Fig.3 Schematic diagrams of vent covers (unit in mm)

    如圖3 所示,實驗中泄爆構(gòu)件采用擊穿壓力均為7.5 kPa 的泄爆板和泄爆膜。泄爆板泄放面積為1.1 m×0.5 m,泄爆膜泄放面積為0.3 m×0.6 m。采用泄爆膜作為泄爆構(gòu)件開展實驗時,在腔體泄爆端安裝尺寸為1.2 m×0.6 m 的隔板,并在隔板中部開設面積為0.6 m×0.3 m 的泄爆口,使用泄爆膜對其進行封閉,如圖3(a)所示;采用泄爆板作為泄爆構(gòu)件開展實驗時,直接在腔體泄爆端安裝尺寸為1.2 m×0.6 m 的泄爆板,如圖3(b)所示。構(gòu)件擊穿壓力通過在泄爆結(jié)構(gòu)背面粘貼應變片以應變信號消失時刻腔體內(nèi)對應的壓力確定,實驗時選取的可燃氣體的體積分數(shù)為4%~8%。

    2.2 結(jié)果對比

    使用分步壓力計算模型對實驗工況進行計算,由于實驗工況為矩形容器端部中心點火,因此需針對火焰的傳播形狀對計算模型進行修正。由高速攝影觀察到火焰由泄爆端流出時其陣面可近似認為是平面,因此為計算簡便假設火焰未接觸壁面時為球形傳播,接觸壁面時為平面?zhèn)鞑?,則端部點火時火焰的體積表達式為:

    式中:2a、2c 為爆炸腔體的長和高,x 為火焰半徑。

    泄爆板和泄爆膜工況的計算值及實驗測量結(jié)果分別如圖4~5 所示。由圖4 可知:采用泄爆板作為泄爆構(gòu)件時,2 個傳感器測得的壓力時程曲線基本重合,符合準靜態(tài)壓力形式[15];此時腔體內(nèi)部為單峰值壓力曲線,分析認為當構(gòu)件擊穿后由于泄放面積較大使得短時間內(nèi)大量氣體外泄,引起腔體內(nèi)部壓力下降至環(huán)境壓力。該情況下模型計算值與實驗結(jié)果吻合較好。

    圖 4 可燃氣體的體積分數(shù)不同、采用泄爆板作為泄爆構(gòu)件時,泄爆壓力的計算值與實驗測量結(jié)果的比較Fig.4 Comparison of gas explosion venting pressure between calculated results by the model and experimental ones in the experiments with vent plates at different volume fractions of combustible gas

    圖 5 可燃氣體的體積分數(shù)不同、采用泄爆膜作為泄爆構(gòu)件時,泄爆壓力的計算值與實驗測量結(jié)果的比較Fig.5 Comparison of gas explosion venting pressures between the model and the tests with vent films at different volume fractions of combustible gas

    觀察圖5 發(fā)現(xiàn),當采用泄爆膜作為泄爆構(gòu)件時腔體內(nèi)壓力曲線均為雙峰值曲線,分析認為由于泄爆膜面積較小,構(gòu)件擊穿后腔體內(nèi)部氣體外流引起壓力下降產(chǎn)生第1 個壓力峰值;氣體外泄同時又會加劇腔體內(nèi)部湍流擾動,扭曲火焰陣面,提高燃燒速率與升壓速率,引起腔體內(nèi)壓力再次上升產(chǎn)生第2 個峰值[1]。同時由圖5 發(fā)現(xiàn)當乙烯體積分數(shù)較高時靠近泄爆口的壓力傳感器T1 測得的第2 個峰值壓力要大于內(nèi)部傳感器T2 相應的測量值,分析原因認為由于泄爆膜開口面積較小,外流氣體流動面積的突然改變、結(jié)合容器壁的限制作用使得近泄爆口處湍流效應明顯,可燃氣體得到充分燃燒,腔體內(nèi)部產(chǎn)生壓力梯度,泄爆口處壓力大于腔體內(nèi)部壓力,且越接近最佳體積分數(shù)該現(xiàn)象越明顯[16-17]。此時理論模型可以較好地描述腔體內(nèi)部的壓力變化情況,但與近泄爆口處壓力變化情況存在誤差。

    為了進一步驗證計算模型的適用性,運用該模型計算其他研究者的實驗工況[5,7,18],其結(jié)果如圖6所示。通過對比發(fā)現(xiàn)計算結(jié)果與實驗結(jié)果吻合較好,該方法具有較好的普適性。

    圖 6 運用本文計算模型計算其他研究者的實驗工況[5, 7, 18]下的瀉爆壓力Fig.6 Gas explosion venting pressures calculated by the model proposed in this paper for the experimental conditions[5, 7, 18] of other researchers

    3 模型修正

    由于傳感器T1 測得的第2 個峰值壓力大于腔體內(nèi)部相應的壓力值是因泄放截面改變引起的火焰湍流傳播引起,因此在此引入湍流加速因數(shù)λ 對火焰?zhèn)鞑ニ俣冗M行修正,以計算靠近泄放口位置的壓力情況。湍流加速作用因泄放截面的改變引起,因此將傳播路徑中的障礙物截面與火焰?zhèn)鞑ソ孛娴谋戎担醋枞? 作為其增強因子的影響因素。研究者們對障礙物對火焰?zhèn)鞑サ募铀僮饔眠M行了大量的實驗研究[16,19-20],雖然實驗環(huán)境不同,但是障礙物阻塞率對于火焰的湍流加速效果表現(xiàn)出一定的規(guī)律性?,F(xiàn)根據(jù)美國防火規(guī)范NFPA68 確定因截面改變引起的湍流加速因數(shù)[9]:

    式中:Aobs為障礙物截面面積,As為火焰?zhèn)鞑ヂ窂浇孛婷娣e, λ0為因自湍流等因素引起的湍流因數(shù)。在本文的泄爆膜實驗工況中Aobs取0.37, As取0.55, λ0取2。

    則未因傳播截面改變引起湍流加速的火焰?zhèn)鞑ニ俣葹閇21]:

    近泄爆口處因截面改變產(chǎn)生湍流加速的火焰?zhèn)鞑ニ俣葹椋?/p>

    式中: S0為基本火焰?zhèn)鞑ニ俣?;p 和p0分別為空間內(nèi)的壓力和初始壓力;T0和Tu分別為初始溫度和未燃氣體溫度; α 和 η為 根據(jù)實驗確定的參數(shù),此處 α取2,η取-0.15。

    現(xiàn)將修正后的火焰?zhèn)鞑ニ俣却雺毫τ嬎隳P?,則計算結(jié)果如圖7 所示。由圖7 可知,經(jīng)湍流修正后的壓力計算模型得到的峰值壓力及升壓時間與傳感器T1 測得的實驗結(jié)果吻合較好,使用湍流修正因數(shù)預測因傳播截面改變引起的湍流區(qū)域的壓力值具有可行性。

    圖 7 針對體積分數(shù)不同的可燃氣體,修正后的壓力計算模型得到的瀉爆壓力與實驗結(jié)果的比較Fig.7 Comparison of gas explosion venting pressures calculated by the modified model and the test data at different volume fractions of combustible gas

    4 結(jié) 論

    (1)通過將爆燃泄放過程進行微分,假設每個微小時間段內(nèi)爆燃泄放過程均按照先燃燒、再泄放、最后壓力平衡的過程獨立分步進行,得到了泄爆壓力分步計算模型。該模型計算簡單,不僅可以計算爆燃泄放過程中的壓力變化情況,還可以得到各時刻空間內(nèi)的氣體密度、體積等物理狀態(tài)參數(shù)。

    (2)在尺寸為2 m×1.2 m×0.6 m 的長方體容器中充滿乙烯-空氣混合氣體進行實驗驗證,發(fā)現(xiàn)當采用大面積泄爆板作為泄爆構(gòu)件時爆燃壓力曲線均為單峰值曲線,該條件下壓力模型計算結(jié)果與實驗結(jié)果吻合較好。

    (3)當采用小面積泄爆膜作為泄爆構(gòu)件時會因外流氣體流動截面改變引起火焰湍流加速,使得腔體內(nèi)產(chǎn)生壓力梯度,近泄爆口處壓力值大于腔體內(nèi)部。通過湍流加速因子對火焰速度進行修正,則修正后的壓力計算模型可以較好地描述近泄爆口處的壓力變化情況。

    猜你喜歡
    腔體湍流燃氣
    擰緊燃氣“安全閥”
    近期實施的燃氣國家標準
    煤氣與熱力(2022年4期)2022-05-23 12:45:00
    探討燃氣工程的造價控制及跟蹤審計
    高鐵復雜腔體鑄造數(shù)值仿真及控制技術(shù)研究
    高鐵制動系統(tǒng)復雜腔體鑄造成形數(shù)值模擬
    重氣瞬時泄漏擴散的湍流模型驗證
    橡膠擠出裝置
    開孔金屬腔體場強增強效應分析
    “青春期”湍流中的智慧引渡(三)
    “青春期”湍流中的智慧引渡(二)
    日本a在线网址| 亚洲在线自拍视频| 欧美+日韩+精品| h日本视频在线播放| 日本黄大片高清| 婷婷亚洲欧美| 亚洲自拍偷在线| 亚洲美女视频黄频| 亚洲精品国产成人久久av| 久久热精品热| 真人做人爱边吃奶动态| 欧美日韩国产亚洲二区| 成人毛片a级毛片在线播放| av天堂在线播放| www.色视频.com| 18禁在线播放成人免费| 亚洲欧美日韩高清专用| 直男gayav资源| 婷婷丁香在线五月| 99久久精品一区二区三区| 欧美激情久久久久久爽电影| 久久热精品热| 日韩强制内射视频| 国产欧美日韩一区二区精品| 久久久精品大字幕| 一个人看的www免费观看视频| 中文字幕av在线有码专区| 亚洲国产精品合色在线| 欧美又色又爽又黄视频| 欧美性感艳星| 人妻少妇偷人精品九色| 国产又黄又爽又无遮挡在线| 免费人成在线观看视频色| 国产视频一区二区在线看| 此物有八面人人有两片| 99在线人妻在线中文字幕| 国产伦人伦偷精品视频| 成年人黄色毛片网站| 美女高潮的动态| 国产综合懂色| 国产精品一区二区三区四区免费观看 | 日本色播在线视频| 狂野欧美激情性xxxx在线观看| 亚洲国产欧美人成| 一区二区三区免费毛片| 国产精品久久久久久久久免| 亚洲精华国产精华液的使用体验 | 超碰av人人做人人爽久久| 男人的好看免费观看在线视频| 精品久久久久久久末码| 性欧美人与动物交配| 国产色爽女视频免费观看| 窝窝影院91人妻| 欧美人与善性xxx| 久久久久久九九精品二区国产| 一个人看的www免费观看视频| 免费观看在线日韩| 久久久成人免费电影| 一个人免费在线观看电影| 国产亚洲精品久久久com| 久久精品国产鲁丝片午夜精品 | 国产蜜桃级精品一区二区三区| 成人av在线播放网站| 国产亚洲精品久久久com| 成年人黄色毛片网站| 国产成人aa在线观看| 一级av片app| 窝窝影院91人妻| 在线观看免费视频日本深夜| 久久99热这里只有精品18| 狂野欧美激情性xxxx在线观看| АⅤ资源中文在线天堂| 97热精品久久久久久| 美女 人体艺术 gogo| 天堂影院成人在线观看| 亚洲va日本ⅴa欧美va伊人久久| av天堂中文字幕网| av国产免费在线观看| 欧美日韩乱码在线| 99久国产av精品| 精品久久久噜噜| 国产乱人伦免费视频| 一进一出抽搐gif免费好疼| 亚洲专区中文字幕在线| 啦啦啦韩国在线观看视频| 日韩一本色道免费dvd| 欧美激情久久久久久爽电影| 国产精品亚洲美女久久久| 亚洲国产精品合色在线| 真人做人爱边吃奶动态| 国产精品一区二区三区四区免费观看 | 国产爱豆传媒在线观看| 一区二区三区激情视频| 色综合婷婷激情| 天堂av国产一区二区熟女人妻| 国产精品久久久久久av不卡| 欧美一区二区亚洲| 乱人视频在线观看| 在线观看一区二区三区| 欧美又色又爽又黄视频| 精品久久久久久,| 婷婷精品国产亚洲av在线| 亚洲国产精品合色在线| 亚洲色图av天堂| 小说图片视频综合网站| 亚洲美女视频黄频| 欧美成人性av电影在线观看| 干丝袜人妻中文字幕| 99久久精品一区二区三区| 日韩中文字幕欧美一区二区| 成人三级黄色视频| 欧美日韩中文字幕国产精品一区二区三区| 欧美性猛交╳xxx乱大交人| 99在线视频只有这里精品首页| 国产美女午夜福利| 免费看日本二区| 欧美精品啪啪一区二区三区| 午夜影院日韩av| 成年女人永久免费观看视频| 久久九九热精品免费| 真实男女啪啪啪动态图| 亚洲中文字幕一区二区三区有码在线看| 12—13女人毛片做爰片一| 99久久中文字幕三级久久日本| 色吧在线观看| 嫩草影院入口| 免费黄网站久久成人精品| 国产免费男女视频| 国产精品综合久久久久久久免费| 色吧在线观看| 久久久国产成人免费| 最近最新中文字幕大全电影3| 亚洲av中文字字幕乱码综合| 一区二区三区高清视频在线| 天堂√8在线中文| 男人的好看免费观看在线视频| 亚洲欧美激情综合另类| 中国美白少妇内射xxxbb| 97碰自拍视频| 亚洲最大成人av| 春色校园在线视频观看| 欧美日本视频| eeuss影院久久| 欧美性感艳星| 一本精品99久久精品77| 国内精品美女久久久久久| 免费不卡的大黄色大毛片视频在线观看 | 国产亚洲精品综合一区在线观看| 欧美三级亚洲精品| 国产一区二区在线av高清观看| 成人午夜高清在线视频| 欧美国产日韩亚洲一区| 看免费成人av毛片| 亚州av有码| 欧美xxxx黑人xx丫x性爽| 国产一区二区三区在线臀色熟女| 国产蜜桃级精品一区二区三区| 欧美性猛交╳xxx乱大交人| 麻豆精品久久久久久蜜桃| 在线国产一区二区在线| av天堂在线播放| 欧美性猛交╳xxx乱大交人| 午夜激情福利司机影院| 非洲黑人性xxxx精品又粗又长| 男插女下体视频免费在线播放| 国产真实乱freesex| 中文在线观看免费www的网站| 国产欧美日韩精品一区二区| 亚洲专区国产一区二区| 日韩欧美三级三区| 最近中文字幕高清免费大全6 | 乱人视频在线观看| 国产伦人伦偷精品视频| 淫秽高清视频在线观看| 91久久精品国产一区二区三区| 最好的美女福利视频网| 亚洲精品色激情综合| 久久久精品欧美日韩精品| 男插女下体视频免费在线播放| 一个人观看的视频www高清免费观看| 一个人观看的视频www高清免费观看| 欧美性猛交黑人性爽| 一个人免费在线观看电影| 欧美性猛交黑人性爽| 深爱激情五月婷婷| videossex国产| 日韩 亚洲 欧美在线| 亚洲男人的天堂狠狠| 黄色丝袜av网址大全| 免费观看精品视频网站| 少妇熟女aⅴ在线视频| 国产高清视频在线播放一区| 91午夜精品亚洲一区二区三区 | 日本精品一区二区三区蜜桃| 亚洲熟妇熟女久久| av女优亚洲男人天堂| 69人妻影院| 亚洲成人久久爱视频| 久久草成人影院| 精品一区二区三区av网在线观看| 99久久精品国产国产毛片| 久久久国产成人免费| 成人国产综合亚洲| 两性午夜刺激爽爽歪歪视频在线观看| netflix在线观看网站| 99久久久亚洲精品蜜臀av| 亚洲精品在线观看二区| 日韩欧美一区二区三区在线观看| 国内揄拍国产精品人妻在线| 亚洲最大成人av| 国产高清有码在线观看视频| 高清毛片免费观看视频网站| 校园春色视频在线观看| 老司机福利观看| 99久久成人亚洲精品观看| 少妇熟女aⅴ在线视频| 男女啪啪激烈高潮av片| 国产精品免费一区二区三区在线| 99久久中文字幕三级久久日本| 国产91精品成人一区二区三区| 日日撸夜夜添| 免费看美女性在线毛片视频| 亚洲va在线va天堂va国产| 精品午夜福利视频在线观看一区| 91精品国产九色| 欧美最黄视频在线播放免费| 99视频精品全部免费 在线| 日韩精品青青久久久久久| avwww免费| 免费大片18禁| 一个人免费在线观看电影| 色精品久久人妻99蜜桃| 日韩欧美免费精品| 国产伦精品一区二区三区视频9| eeuss影院久久| 欧美一级a爱片免费观看看| 久久99热6这里只有精品| 窝窝影院91人妻| 中文字幕免费在线视频6| 色av中文字幕| 在线免费十八禁| 久99久视频精品免费| 可以在线观看的亚洲视频| 亚洲成人久久爱视频| 在线看三级毛片| 亚洲国产精品久久男人天堂| 人妻丰满熟妇av一区二区三区| 亚洲av成人精品一区久久| 少妇熟女aⅴ在线视频| 国产亚洲精品久久久久久毛片| 免费人成在线观看视频色| 婷婷亚洲欧美| 午夜影院日韩av| 在线观看av片永久免费下载| 身体一侧抽搐| 亚洲专区国产一区二区| 亚洲经典国产精华液单| 欧美潮喷喷水| 中文字幕久久专区| 欧美日韩亚洲国产一区二区在线观看| 少妇猛男粗大的猛烈进出视频 | 国产精品福利在线免费观看| 精品一区二区免费观看| 国产精品久久电影中文字幕| 亚洲第一电影网av| 看片在线看免费视频| 女同久久另类99精品国产91| 久久精品国产自在天天线| 久久久久久久精品吃奶| 久久亚洲真实| 国产探花在线观看一区二区| 国产精品福利在线免费观看| 91久久精品国产一区二区三区| 久久久精品欧美日韩精品| 自拍偷自拍亚洲精品老妇| 熟女电影av网| 尤物成人国产欧美一区二区三区| 国内精品一区二区在线观看| 麻豆久久精品国产亚洲av| 99精品在免费线老司机午夜| 国产av在哪里看| 美女 人体艺术 gogo| 欧美zozozo另类| 亚洲五月天丁香| 九九热线精品视视频播放| 精品国内亚洲2022精品成人| 日本免费a在线| 成人无遮挡网站| 国产伦精品一区二区三区视频9| 少妇的逼水好多| 成人一区二区视频在线观看| 国产色爽女视频免费观看| a级一级毛片免费在线观看| 中文字幕高清在线视频| 看片在线看免费视频| 亚洲国产精品合色在线| 天天一区二区日本电影三级| 尾随美女入室| 国产极品精品免费视频能看的| 国产精品亚洲一级av第二区| 日本与韩国留学比较| 久久亚洲精品不卡| 久久久久久国产a免费观看| 美女xxoo啪啪120秒动态图| 欧美区成人在线视频| 91在线精品国自产拍蜜月| 欧美一区二区精品小视频在线| 在线免费观看不下载黄p国产 | 亚洲精品成人久久久久久| 日韩一本色道免费dvd| 免费观看精品视频网站| 午夜a级毛片| 成人永久免费在线观看视频| 99热只有精品国产| 亚洲av中文av极速乱 | 亚洲精品成人久久久久久| 偷拍熟女少妇极品色| 有码 亚洲区| 亚洲国产欧美人成| 国产成年人精品一区二区| 精品乱码久久久久久99久播| АⅤ资源中文在线天堂| 长腿黑丝高跟| 最近视频中文字幕2019在线8| 久久精品国产自在天天线| 午夜免费成人在线视频| 国产精品福利在线免费观看| 亚洲中文字幕一区二区三区有码在线看| 国产又黄又爽又无遮挡在线| 免费看日本二区| 欧美日韩中文字幕国产精品一区二区三区| 在线国产一区二区在线| 日本一本二区三区精品| 国产欧美日韩精品一区二区| 成人一区二区视频在线观看| 亚洲午夜理论影院| 欧美日韩瑟瑟在线播放| 亚洲色图av天堂| 亚洲中文字幕日韩| 我的老师免费观看完整版| 国产精品人妻久久久久久| 亚洲av中文字字幕乱码综合| 免费人成在线观看视频色| 如何舔出高潮| 亚洲欧美激情综合另类| 亚洲精品456在线播放app | 亚洲图色成人| 干丝袜人妻中文字幕| av黄色大香蕉| 窝窝影院91人妻| 国产精品野战在线观看| 国产毛片a区久久久久| 一区二区三区激情视频| 1000部很黄的大片| 一区福利在线观看| 国产精品,欧美在线| 99久久中文字幕三级久久日本| 亚洲成av人片在线播放无| 久久久久久大精品| 久久99热这里只有精品18| 亚洲 国产 在线| 999久久久精品免费观看国产| 久久久久国产精品人妻aⅴ院| 日韩在线高清观看一区二区三区 | 国产精品1区2区在线观看.| a级一级毛片免费在线观看| 亚洲国产精品合色在线| 熟女电影av网| 国产爱豆传媒在线观看| 搡老岳熟女国产| 综合色av麻豆| 男女视频在线观看网站免费| xxxwww97欧美| 99热精品在线国产| 我的老师免费观看完整版| 日韩精品有码人妻一区| 嫩草影院精品99| 精品久久久久久久末码| 国产69精品久久久久777片| 国产在视频线在精品| 91狼人影院| 欧美三级亚洲精品| 无遮挡黄片免费观看| 我要看日韩黄色一级片| 欧美色视频一区免费| 午夜久久久久精精品| 国产精品国产高清国产av| 亚洲精品成人久久久久久| 免费av毛片视频| 五月玫瑰六月丁香| 色综合站精品国产| 观看美女的网站| 午夜福利在线观看吧| 亚洲人成伊人成综合网2020| 国产高清三级在线| 看片在线看免费视频| 中文亚洲av片在线观看爽| 午夜福利在线在线| 午夜a级毛片| 精品久久国产蜜桃| 最近最新免费中文字幕在线| 久久久国产成人精品二区| 一个人看视频在线观看www免费| 日本免费一区二区三区高清不卡| 色综合站精品国产| 99久国产av精品| 国产老妇女一区| 国产精品自产拍在线观看55亚洲| 成年女人毛片免费观看观看9| 老熟妇乱子伦视频在线观看| 男女边吃奶边做爰视频| 亚洲av不卡在线观看| 一进一出抽搐gif免费好疼| 国产精品久久视频播放| 亚洲无线观看免费| 欧美又色又爽又黄视频| 在现免费观看毛片| 国产高清不卡午夜福利| 免费大片18禁| 国产高清三级在线| 精品人妻1区二区| 久久久久久久精品吃奶| 啦啦啦韩国在线观看视频| 国内揄拍国产精品人妻在线| 最近在线观看免费完整版| 可以在线观看的亚洲视频| 亚洲无线观看免费| 丰满乱子伦码专区| 99精品在免费线老司机午夜| 国内揄拍国产精品人妻在线| 亚洲av熟女| 在线播放无遮挡| 国产免费av片在线观看野外av| 男人和女人高潮做爰伦理| 麻豆精品久久久久久蜜桃| 久久精品91蜜桃| a在线观看视频网站| 国产午夜精品论理片| 国产精品爽爽va在线观看网站| 欧美成人一区二区免费高清观看| 美女cb高潮喷水在线观看| 不卡视频在线观看欧美| 在线观看美女被高潮喷水网站| 日本a在线网址| 国产亚洲精品av在线| 深夜精品福利| 尤物成人国产欧美一区二区三区| 亚洲真实伦在线观看| 久99久视频精品免费| 在线天堂最新版资源| 一本久久中文字幕| 国产精品国产高清国产av| 日韩欧美三级三区| 国产综合懂色| 免费电影在线观看免费观看| 国内久久婷婷六月综合欲色啪| 午夜福利在线在线| 狠狠狠狠99中文字幕| 国产午夜精品论理片| 欧美性猛交╳xxx乱大交人| 午夜视频国产福利| 色哟哟·www| 久久久精品大字幕| 色综合亚洲欧美另类图片| 精品久久久久久久人妻蜜臀av| 日本色播在线视频| 熟女人妻精品中文字幕| 身体一侧抽搐| 一本久久中文字幕| 变态另类成人亚洲欧美熟女| 一区二区三区四区激情视频 | 国产一区二区在线av高清观看| 亚洲不卡免费看| 欧美日韩精品成人综合77777| 国产人妻一区二区三区在| 日本-黄色视频高清免费观看| 69av精品久久久久久| 伦精品一区二区三区| 久久久精品欧美日韩精品| 免费一级毛片在线播放高清视频| 香蕉av资源在线| 欧美国产日韩亚洲一区| 亚洲精品日韩av片在线观看| 国产精品美女特级片免费视频播放器| 三级男女做爰猛烈吃奶摸视频| 欧美日韩乱码在线| 91在线观看av| 噜噜噜噜噜久久久久久91| 免费一级毛片在线播放高清视频| 91麻豆精品激情在线观看国产| 91麻豆av在线| 无遮挡黄片免费观看| 欧美色欧美亚洲另类二区| 久久久午夜欧美精品| 午夜激情欧美在线| 神马国产精品三级电影在线观看| 国内精品宾馆在线| 日韩中文字幕欧美一区二区| 一个人观看的视频www高清免费观看| 日日摸夜夜添夜夜添小说| 国产亚洲91精品色在线| 亚洲在线观看片| 高清日韩中文字幕在线| av黄色大香蕉| 丰满的人妻完整版| av黄色大香蕉| 色吧在线观看| 国产高清三级在线| 亚洲av美国av| 一本精品99久久精品77| 少妇丰满av| 亚洲国产精品久久男人天堂| 少妇人妻一区二区三区视频| 成人二区视频| 最好的美女福利视频网| aaaaa片日本免费| 波多野结衣高清无吗| 美女 人体艺术 gogo| 波多野结衣高清无吗| 观看美女的网站| 国内精品久久久久精免费| 日本 av在线| 91麻豆av在线| 最新中文字幕久久久久| 欧美高清性xxxxhd video| 婷婷六月久久综合丁香| 国产精品久久电影中文字幕| 18禁黄网站禁片免费观看直播| 午夜精品久久久久久毛片777| 国产色爽女视频免费观看| 国产私拍福利视频在线观看| 中文在线观看免费www的网站| 亚洲午夜理论影院| 亚洲精品色激情综合| 中文字幕久久专区| 国产精品日韩av在线免费观看| 99热精品在线国产| 国产一区二区亚洲精品在线观看| 美女黄网站色视频| 啦啦啦啦在线视频资源| 国产精品无大码| 久久精品国产99精品国产亚洲性色| 精品人妻熟女av久视频| 午夜福利在线观看免费完整高清在 | av国产免费在线观看| 午夜福利成人在线免费观看| 国内揄拍国产精品人妻在线| 级片在线观看| 国产精品一区二区三区四区久久| 最近最新中文字幕大全电影3| 欧美+亚洲+日韩+国产| 一级黄片播放器| 午夜福利18| 午夜福利成人在线免费观看| 国产伦精品一区二区三区视频9| 蜜桃久久精品国产亚洲av| 久久精品国产亚洲av香蕉五月| 中文字幕免费在线视频6| 久久久成人免费电影| 成人综合一区亚洲| 亚洲无线在线观看| 人妻制服诱惑在线中文字幕| 嫁个100分男人电影在线观看| 日本黄色视频三级网站网址| 18禁黄网站禁片午夜丰满| 婷婷精品国产亚洲av在线| 国产69精品久久久久777片| 五月伊人婷婷丁香| 日韩欧美在线二视频| 国产精品国产三级国产av玫瑰| 又黄又爽又刺激的免费视频.| 国产精品国产三级国产av玫瑰| 观看免费一级毛片| 国产成人福利小说| 国产精品不卡视频一区二区| 日本a在线网址| 国产精品98久久久久久宅男小说| 搡老妇女老女人老熟妇| 免费av不卡在线播放| 日韩亚洲欧美综合| av天堂在线播放| 特大巨黑吊av在线直播| 99热网站在线观看| 免费在线观看日本一区| 噜噜噜噜噜久久久久久91| 人妻制服诱惑在线中文字幕| 久久午夜福利片| 我的老师免费观看完整版| 久久久久性生活片| 男女下面进入的视频免费午夜| 亚洲国产精品久久男人天堂| 国产成年人精品一区二区| 看黄色毛片网站| 狠狠狠狠99中文字幕| 国产精品亚洲一级av第二区| 简卡轻食公司| 精品久久久久久,| 精品人妻一区二区三区麻豆 | 免费av不卡在线播放| 久久午夜亚洲精品久久| 波多野结衣高清无吗| 欧美成人免费av一区二区三区| 国产精品av视频在线免费观看| 日本三级黄在线观看| 波多野结衣高清作品| 韩国av一区二区三区四区| 亚洲国产精品成人综合色| 亚洲成人免费电影在线观看| 啪啪无遮挡十八禁网站| 国产午夜福利久久久久久| 很黄的视频免费| 精品国产三级普通话版| 黄色欧美视频在线观看| 18禁在线播放成人免费| 黄色日韩在线|