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

基于效能評估的空中對抗作戰建模分析

2015-06-19 15:39:18陳長興牛德智唐冬麗
系統工程與電子技術 2015年1期
關鍵詞:飛機模型

陳長興,牛德智,王 卓,符 輝,唐冬麗,陳 強

(1.空軍工程大學理學院,陜西西安710051;2.空軍工程大學防空反導學院,陜西西安710051)

基于效能評估的空中對抗作戰建模分析

陳長興1,牛德智1,王 卓1,符 輝1,唐冬麗2,陳 強1

(1.空軍工程大學理學院,陜西西安710051;2.空軍工程大學防空反導學院,陜西西安710051)

針對效能評估與作戰運用未能有效銜接的問題,研究了空中對抗作戰中的兵力變化問題。給出了數據鏈體制下作戰飛機效能評估計算方法,通過函數映射法將評估結果轉化為平均戰斗力水平,證明了二者的等價性和一致性。以藍徹斯特方程為基礎,建立了空中兵力變化的數學模型和有增援情況時的空戰模型。將用微分方程描述的模型離散化,指出了離散時間步長所應滿足的約束條件,形成了完整的模型求解數值計算方法。仿真實驗結果表明,數據鏈可以改變雙方作戰實力,提高作戰能力,增援作戰時選擇“局部最優增援時刻點”可以獲得更高的增援效率。

效能評估;數據鏈;數值計算;藍徹斯特方程;決策

0 引 言

武器系統的效能評估問題,不僅可以為裝備生產研制過程中提供合理的費效分析,而且可以為衡量武器的整體作戰能力提供量化的參考依據[1-4]。航空武器系統裝備數據鏈終端,是信息化聯合作戰的必然要求。近幾場局部戰爭表明,數據鏈具有獨特的戰場態勢共享、精確指揮控制和武器協同的無縫鏈接優勢,已經越來越成為武器裝備的生命線[5-6]。美軍從20世紀50年代開始研制數據鏈并先后將各種不同類型數據鏈運用到作戰飛機上,其空中作戰能力大幅度提高[7-8]。對裝備數據鏈系統的作戰飛機,進行科學合理的效能評估,并將結果運用到作戰實踐中,完成航空作戰體系循環的全過程[9-12],是一個急需完成的科學命題。

在現有關于數據鏈和作戰效能評估的研究中,一些典型成果如:文獻[13]用一種新方法研究了飛機作戰中的效能評估問題,研究結果僅限于給出不同對象的效能大小;文獻[14]討論了數據鏈支持下的空戰雙方飛機數量的變化問題,分析具有一定的合理性,但是模型中缺少一個直接反映數據鏈作用的量化指標,這樣不利于直接研究數據鏈對空戰雙方格局的影響;文獻[15]引入了戰場感知系數并將其運用到藍徹斯特方程中進行作戰分析,這對于考慮空戰兵力的變化融入了信息因素,但對模型中的作戰效能未進行相關研究。

本文出發點就是給出一種數據鏈體制下的飛機效能評估方法,并將結果運用到空戰對抗的數學建模中,實現評估與作戰的一體化。

1 單機作戰的平均戰斗力水平

1.1 基于數據鏈的效能評估方法

數據鏈下的飛機效能評估思路是:從分析數據鏈的功能,找到其對效能評估的影響因素,進而通過用解析法描述的效能評估模型給出綜合效能值。效能評估的方法較多,如典型的ADC法、AHP法等[3],本文是要選取一種能清晰描述數據鏈對飛機效能影響的方法,數據鏈對飛機效能及空戰的影響定量體現。ADC法用概率量來描述不同指標,但該法難以有效描述出數據鏈的關聯作用,如有效度中的可靠性、可信賴度中的安全性、能力中的精度等都與數據鏈相關,但是這種相關不易模型化體現。AHP法將整個飛機系統按照目標層、準則層、方案層來逐級劃分,建立層級指標體系來進行效能評估,該法雖可把數據鏈的作用考慮進去,但存在權值賦予具有一定的主觀性,指標量化也會因方法而異。所以,本文采用對數法來描述飛機效能,據此考慮數據鏈下的新效能。當然,對數法也有不足之處,此處不再做重點討論。對數法中,飛機的空戰能力指數可表示為[16- 17]

