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

    天然氣爐窯干燥過(guò)程多孔介質(zhì)傳熱傳質(zhì)耦合數(shù)值分析

    2017-07-05 13:22:46同濟(jì)大學(xué)機(jī)械與能源工程學(xué)院秦朝葵
    上海煤氣 2017年3期
    關(guān)鍵詞:磚坯溫度梯度傳質(zhì)

    同濟(jì)大學(xué)機(jī)械與能源工程學(xué)院 周 宇 秦朝葵

    天然氣爐窯干燥過(guò)程多孔介質(zhì)傳熱傳質(zhì)耦合數(shù)值分析

    同濟(jì)大學(xué)機(jī)械與能源工程學(xué)院 周 宇 秦朝葵

    磚坯的干燥是其塑性成型的逆過(guò)程,干燥過(guò)程涉及相變及多孔介質(zhì)傳熱傳質(zhì)的耦合問(wèn)題,文章從非穩(wěn)態(tài)傳熱傳質(zhì)和相平衡的角度出發(fā),基于多孔介質(zhì)傳熱傳質(zhì)和計(jì)算流體動(dòng)力學(xué)基礎(chǔ)理論,建立了磚坯干燥過(guò)程的非線性三維數(shù)學(xué)模型,利用流體計(jì)算軟件FLUENT對(duì)模型進(jìn)行求解,得到磚坯干燥過(guò)程中內(nèi)部含濕量和溫度的變化過(guò)程,討論了干燥介質(zhì)入口速度和溫度對(duì)干燥過(guò)程的影響,結(jié)合所得的干燥數(shù)據(jù)以及臨界干燥溫度梯度,利用Matlab對(duì)干燥介質(zhì)入口速度、溫度與干燥時(shí)間的關(guān)系進(jìn)行擬合,得出在小于本文所設(shè)定的臨界梯度的條件下,較為經(jīng)濟(jì)的干燥條件是溫度為550~573 K,風(fēng)速為0.5~2.5 m/s。

    干燥過(guò)程 多孔介質(zhì) 傳熱傳質(zhì)耦合 數(shù)學(xué)模型

    我國(guó)燒結(jié)磚瓦行業(yè)中,無(wú)論是新建設(shè)的磚廠還是改造擴(kuò)建的磚廠,隧道窯都逐漸成為了首選的磚瓦產(chǎn)品干燥以及焙燒設(shè)備。作為隧道窯所附屬的坯料干燥設(shè)施,干燥窯在整個(gè)生產(chǎn)流程中發(fā)揮著重要的作用。

    干燥過(guò)程中,隧道窯干燥室內(nèi)的坯料與干燥介質(zhì)進(jìn)行強(qiáng)烈的熱濕交換,窯內(nèi)的干燥介質(zhì)沿著爐窯長(zhǎng)度方向具有不同的溫度、速度和濕度。在不同的干燥階段,熱濕交換的強(qiáng)度和速度各不相同,隨著干燥過(guò)程的進(jìn)行,坯料和干燥介質(zhì)的各項(xiàng)參數(shù)都發(fā)生了較大的變化。干燥過(guò)程不僅消耗大量能量,而且受坯料本身物性以及外界條件的影響,干燥效果和熱效率都會(huì)發(fā)生變化,由于影響干燥過(guò)程的因素眾多,干燥中出現(xiàn)的問(wèn)題也比較復(fù)雜,如干燥倒坯、坯料出現(xiàn)的裂紋等,然而目前對(duì)于耐材生產(chǎn)熱力過(guò)程的研究主要集中于燒成階段,對(duì)于干燥階段的研究則相對(duì)較少。

    磚坯的干燥過(guò)程不僅涉及相變過(guò)程,還需要考慮多孔介質(zhì)傳熱傳質(zhì)的相互耦合,按照控制多孔介質(zhì)傳熱傳質(zhì)過(guò)程中動(dòng)量、能量和質(zhì)量傳遞的不同模式,多孔介質(zhì)的干燥過(guò)程一般可視為多孔介質(zhì)中導(dǎo)熱與質(zhì)量擴(kuò)散過(guò)程。多孔介質(zhì)傳熱傳質(zhì)研究發(fā)展到今天己經(jīng)取得了很多的進(jìn)步,研究者們?cè)谶^(guò)去的幾十年內(nèi)建立了數(shù)以百計(jì)的模型以描述多孔介質(zhì)的干燥過(guò)程,并且根據(jù)濕分在在固體物質(zhì)內(nèi)的初始狀態(tài)、驅(qū)動(dòng)力和運(yùn)動(dòng)形式等因素相繼提出液態(tài)擴(kuò)散、蒸發(fā)冷凝理論、毛細(xì)理論、Luikov理論、Whitaker體積平均理論等,由于多孔介質(zhì)本身結(jié)構(gòu)的復(fù)雜性、多變性和隨機(jī)性,使得對(duì)其傳熱傳質(zhì)過(guò)程的研究手段仍然以實(shí)驗(yàn)研究為主,但是實(shí)驗(yàn)研究相對(duì)昂貴,費(fèi)時(shí)較長(zhǎng),相對(duì)于實(shí)驗(yàn)來(lái)說(shuō),數(shù)值計(jì)算的投資要小得多,具有較好的通用性、靈活性、快速性及經(jīng)濟(jì)性。

    本文結(jié)合爐窯實(shí)測(cè)數(shù)據(jù),建立磚坯干燥數(shù)值模型,利用Fluent軟件對(duì)磚坯內(nèi)部水分遷移進(jìn)行數(shù)值計(jì)算,重點(diǎn)分析單塊磚坯在不同干燥條件下的干燥情況,分析干燥介質(zhì)入口速度、溫度對(duì)干燥效果的影響。

    1 模型及邊界條件

    1.1 計(jì)算模型

    利用ICEM繪制單塊磚坯的結(jié)構(gòu)網(wǎng)格,磚坯模型的尺寸為150×150×300 mm,如圖1所示。

    圖1 計(jì)算模型網(wǎng)格示意

    干燥介質(zhì)進(jìn)出口為單進(jìn)單出,設(shè)置左邊面為進(jìn)風(fēng)口,右邊面為出風(fēng)口,墻壁為絕熱面。為了保證計(jì)算精度,模型全部采用六面體結(jié)構(gòu)化網(wǎng)格,網(wǎng)格數(shù)量為49 062。

    1.2 邊界條件

    流場(chǎng)的數(shù)值模擬需在有限區(qū)域內(nèi)進(jìn)行,邊界條件一般是在求解區(qū)域的邊界上,求解的變量隨地點(diǎn)和時(shí)間的變化情況。對(duì)于Fluent計(jì)算,邊界條件的設(shè)置直接影響到計(jì)算結(jié)果的精度。數(shù)值模擬中,初始條件和邊界條件做如下設(shè)定和處理:

    磚坯的初溫設(shè)置為300 K,進(jìn)入干燥窯之前的含水率為6%,出窯含水率要求低于1%。濕坯的密度為2 270 kg/m3,導(dǎo)熱系數(shù)為1.3 W/(m·K),比熱容為910 J/(kg·K),為了簡(jiǎn)化模型以上參數(shù)在模擬干燥過(guò)程中視為常量。

    磚坯設(shè)置為多孔介質(zhì),設(shè)置其多孔率為常數(shù)0.4。需要注意的是,被選為多孔介質(zhì)的區(qū)域必須是Fluid類(lèi)型。

    在直角坐標(biāo)系中定義阻力系數(shù)的辦法是在二維問(wèn)題中定義一個(gè)方向矢量,或在三維問(wèn)題中定義兩個(gè)方向矢量,然后再在每個(gè)方向上定義粘性和慣性阻力系數(shù)。本次計(jì)算中粘性阻力系數(shù)為4e+10,慣性阻力系數(shù)為4e+09。

    計(jì)算過(guò)程中涉及水蒸氣和液態(tài)水共存現(xiàn)象,因此采用兩相流模型,本文對(duì)兩相流的液相、氣相及混合相均采用相應(yīng)的udf進(jìn)行定義。

    除兩相流外,本次計(jì)算開(kāi)啟能量方程,流動(dòng)采用k-epsilon模型。

    2 結(jié)果及分析

    為了分析不同干燥條件下磚坯內(nèi)溫度梯度和水分梯度的變化及其對(duì)干燥效果的影響,對(duì)不同送風(fēng)速度和送風(fēng)溫度下的干燥過(guò)程進(jìn)行了數(shù)值模擬,本節(jié)選取風(fēng)速為5 m/s、2 m/s,溫度為573 K、600 K下的干燥過(guò)程做詳細(xì)分析。

    定義水分飽和度S為坯料中含濕量占總濕分的比例,初始條件水分飽和度為1,完全干燥后為0,取干燥介質(zhì)入口速度垂直方向上的坯料的中心面(以下簡(jiǎn)稱(chēng)中心面),觀察該中心面的平均溫度以及水分飽和度的變化,以得出坯料中心面的干燥曲線。

    為了便于觀察磚坯內(nèi)部溫度、水分梯度的分布,對(duì)模型做如下處理:沿著Z方向和X方向各截取8個(gè)面,如圖2所示。

    圖2 溫度梯度與水分梯度截面示意

    沿Z方向截取的面與干燥介質(zhì)入口風(fēng)速垂直,其坐標(biāo)為-0.14~-0.1 m;沿X方向截取的面與干燥介質(zhì)入口風(fēng)速平行,其坐標(biāo)為0.1~0.14 m。溫度梯度和水分飽和度梯度的算法如下:

    式中:T1、T2——相鄰面的面平均溫度,K;

    S1、S2——相鄰面的面平均水分飽和度;

    △L——相鄰面的面的距離,m。

    2.1 入口風(fēng)速5 m/s,溫度573 K時(shí)干燥過(guò)程

    圖3所示的是干燥介質(zhì)入口風(fēng)速為5 m/s,溫度為573K時(shí)磚坯中心面溫度以及水分飽和度隨著干燥時(shí)間的變化曲線。該干燥曲線與文獻(xiàn)中提及的恒定干燥條件下典型干燥速率曲線基本吻合,間接印證了本模型的正確性。

    圖3 入口5 m/s,573 K時(shí)坯料中心面溫度以及水分飽和度曲線

    從圖中可以看出,在100 s之前屬于升溫階段,磚坯溫度升高但水分飽和度幾乎不變。在100~200 s之間,溫度迅速升高,水分也隨之快速蒸發(fā),該階段屬于加速干燥階段。在干燥時(shí)間t<200 s時(shí),磚的表面飽和度較高,此時(shí)液態(tài)水分的毛細(xì)流起主要作用,隨著干燥的進(jìn)行磚坯的表面脫水使得其表面飽和度比其內(nèi)部飽和度降低要快得多。200~500 s之間干燥平穩(wěn),屬于等速干燥階段,在干燥進(jìn)行約500 s時(shí),磚坯中心面的平均水分飽和度已降至S=0.4。

    圖4和圖5所示分別是垂直風(fēng)速方向和平行風(fēng)速方向截面的溫度以及水分飽和度在不同干燥時(shí)間的梯度。從圖4可以看出,磚坯中心的溫度梯度基本為0,離外表面越近,其溫度梯度越大。隨著干燥時(shí)間的增加,磚坯溫度和水分飽和度梯度呈逐漸減小的趨勢(shì)。從圖5可以看出,初始階段,由于邊界處存在對(duì)流,坯料兩端的溫度梯度高于中間部分,隨著時(shí)間的遞增,內(nèi)部熱傳導(dǎo)的進(jìn)行以及濕度的影響,坯料的內(nèi)部溫差變小,溫度趨于一致,并隨時(shí)間逐漸上升。由于初始時(shí)刻坯料表面的濕勢(shì)遠(yuǎn)高于環(huán)境,所以濕分是以水汽形式由表面向外界傳遞,使得邊界處的濕勢(shì)變化很快。同樣,由于質(zhì)擴(kuò)散系數(shù)較小,中間部分濕度變化相對(duì)比較緩慢,需要較長(zhǎng)時(shí)間能達(dá)到平衡。

    圖4 垂直風(fēng)速方向溫度梯度及水分梯度示意

    在磚坯的邊角處,由于同時(shí)存在兩個(gè)方向上的劇烈傳熱傳質(zhì),其溫度上升和濕度下降的速度都要較其他地方快,因此邊界處也是最易造成干燥不均勻并產(chǎn)生裂縫的地方。磚坯的不同側(cè)面之間溫度有一定的差距,但差別并不大。干燥的開(kāi)始階段側(cè)面干燥得較快,因而溫度首先上升。而正面溫度則上升較晚,但正面在表面干燥后溫度上升得較快,并在干燥的最后階段超過(guò)側(cè)面。傳質(zhì)也有類(lèi)似情況,干燥開(kāi)始階段側(cè)面響應(yīng)快,但是當(dāng)表面干了之后,速率就有明顯的下降。與之相反,正面雖然開(kāi)始干燥速率慢,但是干燥速率變化平穩(wěn)。磚坯在干燥過(guò)程中,由于水分的劇烈蒸發(fā)以及蒸發(fā)面的移動(dòng)會(huì)造成局部溫度梯度和濕度梯度過(guò)于集中,從而使得局部的水分飽和度呈上升狀態(tài)。溫度和濕度梯度過(guò)于集中均會(huì)導(dǎo)致應(yīng)力集中,而且這種集中又特別容易發(fā)生在邊角處,如果這一應(yīng)力超過(guò)磚坯的承受極限就會(huì)產(chǎn)生缺陷。

    圖5 平行風(fēng)速方向溫度梯度及水分梯度示意

    圖6 入口2 m/s,573 K時(shí)坯料中心面溫度以及水分飽和度曲線

    2.2 入口風(fēng)速2 m/s,溫度573 K時(shí)干燥過(guò)程

    圖6所示的是干燥介質(zhì)入口風(fēng)速為2 m/s,溫度為573 K時(shí)磚坯中心面的溫度以及水分飽和度隨著干燥時(shí)間的變化曲線。

    與風(fēng)速為5 m/s時(shí)相比可以看出,溫升顯然小很多,在干燥完全時(shí)的平均溫度為330 K,而5 m/s風(fēng)速下的為440 K。然而,水分的蒸發(fā)速率并沒(méi)有因?yàn)槿肟陲L(fēng)速的減小而減小,相反,與較大的入口風(fēng)速相比,2m/s風(fēng)速下干燥至S=0.1時(shí)所需要的時(shí)間為450 s,而風(fēng)速為5 m/s時(shí)需要600 s,很顯然風(fēng)速并不是決定干燥時(shí)間的唯一因素。

    圖7和圖8所示分別是垂直風(fēng)速方向和平行風(fēng)速方向截面溫度及水分飽和度在不同干燥時(shí)間的梯度。

    圖7 垂直風(fēng)速方向溫度梯度及水分梯度示意

    圖8 平行風(fēng)速方向溫度梯度及水分梯度示意

    從圖中可以看出溫度梯度的變化不是很大,并且水分梯度從一開(kāi)始就很小。

    對(duì)比5 m/s風(fēng)速下迎風(fēng)處和背風(fēng)處的水分梯度相差不大,風(fēng)速為2 m/s時(shí),Z方向背風(fēng)處的水分梯度明顯要大于迎風(fēng)處,原因是風(fēng)速較小時(shí)背風(fēng)面的回流以及擾流使得總的傳熱傳質(zhì)系數(shù)下降。與入口風(fēng)速為5 m/s時(shí)相比,2 m/s風(fēng)速下的平均溫度梯度平均減小了40%,平均水分梯度平均減少了50%,較小的溫度以及水分梯度可以有效防止熱濕應(yīng)力過(guò)大而產(chǎn)生的裂縫。

    2.3 風(fēng)速5 m/s,溫度600 K時(shí)干燥過(guò)程

    圖9所示是入口風(fēng)速為5 m/s,溫度為600 K時(shí)中心面平均溫度以及水分飽和度隨著干燥時(shí)間的變化。

    圖9 入口5 m/s,600 K時(shí)坯料中心面溫度以及水分飽和度曲線

    從圖中可以看出來(lái),風(fēng)速為5 m/s,溫度為600 K時(shí)的溫升很快,水分飽和度的降低也很快,在200 s之內(nèi)水分飽和度S就降至了0.2,從水分飽和度變化曲線可以看出來(lái),200 s之前的干燥屬于快速干燥階段。這種情況下雖然干燥速度很快,但是由于溫度梯度和水分梯度產(chǎn)生的熱濕應(yīng)力也很大,磚坯很容易開(kāi)裂。

    圖10和圖11所示分別是垂直風(fēng)速方向和平行風(fēng)速方向截面溫度以及水分飽和度在不同干燥時(shí)間的梯度。

    與風(fēng)速為5 m/s,溫度為573 K相比,該條件下垂直風(fēng)向和平行風(fēng)向的平均梯度并不是一直都大。在200 s左右是一個(gè)臨界點(diǎn),在此之前,600 K時(shí)的梯度大于573 K時(shí);200 s之后則慢慢變小,原因是600 K的條件下200 s后已經(jīng)基本干燥完畢,這種情況下很容易出現(xiàn)大的熱濕應(yīng)力,從而引起裂縫。

    圖10 垂直風(fēng)速方向溫度梯度及水分梯度示意

    圖11 垂直風(fēng)速方向溫度梯度及水分梯度示意

    3 干燥過(guò)程擬合

    在干燥的幾個(gè)主要階段中,等速干燥過(guò)程會(huì)由于熱質(zhì)傳遞的不匹配而產(chǎn)生大的熱濕應(yīng)力,從而產(chǎn)生裂紋,而降速干燥階段不會(huì)產(chǎn)生裂紋。因此存在一個(gè)干燥的臨界點(diǎn)和干燥的最大溫度梯度。根據(jù)文獻(xiàn)中提及的干燥窯內(nèi)截面水平溫度梯度要求。本節(jié)將結(jié)合前文對(duì)不同干燥條件下水分及溫度梯度的討論,對(duì)特定溫度梯度及水分梯度要求下的干燥進(jìn)行最優(yōu)化分析。本文規(guī)定溫度梯度小于5000 K/m,即50 ℃/cm,時(shí)不會(huì)產(chǎn)生裂紋。因?yàn)楦稍镞^(guò)程中有多個(gè)變量,現(xiàn)在控制一個(gè)量不變,即水分飽和度,觀察干燥至水分飽和度為0.4時(shí)的干燥時(shí)間與干燥介質(zhì)入口溫度以及風(fēng)速的關(guān)系。表1為不同風(fēng)速以及溫度下的最大梯度,灰色部分為臨界干燥條件,即在此干燥條件下其溫度梯度等于50 ℃/cm。

    表1 不同干燥條件下溫度梯度

    根據(jù)表1可以繪制相應(yīng)的二維圖,橫坐標(biāo)為干燥時(shí)間,縱坐標(biāo)為送風(fēng)速度,如圖12所示,圖中灰點(diǎn)代表其溫度梯度大于臨界值50 ℃/cm,小五星點(diǎn)表示其溫度梯度小于50℃/cm,黑點(diǎn)為臨界點(diǎn)。陰影部分的溫度梯度小于50 ℃/cm,該部分所對(duì)應(yīng)的送風(fēng)速度、送風(fēng)溫度為合適的干燥條件。

    圖12 不同干燥條件下的溫度梯度

    為得出干燥時(shí)間較小,并且溫度梯度不大于臨界點(diǎn)的送風(fēng)風(fēng)速以及溫度,利用Matlab繪制三維曲面圖如圖13所示,用來(lái)描述溫度,水分和干燥時(shí)間之間的關(guān)系,圖中所用的干燥條件都滿足其溫度梯度小于臨界點(diǎn)。

    圖13 不同干燥條件下干燥至S=0.4所需的干燥時(shí)間

    從圖13中可以看出,干燥時(shí)間最低的區(qū)域是溫度為550~573 K,風(fēng)速為0.5 ~2.5 m/s。

    4 結(jié)語(yǔ)

    本文針對(duì)多孔介質(zhì)耦合傳熱傳質(zhì),建立相關(guān)的數(shù)學(xué)模型,利用數(shù)值模擬軟件FLUENT模擬了單塊磚坯在不同邊界條件下的干燥效果。

    (1)詳細(xì)考察在送風(fēng)溫度為573K、600K,送風(fēng)風(fēng)速分別為5m/s,2m/s時(shí)磚坯干的燥曲線以及不同方向的溫度梯度和水分梯度。結(jié)果顯示,相同溫度下,隨著干燥介質(zhì)入口風(fēng)速的增大,干燥時(shí)間呈減少的趨勢(shì),但是磚坯內(nèi)部截面溫度和水分梯度變大。同樣,風(fēng)速相同時(shí),隨著溫度的升高,干燥時(shí)間減少,磚坯內(nèi)部截面溫度和水分梯度增大。

    (2)在干燥過(guò)程中由于正面和側(cè)面所處的邊界條件不同,溫度上升和水分下降的過(guò)程也不相同。在磚坯的邊角處,由于同時(shí)存在兩個(gè)方向上與外界的劇烈傳熱傳質(zhì),溫度上升和濕度下降的速度都要比其他地方快,因此邊界處也是最易造成干燥不均勻并產(chǎn)生裂縫的地方。

    (3)根據(jù)所得的干燥數(shù)據(jù)以及所確定的臨界干燥溫度梯度,利用Matlab繪制相關(guān)的三維圖用來(lái)描述送風(fēng)風(fēng)速、溫度與干燥時(shí)間的關(guān)系,并得出在小于本文所設(shè)定的臨界梯度的條件下,即溫度梯度小于50 ℃/cm時(shí),較為經(jīng)濟(jì)的干燥區(qū)域是溫度為550~573 K,風(fēng)速為0.5~2.5 m/s。

    Numerical Analysis of Coupled Heat and Mass Transfer in Porous Media within Adobe Drying Process in Industrial Furnace

    Tongji University College of Mechanical and Energy Engineering Zhou Yu Qin Chaokui

    Adobe drying is the inverse process of plastic molding and involves phase change and coupling of heat and mass transfer. From the perspective of unsteady heat and mass transfer and phase equilibrium, a three dimensional mathematical and computational model based on porous heat and mass transfer and computational fluid dynamics theories was established in this paper. By using the fluid calculation software FLUENT, the changing process of moisture content and temperature in the adobe was obtained. The influence of drying medium inlet velocity and temperature on drying process was discussed. Combining the drying datum and the critical drying temperature gradient, the relationship of drying medium inlet velocity, temperature and drying time was fitted by Matlab. Results show that the more economical drying condition is that the temperature is 550~573 K,and the wind speed is 0.5~2.5 m/s under the condition of less than the critical gradient set in this paper.

    rdrying process, porous media, coupled heat and mass transfer, mathematical model

    猜你喜歡
    磚坯溫度梯度傳質(zhì)
    溫度對(duì)中小跨徑混凝土梁橋應(yīng)力及變形的影響研究
    嚴(yán)寒地區(qū)混凝土箱梁實(shí)測(cè)溫度梯度分析
    山西建筑(2019年21期)2019-12-02 02:35:08
    溫度梯度場(chǎng)對(duì)聲表面波器件影響研究
    電子制作(2018年23期)2018-12-26 01:01:20
    磚坯
    氨基酸鹽吸收二氧化碳過(guò)程的傳質(zhì)特性
    高速鐵路CRTSⅢ型板式無(wú)砟軌道溫度梯度試驗(yàn)研究
    鐵道建筑(2015年3期)2015-12-26 09:41:44
    PTFE膜吸收CO2的工藝參數(shù)對(duì)傳質(zhì)性能的影響
    磚坯
    清潔轉(zhuǎn)向酸H+表面?zhèn)髻|(zhì)行為實(shí)驗(yàn)研究
    煤顆粒熱解的傳熱傳質(zhì)分析
    成人18禁高潮啪啪吃奶动态图| 国产精品不卡视频一区二区| 香蕉丝袜av| 亚洲精品,欧美精品| 美女国产视频在线观看| 国产精品久久久久久av不卡| 国产一级毛片在线| 国产精品不卡视频一区二区| av在线观看视频网站免费| 欧美精品国产亚洲| 不卡av一区二区三区| 黄色 视频免费看| 久久久欧美国产精品| 国产成人精品久久久久久| 汤姆久久久久久久影院中文字幕| 欧美另类一区| 又黄又粗又硬又大视频| 三上悠亚av全集在线观看| 蜜桃在线观看..| 大话2 男鬼变身卡| 国产成人91sexporn| 久久精品久久精品一区二区三区| 999精品在线视频| kizo精华| 只有这里有精品99| 亚洲国产精品一区三区| 欧美成人午夜精品| 亚洲欧美成人精品一区二区| 免费人妻精品一区二区三区视频| 成人18禁高潮啪啪吃奶动态图| 一二三四在线观看免费中文在| 国产成人精品一,二区| 欧美日韩综合久久久久久| 精品国产乱码久久久久久小说| 成人漫画全彩无遮挡| 777久久人妻少妇嫩草av网站| 欧美精品高潮呻吟av久久| 色网站视频免费| 亚洲欧洲国产日韩| 大片免费播放器 马上看| 欧美日韩亚洲高清精品| 国产探花极品一区二区| 哪个播放器可以免费观看大片| 精品第一国产精品| 啦啦啦啦在线视频资源| 在线精品无人区一区二区三| 黄片无遮挡物在线观看| 久久影院123| 少妇 在线观看| 秋霞伦理黄片| 精品福利永久在线观看| 97在线视频观看| 欧美日韩av久久| 精品一区二区三卡| 亚洲国产欧美在线一区| 99九九在线精品视频| 免费黄网站久久成人精品| 国产一区二区三区综合在线观看| 最黄视频免费看| 天美传媒精品一区二区| 天堂俺去俺来也www色官网| 18+在线观看网站| 1024视频免费在线观看| 欧美另类一区| 啦啦啦啦在线视频资源| 日韩熟女老妇一区二区性免费视频| 69精品国产乱码久久久| 男男h啪啪无遮挡| 成年av动漫网址| 岛国毛片在线播放| 日本色播在线视频| 久久韩国三级中文字幕| 国产日韩欧美亚洲二区| 永久免费av网站大全| 亚洲国产欧美在线一区| 丝袜脚勾引网站| 久久精品熟女亚洲av麻豆精品| 亚洲国产精品一区二区三区在线| 国产精品三级大全| 国产乱来视频区| 亚洲第一青青草原| 欧美人与性动交α欧美精品济南到 | 黑丝袜美女国产一区| 免费播放大片免费观看视频在线观看| 久久午夜福利片| 国产在线一区二区三区精| 午夜福利在线免费观看网站| 日日撸夜夜添| 国产老妇伦熟女老妇高清| 久久国产精品男人的天堂亚洲| 国产成人精品久久久久久| 熟女少妇亚洲综合色aaa.| 美女国产高潮福利片在线看| 18+在线观看网站| 国产精品无大码| 看免费av毛片| 男女午夜视频在线观看| 高清黄色对白视频在线免费看| av网站在线播放免费| 亚洲国产精品成人久久小说| 国产精品国产av在线观看| 国产精品成人在线| 久久国产亚洲av麻豆专区| 男人舔女人的私密视频| av.在线天堂| 熟女少妇亚洲综合色aaa.| 一二三四在线观看免费中文在| 黄片小视频在线播放| 中文欧美无线码| 精品少妇久久久久久888优播| 晚上一个人看的免费电影| 美女高潮到喷水免费观看| 日韩不卡一区二区三区视频在线| 中文字幕精品免费在线观看视频| 成年女人在线观看亚洲视频| 久久精品亚洲av国产电影网| 国产精品久久久久久av不卡| 久久精品熟女亚洲av麻豆精品| 亚洲av中文av极速乱| 成人漫画全彩无遮挡| 香蕉精品网在线| 性高湖久久久久久久久免费观看| 精品久久久久久电影网| 亚洲av免费高清在线观看| 成年动漫av网址| 麻豆av在线久日| 久久热在线av| 久久精品国产a三级三级三级| 肉色欧美久久久久久久蜜桃| 久久久精品94久久精品| 在线亚洲精品国产二区图片欧美| 人人妻人人澡人人爽人人夜夜| 久久99蜜桃精品久久| 黑人欧美特级aaaaaa片| 狠狠精品人妻久久久久久综合| 亚洲av电影在线进入| 十分钟在线观看高清视频www| 99久国产av精品国产电影| 在线观看免费视频网站a站| 亚洲精品视频女| 亚洲情色 制服丝袜| 欧美 日韩 精品 国产| 亚洲经典国产精华液单| 2018国产大陆天天弄谢| 国产精品国产三级国产专区5o| 国产探花极品一区二区| 久久久久国产一级毛片高清牌| 9热在线视频观看99| 国产成人午夜福利电影在线观看| 免费播放大片免费观看视频在线观看| 人妻少妇偷人精品九色| 欧美日本中文国产一区发布| 免费不卡的大黄色大毛片视频在线观看| 成年人午夜在线观看视频| 18+在线观看网站| 欧美日韩一级在线毛片| 我要看黄色一级片免费的| h视频一区二区三区| 婷婷色麻豆天堂久久| 久久久久久久精品精品| 国产精品麻豆人妻色哟哟久久| 男的添女的下面高潮视频| 色视频在线一区二区三区| 九九爱精品视频在线观看| 国产有黄有色有爽视频| 人人澡人人妻人| 在线亚洲精品国产二区图片欧美| 成人毛片a级毛片在线播放| 青春草亚洲视频在线观看| 久久久久久人人人人人| 亚洲欧美精品自产自拍| 中文字幕最新亚洲高清| 老女人水多毛片| 亚洲精品日韩在线中文字幕| 嫩草影院入口| 1024香蕉在线观看| 少妇人妻久久综合中文| 99热全是精品| 中文字幕人妻丝袜制服| 人人澡人人妻人| 丁香六月天网| 91国产中文字幕| 99久久综合免费| 美女中出高潮动态图| 亚洲人成电影观看| 性色avwww在线观看| 在线观看免费视频网站a站| 国产熟女午夜一区二区三区| 在线 av 中文字幕| 精品少妇黑人巨大在线播放| 久久久久国产一级毛片高清牌| 丰满饥渴人妻一区二区三| 色94色欧美一区二区| 自拍欧美九色日韩亚洲蝌蚪91| 麻豆乱淫一区二区| 在线天堂最新版资源| 91久久精品国产一区二区三区| 精品国产乱码久久久久久小说| 中文字幕精品免费在线观看视频| 久久鲁丝午夜福利片| 777米奇影视久久| 国产亚洲一区二区精品| 免费女性裸体啪啪无遮挡网站| 国产片内射在线| 日韩制服丝袜自拍偷拍| freevideosex欧美| 最近手机中文字幕大全| 亚洲精品美女久久久久99蜜臀 | 精品国产乱码久久久久久小说| 国产精品嫩草影院av在线观看| 久久精品国产亚洲av涩爱| 爱豆传媒免费全集在线观看| 久久久久久人人人人人| 亚洲精品日韩在线中文字幕| 2021少妇久久久久久久久久久| 黄片小视频在线播放| av免费在线看不卡| 中文字幕亚洲精品专区| 日韩av在线免费看完整版不卡| 丝袜美足系列| 少妇被粗大的猛进出69影院| 桃花免费在线播放| 18禁动态无遮挡网站| 人成视频在线观看免费观看| 色播在线永久视频| 亚洲,欧美精品.| 搡女人真爽免费视频火全软件| av女优亚洲男人天堂| 女人久久www免费人成看片| www.精华液| www.自偷自拍.com| 青草久久国产| 五月天丁香电影| 国产一区二区三区综合在线观看| 亚洲精品国产色婷婷电影| 韩国精品一区二区三区| 日本色播在线视频| 国产乱来视频区| 日本-黄色视频高清免费观看| 国产综合精华液| 久久久a久久爽久久v久久| 如日韩欧美国产精品一区二区三区| 亚洲综合精品二区| 黄色配什么色好看| 久久ye,这里只有精品| 精品亚洲乱码少妇综合久久| 亚洲av男天堂| 精品国产国语对白av| 久久婷婷青草| 亚洲国产av新网站| 久久精品人人爽人人爽视色| 亚洲av电影在线观看一区二区三区| 欧美精品高潮呻吟av久久| 美女主播在线视频| 国产乱来视频区| 久久精品国产亚洲av高清一级| 黄色毛片三级朝国网站| 欧美黄色片欧美黄色片| 国产片内射在线| 久久久a久久爽久久v久久| 日本午夜av视频| 97人妻天天添夜夜摸| 99国产综合亚洲精品| 国产综合精华液| 午夜福利乱码中文字幕| 中文精品一卡2卡3卡4更新| 九九爱精品视频在线观看| 亚洲人成电影观看| 午夜激情av网站| 亚洲av电影在线进入| 在线观看美女被高潮喷水网站| 美女中出高潮动态图| 菩萨蛮人人尽说江南好唐韦庄| 亚洲三区欧美一区| 午夜日本视频在线| 中文乱码字字幕精品一区二区三区| 亚洲中文av在线| 91久久精品国产一区二区三区| 天天影视国产精品| 亚洲av男天堂| 日日啪夜夜爽| 精品一区二区免费观看| 黑人巨大精品欧美一区二区蜜桃| 69精品国产乱码久久久| 妹子高潮喷水视频| 2022亚洲国产成人精品| 国产精品熟女久久久久浪| 久久女婷五月综合色啪小说| 久久青草综合色| 亚洲欧美成人精品一区二区| 亚洲欧洲精品一区二区精品久久久 | videossex国产| 成年人免费黄色播放视频| 国产精品欧美亚洲77777| 日本91视频免费播放| 成人国产麻豆网| 免费看av在线观看网站| 黄片无遮挡物在线观看| 国产精品国产三级国产专区5o| 人成视频在线观看免费观看| 亚洲国产色片| 国产免费福利视频在线观看| 亚洲综合精品二区| 精品国产乱码久久久久久男人| 少妇猛男粗大的猛烈进出视频| 中文字幕人妻丝袜一区二区 | 久久人人爽人人片av| 狠狠婷婷综合久久久久久88av| 亚洲成人手机| 久久久久国产网址| 亚洲精品久久午夜乱码| 精品人妻熟女毛片av久久网站| 搡老乐熟女国产| 男人添女人高潮全过程视频| 少妇人妻久久综合中文| 久久精品国产a三级三级三级| 欧美日韩一级在线毛片| 午夜精品国产一区二区电影| 美女主播在线视频| 国产欧美日韩一区二区三区在线| 69精品国产乱码久久久| 亚洲国产成人一精品久久久| 欧美人与性动交α欧美精品济南到 | 亚洲欧美成人精品一区二区| 有码 亚洲区| 久久国产精品男人的天堂亚洲| 99久久中文字幕三级久久日本| 中文字幕制服av| 中文欧美无线码| 国产乱来视频区| 视频在线观看一区二区三区| 麻豆乱淫一区二区| 街头女战士在线观看网站| 欧美精品亚洲一区二区| 欧美精品人与动牲交sv欧美| 久久久久精品性色| 满18在线观看网站| 久久久久网色| 纵有疾风起免费观看全集完整版| 9热在线视频观看99| 国产成人精品无人区| 日韩免费高清中文字幕av| 大香蕉久久网| 黄色一级大片看看| 日日爽夜夜爽网站| 人人妻人人澡人人爽人人夜夜| 黄色毛片三级朝国网站| 人人妻人人澡人人爽人人夜夜| 卡戴珊不雅视频在线播放| 国产探花极品一区二区| 高清不卡的av网站| 欧美最新免费一区二区三区| 国产在视频线精品| 免费观看av网站的网址| 欧美 亚洲 国产 日韩一| 肉色欧美久久久久久久蜜桃| www日本在线高清视频| 中文字幕另类日韩欧美亚洲嫩草| 欧美国产精品va在线观看不卡| av有码第一页| 毛片一级片免费看久久久久| www.av在线官网国产| 2022亚洲国产成人精品| 麻豆乱淫一区二区| 国产精品无大码| 亚洲av综合色区一区| 91国产中文字幕| 久久精品久久久久久久性| 成人毛片a级毛片在线播放| 啦啦啦在线免费观看视频4| 黄片小视频在线播放| 91久久精品国产一区二区三区| 久久精品国产a三级三级三级| 免费高清在线观看视频在线观看| 蜜桃国产av成人99| 免费观看a级毛片全部| 国产精品国产三级国产专区5o| 只有这里有精品99| 欧美成人精品欧美一级黄| 中文字幕人妻熟女乱码| 久久人人97超碰香蕉20202| 久久精品国产自在天天线| 欧美成人精品欧美一级黄| 久久精品国产自在天天线| 狠狠精品人妻久久久久久综合| 26uuu在线亚洲综合色| 免费黄网站久久成人精品| 精品国产一区二区久久| 男女午夜视频在线观看| 又粗又硬又长又爽又黄的视频| 久久狼人影院| 久久久精品94久久精品| 国产精品二区激情视频| 亚洲 欧美一区二区三区| 女人久久www免费人成看片| 嫩草影院入口| 欧美xxⅹ黑人| 最近的中文字幕免费完整| freevideosex欧美| 黄色 视频免费看| 欧美bdsm另类| 国产在线一区二区三区精| 国产成人免费观看mmmm| 男女无遮挡免费网站观看| 欧美少妇被猛烈插入视频| 欧美变态另类bdsm刘玥| 成人漫画全彩无遮挡| 亚洲在久久综合| 观看美女的网站| 国产在线免费精品| 国产在视频线精品| 满18在线观看网站| 麻豆乱淫一区二区| 久久精品国产自在天天线| 亚洲图色成人| 自拍欧美九色日韩亚洲蝌蚪91| 午夜免费观看性视频| av一本久久久久| 欧美老熟妇乱子伦牲交| 日日摸夜夜添夜夜爱| 久久久国产精品麻豆| 好男人视频免费观看在线| 中文字幕人妻丝袜制服| 狂野欧美激情性bbbbbb| 亚洲精品,欧美精品| 欧美成人精品欧美一级黄| 日本wwww免费看| 新久久久久国产一级毛片| 国产男人的电影天堂91| 亚洲欧美成人精品一区二区| 免费观看无遮挡的男女| 精品国产一区二区三区久久久樱花| 欧美 亚洲 国产 日韩一| 一级片免费观看大全| 国产精品一二三区在线看| 欧美日韩一级在线毛片| 一个人免费看片子| 激情五月婷婷亚洲| 久久久久国产网址| 高清欧美精品videossex| 老鸭窝网址在线观看| 一边摸一边做爽爽视频免费| 日本午夜av视频| 国产乱人偷精品视频| 欧美精品国产亚洲| h视频一区二区三区| 亚洲av电影在线观看一区二区三区| 久久久久精品人妻al黑| 咕卡用的链子| 国产福利在线免费观看视频| 免费黄色在线免费观看| 丝袜人妻中文字幕| 黄色 视频免费看| 日本vs欧美在线观看视频| 日韩视频在线欧美| 亚洲av电影在线观看一区二区三区| 成年人免费黄色播放视频| 日韩欧美精品免费久久| 老鸭窝网址在线观看| 久久久久久免费高清国产稀缺| 欧美日韩一级在线毛片| 亚洲国产毛片av蜜桃av| 中文字幕最新亚洲高清| 女人久久www免费人成看片| 高清欧美精品videossex| 91午夜精品亚洲一区二区三区| 性色av一级| 99热国产这里只有精品6| av在线app专区| 色哟哟·www| 中文字幕人妻丝袜一区二区 | 老熟女久久久| 2018国产大陆天天弄谢| 中文字幕人妻丝袜一区二区 | 老女人水多毛片| 欧美国产精品一级二级三级| 看免费成人av毛片| 不卡av一区二区三区| 亚洲精品久久久久久婷婷小说| 国产精品久久久久久久久免| 各种免费的搞黄视频| 80岁老熟妇乱子伦牲交| 中文字幕av电影在线播放| 熟女电影av网| 国产黄色视频一区二区在线观看| 国产 一区精品| 国产亚洲精品第一综合不卡| 91久久精品国产一区二区三区| www.av在线官网国产| 免费在线观看视频国产中文字幕亚洲 | av国产久精品久网站免费入址| 大话2 男鬼变身卡| 久久久久久免费高清国产稀缺| 精品久久久久久电影网| 人妻 亚洲 视频| av在线播放精品| 午夜福利在线观看免费完整高清在| 欧美成人午夜精品| 亚洲伊人久久精品综合| 午夜福利一区二区在线看| 精品少妇久久久久久888优播| 熟女电影av网| 国产精品一二三区在线看| 色婷婷av一区二区三区视频| 久热这里只有精品99| 日日爽夜夜爽网站| 人成视频在线观看免费观看| 18禁动态无遮挡网站| 国产一区亚洲一区在线观看| 久久亚洲国产成人精品v| 国产成人a∨麻豆精品| 大话2 男鬼变身卡| 丝袜在线中文字幕| 中文精品一卡2卡3卡4更新| 中文字幕人妻丝袜一区二区 | 精品国产一区二区三区四区第35| 精品少妇久久久久久888优播| 免费女性裸体啪啪无遮挡网站| 国产亚洲一区二区精品| 一本色道久久久久久精品综合| 另类亚洲欧美激情| 久久午夜福利片| 亚洲精品视频女| 韩国精品一区二区三区| 丰满乱子伦码专区| 国产精品.久久久| 欧美日韩精品网址| 少妇人妻 视频| 成年女人在线观看亚洲视频| 国产成人免费无遮挡视频| av网站免费在线观看视频| 国产精品久久久久久久久免| 大片电影免费在线观看免费| 老汉色av国产亚洲站长工具| 成年动漫av网址| 久久97久久精品| 天天操日日干夜夜撸| 国产成人欧美| 高清av免费在线| 岛国毛片在线播放| 水蜜桃什么品种好| 极品人妻少妇av视频| kizo精华| 亚洲欧美清纯卡通| 成人手机av| 日本91视频免费播放| 国产综合精华液| 亚洲av日韩在线播放| 亚洲精品中文字幕在线视频| 亚洲av福利一区| 久久久久精品人妻al黑| 曰老女人黄片| 亚洲国产av新网站| 久久亚洲国产成人精品v| 婷婷色综合www| 男人操女人黄网站| 久久精品熟女亚洲av麻豆精品| 国产免费视频播放在线视频| 免费观看在线日韩| 亚洲成人一二三区av| 亚洲综合色网址| 女性生殖器流出的白浆| 一级黄片播放器| 久久久久久伊人网av| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 黄片小视频在线播放| 我要看黄色一级片免费的| 哪个播放器可以免费观看大片| 日韩精品免费视频一区二区三区| 黄色毛片三级朝国网站| 亚洲精品日韩在线中文字幕| 一个人免费看片子| 久久精品久久精品一区二区三区| 欧美日韩一区二区视频在线观看视频在线| 精品国产超薄肉色丝袜足j| 一本—道久久a久久精品蜜桃钙片| 久久av网站| 春色校园在线视频观看| 日韩中字成人| 国产午夜精品一二区理论片| 精品国产国语对白av| 中文欧美无线码| 日韩电影二区| 国产免费一区二区三区四区乱码| 黄片无遮挡物在线观看| 日韩电影二区| av国产久精品久网站免费入址| 欧美精品一区二区免费开放| 日本欧美视频一区| av卡一久久| 黄片无遮挡物在线观看| 精品视频人人做人人爽| 男男h啪啪无遮挡| 亚洲精品,欧美精品| 美女午夜性视频免费| 亚洲婷婷狠狠爱综合网| 国产精品无大码| 秋霞伦理黄片| av电影中文网址| 欧美日韩视频精品一区| 两个人免费观看高清视频| 亚洲国产欧美网| 欧美日韩视频精品一区| 80岁老熟妇乱子伦牲交| 我要看黄色一级片免费的| 久久久久精品久久久久真实原创| 久久亚洲国产成人精品v| 中文乱码字字幕精品一区二区三区| 亚洲国产av新网站| 99国产精品免费福利视频| 国产老妇伦熟女老妇高清| 精品一区二区免费观看| 捣出白浆h1v1| 一边摸一边做爽爽视频免费| 久久久国产一区二区|