• <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)
    国产蜜桃级精品一区二区三区| 免费黄网站久久成人精品| 黄片无遮挡物在线观看| 99热精品在线国产| 国产精品久久久久久av不卡| 丰满的人妻完整版| 三级经典国产精品| 乱码一卡2卡4卡精品| 一边摸一边抽搐一进一小说| 久久精品国产自在天天线| 国产探花极品一区二区| 国产精品伦人一区二区| 欧美精品国产亚洲| 欧美日本亚洲视频在线播放| 一级黄色大片毛片| 伦理电影大哥的女人| 日本一本二区三区精品| 国产视频内射| 日韩欧美精品免费久久| 亚洲真实伦在线观看| 99热精品在线国产| 热99re8久久精品国产| 亚洲婷婷狠狠爱综合网| 国产蜜桃级精品一区二区三区| 亚洲激情五月婷婷啪啪| 亚洲成av人片在线播放无| 成人亚洲精品av一区二区| 亚洲五月天丁香| АⅤ资源中文在线天堂| 亚洲欧洲国产日韩| 亚洲人成网站在线播放欧美日韩| 天堂影院成人在线观看| 国产中年淑女户外野战色| АⅤ资源中文在线天堂| 国产黄片美女视频| 别揉我奶头 嗯啊视频| 精品人妻一区二区三区麻豆| 国产成人aa在线观看| 99久久无色码亚洲精品果冻| 日本爱情动作片www.在线观看| 高清午夜精品一区二区三区 | 一进一出抽搐gif免费好疼| 在线观看免费视频日本深夜| 夫妻性生交免费视频一级片| 精品久久久久久成人av| av女优亚洲男人天堂| 亚洲精品久久久久久婷婷小说 | 久久久久网色| 99久久精品热视频| 国产精品乱码一区二三区的特点| 成人亚洲精品av一区二区| 色播亚洲综合网| 久久久精品94久久精品| 天堂√8在线中文| 国产高清视频在线观看网站| 午夜精品一区二区三区免费看| 欧美激情在线99| 黄色视频,在线免费观看| 亚洲第一区二区三区不卡| 一级毛片aaaaaa免费看小| 18禁在线播放成人免费| 99热只有精品国产| 日本爱情动作片www.在线观看| 日本黄色视频三级网站网址| 国产在线男女| 精品久久久久久成人av| 免费在线观看成人毛片| 亚洲高清免费不卡视频| 99九九线精品视频在线观看视频| 亚洲精品色激情综合| 国产 一区精品| 日韩欧美 国产精品| 国产精品一区二区三区四区免费观看| 精品久久久久久久人妻蜜臀av| 国产91av在线免费观看| 午夜激情福利司机影院| 日韩av在线大香蕉| 中文亚洲av片在线观看爽| 中文资源天堂在线| 国产高清三级在线| 干丝袜人妻中文字幕| 成人无遮挡网站| 亚洲最大成人av| 97超视频在线观看视频| 亚洲国产精品sss在线观看| 边亲边吃奶的免费视频| 波野结衣二区三区在线| 国产午夜精品一二区理论片| 丝袜喷水一区| 一级av片app| 99在线视频只有这里精品首页| 久久韩国三级中文字幕| 中文字幕熟女人妻在线| 亚洲欧洲日产国产| 亚洲综合色惰| 3wmmmm亚洲av在线观看| 亚洲真实伦在线观看| 男女啪啪激烈高潮av片| 麻豆久久精品国产亚洲av| 网址你懂的国产日韩在线| kizo精华| 在线播放国产精品三级| 99久久成人亚洲精品观看| av黄色大香蕉| 久久精品国产亚洲av天美| 少妇熟女欧美另类| 一本久久中文字幕| 18禁在线无遮挡免费观看视频| 亚洲国产精品成人综合色| 免费搜索国产男女视频| 超碰av人人做人人爽久久| 床上黄色一级片| 长腿黑丝高跟| 国产一区二区在线av高清观看| 亚洲av成人av| 亚洲av.av天堂| 久久99精品国语久久久| 国产真实伦视频高清在线观看| 国产av在哪里看| 国产午夜福利久久久久久| 欧美激情国产日韩精品一区| 国产精品野战在线观看| 美女被艹到高潮喷水动态| 国产伦精品一区二区三区四那| 男女那种视频在线观看| 又爽又黄a免费视频| 99久久人妻综合| 亚洲三级黄色毛片| 国内精品久久久久精免费| 中文字幕人妻熟人妻熟丝袜美| 亚洲美女搞黄在线观看| 精品日产1卡2卡| 国产日韩欧美在线精品| 99久久成人亚洲精品观看| 麻豆精品久久久久久蜜桃| 亚洲天堂国产精品一区在线| 成人午夜精彩视频在线观看| 中文资源天堂在线| 亚洲第一区二区三区不卡| 深夜a级毛片| 九九热线精品视视频播放| 色综合亚洲欧美另类图片| 午夜a级毛片| 亚洲自偷自拍三级| 日韩欧美精品v在线| 91aial.com中文字幕在线观看| 亚洲精品成人久久久久久| 久久99精品国语久久久| 亚洲中文字幕一区二区三区有码在线看| h日本视频在线播放| 日韩欧美 国产精品| 中文字幕av成人在线电影| av福利片在线观看| 久久鲁丝午夜福利片| 中文字幕av在线有码专区| 国产精品国产高清国产av| 亚洲av二区三区四区| 老师上课跳d突然被开到最大视频| 亚洲欧美日韩东京热| 国产欧美日韩精品一区二区| 天堂av国产一区二区熟女人妻| 国产成人freesex在线| 国产老妇女一区| 国产伦一二天堂av在线观看| 天天一区二区日本电影三级| 亚洲中文字幕日韩| 白带黄色成豆腐渣| 三级男女做爰猛烈吃奶摸视频| 不卡视频在线观看欧美| 婷婷亚洲欧美| 成人高潮视频无遮挡免费网站| 免费无遮挡裸体视频| 成人亚洲欧美一区二区av| 午夜福利在线观看免费完整高清在 | 日韩av在线大香蕉| 天美传媒精品一区二区| а√天堂www在线а√下载| 日日干狠狠操夜夜爽| 免费大片18禁| 免费人成视频x8x8入口观看| 夜夜爽天天搞| 在线观看午夜福利视频| 一区二区三区免费毛片| 亚洲精品影视一区二区三区av| 能在线免费看毛片的网站| 日本欧美国产在线视频| 能在线免费观看的黄片| 一本一本综合久久| 床上黄色一级片| 亚洲精品成人久久久久久| 色噜噜av男人的天堂激情| 久久久久久久久大av| 日日撸夜夜添| 亚洲av成人av| 国产一区二区激情短视频| 激情 狠狠 欧美| 69人妻影院| 岛国在线免费视频观看| 嫩草影院新地址| 嫩草影院入口| 欧美zozozo另类| 日韩av在线大香蕉| 中国美女看黄片| 18禁在线无遮挡免费观看视频| 亚洲最大成人手机在线| 韩国av在线不卡| 韩国av在线不卡| 精品一区二区三区视频在线| 精品一区二区三区视频在线| 麻豆乱淫一区二区| 桃色一区二区三区在线观看| 99热这里只有是精品在线观看| 直男gayav资源| 亚洲在线观看片| 亚洲国产日韩欧美精品在线观看| 久久精品国产亚洲av涩爱 | 久久中文看片网| 有码 亚洲区| 女人十人毛片免费观看3o分钟| 国产亚洲av嫩草精品影院| 欧美xxxx性猛交bbbb| 亚洲图色成人| 婷婷亚洲欧美| 亚洲精品自拍成人| 国产一区二区三区av在线 | 禁无遮挡网站| 国产精品伦人一区二区| 精品久久久久久久久久免费视频| 欧美变态另类bdsm刘玥| 国语自产精品视频在线第100页| 寂寞人妻少妇视频99o| 亚洲真实伦在线观看| 日韩制服骚丝袜av| 狂野欧美激情性xxxx在线观看| 丝袜美腿在线中文| 最近的中文字幕免费完整| 久久久久久久久久久丰满| 啦啦啦韩国在线观看视频| 九草在线视频观看| 又粗又爽又猛毛片免费看| 日韩欧美一区二区三区在线观看| 国产精品一及| 99在线人妻在线中文字幕| 亚州av有码| 村上凉子中文字幕在线| 国产成人福利小说| 欧美性猛交黑人性爽| 亚洲高清免费不卡视频| 麻豆av噜噜一区二区三区| 亚洲激情五月婷婷啪啪| 成人三级黄色视频| 国产一区二区亚洲精品在线观看| 欧美zozozo另类| 如何舔出高潮| 国产精品人妻久久久久久| 国产精品1区2区在线观看.| 男的添女的下面高潮视频| 亚洲一区二区三区色噜噜| 欧美一级a爱片免费观看看| 国产伦在线观看视频一区| 成年女人看的毛片在线观看| 最好的美女福利视频网| 亚洲内射少妇av| 午夜福利在线观看吧| 国产精品久久久久久久久免| 天天一区二区日本电影三级| 男女下面进入的视频免费午夜| 97热精品久久久久久| 干丝袜人妻中文字幕| 欧美+日韩+精品| 婷婷色综合大香蕉| 日日摸夜夜添夜夜爱| 色尼玛亚洲综合影院| 网址你懂的国产日韩在线| 午夜福利高清视频| avwww免费| 狂野欧美激情性xxxx在线观看| 国产成人a∨麻豆精品| 人妻系列 视频| 国产在线精品亚洲第一网站| 国产伦理片在线播放av一区 | 亚洲国产高清在线一区二区三| 国产午夜精品久久久久久一区二区三区| 久久99热这里只有精品18| or卡值多少钱| 午夜精品在线福利| 国产91av在线免费观看| 成人午夜精彩视频在线观看| 在线观看66精品国产| 蜜臀久久99精品久久宅男| 国产精品久久久久久精品电影| or卡值多少钱| 久久国内精品自在自线图片| 久久精品国产99精品国产亚洲性色| 观看免费一级毛片| 欧美激情国产日韩精品一区| 一区二区三区高清视频在线| 国产伦精品一区二区三区视频9| 最近视频中文字幕2019在线8| av卡一久久| 日韩一区二区三区影片| 久久人妻av系列| 黄色配什么色好看| 日韩av不卡免费在线播放| 日韩欧美三级三区| 国产高清有码在线观看视频| 国产精品一区www在线观看| 精品久久久久久久久久久久久| 人妻制服诱惑在线中文字幕| 国内少妇人妻偷人精品xxx网站| 热99在线观看视频| 综合色丁香网| 国产中年淑女户外野战色| 国产成人精品一,二区 | 长腿黑丝高跟| av视频在线观看入口| 国产视频内射| 99热精品在线国产| 中文亚洲av片在线观看爽| 草草在线视频免费看| 最近中文字幕高清免费大全6| 最近2019中文字幕mv第一页| 精品久久久噜噜| 欧美日韩综合久久久久久| 亚洲精品成人久久久久久| 深爱激情五月婷婷| 国产成人aa在线观看| 日本一本二区三区精品| 国产精品.久久久| 国产av在哪里看| 亚洲欧洲日产国产| 免费观看的影片在线观看| 亚洲丝袜综合中文字幕| 亚洲经典国产精华液单| 久久99精品国语久久久| 床上黄色一级片| 久久九九热精品免费| 日本在线视频免费播放| 久久久久国产网址| 亚洲乱码一区二区免费版| 噜噜噜噜噜久久久久久91| 国产成人aa在线观看| 日韩视频在线欧美| 成人亚洲欧美一区二区av| 在线观看一区二区三区| 91久久精品国产一区二区三区| 1024手机看黄色片| 国产又黄又爽又无遮挡在线| 婷婷亚洲欧美| 五月伊人婷婷丁香| 69人妻影院| 免费看日本二区| 欧美变态另类bdsm刘玥| 久久精品人妻少妇| 日本与韩国留学比较| 黑人高潮一二区| 亚洲精品影视一区二区三区av| 伊人久久精品亚洲午夜| 一区二区三区高清视频在线| 波多野结衣高清无吗| 国产国拍精品亚洲av在线观看| 中文欧美无线码| 精品久久国产蜜桃| 久99久视频精品免费| 亚洲丝袜综合中文字幕| ponron亚洲| 成年女人永久免费观看视频| 在现免费观看毛片| 天堂网av新在线| 欧美日韩一区二区视频在线观看视频在线 | 18禁在线无遮挡免费观看视频| 国产精品久久电影中文字幕| 大香蕉久久网| 亚洲欧洲国产日韩| 久久午夜福利片| 在线播放国产精品三级| 国产精品一区二区在线观看99 | 日韩av在线大香蕉| 99久久成人亚洲精品观看| 精品不卡国产一区二区三区| 久久精品国产亚洲av香蕉五月| 好男人视频免费观看在线| 丝袜喷水一区| 一区二区三区高清视频在线| 亚洲内射少妇av| 少妇熟女欧美另类| 国产精品不卡视频一区二区| 国产真实伦视频高清在线观看| 青青草视频在线视频观看| 女人被狂操c到高潮| 国产免费男女视频| 97人妻精品一区二区三区麻豆| 久久综合国产亚洲精品| 男人狂女人下面高潮的视频| 国产精品乱码一区二三区的特点| 精品久久久久久久久亚洲| а√天堂www在线а√下载| 成人亚洲精品av一区二区| av天堂中文字幕网| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 波多野结衣高清无吗| 我要搜黄色片| 免费黄网站久久成人精品| 三级男女做爰猛烈吃奶摸视频| 高清日韩中文字幕在线| 国产成人影院久久av| 国产成人精品一,二区 | 午夜福利视频1000在线观看| 亚洲在线观看片| 日韩人妻高清精品专区| 狂野欧美激情性xxxx在线观看| 国产精品一二三区在线看| 国产高清不卡午夜福利| 欧美性猛交╳xxx乱大交人| 国产成人aa在线观看| 亚洲,欧美,日韩| 欧美bdsm另类| 六月丁香七月| 又黄又爽又刺激的免费视频.| 久久99精品国语久久久| 午夜亚洲福利在线播放| 国产久久久一区二区三区| 成人特级黄色片久久久久久久| 天堂√8在线中文| 久久精品国产亚洲网站| 国产黄片美女视频| 亚洲av男天堂| 麻豆国产av国片精品| 免费av观看视频| av国产免费在线观看| 精品久久久久久久末码| 国产日韩欧美在线精品| 国产极品天堂在线| 亚洲av第一区精品v没综合| 亚洲国产精品久久男人天堂| 久久精品91蜜桃| 中文字幕熟女人妻在线| 校园春色视频在线观看| 成人毛片60女人毛片免费| 青春草亚洲视频在线观看| 国产精品蜜桃在线观看 | 亚洲成人中文字幕在线播放| 中国美女看黄片| 看片在线看免费视频| 日本与韩国留学比较| 性插视频无遮挡在线免费观看| 波多野结衣巨乳人妻| 日韩欧美精品v在线| 国模一区二区三区四区视频| 亚洲三级黄色毛片| 黄片无遮挡物在线观看| 禁无遮挡网站| 国产精品一区二区三区四区免费观看| 国产成年人精品一区二区| 日本撒尿小便嘘嘘汇集6| 看十八女毛片水多多多| 亚洲最大成人av| 久久精品国产亚洲av天美| 亚洲欧美精品自产自拍| 久久久精品欧美日韩精品| 蜜桃久久精品国产亚洲av| 日韩人妻高清精品专区| 亚洲内射少妇av| 成人二区视频| 黄色日韩在线| 久久99热这里只有精品18| 69av精品久久久久久| 中文在线观看免费www的网站| 日本免费一区二区三区高清不卡| 一个人免费在线观看电影| 久久人人爽人人片av| 日韩一区二区视频免费看| 97人妻精品一区二区三区麻豆| 九九热线精品视视频播放| 久久99热6这里只有精品| 插逼视频在线观看| av天堂中文字幕网| 中文精品一卡2卡3卡4更新| 如何舔出高潮| 国产免费一级a男人的天堂| 国产精品女同一区二区软件| 国产亚洲av片在线观看秒播厂 | 国产一区二区在线观看日韩| 日本一二三区视频观看| av在线亚洲专区| 亚洲熟妇中文字幕五十中出| 尤物成人国产欧美一区二区三区| www日本黄色视频网| 久久精品夜色国产| 身体一侧抽搐| 亚洲精品久久国产高清桃花| 国产精品99久久久久久久久| 亚洲国产欧美人成| 高清在线视频一区二区三区 | 免费一级毛片在线播放高清视频| 男女做爰动态图高潮gif福利片| 麻豆精品久久久久久蜜桃| 国产极品精品免费视频能看的| 99国产精品一区二区蜜桃av| 久久久欧美国产精品| 亚洲av.av天堂| 一级av片app| 久久热精品热| 国产亚洲av嫩草精品影院| 人妻制服诱惑在线中文字幕| 国产精品精品国产色婷婷| 国产日韩欧美在线精品| 国产高清激情床上av| 国产日本99.免费观看| 悠悠久久av| 美女内射精品一级片tv| 日韩成人伦理影院| 99在线人妻在线中文字幕| 高清毛片免费看| 尾随美女入室| 亚洲国产欧美人成| 国产精品麻豆人妻色哟哟久久 | 国产亚洲av嫩草精品影院| 不卡视频在线观看欧美| 精品久久久久久久久av| 精品久久久久久久久久久久久| 久久国产乱子免费精品| 人妻系列 视频| 12—13女人毛片做爰片一| 亚洲五月天丁香| 欧美高清成人免费视频www| 五月伊人婷婷丁香| 日韩高清综合在线| 97人妻精品一区二区三区麻豆| 成人无遮挡网站| kizo精华| 国产精品,欧美在线| 欧美丝袜亚洲另类| 在线a可以看的网站| 国产精品精品国产色婷婷| 亚洲欧美日韩高清在线视频| 男女下面进入的视频免费午夜| 99热全是精品| 国产高潮美女av| 久久国内精品自在自线图片| 欧美激情久久久久久爽电影| 狂野欧美激情性xxxx在线观看| 尤物成人国产欧美一区二区三区| av免费观看日本| 少妇人妻精品综合一区二区 | 国产一级毛片七仙女欲春2| 男女视频在线观看网站免费| 免费一级毛片在线播放高清视频| 丰满人妻一区二区三区视频av| 精品久久久久久成人av| 精品一区二区三区人妻视频| 一夜夜www| 日日摸夜夜添夜夜添av毛片| 麻豆国产97在线/欧美| 一本精品99久久精品77| 久久亚洲精品不卡| 男人狂女人下面高潮的视频| 国产乱人视频| 久久久国产成人免费| 婷婷精品国产亚洲av| 夜夜爽天天搞| 婷婷精品国产亚洲av| 中国美女看黄片| 夫妻性生交免费视频一级片| 国产精品久久视频播放| 美女脱内裤让男人舔精品视频 | 国产精品美女特级片免费视频播放器| 亚洲中文字幕日韩| 亚洲最大成人手机在线| 精品久久久久久成人av| 人人妻人人澡欧美一区二区| 午夜爱爱视频在线播放| 亚洲中文字幕一区二区三区有码在线看| 中国美女看黄片| 免费无遮挡裸体视频| 日韩成人av中文字幕在线观看| av在线亚洲专区| or卡值多少钱| 在线播放无遮挡| 丰满乱子伦码专区| 联通29元200g的流量卡| 长腿黑丝高跟| 成人欧美大片| 色综合站精品国产| 中文资源天堂在线| 国产黄色小视频在线观看| 亚洲性久久影院| 一级毛片久久久久久久久女| 一个人看视频在线观看www免费| 久久韩国三级中文字幕| 欧美日韩综合久久久久久| 日日干狠狠操夜夜爽| 色视频www国产| 男女那种视频在线观看| 亚洲最大成人av| 看十八女毛片水多多多| 插阴视频在线观看视频| 亚洲欧美精品自产自拍| 51国产日韩欧美| 少妇被粗大猛烈的视频| 亚洲在线观看片| 我要搜黄色片| 国产高清不卡午夜福利| 国产真实乱freesex| 国产精品福利在线免费观看| 亚洲av不卡在线观看| 草草在线视频免费看| 啦啦啦观看免费观看视频高清| 亚洲精品国产成人久久av| 97超碰精品成人国产| .国产精品久久| 深爱激情五月婷婷| 熟妇人妻久久中文字幕3abv| 精品免费久久久久久久清纯| 免费看光身美女|