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

    子攔截彈攔截無人機集群防碰撞制導(dǎo)律

    2023-05-30 00:42:12羅瑞寧黃樹彩趙巖張振張飛
    航空兵器 2023年1期

    羅瑞寧 黃樹彩 趙巖 張振 張飛

    引用格式:羅瑞寧,黃樹彩,趙巖,等.子攔截彈攔截無人機集群防碰撞制導(dǎo)律[J].航空兵器,2023,30(1):51-58.

    LuoRuining,HuangShucai,ZhaoYan,etal.CollisionAvoidanceGuidanceLawforSon-InterceptorsInterceptingUAVCluster[J].AeroWeaponry,2023,30(1):51-58.(inChinese)

    摘要:為解決子攔截彈攔截無人機集群過程中可能發(fā)生的碰撞問題,本文將虛擬排斥力與純比例導(dǎo)引律結(jié)合,提出一種防碰撞純比例導(dǎo)引制導(dǎo)律(CollisionAvoidance-PureProportionalGuidanceLaw,CA-PPN)。首先,提出一種子母導(dǎo)彈攔截無人機集群的作戰(zhàn)構(gòu)想,并分析了其中的子攔截彈制導(dǎo)問題。隨后,針對子攔截彈制導(dǎo)中可能發(fā)生的碰撞問題,在子攔截彈之間設(shè)計了一種用于避碰的人工勢場,分析了PPN的特性,在此基礎(chǔ)上將人工勢場下的虛擬排斥力與PPN結(jié)合,并加入視場角與過載約束,提出了CA-PPN。仿真結(jié)果表明,CA-PPN能在制導(dǎo)子攔截彈有效攔截目標的同時避免相互之間的碰撞。

    關(guān)鍵詞:防碰撞制導(dǎo)律;虛擬排斥力;純比例導(dǎo)引;人工勢場;子攔截彈;導(dǎo)彈;無人機集群

    中圖分類號:TJ756;V279

    文獻標識碼:A

    文章編號:1673-5048(2023)01-0051-08

    DOI:10.12132/ISSN.1673-5048.2022.0117

    0引言

    隨著多彈協(xié)同作戰(zhàn)技術(shù)的發(fā)展,多導(dǎo)彈防碰撞已成為多彈編隊飛行、多彈協(xié)同制導(dǎo)過程中不可忽視的一個問題。

    防碰撞問題較早出現(xiàn)于以多機器人為代表的多智能體系統(tǒng),為實現(xiàn)機器人運動中的碰撞規(guī)避,學(xué)者們應(yīng)用最多的是一種基于人工勢場(或稱為虛擬力)的防碰撞方法。人工勢場的方法根據(jù)分子間因距離遠近產(chǎn)生相應(yīng)吸引力與排斥力的物理原理,在機器人之間加入一種隨距離而變化的虛擬力,實現(xiàn)多機器人之間的碰撞規(guī)避與編隊保持。人工勢場的方法在多智能體系統(tǒng)的避障、防碰撞以及編隊保持中均有良好應(yīng)用。文獻[1-3]中將虛擬力與其他控制輸出線性疊加共同控制機器人位移,在實現(xiàn)其各自功能的同時能夠保持隊形且避免碰撞,但其應(yīng)用對象均為自由度較高的機器人,并不能直接應(yīng)用于受到非完整約束的固定翼導(dǎo)彈。文獻[4-5]提出一種旋轉(zhuǎn)的勢場,排斥力與機器人的運動相適應(yīng),且勢場旋轉(zhuǎn)方向由機器人運動方向決定,旋轉(zhuǎn)勢場的設(shè)計更有利于機器人通過障礙且避免了局部最小值的出現(xiàn)。文獻[6]設(shè)計了一種虛擬力,可同時保持多智能體連通與防碰撞,并給出了約束條件的精確表達式與控制力的估計值。文獻[7-8]在無人機之間設(shè)置虛擬力,由虛擬力的大小得到相應(yīng)的位移量,實現(xiàn)規(guī)避障礙及無人機編隊防碰撞,但由于其控制輸出的是位移量,可適用于旋翼無人機,并不能直接應(yīng)用于導(dǎo)彈。文獻[9]在設(shè)計虛擬力時,除考慮相對距離外還加入了相對速度因素。文獻[10-11]針對無人機部署優(yōu)化問題,對典型的虛擬力進行了改進,提高了虛擬力模型精度。

    除基于人工勢場的防碰撞方法外,優(yōu)化的方法也被應(yīng)用在避障與防碰撞中,如文獻[12]用粒子群算法搜索不會發(fā)生碰撞的最短路徑;文獻[13-14]用粒子群與遺傳算法的融合算法求解線性目標函數(shù)的優(yōu)化問題;文獻[15]用線性矩陣表示防碰撞約束,并給出一種用逆動態(tài)學(xué)求可行解的方法;文獻[16]利用凸優(yōu)化的方法將不可微分約束轉(zhuǎn)化為光滑的可微分約束,進而用常規(guī)的函數(shù)求最值方法獲得最優(yōu)解。但優(yōu)化類方法通常求解復(fù)雜,計算量較大,不能確保短時間內(nèi)收斂,難以有效應(yīng)用在飛行速度快、需要快速做出響應(yīng)的導(dǎo)彈類目標的防碰撞設(shè)計中。除優(yōu)化方法外,文獻[17-18]用神經(jīng)網(wǎng)絡(luò)的方法解決無人機的避障與路徑規(guī)劃問題。文獻[19]模仿果蠅根據(jù)視覺特性向障礙物反方向偏轉(zhuǎn)來規(guī)避障礙物的原理,提出一種仿生的避障方法,但僅適用于規(guī)避具有一定體積的固定障礙物,并不能應(yīng)用于彈間防碰撞。

    綜上,人工勢場的方法具有所需信息少、簡單有效,且易與其他控制器結(jié)合的優(yōu)點,但應(yīng)用于多導(dǎo)彈間防碰撞的研究較少,且存在勢場力維數(shù)與導(dǎo)彈制導(dǎo)輸出維數(shù)不匹配的問題。

    針對子攔截彈攔截無人機集群過程中可能出現(xiàn)的彈間碰撞問題,本文將人工勢場的方法引入到多子攔截彈的防碰撞設(shè)計中,提出一種用于防碰撞的勢場力與純比例導(dǎo)引律結(jié)合的方法,得到一種防碰撞制導(dǎo)律,可在多子攔截彈攔截集群目標過程中有效規(guī)避彈間碰撞。

    1問題描述

    子攔截彈攔截無人機集群,是子母導(dǎo)彈攔截無人機集群作戰(zhàn)構(gòu)想中的最后一環(huán)。這里首先對子母導(dǎo)彈攔截無人機集群作戰(zhàn)構(gòu)想進行簡要描述,再重點分析其中的子攔截彈制導(dǎo)問題。

    1.1子母導(dǎo)彈反無人機集群作戰(zhàn)構(gòu)想

    針對中遠程階段反無人機集群問題,提出一種子母導(dǎo)彈攔截無人機集群的作戰(zhàn)構(gòu)想。面向的集群目標是由“小精靈”無人機等用于集群作戰(zhàn)的中/大型無人機組成的集群,此類目標可在防區(qū)外釋放,具有較大的作戰(zhàn)半徑,可執(zhí)行偵察、干擾、打擊任務(wù),具有較大的威脅等級與攔截效益。

    所設(shè)計子母導(dǎo)彈由母艙與運載在艙內(nèi)的子攔截彈構(gòu)成,如圖1所示。母艙搭載渦噴發(fā)動機或固體火箭發(fā)動機,以獲得遠大于無人機的飛行速度;采用復(fù)合式導(dǎo)引頭,在飛至集群目標一定距離后可對集群目標進行識別與截獲,并在適當?shù)奈恢靡詮椛浞绞结尫抛訑r截彈對無人機目標進行攔截。子攔截彈攜帶光學(xué)與紅外導(dǎo)引頭,通過智能識別算法識別并鎖定目標,尾部裝有電動舵系統(tǒng),依靠母艙提供的較大初速度制導(dǎo)攔截目標。

    子母導(dǎo)彈攔截無人機集群的基本作戰(zhàn)過程為:分散布置的綜合情報系統(tǒng)對集群目標進行識別與跟蹤,并將目標信息指示給子母導(dǎo)彈指揮控制系統(tǒng);指揮控制系統(tǒng)引導(dǎo)導(dǎo)彈飛向目標,并在適當?shù)奈恢?,母艙釋放子攔截彈對目標進行攔截。作戰(zhàn)構(gòu)想如圖2所示。

    圖2中的威脅“重心”是用來推導(dǎo)“釋放點”的位置。本文重點分析子攔截彈的制導(dǎo)問題,而不展開討論子母導(dǎo)彈的母艙制導(dǎo)與“釋放點”選取問題。

    1.2子攔截彈制導(dǎo)問題分析

    子攔截彈的速度遠大于無人機的速度,且攔截時距離較近,而無人機目標機動能力較差,加之受限于集群編隊,短時間內(nèi)可做出的機動有限,因此子攔截彈攔截無人機目標過程中可忽略目標機動因素,并假設(shè):子攔截彈攔截無人機目標階段,目標不進行機動。

    在高速攔截器攔截低速目標的情況下,選取純比例導(dǎo)引制導(dǎo)律(PPN)可獲得較好的制導(dǎo)效果,且易于工程實現(xiàn)。

    子母導(dǎo)彈飛抵釋放點時,以彈射的方式釋放子攔截彈,子攔截彈在制導(dǎo)初始時刻已分散在母艙周圍且具有一定間距,通過合理的目標分配,可實現(xiàn)子攔截彈優(yōu)先攔截各自方位目標,避免航跡交叉。但子攔截彈攔截目標的過程中仍可能會產(chǎn)生以下兩種不利情形:

    (1)子攔截彈目標分配不合理、兩目標距離較近或空間位置呈前后關(guān)系,均易導(dǎo)致子攔截彈軌跡發(fā)生交叉。

    (2)由于無人機集群的異構(gòu)特性,目標分配結(jié)果可能存在2個子攔截彈攔截1個高價值無人機目標的情況。

    以上兩種情形下,子攔截彈之間都存在碰撞或相互干擾的風(fēng)險,因此有必要進行防碰撞設(shè)計。

    2防碰撞制導(dǎo)律設(shè)計

    期望在PPN的制導(dǎo)輸出上加入防碰撞的偏置項,使子攔截彈既能在PPN的作用下獲得較好的彈道,又能在即將發(fā)生碰撞的情況下及時作出規(guī)避。

    在導(dǎo)彈之間引入人工勢場是一種簡單且有效的防碰撞手段,為實現(xiàn)PPN與人工勢場的有效結(jié)合,首先對多彈間人工勢場與PPN的性能特點進行分析。

    2.1多導(dǎo)彈間人工勢場

    規(guī)定子攔截彈間的安全距離為rsafe,發(fā)生碰撞或干擾的臨界距離為rmin。根據(jù)人工勢場的原理,子攔截彈i與j之間距離rij>rsafe時,彈間不存在勢場;rmin≤rij≤rsafe時,彈間出現(xiàn)勢場,并產(chǎn)生排斥力;rij

    依據(jù)分子間距離較近時產(chǎn)生相互排斥力的物理原理,基于彈間距離設(shè)計如下人工勢場函數(shù):

    Uij(rij)=0rij>rsafeK2(rsafe-rij)2rmin≤rij≤rsafe

    m∞r(nóng)ij

    式中:m∞為無窮大常數(shù);K為勢場系數(shù),需根據(jù)具體應(yīng)用場景進行設(shè)置。

    對式(1)求偏微分,得到子攔截彈i受到子攔截彈j的排斥力Fij為

    Fij(xi,yi,zi)=-Uij(rij)=0rij>rsafe

    K(rij-rsafe)rijxirijyirijziTrmin≤rij≤rsafe

    ∞r(nóng)ij

    式中:rijxirijyirijziT為單位矢量,由子攔截彈j指向子攔截彈i。

    多彈情形下,子攔截彈i可能同時受到多個來自其他子攔截彈的排斥力,定義由n枚子攔截彈構(gòu)成的多彈系統(tǒng)中,子攔截彈i在t時刻的防碰撞鄰域Ci(t)為

    Ci(t)={j|rij≤rsafe,j≠i,j=1,2,…,n}(3)

    因此,可得到子攔截彈i受到的總排斥力Fi-rep為

    Fi-rep(xi,yi,zi)=∑j∈Ci(t)K(rij-rsafe)rijxirijyirijziT(4)

    觀察式(4)可知,為避免碰撞需對子攔截彈施加一個三維的控制量(二維平面內(nèi)控制量為二維),而常規(guī)戰(zhàn)術(shù)導(dǎo)彈的制導(dǎo)輸出僅為二維(二維平面內(nèi)制導(dǎo)輸出為一維),因此排斥力Fi-rep不能直接附加在制導(dǎo)輸出上。為實現(xiàn)PPN制導(dǎo)律與防碰撞控制量的有效結(jié)合,需進一步分析PPN的特性。

    對于一般戰(zhàn)術(shù)導(dǎo)彈,制導(dǎo)輸出可分解為相互垂直的兩個方向,兩個方向的控制器設(shè)計原理相同,因此這里均針對二維平面制導(dǎo)輸出為一維的情況進行討論。

    2.2PPN特性分析

    縱向平面內(nèi)彈目相對幾何關(guān)系如圖3所示。XOZ為慣性坐標系,VM,aM分別為導(dǎo)彈的速度矢量與加速度矢量,q為導(dǎo)彈M對目標T的視線角,η為導(dǎo)彈速度前置角,er,eθ分別為沿視線方向與垂直于視線方向的單位矢量。規(guī)定視線角q由水平正向起始,逆時針旋轉(zhuǎn)至彈目視線方向為正;導(dǎo)彈速度前置角η由導(dǎo)彈速度方向起始,逆時針旋轉(zhuǎn)至彈目視線方向為正。

    以導(dǎo)彈M位置為原點,以er,eθ矢量方向為兩坐標軸方向,建立二維動坐標系erMeθ[20],彈目相對運動關(guān)系式如下:

    r¨-rω2=aTr-aMr

    rω·+2r·ω=aTθ-aMθ(5)

    式中:r為彈目相對距離;ω為視線角速率;aMr,aTr分別為導(dǎo)彈與目標加速度在視線方向的分量,aMθ,aTθ分別為導(dǎo)彈與目標加速度在垂直于速度方向的分量。

    PN可統(tǒng)一描述為

    aPN=NL×ω(6)

    式中:N為比例系數(shù);L為參考速度矢量;ω為視線角速度矢量,ω=q·ez,ez為垂直于(er,eθ)平面向外的單位矢量,即ez=er×eθ。

    對于PPN而言,式(6)中的L=-VM,因此PPN的制導(dǎo)加速度矢量可表示為

    aPPN=N(-VM)×ω(7)

    由幾何關(guān)系可以得到

    VM=-(r·er+rq·eθ)(8)

    將式(8)代入式(7)中可以得到

    aPPN=Nq·(rq·er-r·eθ)(9)

    在aT=0的條件下,結(jié)合式(9),可將式(5)改寫為

    r¨-rq·2=-Nr·q2rq¨+2r·q·=Nr·q·(10)

    由式(10)中的第2個公式可以得到

    q·=q·0rr0N-2(11)

    在1.2節(jié)的假設(shè)條件下,aT=0,此時有vMsinη=-rq·,結(jié)合式(11)可以得到

    vMsinη=-rq·=-r0q·0rr0N-1(12)

    有初始條件vM0sinη0=-r0q·0,由此可以得到

    sinη=sinη0rr0N-1(13)

    由式(13)可知,PPN使導(dǎo)彈速度前置角朝減小的方向收斂。

    2.3基于人工勢場的防碰撞制導(dǎo)方法

    縱向平面內(nèi),子攔截彈之間的排斥力作用關(guān)系如圖4所示,F(xiàn)為子攔截彈之間的排斥力矢量,ζ為排斥力F與速度V之間的夾角,規(guī)定ζ由排斥力方向起始,逆時針旋轉(zhuǎn)至子攔截彈速度方向為正。

    彈間排斥力F對子攔截彈速度前置角η的影響只有兩種結(jié)果,即促進或抑制η的收斂。由圖4可知,當ζ與η同向(正負相同)時,F(xiàn)抑制η收斂;當ζ與η反向(正負相反)時,F(xiàn)促進η收斂。分析可知,PPN制導(dǎo)下的彈道使速度前置角η不斷減小,因此排斥力F的作用可表現(xiàn)為促進或抑制PPN彈道的收斂。

    基于以上分析,可將彈間排斥力F轉(zhuǎn)化為標量F,作為偏置項直接疊加在PPN制導(dǎo)輸出上,用Jij表示排斥力Fij對子攔截彈i制導(dǎo)輸出的作用方向,則有

    Jij=-1ifsgn(ζij)·sgn(ηi)≥0+1ifsgn(ζij)·sgn(ηi)<0(14)

    據(jù)此,將人工勢場與PPN結(jié)合,得到CA-PPN的表達式為

    ai=aP,i+aU,i(15)

    式中:aP,i為PPN制導(dǎo)輸出;aU,i為人工勢場作用下的防碰撞偏置量,計算公式為

    aU,i=∑j∈CiJij·Fij(16)

    在仿真實驗中發(fā)現(xiàn),當速度方向偏轉(zhuǎn)至與視線方向重合(即導(dǎo)彈速度前置角ηi=0)時,Jij的判定作用失效,易出現(xiàn)Jij的數(shù)值在-1與+1之間不斷切換,控制器輸出產(chǎn)生振蕩的問題。對此,需對Jij進行修正。設(shè)定當η∈(-1°,1°)時,Jij的正負進行如下判定:

    Jij=-1ifsgn(ζij)≥0+1ifsgn(ζij)<0(17)

    在子攔截彈速度前置角較小的情況下,具備較好的攔截適宜性,在式(17)的作用下,不再考慮對比例導(dǎo)引律彈道收斂的影響,僅進行碰撞規(guī)避。

    同時,考慮到子攔截彈的導(dǎo)引頭視場角覆蓋范圍有限,比例導(dǎo)引具有使導(dǎo)彈視場始終追蹤目標的良好特性,但加入防碰撞偏置項后,可能會存在導(dǎo)引頭視場偏離目標的風(fēng)險。為了對子攔截彈視場進行約束,對制導(dǎo)律式(15)進行修正[21],得到最終的CA-PPN的表達式為

    ai=aP,i+cosπ2·ηηmaxaU,i(18)

    式中:ηmax為子攔截彈最大視場角。如此設(shè)計可在視場發(fā)生偏離的時候削弱防碰撞偏置項的作用,目標接近視場邊緣時,偏置項的值趨于零,子攔截彈可在PPN作用下及時偏轉(zhuǎn)重新對準目標。修正后的制導(dǎo)律可確保子攔截彈在不丟失目標的前提下進行碰撞規(guī)避。

    注意到,當出現(xiàn)2個子攔截彈共同攔截1個高價值目標時,在PPN導(dǎo)引下,2枚彈的彈道相似且接近,2枚彈很可能因規(guī)避碰撞而偏離各自的彈道,均不能以最佳的PPN彈道攔截目標。對此,在m枚子攔截彈攔截同1個目標的情況下,將其指定1枚主攔截彈與m-1枚次攔截彈,主攔截彈不受其余m-1枚次攔截彈勢場影響,期望在PPN作用下獲得最佳攔截效果,次攔截彈在勢場作用下主動規(guī)避主攔截彈與其余次攔截彈。另外,假設(shè)在彈目距離較近時,取消主、次攔截彈之間的排斥力作用,各彈因存在初始間距,后續(xù)將實現(xiàn)以不同角度攔截目標而不會發(fā)生碰撞問題。因此,規(guī)定在距離目標500m時取消主、次攔截彈之間的虛擬排斥力,以避免排斥力對次攔截彈攔截效果的影響。

    最后,考慮子攔截彈實際輸出能力有限,需對CA-PPN的制導(dǎo)輸出設(shè)置如下過載約束:

    ai=aiifabs(ai)≤nmax

    sgn(ai)·nmaxifabs(ai)>nmax(19)

    式中:nmax為子攔截彈可輸出的極限過載。

    防碰撞制導(dǎo)律本質(zhì)上是一種協(xié)同制導(dǎo)律,制導(dǎo)律的實現(xiàn)除了目標信息外,還都需要導(dǎo)彈鄰域內(nèi)其他導(dǎo)彈的信息,不同之處在于協(xié)同制導(dǎo)中參與協(xié)同的導(dǎo)彈固定,且須保持實時通信,而本文防碰撞制導(dǎo)中只需要與進入防碰撞鄰域Ci(t)內(nèi)的子攔截彈進行協(xié)同,通信能力只需保證防碰撞鄰域Ci(t)內(nèi)的子攔截彈進行信息交互傳遞。鑒于此,采用文獻[22]中提出的雙層協(xié)同制導(dǎo)架構(gòu)設(shè)計本文防碰撞制導(dǎo)架構(gòu),如圖5所示,協(xié)調(diào)變量為防碰撞鄰域Ci(t)內(nèi)彈間相對位置信息。

    該防碰撞制導(dǎo)架構(gòu)分為上下兩層,上層為協(xié)調(diào)層,用來獲取協(xié)調(diào)變量ζij,ηi,rij,并在防碰撞協(xié)調(diào)策略下得到各子攔截彈的制導(dǎo)律偏置項aU,i傳遞給控制層;下層為控制層,PPN制導(dǎo)輸出aP,i疊加協(xié)調(diào)層輸出的偏置項aU,i形成本地制導(dǎo)律ai。

    3仿真分析

    規(guī)定子攔截彈之間的安全距離rsafe=50m,發(fā)生碰撞的臨界距離rmin=5m,人工勢場系數(shù)K=2,PPN制導(dǎo)律的比例系數(shù)N=3,子攔截彈視場角ηmax=15°,輸出極限過載nmax=20g,g=9.8m/s2。

    構(gòu)造一種會發(fā)生碰撞的彈目初始狀態(tài),如表1所示。其中,目標1為高價值目標,由子攔截彈1和2共同攔截,子攔截彈3攔截目標2。仿真步長設(shè)為0.01s,分為目標不機動與目標做圓弧機動兩種情形對所設(shè)計的防碰撞制導(dǎo)律進行驗證。

    情形一:目標不機動。

    PPN作用下的攔截效果如圖6所示。

    由圖6可知,在1.1~1.4s內(nèi),子攔截彈1和3、子攔截彈2和3的最近距離均在5m附近,隨著接近目標,子攔截彈1和2逐漸接近,攔截末段距離小于5m。因此,目標不機動時,僅在PPN的作用下,子攔截彈之間存在碰撞或相互干擾的風(fēng)險。在CA-PPN作用下攔截效果如圖7所示。

    觀察圖7(a),由于仿真步長的影響,攔截點并未落在目標軌跡上,但觀察彈目軌跡可知,3枚子攔截彈均能以較小的脫靶量攔截目標。由圖7(b)可知,子攔截彈1和3、子攔截彈2和3之間的最近距離均大于10m;子攔截彈1和2的初始距離較近,在虛擬排斥力的作用下,子攔截彈2主動規(guī)避子攔截彈1,使其距離不斷增大,除最后500m共同攔截目標外,其全程距離均遠大于5m。由圖7(c)可知,制導(dǎo)全程子攔截彈的速度前置角均小于15°,可保證目標始終落在視場范圍內(nèi)。由圖7(d)可知,在2.5~3s內(nèi),因子攔截彈1和2先后規(guī)避子攔截彈3,制導(dǎo)輸出出現(xiàn)2次抖動,其余時刻制導(dǎo)輸出較為平滑且均在極限過載之內(nèi)。綜上可知,目標不機動時,在CA-PPN作用下,子攔截彈可有效攔截目標,且相互之間能夠進行碰撞規(guī)避。

    情形二:目標做加速度為5m/s2的圓弧機動。

    在PPN作用下攔截效果如圖8所示。在CA-PPN作用下攔截效果如圖9所示。

    綜合圖8~9可知,目標做圓弧機動時,僅在PPN的作用下,子攔截彈之間存在碰撞或相互干擾的風(fēng)險;而在CA-PPN作用下,子攔截彈仍可有效攔截目標,其相互之間能夠規(guī)避碰撞。攔截效果與分析過程均與目標不機動時接近,證明所設(shè)計的CA-PPN可有效應(yīng)對機動目標。

    對比圖6(a)與7(a)、圖8(a)與9(a)可以看出,3枚子攔截彈均在需要規(guī)避碰撞時偏離了PPN彈道,飛出碰撞鄰域后仍以PPN制導(dǎo),整個彈道體現(xiàn)出PPN彈道的基本特性。同時,為體現(xiàn)CA-PPN的良好性能,仿真部分設(shè)置了需要進行多次碰撞規(guī)避的較為嚴苛的初始條件,在一般條件下僅需進行簡單碰撞規(guī)避時,可認為其彈道特征會更接近于PPN彈道。因此,在比例導(dǎo)引律的基礎(chǔ)上加入防碰撞偏置項的設(shè)計,能在有效規(guī)避碰撞的同時,很大程度上保留PPN的最優(yōu)特性。

    觀察圖7(d)、圖9(d)可發(fā)現(xiàn),CA-PPN作用下的制導(dǎo)輸出并不完全光滑,進行碰撞規(guī)避時刻制導(dǎo)輸出會發(fā)生一次跳躍。產(chǎn)生這種輸出結(jié)果一定程度上歸因于制導(dǎo)仿真的維數(shù)為二維,尤其當2枚子攔截彈的軌跡在某一階段近似對稱時,相互規(guī)避易產(chǎn)生如圖10所示的振蕩。

    此類振蕩現(xiàn)象加大了二維平面內(nèi)碰撞規(guī)避的難度,而在三維攔截場景中子攔截彈可在三個方向機動進行規(guī)避,能很大程度避免軌跡振蕩問題。因此,在實際制導(dǎo)攔截或三維仿真中,本文提出的CA-PPN的性能可表現(xiàn)得更好。

    綜上所述,證明本文提出的CA-PPN制導(dǎo)律可在有效攔截機動與非機動目標的同時,避免子攔截彈之間可能發(fā)生的碰撞。

    4結(jié)論

    本文將用于規(guī)避碰撞的虛擬排斥力與PPN結(jié)合,提出了一種用于規(guī)避彈間碰撞的CA-PPN制導(dǎo)律。所提出的CA-PPN能夠有效規(guī)避子攔截彈之間可能發(fā)生的碰撞,同時還能確保子攔截彈在始終不丟失目標的條件下以較小脫靶量攔截目標,為子攔截彈攔截集群目標提供了一種有效的制導(dǎo)方法。

    注意到,為簡潔直觀進行二維平面上的制導(dǎo)律設(shè)計與仿真驗證,應(yīng)用本文設(shè)計原理分別在縱向通道與偏航通道設(shè)計控制器,即可得到三維空間的CA-PPN,且理論上其防碰撞效果應(yīng)優(yōu)于二維平面的仿真結(jié)果。

    參考文獻:

    [1]ContiR,MeliE,RidolfiA,etal.AnInnovativeDecentralizedStrategyforI-AUVsCooperativeManipulationTasks[J].RoboticsandAutonomousSystems,2015,72:261-276.

    [2]PoonawalaHA,SaticiAC,EckertH,etal.Collision-FreeFormationControlwithDecentralizedConnectivityPreservationforNonholonomic-WheeledMobileRobots[J].IEEETransactionsonControlofNetworkSystems,2015,2(2):122-130.

    [3]RingbackR,WeiJQ,ErstorpES,etal.Multi-AgentFormationTrackingforAutonomousSurfaceVehicles[J].IEEETransactionsonControlSystemsTechnology,2020(99):1-12.

    [4]RezaeeH,AbdollahiF.ADecentralizedCooperativeControlSchemewithObstacleAvoidanceforaTeamofMobileRobots[J].IEEETransactionsonIndustrialElectronics,2014,61(1):347-354.

    [5]HwangJ,LeeJ,ParkC.CollisionAvoidanceControlforFormationFlyingofMultipleSpacecraftUsingArtificialPotentialField[J].AdvancesinSpaceResearch,2022,69(5):2197-2209.

    [6]ChenZY,F(xiàn)anMC,ZhangHT.HowMuchControlisEnoughforNetworkConnectivityPreservationandCollisionAvoidance?[J].IEEETransactionsonCybernetics,2015,45(8):1647-1656.

    [7]WangX,TanGZ,DaiYS,etal.AnOptimalGuidanceStrategyforMoving-TargetInterceptionbyaMultirotorUnmannedAerialVehicleSwarm[J].IEEEAccess,2020,8:121650-121664.

    [8]ZhaoHT,LiuH,LeungYW,etal.Self-AdaptiveCollectiveMotionofSwarmRobots[J].IEEETransactionsonAutomationScienceandEngineering,2018,15(4):1533-1545.

    [9]QiJT,GuoJJ,WangMM,etal.FormationTrackingandObstacleAvoidanceforMultipleQuadrotorswithStaticandDynamicObstacles[J].IEEERoboticsandAutomationLetters,2022,7(2):1713-1720.

    [10]WangX,TanGZ,LiuXY,etal.AMolecularForce-BasedDeploymentAlgorithmforFlightCoverageMaximizationofMulti-RotorUAV[J].JournalofIntelligent&RoboticSystems,2019,95(3/4):1063-1078.

    [11]WangX,TanGZ,LuFL,etal.AMolecularForceField-BasedOptimalDeploymentAlgorithmforUAVSwarmCoverageMaximizationinMobileWirelessSensorNetwork[J].Processes,2020,8(3):369.

    [12]I·nanT,BabaAF.ParticleSwarmOptimization-BasedCollisionAvoidance[J].TurkishJournalofElectricalEngineering&ComputerSciences,2019,27(3):2137-2155.

    [13]DuanHB,LiuSQ.Non-LinearDual-ModeRecedingHorizonControlforMultipleUnmannedAirVehiclesFormationFlightBasedonChaoticParticleSwarmOptimisation[J].IETControlTheory&Applications,2010,4(11):2565-2578.

    [14]DuanHB,LuoQN,ShiYH,etal.HybridParticleSwarmOptimizationandGeneticAlgorithmforMulti-UAVFormationReconfiguration[J].IEEEComputationalIntelligenceMagazine,2013,8(3):16-27.

    [15]GiamouM,Maric′F,RosenDM,etal.ConvexIterationforDistance-GeometricInverseKinematics[J].IEEERoboticsandAutomationLetters,2022,7(2):1952-1959.

    [16]ZhangXJ,LinigerA,BorrelliF.Optimization-BasedCollisionAvoidance[J].IEEETransactionsonControlSystemsTechnology,2021,29(3):972-983.

    [17]WangXH,YadavV,BalakrishnanSN.CooperativeUAVFormationFlyingwithObstacle/CollisionAvoidance[J].IEEETransactionsonControlSystemsTechnology,2007,15(4):672-679.

    [18]HeRK,WeiRX,ZhangQR.UAVAutonomousCollisionAvoidanceApproach[J].Automatika,2017,58(2):195-204.

    [19]FabianST,SumnerME,WardillTJ,etal.AvoidingObstaclesWhileInterceptingaMovingTarget:AMiniatureFlysSolution[J].TheJournalofExperimentalBiology,2022,225(4):jeb243568.

    [20]黎克波,廖選平,梁彥剛,等.基于純比例導(dǎo)引的攔截碰撞角約束制導(dǎo)策略[J].航空學(xué)報,2020,41(S2):724277.

    LiKebo,LiaoXuanping,LiangYangang,etal.GuidanceStrategywithImpactAngleConstraintsBasedonPureProportionalNavigation[J].ActaAeronauticaetAstronauticaSinica,2020,41(S2):724277.(inChinese)

    [21]ZhangYA,WangXL,WuHL.ImpactTimeControlGuidancewithField-of-ViewConstraintAccountingforUncertainSystemLag[J].JournalofAerospaceEngineering,2016,230(3):515-529.

    [22]趙世鈺,周銳.基于協(xié)調(diào)變量的多導(dǎo)彈協(xié)同制導(dǎo)[J].航空學(xué)報,2008,29(6):1605-1611.

    ZhaoShiyu,ZhouRui.Multi-MissileCooperativeGuidanceUsingCoordinationVariables[J].ActaAeronauticaetAstronauticaSinica,2008,29(6):1605-1611.(inChinese)

    CollisionAvoidanceGuidanceLawforSub-Interceptors

    InterceptingUAVCluster

    LuoRuining1*,HuangShucai1,ZhaoYan1,ZhangZhen1,ZhangFei2

    (1.AirandMissileDefenseCollege,AirForceEngineeringUniversity,Xian710051,China;

    2.Unit93942ofPLA,Xianyang713800,China)

    Abstract:Inordertosolvethepossiblecollisionproblemintheprocessofsub-interceptorsinterceptingUAVcluster,thispaperproposesacollisionavoidancepureproportionalguidancelaw(CA-PPN)bycombiningthevirtualrepulsionforcewiththepureproportionalguidancelaw.Firstly,itproposesanoperationalconceptofsubmunitionmissileinterceptingUAVcluster,andanalyzestheguidanceproblemofsub-interceptors.Then,aimingatthepossiblecollisionproblemintheguidanceofsub-interceptors,itdesignsanartificialpotentialfieldforcollisionavoidancebetweensub-interceptors,andanalyzesthecharacteristicsofPPN.Onthisbasis,combiningPPNwiththevirtualrepulsionforceintheartificialpotentialfield,andaddingthefieldofviewangleandoverloadconstraints,theCA-PPNisproposed.ThesimulationresultsshowthatCA-PPNcaneffectivelyinterceptthetargetandavoidthecollisionbetweensub-interceptors.

    Keywords:collisionavoidanceguidancelaw;virtualrepulsionforce;pureproportionalnavigation;artificialpotentialfield;sub-interceptor;missile;UAVcluster

    收稿日期:2022-06-01

    基金項目:國家自然科學(xué)基金項目(61703424)

    *作者簡介:羅瑞寧(1996-),男,河北衡水人,碩士研究生。

    精品不卡国产一区二区三区| 国产探花极品一区二区| 又爽又黄无遮挡网站| 国产av不卡久久| 18+在线观看网站| 国产熟女欧美一区二区| av在线天堂中文字幕| 久久草成人影院| 91在线精品国自产拍蜜月| 九九热线精品视视频播放| 亚洲精品一区av在线观看| 婷婷精品国产亚洲av在线| 欧美+亚洲+日韩+国产| 国产爱豆传媒在线观看| 亚洲经典国产精华液单| 伦理电影大哥的女人| 一个人看的www免费观看视频| 又黄又爽又免费观看的视频| 少妇熟女aⅴ在线视频| 国产欧美日韩精品一区二区| 啦啦啦观看免费观看视频高清| 成年免费大片在线观看| 麻豆成人午夜福利视频| 亚洲精品影视一区二区三区av| 欧美黑人巨大hd| 日本爱情动作片www.在线观看 | 久久久久国产精品人妻aⅴ院| 亚洲av中文字字幕乱码综合| 国内精品久久久久精免费| 国产69精品久久久久777片| 97超级碰碰碰精品色视频在线观看| 国产精品嫩草影院av在线观看 | 91麻豆av在线| 狂野欧美白嫩少妇大欣赏| 成人av在线播放网站| 午夜免费激情av| 99riav亚洲国产免费| 免费大片18禁| 精品国产三级普通话版| eeuss影院久久| 91av网一区二区| 国产精品一区二区三区四区免费观看 | 欧美区成人在线视频| 欧美日韩综合久久久久久 | 午夜福利欧美成人| 亚洲欧美日韩无卡精品| 久久精品91蜜桃| 日本成人三级电影网站| 日韩国内少妇激情av| 国产精品99久久久久久久久| 欧美精品啪啪一区二区三区| 网址你懂的国产日韩在线| 国内精品一区二区在线观看| 免费av观看视频| 搡老熟女国产l中国老女人| 午夜日韩欧美国产| 日日摸夜夜添夜夜添小说| 欧美最黄视频在线播放免费| 久久久久精品国产欧美久久久| 免费黄网站久久成人精品| 亚洲图色成人| 欧洲精品卡2卡3卡4卡5卡区| 99热精品在线国产| 国产伦一二天堂av在线观看| 国产在线精品亚洲第一网站| 国产精品一区二区免费欧美| 99久久成人亚洲精品观看| 成熟少妇高潮喷水视频| 国产v大片淫在线免费观看| 在线免费观看不下载黄p国产 | 变态另类成人亚洲欧美熟女| 欧美日韩综合久久久久久 | 男女下面进入的视频免费午夜| 国产爱豆传媒在线观看| 永久网站在线| 美女被艹到高潮喷水动态| 久久国内精品自在自线图片| 成年人黄色毛片网站| 好男人在线观看高清免费视频| 一夜夜www| 日韩av在线大香蕉| 亚洲av一区综合| 久久午夜福利片| 国产精品久久视频播放| 成人无遮挡网站| www.www免费av| 夜夜爽天天搞| 亚洲精华国产精华精| 一边摸一边抽搐一进一小说| 亚洲五月天丁香| 国产亚洲av嫩草精品影院| 国产黄色小视频在线观看| 极品教师在线免费播放| 极品教师在线免费播放| 国产高清有码在线观看视频| 欧美zozozo另类| 亚洲av第一区精品v没综合| 男女之事视频高清在线观看| 日日干狠狠操夜夜爽| 黄色视频,在线免费观看| 18禁在线播放成人免费| 黄色视频,在线免费观看| 欧美日韩乱码在线| 亚洲aⅴ乱码一区二区在线播放| 久久九九热精品免费| 色播亚洲综合网| 精品国产三级普通话版| 久久精品久久久久久噜噜老黄 | 色哟哟哟哟哟哟| 亚洲欧美日韩东京热| 国产精品1区2区在线观看.| 亚洲av二区三区四区| 国产白丝娇喘喷水9色精品| 搞女人的毛片| 亚洲午夜理论影院| 搞女人的毛片| 亚洲国产精品合色在线| 老司机深夜福利视频在线观看| 啦啦啦韩国在线观看视频| 成人综合一区亚洲| 国产色婷婷99| 成人一区二区视频在线观看| 又黄又爽又刺激的免费视频.| 别揉我奶头 嗯啊视频| 成人精品一区二区免费| 欧美最新免费一区二区三区| 亚洲成人免费电影在线观看| 欧美一区二区精品小视频在线| av.在线天堂| 国产三级中文精品| 久久久久免费精品人妻一区二区| av黄色大香蕉| 午夜激情欧美在线| 欧美3d第一页| 欧美另类亚洲清纯唯美| 国产精品久久久久久久久免| 一卡2卡三卡四卡精品乱码亚洲| 麻豆精品久久久久久蜜桃| 女人被狂操c到高潮| 成人亚洲精品av一区二区| 日韩欧美在线二视频| 国产久久久一区二区三区| 观看免费一级毛片| 日本一本二区三区精品| 久久国产精品人妻蜜桃| 天堂网av新在线| 欧美激情久久久久久爽电影| 又黄又爽又免费观看的视频| 久久香蕉精品热| av中文乱码字幕在线| avwww免费| 久久精品国产自在天天线| 午夜福利欧美成人| 天堂√8在线中文| 听说在线观看完整版免费高清| 久久久久久国产a免费观看| 综合色av麻豆| 中亚洲国语对白在线视频| 午夜激情福利司机影院| 国产精品伦人一区二区| 男女那种视频在线观看| 国产精品野战在线观看| 伦精品一区二区三区| 日本黄色视频三级网站网址| 两个人的视频大全免费| 又粗又爽又猛毛片免费看| 成人三级黄色视频| 亚洲自偷自拍三级| 精品免费久久久久久久清纯| 日韩欧美精品免费久久| 国产日本99.免费观看| 老司机福利观看| 啦啦啦观看免费观看视频高清| 不卡视频在线观看欧美| 国产极品精品免费视频能看的| 亚洲一级一片aⅴ在线观看| 国产精品久久久久久精品电影| 日韩欧美三级三区| 欧美性感艳星| 一本一本综合久久| 亚洲成a人片在线一区二区| 久久国产精品人妻蜜桃| 亚洲,欧美,日韩| 日韩中文字幕欧美一区二区| 亚洲精品在线观看二区| 国产精品福利在线免费观看| 久久精品夜夜夜夜夜久久蜜豆| 久久精品综合一区二区三区| 色综合站精品国产| 在线观看美女被高潮喷水网站| 99久久精品国产国产毛片| aaaaa片日本免费| 精品国内亚洲2022精品成人| 国产精品人妻久久久久久| 成人永久免费在线观看视频| 三级男女做爰猛烈吃奶摸视频| 国产极品精品免费视频能看的| 91av网一区二区| 综合色av麻豆| 亚洲最大成人手机在线| 特大巨黑吊av在线直播| 国产精品久久电影中文字幕| 日韩欧美在线二视频| 欧美潮喷喷水| 三级男女做爰猛烈吃奶摸视频| 网址你懂的国产日韩在线| 免费在线观看影片大全网站| 色吧在线观看| 精品国产三级普通话版| 色综合站精品国产| 丝袜美腿在线中文| 欧美xxxx性猛交bbbb| 91麻豆av在线| 波多野结衣巨乳人妻| 欧美精品啪啪一区二区三区| 精品人妻视频免费看| 麻豆国产97在线/欧美| 亚洲国产精品成人综合色| 国产精品嫩草影院av在线观看 | 中国美白少妇内射xxxbb| 又紧又爽又黄一区二区| 色在线成人网| 国产精品三级大全| www.www免费av| 简卡轻食公司| 国产成人a区在线观看| 欧美日韩综合久久久久久 | 国产一区二区在线av高清观看| 天美传媒精品一区二区| 非洲黑人性xxxx精品又粗又长| 高清在线国产一区| 欧美绝顶高潮抽搐喷水| 中文字幕av成人在线电影| 欧美精品国产亚洲| 97热精品久久久久久| 精品久久久久久成人av| 精品久久久久久久人妻蜜臀av| 韩国av一区二区三区四区| 成人毛片a级毛片在线播放| 久久久成人免费电影| 99热只有精品国产| 亚洲va在线va天堂va国产| 日韩欧美在线二视频| 中文字幕高清在线视频| 国产一区二区三区视频了| 成人国产麻豆网| 亚洲国产精品成人综合色| 国产人妻一区二区三区在| 老师上课跳d突然被开到最大视频| 国产免费男女视频| 深夜a级毛片| 99九九线精品视频在线观看视频| 好男人在线观看高清免费视频| 麻豆精品久久久久久蜜桃| 大又大粗又爽又黄少妇毛片口| 免费无遮挡裸体视频| 国内精品美女久久久久久| 黄色女人牲交| а√天堂www在线а√下载| 国产淫片久久久久久久久| 久久人妻av系列| 日本黄色视频三级网站网址| 一区福利在线观看| 超碰av人人做人人爽久久| x7x7x7水蜜桃| 欧美最黄视频在线播放免费| 一级黄色大片毛片| 男女下面进入的视频免费午夜| 男女之事视频高清在线观看| 国产成人av教育| 成人欧美大片| 国产精品1区2区在线观看.| 成人国产一区最新在线观看| 99热只有精品国产| 亚洲欧美清纯卡通| 十八禁国产超污无遮挡网站| 波野结衣二区三区在线| 他把我摸到了高潮在线观看| 色精品久久人妻99蜜桃| 国产精品久久电影中文字幕| 日韩在线高清观看一区二区三区 | 无遮挡黄片免费观看| 国产伦人伦偷精品视频| 免费看日本二区| 亚洲va日本ⅴa欧美va伊人久久| av中文乱码字幕在线| 成人特级av手机在线观看| 别揉我奶头~嗯~啊~动态视频| 亚洲在线自拍视频| 一级黄色大片毛片| 欧美日韩亚洲国产一区二区在线观看| 亚洲欧美精品综合久久99| 中文字幕久久专区| 免费av不卡在线播放| 日本 av在线| 99九九线精品视频在线观看视频| 国产美女午夜福利| 国产激情偷乱视频一区二区| 一区二区三区激情视频| 日韩欧美三级三区| 日韩亚洲欧美综合| 黄色欧美视频在线观看| 国产午夜福利久久久久久| 午夜a级毛片| 男女视频在线观看网站免费| 中文字幕免费在线视频6| 人妻制服诱惑在线中文字幕| 91av网一区二区| 一区二区三区四区激情视频 | 亚洲av中文av极速乱 | 国产精品久久电影中文字幕| 国产探花极品一区二区| 全区人妻精品视频| 22中文网久久字幕| 亚洲自偷自拍三级| 国产精品女同一区二区软件 | 此物有八面人人有两片| 国产男人的电影天堂91| 1024手机看黄色片| 亚洲精品乱码久久久v下载方式| .国产精品久久| 精品国产三级普通话版| 久久久久九九精品影院| 亚洲第一电影网av| 亚洲第一区二区三区不卡| 五月玫瑰六月丁香| 嫁个100分男人电影在线观看| 人妻丰满熟妇av一区二区三区| 夜夜看夜夜爽夜夜摸| 亚洲自偷自拍三级| 日本撒尿小便嘘嘘汇集6| 少妇高潮的动态图| 国产精品久久久久久精品电影| 日韩高清综合在线| 亚洲在线观看片| ponron亚洲| 嫩草影院精品99| x7x7x7水蜜桃| 又黄又爽又免费观看的视频| 久久精品夜夜夜夜夜久久蜜豆| 亚洲va在线va天堂va国产| 人妻丰满熟妇av一区二区三区| 国产高清视频在线播放一区| 狂野欧美激情性xxxx在线观看| 亚洲一区高清亚洲精品| 国产毛片a区久久久久| 久久99热这里只有精品18| 久久午夜福利片| 99久久精品一区二区三区| 欧美日本亚洲视频在线播放| 超碰av人人做人人爽久久| 别揉我奶头 嗯啊视频| 一进一出抽搐gif免费好疼| 色哟哟·www| 成年女人毛片免费观看观看9| 免费看光身美女| 久久久久久久久大av| 可以在线观看的亚洲视频| 中文字幕高清在线视频| 色在线成人网| 亚洲精品影视一区二区三区av| 九色国产91popny在线| 日本黄色视频三级网站网址| 日韩欧美三级三区| 午夜激情福利司机影院| 国产亚洲精品av在线| 亚洲精品成人久久久久久| 97热精品久久久久久| 麻豆一二三区av精品| 亚洲精品亚洲一区二区| 日本爱情动作片www.在线观看 | 亚洲美女搞黄在线观看 | 亚洲精品一区av在线观看| 欧美日韩乱码在线| 最近最新中文字幕大全电影3| 日本成人三级电影网站| 高清日韩中文字幕在线| 国产中年淑女户外野战色| 国产白丝娇喘喷水9色精品| 美女高潮喷水抽搐中文字幕| 两人在一起打扑克的视频| 三级男女做爰猛烈吃奶摸视频| 波多野结衣高清作品| 国产成人av教育| 老熟妇乱子伦视频在线观看| 日韩欧美三级三区| 免费在线观看成人毛片| 国产黄片美女视频| 国产精品国产三级国产av玫瑰| 国产视频内射| 亚洲av五月六月丁香网| 在线免费观看的www视频| x7x7x7水蜜桃| 黄色日韩在线| 欧美精品国产亚洲| 无人区码免费观看不卡| 波野结衣二区三区在线| 亚洲精品粉嫩美女一区| 欧美极品一区二区三区四区| 色综合色国产| 亚洲熟妇中文字幕五十中出| 久久精品国产亚洲网站| 黄色日韩在线| 女生性感内裤真人,穿戴方法视频| 欧美xxxx性猛交bbbb| 在线a可以看的网站| 国产一区二区激情短视频| 校园春色视频在线观看| 免费在线观看影片大全网站| 日本爱情动作片www.在线观看 | 又爽又黄a免费视频| 精品人妻1区二区| 一区二区三区免费毛片| 国产亚洲精品av在线| 一区福利在线观看| 精品一区二区三区视频在线| 亚洲国产精品合色在线| 国产黄a三级三级三级人| 麻豆av噜噜一区二区三区| 永久网站在线| 亚洲七黄色美女视频| 免费在线观看日本一区| 国产一区二区三区av在线 | 99热这里只有精品一区| 91久久精品电影网| 午夜福利在线观看吧| 在线播放无遮挡| 国产私拍福利视频在线观看| 国产真实乱freesex| 在线播放无遮挡| 国产男人的电影天堂91| 一级av片app| 成人毛片a级毛片在线播放| 99久久无色码亚洲精品果冻| 国产精品一区二区免费欧美| 97人妻精品一区二区三区麻豆| 一本久久中文字幕| 看黄色毛片网站| 国产精品国产三级国产av玫瑰| 欧美丝袜亚洲另类 | 国产精品自产拍在线观看55亚洲| 精品日产1卡2卡| 麻豆久久精品国产亚洲av| 免费大片18禁| 亚洲欧美日韩高清在线视频| 亚洲精品成人久久久久久| 别揉我奶头~嗯~啊~动态视频| 婷婷亚洲欧美| 天堂√8在线中文| 日本在线视频免费播放| 国产乱人视频| 午夜老司机福利剧场| 久久久久久久久久成人| 亚洲熟妇中文字幕五十中出| 久久精品国产亚洲网站| 国产午夜福利久久久久久| 国产熟女欧美一区二区| 国产高清视频在线观看网站| 欧美日韩亚洲国产一区二区在线观看| 又粗又爽又猛毛片免费看| 久久久久久久精品吃奶| 免费搜索国产男女视频| 久久九九热精品免费| 亚洲国产精品合色在线| or卡值多少钱| 中文字幕久久专区| 男女下面进入的视频免费午夜| 成人性生交大片免费视频hd| 婷婷精品国产亚洲av| 18禁在线播放成人免费| 真人做人爱边吃奶动态| 国产欧美日韩精品亚洲av| 国产精品永久免费网站| 亚洲乱码一区二区免费版| 久久午夜亚洲精品久久| 黄色配什么色好看| 舔av片在线| 国产大屁股一区二区在线视频| 国内揄拍国产精品人妻在线| 99精品久久久久人妻精品| 久久香蕉精品热| 日本免费a在线| 男女那种视频在线观看| 亚洲国产精品sss在线观看| 麻豆国产av国片精品| 小说图片视频综合网站| 国内毛片毛片毛片毛片毛片| 中文字幕久久专区| 男人舔女人下体高潮全视频| 夜夜夜夜夜久久久久| 国产精品久久久久久亚洲av鲁大| 男人的好看免费观看在线视频| 两性午夜刺激爽爽歪歪视频在线观看| 久久人人爽人人爽人人片va| 又黄又爽又刺激的免费视频.| 成人国产麻豆网| 久久久国产成人精品二区| 国产女主播在线喷水免费视频网站 | 性色avwww在线观看| 国产精品1区2区在线观看.| 日韩欧美免费精品| 国产av麻豆久久久久久久| 丰满的人妻完整版| 校园春色视频在线观看| 亚洲精品日韩av片在线观看| 美女高潮的动态| 久久欧美精品欧美久久欧美| 桃红色精品国产亚洲av| 又黄又爽又免费观看的视频| 欧美一区二区亚洲| 国产精品不卡视频一区二区| or卡值多少钱| 淫妇啪啪啪对白视频| 国产精品一及| 欧美国产日韩亚洲一区| 无遮挡黄片免费观看| 麻豆av噜噜一区二区三区| 2021天堂中文幕一二区在线观| 在现免费观看毛片| 99视频精品全部免费 在线| 久久午夜福利片| 麻豆成人av在线观看| 一个人观看的视频www高清免费观看| 亚洲,欧美,日韩| 制服丝袜大香蕉在线| 中文字幕高清在线视频| 男女边吃奶边做爰视频| 日日摸夜夜添夜夜添小说| 欧美精品国产亚洲| 久久久久久久亚洲中文字幕| 香蕉av资源在线| 免费黄网站久久成人精品| 国产不卡一卡二| 午夜福利成人在线免费观看| 人人妻,人人澡人人爽秒播| 无人区码免费观看不卡| 国产三级中文精品| 国产精品一区www在线观看 | 欧美成人一区二区免费高清观看| 亚洲va日本ⅴa欧美va伊人久久| 国产高清有码在线观看视频| 欧美日韩瑟瑟在线播放| 免费搜索国产男女视频| 久久久久久久久久成人| 午夜亚洲福利在线播放| 两性午夜刺激爽爽歪歪视频在线观看| 我要看日韩黄色一级片| 人人妻人人看人人澡| 天堂√8在线中文| av专区在线播放| 亚洲av日韩精品久久久久久密| 小蜜桃在线观看免费完整版高清| 丰满乱子伦码专区| 日日摸夜夜添夜夜添小说| 国产成人福利小说| 麻豆一二三区av精品| 免费高清视频大片| 国产又黄又爽又无遮挡在线| 久久国产乱子免费精品| 成熟少妇高潮喷水视频| 久久亚洲精品不卡| 国产一区二区三区在线臀色熟女| 精品一区二区三区av网在线观看| 国产av不卡久久| 国产色爽女视频免费观看| 免费电影在线观看免费观看| 精品久久久久久成人av| 色在线成人网| 成人二区视频| 美女大奶头视频| 亚洲一区高清亚洲精品| 亚洲精品一区av在线观看| 在线看三级毛片| 此物有八面人人有两片| 能在线免费观看的黄片| 久久欧美精品欧美久久欧美| 少妇人妻精品综合一区二区 | 国产国拍精品亚洲av在线观看| 色av中文字幕| 女同久久另类99精品国产91| 久久中文看片网| 国产成人aa在线观看| 免费人成在线观看视频色| 免费看a级黄色片| 亚洲国产精品成人综合色| 极品教师在线视频| av.在线天堂| 欧美日韩乱码在线| 欧美高清性xxxxhd video| 亚洲精品日韩av片在线观看| 国产精品久久久久久久久免| 日韩中字成人| 高清在线国产一区| 亚洲av免费在线观看| 成人综合一区亚洲| 91久久精品国产一区二区三区| 日韩一区二区视频免费看| 淫秽高清视频在线观看| 国产欧美日韩精品一区二区| 久久精品综合一区二区三区| 免费一级毛片在线播放高清视频| 十八禁网站免费在线| 亚洲中文日韩欧美视频| 我要搜黄色片| 亚洲无线观看免费| 日韩av在线大香蕉| 亚洲av二区三区四区| 桃红色精品国产亚洲av| 1000部很黄的大片| 久久久久性生活片| 国产一区二区三区av在线 | 少妇被粗大猛烈的视频| 亚洲图色成人| 深夜精品福利| 国产精品人妻久久久久久|