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

    主動(dòng)防御非奇異終端滑模協(xié)同制導(dǎo)律

    2015-02-28 10:45:00鄒昕光周荻杜潤(rùn)樂劉佳琪
    兵工學(xué)報(bào) 2015年3期
    關(guān)鍵詞:攔截導(dǎo)彈視線制導(dǎo)

    鄒昕光,周荻,杜潤(rùn)樂,劉佳琪

    (1.哈爾濱工業(yè)大學(xué) 電氣工程及自動(dòng)化學(xué)院,黑龍江 哈爾濱150001;2.哈爾濱工業(yè)大學(xué) 航天學(xué)院,黑龍江 哈爾濱150001;3.試驗(yàn)物理與計(jì)算數(shù)學(xué)國(guó)家級(jí)重點(diǎn)實(shí)驗(yàn)室,北京100076)

    0 引言

    先進(jìn)導(dǎo)彈防御系統(tǒng)的研究和開發(fā)給彈道導(dǎo)彈突防帶來了嚴(yán)峻的考驗(yàn)。彈道導(dǎo)彈飛行中段處在大氣層外,此時(shí)推力發(fā)動(dòng)機(jī)關(guān)閉,彈頭和彈體分離,導(dǎo)彈所具有的機(jī)動(dòng)能力很小或者可以忽略不計(jì),導(dǎo)致此階段彈道導(dǎo)彈的飛行軌跡比較固定,易被導(dǎo)彈防御系統(tǒng)攔截。因此研究有效的防御手段以提高彈道導(dǎo)彈突防成功率具有較大的工程實(shí)際意義和理論價(jià)值。

    目前已經(jīng)提出了一些提高彈道導(dǎo)彈突防成功概率的方法,例如電子干擾,偽裝技術(shù),誘餌技術(shù),使用中制導(dǎo)和末制導(dǎo),多彈頭技術(shù),機(jī)動(dòng)變軌技術(shù)等[1]。這些方法本質(zhì)上都可以歸結(jié)為被動(dòng)防御,其出發(fā)點(diǎn)都是尋找有效的躲避手段。與此相對(duì),彈道導(dǎo)彈也可以采取主動(dòng)的防御手段:在受到威脅時(shí)釋放一個(gè)或多個(gè)防御導(dǎo)彈來保護(hù)自身。防御導(dǎo)彈具有和彈道導(dǎo)彈相同的速度,在距離彈道導(dǎo)彈一定距離的范圍內(nèi)伴飛,采用有效的制導(dǎo)律攔截來襲的攔截導(dǎo)彈[1]。為了敘述簡(jiǎn)潔清晰,本文后面將突防的彈道導(dǎo)彈稱為目標(biāo),將導(dǎo)彈防御系統(tǒng)發(fā)射的攔截導(dǎo)彈稱為攔截導(dǎo)彈,將與彈道導(dǎo)彈伴飛的反攔截導(dǎo)彈稱為防御導(dǎo)彈。本文考慮只存在一枚防御導(dǎo)彈的場(chǎng)景。

    由于主動(dòng)防御方式的特點(diǎn)和有效性,該方法在飛機(jī)防御中得到了廣泛的研究。文獻(xiàn)[2]在假設(shè)攔截導(dǎo)彈制導(dǎo)律是已知線性制導(dǎo)律的前提下推導(dǎo)了防御導(dǎo)彈和目標(biāo)的協(xié)同制導(dǎo)律。但假設(shè)攔截導(dǎo)彈的制導(dǎo)律已知并不符合實(shí)際。針對(duì)該問題文獻(xiàn)[3]提出了一種多模自適應(yīng)制導(dǎo)律,該制導(dǎo)律通過一組濾波器來確定攔截導(dǎo)彈的制導(dǎo)律和估計(jì)制導(dǎo)律參數(shù)。但是該方法要求攔截導(dǎo)彈的制導(dǎo)律和制導(dǎo)參數(shù)的組合在一個(gè)有限的已知集合中選取,而集合的大小受限于濾波器的個(gè)數(shù)。這限制了該制導(dǎo)律在工程實(shí)際中的運(yùn)用。文獻(xiàn)[4]根據(jù)防御導(dǎo)彈和目標(biāo)的不同協(xié)同方式(單向協(xié)同或者雙向協(xié)同)設(shè)計(jì)出了3 種不同的最優(yōu)主動(dòng)防御制導(dǎo)律。使用狀態(tài)方程降階方法,將高階運(yùn)動(dòng)方程降階為以零控脫靶量為狀態(tài)的一階動(dòng)態(tài)方程,簡(jiǎn)化了制導(dǎo)律的推導(dǎo)。以上工作都基于最優(yōu)控制理論,最優(yōu)制導(dǎo)律的設(shè)計(jì)需要估計(jì)剩余時(shí)間參數(shù)。而剩余時(shí)間參數(shù)的估計(jì)精度將直接影響制導(dǎo)的精度。攔截導(dǎo)彈一方面需要攔截目標(biāo),另一方面需要逃避防御導(dǎo)彈的攔截,這是一個(gè)微分對(duì)策問題,因此一些文獻(xiàn)使用了微分對(duì)策理論來研究主動(dòng)防御問題[5-6]。和最優(yōu)制導(dǎo)類似,微分對(duì)策制導(dǎo)律需要估計(jì)終端時(shí)刻,而終端時(shí)刻的估計(jì)精度會(huì)直接影響到制導(dǎo)精度。

    如果將目標(biāo),攔截導(dǎo)彈和防御導(dǎo)彈抽象為空間中的點(diǎn),那么它們之間的連線構(gòu)成一個(gè)三角形。一種可行的制導(dǎo)策略是防御導(dǎo)彈控制自身處于目標(biāo)和攔截導(dǎo)彈的視線上且在目標(biāo)和攔截導(dǎo)彈之間,此時(shí)三角形退化為直線,隨著攔截導(dǎo)彈不斷接近目標(biāo)最終將撞上防御導(dǎo)彈[7-9]。這種制導(dǎo)策略稱為視線三角制導(dǎo)策略,在幾何上比較直觀。文獻(xiàn)[7]提出的制導(dǎo)律采用了視線三角制導(dǎo)策略,制導(dǎo)律的表達(dá)式中包含了攔截導(dǎo)彈的加速度項(xiàng),但在實(shí)際中對(duì)于防御導(dǎo)彈而言,攔截導(dǎo)彈的加速度是難以獲得的量。在文獻(xiàn)[8 -9]提出的基于視線三角制導(dǎo)策略的制導(dǎo)律中,防御導(dǎo)彈的制導(dǎo)指令和兩個(gè)參數(shù)成比例。第1 個(gè)參數(shù)是防御導(dǎo)彈-攔截導(dǎo)彈視線角和攔截導(dǎo)彈-目標(biāo)視線角的差。第2 個(gè)參數(shù)是防御導(dǎo)彈和攔截導(dǎo)彈的相對(duì)速度。

    滑??刂朴捎趯?duì)建模誤差和外部未知干擾的良好魯棒性被廣泛應(yīng)用在導(dǎo)彈制導(dǎo)律的設(shè)計(jì)中[10-22]。文獻(xiàn)[13]提出了基于滑模的主動(dòng)防御協(xié)作制導(dǎo)律,研究了大氣層內(nèi)飛機(jī)防御導(dǎo)彈攻擊的問題。有限時(shí)間穩(wěn)定性理論在1986年被首次提出[14],由于具備有保障的收斂性時(shí)間特性,該理論在制導(dǎo)律設(shè)計(jì)問題中被廣泛使用[15-26]。近年來有限時(shí)間收斂滑模控制理論主要用來解決含有角度約束的制導(dǎo)問題,但是一些制導(dǎo)律在被控角度誤差很小時(shí)會(huì)出現(xiàn)奇異問題,導(dǎo)致大的抖振和控制飽和[17-18]。非奇異終端滑模控制方法可以控制滑模變量在有限時(shí)間內(nèi)到達(dá)滑模面并且避免由于奇異問題而造成的控制飽和[19-22]。文獻(xiàn)[23]和文獻(xiàn)[24]針對(duì)該問題分別設(shè)計(jì)了不同的非奇異終端滑模制導(dǎo)律,該制導(dǎo)律在避免奇異問題的同時(shí)保持有限時(shí)間收斂的特性。

    本文試圖設(shè)計(jì)一種大氣層外防御導(dǎo)彈高性能的魯棒制導(dǎo)律。考慮到在制導(dǎo)的大部分時(shí)間內(nèi)目標(biāo)和防御導(dǎo)彈之間的距離和他們與攔截導(dǎo)彈之間的距離相比是一個(gè)小量,將制導(dǎo)分為兩個(gè)階段。在第1 階段對(duì)相對(duì)運(yùn)動(dòng)近似建模,使用基于視線三角制導(dǎo)策略的制導(dǎo)律;第2 階段考慮到近似建模誤差隨著時(shí)間快速增加,因此使用基于零化視線轉(zhuǎn)率策略的制導(dǎo)律。第2 階段的制導(dǎo)律使用更加精確的相對(duì)運(yùn)動(dòng)模型。制導(dǎo)律的設(shè)計(jì)使用非奇異終端滑??刂评碚?。本文在理論上證明基于視線三角制導(dǎo)策略可以顯著地降低防御導(dǎo)彈的需用過載,而基于零化視線轉(zhuǎn)率的制導(dǎo)律使用更加精確的相對(duì)運(yùn)動(dòng)模型,因此分階段的制導(dǎo)方法兼顧了二者的優(yōu)點(diǎn)。仿真驗(yàn)證了制導(dǎo)律的有效性。

    1 問題描述

    目標(biāo)、攔截導(dǎo)彈和防御導(dǎo)彈的相對(duì)運(yùn)動(dòng)關(guān)系如圖1所示,其中目標(biāo)、攔截導(dǎo)彈和防御導(dǎo)彈分別用T、M 和D 來表示。目標(biāo)、攔截導(dǎo)彈和防御導(dǎo)彈的速度分別用vT、vM和vD表示。攔截導(dǎo)彈-目標(biāo)的視線角為qεMT和qβMT;防御導(dǎo)彈-攔截導(dǎo)彈的視線角為qεDM和qβDM. OXYZ 為慣性參考系,其中目標(biāo)和防御導(dǎo)彈的慣性參考系TXTYTZT和DXDYDZD方向和OXYZ 相同,攔截導(dǎo)彈的慣性參考系MXMYMZM方向?yàn)樽鴺?biāo)系OXYZ 繞OY 軸旋轉(zhuǎn)180°.

    圖1 目標(biāo)、攔截導(dǎo)彈和防御導(dǎo)彈相對(duì)運(yùn)動(dòng)幾何關(guān)系Fig.1 Interception geometry in three dimensions

    攔截導(dǎo)彈-目標(biāo)和防御導(dǎo)彈-攔截導(dǎo)彈解耦的相對(duì)運(yùn)動(dòng)方程為

    式中:aTqε和aTqβ代表目標(biāo)加速度在防御導(dǎo)彈-攔截導(dǎo)彈視線系下的縱向平面和側(cè)向平面分量;aDqε和aDqβ代表防御導(dǎo)彈加速度在防御導(dǎo)彈-攔截導(dǎo)彈視線系下的縱向平面和側(cè)向平面分量;aM1qε和aM1qβ代表攔截導(dǎo)彈加速度在攔截導(dǎo)彈-目標(biāo)視線系下的縱向平面和側(cè)向平面分量;aM2qε和aM2qβ代表攔截導(dǎo)彈加速度在防御導(dǎo)彈-攔截導(dǎo)彈視線系下的縱向平面和側(cè)向平面分量。

    2 非奇異終端滑??刂评碚?/h2>

    在推導(dǎo)制導(dǎo)律表達(dá)式之前先給出制導(dǎo)律采用的有限時(shí)間收斂控制理論。

    引理1 對(duì)于非線性系統(tǒng)

    假設(shè)存在一個(gè)定義在原點(diǎn)鄰域U??n上的C1光滑函數(shù)V(x,t),并且存在實(shí)數(shù)γ >0 和0 <λ <1,使得V(x,t)在U 上正定且(x,t)+γVλ(x,t)在U上半負(fù)定,則系統(tǒng)是有限時(shí)間穩(wěn)定的[25-26]。

    證明見文獻(xiàn)[25 -26],為了簡(jiǎn)化計(jì)算,取V(x,t)=,σ 是滑模變量,取根據(jù)引理1,系統(tǒng)有限時(shí)間穩(wěn)定的充分條件有

    (6)式為有限時(shí)間滑模到達(dá)條件。到達(dá)時(shí)間為

    引理2 對(duì)于如下2 階非線性系統(tǒng)

    式中:x=[x1,x2]T∈?2. 將其滑模變量設(shè)置為σ=x1+,β >0,1 <α =p/q <2,p、q 都是正奇數(shù)。引理1 的有限時(shí)間收斂條件,可以保證狀態(tài)x1、在有限時(shí)間內(nèi)收斂到0[19-23].

    狀態(tài)在時(shí)間ts之內(nèi)收斂到0,ts表達(dá)式為

    根據(jù)引理1,tr是有限正值。由于p、q 都是奇數(shù)且p >q >0,β >0,有所以ts≥tr.

    3 制導(dǎo)律設(shè)計(jì)

    在防御導(dǎo)彈制導(dǎo)的大多數(shù)時(shí)間內(nèi)防御導(dǎo)彈和目標(biāo)之間距離相對(duì)于它們與攔截導(dǎo)彈的距離小得多,而且防御導(dǎo)彈和目標(biāo)速度矢量十分接近,因此在大部分時(shí)間里滿足如下關(guān)系:

    將整個(gè)制導(dǎo)分為兩個(gè)階段:階段1 和階段2. 為了簡(jiǎn)化制導(dǎo)律的表達(dá)式,在階段1 中認(rèn)為(9)式和(10)式成立,防御導(dǎo)彈采用視線三角制導(dǎo)策略,保持自身在目標(biāo)和攔截導(dǎo)彈的視線上,這種制導(dǎo)策略可以有效降低防御導(dǎo)彈的需用過載。在階段2 中認(rèn)為(9)式和(10)式不再成立,防御導(dǎo)彈采用零化對(duì)攔截導(dǎo)彈視線轉(zhuǎn)率的制導(dǎo)策略。假設(shè)防御導(dǎo)彈可以測(cè)量相對(duì)于攔截導(dǎo)彈的距離RDM和速度,目標(biāo)可以測(cè)量相對(duì)于攔截導(dǎo)彈的相對(duì)距離RMT和速度并通過通信設(shè)備將RMT和發(fā)送給防御導(dǎo)彈。防御導(dǎo)彈使用一個(gè)門限值μ*來區(qū)分這兩個(gè)階段。并設(shè)

    當(dāng)μ >μ*時(shí),防御導(dǎo)彈使用階段1 制導(dǎo)律,當(dāng)μ≤μ*防御導(dǎo)彈使用階段2 制導(dǎo)律。由于防御導(dǎo)彈總是在目標(biāo)的前下方伴飛,所以0 <RDM<RMT. 當(dāng)μ*=0 時(shí)防御導(dǎo)彈全程使用階段1 制導(dǎo)律;當(dāng)μ*=1 時(shí)防御導(dǎo)彈全程使用階段2 制導(dǎo)律。

    3.1 階段1 制導(dǎo)律設(shè)計(jì)

    為了保證防御導(dǎo)彈能擊中攔截導(dǎo)彈,防御導(dǎo)彈的制導(dǎo)策略是將自己的位置保持在目標(biāo)和導(dǎo)彈連線上,同時(shí)零化視線轉(zhuǎn)率。即

    在該階段中,近似地認(rèn)為(9)式和(10)式成立。目標(biāo)在大氣層外具有很小的機(jī)動(dòng)能力,本文忽略目標(biāo)的機(jī)動(dòng)。因此在該階段中攔截導(dǎo)彈和目標(biāo)的相對(duì)運(yùn)動(dòng)方程(1)式和(2)式可以近似為

    利用非奇異終端滑模設(shè)計(jì)縱向平面制導(dǎo)律,選取滑動(dòng)模態(tài)面為

    式中:β1為常量,β1>0;1 <α1=p1/q1<2,p1、q1都是正奇數(shù)。設(shè)計(jì)目標(biāo)是控制aDqε,保證狀態(tài)xε1和xε2有限時(shí)間內(nèi)趨向0. 對(duì)(16)式取時(shí)間導(dǎo)數(shù)得

    假設(shè)防御導(dǎo)彈-攔截導(dǎo)彈和攔截導(dǎo)彈-目標(biāo)的視線角速率可以分別被防御導(dǎo)彈和目標(biāo)測(cè)量,目標(biāo)通過通信設(shè)備將測(cè)量值發(fā)送給防御導(dǎo)彈,則xε2可以計(jì)算得到。攔截導(dǎo)彈的加速度不能獲得但是可以假定攔截導(dǎo)彈的機(jī)動(dòng)能力是有限的。(17)式中攔截導(dǎo)彈加速度aMqε=aM2qε-aM1qε是攔截導(dǎo)彈加速度在防御導(dǎo)彈-攔截導(dǎo)彈視線系和攔截導(dǎo)彈-目標(biāo)視線系下Y 軸分量的差。按照本制導(dǎo)律設(shè)計(jì)思路,這兩個(gè)視線系的Y 軸接近平行且方向相同,因此aM2qε≈aM1qε,因此攔截導(dǎo)彈在(17)式中引入的加速度很小,這是該制導(dǎo)策略的一個(gè)優(yōu)勢(shì)。假設(shè)|aMqε|≤M1,將攔截導(dǎo)彈的機(jī)動(dòng)看成是系統(tǒng)擾動(dòng)??v向平面制導(dǎo)律設(shè)計(jì)為

    式中:ρ1=M1+1.

    類似的,利用非奇異終端滑模設(shè)計(jì)側(cè)向平面制導(dǎo)律,選取滑動(dòng)模態(tài)面為

    式中:β2為常量都是正奇數(shù)。設(shè)計(jì)目標(biāo)是控制aDqβ,保證狀態(tài)xβ1和xβ2有限時(shí)間內(nèi)趨向0. 對(duì)(16)式取時(shí)間導(dǎo)數(shù)得

    類似地假設(shè)防御導(dǎo)彈-攔截導(dǎo)彈和攔截導(dǎo)彈-目標(biāo)的視線角速率都可以測(cè)量,則xβ2可以計(jì)算得到;假設(shè)防御導(dǎo)彈-攔截導(dǎo)彈的相對(duì)距離和相對(duì)速度可以測(cè)量。攔截導(dǎo)彈的加速度不能獲得但是可以假定攔截導(dǎo)彈的機(jī)動(dòng)能力是有限的。(20)式中攔截導(dǎo)彈加速度aMqβ=aM1qβ+aM2qβ是攔截導(dǎo)彈加速度在防御導(dǎo)彈-攔截導(dǎo)彈視線系和攔截導(dǎo)彈-目標(biāo)視線系下Z 軸分量的和。按照本制導(dǎo)律設(shè)計(jì)思路,這兩個(gè)視線系的Z 軸接近平行且方向相反,因此aM2qβ≈-aM1qβ,因此攔截導(dǎo)彈在(20)式中引入的加速度很小。假設(shè)|aMqβ|≤M2,將攔截導(dǎo)彈的機(jī)動(dòng)看成是系統(tǒng)擾動(dòng)。側(cè)向平面制導(dǎo)律設(shè)計(jì)為

    式中:ρ2=M2+1.

    定理1 制導(dǎo)系統(tǒng)(14)式和(15)式的制導(dǎo)律分別取(18)式和(21)式時(shí),可以保證狀態(tài)xε1、、xβ1和在有限時(shí)間內(nèi)趨向0,即防御導(dǎo)彈能夠在有限時(shí)間內(nèi)將自身位置調(diào)整到目標(biāo)和攔截導(dǎo)彈的視線上并且保證防御導(dǎo)彈-攔截導(dǎo)彈的視線轉(zhuǎn)率漸進(jìn)趨向0.

    證明 根據(jù)(17)式與(18)式有

    根據(jù)(20)式和(21)式有

    由于σ1和σ2采用了非奇異終端滑模形式,根據(jù)引理2 知道xε1、xβ1和在有限時(shí)間內(nèi)收斂到0. 收斂時(shí)間分別為

    考慮到一般情況下攔截導(dǎo)彈都使用零化視線轉(zhuǎn)率策略的制導(dǎo)律對(duì)目標(biāo)制導(dǎo),且因?yàn)楹驮谟邢迺r(shí)間內(nèi)收斂到0,所以防御導(dǎo)彈對(duì)攔截導(dǎo)彈的視線轉(zhuǎn)率會(huì)隨著收斂到0.

    證畢。

    3.2 階段2 制導(dǎo)律設(shè)計(jì)

    在階段2 中,防御導(dǎo)彈采用零化視線轉(zhuǎn)率的制導(dǎo)策略。此時(shí)防御導(dǎo)彈的制導(dǎo)律不再需要考慮目標(biāo)的運(yùn)動(dòng)信息。

    利用非奇異終端滑模設(shè)計(jì)縱向平面制導(dǎo)律,選取滑動(dòng)模態(tài)面為

    設(shè)計(jì)目標(biāo)是控制aDqε,保證 狀 態(tài)有限 時(shí) 間內(nèi)趨向0. 對(duì)(26)式取時(shí)間的導(dǎo)數(shù)得

    攔截導(dǎo)彈加速度難以測(cè)量但是可以假定其是一個(gè)有限值,且門限已知,即

    和階段1 的制導(dǎo)策略相比,(27)式中的攔截導(dǎo)彈加速度項(xiàng)不再具有相互抵消的特點(diǎn),因而這種制導(dǎo)策略下引入的攔截導(dǎo)彈的加速度會(huì)比較大。制導(dǎo)律設(shè)計(jì)為

    式中:ρ3=M3+1.

    類似的,利用非奇異終端滑模設(shè)計(jì)側(cè)向平面制導(dǎo)律,選取滑動(dòng)模態(tài)面為

    設(shè)計(jì)目標(biāo)是控制aDqβ,保證狀態(tài)q·βDM有限時(shí)間內(nèi)趨向0. 對(duì)(30)式取時(shí)間的導(dǎo)數(shù)得

    攔截導(dǎo)彈加速度難以測(cè)量但是可以假定其是一個(gè)有限值,且門限已知,即

    制導(dǎo)律設(shè)計(jì)為

    式中:ρ4=M4+1.

    定理2 制導(dǎo)系統(tǒng)(3)式和(4)式的制導(dǎo)律分別取(29)式和(33)式時(shí),可以保證防御導(dǎo)彈-攔截導(dǎo)彈的視線轉(zhuǎn)率在有限時(shí)間收斂到0.

    證明 根據(jù)(27)式和(29)式有

    設(shè)γ3則

    根據(jù)(31)式和(33)式有

    設(shè)γ4則

    根據(jù)引理1 可知滑模變量σ3和σ4在有限時(shí)間內(nèi)達(dá)到滑模面。收斂時(shí)間分別為

    定理得證。

    4 仿真結(jié)果及分析

    在大氣層外防御導(dǎo)彈通過軌控發(fā)動(dòng)機(jī)產(chǎn)生加速度。而軌控發(fā)動(dòng)機(jī)一般是開關(guān)式發(fā)動(dòng)機(jī)輸出的推力只能是固定值或者0,這不適合制導(dǎo)指令是變化的制導(dǎo)律。針對(duì)該問題,使用了基于脈沖寬度調(diào)制(PWM)的軌控發(fā)動(dòng)機(jī)控制方法,基于沖量等效原理,輸出可變的等效力,來適應(yīng)本文設(shè)計(jì)的制導(dǎo)律。

    下面將本文提出的分階段主動(dòng)防御非奇異終端滑模制導(dǎo)律簡(jiǎn)稱為AD_NTSM 制導(dǎo)律。為了考察AD_NTSM 制導(dǎo)律的性能,先給出一個(gè)仿真算例;接著考察在不同μ*門限值下的制導(dǎo)律性能,以期對(duì)于不同的防御導(dǎo)彈機(jī)動(dòng)能力找到對(duì)應(yīng)最佳的μ*值,兼顧兩個(gè)階段制導(dǎo)律的優(yōu)勢(shì),發(fā)揮最佳性能;最后對(duì)比了AD_NTSM 制導(dǎo)律和比例導(dǎo)引(PN)制導(dǎo)律的性能。

    在仿真中目標(biāo)在慣性坐標(biāo)系下初始位置為(500.0 km,400.0 km,0 km),攔截導(dǎo)彈在慣性系下初始位置為(559.371 km,319.188 km,0.746 km),防御導(dǎo)彈在慣性系下的初始位置為(502.903 km,396.206 km,-1.980 km). 目標(biāo)的初始速度大小為7 071 m/s,初始彈道角為傾角θT0= -35.3°和偏角ψvT0=0°. 攔截導(dǎo)彈初始速度大小為3 000 m/s,初始彈道角為傾角θM0=44°和偏角ψvM0=0°. 防御導(dǎo)彈初始速度大小為7 071 m/s,初始彈道角為傾角θD0=-35°和偏角ψvD0= -2.1°. 目標(biāo)本身不機(jī)動(dòng),只受重力作用。攔截導(dǎo)彈和防御導(dǎo)彈的最大加速度均為3 g. 制導(dǎo)律參數(shù)選擇如下β1=β2=40,M1=M2=20,M3=M4=30,α1=α2=

    目標(biāo)需要通過通信設(shè)備將一些制導(dǎo)信息發(fā)送給防御導(dǎo)彈。目標(biāo)和防御導(dǎo)彈雖然速度很快(7 071 m/s),但是基本上是朝同一方向飛行,大小也大致相同,因此相對(duì)速度很小。而且防御導(dǎo)彈會(huì)在目標(biāo)的前下方5 ~7 km 的位置上伴飛,距離也很短。目標(biāo)和防御導(dǎo)彈之間沒有障礙物和大氣阻隔,通信條件良好。綜上所述制導(dǎo)律對(duì)通信設(shè)備的要求不高,通信設(shè)備比較容易實(shí)現(xiàn)小型化。另外由于通信只存在從目標(biāo)到防御導(dǎo)彈的單向通信,因此彈道導(dǎo)彈上只需要安裝發(fā)射機(jī),防御導(dǎo)彈上只需要安裝接收機(jī)。

    4.1 仿真算例

    在該仿真算例中防御導(dǎo)彈的脫靶量為0.28 m.攔截過程的彈道軌跡如圖2所示。從圖上看目標(biāo)和防御導(dǎo)彈的軌跡十分接近,實(shí)際上在整個(gè)制導(dǎo)階段中目標(biāo)和防御導(dǎo)彈的最小距離不小于5.17 km. 這個(gè)距離對(duì)于目標(biāo)和防御導(dǎo)彈來說是安全的。

    圖2 目標(biāo)、攔截導(dǎo)彈和防御導(dǎo)彈運(yùn)動(dòng)軌跡Fig.2 Motion trajectories of target,interceptor and defending missiles

    防御導(dǎo)彈加速度在彈體系上的分量如圖3所示。前0.5 s 軌控發(fā)動(dòng)機(jī)一直處于關(guān)閉狀態(tài),這是因?yàn)榍?.5 s 防御導(dǎo)彈和攔截導(dǎo)彈的相對(duì)距離超過了防御導(dǎo)彈導(dǎo)引頭的最大探測(cè)距離(本仿真假設(shè)導(dǎo)引頭最大探測(cè)距離為90 km)。在這段時(shí)間內(nèi)制導(dǎo)律不工作,因此不產(chǎn)生制導(dǎo)指令。

    圖3 防御導(dǎo)彈發(fā)動(dòng)機(jī)推力產(chǎn)生的加速度在彈體系下的分量Fig.3 Lateral accelerations of defending missile in body-axis frame

    圖4是防御導(dǎo)彈對(duì)攔截導(dǎo)彈的視線角速率。在大部分時(shí)間內(nèi)該視線角速率都在[- 0.05°/s,0.05°/s]范圍內(nèi),隨著防御導(dǎo)彈接近攔截導(dǎo)彈,視線角速率逐漸發(fā)散。

    圖4 防御導(dǎo)彈-攔截導(dǎo)彈視線角速率Fig.4 Defender-interceptor LOS angular velocities

    圖5描繪了目標(biāo)、攔截導(dǎo)彈和防御導(dǎo)彈三點(diǎn)共線的情況。整個(gè)制導(dǎo)階段防御導(dǎo)彈-攔截導(dǎo)彈和攔截導(dǎo)彈-目標(biāo)視線角之差的范圍分別為:視線高低角之差在[0°, - 0.2°]范圍內(nèi),視線方位角在[-2.2°,-1.5°]范圍內(nèi),共線情況良好。

    圖5 目標(biāo),攔截導(dǎo)彈和防御導(dǎo)彈三點(diǎn)共線情況Fig.5 Three point collinear of the three bodies

    4.2 μ* 對(duì)制導(dǎo)律性能的影響

    我們考察不同的μ*和不同的防御導(dǎo)彈機(jī)動(dòng)能力下整體制導(dǎo)律的性能。μ*從0 變化到1(代表階段1 制導(dǎo)律占據(jù)的時(shí)間越來越短),步長(zhǎng)為0.1. 防御導(dǎo)彈的機(jī)動(dòng)能力從2 g 變化到5 g,步長(zhǎng)為0.5 g.圖6給出了μ*在防御導(dǎo)彈不同機(jī)動(dòng)能力下對(duì)制導(dǎo)性能的影響。

    圖6 μ* 對(duì)制導(dǎo)律性能的影響Fig.6 Choice of μ* for interception phase dividing

    從圖6可以看出,制導(dǎo)全程使用階段1 制導(dǎo)律(μ*=0)的效果要比全程使用階段2 制導(dǎo)律(μ*=1)的效果好。其中原因在階段1 制導(dǎo)律設(shè)計(jì)中討論過,即由于在兩個(gè)視線坐標(biāo)系下的攔截導(dǎo)彈的加速度相互抵消而導(dǎo)致引入的加速度很小,利于防御導(dǎo)彈的攔截。這說明了視線三角制導(dǎo)策略的優(yōu)勢(shì)。在μ*較小時(shí),制導(dǎo)性能隨著μ*增加而提高;而當(dāng)μ*較大時(shí),制導(dǎo)性能隨著μ*增加而下降。這是因?yàn)樵讦?較小時(shí),在制導(dǎo)后期仍舊使用階段1 制導(dǎo)律,此時(shí)近似相對(duì)運(yùn)動(dòng)建模誤差快速變大,階段1 制導(dǎo)律性能下降。如果及時(shí)采用精確建模的階段2 制導(dǎo)律,總體制導(dǎo)性能會(huì)提高。隨著μ*繼續(xù)增大,階段2 制導(dǎo)律越來越早使用。當(dāng)μ*增加到一定程度(例如對(duì)于防御導(dǎo)彈機(jī)動(dòng)能力為2.0 g 或2.5 g 時(shí),μ*=0.4;防御導(dǎo)彈機(jī)動(dòng)能力為3.0 g 或3.5 g 時(shí),μ*=0.3),制導(dǎo)性能進(jìn)入拐點(diǎn),繼續(xù)增加μ*會(huì)造成階段1 制導(dǎo)律使用的時(shí)間過短,其制導(dǎo)優(yōu)勢(shì)難以充分體現(xiàn),從而導(dǎo)致總體制導(dǎo)性能開始下降。

    4.3 AD_NTSM 制導(dǎo)律和PN 制導(dǎo)律比較

    為了橫向比較AD_NTSM 制導(dǎo)律的性能,考察了在相同的防御導(dǎo)彈機(jī)動(dòng)能力下,AD_NTSM 制導(dǎo)律和PN 制導(dǎo)律的性能。碰撞殺傷場(chǎng)景需要的脫靶量很小,一般要求在0.3 ~0.5 m 之間,因此比較了AD_NTSM 制導(dǎo)律和PN 制導(dǎo)律脫靶量小于0.3 m和小于0.5 的概率以及平均脫靶量這3 個(gè)統(tǒng)計(jì)指標(biāo)。仿真考察了防御導(dǎo)彈機(jī)動(dòng)能力從2.0 g 變化到5.0 g,步長(zhǎng)為0.5 g 的7 種不同情況下,每種情況針對(duì)AD_NTSM 制導(dǎo)律和PN 制導(dǎo)律分別進(jìn)行了200 次Monte Carlo 仿真。

    AD_NTSM 制導(dǎo)律門限值取μ*=0.4. 其余仿真參數(shù)和仿真算例的設(shè)置一樣。在Monte Carlo 仿真中攔截導(dǎo)彈和防御導(dǎo)彈的初始位置都分別加上了一個(gè)均勻分布的隨機(jī)向量。這兩個(gè)隨機(jī)向量的2 范數(shù)的平均值分別為3.5 km 和5.0 km.

    仿真結(jié)果如表1所示。從表1中可以看出,在相同機(jī)動(dòng)能力下,AD_NTSM 制導(dǎo)律性能優(yōu)于PN 制導(dǎo)律。這是因?yàn)锳D_NTSM 制導(dǎo)律含有可以抵消攔截導(dǎo)彈機(jī)動(dòng)影響的變結(jié)構(gòu)項(xiàng)。特別的,當(dāng)防御導(dǎo)彈機(jī)動(dòng)能力較小時(shí),AD_NTSM 性能顯著優(yōu)于PN 制導(dǎo)律。除了上面的原因之外,還因?yàn)椴捎靡暰€三角制導(dǎo)策略的AD_NTSM 制導(dǎo)律在階段1 制導(dǎo)中當(dāng)防御導(dǎo)彈-攔截導(dǎo)彈視線系和攔截導(dǎo)彈-目標(biāo)視線系接近平行時(shí)(即目標(biāo),攔截導(dǎo)彈和防御導(dǎo)彈三者接近共線時(shí)),攔截導(dǎo)彈在兩個(gè)視線系下的加速度會(huì)互相抵消,導(dǎo)致引入到(17)式和(20)式中的攔截導(dǎo)彈加速度很小。由于制導(dǎo)律將攔截導(dǎo)彈的加速度當(dāng)做外部干擾處理,小的外部干擾有利于降低防御導(dǎo)彈需用過載。

    表1 AD_NTSM 制導(dǎo)律和PN 制導(dǎo)律性能比較Tab.1 Comparison of performances of AD_NTSM and PN guidance laws

    5 結(jié)論

    本文研究了大氣層外彈道導(dǎo)彈突防的主動(dòng)防御制導(dǎo)問題。將制導(dǎo)過程分為兩個(gè)階段。階段1 設(shè)計(jì)了基于視線三角制導(dǎo)策略的制導(dǎo)律,防御導(dǎo)彈控制自身位置處于目標(biāo)和攔截導(dǎo)彈的視線上。本文首次在理論上證明了該階段制導(dǎo)律能有效地降低防御導(dǎo)彈的需用過載。由于采用近似相對(duì)運(yùn)動(dòng)模型,簡(jiǎn)化了制導(dǎo)律的表達(dá)形式。但是在階段2 中近似建模誤差隨時(shí)間快速增大,因此設(shè)計(jì)了基于零化視線轉(zhuǎn)率的制導(dǎo)律,使用更為精確的相對(duì)運(yùn)動(dòng)模型。兩階段制導(dǎo)律均使用非奇異終端滑模控制方法設(shè)計(jì),保證防御導(dǎo)彈在有限時(shí)間內(nèi)可以將自身調(diào)整到目標(biāo)和攔截導(dǎo)彈的視線上并零化對(duì)攔截導(dǎo)彈的視線轉(zhuǎn)率。最后進(jìn)行了數(shù)值仿真,結(jié)果表明在整個(gè)制導(dǎo)階段防御導(dǎo)彈在較高的精度上達(dá)到了上述2 個(gè)目標(biāo):1)保持在目標(biāo)和攔截導(dǎo)彈的視線上;2)對(duì)攔截導(dǎo)彈的視線角速率控制趨于0. 和PN 律進(jìn)行對(duì)比可以看到,所提制導(dǎo)律在達(dá)到同樣的脫靶量精度情況下,需用過載更小。

    References)

    [1]馮小榮,樊秋林. 對(duì)外大氣層動(dòng)能殺傷攔截器EKV 的突防方法研究[J]. 紅外與激光工程,2006,35(增刊1):52 -55.FENG Xiao-rong,F(xiàn)AN Qiu-Lin. Research on ballistic missile exoatmospheric penertration against EKV[J]. Infrared and Laser Engineering,2006,35(S1):52 -55. (in Chinese)

    [2]Shima T. Optimal cooperative pursuit and evasion strategies against a homing missile[J]. Journal of Guidance,Control,and Dynamics,2011,34(2):414 -425.

    [3]Shaferman V,Shima T. Cooperative multiple-model adaptive guidance for an aircraft defending missile[J]. Journal of Guidance,Control,and Dynamics,2010,33(6):1801 -1813.

    [4]Prokopov O,Shima T. Linear quadratic optimal cooperative strategies for active aircraft protection[J]. Journal of Guidance,Control,and Dynamics,2013,36(3):753 -764.

    [5]Rubinsky S,Gutman S. Three-player pursuit and evasion conflict[J]. Journal of Guidance Control and Dynamics,2014,37(1):98 -110.

    [6]Perelman A,Shima T,Rusnak I. Cooperative differential games strategies for active aircraft protection from a homing missile[J].Journal of Guidance,Control,and Dynamics,2011,34(3):761 -773.

    [7]Ratnoo A,Shima T. Line-of-sight interceptor guidance for defending an aircraft[J]. Journal of Guidance,Control,and Dynamics,2011,34(2):522 -532.

    [8]Yamasaki T,Balakrishnan S N,Takano H. Geometrical approachbased defense-missile intercept guidance for aircraft protection against missile attack[J]. Journal of Aerospace Engineering,2012,226(8):1014 -1028.

    [9]葉繼坤,雷虎民,李炯,等. 基于角度關(guān)系的幾何制導(dǎo)律研究[J]. 系統(tǒng)工程與電子技術(shù),2011,33(12):2721 -2725.YE Ji-kun,LEI Hu-min,LI Jiong,et al. Research on geometric guidance law based on angle relation[J]. Systems Engineering and Electronics,2011,33(12):2721 -2725. (in Chinese)

    [10]Brierly S D,Longchamp R. Application of sliding mode control to air-air interception problem[J]. IEEE Transactions on Aerospace and Electronic Systems,1990,26(2):306 -325.

    [11]Zhou D,Mu C D,Xu W L. Adaptive sliding-mode guidance of a homing missile[J]. Journal of Guidance,Control,and Dynamics,1999,22(4):589 -594.

    [12]Moon J,Kim K,Kim Y. Design of missile guidance law via variable structure control[J]. Journal of Guidance,Control,and Dynamics,2001,24(4):659 -664.

    [13]花文華,陳興林,宋超. 主動(dòng)防御協(xié)同自適應(yīng)滑模制導(dǎo)律[J]. 飛行力學(xué),2012,30(1):66 -70.HUA Wen-hua,CHEN Xing-lin,SONG Chao.Cooperative guidance law for active protection based on adaptive sliding-mode control[J]. Flight Dynamics,2012,30(1):66 -70. (in Chinese)

    [14]Haimo V T. Finite time controllers SIAM[J]. Journal on Control and Optimization,1986,24(4):760 -770.

    [15]Zhou D,Sun S,Teo K L. Guidance Laws with finite time convergence[J]. Journal of Guidance,Control,and Dynamics,2009,32(6):1838 -1846.

    [16]Zhou D,Qu P P,Sun S. A guidance law with terminal impact angle constraint accounting for missile autopilot[J]. Journal of Dynamic Systems,Measurement,and Control,2013,135(5):051009.

    [17]Kumar S R,Rao S,Ghose D. Sliding mode guidance and control for all-aspect interceptor with terminal angle constraints[J].Journal of Guidance,Control,and Dynamics,2012,35(4):1230 -1246.

    [18]Zhang Y,Sun M,Chen Z. Finite-time convergent guidance law with impact angle constraint based on sliding-mode control[J].Nonlinear Dynamics,2012,70(1):619 -625.

    [19]Feng Y,Yu X,Han F L. On nonsingular terminal sliding-mode control of nonlinear systems[J]. Automatica,2013,49(6):1715 -1722.

    [20]Yu S H,Yu X H,Shirinzadeh B,et al. Continuous finite-time control for robotic manipulators with terminal sliding mode[J].Automatica,2005,41(11):1957 -1964.

    [21]Khoo S Y,Xie L H,Man Z H. Robust finite-time consensus tracking algorithm for multirobot systems[J]. IEEE ASME Transactions on Mechatronics,2009,14(2):219 -228.

    [22]馮勇,鮑晟,余星火. 非奇異終端滑??刂葡到y(tǒng)的設(shè)計(jì)方法[J]. 控制與決策,2002,17(2):194 -198.FENG Yong,BAO Sheng,YU Xing-huo. Design method of nonsingular terminal sliding mode control systems[J]. Control and Decision,2002,17(2):194 -198.(in Chinese)

    [23]Kumar S R,Rao S,Ghose D. Nonsingular terminal sliding mode guidance with impact angle constraints[J]. Journal of Guidance,Control,and Dynamics,2014,37(4):1 -17.

    [24]Wang X,Wang J. Partial integrated missile guidance and control with finite time convergence[J]. Journal of Guidance,Control,and Dynamics,2013,36(5):1399 -1409.

    [25]孫勝. 有限時(shí)間收斂尋的導(dǎo)引律[D]. 哈爾濱:哈爾濱工業(yè)大學(xué),2010.SUN Sheng. Guidance laws with finite time convergence for homing missiles[D]. Harbin:Harbin Institute of Technology,2010.(in Chinese)

    [26]Bhat S P,Bernstein D S. Finite-time stability of continuous autonomous systems[J]. SIAM Journal on Control and Optimization,2000,38(3):751 -766.

    猜你喜歡
    攔截導(dǎo)彈視線制導(dǎo)
    要去就去視線盡頭的山
    你吸引了我的視線
    基于MPSC和CPN制導(dǎo)方法的協(xié)同制導(dǎo)律
    基于在線軌跡迭代的自適應(yīng)再入制導(dǎo)
    帶有攻擊角約束的無抖振滑模制導(dǎo)律設(shè)計(jì)
    當(dāng)代視線
    復(fù)合制導(dǎo)方式確保精確入軌
    太空探索(2014年1期)2014-07-10 13:41:49
    雨天戴偏光太陽鏡 視線更清晰
    99精国产麻豆久久婷婷| 国产xxxxx性猛交| 看免费av毛片| 国产成人精品久久二区二区91| 中文亚洲av片在线观看爽| 国产成人一区二区三区免费视频网站| 日韩欧美一区二区三区在线观看| 黄色毛片三级朝国网站| a级毛片在线看网站| 三上悠亚av全集在线观看| 制服人妻中文乱码| 免费高清视频大片| 色综合婷婷激情| 国产亚洲欧美98| 天天添夜夜摸| 欧美日本亚洲视频在线播放| 丝袜美足系列| 黄色片一级片一级黄色片| 国产一区二区激情短视频| 大型黄色视频在线免费观看| 91精品国产国语对白视频| 在线观看日韩欧美| 欧美黄色片欧美黄色片| 欧美亚洲日本最大视频资源| 久久人妻av系列| 午夜老司机福利片| 如日韩欧美国产精品一区二区三区| 天堂俺去俺来也www色官网| 国产成人啪精品午夜网站| 麻豆av在线久日| 日韩成人在线观看一区二区三区| 久久热在线av| 色婷婷av一区二区三区视频| 一个人观看的视频www高清免费观看 | 亚洲国产欧美日韩在线播放| 成人亚洲精品av一区二区 | 人妻久久中文字幕网| 日本五十路高清| 国产一区二区激情短视频| 精品国产乱子伦一区二区三区| 一边摸一边抽搐一进一小说| videosex国产| 精品一区二区三卡| 欧美日韩黄片免| 99国产精品一区二区三区| 91老司机精品| 日韩一卡2卡3卡4卡2021年| 免费女性裸体啪啪无遮挡网站| 男女下面进入的视频免费午夜 | 久久精品91蜜桃| 黑人操中国人逼视频| 精品无人区乱码1区二区| av在线天堂中文字幕 | 亚洲精品国产精品久久久不卡| 乱人伦中国视频| 欧美日韩av久久| 欧美大码av| 亚洲成a人片在线一区二区| 又大又爽又粗| 国产三级在线视频| 国产高清国产精品国产三级| 宅男免费午夜| av免费在线观看网站| 国产欧美日韩一区二区精品| 久久精品国产清高在天天线| 国内久久婷婷六月综合欲色啪| 午夜福利在线观看吧| 91麻豆精品激情在线观看国产 | 国产精品九九99| av在线播放免费不卡| 亚洲精品国产区一区二| 黄色毛片三级朝国网站| e午夜精品久久久久久久| 日韩高清综合在线| 日韩视频一区二区在线观看| 男女做爰动态图高潮gif福利片 | 黄色视频,在线免费观看| 久久国产乱子伦精品免费另类| 激情在线观看视频在线高清| 午夜久久久在线观看| 久99久视频精品免费| 琪琪午夜伦伦电影理论片6080| 亚洲全国av大片| 黄色成人免费大全| 成人影院久久| 久9热在线精品视频| 亚洲性夜色夜夜综合| 天天影视国产精品| 国产一区二区三区视频了| 亚洲欧洲精品一区二区精品久久久| 欧美精品一区二区免费开放| 欧美日韩亚洲综合一区二区三区_| av在线播放免费不卡| 国产精品98久久久久久宅男小说| 国内久久婷婷六月综合欲色啪| 国产av又大| 久久久久久久精品吃奶| 亚洲久久久国产精品| 女生性感内裤真人,穿戴方法视频| 亚洲午夜理论影院| 我的亚洲天堂| 国产精品 欧美亚洲| 国产欧美日韩综合在线一区二区| 在线看a的网站| 欧美黑人欧美精品刺激| 国产又爽黄色视频| 欧美亚洲日本最大视频资源| 9热在线视频观看99| 热re99久久国产66热| 两个人看的免费小视频| 精品欧美一区二区三区在线| 久久久国产成人免费| 欧美精品一区二区免费开放| 电影成人av| 欧美激情极品国产一区二区三区| 天天影视国产精品| 久久久久久免费高清国产稀缺| 少妇的丰满在线观看| 中文字幕另类日韩欧美亚洲嫩草| 丁香欧美五月| 国产无遮挡羞羞视频在线观看| 精品高清国产在线一区| 多毛熟女@视频| 在线观看免费午夜福利视频| 久久天躁狠狠躁夜夜2o2o| 欧美激情极品国产一区二区三区| 亚洲精品在线观看二区| 国产亚洲精品久久久久久毛片| 看免费av毛片| 狂野欧美激情性xxxx| 午夜精品久久久久久毛片777| 免费观看人在逋| 亚洲成人国产一区在线观看| 日韩欧美国产一区二区入口| 999久久久精品免费观看国产| 免费观看精品视频网站| 欧美一级毛片孕妇| 男女下面插进去视频免费观看| 久久精品影院6| 1024香蕉在线观看| av欧美777| 97碰自拍视频| av免费在线观看网站| netflix在线观看网站| 欧美午夜高清在线| 精品少妇一区二区三区视频日本电影| 色播在线永久视频| 丰满迷人的少妇在线观看| 欧美大码av| 可以免费在线观看a视频的电影网站| 99精品久久久久人妻精品| 国产黄色免费在线视频| 国产成人影院久久av| 91成年电影在线观看| 国产不卡一卡二| avwww免费| 国产又色又爽无遮挡免费看| 一级a爱片免费观看的视频| 欧美日韩亚洲高清精品| 亚洲欧美精品综合一区二区三区| 在线播放国产精品三级| 热99re8久久精品国产| 精品久久蜜臀av无| 亚洲精品av麻豆狂野| 日韩免费高清中文字幕av| 日韩中文字幕欧美一区二区| 日日摸夜夜添夜夜添小说| 多毛熟女@视频| 亚洲欧美精品综合久久99| 别揉我奶头~嗯~啊~动态视频| а√天堂www在线а√下载| 欧美激情 高清一区二区三区| 村上凉子中文字幕在线| 日本五十路高清| 最好的美女福利视频网| 久久婷婷成人综合色麻豆| 男女午夜视频在线观看| 999精品在线视频| 亚洲一区二区三区色噜噜 | 亚洲欧美激情在线| 女性被躁到高潮视频| 久久精品亚洲精品国产色婷小说| 男男h啪啪无遮挡| 可以在线观看毛片的网站| 黄色怎么调成土黄色| 亚洲情色 制服丝袜| av天堂久久9| 高清在线国产一区| 黄色丝袜av网址大全| 久久影院123| 欧美中文日本在线观看视频| 国产精品美女特级片免费视频播放器 | 亚洲va日本ⅴa欧美va伊人久久| 精品卡一卡二卡四卡免费| 在线国产一区二区在线| 50天的宝宝边吃奶边哭怎么回事| 国产av一区在线观看免费| 在线播放国产精品三级| 精品欧美一区二区三区在线| 国产亚洲欧美98| 中文字幕最新亚洲高清| 午夜久久久在线观看| 国产片内射在线| 国产片内射在线| 午夜成年电影在线免费观看| 国产人伦9x9x在线观看| 高清av免费在线| 亚洲国产中文字幕在线视频| svipshipincom国产片| 丰满人妻熟妇乱又伦精品不卡| 首页视频小说图片口味搜索| 久久久久九九精品影院| 交换朋友夫妻互换小说| 啪啪无遮挡十八禁网站| 亚洲aⅴ乱码一区二区在线播放 | 97碰自拍视频| 国产野战对白在线观看| 在线观看一区二区三区激情| 欧美av亚洲av综合av国产av| 亚洲美女黄片视频| 91成人精品电影| 国产野战对白在线观看| 午夜影院日韩av| 丁香欧美五月| 亚洲精品国产精品久久久不卡| 天天添夜夜摸| 一夜夜www| 久久精品成人免费网站| 成年人黄色毛片网站| aaaaa片日本免费| 欧美丝袜亚洲另类 | 久99久视频精品免费| 宅男免费午夜| 亚洲少妇的诱惑av| 亚洲专区中文字幕在线| 国产亚洲精品久久久久久毛片| 在线十欧美十亚洲十日本专区| 制服人妻中文乱码| 正在播放国产对白刺激| 人人澡人人妻人| 亚洲在线自拍视频| 亚洲中文日韩欧美视频| 91九色精品人成在线观看| 中文字幕精品免费在线观看视频| 国产精品1区2区在线观看.| 91字幕亚洲| 97碰自拍视频| 午夜视频精品福利| 嫁个100分男人电影在线观看| 99久久精品国产亚洲精品| 久久热在线av| 国产又爽黄色视频| 纯流量卡能插随身wifi吗| 国产男靠女视频免费网站| 人人妻,人人澡人人爽秒播| 日韩精品免费视频一区二区三区| 亚洲av成人av| 黑丝袜美女国产一区| 免费人成视频x8x8入口观看| 一级作爱视频免费观看| 久久久国产一区二区| 精品久久久久久久久久免费视频 | 亚洲专区字幕在线| 777久久人妻少妇嫩草av网站| 午夜视频精品福利| 国产精品香港三级国产av潘金莲| 国产一区二区三区综合在线观看| 丝袜美腿诱惑在线| 丝袜人妻中文字幕| 午夜老司机福利片| 国内久久婷婷六月综合欲色啪| 人妻丰满熟妇av一区二区三区| 欧美乱码精品一区二区三区| 变态另类成人亚洲欧美熟女 | 国产精品国产高清国产av| 国产无遮挡羞羞视频在线观看| 美女大奶头视频| 久久久久久久午夜电影 | 50天的宝宝边吃奶边哭怎么回事| av免费在线观看网站| av视频免费观看在线观看| 中文字幕人妻丝袜制服| 高清毛片免费观看视频网站 | 9色porny在线观看| 日韩精品中文字幕看吧| 一进一出抽搐gif免费好疼 | tocl精华| 91麻豆av在线| 精品一区二区三区av网在线观看| 国产伦一二天堂av在线观看| 激情视频va一区二区三区| 国产精品影院久久| 人人妻人人爽人人添夜夜欢视频| 亚洲av电影在线进入| 日韩大码丰满熟妇| 国产97色在线日韩免费| 9色porny在线观看| 欧美人与性动交α欧美软件| 女人被躁到高潮嗷嗷叫费观| 一a级毛片在线观看| 成在线人永久免费视频| 超色免费av| 中出人妻视频一区二区| 99国产精品99久久久久| 一级a爱视频在线免费观看| 精品电影一区二区在线| 国产91精品成人一区二区三区| 91成年电影在线观看| 在线观看午夜福利视频| 97超级碰碰碰精品色视频在线观看| 国内久久婷婷六月综合欲色啪| 久久久久久久久中文| 精品福利观看| 欧美+亚洲+日韩+国产| 亚洲中文av在线| 久久久久九九精品影院| 狂野欧美激情性xxxx| 窝窝影院91人妻| 不卡av一区二区三区| www.自偷自拍.com| 亚洲久久久国产精品| а√天堂www在线а√下载| 后天国语完整版免费观看| 老司机福利观看| 91麻豆av在线| 黄片播放在线免费| 亚洲欧美精品综合一区二区三区| 大型黄色视频在线免费观看| 精品乱码久久久久久99久播| 亚洲一码二码三码区别大吗| 成人国语在线视频| 99精品在免费线老司机午夜| 日韩中文字幕欧美一区二区| 国产高清激情床上av| 一夜夜www| 午夜成年电影在线免费观看| 久久热在线av| 免费日韩欧美在线观看| 久久精品91蜜桃| 少妇被粗大的猛进出69影院| 日韩视频一区二区在线观看| 操出白浆在线播放| 91老司机精品| 国产亚洲精品久久久久久毛片| 亚洲第一av免费看| av在线天堂中文字幕 | 多毛熟女@视频| 亚洲午夜精品一区,二区,三区| 欧美日本亚洲视频在线播放| 69av精品久久久久久| 成人18禁高潮啪啪吃奶动态图| 国产99白浆流出| 天天添夜夜摸| 91麻豆av在线| 欧美性长视频在线观看| 亚洲成人精品中文字幕电影 | 国产精品久久视频播放| 成人国语在线视频| 在线观看舔阴道视频| www.自偷自拍.com| 亚洲精品美女久久av网站| 丝袜人妻中文字幕| 国内久久婷婷六月综合欲色啪| 女性被躁到高潮视频| 国产xxxxx性猛交| 久久中文字幕一级| 热99re8久久精品国产| 一级毛片高清免费大全| 大型黄色视频在线免费观看| 天堂动漫精品| 一边摸一边抽搐一进一小说| 欧洲精品卡2卡3卡4卡5卡区| 欧美乱色亚洲激情| 成人三级黄色视频| 欧美日韩亚洲综合一区二区三区_| 很黄的视频免费| 亚洲中文日韩欧美视频| 亚洲av成人不卡在线观看播放网| 女生性感内裤真人,穿戴方法视频| 欧美人与性动交α欧美精品济南到| 久久人妻av系列| 亚洲 欧美 日韩 在线 免费| 国产精品一区二区三区四区久久 | 校园春色视频在线观看| 午夜日韩欧美国产| 欧美一区二区精品小视频在线| www.999成人在线观看| 嫁个100分男人电影在线观看| www.熟女人妻精品国产| 国产乱人伦免费视频| 国产一区在线观看成人免费| 天堂中文最新版在线下载| 中文字幕色久视频| 欧美精品亚洲一区二区| 91精品三级在线观看| 欧美色视频一区免费| 国产精品永久免费网站| 少妇 在线观看| 国产一卡二卡三卡精品| 在线看a的网站| 欧美日韩亚洲国产一区二区在线观看| 国产亚洲精品一区二区www| 亚洲欧美日韩无卡精品| 男女之事视频高清在线观看| 亚洲免费av在线视频| 免费观看精品视频网站| 热re99久久精品国产66热6| 欧美黑人欧美精品刺激| 69精品国产乱码久久久| 啦啦啦在线免费观看视频4| 少妇粗大呻吟视频| 国产黄a三级三级三级人| 欧美日韩一级在线毛片| 久久国产亚洲av麻豆专区| av中文乱码字幕在线| 在线天堂中文资源库| 无人区码免费观看不卡| 国产亚洲精品综合一区在线观看 | 激情在线观看视频在线高清| 精品人妻在线不人妻| 老鸭窝网址在线观看| 国产高清videossex| 在线观看www视频免费| 成年人免费黄色播放视频| 精品久久久久久电影网| 757午夜福利合集在线观看| 精品人妻在线不人妻| 制服诱惑二区| 色综合婷婷激情| avwww免费| 欧美精品一区二区免费开放| 91在线观看av| av欧美777| 国产av又大| netflix在线观看网站| 在线观看免费视频网站a站| 久久性视频一级片| 午夜日韩欧美国产| 国产精品香港三级国产av潘金莲| 一边摸一边抽搐一进一小说| 午夜免费鲁丝| 日韩成人在线观看一区二区三区| 老鸭窝网址在线观看| 黄色a级毛片大全视频| 日韩欧美国产一区二区入口| 国产精品一区二区免费欧美| av免费在线观看网站| 亚洲视频免费观看视频| 一进一出抽搐动态| 激情在线观看视频在线高清| 免费在线观看完整版高清| 三级毛片av免费| 午夜福利,免费看| x7x7x7水蜜桃| 9热在线视频观看99| 亚洲精品美女久久久久99蜜臀| 日本撒尿小便嘘嘘汇集6| 欧美不卡视频在线免费观看 | 另类亚洲欧美激情| 真人一进一出gif抽搐免费| 香蕉久久夜色| 亚洲av美国av| 国产区一区二久久| 夜夜看夜夜爽夜夜摸 | 午夜免费激情av| 国产亚洲精品综合一区在线观看 | 色在线成人网| 男人的好看免费观看在线视频 | 色综合婷婷激情| 国产精品香港三级国产av潘金莲| 久久九九热精品免费| 中文字幕精品免费在线观看视频| 国产精品国产av在线观看| 免费在线观看黄色视频的| 欧美精品啪啪一区二区三区| 午夜福利,免费看| 久久精品国产99精品国产亚洲性色 | 美女高潮到喷水免费观看| 成熟少妇高潮喷水视频| 午夜视频精品福利| 性色av乱码一区二区三区2| 成人18禁高潮啪啪吃奶动态图| 国产精品影院久久| www.www免费av| netflix在线观看网站| 亚洲,欧美精品.| 精品欧美一区二区三区在线| 久久久久久人人人人人| 亚洲欧美一区二区三区黑人| 黄色怎么调成土黄色| 啦啦啦免费观看视频1| 99香蕉大伊视频| 国产乱人伦免费视频| 欧美日本中文国产一区发布| 婷婷丁香在线五月| 老熟妇仑乱视频hdxx| 亚洲国产欧美一区二区综合| 免费一级毛片在线播放高清视频 | 国产av在哪里看| 黄片大片在线免费观看| 两性午夜刺激爽爽歪歪视频在线观看 | 成人精品一区二区免费| 一区二区日韩欧美中文字幕| 日本wwww免费看| 久久精品91蜜桃| 狠狠狠狠99中文字幕| 窝窝影院91人妻| 午夜福利一区二区在线看| 亚洲男人的天堂狠狠| 女同久久另类99精品国产91| 少妇裸体淫交视频免费看高清 | 亚洲五月天丁香| 精品高清国产在线一区| 黄片小视频在线播放| 精品久久蜜臀av无| 91字幕亚洲| 成人手机av| 久久精品亚洲熟妇少妇任你| 长腿黑丝高跟| 首页视频小说图片口味搜索| 国产精品久久久久成人av| 麻豆久久精品国产亚洲av | 日韩欧美免费精品| 人人妻人人澡人人看| 免费久久久久久久精品成人欧美视频| 黄片播放在线免费| 精品一区二区三区视频在线观看免费 | 亚洲精品一区av在线观看| 中亚洲国语对白在线视频| 欧美乱码精品一区二区三区| 日本黄色日本黄色录像| 国产av一区在线观看免费| 亚洲国产精品999在线| av福利片在线| 亚洲熟女毛片儿| a级毛片在线看网站| 超碰成人久久| 成人亚洲精品一区在线观看| 怎么达到女性高潮| 亚洲片人在线观看| 夜夜夜夜夜久久久久| 纯流量卡能插随身wifi吗| 国产成人影院久久av| 欧美日韩av久久| 最近最新中文字幕大全免费视频| 看片在线看免费视频| 一二三四社区在线视频社区8| svipshipincom国产片| 村上凉子中文字幕在线| 久久这里只有精品19| 精品电影一区二区在线| 两性夫妻黄色片| 久久国产精品影院| 人成视频在线观看免费观看| 咕卡用的链子| 亚洲av电影在线进入| 一级毛片高清免费大全| 男女床上黄色一级片免费看| 久久人妻av系列| 久久天堂一区二区三区四区| 国产精品久久电影中文字幕| 色在线成人网| 精品久久久久久电影网| 伊人久久大香线蕉亚洲五| 国产欧美日韩综合在线一区二区| 黄网站色视频无遮挡免费观看| 国产精品 国内视频| 两性午夜刺激爽爽歪歪视频在线观看 | 午夜91福利影院| 精品一品国产午夜福利视频| 99在线视频只有这里精品首页| 国产精品亚洲一级av第二区| 亚洲一区二区三区欧美精品| 欧美日韩av久久| 交换朋友夫妻互换小说| 老司机在亚洲福利影院| 亚洲av熟女| 欧美不卡视频在线免费观看 | 国产精品日韩av在线免费观看 | 级片在线观看| 亚洲伊人色综图| 在线观看日韩欧美| 亚洲成人久久性| 国产熟女午夜一区二区三区| 亚洲三区欧美一区| 一区二区三区精品91| 亚洲中文字幕日韩| 91国产中文字幕| 超碰97精品在线观看| 欧美成狂野欧美在线观看| 久久精品91蜜桃| 久久午夜亚洲精品久久| 亚洲精品在线美女| 欧美色视频一区免费| 欧美日韩亚洲高清精品| 在线观看免费午夜福利视频| 一进一出抽搐gif免费好疼 | 97人妻天天添夜夜摸| 免费搜索国产男女视频| 国产精品免费视频内射| 久久精品成人免费网站| 国产亚洲av高清不卡| 桃红色精品国产亚洲av| 日本黄色视频三级网站网址| 成年人黄色毛片网站| 999精品在线视频| 中文字幕人妻熟女乱码| 好男人电影高清在线观看| 精品无人区乱码1区二区| av中文乱码字幕在线| 久久久国产成人精品二区 | 黄色片一级片一级黄色片| 欧美+亚洲+日韩+国产| 亚洲欧美激情综合另类| 91成人精品电影|