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

    基于改進(jìn)人工勢(shì)場(chǎng)法的多機(jī)改航路徑規(guī)劃

    2013-11-04 02:38:50王興隆徐肖豪馮江然
    飛行力學(xué) 2013年4期
    關(guān)鍵詞:勢(shì)場(chǎng)步數(shù)航空器

    王興隆, 徐肖豪, 馮江然

    (中國(guó)民航大學(xué) 空中交通管理研究基地, 天津 300300)

    基于改進(jìn)人工勢(shì)場(chǎng)法的多機(jī)改航路徑規(guī)劃

    王興隆, 徐肖豪, 馮江然

    (中國(guó)民航大學(xué) 空中交通管理研究基地, 天津 300300)

    惡劣天氣是航空運(yùn)輸中造成航班延誤的主要因素,實(shí)施改航策略是減少延誤的一種有效方法,多航空器路徑規(guī)劃是改航策略的關(guān)鍵環(huán)節(jié)。針對(duì)惡劣天氣造成大量航班延誤情況,提出了采用人工勢(shì)場(chǎng)法進(jìn)行多機(jī)改航路徑規(guī)劃。根據(jù)航空器所受引力和斥力建立了傳統(tǒng)人工勢(shì)場(chǎng)模型,為解決在改航期間航空器受力出現(xiàn)局部最小值的問(wèn)題,在航空器之間、航空器與惡劣天氣之間引入不同的旋轉(zhuǎn)控制力。仿真結(jié)果表明,改進(jìn)的人工勢(shì)場(chǎng)法規(guī)劃的路徑更短、更平滑,沒(méi)有大角度轉(zhuǎn)彎和反復(fù)振蕩的情況,且計(jì)算效率提高了33%。

    航空運(yùn)輸; 改航路徑規(guī)劃; 人工勢(shì)場(chǎng)法; 多航空器

    0 引言

    導(dǎo)致航班延誤的因素有很多,其中不利于飛行的天氣,如大霧、雷暴、積冰等是造成航班延誤的主要原因。根據(jù)民航局公布的資料統(tǒng)計(jì),由于天氣原因造成的機(jī)場(chǎng)或空域容量急劇下降而直接或間接導(dǎo)致航班延誤的比例約占不正常航班的三分之一。

    當(dāng)出現(xiàn)惡劣天氣時(shí),實(shí)施改航策略是減少航班延誤的一種有效途徑,改航策略是航班臨時(shí)選擇一條不受影響的航線繞過(guò)飛行受限區(qū)域,航班改航路徑規(guī)劃是改航策略的一個(gè)關(guān)鍵問(wèn)題。國(guó)外對(duì)該領(lǐng)域的研究已有十幾年的時(shí)間,具有代表性的算法主要有:基于網(wǎng)格的改航路徑規(guī)劃法[1];基于多邊形的改航路徑規(guī)劃法[2];基于標(biāo)準(zhǔn)進(jìn)離場(chǎng)程序的改航路徑規(guī)劃法[3];基于自由飛行的改航路徑規(guī)劃法[4]。我國(guó)對(duì)于改航問(wèn)題的研究開始于2000年,目前仍處于研究的初級(jí)階段。我國(guó)學(xué)者從2001年開始對(duì)航班改航路徑規(guī)劃問(wèn)題進(jìn)行研究,提出了基于已有航路點(diǎn)的A*搜索算法[5];文獻(xiàn)[6]提出了基于幾何算法的改航路徑規(guī)劃方法,該方法考慮了最短航段距離、轉(zhuǎn)彎角度和轉(zhuǎn)彎點(diǎn)數(shù)量等因素對(duì)改航路徑的影響;文獻(xiàn)[7]考慮了多個(gè)限制因素,建立了基于人工勢(shì)場(chǎng)算法的動(dòng)態(tài)改航路徑規(guī)劃模型,實(shí)現(xiàn)了改航路徑的快速實(shí)時(shí)規(guī)劃。這些研究只考慮了單航空器(也稱單機(jī))改航路徑的規(guī)劃,在實(shí)際運(yùn)行中,惡劣天氣勢(shì)必影響多個(gè)航空器的正常運(yùn)行,因此,不僅要考慮航空器避讓惡劣天氣,也要分析航空器之間的飛行沖突。針對(duì)這一問(wèn)題,本文研究了惡劣天氣下的多航空器改航路徑規(guī)劃與優(yōu)化。

    1傳統(tǒng)人工勢(shì)場(chǎng)法的多機(jī)改航路徑

    規(guī)劃

    采用傳統(tǒng)人工勢(shì)場(chǎng)法建立人工勢(shì)場(chǎng)模型時(shí),航空器被簡(jiǎn)化為一個(gè)質(zhì)點(diǎn),在二維平面空間中運(yùn)動(dòng)[8]。航空器為了到達(dá)目標(biāo)點(diǎn)需要不斷地朝著目標(biāo)點(diǎn)運(yùn)動(dòng),目標(biāo)點(diǎn)對(duì)航空器產(chǎn)生了引力。同時(shí)航空器在飛行過(guò)程中需要避開受限空域和其他航空器,受限空域和其他各航空器對(duì)航空器產(chǎn)生了斥力,在斥力的作用下航空器始終與受限區(qū)和其他航空器保持一定的間隔,最終完成規(guī)避運(yùn)動(dòng)。航空器在運(yùn)動(dòng)空間中任意位置X=[x,y]T處的引力勢(shì)場(chǎng)函數(shù)為:

    (1)

    則吸引力是目標(biāo)引力勢(shì)函數(shù)的負(fù)梯度:

    Fatt(X) =-grad[UA(X)]=λ(XA-X)

    (2)

    式中,λ為引力位置增益系數(shù);X為航空器當(dāng)前位置;XA=[xA,yA]T為目標(biāo)位置矢量;‖XA-X‖為航空器與目標(biāo)點(diǎn)間的歐氏距離。

    將受限空域和其他航空器看作障礙物,其斥力勢(shì)函數(shù)定義為:

    Urei(X)=

    (3)

    相應(yīng)的斥力為:

    Frei(X) =-grad[Urei]

    (4)

    航空器在空域中所受合力為:

    (5)

    2改進(jìn)人工勢(shì)場(chǎng)法的多機(jī)改航路徑

    規(guī)劃

    傳統(tǒng)人工勢(shì)場(chǎng)法存在一些問(wèn)題:航空器所受的斥力和引力在一條直線上,航空器在改航期間出現(xiàn)合力為零,航空器無(wú)法避開障礙物,只能在一個(gè)范圍內(nèi)來(lái)回振動(dòng),即陷入局部最小值點(diǎn)。為解決該問(wèn)題,采用附加一個(gè)變化的控制力。在航空器之間避讓沖突的情形中,在目標(biāo)航空器周圍構(gòu)筑逆時(shí)針?lè)较虻男D(zhuǎn)力代替斥力,航空器在避讓危險(xiǎn)天氣的過(guò)程中通常是保持在危險(xiǎn)天氣的一定距離之外繞飛,在受限區(qū)域周圍構(gòu)造旋轉(zhuǎn)力場(chǎng)和斥力場(chǎng),航空器在旋轉(zhuǎn)力場(chǎng)和斥力場(chǎng)共同作用形成的復(fù)合人工勢(shì)場(chǎng)中尋找從初始位置到達(dá)目標(biāo)位置的無(wú)碰撞路徑。

    2.1 改進(jìn)人工勢(shì)場(chǎng)法

    本文對(duì)航空器的受力情況進(jìn)行調(diào)整,在目標(biāo)航空器周圍構(gòu)筑逆時(shí)針?lè)较虻男D(zhuǎn)力代替斥力。

    航空器i與j的避讓過(guò)程中,航空器i所受的旋轉(zhuǎn)力Frot為:

    (6)

    式中,krot為旋轉(zhuǎn)力增益系數(shù);rij=[(xi-xj)2+(yi-yj)2]1/2是航空器i和j之間的距離;rj為航空器j的保護(hù)區(qū)半徑;ρroj為航空器j的旋轉(zhuǎn)力勢(shì)場(chǎng)半徑;ρrot為航空器的旋轉(zhuǎn)力場(chǎng)影響的最大范圍。

    假設(shè)惡劣天氣的受限區(qū)為橢圓形,其斥力勢(shì)函數(shù)定義為:

    Ureo(qi,X)=

    (7)

    相應(yīng)的斥力為:

    Freo(qi,X)=

    (8)

    式中,kreo為受限區(qū)斥力的增益系數(shù);XR(x0,y0)為橢圓形危險(xiǎn)天氣區(qū)域的中心;ρ(qi,X)為航空器i到橢圓的距離;rreo為橢圓圓心與航空器i的連線上,圓心到橢圓邊界的長(zhǎng)度;ρreo為受限區(qū)斥力場(chǎng)影響的最大距離范圍。

    受限區(qū)的旋轉(zhuǎn)力為:

    (9)

    式中,ρro為受限區(qū)旋轉(zhuǎn)力場(chǎng)半徑;kro為受限區(qū)的旋轉(zhuǎn)力增益系數(shù)。這樣,在多航空器系統(tǒng)中,航空器i所受的合力為:

    Frot(qi,qj)+Fro(qi,X)]

    (j=1,…,m;i≠j)

    (10)

    以上建立了多航空器避讓危險(xiǎn)天氣時(shí)的模型,需要注意的是:在計(jì)算勢(shì)場(chǎng)力Fatt,Freo,Frot,Fro時(shí),增益系數(shù)katt,kreo,krot,kro的值決定了勢(shì)場(chǎng)力的大小,進(jìn)而決定了航空器運(yùn)動(dòng)的方向。

    2.2 改進(jìn)人工勢(shì)場(chǎng)法的改航路徑規(guī)劃步驟

    步驟1:為規(guī)劃選定模型參數(shù),包括航空器數(shù)量I、航空器i的起始點(diǎn)Oi(x0,y0)、終點(diǎn)Ti(xt,yt)、受限區(qū)中心位置Obs(x0,y0)、受限區(qū)的橢圓方程參數(shù)a,b以及ρreo,ρro,ρrot,katt,kreo,krot,kro,還有移動(dòng)步長(zhǎng)l和最大運(yùn)行步數(shù)J,按運(yùn)動(dòng)空間的信息建立人工勢(shì)場(chǎng)模型;

    步驟2:確定運(yùn)行步數(shù)j(初始步數(shù)為1),將第j步航空器i的當(dāng)前位置坐標(biāo)Xi(xj,yj)存入二維數(shù)組Currenti(xj,yj);

    步驟3:分析航空器、受限區(qū)以及目標(biāo)點(diǎn)的相對(duì)位置,分別計(jì)算每一架航空器所受的引力Fatt,斥力Freo,旋轉(zhuǎn)力Frot,Fro,并表示為分別沿x軸和y軸方向的分量;

    步驟4:計(jì)算每一架航空器所受的勢(shì)場(chǎng)力的合力F(x,y),航空器從點(diǎn)Xi(xj,yj)沿合力方向以固定步長(zhǎng)l移動(dòng)到下一航路點(diǎn)next(xj+1,yj+1),點(diǎn)的計(jì)算如下:

    next(xj+1,yj+1)=Xi(xj,yj)+lF

    (11)

    步驟5:移動(dòng)航空器到點(diǎn)next(xj+1,yj+1),置j+1點(diǎn)為當(dāng)前點(diǎn),j=j+1,xj=xj+1,yj=yj+1;

    步驟6:判斷航空器是否到達(dá)目標(biāo)區(qū)域,如果距離目標(biāo)區(qū)域的距離小于預(yù)設(shè)的閾值(通常為移動(dòng)步長(zhǎng)的整數(shù)倍),航空器到達(dá)目標(biāo)區(qū)域,結(jié)束改航路徑規(guī)劃,輸出數(shù)組Currenti從而生成完整的改航路徑,任務(wù)完成;否則轉(zhuǎn)步驟2繼續(xù)計(jì)算,或者直到步數(shù)j大于設(shè)定的最大步數(shù)J,無(wú)法找到完整路徑,說(shuō)明模型參數(shù)選擇不合適,需要改變參數(shù)重新規(guī)劃。

    3 仿真計(jì)算

    算例為8架航空器,四對(duì)同高度迎頭飛行,航空器的初始狀態(tài)為表1,航空器預(yù)計(jì)飛行航路和狀態(tài)如圖1所示。航空器在飛行中不僅要繞飛天氣,同時(shí)還必須避讓對(duì)向的航空器,且由于改航過(guò)程中需要占用側(cè)向間隔,相鄰兩架航空器的飛行沖突問(wèn)題也要加入考慮之中。

    表1 航空器初始狀態(tài)Table 1 Initial state of the aircraft

    圖1 航空器預(yù)計(jì)飛行航路和狀態(tài)Fig.1 Estimated flight path and state of the aircraft

    首先通過(guò)傳統(tǒng)人工勢(shì)場(chǎng)法作改航路徑的規(guī)劃,各參數(shù)的設(shè)置如下:勢(shì)場(chǎng)范圍ρreo=20 km,ρro=50 km,ρrea=20 km;勢(shì)場(chǎng)力增益系數(shù)katt=15,kreo=150,krea=170,kro=15。經(jīng)MATLAB計(jì)算得到的改航路徑如圖2所示。

    傳統(tǒng)人工勢(shì)場(chǎng)法可以規(guī)劃出符合間隔標(biāo)準(zhǔn)的改航路徑,其中航空器與危險(xiǎn)天氣的間隔保持在10 km以上;航空器之間的最小間隔為7.863 km,大于最低間隔標(biāo)準(zhǔn)的6 km。但是由于多航空器復(fù)雜沖突的存在,傳統(tǒng)人工勢(shì)場(chǎng)法生成的改航路徑非常曲折,包括多處180°以上的大轉(zhuǎn)彎;航空器飛行相當(dāng)長(zhǎng)一段距離后方能到達(dá)目標(biāo)地點(diǎn),特別是迎頭飛行的航空器遭遇之后經(jīng)過(guò)長(zhǎng)時(shí)間的振動(dòng)方得解脫;改航路徑在實(shí)際工作中是不可接受的。此外,整個(gè)仿真的迭代步數(shù)達(dá)到了300步以上,主要原因就是迎頭遭遇的航空器在振動(dòng)過(guò)程中消耗了大量步數(shù),造成計(jì)算資源的浪費(fèi),降低了仿真計(jì)算的效率。

    圖2 傳統(tǒng)人工勢(shì)場(chǎng)法規(guī)劃的改航路徑 Fig.2 Rerouting path of traditional artificial potential field

    采用前面提出的改進(jìn)人工勢(shì)場(chǎng)法解決上述問(wèn)題。由于航空器之間用旋轉(zhuǎn)力代替排斥力,所以將航空器斥力勢(shì)場(chǎng)ρrea變?yōu)樾D(zhuǎn)力勢(shì)場(chǎng)ρrot=20 km,航空器斥力增益系數(shù)krea變?yōu)樾D(zhuǎn)力增益系數(shù)krot=20,其他參數(shù)不變,通過(guò)仿真計(jì)算得到的改航路徑如圖3所示。

    圖3 改進(jìn)人工勢(shì)場(chǎng)法規(guī)劃的改航路徑Fig.3 Rerouting path planned by improved artificial potential field

    改進(jìn)人工勢(shì)場(chǎng)法規(guī)劃出的改航路徑較傳統(tǒng)方法規(guī)劃出的路徑更短,各航空器的航跡也更為平滑,沒(méi)有出現(xiàn)大角度轉(zhuǎn)彎的情況;迎頭相遇的各航空器都能很好地進(jìn)行避讓,沒(méi)有出現(xiàn)反復(fù)振蕩的情況,航空器之間的嚴(yán)重沖突得以避免。航空器與危險(xiǎn)天氣的距離均大于10 km,航空器之間最小間隔為6.729 km,符合最低間隔標(biāo)準(zhǔn)。此外,采用改進(jìn)人工勢(shì)場(chǎng)法的迭代步數(shù)均小于200,說(shuō)明對(duì)人工勢(shì)場(chǎng)法的改進(jìn)節(jié)約了計(jì)算時(shí)間,計(jì)算效率提高了33%,可見本文提出的改進(jìn)人工勢(shì)場(chǎng)法能夠克服傳統(tǒng)方法的一些不足之處,適用于多航空器改航的路徑規(guī)劃問(wèn)題。

    4 結(jié)束語(yǔ)

    本文采用改進(jìn)人工勢(shì)場(chǎng)法對(duì)多航空器改航路徑進(jìn)行了規(guī)劃,解決了傳統(tǒng)方法中航空器受力局部最小值的問(wèn)題,仿真了航空器迎頭沖突的情景。結(jié)果表明,改進(jìn)人工勢(shì)場(chǎng)法規(guī)劃的路徑更短、更平滑,沒(méi)有大角度轉(zhuǎn)彎和反復(fù)振蕩的情況。在此基礎(chǔ)上,利用空管運(yùn)行中的氣象和空域數(shù)據(jù)驗(yàn)證理論方法的實(shí)用性,是下一步研究的主要內(nèi)容。

    [1] Dixon M,Weiner G.Automated aircraft routing through weather-impacted airspace[C]//5th International Conference on Aviation Weather Systems.Vienna,VA,1993:295-298.

    [2] Sridhar B,Chatterji G B,Grabbe S,et al.Integration of traffic flow management decisions[C]//AIAA Guidance,Navigation,and Control Conference.Monterey,Cal,2002:1-9.

    [3] Krozel J,Penny S,Prete J,et al.Automated route generation for avoiding deterministic weather in transition airspace[J].Journal of Guidance,Control,and Dynamics,2007,30(1):144-153.

    [4] Balen C H,Bil C.Optimal re-routing of aircraft around closed airspace in free flight[C]//AIAA Modeling and Simulation Technologies Conference.Honolulu,HI,2008:1-10.

    [5] 宋柯.空中交通流量管理改航策略初步研究[D].南京:南京航空航天大學(xué),2002.

    [6] 李雄,徐肖豪,朱承元,等.基于幾何算法的空中交通改航路徑規(guī)劃研究[J].系統(tǒng)工程,2008,26(8):37-40.

    [7] 徐肖豪,李成功,趙嶷飛,等.基于人工勢(shì)場(chǎng)算法的改航路徑規(guī)劃[J].交通運(yùn)輸工程學(xué)報(bào),2009,9(6):64-68.

    [8] 況菲,王耀南,張輝.動(dòng)態(tài)環(huán)境下基于改進(jìn)人工勢(shì)場(chǎng)的機(jī)器人實(shí)時(shí)路徑規(guī)劃仿真研究[J].計(jì)算機(jī)應(yīng)用,2005,25(10):2415-2417.

    (編輯:方春玲)

    Multi-aircraftreroutingpathplanningbasedonimprovedartificialpotentialfieldalgorithm

    WANG Xing-long, XU Xiao-hao, FENG Jiang-ran

    (Research Base of Air Traffic Management, Civil Aviation University of China, Tianjin 300300, China)

    Severe weather is the main cause for the delay of transport aircraft, rerouting strategy is an effective way to reduce delays, multi-aircraft path planning is the key point for rerouting. For many flights delay caused by severe weather, this paper presents artificial potential field method for multi-aircraft rerouting path planning. Traditional artificial potential field model is established according to the aircraft by gravity force and repulsion force, and different rotation control force between aircraft and aircraft, aircraft and severe weather is introduced to avoid local minimum value during rerouting. Simulation result shows that the path planned by improved artificial potential field method is shorter and smoother, and no large angle turning and repeated shocks, and the computational efficiency is improved by 33%.

    air transportation; rerouting path planning; artificial potential field; multi-aircraft

    V355

    A

    1002-0853(2013)04-0381-04

    2013-01-15;

    2013-04-01; < class="emphasis_bold">網(wǎng)絡(luò)出版時(shí)間

    時(shí)間:2013-06-06 12:25

    國(guó)家自然科學(xué)基金資助(61039001,60972006);中央高??蒲谢鹳Y助(ZXH2011D010)

    王興隆(1979-),男,黑龍江北安人,助理研究員,研究方向?yàn)榭罩薪煌ü芾怼?/p>

    猜你喜歡
    勢(shì)場(chǎng)步數(shù)航空器
    速度和步數(shù),哪個(gè)更重要
    基于Frenet和改進(jìn)人工勢(shì)場(chǎng)的在軌規(guī)避路徑自主規(guī)劃
    基于改進(jìn)人工勢(shì)場(chǎng)方法的多無(wú)人機(jī)編隊(duì)避障算法
    楚國(guó)的探索之旅
    奇妙博物館(2021年4期)2021-05-04 08:59:48
    微信運(yùn)動(dòng)步數(shù)識(shí)人指南
    小演奏家(2018年9期)2018-12-06 08:42:02
    庫(kù)車坳陷南斜坡古流體勢(shì)場(chǎng)對(duì)陸相油氣運(yùn)聚的控制
    論航空器融資租賃出租人的違約取回權(quán)
    基于偶極勢(shì)場(chǎng)的自主水下航行器回塢導(dǎo)引算法
    航空器的順風(fēng)耳——機(jī)載衛(wèi)星通信
    太空探索(2016年3期)2016-07-12 09:58:45
    火星航空器何時(shí)才能首飛
    太空探索(2016年8期)2016-07-10 09:21:58
    亚洲人成电影观看| 午夜老司机福利片| 亚洲欧洲精品一区二区精品久久久| 亚洲国产最新在线播放| 亚洲av综合色区一区| 一区福利在线观看| 亚洲av欧美aⅴ国产| 菩萨蛮人人尽说江南好唐韦庄| 久久青草综合色| 一级黄片播放器| avwww免费| 男女之事视频高清在线观看 | 一级毛片女人18水好多 | 在线观看免费午夜福利视频| 日日爽夜夜爽网站| 久久精品国产a三级三级三级| 欧美人与善性xxx| 亚洲精品久久久久久婷婷小说| 国产激情久久老熟女| 99九九在线精品视频| 久久精品久久久久久久性| 久热这里只有精品99| 波野结衣二区三区在线| 桃花免费在线播放| 王馨瑶露胸无遮挡在线观看| 涩涩av久久男人的天堂| 国产免费福利视频在线观看| 欧美日韩亚洲国产一区二区在线观看 | 高清不卡的av网站| 在线亚洲精品国产二区图片欧美| 亚洲第一青青草原| 成人手机av| 国产人伦9x9x在线观看| 国产亚洲一区二区精品| 深夜精品福利| 欧美日韩福利视频一区二区| 国产成人av激情在线播放| 少妇 在线观看| 久久久精品免费免费高清| 国产成人精品久久二区二区91| 一本—道久久a久久精品蜜桃钙片| 亚洲国产欧美网| 激情视频va一区二区三区| 免费看不卡的av| 色精品久久人妻99蜜桃| 欧美久久黑人一区二区| 欧美日韩亚洲高清精品| 国产精品一区二区在线观看99| 久久精品久久久久久噜噜老黄| av国产精品久久久久影院| 亚洲精品国产区一区二| 日韩熟女老妇一区二区性免费视频| 超色免费av| 自拍欧美九色日韩亚洲蝌蚪91| 日日摸夜夜添夜夜爱| 一二三四在线观看免费中文在| 男女无遮挡免费网站观看| 最近最新中文字幕大全免费视频 | a级片在线免费高清观看视频| 一级,二级,三级黄色视频| 国产男女内射视频| 黄色视频在线播放观看不卡| 欧美黑人精品巨大| 久久性视频一级片| 精品高清国产在线一区| 少妇精品久久久久久久| 一级片免费观看大全| 国产激情久久老熟女| 精品国产国语对白av| 国产片特级美女逼逼视频| 国产又爽黄色视频| 婷婷色麻豆天堂久久| 丰满人妻熟妇乱又伦精品不卡| 水蜜桃什么品种好| 男女床上黄色一级片免费看| 国产成人av教育| 午夜日韩欧美国产| 亚洲少妇的诱惑av| 男的添女的下面高潮视频| 自拍欧美九色日韩亚洲蝌蚪91| 欧美日韩av久久| 国产午夜精品一二区理论片| 一区在线观看完整版| 一级毛片女人18水好多 | 大码成人一级视频| 99热全是精品| 日韩熟女老妇一区二区性免费视频| 精品一区二区三区四区五区乱码 | www.av在线官网国产| 在线观看免费午夜福利视频| 亚洲精品国产av成人精品| 免费黄频网站在线观看国产| 日韩av在线免费看完整版不卡| av一本久久久久| 亚洲熟女毛片儿| 中文字幕人妻丝袜一区二区| 色婷婷av一区二区三区视频| 亚洲欧美精品综合一区二区三区| 国产成人一区二区在线| 国产精品秋霞免费鲁丝片| 精品国产国语对白av| 国产一区有黄有色的免费视频| 大话2 男鬼变身卡| 青春草亚洲视频在线观看| 热re99久久精品国产66热6| 十八禁网站网址无遮挡| 亚洲免费av在线视频| 午夜日韩欧美国产| 一区二区日韩欧美中文字幕| 日本午夜av视频| 久久国产亚洲av麻豆专区| 97在线人人人人妻| 日韩av在线免费看完整版不卡| 久久热在线av| 三上悠亚av全集在线观看| 国产成人精品在线电影| 熟女少妇亚洲综合色aaa.| 亚洲欧美精品综合一区二区三区| 国产在线视频一区二区| 女人精品久久久久毛片| 亚洲欧美精品综合一区二区三区| 国产精品.久久久| 永久免费av网站大全| 久久久久国产一级毛片高清牌| 男男h啪啪无遮挡| 亚洲精品第二区| 国产欧美日韩精品亚洲av| 18禁黄网站禁片午夜丰满| 欧美乱码精品一区二区三区| 少妇被粗大的猛进出69影院| 日本一区二区免费在线视频| a级毛片黄视频| 国产精品国产av在线观看| 久久精品亚洲av国产电影网| 一区二区三区乱码不卡18| 国产xxxxx性猛交| 久久免费观看电影| 国产99久久九九免费精品| 嫩草影视91久久| 免费观看av网站的网址| 国产精品一区二区免费欧美 | 久久精品人人爽人人爽视色| 国产亚洲午夜精品一区二区久久| 国产欧美日韩一区二区三区在线| 天天操日日干夜夜撸| 久久影院123| 国产高清视频在线播放一区 | 各种免费的搞黄视频| 曰老女人黄片| 19禁男女啪啪无遮挡网站| 成人国产av品久久久| 国产精品欧美亚洲77777| 亚洲熟女毛片儿| 视频区欧美日本亚洲| 精品久久蜜臀av无| 美女视频免费永久观看网站| 满18在线观看网站| 脱女人内裤的视频| 久久精品aⅴ一区二区三区四区| 中文欧美无线码| 精品少妇内射三级| 国产精品国产av在线观看| 国产野战对白在线观看| 国产色视频综合| 午夜福利视频精品| 欧美在线黄色| 女警被强在线播放| 秋霞在线观看毛片| 亚洲人成77777在线视频| 高清av免费在线| 黄色视频在线播放观看不卡| 久久国产精品大桥未久av| 国产精品三级大全| 欧美 日韩 精品 国产| 最近中文字幕2019免费版| 国产亚洲午夜精品一区二区久久| 免费在线观看影片大全网站 | 中文字幕高清在线视频| 热99久久久久精品小说推荐| 18禁裸乳无遮挡动漫免费视频| 一级a爱视频在线免费观看| 另类亚洲欧美激情| 日韩大片免费观看网站| 免费在线观看影片大全网站 | 精品卡一卡二卡四卡免费| 久久狼人影院| 久久这里只有精品19| 一边摸一边抽搐一进一出视频| 亚洲人成网站在线观看播放| 国产成人欧美在线观看 | 99热国产这里只有精品6| 青春草亚洲视频在线观看| 飞空精品影院首页| 天堂中文最新版在线下载| 午夜老司机福利片| 啦啦啦啦在线视频资源| 天堂中文最新版在线下载| 黄片播放在线免费| 精品一品国产午夜福利视频| 观看av在线不卡| a级毛片黄视频| 在线天堂中文资源库| xxx大片免费视频| 一区二区三区激情视频| 亚洲国产精品一区二区三区在线| 黑人欧美特级aaaaaa片| 新久久久久国产一级毛片| 欧美日韩成人在线一区二区| 中文乱码字字幕精品一区二区三区| 高清视频免费观看一区二区| 亚洲欧洲日产国产| 亚洲av成人不卡在线观看播放网 | 国产高清videossex| 亚洲av成人不卡在线观看播放网 | 制服人妻中文乱码| 国产免费又黄又爽又色| 一区二区av电影网| 另类精品久久| 天天躁夜夜躁狠狠久久av| 在线观看免费高清a一片| 成人午夜精彩视频在线观看| 日韩大片免费观看网站| 国产一级毛片在线| 亚洲一码二码三码区别大吗| 久久精品国产综合久久久| 久久久欧美国产精品| 菩萨蛮人人尽说江南好唐韦庄| 一边摸一边做爽爽视频免费| 亚洲一码二码三码区别大吗| 人妻人人澡人人爽人人| 女人久久www免费人成看片| 国产激情久久老熟女| 日韩制服丝袜自拍偷拍| 成年人黄色毛片网站| 欧美 亚洲 国产 日韩一| 在线精品无人区一区二区三| 国产亚洲欧美精品永久| 免费在线观看影片大全网站 | 婷婷色综合www| 在线精品无人区一区二区三| 国产人伦9x9x在线观看| 亚洲精品久久成人aⅴ小说| 亚洲欧美成人综合另类久久久| 久久久久国产精品人妻一区二区| 久久人妻熟女aⅴ| 久久精品久久精品一区二区三区| 国产伦理片在线播放av一区| 男女高潮啪啪啪动态图| 免费少妇av软件| 人人妻人人澡人人看| 亚洲熟女毛片儿| 亚洲欧美精品自产自拍| av国产精品久久久久影院| 两个人免费观看高清视频| e午夜精品久久久久久久| 建设人人有责人人尽责人人享有的| 精品国产乱码久久久久久男人| 国产高清国产精品国产三级| 久久国产精品男人的天堂亚洲| 婷婷色综合大香蕉| 亚洲欧美精品综合一区二区三区| 国产成人精品无人区| 午夜日韩欧美国产| 女人久久www免费人成看片| av天堂久久9| 少妇精品久久久久久久| 操美女的视频在线观看| 欧美+亚洲+日韩+国产| 免费人妻精品一区二区三区视频| 欧美日韩av久久| 最近中文字幕2019免费版| 久久99热这里只频精品6学生| 亚洲成人免费av在线播放| 美女主播在线视频| 国产成人欧美| 久久精品亚洲av国产电影网| 一边摸一边做爽爽视频免费| 亚洲国产欧美一区二区综合| 国产精品久久久久久精品电影小说| 成在线人永久免费视频| 大片免费播放器 马上看| 中文字幕人妻丝袜一区二区| 9色porny在线观看| 国产伦理片在线播放av一区| 久久午夜综合久久蜜桃| 赤兔流量卡办理| 少妇的丰满在线观看| 国产欧美日韩一区二区三区在线| 观看av在线不卡| 午夜两性在线视频| 久久99一区二区三区| 王馨瑶露胸无遮挡在线观看| 国产三级黄色录像| 日日爽夜夜爽网站| 欧美激情极品国产一区二区三区| 国产99久久九九免费精品| 亚洲视频免费观看视频| 国产片内射在线| 欧美在线黄色| 色综合欧美亚洲国产小说| av在线app专区| 又大又黄又爽视频免费| 另类精品久久| 国产一区二区三区av在线| 国产精品免费视频内射| 久久免费观看电影| 久热爱精品视频在线9| 精品免费久久久久久久清纯 | 国产精品久久久久久人妻精品电影 | 又大又黄又爽视频免费| 亚洲精品国产av蜜桃| 国产精品一国产av| 啦啦啦啦在线视频资源| 日本av免费视频播放| 一二三四在线观看免费中文在| 国产成人精品久久久久久| 丰满少妇做爰视频| 国产熟女欧美一区二区| 又黄又粗又硬又大视频| 天堂中文最新版在线下载| 日韩大码丰满熟妇| 97人妻天天添夜夜摸| 亚洲七黄色美女视频| 一本久久精品| 丝袜美腿诱惑在线| 好男人视频免费观看在线| 捣出白浆h1v1| 啦啦啦中文免费视频观看日本| 黄片小视频在线播放| 精品一区在线观看国产| 一区二区日韩欧美中文字幕| 黑人猛操日本美女一级片| 亚洲第一av免费看| 亚洲国产欧美在线一区| 精品国产一区二区三区四区第35| xxx大片免费视频| 日本91视频免费播放| 午夜老司机福利片| 十八禁高潮呻吟视频| 欧美日韩一级在线毛片| 亚洲 国产 在线| 日本午夜av视频| 亚洲精品国产区一区二| 99香蕉大伊视频| 欧美亚洲日本最大视频资源| 两性夫妻黄色片| 日韩中文字幕欧美一区二区 | 伊人久久大香线蕉亚洲五| 后天国语完整版免费观看| 欧美精品亚洲一区二区| 国产女主播在线喷水免费视频网站| 一本大道久久a久久精品| 大陆偷拍与自拍| av国产久精品久网站免费入址| 美女高潮到喷水免费观看| 高潮久久久久久久久久久不卡| 欧美老熟妇乱子伦牲交| 亚洲欧洲精品一区二区精品久久久| 人妻 亚洲 视频| 欧美中文综合在线视频| 久久精品熟女亚洲av麻豆精品| av又黄又爽大尺度在线免费看| 国产精品麻豆人妻色哟哟久久| 两人在一起打扑克的视频| 老熟女久久久| 女人精品久久久久毛片| 久久国产精品男人的天堂亚洲| 久久久久国产一级毛片高清牌| 国产在线观看jvid| 熟女少妇亚洲综合色aaa.| av欧美777| 国产精品免费大片| 午夜免费男女啪啪视频观看| 久久国产精品人妻蜜桃| 亚洲国产日韩一区二区| 真人做人爱边吃奶动态| 国产精品av久久久久免费| av国产久精品久网站免费入址| 成人亚洲欧美一区二区av| 好男人视频免费观看在线| 最近中文字幕2019免费版| 日韩免费高清中文字幕av| 国产伦人伦偷精品视频| a 毛片基地| 日韩av在线免费看完整版不卡| 一级黄片播放器| 一区二区三区精品91| 妹子高潮喷水视频| 中文字幕人妻熟女乱码| 丝瓜视频免费看黄片| 国产成人av教育| 人人妻,人人澡人人爽秒播 | 手机成人av网站| 国产精品一区二区在线观看99| 夫妻午夜视频| 自线自在国产av| 美女主播在线视频| 国产精品99久久99久久久不卡| 国产精品av久久久久免费| 在线观看国产h片| tube8黄色片| 免费人妻精品一区二区三区视频| 免费在线观看影片大全网站 | 男人爽女人下面视频在线观看| 天天躁夜夜躁狠狠久久av| 亚洲五月色婷婷综合| 久久精品亚洲av国产电影网| 纯流量卡能插随身wifi吗| 少妇人妻 视频| 国产野战对白在线观看| 2018国产大陆天天弄谢| 亚洲视频免费观看视频| av在线播放精品| 亚洲成色77777| 亚洲欧美一区二区三区久久| 大型av网站在线播放| 一二三四在线观看免费中文在| 超色免费av| 蜜桃国产av成人99| 少妇人妻久久综合中文| 国产又爽黄色视频| 中文字幕最新亚洲高清| 欧美激情极品国产一区二区三区| 一本综合久久免费| 亚洲黑人精品在线| 亚洲中文字幕日韩| 一区在线观看完整版| 嫩草影视91久久| 一本—道久久a久久精品蜜桃钙片| 欧美97在线视频| 国产精品 欧美亚洲| a级片在线免费高清观看视频| 精品少妇久久久久久888优播| 午夜福利免费观看在线| 五月天丁香电影| 日韩av不卡免费在线播放| 日韩中文字幕欧美一区二区 | 国产免费又黄又爽又色| 久久影院123| 看免费成人av毛片| 亚洲精品国产av成人精品| 美女福利国产在线| 免费高清在线观看日韩| 国产片特级美女逼逼视频| 伊人亚洲综合成人网| 狂野欧美激情性bbbbbb| 制服诱惑二区| 亚洲精品国产区一区二| 麻豆av在线久日| 国产精品99久久99久久久不卡| 国产精品成人在线| 一边摸一边做爽爽视频免费| 亚洲一区二区三区欧美精品| 久久精品久久久久久久性| 国产欧美日韩综合在线一区二区| 在线观看一区二区三区激情| 国产精品一区二区在线不卡| 国产一区二区激情短视频 | 亚洲第一av免费看| 久久国产精品人妻蜜桃| 波野结衣二区三区在线| 下体分泌物呈黄色| 国产熟女午夜一区二区三区| 99久久人妻综合| 国产亚洲欧美精品永久| 亚洲成av片中文字幕在线观看| 国产成人精品久久二区二区免费| 国产女主播在线喷水免费视频网站| 亚洲熟女毛片儿| bbb黄色大片| 日本欧美视频一区| 久久影院123| 亚洲av日韩精品久久久久久密 | 女性生殖器流出的白浆| 我要看黄色一级片免费的| 色精品久久人妻99蜜桃| 亚洲精品国产av蜜桃| 国产精品人妻久久久影院| 午夜免费男女啪啪视频观看| 午夜福利视频精品| 亚洲av欧美aⅴ国产| 蜜桃国产av成人99| 少妇的丰满在线观看| 国产不卡av网站在线观看| 国产成人免费无遮挡视频| 中文字幕人妻丝袜制服| 亚洲成人免费电影在线观看 | 在线观看免费视频网站a站| 50天的宝宝边吃奶边哭怎么回事| 亚洲九九香蕉| 少妇 在线观看| 久久精品久久久久久噜噜老黄| 亚洲第一av免费看| 我的亚洲天堂| 国产亚洲精品久久久久5区| 秋霞在线观看毛片| 日韩制服骚丝袜av| 成人午夜精彩视频在线观看| 国产1区2区3区精品| 美女脱内裤让男人舔精品视频| 黄色视频在线播放观看不卡| 免费在线观看黄色视频的| 欧美精品高潮呻吟av久久| 搡老乐熟女国产| 日韩熟女老妇一区二区性免费视频| 老司机影院成人| 国产精品一国产av| 国产日韩欧美在线精品| 久久热在线av| 一级毛片黄色毛片免费观看视频| 日韩中文字幕欧美一区二区 | 考比视频在线观看| 99久久综合免费| 国产亚洲精品第一综合不卡| 青春草亚洲视频在线观看| 色网站视频免费| 亚洲欧美日韩另类电影网站| 少妇人妻久久综合中文| 国产精品国产av在线观看| 99久久综合免费| 一级,二级,三级黄色视频| 成人国产一区最新在线观看 | 欧美变态另类bdsm刘玥| 热re99久久国产66热| 日韩人妻精品一区2区三区| 999久久久国产精品视频| 一边亲一边摸免费视频| 国产三级黄色录像| 国产精品99久久99久久久不卡| 日本一区二区免费在线视频| 日韩精品免费视频一区二区三区| 免费一级毛片在线播放高清视频 | 精品国产乱码久久久久久男人| 亚洲欧美激情在线| 日本vs欧美在线观看视频| 亚洲视频免费观看视频| 一区二区三区乱码不卡18| 只有这里有精品99| 午夜福利视频在线观看免费| 9色porny在线观看| 国产成人精品久久久久久| 悠悠久久av| 成人三级做爰电影| 大香蕉久久网| 97精品久久久久久久久久精品| 午夜老司机福利片| 久久天堂一区二区三区四区| 人妻人人澡人人爽人人| 亚洲美女黄色视频免费看| 国产一区有黄有色的免费视频| 亚洲欧美日韩高清在线视频 | 丝袜美腿诱惑在线| 啦啦啦啦在线视频资源| 国产男人的电影天堂91| 老汉色av国产亚洲站长工具| 国产老妇伦熟女老妇高清| 99香蕉大伊视频| 一区福利在线观看| 亚洲国产欧美网| 免费高清在线观看视频在线观看| 在线观看人妻少妇| 日本av手机在线免费观看| 亚洲国产精品一区二区三区在线| 午夜精品国产一区二区电影| 男女高潮啪啪啪动态图| 亚洲九九香蕉| 国产成人欧美| √禁漫天堂资源中文www| 久久久国产精品麻豆| 一级a爱视频在线免费观看| 少妇人妻久久综合中文| 亚洲av美国av| 不卡av一区二区三区| 婷婷色综合大香蕉| 亚洲精品久久成人aⅴ小说| 9色porny在线观看| 国产在线免费精品| 亚洲av片天天在线观看| 99国产精品一区二区蜜桃av | tube8黄色片| 最新在线观看一区二区三区 | 成年美女黄网站色视频大全免费| 亚洲午夜精品一区,二区,三区| 欧美精品人与动牲交sv欧美| 午夜两性在线视频| 亚洲一卡2卡3卡4卡5卡精品中文| 黑人巨大精品欧美一区二区蜜桃| 亚洲伊人久久精品综合| 国产亚洲精品第一综合不卡| 欧美日韩亚洲国产一区二区在线观看 | 午夜激情av网站| 国产色视频综合| 纯流量卡能插随身wifi吗| 一级毛片 在线播放| 青春草视频在线免费观看| 国产日韩一区二区三区精品不卡| 国产精品.久久久| 亚洲熟女毛片儿| 色综合欧美亚洲国产小说| 午夜福利免费观看在线| 丰满饥渴人妻一区二区三| 狂野欧美激情性xxxx| 国产精品.久久久| 建设人人有责人人尽责人人享有的| 91麻豆av在线| 国产免费现黄频在线看| 亚洲欧美激情在线| 中文精品一卡2卡3卡4更新| 国产精品.久久久| 国产亚洲一区二区精品| 国产男女内射视频| 人体艺术视频欧美日本| 少妇粗大呻吟视频| 搡老乐熟女国产| 纵有疾风起免费观看全集完整版|