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

    基于SA-QPSO算法的多層制造單元內(nèi)部布局方法

    2024-12-04 00:00:00郭軼男姜雪松張宇晨劉森王婧
    關(guān)鍵詞:模擬退火算法

    摘 要:針對(duì)多層制造單元內(nèi)部的設(shè)備布局優(yōu)化問(wèn)題,本文建立考慮單元尺寸、物料搬運(yùn)量、損失時(shí)間以及單元穩(wěn)定性的多目標(biāo)優(yōu)化數(shù)學(xué)模型。為更快速、高效地求解該問(wèn)題,使用模擬退火算法(Simulated Annealing,SA)確定單元內(nèi)設(shè)備所在平面以及層面,使用量子粒子群算法(Quantum Particle Swarm Optimization,QPSO)確定設(shè)備具體坐標(biāo)值和所在高度。以某汽車(chē)零件加工車(chē)間為實(shí)例,運(yùn)用SA-QPSO算法生成直線形、“U”形和環(huán)形3種最優(yōu)空間布局方案,驗(yàn)證了SA-QPSO算法在多層制造單元內(nèi)部布局方法設(shè)計(jì)方面的可行性。

    關(guān)鍵詞:多層制造單元;布局優(yōu)化;模擬退火算法;量子粒子群算法

    中圖分類(lèi)號(hào):TU 318 文獻(xiàn)標(biāo)志碼:A

    設(shè)施布局(Facility Layout Problem,F(xiàn)LP)是智能工廠建設(shè)的重要環(huán)節(jié),科學(xué)合理的設(shè)備布局有重要的意義[1]。隨著使用智能化物流的智能工廠興起,多層制造單元的布局方式越來(lái)越多,許多學(xué)者對(duì)制造單元布局問(wèn)題進(jìn)行深入研究,例如王沙沙等[2]針對(duì)雙向環(huán)形過(guò)道布置問(wèn)題建模并使用混合鯨魚(yú)算法有效求解。丁祥海等[3]以“U”形生產(chǎn)單元為研究對(duì)象,建立了單元生產(chǎn)可重構(gòu)設(shè)施布局決策集成模型,并設(shè)計(jì)了決策問(wèn)題的遺傳算法。目前國(guó)內(nèi)外主要研究單層制造單元,多層制造單元相關(guān)研究較少,其中,單目標(biāo)優(yōu)化模型或雙目標(biāo)優(yōu)化模型的搬運(yùn)成本最低、單元利用率最高?;诖耍疚目紤]了眾多制造車(chē)間的實(shí)際需求,構(gòu)建多目標(biāo)優(yōu)化數(shù)學(xué)模型,研究多層制造單元內(nèi)部設(shè)備布局方法。

    1 多層制造單元車(chē)間內(nèi)部布局問(wèn)題描述以及數(shù)學(xué)模型構(gòu)建

    1.1 問(wèn)題描述及假設(shè)

    多層制造單元內(nèi)部設(shè)備形狀以矩形為主,常見(jiàn)分類(lèi)有直線形制造單元、“U”形制造單元和環(huán)形制造單元。通過(guò)構(gòu)建X、Y、Z坐標(biāo)系,可立體表述制造單元,如圖1所示。其中,Z1為XOY平面,Z2為與XOY平面平行的平面,Z3為YOZ平面,Z4為與YOZ平面平行的平面,Z1、Z2、Z3和Z4均穿過(guò)設(shè)備中心,同層設(shè)備中心在同一條水平線上。

    將單元布局平面抽象化。如圖2所示。其中,直線形制造單元內(nèi)設(shè)備布局僅在平面Z1和Z2內(nèi);“U”形制造單元內(nèi)設(shè)備布局僅在平面Z1、Z2和Z3內(nèi);環(huán)形制造單元內(nèi)設(shè)備布置在所有平面內(nèi)。

    通過(guò)構(gòu)建相同的目標(biāo)函數(shù),采用約束條件選擇不同的布局形式,求解不同的設(shè)備排列組合中最優(yōu)的順序布置,獲得最佳布局。

    1.2 目標(biāo)函數(shù)

    本文以單元體積最小、物流搬運(yùn)強(qiáng)度最低、制造損失時(shí)間最短和單元穩(wěn)定性最高為目標(biāo)構(gòu)建多目標(biāo)數(shù)學(xué)模型。

    1.2.1 單元體積目標(biāo)

    多層制造單元占用空間為單元體積,主要受長(zhǎng)、寬和高3個(gè)要素影響,當(dāng)3個(gè)要素最小時(shí),單元體積就最小,目標(biāo)函數(shù)如公式(1)所示。

    minV=LWH (1)

    式中:V為單元體積;L為單元長(zhǎng)度,在單元內(nèi)設(shè)備的中心坐標(biāo)的x值與設(shè)備投影到X軸上的一半之和或之差組合的集合中,最大值和最小值之差即單元長(zhǎng)度目標(biāo);W為單元寬度,在單元內(nèi)設(shè)備的中心坐標(biāo)的z值與設(shè)施投影到Z軸上的一半之和或之差組合的集合中,最大值和最小值之差即單元寬度;H為單元高度目標(biāo),可將單元內(nèi)每個(gè)設(shè)備的中心在Y軸上的坐標(biāo)值y與該設(shè)備投影到Y(jié)軸上尺寸的一半相加,求得的結(jié)果組成的集合為單元高度集,集合中最大值就是制造單元所需的最低高度,即單元高度目標(biāo)。

    1.2.2 物料搬運(yùn)目標(biāo)

    本文搬運(yùn)距離采用與物流搬運(yùn)實(shí)際距離更接近的曼哈頓距離,如公式(2)所示。

    dij=|xi-xj|+|yi-yj|+|zi-zj| (2)

    式中:dij為設(shè)備i到設(shè)備j的曼哈頓距離,設(shè)備i中心坐標(biāo)為(xi,yi,zi),設(shè)備j中心坐標(biāo)為(xj,yj,zj)。

    第p號(hào)產(chǎn)品從設(shè)備i到設(shè)備j的物料搬運(yùn)批量為ceil[],搬運(yùn)目標(biāo)如公式(3)所示。

    (3)

    式中:D為物料搬運(yùn)強(qiáng)度;M為設(shè)備總數(shù)量;P為產(chǎn)品總數(shù)量;i、j為設(shè)備編號(hào);p為加工產(chǎn)品種類(lèi)編號(hào);ceil為取整;Vp為第p號(hào)產(chǎn)品待加工總數(shù)量;Bp為第p號(hào)產(chǎn)品一次搬運(yùn)數(shù)量。

    1.2.3 損失時(shí)間目標(biāo)

    在生產(chǎn)線平衡問(wèn)題中,損失時(shí)間為衡量制造單元平衡的重要指標(biāo)之一,因此設(shè)損失時(shí)間目標(biāo),如公式(4)所示。

    (4)

    式中:T為損失時(shí)間目標(biāo);np為設(shè)備總數(shù);qi為第i個(gè)設(shè)備加工以及運(yùn)輸?shù)降趇+1個(gè)設(shè)備的時(shí)間。

    1.2.4 穩(wěn)定性目標(biāo)

    在多層制造單元設(shè)備布局中,需要考慮設(shè)備的穩(wěn)定性,穩(wěn)定性目標(biāo)B越小就越穩(wěn)定,目標(biāo)函數(shù)如公式(5)所示。

    (5)

    式中:yi為設(shè)備i的中心高度;wi為設(shè)備i的質(zhì)量。

    1.3 約束條件

    在設(shè)定具體約束前,須設(shè)定的基本約束如下。設(shè)備i的中心坐標(biāo)(xi,yi,zi)為優(yōu)化問(wèn)題的決策變量;(li,wi,hi)為設(shè)備對(duì)應(yīng)的長(zhǎng)、寬、高;設(shè)備單元集均在X、Y、Z軸正方向上,不考慮X、Y、Z軸負(fù)方向,任意設(shè)備不能超過(guò)單元的邊界。

    1.3.1 邊界約束

    邊界約束為設(shè)備對(duì)應(yīng)長(zhǎng)度li的一半與中心點(diǎn)xi坐標(biāo)之和≤單元Z軸方向的長(zhǎng)度邊界約束Lg,設(shè)備所對(duì)應(yīng)寬度wi的一半與中心點(diǎn)yi坐標(biāo)之和≤單元X軸方向的寬度邊界約束Wg,設(shè)備對(duì)應(yīng)高度hi與中心點(diǎn)zi坐標(biāo)之和≤單元Y軸方向的高度約束Hg。

    1.3.2 設(shè)備位置約束

    設(shè)備中心必須位于平面Z1、Z2、Z3和Z4上,zi為非負(fù)區(qū)域。假設(shè)PSi∈I為設(shè)備中心Zi面的歸屬,當(dāng)PSi=1時(shí),zi=0;當(dāng)PSi=2時(shí),zi=Z2_value;當(dāng)PSi=3時(shí),xi=0;當(dāng)PSi=4時(shí),xi=Z4_value。其中,Z2_value為Z2面的z坐標(biāo),Z4_value為Z4面的x坐標(biāo)。

    1.3.3 單元集形狀約束

    假設(shè)Type_choose∈{1,2,3}為單元集形狀,當(dāng)Type_choose=1時(shí),直線形僅選擇平面Z1、Z3;當(dāng)Type_choose=2時(shí),“U”形僅選擇平面Z1、Z2以及Z3;當(dāng)Type_choose=3時(shí),環(huán)形選擇平面Z1、Z2、Z3以及Z4。

    1.3.4 高度層級(jí)約束

    假設(shè)Yk為第k個(gè)高度層級(jí),Yi=0只能從Yk(k=1,2,…,r)中取值,即yi∈{Yk:k=1,2,…,r}。

    2 基于SA-QPSO算法的目標(biāo)函數(shù)求解

    2.1 設(shè)備布局方案的編碼與生成

    在優(yōu)化問(wèn)題中,粒子群優(yōu)化算法(Particle Swarm Optimization,PSO)和模擬退火算法(Simulated Annealing,SA)都是常見(jiàn)且有效的算法,量子粒子群算法(Quantum Particle Swarm Optimization,QPSO)是基于PSO算法提出的一種改進(jìn)算法[4]。本文采用了一種改進(jìn)的混合算法,稱為量子粒子群優(yōu)化算法(SA-QPSO),并將其應(yīng)用于求解單元內(nèi)部設(shè)備布局最優(yōu)化問(wèn)題。流程如圖3所示。

    2.1.1 初步布局方案

    設(shè)備j隨機(jī)生成數(shù)i和數(shù)k,i∈{1,2,3,4},k∈{1,2,…,r},表示將設(shè)備j分配至Zi面的第k層中,由此可得Aik,即Zi面第k層高度的設(shè)備編號(hào)升序向量。設(shè)Aik共有mik個(gè)元素,如果mik>0,則生成一個(gè)mik維的各分量取值為(0,1)的隨機(jī)向量Vik,根據(jù)Vik各分量的大小關(guān)系對(duì)Aik進(jìn)行排序,結(jié)合Aik和Vik可得設(shè)備在zi面第k層上沿x正向或z正向的有序排布,得到初步布局方案。

    2.1.2 Z4的x值生成

    僅考慮x方向的干涉,從平面Z1、Z2找出x坐標(biāo)最大的設(shè)備,使其不與Z4內(nèi)最長(zhǎng)的設(shè)備發(fā)生干涉,又因?yàn)閾Q行約束,所以Z4≥Lz,Z4的下界lb4和上界ub4如公式(6)、公式(7)所示。

    (6)

    式中:Lz為換行約束;Zik,j為Zi的第k層中的第j個(gè)設(shè)備的前向中心間隔;lik,end為平面Zi第k層高度最后1個(gè)設(shè)備的長(zhǎng)度;l4k,j為平面Z4第k層高度設(shè)備j的長(zhǎng)度;G1為給定常數(shù)值。

    (7)

    式中:Lg為單元Z軸方向的長(zhǎng)度邊界約束。

    由于此時(shí)最大余量為正,因此ub4gt;lb4,如公式(8)所示。

    Z4=lb4+(ub4-lb4)·Z4_ratio " " "(8)

    式中:Z4_value為隨機(jī)數(shù)且Z4_ratio∈(0,1)。

    2.1.3 Z2的Z值生成

    Z2的下界lb2和上界ub2如公式(9)、公式(10)所示。

    (9)

    式中:Wz為換行約束;zik,j為Zi的第k層中的第j個(gè)設(shè)備的前向中心間隔;wik,end為平面Zi第k層高度最后1個(gè)設(shè)備的寬度;w2k,j為平面Z2第k層高度設(shè)備j的寬度;G3為給定常數(shù)值。

    (10)

    式中:Wg為單元X軸方向的寬度邊界約束。

    設(shè)隨機(jī)數(shù)Z2-ratio∈(0,1),如公式(11)所示。

    Z2=lb2+(ub2-lb2)·Z2-ratio (11)

    至此,分配至Z1、Z2、Z3和Z4的設(shè)備都已定x坐標(biāo)和z坐標(biāo):平面Z1的z=0;平面Z2的z值已生成;平面Z3的x=0;平面Z4的x值已生成。

    2.1.4 高度層級(jí)Yk的生成

    假設(shè)ηk為第k層設(shè)備的最大高度,那么以層級(jí)r=3為例進(jìn)行分析,如圖4所示。

    對(duì)各層的最高設(shè)備按層級(jí)順序編號(hào),計(jì)算步驟如下。

    步驟一:計(jì)算設(shè)備最大余量,如公式(12)所示。

    (12)

    式中:Marginymax為設(shè)備上側(cè)Y軸坐標(biāo)的最大邊緣值即最大余量;ηk為第k層高度設(shè)備的最大高度;r為高度層;Hg為設(shè)備對(duì)應(yīng)高度hi與中心點(diǎn)zi坐標(biāo)之和≤單元Y軸方向的高度約束;G2為給定常數(shù)值。如果最大余量≤0,則判定無(wú)解Marginymax為最大余量。

    步驟二:生成r-1維(0,1)隨機(jī)向量,如公式(13)、公式(14)所示。

    ra=[ra1,ra2,…,rar-1] (13)

    (14)

    式中:φk為第k個(gè)最高設(shè)備的前向中心間隔,φ1取0,k≥2;βy為占比向量;yratio為生成隨機(jī)數(shù)且yratio∈(0,1)。

    步驟三:生成高度層級(jí),如公式(15)所示。

    (15)

    至此,生成布局方案且各約束均得到了滿足,可得到可行布局方案或判定無(wú)解。該問(wèn)題外層為帶有簡(jiǎn)單約束的組合優(yōu)化問(wèn)題,使用帶有罰函數(shù)的SA算法進(jìn)行優(yōu)化,使各高度層級(jí)均有設(shè)備;內(nèi)層為無(wú)約束優(yōu)化問(wèn)題,使用QPSO算法進(jìn)行求解。

    2.2 SA算法

    2.2.1 解的編碼與初始化

    假設(shè)共有n個(gè)設(shè)備,以2×n的整數(shù)矩陣來(lái)表示設(shè)備組的Zi面屬性和高度層級(jí)屬性,每列對(duì)應(yīng)1個(gè)設(shè)備,第一行為Zi面屬性,第二行為高度層級(jí)屬性。根據(jù)矩陣進(jìn)行。

    編碼下界矩陣如公式(16)所示。

    (16)

    直線形上界矩陣如公式(17)所示。

    (17)

    “U”形上界矩陣如公式(18)所示。

    (18)

    環(huán)形上界矩陣如公式(19)所示。

    (19)

    2.2.2 適應(yīng)度函數(shù)

    對(duì)該嵌套形式的優(yōu)化問(wèn)題來(lái)說(shuō),SA算法的適應(yīng)度為內(nèi)層QPSO算法的最優(yōu)適應(yīng)度值。在最大化目標(biāo)中,令適應(yīng)度為-M;在最小化目標(biāo)中,令適應(yīng)度為M。其中,M為充分大的正數(shù)。

    2.3 QPSO算法

    在QPSO算法中,通過(guò)粒子可以搜索整個(gè)可行解空間,并將粒子群算法中的速度和位置更新公式進(jìn)行了替換[5]。同時(shí)該算法以概率分布機(jī)制替代原有機(jī)制,提高了全局搜索的能力,也提高了算法收斂的精度和速度[6]。QPSO 算法的具體步驟如下。

    步驟一:初始化粒子位置,并設(shè)置初始個(gè)體最優(yōu)值。

    步驟二:計(jì)算各粒子的適應(yīng)度值。

    步驟三:比較粒子適應(yīng)度值的大小,更新個(gè)體最優(yōu)和全局最優(yōu)。

    步驟四:更新粒子位置,如公式(20)所示。

    (20)

    式中:為最新粒子位置信息;pit,j為局部吸引子;Cjt為粒子當(dāng)前個(gè)體最優(yōu)的均值;α為控制參數(shù);u為隨機(jī)變量,由此得到新位置。

    步驟五:按照一定的變異概率對(duì)個(gè)體最優(yōu)進(jìn)行高斯擾動(dòng),即加入服從高斯分布的隨機(jī)數(shù)來(lái)提升算法多樣性。

    步驟六:重復(fù)步驟二~步驟四直到滿足停止條件。

    3 應(yīng)用實(shí)例

    3.1 實(shí)例描述

    采用上述算法對(duì)某汽車(chē)零件單元生產(chǎn)車(chē)間進(jìn)行布局優(yōu)化,本文選擇該車(chē)間的15臺(tái)設(shè)備為研究對(duì)象,共生產(chǎn)12種產(chǎn)品,設(shè)備信息包括設(shè)備編號(hào)、設(shè)備尺寸(長(zhǎng)、寬和高)、設(shè)備質(zhì)量以及每臺(tái)設(shè)備加工1次的加工時(shí)間,具體信息見(jiàn)表1。

    根據(jù)產(chǎn)品的工藝性對(duì)產(chǎn)品進(jìn)行編號(hào),該車(chē)間的加工產(chǎn)品信息包括加工批量、搬運(yùn)批量以及加工產(chǎn)品所經(jīng)過(guò)設(shè)備的工藝路徑,具體信息見(jiàn)表2。

    3.2 結(jié)果比較和分析

    通過(guò)實(shí)例求解得到3種布局形式的最優(yōu)布局方案,解碼后得出了3種形狀的最優(yōu)單元布局方案。

    在直線形單元布局中,平面Z1有7臺(tái)設(shè)備,平面Z2有8臺(tái)設(shè)備。具體布局方案如圖5所示。

    在“U”形單元布局中,平面Z1有8臺(tái)設(shè)備,平面Z3有1臺(tái)設(shè)備,平面Z2有6臺(tái)設(shè)備。具體布局方案如圖6所示。

    在環(huán)形單元布局中,平面Z1有5臺(tái)設(shè)備,平面Z3有2臺(tái)設(shè)備,平面Z2有6臺(tái)設(shè)備,平面Z4有2臺(tái)設(shè)備。具體布局方案如圖7所示。

    SA-QPSO算法當(dāng)更新解時(shí)不會(huì)受目標(biāo)數(shù)值的數(shù)量級(jí)大小影響,而是比較2個(gè)解對(duì)應(yīng)目標(biāo)值的大小關(guān)系,這更有利于迭代尋優(yōu)。綜合實(shí)例求解結(jié)果,可以得出結(jié)論,將SA算法與QPSO算法相結(jié)合,在求解多層制造單元內(nèi)設(shè)備布局問(wèn)題中表現(xiàn)良好,尋優(yōu)能力更強(qiáng)。

    4 結(jié)語(yǔ)

    本文基于SA-QPSO算法提出了一種多層制造單元內(nèi)部布局優(yōu)化設(shè)計(jì)方法。該方法根據(jù)多層制造單元的多個(gè)目標(biāo)函數(shù)構(gòu)建嵌套優(yōu)化模型,進(jìn)行求解,并應(yīng)用于某汽車(chē)零件加工車(chē)間實(shí)際案例,得出直線形、“U”形以及環(huán)形3種最優(yōu)的空間布局方案,驗(yàn)證了當(dāng)較小規(guī)模的多層制造單元布局問(wèn)題求解時(shí),SA算法與QPSO算法相結(jié)合具有有效性。該方法可為相關(guān)的多層制造單元智能工廠提供一定的單元內(nèi)部布局設(shè)計(jì)依據(jù)。SA-QPSO算法的求解思路擴(kuò)展了QPSO算法的應(yīng)用范圍,可應(yīng)用于類(lèi)似問(wèn)題的求解。

    參考文獻(xiàn)

    [1]董峰.建設(shè)智能工廠推進(jìn)智能方案落地[J].中國(guó)工業(yè)和信息化,2022(1):44-47.

    [2]王沙沙,張則強(qiáng),劉俊琦,等.雙向環(huán)形過(guò)道布置問(wèn)題建模及混合鯨魚(yú)算法求解[J].計(jì)算機(jī)集成制造系統(tǒng),2021(10):27.

    [3]丁祥海,查穩(wěn),徐雙燕.“U”形生產(chǎn)單元可重構(gòu)設(shè)施布局方法研究[J].工業(yè)工程與管理,2015,20(1):6.

    [4]遲佳,梁秋艷,張曉玲,等.改進(jìn)量子粒子群自適應(yīng)優(yōu)化算法[J].中國(guó)科技信息,2023(11):86-90.

    [5]胡松琪.量子粒子群算法的研究及應(yīng)用[D].蘭州:蘭州理工大學(xué),2019.

    [6]張婧婧,施亭亭,汪強(qiáng).基于量子粒子群算法的公共空間主功能區(qū)布局優(yōu)化設(shè)計(jì)方法[J].齊齊哈爾大學(xué)學(xué)報(bào)(自然科學(xué)版),2023,39(6):66-71.

    作者簡(jiǎn)介:郭軼男(1998-),女,碩士研究生,主要研究方向?yàn)樯a(chǎn)車(chē)間布局優(yōu)化。

    電子郵箱:2470528100@qq.com。

    通信作者:姜雪松(1979-),男,博士,副教授,研究方向?yàn)楣I(yè)工程與管理、制造系統(tǒng)工程及信息化。

    電子郵箱:xuesongjiang@nefu.edu.cn。

    猜你喜歡
    模擬退火算法
    改進(jìn)模擬退火算法的K—means聚類(lèi)方法在學(xué)生成績(jī)上的應(yīng)用
    道路循環(huán)甩掛運(yùn)輸車(chē)輛調(diào)度研究
    改進(jìn)遺傳模擬退火算法求解TSP
    級(jí)聯(lián)型H橋逆變器的階梯波特定消諧技術(shù)研究
    科技資訊(2017年8期)2017-05-18 09:54:41
    基于圖像特征及改進(jìn)支持向量機(jī)算法的交通標(biāo)志識(shí)別
    模擬退火算法在整車(chē)物流問(wèn)題中的應(yīng)用
    物流科技(2016年12期)2017-04-01 03:12:04
    數(shù)學(xué)建模中的碎紙片拼接復(fù)原要點(diǎn)研究
    智能傳感器中的算法應(yīng)用
    改進(jìn)的模擬退火算法及其在裝填問(wèn)題中的應(yīng)用
    基于BP人工神經(jīng)網(wǎng)絡(luò)的離散型車(chē)間生產(chǎn)調(diào)度指標(biāo)預(yù)測(cè)模型的研究
    科技視界(2016年3期)2016-02-26 09:45:54
    白带黄色成豆腐渣| 精品久久久久久成人av| 国产精品av视频在线免费观看| 国内久久婷婷六月综合欲色啪| or卡值多少钱| 亚洲精品粉嫩美女一区| 国产亚洲精品第一综合不卡| 黑人巨大精品欧美一区二区mp4| 观看免费一级毛片| 桃色一区二区三区在线观看| 成人av在线播放网站| 免费观看人在逋| 亚洲狠狠婷婷综合久久图片| 老汉色av国产亚洲站长工具| 亚洲国产看品久久| 亚洲专区中文字幕在线| 久久久久久免费高清国产稀缺| 午夜日韩欧美国产| 国产精品永久免费网站| 老司机午夜福利在线观看视频| 欧美成人性av电影在线观看| 国产精品九九99| 香蕉丝袜av| 久久久久国内视频| 99国产精品一区二区蜜桃av| 国产91精品成人一区二区三区| 97人妻精品一区二区三区麻豆| 99国产精品99久久久久| 怎么达到女性高潮| 国产在线精品亚洲第一网站| 国产av一区在线观看免费| 国产片内射在线| 级片在线观看| 久久热在线av| 国产av麻豆久久久久久久| 午夜久久久久精精品| 国产aⅴ精品一区二区三区波| √禁漫天堂资源中文www| 国产成人啪精品午夜网站| 亚洲第一欧美日韩一区二区三区| 国产亚洲av嫩草精品影院| 69av精品久久久久久| 久久精品成人免费网站| 在线看三级毛片| 成年女人毛片免费观看观看9| 久久久精品国产亚洲av高清涩受| 中亚洲国语对白在线视频| 一级片免费观看大全| 久久久久久国产a免费观看| 在线永久观看黄色视频| 在线永久观看黄色视频| 免费av毛片视频| 变态另类丝袜制服| 国产亚洲精品久久久久5区| 国产伦人伦偷精品视频| 我要搜黄色片| 国产精品久久视频播放| 高清毛片免费观看视频网站| 久久性视频一级片| 很黄的视频免费| 后天国语完整版免费观看| 很黄的视频免费| 精品无人区乱码1区二区| 99热6这里只有精品| 狂野欧美白嫩少妇大欣赏| 婷婷丁香在线五月| 日韩欧美国产在线观看| 日韩欧美国产在线观看| 一级毛片精品| 久久精品亚洲精品国产色婷小说| 女人爽到高潮嗷嗷叫在线视频| 日本熟妇午夜| 国产精品99久久99久久久不卡| 搞女人的毛片| 成人国产综合亚洲| 精品国产超薄肉色丝袜足j| 999精品在线视频| ponron亚洲| 丁香欧美五月| 精品国产美女av久久久久小说| 嫩草影院精品99| 精品免费久久久久久久清纯| 久久精品成人免费网站| а√天堂www在线а√下载| 国产成人精品无人区| 久久香蕉激情| 中出人妻视频一区二区| 日日夜夜操网爽| 久久久久国产精品人妻aⅴ院| 好看av亚洲va欧美ⅴa在| 黄色成人免费大全| 18禁观看日本| av欧美777| 国产三级黄色录像| 久久精品成人免费网站| 国产单亲对白刺激| 国产精品电影一区二区三区| 亚洲熟女毛片儿| 国产免费av片在线观看野外av| а√天堂www在线а√下载| 亚洲成人中文字幕在线播放| 色播亚洲综合网| 日本在线视频免费播放| av免费在线观看网站| 黄色女人牲交| 久久九九热精品免费| 丰满人妻熟妇乱又伦精品不卡| 一级毛片精品| 日韩免费av在线播放| 久久香蕉国产精品| 97超级碰碰碰精品色视频在线观看| 国产精品电影一区二区三区| 亚洲真实伦在线观看| 丝袜美腿诱惑在线| 欧美国产日韩亚洲一区| 亚洲国产高清在线一区二区三| 1024香蕉在线观看| 欧美成狂野欧美在线观看| 日韩欧美国产在线观看| 老司机深夜福利视频在线观看| 神马国产精品三级电影在线观看 | 五月玫瑰六月丁香| 99热6这里只有精品| 久久精品亚洲精品国产色婷小说| 日韩精品免费视频一区二区三区| av超薄肉色丝袜交足视频| 在线播放国产精品三级| 久久精品91蜜桃| 床上黄色一级片| xxxwww97欧美| 亚洲一区中文字幕在线| 国产精品一及| 中文字幕精品亚洲无线码一区| 亚洲av中文字字幕乱码综合| 精华霜和精华液先用哪个| 一边摸一边抽搐一进一小说| 国产激情偷乱视频一区二区| 国产精品香港三级国产av潘金莲| 99久久无色码亚洲精品果冻| 精品一区二区三区av网在线观看| 亚洲 国产 在线| 国内毛片毛片毛片毛片毛片| 伦理电影免费视频| www.自偷自拍.com| 国产成人精品无人区| 日本撒尿小便嘘嘘汇集6| www国产在线视频色| 99国产极品粉嫩在线观看| 精华霜和精华液先用哪个| 亚洲人成77777在线视频| 亚洲色图 男人天堂 中文字幕| 国产成人精品无人区| 非洲黑人性xxxx精品又粗又长| 成人亚洲精品av一区二区| 久久精品亚洲精品国产色婷小说| 久久婷婷成人综合色麻豆| 黄色女人牲交| 亚洲一码二码三码区别大吗| 精品一区二区三区四区五区乱码| 熟女少妇亚洲综合色aaa.| 国产午夜精品论理片| 97碰自拍视频| 国产一区二区在线av高清观看| 午夜福利在线观看吧| 欧美另类亚洲清纯唯美| 国产精品乱码一区二三区的特点| 麻豆成人av在线观看| 国产精品,欧美在线| 久99久视频精品免费| 又黄又粗又硬又大视频| 国产欧美日韩一区二区三| 亚洲人成电影免费在线| 老熟妇乱子伦视频在线观看| 哪里可以看免费的av片| 色av中文字幕| 一进一出抽搐gif免费好疼| 舔av片在线| 欧洲精品卡2卡3卡4卡5卡区| 身体一侧抽搐| 欧美zozozo另类| 久久热在线av| 欧美在线黄色| 免费在线观看黄色视频的| 国产伦一二天堂av在线观看| 国产av不卡久久| 中出人妻视频一区二区| 18禁观看日本| 在线观看免费日韩欧美大片| 欧美激情久久久久久爽电影| 淫秽高清视频在线观看| 亚洲av成人精品一区久久| 午夜福利在线在线| 午夜福利在线观看吧| 99国产精品99久久久久| 99热只有精品国产| 色综合欧美亚洲国产小说| 一区二区三区激情视频| 久久中文字幕人妻熟女| 国产精品电影一区二区三区| 宅男免费午夜| 人妻久久中文字幕网| 我的老师免费观看完整版| 免费观看人在逋| 日本一本二区三区精品| 久久亚洲精品不卡| 女人被狂操c到高潮| 国产亚洲精品久久久久久毛片| 国产亚洲欧美在线一区二区| 欧美另类亚洲清纯唯美| 人人妻,人人澡人人爽秒播| 亚洲一码二码三码区别大吗| 亚洲精品av麻豆狂野| 97人妻精品一区二区三区麻豆| 日本 欧美在线| 真人做人爱边吃奶动态| 国产欧美日韩一区二区三| x7x7x7水蜜桃| 妹子高潮喷水视频| 亚洲九九香蕉| 丁香六月欧美| 国产一区在线观看成人免费| 免费在线观看黄色视频的| 亚洲国产精品合色在线| 嫩草影院精品99| 国产人伦9x9x在线观看| 久久人人精品亚洲av| 男女视频在线观看网站免费 | 国产精品av视频在线免费观看| 亚洲五月婷婷丁香| 激情在线观看视频在线高清| 午夜福利免费观看在线| 成人18禁高潮啪啪吃奶动态图| 久久久久久九九精品二区国产 | 欧美绝顶高潮抽搐喷水| 波多野结衣高清作品| 夜夜夜夜夜久久久久| 久久久久久免费高清国产稀缺| 俄罗斯特黄特色一大片| 精品国产乱子伦一区二区三区| 91字幕亚洲| 亚洲一码二码三码区别大吗| 亚洲一区二区三区不卡视频| 99在线人妻在线中文字幕| 日韩欧美在线二视频| 可以在线观看毛片的网站| avwww免费| 国产成人一区二区三区免费视频网站| 中文亚洲av片在线观看爽| 欧美性猛交╳xxx乱大交人| 观看免费一级毛片| 久久久国产欧美日韩av| 久久九九热精品免费| 亚洲自偷自拍图片 自拍| 国产精品一区二区三区四区久久| 人妻久久中文字幕网| 一区二区三区高清视频在线| a在线观看视频网站| 国产伦在线观看视频一区| 欧美黑人欧美精品刺激| 蜜桃久久精品国产亚洲av| 又粗又爽又猛毛片免费看| 国产成人影院久久av| 久久精品亚洲精品国产色婷小说| 1024香蕉在线观看| 亚洲成av人片免费观看| 国产亚洲精品综合一区在线观看 | av福利片在线| 特级一级黄色大片| 露出奶头的视频| 男女那种视频在线观看| 18禁黄网站禁片免费观看直播| 亚洲色图 男人天堂 中文字幕| 老司机深夜福利视频在线观看| 美女高潮喷水抽搐中文字幕| 成人18禁在线播放| 国产一区在线观看成人免费| 国产单亲对白刺激| 午夜福利18| 成年人黄色毛片网站| 国产精品98久久久久久宅男小说| 制服丝袜大香蕉在线| 人妻丰满熟妇av一区二区三区| 亚洲乱码一区二区免费版| 国产欧美日韩一区二区精品| 一区二区三区激情视频| 午夜亚洲福利在线播放| 无限看片的www在线观看| 国产一区二区在线av高清观看| 欧美日韩亚洲国产一区二区在线观看| 国产精品久久久久久亚洲av鲁大| 午夜免费激情av| 日韩精品中文字幕看吧| 亚洲欧美一区二区三区黑人| 韩国av一区二区三区四区| 日韩三级视频一区二区三区| 50天的宝宝边吃奶边哭怎么回事| 午夜福利欧美成人| 午夜福利在线观看吧| 精品久久久久久久久久久久久| 香蕉久久夜色| 国产精品影院久久| 听说在线观看完整版免费高清| 在线观看日韩欧美| 亚洲午夜理论影院| av在线播放免费不卡| 一本久久中文字幕| 亚洲中文日韩欧美视频| 在线视频色国产色| www.www免费av| 国产精品野战在线观看| 搡老岳熟女国产| 这个男人来自地球电影免费观看| 夜夜躁狠狠躁天天躁| 免费看a级黄色片| 精品久久久久久久末码| 欧美最黄视频在线播放免费| 搡老岳熟女国产| 欧美乱妇无乱码| 中文资源天堂在线| 99久久精品热视频| 可以在线观看的亚洲视频| 欧美高清成人免费视频www| 美女免费视频网站| 国产亚洲av嫩草精品影院| 天堂av国产一区二区熟女人妻 | 亚洲人成伊人成综合网2020| 国产精品电影一区二区三区| 亚洲国产看品久久| 在线观看免费午夜福利视频| www日本在线高清视频| 五月玫瑰六月丁香| 成人国语在线视频| 国产午夜福利久久久久久| 视频区欧美日本亚洲| 成熟少妇高潮喷水视频| 不卡一级毛片| 狠狠狠狠99中文字幕| 黄色a级毛片大全视频| 欧美av亚洲av综合av国产av| 亚洲五月婷婷丁香| 无限看片的www在线观看| 久久久久久久久久黄片| 国产69精品久久久久777片 | 中文亚洲av片在线观看爽| 天天添夜夜摸| 欧美成人免费av一区二区三区| 全区人妻精品视频| 亚洲性夜色夜夜综合| 老司机福利观看| 99久久久亚洲精品蜜臀av| 69av精品久久久久久| 日日摸夜夜添夜夜添小说| 一a级毛片在线观看| 18禁黄网站禁片午夜丰满| 午夜激情av网站| 99国产综合亚洲精品| 男女那种视频在线观看| 亚洲精品在线观看二区| 在线观看免费日韩欧美大片| 久久久久久久久免费视频了| av中文乱码字幕在线| 亚洲国产精品成人综合色| 脱女人内裤的视频| 国产三级中文精品| 国产亚洲精品久久久久5区| 1024手机看黄色片| 欧美最黄视频在线播放免费| 久久婷婷人人爽人人干人人爱| 给我免费播放毛片高清在线观看| 脱女人内裤的视频| av在线播放免费不卡| 国产精品美女特级片免费视频播放器 | 波多野结衣高清作品| 国产激情久久老熟女| 悠悠久久av| 母亲3免费完整高清在线观看| 欧美最黄视频在线播放免费| 1024香蕉在线观看| 欧美一级毛片孕妇| 久久欧美精品欧美久久欧美| 久久久久国内视频| 国产精品国产高清国产av| 丝袜美腿诱惑在线| 啦啦啦观看免费观看视频高清| 国产精品久久久av美女十八| 欧美日韩福利视频一区二区| 女同久久另类99精品国产91| 精品国产乱码久久久久久男人| 夜夜躁狠狠躁天天躁| 在线观看一区二区三区| 亚洲免费av在线视频| 精品电影一区二区在线| 十八禁人妻一区二区| 99国产精品一区二区三区| 久久久国产精品麻豆| 50天的宝宝边吃奶边哭怎么回事| 男人舔女人的私密视频| 日韩欧美免费精品| 高潮久久久久久久久久久不卡| av国产免费在线观看| 老司机午夜十八禁免费视频| 欧美日韩黄片免| 一级片免费观看大全| 在线十欧美十亚洲十日本专区| 国产亚洲欧美在线一区二区| 国产在线精品亚洲第一网站| 亚洲激情在线av| 欧美另类亚洲清纯唯美| 久久国产乱子伦精品免费另类| 好男人在线观看高清免费视频| 亚洲中文日韩欧美视频| 深夜精品福利| 国产av一区在线观看免费| 亚洲成人精品中文字幕电影| 1024视频免费在线观看| 精品高清国产在线一区| av欧美777| 99久久精品国产亚洲精品| 亚洲 欧美 日韩 在线 免费| 制服丝袜大香蕉在线| 亚洲 欧美一区二区三区| 男女下面进入的视频免费午夜| 欧美成狂野欧美在线观看| 不卡一级毛片| 亚洲欧美激情综合另类| 日韩 欧美 亚洲 中文字幕| 成熟少妇高潮喷水视频| 精品乱码久久久久久99久播| 成人一区二区视频在线观看| 黄色视频,在线免费观看| 亚洲一区二区三区不卡视频| 啪啪无遮挡十八禁网站| 日本黄色视频三级网站网址| 热99re8久久精品国产| 国产男靠女视频免费网站| 成人国语在线视频| 少妇被粗大的猛进出69影院| 中文字幕久久专区| 久久久精品欧美日韩精品| 在线观看舔阴道视频| 欧美成人性av电影在线观看| 全区人妻精品视频| 久久久精品国产亚洲av高清涩受| 午夜福利18| 精品一区二区三区四区五区乱码| 精品福利观看| 亚洲成人国产一区在线观看| 国内精品久久久久久久电影| 欧美日韩精品网址| 90打野战视频偷拍视频| 高清在线国产一区| 亚洲真实伦在线观看| 久久精品成人免费网站| 无限看片的www在线观看| 天堂动漫精品| 日本三级黄在线观看| 国产av不卡久久| 午夜福利在线在线| 999久久久精品免费观看国产| 免费观看人在逋| tocl精华| 99久久99久久久精品蜜桃| 国产野战对白在线观看| 久久天躁狠狠躁夜夜2o2o| 精品一区二区三区av网在线观看| 精品久久久久久久人妻蜜臀av| 亚洲精品在线观看二区| 欧美日韩乱码在线| www国产在线视频色| 欧美色视频一区免费| 黑人欧美特级aaaaaa片| 国产午夜福利久久久久久| 久久香蕉国产精品| 亚洲人成伊人成综合网2020| 国产成+人综合+亚洲专区| 国产免费av片在线观看野外av| 男人的好看免费观看在线视频 | 视频区欧美日本亚洲| 成人18禁在线播放| 三级国产精品欧美在线观看 | 九九热线精品视视频播放| 级片在线观看| 精品久久久久久,| 91国产中文字幕| 日韩欧美 国产精品| 亚洲色图 男人天堂 中文字幕| 不卡一级毛片| 非洲黑人性xxxx精品又粗又长| 又粗又爽又猛毛片免费看| 可以免费在线观看a视频的电影网站| 国产精品久久久av美女十八| 久久伊人香网站| 中文字幕人妻丝袜一区二区| 亚洲午夜理论影院| 99久久无色码亚洲精品果冻| 三级国产精品欧美在线观看 | 男男h啪啪无遮挡| 日本a在线网址| 小说图片视频综合网站| 久99久视频精品免费| 亚洲欧美日韩无卡精品| 国产片内射在线| 黄色视频不卡| 欧美绝顶高潮抽搐喷水| 日本 欧美在线| 黄色 视频免费看| 亚洲成a人片在线一区二区| 一区福利在线观看| 50天的宝宝边吃奶边哭怎么回事| 亚洲国产欧美人成| 欧美精品啪啪一区二区三区| 午夜精品在线福利| 国产亚洲精品第一综合不卡| 欧美国产日韩亚洲一区| 亚洲中文av在线| 欧美日韩瑟瑟在线播放| 欧美成狂野欧美在线观看| 人人妻,人人澡人人爽秒播| 美女扒开内裤让男人捅视频| www.999成人在线观看| 中文资源天堂在线| 免费看十八禁软件| 久久欧美精品欧美久久欧美| 国产精品久久电影中文字幕| 亚洲国产欧洲综合997久久,| 欧美乱色亚洲激情| 男插女下体视频免费在线播放| 91字幕亚洲| 欧美极品一区二区三区四区| 天堂√8在线中文| 特级一级黄色大片| 在线观看免费日韩欧美大片| 国产精品影院久久| 欧美国产日韩亚洲一区| 露出奶头的视频| 午夜激情av网站| 亚洲人成电影免费在线| 国产一区二区在线观看日韩 | 国产av又大| 日韩成人在线观看一区二区三区| 亚洲国产精品久久男人天堂| 国产精品日韩av在线免费观看| 国产伦一二天堂av在线观看| 桃红色精品国产亚洲av| 无遮挡黄片免费观看| 亚洲天堂国产精品一区在线| or卡值多少钱| 在线观看66精品国产| 女警被强在线播放| 99久久精品热视频| 亚洲中文日韩欧美视频| 欧美成人一区二区免费高清观看 | 久久久国产精品麻豆| 一二三四在线观看免费中文在| 毛片女人毛片| 国产爱豆传媒在线观看 | 深夜精品福利| 国产成+人综合+亚洲专区| 午夜精品在线福利| 午夜福利欧美成人| 成人欧美大片| 99久久无色码亚洲精品果冻| 女同久久另类99精品国产91| 亚洲精品久久国产高清桃花| 亚洲一区中文字幕在线| cao死你这个sao货| 热99re8久久精品国产| 日韩中文字幕欧美一区二区| 亚洲精品久久国产高清桃花| 在线a可以看的网站| 久久精品成人免费网站| 国产一区二区三区视频了| 亚洲人成电影免费在线| 国产亚洲精品第一综合不卡| 久久久国产成人精品二区| 18禁黄网站禁片免费观看直播| 午夜激情福利司机影院| 久久国产乱子伦精品免费另类| 美女 人体艺术 gogo| 日本五十路高清| 国产v大片淫在线免费观看| av国产免费在线观看| 麻豆成人av在线观看| 色综合站精品国产| 亚洲欧美日韩高清在线视频| 亚洲激情在线av| 国产在线精品亚洲第一网站| 国产精品一及| 50天的宝宝边吃奶边哭怎么回事| 国产私拍福利视频在线观看| 久久久久久久久免费视频了| 长腿黑丝高跟| 人人妻人人看人人澡| 国产黄色小视频在线观看| 免费无遮挡裸体视频| 日韩大码丰满熟妇| 99久久久亚洲精品蜜臀av| 欧美人与性动交α欧美精品济南到| 午夜视频精品福利| 色精品久久人妻99蜜桃| 亚洲电影在线观看av| 国产麻豆成人av免费视频| 搡老岳熟女国产| 欧美国产日韩亚洲一区| 九色成人免费人妻av| 91在线观看av| 国产成人精品久久二区二区91| 日韩精品中文字幕看吧| 中亚洲国语对白在线视频| 亚洲电影在线观看av| 精品熟女少妇八av免费久了| 人人妻人人澡欧美一区二区| 最新在线观看一区二区三区| 午夜福利18| 亚洲avbb在线观看|