• <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
    強烈動壓巷道支護技術探討
    青草久久国产| 欧美日韩一级在线毛片| 狂野欧美白嫩少妇大欣赏| 欧美色视频一区免费| 一本一本综合久久| 亚洲第一欧美日韩一区二区三区| 免费在线观看日本一区| 久久精品国产亚洲av涩爱 | 亚洲精品国产精品久久久不卡| 成人特级av手机在线观看| 国产伦人伦偷精品视频| 一本一本综合久久| 天堂影院成人在线观看| 中文资源天堂在线| 国产老妇女一区| 18+在线观看网站| 99国产精品一区二区三区| 亚洲人成网站在线播放欧美日韩| 亚洲精品在线观看二区| 国产视频一区二区在线看| 久久久久久九九精品二区国产| 久久精品国产自在天天线| 深爱激情五月婷婷| 1024手机看黄色片| 久久久国产成人精品二区| 国产亚洲欧美98| 国产精品亚洲av一区麻豆| 国产真实乱freesex| 操出白浆在线播放| 毛片女人毛片| 成年女人毛片免费观看观看9| 亚洲在线观看片| 久久精品亚洲精品国产色婷小说| 国产高清视频在线播放一区| 草草在线视频免费看| av天堂在线播放| 久久国产乱子伦精品免费另类| 啦啦啦免费观看视频1| 香蕉av资源在线| 色综合站精品国产| 国产精品久久久久久人妻精品电影| 欧美成狂野欧美在线观看| 欧美黄色淫秽网站| 国内精品久久久久精免费| 搞女人的毛片| 欧美区成人在线视频| 精品电影一区二区在线| 国产一区二区在线观看日韩 | 在线免费观看的www视频| 久久久国产精品麻豆| 欧美在线一区亚洲| 国产麻豆成人av免费视频| 午夜福利成人在线免费观看| 看片在线看免费视频| 中文资源天堂在线| 岛国在线观看网站| 色综合亚洲欧美另类图片| 国产精品 国内视频| 亚洲人成电影免费在线| 99久久精品国产亚洲精品| 特大巨黑吊av在线直播| 亚洲中文日韩欧美视频| 日本黄色视频三级网站网址| 少妇丰满av| 老汉色av国产亚洲站长工具| 欧美成人a在线观看| 欧美日韩福利视频一区二区| 久久人妻av系列| 国产精品乱码一区二三区的特点| 久久久久亚洲av毛片大全| 国产成人系列免费观看| 国产一区二区激情短视频| 亚洲熟妇熟女久久| 高清毛片免费观看视频网站| 亚洲男人的天堂狠狠| 国产激情偷乱视频一区二区| 亚洲成人免费电影在线观看| 国产一区在线观看成人免费| 久久久精品欧美日韩精品| 国产亚洲av嫩草精品影院| 日本一二三区视频观看| tocl精华| 一卡2卡三卡四卡精品乱码亚洲| 精品国内亚洲2022精品成人| 国产老妇女一区| 亚洲乱码一区二区免费版| 美女黄网站色视频| 亚洲激情在线av| 乱人视频在线观看| 在线观看日韩欧美| xxx96com| 变态另类成人亚洲欧美熟女| 午夜福利欧美成人| 天美传媒精品一区二区| 99热6这里只有精品| 尤物成人国产欧美一区二区三区| 亚洲国产精品合色在线| 午夜福利18| 亚洲熟妇熟女久久| 亚洲aⅴ乱码一区二区在线播放| 国产亚洲av嫩草精品影院| 一本久久中文字幕| 熟女少妇亚洲综合色aaa.| av专区在线播放| 国产精品1区2区在线观看.| 日本五十路高清| 丰满人妻一区二区三区视频av | 桃色一区二区三区在线观看| 国产v大片淫在线免费观看| 国产伦精品一区二区三区四那| 老司机在亚洲福利影院| 久久亚洲精品不卡| 69人妻影院| 熟女人妻精品中文字幕| 欧美大码av| 国产精品99久久久久久久久| 国产真实伦视频高清在线观看 | svipshipincom国产片| 精品久久久久久久人妻蜜臀av| 亚洲专区国产一区二区| 一a级毛片在线观看| or卡值多少钱| 国产日本99.免费观看| 久久草成人影院| 国产久久久一区二区三区| 人人妻人人看人人澡| 国产精品99久久99久久久不卡| 亚洲中文字幕一区二区三区有码在线看| 最新中文字幕久久久久| 日韩欧美免费精品| 久99久视频精品免费| 亚洲男人的天堂狠狠| 国产激情欧美一区二区| 亚洲中文字幕日韩| 色综合亚洲欧美另类图片| 女人十人毛片免费观看3o分钟| 亚洲人成伊人成综合网2020| 怎么达到女性高潮| 亚洲成人中文字幕在线播放| 久久久久久久精品吃奶| 一卡2卡三卡四卡精品乱码亚洲| 精品久久久久久久久久久久久| 欧美大码av| 国产午夜精品论理片| 国产探花极品一区二区| 婷婷精品国产亚洲av在线| 国产欧美日韩精品亚洲av| 欧美成人a在线观看| 一级黄片播放器| 丁香六月欧美| 国产一级毛片七仙女欲春2| 一a级毛片在线观看| 欧美日韩乱码在线| 免费看日本二区| 一区二区三区国产精品乱码| 日韩精品青青久久久久久| 亚洲自拍偷在线| 午夜福利免费观看在线| 狂野欧美白嫩少妇大欣赏| 99热精品在线国产| 日本免费一区二区三区高清不卡| 少妇的逼好多水| a在线观看视频网站| 亚洲av电影不卡..在线观看| 老司机在亚洲福利影院| 18禁国产床啪视频网站| 免费搜索国产男女视频| 高清在线国产一区| h日本视频在线播放| 欧美日韩黄片免| 色综合亚洲欧美另类图片| 欧美激情在线99| 久久久色成人| 国产成人a区在线观看| 人人妻,人人澡人人爽秒播| 国产高清视频在线观看网站| 成人av在线播放网站| 国产成人a区在线观看| 亚洲欧美日韩高清在线视频| 国产老妇女一区| 精品一区二区三区视频在线观看免费| 18禁国产床啪视频网站| 免费观看的影片在线观看| 国产精品 欧美亚洲| 欧美+日韩+精品| av中文乱码字幕在线| 亚洲无线在线观看| 中文字幕精品亚洲无线码一区| 高潮久久久久久久久久久不卡| 一个人看的www免费观看视频| 午夜福利欧美成人| 色精品久久人妻99蜜桃| 国产色爽女视频免费观看| 国产黄片美女视频| 欧美在线黄色| 一个人观看的视频www高清免费观看| xxxwww97欧美| 亚洲午夜理论影院| 午夜免费成人在线视频| 亚洲性夜色夜夜综合| 免费在线观看亚洲国产| 岛国在线观看网站| 一区福利在线观看| 搡老妇女老女人老熟妇| 狂野欧美激情性xxxx| 亚洲国产日韩欧美精品在线观看 | 亚洲精品影视一区二区三区av| 亚洲五月婷婷丁香| 成人欧美大片| 久久久久久九九精品二区国产| 琪琪午夜伦伦电影理论片6080| 国产精品美女特级片免费视频播放器| 香蕉久久夜色| 国产精品影院久久| 亚洲七黄色美女视频| 亚洲第一欧美日韩一区二区三区| 婷婷精品国产亚洲av| av欧美777| 欧美黑人巨大hd| 精品人妻一区二区三区麻豆 | 99在线人妻在线中文字幕| 精品熟女少妇八av免费久了| 久久久久久久久中文| 在线观看av片永久免费下载| 精品免费久久久久久久清纯| 国产成人啪精品午夜网站| 成人国产综合亚洲| 日本 av在线| 91久久精品国产一区二区成人 | 国产 一区 欧美 日韩| 国产熟女xx| 亚洲精品在线观看二区| 亚洲欧美一区二区三区黑人| 成人18禁在线播放| 悠悠久久av| 国产真实伦视频高清在线观看 | 无限看片的www在线观看| 国产av在哪里看| 亚洲av电影在线进入| 精品无人区乱码1区二区| 国产成+人综合+亚洲专区| 国产男靠女视频免费网站| 免费看美女性在线毛片视频| 精品午夜福利视频在线观看一区| 亚洲精品一卡2卡三卡4卡5卡| 国产精品av视频在线免费观看| 亚洲男人的天堂狠狠| 亚洲av电影不卡..在线观看| 成人特级黄色片久久久久久久| 日韩欧美三级三区| 在线看三级毛片| 亚洲成av人片在线播放无| 成年版毛片免费区| 久久99热这里只有精品18| 免费大片18禁| 色尼玛亚洲综合影院| 国产熟女xx| 国产精品综合久久久久久久免费| 偷拍熟女少妇极品色| 欧美成人一区二区免费高清观看| 国产单亲对白刺激| 国语自产精品视频在线第100页| 欧美午夜高清在线| 亚洲av免费高清在线观看| 免费人成在线观看视频色| 老司机在亚洲福利影院| 国产精品三级大全| 两人在一起打扑克的视频| 少妇熟女aⅴ在线视频| 久久久久久人人人人人| 国产一区二区亚洲精品在线观看| xxxwww97欧美| av天堂中文字幕网| 久久久久久久久久黄片| 熟女人妻精品中文字幕| 成人永久免费在线观看视频| 一a级毛片在线观看| 日韩高清综合在线| 亚洲中文日韩欧美视频| 午夜福利欧美成人| 男人和女人高潮做爰伦理| 成人一区二区视频在线观看| 99热只有精品国产| 久久久久精品国产欧美久久久| 村上凉子中文字幕在线| 日本撒尿小便嘘嘘汇集6| 999久久久精品免费观看国产| 日韩av在线大香蕉| 欧美日韩国产亚洲二区| 一个人免费在线观看的高清视频| 免费人成在线观看视频色| 国产精品一及| 成人亚洲精品av一区二区| 人人妻人人澡欧美一区二区| 三级毛片av免费| 人人妻,人人澡人人爽秒播| 日韩欧美免费精品| 观看免费一级毛片| 日韩免费av在线播放| 香蕉丝袜av| 内地一区二区视频在线| 成年女人看的毛片在线观看| 欧美性猛交黑人性爽| 中出人妻视频一区二区| 国产精品久久久久久久电影 | 一区二区三区免费毛片| 亚洲精品456在线播放app | 女人高潮潮喷娇喘18禁视频| 女人十人毛片免费观看3o分钟| 亚洲在线观看片| av天堂在线播放| 亚洲欧美日韩卡通动漫| 99精品在免费线老司机午夜| 人人妻人人澡欧美一区二区| 国产精品 国内视频| 久久国产精品人妻蜜桃| 三级男女做爰猛烈吃奶摸视频| 国产一区二区在线av高清观看| 美女高潮的动态| 黄色日韩在线| 亚洲熟妇中文字幕五十中出| 夜夜夜夜夜久久久久| 特级一级黄色大片| АⅤ资源中文在线天堂| 91麻豆av在线| 中文字幕av成人在线电影| 一夜夜www| 丰满人妻熟妇乱又伦精品不卡| 午夜福利高清视频| 制服丝袜大香蕉在线| 男女那种视频在线观看| 国产精品影院久久| 亚洲精华国产精华精| 久久6这里有精品| 日本撒尿小便嘘嘘汇集6| 免费无遮挡裸体视频| 国产欧美日韩精品亚洲av| 国产精品亚洲av一区麻豆| 一区二区三区激情视频| 综合色av麻豆| 色尼玛亚洲综合影院| 伊人久久大香线蕉亚洲五| 男女视频在线观看网站免费| 国产色婷婷99| 色在线成人网| 毛片女人毛片| 国产综合懂色| 久久精品91蜜桃| 一个人免费在线观看的高清视频| 夜夜躁狠狠躁天天躁| 无遮挡黄片免费观看| 性色avwww在线观看| 国产熟女xx| 精品乱码久久久久久99久播| 少妇的丰满在线观看| 色av中文字幕| 国产亚洲精品一区二区www| 欧美黑人欧美精品刺激| 精品日产1卡2卡| 一区二区三区激情视频| 午夜免费观看网址| 国产成人啪精品午夜网站| 好男人电影高清在线观看| 亚洲黑人精品在线| av专区在线播放| 久久久成人免费电影| 熟女人妻精品中文字幕| 三级男女做爰猛烈吃奶摸视频| 国产一区二区三区在线臀色熟女| 看黄色毛片网站| 国产蜜桃级精品一区二区三区| x7x7x7水蜜桃| 变态另类成人亚洲欧美熟女| 听说在线观看完整版免费高清| 狂野欧美白嫩少妇大欣赏| 丰满人妻熟妇乱又伦精品不卡| 99久久精品一区二区三区| 久久精品国产99精品国产亚洲性色| 成人亚洲精品av一区二区| 99久久99久久久精品蜜桃| aaaaa片日本免费| www.熟女人妻精品国产| 一区二区三区激情视频| 悠悠久久av| 一二三四社区在线视频社区8| 俺也久久电影网| 午夜久久久久精精品| 999久久久精品免费观看国产| 露出奶头的视频| 中出人妻视频一区二区| 国产aⅴ精品一区二区三区波| 手机成人av网站| 久久精品夜夜夜夜夜久久蜜豆| 一个人观看的视频www高清免费观看| 在线播放国产精品三级| 中文字幕久久专区| 国产麻豆成人av免费视频| 亚洲欧美日韩高清在线视频| 三级国产精品欧美在线观看| 国产黄a三级三级三级人| 别揉我奶头~嗯~啊~动态视频| 日韩大尺度精品在线看网址| 国产一区二区三区视频了| 亚洲无线在线观看| 国产单亲对白刺激| 国产午夜精品论理片| 亚洲专区国产一区二区| 日韩欧美一区二区三区在线观看| 久久天躁狠狠躁夜夜2o2o| 精品99又大又爽又粗少妇毛片 | 精品一区二区三区视频在线 | 国产成人啪精品午夜网站| 亚洲无线观看免费| 男人舔奶头视频| 99久久久亚洲精品蜜臀av| 欧美一区二区精品小视频在线| 69av精品久久久久久| 真人做人爱边吃奶动态| 国产91精品成人一区二区三区| 亚洲成a人片在线一区二区| 国产午夜福利久久久久久| 丝袜美腿在线中文| 在线观看舔阴道视频| 国产一区二区在线av高清观看| 久久国产精品人妻蜜桃| 18禁黄网站禁片免费观看直播| 欧美乱妇无乱码| 欧美乱色亚洲激情| 老汉色av国产亚洲站长工具| 国产久久久一区二区三区| 欧美av亚洲av综合av国产av| 可以在线观看的亚洲视频| 欧美区成人在线视频| 两个人视频免费观看高清| 91九色精品人成在线观看| 久久精品影院6| 国产精品一区二区免费欧美| 精品电影一区二区在线| 十八禁人妻一区二区| 国产亚洲精品一区二区www| 国产av不卡久久| 国产精品久久电影中文字幕| 亚洲真实伦在线观看| 国产精品嫩草影院av在线观看 | 午夜a级毛片| 在线免费观看不下载黄p国产 | 欧美一级毛片孕妇| 最好的美女福利视频网| 亚洲人成电影免费在线| 精品一区二区三区视频在线观看免费| 日本 欧美在线| 色老头精品视频在线观看| 深爱激情五月婷婷| 亚洲中文日韩欧美视频| 哪里可以看免费的av片| 久久人妻av系列| 他把我摸到了高潮在线观看| 日本精品一区二区三区蜜桃| 两人在一起打扑克的视频| 亚洲专区中文字幕在线| 1024手机看黄色片| 99视频精品全部免费 在线| 日本a在线网址| 国产午夜精品论理片| 国产精品爽爽va在线观看网站| 精品国内亚洲2022精品成人| 国产成年人精品一区二区| 尤物成人国产欧美一区二区三区| 男女下面进入的视频免费午夜| 婷婷丁香在线五月| 免费看a级黄色片| 久99久视频精品免费| 人妻夜夜爽99麻豆av| 国产成人影院久久av| 国语自产精品视频在线第100页| 亚洲七黄色美女视频| 欧美日韩国产亚洲二区| 极品教师在线免费播放| 午夜日韩欧美国产| 免费高清视频大片| 最新美女视频免费是黄的| 国产伦精品一区二区三区视频9 | 小说图片视频综合网站| 69av精品久久久久久| 少妇高潮的动态图| 免费大片18禁| 成人精品一区二区免费| 丰满人妻一区二区三区视频av | 在线观看66精品国产| 女人十人毛片免费观看3o分钟| 搞女人的毛片| 免费观看人在逋| 波多野结衣巨乳人妻| 丰满乱子伦码专区| 高潮久久久久久久久久久不卡| 日韩欧美三级三区| 国产av一区在线观看免费| 婷婷六月久久综合丁香| svipshipincom国产片| av欧美777| 亚洲精品乱码久久久v下载方式 | 国产单亲对白刺激| 久久国产乱子伦精品免费另类| 欧美日韩福利视频一区二区| 身体一侧抽搐| 久久香蕉国产精品| 丁香六月欧美| 琪琪午夜伦伦电影理论片6080| 免费一级毛片在线播放高清视频| 99久久精品热视频| 亚洲第一电影网av| 国产精品,欧美在线| 国产精品自产拍在线观看55亚洲| 一区二区三区国产精品乱码| 欧美日韩中文字幕国产精品一区二区三区| 免费一级毛片在线播放高清视频| 亚洲精品粉嫩美女一区| 国产成人a区在线观看| 老汉色av国产亚洲站长工具| 啦啦啦韩国在线观看视频| 99精品欧美一区二区三区四区| 两个人看的免费小视频| 久久久成人免费电影| 精品人妻一区二区三区麻豆 | 99久久精品热视频| 日韩欧美在线二视频| www.色视频.com| 精品久久久久久久久久免费视频| 在线观看一区二区三区| 特大巨黑吊av在线直播| 18禁黄网站禁片免费观看直播| 成熟少妇高潮喷水视频| 搡老岳熟女国产| 欧美中文日本在线观看视频| 啦啦啦免费观看视频1| 一卡2卡三卡四卡精品乱码亚洲| 日韩 欧美 亚洲 中文字幕| 老司机午夜福利在线观看视频| 精华霜和精华液先用哪个| 亚洲精品456在线播放app | 观看免费一级毛片| 深夜精品福利| 啦啦啦观看免费观看视频高清| 深夜精品福利| 日日夜夜操网爽| 19禁男女啪啪无遮挡网站| 国产激情欧美一区二区| 亚洲av日韩精品久久久久久密| 真人做人爱边吃奶动态| 最近视频中文字幕2019在线8| 一区二区三区国产精品乱码| 久久久久久大精品| 亚洲中文日韩欧美视频| a级毛片a级免费在线| 中文资源天堂在线| 欧美一区二区亚洲| 亚洲成av人片在线播放无| 超碰av人人做人人爽久久 | 波多野结衣高清无吗| 国产成人欧美在线观看| 午夜影院日韩av| 欧美一级毛片孕妇| 国产伦精品一区二区三区视频9 | 老熟妇仑乱视频hdxx| 午夜激情欧美在线| 伊人久久大香线蕉亚洲五| 免费在线观看成人毛片| 午夜福利18| 亚洲av不卡在线观看| 搡女人真爽免费视频火全软件 | 国产在线精品亚洲第一网站| 亚洲熟妇熟女久久| 国产亚洲精品av在线| 国产精品女同一区二区软件 | 亚洲人与动物交配视频| 99国产极品粉嫩在线观看| 久久草成人影院| 久久久成人免费电影| 国产探花在线观看一区二区| 嫩草影院精品99| 色老头精品视频在线观看| 99riav亚洲国产免费| 国产高潮美女av| 欧美丝袜亚洲另类 | 日日摸夜夜添夜夜添小说| 天天添夜夜摸| 在线十欧美十亚洲十日本专区| 一边摸一边抽搐一进一小说| 91久久精品电影网| 一级黄片播放器| 国产一区在线观看成人免费| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 老熟妇乱子伦视频在线观看| 男插女下体视频免费在线播放| 日韩人妻高清精品专区| 很黄的视频免费| 国产不卡一卡二| 制服人妻中文乱码| 热99在线观看视频| 人妻丰满熟妇av一区二区三区| 天天添夜夜摸| 国产黄片美女视频| 午夜精品久久久久久毛片777| 亚洲乱码一区二区免费版| 国产久久久一区二区三区| 中文字幕人妻丝袜一区二区| 小蜜桃在线观看免费完整版高清| 狠狠狠狠99中文字幕| 成人特级黄色片久久久久久久| 无人区码免费观看不卡| 午夜福利欧美成人| 久久人妻av系列| 综合色av麻豆| 老司机在亚洲福利影院| 成人性生交大片免费视频hd|