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

    基于神經(jīng)網(wǎng)絡(luò)和遺傳算法的螺旋折流板換熱器性能預(yù)測

    2016-04-11 11:02:15孫永利王華金肖曉明
    化學(xué)工業(yè)與工程 2016年4期
    關(guān)鍵詞:感知器流板殼程

    孫永利,王華金,郝 麗,肖曉明

    在各種換熱設(shè)備中,管殼式換熱器具有結(jié)構(gòu)簡單、可靠性高、使用壓力范圍廣、使用技術(shù)成熟等優(yōu)點。因此,管殼式換熱器廣泛應(yīng)用于石油、化工、輕工、冶金等過程工業(yè)以及其他工業(yè)部門[1-2]。目前,螺旋折流板換熱器逐步被人們所認(rèn)可,相比于傳統(tǒng)的弓形折流擋板換熱器它有極大的優(yōu)勢:1)增強(qiáng)了殼程換熱;2)通過殼程時壓降更低;3)減少了旁通流;4)降低了殼程污垢熱阻和流激振動[3]。過去數(shù)十年間,為了滿足對螺旋折流板換熱器的精確設(shè)計,人們在實驗和數(shù)據(jù)研究方面做了大量的研究,提出了許多有效數(shù)據(jù)以及殼程傳熱和壓降的關(guān)聯(lián)方法[4-9]。

    近年來,人工神經(jīng)網(wǎng)絡(luò)作為一種新興的優(yōu)選方法,成功地應(yīng)用于許多科學(xué)研究和工程實踐。特定的人工神經(jīng)網(wǎng)絡(luò)被用來處理多種換熱器在穩(wěn)態(tài)傳熱和水力學(xué)特性的分析、性能預(yù)測和動態(tài)控制等領(lǐng)域,并取得了不錯的效果。Diaz等[10]在使用人工神經(jīng)網(wǎng)絡(luò)對單排管翅式換熱器進(jìn)行穩(wěn)態(tài)以及動態(tài)模擬控制上開展了一系列工作。Pacheco-Vega等[11]利用神經(jīng)網(wǎng)絡(luò)對冷凍系統(tǒng)使用非常普遍的以空氣和R22冷凍劑作為工作介質(zhì)的換熱器的總換熱速率進(jìn)行了模擬。Islamoglu等[12]在對管芯式換熱器進(jìn)行換熱速率實驗基礎(chǔ)上,利用神經(jīng)網(wǎng)絡(luò)與反向傳播算法結(jié)合建立了預(yù)測模型。Hao等[13]對5種不同結(jié)構(gòu)參數(shù)的板翅式換熱器一定雷諾數(shù)范圍內(nèi)的換熱系數(shù)j因子和摩擦系數(shù)f因子進(jìn)行了實驗研究,利用所得到的40組有限實驗數(shù)據(jù)建立了板翅式換熱器j因子和f因子的預(yù)測模型。但是大多數(shù)研究都集中于管翅式換熱器,只有Xie等[14]將人工神經(jīng)網(wǎng)絡(luò)應(yīng)用到了對弓形板管殼式換熱器和連續(xù)型螺旋折流板管殼式換熱器換熱性能的預(yù)測分析上。在螺旋折流板換熱器殼程換熱和流體力學(xué)方面的分析,人工神經(jīng)網(wǎng)絡(luò)還沒有大規(guī)模應(yīng)用。

    因此,本研究利用中試試驗所得到的光滑管和橫槽管非連續(xù)螺旋折流板換熱器殼程換熱系數(shù)與壓降數(shù)據(jù),建立了9個輸入?yún)?shù)的采取遺傳算法優(yōu)化的多層感知器神經(jīng)網(wǎng)絡(luò)預(yù)測模型,比較了遺傳算法優(yōu)化與否的預(yù)測誤差以及訓(xùn)練預(yù)測過程中的誤差,并將其泛化能力與回歸關(guān)聯(lián)式進(jìn)行了對比。將試驗與模擬所得到的數(shù)據(jù)結(jié)合,進(jìn)一步地對模型進(jìn)行訓(xùn)練,建立了適用范圍更廣的螺旋折流板換熱器殼程換熱系數(shù)與壓降預(yù)測模型。

    1 螺旋折流板換熱器中試試驗

    螺旋折流板換熱器中試試驗裝置流程圖如圖1所示。系統(tǒng)包括2個獨立的循環(huán):冷卻工作介質(zhì)及加熱工作介質(zhì)循環(huán)系統(tǒng),建立了管殼式換熱器中試試驗平臺對3臺螺旋角分別為7°、13°和25°的螺旋折流板管殼式換熱器的殼程流阻和傳熱特性進(jìn)行研究。此外中試試驗中采用了2種類型的換熱管,分別是普通光管和橫槽換熱管,測試了2種換熱管對于殼程傳熱與流阻性能的影響。螺旋折流板換熱器的結(jié)構(gòu)參數(shù)見表1。利用中試試驗所得到的光滑管和橫槽管螺旋折流板換熱器殼程換熱系數(shù)與壓降數(shù)據(jù),建立采取遺傳算法優(yōu)化的多層感知器神經(jīng)網(wǎng)絡(luò)預(yù)測模型。

    圖1 中試試驗系統(tǒng)流程圖Fig.1 Experimental system schematic

    表1 螺旋折流板換熱器的結(jié)構(gòu)參數(shù)Table 1 Parameters of the tested shell-and-tube heat exchangers with helical baffle

    試驗過程中,殼程入口體積流量取值范圍設(shè)置為50~150 m3/h。管程入口體積流量取值范圍設(shè)置為 50~90 m3/h。

    2 多層感知器優(yōu)化模型的參數(shù)選擇

    本研究建立的多層感知器預(yù)測模型是在利用Matlab2010A軟件中的神經(jīng)網(wǎng)絡(luò)工具包基礎(chǔ)上編程實現(xiàn)的。所建立的網(wǎng)絡(luò)結(jié)構(gòu)中,輸入層中包含9個獨立變量,分別是殼程體積流量Vs、管程體積流量Vt、殼程流體進(jìn)口溫度Ts、管程流體進(jìn)口溫度Tt、螺旋角β、螺距B、螺旋周期數(shù)S、管子槽深e和槽距tp。其中,殼程與管程工作介質(zhì)都是水,殼程入口溫度40~60℃,管程入口溫度70~90℃,其余輸入變量的取值范圍見表2。輸出層包含有2個輸出變量,分別是殼程換熱系數(shù)和殼程壓降,是評價換熱器性能的主要指標(biāo)。

    將換熱器中試試驗中得到的562組實驗數(shù)據(jù)按比例隨機(jī)分成3組,其中的60%用來訓(xùn)練網(wǎng)絡(luò),20%用于驗證網(wǎng)絡(luò)的過擬合性,20%用于測試網(wǎng)絡(luò)性能,即它的泛化能力。多層感知器模型隱含層傳遞函數(shù)設(shè)為Sigmoid函數(shù),輸出層傳遞函數(shù)設(shè)為線性函數(shù)。為了有效的訓(xùn)練網(wǎng)絡(luò),將輸入輸出變量統(tǒng)一歸一化到[0.15,0.85]范圍。訓(xùn)練過程中采用了Levenberg-Marquardt(TRAINLM)函數(shù)作為訓(xùn)練函數(shù)。驗證過程最大迭代失敗次數(shù)設(shè)為20次,學(xué)習(xí)速率設(shè)定為0.01。為了有效地評價多層感知器模型對于殼程換熱系數(shù)以及壓降訓(xùn)練與預(yù)測時的誤差,選擇平均相對誤差作為評價標(biāo)準(zhǔn),按公式(1)計算:

    其中Ae是實驗結(jié)果,Ap是預(yù)測結(jié)果,M是數(shù)據(jù)總數(shù)。

    網(wǎng)絡(luò)的預(yù)測能力,即泛化能力通過預(yù)測輸出值和實驗值之間的平均準(zhǔn)確度R和分散度σ來比較。

    對于用來優(yōu)化多層感知器網(wǎng)絡(luò)的遺傳算法,選擇預(yù)測值和給定輸出值之間的平均相對誤差的絕對值之和作為其適應(yīng)度函數(shù),計算方法如公式(4)。

    在遺傳算法的優(yōu)化過程中,操作選擇輪盤賭算法,種群規(guī)模設(shè)為40,進(jìn)化次數(shù)設(shè)為200,交叉操作和變異操作的概率分別設(shè)置為0.4和0.1。

    表2 結(jié)構(gòu)參數(shù)取值范圍Table 2 Range of structural parameters

    利用遺傳算法優(yōu)化多層感知器神經(jīng)網(wǎng)絡(luò)模型的流程如圖2所示。

    圖2 遺傳算法優(yōu)化多層感知器模型流程Fig.2 Procedures of the MLP network optimization with the genetic algorithm

    3 預(yù)測結(jié)果與討論

    3.1 遺傳算法優(yōu)化多層感知器模型結(jié)果

    圖3顯示的是遺傳算法優(yōu)化多層感知器網(wǎng)絡(luò)過程中的適應(yīng)度值變化。采用的網(wǎng)絡(luò)結(jié)構(gòu)是9-7-5-2,即一層7個節(jié)點的第一隱含層,一層5個節(jié)點的第二隱含層。在初始階段,適應(yīng)度值很大,說明當(dāng)前種群所代表的初始權(quán)值與閥值設(shè)置不適應(yīng)所建立的網(wǎng)絡(luò),會引起較大的訓(xùn)練及預(yù)測誤差。隨著這些種群個體在前期進(jìn)化過程中被不斷淘汰,誤差總和開始迅速下降,到一定進(jìn)化步驟之后逐漸變化為階梯型下降,在大約110步之后,曲線開始變得平緩,并且達(dá)到其最小值,說明優(yōu)化過程結(jié)束。

    圖3 適應(yīng)度值變化曲線Fig.3 The process of evolution

    圖4a)和圖4b)顯示的是經(jīng)過優(yōu)化的多層感知器模型(MLP-GA)與未經(jīng)優(yōu)化的模型(MLP)在對換熱器殼程換熱系數(shù)和壓降進(jìn)行預(yù)測時的相對誤差分布圖。紅色星號代表的是未經(jīng)優(yōu)化的多層感知器神經(jīng)網(wǎng)絡(luò),而藍(lán)色圓圈代表的是經(jīng)過遺傳算法優(yōu)化后的多層感知器模型。

    圖4 對換熱器殼程換熱系數(shù)和壓降進(jìn)行預(yù)測時的相對誤差分布圖Fig.4 Rel ative error scatter of and between the predictions and experiment data

    從圖4中可以看出,經(jīng)過優(yōu)化的網(wǎng)絡(luò)模型顯示出了更高的預(yù)測精確度和穩(wěn)健性,克服了初始權(quán)值與閥值的隨機(jī)選擇對其預(yù)測能力造成的影響。雖然采取遺傳算法進(jìn)行優(yōu)化會減慢多層感知器網(wǎng)絡(luò)的訓(xùn)練速度,但是這種優(yōu)化策略能夠使模型預(yù)測能力得到顯著增強(qiáng)。因此接下來所采用的不同結(jié)構(gòu)的多層感知器網(wǎng)絡(luò)均使用遺傳算法進(jìn)行優(yōu)化。

    3.2 遺傳算法優(yōu)化的多層感知器模型訓(xùn)練預(yù)測結(jié)果

    本研究選擇7種多層感知器模型網(wǎng)絡(luò)結(jié)構(gòu)對實驗數(shù)據(jù)進(jìn)行訓(xùn)練與測試,表3和表4分別顯示的是這7種不同網(wǎng)絡(luò)結(jié)構(gòu)模型對殼程換熱系數(shù)和壓降的預(yù)測結(jié)果。從表3和表4中可以發(fā)現(xiàn),對于單隱含層網(wǎng)絡(luò),當(dāng)節(jié)點數(shù)從5增加到7時,訓(xùn)練過程的平均相對誤差值(MRE)發(fā)生下降,這說明隱含層內(nèi)節(jié)點數(shù)的適當(dāng)增加能夠提高訓(xùn)練過程的準(zhǔn)確度。而對于雙隱含層網(wǎng)絡(luò),當(dāng)?shù)谝浑[含層節(jié)點數(shù)由6增加到8的時候,預(yù)測輸出值和實驗值之間的平均準(zhǔn)確度(R)和分散度(σ)并沒有隨之降低,所以單純地提高節(jié)點數(shù)不一定意味著網(wǎng)絡(luò)預(yù)測性能的提高。盡管9-8-5-2的網(wǎng)絡(luò)結(jié)構(gòu)訓(xùn)練過程中的MRE值最低,但是泛化能力表現(xiàn)卻不是最好的。而9-7-5-2網(wǎng)絡(luò)結(jié)構(gòu)訓(xùn)練過程對于殼程換熱系數(shù)以及壓降預(yù)測的平均相對誤差分別為1.0012%和2.0423%,雖然不是最好的,但是其測試過程 R分別為1.0143和1.0291,σ分別為0.0605和0.0956,是所有網(wǎng)絡(luò)結(jié)構(gòu)中最低的,說明其泛化能力是最好的。因此,在預(yù)測本試驗螺旋折流板換熱器殼程換熱系數(shù)和殼程壓降的多層感知器模型中,9-7-5-2是最佳的網(wǎng)絡(luò)結(jié)構(gòu)。

    表3 優(yōu)化的多層感知器模型不同網(wǎng)絡(luò)結(jié)構(gòu)預(yù)測殼程換熱系數(shù)誤差Table 3 Performance comparison of different MLP topology for heat transfer rate on shell side

    圖5a)和圖5b)顯示的是由經(jīng)過優(yōu)化的多層感知器模型訓(xùn)練過程所預(yù)測的殼程換熱系數(shù)與壓降值同實驗值的對比。

    表4 優(yōu)化的多層感知器模型不同網(wǎng)絡(luò)結(jié)構(gòu)預(yù)測殼程壓降誤差Table 4 Performance comparison of different MLP topology for pressure drop on shell side

    圖5 優(yōu)化的多層感知器模型訓(xùn)練過程預(yù)測值與實驗值對比Fig.5 Prediction of shell-side heat transfer rate and pressure drop by the MLP-GA with training data

    由圖5可以看到,大部分訓(xùn)練過程的預(yù)測輸出值都非常接近于對角線,即訓(xùn)練過程預(yù)測輸出值與實驗值接近一致。

    由于多層感知器模型最重要的是其泛化能力,即對于訓(xùn)練數(shù)據(jù)范圍以外數(shù)據(jù)預(yù)測的準(zhǔn)確度,因此不僅要求其對于訓(xùn)練數(shù)據(jù)預(yù)測準(zhǔn)確度,更要考察其對于測試數(shù)據(jù)預(yù)測的準(zhǔn)確度。

    圖6a)和圖6b)顯示的是經(jīng)過優(yōu)化的多層感知器模型測試階段所預(yù)測的殼程換熱系數(shù)與壓降值同實驗值的對比??梢园l(fā)現(xiàn),基本上預(yù)測輸出值都很接近對角線,說明預(yù)測輸出值與實驗值接近一致,證明了我們所建立的遺傳算法優(yōu)化的多層感知器網(wǎng)絡(luò)在預(yù)測螺旋折流板換熱器殼程換熱系數(shù)和壓降時具有良好的泛化能力。

    圖6 優(yōu)化的多層感知器模型測試過程預(yù)測值與實驗值對比Fig.6 Prediction of shell-side heat transfer rate and pressure drop by the MLP-GA with testing data

    3.3 優(yōu)化的多層感知器模型與關(guān)聯(lián)式預(yù)測精度對比

    將中試試驗所得到的螺旋折流板換熱器水-水試驗數(shù)據(jù)做回歸處理得到努賽爾數(shù)(Nu)及摩擦因子(f),形如 Nu=CRemPr1/3以及 f=φReω的關(guān)聯(lián)式,適用范圍是 8000<Re<12000,2.4<Pr<4.1。殼程進(jìn)口溫度為40℃,管程進(jìn)口溫度為70℃時,3種不同螺旋角的關(guān)聯(lián)式如下所示:

    圖7 優(yōu)化的多層感知器網(wǎng)絡(luò)和關(guān)聯(lián)式預(yù)測精度對比Fig.7 Comparison of shell-side heat transfer rate and pressure drop predictions by the MLP-GA and correlations

    圖7a)和圖7b)顯示的是優(yōu)化的多層感知器預(yù)測模型同殼程進(jìn)口溫度為40℃時的回歸關(guān)聯(lián)式在預(yù)測精度上的對比。由圖7可以看出,不管是對于殼程換熱系數(shù)還是殼程壓降,回歸關(guān)聯(lián)式計算結(jié)果與真實值相差較大。而利用多層感知器預(yù)測模型所得到的預(yù)測結(jié)果具有很好的精確度。

    從表3和表4中也可以看出,回歸關(guān)聯(lián)式所得到的計算結(jié)果精確度要遠(yuǎn)遠(yuǎn)低于我們所測試的所有結(jié)構(gòu)類型的多層感知器預(yù)測網(wǎng)絡(luò)。

    此外,通過多層感知器預(yù)測網(wǎng)絡(luò)能夠直接利用輸入變量的數(shù)值得到對應(yīng)的殼程換熱系數(shù)和壓降,而關(guān)聯(lián)式需要先計算出殼程傳熱Nu數(shù)以及摩擦因子f,然后再將其轉(zhuǎn)化為殼程換熱系數(shù)以及壓降值。因此,不管是從預(yù)測的準(zhǔn)確度,還是方便性來講,多層感知器預(yù)測網(wǎng)絡(luò)與回歸關(guān)聯(lián)式相比,都具有很大優(yōu)勢。

    本研究建立的遺傳算法優(yōu)化的多層感知器模型(MLP-GA)可以用來預(yù)測輸入?yún)?shù)對輸出的影響。圖8a)和圖8b)顯示使用MLP-GA方法在訓(xùn)練數(shù)據(jù)范圍對換熱器殼程不同體積流量下的殼程換熱系數(shù)及殼程壓降的預(yù)測。

    圖8 遺傳算法優(yōu)化的多層感知器模型預(yù)測不同體積流量下?lián)Q熱器殼程換熱系數(shù)及殼程壓降Fig.8 Predictions for shell-side heat transfer rate and pressure drop vs.shell-side volumetric flow rate by MLP-GA

    由圖8可以看出,6種型式換熱器的殼程換熱系數(shù)均隨體積流量的增加而增大,同一體積流量下,螺旋角越小,殼程換熱系數(shù)越高。在相同螺旋角下,橫槽管換熱器的殼程換熱系數(shù)較光滑管高。相同體積流量下,7°螺旋角換熱器殼程壓降明顯高于其他幾種型式。相同螺旋角的情況下,橫槽管換熱器殼程壓降均高于光滑管換熱器。殼程傳熱系數(shù)及殼程壓降的網(wǎng)絡(luò)預(yù)測值與實驗值規(guī)律一致。

    4 結(jié)論

    通過換熱器中試試驗平臺考察了不同螺旋角橫槽管和光滑管螺旋折流板換熱器的殼程性能,利用中試試驗數(shù)據(jù)建立了預(yù)測螺旋折流板換熱器殼程換熱系數(shù)與壓降的多層感知器神經(jīng)網(wǎng)絡(luò)模型,并利用遺傳算法對多層感知器預(yù)測模型進(jìn)行了優(yōu)化,得到以下結(jié)論。

    1)遺傳算法優(yōu)化的多層感知器模型預(yù)測精度要高于未經(jīng)優(yōu)化的模型,遺傳算法優(yōu)化能夠提高多層感知器模型的預(yù)測準(zhǔn)確度和穩(wěn)健性,克服了初始權(quán)值和閥值隨機(jī)選擇的不利影響。

    2)多層感知器模型的網(wǎng)絡(luò)結(jié)構(gòu)要根據(jù)實際應(yīng)用中的預(yù)測精度來進(jìn)行選擇。經(jīng)過比較,我們選擇了訓(xùn)練誤差和泛化性能綜合起來最好的9-7-5-2的網(wǎng)絡(luò)結(jié)構(gòu),此結(jié)構(gòu)訓(xùn)練過程對于殼程換熱系數(shù)和壓降預(yù)測的平均相對誤差分別為1.0012%和2.0432%,測試過程預(yù)測的平均準(zhǔn)確度為1.0143和1.0291,分散度為0.065和0.0956。

    3)遺傳算法優(yōu)化的多層感知器網(wǎng)絡(luò)模型的預(yù)測精度要高于回歸關(guān)聯(lián)式。

    4)采用遺傳算法優(yōu)化的多層感知器模型(MLPGA)預(yù)測體積流量對殼程換熱系數(shù)及殼程壓降的影響。殼程換熱系數(shù)均隨體積流量的增加而增大,且螺旋角越小,殼程換熱系數(shù)越高,橫槽管換熱器的殼程換熱系數(shù)較光滑管高。7°螺旋角換熱器殼程壓降明顯高于13°和25°,橫槽管換熱器殼程壓降高于光滑管換熱器。

    參考文獻(xiàn):

    [1]Master B,Chunangad K,Boxma A,et al.Most frequently used heat exchangers from pioneering research to worldwide applications[J].Heat Transfer Engineering,2006,27(6):4-11

    [2]Stehl K P,Wadekar V V.Different strategies to improve industrial heat exchange[J].Heat Transfer Engineering,2002,23(6):36-48

    [3]Chen G,Wang Q.Experimental and numerical studies of shell-and-tube heat exchangers with helical baffles[J].ASME,2009:601-609

    [4]Jafari N M R,Shafeghat A.Fluid flow analysis and extension of rapid design algorithm for helical baffle heat exchangers[J].Applied Thermal Engineering,2008,28(11/12):1 324-1 332

    [5]Kral D,Stehlik P,Van Der Ploeg H,et al.Helical baffles in shell-and-tube heat exchangers,Part I:Experimental verification[J].Heat Transfer Engineering,1996,17(1):93-101

    [6]Lei Y,He Y,Chu P,et al.Design and optimization of heat exchangers with helical baffles[J].Chemical Engineering Science,2008,63(17):4 386-4 395

    [7]Lutcha J,Nemcansky J.Performance improvement of tubular heat exchangers by helical baffles[J].Chemical Engineering Research& Design,1990,68:263-270

    [8]Peng B,Wang Q,Zhang C,et al.An experimental study of shell-and-tube heat exchangers with continuous helical baffles[J].Journal of Heat Transfer,2007,129(10):1 425-1 431

    [9]Zhang J,Li B,Huang W,et al.Experimental performance comparison of shell-side heat transfer for shell-andtube heat exchangers with middle-overlapped helical baffles and segmental baffles[J].Chemical Engineering Science,2009,64(8):1 643-1 653

    [10]Diaz G,Sen M,Yang K,et al.Dynamic prediction and control of heat exchangersusingartificial neural networks[J].International Journal of Heat and Mass Transfer,2001,45:1 671-1 679

    [11]Pacheco-Vega A,Sen M,Yang K,et al.Neural network analysis of fin-tube refrigerating heat exchanger with limited experimental data[J].International Journal of Heat and Mass Transfer,2001,44(4):763-770

    [12]Islamoglu Y,Kurt A,Parmaksizoglu C.Performance prediction for non-adiabatic capillary tube suction line heat exchanger an artificial neural network[J].Energy Conversion and Management,2005,46:223-232

    [13]Hao P,Xiang L.Neural networks analysis of thermal characteristics on plate-fin heat exchangers with limited experimental data[J].Applied Thermal Engineering,2009,29:2 251-2 256

    [14]Xie G,Wang Q,Zeng M,et al.Heat transfer analysis for shell-and-tube heat exchangers with experimental data by artificial neural networks approach[J].Applied Thermal Engineering,2007,27(5/6):1 096-1 104

    猜你喜歡
    感知器流板殼程
    結(jié)構(gòu)參數(shù)對交錯百葉折流板管殼式換熱器性能影響的研究
    新型錐形孔折流板管殼式換熱器數(shù)值模擬研究
    火箭貯箱噴涂機(jī)器人直接示教運(yùn)動感知器設(shè)計與分析
    管殼式換熱器殼程的傳熱強(qiáng)化探討
    三葉膨脹管換熱器殼程強(qiáng)化傳熱的數(shù)值研究*
    感知器在礦井突水水源識別中的應(yīng)用
    AI超市
    尿濕感知器
    傾斜折流柵式換熱器殼程熱力特性研究
    換熱管束排列對換熱器殼程流動與傳熱影響的研究
    免费一级毛片在线播放高清视频| 免费高清视频大片| 午夜免费激情av| 91麻豆av在线| 老熟妇乱子伦视频在线观看| 精品久久久久久,| 有码 亚洲区| 一区二区三区免费毛片| 亚洲不卡免费看| 内射极品少妇av片p| 日韩欧美国产一区二区入口| 90打野战视频偷拍视频| 日韩欧美精品免费久久 | 非洲黑人性xxxx精品又粗又长| 精品久久久久久久久久免费视频| 日本成人三级电影网站| 一区福利在线观看| 免费看十八禁软件| 久久人妻av系列| 国产精品1区2区在线观看.| 五月伊人婷婷丁香| 69人妻影院| 给我免费播放毛片高清在线观看| 国产 一区 欧美 日韩| 国语自产精品视频在线第100页| 少妇人妻一区二区三区视频| 天美传媒精品一区二区| 成熟少妇高潮喷水视频| 网址你懂的国产日韩在线| 美女免费视频网站| 欧美黑人巨大hd| 色综合欧美亚洲国产小说| 国产精品亚洲美女久久久| 亚洲真实伦在线观看| 淫秽高清视频在线观看| 国产日本99.免费观看| 国产精品野战在线观看| 午夜影院日韩av| 国产一区二区三区在线臀色熟女| e午夜精品久久久久久久| 亚洲午夜理论影院| 亚洲成人中文字幕在线播放| 女人十人毛片免费观看3o分钟| 欧美3d第一页| 婷婷亚洲欧美| 亚洲激情在线av| 在线观看午夜福利视频| 久久久精品大字幕| 欧美+亚洲+日韩+国产| 国产亚洲欧美在线一区二区| 18禁黄网站禁片午夜丰满| 日日摸夜夜添夜夜添小说| 成人av一区二区三区在线看| 国产亚洲欧美在线一区二区| 欧美中文日本在线观看视频| 特级一级黄色大片| 在线观看一区二区三区| 国产精品久久视频播放| 怎么达到女性高潮| 悠悠久久av| 亚洲内射少妇av| 久久九九热精品免费| 国产色爽女视频免费观看| 九九在线视频观看精品| 黄片小视频在线播放| 久久久久九九精品影院| 99精品在免费线老司机午夜| 亚洲 国产 在线| 久久精品综合一区二区三区| 欧美黄色淫秽网站| 国产黄片美女视频| 欧美午夜高清在线| 久久亚洲真实| 国产成人a区在线观看| 日本黄色视频三级网站网址| 亚洲精品粉嫩美女一区| 夜夜夜夜夜久久久久| 亚洲一区二区三区不卡视频| avwww免费| 99久久无色码亚洲精品果冻| 久久久久久久亚洲中文字幕 | 日韩欧美三级三区| 精品欧美国产一区二区三| 精品一区二区三区av网在线观看| 大型黄色视频在线免费观看| 国产伦人伦偷精品视频| 国产精品久久久久久久久免 | 国产99白浆流出| 亚洲av美国av| 成人性生交大片免费视频hd| 中文资源天堂在线| 婷婷六月久久综合丁香| 国产精品嫩草影院av在线观看 | 亚洲成人免费电影在线观看| 久久香蕉国产精品| 国产精品久久久久久久久免 | 一进一出抽搐gif免费好疼| 欧美黄色片欧美黄色片| 亚洲国产色片| 久久精品亚洲精品国产色婷小说| 国产精品久久久人人做人人爽| 嫩草影院入口| 看免费av毛片| 久久国产乱子伦精品免费另类| 午夜亚洲福利在线播放| 欧美大码av| 精品一区二区三区人妻视频| 俄罗斯特黄特色一大片| 亚洲精品乱码久久久v下载方式 | 少妇人妻精品综合一区二区 | 黄色成人免费大全| 成年女人永久免费观看视频| 人妻久久中文字幕网| 亚洲色图av天堂| 国产欧美日韩一区二区三| 高清毛片免费观看视频网站| 波多野结衣巨乳人妻| 人人妻人人澡欧美一区二区| 亚洲精品久久国产高清桃花| 麻豆一二三区av精品| 人人妻,人人澡人人爽秒播| 我要搜黄色片| 亚洲av成人精品一区久久| 日韩欧美 国产精品| 国产一区二区在线av高清观看| xxxwww97欧美| 麻豆成人av在线观看| 亚洲av免费在线观看| 一进一出抽搐动态| 亚洲色图av天堂| 亚洲av免费在线观看| 51国产日韩欧美| 久久精品亚洲精品国产色婷小说| 99热只有精品国产| 国产亚洲欧美98| 国产乱人伦免费视频| 久久中文看片网| 99精品久久久久人妻精品| 午夜日韩欧美国产| 色噜噜av男人的天堂激情| 成人无遮挡网站| 日韩大尺度精品在线看网址| 美女免费视频网站| 午夜精品久久久久久毛片777| 亚洲国产日韩欧美精品在线观看 | 黑人欧美特级aaaaaa片| 久久人人精品亚洲av| 欧美zozozo另类| 国产aⅴ精品一区二区三区波| 悠悠久久av| a级毛片a级免费在线| 欧美在线一区亚洲| 在线视频色国产色| 国产高清视频在线观看网站| 俄罗斯特黄特色一大片| tocl精华| 欧美成人免费av一区二区三区| 国产 一区 欧美 日韩| 国产一区二区三区视频了| 亚洲无线在线观看| 亚洲一区高清亚洲精品| 熟妇人妻久久中文字幕3abv| 亚洲国产精品999在线| 国内毛片毛片毛片毛片毛片| 男插女下体视频免费在线播放| 99久国产av精品| 好看av亚洲va欧美ⅴa在| 午夜福利高清视频| 国内揄拍国产精品人妻在线| 国产精品日韩av在线免费观看| 欧美xxxx黑人xx丫x性爽| 十八禁人妻一区二区| 一级黄片播放器| 免费看日本二区| 小蜜桃在线观看免费完整版高清| 看黄色毛片网站| 激情在线观看视频在线高清| 亚洲中文字幕一区二区三区有码在线看| 精品电影一区二区在线| 日韩 欧美 亚洲 中文字幕| 久久香蕉国产精品| 国产亚洲精品久久久com| 少妇的丰满在线观看| 日韩国内少妇激情av| 欧美黄色片欧美黄色片| 亚洲欧美精品综合久久99| 999久久久精品免费观看国产| 90打野战视频偷拍视频| 精品久久久久久久久久免费视频| 在线观看一区二区三区| 老司机午夜十八禁免费视频| av福利片在线观看| 亚洲av一区综合| 级片在线观看| 丁香欧美五月| 国产中年淑女户外野战色| 手机成人av网站| 亚洲精品国产精品久久久不卡| 欧美一级毛片孕妇| 亚洲av成人不卡在线观看播放网| 日本一二三区视频观看| 国内精品一区二区在线观看| 淫秽高清视频在线观看| 国产aⅴ精品一区二区三区波| 午夜福利在线在线| 亚洲成av人片免费观看| 国语自产精品视频在线第100页| 国产三级黄色录像| 两个人看的免费小视频| 欧美bdsm另类| 免费在线观看亚洲国产| 色哟哟哟哟哟哟| 国产高清有码在线观看视频| 亚洲av中文字字幕乱码综合| 草草在线视频免费看| 熟女人妻精品中文字幕| 在线天堂最新版资源| 偷拍熟女少妇极品色| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 欧美最新免费一区二区三区 | 一个人观看的视频www高清免费观看| 人妻丰满熟妇av一区二区三区| 亚洲成人中文字幕在线播放| 操出白浆在线播放| 精品一区二区三区视频在线观看免费| 欧美成人a在线观看| 欧美高清成人免费视频www| 亚洲精品成人久久久久久| 亚洲精品影视一区二区三区av| 日日摸夜夜添夜夜添小说| 久久精品夜夜夜夜夜久久蜜豆| 国产精品美女特级片免费视频播放器| 中文亚洲av片在线观看爽| 听说在线观看完整版免费高清| 久久国产精品人妻蜜桃| 国产一区二区在线观看日韩 | 一个人看的www免费观看视频| 亚洲精品国产精品久久久不卡| 法律面前人人平等表现在哪些方面| 久久精品人妻少妇| 亚洲精品久久国产高清桃花| 手机成人av网站| or卡值多少钱| 国产激情欧美一区二区| 九九久久精品国产亚洲av麻豆| 美女 人体艺术 gogo| 在线a可以看的网站| or卡值多少钱| 午夜视频国产福利| 亚洲精品色激情综合| 亚洲av不卡在线观看| 精品人妻偷拍中文字幕| 久久九九热精品免费| 国产精品国产高清国产av| 久久久久久久午夜电影| 午夜精品一区二区三区免费看| 最新中文字幕久久久久| 欧美黄色淫秽网站| 人人妻,人人澡人人爽秒播| 男插女下体视频免费在线播放| 亚洲欧美日韩东京热| 国产三级黄色录像| 99热精品在线国产| 亚洲av美国av| 欧洲精品卡2卡3卡4卡5卡区| 亚洲国产精品999在线| 可以在线观看的亚洲视频| 免费av不卡在线播放| 偷拍熟女少妇极品色| 国产探花在线观看一区二区| 日本精品一区二区三区蜜桃| 一本精品99久久精品77| 99久久精品热视频| 欧美中文综合在线视频| 欧美在线黄色| 两人在一起打扑克的视频| 身体一侧抽搐| 嫁个100分男人电影在线观看| 少妇熟女aⅴ在线视频| 久久久久久久午夜电影| 可以在线观看毛片的网站| 香蕉av资源在线| 国产精品久久久久久久久免 | 国产蜜桃级精品一区二区三区| 国产69精品久久久久777片| 国产免费一级a男人的天堂| 天天躁日日操中文字幕| 亚洲国产精品合色在线| 国产伦人伦偷精品视频| 欧美+日韩+精品| 国产亚洲欧美在线一区二区| 99久久精品一区二区三区| 高清日韩中文字幕在线| 天堂影院成人在线观看| 香蕉av资源在线| 欧美另类亚洲清纯唯美| 男人舔女人下体高潮全视频| 色综合欧美亚洲国产小说| 亚洲精品一区av在线观看| 欧美日韩一级在线毛片| 三级国产精品欧美在线观看| 日本一二三区视频观看| 一区二区三区激情视频| 男女视频在线观看网站免费| 精品福利观看| 91久久精品国产一区二区成人 | 欧美性猛交黑人性爽| 国产色爽女视频免费观看| 久久这里只有精品中国| 国产精品 欧美亚洲| 美女免费视频网站| 欧美成人免费av一区二区三区| 国产三级在线视频| 国产97色在线日韩免费| 97碰自拍视频| 亚洲欧美激情综合另类| 女人高潮潮喷娇喘18禁视频| 桃色一区二区三区在线观看| 欧美日本亚洲视频在线播放| 亚洲一区高清亚洲精品| 欧美av亚洲av综合av国产av| 免费看日本二区| 亚洲 国产 在线| 亚洲五月天丁香| 窝窝影院91人妻| av中文乱码字幕在线| 国产精品女同一区二区软件 | 九九久久精品国产亚洲av麻豆| 麻豆成人午夜福利视频| 九九在线视频观看精品| 桃红色精品国产亚洲av| 亚洲 欧美 日韩 在线 免费| 午夜福利视频1000在线观看| 亚洲人成网站在线播放欧美日韩| 男插女下体视频免费在线播放| 三级男女做爰猛烈吃奶摸视频| 国产伦人伦偷精品视频| 少妇熟女aⅴ在线视频| 少妇裸体淫交视频免费看高清| 91在线观看av| 最近视频中文字幕2019在线8| 麻豆一二三区av精品| 国产一区二区在线av高清观看| 国内精品久久久久精免费| 男女床上黄色一级片免费看| 脱女人内裤的视频| 国产亚洲欧美在线一区二区| 99久久99久久久精品蜜桃| 琪琪午夜伦伦电影理论片6080| 日日夜夜操网爽| 少妇的逼好多水| 亚洲在线观看片| 黄片小视频在线播放| 国产aⅴ精品一区二区三区波| 免费av毛片视频| 精品久久久久久久久久久久久| 久久婷婷人人爽人人干人人爱| 欧美极品一区二区三区四区| 婷婷六月久久综合丁香| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国产av不卡久久| 一级毛片女人18水好多| 99在线人妻在线中文字幕| 精品国内亚洲2022精品成人| 亚洲国产欧洲综合997久久,| 丁香六月欧美| 露出奶头的视频| 在线国产一区二区在线| 免费观看人在逋| 国产伦一二天堂av在线观看| 两个人视频免费观看高清| 岛国在线观看网站| 欧美黑人巨大hd| 免费av不卡在线播放| 一区二区三区免费毛片| 精品国产亚洲在线| 久99久视频精品免费| 日韩欧美 国产精品| 亚洲人成电影免费在线| 国产三级中文精品| xxxwww97欧美| av国产免费在线观看| 成熟少妇高潮喷水视频| www.色视频.com| 午夜久久久久精精品| 高清日韩中文字幕在线| 看免费av毛片| 午夜福利免费观看在线| 国产免费男女视频| 91在线精品国自产拍蜜月 | 国产亚洲精品久久久com| 不卡一级毛片| 国产成人aa在线观看| 久久久久久九九精品二区国产| 国产真实乱freesex| 日韩人妻高清精品专区| xxxwww97欧美| 18禁黄网站禁片免费观看直播| 人妻久久中文字幕网| av欧美777| 久久久久精品国产欧美久久久| 久久中文看片网| 国产色婷婷99| 中文字幕熟女人妻在线| 免费看光身美女| 最近最新中文字幕大全免费视频| 又紧又爽又黄一区二区| 精品一区二区三区视频在线 | 亚洲av五月六月丁香网| 午夜免费激情av| 色哟哟哟哟哟哟| 久99久视频精品免费| 欧美又色又爽又黄视频| 国产成人欧美在线观看| 99久久精品热视频| 国内精品久久久久精免费| 国产一区二区激情短视频| 中文字幕高清在线视频| 日本a在线网址| 国产精品国产高清国产av| 一区二区三区免费毛片| 两个人的视频大全免费| 97超级碰碰碰精品色视频在线观看| 免费看a级黄色片| 99久国产av精品| 搡女人真爽免费视频火全软件 | 久久亚洲精品不卡| 男女那种视频在线观看| 欧美一区二区国产精品久久精品| 亚洲无线观看免费| 亚洲av第一区精品v没综合| h日本视频在线播放| 亚洲五月婷婷丁香| 成人国产一区最新在线观看| 久久久色成人| 亚洲电影在线观看av| 久久性视频一级片| 亚洲 国产 在线| 精品欧美国产一区二区三| 黄色成人免费大全| 好男人在线观看高清免费视频| 日韩精品青青久久久久久| 国产亚洲欧美98| 男女下面进入的视频免费午夜| 一本精品99久久精品77| 欧美一区二区国产精品久久精品| 日本熟妇午夜| 日韩欧美精品免费久久 | 女生性感内裤真人,穿戴方法视频| 99久久99久久久精品蜜桃| 国产高清有码在线观看视频| 国产精品久久久久久久电影 | 美女免费视频网站| 国产综合懂色| 露出奶头的视频| 757午夜福利合集在线观看| www.www免费av| 老汉色∧v一级毛片| 亚洲欧美激情综合另类| 精品国产三级普通话版| 欧美一级a爱片免费观看看| 蜜桃久久精品国产亚洲av| 午夜免费观看网址| 一区二区三区国产精品乱码| 男女之事视频高清在线观看| 全区人妻精品视频| 黄色片一级片一级黄色片| 少妇人妻一区二区三区视频| 亚洲成a人片在线一区二区| 亚洲av成人不卡在线观看播放网| 国产av在哪里看| 深爱激情五月婷婷| 精品国内亚洲2022精品成人| 国产精品久久久久久精品电影| 99久久综合精品五月天人人| 国产av在哪里看| 51国产日韩欧美| 每晚都被弄得嗷嗷叫到高潮| 桃红色精品国产亚洲av| 欧美日韩综合久久久久久 | 最近最新中文字幕大全电影3| 一区福利在线观看| 亚洲av不卡在线观看| 亚洲精品日韩av片在线观看 | 国产精品一区二区三区四区免费观看 | 日韩欧美国产一区二区入口| 国产精品美女特级片免费视频播放器| 最好的美女福利视频网| 国产高清videossex| 国产一区二区亚洲精品在线观看| 一个人看视频在线观看www免费 | 人人妻人人澡欧美一区二区| 在线播放国产精品三级| 有码 亚洲区| 国产探花极品一区二区| 久久性视频一级片| 亚洲乱码一区二区免费版| 亚洲精品乱码久久久v下载方式 | 手机成人av网站| 99在线视频只有这里精品首页| 国产午夜精品久久久久久一区二区三区 | 99久久九九国产精品国产免费| 午夜精品在线福利| 日本与韩国留学比较| 级片在线观看| 舔av片在线| 黄色成人免费大全| 美女 人体艺术 gogo| 欧美又色又爽又黄视频| 日本黄色片子视频| 亚洲一区高清亚洲精品| 欧美区成人在线视频| 舔av片在线| 中文字幕高清在线视频| 露出奶头的视频| 一级毛片女人18水好多| 久久久久久久午夜电影| 好看av亚洲va欧美ⅴa在| av视频在线观看入口| 757午夜福利合集在线观看| 一个人免费在线观看电影| 搡老熟女国产l中国老女人| 在线观看美女被高潮喷水网站 | 亚洲精品一卡2卡三卡4卡5卡| 成人av在线播放网站| 中文字幕人成人乱码亚洲影| 啦啦啦韩国在线观看视频| 岛国在线免费视频观看| 91久久精品电影网| 法律面前人人平等表现在哪些方面| 久久香蕉国产精品| 网址你懂的国产日韩在线| 亚洲成人久久爱视频| 日本 av在线| 麻豆一二三区av精品| 国产精品 国内视频| ponron亚洲| 一卡2卡三卡四卡精品乱码亚洲| 变态另类丝袜制服| 欧美绝顶高潮抽搐喷水| 十八禁人妻一区二区| 亚洲乱码一区二区免费版| 国产色爽女视频免费观看| 超碰av人人做人人爽久久 | 久久天躁狠狠躁夜夜2o2o| 日韩欧美国产一区二区入口| 一本精品99久久精品77| 日韩欧美免费精品| 不卡一级毛片| 熟女少妇亚洲综合色aaa.| 欧美一区二区国产精品久久精品| 亚洲专区国产一区二区| 桃色一区二区三区在线观看| 88av欧美| 又黄又爽又免费观看的视频| 狂野欧美白嫩少妇大欣赏| 日韩欧美一区二区三区在线观看| 在线观看舔阴道视频| 久久性视频一级片| 国产真人三级小视频在线观看| 国产激情欧美一区二区| 中亚洲国语对白在线视频| 此物有八面人人有两片| 黄色丝袜av网址大全| 亚洲无线在线观看| 午夜两性在线视频| 久久草成人影院| 99在线人妻在线中文字幕| 网址你懂的国产日韩在线| 国产中年淑女户外野战色| 色吧在线观看| 又黄又粗又硬又大视频| 麻豆国产av国片精品| 一区福利在线观看| 亚洲午夜理论影院| 成年人黄色毛片网站| 亚洲精品成人久久久久久| 99久久精品一区二区三区| 中国美女看黄片| 国产午夜精品久久久久久一区二区三区 | 在线天堂最新版资源| 亚洲精品456在线播放app | 日本精品一区二区三区蜜桃| 日韩精品中文字幕看吧| 国产中年淑女户外野战色| 99热只有精品国产| 少妇高潮的动态图| 三级毛片av免费| 99热只有精品国产| 婷婷亚洲欧美| 国产日本99.免费观看| 中文字幕人妻丝袜一区二区| 老司机在亚洲福利影院| e午夜精品久久久久久久| 久久精品影院6| 亚洲av五月六月丁香网| 亚洲成a人片在线一区二区| 亚洲av成人av| 亚洲专区国产一区二区| 亚洲精品久久国产高清桃花| 日韩成人在线观看一区二区三区| 日韩大尺度精品在线看网址| 日韩av在线大香蕉| 久久欧美精品欧美久久欧美| 精品久久久久久久久久免费视频| 日韩av在线大香蕉| 最好的美女福利视频网| 成人一区二区视频在线观看| 在线观看午夜福利视频| 18禁黄网站禁片免费观看直播| 中文字幕高清在线视频| 亚洲最大成人手机在线| 日本 av在线| 国产精品亚洲美女久久久| 熟女人妻精品中文字幕|