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

    多種因素影響下的船舶航線誘導(dǎo)模型構(gòu)建與應(yīng)用

    2020-10-29 08:44:08閆建鑫邢江豪
    科學(xué)技術(shù)與工程 2020年25期
    關(guān)鍵詞:航行航線路段

    張 赫, 閆建鑫, 邢江豪

    (大連海事大學(xué)交通運(yùn)輸工程學(xué)院, 大連 116026)

    隨著世界經(jīng)濟(jì)和技術(shù)的快速發(fā)展,各種類型船舶的數(shù)量隨之持續(xù)增加,船舶之間交通沖突和擁擠的現(xiàn)象愈發(fā)嚴(yán)重。例如海上船舶在行駛到工程船舶的作業(yè)區(qū)域時(shí),大量的工程船舶根據(jù)作業(yè)地點(diǎn)的不同需要多次穿越航道,這與往來(lái)正常航行的其他船舶形成了一個(gè)復(fù)雜的航線網(wǎng)絡(luò),船舶的航行受到極大的影響。在這種情況下,需要對(duì)船舶的行駛路線進(jìn)行合理規(guī)劃,而誘導(dǎo)技術(shù)是實(shí)現(xiàn)路徑優(yōu)化的有效方法之一。

    路徑誘導(dǎo)技術(shù)自20世紀(jì)產(chǎn)生到現(xiàn)在已經(jīng)幾十年的歷史,相關(guān)的資料也比較完善,目前更多的應(yīng)用于車輛行駛的誘導(dǎo)。對(duì)于汽車的路徑誘導(dǎo)中外相關(guān)領(lǐng)域的專家采用不同方式和方法進(jìn)行了大量研究[1-4],而針對(duì)于船舶的誘導(dǎo)研究較少,僅有喬春福[5]研究了如何根據(jù)船位序列判斷船舶上行和下行的算法,從而為過(guò)往船舶的上下行助航提供誘導(dǎo)幫助。從文獻(xiàn)中可以看出,目前對(duì)船舶的航線誘導(dǎo)還在發(fā)展之中,僅有的研究也只是考慮了行駛水域的特點(diǎn)結(jié)合地圖進(jìn)行上下行的助航,沒(méi)有針對(duì)于某一船舶在復(fù)雜的水域內(nèi)行駛的航線誘導(dǎo)。而陸上汽車的路徑誘導(dǎo)問(wèn)題經(jīng)過(guò)多年的發(fā)展已經(jīng)趨于成熟,相關(guān)理論也相對(duì)完善,這對(duì)實(shí)現(xiàn)船舶的航線誘導(dǎo)具有指導(dǎo)意義。現(xiàn)在考慮與陸上汽車誘導(dǎo)所不同的氣象水文條件影響的基礎(chǔ)上結(jié)合船舶航行的特點(diǎn)將誘導(dǎo)技術(shù)運(yùn)用到船舶的航線規(guī)劃之中,以期為船舶航行提供科學(xué)依據(jù)。

    1 海況對(duì)船舶航行的影響

    船舶在海上航行時(shí),空氣和水對(duì)船舶行駛具有很大的影響,為了便于研究船舶在不同的風(fēng)浪條件下的行駛速度,將船舶阻力分為基本阻力、空氣阻力、洶濤阻力三部分。

    1.1 基本阻力

    船舶基本阻力是船舶在靜水中航行受到的阻力,由摩擦阻力和壓阻力(剩余阻力)組成。結(jié)合傅汝德相似定律和雷諾相似定律可得船舶靜水行駛基本阻力R0的計(jì)算表達(dá)式[6]:

    (1)

    CR=Cf+ΔCf+Cr

    (2)

    式中:CR為基本阻力系數(shù),ρ0為流體密度,V0為船舶在靜水中的行駛速度;S為船舶濕面積;Cf為船舶摩擦阻力系數(shù);ΔCf為粗糙度補(bǔ)貼系數(shù);Cr為船舶剩余阻力系數(shù)。

    摩擦阻力系數(shù)Cf用“1957年國(guó)際船模實(shí)驗(yàn)池實(shí)船-船模換算公式”[式(3)]進(jìn)行求解, 粗糙度補(bǔ)貼系數(shù)ΔCf如表1所示;剩余阻力系數(shù)Cr利用Lap-kelller圖譜進(jìn)行查表計(jì)算獲得[7]。

    (3)

    (4)

    式中:Re為雷諾數(shù);V為船速;Lwl為水線長(zhǎng),近似于船長(zhǎng)L;ν為水運(yùn)動(dòng)黏性系數(shù)。

    表1 粗糙度補(bǔ)貼系數(shù)Table 1 Roughness subsidy coefficients

    1.2 空氣阻力

    船舶在海上航行時(shí),其水面以上的部分在風(fēng)的作用下會(huì)產(chǎn)生空氣阻力影響行駛速度,產(chǎn)生船舶艏艉方向的風(fēng)壓力Fw可按式(5)計(jì)算[8]:

    (5)

    (6)

    (7)

    (8)

    lgY=αw+βwlgX

    (9)

    式中:ρa(bǔ)為空氣密度(1.226 kg/m3);Cw為艏艉方向上風(fēng)力系數(shù);Va為相對(duì)風(fēng)速;Aa為水線以上船體正面受風(fēng)面積;Ba為水線以上船體側(cè)面受風(fēng)面積;Vzf為真風(fēng)風(fēng)速;αa為真風(fēng)方向角,左舷為正;θ為風(fēng)舷角;Ca為風(fēng)力系數(shù);α為風(fēng)力角,是風(fēng)壓力與船舶艏艉線之間的夾角;Y為船舶受風(fēng)面積(包括Aa、Ba);X表示船舶載重量或總質(zhì)量,t;αw、βw為回歸系數(shù),其數(shù)值如表2所示[9]。

    表2 不同船型船舶的受風(fēng)面積回歸系數(shù)Table 2 Regression coefficient of windage area of different types of ships

    風(fēng)力角α和風(fēng)力系數(shù)Ca可根據(jù)巖井聰公式[式(10)和式(11)]進(jìn)行計(jì)算[6-8]:

    (10)

    Ca=1.325-0.05cos(2α)-0.35cos(4α)-0.175cos(6α)

    (11)

    1.3 洶濤阻力

    洶濤阻力分為規(guī)則波干擾和不規(guī)則波干擾兩種,通常只考慮規(guī)則波的影響。對(duì)于規(guī)則波干擾力的計(jì)算比較復(fù)雜,目前主要依靠船模實(shí)驗(yàn)進(jìn)行,Daidola在船模實(shí)驗(yàn)的基礎(chǔ)上提出了洶濤阻力的計(jì)算模型為[6-7]

    (12)

    (13)

    1.4 船舶自由行駛速度

    假定船舶在靜水中行駛的推進(jìn)功率和風(fēng)浪中行駛的推進(jìn)功率保持不變,則船舶在風(fēng)浪中行駛減少的速度用來(lái)抵消風(fēng)浪阻力的影響[6],從而得到船舶在風(fēng)浪下自由航行速度Vf。

    (14)

    由于船舶在不同路段行駛的位置和方向并不一致這就導(dǎo)致相對(duì)風(fēng)速的大小、方向發(fā)生變化,因此船舶在不同路段行駛時(shí)自由航行速度Vf需要根據(jù)行駛方向和位置分別求解。

    2 路段行程時(shí)間計(jì)算

    在對(duì)船舶進(jìn)行航線誘導(dǎo)時(shí),優(yōu)化目標(biāo)通常有距離、時(shí)間和費(fèi)用三種,通過(guò)綜合考慮,現(xiàn)以航線用時(shí)最短作為優(yōu)化目標(biāo)。

    各路段行程時(shí)間分為路段行駛時(shí)間和航線網(wǎng)絡(luò)交叉點(diǎn)延誤時(shí)間兩個(gè)部分。在沒(méi)有航線交叉的路段內(nèi),路段行程時(shí)間就等于路段行駛時(shí)間。船舶在路段的行駛時(shí)間與路段船舶流量密切相關(guān),根據(jù)風(fēng)浪對(duì)船舶航行的影響對(duì)美國(guó)聯(lián)邦公路局函數(shù)即BPR函數(shù)進(jìn)行改進(jìn),建立路段行駛時(shí)間的表達(dá)式:

    (15)

    式(15)中:tx為兩節(jié)點(diǎn)之間路段行駛時(shí)間;ψ為風(fēng)浪對(duì)航速的影響系數(shù),ψ=Vf/V0;t0為兩節(jié)點(diǎn)之間交通量為零時(shí)無(wú)風(fēng)浪影響下的路段行駛時(shí)間;q為路段船舶交通量;c為航路實(shí)際通行能力;γ1、γ2為相關(guān)系數(shù)。

    航線交叉點(diǎn)平均延誤時(shí)間td根據(jù)排隊(duì)論原理進(jìn)行求解。將船舶看作“顧客”,交叉點(diǎn)視為“服務(wù)窗”,船舶在交叉點(diǎn)的排隊(duì)穿越過(guò)程可視為一個(gè)排隊(duì)問(wèn)題。船舶到達(dá)交叉點(diǎn)時(shí),若“服務(wù)窗”被占用,則進(jìn)行排隊(duì)等待,最大排隊(duì)長(zhǎng)度n為航道長(zhǎng)度和單位船舶所需的航道長(zhǎng)度Ld的比值。對(duì)船舶到達(dá)規(guī)律進(jìn)行統(tǒng)計(jì)分析,船舶的到達(dá)規(guī)律可近似服從參數(shù)為λ的泊松分布,船舶通過(guò)交叉點(diǎn)的服務(wù)時(shí)間可看作服從參數(shù)為μ的負(fù)指數(shù)分布(到達(dá)率λ=q,μ可統(tǒng)計(jì)獲得)[10]。因此交叉點(diǎn)的船舶排隊(duì)模型為M/M/1/n/∞模型。

    (16)

    綜上,ij路段總行程時(shí)間tij為

    tij=tx+td

    (17)

    3 模型構(gòu)建

    3.1 影響因素分析

    3.1.1 施工水域影響

    在海上航行的船舶可能會(huì)遇到海上施工區(qū)域,影響船舶航行。由于施工區(qū)域內(nèi)有大量工程船舶進(jìn)行運(yùn)輸作業(yè),使得施工區(qū)內(nèi)航線密集,交叉增多。

    3.1.2 海風(fēng)動(dòng)態(tài)變化

    船舶位置不同,船舶受到風(fēng)力大小也不一致,為此建立船舶風(fēng)速與位置的對(duì)應(yīng)關(guān)系公式Vzf=φ1Vfe,其中φ1為與位置相關(guān)的風(fēng)力作用系數(shù)(0≤φ1≤1),Vfe為基本風(fēng)速。

    3.1.3 波浪動(dòng)態(tài)變化

    波浪振幅的變化與海風(fēng)息息相關(guān),并且影響洶濤阻力的大小。為此建立波浪振幅與位置的對(duì)應(yīng)關(guān)系公式Wa=φ2W0,其中φ2為與位置相關(guān)的波浪作用系數(shù)(0≤φ2≤1),W0為基本波浪振幅。

    3.1.4 沉點(diǎn)的影響

    對(duì)船舶航行產(chǎn)生阻礙的物體所在的位置稱為沉點(diǎn)。在施工區(qū)域內(nèi),進(jìn)行挖泥、疏浚作業(yè)的船舶通常移動(dòng)速度較慢可視為沉點(diǎn),若沉點(diǎn)存在于路段中相應(yīng)路阻為無(wú)窮大。淺灘、島嶼等存在的位置也是沉點(diǎn),為禁止行駛區(qū)域需繞開(kāi)行駛。

    3.2 模型建立

    在對(duì)船舶航線誘導(dǎo)時(shí),將航線網(wǎng)絡(luò)抽象成圖論網(wǎng)絡(luò)G=(P,E,T),其中P是節(jié)點(diǎn)集合表示航道節(jié)點(diǎn);E是弧的集合表示路段;T={tij(Ti)|(i,j)∈E}是各弧的權(quán)值,表示在此路段的行程時(shí)間tij(Ti)為船舶在Ti時(shí)刻從節(jié)點(diǎn)i航行到節(jié)點(diǎn)j所需花費(fèi)的時(shí)間。設(shè)施工區(qū)k(k=1,2,…,l)的節(jié)點(diǎn)集合為Pk存在的時(shí)間為[Tks,Tke],從起點(diǎn)到達(dá)施工區(qū)k的時(shí)間為tk,非施工區(qū)節(jié)點(diǎn)結(jié)合為P0。則建立船舶動(dòng)態(tài)時(shí)變網(wǎng)絡(luò)下的船舶航線誘導(dǎo)優(yōu)化模型為

    (18)

    s.t.

    (19)

    (20)

    (21)

    (22)

    (23)

    k=1,2,…,l

    (24)

    (25)

    Vzf=φ1Vfe

    (26)

    Wa=φ2W0

    (27)

    式(18)為優(yōu)化目標(biāo);式(19)表示路段ij在時(shí)刻被路徑使用則等于1,否則等于0;式(20)表示初始點(diǎn)只能是出發(fā)點(diǎn)不能是到達(dá)點(diǎn),終點(diǎn)只能是到達(dá)點(diǎn)不能是出發(fā)點(diǎn),中間節(jié)點(diǎn)的到達(dá)和出發(fā)次數(shù)相等;式(21)表示所有路段最多只能被選擇一次;式(22)為施工區(qū)內(nèi)沉點(diǎn)對(duì)行程時(shí)間的影響;式(23)、式(24)表示船舶行駛到施工區(qū)域位置時(shí)若在施工區(qū)域存在的時(shí)間范圍內(nèi)則將區(qū)域內(nèi)節(jié)點(diǎn)加入到航線網(wǎng)絡(luò)之中,否則保留原節(jié)點(diǎn);式(25)表示船舶各路段自由行駛速度的計(jì)算;式(26)為風(fēng)速與位置的對(duì)應(yīng)關(guān)系式;式(27)為波浪振幅與位置的對(duì)應(yīng)關(guān)系式;行程時(shí)間計(jì)算過(guò)程中的其他參數(shù)的計(jì)算公式約束第1節(jié)與第2節(jié)已經(jīng)說(shuō)明,此處不再列出。

    對(duì)于各路段船舶流量進(jìn)行實(shí)時(shí)的更新是非常困難的,通常將時(shí)間根據(jù)實(shí)際劃分為多個(gè)時(shí)段,并假定在同一個(gè)時(shí)段航路網(wǎng)各路段的船舶流量不發(fā)生改變。

    4 蟻群算法求解

    (28)

    螞蟻釋放在各節(jié)點(diǎn)之間的信息素隨著時(shí)間的流逝也會(huì)逐漸消失,設(shè)參數(shù)ε3(0<ε3<1)代表信息素的揮發(fā)系數(shù)。當(dāng)m只螞蟻都完成了一次循環(huán)后,需要對(duì)各節(jié)點(diǎn)間t時(shí)刻的τij(t)進(jìn)行迭代更新,公式為

    (29)

    (30)

    式(30)中:Q為螞蟻在一次循環(huán)中所散發(fā)的信息素總量;Lk為第k只螞蟻在本次循環(huán)中所花費(fèi)的總時(shí)間。

    針對(duì)于本文航線誘導(dǎo)問(wèn)題,以船舶出發(fā)點(diǎn)為螞蟻路徑起點(diǎn)o,以船舶目標(biāo)點(diǎn)為螞蟻路徑尋找的終點(diǎn),以螞蟻選擇的節(jié)點(diǎn)集合表示船舶路徑的可行解。假設(shè)初始時(shí)刻T=T0,則用蟻群算法求解的基本過(guò)程如下:

    (1)對(duì)螞蟻總數(shù)目m、信息素影響因子ε1、啟發(fā)函數(shù)影響因子ε2、信息揮發(fā)因子ε3、迭代次數(shù)t、最大迭代次數(shù)tmax等參數(shù)進(jìn)行設(shè)置,令初始時(shí)刻m只螞蟻所在節(jié)點(diǎn)x=o。

    (2)根據(jù)生成的施工區(qū)域和各節(jié)點(diǎn)的坐標(biāo),計(jì)算各節(jié)點(diǎn)間的行駛方向和風(fēng)浪的相對(duì)大小和方向。

    (3)根據(jù)氣象水文條件和船舶本身狀況計(jì)算船舶在各路段的自由行駛速度Vf,并據(jù)此求出對(duì)應(yīng)的風(fēng)浪影響系數(shù)。

    (4)螞蟻根據(jù)與x相連的可行解節(jié)點(diǎn)在T時(shí)刻的船舶流量和風(fēng)浪影響系數(shù)等計(jì)算從節(jié)點(diǎn)x到可行節(jié)點(diǎn)各自的行程時(shí)間,然后按照式(28)依概率選擇下一節(jié)點(diǎn)。

    (5)令T=T+Txr(r為選擇的下一節(jié)點(diǎn),Txr為路段xr行程時(shí)間)、x=r,返回(4)繼續(xù)選擇下一節(jié)點(diǎn)直到到達(dá)目標(biāo)點(diǎn)。

    (6)螞蟻k(k=1,2,…,m)到達(dá)目標(biāo)點(diǎn)則按照實(shí)際用時(shí)記錄此路徑,依次完成所有m只螞蟻循環(huán),記錄此次迭代中的最優(yōu)路徑及用時(shí),并按照式(29)對(duì)信息素進(jìn)行更新。

    (7)令t=t+1,x=o,T=T0若t≤tmax轉(zhuǎn)向(4),進(jìn)行下一次迭代,否則轉(zhuǎn)向(8);

    (8)輸出最優(yōu)解。

    5 算例

    在此給出如圖1所示的從煙臺(tái)到天津的簡(jiǎn)易模擬航線路網(wǎng)(箭頭表示只能單向通行),陰影區(qū)域?yàn)槌咙c(diǎn)存在的區(qū)域,用MATLAB編程求解。利用隨機(jī)算法在航線網(wǎng)絡(luò)中隨機(jī)生成3個(gè)長(zhǎng)度為6~8 km的互不影響矩形施工區(qū)域,施工區(qū)域存在的時(shí)間為[T0+Xk,T0+Yk](0≤Xk≤Yk≤24 h)。在各施工區(qū)域內(nèi)隨機(jī)生成一個(gè)施工船沉點(diǎn),若沉點(diǎn)在航線路段當(dāng)中,則該路段行駛時(shí)間為無(wú)窮大,以15 min為一個(gè)時(shí)段。

    圖1 模擬航線網(wǎng)Fig.1 Analog ship route network

    設(shè)置起點(diǎn)的標(biāo)號(hào)為0,終點(diǎn)的標(biāo)號(hào)為1 000,1號(hào)航線的各點(diǎn)標(biāo)號(hào)為101、102、103、104、105,同理為其他各航線進(jìn)行標(biāo)號(hào)。圖1中在1、3、5航線內(nèi)生成的施工區(qū)航線網(wǎng)絡(luò)及相應(yīng)標(biāo)號(hào)分別如圖2(a)~圖2(c)所示。

    圖2 海上施工區(qū)域航線網(wǎng)絡(luò)Fig.2 Route network in offshore construction area

    現(xiàn)選擇某型號(hào)散貨船進(jìn)行計(jì)算,該船舶船長(zhǎng)180 m,船寬32.2 m,載重量為48 103 t,滿載排水量54 385 t,吃水11.6 m,靜水營(yíng)運(yùn)航速為6.5 m/s,真風(fēng)方向?yàn)槲髌?0°,波長(zhǎng)取0.8倍船長(zhǎng)。

    在用蟻群算法求解時(shí),取m=30,ε1=1,ε2=1,ε3=0.1,tmax=200, 求得時(shí)變網(wǎng)絡(luò)中各代最短用時(shí)和平均用時(shí)隨迭代次數(shù)的變化曲線如圖3所示。

    圖3 時(shí)變網(wǎng)絡(luò)各代最短用時(shí)和平均用時(shí)對(duì)比Fig.3 Comparison of shortest and average time of each generation in time-varying networks

    最終通過(guò)MTALAB編程求得動(dòng)態(tài)網(wǎng)絡(luò)的最佳路徑為0→501→502→503→71→81→91→92→93→94→95→504→505→1 000,用時(shí)86 937 s。在相同情況下,求得1、2、3、4航線的最短用時(shí)分別為90 958、91 704、91 818、92 504 s。最佳路徑比4航線路徑少用時(shí)5 567 s,減少用時(shí)6.4%。

    6 結(jié)論

    通過(guò)建立船舶的航線誘導(dǎo)模型利用MATLAB軟件進(jìn)行算例仿真驗(yàn)證,說(shuō)明了本文提出的船舶航線誘導(dǎo)方法能為船舶在復(fù)雜航線網(wǎng)絡(luò)下行駛提供航行路線決策幫助,可以得出以下結(jié)論。

    (1)減少總用時(shí)時(shí)間。通過(guò)算例將求解的路線與其他路線的最優(yōu)結(jié)果進(jìn)行對(duì)比,結(jié)果表明在此次算例中最佳航線比其他航線最優(yōu)結(jié)果最高可減少用時(shí)6.4%,證明了船舶動(dòng)態(tài)航線誘導(dǎo)的更優(yōu)性和模型的可行性。

    (2)提高模型適用性。充分考慮了海上船舶航行的特點(diǎn),在建模時(shí)考慮施工水域的影響、海風(fēng)和波浪的動(dòng)態(tài)變化影響以及沉點(diǎn)的影響,對(duì)路段行駛時(shí)間求解時(shí)對(duì)BPR模型進(jìn)行改進(jìn)并結(jié)合排隊(duì)論的理論,使得模型比以往更加適用于船舶。

    當(dāng)前模型與實(shí)際情況還有些區(qū)別,缺乏對(duì)突發(fā)情況的考慮,風(fēng)浪對(duì)船舶的影響的計(jì)算也還不夠精確,這是今后改進(jìn)的方向。

    猜你喜歡
    航行航線路段
    冬奧車道都有哪些相關(guān)路段如何正確通行
    到慧骃國(guó)的航行
    部、省、路段監(jiān)測(cè)運(yùn)維聯(lián)動(dòng)協(xié)同探討
    A Survey of Evolutionary Algorithms for Multi-Objective Optimization Problems With Irregular Pareto Fronts
    (21)新航線
    基于XGBOOST算法的擁堵路段短時(shí)交通流量預(yù)測(cè)
    小舟在河上航行
    航行
    青年歌聲(2017年6期)2017-03-13 00:57:56
    太空新航線
    太空探索(2016年5期)2016-07-12 15:17:58
    太空新航線
    太空探索(2016年6期)2016-07-10 12:09:06
    日本一二三区视频观看| 婷婷六月久久综合丁香| 精品久久久久久久人妻蜜臀av| 老女人水多毛片| 高清毛片免费看| 免费看av在线观看网站| 国产精品久久久久久精品电影小说 | 国产高清视频在线观看网站| 草草在线视频免费看| 欧美bdsm另类| 久久99蜜桃精品久久| 亚洲欧美成人综合另类久久久 | 又黄又爽又刺激的免费视频.| 两性午夜刺激爽爽歪歪视频在线观看| 哪个播放器可以免费观看大片| 高清av免费在线| 亚洲欧美清纯卡通| 99久久人妻综合| 麻豆成人av视频| 波多野结衣高清无吗| 久久精品影院6| 欧美激情在线99| 一边摸一边抽搐一进一小说| 观看美女的网站| 免费看美女性在线毛片视频| 偷拍熟女少妇极品色| 一边摸一边抽搐一进一小说| 国产淫片久久久久久久久| 久久99热6这里只有精品| 91精品国产九色| 99热这里只有是精品在线观看| 一级黄片播放器| 在线播放无遮挡| 亚洲图色成人| 亚洲在线自拍视频| 精品人妻一区二区三区麻豆| 免费黄网站久久成人精品| 桃色一区二区三区在线观看| 又爽又黄无遮挡网站| 99热这里只有精品一区| 国产大屁股一区二区在线视频| 又黄又爽又刺激的免费视频.| 青春草视频在线免费观看| 天天一区二区日本电影三级| 久久欧美精品欧美久久欧美| 高清日韩中文字幕在线| 欧美日韩在线观看h| 麻豆国产97在线/欧美| 精品欧美国产一区二区三| 国产乱人偷精品视频| 国产免费男女视频| 国产精品一区二区性色av| 国内少妇人妻偷人精品xxx网站| 国产精品久久久久久久电影| 成人亚洲欧美一区二区av| 日韩av不卡免费在线播放| 亚洲美女视频黄频| 日本午夜av视频| 亚洲av福利一区| 久久99热这里只频精品6学生 | 成人性生交大片免费视频hd| 亚洲国产色片| 久久精品国产99精品国产亚洲性色| 麻豆乱淫一区二区| 久热久热在线精品观看| 亚洲人成网站在线播| 亚洲人成网站高清观看| 成人一区二区视频在线观看| 国产精品不卡视频一区二区| 久久精品国产亚洲av涩爱| 国产精品一区二区在线观看99 | 国产激情偷乱视频一区二区| 麻豆一二三区av精品| 日韩在线高清观看一区二区三区| 人人妻人人澡人人爽人人夜夜 | 亚洲内射少妇av| 乱系列少妇在线播放| 亚洲av中文字字幕乱码综合| 男人舔奶头视频| 欧美性猛交黑人性爽| 一个人观看的视频www高清免费观看| 我要看日韩黄色一级片| 一个人免费在线观看电影| 国产成人a∨麻豆精品| 国产淫片久久久久久久久| 内地一区二区视频在线| 中文字幕久久专区| 久久精品91蜜桃| 美女大奶头视频| 少妇的逼好多水| 性插视频无遮挡在线免费观看| 非洲黑人性xxxx精品又粗又长| 99热这里只有精品一区| 免费在线观看成人毛片| 亚洲精品国产av成人精品| 一个人免费在线观看电影| 亚洲国产最新在线播放| 精品久久久久久久末码| 成人漫画全彩无遮挡| 亚洲伊人久久精品综合 | 神马国产精品三级电影在线观看| 国产精品国产高清国产av| 岛国毛片在线播放| 国产成人freesex在线| 国产伦一二天堂av在线观看| 日韩欧美在线乱码| 免费观看精品视频网站| 观看免费一级毛片| 日本三级黄在线观看| 国产精品蜜桃在线观看| 国产亚洲一区二区精品| 国产精品一二三区在线看| 最近2019中文字幕mv第一页| 少妇熟女欧美另类| 三级毛片av免费| 午夜激情欧美在线| 狂野欧美白嫩少妇大欣赏| 午夜精品在线福利| 一级毛片我不卡| 69人妻影院| 亚洲国产欧美人成| 国产一区二区在线av高清观看| 久久精品综合一区二区三区| 亚洲一级一片aⅴ在线观看| 亚洲av.av天堂| 日韩制服骚丝袜av| 久久久成人免费电影| av免费在线看不卡| 日韩欧美三级三区| 老司机影院成人| 深爱激情五月婷婷| 噜噜噜噜噜久久久久久91| 久久久精品欧美日韩精品| 国国产精品蜜臀av免费| 国产亚洲91精品色在线| 七月丁香在线播放| 日韩精品有码人妻一区| 国产真实伦视频高清在线观看| 亚洲自拍偷在线| 美女大奶头视频| 黄片无遮挡物在线观看| 亚洲成人中文字幕在线播放| 最新中文字幕久久久久| 男女下面进入的视频免费午夜| 深夜a级毛片| 国产精品无大码| 国产精品国产三级专区第一集| 全区人妻精品视频| 免费无遮挡裸体视频| 久久亚洲国产成人精品v| 国产又色又爽无遮挡免| 午夜免费男女啪啪视频观看| 两性午夜刺激爽爽歪歪视频在线观看| 国产精品一区二区三区四区免费观看| 波野结衣二区三区在线| 全区人妻精品视频| 国产爱豆传媒在线观看| 日韩av不卡免费在线播放| 久久99精品国语久久久| 97人妻精品一区二区三区麻豆| 久久久久精品久久久久真实原创| 亚洲av男天堂| 97超视频在线观看视频| 精品久久久久久久久av| 久久综合国产亚洲精品| 日日摸夜夜添夜夜爱| av国产久精品久网站免费入址| 亚洲欧美清纯卡通| 精品欧美国产一区二区三| av女优亚洲男人天堂| 精品一区二区三区视频在线| 国产人妻一区二区三区在| 日韩一区二区三区影片| 18禁动态无遮挡网站| 老师上课跳d突然被开到最大视频| 啦啦啦观看免费观看视频高清| 国产在线男女| 日韩,欧美,国产一区二区三区 | 有码 亚洲区| 日韩强制内射视频| 欧美zozozo另类| 国产视频内射| 日韩在线高清观看一区二区三区| 午夜福利高清视频| 亚洲18禁久久av| 久久久久九九精品影院| 色尼玛亚洲综合影院| 国产毛片a区久久久久| 亚洲精品久久久久久婷婷小说 | 免费在线观看成人毛片| 可以在线观看毛片的网站| 美女脱内裤让男人舔精品视频| 天天躁夜夜躁狠狠久久av| 日韩大片免费观看网站 | 国产一级毛片在线| 一级毛片aaaaaa免费看小| 免费无遮挡裸体视频| 白带黄色成豆腐渣| 欧美三级亚洲精品| 日本免费一区二区三区高清不卡| 国内少妇人妻偷人精品xxx网站| 汤姆久久久久久久影院中文字幕 | 欧美变态另类bdsm刘玥| 在线播放无遮挡| 久久久久久久久中文| 久久鲁丝午夜福利片| 波野结衣二区三区在线| 又粗又爽又猛毛片免费看| www日本黄色视频网| 色综合亚洲欧美另类图片| 国产av不卡久久| 搡老妇女老女人老熟妇| 久久久久九九精品影院| 国产单亲对白刺激| 真实男女啪啪啪动态图| kizo精华| av天堂中文字幕网| 亚洲精品亚洲一区二区| 亚洲伊人久久精品综合 | 2021少妇久久久久久久久久久| 成年免费大片在线观看| 最近最新中文字幕大全电影3| 国产激情偷乱视频一区二区| 三级经典国产精品| 欧美精品一区二区大全| 精华霜和精华液先用哪个| 床上黄色一级片| 国产精品久久久久久av不卡| 国产中年淑女户外野战色| 中文在线观看免费www的网站| 免费一级毛片在线播放高清视频| 两性午夜刺激爽爽歪歪视频在线观看| 国产三级中文精品| 人妻少妇偷人精品九色| 成年av动漫网址| 国产精品乱码一区二三区的特点| 亚洲av中文av极速乱| 亚洲av福利一区| 亚洲欧美清纯卡通| 国产在视频线在精品| 亚洲性久久影院| 永久免费av网站大全| 国产伦精品一区二区三区视频9| 狂野欧美激情性xxxx在线观看| 男人的好看免费观看在线视频| 久久精品久久久久久久性| 日韩视频在线欧美| h日本视频在线播放| 国产亚洲5aaaaa淫片| 亚洲熟妇中文字幕五十中出| 日日摸夜夜添夜夜爱| 99热全是精品| 国产成人a区在线观看| 精品酒店卫生间| 人人妻人人澡欧美一区二区| 国产欧美另类精品又又久久亚洲欧美| 一区二区三区四区激情视频| 免费在线观看成人毛片| 欧美97在线视频| 插逼视频在线观看| 国产成人免费观看mmmm| 能在线免费看毛片的网站| 在线a可以看的网站| 国产成人免费观看mmmm| 国产精品熟女久久久久浪| 中文字幕免费在线视频6| 国产淫语在线视频| 国产精品一区二区三区四区免费观看| 永久免费av网站大全| 国内揄拍国产精品人妻在线| 欧美性猛交黑人性爽| 能在线免费看毛片的网站| 99热这里只有是精品在线观看| 观看免费一级毛片| 麻豆成人av视频| 国产午夜福利久久久久久| 久久精品综合一区二区三区| 建设人人有责人人尽责人人享有的 | 亚洲欧美日韩高清专用| 午夜免费男女啪啪视频观看| 国产不卡一卡二| 性色avwww在线观看| 99久久精品一区二区三区| 亚洲自拍偷在线| av卡一久久| 成人综合一区亚洲| 国产精品国产三级专区第一集| 99久久无色码亚洲精品果冻| 欧美97在线视频| 亚洲欧美精品专区久久| 久久韩国三级中文字幕| 中文在线观看免费www的网站| 中文天堂在线官网| 超碰av人人做人人爽久久| 欧美极品一区二区三区四区| 国产成人精品久久久久久| 深爱激情五月婷婷| 观看免费一级毛片| 亚洲成人精品中文字幕电影| 嫩草影院新地址| 亚洲久久久久久中文字幕| 在线天堂最新版资源| 一级毛片电影观看 | 久久亚洲国产成人精品v| 最近最新中文字幕免费大全7| 午夜激情福利司机影院| 波野结衣二区三区在线| 亚洲自偷自拍三级| 成人性生交大片免费视频hd| 欧美97在线视频| 高清午夜精品一区二区三区| 久久人人爽人人爽人人片va| 欧美丝袜亚洲另类| 成人高潮视频无遮挡免费网站| 美女黄网站色视频| 人妻少妇偷人精品九色| 在线播放无遮挡| 日韩欧美精品v在线| 最近中文字幕高清免费大全6| 午夜福利成人在线免费观看| 精品国产一区二区三区久久久樱花 | 性色avwww在线观看| 欧美不卡视频在线免费观看| 一二三四中文在线观看免费高清| 看免费成人av毛片| av在线播放精品| 麻豆成人av视频| 淫秽高清视频在线观看| 人人妻人人看人人澡| 变态另类丝袜制服| 日韩一区二区视频免费看| 国产精品国产三级国产av玫瑰| 亚洲图色成人| 国产成人91sexporn| 一卡2卡三卡四卡精品乱码亚洲| 亚洲成人精品中文字幕电影| 国产单亲对白刺激| 免费观看a级毛片全部| 日韩一区二区三区影片| 亚洲av中文av极速乱| 人妻少妇偷人精品九色| 边亲边吃奶的免费视频| 国产人妻一区二区三区在| 亚洲中文字幕一区二区三区有码在线看| 一区二区三区免费毛片| av福利片在线观看| 能在线免费看毛片的网站| 人人妻人人澡人人爽人人夜夜 | 久久久色成人| 亚洲欧美精品自产自拍| 99在线视频只有这里精品首页| 日本黄色视频三级网站网址| 男人的好看免费观看在线视频| 成人亚洲欧美一区二区av| 国产精品久久久久久av不卡| 夫妻性生交免费视频一级片| 亚洲高清免费不卡视频| 夫妻性生交免费视频一级片| 亚洲va在线va天堂va国产| 国语自产精品视频在线第100页| 99国产精品一区二区蜜桃av| 国产成人freesex在线| 少妇裸体淫交视频免费看高清| 中国美白少妇内射xxxbb| 日韩强制内射视频| 国产伦在线观看视频一区| 欧美性感艳星| 不卡视频在线观看欧美| 最近中文字幕高清免费大全6| 熟妇人妻久久中文字幕3abv| 岛国在线免费视频观看| 国语自产精品视频在线第100页| 91久久精品国产一区二区成人| 久久精品国产亚洲av涩爱| 最近视频中文字幕2019在线8| 国产 一区 欧美 日韩| 夜夜看夜夜爽夜夜摸| 黄片无遮挡物在线观看| 老司机影院毛片| 日本色播在线视频| 99在线视频只有这里精品首页| 国产又色又爽无遮挡免| 天天躁夜夜躁狠狠久久av| 久久国产乱子免费精品| 精品一区二区三区视频在线| 小说图片视频综合网站| 一级二级三级毛片免费看| 别揉我奶头 嗯啊视频| 亚洲精品aⅴ在线观看| 欧美xxxx黑人xx丫x性爽| 麻豆成人av视频| 色综合站精品国产| 亚洲不卡免费看| 26uuu在线亚洲综合色| 午夜激情福利司机影院| 国产熟女欧美一区二区| a级毛片免费高清观看在线播放| 亚洲欧美一区二区三区国产| 白带黄色成豆腐渣| 亚洲精品国产成人久久av| 成人鲁丝片一二三区免费| 边亲边吃奶的免费视频| 久久精品综合一区二区三区| 18禁动态无遮挡网站| 午夜福利高清视频| 色播亚洲综合网| 人人妻人人看人人澡| 亚洲国产精品专区欧美| a级毛色黄片| 国产一级毛片七仙女欲春2| 成人午夜高清在线视频| 国产探花在线观看一区二区| 97人妻精品一区二区三区麻豆| 能在线免费观看的黄片| 色网站视频免费| 亚洲av一区综合| 国产在视频线在精品| 日日啪夜夜撸| 亚洲国产色片| 精品不卡国产一区二区三区| 噜噜噜噜噜久久久久久91| 日本三级黄在线观看| av播播在线观看一区| 亚洲四区av| 国产精品综合久久久久久久免费| 久久人人爽人人片av| 少妇人妻一区二区三区视频| 亚洲精品一区蜜桃| 中文精品一卡2卡3卡4更新| 一级毛片aaaaaa免费看小| 国产亚洲最大av| 国产成人一区二区在线| 嫩草影院入口| 人妻夜夜爽99麻豆av| 中文字幕精品亚洲无线码一区| av播播在线观看一区| 亚洲av电影不卡..在线观看| 美女国产视频在线观看| 秋霞伦理黄片| 一夜夜www| 国产精品一区二区在线观看99 | 美女内射精品一级片tv| av在线蜜桃| 精品无人区乱码1区二区| 91av网一区二区| 精品99又大又爽又粗少妇毛片| 国产精品一及| 精品久久久久久久人妻蜜臀av| 免费看美女性在线毛片视频| 夫妻性生交免费视频一级片| 国产真实乱freesex| 亚洲精品,欧美精品| 极品教师在线视频| 丰满人妻一区二区三区视频av| 精品久久国产蜜桃| 国产成人一区二区在线| 18禁在线无遮挡免费观看视频| 亚洲中文字幕一区二区三区有码在线看| 水蜜桃什么品种好| 亚洲精品成人久久久久久| 精品久久久久久久久久久久久| 国产av不卡久久| 中文资源天堂在线| 天天躁日日操中文字幕| 大话2 男鬼变身卡| 日本黄色片子视频| 色综合亚洲欧美另类图片| 中文字幕久久专区| 久久综合国产亚洲精品| 欧美日韩综合久久久久久| 1024手机看黄色片| 亚洲乱码一区二区免费版| 高清毛片免费看| 国产亚洲91精品色在线| 国产精品99久久久久久久久| 国语对白做爰xxxⅹ性视频网站| 亚洲精品aⅴ在线观看| 欧美丝袜亚洲另类| 亚洲怡红院男人天堂| 亚洲人成网站高清观看| 欧美性感艳星| 欧美区成人在线视频| 日韩 亚洲 欧美在线| 97人妻精品一区二区三区麻豆| 日韩欧美国产在线观看| 国产视频首页在线观看| 亚洲经典国产精华液单| 午夜爱爱视频在线播放| 亚洲国产精品sss在线观看| 中文字幕人妻熟人妻熟丝袜美| 亚洲欧洲日产国产| 伊人久久精品亚洲午夜| 午夜免费激情av| 久久欧美精品欧美久久欧美| 欧美成人午夜免费资源| АⅤ资源中文在线天堂| 日本熟妇午夜| 男人和女人高潮做爰伦理| 欧美日本亚洲视频在线播放| 99久久人妻综合| 男女边吃奶边做爰视频| 色综合色国产| 人妻少妇偷人精品九色| 欧美bdsm另类| 欧美潮喷喷水| 丝袜美腿在线中文| 国产精品综合久久久久久久免费| 女的被弄到高潮叫床怎么办| 国产亚洲91精品色在线| 国产精品乱码一区二三区的特点| 99久国产av精品国产电影| 午夜福利在线在线| 中文字幕人妻熟人妻熟丝袜美| 成人特级av手机在线观看| 国产精品.久久久| 亚洲国产精品成人久久小说| 国产一区二区在线观看日韩| 精品国内亚洲2022精品成人| 国产黄片视频在线免费观看| 亚洲欧美精品自产自拍| 国产一级毛片在线| 中文精品一卡2卡3卡4更新| 特级一级黄色大片| 午夜精品在线福利| 美女大奶头视频| 一级爰片在线观看| 一级二级三级毛片免费看| 欧美区成人在线视频| 日韩一区二区三区影片| 我的女老师完整版在线观看| 中文资源天堂在线| 亚洲欧美清纯卡通| 国产精品无大码| 亚洲人成网站高清观看| 日本av手机在线免费观看| 你懂的网址亚洲精品在线观看 | 亚洲国产精品专区欧美| 97超碰精品成人国产| 日韩欧美精品免费久久| 国产三级在线视频| 亚洲av不卡在线观看| 国产伦精品一区二区三区视频9| 精品久久久久久久人妻蜜臀av| 国产精品久久久久久av不卡| 久久久亚洲精品成人影院| 亚洲国产精品久久男人天堂| 91久久精品国产一区二区三区| 一个人看的www免费观看视频| 熟女人妻精品中文字幕| 国产伦精品一区二区三区视频9| 少妇的逼水好多| av国产久精品久网站免费入址| 久久久午夜欧美精品| 国产av一区在线观看免费| 舔av片在线| 国产精品久久久久久久电影| 国产单亲对白刺激| 高清av免费在线| 99久久九九国产精品国产免费| 简卡轻食公司| 草草在线视频免费看| 18+在线观看网站| 一区二区三区乱码不卡18| 亚洲真实伦在线观看| 十八禁国产超污无遮挡网站| 在线观看美女被高潮喷水网站| 大话2 男鬼变身卡| 一卡2卡三卡四卡精品乱码亚洲| 国国产精品蜜臀av免费| 成人亚洲欧美一区二区av| 国产亚洲一区二区精品| 国产精品一区二区三区四区久久| 不卡视频在线观看欧美| 能在线免费观看的黄片| 亚洲伊人久久精品综合 | 少妇的逼好多水| 国产黄片美女视频| 三级国产精品片| 国产精品人妻久久久影院| 两性午夜刺激爽爽歪歪视频在线观看| 国产淫片久久久久久久久| 久久亚洲精品不卡| 亚洲国产日韩欧美精品在线观看| 亚洲四区av| 午夜a级毛片| 欧美性感艳星| 国产一级毛片七仙女欲春2| 欧美极品一区二区三区四区| 91久久精品国产一区二区三区| 中文亚洲av片在线观看爽| 久热久热在线精品观看| 一级毛片电影观看 | 一二三四中文在线观看免费高清| 久久精品人妻少妇| 国产色爽女视频免费观看| 一区二区三区四区激情视频| 国产三级中文精品| 国产精品永久免费网站| 综合色丁香网| 22中文网久久字幕| 老女人水多毛片| 岛国在线免费视频观看| 波多野结衣高清无吗| 一级爰片在线观看| 国产精品久久久久久av不卡| 人人妻人人澡欧美一区二区| 青春草视频在线免费观看| 午夜福利在线观看吧| 深夜a级毛片| 国产在线男女| 午夜福利视频1000在线观看| 久久久午夜欧美精品| 国产精品久久电影中文字幕| 五月玫瑰六月丁香| 99热全是精品|