999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

紅外空空導彈抗干擾性能評估指標體系研究

2020-12-14 09:30:46張喜濤白曉東王煒強
紅外技術 2020年11期

張喜濤,白曉東,閆 琳,王煒強

〈制導與對抗〉

紅外空空導彈抗干擾性能評估指標體系研究

張喜濤1,2,白曉東1,2,閆 琳1,2,王煒強1,2

(1. 中國空空導彈研究院,河南 洛陽 471009;2. 航空制導武器航空科技重點實驗室,河南 洛陽 471009)

評估紅外空空導彈抗干擾能力強弱的傳統指標較為單一,多是利用綜合抗干擾概率進行判定。針對該問題,建立了包含導彈總體、制導系統、導引頭3個層次的抗干擾性能評估指標體系,實現抗干擾性能評估指標的分解,提高了對制導系統和導引頭抗干擾性能的評估能力,提升了利用導引頭、制導系統評估結果對導彈總體抗干擾性能進行預估的能力。

紅外導彈;抗干擾;指標體系;指標分解;性能評估;性能預估

0 引言

現代戰場環境中,沒有人工干擾的凈空作戰條件已經不復存在,對敵方的紅外制導武器實施人工干擾已經成為一種常規的作戰手段,為了能夠適應現代化作戰環境,紅外制導武器在研制過程中,必須準確評估其在復雜戰場環境中的抗干擾性能[1]。國內抗干擾性能評估主要集中在導彈總體層面,往往以抗干擾概率來評估導彈總體、制導系統乃至導引頭的抗干擾能力[2-3],對抗干擾評估指標體系缺乏系統研究,從而無法在研制及鑒定過程中對抗干擾性能進行充分評估。由此,一方面使得軍方對導彈抗干擾能力認識不夠充分,影響導彈實際作戰效能的發揮;另一方面,使研制方對導彈的抗干擾性能掌控不深入,嚴重制約了研制進度。

上海機電工程研究所唐善軍等[4],提出識別概率、最大最小能量壓制、視線角速度誤差因子、干擾結束時間占比因子和干擾占空比等6項指標系統,對抗干擾過程中目標的識別能力和跟蹤能力進行評價,但其無法區分導彈不同層次的評估。南京航空航天大學韓培駿等[5],詳細分析了紅外光學系統、調制器、探測器、圖像處理系統、陀螺伺服系統的性能參數,將指標設計為發射干擾時的彈目距離、抗干擾檢測時間、抗干擾軟件識別概率和干擾/目標能量比,無法區分對抗條件與抗干擾性能之間的關系,其采用層次分析法建立了紅外成像導引系統抗干擾性能計算模型,該模型具有主觀性,導彈總體和制導系統無有效評估途徑。南京航空航天大學許友平等[6]綜合導引頭固有性能指標和性能改善指標,建立了導引頭抗干擾評估指標集合,無法對導彈抗干擾性能進行多層次的評估。

本文在對抗環境影響機理分析的基礎上,從紅外成像空空導彈總體、制導系統、導引頭3個層次詳細分析了干擾的影響機理,闡明了傳統抗干擾性能評估指標的不足,提出了導彈總體、制導系統和導引頭3個層次的抗干擾性能評估指標體系,描述了抗干擾性能評估指標的分解過程,為研制過程中抗干擾改進設計及抗干擾性能快速評估提供指導。

1 紅外對抗過程影響機理分析

紅外制導空空導彈常用于近距格斗,其作戰任務是對目標進行殺傷,作戰過程中面臨人工干擾的對抗場景。目標與導彈博弈對抗過程中,遵循以下干擾投放原則:導彈未發射時,干擾以預防為主進行稀疏投放,主要影響導引頭正確截獲目標;導彈發射后,干擾以誘偏為目的進行密集投放,主要影響導引頭對目標的正確測量,進而對導彈制導系統產生影響,降低導彈的殺傷能力。干擾的兩種投放策略分別影響導引頭對目標的截獲和測量,作用機理差異較大。本文按照對抗環境、導彈總體、制導系統及導引頭多個層次對紅外對抗機理進行研究,提出用于近距格斗導彈截獲目標后的紅外空空導彈抗干擾性能評估指標體系。

