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

    權期望填充準則作用下的Kriging自適應建模及全局優(yōu)化算法

    2022-04-04 06:11:40彭行坤林成龍馬義中
    計算機集成制造系統(tǒng) 2022年3期
    關鍵詞:優(yōu)化模型

    彭行坤,林成龍,馬義中

    (南京理工大學 經(jīng)濟管理學院,江蘇 南京 210094)

    0 引言

    隨著定制化產(chǎn)品等營銷策略的出現(xiàn),消費者個性化需求使產(chǎn)品設計和生產(chǎn)具備小批量、定制化的特點,傳統(tǒng)實物試驗設計已不能滿足產(chǎn)品復雜程度高,更新速度快的系統(tǒng)設計要求[1]?;谟嬎銠C學科建模,數(shù)值計算和計算機軟硬件技術發(fā)展而來的高精度仿真模型被廣泛應用于產(chǎn)品和系統(tǒng)優(yōu)化設計,有效提升了傳統(tǒng)實物試驗設計的質(zhì)量和可靠度[2-3]。高精度仿真建模過程中的計算流體力學(Computational Fluid Dynamics, CFD)和計算固體力學(Computational Solid Mechanics, CSM)技術極大地減少了研發(fā)成本,但建模需耗費大量時間且因?qū)嶋H優(yōu)化問題中往往無明確的函數(shù)關系,給工程建模和優(yōu)化帶來了極大困難[4-5]。近些年來,僅依靠初始樣本就可近似模擬輸入變量及響應變量二者之間函數(shù)關系的代理模型受到關注,先后產(chǎn)生了諸如響應面建模、KRIGING模型、徑向基函數(shù)、支持向量回歸及組合建模等代理模型[2,5]。上述模型因有效減少優(yōu)化過程的計算量且具備較好的精度及穩(wěn)健性在工程領域中得廣泛應用[5-6]。其中,基于統(tǒng)計技術的Kriging代理模型因其能提供對未知試驗點預測不確定性的度量且在小樣本數(shù)據(jù)條件下具有高的預測精度受到眾多科研工作者的青睞。

    Kriging代理模型自1951年Krige博士提出以來,衍生出了梯度增強Kriging模型、Co Kriging模型、Blind Kriging模型等眾多改進版本[7-8],尤其在Jones等[9-10]基于統(tǒng)計假設提出最大化期望改進(Expected Improvement, EI)填充策略,對Kriging模型進行更新并尋找最優(yōu)解的高效全局優(yōu)化(Efficient Global Optimization,EGO)算法后,大大提升了Kriging模型在工程優(yōu)化問題中的使用頻率及求解效率[3,6-7]。其中,SCHONLAU等[10]引入?yún)?shù)來調(diào)整EI準則的區(qū)域探索能力,提出了適用于包含黑箱約束問題的處理方法;ZHAN等[11]提出了通過設定閾值來調(diào)整和改進EI探索能力的并行EGO算法;SUPRAYITNO等[12]提出的可靠域Kriging代理進化優(yōu)化(Evolutionary Optimization using Reliable regional Kriging Surrogate, EORKS)算法,通過混合填充策略防止過早局部收斂。此外,基于代理模型化產(chǎn)生了諸多有代表性的優(yōu)化算法,例如基于置信下線準則(Lower Confidence Boundary, LCB)、可行性概率(Probability of Feasibility, PoF)準則的眾多算法[7,13-14]。隨著工程研究對象的復雜化,在計算機試驗建模技術的快速發(fā)展推動下,基于Kriging模型的代理優(yōu)化算法已成為解決高度非線性,且具備昂貴仿真黑盒模型特點問題的重要解決方案之一[6,13]。

    計算機試驗建模過程中發(fā)現(xiàn),EI填充準則新增試驗樣本點具有較好的全局探索能力,有效地改善了Kriging模型預測精度及建模效率[2,11]。在EI準則及PoF準則的基礎上產(chǎn)生了諸多填充準則,如張建俠等[15]提出了多目標策略和聚類分析的代理優(yōu)化算法,并將其應用于壓力容器的優(yōu)化設計中;JIAO等[16]在約束條件輔助下提出的約束完備EI準則,在探索當前可行區(qū)域的同時,還可開發(fā)不可行但有發(fā)展前途的區(qū)域;SUN等[17]通過Pareto前沿點集實現(xiàn)多點選取,改善新樣本點探索可能性的同時提出了并行多目標優(yōu)化算法;CINQUEGRANA等[18]通過給定逐步下降的權值來調(diào)整EI的探索能力,并將其應用于翼型優(yōu)化設計中。張鵬等[19]依據(jù)復相關系數(shù)R2實現(xiàn)Kriging模型的自適應設計,并將其應用于注塑機模板應力、變形測試及優(yōu)化。研究發(fā)現(xiàn),上述準則多集中于PoF準則,改變可靠性程度或者提升可行性區(qū)間的變化,鮮少有文獻以EI新增試驗點作為參考點,從而減少EI準則貪婪特性導致的陷入局部優(yōu)化問題[20]。因此,研究具有距離特點的權函數(shù)來建立新的期望填充準則顯然是一種大膽嘗試,但就目前研究來看尚有欠缺及不足[21-22]。

    本文對EI準則進行修改,提出兩種基于參考點的權期望調(diào)整準則,并將其與代理優(yōu)化算法相結(jié)合,綜合運用Kriging建模技術及進化算法實現(xiàn)對包含約束條件的優(yōu)化問題進行求解,提出了基于權期望填充準則的全局優(yōu)化算法。該算法通過具有距離特性的權函數(shù)實現(xiàn)對EI填充準則新增試驗點進行調(diào)整,使其在具備平衡全局探索和局部探索能力的同時,可以使新增樣本點具備距離特性,以提升建模效率及模型預測精度。最后,通過數(shù)值試驗算例及工程實例進行比較分析,驗證了所提權期望填充準則作用下的全局優(yōu)化算法的有效性。

    1 權期望加點準則作用下的全局優(yōu)化算法

    受制于工況、機械性能等因素,工程中的全局優(yōu)化問題一般均具有約束條件,其可用一般意義上的數(shù)學公式進行表示,具體形式如下:

    minf(x);

    s.t.

    gi(x)≤0,i=1,2,…,r,

    x∈D。

    (1)

    其中:f(x)是目標函數(shù);r為整數(shù),指共具有r個約束;gi(x)表示第i個約束條件;x∈D?Rd表示d維設計變量;D=[xl,xr]表示變量的矢量矩陣,xl,xr分別表示變量的左、右邊界。

    1.1 Kriging代理模型

    假定樣本集Ω中包含m個試驗點,定義樣本觀測變量矩陣為X=[x1,x2,…,xm]T,其對應的響應值矩陣為Y=[y1,y2,…,ym]T,Kriging代理模型可表示為:

    y(x)=μ+z(x)。

    (2)

    (3)

    Kriging模型中的相關函數(shù)通常選擇常用的核函數(shù),如高斯函數(shù)、樣條函數(shù)、Matérn函數(shù)等。具體可參考文獻[2-3]。

    1.2 權期望填充準則

    (4)

    EI(x)=EI1+EI2,

    為了平衡全局和局部探索能力,一般會對EI1和EI2分別賦予不同的權重,對兩者進行簡單的加權或采用熵權準則來依據(jù)EI1和EI2的重要程度確定權重[18],將EI函數(shù)修改為如下格式:

    EI(x)=w1EI1+w2EI2。

    式中w1,w2∈R,且w1+w2=1。

    xR=argmaxEI(x)。

    ZHAN等[21-22]在研究多點填充準則時提出了影響力函數(shù)(Influence Function,IF),來調(diào)整和近似EI改進函數(shù),稱為偽期望改進(Pseudo Expected Improvement, PEI)準則,其本質(zhì)思想是構造權函數(shù)(Weighted Function, WF)。以經(jīng)典的高斯相關核函數(shù)為例,權函數(shù)表達形式如下:

    權函數(shù)滿足權的一般性質(zhì)且具備距離特性,同時這也是具有空間距離特性相關函數(shù)對于Kriging模型特有的校正和平滑功能,但上述PEI準則僅被用于并行填充設計,未進行單個加點嘗試。為了避免EI準則過于貪婪而陷入局部最優(yōu),通過權函數(shù)對其進行調(diào)整的同時設定變量維度d的倒數(shù)作為冪指數(shù)來減緩其收斂速度?;谏鲜隹紤],將權函數(shù)分別應用于EI函數(shù)的不同位置,獲得兩種新的權期望填充準則,其表現(xiàn)形式如下:

    其中d為設計變量維數(shù)。

    1.3 可行性概率準則及約束處理方法

    EI填充準則是為解決無約束優(yōu)化問題設計的,對于具有復雜非線性黑箱約束的工程優(yōu)化問題,考慮新增試驗點在約束范圍內(nèi)的可行性概率準則非常適用[10]。

    現(xiàn)有文獻大多將PoF作為EI準則權重的約束應對方法應用在實際問題中,即以PoF為權重的約束期望改進(Constrained Expected Improvement, CEI)[2,10],可表示為:

    CEI(x)=E(1gi(x)≤0·I(x))=EI(x)·PoF(x)。

    同理,依據(jù)新的填充準則可獲得約束條件下的權期望改進準則:

    約束權期望準則偽實現(xiàn)代碼如下:

    步驟1計算滿足最大期望約束條件的參考點xR:

    步驟2計算新增試驗樣本點x:

    1.4 EGO算法及其實現(xiàn)流程

    EGO算法一般是指利用少量初始樣本構建初始Kriging模型,采用最大化期望改進填充準則選取增加新試驗點進行樣本填充設計,通過更新迭代過程不斷自適應修正模型,直至滿足終止條件獲取最優(yōu)解的算法實現(xiàn)過程。其算法流程圖1所示。

    上述流程實現(xiàn)具體偽代碼如下:

    步驟1最大最小拉丁抽取5d+1個觀測點,獲取樣本集(X,y)和約束矩陣G=[g1,g2,…,gr]T,設定終止條件,獲得初始最優(yōu)可行解ymin(x)。

    步驟2While終止條件不滿足do

    依據(jù)(X,y)建立目標函數(shù)初始Kriging模型

    fori=1 tor

    依據(jù)(X,G)構建約束初始Kriging模型

    end for

    步驟3判斷gi(x)≤0,樣本集是否有可行解。

    if無可行解,then采用可行性概率準則增加可行試驗點,

    else

    end

    步驟4更新樣本點x,y(x),gi(x),更新樣本集

    X=X∪x,y=y∪y(x),

    G=G∪[g1,g2,…,gr]T。

    步驟5更新當前最優(yōu)可行解ymin(x)

    end while

    2 測試算例與分析

    2.1 數(shù)值算例

    為了更好地說明新提加點準則的有效性,選取分別包含5個和10個變量的G4、G7數(shù)學測試函數(shù)[23],以及包含7個變量的減速器設計(Speed reduce design, SRD)案例[24]、11個變量的汽車側(cè)面碰撞(Car Side Impact, CSI)案例[25]、進行試驗驗證,并將結(jié)果與經(jīng)典EI準則進行對比,驗證所提算法的有效性、高效性和穩(wěn)健性,算例具體信息如表1所示。

    表1 算例相關信息

    2.2 優(yōu)化結(jié)果與分析

    所有測試均在MATLAB 2017a,Think Pad環(huán)境下運行;為保證結(jié)構穩(wěn)定性,采用最大最小拉丁設計抽取5d+1個樣本作為輸入樣本矩陣X,仿真獲取相應的響應矩陣Y及約束矩陣G。為方便分析,本文以仿真試驗的總次數(shù)T作為優(yōu)化算法的終止條件,當T=100時,算法停止,并返回最優(yōu)解。

    (1)算法收斂圖

    首先比較各加點準則下算法的收斂效果,為消除初始試驗設計隨機性對優(yōu)化結(jié)果的影響,用30組初始試驗設計下優(yōu)化結(jié)果的平均值隨迭代次數(shù)變化的曲線來評估各算法的尋優(yōu)能力和收斂性。3種加點方法的收斂效果如圖2所示。

    如圖2所示,迭代次數(shù)越少,結(jié)果越接近已知最優(yōu)值,則表明該準則下的算法收斂速度更快;最終解接近或小于已知最優(yōu)值則表明該準則下的算法精度更好。從圖2可以看出,4個算例在3種加點準則的作用下隨著迭代次數(shù)的增加均可收斂至近似最優(yōu)解,說明了兩種新權期望填充準則作用下的代理優(yōu)化算法的有效性。對于G7、CSI和SRD算例,在到達相同精度的情況下,CWEI準則相較于EI準則的迭代次數(shù)更少,說明CWEI準則作用下的優(yōu)化算法是加快收斂速度的有效方法,而對于CPEI準則,其收斂速度雖然相較于經(jīng)典EI準則和CWEI準則更慢,但其近似最優(yōu)解具有更好的精度。對于G4算例,雖然CPEI和CWEI準則收斂速度相較于EI準則更慢,但整體差距并不是很大,且能夠很好地收斂到近似最優(yōu)解。

    (2)算法性能表

    表2 不同優(yōu)化算法結(jié)果比較

    續(xù)表2

    對于G4和SRD算例,CPEI準則的最大值、平均值和標準差都優(yōu)于EI準則,說明CPEI準則作用下的優(yōu)化算法能夠更好地對最優(yōu)解區(qū)域進行探索,且該準則受初始試驗設計隨機性的影響最小;而CWEI準則在標準差方面優(yōu)于EI準則,雖然最大值和平均值沒有比EI準則更優(yōu),但是其結(jié)果接近已知最優(yōu)值,說明該準則是有效的。對于CSI算例,CWEI準則的結(jié)果更好,CPEI準則也具有較好的表現(xiàn)。對于G7算例,雖然EI準則的結(jié)果更優(yōu),但CWEI和CPEI準則的平均值相較于EI準則差距并不大,說明兩種準則作用下的優(yōu)化算法也具有較好的精度??傮w來說,新提兩種加點準則具有較好的有效性和尋優(yōu)穩(wěn)健性。

    (3)算法箱線圖

    對于每個算例中每種算法的30次結(jié)果,繪制結(jié)果的箱線圖,縱坐標是30次結(jié)果的近似最優(yōu)值,結(jié)果如圖3所示。

    對于G4算例,CWEI和CPEI準則的結(jié)果均比EI準則的結(jié)果更加集中,說明兩種準則受實驗隨機性影響較小,對于G7、CSI和SRD算例,CWEI和CPEI準則的結(jié)果均和EI準則的結(jié)果相近,且波動范圍不大。總體來說,CWEI準則及CPEI準則均可獲得較好的近似最優(yōu)解,且穩(wěn)健性較好。

    2.3 工程實例

    該實例是一個圓柱形容器設計,由圓柱體和半球體組成,結(jié)構如圖4所示。該實例有4個設計變量:Ts(圓柱體殼厚度)、Th(半球體厚度)、R(內(nèi)半徑)和L(容器的圓柱形部分的長度);優(yōu)化目的是使總成本降到最低;包括材料,成型和焊接的成本。工程實例的數(shù)學模型如下:

    minf(x)=0.622 4TsRL+1.778 1ThR2+

    s.t.

    g1(x)=-TS+0.019 3R≤0;g2(x)=-Th+0.009 54R≤0;

    g4(x)=L-240≤0,

    0.062 5≤Ts,Th≤0.062 5×99,10≤R,L≤200。

    運用本文提出的CWEI和CPEI準則,對此工程案例進行求解,并與EI準則以及其他文獻中的結(jié)果進行對比,如表3所示。對比結(jié)果發(fā)現(xiàn),近似解接近最優(yōu)值,說明了新提兩種加點準則的有效性,且兩種準則作用下的算法能夠很好地解決文獻中的實際工程問題。

    表3 結(jié)果對比表

    3 結(jié)束語

    算例試驗以及工程實例結(jié)果表明,基于Kriging代理優(yōu)化模型和兩種權期望填充準則的EGO算法可有效地處理和優(yōu)化復雜的黑箱過程。本文基于經(jīng)典EI準則和Kriging模型預測不確定性提出的兩種權期望填充準則具有良好的空間探索能力,所提兩種新EGO算法在解決單目標約束優(yōu)化問題時有效,且在算法的可行性、建模效率和解的穩(wěn)健性方面具有一定優(yōu)勢。具體結(jié)論如下:①3種填充準則作用下的EGO算法均可快速地收斂到近似最優(yōu)解,但從收斂速度來看,CWEI和CPEI準則略優(yōu)于EI準則;②3種填充準則作用下,CWEI準則受試驗隨機性影響最小,穩(wěn)健性最好,可有效減少試驗次數(shù);③相同條件下,CPEI和CWEI準則相比EI準則近似最優(yōu)解的標準差更小,解的質(zhì)量更高;④兩種新的加點準則,只需要較少的試驗就能求得優(yōu)化問題的最優(yōu)解,是提高工程優(yōu)化設計效率和質(zhì)量的有效途徑。

    試驗結(jié)果表明,通過權函數(shù)可有效調(diào)整EI準則的全局和局部探索能力,注意到3種不同的填充準則可實現(xiàn)不同位置的加點,故3種填充方法也可實現(xiàn)多點同時加點。未來,實現(xiàn)現(xiàn)有準則基礎上的一次多點加點是可以研究的領域;將該方法應用到解決多目標優(yōu)化問題也是未來研究的方向之一。

    猜你喜歡
    優(yōu)化模型
    一半模型
    超限高層建筑結(jié)構設計與優(yōu)化思考
    民用建筑防煙排煙設計優(yōu)化探討
    關于優(yōu)化消防安全告知承諾的一些思考
    一道優(yōu)化題的幾何解法
    由“形”啟“數(shù)”優(yōu)化運算——以2021年解析幾何高考題為例
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權M-估計的漸近分布
    3D打印中的模型分割與打包
    FLUKA幾何模型到CAD幾何模型轉(zhuǎn)換方法初步研究
    久久香蕉激情| 久久婷婷成人综合色麻豆| 亚洲美女黄片视频| 在线天堂中文资源库| 久久精品91无色码中文字幕| 日本三级黄在线观看| 午夜福利影视在线免费观看| 亚洲专区中文字幕在线| 成人免费观看视频高清| 日本a在线网址| av有码第一页| 一区二区日韩欧美中文字幕| 国产又色又爽无遮挡免费看| 波多野结衣高清无吗| 欧美精品亚洲一区二区| 黑人巨大精品欧美一区二区mp4| 国产亚洲av高清不卡| ponron亚洲| 天天影视国产精品| 一进一出好大好爽视频| 国产精品九九99| 国内毛片毛片毛片毛片毛片| 满18在线观看网站| 真人一进一出gif抽搐免费| 亚洲欧美激情综合另类| 国产精品久久电影中文字幕| 亚洲欧美日韩高清在线视频| 国产成人欧美| 免费高清视频大片| 欧美色视频一区免费| 久久热在线av| 国产高清视频在线播放一区| 国产精品二区激情视频| 如日韩欧美国产精品一区二区三区| 在线播放国产精品三级| 动漫黄色视频在线观看| 久久久水蜜桃国产精品网| 国产精品成人在线| 国产精品一区二区在线不卡| 久久婷婷成人综合色麻豆| 人人妻人人爽人人添夜夜欢视频| 久久精品亚洲av国产电影网| 国产在线精品亚洲第一网站| 少妇粗大呻吟视频| 亚洲av电影在线进入| 视频区图区小说| 久久国产亚洲av麻豆专区| 黄色 视频免费看| 亚洲久久久国产精品| 亚洲五月婷婷丁香| 热99国产精品久久久久久7| 久久精品91蜜桃| 国产91精品成人一区二区三区| 久久天躁狠狠躁夜夜2o2o| 亚洲avbb在线观看| 国产视频一区二区在线看| 日本 av在线| 国产精品久久久久久人妻精品电影| 精品久久蜜臀av无| 亚洲专区字幕在线| 国产主播在线观看一区二区| 青草久久国产| 在线观看免费高清a一片| 午夜两性在线视频| 欧美日韩亚洲高清精品| 日韩欧美一区视频在线观看| 女人高潮潮喷娇喘18禁视频| 波多野结衣高清无吗| 亚洲自拍偷在线| 午夜影院日韩av| 久久狼人影院| 99在线人妻在线中文字幕| 一边摸一边抽搐一进一小说| 99久久综合精品五月天人人| 色播在线永久视频| 午夜两性在线视频| 神马国产精品三级电影在线观看 | 9热在线视频观看99| 丝袜美足系列| 搡老熟女国产l中国老女人| 一边摸一边抽搐一进一出视频| 国产一区二区三区视频了| 亚洲欧美日韩无卡精品| 欧美激情 高清一区二区三区| 午夜福利在线免费观看网站| 视频在线观看一区二区三区| 人人妻人人澡人人看| svipshipincom国产片| 在线国产一区二区在线| 侵犯人妻中文字幕一二三四区| 亚洲欧美日韩另类电影网站| 免费在线观看完整版高清| 国产精华一区二区三区| 丝袜美足系列| 久久亚洲精品不卡| 69精品国产乱码久久久| 无限看片的www在线观看| 国产一区在线观看成人免费| 日韩人妻精品一区2区三区| 伦理电影免费视频| 日韩免费av在线播放| 国产免费现黄频在线看| 精品久久蜜臀av无| 丰满饥渴人妻一区二区三| 国产又色又爽无遮挡免费看| 成人亚洲精品一区在线观看| 两性夫妻黄色片| svipshipincom国产片| 黄色成人免费大全| 成人免费观看视频高清| 老司机午夜福利在线观看视频| 欧美日韩av久久| 亚洲成人国产一区在线观看| 中文字幕人妻熟女乱码| 中文字幕精品免费在线观看视频| 国产91精品成人一区二区三区| 国产av又大| 啦啦啦在线免费观看视频4| 日韩欧美国产一区二区入口| 亚洲国产欧美一区二区综合| 久久青草综合色| 久久久久亚洲av毛片大全| 亚洲欧美一区二区三区黑人| 91老司机精品| 亚洲成人免费av在线播放| 性色av乱码一区二区三区2| 精品久久久久久成人av| 久久久久久人人人人人| 久久久久久久精品吃奶| 免费观看精品视频网站| 亚洲精品一区av在线观看| 日日干狠狠操夜夜爽| svipshipincom国产片| 热re99久久精品国产66热6| 美女高潮到喷水免费观看| 免费在线观看亚洲国产| 在线观看66精品国产| 国产精品一区二区三区四区久久 | 精品乱码久久久久久99久播| 日本一区二区免费在线视频| 中文字幕精品免费在线观看视频| 午夜福利一区二区在线看| 一夜夜www| 亚洲九九香蕉| 美女国产高潮福利片在线看| 在线十欧美十亚洲十日本专区| 黑人巨大精品欧美一区二区mp4| 少妇粗大呻吟视频| 久热这里只有精品99| 免费女性裸体啪啪无遮挡网站| 一区在线观看完整版| 国产高清videossex| 高清欧美精品videossex| 神马国产精品三级电影在线观看 | 俄罗斯特黄特色一大片| 国产国语露脸激情在线看| 50天的宝宝边吃奶边哭怎么回事| 国产精品一区二区精品视频观看| www.自偷自拍.com| 日日干狠狠操夜夜爽| 精品第一国产精品| 在线视频色国产色| 两性夫妻黄色片| 在线观看舔阴道视频| 99re在线观看精品视频| 三级毛片av免费| www日本在线高清视频| 久久天堂一区二区三区四区| 夜夜夜夜夜久久久久| 黄色成人免费大全| 俄罗斯特黄特色一大片| 精品久久久久久久久久免费视频 | 欧美午夜高清在线| 日韩高清综合在线| 9191精品国产免费久久| 精品一区二区三卡| 国产av又大| 国产亚洲精品一区二区www| 国产主播在线观看一区二区| 热99re8久久精品国产| 亚洲国产精品合色在线| 国产精品日韩av在线免费观看 | 9191精品国产免费久久| 深夜精品福利| 熟女少妇亚洲综合色aaa.| 黄色 视频免费看| 黄频高清免费视频| 国产精品野战在线观看 | 91老司机精品| 精品一品国产午夜福利视频| 人人妻人人澡人人看| 激情在线观看视频在线高清| 成人亚洲精品一区在线观看| 黑人欧美特级aaaaaa片| 国产精品香港三级国产av潘金莲| 三级毛片av免费| 亚洲色图av天堂| 国产欧美日韩精品亚洲av| 国产成人欧美在线观看| 18禁观看日本| 在线观看免费日韩欧美大片| 老熟妇乱子伦视频在线观看| 国产精品自产拍在线观看55亚洲| 精品一区二区三区av网在线观看| 国产又色又爽无遮挡免费看| 日本撒尿小便嘘嘘汇集6| 侵犯人妻中文字幕一二三四区| 亚洲av日韩精品久久久久久密| √禁漫天堂资源中文www| 成人av一区二区三区在线看| 国产精品久久视频播放| 亚洲一区二区三区不卡视频| 99在线人妻在线中文字幕| av网站在线播放免费| 亚洲国产欧美日韩在线播放| 1024视频免费在线观看| 成人永久免费在线观看视频| 欧美激情极品国产一区二区三区| 丰满人妻熟妇乱又伦精品不卡| 97碰自拍视频| 两性夫妻黄色片| 久久精品亚洲熟妇少妇任你| 午夜福利免费观看在线| 精品久久久久久久久久免费视频 | 精品国产一区二区久久| 身体一侧抽搐| 亚洲一区二区三区欧美精品| 老司机深夜福利视频在线观看| 国产成年人精品一区二区 | 99热只有精品国产| 法律面前人人平等表现在哪些方面| 老熟妇仑乱视频hdxx| av超薄肉色丝袜交足视频| 可以在线观看毛片的网站| avwww免费| 午夜老司机福利片| 麻豆久久精品国产亚洲av | 丝袜在线中文字幕| 亚洲人成电影观看| 欧美日韩精品网址| 亚洲人成77777在线视频| 中文字幕最新亚洲高清| 欧美日本亚洲视频在线播放| 人人妻人人澡人人看| 午夜两性在线视频| 色精品久久人妻99蜜桃| 不卡av一区二区三区| 欧美日本中文国产一区发布| 精品第一国产精品| 他把我摸到了高潮在线观看| 亚洲avbb在线观看| 亚洲第一av免费看| 极品人妻少妇av视频| 国产亚洲av高清不卡| 国内毛片毛片毛片毛片毛片| 色综合婷婷激情| 国产91精品成人一区二区三区| 狂野欧美激情性xxxx| 老熟妇仑乱视频hdxx| 这个男人来自地球电影免费观看| 欧美不卡视频在线免费观看 | 女性被躁到高潮视频| 国产三级黄色录像| 黄色a级毛片大全视频| 巨乳人妻的诱惑在线观看| 在线观看免费午夜福利视频| 妹子高潮喷水视频| 久久精品成人免费网站| 曰老女人黄片| 一夜夜www| 一级毛片高清免费大全| 免费高清视频大片| 精品一区二区三区av网在线观看| 亚洲精品久久成人aⅴ小说| 免费不卡黄色视频| 97碰自拍视频| bbb黄色大片| 亚洲熟女毛片儿| а√天堂www在线а√下载| 真人做人爱边吃奶动态| 一二三四社区在线视频社区8| 久久久国产一区二区| 国产深夜福利视频在线观看| av福利片在线| 女人被躁到高潮嗷嗷叫费观| 两性夫妻黄色片| 黄色成人免费大全| 一级片免费观看大全| 搡老熟女国产l中国老女人| 亚洲熟妇熟女久久| 91老司机精品| 88av欧美| xxx96com| 他把我摸到了高潮在线观看| 国产精品永久免费网站| 国产高清激情床上av| 狂野欧美激情性xxxx| 97碰自拍视频| 国产成人一区二区三区免费视频网站| 精品电影一区二区在线| 操出白浆在线播放| 老熟妇仑乱视频hdxx| 一区二区三区精品91| 亚洲av五月六月丁香网| 亚洲黑人精品在线| 亚洲色图av天堂| 一夜夜www| 在线观看舔阴道视频| 国产精品二区激情视频| 性欧美人与动物交配| 国产精品国产高清国产av| 在线观看一区二区三区| 天堂中文最新版在线下载| 日韩有码中文字幕| 热99re8久久精品国产| 午夜成年电影在线免费观看| 欧美精品亚洲一区二区| 日本精品一区二区三区蜜桃| 99久久综合精品五月天人人| 国产在线观看jvid| 91麻豆av在线| 精品国产一区二区三区四区第35| 少妇的丰满在线观看| 久久久久亚洲av毛片大全| 欧美日韩一级在线毛片| 亚洲精品国产区一区二| 久久午夜亚洲精品久久| 国产一区二区三区视频了| a在线观看视频网站| 97人妻天天添夜夜摸| 狠狠狠狠99中文字幕| 午夜福利影视在线免费观看| 日本免费a在线| 丝袜美足系列| 亚洲性夜色夜夜综合| 亚洲中文日韩欧美视频| 国产免费现黄频在线看| 亚洲午夜理论影院| 日韩欧美一区二区三区在线观看| 男女高潮啪啪啪动态图| 国产欧美日韩一区二区三| www.精华液| 美女大奶头视频| 操美女的视频在线观看| 欧美一区二区精品小视频在线| 午夜精品在线福利| 亚洲国产欧美日韩在线播放| 操出白浆在线播放| 99久久人妻综合| 在线观看免费午夜福利视频| 欧美+亚洲+日韩+国产| 日本一区二区免费在线视频| 黑人巨大精品欧美一区二区蜜桃| 97超级碰碰碰精品色视频在线观看| 国产深夜福利视频在线观看| e午夜精品久久久久久久| 好男人电影高清在线观看| 欧美日韩福利视频一区二区| 国产深夜福利视频在线观看| 夜夜夜夜夜久久久久| 日本五十路高清| 久久这里只有精品19| 日韩免费av在线播放| av福利片在线| 久久久久亚洲av毛片大全| 国产精品秋霞免费鲁丝片| 亚洲专区中文字幕在线| 性少妇av在线| 91成人精品电影| 亚洲免费av在线视频| 亚洲精品在线观看二区| 国产成人欧美在线观看| 在线观看一区二区三区| 久久久久精品国产欧美久久久| 日韩欧美一区二区三区在线观看| 精品第一国产精品| 99在线人妻在线中文字幕| 亚洲色图综合在线观看| 国产免费现黄频在线看| 国产伦一二天堂av在线观看| 亚洲欧美一区二区三区黑人| 欧美精品亚洲一区二区| 黄片大片在线免费观看| 麻豆久久精品国产亚洲av | 悠悠久久av| 亚洲一卡2卡3卡4卡5卡精品中文| 999精品在线视频| 久久香蕉激情| 乱人伦中国视频| 日本wwww免费看| av天堂在线播放| 操美女的视频在线观看| 日韩成人在线观看一区二区三区| 熟女少妇亚洲综合色aaa.| 久热这里只有精品99| 国产av一区在线观看免费| 男人舔女人下体高潮全视频| 国产高清videossex| 日韩欧美三级三区| 国产精品电影一区二区三区| 久久久精品欧美日韩精品| 搡老熟女国产l中国老女人| 91av网站免费观看| 欧美不卡视频在线免费观看 | 女人被躁到高潮嗷嗷叫费观| 极品人妻少妇av视频| 91成人精品电影| 成年女人毛片免费观看观看9| 国产一区在线观看成人免费| 国产片内射在线| 欧美黄色片欧美黄色片| 可以在线观看毛片的网站| videosex国产| 无限看片的www在线观看| 欧美成人性av电影在线观看| 国产精品98久久久久久宅男小说| 国产精华一区二区三区| 天天躁夜夜躁狠狠躁躁| 性欧美人与动物交配| 99re在线观看精品视频| 91大片在线观看| 88av欧美| 亚洲av片天天在线观看| 满18在线观看网站| 国产亚洲精品久久久久久毛片| 国产精华一区二区三区| 嫩草影视91久久| 欧美黄色淫秽网站| 午夜福利欧美成人| 成人18禁高潮啪啪吃奶动态图| 亚洲av美国av| 美女福利国产在线| 性色av乱码一区二区三区2| 欧美日韩av久久| 超色免费av| 亚洲国产欧美日韩在线播放| av片东京热男人的天堂| 在线观看免费视频网站a站| 国产aⅴ精品一区二区三区波| 日本精品一区二区三区蜜桃| 亚洲av成人一区二区三| 男女下面插进去视频免费观看| 又紧又爽又黄一区二区| 波多野结衣一区麻豆| 日韩欧美国产一区二区入口| 大码成人一级视频| 亚洲精品一二三| 亚洲 欧美一区二区三区| 天天添夜夜摸| 一边摸一边抽搐一进一出视频| 国产区一区二久久| 久久狼人影院| 淫妇啪啪啪对白视频| 国产精品九九99| av电影中文网址| 日韩av在线大香蕉| 999精品在线视频| 很黄的视频免费| 天天影视国产精品| 人妻久久中文字幕网| 免费av毛片视频| av国产精品久久久久影院| 黑丝袜美女国产一区| 色老头精品视频在线观看| 18美女黄网站色大片免费观看| 动漫黄色视频在线观看| 日韩一卡2卡3卡4卡2021年| 国产精品久久久av美女十八| 高清在线国产一区| 国产一区二区在线av高清观看| netflix在线观看网站| 男女下面插进去视频免费观看| 亚洲成人免费电影在线观看| 一进一出抽搐动态| 美女午夜性视频免费| 久久久久精品国产欧美久久久| 国产一区二区三区综合在线观看| av电影中文网址| 色精品久久人妻99蜜桃| 午夜精品国产一区二区电影| 国产成人欧美在线观看| 亚洲第一青青草原| 一区二区三区国产精品乱码| 久久99一区二区三区| 黄网站色视频无遮挡免费观看| 视频区图区小说| 国产一区在线观看成人免费| 夜夜躁狠狠躁天天躁| 亚洲性夜色夜夜综合| 午夜激情av网站| 麻豆国产av国片精品| 午夜免费成人在线视频| 久久伊人香网站| 在线天堂中文资源库| 国产精品久久久久成人av| 日韩高清综合在线| 亚洲伊人色综图| √禁漫天堂资源中文www| 欧美日本中文国产一区发布| 亚洲va日本ⅴa欧美va伊人久久| 日韩高清综合在线| 人人妻,人人澡人人爽秒播| 免费久久久久久久精品成人欧美视频| 91精品三级在线观看| 搡老熟女国产l中国老女人| 中亚洲国语对白在线视频| 99久久久亚洲精品蜜臀av| 首页视频小说图片口味搜索| 亚洲视频免费观看视频| 男男h啪啪无遮挡| 精品福利永久在线观看| 欧美亚洲日本最大视频资源| 大码成人一级视频| 91精品三级在线观看| 无限看片的www在线观看| 午夜免费激情av| 人人妻,人人澡人人爽秒播| svipshipincom国产片| 欧美日本中文国产一区发布| 欧美国产精品va在线观看不卡| a在线观看视频网站| 99国产极品粉嫩在线观看| 色尼玛亚洲综合影院| 女性生殖器流出的白浆| a级毛片黄视频| 午夜精品久久久久久毛片777| 日韩成人在线观看一区二区三区| 国产av一区在线观看免费| 亚洲第一青青草原| 我的亚洲天堂| 久久中文看片网| 中文字幕最新亚洲高清| 一进一出好大好爽视频| 午夜成年电影在线免费观看| 精品一品国产午夜福利视频| 日韩欧美在线二视频| 一区二区日韩欧美中文字幕| 麻豆一二三区av精品| 精品福利永久在线观看| 天天躁夜夜躁狠狠躁躁| 黑人猛操日本美女一级片| 国产一区二区激情短视频| 欧美日韩av久久| 无人区码免费观看不卡| 一区二区三区激情视频| aaaaa片日本免费| 1024香蕉在线观看| 成人永久免费在线观看视频| 成人手机av| 国产三级在线视频| 久久人人爽av亚洲精品天堂| 日本免费a在线| 国产一区二区三区视频了| 欧美一级毛片孕妇| 啦啦啦 在线观看视频| 手机成人av网站| 欧美日韩一级在线毛片| 男女做爰动态图高潮gif福利片 | 色婷婷久久久亚洲欧美| 波多野结衣一区麻豆| 久久狼人影院| 中文字幕高清在线视频| 又大又爽又粗| 涩涩av久久男人的天堂| 欧美久久黑人一区二区| 国产亚洲精品久久久久5区| 老汉色∧v一级毛片| 成人亚洲精品av一区二区 | 亚洲一区中文字幕在线| 中文字幕人妻熟女乱码| 黄色女人牲交| 好看av亚洲va欧美ⅴa在| 成年人黄色毛片网站| 国产又爽黄色视频| 丁香六月欧美| 在线观看午夜福利视频| 日韩成人在线观看一区二区三区| 一级毛片女人18水好多| 久久久久久久精品吃奶| 99久久人妻综合| 精品第一国产精品| 久久精品国产亚洲av高清一级| 纯流量卡能插随身wifi吗| 多毛熟女@视频| 国产精品爽爽va在线观看网站 | 精品第一国产精品| 狠狠狠狠99中文字幕| 男人舔女人的私密视频| 国产人伦9x9x在线观看| 性色av乱码一区二区三区2| 老熟妇仑乱视频hdxx| 高清在线国产一区| 少妇裸体淫交视频免费看高清 | 免费看十八禁软件| 精品久久久久久,| 国产激情欧美一区二区| 久久精品91无色码中文字幕| 久久久国产欧美日韩av| 可以在线观看毛片的网站| 国产黄色免费在线视频| 久久欧美精品欧美久久欧美| 亚洲在线自拍视频| 人成视频在线观看免费观看| 午夜免费成人在线视频| 午夜免费鲁丝| 日韩人妻精品一区2区三区| 亚洲美女黄片视频| 亚洲精品av麻豆狂野| 一边摸一边抽搐一进一小说| 亚洲美女黄片视频| 亚洲精品av麻豆狂野| 十八禁人妻一区二区| 欧美在线黄色| 亚洲午夜理论影院| 99在线视频只有这里精品首页|