式中,I為空戰能力指數;C表示機動性;D表示火力;E表示探測目標能力;ε1~ε4分別表示生存力、電子對抗能力、航程和飛行員操縱效能系數。

由于具有態勢共享優勢,數據鏈使得超視距打擊成為可能。另外,數據鏈在傳感器、指控單元和武器平臺之間實時傳輸戰術信息,可以遂行發現即打擊的任務,實現了從傳感器到射手的無縫鏈接。分析數據鏈的這些作用及對應對數法模型中的參數可知,探測目標和電子干擾能力會在數據鏈作用下得到提高,這樣,數據鏈環境下的空戰能力可以表示為

式中,E′和ε′2分別表示使用數據鏈后的探測和電子對抗能力,且E′=E+λ1E,ε′2=ε2+λ2ε2,λ1和λ2表示數據鏈對這兩項參數的提升率。

對于λ1的確定,需要根據不同作戰模式和數據鏈體制有所不同,這里僅考慮一種簡單的雙機協同作戰時數據鏈對探測能力提升率的計算方式。假定雙機為同向直線飛行的情況,則提升率可以用探測面積的變化比來表示:

式中,φ為雷達搜索的方位角;r為雷達探測距離;d為雙機之間距離(d<r)。

對于λ2的確定,可以在現有處理方法基礎上對其統一賦值研究。文獻[16]給出了機載電子電子對抗設備的電子對抗能力,全向雷達告警系統為1.05,全向雷達告警+消極干擾投放+導彈逼近警告系統的綜合系統為1.16~1.20,有數據鏈支撐下的電子對抗設備,還會增加積極干擾投放功能,所以對其賦值時可以使其適當大于無數據鏈支持下的電子對抗能力。

1.2 空戰能力指數轉化為平均戰斗力水平

由式(1)及文獻研究的結果知,空戰能力指數一般遠大于1。在空戰中,考慮單機一次戰斗只打擊一架敵機的情況,則應該用單機的平均戰斗力水平來衡量殺傷效率或作戰威力。

作如下變換:

將I變換到在區間(0,1)上的平均戰斗力水平γ(其中k>0,作用是避免γ在數值上過于集中)。為此證明γ和I是等價的。

證明 令A={x|x∈Z,x>0},B={x|x∈Z,0<x<1},作一映射f:A→B,有

因為x∈A,則ln(kx)∈(-∞,+∞)。且arctan(ln(kx))∈(-π/2,π/2),所以1/π·arctan(ln(kx))∈(-1/2,1/2)。即有1/π·arctan(ln(kx))+1/2∈(0,1)。

說明f∶A→B是由A到B的一一映射,這樣A和B是等勢(等價)的,即A~B。

因為x>0,k>0,所以f′(x)>0。且f(x)為增函數,即在式(4)中,γ隨I的增大而增大,隨γ的減小而減小,二者變化關系一致。證畢

2 空戰中的兵力變化建模

2.1 藍徹斯特方程分析方法

在對作戰兵力變化規律的研究中,本文以藍徹斯特方程為基礎,并分析不同作戰場景下的動態兵力變化過程。此處給出藍徹斯特方程(平方律)的表達式為

式中,R,B為紅、藍方數量;t為時間;α,β為藍、紅方每個戰斗單位的平均戰斗力水平。其求解可通過式(1)~式(4)來具體確定。

對式(7)求解微分方程組的狀態解為