1.1 復雜對抗環境影響分析

目標、干擾在導彈攻擊過程中的相互關系構成了復雜多變的作戰對抗環境。在導彈攻擊目標過程中,目標機動使得目標的輻射特性和運動特性發生持續改變,在干擾的遮擋情況下,導引頭無法完成對目標的正確測量,進而影響導彈的殺傷能力。干擾投放模式主要包含干擾投放方向、組數、組間隔、每組彈數、彈間隔等要素,上述要素不僅能夠影響干擾在空間維度上的分布,同時也能夠影響干擾在時間維度上的分布。不同的干擾投放時機對導彈的擾動效果不同,近距投放干擾,對導彈跟蹤造成的擾動更大,更容易導致導彈丟失目標或脫靶。

對抗態勢是對目標、干擾與導彈在相對位置、速度等方面的描述,包含了對抗環境的當前狀態以及過程發展趨勢,對抗態勢的構成要素包括發射時刻的彈目距離、目標速度、導彈速度、進入角、離軸角、目標高度、導彈高度。發射時刻導彈和目標的相對位置、速度狀態決定了導彈發射后攻擊目標的動態變化過程,也直接影響干擾和導彈的對抗效果。

依據空戰復雜對抗環境的影響分析,抽取作戰環境的重要影響因素,將其定量化描述,并形成若干實例樣本,是空空導彈抗干擾性能評估的輸入和基礎。實例樣本由導彈發射條件、攻擊目標機動方式、干擾投放模式3個維度確定。在設計樣本時,依據攻擊目標和對抗態勢,確定導彈發射條件,在發射條件上疊加干擾投放模式形成干擾對抗條件,選取真實戰場中發生概率較高的樣本條件,最終形成戰術級紅外對抗作戰樣本空間。

1.2 干擾對導彈總體性能影響分析

紅外空空導彈的作戰任務是對目標進行殺傷,主要通過作戰條件下的殺傷概率進行評估。按照全概率公式,導彈的殺傷概率為[7]:

式中:(,)為遇靶平面制導誤差的分布;1(|,)為給定制導誤差(,)時,引信引爆點沿軸向的分布,2(,)為引信啟動概率;(,,)為戰斗部在點(,,)起爆條件下的分布。

制導誤差是在導彈外部和內部各種擾動作用下,形成的導彈實際彈道相對于理想彈道的誤差。導彈內部擾動是非均勻分布的,制導誤差的分布也具有一定的偏差;在外部干擾的作用下,制導誤差被進一步放大,該放大過程與干擾的投放模式有關,是非線性的放大過程[8]。

引信的啟動概率主要受引信的探測性能影響,與引信的靈敏度和目標的有效投影面積有關,因此,引信的啟動概率與目標的尺寸和目標的姿態關系密切相關,此外,干擾在引信工作范圍內投放,其核心區域也有一定的概率觸發引信啟動信號,對引信的啟動概率有一定的影響。在彈目交會階段,引信依據交會方向和交會速度輸出引爆信號[9],以實現對目標的最大殺傷。引信引爆點分布的主要影響因素為彈目交會方向和交會速度。

戰斗部起爆后的分布與交會速度關系密切[10],交會速度越大,破片動態飛散區越前傾;交會速度與彈軸夾角越大,破片動態飛散區不對稱性越嚴重。

導彈的殺傷性能不僅與導彈性能有關,也與飛機的易損性有關。空中目標往往有多個要害點,各個要害點有不同的易損性。如果認為殺傷目標上任意要害部位是獨立的,都將以一定的概率構成對整體目標的殺傷,則導彈的殺傷概率可改寫為:

式中:()為導彈對目標第個要害部位的殺傷概率;0()為目標第個要害部位被命中后的易損概率;為目標的要害部位個數。

