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

    基于數(shù)字地形平滑和人工勢(shì)場(chǎng)的無(wú)人機(jī)近地三維路徑規(guī)劃

    2024-02-02 14:00:04郝向宇史涵意馬文峰
    兵器裝備工程學(xué)報(bào) 2024年1期
    關(guān)鍵詞:曲率坡度威脅

    朱 熠,田 輝,郝向宇,史涵意,馬文峰,王 聰

    (陸軍工程大學(xué) 野戰(zhàn)工程學(xué)院, 南京 210007)

    0 引言

    通常來(lái)說(shuō),路徑規(guī)劃的目的是生成一條最快、最簡(jiǎn)單的路徑,擺脫復(fù)雜地形、障礙物和威脅源等所有不利因素[1-3]。然而,當(dāng)無(wú)人機(jī)在近地區(qū)域飛行時(shí),由于飛行任務(wù)可能對(duì)地形有所限制和要求,如地形跟隨(terrain following,TF)和地形回避(terrain avoidance,TA)。因此,威脅因素仍在考慮之中,此時(shí)三維路徑規(guī)劃的任務(wù)就變得艱巨而富有挑戰(zhàn)性。

    在處理近地區(qū)域三維路徑規(guī)劃問題時(shí),首先需要一個(gè)數(shù)字地形模型。以前,地形模型通常由路徑設(shè)計(jì)者自己創(chuàng)建,現(xiàn)在有許多真實(shí)地圖的數(shù)字模型可用。例如,航天飛機(jī)地貌雷達(dá)任務(wù)(shuttle radar topography mission,SRTM)模型是最常用的模型之一[4-6],其所有數(shù)據(jù)覆蓋了地球80%以上的陸地表面[7-8]。通過對(duì)SRTM地形數(shù)據(jù)進(jìn)行一些輕微的轉(zhuǎn)換,可以使用Matlab軟件進(jìn)行編輯和處理[9-10]。通過地形平滑方法,可以生成安全飛行曲面。不少文獻(xiàn)通過對(duì)飛行軌跡的優(yōu)化,以滿足飛機(jī)的機(jī)動(dòng)性和地形跟隨的要求[11-12]。

    在通過地形平滑生成安全飛行曲面后,可以通過搜索算法獲得最優(yōu)飛行路徑。關(guān)于無(wú)人機(jī)路徑規(guī)劃算法大多是基于搜索算法的改進(jìn)和評(píng)估指標(biāo)的優(yōu)化[13-14]。本文中提出一種基于數(shù)字地形平滑和人工勢(shì)場(chǎng)法的三維路徑規(guī)劃方法,該方法將地形平滑技術(shù)和向量場(chǎng)搜索相結(jié)合。仿真結(jié)果表明,所設(shè)計(jì)的方法可以在安全飛行曲面上獲得最優(yōu)的三維路徑。

    1 數(shù)字地形平滑

    1.1 數(shù)字地形模型

    為獲得數(shù)字地形模型,采用SRTM數(shù)字高程模型。下載SRTM的59_06地形數(shù)據(jù),如圖1(a)所示。由于原始59_06地形非常大,因此只選擇了其中的一小部分,如圖1(b)所示。

    圖1 59_06地形數(shù)據(jù)和所選定區(qū)域

    使用Global Mapper或Surfer軟件可以方便地編輯SRTM地形數(shù)據(jù)。通過Surfer,SRTM地形數(shù)據(jù)可以從.grd格式轉(zhuǎn)換為.dat格式,.dat格式的文件可以在Matlab軟件中進(jìn)行編輯。

    1.2 地形平滑處理

    無(wú)人機(jī)的飛行狀態(tài)會(huì)受到其機(jī)動(dòng)性能的影響,這些性能主要包括最大爬升角、最小轉(zhuǎn)彎半徑、最大機(jī)動(dòng)過載、飛行速度和飛行高度的限制。通過平滑地形模型并生成安全飛行面可以滿足無(wú)人機(jī)機(jī)動(dòng)性能的這些要求。在本節(jié)中,提出了設(shè)置地面以上安全高度、減小坡度和減小曲率的設(shè)計(jì)思路。

    1) 設(shè)置高于地面的安全高度。地面的隱蔽要求和防撞要求之間存在矛盾。為了實(shí)現(xiàn)平衡,設(shè)置了地面以上的最低安全高度。當(dāng)無(wú)人機(jī)飛得高于安全高度時(shí),飛機(jī)可以避免撞到地面。在飛行高度高于安全高度的情況下,無(wú)人機(jī)可以最大限度地降低航跡的高度。在不同的地形中,最低安全海拔是不同的:在海洋和平原地區(qū),安全海拔較低;而在丘陵和山區(qū),安全海拔較高。一般情況下,適宜的安全海拔為:15~30 m的海域、60~80 m的平原區(qū)和約120 m的丘陵區(qū)。

    2) 降低坡度。在低空突防任務(wù)中,無(wú)人機(jī)有一個(gè)最大爬升角的限制。如果前方有坡度大于最大爬升角的山峰或障礙物,無(wú)人機(jī)需要提前爬升。否則,將導(dǎo)致與山體或障礙物發(fā)生碰撞。因此,需要對(duì)坡度較大的地形進(jìn)行修改,以滿足無(wú)人機(jī)縱向飛行的約束。坡度約束由無(wú)人機(jī)的最大爬升角決定。最常用的方法是用無(wú)人機(jī)的最大爬升角代替原始地形的最大坡度。坡度降低方法如圖2所示。

    圖2 減小坡度

    圖2中,i為描述高度信息的離散網(wǎng)格點(diǎn),h(i)為對(duì)應(yīng)的高程值,d為網(wǎng)格間距。記坡度為k=(h(i+1)-h(i))/d,最大可飛坡度為kmax,坡度處理邏輯為:

    ifk>kmax,h(i)=h(i+1)-d·kmax,

    i=1,…,n

    else ifk<-kmax,h(i+1)=h(i)-l·kmax

    elseh(i)=h(i),h(i+1)=h(i+1)

    上述邏輯只考慮了相鄰2個(gè)點(diǎn)之間的坡度問題,沒有考慮坡度的整體變化趨勢(shì),需要多次迭代才能得到理想的效果,甚至可能會(huì)出現(xiàn)峰巔相互拉高而不收斂的現(xiàn)象。

    綜合來(lái)說(shuō)主要的改進(jìn)思路包括2個(gè)方面:

    ① 以最高山峰點(diǎn)為起點(diǎn),依次對(duì)其兩側(cè)地形坡度進(jìn)行處理。首先尋找每行每列最高的山峰,標(biāo)記它們的位置;然后從最高山峰開始分別從正向和負(fù)向2個(gè)方向進(jìn)行坡度約束處理。

    ② 對(duì)山峰進(jìn)行保護(hù),對(duì)小于最大可飛坡度的山峰盡量少做處理。

    采用上述改進(jìn)邏輯,可以提高程序處理的速度,山峰不會(huì)互相削平或者相互拔高。

    3) 降低曲率。無(wú)人機(jī)在鉛垂平面內(nèi)運(yùn)動(dòng)軌跡的曲率ρa(bǔ)和法向加速度nz的關(guān)系為

    (1)

    根據(jù)無(wú)人機(jī)的速度和最大法向加速度計(jì)算出地形曲率的最大值ρa(bǔ)max,然后對(duì)地形的曲率ρx、ρy進(jìn)行約束處理,使它們的絕對(duì)值小于ρa(bǔ)max。由于相互正交的x和y方向的處理方法相同,下面只介紹x方向的處理方法。

    設(shè)地形高程函數(shù)為h=h(x,y),則曲率ρx為

    (2)

    圖3 地形曲率約束處理

    為求柵格點(diǎn)的曲率必須已知其一階偏導(dǎo)數(shù)和二階偏導(dǎo)數(shù)。利用離散數(shù)據(jù)求近似偏導(dǎo)數(shù)的方法較多??紤]到在實(shí)時(shí)軌跡優(yōu)化時(shí),采用3次卷積插值算法對(duì)數(shù)字地形進(jìn)行內(nèi)插。用上述方法內(nèi)插出的地形來(lái)求柵格點(diǎn)的一階偏導(dǎo)數(shù)和二階偏導(dǎo)數(shù),然后求出柵格點(diǎn)的曲率。曲率約束處理采用與坡度約束處理同樣的邏輯,具體步驟如下:

    步驟1在坡度約束處理完成的基礎(chǔ)上,從最高山峰點(diǎn)兩側(cè)分別開始進(jìn)行曲率約束處理;

    步驟2網(wǎng)格點(diǎn)(xi,yj)的一階偏導(dǎo)數(shù)為

    h′x=(hi+1, j-hi-1, j)/2l

    (3)

    由于3次卷積插值在柵格點(diǎn)處的二階偏導(dǎo)數(shù)不連續(xù),因此網(wǎng)格點(diǎn)(xi,yj)處的二階偏導(dǎo)數(shù)需要分左右2種情況來(lái)求。網(wǎng)格點(diǎn)(xi,yj)處的二階右偏導(dǎo)數(shù)為

    h″xr=(2hi-1, j-5hi, j+4hi+1, j-hi+2, j)/l2

    (4)

    網(wǎng)格點(diǎn)(xi,yj)處的二階左偏導(dǎo)數(shù)須利用網(wǎng)格點(diǎn)(xi-1,yj)處展開的插值公式求得,表達(dá)式為

    h″xl=(-hi-2, j+4hi-1, j-5hi, j+2hi+1, j)/l2

    (5)

    把式(3)代入式(5)中可得網(wǎng)格點(diǎn)處的左右曲率。

    步驟3網(wǎng)格曲率約束的處理會(huì)改變?cè)瓉?lái)已進(jìn)行過坡度修正的地形,因而每次計(jì)算網(wǎng)格曲率之前都要進(jìn)行網(wǎng)格坡度的重新判斷和調(diào)整;

    步驟4措施3已確保了地形坡度約束是滿足的,因而山峰和山谷之間的網(wǎng)格點(diǎn)是不需要進(jìn)行曲率修正的,曲率的約束處理只存在于山峰和山谷。為了對(duì)山峰進(jìn)行保護(hù),山峰和山谷須分別設(shè)計(jì)處理方法。

    1.3 安全飛行曲面生成

    在本節(jié)中,在Matlab中分析了SRTM高程地形模型。Matlab中的“l(fā)oad”函數(shù)可用于讀取SRTM的.dat文件,該文件已轉(zhuǎn)換為包含經(jīng)度、緯度和高度分量的三維矢量?!熬W(wǎng)格”函數(shù)用于展開平面中的點(diǎn)。插值函數(shù)“gridata”用于完成高度方向上的點(diǎn)。然后可以在Matlab中使用“網(wǎng)格”函數(shù)生成數(shù)字地形?!癱ontour3”函數(shù)可用于生成等高線,如圖4(a)和圖4(b)所示。

    圖4 地形數(shù)據(jù)的等高線

    然后,上述方法可以用于平滑數(shù)字地形,主要考慮的因素是海拔、坡度和曲率的安全系數(shù)。平滑的結(jié)果如圖5所示。

    通過比較圖5和圖4(a),可以看到地形平滑的結(jié)果,即生成的安全飛行曲面。可見,總體趨勢(shì)是完全保留的,最低高度提高了約150 m,最高山峰的高度保持不變,相鄰的陡峭山峰平滑地連接在一起,因此為了安全飛行,坡度和曲率都減小了。

    圖5 平滑處理后的安全飛行曲面

    2 三維最優(yōu)路徑規(guī)劃

    2.1 機(jī)動(dòng)性約束指標(biāo)

    基于戰(zhàn)術(shù)運(yùn)用和飛行器性能的要求,無(wú)人機(jī)的機(jī)動(dòng)性約束一般情況下包括以下5點(diǎn):

    1) 最短直飛距離:即飛行器在姿態(tài)改變前必須保持的最短直飛距離,設(shè)航跡由{Li|i=1,2,…,n}組成,最短直飛距離為L(zhǎng)min,該約束可表示為

    Li>Lmin(i=1,…,n)

    (6)

    2) 最小轉(zhuǎn)彎半徑:描述的是無(wú)人機(jī)過載約束的能力,在轉(zhuǎn)彎過程中不考慮切向加速度因素,則有:

    (7)

    式(7)中:R為轉(zhuǎn)彎半徑,V為無(wú)人機(jī)飛行速度,φ為滾轉(zhuǎn)角。

    由于無(wú)人機(jī)有最小飛行速度Vmin,而且有最大滾轉(zhuǎn)角φmax限制,所以R有最小限制,當(dāng)規(guī)劃的航跡的曲率大于該曲率時(shí),無(wú)人機(jī)受最大過載限制將不能準(zhǔn)確跟蹤航跡飛行。因此,路徑的曲率半徑必須大于無(wú)人機(jī)最小轉(zhuǎn)彎半徑。即約束表達(dá)式如下:

    (8)

    3) 最大爬升/俯沖角:最大爬升/俯沖角的限制是產(chǎn)生的航跡在飛行控制系統(tǒng)能夠允許的最大爬升/俯沖角內(nèi)。假設(shè)最大爬升/俯沖角為γmax,該約束可以這樣表示:

    (9)

    4) 最大飛行路徑:對(duì)于燃油和飛行時(shí)間的限制,飛行路徑的總長(zhǎng)度必須小于預(yù)設(shè)的最大距離。將飛行路徑的最大長(zhǎng)度設(shè)置為L(zhǎng)max,這樣該約束可以表示為

    (10)

    5) 目標(biāo)進(jìn)入方向:對(duì)于困難任務(wù)的要求,有時(shí)無(wú)人機(jī)需要從某個(gè)方向接近目標(biāo)。這就是進(jìn)入方向的約束。

    2.2 威脅評(píng)估指標(biāo)

    飛行路徑規(guī)劃的評(píng)估指標(biāo)設(shè)計(jì)為以下2種:

    1) 代價(jià)指標(biāo)

    (11)

    2) 威脅指標(biāo)

    (12)

    (13)

    式(13)中:Kj為第j個(gè)威脅源的強(qiáng)度,Rsj為飛行器距第j個(gè)威脅源的斜距,Nt為已知威脅的個(gè)數(shù)。

    2.3 最優(yōu)路徑搜索算法

    航跡生成的原理分成2步:首先在xy平面內(nèi)用路徑搜索算法規(guī)劃出能夠規(guī)避威脅的最優(yōu)二維路徑;然后根據(jù)第1節(jié)所得地形平滑后的安全曲面,擴(kuò)展出三維航路。故本節(jié)重點(diǎn)介紹xy平面內(nèi)的路徑搜索算法。

    1) 算法描述。算法采用人工勢(shì)場(chǎng)法的基本思想,其基本原理就是對(duì)所處障礙環(huán)境中每一點(diǎn)都賦一個(gè)勢(shì)場(chǎng)值,其值是目標(biāo)點(diǎn)的引力和障礙點(diǎn)的斥力疊加之和。因此路徑規(guī)劃也就是從起始點(diǎn)沿著勢(shì)場(chǎng)下降最快的方向到達(dá)目標(biāo)點(diǎn)。

    將威脅源看成“規(guī)避場(chǎng)”,即每一個(gè)威脅障礙的位置原點(diǎn)都看作在產(chǎn)生一個(gè)排斥性的斥力場(chǎng)作用在即將路過的無(wú)人機(jī)上面。產(chǎn)生的斥力大小與它們之間距離的平方成反比,平面內(nèi)任意一點(diǎn)受到第i個(gè)威脅障礙的斥力計(jì)算公式為

    (14)

    對(duì)于路徑的目標(biāo)點(diǎn)同樣規(guī)劃一個(gè)場(chǎng),這是一個(gè)吸引力場(chǎng),吸引力的大小正比于無(wú)人機(jī)和目標(biāo)的距離,其計(jì)算公式為

    (15)

    式(15)中:xaim、yaim為目標(biāo)點(diǎn)的位置坐標(biāo)。這個(gè)力同樣對(duì)平面中的點(diǎn)生成表格,是一張二維引力表。設(shè)當(dāng)無(wú)人機(jī)處于初始位置時(shí),引力達(dá)到最大值;當(dāng)無(wú)人機(jī)運(yùn)行到目標(biāo)點(diǎn)時(shí),引力達(dá)到最小值。這樣的一個(gè)10×10網(wǎng)格點(diǎn)的引力表格如表2所示,表中數(shù)值為900的是無(wú)人機(jī)初始位置即引力最大的點(diǎn),數(shù)值為300的是目標(biāo)位置即引力最小的點(diǎn)。

    2) 搜索過程。對(duì)斥力表和引力表進(jìn)行合成,可以生成一張合力表,通過合力引導(dǎo)無(wú)人機(jī)運(yùn)行出最優(yōu)軌跡。具體步驟如下:

    步驟1將整個(gè)平面區(qū)域離散化成100×100的網(wǎng)格地圖,然后從第一個(gè)網(wǎng)格點(diǎn)P1,1開始,選取一個(gè)3×3的子圖,共包含9個(gè)網(wǎng)格元素;

    步驟2計(jì)算這個(gè)3×3的子圖中的9個(gè)元素相互之間的受力值,選出所受合力最小的網(wǎng)格點(diǎn);

    步驟3進(jìn)行解的攝動(dòng)操作,即以步驟2中所得的合力最小網(wǎng)格點(diǎn)為中心,重新選取一個(gè)3×3的子圖;

    步驟4比較當(dāng)前最小網(wǎng)格點(diǎn)的坐標(biāo)與目標(biāo)點(diǎn)是否相等,如果相等則終止程序,否則繼續(xù)程序;

    步驟5重復(fù)步驟2—步驟4直到完成所有100×100網(wǎng)格點(diǎn)。

    3) 局部最優(yōu)的處理。如果在第2)步搜索過程中無(wú)人機(jī)進(jìn)入了一個(gè)區(qū)域,算出的當(dāng)前最小合力網(wǎng)格點(diǎn)與上一步得出的網(wǎng)格點(diǎn)相同,而且這個(gè)點(diǎn)不是目標(biāo)點(diǎn),則搜索算法進(jìn)入了一個(gè)局部最優(yōu)循環(huán)中。

    處理這個(gè)局部最優(yōu)循環(huán)的方法是在這個(gè)點(diǎn)的位置增加一個(gè)高能量的威脅場(chǎng),然后重新開始第2)步,這樣無(wú)人機(jī)就會(huì)避開這個(gè)區(qū)域。

    表1 10×10網(wǎng)格點(diǎn)的單個(gè)威脅場(chǎng)斥力表

    表2 10×10網(wǎng)格點(diǎn)的目標(biāo)點(diǎn)引力表

    3 三維路徑規(guī)劃仿真

    根據(jù)上節(jié)描述的路徑搜索算法,并參考文獻(xiàn)[15]中的仿真框架設(shè)計(jì)思路,采用Matlab軟件設(shè)計(jì)了GUI(graphical user interface,圖形用戶界面)程序,優(yōu)化了程序界面布局和回調(diào)函數(shù)中的搜索算法,改進(jìn)后的程序能夠直接導(dǎo)入1.3節(jié)生成的平滑后的安全飛行曲面。其中相關(guān)控件的名稱和功能如表3所示。

    表3 仿真界面中各控件功能說(shuō)明

    平滑后的數(shù)字地形及其等高線就顯示在坐標(biāo)軸系1和坐標(biāo)軸系2中。在表3中所有控件的回調(diào)函數(shù)都已編碼之后,GUI程序就設(shè)計(jì)完成了。

    路徑規(guī)劃結(jié)果分4個(gè)步驟如下:

    步驟1加載地圖?!拜d入地形”按鈕是指加載平滑處理后的三維地形的地圖,如圖6所示,即包括平滑處理后的安全飛行曲面及其相應(yīng)的等高線。在圖6中,3D地形的軸的單位由之前的緯度(°)、經(jīng)度(°)和米(m)轉(zhuǎn)換而來(lái)。

    圖6 在GUI界面中載入三維平滑地形

    步驟2添加新威脅。點(diǎn)擊界面中的“新增威脅障礙”按鈕,即可添加威脅源。在等高線的任意位置左鍵單擊1次,再單擊鼠標(biāo)右鍵即可完成一處威脅源的添加,柱狀威脅源將出現(xiàn)在三維地形和等高線中。在任意位置先后隨機(jī)設(shè)置5個(gè)威脅源后,得到如圖7中所示的5排柱形威脅源。

    圖7 設(shè)置威脅障礙

    步驟3確定目標(biāo)點(diǎn)。點(diǎn)擊界面中的“輸入目標(biāo)位置”按鈕。然后左鍵單擊等高線圖中威脅屏障外的任何位置,就會(huì)顯示一個(gè)紅色的小圓圈作為標(biāo)記。

    步驟4開始搜索。點(diǎn)擊界面中的“開始生成航路”按鈕。程序開始運(yùn)行,依電腦計(jì)算性能不同,等待約1~2 s后,最佳飛行路徑就出現(xiàn)在3D地形圖和等高線圖中,如圖8(a)所示。

    在圖8(a)中,在右側(cè)的等高線圖中,繪制了一條藍(lán)色星形線作為最佳路徑。同時(shí),在左側(cè)的3D地形圖中,繪制了一條紅線作為最佳路徑??梢钥闯?所設(shè)計(jì)的路徑規(guī)劃算法可以有效地避免威脅,快速找到到達(dá)目標(biāo)的最短路徑,并緊貼3D安全飛行曲面之上飛行。

    在圖8(b)中單獨(dú)繪制出了規(guī)劃好最優(yōu)軌跡(圖8中紅色曲線)和搜索過程中形成的最優(yōu)航點(diǎn)的連線(圖8中藍(lán)色曲線),顯然,最終的最優(yōu)軌跡是在最優(yōu)航點(diǎn)連線的基礎(chǔ)上平滑處理之后得到的。

    圖8 三維最優(yōu)路徑規(guī)劃結(jié)果

    4 結(jié)論

    本文中設(shè)計(jì)了一種無(wú)人機(jī)近地三維最優(yōu)路徑規(guī)劃方法,通過仿真實(shí)驗(yàn),可以得到以下結(jié)論:

    1) 通過真實(shí)地圖的數(shù)字高程模型作為無(wú)人機(jī)近地路徑規(guī)劃的地形數(shù)據(jù)來(lái)源是一種簡(jiǎn)潔、可行的地形建模方法。

    2) 通過考慮安全高度、減小坡度和曲率等因素對(duì)原始數(shù)字地形進(jìn)行預(yù)處理,可以有效的生成無(wú)人機(jī)的安全飛行曲面,形成對(duì)復(fù)雜三維地形的路徑預(yù)規(guī)劃。

    3) 所設(shè)計(jì)的GUI程序可以加載優(yōu)化后的安全飛行曲面,并通過最優(yōu)算法避開威脅區(qū)域到達(dá)目標(biāo)位置(其中威脅數(shù)量、威脅位置和目標(biāo)位置可自由指定),獲得最佳的飛行路徑。

    本文中設(shè)計(jì)的規(guī)劃方法,在安全飛行曲面的基礎(chǔ)上,可以避開特定的威脅,獲得最佳的飛行路徑,且滿足無(wú)人機(jī)操縱性能的約束要求,為無(wú)人機(jī)近地三維路徑規(guī)劃的設(shè)計(jì)與仿真提供了一種新的思路。

    猜你喜歡
    曲率坡度威脅
    大曲率沉管安裝關(guān)鍵技術(shù)研究
    一類雙曲平均曲率流的對(duì)稱與整體解
    人類的威脅
    半正迷向曲率的四維Shrinking Gradient Ricci Solitons
    關(guān)于公路超高漸變段合成坡度解析與應(yīng)用
    受到威脅的生命
    面對(duì)孩子的“威脅”,我們要會(huì)說(shuō)“不”
    家教世界(2017年11期)2018-01-03 01:28:49
    基于圖像處理的定位器坡度計(jì)算
    電氣化鐵道(2016年4期)2016-04-16 05:59:46
    Why Does Sleeping in Just Make Us More Tired?
    坡度在巖石風(fēng)化層解譯中的應(yīng)用
    河北遙感(2015年2期)2015-07-18 11:11:14
    亚洲中文字幕一区二区三区有码在线看| 国产伦在线观看视频一区| 久久99热6这里只有精品| 久久久久性生活片| 91麻豆av在线| xxxwww97欧美| aaaaa片日本免费| 又爽又黄无遮挡网站| 亚洲精品色激情综合| 欧美黑人欧美精品刺激| 久久精品国产自在天天线| 国产不卡一卡二| 精华霜和精华液先用哪个| 日本 av在线| 美女大奶头视频| 国产高潮美女av| 国产老妇女一区| 伊人久久精品亚洲午夜| 老熟妇乱子伦视频在线观看| 伊人久久精品亚洲午夜| 国产视频一区二区在线看| 啦啦啦观看免费观看视频高清| 免费一级毛片在线播放高清视频| 国产又黄又爽又无遮挡在线| 99热6这里只有精品| 最后的刺客免费高清国语| 少妇丰满av| 国产精品综合久久久久久久免费| 亚洲成人久久性| 麻豆成人av在线观看| 99久久九九国产精品国产免费| 日韩大尺度精品在线看网址| 老司机午夜福利在线观看视频| 日本撒尿小便嘘嘘汇集6| 国产精品久久久久久av不卡| 国产一区二区在线观看日韩| 国产精品亚洲一级av第二区| 欧美高清成人免费视频www| 黄色欧美视频在线观看| 成人亚洲精品av一区二区| videossex国产| 99riav亚洲国产免费| 成人特级av手机在线观看| 国产精品一区二区三区四区久久| 久久人人精品亚洲av| 网址你懂的国产日韩在线| 日韩在线高清观看一区二区三区 | 欧美极品一区二区三区四区| 亚洲五月天丁香| 变态另类丝袜制服| 男女下面进入的视频免费午夜| 国产女主播在线喷水免费视频网站 | 亚洲av熟女| 在线免费观看的www视频| 久久精品久久久久久噜噜老黄 | 少妇熟女aⅴ在线视频| 91狼人影院| 久久久久国内视频| 五月伊人婷婷丁香| 老熟妇乱子伦视频在线观看| av黄色大香蕉| 性色avwww在线观看| 99热只有精品国产| 91午夜精品亚洲一区二区三区 | 免费高清视频大片| 我的老师免费观看完整版| 国产激情偷乱视频一区二区| 国产高清三级在线| 久99久视频精品免费| 国产亚洲av嫩草精品影院| 观看美女的网站| 亚洲精品粉嫩美女一区| 色吧在线观看| 18禁黄网站禁片免费观看直播| 久久欧美精品欧美久久欧美| 在线观看舔阴道视频| 国产av麻豆久久久久久久| 成人二区视频| 18禁黄网站禁片午夜丰满| 1024手机看黄色片| 国产亚洲精品久久久com| 嫁个100分男人电影在线观看| 国产精品女同一区二区软件 | 在线看三级毛片| 男人和女人高潮做爰伦理| 制服丝袜大香蕉在线| 欧美黑人欧美精品刺激| 久久精品久久久久久噜噜老黄 | 午夜福利在线观看免费完整高清在 | 女同久久另类99精品国产91| 人人妻人人看人人澡| 成人亚洲精品av一区二区| 人人妻人人看人人澡| 中国美女看黄片| 九色成人免费人妻av| 亚洲精品色激情综合| 婷婷六月久久综合丁香| 亚洲精品在线观看二区| 国内精品美女久久久久久| 中文字幕久久专区| 免费看日本二区| 久久久久性生活片| 午夜福利在线观看吧| 一个人免费在线观看电影| 女人十人毛片免费观看3o分钟| 亚洲av一区综合| 亚洲熟妇中文字幕五十中出| 午夜影院日韩av| 亚洲av中文av极速乱 | 亚洲av熟女| 啪啪无遮挡十八禁网站| 成人无遮挡网站| 久久午夜福利片| 男女视频在线观看网站免费| 中亚洲国语对白在线视频| 18禁在线播放成人免费| 我的女老师完整版在线观看| 国产精品久久久久久亚洲av鲁大| 欧美激情久久久久久爽电影| 国产探花极品一区二区| 天天一区二区日本电影三级| 天天一区二区日本电影三级| av福利片在线观看| 国内精品一区二区在线观看| 成人国产综合亚洲| 网址你懂的国产日韩在线| 看黄色毛片网站| 久久国产精品人妻蜜桃| 99久久成人亚洲精品观看| 成人二区视频| 天堂网av新在线| 联通29元200g的流量卡| 日日夜夜操网爽| 午夜视频国产福利| 啦啦啦观看免费观看视频高清| 精品国内亚洲2022精品成人| 黄色欧美视频在线观看| a级一级毛片免费在线观看| 久久精品国产亚洲av涩爱 | 婷婷精品国产亚洲av| 欧美日韩亚洲国产一区二区在线观看| 嫩草影院入口| 久久久久久大精品| 国产精品av视频在线免费观看| 搡女人真爽免费视频火全软件 | 乱码一卡2卡4卡精品| 老女人水多毛片| 国产精品永久免费网站| av黄色大香蕉| 人妻丰满熟妇av一区二区三区| 久久国产精品人妻蜜桃| 国产激情偷乱视频一区二区| 中文资源天堂在线| 91狼人影院| 国产高清不卡午夜福利| 三级国产精品欧美在线观看| 国产乱人视频| 精品人妻视频免费看| 热99在线观看视频| 日本 av在线| 黄色女人牲交| 亚洲精品色激情综合| 国产一区二区激情短视频| 一个人免费在线观看电影| 久久亚洲真实| 国模一区二区三区四区视频| www.www免费av| 岛国在线免费视频观看| 亚洲一级一片aⅴ在线观看| 欧美又色又爽又黄视频| 在线a可以看的网站| 动漫黄色视频在线观看| 日本-黄色视频高清免费观看| 全区人妻精品视频| 变态另类成人亚洲欧美熟女| 一个人看的www免费观看视频| av福利片在线观看| 国产视频一区二区在线看| 亚洲国产高清在线一区二区三| 亚洲国产高清在线一区二区三| av.在线天堂| 欧美日韩黄片免| 国产高清视频在线播放一区| 国产高潮美女av| 亚洲专区中文字幕在线| 男人舔奶头视频| 97超级碰碰碰精品色视频在线观看| 99riav亚洲国产免费| 欧美成人一区二区免费高清观看| 99久久无色码亚洲精品果冻| 免费搜索国产男女视频| 99国产精品一区二区蜜桃av| 日日干狠狠操夜夜爽| 成人性生交大片免费视频hd| 在现免费观看毛片| 在线免费十八禁| 国产视频内射| 精品久久久噜噜| 赤兔流量卡办理| 搡老熟女国产l中国老女人| 韩国av一区二区三区四区| 九九在线视频观看精品| 在线a可以看的网站| 国产高清视频在线播放一区| 美女免费视频网站| 国产成人福利小说| 欧美成人性av电影在线观看| 男插女下体视频免费在线播放| 一边摸一边抽搐一进一小说| 简卡轻食公司| 国产在线精品亚洲第一网站| 小蜜桃在线观看免费完整版高清| 18+在线观看网站| h日本视频在线播放| av.在线天堂| 午夜福利在线在线| 黄色女人牲交| 露出奶头的视频| 国产精品久久视频播放| 国产麻豆成人av免费视频| 中文字幕精品亚洲无线码一区| 中出人妻视频一区二区| 欧美高清性xxxxhd video| 成人永久免费在线观看视频| 亚洲精品一区av在线观看| 国产主播在线观看一区二区| 亚洲乱码一区二区免费版| 伊人久久精品亚洲午夜| 欧美一区二区精品小视频在线| 赤兔流量卡办理| 搞女人的毛片| 国产毛片a区久久久久| 韩国av一区二区三区四区| 色综合色国产| 欧美日韩黄片免| 免费av毛片视频| 国产一区二区三区在线臀色熟女| 69av精品久久久久久| 国内精品宾馆在线| 精品午夜福利视频在线观看一区| 亚洲三级黄色毛片| 亚洲va日本ⅴa欧美va伊人久久| 亚洲国产高清在线一区二区三| 99国产极品粉嫩在线观看| av国产免费在线观看| 久久精品国产自在天天线| 国产黄a三级三级三级人| 91久久精品国产一区二区成人| 免费观看的影片在线观看| 日本撒尿小便嘘嘘汇集6| 一本精品99久久精品77| 少妇的逼水好多| 国产精品一区二区三区四区免费观看 | 亚洲内射少妇av| 亚洲一级一片aⅴ在线观看| 狠狠狠狠99中文字幕| 国产成人影院久久av| 51国产日韩欧美| 中文字幕人妻熟人妻熟丝袜美| 在线播放无遮挡| 1000部很黄的大片| 国产免费av片在线观看野外av| 欧美一区二区亚洲| 午夜福利视频1000在线观看| 女的被弄到高潮叫床怎么办 | 日本与韩国留学比较| 亚洲精品色激情综合| 国产一区二区三区在线臀色熟女| 欧美xxxx黑人xx丫x性爽| 亚洲四区av| 一级黄片播放器| 欧美色欧美亚洲另类二区| 人妻丰满熟妇av一区二区三区| 全区人妻精品视频| 成人三级黄色视频| 亚洲最大成人av| 久久九九热精品免费| 久久草成人影院| 国产亚洲精品综合一区在线观看| 国产午夜福利久久久久久| 久久亚洲精品不卡| 99在线人妻在线中文字幕| 午夜亚洲福利在线播放| 久久久久国内视频| 欧美日韩亚洲国产一区二区在线观看| 亚洲国产日韩欧美精品在线观看| 国产白丝娇喘喷水9色精品| 人人妻,人人澡人人爽秒播| 国产精品人妻久久久影院| 性插视频无遮挡在线免费观看| 国产69精品久久久久777片| 窝窝影院91人妻| 一个人看视频在线观看www免费| 亚洲av二区三区四区| 搞女人的毛片| 国产欧美日韩精品一区二区| 久久久久久久久久成人| 国产乱人视频| 国产免费av片在线观看野外av| 日本精品一区二区三区蜜桃| 国产黄色小视频在线观看| 男人和女人高潮做爰伦理| 亚洲男人的天堂狠狠| 久久久久国产精品人妻aⅴ院| 色在线成人网| 少妇裸体淫交视频免费看高清| 可以在线观看的亚洲视频| 乱码一卡2卡4卡精品| .国产精品久久| 草草在线视频免费看| 一个人观看的视频www高清免费观看| 国产精品久久久久久久久免| 熟妇人妻久久中文字幕3abv| 国产真实伦视频高清在线观看 | 日韩欧美国产在线观看| 永久网站在线| 国国产精品蜜臀av免费| 久99久视频精品免费| 级片在线观看| h日本视频在线播放| 亚洲国产精品sss在线观看| avwww免费| 国产精品三级大全| 日本 欧美在线| 日本欧美国产在线视频| 可以在线观看的亚洲视频| 亚洲va在线va天堂va国产| 国产精品嫩草影院av在线观看 | 国产精华一区二区三区| 美女cb高潮喷水在线观看| 久久精品国产99精品国产亚洲性色| 性欧美人与动物交配| 观看免费一级毛片| 三级国产精品欧美在线观看| 欧美又色又爽又黄视频| 99国产极品粉嫩在线观看| 欧美日韩精品成人综合77777| 成人午夜高清在线视频| 国产美女午夜福利| 欧美zozozo另类| 能在线免费观看的黄片| 亚洲成人久久爱视频| 欧美成人一区二区免费高清观看| 亚洲精品久久国产高清桃花| 男人的好看免费观看在线视频| 精品人妻熟女av久视频| 日本一本二区三区精品| 日韩在线高清观看一区二区三区 | 欧美潮喷喷水| 亚洲成av人片在线播放无| 大型黄色视频在线免费观看| 亚洲av美国av| 一个人免费在线观看电影| 亚洲精品久久国产高清桃花| 最好的美女福利视频网| 日本黄大片高清| 深夜a级毛片| 黄色女人牲交| 91久久精品国产一区二区三区| 精品一区二区三区视频在线| av黄色大香蕉| 国产亚洲精品综合一区在线观看| 亚洲成人久久性| 少妇丰满av| 国产一区二区亚洲精品在线观看| 欧美色视频一区免费| 99精品久久久久人妻精品| 如何舔出高潮| 精品一区二区三区av网在线观看| 在线a可以看的网站| 自拍偷自拍亚洲精品老妇| 久久中文看片网| 成年女人永久免费观看视频| 男人和女人高潮做爰伦理| 亚洲人成网站在线播| 麻豆av噜噜一区二区三区| 赤兔流量卡办理| 国产老妇女一区| 男女下面进入的视频免费午夜| 不卡视频在线观看欧美| 亚洲人成网站在线播| 国产精品人妻久久久久久| 大又大粗又爽又黄少妇毛片口| 看免费成人av毛片| 免费看a级黄色片| 成人国产麻豆网| 少妇的逼好多水| 色精品久久人妻99蜜桃| 麻豆成人午夜福利视频| 国产欧美日韩精品一区二区| 国产亚洲精品久久久久久毛片| 亚洲精品亚洲一区二区| 久久精品综合一区二区三区| 亚洲va在线va天堂va国产| 精品无人区乱码1区二区| 久久热精品热| 嫩草影视91久久| 国产精品三级大全| 国产精品久久久久久av不卡| 午夜激情欧美在线| 九色国产91popny在线| 熟妇人妻久久中文字幕3abv| 日韩欧美国产一区二区入口| 精品久久久久久久久亚洲 | 日韩欧美精品v在线| 香蕉av资源在线| 亚洲性夜色夜夜综合| 国产精华一区二区三区| 亚洲在线观看片| 亚洲精品一卡2卡三卡4卡5卡| 亚洲成人精品中文字幕电影| 深夜精品福利| 国产人妻一区二区三区在| 国产aⅴ精品一区二区三区波| 特级一级黄色大片| 美女黄网站色视频| 女人十人毛片免费观看3o分钟| 一级av片app| .国产精品久久| www.色视频.com| 亚洲四区av| 亚洲图色成人| 人妻制服诱惑在线中文字幕| 国产精品一区二区三区四区久久| 给我免费播放毛片高清在线观看| 欧美高清性xxxxhd video| 亚洲成人久久性| 国产私拍福利视频在线观看| 国产一区二区在线av高清观看| 日韩 亚洲 欧美在线| 日韩一区二区视频免费看| 午夜福利高清视频| 岛国在线免费视频观看| 成人av在线播放网站| 直男gayav资源| 身体一侧抽搐| 免费av毛片视频| а√天堂www在线а√下载| 九九爱精品视频在线观看| netflix在线观看网站| 日日啪夜夜撸| 久久精品国产自在天天线| 少妇人妻精品综合一区二区 | 亚洲美女视频黄频| 久9热在线精品视频| 亚洲三级黄色毛片| 波多野结衣高清无吗| 毛片一级片免费看久久久久 | 国产69精品久久久久777片| 我要搜黄色片| 露出奶头的视频| 免费黄网站久久成人精品| 婷婷精品国产亚洲av| 老师上课跳d突然被开到最大视频| 看片在线看免费视频| 婷婷亚洲欧美| 国产伦精品一区二区三区视频9| 日韩亚洲欧美综合| 精品久久久久久成人av| 啦啦啦韩国在线观看视频| 1000部很黄的大片| 久久久久久久午夜电影| 国产日本99.免费观看| 国产精品乱码一区二三区的特点| 亚洲一级一片aⅴ在线观看| 亚洲精华国产精华精| 亚洲久久久久久中文字幕| 悠悠久久av| 国产精品爽爽va在线观看网站| 久久久久精品国产欧美久久久| 国产精品国产高清国产av| 国产91精品成人一区二区三区| 免费在线观看成人毛片| 尾随美女入室| 亚洲精品色激情综合| 禁无遮挡网站| 欧美国产日韩亚洲一区| 国产精品一区二区性色av| 日韩欧美精品v在线| 69人妻影院| 国产精品一区二区免费欧美| 中文在线观看免费www的网站| 久久精品91蜜桃| 国产精品av视频在线免费观看| 丝袜美腿在线中文| 国产美女午夜福利| 国产精品人妻久久久久久| 婷婷精品国产亚洲av在线| 免费大片18禁| 久久精品国产亚洲网站| 国产伦一二天堂av在线观看| 男女做爰动态图高潮gif福利片| videossex国产| 最近最新中文字幕大全电影3| 自拍偷自拍亚洲精品老妇| 久久久久久久午夜电影| 成年免费大片在线观看| 日本熟妇午夜| 黄色女人牲交| 午夜精品久久久久久毛片777| 免费无遮挡裸体视频| 欧美黑人欧美精品刺激| 特大巨黑吊av在线直播| 韩国av一区二区三区四区| 人人妻,人人澡人人爽秒播| 乱系列少妇在线播放| 国产淫片久久久久久久久| 亚洲第一电影网av| 人人妻,人人澡人人爽秒播| 欧美一区二区亚洲| 中文字幕av在线有码专区| 久久亚洲精品不卡| or卡值多少钱| 在线播放国产精品三级| 日韩精品青青久久久久久| 久久精品国产亚洲av香蕉五月| 国产精品久久视频播放| 成人综合一区亚洲| 国产一区二区三区av在线 | 久久亚洲精品不卡| 长腿黑丝高跟| 免费黄网站久久成人精品| 亚洲中文字幕日韩| 91在线精品国自产拍蜜月| 91精品国产九色| 国产视频一区二区在线看| 嫩草影院新地址| 成人美女网站在线观看视频| 男女做爰动态图高潮gif福利片| 亚洲最大成人手机在线| 亚洲av成人精品一区久久| 久久人人精品亚洲av| 成人午夜高清在线视频| 欧美+日韩+精品| 99热网站在线观看| 国产成人影院久久av| 少妇高潮的动态图| 国产不卡一卡二| 日本黄色片子视频| 国产精品国产高清国产av| 欧美bdsm另类| 亚洲精品久久国产高清桃花| 国产aⅴ精品一区二区三区波| 亚洲中文字幕一区二区三区有码在线看| 欧美成人性av电影在线观看| 人妻丰满熟妇av一区二区三区| 最新中文字幕久久久久| 成年人黄色毛片网站| 国产精品电影一区二区三区| 久久久久久伊人网av| 日本一二三区视频观看| 精品福利观看| 欧美bdsm另类| 亚洲七黄色美女视频| 91在线观看av| 久久精品夜夜夜夜夜久久蜜豆| 熟女人妻精品中文字幕| 亚洲图色成人| 91午夜精品亚洲一区二区三区 | 18禁黄网站禁片午夜丰满| 欧美日韩国产亚洲二区| 亚洲欧美日韩无卡精品| 非洲黑人性xxxx精品又粗又长| 婷婷精品国产亚洲av| 日本爱情动作片www.在线观看 | 国产精品综合久久久久久久免费| 亚洲最大成人av| 少妇的逼水好多| 成人国产麻豆网| 如何舔出高潮| 欧美最新免费一区二区三区| 少妇的逼好多水| 国产大屁股一区二区在线视频| 亚洲成人久久性| 高清毛片免费观看视频网站| 真人做人爱边吃奶动态| 大型黄色视频在线免费观看| 俺也久久电影网| 午夜免费男女啪啪视频观看 | av黄色大香蕉| 88av欧美| 亚洲美女搞黄在线观看 | 伊人久久精品亚洲午夜| 国产av不卡久久| 国产亚洲精品综合一区在线观看| 久久精品影院6| 国产亚洲精品av在线| 国国产精品蜜臀av免费| av在线亚洲专区| 国产精品三级大全| 亚洲无线观看免费| 国产一区二区亚洲精品在线观看| 日韩在线高清观看一区二区三区 | 日日撸夜夜添| 久久久久久久亚洲中文字幕| 少妇人妻一区二区三区视频| 乱系列少妇在线播放| 人妻制服诱惑在线中文字幕| 男人舔奶头视频| 国模一区二区三区四区视频| 中文资源天堂在线| 少妇裸体淫交视频免费看高清| 可以在线观看毛片的网站| 女生性感内裤真人,穿戴方法视频| 乱码一卡2卡4卡精品| 亚洲av成人av| 中文字幕人妻熟人妻熟丝袜美| 久久人人爽人人爽人人片va| 亚洲天堂国产精品一区在线| 亚洲国产欧洲综合997久久,| 国产精品永久免费网站| 亚洲av免费在线观看| 国产美女午夜福利| av天堂在线播放| 久久亚洲真实|