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

    通用航空低空避撞研究

    2018-11-22 11:17:56郭文豪黃俊賢張書煜張潤(rùn)澤江波
    科技創(chuàng)新與應(yīng)用 2018年31期
    關(guān)鍵詞:通用航空

    郭文豪 黃俊賢 張書煜 張潤(rùn)澤 江波

    摘 要:為了避免通航飛機(jī)在低空飛行時(shí),與建筑物、地形等發(fā)生碰撞,根據(jù)不同障礙物的特點(diǎn)和通用航空器的性能,利用矢量分解的方法,將通航飛機(jī)在三維空間中的避撞分解為飛機(jī)投影在水平方向和垂直方向上的避撞。算法利用障礙物的信息,判斷碰撞的風(fēng)險(xiǎn),從而給出高度調(diào)整或航向調(diào)整的避撞方法。運(yùn)用該算法通過Matlab進(jìn)行障礙物避撞的仿真實(shí)驗(yàn),結(jié)果表明本算法能有效避免通航飛機(jī)與障礙物的相撞。

    關(guān)鍵詞:低空開放;通用航空;幾何算法;碰撞風(fēng)險(xiǎn);避撞路徑

    中圖分類號(hào):V328 文獻(xiàn)標(biāo)志碼:A 文章編號(hào):2095-2945(2018)31-0001-05

    Abatract: In order to avoid the collision of navigable aircraft with buildings and terrain while flying at low altitude, according to the characteristics of different obstacles and the performance of general aircraft, the method of vector decomposition is used. In this paper, the collision avoidance of navigable aircraft in three-dimensional space is divided into the collision avoidance of aircraft projected in the horizontal direction and the vertical direction. The algorithm uses the information of obstacles to judge the risk of collision, so as to give a collision avoidance method for height adjustment or course adjustment. The algorithm is applied to the simulation of obstacle avoidance by Matlab. The results show that the algorithm can effectively avoid the collision between navigable aircraft and obstacles.

    Keymords: low altitude open; general aviation; geometric algorithm; collision risk; collision avoidance path

    1 概述

    近年來,中國經(jīng)濟(jì)的快速發(fā)展,帶來通用航空業(yè)務(wù)的巨大需求,尤其是公務(wù)航空、私人航空、緊急救護(hù)、工業(yè)航空和農(nóng)林航空存在巨大增長(zhǎng)潛力[1]。預(yù)測(cè)顯示到2020年中國通用航空飛機(jī)數(shù)量可達(dá)12610架,產(chǎn)業(yè)增加值將達(dá)50534百萬元[2]。2017年2月提出的《通用航空“十三五”發(fā)展規(guī)劃》將進(jìn)一步擴(kuò)大低空空域開放[3],預(yù)示著未來幾年我國將迎來低空開放的新局面。但目前我國通用航空尚處于初步發(fā)展階段,事故萬時(shí)率遠(yuǎn)遠(yuǎn)高于運(yùn)輸航空。在2001年至2010年間,運(yùn)輸航發(fā)生了三起航空事故,而通用航空發(fā)生了22起飛行事故。其中直接原因?yàn)闄C(jī)組原因的事故占事故總數(shù)的72.7%[4]。這嚴(yán)重限制了通航產(chǎn)業(yè)的健康發(fā)展。為提高通航飛行安全系數(shù),減少事故悲劇,急需建立一套可靠有效的低空飛行對(duì)地避撞算法。而讓航空器避免與高大建筑物或地形等障礙物發(fā)生碰撞,是我們著重要解決的問題。

    鑒于此,本項(xiàng)目將對(duì)低空開放后通用航空器對(duì)地面障礙物在三維空間內(nèi)的幾何避撞算法展開研究,綜合考慮航空器性能和障礙物的不同特征,利用矢量分解,將三維空間的避撞分解為在水平方向和垂直方向的避撞,并通過算例分析驗(yàn)證該算法的作用。

    2 幾何碰撞算法

    因?yàn)轱w行器尺寸遠(yuǎn)小于安全間隔,所以飛行器任何時(shí)候均可等效為一個(gè)點(diǎn)。本算法中將飛機(jī)保護(hù)區(qū)模型視為三維空間內(nèi)圓心為(x0,y0,H0),半徑為r米的球體[5]。通過將飛機(jī)模型投影到水平和垂直方向上,再把水平和垂直兩個(gè)方向的結(jié)果相結(jié)合,使復(fù)雜的三維防撞問題簡(jiǎn)化為兩個(gè)二維防撞問題[6]。

    該避撞算法會(huì)在飛機(jī)按照當(dāng)前速度到達(dá)障礙物前60s(參照增強(qiáng)近地告警系統(tǒng)的設(shè)置),根據(jù)飛機(jī)實(shí)時(shí)的速度、高度、加速度判斷是否存在避撞風(fēng)險(xiǎn),若存在則提供避撞路徑;若不存在,則可保持現(xiàn)有飛行狀態(tài)繼續(xù)飛行。

    2.1 對(duì)高大建筑物避撞

    2.1.1 垂直方向避撞

    建筑物在垂直面上的投影為若干定位點(diǎn)組成的平面圖形,這些定位點(diǎn)的坐標(biāo)為(xn,Hn),飛機(jī)垂直投影實(shí)時(shí)位置坐標(biāo)為小圓圓心坐標(biāo)(x0,H0)。

    首先,判斷飛行高度是否不低于建筑物高度之上r米:若不低于,則無碰撞風(fēng)險(xiǎn);若低于,則存在碰撞風(fēng)險(xiǎn),進(jìn)行路徑計(jì)算。

    避撞路徑的計(jì)算:選取障礙物最高點(diǎn)中與航空器直線距離最短的點(diǎn)(Hn的值最大的前提下,xn-x0值最?。┯?jì)算避撞路徑。設(shè)飛機(jī)水平速度為Vx,水平位移為X,爬升加速度為a,爬升時(shí)間設(shè)定為t(60s)。根據(jù)拋物線運(yùn)動(dòng)的分解

    X=Vxt (1)

    H=at2 (2)

    完成加速上升階段飛機(jī)坐標(biāo)為(X'=X0+Vxt,H'=H0+at2)

    上升軌跡為

    H=X2 (3)

    (X',H')處的切線斜率為

    k=X' (4)

    根據(jù)點(diǎn)到直線的距離方程

    d= (5)

    當(dāng)定位點(diǎn)與切線的距離等于r即d=r時(shí):

    r= (6)

    C=H0+at2-kX'(7)

    則飛機(jī)剛好避開障礙物。再由

    x= (8)

    得到經(jīng)過爬升時(shí)間t剛好避開障礙物的爬升加速度amin的解。又因?yàn)樗帽茏猜窂讲荒艽┻^建筑物投影,可排除不符合要求的解,進(jìn)而得到避撞所需最小上爬升率Vy=at。

    若amin大于amax(飛機(jī)性能允許的最大上升加速度),則飛機(jī)無法利用爬升避開障礙物,需要檢索飛機(jī)前方120°范圍內(nèi)的其他高點(diǎn),是否存在飛機(jī)采用不超過最大可用爬升加速度能夠飛越的高點(diǎn),若存在這樣的高點(diǎn),則選擇滿足要求的高點(diǎn)中高度最低的高點(diǎn)(x,y,H),進(jìn)行垂直避撞;若不存在,則對(duì)當(dāng)前障礙物采取水平避撞。

    2.1.2 水平方向避撞

    建筑物在水平面上的投影為若干定位點(diǎn)組成的平面圖形,定位點(diǎn)的坐標(biāo)為(xn,yn),飛機(jī)水平投影實(shí)時(shí)位置坐標(biāo)為小圓圓心坐標(biāo)(x0,y0)。

    首先通過實(shí)時(shí)速度(由慣導(dǎo)系統(tǒng)得到)得出通航飛機(jī)的預(yù)計(jì)飛行軌跡(直線),計(jì)算預(yù)計(jì)軌跡是否經(jīng)過建筑物的投影,若不經(jīng)過,則無碰撞危險(xiǎn);若經(jīng)過則發(fā)出碰撞警告,提供避撞決策。

    避撞路徑的計(jì)算:計(jì)算所有定位點(diǎn)與圓心連線的斜率。

    K= (9)

    選出斜率最小值Kmin和最大值Kmax對(duì)應(yīng)的定位點(diǎn)計(jì)算避撞路徑。飛機(jī)的預(yù)計(jì)路徑為kx-y+(y0-kx0)=0,根據(jù)點(diǎn)到直線距離的方程,當(dāng)d=r時(shí),解得k的兩個(gè)解表示為k1=m,k2=n,對(duì)應(yīng)直線為L(zhǎng)1:m(x-x0)+y0=y;直線L2:n(x-x0)+y0=y。通過驗(yàn)算L1,L2不通過建筑物投影,去掉不符合要求的直線,從而得到避撞路徑。最后選取轉(zhuǎn)彎角度較小的作為避撞路徑,若兩轉(zhuǎn)彎角度相等,則都可作為避撞路徑。

    a?燮a?燮aVmca

    amax:飛機(jī)最大可用爬升加速度;Vmca:空中最小操縱速度;Vmax:飛機(jī)最大可用速度;ɑ:飛機(jī)轉(zhuǎn)彎角度。

    2.2 對(duì)地形的避撞

    2.2.1 垂直方向避撞

    地形障礙在垂直面上的投影為若干定位點(diǎn)組成的平面圖形,這些定位點(diǎn)的坐標(biāo)設(shè)為(xn,Hn),飛機(jī)垂直投影實(shí)時(shí)位置坐標(biāo)為小圓圓心坐標(biāo)(x0,H0)。

    首先判斷飛行高度是否高于高點(diǎn)高度之上r米:若不低于,則無碰撞風(fēng)險(xiǎn);若低于,則存在碰撞風(fēng)險(xiǎn),進(jìn)行路徑計(jì)算。

    避撞路徑的計(jì)算:設(shè)飛機(jī)水平速度為Vx,爬升加速度為a(a>0),爬升時(shí)間設(shè)定為t(60s)。同樣根據(jù)點(diǎn)到直線的距離方程將高點(diǎn)代入計(jì)算,這里分d≥r和d

    當(dāng)d≥r時(shí),需驗(yàn)證加速爬升階段與地形障礙的最小距離是否小于r:對(duì)于投影上的一個(gè)定位點(diǎn),拋物線上都存在與之對(duì)應(yīng)的一點(diǎn)使兩點(diǎn)間的距離最短,這個(gè)距離即定位點(diǎn)到拋物線的最短距離,且有過該兩點(diǎn)的直線與拋物線上過對(duì)應(yīng)點(diǎn)的切線相垂直。得到拋物線切線方程:

    x3+(Vx2-aHn)x-xnVx2=0 (10)

    可利用盛金公式法[7]對(duì)拋物線切線方程求解得到拋物線上與定位點(diǎn)對(duì)應(yīng)的點(diǎn)的坐標(biāo)(xs,Hs)(盛金公式:一元三次方程aX3+bX2+cX+d=0,(a,b,c,d∈R,且a≠0),重根判別式

    A=b2-3acB=bc-9adC=c2-3bd (11)

    總判別式

    ?駐=B2-4ac (12)

    當(dāng)A=B=0時(shí)

    x=x2=x3===(13)

    當(dāng)?駐=B2-4AC>0時(shí)

    x1=(14)

    x2,x3=(15)

    其中

    Y1,2=Ab+3a

    i2=-1

    當(dāng)?駐=B2-4AC=0時(shí)

    x1=+K (17)

    x2=x3= (18)

    其中K=,(A≠0);

    當(dāng)?駐=B2-4AC<0時(shí),

    x1=(19)

    x2,x3= (20)

    ?茲=arccosT,T=(A>0,-1

    求出拋物線上與定位點(diǎn)相對(duì)應(yīng)的點(diǎn)的坐標(biāo)后,分別計(jì)算定位點(diǎn)到拋物線的最短距離即兩點(diǎn)間的距離,取其中的最小值記為L(zhǎng)min,當(dāng)Lmin≥r時(shí),保持當(dāng)前爬升加速度飛機(jī)無碰撞風(fēng)險(xiǎn),避撞所需爬升率Vy=at;當(dāng)Lmin

    = (21)

    r= (22)

    解出(xb,Hb),可知

    T=(23)

    Hb=a′T2+H0 (24)

    得出a′,即可得到Vy=a′T

    若a′>amax(飛機(jī)性能的最大可用爬升加速度),則需要在采用amax的前提下減小Vx,有

    H2=amax+H0 (25)

    求出Vx′(Vx′≥空中最小操縱速度Vmca)

    Vy=amaxt (26)

    若Vx′Vmca)能夠飛越的高點(diǎn),若存在這樣的高點(diǎn),則選擇滿足要求的高度最低的高點(diǎn)進(jìn)行垂直避撞;若不存在,則采取水平避撞。

    當(dāng)d

    amin?燮a?燮amaxVmca

    amax:飛機(jī)最大可用爬升加速度;Vmca:空中最小操縱速度;Vmax:飛機(jī)最大可用速度;ɑ:飛機(jī)轉(zhuǎn)彎角度。

    2.2.2 水平方向避撞

    高點(diǎn)的坐標(biāo)為(xn,Hn),飛機(jī)實(shí)時(shí)位置坐標(biāo)為小圓圓心坐標(biāo)(x0,H0)。

    通過實(shí)時(shí)速度矢量(由慣導(dǎo)系統(tǒng)得到)得出通航飛機(jī)的預(yù)計(jì)飛行軌跡(直線),計(jì)算預(yù)計(jì)軌跡是否在安全余度內(nèi)經(jīng)過高點(diǎn),若不經(jīng)過,則無碰撞危險(xiǎn);若經(jīng)過則發(fā)出碰撞警告,提供避撞決策。

    避撞路徑的計(jì)算:計(jì)算飛機(jī)前方180°扇形區(qū)域內(nèi)所有定位點(diǎn)與圓心連線的斜率。選出斜率最小值Kmin和最大值Kmax對(duì)應(yīng)的定位點(diǎn)計(jì)算避撞路徑。根據(jù)點(diǎn)到直線的距離為r時(shí)求出直線L1:

    mx-y+Cm=0(Cm=y0-mx0)

    直線L2:

    nx-y+Cn=0(Cn=y0-nx0)

    對(duì)最小值Kmin定位點(diǎn)的避撞路徑選擇m和n中較小的;對(duì)最大值Kmax定位點(diǎn)的避撞路徑選擇m和n中斜率較大的。兩避撞路徑中取轉(zhuǎn)彎角度較小的作為避撞路徑,若兩轉(zhuǎn)彎角度相等,則都可作為避撞路徑。

    3 算法仿真

    3.1 避撞仿真流程及其數(shù)據(jù)設(shè)置

    通過獲取飛機(jī)飛行前方的障礙物高度和位置信息,可判斷飛機(jī)與最高障礙物之間的相對(duì)運(yùn)動(dòng)關(guān)系和位置關(guān)系。如飛機(jī)此時(shí)存在與障礙物相撞的危險(xiǎn),則通過以上避撞算法進(jìn)行避撞,其中包括垂直方向上的避撞和水平方向上的避撞。本節(jié)應(yīng)用Matlab進(jìn)行算法的仿真實(shí)驗(yàn),以驗(yàn)證算法的正確性。仿真流程圖如圖1。參考文獻(xiàn)[8]的設(shè)置,本文將飛機(jī)保護(hù)區(qū)半徑設(shè)為161m,設(shè)飛機(jī)以50m/s的水平速度勻速飛行。其余參數(shù)由Matlab隨機(jī)生成。

    3.2 高大建筑物避撞仿真

    3.2.1 垂直避撞仿真

    對(duì)于規(guī)則建筑物的垂直避撞仿真實(shí)驗(yàn),現(xiàn)生成垂直面下的位置數(shù)據(jù):飛機(jī)(0,0),障礙物點(diǎn)1(1450,75),障礙物點(diǎn)2(1510,75)。根據(jù)以上數(shù)據(jù),由垂直避撞算法可得,當(dāng)飛機(jī)以1.02m/s2的加速度上升至43m時(shí),飛機(jī)將不會(huì)與障礙物相撞。如圖2和圖3所示,飛機(jī)進(jìn)行調(diào)整后,能保持與障礙物之間的最小距離等于161m,即成功避免了相撞問題。

    3.2.2 水平避撞仿真

    對(duì)于規(guī)則建筑物的水平避撞仿真實(shí)驗(yàn),現(xiàn)生成水平面下的位置數(shù)據(jù):飛機(jī)(0,0),障礙物點(diǎn)1(1000,100),障礙物點(diǎn)2(1000,200),障礙物點(diǎn)3(1060,100),障礙物點(diǎn)4(1060,200)。根據(jù)以上數(shù)據(jù),由水平避撞算法可得,當(dāng)飛機(jī)左偏至20.4°或右偏至3.5°時(shí),飛機(jī)在未來將不會(huì)與障礙物相撞。如圖4和圖5所示,飛機(jī)進(jìn)行調(diào)整后,能保持與障礙物之間的最小距離大于等于161m,即成功避免了相撞問題。

    3.3 地形避撞仿真

    3.3.1 垂直避撞仿真

    對(duì)于規(guī)則地形的垂直避撞仿真實(shí)驗(yàn),現(xiàn)生成垂直面下的位置數(shù)據(jù):飛機(jī)(-7,243),障礙物點(diǎn)1(1425,325),障礙物點(diǎn)2(1660,273),障礙物點(diǎn)3(1310,281),障礙物點(diǎn)4(1612,263),障礙物點(diǎn)5(1184,243),障礙物點(diǎn)6(1124,222),障礙物點(diǎn)7(1053,188)。求得當(dāng)飛機(jī)以1.25m/s2的加速度上升至40m并維持此狀態(tài)爬升時(shí),可安全飛越地形。飛機(jī)上升軌跡如圖6,規(guī)避期間與最高點(diǎn)的距離變化如圖7。

    3.3.2 水平避撞仿真

    對(duì)于不規(guī)則地形物的水平避撞仿真實(shí)驗(yàn),現(xiàn)生成了如下位置數(shù)據(jù):飛機(jī)(0,0),障礙物點(diǎn)1(1000,100),障礙物點(diǎn)2(900,200),障礙物點(diǎn)3(1200,80),障礙物點(diǎn)4(1100,300),障礙物點(diǎn)5(1300,150)。根據(jù)以上數(shù)據(jù),由水平避撞算法可得,當(dāng)飛機(jī)左偏至23.4°或右偏至3.9°時(shí),飛機(jī)在未來將不會(huì)與障礙物相撞。如圖8和圖9所示,飛機(jī)進(jìn)行調(diào)整后,就能保持與障礙物之間的最小距離大于等于161m,即成功避免了相撞問題。

    3.4 避撞仿真結(jié)果

    現(xiàn)通過Matlab隨機(jī)產(chǎn)生飛機(jī)和障礙物點(diǎn)的位置,并進(jìn)行100000次的仿真實(shí)驗(yàn)驗(yàn)結(jié)果如表1所示??梢姳舅惴軌蛴行П苊怙w機(jī)與障礙物的相撞,且避撞成功率達(dá)到了100%。

    4 結(jié)束語

    本文提出了基于幾何分析的在三維空間中通航低空避撞算法,其中包括垂直避撞和水平避撞兩部分,旨在為通航飛行避撞提供有效的避撞策略。為使算法更具有針對(duì)性,本文將地面障礙物劃分為規(guī)則的建筑物和不規(guī)則的地形障礙物,并分別進(jìn)行了避撞算法的研究。最后將各部分的算法進(jìn)行了仿真實(shí)驗(yàn),實(shí)驗(yàn)表明,算法可給出正確的避撞方法,使飛機(jī)避開前方的地面障礙物。為使算法更加貼近實(shí)際情況,未來將對(duì)算法的精密度和一些特殊情況如高壓線等做進(jìn)一步的研究。

    參考文獻(xiàn):

    [1]于光妍,杜青蕓.從產(chǎn)業(yè)角度看通用航空布局模式[J].中國民用航空,2013(11):44-45.

    [2]王成,鄭海濤,王惠文,等.中國通用航空產(chǎn)業(yè)規(guī)模的預(yù)測(cè)研究[J].生產(chǎn)力研究,2013(04):152-155.

    [3]中國人民共和國國家發(fā)展和改革委員會(huì).通用航空“十三五”發(fā)展規(guī)劃[EB/OL].http://www.caac.gov.cn/XXGK/XXGK/ZCFBJD/2017

    02/t20170217_42570.html,2017-02-17.

    [4]田磊,常松濤.通用航空事故高發(fā)原因分析[J].科技視界,2017(11):119+134.

    [5]宮峰勛.基于球形邊界的ATC飛行沖突概率預(yù)測(cè)[J].計(jì)算機(jī)工程與應(yīng)用,2005(28):177-178.

    [6]彭良福,顏悅,石佳奇,等.通航飛機(jī)避碰區(qū)域防撞模型的等效性研究[J].中國安全科學(xué)學(xué)報(bào),2016,26(08):95-99.

    [7]范盛金.一元三次方程的新求根公式與新判別法[J].海南師范學(xué)院學(xué)報(bào)(自然科學(xué)版),1989,2(2):91-98.

    [8]Gafiel M, Frazzoli E, Hansman R.J. Impact of GPS and ADS-B reported accuracy on conflict detection performance in dense.

    猜你喜歡
    通用航空
    通用航空飛行服務(wù)站運(yùn)營(yíng)模式研究
    科技資訊(2017年2期)2017-03-27 15:47:48
    SWOT矩陣在通航產(chǎn)業(yè)發(fā)展戰(zhàn)略研究中的應(yīng)用
    科技視界(2016年27期)2017-03-14 18:40:09
    大數(shù)據(jù)時(shí)代跨境電商與物流協(xié)同發(fā)展研究
    商情(2016年48期)2017-03-12 20:34:20
    我國通用航空發(fā)展問題與對(duì)策分析
    通用航空醫(yī)療運(yùn)營(yíng)模式分析
    商情(2016年40期)2016-11-28 10:25:29
    通用航空制造業(yè)及其發(fā)展
    通用航空維修保障能力的評(píng)估體系研究
    山東通用航空發(fā)展研究分析
    興趣教學(xué)模式在“通用航空運(yùn)行與保障”課程中的應(yīng)用探討
    通用航空典型機(jī)型CBT系統(tǒng)智能生成技術(shù)研究
    科技視界(2016年2期)2016-03-30 21:06:57
    av在线播放精品| 亚洲国产最新在线播放| 卡戴珊不雅视频在线播放| 欧美人与性动交α欧美软件 | 五月伊人婷婷丁香| 国产亚洲午夜精品一区二区久久| 成年人午夜在线观看视频| 久久人人97超碰香蕉20202| 免费看不卡的av| 大陆偷拍与自拍| 成人黄色视频免费在线看| 高清av免费在线| 国产免费视频播放在线视频| 九九爱精品视频在线观看| 欧美+日韩+精品| 男的添女的下面高潮视频| 成人国产麻豆网| 精品第一国产精品| 免费看av在线观看网站| 尾随美女入室| 免费人妻精品一区二区三区视频| 精品国产一区二区久久| 亚洲av电影在线进入| 亚洲国产成人一精品久久久| 美女视频免费永久观看网站| 成人二区视频| 婷婷成人精品国产| 亚洲性久久影院| 亚洲精品日韩在线中文字幕| 97在线视频观看| 国产片特级美女逼逼视频| 乱码一卡2卡4卡精品| 国产白丝娇喘喷水9色精品| 五月伊人婷婷丁香| 丰满少妇做爰视频| 精品亚洲成国产av| 国产av精品麻豆| 视频区图区小说| 久久 成人 亚洲| 999精品在线视频| 尾随美女入室| 成人国产麻豆网| 美女福利国产在线| 涩涩av久久男人的天堂| 成人亚洲欧美一区二区av| 又大又黄又爽视频免费| 大片免费播放器 马上看| 久久久久久久国产电影| 99精国产麻豆久久婷婷| 在线观看三级黄色| 国产精品99久久99久久久不卡 | 少妇人妻精品综合一区二区| 亚洲精品久久成人aⅴ小说| 日本黄大片高清| 国产一区亚洲一区在线观看| 国产精品一区www在线观看| 日韩一区二区视频免费看| 国产成人精品久久久久久| 国产欧美另类精品又又久久亚洲欧美| 色婷婷久久久亚洲欧美| 亚洲精品久久成人aⅴ小说| av网站免费在线观看视频| 一边亲一边摸免费视频| 韩国av在线不卡| 高清在线视频一区二区三区| 成人亚洲欧美一区二区av| 日日撸夜夜添| 国产老妇伦熟女老妇高清| 日本欧美视频一区| 亚洲精品乱码久久久久久按摩| 男女边吃奶边做爰视频| 热99国产精品久久久久久7| 精品人妻一区二区三区麻豆| 亚洲欧美一区二区三区黑人 | 中文乱码字字幕精品一区二区三区| 高清在线视频一区二区三区| 男的添女的下面高潮视频| 多毛熟女@视频| 国产成人aa在线观看| 久久久精品免费免费高清| 久久久久久久久久久久大奶| 制服人妻中文乱码| 国产深夜福利视频在线观看| 午夜福利在线观看免费完整高清在| 国产一区亚洲一区在线观看| 一二三四在线观看免费中文在 | 桃花免费在线播放| 曰老女人黄片| 亚洲久久久国产精品| 欧美人与性动交α欧美精品济南到 | 街头女战士在线观看网站| 成人漫画全彩无遮挡| 丰满饥渴人妻一区二区三| 免费看不卡的av| 国产精品无大码| 人人妻人人添人人爽欧美一区卜| 久久久精品区二区三区| 在线观看人妻少妇| 五月开心婷婷网| 1024视频免费在线观看| 亚洲国产欧美在线一区| 亚洲国产最新在线播放| 内地一区二区视频在线| 亚洲,欧美,日韩| 亚洲国产精品一区二区三区在线| 99香蕉大伊视频| 视频中文字幕在线观看| 又黄又粗又硬又大视频| 青春草国产在线视频| 国产又爽黄色视频| 久久女婷五月综合色啪小说| 观看av在线不卡| 中国三级夫妇交换| 一区二区三区四区激情视频| 亚洲精品日本国产第一区| 建设人人有责人人尽责人人享有的| 丁香六月天网| av线在线观看网站| 午夜视频国产福利| 亚洲久久久国产精品| 中文字幕另类日韩欧美亚洲嫩草| 亚洲精品视频女| 蜜臀久久99精品久久宅男| 久久这里只有精品19| 国产亚洲欧美精品永久| 欧美精品国产亚洲| 黄色怎么调成土黄色| 丰满少妇做爰视频| 一二三四在线观看免费中文在 | 人人妻人人澡人人爽人人夜夜| 亚洲欧美成人精品一区二区| 美女脱内裤让男人舔精品视频| 啦啦啦在线观看免费高清www| www日本在线高清视频| av黄色大香蕉| 男女边摸边吃奶| 制服人妻中文乱码| 丁香六月天网| 亚洲一区二区三区欧美精品| 男女边吃奶边做爰视频| 老熟女久久久| freevideosex欧美| 人妻人人澡人人爽人人| 午夜福利,免费看| 波多野结衣一区麻豆| 免费观看无遮挡的男女| 成年女人在线观看亚洲视频| 2022亚洲国产成人精品| 亚洲成人手机| 一级爰片在线观看| 日韩人妻精品一区2区三区| 中国国产av一级| 熟妇人妻不卡中文字幕| 精品福利永久在线观看| 亚洲,一卡二卡三卡| 国产乱人偷精品视频| 日本免费在线观看一区| 国产国拍精品亚洲av在线观看| 欧美日韩av久久| 高清视频免费观看一区二区| 亚洲av福利一区| 在线免费观看不下载黄p国产| 欧美+日韩+精品| 日日撸夜夜添| 成人二区视频| 精品福利永久在线观看| 丝袜美足系列| 亚洲第一av免费看| 免费人妻精品一区二区三区视频| 中文字幕精品免费在线观看视频 | 校园人妻丝袜中文字幕| 精品亚洲成国产av| 国产精品国产三级专区第一集| 久久青草综合色| 亚洲伊人色综图| 欧美精品人与动牲交sv欧美| 美女福利国产在线| 久久亚洲国产成人精品v| 国产麻豆69| 日韩免费高清中文字幕av| 男女无遮挡免费网站观看| 街头女战士在线观看网站| 咕卡用的链子| 高清视频免费观看一区二区| 国产亚洲欧美精品永久| 国产精品蜜桃在线观看| 亚洲一区二区三区欧美精品| 欧美人与性动交α欧美软件 | 精品国产国语对白av| 国产精品久久久av美女十八| 日本色播在线视频| 在现免费观看毛片| 亚洲精品美女久久久久99蜜臀 | 18+在线观看网站| 美女主播在线视频| 在线 av 中文字幕| 少妇高潮的动态图| 中文精品一卡2卡3卡4更新| 免费黄色在线免费观看| 欧美xxⅹ黑人| av卡一久久| 天天影视国产精品| 日本与韩国留学比较| 嫩草影院入口| 中国美白少妇内射xxxbb| 日韩制服骚丝袜av| a级毛片黄视频| 国产成人精品福利久久| 少妇的丰满在线观看| 少妇 在线观看| www.av在线官网国产| xxx大片免费视频| 三级国产精品片| 国产精品.久久久| 黑人猛操日本美女一级片| 中文字幕免费在线视频6| 久久人妻熟女aⅴ| 视频中文字幕在线观看| 国产xxxxx性猛交| 欧美日韩精品成人综合77777| 国产精品久久久久成人av| 在线精品无人区一区二区三| 欧美少妇被猛烈插入视频| 香蕉精品网在线| 色94色欧美一区二区| 一级片'在线观看视频| 久久久久久久国产电影| 香蕉国产在线看| 黄色毛片三级朝国网站| 国产精品久久久av美女十八| 亚洲国产av影院在线观看| 亚洲国产精品成人久久小说| 久久国产精品大桥未久av| 成人漫画全彩无遮挡| 飞空精品影院首页| 精品福利永久在线观看| 又黄又粗又硬又大视频| 夫妻性生交免费视频一级片| 黄色怎么调成土黄色| 国产色爽女视频免费观看| 欧美亚洲 丝袜 人妻 在线| 欧美日韩国产mv在线观看视频| 宅男免费午夜| 亚洲国产欧美在线一区| 国产亚洲精品第一综合不卡 | 乱码一卡2卡4卡精品| 国产 一区精品| av不卡在线播放| 亚洲,欧美,日韩| 亚洲内射少妇av| 久久久久久久国产电影| 亚洲熟女精品中文字幕| 亚洲五月色婷婷综合| 日本免费在线观看一区| 精品国产国语对白av| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲国产精品专区欧美| 欧美亚洲日本最大视频资源| 天天操日日干夜夜撸| 大片免费播放器 马上看| 免费在线观看完整版高清| 亚洲国产欧美在线一区| 一个人免费看片子| 欧美日韩精品成人综合77777| 51国产日韩欧美| 99精国产麻豆久久婷婷| 大片免费播放器 马上看| 午夜福利,免费看| 高清av免费在线| 高清在线视频一区二区三区| 日韩不卡一区二区三区视频在线| 成人二区视频| av播播在线观看一区| 丝袜脚勾引网站| 黑丝袜美女国产一区| 久久热在线av| 天天影视国产精品| 各种免费的搞黄视频| 在线观看www视频免费| 成年美女黄网站色视频大全免费| 热re99久久国产66热| 亚洲精品一二三| 国产国拍精品亚洲av在线观看| 亚洲伊人色综图| 国产 一区精品| 在线免费观看不下载黄p国产| 欧美xxⅹ黑人| 大香蕉97超碰在线| 亚洲成色77777| 精品少妇久久久久久888优播| av视频免费观看在线观看| 久久99精品国语久久久| 国产乱来视频区| 国产免费福利视频在线观看| 最近手机中文字幕大全| 最近中文字幕2019免费版| 日日爽夜夜爽网站| 黄色一级大片看看| 欧美性感艳星| 亚洲性久久影院| www.熟女人妻精品国产 | 9色porny在线观看| 在线观看人妻少妇| 日韩大片免费观看网站| 精品亚洲成国产av| 女人精品久久久久毛片| videos熟女内射| 亚洲欧美成人综合另类久久久| 最近的中文字幕免费完整| 亚洲欧美精品自产自拍| 亚洲四区av| 精品酒店卫生间| www.熟女人妻精品国产 | 日本爱情动作片www.在线观看| 赤兔流量卡办理| 一级片'在线观看视频| 老女人水多毛片| 国产熟女午夜一区二区三区| 婷婷色综合www| 国产一区二区三区综合在线观看 | 精品人妻一区二区三区麻豆| 亚洲成人av在线免费| 桃花免费在线播放| 欧美人与性动交α欧美软件 | 狂野欧美激情性xxxx在线观看| 久久国产精品男人的天堂亚洲 | 亚洲精华国产精华液的使用体验| 成人影院久久| 两个人免费观看高清视频| 一级片'在线观看视频| 久久久a久久爽久久v久久| 国产成人91sexporn| 另类精品久久| 日本午夜av视频| 欧美成人午夜精品| 18在线观看网站| 国产成人aa在线观看| 国产免费一级a男人的天堂| 精品久久久精品久久久| 欧美日韩一区二区视频在线观看视频在线| 又黄又粗又硬又大视频| 在线看a的网站| 亚洲国产最新在线播放| 你懂的网址亚洲精品在线观看| 精品少妇内射三级| 女性生殖器流出的白浆| 好男人视频免费观看在线| 成年人免费黄色播放视频| 狂野欧美激情性xxxx在线观看| 欧美变态另类bdsm刘玥| 美女国产视频在线观看| av线在线观看网站| 桃花免费在线播放| 亚洲欧美色中文字幕在线| 亚洲一级一片aⅴ在线观看| 久久久久久人妻| 国产 一区精品| 亚洲av综合色区一区| 在线观看免费日韩欧美大片| 婷婷成人精品国产| 久久精品久久久久久噜噜老黄| 久久久久久人人人人人| 99香蕉大伊视频| 国产精品99久久99久久久不卡 | 国产福利在线免费观看视频| 国产精品欧美亚洲77777| 国产精品蜜桃在线观看| 久久人妻熟女aⅴ| 视频中文字幕在线观看| 精品一区二区三卡| 亚洲色图 男人天堂 中文字幕 | 亚洲国产看品久久| 看非洲黑人一级黄片| a 毛片基地| 国产男女超爽视频在线观看| 狂野欧美激情性xxxx在线观看| 亚洲伊人色综图| 大香蕉久久成人网| 亚洲综合精品二区| 午夜福利,免费看| av又黄又爽大尺度在线免费看| 欧美3d第一页| 久久99一区二区三区| 久久精品久久精品一区二区三区| 国产又爽黄色视频| 黑人高潮一二区| 精品一区二区三卡| 激情五月婷婷亚洲| 人人妻人人爽人人添夜夜欢视频| 日本av手机在线免费观看| 日本黄色日本黄色录像| 久久国产亚洲av麻豆专区| 欧美激情极品国产一区二区三区 | 最近最新中文字幕大全免费视频 | 久久精品国产自在天天线| 成年av动漫网址| 久久久久精品久久久久真实原创| 新久久久久国产一级毛片| 桃花免费在线播放| 日韩精品免费视频一区二区三区 | 久久久久久久国产电影| 91aial.com中文字幕在线观看| 美女视频免费永久观看网站| 久久这里只有精品19| 新久久久久国产一级毛片| 一本—道久久a久久精品蜜桃钙片| 亚洲国产毛片av蜜桃av| 国产亚洲欧美精品永久| 97超碰精品成人国产| 久久影院123| 欧美变态另类bdsm刘玥| 中文精品一卡2卡3卡4更新| 久久精品国产a三级三级三级| 久久99精品国语久久久| 国产成人aa在线观看| 欧美日韩综合久久久久久| 五月玫瑰六月丁香| 中文字幕最新亚洲高清| tube8黄色片| 草草在线视频免费看| a级毛色黄片| 成人毛片60女人毛片免费| 精品少妇内射三级| 亚洲,欧美,日韩| 9色porny在线观看| 尾随美女入室| 精品一区在线观看国产| 精品人妻在线不人妻| 多毛熟女@视频| 日韩精品免费视频一区二区三区 | av国产久精品久网站免费入址| kizo精华| 成人国产av品久久久| 亚洲av中文av极速乱| 亚洲精品日本国产第一区| 黑丝袜美女国产一区| 大片免费播放器 马上看| 亚洲五月色婷婷综合| 永久免费av网站大全| 亚洲国产精品一区三区| 亚洲第一av免费看| 亚洲成国产人片在线观看| 久久精品国产鲁丝片午夜精品| 日韩制服骚丝袜av| 精品人妻一区二区三区麻豆| 美女国产视频在线观看| 国产乱人偷精品视频| 国产欧美日韩综合在线一区二区| 高清欧美精品videossex| 欧美最新免费一区二区三区| 亚洲,一卡二卡三卡| 国产在线视频一区二区| 免费在线观看黄色视频的| 国产69精品久久久久777片| 亚洲精品视频女| 熟女人妻精品中文字幕| 精品人妻偷拍中文字幕| 国产成人aa在线观看| 日本-黄色视频高清免费观看| 黄色毛片三级朝国网站| 国产成人精品福利久久| 高清不卡的av网站| 欧美97在线视频| 各种免费的搞黄视频| 日日啪夜夜爽| 26uuu在线亚洲综合色| 建设人人有责人人尽责人人享有的| 天天操日日干夜夜撸| 99久国产av精品国产电影| 天堂8中文在线网| 久久久久精品性色| 亚洲人成77777在线视频| 亚洲色图综合在线观看| 捣出白浆h1v1| 最近最新中文字幕免费大全7| 国产精品久久久久久久久免| 婷婷色麻豆天堂久久| 亚洲精品日本国产第一区| 亚洲人成77777在线视频| 久久99热这里只频精品6学生| 高清视频免费观看一区二区| 黑人欧美特级aaaaaa片| 久久国产精品男人的天堂亚洲 | 麻豆乱淫一区二区| a级毛片在线看网站| 男男h啪啪无遮挡| 蜜桃在线观看..| 18禁裸乳无遮挡动漫免费视频| 亚洲国产精品国产精品| 搡女人真爽免费视频火全软件| 久久久久久人人人人人| 免费在线观看黄色视频的| 另类精品久久| av在线观看视频网站免费| 国产日韩欧美视频二区| 国产成人精品无人区| 国产在视频线精品| 在线天堂中文资源库| 国产片内射在线| 国产色婷婷99| 亚洲国产欧美在线一区| 一级爰片在线观看| 成人国产av品久久久| av视频免费观看在线观看| 久久久久视频综合| 亚洲第一av免费看| 精品国产露脸久久av麻豆| 成人无遮挡网站| 不卡视频在线观看欧美| 亚洲av电影在线观看一区二区三区| 天天影视国产精品| 国产精品人妻久久久影院| 久久久a久久爽久久v久久| 成人国语在线视频| 亚洲国产毛片av蜜桃av| 女人精品久久久久毛片| 久久久国产欧美日韩av| 亚洲av.av天堂| 超碰97精品在线观看| 男人操女人黄网站| 建设人人有责人人尽责人人享有的| av在线观看视频网站免费| 国产亚洲一区二区精品| 18在线观看网站| 丰满饥渴人妻一区二区三| 久久免费观看电影| 丝袜脚勾引网站| freevideosex欧美| 欧美丝袜亚洲另类| av福利片在线| 极品少妇高潮喷水抽搐| 99久久中文字幕三级久久日本| 国产色爽女视频免费观看| 十分钟在线观看高清视频www| 精品国产一区二区三区四区第35| 一级毛片电影观看| 精品少妇黑人巨大在线播放| 丰满乱子伦码专区| 深夜精品福利| 亚洲高清免费不卡视频| 90打野战视频偷拍视频| 国产成人免费无遮挡视频| 九九在线视频观看精品| 天美传媒精品一区二区| 免费看光身美女| 香蕉国产在线看| 久久久久久久大尺度免费视频| 亚洲伊人久久精品综合| 三级国产精品片| 久久精品久久久久久噜噜老黄| 内地一区二区视频在线| 精品久久蜜臀av无| av在线老鸭窝| 国产精品免费大片| 欧美人与善性xxx| 好男人视频免费观看在线| 男女免费视频国产| 一二三四在线观看免费中文在 | 青青草视频在线视频观看| 久久久国产精品麻豆| 大话2 男鬼变身卡| 又粗又硬又长又爽又黄的视频| 秋霞伦理黄片| 亚洲精品乱久久久久久| 十八禁高潮呻吟视频| 亚洲综合色惰| 亚洲少妇的诱惑av| 国产精品偷伦视频观看了| 亚洲精品456在线播放app| 国产国语露脸激情在线看| 亚洲图色成人| 亚洲精品乱久久久久久| 久久国产精品男人的天堂亚洲 | 精品一区在线观看国产| 内地一区二区视频在线| 成人国语在线视频| 日韩欧美一区视频在线观看| 国产精品人妻久久久影院| 热99久久久久精品小说推荐| 黄片播放在线免费| 免费大片黄手机在线观看| 国产在视频线精品| 久久免费观看电影| 美女视频免费永久观看网站| www.色视频.com| 日韩精品免费视频一区二区三区 | 久久久久久久大尺度免费视频| 美女国产视频在线观看| 精品一品国产午夜福利视频| 9色porny在线观看| 免费av不卡在线播放| a级毛片黄视频| 超色免费av| 国产一级毛片在线| av片东京热男人的天堂| 亚洲欧美色中文字幕在线| 伦理电影免费视频| 在线天堂中文资源库| 亚洲精品国产av成人精品| 亚洲一区二区三区欧美精品| 在线亚洲精品国产二区图片欧美| 亚洲精品乱久久久久久| 欧美老熟妇乱子伦牲交| 99久久综合免费| 母亲3免费完整高清在线观看 | 欧美日韩av久久| 哪个播放器可以免费观看大片| 天天躁夜夜躁狠狠躁躁| 日日撸夜夜添| 熟女av电影| 美国免费a级毛片| 欧美日韩成人在线一区二区| 亚洲国产精品国产精品| 丰满乱子伦码专区| 国产白丝娇喘喷水9色精品|