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

采用遺傳算法對(duì)Jaumann吸波結(jié)構(gòu)材料優(yōu)化設(shè)計(jì)*

2017-06-27 08:14:35鄭日升張慶兵劉恒軍肖志河戚開(kāi)南
現(xiàn)代防御技術(shù) 2017年3期
關(guān)鍵詞:優(yōu)化結(jié)構(gòu)

鄭日升,張慶兵,劉恒軍,肖志河,戚開(kāi)南

(1.北京電子工程總體研究所,北京 100854;2.電磁散射重點(diǎn)實(shí)驗(yàn)室,北京 100854)

采用遺傳算法對(duì)Jaumann吸波結(jié)構(gòu)材料優(yōu)化設(shè)計(jì)*

鄭日升1,張慶兵1,劉恒軍1,肖志河2,戚開(kāi)南2

(1.北京電子工程總體研究所,北京 100854;2.電磁散射重點(diǎn)實(shí)驗(yàn)室,北京 100854)

以頻率帶寬為優(yōu)化目標(biāo)函數(shù),采用遺傳算法針對(duì)3層Jaumann吸波結(jié)構(gòu)材料進(jìn)行優(yōu)化設(shè)計(jì)。分析了不同介電常數(shù),不同隔離層層數(shù),不同隔離層厚度以及不同電阻參數(shù)對(duì)吸波材料反射率的影響;在吸波結(jié)構(gòu)材料反射率T<-15 dB條件下,獲得了較寬的頻率帶寬,并得到了該頻率帶寬內(nèi)相應(yīng)介電常數(shù),隔離層厚度和電阻的優(yōu)化參數(shù)。

材料隱身;雷達(dá)吸波結(jié)構(gòu);遺傳算法;反射率;阻抗匹配;頻率帶寬

0 引言

降低飛行器雷達(dá)散射截面通常有2種方法:外形隱身和材料隱身。由于飛行器外形隱身具有較大的局限性,比如影響飛行器的氣動(dòng)特性。因此通常結(jié)合材料隱身技術(shù)提高飛行器的生存能力和突防能力。而吸波結(jié)構(gòu)材料隱身技術(shù)是降低飛行器反射率的一種重要方法[1],也是目前研究隱身材料技術(shù)的熱門(mén)課題之一。

迄今為止,隱身材料已經(jīng)進(jìn)行了接近60年的研究,被廣泛應(yīng)用于各國(guó)的武器裝備和尖端設(shè)施上。美國(guó)是最早將材料技術(shù)應(yīng)用于武器裝備,并得到了較好的效果,其中最早應(yīng)用飛行器是F-117,使得該飛行器雷達(dá)散射截面平均值0.1~0.01 m2,該隱身飛機(jī)的尾翼和機(jī)身等部件采用了復(fù)合材料,但是給飛行器增加了額外10%的質(zhì)量,并且該飛行器結(jié)合了氣動(dòng)外形隱身技術(shù),對(duì)發(fā)動(dòng)機(jī)艙內(nèi)采用了蜂窩夾芯結(jié)構(gòu)外形,降低了在某一頻段內(nèi)的雷達(dá)散射截面[2-3]。F-22是美國(guó)一種最為典型的隱身飛行器,其中在進(jìn)氣道、舵面前緣等部件都安裝了雷達(dá)吸波結(jié)構(gòu)復(fù)合材料。由于受到空氣動(dòng)力性能等方面的限制,雷達(dá)罩采用了特殊外形設(shè)計(jì)與布局,2種技術(shù)的結(jié)合,使得該飛行器雷達(dá)散射截面強(qiáng)散射源得到了有效的控制[4-6]。材料隱身技術(shù)的研究主要是針對(duì)材料本身,即通過(guò)研究不同吸收劑的電磁特性,以及選擇吸收劑、粘結(jié)劑和基本材料的搭配比例,研制出具有高吸收性的新型材料,然后將其整體涂敷在飛行器表面[7-9]。吸波材料隱身技術(shù)越來(lái)越受到國(guó)內(nèi)外專家們的重視。目前,吸波材料隱身技術(shù)也在不斷的改進(jìn),主要針對(duì)材料反射率低、吸收頻帶寬、質(zhì)量輕、隔離層厚度小為優(yōu)化目標(biāo)。國(guó)內(nèi)外都在研制一種新型的蜂窩夾心結(jié)構(gòu)吸波材料[10-12]:由玻璃鋼、蜂窩夾芯和碳纖維材料3部分組成。目前,采用優(yōu)化算法對(duì)降低吸波結(jié)構(gòu)材料的反射率也得到了較好的效果[13-15],為材料隱身實(shí)驗(yàn)測(cè)試提供基礎(chǔ)。

