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

    面向臨近空間目標(biāo)攔截的預(yù)測命中點設(shè)計方法

    2020-10-12 06:27:10尹中杰朱柏羊韋文書安帥斌
    宇航學(xué)報 2020年9期
    關(guān)鍵詞:航程超聲速彈道

    尹中杰,劉 凱,朱柏羊,韋文書,安帥斌

    (1. 大連理工大學(xué)航空航天學(xué)院,大連 116024;2. 遼寧省空天飛行器前沿技術(shù)重點實驗室,大連 116024; 3. 北京空天技術(shù)研究所,北京 100074;4.中國運載火箭技術(shù)研究院,北京 100076)

    0 引 言

    近年來,隨著高超聲速技術(shù)的發(fā)展,以美俄為代表的軍事大國正在大力發(fā)展臨近空間高超聲速武器[1]。由于高超武器具有飛行速度快、機(jī)動能力強(qiáng)、穩(wěn)定跟蹤難度大的特點,導(dǎo)致防空反導(dǎo)作戰(zhàn)面臨巨大挑戰(zhàn)。因此,以美國為代表的世界大國開始聚焦高超聲速武器防御技術(shù)。2018年美國導(dǎo)彈防御局(MDA)授出21份“高超聲速防御武器系統(tǒng)概念”研究合同用于提高臨近空間防御能力[2],其中較為引人注意的是3個基于高超聲速平臺的臨近空間目標(biāo)攔截方案,分別是波音公司的“針對高超聲速武器的超高聲速攔截器概念(HYVINT)”、通用原子電磁系統(tǒng)公司的“超高速攔截器系統(tǒng)”、德雷伯實驗室的“超高速吸氣者(HSAB)”。由于高超聲速平臺具有快速響應(yīng)、射程遠(yuǎn)以及中段機(jī)動修正能力強(qiáng)的優(yōu)勢,可有效解決傳統(tǒng)火箭動力攔截系統(tǒng)反臨近空間高超聲速目標(biāo)的瓶頸問題,是后續(xù)發(fā)展中具有潛力的一類攔截方案。

    在臨近空間目標(biāo)攔截方案中,預(yù)測命中點設(shè)計問題由于因為存在可行解搜索耗時長、攔截誤差大的問題被廣泛的研究。為進(jìn)一步發(fā)揮高超攔截平臺的射程遠(yuǎn)的優(yōu)勢,需要針對作戰(zhàn)需求基于攔截平臺的動力學(xué)特性重新設(shè)計預(yù)測命中點規(guī)劃方法。

    預(yù)測命中點是攻擊彈、攔截彈同時間到達(dá)的一點,預(yù)測命中規(guī)劃需要在攔截彈發(fā)射前,通過算法計算出合理的攔截點,并為攔截彈裝訂相應(yīng)的諸元,使得攔截彈準(zhǔn)時精確的飛抵預(yù)測命中點。在該方面,國內(nèi)外學(xué)者普遍針對火箭動力攔截彈設(shè)計基于標(biāo)準(zhǔn)彈道族的預(yù)測命中點規(guī)劃方法。

    基于標(biāo)準(zhǔn)彈道族的預(yù)測命中點設(shè)計是指攔截方離線獲得攔截覆蓋區(qū)域得到標(biāo)準(zhǔn)彈道族,并通過調(diào)整攔截彈發(fā)射方位角形成等時間攔截曲面,并結(jié)合目標(biāo)預(yù)報彈道篩選時間相同、位置重合相遇點的方法。該計算方法可以針對不同類型的目標(biāo)飛行器快速有效地尋找準(zhǔn)確的預(yù)測命中點。

    在標(biāo)準(zhǔn)彈道族設(shè)計方面,文獻(xiàn)[3-4]通過設(shè)計攻角指令設(shè)計完成了基于火箭動力的拋物線、直線標(biāo)準(zhǔn)彈道族,上述兩種彈道設(shè)計方式具有方法簡單、可行性強(qiáng)的特點。但是由于彈道的物理特性的因素,上述兩種彈道并不能很完整的囊括攔截平臺的所有攔截覆蓋范圍,存在覆蓋區(qū)域有缺失的問題。文獻(xiàn)[5]首先對攔截需求進(jìn)行分析,引入性能指標(biāo),優(yōu)化出指標(biāo)最優(yōu)攔截彈道,組成針對不同優(yōu)化指標(biāo)的多個標(biāo)準(zhǔn)彈道族。這種設(shè)計方法可以針對不同攔截狀況,針對攔截需求選取相對應(yīng)的標(biāo)準(zhǔn)彈道族進(jìn)行更有針對性的預(yù)測命中點搜索。同時由于該方法有多個標(biāo)準(zhǔn)彈道族,可以在某一彈道族搜索失敗后,切換至另一彈道族進(jìn)行搜索,從而提高攔截成功率。并且通過彈道族疊加,在攔截覆蓋范圍上相較于文獻(xiàn)[3-4]具有一定優(yōu)勢,但是通過彈道優(yōu)化算法并不能優(yōu)化出攔截空域最大的標(biāo)準(zhǔn)彈道族,所以同樣無法克服無法囊括所有攔截覆蓋范圍的問題。

    在搜索算法方面,文獻(xiàn)[6]采用迭代發(fā)射諸元的方式解算預(yù)測命中點。通過假設(shè)初始的攻角、攔截時間獲得彈道,通過位置偏差迭代攻角、時間修正量,從而逼近發(fā)射諸元正確值。這種方法采用在線求解彈道并修正誤差的方式,具有攔截誤差小的優(yōu)勢,但因為解算諸元較多存在遺漏可行解的問題。文獻(xiàn)[7-8]針對彈道導(dǎo)彈被動段攔截,假設(shè)彈目距離是單調(diào)遞減的,基于斜距迭代的方法,獲得相應(yīng)的發(fā)射諸元從而獲得攔截彈彈道。但是臨近空間高超目標(biāo)相較于彈道導(dǎo)彈具有極強(qiáng)的機(jī)動能力,能夠進(jìn)行大范圍的機(jī)動,攔截陣地與目標(biāo)的斜距可能并不單調(diào)。預(yù)報彈道上可能存在有多個斜距相同的點,會導(dǎo)致迭代出的發(fā)射諸元存在偏差,致使攔截失敗。同時由于搜索算法并未考慮標(biāo)準(zhǔn)彈道族內(nèi)各彈道之間存在的時間、空間規(guī)律,沒有設(shè)置對應(yīng)的彈道搜索區(qū)間篩選、導(dǎo)致需要對標(biāo)準(zhǔn)彈道族內(nèi)的彈道逐一搜索,導(dǎo)致計算耗時長,方法存在局限性。

    針對以上問題,本文瞄準(zhǔn)基于吸氣式高超聲速平臺的臨近空間目標(biāo)攔截概念方案,開展預(yù)測命中點設(shè)計問題研究,基于攔截空域最大原則,結(jié)合吸氣式動力攔截平臺的動力學(xué)特性構(gòu)建標(biāo)準(zhǔn)彈道族,在對目標(biāo)預(yù)報彈道搜索區(qū)間預(yù)處理的基礎(chǔ)上,給出高超聲速攔截平臺的預(yù)測命中點高效搜索算法,快速完成預(yù)測命中點規(guī)劃。

    1 方案構(gòu)想與吸氣式攔截平臺質(zhì)點動力學(xué)建模

    1.1 方案構(gòu)想

    新式反臨近空間目標(biāo)高超聲速導(dǎo)彈由助推器、攔截平臺、動能攔截器組成。助推器用于爬升段的助推工作,由固體火箭發(fā)動機(jī)提供動力。攔截平臺為吸氣式高超聲速飛行器,其任務(wù)是搭載一枚或多枚動能攔截器,并將其運抵至導(dǎo)引頭開機(jī)距離。攔截平臺采用軸對稱升力體構(gòu)型,動力方面采用亞燃/超燃雙模態(tài)沖壓發(fā)動機(jī),動能攔截器采用紅外式側(cè)窗導(dǎo)引頭,具有不小于50 km的探測距離。

    在裝備攔截平臺后,防御系統(tǒng)將大體按照如下的方式實施防御作戰(zhàn):

    1)敵方的高超聲速武器放射后,通過天基預(yù)警衛(wèi)星探測目標(biāo),將預(yù)警信息傳送給作戰(zhàn)指揮系統(tǒng)。作戰(zhàn)指揮系統(tǒng)將對目標(biāo)進(jìn)行威脅評估、資源分配,并引導(dǎo)地基雷達(dá)搜索、捕獲和跟蹤目標(biāo)。作戰(zhàn)指揮系統(tǒng)制定交戰(zhàn)計劃,向作戰(zhàn)平臺發(fā)送目標(biāo)數(shù)據(jù)、彈道信息和發(fā)射指令。作戰(zhàn)平臺按照初始發(fā)射方位需求調(diào)整平臺姿態(tài),然后實施發(fā)射高超聲速導(dǎo)彈指令,發(fā)射高超聲速導(dǎo)彈。

    2)高超聲速導(dǎo)彈發(fā)射后,首先由助推火箭加速,助推到Ma4以上后,助推器分離,攔截平臺啟動吸氣式動力發(fā)動機(jī)加速飛行至指定巡航高度/速度,然后依靠吸氣式發(fā)動機(jī)高比沖優(yōu)勢的高超聲速巡航飛行。飛行過程中,地面雷達(dá)站仍然實時跟蹤來襲目標(biāo),并據(jù)攔截平臺和來襲目標(biāo)信息實時解算新的預(yù)測命中點信息和預(yù)設(shè)彈道信息并實時發(fā)送給攔截平臺,攔截平臺按照接收信息對飛行彈道進(jìn)行修正,導(dǎo)引向新的預(yù)測命中點。

    3)當(dāng)攔截平臺飛到預(yù)測命中點附近時,動能攔截器與攔截平臺分離,攔截器立即用紅外感應(yīng)、跟蹤、識別目標(biāo),確定瞄準(zhǔn)點;在姿軌控發(fā)動機(jī)的控制下,自行接近目標(biāo),最后通過直接碰撞攔截并摧毀目標(biāo)。

    圖1 攔截流程示意圖Fig.1 Operational flow diagram

    由于吸氣式攔截平臺所采用的超燃沖壓發(fā)動機(jī)與固體火箭發(fā)動機(jī)存在很大的差異,而現(xiàn)有的標(biāo)準(zhǔn)彈道族設(shè)計方法并不能有效利用超燃沖壓發(fā)動機(jī)工作特性,導(dǎo)致其無法為吸氣式攔截平臺設(shè)計出攔截空域最大的標(biāo)準(zhǔn)彈道族。除此以外,臨近空間高機(jī)動目標(biāo)相較于彈道導(dǎo)彈具有機(jī)動能力強(qiáng)的特點,傳統(tǒng)的預(yù)測命中點解算流程在攔截臨近空間高機(jī)動目標(biāo)時存在解算速度慢,預(yù)測命中點實際攔截誤差較大的問題?;诖?,這里針對攔截平臺動力學(xué)特性重新設(shè)計了標(biāo)準(zhǔn)彈道族與預(yù)測命中點規(guī)劃方法,從而為基于吸氣式高超攔截平臺的預(yù)臨近空間高機(jī)動目標(biāo)攔截任務(wù)提供相關(guān)技術(shù)途徑。

    1.2 攔截平臺質(zhì)點動力學(xué)模型

    模型在球形大地假設(shè)下,不考慮地球自轉(zhuǎn)、在發(fā)射坐標(biāo)系下建立質(zhì)點動力學(xué)模型,根據(jù)攔截平臺的特性完成坐標(biāo)系的轉(zhuǎn)化、受力分析并建立動力學(xué)方程[9]:

    (1)

    式中:x,y,z是發(fā)射系ox,oy,oz軸位置坐標(biāo);Vx,Vy,Vz發(fā)射系ox,oy,oz軸速度;θ,σ是發(fā)射系下彈道傾角、彈道偏角;D,L,Z,T分別是阻力、升力、側(cè)向力、推力;Isp是比沖;Mfv是發(fā)射系與彈體坐標(biāo)系的轉(zhuǎn)換矩陣;Mfd是發(fā)射系與彈道坐標(biāo)系的轉(zhuǎn)換矩陣;h,r分別是高度、地心距;m為攔截彈質(zhì)量;g為地表重力加速度;Re是地球半徑。

    轉(zhuǎn)換矩陣具體形式如下:

    (2)

    (3)

    式中:υ為傾側(cè)角,φ為俯仰角,ψ為偏航角,γ為滾轉(zhuǎn)角。

    1.3 預(yù)測命中點規(guī)劃問題描述

    預(yù)測命中點是預(yù)先計算來襲目標(biāo)和攔截彈的飛行彈道,確定出瞬時遭遇點,即預(yù)測命中點,當(dāng)預(yù)測命中點確定后,就導(dǎo)引攔截彈向其接近。基于此,預(yù)測命中點規(guī)劃問題可以描述為:結(jié)合預(yù)報彈道與標(biāo)準(zhǔn)彈道族,通過搜索算法尋找滿足時間約束、位置約束的攔截點。并為攔截彈裝訂相應(yīng)的諸元,確保其準(zhǔn)時精確到達(dá)攔截點。

    約束條件如下:

    (4)

    式中:tf為預(yù)測飛行時間,xtf,ytf,ztf為預(yù)測命中點在發(fā)射坐標(biāo)系下的位置。

    2 基于吸氣式攔截平臺的預(yù)測命中點計算方法

    2.1 吸氣式攔截平臺飛行策略設(shè)計

    攔截平臺發(fā)射后,共經(jīng)歷了固體火箭助推段與超燃沖壓動力巡航段兩種飛行狀態(tài),這里重點研究巡航階段的飛行軌跡策略。

    助推段設(shè)計采用單級運載火箭上升段指令生成規(guī)律進(jìn)行設(shè)計。該方法將火箭上升過程分為若干飛行階段,每段設(shè)定固定的程序角指令,這些指令由有限的參數(shù)決定,根據(jù)任務(wù)需求離線進(jìn)行優(yōu)化得到這些參數(shù)及其對應(yīng)程序指令,裝訂在助推器系統(tǒng)中,上升段飛行過程中利用姿態(tài)控制系統(tǒng)實現(xiàn)程序指令即可。為尋求彈道軌跡的快速生成,將整個上升段分為3段,包括垂直上升段(0~T1),負(fù)攻角轉(zhuǎn)彎段(T1~T2),重力轉(zhuǎn)彎段(T2~T3)。

    攔截彈全程飛行攻角近似公式為:

    (5)

    (6)

    (7)

    式中: 1/v0為推重比,αm是最大負(fù)攻角幅值,tm為負(fù)攻角值最大時刻,T2為負(fù)攻角轉(zhuǎn)彎段,可以自由設(shè)置;T3為固體火箭燃料耗盡時間。

    巡航飛行初始狀態(tài)為助推段的末狀態(tài),在巡航段動力由超燃沖壓發(fā)動機(jī)提供,推力通過數(shù)據(jù)插值得到,即T=T(q,kr,Ma),其中,q為動壓、kr為推力調(diào)節(jié)閥門開度、Ma為馬赫數(shù)。巡航段考慮攔截平臺的飛行特性,要求攔截平臺盡可能處于等高等速的飛行狀態(tài)。為滿足攔截需求,可以基于飛行高度將彈道分為平飛段與爬升段。

    由于巡航速度大于助推段結(jié)束時的速度,這要求攔截平臺首先要加速至規(guī)定速度,再轉(zhuǎn)入等速飛行狀態(tài)。在加速飛行過程推力為超燃沖壓發(fā)動機(jī)所能提供的最大值,即T=Tmax。

    等速飛行狀態(tài),為使速度保持恒定,因此有

    (8)

    1)平飛段

    h=r-Re

    (9)

    可以對x,y,z求導(dǎo)并再對Vx,Vy,Vz求導(dǎo),獲得高度對攻角的二階導(dǎo)數(shù)變化率:

    (10)

    (11)

    (12)

    式中:dk是搜索方向,則令

    (13)

    即:

    (14)

    這里:

    2)爬升段

    爬升段是指攔截平臺爬升至指定高度的一段飛行過程,可以細(xì)分為定攻角爬升段、直線爬升段、平滑過渡段三個部分。

    首先,攔截平臺在定攻角爬升段通過大攻角增加升力改變彈道傾角,從而使得攔截平臺具有爬升能力。當(dāng)彈道傾角達(dá)到預(yù)設(shè)值后,通過直線爬升的方式爬升至一定高度。為使攔截彈道保持直線爬升,需要彈道傾角變化率恒為0°/s,彈道傾角變化率表達(dá)式為:

    (15)

    可以設(shè)置攻角的迭代初值α0,通過牛頓迭代,即:

    (16)

    最后通過平滑過渡段,通過以飛行高度為自變量,彈道傾角為因變量的過渡函數(shù),獲得彈道傾角的程序指令。

    φ(h)=

    (17)

    θd(h)=φ(h)θ1+(1-φ(h))θ0

    (18)

    式中:hmax為攔截高度,hmin為指定的平滑段起始高度,h為實際高度,θ0為平滑段起始彈道傾角,θ1為期望彈道傾角,由于等高等速的設(shè)計需求,在平滑末端在彈道坐標(biāo)系下期望彈道傾角θ1=0,彈道傾角指令可以簡化為:

    θd(h)=(1-φ(h))θ0

    (19)

    為了將彈道坐標(biāo)系下的彈道傾角指令轉(zhuǎn)化為發(fā)射系下的彈道傾角指令,需要通過坐標(biāo)轉(zhuǎn)換矩陣,通過轉(zhuǎn)化獲得:

    (20)

    (21)

    式中:δd為彈道系下彈道偏角,由于攔截空域是由一束存在于同一平面內(nèi)的標(biāo)準(zhǔn)彈道族通過旋轉(zhuǎn)發(fā)射方位角的方式獲得,所以υ≡0°。獲得在發(fā)射坐標(biāo)系下的彈道傾角θ′后,通過與當(dāng)前實際的發(fā)射坐標(biāo)系彈道傾角θ相減做差,基于迭代步長Δt,可以得到彈道傾角的變化率:

    (22)

    發(fā)射坐標(biāo)系下彈道傾角變化率還可以表示為:

    (23)

    即可通過彈道傾角變化指令完成攻角的迭代。

    2.2 攔截平臺標(biāo)準(zhǔn)彈道族設(shè)計

    高超攔截平臺采用火箭助推垂直發(fā)射,按照相應(yīng)的攻角指令完成助推段飛行。助推段結(jié)束后,通過超燃沖壓發(fā)動機(jī)加速至指定速度并達(dá)到指定高度。

    巡航式攔截器射表由兩種不同的彈道族組成,分別是巡航彈道族與爬升彈道族。

    巡航彈道族是指助推段結(jié)束后,通過超燃沖壓動力直接加速至指定速度,爬升至攔截高度,并在攔截高度上等高等速飛行,并由不同攔截高度的巡航彈道完成巡航彈道族。該彈道族的初始高度變化率為攔截平臺允許的最大高度變化率,所以能在最短的時間到達(dá)指定攔截高度,所以在初期有較大的攔截空域,但是在后期,由于爬升時間過早,導(dǎo)致部分彈道高度過高,比沖較低,存在攔截航程較小的問題,所以單一由巡航彈道所組成的標(biāo)準(zhǔn)彈道族不能很好的覆蓋攔截平臺所有的攔截空域,可以通過引入采用其他飛行策略的標(biāo)準(zhǔn)彈道族組成標(biāo)準(zhǔn)彈道族群,從而彌補(bǔ)巡航彈道族未覆蓋的攔截空域。

    為彌補(bǔ)巡航彈道族在攔截高度較高時攔截航程較小,無法覆蓋全部的攔截區(qū)域的問題,引入了爬升彈道族。爬升彈道在助推段結(jié)束后,首先在燃料最省的基準(zhǔn)彈道飛行,在某一時刻爬升,當(dāng)高度滿足攔截高度再進(jìn)行等高飛行。并由不同爬升時刻、不同攔截高度的爬升彈道組成爬升彈道族。爬升彈道利用基準(zhǔn)彈道比沖大的特點,可以獲得較大的航程。

    由于巡航式高超攔截平臺具有寬速域特性,可以選擇不同的速度進(jìn)行飛行。那么在同一個攔截點,采用不同飛行速度的攔截平臺將有不相同的攔截時間?;诖丝梢栽O(shè)計出多個不同飛行速度的標(biāo)準(zhǔn)彈道族組成標(biāo)準(zhǔn)彈道族群,通過擴(kuò)大搜索范圍的方式提高攔截成功率。

    圖2 攔截示意圖Fig.2 Interception diagram

    2.3 攔截搜索區(qū)間計算

    為提高搜索算法的運算效率,需要基于吸氣式攔截平臺動力學(xué)特性進(jìn)行有針對的減小搜索范圍,同時為保證攔截器處于逆軌攔截狀態(tài),且不存在不同時間、航程相同的問題,需要對目標(biāo)飛行器進(jìn)行攔截區(qū)域劃分,并分段討論。

    首先,獲得航程與時間的關(guān)系,篩選出目標(biāo)—陣地航程單調(diào)遞減區(qū)間,確定一定時間內(nèi)的航程極值點個數(shù)n與相對應(yīng)的時間tn,由此將攔截時間分為n+1段,單調(diào)時間段為t∈[ti-1,ti],i∈[0,n+1],由此確保在規(guī)定時間段內(nèi),航程變化是單調(diào)的,不存在航程相同的可能。

    由于預(yù)測命中點存在時間與空間上的約束,結(jié)合預(yù)測命中點規(guī)劃需求,可以將其轉(zhuǎn)化為帶有航程約束、高度約束的攔截時間與發(fā)射方位角求解。則可以尋找在攔截時間段內(nèi)滿足航程約束的彈道集合,再從彈道集合中篩選出滿足高度約束的彈道,從而根據(jù)滿足航程、高度約束的攔截點位置計算出攔截時間與攔截平臺發(fā)射方位角。

    圖3 時間預(yù)劃分示意圖Fig.3 Time division diagram

    由于攔截平臺動力學(xué)特性與標(biāo)準(zhǔn)彈道族構(gòu)成方式,同一速度下標(biāo)準(zhǔn)彈道族內(nèi)各條彈道,在相同時間航程差較小。而基準(zhǔn)彈道(燃料最省彈道)又具有最遠(yuǎn)的航程,可以覆蓋本彈道族內(nèi)在規(guī)定時間區(qū)間所有彈道的航程區(qū)間。

    所以可以判斷在某一區(qū)間段內(nèi)目標(biāo)航程—時間與基準(zhǔn)彈道航程—時間曲線是否存在交點,若存在交點,可以認(rèn)為在標(biāo)準(zhǔn)彈道族中存在滿足航程約束的多條彈道,具有篩選到滿足高度約束的預(yù)測命中點的可能性。

    交點解算流程如下:

    1)判斷基準(zhǔn)彈道與目標(biāo)彈道在某一時間區(qū)間t∈[ti-1,ti],i∈[0,n+1]有沒有航程重疊的區(qū)域。

    如果:

    R1

    (24)

    則兩個區(qū)間沒有交點,搜索結(jié)束。其中R1指基準(zhǔn)彈道的航程最大值,R2指目標(biāo)的距離陣地航程最小值。

    否則,找出兩者的航程重疊區(qū)域

    tqj=[t2(R1),t2(R2)]∩[t1(R2),t1(R1)]

    (25)

    式中:t2(R)是指目標(biāo)航程為R的時間,t1(R)是指攔截平臺航程為R的時間,tqj為篩選后精簡的時間區(qū)間。

    ΔR=Rgj(tqj)-Rlj(tqj)

    (26)

    式中:Rlj(tqj)指攔截平臺在時間區(qū)間tqj內(nèi)的航程區(qū)間,Rgj(tqj)指目標(biāo)在時間區(qū)間tqj內(nèi)的航程區(qū)間,ΔR為相對距離。

    2)計算交點處飛行時間,令:

    a=sgn(ΔR)′

    (27)

    由于ΔR是彈目的相對距離,sgn(ΔR)=1表示該時刻彈目相對接近,sgn(ΔR)=-1表示該時刻彈目相對遠(yuǎn)離,則在預(yù)測命中點ΔR的符號應(yīng)該發(fā)生切換,sgn(ΔR)導(dǎo)數(shù)的幅值變化±2為即:

    (28)

    (29)

    (30)

    [t0,R]=(A/B)′

    (31)

    式中:t0為交點處的時間,R為交點處的航程。

    通過以交點處的時間t0作為時間迭代的初始值,并通過將t∈[t0-Δtm,t0+Δtm]與單調(diào)區(qū)間取交集的方式獲得搜索區(qū)間,Δtm為可設(shè)置的預(yù)估時間。

    2.4 基于射程迭代的發(fā)射諸元解算

    標(biāo)準(zhǔn)彈道族由兩種彈道組成,其中全部的等高巡航彈道可以看作一束彈道,二次爬升時刻tmin為0 s;而爬升彈道中,相同爬升時刻tmin的彈道亦可以看為一束彈道。

    在一束彈道中,由于爬升時刻相同,速度相同,在相同時間飛行高度最低的彈道航程最遠(yuǎn),飛行高度最高的彈道航程最近。則基于這一特性,由最低端的彈道開始,向最高彈道逐一搜索,以爬升起始時間tminlj,燃料耗盡時間tmaxlj_i為攔截區(qū)間ti∈[tminlj,tmaxlj_i],將攔截區(qū)間ti∈[tminlj,tmaxlj_i]與搜索區(qū)間取交集為t∈[tmin,tmax_i],如果t∈[tmin,tmax_i]為非空集,并且:

    (32)

    式中:Rgj為目標(biāo)與攔截陣地的航程,Rlj為攔截平臺與攔截陣地的航程。通過迭代的方式搜索Rgj=Rlj的時間,即:

    F0=Rgj(tmin)-Rlj(tmin)

    (33)

    F1=Rgj(tmin+Δt)-Rlj(tmin+Δt)

    (34)

    (35)

    滿足|(Rgj(t0)-Rlj(t0))|><ΔL,可獲得攔截時間t0。

    由于目標(biāo)的航程在搜索區(qū)間內(nèi)單調(diào)遞減,而一束彈道內(nèi)相同時間高度越高,航程越短,則在下一搜索過程tmin=t,并重復(fù)上述過程,若:

    Rgj(tmin)

    Rlj(tmax)

    則結(jié)束在這一束彈道的搜索。

    獲得滿足航程約束的攔截彈道后,可基于目標(biāo)預(yù)報彈道中的高度信息篩選出符合攔截要求,并且距離誤差最小的攔截點作為預(yù)測命中點,并計算發(fā)射諸元,發(fā)射諸元包括飛行速度V,發(fā)射方位角Af,攻角α。發(fā)射方位角解算流程如下:

    B=arccos(sinλ0sinφf+cosφ0cosφfcos(λf-λ0))

    (36)

    A1=(sinφf-cosBsinφ0)/(sinBcosφ0)

    (37)

    A2=(sin(λf-λ0)cosφf)/sinB

    (38)

    (39)

    式中:λ0,φ0是攔截陣地經(jīng)緯度,λf,φf是預(yù)測命中點經(jīng)緯度。

    除了發(fā)射方位角Af,為確保攔截平臺能夠滿足約束,還需要基于離線獲得的標(biāo)準(zhǔn)彈道族中,距離預(yù)測命中點最近的攔截彈道獲取攔截爬升時間、攔截高度信息。從而獲得每一時刻的攻角信息,從而為攔截平臺裝訂發(fā)發(fā)射諸元完成預(yù)測命中點規(guī)劃。

    圖4 程序框圖Fig.4 Program chart

    3 仿真分析

    3.1 仿真輸入

    考慮攔截方利用探測系統(tǒng)提供的彈道數(shù)據(jù)預(yù)報了未來600 s的攻擊彈彈道,結(jié)果如圖5所示。

    表1 仿真輸入數(shù)據(jù)Table 1 Chemical composition of experimental alloy

    圖5 目標(biāo)預(yù)報軌跡Fig.5 Target prediction trajectory

    圖6 目標(biāo)距離攔截陣地航程Fig.6 Target and intercept position range

    由圖6可知目標(biāo)距離攔截陣地的斜距先變小,再變大,再變小,考慮到攔截器僅具有逆軌攔截能力,所以搜索區(qū)間為兩個航程—時間成反比的區(qū)間,即t1∈[0,440],t2∈[459,600]。

    3.2 標(biāo)準(zhǔn)彈道族設(shè)計仿真分析

    針對高超攔截平臺,這里基于沖壓發(fā)動機(jī)特性設(shè)計了包含巡航彈道與爬升彈道的標(biāo)準(zhǔn)彈道族。

    為驗證包含巡航彈道與爬升彈道的標(biāo)準(zhǔn)彈道族是否具有覆蓋所有攔截區(qū)域的能力,這里引入直線彈道。通過設(shè)置固定攔截高度,對比航程與爬升至指定高度的時間,驗證標(biāo)準(zhǔn)彈道族設(shè)計的正確性。

    假設(shè)攔截高度為29.5 km,攔截平臺飛行速度1800 m/s,直線爬升彈道的彈道傾角為0.4°。

    圖7 航程—時間曲線Fig.7 Range VS time

    圖8 高度—時間曲線Fig.8 Height VS time

    通過對比圖7和圖8發(fā)現(xiàn),巡航彈道在攔截初期能夠較早的爬升至預(yù)設(shè)攔截高度,而爬升彈道具有航程上的優(yōu)勢。直線彈道在航程與爬升至攔截高度的時間上均無優(yōu)勢,這也從側(cè)面說明了采用兩種彈道相組合的標(biāo)準(zhǔn)彈道族能夠有效的擴(kuò)展攔截空域。

    圖9 覆蓋空域示意圖Fig.9 Dchematic diagram of covered airspace

    3.3 預(yù)測命中點規(guī)劃仿真分析

    攔截平臺為巡航式飛行器,由超燃沖壓發(fā)動機(jī)提供動力,飛行攻角范圍為α∈[-4°,6°],為確保發(fā)動機(jī)正常工作,需要動壓滿足約束q∈[30 kPa, 70 kPa]。結(jié)合攔截平臺的飛行速度約束V∈[1800 m/s, 2100 m/s],攔截高度被限制為H∈[26 km, 30 km]。

    基于以上限制,假設(shè)前期獲得4個攔截射表,飛行速度分別為1800 m/s,1900 m/s,2000 m/s,2100 m/s,最大攔截允許誤差為1 km,并設(shè)置搜索區(qū)間為t1∈[0 s,440 s],t2∈[459 s,600 s]。

    在t1∈[0 s,440 s]區(qū)間,沒有搜索到可行解。在t2∈[459 s,600 s]區(qū)間,搜索到可行解。采用速度1800 m/s標(biāo)準(zhǔn)彈道族,攔截時間為503 s,發(fā)射方位角-19°,采用爬升彈道,切換時間為400 s,距離誤差146 m,搜索時間小于0.1 s,如圖10所示。

    圖10 軌跡示意圖Fig.10 Track diagram

    4 結(jié) 論

    本文瞄準(zhǔn)基于吸氣式高超聲速平臺的臨近空間目標(biāo)攔截概念方案。結(jié)合吸氣式動力飛行策略構(gòu)建高超聲速攔截平臺的標(biāo)準(zhǔn)彈道族;通過將多約束攔截預(yù)測命中點設(shè)計問題分解為基于航程約束的攔截時間迭代計算問題和滿足高度約束的攔截點的高效篩選問題,從而完成預(yù)測命中點規(guī)劃,為下一代反導(dǎo)攔截系統(tǒng)提供潛在的技術(shù)途徑。

    猜你喜歡
    航程超聲速彈道
    殲-16挑戰(zhàn)更大航程
    彈道——打勝仗的奧秘
    高超聲速出版工程
    高超聲速飛行器
    西進(jìn)執(zhí)教 一段人生的奇異航程
    海峽姐妹(2019年5期)2019-06-18 10:40:34
    一維彈道修正彈無線通信系統(tǒng)研制
    電子制作(2019年7期)2019-04-25 13:17:48
    超聲速旅行
    飛越北極的航程
    人生航程 “漫”條“思”理
    航海(2016年2期)2016-05-19 03:57:11
    基于PID控制的二維彈道修正彈仿真
    欧美一区二区精品小视频在线| 免费在线观看影片大全网站| 在线a可以看的网站| 亚洲无线观看免费| 亚洲av成人精品一区久久| 国产在线精品亚洲第一网站| 级片在线观看| 窝窝影院91人妻| 久久久久久久久中文| 亚洲va日本ⅴa欧美va伊人久久| av福利片在线观看| 一区二区三区激情视频| 身体一侧抽搐| 婷婷精品国产亚洲av在线| 久久精品国产自在天天线| 日韩欧美三级三区| 狂野欧美激情性xxxx在线观看| 观看美女的网站| 一区二区三区免费毛片| 欧美日韩乱码在线| 观看美女的网站| 久久久久九九精品影院| 国产美女午夜福利| 欧美中文日本在线观看视频| 精品一区二区免费观看| 免费av观看视频| 成人精品一区二区免费| 国产精品一区二区免费欧美| 亚洲人成网站高清观看| 韩国av在线不卡| 毛片女人毛片| 欧美3d第一页| 久久精品国产亚洲av涩爱 | 最新在线观看一区二区三区| 亚洲性久久影院| 午夜免费成人在线视频| 噜噜噜噜噜久久久久久91| 精品人妻熟女av久视频| 久久久久免费精品人妻一区二区| 18禁黄网站禁片免费观看直播| 少妇裸体淫交视频免费看高清| 麻豆精品久久久久久蜜桃| 午夜视频国产福利| 日韩欧美精品v在线| 免费大片18禁| 亚洲久久久久久中文字幕| 国产精品乱码一区二三区的特点| 成人国产麻豆网| 久久久久久久精品吃奶| av黄色大香蕉| 亚洲成人久久性| 亚洲熟妇熟女久久| 中文字幕精品亚洲无线码一区| 国产精品野战在线观看| 午夜影院日韩av| 极品教师在线免费播放| 国产精品久久电影中文字幕| 亚洲人成网站在线播放欧美日韩| 天堂√8在线中文| 久久久国产成人免费| 免费观看精品视频网站| 精品人妻1区二区| 欧美性感艳星| 欧美高清性xxxxhd video| 亚洲av五月六月丁香网| 亚洲自偷自拍三级| 春色校园在线视频观看| 观看免费一级毛片| 日韩中文字幕欧美一区二区| 有码 亚洲区| 可以在线观看的亚洲视频| 看片在线看免费视频| 午夜福利高清视频| 我要搜黄色片| 色尼玛亚洲综合影院| 麻豆成人av在线观看| 亚洲综合色惰| 国产真实乱freesex| 国产精品电影一区二区三区| 久久欧美精品欧美久久欧美| 久久久久久久精品吃奶| 亚洲在线观看片| 99久久久亚洲精品蜜臀av| 九色成人免费人妻av| 久99久视频精品免费| 天堂网av新在线| 国产精品综合久久久久久久免费| a级毛片免费高清观看在线播放| 在线免费观看不下载黄p国产 | 亚洲av熟女| 老司机福利观看| 免费无遮挡裸体视频| 久久99热这里只有精品18| 日韩高清综合在线| 黄色欧美视频在线观看| 免费看a级黄色片| 在线观看美女被高潮喷水网站| 久久草成人影院| 亚洲国产精品合色在线| 亚洲一级一片aⅴ在线观看| 在现免费观看毛片| 成人一区二区视频在线观看| 国产 一区精品| 欧美激情在线99| av黄色大香蕉| 波多野结衣高清无吗| 国内毛片毛片毛片毛片毛片| 国产 一区 欧美 日韩| 动漫黄色视频在线观看| 免费看美女性在线毛片视频| 亚洲天堂国产精品一区在线| 国产免费一级a男人的天堂| 色哟哟·www| 精品无人区乱码1区二区| av在线亚洲专区| 精品久久久久久久久亚洲 | 在线观看av片永久免费下载| 一夜夜www| 亚洲精品乱码久久久v下载方式| 三级男女做爰猛烈吃奶摸视频| 国产精品国产三级国产av玫瑰| 久久久国产成人免费| 日韩在线高清观看一区二区三区 | 热99re8久久精品国产| 韩国av一区二区三区四区| 日本色播在线视频| 男女视频在线观看网站免费| 啦啦啦韩国在线观看视频| 波多野结衣高清无吗| 最近中文字幕高清免费大全6 | 美女免费视频网站| 免费观看的影片在线观看| 免费观看精品视频网站| 欧美日韩精品成人综合77777| 男女视频在线观看网站免费| 欧美潮喷喷水| 国产av在哪里看| 亚洲欧美日韩卡通动漫| 啦啦啦观看免费观看视频高清| 久久精品人妻少妇| 日本免费a在线| 九九爱精品视频在线观看| 全区人妻精品视频| 美女xxoo啪啪120秒动态图| 人妻久久中文字幕网| 日本成人三级电影网站| 亚洲狠狠婷婷综合久久图片| 又爽又黄a免费视频| 国产精品一区二区性色av| 麻豆国产97在线/欧美| 欧美色欧美亚洲另类二区| 少妇的逼水好多| 亚洲电影在线观看av| 男女边吃奶边做爰视频| 好男人在线观看高清免费视频| 国产在线精品亚洲第一网站| www日本黄色视频网| 性色avwww在线观看| 看片在线看免费视频| 五月伊人婷婷丁香| 亚洲四区av| 中国美白少妇内射xxxbb| 少妇人妻精品综合一区二区 | 一级av片app| 成年版毛片免费区| 自拍偷自拍亚洲精品老妇| 97人妻精品一区二区三区麻豆| 村上凉子中文字幕在线| 欧美bdsm另类| 亚洲自偷自拍三级| av视频在线观看入口| 波多野结衣高清无吗| 亚洲真实伦在线观看| 国产精品,欧美在线| 国产精品一区www在线观看 | 桃红色精品国产亚洲av| 欧美人与善性xxx| 国内揄拍国产精品人妻在线| 国产成人av教育| 伊人久久精品亚洲午夜| 国产精品三级大全| 一级av片app| 熟妇人妻久久中文字幕3abv| 在线天堂最新版资源| 尤物成人国产欧美一区二区三区| 国产色婷婷99| 成年版毛片免费区| 亚洲美女黄片视频| 国产免费av片在线观看野外av| 亚洲av免费在线观看| 99久久无色码亚洲精品果冻| 最新在线观看一区二区三区| av在线亚洲专区| 日韩欧美国产一区二区入口| 琪琪午夜伦伦电影理论片6080| 别揉我奶头~嗯~啊~动态视频| 亚洲在线观看片| 日本黄大片高清| 天堂动漫精品| 久久国内精品自在自线图片| 精品久久久久久久久久久久久| 国产白丝娇喘喷水9色精品| 国产一区二区在线观看日韩| 91在线精品国自产拍蜜月| 欧美丝袜亚洲另类 | 亚洲成人精品中文字幕电影| 欧美成人a在线观看| 中文字幕久久专区| 久久久成人免费电影| 无人区码免费观看不卡| 国产精品久久久久久久电影| 成人午夜高清在线视频| 亚洲国产日韩欧美精品在线观看| 亚洲av成人av| 一本精品99久久精品77| 久久精品国产自在天天线| 简卡轻食公司| 国产精品久久久久久久电影| 夜夜爽天天搞| 两性午夜刺激爽爽歪歪视频在线观看| 男女边吃奶边做爰视频| 日本一本二区三区精品| 黄片wwwwww| 午夜爱爱视频在线播放| 日韩中文字幕欧美一区二区| 啪啪无遮挡十八禁网站| 成人永久免费在线观看视频| 国产黄a三级三级三级人| 特大巨黑吊av在线直播| 男人舔奶头视频| 一边摸一边抽搐一进一小说| 成人国产综合亚洲| 免费看美女性在线毛片视频| 国产在线男女| 天堂动漫精品| 国产综合懂色| 男人舔女人下体高潮全视频| 亚洲一区高清亚洲精品| 97碰自拍视频| 欧美色视频一区免费| 欧美日韩瑟瑟在线播放| 亚洲性夜色夜夜综合| 国产aⅴ精品一区二区三区波| 日韩欧美精品v在线| 久9热在线精品视频| 伦理电影大哥的女人| 色综合站精品国产| 亚洲三级黄色毛片| 国产精品免费一区二区三区在线| 一卡2卡三卡四卡精品乱码亚洲| 亚洲专区中文字幕在线| 日本欧美国产在线视频| 九九在线视频观看精品| 最近中文字幕高清免费大全6 | 丰满人妻一区二区三区视频av| 干丝袜人妻中文字幕| 精品一区二区三区人妻视频| a在线观看视频网站| 久久精品国产亚洲av涩爱 | 天堂动漫精品| 最近最新中文字幕大全电影3| 无遮挡黄片免费观看| 精品一区二区三区视频在线| 亚洲成人免费电影在线观看| 天天躁日日操中文字幕| 日本熟妇午夜| 搡老妇女老女人老熟妇| 人人妻人人看人人澡| 国产91精品成人一区二区三区| 成人鲁丝片一二三区免费| 国产精品综合久久久久久久免费| 在现免费观看毛片| 久久香蕉精品热| 国产黄色小视频在线观看| 成人欧美大片| 国产男人的电影天堂91| 国产精品一区二区免费欧美| 99久久中文字幕三级久久日本| 极品教师在线视频| 色在线成人网| 国产成人av教育| 永久网站在线| 毛片女人毛片| 久久精品91蜜桃| 久久久久久久精品吃奶| 国产成人一区二区在线| 我要看日韩黄色一级片| 赤兔流量卡办理| 国产精品人妻久久久影院| 亚洲avbb在线观看| 色噜噜av男人的天堂激情| 久久久久性生活片| 亚洲精品亚洲一区二区| 能在线免费观看的黄片| 91久久精品电影网| 好男人在线观看高清免费视频| 在线a可以看的网站| 日日干狠狠操夜夜爽| 亚洲国产日韩欧美精品在线观看| or卡值多少钱| 中文字幕av在线有码专区| 亚洲自偷自拍三级| 亚洲在线自拍视频| 天堂动漫精品| 国产色婷婷99| 亚洲av成人av| 国产探花在线观看一区二区| 国产亚洲91精品色在线| 国产成人a区在线观看| 他把我摸到了高潮在线观看| or卡值多少钱| 老司机午夜福利在线观看视频| 国产三级中文精品| 亚洲四区av| av天堂中文字幕网| 欧美成人一区二区免费高清观看| 日韩国内少妇激情av| 综合色av麻豆| 色综合婷婷激情| 女生性感内裤真人,穿戴方法视频| 精品久久久久久久久久免费视频| 久久中文看片网| 老司机福利观看| 亚洲欧美日韩无卡精品| 干丝袜人妻中文字幕| 成人毛片a级毛片在线播放| a在线观看视频网站| 99久久九九国产精品国产免费| 日本精品一区二区三区蜜桃| 少妇裸体淫交视频免费看高清| 最新中文字幕久久久久| а√天堂www在线а√下载| 一级av片app| 国产精品久久久久久久电影| 久久久久久大精品| 别揉我奶头 嗯啊视频| av天堂在线播放| 精品久久久久久久末码| 国产av麻豆久久久久久久| 美女大奶头视频| 在线观看午夜福利视频| 成人鲁丝片一二三区免费| 99久久精品一区二区三区| 国产极品精品免费视频能看的| 亚洲18禁久久av| 最近在线观看免费完整版| 在线播放国产精品三级| 亚洲成人中文字幕在线播放| 成人亚洲精品av一区二区| 精华霜和精华液先用哪个| 欧美3d第一页| 黄片wwwwww| 精品午夜福利视频在线观看一区| 午夜a级毛片| 亚洲av不卡在线观看| 一级a爱片免费观看的视频| 色视频www国产| 中文字幕av成人在线电影| 麻豆成人午夜福利视频| 22中文网久久字幕| 欧美色视频一区免费| 不卡一级毛片| 熟女人妻精品中文字幕| 欧美又色又爽又黄视频| 麻豆av噜噜一区二区三区| 日日啪夜夜撸| 欧美3d第一页| 色噜噜av男人的天堂激情| 老司机福利观看| 3wmmmm亚洲av在线观看| 内地一区二区视频在线| 免费av毛片视频| 高清在线国产一区| 免费不卡的大黄色大毛片视频在线观看 | 日日干狠狠操夜夜爽| 日本熟妇午夜| 亚洲精华国产精华精| 老熟妇乱子伦视频在线观看| 97热精品久久久久久| .国产精品久久| 国模一区二区三区四区视频| 午夜福利在线观看免费完整高清在 | 久久久久久伊人网av| 直男gayav资源| 亚洲综合色惰| 亚洲欧美日韩高清在线视频| 精品久久久久久久久av| 免费一级毛片在线播放高清视频| 91久久精品国产一区二区三区| 中文字幕av在线有码专区| 嫩草影院精品99| 国产亚洲精品av在线| 欧美日韩乱码在线| 网址你懂的国产日韩在线| 亚洲精品影视一区二区三区av| 国产高清不卡午夜福利| 一本一本综合久久| 日韩精品青青久久久久久| 非洲黑人性xxxx精品又粗又长| 午夜免费成人在线视频| 极品教师在线视频| 他把我摸到了高潮在线观看| 国内精品久久久久精免费| 亚洲国产欧美人成| 国产精品99久久久久久久久| 女人十人毛片免费观看3o分钟| 黄色一级大片看看| 亚洲va在线va天堂va国产| 精品一区二区三区视频在线观看免费| 亚洲成人久久爱视频| 久久久午夜欧美精品| 夜夜夜夜夜久久久久| 日本色播在线视频| 99久久成人亚洲精品观看| 免费不卡的大黄色大毛片视频在线观看 | 亚洲精华国产精华液的使用体验 | 欧美日韩中文字幕国产精品一区二区三区| 欧美日本视频| 色播亚洲综合网| 波野结衣二区三区在线| 春色校园在线视频观看| 亚洲 国产 在线| 日本黄色片子视频| 亚洲美女视频黄频| 黄色欧美视频在线观看| 嫩草影院精品99| 天天一区二区日本电影三级| 午夜福利在线观看吧| 永久网站在线| 久久精品国产亚洲网站| 国产不卡一卡二| 国产在视频线在精品| 波多野结衣高清作品| 国产在线男女| 日本爱情动作片www.在线观看 | 免费av毛片视频| 天堂网av新在线| 九九热线精品视视频播放| 欧美激情国产日韩精品一区| 欧美黑人巨大hd| 熟女电影av网| 久99久视频精品免费| 自拍偷自拍亚洲精品老妇| 性欧美人与动物交配| 欧美又色又爽又黄视频| 好男人在线观看高清免费视频| 亚洲成av人片在线播放无| 亚洲成人中文字幕在线播放| 精品欧美国产一区二区三| 精品不卡国产一区二区三区| 偷拍熟女少妇极品色| 日日夜夜操网爽| 久久人人爽人人爽人人片va| 桃色一区二区三区在线观看| 日本-黄色视频高清免费观看| 欧美zozozo另类| 少妇高潮的动态图| 久久亚洲真实| 小说图片视频综合网站| 一边摸一边抽搐一进一小说| 少妇高潮的动态图| 男女视频在线观看网站免费| 三级男女做爰猛烈吃奶摸视频| 亚洲五月天丁香| 免费av不卡在线播放| 国内毛片毛片毛片毛片毛片| 久久久久久久亚洲中文字幕| 亚洲综合色惰| 国产综合懂色| 精品一区二区免费观看| 中文字幕免费在线视频6| 国产精品亚洲一级av第二区| 女人十人毛片免费观看3o分钟| 赤兔流量卡办理| 长腿黑丝高跟| 真实男女啪啪啪动态图| www.色视频.com| 国产毛片a区久久久久| 变态另类丝袜制服| 真实男女啪啪啪动态图| 色噜噜av男人的天堂激情| 亚洲一区高清亚洲精品| 亚洲专区中文字幕在线| 日本一二三区视频观看| 99视频精品全部免费 在线| 毛片一级片免费看久久久久 | 老熟妇仑乱视频hdxx| 国产亚洲av嫩草精品影院| 欧美区成人在线视频| 国模一区二区三区四区视频| 亚洲av五月六月丁香网| 亚洲av第一区精品v没综合| 韩国av一区二区三区四区| 国产精品日韩av在线免费观看| 黄色欧美视频在线观看| 亚洲精品色激情综合| 淫妇啪啪啪对白视频| 九九久久精品国产亚洲av麻豆| 别揉我奶头 嗯啊视频| 香蕉av资源在线| 99热这里只有是精品50| 少妇被粗大猛烈的视频| 又爽又黄无遮挡网站| 久久国产乱子免费精品| 老司机福利观看| 日韩中字成人| 午夜激情欧美在线| 日韩精品青青久久久久久| 欧美成人a在线观看| 国产老妇女一区| 久久国内精品自在自线图片| 免费观看在线日韩| 99在线视频只有这里精品首页| 国产精品一区二区性色av| 精品欧美国产一区二区三| 变态另类成人亚洲欧美熟女| 丝袜美腿在线中文| 国内精品一区二区在线观看| 一区二区三区免费毛片| 久久欧美精品欧美久久欧美| 国内久久婷婷六月综合欲色啪| 国产一区二区亚洲精品在线观看| 伦理电影大哥的女人| 小蜜桃在线观看免费完整版高清| 岛国在线免费视频观看| 日韩欧美 国产精品| 男人狂女人下面高潮的视频| 亚州av有码| 搞女人的毛片| 给我免费播放毛片高清在线观看| 欧美日本亚洲视频在线播放| 日本一本二区三区精品| 一夜夜www| 三级国产精品欧美在线观看| 女的被弄到高潮叫床怎么办 | 亚洲无线观看免费| 毛片女人毛片| 美女免费视频网站| 日日夜夜操网爽| 亚洲国产精品成人综合色| 国产高潮美女av| 精品99又大又爽又粗少妇毛片 | 亚洲国产精品合色在线| 简卡轻食公司| 久久精品国产亚洲av香蕉五月| 国产伦精品一区二区三区四那| 一个人观看的视频www高清免费观看| 22中文网久久字幕| 直男gayav资源| 国产一区二区三区在线臀色熟女| av福利片在线观看| 99热这里只有是精品在线观看| 美女黄网站色视频| 日韩国内少妇激情av| 久久久国产成人免费| 一区福利在线观看| 看片在线看免费视频| 久久午夜福利片| 神马国产精品三级电影在线观看| 天堂影院成人在线观看| 在线观看美女被高潮喷水网站| 国产一区二区三区视频了| 久久精品夜夜夜夜夜久久蜜豆| 黄色丝袜av网址大全| 听说在线观看完整版免费高清| 成人特级黄色片久久久久久久| 亚洲一区高清亚洲精品| 久久精品影院6| h日本视频在线播放| 国内精品久久久久久久电影| 97超视频在线观看视频| 少妇熟女aⅴ在线视频| 午夜精品久久久久久毛片777| 国产一区二区在线观看日韩| 波多野结衣巨乳人妻| 亚洲精品久久国产高清桃花| 国产一区二区三区av在线 | 国内久久婷婷六月综合欲色啪| 国产精品自产拍在线观看55亚洲| 动漫黄色视频在线观看| 成熟少妇高潮喷水视频| a级毛片免费高清观看在线播放| 国产老妇女一区| 久久精品国产99精品国产亚洲性色| 国产欧美日韩精品亚洲av| 国产乱人伦免费视频| 欧美色欧美亚洲另类二区| 国产精品国产三级国产av玫瑰| 天天一区二区日本电影三级| 国产美女午夜福利| 久久中文看片网| 看免费成人av毛片| 国产精品美女特级片免费视频播放器| 亚洲综合色惰| 天天一区二区日本电影三级| 狂野欧美白嫩少妇大欣赏| 香蕉av资源在线| 天天一区二区日本电影三级| 一区二区三区免费毛片| 悠悠久久av| 一夜夜www| x7x7x7水蜜桃| 少妇人妻一区二区三区视频| 中出人妻视频一区二区| 男女做爰动态图高潮gif福利片| 国产aⅴ精品一区二区三区波| 天堂√8在线中文| 又黄又爽又免费观看的视频| 久久精品国产亚洲网站| 亚洲熟妇中文字幕五十中出| 婷婷精品国产亚洲av| 成人鲁丝片一二三区免费| 亚洲精品粉嫩美女一区| 变态另类成人亚洲欧美熟女| 91精品国产九色|