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

    基于改進(jìn)鯨魚優(yōu)化算法的裝配線再平衡方法*

    2023-08-02 06:32:06閆炳南
    關(guān)鍵詞:裝配線鯨魚工作站

    郭 宏,閆炳南,羅 雷,冒 源,原 超,徐 壯

    (1.太原科技大學(xué)機械工程學(xué)院,太原 030024;2.山西平陽重工機械有限責(zé)任公司,侯馬 043003)

    0 引言

    裝配線平衡問題(assembly line balancing problem,ALBP)主要是研究在滿足各類約束的前提下,將各裝配任務(wù)分配到工作站,來實現(xiàn)減少成本和提升效率的目的[1]。

    當(dāng)前,已有許多學(xué)者對裝配線平衡問題進(jìn)行了較深入地研究。ALBUS等[2]采用一種基于混合整數(shù)規(guī)劃的線性優(yōu)化方法來解決裝配線的動態(tài)資源分配問題。KATIRAEE等[3]提出了一個雙目標(biāo)線性規(guī)劃模型來解決裝配線平衡問題。黃輝等[4]提出一種混合算法來解決多目標(biāo)雙邊裝配線平衡問題。趙聯(lián)鵬等[5]提出了一種新的搜索算法去求解裝配線平衡問題,以提高裝配線效率。

    現(xiàn)有的大多數(shù)模型都建立在確定需求的情況下,雖然一些模型考慮了需求的變化,但很少考慮再平衡成本的問題。李金霖等[6]使用多能工應(yīng)對需求變化的混裝線平衡問題,以最小化人工成本并滿足不同情境的需求。CHICA等[7]提出了一個用于裝配線平衡的多目標(biāo)模型,并使用自適應(yīng)進(jìn)化算法提高效率。張于賢等[8]為了解決算法的最優(yōu)解存在著達(dá)不到裝配線平衡最優(yōu)準(zhǔn)則且無法對瓶頸工序進(jìn)行操作分析的情形,提出了算法和IE技術(shù)兩階段組合求解法,求出了一個能滿足生產(chǎn)需要的有效可行解。

    針對不確定需求下的再平衡問題,本文采用只改變工作站數(shù)量以及重新規(guī)劃邊界的方式來建立裝配線模型,并提出了一種改進(jìn)多目標(biāo)鯨魚的優(yōu)化算法來求解該模型。最后,與其他多目標(biāo)算法進(jìn)行對比驗證了此算法的可行性。

    1 問題描述與數(shù)學(xué)建模

    1.1 可調(diào)整產(chǎn)能裝配線平衡問題描述

    當(dāng)訂單需求改變時,裝配線的生產(chǎn)周期也需要隨之變化。此時,有必要重新平衡裝配線。但為了減少成本,設(shè)施布局應(yīng)避免改動。在不改變裝配線上工位設(shè)備順序的前提下,可以通過增加或減少工作臺的數(shù)量以及重新規(guī)劃工作臺的邊界來實現(xiàn)不同的任務(wù)分配。如圖1所示,分為兩種情境,生產(chǎn)周期C=30和生產(chǎn)周期C=50,已知每個設(shè)備所分配任務(wù)時間,要求每個工作站所有設(shè)備任務(wù)時間之和不大于周期時間。當(dāng)生產(chǎn)周期較小時(圖中舉例為C=30),設(shè)備1和2分配給工作站1 的工人,設(shè)備3,4,5分配給工作站2的工人;當(dāng)需求減少時,生產(chǎn)周期增加(圖中示例為C=50),工作站1的工人操作設(shè)備1和2仍有余力,此時將設(shè)備3也交給此工人操作。

    本文所述的裝配線平衡問題可描述為:在已知兩個或多個不同需求情境的條件下,如何完成任務(wù)分配,以最小化多個需求下平均工作站數(shù)量和平均時間均衡指數(shù)。

    為有效描述該問題,做以下基本假設(shè):①裝配圖已給定且其中不存在回路;②為了減少成本,同一個任務(wù)只能分配給同一個設(shè)備;③每個任務(wù)都對應(yīng)一個設(shè)備;④可調(diào)整產(chǎn)能裝配線全部使用多能工進(jìn)行裝配;⑤普通工人的每天工作時間為T小時,多能工每天可用裝配時間設(shè)為T/2小時,多能工每天的工資水平是普通工人的2倍。

    1.2 符號說明

    N為裝配順序圖中任務(wù)的數(shù)量,N*為所有任務(wù)所組成的集合,S*為所有需求情境組成的集合,Ps為需求情境s發(fā)生的比例,Ds為情境s下每天的需求數(shù)量(單位:件/天),T為1個工人1天可使用的工作時間,Cs為需求情境s的生產(chǎn)節(jié)拍,Cs=T/Ds;Tsj為需求情境s下第j個工作站的總體任務(wù)時間,Lsj為情境s下第j個工作站的時間均衡指數(shù),Lsj=(Cs-Tsj)/Cs;H*為已使用設(shè)備所組成的集合,xih為任務(wù)i被分到第h個設(shè)備,則為1,否則為0;yhjs為設(shè)備h在情境s被分到第j個工作臺,則為1,否則為0。

    1.3 建立數(shù)學(xué)模型

    建立變產(chǎn)能裝配線數(shù)學(xué)模型,其目標(biāo)函數(shù)為:

    minF={f1,f2}

    (1)

    式中:

    (2)

    (3)

    約束條件如下:

    (4)

    (5)

    (6)

    (7)

    (8)

    (9)

    xih∈(0,1),?i∈N*,?h∈H*

    (10)

    yhjs∈(0,1),?h∈H*,?s∈S*

    (11)

    式(1)~式(3)表示最小化每個場景中使用的工作站數(shù)量的平均值,以及每個場景下每個工作站的時間平衡指數(shù)。式(4)表示每個任務(wù)都對應(yīng)一個設(shè)備。式(5)表示每個站點上都只有一個設(shè)備。式(6)表示任務(wù)優(yōu)先順序約束,每個任務(wù)在生產(chǎn)線上的工位必須在它所有前置任務(wù)完成之后。式(7)表示在給定情境下每個位置只能分配到一個工作站,不同的工作站之間不能共享同一個設(shè)備。式8)表示每個工作站的所包含的區(qū)域必須是連續(xù)。式(9)表示生產(chǎn)節(jié)拍約束,這意味著必須完成每種情況下所需的任務(wù)。式(10)和式(11)分別表示xij和yhjs的值范圍。

    2 算法構(gòu)建

    本文研究的裝配線平衡問題屬于NP難問題,而群智能算法很適合解決這一問題。鯨魚優(yōu)化算法(WOA)是一種群智能算法,具有結(jié)構(gòu)簡單、計算方便的優(yōu)點,在工程領(lǐng)域已廣泛應(yīng)用[9]。

    為更好地求解變產(chǎn)能裝配線模型,本文提出了一種改進(jìn)多目標(biāo)鯨魚優(yōu)化算法。首先,使用混沌映射初始化種群;其次,在位置更新時對搜索步長進(jìn)行動態(tài)調(diào)整來提高算法的效率,并使用新的頭鯨選擇方式來避免局部最優(yōu);最后,采用改進(jìn)擁擠距離來更新外部檔案以得到最優(yōu)解集。

    2.1 鯨魚優(yōu)化算法(WOA)

    鯨魚優(yōu)化算法(WOA)是一種新的元啟發(fā)式優(yōu)化算法,WOA的主要靈感來自座頭鯨的特殊捕食行為。鯨魚在獵物周圍以螺旋形移動,在覓食時使用縮小的圓圈向獵物移動,這種行為被稱為泡泡網(wǎng)覓食。

    2.1.1 包圍捕食階段

    WOA通過在螺旋模型和收縮包圍獵物之間進(jìn)行50%的選擇概率來模擬這種狩獵機制,以生成當(dāng)前鯨魚的新位置。通常,包圍機制的數(shù)學(xué)表達(dá)式為:

    D=|C·X*(t)-X(t)|

    (12)

    X(t+1)=X*(t)-A·D

    (13)

    A=2a×r-a

    (14)

    C=2×r

    (15)

    (16)

    式中:X為位置矢量,X*(t)為一個向量,包含迄今為止找到的最佳位置;t為當(dāng)前迭代次數(shù),||為絕對值運算,·為逐個元素相乘,A和C為系數(shù)變量,a為一個隨著迭代從2不斷線性減小到0的參數(shù),r為[0,1]范圍內(nèi)隨機生成的數(shù)。

    2.1.2 螺旋捕食階段

    WOA試圖用螺旋模型來模擬鯨魚的螺旋形運動,以得出獵物和鯨魚之間的位置。螺旋形狀的數(shù)學(xué)表達(dá)式為:

    D=|X*(t)-X(t)|

    (17)

    X(t+1)=D*·ebl·cos(2πl(wèi))+X*(t)

    (18)

    式中:b為一個常數(shù),用來描述對數(shù)螺旋運動;l為一個1和-1范圍內(nèi)的隨機數(shù)。

    2.1.3 隨機搜索階段

    為了使WOA能夠探索優(yōu)化問題搜索空間內(nèi)的大部分區(qū)域,它從當(dāng)前種群中隨機選擇一個個體,將當(dāng)前的鯨魚移向它。如果|A|>1,則當(dāng)前鯨魚將移向從種群中隨機選擇的鯨魚。該階段的數(shù)學(xué)公式為:

    X(t+1)=X*(t)-A·D

    (19)

    D=|C·Xr(t)-X(t)|

    (20)

    式中:Xr(t)為從當(dāng)前種群中隨機選擇的個體。

    2.2 頭鯨選擇方式

    在WOA的進(jìn)化過程中,頭鯨在逃離局部最優(yōu)解方面發(fā)揮著重要作用。所以本文通過提出一種成對競爭機制來選擇更好的頭鯨,以其來指導(dǎo)鯨魚種群的更新。

    在成對競爭機制中,鯨魚個體通過競爭機制更新,而不是使用全局和單個最優(yōu)個體,從而大大提高了種群的多樣性,避免過早收斂。具體來說,在成對競爭機制中,從當(dāng)前種群中隨機配對選擇個體進(jìn)行競爭,并通過向優(yōu)勝者學(xué)習(xí)來更新競爭中的失敗者,而優(yōu)勝者則直接傳遞給下一代種群。

    假設(shè)種群P(t)中有N個個體,每個個體都有一個n維位置Xi(t)=(xi,1(t),xi,2(t),…,xi,n(t))和n維速度向量Vi(t)=(vi,1(t),vi,2(t),…,vi,n(t))。在第k輪比賽中,勝利者和失敗者的位置和速度分別表示為Xw,k(t)、Xl,k(t)、Vw,k(t)和Vl,k(t),其中k=1,2,…,N/2。在第k輪比賽后,失敗者的速度將使用以下策略進(jìn)行更新:

    (21)

    而失敗者的位置可以根據(jù)新的速度更新:

    Xl,k(t+1)=Xl,k(t)+Vl,k(t+1)

    (22)

    在創(chuàng)建精英種群之后,將在其中的鯨魚個體之間進(jìn)行成對競爭,獲勝者將被用來指導(dǎo)當(dāng)前種群中鯨魚個體的移動方向。對于每一次競爭,給定要更新的種群中的一個鯨魚個體p,從精英種群中隨機選擇兩個精英鯨魚個體a和b。分別計算a、b和p之間的角度,角度較小的精英個體獲勝,這樣鯨魚個體p將從更接近其收斂方向的精英個體學(xué)習(xí)。具體示例如圖2所示。

    圖2 成對競爭機制示例

    圖2中Xi是要更新的解決方案;X1和X2是從精英列表中隨機選擇的精英個體;θ1是X1和Xi之間的角度;θ2是X2和Xi之間的角度。因為θ2<θ1;X2是獲勝者,被選為領(lǐng)導(dǎo)者。

    2.3 外部檔案更新策略

    外部歸檔是多目標(biāo)鯨魚算法(MOWOA)中常用的模塊,在更新外部檔案的策略方面,有許多選擇,如網(wǎng)格選擇、非支配排序、高效非支配排序(ENS)和擁擠距離法等。

    由于擁擠距離法可以更好地維持種群中個體的多樣性,所以本文使用擁擠距離法更新外部檔案。然而,使用傳統(tǒng)的擁擠距離來反映密度在某種程度上是有缺陷的,所以本文采用一種新的方法來計算距離。其具體方法如圖3所示。

    圖3 擁擠距離

    圖3中有3種解決方案,用x1,x2,x3表示;目標(biāo)數(shù)為2。解x2的傳統(tǒng)擁擠距離為CD(x2)=R1+R2+S1+S2。事實上,只要x2位于邊界框中,其中x1和x3是對角頂點,CD(x2)的值就不會改變。為了解決這個問題,本文使用一種改進(jìn)的擁擠距離,其表達(dá)式定義為[10]:

    (23)

    式中:M為目標(biāo)數(shù),Si和Ri為圖3所示的距離。

    2.4 多策略改進(jìn)多目標(biāo)鯨魚優(yōu)化算法

    原始WOA算法通常使用隨機生成初始種群的方法。然而根據(jù)現(xiàn)有研究,初始種群將直接影響算法的計算效率。因為混沌序列有較強的隨機性以及整體穩(wěn)定的特性,所以本文使用混沌序列初始化。其中,Tent映射公式為:

    (24)

    式中:X(t)∈[0,1]。

    在WOA中,因子a是一個距離控制參數(shù),從2線性減少到0。因子a隨固定步長變化,因此當(dāng)?shù)螖?shù)較大時,大多數(shù)迭代將在找到最佳解之前結(jié)束,否則將消耗大量迭代次數(shù)而沒有任何進(jìn)展。所以,本文設(shè)計了新的動態(tài)控制系數(shù),使其與迭代次數(shù)成比例減少。其數(shù)學(xué)表達(dá)式如下:

    (25)

    (26)

    式中:r是一個介于0和1之間的隨機數(shù)。在前面的方程中,我們試圖控制距離控制參數(shù)a以值2開始,該值隨著迭代逐漸減小,并乘以0和1之間的隨機數(shù)r,直到達(dá)到標(biāo)準(zhǔn)參數(shù)無法達(dá)到的區(qū)域。

    步驟1:初始化。首先初始化外部檔案,設(shè)置初始容量為C,并在滿足優(yōu)先關(guān)系前提下使用式(24)生成N個初始個體X1,X2,…,XN;

    步驟2:確定非支配解集。對每一個鯨魚個體位置Pi計算適應(yīng)度函數(shù),并對其進(jìn)行Pareto非支配排序,得到非支配解集;

    步驟3:更新外部檔案。使用式(23)計算前C個個體,并將其加入外部檔案;

    步驟4:頭鯨選擇。使用競爭機制進(jìn)行頭鯨選擇, 首先,根據(jù)擁擠距離排序從外部檔案中選擇出前10個個體并組成集合E,從E中隨機選擇兩個個體X1和X2;然后,從當(dāng)前種群中取出一個個體Xi,計算X1和Xi之間的角度θ1以及X2和Xi之間的角度θ2,如果θ1<θ2,X1為領(lǐng)導(dǎo)者X*,否則,X2為領(lǐng)導(dǎo)者X*;之后,將X*代入到式(13)或式(18)或式(19)中得到X(t+1),并將其添加到集合X′,最后,判斷集合X是否為空,為空則重新從E中隨機選擇兩個個體循環(huán)執(zhí)行,否則返回新種群X′;

    步驟5:鯨魚位置更新。根據(jù)式(25)和式(26)更新參數(shù)a和C,并相應(yīng)地更新算法中A、l和r等相關(guān)參數(shù)。然后根據(jù)p和|A|的大小來選擇相應(yīng)的更新方案,并使用步驟4選出的頭鯨來指導(dǎo)更新。

    MIMOWOA算法流程圖如圖4所示。

    圖4 MIMOWOA算法流程圖

    3 仿真實驗分析

    3.1 參數(shù)分析

    使用python實現(xiàn)上述算法程序,并在CPU為CMD 6-Core,內(nèi)存為8 GB的PC上進(jìn)行。實驗參數(shù)設(shè)定為種群規(guī)模30,迭代次數(shù)1000。實驗選用SALBP-1標(biāo)準(zhǔn)算例庫中的9個不同規(guī)模的數(shù)據(jù)進(jìn)行實驗。

    3.2 實驗對比

    為了驗證此算法的優(yōu)越性,將它與原有的多目標(biāo)鯨魚優(yōu)化算法進(jìn)行比較,求解結(jié)果如表1所示。

    表1 不同問題規(guī)模下的結(jié)果對比

    表1中分為小,中,大3種不同任務(wù)數(shù)量的任務(wù),給定相應(yīng)的周期時間,對比MOWOA和MIMOWOA 兩種算法對同一問題的計算結(jié)果。從表中可以明顯看出,MIMOWOA算法得出的站點數(shù)量和均衡指數(shù)都低于MOWOA算法。如圖5所示,為了更加直觀地體現(xiàn)出優(yōu)越性,對比任務(wù)數(shù)量為101的帕累托前沿。其中,MIMOWOA在兩個目標(biāo)函數(shù)方向都能實現(xiàn)帕累托最優(yōu)。這主要由于MIMOWOA采用動態(tài)距離控制系數(shù)以及成對競爭機制后,在一定程度上提高了裝配線的均衡性。

    圖5 任務(wù)數(shù)量為101的帕累托前沿

    為了進(jìn)一步驗證算法的有效性,將改進(jìn)算法與其他多目標(biāo)優(yōu)化算法進(jìn)行比較。其中將算法性能的評價指標(biāo)設(shè)定為帕累托解的個數(shù)(NF)、解的均勻性(ES)和非支配解的最優(yōu)度(DPS)[11]。各項評價指標(biāo)計算為:

    (1)非支配解的數(shù)量:

    NF=|Np|

    (27)

    式中:NF為算法得到的帕累托解的數(shù)量,NF值越大,算法越優(yōu)。

    (2)非支配解最優(yōu)度:

    (28)

    (3)解的均勻性指標(biāo):

    (29)

    (30)

    式中:Di為xi與其他相鄰解的距離。ES越小,說明解越分散。

    實驗包含3種不同規(guī)模的任務(wù)數(shù)量,每個任務(wù)重復(fù)15次,選取平均值作為最終結(jié)果如表2所示。MW為改進(jìn)多目標(biāo)鯨魚算法,NS為遺傳算法,MP為粒子群算法,加粗字體表示不同問題規(guī)模下各項評價指標(biāo)的最優(yōu)值。

    表2 3種算法的數(shù)值結(jié)果對比

    分析表2可得:

    (1)分析評價性能NF的值可知,MIMOWOA得到的帕累托解的數(shù)量明顯更多。

    (2)分析評價性能DPS的值可知,MIMOWOA得到的DPS值較其余算法更大一些,表明MIMOWOA的全局搜索性更好。

    (3)分析評價性能ES的值可知,MIMOWOA得出的解的分布與其他算法相比較均勻。這是因為MIMWOA采用了擁擠距離排序的外部檔案更新策略,使解的分布更均勻。

    總而言之,MIMOWOA算法在解的數(shù)量或者質(zhì)量方面都明顯優(yōu)于其余兩個算法,這也驗證說明了其求解ALBP問題時的可行性。

    4 結(jié)束語

    本文研究了可調(diào)整產(chǎn)能裝配線平衡問題,采用只改變工作站數(shù)量以及重新規(guī)劃邊界的方式來建立裝配線模型,降低了裝配線再平衡成本。并針對該模型構(gòu)建了一種改進(jìn)多目標(biāo)鯨魚優(yōu)化算法,在原始鯨魚算法中使用混沌映射初始化種群,并在鯨魚搜索過程中引入動態(tài)距離控制系數(shù)和成對競爭機制來提高算法的收斂速度。最后用實驗驗證了所提出的算法在解決可調(diào)整產(chǎn)能裝配線平衡問題上的可行性和有效性。

    當(dāng)然本文只考慮了在兩種周期下進(jìn)行轉(zhuǎn)換的裝配線方案,后續(xù)可以進(jìn)一步開展多周期轉(zhuǎn)換的研究。

    猜你喜歡
    裝配線鯨魚工作站
    小鯨魚
    幼兒100(2022年41期)2022-11-24 03:20:20
    左權(quán)浙理大 共建工作站
    汽車零部件自動化裝配線防錯設(shè)計
    汽車工藝師(2021年7期)2021-07-30 08:03:26
    迷途鯨魚
    戴爾Precision 5750移動工作站
    電腦報(2020年32期)2020-09-06 13:55:22
    基于SPS模式的轉(zhuǎn)向架軸箱裝配線仿真研究
    鯨魚
    鯨魚島——拖延癥
    動漫星空(2018年4期)2018-10-26 02:11:54
    混流裝配線第二類平衡問題優(yōu)化研究
    基于Flexsim的隨機混流裝配線平衡設(shè)計與仿真
    日韩精品青青久久久久久| 人成视频在线观看免费观看| 国产精品爽爽va在线观看网站 | 国产精品,欧美在线| 日韩三级视频一区二区三区| 国产蜜桃级精品一区二区三区| 免费一级毛片在线播放高清视频| 欧美日韩乱码在线| 国产午夜精品久久久久久| 亚洲男人天堂网一区| 桃红色精品国产亚洲av| 午夜激情福利司机影院| 悠悠久久av| 波多野结衣av一区二区av| 亚洲成人精品中文字幕电影| 男女做爰动态图高潮gif福利片| 国产精品自产拍在线观看55亚洲| 性色av乱码一区二区三区2| 精品免费久久久久久久清纯| 19禁男女啪啪无遮挡网站| 国产91精品成人一区二区三区| 国产一卡二卡三卡精品| 亚洲熟妇熟女久久| 51午夜福利影视在线观看| 日韩av在线大香蕉| 午夜亚洲福利在线播放| 两个人视频免费观看高清| 免费搜索国产男女视频| 久久中文字幕人妻熟女| 国产精品 国内视频| 日日干狠狠操夜夜爽| 香蕉丝袜av| 国产亚洲精品久久久久5区| 成人午夜高清在线视频 | 黄色毛片三级朝国网站| 中文字幕精品免费在线观看视频| 变态另类成人亚洲欧美熟女| 成人av一区二区三区在线看| 99在线人妻在线中文字幕| 成年人黄色毛片网站| 久久国产精品人妻蜜桃| 欧美日韩黄片免| 久久精品国产清高在天天线| 久久久久国内视频| 最近最新中文字幕大全电影3 | 国产精品自产拍在线观看55亚洲| 日本撒尿小便嘘嘘汇集6| 亚洲av片天天在线观看| 国产成人系列免费观看| 久久精品91无色码中文字幕| 99热6这里只有精品| 法律面前人人平等表现在哪些方面| av片东京热男人的天堂| 夜夜夜夜夜久久久久| 一本精品99久久精品77| 法律面前人人平等表现在哪些方面| 午夜成年电影在线免费观看| 精品国内亚洲2022精品成人| 国产激情偷乱视频一区二区| 黄色a级毛片大全视频| 国产精品亚洲av一区麻豆| 老司机福利观看| 两个人看的免费小视频| 麻豆一二三区av精品| 曰老女人黄片| ponron亚洲| 少妇被粗大的猛进出69影院| 亚洲久久久国产精品| 香蕉丝袜av| 99国产精品一区二区蜜桃av| 又大又爽又粗| 十分钟在线观看高清视频www| 一卡2卡三卡四卡精品乱码亚洲| 最新美女视频免费是黄的| 一级毛片精品| 国产乱人伦免费视频| 男女做爰动态图高潮gif福利片| 亚洲精华国产精华精| svipshipincom国产片| 欧美精品亚洲一区二区| 十八禁网站免费在线| 久9热在线精品视频| 脱女人内裤的视频| 午夜福利18| 欧美日韩乱码在线| √禁漫天堂资源中文www| 亚洲欧美精品综合一区二区三区| 老汉色av国产亚洲站长工具| 精品熟女少妇八av免费久了| 久久精品亚洲精品国产色婷小说| 国产在线观看jvid| 午夜免费激情av| 成人18禁高潮啪啪吃奶动态图| 身体一侧抽搐| 一个人免费在线观看的高清视频| 精品日产1卡2卡| 手机成人av网站| 国产免费av片在线观看野外av| 黄片大片在线免费观看| 中文字幕高清在线视频| 99re在线观看精品视频| 欧美黑人巨大hd| 日韩欧美在线二视频| 国产精品亚洲一级av第二区| 黄色视频不卡| 久久久国产成人精品二区| www日本在线高清视频| 日韩 欧美 亚洲 中文字幕| 亚洲五月天丁香| xxx96com| 午夜免费成人在线视频| 国产精品 国内视频| xxx96com| 亚洲真实伦在线观看| 两个人视频免费观看高清| 又紧又爽又黄一区二区| av有码第一页| 国产单亲对白刺激| 国产精品99久久99久久久不卡| 亚洲成人国产一区在线观看| 亚洲欧美激情综合另类| 国产精品98久久久久久宅男小说| 欧美激情 高清一区二区三区| 真人一进一出gif抽搐免费| 看黄色毛片网站| 国产野战对白在线观看| 丝袜人妻中文字幕| 看片在线看免费视频| 无遮挡黄片免费观看| 国产精品 国内视频| 人人澡人人妻人| 亚洲全国av大片| 国内毛片毛片毛片毛片毛片| 正在播放国产对白刺激| 精品人妻1区二区| 亚洲,欧美精品.| 久久精品成人免费网站| 午夜成年电影在线免费观看| 欧美国产精品va在线观看不卡| 欧美av亚洲av综合av国产av| 国产色视频综合| 看片在线看免费视频| 美女大奶头视频| 国产精品久久视频播放| 在线国产一区二区在线| 精华霜和精华液先用哪个| 欧美成人一区二区免费高清观看 | 欧美成人免费av一区二区三区| 精品少妇一区二区三区视频日本电影| 啦啦啦 在线观看视频| 亚洲欧美激情综合另类| 国产精品一区二区精品视频观看| 欧美日韩中文字幕国产精品一区二区三区| 国产精品,欧美在线| 午夜精品在线福利| 国产伦在线观看视频一区| avwww免费| 久久久久九九精品影院| 一本综合久久免费| 中文字幕人妻丝袜一区二区| 成年免费大片在线观看| 一区二区三区激情视频| 国产亚洲欧美在线一区二区| 国产久久久一区二区三区| 黄色成人免费大全| 国产真人三级小视频在线观看| 免费电影在线观看免费观看| 香蕉久久夜色| 一级黄色大片毛片| 日韩 欧美 亚洲 中文字幕| 制服丝袜大香蕉在线| 成年免费大片在线观看| 成人一区二区视频在线观看| 亚洲人成网站在线播放欧美日韩| 精品人妻1区二区| 男女视频在线观看网站免费 | av电影中文网址| 婷婷六月久久综合丁香| 日韩一卡2卡3卡4卡2021年| 手机成人av网站| 亚洲熟妇熟女久久| 久久久久精品国产欧美久久久| 一进一出抽搐gif免费好疼| 亚洲精品中文字幕一二三四区| 观看免费一级毛片| 亚洲精品一区av在线观看| cao死你这个sao货| netflix在线观看网站| 久久99热这里只有精品18| 久99久视频精品免费| 岛国在线观看网站| 一个人免费在线观看的高清视频| 法律面前人人平等表现在哪些方面| 少妇粗大呻吟视频| 日韩大码丰满熟妇| 老司机深夜福利视频在线观看| 宅男免费午夜| 成年女人毛片免费观看观看9| 黄网站色视频无遮挡免费观看| 免费高清在线观看日韩| 精品久久久久久久久久免费视频| 欧美av亚洲av综合av国产av| 国产成人av教育| 欧美日韩精品网址| 丝袜美腿诱惑在线| 脱女人内裤的视频| 日本 av在线| 女警被强在线播放| 妹子高潮喷水视频| 老鸭窝网址在线观看| 一级毛片女人18水好多| 一进一出好大好爽视频| 最新在线观看一区二区三区| 久久久久久人人人人人| 免费看十八禁软件| 欧美日韩瑟瑟在线播放| 啦啦啦 在线观看视频| 亚洲男人的天堂狠狠| ponron亚洲| av超薄肉色丝袜交足视频| 成人精品一区二区免费| 亚洲国产欧洲综合997久久, | 岛国视频午夜一区免费看| 国产区一区二久久| 亚洲专区中文字幕在线| 国产色视频综合| 欧美精品啪啪一区二区三区| 亚洲精品国产精品久久久不卡| 国产精品久久电影中文字幕| 亚洲欧美日韩无卡精品| x7x7x7水蜜桃| 99久久久亚洲精品蜜臀av| 亚洲一码二码三码区别大吗| 在线观看舔阴道视频| 免费电影在线观看免费观看| 欧美中文综合在线视频| 18禁黄网站禁片免费观看直播| 十八禁人妻一区二区| 看免费av毛片| 亚洲人成网站在线播放欧美日韩| 女人高潮潮喷娇喘18禁视频| 成在线人永久免费视频| 在线看三级毛片| 日日摸夜夜添夜夜添小说| 日韩欧美 国产精品| 男女床上黄色一级片免费看| 亚洲精华国产精华精| 日本三级黄在线观看| 免费无遮挡裸体视频| 精华霜和精华液先用哪个| 一夜夜www| 久久久国产成人精品二区| 国产免费男女视频| 国产又色又爽无遮挡免费看| xxxwww97欧美| 18禁观看日本| 国产欧美日韩一区二区精品| 亚洲 欧美 日韩 在线 免费| 一边摸一边抽搐一进一小说| 国产精品国产高清国产av| 国产黄a三级三级三级人| 欧美中文综合在线视频| 久久国产精品男人的天堂亚洲| 精品不卡国产一区二区三区| 国产精品av久久久久免费| 级片在线观看| 久热这里只有精品99| 国产区一区二久久| 婷婷精品国产亚洲av在线| 两性夫妻黄色片| 亚洲国产精品合色在线| 精品免费久久久久久久清纯| 一级毛片高清免费大全| 亚洲熟妇熟女久久| 国内精品久久久久精免费| 国产av一区二区精品久久| 成年免费大片在线观看| 神马国产精品三级电影在线观看 | 最近最新中文字幕大全电影3 | av片东京热男人的天堂| 日本免费一区二区三区高清不卡| 久久午夜亚洲精品久久| 国产亚洲精品久久久久久毛片| 国产黄色小视频在线观看| 亚洲av中文字字幕乱码综合 | 脱女人内裤的视频| 天堂动漫精品| 欧美久久黑人一区二区| 国产免费av片在线观看野外av| 女人被狂操c到高潮| 精品电影一区二区在线| 男人操女人黄网站| 精品一区二区三区视频在线观看免费| 日韩高清综合在线| 777久久人妻少妇嫩草av网站| 欧美日韩亚洲综合一区二区三区_| 久久国产亚洲av麻豆专区| 深夜精品福利| 精品久久久久久久久久免费视频| 久久久久久久久久黄片| 久久久久精品国产欧美久久久| 国产真人三级小视频在线观看| xxxwww97欧美| 最好的美女福利视频网| xxxwww97欧美| 欧美日本亚洲视频在线播放| 国产区一区二久久| 无遮挡黄片免费观看| 久久久精品国产亚洲av高清涩受| 日韩精品青青久久久久久| 国产av又大| 又紧又爽又黄一区二区| 欧美成人一区二区免费高清观看 | 亚洲一区二区三区不卡视频| 嫩草影院精品99| 国产激情欧美一区二区| 久久香蕉国产精品| 99精品在免费线老司机午夜| 国产一区二区在线av高清观看| 午夜久久久久精精品| 久久狼人影院| 亚洲精品国产一区二区精华液| 丁香欧美五月| 成熟少妇高潮喷水视频| 99re在线观看精品视频| 国产激情久久老熟女| 美女大奶头视频| 久久久久久亚洲精品国产蜜桃av| 久久人人精品亚洲av| bbb黄色大片| 亚洲国产欧美日韩在线播放| 国产熟女午夜一区二区三区| 亚洲国产欧美日韩在线播放| 日韩有码中文字幕| 好男人电影高清在线观看| 好男人在线观看高清免费视频 | 亚洲成a人片在线一区二区| 久久久久久九九精品二区国产 | 丁香六月欧美| av有码第一页| 国产亚洲精品综合一区在线观看 | x7x7x7水蜜桃| 91成人精品电影| 国产视频内射| 日韩成人在线观看一区二区三区| 久久久久国产一级毛片高清牌| 久久人人精品亚洲av| xxxwww97欧美| 成人三级黄色视频| 熟女少妇亚洲综合色aaa.| 一级a爱片免费观看的视频| e午夜精品久久久久久久| 自线自在国产av| 日本一区二区免费在线视频| 日韩精品青青久久久久久| 无人区码免费观看不卡| 男女午夜视频在线观看| www国产在线视频色| 波多野结衣av一区二区av| 成人国语在线视频| 日韩 欧美 亚洲 中文字幕| 成熟少妇高潮喷水视频| 国产蜜桃级精品一区二区三区| 日韩欧美在线二视频| 亚洲精品久久成人aⅴ小说| 欧美另类亚洲清纯唯美| 国产免费av片在线观看野外av| 欧美国产日韩亚洲一区| 亚洲激情在线av| 亚洲第一电影网av| 国产私拍福利视频在线观看| 日韩大尺度精品在线看网址| 成人免费观看视频高清| 国产亚洲精品一区二区www| 极品教师在线免费播放| www日本在线高清视频| 国产高清有码在线观看视频 | 黄色 视频免费看| 精品欧美国产一区二区三| 麻豆久久精品国产亚洲av| 麻豆成人午夜福利视频| 欧美+亚洲+日韩+国产| 精品无人区乱码1区二区| 午夜福利在线在线| 欧美色欧美亚洲另类二区| av中文乱码字幕在线| 亚洲人成伊人成综合网2020| 天天躁狠狠躁夜夜躁狠狠躁| 91麻豆av在线| 亚洲精品久久成人aⅴ小说| 国内久久婷婷六月综合欲色啪| 色综合亚洲欧美另类图片| 久久久久精品国产欧美久久久| 国产欧美日韩一区二区三| 中文在线观看免费www的网站 | 1024手机看黄色片| 此物有八面人人有两片| 中出人妻视频一区二区| 日韩 欧美 亚洲 中文字幕| 午夜久久久在线观看| 两个人免费观看高清视频| 色尼玛亚洲综合影院| 国产私拍福利视频在线观看| 亚洲 欧美一区二区三区| 成人免费观看视频高清| 又黄又粗又硬又大视频| 日本 欧美在线| 国产亚洲av高清不卡| 91国产中文字幕| 老司机靠b影院| 可以在线观看的亚洲视频| 午夜福利在线在线| 特大巨黑吊av在线直播 | 一级毛片女人18水好多| 精品久久久久久久毛片微露脸| 国产熟女xx| 日本a在线网址| 免费无遮挡裸体视频| 国产成人影院久久av| 亚洲精品中文字幕在线视频| 很黄的视频免费| 亚洲国产精品sss在线观看| 亚洲男人的天堂狠狠| 亚洲激情在线av| 精品国产乱子伦一区二区三区| 女性生殖器流出的白浆| 一本综合久久免费| 亚洲五月天丁香| 亚洲电影在线观看av| 在线观看日韩欧美| 2021天堂中文幕一二区在线观 | 一二三四社区在线视频社区8| 久久青草综合色| 欧美性猛交黑人性爽| 一级毛片高清免费大全| а√天堂www在线а√下载| 国产黄片美女视频| 国产色视频综合| 亚洲 国产 在线| 首页视频小说图片口味搜索| 国产主播在线观看一区二区| 国内精品久久久久精免费| a级毛片a级免费在线| 亚洲自偷自拍图片 自拍| 日本a在线网址| 变态另类成人亚洲欧美熟女| 中文字幕人妻熟女乱码| 久久人人精品亚洲av| 99热这里只有精品一区 | 老司机靠b影院| 国产成人av教育| 亚洲国产精品sss在线观看| 亚洲国产精品久久男人天堂| 一级毛片精品| 亚洲第一青青草原| 国产精品野战在线观看| 国产黄色小视频在线观看| 一级a爱视频在线免费观看| 久久九九热精品免费| 自线自在国产av| 欧美一级毛片孕妇| 亚洲av美国av| 国产熟女xx| 国产精品 国内视频| 国产一级毛片七仙女欲春2 | 日韩欧美一区二区三区在线观看| 国产精品电影一区二区三区| 国产黄色小视频在线观看| 麻豆成人av在线观看| 男人操女人黄网站| 亚洲欧美日韩无卡精品| 国产aⅴ精品一区二区三区波| 日韩三级视频一区二区三区| 99国产精品一区二区三区| 在线观看午夜福利视频| 日韩精品青青久久久久久| 人妻久久中文字幕网| 免费在线观看完整版高清| 一级毛片高清免费大全| 99国产精品一区二区蜜桃av| 999精品在线视频| 国产区一区二久久| 国产一卡二卡三卡精品| 国产亚洲精品综合一区在线观看 | 国产av一区在线观看免费| 黄色视频不卡| 欧美激情久久久久久爽电影| 国产亚洲av嫩草精品影院| 午夜福利在线在线| 高清在线国产一区| 午夜福利免费观看在线| 欧美精品亚洲一区二区| 国产精品美女特级片免费视频播放器 | 国产精品永久免费网站| 亚洲中文av在线| 色老头精品视频在线观看| 久久午夜亚洲精品久久| 亚洲国产毛片av蜜桃av| 久久久久精品国产欧美久久久| 久久久久久国产a免费观看| 岛国在线观看网站| 亚洲一卡2卡3卡4卡5卡精品中文| 国产成人系列免费观看| 久久久久国产一级毛片高清牌| 成人三级黄色视频| 十分钟在线观看高清视频www| av中文乱码字幕在线| 18禁黄网站禁片午夜丰满| 亚洲真实伦在线观看| 欧美一级毛片孕妇| 国产免费av片在线观看野外av| 亚洲在线自拍视频| www.999成人在线观看| 一区二区三区国产精品乱码| 免费观看人在逋| 欧美乱码精品一区二区三区| 黄色视频不卡| 亚洲 欧美一区二区三区| 亚洲,欧美精品.| 淫妇啪啪啪对白视频| 免费电影在线观看免费观看| 国产精品久久电影中文字幕| 国产精品香港三级国产av潘金莲| 久久伊人香网站| 欧美黄色片欧美黄色片| 黄色毛片三级朝国网站| 啦啦啦 在线观看视频| 窝窝影院91人妻| 国产三级黄色录像| 亚洲av成人一区二区三| 国产成人精品无人区| 黄色成人免费大全| 大型av网站在线播放| 又黄又爽又免费观看的视频| 大型av网站在线播放| 国产精品99久久99久久久不卡| 中国美女看黄片| 亚洲熟女毛片儿| 一本一本综合久久| ponron亚洲| 日韩一卡2卡3卡4卡2021年| 精品久久久久久,| 婷婷丁香在线五月| 久热这里只有精品99| 制服人妻中文乱码| 一二三四社区在线视频社区8| 看片在线看免费视频| 嫩草影院精品99| 免费在线观看视频国产中文字幕亚洲| а√天堂www在线а√下载| 国产黄片美女视频| a在线观看视频网站| 老汉色∧v一级毛片| 欧美色视频一区免费| 久久精品人妻少妇| 国产区一区二久久| 精品国产一区二区三区四区第35| 精品午夜福利视频在线观看一区| 日韩精品免费视频一区二区三区| 精品欧美一区二区三区在线| 91成人精品电影| 免费看日本二区| 亚洲第一av免费看| 特大巨黑吊av在线直播 | 91麻豆av在线| 亚洲精品粉嫩美女一区| 亚洲欧美精品综合一区二区三区| 欧美色欧美亚洲另类二区| 日韩欧美三级三区| www日本在线高清视频| 亚洲人成电影免费在线| 草草在线视频免费看| 国产成人精品久久二区二区免费| 久久久久九九精品影院| 97人妻精品一区二区三区麻豆 | 国产91精品成人一区二区三区| 国产伦一二天堂av在线观看| 麻豆国产av国片精品| 一a级毛片在线观看| 色综合婷婷激情| 在线观看免费日韩欧美大片| 美女午夜性视频免费| 亚洲中文日韩欧美视频| 一区二区三区精品91| 美女午夜性视频免费| 看黄色毛片网站| 在线观看免费日韩欧美大片| 中文字幕人成人乱码亚洲影| 麻豆一二三区av精品| 高清毛片免费观看视频网站| 国产在线观看jvid| 亚洲av中文字字幕乱码综合 | 欧美又色又爽又黄视频| 人人澡人人妻人| 亚洲aⅴ乱码一区二区在线播放 | 丝袜美腿诱惑在线| 每晚都被弄得嗷嗷叫到高潮| 日本三级黄在线观看| 香蕉丝袜av| 午夜福利高清视频| 最新美女视频免费是黄的| 欧美成人性av电影在线观看| 日韩精品青青久久久久久| 亚洲人成网站高清观看| 免费搜索国产男女视频| 最好的美女福利视频网| 亚洲 欧美一区二区三区| 在线观看午夜福利视频| 亚洲第一av免费看| 在线国产一区二区在线| 精品欧美国产一区二区三| 久久婷婷人人爽人人干人人爱| 成人永久免费在线观看视频| 可以在线观看的亚洲视频| 在线免费观看的www视频|