1 吸波結(jié)構(gòu)材料反射率分析計(jì)算方法

1.1 反射率計(jì)算方法

吸波結(jié)構(gòu)材料的平板結(jié)構(gòu)通常采用電阻抗加載多層介質(zhì),通過(guò)多層介質(zhì)之間的阻抗匹配組合,使得在某一反射率條件下達(dá)到較寬頻帶,并且達(dá)到最佳的電磁波吸收能力。吸波結(jié)構(gòu)材料一般包括3層,外層是透波層,通常是吸波材料制成,中間是電磁吸收層,一般遵循阻抗?jié)u變?cè)瓌t,內(nèi)層是反射層,通常是用金屬基底或呈金屬反射特性的碳纖維復(fù)合材料構(gòu)成。此層具有反射電磁波的特性,吸波結(jié)構(gòu)材料是基于相位相互抵消的原理進(jìn)行研制的。

首先將遺傳算法應(yīng)用到多層吸波材料優(yōu)化設(shè)計(jì),將所要求頻帶內(nèi)整體反射率最小作為優(yōu)化目標(biāo),本文將目標(biāo)函數(shù)定為反射率T<-15 dB條件下的頻率帶寬,在指定吸波材料厚度d/mm和電阻R/Ω以及介電常數(shù)e約束的情況下,準(zhǔn)確計(jì)算出最優(yōu)的多層吸波材料組合。由于多層吸波結(jié)構(gòu)材料目標(biāo)函數(shù)是求解一個(gè)極值問(wèn)題,因此可以用遺傳算法來(lái)優(yōu)化。如圖1所示,是典型的Jaumann 吸波結(jié)構(gòu),金屬板在最底層,在底板上層是有一定厚度和阻抗的材料,為計(jì)算方便,本文只針對(duì)吸波結(jié)構(gòu)型材料其中3個(gè)參數(shù)進(jìn)行優(yōu)化研究。

圖1 多層吸波結(jié)構(gòu)材料示意圖Fig.1 Sketch of multi-layer absorbing structure material

其中,各層吸波材料的輸入阻抗為

(1)

式中:R為電阻;Zin1為輸入阻抗;Zl12為負(fù)載阻抗,即任何一個(gè)電氣元件都有一定的阻抗,相對(duì)于電源而言,該阻抗就是負(fù)載阻抗。

當(dāng)電磁波從初始阻抗的自由空間入射到輸入阻抗為Zin1的界面上時(shí),一部分經(jīng)飛行器表面直接反射,另一部分電磁波進(jìn)入吸波材料中,多層吸波材料對(duì)電磁波的反射特性可由反射率T計(jì)算公式為

(2)

P為反射系數(shù),轉(zhuǎn)化為反射率為

T=201g|P|.

(3)

1.2 目標(biāo)函數(shù)構(gòu)造

在實(shí)際工程需求方面考慮,吸波結(jié)構(gòu)材料通常由3層以上的材料組成,本文初步探討3層吸波材料的反射率進(jìn)行優(yōu)化。反射率與入射頻率、隔離層厚度和相對(duì)介電常數(shù)等參數(shù)密切相關(guān),要想得到較好的隱身效果,即得到不同參數(shù)的一個(gè)最佳排列組合。所以反射率就是關(guān)于相對(duì)介電常數(shù)、材料電阻和各層厚度的函數(shù),記為

(4)

式中:T為材料的不同組合方式得到的反射率。

本文中目標(biāo)函數(shù)構(gòu)造考慮反射率T<-15 dB的面積S最小為目標(biāo),如圖2所示,陰影部分為頻率區(qū)間[fa,fb]上的反射率曲線,當(dāng)反射率T<-15 dB時(shí),就不考慮,也就是說(shuō)當(dāng)T>-15 dB才增加不同頻率下的反射率,當(dāng)多層吸波材料在指定頻率區(qū)間內(nèi)達(dá)到最佳吸波效果時(shí),本文取反射率T的總和為極小值,即反射率曲線與頻率軸圍成的區(qū)域的面積應(yīng)是最小值,如圖2所示。

圖2 吸波材料的反射率面積示意圖Fig.2 Sketch of reflectivity area of absorbing structure material

(5)

