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

一種毫米波UWB多模態(tài)OAM介質(zhì)諧振器陣列天線

2023-03-07 01:26:12王艷妮孫學(xué)宏劉麗萍
無線電工程 2023年2期
關(guān)鍵詞:模態(tài)

王艷妮,孫學(xué)宏,劉麗萍,2,方 亮

(1.寧夏大學(xué) 物理與電子電氣工程學(xué)院,寧夏 銀川 750021;2.寧夏沙漠信息智能感知重點(diǎn)實(shí)驗(yàn)室,寧夏 銀川 750021)

0 引言

隨著移動(dòng)互聯(lián)網(wǎng)應(yīng)用的蓬勃發(fā)展,頻譜資源越來越匱乏,傳統(tǒng)的提高頻譜利用率的技術(shù),如正交頻分復(fù)用、時(shí)分復(fù)用和碼分復(fù)用等,已經(jīng)無法滿足高速?gòu)?fù)雜的信號(hào)傳輸。軌道角動(dòng)量(Orbital Angular Momentum,OAM)技術(shù)作為一種新型的無線通信技術(shù),能夠有效增大系統(tǒng)容量、提高頻譜利用率[1-2],攜帶OAM的渦旋電磁波可包含無限多正交模態(tài),并且不同模態(tài)之間都相互正交且互不干擾,能夠?qū)崿F(xiàn)不同模態(tài)的拓展OAM信號(hào)共用同一載頻進(jìn)行傳輸[3],為現(xiàn)有的無線通信復(fù)用技術(shù)提供了一種新的思路。

目前,國(guó)內(nèi)外學(xué)者已經(jīng)應(yīng)用了多種方法來產(chǎn)生OAM無線電波,例如螺旋反射陣天線[4-6]、貼片天線[7-8]和介質(zhì)諧振器天線[9-12]。其中,介質(zhì)諧振器天線由于其輻射單元無金屬損耗,天線輻射效率較高,受到了廣泛關(guān)注。2017年,文獻(xiàn)[9]設(shè)計(jì)了一款工作在17.38 GHz的介質(zhì)諧振器OAM天線,該天線通過激發(fā)回音壁模式產(chǎn)生OAM波束[9]。文獻(xiàn)[10]利用縫隙耦合圓柱形介質(zhì)諧振器設(shè)計(jì)了一種攜帶渦旋電磁波的OAM天線,該天線通過每個(gè)天線單元的延遲輻射相互作用產(chǎn)生2種模態(tài)的OAM波束[10]。2020年,文獻(xiàn)[11]設(shè)計(jì)了一種半球形介質(zhì)諧振器天線,可以在5.8和7.4 GHz處分別產(chǎn)生l為2和3的OAM波束[11];同年,文獻(xiàn)[12]以水作為介質(zhì)設(shè)計(jì)了一種OAM天線,該天線能夠在不同諧振頻率處分別生成2種不同模態(tài)的OAM波束[12],但是目前設(shè)計(jì)的OAM天線大多分布在C波段、X波段和Ku波段,仍然存在頻帶擁擠的問題。工信部將30~300 GHz用于第五代移動(dòng)通信(5G)毫米波(Millimeter Wave,MMW)頻段的技術(shù)研發(fā)試驗(yàn)[13],MMW通信不僅可以實(shí)現(xiàn)大容量信息傳輸[14],還可以有效減小天線尺寸,因此在高速發(fā)展的信息時(shí)代,開發(fā)MMW通信技術(shù)將有效緩解頻譜資源緊張的問題[15-16]。本文設(shè)計(jì)了一種MMW多模態(tài)OAM介質(zhì)諧振器陣列天線,該天線中心頻率在35 GHz,仿真結(jié)果表明該天線能夠在MMW頻段穩(wěn)定產(chǎn)生多種模態(tài)的OAM波束,相對(duì)帶寬達(dá)到74.3%,將其應(yīng)用在OAM-MIMO系統(tǒng)當(dāng)中能夠有效增加系統(tǒng)容量,在MMW無線通信領(lǐng)域?qū)⒂袕V闊的應(yīng)用前景。

