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

    考慮材料參數(shù)的二維固/固聲子晶體拓?fù)鋬?yōu)化設(shè)計

    2019-12-10 02:44:20許衛(wèi)鍇寧金英韓世超余欣宇祁武超
    關(guān)鍵詞:散射體帶隙聲子

    許衛(wèi)鍇,寧金英,韓世超,余欣宇,祁武超

    (沈陽航空航天大學(xué) 航空宇航學(xué)院,沈陽 110136)

    聲子晶體[1]是在光子晶體的基礎(chǔ)上發(fā)展起來的,由不同彈性性質(zhì)的材料按結(jié)構(gòu)周期性復(fù)合在一起形成的人工復(fù)合材料。這種周期性結(jié)構(gòu)由于帶隙的存在及結(jié)構(gòu)的可操縱性,使其廣泛應(yīng)用于頻率濾波器、波導(dǎo)和振動保護裝置等方面[2,3]。近年來,對聲子晶體帶隙的拓?fù)鋬?yōu)化已成為該領(lǐng)域的研究熱點。Huang等[4]提出基于梯度變換的優(yōu)化方法,該方法以波數(shù)的余弦代替衰減系數(shù)作為目標(biāo)函數(shù)以分析彈性波在特定頻率范圍內(nèi)的衰減特性;Liu等[5]利用兩階段的遺傳算法來尋找二維方形晶格聲子晶體規(guī)定頻帶之間的最大相對帶隙;Lin等[6]提出一種用于獲得聲子晶體形狀拓?fù)鋬?yōu)化的靈敏度導(dǎo)數(shù)的時間伴隨方法,實現(xiàn)了有效的降噪結(jié)果;Dong等利用遺傳算法和NSGA-II分別研究了二維正方形固/固聲子晶體的無約束和有約束的單目標(biāo)優(yōu)化設(shè)計[7]以及多孔單相聲子晶體的多目標(biāo)優(yōu)化設(shè)計[8],并獲得了較寬的相對帶隙;Li等[9]提出一種新的基于雙向進(jìn)化的結(jié)構(gòu)優(yōu)化(BESO)方法和有限元分析的拓?fù)鋬?yōu)化方法,并獲得了新的聲子晶體拓?fù)錁?gòu)型;Zhang等[10]提出一種獨立散射模型來評估材料的反平面剪切波衰減,并利用該模型設(shè)計了一種由平行排列的PMMA薄板組成的高質(zhì)量濾波器[11];He等[12]基于聲子晶體方向帶隙的拓?fù)鋬?yōu)化微結(jié)構(gòu)設(shè)計,實現(xiàn)了彈性波的單向傳輸特性;Xu等[13]研究了多相固態(tài)聲子晶體的拓?fù)鋬?yōu)化設(shè)計,并得到了更寬的最大化相對帶隙;劉堅等[14]將不確定模型引入局域共振聲子晶體,并采用Chebyshev展開法以帶隙最大化為目標(biāo)進(jìn)行優(yōu)化,數(shù)值分析和優(yōu)化后的結(jié)果都表明帶隙有大幅度的拓寬。

    聲子晶體作為一種復(fù)合介質(zhì),影響其帶隙特性的因素主要分為結(jié)構(gòu)參數(shù)和材料參數(shù)兩類。結(jié)構(gòu)參數(shù)主要包括散射體的形狀、晶格形式、散射體填充率等,體現(xiàn)為聲子晶體的拓?fù)湫问?;另一方面,材料參?shù)包括各組元材料的密度、彈性模量、聲阻抗以及波速等[1,15],同樣會對聲子晶體的帶隙特征產(chǎn)生影響。趙浩江等[16]分析了材料參數(shù)對薄板型聲子晶體的影響,得出散射體和基體材料的楊氏模量比和質(zhì)量密度比是影響薄板彎曲振動帶隙的關(guān)鍵材料參數(shù);張昭等[17]通過平面波展開法討論了薄板型聲子晶體的帶隙特性,并討論了散射體彈性模量比、密度比等對聲子晶體帶隙特性和減振的影響。趙元年等[18]利用有限元法研究了超晶格下材料參數(shù)對二維聲子晶體振動缺陷態(tài)的影響。Han等[19]通過選取4種組份材料,討論了材料參數(shù)對聲子晶體帶隙的影響,發(fā)現(xiàn)在改變材料參數(shù)時,得到的每兩個模型之間都有明顯差別。Li等[20]研究了材料參數(shù)對三相聲子晶體帶隙的影響,并討論了布拉格散射帶隙與局域共振帶隙的不同。

    然而,目前對聲子晶體帶隙特征的研究與優(yōu)化多是從一種影響因素出發(fā),或者是討論特定材料組合的拓?fù)錁?gòu)型,或者是固定構(gòu)型后討論材料參數(shù)的影響,并未考慮聲子晶體在兩個因素共同影響下的帶隙特征。本文基于混合變量的遺傳算法,研究同時考慮材料參數(shù)及拓?fù)錁?gòu)型對固/固聲子晶體帶隙的影響,從而達(dá)到材料參數(shù)及拓?fù)錁?gòu)型的一體化設(shè)計。

    1 理論方法

    1.1 聲子晶體彈性波方程

    考慮由兩種各向同性彈性固體材料組成的正方晶格的二維聲子晶體,其彈性波動方程[21]為

    (1)

    式(1)中,λ和μ是Lame常數(shù),r為位置矢量,ρ、u和ω分別是密度、位移矢量和角頻率。▽表示梯度算子。

    假設(shè)彈性波在xy平面內(nèi)傳播,則式(1)可以分解為平面混合模態(tài)和反平面剪切模式兩個方程

    (2)

    (3)

    由Bloch定理,位移場可表示為

    u(r)=ei(k·r)uk(r)

    (4)

    式(4)中uk(r)是與單胞結(jié)構(gòu)具有相同周期性的函數(shù),且k=(kx,ky)是限于第一布里淵區(qū)域的平面矢量。因此,單胞內(nèi)離散的特征值方程為

    [K(k)-ω2M]U=0

    (5)

    式(5)中K,U,M分別表示剛度矩陣、位移矢量和質(zhì)量矩陣。利用有限元軟件COMSOL,可方便有效地求解該本征頻率問題,并得到聲子晶體的帶隙結(jié)果。

    1.2 基于混合變量的遺傳算法

    遺傳算法[22]是一種全局、隨機和通用的智能算法,被廣泛地應(yīng)用于各種連續(xù)體/離散體的拓?fù)鋬?yōu)化領(lǐng)域。與其他常用的拓?fù)鋬?yōu)化方法相比,遺傳算法不依賴于目標(biāo)函數(shù)的梯度及敏度,具有非常強的魯棒性,因此特別適合求解非常復(fù)雜的波動問題,并在聲子晶體的優(yōu)化設(shè)計領(lǐng)域得到了廣泛地應(yīng)用[5,7-9]。

    目前應(yīng)用遺傳算法進(jìn)行的聲子晶體拓?fù)鋬?yōu)化更多地體現(xiàn)在拓?fù)錁?gòu)型的設(shè)計上,即以二進(jìn)制編碼0和1分別代表不同的材料,并通過編碼、解碼、選擇、交叉和變異等操作實現(xiàn)最優(yōu)值的求解。這種方法的優(yōu)點是編碼方式簡單、相對獨立,使得適應(yīng)度評價過程易于實現(xiàn)。然而,一般情況下0和1只能代表固定的材料,并沒有考慮材料參數(shù)對帶隙的影響。

    事實上,材料參數(shù)對聲子晶體的帶隙特性也具有重要的影響[15-20]。為此,本文引入密度、彈性模量等材料參數(shù)作為連續(xù)設(shè)計變量,使得該優(yōu)化問題變?yōu)橥瑫r考慮離散拓?fù)?二進(jìn)制變量)和材料參數(shù)(連續(xù)變量)的混合變量問題。求解此問題需要對變量進(jìn)行重新編碼:染色體分為兩部分,第一部分仍然是代表材料分布的離散變量,第二部分則是將代表材料參數(shù)的連續(xù)變量進(jìn)行轉(zhuǎn)化得到。

    例如,假定給定參數(shù)區(qū)間為[a,b],可將二進(jìn)制串與區(qū)間[a,b]內(nèi)的十進(jìn)制數(shù)建立映射關(guān)系。首先,可用式(6)將二進(jìn)制串(bnbn-1b0)2代表的二進(jìn)制數(shù)轉(zhuǎn)化為十進(jìn)制數(shù)x′;其次,應(yīng)用式(7)將x′轉(zhuǎn)化為對應(yīng)區(qū)間[a,b]內(nèi)的實數(shù)。

    (6)

    (7)

    式(7)中n+1為代表材料參數(shù)二進(jìn)制串的長度。

    鑒于前人的研究結(jié)果[15-20],不失一般性,可取散射體和基體的參數(shù)比值作為設(shè)計變量,分別討論使它們單獨變化及同時變化時對聲子晶體帶隙特性的影響。

    圖1 二維正方晶格不可約布里淵區(qū)

    2 算例分析

    本文選取四方對稱的聲子晶體進(jìn)行優(yōu)化,其不可約布里淵區(qū)示意圖如圖1所示。晶胞尺寸為0.02 m。為簡單起見,選擇目前聲子晶體中常用的環(huán)氧樹脂作為基體材料,其材料參數(shù)為密度ρ=1 200 kg/m3,楊氏模量E=4.51 GPa,泊松比υ=0.42。為了與文獻(xiàn)[7]、[23]的結(jié)果進(jìn)行對比,初始變量選擇鉛作為散射體,其材料參數(shù)為密度ρ=11 350 kg/m3,楊氏模量E=16 GPa,泊松比υ=0.42。為了減少計算成本,本文采用兩階段遺傳算法。在第一階段,晶胞首先被劃分為16×16的粗網(wǎng)格,得到的結(jié)果在第二階段被重新劃分為32×32的細(xì)網(wǎng)格并作為種子結(jié)構(gòu)參與第二階段的優(yōu)化。針對四方對稱的聲子晶體來說,只有1/8的區(qū)域需要進(jìn)行0、1編碼,因此在第一階段和第二階段所對應(yīng)的用于材料分布編碼的染色體長度分別為36和136,而對應(yīng)于連續(xù)變量的染色體長度要依據(jù)所選變量的范圍和求解精度來確定,總體染色體的長度為二者之和。遺傳算法相關(guān)參數(shù)為:種群大小為30,交叉率為0.9,突變率的初始值設(shè)為0.1。為簡單起見,本文僅以平面內(nèi)混合波模態(tài)為例展開討論,并選取相對帶寬(帶隙寬度/帶隙中間頻率)為優(yōu)化目標(biāo),即

    (8)

    2.1 考慮密度變化對帶隙的影響

    以散射體的密度作為連續(xù)變量進(jìn)行優(yōu)化時,根據(jù)常見材料的密度值,選取密度范圍為1 135~19 500 kg/m3,其中密度的上限與金相同。優(yōu)化結(jié)果如圖2所示,其中能帶圖的縱坐標(biāo)采用歸一化頻率Ω=ωα/2πCt,Epoxy19 500 kg/m3。

    從圖2中可以看出,優(yōu)化后的材料構(gòu)型與文獻(xiàn)[23]類似,為含4個圓角的近似方形的散射體。計算得到的最大相對帶寬為0.36,略低于文獻(xiàn)[7]、[23]的結(jié)果,這是由于我們的計算條件有限,所選取的離散單元數(shù)較少而引起的。另一方面,優(yōu)化后散射體密度為設(shè)置的密度范圍的上限,即19 500 kg/m3,同時帶隙位置有所降低,這表明散射體與基體的密度差越大越易產(chǎn)生較寬的低頻帶隙[16]。然而,當(dāng)彈性模量的差值未達(dá)到一定量級的時候,僅僅改變材料的密度并不能顯著提高聲子晶體的最大相對帶寬。例如,針對圖2的拓?fù)錁?gòu)型,將散射體重新設(shè)置為鉛,計算得到的相對帶寬為0.3。

    圖2 考慮密度變化時的優(yōu)化結(jié)果

    2.2 考慮彈性模量變化對帶隙的影響

    當(dāng)以散射體的彈性模量作為連續(xù)變量進(jìn)行優(yōu)化時,同樣根據(jù)常見材料的模量值,給定彈性模量范圍為4×109~3.45×1011Pa。其中4×109為所選鉛材料彈性模量的1/4,小于基體的彈性模量,3.45×1011與材料參數(shù)表中常用材料鎢的彈性模量相同,散射體的密度與鉛一致,優(yōu)化結(jié)果如圖3所示。

    相比于文獻(xiàn)[7]、[23],我們得到了截然不同的構(gòu)型,且相對帶寬達(dá)到了1.07,不但遠(yuǎn)大于散射體為鉛的情況,甚至大于文獻(xiàn)[23]中以銅為散射體的情況(0.902)。這似乎表明彈性模量對聲子晶體帶隙有著更為重要的影響。事實上,在文獻(xiàn)[23]中,以銅為散射體的帶隙結(jié)果(0.902)也遠(yuǎn)大于以鉛為散射體(0.455)的情況,這應(yīng)該是由于銅的彈性模量遠(yuǎn)大于鉛所導(dǎo)致的。

    圖3 考慮彈性模量變化時的優(yōu)化結(jié)果

    再者,當(dāng)把算例1和算例2的結(jié)果同時與文獻(xiàn)[23]的結(jié)果比較時可以發(fā)現(xiàn),除禁帶位置降低外,當(dāng)僅以密度為變量優(yōu)化時,能帶圖第三階頻率下移,當(dāng)僅以彈性模量為變量優(yōu)化時,能帶圖第四階頻率上移。這表明密度對禁帶的下邊界影響較大,而彈性模量對禁帶上邊界影響較大,這也與文獻(xiàn)[24]中的結(jié)論一致。

    2.3 同時考慮密度和彈性模量對帶隙的影響

    以密度和彈性模量同時為變量時,兩種參數(shù)的范圍選擇與算例1、2相同,優(yōu)化后密度結(jié)果為19 320 kg/m3,彈性模量為3.45×1011Pa。優(yōu)化結(jié)果如圖4所示??梢钥闯觯?dāng)彈性模量保持較大的比值時,增大密度差能進(jìn)一步的提高相對帶寬,且可實現(xiàn)禁帶寬度增大和頻率位置降低兩個目標(biāo),并得到了遠(yuǎn)大于先前工作的相對帶隙,結(jié)果再一次證明了彈性模量在聲子晶體設(shè)計中的重要性。

    圖4 密度和彈性模量同時變化時的優(yōu)化結(jié)果

    3 結(jié)論

    本文應(yīng)用遺傳算法對聲子晶體拓?fù)錁?gòu)型及影響聲子晶體帶隙特性的密度和彈性模量兩種材料參數(shù)進(jìn)行優(yōu)化。通過引入混合變量,分別考慮了密度變化、彈性模量變化以及同時考慮兩種參數(shù)作為連續(xù)變量對聲子晶體帶隙的影響,結(jié)果表明兩種參數(shù)對帶隙特性都有顯著影響,其中彈性模量具有更重要的地位。優(yōu)化后密度和彈性模量都接近甚至達(dá)到給定范圍的上限,這表明散射體的密度和彈性模量與基體材料對應(yīng)參數(shù)的差值越大越易產(chǎn)生較寬的帶隙,其中密度對帶隙的下邊界影響較大,彈性模量對帶隙的上邊界影響較大。通過對比3個算例可以發(fā)現(xiàn),當(dāng)彈性模量具有更大的差值時,更易產(chǎn)生較寬的帶隙,這為設(shè)計聲子晶體時的材料選擇和帶隙調(diào)控提供了一定的理論依據(jù)。

    猜你喜歡
    散射體帶隙聲子
    半無限板類聲子晶體帶隙仿真的PWE/NS-FEM方法
    一種基于單次散射體定位的TOA/AOA混合定位算法*
    密度泛函理論計算半導(dǎo)體材料的帶隙誤差研究
    納米表面聲子 首次實現(xiàn)三維成像
    聲子晶體覆蓋層吸聲機理研究
    二維結(jié)構(gòu)中亞波長缺陷的超聲特征
    無損檢測(2019年11期)2019-11-20 07:07:50
    一種基于BJT工藝的無運放低溫度系數(shù)的帶隙基準(zhǔn)源
    間距比對雙振子局域共振軸縱振帶隙的影響
    一款高PSRR低溫度系數(shù)的帶隙基準(zhǔn)電壓源的設(shè)計
    電子制作(2018年1期)2018-04-04 01:48:38
    高斯波包散射體成像方法
    老司机深夜福利视频在线观看| 国产成人av激情在线播放| 精品无人区乱码1区二区| 中文亚洲av片在线观看爽| 在线观看午夜福利视频| 精品一区二区三区视频在线观看免费| 一卡2卡三卡四卡精品乱码亚洲| 男女下面进入的视频免费午夜 | 高清在线国产一区| 午夜亚洲福利在线播放| 中文字幕精品免费在线观看视频| 国产激情久久老熟女| 美女 人体艺术 gogo| 久久天躁狠狠躁夜夜2o2o| 亚洲av熟女| www.999成人在线观看| 国产亚洲精品久久久久久毛片| 成年女人毛片免费观看观看9| www国产在线视频色| 两性夫妻黄色片| 日本 欧美在线| 午夜亚洲福利在线播放| a级毛片在线看网站| 欧美日韩瑟瑟在线播放| 精品国产乱子伦一区二区三区| 色综合亚洲欧美另类图片| 亚洲国产精品sss在线观看| 亚洲精品在线观看二区| 黄网站色视频无遮挡免费观看| 又黄又粗又硬又大视频| 免费无遮挡裸体视频| 成年版毛片免费区| 久久久久久久久免费视频了| 9色porny在线观看| 丁香六月欧美| 亚洲专区字幕在线| 搞女人的毛片| 国产激情欧美一区二区| 欧美日韩瑟瑟在线播放| 久久国产精品影院| 无限看片的www在线观看| 色婷婷久久久亚洲欧美| 中文亚洲av片在线观看爽| 日韩精品青青久久久久久| 久久国产乱子伦精品免费另类| 国产男靠女视频免费网站| 人人妻,人人澡人人爽秒播| 国产主播在线观看一区二区| 亚洲精品一区av在线观看| 男人舔女人下体高潮全视频| 久久天堂一区二区三区四区| 国产精品香港三级国产av潘金莲| 大型黄色视频在线免费观看| 此物有八面人人有两片| 日日爽夜夜爽网站| 欧美最黄视频在线播放免费| 丝袜美足系列| 99国产精品99久久久久| 美女国产高潮福利片在线看| 免费无遮挡裸体视频| 国产一级毛片七仙女欲春2 | 欧美日本视频| 夜夜夜夜夜久久久久| 夜夜夜夜夜久久久久| 成熟少妇高潮喷水视频| 国产精品久久电影中文字幕| 老鸭窝网址在线观看| 日韩欧美免费精品| 亚洲欧美一区二区三区黑人| 午夜免费激情av| 午夜老司机福利片| 国产精品久久电影中文字幕| 日日爽夜夜爽网站| 国产精品1区2区在线观看.| 亚洲 欧美 日韩 在线 免费| 欧美日本亚洲视频在线播放| 国产高清激情床上av| 午夜福利成人在线免费观看| 亚洲欧美精品综合久久99| 成人18禁在线播放| www.www免费av| 午夜免费激情av| 久久中文字幕人妻熟女| 亚洲精品粉嫩美女一区| 久久婷婷人人爽人人干人人爱 | 亚洲人成77777在线视频| 国产精品爽爽va在线观看网站 | 国产亚洲精品第一综合不卡| 欧美av亚洲av综合av国产av| 美女午夜性视频免费| 国语自产精品视频在线第100页| 精品不卡国产一区二区三区| 日本免费一区二区三区高清不卡 | 极品人妻少妇av视频| 亚洲欧美激情综合另类| 国产精品久久电影中文字幕| 欧美黄色淫秽网站| 亚洲第一电影网av| 怎么达到女性高潮| 欧美不卡视频在线免费观看 | 嫩草影视91久久| 波多野结衣av一区二区av| 亚洲国产看品久久| 亚洲 国产 在线| 亚洲精品久久成人aⅴ小说| 99精品欧美一区二区三区四区| 国产又色又爽无遮挡免费看| 一区福利在线观看| 日韩视频一区二区在线观看| 亚洲国产精品成人综合色| 黄片大片在线免费观看| 可以免费在线观看a视频的电影网站| 国产成人啪精品午夜网站| 国产99久久九九免费精品| 啦啦啦观看免费观看视频高清 | 女性生殖器流出的白浆| 午夜福利18| 亚洲电影在线观看av| 国产精品影院久久| 他把我摸到了高潮在线观看| 91麻豆精品激情在线观看国产| 久久久水蜜桃国产精品网| 18禁美女被吸乳视频| 久久亚洲真实| 精品国产一区二区久久| 精品久久蜜臀av无| √禁漫天堂资源中文www| 超碰成人久久| 动漫黄色视频在线观看| 18禁观看日本| 1024香蕉在线观看| 制服丝袜大香蕉在线| 亚洲人成伊人成综合网2020| 日本vs欧美在线观看视频| 欧美中文综合在线视频| 欧美日本视频| 亚洲中文字幕一区二区三区有码在线看 | 欧美成人一区二区免费高清观看 | 午夜亚洲福利在线播放| 美女高潮到喷水免费观看| 一个人观看的视频www高清免费观看 | 99re在线观看精品视频| 人人妻人人爽人人添夜夜欢视频| 久久婷婷人人爽人人干人人爱 | 在线国产一区二区在线| 亚洲最大成人中文| 咕卡用的链子| 天天一区二区日本电影三级 | 精品高清国产在线一区| 国产成人精品久久二区二区免费| 国产在线精品亚洲第一网站| 亚洲国产精品合色在线| 99久久久亚洲精品蜜臀av| 亚洲中文字幕日韩| 91成人精品电影| 国产一区二区在线av高清观看| 久热这里只有精品99| 天堂√8在线中文| 亚洲av熟女| 精品无人区乱码1区二区| 级片在线观看| 久久精品亚洲精品国产色婷小说| 久久 成人 亚洲| 又紧又爽又黄一区二区| 久久欧美精品欧美久久欧美| 国产欧美日韩一区二区三| 在线免费观看的www视频| 国产精品亚洲美女久久久| 色婷婷久久久亚洲欧美| 亚洲 欧美 日韩 在线 免费| av视频免费观看在线观看| 国产一区二区三区在线臀色熟女| 午夜免费激情av| 国产亚洲欧美精品永久| 色综合婷婷激情| 不卡一级毛片| 国产精品,欧美在线| 亚洲aⅴ乱码一区二区在线播放 | 一a级毛片在线观看| 亚洲国产毛片av蜜桃av| 午夜精品在线福利| 欧美乱码精品一区二区三区| 美女高潮喷水抽搐中文字幕| 精品欧美一区二区三区在线| 18禁观看日本| 嫁个100分男人电影在线观看| 啦啦啦韩国在线观看视频| 人人妻人人澡欧美一区二区 | 亚洲一区二区三区色噜噜| 亚洲熟妇熟女久久| 国产精品久久电影中文字幕| 免费看a级黄色片| 亚洲国产欧美一区二区综合| 久久精品国产综合久久久| 国产免费av片在线观看野外av| 禁无遮挡网站| 老司机深夜福利视频在线观看| 一a级毛片在线观看| 欧美色视频一区免费| 一本综合久久免费| 亚洲av美国av| 欧美黄色淫秽网站| av免费在线观看网站| 欧美日韩中文字幕国产精品一区二区三区 | 丁香六月欧美| 国产私拍福利视频在线观看| 国产成人欧美在线观看| 亚洲精品在线美女| 最近最新免费中文字幕在线| 精品国产国语对白av| 一进一出抽搐gif免费好疼| 亚洲av成人av| 婷婷丁香在线五月| 欧美日本亚洲视频在线播放| 一本大道久久a久久精品| 国产欧美日韩综合在线一区二区| 一本久久中文字幕| 老司机午夜十八禁免费视频| 可以在线观看毛片的网站| 九色亚洲精品在线播放| 日韩 欧美 亚洲 中文字幕| 男女之事视频高清在线观看| a在线观看视频网站| 黄片播放在线免费| 97人妻精品一区二区三区麻豆 | 啦啦啦 在线观看视频| 自拍欧美九色日韩亚洲蝌蚪91| 97人妻天天添夜夜摸| 国产一区二区三区综合在线观看| 熟女少妇亚洲综合色aaa.| 女人被躁到高潮嗷嗷叫费观| av视频在线观看入口| 免费观看精品视频网站| 少妇被粗大的猛进出69影院| 免费在线观看黄色视频的| 亚洲色图av天堂| 亚洲成人国产一区在线观看| 真人一进一出gif抽搐免费| 99久久99久久久精品蜜桃| 波多野结衣av一区二区av| 久久中文字幕一级| 看黄色毛片网站| 他把我摸到了高潮在线观看| 亚洲avbb在线观看| 成人三级黄色视频| 人人妻人人澡人人看| 精品不卡国产一区二区三区| 亚洲国产精品成人综合色| 国产成人精品久久二区二区免费| 午夜免费观看网址| 亚洲av电影不卡..在线观看| 村上凉子中文字幕在线| 国产极品粉嫩免费观看在线| 黑人巨大精品欧美一区二区mp4| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲国产日韩欧美精品在线观看 | 亚洲黑人精品在线| 女人爽到高潮嗷嗷叫在线视频| 91麻豆精品激情在线观看国产| 国产麻豆69| 又紧又爽又黄一区二区| 色播亚洲综合网| 91在线观看av| 亚洲无线在线观看| 99精品在免费线老司机午夜| 国产91精品成人一区二区三区| 97超级碰碰碰精品色视频在线观看| 亚洲欧美日韩另类电影网站| 日日干狠狠操夜夜爽| 黄色毛片三级朝国网站| 99国产精品免费福利视频| 美女扒开内裤让男人捅视频| 欧美av亚洲av综合av国产av| 操美女的视频在线观看| 天天躁夜夜躁狠狠躁躁| 亚洲成av人片免费观看| 国产精品一区二区精品视频观看| 99国产综合亚洲精品| 啦啦啦观看免费观看视频高清 | 一级毛片女人18水好多| 久久中文看片网| 久久精品国产99精品国产亚洲性色 | 国产一区二区三区视频了| 90打野战视频偷拍视频| 成人三级做爰电影| 咕卡用的链子| 国产亚洲精品第一综合不卡| 无人区码免费观看不卡| 国产精品综合久久久久久久免费 | av电影中文网址| 99国产精品免费福利视频| 午夜成年电影在线免费观看| 国产一区在线观看成人免费| 黄色片一级片一级黄色片| 欧美激情久久久久久爽电影 | 日韩有码中文字幕| 久久久久久久久中文| 午夜福利18| 中文字幕人成人乱码亚洲影| or卡值多少钱| 免费在线观看亚洲国产| 亚洲avbb在线观看| 九色国产91popny在线| 久久精品91蜜桃| 手机成人av网站| 麻豆久久精品国产亚洲av| av片东京热男人的天堂| 淫妇啪啪啪对白视频| 成年女人毛片免费观看观看9| 精品国产一区二区久久| 精品国产国语对白av| 巨乳人妻的诱惑在线观看| 亚洲精华国产精华精| 欧美中文综合在线视频| 激情视频va一区二区三区| 久久久久久大精品| 久久精品国产99精品国产亚洲性色 | 丰满人妻熟妇乱又伦精品不卡| 九色亚洲精品在线播放| 天堂动漫精品| 99久久精品国产亚洲精品| 色老头精品视频在线观看| 中文字幕人妻丝袜一区二区| 国产av又大| 中文字幕精品免费在线观看视频| av在线播放免费不卡| 中文字幕人妻熟女乱码| 日韩国内少妇激情av| 日韩欧美三级三区| 日本vs欧美在线观看视频| av电影中文网址| 99re在线观看精品视频| 免费女性裸体啪啪无遮挡网站| 麻豆国产av国片精品| 成人精品一区二区免费| 亚洲av成人不卡在线观看播放网| 久久中文字幕人妻熟女| 狂野欧美激情性xxxx| 一二三四在线观看免费中文在| 国产欧美日韩一区二区三| 女警被强在线播放| 老司机深夜福利视频在线观看| av在线播放免费不卡| 精品欧美国产一区二区三| 又黄又粗又硬又大视频| 国产精品久久久人人做人人爽| 成人手机av| 国产亚洲精品综合一区在线观看 | 日日夜夜操网爽| 97碰自拍视频| 嫩草影院精品99| 天堂影院成人在线观看| 亚洲中文字幕一区二区三区有码在线看 | 美女高潮到喷水免费观看| 国产精品永久免费网站| 午夜免费成人在线视频| 人成视频在线观看免费观看| 99国产精品免费福利视频| 国产精品美女特级片免费视频播放器 | 99久久国产精品久久久| 在线十欧美十亚洲十日本专区| 男女做爰动态图高潮gif福利片 | 一级毛片精品| 99re在线观看精品视频| 黄色片一级片一级黄色片| 亚洲国产高清在线一区二区三 | 婷婷六月久久综合丁香| 十八禁人妻一区二区| 亚洲最大成人中文| 亚洲久久久国产精品| 亚洲欧美精品综合久久99| 午夜福利高清视频| 久久天躁狠狠躁夜夜2o2o| 亚洲专区字幕在线| 天堂√8在线中文| 精品日产1卡2卡| 欧美+亚洲+日韩+国产| 日韩大码丰满熟妇| 亚洲va日本ⅴa欧美va伊人久久| 国产av一区二区精品久久| 精品一区二区三区四区五区乱码| 国产精品98久久久久久宅男小说| 黄色视频不卡| 婷婷丁香在线五月| 黑人巨大精品欧美一区二区蜜桃| 国产成人精品无人区| 99久久99久久久精品蜜桃| 在线天堂中文资源库| 日本在线视频免费播放| 国产精品久久久av美女十八| 操出白浆在线播放| 精品国产亚洲在线| 精品免费久久久久久久清纯| 一级毛片精品| 欧美黑人欧美精品刺激| 大码成人一级视频| 99久久久亚洲精品蜜臀av| 中国美女看黄片| 美女国产高潮福利片在线看| 欧美日韩亚洲综合一区二区三区_| 正在播放国产对白刺激| 免费少妇av软件| 操出白浆在线播放| 一边摸一边做爽爽视频免费| 大型av网站在线播放| 成人永久免费在线观看视频| 久久久精品欧美日韩精品| 久久 成人 亚洲| 国产一区二区三区综合在线观看| 美女扒开内裤让男人捅视频| 欧美av亚洲av综合av国产av| 久久久精品国产亚洲av高清涩受| 手机成人av网站| 中文字幕另类日韩欧美亚洲嫩草| 两个人看的免费小视频| 亚洲国产毛片av蜜桃av| 在线观看免费午夜福利视频| 在线观看www视频免费| 国产精品久久久久久人妻精品电影| 不卡av一区二区三区| 久久国产精品人妻蜜桃| 午夜影院日韩av| or卡值多少钱| 欧美久久黑人一区二区| 大型黄色视频在线免费观看| 免费搜索国产男女视频| 欧美最黄视频在线播放免费| 国产又爽黄色视频| 在线观看免费日韩欧美大片| 不卡av一区二区三区| 88av欧美| 久久久久久亚洲精品国产蜜桃av| 亚洲人成网站在线播放欧美日韩| 嫩草影院精品99| 久久青草综合色| 久久天堂一区二区三区四区| 大码成人一级视频| av中文乱码字幕在线| 18禁黄网站禁片午夜丰满| 久久久久国产精品人妻aⅴ院| 久久精品国产亚洲av香蕉五月| 无人区码免费观看不卡| 亚洲中文字幕日韩| АⅤ资源中文在线天堂| 国产区一区二久久| 精品福利观看| 99国产综合亚洲精品| 看免费av毛片| 丁香欧美五月| 嫁个100分男人电影在线观看| 99国产精品99久久久久| 久久国产精品男人的天堂亚洲| 99riav亚洲国产免费| 午夜福利影视在线免费观看| 亚洲精品美女久久av网站| 欧美成人一区二区免费高清观看 | 亚洲国产欧美网| 精品国产美女av久久久久小说| 国产麻豆69| 少妇粗大呻吟视频| a级毛片在线看网站| netflix在线观看网站| 日韩欧美三级三区| 午夜久久久在线观看| 欧美日本亚洲视频在线播放| 亚洲一卡2卡3卡4卡5卡精品中文| 免费观看精品视频网站| 日韩成人在线观看一区二区三区| 国产精品日韩av在线免费观看 | 欧美不卡视频在线免费观看 | 他把我摸到了高潮在线观看| 国产精品美女特级片免费视频播放器 | 亚洲欧美激情在线| 最近最新免费中文字幕在线| 午夜福利影视在线免费观看| 一个人免费在线观看的高清视频| 一级毛片女人18水好多| 少妇被粗大的猛进出69影院| 欧美成人免费av一区二区三区| 午夜精品国产一区二区电影| 午夜福利在线观看吧| 久久国产亚洲av麻豆专区| 中文字幕人成人乱码亚洲影| 久久精品国产99精品国产亚洲性色 | 在线播放国产精品三级| 日韩 欧美 亚洲 中文字幕| 精品久久蜜臀av无| 乱人伦中国视频| 精品欧美国产一区二区三| 亚洲成人免费电影在线观看| 电影成人av| 亚洲午夜精品一区,二区,三区| 免费高清在线观看日韩| 91精品国产国语对白视频| 亚洲无线在线观看| 国产熟女xx| 欧美日韩乱码在线| 老司机深夜福利视频在线观看| 91成年电影在线观看| 在线观看舔阴道视频| 俄罗斯特黄特色一大片| 欧美性长视频在线观看| 男女下面插进去视频免费观看| 国产免费男女视频| 亚洲精品中文字幕一二三四区| 精品久久久精品久久久| 免费高清在线观看日韩| 国产亚洲精品av在线| 91av网站免费观看| 欧美乱色亚洲激情| 国产一区二区三区在线臀色熟女| 精品第一国产精品| 99精品欧美一区二区三区四区| xxx96com| 91精品三级在线观看| 91老司机精品| 天堂√8在线中文| 国产成人免费无遮挡视频| 国产三级黄色录像| 国产精品野战在线观看| 女生性感内裤真人,穿戴方法视频| 亚洲欧洲精品一区二区精品久久久| 日本五十路高清| 99在线视频只有这里精品首页| 97人妻天天添夜夜摸| 免费女性裸体啪啪无遮挡网站| 国产激情欧美一区二区| 日韩三级视频一区二区三区| 看片在线看免费视频| 99在线人妻在线中文字幕| 国产精品一区二区免费欧美| 精品人妻在线不人妻| 日韩一卡2卡3卡4卡2021年| 老汉色∧v一级毛片| 最近最新免费中文字幕在线| 国产精品1区2区在线观看.| 最近最新免费中文字幕在线| 97人妻天天添夜夜摸| 国产精品久久久久久亚洲av鲁大| 麻豆成人av在线观看| 日本精品一区二区三区蜜桃| 免费久久久久久久精品成人欧美视频| 动漫黄色视频在线观看| 嫁个100分男人电影在线观看| 久久精品亚洲熟妇少妇任你| 精品乱码久久久久久99久播| 亚洲国产中文字幕在线视频| 亚洲国产精品999在线| av片东京热男人的天堂| 男女午夜视频在线观看| 操美女的视频在线观看| 此物有八面人人有两片| 91国产中文字幕| 久久久久久免费高清国产稀缺| 69av精品久久久久久| 亚洲va日本ⅴa欧美va伊人久久| 99久久国产精品久久久| av免费在线观看网站| 亚洲欧美一区二区三区黑人| 手机成人av网站| 亚洲人成电影观看| 精品国产亚洲在线| tocl精华| 每晚都被弄得嗷嗷叫到高潮| 桃红色精品国产亚洲av| 美女国产高潮福利片在线看| 男女之事视频高清在线观看| 国产伦一二天堂av在线观看| 国产麻豆69| 天天躁夜夜躁狠狠躁躁| 日本vs欧美在线观看视频| 久久亚洲精品不卡| 久久人人爽av亚洲精品天堂| videosex国产| 无限看片的www在线观看| 男女下面插进去视频免费观看| av天堂久久9| 精品国产美女av久久久久小说| 老司机午夜福利在线观看视频| 波多野结衣高清无吗| 国产麻豆成人av免费视频| 亚洲av成人不卡在线观看播放网| 亚洲人成电影免费在线| 男人舔女人的私密视频| 亚洲人成伊人成综合网2020| 久久欧美精品欧美久久欧美| 欧美日韩福利视频一区二区| 亚洲一区二区三区色噜噜| 亚洲一卡2卡3卡4卡5卡精品中文| 国产99白浆流出| 国内精品久久久久久久电影| 啦啦啦免费观看视频1| 欧美另类亚洲清纯唯美| 欧美在线黄色| www.熟女人妻精品国产| 99精品欧美一区二区三区四区| 久久热在线av| 国产一区二区三区在线臀色熟女| 免费在线观看视频国产中文字幕亚洲| 久久中文字幕人妻熟女| 97人妻天天添夜夜摸| 后天国语完整版免费观看| 亚洲av成人不卡在线观看播放网| 亚洲人成伊人成综合网2020| 国产又色又爽无遮挡免费看| 精品久久久久久久久久免费视频| 国产一区二区在线av高清观看| 91在线观看av| 亚洲精华国产精华精| 亚洲精品国产区一区二| 国产成人一区二区三区免费视频网站| 日本免费一区二区三区高清不卡 | 色综合婷婷激情| 黄色视频,在线免费观看| 手机成人av网站|