則求面積S的極小值內(nèi)的頻率帶寬為目標(biāo)函數(shù),即確保T<-15 dB條件下,在頻帶寬內(nèi)的反射率面積最小,從而可以在確保整個(gè)頻帶內(nèi)得到反射率的最優(yōu)解,而不只是一個(gè)頻率點(diǎn)上得到最好的結(jié)果,因此這是求解全局最優(yōu)解問(wèn)題,更加符合實(shí)際工程中需要在一定的頻帶內(nèi)達(dá)到隱身效果的要求,采用遺傳算法能夠針對(duì)該目標(biāo)函數(shù)進(jìn)行優(yōu)化。

2 遺傳算法優(yōu)化設(shè)計(jì)

2.1 編碼方案設(shè)計(jì)

基于以上分析,以頻率帶寬為目標(biāo)函數(shù),首先針對(duì)優(yōu)化編碼進(jìn)行設(shè)計(jì),采用二進(jìn)制進(jìn)行染色體編碼。以Ti表示第i層的材料反射率編碼,則第i層材料的編碼是形如MiTi一個(gè)進(jìn)制串。而對(duì)于層吸波材料整體,其染色體編碼形式為

M1T1M2T2…MiTi.

利用遺傳算法進(jìn)行優(yōu)化設(shè)計(jì)時(shí),所采用的編碼方案為:用3位二進(jìn)制碼串表示吸波材料相對(duì)介電常數(shù),用7位二進(jìn)制碼串表示材料的隔離層厚度,則i層吸波材料的二進(jìn)制串長(zhǎng)度應(yīng)為10i。

2.2 適應(yīng)度函數(shù)設(shè)計(jì)

設(shè)種群中某個(gè)個(gè)體的反射率為T(mén),由式(3)可以計(jì)算得出Ti,Ti就是所要求的種群中某個(gè)個(gè)體的反射率,并且需要求解目標(biāo)函數(shù)的最小值。由于反射率尺取對(duì)數(shù)后總是負(fù)值,假設(shè)個(gè)體的適應(yīng)度函數(shù)為E,令Fi=-Ti,因此,反射率目標(biāo)函數(shù)面積最小值等價(jià)于適應(yīng)度函數(shù)的最小值,對(duì)于遺傳算法來(lái)說(shuō),這是一個(gè)求解全局最小值的問(wèn)題。

2.3 選擇、交叉、變異的設(shè)計(jì)

遺傳優(yōu)化算法采用比例選擇方法,采用單點(diǎn)交叉算子作為交叉算子,選用基本位變異算子作為變異算子?;谶z傳算法的多層吸波材料優(yōu)化初始參數(shù)設(shè)置為:種群大小M為600,運(yùn)行代數(shù)迭代范圍為0~600,交叉概率G為0.1~0.99,變異概率G1為0.001~0.1。

2.4 基于遺傳優(yōu)化算法吸波結(jié)構(gòu)材料頻帶優(yōu)化設(shè)計(jì)

下面首先分析3個(gè)典型不同參數(shù)對(duì)反射率的影響,研究3層吸波結(jié)構(gòu)材料在頻率2~18 GHz范圍內(nèi)反射率與頻率關(guān)系。圖3是未經(jīng)優(yōu)化得到的反射率T與頻率之間的關(guān)系,從圖中可以看出,在反射率T<-15 dB時(shí),在頻率5~7 GHz范圍內(nèi)反射率T較小,其次是頻寬14~16 GHz范圍內(nèi)的反射率T略小于-15 dB。而在其余的頻寬內(nèi)反射率T都大于-15 dB。即陰影部分的面積。本文目的就是求陰影部分面積最小條件下的頻寬。

圖3 頻率帶寬2~18 GHz條件下優(yōu)化前反射率T<-15 dBFig.3 Reflectivity T<-15 dB under the conditionof frequency bandwidth at 2 GHz~18 GHzbefore optimization

下面初步分析不同介電常數(shù)e,不同厚度d及不同隔離層層數(shù)N對(duì)反射率的影響,如圖4所示,比較了e=1,e=3,e=5 3種介電常數(shù)的反射率。當(dāng)e=5時(shí),反射率相對(duì)較小,約在3.7 GHz

圖4 介電常數(shù)e對(duì)吸波材料反射率T的影響Fig.4 Effect of dielectric constant e on the reflectivity T of absorbing material