1 相關(guān)原理

1.1 環(huán)形相控陣原理

環(huán)形相控陣列天線模型如圖1所示。

圖1 環(huán)形相控陣列天線模型Fig.1 Ring phased array antenna model

將矩形介質(zhì)諧振器按圖1所示放置,利用圓形相控陣列天線原理[17],將N個(gè)輻射單元等間距地設(shè)置在半徑為d的圓環(huán)上,第n個(gè)輻射單元上的電流分布為:

In=Ie-jlφne-βφn,

(1)

式中,I為電流初始強(qiáng)度大小。當(dāng)激勵(lì)第n個(gè)陣元時(shí),它的相位角應(yīng)為φn=n·2πl(wèi)/N,其中φn表示第n個(gè)陣元的空間相位角;l為OAM的模態(tài)值;β為電流強(qiáng)度沿圓環(huán)傳遞時(shí)的衰減系數(shù),在理想情況下,β=0,環(huán)上的電流分布為In=Ie-jlφn。

每個(gè)單個(gè)介質(zhì)諧振器為一個(gè)天線原件,其輻射矢量位是:

(2)

(3)

通過Maxwell方程變形可以得到含有相位因子e-jlφn的電場(chǎng)和磁場(chǎng)分量,該相位因子就是產(chǎn)生旋波所需要的相位因子[18]。

1.2 介質(zhì)諧振器相關(guān)理論

介質(zhì)諧振器天線是指能夠限定電磁能量在一定區(qū)域內(nèi)振蕩的結(jié)構(gòu)。當(dāng)能量在諧振腔中周期性振蕩時(shí),金屬側(cè)壁可以視為理想電壁,在能量振蕩過程中,電磁波發(fā)生完全反射,能量不會(huì)輻射出腔體;在介質(zhì)諧振器中,側(cè)壁只能近似于磁壁,能量在振蕩過程中部分會(huì)輻射出去。電磁波在介質(zhì)界面上的反射與透射如圖2所示,一束平面電磁波Et斜入射向介質(zhì)諧振器側(cè)壁,在介質(zhì)諧振器與空氣交界面處發(fā)生透射和反射。

圖2 電磁波在介質(zhì)界面上的反射與透射Fig.2 Reflection and transmission of electromagnetic waves at dielectric interfaces

按照折射定律:

(4)

式中,εr為介質(zhì)諧振器的相對(duì)介電常數(shù)。

根據(jù)式(4),存在臨界角θ0,當(dāng)θi>θ0時(shí),發(fā)生全反射,不發(fā)生透射。由折射定律,當(dāng)εr趨于無窮大時(shí),θi趨近于0,即電磁波在諧振器內(nèi)振蕩時(shí),任意入射角度的波均可以看作全反射,不向外輻射能量,此時(shí)諧振器側(cè)壁可以看作理想磁壁。但實(shí)際應(yīng)用中介質(zhì)材料的相對(duì)介電常數(shù)有確定值,故在電磁波振蕩過程中,部分能量會(huì)透射出去,從而使介質(zhì)諧振器形成輻射體。