綜上,假設干擾觸發引信啟動的概率較小,干擾主要通過放大制導誤差的分布,影響導彈的殺傷概率。假設制導誤差均勻分布,干擾對制導誤差分布的放大系數為,則通過公式(1)、(2)可得干擾對導彈的殺傷概率為:

2=×1(3)

如果考慮干擾對引起啟動概率的影響及制導誤差的非均勻分布,則需要建立制導系統、引信、戰斗部及目標的仿真模型,利用對抗樣本空間進行干擾對抗下的導彈殺傷概率仿真。

1.3 干擾對制導系統性能影響分析

制導系統主要由導引頭和飛控系統組成。導引頭負責測量導彈與目標之間的相對運動信息,將視線角速度提供給飛控系統,飛控系統按照一定的導引規律輸出控制指令,控制導彈修正偏差,準確飛向目標。假設目標形心為唯一要害點,則制導誤差可以用脫靶量代替。

式中:為比例導引系數;c為彈目相對速度;為視線角;F為遇靶時刻;為當前時間。

在不考慮導彈過載限制的條件下,對導彈過載公式進行積分,并將視線角代入,可得脫靶量、比例導引系數及時間的關系為:

式中:為常數;為中間變量。

由上式可知,(F)為0,表示制導系統在機動能力足夠,信號無延時的理想狀態時,脫靶量必將衰減為零。

導彈引入延時和擾動[11],假設延時時間常數為,利用伴隨分析方法[12],計算時刻引入過載為的擾動,則脫靶量、飛行時間及時間常數的關系為:

圖1為比例導引系數等于4時,擾動幅度相同,時間常數不同時,擾動時機與脫靶量的歸一化曲線。

從仿真曲線可以看出,在剩余飛行時間是制導時間常數5倍以上進行擾動時,脫靶量最終將收斂為0;制導系統時間常數越大,脫靶量發散區間越大;脫靶量對擾動的時機較為敏感,當剩余飛行時間等于制導時間常數時進行擾動,最終將產生最大的脫靶量。因此,縮小制導系統時間常數,可以在一定程度上降低脫靶量。

圖2為比例導引系數等于4時,時間常數相同,擾動幅度不同時,擾動時機與脫靶量的歸一化曲線。從仿真曲線可以看出,擾動幅度越大,脫靶量越大。

圖2 脫靶量歸一化曲線

考察0時刻擾動引起的脫靶量變化,將公式(6)記作:

將上式進行泰勒展開,得到:

式中:0為常數;G為函數相對于變量的導數,在導彈時間常數變化量D、比例導引系數變化量D、干擾投放時刻變化量D0為0,干擾投放時彈目接近速度變化量Dc不大的條件下,令D0=D=D=Dc=0,得到時刻干擾引起視線角速度變化與脫靶量的關系為:

則公式(3)中的干擾對制導誤差分布的放大系數可表示為:

=1+()

()為變量的均值。假設脫靶量服從瑞利分布,即:

式中:2為脫靶量的均方,當殺傷半徑為,脫靶量滿足殺傷半徑要求的概率為:

綜上,無干擾條件下,脫靶量滿足殺傷半徑要求的概率為:

()為變量的方差。干擾條件下,脫靶量滿足殺傷半徑要求的概率為:

抗干擾成功概率為:

6=5/4

1.4 干擾對導引頭性能影響分析

導引頭主要由探測系統、信息處理軟件、控制與跟蹤系統組成。探測系統接收來自外部的紅外輻射,形成數字灰度圖像,信息處理軟件對圖像進行處理,實現目標的檢測、識別和跟蹤,提取目標在視場中的位置信號給控制與跟蹤系統,控制與跟蹤系統利用該信號完成導引頭的控制與跟蹤功能,并輸出視線角速度給飛控系統[13]。

導引頭抗干擾性能傳統指標為導引頭抗干擾概率[14],該指標混淆了導引頭和制導系統的功能要求,導引頭無法輸出脫靶量,無法形成抗干擾概率。導引頭的功能是完成對目標的測量,形成視線角速度。導引頭對目標的測量可以用識別概率進行評估,傳統的識別概率常選擇單幀識別概率或者片段識別概率[15],在彈道過程中,存在干擾對目標遮擋,目標不存在的場景,且該場景隨機出現,無法利用單幀識別概率進行評估,導引頭的識別是一個連續判定過程,無法合理進行片段劃分,因此,傳統的識別概率不適用于導引頭識別概率的評價。本文通過干擾對導引頭的機理影響分析,提出導引頭識別概率的評估指標。

