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

    遵循COLREGs的USV多階段自適應(yīng)ph-Radau偽譜避障軌跡規(guī)劃

    2018-10-15 11:27:18戴懷志蒲紅紅
    指揮與控制學(xué)報 2018年3期
    關(guān)鍵詞:偽譜航路障礙物

    戴懷志 蒲紅紅

    障礙物規(guī)避策略對無人水面航行器(Unmanned Surface Vehicle,USV)順利完成自主航行與作業(yè)十分重要,而為實現(xiàn)有效且協(xié)調(diào)的避障,將USV作為洋流中航行船舶的一種,遵循國際海上避碰規(guī)則公約(International Regulations for Preventing Collisions at Sea,COLREGs)是必要的[1].

    針對USV避障問題,多數(shù)學(xué)者均采用遵循COLREGs的速度避障法[2?4].文獻[2]針對未知動態(tài)環(huán)境,線性化處理相對速度,并將其作為約束條件整合到混合整數(shù)線性規(guī)劃器中,選擇多目標(biāo)函數(shù)作為優(yōu)化函數(shù),提出一種遵循COLREGs的在線路徑規(guī)劃算法.文獻[3]則利用動態(tài)障礙物與USV的相對運動狀態(tài),提出一種遵循COLREGs的粒子群優(yōu)化(Particle Swarm Optimization,PSO)障礙物規(guī)避策略,實時解算規(guī)避障礙物的最優(yōu)航行和航速改變量.而文獻[4]在遵循COLREGs的前提下,靜態(tài)和動態(tài)障礙物的規(guī)避問題,提出基于PSO的動、靜態(tài)障礙物規(guī)避分層路徑規(guī)劃方法.

    然而,上述文獻中USV避障策略無法從USV運動模型出發(fā)且未充分考慮航行機動性能限制等約束條件,難以保證所規(guī)劃路徑的可行性,本質(zhì)上歸屬于路徑規(guī)劃范疇.而基于偽譜法的軌跡優(yōu)化策略是一類活躍在航天器領(lǐng)域的軌跡規(guī)劃數(shù)值解法,以最優(yōu)控制理論為通用框架,可方便地加入各種約束條件,且求解精度高,收斂速度快.此外,航天器軌跡優(yōu)化策略可分為偽譜法和多階段偽譜法[5?7]求解兩類軌跡優(yōu)化策略.文獻[8]從六自由度運動模型出發(fā),充分考慮飛行機動性能限制和三維空間靜態(tài)障礙物約束,建立連續(xù)時間避障最優(yōu)軌跡規(guī)劃模型,給出基于Gauss偽譜法的無人直升機避障軌跡規(guī)劃策略.文獻[9]則引入內(nèi)點約束將突防軌跡分段,并綜合考慮求解精度和速率,結(jié)合Radau偽譜法和自適應(yīng)hp節(jié)點配置法,求解含靜態(tài)禁飛區(qū)和必經(jīng)航路點約束的多階段再入飛行器軌跡快速優(yōu)化問題.值得注意的是,以上基于偽譜法的航天器軌跡優(yōu)化策略均為針對靜態(tài)障礙物和靜態(tài)必經(jīng)路徑點的處理.

    因此,基于以上分析,本文針對USV對動、靜態(tài)障礙物的機動避障問題,在遵循COLREGs的前提下,提出一種基于自適應(yīng)ph-Radau偽譜法的USV多階段避障軌跡規(guī)劃算法.該算法充分考慮了機動性能、狀態(tài)方程以及障礙物(包括動態(tài)和靜態(tài))等約束條件,能夠保證所規(guī)劃軌跡的可行性;COLREGs的遵守和內(nèi)點約束的引入,使得USV能同時對動態(tài)和靜態(tài)障礙物有效且協(xié)調(diào)地機動避讓;相比Gauss偽譜法,所設(shè)計算法具備更高的求解精度和速率.

    1 國際海上避碰規(guī)則公約與避讓策略

    1.1 國際海上避碰規(guī)則公約

    目前還未有相關(guān)法律對USV的航行行為作出明確規(guī)定,但較為合理的方案就是將USV作為船舶的一種,遵循國際海事組織制定的COLREGs.

    該公約第13~15條規(guī)定了碰撞局面及相關(guān)避讓規(guī)則.以USV航向為基準(zhǔn),根據(jù)障礙物與USV的航向角度差值,定義圖1所示的碰撞局面和如下避讓規(guī)則[3]:

    1)追越:若障礙物與USV的航向角度差值在[0?,45?)和[315?,360?)之間,當(dāng)接近危險距離時,要求USV從被超越障礙物左舷處安全通過.

    2)正面相遇:若障礙物與USV的航向角度差值在(165?,195?]之間,當(dāng)接近危險距離時,要求USV從障礙物右舷處安全通過.

    3)交叉相遇:若障礙物與USV的航向角度差值在[45?,165?]之間,障礙物與USV 右交叉相遇,當(dāng)接近危險距離時,要求USV從障礙物右舷處安全通過;若航向角度差值在(195?,315?)之間,障礙物與USV左交叉相遇,當(dāng)接近危險距離時,要求USV從障礙物左舷處安全通過.

    1.2 避讓策略

    當(dāng)USV在攜帶傳感器的探測范圍內(nèi)檢測到未知的動、靜態(tài)障礙物并存在碰撞危險時,依據(jù)COLREGs判定碰撞局面,在障礙物和USV保持當(dāng)前位姿不變的假定條件下估算相撞時刻,設(shè)定合適的避讓航路點和禁航區(qū)約束.

    其中,避讓航路點可看作是為了提高海事交通安全性,即有效且協(xié)調(diào)地安全避障,讓USV遵循COLREGs而設(shè)定的必經(jīng)地理位置.而為避免避讓策略的復(fù)雜性,僅在動態(tài)障礙物規(guī)避時設(shè)置避讓航路點即可.此外,為保證USV在安全避讓障礙物后能“優(yōu)雅”且快速地返回原始機動航行路線,需限制避讓航路點處的其他條件(如時間、姿態(tài)),避讓航路點約束可設(shè)置如下:

    式中,tl為要求USV到達避讓航路點處的估算時刻,即時間約束;(xl,yl,ψl)為第l個避讓航路點的位置和航向;ψlmin和ψlmax為航向ψl的邊界值;lN為避讓航路點數(shù)目.

    為保障避讓策略的可靠性,將未知靜態(tài)障礙物看作是一類特殊的動態(tài)障礙物并考慮避讓過程中仍存在碰撞危險,要求USV始終與障礙物保持一定的安全距離,即設(shè)置禁航區(qū)約束.將障礙物近似為一個速度和航向均為常值,半徑為Rj的圓域,j=1,2,···jN,jN為障礙物數(shù)目.假設(shè)初始時刻t0=0圓心坐標(biāo)為(xj0,yj0),移動速度為Vj,航向為ψj,t時刻的障礙物可描述為

    則禁航區(qū)約束為

    2 避障最優(yōu)軌跡規(guī)劃模型

    2.1 約束條件

    1)狀態(tài)方程約束

    水平面三自由度欠驅(qū)動USV運動模型,即狀態(tài)方程[10]約束如下:

    式中,m11=m?X˙u,m22=m?Y˙v,m33=Iz?N˙r;C13=?m22v,C23=m11u;d11=?Xu,d22=?Xv,d33=?Nr;τ=[Tu0Tr]T為控制輸入向量,τu為前進推力,τr為偏航力矩.注意到,式(4)中僅存在作用在縱蕩和艏搖兩個運動自由度上的控制力/力矩τu和τr,橫蕩運動自由度上不存在控制力/力矩,系統(tǒng)的控制輸入少于運動自由度個數(shù),故該三自由度USV為典型的欠驅(qū)動系統(tǒng).

    2)機動性能約束

    為保證所使用運動模型的可靠性,根據(jù)航行器的機動性能對狀態(tài)變量作如下約束:

    式中,umax、vmax和ψmax均為非負常值.此外,考慮到控制變量受其執(zhí)行機構(gòu)性能與姿態(tài)控制系統(tǒng)的限制,為避免USV在航行過程中失控,對其前進推力和偏航力矩作如下約束:

    式中,Tumax和Trmax均為非負常值.

    3)禁航區(qū)約束

    遵循COLREGs設(shè)置禁航區(qū)和避讓航路點約束,其中禁航區(qū)約束為式(3).

    4)邊界條件約束

    初始時刻t0=0的USV位姿為x x x(t0),指定終端時刻tf的避障結(jié)束位置(xf,yf)為終端目標(biāo),即

    為提高避障可靠性,引入終端目標(biāo)處的其他條件約束(如速度、姿態(tài)),即

    式中,ufmin、ufmax、ψfmin和ψfmax均為常值.

    2.2 性能指標(biāo)

    考慮到大洋中自主航行的USV攜帶能量有限,且能量補給困難,可選擇已到達終端目標(biāo)處的能耗最少為性能指標(biāo),即

    當(dāng)然,依據(jù)USV的具體作業(yè)與任務(wù)需求,該性能指標(biāo)也可選擇為航行時間最短、航程最短等.

    2.3 多階段避障最優(yōu)軌跡規(guī)劃模型

    將避障軌跡在避讓航路點處分段,并引入內(nèi)點約束,將避讓航路點轉(zhuǎn)化為內(nèi)點約束,進而可借鑒多階段偽譜法軌跡優(yōu)化策略,建立USV多階段跡規(guī)劃模型.

    給定避讓航路點數(shù)目lN,將原始時間區(qū)域I=[t0,tf]劃分為P段,P=lN+1,階段Ip=[tp?1,tp],p=1,2,···,P滿足:且于是,含有避讓航路點約束的USV避障軌跡規(guī)劃問題可描述為多階段軌跡規(guī)劃問題,即尋找控制變量u u u(p)(t)∈Rnu,最小化Bolza型性能指標(biāo):

    滿足約束條件:

    本文避讓航路點約束式(1)轉(zhuǎn)化為如下內(nèi)點約束:

    式中,ψl滿足ψlmin≤ψl≤ψlmax.

    3 自適應(yīng)ph-Radau偽譜法

    在階段Ip=[tp?1,tp],p=1,2,···,P內(nèi),以偽譜法,如Gauss、Radau偽譜法離散化求解策略,可將連續(xù)時間最優(yōu)控制式(11)~式(14)轉(zhuǎn)化為含有帶有代數(shù)方程約束的大規(guī)模非線性規(guī)劃(Nonlinear Porgram,NLP)問題,借助序列二次規(guī)劃等數(shù)值方法,可獲得原連續(xù)時間最優(yōu)控制問題的近似數(shù)值解[11?12].

    然而,多數(shù)實際問題均為非光滑問題,固定網(wǎng)格個數(shù)和配點數(shù)配置方式的偽譜法無法捕捉非光滑最優(yōu)控制問題的狀態(tài)或控制的不連續(xù)性;且求解精度與全局插值多項式維數(shù)(即配點數(shù))有關(guān),配點數(shù)的增加在提高求解精度的同時也增加計算代價[13].h法和p法是兩典型的網(wǎng)格細化算法,可有效解決NLP維度高與求解速率低的矛盾問題.因此,本文在Radau偽譜法的基礎(chǔ)上,結(jié)合自適應(yīng)網(wǎng)格劃分與節(jié)點個數(shù)分配[14?15],即基于自適應(yīng)ph-Rddau偽譜法求解多階段避障軌跡規(guī)劃問題.

    其中,自適應(yīng)ph配置方法為:當(dāng)任意網(wǎng)格區(qū)間的最大相對誤差,k=1,2,···,K大于設(shè)定的誤差門限值ε時,自適應(yīng)進行網(wǎng)格劃分或增加Legendre Gauss Radau(LGR)節(jié)點個數(shù).

    1)誤差判定準(zhǔn)則

    假設(shè)在階段Ip的第k個網(wǎng)格區(qū)間[tk?1,tk]內(nèi),Nk個 LGR 配點和末端點=1構(gòu)成了該網(wǎng)格區(qū)間的插值節(jié)點,而連續(xù)時間最優(yōu)控制式(11)~式(14)及其狀態(tài)變量、控制變量的Radau偽譜法離散化過程見文獻[12?14].考慮到增加配點個數(shù)可提高求解精度,因而設(shè)定為網(wǎng)格區(qū)間[tk?1,tk]內(nèi)的Mk個LGR配點,且Mk=Nk+1,=?1,=1.取Mk個點作為誤差判定準(zhǔn)則的采樣點,則近似狀態(tài)?變量·和近似控?制變·量在采樣點上可表示為、,i=1,2,···,Mk.以? 為插·值節(jié)點,?得到·近似狀態(tài)變量和控制變量和,則狀態(tài)方程在LGR配點=1,2,···,Mk上可近似為

    式 中,=j=1,2,···,Mk為網(wǎng)格區(qū)間[tk?1,tk]的Mk×Mk維LGR積分權(quán)重矩陣.則第i個狀態(tài)在采樣點上的絕對誤差和相對誤差可表示為

    網(wǎng)格區(qū)間[tk?1,tk]中最大相對誤差為

    2)自適應(yīng)ph網(wǎng)格配置方法

    令Nmin和Nmax分別為每個網(wǎng)格區(qū)間內(nèi)的最小和最大LGR配點個數(shù).假定網(wǎng)格區(qū)間[tk?1,tk]內(nèi)LGR配點個數(shù)為Nk,Nk≤Nmax,最大相對誤差大于誤差門限值ε,則令

    增加配點個數(shù),并取整為

    當(dāng)進行網(wǎng)格細化時,要求細化后所有網(wǎng)格區(qū)間的配點個數(shù)之和與一致,且每個網(wǎng)格區(qū)間的配點個數(shù)均為Nmin,即將網(wǎng)格區(qū)間[tk?1,tk]進行Bq等分,等分子區(qū)間配點個數(shù)之和為BqNmin.等分子區(qū)間的個數(shù)Bq取

    值得注意的是,在該網(wǎng)格更新策略中,無論是增加配點個數(shù)還是進一步網(wǎng)格細化,配點總數(shù)均為?Nq.該網(wǎng)格更新策略可稱為自適應(yīng)ph網(wǎng)格更新策略,即先p后h的網(wǎng)格更新策略,當(dāng)區(qū)間的配點個數(shù)超過最大配點個數(shù)Nmax時,才對網(wǎng)格進行細化處理.而結(jié)合Radau偽譜法與自適應(yīng)ph網(wǎng)格更新策略的自適應(yīng)ph-Radau偽譜法具體迭代策略可參見文獻[15].

    4 仿真實驗與分析

    采用挪威科技大學(xué)試驗?zāi)P痛珻yberShipΠ[16]模型進行仿真實驗,驗證所設(shè)計算法的可行性與有效性,即是否能遵循COLREGs制定的避讓航路點和禁航區(qū)約束,成功規(guī)避動、靜態(tài)障礙物的同時,也能兼顧算法的求解精度與速率.

    4.1 避讓航路點約束驗證實驗

    選擇靜態(tài)和動態(tài)障礙物各一個進行避障實例驗證.為提高避障安全度,對靜、動態(tài)障礙物進行膨脹處理.障礙物運動參數(shù)如表1所示.

    表1 障礙物運動參數(shù)

    遵循COLREGs判定存在碰撞危險,且與動態(tài)障礙物的碰撞局面為正面相遇,避讓航路點數(shù)目lN=1,障礙物數(shù)目jN=2.其中,估算避讓時刻為tl=40 s,避讓航路點約束為

    而帶入障礙物運動參數(shù)表1和式(2)~式(3),即可獲得禁航區(qū)約束.

    此外,給定USV的邊界條件約束:

    USV的機動性能約束:

    給定自適應(yīng)ph-Radau偽譜法離散化求解策略中的參數(shù):

    式中,Mmax為最大迭代次數(shù);K0和N0分別表示階段內(nèi)網(wǎng)格區(qū)間初始個數(shù)和網(wǎng)格區(qū)間內(nèi)LGR配點個數(shù);Nmin和Nmax分別為網(wǎng)格區(qū)間內(nèi)最小和最大配點個數(shù).

    仿真結(jié)果如圖2~圖5所示.圖2為在最優(yōu)控制變量作用下USV最優(yōu)航行軌跡變化曲線,其中,黑色方框?qū)?yīng)初始時刻t0、避讓時刻tl和終端時刻tf對應(yīng)的USV具體位置;橙色實心圓對應(yīng)初始時刻t0的動態(tài)障礙物具體位置,灰色實心圓分別對應(yīng)避讓時刻tl和終端時刻tf的動態(tài)障礙物具體位置,黑色箭頭表示動態(tài)障礙物行進過程;黑色實心圓為靜態(tài)障礙物.圖3為禁航區(qū)變化曲線,即避障過程中USV與靜、動態(tài)障礙物之間距離的時間變化曲線.注意到,在禁航區(qū)約束中取膨脹2m后的障礙物半徑,而在實際繪制USV最優(yōu)航行軌跡變化曲線時取原障礙物半徑.圖4為最優(yōu)控制輸入變化曲線.圖5為USV最優(yōu)狀態(tài)變化曲線.可以看出:1)控制變量在整個避障過程中較為平滑,無較大跳變出現(xiàn),USV的初始狀態(tài)和終端狀態(tài)也均在約束范圍內(nèi),即機動性能約束條件得以滿足;2)USV準(zhǔn)確地經(jīng)過了避讓航路點,且禁航區(qū)約束的設(shè)置使得USV在整個航行過程中與動、靜態(tài)障礙物保持一定安全距離,即避讓航路點約束和禁航區(qū)約束條件得以滿足;3)USV的初始狀態(tài)和終端狀態(tài)均在約束范圍內(nèi),邊界條件約束得以滿足;4)由避讓航路點造成的兩階段避障軌跡和各狀態(tài)變量銜接處的變化均較為平滑,表明本章所設(shè)計的多階段避障軌跡規(guī)劃策略中內(nèi)點約束的有效性.

    此外,該自適應(yīng)ph-Radau偽譜法多階段避障軌跡規(guī)劃算法得到的網(wǎng)格區(qū)間最大相對誤差為5.2344×10?4,滿足誤差門限值.文獻[4]中的避障優(yōu)化軌跡與本文的優(yōu)化軌跡差別較大,這是由于基于PSO算法優(yōu)化避障軌跡過程中未從USV無人運載平臺的狀態(tài)方程出發(fā),而在實際航行過程中各運動自由度存在一定的耦合關(guān)系.

    4.2 禁航區(qū)約束驗證實驗

    保持上述仿真參數(shù)不變,驗證本章所設(shè)計算法是否能夠有效處理制定的禁航區(qū)約束.圖6給出了同時包含避讓航路點約束、禁航區(qū)約束和僅包含避讓航路點約束的USV航行軌跡變化曲線.圖7對應(yīng)兩種情況的禁航區(qū)變化曲線.

    其中,同時包含避讓航路點和禁航區(qū)約束的情況存在禁航區(qū)變化曲線,僅包含避讓航路點約束的情況不存在禁航區(qū)約束,但考慮到禁航區(qū)實質(zhì)即為USV與動、靜態(tài)障礙物之間的距離變化,故為了描述方便,仍用“禁航區(qū)變化曲線”描述僅包含避讓航路點約束情況下USV與障礙物之間距離的變化.此外,取原障礙物半徑進行USV最優(yōu)航行軌跡變化曲線的繪制,而在禁航區(qū)約束中取膨脹半徑為2 m后的障礙物半徑.

    可以看出:1)兩種情況下,USV均精確通過避讓航路點;2)僅考慮避讓航路點約束,即忽略避讓過程中存在的碰撞危險,與同時考慮避讓航路點、禁航區(qū)約束的情況下,USV的最優(yōu)航行軌跡是不同的.對于存在禁航區(qū)的情況,USV需要做“S”形機動繞過禁航區(qū),以達到指定目標(biāo)位置.

    因此,圖2~圖7的分析表明了將避讓航路點轉(zhuǎn)化為內(nèi)點,并考慮禁航區(qū)約束的多階段避障軌跡規(guī)劃策略的可行性.

    4.3 求解精度和速度驗證實驗

    為了更好地說明本章所設(shè)計算法的有效性,給出采用Gauss偽譜法離散化求解策略和本文求解策略的USV最優(yōu)航行軌跡對比,如圖8所示.表2對應(yīng)兩種情況下的部分仿真結(jié)果,其中LG點為Gauss偽譜法離散化求解策略所需的Lagendre-Gauss離散點.為便于比較兩種求解策略的求解精度,給出圖8中障礙物膨脹2 m后的USV最優(yōu)航跡局部放大圖,如圖9所示.其中,粉色區(qū)域?qū)?yīng)障礙物的膨脹部分.

    從圖8、圖9和表2可明顯地看出:1)在離散點個數(shù)近乎相同時,Gauss偽譜離散化求解策略下的靜態(tài)障礙物禁航區(qū)約束未被滿足,得到的最優(yōu)航行軌跡穿過了禁航區(qū);而自適應(yīng)ph-Radau偽譜法則根據(jù)設(shè)定的解的誤差判定準(zhǔn)則,自動檢測到軌跡精度的滿足情況,自適應(yīng)增加離散點個數(shù)或者進一步網(wǎng)格細化.3)自適應(yīng)ph-Radau偽譜法的仿真運行時間遠小于Gauss偽譜法.

    表2 部分仿真結(jié)果

    因此,本章自適應(yīng)ph-Radau偽譜法多階段避障軌跡規(guī)劃算法能夠較好處理避讓航路點和禁航區(qū)約束的同時,也能兼顧解的求解精度和速度.

    5 結(jié)論

    本文提出一種基于自適應(yīng)ph-Radau偽譜法的USV多階段避障軌跡規(guī)劃算法.相比基于PSO等速度避障法,該算法從USV運動模型出發(fā),充分考慮航行機動性能限制等約束條件,保證了所規(guī)劃避障軌跡的可行性.而遵循COLREGs設(shè)定的必然航路點和禁航區(qū)約束,以及內(nèi)點約束的引入,保證了USV能同時有效且協(xié)調(diào)地規(guī)避動態(tài)和靜態(tài)障礙物.此外,相比Gausss偽譜法,該算法能兼顧求解精度和速率.

    如何將該方法應(yīng)用于外界環(huán)境,如海風(fēng)、海浪、海流擾動下的USV避障軌跡規(guī)劃問題可作為進一步的研究方向.

    猜你喜歡
    偽譜航路障礙物
    基于實時航路的PFD和ND的仿真研究
    矩陣偽譜的新定位集及其在土壤生態(tài)系統(tǒng)的應(yīng)用
    高低翻越
    SelTrac?CBTC系統(tǒng)中非通信障礙物的設(shè)計和處理
    紊流環(huán)境下四維軌跡優(yōu)化的偽譜方法研究
    應(yīng)召反潛時無人機監(jiān)聽航路的規(guī)劃
    托勒密世界地圖與新航路的開辟
    偽譜法及其在飛行器軌跡優(yōu)化設(shè)計領(lǐng)域的應(yīng)用綜述*
    基于Event改進模型的交叉航路碰撞風(fēng)險評估
    基于Radau偽譜方法的軌跡優(yōu)化
    自拍偷自拍亚洲精品老妇| 我要看日韩黄色一级片| 亚洲国产最新在线播放| 在线亚洲精品国产二区图片欧美 | 亚洲国产精品专区欧美| 国产精品一区二区三区四区免费观看| 亚洲人与动物交配视频| 亚洲国产最新在线播放| 乱系列少妇在线播放| 美女cb高潮喷水在线观看| 能在线免费看毛片的网站| 久久精品久久精品一区二区三区| 亚洲国产欧美在线一区| 成人国产av品久久久| 少妇熟女欧美另类| 国产精品99久久99久久久不卡 | 精品久久久久久久久亚洲| 亚洲欧美一区二区三区黑人 | 91狼人影院| 晚上一个人看的免费电影| 97超视频在线观看视频| 波野结衣二区三区在线| 欧美成人午夜免费资源| av国产精品久久久久影院| 99热网站在线观看| 一级毛片黄色毛片免费观看视频| 夜夜爽夜夜爽视频| 国产成人aa在线观看| 18禁在线播放成人免费| 国产一区二区在线观看日韩| 欧美最新免费一区二区三区| 在线观看一区二区三区激情| 青春草国产在线视频| 热99国产精品久久久久久7| 精品久久久精品久久久| 国产 精品1| 亚洲国产色片| 亚洲成人一二三区av| 日本黄色片子视频| 麻豆精品久久久久久蜜桃| 一级黄片播放器| 在线免费观看不下载黄p国产| 国产69精品久久久久777片| 久久婷婷青草| 午夜福利高清视频| 97在线视频观看| 女性被躁到高潮视频| 日日啪夜夜爽| 全区人妻精品视频| 国产欧美另类精品又又久久亚洲欧美| 毛片女人毛片| 一级爰片在线观看| av天堂中文字幕网| av免费在线看不卡| 99久久精品国产国产毛片| 99久久人妻综合| 看免费成人av毛片| 草草在线视频免费看| 天美传媒精品一区二区| 极品教师在线视频| 伦理电影大哥的女人| 少妇人妻一区二区三区视频| 最近中文字幕高清免费大全6| 天堂俺去俺来也www色官网| 久久久久久久久久久丰满| 99热这里只有精品一区| 在线观看美女被高潮喷水网站| 国产精品三级大全| 女人十人毛片免费观看3o分钟| 一本一本综合久久| 亚洲国产成人一精品久久久| 亚洲图色成人| 80岁老熟妇乱子伦牲交| 大陆偷拍与自拍| 亚洲精品视频女| 久久99精品国语久久久| 一本—道久久a久久精品蜜桃钙片| 久久久久网色| av一本久久久久| 偷拍熟女少妇极品色| 97热精品久久久久久| 最近2019中文字幕mv第一页| 亚洲美女搞黄在线观看| 成人综合一区亚洲| 一个人看视频在线观看www免费| 免费观看在线日韩| 边亲边吃奶的免费视频| 国产真实伦视频高清在线观看| 色婷婷av一区二区三区视频| 国产国拍精品亚洲av在线观看| 国产精品一区二区在线不卡| 美女福利国产在线 | 成年人午夜在线观看视频| 99热这里只有是精品在线观看| av线在线观看网站| 久久精品熟女亚洲av麻豆精品| 一个人看视频在线观看www免费| 国产免费一级a男人的天堂| 国产精品人妻久久久影院| 婷婷色综合www| 最近最新中文字幕大全电影3| 免费av不卡在线播放| 亚洲美女视频黄频| 免费观看的影片在线观看| 最近最新中文字幕大全电影3| 国产精品久久久久久久电影| av免费在线看不卡| 欧美精品亚洲一区二区| 国产精品人妻久久久久久| 亚洲av欧美aⅴ国产| 亚洲av成人精品一二三区| av又黄又爽大尺度在线免费看| 国国产精品蜜臀av免费| 一二三四中文在线观看免费高清| 最近2019中文字幕mv第一页| 如何舔出高潮| 久久毛片免费看一区二区三区| 看免费成人av毛片| 亚洲av国产av综合av卡| 九九久久精品国产亚洲av麻豆| 欧美精品国产亚洲| 嫩草影院入口| 成人二区视频| 亚洲第一av免费看| 网址你懂的国产日韩在线| 午夜福利网站1000一区二区三区| 国产伦精品一区二区三区视频9| 亚洲欧美成人综合另类久久久| 亚洲欧美日韩无卡精品| av不卡在线播放| 超碰av人人做人人爽久久| 亚洲欧美清纯卡通| 国产精品久久久久久久久免| a级一级毛片免费在线观看| 精品国产一区二区三区久久久樱花 | 国产一区二区三区综合在线观看 | 男女下面进入的视频免费午夜| 中国三级夫妇交换| 激情 狠狠 欧美| 性色avwww在线观看| 国产亚洲一区二区精品| 欧美97在线视频| 成人亚洲精品一区在线观看 | 色婷婷av一区二区三区视频| 看十八女毛片水多多多| 日韩av不卡免费在线播放| 一级毛片 在线播放| 91在线精品国自产拍蜜月| 水蜜桃什么品种好| 国产深夜福利视频在线观看| 日日摸夜夜添夜夜添av毛片| 成年美女黄网站色视频大全免费 | 嘟嘟电影网在线观看| 成人一区二区视频在线观看| 久久久久久久久久人人人人人人| 欧美精品一区二区免费开放| 国产国拍精品亚洲av在线观看| 午夜老司机福利剧场| 香蕉精品网在线| 国产伦理片在线播放av一区| 人妻夜夜爽99麻豆av| a级毛片免费高清观看在线播放| 日本黄色片子视频| 亚洲电影在线观看av| 一级片'在线观看视频| 久久久久性生活片| av福利片在线观看| 一级片'在线观看视频| 久久精品国产亚洲网站| 日韩,欧美,国产一区二区三区| 亚洲av成人精品一二三区| 国语对白做爰xxxⅹ性视频网站| av在线老鸭窝| 亚洲真实伦在线观看| 欧美zozozo另类| 国产亚洲av片在线观看秒播厂| 日韩欧美 国产精品| 舔av片在线| 国产爽快片一区二区三区| 肉色欧美久久久久久久蜜桃| 亚洲精品久久午夜乱码| 欧美高清性xxxxhd video| 亚洲国产高清在线一区二区三| a级毛片免费高清观看在线播放| a 毛片基地| 99热这里只有是精品在线观看| 深爱激情五月婷婷| 97在线视频观看| 国产老妇伦熟女老妇高清| 国产淫语在线视频| 日韩国内少妇激情av| 人体艺术视频欧美日本| 亚洲一区二区三区欧美精品| 自拍欧美九色日韩亚洲蝌蚪91 | 美女福利国产在线 | 国产欧美亚洲国产| 亚洲真实伦在线观看| 欧美少妇被猛烈插入视频| 成年人午夜在线观看视频| 久久久a久久爽久久v久久| 国产av精品麻豆| 国产伦精品一区二区三区视频9| 男女国产视频网站| 一级毛片我不卡| 欧美3d第一页| 成年女人在线观看亚洲视频| 精品熟女少妇av免费看| 菩萨蛮人人尽说江南好唐韦庄| 亚洲,一卡二卡三卡| 亚洲精品456在线播放app| 国产在线男女| 日本-黄色视频高清免费观看| 亚洲国产精品专区欧美| 女人十人毛片免费观看3o分钟| 亚洲av综合色区一区| 91精品一卡2卡3卡4卡| 男的添女的下面高潮视频| 又粗又硬又长又爽又黄的视频| 欧美精品亚洲一区二区| 波野结衣二区三区在线| 视频中文字幕在线观看| a级一级毛片免费在线观看| 少妇丰满av| 五月天丁香电影| 久久久久久久久久成人| 夫妻午夜视频| 老熟女久久久| 精品人妻视频免费看| 亚洲综合精品二区| 亚洲国产日韩一区二区| 欧美另类一区| av在线老鸭窝| 日韩伦理黄色片| 国产精品久久久久久精品电影小说 | 日本vs欧美在线观看视频 | 国产av一区二区精品久久 | 妹子高潮喷水视频| 高清黄色对白视频在线免费看 | 晚上一个人看的免费电影| 在线看a的网站| 最近中文字幕高清免费大全6| 男人添女人高潮全过程视频| 亚洲av男天堂| 如何舔出高潮| 在线观看免费日韩欧美大片 | 久久久久久久精品精品| 亚洲人成网站在线观看播放| 中国三级夫妇交换| 纯流量卡能插随身wifi吗| 1000部很黄的大片| 一级毛片我不卡| av视频免费观看在线观看| 日产精品乱码卡一卡2卡三| 久久av网站| 一个人看的www免费观看视频| 亚洲美女视频黄频| 日本欧美国产在线视频| 亚洲精华国产精华液的使用体验| 亚洲天堂av无毛| 女人十人毛片免费观看3o分钟| 精品亚洲成国产av| 国产成人精品久久久久久| 欧美国产精品一级二级三级 | 中国国产av一级| 99热6这里只有精品| 免费人成在线观看视频色| 永久免费av网站大全| 2018国产大陆天天弄谢| 在线 av 中文字幕| 亚洲国产毛片av蜜桃av| 成人亚洲欧美一区二区av| 亚洲伊人久久精品综合| 日本一二三区视频观看| 新久久久久国产一级毛片| 肉色欧美久久久久久久蜜桃| 哪个播放器可以免费观看大片| 日本爱情动作片www.在线观看| 男人狂女人下面高潮的视频| 亚洲精品一二三| 两个人的视频大全免费| 黄片无遮挡物在线观看| 精品久久国产蜜桃| 在现免费观看毛片| 亚洲人成网站高清观看| 亚洲中文av在线| 国产一区二区在线观看日韩| 91精品伊人久久大香线蕉| av线在线观看网站| 天堂中文最新版在线下载| 国国产精品蜜臀av免费| 建设人人有责人人尽责人人享有的 | 久久精品熟女亚洲av麻豆精品| 国产黄色免费在线视频| av国产精品久久久久影院| av线在线观看网站| 高清黄色对白视频在线免费看 | 免费观看性生交大片5| 欧美成人午夜免费资源| 亚洲国产高清在线一区二区三| 丰满乱子伦码专区| 国产无遮挡羞羞视频在线观看| 色综合色国产| 久久人人爽人人爽人人片va| 久久久久性生活片| 精品酒店卫生间| 最近最新中文字幕免费大全7| 欧美精品亚洲一区二区| 人妻制服诱惑在线中文字幕| 十分钟在线观看高清视频www | 亚洲aⅴ乱码一区二区在线播放| 中文字幕人妻熟人妻熟丝袜美| 亚洲精品aⅴ在线观看| 天堂8中文在线网| 国产日韩欧美亚洲二区| 欧美精品一区二区免费开放| 欧美日韩国产mv在线观看视频 | av视频免费观看在线观看| 99re6热这里在线精品视频| 国产黄频视频在线观看| 王馨瑶露胸无遮挡在线观看| 久久99热6这里只有精品| 91狼人影院| 91久久精品电影网| 18禁裸乳无遮挡动漫免费视频| 99热6这里只有精品| 伦理电影大哥的女人| 日韩大片免费观看网站| 五月天丁香电影| 成年女人在线观看亚洲视频| 亚洲精品视频女| 亚洲色图综合在线观看| 精品亚洲成国产av| 国产爱豆传媒在线观看| 国产综合精华液| 亚洲精品一区蜜桃| 国内少妇人妻偷人精品xxx网站| 国产欧美另类精品又又久久亚洲欧美| 大香蕉97超碰在线| 国产成人a∨麻豆精品| 中文在线观看免费www的网站| 高清黄色对白视频在线免费看 | 最近中文字幕2019免费版| 日本av免费视频播放| 美女xxoo啪啪120秒动态图| 欧美日韩国产mv在线观看视频 | 亚州av有码| 亚洲欧美中文字幕日韩二区| 亚洲最大成人中文| 久久99热这里只有精品18| 性色avwww在线观看| 亚洲图色成人| 另类亚洲欧美激情| 亚洲精品日韩av片在线观看| 国模一区二区三区四区视频| 久久这里有精品视频免费| av免费观看日本| 亚洲av免费高清在线观看| 国产欧美亚洲国产| 在线观看三级黄色| 热re99久久精品国产66热6| 99视频精品全部免费 在线| 伦理电影免费视频| 欧美成人午夜免费资源| 国产爱豆传媒在线观看| 久久97久久精品| 美女中出高潮动态图| 日本av免费视频播放| 日韩亚洲欧美综合| 99国产精品免费福利视频| 99久久中文字幕三级久久日本| 欧美极品一区二区三区四区| 搡女人真爽免费视频火全软件| 女的被弄到高潮叫床怎么办| 亚洲成人中文字幕在线播放| 国产亚洲一区二区精品| 久久国产乱子免费精品| 午夜精品国产一区二区电影| 九九久久精品国产亚洲av麻豆| 免费黄网站久久成人精品| 久久久久精品久久久久真实原创| 高清毛片免费看| 搡女人真爽免费视频火全软件| 免费黄网站久久成人精品| 久久精品久久久久久噜噜老黄| 18禁动态无遮挡网站| 五月伊人婷婷丁香| 特大巨黑吊av在线直播| 亚洲美女视频黄频| 乱码一卡2卡4卡精品| 国产成人91sexporn| 99久久人妻综合| 日本vs欧美在线观看视频 | 女人久久www免费人成看片| 免费观看无遮挡的男女| 久久国产精品大桥未久av | 涩涩av久久男人的天堂| 日韩中文字幕视频在线看片 | 99九九线精品视频在线观看视频| 伦精品一区二区三区| 天美传媒精品一区二区| 看十八女毛片水多多多| 国语对白做爰xxxⅹ性视频网站| 欧美区成人在线视频| 一本久久精品| 男女免费视频国产| 日韩 亚洲 欧美在线| 日韩av在线免费看完整版不卡| 国产精品99久久99久久久不卡 | 99精国产麻豆久久婷婷| 免费观看的影片在线观看| 国产免费福利视频在线观看| 欧美精品亚洲一区二区| 在线观看免费日韩欧美大片 | 下体分泌物呈黄色| 久久久午夜欧美精品| 色视频www国产| 欧美日韩在线观看h| 美女视频免费永久观看网站| 亚洲美女视频黄频| 亚洲国产毛片av蜜桃av| 欧美国产精品一级二级三级 | 亚洲精品久久久久久婷婷小说| 少妇熟女欧美另类| av免费观看日本| 99热这里只有精品一区| 18禁裸乳无遮挡动漫免费视频| 欧美97在线视频| 黑人猛操日本美女一级片| 欧美老熟妇乱子伦牲交| 日日摸夜夜添夜夜爱| 国产高清不卡午夜福利| 国产老妇伦熟女老妇高清| 日韩av不卡免费在线播放| 建设人人有责人人尽责人人享有的 | 免费久久久久久久精品成人欧美视频 | av.在线天堂| 精品一区二区三区视频在线| 蜜臀久久99精品久久宅男| 男人狂女人下面高潮的视频| 日韩av不卡免费在线播放| 久久综合国产亚洲精品| 乱码一卡2卡4卡精品| 身体一侧抽搐| 色视频在线一区二区三区| 欧美97在线视频| 亚洲av.av天堂| 极品教师在线视频| 一区在线观看完整版| 亚洲不卡免费看| 色婷婷av一区二区三区视频| 国产淫语在线视频| 免费大片18禁| 伦理电影大哥的女人| 国产在视频线精品| 国产爱豆传媒在线观看| 亚洲精品中文字幕在线视频 | 亚洲av.av天堂| 青春草视频在线免费观看| 亚洲国产欧美在线一区| 国产乱人偷精品视频| 国产亚洲av片在线观看秒播厂| 精品国产一区二区三区久久久樱花 | 一区二区av电影网| 国产精品国产三级国产专区5o| 天堂中文最新版在线下载| 卡戴珊不雅视频在线播放| 久久毛片免费看一区二区三区| 国产精品国产三级国产av玫瑰| 亚洲av在线观看美女高潮| 久久99蜜桃精品久久| 一本久久精品| 夫妻性生交免费视频一级片| 日本av手机在线免费观看| 尾随美女入室| 成人影院久久| 亚洲丝袜综合中文字幕| 亚洲av男天堂| 在线亚洲精品国产二区图片欧美 | 亚洲av成人精品一区久久| av在线观看视频网站免费| 亚洲精品国产av成人精品| 精品少妇黑人巨大在线播放| 久久影院123| 国产高清三级在线| 国产精品成人在线| 一区二区三区乱码不卡18| 日韩国内少妇激情av| 噜噜噜噜噜久久久久久91| 国产午夜精品一二区理论片| 中文精品一卡2卡3卡4更新| 九九爱精品视频在线观看| 一区在线观看完整版| 美女视频免费永久观看网站| 免费观看a级毛片全部| 嘟嘟电影网在线观看| 一区在线观看完整版| 久热久热在线精品观看| 五月伊人婷婷丁香| 特大巨黑吊av在线直播| 久久精品国产自在天天线| 中文字幕亚洲精品专区| 精品国产乱码久久久久久小说| 午夜精品国产一区二区电影| 日韩精品有码人妻一区| 熟妇人妻不卡中文字幕| 久久久久久久久大av| 又黄又爽又刺激的免费视频.| 一区二区三区精品91| 一本久久精品| 国产av码专区亚洲av| 国产精品久久久久久av不卡| 黄色日韩在线| 好男人视频免费观看在线| 成人亚洲欧美一区二区av| 国语对白做爰xxxⅹ性视频网站| 看免费成人av毛片| 黄色怎么调成土黄色| 成人高潮视频无遮挡免费网站| 国国产精品蜜臀av免费| 亚洲av中文字字幕乱码综合| 国产亚洲91精品色在线| 亚洲美女搞黄在线观看| 视频区图区小说| 中文字幕制服av| 日本午夜av视频| 老熟女久久久| av免费在线看不卡| 国产成人a区在线观看| 一边亲一边摸免费视频| 亚洲第一av免费看| 国产一区亚洲一区在线观看| 成年女人在线观看亚洲视频| 亚洲精品色激情综合| 免费人成在线观看视频色| 啦啦啦在线观看免费高清www| 欧美一级a爱片免费观看看| 一级a做视频免费观看| 国产午夜精品一二区理论片| 久久久久人妻精品一区果冻| kizo精华| 青青草视频在线视频观看| 小蜜桃在线观看免费完整版高清| 精品午夜福利在线看| 男人添女人高潮全过程视频| 菩萨蛮人人尽说江南好唐韦庄| av在线app专区| 国产精品av视频在线免费观看| 国产在线一区二区三区精| 各种免费的搞黄视频| 美女中出高潮动态图| 男女边吃奶边做爰视频| 国产在线男女| 青青草视频在线视频观看| 看十八女毛片水多多多| 一区二区三区乱码不卡18| 最近中文字幕2019免费版| 国产亚洲精品久久久com| 纵有疾风起免费观看全集完整版| 色婷婷久久久亚洲欧美| 亚洲第一区二区三区不卡| 一本一本综合久久| 亚洲精品aⅴ在线观看| 日韩 亚洲 欧美在线| 午夜福利高清视频| 建设人人有责人人尽责人人享有的 | 日本色播在线视频| 99久久人妻综合| 国产精品国产三级国产av玫瑰| 男人添女人高潮全过程视频| 777米奇影视久久| 国产成人午夜福利电影在线观看| 啦啦啦在线观看免费高清www| 99热这里只有是精品50| 日韩大片免费观看网站| 一个人免费看片子| 七月丁香在线播放| 欧美97在线视频| 伦理电影免费视频| 九草在线视频观看| 亚洲,一卡二卡三卡| 女性生殖器流出的白浆| 免费观看的影片在线观看| 午夜福利在线观看免费完整高清在| 五月天丁香电影| 亚洲精品日韩在线中文字幕| 欧美极品一区二区三区四区| 看非洲黑人一级黄片| 狠狠精品人妻久久久久久综合| 国产一区二区在线观看日韩| 伊人久久国产一区二区| 草草在线视频免费看| 久久精品熟女亚洲av麻豆精品| av卡一久久| 自拍偷自拍亚洲精品老妇| 精品久久久噜噜| 91久久精品国产一区二区成人| 日日啪夜夜撸| 国产成人精品一,二区| 国产成人免费无遮挡视频| 久久久午夜欧美精品| 国产成人a区在线观看| 国产成人免费无遮挡视频| 这个男人来自地球电影免费观看 | 成人二区视频| 色哟哟·www| 国产伦在线观看视频一区| 91久久精品国产一区二区成人| 国产有黄有色有爽视频| 日产精品乱码卡一卡2卡三| 国产精品久久久久久精品电影小说 | 深夜a级毛片| 中文字幕人妻熟人妻熟丝袜美| 欧美另类一区| 国产精品国产av在线观看| av卡一久久|