本文采用的是矩形介質(zhì)諧振器,其結(jié)構(gòu)如圖3所示,擁有x,y,z,εr四個(gè)維度的可調(diào)參數(shù),設(shè)計(jì)自由度較高,可以通過選擇適當(dāng)?shù)拈L(zhǎng)寬比來克服模式退化,從而限制了非期望模式的激發(fā),并可以降低交叉極化水平。基于Bladel的理論,Mongia評(píng)估了由Okaya和Barash首先提出的關(guān)于矩形介質(zhì)諧振器中TE和TM模式的分類方法。經(jīng)過理論推導(dǎo)發(fā)現(xiàn),矩形介質(zhì)諧振器中的TM模式不滿足矩形介質(zhì)諧振器中的電磁場(chǎng)分布條件,同時(shí)在實(shí)驗(yàn)中也沒有發(fā)現(xiàn)TM模式的存在。然而,對(duì)于TE模式來說,無論是基于理論推演還是實(shí)驗(yàn)測(cè)量都論證了其真實(shí)存在,因此,目前關(guān)于矩形介質(zhì)諧振器的模式研究也主要集中在對(duì)TE模式的開發(fā)。根據(jù)極化方向分為TEx,TEy,TEz模式,其基模為TEδ11模。介質(zhì)諧振器工作時(shí)直接置于導(dǎo)體面上,故其偶次模被短路,僅存在奇次模。由于其內(nèi)部電磁場(chǎng)分布十分復(fù)雜,選用Mongia提出的矩形介質(zhì)諧振器的介質(zhì)波導(dǎo)模型進(jìn)行求解[19]。

圖3 矩形介質(zhì)諧振器結(jié)構(gòu)示意Fig.3 Rectangular dielectric resonator structure

(5)

根據(jù)相應(yīng)的邊界條件,有:

(6)

(7)

(8)

(9)

(10)

式中,α,γ為x,y方向的衰減常量;k0為自由空間的波數(shù)。理想模式下,場(chǎng)僅分布在諧振器內(nèi)部,近似得kx=mπ/a,ky=nπ/b。結(jié)合超越方程[20]:

(11)

得到介質(zhì)諧振器中各個(gè)方向電場(chǎng)和磁場(chǎng)值為:

Ex=0,

(12)

(13)

(14)

(15)

(16)

(17)

矩形介質(zhì)諧振器基模TEδ11的諧振頻率為[21]:

(18)

矩形介質(zhì)諧振器天線的輻射Q值可以表示為:

(19)

式中,We為介質(zhì)諧振器天線儲(chǔ)存的電能;Pm為磁偶極子的等效磁偶極矩。

(20)

(21)

1.3 天線陣列耦合分析

互耦是陣列天線固有的特征,各單元之間的耦合作用與單元之間的距離成反比[22]。由于互耦的存在,陣列中各天線單元的輻射特性和阻抗特性與自由空間中的天線單元有所不同,會(huì)導(dǎo)致陣列整體性能的改變[23]。每一個(gè)天線陣元都可以等效為存在電磁耦合的開放型電路[24],因而,在發(fā)射信號(hào)時(shí),單個(gè)天線陣元的表面電流不僅有饋電電流,還會(huì)有相鄰陣元散射激勵(lì)引起的電流。陣列天線陣元間的互耦示意如圖4所示,在接收信號(hào)時(shí),a0號(hào)陣元接收的信號(hào)是入射波以及周圍天線陣元響應(yīng)的疊加。簡(jiǎn)言之,各陣元之間存在由能量耦合引起的互耦效應(yīng)的存在會(huì)導(dǎo)致天線與理想情況出現(xiàn)偏差,而這種偏差往往就是增益下降以及旁瓣電平值增大等,最終導(dǎo)致陣列天線方向性變差。

圖4 陣列天線陣元間的互耦示意Fig.4 Schematic of mutual coupling between array antenna elements

基于上述模型,可以考慮將陣列天線用端口網(wǎng)絡(luò)等效表示,陣元天線由N個(gè)單元組成,且天線單元均可等效為阻抗ZL,其內(nèi)部阻抗用Zg表示,激勵(lì)源的開路電壓表示為Vg,陣列天線等效端口網(wǎng)絡(luò)如圖5所示。

圖5 陣列天線等效端口網(wǎng)絡(luò)Fig.5 Equivalent port network of array antenna

