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

    基于Kriging模型的快速拆裝葉輪輪盤(pán)多目標(biāo)優(yōu)化設(shè)計(jì)

    2017-07-21 01:45:03冀春俊徐浩孫琦孫淼
    風(fēng)機(jī)技術(shù) 2017年3期
    關(guān)鍵詞:輪盤(pán)端面葉輪

    冀春俊 徐浩 孫琦 孫淼

    (大連理工大學(xué)能源與動(dòng)力學(xué)院)

    基于Kriging模型的快速拆裝葉輪輪盤(pán)多目標(biāo)優(yōu)化設(shè)計(jì)

    冀春俊 徐浩 孫琦 孫淼

    (大連理工大學(xué)能源與動(dòng)力學(xué)院)

    為了降低葉輪安裝端面變形量以及減輕葉輪質(zhì)量,對(duì)快速拆裝形式的葉輪輪盤(pán)結(jié)構(gòu)進(jìn)行多目標(biāo)優(yōu)化設(shè)計(jì),使用拉丁超立方試驗(yàn)設(shè)計(jì)和有限元分析生成初始樣本數(shù)據(jù)并構(gòu)造出Kriging近似模型,運(yùn)用遺傳算法對(duì)近似模型進(jìn)行尋優(yōu)。對(duì)優(yōu)化后的參數(shù)值進(jìn)行工程取整并進(jìn)行有限元分析得出葉輪輪盤(pán)結(jié)構(gòu)設(shè)計(jì)模型,最終優(yōu)化結(jié)果顯示葉輪質(zhì)量減小了15%,葉輪摩擦端面軸向變形量降低了29%。優(yōu)化后的葉輪模型可提高拉桿螺栓連接轉(zhuǎn)子結(jié)構(gòu)的可靠性。

    快速拆裝葉輪;優(yōu)化設(shè)計(jì);Kriging模型

    0 引言

    離心壓縮機(jī)在石油、化工和造紙等國(guó)民生產(chǎn)領(lǐng)域中有著廣泛應(yīng)用。隨著離心壓縮機(jī)向著小流量、高速化方向發(fā)展,對(duì)葉輪提出越來(lái)越高的設(shè)計(jì)要求,葉輪結(jié)構(gòu)設(shè)計(jì)的好壞直接影響著整個(gè)機(jī)組的可靠性[1]。在現(xiàn)代高性能葉輪設(shè)計(jì)中,優(yōu)化設(shè)計(jì)是必不可少的工作,由于離心葉輪的葉片形狀設(shè)計(jì)是在氣動(dòng)設(shè)計(jì)中完成,因此對(duì)葉輪形狀設(shè)計(jì)主要針對(duì)輪盤(pán)。目前對(duì)離心葉輪進(jìn)行結(jié)構(gòu)優(yōu)化主要是降低應(yīng)力、減輕質(zhì)量和整體變形量,此外對(duì)不同安裝形式的葉輪其裝配工藝的優(yōu)化同樣重要[2],但是這些方面又存在相互制約的關(guān)系,因此有必要使用優(yōu)化設(shè)計(jì)的方法來(lái)平衡設(shè)計(jì)之間的沖突。

    Kriging模型的數(shù)學(xué)思想是在20世紀(jì)初隨著經(jīng)典統(tǒng)計(jì)學(xué)理論逐漸形成的,第一次成功運(yùn)用是地質(zhì)學(xué)家DanieKrige進(jìn)行地質(zhì)情況模擬及礦產(chǎn)分析。Kriging代理模型由線性回歸部分和非參數(shù)部分組成,它較單個(gè)的參數(shù)化模型更具有靈活性,同時(shí)克服了非參數(shù)化模型處理高維數(shù)據(jù)而存在的局限性[3]。目前,Kriging模型作為一種有效的近似技術(shù)已經(jīng)廣泛應(yīng)用在結(jié)構(gòu)優(yōu)化、多學(xué)科優(yōu)化設(shè)計(jì)、航空設(shè)計(jì)等工程優(yōu)化設(shè)計(jì)領(lǐng)域。

    本文運(yùn)用Kriging模型-遺傳算法對(duì)某壓縮機(jī)離心葉輪輪盤(pán)進(jìn)行形狀優(yōu)化設(shè)計(jì),選擇葉輪質(zhì)量和摩擦端面的軸向變形量最小為優(yōu)化目標(biāo)[4],約束條件為葉輪最大應(yīng)力和葉頂變形量。優(yōu)化后的葉輪在葉輪應(yīng)力和葉頂變形量符合要求的情況下質(zhì)量得到減輕,摩擦端面變形量大幅減少,對(duì)整個(gè)轉(zhuǎn)子運(yùn)轉(zhuǎn)的可靠性有著積極的影響。

    1葉輪結(jié)構(gòu)描述

    快速拆裝葉輪是利用一個(gè)兩端帶螺紋的拉桿來(lái)連接葉輪與軸,連接結(jié)構(gòu)如圖1所示。在安裝過(guò)程中,先用液壓工具將拉桿產(chǎn)生彈性變形拉長(zhǎng),然后擰緊鎖緊螺母,在液壓卸除后,將葉輪緊緊壓在軸的端面上,從而產(chǎn)生非常大的摩擦力,這個(gè)摩擦力在整個(gè)轉(zhuǎn)速范圍內(nèi)都能保證葉輪與軸不會(huì)發(fā)生松動(dòng)[5]。當(dāng)葉輪需要拆卸時(shí),先用液壓工具將拉桿拉長(zhǎng),然后松開(kāi)鎖緊螺母,卸除液壓后葉輪可與軸分離,整個(gè)過(guò)程實(shí)現(xiàn)快速拆裝。

    圖1葉輪連接結(jié)構(gòu)圖Fig.1Impeller connecting structure

    本文研究的葉輪為某單級(jí)高速直驅(qū)壓縮機(jī)主要零件,葉輪通過(guò)拉桿螺栓直接與高速電機(jī)軸相連接,通過(guò)端面摩擦將電機(jī)扭矩傳遞給葉輪,葉輪由7A09鋁合金加工而成,工作轉(zhuǎn)速20 000r/min。在高速旋轉(zhuǎn)時(shí)葉輪受到離心力的作用會(huì)產(chǎn)生一定的變形量,葉輪整體變形趨勢(shì)為朝向進(jìn)口方向。對(duì)于快速拆裝形式的葉輪,摩擦端面的軸向變形會(huì)減弱螺栓預(yù)緊力造成接觸面的分離,因此對(duì)此類結(jié)構(gòu)形式的轉(zhuǎn)子,葉輪合理的結(jié)構(gòu)設(shè)計(jì)和拉桿螺栓預(yù)緊力的準(zhǔn)確計(jì)算是保證整體結(jié)構(gòu)安全可靠的關(guān)鍵。

    2優(yōu)化設(shè)計(jì)

    2.1Kriging近似模型

    Kriging模型是一種估計(jì)方差最小的無(wú)偏估計(jì)模型,能夠很好地逼近具有非線性特征的復(fù)雜問(wèn)題。在計(jì)算過(guò)程中,Kriging代理模型計(jì)算量小同時(shí)計(jì)算結(jié)果與高精度復(fù)雜模型接近,在優(yōu)化設(shè)計(jì)中使用代理模型可以提高優(yōu)化效率,降低計(jì)算成本[6]。

    對(duì)選定p個(gè)訓(xùn)練樣本空間

    訓(xùn)練樣本所對(duì)應(yīng)的p個(gè)實(shí)際響應(yīng)值向量

    那么用Kriging代理模型來(lái)描述第i個(gè)訓(xùn)練樣本的模擬響應(yīng)向量y?i(x)

    其中,F(xiàn)(β,x)為回歸模型部分;z(x)是模型非參數(shù)部分相關(guān)函數(shù);z(x)服從正態(tài)分布,其均值為0,方差σ2,訓(xùn)練樣本點(diǎn)之間協(xié)方差

    θ是高斯相關(guān)函數(shù)的關(guān)鍵參數(shù),通過(guò)優(yōu)化θ能夠自適應(yīng)調(diào)節(jié)設(shè)計(jì)點(diǎn)之間的相關(guān)性

    相關(guān)函數(shù)確定以后,就可以建立觀測(cè)點(diǎn)的模擬響應(yīng)值y?,關(guān)于待測(cè)點(diǎn)的表達(dá)式

    建立Kriging近似方程后,常用的近似模型檢驗(yàn)方法有復(fù)相關(guān)系數(shù)(R2)、修正的復(fù)相關(guān)系數(shù)(adjustedR2)、均方根誤差(RMSE)等。

    2.2設(shè)計(jì)變量的選擇

    離心葉輪在不改變?nèi)~片形狀只對(duì)輪盤(pán)進(jìn)行優(yōu)化時(shí)

    的輪盤(pán)截面形狀及設(shè)計(jì)參數(shù)如圖2所示。

    圖2 輪盤(pán)結(jié)構(gòu)優(yōu)化設(shè)計(jì)參數(shù)Fig.2Design parameters of wheel structure

    設(shè)計(jì)參數(shù)P1為輪盤(pán)靠近葉輪出口處厚度值,P1減小可降低葉輪慣性載荷,P1增大可改善葉根應(yīng)力值。

    設(shè)計(jì)參數(shù)P2為錐部斜線與垂直面夾角,P2對(duì)減輕葉輪質(zhì)量,降低葉輪離心載荷,保證葉輪出口處剛度起到較重要作用。

    設(shè)計(jì)參數(shù)P3為控制錐部位置的參數(shù)值,P3對(duì)葉輪應(yīng)力以及質(zhì)量有重要影響。

    設(shè)計(jì)參數(shù)P4為輪盤(pán)背面錐部倒圓角,P4對(duì)改善葉輪應(yīng)力分布,減輕質(zhì)量,改善輪心處應(yīng)力狀況有一定作用,但是如果此處選擇不當(dāng)易形成高應(yīng)力區(qū)。

    設(shè)計(jì)參數(shù)P5為葉片根部倒圓角值,P5對(duì)葉根應(yīng)力值有較大影響。

    各設(shè)計(jì)參數(shù)相應(yīng)的初值及取值范圍如表1所示。

    表1 設(shè)計(jì)參數(shù)變化表Tab.1Design parameter variation

    2.3 優(yōu)化過(guò)程

    在對(duì)葉輪進(jìn)行多目標(biāo)優(yōu)化時(shí)利用三維建模軟件SolidWorks建立起葉輪關(guān)于設(shè)計(jì)變量的參數(shù)化模型,通過(guò)ANSYS Workbench平臺(tái)與SolidWorks之間的連接接口實(shí)現(xiàn)數(shù)據(jù)傳遞。采用拉丁超立方進(jìn)行試驗(yàn)設(shè)計(jì)并對(duì)試驗(yàn)點(diǎn)進(jìn)行有限元分析從而得到一組樣本點(diǎn),利用樣本數(shù)據(jù)建立起Kriging近似響應(yīng)面模型。對(duì)模型精度進(jìn)行分析,當(dāng)構(gòu)建的近似模型滿足精度要求時(shí)采用遺傳算法進(jìn)行尋優(yōu)得到Pareto最優(yōu)解,對(duì)得到的最優(yōu)解進(jìn)行數(shù)值圓整處理,從而得到最終的葉輪優(yōu)化形狀??焖俨鹧b葉輪輪盤(pán)結(jié)構(gòu)多目標(biāo)優(yōu)化流程如圖3所示。

    圖3 葉輪輪盤(pán)多目標(biāo)優(yōu)化流程圖Fig.3Multi objective optimization process of impeller disc

    由于葉輪葉片形狀已經(jīng)在氣動(dòng)設(shè)計(jì)中完成,因此,本文優(yōu)化主要針對(duì)輪盤(pán)背部形狀進(jìn)行。對(duì)于拉桿裝配連接的葉輪在高速狀態(tài)下工作時(shí),會(huì)出現(xiàn)摩擦端面距離變短的情況,這將減弱拉桿預(yù)緊力,造成葉輪與軸脫離,計(jì)算中應(yīng)考慮葉輪摩擦端面的軸向變形量[7]。本文優(yōu)化設(shè)計(jì)目標(biāo)確定為葉輪整體質(zhì)量最輕及葉輪摩擦端面的軸向變形量最小。

    根據(jù)本文葉輪的安裝形式,進(jìn)行有限元分析時(shí)對(duì)葉輪與鎖緊螺母連接的部位進(jìn)行全約束,對(duì)葉輪摩擦端面處約束周向位移,對(duì)葉輪整體施加慣性載荷。由于葉輪形狀具有周期對(duì)稱性,本文葉輪葉片數(shù)13個(gè),因此選取葉輪1/13扇區(qū)作為計(jì)算對(duì)象,葉輪單扇區(qū)三維模型及邊界條件設(shè)置如圖4所示。

    圖4 葉輪三維計(jì)算模型及邊界條件設(shè)置Fig.43D calculation model and boundary conditions of impeller

    2.4 優(yōu)化過(guò)程數(shù)學(xué)模型

    優(yōu)化設(shè)計(jì)問(wèn)題的數(shù)學(xué)模型:

    式中,f和y為目標(biāo)函數(shù);X為設(shè)計(jì)變量;gi為約束條件。

    在本文中目標(biāo)函數(shù)為葉輪質(zhì)量和葉輪摩擦端面軸向變形量最小,約束函數(shù)為葉輪應(yīng)力不超過(guò)材料許用應(yīng)力值350MPa,葉頂軸向變形量不超過(guò)0.5mm。

    3 葉輪輪盤(pán)多目標(biāo)優(yōu)化設(shè)計(jì)及結(jié)果分析

    3.1 葉輪輪盤(pán)結(jié)構(gòu)多目標(biāo)優(yōu)化

    優(yōu)化過(guò)程中常用的試驗(yàn)設(shè)計(jì)方法有中心復(fù)合設(shè)計(jì)、均與設(shè)計(jì)、正交設(shè)計(jì)和拉丁超立方設(shè)計(jì)等。本文采用拉丁超立方試驗(yàn)設(shè)計(jì)的方法對(duì)葉輪結(jié)構(gòu)5個(gè)設(shè)計(jì)變量進(jìn)行20次試驗(yàn)設(shè)計(jì),在有限元分析模塊中對(duì)這些試驗(yàn)樣本進(jìn)行計(jì)算,得到初始樣本數(shù)據(jù)點(diǎn)。選擇Kriging近似模型方法擬合樣本數(shù)據(jù),建立起近似響應(yīng)面模型。

    對(duì)于多目標(biāo)優(yōu)化問(wèn)題在優(yōu)化過(guò)程中會(huì)根據(jù)優(yōu)化問(wèn)題的目標(biāo)函數(shù)構(gòu)成一個(gè)適值函數(shù),對(duì)一個(gè)或多個(gè)解構(gòu)成種群進(jìn)行評(píng)估分析,因此多目標(biāo)優(yōu)化得出的結(jié)果是多個(gè)最優(yōu)解組成的集合,一般稱為Pareto最優(yōu)解。在Pareto最優(yōu)解集中任意兩個(gè)解不能直接對(duì)比其優(yōu)劣性,如何再進(jìn)行選擇要看實(shí)際問(wèn)題的需要[8]。

    本文采用多目標(biāo)遺傳算法MOGA進(jìn)行尋優(yōu),MOGA是基于NSGA-II(非支配排序遺傳算法-Ⅱ)的一個(gè)變種。它支持多種目標(biāo)和約束找到全局最優(yōu)。一個(gè)迭代的多目標(biāo)遺傳算法,提供比篩選一個(gè)更精確的方法,非常適合用于計(jì)算全局最大最小值,同時(shí)可以規(guī)避局部最優(yōu)的陷阱[9]。優(yōu)化設(shè)計(jì)時(shí)設(shè)定初始種群數(shù)為1 000,最大迭代數(shù)為20,優(yōu)化計(jì)算獲取3個(gè)最優(yōu)方案。

    3.2 優(yōu)化結(jié)果及分析

    利用上述優(yōu)化方法進(jìn)行優(yōu)化設(shè)計(jì)后得到3組優(yōu)化方案,各方案參數(shù)結(jié)果及葉輪力學(xué)性能見(jiàn)表2和表格3。

    表3 葉輪輪盤(pán)結(jié)構(gòu)多目標(biāo)優(yōu)化葉輪性能結(jié)果Tab.3Performance of multi objective optimization impeller disc structure

    從表2和表3中可以看出3個(gè)優(yōu)化方案數(shù)值比較接近,但優(yōu)化后的參數(shù)值不能為工程所應(yīng)用,從生產(chǎn)加工角度對(duì)優(yōu)化方案中各參數(shù)進(jìn)行數(shù)值圓整處理,并進(jìn)行有限元分析從而得到最終優(yōu)化方案的結(jié)果。對(duì)比最終優(yōu)化方案與葉輪原始方案的設(shè)計(jì)參數(shù)和有限元計(jì)算結(jié)果如表4和表5所示。從表5可以看出最終方案與優(yōu)化得到的3個(gè)方案結(jié)果相差不大,這種再次進(jìn)行有限元分析的方法可以檢驗(yàn)出本文構(gòu)建的Kriging近似模型精度較好,滿足真實(shí)情況。

    表4 優(yōu)化前后葉輪結(jié)構(gòu)尺寸對(duì)比Tab.4Comparison of impeller structure size before and after optimization

    表5 葉輪優(yōu)化前后結(jié)果對(duì)比Tab.5Comparison of results before and after optimization

    由表5可以看出通過(guò)多目標(biāo)優(yōu)化設(shè)計(jì)后的葉輪單個(gè)通道的質(zhì)量較原始葉輪單通道質(zhì)量減輕了15.71%,實(shí)現(xiàn)了葉輪減重的目標(biāo)。葉輪摩擦端面軸向變形量如圖5,最大軸向變形量為0.032 8mm比原始結(jié)構(gòu)減少29%,提高了工藝裝配的穩(wěn)定性,同時(shí)對(duì)拉桿預(yù)緊力設(shè)計(jì)提供了數(shù)值依據(jù)。

    圖5 摩擦端面軸向位移云圖Fig.5Axial displacement contour of friction end-face

    葉輪應(yīng)力結(jié)果如圖6,最大等效應(yīng)力值為344MPa出現(xiàn)在葉片根部靠近出口位置,葉頂軸向變形量如圖7,葉片頂部最大軸向變形量為0.448 9mm。葉輪最大

    圖6 優(yōu)化后葉輪應(yīng)力云圖Fig.6Stress contour of optimized impeller

    4 結(jié)論

    本文對(duì)快速拆裝這一結(jié)構(gòu)形式的葉輪作為研究對(duì)象,提出葉輪質(zhì)量和葉輪摩擦端面變形量為目標(biāo)函數(shù),葉輪最大應(yīng)力和葉頂變形量為約束條件的多目標(biāo)優(yōu)化。優(yōu)化過(guò)程為采用拉丁超立方進(jìn)行試驗(yàn)設(shè)計(jì)獲得樣本點(diǎn),通過(guò)樣本點(diǎn)構(gòu)造出Kriging近似響應(yīng)面模型,對(duì)Kriging模型運(yùn)用遺傳算法進(jìn)行優(yōu)化。

    通過(guò)Kriging模型-遺傳算法優(yōu)化的方法能有效減少設(shè)計(jì)樣本點(diǎn)的選取數(shù)量,提高優(yōu)化設(shè)計(jì)效率。對(duì)優(yōu)化后的數(shù)值進(jìn)行取整后得到最終優(yōu)化方案,與原始結(jié)構(gòu)相比優(yōu)化后葉輪質(zhì)量減小了15.71%,葉輪摩擦端面軸向變形量減小了29%。葉輪質(zhì)量減小后在高速旋轉(zhuǎn)狀態(tài)下離心力減小,結(jié)構(gòu)安全性得到提高。葉輪摩擦端面變形量減小對(duì)拉桿螺栓預(yù)緊力的影響變小,葉輪與軸接觸面連接更加可靠。應(yīng)力值和葉頂軸向變形量較葉輪原始結(jié)構(gòu)均有所增大但都在約束范圍內(nèi)。

    圖7 葉頂軸向位移云圖Fig.7Axial displacement contour of tip

    [1]James M.Sorokes,MarkJ.Kuzdzal,張海界.離心壓縮機(jī)的發(fā)展歷程[J].風(fēng)機(jī)技術(shù),2011(3):61-71.

    [2]王睿,王躍方.通風(fēng)機(jī)轉(zhuǎn)子過(guò)盈裝配松脫轉(zhuǎn)速數(shù)值分析[J].風(fēng)機(jī)技術(shù),2017(1):38-42.

    [3]賴喜濤,溫衛(wèi)東,馮大俊.基于Kriging模型的離心葉輪結(jié)構(gòu)優(yōu)化設(shè)計(jì)[J].沈陽(yáng)航空航天大學(xué)學(xué)報(bào),2014(4):17-22.

    [4]冀春俊,高亞威,孫琦,等.基于ANSYS的離心風(fēng)機(jī)葉輪結(jié)構(gòu)優(yōu)化研究[J].風(fēng)機(jī)技術(shù),2016(6):49-52.

    [5]冀春俊,仵忠浩,張晨陽(yáng),等.一種離心式風(fēng)機(jī)及用于連接其葉輪和轉(zhuǎn)軸的連接裝置[P].黑龍江:CN104929975A,2015.

    [6]李飛.代理模型預(yù)測(cè)研究及其在葉輪機(jī)械中的應(yīng)用[D].大連理工大學(xué),2014.

    [7]關(guān)振群,羅陽(yáng)軍,王學(xué)軍,等.離心式葉輪三維參數(shù)化形狀優(yōu)化設(shè)計(jì)方法[J].機(jī)械強(qiáng)度,2006(6):833-838.

    [8]唐明裴,閻貴平.結(jié)構(gòu)靈敏度分析及計(jì)算方法概述[J].中國(guó)鐵道科學(xué),2003(1):76-81.

    [9]陳國(guó)棟.基于代理模型的多目標(biāo)優(yōu)化方法及其在車身設(shè)計(jì)中的應(yīng)用[D].湖南大學(xué),2012.

    Multi-objective Structure Optimization of a Fast Disassembling Connection Impeller Disk Based on Kriging Model

    Chun-jun JiHao XuQi SunMiao Sun
    (School of energy and Power Engineering,Dalian University of Technology)

    In order to improve the installation reliability and reduce the weight of impeller,a multi-objective optimization method was used for the fast disassembling connection impeller disk structure.The initial sampling data was obtained and a Kriging model was constructed by a latin hypercube sampling design and a finite element method.A genetic algorithm was applied for global optimization.The optimized design parameters were obtained by a finite element analysis for the impeller wheel structure.After the optimization,the weight of impeller was decreased by 15%,the axial deformation of the impeller friction section was decreased by 49%.The optimized impeller model can improve the reliability of the rod bolt which connects the rotor structure.

    fast disassembling connection impeller,optimal design,Kriging model

    TH452;TK05

    1006-8155-(2017)03-0044-05

    A

    10.16492/j.fjjs.2017.03.0008

    2016-12-20遼寧大連116024

    猜你喜歡
    輪盤(pán)端面葉輪
    KDF3E成型機(jī)濾棒端面觸頭的原因及排除方法
    1.4317 QT2鋼在高能泵葉輪上的應(yīng)用
    某型航空發(fā)動(dòng)機(jī)鈦合金輪盤(pán)模擬疲勞試驗(yàn)件設(shè)計(jì)
    銅基合金襯套端面鍍鉻質(zhì)量的改善
    應(yīng)用石膏型快速精密鑄造技術(shù)制造葉輪
    基于ANSYS的輪盤(pán)轉(zhuǎn)子模態(tài)影響因素分析
    離心泵葉輪切割方法
    老舊端面磨齒機(jī)故障處理
    基于CFD/CSD耦合的葉輪機(jī)葉片失速顫振計(jì)算
    貫通軸端面花鍵拉床的設(shè)計(jì)
    久久中文字幕人妻熟女| 99国产精品99久久久久| 国产午夜福利久久久久久| 久久久国产欧美日韩av| 俄罗斯特黄特色一大片| xxxwww97欧美| 国产激情欧美一区二区| 亚洲人成网站在线播放欧美日韩| 男人舔女人下体高潮全视频| 欧美成人性av电影在线观看| 少妇的丰满在线观看| 亚洲中文日韩欧美视频| 自线自在国产av| 欧美午夜高清在线| 大香蕉久久成人网| 中文字幕最新亚洲高清| 一边摸一边做爽爽视频免费| 久久天堂一区二区三区四区| 日韩av在线大香蕉| 午夜久久久在线观看| 一本大道久久a久久精品| 嫩草影视91久久| 免费高清视频大片| 在线观看日韩欧美| 免费在线观看日本一区| 日韩中文字幕欧美一区二区| 国产成年人精品一区二区| 亚洲欧美一区二区三区黑人| cao死你这个sao货| 成人午夜高清在线视频 | 首页视频小说图片口味搜索| 国产成人一区二区三区免费视频网站| 精品第一国产精品| 无限看片的www在线观看| 精品人妻1区二区| 中文字幕精品免费在线观看视频| 久久伊人香网站| 后天国语完整版免费观看| 长腿黑丝高跟| 大型av网站在线播放| 欧美丝袜亚洲另类 | av福利片在线| 成人一区二区视频在线观看| 日本免费一区二区三区高清不卡| 国产精品自产拍在线观看55亚洲| 久久草成人影院| 中文资源天堂在线| 国产成人系列免费观看| 99国产精品一区二区蜜桃av| 高清在线国产一区| 亚洲中文字幕日韩| 在线天堂中文资源库| 午夜精品在线福利| 看黄色毛片网站| 成人国产一区最新在线观看| 欧美激情极品国产一区二区三区| xxxwww97欧美| 手机成人av网站| 日韩成人在线观看一区二区三区| 日韩av在线大香蕉| 两个人看的免费小视频| 12—13女人毛片做爰片一| 中文字幕人妻丝袜一区二区| 男人操女人黄网站| 人人妻人人看人人澡| 日韩欧美国产一区二区入口| 狠狠狠狠99中文字幕| 狂野欧美激情性xxxx| 国产亚洲精品综合一区在线观看 | 夜夜夜夜夜久久久久| 男女视频在线观看网站免费 | 69av精品久久久久久| e午夜精品久久久久久久| 真人做人爱边吃奶动态| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品久久久久久人妻精品电影| 91麻豆av在线| 亚洲熟女毛片儿| 国产精品免费一区二区三区在线| 久久久久九九精品影院| 中文在线观看免费www的网站 | 香蕉av资源在线| 免费搜索国产男女视频| 国产在线精品亚洲第一网站| 久久欧美精品欧美久久欧美| 天天添夜夜摸| www国产在线视频色| 搡老熟女国产l中国老女人| 精品人妻1区二区| 国内毛片毛片毛片毛片毛片| 久久伊人香网站| 亚洲人成网站在线播放欧美日韩| 亚洲精品美女久久久久99蜜臀| 制服诱惑二区| 亚洲av片天天在线观看| 亚洲国产看品久久| 亚洲五月婷婷丁香| 91在线观看av| avwww免费| a级毛片在线看网站| 国产片内射在线| 亚洲国产高清在线一区二区三 | √禁漫天堂资源中文www| 美女国产高潮福利片在线看| 久久精品国产清高在天天线| 久久久久国内视频| av在线天堂中文字幕| 制服诱惑二区| 麻豆一二三区av精品| 此物有八面人人有两片| 国产高清有码在线观看视频 | 高潮久久久久久久久久久不卡| 国产一级毛片七仙女欲春2 | 一级黄色大片毛片| 久久久久久国产a免费观看| bbb黄色大片| 国产亚洲精品综合一区在线观看 | 国产伦人伦偷精品视频| 久久久久精品国产欧美久久久| 国产亚洲精品久久久久久毛片| 在线免费观看的www视频| 久9热在线精品视频| 久久国产亚洲av麻豆专区| 婷婷精品国产亚洲av在线| 国内毛片毛片毛片毛片毛片| 天堂影院成人在线观看| 丝袜在线中文字幕| 性色av乱码一区二区三区2| 免费看日本二区| 免费在线观看视频国产中文字幕亚洲| 久久狼人影院| 桃红色精品国产亚洲av| 亚洲国产欧美网| 天天躁狠狠躁夜夜躁狠狠躁| 久久九九热精品免费| 国产在线精品亚洲第一网站| 2021天堂中文幕一二区在线观 | 精品国产超薄肉色丝袜足j| 美女高潮到喷水免费观看| 在线av久久热| 夜夜看夜夜爽夜夜摸| 长腿黑丝高跟| 成人亚洲精品av一区二区| 久热爱精品视频在线9| 久久99热这里只有精品18| 成人免费观看视频高清| 精品福利观看| 夜夜夜夜夜久久久久| 国产单亲对白刺激| 亚洲五月色婷婷综合| 亚洲久久久国产精品| 三级毛片av免费| 中文字幕精品免费在线观看视频| 天堂动漫精品| 久久久久国产一级毛片高清牌| 国产精品99久久99久久久不卡| 欧美黑人欧美精品刺激| 亚洲成国产人片在线观看| 变态另类成人亚洲欧美熟女| 国产亚洲欧美在线一区二区| 亚洲国产中文字幕在线视频| 在线观看www视频免费| 亚洲中文av在线| 欧美黄色淫秽网站| 亚洲精品在线美女| 男女之事视频高清在线观看| 国产激情久久老熟女| 在线观看www视频免费| 精品久久久久久久毛片微露脸| 久久草成人影院| 久久香蕉精品热| 黄色视频不卡| 国产单亲对白刺激| 国产爱豆传媒在线观看 | 欧美成狂野欧美在线观看| 精品国内亚洲2022精品成人| 麻豆一二三区av精品| 97人妻精品一区二区三区麻豆 | 九色国产91popny在线| 香蕉丝袜av| 亚洲精品粉嫩美女一区| 一个人观看的视频www高清免费观看 | 香蕉国产在线看| 国产人伦9x9x在线观看| 757午夜福利合集在线观看| 久久欧美精品欧美久久欧美| 国产三级黄色录像| tocl精华| 制服人妻中文乱码| 少妇 在线观看| 欧美av亚洲av综合av国产av| 国产亚洲欧美在线一区二区| a级毛片a级免费在线| 91麻豆av在线| 韩国精品一区二区三区| 久久久久九九精品影院| 欧美国产精品va在线观看不卡| 香蕉av资源在线| 欧美日韩一级在线毛片| 国产精品亚洲一级av第二区| 狠狠狠狠99中文字幕| 午夜免费鲁丝| 日本 av在线| 久9热在线精品视频| 亚洲最大成人中文| 国产三级黄色录像| 男女那种视频在线观看| 丁香欧美五月| 黄色片一级片一级黄色片| 久久久久国产精品人妻aⅴ院| 天堂√8在线中文| 精品午夜福利视频在线观看一区| 麻豆久久精品国产亚洲av| 性色av乱码一区二区三区2| 人人澡人人妻人| 欧美av亚洲av综合av国产av| bbb黄色大片| 欧美精品亚洲一区二区| 亚洲精品国产精品久久久不卡| 久久久久国产一级毛片高清牌| 丝袜在线中文字幕| 午夜福利免费观看在线| 亚洲成a人片在线一区二区| 一本精品99久久精品77| 亚洲男人的天堂狠狠| 亚洲欧美日韩高清在线视频| 国产精品 欧美亚洲| 亚洲av成人一区二区三| 欧美午夜高清在线| 特大巨黑吊av在线直播 | 日本撒尿小便嘘嘘汇集6| 黑人巨大精品欧美一区二区mp4| 亚洲无线在线观看| 国产精品一区二区免费欧美| 男女下面进入的视频免费午夜 | 在线观看免费日韩欧美大片| 亚洲 欧美一区二区三区| 日韩精品免费视频一区二区三区| 国产日本99.免费观看| 午夜两性在线视频| 亚洲一码二码三码区别大吗| 一级毛片高清免费大全| 一区二区日韩欧美中文字幕| 欧美黑人巨大hd| 啦啦啦观看免费观看视频高清| 久久婷婷成人综合色麻豆| 免费在线观看完整版高清| 可以免费在线观看a视频的电影网站| 欧美三级亚洲精品| 国产精品久久久久久人妻精品电影| 麻豆一二三区av精品| 国产免费av片在线观看野外av| av片东京热男人的天堂| ponron亚洲| 免费看十八禁软件| 一本大道久久a久久精品| 久久久久久久久久黄片| 少妇的丰满在线观看| 12—13女人毛片做爰片一| 最近最新免费中文字幕在线| 在线十欧美十亚洲十日本专区| 国产精品九九99| 日本在线视频免费播放| 日韩 欧美 亚洲 中文字幕| 欧美成狂野欧美在线观看| 国产亚洲精品综合一区在线观看 | 两个人看的免费小视频| 老司机靠b影院| 琪琪午夜伦伦电影理论片6080| 美女 人体艺术 gogo| 久久精品亚洲精品国产色婷小说| 男人舔女人下体高潮全视频| 国产伦人伦偷精品视频| 亚洲avbb在线观看| 欧美日本视频| 一边摸一边做爽爽视频免费| 国产一区二区在线av高清观看| 亚洲avbb在线观看| 啦啦啦免费观看视频1| 国产成人一区二区三区免费视频网站| 12—13女人毛片做爰片一| 麻豆av在线久日| 99国产精品一区二区三区| 啦啦啦 在线观看视频| 在线播放国产精品三级| 一区二区日韩欧美中文字幕| 午夜福利成人在线免费观看| 亚洲成人精品中文字幕电影| 中文字幕精品亚洲无线码一区 | 99国产精品一区二区蜜桃av| 一本综合久久免费| 亚洲av电影在线进入| 国产精品一区二区免费欧美| 国产一卡二卡三卡精品| 国产久久久一区二区三区| 2021天堂中文幕一二区在线观 | 国产在线观看jvid| 欧美丝袜亚洲另类 | 麻豆久久精品国产亚洲av| 国产精品影院久久| 性色av乱码一区二区三区2| 国产成人欧美在线观看| 国产亚洲精品一区二区www| 午夜日韩欧美国产| 在线观看www视频免费| 成人午夜高清在线视频 | 日韩av在线大香蕉| 亚洲自偷自拍图片 自拍| 亚洲中文av在线| 50天的宝宝边吃奶边哭怎么回事| 国产91精品成人一区二区三区| 国产亚洲欧美精品永久| 美女免费视频网站| 成人国语在线视频| 国产三级在线视频| 国产97色在线日韩免费| 亚洲av成人不卡在线观看播放网| 午夜老司机福利片| 日韩欧美国产一区二区入口| 非洲黑人性xxxx精品又粗又长| 欧美日韩黄片免| 国产精品久久视频播放| 在线观看66精品国产| 51午夜福利影视在线观看| 色综合婷婷激情| 日韩三级视频一区二区三区| 脱女人内裤的视频| 黄色视频不卡| 真人做人爱边吃奶动态| 久久久国产精品麻豆| 亚洲欧美日韩高清在线视频| www.999成人在线观看| 级片在线观看| 国产蜜桃级精品一区二区三区| 中出人妻视频一区二区| 国产一区在线观看成人免费| 国产v大片淫在线免费观看| 99国产极品粉嫩在线观看| 国产精品美女特级片免费视频播放器 | aaaaa片日本免费| av免费在线观看网站| 日韩精品青青久久久久久| 国产主播在线观看一区二区| 午夜福利免费观看在线| 别揉我奶头~嗯~啊~动态视频| 精品高清国产在线一区| 欧美三级亚洲精品| 国产蜜桃级精品一区二区三区| 操出白浆在线播放| 美女国产高潮福利片在线看| 天堂√8在线中文| 中出人妻视频一区二区| 欧美亚洲日本最大视频资源| 免费观看人在逋| 午夜久久久久精精品| 久久精品人妻少妇| 99国产精品一区二区蜜桃av| 国产激情欧美一区二区| 美女午夜性视频免费| 国产色视频综合| 免费在线观看视频国产中文字幕亚洲| 色综合亚洲欧美另类图片| 精品久久久久久成人av| 草草在线视频免费看| 欧美日韩亚洲国产一区二区在线观看| 男人舔女人下体高潮全视频| 在线av久久热| 午夜视频精品福利| 欧美另类亚洲清纯唯美| av欧美777| 大型av网站在线播放| 欧美大码av| 999久久久国产精品视频| 久久久久久久午夜电影| 欧美黄色淫秽网站| 国产视频一区二区在线看| 久久久久久大精品| 欧美一级毛片孕妇| 亚洲人成网站在线播放欧美日韩| 欧美性猛交╳xxx乱大交人| 色尼玛亚洲综合影院| 一区二区三区激情视频| 免费av毛片视频| 日韩大尺度精品在线看网址| 午夜免费成人在线视频| 午夜福利在线观看吧| 一区二区三区国产精品乱码| 午夜免费观看网址| 欧美丝袜亚洲另类 | 一区二区三区国产精品乱码| 国产在线观看jvid| 给我免费播放毛片高清在线观看| 搡老熟女国产l中国老女人| 18禁国产床啪视频网站| 操出白浆在线播放| 久久欧美精品欧美久久欧美| 亚洲成国产人片在线观看| 国产精品久久久av美女十八| 国产成人啪精品午夜网站| 国产av一区二区精品久久| 天堂影院成人在线观看| 十八禁网站免费在线| 高潮久久久久久久久久久不卡| 99在线人妻在线中文字幕| 国产激情偷乱视频一区二区| 精品久久久久久久毛片微露脸| 欧美日本视频| 欧美性猛交╳xxx乱大交人| 又黄又粗又硬又大视频| 日日干狠狠操夜夜爽| www日本在线高清视频| 啪啪无遮挡十八禁网站| 成人特级黄色片久久久久久久| tocl精华| 神马国产精品三级电影在线观看 | 男男h啪啪无遮挡| 免费高清在线观看日韩| 欧美精品啪啪一区二区三区| 久久精品91无色码中文字幕| 亚洲中文字幕一区二区三区有码在线看 | 黄色女人牲交| 50天的宝宝边吃奶边哭怎么回事| 视频区欧美日本亚洲| 久久草成人影院| 在线视频色国产色| 色尼玛亚洲综合影院| 亚洲,欧美精品.| 91老司机精品| 精品久久久久久久毛片微露脸| 国产伦一二天堂av在线观看| 一个人观看的视频www高清免费观看 | 两性夫妻黄色片| 男男h啪啪无遮挡| 亚洲人成电影免费在线| 少妇熟女aⅴ在线视频| 国产一区二区激情短视频| 亚洲av电影不卡..在线观看| 久久久久久九九精品二区国产 | 久久人妻福利社区极品人妻图片| 日本在线视频免费播放| 午夜福利18| 麻豆成人午夜福利视频| 国产黄片美女视频| 国产一区二区三区在线臀色熟女| 伦理电影免费视频| 亚洲美女黄片视频| 亚洲色图 男人天堂 中文字幕| 国产乱人伦免费视频| 精品久久久久久久毛片微露脸| 高潮久久久久久久久久久不卡| 草草在线视频免费看| 欧美黄色淫秽网站| 999久久久国产精品视频| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲 欧美 日韩 在线 免费| 搡老妇女老女人老熟妇| 男男h啪啪无遮挡| 色综合亚洲欧美另类图片| 黑丝袜美女国产一区| 老鸭窝网址在线观看| 女性被躁到高潮视频| 欧美黑人欧美精品刺激| 国产亚洲欧美98| 此物有八面人人有两片| 欧美日韩一级在线毛片| 麻豆成人av在线观看| 日韩大码丰满熟妇| 久热爱精品视频在线9| 狠狠狠狠99中文字幕| 搡老熟女国产l中国老女人| 亚洲人成77777在线视频| 精品国产国语对白av| 国产一区二区激情短视频| 一级毛片女人18水好多| 久久久久久久久久黄片| 中文字幕最新亚洲高清| 一进一出好大好爽视频| 国产精品永久免费网站| 日韩精品免费视频一区二区三区| 精品欧美一区二区三区在线| 午夜福利18| 欧美中文日本在线观看视频| 亚洲自拍偷在线| 97超级碰碰碰精品色视频在线观看| 波多野结衣高清无吗| 在线观看免费午夜福利视频| 国产爱豆传媒在线观看 | 久久久久久久精品吃奶| 国产免费男女视频| 99国产综合亚洲精品| 国产伦在线观看视频一区| 免费看美女性在线毛片视频| 91在线观看av| 听说在线观看完整版免费高清| 韩国av一区二区三区四区| 天天添夜夜摸| 国产亚洲精品一区二区www| 欧美激情高清一区二区三区| 夜夜躁狠狠躁天天躁| 亚洲成人国产一区在线观看| 亚洲欧洲精品一区二区精品久久久| avwww免费| 丰满人妻熟妇乱又伦精品不卡| 欧美黑人精品巨大| 午夜成年电影在线免费观看| 91成人精品电影| 成人国语在线视频| 美女免费视频网站| 午夜福利成人在线免费观看| 国产不卡一卡二| 老司机午夜福利在线观看视频| 美女免费视频网站| 国产亚洲欧美精品永久| 日韩精品免费视频一区二区三区| 午夜久久久在线观看| 日韩欧美一区视频在线观看| 日本三级黄在线观看| 国产私拍福利视频在线观看| 欧美另类亚洲清纯唯美| avwww免费| 伦理电影免费视频| 国产成年人精品一区二区| 国内精品久久久久精免费| 19禁男女啪啪无遮挡网站| 一级毛片精品| 成在线人永久免费视频| 久久久久免费精品人妻一区二区 | 精品午夜福利视频在线观看一区| 啦啦啦免费观看视频1| 精品欧美一区二区三区在线| 欧美在线黄色| 日本在线视频免费播放| 久久精品国产亚洲av香蕉五月| 在线播放国产精品三级| 亚洲,欧美精品.| 久久青草综合色| 天堂动漫精品| 亚洲男人的天堂狠狠| 国产精品久久久久久亚洲av鲁大| 午夜激情福利司机影院| 亚洲三区欧美一区| 99国产综合亚洲精品| 一本综合久久免费| 69av精品久久久久久| 老熟妇乱子伦视频在线观看| 国产精品久久久久久人妻精品电影| 老司机在亚洲福利影院| 午夜福利在线观看吧| 欧美三级亚洲精品| 中文字幕人妻熟女乱码| 久久欧美精品欧美久久欧美| 亚洲成国产人片在线观看| av在线播放免费不卡| svipshipincom国产片| 久久中文看片网| 免费观看精品视频网站| 成人免费观看视频高清| 91在线观看av| 男女之事视频高清在线观看| 不卡av一区二区三区| 免费人成视频x8x8入口观看| 一边摸一边做爽爽视频免费| 又大又爽又粗| 亚洲精品中文字幕在线视频| 国产成人av教育| 男女午夜视频在线观看| 婷婷丁香在线五月| 一区二区日韩欧美中文字幕| 国产aⅴ精品一区二区三区波| 日本精品一区二区三区蜜桃| 日日爽夜夜爽网站| 免费高清在线观看日韩| 欧美绝顶高潮抽搐喷水| 亚洲黑人精品在线| 亚洲三区欧美一区| 麻豆国产av国片精品| 自线自在国产av| 国产区一区二久久| 亚洲第一av免费看| 午夜福利成人在线免费观看| 日韩欧美一区二区三区在线观看| 国产aⅴ精品一区二区三区波| 在线天堂中文资源库| 亚洲av片天天在线观看| 男女之事视频高清在线观看| 国产精品99久久99久久久不卡| 黄片大片在线免费观看| 999久久久国产精品视频| 首页视频小说图片口味搜索| 一区二区日韩欧美中文字幕| 91老司机精品| 亚洲精品国产精品久久久不卡| 久久精品亚洲精品国产色婷小说| 一个人观看的视频www高清免费观看 | 午夜免费成人在线视频| 少妇的丰满在线观看| 色综合亚洲欧美另类图片| 女人被狂操c到高潮| 少妇被粗大的猛进出69影院| 夜夜躁狠狠躁天天躁| 午夜福利免费观看在线| 99国产精品一区二区三区| 在线观看午夜福利视频| 不卡一级毛片| x7x7x7水蜜桃| 精华霜和精华液先用哪个| 亚洲第一电影网av| 深夜精品福利| 真人做人爱边吃奶动态| 亚洲av第一区精品v没综合| 侵犯人妻中文字幕一二三四区| 后天国语完整版免费观看| av中文乱码字幕在线| 99热这里只有精品一区 | 美女 人体艺术 gogo|