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

    一種預(yù)變形周期剛架結(jié)構(gòu)振動帶隙特性分析

    2022-10-09 09:02:58郭振坤溫佳琦張軍董挺
    關(guān)鍵詞:單胞架結(jié)構(gòu)端點(diǎn)

    郭振坤, 溫佳琦, 張軍, 董挺

    (1.北京建筑大學(xué) 機(jī)電與車輛工程學(xué)院, 北京 102616; 2.北京理工大學(xué) 機(jī)械與車輛學(xué)院, 北京 100081)

    工程減振問題在現(xiàn)代工程研究中十分常見[1],為滿足工業(yè)和運(yùn)輸業(yè)高速發(fā)展的需求,各種機(jī)械裝置、精密儀器受到振動的影響愈加突出[2],同時(shí)為加快國防現(xiàn)代化的腳步,要求飛機(jī)、坦克、艦船等能夠在最為惡劣的環(huán)境中工作[3],其工作性能都與減振技術(shù)密切相關(guān)。因此,關(guān)于如何減振的研究一直在持續(xù),最行之有效的方法是究其原因?qū)ΠY下藥,學(xué)者們做出了很多有價(jià)值的研究[4-5]。

    周期結(jié)構(gòu)因?yàn)槠洫?dú)特的減振優(yōu)勢,在許多工程領(lǐng)域越來越受到關(guān)注,例如聲子晶體和光子晶體。Kushwaha等[6]利用平面波展開法研究了聲子晶體帶隙,即可以在某些頻率范圍內(nèi),阻止彈性波或聲波的傳播。這些阻止波傳播的頻率范圍稱為頻率禁帶,而那些允許波傳播的頻率范圍稱為頻率通帶;當(dāng)振動或波的頻率處于周期結(jié)構(gòu)的頻率帶隙內(nèi)時(shí),將無法被傳遞。這樣的禁帶特性吸引了國內(nèi)外眾多學(xué)者進(jìn)行研究,溫激鴻等[7]利用周期彈簧振子結(jié)構(gòu)簡化了聲子晶體,研究了聲子晶體的振動特性,此后眾多學(xué)者在此基礎(chǔ)上研究了多種結(jié)構(gòu)的振動帶隙并將其應(yīng)用于振動抑制和減振降噪等工程領(lǐng)域。

    對于周期性聲子晶體結(jié)構(gòu),通??梢酝ㄟ^改變部件的拓?fù)浣Y(jié)構(gòu)、剛度、密度或在某些局部增加共振器來優(yōu)化振動和彈性波的傳播特性。Yu等[8]利用傳遞矩陣法(TMM)研究了局域共振超材料梁的彎曲波衰減特性,其中局域振子是由軟橡膠環(huán)和銅環(huán)設(shè)計(jì)而成。在此基礎(chǔ)上,為提高工程適用能力和減振效果,針對如何根據(jù)實(shí)際工程環(huán)境靈活調(diào)整帶隙位置和寬度,學(xué)者們完成了大量的設(shè)計(jì)和研究工作[9-11]。

    周期剛架結(jié)構(gòu)作為在工程建筑、機(jī)械中常見的結(jié)構(gòu)也受到了廣泛的關(guān)注。吳志靜等[12]采用譜元法研究了周期桁架結(jié)構(gòu)的振動帶隙特性,這種方法具有在中高頻區(qū)域計(jì)算誤差小、劃分單元少且計(jì)算效率高等優(yōu)點(diǎn)?;谠摲椒?,Zuo等[13]通過進(jìn)一步改變結(jié)構(gòu)構(gòu)型和幾何參數(shù)等對周期剛架結(jié)構(gòu)進(jìn)行設(shè)計(jì)優(yōu)化,并討論了影響帶隙特性的因素。一些學(xué)者還提出了較為新穎的超材料結(jié)構(gòu)。蔡昌琦等[14]利用準(zhǔn)零剛度理論設(shè)計(jì)出一種具有低頻帶隙區(qū)間嵌入式準(zhǔn)零剛度振子的新型超材料梁。Muhammad等[15]提出了一種具有超寬帶隙的周期性和非周期性3D復(fù)合超結(jié)構(gòu),采用聚合物外殼的設(shè)計(jì)和球形/圓柱形鋼塊的嵌入,提高共振的動態(tài)特性和機(jī)械性能,使其擁有非常明顯的低頻帶隙。進(jìn)一步,一些研究者發(fā)現(xiàn)還可以通過結(jié)構(gòu)元件的特定預(yù)變形來控制周期結(jié)構(gòu)的振動和波傳播特性。Mellmann等[16]提出一種預(yù)變形網(wǎng)格結(jié)構(gòu),計(jì)算了波在這種網(wǎng)格結(jié)構(gòu)中的傳播特性,發(fā)現(xiàn)可以通過結(jié)構(gòu)預(yù)變形靈活地調(diào)控波帶隙的位置和寬度。受此啟發(fā),本文設(shè)計(jì)了一種預(yù)變形周期剛架結(jié)構(gòu),通過內(nèi)部構(gòu)件的預(yù)變形實(shí)現(xiàn)剛度的調(diào)節(jié),進(jìn)而控制波帶隙特性。此外,利用有限元方法進(jìn)行了模態(tài)分析和振動透射率的計(jì)算,驗(yàn)證譜元法的正確性。

    1 動力學(xué)模型

    圖1(a)所示為左端固支的預(yù)變形周期剛架結(jié)構(gòu),該結(jié)構(gòu)包含20個(gè)單胞,每個(gè)單胞由11根梁/桿結(jié)構(gòu)組成,由圖1(b)看出,該單胞結(jié)構(gòu)是通過將X型結(jié)構(gòu)進(jìn)行折斷所形成,變形角用α表示。

    圖1 預(yù)變形周期剛架結(jié)構(gòu)示意Fig.1 Schematic diagram of the pre-deformed periodic rigid frame structure

    1.1 桿單元求解

    軸向振動桿件的偏微分運(yùn)動方程表示為:

    (1)

    式中:u(x,t)為軸向位移;E和ρ分別是材料彈性模量和密度;A是桿件橫截面積。式(1)的解可以假設(shè)成:

    (2)

    式中:Un(x,ωn)為微元的軸向位移;ωn為桿的圓頻率。

    將式(2)代入(1)中可以得到對于每個(gè)離散的頻率ωn的特征值問題:

    EAU″+ω2ρAU=0

    (3)

    假設(shè)其通解為:

    U(x)=be-ik(ω)x

    (4)

    將式(4)代入式(3)中得到色散關(guān)系:

    (5)

    即可得到2個(gè)實(shí)根k1和k2:

    k1=-k2=kL

    (6)

    將式(6)代回式(4)可得:

    U(x)=b1e-ikLx+b2eikLx

    (7)

    則對于長度為L的桿可以用矩陣形式來表示:

    U(x)=e(x,ω)b

    (8)

    用U1和U2表示桿件兩端的節(jié)點(diǎn)位移,則有U1=U(0)和U2=U(L),可以記為:

    (9)

    (10)

    令dR=HR(ω)b,則有:

    (11)

    用N1和N2表示桿件兩端節(jié)點(diǎn)力,桿的節(jié)點(diǎn)力函數(shù)為N(x)。

    整理可得:

    (12)

    式中:NR1(x,ω)=csc(kLL)sin(kL(L-x));NR2(x,ω)=csc(kLL)sin(kLx)。

    (13)

    軸向力分量與U(x)有關(guān),即:

    N(x)=EAU′(x)

    (14)

    將節(jié)點(diǎn)位移代入式(4),可得:

    U(x)=NR1U1+NR2U2

    (15)

    類似地,軸向力也可以寫成譜表達(dá)式,表示為:

    (16)

    節(jié)點(diǎn)軸向力和由材料強(qiáng)度定義的內(nèi)部力之間的平衡條件有:

    (17)

    則有:

    SR(ω)dR=fc(ω)

    (18)

    根據(jù)軸向力與位移的關(guān)系,可以得到頻域上位移與節(jié)點(diǎn)力的關(guān)系:

    (19)

    式中SR(ω)為桿單元的剛度矩陣,可以表示為:

    (20)

    (21)

    整理可得:

    (22)

    1.2 梁單元求解

    類似地,可以推導(dǎo)出梁單元位移與節(jié)點(diǎn)力的關(guān)系。Timoshenko梁的橫向振動需要考慮轉(zhuǎn)動慣量和剪切變形,用譜表達(dá)式對橫向位移w(x,t)和轉(zhuǎn)角θ(x,t)分別進(jìn)行表示:

    (23)

    (24)

    梁的橫向振動偏微分方程為:

    (25)

    (26)

    式中:κ的取值與截面形狀有關(guān);G是剪切模量;A是橫截面積;I是截面慣性矩。

    將式(23)和(24)代入偏微分方程(25)、(26)中得:

    (27)

    (28)

    式(27)、(28)的通解分別假設(shè)為:W(x)=ae-ik(ω)x和Θ(x)=βeik(ω)x,可得:

    (29)

    則有:

    (30)

    求解方程(30)可以得到4個(gè)根:

    k1=-k2=

    (31)

    k2=-k4=

    (32)

    代入式(29)中可得:

    (33)

    (34)

    由此,式(27)和(28)的通解可得:

    W(x,ω)=a1e-iktx+a2eiktx+a3e-ikex+a4eikex=

    ew(x,ω)a

    (35)

    Θ(x,ω)=β1a1e-iktx+β2a2eiktx+β3a3e-ikex+

    β4a4eikex=eθ(x,ω)a

    (36)

    其中:

    (37)

    (38)

    eθ(x,ω)=ew(x,ω)B(ω)

    (39)

    梁的兩端節(jié)點(diǎn)位移和轉(zhuǎn)角可以表示為:

    (40)

    將式(35)和(36)代入式(40)可得:

    (41)

    式中:

    (42)

    (43)

    將式(41)代入到式(35)和(36)中可得:

    (44)

    (45)

    梁的節(jié)點(diǎn)力為:

    (46)

    將節(jié)點(diǎn)力表示為譜表達(dá)形式:

    (47)

    (48)

    es(x,ω)=[-ktrte-iktx-ktrteiktx-keree-ikex-kereeikex]。

    則有端點(diǎn)位移與力的關(guān)系為:

    (49)

    即:

    SB(ω)dB=f

    (50)

    式中SB(ω)為Timoshenko梁的單元剛度矩陣:

    (51)

    其中:

    最終推導(dǎo)出總的力與位移的關(guān)系為:

    FB=SB(ω)UB

    (52)

    式中:FB和UB為節(jié)點(diǎn)力矩陣和節(jié)點(diǎn)位移矩陣;SB(ω)為剛度矩陣。

    1.3 整體結(jié)構(gòu)組裝

    利用坐標(biāo)轉(zhuǎn)換矩陣,將上述梁、桿結(jié)構(gòu)的節(jié)點(diǎn)力和節(jié)點(diǎn)位移從局部坐標(biāo)系轉(zhuǎn)換為整體坐標(biāo)系,并采用有限元組裝和矩陣組合的方法,可得到每一個(gè)子結(jié)構(gòu)在整體坐標(biāo)系下的動力學(xué)剛度矩陣:

    (53)

    (54)

    式中n=1,2,3,4,5,6,7,8,9,10,11。進(jìn)而可以得到圖1(a)所示周期剛架結(jié)構(gòu)的節(jié)點(diǎn)位移與節(jié)點(diǎn)力的關(guān)系:

    F=S(ω)U

    (55)

    式中:S(ω)為3(7m+2)×3(7m+2)的動力學(xué)剛度矩陣;m是單胞周期數(shù)。對式(55)進(jìn)行求解,即可完成結(jié)構(gòu)動力學(xué)計(jì)算。

    2 數(shù)值模擬

    為分析預(yù)變形周期剛架結(jié)構(gòu)振動帶隙特性,取20個(gè)單胞組成整體結(jié)構(gòu),一端施加固定約束,每個(gè)單胞均由實(shí)心圓截面桿組成,其中橫桿和豎桿長度分別用l1和l2來表示,且l1=l2=1 m,圓截面半徑r=0.04 m,變形角α=π/12。材料參數(shù)為:ρ=7 800 kg/m3,E=210 GPa。通過求解式(10)特征值問題,可以得到結(jié)構(gòu)的固有頻率。

    為驗(yàn)證SEM建模和編程求解的正確性,采用有限元商業(yè)軟件ANSYS進(jìn)行仿真計(jì)算。此外,還計(jì)算了無預(yù)變形的X型周期剛架結(jié)構(gòu)(如圖1(b)虛線所示,即α=0)固有頻率。為保證合理性,基于質(zhì)量守恒定律,在確定2種結(jié)構(gòu)單胞跨度相等的前提下,調(diào)整X型結(jié)構(gòu)桿件的圓截面半徑,令其等于0.040 7 m。

    表1 2種不同剛架結(jié)構(gòu)固有頻率

    圖2 預(yù)變形周期剛架結(jié)構(gòu)前四階振型Fig.2 The first fourth vibration modes of the pre-deformed periodic rigid frame structure

    圖4表示的是周期剛架結(jié)構(gòu)的頻域分析,在固定端施加豎直向下、大小為1 N的簡諧力,通過求解式(55),取自由端點(diǎn)響應(yīng)與固定端位置響應(yīng)的比值,可以得到振動透射率??梢园l(fā)現(xiàn)該結(jié)構(gòu)的頻響曲線具有多處低谷,即振動帶隙,F(xiàn)EM結(jié)果較好地驗(yàn)證了SEM計(jì)算得到的振動帶隙結(jié)果的正確性。

    圖3 2種不同單胞結(jié)構(gòu)Fig.3 Two different unit cells

    圖4 2種不同周期剛架結(jié)構(gòu)的頻響透射率曲線Fig.4 Frequency response transmittance of two different rigid frame structures

    此外,通過比較可知,預(yù)變形剛架結(jié)構(gòu)第1個(gè)振動帶隙位置比X型剛架結(jié)構(gòu)出現(xiàn)在更低頻的位置,前者的頻率范圍大致是210~258 Hz,而后者的頻率范圍是270~338 Hz。并且預(yù)變形剛架結(jié)構(gòu)在1 576~1 753 Hz頻率范圍內(nèi)有一個(gè)更寬的振動帶隙。因此,可以得出結(jié)論,預(yù)變形周期結(jié)構(gòu)可以實(shí)現(xiàn)更加低頻的振動帶隙。

    圖4中通過透射率,即自由端點(diǎn)位置響應(yīng)與固定端位置響應(yīng)的比值來判斷結(jié)構(gòu)的振動帶隙范圍,當(dāng)自由端點(diǎn)的位置響應(yīng)較小時(shí),頻響曲線出現(xiàn)低谷,且響應(yīng)越小對應(yīng)的透射率也就越小,這就會出現(xiàn)結(jié)構(gòu)恰好不處于帶隙區(qū)間范圍內(nèi)而選擇的自由端點(diǎn)因受迫振動產(chǎn)生的響應(yīng)很小的情況,這種情況可以通過觀察帶隙區(qū)間外和區(qū)間內(nèi)對應(yīng)的模態(tài)進(jìn)行比較來辨別。比如從圖4中發(fā)現(xiàn)預(yù)變形結(jié)構(gòu)的第1個(gè)振動帶隙區(qū)間為210~258 Hz,所以選取了120 Hz和222 Hz 2個(gè)頻率,如圖5(a)、(b)分別為對應(yīng)頻率的模態(tài),可以看出120 Hz頻率下自由端點(diǎn)的位置響應(yīng)較大,且這種模態(tài)不屬于帶隙區(qū)間內(nèi),因?yàn)闆]有抑制振動波在結(jié)構(gòu)內(nèi)傳播;相對而言222 Hz頻率下自由端點(diǎn)的位置響應(yīng)很小,幾乎沒有,對應(yīng)的模態(tài)表明振動受到了抑制,后面幾個(gè)周期的單胞幾乎沒有受到影響,符合了結(jié)構(gòu)的帶隙特性,可以進(jìn)一步證明帶隙的存在以及結(jié)果正確性。

    圖5 預(yù)變形周期剛架結(jié)構(gòu)的振動模態(tài)Fig.5 Vibration modes of the pre-deformed rigid frame structures

    3 結(jié)論

    1)分析確定了結(jié)構(gòu)的振動帶隙,預(yù)變形結(jié)構(gòu)的剛度比X型結(jié)構(gòu)的剛度更低,證明了預(yù)變形結(jié)構(gòu)可以實(shí)現(xiàn)更低頻的振動帶隙的同時(shí),在高頻區(qū)域產(chǎn)生更寬頻的帶隙;

    2)在仿真計(jì)算中取自由端點(diǎn)響應(yīng)與固定端位置響應(yīng)的比值,得到振動透射率,通過FEM驗(yàn)證了上述理論結(jié)果,并定性分析了預(yù)變形結(jié)構(gòu)固有頻率更低的原因;

    3)通過自由端點(diǎn)位置響應(yīng)與固定端位置響應(yīng)的比值來判斷結(jié)構(gòu)的振動帶隙范圍,會出現(xiàn)結(jié)構(gòu)恰好不處于帶隙區(qū)間范圍內(nèi)而選擇的自由端點(diǎn)因受迫振動產(chǎn)生的響應(yīng)很小的情況,這種情況可以通過觀察模態(tài)進(jìn)行辨別,禁帶范圍內(nèi)振動波在經(jīng)過幾個(gè)周期內(nèi)就會得到充分抑制,可以進(jìn)一步證明帶隙的存在以及結(jié)果正確性;

    4)在此基礎(chǔ)上,可以在工程結(jié)構(gòu)設(shè)計(jì)過程中考慮調(diào)節(jié)變形角α的大小,來改變結(jié)構(gòu)的剛度,以獲得更符合條件的振動帶隙特性,但要兼顧結(jié)構(gòu)穩(wěn)定性的問題。

    猜你喜歡
    單胞架結(jié)構(gòu)端點(diǎn)
    海洋平臺板架結(jié)構(gòu)簡化評估對比研究
    某輕卡車架結(jié)構(gòu)強(qiáng)度性能優(yōu)化
    非特征端點(diǎn)條件下PM函數(shù)的迭代根
    基于NURBS的點(diǎn)陣材料參數(shù)化建模方法
    復(fù)合材料周期結(jié)構(gòu)數(shù)學(xué)均勻化方法的一種新型單胞邊界條件
    不等式求解過程中端點(diǎn)的確定
    叉車門架結(jié)構(gòu)分析和定期檢驗(yàn)探討
    考慮界面層影響的三維機(jī)織復(fù)合材料單胞模型研究
    參數(shù)型Marcinkiewicz積分算子及其交換子的加權(quán)端點(diǎn)估計(jì)
    底部鋼框架頂層門式剛架結(jié)構(gòu)建模研究
    岛国在线免费视频观看| 90打野战视频偷拍视频| 亚洲成a人片在线一区二区| 99国产综合亚洲精品| 午夜福利欧美成人| 国产爱豆传媒在线观看| 欧美乱色亚洲激情| 午夜视频国产福利| 国产午夜精品论理片| 看十八女毛片水多多多| 在线十欧美十亚洲十日本专区| 日日摸夜夜添夜夜添小说| 日韩欧美在线二视频| 免费观看的影片在线观看| 长腿黑丝高跟| 亚洲真实伦在线观看| 国产精品一区二区三区四区免费观看 | 青草久久国产| 色综合欧美亚洲国产小说| 脱女人内裤的视频| 午夜精品一区二区三区免费看| 久久精品人妻少妇| 身体一侧抽搐| 欧美一区二区精品小视频在线| 国产三级中文精品| 丁香六月欧美| 成人美女网站在线观看视频| 婷婷精品国产亚洲av| 99国产极品粉嫩在线观看| 亚洲激情在线av| 99国产综合亚洲精品| 欧美成狂野欧美在线观看| 久久久久性生活片| 免费一级毛片在线播放高清视频| 三级国产精品欧美在线观看| 97超视频在线观看视频| 亚洲欧美清纯卡通| 天堂影院成人在线观看| av天堂中文字幕网| 国产主播在线观看一区二区| 一进一出好大好爽视频| 久久久久久久精品吃奶| 亚洲av不卡在线观看| 国产伦精品一区二区三区四那| av福利片在线观看| 久久亚洲精品不卡| 91久久精品国产一区二区成人| 国产成人aa在线观看| 级片在线观看| 在线观看66精品国产| 亚洲欧美日韩无卡精品| 成人特级av手机在线观看| 久久久久国内视频| 亚洲国产精品久久男人天堂| 国产毛片a区久久久久| 欧美一区二区精品小视频在线| 亚洲激情在线av| 69人妻影院| 国产美女午夜福利| 天堂av国产一区二区熟女人妻| 在线a可以看的网站| 97人妻精品一区二区三区麻豆| 亚洲内射少妇av| 欧美xxxx黑人xx丫x性爽| 亚洲在线观看片| 天天一区二区日本电影三级| 亚洲成人精品中文字幕电影| 99国产综合亚洲精品| 久久久国产成人免费| 久久精品国产亚洲av天美| 久久久久免费精品人妻一区二区| 最好的美女福利视频网| 亚洲狠狠婷婷综合久久图片| 免费在线观看成人毛片| 久久久精品大字幕| 国产毛片a区久久久久| 亚洲精品一区av在线观看| 精品久久国产蜜桃| 男人的好看免费观看在线视频| 成熟少妇高潮喷水视频| 少妇高潮的动态图| 自拍偷自拍亚洲精品老妇| 国产欧美日韩精品亚洲av| 欧美成人a在线观看| 久久久久性生活片| 99精品在免费线老司机午夜| 久久香蕉精品热| 最新中文字幕久久久久| 亚洲国产精品成人综合色| 亚洲成a人片在线一区二区| 欧美潮喷喷水| 久久久久久九九精品二区国产| 高清日韩中文字幕在线| 日本 欧美在线| 精品午夜福利视频在线观看一区| 少妇熟女aⅴ在线视频| 欧美+亚洲+日韩+国产| 在线免费观看的www视频| 亚洲国产精品成人综合色| 又爽又黄无遮挡网站| 18+在线观看网站| 麻豆久久精品国产亚洲av| 免费av毛片视频| 日日干狠狠操夜夜爽| 欧美色欧美亚洲另类二区| 欧美一级a爱片免费观看看| 亚洲av二区三区四区| 五月伊人婷婷丁香| а√天堂www在线а√下载| 能在线免费观看的黄片| 级片在线观看| 国产精品一区二区性色av| 中文字幕免费在线视频6| 亚洲,欧美,日韩| av中文乱码字幕在线| 欧美日韩中文字幕国产精品一区二区三区| 18禁黄网站禁片午夜丰满| 久久久久久久亚洲中文字幕 | 欧美+亚洲+日韩+国产| 国产真实伦视频高清在线观看 | 舔av片在线| 夜夜夜夜夜久久久久| 亚洲在线自拍视频| 欧美bdsm另类| 性欧美人与动物交配| 久久亚洲精品不卡| 国产精品亚洲美女久久久| 男插女下体视频免费在线播放| 亚洲美女视频黄频| 女同久久另类99精品国产91| 最近最新免费中文字幕在线| 亚洲国产欧美人成| 久久国产乱子免费精品| 国产精品一区二区三区四区免费观看 | 亚洲最大成人av| av在线观看视频网站免费| 99国产精品一区二区三区| 国产精品嫩草影院av在线观看 | 日本黄色片子视频| 极品教师在线免费播放| 亚洲国产精品999在线| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲成a人片在线一区二区| 国产精品久久久久久人妻精品电影| 国产国拍精品亚洲av在线观看| 大型黄色视频在线免费观看| 亚洲经典国产精华液单 | 中文字幕人成人乱码亚洲影| 97碰自拍视频| 国产白丝娇喘喷水9色精品| 在线播放无遮挡| 亚洲av成人精品一区久久| 女人被狂操c到高潮| 久久久国产成人精品二区| 精品一区二区三区视频在线| 91在线精品国自产拍蜜月| 国产91精品成人一区二区三区| 深夜a级毛片| 人妻夜夜爽99麻豆av| 欧美最新免费一区二区三区 | 国产午夜精品论理片| 国产精品一区二区免费欧美| 精品无人区乱码1区二区| 久久久久九九精品影院| 亚洲av第一区精品v没综合| 亚洲av日韩精品久久久久久密| 他把我摸到了高潮在线观看| 在线天堂最新版资源| 亚洲三级黄色毛片| 桃色一区二区三区在线观看| 日韩欧美精品免费久久 | 亚洲天堂国产精品一区在线| 午夜精品久久久久久毛片777| 99国产精品一区二区蜜桃av| 我要看日韩黄色一级片| 我要看日韩黄色一级片| 亚洲av免费在线观看| 村上凉子中文字幕在线| 免费av不卡在线播放| 国产久久久一区二区三区| 真人一进一出gif抽搐免费| 淫秽高清视频在线观看| 男人舔女人下体高潮全视频| 人妻夜夜爽99麻豆av| 一进一出好大好爽视频| 国产乱人伦免费视频| 美女被艹到高潮喷水动态| 国产久久久一区二区三区| 一本一本综合久久| 在线观看午夜福利视频| 亚洲aⅴ乱码一区二区在线播放| 琪琪午夜伦伦电影理论片6080| 热99在线观看视频| 精品人妻视频免费看| 男人舔奶头视频| 成人国产一区最新在线观看| 午夜两性在线视频| 757午夜福利合集在线观看| 亚洲一区二区三区不卡视频| 看片在线看免费视频| 尤物成人国产欧美一区二区三区| 欧美国产日韩亚洲一区| 久久草成人影院| 国产精品伦人一区二区| 久久九九热精品免费| 狂野欧美白嫩少妇大欣赏| 亚洲精华国产精华精| 好男人在线观看高清免费视频| 国产精品美女特级片免费视频播放器| 日韩欧美精品免费久久 | 亚洲欧美精品综合久久99| 国产亚洲av嫩草精品影院| 精品久久国产蜜桃| 男女下面进入的视频免费午夜| 亚洲中文字幕一区二区三区有码在线看| 日本成人三级电影网站| 国产精品亚洲一级av第二区| 国产精品亚洲美女久久久| 夜夜躁狠狠躁天天躁| 日本三级黄在线观看| 久久精品国产99精品国产亚洲性色| 十八禁人妻一区二区| 禁无遮挡网站| 久久久久国内视频| 男女之事视频高清在线观看| 俄罗斯特黄特色一大片| 国产欧美日韩精品亚洲av| 欧美精品国产亚洲| 国产精品一及| 婷婷六月久久综合丁香| eeuss影院久久| 男人舔女人下体高潮全视频| 村上凉子中文字幕在线| 成人国产一区最新在线观看| 国产一级毛片七仙女欲春2| 少妇熟女aⅴ在线视频| 特级一级黄色大片| 婷婷精品国产亚洲av| 女同久久另类99精品国产91| 成人国产一区最新在线观看| 国产三级在线视频| 久久国产乱子伦精品免费另类| av国产免费在线观看| 无遮挡黄片免费观看| 最好的美女福利视频网| 美女大奶头视频| 国产91精品成人一区二区三区| 波多野结衣高清作品| 欧美性感艳星| 偷拍熟女少妇极品色| 午夜亚洲福利在线播放| 成人午夜高清在线视频| 成人无遮挡网站| 在现免费观看毛片| 欧美xxxx性猛交bbbb| 成人鲁丝片一二三区免费| 国产精品伦人一区二区| 亚洲精品乱码久久久v下载方式| 亚洲午夜理论影院| 免费av不卡在线播放| 成人国产一区最新在线观看| 国内揄拍国产精品人妻在线| 国产国拍精品亚洲av在线观看| 少妇裸体淫交视频免费看高清| 久久天躁狠狠躁夜夜2o2o| 国产精品美女特级片免费视频播放器| 亚洲欧美精品综合久久99| 欧美日韩亚洲国产一区二区在线观看| 午夜久久久久精精品| 日韩欧美在线乱码| 国内精品一区二区在线观看| 不卡一级毛片| 亚洲美女黄片视频| 精品一区二区免费观看| 久久久久久大精品| 69人妻影院| 亚洲五月婷婷丁香| 亚洲第一欧美日韩一区二区三区| 色5月婷婷丁香| 精品一区二区三区人妻视频| 91九色精品人成在线观看| 欧美区成人在线视频| 麻豆国产97在线/欧美| 国内少妇人妻偷人精品xxx网站| 欧美性猛交黑人性爽| 国产视频内射| 国产乱人伦免费视频| 男女床上黄色一级片免费看| 成人av在线播放网站| 在线国产一区二区在线| 午夜福利欧美成人| 亚洲最大成人手机在线| 国产精品野战在线观看| 99国产极品粉嫩在线观看| 偷拍熟女少妇极品色| 日韩成人在线观看一区二区三区| 日韩人妻高清精品专区| 又粗又爽又猛毛片免费看| 国产蜜桃级精品一区二区三区| 久久精品综合一区二区三区| 亚洲黑人精品在线| 热99re8久久精品国产| 亚洲欧美激情综合另类| 最近最新中文字幕大全电影3| 国产激情偷乱视频一区二区| 色5月婷婷丁香| 九九在线视频观看精品| 精品人妻1区二区| 天堂动漫精品| 男女之事视频高清在线观看| 精品一区二区三区视频在线| 亚洲欧美精品综合久久99| 丁香六月欧美| 小说图片视频综合网站| 国产伦精品一区二区三区视频9| 久久6这里有精品| 91麻豆av在线| 看片在线看免费视频| 成人欧美大片| 波野结衣二区三区在线| 免费看a级黄色片| 变态另类丝袜制服| 偷拍熟女少妇极品色| 在现免费观看毛片| 亚洲第一电影网av| 国产亚洲精品综合一区在线观看| 日本黄色视频三级网站网址| 精华霜和精华液先用哪个| 国产av不卡久久| 国产黄色小视频在线观看| 能在线免费观看的黄片| 日本一本二区三区精品| 看十八女毛片水多多多| 日韩欧美精品v在线| 国产精品精品国产色婷婷| 最近在线观看免费完整版| 国产在视频线在精品| 国产精品乱码一区二三区的特点| 亚洲美女黄片视频| 午夜福利免费观看在线| 国产欧美日韩一区二区精品| 天堂影院成人在线观看| 国产伦一二天堂av在线观看| 亚洲专区中文字幕在线| 97超视频在线观看视频| 男女视频在线观看网站免费| 国产精品伦人一区二区| 小蜜桃在线观看免费完整版高清| 国内精品美女久久久久久| 国产欧美日韩精品亚洲av| 欧美成人a在线观看| 欧美中文日本在线观看视频| 脱女人内裤的视频| 白带黄色成豆腐渣| 亚洲国产精品999在线| 午夜福利免费观看在线| 一本久久中文字幕| 日本三级黄在线观看| 亚洲国产日韩欧美精品在线观看| 三级国产精品欧美在线观看| 国产高清激情床上av| 搡老妇女老女人老熟妇| 性欧美人与动物交配| 亚洲午夜理论影院| 亚洲电影在线观看av| www日本黄色视频网| 嫩草影院新地址| 少妇人妻精品综合一区二区 | 亚洲av二区三区四区| 国内少妇人妻偷人精品xxx网站| 国产精品av视频在线免费观看| 波多野结衣高清作品| 国内毛片毛片毛片毛片毛片| 国产麻豆成人av免费视频| 国产欧美日韩精品一区二区| 久久久久久九九精品二区国产| 少妇被粗大猛烈的视频| 国产精品久久久久久精品电影| 国产精品av视频在线免费观看| 91午夜精品亚洲一区二区三区 | 三级国产精品欧美在线观看| 一级av片app| 欧美一级a爱片免费观看看| 波多野结衣巨乳人妻| 亚洲 欧美 日韩 在线 免费| www日本黄色视频网| 亚洲av第一区精品v没综合| 国产中年淑女户外野战色| 韩国av一区二区三区四区| 国产男靠女视频免费网站| 久久精品影院6| 毛片女人毛片| 成人精品一区二区免费| 亚洲成人免费电影在线观看| 丰满的人妻完整版| 亚洲精品久久国产高清桃花| 国产在视频线在精品| 国产大屁股一区二区在线视频| 美女被艹到高潮喷水动态| 全区人妻精品视频| 美女黄网站色视频| 1024手机看黄色片| 欧美一级a爱片免费观看看| 一个人看的www免费观看视频| 色综合亚洲欧美另类图片| 久久亚洲真实| 国产精品永久免费网站| 久久久久国产精品人妻aⅴ院| 18禁黄网站禁片免费观看直播| 男女之事视频高清在线观看| 久久精品国产99精品国产亚洲性色| 精品福利观看| 99久久精品热视频| 久久久久国内视频| 精品一区二区三区av网在线观看| 久久香蕉精品热| 欧美精品国产亚洲| 亚洲人成电影免费在线| 国产免费男女视频| 久久精品国产亚洲av香蕉五月| 国产高清有码在线观看视频| 老司机午夜福利在线观看视频| 波野结衣二区三区在线| 久久中文看片网| aaaaa片日本免费| 国产精品国产高清国产av| 免费av毛片视频| 久久午夜亚洲精品久久| 国产探花在线观看一区二区| 日韩欧美国产一区二区入口| 中文字幕av在线有码专区| 99国产精品一区二区三区| 97人妻精品一区二区三区麻豆| 国产精品久久视频播放| 日韩欧美国产在线观看| 日本免费一区二区三区高清不卡| 99国产综合亚洲精品| 亚洲av不卡在线观看| 国产老妇女一区| www.999成人在线观看| 日本黄色片子视频| 波野结衣二区三区在线| 亚洲国产欧美人成| 精品午夜福利视频在线观看一区| 又紧又爽又黄一区二区| 国产亚洲欧美在线一区二区| 成人鲁丝片一二三区免费| 波野结衣二区三区在线| 琪琪午夜伦伦电影理论片6080| 欧美日韩亚洲国产一区二区在线观看| 三级国产精品欧美在线观看| 制服丝袜大香蕉在线| 国内久久婷婷六月综合欲色啪| 国产免费一级a男人的天堂| 日日干狠狠操夜夜爽| 免费在线观看日本一区| 成人三级黄色视频| 欧美xxxx黑人xx丫x性爽| 国产成人a区在线观看| 午夜福利18| 亚洲精品久久国产高清桃花| 99久久99久久久精品蜜桃| 宅男免费午夜| 国产日本99.免费观看| 精品福利观看| 亚洲av免费在线观看| 日本撒尿小便嘘嘘汇集6| 久久国产精品影院| 九九热线精品视视频播放| 欧美丝袜亚洲另类 | 高潮久久久久久久久久久不卡| 亚洲精品亚洲一区二区| 欧美xxxx黑人xx丫x性爽| 色综合婷婷激情| 国产精品,欧美在线| 首页视频小说图片口味搜索| 成人av一区二区三区在线看| 一个人观看的视频www高清免费观看| av在线老鸭窝| 国产三级在线视频| 成年人黄色毛片网站| 色5月婷婷丁香| 国产成年人精品一区二区| 中出人妻视频一区二区| 一区二区三区高清视频在线| 亚洲av熟女| 国产精品野战在线观看| 97人妻精品一区二区三区麻豆| 男人和女人高潮做爰伦理| 又黄又爽又免费观看的视频| 床上黄色一级片| 丰满人妻熟妇乱又伦精品不卡| 亚洲精品一区av在线观看| 日韩亚洲欧美综合| 国产精品1区2区在线观看.| 国产一区二区三区在线臀色熟女| 国产 一区 欧美 日韩| 欧美日本亚洲视频在线播放| 国产成人影院久久av| 成人高潮视频无遮挡免费网站| 欧美一区二区国产精品久久精品| 国产精品久久电影中文字幕| 日本免费一区二区三区高清不卡| 精品久久久久久久末码| 在线播放无遮挡| 一二三四社区在线视频社区8| 色av中文字幕| 成人国产一区最新在线观看| 久久久久九九精品影院| 国产毛片a区久久久久| 我要看日韩黄色一级片| 精品午夜福利视频在线观看一区| 熟女电影av网| 欧美日韩乱码在线| 麻豆一二三区av精品| 欧美黑人欧美精品刺激| 日本在线视频免费播放| 久久人人精品亚洲av| 好男人在线观看高清免费视频| 变态另类成人亚洲欧美熟女| 国内毛片毛片毛片毛片毛片| 99热精品在线国产| 在线免费观看的www视频| 亚洲av电影不卡..在线观看| 午夜福利在线在线| 欧美色视频一区免费| 久久国产精品影院| 日本a在线网址| 小蜜桃在线观看免费完整版高清| 国产91精品成人一区二区三区| 婷婷丁香在线五月| 亚洲美女黄片视频| 亚洲人成伊人成综合网2020| 亚洲第一区二区三区不卡| 少妇人妻一区二区三区视频| 国产亚洲av嫩草精品影院| 国产在线男女| 内地一区二区视频在线| 亚洲专区中文字幕在线| 禁无遮挡网站| 国语自产精品视频在线第100页| 俄罗斯特黄特色一大片| 免费av观看视频| 88av欧美| 高清在线国产一区| 亚洲自拍偷在线| 老鸭窝网址在线观看| 精品久久久久久久人妻蜜臀av| eeuss影院久久| 午夜福利视频1000在线观看| 老司机深夜福利视频在线观看| 亚洲无线在线观看| 欧美最黄视频在线播放免费| 97超视频在线观看视频| 直男gayav资源| av欧美777| 一级av片app| 国产野战对白在线观看| 国产中年淑女户外野战色| 国产精品综合久久久久久久免费| 亚洲成av人片免费观看| 成人特级黄色片久久久久久久| 一区二区三区免费毛片| 熟妇人妻久久中文字幕3abv| 我要看日韩黄色一级片| 18+在线观看网站| 久久久久久大精品| 赤兔流量卡办理| 久久久久久大精品| 怎么达到女性高潮| 一二三四社区在线视频社区8| 国产成人a区在线观看| 国产私拍福利视频在线观看| 精品久久久久久成人av| 男人舔女人下体高潮全视频| 男人狂女人下面高潮的视频| 国产野战对白在线观看| 少妇人妻精品综合一区二区 | 国产成人啪精品午夜网站| 97人妻精品一区二区三区麻豆| 97超级碰碰碰精品色视频在线观看| 日韩欧美 国产精品| 日本一二三区视频观看| 在线观看免费视频日本深夜| 在线观看午夜福利视频| 婷婷精品国产亚洲av| 日韩av在线大香蕉| 婷婷亚洲欧美| 欧美成人性av电影在线观看| 日韩国内少妇激情av| 嫩草影院新地址| 搞女人的毛片| 91在线精品国自产拍蜜月| 黄色一级大片看看| 欧美黑人巨大hd| 欧美日韩综合久久久久久 | 热99在线观看视频| 直男gayav资源| 欧美午夜高清在线| 亚洲天堂国产精品一区在线| 精品乱码久久久久久99久播| 国产aⅴ精品一区二区三区波| 97超视频在线观看视频| 99热这里只有是精品50| 久久精品影院6| 欧美潮喷喷水| 欧美日韩中文字幕国产精品一区二区三区| 午夜视频国产福利| 露出奶头的视频| 欧美+日韩+精品| 91在线精品国自产拍蜜月| 日韩精品中文字幕看吧| 欧美成人a在线观看| .国产精品久久| 高清毛片免费观看视频网站|