令x=1-R/R0表示紅方兵力的消耗率,y=B/B0表示藍方兵力的剩余率。L=β·R20/α·B20表示紅藍雙方的戰斗實力比,L>1時表示紅方戰斗實力強,L<1時表示藍方戰斗實力強,L=1時雙方實力相等。將x,y,L代入式(8),得

式(9)描述了在雙方實力比L約束下的消耗率與剩余率之間關系,在L確定的情況下,可以分析出戰斗結局。

2.2 數據鏈體制下的空中力量對抗建模

考慮在數據鏈情況下的空戰模型。首先求解平均戰斗力水平的提升率ηR和ηB,即有ηR=(γR′-γR)/γR,ηB=(γ′B-γB)/γB。

在空戰過程中,綜合性能更優的數據鏈一方擁有更大優勢,不僅可以實施更加積極主動的電子干擾,而且可以在有目標來襲時在更遠距離上發出預警提示等,這樣,定義一數據鏈水平占優因子為

則數據鏈體制下的空戰兵力變化規律建模為

從式(10)中可知,η實質是個可正可負的量,為正時對紅方有增進作用,反之對藍方。

類似式(7)~式(9)的分析可得

式中,(1+η)/(1-η))·L定義為等效實力比。

另外,考慮有增援時的模型為

式中,μ1(t),μ2(t)分別為雙方的增援率(可以恒定,也可以變化);t1,t2為增援開始時刻;ε(t)為階躍函數,用來限制增援兵力加入時間。

2.3 模型的數值計算方法

用微分方程組描述的數學模型,意義比較明確,體系也比較完整,但是微分方程組的求解一般比較困難。若在式(13)的基礎上在綜合考慮別的因素(如非戰斗減員等),空戰對抗模型將更加復雜,雖然在某些情況下可以借助特定軟件求解,但是模型中各變量的變化規律及內在的影響關系不易分析和掌握。所以,需要將微分方程組描述的數學模型離散化,研究其數值計算方法。

首先討論最基本的藍徹斯特方程,由微分方程的定義入手,可以將式(7)中的第1個方程(第2個方程的分析方法類似)離散化為

式中,Δt為一個一定小的時間步長,在該時間步長的Δt下,可以基本保持原微分方程的變量變化規律。在離散取值時,依次令t=0,Δt,2Δt,…,nΔt(考慮在一場空戰中,它是由若干個戰斗波次構成的,而每個戰斗波次又可以細分為一定數量的戰斗回合,則此處的n可以理解為總的戰斗回合數),且將RnΔt簡記為Rn,則式(14)可化為

在數值計算時,需要確定Δt取多大的時間步長才合適,文獻[15]中將藍徹斯特方程直接轉化為差分方程形式,其實質是直接令Δt=1,這樣處理不僅會失去建立在微分方程基礎上的作戰模型中變量的變化規律(主要是與Δt足夠小的要求相違背),而且所得的計算數據太過片面性,對研究作戰規律不利。為此,寫出式(15)的完整形式為

將式(16)中方程組的左右兩邊分別相加,得到

當藍方最后失敗時,有Bn-1為0(或近似為0),且可以將Bn-1的變化近似用等差數列來描述,則式(17)化為

進一步,可以求出

由式(8)可知,令B=0,可以得到

至此,得到了離散取樣間隔與戰斗回合數(或描述戰斗兵力變化的細分過程數)之間的關系,在數值計算(包括式(11)和式(13))時,可以取時間間隔小于或者等于式(18)中的時間間隔,從而研究戰斗中的雙方兵力變化過程。類似地,若采用式(7)中的第2個方程進行討論,考慮最終失敗的為紅方,則時間間隔可取為

這樣,數據鏈體制下的空戰兵力變化離散模型可以表述為

式中,m1,m2=0,1,2,…為正整數,表示某一增援時刻。在戰斗結束時(n=N),紅方(假定為獲勝方)共增援的飛機數可以表示為Rs=μ1·Δt·(N-nm1+1),共損失的飛機數可以表示為Rd=R0+Rs-RN。

