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

    基于Kriging模型與NSGA-II算法的堆垛機結(jié)構(gòu)優(yōu)化設計

    2022-05-19 03:34:56楊瑞剛徐格寧
    機械設計與制造 2022年5期
    關鍵詞:堆垛雙高橫梁

    楊瑞剛,高 林,董 青,徐格寧

    (太原科技大學機械工程學院,山西 太原 030024)

    1 引言

    隨著物流倉儲的快速發(fā)展,巷道式堆垛機應時代而生,巷道式堆垛機作為物流倉儲的專用堆垛機,是物流倉儲中最要的搬運設備。隨著用地面積的緊張,堆垛機正向高速度和大高度方向發(fā)展,這對堆垛機的設計提出更高的要求,針對國內(nèi)的堆垛機性能指標落后,結(jié)構(gòu)不合理、笨重等特點[1],提出了對現(xiàn)有堆垛機結(jié)構(gòu)在滿足結(jié)構(gòu)強度、剛度、穩(wěn)定性條件下的快速優(yōu)化設計。

    雙高雙立柱式堆垛機在啟動,加速運行,停止時,堆垛機金屬結(jié)構(gòu)受到各種載荷作用,易發(fā)生撓曲變形,產(chǎn)生局部應力集中[2],故對堆垛機結(jié)構(gòu)優(yōu)化就顯得異常重要。文獻[2]通過運用參數(shù)化建模運用ANSYS求解,以立柱頂部撓度和立柱根部應力為優(yōu)化目標,通過改變單個設計參數(shù)尋優(yōu),具有一定局限性,容易陷入局部最優(yōu)。文獻[1]通過運用MOGA遺傳算法對單立柱堆垛機設置不同的約束條件進行靜態(tài)優(yōu)化。但由于設計參數(shù)較少,當設計參數(shù)關系復雜時,對其進行直接優(yōu)化,工作量巨大,效率極低。而代理模型可以減少調(diào)用真實模型次數(shù),提高優(yōu)化效率,文獻[3]運用Kriging模型和穩(wěn)健性設計相結(jié)合對乘務員約束系統(tǒng)進行穩(wěn)健性設計,提高了設計質(zhì)量。文獻[4]通過kriging模型和遺傳算法對泵葉輪進行了尺寸優(yōu)化,優(yōu)化后效率提高5.53%。文獻[5]通過Krig?ing模型和模擬退火粒子群算法對桁架模型進行了模型修正和損傷識別等驗證。文獻[6]將Kriging模型與傳統(tǒng)優(yōu)化算法和遺傳算法進行對比,確定Kriging模型的有效性。并對翼型進行優(yōu)化,翼型升阻比提高了81%,優(yōu)化效果明顯。文獻[7]通過將代理模型應用到飛機旋翼上進行優(yōu)化設計。文獻[8]通過運用Kriging模型可以對非線性模型有很好的逼近能力。Kriging廣泛應用于工程結(jié)構(gòu)和飛行器多學科優(yōu)化中,但將Kriging模型用于堆垛機整機優(yōu)化還很少。

    鑒于此,這里提出一種適合雙高(速度高與高度高)堆垛機結(jié)構(gòu)特點的方案,即基于Kriging 模型與非支配排序遺傳算法(NSGA-II)的相結(jié)合的優(yōu)化方法,將雙高雙立柱堆垛機的位移和質(zhì)量作為待優(yōu)化目標,綜合考慮雙高雙立柱整機的約束情況,采用參數(shù)化三維建模,通過最佳填充樣本空間(OSF)進行試驗設計,考慮不同設計參數(shù)進行組合,避免陷入單個設計參數(shù)局部最優(yōu)。通過試驗設計構(gòu)建克里金代理模型代替有限元模型(真實模型),可以最大程度避免在進行優(yōu)化時調(diào)用真實模型次數(shù),提高了優(yōu)化效率。最后通過精英保留策略快速非支配排序多目標遺傳(NSGA-II)算法對Kriging代理模型尋找全局最優(yōu)Pareto解。

    2 雙高雙立柱式堆垛機優(yōu)化設計流程

    由于雙高雙立柱式堆垛機結(jié)構(gòu)是一個非線性問題,對其進行直接優(yōu)化設計,將不斷調(diào)用真實模型進行尋優(yōu),將大大降低了優(yōu)化效率,為了解決此類優(yōu)化設計問題,通過構(gòu)造代理模型將成為關鍵。針對雙高雙立柱式堆垛機的特性,故需要找到一個精度高和能描述高度非線性過程的代理模型,文獻[3]中指出通過構(gòu)造Kriging 構(gòu)造代理模型不僅可以描述高度非線性過程,還可以去除噪聲。被廣泛應用于高度非線性問題。

    基于Kriging模型的雙高雙立柱式堆垛機靜態(tài)多目標優(yōu)化設計流程,如圖1所示。

    圖1 堆垛機結(jié)構(gòu)優(yōu)化整體流程圖Fig.1 Overall Flowchart of Stacker Structure Optimization

    具體過程如下:

    (1)對雙高雙立柱式堆垛機進行理論分析并通過構(gòu)建有限元模型進行驗證,確保優(yōu)化設計模型的正確性;

    (2)基于最佳樣本空間(OSF)的試驗設計方法和Kriging理論建立代理模型;

    (3)通過NSGA-II 多目標優(yōu)化算法對Kriging 代理模型進行求解,獲得多個Pareto非劣解集;

    (4)通過評價準則確定最佳Pareto非劣解。

    3 雙高雙立柱堆垛機結(jié)構(gòu)靜力學分析

    這里研究對象為雙高雙立柱式堆垛機,其基本參數(shù),如表1所示。雙高雙立柱式堆垛機模型圖,如圖2所示。

    表1 堆垛機基本參數(shù)Tab.1 Stacker Basic Parameters

    圖2 堆垛機模型Fig.2 Stacker Model

    其中為主立柱、輔立柱、上橫梁、下橫梁各截面形式,如圖3所示。其中L1為主立柱下端腹板長度,T1為主立柱下端翼板厚度,其中L2為主立柱上端腹板長度,T2為主立柱上端翼板厚度,T3為輔立柱腹板厚度,T4為輔立柱翼板厚度,T5為上橫梁腹板厚度,T4為上橫梁翼板厚度,T7為下橫梁腹板厚度,T8為下橫梁翼板厚度。

    圖3 各截面形式Fig.3 Each Section Form

    3.1 堆垛機力學模型與有限元模型

    當雙高雙立柱式堆垛機載貨臺運行到最高位置時,貨叉全伸叉取滿載貨物。為最危險工況[9],其堆垛機門架受力示意圖,如圖4所示。堆垛機門架結(jié)構(gòu)力學模型,如圖5所示。

    圖4 堆垛機門架受力示意圖Fig.4 Stacking Machine Frame Force Diagram

    圖5 堆垛機門架結(jié)構(gòu)力學模型Fig.5 Mechanical Model of Door Frame Structure of Stacker

    其中,自重載荷為:G1為上橫梁及附件重力,G2為載貨臺、貨物、司機室及附件的重力,G3為起升裝置的重力,G4為電氣控制柜的重力,G5為下橫梁及附件重力,G6、G7為立柱重力,G8為運行機構(gòu)自重,GZ1和GZ2分別為主立柱和輔立柱單位長度平均重力,qs和qx分別為上、下橫梁單位長度平均重力,而h1~h5和h8分別對應于G1~G5和G8分別到下橫梁中心線的質(zhì)心高度。堆垛機在加速或者減速行走時產(chǎn)生的水平力為:H1~H5和H8對應于G1~G5和G8的水平慣性力,q1、q2分別為立柱的單位長度平均質(zhì)量產(chǎn)生的水平慣性力。F1和F2分別為主立柱和輔立柱頂部作用力。

    根據(jù)雙高雙立柱式堆垛機的工作特點、不同工況、考慮各項載荷實際出現(xiàn)的機率,按對結(jié)構(gòu)最不利的作用情況,將可能同時出現(xiàn)的載荷分別按A1和A4進行組合,通過對比得到合理的載荷組合[10]。

    按A1進行載荷組合時,取?1=1.1,起升動載系數(shù)?2按下式計算:

    式中:?2min—與起升狀態(tài)級別相對應的起升動載系數(shù)的最小值;

    β2—按起升狀態(tài)級別設定的系數(shù);

    Vq—穩(wěn)定(額定)起升速度,單位為m/s。

    載荷組合A1下的總輪壓:

    按A4進行載荷組合時運行沖擊系數(shù)按式計算:

    載荷組合A4下的總輪壓:

    由于LUy4>LUy1因此因按載荷組合A4進行計算。

    3.2 雙高雙立柱堆垛機有限元模型

    雙高雙立柱式堆垛機的金屬結(jié)構(gòu)主要有,主立柱(變截面),輔立柱,上橫梁,下橫梁組成。建立雙高雙立柱式堆垛機的有限元的模型時要簡化實際模型的中的一些不重要部件,使建立的有限元模型既可以反應實際模型的真實狀況,又可以得到較為精確的結(jié)果。

    雙高雙立柱式堆垛機的模型采用Solid works參數(shù)化建模,導入ANSYS workbench中,其中單元采用solid45單元,因為實體單元可以更好的模擬實際情況。主立柱,輔立柱,上橫梁,下橫梁之間均采用ANSYS workbench中的bonded接觸。采用自由網(wǎng)格劃分,建立有限元模型。其單元總數(shù)為17597,節(jié)點總數(shù)為109516,如圖6所示。

    圖6 有限元模型Fig.6 Finite Element Model

    根據(jù)雙高雙立柱式堆垛機的實際工作情況,在下橫梁底部兩端面與上橫梁兩側(cè)面施加面約束,約束類型為:在A面施加固定約束,約束X,Y,Z三個方向自由度,在B面施加位移約束,限制住Y方向和Z方向的自由度,釋放掉X方向的自由度。在C面和D面施加位移約束,限制Z方向的自由度,防止堆垛機前后擺動,使其始終保為平面,釋放掉X和Y方向自由度。將堆垛機力學模型上面的力等效施加到有限元模型上,如圖7所示。

    圖7 載荷施加Fig.7 Load Application

    3.3 有限元求解

    通過云圖9可知,在載貨臺移動到最高位置,且滿載貨物貨叉伸出時的工況下,最大等效應力位移輔立柱下端,應力值為103.21MPa,且整個堆垛機金屬結(jié)構(gòu)在x方向產(chǎn)生的位移為41.157mm,小于許用位移值。由于堆垛機主立柱材料采用Q345,許用應力為[σ]=0.5*σs+0.35*σb ni=227MPa ,輔立柱材料采用Q235,許用應力為158MPa,由云圖9可知,應力均小于許用應力。

    圖8 堆垛機位移云圖Fig.8 Stacker Displacement Cloud Diagram

    圖9 堆垛機應力云圖Fig.9 Stacker Stress Cloud Diagram

    4 Kriging代理模型構(gòu)建

    由于雙高雙立柱式堆垛機在結(jié)構(gòu)分析時為非線性過程,而Kriging 模型可以很好的描述非線性過程,同時也能光滑目標響應、去除噪聲和提高優(yōu)化設計效率[3]。Kriging模型的主要是有線性回歸部分和隨機分布部分所組成[11],其樣本點X=[x1,x2,x3…xn]T和響應值Y=[y1,y2,y3…yn]T,其中,樣本點和響應點通過試驗設計法中的最佳空間填充(OSF)來獲得。其Krig?ing模型的數(shù)學表達式為:

    式中:Y(x)—未知Kriging 模型;fi(x)—樣本點的多項式回歸組合,提供了設計空間的全局近似模擬;βi—回歸系數(shù)的最小二乘估計;z(x)—從正態(tài)分布的隨機項,是在全局近似模擬的基礎上創(chuàng)建的局部偏差,其協(xié)方差方程表達式為:

    式中:R—N*N的協(xié)方差矩陣,沿對角線正定;r[xi,xj]—任意N個樣本點中xi和xj空間函數(shù)相關,為高斯的相關函數(shù):

    式中:θK—用于擬合模型的不確定參數(shù);M—設計參數(shù)的數(shù)量。

    通過確定相關函數(shù),則樣本點x處的響應估計為:

    式中:y—樣本點數(shù)據(jù)的響應值所組成的列向量;f—單位列矢量;

    rT(x)—樣本點和預測點之間得相關矢量,如式(9)所示:

    其中式(8)中β估計和式(6)全局模型方差估計值σ2如下式所示:

    通過最大似然估計求解相關參數(shù)θK得到:

    通過式(12)求解出θK后,可得到Kriging擬合精度最優(yōu)的代理模型。

    由于雙高雙立柱式堆垛機為非線性問題,因此在構(gòu)建Krig?ing代理模型時,試驗點的選取就顯得異常重要。這里通過采用最佳填充空間(OSF)設計得到試驗設計方法,最佳填充樣本空間(OSF)相對于中心復合設計(CCD)可在整個設計空間均勻的分配設計參數(shù),以最小的數(shù)量獲得對設計點的最大洞察。根據(jù)設計點各個參數(shù)的取值范圍,生成85組的樣本點,進行計算,通過最佳填充空間(OSF)得到的樣本點和響應值帶入Kriging模型中構(gòu)建響應面模型。通過擬合優(yōu)度散點圖可以觀測Kriging模型擬合出的響應面精度好壞,其中,樣本點越接近對角線,則表示Kriging模型擬合出的響應面精度越高,如圖10所示。

    圖10 擬合優(yōu)度散點圖Fig.10 Goodness of Fit Scatter Plot

    設計參數(shù)與響應值通過Kriging模型擬合出的響應面圖,如圖11~圖13所示??梢酝ㄟ^對不同的設計參數(shù)進行組合,得到設計參數(shù)與響應值不同的響應面圖。通過響應面圖還可以看出各個設計參數(shù)對響應值之間得影響。

    圖11 設計參數(shù)與應力的擬合響應面Fig.11 Fitting Response Surface of Design Parameters and Stress

    圖12 設計參數(shù)與位移的擬合響應面圖Fig.12 Fitting Response Surface of Design Parameters and Displacement

    圖13 設計參數(shù)與質(zhì)量的擬合響應面圖Fig.13 Fitting Response Surface of Design Parameters and Mass

    5 堆垛機結(jié)構(gòu)多目標優(yōu)化數(shù)學模型

    基于雙高雙立柱堆垛機金屬結(jié)構(gòu)模型,建立以主立柱下端腹板長度L1、主立柱下端翼板厚度T1、主立柱上端腹板長度L2、主立柱上端翼板厚度、T2輔立柱腹板厚度T3、輔立柱翼板厚度T4、上橫梁腹板厚度T5、上橫梁翼板厚度T4、下橫梁腹板厚度T7、下橫梁翼板厚度T8為設計參數(shù),堆垛機結(jié)構(gòu)強度為約束條件,堆垛機結(jié)構(gòu)剛度和質(zhì)量為優(yōu)化目標的多目標數(shù)學優(yōu)化模型如下所示:

    式中:m—堆垛機結(jié)構(gòu)質(zhì)量;y—堆垛機結(jié)構(gòu)最大位移;σmax—堆垛機結(jié)構(gòu)最大應力值;[σ]—鋼材許用值;xi—設計參數(shù)。

    5.1 優(yōu)化算法

    多目標遺傳算法(NSGA-II)是通過精英保留策略快速非支配排序的算法,是一種基本Pareto最優(yōu)解集的多目標算法[13-14]。其分析流程圖,如圖14所示。其具體步驟為:(1)在變量范圍內(nèi)隨機生成初始化種群,根據(jù)Pareto 的等級關系對生成的種群進行非支配排序后,帶入遺傳算法中進行選擇、交叉、變異后,生成第一代子代種群。(2)第二代開始,將父代與子代合并,并進行非支配排序,同時對非支配關系層中的個體進行擁擠度計算,使解在目標空間更加均勻。根據(jù)個體擁擠度和非支配關系選擇合適的種群組成新的父代種群,然后在帶入遺傳算法進行選擇、交叉、變異生成新的子代種群。循環(huán)進行,直到滿足進化代數(shù)要求[15]。

    圖14 NSGA-II算法流程圖Fig.14 NSGA-II Algorithm Flowchart

    5.2 模型優(yōu)化對比

    由于雙高雙立柱式堆垛機有限元模型單元數(shù)量眾多,為了提高雙高雙立柱式堆垛機的優(yōu)化效率。分別采用NSGA-II多目標優(yōu)化算法對堆垛機有限元模型進行優(yōu)化和基于最佳樣本空間(OSF)試驗設計構(gòu)建的Kriging代理模型進行優(yōu)化對比,由于直接優(yōu)化與初始種群規(guī)模有直接影響,為了減少仿真時間,NSGA-II參數(shù)設置:初始種群規(guī)模為100,最大允許的Pareto 百分比為70%,進化代數(shù)50,突變因子概率為0.02,交叉因子概率為0.98。其對比結(jié)果,如表2所示。

    表2 不同模型的優(yōu)化對比Tab.2 Comparison of Model Optimization

    通過表2 表明:基于最佳樣本空間(OSF)試驗設計構(gòu)建的Kriging 代理模型進行優(yōu)化相對于NSGA-II 對真實模型優(yōu)化調(diào)用次數(shù)較少65次,優(yōu)化時間減少了12.992h,優(yōu)化效率提升明顯。NSGA-II 多目標優(yōu)化算法對參數(shù)進行設置[13],初始種群規(guī)模為300,最大允許的Pareto百分比為70%,進化代數(shù)300,突變因子概率為0.02,交叉因子概率為0.98。通過對Kriging擬合代理模型進行自適應多目標優(yōu)化,得到的為Pareto解非支配排序散點圖,如圖15所示。

    圖15 Pareto解集散點圖Fig.15 Pareto Solution Scatter Plot

    通過在Pareto非支配排序最優(yōu)解集中選出5個優(yōu)化設計方案中,通過對其圓整后,如表2所示。

    由表3可知在Pareto最優(yōu)解集中的5個優(yōu)化設計方案中可以看出目標函數(shù)δmax,σmax,與mass之間存在某種制約關系,因而選擇出一個折中方案。方案1 相對于剩余方案質(zhì)量較少了(5.4~8.7)%,而方案一的應力和位移相對于其他方案相差不大,因此選擇方案1作為堆垛機多目標優(yōu)化的最優(yōu)解。通過表四最優(yōu)解與初始解進行對比可知,其堆垛機整機質(zhì)量下降了5%,最大位移減小了2.7%,最大應力增加了14%,卻依然滿足靜強度和靜剛度要求。

    表3 5個最優(yōu)Pareto解Tab.3 5 Optimal Pareto Solutions

    表4 最優(yōu)方案和原始方案對比Tab.4 Comparison of Optimal and Original Plans

    6 結(jié)論

    這里以雙高雙立柱堆垛機整體結(jié)構(gòu)為工程實例,提出以三維模型的參數(shù),最佳填充樣本空間(OSF)、Kriging模型非支配排序多目標遺傳算法(NSGA-II)相結(jié)合的多目標靜態(tài)優(yōu)化設計方法,其特點為:(1)在雙高雙立柱堆垛機多目標優(yōu)化過程中,通過采用最佳樣本填充空間(OSF)進行試驗設計,在通過Kriging模型構(gòu)建近似響應面模型代替復雜真實的模型。從而避免優(yōu)化時調(diào)用真實模型的次數(shù),提高優(yōu)化效率。(2)通過考慮最大應力、最大位移和整機質(zhì)量等要求,在雙高雙立柱堆垛機結(jié)構(gòu)多目標優(yōu)化過程中,通過運用非支配排序多目標遺傳算法(NSGA-II)進行尋優(yōu)獲得多個最優(yōu)的非支配排序Pareto解。(3)通過上述方法能夠有效的提高結(jié)構(gòu)設計的效率和質(zhì)量,對雙高雙立柱式堆垛機結(jié)構(gòu)設計提供了參考。

    猜你喜歡
    堆垛雙高橫梁
    同頻共振四鏈聯(lián)動,彰顯“雙高”院??蒲袚?/a>
    搬易通推出MCC系列人上型三向堆垛車
    我校2020年“雙高計劃”建設研究成果展
    軟橫跨改軟索式硬橫梁研究
    我校隆重召開“雙高計劃”建設啟動大會
    牽手教育: 中等職業(yè)學?!八脑⑴e, 德技雙高” 的校本實踐
    廣東教育(2020年3期)2020-04-03 02:03:18
    基于鋼橫梁剛度折減的鋼板組合梁地震響應分析
    工程與建設(2019年4期)2019-10-10 01:45:46
    自動化立體倉庫用堆垛機的幾種換軌方式及應用案例
    立式車床高精度、高剛度進給式橫梁設計
    淺談框架式玻璃幕墻閉腔橫梁系統(tǒng)
    在线观看免费视频网站a站| 欧美变态另类bdsm刘玥| 在线精品无人区一区二区三 | 欧美日韩在线观看h| 欧美bdsm另类| 有码 亚洲区| 男女国产视频网站| 国产精品熟女久久久久浪| 国产亚洲一区二区精品| 久久久久久久久久久丰满| 国产在线男女| videos熟女内射| 国产精品不卡视频一区二区| 人体艺术视频欧美日本| 精华霜和精华液先用哪个| 一级片'在线观看视频| 天美传媒精品一区二区| 人妻系列 视频| 欧美极品一区二区三区四区| 精品一品国产午夜福利视频| 一个人看视频在线观看www免费| 人妻制服诱惑在线中文字幕| 亚洲国产精品999| 国产成人精品福利久久| 国产成人一区二区在线| 免费av不卡在线播放| 午夜福利高清视频| 亚洲自偷自拍三级| 日本一二三区视频观看| 天堂俺去俺来也www色官网| 乱码一卡2卡4卡精品| 亚洲精华国产精华液的使用体验| 欧美日韩综合久久久久久| freevideosex欧美| 亚洲一级一片aⅴ在线观看| 欧美一区二区亚洲| 水蜜桃什么品种好| 最近最新中文字幕大全电影3| 亚洲性久久影院| 久久人妻熟女aⅴ| 日本色播在线视频| 搡老乐熟女国产| 美女cb高潮喷水在线观看| 国产一级毛片在线| 免费看不卡的av| 大陆偷拍与自拍| 一区二区三区精品91| 国产精品人妻久久久影院| 欧美一区二区亚洲| 狂野欧美激情性bbbbbb| 成人一区二区视频在线观看| 熟女电影av网| 国产精品一及| 黄色一级大片看看| 日韩一本色道免费dvd| 最近最新中文字幕免费大全7| 午夜免费鲁丝| 王馨瑶露胸无遮挡在线观看| 欧美极品一区二区三区四区| 午夜免费观看性视频| .国产精品久久| 国产真实伦视频高清在线观看| 又粗又硬又长又爽又黄的视频| 男男h啪啪无遮挡| 亚洲精品日本国产第一区| 九九爱精品视频在线观看| 国产伦在线观看视频一区| 男女免费视频国产| 精品午夜福利在线看| 天天躁日日操中文字幕| 国产69精品久久久久777片| 欧美极品一区二区三区四区| 国产欧美亚洲国产| 美女主播在线视频| 啦啦啦中文免费视频观看日本| 最近2019中文字幕mv第一页| 91精品国产九色| 国产精品国产av在线观看| 国产有黄有色有爽视频| 亚洲国产精品专区欧美| 狂野欧美激情性xxxx在线观看| 国产精品久久久久成人av| 久久久久久人妻| 国产伦在线观看视频一区| 亚洲激情五月婷婷啪啪| 黑人高潮一二区| 成人影院久久| 国产视频内射| 熟女av电影| 免费人妻精品一区二区三区视频| 欧美日韩视频高清一区二区三区二| 美女高潮的动态| 国产成人freesex在线| 日韩伦理黄色片| 亚洲成人手机| 美女中出高潮动态图| 午夜福利高清视频| 女人十人毛片免费观看3o分钟| 极品教师在线视频| 精品国产一区二区三区久久久樱花 | 亚洲欧美一区二区三区黑人 | 91久久精品国产一区二区三区| 亚洲va在线va天堂va国产| 亚洲精品456在线播放app| 午夜免费男女啪啪视频观看| 中文字幕av成人在线电影| 丰满人妻一区二区三区视频av| 毛片一级片免费看久久久久| 成年av动漫网址| 精品少妇黑人巨大在线播放| 哪个播放器可以免费观看大片| 噜噜噜噜噜久久久久久91| 在线免费十八禁| 成人综合一区亚洲| 永久免费av网站大全| 777米奇影视久久| 搡老乐熟女国产| 在线 av 中文字幕| 午夜视频国产福利| 永久免费av网站大全| 亚洲美女视频黄频| 亚洲综合精品二区| 精品一区二区三卡| 久久久亚洲精品成人影院| 十分钟在线观看高清视频www | 欧美日韩视频高清一区二区三区二| 韩国av在线不卡| 免费看av在线观看网站| 精品酒店卫生间| 哪个播放器可以免费观看大片| 少妇人妻久久综合中文| 亚洲欧美清纯卡通| 久热这里只有精品99| 男女免费视频国产| 免费播放大片免费观看视频在线观看| 国模一区二区三区四区视频| 久久国产精品男人的天堂亚洲 | 婷婷色综合大香蕉| 一区二区三区四区激情视频| 国产一级毛片在线| 老司机影院毛片| 国产伦精品一区二区三区视频9| 人人妻人人添人人爽欧美一区卜 | 看十八女毛片水多多多| 精品久久久久久久久亚洲| 免费大片18禁| 纯流量卡能插随身wifi吗| 老师上课跳d突然被开到最大视频| 午夜免费男女啪啪视频观看| 亚洲美女黄色视频免费看| 国产精品爽爽va在线观看网站| 最近的中文字幕免费完整| 在线观看免费视频网站a站| 插逼视频在线观看| 97超视频在线观看视频| 久久久久网色| 赤兔流量卡办理| 亚洲精品乱码久久久久久按摩| 久久人人爽人人爽人人片va| 久久精品夜色国产| 国产毛片在线视频| 黄片wwwwww| 免费观看无遮挡的男女| 蜜桃亚洲精品一区二区三区| 老司机影院毛片| 99热网站在线观看| 五月天丁香电影| 在线观看一区二区三区激情| 亚洲欧美精品自产自拍| 国产男女内射视频| 最近的中文字幕免费完整| 国产午夜精品一二区理论片| 亚洲成人中文字幕在线播放| 亚洲成人手机| av在线老鸭窝| 国产精品福利在线免费观看| 亚洲婷婷狠狠爱综合网| 97在线人人人人妻| 夜夜看夜夜爽夜夜摸| 在线观看三级黄色| 国产中年淑女户外野战色| 亚洲综合精品二区| 国产精品熟女久久久久浪| 人妻 亚洲 视频| 久久久久久九九精品二区国产| 多毛熟女@视频| 欧美xxxx性猛交bbbb| 一级毛片aaaaaa免费看小| 91久久精品电影网| av一本久久久久| a级一级毛片免费在线观看| 欧美三级亚洲精品| 免费观看无遮挡的男女| 欧美高清成人免费视频www| 久久精品国产自在天天线| 国产精品久久久久成人av| 中文资源天堂在线| 国产国拍精品亚洲av在线观看| 韩国av在线不卡| 国产男女超爽视频在线观看| 青青草视频在线视频观看| 极品少妇高潮喷水抽搐| 自拍偷自拍亚洲精品老妇| 午夜老司机福利剧场| 另类亚洲欧美激情| 久久久久久久久大av| 97精品久久久久久久久久精品| 少妇高潮的动态图| 搡老乐熟女国产| 干丝袜人妻中文字幕| 亚洲欧美成人精品一区二区| 亚洲性久久影院| 国产精品嫩草影院av在线观看| 尤物成人国产欧美一区二区三区| 日韩在线高清观看一区二区三区| 午夜福利在线在线| 免费大片18禁| 欧美成人精品欧美一级黄| 国产淫语在线视频| 亚洲久久久国产精品| 精品国产一区二区三区久久久樱花 | 最近中文字幕2019免费版| 国产成人免费观看mmmm| 国产精品熟女久久久久浪| av在线老鸭窝| 国产av精品麻豆| 九九久久精品国产亚洲av麻豆| 少妇猛男粗大的猛烈进出视频| 蜜桃亚洲精品一区二区三区| av福利片在线观看| 成人国产av品久久久| 精华霜和精华液先用哪个| av又黄又爽大尺度在线免费看| 久久青草综合色| 建设人人有责人人尽责人人享有的 | 小蜜桃在线观看免费完整版高清| 午夜福利网站1000一区二区三区| 大片电影免费在线观看免费| 夜夜骑夜夜射夜夜干| 老女人水多毛片| 色网站视频免费| 一区二区av电影网| 天美传媒精品一区二区| 国产精品伦人一区二区| 国产欧美亚洲国产| 在线观看一区二区三区| 免费观看无遮挡的男女| 欧美高清性xxxxhd video| 亚洲经典国产精华液单| 高清av免费在线| 晚上一个人看的免费电影| 午夜福利网站1000一区二区三区| 午夜福利影视在线免费观看| 亚洲经典国产精华液单| 亚洲,一卡二卡三卡| 成年女人在线观看亚洲视频| 国产精品人妻久久久影院| 欧美丝袜亚洲另类| 亚洲精品自拍成人| av不卡在线播放| 免费看光身美女| 2018国产大陆天天弄谢| 日日啪夜夜爽| 久久久久精品久久久久真实原创| 亚洲精品中文字幕在线视频 | 精品久久久久久久久亚洲| 乱系列少妇在线播放| 99精国产麻豆久久婷婷| 国产成人免费无遮挡视频| 在线播放无遮挡| 亚州av有码| 久久久久久久大尺度免费视频| 午夜福利高清视频| 麻豆成人午夜福利视频| 亚洲第一区二区三区不卡| 高清不卡的av网站| 欧美激情国产日韩精品一区| 国产精品精品国产色婷婷| 国产一区有黄有色的免费视频| 五月伊人婷婷丁香| 国产精品伦人一区二区| 黄片无遮挡物在线观看| 亚洲精品成人av观看孕妇| 亚洲欧美精品自产自拍| 美女cb高潮喷水在线观看| 国产精品.久久久| 亚洲av电影在线观看一区二区三区| 中文字幕久久专区| 亚洲自偷自拍三级| 一级爰片在线观看| 我的女老师完整版在线观看| av线在线观看网站| 免费人妻精品一区二区三区视频| 久久人人爽人人片av| 高清黄色对白视频在线免费看 | 乱码一卡2卡4卡精品| 精品国产乱码久久久久久小说| 国产乱来视频区| 成人毛片a级毛片在线播放| 久久午夜福利片| 欧美激情极品国产一区二区三区 | av福利片在线观看| 欧美日韩视频高清一区二区三区二| 久久久精品免费免费高清| 最近最新中文字幕免费大全7| 国产在视频线精品| 国产一区二区三区综合在线观看 | 国产精品一区二区在线观看99| 天美传媒精品一区二区| 精品少妇黑人巨大在线播放| 亚洲精品国产色婷婷电影| 久久久a久久爽久久v久久| 又大又黄又爽视频免费| 一级av片app| 日韩av在线免费看完整版不卡| 免费观看的影片在线观看| h日本视频在线播放| 天美传媒精品一区二区| 人人妻人人爽人人添夜夜欢视频 | 美女xxoo啪啪120秒动态图| 国产免费又黄又爽又色| 一本—道久久a久久精品蜜桃钙片| 久久久久久人妻| 国产免费一级a男人的天堂| 黄色怎么调成土黄色| 免费看av在线观看网站| 日韩大片免费观看网站| 51国产日韩欧美| 国产探花极品一区二区| 日韩三级伦理在线观看| 在线亚洲精品国产二区图片欧美 | 欧美丝袜亚洲另类| 久久国产精品大桥未久av | 国产av国产精品国产| 中文在线观看免费www的网站| 一区在线观看完整版| 日韩 亚洲 欧美在线| 肉色欧美久久久久久久蜜桃| 久久久久久久精品精品| 精品一区二区三卡| 又大又黄又爽视频免费| 少妇人妻精品综合一区二区| 欧美最新免费一区二区三区| 又粗又硬又长又爽又黄的视频| 啦啦啦啦在线视频资源| 亚洲电影在线观看av| 欧美日韩视频高清一区二区三区二| 国产亚洲av片在线观看秒播厂| 久久久午夜欧美精品| 国产成人a∨麻豆精品| 又爽又黄a免费视频| 国产有黄有色有爽视频| 亚洲精品国产成人久久av| 三级国产精品片| 欧美日韩综合久久久久久| 欧美另类一区| 中文字幕亚洲精品专区| 制服丝袜香蕉在线| 成人影院久久| 精品亚洲乱码少妇综合久久| 成人毛片a级毛片在线播放| 天美传媒精品一区二区| 在线观看人妻少妇| 国产黄色免费在线视频| 亚洲成人手机| 亚洲精华国产精华液的使用体验| 美女xxoo啪啪120秒动态图| 熟女av电影| 99九九线精品视频在线观看视频| 国产毛片在线视频| 国产综合精华液| 色5月婷婷丁香| 一级二级三级毛片免费看| 人妻 亚洲 视频| 丰满乱子伦码专区| 麻豆成人av视频| 美女xxoo啪啪120秒动态图| 人妻一区二区av| 午夜免费鲁丝| 性高湖久久久久久久久免费观看| 亚洲高清免费不卡视频| 亚洲电影在线观看av| 久久国内精品自在自线图片| 身体一侧抽搐| 国产精品一区www在线观看| 男的添女的下面高潮视频| 欧美精品亚洲一区二区| 亚洲精品,欧美精品| 干丝袜人妻中文字幕| 免费观看a级毛片全部| 在线亚洲精品国产二区图片欧美 | 这个男人来自地球电影免费观看 | 日韩 亚洲 欧美在线| 亚洲精品亚洲一区二区| 美女高潮的动态| 一区在线观看完整版| 中文字幕免费在线视频6| 99热国产这里只有精品6| 亚洲不卡免费看| 日日撸夜夜添| 成年女人在线观看亚洲视频| 国产精品一区二区在线不卡| 亚洲欧美日韩东京热| 在线观看免费高清a一片| 老司机影院毛片| 我要看黄色一级片免费的| 成人毛片60女人毛片免费| 日韩 亚洲 欧美在线| 国产精品国产三级专区第一集| 热re99久久精品国产66热6| 成人黄色视频免费在线看| 亚洲国产成人一精品久久久| 日日撸夜夜添| 国产欧美日韩精品一区二区| 一级毛片久久久久久久久女| 伦理电影大哥的女人| 亚洲婷婷狠狠爱综合网| 亚洲av日韩在线播放| 亚洲欧美一区二区三区国产| 只有这里有精品99| 蜜桃久久精品国产亚洲av| 精品久久久久久电影网| 免费少妇av软件| 在线观看三级黄色| 欧美日本视频| 高清不卡的av网站| 老司机影院毛片| 色婷婷久久久亚洲欧美| 日韩av在线免费看完整版不卡| 亚洲欧美中文字幕日韩二区| freevideosex欧美| 在线亚洲精品国产二区图片欧美 | 日韩伦理黄色片| 精品人妻一区二区三区麻豆| 国产男人的电影天堂91| 另类亚洲欧美激情| 国产高清有码在线观看视频| 国产男女超爽视频在线观看| 久久久久精品性色| 汤姆久久久久久久影院中文字幕| 国产美女午夜福利| 成人毛片a级毛片在线播放| 一级黄片播放器| 国产精品成人在线| 久久精品国产a三级三级三级| 黑人高潮一二区| 国产成人精品一,二区| 亚洲国产av新网站| 直男gayav资源| 亚洲av福利一区| 日本欧美视频一区| 日韩三级伦理在线观看| 在线观看一区二区三区激情| 高清av免费在线| 中国美白少妇内射xxxbb| 女人十人毛片免费观看3o分钟| 美女脱内裤让男人舔精品视频| 国产淫片久久久久久久久| 精品熟女少妇av免费看| 黄色配什么色好看| 国精品久久久久久国模美| 亚洲在久久综合| 91在线精品国自产拍蜜月| 国产伦理片在线播放av一区| 色哟哟·www| 精品视频人人做人人爽| 2018国产大陆天天弄谢| 欧美变态另类bdsm刘玥| 激情 狠狠 欧美| 99国产精品免费福利视频| 久久97久久精品| 三级经典国产精品| 视频中文字幕在线观看| 久久鲁丝午夜福利片| av在线播放精品| 在线精品无人区一区二区三 | 熟女人妻精品中文字幕| 在线天堂最新版资源| 亚洲aⅴ乱码一区二区在线播放| 亚洲欧美日韩东京热| 亚洲aⅴ乱码一区二区在线播放| 99re6热这里在线精品视频| 我要看日韩黄色一级片| 免费人成在线观看视频色| 色婷婷av一区二区三区视频| 亚洲人成网站在线播| 亚洲成人中文字幕在线播放| 天堂俺去俺来也www色官网| 日韩一本色道免费dvd| 在线免费十八禁| 男女免费视频国产| 日韩亚洲欧美综合| 精品国产一区二区三区久久久樱花 | 国产黄片视频在线免费观看| 色5月婷婷丁香| 美女xxoo啪啪120秒动态图| 亚洲精品一二三| 激情 狠狠 欧美| 亚洲欧美成人精品一区二区| 国产精品爽爽va在线观看网站| 久久久精品免费免费高清| 国产亚洲午夜精品一区二区久久| 大又大粗又爽又黄少妇毛片口| 亚洲av福利一区| 成人美女网站在线观看视频| 国产精品爽爽va在线观看网站| 五月开心婷婷网| 老司机影院成人| 亚洲熟女精品中文字幕| 免费少妇av软件| 观看美女的网站| 亚洲国产高清在线一区二区三| 日韩一区二区三区影片| 日本wwww免费看| 成人午夜精彩视频在线观看| 国精品久久久久久国模美| 国产精品一区二区性色av| 亚洲精品中文字幕在线视频 | 久久久久人妻精品一区果冻| 亚洲成人av在线免费| 熟女人妻精品中文字幕| 精品国产三级普通话版| 国产av一区二区精品久久 | 久久久久久久久久成人| 毛片女人毛片| 街头女战士在线观看网站| 久久久a久久爽久久v久久| 欧美性感艳星| 久久人妻熟女aⅴ| 免费看av在线观看网站| av又黄又爽大尺度在线免费看| 麻豆成人av视频| 国产精品久久久久成人av| 黄片无遮挡物在线观看| 一级爰片在线观看| 久久精品久久精品一区二区三区| 亚洲欧美精品自产自拍| 欧美成人精品欧美一级黄| 又粗又硬又长又爽又黄的视频| 亚洲欧美日韩东京热| 国产老妇伦熟女老妇高清| 亚洲天堂av无毛| 大又大粗又爽又黄少妇毛片口| 一区二区三区四区激情视频| av网站免费在线观看视频| 久久久久久人妻| av一本久久久久| 亚洲成人av在线免费| 国产乱来视频区| 中国三级夫妇交换| 人人妻人人爽人人添夜夜欢视频 | 成年免费大片在线观看| 亚州av有码| 国产成人freesex在线| 男人爽女人下面视频在线观看| 精品国产乱码久久久久久小说| 在线观看免费视频网站a站| 午夜福利视频精品| 这个男人来自地球电影免费观看 | 精华霜和精华液先用哪个| 日本爱情动作片www.在线观看| 不卡视频在线观看欧美| 亚洲婷婷狠狠爱综合网| 精品熟女少妇av免费看| 欧美精品亚洲一区二区| 国产精品成人在线| 蜜桃久久精品国产亚洲av| 久久久久视频综合| 亚洲av中文字字幕乱码综合| 美女xxoo啪啪120秒动态图| 国产一区二区三区综合在线观看 | 97热精品久久久久久| 亚洲欧美精品自产自拍| 日韩 亚洲 欧美在线| 国产精品成人在线| 国产爽快片一区二区三区| 亚洲av男天堂| 亚洲丝袜综合中文字幕| 丰满迷人的少妇在线观看| 岛国毛片在线播放| 97热精品久久久久久| 日本黄大片高清| 日韩 亚洲 欧美在线| 日日啪夜夜撸| 超碰97精品在线观看| 国产伦精品一区二区三区视频9| 国产高清三级在线| 男女国产视频网站| 国产精品不卡视频一区二区| 啦啦啦中文免费视频观看日本| 久久久久久久久久人人人人人人| 99久久精品国产国产毛片| 国产精品成人在线| 在线观看免费视频网站a站| 能在线免费看毛片的网站| 国产 一区精品| 国产欧美日韩精品一区二区| 天天躁日日操中文字幕| 亚洲欧美精品专区久久| 日韩视频在线欧美| 久久午夜福利片| 久久国产乱子免费精品| av专区在线播放| 性色avwww在线观看| 成人黄色视频免费在线看| 国产av一区二区精品久久 | 性高湖久久久久久久久免费观看| 久久久久性生活片| av又黄又爽大尺度在线免费看| 成人免费观看视频高清| 亚洲欧美日韩另类电影网站 | 亚洲精品一区蜜桃| 伦理电影免费视频| 日本黄大片高清| 久久精品国产鲁丝片午夜精品|