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

    發(fā)動機冷卻風扇性能的優(yōu)化設計研究*

    2015-03-02 06:25:18王天宇劉會霞馮士琳蘇占龍
    機電工程 2015年6期
    關鍵詞:動壓靜壓風扇

    王天宇,王 霄 ,劉會霞,馮士琳,蘇占龍

    (江蘇大學機械工程學院,江蘇鎮(zhèn)江212013)

    0 引 言

    隨著人們對汽車動力性能的要求越來越高,發(fā)動機的工作轉速也在不斷提高,這對汽車冷卻系統(tǒng)提出了更高的要求。因此,在提高發(fā)動機冷卻風扇的整體性能同時縮短設計周期和降低設計成本,是風扇設計的一個重要目標。

    隨著計算流體力學(computational fluid dynamics,CFD)與計算機技術的發(fā)展,數值模擬已經作為一種強大的工具廣泛應用于風扇設計中。數值模擬與試驗研究相比,優(yōu)點在于模型得到簡化,條件容易控制,修改方便,并且具有很好的可重復性,所以設計和研究的成本和所需時間都比試驗研究大大降低。因此,數值模擬已經成為風扇研究不可或缺的重要工具。

    為了提高風扇性能,國內外學者對風扇葉片的翼型[1]、彎掠技術[2]和葉片扭曲技術[3]進行了大量的研究。但多數對風扇靜特性進行對比和研究,未對風扇內部流場進行詳細分析。在風扇優(yōu)化方面,一些研究人員單獨討論不同參數對風扇性能的影響[4],然后綜合選擇最佳參數達到優(yōu)化目的;也有一些學者使用優(yōu)化算法對結構參數進行綜合考慮,如正交試驗法[5]等,對風扇性能進行優(yōu)化。但這些研究多是單目標優(yōu)化,對風扇多個性能指標進行多目標優(yōu)化尚沒有深入研究。

    本研究將參數化建模、CFD 流場分析與多目標優(yōu)化技術相結合,集成在Isight 優(yōu)化平臺中,以提高風扇靜壓和動壓性能為優(yōu)化目標,效率和軸向尺寸為約束條件,對風扇進行多目標優(yōu)化設計。

    1 風扇參數化造型

    風扇的幾何表達方式有:CAD 模型、離散點模型和參數化模型等3 種形式。

    本研究通過參數化造型和擬合手段,將復雜的風扇結構用若干個簡單的控制參數表達出來,利用改變控制參數達到改變整個風扇造型的目的。優(yōu)化中只要通過鎖定部分控制點及參數,僅對幾個重要參數進行調整就能方便、直觀地改變模型。

    1.1 設計參數的確定

    一般來說,研究者沿徑向葉高選取若干個斷面,建立各個斷面的二維翼型,建立一定的積疊線,通過坐標轉換形成三維翼型,將各斷面曲線進行蒙皮操作,得到三維葉片[6]。

    影響風扇氣動性能的參數主要有:葉片數Z、葉輪直徑D、輪轂比Xb、翼型、弦長l、各斷面翼型安裝角β等[7-8]。在優(yōu)化設計之前需選定一些合理的參數。本研究以現有的產品風扇模型為基礎,由于客觀條件限定,葉片數、輪轂比及翼型參數等參數不作修改。風扇葉片的翼型安裝角對風扇性能有十分重要的影響。翼型速度三角形如圖1 所示。軸向速度va沿葉高保持不變時,氣流流經葉片的合速度vm與旋轉平面的夾角ψ 從葉根到葉尖逐漸減小。所以,當葉片翼型的安裝角β 一定時,氣流攻角α=β-ψ,則從葉根到葉尖逐漸增大。所以為使葉片各個斷面均在最有利的攻角下工作,需要對葉片進行扭曲[9]。另一方面,在一定范圍內增加安裝角可以使風扇流量和全壓增加,從而提高風扇冷卻能力。所以研究者通過對葉片安裝角進行優(yōu)化可以達到優(yōu)化風扇性能的目的。

    圖1 翼型速度三角形

    1.2 安裝角的參數化描述

    本研究沿葉高方向選擇12 個斷面的安裝角βi為設計參數,對其進行優(yōu)化設計,風扇的參數化模型如圖2 所示。設計變量過多會影響模型修改和優(yōu)化的速度和效率。因此需要對各斷面的翼型安裝角進行描述,盡量用較少的參數控制斷面翼型安裝角的變化。原模型沿葉高的12 個斷面翼型安裝角為22.977°、21.107°、19.553°、18.300°、17.417°、16.864°、16.869°、17.618°、19.211°、20.396°、20.618°和21.000°。

    圖2 風扇模型

    由于風扇廠家提供的模型葉頂斷面的翼型弦長和安裝角與其他斷面變化規(guī)律不同,葉頂斷面的翼型安裝角被單獨作為一個參數進行優(yōu)化,記為βt。本研究利用Matlab 對其余斷面的翼型安裝角進行二次多項式擬合,得到安裝角β 關于r 的關系曲線如下:

    式中:βi—翼型安裝角,r—各斷面的相對葉高值。

    本研究定義該二次多項式系數為a0,a1,a2與βt一起作為優(yōu)化的設計變量,通過控制4 個參數的變化改變翼型安裝角的分布,從而實現風扇模型的參數化。

    2 數值模型與計算方法

    2.1 流場數值模型建立

    本研究利用三維建模軟件CATIA 根據風扇原始模型建立參數化模型。保持各斷面翼型和葉形積疊線不變,僅改變翼型安裝角。原始模型風扇相關參數如表1所示。

    表1 風扇相關參數

    本研究將模型保存為model 模型,導入前處理軟件Gambit 中進行網格劃分。風扇流場區(qū)域劃分為進口區(qū)、出口區(qū)、過渡區(qū)和旋轉流體區(qū)。本研究綜合考慮試驗設備情況,將進口區(qū)域長度設置為2 000 mm,出口區(qū)域長度設置為3 000 mm。過渡區(qū)域是連接旋轉流體區(qū)和進出口區(qū)域的部分。旋轉流體區(qū)包含風扇,是流場中與風扇一起旋轉起來的區(qū)域。風扇流場區(qū)域模型如圖3 所示。

    圖3 風扇流場區(qū)域模型簡圖(單位:mm)

    模型旋轉軸為x 軸,位于圖3 中心線上,坐標原點位于風扇輪轂中心,x 軸正方向指向出口。本研究在Gambit 中對流場模型進行網格劃分。對于穩(wěn)態(tài)求解問題,需進行網格無關性驗證。網格劃分方案與模擬結果如表2 所示,筆者采取4 種方案對流場進行網格劃分和數值模擬,將模擬結果進行對比,從模型中提取轉速為2 800 r/min、入口流量為0.921 779 kg/s 條件下的風扇靜壓和靜壓效率。結果表明,1.99 ×106網格時風扇靜壓和靜壓效率與試驗結果(靜壓159.6 Pa,靜壓效率31.6%)相差較大;當網格數量從2.8 ×106增大到4.1×106,隨著網格數量增加風扇靜壓和靜壓效率變化不大,可以認為2.8 ×106的網格數量已達到網格無關。2.8×106網格劃分方案具體網格尺寸為:扇葉和頂圈表面為2 mm的tri 面網格,旋轉流體區(qū)其他表面為5 mm 的tri 面網格,旋轉流體區(qū)為3 mm 的Tet/Hybrid 體網格,過渡區(qū)為6 mm的Tet/Hybrid 體網格,進、出口區(qū)為20 mm的Hex/Wedge 體網格。流場網格劃分結果如圖4 所示。

    表2 網格劃分方案與模擬結果

    圖4 流場網格

    2.2 邊界條件與計算方法

    本研究采用質量流量進口(mass flow rate)和壓力出口(pressure outlet),設置出口相對壓力為0,過渡區(qū)與進、出口區(qū)的交界面為interior 邊界。風扇表面為wall,設置旋轉表面,相對速度為0,其他壁面設置為無滑移固壁條件。旋轉流體區(qū)使用動參考系(moving reference frame,MRF),使旋轉流體區(qū)轉動起來,轉速為2 800 r/min。過渡區(qū)和進、出口區(qū)則設置為靜止流體。

    筆者在研究風扇的內部流場時采用定常模擬,選擇RNG k-ε 模型進行湍流模擬。使用SIMPLE 算法進行壓力-速度耦合,控制方程使用二階迎風格式進行離散。進行湍流模擬時進出口邊界條件需要給出湍流參數,本研究選擇給定湍流強度I 和水力直徑DH(intensity and hydraulic diameter)。由于進出口區(qū)域截面為非圓管,故使用等效水力直徑,如下式計算:

    式中:DH—水力直徑,A—浸潤截面面積,P—橫截面濕周周長。

    計算湍流強度I 需要由等效水力直徑計算出雷諾數Re,計算公式如下:

    式中:Re—雷諾數,v—空氣相對流速,μ—空氣動力粘度。

    湍流強度I 計算公式為:

    3 基于Isight 平臺的風扇參數化優(yōu)化及分析

    Isight 是廣泛應用于多領域的功能強大的多學科設計優(yōu)化平臺。該軟件可以通過集成和驅動其他軟件實現和管理復雜的仿真過程,運用多種優(yōu)化算法尋得優(yōu)化方案,從而達到縮短設計周期、提高設計效率和降低研發(fā)成本的目的。Isight 能夠集成廣泛的商業(yè)CAD/CAE 軟件,可以快速建立復雜的仿真優(yōu)化流程。仿真流程可以實現自動化和可視化,方便設計人員控制和修改優(yōu)化過程。

    3.1 試驗設計與近似模型

    試驗設計(design of experiment,DOE)是一種應用統(tǒng)計學思想來處理變異從而達到改進產品質量和工藝的方法。本研究試驗設計采用最優(yōu)拉丁超立方設計方法,使所有的試驗點盡量均勻地分布在設計空間,具有非常好的空間填補性和均衡性。4 個設計變量a0,a1,a2,βt的設計空間分別為:13~33,0~5,0~0.5,15°~28°。本研究由最優(yōu)拉丁超立方法得到15 組不同的參數樣本點,以供建立輸出響應的近似模型。

    近似模型方法(approximation models)是通過已知點的輸入變量與響應值預測未知點的響應值的方法。常用的近似模型有:響應面近似模型、克里格近似模型、徑向基神經網絡近似模型、正交多項式近似模型等。由于克里格近似模型無論擬合低階或者是高階非線性都有著廣泛的應用,本研究利用試驗設計得到的樣本數據建立克里格近似模型(kriging)。該模型用于后續(xù)的優(yōu)化設計中,可以大幅度提高優(yōu)化效率。

    3.2 優(yōu)化算法

    近年來,多目標遺傳算法在實際工程優(yōu)化問題的解決中得到越來越多的應用。多目標遺傳算法是一種能夠不斷優(yōu)化特定問題無支配前沿的優(yōu)化算法,非支配排序遺傳算法(non-donminated sorting genetic algorithms,NSGA Ⅱ)是一種典型的多目標遺傳算法[10]。NSGAⅡ是由K.Ded 等于2000年在原有的NSGA 算法基礎上提出的改良版。NSGA 算法,缺少保護最優(yōu)個體的策略,計算復雜度較大,且需要設置共享參數,而NSGAⅡ減少了復雜度,提出的擁擠度算子無須參數指定,可以保存最優(yōu)個體[11]。擁擠度算子的提出不僅保持了種群的多樣性,也使得優(yōu)化結果在目標空間中均勻分布。

    3.3 優(yōu)化模型

    風扇優(yōu)化目的就是尋求一定的結構參數組合,使其達到最佳的氣動性能。靜壓是評價風扇性能的重要參數之一,它代表了風扇推動空氣在流道中流動克服阻力的能力。靜壓效率體現了風扇內部流動品質的好壞,在軸功率一定的情況下,靜壓效率越高,風扇能量損失越少,冷卻性能越好。靜壓效率正比于靜壓,因此,只需優(yōu)化靜壓就可以達到優(yōu)化靜壓效率的目的。風扇動壓代表氣體對葉片的氣動力,該氣動力產生的轉矩為風扇旋轉提供動力,這樣將降低風扇消耗的功率,節(jié)省燃油[12]。

    建立多目標優(yōu)化模型有:

    式中:Ps—風扇靜壓;Pd—風扇動壓。

    為防止出現偏差較大的子代,還需對靜壓效率進行約束,同時考慮到發(fā)動機艙的空間,也需控制風扇的軸向尺寸,故約束條件如下:

    (1)靜壓效率約束:ηmin≤η≤ηmax;

    (2)軸向尺寸約束:Hmin≤H≤Hmax。

    其中:η—靜壓效率,H—軸向尺寸。

    設計變量:x={a0,a1,a2,βt}。

    4 分析結果與試驗結果

    4.1 優(yōu)化后氣動性能與原型對比

    進口流量取0.921 779 kg/s,遺傳迭代次數為20代,種群大小為15,交叉概率取0.9,變異分布指數取10。優(yōu)化后Pareto 解個數為10 個,綜合考慮效率及尺寸問題,選擇最優(yōu)解如表3 所示。

    表3 優(yōu)化結果

    風扇優(yōu)化前和優(yōu)化后的葉片參數整理如表4 所示。優(yōu)化后的風扇葉片葉根和葉尖安裝角有所增大,而中部安裝角增大不多,增加了葉片扭曲。

    表4 優(yōu)化前、后參數對比

    試驗結果與優(yōu)化前、后模擬結果性能參數對比如表5 所示。

    表5 試驗結果與優(yōu)化前、后模擬結果性能參數對比

    由表5 可以看出,優(yōu)化前的數值模擬結果與試驗結果誤差基本不大于5%,說明數值模擬所設置的邊界條件是合理的,流場分析具有一定的可靠性。優(yōu)化后風扇靜壓為175.923 1 Pa,比優(yōu)化前提高了12.840 6%,說明風扇推動空氣克服流道阻力的能力有所提高,使得風扇冷卻能力提高。風扇動壓有一定增加,說明優(yōu)化后可以減少發(fā)動機的功率消耗。而且風扇靜壓效率提高到了38.291 3%。

    4.2 壓力場計算結果分析

    風扇吸力面靜壓云圖如圖5 所示。優(yōu)化前葉片低壓區(qū)主要集中在葉片上部1/3 區(qū)域,壓力由葉片前緣向葉片尾緣逐漸升高。優(yōu)化后葉片低壓區(qū)由上部1/3區(qū)域擴展到葉片2/3 區(qū)域,靜壓沿葉高方向分布更加均勻,葉片前緣到葉片后緣壓力梯度減小。優(yōu)化后風扇進口靜壓最大值為- 192. 972 52 Pa,最小值為-355.741 91 Pa,靜壓最大值比優(yōu)化前降低18.470 11 Pa,靜壓最小值比優(yōu)化前降低了6.155 3 Pa。進口平均靜壓為-246.955 75 Pa,比優(yōu)化前的-222.159 21 Pa 有所降低。

    圖5 葉片吸力面靜壓云圖

    風扇壓力面靜壓云圖如圖6 所示。由圖6 可知,靜壓最大值出現在葉片前緣中部區(qū)域,優(yōu)化前高壓區(qū)較小,壓力梯度大。優(yōu)化后,高壓區(qū)域明顯增大,沿葉高方向分布變廣,且向尾緣擴展。出口壓力的高壓區(qū)分布更廣,邊緣低壓區(qū)域明顯減少。優(yōu)化后出口的靜壓平均值為-28.165 892 Pa,比優(yōu)化前的平均靜壓值所提高,最大值為21.405 825 Pa,最小值為-242.646 8 Pa,優(yōu)化后壓力范圍增大。

    圖6 葉片壓力面靜壓云圖

    優(yōu)化后風扇壓力面靜壓提高,而吸力面靜壓降低,故風扇靜壓增大,還增大了風扇靜壓分布范圍,增強了風扇推動空氣克服流道阻力的能力。如此,風扇可以使空氣順利通過散熱器將更多熱量及時散發(fā)出去,提高風扇的靜壓效率,從而提高冷卻效率。

    風扇子午面動壓云圖如圖7 所示。圖7 表明,優(yōu)化后風扇葉根處高動壓區(qū)明顯減小,輪轂附近的回流減少,使風扇功率損失減少。由于風扇進出口靜壓差增大,頂圈出口處動壓增大。葉片壓力面附近動壓梯度減小,葉高方向氣壓更加穩(wěn)定。風扇進口動壓平均值變化不大,出口動壓平均值升高到84.171 608 Pa,對葉片提供的氣動力增大,減少了風扇消耗的功率。

    圖7 風扇子午面動壓云圖

    4.3 速度場計算結果分析

    子午面的速度分布圖如圖8 所示。葉高方向上氣流軸向速度有所提高,且軸向速度變化更小,氣流比優(yōu)化前更加平穩(wěn)。葉根處的速度回流明顯減小,頂圈進口處回流減小而中部回流增大,出口處回流速度增大,氣流能量因為回流產生的損失減小。風扇入口、出口速度云圖如圖9、圖10 所示,氣流進入風扇由于扇葉轉動產生速度變化,從進口開始外圈氣流速度開始增大,隨著軸向推進,葉片中上部附近速度越來越大,葉片之間的低速區(qū)速度逐漸增大,低速區(qū)逐漸減小。優(yōu)化后葉片之間的低速區(qū)減小,葉片之間的速度更加均勻,減少了氣流能量的摻混,可以提高風扇的通風量。入口平均速度為7.731 21 m/s,出口速度平均速度為10.307 855 m/s,均比優(yōu)化前模型有所提高。

    圖8 風扇子午面速度分布圖

    圖9 風扇入口速度云圖

    圖10 風扇出口速度云圖

    5 結束語

    本研究對發(fā)動機風扇進行了多目標優(yōu)化。

    (1)經優(yōu)化前、后靜態(tài)性能參數對比顯示,改變安裝角葉高方向分布規(guī)律能夠達到改善風扇性能的目的,通過利用擬合方法減少優(yōu)化參數可以提高優(yōu)化效率,利用Isight 進行優(yōu)化是有效可行的。

    (2)靜壓計算結果表明,優(yōu)化前模型壓力沿葉高方向分布不夠均勻,壓力面做功區(qū)域較小;優(yōu)化后模型吸力面低壓區(qū)域沿葉高方向擴展,壓力面高壓區(qū)域增大,吸力面和壓力面壓差增大,風扇靜壓有所提高。動壓計算結果表明,優(yōu)化后葉根動壓減小,頂圈處動壓增大,葉片附近動壓梯度減小,平均動壓增大,提供給葉片的氣動力增大。

    (3)速度計算結果表明,葉根回流減小明顯,頂圈附近進出口回流增大,中部回流減小。優(yōu)化后風扇進出口速度有所提高,葉片間的速度更加均勻,低速區(qū)范圍減小,提高風扇的通風量。

    [1]張 旭,李 偉,邢靜忠.相對彎度對鈍尾緣改型提升翼型氣動性能的影響[J].農業(yè)機械學報,2014,45(8):214-219.

    [2]HURAULT J,KOUIDRI S,BAKIR F. Experimental investigations on the wall pressure measurement on the blade of axial flow fans[J]. Experimental Thermal and Fluid Science,2012(40):29-37.

    [3]田 夫,袁國凱,趙真真,等.單葉輪軸流通風機的葉片扭曲規(guī)律研究及薄葉型設計[J]. 風機技術,2015(5):49-54.

    [4]習 羽.工程車輛冷卻風扇流體特性研究[D]. 長春:吉林大學機械科學與工程學院,2013.

    [5]唐 釗.發(fā)動機冷卻風扇葉片參數的研究和優(yōu)化[D].廣州:華南理工大學汽車工程學院,2012.

    [6]鐘守山.發(fā)動機冷卻風扇造型設計與性能計算方法的研究[D].廣州:華南理工大學汽車工程學院,2011.

    [7]湯黎明. 工程機械冷卻風扇流場特性與氣動噪聲研究[D].長春:吉林大學機械學院,2014.

    [8]霍進剛.風扇冷卻系統(tǒng)在閉式回路中的應用[J].液壓氣動與密封,2011(10):52-55.

    [9]吳秉禮,高延福.空氣冷卻軸流通風機[M].哈爾濱:哈爾濱工業(yè)大學出版社,2006.

    [10]DEHGHANIAN P,HOSSEINI S H,MOEIN,et al. Optimal siting of DG units in power systems from a probabilistic multi-objective optimizati-on perspective[J]. Electrical Power and Energy Systems,2013(51):14-26.

    [11]CHEN Zhi-huan,YUAN Xiao-hui,JI Bin,et al. Design of a fractional order PID controller for hydraulic turbine regulating system using chaotic non-dominated sorting genetic algorithm II[J]. Energy Conversion and Management,2014(84):390-404.

    [12]宮春梅.基于逆向工程的汽車冷卻風扇葉片的研究與改型設計[D].秦皇島:燕山大學機械工程學院,2012.

    猜你喜歡
    動壓靜壓風扇
    靜壓法沉樁對周邊環(huán)境影響及質量控制
    國內首個現代箔片氣體動壓軸承技術培訓班在長沙成功舉辦
    風機技術(2021年3期)2021-08-05 07:41:36
    靜壓托換樁在某濕陷性黃土場地地基加固中的應用
    超精密液體靜壓轉臺裝配技術
    一種基于空氣靜壓支承的自調心裝置
    電風扇
    故事大王(2017年4期)2017-05-08 07:53:40
    基于智能手機控制風扇運行的實現
    新蒙迪歐車冷卻風扇常高速運轉
    南屯煤礦深部泵房硐室群動壓失穩(wěn)機理及控制對策
    煤炭學報(2015年10期)2015-12-21 01:55:27
    強烈動壓巷道支護技術探討
    亚洲高清免费不卡视频| 少妇被粗大猛烈的视频| 久久久欧美国产精品| 麻豆乱淫一区二区| 国产成人福利小说| 人妻夜夜爽99麻豆av| 国产在线一区二区三区精| 激情五月婷婷亚洲| 成人毛片60女人毛片免费| 成人美女网站在线观看视频| 国产免费视频播放在线视频| 国产亚洲5aaaaa淫片| 国产精品一二三区在线看| 日韩成人av中文字幕在线观看| 欧美精品一区二区大全| 亚洲美女搞黄在线观看| 国产精品嫩草影院av在线观看| 日本猛色少妇xxxxx猛交久久| 中文字幕av成人在线电影| 亚洲美女搞黄在线观看| 偷拍熟女少妇极品色| 欧美日韩精品成人综合77777| 亚洲精品一二三| 久久精品国产亚洲网站| av一本久久久久| 亚洲精品视频女| 成年女人在线观看亚洲视频 | 婷婷色综合www| 久久国产乱子免费精品| 亚洲国产精品999| 春色校园在线视频观看| 一级毛片 在线播放| 亚洲精品aⅴ在线观看| 男女那种视频在线观看| 久久精品夜色国产| 一级二级三级毛片免费看| 免费看a级黄色片| 国产探花极品一区二区| 亚洲天堂国产精品一区在线| 亚洲精品日本国产第一区| a级毛片免费高清观看在线播放| 国产成人91sexporn| 自拍偷自拍亚洲精品老妇| 三级经典国产精品| 国产一区有黄有色的免费视频| 日本黄色片子视频| 又爽又黄无遮挡网站| 一二三四中文在线观看免费高清| 亚洲精品成人久久久久久| 久久久久久久久久久免费av| 精品久久久噜噜| 欧美成人精品欧美一级黄| 亚洲国产精品成人综合色| 草草在线视频免费看| 1000部很黄的大片| 黄色一级大片看看| 97在线人人人人妻| 黄色视频在线播放观看不卡| 一级毛片久久久久久久久女| 国产欧美另类精品又又久久亚洲欧美| 久久热精品热| 久久精品国产亚洲av涩爱| 国产高清三级在线| 91精品国产九色| 禁无遮挡网站| 欧美国产精品一级二级三级 | 欧美精品国产亚洲| 亚洲欧美日韩无卡精品| 黄色配什么色好看| 涩涩av久久男人的天堂| av国产久精品久网站免费入址| 国产精品一及| 91aial.com中文字幕在线观看| 99久久九九国产精品国产免费| 夜夜爽夜夜爽视频| 日韩亚洲欧美综合| 国精品久久久久久国模美| 国产一区亚洲一区在线观看| 日本-黄色视频高清免费观看| 国产高清三级在线| 丝袜脚勾引网站| 啦啦啦在线观看免费高清www| 一级毛片黄色毛片免费观看视频| 国产在线男女| 99热这里只有精品一区| 一本一本综合久久| 国产亚洲91精品色在线| 最近中文字幕2019免费版| 一级毛片aaaaaa免费看小| 国产探花在线观看一区二区| 直男gayav资源| 久久精品国产亚洲网站| 91久久精品国产一区二区三区| 日本午夜av视频| 午夜福利网站1000一区二区三区| 国产午夜精品一二区理论片| 成年免费大片在线观看| 毛片一级片免费看久久久久| 亚洲精品国产成人久久av| 婷婷色综合大香蕉| 国产一区二区在线观看日韩| av免费观看日本| 日韩伦理黄色片| 少妇人妻久久综合中文| 久久国产乱子免费精品| 深夜a级毛片| 午夜日本视频在线| 国产黄色免费在线视频| 免费播放大片免费观看视频在线观看| 99久久中文字幕三级久久日本| 午夜激情久久久久久久| 欧美性感艳星| 欧美xxxx黑人xx丫x性爽| 日本熟妇午夜| 天天一区二区日本电影三级| 97超碰精品成人国产| 99re6热这里在线精品视频| 亚洲国产高清在线一区二区三| 男女边摸边吃奶| 成人毛片a级毛片在线播放| 欧美精品一区二区大全| 国产一区二区在线观看日韩| 国产欧美另类精品又又久久亚洲欧美| 国产国拍精品亚洲av在线观看| 国产黄a三级三级三级人| 国精品久久久久久国模美| 久久久久九九精品影院| 精品一区二区三区视频在线| 色视频www国产| 日日啪夜夜爽| 亚洲在线观看片| 男女那种视频在线观看| 91狼人影院| 久久久久久久久久人人人人人人| 亚洲,欧美,日韩| 国产成人午夜福利电影在线观看| 赤兔流量卡办理| 欧美日本视频| 国产成人午夜福利电影在线观看| 午夜免费男女啪啪视频观看| 亚洲人成网站在线观看播放| 中文资源天堂在线| 欧美老熟妇乱子伦牲交| 国产伦理片在线播放av一区| 国产精品人妻久久久久久| 久久久精品免费免费高清| 一级毛片aaaaaa免费看小| 成人一区二区视频在线观看| 午夜免费观看性视频| 18+在线观看网站| 亚洲av中文字字幕乱码综合| 国产在线男女| 欧美+日韩+精品| 99久久九九国产精品国产免费| 国产精品.久久久| 最近中文字幕2019免费版| 国产av码专区亚洲av| 国产有黄有色有爽视频| 人妻夜夜爽99麻豆av| 日韩av不卡免费在线播放| 国产黄色免费在线视频| 国产一区二区三区av在线| 国产黄频视频在线观看| 国产精品成人在线| 日韩中字成人| 国内揄拍国产精品人妻在线| 亚洲伊人久久精品综合| 国产成人免费观看mmmm| 国产一区亚洲一区在线观看| 黑人高潮一二区| 久久久国产一区二区| 97在线视频观看| 亚洲最大成人av| 欧美zozozo另类| 91精品伊人久久大香线蕉| av卡一久久| 日本一二三区视频观看| 亚洲欧美中文字幕日韩二区| 狂野欧美白嫩少妇大欣赏| 久久99热6这里只有精品| 欧美 日韩 精品 国产| 中文字幕制服av| 乱码一卡2卡4卡精品| 亚洲在线观看片| 一级毛片aaaaaa免费看小| 国产伦在线观看视频一区| 白带黄色成豆腐渣| 精品国产露脸久久av麻豆| 午夜福利在线观看免费完整高清在| 久久热精品热| 91狼人影院| 国产精品一区二区在线观看99| 麻豆久久精品国产亚洲av| 亚洲av中文av极速乱| 亚洲成人中文字幕在线播放| 一区二区av电影网| 各种免费的搞黄视频| 国产成人福利小说| 在线观看国产h片| 日本爱情动作片www.在线观看| 在线观看免费高清a一片| 毛片女人毛片| 麻豆乱淫一区二区| 一级av片app| 高清在线视频一区二区三区| 国产成人a区在线观看| 看黄色毛片网站| 十八禁网站网址无遮挡 | 午夜福利高清视频| 中文欧美无线码| 亚洲av男天堂| 国产亚洲最大av| 综合色丁香网| 久久久久九九精品影院| 免费黄色在线免费观看| 国产欧美日韩一区二区三区在线 | 综合色av麻豆| 国产免费视频播放在线视频| 日韩亚洲欧美综合| 九九爱精品视频在线观看| 成人午夜精彩视频在线观看| 黄色日韩在线| 国产 一区精品| 国国产精品蜜臀av免费| 视频中文字幕在线观看| 午夜爱爱视频在线播放| 精品99又大又爽又粗少妇毛片| av专区在线播放| 黄色日韩在线| 自拍偷自拍亚洲精品老妇| 亚洲自偷自拍三级| 免费大片18禁| 成年人午夜在线观看视频| 伊人久久国产一区二区| 国产一区二区三区综合在线观看 | 一区二区三区乱码不卡18| 国产乱来视频区| 国产成人精品一,二区| 人妻一区二区av| 亚洲国产高清在线一区二区三| 免费av不卡在线播放| 国内精品美女久久久久久| 亚洲成人中文字幕在线播放| 自拍欧美九色日韩亚洲蝌蚪91 | 日本熟妇午夜| 三级国产精品欧美在线观看| 菩萨蛮人人尽说江南好唐韦庄| 91久久精品国产一区二区三区| 婷婷色综合www| 免费看日本二区| 99久久精品一区二区三区| 亚洲经典国产精华液单| 国产极品天堂在线| 色吧在线观看| 18禁在线播放成人免费| 日韩亚洲欧美综合| 久久久久久久国产电影| 人妻一区二区av| 午夜福利网站1000一区二区三区| 美女国产视频在线观看| 18禁在线播放成人免费| 下体分泌物呈黄色| 免费观看的影片在线观看| 99热这里只有是精品在线观看| 伊人久久国产一区二区| 国产精品人妻久久久久久| 人体艺术视频欧美日本| 2022亚洲国产成人精品| 最近最新中文字幕免费大全7| 嫩草影院精品99| 亚洲欧美清纯卡通| 干丝袜人妻中文字幕| 欧美97在线视频| 黄色日韩在线| 日韩欧美精品v在线| 又粗又硬又长又爽又黄的视频| 校园人妻丝袜中文字幕| 在现免费观看毛片| 欧美高清性xxxxhd video| 亚洲无线观看免费| 天堂中文最新版在线下载 | av又黄又爽大尺度在线免费看| 偷拍熟女少妇极品色| 在现免费观看毛片| 免费黄网站久久成人精品| 中国美白少妇内射xxxbb| 麻豆久久精品国产亚洲av| 又大又黄又爽视频免费| 国产一区二区亚洲精品在线观看| 精品久久久久久久末码| 男人狂女人下面高潮的视频| 成年女人在线观看亚洲视频 | 亚洲av一区综合| 午夜爱爱视频在线播放| 亚洲天堂国产精品一区在线| 日韩国内少妇激情av| 男女啪啪激烈高潮av片| 中国三级夫妇交换| 日日摸夜夜添夜夜爱| 秋霞在线观看毛片| 国产色爽女视频免费观看| 国产白丝娇喘喷水9色精品| 中文天堂在线官网| 18禁裸乳无遮挡动漫免费视频 | 在线播放无遮挡| 免费观看av网站的网址| 久久精品国产自在天天线| 成人漫画全彩无遮挡| 简卡轻食公司| 男女那种视频在线观看| 国产黄片美女视频| 日韩欧美精品v在线| 激情五月婷婷亚洲| 亚洲国产色片| 尤物成人国产欧美一区二区三区| 伊人久久精品亚洲午夜| 欧美亚洲 丝袜 人妻 在线| 国产成人福利小说| 久久久久九九精品影院| 亚洲不卡免费看| 十八禁网站网址无遮挡 | 国产精品伦人一区二区| 国产黄片美女视频| 夜夜看夜夜爽夜夜摸| 亚洲怡红院男人天堂| 丰满少妇做爰视频| 国产高清三级在线| 22中文网久久字幕| 日本猛色少妇xxxxx猛交久久| 简卡轻食公司| 亚洲三级黄色毛片| 两个人的视频大全免费| 日韩av免费高清视频| 亚洲精品日本国产第一区| 国产日韩欧美在线精品| 纵有疾风起免费观看全集完整版| 免费av毛片视频| 久久久久九九精品影院| videos熟女内射| 美女脱内裤让男人舔精品视频| 亚洲国产av新网站| 好男人视频免费观看在线| 蜜桃久久精品国产亚洲av| 亚洲色图av天堂| 视频区图区小说| 嫩草影院入口| 在线观看一区二区三区激情| 亚洲av日韩在线播放| 美女xxoo啪啪120秒动态图| 国产精品99久久99久久久不卡 | 欧美日韩在线观看h| 亚洲精品国产成人久久av| 少妇人妻精品综合一区二区| 麻豆成人av视频| 2021天堂中文幕一二区在线观| 亚洲精品国产色婷婷电影| 最近的中文字幕免费完整| 日日摸夜夜添夜夜添av毛片| 国产黄色免费在线视频| 一级二级三级毛片免费看| 国产精品久久久久久久久免| 亚洲无线观看免费| 日本色播在线视频| 一区二区三区四区激情视频| 久久久久久久国产电影| 女人久久www免费人成看片| 国内揄拍国产精品人妻在线| 好男人视频免费观看在线| 久久精品国产亚洲av天美| 高清午夜精品一区二区三区| 菩萨蛮人人尽说江南好唐韦庄| 日本欧美国产在线视频| 国产欧美日韩精品一区二区| 亚洲欧美日韩无卡精品| 一区二区三区免费毛片| 国产综合懂色| 国产午夜福利久久久久久| 精品久久久久久久久av| 五月天丁香电影| 久久亚洲国产成人精品v| 久久久久精品久久久久真实原创| 一个人观看的视频www高清免费观看| 欧美成人午夜免费资源| 国产男女超爽视频在线观看| 一本色道久久久久久精品综合| av女优亚洲男人天堂| 亚洲精华国产精华液的使用体验| 欧美亚洲 丝袜 人妻 在线| 在线观看美女被高潮喷水网站| 亚洲国产成人一精品久久久| 国产高清不卡午夜福利| 一级片'在线观看视频| 日本黄色片子视频| 岛国毛片在线播放| 少妇人妻久久综合中文| 春色校园在线视频观看| 国产午夜精品一二区理论片| 亚洲成人精品中文字幕电影| 日日摸夜夜添夜夜爱| 丝袜喷水一区| 亚洲精品乱码久久久久久按摩| 久久综合国产亚洲精品| 精品99又大又爽又粗少妇毛片| 免费观看av网站的网址| a级毛片免费高清观看在线播放| 欧美丝袜亚洲另类| 国产成人一区二区在线| 国产精品秋霞免费鲁丝片| 国产亚洲精品久久久com| 国产精品国产av在线观看| 国产真实伦视频高清在线观看| av国产久精品久网站免费入址| 赤兔流量卡办理| 日韩中字成人| 久久久久网色| 高清视频免费观看一区二区| 色吧在线观看| 国产美女午夜福利| 亚洲最大成人中文| 少妇人妻久久综合中文| 在现免费观看毛片| 国产精品国产av在线观看| 久久精品国产自在天天线| 日韩免费高清中文字幕av| 18+在线观看网站| 免费av毛片视频| 简卡轻食公司| 国产欧美亚洲国产| 亚洲精品乱久久久久久| 51国产日韩欧美| 亚洲国产精品999| 精品久久久久久久久亚洲| 欧美日韩国产mv在线观看视频 | 国产91av在线免费观看| 欧美日韩一区二区视频在线观看视频在线 | 九九爱精品视频在线观看| 日韩欧美 国产精品| 综合色丁香网| 免费黄频网站在线观看国产| 纵有疾风起免费观看全集完整版| 听说在线观看完整版免费高清| 熟妇人妻不卡中文字幕| 欧美日韩一区二区视频在线观看视频在线 | 黄片无遮挡物在线观看| 好男人视频免费观看在线| 成年女人在线观看亚洲视频 | 国产淫片久久久久久久久| 特级一级黄色大片| 丰满人妻一区二区三区视频av| 亚洲人与动物交配视频| 精品酒店卫生间| 国产欧美另类精品又又久久亚洲欧美| 国产乱来视频区| 黄色日韩在线| 五月玫瑰六月丁香| 中文天堂在线官网| 最近中文字幕高清免费大全6| 真实男女啪啪啪动态图| 午夜福利在线在线| 亚洲精品国产av成人精品| 熟女av电影| av免费在线看不卡| 日韩av在线免费看完整版不卡| 国产精品秋霞免费鲁丝片| 干丝袜人妻中文字幕| 亚洲国产高清在线一区二区三| a级毛片免费高清观看在线播放| 久久久久久久大尺度免费视频| 国产成人精品久久久久久| 在线观看国产h片| 又黄又爽又刺激的免费视频.| 只有这里有精品99| 18禁裸乳无遮挡免费网站照片| 18+在线观看网站| 一级毛片 在线播放| 日本黄色片子视频| 国产av码专区亚洲av| 听说在线观看完整版免费高清| 久久久久久久亚洲中文字幕| 日韩伦理黄色片| 久久99热6这里只有精品| 我的老师免费观看完整版| 视频中文字幕在线观看| 丝袜喷水一区| 国产伦精品一区二区三区视频9| 一级黄片播放器| 交换朋友夫妻互换小说| 嫩草影院入口| 亚洲精华国产精华液的使用体验| 免费av毛片视频| 麻豆精品久久久久久蜜桃| 大话2 男鬼变身卡| 欧美激情在线99| 少妇的逼水好多| 搡女人真爽免费视频火全软件| 色视频www国产| 麻豆久久精品国产亚洲av| 一个人观看的视频www高清免费观看| 亚洲三级黄色毛片| 亚洲欧洲国产日韩| 啦啦啦啦在线视频资源| 永久网站在线| 国精品久久久久久国模美| 午夜免费男女啪啪视频观看| 亚洲精品色激情综合| 亚洲久久久久久中文字幕| 2021天堂中文幕一二区在线观| 精品一区二区三区视频在线| 夫妻性生交免费视频一级片| 69人妻影院| 国产 一区 欧美 日韩| 精华霜和精华液先用哪个| 三级国产精品片| 中文天堂在线官网| 久久久久九九精品影院| 久久久久久久久久久免费av| 在线观看美女被高潮喷水网站| 在线观看一区二区三区激情| 99久久精品热视频| 又爽又黄无遮挡网站| 亚洲一区二区三区欧美精品 | 日韩在线高清观看一区二区三区| 国产亚洲91精品色在线| 大又大粗又爽又黄少妇毛片口| 爱豆传媒免费全集在线观看| 亚洲av在线观看美女高潮| 少妇被粗大猛烈的视频| 性色av一级| av网站免费在线观看视频| av女优亚洲男人天堂| 国产成人精品久久久久久| 亚洲国产最新在线播放| 搞女人的毛片| 我要看日韩黄色一级片| a级一级毛片免费在线观看| 国产成人福利小说| 80岁老熟妇乱子伦牲交| 国产成人精品婷婷| 亚洲国产精品国产精品| 亚洲美女视频黄频| 成人亚洲精品av一区二区| 亚洲成人中文字幕在线播放| 久久精品国产亚洲网站| 日韩一本色道免费dvd| 99久久精品国产国产毛片| 欧美最新免费一区二区三区| 国产精品女同一区二区软件| 人妻 亚洲 视频| 精品久久久久久久久av| 在线观看人妻少妇| 99久久中文字幕三级久久日本| 久久精品国产亚洲av天美| 亚洲精品成人久久久久久| 色视频www国产| 国内少妇人妻偷人精品xxx网站| 国产美女午夜福利| 亚洲精品自拍成人| 亚洲国产色片| 欧美精品人与动牲交sv欧美| 精品一区二区三区视频在线| 免费观看av网站的网址| 美女内射精品一级片tv| 一个人看的www免费观看视频| 免费观看无遮挡的男女| 不卡视频在线观看欧美| av在线天堂中文字幕| kizo精华| 少妇裸体淫交视频免费看高清| 伦精品一区二区三区| av一本久久久久| 午夜视频国产福利| 九色成人免费人妻av| 美女被艹到高潮喷水动态| av又黄又爽大尺度在线免费看| 国产 一区 欧美 日韩| 亚洲国产av新网站| 涩涩av久久男人的天堂| 国产精品国产三级国产av玫瑰| 别揉我奶头 嗯啊视频| 一区二区av电影网| 国产91av在线免费观看| 七月丁香在线播放| 国产成人a区在线观看| 国产男人的电影天堂91| 日韩视频在线欧美| 麻豆国产97在线/欧美| 大又大粗又爽又黄少妇毛片口| 欧美另类一区| 18+在线观看网站| 嫩草影院精品99| 人妻一区二区av| 国产精品精品国产色婷婷| 欧美日本视频| 欧美性感艳星| 精品少妇黑人巨大在线播放| 欧美日韩视频高清一区二区三区二| 十八禁网站网址无遮挡 | 欧美97在线视频| 少妇的逼好多水| 国产欧美日韩精品一区二区| a级毛片免费高清观看在线播放| 看黄色毛片网站| 在线观看国产h片| 成人高潮视频无遮挡免费网站| 18禁在线无遮挡免费观看视频| 亚洲av二区三区四区| 丰满少妇做爰视频| 大话2 男鬼变身卡| 国产人妻一区二区三区在| 日日啪夜夜撸| 韩国av在线不卡| 中文字幕av成人在线电影| 亚洲精品成人av观看孕妇| 久久国产乱子免费精品| 99久久人妻综合|