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

    不確定環(huán)境下的多無人機協(xié)同搜索航路規(guī)劃

    2011-02-22 07:31:32吳文超黃長強宋磊唐上欽白壬潮
    兵工學(xué)報 2011年11期
    關(guān)鍵詞:禁飛區(qū)抗力航路

    吳文超,黃長強,宋磊,唐上欽,白壬潮

    (1.空軍工程大學(xué) 工程學(xué)院,陜西 西安710038;2.空軍裝備部,北京100843)

    利用無人機(UAV)執(zhí)行偵察任務(wù)在眾多領(lǐng)域都得到了廣泛應(yīng)用,如對敵區(qū)目標(biāo)情況進(jìn)行偵察監(jiān)視,在公海、沙漠或山區(qū)展開搜索營救以及對礦藏的資源勘察等。相對于有人駕駛飛機,UAV 成本低,可操縱性好,尤其適合于執(zhí)行危險的任務(wù)[1]。由于單架UAV 在信息收集方面的局限性,多架UAV 協(xié)同執(zhí)行復(fù)雜搜索任務(wù)吸引了越來越多的關(guān)注[2-3]。一些窮舉覆蓋航路規(guī)劃模型,如Zamboni 搜索方法已經(jīng)得到了發(fā)展應(yīng)用[4-6]。然而該方法缺乏靈活性,如在設(shè)計某地區(qū)的UAV 森林火災(zāi)偵察路線時,不能對該區(qū)域的湖泊以及其他無林地帶進(jìn)行規(guī)避,在載油有限和時間緊張的情況下,將影響到對其他更有價值區(qū)域的偵察活動。文獻(xiàn)[7]中提出的搜索算法可以對整個環(huán)境展開遍歷搜索,但沒有考慮UAV 的機動限制和搜索環(huán)境的動態(tài)特性。如UAV的物理機動性能會限制其最小轉(zhuǎn)彎半徑,而在大部分窮舉搜索算法中,允許UAV 向任何方向運動[8]。此外,任務(wù)時間和機載燃油的限制可能不允許UAV完全覆蓋整個環(huán)境,從而需要對目標(biāo)存在概率較高的特定區(qū)域進(jìn)行重點偵察。文獻(xiàn)[9]對UAV 的協(xié)同搜索展開了研究,UAV 在搜索過程中一定程度上可以避開已知區(qū)域,但是其算法無法保證對禁飛區(qū)實現(xiàn)完全回避。文獻(xiàn)[10]基于遺傳算法設(shè)計了UAV 的偵察航路,其研究是在假設(shè)目標(biāo)位置已知情況下進(jìn)行的,在環(huán)境中的目標(biāo)分布信息很少的情況下,應(yīng)用受到了限制。本文主要針對不確定環(huán)境下的多UAV 協(xié)同搜索問題展開研究,所謂不確定環(huán)境指的是在該環(huán)境中關(guān)于目標(biāo)的分布是有限的(甚至沒有)先驗信息。

    1 協(xié)同搜索的問題描述

    一種典型的UAV 協(xié)同搜索在線規(guī)劃和飛行控制框架如圖1所示。它包括上層決策(任務(wù)規(guī)劃)和由傳感器、跟蹤控制器和飛行器本體組成的底層控制2 部分。UAV 接收地面站發(fā)送的任務(wù),同時借助本機傳感器或它機傳感器(通過數(shù)據(jù)鏈)實時更新環(huán)境信息,進(jìn)行航路規(guī)劃產(chǎn)生參考飛行航跡,而后跟蹤控制器根據(jù)UAV 的動力學(xué)模型產(chǎn)生激勵輸入信號,控制UAV 沿參考軌跡飛行。

    與相互獨立控制的單架UAV 相比,多架UAV協(xié)同搜索能夠提供更有效的偵察能力。多架UAV在不確定環(huán)境中進(jìn)行協(xié)同搜索時,一般需要滿足以下原則[11-12]:

    1)在條件允許的情況下覆蓋盡可能多的區(qū)域。

    2)為了提高搜索效率,應(yīng)對環(huán)境中的無回報區(qū)域(環(huán)境信息完全已知區(qū)域)盡可能進(jìn)行規(guī)避,而對環(huán)境中目標(biāo)存在概率較大的區(qū)域進(jìn)行重點偵察。

    3)滿足UAV 機動限制要求,在數(shù)據(jù)通訊延遲和編隊中的某架UAV 出現(xiàn)故障時不影響其他UAV的協(xié)同搜索任務(wù)。

    4)保證UAV 飛行安全。

    需要說明:當(dāng)對環(huán)境信息完全不確定的情況下應(yīng)當(dāng)盡可能的進(jìn)行區(qū)域覆蓋,以求獲得盡可能多的目標(biāo)信息,而當(dāng)存在部分先驗已知信息時,需要對目標(biāo)存在概率很小的區(qū)域進(jìn)行規(guī)避,對目標(biāo)存在概率較大的區(qū)域進(jìn)行重點偵察。所以第1)、2)條只是應(yīng)用的前提條件不同,并不矛盾。第3)條是從應(yīng)用角度提出的,要求得到的路徑具有物理可行性,并能夠適應(yīng)突發(fā)事件的發(fā)生。第4)條是從安全角度提出的,由于政治、軍事或地理方面的原因,UAV 在搜索過程中要保證徹底避免進(jìn)入禁飛區(qū)。另外由于UAV 可以在不同高度下航行,因此本文不考慮UAV之間的碰撞規(guī)避問題。

    根據(jù)上述原則,本文將環(huán)境劃分為目標(biāo)信息未知環(huán)境,已知環(huán)境和禁飛區(qū),并對這3 種區(qū)域區(qū)別對待:UAV 主要針對未知環(huán)境展開搜索;UAV 應(yīng)當(dāng)盡量避免對已知環(huán)境進(jìn)行搜索(但在受到轉(zhuǎn)彎限制的情況下可以經(jīng)過);UAV 必須完全避開禁飛區(qū)。下面就路徑?jīng)Q策算法設(shè)計以及針對3 種區(qū)域的不同對策設(shè)計展開研究。

    2 協(xié)同搜索航路規(guī)劃

    2.1 環(huán)境模型

    本文使用笛卡爾柵格識別地圖描述環(huán)境,每一柵格被賦予一個值代表UAV 關(guān)于該區(qū)域目標(biāo)分布的感知能力。設(shè)環(huán)境E 為一個Lx×Ly的有界區(qū)域,將UAV 在t=kT 時刻對于(x,y)坐標(biāo)處的目標(biāo)分布的不確定性能力表示為η(x,y,k),且滿足η(x,y,k)∈[0,1].η(x,y,k)值越大,表示的不確定性越高,如果η(x,y,k)=1,表示UAV 在t 時刻完全不知道目標(biāo)在單元(x,y)的分布情況。

    2.2 UAV 模型

    假定所有UAV 以定常速度v 運動,θmax為UAV物理操縱性能限制最大轉(zhuǎn)彎角,UAV 的運動數(shù)學(xué)模型如下:

    當(dāng)某架UAV 經(jīng)過某區(qū)域(x,y)時,有

    式中:p 為機載傳感器對目標(biāo)的探測能力。如果對某單元多次搜索,η(x,y,k)趨向于0.

    2.3 路徑?jīng)Q策算法

    UAV 使用識別地圖儲存關(guān)于環(huán)境的不確定知識,通過數(shù)據(jù)融合算法載入UAV 傳感器,隨著識別地圖不斷更新,航路點決策函數(shù)利用地圖確定搜索特定環(huán)境區(qū)域的值,并且產(chǎn)生航路導(dǎo)引每架UAV 進(jìn)行搜索。這種航路應(yīng)該考慮到UAV 的操縱限制,即首先要求產(chǎn)生的航路物理可行,然后利用回報函數(shù)定義一種控制問題(該回報函數(shù)表征探索目標(biāo)的收益),最終選擇一條可產(chǎn)生最大回報的航路。算法設(shè)計如下:

    首先將路徑規(guī)劃問題時間離散化,即只允許UAV 在離散時間間隔T 內(nèi)做出決策。假定UAVi 在時刻t=kT 的位置為pi(k),將pi(k+1)的可能位置離散化為m 個點,每個點的位置可以表示為(k+1),j=1,2,…,m,則UAVi 在時刻t=kT 的所有可能點的位置可以表示為如下集合:(k+1)={p1i(k+1),…(k+1),…(k+1)}UAV 從pi(k)運動到下一位置pi(k+1)的距離為dt,并且應(yīng)滿足±θm的角度范圍限制??紤]到UAV 之間的數(shù)據(jù)傳輸延遲,UAVi 的開始q 個位置,pi(1),pi(2),…,pi(q)需要提前做出選擇。當(dāng)UAVi 在時刻t 位于位置pi(k)時,它已經(jīng)決定了下面的q 個位置:pi(k+1),pi(k+2),…,pi(k+q).當(dāng)UAV 到達(dá)pi(k)后,它開始選擇位置pi(k+q+1),將在時刻t=k+q+1 到達(dá)。只考慮UAV 3 種航路選擇狀態(tài):左轉(zhuǎn)θmax,直飛,右轉(zhuǎn)θmax,即m 值取3,當(dāng)q=3 時的航路搜索樹如圖2所示。

    3 搜索回報函數(shù)

    協(xié)同航路搜索的關(guān)鍵是設(shè)計一個搜索回報函數(shù),對每條可行路徑的回報值進(jìn)行評估。根據(jù)前面提出的協(xié)同航路規(guī)劃方法,每架UAV 使用回報函數(shù)J 選擇它的搜索航路:

    式中:J1,J2,J3為航路選擇點的對應(yīng)回報標(biāo)準(zhǔn);ωi為相應(yīng)的權(quán)重,0≤ωi≤1,且為重要性因子,γ≥1.收益標(biāo)準(zhǔn)的優(yōu)先性通過調(diào)整權(quán)重ωi實現(xiàn),航路選擇點的重要性通過調(diào)整γ 來體現(xiàn)。如果UAVi 在時刻t=kT 選擇第j∈[1,2,…,m]條路徑,則J1(i,j,k),J2(i,j,k)分別計算如下:

    圖2 3 步航路搜索樹示意圖Fig.2 Illustration of recursive 3-step planning tree

    式中:(x,y)為UAVi 沿第j 條航路飛行的航路點;η(x,y;k)為UAVi 在時刻k 關(guān)于點(x,y)的不確定值。從(7)式和(8)式可以看出,J1(i,j,k)表示當(dāng)前位置pi(k)與(k-1)之間不確定性減少的程度,J2(i,j,k)表示(k),(k),…(k)的平均不確定度。

    UAV 之間通過數(shù)據(jù)鏈傳輸實現(xiàn)信息共享,在UAV 相互之間距離很近,運動方向大致相同或相反的情況下容易出現(xiàn)航路重疊。為避免出現(xiàn)這種情況,將其他UAV 視為軟障礙,運用一種人工勢場算法使一定距離范圍內(nèi)的其他UAV 對UAVi 施加綜合抗力Fi(k),使UAVi 從(k+q+1)中選擇航路點pi(k+q+1),實現(xiàn)航路規(guī)避。

    在時刻k 由其他所有UAV 作用到UAVi 上的綜合抗力可以表示為

    其中:

    式中:Fi(k)為在時刻t=kT 由相鄰UAV 作用到UAVi 上的綜合抗力;Rij為UAVi 與UAVj 的距離;Rij為對應(yīng)的單位向量。k1取值對應(yīng)于距離Rij為0時的抗力大小。設(shè)計參數(shù)μ >0,對應(yīng)于抗力隨距離增加而減少的速率。Rmax、βmax、φmax分別為UAVj 與UAVi 之間產(chǎn)生抗力作用的最大允許距離、最大允許角度和最大允許方位角之差。如圖3所示。

    式中:k2為非負(fù)參數(shù);αi(j,k)為綜合抗力Fi(k)與(k+q+1)中每一可行路徑的方位差。從協(xié)同的觀點來看,顯然應(yīng)該選擇具有最小αi(j,k)的航路。

    將(6)式標(biāo)準(zhǔn)化,有

    Jn為標(biāo)準(zhǔn)化代價函數(shù),計算如下

    圖3 UAV 抗力產(chǎn)生條件示意圖Fig.3 Sketch of rivaling force conditions between UAVs

    4 禁飛區(qū)回避決策設(shè)計

    針對禁飛區(qū)回避問題一種可能的方法是將禁飛區(qū)內(nèi)的搜索回報函數(shù)設(shè)定為最小值。但是這并不能保證UAV 不會進(jìn)入禁飛區(qū)。如圖4中當(dāng)UAV 以航向角θ=90°飛到節(jié)點b2時,受最大轉(zhuǎn)彎角的限制,當(dāng)可選節(jié)點均位于禁飛區(qū)時,UAV 將無法避免進(jìn)入禁飛區(qū)。下面介紹一種新的解決辦法,以θmax=45°為例進(jìn)行說明。

    解決的思路是在UAV 要選擇的位置還沒有進(jìn)入禁飛區(qū)時提前做出選擇,當(dāng)pi(k+q+1)與禁飛區(qū)的邊界距離小于1 個步長dt 的距離時,如圖4陰影部分所示,UAV 開始從b1,b2,b3三個位置中做出選擇。圖4中UAV 從a1進(jìn)入b2時,由于受到轉(zhuǎn)彎限制,UAV 無法避開禁飛區(qū),因此b2為不可行點,即使b2點具有比b1和b3更高的收益也必須放棄,而此時b1和b3為可行點,UAV 可分別到達(dá)點c1和c7,展開下一步搜索。需要注意的是當(dāng)UAV 從b1和b3的正下方進(jìn)入時,這2 點又都變成了不可行點,因此判斷某點是否可行點取決于該點的位置和UAV 進(jìn)入該點的角度。針對圖4中的矩形禁飛區(qū)進(jìn)一步分析可發(fā)現(xiàn),當(dāng)x1≤xb≤x2并且y1- dt≤yb<y1時,如果(xb,yb)處UAV 的航向角為90°,則該點為不可行點,航向角為45°時則該點為可行點。同理可對陰影框的其他部分進(jìn)行分析。為了保證下一步搜索時UAV 不會進(jìn)入禁飛區(qū),將所有位于禁飛區(qū)內(nèi)的選擇點不論航向角大小均列為不可行點。

    設(shè)M={(xi,yi),θi}為陰影區(qū)內(nèi)的所有點和所有航向角的集合;Mf為陰影區(qū)內(nèi)所有可行點和可行航向角的集合,則在圖4中Mf應(yīng)滿足(14)式~(21)式:

    需要說明的是,(14)式~(21)式只針對圖4.隨著禁飛區(qū)位置和形狀的不同,Mf會有所變化,需要根據(jù)具體情況分析,但分析策略類似。

    UAV 整個搜索程序如下:

    1)任務(wù)參數(shù)初始化。

    2)for k=1 to n.

    3)if pi(k+q+1)∈E.

    4)if (pi(k+q+1),θi(k+q+1))?M.

    5)計算pi(k+q+1)的m 個可能位置處的回報函數(shù)J,選擇J 值最大的航路。

    6)else 根據(jù)pi(k+q+1)的m 個可能位置是否滿足Mf中可行點的條件從中選取m'個可行位置和對應(yīng)的可行航向角,選擇J 值最大的航路。

    7)end if.

    8)else 執(zhí)行轉(zhuǎn)彎程序,使UAV 回到環(huán)境中。

    9)end if.

    10)k=k+1.

    11)end for.

    5 仿真實驗分析

    設(shè)E(x,y)為待偵察區(qū)域,50 km≤x≤200 km,50 km≤y≤200 km,A 為根據(jù)先驗知識不需要進(jìn)行偵察的區(qū)域,B 為禁飛區(qū)。假定在每一樣本時間內(nèi),UAV 只允許在-45°~45°之間的角度進(jìn)行直飛,左轉(zhuǎn)或右轉(zhuǎn)。

    實驗1 設(shè)需要進(jìn)行重點偵察的區(qū)域F(x,y)坐標(biāo)為:130 km≤x≤180 km,60 km≤y≤110 km.剩余部分為完全不確定區(qū)域,即η(x,y)=1,ω1=1/3,ω2=1/3,ω3=1/3,(x,y)∈F 時,γ=3,(x,y)?F時,γ=1.機載傳感器對目標(biāo)的探測能力p=0.6.5 架UAV 從不同基地起飛開始執(zhí)行偵察任務(wù),其起飛點坐標(biāo)分別為(50,50),(80,50),(110,50),(140,50),(170,50),運行60 步長的仿真結(jié)果如圖5所示。

    圖5 設(shè)置重點搜索區(qū)域的5 架UAV 協(xié)同搜索仿真結(jié)果Fig.5 Simulation results of cooperative search with a key search zone for five UAVs

    由圖5可知,UAV1 有一段航路經(jīng)過A 區(qū),而UAV1 和UAV2 盡管一度接近禁飛區(qū)但最后都順利避開。由于對F(x,y)搜索能獲得更高的收益,該區(qū)航路明顯比較密集。

    實驗2 為了便于與其他算法進(jìn)行比較,設(shè)γ為常數(shù),取值1,其他參數(shù)設(shè)定同試驗1.試驗分2部分進(jìn)行:首先運行80 步長,然后假定UAV1、UAV4 和UAV5 設(shè)備故障返航,分別對UAV2 和UAV3 運行20 步長,分別運用Zamboni 搜索、隨機搜索和本文的協(xié)同搜索算法進(jìn)行仿真,仿真結(jié)果分別如圖6~圖8所示(UAV2 和UAV3 的后半段航跡用星號線表示)。

    圖6 5 架UAV 的Zamboni 搜索仿真結(jié)果Fig.6 Simulation results of Zamboni search for five UAVs

    圖7 5 架UAV 的隨機搜索仿真結(jié)果Fig.7 Simulation results of random search for five UAVs

    由圖6和圖7可以看出,Zamboni 搜索和隨機搜索的航跡不會受到區(qū)域A 和B 的影響,而圖8中的協(xié)同搜索只有一架UAV 的航跡經(jīng)過區(qū)域A,所有UAV 都避開了區(qū)域B.另外,Zamboni 覆蓋算法和協(xié)同搜索算法幾乎沒有路徑重疊,而隨機搜索出現(xiàn)了多處多架UAV 航跡段重疊的情況。在3 架UAV 故障后,Zamboni 算法中的UAV2 和UAV3 仍然按照既定航線搜索,故障UAV 未完成的任務(wù)將得不到執(zhí)行,算法不具有魯棒性。圖7中UAV3 后20 步長中的航跡與UAV1 出現(xiàn)了重疊,這也降低了搜索收益。而從圖8可以看出UAV2 和UAV3 后20 步長中仍能保持之前的協(xié)同搜索策略,并且依然能夠避開禁飛區(qū)。

    圖9為搜索效率隨仿真步長的變化情況。由于多處航路經(jīng)過區(qū)域A 和B,造成Zamboni 算法的偵察效率比協(xié)同搜索低,而部分UAV 故障導(dǎo)致在仿真后半段搜索效率幾乎沒有提高;隨機搜索不僅多處航路經(jīng)過區(qū)域A 和B,而且出現(xiàn)多處航路重疊,所以搜索效率也不高;與前2 種非協(xié)同方法相比,協(xié)同搜索可以使搜索環(huán)境的不確定度減少的更快。

    圖8 5 架UAV 的協(xié)同搜索仿真結(jié)果Fig.8 Simulation results of cooperative search for five UAVs

    圖9 5 架UAV 在3 種搜索樣式下的搜索效率Fig.9 Search efficiencies in three search patterns for five UAVs

    6 結(jié)論

    通過對目標(biāo)信息未知環(huán)境,已知環(huán)境和禁飛區(qū)實施不同對策,解決了UAV 在不確定環(huán)境中的協(xié)同搜索問題。仿真實驗結(jié)果表明,本文提出的協(xié)同搜索算法對未知環(huán)境可以實施重點偵察,在機載燃油充足的情況下也可以展開覆蓋搜索。對已知環(huán)境能夠?qū)崿F(xiàn)較好的規(guī)避,并完全避免了UAV 進(jìn)入禁飛區(qū)。UAV 在整個仿真過程中都能夠滿足機動限制,并能夠在友機故障的情況下繼續(xù)實現(xiàn)協(xié)同。與3 種非協(xié)同搜索算法相比,本文提出的協(xié)同搜索算法具有更高的搜索效率。

    References)

    [1] Schoenwald D A.AUVs:in space,air,water,and on the ground[J].IEEE Control Syst,2000,20(6):15-18.

    [2] Chandler P,Rasmussen S,Pachter M.UAV cooperative path planning[C]∥Proceedings of AIAA Guidance,Navigation,and Control Conference and Exhibit.Denver:American Institute of Aeronautics and Astronautics,2000:1255-1265.

    [3] Chandler P,Pachter M,Rasmussen S.Hierarchical control for autonomous teams[C]∥Proceedings of AIAA Guidance,Navigation,and Control Conference and Exhibit.Montreal:American Institute of Aeronautics and Astronautics,2001:632-642.

    [4] Svennebring J,Koenig S.Building terrain-covering ant robots:a feasibility study[J].Autonomous Robots,2004,16(3):313-332.

    [5] Wagner I A,Lindenbaum M,Bruckstein A M.MAC versus PC:determinism and randomness as complementary approaches to robotic exploration of continuous unknown domains[J].International Journal of Robotics Research,2000,19(1):12-31.

    [6] Yang S,Luo C.A neural network approach to complete coverage path planning[J].IEEE Transactions on Systems,Man and Cybernetics,2004,34(1):718-725.

    [7] Choset H.Coverage for robotics:a survey of recent results[J].Annals of Mathematics and Artificial Intelligence,2001,31:113-126.

    [8] Sujit P B,Ghose D.Search using multiple UAVs with flight time constraints[J].IEEE Transactions on Aerospace and Electronic Systems,2004,40(2):491-510.

    [9] Yang Y.Cooperative search by uninhabited air vehicles in dynamic environment[D].Cincinnati:University of Cincinnati,2005.

    [10] 柳長安,王和平,李為吉.無人機的偵察航路規(guī)劃[J].西北工業(yè)大學(xué)學(xué)報,2003,21(4):490-493.LIU Chang-an,WANG He-ping,LI Wei-ji.On path planning for more efficient reconnaissance of UAV[J].Journal of Northwestern Polytechnical University,2003,21(4):490-493.(in Chinese)

    [11] Polycarpou M,Yang Y,Passino K.A cooperative search framework for distributed agents[C]∥Proceedings of the 2001 IEEE International Symposium on Intelligent Control.Mexico:Institute of Electrical and Electronics Engineers,2001:1-6.

    [12] Butenko S,Murphey R,Pardalos P.Cooperative control:models,applications and algorithms[M].Dordrecht:Kluwer Academic Publishers,2003:283-321.

    猜你喜歡
    禁飛區(qū)抗力航路
    橋梁低樁承臺基底豎向抗力作用效應(yīng)的計算研究
    交通科技(2021年4期)2021-09-03 09:47:22
    基于實時航路的PFD和ND的仿真研究
    大疆更新多邊形禁飛區(qū)策略
    航空模型(2017年7期)2017-08-15 10:24:44
    巖塊的彈性模量及巖體單位彈性抗力系數(shù)的確定方法
    應(yīng)召反潛時無人機監(jiān)聽航路的規(guī)劃
    托勒密世界地圖與新航路的開辟
    基于Event改進(jìn)模型的交叉航路碰撞風(fēng)險評估
    30MnSi熱變形抗力的數(shù)學(xué)模型
    上海金屬(2014年2期)2014-12-18 06:52:51
    人臉禁飛區(qū)
    自然與科技(2013年6期)2013-04-29 00:44:03
    驗證荷載對結(jié)構(gòu)抗力分布的時效影響
    免费高清在线观看日韩| 青春草亚洲视频在线观看| 精品卡一卡二卡四卡免费| 麻豆乱淫一区二区| 久久久精品免费免费高清| 男女之事视频高清在线观看 | 久久热在线av| 99热全是精品| 亚洲久久久国产精品| 色精品久久人妻99蜜桃| 制服人妻中文乱码| 久久久精品区二区三区| 一本一本久久a久久精品综合妖精| 最近最新中文字幕大全免费视频 | 国产黄频视频在线观看| 激情五月婷婷亚洲| 国产日韩欧美视频二区| 午夜福利,免费看| 麻豆av在线久日| 亚洲精品国产一区二区精华液| 好男人视频免费观看在线| 一区二区三区激情视频| 亚洲男人天堂网一区| 亚洲一区二区三区欧美精品| 黑人欧美特级aaaaaa片| 悠悠久久av| 欧美黄色片欧美黄色片| 亚洲色图 男人天堂 中文字幕| av天堂久久9| 亚洲精品一卡2卡三卡4卡5卡 | 亚洲成人国产一区在线观看 | 成年av动漫网址| 蜜桃在线观看..| 国产亚洲欧美在线一区二区| 18禁国产床啪视频网站| 国产99久久九九免费精品| 9热在线视频观看99| 亚洲精品国产区一区二| av在线老鸭窝| 九草在线视频观看| 欧美在线一区亚洲| 亚洲欧美成人综合另类久久久| 一级毛片电影观看| 日本av手机在线免费观看| 王馨瑶露胸无遮挡在线观看| 免费在线观看日本一区| 9色porny在线观看| 国产成人a∨麻豆精品| 2021少妇久久久久久久久久久| 狂野欧美激情性bbbbbb| 国产无遮挡羞羞视频在线观看| 日本五十路高清| 成人午夜精彩视频在线观看| 在线天堂中文资源库| 男女之事视频高清在线观看 | av又黄又爽大尺度在线免费看| 国产精品欧美亚洲77777| 91精品伊人久久大香线蕉| 999久久久国产精品视频| 19禁男女啪啪无遮挡网站| 最黄视频免费看| 啦啦啦在线免费观看视频4| 欧美成人精品欧美一级黄| 一二三四社区在线视频社区8| 国产色视频综合| 我的亚洲天堂| 日韩制服丝袜自拍偷拍| 亚洲国产精品一区二区三区在线| 午夜久久久在线观看| 97人妻天天添夜夜摸| 精品人妻在线不人妻| 久久天躁狠狠躁夜夜2o2o | 视频区欧美日本亚洲| 亚洲欧洲国产日韩| 国产精品久久久久成人av| 十八禁高潮呻吟视频| 丝袜美腿诱惑在线| 伦理电影免费视频| 亚洲精品国产av蜜桃| 国产三级黄色录像| 首页视频小说图片口味搜索 | 婷婷色麻豆天堂久久| 国产高清国产精品国产三级| 欧美成人精品欧美一级黄| 老汉色av国产亚洲站长工具| 欧美人与性动交α欧美精品济南到| xxxhd国产人妻xxx| 久久久久久久国产电影| 免费女性裸体啪啪无遮挡网站| av又黄又爽大尺度在线免费看| 亚洲精品日本国产第一区| 久久精品熟女亚洲av麻豆精品| cao死你这个sao货| 久久免费观看电影| 999久久久国产精品视频| 9热在线视频观看99| 国产精品一区二区精品视频观看| 国产视频一区二区在线看| 久久久久久免费高清国产稀缺| 亚洲精品国产av蜜桃| 国产精品香港三级国产av潘金莲 | 亚洲免费av在线视频| 成人三级做爰电影| 精品一区二区三区av网在线观看 | 国产精品.久久久| 亚洲欧美精品自产自拍| 69精品国产乱码久久久| 亚洲欧洲日产国产| 国语对白做爰xxxⅹ性视频网站| 婷婷成人精品国产| 久久精品熟女亚洲av麻豆精品| 久久久久久久久久久久大奶| 亚洲 欧美一区二区三区| 亚洲 国产 在线| 成人亚洲欧美一区二区av| 一本色道久久久久久精品综合| 国产片特级美女逼逼视频| 久久久精品免费免费高清| 在线观看免费日韩欧美大片| 大片免费播放器 马上看| 国产精品久久久av美女十八| 美女大奶头黄色视频| 精品第一国产精品| 亚洲国产精品国产精品| 中国美女看黄片| 狠狠婷婷综合久久久久久88av| bbb黄色大片| 一级毛片我不卡| 美女中出高潮动态图| 男的添女的下面高潮视频| 制服人妻中文乱码| 久久久国产欧美日韩av| 最新的欧美精品一区二区| xxxhd国产人妻xxx| 巨乳人妻的诱惑在线观看| 亚洲av国产av综合av卡| 一区二区三区四区激情视频| 国产亚洲精品第一综合不卡| 国产成人影院久久av| 啦啦啦在线免费观看视频4| 国产野战对白在线观看| 91国产中文字幕| 国产成人免费无遮挡视频| 黑人巨大精品欧美一区二区蜜桃| 天天躁夜夜躁狠狠久久av| 久久综合国产亚洲精品| 精品亚洲成国产av| 欧美精品人与动牲交sv欧美| 黄色怎么调成土黄色| 精品亚洲乱码少妇综合久久| 欧美亚洲 丝袜 人妻 在线| 国产成人免费无遮挡视频| 国产精品麻豆人妻色哟哟久久| 美女大奶头黄色视频| 国产成人av激情在线播放| 韩国精品一区二区三区| 18在线观看网站| 少妇猛男粗大的猛烈进出视频| 国产一区有黄有色的免费视频| 这个男人来自地球电影免费观看| av国产久精品久网站免费入址| 国产精品香港三级国产av潘金莲 | 午夜福利视频在线观看免费| 成人手机av| 亚洲成人手机| kizo精华| 18禁裸乳无遮挡动漫免费视频| 久久久久国产精品人妻一区二区| 久久久国产精品麻豆| 婷婷丁香在线五月| 不卡av一区二区三区| 亚洲国产精品国产精品| 九草在线视频观看| 久久精品亚洲av国产电影网| 日韩伦理黄色片| 操美女的视频在线观看| 黑人欧美特级aaaaaa片| 国产精品99久久99久久久不卡| 女人高潮潮喷娇喘18禁视频| 国产一区二区激情短视频 | 后天国语完整版免费观看| 国产精品三级大全| 美女国产高潮福利片在线看| 亚洲精品一区蜜桃| 黄色视频不卡| 成人手机av| 欧美久久黑人一区二区| 一边摸一边抽搐一进一出视频| 国产精品香港三级国产av潘金莲 | 欧美国产精品va在线观看不卡| 性色av一级| 啦啦啦在线观看免费高清www| 成人影院久久| 亚洲精品第二区| 久久久精品免费免费高清| 18在线观看网站| 中国国产av一级| 午夜两性在线视频| 国产又色又爽无遮挡免| 一级片免费观看大全| 免费久久久久久久精品成人欧美视频| 久久国产精品影院| 国产午夜精品一二区理论片| 亚洲人成电影观看| 中文字幕最新亚洲高清| 热re99久久精品国产66热6| 亚洲精品美女久久久久99蜜臀 | 女人爽到高潮嗷嗷叫在线视频| 啦啦啦啦在线视频资源| 精品一区二区三卡| 色网站视频免费| 亚洲av男天堂| 老司机深夜福利视频在线观看 | 亚洲精品第二区| 人妻 亚洲 视频| 大话2 男鬼变身卡| 免费在线观看黄色视频的| 少妇人妻久久综合中文| av天堂久久9| 亚洲 欧美一区二区三区| 999久久久国产精品视频| 欧美日本中文国产一区发布| 欧美精品一区二区免费开放| 9色porny在线观看| 韩国高清视频一区二区三区| 亚洲视频免费观看视频| 亚洲av国产av综合av卡| 男女边吃奶边做爰视频| 又紧又爽又黄一区二区| 国产成人一区二区在线| 婷婷色综合大香蕉| 日韩一区二区三区影片| 国产福利在线免费观看视频| 国产免费一区二区三区四区乱码| a级片在线免费高清观看视频| 99热网站在线观看| 女人被躁到高潮嗷嗷叫费观| 欧美精品一区二区大全| 久久午夜综合久久蜜桃| 老汉色∧v一级毛片| 欧美激情极品国产一区二区三区| 国产一区二区在线观看av| 精品卡一卡二卡四卡免费| 日本五十路高清| 大陆偷拍与自拍| 国产成人免费观看mmmm| 日韩免费高清中文字幕av| 免费黄频网站在线观看国产| 国产亚洲精品久久久久5区| 七月丁香在线播放| netflix在线观看网站| 丁香六月天网| 最近最新中文字幕大全免费视频 | 免费少妇av软件| 男女边摸边吃奶| 我的亚洲天堂| 成人国产一区最新在线观看 | 欧美黄色片欧美黄色片| 两性夫妻黄色片| 久久青草综合色| 成人国产一区最新在线观看 | 热99国产精品久久久久久7| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲欧美一区二区三区黑人| 中文精品一卡2卡3卡4更新| 手机成人av网站| 国产午夜精品一二区理论片| 免费看十八禁软件| 欧美日韩视频高清一区二区三区二| 欧美亚洲日本最大视频资源| 亚洲精品美女久久av网站| 日韩制服骚丝袜av| 夜夜骑夜夜射夜夜干| 国产xxxxx性猛交| 免费高清在线观看视频在线观看| 国产精品一区二区免费欧美 | 97在线人人人人妻| 99国产综合亚洲精品| 久久久久久亚洲精品国产蜜桃av| 精品久久久精品久久久| 午夜影院在线不卡| 久久精品国产亚洲av高清一级| 国产亚洲午夜精品一区二区久久| 久热爱精品视频在线9| 午夜福利,免费看| 精品熟女少妇八av免费久了| 色播在线永久视频| 国产日韩欧美亚洲二区| 成年美女黄网站色视频大全免费| 国产精品 欧美亚洲| 妹子高潮喷水视频| 性高湖久久久久久久久免费观看| 成人影院久久| 国产欧美日韩一区二区三 | 国产成人av激情在线播放| 人妻 亚洲 视频| xxx大片免费视频| 亚洲图色成人| 一区二区av电影网| 狠狠精品人妻久久久久久综合| 汤姆久久久久久久影院中文字幕| 亚洲欧美日韩另类电影网站| 天天操日日干夜夜撸| 女人爽到高潮嗷嗷叫在线视频| 欧美激情极品国产一区二区三区| 两个人看的免费小视频| 亚洲人成网站在线观看播放| 丰满人妻熟妇乱又伦精品不卡| 国产日韩欧美在线精品| 少妇的丰满在线观看| 丝瓜视频免费看黄片| 国产黄频视频在线观看| 一本久久精品| www.熟女人妻精品国产| 成年人午夜在线观看视频| 亚洲成人免费av在线播放| 成人免费观看视频高清| 熟女少妇亚洲综合色aaa.| 妹子高潮喷水视频| 热99久久久久精品小说推荐| 男人操女人黄网站| 精品少妇久久久久久888优播| 久久人人97超碰香蕉20202| 51午夜福利影视在线观看| 国产人伦9x9x在线观看| www.自偷自拍.com| 亚洲av欧美aⅴ国产| 青春草亚洲视频在线观看| 亚洲专区中文字幕在线| 久久人妻熟女aⅴ| 人人澡人人妻人| 韩国精品一区二区三区| 亚洲七黄色美女视频| 欧美久久黑人一区二区| av天堂在线播放| 国产黄色免费在线视频| 少妇猛男粗大的猛烈进出视频| 国产又爽黄色视频| 亚洲熟女精品中文字幕| 欧美人与性动交α欧美软件| 九草在线视频观看| 精品第一国产精品| 精品视频人人做人人爽| 国产国语露脸激情在线看| 久久精品亚洲熟妇少妇任你| 日韩熟女老妇一区二区性免费视频| 久久久久久久精品精品| 精品少妇黑人巨大在线播放| 一级毛片女人18水好多 | 国产高清国产精品国产三级| 久久av网站| 青春草亚洲视频在线观看| 1024香蕉在线观看| 国产成人av激情在线播放| 尾随美女入室| xxxhd国产人妻xxx| 国产亚洲欧美在线一区二区| 99久久99久久久精品蜜桃| 色94色欧美一区二区| 99国产精品99久久久久| 99精国产麻豆久久婷婷| 免费在线观看视频国产中文字幕亚洲 | 午夜福利乱码中文字幕| 男女国产视频网站| 黄片小视频在线播放| 亚洲国产精品一区三区| 亚洲,欧美,日韩| 精品第一国产精品| 亚洲国产日韩一区二区| 在线观看免费高清a一片| 国产三级黄色录像| 丰满少妇做爰视频| 男女免费视频国产| 97人妻天天添夜夜摸| a级毛片在线看网站| 黄色视频在线播放观看不卡| 男女边摸边吃奶| 亚洲欧美日韩另类电影网站| 如日韩欧美国产精品一区二区三区| 亚洲精品国产区一区二| 黄色片一级片一级黄色片| 天天躁夜夜躁狠狠躁躁| 国产国语露脸激情在线看| 日韩制服丝袜自拍偷拍| 亚洲精品乱久久久久久| 精品一品国产午夜福利视频| 欧美大码av| 一本一本久久a久久精品综合妖精| 青青草视频在线视频观看| 国产精品二区激情视频| 亚洲精品av麻豆狂野| www.av在线官网国产| 免费一级毛片在线播放高清视频 | 在线精品无人区一区二区三| 日日夜夜操网爽| 亚洲中文av在线| 国产成人精品久久二区二区免费| 国产伦人伦偷精品视频| 亚洲午夜精品一区,二区,三区| 亚洲av欧美aⅴ国产| 97精品久久久久久久久久精品| 免费看十八禁软件| 丝袜在线中文字幕| 亚洲欧洲国产日韩| 国产一区有黄有色的免费视频| 啦啦啦在线观看免费高清www| 在线观看免费视频网站a站| www.熟女人妻精品国产| 99国产精品免费福利视频| 啦啦啦在线免费观看视频4| 人人妻人人澡人人爽人人夜夜| 天天躁日日躁夜夜躁夜夜| 男女之事视频高清在线观看 | 一级毛片黄色毛片免费观看视频| 亚洲欧美中文字幕日韩二区| 国产成人影院久久av| 另类精品久久| 久久这里只有精品19| 搡老乐熟女国产| 宅男免费午夜| 久久 成人 亚洲| 国产精品一区二区在线观看99| 99久久99久久久精品蜜桃| 肉色欧美久久久久久久蜜桃| 亚洲欧美精品综合一区二区三区| 一区二区三区精品91| 少妇 在线观看| 香蕉国产在线看| 国产爽快片一区二区三区| 七月丁香在线播放| 午夜免费观看性视频| 成人影院久久| 欧美激情极品国产一区二区三区| 一本久久精品| 一级,二级,三级黄色视频| 啦啦啦啦在线视频资源| 国产麻豆69| 国产深夜福利视频在线观看| 国产精品人妻久久久影院| 一区福利在线观看| 亚洲情色 制服丝袜| 美女高潮到喷水免费观看| av国产久精品久网站免费入址| 男女免费视频国产| 国产精品一二三区在线看| 欧美日韩福利视频一区二区| 午夜福利视频在线观看免费| 国产片内射在线| 99re6热这里在线精品视频| 日日摸夜夜添夜夜爱| 天天添夜夜摸| 日韩 亚洲 欧美在线| 制服诱惑二区| 国产亚洲精品第一综合不卡| xxxhd国产人妻xxx| 国产精品一区二区在线观看99| 热99国产精品久久久久久7| 麻豆av在线久日| 国产一级毛片在线| 久久久久久亚洲精品国产蜜桃av| 午夜福利一区二区在线看| 天堂8中文在线网| 日韩熟女老妇一区二区性免费视频| 无限看片的www在线观看| 久久久亚洲精品成人影院| 热re99久久精品国产66热6| 国产成人av教育| 无遮挡黄片免费观看| √禁漫天堂资源中文www| 中文字幕人妻丝袜一区二区| 女人爽到高潮嗷嗷叫在线视频| 欧美精品人与动牲交sv欧美| 黄色 视频免费看| 国产成人免费观看mmmm| 赤兔流量卡办理| 啦啦啦在线免费观看视频4| 中文字幕人妻丝袜制服| 亚洲欧美精品综合一区二区三区| 中文字幕人妻熟女乱码| 亚洲国产毛片av蜜桃av| 国产一区二区 视频在线| 国产精品一区二区在线观看99| 老司机影院成人| 国产精品国产三级专区第一集| 91字幕亚洲| 一级片'在线观看视频| 免费在线观看完整版高清| 大片免费播放器 马上看| 人成视频在线观看免费观看| 亚洲精品国产区一区二| 一本—道久久a久久精品蜜桃钙片| 国产av一区二区精品久久| 国产99久久九九免费精品| 免费女性裸体啪啪无遮挡网站| 秋霞在线观看毛片| 亚洲av成人精品一二三区| 日韩电影二区| 啦啦啦中文免费视频观看日本| 你懂的网址亚洲精品在线观看| 国产亚洲av高清不卡| 亚洲av成人精品一二三区| 日韩精品免费视频一区二区三区| 校园人妻丝袜中文字幕| 亚洲欧美精品自产自拍| 婷婷成人精品国产| 国产一区亚洲一区在线观看| 亚洲人成网站在线观看播放| 亚洲中文日韩欧美视频| a级毛片在线看网站| 天天躁夜夜躁狠狠躁躁| 欧美中文综合在线视频| 女性生殖器流出的白浆| 国产精品一区二区精品视频观看| 精品国产乱码久久久久久小说| 久久国产精品男人的天堂亚洲| 黑人欧美特级aaaaaa片| 你懂的网址亚洲精品在线观看| 日韩一区二区三区影片| 一级毛片 在线播放| a级片在线免费高清观看视频| 热re99久久国产66热| 丝袜美腿诱惑在线| 欧美性长视频在线观看| 大型av网站在线播放| 99国产精品一区二区三区| 最近最新中文字幕大全免费视频 | 亚洲av电影在线进入| 久久 成人 亚洲| 国产男女超爽视频在线观看| 精品国产国语对白av| 欧美日韩综合久久久久久| 午夜日韩欧美国产| 亚洲九九香蕉| 老司机亚洲免费影院| 天天操日日干夜夜撸| 亚洲精品一卡2卡三卡4卡5卡 | 国产一区二区激情短视频 | 国产成人a∨麻豆精品| 色94色欧美一区二区| 91九色精品人成在线观看| 免费观看人在逋| 久久久久网色| 老鸭窝网址在线观看| 亚洲精品成人av观看孕妇| 久久99精品国语久久久| 国产亚洲av高清不卡| 午夜av观看不卡| 成年av动漫网址| 女人精品久久久久毛片| 亚洲第一av免费看| 熟女少妇亚洲综合色aaa.| 欧美变态另类bdsm刘玥| 精品第一国产精品| 岛国毛片在线播放| 丰满少妇做爰视频| 亚洲欧美一区二区三区黑人| 国产精品99久久99久久久不卡| 欧美亚洲 丝袜 人妻 在线| 欧美亚洲日本最大视频资源| 视频区图区小说| 日本91视频免费播放| 91麻豆av在线| 久久毛片免费看一区二区三区| 精品一区二区三区av网在线观看 | 一区福利在线观看| 热re99久久国产66热| 18在线观看网站| tube8黄色片| 脱女人内裤的视频| 伊人亚洲综合成人网| 韩国精品一区二区三区| 成人18禁高潮啪啪吃奶动态图| 免费女性裸体啪啪无遮挡网站| 老司机靠b影院| 高清av免费在线| 国产成人精品在线电影| 日韩一区二区三区影片| 国产精品人妻久久久影院| 国产在线免费精品| 久久久久国产精品人妻一区二区| 女人精品久久久久毛片| 一二三四社区在线视频社区8| 亚洲欧美中文字幕日韩二区| 菩萨蛮人人尽说江南好唐韦庄| 亚洲,欧美精品.| 久久精品亚洲熟妇少妇任你| 日本一区二区免费在线视频| 国产一卡二卡三卡精品| 欧美国产精品va在线观看不卡| 国产欧美日韩综合在线一区二区| 午夜免费观看性视频| 欧美乱码精品一区二区三区| 成人黄色视频免费在线看| 国产精品香港三级国产av潘金莲 | a级毛片在线看网站| videos熟女内射| 男女高潮啪啪啪动态图| 久久久久精品人妻al黑| 国产在视频线精品| 午夜福利,免费看| 午夜av观看不卡| 晚上一个人看的免费电影| 亚洲伊人久久精品综合| 天天影视国产精品| 青春草亚洲视频在线观看| 九草在线视频观看| 丰满少妇做爰视频| 男的添女的下面高潮视频| 妹子高潮喷水视频| 少妇人妻久久综合中文| 一区福利在线观看| 精品少妇久久久久久888优播| 99精国产麻豆久久婷婷| 欧美精品亚洲一区二区| 亚洲一卡2卡3卡4卡5卡精品中文|