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

    材料粘彈性對于一維局域共振聲子晶體帶隙的影響*

    2016-11-12 05:21:26左曙光韋錫晉倪天心吳旭東
    功能材料 2016年10期
    關(guān)鍵詞:粘彈性帶隙聲子

    左曙光,韋錫晉,倪天心,吳旭東

    (同濟大學 新能源汽車工程中心, 上海 201804)

    ?

    材料粘彈性對于一維局域共振聲子晶體帶隙的影響*

    左曙光,韋錫晉,倪天心,吳旭東

    (同濟大學 新能源汽車工程中心, 上海 201804)

    聲子晶體是具有聲子帶隙的人工周期性復合結(jié)構(gòu)。作為常見的聲子晶體的組成材料,橡膠材料的粘彈特性很少被考慮到,但其會對聲子晶體帶隙特性有重要影響。研究聲子晶體考慮橡膠的粘彈性后的帶隙特性,目的在于研究聲子晶體在實際應用中的特性,以更好的應用聲子晶體。建立了一維聲子晶體梁的有限元模型。使用微分型粘彈性本構(gòu)模型中的三參數(shù)Maxwell模型,推導了它的松弛剪切模量,將模型中的粘彈性轉(zhuǎn)化為有限元模型中的材料參數(shù)。并通過進行有限元計算,比較和分析了考慮粘彈性之后對于不同形式聲子晶體帶隙的影響,發(fā)現(xiàn)粘彈性會影響帶隙的位置和寬度,同時減小了振動的衰減。

    局域共振;聲子晶體;帶隙;粘彈性;有限元

    0 引 言

    從20世紀90年代開始,對聲子晶體有了比較廣泛的研究[1-3]。在2000年以前,所研究的聲子晶體的彈性波帶隙頻率對應的彈性波波長同聲子晶體晶格尺寸相當,即彈性波帶隙的形成遵循Bragg散射機理。2000年,劉正猷教授發(fā)現(xiàn)使用軟材料包覆的鉛球埋在環(huán)氧樹脂中形成的聲子晶體的帶隙頻率對應的波長遠遠大于晶格尺寸,由此提出了聲子晶體的局域共振帶隙機理[4]。

    在實際應用中,材料具有很多復雜的屬性,比如經(jīng)常被用到作包覆層材料的橡膠具有粘彈性,會對帶隙特性產(chǎn)生一定影響。粘彈材料是同時具有流體粘性和固體彈性的材料,如塑料、橡膠等聚合物等。對于粘彈性材料,彈性模量不僅為復數(shù),而且是頻率的函數(shù)。將粘彈材料引入聲子晶體的材料中,會帶來一些特殊的影響,引起了學者的關(guān)注。在處理粘彈材料時,有的學者將其視為具有特定阻尼的彈性材料,沒有考慮阻尼和彈性常數(shù)隨頻率的變化[5-6];有的學者用復數(shù)來表示彈性模量,并認為其虛部也就是損耗模量表征為對能量的耗散,而決定彈性波傳播行為的是其實部也就是儲能模量,因此只考慮儲能模量對聲子晶體色散關(guān)系的影響[7]。王剛、溫激鴻等還對比了理論計算和實驗測試的結(jié)果,認為粘彈阻尼不改變結(jié)構(gòu)的色散關(guān)系,只表現(xiàn)為對高頻能量的耗散作用,從而增大了帶隙的寬度[8]。上述研究中都認為損耗模量只表現(xiàn)為對能量的耗散而對色散關(guān)系沒有影響。Merheb等同時考慮了材料粘彈性對聲子晶體結(jié)構(gòu)色散關(guān)系和耗散的影響[9]。這些人的研究都只針對了簡單的布拉格形式的聲子晶體,并沒有研究材料粘彈性對于局域共振聲子晶體帶隙特性的影響,這兩者因為機理不同,粘彈性的影響也會有所不同。基于此,本文主要針對粘彈性材料屬性對于局域共振聲子晶體帶隙特性的影響進行研究。

    目前的理論算法都無法考慮到粘彈性材料的彈性模量隨頻率變化的特性。所以本文使用有限元算法進行計算。本文首先建立了一維聲子晶體的有限元模型,通過引入三參數(shù)微分型粘彈本構(gòu)模型來模擬橡膠材料的粘彈性,通過推導將理論模型中的參數(shù)轉(zhuǎn)化為有限元模型中的材料參數(shù),給有限元模型賦予了粘彈屬性,最后使用有限元方法研究材料粘彈性對聲子晶體能帶結(jié)構(gòu)的影響,分析和討論了損耗模量和儲能模量對于聲子晶體帶隙特性的影響。

    1 一維聲子晶體有限元模型的建立

    如圖1所示,為本文所使用的典型一維聲子晶體梁結(jié)構(gòu)的示意圖[10]。梁使用的是有機玻璃,界面尺寸為b=0.01 m和h=0.005 m。振子由橡膠(長0.01 m、寬0.01 m、高0.002 m)粘結(jié)銅塊構(gòu)成,晶格常數(shù)為0.08 m,銅塊質(zhì)量為0.013 kg。材料參數(shù)如表1所示。

    因為該結(jié)構(gòu)簡單,組成部分俱為規(guī)則的立方體單元,選擇使用六面體實體單元。在Abaqus中典型的六面體實體單元主要有線性單元和二次單元,它們又分為完全積分和縮減積分兩種。因為分析的振動主要是彎曲振動,使用完全積分單元會導致剪力自鎖,使得單元在彎曲時剛度過大,所以選擇縮減積分單元。為糾正使用線性縮減積分單元(C3DR8)時可能產(chǎn)生的剛度“泄露”現(xiàn)象,仿真中同時使用了二次精度控制和沙漏控制。

    圖1 一維局域共振聲子晶體示意圖

    表1 材料屬性

    該聲子晶體結(jié)構(gòu)由3種材料組成,其接觸面在實際中需要用膠合物連接。在本文中是考慮材料的粘彈性,所以對膠合物的性質(zhì)不作討論,在有限元軟件中使用共節(jié)點的方式進行連接。即3種材料在接觸面處不會發(fā)生滑移。繪制完網(wǎng)格后使用Abaqus中的check element功能,用網(wǎng)格大小角度等參數(shù)來判斷繪制的網(wǎng)格是否有比較高的質(zhì)量。

    給模型賦予材料參數(shù)后即可開始進行有限元計算。給梁的一端加上垂直于板面幅值為1 mm/s2的簡諧加速度激勵,進行掃頻分析,在梁的另一端讀取加速度響應即為傳輸頻響函數(shù)。加速度出現(xiàn)衰減的頻率范圍即為所求的帶隙。綜合考慮計算時間和所需求的精度,選擇計算步長為5 Hz。在理論計算中所求的聲子晶體是無限周期的,而在有限結(jié)構(gòu)中,帶隙的衰減大小與位置和結(jié)構(gòu)的周期數(shù)有明顯的關(guān)系,為了能夠清晰得出帶隙的位置從而得到準確的結(jié)論必須選擇合適的周期數(shù)。圖2所示為不同周期數(shù)的聲子晶體的局域共振帶隙的計算值??梢娭芷跀?shù)達到5個之后計算得到的帶隙值已經(jīng)基本一致,鑒于周期數(shù)越多計算時間越久,選擇5個周期的有限元模型進行之后的計算。有限元模型如圖3所示。

    2 粘彈性材料模型及其向有限元軟件的轉(zhuǎn)化

    粘彈性是材料的一種本構(gòu)性質(zhì),粘彈性固體材料在性質(zhì)上介于彈性體和粘性流體之間:它對突加載荷像彈性體一樣產(chǎn)生瞬時響應;對動載荷響應又像粘性流體那樣滯后,而滯后的時間介于彈性體和粘性流體之間。組成聲子晶體結(jié)構(gòu)中橡膠材料的粘彈性對于帶隙性質(zhì)有很大的影響,所以引進粘彈性模型來對材料粘彈性的影響進行研究。

    圖2 不同周期數(shù)有限元計算得到的帶隙值

    圖3 一維局域共振聲子晶體模型

    微分型粘彈性本構(gòu)關(guān)系是基于粘彈性材料的機械模型(彈簧和粘壺組成的模型)所建立的本構(gòu)方程。目前常見幾種粘彈性模型如圖4所示[12]。圖4(a)的K-V模型能表現(xiàn)粘彈性材料的蠕變特性,不能表現(xiàn)其松弛特性。圖4(b)的二參數(shù)Maxwell模型能表現(xiàn)松弛特性,不能表現(xiàn)蠕變特性。而對于三參數(shù)Maxwell模型和Maxwell并聯(lián)模型,既能描述橡膠的應力松弛,又能描述橡膠的蠕變特性??紤]到本文是為了探討粘彈性對于聲子晶體帶隙特性的影響,而不是需要并聯(lián)模型的多個參數(shù)能夠很精確的描述實際橡膠的運動特性,所以在研究中選擇三參數(shù)Maxwell模型。圖4(c)的三參數(shù)Maxwell由3個部分組成,粘性系數(shù)為η的粘壺,和粘壺串聯(lián)的彈性模量為E2的彈簧,以及和他們并聯(lián)的彈性模量為E1的彈簧。

    圖4 K-V模型、二參數(shù)Maxwell模型、三參數(shù)Maxwell模型和為Maxwell并聯(lián)模型

    橡膠材料粘彈特性可以分別在時域和頻域中描述。ABAQUS時域中提供的粘彈本構(gòu)模型是基于松弛剪切模量Prony級數(shù)展開形式的積分型本構(gòu)方程。可以通過Laplace變換,將線性粘彈性微分本構(gòu)關(guān)系轉(zhuǎn)化為Prony級數(shù)的形式。松弛剪切模量表示的是一個應力與應變速度的關(guān)系。對于三參數(shù)Maxwell模型,其線性粘彈微分本構(gòu)方程可以表示為

    (1)

    三參數(shù)Maxwell模型的松弛剪切模量為

    (2)

    再將其做Laplace逆變換,得到時域松弛剪切模量

    (3)

    模型突然加載,粘壺來不及反應,力經(jīng)過粘壺沒有幅值和相位的改變。所以瞬時剪切模量

    (4)

    把剪切模量去量綱化

    (5)

    將式(5)展開成Prony級數(shù)的形式

    (6)

    在Abaqus中,時域和頻域的參數(shù)值可以自動轉(zhuǎn)化,所以根據(jù)上述分析即可得到Abaqus中需要輸入的材料參數(shù)

    (7)

    (8)

    對于粘彈性材料,在研究其動態(tài)問題時,往往根據(jù)類比彈簧的彈性模量引入復模量模型。

    三參數(shù)Maxwell模型的微分本構(gòu)模型,將其寫成拉氏域動態(tài)模量為

    (9)

    將拉氏域變換到頻域可得

    (10)

    其中

    表示單元的相對粘度。

    (11)

    式(11)為復模量的實部,為儲能模量。

    (12)

    式(12)為復模量的虛部,為損失模量。

    3 粘彈性對于一維局域共振聲子晶體帶隙特性的影響

    在分析材料粘彈性影響之前,先對與之類似的材料阻尼特性進行考慮。阻尼作為粘彈性的一種特殊情況,采用復阻尼形式可以表示為

    μ為阻尼因子。如圖5所示為不同阻尼因子的計算所得的虛波矢圖。虛波矢即為帶隙位置,虛波矢越大,代表衰減越大??梢钥闯鱿鹉z的阻尼因子越大對帶隙的形成越不利,嚴重影響帶隙內(nèi)的振動衰減能力。當橡膠阻尼因子達到0.2時,帶隙的形狀已經(jīng)發(fā)生了很大的變化,難以對振動形成有效的衰減。

    圖5 不同阻尼因子帶隙對比

    考慮完了特殊情況,對粘彈性進行分析。如圖6所示,為考慮粘彈性的一維局域共振型聲子晶體和不考慮粘彈性的一維局域共振型聲子晶體彎曲振動頻率響應曲線的比較。計算中,材料的具體的數(shù)據(jù)與理論計算模型參數(shù)相同,粘彈性參數(shù)E2=0.5 MPa,η=0.0005。

    從圖6可以看出,布拉格帶隙的帶隙位置沒有因為橡膠的粘彈性而發(fā)生變化,這是因為布拉格帶隙的帶隙機理與局域共振帶隙不同,橡膠塊對于布拉格帶隙機理起到的只是質(zhì)量塊的作用。而對于局域共振帶隙可以看到帶隙起始頻率和帶隙截止頻率都出現(xiàn)了增大。這可以由橡膠復模量的實部隨著頻率增大解釋。橡膠復模量的實部隨著頻率增大而單調(diào)遞增,并且趨近于E1+E2。同時由于復模量實部的增大也使共振峰向高頻移動。這種特性對于擴大聲子晶體帶隙范圍有著明顯的效果。根據(jù)文獻[12]局域共振帶隙具有簡化模型,根據(jù)簡化模型可以得到帶隙的起始頻率和截止頻率的計算公式,分別為

    f1和f2分別為起始帶隙和截止帶隙,k為橡膠剛度,m1為振子質(zhì)量,m2為單個單元基體質(zhì)量。由公式可以看出,起始帶隙和截止帶隙均會隨著剛度的增大而增大,但是剛度在頻率低時增加小,在頻率高時增加多,也就是可以在保持起始帶隙變化很小的情況下,使得截止帶隙有較大的增長,這也就擴大了帶隙的范圍。同時從圖6可以看出,曲線的形狀并沒有發(fā)生明顯改變,只是出現(xiàn)的一些共振峰值都向高頻移動,并且幅值減小了,并且這種現(xiàn)象越接近高頻越明顯。這些明顯是由粘彈性模型中的損耗模量引起的,但是這部分變化對帶隙機理不會起影響,對于有限周期的聲子晶體,損耗模量的存在會使得衰減振動的能力減弱,但是引起的帶隙外的振動也會減小。

    圖6 振子梁彎曲振動考慮粘彈性與不考慮粘彈性頻率響應函數(shù)的對比

    下面針對三參數(shù)Maxwell模型中的參數(shù)變化具體分析。E1屬于不變的量,是橡膠的長效模量,它的變化相當于增大了橡膠的剛度,不影響損耗模量,所以不做具體討論。

    如圖7所示為E2增大后對于布拉格帶隙的影響,中心頻率也就是帶隙位置基本沒有發(fā)生變化,與之前的分析相同。在所取的參數(shù)范圍,由式(11)和(12)可得,儲能模量和損耗模量是遞增的。因為E2增大使得儲能模量和損耗模量都明顯增加,這使得帶隙的范圍擴大了。由理論計算可以分析彈性模量變大的影響,如圖8所示??梢钥闯?,布拉格帶隙的變化與橡膠剛度的變化并沒有關(guān)系。這證明是復模量中的損耗模量造成了帶隙的擴大。對此的理論解釋,因為損耗模量的存在,使得原本不處于帶隙頻率范圍內(nèi)的頻率段也產(chǎn)生了衰減,這種衰減并不是因為振子的振動和基體中波的耦合作用而產(chǎn)生的,而是因為橡膠中的能量衰減所造成的。

    圖7 E2變化產(chǎn)生的第一帶隙各參數(shù)的變化

    圖8 不同剛度橡膠帶隙的對比

    而對于第二帶隙也就是局域共振帶隙的影響就比較明顯。對于起始頻率影響比較小,對于截止頻率影響比較大,這也是式(11)所決定的。低頻影響小,高頻影響大。變化對于局域共振帶隙影響如圖9所示。起始頻率和截止頻率都會增大,截止頻率增加的更多,這就導致帶隙位置向高頻移動同時寬度增加。

    隨著η的增大儲能模量是單調(diào)遞增的,而損耗能量會先增大后減小。所以得到的帶隙的變化如圖10所示,起始頻率和截止頻率會隨著η的增加而增加,截止頻率增加的更多。損耗能量的變化在正常取值范圍不會明顯的體現(xiàn)出來。

    圖9 E2變化產(chǎn)生的第二帶隙各參數(shù)的變化

    圖10 η變化產(chǎn)生的第二帶隙各參數(shù)的變化

    4 結(jié) 論

    根據(jù)粘彈性材料的三參數(shù)Maxwell模型,利用有限元法對由粘彈性材料組成的一維聲子晶體的能帶結(jié)構(gòu)及帶隙特性進行了計算和理論分析。三參數(shù)模型中的3個參數(shù),除了不屬于粘彈性的考慮范疇長效彈性模量,其它兩個參數(shù)可以由上述分析得,它們都在不同程度上改變了粘彈性橡膠的儲能模量和損耗能量,并且他們的影響也可以通過儲能模量和損耗模量的改變來進行分析。

    由于材料粘彈特性的頻率相關(guān)性,儲能模量和損耗模量都是頻率的函數(shù),從而對聲子晶體的能帶結(jié)構(gòu)產(chǎn)生了重大影響。對比結(jié)果表明,對于局域共振類型的聲子晶體材料,其中作為彈簧作用的橡膠材料的粘彈性對于布拉格機理的帶隙基本不起作用,僅是損耗模量增大會對帶隙有微小的放大作用。而對于局域共振機理的帶隙,儲能模量隨頻率的變化會顯著的影響帶隙的位置和大小,而損耗模量會改變共振峰的位置,不影響帶隙的位置和大小,但會對衰減的幅度起到一定影響。

    材料粘彈性對聲子晶體能帶結(jié)構(gòu)的這些影響歸因于材料粘彈阻尼的頻率相關(guān)性。研究結(jié)果可以為粘彈性阻尼材料應用于聲子晶體提供一定的理論依據(jù),也為調(diào)節(jié)聲子晶體能帶結(jié)構(gòu)和獲得寬頻帶隙提供了一種新的方法。

    [1] Sigalas M M, Economou E N. Elastic and acoustic wave band structure[J]. Journal of Sound and Vibration, 1992, 158(2): 377-382.

    [2] Kushwaha M S, Halevi P, Dobrzynski L, et al. Acoustic band structure of periodic elastic composites[J]. Physical Review Letters, 1993, 71(13): 2022.

    [3] Sánchez-Pérez J V, Caballero D, Martinez-Sala R, et al. Sound attenuation by a two-dimensional array of rigid cylinders[J]. Physical Review Letters, 1998, 80(24): 5325.

    [4] Liu Z, Zhang X, Mao Y, et al. Locally resonant sonic materials[J]. Science, 2000, 289(5485): 1734-1736.

    [5] Psarobas I E. Viscoelastic response of sonic band-gap materials[J]. Physical Review B, 2001, 64(1): 012303.

    [6] Hussein M I. Theory of damped bloch waves in elastic media[J]. Physical Review B, 2009, 80(21): 212301.

    [7] Wang Gang, Wen Jihong, Liu Yaozong, et al. Vibration band gap of one-dimensional periodic structure with viscoelastic material[J]. Journal of Mechanical Engineering, 2004, 40(7): 47-50.

    王 剛, 溫激鴻, 劉耀宗, 等. 一維粘彈材料周期結(jié)構(gòu)的振動帶隙研究[J]. 機械工程學報, 2004, 40(7): 47-50.

    [8] Wen Jihong, Wang Gang, Liu Yaozong,et al. Vibration band gaps of metal/nitrile rubber rod-like phononic crystals structure[J]. Journal of Vibration Engineering, 2005, 18(1): 1-7.

    溫激鴻, 王 剛, 劉耀宗, 等. 金屬/丁腈橡膠桿狀結(jié)構(gòu)聲子晶體振動帶隙研究[J]. 振動工程學報, 2005, 18(1): 1-7.

    [9] Merheb B, Deymier P A, Jain M, et al. Elastic and viscoelastic effects in rubber/air acoustic band gap structures: A theoretical and experimental study[J]. Journal of Applied Physics, 2008, 104(6): 064913.

    [10] Wen Xisen, Wen Jihong, Yu Dianlong, et al. Phononic crystals[M]. Beijing: National Defense Industry Press, 2009: 2-15.

    溫熙森,溫激鴻,郁殿龍,等.聲子晶體[M].北京:國防工業(yè)出版社,2009:2-15.

    [11] Hu Jing. Analysis and optimization of EV suspension rubber bushing vibration characteristics[D]. Shanghai: Tongji University, 2014.

    胡 競. 電動車懸架橡膠襯套隔振特性分析與優(yōu)化設(shè)計[D]. 上海: 同濟大學, 2014.

    [12] Wang Gang. Research of band gap mechanism and damping characteristics of locally resonant phononic crystals[D]. Changsha: National Defense Scientific and Technical University, 2005.

    王 剛. 聲子晶體局域共振帶隙機理及減振特性研究[D]. 長沙: 國防科技大學, 2005.

    The effects of viscoelastic material on the band gap of one-dimensional locally resonant phononic crystal

    ZUO Shuguang, WEI Xijin, NI Tianxin, WU Xudong

    (Clean Energy Automotive Engineering Center, Tongji University, Shanghai 201804,China)

    Phononic crystals are artificial periodic composite structure with phononic band gap. As common phononic crystals constituent materials, the viscoelastic property of the rubber material is rarely taken into account, but it has a major impact on phononic crystal band gap characteristics. This paper studies the band gap properties of phononic crystal considering the viscoelasticity of the rubber materials for the purpose of studying phononic crystal features in the practical application so that we can have a better application of its vibration control characteristics. The shear modulus relaxation parameter is derived and transformed into the finite element model in Abaqus. The finite element model considering viscoelastic property is thus proposed and employed to analysis. According to the comparison and analysis of the influence on different forms of phononic band gaps after considering the viscoelasticity, it is found that the material viscoelasticity has the influence on the location and width of band gap and will weaken the attenuation ability of the band gap.

    local resonance; phononic crystal; band gap; viscoelasticity; finite element

    1001-9731(2016)10-10162-06

    國家自然科學基金青年基金資助項目(51305303)

    2015-09-10

    2016-04-22 通訊作者:吳旭東,E-mail: wuxudong@#edu.cn

    左曙光 (1968-),男,湖南沅江人,教授,博士生導師,主要從事新能源電動汽車振動與噪聲的研究。

    O328

    A

    10.3969/j.issn.1001-9731.2016.10.030

    猜你喜歡
    粘彈性帶隙聲子
    半無限板類聲子晶體帶隙仿真的PWE/NS-FEM方法
    二維粘彈性棒和板問題ADI有限差分法
    密度泛函理論計算半導體材料的帶隙誤差研究
    納米表面聲子 首次實現(xiàn)三維成像
    聲子晶體覆蓋層吸聲機理研究
    時變時滯粘彈性板方程的整體吸引子
    一種基于BJT工藝的無運放低溫度系數(shù)的帶隙基準源
    間距比對雙振子局域共振軸縱振帶隙的影響
    一款高PSRR低溫度系數(shù)的帶隙基準電壓源的設(shè)計
    電子制作(2018年1期)2018-04-04 01:48:38
    不可壓粘彈性流體的Leray-α-Oldroyd模型整體解的存在性
    又大又黄又爽视频免费| 国产91av在线免费观看| 国产白丝娇喘喷水9色精品| 91精品国产国语对白视频| 国产成人免费观看mmmm| 亚洲精品亚洲一区二区| 精品少妇黑人巨大在线播放| 午夜免费男女啪啪视频观看| 啦啦啦啦在线视频资源| 一区二区三区免费毛片| 九色成人免费人妻av| 一区二区三区精品91| 精品国产一区二区三区久久久樱花 | 激情 狠狠 欧美| 观看av在线不卡| 人妻夜夜爽99麻豆av| 国产一区二区三区综合在线观看 | 乱系列少妇在线播放| 日韩成人av中文字幕在线观看| 各种免费的搞黄视频| 国产乱人视频| 人体艺术视频欧美日本| 十八禁网站网址无遮挡 | 韩国av在线不卡| 18禁裸乳无遮挡免费网站照片| 久久ye,这里只有精品| 日韩中文字幕视频在线看片 | av女优亚洲男人天堂| 久久6这里有精品| 亚洲美女视频黄频| 亚洲精品日韩在线中文字幕| 成年美女黄网站色视频大全免费 | 久久精品国产自在天天线| 男的添女的下面高潮视频| 国产高清三级在线| 亚洲精品成人av观看孕妇| 精品久久久久久久末码| 搡老乐熟女国产| av又黄又爽大尺度在线免费看| 最黄视频免费看| 黄片wwwwww| 欧美精品人与动牲交sv欧美| 十八禁网站网址无遮挡 | 亚洲av综合色区一区| 久久久久网色| 国产精品熟女久久久久浪| 日韩成人av中文字幕在线观看| 自拍欧美九色日韩亚洲蝌蚪91 | 亚洲av国产av综合av卡| 欧美日韩综合久久久久久| 黄色一级大片看看| 在线播放无遮挡| 精品亚洲成国产av| 蜜桃亚洲精品一区二区三区| 亚洲欧美中文字幕日韩二区| 91久久精品国产一区二区成人| 国产精品精品国产色婷婷| 国产亚洲91精品色在线| 七月丁香在线播放| 啦啦啦在线观看免费高清www| 丰满人妻一区二区三区视频av| 久久久久久伊人网av| 成人黄色视频免费在线看| 这个男人来自地球电影免费观看 | av天堂中文字幕网| 免费看光身美女| 亚洲欧美成人精品一区二区| 中文字幕久久专区| 国产成人a区在线观看| 国产爽快片一区二区三区| 美女视频免费永久观看网站| 熟妇人妻不卡中文字幕| 晚上一个人看的免费电影| 成年免费大片在线观看| 日本欧美国产在线视频| 春色校园在线视频观看| 亚洲高清免费不卡视频| 卡戴珊不雅视频在线播放| 日韩成人av中文字幕在线观看| 欧美丝袜亚洲另类| 18禁在线播放成人免费| 777米奇影视久久| 国产精品人妻久久久影院| 欧美精品一区二区免费开放| 日韩av不卡免费在线播放| 肉色欧美久久久久久久蜜桃| 欧美最新免费一区二区三区| 成人亚洲精品一区在线观看 | 亚洲图色成人| 久久影院123| 搡女人真爽免费视频火全软件| 在线看a的网站| 人妻制服诱惑在线中文字幕| 国产伦在线观看视频一区| 男人添女人高潮全过程视频| 精品99又大又爽又粗少妇毛片| 亚洲人与动物交配视频| 舔av片在线| 色婷婷av一区二区三区视频| 国产极品天堂在线| 交换朋友夫妻互换小说| 夫妻性生交免费视频一级片| 国产精品国产av在线观看| 欧美日韩亚洲高清精品| 干丝袜人妻中文字幕| 久久99热这里只频精品6学生| 久久久国产一区二区| 午夜福利在线观看免费完整高清在| 国产伦在线观看视频一区| av专区在线播放| 99热这里只有是精品50| 如何舔出高潮| 99热6这里只有精品| 亚洲国产最新在线播放| 少妇裸体淫交视频免费看高清| 97超视频在线观看视频| 亚洲成人av在线免费| 国产大屁股一区二区在线视频| 少妇熟女欧美另类| 色婷婷久久久亚洲欧美| 免费大片18禁| 亚洲内射少妇av| 少妇被粗大猛烈的视频| 国产黄频视频在线观看| 亚洲国产高清在线一区二区三| 成人午夜精彩视频在线观看| 又爽又黄a免费视频| 18禁动态无遮挡网站| 国产无遮挡羞羞视频在线观看| av不卡在线播放| 欧美zozozo另类| 成人18禁高潮啪啪吃奶动态图 | 久久女婷五月综合色啪小说| 妹子高潮喷水视频| 亚洲av综合色区一区| 在线免费观看不下载黄p国产| 欧美成人一区二区免费高清观看| 久久久久久九九精品二区国产| 亚洲成色77777| 欧美区成人在线视频| 国产高清不卡午夜福利| 久久人人爽人人片av| 在现免费观看毛片| 久久韩国三级中文字幕| 人妻少妇偷人精品九色| 久久久久国产精品人妻一区二区| 亚洲欧美精品专区久久| 免费人成在线观看视频色| 国产淫片久久久久久久久| 一本一本综合久久| 你懂的网址亚洲精品在线观看| 日本欧美视频一区| 超碰97精品在线观看| 日韩,欧美,国产一区二区三区| 精品亚洲成a人片在线观看 | 国产一区有黄有色的免费视频| av在线app专区| 男人狂女人下面高潮的视频| 在线免费观看不下载黄p国产| 国产精品嫩草影院av在线观看| 久久97久久精品| av在线播放精品| av又黄又爽大尺度在线免费看| 中文天堂在线官网| 美女主播在线视频| 国产精品无大码| 欧美精品一区二区免费开放| 国产av码专区亚洲av| 国产成人午夜福利电影在线观看| 亚洲av中文字字幕乱码综合| 国产淫语在线视频| a级一级毛片免费在线观看| 国产欧美日韩精品一区二区| 中文字幕人妻熟人妻熟丝袜美| 新久久久久国产一级毛片| 日日啪夜夜爽| 一级黄片播放器| 国产精品99久久99久久久不卡 | 午夜免费观看性视频| 久久精品国产亚洲av天美| 嘟嘟电影网在线观看| 极品教师在线视频| 秋霞在线观看毛片| 美女高潮的动态| 狂野欧美激情性bbbbbb| 久久久久性生活片| 51国产日韩欧美| 久久久久国产精品人妻一区二区| 人人妻人人添人人爽欧美一区卜 | 日韩成人av中文字幕在线观看| 日韩av免费高清视频| 国产精品精品国产色婷婷| 国产精品人妻久久久影院| 久久人人爽人人爽人人片va| 亚洲国产欧美人成| 一本—道久久a久久精品蜜桃钙片| 一区二区三区精品91| 天天躁夜夜躁狠狠久久av| 91aial.com中文字幕在线观看| 精品人妻偷拍中文字幕| 久久久亚洲精品成人影院| av在线观看视频网站免费| 久久人妻熟女aⅴ| 男女啪啪激烈高潮av片| 中文天堂在线官网| 日日啪夜夜撸| 久久久午夜欧美精品| 午夜日本视频在线| 欧美精品一区二区大全| 国产精品.久久久| 99热这里只有是精品50| 亚洲人成网站高清观看| 国产免费又黄又爽又色| 成人一区二区视频在线观看| 亚洲国产精品国产精品| 一级爰片在线观看| 国内精品宾馆在线| 欧美日韩国产mv在线观看视频 | 国产v大片淫在线免费观看| 久久人人爽人人片av| 韩国av在线不卡| 久久精品久久久久久久性| 久久精品久久久久久久性| 啦啦啦啦在线视频资源| 天堂俺去俺来也www色官网| 熟妇人妻不卡中文字幕| 免费黄频网站在线观看国产| 久久久久久久久久久丰满| 丝袜喷水一区| 久久婷婷青草| 男女无遮挡免费网站观看| 久久国产精品男人的天堂亚洲 | 亚洲精品aⅴ在线观看| 蜜臀久久99精品久久宅男| 全区人妻精品视频| 午夜日本视频在线| 日本色播在线视频| 美女国产视频在线观看| 激情五月婷婷亚洲| 亚洲av电影在线观看一区二区三区| 国产精品一区二区在线不卡| 亚洲欧美日韩卡通动漫| 卡戴珊不雅视频在线播放| 插阴视频在线观看视频| 2021少妇久久久久久久久久久| 国产一级毛片在线| 午夜精品国产一区二区电影| 亚洲第一av免费看| 人妻一区二区av| 少妇高潮的动态图| 免费观看av网站的网址| freevideosex欧美| 五月伊人婷婷丁香| 日日撸夜夜添| 一级毛片 在线播放| 各种免费的搞黄视频| 一级毛片电影观看| 国精品久久久久久国模美| 啦啦啦视频在线资源免费观看| 王馨瑶露胸无遮挡在线观看| 久久青草综合色| 插阴视频在线观看视频| 99热这里只有是精品在线观看| 色婷婷久久久亚洲欧美| 极品教师在线视频| 亚洲国产毛片av蜜桃av| 国产色婷婷99| 尾随美女入室| 最近中文字幕2019免费版| 在线看a的网站| 婷婷色综合大香蕉| 97超视频在线观看视频| 天堂8中文在线网| 亚洲精品久久午夜乱码| 人人妻人人添人人爽欧美一区卜 | 日本vs欧美在线观看视频 | 黑丝袜美女国产一区| 一级毛片黄色毛片免费观看视频| 人妻夜夜爽99麻豆av| 纵有疾风起免费观看全集完整版| 免费人妻精品一区二区三区视频| 国产成人91sexporn| 亚洲综合色惰| 欧美少妇被猛烈插入视频| 91精品伊人久久大香线蕉| 我要看黄色一级片免费的| 少妇熟女欧美另类| 特大巨黑吊av在线直播| 久久精品熟女亚洲av麻豆精品| 亚洲久久久国产精品| 高清黄色对白视频在线免费看 | 国产一级毛片在线| h视频一区二区三区| 国产在线免费精品| 成人毛片a级毛片在线播放| 日韩大片免费观看网站| 欧美性感艳星| 一本一本综合久久| 狠狠精品人妻久久久久久综合| 亚洲人成网站在线观看播放| 小蜜桃在线观看免费完整版高清| 大香蕉久久网| 国产乱人视频| 国产成人91sexporn| 激情 狠狠 欧美| 成人无遮挡网站| 国产淫片久久久久久久久| 丰满少妇做爰视频| 国产精品久久久久久精品古装| 一个人免费看片子| 精品久久久噜噜| 国产成人一区二区在线| 亚洲欧美精品自产自拍| 国产精品偷伦视频观看了| 国产一区二区三区av在线| av免费在线看不卡| 一二三四中文在线观看免费高清| 日本欧美国产在线视频| 亚洲人与动物交配视频| 国产亚洲一区二区精品| 日韩视频在线欧美| 国产白丝娇喘喷水9色精品| 国产av码专区亚洲av| 哪个播放器可以免费观看大片| 777米奇影视久久| 久久人人爽av亚洲精品天堂 | 极品少妇高潮喷水抽搐| 老女人水多毛片| 欧美人与善性xxx| 久久久久久九九精品二区国产| 精品久久久久久电影网| 色吧在线观看| 久久精品久久精品一区二区三区| 色婷婷久久久亚洲欧美| 国产大屁股一区二区在线视频| 国产欧美日韩精品一区二区| 成人免费观看视频高清| 国产永久视频网站| 美女国产视频在线观看| 亚洲欧美日韩卡通动漫| 性色av一级| 日日啪夜夜爽| 内射极品少妇av片p| 男女边摸边吃奶| 少妇被粗大猛烈的视频| 日本av手机在线免费观看| 六月丁香七月| 亚洲成人手机| 91在线精品国自产拍蜜月| 亚洲欧美清纯卡通| 国产91av在线免费观看| 不卡视频在线观看欧美| 熟女电影av网| 91久久精品电影网| 国产精品人妻久久久影院| 午夜老司机福利剧场| 免费看日本二区| 蜜桃久久精品国产亚洲av| 中文资源天堂在线| 国产av精品麻豆| 国产真实伦视频高清在线观看| 99九九线精品视频在线观看视频| 男人和女人高潮做爰伦理| 国产乱来视频区| 国产欧美日韩精品一区二区| 亚洲成人av在线免费| 亚洲欧美精品专区久久| 日韩强制内射视频| 一区二区三区乱码不卡18| 国产久久久一区二区三区| 亚洲久久久国产精品| 人人妻人人澡人人爽人人夜夜| 自拍偷自拍亚洲精品老妇| 日韩中文字幕视频在线看片 | 亚洲av不卡在线观看| 久久久久久久久久久丰满| 亚洲国产日韩一区二区| 亚洲激情五月婷婷啪啪| 国产成人91sexporn| 啦啦啦视频在线资源免费观看| 狂野欧美激情性bbbbbb| 夜夜爽夜夜爽视频| 中文资源天堂在线| 国产午夜精品久久久久久一区二区三区| 色哟哟·www| 色婷婷久久久亚洲欧美| 三级国产精品片| 亚洲精华国产精华液的使用体验| 亚洲精品亚洲一区二区| 五月玫瑰六月丁香| 高清不卡的av网站| 欧美3d第一页| 国产欧美日韩精品一区二区| 国产免费福利视频在线观看| 777米奇影视久久| 美女xxoo啪啪120秒动态图| 亚洲国产精品成人久久小说| av在线app专区| 人人妻人人澡人人爽人人夜夜| 男女下面进入的视频免费午夜| 又大又黄又爽视频免费| 久久综合国产亚洲精品| 国产成人aa在线观看| 精品视频人人做人人爽| 你懂的网址亚洲精品在线观看| 一区二区av电影网| 亚洲精品国产色婷婷电影| 国产成人aa在线观看| 国产av一区二区精品久久 | 九草在线视频观看| 日产精品乱码卡一卡2卡三| 国产精品爽爽va在线观看网站| 国产精品嫩草影院av在线观看| av在线app专区| 一个人看的www免费观看视频| 美女主播在线视频| 色视频在线一区二区三区| 如何舔出高潮| 大香蕉久久网| 亚洲欧洲国产日韩| 91精品国产国语对白视频| 亚洲精品日本国产第一区| av在线app专区| 国产成人一区二区在线| 亚洲av国产av综合av卡| 日本猛色少妇xxxxx猛交久久| 婷婷色综合www| 新久久久久国产一级毛片| 小蜜桃在线观看免费完整版高清| av在线蜜桃| 亚洲自偷自拍三级| 午夜福利影视在线免费观看| 久久精品国产亚洲av天美| 国产精品99久久久久久久久| 亚洲成色77777| 哪个播放器可以免费观看大片| 国产视频首页在线观看| 少妇人妻精品综合一区二区| 夜夜爽夜夜爽视频| 色5月婷婷丁香| 人人妻人人澡人人爽人人夜夜| 国产一区二区在线观看日韩| 啦啦啦中文免费视频观看日本| 国产高清三级在线| 日日摸夜夜添夜夜添av毛片| 亚洲av男天堂| 高清欧美精品videossex| 国产成人freesex在线| 免费久久久久久久精品成人欧美视频 | 亚洲国产精品国产精品| 色综合色国产| 亚洲天堂av无毛| 亚洲美女搞黄在线观看| 亚洲久久久国产精品| 欧美亚洲 丝袜 人妻 在线| 最新中文字幕久久久久| 国产综合精华液| 一区二区三区乱码不卡18| 国产v大片淫在线免费观看| 成人综合一区亚洲| 久久久久久久久久久免费av| 亚洲av成人精品一二三区| 国国产精品蜜臀av免费| 亚洲婷婷狠狠爱综合网| 久久久久久久久久久丰满| 成人亚洲欧美一区二区av| 亚洲国产欧美在线一区| 亚洲欧美日韩卡通动漫| 中文字幕人妻熟人妻熟丝袜美| 欧美日韩综合久久久久久| 中国国产av一级| 国产在线一区二区三区精| 国产爱豆传媒在线观看| 日本爱情动作片www.在线观看| 久久精品夜色国产| 久久精品国产亚洲av涩爱| 午夜福利高清视频| 国产精品精品国产色婷婷| 亚洲精品第二区| 国产精品99久久久久久久久| 亚洲精品久久久久久婷婷小说| 国产成人免费无遮挡视频| 国产欧美日韩一区二区三区在线 | tube8黄色片| 欧美精品一区二区大全| 中国国产av一级| 99热6这里只有精品| 国产探花极品一区二区| 七月丁香在线播放| 中文乱码字字幕精品一区二区三区| 日韩在线高清观看一区二区三区| 亚洲成人中文字幕在线播放| 91久久精品国产一区二区成人| 国产在线一区二区三区精| 波野结衣二区三区在线| 高清日韩中文字幕在线| 中文乱码字字幕精品一区二区三区| 国产毛片在线视频| 国产精品一及| 精品一区二区三区视频在线| 丝瓜视频免费看黄片| 亚洲国产日韩一区二区| 日韩制服骚丝袜av| 插阴视频在线观看视频| 国产av一区二区精品久久 | 最近中文字幕2019免费版| 黑人高潮一二区| 亚洲精品第二区| 国产在线免费精品| 最新中文字幕久久久久| 夜夜骑夜夜射夜夜干| 欧美人与善性xxx| 亚洲国产精品专区欧美| 人妻少妇偷人精品九色| 麻豆成人av视频| 亚洲av成人精品一二三区| 伦理电影免费视频| 熟女人妻精品中文字幕| 欧美老熟妇乱子伦牲交| 亚洲激情五月婷婷啪啪| 91久久精品国产一区二区三区| 久久人人爽人人片av| 日本黄色日本黄色录像| av国产免费在线观看| 自拍欧美九色日韩亚洲蝌蚪91 | 国产深夜福利视频在线观看| 亚洲欧美日韩东京热| 久久久久久久久大av| 天天躁日日操中文字幕| 亚洲av在线观看美女高潮| 国产高潮美女av| 男女国产视频网站| 国产深夜福利视频在线观看| 欧美精品一区二区大全| av国产精品久久久久影院| 最近最新中文字幕大全电影3| 久久久久久久久久人人人人人人| 免费不卡的大黄色大毛片视频在线观看| 久久久久国产精品人妻一区二区| 人人妻人人澡人人爽人人夜夜| 午夜激情久久久久久久| 在线亚洲精品国产二区图片欧美 | 黑人猛操日本美女一级片| 日韩成人伦理影院| 少妇丰满av| 久久精品国产亚洲网站| 亚洲精品自拍成人| 欧美xxxx性猛交bbbb| 纵有疾风起免费观看全集完整版| 国产精品久久久久成人av| 国产av精品麻豆| 国产免费一区二区三区四区乱码| 五月开心婷婷网| 亚洲色图综合在线观看| 欧美+日韩+精品| 视频中文字幕在线观看| 亚洲无线观看免费| 成人影院久久| 国产又色又爽无遮挡免| 激情五月婷婷亚洲| 国产av一区二区精品久久 | 校园人妻丝袜中文字幕| 99re6热这里在线精品视频| 噜噜噜噜噜久久久久久91| 亚洲人成网站在线观看播放| 欧美少妇被猛烈插入视频| 日韩电影二区| xxx大片免费视频| 欧美一级a爱片免费观看看| 大片电影免费在线观看免费| 成人一区二区视频在线观看| 日本猛色少妇xxxxx猛交久久| 国产成人a区在线观看| 亚洲精品色激情综合| 成人综合一区亚洲| 精品久久国产蜜桃| 人妻 亚洲 视频| 成人二区视频| 亚洲精品视频女| 在线免费观看不下载黄p国产| 亚洲一区二区三区欧美精品| 一级黄片播放器| 蜜臀久久99精品久久宅男| 在线免费观看不下载黄p国产| 婷婷色麻豆天堂久久| 韩国av在线不卡| 日本与韩国留学比较| 国产精品久久久久久久久免| 中文欧美无线码| 男人狂女人下面高潮的视频| 国产一区二区在线观看日韩| 纯流量卡能插随身wifi吗| 最新中文字幕久久久久| 国产爽快片一区二区三区| 久久99热这里只有精品18| 亚洲国产精品国产精品| 日本vs欧美在线观看视频 | 国产欧美日韩一区二区三区在线 | 蜜桃久久精品国产亚洲av| 男女下面进入的视频免费午夜| 亚洲精品乱码久久久v下载方式| 日日撸夜夜添| 亚洲国产精品国产精品| 黑人高潮一二区| 男女无遮挡免费网站观看| 狂野欧美激情性xxxx在线观看| 国产精品人妻久久久久久| 午夜免费鲁丝| 亚洲欧美日韩东京热| 国产成人免费无遮挡视频| av福利片在线观看| 一本一本综合久久| 国产乱人视频| 天堂俺去俺来也www色官网| 免费观看a级毛片全部|