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

    隧道窯冷卻帶流-固耦合散熱數(shù)值模擬

    2016-05-31 01:06:32許寧王雪峰張柱銀黃建軍鹽城工學院機械工程學院江蘇鹽城05湖南工業(yè)大學機械工程學院湖南株洲007中國建材西安墻體材料研究設計院陜西西安7006中國建筑材料集團有限公司北京0006
    中國陶瓷工業(yè) 2016年2期
    關(guān)鍵詞:散熱數(shù)值模擬

    許寧,王雪峰,張柱銀,王 君,黃建軍,肖 慧,郅 曉(. 鹽城工學院機械工程學院,江蘇 鹽城 05;. 湖南工業(yè)大學機械工程學院,湖南 株洲 007;. 中國建材西安墻體材料研究設計院,陜西 西安 7006;. 中國建筑材料集團有限公司,北京 0006)

    ?

    隧道窯冷卻帶流-固耦合散熱數(shù)值模擬

    許寧1, 2,王雪峰1,2,張柱銀2,王 君2,黃建軍1,肖 慧3,郅 曉4
    (1. 鹽城工學院機械工程學院,江蘇 鹽城 224051;2. 湖南工業(yè)大學機械工程學院,湖南 株洲 412007;3. 中國建材西安墻體材料研究設計院,陜西 西安 710061;4. 中國建筑材料集團有限公司,北京 100036)

    摘 要:以某燒結(jié)磚隧道窯為研究對象,為獲得隧道窯冷卻帶散熱溫度場數(shù)據(jù),建立了隧道窯窯壁-磚塊-流體耦合傳熱物理模型與數(shù)學模型,采用FLUENT流體計算軟件對其進行了隧道窯窯壁-磚塊-流體耦合傳熱仿真,獲得了窯內(nèi)整體比較精確直觀的溫度場分布。仿真的結(jié)果表明:風速為5 m/s、孔隙率為0.5以及磚塊適當稀碼時,隧道窯冷卻帶降溫效果較好,其與隧道窯燒成曲線冷卻段吻合也較好,為風速調(diào)整、碼磚方式和孔隙率的選擇提供了理論依據(jù),提高了燒成質(zhì)量和生產(chǎn)效率。

    關(guān)鍵詞:隧道窯冷卻帶;流-固耦合;散熱;數(shù)值模擬

    E-mail:xuning196402@163.com

    0 引 言

    隧道窯是現(xiàn)代化的連續(xù)式燒成的熱工設備,廣泛用于磚塊和陶瓷產(chǎn)品的焙燒生產(chǎn),燒結(jié)磚隧道窯冷卻帶在燒成帶之后,主要的作用是使焙燒好的磚塊迅速冷卻。經(jīng)常在窯尾安裝風機,通過改變風速的大小從而控制磚塊冷卻的速率與溫度。磚塊內(nèi)部傳熱機理是比較復雜的,涉及流相及固相的熱傳導和對流換熱[1]。磚塊上下溫差受冷卻風影響較大,出磚時磚塊表面容易產(chǎn)生大的裂紋,影響磚塊的質(zhì)量。

    隨著微型計算機技術(shù)的發(fā)展以及計算方法的不斷更新進步,計算機應用于燒結(jié)磚隧道窯冷卻帶的三維流-固耦合數(shù)值模擬成為現(xiàn)實。有關(guān)流固耦合的計算方法大致可以分為兩類:一類是解析解;一類是數(shù)值解。數(shù)值解法通常采用的是有限單元法或有限差分法。Anton is Antoniou等人計算了二維流-固耦合流場[2]。梅國棟等采用有限差分法分析三維流-固耦合問題[3]。浙江大學李迎博士建立了內(nèi)燃機活塞組-缸套-冷卻水-機體流固耦合系統(tǒng)。對柴油機穩(wěn)態(tài)工況下的耦合傳熱、周期性瞬態(tài)傳熱以及冷啟動、外特性的瞬態(tài)傳熱與流動進行了數(shù)值仿真[4]。Ed Bendel主要通過測量排氣管道表面的溫度值作為邊界條件進行數(shù)值模擬,并且利用Rad Therm和Star-CD的方法來分析對流、傳導的換熱特性[5]。劉上在FLUENT軟件平臺上,采用氣固耦合傳熱計算的方法模擬了二維擴散燃燒實驗器中的溫度、組分分布[6]。大連理工大學周禛建立了單缸的缸蓋- 缸套 - 冷卻水耦合傳熱模型并且進行了耦合傳熱仿真[7]。

    本文運用FLUENT軟件對隧道窯冷卻帶流-固耦合散熱進行數(shù)值模擬。模擬不同的風速、孔隙率、氣體流動的數(shù)學模型以及磚塊碼放方式對隧道窯散熱的影響。為磚塊的散熱找到最佳的風速、孔隙率、數(shù)學模型以及碼磚的方式。為降低磚塊的裂縫,提高磚塊的質(zhì)量尋找有效途徑。

    1 隧道窯冷卻帶的物理模型

    在GAMBIT中建立了隧道窯冷卻帶的物理模型。為了減少計算量,提高運算效率對模型進行適當?shù)暮喕?,選取六節(jié)窯車進行流-固耦合散熱的數(shù)值模擬。隧道窯冷卻帶的截面寬度為Y=3600 mm,窯車頂面到窯頂之間的距離Z=1200 mm,冷卻帶的長度為X=8000 mm。對模型進行假設:冷卻帶要求保溫性良好;窯壁可以假設為絕熱;窯尾的風機可以簡化為速度入口;出口設為壓力出口;不考慮磚塊之間的傳熱影響以及磚垛之間的輻射傳熱影響。

    圖1 隧道窯冷卻帶的物理模型Fig.1 Physical model of the tunnel kiln cooling zone

    隧道窯冷卻帶的物理模型分為磚塊區(qū)、流體區(qū)。根據(jù)數(shù)值模擬要求采用Tet/Hybrid混合網(wǎng)格劃分方法,即結(jié)構(gòu)網(wǎng)格和非結(jié)構(gòu)網(wǎng)格相結(jié)合的劃分方法。對磚塊區(qū)的網(wǎng)格劃分,格間距為4,流體區(qū)的網(wǎng)格劃分格間距為6??傆嬀W(wǎng)格總數(shù)為670155。檢查網(wǎng)格質(zhì)量確保minimum volume不為負數(shù)。采用有限體積法對控制方程進行離散化,動量方程和能量方程用一階迎風格式,迭代直至收斂。合理的劃分網(wǎng)格有利于減少計算時間節(jié)約計算成本。

    2 隧道窯冷卻帶的數(shù)學模型

    本文計算中運用的控制方程為三維常物性不可壓縮流動的連續(xù)性方程及能量方程。計算域為磚塊與空氣、粉塵流動的區(qū)域。

    2.1標準的k-ε模型

    式中,Gk是由于平均速度梯度引起的湍動能產(chǎn)生;Gb是由于浮力影響引起的湍動能產(chǎn)生;YM為可壓縮湍流脈動膨脹對總的耗散率的影響;C1ε、C2ε、C3ε為經(jīng)驗常數(shù),F(xiàn)LUENT中默認值為C1ε=1.44、C2ε=1.92、C3ε=0.09;σk=1.0、σε=1.3。

    2.2流固耦合控制方程簡化

    隧道窯冷卻帶內(nèi)氣體的流動看做粘性不可壓縮的湍流流動。它在流動與傳熱的過程中遵循三大守恒方程:質(zhì)量守恒方程、動量守恒方程(N-S方程)、能量守恒方程。

    質(zhì)量守恒方程

    最后在求解過程中采用非耦合隱式格式(即SIMPLE算法)。該算法占用的內(nèi)存小,計算快捷靈活。用有限容積法導出的離散方程可以保證具有守恒特性,而且離散方程系數(shù)的物理意義明確,是目前流動與傳熱問題的數(shù)值計算中應用最廣的一種方法。

    2.3模型邊界條件及參數(shù)設定

    在對隧道窯冷卻帶三維流固耦合散熱數(shù)值模擬的研究中,將流體的進口設置為速度入口,流體初始速度方向要與源面垂直;流體的出口設定為壓力出口[8-10];隧道窯的窯壁在整個計算的過程中設為絕熱壁;內(nèi)部窯壁與磚塊的區(qū)域設為流體區(qū)域[11],流體不可壓縮,流體中有氣相和固相,將固相中的顆粒假設為密度相同直徑均勻的顆粒[12],忽略粒子之間的相互作用。磚塊的密度為2500 kg/m3,導熱系數(shù)為80 w/m.k,熱生成率為0.57。

    3 數(shù)值模擬結(jié)果與分析

    3.1隧道窯冷卻帶流-固耦合溫度場分布

    根據(jù)FLUENT軟件數(shù)值模擬計算經(jīng)驗,流體從隧道窯窯尾以5 m/s左右的速度進入,進入窯內(nèi)流體的溫度為常溫。

    圖2是從窯長方向取得XY截面,從圖中可以看出,顏色越深代表溫度越高,氣體從窯尾進入對磚塊進行冷卻,磚塊運動的方向與氣體運動方向相反,散熱溫度隨窯車運動的方向逐漸降低。截面內(nèi)溫度分層情況比較均勻,磚塊的左右平均溫差不超過42 ℃。且在磚垛與窯壁之間或者磚垛與磚垛之間的間隙冷卻的速度較快,主要因為磚垛與窯壁或者磚垛與磚垛之間的阻力較小氣體流速較快。

    圖2 Z=0.5m截面溫度場分布圖Fig.2 Temperature feld distribution on the cross-section of Z = 0.5 m

    圖3 X=3.6m截面溫度場分布圖Fig.3 Temperature feld distribution on the cross-section of X = 3.6 m

    圖4 X=0.8m 截面溫度場分布圖Fig.4 Temperature feld distribution on the cross-section of X = 0.8 m

    圖3、圖4分別是沿著窯寬方向取得兩個YZ截面。圖5是沿著窯長方向取得XZ截面。從圖中都可以看到窯頂?shù)臏囟纫雀G車臺面的溫度低,溫度的變化梯度均勻。氣體從窯體進入磚垛時速度會加大,這是由于進入窯體的風速為5 m/s,馬赫數(shù)小于0.2可以認為氣體的流動是低速,不考慮氣體的壓縮性,氣體的速度與截面積成反比。氣體進入磚垛時截面積變小則速度增大。氣體與固相邊界(即:磚垛)相互流動時還會產(chǎn)生湍流。在整個流-固耦合散熱的過程中,磚垛附近產(chǎn)生較強的湍流影響冷卻速度,在磚垛與窯頂之間主要是層流,冷卻的速度較快,則出現(xiàn)窯頂?shù)臏囟纫雀G車臺面附近區(qū)域的溫度低。

    3.2計算結(jié)果分析

    為了研究不同風速對磚塊的平均溫度和上下溫差影響。選擇了10組不同風速進行數(shù)值模擬,計算殘差取為10-3,在計算迭代過程中,當監(jiān)測點的溫度呈現(xiàn)穩(wěn)定不變或者近似穩(wěn)定不變時,即認為收斂。對比結(jié)果如下:

    從圖6可以看出,要保證隧道窯冷卻帶磚塊冷卻的溫度與速度,同時還要保持窯內(nèi)截面的溫差最小,就要控制進入窯尾冷卻風的風速。風速過低(風速小于5 m/s)時,隨著氣體流速升高溫度逐漸降低,截面上下溫差也在逐漸降低。風速達到5 m/s時,氣體與磚垛達到產(chǎn)生湍流的條件,此時截面的上下溫差達到最低。風速繼續(xù)增加時,湍流的強度也隨著風速增加逐漸降低,此時磚垛附近區(qū)域的散熱平均溫度下降逐漸變緩,截面的上下溫差也在不斷的增大。

    從圖7可以看出,孔隙率不斷增大,磚塊碼放逐漸變得稀疏。Z=0.5 m時,XY截面的平均溫度逐漸降低。XY截面的上下溫差先減小后增大再減小的過程,主要是因為孔隙率小于0.5時,隨著孔隙率的逐漸增加流體與磚垛的相互流動時還沒有產(chǎn)生湍流,溫差還在逐漸減小。當孔隙率在0.5-0.7之間時,湍流逐漸形成并且逐漸加強,影響散熱溫差也在逐漸增大。當孔隙率大于0.7時,磚塊已經(jīng)很稀疏,湍流現(xiàn)象逐漸消失,截面的溫差也在逐漸降低。但考慮生產(chǎn)的產(chǎn)量與成本選取孔隙率為0.5較為合適。

    從圖8可以看出,二種不同氣體的流動模型冷卻曲線與隧道窯燒成曲線冷卻段曲線對比結(jié)果,兩條冷卻曲線都與隧道窯燒成曲線冷卻段曲線接近,但采用模型方式對磚塊進行冷卻曲線更接近隧道窯冷卻曲線。主要因為氣體在磚垛內(nèi)流動的方式還是湍流,同時在做數(shù)值模擬時對物理模型和數(shù)學模型都進行了適當?shù)暮喕吔鐥l件以及參數(shù)設定與實際有點差別,計算的方法也不夠完善。

    圖5 Y=1.7m 截面溫度場分布圖Fig.5 Temperature feld distribution on the cross-section of Y = 1.7 m

    圖6 Z=0.5m時不同風速對XY截面平均溫度和溫差的影響Fig.6 Effect of different draft speed on average temperature and temperature difference on XY cross-section when Z = 0.5 m

    圖7 Z=0.5m時不同孔隙率對XY截面平均溫度和溫差的影響Fig. 7 Effect of different pore rate on average temperature and temperature difference on XY cross-section when Z = 0.5 m

    圖8 不同氣體流動模型冷卻曲線與隧道窯燒成曲線冷卻段對比Fig.8 Contrast of the cooling curves of different fuid-fow models with the cooling segment of the tunnel kiln fring curve

    圖9 不同碼放方式對磚塊散熱溫度變化的影響Fig.9 Effect of different bricks arrangement on brick heat transfer

    從圖9可以看出,隨著風速加快Z=0.5 m的XY截面平均溫度下降,在相同風速、孔隙率下磚塊稀碼比磚塊密碼的平均溫度低。從圖中曲線下降的斜率可以看出:磚塊稀碼溫度隨著風速下降的速度快于磚塊密碼。主要因為磚塊稀碼時,氣體與磚垛接觸面積加大,溫度降低速度也會加快,平均溫度也會相對較低。磚塊稀碼時氣體在磚垛內(nèi)部的湍流強度小于磚塊密碼,其散熱效果也要優(yōu)于磚塊密碼。

    4 結(jié) 論

    本文運用FLUENT對隧道窯冷卻帶六節(jié)窯車段進行流-固耦合散熱數(shù)值模擬,結(jié)論如下:

    (1)用數(shù)值模擬的方法對隧道窯冷卻帶磚塊的散熱進行模擬,可以代替大量的反復實驗,節(jié)約大量的人力物力縮短研究時間。

    (2)文中分別從風速、孔隙率、碼磚方式的角度展開研究。研究表明:風速5 m/s,孔隙率0.5時磚塊散熱時上下溫差較小,平均溫度較低;在相同風速與孔隙率的條件下磚塊稀碼比磚塊密碼的散熱條件好。

    (3)選擇合適的數(shù)學計算模型可以有效的模擬隧道窯冷卻帶散熱的過程,為以后的實驗研究提供理論依據(jù)。

    (4)文中的散熱數(shù)值模擬均基于實際生產(chǎn)條件,模擬的結(jié)果與實際數(shù)據(jù)相吻合,表明文章中數(shù)值模擬的結(jié)果是可行的,接下來進一步的工作將模擬的結(jié)果應用于實際生產(chǎn)中,針對生產(chǎn)的實際情況做適當?shù)恼{(diào)整。

    參考文獻:

    [1] 百敏麗, 薛瑩瑩, 呂繼祖, 等. 隔熱材料導熱系數(shù)的數(shù)值模擬預測[J]. 熱科學與技術(shù), 2007, 6(2): 136-140.

    [2] AKYUZLU K M, ANTONIOU A. Determination of regression rate in an ablation hybrid rocket solid fuel using a physics based comprehend div mathematical model. AIAA, 2002: 2002-3577. [3] 梅國棟, 王云海. 三維流固耦合數(shù)值模擬在銅鑼山隧道安全性評價中的應用[J]. 中國安全生產(chǎn)科學技術(shù), 2009, 5(6): 57-61. [4] 李迎. 內(nèi)燃機流固耦合傳熱問題數(shù)值仿真與應用研究[D]. 浙江大學, 2006.

    [5] BENDELL E. Investigation of a coupled CFD and thermal modeling methodology for prediction of vehicle underbody temperature. SAE, 2005-01-2044.

    [6] 劉上, 孫得川. 固液發(fā)動機實驗燃燒器的氣固耦合傳熱計算[J]. 推進技術(shù), 2008, 29(3): 258-261.

    [7] 周禛, 薛冬新, 唐斌, 等. 某中速柴油機冷卻液流動及流固耦合傳熱計算分析[J]. 柴油機, 2013, 35(1): 17-32.

    [8] KAYA S, MANCHAN E. Modelling and optimization of the firing zone of a tunnel kiln to predict the optimal feed locations and mass fluxes of the fuel and secondary air. Applied Energy,2009, 86: 325-332.

    [9] NICOLAU V P, DADAM A P. Numerical and experimental thermal analysis of a tunnel kiln used in ceramic production. Journal of the Brazilian Society of Mechanical Sciences & Engineering 2009, 31 (4): 297-304.

    [10] SANTOS G M. Study of thermal behavior of a tunnel kiln used in red ceramic industry. Master Thesis, Federal University of Santa Catarina, Santa Catarina, Brazil, 2001.

    [11] NACCACHE M F, GOMES M S P. Numerical simulation of flow and heat transfer through a tunnel kiln. 18th International Congress of Mechanical Engineering, November 6-11, 2005.

    [12] OBA R, POSSAMAI T S. Thermal analysis of tunnel kiln used to produce roof tiles. Applied Thermal Engineering, 2013. 63(1): 59-65.

    通信聯(lián)系人:許寧,男,博士,教授。

    Received date:2015-12-18. Revised date: 2015-12-20.

    Correspondent author:XU Ning, male, Ph.D., Professor.

    Numerical Simulation of Fluid - Solid Coupling Heat Transfer in Tunnel Kiln Cooling Zone

    XU Ning1, 2, Wang Xuefeng1, 2, ZHANG Zhuyin2, WANG Jun2, HUANG Janjun1, XIAO Hui3, ZHI Xiao4
    (1. School of Mechanical Engineering, Yancheng Institute of Technology, Yancheng 224051, Jiangsu, China; 2. School of Mechanical Engineering, Hunan University of Technology, Zhuzhou 412007, Hunan, China; 3. Xi'an Research & Design Institute of Wall & Roof Materials, Xi'an 710061, Shaanxi, China; 4. China National Materials Group Corporation, Beijing 100036, China)

    Abstract:In order to obtain the cooling temperature feld data in the tunnel kiln cooling zone, the tunnel kiln wall bricks-fuid coupling heat transfer physical model and mathematic model were established, taking a sintered brick tunnel kiln as the research object. The tunnel kiln wall bricks-fuid coupling heat transfer simulations were performed by commercial software FLUENT, and the overall accurate intuitive temperature feld distribution was obtained. The simulation results show: the cooling effect on the tunnel kiln cooling zone is better when the draft speed is 5 m/s, the pore rate is 0.5, and the bricks arrangement is sparse, and that is also agreeable with the tunnel kiln fring curve's cooling segment. This information provides a theoretical basis for the adjustment of the draft speed, pore rate and bricks arrangement to improve the quality of bricks fring and productivity.

    Key words:tunnel kiln cooling zone; fuid - solid coupling; heat transfer; numerical simulation

    基金項目:國家科技支撐項目水泥、墻材窯爐大氣污染治理關(guān)鍵技術(shù)研究與示范(2013BCA13B01);湖南省自然科學基金資助項目(13JJ9034);江蘇省新型環(huán)保重點實驗室資助項目(AE201109)。

    收稿日期:2015-12-18。

    修訂日期:2015-12-20。

    DOI:10.13958/j.cnki.ztcg.2016.02.006

    中圖分類號:TQ174.6+53

    文獻標志碼:A

    文章編號:1006-2874(2016)02-0026-05

    猜你喜歡
    散熱數(shù)值模擬
    全密封式大功率硅整流柜的設計思考與應用
    商情(2016年43期)2016-12-26 00:00:00
    功率芯片表面絕緣層厚度對石墨烯散熱效果的影響
    IGBT模塊壓接技術(shù)研究
    LED光源散熱不良的原因及對策
    UPS儲能電源車蓄電池組散熱的設計與實現(xiàn)
    張家灣煤礦巷道無支護條件下位移的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 23:14:27
    張家灣煤礦開切眼錨桿支護參數(shù)確定的數(shù)值模擬
    科技視界(2016年18期)2016-11-03 22:57:21
    跨音速飛行中機翼水汽凝結(jié)的數(shù)值模擬研究
    科技視界(2016年18期)2016-11-03 20:38:17
    姚橋煤礦采空區(qū)CO2防滅火的數(shù)值模擬分析
    雙螺桿膨脹機的流場數(shù)值模擬研究
    科技視界(2016年22期)2016-10-18 14:53:19
    成人综合一区亚洲| 国产亚洲欧美98| 啦啦啦观看免费观看视频高清| 国产伦一二天堂av在线观看| 亚洲美女黄片视频| 国产精品日韩av在线免费观看| 午夜福利视频1000在线观看| 男女那种视频在线观看| 欧美三级亚洲精品| 国产蜜桃级精品一区二区三区| 校园人妻丝袜中文字幕| 国内精品久久久久久久电影| 免费一级毛片在线播放高清视频| 99久久久亚洲精品蜜臀av| 美女黄网站色视频| 禁无遮挡网站| 99热这里只有是精品50| 亚洲中文字幕日韩| АⅤ资源中文在线天堂| 欧美xxxx性猛交bbbb| 欧美一级a爱片免费观看看| 欧美最新免费一区二区三区| 香蕉av资源在线| 毛片一级片免费看久久久久 | 人妻制服诱惑在线中文字幕| 婷婷色综合大香蕉| 日本免费a在线| 亚洲欧美激情综合另类| 亚洲,欧美,日韩| 大型黄色视频在线免费观看| 有码 亚洲区| 日本欧美国产在线视频| 可以在线观看毛片的网站| 不卡一级毛片| 亚州av有码| a在线观看视频网站| 精品久久久久久久久久久久久| 成人国产综合亚洲| 日韩欧美一区二区三区在线观看| 日本一二三区视频观看| 国产色婷婷99| 麻豆久久精品国产亚洲av| 日日撸夜夜添| 久久精品综合一区二区三区| 欧美性感艳星| 国产私拍福利视频在线观看| 在线观看一区二区三区| 99久久中文字幕三级久久日本| 精品久久国产蜜桃| 噜噜噜噜噜久久久久久91| 麻豆国产97在线/欧美| 国产免费av片在线观看野外av| 男人舔奶头视频| 亚洲人与动物交配视频| 中出人妻视频一区二区| 国产在线男女| 婷婷精品国产亚洲av| 午夜免费成人在线视频| 香蕉av资源在线| 91久久精品国产一区二区三区| 黄色视频,在线免费观看| 国产一区二区激情短视频| 女生性感内裤真人,穿戴方法视频| 综合色av麻豆| 欧美激情在线99| 黄片wwwwww| 高清日韩中文字幕在线| 十八禁网站免费在线| 99九九线精品视频在线观看视频| 日本a在线网址| 日韩亚洲欧美综合| 亚洲四区av| 99精品久久久久人妻精品| 日韩亚洲欧美综合| 级片在线观看| 欧美一区二区精品小视频在线| 女生性感内裤真人,穿戴方法视频| 国产美女午夜福利| 又黄又爽又免费观看的视频| 国产真实伦视频高清在线观看 | 国产高清三级在线| 亚洲一区二区三区色噜噜| 最新中文字幕久久久久| 精品久久久久久久久亚洲 | 日本免费a在线| 最近中文字幕高清免费大全6 | 嫁个100分男人电影在线观看| 亚洲最大成人手机在线| 欧美最黄视频在线播放免费| 日韩人妻高清精品专区| 免费高清视频大片| 岛国在线免费视频观看| 国产91精品成人一区二区三区| 亚洲成人精品中文字幕电影| ponron亚洲| 日韩精品有码人妻一区| 我要搜黄色片| 禁无遮挡网站| 国内久久婷婷六月综合欲色啪| 国产综合懂色| 国产高清有码在线观看视频| av天堂在线播放| 国产国拍精品亚洲av在线观看| 我要搜黄色片| x7x7x7水蜜桃| av在线老鸭窝| 亚洲国产精品合色在线| 老熟妇仑乱视频hdxx| www.色视频.com| 少妇裸体淫交视频免费看高清| 亚洲最大成人手机在线| 亚洲人成网站在线播| 国产亚洲精品综合一区在线观看| 国产一级毛片七仙女欲春2| 久久精品夜夜夜夜夜久久蜜豆| 最近视频中文字幕2019在线8| 午夜精品在线福利| 天堂影院成人在线观看| 亚洲欧美日韩东京热| 国产三级中文精品| 亚洲18禁久久av| 中亚洲国语对白在线视频| 桃色一区二区三区在线观看| 欧美最黄视频在线播放免费| 亚洲色图av天堂| 午夜福利18| 女生性感内裤真人,穿戴方法视频| 亚洲精品久久国产高清桃花| 国产日本99.免费观看| 精品人妻一区二区三区麻豆 | 性插视频无遮挡在线免费观看| 一级黄片播放器| 成人一区二区视频在线观看| 精品一区二区三区视频在线观看免费| 欧美日韩精品成人综合77777| 亚洲精品国产成人久久av| 久久久色成人| 九色国产91popny在线| 欧美激情久久久久久爽电影| 日本与韩国留学比较| 精品久久国产蜜桃| 国产精品嫩草影院av在线观看 | 欧美最新免费一区二区三区| 亚洲成人中文字幕在线播放| 91久久精品国产一区二区三区| 99久久精品国产国产毛片| 直男gayav资源| 成人永久免费在线观看视频| 久久精品人妻少妇| 久久久久久国产a免费观看| 亚洲欧美日韩高清专用| 国产大屁股一区二区在线视频| 国语自产精品视频在线第100页| 黄片wwwwww| 嫩草影院精品99| 亚洲黑人精品在线| 欧美色欧美亚洲另类二区| 免费观看精品视频网站| 日韩中文字幕欧美一区二区| 男女做爰动态图高潮gif福利片| 日本免费a在线| 99热这里只有是精品在线观看| 啦啦啦观看免费观看视频高清| 亚洲中文字幕日韩| 香蕉av资源在线| 国产精品久久视频播放| 91麻豆av在线| 亚洲va日本ⅴa欧美va伊人久久| 国产三级在线视频| 我要看日韩黄色一级片| 少妇人妻精品综合一区二区 | 久久久久久久久久黄片| 美女cb高潮喷水在线观看| 别揉我奶头 嗯啊视频| 色综合色国产| 午夜福利18| 嫩草影院精品99| 日本色播在线视频| 国产91精品成人一区二区三区| 极品教师在线免费播放| 久久精品国产亚洲网站| 免费人成在线观看视频色| 中文字幕av成人在线电影| 亚洲熟妇中文字幕五十中出| 国产精华一区二区三区| 黄色丝袜av网址大全| 国产高清不卡午夜福利| 亚洲综合色惰| 一边摸一边抽搐一进一小说| 国产av一区在线观看免费| 欧美绝顶高潮抽搐喷水| 无遮挡黄片免费观看| 亚洲无线在线观看| 精品无人区乱码1区二区| 亚洲精华国产精华液的使用体验 | 91麻豆av在线| 婷婷亚洲欧美| 国产蜜桃级精品一区二区三区| 国产精品久久久久久av不卡| 日韩欧美在线二视频| 波野结衣二区三区在线| 99热网站在线观看| 国产精华一区二区三区| 国产色婷婷99| 国产精品一区二区性色av| 长腿黑丝高跟| 国产淫片久久久久久久久| a级一级毛片免费在线观看| 国内少妇人妻偷人精品xxx网站| 精品午夜福利在线看| 国产欧美日韩一区二区精品| 成人欧美大片| 色哟哟哟哟哟哟| 成人特级av手机在线观看| xxxwww97欧美| 欧美激情久久久久久爽电影| 日韩精品有码人妻一区| 精品久久久久久久末码| 国产在视频线在精品| 九九在线视频观看精品| 亚洲国产欧洲综合997久久,| 能在线免费观看的黄片| 精品国内亚洲2022精品成人| 日本免费一区二区三区高清不卡| 成年女人永久免费观看视频| 别揉我奶头~嗯~啊~动态视频| 赤兔流量卡办理| 色在线成人网| 岛国在线免费视频观看| 亚洲图色成人| 亚洲国产高清在线一区二区三| 丰满人妻一区二区三区视频av| www.www免费av| av黄色大香蕉| 国产女主播在线喷水免费视频网站 | x7x7x7水蜜桃| 日韩在线高清观看一区二区三区 | 99九九线精品视频在线观看视频| 国内精品久久久久久久电影| 日本熟妇午夜| 色噜噜av男人的天堂激情| 日本 欧美在线| 乱码一卡2卡4卡精品| 两个人的视频大全免费| 有码 亚洲区| 久9热在线精品视频| 人妻丰满熟妇av一区二区三区| 成年女人看的毛片在线观看| 亚洲自偷自拍三级| 国产淫片久久久久久久久| 免费观看精品视频网站| 日日啪夜夜撸| 搡老妇女老女人老熟妇| 天堂影院成人在线观看| 我的老师免费观看完整版| 亚洲内射少妇av| 国产中年淑女户外野战色| 欧美日本视频| eeuss影院久久| 亚洲精品影视一区二区三区av| 国产大屁股一区二区在线视频| 人妻夜夜爽99麻豆av| 久久精品综合一区二区三区| 欧美3d第一页| av福利片在线观看| 国产男人的电影天堂91| 日本撒尿小便嘘嘘汇集6| 日本熟妇午夜| 91麻豆精品激情在线观看国产| 真人一进一出gif抽搐免费| 小蜜桃在线观看免费完整版高清| 一级a爱片免费观看的视频| 久久午夜福利片| 国产又黄又爽又无遮挡在线| 欧美日韩黄片免| 18禁在线播放成人免费| 欧美另类亚洲清纯唯美| 黄色日韩在线| 免费高清视频大片| 桃色一区二区三区在线观看| 成年版毛片免费区| 美女高潮喷水抽搐中文字幕| 亚洲人与动物交配视频| 国产免费男女视频| 亚洲美女视频黄频| 真人一进一出gif抽搐免费| 日日摸夜夜添夜夜添小说| 欧美日本亚洲视频在线播放| 最近中文字幕高清免费大全6 | 亚洲成人免费电影在线观看| 国产欧美日韩精品亚洲av| 成人鲁丝片一二三区免费| 在线a可以看的网站| 色哟哟·www| 97超级碰碰碰精品色视频在线观看| 亚洲av五月六月丁香网| 国产精品一区二区三区四区久久| 狂野欧美白嫩少妇大欣赏| 人妻制服诱惑在线中文字幕| 国产私拍福利视频在线观看| 午夜激情欧美在线| 亚洲国产精品合色在线| 免费av毛片视频| 99精品久久久久人妻精品| 成人高潮视频无遮挡免费网站| 国产免费av片在线观看野外av| 能在线免费观看的黄片| 熟妇人妻久久中文字幕3abv| 99热网站在线观看| 精品欧美国产一区二区三| 亚洲欧美日韩高清在线视频| 国产三级中文精品| 99九九线精品视频在线观看视频| 麻豆av噜噜一区二区三区| 日韩,欧美,国产一区二区三区 | 欧美性猛交╳xxx乱大交人| 99久国产av精品| 全区人妻精品视频| 亚洲最大成人av| 欧洲精品卡2卡3卡4卡5卡区| 国产一级毛片七仙女欲春2| 亚洲最大成人av| 午夜福利成人在线免费观看| 人妻夜夜爽99麻豆av| 免费看美女性在线毛片视频| xxxwww97欧美| 小蜜桃在线观看免费完整版高清| 久久精品国产自在天天线| 黄色视频,在线免费观看| 久久久国产成人精品二区| 999久久久精品免费观看国产| 一级毛片久久久久久久久女| 国产av一区在线观看免费| 日韩人妻高清精品专区| 亚洲aⅴ乱码一区二区在线播放| 精品一区二区三区av网在线观看| 国产精品福利在线免费观看| 99riav亚洲国产免费| 欧美日本视频| 国产高清视频在线观看网站| 久久草成人影院| 国产一级毛片七仙女欲春2| 欧美性感艳星| 淫妇啪啪啪对白视频| 最好的美女福利视频网| 日韩欧美精品免费久久| 国产伦精品一区二区三区四那| 亚洲精品成人久久久久久| 男人和女人高潮做爰伦理| 少妇熟女aⅴ在线视频| 亚洲18禁久久av| 亚洲欧美精品综合久久99| 1000部很黄的大片| 亚洲欧美日韩卡通动漫| 日日摸夜夜添夜夜添小说| 99久久精品热视频| 麻豆精品久久久久久蜜桃| 女同久久另类99精品国产91| 国产精品自产拍在线观看55亚洲| 一本精品99久久精品77| 色在线成人网| 欧美色视频一区免费| 色综合色国产| 听说在线观看完整版免费高清| 亚洲图色成人| 免费看a级黄色片| 国产精华一区二区三区| 欧美人与善性xxx| 国产精品三级大全| 日本三级黄在线观看| 亚洲精品在线观看二区| а√天堂www在线а√下载| 日韩一区二区视频免费看| 久久精品国产亚洲av香蕉五月| 国产精品一区www在线观看 | 91久久精品国产一区二区成人| 国产精品电影一区二区三区| 免费av毛片视频| 性欧美人与动物交配| 国产爱豆传媒在线观看| 女人十人毛片免费观看3o分钟| 国产av麻豆久久久久久久| 成人国产麻豆网| 无人区码免费观看不卡| 特级一级黄色大片| 精品午夜福利在线看| 欧美日韩亚洲国产一区二区在线观看| 97碰自拍视频| 久久九九热精品免费| 免费无遮挡裸体视频| 淫秽高清视频在线观看| 又爽又黄a免费视频| 深夜a级毛片| 国产精品一区www在线观看 | 亚洲精华国产精华精| av在线亚洲专区| 99热只有精品国产| 蜜桃亚洲精品一区二区三区| 久久久久久九九精品二区国产| 国产高清有码在线观看视频| 亚洲成av人片在线播放无| 免费电影在线观看免费观看| 久久久色成人| 国产精品99久久久久久久久| 亚洲av美国av| 午夜精品久久久久久毛片777| 亚洲va日本ⅴa欧美va伊人久久| 有码 亚洲区| 老女人水多毛片| 久久精品国产自在天天线| 麻豆精品久久久久久蜜桃| 精品一区二区三区视频在线观看免费| 日本a在线网址| 国产在线男女| www日本黄色视频网| 国产美女午夜福利| 成人特级黄色片久久久久久久| 免费观看人在逋| 国产在线男女| 天堂动漫精品| 精品人妻一区二区三区麻豆 | 变态另类成人亚洲欧美熟女| 天美传媒精品一区二区| 99久久精品一区二区三区| 久久久久久伊人网av| 精品日产1卡2卡| 无遮挡黄片免费观看| 午夜激情福利司机影院| 悠悠久久av| 99riav亚洲国产免费| 欧美激情久久久久久爽电影| 我要看日韩黄色一级片| 可以在线观看毛片的网站| 一进一出好大好爽视频| 成年女人毛片免费观看观看9| 午夜福利高清视频| 免费观看的影片在线观看| 网址你懂的国产日韩在线| 日本免费a在线| 中文字幕av在线有码专区| 日本精品一区二区三区蜜桃| 亚洲美女搞黄在线观看 | 久久精品国产亚洲av香蕉五月| 亚洲精品456在线播放app | 精品无人区乱码1区二区| 1000部很黄的大片| 国产亚洲av嫩草精品影院| 婷婷亚洲欧美| 超碰av人人做人人爽久久| 可以在线观看毛片的网站| 美女黄网站色视频| 亚洲精品乱码久久久v下载方式| 国产真实伦视频高清在线观看 | 国产高清视频在线播放一区| 国产亚洲av嫩草精品影院| 亚洲av二区三区四区| 丝袜美腿在线中文| 成人高潮视频无遮挡免费网站| 日韩欧美免费精品| 亚洲性久久影院| 久久久精品欧美日韩精品| 精品午夜福利视频在线观看一区| 久久欧美精品欧美久久欧美| 午夜精品久久久久久毛片777| 国产一级毛片七仙女欲春2| 午夜免费激情av| 在线观看66精品国产| 一区二区三区四区激情视频 | 亚洲精华国产精华精| 精品福利观看| 搡老妇女老女人老熟妇| 中文资源天堂在线| 国产精品日韩av在线免费观看| 成人一区二区视频在线观看| 久久人人爽人人爽人人片va| 成人美女网站在线观看视频| 他把我摸到了高潮在线观看| 欧美高清性xxxxhd video| eeuss影院久久| 91精品国产九色| 日韩欧美精品v在线| 国产69精品久久久久777片| 自拍偷自拍亚洲精品老妇| 国产日本99.免费观看| 91午夜精品亚洲一区二区三区 | 熟女电影av网| 俺也久久电影网| 国产精品1区2区在线观看.| 少妇猛男粗大的猛烈进出视频 | 久久人人爽人人爽人人片va| 久久久久久久亚洲中文字幕| 97超视频在线观看视频| 久久国产精品人妻蜜桃| 他把我摸到了高潮在线观看| 日韩欧美在线二视频| 国产精品av视频在线免费观看| 蜜桃久久精品国产亚洲av| 精品一区二区免费观看| 高清在线国产一区| 免费av毛片视频| 老师上课跳d突然被开到最大视频| 久久久久久久久大av| 国产亚洲av嫩草精品影院| 日韩中文字幕欧美一区二区| 国产精华一区二区三区| 日本黄色视频三级网站网址| 久久久色成人| 国产乱人视频| 在线免费观看不下载黄p国产 | 中文字幕精品亚洲无线码一区| 此物有八面人人有两片| 欧美色欧美亚洲另类二区| 午夜福利18| 成年女人毛片免费观看观看9| 国产真实伦视频高清在线观看 | 久久亚洲精品不卡| 欧美国产日韩亚洲一区| 黄色女人牲交| 亚洲在线观看片| 国产一区二区在线av高清观看| 成人欧美大片| 国产成人福利小说| 亚洲av中文字字幕乱码综合| 床上黄色一级片| 99热精品在线国产| 成人特级av手机在线观看| 午夜日韩欧美国产| 免费观看人在逋| 亚洲av五月六月丁香网| 高清毛片免费观看视频网站| 国产一区二区三区视频了| 99在线视频只有这里精品首页| 欧美+亚洲+日韩+国产| 麻豆av噜噜一区二区三区| 午夜福利18| 一卡2卡三卡四卡精品乱码亚洲| 一a级毛片在线观看| 国内久久婷婷六月综合欲色啪| 美女黄网站色视频| 亚洲三级黄色毛片| 日韩强制内射视频| 国产成人一区二区在线| 色播亚洲综合网| 日韩中文字幕欧美一区二区| 熟女人妻精品中文字幕| av天堂在线播放| 女人被狂操c到高潮| 日本 av在线| 乱码一卡2卡4卡精品| 国产国拍精品亚洲av在线观看| 免费av观看视频| 免费看美女性在线毛片视频| 免费无遮挡裸体视频| 成年人黄色毛片网站| 大型黄色视频在线免费观看| 性欧美人与动物交配| 日韩精品有码人妻一区| 男女那种视频在线观看| 国产精品日韩av在线免费观看| 久久久久免费精品人妻一区二区| 亚洲国产精品成人综合色| 一级黄色大片毛片| 美女cb高潮喷水在线观看| 国产高清视频在线播放一区| 日本撒尿小便嘘嘘汇集6| 国产一区二区三区av在线 | 黄色视频,在线免费观看| 亚洲,欧美,日韩| 国产精品久久久久久av不卡| 日本-黄色视频高清免费观看| 一进一出抽搐动态| 亚洲美女黄片视频| 蜜桃久久精品国产亚洲av| 亚洲18禁久久av| 国产伦精品一区二区三区四那| 长腿黑丝高跟| 老女人水多毛片| 亚洲精品国产成人久久av| 国产精品乱码一区二三区的特点| 欧美xxxx黑人xx丫x性爽| 三级国产精品欧美在线观看| x7x7x7水蜜桃| 男女下面进入的视频免费午夜| 国产高清三级在线| 丰满乱子伦码专区| 久久久久国产精品人妻aⅴ院| 久久久久国内视频| 精品人妻偷拍中文字幕| 女人被狂操c到高潮| 精品一区二区免费观看| av专区在线播放| 免费观看人在逋| 国产老妇女一区| 国产成人av教育| 午夜福利欧美成人| 黄色欧美视频在线观看| 深夜精品福利| 在线观看舔阴道视频| 色精品久久人妻99蜜桃| 午夜福利在线观看免费完整高清在 | 999久久久精品免费观看国产| x7x7x7水蜜桃| 男女做爰动态图高潮gif福利片| 中文字幕久久专区| 听说在线观看完整版免费高清| 国产毛片a区久久久久| 成年人黄色毛片网站| 自拍偷自拍亚洲精品老妇| 中文字幕人妻熟人妻熟丝袜美| 综合色av麻豆| 欧美激情久久久久久爽电影| 好男人在线观看高清免费视频| 亚洲真实伦在线观看| 色尼玛亚洲综合影院| 变态另类丝袜制服|