• <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不卡免费在线播放| 日日啪夜夜撸| or卡值多少钱| 国产淫片久久久久久久久| 内地一区二区视频在线| 黄色日韩在线| 黄片wwwwww| 18禁在线播放成人免费| 国产精品一区www在线观看| 一卡2卡三卡四卡精品乱码亚洲| 人妻久久中文字幕网| 波多野结衣高清作品| 国产精品麻豆人妻色哟哟久久 | 国产精品.久久久| 亚洲av二区三区四区| 夜夜看夜夜爽夜夜摸| 欧美精品一区二区大全| 欧美+亚洲+日韩+国产| 91aial.com中文字幕在线观看| 成人三级黄色视频| 菩萨蛮人人尽说江南好唐韦庄 | 在线观看免费视频日本深夜| 少妇人妻一区二区三区视频| 激情 狠狠 欧美| 国产成人aa在线观看| 联通29元200g的流量卡| 日韩欧美 国产精品| 亚洲综合色惰| videossex国产| 午夜福利成人在线免费观看| 一级黄色大片毛片| 亚洲av男天堂| 18禁在线无遮挡免费观看视频| 干丝袜人妻中文字幕| 日韩在线高清观看一区二区三区| 国产三级中文精品| 亚洲内射少妇av| 精品一区二区三区视频在线| 赤兔流量卡办理| 国产乱人视频| 免费看美女性在线毛片视频| 精品日产1卡2卡| 精品少妇黑人巨大在线播放 | 亚洲精品国产成人久久av| 国产真实乱freesex| 免费观看精品视频网站| 亚洲精品自拍成人| 久久综合国产亚洲精品| 一卡2卡三卡四卡精品乱码亚洲| 精品久久久噜噜| 老女人水多毛片| 精品不卡国产一区二区三区| 亚洲综合色惰| 国产亚洲精品久久久com| 中文在线观看免费www的网站| 欧美不卡视频在线免费观看| 久久精品国产自在天天线| 欧美zozozo另类| 干丝袜人妻中文字幕| 久久久久久久久久成人| 久久亚洲精品不卡| 毛片女人毛片| 联通29元200g的流量卡| 国产午夜精品久久久久久一区二区三区| 国产一区二区激情短视频| 国产真实乱freesex| 国产精品野战在线观看| 日韩一本色道免费dvd| 欧美一区二区国产精品久久精品| 国产精品久久久久久精品电影小说 | 美女 人体艺术 gogo| 99久久精品一区二区三区| 99久久久亚洲精品蜜臀av| 国产亚洲91精品色在线| 亚洲欧美清纯卡通| 在线免费观看不下载黄p国产| 日韩大尺度精品在线看网址| 久久久午夜欧美精品| 国产精品精品国产色婷婷| 99热6这里只有精品| 国产在线男女| eeuss影院久久| 久久6这里有精品| 亚洲人成网站在线播| 人人妻人人澡人人爽人人夜夜 | 婷婷精品国产亚洲av| 日韩制服骚丝袜av| 国产真实乱freesex| 少妇的逼水好多| 自拍偷自拍亚洲精品老妇| 偷拍熟女少妇极品色| 桃色一区二区三区在线观看| www.av在线官网国产| 一区福利在线观看| 69av精品久久久久久| 免费黄网站久久成人精品| 人人妻人人澡人人爽人人夜夜 | 亚洲国产高清在线一区二区三| 色吧在线观看| 国产麻豆成人av免费视频| 久久久久久久久中文| 亚洲真实伦在线观看| 国产精品女同一区二区软件| 亚洲不卡免费看| 中文亚洲av片在线观看爽| 黑人高潮一二区| 麻豆国产av国片精品| 国产精品野战在线观看| 久久中文看片网| 精品久久久久久久末码| av黄色大香蕉| 日韩人妻高清精品专区| 久久综合国产亚洲精品| 精品久久久久久久久久久久久| 人人妻人人看人人澡| 我的老师免费观看完整版| 久久久久久久午夜电影| 99热这里只有是精品在线观看| 午夜免费男女啪啪视频观看| 亚洲精品成人久久久久久| 免费人成视频x8x8入口观看| 国产乱人偷精品视频| 高清午夜精品一区二区三区 | 精品国内亚洲2022精品成人| 色视频www国产| 非洲黑人性xxxx精品又粗又长| 女的被弄到高潮叫床怎么办| 日本黄色视频三级网站网址| 成年女人永久免费观看视频| 欧美成人免费av一区二区三区| 国产在线男女| 草草在线视频免费看| 男人的好看免费观看在线视频| 又爽又黄a免费视频| 成人毛片60女人毛片免费| 97在线视频观看| 99热全是精品| 女同久久另类99精品国产91| 久久久久久国产a免费观看| 免费看av在线观看网站| 国产精品.久久久| 91av网一区二区| 日本免费a在线| 成熟少妇高潮喷水视频| 精品熟女少妇av免费看| 青青草视频在线视频观看| 蜜桃久久精品国产亚洲av| 亚洲成人av在线免费| 国产高清视频在线观看网站| 一级黄片播放器| 亚洲国产精品合色在线| 久久人人爽人人爽人人片va| 特大巨黑吊av在线直播| 久久久久久九九精品二区国产| 色综合亚洲欧美另类图片| 99久国产av精品国产电影| kizo精华| 99久久中文字幕三级久久日本| 国产精品嫩草影院av在线观看| 精品久久久噜噜| 自拍偷自拍亚洲精品老妇| 国产精品久久久久久av不卡| 国产在线男女| 国产精品一及| 久久久久性生活片| 在线天堂最新版资源| 亚洲国产欧美人成| 国产在视频线在精品| 日韩人妻高清精品专区| 日韩欧美精品免费久久| 亚洲,欧美,日韩| 精品一区二区免费观看| 女人被狂操c到高潮| 三级男女做爰猛烈吃奶摸视频| 久久中文看片网| 日韩一区二区三区影片| 人人妻人人澡人人爽人人夜夜 | 亚洲第一区二区三区不卡| 亚洲精品国产成人久久av| 99热这里只有精品一区| 此物有八面人人有两片| 男插女下体视频免费在线播放| 大又大粗又爽又黄少妇毛片口| 男人舔奶头视频| 女人十人毛片免费观看3o分钟| 99久久人妻综合| 悠悠久久av| 少妇裸体淫交视频免费看高清| 亚洲精品456在线播放app| 国产亚洲精品久久久久久毛片| 久久人人爽人人爽人人片va| 日本三级黄在线观看| 99久久久亚洲精品蜜臀av| 搡老妇女老女人老熟妇| 美女被艹到高潮喷水动态| 哪个播放器可以免费观看大片| 国产精品久久久久久av不卡| 欧美日韩国产亚洲二区| av免费在线看不卡| 特大巨黑吊av在线直播| 2022亚洲国产成人精品| 老司机影院成人| 亚洲在久久综合| 丰满人妻一区二区三区视频av| 日韩 亚洲 欧美在线| 婷婷亚洲欧美| 一级二级三级毛片免费看| 亚洲国产精品国产精品| 97人妻精品一区二区三区麻豆| av卡一久久| 女的被弄到高潮叫床怎么办| 欧美又色又爽又黄视频| videossex国产| 日韩在线高清观看一区二区三区| 亚洲综合色惰| 国产精品一二三区在线看| 人人妻人人澡人人爽人人夜夜 | 日韩精品青青久久久久久| 久久韩国三级中文字幕| 欧美+日韩+精品| 爱豆传媒免费全集在线观看| 精品99又大又爽又粗少妇毛片| 不卡一级毛片| 97热精品久久久久久| 成人午夜精彩视频在线观看| 免费不卡的大黄色大毛片视频在线观看 | 深夜精品福利| 亚洲美女搞黄在线观看| 日韩精品青青久久久久久| videossex国产| 日韩成人伦理影院| 在线免费观看的www视频| 女人被狂操c到高潮| 精品一区二区三区人妻视频| 亚洲高清免费不卡视频| 国产成人freesex在线| 久久精品影院6| 成年av动漫网址| 国产精品久久久久久亚洲av鲁大| 国产精品久久久久久av不卡| 一个人看视频在线观看www免费| 美女黄网站色视频| 色尼玛亚洲综合影院| 只有这里有精品99| 亚洲欧美中文字幕日韩二区| 国产熟女欧美一区二区| 亚洲欧美日韩卡通动漫| 国产免费男女视频| 91久久精品国产一区二区三区| 蜜桃久久精品国产亚洲av| 久久欧美精品欧美久久欧美| 在线免费十八禁| 人人妻人人澡人人爽人人夜夜 | av女优亚洲男人天堂| 伦理电影大哥的女人| 色综合站精品国产| 亚洲在线观看片| av国产免费在线观看| 亚洲欧美中文字幕日韩二区| av天堂中文字幕网| 免费观看人在逋| 免费观看精品视频网站| 日韩三级伦理在线观看| 尤物成人国产欧美一区二区三区| 国产精品无大码| 亚洲欧美清纯卡通| 夫妻性生交免费视频一级片| 亚洲七黄色美女视频| 赤兔流量卡办理| 免费电影在线观看免费观看| 国产av一区在线观看免费| 2021天堂中文幕一二区在线观| 天天躁日日操中文字幕| 天堂网av新在线| 在线观看免费视频日本深夜| 黄片无遮挡物在线观看| 搡老妇女老女人老熟妇| 精品99又大又爽又粗少妇毛片| 久久久久免费精品人妻一区二区| av天堂在线播放| 日日摸夜夜添夜夜爱| 国产精品永久免费网站| 免费观看在线日韩| 国产精品女同一区二区软件| 国产高清三级在线| 哪个播放器可以免费观看大片| 午夜福利视频1000在线观看| 不卡视频在线观看欧美| 国产亚洲91精品色在线| 内地一区二区视频在线| 久久欧美精品欧美久久欧美| 天天躁日日操中文字幕| 在线播放无遮挡| 亚州av有码| 99久国产av精品| 欧美一区二区精品小视频在线| 久久6这里有精品| 亚洲人成网站在线观看播放| 亚洲综合色惰| 欧美3d第一页| 日本-黄色视频高清免费观看| 欧美性猛交黑人性爽| 国产亚洲av片在线观看秒播厂 | 如何舔出高潮| ponron亚洲| 18禁在线无遮挡免费观看视频| 黄色欧美视频在线观看| av天堂在线播放| 夜夜爽天天搞| 国产精品麻豆人妻色哟哟久久 | 国产精品乱码一区二三区的特点| 亚洲av男天堂| 99久国产av精品| 男人舔奶头视频| 毛片女人毛片| 网址你懂的国产日韩在线| 色综合站精品国产| 国产又黄又爽又无遮挡在线| 少妇人妻一区二区三区视频| 可以在线观看的亚洲视频| 一级黄色大片毛片| 天美传媒精品一区二区| 天堂av国产一区二区熟女人妻| 男女做爰动态图高潮gif福利片| 精品人妻熟女av久视频| 在线天堂最新版资源| 女人十人毛片免费观看3o分钟| 精品久久久噜噜| 久久精品国产亚洲av香蕉五月| 精品无人区乱码1区二区| 亚洲国产高清在线一区二区三| 国产精品女同一区二区软件| 成人美女网站在线观看视频| 久久99蜜桃精品久久| 禁无遮挡网站| 免费看光身美女| 国产成人精品婷婷| 久久精品国产99精品国产亚洲性色| 国产精品1区2区在线观看.| 亚洲国产欧美在线一区| 亚洲欧美精品综合久久99| 日本黄大片高清| 日本免费a在线| 欧美最黄视频在线播放免费| 亚洲精品日韩av片在线观看| 男人和女人高潮做爰伦理| 天天一区二区日本电影三级| 国产成人精品婷婷| 精品99又大又爽又粗少妇毛片| av在线蜜桃| 日韩av在线大香蕉| 99精品在免费线老司机午夜| 久久精品夜夜夜夜夜久久蜜豆| 男人和女人高潮做爰伦理| 成人午夜精彩视频在线观看| 久久精品影院6| 亚洲在线观看片| 亚州av有码| 日本一本二区三区精品| 久久精品影院6| 啦啦啦观看免费观看视频高清| 久久精品国产清高在天天线| av在线老鸭窝| 国产成人精品婷婷| 麻豆一二三区av精品| 色哟哟·www| 99热精品在线国产| 一级毛片电影观看 | 日韩国内少妇激情av| 超碰av人人做人人爽久久| 亚洲,欧美,日韩| 日韩成人av中文字幕在线观看| 欧美高清成人免费视频www| 日韩强制内射视频| av在线亚洲专区| 偷拍熟女少妇极品色| 麻豆成人av视频| 校园人妻丝袜中文字幕| 老熟妇乱子伦视频在线观看| 国产精品久久视频播放| 免费观看在线日韩| 久99久视频精品免费| 午夜福利高清视频| 精品无人区乱码1区二区| 日本五十路高清| 久久精品夜色国产| 2022亚洲国产成人精品| 久久草成人影院| 美女大奶头视频| 亚洲欧美精品综合久久99| 欧美bdsm另类| 在线播放国产精品三级| 亚州av有码| 大又大粗又爽又黄少妇毛片口| 亚洲18禁久久av| 久久久精品大字幕| 精品久久久久久久久av| 欧美高清成人免费视频www| 美女xxoo啪啪120秒动态图| 国产免费男女视频| 夜夜爽天天搞| 亚洲电影在线观看av| 在现免费观看毛片| 男人舔奶头视频| 亚洲成a人片在线一区二区| 亚洲欧美日韩无卡精品| 欧美极品一区二区三区四区| 我的女老师完整版在线观看| 99精品在免费线老司机午夜| 青春草亚洲视频在线观看| 亚洲av中文字字幕乱码综合| 国产精品免费一区二区三区在线| 亚洲最大成人手机在线| 久久久久久大精品| 热99re8久久精品国产| 欧美激情久久久久久爽电影| 国产av不卡久久| 级片在线观看| 变态另类丝袜制服| 亚洲欧美成人精品一区二区| 岛国在线免费视频观看| 超碰av人人做人人爽久久| 特级一级黄色大片| 18禁在线无遮挡免费观看视频| 九草在线视频观看| 一边摸一边抽搐一进一小说| 我的老师免费观看完整版| 国产麻豆成人av免费视频| av在线观看视频网站免费| 欧美成人精品欧美一级黄| 国内揄拍国产精品人妻在线| 国产亚洲av片在线观看秒播厂 | 不卡视频在线观看欧美| 成人无遮挡网站| 亚洲欧美精品自产自拍| 中文字幕精品亚洲无线码一区| 神马国产精品三级电影在线观看| 国产亚洲欧美98| 美女大奶头视频| 久久精品国产亚洲av天美| 国产精品国产三级国产av玫瑰| 国产成人a区在线观看| 又粗又硬又长又爽又黄的视频 | 美女黄网站色视频| 精品久久久久久久人妻蜜臀av| 久久久久久久久久黄片| 级片在线观看| 日本黄大片高清| 国产女主播在线喷水免费视频网站 | 国产av一区在线观看免费| 国产精品国产三级国产av玫瑰| 国产在视频线在精品| 长腿黑丝高跟| 国产成人一区二区在线| 欧美在线一区亚洲| 欧美三级亚洲精品| 秋霞在线观看毛片| 日韩精品有码人妻一区| 亚洲人成网站高清观看| 久久久精品欧美日韩精品| 午夜免费男女啪啪视频观看| 69人妻影院| 欧美日韩综合久久久久久| av女优亚洲男人天堂| 最好的美女福利视频网| 免费黄网站久久成人精品| 国产成年人精品一区二区| 免费一级毛片在线播放高清视频| 欧美一区二区国产精品久久精品| 欧美日韩国产亚洲二区| 欧美成人一区二区免费高清观看| 日本黄大片高清| 在线免费观看不下载黄p国产| 成人毛片a级毛片在线播放| 麻豆精品久久久久久蜜桃| 国内揄拍国产精品人妻在线| 99热精品在线国产| 亚洲av第一区精品v没综合| 国产高清不卡午夜福利| 一进一出抽搐动态| 3wmmmm亚洲av在线观看| 黄色欧美视频在线观看| 在线观看美女被高潮喷水网站| 少妇高潮的动态图| 麻豆国产av国片精品| 黄片无遮挡物在线观看| 亚洲欧美精品综合久久99| 欧美bdsm另类| 别揉我奶头 嗯啊视频| 久久久久久久久久黄片| 亚洲欧美成人精品一区二区| 久久午夜福利片| 一区二区三区高清视频在线| 中出人妻视频一区二区| 五月玫瑰六月丁香| 内射极品少妇av片p| 最近手机中文字幕大全| 久久久午夜欧美精品| 高清毛片免费看| 床上黄色一级片| 国产探花极品一区二区| 日本欧美国产在线视频| 亚洲精品日韩在线中文字幕 | 大香蕉久久网| 亚洲av免费高清在线观看| 人人妻人人澡欧美一区二区| 亚洲国产精品国产精品| 亚洲精品国产成人久久av| 免费无遮挡裸体视频| 亚洲av免费高清在线观看| 亚洲乱码一区二区免费版| 久久精品国产鲁丝片午夜精品| 97超视频在线观看视频| 麻豆一二三区av精品| 亚洲aⅴ乱码一区二区在线播放| 全区人妻精品视频| 日本黄色片子视频| 国产精品国产三级国产av玫瑰| 久久这里有精品视频免费| 欧美激情在线99| 日韩欧美 国产精品| 老司机福利观看| 哪里可以看免费的av片| 人妻制服诱惑在线中文字幕| 麻豆成人午夜福利视频| 在线观看一区二区三区| 联通29元200g的流量卡| 99国产精品一区二区蜜桃av| 悠悠久久av| 成人亚洲精品av一区二区| 亚洲内射少妇av| 少妇高潮的动态图| 一进一出抽搐动态| 久久这里只有精品中国| 国产成人午夜福利电影在线观看| 亚洲aⅴ乱码一区二区在线播放| 日本-黄色视频高清免费观看| 一级黄片播放器| 高清日韩中文字幕在线| 麻豆成人av视频| 久久九九热精品免费| 91av网一区二区| 中文字幕久久专区| 69人妻影院| 一进一出抽搐gif免费好疼| 99热这里只有精品一区| 成人特级av手机在线观看| 久久久久国产网址| 亚洲成人久久爱视频| 国产精品av视频在线免费观看| 国内精品久久久久精免费| 深爱激情五月婷婷| 亚洲中文字幕一区二区三区有码在线看| 狂野欧美白嫩少妇大欣赏| 搡女人真爽免费视频火全软件| 欧美最黄视频在线播放免费| 99国产精品一区二区蜜桃av| 大型黄色视频在线免费观看| 精品久久久久久久久久久久久| a级毛片免费高清观看在线播放| 大香蕉久久网| 美女 人体艺术 gogo| 99热网站在线观看| 欧美日韩国产亚洲二区| 老司机影院成人| 男女下面进入的视频免费午夜| 欧美最黄视频在线播放免费| 国产成人a∨麻豆精品| 亚洲精品久久久久久婷婷小说 | 亚洲av中文字字幕乱码综合| 一区二区三区免费毛片| 91aial.com中文字幕在线观看| 国产综合懂色| 观看美女的网站| 三级国产精品欧美在线观看| 国产精品久久电影中文字幕| 国产视频首页在线观看| 欧美三级亚洲精品| 欧美极品一区二区三区四区| 国产精品爽爽va在线观看网站| 国产黄片美女视频| 桃色一区二区三区在线观看| 国产高清视频在线观看网站| 人人妻人人澡人人爽人人夜夜| 国产黄色视频一区二区在线观看| 亚洲欧洲国产日韩| 免费看不卡的av| 我要看黄色一级片免费的| 日韩中字成人| 午夜免费男女啪啪视频观看| 我要看黄色一级片免费的| 欧美bdsm另类| 亚洲欧美成人精品一区二区| 女性被躁到高潮视频| 日本wwww免费看| 久久久久久久久久成人| 精品人妻熟女av久视频| 十分钟在线观看高清视频www| 人人妻人人爽人人添夜夜欢视频| 一本久久精品| 日韩伦理黄色片| 国产精品不卡视频一区二区| 老司机亚洲免费影院| 最新中文字幕久久久久| 中国美白少妇内射xxxbb| 精品一区二区三区视频在线| 国产成人91sexporn| 日本黄色日本黄色录像| 最近的中文字幕免费完整| 亚州av有码| 国产免费现黄频在线看| 国产日韩欧美视频二区| 久久久国产一区二区| 最近中文字幕2019免费版|