導引頭輸出的視線角速度誤差源主要有以下3種類型:

1)跟蹤點在目標本體上一點跳到另外一點,是非干擾條件下導引頭跟蹤目標的主要誤差源;

2)跟蹤點被干擾誘偏,干擾與目標分離后識別出目標,跟蹤點回到目標,是干擾條件下導引頭跟蹤目標的典型過程;

3)跟蹤點被干擾誘偏,干擾與目標分離后,跟蹤點鎖定在干擾上,是干擾條件下導引頭跟蹤干擾的典型過程。

類型1將導致視線角速度的一個跳動,類型2將導致視線角速度的兩個跳動,第一個跳動為跟蹤目標干擾重合體,第二個跳動為跳回目標,兩個跳動的時間間隔為導引頭識別目標的時間間隔。類型3將導致視線角速度的發散。因此,除了上文分析的擾動時機和擾動幅度,導引頭識別時間對視線角速度誤差也產生影響。如圖3所示,同樣擾動幅度和擾動時機下,識別時間越長,導引頭視線角速度誤差越大。

圖3 不同識別時間的視線角速度精度歸一化曲線

根據導引頭跟蹤目標的典型過程,假設導彈系統為線性系統,則視線角速度引起的兩個跳動所導致的脫靶量可以相加,則將公式(4)更改為:

由上式可以計算不同剩余飛行時間下脫靶量滿足殺傷半徑,單顆干擾下對視線角速度的要求,如圖4所示。

圖4以制導時間常數0.5s為例,在剩余飛行時間大于5倍制導時間常數時,主要考慮導彈的過載能力和導引頭的識別能力,視線角速度誤差的要求較寬泛,隨著剩余飛行時間的縮小,視線角速度誤差的要求迅速提高,幅度與彈目接近速度等彈道條件相關,需要進行仿真確定;剩余飛行時間在時間常數的2倍至3倍區間內,允許視線角速度誤差有一定偏差,對脫靶量影響不大;剩余飛行時間在1倍時間常數附近時,視線角速度誤差要求最嚴格,此時的跟蹤點跳動將導致脫靶量大幅提高;剩余飛行時間常數小于0.5倍時間常數時,則飛控系統不響應導引頭信號,視線角速度誤差要求較為寬泛。

對于多顆干擾投放,利用線性系統相加的方式簡化計算多個擾動引起的總脫靶量,在總脫靶量滿足殺傷半徑的前提下,對視線角速度誤差在整個彈道過程中進行分配。導引頭在彈道過程中輸出的視線角速度滿足誤差要求,則認為導引頭正確識別目標。

2 多層次評估指標體系及評估過程

紅外空空導彈抗干擾評估指標體系包含了能夠反映導彈總體、制導系統及導引頭抗干擾能力的指標,如圖5所示,分別為所有要害部位的聯合殺傷概率、滿足制導誤差要求的抗干擾概率和視線角速度誤差要求的識別概率。

圖5 抗干擾評估指標體系

導彈總體的抗干擾性能評估指標為所有要害部位的聯合殺傷概率,以紅外對抗作戰樣本空間為輸入,依據目標環境模型、制導系統模型、引信模型及戰斗部模型進行聯合仿真得到,同時,依據各樣本條件下制導誤差的分布,給出滿足殺傷概率要求的各樣本制導誤差要求,為制導系統抗干擾性能評估奠定基礎。

制導系統的抗干擾性能評估指標為滿足制導誤差要求的抗干擾概率,以紅外對抗作戰樣本空間為輸入,依據目標環境模型、導引頭及飛控系統模型進行聯合仿真得到,同時,依據導彈總體對制導誤差的要求及制導時間常數的影響,給出滿足抗干擾概率要求的各樣本視線角速度誤差要求,為導引頭抗干擾性能評估奠定基礎。