3 仿真算例

假定在空中對抗作戰中,紅、藍方所用分別為米格-29和F-15E機型。仿真中假定僅紅方有數據鏈支持。米格-29的雷達型號為NO-9,其最大探測距離為90 km,這樣協同作戰距離為70 km時,根據式(3)可求λ1=0.95。設定紅、藍方的電子對抗能力分別為1.05和1.2,紅方使用數據鏈后子對抗能力變為1.8。評估指標中的其余參數參照文獻[17],將飛機效能的每項參數列表顯示并求解I(k取1/10)和γ如表1所示。

表1 飛機作戰效能指標及效能評估值

從表中可以看出,數據鏈的引入對空戰能力指數和平均戰斗力水平都有所提升。

3.1 戰斗實力比對空戰結局的影響

取L=[0.6,0.8,1,1.2,1.4],仿真藍方剩余率與紅方消耗率的變化曲線如圖1所示。

圖1 藍方剩余率和紅方消耗率的變化關系

圖1 的結果驗證了第2.1節中對L不同空戰結局不同的分析。

由表1中數據可知ηR=0.21,則η=0.21。這樣在原來的L下仿真可得到剩余率與消耗率曲線(見圖2)。分析圖像知,在數據鏈情況下,當紅、藍方的實力比介于0.6~0.8之間時,雙方實力平衡。令(1+η)/(1-η))·L=1,可求L=0.65。這說明,采用數據鏈后,紅方只要大于藍方實力的65%時就可獲勝,反之失敗。圖2中也標出了雙方實力平衡的分界線。

圖2 引入數據鏈后剩余率和消耗率的變化關系

3.2 不同的增援時刻點對空戰的影響

設定R=10,B=15,由表1結果知,α=0.7和β=0.6, η=0.21。由實力比(或等效實力比)計算知,紅方整體實力仍然弱于藍方,這樣實驗中僅假定紅方有增援(首先考慮nm1=0的情況)且增援率為2.5。通過式(15)~式(22),仿真4種模式(分別為原始實力下雙方作戰、紅方僅引入數據鏈雙方作戰、紅方僅擁有一定速率增援力量雙方作戰、紅方同時擁有數據鏈和增援力量雙方作戰)下雙方的動態損耗見圖3。

圖3 不同作戰條件下的雙方動態損耗圖

從圖3可以看出,由于紅方實力較弱,原始條件下的作戰結果為失敗;在無數據鏈支持下,一定的增援速率下并不能改變最終戰局(當增援率一定大時可以改變)。在數據鏈體制下,紅方有效延長了戰斗進程,這在增援力量到來的情況下最終扭轉戰局,獲得勝利。

接下來確定紅方實施增援的有效增援時間范圍。令增援開始時刻m1分別為0、15和30,取離散時間步長為Δt=0.025,仿真得到雙方兵力變化的動態損耗曲線見圖4。

圖4 紅方不同的增援起始時刻時雙方的動態損耗圖

由圖4知,當m1=30時,增援不能改變紅方失敗的結果。這說明,空戰增援中存在一個時間分界點,在此之前增援,可以扭轉雙方戰局;反之不能。經過多次仿真可以確定有效增援時刻范圍為[0,28],換算到時間范圍即為[0,0.7]。

進一步,讓m1在[0,28]之間取值,步進量為1,其余仿真參數不變,得到的雙方動態損耗曲線見圖5。

圖5 有效增援時刻范圍內m1變化時的雙方動態損耗

當m1增大時(即增援開始時刻越晚),雙方的空戰進程維持時間變長,紅方勝利時剩余的飛機數量越少,由于相鄰時間上的飛機數量變化規律比較接近,為進一步分析增援效率的高低,圖6給出了m1不同時紅方共增援和損失的飛機數量變化曲線。

圖6 m1值不同時的紅方增援和損失飛機數量變化

