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

    基于全局敏感性分析和PE-HMGE的軸承潤(rùn)滑特性優(yōu)化*

    2019-01-29 06:46:42李林潔韋靜思何志良占文鋒
    汽車工程 2018年12期
    關(guān)鍵詞:曲軸全局敏感性

    李林潔,韋靜思,何志良,陳 宇,占文鋒

    (廣州汽車集團(tuán)有限公司汽車工程研究院,廣州 511434)

    前言

    摩擦損失是發(fā)動(dòng)機(jī)能量損耗的主要原因之一,降低發(fā)動(dòng)機(jī)摩擦損失對(duì)汽車節(jié)能減排具有重要意義。其中,軸承摩擦副是發(fā)動(dòng)機(jī)的重要組成部分,其潤(rùn)滑性能的好壞直接影響發(fā)動(dòng)機(jī)的性能及使用壽命。近年來,采用彈性流體動(dòng)力學(xué)(elasto-hydrodynamic,EHD)對(duì)軸承潤(rùn)滑特性進(jìn)行分析的方法逐漸代替了傳統(tǒng)經(jīng)典流體動(dòng)力潤(rùn)滑方法,成為主流。申立中等[1]對(duì)柴油機(jī)連桿大頭軸承進(jìn)行EHD模擬計(jì)算,重點(diǎn)研究了軸承相對(duì)平均間隙、曲柄銷油孔直徑和軸瓦寬度對(duì)大頭軸承潤(rùn)滑特性的影響。

    隨著油耗和排放法規(guī)的日益嚴(yán)格,對(duì)于軸承的要求不能僅局限于其承載和潤(rùn)滑能力滿足需求,而是在上述前提下,盡可能降低摩擦功耗和消耗的機(jī)油量,以提高發(fā)動(dòng)機(jī)有效功率,改善燃油經(jīng)濟(jì)性。這需要對(duì)軸承進(jìn)行以多目標(biāo)優(yōu)化為基礎(chǔ)的分析與設(shè)計(jì)。文獻(xiàn)中也有這方面相應(yīng)的研究成果,張俊紅等[2]采用極差分析的方法確定了影響最小油膜厚度、最大油膜壓力和摩擦損失功率因素的主次關(guān)系,利用BP神經(jīng)網(wǎng)絡(luò)理論建立了主軸承潤(rùn)滑特性模型,并對(duì)影響潤(rùn)滑特性的主要參數(shù)進(jìn)行了優(yōu)化。

    由于EHD采用有限元法求解軸頸、軸瓦和軸承壁的彈性變形,用有限差分法求解雷諾方程,整個(gè)過程涉及有限單元網(wǎng)格與有限差分網(wǎng)格之間的相互映射,因此計(jì)算成本過高,很難滿足優(yōu)化設(shè)計(jì)的計(jì)算要求。所以,建立準(zhǔn)確表征軸承潤(rùn)滑特性數(shù)值模型是對(duì)軸承進(jìn)行準(zhǔn)確多目標(biāo)優(yōu)化設(shè)計(jì)的關(guān)鍵。本文中從多體動(dòng)力學(xué)角度對(duì)某汽油發(fā)動(dòng)機(jī)主軸承進(jìn)行EHD分析,并通過拆除倒拖摩擦功試驗(yàn)方法校驗(yàn)仿真模型。針對(duì)主軸承潤(rùn)滑系統(tǒng)非線性的特點(diǎn),利于基于方差的全局敏感性分析方法——Sobol法[3],定量分析輸入變量對(duì)系統(tǒng)輸出的影響程度,確定敏感性參數(shù)。運(yùn)用克里格方法(Kriging)[4]建立可以準(zhǔn)確表征主軸承潤(rùn)滑特性的數(shù)學(xué)模型,采用具有良好的探索性能的全局多目標(biāo)梯度算法——PE-HMGE優(yōu)化算法[5],計(jì)算目標(biāo)函數(shù)Pareto前沿個(gè)數(shù),并重構(gòu)Kriging模型,從而實(shí)現(xiàn)主軸承潤(rùn)滑特性多目標(biāo)優(yōu)化設(shè)計(jì)。

    1 理論介紹

    1.1 全局敏感性分析理論

    基于方差的Sobol全局敏感性分析[3]是一種定量分析方法,其核心思想是將模型分解為單個(gè)參數(shù)及參數(shù)之間相互組合的函數(shù),采用蒙特卡洛估計(jì)[6]得到參數(shù)的各階次敏感度,具體方法如下。

    假設(shè)系統(tǒng)的設(shè)計(jì)參數(shù)個(gè)數(shù)為k,在所有設(shè)計(jì)參數(shù)的初始設(shè)計(jì)空間內(nèi)隨機(jī)抽樣n次,共抽取A,B兩個(gè)矩陣,如式(1)所示,其中,矩陣的每一行表示k個(gè)參數(shù)的組合。

    將式(1)中矩陣A的第i列換成矩陣B的第i列,其余列保持不變,得矩陣記為Ci;再將矩陣B的第i列換成矩陣A的第i列,其余列保持不變,得到矩陣記為 C-i,如式(2)所示:

    將每組設(shè)計(jì)參數(shù)分別代入到系統(tǒng)模型中進(jìn)行計(jì)算,可得到對(duì)應(yīng)每組設(shè)計(jì)參數(shù)的輸出響應(yīng)值,則系統(tǒng)響應(yīng)的方差估計(jì)和參數(shù)的敏感性指標(biāo)計(jì)算如式(3)~(6)所示:

    輸入?yún)?shù) xi的主效應(yīng)指標(biāo)的估計(jì)為

    輸入?yún)?shù) xi的全效應(yīng)指標(biāo)的估計(jì)為

    1.2 K riging模型及優(yōu)化方法

    Kriging模型采用的是最優(yōu)內(nèi)插法,是一種基于統(tǒng)計(jì)理論的估計(jì)方差最小的無偏估計(jì)模型,其數(shù)學(xué)表達(dá)式為

    式中:β為回歸系數(shù);fT(X)為回歸模型;z(X)為隨機(jī)分布誤差。

    在初始設(shè)計(jì)變量及設(shè)計(jì)空間S確定的條件下,采用優(yōu)化的拉丁方抽樣方法對(duì)設(shè)計(jì)空間進(jìn)行抽樣。

    式中:n為設(shè)計(jì)參數(shù)的個(gè)數(shù);N1為構(gòu)造Kriging代理模型所需的最少點(diǎn)。

    利用初始樣本點(diǎn)和響應(yīng)值建立Kriging初始模型。為保證初始模型在局部的精度,優(yōu)化迭代過程中,模型應(yīng)不斷使用新增樣本點(diǎn)進(jìn)行更新,新增樣本點(diǎn)個(gè)數(shù)的確定方法如下所述。

    引入適應(yīng)度函數(shù)[8]Gi,其表達(dá)式為

    式中:Gi為第i個(gè)點(diǎn)的適應(yīng)度為第 i個(gè)點(diǎn)第 k個(gè)目標(biāo)值的歸一化值,k=1,2,…,m;l為前沿指數(shù),這里取值為1。當(dāng)Pareto前沿點(diǎn)密集和均勻分布時(shí),所有的Pareto點(diǎn)的適應(yīng)度值都趨近1。

    首次迭代時(shí),利用式(11)計(jì)算得出的Pareto前沿點(diǎn)個(gè)數(shù),記為M1;后續(xù)迭代時(shí),提取第一次樣本點(diǎn)和本次迭代之前所有輪迭代的新增樣本點(diǎn),重構(gòu)Kriging模型,求解Pareto解個(gè)數(shù),記為Mq,其中q為迭代次數(shù)。將M1和Mq合并在一起,再進(jìn)行適應(yīng)度計(jì)算,重新確定Pareto前沿點(diǎn)個(gè)數(shù),記為Qq。新增樣本點(diǎn)數(shù)確定之后,確定新增樣本點(diǎn)的設(shè)計(jì)空間,其確定方法如下所述。

    初步確定新增樣本點(diǎn)設(shè)計(jì)空間大小。設(shè)第q次迭代新增樣本點(diǎn)空間為

    在新的設(shè)計(jì)空間內(nèi),采用具有良好的探索性能的全局多目標(biāo)梯度算法——PE-HMGE優(yōu)化算法[5],選擇接近Pareto前沿的個(gè)體,將第q次迭代中適應(yīng)度值最接近1的那組數(shù)據(jù)作為第k個(gè)目標(biāo)的最優(yōu)解xkq*,將其代入到分析模型中得到最優(yōu)目標(biāo)值f(xkq*)。

    迭代過程中,判斷第k個(gè)目標(biāo)在第q次迭代和第q-1次迭代產(chǎn)生的最優(yōu)目標(biāo)函數(shù)值的相對(duì)誤差是否滿足收斂值,若滿足,則f(xkq*)即為最優(yōu)值,優(yōu)化流程結(jié)束;如不滿足則繼續(xù)迭代,直到滿足收斂準(zhǔn)則,整個(gè)分析過程的流程如圖1所示。

    圖1 多參數(shù)全局敏感性分析和優(yōu)化流程

    2 模型建立

    2.1 曲軸系 -缸體多體動(dòng)力學(xué)模型

    本文中運(yùn)用AVL的Excite Power Unit平臺(tái)建立曲軸系 -缸體多體動(dòng)力學(xué)模型,主要包括了曲軸、連桿、缸體、軸瓦等部件模塊,其中主軸承主要設(shè)計(jì)參數(shù)如表1所示。

    表1 主軸承主要參數(shù)

    2.2 主軸承EHD仿真結(jié)果分析

    為更準(zhǔn)確計(jì)算軸承的潤(rùn)滑特性,主軸承采用EHD(elastic hydro-dynamic)軸承模型。本文軸承潤(rùn)滑特性分析中,采用峰值油膜壓力(peak oil film pressure,POFP)和最小油膜厚度(minimum oil film thickness,MOFT)表示軸承的承載能力,峰值粗糙接觸壓力(peak asperity contact pressure,PASP)和總摩擦損失(total friction power loss,TFPL)表示軸承的潤(rùn)滑性能。

    POFP是軸承油膜壓力的最大值,一般POFP限值為:Plimit=HB/0.102,其中HB為材料的布氏硬度;若軸瓦和軸頸的表面粗糙度均方差分別為σ1和σ2,則潤(rùn)滑表面的合成均方差即峰值粗糙度為:σa=建立流體動(dòng)壓潤(rùn)滑的最低要求是MOFT不小于σa;PASP越小,表明潤(rùn)滑性能越好,通常不大于100MPa;TFPL是液動(dòng)摩擦損失和粗糙摩擦損失之和。

    2.3 曲軸摩擦功試驗(yàn)

    本文中采用倒拖法分解摩擦功,發(fā)動(dòng)機(jī)僅裝配曲軸,去除其他附件系統(tǒng),在不影響曲軸平衡率的前提下,堵塞油道。測(cè)試臺(tái)架外界冷卻液和機(jī)油溫度及壓力控制系統(tǒng),調(diào)節(jié)冷卻液和機(jī)油溫度恒定,同時(shí)保證主油道油壓、進(jìn)出水壓與發(fā)動(dòng)機(jī)正常運(yùn)行時(shí)一致,調(diào)節(jié)轉(zhuǎn)速不點(diǎn)火倒拖,測(cè)量轉(zhuǎn)矩值。圖2為曲軸摩擦功測(cè)試結(jié)果和曲軸EHD仿真結(jié)果對(duì)比,試驗(yàn)和仿真結(jié)果在全轉(zhuǎn)速工況下變化趨勢(shì)基本一致,利用試驗(yàn)結(jié)果校驗(yàn)仿真模型,通過仿真手段進(jìn)行參數(shù)研究,優(yōu)化結(jié)構(gòu)設(shè)計(jì)。

    圖2 曲軸主軸承摩擦損失

    2.4 基于方差的全局敏感性分析

    以第2主軸承為例,對(duì)發(fā)動(dòng)機(jī)主軸承各參數(shù)進(jìn)行敏感度分析。主軸承可考慮的參數(shù)包括軸承寬度、軸承間隙、軸瓦磨損量、軸頸鼓型量、機(jī)油動(dòng)力黏度、供油壓力等,如果全部參與計(jì)算,設(shè)計(jì)空間維數(shù)高達(dá)12維,必須增加采樣點(diǎn)數(shù)來描述設(shè)計(jì)空間的信息,從而增加計(jì)算規(guī)模。通過靈敏度分析,摒棄非敏感參數(shù),精簡(jiǎn)模型,降低計(jì)算規(guī)模。采用蒙特卡羅方法抽取104個(gè)樣本點(diǎn),分別代表A,B兩個(gè)矩陣,通過Sobol全局敏感性分析計(jì)算各參數(shù)的主效應(yīng),結(jié)果如表2所示。

    表2 全局敏感性分析結(jié)果

    2.5 主軸承K riging潤(rùn)滑特性數(shù)值模型

    根據(jù)全局敏感性分析的結(jié)果,選擇軸承寬度、軸承間隙、軸瓦磨損量、軸頸油孔直徑和軸頸油孔位置5個(gè)參數(shù)為設(shè)計(jì)變量。數(shù)值模型如下,

    max MOFT

    min POFP,PASP,TFPL

    s.t.MOFT≥0.8μm

    POFP≤200MPa

    PASP≤60MPa

    應(yīng)用前文中提出的基于PE-HMGE優(yōu)化算法的模型優(yōu)化方法,重構(gòu)Kriging潤(rùn)滑特性數(shù)值模型。隨機(jī)選取5組設(shè)計(jì)變量,對(duì)比仿真模型和優(yōu)化模型計(jì)算結(jié)果,模型精度如表3~表6所示。

    表3 POFP模型

    表4 MOFT模型

    表5 PASP模型

    從上面列表中可以看出,各優(yōu)化模型的預(yù)測(cè)值與實(shí)際仿真結(jié)果之間的相對(duì)誤差較小,說明Kriging方法在較少的采樣點(diǎn)的情況下,可以建立較為準(zhǔn)確的模型。在優(yōu)化迭代的過程中,數(shù)值模型最優(yōu)解附近的樣本點(diǎn)不斷地增加,使得最優(yōu)解附近的精度得到了很大的提高,因此得到的解更加具有可信性,與理論最優(yōu)解更接近。

    表6 TFPL模型

    基于優(yōu)化模型得到的最優(yōu)化設(shè)計(jì)結(jié)果如表7所示。

    表7 最優(yōu)設(shè)計(jì)結(jié)果

    將優(yōu)化后的設(shè)計(jì)變量輸入到發(fā)動(dòng)機(jī)軸承的實(shí)際物理模型中進(jìn)行計(jì)算,以第2主軸承為例,優(yōu)化前后軸承潤(rùn)滑特性如圖3所示。

    從圖3中可以看出,經(jīng)過優(yōu)化后,軸承的最大峰值油膜壓力(POPF)由原機(jī)的 79.6提高至85.2MPa,但仍滿足軸瓦的設(shè)計(jì)標(biāo)準(zhǔn);最小油膜厚度(MOFT)由原機(jī)的 0.77提高到 1.07μm,提高了39%;最大峰值粗糙接觸壓力(PASP)由原機(jī)的65.8降低至36.2MPa,降低了45%;通過計(jì)算得到的平均有效摩擦損失由原機(jī)的242.8降低至195.6W,降低了19%,軸承的潤(rùn)滑特性得到明顯改善。全轉(zhuǎn)速工況下,曲軸主軸承的總摩擦損失也有明顯的下降,如圖4所示。

    圖3 軸承潤(rùn)滑特性優(yōu)化前后對(duì)比

    將優(yōu)化方案應(yīng)用于實(shí)際樣機(jī)中,進(jìn)行曲軸摩擦損失測(cè)試,結(jié)果如圖5所示。在測(cè)試轉(zhuǎn)速范圍內(nèi),摩擦損失降低了17.6%,與仿真結(jié)果相符。

    3 結(jié)論

    圖4 曲軸主軸承摩擦損失仿真結(jié)果

    圖5 曲軸主軸承摩擦損失試驗(yàn)結(jié)果

    (1)本文中將全局敏感性分析和PE-HMGE優(yōu)化算法相結(jié)合,研究復(fù)雜非線性系統(tǒng)多參數(shù)的優(yōu)化求解問題,提出了多參數(shù)多目標(biāo)優(yōu)化策略,并將該方法應(yīng)用到軸承的潤(rùn)滑特性優(yōu)化設(shè)計(jì)中。

    (2)通過全局敏感性分析方法,定量確定了每個(gè)設(shè)計(jì)參數(shù)對(duì)潤(rùn)滑特性的影響,從而篩選出重要參數(shù),大大減少了計(jì)算量。

    (3)本文中提出的模型優(yōu)化方法在軸承設(shè)計(jì)中取得良好的優(yōu)化效果。對(duì)比仿真結(jié)果,最小油膜厚度提高了39%;最大峰值粗糙接觸壓力降低了45%;平均有效摩擦損失降低了19%。測(cè)試結(jié)果也表明軸承的潤(rùn)滑特性得到明顯改善,證明了該研究方法的有效性。

    猜你喜歡
    曲軸全局敏感性
    Cahn-Hilliard-Brinkman系統(tǒng)的全局吸引子
    量子Navier-Stokes方程弱解的全局存在性
    淺析曲軸磨床的可靠性改進(jìn)
    一種曲軸豎置熱處理裝置設(shè)計(jì)
    落子山東,意在全局
    金橋(2018年4期)2018-09-26 02:24:54
    釔對(duì)Mg-Zn-Y-Zr合金熱裂敏感性影響
    基于支持向量機(jī)的曲軸動(dòng)平衡控制方法
    AH70DB鋼焊接熱影響區(qū)組織及其冷裂敏感性
    焊接(2016年1期)2016-02-27 12:55:37
    如何培養(yǎng)和提高新聞敏感性
    新聞傳播(2015年8期)2015-07-18 11:08:24
    微小RNA與食管癌放射敏感性的相關(guān)研究
    免费看a级黄色片| 丰满人妻熟妇乱又伦精品不卡| 欧美bdsm另类| 91字幕亚洲| 搞女人的毛片| 精品一区二区三区av网在线观看| 国产精品永久免费网站| 午夜久久久久精精品| 国产精品电影一区二区三区| www.熟女人妻精品国产| 两人在一起打扑克的视频| 国内毛片毛片毛片毛片毛片| 波多野结衣高清作品| 国产三级黄色录像| 午夜福利免费观看在线| 天堂动漫精品| 亚洲成人精品中文字幕电影| 网址你懂的国产日韩在线| 级片在线观看| 国产在视频线在精品| 真实男女啪啪啪动态图| 真实男女啪啪啪动态图| 国产av麻豆久久久久久久| netflix在线观看网站| 欧美日韩国产亚洲二区| 变态另类成人亚洲欧美熟女| 99久久久亚洲精品蜜臀av| 人人妻,人人澡人人爽秒播| 九九久久精品国产亚洲av麻豆| 久久久久久国产a免费观看| 嫩草影院精品99| 久久久色成人| 麻豆一二三区av精品| 国产高潮美女av| 久久精品久久久久久噜噜老黄 | 18美女黄网站色大片免费观看| 老司机福利观看| 亚洲,欧美,日韩| 日本五十路高清| 国产一区二区三区在线臀色熟女| 蜜桃亚洲精品一区二区三区| 伦理电影大哥的女人| 18禁黄网站禁片免费观看直播| 99国产综合亚洲精品| 日本黄色片子视频| 少妇熟女aⅴ在线视频| 精品久久久久久,| 舔av片在线| 在线免费观看不下载黄p国产 | 精品人妻偷拍中文字幕| 亚洲av电影在线进入| 51国产日韩欧美| 成人欧美大片| 国产精品久久久久久亚洲av鲁大| 国产一区二区在线av高清观看| 真人一进一出gif抽搐免费| 又黄又爽又刺激的免费视频.| 国产老妇女一区| 欧美午夜高清在线| 真人一进一出gif抽搐免费| 精品久久久久久久久久免费视频| 两个人的视频大全免费| 一个人免费在线观看的高清视频| 亚洲av不卡在线观看| 亚洲va日本ⅴa欧美va伊人久久| 搡老岳熟女国产| 欧美成人性av电影在线观看| 国产在视频线在精品| 欧美日韩综合久久久久久 | 国产一区二区亚洲精品在线观看| 桃色一区二区三区在线观看| 黄色一级大片看看| 高清日韩中文字幕在线| 亚洲国产色片| 搡老熟女国产l中国老女人| 最好的美女福利视频网| 91久久精品国产一区二区成人| 亚洲人成网站在线播| 人妻制服诱惑在线中文字幕| 香蕉av资源在线| 国产亚洲精品久久久com| 真实男女啪啪啪动态图| 欧美又色又爽又黄视频| 丝袜美腿在线中文| 亚洲欧美日韩高清在线视频| 在线观看av片永久免费下载| 午夜亚洲福利在线播放| 国产av一区在线观看免费| 欧美区成人在线视频| 美女高潮喷水抽搐中文字幕| 99热6这里只有精品| 国产亚洲精品综合一区在线观看| 婷婷精品国产亚洲av在线| 黄色视频,在线免费观看| 国产成人欧美在线观看| 熟女电影av网| 90打野战视频偷拍视频| 国产精品伦人一区二区| 99视频精品全部免费 在线| 免费电影在线观看免费观看| 能在线免费观看的黄片| 人妻制服诱惑在线中文字幕| 国产成年人精品一区二区| 在现免费观看毛片| 久久久久久久久中文| 免费看a级黄色片| 两个人视频免费观看高清| 欧洲精品卡2卡3卡4卡5卡区| 国内精品美女久久久久久| 亚洲熟妇熟女久久| 一级av片app| 国产黄片美女视频| 精品久久久久久久久av| 老司机福利观看| 麻豆成人av在线观看| 99热只有精品国产| 国产日本99.免费观看| 好看av亚洲va欧美ⅴa在| 国产不卡一卡二| 一边摸一边抽搐一进一小说| 人人妻人人看人人澡| 亚洲乱码一区二区免费版| 精品免费久久久久久久清纯| 亚洲人成网站高清观看| 2021天堂中文幕一二区在线观| 制服丝袜大香蕉在线| 一卡2卡三卡四卡精品乱码亚洲| 亚洲经典国产精华液单 | 91九色精品人成在线观看| 婷婷色综合大香蕉| 久久天躁狠狠躁夜夜2o2o| 又黄又爽又刺激的免费视频.| 成人毛片a级毛片在线播放| 黄色视频,在线免费观看| 亚洲avbb在线观看| 日韩人妻高清精品专区| 亚洲精华国产精华精| 男女那种视频在线观看| 人人妻人人澡欧美一区二区| 日韩欧美精品免费久久 | 一个人免费在线观看电影| 特大巨黑吊av在线直播| 麻豆久久精品国产亚洲av| a级毛片a级免费在线| 国产单亲对白刺激| 国产极品精品免费视频能看的| 欧美一区二区国产精品久久精品| 99热这里只有精品一区| 国产高清视频在线播放一区| 亚洲最大成人av| 一卡2卡三卡四卡精品乱码亚洲| 夜夜看夜夜爽夜夜摸| 欧美性猛交黑人性爽| 99久国产av精品| 国产老妇女一区| 女生性感内裤真人,穿戴方法视频| 亚洲一区高清亚洲精品| 不卡一级毛片| 欧美激情在线99| 亚洲人成网站高清观看| 国产av在哪里看| 国产精品久久久久久人妻精品电影| 又黄又爽又刺激的免费视频.| 97碰自拍视频| 桃红色精品国产亚洲av| 国产成+人综合+亚洲专区| 草草在线视频免费看| 午夜精品久久久久久毛片777| 国产熟女xx| 国产高清激情床上av| 成年人黄色毛片网站| 亚洲午夜理论影院| 国产精品免费一区二区三区在线| 国产单亲对白刺激| 欧美一区二区精品小视频在线| 亚洲av一区综合| www.熟女人妻精品国产| 村上凉子中文字幕在线| 午夜精品在线福利| 欧美不卡视频在线免费观看| 嫩草影院精品99| 国产精品永久免费网站| 午夜影院日韩av| 国产精品久久视频播放| 国产在线男女| 国产精品一及| 免费在线观看亚洲国产| av福利片在线观看| 日韩欧美免费精品| 国产av麻豆久久久久久久| 少妇被粗大猛烈的视频| 51国产日韩欧美| 网址你懂的国产日韩在线| 亚洲精品粉嫩美女一区| 别揉我奶头 嗯啊视频| 午夜福利高清视频| 国产精品久久视频播放| 一个人看视频在线观看www免费| 91麻豆精品激情在线观看国产| 亚洲成人久久爱视频| 国产精品国产高清国产av| 日本精品一区二区三区蜜桃| 成人亚洲精品av一区二区| 国产一区二区三区视频了| 我要搜黄色片| 国内精品久久久久精免费| 日韩欧美精品v在线| 免费观看的影片在线观看| 国产伦一二天堂av在线观看| 精品久久久久久久久av| 国产精品久久电影中文字幕| 日韩欧美国产一区二区入口| 十八禁人妻一区二区| 日日摸夜夜添夜夜添av毛片 | a在线观看视频网站| 亚洲成av人片在线播放无| 九九久久精品国产亚洲av麻豆| 国产成人av教育| 国产白丝娇喘喷水9色精品| 女生性感内裤真人,穿戴方法视频| 亚洲久久久久久中文字幕| 国产乱人视频| 美女cb高潮喷水在线观看| 一本综合久久免费| 亚洲av日韩精品久久久久久密| 亚洲第一电影网av| 自拍偷自拍亚洲精品老妇| 18禁裸乳无遮挡免费网站照片| 国产激情偷乱视频一区二区| 欧美最新免费一区二区三区 | 久久伊人香网站| 久久精品人妻少妇| 青草久久国产| 国语自产精品视频在线第100页| 欧美一区二区国产精品久久精品| 淫妇啪啪啪对白视频| 精品人妻视频免费看| 99在线视频只有这里精品首页| 又黄又爽又免费观看的视频| 久久久精品欧美日韩精品| 一本精品99久久精品77| 久久精品国产自在天天线| 能在线免费观看的黄片| 永久网站在线| 97热精品久久久久久| 在线观看舔阴道视频| 美女高潮喷水抽搐中文字幕| 高清在线国产一区| 国产一区二区在线观看日韩| 亚洲欧美激情综合另类| 婷婷精品国产亚洲av| 国产欧美日韩一区二区三| 久久久久国内视频| 人妻丰满熟妇av一区二区三区| av天堂在线播放| 男插女下体视频免费在线播放| 日本熟妇午夜| 免费看美女性在线毛片视频| 精品久久久久久久人妻蜜臀av| 精品一区二区免费观看| 精品午夜福利视频在线观看一区| av在线观看视频网站免费| 美女cb高潮喷水在线观看| 亚洲欧美日韩东京热| 有码 亚洲区| 一本精品99久久精品77| 黄色女人牲交| 中文字幕免费在线视频6| 99国产极品粉嫩在线观看| 午夜精品在线福利| 亚洲国产高清在线一区二区三| 日日干狠狠操夜夜爽| 赤兔流量卡办理| 亚洲,欧美精品.| 一卡2卡三卡四卡精品乱码亚洲| 午夜福利在线在线| 51国产日韩欧美| 久久精品91蜜桃| 亚洲欧美日韩卡通动漫| 亚洲va日本ⅴa欧美va伊人久久| 丰满乱子伦码专区| 日日夜夜操网爽| 精品一区二区免费观看| 亚洲美女搞黄在线观看 | 国产高清视频在线播放一区| 男人和女人高潮做爰伦理| .国产精品久久| 精品免费久久久久久久清纯| 熟妇人妻久久中文字幕3abv| 成人特级av手机在线观看| 国产高清三级在线| 久久精品夜夜夜夜夜久久蜜豆| 最新在线观看一区二区三区| 国产乱人视频| 日韩 亚洲 欧美在线| 成年女人看的毛片在线观看| 最后的刺客免费高清国语| 成人国产综合亚洲| 51午夜福利影视在线观看| 免费av观看视频| 长腿黑丝高跟| 欧美一级a爱片免费观看看| av欧美777| 色视频www国产| 亚洲 国产 在线| 男人和女人高潮做爰伦理| 在线a可以看的网站| 精品免费久久久久久久清纯| 国产伦一二天堂av在线观看| 我要搜黄色片| 熟妇人妻久久中文字幕3abv| 国产伦人伦偷精品视频| 91狼人影院| 成人特级av手机在线观看| 免费观看精品视频网站| 波多野结衣高清作品| 国产aⅴ精品一区二区三区波| 精品久久久久久久久亚洲 | 99国产精品一区二区三区| 中文字幕熟女人妻在线| 看片在线看免费视频| 欧美黑人欧美精品刺激| 十八禁国产超污无遮挡网站| 91麻豆精品激情在线观看国产| 欧美高清性xxxxhd video| 亚洲成人精品中文字幕电影| 国产单亲对白刺激| 精华霜和精华液先用哪个| 欧美色欧美亚洲另类二区| 91久久精品电影网| 麻豆成人午夜福利视频| 两个人的视频大全免费| 欧美日韩福利视频一区二区| 国产乱人视频| 国产中年淑女户外野战色| 每晚都被弄得嗷嗷叫到高潮| 麻豆国产av国片精品| av福利片在线观看| 十八禁人妻一区二区| 一级黄色大片毛片| 18美女黄网站色大片免费观看| 老熟妇乱子伦视频在线观看| 亚洲欧美日韩卡通动漫| 亚洲va日本ⅴa欧美va伊人久久| 欧美一区二区亚洲| 国产精品亚洲美女久久久| 亚洲av日韩精品久久久久久密| 窝窝影院91人妻| 国产高清视频在线观看网站| 最近最新中文字幕大全电影3| 久久99热这里只有精品18| 日韩欧美在线二视频| av国产免费在线观看| 一级av片app| 久久国产乱子免费精品| 国产v大片淫在线免费观看| 成人三级黄色视频| 亚洲av电影不卡..在线观看| 丁香欧美五月| 色av中文字幕| 国产伦精品一区二区三区视频9| 老司机福利观看| eeuss影院久久| 91麻豆av在线| avwww免费| 蜜桃久久精品国产亚洲av| 一区福利在线观看| 欧美国产日韩亚洲一区| 日韩欧美三级三区| 人妻丰满熟妇av一区二区三区| 国产精品一区二区三区四区久久| 可以在线观看毛片的网站| 国产欧美日韩精品亚洲av| 看免费av毛片| 淫秽高清视频在线观看| 欧美日韩国产亚洲二区| 女人十人毛片免费观看3o分钟| 欧美绝顶高潮抽搐喷水| 久久6这里有精品| 在线十欧美十亚洲十日本专区| 听说在线观看完整版免费高清| 欧美高清成人免费视频www| 亚洲av电影在线进入| 国产亚洲精品av在线| 国产精品女同一区二区软件 | 国产真实伦视频高清在线观看 | 黄色一级大片看看| 中国美女看黄片| 麻豆成人av在线观看| 俺也久久电影网| 午夜福利在线观看免费完整高清在 | 色尼玛亚洲综合影院| 日韩中字成人| 国内久久婷婷六月综合欲色啪| www.色视频.com| 黄色丝袜av网址大全| 我的女老师完整版在线观看| 99热这里只有是精品50| 九色成人免费人妻av| 窝窝影院91人妻| 久久久久久大精品| 欧美绝顶高潮抽搐喷水| 嫩草影院入口| 九九在线视频观看精品| 成年人黄色毛片网站| 我的老师免费观看完整版| 能在线免费观看的黄片| 村上凉子中文字幕在线| 国产精品久久电影中文字幕| 国产av在哪里看| 欧美中文日本在线观看视频| 一区福利在线观看| 十八禁国产超污无遮挡网站| 成年女人永久免费观看视频| netflix在线观看网站| 少妇裸体淫交视频免费看高清| 欧美激情国产日韩精品一区| 99国产精品一区二区蜜桃av| 久久国产精品影院| 欧美色视频一区免费| 中文字幕av成人在线电影| 99热这里只有是精品在线观看 | 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国产精品伦人一区二区| 久久99热这里只有精品18| 观看免费一级毛片| 哪里可以看免费的av片| 国产午夜精品论理片| 亚洲欧美日韩卡通动漫| 在线天堂最新版资源| 一级作爱视频免费观看| 亚洲精品在线观看二区| 欧美极品一区二区三区四区| 国产91精品成人一区二区三区| 啦啦啦观看免费观看视频高清| 国产一区二区亚洲精品在线观看| 国产精品99久久久久久久久| 男女视频在线观看网站免费| 一区二区三区激情视频| 精品人妻熟女av久视频| 亚洲国产精品999在线| 日本三级黄在线观看| 亚洲av美国av| 亚洲中文日韩欧美视频| 一区二区三区四区激情视频 | 99热6这里只有精品| 国产真实伦视频高清在线观看 | 精品不卡国产一区二区三区| 欧美色视频一区免费| 国产精品影院久久| 内地一区二区视频在线| 美女黄网站色视频| 国产伦精品一区二区三区四那| 亚洲狠狠婷婷综合久久图片| 少妇高潮的动态图| 男插女下体视频免费在线播放| 国产爱豆传媒在线观看| 成年女人毛片免费观看观看9| 久久精品国产自在天天线| 国产大屁股一区二区在线视频| 亚洲内射少妇av| 亚洲无线观看免费| 国产精品久久视频播放| 久久久久久九九精品二区国产| 两人在一起打扑克的视频| 男女下面进入的视频免费午夜| 久久国产精品人妻蜜桃| 国产一区二区在线av高清观看| 国产伦在线观看视频一区| 毛片一级片免费看久久久久 | 免费看a级黄色片| 精品99又大又爽又粗少妇毛片 | av女优亚洲男人天堂| 12—13女人毛片做爰片一| 成年人黄色毛片网站| 熟女电影av网| 在线观看午夜福利视频| 99久久九九国产精品国产免费| 最近中文字幕高清免费大全6 | 日本一本二区三区精品| 我要看日韩黄色一级片| 亚洲精品一卡2卡三卡4卡5卡| 毛片女人毛片| 欧美日本亚洲视频在线播放| 又爽又黄a免费视频| 久久久久久国产a免费观看| 免费看光身美女| 日本一二三区视频观看| 日本五十路高清| 国产久久久一区二区三区| 天堂动漫精品| 久久精品国产亚洲av涩爱 | 久久热精品热| 99久国产av精品| 欧美黄色片欧美黄色片| 久久久久免费精品人妻一区二区| 欧美成人一区二区免费高清观看| 午夜福利高清视频| 中文字幕免费在线视频6| 天堂动漫精品| 欧美日韩国产亚洲二区| 十八禁人妻一区二区| 波多野结衣高清无吗| 内地一区二区视频在线| www.999成人在线观看| 欧美乱色亚洲激情| 午夜福利视频1000在线观看| 人人妻,人人澡人人爽秒播| 欧洲精品卡2卡3卡4卡5卡区| 精品久久国产蜜桃| 欧美不卡视频在线免费观看| 国产人妻一区二区三区在| 村上凉子中文字幕在线| 男女那种视频在线观看| 亚洲av一区综合| 国产精品久久久久久久电影| 两人在一起打扑克的视频| 一进一出好大好爽视频| 少妇的逼水好多| 嫩草影院入口| 久久欧美精品欧美久久欧美| 久久精品国产亚洲av天美| 亚洲熟妇熟女久久| 国产精品国产高清国产av| 国产淫片久久久久久久久 | 国产视频内射| 国产精品,欧美在线| 欧美高清成人免费视频www| 中亚洲国语对白在线视频| 97碰自拍视频| 级片在线观看| 极品教师在线免费播放| 一区二区三区激情视频| 国产成人欧美在线观看| 亚洲国产欧洲综合997久久,| 亚洲最大成人中文| 一级作爱视频免费观看| 一区福利在线观看| 赤兔流量卡办理| 一区二区三区免费毛片| av中文乱码字幕在线| 中文字幕人妻熟人妻熟丝袜美| 90打野战视频偷拍视频| 女同久久另类99精品国产91| 日本撒尿小便嘘嘘汇集6| 日韩欧美一区二区三区在线观看| netflix在线观看网站| a在线观看视频网站| 欧美日本亚洲视频在线播放| 99热这里只有是精品50| 两个人视频免费观看高清| 国产高清视频在线播放一区| 内地一区二区视频在线| 亚洲av电影在线进入| 国产亚洲精品久久久久久毛片| aaaaa片日本免费| 久久久久久久久大av| 久久久久免费精品人妻一区二区| 青草久久国产| 少妇被粗大猛烈的视频| 国产在线男女| 俺也久久电影网| 日韩欧美精品v在线| 啦啦啦观看免费观看视频高清| 亚洲七黄色美女视频| 国产精品综合久久久久久久免费| 老司机福利观看| 琪琪午夜伦伦电影理论片6080| 三级男女做爰猛烈吃奶摸视频| x7x7x7水蜜桃| 亚洲精品色激情综合| 欧美精品啪啪一区二区三区| 国产精品美女特级片免费视频播放器| 九九久久精品国产亚洲av麻豆| 91在线观看av| 久久午夜亚洲精品久久| 午夜福利欧美成人| 精品久久久久久久人妻蜜臀av| 午夜老司机福利剧场| 精品人妻视频免费看| 欧美潮喷喷水| 日本 欧美在线| 日韩欧美国产一区二区入口| 九色国产91popny在线| 69人妻影院| 亚洲专区中文字幕在线| 乱码一卡2卡4卡精品| 亚洲av第一区精品v没综合| 久久久久久九九精品二区国产| www.熟女人妻精品国产| 免费观看人在逋| 人妻夜夜爽99麻豆av| 中文字幕免费在线视频6| 国产成人福利小说| 激情在线观看视频在线高清| bbb黄色大片| 亚洲,欧美精品.| 黄色一级大片看看| 免费在线观看亚洲国产| 丰满人妻一区二区三区视频av| 欧美精品啪啪一区二区三区| 在线免费观看不下载黄p国产 | 亚洲av免费高清在线观看| 少妇被粗大猛烈的视频| 免费av观看视频| 国产精品嫩草影院av在线观看 | 长腿黑丝高跟| 亚洲专区国产一区二区| 麻豆av噜噜一区二区三区| 床上黄色一级片| 三级男女做爰猛烈吃奶摸视频| 国产亚洲欧美98| 久久精品影院6| 99久久九九国产精品国产免费| 51国产日韩欧美|