用Zij(i≠j)表示陣元i和陣元j之間的互阻抗,Zii表示陣元i的自阻抗,vj表示端口j的電壓,voj表示j的開路電壓。那么,以上的網(wǎng)絡(luò)關(guān)系可以表示為[25]:

(22)

式(22)用矩陣的形式可以表示為:Z0V=V0,其中,V0為陣列天線端口的開路電壓,Z0是用ZL歸一化處理后的陣列天線廣義阻抗矩陣,矩陣V代表陣列天線的輸出電壓組成的矩陣。

在陣列天線中,互耦的大小通常還可以用互耦系數(shù)Sij更直觀地度量:

(23)

S=(Z+Z0)-1(Z-Z0),

(24)

式中,Z0為對(duì)角矩陣,一般對(duì)角元素為50 Ω;Z為天線陣的阻抗矩陣。將天線陣列等效為端口網(wǎng)絡(luò)后,以上2個(gè)量就包含了天線間的互耦信息,從而可以將互耦的影響量化分析。

2 天線結(jié)構(gòu)設(shè)計(jì)

介質(zhì)諧振器具有多饋電方式、高輻射率和低損耗率的特點(diǎn),更適用于超寬帶(Ultra-wide Band,UWB)天線的設(shè)計(jì),本文將矩形介質(zhì)諧振器作為OAM陣列天線的輻射陣元,為了減小饋電網(wǎng)絡(luò)對(duì)輻射單元的影響,降低陣元間耦合,所以采用同軸探針對(duì)天線進(jìn)行饋電。陣元的三維結(jié)構(gòu)示意如圖6(a)所示,其中W是介質(zhì)諧振器的寬度,L是介質(zhì)諧振器的長(zhǎng)度,H是諧振器的高度,諧振器采用介電常數(shù)為6的TP材料。陣列天線的俯視圖如圖6(b)所示,陣列天線的仰視圖如圖6(c)所示。天線采用介電常數(shù)為2.2的RT5880作為介質(zhì)基板,H1是探針的高度,R是介質(zhì)基板的半徑,H2是介質(zhì)基板的厚度,xf是探針距離原點(diǎn)的距離,d是圓環(huán)半徑(各個(gè)陣元與原點(diǎn)的距離)。天線結(jié)構(gòu)的各個(gè)參數(shù)如表1所示。

圖6 天線結(jié)構(gòu)示意Fig.6 Antenna structure schematic

表1 天線結(jié)構(gòu)參數(shù)Tab.1 Antenna structure parameters 單位:mm

用相控陣方法產(chǎn)生OAM波束,需要N個(gè)相同的陣元均勻放置在同一半徑的圓上,給每個(gè)陣元饋以等幅、等相位差的信號(hào)。第n個(gè)陣列單元激勵(lì)時(shí)的相位角用φn=n·2πl(wèi)/N表示,其中,OAM波束的模態(tài)l取值與陣元的個(gè)數(shù)N有關(guān),需滿足-N/2

本文設(shè)計(jì)的天線陣元數(shù)為12,即N=12,通過迭代優(yōu)化,將陣元分布在半徑為11 mm的圓環(huán)上,陣元沿陣列中心順時(shí)針旋轉(zhuǎn)30°組成一個(gè)陣列。根據(jù)-N/2

3 仿真結(jié)果與分析

3.1 參數(shù)分析

基于圖6的模型,通過仿真分析,得到如圖7所示的回波損耗S11隨圓環(huán)半徑d=0.8λ,d=0.9λ,d=λ,d=1.1λ的變化曲線,其中λ是波長(zhǎng),可以看出S11在圓環(huán)半徑為λ時(shí)最小,并且隨著圓環(huán)半徑的減小,頻率明顯上移,說明圓環(huán)半徑d對(duì)于天線的匹配情況影響極大,當(dāng)d=λ時(shí),S11小于-30 dB,帶寬范圍達(dá)到最大,匹配良好。