分析圖6可知,越在增援的后期,需增援的飛機總數越多,同時損失的飛機數量也急劇攀升,這說明增援力量越早加入戰斗,越能以較少的飛機數量贏得戰斗。由此產生的另外一個問題是,若增援只能在更小的一個時間范圍內實施(如圖6中第10~15個時刻點之間進行增援),則需要給出一個科學有效的增援時刻點選取準則。

為此,引入“損增比”的概念,表示整個戰斗中單方損失的飛機數量與增援的飛機數量之比,其深層含義是增援單位數量的飛機時損失的飛機數量,可給出其定義式為

根據式(24)計算“損增比”ρ,ρ值比較小的時刻點可以作為增援時刻的選擇點。圖7給出了與圖6對應的“損增比”變化曲線。

圖7 不同增援時刻時的損增比曲線

從圖7中可以看出,隨m1值變化時,“損增比”曲線并不是簡單的單調遞增關系,而是存在一系列極小值點,在這些極小值點周圍,損增比數值都相對較大,說明增援獲得的效率并不是最高(在有效增援時間范圍的末期,雖然損增比數值較小,但是其需要的增援飛機數量較多,即需要較大的增援規模才能實現,實際作戰中并不一定總能滿足這種規模要求)。據此可得出結論,在不同的增援時刻周圍,存在使增援效率最高的局部最優點,即“局部最優增援時刻點”。掌握了該規律,便于指揮員對整場空戰中在何時進行空中力量增援,以獲得最高增援效率提供參考依據。

4 結 論

本文對數據鏈體系下飛機空中對抗作戰的兵力變化問題進行了研究。對數據鏈情形下的飛機效能評估模型進行改進,并將評估結果轉化為平均戰斗力水平運用到空戰模型中。以藍徹斯特方程為基礎進行分析空中對抗過程中的兵力變化規律,建立了具有反映數據鏈效應的空中力量體系對抗數學模型和有增援時的空戰模型。為更好地分析和把握空戰規律,將微分方程模型離散化,給出了其數值計算方法,為仿真實現空中力量的變化過程提供了一種有效的實現途徑。仿真算例說明了數據鏈對戰斗能力的提升作用,增援時應選擇“局部最優增援時刻點”以獲得最大的增援效率,實驗結果為運籌和指揮空戰全局提供了重要參考。

[1]Yang F,Yang C C,Liang L,et al.New approach to determine common weights in DEA efficiency evaluation model[J].Journal of Systems Engineering and Electronics,2010,21(4):609- 615.

[2]Liu H,Xu M,Chen N Y,et al.Task effectiveness evaluation on anti-ship attack of fighter aircraft[J].Transactions of Nanjing University of Aeronautics&Astronautics,2012,29(2):164- 170.

[3]Yang Y L,Tai H X,Shi T.Weighting indicators of building energy efficiency assessment taking account of experts'priority[J].Journal of Central South University of Technology,2012,19(3):803- 808.

[4]Gray F B,Hommedieu W.A combat-effectiveness approach to information assurance[EB/OL].[2013- 10- 30].http:∥enu.kz/repository/2010/AIAA- 2010- 1748.pdf.

[5]Xu Y L,Li Q Z,Xie S F,et al.Study on algorithm and communication protocol of differential GPS positioning based on pseudo range[J].Information Technology and Application,2009,1(3):606- 609.

[6]Tang L,Zhang R R,Chen Q B.Cross-layer design of combining AMC with H ARQ in cooperative relay system with perfect and imperfect CSI[J].Journal of Harbin Institute of Technology(New Series),2012,19(4):118- 128.

[7]Enel L,Guillem F.Improvements in navy data networks and tactical communication systems[J].World Academy of Science,Engineering and Technology,2008,2(6):97- 101.

[8]Zhou L,Zhang H Y,Wang T,et al.Static check of WS-CDL documents[C]∥Proc.of the IEEE International Symposium on Service Oriented System Engineering,2008:142- 147.