下面分析不同電阻R對(duì)反射率T的影響。其中R分別為:500,700,900 Ω。從圖5中可以發(fā)現(xiàn),隨著電阻R增大,反射率T減小,當(dāng)R=900 Ω時(shí),反射率T約-22 dB,R=700 Ω時(shí)。反射率T達(dá)到約-16 dB,而當(dāng)R=500 Ω時(shí),反射率T增大為-7 dB。當(dāng)R低于700 Ω時(shí),反射率T大于-15 dB,如圖5所示。

圖5 阻抗R對(duì)吸波材料反射率T的影響Fig.5 Effect of impedance R on reflectivity T of absorbing material

圖6是分析隔離層厚度d對(duì)反射率T的影響。從圖中發(fā)現(xiàn),不同的隔離層厚度對(duì)反射率T也有一定的影響。當(dāng)d=5 mm時(shí),在頻率范圍3.2 GHz

圖6 隔離層厚度對(duì)反射率的影響Fig.6 Effect of the thickness of the isolation layer on the reflectivity

從圖7中發(fā)現(xiàn),當(dāng)N=1時(shí),反射率T都大于-15 dB,最小反射率達(dá)到約-7 dB。當(dāng)N=2時(shí),反射率T有一定程度的減小,最小反射率約-11 dB,而當(dāng)N=3時(shí),有2個(gè)頻帶范圍內(nèi)T小于-15 dB。通過(guò)不同隔離層層數(shù)分析發(fā)現(xiàn),當(dāng)層數(shù)增加時(shí),反射率T逐漸減小,頻帶也隨著增加。但是具體層數(shù)需要根據(jù)實(shí)際需要多方面衡量。

圖7 隔離層層數(shù)N對(duì)反射率的影響Fig.7 Effect of the number of isolation layer N on the reflectivity

下面是針對(duì)3層吸波結(jié)構(gòu)型材料9個(gè)參數(shù)變量進(jìn)行優(yōu)化設(shè)計(jì),圖8是遺傳算法優(yōu)化計(jì)算界面。

