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

    二維蜂窩狀光子晶體能帶的數(shù)學(xué)模型與仿真研究

    2023-04-04 08:12:16
    科技創(chuàng)新與應(yīng)用 2023年9期
    關(guān)鍵詞:蜂窩狀帶隙能帶

    楊 之

    (北京郵電大學(xué),北京 100876)

    光子晶體是由宏觀光學(xué)介質(zhì)代替微觀粒子按一定規(guī)律周期性排列的人工微觀結(jié)構(gòu),其周期陣列限制了電磁波在光子晶體中的傳播,對特定波長的光子具有阻擋作用,從而形成了光子帶隙[1]。根據(jù)應(yīng)用需要設(shè)計光子晶體的結(jié)構(gòu),由于光子帶隙的存在可以禁止光在一定頻率范圍內(nèi)的傳播,使得光子晶體有著廣泛的應(yīng)用,其應(yīng)用包括窄帶濾波器、分復(fù)用器等[2-5]。獲得寬的光子晶體帶隙是設(shè)計高性能光子晶體器件的基礎(chǔ)[6],所以如何設(shè)計光子晶體結(jié)構(gòu)來獲得更寬的帶隙一直是光子晶體領(lǐng)域的研究重點(diǎn)。

    已有研究結(jié)果表明,二維蜂窩狀光子晶體具有較大的帶隙[7],且?guī)洞笮『徒橘|(zhì)柱的半徑大小與介電常數(shù)密切相關(guān)[8-9]。改變光子晶體的晶格結(jié)構(gòu)、降低結(jié)構(gòu)對稱性和引入新材料等措施都可以增大光子晶體的帶隙,但是由于直接進(jìn)行實驗來設(shè)計新型光子晶體結(jié)構(gòu)的成本巨大,現(xiàn)在一般先采用數(shù)值模擬方法進(jìn)行結(jié)構(gòu)設(shè)計,以期獲得結(jié)構(gòu)設(shè)計的指導(dǎo)。目前,數(shù)值模擬二維光子晶體帶隙的常用方法有平面波展開法[10-11]、時域有限差分法[2,12]和有限元方法[13],其中有限元方法因其計算效率較高,常用于三角晶格光子晶體和正方晶格光子晶體的數(shù)值模擬,但目前在蜂窩狀光子晶體領(lǐng)域研究較少。

    本文采用有限元法研究二維蜂窩狀光子晶體的能帶結(jié)構(gòu),以往在對蜂窩晶格光子晶體進(jìn)行建模時通常選擇正六邊形結(jié)構(gòu)作為單元晶格,本文在數(shù)值模擬過程中直接采用了易于操作的三角晶格胞元,使建模仿真和邊界處理更加容易。并利用COMSOL 軟件計算了二維蜂窩狀光子晶體的能帶結(jié)構(gòu),通過對比文獻(xiàn)[14]的結(jié)果,驗證了方法的正確性和有效性。進(jìn)一步,本文采用數(shù)值模擬的方法探究了二維蜂窩狀光子晶體的參數(shù)與帶隙的關(guān)系,并通過在蜂窩狀光子晶體中應(yīng)用新型材料構(gòu)造了寬帶隙的新型二維蜂窩狀光子晶體結(jié)構(gòu)。

    1 能帶結(jié)構(gòu)計算的數(shù)學(xué)模型

    在本節(jié)中,將介紹二維蜂窩狀光子晶體能帶求解的數(shù)學(xué)模型??紤]由周期介電圓柱組成的二維蜂窩狀光子晶體陣列,如圖1(a)所示,a0為晶格常數(shù),即相鄰2 個正六邊形中心的距離,a1和a2是長度a0的向量。放大的單個六角蜂窩晶格如圖1(b)所示,其中R 為六邊形的邊長,d 為圓柱體的直徑,εa和εb分別是周圍背景和介電圓柱的介電常數(shù)。本文采用三角晶格胞元進(jìn)行數(shù)值計算,光子晶體等效至xoy 平面后選取的三角形晶格胞元,如圖1(c)所示,其邊長為a0。

    圖1 蜂窩狀光子晶體示意圖

    圖2 三角晶格第一布里淵區(qū)

    下面根據(jù)圖1(c)所示的三角單元晶格,推導(dǎo)光子晶體特征值問題。在經(jīng)典電磁理論中,電磁現(xiàn)象可由Maxwell 方程組(1)描述。在光子晶體中,研究的對象是電磁波,其在光子晶體中傳播時滿足Maxwell 方程組

    式中:E 是電場強(qiáng)度;H 是磁場強(qiáng)度。對于二維光子晶體,直接對Maxwell 方程進(jìn)行求解比較困難,由于電磁波可分解為簡諧波,即有E(r,t)=E(r)e-iωt,H(r,t)=H(r)e-iωt,這里ω 是頻率。代入式(1)可得到與時間變量無關(guān)的Helmholtz 方程,即電場方程(2)和磁場方程(3)

    考慮二維光子晶體帶隙問題時,可簡化為TM 和TE 兩種偏振情形。在TM 模式下,有E(r)=(0,0,ER());在TE 模式下,有H(r)=(0,0,HR())。這里是二維坐標(biāo)(x,y)。代入式(2)和(3),分別可以得到TM 模的電場波動方程(4)和TE 模的磁場波動方程(5)

    由于二維光子晶體的周期性,滿足Bloch 定理,可將上述波動方程化至單位原胞上再進(jìn)行數(shù)值實驗。在二維情況下,記,k=(a,b)表示的是二維平面空間內(nèi)的波矢矢量,其取值范圍為不可約布里淵區(qū)B。首先考慮TM 模

    Ω 是求解的二維集合區(qū)域,在單位原胞邊界∏=?Ω上,滿足邊界條件

    式中:Γ1和Γ3是單位原胞的一組對邊;Γ2和Γ4是單位原胞的另一組對邊。求解區(qū)域如圖1(c)三角單元晶格所示。

    求解TE 模的過程類似,為方便討論變分過程,可將2 種模式統(tǒng)一記成如下特征值形式,給定矢量k,求λ 和uR使得

    uR表示待求場量即特征向量。只要矢量k 遍歷第一不可約布里淵區(qū)域的邊界,就能夠繪制出能帶結(jié)構(gòu)圖。

    下面針對式(8),給出其有限元變分過程和離散過程。首先,定義Sobolev 空間

    對于固定的波矢k,當(dāng)求解區(qū)域是線性三角元時,u、v在單元e 內(nèi)的近似ue、ve,取線性三角元的基函數(shù)和的線性組合,這里所說的基函數(shù)是關(guān)于(x,y)的勢函數(shù),不妨取試探函數(shù)ve為單元基函數(shù),則有

    這里x 取值{i,j,k},而ui,uj,uk是待定常數(shù)。

    根據(jù)SRSS評分,擬定療效評定標(biāo)準(zhǔn),分為顯效、有效、無效3個指標(biāo):(1)顯效:治療后,SRSS評分降低>50%。(2)有效:治療后,SRSS評分降低25%-50%。(3)無效:治療后,SRSS評分降低<25%[3]。

    由此得到單元剛度矩陣,然后根據(jù)線性三角元局部編號和整體編號關(guān)系,生成總剛度矩陣A 和B,得到變分問題對應(yīng)的離散特征值問題:

    在求解廣義特征值問題時,還需要考慮周期邊界條件(7),從而得到最終的待求解廣義特征值問題

    對于形如(14)的廣義特征值問題,在二維蜂窩狀光子晶體帶隙模型中,矩陣A 和B 常具有正定性,則有Cholesky 三角分解

    則對式(14)兩側(cè)的左側(cè)乘L-1,并在右側(cè)整理提出因子LT得

    式(16)可以看成以λ 為特征值,LT為特征向量的特征值問題,并且式(16)和(14)具有相同的特征值,系數(shù)矩陣仍滿足正定性,這樣就將廣義特征值問題(14)轉(zhuǎn)化為特征值問題(16),可利用直接法或迭代法進(jìn)行求解特征值問題(16)。

    2 數(shù)值仿真結(jié)果

    2.1 算例1:二維蜂窩狀光子晶體仿真

    本節(jié)考慮圖1 所示的二維蜂窩狀光子晶體,周圍背景的介電常數(shù)是εa=1,圓柱體的介電常數(shù)是εb=11.7,六邊形的長度為R,圓柱體截面圓形的直徑為d且d=2/3R。選擇三角單元晶格模擬得到的能帶圖如圖3 所示,在其他參數(shù)不變的情況下,僅改變晶格參數(shù)的大小,此時三角晶格中圓形的大小和位置會發(fā)生改變。對不同參數(shù)取值下的能帶結(jié)構(gòu)進(jìn)行仿真,結(jié)果如圖3所示,當(dāng)a0/R=3.125 時,帶隙歸一化頻率在0.48~0.51之間;當(dāng)a0/R=3 時,帶隙完全消失;當(dāng)a0/R=2.8 時,帶隙歸一化頻率在0.44~0.48 之間。3 種情況下的計算結(jié)果均和文獻(xiàn)[14]中的仿真結(jié)果一致,并且仿真所用時間比采用常見三角晶格進(jìn)行仿真時間更短、效率更高。

    圖3 二維蜂窩狀光子晶體的TM 模式的能帶圖

    圖4 TM 模式下改變晶格參數(shù)的Gap map 變化

    2.2 算例2:二維新結(jié)構(gòu)蜂窩狀光子晶體仿真

    下面,研究當(dāng)改變2.1 節(jié)中晶格結(jié)構(gòu)時光子晶體的能帶結(jié)構(gòu)。在每個六角蜂窩晶格的中心加入材質(zhì)為Si(介電常數(shù)為17.64),半徑r 同樣為d/2 的介電圓柱,周圍背景為空氣介電常數(shù)εa=1,周圍的圓柱介電常數(shù)仍為εb=11.7,每個晶格中中心圓柱的介電常數(shù)為εc=17.64,其余參數(shù)不變,對這種新結(jié)構(gòu)進(jìn)行仿真。新結(jié)構(gòu)蜂窩狀光子晶體示意圖如圖5 所示。

    圖5 新結(jié)構(gòu)蜂窩狀光子晶體示意圖

    取晶格參數(shù)a0/R=3 時為例,由2.1 節(jié)可知圖1 所示的蜂窩狀光子晶體在此晶格參數(shù)條件下并無任何帶隙產(chǎn)生,而新結(jié)構(gòu)的能帶結(jié)構(gòu)如圖6 所示,此時會產(chǎn)生1 條新的帶隙,與2.1 節(jié)中結(jié)構(gòu)相比新結(jié)構(gòu)光子晶體產(chǎn)生的帶隙位置向低頻方向移動。此時產(chǎn)生的帶隙歸一化頻率范圍為[0.264 954 9,0.279 614 2],寬度為0.014 659 2。

    圖6 新結(jié)構(gòu)蜂窩狀光子晶體得到的能帶圖

    為探究新結(jié)構(gòu)的帶隙隨參數(shù)的變化情況,類似2.1節(jié)中保持其他參數(shù)不變僅改變晶格參數(shù)a0/R,a0/R 的取值從2.7~7.7,得到新結(jié)構(gòu)TM 模式下蜂窩狀光子晶體的帶隙變化圖,如圖7 所示。由圖7 可知,當(dāng)晶格參數(shù)a0/R=5.6 時可產(chǎn)生最寬帶隙,此時有3 條帶隙,寬度分別為0.037 548,0.002 096 和0.134 911。

    圖7 改變晶格參數(shù)a0/R 得到的Gap map

    為進(jìn)一步探究蜂窩狀光子晶體帶隙的產(chǎn)生和變化情況,改變蜂窩晶格中心處介質(zhì)圓柱的截面半徑r,得到帶隙變化如圖8 所示,當(dāng)r 取0.75d 時可獲得最寬的帶隙,此時有4 條帶隙存在。

    圖8 改變中心圓柱截面半徑r 得到的Gap map

    2.3 算例3:二維蜂窩狀復(fù)合膠體光子晶體仿真

    當(dāng)前的光子晶體能帶研究主要是針對人為可控的物理參數(shù)的對比研究,但是對于不同尺寸粒子混合形成的膠體光子晶體的光子帶隙則少有研究。膠體光子晶體是一類利用交替納米顆粒通過自組裝而成的周期有序結(jié)構(gòu)的材料,可以保留2 種材料的光學(xué)特性,具有制備簡單、易于實現(xiàn)功能化和成本低廉等優(yōu)勢,可應(yīng)用于傳感器、防偽圖案和顯示器件等領(lǐng)域。這種光子晶體由于其獨(dú)特的結(jié)構(gòu)特征及光學(xué)特性的廣泛應(yīng)用,近年來已成為光子晶體領(lǐng)域的研究熱點(diǎn)。

    采用這種新型材料,由大粒徑聚苯乙烯(PS)和小粒徑二氧化硅(SiO2)組成二維復(fù)合膠體光子晶體,其結(jié)構(gòu)如圖9 所示,背景介質(zhì)為空氣εa=1;在此種結(jié)構(gòu)中大的微球(PS)介電常數(shù)為2.53,半徑為d/2;小的微球(SiO2)介電常數(shù)為2.13,其半徑為r=d/4。本文基于有限元法,利用COMSOL 軟件的波動光學(xué)模塊建立光子晶體模型,并使用Floquet 周期性邊界條件,探究不同參數(shù)下膠體光子晶體的帶隙寬度。

    利用COMSOL 軟件對圖9 所示的膠體光子晶體的光子帶隙進(jìn)行了仿真模擬,結(jié)果如圖10 所示,此時在[0.847 558 0,0.840 933 8]的歸一化頻率的區(qū)間產(chǎn)生一條寬度為0.006 624 2 的帶隙。

    圖9 由PS 和SiO2 組成的二維復(fù)合膠體光子晶體結(jié)構(gòu)

    圖10 二維復(fù)合膠體光子晶體仿真結(jié)果圖

    采用二維復(fù)合膠體光子晶體可以產(chǎn)生如上所述的帶隙,而如果僅使用一種材料PS,此時能帶結(jié)構(gòu)無任何帶隙產(chǎn)生,如圖11 所示。由此可見,引入二維復(fù)合膠體材料能有效地增大蜂窩狀光子晶體的帶隙寬度,同時還可以保留2 種材料的光學(xué)特性。

    圖11 僅由PS 材料構(gòu)成的光子晶體能帶圖

    取上面討論的最佳參數(shù):晶格常數(shù)選取a0/R=5.6,中心圓柱界面半徑r 取0.75d 時,可獲得寬帶隙,此時應(yīng)用新材料的蜂窩狀光子晶體的能帶圖如圖12 所示,此時會產(chǎn)生3 條帶隙。第一條帶隙的歸一化頻率范圍為[0.467 225,0.487 643],帶寬為0.020 717;第二條帶隙的歸一化頻率范圍為[0.827 783,0.831 438],帶寬為0.003 655;第三條帶隙的歸一化頻率范圍為[0.929 49,0.953 844],帶寬為0.024 347。由圖10 和圖12 比較可知,選取合適的晶格參數(shù)a0/R 和圓形半徑r來獲得寬帶隙的方法是有效的。

    圖12 最優(yōu)參數(shù)下二維復(fù)合膠體光子晶體仿真結(jié)果圖

    3 結(jié)論

    本文采用了有限元方法,研究了二維蜂窩狀光子晶體的能帶結(jié)構(gòu)數(shù)值模擬和寬帶隙的結(jié)構(gòu)設(shè)計問題。對有限元方法和COMSOL 軟件的結(jié)合運(yùn)用,使得光子晶體能帶結(jié)構(gòu)的求解變得簡單高效。數(shù)值模擬方法是新材料運(yùn)用、新結(jié)構(gòu)設(shè)計的一種非常高效、便捷的方法,通過優(yōu)化蜂窩狀光子晶體的晶格參數(shù)及應(yīng)用新材料,基于數(shù)值模擬方法探究了二維蜂窩狀光子晶體的參數(shù)與帶隙的關(guān)系,通過選取最佳參數(shù)構(gòu)造了寬帶隙的新型二維蜂窩狀光子晶體結(jié)構(gòu)。本文的研究為二維蜂窩狀光子晶體的能帶結(jié)構(gòu)求解提供了一種有效求解途徑,有望為光子晶體設(shè)備的制備提供指導(dǎo)和數(shù)值算法支撐。

    猜你喜歡
    蜂窩狀帶隙能帶
    酷颯主場
    密度泛函理論計算半導(dǎo)體材料的帶隙誤差研究
    吃東西時注意多
    顆粒和蜂窩狀廢棄活性炭再生方法的探究
    硫碘循環(huán)制氫系統(tǒng)中蜂窩狀活性炭成型工藝研究及活性評價
    能源工程(2020年6期)2021-01-26 00:55:10
    汽車轉(zhuǎn)向管柱吸能帶變形研究和仿真優(yōu)化
    一種基于BJT工藝的無運(yùn)放低溫度系數(shù)的帶隙基準(zhǔn)源
    Life OR Death Decision
    間距比對雙振子局域共振軸縱振帶隙的影響
    一款高PSRR低溫度系數(shù)的帶隙基準(zhǔn)電壓源的設(shè)計
    電子制作(2018年1期)2018-04-04 01:48:38
    亚洲欧美成人精品一区二区| 亚洲三级黄色毛片| 色哟哟·www| 亚洲中文字幕日韩| 亚洲国产精品国产精品| 性色avwww在线观看| 国产成年人精品一区二区| av.在线天堂| 99久久精品一区二区三区| 日本欧美国产在线视频| 丰满人妻一区二区三区视频av| 不卡视频在线观看欧美| 成人鲁丝片一二三区免费| 欧美激情在线99| 国产一区二区在线观看日韩| 欧美激情久久久久久爽电影| 99热这里只有精品一区| 国产麻豆成人av免费视频| av国产久精品久网站免费入址| 六月丁香七月| 69av精品久久久久久| 日韩av在线大香蕉| 91在线精品国自产拍蜜月| 国产黄a三级三级三级人| 麻豆久久精品国产亚洲av| 午夜福利成人在线免费观看| 久久久久久久久大av| 亚洲国产精品sss在线观看| 人人妻人人澡欧美一区二区| 91午夜精品亚洲一区二区三区| av福利片在线观看| 欧美最新免费一区二区三区| 欧美激情国产日韩精品一区| 国产午夜精品论理片| 欧美变态另类bdsm刘玥| 男女那种视频在线观看| 久久久a久久爽久久v久久| 国产又黄又爽又无遮挡在线| 亚洲天堂国产精品一区在线| 免费黄网站久久成人精品| 插阴视频在线观看视频| 亚洲最大成人av| 欧美日韩精品成人综合77777| 午夜福利在线观看免费完整高清在| 亚洲在线自拍视频| 午夜激情福利司机影院| 听说在线观看完整版免费高清| 97人妻精品一区二区三区麻豆| 黄色日韩在线| 欧美97在线视频| 麻豆乱淫一区二区| 久久久久九九精品影院| av视频在线观看入口| 我要看日韩黄色一级片| av又黄又爽大尺度在线免费看 | 男女边吃奶边做爰视频| 搡女人真爽免费视频火全软件| 成人三级黄色视频| 美女黄网站色视频| 变态另类丝袜制服| 久久人妻av系列| 菩萨蛮人人尽说江南好唐韦庄 | 国产亚洲精品av在线| 色尼玛亚洲综合影院| 国产精品国产三级国产专区5o | 白带黄色成豆腐渣| 欧美变态另类bdsm刘玥| 99热这里只有精品一区| 国产精品一二三区在线看| 韩国高清视频一区二区三区| 直男gayav资源| 日日摸夜夜添夜夜爱| 两性午夜刺激爽爽歪歪视频在线观看| 69av精品久久久久久| 欧美三级亚洲精品| 欧美激情国产日韩精品一区| 中文字幕久久专区| 久久精品国产鲁丝片午夜精品| 日本黄色片子视频| 一区二区三区免费毛片| 日本免费a在线| 色综合站精品国产| 亚洲av免费在线观看| 免费播放大片免费观看视频在线观看 | 伦精品一区二区三区| 啦啦啦韩国在线观看视频| 亚洲在久久综合| 熟妇人妻久久中文字幕3abv| 精品国内亚洲2022精品成人| 精品久久久久久久久久久久久| 久热久热在线精品观看| 一级黄片播放器| 国产真实乱freesex| 欧美成人一区二区免费高清观看| 嫩草影院新地址| 青春草亚洲视频在线观看| 日本猛色少妇xxxxx猛交久久| 美女黄网站色视频| 亚洲美女视频黄频| 日本午夜av视频| 欧美日本亚洲视频在线播放| 精品人妻偷拍中文字幕| 国产免费福利视频在线观看| 亚洲成人av在线免费| 伦理电影大哥的女人| 日韩成人伦理影院| 三级男女做爰猛烈吃奶摸视频| 国产成人精品久久久久久| 国产美女午夜福利| 国产 一区 欧美 日韩| 日韩,欧美,国产一区二区三区 | 亚洲综合精品二区| 简卡轻食公司| 国产高潮美女av| 身体一侧抽搐| 麻豆av噜噜一区二区三区| 亚洲伊人久久精品综合 | av在线蜜桃| 精品一区二区三区人妻视频| 国产 一区精品| 国产精品综合久久久久久久免费| 国产精品人妻久久久久久| 久久久久久久久久成人| 国产亚洲av片在线观看秒播厂 | 欧美成人精品欧美一级黄| 最近的中文字幕免费完整| 毛片一级片免费看久久久久| 日韩成人av中文字幕在线观看| 精品午夜福利在线看| 最近最新中文字幕免费大全7| 日本五十路高清| 日本免费在线观看一区| 特级一级黄色大片| 人人妻人人看人人澡| 美女被艹到高潮喷水动态| 亚洲成av人片在线播放无| 联通29元200g的流量卡| 国产精品伦人一区二区| 国模一区二区三区四区视频| 成人av在线播放网站| 一本久久精品| 免费看av在线观看网站| 成人综合一区亚洲| 熟女人妻精品中文字幕| 日韩欧美三级三区| 少妇的逼水好多| 18禁动态无遮挡网站| 最近的中文字幕免费完整| 乱系列少妇在线播放| 乱码一卡2卡4卡精品| 三级经典国产精品| 国产伦精品一区二区三区视频9| 国产精品久久视频播放| 内地一区二区视频在线| 国产精品一区二区在线观看99 | 国产精品1区2区在线观看.| 欧美bdsm另类| 久久精品久久精品一区二区三区| 秋霞在线观看毛片| 99热这里只有精品一区| h日本视频在线播放| 久久精品久久久久久噜噜老黄 | 久久人妻av系列| 看免费成人av毛片| 亚洲乱码一区二区免费版| 亚洲伊人久久精品综合 | 一边摸一边抽搐一进一小说| 成人二区视频| 久久亚洲国产成人精品v| 久久久久久久久中文| 人体艺术视频欧美日本| 汤姆久久久久久久影院中文字幕 | 寂寞人妻少妇视频99o| 国内精品美女久久久久久| 日韩中字成人| 亚洲精品乱久久久久久| 亚洲内射少妇av| 久久99热这里只频精品6学生 | 又粗又硬又长又爽又黄的视频| 亚洲欧美精品专区久久| 国产av码专区亚洲av| 97超碰精品成人国产| 精品99又大又爽又粗少妇毛片| 久久久色成人| 最近中文字幕高清免费大全6| 日韩三级伦理在线观看| 爱豆传媒免费全集在线观看| 只有这里有精品99| 欧美精品国产亚洲| av视频在线观看入口| 美女黄网站色视频| 哪个播放器可以免费观看大片| 中文亚洲av片在线观看爽| 久久婷婷人人爽人人干人人爱| 久久亚洲精品不卡| videossex国产| 国产在视频线在精品| 精品国产露脸久久av麻豆 | 国产极品天堂在线| 免费不卡的大黄色大毛片视频在线观看 | 亚洲经典国产精华液单| 亚洲中文字幕日韩| 久久久久国产网址| 亚洲欧美日韩无卡精品| 久久99热这里只频精品6学生 | 亚洲四区av| 国产日韩欧美在线精品| 亚洲精品aⅴ在线观看| 亚洲欧洲日产国产| 51国产日韩欧美| 亚洲成av人片在线播放无| 亚洲欧美精品自产自拍| 国产成人免费观看mmmm| 成人午夜高清在线视频| 亚洲一区高清亚洲精品| 波多野结衣巨乳人妻| 亚洲中文字幕日韩| 日本与韩国留学比较| 噜噜噜噜噜久久久久久91| 一本久久精品| 久久久a久久爽久久v久久| 欧美最新免费一区二区三区| 国产一区二区在线观看日韩| 欧美性感艳星| 在线免费十八禁| 秋霞伦理黄片| 成人毛片60女人毛片免费| 亚洲av成人av| 熟女人妻精品中文字幕| 午夜视频国产福利| 中文天堂在线官网| 男女边吃奶边做爰视频| 久久精品夜色国产| 搞女人的毛片| 国产av在哪里看| 丰满乱子伦码专区| 最近2019中文字幕mv第一页| 亚洲在线自拍视频| 在线a可以看的网站| 日日撸夜夜添| 精品久久久久久久末码| 中文字幕免费在线视频6| 91精品一卡2卡3卡4卡| 国产极品天堂在线| 亚洲av电影在线观看一区二区三区 | 欧美日韩综合久久久久久| 久久99热这里只有精品18| 能在线免费看毛片的网站| 级片在线观看| 亚洲av日韩在线播放| 日韩欧美三级三区| 免费观看的影片在线观看| 久久久精品大字幕| 国产免费一级a男人的天堂| 五月玫瑰六月丁香| 国产黄色视频一区二区在线观看 | 亚洲国产欧美人成| 51国产日韩欧美| 欧美97在线视频| 人妻夜夜爽99麻豆av| 国产黄a三级三级三级人| 国产精品一区www在线观看| 人人妻人人看人人澡| 久久99蜜桃精品久久| 久久久精品94久久精品| 成年免费大片在线观看| 亚洲欧美日韩东京热| 免费av毛片视频| 91午夜精品亚洲一区二区三区| 日本午夜av视频| 我要搜黄色片| 非洲黑人性xxxx精品又粗又长| 欧美+日韩+精品| 久久这里只有精品中国| 91精品伊人久久大香线蕉| 最近最新中文字幕大全电影3| 国产极品精品免费视频能看的| 97在线视频观看| 免费看日本二区| 在线观看66精品国产| 长腿黑丝高跟| 精品无人区乱码1区二区| 欧美激情国产日韩精品一区| 九九久久精品国产亚洲av麻豆| 乱系列少妇在线播放| 国产人妻一区二区三区在| 欧美一级a爱片免费观看看| 日本五十路高清| 国产高清有码在线观看视频| 免费黄网站久久成人精品| 午夜福利在线观看免费完整高清在| 久久久久久久午夜电影| 又爽又黄无遮挡网站| 精品国产露脸久久av麻豆 | 一区二区三区高清视频在线| 欧美性感艳星| 直男gayav资源| 欧美精品一区二区大全| 国产免费视频播放在线视频 | 久久久精品大字幕| 成人毛片60女人毛片免费| 性插视频无遮挡在线免费观看| 老女人水多毛片| 亚洲精品aⅴ在线观看| 色综合站精品国产| 亚洲av不卡在线观看| 亚洲国产高清在线一区二区三| or卡值多少钱| 18禁动态无遮挡网站| 国产精品久久久久久久电影| a级毛片免费高清观看在线播放| 女人被狂操c到高潮| 麻豆成人av视频| 最新中文字幕久久久久| 男人和女人高潮做爰伦理| 欧美日本视频| 精品一区二区免费观看| 在线免费十八禁| 国产探花极品一区二区| 日日撸夜夜添| 99热全是精品| 女人久久www免费人成看片 | 中文字幕久久专区| 黄片无遮挡物在线观看| 日韩av在线大香蕉| 国产色爽女视频免费观看| 国产三级中文精品| 午夜老司机福利剧场| 小蜜桃在线观看免费完整版高清| 久久精品久久久久久噜噜老黄 | 两性午夜刺激爽爽歪歪视频在线观看| 久久久久久伊人网av| 国产淫片久久久久久久久| 国产探花极品一区二区| 日韩欧美三级三区| 美女高潮的动态| 欧美日韩一区二区视频在线观看视频在线 | 性色avwww在线观看| 国产伦精品一区二区三区四那| 五月伊人婷婷丁香| 亚洲欧美日韩高清专用| 国产伦在线观看视频一区| 国产私拍福利视频在线观看| 国产乱来视频区| 亚洲成人av在线免费| 中文字幕熟女人妻在线| av黄色大香蕉| 免费看a级黄色片| 亚洲精品国产av成人精品| 中文资源天堂在线| 美女黄网站色视频| 桃色一区二区三区在线观看| 99视频精品全部免费 在线| 中文字幕熟女人妻在线| 又爽又黄a免费视频| av在线观看视频网站免费| 中文乱码字字幕精品一区二区三区 | 最近2019中文字幕mv第一页| 免费无遮挡裸体视频| 国国产精品蜜臀av免费| videos熟女内射| 91av网一区二区| 熟女人妻精品中文字幕| h日本视频在线播放| 国产精品三级大全| 国产又黄又爽又无遮挡在线| 看非洲黑人一级黄片| 精品一区二区三区视频在线| 岛国毛片在线播放| 成年女人永久免费观看视频| 国产久久久一区二区三区| 亚洲精品日韩在线中文字幕| 丰满少妇做爰视频| 欧美人与善性xxx| 国国产精品蜜臀av免费| 内地一区二区视频在线| 日韩一本色道免费dvd| 欧美xxxx黑人xx丫x性爽| 久久久久免费精品人妻一区二区| 国产成人免费观看mmmm| 干丝袜人妻中文字幕| 1000部很黄的大片| 久久6这里有精品| 中文字幕av在线有码专区| 久久久成人免费电影| 国产亚洲一区二区精品| 日本黄大片高清| 一二三四中文在线观看免费高清| 国产探花在线观看一区二区| av国产免费在线观看| 精品人妻熟女av久视频| 欧美日韩综合久久久久久| 午夜亚洲福利在线播放| 成人欧美大片| av免费在线看不卡| 国产单亲对白刺激| 纵有疾风起免费观看全集完整版 | 亚洲av免费在线观看| 亚洲国产成人一精品久久久| 久久精品人妻少妇| 男插女下体视频免费在线播放| 亚洲欧美日韩东京热| 91久久精品国产一区二区三区| 嫩草影院新地址| 直男gayav资源| 日韩av不卡免费在线播放| 嫩草影院入口| 欧美日本视频| 亚洲欧洲日产国产| 边亲边吃奶的免费视频| 97热精品久久久久久| 国产成人一区二区在线| 床上黄色一级片| 亚洲精品乱码久久久v下载方式| 2021天堂中文幕一二区在线观| 国产亚洲精品久久久com| 久久久久久久久久久免费av| av在线天堂中文字幕| 一级爰片在线观看| 草草在线视频免费看| 久久久国产成人精品二区| 日韩人妻高清精品专区| 亚洲综合精品二区| 成年版毛片免费区| 视频中文字幕在线观看| 欧美高清成人免费视频www| 国产精品国产三级国产av玫瑰| 中文字幕制服av| 欧美又色又爽又黄视频| 99久国产av精品国产电影| 国内精品一区二区在线观看| 色哟哟·www| 亚洲天堂国产精品一区在线| 91aial.com中文字幕在线观看| 国产老妇伦熟女老妇高清| 白带黄色成豆腐渣| 嫩草影院精品99| 精品国内亚洲2022精品成人| 国产私拍福利视频在线观看| 欧美成人免费av一区二区三区| 亚洲欧美精品专区久久| 精品无人区乱码1区二区| 麻豆久久精品国产亚洲av| 午夜福利在线观看吧| 免费看光身美女| 欧美一区二区亚洲| 天堂av国产一区二区熟女人妻| 国产男人的电影天堂91| 丝袜喷水一区| 免费不卡的大黄色大毛片视频在线观看 | 亚洲电影在线观看av| 有码 亚洲区| 亚洲在线观看片| 久久久久久久久大av| 99热这里只有精品一区| 哪个播放器可以免费观看大片| 人妻系列 视频| 日韩人妻高清精品专区| 亚洲欧美日韩高清专用| 你懂的网址亚洲精品在线观看 | 在线a可以看的网站| 一二三四中文在线观看免费高清| 97超碰精品成人国产| 一级毛片电影观看 | 亚洲精品日韩在线中文字幕| 亚洲国产精品成人久久小说| 熟女电影av网| 97热精品久久久久久| 久久久久久久久久久免费av| av卡一久久| 国产成人freesex在线| 深爱激情五月婷婷| 亚洲精华国产精华液的使用体验| 国产精品三级大全| 国产淫语在线视频| 亚洲国产欧美在线一区| 十八禁国产超污无遮挡网站| 七月丁香在线播放| 又爽又黄a免费视频| 精品久久久久久电影网 | 亚洲伊人久久精品综合 | 天堂中文最新版在线下载 | 九九在线视频观看精品| 亚洲aⅴ乱码一区二区在线播放| 搞女人的毛片| 日韩成人av中文字幕在线观看| 国产午夜精品论理片| 一级二级三级毛片免费看| 18+在线观看网站| 九九热线精品视视频播放| a级一级毛片免费在线观看| 日本猛色少妇xxxxx猛交久久| 亚洲三级黄色毛片| 长腿黑丝高跟| 啦啦啦观看免费观看视频高清| 爱豆传媒免费全集在线观看| 亚洲av成人精品一区久久| 国产成人91sexporn| 国产av一区在线观看免费| 美女国产视频在线观看| 亚洲av不卡在线观看| 婷婷色综合大香蕉| 2021天堂中文幕一二区在线观| 国产乱来视频区| 亚洲乱码一区二区免费版| 国产精品一二三区在线看| 久久99蜜桃精品久久| 女人被狂操c到高潮| 成人三级黄色视频| 国产精品蜜桃在线观看| 国产一区二区在线观看日韩| 午夜精品国产一区二区电影 | 午夜福利成人在线免费观看| 黄片wwwwww| 高清毛片免费看| 99热网站在线观看| 国产综合懂色| 国产精品国产三级国产av玫瑰| av在线老鸭窝| 亚洲av.av天堂| 男女那种视频在线观看| 久久久久久久久久黄片| 日韩视频在线欧美| 国产一区二区在线av高清观看| 国产欧美日韩精品一区二区| 久久久久久久午夜电影| 18+在线观看网站| 少妇裸体淫交视频免费看高清| 黄色一级大片看看| 成人无遮挡网站| 欧美三级亚洲精品| 99在线视频只有这里精品首页| 国产精品综合久久久久久久免费| 精品人妻熟女av久视频| 久久欧美精品欧美久久欧美| 韩国高清视频一区二区三区| 深夜a级毛片| 成年av动漫网址| 亚洲怡红院男人天堂| 国产综合懂色| 久久精品影院6| 久久久久久久久中文| 超碰av人人做人人爽久久| 亚洲在久久综合| 男女视频在线观看网站免费| 欧美日韩综合久久久久久| 大香蕉97超碰在线| 日韩欧美 国产精品| 一级黄片播放器| 国产一区二区三区av在线| 久久久久久久国产电影| 菩萨蛮人人尽说江南好唐韦庄 | 又粗又爽又猛毛片免费看| 韩国高清视频一区二区三区| 亚洲成人中文字幕在线播放| 中文字幕亚洲精品专区| 婷婷六月久久综合丁香| 男女那种视频在线观看| 中文精品一卡2卡3卡4更新| 看片在线看免费视频| 美女xxoo啪啪120秒动态图| 熟妇人妻久久中文字幕3abv| 少妇猛男粗大的猛烈进出视频 | 午夜福利视频1000在线观看| 午夜福利在线观看免费完整高清在| 日韩亚洲欧美综合| 中文天堂在线官网| 国产精品熟女久久久久浪| 欧美日韩一区二区视频在线观看视频在线 | 好男人视频免费观看在线| 插逼视频在线观看| 久久久国产成人免费| 欧美又色又爽又黄视频| 一级av片app| 午夜a级毛片| 久久久久久久久久久丰满| 亚洲自偷自拍三级| 亚洲三级黄色毛片| 亚洲精品国产成人久久av| 观看免费一级毛片| 天堂√8在线中文| 天堂av国产一区二区熟女人妻| 国产精品日韩av在线免费观看| 亚洲精华国产精华液的使用体验| 偷拍熟女少妇极品色| 大香蕉久久网| 亚洲国产最新在线播放| 大话2 男鬼变身卡| 国产麻豆成人av免费视频| 人人妻人人澡欧美一区二区| 国产精品久久久久久久久免| 免费播放大片免费观看视频在线观看 | 极品教师在线视频| 草草在线视频免费看| 男女那种视频在线观看| 国产黄a三级三级三级人| 国产欧美日韩精品一区二区| 一级二级三级毛片免费看| 亚洲欧美日韩东京热| videossex国产| 国内精品美女久久久久久| 日韩 亚洲 欧美在线| 国产精品福利在线免费观看| 亚洲人成网站在线播| 亚洲av免费在线观看| 久久精品熟女亚洲av麻豆精品 | 亚洲真实伦在线观看| 啦啦啦啦在线视频资源| 亚洲欧美成人综合另类久久久 | 久久久成人免费电影| 亚洲在线观看片| 乱人视频在线观看| 国产精品国产高清国产av| 午夜免费激情av| 91久久精品电影网|