[9]Lee J,Kang S H,Rosenberger J,et al.A hybrid approach of goal programming for weapon systems selection[J].Computers &Industrial Engineering,2010,58(3):521- 527.

[10]Chen J,Mu Y P,Li D F.Study on effect evaluation of weapon coordinated use of fleet[J].Military Operations Research and Systems Engineering,2009,23(3):66- 71.

[11]Johnson I R,Mackay N J.Lanchester models and the battle of Britain[J].Naval Research Logistics,2011,58(3):210- 222.

[12]Yan J J,Wang Y,Wang X Y.Modeling and simulation of Lanchester equation based on game theory[C]∥P roc.of International Conference on the Business Management and Electronic Information,2011:212- 214.

[13]Xiang G D,Fu X D,Jiao L M.A method for efficiency evaluation under complex circumstance[J].Electronics Optics&Control,2013,20(1):19- 23.(相國棟,傅曉冬,焦利明.復雜環境下效能評估方法[J].電光與控制,2013,20(1):19- 23.)

[14]Xia B H,Zhou Y Y.Influence analysis of data link on modern air-air combat effectiveness[J].Computer Engineering,2011,37(9):257- 260.(夏白樺,周于義.數據鏈對現代空-空作戰效能的影響分析[J].計算機工程,2011,37(9):257- 260.)

[15]Wu J,Yang F,Liang Y,et al.Generalized lanchester combat model for information warfare[J].Fire Control&Command Control,2010,35(2):50- 53.(吳俊,楊峰,梁彥,等.面向信息化戰爭的廣義蘭切斯特作戰模型[J].火力與指揮控制,2010,35(2):50- 53.)

[16]Zhu B L,Zhu R C,Xiong X F.Aircraft effectiveness evaluation[M].2nd ed.Beijing:Aviation Industry Press,2006.(朱寶鎏,朱榮昌,熊笑非.作戰飛機效能評估[M].2版.北京:航空工業出版社,2006.)

[17]Huang J C,Zhang Y,Yang L,et al.A quantitative evaluation method for operation efficiency of tactical data link[J].Fire Control&Command Control,2010,35(12):66- 69.(黃金才,張勇,楊磊,等.戰術數據鏈作戰效能定量評估方法[J].火力與指揮控制,2010,35(12):66- 69.)

Modeling of air combat based on effectiveness evaluation

CHEN Chang-xing1,NIU De-zhi1,WANG Zhuo1,FU Hui1,TANG Dong-li2,CHEN Qiang1
(1.College of Science,Air Force Engineering University,Xi’an 710051,China;
2.College of Air and Missile Defense,Air Force Engineering University,Xi’an 710051,China)

Aiming at aircraft effectiveness evaluation’s connection with combat application,force change problem in air combat is studied.Effectiveness evaluation computation method under data link mechanism is given,and value of effectiveness evaluation is converted to average combat level by function mapping means,which is improved to equivalent and consistent with the value.On basis of lanchester equation,mathematic model of force change in air combat is set up,including the case of being reinforced.The model described by differential equation is dispersed with a subject in discrete time internal,thus the whole numerical calculation method about solving the model is formulated.Simulation shows that data link could change combat force of each other and improve combat ability,and greater efficiency could be obtained by one’s choosing local optimum reinforcement time to join the combat.

effectiveness evaluation;data link;numerical computation;lanchester equation;decision making

TN 957.51

A

10.3969/j.issn.1001-506X.2015.01.14

陳長興(1964-),男,教授,博士,主要研究方向為現代通信理論、信息系統建模、效能評估。

E-mail:xachenchangxing@yahoo.com

牛德智(1984-),通信作者,男,博士研究生,主要研究方向為通信系統與雷達信息處理、系統工程與效能評估。

E-mail:niudezhi_001@163.com

王 卓(1984-),男,講師,碩士,主要研究方向為信號與信息處理、作戰理論與應用。

