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

    全應(yīng)變梯度撓曲電納米梁有限單元法研究

    2020-10-27 03:16:18陳玲玲王炳雷
    關(guān)鍵詞:效應(yīng)模型

    陳玲玲, 楊 旭, 劉 洋, 王炳雷

    (山東大學(xué) 土建與水利學(xué)院工程力學(xué)系,濟(jì)南 250061)

    1 引 言

    壓電效應(yīng)是最常見的力電耦合效應(yīng),是應(yīng)變?chǔ)藕碗姌O化P之間的線性耦合[1]。雖然在生活中有著廣泛的應(yīng)用,但壓電效應(yīng)不是普遍存在的。因?yàn)樗淮嬖谟诜侵行膶?duì)稱介電材料中,并且其激發(fā)還受限于居里溫度Tc。

    不同于壓電效應(yīng),撓曲電效應(yīng)是一種普遍存在的力電耦合效應(yīng),是應(yīng)變梯度?ε/?x與電極化P之間的線性耦合。理論上,撓曲電效應(yīng)存在于所有的電介質(zhì)材料中,如液晶、聚合物和脂質(zhì)雙層膜等[2-4]。它在能量收集、傳感和驅(qū)動(dòng)、高級(jí)顯微鏡、人工肌肉和微創(chuàng)手術(shù)方面有廣泛的應(yīng)用前景[5,6]。撓曲電效應(yīng)的激發(fā)沒有溫度限制,并且其在納米尺度下會(huì)變得更加顯著(即存在尺寸依賴性),這是因?yàn)殡S著結(jié)構(gòu)尺寸的減小,應(yīng)變梯度會(huì)相應(yīng)地增大。

    近些年,撓曲電效應(yīng)引起了人們的廣泛關(guān)注,主要原因有兩個(gè)。(1) 鐵電陶瓷的撓曲電效應(yīng)試驗(yàn)表明,一些高介電常數(shù)的材料可以獲得比預(yù)期更強(qiáng)的撓曲電特性[7];(2) 隨著先進(jìn)的納米技術(shù)和小型化趨勢的發(fā)展,納米結(jié)構(gòu)材料得到了廣泛的開發(fā)。撓曲電效應(yīng)最重要的應(yīng)用之一就是利用其尺寸依賴性而制成的俘能器[8]。同時(shí),撓曲電效應(yīng)還可以顯著提高納米梁的壓電性,這為其在微納機(jī)電系統(tǒng)中的應(yīng)用提供了新的途徑。

    Kogan[9]為了在理論上解釋電介質(zhì)中的撓曲電效應(yīng),首次提出了撓曲電理論,并估計(jì)了撓曲電系數(shù)的取值范圍。申勝平等[10]從電焓學(xué)的角度出發(fā),同時(shí)考慮了撓曲電效應(yīng)、表面效應(yīng)和靜電效應(yīng),嚴(yán)格推導(dǎo)出了納米電介質(zhì)的控制方程。周亞榮等[11]給出了考慮撓曲電效應(yīng)的納米歐拉梁的理論框架。王炳雷等[12]研究了考慮撓曲電效應(yīng)的揉皺電介質(zhì)薄膜的能量俘獲機(jī)制。

    應(yīng)變梯度是引起撓曲電效應(yīng)的直接原因,為了準(zhǔn)確地描述微構(gòu)件中的撓曲電效應(yīng),必然要引入應(yīng)變梯度理論。應(yīng)變梯度理論的發(fā)展已有半個(gè)多世紀(jì),應(yīng)變梯度理論最先由Mindlin[13]提出,除了 2個(gè) Lamé常數(shù)之外,對(duì)于各向同性材料還需要額外的16個(gè)獨(dú)立的材料尺度參數(shù)。Mindlin等[14]對(duì)其之前的理論進(jìn)行了簡化,使材料尺度參數(shù)從 16個(gè) 減少到了5個(gè),但仍很難通過實(shí)驗(yàn)確定。周慎杰等[15]提出了全應(yīng)變梯度彈性理論,通過對(duì)兩組應(yīng)變梯度張量進(jìn)行正交分解,嚴(yán)格證明了只有3個(gè)材料尺度參數(shù)是獨(dú)立的。這3個(gè)材料尺度參數(shù)l0,l1和l2分別與膨脹梯度張量、拉伸梯度張量和轉(zhuǎn)動(dòng)梯度張量相關(guān)。

    另一方面,在撓曲電理論模型中必須考慮電極化P和應(yīng)變梯度u之間的耦合,然而對(duì)于總能量來說,應(yīng)變梯度彈性項(xiàng)(U~u·gu/2)也是至關(guān)重要的一項(xiàng)。但有些研究并未考慮這一項(xiàng)的影響。鄧謙[16]的研究表明,以往的研究低估了應(yīng)變梯度彈性項(xiàng)的影響。納米尺度下,其影響是不可忽視的。撓曲電的尺寸效應(yīng)只表明了隨著結(jié)構(gòu)的特征尺寸減小,其撓曲電響應(yīng)會(huì)增強(qiáng),卻沒有給出在什么條件下這種增強(qiáng)現(xiàn)象會(huì)減弱,而應(yīng)變梯度彈性項(xiàng)可以抑制這種現(xiàn)象。因此,為了考慮應(yīng)變梯度彈性項(xiàng)的影響,本文采用全應(yīng)變梯度彈性理論[15]。

    然而,撓曲電問題中的高階偏微分方程常給理論求解帶來困難。為此,人們開始尋求有效的數(shù)值方法。其中,有限單元法一直是解決邊值問題的有效方法。王凱等[17]建立了考慮應(yīng)變梯度理論的二維及三維問題的自然鄰近混合伽遼金法。鄧峰等[18]提出了四種二維混合有限單元來研究撓曲電效應(yīng),隨后又創(chuàng)建了三維混合有限單元法,用來模擬撓曲電球殼和截棱錐的力電耦合行為[19]。與理論研究和實(shí)驗(yàn)分析相比,有關(guān)撓曲電效應(yīng)的數(shù)值研究相對(duì)較少,并且撓曲電效應(yīng)的有限單元法也需要高階的有限單元來滿足其連續(xù)性。

    因此,為了求解納米歐拉梁的撓曲電問題,本文提出了一種基于全應(yīng)變梯度彈性理論的包含3個(gè)獨(dú)立材料尺度參數(shù)的有限單元法,并通過一個(gè)簡單的歐拉簡支梁算例,來研究撓曲電效應(yīng)和應(yīng)變梯度彈性項(xiàng)在撓曲電問題中的影響。

    2 基本理論

    根據(jù)全應(yīng)變梯度彈性理論,在本構(gòu)關(guān)系中,除了2個(gè)經(jīng)典的拉梅常數(shù)外,還引入了3個(gè)獨(dú)立的材料尺度參數(shù)[15]。因此,考慮撓曲電效應(yīng)的線彈性材料的內(nèi)能密度函數(shù)可以表示為

    (1)

    在傳統(tǒng)的電焓密度中引入撓曲電理論,則得到新的電焓密度H的表達(dá)式為[20]

    (2)

    式中e0=8.85×10- 12F/m為真空介電常數(shù),Φ為電勢,z為梁的厚度方向。

    在整個(gè)電介質(zhì)結(jié)構(gòu)中,考慮撓曲電效應(yīng)的哈密頓變分原理為

    (3)

    式中δW為外力功的變分。

    3 納米歐拉梁的理論模型

    考慮一個(gè)長度為L、厚度為h、寬度為b的壓電納米簡支梁中的力電耦合響應(yīng),如圖1所示。以梁未變形前的中心軸為x軸,建立笛卡爾坐標(biāo)系,z軸沿著梁的厚度方向。在梁的上下表面(z=±h/2)施加一個(gè)外加電壓ΔV,在梁的中點(diǎn)處(x=L/2)施加一個(gè)集中力F。當(dāng)施加電壓ΔV時(shí),簡支梁會(huì)在z方向產(chǎn)生機(jī)械變形;當(dāng)施加集中力F時(shí),簡支梁會(huì)在z方向產(chǎn)生電極化。在下面的推導(dǎo)中,坐標(biāo)(x,y,z)由(1,2,3)表示。

    設(shè)彎曲梁的z向位移為w(x),則歐拉梁中任意一點(diǎn)的位移可以表示為

    u1(x,z)=-z(?w(x)/?x)u2(x,z)=0u3(x,z)=w(x)

    (4)

    圖1 簡支梁力電耦合響應(yīng)

    the simply supported beam

    從圖1可以看出,簡支梁的模型和載荷關(guān)于x=L/2對(duì)稱,故可取梁的一半(0

    (5)

    (6)

    式中I和A分別為梁的極慣性矩和橫截面積。

    考慮到電極化和電勢只存在于z方向,并結(jié)合式(5)的后兩個(gè)式子以及式(6)的最后一個(gè)邊界條件,可得電極化和電勢的表達(dá)式為

    (7)

    (8)

    將式(7)代入式(5)的第一個(gè)式子,可得壓電納米簡支梁的撓度控制方程以及在x=0和x=L/2處的邊界條件

    (9)

    (10)

    (11)

    由撓度控制方程及邊界條件求解得到撓度的表達(dá)式為

    w=C1+C2x+C3x2+C4x3+C5eS0x+C6e-S0x

    (12)

    式中

    (13)

    令控制方程(9)和邊界條件(10)有關(guān)電學(xué)的項(xiàng)為0,即令f1133=f1113=0,d311=0,ΔV=0,再令λ1=2,λ3=2,λ5=1,便可以得到基于Lam等[22]的應(yīng)變梯度理論的簡支梁撓度表達(dá)式

    (14)

    4 納米歐拉梁的有限元模型

    基于有限單元法理論,建立納米歐拉梁的有限元模型以研究納米梁的力電耦合特性。首先根據(jù)有限元理論,將簡支梁的控制方程(9)變換為弱形式方程

    (15)

    式中M*=-[f1133ΔV/(a33h)]A相當(dāng)于作用在梁兩端的等效彎矩。

    由式(15)可以看出,弱形式方程包含w的三階導(dǎo)數(shù),因此,本文將歐拉納米梁離散成兩節(jié)點(diǎn)六自由度,滿足C2弱連續(xù)的新歐拉梁單元,該單元包含三個(gè)獨(dú)立的材料尺度參數(shù),可以表征尺寸效應(yīng)。

    歐拉納米梁的撓度w可以由一組含有待定參數(shù)的試函數(shù)表示為

    w=Na

    (16)

    N1=1-10(x/L)3+15(x/L)4-6(x/L)5

    N2=[x/L-6(x/L)3+8(x/L)4-3(x/L)5]L

    N4=10(x/L)3-15(x/L)4+6(x/L)5

    N5=[-4(x/L)3+7(x/L)4-3(x/L)5]L

    (17)

    將撓度試函數(shù)(16)代入弱形式方程(15),就可以得到新的納米歐拉梁單元的單元?jiǎng)偠染仃嘖e,

    (18)

    式中l(wèi)e為每個(gè)單元的長度。

    集合單元?jiǎng)偠染仃嘖e形成結(jié)構(gòu)的總體剛度矩陣K,最終得到有限元的求解方程

    Kw=F

    (19)

    式中K,F(xiàn)和w分別為整體剛度矩陣、整體載荷向量和整體位移向量。

    由式(18)可知,不同于傳統(tǒng)的歐拉梁單元,新的歐拉梁單元包含3個(gè)獨(dú)立的材料尺度參數(shù),因此它可以表征尺寸效應(yīng)。并且,傳統(tǒng)的歐拉梁單元只滿足C0連續(xù),而新的歐拉梁單元可以滿足C1連續(xù)和C2弱連續(xù)。

    5 算例分析

    5.1 參數(shù)設(shè)置

    以納米歐拉簡支梁為例,分析和討論其撓曲電響應(yīng)。所選材料為BaTiO3。梁的幾何形狀設(shè)定為L=20h和b=3h。懸臂梁自由端的集中力為F=3 nN,梁上下表面的外加電壓為ΔV=-0.2 V。3個(gè)獨(dú)立的材料尺度參數(shù)設(shè)定為l0=l1=l2=10 nm。其他詳細(xì)參數(shù)設(shè)置列入表1。

    5.2 收斂性分析

    無量綱化的撓度、轉(zhuǎn)角和彎矩表示如下

    (20)

    為確保本文研究方法的收斂性,下面利用新的歐拉梁單元進(jìn)行有限元計(jì)算。在x=L/3處的無量綱化的撓度、轉(zhuǎn)角和彎矩的理論計(jì)算結(jié)果分別為3.279707,5.226900和-29.795157,表2給出了相應(yīng)的有限元計(jì)算結(jié)果??梢钥闯?本文計(jì)算方法收斂且收斂速度非常快。再對(duì)比理論計(jì)算結(jié)果和有限元計(jì)算結(jié)果可以看出,有限元計(jì)算結(jié)果最大誤差為0.09%,說明本有限單元法的計(jì)算精度足夠高。由表2可知,30個(gè)單元已經(jīng)足以滿足本文的精度要求,因此,在本文將采用30個(gè)單元進(jìn)行計(jì)算以及結(jié)果分析。

    5.3 方法驗(yàn)證

    為了驗(yàn)證本文有限單元法的精確度,圖2給出兩種模型的理論與數(shù)值結(jié)果對(duì)比。一種是基于Lam等的應(yīng)變梯度理論的簡支梁模型,式(14)給出了其理論結(jié)果。一種是基于全應(yīng)變梯度理論的簡支梁模型,即為不考慮電學(xué)相關(guān)項(xiàng)(f1133=f1113=0,d311=0,ΔV=0)的本文研究模型。從圖2可以看出,數(shù)值結(jié)果與理論結(jié)果相擬合,本文有限單元法具有很高的精度。對(duì)比兩種模型的數(shù)值發(fā)現(xiàn),以往的應(yīng)變梯度理論模型可能低估了梁的等效剛度。

    表1 梁的材料參數(shù)和幾何尺寸

    表2 x =L /3處的無量綱化撓度和轉(zhuǎn)角

    5.4 結(jié)果討論

    為了方便討論,在接下來的結(jié)果分析中將給出4種模型。(1) SF模型,同時(shí)考慮撓曲電效應(yīng)和應(yīng)變梯度彈性項(xiàng),此為當(dāng)前的研究模型;(2) FL模型,只考慮撓曲電效應(yīng);(3) SG模型,只考慮應(yīng)變梯度彈性項(xiàng)。(4)NF模型:既不考慮撓曲電效應(yīng),也不考慮應(yīng)變梯度彈性項(xiàng)。

    圖3為四種模型的撓度曲線。首先對(duì)比FL,SG和NF模型的撓度值,可以看出撓曲電效應(yīng)可以明顯地增大梁的撓度,即考慮撓曲電效應(yīng)后梁的等效彎曲剛度減小,而考慮應(yīng)變梯度彈性項(xiàng)后,梁的撓度減小,即梁的等效剛度增加。在當(dāng)前結(jié)構(gòu)尺寸設(shè)定下,撓曲電效應(yīng)和應(yīng)變梯度彈性項(xiàng)對(duì)梁的變形的影響是相反的。這一點(diǎn)可以從梁的等效彎曲剛度k2看出,撓曲電效應(yīng)會(huì)減小梁的等效彎曲剛度,而應(yīng)變梯度彈性項(xiàng)則會(huì)增大梁的等效彎曲剛度。同時(shí),因?yàn)閾锨娦?yīng),外加電壓相當(dāng)于在梁的兩端施加了等效彎矩M*,這相當(dāng)于增加了外加載荷,從而使梁的撓度增大。再來看SF模型,相比撓曲電效應(yīng)和應(yīng)變梯度彈性項(xiàng)皆不考慮的NF模型,SF模型梁的撓度有所增大,但撓度值還是小于FL模型,這說明應(yīng)變梯度彈性項(xiàng)對(duì)撓曲電效應(yīng)存在一定的抑制效果,但在當(dāng)前模型尺寸下?lián)锨娦?yīng)還是起主導(dǎo)作用的。最后,圖3的數(shù)值結(jié)果與理論結(jié)果擬合良好,再次驗(yàn)證了本文有限單元法的有效性。

    圖2 兩種模型的理論與數(shù)值結(jié)果對(duì)比

    圖3 四種模型的簡支梁的撓度

    令簡支梁的材料參數(shù)不變,幾何尺寸和外加集中力F進(jìn)行等比例縮放,計(jì)算可得四種模型歸一化的最大撓度,如圖4所示。wN F為NF模型下梁的最大撓度??梢钥闯?SF模型歸一化的最大撓度值隨著結(jié)構(gòu)尺寸的減小先增大后減小,存在極值點(diǎn)。而FL模型和SG模型歸一化的最大撓度值則分別隨著結(jié)構(gòu)尺寸的減小單調(diào)遞增和單調(diào)遞減。對(duì)于傳統(tǒng)的NF模型,隨著尺寸的減小,其歸一化的最大撓度值不變。從圖4可以看出與圖3一致的結(jié)論,即考慮撓曲電效應(yīng)之后,梁的撓度會(huì)增大,而考慮應(yīng)變梯度彈性項(xiàng)后,梁的撓度會(huì)減小。同時(shí)圖4還可以看出,撓曲電效應(yīng)和應(yīng)變梯度效應(yīng)都存在尺寸效應(yīng),結(jié)構(gòu)尺寸越小,其作用越明顯。而對(duì)于SF和FL模型,當(dāng)結(jié)構(gòu)尺寸比較大時(shí),SF模型的歸一化的最大撓度值是趨向于FL模型的,說明此時(shí)撓曲電效應(yīng)起主導(dǎo)作用,應(yīng)變梯度彈性項(xiàng)的作用不明顯。隨著幾何尺寸減小,FL模型的撓度值會(huì)迅速增大,而SF模型會(huì)出現(xiàn)拐點(diǎn),說明此時(shí)應(yīng)變梯度彈性項(xiàng)對(duì)撓曲電效應(yīng)的抑制作用增強(qiáng)。當(dāng)幾何尺寸與材料尺度參數(shù)相當(dāng)時(shí),兩者的作用幾乎相互抵消,所以此時(shí)SF模型梁的最大撓度趨近NF模型。

    圖4 四種模型簡支梁的最大歸一化撓度

    圖3~圖5為納米簡支梁的力學(xué)響應(yīng),接下來將討論其電學(xué)響應(yīng)。圖6為外加電壓為0,僅有外加集中力作用下時(shí),SF模型在x=L/3處沿厚度方向的電勢。由電勢的表達(dá)式(8)可知,沿梁厚度方向的電勢為二次函數(shù),外加電壓為0時(shí),梁的上下表面的電勢為0,與圖6曲線一致。如圖6所示,隨著撓曲電系數(shù)的增大,梁中產(chǎn)生的電勢也會(huì)隨之增大。f1133=5 V時(shí)的電勢相較于f1133=0時(shí)增大了2.6%,f1133=15 V時(shí)的電勢相較于f1133=10 V時(shí)增大了17.4%,這說明撓曲電系數(shù)越大,電勢增長越快。因此在撓曲電效應(yīng)的應(yīng)用中,應(yīng)當(dāng)優(yōu)先選擇撓曲電系數(shù)大的材料,以保證獲得更高的能量轉(zhuǎn)化效率。

    最后,圖7給出了在只施加集中力時(shí),SF,F(xiàn)L和NF三種模型在不同撓曲電系數(shù)情況下的歸一化的能量效率deff/dB。能量效率定義為施加在梁上的每單位力載荷產(chǎn)生的電荷量,表示為deff=Q/F,其中Q=?P(x)dxdy為產(chǎn)生的電荷量,F為外加的集中力。dB=1.05×10-12為NF模型的能量效率。如圖7所示,隨著梁幾何尺寸的減小,FL模型的能量效率持續(xù)增高,這驗(yàn)證了撓曲電效應(yīng)存在尺寸依賴性,而SF模型的能量效率則是先增大后減小,也存在與圖4和圖5中相似的極值點(diǎn),這再次說明應(yīng)變梯度彈性項(xiàng)對(duì)撓曲電效應(yīng)具有抑制作用,并且尺寸越小,其抑制作用越明顯。因此,在納米尺度的研究中,應(yīng)變梯度彈性項(xiàng)的影響非常顯著,不可忽略。當(dāng)梁的尺寸足夠大時(shí),SF和FL模型的能量效率趨于相同。此外,對(duì)比不同撓曲電系數(shù)的能量效率曲線可以發(fā)現(xiàn),更大的撓曲電系數(shù)有助于獲得更高的能量效率,與圖6的結(jié)論一致。

    圖5 不同電壓下SF模型的最大撓度

    圖6 x =L /3處不同撓曲電系數(shù)的電勢

    圖7 不同模型的能量效率

    6 結(jié) 論

    本文基于全應(yīng)變梯度彈性理論,利用能量法推得了考慮高階項(xiàng)的納米歐拉梁的控制方程和邊界條件。并且在傳統(tǒng)歐拉梁單元的基礎(chǔ)上,創(chuàng)建了一種新的滿足C2弱連續(xù)的歐拉納米梁單元。這種新的梁單元包含三個(gè)獨(dú)立的材料尺度參數(shù),有2個(gè)節(jié)點(diǎn),6個(gè)自由度。而后,從新的歐拉梁控制方程出發(fā),根據(jù)有限元理論推得弱形式方程,繼而推得梁的剛度方程。最后,用有限單元法求解本文算例簡支梁的力電響應(yīng)。通過分析撓度、電勢和能量效率,得出如下結(jié)論。(1) 撓曲電效應(yīng)和應(yīng)變梯度彈性都存在尺寸依賴性,尤其是當(dāng)結(jié)構(gòu)的特征尺寸接近材料尺度參數(shù)時(shí),二者的作用效果都將變得非常顯著。(2) 外加電壓相當(dāng)于在梁的兩端施加等效彎矩,電壓的正負(fù)決定等效彎矩的方向。(3) 應(yīng)變梯度彈性項(xiàng)對(duì)撓曲電效應(yīng)以及能量的轉(zhuǎn)化存在一定的抑制作用,因此在撓曲電問題中,其影響不可忽略。

    猜你喜歡
    效應(yīng)模型
    一半模型
    鈾對(duì)大型溞的急性毒性效應(yīng)
    懶馬效應(yīng)
    場景效應(yīng)
    重要模型『一線三等角』
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    應(yīng)變效應(yīng)及其應(yīng)用
    3D打印中的模型分割與打包
    FLUKA幾何模型到CAD幾何模型轉(zhuǎn)換方法初步研究
    偶像效應(yīng)
    国产在线视频一区二区| 黄色配什么色好看| 看十八女毛片水多多多| 国产成人精品久久久久久| 丝袜人妻中文字幕| 国产成人91sexporn| 久久久久久久精品精品| 咕卡用的链子| 青春草国产在线视频| 日韩制服丝袜自拍偷拍| 亚洲欧美成人精品一区二区| 久久久久久久久久久久大奶| 美女国产视频在线观看| 制服人妻中文乱码| 热99国产精品久久久久久7| 中文字幕精品免费在线观看视频 | 夫妻性生交免费视频一级片| 人人妻人人爽人人添夜夜欢视频| 久久狼人影院| 亚洲精品乱久久久久久| 成人亚洲欧美一区二区av| 久久亚洲国产成人精品v| av视频免费观看在线观看| 中国美白少妇内射xxxbb| h视频一区二区三区| av在线老鸭窝| av线在线观看网站| 香蕉精品网在线| 性色av一级| 亚洲国产色片| 97人妻天天添夜夜摸| 国语对白做爰xxxⅹ性视频网站| 亚洲精品久久久久久婷婷小说| 午夜福利在线观看免费完整高清在| 亚洲成色77777| 亚洲成av片中文字幕在线观看 | 又大又黄又爽视频免费| 在线天堂最新版资源| 亚洲国产日韩一区二区| 国产深夜福利视频在线观看| 两个人看的免费小视频| 全区人妻精品视频| 久久亚洲国产成人精品v| 秋霞伦理黄片| 精品人妻一区二区三区麻豆| 成人二区视频| 国产亚洲精品久久久com| 热re99久久国产66热| 青青草视频在线视频观看| 制服人妻中文乱码| 一级a做视频免费观看| 草草在线视频免费看| 中文欧美无线码| 国产一区二区三区av在线| 国产 一区精品| 大片电影免费在线观看免费| 婷婷色综合www| 黄网站色视频无遮挡免费观看| 三上悠亚av全集在线观看| 国产高清三级在线| 99国产综合亚洲精品| 最新中文字幕久久久久| 日日爽夜夜爽网站| 婷婷成人精品国产| 亚洲三级黄色毛片| 国产精品久久久久久久电影| 国产麻豆69| 色视频在线一区二区三区| 女性被躁到高潮视频| 飞空精品影院首页| 日韩在线高清观看一区二区三区| 观看av在线不卡| 夫妻午夜视频| 美女国产高潮福利片在线看| 免费黄频网站在线观看国产| 美女脱内裤让男人舔精品视频| 少妇被粗大猛烈的视频| 中文字幕最新亚洲高清| 韩国av在线不卡| 各种免费的搞黄视频| 午夜老司机福利剧场| 国产精品久久久久久久久免| 国产高清不卡午夜福利| 日本色播在线视频| 天堂8中文在线网| 看免费av毛片| 精品亚洲成a人片在线观看| 中国美白少妇内射xxxbb| 纵有疾风起免费观看全集完整版| 精品久久蜜臀av无| 一本大道久久a久久精品| 精品国产一区二区三区四区第35| 国产一区二区激情短视频 | 国产免费现黄频在线看| 波多野结衣一区麻豆| 国产精品三级大全| 下体分泌物呈黄色| 精品一品国产午夜福利视频| 天堂中文最新版在线下载| 人人妻人人澡人人看| 9色porny在线观看| 色婷婷久久久亚洲欧美| 51国产日韩欧美| 国产成人a∨麻豆精品| 高清不卡的av网站| 久久久久久久精品精品| 久久 成人 亚洲| 好男人视频免费观看在线| 亚洲欧美中文字幕日韩二区| 久久精品久久久久久久性| 国产精品国产三级国产专区5o| 欧美 亚洲 国产 日韩一| 国产麻豆69| 美女主播在线视频| 精品国产乱码久久久久久小说| 制服诱惑二区| 少妇人妻久久综合中文| 十分钟在线观看高清视频www| 成人国产av品久久久| 少妇高潮的动态图| 国产精品不卡视频一区二区| 色5月婷婷丁香| 国产不卡av网站在线观看| 久久毛片免费看一区二区三区| 国产免费又黄又爽又色| 男女无遮挡免费网站观看| 少妇的逼水好多| 人妻一区二区av| 国产探花极品一区二区| videosex国产| 亚洲成色77777| 夫妻性生交免费视频一级片| 全区人妻精品视频| 插逼视频在线观看| 狂野欧美激情性bbbbbb| 丝袜美足系列| av电影中文网址| 黄色毛片三级朝国网站| 成人黄色视频免费在线看| 欧美日韩一区二区视频在线观看视频在线| 国产欧美亚洲国产| 欧美激情极品国产一区二区三区 | 国产 一区精品| 国产欧美日韩综合在线一区二区| 国产成人91sexporn| 精品午夜福利在线看| 大片电影免费在线观看免费| 丰满迷人的少妇在线观看| 久久青草综合色| 人人澡人人妻人| 一区二区三区四区激情视频| 热re99久久国产66热| 亚洲av免费高清在线观看| 青春草国产在线视频| 欧美日韩亚洲高清精品| 自线自在国产av| 国产精品国产三级国产av玫瑰| 欧美性感艳星| 国产免费一区二区三区四区乱码| 久久午夜福利片| videossex国产| 一本久久精品| 成年人午夜在线观看视频| 麻豆乱淫一区二区| 国产无遮挡羞羞视频在线观看| 久久久久久人妻| av线在线观看网站| 又黄又爽又刺激的免费视频.| 久久久欧美国产精品| 男女午夜视频在线观看 | 国产av码专区亚洲av| 亚洲精品视频女| 久久精品国产自在天天线| 午夜久久久在线观看| 亚洲精品日本国产第一区| 青春草视频在线免费观看| 一级,二级,三级黄色视频| 国产精品一国产av| 久久 成人 亚洲| 亚洲美女搞黄在线观看| 女性生殖器流出的白浆| 天天躁夜夜躁狠狠久久av| 久久精品aⅴ一区二区三区四区 | 国产 精品1| 丝袜在线中文字幕| 丝瓜视频免费看黄片| 伦精品一区二区三区| 男女免费视频国产| 丝袜人妻中文字幕| 五月伊人婷婷丁香| 日韩中文字幕视频在线看片| 亚洲内射少妇av| 黄色怎么调成土黄色| 在线 av 中文字幕| 在线观看www视频免费| 亚洲国产毛片av蜜桃av| 午夜日本视频在线| 日本av手机在线免费观看| 国产乱来视频区| 九九爱精品视频在线观看| av卡一久久| 七月丁香在线播放| 一级毛片 在线播放| 亚洲一区二区三区欧美精品| 国产黄色视频一区二区在线观看| 国产片特级美女逼逼视频| 一区二区三区乱码不卡18| 国产精品嫩草影院av在线观看| 国产日韩欧美视频二区| 在线亚洲精品国产二区图片欧美| 十八禁高潮呻吟视频| 成年女人在线观看亚洲视频| 两个人免费观看高清视频| 免费大片黄手机在线观看| 久久99一区二区三区| 一区二区三区精品91| 亚洲综合色惰| 天美传媒精品一区二区| 午夜福利影视在线免费观看| 久久精品熟女亚洲av麻豆精品| 久久久久久久久久久久大奶| 天堂8中文在线网| 亚洲精品国产av蜜桃| 伦理电影大哥的女人| 色婷婷av一区二区三区视频| 老司机亚洲免费影院| 午夜免费观看性视频| 亚洲av国产av综合av卡| 久久人人爽av亚洲精品天堂| 狠狠精品人妻久久久久久综合| 国产高清不卡午夜福利| 多毛熟女@视频| 久久精品久久久久久久性| 日韩在线高清观看一区二区三区| 亚洲人成77777在线视频| 国产精品久久久久久精品古装| 黄色 视频免费看| 亚洲精品中文字幕在线视频| 久久久久久久亚洲中文字幕| 亚洲,欧美,日韩| 91精品三级在线观看| 中文字幕另类日韩欧美亚洲嫩草| 欧美日韩成人在线一区二区| 26uuu在线亚洲综合色| 色哟哟·www| 18禁在线无遮挡免费观看视频| 久久久久网色| 国产精品久久久久久精品电影小说| 欧美丝袜亚洲另类| 久久国产精品男人的天堂亚洲 | 国产福利在线免费观看视频| 欧美日韩精品成人综合77777| 制服人妻中文乱码| av黄色大香蕉| 丝袜人妻中文字幕| 女性被躁到高潮视频| 国产毛片在线视频| 精品一区二区免费观看| 国产精品99久久99久久久不卡 | 午夜91福利影院| 日本91视频免费播放| 满18在线观看网站| 国产综合精华液| 久久狼人影院| 两性夫妻黄色片 | 狠狠精品人妻久久久久久综合| 成人亚洲欧美一区二区av| 亚洲av福利一区| 韩国av在线不卡| 国产一区二区激情短视频 | 少妇被粗大猛烈的视频| 人妻系列 视频| 90打野战视频偷拍视频| 少妇的逼好多水| 另类精品久久| 夫妻性生交免费视频一级片| 狂野欧美激情性xxxx在线观看| 婷婷成人精品国产| 国产熟女欧美一区二区| 在线观看国产h片| 丝瓜视频免费看黄片| 亚洲国产精品一区二区三区在线| 欧美日韩成人在线一区二区| 亚洲欧美色中文字幕在线| 久久这里有精品视频免费| 亚洲国产精品国产精品| 赤兔流量卡办理| 纵有疾风起免费观看全集完整版| 日韩av不卡免费在线播放| 亚洲美女视频黄频| 国产精品一区二区在线不卡| 亚洲国产色片| 亚洲,欧美精品.| 日韩不卡一区二区三区视频在线| 免费观看无遮挡的男女| 亚洲国产精品国产精品| 中文字幕人妻丝袜制服| 国产成人精品婷婷| 97人妻天天添夜夜摸| 精品人妻在线不人妻| 18在线观看网站| 一级a做视频免费观看| a级毛片黄视频| 天天躁夜夜躁狠狠躁躁| 亚洲人成网站在线观看播放| 插逼视频在线观看| 少妇人妻 视频| 两个人看的免费小视频| 国产成人a∨麻豆精品| 国产在视频线精品| 激情视频va一区二区三区| 熟妇人妻不卡中文字幕| 熟妇人妻不卡中文字幕| 国产探花极品一区二区| 久久99蜜桃精品久久| 欧美另类一区| 97精品久久久久久久久久精品| 男女午夜视频在线观看 | 亚洲精品久久午夜乱码| 最近手机中文字幕大全| 欧美精品亚洲一区二区| 久久久精品94久久精品| 王馨瑶露胸无遮挡在线观看| 一边摸一边做爽爽视频免费| tube8黄色片| 精品熟女少妇av免费看| 精品一区二区三区四区五区乱码 | 人妻一区二区av| 亚洲一码二码三码区别大吗| 18禁裸乳无遮挡动漫免费视频| 美女内射精品一级片tv| 亚洲精品国产av蜜桃| 亚洲国产日韩一区二区| 超色免费av| 国产片内射在线| 51国产日韩欧美| 久久毛片免费看一区二区三区| 欧美日韩视频精品一区| av女优亚洲男人天堂| 少妇高潮的动态图| 国产精品一区二区在线不卡| 成人综合一区亚洲| 日韩不卡一区二区三区视频在线| a级毛色黄片| 美女脱内裤让男人舔精品视频| 黑人欧美特级aaaaaa片| 日本欧美视频一区| 中文字幕人妻熟女乱码| 久久久精品区二区三区| 国产成人精品久久久久久| 久久青草综合色| 亚洲第一av免费看| 亚洲成国产人片在线观看| 婷婷色av中文字幕| 亚洲精品日韩在线中文字幕| 老司机亚洲免费影院| 蜜桃在线观看..| 久久久国产一区二区| 久久人人97超碰香蕉20202| 亚洲精品日韩在线中文字幕| 大香蕉97超碰在线| 波野结衣二区三区在线| 久久综合国产亚洲精品| 久久久精品免费免费高清| 777米奇影视久久| 狠狠精品人妻久久久久久综合| 亚洲精品乱码久久久久久按摩| 国语对白做爰xxxⅹ性视频网站| 亚洲成av片中文字幕在线观看 | 男的添女的下面高潮视频| 我要看黄色一级片免费的| av有码第一页| 丰满迷人的少妇在线观看| 国产成人精品一,二区| 日本av免费视频播放| 99久国产av精品国产电影| 插逼视频在线观看| 免费观看无遮挡的男女| 亚洲,欧美精品.| 国产男人的电影天堂91| av线在线观看网站| 国产 一区精品| 久久久久网色| 99视频精品全部免费 在线| 日韩成人av中文字幕在线观看| 秋霞伦理黄片| 精品一区二区免费观看| 午夜精品国产一区二区电影| 黄片无遮挡物在线观看| 亚洲精品日本国产第一区| 国产日韩欧美在线精品| 国产精品秋霞免费鲁丝片| 久久久久久人人人人人| 亚洲精品久久久久久婷婷小说| 国产精品熟女久久久久浪| 90打野战视频偷拍视频| 亚洲 欧美一区二区三区| 日本欧美视频一区| 一边亲一边摸免费视频| 久久久久久久亚洲中文字幕| 国产精品一国产av| 久久精品夜色国产| 亚洲欧美一区二区三区黑人 | 一级毛片我不卡| 欧美激情 高清一区二区三区| 国产精品人妻久久久影院| 天天躁夜夜躁狠狠躁躁| 如日韩欧美国产精品一区二区三区| av女优亚洲男人天堂| 国产av码专区亚洲av| 一二三四中文在线观看免费高清| 国产毛片在线视频| 国产高清国产精品国产三级| h视频一区二区三区| 成人国产av品久久久| 中国美白少妇内射xxxbb| av在线观看视频网站免费| 夫妻午夜视频| 91国产中文字幕| 岛国毛片在线播放| 视频区图区小说| 午夜福利视频精品| 国产高清三级在线| 一区二区三区精品91| 亚洲婷婷狠狠爱综合网| 欧美精品高潮呻吟av久久| 久久久久久久久久人人人人人人| 精品久久久久久电影网| 好男人视频免费观看在线| 免费人成在线观看视频色| 亚洲精品日韩在线中文字幕| 免费观看性生交大片5| 亚洲美女搞黄在线观看| 天天操日日干夜夜撸| 欧美少妇被猛烈插入视频| 欧美成人精品欧美一级黄| 成人亚洲欧美一区二区av| 国产国语露脸激情在线看| 免费高清在线观看视频在线观看| 中国三级夫妇交换| 一级毛片我不卡| 精品一品国产午夜福利视频| 久久久久精品久久久久真实原创| 国产精品国产三级国产专区5o| 满18在线观看网站| 国产在视频线精品| 五月开心婷婷网| 亚洲精品国产色婷婷电影| 巨乳人妻的诱惑在线观看| 美女国产视频在线观看| 纯流量卡能插随身wifi吗| 国产极品天堂在线| 国产黄色视频一区二区在线观看| 免费观看在线日韩| 午夜老司机福利剧场| 精品酒店卫生间| 色视频在线一区二区三区| 少妇被粗大的猛进出69影院 | 国产麻豆69| 韩国高清视频一区二区三区| 亚洲一区二区三区欧美精品| www日本在线高清视频| 全区人妻精品视频| 嫩草影院入口| 欧美激情极品国产一区二区三区 | 日韩欧美一区视频在线观看| 亚洲国产精品专区欧美| 婷婷色av中文字幕| 香蕉丝袜av| 午夜免费鲁丝| 午夜视频国产福利| 一本—道久久a久久精品蜜桃钙片| 欧美少妇被猛烈插入视频| 亚洲精品456在线播放app| 只有这里有精品99| 国产精品99久久99久久久不卡 | 99视频精品全部免费 在线| 老司机影院成人| 日产精品乱码卡一卡2卡三| 天天躁夜夜躁狠狠躁躁| 五月玫瑰六月丁香| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 亚洲 欧美一区二区三区| 色网站视频免费| 中文乱码字字幕精品一区二区三区| 2022亚洲国产成人精品| 老熟女久久久| 男女边吃奶边做爰视频| 欧美激情极品国产一区二区三区 | 最新的欧美精品一区二区| 女人被躁到高潮嗷嗷叫费观| 各种免费的搞黄视频| 亚洲精品一二三| 91在线精品国自产拍蜜月| 狠狠婷婷综合久久久久久88av| a 毛片基地| 美女主播在线视频| 午夜影院在线不卡| 观看美女的网站| 99国产综合亚洲精品| 少妇人妻精品综合一区二区| 久久精品久久久久久噜噜老黄| 国产免费又黄又爽又色| 国产综合精华液| 在线天堂最新版资源| 久久精品夜色国产| 搡女人真爽免费视频火全软件| 亚洲国产毛片av蜜桃av| 美女xxoo啪啪120秒动态图| 人人妻人人澡人人爽人人夜夜| 亚洲精品日本国产第一区| av播播在线观看一区| 在线观看人妻少妇| 成人影院久久| 免费av不卡在线播放| 下体分泌物呈黄色| 最新中文字幕久久久久| 男女无遮挡免费网站观看| 久久精品久久精品一区二区三区| 国产福利在线免费观看视频| 婷婷成人精品国产| 国产男人的电影天堂91| 亚洲久久久国产精品| 美国免费a级毛片| 午夜日本视频在线| 日日啪夜夜爽| 亚洲欧洲精品一区二区精品久久久 | 各种免费的搞黄视频| 97超碰精品成人国产| 欧美日韩国产mv在线观看视频| 久久久亚洲精品成人影院| 中文字幕亚洲精品专区| 国产欧美日韩综合在线一区二区| 女人被躁到高潮嗷嗷叫费观| 免费播放大片免费观看视频在线观看| 国产精品久久久久久精品古装| 国产精品偷伦视频观看了| 国产亚洲一区二区精品| 国产熟女午夜一区二区三区| 韩国高清视频一区二区三区| 日韩av不卡免费在线播放| 最近的中文字幕免费完整| 天天操日日干夜夜撸| 中文字幕精品免费在线观看视频 | 春色校园在线视频观看| 精品熟女少妇av免费看| 精品一区二区三区视频在线| 成人毛片a级毛片在线播放| 中文乱码字字幕精品一区二区三区| 欧美精品一区二区大全| 亚洲欧美中文字幕日韩二区| 中文天堂在线官网| 少妇人妻精品综合一区二区| 国产精品国产三级专区第一集| 亚洲精品中文字幕在线视频| 99香蕉大伊视频| 精品99又大又爽又粗少妇毛片| av片东京热男人的天堂| 久久久久久久亚洲中文字幕| 男女边吃奶边做爰视频| 97超碰精品成人国产| 综合色丁香网| 色94色欧美一区二区| 美国免费a级毛片| 一级片免费观看大全| 99热全是精品| 国产 一区精品| 欧美97在线视频| 卡戴珊不雅视频在线播放| 母亲3免费完整高清在线观看 | 在线观看免费视频网站a站| 中文字幕av电影在线播放| 日韩av不卡免费在线播放| 乱人伦中国视频| 久久综合国产亚洲精品| 制服丝袜香蕉在线| 久久久久久人妻| 男女国产视频网站| 你懂的网址亚洲精品在线观看| 免费看av在线观看网站| 2018国产大陆天天弄谢| 国产精品一区二区在线观看99| 亚洲欧美一区二区三区国产| 精品久久久久久电影网| 国产精品不卡视频一区二区| 热99国产精品久久久久久7| 99香蕉大伊视频| 国产精品国产av在线观看| 一区二区三区乱码不卡18| 一区在线观看完整版| 人妻 亚洲 视频| 水蜜桃什么品种好| 日本vs欧美在线观看视频| 在线精品无人区一区二区三| 亚洲av免费高清在线观看| 视频区图区小说| 久久av网站| 婷婷色综合www| av福利片在线| 赤兔流量卡办理| 欧美 亚洲 国产 日韩一| 婷婷色综合大香蕉| 免费看不卡的av| 欧美国产精品va在线观看不卡| 街头女战士在线观看网站| 在线观看免费视频网站a站| 美女内射精品一级片tv| 如日韩欧美国产精品一区二区三区| 在线观看免费视频网站a站| 好男人视频免费观看在线| 欧美日韩视频精品一区| 午夜福利乱码中文字幕| 午夜福利影视在线免费观看| 欧美bdsm另类| 国产不卡av网站在线观看|