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

    基于DOE的二沖程汽油機Wiebe燃燒參數(shù)優(yōu)化研究

    2018-03-01 03:24:24王在良楊海青王思奇
    關鍵詞:缸內(nèi)扭矩數(shù)值

    王在良,楊海青,王思奇

    (南京航空航天大學 能源與動力學院, 南京 210016)

    汽油機排量小、轉(zhuǎn)速高,進、排氣管中的非定常流對工作過程有強烈的影響,致使數(shù)值計算研究相對比較困難,但汽油機為大批量生產(chǎn)、成本低,容易獲得試驗樣機,因此,長期以來,主要是通過臺架試驗的方式對其進行研究[1]。近年來,隨著計算機輔助工程(CAE)的巨大進步,其預測能力已經(jīng)發(fā)展到相當成熟的階段,許多設計驗證可通過計算機完成。對于燃燒模型的探究,現(xiàn)有的數(shù)學模型包括零維、準維和多維燃燒模型[2-3]。零維燃燒模型由于具有參數(shù)間關系簡單、容易求解等諸多優(yōu)點而被廣泛應用,通常用于計算速度特性、負荷特性等汽油機的各項性能參數(shù)。零維燃燒模型計算中較常采用韋伯半經(jīng)驗公式,其主要是通過參考所研究發(fā)動機的運行數(shù)據(jù),選定Wiebe公式中的一些經(jīng)驗參數(shù)后,用Wiebe公式模擬實際發(fā)動機的燃燒放熱率,最終實現(xiàn)發(fā)動機的數(shù)值模擬研究。然而,這些經(jīng)驗參數(shù)根據(jù)機型的不同變化范圍較大,經(jīng)驗參數(shù)的選取對模型的正確性有至關重要的影響。以往通過經(jīng)驗或多次驗證調(diào)整的方式選取的方法具有很大的盲目性,且工作量大,計算結果與實測值出入較大。在GT-Power中,基于DOE(design of experiment)試驗設計研究方法在發(fā)動機性能優(yōu)化方面是一種有效的研究方法。胡磊等[4]將發(fā)動機噴油系統(tǒng)參數(shù)作為DOE試驗設計變量,通過優(yōu)化得到最佳噴油參數(shù),使得功率、排放達到設計要求指標。Rahul C Chikurde等[5]利用DOE對排氣噪聲做了優(yōu)化研究,詳細討論了排氣系統(tǒng)各個結構參數(shù)對發(fā)動機的影響。因此,DOE方法對模型的校核與優(yōu)化有很大幫助。

    本文針對1臺直列兩缸風冷進氣管電噴式二沖程汽油機,借助于發(fā)動機性能仿真軟件 GT-Power,應用DOE(design of experiment)研究方法,對Wiebe燃燒模型特性參數(shù)進行了擬合優(yōu)化,得到了發(fā)動機在各個工況下的燃燒模型參數(shù)變化MAP圖,最終建立了準確、完整的發(fā)動機一維仿真模型。

    1 建模與試驗驗證

    本文研究對象是1臺Hirth公司的直列兩缸風冷二沖程汽油機,其進氣方式為簧片閥進氣,采用進氣管電控燃油噴射,主要技術參數(shù)如表1所示。

    表1 主要技術參數(shù)

    1.1 一維仿真模型的建立

    根據(jù)發(fā)動機結構,將發(fā)動機按照進氣系統(tǒng)、缸體、排氣系統(tǒng)三大部分進行搭建。為使數(shù)值模擬盡可能地反映發(fā)動機運行情況,需根據(jù)發(fā)動機的結構,準確地輸入其結構參數(shù),同時選擇合理的工作過程計算模型。本研究中傳熱模型選擇Woschni傳熱模型。燃燒模型選擇非預測型的SI Wiebe燃燒模型,該模型通過參考所研究發(fā)動機的運行數(shù)據(jù),選定韋伯公式中的一些經(jīng)驗參數(shù)后,用韋伯公式模擬實際發(fā)動機的燃燒放熱率,SI Wiebe燃燒模型計算方法如式(1)所示[6]。

    Combustion(θ)=(CE)[1-e-(WC)(θ-SOC)(E+1)]

    (1)

    式中θ為瞬時曲軸角度。

    燃燒模型中的計算常數(shù):

    (2)

    WC為Wiebe常數(shù);

    (3)

    SOC為燃燒起始常數(shù);

    BMC=-ln(1-BM)

    (4)

    BMC為燃燒中間常數(shù);

    BSC=-ln(1-BS)

    (5)

    BSC為燃燒起始常數(shù);

    BEC=-ln(1-BE)

    (6)

    BEC為燃燒結束常數(shù);

    Wiebe方程的輸入量中:AA為錨角(本研究中為Wiebe曲線中燃料燃50%時刻);D為燃燒持續(xù)轉(zhuǎn)角;E為燃燒品質(zhì)指數(shù);CE為燃燒效率;BM為AA角度下的已燃燃油比例,取50%;BS為燃燒持續(xù)開始時的已燃燃油比例,取10%;BE為燃燒持續(xù)結束時的已燃燃油比例,取90%。

    在GT-Power中,Wiebe函數(shù)中決定燃料燃燒特性的關鍵參數(shù)為AA、D和E,而AA與D反映的是燃燒的快慢以及持續(xù)時間,E決定放熱率曲線的形狀,因此燃燒特性參數(shù)的準確性是建立仿真模型的關鍵。所建立的仿真模型結構如圖1所示。

    圖1 仿真模型結構

    1.2 模型的驗證

    圖2為發(fā)動機在全負荷工況下,仿真模擬運行得到的進氣流量、扭矩及功率與臺架試驗值的結果對比。由圖2可知:發(fā)動機進氣流量在各轉(zhuǎn)速下達到了很好的吻合,說明建模中發(fā)動機結構幾何參數(shù)設置滿足實際發(fā)動機一致性要求,包括進氣、排氣系統(tǒng)管路、簧片閥及其缸內(nèi)幾何結構等的設置。而扭矩和功率在有些轉(zhuǎn)速工況下偏差較大,最大達到了10%。其主要原因是建模過程中在選擇燃燒參數(shù)時對Wiebe參數(shù)是以經(jīng)驗選取的,難以準確地反映在各個工況下的發(fā)動機缸內(nèi)混合氣燃燒特性,故而導致較大的建模誤差。因此,本文在建立一維性能仿真模型中,采用DOE方法對SI Wiebe燃燒模型燃燒控制參數(shù)做系統(tǒng)的試驗選取研究,以達到數(shù)值模擬計算要求。

    圖2 試驗值與模擬值對比

    2 DOE方法介紹與設置

    試驗設計(design of experiment)簡稱 DOE,是一種結構性的系統(tǒng)研究自變量與因變量之間關系的研究方法。該方法是通過GT-Power 中DOE模塊設置和運行來研究一個或多個輸入?yún)?shù)(因子)對多個輸出結果(響應)的影響。它不僅可以呈現(xiàn)因子間的相互關系,還能夠通過合理減少試驗數(shù)量,縮短產(chǎn)品開發(fā)周期,節(jié)省研究成本[7]。在 DOE 中的多個試驗運行計算完成后,通過 DOE 后處理器可以對試驗結果進行響應面擬合分析和優(yōu)化設計。在DOE設置中,可以選擇多種不同形式的DOE抽樣方法,包括全因子法、D-Optimun以及拉丁超立方法。本文使用拉丁超立方法。拉丁超立方法是一種受約束的多維分層抽樣方法,對給定的試驗次數(shù),可將水平區(qū)間等分成互不重疊的子區(qū)間,隨后在各個子區(qū)間上進行等概率抽樣?;纠碚撊缦耓8-9]:假設K維隨機變量X的各個元素的概率分布函數(shù)為Fi(i=1,2,…,K)。向量X的各元素相互獨立,每個元素進行N次抽樣,xjk為第k(k=1,2,…,K)個元素的第j(j=1,2,…,N)次抽樣的值,定義N×K維矩陣P。P的每一列由數(shù)列{1,2,…,N}中各元素的隨機排列組成。 令隨機變量ξjk服從區(qū)間[0,1]上的均勻分布,則抽樣后得到的結果為:

    (7)

    式中pjk為N×K維矩陣P的j行k列元素。

    表2 Wiebe公式中參數(shù)的取值

    拉丁超立方抽樣能滿足試驗設計要求的均勻性和正交性的指標,通過少量的試驗就能對響應面進行精確的擬合。本研究分別將SI Wiebe燃燒模型中決定燃燒特性的關鍵參數(shù)AA(燃料燃50%時刻)、D(燃燒10%~90%持續(xù)轉(zhuǎn)角)及E(燃燒品質(zhì)指數(shù))作為DOE輸入因子進行運算,優(yōu)化范圍如表2所示,每個工況點設置100個試驗設計計算樣本。

    3 基于DOE的擬合分析

    3.1 響應面的擬合

    建立響應面的主要目的是令模型能針對范圍內(nèi)的不同因子輸入?yún)?shù)進行插值得出所需響應,通過在響應面搜索合適的插值,得出所需最優(yōu)化的最佳響應[10]。對每個轉(zhuǎn)速工況中的因子和所對應的響應,采用簡單最小二次多項式擬合法(OLS)進行響應面的擬合,擬合時選取擬合階數(shù)為2階。圖3為擬合結束后在轉(zhuǎn)速為4 000 r/min時AA與D對扭矩的面響應。

    圖3 扭矩的面響應

    3.2 質(zhì)量評價

    DOE優(yōu)化是依靠對面響應的插值得到優(yōu)化結果,因此面響應擬合質(zhì)量對優(yōu)化選擇結果影響很大。響應面的擬合質(zhì)量評價主要由3個評價指標組成,分別為 R-Sqr 指標、Adj.R-sqr 指標和 Q-Sqr 指標。其中,R-Sqr 指標計算公式如式(8)所示。

    (8)

    式中:YP為預測響應值;YO為觀測響應值;n為試驗次數(shù)。扭矩、功率和燃油消耗率的R2均在0.99以上,滿足優(yōu)化計算要求。

    同時,選取油門開度為60%、轉(zhuǎn)速為4 000 r/min工況下扭矩和功率的預測數(shù)據(jù)和模擬數(shù)據(jù)進行分析,以檢驗抽樣數(shù)據(jù)和所建模型的合理性。圖4為扭矩和功率響應擬合殘差,x軸為模擬所得的觀測值,y軸為抽樣預測值,兩者的數(shù)值越向中間對角線靠攏,表示誤差越小,該響應擬合精度越高。

    圖4 響應擬合殘差

    3.3 顯著性分析

    圖5為不同因子之間對扭矩的影響關系,分為主效應與相互效應。主效應主要是對各個因子進行單獨的分析,根據(jù)其對所選響應的影響程度,以量化的數(shù)值來表征對該響應的相關影響。如圖5所示,E值為正,說明E值與響應正相關,E值增大,響應數(shù)值也增大;AA和D值的幅值為負,說明與響應負相關,因子數(shù)值增大,響應數(shù)值會減少。相互效應主要表征各個因子之間相乘組成的項對不同響應的影響。由于擬合的響應面為一組多項式,因此,相互響應能探究多項式中不同因子的組合對響應變量的相關性影響。圖中的AA*D、AA*E均為正值,說明其乘積越大,對響應的影響就越大;D*E為負值,且數(shù)值很小,由此可知該相乘項對響應的影響很小。

    3.4 燃燒特性分析

    利用DOE計算后的結果數(shù)據(jù),對Wiebe燃燒模型中關鍵參數(shù)以扭矩變化為例做詳細的分析。

    圖5 扭矩主效應

    由圖6可知:當在全負荷、轉(zhuǎn)速為3 500 r/min時,隨著AA值增大,發(fā)動機扭矩呈現(xiàn)增大后減小趨勢。這是因為在轉(zhuǎn)速相對較小時,缸內(nèi)燃燒時間相對充足,再加上缸內(nèi)新鮮充量比較多,若AA值過小,說明過早開始燃燒,在還未達到上止點前就已經(jīng)有大量的已燃燃油,產(chǎn)生燃燒負功,導致扭矩輸出低;當AA值逐步增大到7 ℃A BTDC時,形成負功現(xiàn)象不再產(chǎn)生,扭矩最大達到40.27 N·m;但當AA值從7 ℃A BTDC進一步增大到15 ℃A BTDC時,意味著燃燒過程后移,更多的燃油在活塞下行時燃燒,此時缸內(nèi)容積變大,導致缸內(nèi)燃燒爆發(fā)壓力下降,致使扭矩輸出下降了2.5%。在6 000 r/min時,在相同的燃燒時間內(nèi)所占的曲軸轉(zhuǎn)角變小,因此,AA值與扭矩從3 ℃A BTDC開始已經(jīng)呈負相關趨勢,在樣本AA值為3~15 ℃A BTDC范圍內(nèi)扭矩下降了6%。

    圖6 AA值與扭矩變化關系

    從圖7可以看出:隨著D值的增大,燃燒持續(xù)時間變長,燃燒平緩,缸內(nèi)爆發(fā)壓力升高率變小,峰值壓力下降,發(fā)動機扭矩輸出則會相應減小。

    燃燒品質(zhì)指數(shù)E是表征放熱率分布的一個參數(shù),其大小影響放熱率曲線的形狀,E值較小,初期放熱量多,壓力升高率大,燃燒粗暴:反之E值增大,初期放熱量小,放熱率曲線峰值向后移,壓力升高率小,燃燒柔和。由圖8可知,E值的變化對發(fā)動機動力輸出特性有重要的影響。

    圖7 D與扭矩變化關系

    4 GA優(yōu)化與結果驗證

    在采用擬合精度較好的響應面的基礎下,使用基于GA Standard遺傳算法進行優(yōu)化。為了達到與臺架試驗數(shù)據(jù)相對一致的仿真計算模型,優(yōu)化目標變量設置為扭矩、功率和燃油消耗率的目標值(來自于發(fā)動機臺架試驗)。設置最大迭代次數(shù)為100次,在20個迭代后沒有變化時判斷為收斂,取樣個體為40個,突變率為10%,采用單點交叉方式進行雜交,最優(yōu)淘汰制進行選擇,統(tǒng)一變異方式進行突變。優(yōu)化結果如圖9所示。最終得到完整的在各個工況下的Wiebe參數(shù)MAP。

    在半經(jīng)驗的Wiebe公式中,掌握決定燃燒放熱率的特性參數(shù)變化規(guī)律是建立仿真模型的首要條件。如圖9(a)(b)所示,AA值和D值總體呈現(xiàn)隨著負荷的增加而減小,原因是在小負荷工況下,發(fā)動機缸內(nèi)進氣量少,廢氣稀釋現(xiàn)象嚴重,導致混合氣濃度偏低,火焰?zhèn)鞑ニ俣嚷?,因此燃燒持續(xù)期會變長,AA值也比較大;而在大負荷工況時,缸內(nèi)進氣量充足,燃燒速率變快,則會使燃燒時間縮短。

    同時,AA值反映著點火時刻的大小,圖9(a)中AA值在隨轉(zhuǎn)速變化的過程中,總體呈現(xiàn)先增大后減小的趨勢,因為當在小轉(zhuǎn)速時,轉(zhuǎn)過相同的轉(zhuǎn)角時所用的時間較長,缸內(nèi)混合氣可在較小的曲軸轉(zhuǎn)角內(nèi)完成燃燒過程,故AA值相對較??;當在中間轉(zhuǎn)速時,AA值變大,因為此時燃燒條件逐漸變好,燃燒效率變高,為避免過早燃燒,應適當?shù)赝七t點火,AA值較大;而在高轉(zhuǎn)速時,轉(zhuǎn)過相同曲軸轉(zhuǎn)角的時間變短,無法滿足燃燒所需時間要求,會產(chǎn)生后燃現(xiàn)象,為使達到更高的指示功率,應適當?shù)靥崆包c火,則AA值變小。

    由圖9(b)可知:在高速大負荷工況時,為發(fā)動機標定工況,燃燒條件適宜,燃燒持續(xù)期D值也將最短;在低速大負荷時,雖然進氣量充足,但缸內(nèi)湍流強度弱,火焰?zhèn)鞑ニ俣嚷?,故相較于高速工況D值會稍有變長;當發(fā)動機處于低速小負荷時,缸內(nèi)進氣量少,火焰?zhèn)鞑ニ俣冗M一步降低,致使D值明顯變大;在低速大負荷時,進氣損失進一步增加,混合氣變稀,且相同燃燒時間內(nèi)所轉(zhuǎn)過的曲軸轉(zhuǎn)角也會變大,故D值也會在整個發(fā)動機運行工況內(nèi)最大。

    圖9(c)中E值在發(fā)動機運行整個工況內(nèi)的變動范圍為1.5~3,符合汽油機的燃燒特點。

    通過優(yōu)化后,將發(fā)動機的輸出扭矩、功率與試驗數(shù)值進行對比,誤差值都控制在4%左右,較初始建模誤差有很大改善,可為發(fā)動機模擬各項性能研究做良好的保障。圖10為優(yōu)化后在節(jié)氣門開度為60%和100%的扭矩、功率對比結果,表3、4分別為扭矩和功率在各個工況下誤差對比。

    圖9 Wiebe參數(shù)MAP

    圖10 優(yōu)化后試驗與模擬對比

    表3 扭矩誤差對比

    表4 功率誤差對比

    5 結束語

    基于DOE的Wiebe燃燒特性參數(shù)優(yōu)化,運用GA Standard遺傳算法能夠?qū)崿F(xiàn)燃燒參數(shù)優(yōu)化的目的。本文在保證所建模型試驗環(huán)境條件、進氣流量、摩擦功率等與發(fā)動機試驗工況具有很高一致性的前提下,通過對AA、D、E值采用拉丁超立方抽樣選取,選用遺傳算法優(yōu)化,使一維整機性能模型的性能輸出與臺架試驗達到很好的吻合。

    發(fā)動機缸內(nèi)燃燒狀態(tài)非常惡劣,很難利用儀器直接測試研究。本研究從燃燒特性角度對建模做出了詳細的優(yōu)化分析研究,揭示了燃燒參數(shù)在各個工況下的差異性規(guī)律,為數(shù)值建模提供一定的理論指導。

    [1] 朱訪君,吳堅.內(nèi)燃機工作過程數(shù)值計算及其優(yōu)化[M].北京:國防工業(yè)出版社,1997:215-222.

    [2] 周松,王銀燕,明平劍.內(nèi)燃機工作過程仿真技術[M].北京:北京航空航天大學出版社,2012.

    [3] 陳明,丁寧,張小矛.車用內(nèi)燃機仿真[M].北京:人民交通出版社股份有限公司,2016:5.

    [4] 胡磊,石君明,孫躍東,等.基于DOE的柴油機噴油系統(tǒng)參數(shù)優(yōu)化研究 [J].合肥工業(yè)大學學報,2015,38(6):744-747,832.

    [5] RAHUL C CHIKURDE,MOHNISH KUMAR,TEJINDERPAL SINGH.Optimization and validation of exhaust tailpipe noise for passenger car[J].SAE International,2013(26):101.

    [6] Gamma Technologies.GT-SUITE Engine Performance Application Manual VERSION7.4[Z].2013:35-36.

    [7] 倪計民,杜倩穎,周英杰,等.DOE 在高壓共軌柴油機優(yōu)化設計中的應用[J].內(nèi)燃機學報,2009,27(3):231-236.

    [8] 鄧乾旺,文文.基于拉丁超立方抽樣的薄板裝配誤差分析[J].中國機械工程,2012,23(8):947-951.

    [9] MICHAEL D SHIELDS,ZHANG JIAXIN.The generalization of latin hypercube sampling[J].Reliability Engineering and Safety,2016,148:96-108.

    [10] 劉德輝.基于GT-Power的490柴油機進排氣系統(tǒng)的數(shù)值仿真與優(yōu)化研究[D].廣州:華南理工大學,2016.

    猜你喜歡
    缸內(nèi)扭矩數(shù)值
    電動汽車扭矩濾波控制策略研究
    用固定數(shù)值計算
    數(shù)值大小比較“招招鮮”
    水垂比對管柱摩阻和扭矩的影響實驗
    中國煤層氣(2021年5期)2021-03-02 05:53:14
    缸內(nèi)直噴發(fā)動機的燃燒和爆震仿真分析
    北京汽車(2018年6期)2019-01-03 11:27:38
    液壓扭矩扳手的不確定度評定
    基于Fluent的GTAW數(shù)值模擬
    焊接(2016年2期)2016-02-27 13:01:02
    氣密封接頭最佳上扣扭矩計算與分析
    支持直接噴射汽油機歐6開發(fā)目標的缸內(nèi)和循環(huán)可辨顆粒生成的評估
    淺析汽油機缸內(nèi)直噴技術
    欧美在线黄色| 国产精品免费一区二区三区在线| 一级黄色大片毛片| 亚洲国产高清在线一区二区三 | 人人妻,人人澡人人爽秒播| 欧美午夜高清在线| 淫秽高清视频在线观看| 午夜亚洲福利在线播放| 免费在线观看黄色视频的| 亚洲男人的天堂狠狠| 看黄色毛片网站| 国产精品99久久99久久久不卡| 妹子高潮喷水视频| 国产男靠女视频免费网站| 午夜a级毛片| 亚洲男人天堂网一区| 在线天堂中文资源库| 大型av网站在线播放| 国产成人一区二区三区免费视频网站| 欧美日韩亚洲国产一区二区在线观看| 亚洲avbb在线观看| 97超级碰碰碰精品色视频在线观看| 亚洲一区高清亚洲精品| 亚洲天堂国产精品一区在线| 日日夜夜操网爽| 日韩国内少妇激情av| 侵犯人妻中文字幕一二三四区| 欧美日韩黄片免| 99国产极品粉嫩在线观看| 日韩 欧美 亚洲 中文字幕| 香蕉av资源在线| 国产三级黄色录像| 日日干狠狠操夜夜爽| 久久久久九九精品影院| 久久这里只有精品19| 午夜福利成人在线免费观看| 国产久久久一区二区三区| 怎么达到女性高潮| 欧美绝顶高潮抽搐喷水| 99久久久亚洲精品蜜臀av| 亚洲第一青青草原| 视频在线观看一区二区三区| 欧美中文综合在线视频| 手机成人av网站| www日本黄色视频网| 日本a在线网址| 少妇的丰满在线观看| av电影中文网址| 午夜免费成人在线视频| 90打野战视频偷拍视频| 国产极品粉嫩免费观看在线| 国产精品野战在线观看| 久久国产乱子伦精品免费另类| 亚洲第一av免费看| 亚洲精品粉嫩美女一区| 午夜福利在线观看吧| 性欧美人与动物交配| 久久 成人 亚洲| 黄网站色视频无遮挡免费观看| 97碰自拍视频| 在线免费观看的www视频| 日本免费a在线| 日本 欧美在线| 18禁黄网站禁片午夜丰满| 亚洲avbb在线观看| 欧美日本视频| 日本免费一区二区三区高清不卡| 亚洲第一av免费看| 欧美激情 高清一区二区三区| 99在线人妻在线中文字幕| 国产av不卡久久| 99国产精品一区二区蜜桃av| 久久精品成人免费网站| x7x7x7水蜜桃| 久久久久精品国产欧美久久久| 热re99久久国产66热| 国产乱人伦免费视频| 成人国产一区最新在线观看| 男女做爰动态图高潮gif福利片| 欧美激情极品国产一区二区三区| 免费人成视频x8x8入口观看| 中文字幕人妻熟女乱码| 亚洲avbb在线观看| 国产激情偷乱视频一区二区| 侵犯人妻中文字幕一二三四区| 久久国产精品男人的天堂亚洲| 长腿黑丝高跟| av在线天堂中文字幕| 女人被狂操c到高潮| 欧美性长视频在线观看| 国产精品香港三级国产av潘金莲| 久久久久久久精品吃奶| 亚洲第一欧美日韩一区二区三区| 一级毛片精品| 欧美乱色亚洲激情| 亚洲国产欧美网| 亚洲av五月六月丁香网| 日韩免费av在线播放| 国产又色又爽无遮挡免费看| 亚洲精品粉嫩美女一区| 日本精品一区二区三区蜜桃| 亚洲激情在线av| 亚洲人成77777在线视频| 又大又爽又粗| 在线av久久热| 免费观看精品视频网站| 国产精品国产高清国产av| 亚洲成人久久性| 一区福利在线观看| 黄色视频,在线免费观看| 91国产中文字幕| 成人亚洲精品av一区二区| 亚洲中文日韩欧美视频| 欧美三级亚洲精品| 精品国产一区二区三区四区第35| 女性生殖器流出的白浆| 亚洲色图 男人天堂 中文字幕| 精品国内亚洲2022精品成人| 自线自在国产av| 91av网站免费观看| 国产97色在线日韩免费| www日本在线高清视频| 国内揄拍国产精品人妻在线 | 精品高清国产在线一区| 最好的美女福利视频网| 欧美日韩乱码在线| √禁漫天堂资源中文www| 在线观看免费视频日本深夜| www.999成人在线观看| 国产av又大| 久久 成人 亚洲| 村上凉子中文字幕在线| 男女视频在线观看网站免费 | 在线观看66精品国产| 久久久久久久久久黄片| 亚洲人成电影免费在线| 香蕉久久夜色| 国产精品香港三级国产av潘金莲| 十八禁人妻一区二区| 中文字幕av电影在线播放| 午夜免费成人在线视频| 99国产精品一区二区蜜桃av| 制服人妻中文乱码| 国产成人欧美| or卡值多少钱| 亚洲第一欧美日韩一区二区三区| 麻豆成人av在线观看| 一二三四社区在线视频社区8| bbb黄色大片| 成人永久免费在线观看视频| 国产av一区二区精品久久| 一本一本综合久久| 中文字幕高清在线视频| 欧美av亚洲av综合av国产av| 69av精品久久久久久| 成人18禁高潮啪啪吃奶动态图| 亚洲国产欧美日韩在线播放| 国产精品自产拍在线观看55亚洲| 欧美一级毛片孕妇| 在线观看免费视频日本深夜| 久久香蕉国产精品| 日韩一卡2卡3卡4卡2021年| 国产激情偷乱视频一区二区| 亚洲人成伊人成综合网2020| 99久久精品国产亚洲精品| 可以在线观看毛片的网站| 在线观看66精品国产| 久久99热这里只有精品18| 免费女性裸体啪啪无遮挡网站| 中文字幕精品亚洲无线码一区 | 每晚都被弄得嗷嗷叫到高潮| 久久久久久久久久黄片| 午夜福利在线观看吧| 三级毛片av免费| 久久久久国产一级毛片高清牌| 一边摸一边做爽爽视频免费| 国产精品香港三级国产av潘金莲| 少妇裸体淫交视频免费看高清 | 国产成人精品无人区| 国产亚洲欧美98| 国产激情偷乱视频一区二区| 亚洲国产欧美日韩在线播放| 亚洲av美国av| 久久久久九九精品影院| 国产精品一区二区免费欧美| 久久国产精品影院| 男女做爰动态图高潮gif福利片| 久久久国产成人免费| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲三区欧美一区| 亚洲无线在线观看| 国产亚洲欧美精品永久| 老司机靠b影院| 日本a在线网址| 人人澡人人妻人| 亚洲精品一区av在线观看| av超薄肉色丝袜交足视频| 99精品在免费线老司机午夜| 国产又爽黄色视频| 91国产中文字幕| 一本综合久久免费| 国产一卡二卡三卡精品| 老司机午夜十八禁免费视频| 精品福利观看| 国产亚洲欧美在线一区二区| 校园春色视频在线观看| 成人av一区二区三区在线看| 51午夜福利影视在线观看| 亚洲在线自拍视频| 在线视频色国产色| 大香蕉久久成人网| 好看av亚洲va欧美ⅴa在| 18禁国产床啪视频网站| 免费在线观看亚洲国产| 午夜两性在线视频| 俄罗斯特黄特色一大片| 亚洲天堂国产精品一区在线| 最好的美女福利视频网| 成人手机av| 两性午夜刺激爽爽歪歪视频在线观看 | 女生性感内裤真人,穿戴方法视频| 精品久久蜜臀av无| 欧美日韩黄片免| 国产高清有码在线观看视频 | 亚洲一区二区三区不卡视频| 久久中文看片网| av片东京热男人的天堂| 亚洲成人免费电影在线观看| 少妇裸体淫交视频免费看高清 | 女人爽到高潮嗷嗷叫在线视频| 91麻豆精品激情在线观看国产| 成年免费大片在线观看| 国产一级毛片七仙女欲春2 | 亚洲精品中文字幕在线视频| 国产高清有码在线观看视频 | 欧美精品啪啪一区二区三区| 亚洲熟妇熟女久久| 校园春色视频在线观看| 伊人久久大香线蕉亚洲五| 1024视频免费在线观看| 亚洲av成人一区二区三| 国产成年人精品一区二区| 精品午夜福利视频在线观看一区| 每晚都被弄得嗷嗷叫到高潮| 亚洲五月天丁香| 人人妻人人看人人澡| 少妇被粗大的猛进出69影院| 国产区一区二久久| 精品一区二区三区视频在线观看免费| 性色av乱码一区二区三区2| 国内精品久久久久精免费| 午夜a级毛片| 免费在线观看亚洲国产| 妹子高潮喷水视频| cao死你这个sao货| 99久久99久久久精品蜜桃| 亚洲av美国av| 亚洲精品国产一区二区精华液| 久久欧美精品欧美久久欧美| 禁无遮挡网站| 亚洲va日本ⅴa欧美va伊人久久| 波多野结衣高清作品| 国产亚洲欧美精品永久| 午夜免费成人在线视频| 后天国语完整版免费观看| 白带黄色成豆腐渣| 嫩草影院精品99| 黄色a级毛片大全视频| 一区福利在线观看| 国产亚洲欧美98| 成人av一区二区三区在线看| 视频在线观看一区二区三区| 超碰成人久久| 成人免费观看视频高清| 黄片播放在线免费| 国产精品香港三级国产av潘金莲| 在线看三级毛片| 亚洲精品国产区一区二| 老鸭窝网址在线观看| av电影中文网址| 日韩三级视频一区二区三区| 国产片内射在线| 在线视频色国产色| 久久久精品国产亚洲av高清涩受| 黄网站色视频无遮挡免费观看| 国产成人精品久久二区二区91| 精品高清国产在线一区| 中文字幕人妻丝袜一区二区| 丰满人妻熟妇乱又伦精品不卡| 久久国产乱子伦精品免费另类| 欧美日韩福利视频一区二区| 国产亚洲精品久久久久久毛片| 搡老妇女老女人老熟妇| netflix在线观看网站| 亚洲国产欧洲综合997久久, | 国产一区二区激情短视频| 伦理电影免费视频| www日本黄色视频网| 国内精品久久久久久久电影| 免费在线观看影片大全网站| 高清毛片免费观看视频网站| 成年版毛片免费区| 免费看美女性在线毛片视频| 丁香六月欧美| 88av欧美| www.熟女人妻精品国产| 欧美在线黄色| 欧美黄色淫秽网站| 神马国产精品三级电影在线观看 | 精品国产乱子伦一区二区三区| 草草在线视频免费看| 国产精品亚洲av一区麻豆| 大型黄色视频在线免费观看| 欧美乱妇无乱码| 少妇 在线观看| 亚洲aⅴ乱码一区二区在线播放 | 色在线成人网| 很黄的视频免费| 国产精品综合久久久久久久免费| 久久性视频一级片| 少妇 在线观看| 日韩av在线大香蕉| 少妇熟女aⅴ在线视频| 叶爱在线成人免费视频播放| 亚洲av中文字字幕乱码综合 | 99国产极品粉嫩在线观看| 麻豆国产av国片精品| 久久香蕉精品热| 亚洲五月天丁香| 亚洲一区中文字幕在线| 免费女性裸体啪啪无遮挡网站| 久久精品国产综合久久久| 欧美激情极品国产一区二区三区| 啦啦啦观看免费观看视频高清| 欧美午夜高清在线| www.熟女人妻精品国产| 草草在线视频免费看| 成人亚洲精品av一区二区| 搡老岳熟女国产| 欧美丝袜亚洲另类 | 亚洲国产精品合色在线| 99国产综合亚洲精品| 热re99久久国产66热| 亚洲精品美女久久av网站| 观看免费一级毛片| 亚洲激情在线av| 亚洲第一电影网av| 一区福利在线观看| 美女扒开内裤让男人捅视频| 美女免费视频网站| 国产成人欧美在线观看| 波多野结衣av一区二区av| 99久久精品国产亚洲精品| 亚洲av成人不卡在线观看播放网| 久99久视频精品免费| 亚洲国产欧美日韩在线播放| 国产熟女午夜一区二区三区| 夜夜看夜夜爽夜夜摸| 男女视频在线观看网站免费 | 久热这里只有精品99| 亚洲男人的天堂狠狠| 一a级毛片在线观看| 精品午夜福利视频在线观看一区| 精华霜和精华液先用哪个| 欧美中文日本在线观看视频| 午夜免费成人在线视频| av片东京热男人的天堂| 亚洲美女黄片视频| 亚洲成国产人片在线观看| 欧美日本亚洲视频在线播放| 97超级碰碰碰精品色视频在线观看| 婷婷精品国产亚洲av在线| 精品一区二区三区四区五区乱码| 此物有八面人人有两片| 午夜激情福利司机影院| 成人国产一区最新在线观看| 非洲黑人性xxxx精品又粗又长| 老司机靠b影院| 欧美人与性动交α欧美精品济南到| 亚洲精品在线美女| 午夜亚洲福利在线播放| 亚洲欧美日韩高清在线视频| 天堂动漫精品| 久久伊人香网站| 韩国精品一区二区三区| 99久久国产精品久久久| 日韩精品青青久久久久久| 无人区码免费观看不卡| 神马国产精品三级电影在线观看 | av电影中文网址| 免费观看人在逋| 99国产精品一区二区三区| 一夜夜www| 国产精品乱码一区二三区的特点| cao死你这个sao货| 变态另类成人亚洲欧美熟女| 久久精品成人免费网站| 午夜两性在线视频| 国产av又大| 老熟妇仑乱视频hdxx| 国产精品美女特级片免费视频播放器 | 欧美久久黑人一区二区| 深夜精品福利| 人妻丰满熟妇av一区二区三区| 亚洲国产欧洲综合997久久, | 久久精品国产亚洲av高清一级| 日韩精品中文字幕看吧| 亚洲精华国产精华精| 国产精品亚洲一级av第二区| 亚洲成av人片免费观看| 午夜视频精品福利| 国产国语露脸激情在线看| 韩国精品一区二区三区| 首页视频小说图片口味搜索| 搡老岳熟女国产| 99热这里只有精品一区 | 欧美精品亚洲一区二区| 又大又爽又粗| 色哟哟哟哟哟哟| 国产精品美女特级片免费视频播放器 | 成人欧美大片| 午夜福利免费观看在线| 国产免费男女视频| 日韩av在线大香蕉| 亚洲熟女毛片儿| 午夜精品在线福利| 亚洲全国av大片| 免费在线观看成人毛片| 18禁黄网站禁片免费观看直播| 欧美中文综合在线视频| 亚洲,欧美精品.| 99精品在免费线老司机午夜| 亚洲性夜色夜夜综合| 国产免费男女视频| 亚洲精华国产精华精| 久久精品91蜜桃| 精品熟女少妇八av免费久了| 男女午夜视频在线观看| 美女扒开内裤让男人捅视频| 免费在线观看完整版高清| 一级毛片精品| 欧美日韩乱码在线| 国产精品免费视频内射| 亚洲色图 男人天堂 中文字幕| 自线自在国产av| 久久久国产成人精品二区| av在线播放免费不卡| 国产人伦9x9x在线观看| 亚洲av成人av| 亚洲国产欧美网| 亚洲七黄色美女视频| 国产蜜桃级精品一区二区三区| 亚洲人成伊人成综合网2020| 男女下面进入的视频免费午夜 | 久久狼人影院| 精品第一国产精品| 欧美性猛交黑人性爽| 国产免费男女视频| 这个男人来自地球电影免费观看| 国产爱豆传媒在线观看 | 韩国av一区二区三区四区| 国产精品精品国产色婷婷| 日韩精品中文字幕看吧| 久久久久九九精品影院| 久久欧美精品欧美久久欧美| 亚洲男人的天堂狠狠| 久久久久久久久久黄片| 亚洲国产欧美网| av视频在线观看入口| 一级黄色大片毛片| 别揉我奶头~嗯~啊~动态视频| 午夜福利在线在线| 亚洲男人的天堂狠狠| 亚洲成av片中文字幕在线观看| 色播在线永久视频| 夜夜爽天天搞| 国产高清视频在线播放一区| 十分钟在线观看高清视频www| 男女做爰动态图高潮gif福利片| 国内精品久久久久精免费| 日韩视频一区二区在线观看| 精品福利观看| 在线观看免费午夜福利视频| 啦啦啦免费观看视频1| 国产aⅴ精品一区二区三区波| 久久久久久国产a免费观看| 日韩国内少妇激情av| 欧美成人一区二区免费高清观看 | 在线观看66精品国产| 1024手机看黄色片| 成人三级黄色视频| 91大片在线观看| 亚洲中文av在线| 少妇熟女aⅴ在线视频| 大香蕉久久成人网| 国产三级在线视频| 女人被狂操c到高潮| 精品国产美女av久久久久小说| 少妇的丰满在线观看| 1024香蕉在线观看| 精品国产国语对白av| 一级片免费观看大全| 国产高清视频在线播放一区| 婷婷丁香在线五月| 长腿黑丝高跟| www日本在线高清视频| 成年人黄色毛片网站| 亚洲精品一区av在线观看| 狂野欧美激情性xxxx| 久久国产精品男人的天堂亚洲| 欧美日韩亚洲国产一区二区在线观看| 在线观看一区二区三区| 欧美成人免费av一区二区三区| 色在线成人网| 久久久精品国产亚洲av高清涩受| 国产精品影院久久| 午夜老司机福利片| 亚洲av成人不卡在线观看播放网| aaaaa片日本免费| 中文字幕人成人乱码亚洲影| 国产精品永久免费网站| 国产精品,欧美在线| 正在播放国产对白刺激| 久久中文看片网| 黑人欧美特级aaaaaa片| 波多野结衣高清作品| 色播在线永久视频| 免费女性裸体啪啪无遮挡网站| 日日摸夜夜添夜夜添小说| 国产精品自产拍在线观看55亚洲| 天天一区二区日本电影三级| 母亲3免费完整高清在线观看| 国产精品一区二区精品视频观看| 一边摸一边抽搐一进一小说| 精品少妇一区二区三区视频日本电影| 亚洲国产精品成人综合色| 国产又黄又爽又无遮挡在线| 此物有八面人人有两片| 麻豆成人午夜福利视频| 少妇粗大呻吟视频| 欧美色欧美亚洲另类二区| 日韩欧美国产一区二区入口| 日韩一卡2卡3卡4卡2021年| 搡老岳熟女国产| 中文字幕最新亚洲高清| 日韩三级视频一区二区三区| 搡老熟女国产l中国老女人| 日韩av在线大香蕉| 亚洲中文字幕一区二区三区有码在线看 | 欧美精品亚洲一区二区| 又大又爽又粗| 亚洲精品中文字幕在线视频| 国产精品美女特级片免费视频播放器 | 久久精品国产亚洲av高清一级| 免费无遮挡裸体视频| 两性夫妻黄色片| 黄色成人免费大全| 母亲3免费完整高清在线观看| 精品国产亚洲在线| 国产精品美女特级片免费视频播放器 | 99精品久久久久人妻精品| 国产av在哪里看| 夜夜夜夜夜久久久久| 欧美在线黄色| 中文在线观看免费www的网站 | 国产三级黄色录像| 成人精品一区二区免费| 亚洲无线在线观看| 91九色精品人成在线观看| 久久中文字幕一级| 男男h啪啪无遮挡| 俄罗斯特黄特色一大片| av电影中文网址| 色播亚洲综合网| 身体一侧抽搐| 热99re8久久精品国产| 精华霜和精华液先用哪个| 久久人妻av系列| 久久久久久国产a免费观看| av电影中文网址| 亚洲精华国产精华精| 精品一区二区三区视频在线观看免费| 美女 人体艺术 gogo| 亚洲精品在线观看二区| 欧美丝袜亚洲另类 | 精品久久久久久久久久久久久 | 欧美乱妇无乱码| 成人特级黄色片久久久久久久| 热99re8久久精品国产| 亚洲国产看品久久| 最新美女视频免费是黄的| 宅男免费午夜| 看免费av毛片| 欧美性猛交╳xxx乱大交人| 亚洲国产精品合色在线| 亚洲一码二码三码区别大吗| 在线观看www视频免费| 国内精品久久久久久久电影| 巨乳人妻的诱惑在线观看| 正在播放国产对白刺激| 天堂动漫精品| 国产亚洲精品av在线| 人成视频在线观看免费观看| 久久国产精品影院| 亚洲人成电影免费在线| 欧美日本亚洲视频在线播放| 国产精品久久电影中文字幕| 国产精品久久久久久精品电影 | 欧美性猛交黑人性爽| 欧美日韩瑟瑟在线播放| 18美女黄网站色大片免费观看| 不卡一级毛片| 人人妻人人澡人人看| www日本黄色视频网| 一级毛片女人18水好多|