圖7 S11隨圓環(huán)半徑d的變化曲線Fig.7 Variation curve of S11 with the radius of the circle d

當(dāng)d=λ時(shí),本文進(jìn)一步研究了探針高度對(duì)S11的影響。S11隨探針高度H1變化(H1為0.6,0.7,0.8,0.9,1.0,1.1,1.2)的變化曲線如圖8所示,可以看出H1的增加能夠引起S11減小,當(dāng)H1到達(dá)1 mm時(shí),S11最小,再增加H1的高度反而會(huì)使S11增大。

圖8 S11隨探針高度H1的變化曲線Fig.8 Variation curve of S11 with probe height H1

基于以上仿真結(jié)果分析,最終確定設(shè)計(jì)的天線圓環(huán)半徑為d=λ,H1=1.0 mm。此時(shí),中心頻率為35 GHz時(shí),工作頻率在30~56 GHz,S11小于-10 dB,匹配良好,相對(duì)帶寬達(dá)到74.3%,遠(yuǎn)遠(yuǎn)超過UWB所定義的25%,進(jìn)一步提高了該天線抗多徑衰落、數(shù)據(jù)傳輸安全等優(yōu)勢(shì)。

3.2 仿真結(jié)果分析

在40 mm×40 mm×20 mm處觀察到的不同模態(tài)下的相位分布如圖9所示。圖像呈現(xiàn)螺旋分布狀態(tài),符合OAM渦旋電磁波的特征。當(dāng)OAM模式值為正時(shí),得到的是逆時(shí)針旋轉(zhuǎn)的渦旋相位波前;當(dāng)OAM模式值為負(fù)時(shí),得到的是順時(shí)針旋轉(zhuǎn)的渦旋相位波前。由仿真結(jié)果可以看出,本文設(shè)計(jì)的天線產(chǎn)生的OAM渦旋電磁波具有良好的旋轉(zhuǎn)性,即說明本文提出的介質(zhì)諧振器天線可以很好地產(chǎn)生5種獨(dú)立的OAM波束,可分別獨(dú)立用于多路復(fù)用、信道調(diào)制和信息傳輸。另外,若將陣元的饋電相位差由正變負(fù),又可產(chǎn)生l為-1,-2,-3,-4,-5的渦旋波束,此時(shí)的相位分布如圖9(f)~(j)所示。當(dāng)l=±5時(shí),天線的遠(yuǎn)場(chǎng)輻射圖略有失真,這是因?yàn)樘炀€所能生成的模態(tài)數(shù)受陣元個(gè)數(shù)的嚴(yán)格控制-N/2

圖9 相位分布Fig.9 Phase distribution

圖10是l為±1,±2,±3,±4,±5的E面(Phi=0°),H面方向圖(Phi=90°)。可以看到,方向圖覆蓋面形狀較好,且因?yàn)榫哂辛己玫膶?duì)稱性,所以E面和H面的增益方向圖幾乎相同。在0°和180°處,陣列增益急劇下降,電磁波輻射能量很低,并且隨模態(tài)l增加,呈現(xiàn)主瓣逐漸減少,旁瓣不斷增大的趨勢(shì),這是OAM波束特有的現(xiàn)象。

(a)E面

(b)H面圖10 天線方向圖Fig.10 Radiation pattern of the proposed antenna

因此,OAM在應(yīng)用于無線通信時(shí),會(huì)出現(xiàn)主軸方向上的輻射能量為零的現(xiàn)象。

不同模態(tài)下的3D增益圖如圖11所示。

