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

    關(guān)于易腐產(chǎn)品的M/M/1 生產(chǎn)服務(wù)庫存模型的最優(yōu)控制

    2021-01-09 02:44:20岳德權(quán)張媛媛田瑞玲
    工程數(shù)學(xué)學(xué)報 2020年6期
    關(guān)鍵詞:易腐穩(wěn)態(tài)庫存

    王 賽, 岳德權(quán), 張媛媛, 田瑞玲

    (燕山大學(xué)理學(xué)院,秦皇島 066004)

    1 引言

    隨著經(jīng)濟的迅速發(fā)展,企業(yè)面臨的市場競爭越來越激烈,生產(chǎn)制造企業(yè)需要盡可能的降低庫存管理成本才能有可觀的經(jīng)濟效益.對于生產(chǎn)制造企業(yè)來說,庫存不僅用來減少顧客需求,也可以用來調(diào)節(jié)生產(chǎn),使得生產(chǎn)能夠平穩(wěn)的進行.因此生產(chǎn)庫存模型及其庫存的控制管理一直是庫存理論的重要研究內(nèi)容.

    生產(chǎn)庫存策略是生產(chǎn)庫存系統(tǒng)的重要研究問題.施文武等[1]研究了一種多周期隨機需求的生產(chǎn)庫存控制系統(tǒng),給出了費用函數(shù),并設(shè)計算法求得最優(yōu)控制策略.Benjaafar 等[2]研究了帶有不耐煩顧客的生產(chǎn)庫存系統(tǒng),確定了最佳生產(chǎn)策略,分析了最優(yōu)庫存控制策略對于各項性能指標的影響.Krishnamoorthy 和Viswanath[3]研究了具有(s,S)策略的帶有隨機服務(wù)時間的生產(chǎn)服務(wù)庫存系統(tǒng),利用擬生滅過程方法得到了系統(tǒng)的穩(wěn)態(tài)分布和費用函數(shù).Baek 和Moon[4]研究了具有(s,S)策略的多個服務(wù)員的生產(chǎn)服務(wù)庫存系統(tǒng),給出了系統(tǒng)的性能指標和費用函數(shù),并通過數(shù)值例子確定了系統(tǒng)最優(yōu)控制策略.Axs¨ater[5]研究了具有(s,S)策略且容量有限的生產(chǎn)庫存系統(tǒng),假設(shè)系統(tǒng)生產(chǎn)時間服從伽馬分布,利用M/G/1 的相關(guān)理論,得到了最優(yōu)的生產(chǎn)策略.Nair 和Jose[6]研究了具有(s,S)策略,兩種服務(wù)模式和重試需求的生產(chǎn)系統(tǒng),求出了系統(tǒng)的性能指標,并建立了費用函數(shù).

    在上述關(guān)于生產(chǎn)庫存系統(tǒng)的研究中,都沒有考慮產(chǎn)品的易腐性.我國每年因為易腐性產(chǎn)品造成的損失之和高達千億[7],所以易腐庫存系統(tǒng)的分析是很重要的.易腐性產(chǎn)品是指那些必須在有限時間內(nèi)售出,否則將發(fā)生變質(zhì)、損壞、揮發(fā)、過期且必須進行清倉處理的商品,如生鮮食品、水果、蔬菜、海鮮等.其顯著特點是在儲存和流通的過程中產(chǎn)品的數(shù)量會因腐爛變質(zhì)、揮發(fā)、失效等逐漸減少.

    易腐品的庫存問題已經(jīng)引起了學(xué)者的廣泛關(guān)注.Kouki 等[8]研究了具有多種類型產(chǎn)品的易腐庫存問題,并且給出了系統(tǒng)的最優(yōu)庫存控制策略.Sivakumar[9]研究了具有(s,S)策略,帶有有限個顧客的易腐品庫存系統(tǒng),根據(jù)穩(wěn)態(tài)方程求出了穩(wěn)態(tài)概率向量,并建立了費用函數(shù).Vijaya Laxmi 和Soujanya[10]分析了帶有服務(wù)中斷、重試需求和負顧客的易腐品庫存模型,并利用矩陣分析法求出了模型在穩(wěn)態(tài)下的聯(lián)合概率分布,討論了系統(tǒng)的各項性能指標,確定了費用函數(shù)并得到了庫存的最優(yōu)控制方案.Shophia Lawrence 等[11]研究了顧客有限的易腐產(chǎn)品的庫存系統(tǒng).假設(shè)服務(wù)時間服從位相分布,產(chǎn)品壽命服從指數(shù)分布,利用拉普拉斯變換方法推導(dǎo)出了系統(tǒng)的性能指標,構(gòu)建了費用函數(shù).Manuel 等[12]研究了具有(s,S)策略的兩類客戶的易腐產(chǎn)品庫存系統(tǒng),求得了穩(wěn)態(tài)概率向量,并給出了費用函數(shù)的表達式.

    在上述關(guān)于易腐產(chǎn)品的庫存控制文獻中,都是關(guān)注庫存系統(tǒng)的性能分析和庫存策略問題.目前在生產(chǎn)服務(wù)庫存系統(tǒng)中,關(guān)于易腐產(chǎn)品的生產(chǎn)庫存策略的研究工作較少.Sangeetha 等[13]研究了帶有重試需求的易腐品生產(chǎn)庫存系統(tǒng).假設(shè)系統(tǒng)的顧客數(shù)量有限,基于(s,S)策略建立半馬爾可夫過程,利用線性規(guī)劃方法求出系統(tǒng)的最優(yōu)生產(chǎn)速率,并給出數(shù)值例子進行分析.本文討論了易腐產(chǎn)品的M/M/1 生產(chǎn)服務(wù)庫存模型.假設(shè)顧客數(shù)量是無限的,規(guī)定當系統(tǒng)庫存為零時,允許顧客進入系統(tǒng).根據(jù)擬生滅過程理論,我們求得了系統(tǒng)的平衡條件,進一步得到了穩(wěn)態(tài)概率向量矩陣幾何解,性能指標及成本函數(shù).最后給出數(shù)值例子,利用遺傳算法實現(xiàn)最優(yōu)解的有效搜索,得到了系統(tǒng)的最優(yōu)控制策略.

    2 模型描述

    我們考慮具有(s,S)生產(chǎn)策略的易腐產(chǎn)品的馬爾可夫生產(chǎn)服務(wù)庫存模型,可描述為:

    顧客需求的到達服從參數(shù)為λ 的泊松過程,顧客到達系統(tǒng)后形成一個隊列,每位顧客的庫存需求量為一個單位庫存.系統(tǒng)中只有一個服務(wù)員,服務(wù)過程需要一定的時間,系統(tǒng)服務(wù)時間服從參數(shù)為μ的指數(shù)分布,服務(wù)規(guī)則為FCFS.

    系統(tǒng)采用(s,S)生產(chǎn)策略,即當系統(tǒng)的庫存水平下降到s 時,生產(chǎn)系統(tǒng)立即啟動,并且每次只生產(chǎn)一個,直至庫存水平達到最大值S 時,才停止生產(chǎn),其中s < S.生產(chǎn)時間服從參數(shù)為η 的指數(shù)分布.

    產(chǎn)品壽命服從參數(shù)為γ 的指數(shù)分布,變質(zhì)后的產(chǎn)品將不再具有價值.系統(tǒng)是無損失的,即當系統(tǒng)庫存為零時,顧客可以進入系統(tǒng)等待.假設(shè)需求到達過程,服務(wù)過程和生產(chǎn)過程是相互獨立的.

    3 系統(tǒng)平衡條件

    3.1 狀態(tài)過程

    定義系統(tǒng)的狀態(tài)過程為ψ = {(N(t),I(t),C(t)); t ≥0},其中N(t)表示時刻t 系統(tǒng)中的顧客數(shù),I(t)表示時刻t 系統(tǒng)的庫存水平,C(t)表示時刻t 系統(tǒng)的生產(chǎn)狀態(tài),其中C(t) = 0 表示生產(chǎn)系統(tǒng)關(guān)閉,C(t) = 1 表示生產(chǎn)系統(tǒng)開啟.根據(jù)(s,S)策略,C(t)應(yīng)滿足如下關(guān)系

    過程ψ ={(N(t),I(t),C(t)); t ≥0}的狀態(tài)空間為

    將狀態(tài)按字典序排列可得過程ψ ={(N(t),I(t),C(t)); t ≥0}的無窮小生成元如下

    其中

    其中

    D1是(2S ?s)×S 維矩陣,D2是(2S ?s)×(S ?s)維矩陣,C 是(2S ?s)×(2S ?s)維矩陣,I 是(2S ?s)維單位陣,I1是(2S ?s ?1)維單位陣.

    3.2 系統(tǒng)平衡條件

    由矩陣Q 的結(jié)構(gòu)可知過程ψ = {(N(t),I(t),C(t)); t ≥0}是擬生滅過程[14].令H =C+B+A,則有

    H =(H1H2),

    其中

    L1(i)=?(η+μ+iγ), 1 ≤i ≤S ?1,

    L2(i)=?(μ+iγ), s+1 ≤i ≤S, G(i)=iγ+μ, 1 ≤i ≤S,

    H1是(2S ?s)×S 維矩陣,H2是(2S ?s)×(S ?s)維矩陣.

    定義H 的平穩(wěn)概率向量為

    π =(π(0,1),π(1,1),··· ,π(S ?1,1),π(s+1,0),π(s+2,0),··· ,π(S,0)),

    平穩(wěn)概率向量π 滿足如下條件

    其中e1是所有元素為1 的(2S ?s)維列向量.通過計算整理,概率向量π 的各個分量可表示為

    其中

    由文獻[14]知,過程ψ ={(N(t),I(t),C(t)); t ≥0}正常返的充分必要條件為πAe1<πCe1.經(jīng)計算整理可得

    通過計算,πAe1<πCe1等價于

    所以,公式(1)就是系統(tǒng)穩(wěn)態(tài)平衡的充分必要條件.

    注1 根據(jù)給出的πAe1和πCe1的表達式,可以看出πAe1表示系統(tǒng)的到達率,πCe1表示當庫存非零時系統(tǒng)的服務(wù)率,即系統(tǒng)的有效服務(wù)率.根據(jù)排隊論的知識,當系統(tǒng)的到達率λ 小于系統(tǒng)的有效服務(wù)率時,系統(tǒng)達到平衡狀態(tài).

    4 矩陣幾何解

    本節(jié)利用擬生滅過程理論,首先給出系統(tǒng)的穩(wěn)態(tài)概率分布的矩陣幾何解,然后討論率陣R 和邊界狀態(tài)概率的計算問題.定義穩(wěn)態(tài)概率為

    在穩(wěn)態(tài)條件ρ < 1 下,Q 的平穩(wěn)概率向量P 存在.相應(yīng)于Q 的分塊結(jié)構(gòu),穩(wěn)態(tài)概率向量P 分塊如下

    P =(P0,P1,··· ,Pi,···),

    其中

    Pi={P(i,0,1),P(i,1,1),··· ,P(i,S ?1,1),P(i,s+1,0),··· ,P(i,S,0)}, i ≥0.

    穩(wěn)態(tài)向量P 滿足如下平衡方程

    其中e 是元素都為1 的適當維數(shù)的列向量.根據(jù)公式(2),我們可以得到如下方程組

    由文獻[15]知,系統(tǒng)的穩(wěn)態(tài)概率向量有矩陣幾何解

    其中R 是矩陣二次方程

    的最小非負解,其譜半徑sp(R)<1.P0是方程組的唯一正解,其中0 矩陣是元素都是0 的(2S ?s)維行向量,I 是(2S ?s)維單位矩陣.

    為了計算系統(tǒng)的平穩(wěn)向量,我們需要求出方程(7)的最小非負解R.在計算R 時,我們采用改進的循環(huán)簡約算法.循環(huán)簡約算法在文獻[16,17]中有詳細的介紹,循環(huán)簡約算法步驟如下:

    步驟1 令

    步驟2 矩陣^B(j)和率陣R 由以下循環(huán)迭代方法求得

    直至迭代得到R 值的差的范數(shù)滿足‖R(j+1)?R(j)‖≤ε,從而求得率陣R;

    步驟3 由方程(8)及Pi=P0Ri, i ≥0,即可求得平穩(wěn)分布.

    注2 為計算系統(tǒng)的平穩(wěn)概率向量,我們需要求解矩陣二次方程(7)的最小非負解.在本文中,矩陣B、矩陣C 和矩陣A 均不是特殊矩陣,若直接利用方程(7)求解率陣R 則計算過程過于復(fù)雜,不易直接求出.例如,當采用(1,2)策略時,為求解R 需要求解一個九元二次方程組,求解過程相當復(fù)雜.因此我們使用上述算法,可近似求出率陣R 的數(shù)值計算結(jié)果.

    例1 給定系統(tǒng)的參數(shù)值s = 1, S = 2, λ = 3.5, η = 5, γ = 0.05, μ = 8,根據(jù)率陣R 的計算方法,可近似求出R 如下

    利用方程組(8)求解向量P0,結(jié)果如下

    P0=(0.04132849666420,0.07048216383424,0.09789189421422),

    再由公式(6)可以計算出系統(tǒng)的穩(wěn)態(tài)概率向量.

    5 穩(wěn)態(tài)性能指標及成本分析

    5.1 穩(wěn)態(tài)指標

    根據(jù)上述穩(wěn)態(tài)概率向量的表達式,易得系統(tǒng)穩(wěn)態(tài)性能指標,具體如下:

    1) 平均隊長

    其中δ1是元素為1 的(2S ?s)維列向量.

    2) 平均等待隊長

    3) 平均庫存

    其中δ2為(2S ?s)維列向量,且δ2=(0,1,2,··· ,S ?1,s+1,s+2,··· ,S)T.

    4) 平均易腐率

    5) 平均生產(chǎn)率

    其中δ3=(a b)T,a 是元素都為1 的S 維行向量,b 是元素都為0 的(S ?s)維行向量.

    6) 平均生產(chǎn)啟動率

    其中δ4為第(S+1)個元素為1,其余元素為0 的(2S ?s)維列向量.

    5.2 成本函數(shù)分析

    假設(shè)該庫存系統(tǒng)的成本主要由庫存保管成本、生產(chǎn)成本、設(shè)備啟動所需要的固定成本、產(chǎn)品腐爛的損失成本以及系統(tǒng)中每個顧客的等待成本組成.令單位時間單位庫存的保管成本為Cinv,單位時間的生產(chǎn)成本為Crp,每單位產(chǎn)品腐爛的成本損失為Cp,系統(tǒng)中每個顧客的等待成本Clq,每次重啟設(shè)備所需要的固定成本為K,則系統(tǒng)的成本函數(shù)C(s,S)為

    顯然,系統(tǒng)的成本函數(shù)是關(guān)于庫存控制決策變量s, S 的非線性函數(shù),且決策變量是離散的整數(shù)型變量.因為成本函數(shù)的表達式具有高復(fù)雜性,所以對于最小值的計算有一定難度.對于具體的數(shù)值例子,本文采用遺傳算法求解最優(yōu)值.

    5.3 遺傳算法

    遺傳算法(GA)是模擬自然界遺傳、變異、適者生存的進化思想來求解優(yōu)化問題,對目標函數(shù)和約束條件的要求較低.對于存在有限的可行解空間的純整數(shù)規(guī)劃問題,可采用適當?shù)倪z傳算法得到最優(yōu)解,本文采用了豐建榮等[18]提出的利用二進制編碼的遺傳算法,具體步驟如下:

    步驟1 隨機產(chǎn)生初始種群,個體數(shù)目一定.如果個體無效,那么重新隨機生成初始個體,直至有效;

    步驟2 評價群體的適應(yīng)度,找到最好的染色體;

    步驟3 用輪盤賭策略對每一代種群中的染色體進行選擇;

    步驟4 進行個體交叉操作,判斷個體的有效性.如果是無效個體,則隨機生成一個交叉位置進行交叉,直至有效;

    步驟5 進行個體變異操作.如果產(chǎn)生無效個體,則重新進行變異操作,直至有效個體產(chǎn)生;

    步驟6 由選擇、交叉和變異產(chǎn)生新一代種群;

    步驟7 對新種群適應(yīng)度評價,找到最好的染色體,將它與上一次進化中最好的染色體比較,記錄每一代進化中最好的適應(yīng)度和平均適應(yīng)度;

    步驟8 如果不滿足算法終止的條件,轉(zhuǎn)到步驟3.否則,輸出當前最優(yōu)個體,算法結(jié)束.

    6 數(shù)值分析

    本節(jié)通過數(shù)值例子來考察系統(tǒng)參數(shù)的敏感性.根據(jù)上一節(jié)設(shè)置的遺傳算法,設(shè)置種群大小為50,交叉率為0.6,變異率為0.01,算法終止條件為100 次(通過反復(fù)實驗,程序迭代100 次結(jié)束,穩(wěn)定收斂到最優(yōu)解),并用Matlab R2014a 編程進行數(shù)值實驗.

    表1 給出了參數(shù)λ 取不同的值時,系統(tǒng)的最優(yōu)控制策略和一些性能指標的取值變化,其中參數(shù)設(shè)置為η = 5, γ = 0.05, μ = 8, Cinv= 10, Crp= 20, Cp= 50, Clq=200, K =2500.

    表1 參數(shù)λ 的敏感性分析

    由表1 可見,隨著參數(shù)λ 的增大,系統(tǒng)的安全庫存水平和最大庫存水平都呈上升趨勢.系統(tǒng)的平均庫存、平均易腐率、平均生產(chǎn)速率、平均等待隊長逐漸增加,平均生產(chǎn)啟動率逐漸減小,最終成本函數(shù)逐漸增加.對于企業(yè)來說,在系統(tǒng)成本參數(shù)確定的條件下,市場需求增加時,庫存控制的最優(yōu)策略應(yīng)相應(yīng)的上調(diào)安全庫存水平和最大庫存水平.

    表2 給出了參數(shù)η 取不同的值時,系統(tǒng)的最優(yōu)控制策略和一些性能指標的取值變化,其中參數(shù)為λ=3.5, γ =0.05, μ=8, Cinv=10, Crp=20, Cp=50, Clq=200, K =2500.

    表2 參數(shù)η 的敏感性分析

    由表2 可見,隨著參數(shù)η 的增加,系統(tǒng)的安全庫存水平,最大庫存水平均呈下降趨勢.系統(tǒng)的平均庫存、平均生產(chǎn)率、平均易腐率先減后增,平均等待隊長逐漸減小,平均生產(chǎn)啟動率逐漸增大,最終導(dǎo)致系統(tǒng)的成本逐漸增大.對于企業(yè)來說,可以通過技術(shù)改進,更新設(shè)備等方法提高生產(chǎn)率.同時當生產(chǎn)效率提高時,庫存控制的最優(yōu)策略需要下調(diào)安全庫存水平和最大庫存水平,企業(yè)在調(diào)整時還需要注意生產(chǎn)效率過高會造成產(chǎn)品堆積.

    表3 給出了參數(shù)γ 取不同的值時,系統(tǒng)的最優(yōu)控制策略和一些性能指標的取值變化,參數(shù)為λ=3.5, η =5, μ=8, Cinv=10, Crp=20, Cp=50, Clq=200, K =2500.

    表3 參數(shù)γ 的敏感性分析

    由表3 可見,隨著γ 的增加,系統(tǒng)最大庫存水平逐漸減少.系統(tǒng)的平均庫存、平均生產(chǎn)啟動率逐漸減少,平均生產(chǎn)率、產(chǎn)品的平均易腐率逐漸增加,平均等待隊長先增后減,最終導(dǎo)致系統(tǒng)費用增加.易腐產(chǎn)品的壽命參數(shù)γ 越小,產(chǎn)品的腐爛速度則越慢,系統(tǒng)成本越小.對于企業(yè)來說,可優(yōu)化易腐品的保存方法,延長產(chǎn)品壽命,進而降低庫存成本.

    表4 給出了參數(shù)μ取不同的值時,系統(tǒng)的最優(yōu)控制策略和一些性能指標的取值變化,參數(shù)設(shè)置為λ=3.5, η =5, γ =0.05, Cinv=10, Crp=20, Cp=50, Clq=200, K =2500.

    表4 參數(shù)μ的敏感性分析

    由表4 可見,隨著參數(shù)μ的增加,系統(tǒng)的安全庫存和最大庫存水平保持穩(wěn)定.系統(tǒng)平均生產(chǎn)啟動率逐漸增加,但增幅很小.平均生產(chǎn)率、平均庫存和平均易腐率、平均等待隊長逐漸減小,系統(tǒng)成本逐漸減小.在系統(tǒng)成本參數(shù)確定的條件下,服務(wù)速率越高,系統(tǒng)的成本則越低.所以對于企業(yè)而言,應(yīng)該提高員工的服務(wù)效率,進而降低系統(tǒng)的庫存成本.

    7 結(jié)論

    本文研究了有關(guān)易腐產(chǎn)品的M/M/1 生產(chǎn)服務(wù)庫存模型,建立了關(guān)于系統(tǒng)中的顧客數(shù)、庫存水平和生產(chǎn)狀態(tài)三個隨機變量的三維擬生滅過程.利用擬生滅過程理論,求出了系統(tǒng)的穩(wěn)態(tài)平衡條件,給出了系統(tǒng)的穩(wěn)態(tài)概率的矩陣幾何解,并求出一些系統(tǒng)穩(wěn)態(tài)下的性能指標.進一步構(gòu)建了成本函數(shù),針對系統(tǒng)成本函數(shù)的非線性、整數(shù)型變量的特點,運用遺傳算法,實現(xiàn)了最優(yōu)解的有效搜索.通過數(shù)值實驗分析了系統(tǒng)參數(shù)的敏感性,這些分析結(jié)果,有助于生產(chǎn)型企業(yè)在外部參數(shù)發(fā)生變化時,根據(jù)實際情況及時調(diào)節(jié)庫存控制策略,對庫存管理實踐具有指導(dǎo)借鑒意義.

    猜你喜歡
    易腐穩(wěn)態(tài)庫存
    易腐果蔬動態(tài)保質(zhì)期評估和庫存管理策略探討
    ——基于集成射頻識別技術(shù)
    可變速抽水蓄能機組穩(wěn)態(tài)運行特性研究
    阿U漫說垃圾分類
    智慧少年(2022年2期)2022-06-23 15:03:57
    易腐垃圾處理技術(shù)及其效果研究進展
    碳化硅復(fù)合包殼穩(wěn)態(tài)應(yīng)力與失效概率分析
    電廠熱力系統(tǒng)穩(wěn)態(tài)仿真軟件開發(fā)
    煤氣與熱力(2021年4期)2021-06-09 06:16:54
    元中期歷史劇對社會穩(wěn)態(tài)的皈依與維護
    中華戲曲(2020年1期)2020-02-12 02:28:18
    家庭易腐垃圾處理現(xiàn)狀分析與建議
    一二線城市庫存減少5.2%
    營銷4C與房產(chǎn)去庫存
    中國市場(2016年45期)2016-05-17 05:15:23
    亚洲专区国产一区二区| 国产亚洲av嫩草精品影院| 久久久久性生活片| 真人做人爱边吃奶动态| 好男人电影高清在线观看| 久久中文看片网| 日本免费a在线| 熟女人妻精品中文字幕| 国产探花在线观看一区二区| 日本黄色视频三级网站网址| 丰满人妻一区二区三区视频av | 久久久精品欧美日韩精品| 精品欧美国产一区二区三| 高清在线国产一区| 精品不卡国产一区二区三区| 国产高清激情床上av| 性色av乱码一区二区三区2| 三级毛片av免费| 亚洲欧美日韩高清在线视频| 亚洲成人久久性| 18禁国产床啪视频网站| 久久性视频一级片| 精品国产乱子伦一区二区三区| 变态另类成人亚洲欧美熟女| 国产高清视频在线观看网站| 在线观看免费午夜福利视频| 亚洲中文日韩欧美视频| 免费在线观看影片大全网站| 久久久水蜜桃国产精品网| 俺也久久电影网| 精品99又大又爽又粗少妇毛片 | 男女视频在线观看网站免费| 国产97色在线日韩免费| 国产亚洲精品一区二区www| 成人国产一区最新在线观看| 变态另类丝袜制服| 久久精品91无色码中文字幕| 99久久成人亚洲精品观看| 757午夜福利合集在线观看| 2021天堂中文幕一二区在线观| 女生性感内裤真人,穿戴方法视频| 黄片大片在线免费观看| 久久久久久久久免费视频了| 亚洲一区二区三区色噜噜| 丰满人妻一区二区三区视频av | 国产麻豆成人av免费视频| 中文字幕久久专区| 此物有八面人人有两片| 国产一区二区激情短视频| 久9热在线精品视频| 最新美女视频免费是黄的| 三级男女做爰猛烈吃奶摸视频| 一区福利在线观看| 丝袜人妻中文字幕| 日本一本二区三区精品| ponron亚洲| 19禁男女啪啪无遮挡网站| 后天国语完整版免费观看| 色尼玛亚洲综合影院| 白带黄色成豆腐渣| 久久精品国产清高在天天线| 久久午夜亚洲精品久久| 国产精品野战在线观看| 亚洲精品一区av在线观看| 欧美性猛交╳xxx乱大交人| 亚洲人成伊人成综合网2020| 老司机深夜福利视频在线观看| 久久人人精品亚洲av| 国产精品精品国产色婷婷| 精品无人区乱码1区二区| 欧美乱码精品一区二区三区| 日韩 欧美 亚洲 中文字幕| av视频在线观看入口| 婷婷丁香在线五月| 国语自产精品视频在线第100页| 午夜a级毛片| 亚洲黑人精品在线| 两个人看的免费小视频| 大型黄色视频在线免费观看| 久久中文字幕一级| 欧美中文综合在线视频| 国产极品精品免费视频能看的| 久久久久久久久久黄片| 国内精品久久久久久久电影| 久久久久国产精品人妻aⅴ院| 中文字幕人成人乱码亚洲影| a在线观看视频网站| 观看免费一级毛片| 人人妻,人人澡人人爽秒播| 噜噜噜噜噜久久久久久91| 在线看三级毛片| 国产午夜精品论理片| 18禁黄网站禁片免费观看直播| 狠狠狠狠99中文字幕| 99热这里只有精品一区 | av在线天堂中文字幕| 国产69精品久久久久777片 | 精品一区二区三区视频在线 | 成人国产一区最新在线观看| 久久久久久久精品吃奶| 人妻丰满熟妇av一区二区三区| 每晚都被弄得嗷嗷叫到高潮| 亚洲av第一区精品v没综合| 特大巨黑吊av在线直播| 一区二区三区激情视频| 九九热线精品视视频播放| 美女高潮的动态| 国产免费男女视频| 亚洲激情在线av| 美女cb高潮喷水在线观看 | 99热只有精品国产| 久久久久久国产a免费观看| 99国产精品一区二区蜜桃av| 高清在线国产一区| 亚洲av成人不卡在线观看播放网| 国产成人精品无人区| 欧洲精品卡2卡3卡4卡5卡区| 悠悠久久av| 国产精品爽爽va在线观看网站| 日本撒尿小便嘘嘘汇集6| 特级一级黄色大片| 天堂动漫精品| 首页视频小说图片口味搜索| 亚洲色图av天堂| 欧美午夜高清在线| 三级毛片av免费| 搞女人的毛片| 97超视频在线观看视频| 亚洲欧美激情综合另类| 老熟妇乱子伦视频在线观看| 亚洲精品美女久久av网站| 中文字幕精品亚洲无线码一区| 久久久久免费精品人妻一区二区| 国产精品久久久人人做人人爽| 嫁个100分男人电影在线观看| 亚洲在线观看片| 欧美又色又爽又黄视频| 中亚洲国语对白在线视频| 国产 一区 欧美 日韩| 我的老师免费观看完整版| 成人特级av手机在线观看| 国产亚洲欧美98| 757午夜福利合集在线观看| 18禁国产床啪视频网站| 国产精品亚洲av一区麻豆| 免费av不卡在线播放| 精品乱码久久久久久99久播| 成年免费大片在线观看| 亚洲人成网站在线播放欧美日韩| 国产精品爽爽va在线观看网站| 色综合欧美亚洲国产小说| 曰老女人黄片| 成年女人永久免费观看视频| 法律面前人人平等表现在哪些方面| 熟女电影av网| 精品久久久久久久毛片微露脸| 亚洲精品456在线播放app | 成人av在线播放网站| 亚洲午夜理论影院| 深夜精品福利| 成年女人毛片免费观看观看9| 国产黄片美女视频| 亚洲av日韩精品久久久久久密| 久久精品国产99精品国产亚洲性色| 精品福利观看| 九九热线精品视视频播放| 老司机福利观看| 亚洲精品在线美女| 日韩av在线大香蕉| 国产亚洲精品av在线| 国产不卡一卡二| 好男人在线观看高清免费视频| 91九色精品人成在线观看| avwww免费| 国产精华一区二区三区| 亚洲自偷自拍图片 自拍| 久久精品国产清高在天天线| 成年女人毛片免费观看观看9| 亚洲成人免费电影在线观看| 欧美日韩乱码在线| 国产探花在线观看一区二区| 五月伊人婷婷丁香| 亚洲成人久久爱视频| 欧美绝顶高潮抽搐喷水| 成年女人永久免费观看视频| 级片在线观看| 亚洲自拍偷在线| 我要搜黄色片| 首页视频小说图片口味搜索| 精品国产超薄肉色丝袜足j| 精品乱码久久久久久99久播| 亚洲欧美日韩高清专用| 又黄又爽又免费观看的视频| 男女床上黄色一级片免费看| 久久久久国产一级毛片高清牌| 午夜两性在线视频| 91av网站免费观看| 脱女人内裤的视频| 国产欧美日韩一区二区三| 在线观看免费午夜福利视频| 小说图片视频综合网站| 国内毛片毛片毛片毛片毛片| 淫妇啪啪啪对白视频| 国产av在哪里看| 一个人观看的视频www高清免费观看 | 国语自产精品视频在线第100页| 亚洲精品粉嫩美女一区| 国产精品 欧美亚洲| 好男人电影高清在线观看| 夜夜爽天天搞| 在线观看舔阴道视频| 他把我摸到了高潮在线观看| 免费看日本二区| 精品国产亚洲在线| 嫩草影院入口| 黑人操中国人逼视频| 1024香蕉在线观看| 欧美日韩亚洲国产一区二区在线观看| 国产男靠女视频免费网站| 欧美xxxx黑人xx丫x性爽| 99国产精品99久久久久| 亚洲,欧美精品.| 国产成人欧美在线观看| 亚洲专区中文字幕在线| 久久午夜亚洲精品久久| 在线a可以看的网站| 天天添夜夜摸| 精品国内亚洲2022精品成人| 他把我摸到了高潮在线观看| 精品午夜福利视频在线观看一区| 99国产极品粉嫩在线观看| 午夜福利在线观看吧| 国产亚洲精品av在线| 美女高潮喷水抽搐中文字幕| 亚洲avbb在线观看| 国产av在哪里看| 久久九九热精品免费| 亚洲专区字幕在线| 亚洲无线在线观看| 欧美一级a爱片免费观看看| 亚洲国产精品999在线| 俺也久久电影网| 老司机深夜福利视频在线观看| 国产av不卡久久| 校园春色视频在线观看| 国产男靠女视频免费网站| 男人和女人高潮做爰伦理| 欧美乱色亚洲激情| 日韩欧美精品v在线| 免费高清视频大片| 国产极品精品免费视频能看的| 国产精品一及| 很黄的视频免费| 日日摸夜夜添夜夜添小说| 美女午夜性视频免费| 亚洲黑人精品在线| 级片在线观看| 国产亚洲精品久久久久久毛片| 又紧又爽又黄一区二区| 18禁国产床啪视频网站| xxxwww97欧美| 一个人看视频在线观看www免费 | 午夜福利在线观看免费完整高清在 | 一夜夜www| 在线看三级毛片| 少妇熟女aⅴ在线视频| 啦啦啦免费观看视频1| a在线观看视频网站| 日本一本二区三区精品| 999精品在线视频| 日本精品一区二区三区蜜桃| 亚洲成人久久爱视频| 91九色精品人成在线观看| 啦啦啦免费观看视频1| 桃红色精品国产亚洲av| 欧美最黄视频在线播放免费| av天堂在线播放| 69av精品久久久久久| 最新美女视频免费是黄的| 成人av一区二区三区在线看| 精品久久久久久久毛片微露脸| 精品午夜福利视频在线观看一区| 黄色成人免费大全| АⅤ资源中文在线天堂| 99久久精品国产亚洲精品| 久久久久国产精品人妻aⅴ院| 日韩欧美一区二区三区在线观看| 久久这里只有精品19| 我的老师免费观看完整版| 一边摸一边抽搐一进一小说| 成年版毛片免费区| 岛国在线免费视频观看| 国内精品一区二区在线观看| 91在线精品国自产拍蜜月 | 免费在线观看影片大全网站| 男人和女人高潮做爰伦理| 三级男女做爰猛烈吃奶摸视频| 亚洲欧美一区二区三区黑人| 一本综合久久免费| 亚洲av日韩精品久久久久久密| 久久国产精品人妻蜜桃| 一进一出好大好爽视频| av黄色大香蕉| 99久久精品热视频| 欧美国产日韩亚洲一区| 成人三级黄色视频| 性欧美人与动物交配| 国产乱人视频| 欧美日韩福利视频一区二区| 免费在线观看影片大全网站| 久久久国产成人免费| 亚洲中文字幕日韩| 性色avwww在线观看| 国产毛片a区久久久久| av在线天堂中文字幕| 波多野结衣巨乳人妻| 天堂√8在线中文| 中文字幕最新亚洲高清| 国产欧美日韩一区二区三| 国产精品久久电影中文字幕| 国产精品av视频在线免费观看| 精品日产1卡2卡| 午夜福利在线在线| 精品一区二区三区四区五区乱码| 久久精品人妻少妇| 99久久国产精品久久久| 日韩人妻高清精品专区| 成人精品一区二区免费| 亚洲av成人av| 美女被艹到高潮喷水动态| 十八禁网站免费在线| 国产人伦9x9x在线观看| 国产欧美日韩一区二区三| 国产精品 欧美亚洲| 男女下面进入的视频免费午夜| 欧美3d第一页| 久久国产乱子伦精品免费另类| 亚洲人与动物交配视频| 好男人电影高清在线观看| 亚洲精品久久国产高清桃花| 国语自产精品视频在线第100页| 国产成人精品久久二区二区91| 99国产精品一区二区三区| 听说在线观看完整版免费高清| 国产成人精品久久二区二区免费| 久久久久久国产a免费观看| 窝窝影院91人妻| 哪里可以看免费的av片| 啦啦啦观看免费观看视频高清| 久久久久久久久久黄片| 国产探花在线观看一区二区| 久久久久久国产a免费观看| 可以在线观看毛片的网站| 五月玫瑰六月丁香| 一级毛片精品| 成人特级av手机在线观看| 麻豆一二三区av精品| 国产精品日韩av在线免费观看| 脱女人内裤的视频| 亚洲无线在线观看| 国产午夜精品论理片| 国产精品女同一区二区软件 | 国产精品久久久久久人妻精品电影| 国产免费男女视频| 欧美在线一区亚洲| 婷婷亚洲欧美| 嫩草影院入口| 一区二区三区激情视频| 午夜精品久久久久久毛片777| 免费在线观看影片大全网站| 一进一出抽搐动态| 久久久久亚洲av毛片大全| 国产成人av教育| 日韩人妻高清精品专区| e午夜精品久久久久久久| 精品人妻1区二区| e午夜精品久久久久久久| 久久人人精品亚洲av| 国产成人欧美在线观看| 一区福利在线观看| 国产综合懂色| 日韩三级视频一区二区三区| 麻豆一二三区av精品| 午夜福利免费观看在线| 久久精品亚洲精品国产色婷小说| 99国产精品一区二区蜜桃av| 黑人巨大精品欧美一区二区mp4| 性色avwww在线观看| 久久精品91蜜桃| 毛片女人毛片| 又粗又爽又猛毛片免费看| 亚洲人成伊人成综合网2020| 97人妻精品一区二区三区麻豆| 婷婷丁香在线五月| 国产精品久久久久久精品电影| 中文字幕最新亚洲高清| 天堂√8在线中文| 看片在线看免费视频| 成年女人永久免费观看视频| 欧美高清成人免费视频www| 久久中文字幕人妻熟女| 国产伦一二天堂av在线观看| 岛国在线免费视频观看| 精品不卡国产一区二区三区| 精品一区二区三区视频在线 | 国产精品av久久久久免费| 特级一级黄色大片| 久久亚洲真实| 亚洲中文字幕日韩| 可以在线观看的亚洲视频| 久久久久亚洲av毛片大全| 成人鲁丝片一二三区免费| 美女午夜性视频免费| 亚洲精品乱码久久久v下载方式 | 最近在线观看免费完整版| 18禁观看日本| 国产亚洲欧美98| 欧美在线黄色| 欧美国产日韩亚洲一区| 免费高清视频大片| 91九色精品人成在线观看| 91在线精品国自产拍蜜月 | 国产91精品成人一区二区三区| 制服人妻中文乱码| 亚洲中文av在线| 午夜日韩欧美国产| 亚洲av熟女| 精品不卡国产一区二区三区| 欧美日韩一级在线毛片| 亚洲真实伦在线观看| 国产精品亚洲av一区麻豆| 男女下面进入的视频免费午夜| 亚洲欧洲精品一区二区精品久久久| 国产高清视频在线观看网站| netflix在线观看网站| 国内毛片毛片毛片毛片毛片| 一个人看的www免费观看视频| 久久久久国内视频| 18禁观看日本| 婷婷丁香在线五月| 亚洲,欧美精品.| 丁香六月欧美| 91av网站免费观看| 精品不卡国产一区二区三区| 欧美日韩一级在线毛片| av中文乱码字幕在线| 成年免费大片在线观看| 九九久久精品国产亚洲av麻豆 | 一a级毛片在线观看| 午夜久久久久精精品| 国产黄a三级三级三级人| 国产熟女xx| 亚洲黑人精品在线| 999精品在线视频| bbb黄色大片| 国产成人系列免费观看| 制服人妻中文乱码| a级毛片a级免费在线| 亚洲自偷自拍图片 自拍| 午夜福利在线观看吧| 1000部很黄的大片| 首页视频小说图片口味搜索| 亚洲激情在线av| 国产精品野战在线观看| 国产高清有码在线观看视频| 久久香蕉精品热| 可以在线观看的亚洲视频| 亚洲欧美日韩无卡精品| 制服人妻中文乱码| 亚洲av电影在线进入| 色老头精品视频在线观看| 91久久精品国产一区二区成人 | 女生性感内裤真人,穿戴方法视频| 国产成人精品无人区| 18禁黄网站禁片午夜丰满| 亚洲电影在线观看av| 免费看a级黄色片| 国内毛片毛片毛片毛片毛片| 色尼玛亚洲综合影院| 少妇人妻一区二区三区视频| 国产精华一区二区三区| 老熟妇乱子伦视频在线观看| 听说在线观看完整版免费高清| 99精品在免费线老司机午夜| 免费电影在线观看免费观看| 人妻久久中文字幕网| 在线视频色国产色| 亚洲人成电影免费在线| 亚洲中文字幕日韩| 欧美一区二区国产精品久久精品| 人妻丰满熟妇av一区二区三区| 美女高潮喷水抽搐中文字幕| 国产av一区在线观看免费| av国产免费在线观看| 老司机午夜十八禁免费视频| 三级国产精品欧美在线观看 | 狂野欧美激情性xxxx| 好男人电影高清在线观看| 免费观看精品视频网站| 又粗又爽又猛毛片免费看| 桃红色精品国产亚洲av| 网址你懂的国产日韩在线| 久久精品综合一区二区三区| 高潮久久久久久久久久久不卡| 在线免费观看的www视频| 欧美另类亚洲清纯唯美| 亚洲自拍偷在线| 欧美色欧美亚洲另类二区| 亚洲中文字幕日韩| 国产精品爽爽va在线观看网站| 99久久精品一区二区三区| 深夜精品福利| 午夜久久久久精精品| 九九久久精品国产亚洲av麻豆 | 精品国产乱子伦一区二区三区| 在线永久观看黄色视频| 国产v大片淫在线免费观看| 此物有八面人人有两片| 一本一本综合久久| 国产精品一区二区精品视频观看| 亚洲av熟女| 变态另类成人亚洲欧美熟女| 久久久国产成人免费| 国产av在哪里看| 亚洲男人的天堂狠狠| 国产私拍福利视频在线观看| 看片在线看免费视频| 国产伦人伦偷精品视频| 国产视频一区二区在线看| 日韩有码中文字幕| 最新中文字幕久久久久 | 三级国产精品欧美在线观看 | 一级a爱片免费观看的视频| 亚洲第一电影网av| 亚洲精品国产精品久久久不卡| 亚洲va日本ⅴa欧美va伊人久久| tocl精华| 色噜噜av男人的天堂激情| 可以在线观看毛片的网站| 亚洲国产精品成人综合色| 每晚都被弄得嗷嗷叫到高潮| 精品99又大又爽又粗少妇毛片 | 网址你懂的国产日韩在线| 国产麻豆成人av免费视频| 精华霜和精华液先用哪个| 人妻久久中文字幕网| www.熟女人妻精品国产| 在线观看一区二区三区| 亚洲国产中文字幕在线视频| 亚洲美女视频黄频| 亚洲avbb在线观看| 久久精品国产清高在天天线| 国产综合懂色| 美女黄网站色视频| 久久久精品大字幕| 一个人免费在线观看的高清视频| 国产精品亚洲一级av第二区| 日本撒尿小便嘘嘘汇集6| 丰满人妻一区二区三区视频av | 亚洲av成人一区二区三| 男人的好看免费观看在线视频| 日韩成人在线观看一区二区三区| 中文字幕久久专区| 中文字幕熟女人妻在线| 亚洲专区字幕在线| 亚洲aⅴ乱码一区二区在线播放| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲av第一区精品v没综合| 亚洲av五月六月丁香网| 亚洲aⅴ乱码一区二区在线播放| 久久久久久九九精品二区国产| 亚洲精品乱码久久久v下载方式 | 亚洲成a人片在线一区二区| 1024香蕉在线观看| 不卡av一区二区三区| 国产欧美日韩一区二区精品| 亚洲片人在线观看| 中文资源天堂在线| 好男人在线观看高清免费视频| 久久精品夜夜夜夜夜久久蜜豆| 日韩欧美在线乱码| 给我免费播放毛片高清在线观看| 97人妻精品一区二区三区麻豆| 一级毛片精品| 又大又爽又粗| av国产免费在线观看| 99国产精品99久久久久| 久久久久九九精品影院| 欧美在线黄色| 久久久久久久久中文| 免费观看的影片在线观看| 少妇丰满av| 韩国av一区二区三区四区| 真人一进一出gif抽搐免费| 草草在线视频免费看| 久久精品国产综合久久久| 免费观看的影片在线观看| 亚洲国产看品久久| 久久精品国产综合久久久| 国产欧美日韩一区二区精品| 久久99热这里只有精品18| 搞女人的毛片| 99久久国产精品久久久| 亚洲av免费在线观看| 一级黄色大片毛片| 欧美日韩福利视频一区二区| 91麻豆精品激情在线观看国产| 在线十欧美十亚洲十日本专区| 国产麻豆成人av免费视频| 亚洲av免费在线观看| 老司机在亚洲福利影院| 最新在线观看一区二区三区| 一二三四社区在线视频社区8| 一二三四在线观看免费中文在| 欧美中文日本在线观看视频|