導引頭的抗干擾性能評估指標為滿足視線角速度誤差要求的識別概率,以紅外對抗作戰樣本空間為輸入,依據目標環境模型、導引頭及飛控系統模型進行聯合仿真得到。

3 仿真分析

文獻[16]以國外某型導彈為例,模擬導彈相對飛機以不同的方位角和距離發射,使用仿真評估干擾彈的影響,圖6為對抗態勢示意圖,圖7為對抗結果,圖中的每一個點表示從該位置處發射的一枚導彈,干擾的投放策略是固定的,如果對抗是有效的,能夠欺騙導彈遠離飛行器,脫靶距離很大,如果導彈飛行過程中離飛行器很近,飛行器就是易受攻擊的,圖中淺色表示脫靶量小,深色表示脫靶量大。

圖6 文獻[16]干擾對抗場景

圖7 文獻[16]干擾對抗結果

從圖7中可以看出,導彈尾后攻擊較迎頭攻擊命中概率高,導彈遠距離攻擊較近距離攻擊命中概率高。

從本文對制導系統的分析可以看出,干擾擾動時的剩余飛行時間對脫靶量的影響較大,由于尾后彈道時間長,同樣干擾條件下,其剩余飛行時間較迎頭要大,干擾對視線角速度的擾動小于制導系統的要求,制導系統的脫靶量小,導彈的命中概率高。同理,相對于近距離攻擊,遠距離攻擊時,剩余飛行時間長,干擾對視線角速度的擾動小于制導系統的要求,制導系統的脫靶量小,導彈的命中概率高。

4 結論

本文在對抗環境影響機理分析的基礎上,提出紅外對抗作戰樣本空間,是抗干擾性能評估的基礎。通過干擾對紅外成像空空導彈總體、制導系統及導引頭本質能力影響機理的分析,闡述了傳統抗干擾性能評估指標的不足,并建立了包含導彈總體、制導系統和導引頭3個層次的抗干擾性能評估指標體系,分別為所有要害部位的聯合殺傷概率、滿足制導誤差要求的抗干擾概率和滿足視線角速度誤差要求的識別概率。仿真結果與本文分析結果一致,說明本文提出的指標體系具有一定的合理性。多層次的評估指標體系能夠提高對制導系統和導引頭的抗干擾性能的評估能力,提升利用導引頭、制導系統評估結果對導彈總體抗干擾性能進行預估的能力。

[1] 王煒強, 賈曉洪, 楊東升, 等. 制導武器效能評估試驗設計方法綜述與應用探討[J].航空兵器, 2015(6): 46-48.

WANG Weiqiang, JIA Xiaohong, YANG Dongsheng,et al. An Overview and Application of DOE for Guided Weapon Performance Evaluation[J]., 2015(6): 46-48.

[2] 張凱, 張邵宇, 楊東升. 紅外制導導彈抗干擾能力仿真鑒定評估技術[J]. 上海航天, 2019, 36(4): 83-89.

ZHANG Kai, ZHANG Shaoyu, YANG Dongsheng. Simulation Test and Evaluation of Anti-Jamming Capability of Infrared Guided Missile[J]., 2019, 36(4): 83-89.

[3] 王泉, 董維浩, 劉新愛, 等. 新型紅外空空導彈抗干擾能力評估分析[J]. 航天電子對抗, 2019(3): 16-19.

WANG Quan, DONG Weihao, LIU Xin’ai, et al. Analysis of evaluation for new infrared air-to-air missile countermeasures performance[J]., 2019(3): 16-19.

[4] 唐善軍, 王楓, 陳曉東. 紅外導彈抗干擾能力指標體系和評估研究[J]. 上海航天, 2017, 34(4): 144-149.

TANG Shanjun, WANG Feng, CHEN Xiaodong. Anti-Jamming Performance Index System and Evaluation of Infrared Missile[J]., 2017, 34(4): 144-149.