E-mail:wz93861@163.com

符 輝(1986-),男,碩士,主要研究方向為信息處理技術。

E-mail:43349291@qq.com

唐冬麗(1987-),男,碩士,主要研究方向為測控技術。

E-mail:tangdongli123@126.com

陳 強(1989-),男,碩士研究生,主要研究方向為信息系統建模與仿真。

E-mail:chenqiang2013@163.com

1001-506X(2015)01-0079-06

網址:www.sys-ele.com

2013- 11- 04;

2014- 03- 29;網絡優先出版日期:2014- 05- 30。

網絡優先出版地址:http:∥w ww.cnki.net/kcms/detail/11.2422.TN.20140530.1441.006.html

陜西省電子信息綜合集成重點實驗室項目(201107Y16)資助課題

猜你喜歡
飛機模型
一半模型
鷹醬想要“小飛機”
飛機失蹤
環球時報(2022-05-30)2022-05-30 15:16:57
國航引進第二架ARJ21飛機
重要模型『一線三等角』
重尾非線性自回歸模型自加權M-估計的漸近分布
“拼座飛機”迎風飛揚
當代陜西(2019年11期)2019-06-24 03:40:28
乘坐飛機
3D打印中的模型分割與打包
神奇飛機變變變
主站蜘蛛池模板: 免费Aⅴ片在线观看蜜芽Tⅴ| 欧美一级大片在线观看| 99精品在线看| 欧洲免费精品视频在线| 制服丝袜亚洲| 狠狠综合久久| 好紧好深好大乳无码中文字幕| 自拍偷拍一区| 欧美一区二区人人喊爽| 国产迷奸在线看| 欧美一级黄片一区2区| 国产亚洲欧美日韩在线一区二区三区| 丝袜久久剧情精品国产| 影音先锋丝袜制服| 久久久久久高潮白浆| 国国产a国产片免费麻豆| 无码啪啪精品天堂浪潮av| 91精品国产福利| 欧美成人精品欧美一级乱黄| 天堂av综合网| 久久久精品国产亚洲AV日韩| 亚洲精品另类| 99精品在线视频观看| 综1合AV在线播放| 一级毛片基地| 91精品网站| 激情综合五月网| 久久国产毛片| 中文字幕无码电影| 91色在线观看| 中文无码毛片又爽又刺激| 99国产在线视频| 亚洲成综合人影院在院播放| 99精品久久精品| 久久久久久尹人网香蕉| 干中文字幕| 婷婷五月在线| 亚洲福利片无码最新在线播放| 91偷拍一区| 在线欧美国产| 国产日本一区二区三区| 亚洲天堂视频网站| 毛片网站免费在线观看| 香蕉在线视频网站| 2021最新国产精品网站| 麻豆精品在线| 91九色国产porny| 一级毛片在线免费看| 亚洲一区二区日韩欧美gif| 亚洲第一视频网| 亚洲国产精品VA在线看黑人| 亚洲av无码专区久久蜜芽| 91九色最新地址| 国产你懂得| 日韩色图区| 国产精品流白浆在线观看| 久久精品无码一区二区国产区| 风韵丰满熟妇啪啪区老熟熟女| 午夜视频www| 天天综合天天综合| 二级毛片免费观看全程| 成人日韩欧美| 国产经典在线观看一区| 不卡网亚洲无码| 久久女人网| 亚洲AV人人澡人人双人| 亚洲女同欧美在线| 五月婷婷综合网| 国产91高跟丝袜| 日韩国产精品无码一区二区三区| 午夜三级在线| 国产男女免费视频| 国产一级二级在线观看| 国产菊爆视频在线观看| 国产不卡网| 色哟哟色院91精品网站| 亚洲人成网站18禁动漫无码| 三级欧美在线| 日韩大乳视频中文字幕| 一级毛片中文字幕| 国产丰满大乳无码免费播放 | 99久久精品国产综合婷婷|