两个人的电影免费视频_国产精品久久久久久久久成人_97视频在线观看播放_久久这里只有精品777_亚洲熟女少妇二三区_4438x8成人网亚洲av_内谢国产内射夫妻免费视频_人妻精品久久久久中国字幕

電磁波與孔縫腔體內(nèi)復(fù)雜傳輸線耦合預(yù)測(cè)方法與流程

文檔序號(hào):11951040閱讀:564來源:國知局
電磁波與孔縫腔體內(nèi)復(fù)雜傳輸線耦合預(yù)測(cè)方法與流程
本發(fā)明針對(duì)雷達(dá)電磁波與孔縫腔體內(nèi)復(fù)雜傳輸線(TL)耦合預(yù)測(cè)問題,提出一種基于矩量法(MOM)、并矢格林函數(shù)(GF)和傳輸線PI等效電路與Agrawal混合模型的預(yù)測(cè)方法,屬于微波和電磁兼容領(lǐng)域。
背景技術(shù)
:隨著電子科學(xué)技術(shù)的發(fā)展,海戰(zhàn)場上的電磁環(huán)境逐漸變得惡劣,敏感設(shè)備越來越容易受電磁波干擾,房間內(nèi)雷達(dá)設(shè)施或電子箱內(nèi)設(shè)備也易受到干擾。房間或箱體等腔體由于通風(fēng)散熱口等需要,不可避免存在孔縫,這使得電磁波可通過腔體的孔縫耦合到腔體內(nèi)。故電磁場通過孔縫對(duì)腔體內(nèi)線纜的耦合響應(yīng)已經(jīng)是EMC領(lǐng)域的重要問題。很多學(xué)者已經(jīng)對(duì)此類問題作出大量的研究,主要分為數(shù)值方法、混合算法和解析算法。數(shù)值算法雖然準(zhǔn)確,但是耗費(fèi)很大,計(jì)算時(shí)間較長,不能滿足計(jì)算的快速性。有人提出一種基于電磁場下BLT方程求解傳輸線的負(fù)載響應(yīng)的混合算法,但是BLT混合算法有時(shí)不能求解復(fù)雜傳輸線的耦合響應(yīng)。本發(fā)明分析電磁波對(duì)房間這一大腔體內(nèi)的復(fù)雜分布的傳輸線耦合,用MOM求得雷達(dá)天線產(chǎn)生的入射電磁場;采用MOM-GF混合算法計(jì)算孔縫處等效磁流并求得腔體內(nèi)的電磁場;接著運(yùn)用所得電磁場與傳輸線混合模型計(jì)算傳輸線的負(fù)載電流響應(yīng)。技術(shù)實(shí)現(xiàn)要素:本發(fā)明針對(duì)雷達(dá)電磁波與孔縫腔體內(nèi)復(fù)雜傳輸線耦合預(yù)測(cè)問題,提出一種基于矩量法(MOM)、并矢格林函數(shù)和傳輸線等效電路與Agrawal混合模型的預(yù)測(cè)方法。本發(fā)明的目的是這樣實(shí)現(xiàn)的:電磁波與孔縫腔體內(nèi)復(fù)雜傳輸線耦合預(yù)測(cè)方法包括如下步驟:步驟1.運(yùn)用MOM計(jì)算雷達(dá)天線入射電磁場;步驟2.利用MOM-GF分析孔縫腔體內(nèi)電磁場;步驟3.對(duì)步驟2所述孔縫腔體內(nèi)任意截面形狀的多導(dǎo)體傳輸線電容、電感參數(shù)計(jì)算;步驟4.建立PI等效電路與Agrawal混合模型對(duì)傳輸線分析。步驟1所述計(jì)算電磁場的方法如下:(1):利用雷達(dá)天線表面積分方程建立雷達(dá)天線的電場模型;(2):用RWG基函數(shù)將目標(biāo)電流展開得到線性方程組:[Z][I]=[V]其中,Z為阻抗矩陣,I為待求電流系數(shù),V為激勵(lì)向量;(3):根據(jù)步驟(1)所述的電場模型采用步驟(2)所述的電流系數(shù)I得到雷達(dá)天線的輻射電場E;(4):根據(jù)步驟(3)所述的雷達(dá)天線輻射場電場求得磁場強(qiáng)度H:H=1η·r^×E]]>其中,為雷達(dá)天線遠(yuǎn)區(qū)某點(diǎn)的坐標(biāo)向量;η≈377。步驟2所述孔縫腔體內(nèi)電磁場的分析方法如下:1):用脈沖基函數(shù)將孔縫分解為P(x)×Q(y)塊面元,得到孔縫每個(gè)面元中心點(diǎn)的等效磁流;2):根據(jù)步驟1)所述的孔縫每個(gè)面元中心點(diǎn)的等效磁流得到孔縫腔體內(nèi)任意一點(diǎn)的電場:Exb=Σm,n[jZ0k0kxkzAmn+kyBmn]cos(kxx)sin(kyy)sin(kz·(z+c))]]>Eyb=Σm,n[jZ0k0kxkzAmn-kyBmn]sin(kxx)cos(kyy)sin(kz·(z+c))]]>Ezb=-Σm,njZ0kc2k0Amnsin(kxx)sin(kyy)cos(kz·(z+c))]]>其中,ϵn=1,n=02,n>0]]>Ixmn=∫sMxsin(kxx)cos(kyy)ds]]>Iymn=∫sMycos(kxx)sin(kyy)ds]]>kx=mπ/a;ky=nπ/b;kc2=kx2+ky2;]]>Z0為自由空間波阻抗;j為虛數(shù)單位;k0為自由空間波數(shù);Y0為自由空間導(dǎo)納;a,b和c分別為腔體的長、高和寬;Mx和My分別孔縫處x和y方向的等效磁流。步驟3所述多導(dǎo)體傳輸線電容、電感參數(shù)計(jì)算方法如下:運(yùn)用傳輸線PI等效電路模型將n+1傳輸線多導(dǎo)體系統(tǒng)第i個(gè)導(dǎo)體傳輸線劃分為Ni段,每一段的傳輸線通過矩量法求得電容、互容、自感和互感,其中,0<i≤n+1,i和n均為大于1的正整數(shù)。步驟4所述傳輸線分析方法如下:首先,運(yùn)用傳輸線PI等效電路模型將傳輸線劃分為N段,N為大于1的自然數(shù);其次,將傳輸線網(wǎng)絡(luò)分為多個(gè)多導(dǎo)體和單導(dǎo)體的集合,而電磁場對(duì)多導(dǎo)體和單導(dǎo)體傳輸線的耦合電壓采用Agrawal模型求出;最后,根據(jù)耦合電壓輸出最后的電流響應(yīng)。電磁波波與孔縫內(nèi)復(fù)雜傳輸線耦合計(jì)算問題,單獨(dú)運(yùn)用全波算法耗費(fèi)大、時(shí)間長,本發(fā)明提出一種基于矩量法、并矢格林函數(shù)和傳輸線PI等效電路與Agrawal混合模型的計(jì)算方法。該算法適用多導(dǎo)體的情況,對(duì)于解決電磁波通過孔縫耦合到其他敏感設(shè)備的分析具有很強(qiáng)的優(yōu)勢(shì)。附圖說明圖1:電磁波與孔縫腔體內(nèi)多導(dǎo)體束傳輸線耦合模型;其中,A為雷達(dá)天線簡化模型。圖2:傳輸線第i個(gè)導(dǎo)體的截面圖;圖3:傳輸線PI等效電路模型。具體實(shí)施方式下面該方法作進(jìn)一步的詳細(xì)說明。步驟1.運(yùn)用MOM計(jì)算雷達(dá)天線入射電磁場圖1所示為一個(gè)電磁波與孔縫腔體內(nèi)傳輸線耦合模型圖。為便于計(jì)算,將雷達(dá)天線A用線天線代替,其長度L0、半徑r0,激勵(lì)源在天線中心,位置為P0(x0,y0,z0)。腔體是一個(gè)體積為a(x)×b(y)×c(z)的矩形腔體且默認(rèn)腔體內(nèi)外都為自由空間;孔縫也為矩形,面積等于l(x)×w(y),其近原點(diǎn)坐標(biāo)為(xa,yb,0)。傳輸線的長L、半徑r、首末段負(fù)載Z1Z2。首先,運(yùn)用矩量法計(jì)算位于P(x,y,0)的孔縫中心處的雷達(dá)電磁場。理想導(dǎo)體構(gòu)成的天線表面電場積分方程:其中,為導(dǎo)體表面在點(diǎn)的單位矢量;為輻射電場;μ0為真空磁導(dǎo)率,ε0為真空介電值常數(shù);為點(diǎn)的表面電流;k為波數(shù);為并矢和為單位坐標(biāo)向量;g(r,r′)=e-ik|r-r′|/|r-r′|;j為虛數(shù)單位。用RWG基函數(shù)將目標(biāo)電流展開得到線性方程組:[Z][I]=[V](2)其中,Z為阻抗矩陣,I為待求電流系數(shù),V為激勵(lì)向量。阻抗矩陣的表達(dá)式為:其中,為第m條邊對(duì)應(yīng)的RWG基函數(shù);ω=2πf;ε、μ分別為電介質(zhì)常數(shù)和磁導(dǎo)率。將求得的電流帶入電場積分方程(1)可得入射電場根據(jù)理想導(dǎo)體的邊界條件:可得天線的輻射電場。運(yùn)用線天線遠(yuǎn)區(qū)輻射場電場強(qiáng)度與磁場強(qiáng)度關(guān)系可以求得磁場強(qiáng)度:H=1η·r^×E---(4)]]>其中,為遠(yuǎn)區(qū)某點(diǎn)的坐標(biāo)向量。步驟2.利用MOM-GF分析腔體內(nèi)電磁場為求孔縫腔體內(nèi)的電磁場,根據(jù)Schelkunoff場等效原理,滿足內(nèi)外等效電流切向連續(xù),在孔縫處引入等效磁流M(忽略傳輸線對(duì)其的影響)。腔體外的a區(qū)和腔體內(nèi)b區(qū)等效磁流分別滿足:其中,為入射電場;為磁流M在腔體外產(chǎn)生的電場;為腔體內(nèi)產(chǎn)生的電場。根據(jù)公式(5)、(6)可知:其中,為與對(duì)應(yīng)的磁場;為腔體內(nèi)部磁場。其中,k0是自由空間波數(shù);Y0為自由空間導(dǎo)納;此處s為孔縫范圍,而且為二維格林函數(shù),只與x、y有關(guān)。其中,依據(jù)公式(7)-(9)可得:為求解此公式,用脈沖基函數(shù)將孔縫分解為P(x)×Q(y)塊面元,即展開M,得到線性方程組:[Ya+Yb]·[M]=[C](14)其中,[M]的元素為包含位置編號(hào)為(p,q)的面元上磁流的x、y分量,即Mxpq和Mypq[12,15];[C]中元素為:其中,spq為編號(hào)(p,q)的面元范圍。Ya、Yb分別為a和b區(qū)的廣義導(dǎo)納,由于等效磁流有x、y分量,故Ya、Yb有四個(gè)分量,即和求得每個(gè)面元中心點(diǎn)的等效磁流后,近似認(rèn)為每個(gè)小面元上的等效磁流相等。運(yùn)用下列公式可得腔體內(nèi)任意一點(diǎn)的電場:Exb=Σm,n[jZ0k0kxkzAmn+kyBmn]cos(kxx)sin(kyy)sin(kz·(z+c))---(16)]]>Eyb=Σm,n[jZ0k0kxkzAmn-kyBmn]sin(kxx)cos(kyy)sin(kz·(z+c))---(17)]]>Ezb=-Σm,njZ0kc2k0Amnsin(kxx)sin(kyy)cos(kz·(z+c))---(18)]]>其中,ϵn=1,n=02,n>0---(20)]]>Ixmn=∫sMxsin(kxx)cos(kyy)dsIymn=∫sMycos(kxx)sin(kyy)ds---(21)]]>Z0為自由空間波阻抗;j為虛數(shù)單位;k0為自由空間波數(shù);Y0為自由空間導(dǎo)納;Mx和My分別孔縫處x和y方向的等效磁流。步驟3.任意截面形狀的多導(dǎo)體傳輸線電容、電感參數(shù)計(jì)算對(duì)于平行均勻的傳輸線,每Δx長度段的單位長度電感和電容相同;對(duì)于非平行或非均勻的傳輸線每段是不同的,需要分別計(jì)算。傳輸線第i個(gè)導(dǎo)體的截面圖如圖2所示,是非圓的非規(guī)則截面,其沿著x軸方向放置。其電勢(shì)與電荷密度滿足積分形式的泊松方程:φi(y,z)=14πϵ∫∫siρi(y,z)r---(22)]]>式中,φi(y,z)為某點(diǎn)產(chǎn)生的電勢(shì);r為導(dǎo)體上電荷分布的某點(diǎn)與產(chǎn)生電勢(shì)點(diǎn)的距離;ρi(y,z)C/m2為電荷密度;si為導(dǎo)體截面;ε為周圍介質(zhì)介電常數(shù)。對(duì)于非規(guī)則截面的傳輸線,無法直接求單位長度電容的解析解,在這里可以利用矩量法(MOM)求近似解。首先將導(dǎo)體表面電荷分布近似成沿著x軸劃分為Ni個(gè)細(xì)絲電荷,每段電荷密度為ρik。根據(jù)傳輸線導(dǎo)體之間間距與導(dǎo)體半徑之比大于4及以上時(shí),導(dǎo)體表面電荷分布可視為近似均勻分布(對(duì)于不規(guī)則截面,可認(rèn)為間距與最大半徑之比大于4及以上),假設(shè)在導(dǎo)體表面存在均勻分布的未知的電荷密度ρi(y,z)(單位C/m2),那么將導(dǎo)體全部電荷分布表現(xiàn)為基函數(shù)ρik的線性組合:ρi(y,z)=αi1ρi1(y,z)+αi2ρi2(y,z)+...+αiNiρiNi(y,z)=Σk=1Niαikρik(y,z)---(23)]]>式中,αik為電荷分布的具體水平,可以由第i個(gè)導(dǎo)體上的電勢(shì)為φi的邊界條件確定。x方向單位長度總電荷可以由該導(dǎo)體上各個(gè)子部分的電荷之和得到:qi=Σk=1Niαik∫ciρikdci---(24)]]>式中,ci是圖2中導(dǎo)體截面邊界。用脈沖法展開公式(12),簡化為:式中,是第i個(gè)導(dǎo)體中第k個(gè)子段的寬度。在第i個(gè)導(dǎo)體中,任意一點(diǎn)的電勢(shì)都是所有未知電荷分布的貢獻(xiàn),可以將它們的關(guān)系寫成如下的公式:φi(y,z)=ζi1αi1+ζi2αi2+...+ζiNiαiNi=Σk=1Niζikαik---(26)]]>式中,ζik是每個(gè)基函數(shù)單獨(dú)對(duì)電勢(shì)的貢獻(xiàn),即其值大小可由下列公式計(jì)算可得:dφi=-ρiln(r)2πϵdci⇒φi=-∫ciρiln(r)2πϵdci---(27)]]>式中,dφi為每個(gè)字段中心處的電勢(shì)。對(duì)于一個(gè)n+1傳輸線多導(dǎo)體系統(tǒng),每個(gè)導(dǎo)體上的電勢(shì)是系統(tǒng)中全部電荷產(chǎn)生的電勢(shì),因此可得導(dǎo)體i的電勢(shì)為:φi=Σi=0nΣk=1Niζikαik---(28)]]>其中,0<i≤n+1,n為大于1的正整數(shù)。根據(jù)展開系數(shù),可以得到聯(lián)立方程:Φ0...Φi...Φn=D00...D0j...D0n...............Di0...Dij...Din...............Dn0...Dnj...DnnA0...Ai...An---(29)]]>⇒A0...Ai...An=B00...B0j...B0n...............Bi0...Bij...Bin...............Bn0...Bnj...BnnΦ0...Φi...Φn---(30)]]>式中,Φi為第i個(gè)導(dǎo)體上的匹配點(diǎn)處的電勢(shì)向量,值為(φi…φi…φi)T是一個(gè)Ni×1矩陣;Ai第i個(gè)導(dǎo)體上電荷分布的展開系數(shù)向量,其也是一個(gè)Ni×1的矩陣,值為B=D-1。由公式(18)得到:Ai=BijΦj(31)電荷分布qi、電勢(shì)φi與廣義電容矩陣C的關(guān)系,可以計(jì)算傳輸線之間的自容和互容:由公式(13)、(19)和(20)可以得到廣義電容矩陣中元素的大?。菏街校绻磳?dǎo)體i的各子段的寬度相同,記為那么公式(21)簡化為:式中,∑Bij的值是子矩陣Bij所有元素之和。假設(shè)導(dǎo)體之間充滿介電常數(shù)為εr的介質(zhì)或非均勻介質(zhì):含介質(zhì)的單位長度電容為C;不含介質(zhì)的單位長度電容為C0;單位長度電感為L。它們之間的關(guān)系:均勻介質(zhì):非均勻介質(zhì):式中:εr,eff為有效介電常數(shù)。步驟4.建立PI等效電路與Agrawal混合模型對(duì)傳輸線分析首先,運(yùn)用傳輸線PI等效電路模型將傳輸線劃分為N段,如圖3所示,N=Nl+NL+Nr,Nl、Nr為首端和末端接地的分段數(shù),NL為傳輸線的分段數(shù);ΔL為每個(gè)分段的長度;li1,i1,i、ci1,i1,i分別為第i1個(gè)傳輸線第i段的單位自感和自容,li1,i2,i、ci1,i2,i分別為第i1個(gè)i2和傳輸線第i段的單位互感和互容,其值由步驟3可得,由此可得第i段的阻抗Zli和導(dǎo)納Yci:Zli=j(luò)ωliΔL;Yci=j(luò)ωciΔL。由此可得所有傳輸線第i段阻抗Zli的元素:Zi1,i1,i=j(luò)ωli1,i1,iΔL;Zi1,i2,i=j(luò)ωli1,i2,iΔL。所有傳輸線第i段阻抗Yci的元素:Yi1,i1,i=1/jωci1,i1,iΔL;Yi1,i2,i=1/jωci1,i2,iΔL。將復(fù)雜的傳輸線網(wǎng)絡(luò)分為多個(gè)多導(dǎo)體和單導(dǎo)體的集合,而電磁場對(duì)多導(dǎo)體和單導(dǎo)體傳輸線的耦合電壓可用Agrawal模型求出[13]。在首端、末端的電壓受電場Ey激勵(lì),其值可由下列公式求出:Vli=∫yli-1yliEy(xli,y,zli)dy;Vri=-∫yri-1yriEy(xri,y,zri)dy---(37)]]>Vxl(i1,i2,i)=∫xli1xli2Ex(x,yli,zli)dxVzl(i1,i2,i)=∫zli1zli2Ez(x,yli,zli)dx⇒Vdl(i1,i2,i)=V2xl(i1,i2,i)+V2zl(i1,i2,)---(38)]]>其中,Vli、Vri為首末端第i段的自感電壓;xliylizlixriyri和zri分別為首末段第i段的位置坐標(biāo);為首端第i1和第i2個(gè)傳輸線第i段的互感電壓,末端同首端方法可求。為求傳輸線上的耦合電壓,首先將傳輸線上每一段分解為x、y和z軸的平方和的二次根,即:分別求x、y和z軸方向電壓,最后再合成總電壓。VxLi=∫xLi-1xLiEx(x,yLi,zLi)dx≈Ex(x,yLi,zLi)·ΔxLi---(39)]]>VyLi=∫yLi-1yLiEy(xLi,y,zLi)dy≈Ey(xLi,y,zLi)·ΔyLi---(40)]]>VzLi=∫zLi-1zLiEz(xLi,yLi,z)dz≈Ez(xLi,yLi,z)·ΔzLi---(41)]]>VLi=VxLi2+VyLi2+VzLi2---(42)]]>其中,VLi分別為第i段的x、y和z軸自感方向電壓與自感總電壓?;ジ须妷捍笮∏蠼夥椒ㄍ?5.34),將x、y和z軸方向電壓合成總電壓如同式(5.35)-(5.38)。為方便表示,將首段到末端每段電壓用向量V1,V2…Vi…VN-1,VN編號(hào),每段電流用I1、I2…Ii…IN-1、IN編號(hào),Vi、Ii由n個(gè)傳輸線的第i段自電壓電流和互電壓電流構(gòu)成。通過電壓與電流的關(guān)系可得電流的大?。浩渲校珽i為單位向量,如果第i段時(shí)有傳輸線i1到達(dá)末端,則Ei第i段之后第i1個(gè)位置值為0。當(dāng)前第1頁1 2 3 
當(dāng)前第1頁1 2 3 
網(wǎng)友詢問留言 已有0條留言
  • 還沒有人留言評(píng)論。精彩留言會(huì)獲得點(diǎn)贊!