通過(guò)遺傳優(yōu)化算法得到在頻率2~18 GHz范圍內(nèi),在T<-15 dB條件下較寬的頻帶(5.5 GHz

圖8 遺傳算法運(yùn)行迭代界面示意圖Fig.8 Sketch of genetic algorithm running iterative interface

圖9 反射率T<-15 dB條件下遺傳算法優(yōu)化頻帶結(jié)果Fig.9 Optimized bandwidth results under the condition of reflectivity at T<-15 dB based on genetic algorithm

表1~3是表示優(yōu)化后吸波結(jié)構(gòu)材料組合得到相應(yīng)的參數(shù)值,并獲得了在該反射率面積內(nèi)反射率總和為-233.82 dB。

表1 隔離層厚度dTable 1 Thickness of isolation layer d mm

表2 金屬片電阻RTable 2 Metal resistance R Ω

表3 介電常數(shù)e及反射率T總和Table 3 Dielectric constant e and the sum of the reflectivity T

3 結(jié)束語(yǔ)

本文通過(guò)構(gòu)造了用遺傳算法設(shè)計(jì)寬帶吸收結(jié)構(gòu)材料的合理目標(biāo)函數(shù),獲得了頻率范圍2~18 GHz條件下反射率T<-15 dB時(shí)較寬的頻帶,并得到了相應(yīng)的相對(duì)介電常數(shù),隔離層厚度以及金屬片電阻的優(yōu)化值。

首先根據(jù)典型的3層吸波結(jié)構(gòu)材料反射率的分布,構(gòu)造了反射率T<-15 dB時(shí),以反射率曲線與頻率軸圍成的區(qū)域面積最小值為目標(biāo),確定目標(biāo)函數(shù)為頻帶寬。

其次分析了不同介電常數(shù)和不同金屬片電阻對(duì)反射率的影響。發(fā)現(xiàn)隨著介電常數(shù)的增加,吸波結(jié)構(gòu)材料的反射率T增加,頻率帶寬也相應(yīng)減小。獲得了隨著電阻R增大,反射率T反而減小,而且頻率帶寬也較大。

接著研究了不同隔離層厚度以及隔離層層數(shù)對(duì)反射率的影響,研究表明,當(dāng)隔離層厚度增加時(shí),反射率減小,頻率帶寬增加;隔離層層數(shù)增加時(shí),反射率T逐漸減小,頻帶相應(yīng)的增加。

最后通過(guò)遺傳算法針對(duì)反射率T<-15 dB進(jìn)行頻帶寬優(yōu)化設(shè)計(jì),得到了在3層吸波結(jié)構(gòu)材料在頻率5.5 GHz

[1] 孫敏.隱身材料技術(shù)[M].北京:國(guó)防大學(xué)出版社,2013. SUN Min.Stealth Material Technology[M].Beijing:The National Defense University Press,2013.

[2] 趙東林,周萬(wàn)城.結(jié)構(gòu)吸波材料及其結(jié)構(gòu)型式設(shè)計(jì)[J].兵器材料科學(xué)與工程,1997,20(6):53-57. ZHAO Dong-lin,ZHOU Wan-cheng.Structural Absorbing Materials and Structure Design[J].Weapon of Materials Science and Engineering,1997,20(6):53-57.

[3] 李江海,孫秦.結(jié)構(gòu)型吸波材料隱身特性的優(yōu)化設(shè)計(jì)算法研究[J].兵器材料科學(xué)與工程,2004,27(2):9-11. LI Jiang-hai,SUN Qin.Research on Characteristic of Absorbing Material Stealth by Structural Design Optimization Algorithm[J].Weapon of Materials Science and Engineering,2004,27(2):9-11.

[4] FARMER S J,LAIGHT A,GRANT I D.Structural Absorbent Materials[C]∥IEE Colloquium on Low Profile Absorbers and Scatter(Digest No.125),1992.London:UK,May 1992,1-4

[5] J.瓊斯.隱身技術(shù)——黑色魔力的藝術(shù)[M].北京:航空工業(yè)出版社,1991:33-70. JONES J.Stealth Technology——The Art of Black Magic[M].Beijing:Aviation Industry Press,1991:33-70.

[6] 趙云峰.先進(jìn)復(fù)合材料在結(jié)構(gòu)隱身材料中的應(yīng)用[J].宇航材料工藝,1989(Z1):4-5. ZHAO Yun-feng.The Application of Advanced Composite Materials in the Structure Stealth Materials[J].Aerospace Material Process,1989(Z1):4-5.

[7] 夏新仁.隱身技術(shù)發(fā)展現(xiàn)狀與趨勢(shì)[J].中國(guó)航天,2002(1):40-44. XIA Xin-ren.Stealth Technology Development Status and Trends[J].China’s Space,2002(1):40~44.

[8] 同武勤,凌永順,蔣金水.雷達(dá)材料隱身及效果測(cè)試技術(shù)分析[J].艦船電子對(duì)抗,2005,28(1):25-28. TONG Wu-qin,LING Yong-shun,JIANG Jin-shui.Radar Stealth Material and Effect Test Technical Analysis[J].Shipboard Electronic Countermeasure,2005,28(1):25-28.

[9] SMITH F C.Effective Permittivity of Dielectric Honeycombs[J].IEE Proceedings Microwaves Antennas and Propagation,1999,146(1):55-59.

[10] 王略,章仲安.低RCS飛行器外形設(shè)計(jì)實(shí)踐[J].航空學(xué)報(bào),1995,16(5):692-695. WANG Lue,ZHANG Zhong-an.The Design of Low RCS Aircraft Shape[J].Journal of Aviation,1995,16(5):692-695.

[11] BUCCI O M.Scattering Form Wedge Tapered Absorbers[J].IEEE Tram,1994,19(7):96-104.

[12] MICHAEL A.Pegasus to Test Naval UACV Concept[J].Aviation Week & Space Technology,2001(8):20-2.[13] MICHIELSSEN E,SAJER J M,RANJITHAN S,et al.Design of Light Weight,Broadband Microwave Absorbers Using Genetic Algorithm[J].IEEE Transactionson Microwave Theory and Techniques,1993,41(7):1024-1031.

[14] MOSALLAEI H,RAHMAT S Y.RCS Reduction of Canonical Targets Using Genetic Algorithm Synthesized RAM [J].IEEE Transactions on Antennas and Propagation,2000,48(10):1594-1606.

[15] 張秀成,何華輝.遺傳算法在涂層吸波材料設(shè)計(jì)中的應(yīng)用[J].磁性材料及器件,1999,30(5):11-14. ZHANG Xiu-cheng,HE Hua-hui.Application of Genetic Algorithm in the Coating Absorbing Material Design[J].Magnetic Materials and Devices,1999,30(5):11-14.

Optimized Design of Jaumann Radar Absorbing Materials Using Genetic Algorithm

ZHENG Ri-sheng1,ZHANG Qing-bing1,LIU Heng-jun1,XIAO Zhi-he2,QI Kai-nan2

(1.Beijing Institute of Electronic System Engineering,Beijing 100854,China;2.Electromagnetic Scattring Laboratory,Beijing 100854,China)

An optimized design of Jaumann radar absorbing materials using a genetic algorithm is proposed. The critical parameters are optimally designed for the frequency bandwidth. The phenomenon analysis of radar absorbing materials is based on different relative permittivity, the thickness of isolation layer and the resistance. The wide frequency bandwidth, the optimized parameter value of relative permittivity, the isolation layer and the resistance under the reflectivity (T<-15 dB) are obtained.

material stealth;radar absorbing structure;genetic algorithm;reflectivity;impedance matching;frequency bandwidth

2016-09-30;

2016-11-30

鄭日升(1981-),男,湖北崇陽(yáng)人。工程師,博士,研究方向?yàn)轱w行器設(shè)計(jì)。

