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

    基于流固耦合的脫硫攪拌器振動(dòng)特性分析

    2022-09-21 08:49:02蔣麒麟但斌斌龔昌運(yùn)王偉劉洋都李平
    機(jī)床與液壓 2022年17期
    關(guān)鍵詞:攪拌器鐵水流體

    蔣麒麟,但斌斌,龔昌運(yùn),王偉,劉洋,都李平

    (1.武漢科技大學(xué)冶金裝備及其控制教育部重點(diǎn)實(shí)驗(yàn)室,湖北武漢 430081;2.武漢科技大學(xué)機(jī)械傳動(dòng)與制造工程湖北省重點(diǎn)實(shí)驗(yàn)室,湖北武漢 430081;3.武漢鋼鐵有限公司技術(shù)中心,湖北武漢 430080)

    0 前言

    現(xiàn)代工業(yè)中使用的特種鋼由于其惡劣的服役環(huán)境,對(duì)鋼中的硫含量要求更嚴(yán)格,需要低于0.005%。鐵水脫硫過程是指轉(zhuǎn)爐冶煉前的鐵水預(yù)處理工藝。其中,KR脫硫由于其優(yōu)良的動(dòng)態(tài)條件和高脫硫率,在提高鋼材性能、保證連鑄和鋼坯質(zhì)量、降低預(yù)處理成本等方面具有明顯優(yōu)勢(shì)。鐵水脫硫攪拌器是實(shí)施KR脫硫的重要設(shè)備。關(guān)于KR脫硫攪拌器,劉鵬、李文明、畢學(xué)工、肖林偉等都進(jìn)行了相關(guān)研究。劉鵬應(yīng)用Fluent軟件,結(jié)合多重參照系方法,采用Eulerian多流體模型和標(biāo)準(zhǔn)-湍流模型,研究了KR機(jī)械攪拌脫硫混合器中流體的流動(dòng)特性。李文明應(yīng)用數(shù)值方法研究了KR機(jī)械攪拌和噴吹攪拌的動(dòng)力學(xué)特性,分析了KR機(jī)械攪拌流場(chǎng)和鐵水罐內(nèi)自由液面凹進(jìn)深度的變化規(guī)律以及噴吹攪拌時(shí)氣泡對(duì)流場(chǎng)的影響。畢學(xué)工等通過建立不同工況下的混合模型,利用軟件Fluent仿真得到影響混合效果的最佳工況。肖林偉等通過模擬攪拌器的中后期狀態(tài),根據(jù)脫硫劑的運(yùn)動(dòng)和分布,分析出攪拌器中后期所需的轉(zhuǎn)速。

    上述研究人員主要研究了攪拌器對(duì)脫硫和攪拌的影響,但忽略了KR脫硫攪拌器在脫硫過程中的振動(dòng)問題,沒有分析攪拌器自身在脫硫過程中的振動(dòng)特性。本文作者利用有限元軟件對(duì)基于流固耦合的脫硫攪拌器振動(dòng)特性進(jìn)行數(shù)值模擬,旨在為攪拌器結(jié)構(gòu)的設(shè)計(jì)和優(yōu)化提供參考。

    1 攪拌器數(shù)學(xué)模型的建立及求解

    1.1 流場(chǎng)控制方程

    在實(shí)際生產(chǎn)過程中,鐵水在罐內(nèi)的流動(dòng)非常復(fù)雜,為方便研究,對(duì)模型做以下簡(jiǎn)化和假設(shè):(1)罐內(nèi)鐵水的黏度恒定;(2)忽略溫度場(chǎng)和濃度場(chǎng)對(duì)鐵水物性參數(shù)的影響;(3)考慮重力加速度沿軸向下,大小為9.8 m/s。

    在非定常條件下,罐內(nèi)鐵水的基本方程包括質(zhì)量控制方程與動(dòng)量控制方程。

    質(zhì)量控制方程(連續(xù)方程):

    (1)

    式中:第2、3、4項(xiàng)為質(zhì)量流密度的散度,表示單位時(shí)間內(nèi)通過單位面積的流體質(zhì)量;為流體密度,kg/m;為時(shí)間,s;為流場(chǎng)速度,m/s;、分別為速度矢量在、、方向的分量。

    動(dòng)量控制方程:

    (2)

    式中:為壓力,Pa;為應(yīng)力張量;為包含其他與模型相關(guān)源項(xiàng),如自定義源項(xiàng)。

    1.2 流固耦合方程

    流體與固體模型的網(wǎng)格單元類型可以不同,但在耦合面上必須滿足應(yīng)力、位移、熱流量和溫度等變量守恒,即滿足以下4個(gè)方程:

    =

    (3)

    =

    (4)

    =

    (5)

    =

    (6)

    式中:、分別為流體與固體在耦合面上的位移,m;、分別為流體與固體在耦合面上的剪切應(yīng)力,Pa;、分別為流體和固體在耦合面上的熱流量,W;、分別為流體和固體在耦合面上的溫度,℃。當(dāng)不考慮溫度變化時(shí),應(yīng)同時(shí)滿足位移平衡方程和應(yīng)力平衡方程,即通過滿足這兩個(gè)變量的守恒實(shí)現(xiàn)流體與固體之間分析參數(shù)的傳遞。

    2 攪拌器仿真模型

    2.1 攪拌器三維模型的建立

    根據(jù)我國(guó)某鋼廠KR脫硫工藝,結(jié)合水模型試驗(yàn)中的分析方法,應(yīng)用相似原理,按1∶6建立鐵水罐和攪拌器的數(shù)值模型。攪拌器和鐵水罐所用材料為Q235,密度7 850 kg/m,彈性模量206 GPa,泊松比0.3,質(zhì)量4 522 kg。

    圖1 攪拌器模型

    2.2 邊界條件及計(jì)算

    Gambit軟件作為Fluent的前處理器被廣泛地用于數(shù)值模擬網(wǎng)格劃分,將三維軟件中建立好的三葉常規(guī)攪拌器三維模型導(dǎo)入Gambit中進(jìn)行前處理網(wǎng)格劃分,并在Gambit中將它劃分為3個(gè)區(qū)域:動(dòng)區(qū)域、靜區(qū)域和空氣區(qū)。其中,動(dòng)區(qū)域存在斜邊和倒角,幾何形狀不規(guī)則,使用結(jié)構(gòu)四面體網(wǎng)格,采用Tet/Hybrid網(wǎng)格類型;靜區(qū)域與空氣區(qū)模型幾何形狀規(guī)則,使用結(jié)構(gòu)六面體網(wǎng)格,采用Hex/Wedge網(wǎng)格類型,網(wǎng)格數(shù)目共348 469個(gè)。劃分好的網(wǎng)格如圖2所示。

    圖2 常規(guī)三葉模型網(wǎng)格

    求解模式為有限體積法,流體流動(dòng)為定常流動(dòng)。圖3所示為攪拌器邊界條件設(shè)置,采用VOF模型模擬流場(chǎng)中的氣液兩相流,使用多重參考系法處理動(dòng)區(qū)域與靜區(qū)域的能量傳輸現(xiàn)象,兩相交面設(shè)定為交界面(interface),容器固體壁面為無滑移邊界條件(no slip)。攪拌器葉片為運(yùn)動(dòng)壁面,動(dòng)網(wǎng)格與攪拌器同步轉(zhuǎn)動(dòng),轉(zhuǎn)速為150 r/min;罐壁及下底面設(shè)置為(wall);自由液面初始為靜止?fàn)顟B(tài),液面初始高度設(shè)定為380 mm,攪拌器潛入液面深度為300 mm。采用Simple算法對(duì)壓力-速度耦合進(jìn)行求解,離散格式采用一階迎風(fēng),所有項(xiàng)的殘差收斂范圍均為1×10。時(shí)間步長(zhǎng)設(shè)定為0.005 s,時(shí)間步設(shè)定為2 000步,共分析時(shí)長(zhǎng)為10 s。

    圖3 攪拌器邊界條件設(shè)置

    3 仿真計(jì)算結(jié)果分析

    3.1 流固耦合結(jié)果分析

    由圖4可以看出:攪拌器附近的區(qū)域流速較大,鐵水罐壁及攪拌器上方區(qū)域流速較小。這是由于鐵水罐內(nèi)流體的速度來源于攪拌器的轉(zhuǎn)動(dòng),攪拌器的轉(zhuǎn)動(dòng)帶動(dòng)動(dòng)區(qū)域鐵水作旋轉(zhuǎn)運(yùn)動(dòng),使得最先獲得速度的是動(dòng)區(qū)域附近的流體,該區(qū)域的流體最大流速為1.2 m/s。隨著連續(xù)攪拌,動(dòng)區(qū)域流體的動(dòng)能向外擴(kuò)散以傳遞能量,當(dāng)?shù)竭_(dá)鐵水罐壁時(shí),由于鐵水罐壁的阻礙和吸附,使得在速度傳遞過程中損失了大量能量,在鐵水罐壁處速度為0.1 m/s。由圖5可以看出:攪拌頭底部的流體壓力較大,而攪拌軸附近的流體壓力較小,底部最大流體壓力約為1 350 Pa,攪拌軸附近最小流體壓力為750 Pa。這是因?yàn)楦鶕?jù)能量守恒原理在攪拌流場(chǎng)中的應(yīng)用,流體速度大的區(qū)域壓力較小,速度小的區(qū)域壓力較大。在攪拌過程中,由于鐵水罐底部流體速度小,底部附近的壓力大,攪拌頭動(dòng)區(qū)域附近流體流速大,攪拌頭動(dòng)區(qū)域附近壓力小。

    圖4 x截面速度分布

    圖5 x-y平面壓力分布

    由圖6可以看出:攪拌器的最大變形發(fā)生在攪拌葉片底部區(qū)域,而最小變形發(fā)生攪拌器軸的上部。葉片底部區(qū)域壓力大流速小,以至于變形量大,磨損嚴(yán)重。攪拌軸上部放置在空氣中,沒有流體與之接觸,因此變形最小。由圖7可知:最大等效應(yīng)力發(fā)生在攪拌軸與攪拌葉片的連接處,數(shù)值為2 299.6 Pa,此處極易產(chǎn)生應(yīng)力集中;最小等效應(yīng)力出現(xiàn)在攪拌軸上部,數(shù)值為14.219 Pa,此處攪拌軸置于空氣,不受流場(chǎng)的壓力作用,所以該處的等效應(yīng)力最小。

    圖6 攪拌器變形分布云圖 圖7 攪拌器應(yīng)力分布云圖

    3.2 預(yù)應(yīng)力模態(tài)分析

    在結(jié)構(gòu)振動(dòng)中,低階模態(tài)占主導(dǎo)地位,因此選擇攪拌器的前4階模態(tài)進(jìn)行觀測(cè),結(jié)果如圖8和表1所示??梢钥闯觯涸?~2 600 Hz頻率范圍內(nèi),最先發(fā)生變形的部位主要在攪拌器葉片處,該部位由于強(qiáng)大的流體壓力而產(chǎn)生較大的變形。

    由圖8和表1可知:1階模態(tài)的固有頻率為253.19 Hz,其振型主要是攪拌器沿軸正向的彎曲變形,變形從攪拌器頂部至葉片底部逐漸增加,在葉片底部達(dá)到最大變形為12.689 mm;2階模態(tài)的固有頻率為901.62 Hz,其振型主要是攪拌器繞軸的扭轉(zhuǎn)變形,攪拌器整體變形大,變形從攪拌器頂部至底部逐漸增大,在葉片端面處達(dá)到最大變形17.73 mm;3階模態(tài)的固有頻率為1 818.90 Hz,振型主要是攪拌器沿軸負(fù)向有明顯的彎曲變形,由中間向兩端變形逐漸減?。?階模態(tài)的固有頻率為2 552.10 Hz,振型主要是攪拌器以底端方向?yàn)橹行妮S線產(chǎn)生扭轉(zhuǎn)變形,攪拌器葉片處達(dá)到最大變形9.661 7 mm。

    圖8 攪拌器前4階模態(tài)振型

    表1 攪拌器前4階固有頻率

    3.3 攪拌器諧響應(yīng)分析

    攪拌過程中,攪拌器處于溫度場(chǎng)、流場(chǎng)、力場(chǎng)等多物理場(chǎng)耦合作用的復(fù)雜服役環(huán)境下,會(huì)受到非定常壓力場(chǎng)及不平衡的激勵(lì),從而引起不平衡響應(yīng)。因此,通過諧響應(yīng)分析來計(jì)算攪拌器的不平衡響應(yīng)。

    根據(jù)動(dòng)平衡精度理論:

    (7)

    代入攪拌器設(shè)計(jì)參數(shù):動(dòng)平衡精度=2.5 mm/s,設(shè)計(jì)工作轉(zhuǎn)速=150 r/min,系統(tǒng)總質(zhì)量=4 522 kg,校正半徑=208 mm,得攪拌器的最大剩余不平衡量為3 460.08 g。

    偏心距為1.45 μm,攪拌器因質(zhì)心偏心所產(chǎn)生的力為

    =

    (8)

    式中:為不平衡質(zhì)量;為轉(zhuǎn)子的轉(zhuǎn)速。

    =cos

    (9)

    =sin

    (10)

    將這個(gè)力近似為2個(gè)相互垂直方向上的簡(jiǎn)諧力,兩者之間的相位相差為90°。采用模態(tài)疊加法,選取攪拌器軸與葉片結(jié)合處點(diǎn)作為關(guān)鍵監(jiān)測(cè)點(diǎn),分析攪拌器的位移、加速度、應(yīng)力和應(yīng)變與激勵(lì)頻率之間的響應(yīng)關(guān)系。

    在模態(tài)分析得到的邊界條件下,攪拌器激勵(lì)頻率求解范圍設(shè)置為0~2 600 Hz,運(yùn)算次數(shù)為260,即每次求解間隔10 Hz。在攪拌器點(diǎn)處2個(gè)垂直方向上添加簡(jiǎn)諧激振力,幅值為10 N、相位分別為0°和90°。通過對(duì)交匯點(diǎn)的振動(dòng)諧響應(yīng)分析,得到、、方向上位移、加速度、應(yīng)力及應(yīng)變與激勵(lì)頻率之間的響應(yīng)曲線,如圖9所示。

    圖9 振動(dòng)特性響應(yīng)曲線

    由圖9可知:激振力在攪拌器、、3個(gè)方向引起的響應(yīng)趨勢(shì)基本一致,隨著激勵(lì)頻率的增加,交匯處點(diǎn)的位移、應(yīng)力與應(yīng)變都會(huì)在1階(約260 Hz)和3階(約1 800 Hz)固有頻率附近出現(xiàn)峰值。由圖9(a)可知:激勵(lì)頻率為260 Hz時(shí),方向和方向的位移峰值分別達(dá)到0.001 6、0.000 3 mm;激勵(lì)頻率為1 800 Hz時(shí),方向和方向的位移峰值分別達(dá)到0.006 3 mm和0.002 2 mm。由圖9(b)可知:激勵(lì)頻率為1 800 Hz時(shí),方向上的加速度峰值為8×10mm/s,方向上的加速度峰值為2.3×10mm/s。在攪拌器系統(tǒng)設(shè)計(jì)和結(jié)構(gòu)優(yōu)化中,可以在和方向上進(jìn)行約束,以減小振動(dòng)。由圖9(c)(d)可知:由于存在、方向上的流體沖擊,點(diǎn)位置處方向上應(yīng)力響應(yīng)最大、方向上的應(yīng)變響應(yīng)最大,激勵(lì)頻率為1 800 Hz時(shí),方向上的應(yīng)力達(dá)到9.5 MPa,方向上的應(yīng)變達(dá)到1.5×10mm/mm。攪拌器軸和葉片結(jié)合處附近可以作為振動(dòng)傳感器安裝位置,對(duì)攪拌器進(jìn)行壓力與速度的振動(dòng)數(shù)據(jù)監(jiān)測(cè),以期為攪拌器結(jié)構(gòu)設(shè)計(jì)與優(yōu)化提供參考。

    4 結(jié)論

    本文作者采用單向流固耦合的方法研究了攪拌器的振動(dòng)特性,通過分析數(shù)值結(jié)果發(fā)現(xiàn):

    (1)動(dòng)區(qū)域附近的流體流速最大,數(shù)值為1.2 m/s,在混合罐的底部死區(qū)流速最小,數(shù)值為0.1 m/s;攪拌器底部的流體壓力較大,而攪拌軸附近的流體壓力較?。蛔畲罅黧w壓力為1 350 Pa,最小流體壓力為750 Pa;

    (2)流固耦合作用下最大等效應(yīng)力出現(xiàn)在攪拌軸與攪拌葉片連接處,數(shù)值為2 299.6 Pa;最小等效應(yīng)力出現(xiàn)在攪拌軸上部置于空氣部分,數(shù)值為14.219 Pa;

    (3)通過研究交匯點(diǎn)的振動(dòng)特性,發(fā)現(xiàn)由于內(nèi)部流體的脈動(dòng),當(dāng)攪拌器的激勵(lì)頻率為260 Hz時(shí),和方向的位移峰值分別達(dá)到0.001 6、0.000 3 mm;當(dāng)攪拌器的激勵(lì)頻率為1 800 Hz時(shí),和方向的位移峰值分別達(dá)到0.006 3、0.002 2 mm,和方向的加速度峰值分別達(dá)到8×10、2.3×10mm/s。

    猜你喜歡
    攪拌器鐵水流體
    流體壓強(qiáng)知多少
    山鋼5100m3高爐提高鐵水裝準(zhǔn)率的實(shí)踐
    山東冶金(2022年1期)2022-04-19 13:40:16
    山雨欲來風(fēng)滿樓之流體壓強(qiáng)與流速
    低鐵比條件下低硅、低溫鐵水加煤塊冶煉生產(chǎn)實(shí)踐
    山東冶金(2019年5期)2019-11-16 09:09:06
    解析攪拌器在原油儲(chǔ)罐中的應(yīng)用及發(fā)展
    靈活變化的攪拌器
    等效流體體積模量直接反演的流體識(shí)別方法
    基于Bootstrap的高爐鐵水硅含量二維預(yù)報(bào)
    盧成:發(fā)展集裝箱鐵水聯(lián)運(yùn)大有可為
    專用汽車(2016年8期)2016-03-01 04:15:13
    賽萊默發(fā)布一款全新的高效低速潛污攪拌器
    亚洲精华国产精华液的使用体验 | 麻豆国产97在线/欧美| 亚洲天堂国产精品一区在线| 亚洲精品国产成人久久av| 久久久久久久久久黄片| 国产极品精品免费视频能看的| 国产主播在线观看一区二区| 搡女人真爽免费视频火全软件 | 美女xxoo啪啪120秒动态图| 久久精品国产亚洲网站| 精品一区二区免费观看| 神马国产精品三级电影在线观看| 国产白丝娇喘喷水9色精品| 内地一区二区视频在线| 99热6这里只有精品| 国产精品爽爽va在线观看网站| 久久久久久九九精品二区国产| 久久久精品大字幕| 国产日本99.免费观看| 亚洲av熟女| 亚洲最大成人av| 国产国拍精品亚洲av在线观看| 国产精品久久久久久亚洲av鲁大| 观看免费一级毛片| 国产亚洲av嫩草精品影院| 欧美日韩精品成人综合77777| 免费av不卡在线播放| 亚洲国产精品合色在线| 国模一区二区三区四区视频| 欧美3d第一页| 国产综合懂色| 少妇熟女aⅴ在线视频| 国产一区二区三区视频了| 国产午夜精品久久久久久一区二区三区 | 男女做爰动态图高潮gif福利片| 国产一区二区亚洲精品在线观看| 免费av毛片视频| 日本与韩国留学比较| 欧美日本亚洲视频在线播放| 亚洲精品在线观看二区| 欧美区成人在线视频| 九色国产91popny在线| 内地一区二区视频在线| 夜夜夜夜夜久久久久| 欧美xxxx性猛交bbbb| 啦啦啦观看免费观看视频高清| 国产亚洲精品久久久久久毛片| 男女做爰动态图高潮gif福利片| 亚洲黑人精品在线| 中文在线观看免费www的网站| 国语自产精品视频在线第100页| 3wmmmm亚洲av在线观看| 一区二区三区免费毛片| 又黄又爽又刺激的免费视频.| 麻豆成人av在线观看| 日韩欧美一区二区三区在线观看| 99久久精品国产国产毛片| 国产在视频线在精品| 国产精品1区2区在线观看.| 少妇丰满av| 69人妻影院| 国产69精品久久久久777片| 久久久久久久久久久丰满 | 亚洲欧美日韩高清在线视频| 人妻丰满熟妇av一区二区三区| 免费看a级黄色片| 在线观看66精品国产| 亚洲精品影视一区二区三区av| 一进一出抽搐gif免费好疼| 中文字幕久久专区| 又爽又黄a免费视频| 亚洲av成人精品一区久久| 国产探花极品一区二区| 看十八女毛片水多多多| 亚洲国产精品sss在线观看| 色哟哟·www| 亚洲精品粉嫩美女一区| 毛片一级片免费看久久久久 | 天堂√8在线中文| 我要搜黄色片| 在线a可以看的网站| 亚洲精品乱码久久久v下载方式| 欧美成人a在线观看| 亚洲av美国av| 十八禁网站免费在线| 国产高清激情床上av| 国产一区二区三区在线臀色熟女| 国产亚洲精品久久久com| x7x7x7水蜜桃| a级一级毛片免费在线观看| 国产白丝娇喘喷水9色精品| 国产午夜精品久久久久久一区二区三区 | 久久久久久久午夜电影| 乱系列少妇在线播放| 亚洲人成网站在线播| 日本在线视频免费播放| 国产精品一区二区三区四区免费观看 | 成人无遮挡网站| 一进一出抽搐动态| 精品人妻一区二区三区麻豆 | 国产亚洲精品av在线| 日日摸夜夜添夜夜添小说| 久久久久久久久久久丰满 | 国产成人aa在线观看| 黄色丝袜av网址大全| 老司机深夜福利视频在线观看| 久久热精品热| 免费不卡的大黄色大毛片视频在线观看 | 亚洲精品色激情综合| 一夜夜www| 亚洲成人久久爱视频| 一本久久中文字幕| 有码 亚洲区| 久久欧美精品欧美久久欧美| 一本精品99久久精品77| 欧美国产日韩亚洲一区| 精品久久久久久久久久久久久| 国产精品久久久久久久久免| 午夜福利成人在线免费观看| netflix在线观看网站| netflix在线观看网站| 午夜福利高清视频| 亚洲真实伦在线观看| 国产精品人妻久久久影院| 欧美成人a在线观看| 男人狂女人下面高潮的视频| 久久精品国产亚洲av涩爱 | 免费观看的影片在线观看| 成人亚洲精品av一区二区| 国产免费男女视频| 亚洲成人精品中文字幕电影| 久久精品国产自在天天线| 国产精品爽爽va在线观看网站| 国产人妻一区二区三区在| 国产大屁股一区二区在线视频| 亚洲国产精品合色在线| 久久精品国产自在天天线| 男女之事视频高清在线观看| 国产亚洲av嫩草精品影院| 亚洲不卡免费看| 国产精品一区二区三区四区免费观看 | 联通29元200g的流量卡| 国产一级毛片七仙女欲春2| 身体一侧抽搐| 亚洲av免费在线观看| 少妇裸体淫交视频免费看高清| 欧美zozozo另类| 干丝袜人妻中文字幕| 午夜日韩欧美国产| 一区二区三区高清视频在线| 国产精品一区二区三区四区久久| 国产亚洲精品久久久com| 久久久久精品国产欧美久久久| 少妇高潮的动态图| 欧美+亚洲+日韩+国产| 亚洲内射少妇av| 久久久色成人| 日韩欧美在线乱码| 91午夜精品亚洲一区二区三区 | 中文字幕av成人在线电影| 2021天堂中文幕一二区在线观| 国产精品日韩av在线免费观看| 国产淫片久久久久久久久| 久久久久久久亚洲中文字幕| 一区二区三区激情视频| 哪里可以看免费的av片| 三级毛片av免费| 国产欧美日韩精品一区二区| 久久热精品热| 成人无遮挡网站| 成人综合一区亚洲| 国产精品99久久久久久久久| 丰满的人妻完整版| 很黄的视频免费| 少妇高潮的动态图| 淫秽高清视频在线观看| 色尼玛亚洲综合影院| 国产免费一级a男人的天堂| 国产一区二区三区在线臀色熟女| 真实男女啪啪啪动态图| 2021天堂中文幕一二区在线观| 久久天躁狠狠躁夜夜2o2o| 女人被狂操c到高潮| 久久久久精品国产欧美久久久| 91在线观看av| 国产精品无大码| 色精品久久人妻99蜜桃| 男女那种视频在线观看| 国产高潮美女av| 嫩草影院精品99| 中出人妻视频一区二区| 亚洲精品一卡2卡三卡4卡5卡| 日日撸夜夜添| 精品一区二区三区视频在线| 18+在线观看网站| 日本色播在线视频| 国产精品久久久久久av不卡| 日韩欧美国产在线观看| 国产探花极品一区二区| .国产精品久久| 精品不卡国产一区二区三区| 99国产极品粉嫩在线观看| 国产成人福利小说| 成年人黄色毛片网站| 又黄又爽又刺激的免费视频.| 国产亚洲精品久久久久久毛片| 丰满人妻一区二区三区视频av| 在线观看美女被高潮喷水网站| 久久久久久久久久久丰满 | а√天堂www在线а√下载| 国产高清三级在线| 男女做爰动态图高潮gif福利片| 国产成人av教育| 亚洲av五月六月丁香网| 乱人视频在线观看| 老女人水多毛片| 欧美+亚洲+日韩+国产| 一区二区三区免费毛片| 日本熟妇午夜| 欧美一区二区国产精品久久精品| 国产在线精品亚洲第一网站| 午夜a级毛片| 狂野欧美白嫩少妇大欣赏| 欧美不卡视频在线免费观看| 九九热线精品视视频播放| 国产av一区在线观看免费| 日日啪夜夜撸| 精品欧美国产一区二区三| 欧美精品啪啪一区二区三区| 九九热线精品视视频播放| 国产av一区在线观看免费| 日日啪夜夜撸| 午夜爱爱视频在线播放| 九九在线视频观看精品| 欧美成人性av电影在线观看| 亚洲天堂国产精品一区在线| 成人国产麻豆网| bbb黄色大片| 国产乱人伦免费视频| 久久久久久久午夜电影| 亚洲人成网站高清观看| 草草在线视频免费看| 午夜激情欧美在线| 成人av一区二区三区在线看| 日韩av在线大香蕉| 色5月婷婷丁香| 波多野结衣高清作品| 成人亚洲精品av一区二区| 大型黄色视频在线免费观看| 禁无遮挡网站| 精品不卡国产一区二区三区| 夜夜爽天天搞| 看片在线看免费视频| 久久这里只有精品中国| 在线看三级毛片| 日本-黄色视频高清免费观看| 亚洲七黄色美女视频| 看片在线看免费视频| 国产私拍福利视频在线观看| 国产视频内射| 久久久久久九九精品二区国产| 久久天躁狠狠躁夜夜2o2o| 91狼人影院| 精品久久久久久久久av| av在线老鸭窝| 黄色女人牲交| 内地一区二区视频在线| 日韩中字成人| 很黄的视频免费| 亚洲第一区二区三区不卡| 99国产精品一区二区蜜桃av| 中文资源天堂在线| 天堂√8在线中文| 在线播放国产精品三级| 亚洲美女视频黄频| 麻豆国产av国片精品| av专区在线播放| 亚洲最大成人中文| 在线观看av片永久免费下载| 国产高清视频在线观看网站| 国内精品宾馆在线| 日日摸夜夜添夜夜添小说| 亚洲美女黄片视频| 夜夜爽天天搞| 久久久久久久亚洲中文字幕| 动漫黄色视频在线观看| 哪里可以看免费的av片| 99热这里只有是精品在线观看| 一个人看视频在线观看www免费| 欧美bdsm另类| 最近最新免费中文字幕在线| 国产色婷婷99| 男人舔奶头视频| 一边摸一边抽搐一进一小说| 看免费成人av毛片| 亚洲七黄色美女视频| 国内精品久久久久精免费| 久久亚洲精品不卡| 免费黄网站久久成人精品| 国产午夜精品久久久久久一区二区三区 | a级毛片免费高清观看在线播放| 乱码一卡2卡4卡精品| 91在线观看av| 最近视频中文字幕2019在线8| 99精品在免费线老司机午夜| 哪里可以看免费的av片| 中文字幕久久专区| 久久久久久伊人网av| 午夜福利成人在线免费观看| 欧美bdsm另类| 亚洲男人的天堂狠狠| 国产真实伦视频高清在线观看 | 欧美bdsm另类| 亚洲avbb在线观看| 热99在线观看视频| 亚洲图色成人| 黄色女人牲交| av黄色大香蕉| 色噜噜av男人的天堂激情| 欧美3d第一页| 亚洲精品日韩av片在线观看| 老司机午夜福利在线观看视频| 亚洲 国产 在线| 男人和女人高潮做爰伦理| 国产亚洲精品久久久久久毛片| 免费搜索国产男女视频| 国产精品一区二区三区四区免费观看 | 欧美在线一区亚洲| 欧美精品啪啪一区二区三区| 亚洲七黄色美女视频| 国产伦一二天堂av在线观看| 午夜亚洲福利在线播放| 特级一级黄色大片| 91在线精品国自产拍蜜月| 日日摸夜夜添夜夜添av毛片 | 亚洲国产精品合色在线| 欧美在线一区亚洲| 国语自产精品视频在线第100页| 成人午夜高清在线视频| 人妻制服诱惑在线中文字幕| 美女xxoo啪啪120秒动态图| 中文字幕av在线有码专区| 亚洲av成人av| 国产在线精品亚洲第一网站| 国产亚洲91精品色在线| 国内精品久久久久精免费| 美女高潮喷水抽搐中文字幕| 亚洲av成人精品一区久久| 丰满人妻一区二区三区视频av| 尤物成人国产欧美一区二区三区| 国内揄拍国产精品人妻在线| 内地一区二区视频在线| 免费av毛片视频| 高清在线国产一区| 国产精品国产三级国产av玫瑰| 精品人妻一区二区三区麻豆 | 成人特级av手机在线观看| 啦啦啦观看免费观看视频高清| 五月玫瑰六月丁香| 国产精品一区二区免费欧美| 久久人人爽人人爽人人片va| 免费一级毛片在线播放高清视频| 噜噜噜噜噜久久久久久91| 看片在线看免费视频| 国国产精品蜜臀av免费| 中文字幕av在线有码专区| 亚洲av日韩精品久久久久久密| 麻豆一二三区av精品| 久久人人爽人人爽人人片va| 精品久久久久久久人妻蜜臀av| 久久久久久国产a免费观看| 国产精品野战在线观看| 中文字幕熟女人妻在线| 特级一级黄色大片| 91久久精品国产一区二区成人| 少妇高潮的动态图| 亚洲精华国产精华精| 成人高潮视频无遮挡免费网站| 一个人免费在线观看电影| 国产成人一区二区在线| 99国产精品一区二区蜜桃av| 他把我摸到了高潮在线观看| 亚洲天堂国产精品一区在线| 亚洲av美国av| 国产伦精品一区二区三区四那| 色综合站精品国产| 国产成人一区二区在线| 99国产精品一区二区蜜桃av| 精品不卡国产一区二区三区| 日韩亚洲欧美综合| 成人毛片a级毛片在线播放| 久久久国产成人免费| 我的老师免费观看完整版| 在现免费观看毛片| 1000部很黄的大片| 亚洲人成伊人成综合网2020| 美女大奶头视频| 又爽又黄无遮挡网站| 色在线成人网| 人妻久久中文字幕网| 欧美bdsm另类| 国产 一区 欧美 日韩| 亚洲精品色激情综合| 啦啦啦韩国在线观看视频| 变态另类丝袜制服| 尾随美女入室| 狠狠狠狠99中文字幕| 国产白丝娇喘喷水9色精品| 一区二区三区免费毛片| 三级国产精品欧美在线观看| 18禁黄网站禁片午夜丰满| 人人妻,人人澡人人爽秒播| 日本撒尿小便嘘嘘汇集6| 亚洲精品一区av在线观看| 国产亚洲欧美98| 国产真实伦视频高清在线观看 | 欧美精品国产亚洲| 中文字幕av在线有码专区| 深爱激情五月婷婷| 亚洲国产精品合色在线| 国产成人影院久久av| 精品一区二区三区人妻视频| 无人区码免费观看不卡| 一级a爱片免费观看的视频| 人妻久久中文字幕网| 久久久成人免费电影| 日本一本二区三区精品| 亚洲熟妇中文字幕五十中出| 麻豆一二三区av精品| 18禁在线播放成人免费| 国产白丝娇喘喷水9色精品| 俺也久久电影网| 亚洲中文字幕一区二区三区有码在线看| 麻豆成人午夜福利视频| 国产精品福利在线免费观看| 国产久久久一区二区三区| 九九久久精品国产亚洲av麻豆| 亚洲色图av天堂| 俄罗斯特黄特色一大片| 日韩国内少妇激情av| 久久亚洲真实| 国产午夜福利久久久久久| 国产美女午夜福利| 色在线成人网| 99久久成人亚洲精品观看| 99热只有精品国产| 免费一级毛片在线播放高清视频| 国产一区二区亚洲精品在线观看| 啦啦啦观看免费观看视频高清| 我的女老师完整版在线观看| 国产日本99.免费观看| 夜夜看夜夜爽夜夜摸| 免费看美女性在线毛片视频| 中文字幕高清在线视频| 哪里可以看免费的av片| 国产在线精品亚洲第一网站| 国产精品一及| 亚洲精品乱码久久久v下载方式| 亚洲乱码一区二区免费版| 人妻少妇偷人精品九色| 国产视频一区二区在线看| 久久6这里有精品| av在线亚洲专区| 变态另类丝袜制服| 久久精品国产鲁丝片午夜精品 | 香蕉av资源在线| 人妻久久中文字幕网| 村上凉子中文字幕在线| 少妇人妻精品综合一区二区 | 欧美成人a在线观看| 欧美日韩瑟瑟在线播放| 国产一区二区三区在线臀色熟女| 窝窝影院91人妻| 久久人人爽人人爽人人片va| 日本一本二区三区精品| 亚洲一级一片aⅴ在线观看| 久久草成人影院| 村上凉子中文字幕在线| 又爽又黄无遮挡网站| 国产蜜桃级精品一区二区三区| 国产成人a区在线观看| 国产一区二区三区在线臀色熟女| 亚洲色图av天堂| 亚洲av免费高清在线观看| 日本五十路高清| 校园春色视频在线观看| 国产高清视频在线观看网站| 亚洲,欧美,日韩| 一区二区三区免费毛片| 男女啪啪激烈高潮av片| 免费黄网站久久成人精品| 国产男人的电影天堂91| 一个人看视频在线观看www免费| 18禁裸乳无遮挡免费网站照片| 亚洲 国产 在线| 午夜福利成人在线免费观看| 亚洲av成人av| 久久精品国产亚洲av涩爱 | 久久亚洲真实| 国内少妇人妻偷人精品xxx网站| 老司机午夜福利在线观看视频| 久久中文看片网| 赤兔流量卡办理| 啪啪无遮挡十八禁网站| 22中文网久久字幕| 亚洲久久久久久中文字幕| 99在线人妻在线中文字幕| 欧美绝顶高潮抽搐喷水| 成年人黄色毛片网站| 亚洲自拍偷在线| 一本精品99久久精品77| 欧美性猛交╳xxx乱大交人| 免费无遮挡裸体视频| 国产真实乱freesex| 无遮挡黄片免费观看| 国产色爽女视频免费观看| 黄色视频,在线免费观看| 麻豆成人av在线观看| 中文字幕精品亚洲无线码一区| 亚洲美女视频黄频| 中文字幕高清在线视频| 欧美xxxx黑人xx丫x性爽| 成人毛片a级毛片在线播放| 女的被弄到高潮叫床怎么办 | 久久久久免费精品人妻一区二区| 亚洲精品影视一区二区三区av| 日韩人妻高清精品专区| 国产久久久一区二区三区| 国国产精品蜜臀av免费| 国产精品98久久久久久宅男小说| 三级毛片av免费| 九九久久精品国产亚洲av麻豆| 波多野结衣巨乳人妻| 中亚洲国语对白在线视频| 亚洲精品亚洲一区二区| 校园春色视频在线观看| 琪琪午夜伦伦电影理论片6080| 看十八女毛片水多多多| 国产主播在线观看一区二区| 波多野结衣巨乳人妻| 嫩草影院入口| 欧美又色又爽又黄视频| 老师上课跳d突然被开到最大视频| 国产成人福利小说| h日本视频在线播放| 少妇人妻精品综合一区二区 | 日韩中文字幕欧美一区二区| 日本撒尿小便嘘嘘汇集6| 美女黄网站色视频| 搞女人的毛片| 亚洲国产色片| 一级a爱片免费观看的视频| 精品国内亚洲2022精品成人| 狠狠狠狠99中文字幕| 一进一出抽搐gif免费好疼| 99热这里只有精品一区| 亚洲精品国产成人久久av| 国产精品亚洲一级av第二区| 日本撒尿小便嘘嘘汇集6| 97超视频在线观看视频| 91狼人影院| 久久午夜亚洲精品久久| 18禁黄网站禁片免费观看直播| 免费人成在线观看视频色| 亚洲欧美日韩卡通动漫| 又粗又爽又猛毛片免费看| 亚洲图色成人| 俄罗斯特黄特色一大片| 色哟哟哟哟哟哟| 精品99又大又爽又粗少妇毛片 | 国产黄片美女视频| 无人区码免费观看不卡| 我的女老师完整版在线观看| 欧美绝顶高潮抽搐喷水| 国产人妻一区二区三区在| 99久久精品热视频| 亚洲性久久影院| 最近在线观看免费完整版| 国产主播在线观看一区二区| 国产精品不卡视频一区二区| 俺也久久电影网| av国产免费在线观看| 午夜视频国产福利| 久久久久九九精品影院| 深夜a级毛片| av专区在线播放| 久久久久国内视频| 日本色播在线视频| 真人做人爱边吃奶动态| 久久久久九九精品影院| 亚洲成人中文字幕在线播放| 免费一级毛片在线播放高清视频| 在线观看午夜福利视频| 亚洲美女黄片视频| 可以在线观看毛片的网站| 一边摸一边抽搐一进一小说| 亚洲精品日韩av片在线观看| 免费一级毛片在线播放高清视频| 午夜老司机福利剧场| 一级a爱片免费观看的视频| 成人精品一区二区免费| 毛片女人毛片| 国产主播在线观看一区二区| 国产午夜精品论理片| 99精品久久久久人妻精品| 国产精品久久久久久精品电影| 搞女人的毛片| bbb黄色大片| 亚洲av免费在线观看| 亚洲色图av天堂| 麻豆av噜噜一区二区三区| 九色成人免费人妻av| 亚洲人成网站在线播放欧美日韩| 欧美成人免费av一区二区三区| 狂野欧美激情性xxxx在线观看|