1
一区二区日韩欧美中文字幕 | 大香蕉久久成人网| 国产 一区精品| 国产 一区精品| 中国三级夫妇交换| 欧美性感艳星| 精品国产国语对白av| 人妻系列 视频| 如何舔出高潮| 免费人成在线观看视频色| 久久狼人影院| 亚洲精品久久成人aⅴ小说 | 亚洲美女黄色视频免费看| www.色视频.com| 欧美另类一区| 国产免费视频播放在线视频| 九色亚洲精品在线播放| 亚洲国产精品国产精品| 欧美激情 高清一区二区三区| av.在线天堂| 免费人妻精品一区二区三区视频| 国产极品粉嫩免费观看在线 | 韩国av在线不卡| 在线亚洲精品国产二区图片欧美 | 在线观看国产h片| 久久 成人 亚洲| 校园人妻丝袜中文字幕| 丰满乱子伦码专区| 久久久国产一区二区| 亚洲国产精品一区三区| 久久97久久精品| 国产亚洲最大av| 久久人人爽人人爽人人片va| 久久久久久久久大av| 精品久久久久久电影网| 自拍欧美九色日韩亚洲蝌蚪91| 国产亚洲最大av| 啦啦啦视频在线资源免费观看| 免费日韩欧美在线观看| 国产欧美日韩一区二区三区在线 | 如日韩欧美国产精品一区二区三区 | 久久久久久久国产电影| 国产 精品1| 性色av一级| 亚洲不卡免费看| 日本欧美视频一区| 精品少妇黑人巨大在线播放| 男女边摸边吃奶| 97在线视频观看| 亚洲欧美一区二区三区国产| 麻豆成人av视频| 黄色配什么色好看| 久久99蜜桃精品久久| 纵有疾风起免费观看全集完整版| 人妻夜夜爽99麻豆av| 婷婷色综合大香蕉| 美女视频免费永久观看网站| 亚洲国产精品999| 国产精品国产三级专区第一集| 国产精品免费大片| 97精品久久久久久久久久精品| 在线观看免费视频网站a站| 国产免费福利视频在线观看| 欧美另类一区| 热re99久久精品国产66热6| 久久亚洲国产成人精品v| 婷婷色麻豆天堂久久| 又粗又硬又长又爽又黄的视频| 国语对白做爰xxxⅹ性视频网站| 三上悠亚av全集在线观看| 久久久久久久久久成人| 国产国语露脸激情在线看| 亚洲,一卡二卡三卡| 中文字幕制服av| 人妻制服诱惑在线中文字幕| 日本wwww免费看| 亚洲高清免费不卡视频| av在线app专区| 日韩,欧美,国产一区二区三区| 成人黄色视频免费在线看| 久久青草综合色| 亚洲国产欧美日韩在线播放| 十八禁网站网址无遮挡| 七月丁香在线播放| 51国产日韩欧美| 日韩成人伦理影院| 日韩不卡一区二区三区视频在线| 在现免费观看毛片| 麻豆成人av视频| 美女主播在线视频| 日韩成人av中文字幕在线观看| 欧美3d第一页| 国产av码专区亚洲av| 久久久久视频综合| 麻豆精品久久久久久蜜桃| 中国美白少妇内射xxxbb| a级毛片在线看网站| 国产视频首页在线观看| 中文字幕免费在线视频6| 久久女婷五月综合色啪小说| 日韩av免费高清视频| 久久久久国产精品人妻一区二区| 欧美激情 高清一区二区三区| 男人操女人黄网站| 一本久久精品| 日日摸夜夜添夜夜添av毛片| 欧美bdsm另类| 亚洲四区av| 99久久精品一区二区三区| 天堂中文最新版在线下载| 18禁裸乳无遮挡动漫免费视频| 国产亚洲精品久久久com| 日韩精品免费视频一区二区三区 | 熟女电影av网| 丰满乱子伦码专区| 中文字幕免费在线视频6| 99精国产麻豆久久婷婷| 国产一区亚洲一区在线观看| 中文字幕av电影在线播放| 国产深夜福利视频在线观看| 亚洲在久久综合| 最新的欧美精品一区二区| 久久综合国产亚洲精品| 日日啪夜夜爽| 国产精品国产三级专区第一集| 国产综合精华液| 国产爽快片一区二区三区| 一级毛片黄色毛片免费观看视频| 精品久久久久久久久av| 如何舔出高潮| 两个人的视频大全免费| av天堂久久9| 欧美日韩亚洲高清精品| 女人久久www免费人成看片| 精品一区二区免费观看| 久久人人爽人人片av| 不卡视频在线观看欧美| 国产精品人妻久久久影院| 下体分泌物呈黄色| 成人免费观看视频高清| 欧美国产精品一级二级三级| 18禁观看日本| a级毛片在线看网站| 久久久久久久精品精品| 高清在线视频一区二区三区| 日韩av在线免费看完整版不卡| 日本爱情动作片www.在线观看| 新久久久久国产一级毛片| 亚洲精品aⅴ在线观看| 丰满迷人的少妇在线观看| 熟妇人妻不卡中文字幕| 亚洲av国产av综合av卡| 日韩免费高清中文字幕av| 一区二区日韩欧美中文字幕 | 亚洲天堂av无毛| 久久国产亚洲av麻豆专区| 欧美精品高潮呻吟av久久| 一级,二级,三级黄色视频| 久久99热这里只频精品6学生| 国产乱人偷精品视频| 亚洲av日韩在线播放| 亚洲情色 制服丝袜| 亚洲不卡免费看| 人成视频在线观看免费观看| 欧美性感艳星| 99热这里只有精品一区| 午夜福利网站1000一区二区三区| 久久人人爽av亚洲精品天堂| av线在线观看网站| 免费人妻精品一区二区三区视频| 2021少妇久久久久久久久久久| 亚洲欧美成人精品一区二区| 美女国产视频在线观看| 久久精品久久久久久久性| 人人妻人人澡人人爽人人夜夜| 在线观看三级黄色| 99热这里只有精品一区| 亚洲精品自拍成人| 我要看黄色一级片免费的| 亚洲精华国产精华液的使用体验| 成人国产av品久久久| 高清黄色对白视频在线免费看| 搡老乐熟女国产| 亚洲不卡免费看| 美女大奶头黄色视频| 最近中文字幕高清免费大全6| 国产亚洲av片在线观看秒播厂| 午夜日本视频在线| 我要看黄色一级片免费的| 亚洲国产精品999| 美女国产视频在线观看| 飞空精品影院首页| 伊人亚洲综合成人网| 啦啦啦中文免费视频观看日本| 老女人水多毛片| freevideosex欧美| 在线 av 中文字幕| 国产在线免费精品| 亚洲精品国产色婷婷电影| 久久久久久久大尺度免费视频| 大话2 男鬼变身卡| 王馨瑶露胸无遮挡在线观看| 欧美+日韩+精品| 亚洲国产欧美在线一区| 国产欧美日韩一区二区三区在线 | 如何舔出高潮| 天天操日日干夜夜撸| av在线播放精品| 久久婷婷青草| 亚洲一级一片aⅴ在线观看| 国产无遮挡羞羞视频在线观看| 夜夜看夜夜爽夜夜摸| 丰满迷人的少妇在线观看| 91成人精品电影| 最新的欧美精品一区二区| 日韩欧美精品免费久久| 久久久久国产网址| 亚洲情色 制服丝袜| 久久毛片免费看一区二区三区| 久久精品久久久久久噜噜老黄| 熟女人妻精品中文字幕| 久久国产精品大桥未久av| 日韩强制内射视频| 欧美97在线视频| 少妇被粗大猛烈的视频| 一级黄片播放器| 国产有黄有色有爽视频| 日日摸夜夜添夜夜爱| av专区在线播放| 国产一区二区在线观看av| 97精品久久久久久久久久精品| 久久久久久久久久人人人人人人| 啦啦啦中文免费视频观看日本| 久久狼人影院| 久久久久视频综合| 亚洲少妇的诱惑av| 99热国产这里只有精品6| 人妻少妇偷人精品九色| 久久久国产精品麻豆| 精品一区在线观看国产| 免费少妇av软件| 少妇被粗大猛烈的视频| 在线观看www视频免费| 狠狠婷婷综合久久久久久88av| 中文字幕最新亚洲高清| 日韩一区二区视频免费看| 草草在线视频免费看| 国产日韩欧美视频二区| 久久热精品热| 亚洲精品日本国产第一区| 久久精品人人爽人人爽视色| 亚洲色图综合在线观看| 97超碰精品成人国产| 国产精品人妻久久久影院| 中文字幕免费在线视频6| 2022亚洲国产成人精品| 国产黄色视频一区二区在线观看| 中文字幕人妻丝袜制服| 国模一区二区三区四区视频| 成人亚洲精品一区在线观看| 黄色怎么调成土黄色| 青春草视频在线免费观看| 日韩制服骚丝袜av| 婷婷色av中文字幕| 在线观看免费日韩欧美大片 | 91国产中文字幕| 女性被躁到高潮视频| 国产精品一区二区三区四区免费观看| 国产片内射在线| 精品久久久噜噜| 久久久久精品久久久久真实原创| 欧美3d第一页| 日本与韩国留学比较| 久久久久网色| 制服丝袜香蕉在线| 国产亚洲精品久久久com| 亚洲中文av在线| 亚洲情色 制服丝袜| h视频一区二区三区| 亚洲怡红院男人天堂| 久久99热6这里只有精品| 永久免费av网站大全| 一区二区三区乱码不卡18| 欧美激情极品国产一区二区三区 | 日本黄色片子视频| 午夜日本视频在线| 少妇 在线观看| 国产av码专区亚洲av| 国产成人免费无遮挡视频| 午夜福利,免费看| 欧美另类一区| 少妇的逼水好多| 亚洲国产欧美在线一区| 国产日韩欧美在线精品| 18禁裸乳无遮挡动漫免费视频| 内地一区二区视频在线| 日韩伦理黄色片| 国产69精品久久久久777片| 纵有疾风起免费观看全集完整版| 老司机亚洲免费影院| 成人免费观看视频高清| 美女福利国产在线| 女性生殖器流出的白浆| 久久国产亚洲av麻豆专区| 99re6热这里在线精品视频| 激情五月婷婷亚洲| 另类亚洲欧美激情| 精品一区在线观看国产| 18+在线观看网站| 有码 亚洲区| 亚洲av成人精品一区久久| 日韩不卡一区二区三区视频在线| 自线自在国产av| 精品久久久久久久久亚洲| 韩国av在线不卡| 天天躁夜夜躁狠狠久久av| 一本大道久久a久久精品| 国产男女超爽视频在线观看| 一区二区三区免费毛片| 美女脱内裤让男人舔精品视频| 国产黄片视频在线免费观看| 午夜福利,免费看| 欧美少妇被猛烈插入视频| 免费人成在线观看视频色| 国产一区亚洲一区在线观看| 黄色欧美视频在线观看| 一级片'在线观看视频| 日本91视频免费播放| 国产一区二区在线观看av| 亚洲人成网站在线播| 国产 一区精品| 性色av一级| 欧美日韩成人在线一区二区| 啦啦啦啦在线视频资源| 成人黄色视频免费在线看| 妹子高潮喷水视频| 婷婷色综合大香蕉| 亚洲美女搞黄在线观看| 色婷婷av一区二区三区视频| 高清在线视频一区二区三区| 亚洲一级一片aⅴ在线观看| 十八禁高潮呻吟视频| av免费观看日本| 成人二区视频| 欧美日韩精品成人综合77777| 久久久精品区二区三区| 久久99蜜桃精品久久| 国产视频内射| 国产国拍精品亚洲av在线观看| 高清在线视频一区二区三区| 日韩大片免费观看网站| 少妇精品久久久久久久| 亚州av有码| 亚洲av成人精品一区久久| 亚州av有码| 亚洲av成人精品一区久久| 日韩中文字幕视频在线看片| 男人添女人高潮全过程视频| 九色成人免费人妻av| 国产色婷婷99| 黑人猛操日本美女一级片| 日韩一区二区视频免费看| 免费久久久久久久精品成人欧美视频 | a级毛色黄片| 蜜臀久久99精品久久宅男| 日本黄色视频三级网站网址 | 精品国内亚洲2022精品成人 | 亚洲一码二码三码区别大吗| 手机成人av网站| 国产成人精品久久二区二区91| 国产精品熟女久久久久浪| 久久人妻熟女aⅴ| 国产精品熟女久久久久浪| 男女免费视频国产| 亚洲午夜精品一区,二区,三区| 最近最新中文字幕大全免费视频| a级毛片黄视频| 夜夜爽天天搞| 日本黄色视频三级网站网址 | 中文字幕av电影在线播放| 国产一区二区激情短视频| 国产免费视频播放在线视频| 午夜福利在线免费观看网站| 中文字幕高清在线视频| 国产精品av久久久久免费| 高清欧美精品videossex| 天堂中文最新版在线下载| 在线 av 中文字幕| 又大又爽又粗| 飞空精品影院首页| 视频区图区小说| 国产精品影院久久| 老司机午夜福利在线观看视频 | 日本vs欧美在线观看视频| 十八禁高潮呻吟视频| 18禁裸乳无遮挡动漫免费视频| 亚洲第一av免费看| 国产精品欧美亚洲77777| 18禁国产床啪视频网站| 91大片在线观看| 亚洲成av片中文字幕在线观看| 日本av手机在线免费观看| 19禁男女啪啪无遮挡网站| 成人精品一区二区免费| 99国产综合亚洲精品| 久久精品亚洲精品国产色婷小说| 男女下面插进去视频免费观看| 国产精品98久久久久久宅男小说| 国产色视频综合| 飞空精品影院首页| 国产成人免费观看mmmm| 久久精品熟女亚洲av麻豆精品| 高清在线国产一区| 亚洲熟妇熟女久久| 91九色精品人成在线观看| 黄频高清免费视频| 又紧又爽又黄一区二区| 久9热在线精品视频| 亚洲专区中文字幕在线| 国产单亲对白刺激| 国产主播在线观看一区二区| svipshipincom国产片| 精品乱码久久久久久99久播| 丁香六月欧美| 免费观看av网站的网址| 水蜜桃什么品种好| 国产精品成人在线| 丝瓜视频免费看黄片| 色婷婷av一区二区三区视频| av欧美777| 最近最新中文字幕大全电影3 | 满18在线观看网站| 九色亚洲精品在线播放| 一区在线观看完整版| 乱人伦中国视频| 91麻豆av在线| 国产伦理片在线播放av一区| 午夜福利视频精品| 亚洲欧美一区二区三区黑人| 中文亚洲av片在线观看爽 | 美女福利国产在线| 亚洲av日韩在线播放| 亚洲精品成人av观看孕妇| tocl精华| 男女午夜视频在线观看| 久久久精品94久久精品| 久久久水蜜桃国产精品网| 99久久99久久久精品蜜桃| 最近最新中文字幕大全免费视频| 国产亚洲一区二区精品| 极品人妻少妇av视频| 日韩 欧美 亚洲 中文字幕| 欧美大码av| 国产在线免费精品| 亚洲va日本ⅴa欧美va伊人久久| 最近最新免费中文字幕在线| 老汉色∧v一级毛片| 亚洲熟女精品中文字幕| 十八禁高潮呻吟视频| 国产男女内射视频| 老汉色av国产亚洲站长工具| 夜夜骑夜夜射夜夜干| av超薄肉色丝袜交足视频| 亚洲三区欧美一区| 最新的欧美精品一区二区| 欧美精品亚洲一区二区| 久久九九热精品免费| 老熟妇仑乱视频hdxx| 欧美精品一区二区免费开放| 一区二区三区激情视频| 日本vs欧美在线观看视频| 精品一品国产午夜福利视频| 看免费av毛片| 一夜夜www| 亚洲成人免费电影在线观看| 久久久精品国产亚洲av高清涩受| avwww免费| 高潮久久久久久久久久久不卡| 国产精品自产拍在线观看55亚洲 | 国产视频一区二区在线看| 欧美成狂野欧美在线观看| 国产成人欧美| 欧美人与性动交α欧美精品济南到| 老司机靠b影院| 桃红色精品国产亚洲av| 成年版毛片免费区| 纵有疾风起免费观看全集完整版| 老熟妇乱子伦视频在线观看| 成年人黄色毛片网站| 午夜免费成人在线视频| 正在播放国产对白刺激| 国产一区二区三区在线臀色熟女 | 在线观看免费视频网站a站| 午夜免费鲁丝| 俄罗斯特黄特色一大片| 少妇 在线观看| bbb黄色大片| 国产淫语在线视频| 黑人猛操日本美女一级片| kizo精华| 99九九在线精品视频| cao死你这个sao货| 在线亚洲精品国产二区图片欧美| 大片电影免费在线观看免费| 一二三四社区在线视频社区8| av电影中文网址| 久久99一区二区三区| 妹子高潮喷水视频| 老司机午夜福利在线观看视频 | 人妻久久中文字幕网| 男女床上黄色一级片免费看| 我要看黄色一级片免费的| 一区二区三区精品91| 国产日韩欧美亚洲二区| 久久午夜综合久久蜜桃| 99国产综合亚洲精品| 国产精品亚洲一级av第二区| 免费不卡黄色视频| 亚洲第一青青草原| 又黄又粗又硬又大视频| a在线观看视频网站| 国产日韩欧美亚洲二区| 国产一区二区 视频在线| 久久中文字幕一级| 热re99久久精品国产66热6| 夜夜夜夜夜久久久久| 亚洲精品一二三| 亚洲av成人一区二区三| 国产黄色免费在线视频| 国产亚洲精品第一综合不卡| 久久国产精品人妻蜜桃| 免费观看人在逋| 精品国产亚洲在线| 大码成人一级视频| 国产欧美亚洲国产| 成年版毛片免费区| 人人妻人人澡人人看| 久久亚洲精品不卡| 亚洲av成人一区二区三| 超色免费av| 女人高潮潮喷娇喘18禁视频| 国产人伦9x9x在线观看| 日韩欧美免费精品| 久久久久久免费高清国产稀缺| 亚洲成国产人片在线观看| 美女高潮喷水抽搐中文字幕| 一本综合久久免费| 亚洲天堂av无毛| 精品视频人人做人人爽| 成人精品一区二区免费| 99riav亚洲国产免费| 色在线成人网| 看免费av毛片| 日韩人妻精品一区2区三区| 大型黄色视频在线免费观看| av一本久久久久| 国产精品av久久久久免费| 日韩中文字幕欧美一区二区| 一本—道久久a久久精品蜜桃钙片| 国产精品二区激情视频| 亚洲一区中文字幕在线| 一进一出抽搐动态| 日韩欧美一区视频在线观看| 亚洲av成人一区二区三| 国产色视频综合| 免费少妇av软件| 窝窝影院91人妻| 亚洲国产欧美一区二区综合| 99国产极品粉嫩在线观看| 久久热在线av| 亚洲七黄色美女视频| 少妇粗大呻吟视频| 高清在线国产一区| 中文字幕另类日韩欧美亚洲嫩草| 悠悠久久av| 国产精品久久电影中文字幕 | 在线看a的网站| 老司机靠b影院| svipshipincom国产片| 国产精品av久久久久免费| 999精品在线视频| √禁漫天堂资源中文www| 91成年电影在线观看| 丰满饥渴人妻一区二区三| 午夜激情久久久久久久| 成人亚洲精品一区在线观看| 日本黄色视频三级网站网址 | 黄色成人免费大全| 夜夜夜夜夜久久久久| 日本黄色日本黄色录像| 亚洲专区国产一区二区| 亚洲精品乱久久久久久| av有码第一页| 女同久久另类99精品国产91| 一本—道久久a久久精品蜜桃钙片| 久久久水蜜桃国产精品网| 少妇 在线观看| 成年女人毛片免费观看观看9 | 一级片免费观看大全| videos熟女内射| 精品国内亚洲2022精品成人 | 老司机午夜十八禁免费视频| 国产精品二区激情视频| 欧美精品亚洲一区二区| 成人国产av品久久久| 啦啦啦视频在线资源免费观看| 18在线观看网站| 天天影视国产精品| 丝袜美足系列| 男人操女人黄网站| 在线永久观看黄色视频| 1024视频免费在线观看| 涩涩av久久男人的天堂| av国产精品久久久久影院| 日韩大码丰满熟妇|