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

    壓電材料雙曲殼熱彈耦合作用下的混沌運(yùn)動?

    2019-01-25 09:53:26李林利薛春霞
    物理學(xué)報(bào) 2019年1期
    關(guān)鍵詞:雙曲壓電溫度場

    李林利 薛春霞

    (中北大學(xué)理學(xué)院,太原 030051)

    (2018年9月15日收到;2018年11月14日收到修改稿)

    運(yùn)用彈性力學(xué)有限變形基本理論推導(dǎo)出了壓電材料雙曲殼在外激力和溫度場作用下的非線性振動方程和協(xié)調(diào)方程.通過Bubnov-Galerkin原理,得到該結(jié)構(gòu)的非線性動力學(xué)方程.利用Melnikov方法,得到系統(tǒng)產(chǎn)生Smale馬蹄變換意義下混沌的條件,用四階Runge-Kutta法編寫程序?qū)ο到y(tǒng)進(jìn)行數(shù)值求解,并繪制出相應(yīng)的分岔圖、Lyapunov指數(shù)圖、相軌跡圖以及Poincaré截面圖,分析了溫度場對壓電材料雙曲殼系統(tǒng)的非線性特性的影響.仿真結(jié)果表明,隨著溫度的升高,系統(tǒng)的混沌與周期區(qū)交替出現(xiàn),溫度場的改變可影響和控制系統(tǒng)的振動特性.

    1 引 言

    隨著現(xiàn)代高新科學(xué)技術(shù)的迅猛發(fā)展,壓電材料作為智能材料的一種,具有薄而輕、成本低、響應(yīng)速度快、轉(zhuǎn)換效率高、壽命長等優(yōu)點(diǎn),被廣泛應(yīng)用于眾多高科技設(shè)備和壓電智能器件中.而雙曲殼結(jié)構(gòu)作為工程中常見的結(jié)構(gòu)之一,常在機(jī)械場和溫度場等多場作用的環(huán)境下工作,其動力特性對系統(tǒng)的結(jié)構(gòu)安全有重要影響.因此對壓電材料雙曲殼在熱彈耦合作用下的動力學(xué)分析具有重要的理論意義和實(shí)用價(jià)值.國內(nèi)外學(xué)者針對普通彈性梁板結(jié)構(gòu)在熱、電、磁多物理場耦合作用下混沌非線性動力學(xué)行為的研究,已取得了一定的成果[1-5].薛春霞等[6]分析了磁感應(yīng)強(qiáng)度和外載荷對金屬薄板振動特性的影響.陳趙江等[7]通過實(shí)驗(yàn)方法和理論推導(dǎo)對高功率超聲脈沖激勵下金屬薄板的振動特性進(jìn)行了分析研究.與此同時(shí),對一些特殊功能材料板振動問題的研究也取得了顯著成就[8,9].Hao等[10]首次在熱環(huán)境下分析了簡支功能梯度材料(FGMS)矩形板在橫向和平面激勵下的非線性動力學(xué)問題.栗蕾[11]系統(tǒng)地研究了考慮初始缺陷的網(wǎng)殼結(jié)構(gòu)的非線性動力學(xué)行為和系統(tǒng)整體穩(wěn)定性等問題.周運(yùn)朱[12]通過幾何非線性和幾何材料雙重非線性,分析和考察了復(fù)雜曲面單層網(wǎng)殼結(jié)構(gòu)的穩(wěn)定性能.曹洲等[13]分析了各物理參數(shù)和溫度場對功能梯度夾層雙曲殼自由振動固有頻率的影響.此外,數(shù)值模擬驗(yàn)證非線性系統(tǒng)的混沌動力學(xué)特性方面的研究對理論的分析奠定了必要的基礎(chǔ)[14-16].

    綜上所述,對于雙曲殼結(jié)構(gòu)在熱彈耦合作用下的混沌運(yùn)動研究的較少.本文以壓電材料雙曲殼結(jié)構(gòu)為研究對象,通過MATLAB軟件建立數(shù)學(xué)模型并進(jìn)行數(shù)值計(jì)算,對在溫度場和簡諧外激勵共同作用下的非線性振動行為進(jìn)行分析.結(jié)果表明,溫度場的改變可影響和控制系統(tǒng)的振動特性.

    2 壓電材料雙曲殼的基本方程

    薄殼振動理論的基本假設(shè):

    1)變形前垂直中曲面的直線在變形后仍保持直線,并垂直中曲面;

    2)沿中曲面垂直方向的法向應(yīng)力忽略不計(jì);

    3)忽略不計(jì)轉(zhuǎn)動慣性力矩和面內(nèi)慣性力.

    四邊簡支壓電材料雙曲殼模型如圖1所示,在雙曲殼模型的中面建立正交曲線坐標(biāo)系(x1,x2,x3),坐標(biāo)原點(diǎn)為o點(diǎn).圖中x1,x2分別是中面曲率線方向,x3表示沿著厚度的方向,向上為正方向;a表示沿x1坐標(biāo)方向的中面曲線長度,b表示沿x2坐標(biāo)方向的中面曲線長度;R1,R2分別表示x1,x2坐標(biāo)方向上的曲率半徑;h表示雙曲殼的厚度.

    圖1 壓電材料雙曲殼模型Fig.1.Hyperbolic shell model for piezoelectric materials.

    對于橫觀各向同性壓電材料的雙曲殼,在考慮溫度效應(yīng)時(shí)的壓電方程[17]如下:

    其中t31=Eα/(1-2υ)為熱-機(jī)械耦合常數(shù),E和υ分別為楊氏模量和泊松比,α為線性膨脹系數(shù);C為比熱,c11,c12,c21,c22,c66表示壓電材料的剛度;e31為壓電系數(shù),ε3為介電系數(shù),m3為熱-壓電耦合常數(shù),E3為電場,T為溫度增量;而σ1,σ2和τ12分別表示x1方向和x2方向的正應(yīng)力和切應(yīng)力;DE表示電位移,S表示熵,α0=C/T0,T0是初始溫度.

    應(yīng)用已知的扁殼大撓度幾何方程,得到雙曲殼的位移應(yīng)變關(guān)系為[18]:

    其中,k1=1/R1,k2=1/R2.

    設(shè)變溫T(x1,x2,x3)=T0(x1,x2)+x3(x1,x2),其中,T0(x1,x2)是初始狀態(tài)時(shí),沿殼體厚度方向分布的平均溫度,(x1,x2)是溫度梯度.設(shè)f1(ξ,η)和f2(ξ,η)是殼體的內(nèi)、外表面的溫度,而且

    其中,

    薄膜內(nèi)力

    彎曲內(nèi)力

    雙曲殼的運(yùn)動方程[19]

    其中,Q1,Q2分別是垂直于x1軸及垂直于x2軸的截面上單位寬度的橫向剪力.通過引入Airy應(yīng)力函數(shù)φ可以簡化求解過程:

    將(7)式代入(5)式,再代入方程(6),由于c11=c22,注意到c11=c12+2c66,可得到壓電材料矩形板的橫向熱振動方程如下:

    協(xié)調(diào)方程為:

    ?Γ為高斯曲率的改變,

    由于雙曲殼是扁殼,可近似地由撓度w所產(chǎn)生的中曲面的主曲率和扭率表示

    結(jié)合(1),(2)和(7)式,并代入方程(9)中,協(xié)調(diào)方程變型為

    方程(8)和(11)是熱彈耦合作用下的薄殼大撓度問題的非線性振動方程組.不考慮材料壓電特性時(shí),對比文獻(xiàn)[20,21],對于熱彈耦合作用下的普通彈性板殼結(jié)構(gòu)該方程組仍適用,故本文的方法和結(jié)論具有可靠性.

    四邊簡支的雙曲殼的邊界條件為:

    為簡化計(jì)算,令R1=R2=R,引入無量綱參數(shù)ξ=x1/R,η=x2/R,橫向位移w設(shè)為

    式中α1= πR/a,β1= πR/b.

    由方程(8)得:

    將方程(12)代入方程(11),解出應(yīng)力函數(shù)φ再代入方程(8),利用Bubnov-Galerkin原理,由

    其中

    求得

    其中,

    引入:w0(t)=f(t)+并代入方程(16),得到:

    3 系統(tǒng)非線性分岔與混沌分析

    方程(17)可改寫為我們熟悉的非線性動力系統(tǒng)的形式:

    (18)式是熱彈耦合作用下的非線性動力學(xué)方程.其中,

    λ1,λ2均大于0,便是我們要求的熱彈耦合作用下雙曲殼的大撓度問題的非線性動力系統(tǒng).

    雙曲殼在強(qiáng)迫振動作用下的非線性動力系統(tǒng)的無擾動力系統(tǒng)為

    系統(tǒng)(19)的Hamilton量為

    H值由運(yùn)動的初始條件確定.從物理意義上看,方程(19)可看成是描述單位質(zhì)量的質(zhì)點(diǎn)在外力-λ1x+λ2x2-λ3x3作用下所做的運(yùn)動,而Hamilton量中表示質(zhì)點(diǎn)的動能,表示外力對質(zhì)點(diǎn)所做功的負(fù)值,即為勢能,H代表系統(tǒng)總的機(jī)械能.不同的H值,(20)式代表相平面上不同的曲線.

    無擾動系統(tǒng)的Jocabi矩陣為

    特征方程為

    圖2 同宿軌道Fig.2.Homoclinic orbit.

    可解出同宿軌道的參數(shù)方程為

    對動力系統(tǒng)定義Melnikov函數(shù)

    由Melnikov理論可知,若Melnikov函數(shù)存在零點(diǎn),則穩(wěn)定流形和不穩(wěn)定流形相交,存在不變集合,系統(tǒng)有可能會出現(xiàn)混沌.(24)式中等式右端第一項(xiàng)是阻尼力δ˙x沿軌道所做的功,第二項(xiàng)是以τ所處位置為出發(fā)點(diǎn)外力沿軌道所做的功.

    使|sinωτ|<1,得到

    且僅當(dāng)(26)式成立時(shí),動力系統(tǒng)(18)具有Smale馬蹄變換意義下的混沌運(yùn)動的閾值條件.如圖3所示,曲線上方表示可能出現(xiàn)混沌的參數(shù)范圍.

    圖3 混沌閾值曲線Fig.3.Chaotic threshold curve.

    4 數(shù)值計(jì)算

    本文通過MATLAB中simulink仿真模塊對系統(tǒng)建模,具體的仿真模型如圖4所示.取雙曲殼結(jié)構(gòu)的材料為BaTiO3陶瓷,材料參數(shù)列于表1.雙曲殼長a=400 mm,寬b=400 mm,厚h=70 mm,曲率半徑R1=R2=500 mm.將這些參數(shù)代入(26)式計(jì)算,通過改變溫度T0的大小,使方程(26)成立時(shí),系統(tǒng)就可能出現(xiàn)Smale馬蹄變換意義下的混沌.

    表1 鈦酸鋇材料參數(shù)[22]Table 1.Barium titanate material parameters[22].

    利用四階Runge-Kutta法編制數(shù)值仿真程序,通過改變溫度T0得到系統(tǒng)的分岔圖和相對應(yīng)的Lyapunov指數(shù)圖,如圖5和圖6所示.可以看出隨著溫度的增加,系統(tǒng)的混沌與周期區(qū)交替出現(xiàn).由圖6可看出,溫度T0在32?C和41?C附近以及36?C—37?C的范圍內(nèi),Lyapunov指數(shù)小于0,對應(yīng)圖5中的周期區(qū);而當(dāng)圖6中的Lyapunov指數(shù)大于0時(shí),對應(yīng)圖5中的混沌區(qū).

    圖4 仿真模型Fig.4.Simulation model.

    圖5 隨溫度T0變化的分岔圖Fig.5.Bifurcation diagram with temperature T0change.

    圖6 隨溫度T0變化的Lyapunov指數(shù)圖Fig.6.Lyapunov exponent diagram varying with temperature T0.

    當(dāng)ω0=17 Hz,q0=200 N/m2時(shí),通過改變溫度T0的大小,分別做出熱彈耦合作用下雙曲殼系統(tǒng)的相圖和Poincaré截面圖.

    由上面的分析可知,仿真結(jié)果和Melnikov方法的判定一致.此外,仿真結(jié)果反映出的混沌運(yùn)動更加復(fù)雜,隨著溫度的增加,系統(tǒng)的混沌與周期區(qū)交替出現(xiàn),而Melnikov方法是無法得到這樣的結(jié)論的.

    由圖7—圖9可以看出,圖7相圖和時(shí)程曲線很規(guī)整,相圖繞了2圈,Poincaré映射為兩個點(diǎn),故代表周期為2的運(yùn)動;同理可知,圖8代表周期為三的運(yùn)動.由圖7和圖8可看出,k周期運(yùn)動有明顯的周期重復(fù)性,對應(yīng)的相圖繞k圈,Poincaré截面圖顯示k個點(diǎn).而圖9的相圖軌跡互不重疊,相互纏繞,Poincaré截面圖也雜亂無章,結(jié)合這組參數(shù)已落在Melnikov函數(shù)和系統(tǒng)分岔圖給定的混沌區(qū),且對應(yīng)的Lyapunov指數(shù)大于0,因此可判斷圖9代表著混沌運(yùn)動.

    從計(jì)算結(jié)果可以看出,壓電材料雙曲殼在四邊簡支的條件下,隨著溫度場的變化,會改變系統(tǒng)的運(yùn)動狀態(tài).就物理意義而言,溫度的變化給系統(tǒng)的剛度造成了附加效應(yīng).由(16)式可看出,溫度的變化使l1發(fā)生變化,相當(dāng)于系統(tǒng)的剛度發(fā)生變化,從而影響和改變了系統(tǒng)的振動特性.

    圖7 T0=22?C時(shí)的相圖和Poincaré截面圖Fig.7.Phase diagrams and Poincaré section diagrams at T0=22 ?C.

    圖8 T0=36?C時(shí)的相圖和Poincaré截面圖Fig.8.Phase diagrams and Poincaré section diagrams at T0=36 ?C.

    圖9 T0=45?C時(shí)的相圖和Poincaré截面圖Fig.9.Phase diagrams and Poincaré section diagrams at T0=45 ?C.

    5 結(jié) 論

    本文研究了熱彈耦合效應(yīng)下簡支的壓電材料雙曲殼的動力學(xué)行為.利用Melnikov理論,推導(dǎo)出系統(tǒng)出現(xiàn)混沌運(yùn)動的判據(jù),通過MATLAB軟件進(jìn)行數(shù)值模擬,在熱彈耦合作用下系統(tǒng)的振動方程具有明顯的非線性,隨著溫度的增加,運(yùn)動行為較復(fù)雜,系統(tǒng)的混沌與周期區(qū)交替出現(xiàn).因此,控制和調(diào)整溫度場的變化,可實(shí)現(xiàn)控制系統(tǒng)的運(yùn)動狀態(tài),為系統(tǒng)安全性和可靠性提供參考.

    猜你喜歡
    雙曲壓電溫度場
    中國科學(xué)技術(shù)館之“雙曲隧道”
    軍事文摘(2021年22期)2022-01-18 06:22:48
    鋁合金加筋板焊接溫度場和殘余應(yīng)力數(shù)值模擬
    雙曲型交換四元數(shù)的極表示
    《壓電與聲光》征稿啟事
    壓電與聲光(2019年1期)2019-02-22 09:46:06
    基于紋影法的溫度場分布測量方法
    MJS工法與凍結(jié)法結(jié)合加固區(qū)溫度場研究
    建筑科技(2018年6期)2018-08-30 03:41:08
    新型壓電疊堆泵設(shè)計(jì)及仿真
    一階雙曲型偏微分方程的模糊邊界控制
    基于雙曲和代數(shù)多項(xiàng)式的HC-Bézier曲線
    X80鋼層流冷卻溫度場的有限元模擬
    日韩一卡2卡3卡4卡2021年| 成人av一区二区三区在线看| 两个人免费观看高清视频| 亚洲无线在线观看| 美女 人体艺术 gogo| 国产成人av激情在线播放| 少妇熟女aⅴ在线视频| 亚洲伊人色综图| АⅤ资源中文在线天堂| 麻豆国产av国片精品| 亚洲欧洲精品一区二区精品久久久| 美国免费a级毛片| 国产一级毛片七仙女欲春2 | 久久香蕉国产精品| 久久午夜亚洲精品久久| 国产99久久九九免费精品| 黄网站色视频无遮挡免费观看| 久久国产精品男人的天堂亚洲| 午夜福利免费观看在线| 99久久精品国产亚洲精品| 一级毛片精品| 天天躁夜夜躁狠狠躁躁| 久久精品人人爽人人爽视色| 老司机午夜十八禁免费视频| 精品国产国语对白av| 日本免费一区二区三区高清不卡 | 热re99久久国产66热| 亚洲一区二区三区色噜噜| 国产欧美日韩一区二区三| 亚洲一码二码三码区别大吗| 免费久久久久久久精品成人欧美视频| 97人妻天天添夜夜摸| 中文字幕人成人乱码亚洲影| av在线天堂中文字幕| 国产麻豆成人av免费视频| 国内精品久久久久精免费| 日本黄色视频三级网站网址| 神马国产精品三级电影在线观看 | 亚洲成人久久性| 国产99白浆流出| 变态另类成人亚洲欧美熟女 | av网站免费在线观看视频| 男人的好看免费观看在线视频 | 久久狼人影院| 欧美老熟妇乱子伦牲交| 亚洲精品美女久久av网站| 免费看美女性在线毛片视频| 欧美成人性av电影在线观看| 国产成人影院久久av| 在线观看午夜福利视频| 久久久国产成人精品二区| 91九色精品人成在线观看| 日韩欧美国产一区二区入口| 国产三级黄色录像| 日韩一卡2卡3卡4卡2021年| 国产成人精品久久二区二区免费| 村上凉子中文字幕在线| 亚洲国产精品成人综合色| 亚洲第一欧美日韩一区二区三区| 色婷婷久久久亚洲欧美| 久久天堂一区二区三区四区| 午夜精品在线福利| 丝袜美足系列| 国产在线观看jvid| 大型av网站在线播放| 成人三级做爰电影| 精品不卡国产一区二区三区| 久久久水蜜桃国产精品网| 91av网站免费观看| 午夜久久久在线观看| 女性被躁到高潮视频| 久久国产精品男人的天堂亚洲| 免费女性裸体啪啪无遮挡网站| 性欧美人与动物交配| 精品国产国语对白av| 最好的美女福利视频网| 免费看美女性在线毛片视频| 国产精品 欧美亚洲| 亚洲专区国产一区二区| АⅤ资源中文在线天堂| 精品无人区乱码1区二区| 精品卡一卡二卡四卡免费| 成人国产一区最新在线观看| 久久人妻熟女aⅴ| 成人三级做爰电影| 国产午夜福利久久久久久| 日本精品一区二区三区蜜桃| 精品国产一区二区久久| 99国产精品一区二区三区| 国产一区在线观看成人免费| 国产精品久久久久久人妻精品电影| 亚洲国产精品久久男人天堂| 亚洲成人免费电影在线观看| 亚洲狠狠婷婷综合久久图片| 亚洲精品国产一区二区精华液| 成年版毛片免费区| 人人妻人人爽人人添夜夜欢视频| 免费观看人在逋| 国产av精品麻豆| 久久精品91蜜桃| 午夜视频精品福利| 亚洲中文字幕日韩| 国产真人三级小视频在线观看| 日本精品一区二区三区蜜桃| 日韩欧美免费精品| 久久精品国产亚洲av高清一级| 又大又爽又粗| 国产国语露脸激情在线看| 日本在线视频免费播放| 看免费av毛片| 成年版毛片免费区| 日本五十路高清| 国产欧美日韩一区二区三| 久久中文字幕人妻熟女| 亚洲人成电影观看| 久久久久久亚洲精品国产蜜桃av| 黄片播放在线免费| 99国产精品一区二区三区| 一进一出抽搐gif免费好疼| 色综合亚洲欧美另类图片| 亚洲人成电影观看| 99久久99久久久精品蜜桃| 国产三级在线视频| 嫁个100分男人电影在线观看| 50天的宝宝边吃奶边哭怎么回事| 日本五十路高清| 男女下面插进去视频免费观看| 亚洲国产精品久久男人天堂| 久久人人97超碰香蕉20202| 叶爱在线成人免费视频播放| 国产高清videossex| 亚洲国产高清在线一区二区三 | 国产私拍福利视频在线观看| 无人区码免费观看不卡| 在线观看66精品国产| 99热只有精品国产| 免费不卡黄色视频| 亚洲欧美日韩另类电影网站| 在线天堂中文资源库| 国产精品秋霞免费鲁丝片| 热re99久久国产66热| 又黄又爽又免费观看的视频| 热99re8久久精品国产| 色尼玛亚洲综合影院| 91精品国产国语对白视频| 国产精品亚洲美女久久久| 亚洲欧美精品综合一区二区三区| 亚洲国产毛片av蜜桃av| 国产精品爽爽va在线观看网站 | 久久婷婷成人综合色麻豆| 夜夜看夜夜爽夜夜摸| 黑人巨大精品欧美一区二区蜜桃| www.999成人在线观看| 精品免费久久久久久久清纯| 99久久精品国产亚洲精品| 亚洲色图综合在线观看| 激情在线观看视频在线高清| 好男人电影高清在线观看| 最好的美女福利视频网| 久久亚洲精品不卡| 亚洲av日韩精品久久久久久密| 丝袜人妻中文字幕| 精品欧美一区二区三区在线| 色播在线永久视频| 丝袜在线中文字幕| 岛国视频午夜一区免费看| 亚洲国产欧美一区二区综合| 在线av久久热| 国产亚洲精品av在线| 好男人电影高清在线观看| 老汉色av国产亚洲站长工具| 久久精品成人免费网站| 多毛熟女@视频| 精品久久久精品久久久| 大型黄色视频在线免费观看| 欧美成人一区二区免费高清观看 | 欧美精品亚洲一区二区| 亚洲精品中文字幕在线视频| 亚洲国产精品999在线| 天天添夜夜摸| 真人一进一出gif抽搐免费| 欧美一级毛片孕妇| 国产片内射在线| 国产精品99久久99久久久不卡| 日韩免费av在线播放| 国产私拍福利视频在线观看| 日本免费a在线| 欧美成狂野欧美在线观看| 99精品久久久久人妻精品| 成人国语在线视频| 正在播放国产对白刺激| 亚洲欧美精品综合一区二区三区| 精品国产一区二区三区四区第35| 黄色视频,在线免费观看| 中文字幕人妻熟女乱码| 久久人人爽av亚洲精品天堂| 一边摸一边做爽爽视频免费| 久热这里只有精品99| 国产成人精品无人区| 大码成人一级视频| 亚洲欧美精品综合久久99| xxx96com| 亚洲av成人一区二区三| 在线观看一区二区三区| 97超级碰碰碰精品色视频在线观看| 变态另类成人亚洲欧美熟女 | 日日摸夜夜添夜夜添小说| 国产又爽黄色视频| 午夜老司机福利片| 亚洲精品在线美女| av视频在线观看入口| 中文字幕久久专区| 中文字幕人成人乱码亚洲影| 极品人妻少妇av视频| 又紧又爽又黄一区二区| 欧美日韩亚洲综合一区二区三区_| 国产精品久久电影中文字幕| 欧美一区二区精品小视频在线| 免费女性裸体啪啪无遮挡网站| 91大片在线观看| 97碰自拍视频| 午夜福利视频1000在线观看 | 一级片免费观看大全| 亚洲色图综合在线观看| 在线永久观看黄色视频| 在线观看日韩欧美| 亚洲第一欧美日韩一区二区三区| 一夜夜www| 满18在线观看网站| 亚洲欧美精品综合久久99| 日韩国内少妇激情av| 午夜福利高清视频| 老汉色av国产亚洲站长工具| 久久久精品国产亚洲av高清涩受| 精品卡一卡二卡四卡免费| 国产视频一区二区在线看| 99re在线观看精品视频| 日本a在线网址| 成年女人毛片免费观看观看9| 一本大道久久a久久精品| 亚洲国产毛片av蜜桃av| 中文字幕另类日韩欧美亚洲嫩草| 三级毛片av免费| 亚洲中文日韩欧美视频| 黑人巨大精品欧美一区二区蜜桃| 美女高潮到喷水免费观看| 精品熟女少妇八av免费久了| 精品第一国产精品| 亚洲精品在线美女| АⅤ资源中文在线天堂| 一级作爱视频免费观看| 99国产精品99久久久久| 999久久久国产精品视频| 满18在线观看网站| 精品国产乱子伦一区二区三区| 欧美在线黄色| 老司机午夜十八禁免费视频| 久久婷婷成人综合色麻豆| av视频在线观看入口| 久久青草综合色| 国产一区在线观看成人免费| 国产国语露脸激情在线看| 精品少妇一区二区三区视频日本电影| 欧美激情高清一区二区三区| 波多野结衣巨乳人妻| 欧美午夜高清在线| √禁漫天堂资源中文www| 久久久久亚洲av毛片大全| 嫁个100分男人电影在线观看| 免费在线观看视频国产中文字幕亚洲| 成人免费观看视频高清| 亚洲 欧美一区二区三区| 午夜福利一区二区在线看| 免费不卡黄色视频| 成人国产一区最新在线观看| 久久午夜亚洲精品久久| 日韩中文字幕欧美一区二区| 69av精品久久久久久| 久久国产亚洲av麻豆专区| 久久精品国产清高在天天线| 国产激情久久老熟女| 亚洲精品一区av在线观看| 国产精华一区二区三区| 午夜激情av网站| 亚洲人成77777在线视频| 国产精品一区二区免费欧美| 国产精品久久久久久精品电影 | 亚洲精品国产一区二区精华液| 成人国产一区最新在线观看| 久9热在线精品视频| 激情在线观看视频在线高清| 日日干狠狠操夜夜爽| 成人精品一区二区免费| 99国产精品免费福利视频| 校园春色视频在线观看| 午夜福利18| 在线观看日韩欧美| 麻豆成人av在线观看| 一夜夜www| 男男h啪啪无遮挡| 99热只有精品国产| 国产人伦9x9x在线观看| 亚洲精品久久成人aⅴ小说| 性少妇av在线| 纯流量卡能插随身wifi吗| 中国美女看黄片| 久久久久国产精品人妻aⅴ院| 91成年电影在线观看| 久久人妻av系列| 精品国产一区二区三区四区第35| 国产单亲对白刺激| 亚洲美女黄片视频| 精品无人区乱码1区二区| 欧美日韩黄片免| 日本 欧美在线| 搞女人的毛片| 久久人妻福利社区极品人妻图片| 中文字幕久久专区| 国产精品二区激情视频| 欧美大码av| 久久狼人影院| 久久精品国产99精品国产亚洲性色 | 丁香欧美五月| 美女 人体艺术 gogo| 每晚都被弄得嗷嗷叫到高潮| 禁无遮挡网站| 欧美成人午夜精品| 日韩高清综合在线| 成人特级黄色片久久久久久久| 久久亚洲真实| 欧美 亚洲 国产 日韩一| 国产三级黄色录像| 一级毛片精品| 亚洲精品粉嫩美女一区| 国产1区2区3区精品| 88av欧美| 亚洲天堂国产精品一区在线| 亚洲成人国产一区在线观看| 两个人看的免费小视频| 国产成年人精品一区二区| 国产日韩一区二区三区精品不卡| 欧美老熟妇乱子伦牲交| 天堂影院成人在线观看| 国产不卡一卡二| 女人精品久久久久毛片| 国产伦一二天堂av在线观看| 亚洲九九香蕉| 精品久久久精品久久久| 日本免费a在线| 一本综合久久免费| 中文字幕另类日韩欧美亚洲嫩草| 两性午夜刺激爽爽歪歪视频在线观看 | 夜夜看夜夜爽夜夜摸| 亚洲五月天丁香| 757午夜福利合集在线观看| 久久中文字幕一级| 韩国av一区二区三区四区| 国产精品九九99| 亚洲情色 制服丝袜| 欧美午夜高清在线| 欧洲精品卡2卡3卡4卡5卡区| 日韩欧美一区视频在线观看| 色综合欧美亚洲国产小说| 国产男靠女视频免费网站| 国产成人一区二区三区免费视频网站| 91字幕亚洲| bbb黄色大片| 大型黄色视频在线免费观看| 日本撒尿小便嘘嘘汇集6| 亚洲伊人色综图| 国产麻豆69| 久久人人97超碰香蕉20202| 国产麻豆成人av免费视频| 亚洲全国av大片| 亚洲人成电影免费在线| 久久欧美精品欧美久久欧美| 女人被狂操c到高潮| 无遮挡黄片免费观看| 国产97色在线日韩免费| 国产又爽黄色视频| 国内久久婷婷六月综合欲色啪| 日韩欧美一区视频在线观看| 久久影院123| 1024香蕉在线观看| 日本三级黄在线观看| 在线观看免费视频网站a站| 色综合欧美亚洲国产小说| 一级a爱片免费观看的视频| 国产成人欧美在线观看| 一级a爱片免费观看的视频| www国产在线视频色| 日韩成人在线观看一区二区三区| 国产亚洲欧美精品永久| 午夜a级毛片| 亚洲欧美精品综合久久99| 国产不卡一卡二| 欧洲精品卡2卡3卡4卡5卡区| 国产精品免费一区二区三区在线| 一级毛片女人18水好多| 欧美日韩中文字幕国产精品一区二区三区 | 国产成人欧美| 国产欧美日韩综合在线一区二区| 国产精品久久久久久亚洲av鲁大| 亚洲人成伊人成综合网2020| 最好的美女福利视频网| 中文字幕色久视频| 亚洲第一av免费看| 国产成人欧美在线观看| 国产高清激情床上av| 欧美国产日韩亚洲一区| 一进一出好大好爽视频| 精品久久久久久久久久免费视频| 99在线人妻在线中文字幕| 99精品在免费线老司机午夜| 两个人免费观看高清视频| 桃色一区二区三区在线观看| 午夜激情av网站| 国产一区二区三区在线臀色熟女| 级片在线观看| 夜夜看夜夜爽夜夜摸| 啪啪无遮挡十八禁网站| 我的亚洲天堂| 亚洲av成人不卡在线观看播放网| 99riav亚洲国产免费| 在线观看舔阴道视频| 亚洲熟妇中文字幕五十中出| 亚洲人成伊人成综合网2020| 韩国av一区二区三区四区| 天堂√8在线中文| 久久人妻av系列| 丰满的人妻完整版| 色综合欧美亚洲国产小说| 人人妻人人澡欧美一区二区 | 亚洲国产欧美日韩在线播放| 禁无遮挡网站| 精品福利观看| 久久久久久国产a免费观看| 两个人看的免费小视频| 婷婷六月久久综合丁香| 久久久久国产一级毛片高清牌| 亚洲伊人色综图| 99久久国产精品久久久| 黄色成人免费大全| 满18在线观看网站| 一本综合久久免费| 久久久久久国产a免费观看| 级片在线观看| 日韩大码丰满熟妇| 亚洲中文日韩欧美视频| 亚洲激情在线av| www.熟女人妻精品国产| 亚洲久久久国产精品| 久久午夜综合久久蜜桃| 女生性感内裤真人,穿戴方法视频| 国产精品电影一区二区三区| 19禁男女啪啪无遮挡网站| 精品久久久久久久人妻蜜臀av | 制服诱惑二区| 91字幕亚洲| 三级毛片av免费| 亚洲欧美日韩无卡精品| 校园春色视频在线观看| 精品乱码久久久久久99久播| 精品久久久久久久人妻蜜臀av | 国产精品电影一区二区三区| 免费看美女性在线毛片视频| 97人妻天天添夜夜摸| 久久性视频一级片| 久久欧美精品欧美久久欧美| 一边摸一边抽搐一进一小说| 老司机深夜福利视频在线观看| 中文字幕人妻丝袜一区二区| 久久伊人香网站| 色在线成人网| 操出白浆在线播放| av免费在线观看网站| 国产精品自产拍在线观看55亚洲| 欧美国产日韩亚洲一区| 亚洲av美国av| 真人做人爱边吃奶动态| 亚洲av第一区精品v没综合| 一级黄色大片毛片| 欧美日韩一级在线毛片| 久久国产乱子伦精品免费另类| 欧洲精品卡2卡3卡4卡5卡区| 亚洲久久久国产精品| 韩国av一区二区三区四区| 国产成人精品久久二区二区91| 亚洲欧美激情在线| 黄色片一级片一级黄色片| 亚洲国产精品sss在线观看| 日本在线视频免费播放| 亚洲中文av在线| 国产亚洲精品一区二区www| 亚洲国产精品成人综合色| 亚洲国产高清在线一区二区三 | 久久草成人影院| 国产精品久久久久久精品电影 | 欧美一区二区精品小视频在线| 国产97色在线日韩免费| 精品国产美女av久久久久小说| 亚洲 欧美一区二区三区| 叶爱在线成人免费视频播放| 国产精品免费视频内射| 免费人成视频x8x8入口观看| 91在线观看av| 国产午夜精品久久久久久| 亚洲第一电影网av| 亚洲欧美日韩无卡精品| 日本一区二区免费在线视频| 夜夜看夜夜爽夜夜摸| 男女之事视频高清在线观看| 怎么达到女性高潮| 国产主播在线观看一区二区| 99久久精品国产亚洲精品| 少妇 在线观看| 国产亚洲精品久久久久5区| 窝窝影院91人妻| 88av欧美| 日本在线视频免费播放| 国产精品香港三级国产av潘金莲| 桃色一区二区三区在线观看| 悠悠久久av| cao死你这个sao货| 中文字幕av电影在线播放| 精品乱码久久久久久99久播| 国内久久婷婷六月综合欲色啪| 丁香欧美五月| 老司机午夜十八禁免费视频| 免费在线观看影片大全网站| 99riav亚洲国产免费| 777久久人妻少妇嫩草av网站| 国产伦人伦偷精品视频| 亚洲人成伊人成综合网2020| 午夜福利高清视频| 男人操女人黄网站| 黄网站色视频无遮挡免费观看| 又大又爽又粗| 精品福利观看| 色播亚洲综合网| 手机成人av网站| 久久亚洲真实| 亚洲欧美精品综合一区二区三区| 啪啪无遮挡十八禁网站| 少妇粗大呻吟视频| 成熟少妇高潮喷水视频| 女生性感内裤真人,穿戴方法视频| 成人国产综合亚洲| 国产麻豆成人av免费视频| 在线观看免费视频网站a站| 国产精品野战在线观看| 母亲3免费完整高清在线观看| 一区二区三区高清视频在线| 亚洲国产中文字幕在线视频| 成人亚洲精品一区在线观看| 免费无遮挡裸体视频| av福利片在线| avwww免费| 一级毛片精品| 国产精品 国内视频| 久久欧美精品欧美久久欧美| 久久午夜综合久久蜜桃| 国产亚洲精品久久久久5区| 亚洲精品国产区一区二| 欧美成人免费av一区二区三区| 十八禁网站免费在线| www日本在线高清视频| 欧美色视频一区免费| 一级a爱片免费观看的视频| 夜夜躁狠狠躁天天躁| 日日干狠狠操夜夜爽| 亚洲国产欧美日韩在线播放| 熟女少妇亚洲综合色aaa.| 日韩三级视频一区二区三区| 美女国产高潮福利片在线看| 国产国语露脸激情在线看| 久久精品aⅴ一区二区三区四区| 在线国产一区二区在线| 亚洲精品久久国产高清桃花| 欧美日本亚洲视频在线播放| 黑人欧美特级aaaaaa片| 欧美激情高清一区二区三区| 黄色 视频免费看| 色综合欧美亚洲国产小说| 国产亚洲欧美在线一区二区| 97人妻天天添夜夜摸| 美女 人体艺术 gogo| 欧美午夜高清在线| 又黄又粗又硬又大视频| 一区在线观看完整版| avwww免费| 国产成人精品久久二区二区免费| 中亚洲国语对白在线视频| 久久久久亚洲av毛片大全| 久久人妻av系列| 日本免费a在线| 欧美一区二区精品小视频在线| 欧美国产日韩亚洲一区| 欧美日韩亚洲综合一区二区三区_| 久久婷婷成人综合色麻豆| 美女扒开内裤让男人捅视频| 琪琪午夜伦伦电影理论片6080| 国产乱人伦免费视频| 在线观看午夜福利视频| 亚洲 国产 在线| 亚洲第一电影网av| 午夜亚洲福利在线播放| 激情视频va一区二区三区| 91麻豆av在线| videosex国产| 黄网站色视频无遮挡免费观看| 人妻久久中文字幕网| 亚洲av日韩精品久久久久久密| 亚洲成人免费电影在线观看| 人妻久久中文字幕网|