[5] 韓培駿, 姜長生, 許友平. 紅外成像導引系統抗干擾性能評估研究[J]. 電光與控制, 2012, 19(5): 82-86.

HAN Peijun, JIANG Changsheng, XU Youping. Anti-Jamming Performance Evaluation of Infrared Imaging Guidance Systems[J]., 2012, 19(5): 82-86.

[6] 許友平, 吳慶憲, 姜長生, 等. 基于支持向量機的紅外成像導引頭抗干擾性能評估方法[J].電光與控制, 2013, 20(12): 6-9.

XU Youping, WU Qingxian, JIANG Changsheng, et al. SVM Based Evaluation on Anti-Jamming Performance of Infrared Imaging Seekers[J]., 2013, 20(12): 6-9.

[7] 朱逸冰. 防空導彈武器系統射擊效能分析[J]. 艦船電子工程, 2019, 39(7): 154-156.

ZHU Yibing. Analysis of Firing Efficency of Antiaircraft Missile System[J]., 2019, 39(7): 154-156.

[8] 李權成, 朱傳祥, 凡永華, 等. 復雜環境對紅外空空導彈制導誤差影響規律研究[J]. 西北工業大學學報, 2019, 37(3): 457-464.

LI Quancheng, ZHU Chuanxiang, FAN Yonghua, et al. Study on Infrared Air-to-Air Missile Guidance Accuracy Affected by Complicated Environment[J]., 2019, 37(3): 457-464.

[9] 陳銘, 史志中, 蔡克榮. 彈目相對速度對防空導彈引戰配合的影響[J]. 制導與引信, 2015, 36(3): 7-11.

CHEN Ming, SHI Zhizhong, CAI Kerong. Influence of Relative Velocity of Missile and Target on Fuze Warhead Matching for Antiaircraft Missile [J]., 2015, 36(3): 7-11.

[10] 姚志敏, 史連艷, 禹海全, 等. 某型防空導彈殺傷概率評估[J]. 火力與指揮控制, 2005, 30(5): 32-37.

YAO Zhimin, SHI Lianyan, YU Haiquan, et al. The Evaluation of Anti-air Missiles' killing Probability[J]., 2005, 30(5): 32-37.

[11] 李炯, 張濤, 張金鵬. 基于零效脫靶量的制導估計一體化方法[J].彈道學報, 2017, 29(4): 35-39.

LI Jiong, ZHANG Tao, ZHANG Jinpeng. Unitization Method of Guidance and Estimation Based on Zero-Effort-Miss[J]., 2017, 29(4): 35-39.

[12] 王煒強, 賈曉洪, 付奎生, 等. 基于空中紅外對抗隨機過程的制導誤差分析[J]. 紅外技術, 2019, 41(2): 163-170.

WANG Weiqiang, JIA Xiaohong, FU Kuisheng, et al. Guidance Precision Analysis Based on Airborne IRCM Stochastic Process [J]., 2019, 41(2): 163-170.

[13] 周永恒, 崔少輝, 方丹. 紅外導引頭視線角速率誤差測試方法研究[J]. 中國測試, 2018, 44(9): 39-43.

ZHOU Yongheng, CUI Shaohui, FANG Dan. Research on accuracy test method of line of sight angle rate for infrared seeker [J]., 2018, 44(9): 39-43.

[14] 馬曉平, 趙良玉. 紅外導引頭關鍵技術國內外研究現狀綜述[J]. 航空兵器, 2018(3): 3-10.

MA Xiaoping, ZHAO Liangyu. An Overview of Infrared Seeker Key Technologies at Home and Abroad [J]., 2018(3): 3-10.

[15] 李萍, 張波, 尚怡君. 基于紅外圖像和特征融合的飛機目標識別方法[J]. 電光與控制, 2016, 23(8): 92-96.

LI Ping, ZHANG Bo, SHANG Yijun. Aircraft Target Identification Based on Infrared Image and Feature Fusion[J]., 2016, 23(8): 92-96.

[16] 滕小虎. OSSIM 系統在導彈對抗中的紅外場景仿真應用[J].艦船電子工程, 2018(8): 93-97.