出現(xiàn)中心上方下凹的現(xiàn)象,符合渦旋電磁波中心能量最低、邊緣能量高的物理特性。與傳統(tǒng)的OAM天線相比,在本文中,當(dāng)l=±1時(shí)OAM的狀態(tài)有所不同,這種差異是由于陣元差異造成的,使用的矩形介質(zhì)諧振器本身輻射能力強(qiáng),互偶性強(qiáng),因此造成中心能量的指向性增強(qiáng)。此外,由圖11可以看到,隨著模態(tài)數(shù)的增加,增益不斷降低,中心空域面積增大,圖10也顯示了該特征,所以應(yīng)該正確選擇接收天線的半徑,并且將其垂直分布在陣列軸的波束周長(zhǎng)上,以確保接收位置不會(huì)處在OAM的空域范圍。接收天線示意如圖12所示,左邊是陣列天線,右邊的紅點(diǎn)表示2個(gè)相鄰的接收天線,通過二者之間的相位差確定OAM狀態(tài)。

圖12 接收天線示意Fig.12 Receiving antenna schematic

近幾年的參考文獻(xiàn)與本文的對(duì)比如表2所示。相比較而言,本文設(shè)計(jì)不僅結(jié)構(gòu)緊湊、工作頻段高,而且實(shí)現(xiàn)了UWB。此外,本文應(yīng)用多陣元陣列天線,實(shí)現(xiàn)了OAM多模態(tài)的靈活切換,有利于OAM在模態(tài)復(fù)用方面的進(jìn)一步研究。

表2 參考文獻(xiàn)與本文對(duì)比Tab.2 Comparison of references with this paper

4 OAM-MIMO信道容量驗(yàn)證

設(shè)計(jì)了一個(gè)毫米波多模態(tài)OAM通信系統(tǒng),如圖13所示。根據(jù)發(fā)射天線和接收天線的互易性,采用相同的介質(zhì)諧振器陣列天線作為發(fā)射天線和接收天線,利用Matlab仿真平臺(tái)進(jìn)行OAM調(diào)制解調(diào),從而完成OAM信道搭建,利用計(jì)算的信道數(shù)據(jù)分析了基于OAM的MIMO通信與普通MIMO通信對(duì)信道容量的影響情況。

圖13 OAM通信系統(tǒng)Fig.13 OAM communication system

利用有耗環(huán)形理論建立了介質(zhì)諧振器的等效模型,設(shè)定N=12,則信道矩陣可表示為[26]:

(25)

(26)

式中,a∈N×1;A表示激勵(lì)信號(hào)幅度矩陣;A∈L×L是數(shù)模矩陣A=diag{A0,A1,…,AL},L=10,表示模態(tài)總數(shù)量;X表示激勵(lì)信號(hào)矩陣;G表示調(diào)制矩陣:

(27)

接收端的信號(hào)為:

b=Ha+n,

(28)

式中,n表示零均值方差為σ2的加性高斯白噪聲矩陣。OAM信道矩陣可表示為:HOAM=MHG,M是解調(diào)矩陣,M=GH,OAM系統(tǒng)的信道容量可表示為:

(29)

圖14 OAM-MIMO與傳統(tǒng)MIMO系統(tǒng)信道容量對(duì)比Fig.14 Channel capacity comparison of OAM-MIMO and conventional MIMO system

5 結(jié)束語

本文基于環(huán)形相控陣原理提出了一種工作在MMW頻段的UWB介質(zhì)諧振器天線,用以產(chǎn)生OAM波束。討論了圓環(huán)半徑和探針高度對(duì)回波損耗以及天線帶寬的影響,天線的工作頻段為30~56 GHz,相對(duì)帶寬達(dá)到74.3%,具有工作頻率高、電尺寸小的優(yōu)點(diǎn),因此對(duì)于加工工藝具有一定的挑戰(zhàn)性。仿真結(jié)果表明,該天線在等幅、等相位差的條件下可以穩(wěn)定產(chǎn)生l為±1,±2,±3,±4,±5的OAM波束,并且每種波束都具有良好的旋轉(zhuǎn)性和對(duì)稱性。此外,將該天線應(yīng)用于OAM-MIMO系統(tǒng)將有效提高系統(tǒng)容量,使OAM在MMW頻段的模態(tài)復(fù)用具有一定的現(xiàn)實(shí)意義。后期將對(duì)本文設(shè)計(jì)的天線進(jìn)行加工試驗(yàn),從而在未來進(jìn)行推廣。