通信地址:100854 北京市142信箱30分箱 E-mail:zhengrisheng2000@163.com

10.3969/j.issn.1009-086x.2017.03.005

TP301.6;TJ765.5

A

1009-086X(2017)-03-0028-06

猜你喜歡
優(yōu)化結(jié)構(gòu)
超限高層建筑結(jié)構(gòu)設(shè)計(jì)與優(yōu)化思考
《形而上學(xué)》△卷的結(jié)構(gòu)和位置
民用建筑防煙排煙設(shè)計(jì)優(yōu)化探討
關(guān)于優(yōu)化消防安全告知承諾的一些思考
一道優(yōu)化題的幾何解法
由“形”啟“數(shù)”優(yōu)化運(yùn)算——以2021年解析幾何高考題為例
論結(jié)構(gòu)
新型平衡塊結(jié)構(gòu)的應(yīng)用
模具制造(2019年3期)2019-06-06 02:10:54
論《日出》的結(jié)構(gòu)
創(chuàng)新治理結(jié)構(gòu)促進(jìn)中小企業(yè)持續(xù)成長(zhǎng)
主站蜘蛛池模板: 久久综合婷婷| 思思99思思久久最新精品| 日本免费a视频| 免费人成黄页在线观看国产| 国产黑人在线| 欧洲成人在线观看| 亚洲第一成网站| 免费观看国产小粉嫩喷水| 婷婷色一二三区波多野衣| 99久久精品免费看国产电影| 成人福利在线观看| 好吊妞欧美视频免费| 99在线观看免费视频| 国产91在线|中文| 无码人妻免费| 日本在线欧美在线| 久久国产精品影院| 波多野结衣一区二区三区AV| 国产91丝袜在线播放动漫 | 中文字幕无码制服中字| 99re在线视频观看| 亚洲高清在线天堂精品| 福利姬国产精品一区在线| 欧美笫一页| 欧美日韩中文国产| 欧美日韩福利| 91精品国产自产在线老师啪l| 女人天堂av免费| 99久久精品免费观看国产| 三上悠亚在线精品二区| 亚洲欧美日韩另类| 国产在线精彩视频二区| 亚洲最大福利网站| 国产91色| 国产色网站| AV老司机AV天堂| 中文一区二区视频| 欧美a在线看| 高清无码不卡视频| 欧美视频免费一区二区三区| 99热这里只有免费国产精品| 日韩黄色大片免费看| 91精品久久久久久无码人妻| 无码免费视频| 久久精品丝袜高跟鞋| 666精品国产精品亚洲| 欧美色香蕉| 不卡网亚洲无码| 亚洲看片网| 91成人在线免费视频| 日本亚洲成高清一区二区三区| a亚洲视频| 亚洲无码精品在线播放| 欧美成人怡春院在线激情| 中国丰满人妻无码束缚啪啪| 成人免费网站在线观看| 国产美女免费网站| 国产在线八区| 性欧美久久| 国产丝袜啪啪| 亚洲国产精品成人久久综合影院| 国产在线无码一区二区三区| 久久国产精品波多野结衣| 亚洲自偷自拍另类小说| 日韩欧美中文亚洲高清在线| 成人综合在线观看| 99伊人精品| 99精品视频播放| 日本91视频| 在线精品视频成人网| 亚洲日韩欧美在线观看| 日韩国产 在线| 欧美a级完整在线观看| aa级毛片毛片免费观看久| 亚洲第一极品精品无码| 992Tv视频国产精品| 性色一区| 国产一区成人| 婷婷六月综合网| 国产午夜无码专区喷水| 国产一区二区人大臿蕉香蕉| 欧美日韩在线第一页|