• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于Dirichlet-to-Neumann映射計算旋磁光子晶體的帶隙結(jié)構(gòu)*

    2020-05-29 11:56:44楊華嬌胡真
    關(guān)鍵詞:帶隙晶格光子

    楊華嬌, 胡真

    (河海大學(xué) 理學(xué)院,江蘇 南京 211100)

    光子晶體(Photonic Crystals)是一種新型的周期性人造材料,周期為光波長量級[1].光子晶體最重要的性質(zhì)之一是具有光子帶隙(Photonic Band Gap),頻率落在光子帶隙范圍內(nèi)的光波將無法在光子晶體內(nèi)部進行傳播,從而光子晶體具有了控制光的傳播的能力,被制成各種光子晶體元件應(yīng)用于各個領(lǐng)域,如無損耗光波導(dǎo)[2]、高輻射亮度超微激光器[3]、高光電轉(zhuǎn)化效率薄膜光伏電池[4]等.

    在過去的二十年中,對于二維光子晶體的研究主要集中在各向同性材料制成的光子晶體.近年來,人們發(fā)現(xiàn)基于旋磁各向異性材料制作的二維旋磁光子晶體具有更多特殊的性質(zhì).旋磁各向異性材料是指該材料的相對磁導(dǎo)率會隨著方向的變化而改變.利用二維旋磁光子晶體,Wang等設(shè)計并制作出了只能單向傳播的波導(dǎo)管[5-6];He等設(shè)計了一種單向交叉波導(dǎo)分路器[7];Chen等設(shè)計了一種無色散的慢光波導(dǎo)管[8].

    分析設(shè)計二維旋磁光子晶體最基本的問題是計算旋磁光子晶體的帶隙結(jié)構(gòu).計算帶隙結(jié)構(gòu)需要有效的數(shù)值方法,傳統(tǒng)的數(shù)值方法有:平面波展開法[9]、有限元法[10]、多重散射法[11]等等.但是這些方法各有不足:平面波展開法和有限元方法對于非色散介質(zhì)而言,建立的是線性特征值問題,但對于色散介質(zhì)建立的卻是非線性特征值問題;多重散射法需要格點求和技術(shù),計算過程復(fù)雜.

    近年來Lu等[15]提出了一種能有效分析各種光子晶體結(jié)構(gòu)的新的數(shù)值方法,叫作Dirichlet-to-Neumann(DtN)映射方法.單元晶格的DtN映射將單元晶格邊界上的波動場映射成邊界上的法向?qū)?shù),通過DtN映射避免了在晶格內(nèi)部的離散,從而達(dá)到簡化運算的效果.DtN映射方法已經(jīng)被廣泛用于分析二維各向同性光子晶體結(jié)構(gòu),

    本文將DtN映射方法應(yīng)用于計算二維旋磁光子晶體的帶隙結(jié)構(gòu).旋磁光子晶體的相對磁導(dǎo)率μ為一個矩陣,控制方程及不同介質(zhì)邊界面上的連續(xù)性條件均比較復(fù)雜,因此需要重新在旋磁光子晶體的單元晶格上構(gòu)造出相應(yīng)的DtN映射,再基于該DtN映射建立特征值問題來求解旋磁光子晶體的帶隙結(jié)構(gòu).利用擴展的DtN映射建立的特征值問題以波數(shù)β的函數(shù)為特征值,以頻率ω為給定的參數(shù).通過旋磁光子晶體單元晶格的DtN映射,得到較小矩陣的特征值問題,且無論是色散介質(zhì)還是非色散介質(zhì),該方法建立的特征值問題均為線性的.

    1 基本方程

    旋磁光子晶體的相對磁導(dǎo)率μ為3×3維的對稱矩陣

    (1)

    本文僅考慮二維問題,且只研究TE模式,利用Maxwell方程組,可得TE模式下的控制方程為

    (2)

    其中u為電場在z方向上的分量,是二維梯度算子,

    2 二維旋磁光子晶體單元晶格的DtN映射

    為了建立特征值問題來求解二維旋磁光子晶體的帶隙結(jié)構(gòu),需要首先構(gòu)造出其單元晶格的DtN映射.單元晶格的DtN映射Λ將單元晶格邊界Γ上的波動場u映射成u在Γ上的法向?qū)?shù),即

    Λu|Γ=?vu|Γ

    (3)

    其中v表示邊界Γ上的單位法向量.旋磁光子晶體是由旋磁各向異性材料在空間中周期排列而構(gòu)成的人造材料,常用的周期排列方式有正方形、三角形和蜂巢狀三種(如圖1所示);圖中線段框出的部分為對應(yīng)結(jié)構(gòu)的單元晶格,a為晶格常數(shù),其中正方形和三角形周期排列的磁光子晶體的單元晶格都是只包含一個圓柱的晶格,他們的DtN映射可以直接構(gòu)造;蜂巢狀周期排列的旋磁光子晶體的單元晶格由三個小正六邊形晶格組成(兩個包含圓柱,一個不包含圓柱),其DtN映射可以利用這三個小晶格的DtN映射構(gòu)造得到.

    (a)正方形周期排列介質(zhì)柱型 (b)三角形周期排列薄板型 (c)蜂巢狀周期排列介質(zhì)柱型

    圖1 旋磁光子晶體

    Fig.1 Gyromagnetic photonic crystals

    2.1 只包含一個圓柱的晶格的DtN映射

    控制方程(4)在只含有一個圓柱的結(jié)構(gòu)上的通解u(x,y)可以用特解φj(r)eijθ的線性組合來表示,即

    (4)

    其中,cj為未知系數(shù),(r,θ)是點(x,y)在極坐標(biāo)下的坐標(biāo),φj(r)與Bessel函數(shù)Jj、Yj有關(guān),表示為

    (5)

    Jj(q)Aj+Yj(q)Bj=Jj(s)

    (6)

    其中,s=k0n1a,q=k0n2a.

    (7)

    為了構(gòu)造晶格的DtN映射,還需要將通解u(x,y)的法向?qū)?shù)用特解法向?qū)?shù)的線性組合表示出來:

    (8)

    由(4)式和(8)式消去未知系數(shù)cj,從而得到只包含一個圓柱的晶格的DtN映射Λ.例如,假設(shè)在邊界Γ上共離散p個點(p為偶數(shù)),離散點的坐標(biāo)為(xk,yk),(rk,θk)為該點的極坐標(biāo),其中k=1,2,…,p,則由(4)式得到

    (9)

    再由(8)式得到

    (10)

    將(9)和(10)改寫為矩陣形式

    (11)

    Λ=BAV-1

    若在正邊形晶格的每條邊上離散N個點,則p=4N,Λ近似為(4N)×(4N)的矩陣;若在六邊形晶格的每條邊上離散N個點,則p=6N,Λ近似為(6N)×(6N)的矩陣.

    2.2 不包含圓柱的晶格的DtN映射

    對不包含圓柱的晶格而言,其DtN映射的構(gòu)造過程和包含圓柱的晶格的DtN映射構(gòu)造過程類似,唯一的區(qū)別在于方程(4)的特解更加簡單,為

    φj(r)=Jj(k0nr)

    2.3 蜂巢狀單元晶格的DtN映射

    蜂巢狀單元晶格由兩個只包含一個圓柱的小晶格和一個不包含圓柱的小晶格組成.為了構(gòu)造出蜂巢狀單元晶格的DtN映射,需要先建立Ω2、Ω3(如圖2(c)所示)這兩個小晶格的DtN映射Λ(2)、Λ(3),再利用2.2節(jié)中介紹的方法構(gòu)造不包含圓柱的晶格Ω1的映射Λ(1).利用這三個小晶格的DtN映射可以在公共邊上建立如下方程:

    A1V=A2U

    (12)

    其中,U為蜂巢狀單元晶邊界上的波動場,V為三條公共邊上的波動場.

    蜂巢狀單元晶格邊界上波動場的法向?qū)?shù)可以由對應(yīng)小晶格的DtN映射得到

    ?vU=B1U+B2V

    (13)

    通過方程(12)和(13)消去公共邊上的波動場V,從而得到蜂巢狀單元晶格的DtN映射

    3 特征值問題

    構(gòu)造了二維旋磁光子晶體單元晶格的DtN映射之后,便可以建立求解旋磁光子晶體帶隙結(jié)構(gòu)的特征值問題.這一建立過程與求解各向同性光子晶體帶隙結(jié)構(gòu)的建立過程類似,其中求解正方形、三角形和蜂巢狀周期排列的各向同性光子晶體帶隙結(jié)構(gòu)的特征值問題分別在Yuan[16]、Yuan[17]和李靜[18]的文章中有詳細(xì)介紹.下文簡要介紹這一過程.

    分析旋磁光子晶體的帶隙結(jié)構(gòu),需要求解方程(4)在相應(yīng)結(jié)構(gòu)中的Bloch模式解:

    u(x,y)=ei(αx+βy)Ψ(x,y)

    其中,(α,β)為Bloch波矢,Ψ為周期函數(shù).

    將在第2節(jié)中得到的單元晶格的DtN映射Λ按照單元晶格的邊數(shù)m進行分塊,使其每一個分塊均為N×N的矩陣,且記單元晶格每條邊上的波動場分別為u1,u2,…,um(正方形m=4;六邊形m=6,蜂巢狀m=12),從而可將(3)式改寫為

    (14)

    再分別利用正方形、三角形和蜂巢狀周期排列的旋磁光子晶體的擬周期邊界條件和(14)式,便可以得到求解他們帶隙結(jié)構(gòu)的特征值問題.由于這三種結(jié)構(gòu)建立特征值問題的過程類似,只以正方形周期排列的旋磁光子晶體結(jié)構(gòu)為例說明此過程.

    圖2 正方形周期排列旋磁光子晶體的布里淵區(qū)Fig.2 Brillouin zone of a square photonic crystals

    正方形單元晶格的擬周期邊界條件為:

    u3=ρβu1,?vu3=ρβ?vu1;u4=ραu2,?vu4=ρα?vu2

    (15)

    其中,ρα=eiαL,ρβ=eiβL.

    將(13)式代入(12)式,便可得到求解正方形周期排列旋磁光子晶體帶隙結(jié)構(gòu)的特征值問題:

    其中,I為單位矩陣,V=[λUT,UT]T,U=[u1,u2]T,特征值λ和矩陣M0、M1、M2的取值與正方形周期排列旋磁光子晶體的布里淵區(qū)有關(guān).

    在ΓX、XM、MΓ邊上特征值λ與矩陣M0、M1、M2的取值分別如下所示:

    (1)從Γ到X,λ=eiαa(0<α≤π/a)且

    (2)從X到M,λ=eiβa(0<β≤π/a)且

    (3)從M到Γ,λ=eiαa=eiβa(0<α=β≤π/a)且

    建立求解三角形和蜂巢狀周期排列的旋磁光子晶體帶隙結(jié)構(gòu)的特征值問題的過程與此類似,差別僅在于擬周期邊界條件以及相對應(yīng)的布里淵區(qū)略有不同.

    4 數(shù)值算例

    分別以正方形、三角形和蜂巢狀周期排列的二維旋磁光子晶體為例驗證方法的有效性,并且基于該方法優(yōu)化設(shè)計出了帶隙更寬的蜂巢狀周期排列的旋磁光子晶體.

    4.1 正方形周期排列的介質(zhì)柱型旋磁光子晶體

    首先計算二維正方形周期排列的介質(zhì)柱型旋磁光子晶體的帶隙結(jié)構(gòu).Wang等用有限元法也計算過這一旋磁光子晶體的帶隙結(jié)構(gòu),旋磁各向異性材料采用的是YIG材料[5].YIG的相對介電常數(shù)為εr=15,介質(zhì)柱半徑為r=0.11a,a為晶格常數(shù).當(dāng)外加磁場的作用下引發(fā)YIG強烈的各向異性,此時μr=14,μk=-12.4.計算結(jié)果如圖3所示,陰影部分表示光子帶隙,標(biāo)準(zhǔn)化頻率的帶隙范圍為[0.324 2,0.448]、[0.527,0.576]和[0.611,0.647 5],結(jié)果與Wang等一致.

    圖3 計算結(jié)果 圖4 計算結(jié)果

    4.2 三角形周期排列的薄板型旋磁光子晶體

    接下來,計算二維三角形周期排列的薄板型旋磁光子晶體的帶隙結(jié)構(gòu).Lu等用有限元法和平面波展開法計算過這一旋磁光子晶體的帶隙結(jié)構(gòu),采用的旋磁材料也為YIG材料,材料參數(shù)與4.1中相同,且空氣洞半徑為r=0.44a[14].同樣,在外加磁場的作用下光子晶體磁化.計算結(jié)果如圖4所示,陰影部分表示光子帶隙,標(biāo)準(zhǔn)化頻率的帶隙范圍為[0.341 5,0.361 1],與Lu等一致.

    4.3 蜂巢狀周期排列的介質(zhì)柱型旋磁光子晶體

    最后,計算二維蜂巢狀周期排列的介質(zhì)柱型旋磁光子晶體的帶隙結(jié)構(gòu).Ao等用多重散射法計算過這一旋磁光子晶體的帶隙結(jié)構(gòu),采用的旋磁各向異性材料為鐵氧體,相對介電常數(shù)εm=15,介質(zhì)柱半徑r=0.2a(其中a=10 mm)[13].鐵氧體光子晶體被外加的恒定磁場沿+z方向飽和磁化后便會呈現(xiàn)出各向異性的性質(zhì),其中H0=500 Oe,H0是外加磁場,4πMs=1 750 G(4πMs是飽和磁化強度

    ω0=γH0是進動頻率(旋磁比γ=2.8×106rad·s-1·G-1).計算結(jié)果如圖5所示,第三條能帶和第四條能帶間的標(biāo)準(zhǔn)化的帶隙范圍為[0.252 9,0.272 5],與Ao等一致.

    由于蜂巢狀旋磁光子晶體的單元晶格由三個小正六邊形晶格構(gòu)成,所以有更多的結(jié)構(gòu)參數(shù)可供調(diào)整,特別適合進行優(yōu)化設(shè)計.下面嘗試對圖1(c)所示的蜂巢狀旋磁光子晶體結(jié)構(gòu)進行優(yōu)化,僅通過調(diào)整小晶格Ω2、Ω3中圓柱的半徑大小,來擴大其帶隙范圍.當(dāng)r2=r3=0.27 a時,得到了較寬的帶隙結(jié)構(gòu),計算結(jié)果如圖6所示,第二條能帶和第三條能帶間的帶隙范圍擴大到了[0.221 3,0.257 5].

    圖5 計算結(jié)果 圖6計算結(jié)果

    5 結(jié) 語

    擴展了DtN映射方法,建立了二維旋磁光子晶體單元晶格的DtN映射,并將其用于計算正方形、三角形及蜂巢狀周期排列的旋磁光子晶體的帶隙結(jié)構(gòu),方法避免了在單元晶格內(nèi)部的離散,只需要在邊界上進行離散,從而得到矩陣較小的特征值問題,且計算快速并易于調(diào)整參數(shù),非常有利于優(yōu)化設(shè)計出具有更寬帶隙的二維旋磁光子晶體.目前考慮的旋磁光子晶體都是由圓柱周期排列構(gòu)成的,柱體的截面均為圓形,所以可以使用柱面波展開法來構(gòu)造DtN映射,接下來將進一步研究非圓柱構(gòu)成的旋磁光子晶體的帶隙結(jié)構(gòu),如柱體的截面為三角形、橢圓等其他形狀,這時需要重新設(shè)法來構(gòu)造單元晶格的DtN映射.此外,還可以研究三維旋磁光子晶體的帶隙結(jié)構(gòu),比如在有限厚度的板上周期性的挖取空氣洞形成的結(jié)構(gòu)等.

    猜你喜歡
    帶隙晶格光子
    《光子學(xué)報》征稿簡則
    密度泛函理論計算半導(dǎo)體材料的帶隙誤差研究
    非線性光學(xué)晶格中的梯度流方法
    一種基于BJT工藝的無運放低溫度系數(shù)的帶隙基準(zhǔn)源
    間距比對雙振子局域共振軸縱振帶隙的影響
    一款高PSRR低溫度系數(shù)的帶隙基準(zhǔn)電壓源的設(shè)計
    電子制作(2018年1期)2018-04-04 01:48:38
    一個新非線性可積晶格族和它們的可積辛映射
    一族拉克斯可積晶格方程
    在光子帶隙中原子的自發(fā)衰減
    三維復(fù)式晶格的熱容研究
    亚洲欧洲精品一区二区精品久久久| 色尼玛亚洲综合影院| 久久久久免费精品人妻一区二区| 久久精品国产亚洲av香蕉五月| 国产高清三级在线| 国产精品久久久人人做人人爽| 亚洲成人久久爱视频| 丁香欧美五月| 午夜久久久久精精品| 国产真人三级小视频在线观看| 一边摸一边抽搐一进一小说| 国产午夜福利久久久久久| 日本在线视频免费播放| 免费大片18禁| 午夜福利在线观看免费完整高清在 | 亚洲国产精品sss在线观看| 老司机深夜福利视频在线观看| 日本一本二区三区精品| 中文字幕精品亚洲无线码一区| 91av网站免费观看| 在线永久观看黄色视频| 日本黄大片高清| 久久久久免费精品人妻一区二区| 一本精品99久久精品77| 午夜福利在线观看吧| 成年女人永久免费观看视频| 国产不卡一卡二| 午夜视频精品福利| 无限看片的www在线观看| 久久人妻av系列| 噜噜噜噜噜久久久久久91| 国产日本99.免费观看| 岛国视频午夜一区免费看| 午夜影院日韩av| 国产毛片a区久久久久| 国产成人精品久久二区二区免费| 狂野欧美激情性xxxx| 欧美色视频一区免费| 国产精品九九99| 免费看a级黄色片| 天天一区二区日本电影三级| 久久热在线av| 久久这里只有精品19| 性色avwww在线观看| 亚洲熟妇中文字幕五十中出| 少妇丰满av| 人人妻,人人澡人人爽秒播| 97碰自拍视频| 一区二区三区激情视频| 岛国在线观看网站| 欧美一级a爱片免费观看看| 国产精品永久免费网站| 老司机午夜福利在线观看视频| 国产淫片久久久久久久久 | 99re在线观看精品视频| h日本视频在线播放| 99久久精品一区二区三区| 日韩欧美在线二视频| 亚洲色图 男人天堂 中文字幕| 欧美一区二区国产精品久久精品| 巨乳人妻的诱惑在线观看| 日本成人三级电影网站| 怎么达到女性高潮| 精品日产1卡2卡| 禁无遮挡网站| 欧美性猛交黑人性爽| 草草在线视频免费看| 午夜视频精品福利| 免费看光身美女| 国产精品电影一区二区三区| 久久久久性生活片| 欧美zozozo另类| 亚洲av成人av| 亚洲成人中文字幕在线播放| 国产黄片美女视频| 免费一级毛片在线播放高清视频| 成人鲁丝片一二三区免费| 亚洲精华国产精华精| 亚洲中文av在线| 日本 av在线| 日韩人妻高清精品专区| 亚洲国产精品合色在线| 欧美一级毛片孕妇| 十八禁人妻一区二区| 精品国产美女av久久久久小说| 亚洲精品在线美女| 亚洲欧美日韩高清在线视频| bbb黄色大片| 91在线精品国自产拍蜜月 | 亚洲 欧美 日韩 在线 免费| 国产精品99久久久久久久久| 亚洲人成网站高清观看| 久久草成人影院| 狂野欧美白嫩少妇大欣赏| 老鸭窝网址在线观看| 一个人免费在线观看的高清视频| 亚洲国产看品久久| 亚洲一区二区三区色噜噜| 国产精品乱码一区二三区的特点| 久久热在线av| 国产v大片淫在线免费观看| 男插女下体视频免费在线播放| 少妇裸体淫交视频免费看高清| 久久中文看片网| 美女大奶头视频| 日韩人妻高清精品专区| 国产午夜福利久久久久久| 国产黄片美女视频| 久久精品亚洲精品国产色婷小说| 国产高清三级在线| 精品久久久久久久久久久久久| 99热6这里只有精品| 两性夫妻黄色片| 女警被强在线播放| 日韩av在线大香蕉| 99在线人妻在线中文字幕| 91麻豆精品激情在线观看国产| 巨乳人妻的诱惑在线观看| 亚洲欧美日韩无卡精品| 97碰自拍视频| 久久久久久九九精品二区国产| 国产极品精品免费视频能看的| 曰老女人黄片| 久久国产精品影院| 欧美黄色淫秽网站| 桃红色精品国产亚洲av| 国产伦人伦偷精品视频| 国产一区二区在线av高清观看| 国产伦精品一区二区三区视频9 | 国产成人啪精品午夜网站| 亚洲一区高清亚洲精品| 在线免费观看不下载黄p国产 | 日本免费a在线| 久久精品综合一区二区三区| 久久精品影院6| 国产精品,欧美在线| 黄色女人牲交| 99re在线观看精品视频| 18禁观看日本| 日本黄大片高清| 老司机在亚洲福利影院| 搡老熟女国产l中国老女人| 免费观看精品视频网站| 一级毛片女人18水好多| 成年女人毛片免费观看观看9| 99精品久久久久人妻精品| 亚洲欧美日韩卡通动漫| 香蕉丝袜av| 99国产极品粉嫩在线观看| www日本在线高清视频| 国内精品久久久久久久电影| 日本 欧美在线| 叶爱在线成人免费视频播放| 成人特级av手机在线观看| 国产精品久久久久久精品电影| 日韩欧美精品v在线| www.自偷自拍.com| 亚洲精品在线美女| 中文亚洲av片在线观看爽| 老司机在亚洲福利影院| 在线国产一区二区在线| 亚洲天堂国产精品一区在线| 成人亚洲精品av一区二区| 精品久久蜜臀av无| 精品无人区乱码1区二区| 午夜免费激情av| 一边摸一边抽搐一进一小说| 国产成人精品久久二区二区91| 欧美性猛交黑人性爽| 久久久精品欧美日韩精品| 亚洲国产欧美人成| 亚洲精品粉嫩美女一区| 美女黄网站色视频| 久久中文看片网| 欧美午夜高清在线| 男女床上黄色一级片免费看| 国产高清激情床上av| 97人妻精品一区二区三区麻豆| 男女午夜视频在线观看| 国产av在哪里看| 又粗又爽又猛毛片免费看| 色哟哟哟哟哟哟| 国产午夜福利久久久久久| 免费av不卡在线播放| 国内精品一区二区在线观看| 又黄又粗又硬又大视频| 无限看片的www在线观看| 一本综合久久免费| 国产午夜精品论理片| 成人精品一区二区免费| netflix在线观看网站| 一级毛片精品| 国产1区2区3区精品| 成人欧美大片| 熟女少妇亚洲综合色aaa.| 成年免费大片在线观看| 欧美丝袜亚洲另类 | 亚洲aⅴ乱码一区二区在线播放| 最近视频中文字幕2019在线8| 久久中文字幕人妻熟女| 午夜亚洲福利在线播放| 免费大片18禁| 天堂网av新在线| 久久久久精品国产欧美久久久| 国产成人啪精品午夜网站| 国内久久婷婷六月综合欲色啪| 在线播放国产精品三级| 一本精品99久久精品77| 国产在线精品亚洲第一网站| 脱女人内裤的视频| 日日夜夜操网爽| 看黄色毛片网站| 岛国在线观看网站| 一级毛片高清免费大全| 天堂网av新在线| 精品国内亚洲2022精品成人| 99久久精品一区二区三区| 国产熟女xx| 国产精品国产高清国产av| 国产乱人伦免费视频| 在线观看免费视频日本深夜| 全区人妻精品视频| 每晚都被弄得嗷嗷叫到高潮| 日韩精品青青久久久久久| 久久人人精品亚洲av| 亚洲国产看品久久| 午夜免费观看网址| 九色成人免费人妻av| 露出奶头的视频| 亚洲七黄色美女视频| 在线观看午夜福利视频| 18禁裸乳无遮挡免费网站照片| 午夜精品久久久久久毛片777| 日韩有码中文字幕| 国产麻豆成人av免费视频| 非洲黑人性xxxx精品又粗又长| 国产成人欧美在线观看| 日本在线视频免费播放| 人妻丰满熟妇av一区二区三区| 一级毛片高清免费大全| 无限看片的www在线观看| 国产毛片a区久久久久| 亚洲精品在线美女| 波多野结衣高清无吗| 国产又色又爽无遮挡免费看| 国产高清激情床上av| 日本精品一区二区三区蜜桃| 久久久久久国产a免费观看| 90打野战视频偷拍视频| 国产伦人伦偷精品视频| 国产黄a三级三级三级人| 欧美日韩中文字幕国产精品一区二区三区| 99热6这里只有精品| 成人亚洲精品av一区二区| 精品久久久久久久久久免费视频| 成人精品一区二区免费| 国产精品一区二区三区四区免费观看 | 真人做人爱边吃奶动态| 婷婷精品国产亚洲av| 免费在线观看影片大全网站| 精品久久久久久,| 久久中文看片网| 偷拍熟女少妇极品色| a级毛片在线看网站| 国产真实乱freesex| 99久久精品国产亚洲精品| 国产精品久久电影中文字幕| 亚洲av成人精品一区久久| 黑人操中国人逼视频| 99国产精品一区二区三区| 999久久久精品免费观看国产| 亚洲精品一区av在线观看| 亚洲美女视频黄频| 99在线人妻在线中文字幕| 两个人视频免费观看高清| 午夜精品在线福利| 亚洲精品456在线播放app | 一卡2卡三卡四卡精品乱码亚洲| 国产高清三级在线| 欧美日韩国产亚洲二区| 亚洲中文字幕日韩| 亚洲成人中文字幕在线播放| 亚洲国产欧洲综合997久久,| 狠狠狠狠99中文字幕| 热99re8久久精品国产| 精品久久久久久久久久久久久| 亚洲成av人片在线播放无| 三级男女做爰猛烈吃奶摸视频| 婷婷精品国产亚洲av在线| 久久中文看片网| 国产黄色小视频在线观看| 亚洲成av人片免费观看| 99久国产av精品| 欧美一级a爱片免费观看看| 免费看光身美女| 国产高清视频在线播放一区| 啦啦啦免费观看视频1| 床上黄色一级片| 国产精华一区二区三区| 美女高潮的动态| 国产精品一区二区免费欧美| 老司机深夜福利视频在线观看| 黄色 视频免费看| 国产精品av视频在线免费观看| 少妇的丰满在线观看| 99久国产av精品| 啦啦啦韩国在线观看视频| 黄色日韩在线| 午夜影院日韩av| 一二三四社区在线视频社区8| 亚洲aⅴ乱码一区二区在线播放| 国产麻豆成人av免费视频| www.熟女人妻精品国产| netflix在线观看网站| 全区人妻精品视频| 深夜精品福利| 好男人在线观看高清免费视频| 夜夜躁狠狠躁天天躁| 日本熟妇午夜| 日本五十路高清| 长腿黑丝高跟| 国产精品久久视频播放| 成人鲁丝片一二三区免费| 精品国内亚洲2022精品成人| 国产精品99久久99久久久不卡| 国产男靠女视频免费网站| 天堂√8在线中文| 日韩高清综合在线| 老司机福利观看| 欧美日本视频| 我的老师免费观看完整版| 国产蜜桃级精品一区二区三区| 午夜亚洲福利在线播放| 国产美女午夜福利| xxx96com| 国产精品98久久久久久宅男小说| 久久人人精品亚洲av| 国产精品爽爽va在线观看网站| a在线观看视频网站| 国产激情久久老熟女| 成人永久免费在线观看视频| 日本免费一区二区三区高清不卡| 国产91精品成人一区二区三区| av视频在线观看入口| 午夜激情福利司机影院| 人妻丰满熟妇av一区二区三区| av天堂中文字幕网| 美女午夜性视频免费| 叶爱在线成人免费视频播放| 亚洲欧美日韩高清专用| 亚洲天堂国产精品一区在线| 成人国产综合亚洲| 欧美日本亚洲视频在线播放| 国产精品日韩av在线免费观看| 禁无遮挡网站| 天堂影院成人在线观看| h日本视频在线播放| 美女午夜性视频免费| 少妇裸体淫交视频免费看高清| 美女 人体艺术 gogo| 国产激情欧美一区二区| 国产午夜福利久久久久久| 国产三级中文精品| av欧美777| 日韩欧美免费精品| 日本撒尿小便嘘嘘汇集6| 国产69精品久久久久777片 | av福利片在线观看| 麻豆久久精品国产亚洲av| 精品乱码久久久久久99久播| 一个人观看的视频www高清免费观看 | a级毛片在线看网站| 精品午夜福利视频在线观看一区| 欧美乱色亚洲激情| 2021天堂中文幕一二区在线观| 成人特级av手机在线观看| 欧美日韩一级在线毛片| 怎么达到女性高潮| 人妻夜夜爽99麻豆av| 国产高清视频在线观看网站| www.www免费av| 高清在线国产一区| 久久国产精品影院| 男女做爰动态图高潮gif福利片| cao死你这个sao货| 一卡2卡三卡四卡精品乱码亚洲| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲乱码一区二区免费版| 免费大片18禁| 国产97色在线日韩免费| 国产亚洲精品av在线| 精品人妻1区二区| 亚洲狠狠婷婷综合久久图片| 两人在一起打扑克的视频| 亚洲av片天天在线观看| 神马国产精品三级电影在线观看| 亚洲精品国产精品久久久不卡| 亚洲欧美日韩高清专用| 真实男女啪啪啪动态图| 久久精品91蜜桃| 岛国视频午夜一区免费看| 麻豆国产av国片精品| 中文字幕久久专区| ponron亚洲| 欧美一区二区国产精品久久精品| 亚洲人与动物交配视频| 天堂av国产一区二区熟女人妻| 午夜视频精品福利| 天天添夜夜摸| 成人鲁丝片一二三区免费| 一个人看视频在线观看www免费 | 国产一级毛片七仙女欲春2| 亚洲,欧美精品.| 香蕉国产在线看| 美女被艹到高潮喷水动态| www.自偷自拍.com| 中文字幕熟女人妻在线| 欧美黄色淫秽网站| 少妇熟女aⅴ在线视频| 国产又色又爽无遮挡免费看| 啪啪无遮挡十八禁网站| 高清在线国产一区| 国产精品久久久久久久电影 | www日本黄色视频网| 国产高清有码在线观看视频| 午夜福利在线在线| 午夜亚洲福利在线播放| 欧美另类亚洲清纯唯美| 亚洲精品一卡2卡三卡4卡5卡| 18禁国产床啪视频网站| 欧美成狂野欧美在线观看| 丁香六月欧美| 日韩欧美国产一区二区入口| 国产精品女同一区二区软件 | tocl精华| 少妇的逼水好多| 日韩欧美一区二区三区在线观看| 精品无人区乱码1区二区| 亚洲精品一卡2卡三卡4卡5卡| 18禁国产床啪视频网站| 别揉我奶头~嗯~啊~动态视频| 国产激情偷乱视频一区二区| 亚洲真实伦在线观看| av视频在线观看入口| 精品久久久久久久久久久久久| 精品一区二区三区视频在线观看免费| 国产成人福利小说| 亚洲成a人片在线一区二区| 天天一区二区日本电影三级| 草草在线视频免费看| 男人和女人高潮做爰伦理| 人人妻人人看人人澡| 中文字幕人成人乱码亚洲影| 这个男人来自地球电影免费观看| 91老司机精品| 99国产极品粉嫩在线观看| 岛国在线免费视频观看| 日韩欧美国产在线观看| 日韩免费av在线播放| 观看美女的网站| 国产精品99久久久久久久久| 天天一区二区日本电影三级| 香蕉av资源在线| 亚洲人与动物交配视频| 欧美绝顶高潮抽搐喷水| 一进一出好大好爽视频| 亚洲成人久久爱视频| 欧美精品啪啪一区二区三区| 日韩欧美三级三区| 精品无人区乱码1区二区| 色老头精品视频在线观看| 99精品欧美一区二区三区四区| 午夜激情欧美在线| 国产精品一区二区三区四区免费观看 | 亚洲熟妇中文字幕五十中出| 老司机午夜十八禁免费视频| 曰老女人黄片| 天堂√8在线中文| 少妇的丰满在线观看| 中文亚洲av片在线观看爽| 综合色av麻豆| 久久伊人香网站| 日韩人妻高清精品专区| 国产成人精品无人区| 男女那种视频在线观看| 日本与韩国留学比较| 亚洲精华国产精华精| 免费一级毛片在线播放高清视频| 美女免费视频网站| 午夜精品在线福利| 特级一级黄色大片| 88av欧美| 亚洲第一欧美日韩一区二区三区| 欧美三级亚洲精品| 欧美日韩乱码在线| 男女那种视频在线观看| 色在线成人网| 给我免费播放毛片高清在线观看| 国产伦在线观看视频一区| 黄片小视频在线播放| 日本撒尿小便嘘嘘汇集6| 在线观看免费视频日本深夜| 噜噜噜噜噜久久久久久91| 精品久久久久久久久久久久久| 97碰自拍视频| 亚洲精品久久国产高清桃花| 亚洲国产日韩欧美精品在线观看 | 村上凉子中文字幕在线| 在线观看午夜福利视频| 亚洲无线观看免费| 亚洲av美国av| 99久国产av精品| 999久久久国产精品视频| 夜夜爽天天搞| 在线国产一区二区在线| 久久天堂一区二区三区四区| 夜夜夜夜夜久久久久| 好男人电影高清在线观看| 极品教师在线免费播放| 99热这里只有是精品50| 精华霜和精华液先用哪个| 一级作爱视频免费观看| 老司机在亚洲福利影院| 亚洲av电影在线进入| 男女那种视频在线观看| 亚洲成人中文字幕在线播放| 一级作爱视频免费观看| 黄色日韩在线| 欧美日韩乱码在线| 美女被艹到高潮喷水动态| 黄色女人牲交| 久久精品国产99精品国产亚洲性色| 老汉色∧v一级毛片| 亚洲国产欧美一区二区综合| 韩国av一区二区三区四区| 长腿黑丝高跟| 亚洲第一电影网av| 无遮挡黄片免费观看| 午夜精品久久久久久毛片777| 国语自产精品视频在线第100页| 中文字幕熟女人妻在线| 精品一区二区三区av网在线观看| avwww免费| 成人鲁丝片一二三区免费| 午夜视频精品福利| 国模一区二区三区四区视频 | 久久久久久久久久黄片| 欧美日韩亚洲国产一区二区在线观看| 午夜福利在线在线| 国产欧美日韩一区二区精品| av福利片在线观看| 搡老妇女老女人老熟妇| 真人做人爱边吃奶动态| 天堂动漫精品| 欧美大码av| 久久精品夜夜夜夜夜久久蜜豆| 国产亚洲欧美98| 日本免费一区二区三区高清不卡| 桃红色精品国产亚洲av| 网址你懂的国产日韩在线| 在线观看美女被高潮喷水网站 | 午夜两性在线视频| 俺也久久电影网| 日韩中文字幕欧美一区二区| 一个人观看的视频www高清免费观看 | svipshipincom国产片| 亚洲国产精品999在线| 很黄的视频免费| 日本免费一区二区三区高清不卡| 免费观看的影片在线观看| 又紧又爽又黄一区二区| 91久久精品国产一区二区成人 | 午夜精品久久久久久毛片777| 久久香蕉精品热| 久久天躁狠狠躁夜夜2o2o| 岛国在线观看网站| 国内毛片毛片毛片毛片毛片| 亚洲中文日韩欧美视频| 久久香蕉国产精品| 怎么达到女性高潮| 深夜精品福利| 国产aⅴ精品一区二区三区波| 草草在线视频免费看| 日韩欧美在线二视频| 色噜噜av男人的天堂激情| 国产亚洲欧美98| 婷婷精品国产亚洲av在线| 美女高潮的动态| 日韩欧美国产一区二区入口| 精品熟女少妇八av免费久了| 精品国产超薄肉色丝袜足j| 俺也久久电影网| 国产视频一区二区在线看| 久久精品91无色码中文字幕| 日韩精品中文字幕看吧| 亚洲精品中文字幕一二三四区| 香蕉久久夜色| 叶爱在线成人免费视频播放| 欧美日本视频| 成在线人永久免费视频| 欧美高清成人免费视频www| 久久久国产精品麻豆| 一区二区三区国产精品乱码| 久久久久久国产a免费观看| 一级毛片精品| 国产亚洲欧美在线一区二区| 国产精品野战在线观看| 国产三级黄色录像| 看片在线看免费视频| 久久久久久久午夜电影| 久久久久久久久中文| 露出奶头的视频| 国产亚洲精品av在线| 亚洲九九香蕉| 两性午夜刺激爽爽歪歪视频在线观看|