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

    平直翅片管式換熱器結(jié)構(gòu)參數(shù)的優(yōu)化

    2019-01-03 02:15:46費(fèi)繼友田士博王楓李花
    關(guān)鍵詞:翅片管徑換熱器

    費(fèi)繼友, 田士博, 王楓, 李花

    (大連交通大學(xué) 機(jī)車(chē)車(chē)輛工程學(xué)院,遼寧 大連 116028)*

    0 引言

    平直翅片式換熱器作為一種能源工業(yè)的重要設(shè)備,以其結(jié)構(gòu)簡(jiǎn)單、工藝成熟、互換性高等特點(diǎn)得到廣泛應(yīng)用.盡管此翅片結(jié)構(gòu)可以降低空氣流阻,增加傳熱面積進(jìn)而強(qiáng)化換熱效果,但往往存在著擾流強(qiáng)度弱,換熱效率低,整機(jī)體積大的缺點(diǎn)[1].因此對(duì)平直翅片式換熱器進(jìn)行結(jié)構(gòu)優(yōu)化研究,有助于提高換熱效率及增強(qiáng)結(jié)構(gòu)緊湊性.

    目前,平直式換熱器的優(yōu)化多以單目標(biāo)函數(shù)進(jìn)行研究[2-5],例如:換熱因子j、阻力因子f、換熱量等,且對(duì)換熱器體積的研究較少.康海軍[2]、陰繼翔等[3]以阻力因子f,PEC(綜合性能評(píng)價(jià)指標(biāo))分別為單目標(biāo)函數(shù)對(duì)平直翅片結(jié)構(gòu)進(jìn)行優(yōu)化,僅提高了換熱器效率.馬士偉[4]、葉立平等[5]以換熱量為單目標(biāo)函數(shù)對(duì)平直式翅片結(jié)構(gòu)進(jìn)行分析,僅提高了其換熱量.然而對(duì)于換熱器而言,評(píng)價(jià)其換熱效率的標(biāo)準(zhǔn)絕非單一目標(biāo)即可,由于這些目標(biāo)互相耦合,互相關(guān)聯(lián),所以需要多目標(biāo)同時(shí)參與并行分析,往往要進(jìn)一步對(duì)翅片結(jié)構(gòu)進(jìn)行多目標(biāo)優(yōu)化設(shè)計(jì).而遺傳算法可以較完美的解決多目標(biāo)優(yōu)化問(wèn)題,楊輝著等[6]利用遺傳算法完成了板式換熱器的多目標(biāo)結(jié)構(gòu)優(yōu)化,徐東等[7]利用遺傳算法完成鋸齒型翅片的多目標(biāo)結(jié)構(gòu)優(yōu)化.另外,由于換熱器在制冷設(shè)備中所占體積較大,其體積的優(yōu)化對(duì)整機(jī)的緊湊性有著重要影響,在平直式翅片換熱器的優(yōu)化設(shè)計(jì)中,還缺乏多目標(biāo)優(yōu)化下對(duì)體積的研究.

    因此,本文以平直翅片管換熱器為研究對(duì)象,在維持入口雷諾數(shù)恒定時(shí),基于遺傳算法以最大換熱因子j,最小流動(dòng)阻力f,最小換熱器體積Vhe為目標(biāo)函數(shù)對(duì)翅片間距,管徑,橫向中心距,縱向中心距進(jìn)行優(yōu)化設(shè)計(jì).

    1 換熱器數(shù)值模擬方法及計(jì)算模型

    本文計(jì)算模型中的對(duì)流換熱過(guò)程遵循質(zhì)量守恒、動(dòng)量守恒、能量守恒,并采用有限體積法對(duì)其進(jìn)行離散,模擬計(jì)算采用SIMPLE算法,動(dòng)量及能量方程使用二階迎風(fēng)格式[8].

    1.1 控制方程及湍流模型

    本文所用方程包含質(zhì)量、動(dòng)量、能量控制方程如下:

    質(zhì)量方程為

    (1)

    動(dòng)量方程為

    (2)

    能量方程為

    (3)

    流體在翅片間圓柱繞流的情況十分復(fù)雜,為模擬這一情況,工程上常采用半經(jīng)驗(yàn)的標(biāo)準(zhǔn)k-ε方程[8]如下:

    (4)

    (5)

    1.2 目標(biāo)因子定義

    (1)阻力因子

    (6)

    式中,ΔP為進(jìn)出口空氣壓降,Pa,ρ為流道中空氣的密度,kg/m3,μmax為流道最小截面處空氣流速,m/s,De為圓管當(dāng)量直徑,m,L為翅片沿空氣流通方向的長(zhǎng)度,m.

    (2)換熱因子

    (7)

    式中,Nu努塞爾系數(shù),Re為雷諾系數(shù),Pr為普朗特系數(shù);

    (3)換熱器體積Vhe,本模型中所計(jì)算的體積為圖1中翅片,管壁以及計(jì)算域體積,不包含進(jìn)出口的延長(zhǎng)部分,單位為m3;

    1.3 翅片管幾何模型

    換熱器翅片結(jié)構(gòu)數(shù)量龐大,為減少CFD模擬計(jì)算成本,需將模型做必要的簡(jiǎn)化處理.如圖1,

    (a)

    (b)

    針對(duì)目前平直式翅片換熱器常用結(jié)構(gòu)參數(shù)范圍,以H取1.4~2.6 mm,R取5~9 mm,S1取16~21 mm,S2取8~13 mm為優(yōu)化范圍.

    1.4 網(wǎng)格劃分及邊界條件

    利用ICEM模塊生成精度更高的參數(shù)結(jié)構(gòu)化網(wǎng)格可以很好的與Fluent進(jìn)行實(shí)時(shí)關(guān)聯(lián),同時(shí)為保證捕捉避免流動(dòng)特性,對(duì)近壁面網(wǎng)格進(jìn)行局部加密處理,如圖2所示.

    圖2 網(wǎng)格邊界層加密示意圖

    為保證模型計(jì)算精度,對(duì)模型進(jìn)行網(wǎng)格無(wú)關(guān)性驗(yàn)證,圖3為j和f隨網(wǎng)格數(shù)量從17 320加密至661680變化的曲線,可以觀察到網(wǎng)格數(shù)量在320000時(shí),變化趨于平緩,考慮計(jì)算的精度以及計(jì)算的成本和周期,選擇網(wǎng)格數(shù)量為320 000進(jìn)行計(jì)算.

    圖3 網(wǎng)格無(wú)關(guān)性驗(yàn)證

    模擬計(jì)算時(shí),邊界條件設(shè)置如下:

    (1)入口:流體域空氣進(jìn)口為速度入口,空氣流速均勻分布;

    (2)出口:流動(dòng)狀態(tài)未知,故為自由出流;

    (3)翅片:翅片與流體接觸表面采用自身導(dǎo)熱和表面對(duì)流換熱的耦合方式;

    (4)熱管:忽略圓管的厚度,認(rèn)為管內(nèi)外壁面的溫度相同,設(shè)定熱管的壁面為恒溫壁面,溫度為316 K,其余幾何體模型表面均設(shè)置為對(duì)稱性邊界條件;

    (5)翅片材質(zhì)為鋁,管壁材質(zhì)為銅,且忽略外界輻射換熱及由于銅管與翅片接觸不良產(chǎn)生的熱損失,通道流體為空氣,視為常物性.

    1.5 數(shù)值模擬的正確性驗(yàn)證

    仿真算法驗(yàn)證所用結(jié)構(gòu)參數(shù)如下:H為2.00mm;F為0.15 mm;R為7.00mm;S1為18.20 mm;S210.50 mm.

    為保證仿真計(jì)算結(jié)果準(zhǔn)確性,以上述結(jié)構(gòu)參數(shù),改變?nèi)肟陲L(fēng)速(2.5~6.5 m/s),得到衡量空氣阻力的系數(shù)f,以及衡量換熱的系數(shù)j.將仿真計(jì)算結(jié)果與文獻(xiàn)[10]中實(shí)驗(yàn)值進(jìn)行對(duì)比,如圖4可知,f因子與實(shí)驗(yàn)對(duì)比,均方根誤差為8.07%,j因子與實(shí)驗(yàn)對(duì)比,均方根誤差為1.66%,吻合度較好,雖然在低雷諾數(shù)最大誤差達(dá)到20%,但實(shí)際工況雷諾數(shù)要遠(yuǎn)大于此,因此可以認(rèn)為數(shù)值算法切實(shí)可靠.

    圖4 仿真模擬與實(shí)驗(yàn)數(shù)值結(jié)果對(duì)比

    2 換熱器數(shù)值模擬結(jié)果分析

    2.1 翅片間距H的影響

    如圖5,R=7 mm,S1=18.5 mm,S2=10.5mm時(shí),不同翅片間距H(1.4~2.6 mm)對(duì)f因子、j因子、體積Vhe的影響.如圖5可知,隨著翅片間距H的增加,阻力系數(shù)f下降,換熱因子j上升,體積Vhe增加.隨著H的增加,導(dǎo)致空氣流道的擴(kuò)大,流速下降必然使流動(dòng)阻力下降,進(jìn)而促進(jìn)湍流充分發(fā)展,帶走更多熱量最終導(dǎo)致?lián)Q熱因子f的增加,但由于間距的增加,為保持相同的換熱面積,必然使整機(jī)體積增加.增加H有利于強(qiáng)化換熱,但考慮到整機(jī)體積,應(yīng)在允許的空間條件下適當(dāng)增加翅片間距.

    圖5 翅片間距 H對(duì)換熱及體積的影響

    2.2 管徑R的影響

    如圖6,H=2.0 mm,S1=18.5 mm,S2=10.5 mm時(shí),不同管徑R(5~9 mm)對(duì)f因子、j因子、體積Vhe的影響.如圖6可知,隨著R增加,阻力系數(shù)f增大,換熱因子j減小,體積Vhe是不變的.隨著管徑R增加,流道中管路阻擋面積擴(kuò)大,進(jìn)而使流動(dòng)受阻加劇,壓降升高,使阻力系數(shù)f快速增加,由于流動(dòng)受阻嚴(yán)重,湍流得不到有效的擴(kuò)展,熱量不能及時(shí)的隨主流流出,最終導(dǎo)致?lián)Q熱因子的下降.由于管徑的改變不影響相鄰管路的中心距離,翅片所在的整體空間不發(fā)生尺寸變化,所以體積是不變的.不論從換熱效率還是工程成本考慮,小管徑已成為必要的發(fā)展趨勢(shì),但考慮管內(nèi)制冷劑的流動(dòng)阻力,應(yīng)結(jié)合使用工況下選取適當(dāng)?shù)墓軓?

    圖6 管徑 R對(duì)換熱及體積的影響

    2.3 管橫向間距S1的影響

    如圖7,H=2.0 mm,R=7.0 mm,S2=10.5mm時(shí),不同橫向中心距S1(16~21 mm)對(duì)f因子、j因子、體積Vhe的影響.隨著S1的增加,換熱因子j下降,阻力系數(shù)不斷上升后趨于平緩,體積Vhe線性增加.隨著S1增加,管路之間橫向距離的增加,降低了最小截面流速,使湍流擾動(dòng)程度下降,且流體流動(dòng)長(zhǎng)度增加使其流動(dòng)阻力增加,長(zhǎng)度的增加導(dǎo)致熱量不能及時(shí)的導(dǎo)出,進(jìn)而使換熱因子j呈下降趨勢(shì),橫向間距的增加必然導(dǎo)致翅片寬度的增加,導(dǎo)致其體積線性增長(zhǎng).橫向間距對(duì)f,j,Vhe均有顯著的影響,不論從換熱效率還是從體積上分析,可見(jiàn)S1應(yīng)越小越好.

    圖7 橫向間距 S1對(duì)換熱及體積的影響

    2.4 管縱向間距S2的影響

    如圖8,H=2.0 mm,R=7.0 mm,S1=18.5mm時(shí),不同縱向中心距S2(8~13 mm)對(duì)f因子、j因子、體積Vhe的影響.隨著縱向間距S2的增加,換熱因子j,阻力系數(shù)f均下降,體積Vhe線性增加.隨著S2的增加,迎風(fēng)流道面積增加,空氣流通速度下降,進(jìn)而流動(dòng)阻力下降,由于翅片間距不變,而S2的增加相當(dāng)于翅片高度,進(jìn)而導(dǎo)致湍流進(jìn)一步發(fā)展惡化,并不能起到促進(jìn)換熱的作用,所以換熱因子j呈下降趨勢(shì),同時(shí)S2的增加,必然導(dǎo)致計(jì)算域體積擴(kuò)大,所以導(dǎo)致體積Vhe的增加.阻力因子f下降幅度相對(duì)較大,可以認(rèn)為S2增加對(duì)換熱效率有利,但要考慮其帶來(lái)的體積增長(zhǎng),所以應(yīng)該適當(dāng)選取合適的S2.

    圖8 縱向間距 S2對(duì)換熱及體積的影響

    3 多目標(biāo)優(yōu)化結(jié)果分析

    3.1 多目標(biāo)優(yōu)化原理

    本文所用多目標(biāo)遺傳算法(MOGA)的優(yōu)化過(guò)程如圖9.對(duì)于換熱器結(jié)構(gòu)的優(yōu)化,不僅要求具有充分的換熱效果,較低的空氣流動(dòng)阻力及更小的體積結(jié)構(gòu).為兼顧三者,首先在試驗(yàn)設(shè)計(jì)中使用精度高,預(yù)測(cè)性好的中心組合設(shè)計(jì),是響應(yīng)曲面中最常用的二階實(shí)驗(yàn)設(shè)計(jì)[11].然后根據(jù)試驗(yàn)點(diǎn)建立響應(yīng)面方法構(gòu)建目標(biāo)函數(shù),響應(yīng)面方法是一種通過(guò)已知點(diǎn)來(lái)預(yù)測(cè)未知或需要觀察點(diǎn)的插值方法,其利用方差的變化來(lái)表達(dá)空間的變化,可以保證由空間分布得到的預(yù)測(cè)值的誤差最小,使結(jié)果良好的全局性及預(yù)測(cè)性.最后利用遺傳算法對(duì)目標(biāo)函數(shù)進(jìn)行迭代求解并利用Pareto機(jī)制對(duì)解集進(jìn)行最優(yōu)選取[12].

    圖9 遺傳算法優(yōu)化流程

    在多目標(biāo)優(yōu)化中,各個(gè)目標(biāo)函數(shù)往往是互相矛盾的,不存在一組解使所有目標(biāo)函數(shù)達(dá)到最優(yōu)值,經(jīng)多目標(biāo)遺傳算法求近似解的分析方法,克服了以往單目標(biāo)分析的不足,在對(duì)換熱器結(jié)構(gòu)優(yōu)化的處理上非常有效[13].

    3.2 多目標(biāo)優(yōu)化處理

    基于Workbench的目標(biāo)驅(qū)動(dòng)優(yōu)化系統(tǒng),對(duì)模型進(jìn)行參數(shù)化建模,并使全局網(wǎng)格參數(shù)化,設(shè)定換熱因子j,阻力系數(shù)f,換熱器體積Vhe為目標(biāo)函數(shù),維持雷諾數(shù)為2219時(shí),根據(jù)各結(jié)構(gòu)常見(jiàn)的變化范圍,給予結(jié)構(gòu)參數(shù)優(yōu)化尺寸如表1.在優(yōu)化模塊中,選取精度高的中心組合設(shè)計(jì)并生成正交試驗(yàn)表格.優(yōu)化算法選用多目標(biāo)遺傳算法,設(shè)置響應(yīng)面生成樣本點(diǎn)1000,迭代樣本200,允許max pareto為70%.換熱因子j最大化,阻力因子f最小化,翅片體積Vhe最小化.

    表1 結(jié)構(gòu)參數(shù)變化范圍 mm

    3.3 多目標(biāo)優(yōu)化結(jié)果分析

    表2為以最大換熱因子j,最小阻力系數(shù)f,和最小體積Vhe為目標(biāo)函數(shù)對(duì)翅片間距、管徑、熱管橫向中心距、熱管縱向中心距的優(yōu)化結(jié)果.維持雷諾數(shù)為2 219時(shí),有必要考慮實(shí)際生產(chǎn)當(dāng)H=2.6mm,R=7.00 mm,S1=17.20 mm,S2=8.50mm時(shí),阻力系數(shù)f減少26.30%,換熱因子j提高1.11%,體積Vhe減少13.93%.優(yōu)化結(jié)果為減小體積而犧牲了基管的衡縱間距,進(jìn)而使換熱面積減小,最終導(dǎo)致?lián)Q熱因子下降.為彌補(bǔ)j因子的下降,采用維持小管徑而略增大翅片間距的措施來(lái)保證原有的換熱效率,較大的間距不僅補(bǔ)足了原有換熱因子的下降,而且降低了阻力因子f,進(jìn)而降低設(shè)備運(yùn)行成本.綜合認(rèn)為換熱器在維持原有的換熱效率時(shí),體積及空氣流動(dòng)阻力均有下降,由此可見(jiàn),利用多目標(biāo)遺傳算法,可以有效的對(duì)平直式翅片換熱器進(jìn)行優(yōu)化設(shè)計(jì).

    4 結(jié)論

    利用Workbench目標(biāo)驅(qū)動(dòng)優(yōu)化系統(tǒng),以最小f因子,最大j因子,最小Vhe為目標(biāo)函數(shù),應(yīng)用多目標(biāo)遺傳算法對(duì)平直翅片式換熱器的翅片結(jié)構(gòu)參數(shù)進(jìn)行優(yōu)化計(jì)算,得出以下結(jié)論:

    (1)維持入口雷諾數(shù)為2 219時(shí),翅片間距H=2.60 mm、管徑R=7.00 mm、 管橫向間距S1=17.20 mm,管縱向間距S2=8.50 mm時(shí)性能較佳,其阻力系數(shù)下降26.30%,換熱因子提高1.11%,換熱器體積減少13.93%;

    (2)基于Workbench多目標(biāo)構(gòu)優(yōu)化技術(shù),可以有效的克服傳統(tǒng)優(yōu)化方法不全面的缺點(diǎn),以及對(duì)實(shí)驗(yàn)關(guān)聯(lián)式的依賴,有助于提高實(shí)驗(yàn)效率,保證結(jié)果精度,具有實(shí)際的參考意義.

    表2結(jié)構(gòu)優(yōu)化結(jié)果

    H/mmR/mmS1/mmS2/mmfjVhe/m3優(yōu)化前2.007.0018.5010.500.890.01380.28優(yōu)化后2.587.4217.238.690.670.01440.22工程化分析2.607.0017.208.500.650.01390.24變化率-26.30%1.11%-13.93%

    猜你喜歡
    翅片管徑換熱器
    垂直翅片管自然對(duì)流傳熱特性的數(shù)值研究
    ASM-600油站換熱器的國(guó)產(chǎn)化改進(jìn)
    能源工程(2021年3期)2021-08-05 07:26:14
    大管徑預(yù)微導(dǎo)糾偏防護(hù)窗頂管施工技術(shù)研究
    大功率COB-LED的翅片散熱器優(yōu)化研究
    翅片管式換熱器的傳熱研究進(jìn)展
    制冷(2019年2期)2019-12-09 08:10:36
    超硬翅片滾刀加工在CNC磨床上的實(shí)現(xiàn)
    大功率LED燈的新型底部開(kāi)縫翅片散熱性能研究
    寒區(qū)某水電站大管徑PCCP管施工過(guò)程
    310S全焊式板殼型換熱器的焊接
    焊接(2016年10期)2016-02-27 13:05:35
    拓?fù)淙毕輰?duì)Armchair型小管徑多壁碳納米管輸運(yùn)性質(zhì)的影響
    亚洲片人在线观看| 国产精华一区二区三区| 伦理电影大哥的女人| 欧美性猛交╳xxx乱大交人| 怎么达到女性高潮| 亚洲av.av天堂| 中文字幕免费在线视频6| 精品人妻视频免费看| xxxwww97欧美| 91午夜精品亚洲一区二区三区 | АⅤ资源中文在线天堂| 九色国产91popny在线| 搡女人真爽免费视频火全软件 | 校园春色视频在线观看| 日韩中字成人| 欧美另类亚洲清纯唯美| avwww免费| 91麻豆av在线| 国产精品av视频在线免费观看| 亚洲欧美日韩高清专用| 亚洲av日韩精品久久久久久密| 日本一本二区三区精品| 亚洲自偷自拍三级| av专区在线播放| 深夜精品福利| 精品人妻一区二区三区麻豆 | 亚洲欧美精品综合久久99| 美女免费视频网站| 精品午夜福利视频在线观看一区| 色在线成人网| 人妻久久中文字幕网| 国产免费av片在线观看野外av| 国产免费男女视频| 欧美日韩瑟瑟在线播放| 天堂影院成人在线观看| 午夜激情欧美在线| 久久午夜福利片| 90打野战视频偷拍视频| 国产伦在线观看视频一区| 国产一级毛片七仙女欲春2| 给我免费播放毛片高清在线观看| 永久网站在线| 很黄的视频免费| 少妇丰满av| 精品久久久久久久久久免费视频| 成人无遮挡网站| 免费搜索国产男女视频| 天美传媒精品一区二区| 国产精品久久视频播放| 午夜激情福利司机影院| 18禁裸乳无遮挡免费网站照片| 欧美性感艳星| 哪里可以看免费的av片| 欧美性猛交╳xxx乱大交人| 露出奶头的视频| 最新中文字幕久久久久| 在线观看av片永久免费下载| 国产精品国产高清国产av| 淫妇啪啪啪对白视频| 人人妻人人看人人澡| 亚洲七黄色美女视频| 国内精品美女久久久久久| 首页视频小说图片口味搜索| 丰满人妻一区二区三区视频av| 9191精品国产免费久久| 午夜福利在线观看吧| 国产乱人视频| 亚洲avbb在线观看| 亚洲精品粉嫩美女一区| 欧美乱色亚洲激情| 久久性视频一级片| 2021天堂中文幕一二区在线观| 高清在线国产一区| 欧美日韩中文字幕国产精品一区二区三区| 久久性视频一级片| 91久久精品电影网| 国产精品久久久久久久久免 | 精品乱码久久久久久99久播| 国产亚洲精品综合一区在线观看| 欧美黑人巨大hd| 亚洲专区中文字幕在线| 欧美极品一区二区三区四区| 亚洲熟妇中文字幕五十中出| 内射极品少妇av片p| 给我免费播放毛片高清在线观看| 午夜福利18| 怎么达到女性高潮| 亚洲国产精品成人综合色| 99久久久亚洲精品蜜臀av| 成人亚洲精品av一区二区| 亚洲国产欧美人成| 两性午夜刺激爽爽歪歪视频在线观看| 久久久久九九精品影院| 亚洲精品日韩av片在线观看| 高潮久久久久久久久久久不卡| 国产黄色小视频在线观看| 亚洲电影在线观看av| 亚洲第一欧美日韩一区二区三区| 哪里可以看免费的av片| 一个人看的www免费观看视频| 热99re8久久精品国产| 香蕉av资源在线| 特级一级黄色大片| 亚洲国产色片| 男女做爰动态图高潮gif福利片| 亚洲国产日韩欧美精品在线观看| 嫁个100分男人电影在线观看| 精品人妻视频免费看| 国产伦精品一区二区三区四那| 一边摸一边抽搐一进一小说| 少妇人妻一区二区三区视频| 国产精品电影一区二区三区| 久久午夜福利片| 一级a爱片免费观看的视频| 国产伦精品一区二区三区四那| 日韩高清综合在线| 欧美在线黄色| 亚洲中文日韩欧美视频| 亚洲成av人片免费观看| 热99在线观看视频| 内射极品少妇av片p| 狂野欧美白嫩少妇大欣赏| 国产91精品成人一区二区三区| 一个人看视频在线观看www免费| 欧美日韩福利视频一区二区| 亚洲激情在线av| 黄色配什么色好看| aaaaa片日本免费| 精品久久国产蜜桃| 香蕉av资源在线| or卡值多少钱| 色哟哟·www| 国产av不卡久久| 十八禁网站免费在线| 亚州av有码| 12—13女人毛片做爰片一| 国内精品久久久久精免费| 欧美成狂野欧美在线观看| 久久久久久久久久黄片| 在线天堂最新版资源| www.熟女人妻精品国产| 亚洲国产精品成人综合色| 精品久久久久久久久亚洲 | 国产精品久久电影中文字幕| 制服丝袜大香蕉在线| 日韩高清综合在线| 男插女下体视频免费在线播放| 男人的好看免费观看在线视频| 中文字幕精品亚洲无线码一区| 欧美3d第一页| 成年免费大片在线观看| 97超级碰碰碰精品色视频在线观看| 18禁黄网站禁片午夜丰满| 亚洲成人久久爱视频| 狂野欧美白嫩少妇大欣赏| 国产精品嫩草影院av在线观看 | 国产av麻豆久久久久久久| 可以在线观看的亚洲视频| 国产亚洲欧美在线一区二区| 国内少妇人妻偷人精品xxx网站| 国产探花在线观看一区二区| 禁无遮挡网站| 国产av不卡久久| ponron亚洲| 国产亚洲精品综合一区在线观看| aaaaa片日本免费| 乱人视频在线观看| 国产野战对白在线观看| 国产伦在线观看视频一区| 欧美中文日本在线观看视频| 3wmmmm亚洲av在线观看| 欧美精品啪啪一区二区三区| 午夜福利成人在线免费观看| 免费人成视频x8x8入口观看| av欧美777| 九九在线视频观看精品| 亚洲av成人精品一区久久| 久9热在线精品视频| 久久久久久国产a免费观看| ponron亚洲| 韩国av一区二区三区四区| 性欧美人与动物交配| av视频在线观看入口| 黄色女人牲交| 亚洲精品在线美女| 国产精品野战在线观看| 99久久精品国产亚洲精品| 免费大片18禁| 国内揄拍国产精品人妻在线| 女同久久另类99精品国产91| 日本一本二区三区精品| 18禁黄网站禁片午夜丰满| 欧美性感艳星| 午夜a级毛片| 夜夜看夜夜爽夜夜摸| 脱女人内裤的视频| 亚洲中文字幕日韩| 亚洲av美国av| 青草久久国产| 色哟哟·www| 亚洲精品日韩av片在线观看| av福利片在线观看| 观看免费一级毛片| 日韩有码中文字幕| 午夜亚洲福利在线播放| 少妇丰满av| 日本黄色视频三级网站网址| 99精品在免费线老司机午夜| 两人在一起打扑克的视频| 91久久精品国产一区二区成人| 宅男免费午夜| 日韩免费av在线播放| 欧美+日韩+精品| 色尼玛亚洲综合影院| 久久久色成人| 一个人免费在线观看电影| 搡女人真爽免费视频火全软件 | 精品一区二区三区人妻视频| 亚洲午夜理论影院| 日韩欧美在线二视频| 好男人在线观看高清免费视频| 亚洲最大成人av| 日本 av在线| ponron亚洲| 国产白丝娇喘喷水9色精品| 国产探花在线观看一区二区| 国产单亲对白刺激| 男人的好看免费观看在线视频| 最近视频中文字幕2019在线8| avwww免费| 色哟哟哟哟哟哟| 久久精品国产亚洲av涩爱 | 中国美女看黄片| 女生性感内裤真人,穿戴方法视频| 午夜福利18| 高清日韩中文字幕在线| 国产国拍精品亚洲av在线观看| 亚洲性夜色夜夜综合| 色精品久久人妻99蜜桃| 成人性生交大片免费视频hd| 久久精品国产99精品国产亚洲性色| 国产欧美日韩一区二区精品| 亚洲欧美日韩卡通动漫| 色哟哟·www| 黄色女人牲交| 无遮挡黄片免费观看| 999久久久精品免费观看国产| 国产伦在线观看视频一区| 久久久精品大字幕| 色综合欧美亚洲国产小说| 简卡轻食公司| 性插视频无遮挡在线免费观看| АⅤ资源中文在线天堂| 亚洲狠狠婷婷综合久久图片| 天天一区二区日本电影三级| 精华霜和精华液先用哪个| 久久天躁狠狠躁夜夜2o2o| 精品一区二区三区av网在线观看| 99久久久亚洲精品蜜臀av| 一个人免费在线观看的高清视频| 国产69精品久久久久777片| 婷婷亚洲欧美| 黄片小视频在线播放| 亚洲国产色片| 国产亚洲精品综合一区在线观看| 日韩欧美精品v在线| 啪啪无遮挡十八禁网站| 日韩欧美国产在线观看| 国内精品美女久久久久久| 免费高清视频大片| 国产亚洲欧美98| 亚洲欧美日韩东京热| 日本五十路高清| 51午夜福利影视在线观看| 色尼玛亚洲综合影院| 国产一级毛片七仙女欲春2| 成人特级av手机在线观看| 国产成人aa在线观看| 久久久久免费精品人妻一区二区| 波多野结衣高清无吗| 天堂√8在线中文| 国产成人av教育| netflix在线观看网站| eeuss影院久久| 亚洲黑人精品在线| 国内精品久久久久精免费| 国产三级中文精品| 成人特级av手机在线观看| 国产成人啪精品午夜网站| 亚洲成人久久性| 一夜夜www| 黄色丝袜av网址大全| ponron亚洲| 欧美成人性av电影在线观看| 亚洲人成网站高清观看| 嫩草影视91久久| 嫩草影院精品99| 黄色丝袜av网址大全| 欧美绝顶高潮抽搐喷水| 日韩欧美精品v在线| 亚洲精品影视一区二区三区av| 夜夜躁狠狠躁天天躁| 久久久久久九九精品二区国产| 免费高清视频大片| 一区二区三区免费毛片| 国产精品野战在线观看| 国产真实伦视频高清在线观看 | 精品久久久久久,| 免费人成视频x8x8入口观看| 中文字幕免费在线视频6| 亚洲国产日韩欧美精品在线观看| 一进一出抽搐gif免费好疼| 99国产综合亚洲精品| 亚洲一区高清亚洲精品| 白带黄色成豆腐渣| 亚洲熟妇中文字幕五十中出| 在线免费观看不下载黄p国产 | 亚洲中文日韩欧美视频| 国产极品精品免费视频能看的| 天堂动漫精品| 日韩av在线大香蕉| 中文在线观看免费www的网站| 欧美+亚洲+日韩+国产| 亚洲欧美日韩卡通动漫| 亚洲在线观看片| 久久国产乱子伦精品免费另类| 青草久久国产| 两性午夜刺激爽爽歪歪视频在线观看| 美女xxoo啪啪120秒动态图 | 亚洲精品影视一区二区三区av| 亚洲va日本ⅴa欧美va伊人久久| 极品教师在线免费播放| 韩国av一区二区三区四区| 欧美日韩中文字幕国产精品一区二区三区| 天堂av国产一区二区熟女人妻| 97超视频在线观看视频| 最近最新中文字幕大全电影3| 久久亚洲精品不卡| 日本三级黄在线观看| 亚洲国产精品999在线| 一级黄片播放器| 日韩欧美 国产精品| 嫁个100分男人电影在线观看| 变态另类成人亚洲欧美熟女| 久久久久国内视频| 在线观看免费视频日本深夜| 日本一二三区视频观看| 成年女人看的毛片在线观看| 亚洲欧美日韩高清专用| 成人av一区二区三区在线看| 久久性视频一级片| 永久网站在线| 精品欧美国产一区二区三| 免费av不卡在线播放| 精品免费久久久久久久清纯| 国内精品一区二区在线观看| 少妇裸体淫交视频免费看高清| 91狼人影院| 黄色视频,在线免费观看| 日韩欧美在线乱码| 18禁在线播放成人免费| 一本综合久久免费| 午夜亚洲福利在线播放| 国内精品一区二区在线观看| 99久久无色码亚洲精品果冻| 亚洲专区国产一区二区| a在线观看视频网站| 麻豆成人av在线观看| 欧美成人一区二区免费高清观看| 女同久久另类99精品国产91| 欧美乱妇无乱码| 永久网站在线| 精品人妻熟女av久视频| 波多野结衣高清无吗| 美女 人体艺术 gogo| 国产三级中文精品| 美女高潮的动态| 亚洲精品一区av在线观看| 在线看三级毛片| 国产精品久久久久久亚洲av鲁大| 亚洲一区高清亚洲精品| 日韩中字成人| 亚洲第一电影网av| 亚洲av免费在线观看| 久久久精品大字幕| 久久久久久国产a免费观看| 麻豆一二三区av精品| 高清日韩中文字幕在线| 亚洲精品粉嫩美女一区| 成人三级黄色视频| 18禁裸乳无遮挡免费网站照片| 色综合婷婷激情| 亚洲激情在线av| 别揉我奶头~嗯~啊~动态视频| 国产黄片美女视频| 久久久久九九精品影院| 亚洲成人久久性| 简卡轻食公司| 最好的美女福利视频网| 老司机午夜福利在线观看视频| 久久国产乱子免费精品| 成人亚洲精品av一区二区| 久久亚洲精品不卡| 老女人水多毛片| 天堂动漫精品| 国产精品人妻久久久久久| 国内毛片毛片毛片毛片毛片| 天堂网av新在线| 丁香欧美五月| 最近在线观看免费完整版| 在线十欧美十亚洲十日本专区| 国产精品爽爽va在线观看网站| 国产精品不卡视频一区二区 | 黄片小视频在线播放| 亚洲精品乱码久久久v下载方式| 日韩欧美一区二区三区在线观看| 午夜久久久久精精品| 日本精品一区二区三区蜜桃| 小蜜桃在线观看免费完整版高清| 国产精品女同一区二区软件 | 欧美精品啪啪一区二区三区| 黄片小视频在线播放| 欧美另类亚洲清纯唯美| 97碰自拍视频| 久久性视频一级片| av专区在线播放| 欧美在线黄色| x7x7x7水蜜桃| 99久久久亚洲精品蜜臀av| 欧美成人一区二区免费高清观看| 成人永久免费在线观看视频| 国产一区二区激情短视频| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国产一区二区激情短视频| 国产高清三级在线| 搡女人真爽免费视频火全软件 | 男人舔奶头视频| 午夜精品在线福利| 日本成人三级电影网站| 69人妻影院| 午夜日韩欧美国产| 精品一区二区免费观看| 成年免费大片在线观看| 91久久精品电影网| 18禁裸乳无遮挡免费网站照片| 久久久久免费精品人妻一区二区| 国产国拍精品亚洲av在线观看| 国产精品久久电影中文字幕| 国产av在哪里看| 少妇高潮的动态图| av在线蜜桃| 99国产综合亚洲精品| 1024手机看黄色片| 精品人妻一区二区三区麻豆 | 久久亚洲精品不卡| 18禁黄网站禁片午夜丰满| 欧美区成人在线视频| 欧美精品国产亚洲| 色av中文字幕| 少妇人妻一区二区三区视频| av女优亚洲男人天堂| 久久国产精品影院| 亚洲av第一区精品v没综合| 真人做人爱边吃奶动态| 宅男免费午夜| 色av中文字幕| 久久国产乱子伦精品免费另类| 久9热在线精品视频| 在线观看免费视频日本深夜| 久久久久久久精品吃奶| 国产免费男女视频| 变态另类丝袜制服| 观看免费一级毛片| 亚洲激情在线av| 亚洲国产精品久久男人天堂| 午夜福利18| 日韩 亚洲 欧美在线| 日本黄色视频三级网站网址| 日日干狠狠操夜夜爽| 18+在线观看网站| 久久久久久久久久黄片| 别揉我奶头 嗯啊视频| 精品久久久久久,| 国产亚洲精品久久久com| 少妇熟女aⅴ在线视频| 少妇的逼好多水| 国产又黄又爽又无遮挡在线| 免费看a级黄色片| 成人av在线播放网站| 国产亚洲精品久久久久久毛片| 丝袜美腿在线中文| 中文资源天堂在线| 免费无遮挡裸体视频| 国产精品,欧美在线| 精品久久久久久久久久久久久| 亚洲美女视频黄频| 老司机午夜福利在线观看视频| 午夜福利高清视频| 欧美xxxx黑人xx丫x性爽| 日韩欧美在线乱码| a级毛片a级免费在线| 久久久久精品国产欧美久久久| 91av网一区二区| 特大巨黑吊av在线直播| 日韩欧美国产一区二区入口| 51国产日韩欧美| 亚洲 国产 在线| bbb黄色大片| 国产爱豆传媒在线观看| 久久九九热精品免费| 网址你懂的国产日韩在线| 男女视频在线观看网站免费| 变态另类成人亚洲欧美熟女| 久久久久九九精品影院| 999久久久精品免费观看国产| 精品欧美国产一区二区三| 日本成人三级电影网站| 日韩av在线大香蕉| 小说图片视频综合网站| 精品久久久久久久末码| 天天一区二区日本电影三级| 色播亚洲综合网| 国产精品一区二区三区四区免费观看 | 一进一出抽搐动态| 亚洲精品亚洲一区二区| 少妇人妻精品综合一区二区 | 丰满的人妻完整版| 精品人妻一区二区三区麻豆 | 高潮久久久久久久久久久不卡| 99久久精品热视频| 又黄又爽又刺激的免费视频.| 美女黄网站色视频| 中文字幕久久专区| 成人特级av手机在线观看| 免费人成在线观看视频色| 欧美一区二区国产精品久久精品| 亚洲国产日韩欧美精品在线观看| 亚洲狠狠婷婷综合久久图片| 国产精品不卡视频一区二区 | 亚洲,欧美,日韩| 51午夜福利影视在线观看| 天天一区二区日本电影三级| 嫩草影院新地址| 欧美精品啪啪一区二区三区| 天美传媒精品一区二区| 国产精品国产高清国产av| 丝袜美腿在线中文| 在线观看免费视频日本深夜| 欧美成狂野欧美在线观看| 三级国产精品欧美在线观看| 麻豆成人午夜福利视频| av福利片在线观看| 女人十人毛片免费观看3o分钟| 国产成+人综合+亚洲专区| 黄色视频,在线免费观看| 真实男女啪啪啪动态图| 久久精品久久久久久噜噜老黄 | netflix在线观看网站| 国产精品精品国产色婷婷| 伊人久久精品亚洲午夜| 久久亚洲真实| 国产乱人伦免费视频| 夜夜爽天天搞| 国产黄片美女视频| 最近最新免费中文字幕在线| 亚洲专区国产一区二区| 一区二区三区四区激情视频 | ponron亚洲| 夜夜夜夜夜久久久久| 国产三级在线视频| 成熟少妇高潮喷水视频| 少妇人妻一区二区三区视频| 狂野欧美白嫩少妇大欣赏| 国产一级毛片七仙女欲春2| 国产成人影院久久av| 国产视频一区二区在线看| 18美女黄网站色大片免费观看| 黄色配什么色好看| 欧美激情在线99| 国产黄片美女视频| 俺也久久电影网| 三级国产精品欧美在线观看| 国产美女午夜福利| 午夜免费成人在线视频| 夜夜夜夜夜久久久久| 3wmmmm亚洲av在线观看| 亚洲不卡免费看| 国产精品一及| 国产色婷婷99| 亚洲国产欧洲综合997久久,| 可以在线观看毛片的网站| 午夜福利免费观看在线| 韩国av一区二区三区四区| 国产伦人伦偷精品视频| 亚洲国产精品合色在线| 一个人观看的视频www高清免费观看| 又紧又爽又黄一区二区| 在现免费观看毛片| 中文字幕av成人在线电影| 99久久精品国产亚洲精品| 2021天堂中文幕一二区在线观| 级片在线观看| а√天堂www在线а√下载| 国内少妇人妻偷人精品xxx网站| 国产av一区在线观看免费| 美女xxoo啪啪120秒动态图 | 欧美成人a在线观看| 日本成人三级电影网站| 免费大片18禁| 日日干狠狠操夜夜爽| 国产视频一区二区在线看| 国产精品美女特级片免费视频播放器| 亚洲精品一区av在线观看| 欧美另类亚洲清纯唯美| 我的女老师完整版在线观看|