猜你喜歡
模態(tài)
基于BERT-VGG16的多模態(tài)情感分析模型
跨模態(tài)通信理論及關(guān)鍵技術(shù)初探
一種新的基于模態(tài)信息的梁結(jié)構(gòu)損傷識(shí)別方法
多跨彈性支撐Timoshenko梁的模態(tài)分析
車輛CAE分析中自由模態(tài)和約束模態(tài)的應(yīng)用與對(duì)比
國(guó)內(nèi)多模態(tài)教學(xué)研究回顧與展望
基于HHT和Prony算法的電力系統(tǒng)低頻振蕩模態(tài)識(shí)別
由單個(gè)模態(tài)構(gòu)造對(duì)稱簡(jiǎn)支梁的抗彎剛度
利用源強(qiáng)聲輻射模態(tài)識(shí)別噪聲源
日版《午夜兇鈴》多模態(tài)隱喻的認(rèn)知研究
電影新作(2014年1期)2014-02-27 09:07:36
主站蜘蛛池模板: 久久综合色天堂av| 99er这里只有精品| 亚洲日本韩在线观看| 亚洲免费三区| 热99re99首页精品亚洲五月天| 国产sm重味一区二区三区| 亚洲男人天堂2020| 小说 亚洲 无码 精品| 亚洲精品色AV无码看| 97se亚洲综合不卡| 国产00高中生在线播放| 88国产经典欧美一区二区三区| 无码日韩人妻精品久久蜜桃| 亚洲色无码专线精品观看| 精品国产91爱| 日韩中文精品亚洲第三区| 亚洲第一黄色网址| 国产情精品嫩草影院88av| 99热这里只有免费国产精品| 久久久精品久久久久三级| 欧美午夜性视频| 91精品免费久久久| 欧美色伊人| 亚洲美女操| 国产综合无码一区二区色蜜蜜| 国产不卡在线看| 久久这里只有精品国产99| 婷婷成人综合| 强乱中文字幕在线播放不卡| 99久久亚洲精品影院| 在线免费亚洲无码视频| 啪啪永久免费av| 美女黄网十八禁免费看| 亚洲bt欧美bt精品| 国产丝袜第一页| 国产精品.com| 亚洲国产欧美国产综合久久| 国产精品久久久久久影院| 亚洲精选无码久久久| 97视频在线观看免费视频| 精品国产免费第一区二区三区日韩| 色婷婷视频在线| 亚洲综合色婷婷| 欧美影院久久| 亚洲欧州色色免费AV| 成人伊人色一区二区三区| 久久久久久久久18禁秘| 亚洲日韩精品无码专区| 亚洲第一成年网| 久久五月视频| 国产91全国探花系列在线播放| 久久国语对白| 色九九视频| 香蕉eeww99国产精选播放| 被公侵犯人妻少妇一区二区三区| 超碰91免费人妻| 高清欧美性猛交XXXX黑人猛交| 在线看片中文字幕| 第九色区aⅴ天堂久久香| 人妻精品久久无码区| 国产综合网站| 国产区人妖精品人妖精品视频| 99ri国产在线| 91黄视频在线观看| 一本综合久久| 色综合中文综合网| 婷五月综合| 人妻丰满熟妇αv无码| 欧洲欧美人成免费全部视频| 国产成人欧美| 丝袜国产一区| 丁香六月综合网| 97人人做人人爽香蕉精品| 欧美另类视频一区二区三区| 国产导航在线| 亚洲伊人久久精品影院| 欧美区日韩区| 99无码熟妇丰满人妻啪啪| AV老司机AV天堂| 四虎亚洲精品| 91久久精品日日躁夜夜躁欧美| AV老司机AV天堂|