亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        雙信道應(yīng)答器多天線間的互耦分析及降耦研究

        2017-04-09 10:15:56朱林富趙會兵鐘志旺陳建譯
        中國鐵道科學(xué) 2017年3期

        朱林富,趙會兵,鐘志旺,陳建譯

        (1.北京交通大學(xué) 電子信息工程學(xué)院,北京 100044; 2.北京交通大學(xué) 軌道交通控制與安全國家重點(diǎn)實驗室,北京 100044;3.廣州鐵路(集團(tuán))公司 電務(wù)處,廣東 廣州 510088)

        應(yīng)答器是射頻識別技術(shù)(Radio Frequency Identification,RFID)在鐵路信號領(lǐng)域的特殊應(yīng)用,其有1個上行鏈路傳輸信道,工作頻率為4.234 MHz,采用頻移鍵控(Frequency Shift Keying,F(xiàn)SK)調(diào)制方式。雙信道應(yīng)答器是再增加1個上行鏈路傳輸信道,工作頻率為9.032 MHz,采用相移鍵控(Phase Shift Keying,PSK)調(diào)制方式。雙信道同時傳輸應(yīng)答器報文,將應(yīng)答器數(shù)據(jù)容量從1 023位擴(kuò)大為2 046位,可滿足既有線對鐵路信號大容量數(shù)據(jù)的需求。

        目前對應(yīng)答器的研究分別在應(yīng)用層、傳輸層和物理層3個層面上進(jìn)行。在應(yīng)用層面上的研究主要包括:應(yīng)答器的報文解析、測試[1-2]和布置[3]。在傳輸層面上的研究主要包括:電磁信號轉(zhuǎn)換為數(shù)字信號后的校驗、編碼、譯碼;定量評價應(yīng)答器傳輸模塊(Balise Transmission Module,BTM)的動態(tài)特性和高速適用性[4-5]。在物理層面上的研究主要包括:應(yīng)答器電磁特性和場強(qiáng)的分析[6];車載天線與地面應(yīng)答器之間射頻能量和數(shù)據(jù)傳輸過程的分析[7];從電磁場角度對上行鏈路和射頻能量傳輸過程的建模分析[8-9];周圍復(fù)雜空間電磁環(huán)境對應(yīng)答器“A”接口性能會產(chǎn)生影響,因此需要對周圍空間介質(zhì)進(jìn)行約束[10],例如應(yīng)答器附近的護(hù)軌侵入其無金屬區(qū)后需對其進(jìn)行截斷[11]。物理層面的研究是傳輸層面和應(yīng)用層面研究的基礎(chǔ)。但目前物理層面的研究很少關(guān)注應(yīng)答器的內(nèi)部結(jié)構(gòu)。

        彼此靠近的多個天線之間以一種復(fù)雜的方式相互作用,這種現(xiàn)象稱為互耦[12-14]。天線間的互耦,一方面會引起阻抗變化,從而破壞天線的調(diào)諧狀態(tài),使調(diào)諧頻率偏移;另一方面敏感天線收到源天線發(fā)射的高次諧波,會產(chǎn)生感應(yīng)電壓,成為帶內(nèi)干擾信號,影響敏感天線的正常工作。雙信道應(yīng)答器內(nèi)部4個天線間的互耦,使天線加工調(diào)試?yán)щy,產(chǎn)品合格率低。對于多天線間的去耦,Chi-Yuk Chiu提出在接地面刻槽,阻止平面波的傳播[15],但應(yīng)答器的信號波長很大,這種刻槽阻波的方法顯然不適合。Won-Kyu Choi提出的去耦方法是:在1個天線工作時,加載1個開關(guān),令相鄰天線均停止工作[16],但這個方法顯然也不適用。雙信道應(yīng)答器多天線間的降耦方法包括增加濾波器、優(yōu)化各天線的匹配電路、增加天線間距、調(diào)整天線間的排列順序、減小天線間的重合面積等。本文針對增加濾波器的方法,在物理層面上,采用HFSS(High Frequency Structure Simulator,高頻結(jié)構(gòu)仿真器)和CST(Computer Simulation Technology,計算機(jī)仿真)軟件,建立雙信道應(yīng)答器天線的三維仿真模型,分析應(yīng)答器內(nèi)部多天線間的互耦作用,以及增加濾波器后的降耦效果,從而為雙信道應(yīng)答器能夠大規(guī)模批量生產(chǎn)提供理論支持。

        1 雙信道應(yīng)答器天線的建模

        1.1 天線的等效電路及其帶寬計算公式

        天線有串聯(lián)和并聯(lián)2種等效電路,如圖1所示。圖中:L1,R2,C3分別為天線的串聯(lián)等效電感、串聯(lián)損耗電阻、串聯(lián)寄生電容;L2,R4,C6分別為天線的并聯(lián)等效電感、并聯(lián)損耗電阻、并聯(lián)寄生電容。寄生電容的形成原因是:線圈中流經(jīng)有高頻電流,在兩段平行的導(dǎo)線之間存在位移電流,兩段平行導(dǎo)線如同電容的2個極板。

        對于同一天線,串聯(lián)等效電路的阻抗Zws與并聯(lián)等效電路的阻抗Zwp應(yīng)該相等,即

        Zws=Zwp

        (1)

        由于天線主要呈現(xiàn)電感性質(zhì),寄生電容對天線的阻抗貢獻(xiàn)很小,所以可忽略寄生電容C3和C6。設(shè)ω為角頻率,故串聯(lián)、并聯(lián)等效電路的阻抗計算公式分別為

        Zws=R2+jωL1

        (2)

        (3)

        對于同一電路,串聯(lián)等效電路的質(zhì)量因子Qws與并聯(lián)等效電路的質(zhì)量因子Qwp應(yīng)該相等,即

        Qws=Qwp

        (4)

        串聯(lián)、并聯(lián)等效電路質(zhì)量因子的計算式分別為

        (5)

        (6)

        由式(3)和式(6)可得

        (7)

        由式(1)、式(2)和式(7)可得

        (8)

        對于式(8),根據(jù)實部和虛部分別對應(yīng)相等,得

        (9)

        (10)

        由式(9)得

        (11)

        由式(4)—式(6)和式(10)可得

        (12)

        可見:式(11)即為并聯(lián)損耗電阻R4與串聯(lián)損耗電阻R2的轉(zhuǎn)換關(guān)系式;式(12)即為并聯(lián)等效電感L2與串聯(lián)等效電感L1的轉(zhuǎn)換關(guān)系式。

        圖1 天線的2種等效電路

        線圈帶寬Bw的表達(dá)式為

        (13)

        式中:ω0為工作角頻率。

        由式(4)—式(6)和式(13)可得

        (14)

        (15)

        由式(15)可知,通過調(diào)整R1和R3,就可調(diào)整天線的帶寬BAntenna。

        同時,匹配電路除了可以用于調(diào)整天線帶寬,還有1個重要作用就是將天線調(diào)諧在工作頻率上,即調(diào)整其中的調(diào)諧電容C1和C2,C4和C5,使天線阻抗處于共軛匹配狀態(tài)。

        1.2 天線的仿真模型

        雙信道應(yīng)答器內(nèi)部有4個天線,自上而下分別是27.095 MHz射頻能量接收天線(簡稱為27.095 MHz接收天線)、4.234 MHz FSK上行鏈路發(fā)射天線(簡稱為4.234 MHz FSK發(fā)射天線)、9.032 MHz PSK上行鏈路發(fā)射天線(簡稱為9.032 MHz PSK發(fā)射天線)、9.032 MHz編程接收天線(簡稱為編程天線),其中編程天線與27.095 MHz接收天線位于同一平面內(nèi),如圖2所示。

        圖2 雙信道應(yīng)答器內(nèi)部天線位置示意圖

        采用HFSS軟件建立雙信道應(yīng)答器天線的三維仿真模型,如圖3所示,基于有限元法對多天線間的互耦和降耦進(jìn)行仿真分析。為了驗證降耦方法的有效性和合理性,采用CST軟件,基于有限積分法進(jìn)行交叉驗證。

        圖3 雙信道應(yīng)答器天線的仿真模型

        根據(jù)SUBSET-036《FFFIS for Eurobalise》規(guī)范[17]要求,雙信道應(yīng)答器的激活參考區(qū)域為200 mm×390 mm。 為了生成該激活參考區(qū)域,設(shè)置天線的模型參數(shù)為:除編程天線外,其他3個天線的長度均為420 mm,寬度均為230 mm;導(dǎo)線的寬度為1.5 mm;天線間的垂直距離為1 mm;編程天線為2個相同矩形組成的8字形天線,其中每個矩形的長度為90 mm,寬度為62 mm;電路板厚度為5 mm。仿真時采用圖1(a)中的天線串聯(lián)等效電路進(jìn)行仿真計算,所以下面的分析中,所有參數(shù)的名稱中均省略“串聯(lián)”二字。

        2 雙信道應(yīng)答器多天線間的互耦分析

        雙信道應(yīng)答器內(nèi)部多天線間的互耦,一方面會影響各天線的等效參數(shù),包括等效電感、損耗電阻和自諧振頻率;另一方面會使天線間產(chǎn)生轉(zhuǎn)移阻抗。等效參數(shù)的變化和轉(zhuǎn)移阻抗的產(chǎn)生都會導(dǎo)致天線間失諧,調(diào)諧頻率偏移。雙信道應(yīng)答器內(nèi)部多天線間互耦的程度采用S參數(shù)中的傳輸系數(shù)Sij(i,j=1,2,3,4;且i≠j)和反射系數(shù)Sij(i,j=1,2,3,4;且i=j)表征。

        2.1 天線的等效參數(shù)仿真

        在獨(dú)立空間中只有1個天線時,就沒有天線間的耦合作用。因此,首先對應(yīng)答器內(nèi)部只有4個天線中的1個天線時的情況分別進(jìn)行仿真,然后對有4個天線時的情況進(jìn)行仿真,對比2種情況下的等效參數(shù),分析其變化趨勢,由此得到多天線間的互耦影響。

        2.1.1應(yīng)答器內(nèi)只有1個天線時

        假設(shè)應(yīng)答器內(nèi)只有1個天線,分別仿真計算得到4個天線的阻抗圓圖;根據(jù)該阻抗圓圖得到天線沒有加調(diào)諧和濾波電路時在工作頻率下的50 Ω歸一化阻抗,將歸一化阻抗乘以50 Ω得到天線的阻抗,結(jié)果均見表1。

        表1 應(yīng)答器內(nèi)只有1個天線時的天線阻抗

        分析式(2)的組成可知:阻抗=損耗電阻+j電抗。因此,由表1對應(yīng)地可以得到天線的損耗電阻R2和電抗X。又因為

        X=ωL1=2πf0L1

        (16)

        式中:f0為工作頻率;L1為等效電感。

        由式(16)可得

        (17)

        根據(jù)表1得到天線的損耗電阻R2和電抗X,根據(jù)式(17)計算得到天線的等效電感L1,根據(jù)天線回波損耗得到自諧振頻率fsr。獨(dú)立空間中4個天線的等效電路參數(shù)見表2。

        2.1.2應(yīng)答器內(nèi)有4個天線時

        采用同樣的方法,當(dāng)應(yīng)答器內(nèi)有4個天線時,通過仿真得到多天線的阻抗和等效參數(shù),見表3和表4。

        表2 應(yīng)答器內(nèi)只有1個天線時的天線等效參數(shù)

        表3 應(yīng)答器內(nèi)有4個天線時的天線阻抗

        表4 應(yīng)答器內(nèi)有4個天線時的等效參數(shù)

        比較表2和表4可知:有4個天線時與只有1個天線時相比,因多天線間的互感,導(dǎo)致?lián)p耗電阻增加、等效電感和自諧振頻率減小??梢?,多天線間的互耦會引起天線間失諧、調(diào)諧頻率偏移。

        2.2 應(yīng)答器內(nèi)有4個天線時的轉(zhuǎn)移阻抗

        雙信道應(yīng)答器內(nèi)部4個天線可以等效為4端口網(wǎng)絡(luò),如圖4所示。圖4中:Z為端口阻抗,下角標(biāo)中的數(shù)字對應(yīng)端口編號。

        圖4 4端口網(wǎng)絡(luò)等效模型

        各天線間產(chǎn)生的轉(zhuǎn)移阻抗用矩陣表示,為

        (18)

        除j端口外其他3個端口均為開路,Zij(i,j=1,2,3,4;且i≠j)為j端口到i端口產(chǎn)生的轉(zhuǎn)移阻抗;除i端口外其他3個端口均為開路,Zij(i,j=1,2,3,4;且i=j)為i端口的自阻抗。采用圖3所示的仿真模型,仿真得到歸一化的轉(zhuǎn)移阻抗矩陣矩陣(對角線上為自阻抗)為

        由該轉(zhuǎn)移阻抗矩陣可得如下結(jié)論。

        (1)第1列元素中轉(zhuǎn)移阻抗最大的是Z21。這是因為對于端口1,端口2與其之間的距離為1 mm, 端口3與其之間的距離為2 mm,Z21大于Z31,說明天線間的距離影響了轉(zhuǎn)移阻抗。

        (2)第2列元素中最大的轉(zhuǎn)移阻抗是Z12。這是因為,端口1的工作頻率位于端口2工作頻率的高次諧波頻帶內(nèi),說明高次諧波會增加轉(zhuǎn)移阻抗。

        (3)第3列元素中最大的轉(zhuǎn)移阻抗是Z13。這是因為,端口1的工作頻率為27.095 MHz,位于端口3的工作頻率9.032 MHz的高次諧波頻帶內(nèi),因此其對端口1的轉(zhuǎn)移阻抗最大。

        (4)第4行元素中最大的轉(zhuǎn)移阻抗是Z34。這是因為端口4與端口3的工作頻率相同,所以4端口對3端口產(chǎn)生的轉(zhuǎn)移阻抗最大。

        2.3 應(yīng)答器內(nèi)有4個天線時的傳輸系數(shù)

        將27.095 MHz接收天線、4.234 MHz FSK發(fā)射天線、9.032 MHz PSK發(fā)射天線、編程天線分別編號為1,2,3,4。對于圖4所示的仿真模型,基于有限元法,采用HFSS軟件仿真計算應(yīng)答器內(nèi)4個天線間的S參數(shù),結(jié)果見表5。

        表5 多天線的S參數(shù) dB

        由表5可知:其中傳輸系數(shù)S12,S13,S23和S32的值均大于-10 dB[15]。這是因為:4.234 MHz FSK發(fā)射天線接收信號的上邊頻為4.516 MHz,其二次諧波頻率為9.032 MHz,六次諧波頻率為27.096 MHz。設(shè)9.032 MHz PSK發(fā)射天線和27.095 MHz接收天線的帶寬都為1 MHz,則其帶寬對應(yīng)的頻率范圍分別為(8.532 MHz,9.532 MHz)和(26.595 MHz,27.595 MHz)。因此,4.234 MHz FSK發(fā)射天線的信號高次諧波會落于9.032 MHz PSK發(fā)射天線和27.095 MHz接收天線的帶寬內(nèi),9.032 MHz PSK發(fā)射天線的信號高次諧波會落于27.095 MHz接收天線的帶寬內(nèi),成為帶內(nèi)高次諧波干擾,帶內(nèi)高次諧波引起天線間的耦合增大,使得傳輸系數(shù)S12,S13和S32均大于-10 dB;4.234 MHz FSK發(fā)射天線與9.032 MHz PSK發(fā)射天線之間距離僅為1 mm,且2個天線的工作頻率差較小,因此形成緊耦合,使得傳輸系數(shù)S23大于-10 dB。

        3 雙信道應(yīng)答器多天線間的降耦

        為了從源頭降低諧波干擾,在應(yīng)答器天線信號源端加載低通濾波器。由于編程天線和發(fā)射天線均沒有信號源,所以不用加載濾波器,4.234 MHz FSK發(fā)射天線和9.032 MHz PSK僅在發(fā)射天線的信號源處加載RC低通濾波器,濾除高次諧波。對于圖4所示的仿真模型,仍基于有限元法,采用HFSS軟件仿真增加濾波器后應(yīng)答器內(nèi)4個天線間的S參數(shù);為了醒目,僅列出S12,S13,S23和S32這4個參數(shù)濾波前、后的值,見表6。圖5顯示了濾波前后S參數(shù)隨頻率的變化關(guān)系,以圖5(a)和(a′)為例,在0~35 MHz頻率范圍內(nèi),S11在27.095 MHz處取得最小值,當(dāng)4個天線都采用27.095 MHz的激勵信號時,S12和S13在濾波后明顯下降。

        為了驗證采用濾波器降耦的有效性,仍然采用圖3所示的仿真模型,基于有限積分法,采用CST軟件仿真增加濾波器前、后應(yīng)答器內(nèi)4個天線間的S參數(shù);同樣,也僅列出S12,S13,S23和S32這4個傳輸系數(shù)濾波前、后的值,見表7。

        圖5 濾波前后的S參數(shù)

        表6基于有限元法仿真濾波前后其中的4個傳輸系數(shù)dB

        表7基于有限積分法仿真濾波前后其中的4個傳輸系數(shù)dB

        由表6和表7可知,濾波后這4個傳輸系數(shù)的值均降低到-10 dB以下,從而證明了降耦方法的有效性和合理性。說明增加RC低通濾波器后,位于敏感天線帶寬內(nèi)的高次諧波分量被濾除,干擾降低,天線間傳輸系數(shù)的值下降,從而提高了天線工作的可靠性。

        4 結(jié) 語

        雙信道應(yīng)答器多天線間互耦引起損耗電阻增加、等效電感和自諧振頻率下降;同時,互耦在多天線間產(chǎn)生轉(zhuǎn)移阻抗。等效參數(shù)變化和轉(zhuǎn)移阻抗引起天線阻抗變化,進(jìn)而使天線調(diào)諧頻率發(fā)生偏移。使用S參數(shù)中的傳輸系數(shù)表征互耦程度?;谟邢拊ㄓ嬎愕玫蕉嗵炀€間的傳輸系數(shù)S12,S13,S23和S32均大于-10 dB。為了降低多天線間的互耦,在4.234 MHz FSK發(fā)射天線和9.032 MHz PSK發(fā)射天線的信號源輸出端加載RC低通濾波器,濾除了天線帶寬內(nèi)的高次諧波,濾波后多天線間的傳輸系數(shù)均低于-10 dB。采用有限積分法仿真計算了增加濾波器前、后多天線間的傳輸系數(shù),進(jìn)一步驗證了增加濾波器降耦的有效性和合理性。

        [1]謝立鵬.應(yīng)答器測試關(guān)鍵技術(shù)與測試管理系統(tǒng)的研究 [D]. 北京:北京交通大學(xué), 2008.

        (XIE Lipeng. The Research of Crucial Technology and Management System for Balise Testing [D]. Beijing: Beijing Jiaotong University, 2008. in Chinese)

        [2]徐寧, 張季良, 王財進(jìn), 等. 列車高速運(yùn)行條件下應(yīng)答器車載測試設(shè)備關(guān)鍵技術(shù)[J]. 中國鐵道科學(xué), 2010, 31(4):131-137.

        (XU Ning, ZHANG Jiliang, WANG Caijin, et al. The Crucial Techniques of Balise Onboard Test Equipment for High Speed Train[J]. China Railway Science, 2010, 31(4):131-137. in Chinese)

        [3]CHEN Dewang, CHEN Rong, LI Yidong, et al. Online Learning Algorithms for Train Automatic Stop Control Using Precise Location Data of Balises[J]. IEEE Transactions on Intelligent Transportation Systems, 2013, 14(3):1526-1535.

        [4]趙會兵, 唐抗尼, 李偉, 等. 應(yīng)答器傳輸模塊的動態(tài)特性及高速條件下的適用性評價[J]. 中國鐵道科學(xué), 2010, 31(3):93-98.

        (ZHAO Huibing, TANG Kangni, LI Wei, et al. The Dynamic Characteristics of Balise Transmission Module and the Assessment on Its Adaptability under High Speed Condition[J]. China Railway Science, 2010, 31(3):93-98. in Chinese)

        [5]楊志杰, 范浦輝, 薛瑞民, 等. 適應(yīng)于高速運(yùn)營與提速的查詢應(yīng)答器系統(tǒng)[J]. 中國鐵道科學(xué), 2002, 23(2):42-47.

        (YANG Zhijie, FAN Puhui, XUE Ruimin, et al. Balise System Used in High Speed or Speed-Increase Line [J]. China Railway Science, 2002, 23(2):42-47. in Chinese)

        [6]SEMET C, MEGANCK P, GABILLARD R, et al. Feasibility and Security for a New Automatic Train Localization System by Electronic Beacon [C]//IEEE 46th Vehicular Technology Conference. Atlanta: IEEE, 1996:1604-1608.

        [7]梁迪, 趙會兵, 全宏宇, 等. 應(yīng)答器傳輸系統(tǒng)的電磁耦合機(jī)理及工程安裝優(yōu)化研究[J]. 鐵道學(xué)報, 2014, 36(5):64-70.

        (LIANG Di, ZHAO Huibing, QUAN Hongyu, et al. Research on Electromagnetic Coupling Mechanism and Mounting Parameter Optimization of Balise Transmission System[J]. Journal of the China Railway Society, 2014, 36(5):64-70. in Chinese)

        [8]ZHAO Linhai, JIANG Ying. Modenliing and Simulation of Balise Up-Link Data Transmission Based on Finite Element Method [J]. Journal of Theoretical & Applied Information Technology, 2012, 46(2):867-874.

        [9]ZHAO Linhai, JIANG Ying. Modeling and Optimization Research for Dynamic Transmission Process of Balise Tele-Powering Signal in High-Speed Railways[J]. Progress in Electromagnetics Research, 2013, 140:563-588.

        [10]趙會兵, 孟亞松, 梁迪. 空間環(huán)境介質(zhì)對應(yīng)答器 “A” 接口性能的影響研究[J]. 鐵道學(xué)報, 2015, 37(3):54-58.

        (ZHAO Huibing, MENG Yasong, LIANG Di. Study on Influence of Environmental Mediums on Performance of Interface ‘A’ of Balise System[J]. Journal of the China Railway Society, 2015, 37(3):54-58. in Chinese)

        [11]羅麗燕, 劉中田, 周果. 護(hù)軌對應(yīng)答器干擾的仿真研究[J]. 鐵道學(xué)報, 2014, 36(10):59-64.

        (LUO Liyan, LIU Zhongtian, ZHOU Guo. Simulation of Interference of Guardrail in Balise[J]. Journal of the China Railway Society, 2014, 36(10):59-64. in Chinese)

        [12]伍裕江, 聶在平. 一種新的互耦補(bǔ)償方法及其理論分析[J]. 電子學(xué)報, 2007, 35(3):492-496.

        (WU Yujiang, NIE Zaiping. A New Mutual Coupling Compensation Method and Its Theoretical Analysis[J]. Acta Electronica Sinica, 2007, 35(3):492-496. in Chinese)

        [13]謝俊, 胡皓全. 任意多邊形狀天線之間的互耦計算[J]. 電子科技大學(xué)學(xué)報, 1995, 24(增1):37-41.

        (XIE Jun, HU Haoquan. Research on the Mutual Coupling between Arbitrary Polygon Thin Linear Antenna[J]. Journal of UEST of China, 1995, 24(Supplement 1):37-41. in Chinese)

        [14]YAMADA Hiroyoshi, MORISHITA Masayuki, YAMAGUCHI Yoshio. A Simple Mutual Coupling Compensation Technique in Array of Single-Mode Elements by a Weighted Mutual Coupling Matrix Based on the Impedance Matrix[J]. IEICE Transactions on Communications, 2007, 90(9):2288-2296.

        [15]CHIU Chiyuk, CHENG Chiho, MURCH Ross D, et al. Reduction of Mutual Coupling between Closely-Packed Antenna Elements[J]. IEEE Transactions on Antennas and Propagation, 2007, 55(6):1732-1738.

        [16]CHOI Wonkyu, JEONG Seunghwan, PARK Chanwon, et al. HF RFID Reader Antenna with Loop Switch for Avoiding Magnetic Coupling [J]. ETRI Journal, 2014, 36(6):881-888.

        [17]Union Industry of Signalling. UNISIG SUBSET-036 Form Fit Function Interface Specification for Eurobalise[S]. Brussels: UNISIG Consortium, 2012.

        精品一区二区三区免费视频| 亚州韩国日本区一区二区片| 超碰青青草手机在线免费观看| 超级乱淫片国语对白免费视频| 国产免费艾彩sm调教视频| 精品熟女日韩中文十区| 中文字幕久久精品波多野结百度 | 午夜精品久久久久久久99老熟妇| 成 人 免费 黄 色 视频| 国产原创精品视频| 91九色精品日韩内射无| 国产人成视频在线视频| 中国国语毛片免费观看视频| 久久久亚洲经典视频| 美女人妻中文字幕av| 色综合久久网| 国产无遮挡裸体免费视频| 精品国产一区二区三区香蕉| 蜜桃在线视频一区二区| 亚洲av成人无码一区二区三区在线观看| 日韩无套内射视频6| 高清高速无码一区二区| 国产夫妻精品自拍视频| 久久av一区二区三区下 | 久久夜色精品国产亚洲av动态图| 久激情内射婷内射蜜桃人妖| av一区无码不卡毛片| 日韩精品有码中文字幕| av一区二区三区在线| 国产色在线 | 日韩| 国产熟妇搡bbbb搡bbbb搡| 手机在线中文字幕国产| 男性av天堂一区二区| 乱人妻中文字幕| 久久天天爽夜夜摸| 精品在线亚洲一区二区三区| 精品亚洲a∨无码一区二区三区| 欧美第一黄网免费网站| A亚洲VA欧美VA国产综合| 色和尚色视频在线看网站| s级爆乳玩具酱国产vip皮裤|