• <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模型整體解的存在性
    亚洲aⅴ乱码一区二区在线播放| 两个人的视频大全免费| 嘟嘟电影网在线观看| 精品久久久久久久末码| 亚洲精品日韩在线中文字幕| 中国美白少妇内射xxxbb| 99久久成人亚洲精品观看| 亚洲一区高清亚洲精品| 桃色一区二区三区在线观看| 狂野欧美白嫩少妇大欣赏| 91久久精品国产一区二区成人| 97超视频在线观看视频| 日日啪夜夜撸| 亚洲精品日韩在线中文字幕| 日本与韩国留学比较| 久久精品国产自在天天线| 国产高潮美女av| 熟妇人妻久久中文字幕3abv| 亚洲国产精品成人久久小说| 天堂√8在线中文| 婷婷色综合大香蕉| 草草在线视频免费看| 欧美日韩精品成人综合77777| 成人av在线播放网站| 亚洲av.av天堂| 欧美一区二区亚洲| 禁无遮挡网站| 国产精品久久电影中文字幕| 成人午夜精彩视频在线观看| 日本黄色视频三级网站网址| 国产私拍福利视频在线观看| 能在线免费看毛片的网站| 青春草国产在线视频| 国产成人a区在线观看| 免费无遮挡裸体视频| 久久99精品国语久久久| 又粗又硬又长又爽又黄的视频| 国产黄片视频在线免费观看| 免费看av在线观看网站| 精品少妇黑人巨大在线播放 | 亚洲四区av| 久久久久久久国产电影| 国产黄片美女视频| 国产私拍福利视频在线观看| 精品欧美国产一区二区三| 日韩欧美国产在线观看| 天堂中文最新版在线下载 | 老司机影院毛片| 在线观看一区二区三区| 日日干狠狠操夜夜爽| 久久久国产成人精品二区| 久久亚洲国产成人精品v| 在线观看av片永久免费下载| 麻豆国产97在线/欧美| 国产亚洲精品久久久com| 欧美极品一区二区三区四区| 一个人看的www免费观看视频| 五月伊人婷婷丁香| av卡一久久| 成人二区视频| 亚洲经典国产精华液单| 亚州av有码| 97热精品久久久久久| 别揉我奶头 嗯啊视频| 精品酒店卫生间| 亚洲成人中文字幕在线播放| 国产成人a∨麻豆精品| 欧美成人免费av一区二区三区| 久久久久久久久久成人| 亚洲欧美成人综合另类久久久 | 亚洲av日韩在线播放| 亚洲av免费在线观看| 国产在视频线精品| 中文字幕久久专区| 久久精品人妻少妇| 日本av手机在线免费观看| 一级黄片播放器| a级一级毛片免费在线观看| 午夜福利高清视频| 久久综合国产亚洲精品| 中文字幕亚洲精品专区| 免费无遮挡裸体视频| 亚洲美女视频黄频| 一区二区三区乱码不卡18| 高清在线视频一区二区三区 | 亚洲欧美日韩卡通动漫| 级片在线观看| 永久免费av网站大全| 亚洲乱码一区二区免费版| 免费观看a级毛片全部| 亚洲av福利一区| 99在线人妻在线中文字幕| 欧美一区二区国产精品久久精品| 国产精品人妻久久久影院| 丰满人妻一区二区三区视频av| 精品久久久久久久久亚洲| 久久热精品热| 国产黄片美女视频| 亚洲国产高清在线一区二区三| 精品少妇黑人巨大在线播放 | 91在线精品国自产拍蜜月| 青青草视频在线视频观看| 国产精品三级大全| 最近中文字幕高清免费大全6| 国产高清有码在线观看视频| 最近手机中文字幕大全| 联通29元200g的流量卡| av在线观看视频网站免费| 一级毛片久久久久久久久女| 亚洲美女搞黄在线观看| av在线蜜桃| 国产av不卡久久| 国产淫片久久久久久久久| 可以在线观看毛片的网站| 91午夜精品亚洲一区二区三区| 丰满乱子伦码专区| 欧美成人精品欧美一级黄| 黄色配什么色好看| 亚洲一级一片aⅴ在线观看| 亚洲精品国产av成人精品| 亚洲人成网站高清观看| 91午夜精品亚洲一区二区三区| 一级二级三级毛片免费看| 久久久国产成人精品二区| 国内精品一区二区在线观看| 免费电影在线观看免费观看| 国产美女午夜福利| 少妇熟女aⅴ在线视频| 看非洲黑人一级黄片| 久久韩国三级中文字幕| 夜夜爽夜夜爽视频| 久久久午夜欧美精品| 亚洲av一区综合| 韩国高清视频一区二区三区| 三级经典国产精品| 99在线人妻在线中文字幕| 少妇人妻精品综合一区二区| 国内精品一区二区在线观看| 国产亚洲91精品色在线| 免费看av在线观看网站| 久久精品91蜜桃| 在线观看美女被高潮喷水网站| 国产精品久久久久久久久免| 欧美三级亚洲精品| av在线天堂中文字幕| 婷婷色av中文字幕| 亚洲中文字幕日韩| 免费大片18禁| 国产综合懂色| 久久99精品国语久久久| 亚洲精品国产av成人精品| 国产欧美日韩精品一区二区| 婷婷色综合大香蕉| 精品欧美国产一区二区三| 精品久久久久久久人妻蜜臀av| 波多野结衣高清无吗| 观看美女的网站| 日本三级黄在线观看| 亚洲电影在线观看av| 国产一区二区亚洲精品在线观看| 午夜精品国产一区二区电影 | 亚洲在线自拍视频| 久久99精品国语久久久| 亚洲综合精品二区| 尾随美女入室| 国产精品一及| 亚洲一区高清亚洲精品| 欧美日韩一区二区视频在线观看视频在线 | 一夜夜www| 亚洲久久久久久中文字幕| 熟女电影av网| 免费看美女性在线毛片视频| 国产色爽女视频免费观看| 日韩av在线免费看完整版不卡| 少妇裸体淫交视频免费看高清| 国产精品久久电影中文字幕| 亚洲精品日韩在线中文字幕| 欧美xxxx黑人xx丫x性爽| 亚洲av成人精品一区久久| 亚洲精品亚洲一区二区| 国模一区二区三区四区视频| 在线播放无遮挡| 男女下面进入的视频免费午夜| 男插女下体视频免费在线播放| 美女被艹到高潮喷水动态| 亚洲自拍偷在线| 日韩高清综合在线| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲国产欧美在线一区| 18禁动态无遮挡网站| av.在线天堂| 综合色av麻豆| 老女人水多毛片| 又粗又爽又猛毛片免费看| 一夜夜www| 久久久久久久亚洲中文字幕| 国产 一区精品| 亚洲国产最新在线播放| 亚洲在久久综合| 日日撸夜夜添| 久久久国产成人精品二区| 国产精品久久久久久久电影| 在线观看美女被高潮喷水网站| 亚洲天堂国产精品一区在线| 美女大奶头视频| 国产精品久久久久久av不卡| 亚洲精华国产精华液的使用体验| 美女国产视频在线观看| 免费看av在线观看网站| 22中文网久久字幕| 美女cb高潮喷水在线观看| 在线播放无遮挡| 国产精品久久久久久精品电影| 国产精品一区www在线观看| 日韩,欧美,国产一区二区三区 | 日日摸夜夜添夜夜添av毛片| 亚州av有码| 婷婷色av中文字幕| 亚洲精品乱码久久久久久按摩| 国产成人a∨麻豆精品| 免费一级毛片在线播放高清视频| 美女黄网站色视频| 国产精品伦人一区二区| 久久久久久九九精品二区国产| av在线播放精品| 国产成人a区在线观看| 亚洲欧美日韩无卡精品| 能在线免费观看的黄片| 你懂的网址亚洲精品在线观看 | 在线a可以看的网站| 久久久久久国产a免费观看| 亚洲不卡免费看| 午夜福利高清视频| 久久国内精品自在自线图片| 亚洲av二区三区四区| 男女边吃奶边做爰视频| 男女那种视频在线观看| 午夜a级毛片| 中文乱码字字幕精品一区二区三区 | 国产精品一区二区三区四区免费观看| 乱人视频在线观看| 青春草视频在线免费观看| 日本色播在线视频| 99久国产av精品国产电影| 国产成人一区二区在线| 亚洲av免费在线观看| 日本一二三区视频观看| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲av免费在线观看| 在现免费观看毛片| 欧美最新免费一区二区三区| 日本av手机在线免费观看| 国产探花在线观看一区二区| 久久久久久久久大av| 99热这里只有是精品在线观看| 成年女人看的毛片在线观看| 性插视频无遮挡在线免费观看| 亚洲国产日韩欧美精品在线观看| 国产精品,欧美在线| 久久久国产成人精品二区| 91久久精品电影网| 嫩草影院入口| 日韩,欧美,国产一区二区三区 | 少妇高潮的动态图| 国产亚洲最大av| 国产在线一区二区三区精 | 最后的刺客免费高清国语| 精品少妇黑人巨大在线播放 | 少妇高潮的动态图| 免费观看性生交大片5| 国产欧美日韩精品一区二区| 国产亚洲91精品色在线| 高清在线视频一区二区三区 | 秋霞在线观看毛片| 边亲边吃奶的免费视频| 在现免费观看毛片| 永久网站在线| 亚洲va在线va天堂va国产| 国产亚洲精品av在线| 久久精品影院6| 久久精品国产亚洲av天美| 国产精品福利在线免费观看| 欧美xxxx性猛交bbbb| 中文亚洲av片在线观看爽| 天美传媒精品一区二区| 亚洲国产日韩欧美精品在线观看| 国产精品爽爽va在线观看网站| 久久久a久久爽久久v久久| 少妇丰满av| 久久久久性生活片| 日韩人妻高清精品专区| 亚洲国产色片| 久久精品91蜜桃| 国产精品久久久久久av不卡| 日本免费a在线| 久久这里有精品视频免费| 亚洲av.av天堂| 99久国产av精品| 久久久午夜欧美精品| 国产片特级美女逼逼视频| 精品国产一区二区三区久久久樱花 | 乱系列少妇在线播放| 人体艺术视频欧美日本| 一个人观看的视频www高清免费观看| 亚洲国产日韩欧美精品在线观看| 中文字幕熟女人妻在线| 日韩视频在线欧美| 久久精品影院6| 真实男女啪啪啪动态图| 99热网站在线观看| 亚洲国产精品久久男人天堂| 久久精品国产自在天天线| 超碰97精品在线观看| 久久精品人妻少妇| 免费黄色在线免费观看| 婷婷色麻豆天堂久久 | 欧美日韩精品成人综合77777| 免费在线观看成人毛片| 日韩精品青青久久久久久| 亚洲欧美中文字幕日韩二区| 国产极品精品免费视频能看的| 亚洲精品影视一区二区三区av| 欧美激情国产日韩精品一区| 嫩草影院新地址| 91狼人影院| a级毛片免费高清观看在线播放| 亚洲人成网站高清观看| 精品午夜福利在线看| 最近最新中文字幕大全电影3| 国内揄拍国产精品人妻在线| 国产91av在线免费观看| 搡女人真爽免费视频火全软件| 嘟嘟电影网在线观看| 2022亚洲国产成人精品| 99热全是精品| 蜜臀久久99精品久久宅男| 国产麻豆成人av免费视频| 午夜福利在线观看免费完整高清在| 亚洲欧美日韩卡通动漫| 中文字幕久久专区| 久久精品久久久久久噜噜老黄 | 日韩 亚洲 欧美在线| 欧美成人午夜免费资源| 久久精品人妻少妇| 一区二区三区免费毛片| 色5月婷婷丁香| 国产黄a三级三级三级人| 直男gayav资源| 日本五十路高清| 一区二区三区四区激情视频| 夜夜爽夜夜爽视频| 黄片wwwwww| 亚洲人与动物交配视频| 老司机影院毛片| 国国产精品蜜臀av免费| 男女那种视频在线观看| 精品国产三级普通话版| 美女高潮的动态| 成人毛片a级毛片在线播放| 国产精品久久久久久av不卡| av在线亚洲专区| 男女边吃奶边做爰视频| 别揉我奶头 嗯啊视频| 国产一区二区三区av在线| 国产黄色视频一区二区在线观看 | 国产精品人妻久久久影院| 蜜桃久久精品国产亚洲av| 国产三级在线视频| 亚洲国产精品成人久久小说| 免费无遮挡裸体视频| 七月丁香在线播放| 麻豆精品久久久久久蜜桃| 三级国产精品欧美在线观看| 成人无遮挡网站| 午夜免费男女啪啪视频观看| 亚洲精品日韩在线中文字幕| 插阴视频在线观看视频| 免费观看精品视频网站| 深夜a级毛片| 欧美色视频一区免费| 亚洲最大成人中文| 久久精品夜夜夜夜夜久久蜜豆| 精品一区二区三区视频在线| 18+在线观看网站| 久久这里只有精品中国| 一本一本综合久久| 中文资源天堂在线| 日韩欧美三级三区| 亚洲成人久久爱视频| 国产成人精品一,二区| 中文天堂在线官网| 两个人视频免费观看高清| 男女那种视频在线观看| 日本免费一区二区三区高清不卡| 欧美成人a在线观看| 看片在线看免费视频| 亚洲成人中文字幕在线播放| 精品人妻熟女av久视频| 国产老妇女一区| 日韩欧美精品免费久久| 亚洲电影在线观看av| 国产老妇伦熟女老妇高清| 国产女主播在线喷水免费视频网站 | 亚洲国产精品成人久久小说| 久久久午夜欧美精品| 美女高潮的动态| 精华霜和精华液先用哪个| 精品久久久久久成人av| 99热网站在线观看| 国产淫语在线视频| or卡值多少钱| 国产探花极品一区二区| 免费搜索国产男女视频| 国产av码专区亚洲av| 免费观看的影片在线观看| 寂寞人妻少妇视频99o| 精品一区二区三区视频在线| 日本五十路高清| 嘟嘟电影网在线观看| 免费播放大片免费观看视频在线观看 | 天天一区二区日本电影三级| 91精品伊人久久大香线蕉| 插逼视频在线观看| 亚洲国产色片| av播播在线观看一区| 18禁动态无遮挡网站| 久久99热这里只有精品18| 中文在线观看免费www的网站| 亚洲18禁久久av| 日韩一区二区视频免费看| 欧美性猛交黑人性爽| 欧美人与善性xxx| 色播亚洲综合网| 夫妻性生交免费视频一级片| 日本欧美国产在线视频| 日韩一本色道免费dvd| 久久精品国产亚洲网站| 久久久久精品久久久久真实原创| 国产精品麻豆人妻色哟哟久久 | 国产黄色视频一区二区在线观看 | 中文在线观看免费www的网站| 亚洲av福利一区| 久久久欧美国产精品| 嫩草影院新地址| 国产精品1区2区在线观看.| 毛片女人毛片| 亚洲精品乱码久久久久久按摩| 亚洲一区高清亚洲精品| av国产免费在线观看| 久久99热6这里只有精品| 男人舔女人下体高潮全视频| 国产成人精品一,二区| 成人特级av手机在线观看| 搞女人的毛片| 欧美xxxx黑人xx丫x性爽| 国产免费一级a男人的天堂| 亚洲人与动物交配视频| 免费观看在线日韩| av国产久精品久网站免费入址| 女人久久www免费人成看片 | 九草在线视频观看| 欧美激情在线99| 国产免费男女视频| 免费观看人在逋| av在线亚洲专区| 国产精品精品国产色婷婷| 麻豆一二三区av精品| 国国产精品蜜臀av免费| 成人午夜精彩视频在线观看| 成年女人永久免费观看视频| 99热精品在线国产| 男女边吃奶边做爰视频| 亚洲丝袜综合中文字幕| 国产成人福利小说| av专区在线播放| 欧美一区二区国产精品久久精品| 久久久精品大字幕| 国产成人91sexporn| 免费黄网站久久成人精品| 日韩欧美精品v在线| 能在线免费观看的黄片| 汤姆久久久久久久影院中文字幕 | 精品久久久噜噜| 我要看日韩黄色一级片| 老司机影院毛片| 成人av在线播放网站| av播播在线观看一区| 在线播放国产精品三级| 国产一区二区在线av高清观看| 国产午夜精品论理片| 久久午夜福利片| 22中文网久久字幕| 极品教师在线视频| 国产成年人精品一区二区| 国产成人午夜福利电影在线观看| av卡一久久| 一区二区三区高清视频在线| 男人狂女人下面高潮的视频| 久久久久久久久久黄片| 国产精品不卡视频一区二区| 久久午夜福利片| 国产乱人视频| 亚洲美女搞黄在线观看| 国产一区二区亚洲精品在线观看| 99久久九九国产精品国产免费| 好男人在线观看高清免费视频| 国产午夜精品久久久久久一区二区三区| 免费黄网站久久成人精品| 99在线人妻在线中文字幕| 久久精品国产亚洲av涩爱| 春色校园在线视频观看| 蜜臀久久99精品久久宅男| 99久久人妻综合| 97在线视频观看| 成年版毛片免费区| 国产美女午夜福利| 亚洲欧洲国产日韩| 成人欧美大片| 一个人观看的视频www高清免费观看| 永久免费av网站大全| 国产v大片淫在线免费观看| 国产一区亚洲一区在线观看| 内射极品少妇av片p| 久久国内精品自在自线图片| 亚洲欧美日韩无卡精品| 欧美一级a爱片免费观看看| 一个人免费在线观看电影| 国产av码专区亚洲av| 国产精品国产三级国产专区5o | 久久久欧美国产精品| 亚洲五月天丁香| 午夜福利高清视频| 亚洲国产精品久久男人天堂| 亚洲精品国产成人久久av| av天堂中文字幕网| 人妻系列 视频| 国产私拍福利视频在线观看| 亚洲性久久影院| 99久久精品国产国产毛片| 国产精品无大码| 成人美女网站在线观看视频| 91在线精品国自产拍蜜月| 成人毛片60女人毛片免费| 亚洲国产成人一精品久久久| 国产亚洲午夜精品一区二区久久 | 不卡视频在线观看欧美| 精品欧美国产一区二区三| 亚洲精品成人久久久久久| 国产淫语在线视频| 国产精品久久久久久av不卡| 婷婷色麻豆天堂久久 | 久久久精品欧美日韩精品| 99视频精品全部免费 在线| 麻豆av噜噜一区二区三区| 我的女老师完整版在线观看| 国产精品99久久久久久久久| 国产老妇伦熟女老妇高清| 26uuu在线亚洲综合色| 亚洲丝袜综合中文字幕| 六月丁香七月| 人人妻人人澡欧美一区二区| 国产 一区 欧美 日韩| 欧美性猛交╳xxx乱大交人| 亚洲电影在线观看av| 99久国产av精品| 日本猛色少妇xxxxx猛交久久| 三级经典国产精品| 国产精品国产三级国产专区5o | 91狼人影院| 亚洲av免费高清在线观看| 亚洲精品日韩av片在线观看| 亚洲三级黄色毛片| 毛片一级片免费看久久久久| 汤姆久久久久久久影院中文字幕 | 色5月婷婷丁香| 午夜精品在线福利| 午夜a级毛片| 日韩欧美在线乱码| 国语对白做爰xxxⅹ性视频网站| 国产一区二区三区av在线| 99热精品在线国产| 一区二区三区免费毛片| 国产一区二区三区av在线| 联通29元200g的流量卡| 18禁在线无遮挡免费观看视频| 久久久亚洲精品成人影院| 乱码一卡2卡4卡精品| 久久99热这里只频精品6学生 | 国内精品宾馆在线| www日本黄色视频网| 国产欧美日韩精品一区二区| 成人无遮挡网站| 久久国产乱子免费精品| 免费人成在线观看视频色| 好男人视频免费观看在线| 黄片无遮挡物在线观看| 一级毛片aaaaaa免费看小| 97热精品久久久久久| 少妇人妻一区二区三区视频| 成人一区二区视频在线观看| 亚洲国产最新在线播放| 少妇人妻一区二区三区视频| 免费一级毛片在线播放高清视频| 亚洲va在线va天堂va国产| av国产免费在线观看| 午夜免费男女啪啪视频观看| 亚洲最大成人手机在线| 久久精品国产自在天天线| 亚洲欧美成人精品一区二区| 18禁裸乳无遮挡免费网站照片| 亚洲精品成人久久久久久| 国产精品人妻久久久久久| 老司机影院成人| 久久这里只有精品中国| 高清视频免费观看一区二区 |