TENG Xiaohu. IR Scene Simulation Application of OSSIM System in Missile Countermine[J]., 2018(8): 93-97.

Anti-Jamming Performance Evaluation Index System for Infrared Air-to-Air Missiles

ZHANG Xitao1,2,BAI Xiaodong1,2,YAN Lin1,2,WANG Weiqiang1,2

(1.,471009,;2.,471009,)

The traditional index designed to evaluate infrared anti-jamming ability is relatively simpleand is mostly determined by the comprehensive anti-jamming probability. To solve this problem, we establishan anti-jamming performance evaluation index system including three levels of missiles, guidance systems, and seekers; construct the anti-jamming performance evaluation function; decompose anti-jamming performance evaluation indexes; and improve the anti-jamming performance of the guidance system and seeker. The evaluation ability of the system improves the ability touse the evaluation results of the seeker and guidance system to predict the anti-jamming performance of the guidance system.

infrared missile, anti-jamming, index system, index refinement, performance evaluation, performance prediction

TJ761.1

A

1001-8891(2020)11-1089-06

2020-05-12;

2020-11-06.

張喜濤(1986-),男,河南洛陽,高工/碩士,研究方向是紅外制導技術。E-mail: zhangxitao86@126.com。

主站蜘蛛池模板: 无码精品国产dvd在线观看9久| 欧美亚洲欧美| 性做久久久久久久免费看| 国产成人永久免费视频| 日本国产精品一区久久久| a毛片在线| 四虎永久在线| 人妻夜夜爽天天爽| 日韩AV无码免费一二三区| 日韩精品久久久久久久电影蜜臀| 国产在线视频欧美亚综合| 青草视频在线观看国产| 亚洲成人www| 波多野结衣AV无码久久一区| 亚洲国产欧洲精品路线久久| 青青热久免费精品视频6| 国产网友愉拍精品| 亚洲日本中文字幕天堂网| 中文字幕亚洲综久久2021| 欧美人与牲动交a欧美精品| www.99精品视频在线播放| 一区二区三区四区日韩| 日韩精品亚洲精品第一页| 亚洲无码高清免费视频亚洲| 精品国产成人av免费| 国产福利一区视频| 欧美色香蕉| 国产久操视频| 激情网址在线观看| 香蕉国产精品视频| 狠狠v日韩v欧美v| 国产精品免费入口视频| 亚洲Av综合日韩精品久久久| 国产成人精品高清不卡在线| 精品成人一区二区三区电影| 动漫精品中文字幕无码| 在线无码av一区二区三区| 国产亚洲成AⅤ人片在线观看| 最新日韩AV网址在线观看| 精品国产www| 午夜视频免费一区二区在线看| 欧美亚洲国产视频| 不卡的在线视频免费观看| 欧美在线国产| 拍国产真实乱人偷精品| 人人看人人鲁狠狠高清| a级毛片视频免费观看| 国产精品嫩草影院av| 欧美日韩国产成人高清视频| 青青草原国产| 1769国产精品视频免费观看| 久久久无码人妻精品无码| 无码精品福利一区二区三区| 日本黄网在线观看| 日韩精品一区二区深田咏美| 国产高颜值露脸在线观看| 色视频国产| 在线不卡免费视频| 亚洲一道AV无码午夜福利| 日本不卡免费高清视频| 无码区日韩专区免费系列 | 狠狠色综合久久狠狠色综合| 亚洲人成日本在线观看| 精品国产女同疯狂摩擦2| 大陆精大陆国产国语精品1024 | 综合色婷婷| 一边摸一边做爽的视频17国产| 美女国产在线| 日韩精品毛片| 日韩av电影一区二区三区四区| 中文字幕在线一区二区在线| 欧美无专区| 在线免费a视频| 国产97视频在线观看| 欧美三级自拍| 国产自视频| 亚洲日韩日本中文在线| 欧美精品影院| 伊人久久婷婷五月综合97色| 婷婷色一区二区三区| 久久国产免费观看| 美女黄网十八禁免费看|