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

    強非線性杜芬系統(tǒng)的周期解及其分岔

    2017-10-19 09:16:07王旭李欣業(yè)王振靜韓善凱
    關(guān)鍵詞:受迫振動項數(shù)平均法

    王旭,李欣業(yè),王振靜,韓善凱

    (河北工業(yè)大學(xué) 機械工程學(xué)院,天津 300130)

    強非線性杜芬系統(tǒng)的周期解及其分岔

    王旭,李欣業(yè),王振靜,韓善凱

    (河北工業(yè)大學(xué) 機械工程學(xué)院,天津 300130)

    基于廣義諧波平衡法,求解了強非線性杜芬振子自由振動和簡諧激勵下受迫振動的周期-m解,并與數(shù)值解進行了比較,從而討論非線性項的系數(shù)以及激勵參數(shù)對系統(tǒng)周期解的影響.對自由振動而言,倍周期響應(yīng)的周期是派生系統(tǒng)固有周期的整倍數(shù);對受迫振動而言,倍周期響應(yīng)的周期是外激勵周期的整倍數(shù).結(jié)果表明,為使近似解析諧波解與數(shù)值解比較接近,系統(tǒng)的非線性越強,所需的諧波項數(shù)越多;所設(shè)倍周期分岔解的周期越大,所需的項數(shù)也越多.

    廣義諧波平衡法;強非線性;杜芬系統(tǒng);周期-m解

    眾所周知,在傳統(tǒng)的求解非線性振動系統(tǒng)的各種近似解析方法中,諧波平衡法是概念最明了、使用最簡便的,且不僅限于弱非線性系統(tǒng)[1-3].也正因如此,許多學(xué)者將諧波平衡法與其他方法或理論相結(jié)合,提出了新的改進方法,使諧波平衡法的理論有了十分豐富的發(fā)展[4].

    1981年Lau和Cheung[5]最先將諧波平衡法與增量法結(jié)合,并將該方法稱為增量諧波平衡法(IHB法).Yuste[6]將諧波平衡法與雅可比橢圓函數(shù)結(jié)合起來,得到橢圓函數(shù)諧波平衡法.Huseyin等[7]提出了一種把諧波平衡法與多尺度法結(jié)合起來的內(nèi)在的時間多尺度諧波平衡法來分析非線性周期振蕩和分岔問題,在此基礎(chǔ)上,Summers等[8]提出了雙時間尺度諧波平衡法.

    早在Lagrange時期,在確定行星運行軌道的演化規(guī)律時就提出了平均的概念,后經(jīng)Krylov和Bogolyubov提出著名的KB變換后,逐步形成了具有嚴(yán)格數(shù)學(xué)理論基礎(chǔ)的平均法[9].Chow[10]首次在分岔問題中使用了平均法,研究了二維中心流Hopf分岔的穩(wěn)定性交換公式.平均法經(jīng)過徐兆和張佑啟等多位學(xué)者的不斷完善,逐步發(fā)展為廣義平均法[11].在此基礎(chǔ)上又加入廣義諧波函數(shù)的概念,通過這種處理平均法不但可以用于求弱非線性問題,對強非線性系統(tǒng)也同樣適用.

    廣義諧波平衡法[12]綜合了諧波平衡法和平均法的優(yōu)勢,既是諧波平衡法的發(fā)展,又是平均法的發(fā)展.本文將基于此法構(gòu)造杜芬系統(tǒng)自由振動和受迫振動的近似解析周期解.由于篇幅限制,略去了對基于此法構(gòu)造的周期解的穩(wěn)定性分析.

    1 廣義諧波平衡法

    考慮如下形式的運動微分方程

    (1)

    與諧波平衡法類似,設(shè)其解的形式如下:

    (2)

    當(dāng)該方程為自治系統(tǒng)時,其基頻Ω是對應(yīng)的線性系統(tǒng)的固有頻率.當(dāng)該方程為非自治系統(tǒng)時,對應(yīng)的是外激勵周期的整數(shù)倍,與平均法類似,式(2)中的諸系數(shù)是時變的.對式(2)求二階導(dǎo)數(shù),可得

    (3)

    將式(3)代入方程(1),再利用三角函數(shù)的正交性可得

    (4)

    其中

    式(4)又可以改寫成如下一階微分方程組的形式:

    (5)

    ?

    ?

    式(5)是用于確定近似解(2)的平均方程,若令z=(a0,b1,…,bn,c1,…,cn)T,又可寫成如下形式:

    (6)

    與平均法類似,為求形如式(2)的穩(wěn)態(tài)解,令各平均方程等號右邊的表達式為零,可得到關(guān)于諸系數(shù)的2n+1個非線性代數(shù)方程組

    (7)

    一般地,式(7)可以通過Newton-Raphson方法進行數(shù)值求解.根據(jù)李雅普諾夫一次近似理論,可以通過求相應(yīng)雅克比矩陣的特征值來判斷解的穩(wěn)定性.

    2 杜芬系統(tǒng)自由振動的周期解

    在不考慮阻尼的情況下,杜芬系統(tǒng)的自由振動微分方程可以表示為

    (8)

    相應(yīng)的平均方程為

    (9)

    圖1和圖2分別表示在α、β取不同值時,x的近似解析解和數(shù)值解的時間歷程曲線,其中圖例xn表示數(shù)值解,xa表示近似解析解,下角標(biāo)中的數(shù)字代表諧波項的項數(shù).

    圖1 α=1,β=0.1時的時間歷程Fig.1 Time history with α=1,β=0.1

    圖2 α=4,β=5時的時間歷程Fig.2 Time history with α=4,β=5

    圖1表明,當(dāng)非線性較弱時,用3項諧波的和作為近似解就能較好地與數(shù)值解吻合,繼續(xù)增加諧波項的項數(shù),對提高解的精度意義不大.圖2表明,當(dāng)系統(tǒng)的非線性較強時,由于派生系統(tǒng)的固有周期變小,用10項諧波的和作為近似解也能較好地與數(shù)值解吻合.所以在接下來討論杜芬系統(tǒng)的受迫振動時,近似解是按至少包含10項諧波給出的.

    3 杜芬系統(tǒng)受迫振動的周期解

    在自由振動基礎(chǔ)上,以杜芬系統(tǒng)的受迫振動為例,驗證廣義諧波平衡法對周期-m解的適用性.

    3.1 周期-m解

    杜芬系統(tǒng)的受迫振動微分方程可以表示為

    (10)

    設(shè)其解為

    (11)

    此時的平均方程為

    n.

    (12)

    將式(12)簡寫成

    (13)

    g(m)(z(m))=06×1.

    (14)

    此時得到的平衡點即為近似解析解,即方程(10)中動力學(xué)系統(tǒng)的穩(wěn)態(tài)解.

    3.2 激勵參數(shù)的影響

    強非線性杜芬系統(tǒng)在作受迫振動時可能出現(xiàn)倍周期分岔的情況,以隨外激勵幅值和外激勵頻率變化而產(chǎn)生的倍周期分岔為例,對廣義諧波平衡法求解周期-m解的適用性進行分析.

    3.2.1 激勵幅值的影響

    圖3和圖4分別為μ=0.3,α=-1,β=1,Ω=1.2時,不同的F值對應(yīng)的x的時間歷程曲線.

    圖3 F=0.2時的時間歷程Fig.3 Time history with F=0.2

    圖4 F=0.287時的時間歷程Fig.4 Time history with F=0.287

    與圖3對比,圖4表明系統(tǒng)的受迫振動響應(yīng)周期為圖3所示響應(yīng)周期的2倍.

    3.2.2 激勵頻率的影響

    圖5、6、7分別表示方程(9)中選取μ=0.25、α=-1、β=1、F=0.26時隨Ω值的變化,響應(yīng)出現(xiàn)倍周期分岔時,x的時間歷程曲線.

    圖5—圖7表明,隨著外激勵頻率的減小,受迫振動響應(yīng)出現(xiàn)了倍周期分岔,即響應(yīng)周期倍增的現(xiàn)象.由于周期不斷變大,所以近似解中的諧波項數(shù)越來越多,以期能與數(shù)值解有比較好地吻合.

    圖5 Ω=1.3時的時間歷程Fig.5 Time history with Ω=1.3

    圖6 Ω=1.24時的時間歷程Fig.6 Time history with Ω=1.24

    圖7 Ω=1.2時的時間歷程Fig.7 Time history with Ω=1.2

    4 討論與結(jié)論

    廣義諧波平衡法結(jié)合了諧波平衡法和平均法的優(yōu)點,即不僅可以進行解的穩(wěn)定性分析,又適用于強非線性系統(tǒng),亦適用于多自由度系統(tǒng).本文利用此法求解強非線性杜芬系統(tǒng)的近似周期解時,僅考慮了倍周期分岔的情況,即對自由振動,僅考慮了周期為派生系統(tǒng)固有周期整倍數(shù)的情況,對受迫振動僅考慮了周期為外激勵周期整倍數(shù)的情況.對于受迫振動而言,這種倍周期運動相當(dāng)于亞諧振動.事實上,也可用此法求強非線性系統(tǒng)的超諧振動解,即解的周期為派生系統(tǒng)固有周期(對自由振動)或激勵周期(對受迫振動)的幾分之一倍的情況.本文的結(jié)果表明,不論是自由振動還是受迫振動,只要所需諧波項數(shù)足夠多,基于此法構(gòu)造的近似解析解都能很好地與數(shù)值解吻合.為使近似解析諧波解與數(shù)值解比較接近,則系統(tǒng)的非線性越強,所需的諧波項數(shù)越多;所設(shè)倍周期解的周期越大,所需的項數(shù)也越多.

    [1] 陳予恕.非線性振動[M].北京: 高等教育出版社,2002.

    [2] 陳予恕.非線性振動、分岔和混沌理論及其應(yīng)用[J].振動工程學(xué)報,1992,5(3): 235-250.DOI:10.16385/j.cnki.issn.1004-4523.1992.03.009.

    [3] 劉延柱,陳立群.非線性振動[M].北京: 高等教育出版社,2001.

    [4] CHEN Y M,LIU J K.A new method based on the harmonic balance method for nonlinear oscillator [J].Physics Letters A,2007,368(5): 371-378.DOI:10.1016/j.physleta.2007.04.025.

    [5] LAU S L,CHEUNG Y K.Amplitude incremental variational principle for nonlinear of elastic systems[J].Journal of Applied Mechanies,1981,48: 959-964.DOI: 10.1115/1.3157762.

    [6] YUSTE S B.Comments on the method of harmonic-balance in which Jacobi elliptic function are used[J].Journal of Sound and Vibration,1991,145(3): 381-390.DOI: 10.1016/0022-460X(91)90109-W.

    [7] HUSEYIN K,LIN R.An intrinsic multiple-scale harmonic balance method for nonlinear vibration and bifurcation problems[J].International Journal of Non-linear Mechanics,1991,26: 727-740.DOI: 10.1016/0020-7462(91)90023-M.

    [8] SUMMERS J L,SAVAGE M D.2 timescale harmonic-balance.I.Application to autonomous one-dimensional nonlinear oscillators[J].Philosophical Transactions of the Royal Society A,1992,3: 340-473.DOI: 10.1098/rsta.1992.0077.

    [9] 張琪昌,王洪禮.分岔與混沌理論及應(yīng)用[M].天津: 天津大學(xué)出版社,2005.

    [10] CHOW,S N,MALLER-PARET.Integral averaging and bifurcations[J].Differential Equations,1977,26: 112- 159.DOI: 10.1016/0022-0396(77)90101-2.

    [11] CHEUNG Y K,XU Z.Intemal resonance of strongly non-linear autonomous vibrating systems with many degrees of freedom[J].Journal of Sound and Vibration,1995,180(2): 229-238.DOI:10.1006/jsvi.1995.0076.

    [12] ALBERT C,LUO J.ARASH BAGHAEI LAKEH.Period-m motions and bifurcation trees in a periodically forced van der Pol-Duffing oscillator[J].International Journal of Dynamics and Control,2014,2(4):474-493.DOI: 10.1007/s40435-014-0058-9.

    (責(zé)任編輯:王蘭英)

    PeriodicsolutionsandtheirbifurcationofDuffingsystemwithstrongnon-linearity

    WANGXu,LIXinye,WANGZhenjing,HANShankai

    (School of Mechanical Engineering,Hebei University of Technology,Tianjin 300130,China)

    Based on the generalized harmonic balance method,the period-m responses of free and forced Duffing oscillators are constructed and compared with numerical solutions from which the effects of nonlinear term and excitation parameters can be observed.For free vibrations,the response periods are assumed to be integer times of the natural period of the corresponding linear system and the excitation period respectively.It is shown that the number of harmonic terms should be large enough for both strongly nonlinear systems and period-doubling bifurcation solutions with large period so that the approximate analytic solutions are comparative to numerical results.

    generalized harmonic balance method;strong non-linearity;Duffing systems;period-msolution

    O322

    A

    1000-1565(2017)05-0457-07

    10.3969/j.issn.1000-1565.2017.05.003

    2017-03-27

    河北省高層次人才資助項目(C201400309)

    王旭(1990—),男,河北保定人,河北工業(yè)大學(xué)在讀碩士研究生.E-mail:1216380802@qq.com

    李欣業(yè)(1966—),男,河北遷安人,河北工業(yè)大學(xué)教授,主要從事復(fù)雜系統(tǒng)的動力學(xué)建模與分析和非線性動力學(xué)與控制研究.E-mail:xylihebut@163.com

    猜你喜歡
    受迫振動項數(shù)平均法
    基于平均法的金屬橡膠隔振器非線性振動特性研究
    等比數(shù)列的性質(zhì)、推論和應(yīng)用
    具有初邊值條件的集值脈沖微分方程的平均法
    利用自制教具探究受迫振動與共振
    平均法處理自由落體頻閃數(shù)據(jù)的研究
    物理教師(2017年5期)2017-06-09 11:21:18
    求 和
    論高次方程
    《推理與證明》必考題型賞析
    永磁行星齒輪傳動系統(tǒng)受迫振動分析
    有限差分法模擬電梯懸掛系統(tǒng)橫向受迫振動
    美女免费视频网站| 真实男女啪啪啪动态图| 91麻豆av在线| 老师上课跳d突然被开到最大视频| 国内精品久久久久精免费| 精品国内亚洲2022精品成人| 色综合婷婷激情| 欧美国产日韩亚洲一区| 国产av麻豆久久久久久久| 在线观看舔阴道视频| 欧美精品啪啪一区二区三区| 亚洲专区国产一区二区| 国产女主播在线喷水免费视频网站 | 国产精品久久久久久亚洲av鲁大| 午夜a级毛片| 男女做爰动态图高潮gif福利片| 国产高潮美女av| 日本五十路高清| 国产午夜精品久久久久久一区二区三区 | 日本三级黄在线观看| 亚洲avbb在线观看| 欧美潮喷喷水| 内射极品少妇av片p| 国产黄片美女视频| 变态另类成人亚洲欧美熟女| 12—13女人毛片做爰片一| 村上凉子中文字幕在线| 精品午夜福利在线看| 国产精品日韩av在线免费观看| 欧美性感艳星| 亚洲精品久久国产高清桃花| 亚洲一区高清亚洲精品| 久久精品国产亚洲av天美| 国产亚洲精品综合一区在线观看| 国语自产精品视频在线第100页| 国产视频一区二区在线看| 精品一区二区三区人妻视频| 美女大奶头视频| 免费av毛片视频| 欧美日韩中文字幕国产精品一区二区三区| 一个人观看的视频www高清免费观看| 在线免费十八禁| 成年免费大片在线观看| 亚洲欧美日韩高清在线视频| 欧美日韩瑟瑟在线播放| 久久热精品热| 男女做爰动态图高潮gif福利片| 欧美区成人在线视频| 老熟妇仑乱视频hdxx| 人妻夜夜爽99麻豆av| 久久久久久伊人网av| 欧美日韩中文字幕国产精品一区二区三区| 久久午夜福利片| 成人av在线播放网站| 一级av片app| 中文字幕人妻熟人妻熟丝袜美| 国产精品电影一区二区三区| 一个人看的www免费观看视频| 国产一级毛片七仙女欲春2| 亚洲专区中文字幕在线| 日韩精品有码人妻一区| 长腿黑丝高跟| 观看美女的网站| 长腿黑丝高跟| 国产成人福利小说| 欧美精品国产亚洲| 亚洲av美国av| 亚洲av.av天堂| 国产男靠女视频免费网站| 国产高清不卡午夜福利| 国产乱人伦免费视频| 中文字幕av在线有码专区| 91麻豆av在线| 成熟少妇高潮喷水视频| 欧美极品一区二区三区四区| 日本撒尿小便嘘嘘汇集6| 毛片女人毛片| 春色校园在线视频观看| 日韩高清综合在线| 免费人成视频x8x8入口观看| 成人午夜高清在线视频| 精品免费久久久久久久清纯| 欧美日韩乱码在线| 极品教师在线免费播放| 舔av片在线| 成人国产一区最新在线观看| 免费电影在线观看免费观看| 亚洲va在线va天堂va国产| 三级毛片av免费| 国产麻豆成人av免费视频| 蜜桃亚洲精品一区二区三区| 国产激情偷乱视频一区二区| 国产熟女欧美一区二区| 欧美bdsm另类| 亚洲精品色激情综合| 亚洲狠狠婷婷综合久久图片| 欧美成人一区二区免费高清观看| 成人二区视频| 中文字幕av在线有码专区| 亚洲专区国产一区二区| 2021天堂中文幕一二区在线观| 国产成年人精品一区二区| 亚洲av第一区精品v没综合| 亚洲av免费高清在线观看| 久久久久性生活片| 欧美另类亚洲清纯唯美| 波野结衣二区三区在线| 黄片wwwwww| 麻豆国产97在线/欧美| 色播亚洲综合网| 午夜激情福利司机影院| 国产精品电影一区二区三区| 变态另类成人亚洲欧美熟女| eeuss影院久久| 色综合站精品国产| 午夜精品在线福利| 午夜日韩欧美国产| 自拍偷自拍亚洲精品老妇| 久久久久久久精品精品| 欧美成人午夜免费资源| 成人亚洲精品一区在线观看 | 黑人高潮一二区| av国产精品久久久久影院| 精品一区二区三区视频在线| 欧美成人精品欧美一级黄| 秋霞伦理黄片| 午夜福利网站1000一区二区三区| 亚洲av男天堂| av在线播放精品| 少妇人妻精品综合一区二区| 中文字幕人妻熟人妻熟丝袜美| 十八禁网站网址无遮挡 | 欧美变态另类bdsm刘玥| 亚洲欧美精品自产自拍| 亚洲欧美成人综合另类久久久| 搡老乐熟女国产| 国产亚洲精品久久久com| 午夜精品国产一区二区电影| 熟女av电影| 成人毛片60女人毛片免费| 久久久久久伊人网av| 国精品久久久久久国模美| 少妇人妻精品综合一区二区| 成人综合一区亚洲| 我的老师免费观看完整版| 日韩制服骚丝袜av| 欧美成人a在线观看| 最近的中文字幕免费完整| 亚洲av成人精品一区久久| av福利片在线观看| 精华霜和精华液先用哪个| av一本久久久久| 毛片女人毛片| 夜夜爽夜夜爽视频| 久久人妻熟女aⅴ| 91aial.com中文字幕在线观看| 亚洲综合色惰| 国产伦精品一区二区三区视频9| 最新中文字幕久久久久| 丰满人妻一区二区三区视频av| 国产乱人视频| 国产精品福利在线免费观看| 日韩中字成人| 国产成人91sexporn| 97超视频在线观看视频| 亚洲自偷自拍三级| 久久人人爽av亚洲精品天堂 | 久久久久久人妻| 亚洲av综合色区一区| 国产亚洲精品久久久com| 日本欧美国产在线视频| 水蜜桃什么品种好| 国产视频内射| 国产亚洲一区二区精品| 高清在线视频一区二区三区| 国产 一区 欧美 日韩| 久久精品国产亚洲av涩爱| 日韩人妻高清精品专区| 2018国产大陆天天弄谢| 熟女人妻精品中文字幕| 在线观看人妻少妇| av天堂中文字幕网| 最新中文字幕久久久久| 一区二区三区免费毛片| 国产欧美日韩一区二区三区在线 | av又黄又爽大尺度在线免费看| 亚洲内射少妇av| 久久午夜福利片| 国产精品久久久久久久久免| 熟女av电影| 中文字幕av成人在线电影| 日本色播在线视频| 亚洲一级一片aⅴ在线观看| 国产亚洲欧美精品永久| 色婷婷久久久亚洲欧美| 成人美女网站在线观看视频| 久久精品人妻少妇| 黑人高潮一二区| 亚洲四区av| 高清在线视频一区二区三区| av播播在线观看一区| 国产黄片视频在线免费观看| 晚上一个人看的免费电影| 久久毛片免费看一区二区三区| 久久热精品热| 99九九线精品视频在线观看视频| a级毛色黄片| 人妻 亚洲 视频| 大话2 男鬼变身卡| 免费久久久久久久精品成人欧美视频 | 99精国产麻豆久久婷婷| 精品人妻偷拍中文字幕| 日本av免费视频播放| 色婷婷av一区二区三区视频| 日韩制服骚丝袜av| 七月丁香在线播放| 色哟哟·www| 精品久久国产蜜桃| 各种免费的搞黄视频| 九色成人免费人妻av| 男女国产视频网站| 91久久精品国产一区二区三区| 久久久精品94久久精品| 欧美少妇被猛烈插入视频| 丰满迷人的少妇在线观看| 国内精品宾馆在线| 日韩一区二区视频免费看| 男女边摸边吃奶| 搡老乐熟女国产| 欧美丝袜亚洲另类| 成人美女网站在线观看视频| 亚洲av中文av极速乱| 多毛熟女@视频| 一本一本综合久久| 联通29元200g的流量卡| 男人舔奶头视频| 一本色道久久久久久精品综合| 免费观看无遮挡的男女| 久久精品国产亚洲av天美| 欧美最新免费一区二区三区| 亚洲国产av新网站| 六月丁香七月| 中文在线观看免费www的网站| 国产在线视频一区二区| 日韩不卡一区二区三区视频在线| 美女cb高潮喷水在线观看| 国产一区有黄有色的免费视频| 精品国产露脸久久av麻豆| 国产成人免费观看mmmm| 色婷婷久久久亚洲欧美| 赤兔流量卡办理| 91午夜精品亚洲一区二区三区| 国产精品免费大片| 亚洲欧美一区二区三区国产| 国产成人午夜福利电影在线观看| 国产亚洲5aaaaa淫片| 久久影院123| 麻豆精品久久久久久蜜桃| 熟女av电影| 久久国产乱子免费精品| 美女福利国产在线 | 日韩欧美一区视频在线观看 | 男女边摸边吃奶| 91久久精品国产一区二区三区| 狂野欧美激情性bbbbbb| 亚洲性久久影院| 午夜日本视频在线| 国产精品欧美亚洲77777| 99热全是精品| 亚洲无线观看免费| 亚洲怡红院男人天堂| 91精品国产国语对白视频| 一级爰片在线观看| 亚洲真实伦在线观看| 黄色怎么调成土黄色| 丝瓜视频免费看黄片| 看十八女毛片水多多多| 亚洲色图av天堂| 亚洲av欧美aⅴ国产| 多毛熟女@视频| 国产精品国产三级国产专区5o| 精品少妇黑人巨大在线播放| 中文欧美无线码| 国产黄片视频在线免费观看| 蜜桃久久精品国产亚洲av| 美女脱内裤让男人舔精品视频| 亚洲国产欧美在线一区| 亚洲欧美一区二区三区国产| 国产精品一区二区三区四区免费观看| 久久99热这里只有精品18| 日韩欧美精品免费久久| 日韩电影二区| 亚洲熟女精品中文字幕| 在线观看国产h片| 少妇熟女欧美另类| 性高湖久久久久久久久免费观看| 久久久亚洲精品成人影院| 国产国拍精品亚洲av在线观看| 国产黄色视频一区二区在线观看| 免费看光身美女| 久久久久久伊人网av| 中国国产av一级| 国产免费又黄又爽又色| 国产黄频视频在线观看| 欧美成人一区二区免费高清观看| 18+在线观看网站| 国产精品福利在线免费观看| 国产精品久久久久久av不卡| 欧美日韩亚洲高清精品| 男女免费视频国产| 搡老乐熟女国产| 九九爱精品视频在线观看| 国产成人精品久久久久久| 丰满乱子伦码专区| 国产91av在线免费观看| 99热国产这里只有精品6| 性色av一级| 亚洲激情五月婷婷啪啪| 最近的中文字幕免费完整| 亚洲av日韩在线播放| 一级毛片我不卡| 少妇高潮的动态图| 久久久久久久精品精品| 久久精品国产亚洲av涩爱| 狂野欧美激情性xxxx在线观看| 久久久久久久亚洲中文字幕| 久久久久视频综合| 一个人免费看片子| 国产 一区精品| 久久午夜福利片| 蜜桃亚洲精品一区二区三区| 国产一区有黄有色的免费视频| 久久精品人妻少妇| 久久影院123| av视频免费观看在线观看| 欧美日韩精品成人综合77777| 亚洲va在线va天堂va国产| 国产精品国产av在线观看| 免费播放大片免费观看视频在线观看| 欧美成人一区二区免费高清观看| 欧美老熟妇乱子伦牲交| 97超视频在线观看视频| 日韩 亚洲 欧美在线| 欧美极品一区二区三区四区| 两个人的视频大全免费| 久久久精品94久久精品| 一级av片app| 中文资源天堂在线| 欧美精品一区二区大全| 亚洲精品亚洲一区二区| 综合色丁香网| 久久精品熟女亚洲av麻豆精品| 国产一区二区三区av在线| 一级毛片我不卡| 18+在线观看网站| 如何舔出高潮| 身体一侧抽搐| 精品少妇久久久久久888优播| 涩涩av久久男人的天堂| 人人妻人人澡人人爽人人夜夜| 国产一区亚洲一区在线观看| 成人二区视频| 我的老师免费观看完整版| 欧美性感艳星| 精品99又大又爽又粗少妇毛片| kizo精华| 日韩三级伦理在线观看| 久热久热在线精品观看| 三级国产精品欧美在线观看| 国产高清有码在线观看视频| h视频一区二区三区| 中国美白少妇内射xxxbb| 中国三级夫妇交换| 最近的中文字幕免费完整| 国产精品一区二区在线观看99| 麻豆国产97在线/欧美| 欧美高清性xxxxhd video| 日韩av在线免费看完整版不卡| 少妇猛男粗大的猛烈进出视频| 成人午夜精彩视频在线观看| av免费在线看不卡| 男人狂女人下面高潮的视频| 三级经典国产精品| 国产精品欧美亚洲77777| 女性被躁到高潮视频| 日本黄色日本黄色录像| 麻豆国产97在线/欧美| 久久久久久久大尺度免费视频| 肉色欧美久久久久久久蜜桃| 在线观看美女被高潮喷水网站| 丝袜喷水一区| 国产 一区精品| 天堂8中文在线网| 亚洲av综合色区一区| 成人毛片60女人毛片免费| av国产精品久久久久影院| 日韩欧美一区视频在线观看 | 日韩国内少妇激情av| xxx大片免费视频| 六月丁香七月| 国产毛片在线视频| 亚洲欧洲国产日韩| 国产精品一二三区在线看| 久久鲁丝午夜福利片| 免费久久久久久久精品成人欧美视频 | 国语对白做爰xxxⅹ性视频网站| 国产久久久一区二区三区| 国产乱人视频| 99久久精品国产国产毛片| 国产欧美亚洲国产| 99九九线精品视频在线观看视频| 亚洲国产精品国产精品| 国产精品不卡视频一区二区| 只有这里有精品99| 久久99热这里只频精品6学生| 成人综合一区亚洲| 午夜福利在线在线| 国产免费一区二区三区四区乱码| 久久亚洲国产成人精品v| 国产成人91sexporn| 日本黄色片子视频| 亚洲精品aⅴ在线观看| 久久久国产一区二区| 尤物成人国产欧美一区二区三区| 国产极品天堂在线| 中文精品一卡2卡3卡4更新| 亚洲精品乱码久久久v下载方式| 男女边摸边吃奶| 天堂8中文在线网| 日韩欧美一区视频在线观看 | 国产欧美日韩一区二区三区在线 | 国产一区二区在线观看日韩| 麻豆乱淫一区二区| 国产女主播在线喷水免费视频网站| 色视频在线一区二区三区| 精品视频人人做人人爽| 成人国产麻豆网| 妹子高潮喷水视频| 国产在线视频一区二区| 久久久久人妻精品一区果冻| 在线观看一区二区三区激情| 日韩精品有码人妻一区| 国产精品伦人一区二区| 亚洲精品日韩在线中文字幕| 日韩亚洲欧美综合| 欧美日韩综合久久久久久| 精品一区在线观看国产| 性色av一级| 男女啪啪激烈高潮av片| 国内精品宾馆在线| 赤兔流量卡办理| 18禁裸乳无遮挡动漫免费视频| 亚洲av综合色区一区| 尾随美女入室| 欧美国产精品一级二级三级 | 中国三级夫妇交换| 高清毛片免费看| 中文欧美无线码| 久久精品夜色国产| 91久久精品电影网| 亚洲四区av| 中文资源天堂在线| 亚洲欧洲日产国产| 涩涩av久久男人的天堂| 国产成人免费观看mmmm| 如何舔出高潮| 婷婷色综合大香蕉| 在线观看人妻少妇| 日韩av免费高清视频| 男女无遮挡免费网站观看| 亚洲激情五月婷婷啪啪| 天堂8中文在线网| 日韩一本色道免费dvd| 国产视频内射| 99热国产这里只有精品6| 久热久热在线精品观看| 欧美国产精品一级二级三级 | 三级国产精品片| 日韩伦理黄色片| 精品人妻一区二区三区麻豆| 一级a做视频免费观看| 国产精品99久久久久久久久| 王馨瑶露胸无遮挡在线观看| 日韩欧美 国产精品| 最近手机中文字幕大全| freevideosex欧美| 亚洲精品国产av蜜桃| 国产精品成人在线| 亚洲精品成人av观看孕妇| 国产爱豆传媒在线观看| 在线免费十八禁| 91aial.com中文字幕在线观看| 一本久久精品| 久久人妻熟女aⅴ| 王馨瑶露胸无遮挡在线观看| tube8黄色片| 日韩制服骚丝袜av| 亚洲综合精品二区| 日韩av不卡免费在线播放| 夜夜骑夜夜射夜夜干| 国内少妇人妻偷人精品xxx网站| 五月天丁香电影| 国产午夜精品久久久久久一区二区三区| 狂野欧美激情性bbbbbb| 成人毛片60女人毛片免费| 夜夜骑夜夜射夜夜干| 国产一区二区三区综合在线观看 | 老熟女久久久| 国产亚洲欧美精品永久| 国产色爽女视频免费观看| 噜噜噜噜噜久久久久久91| 另类亚洲欧美激情| 国产真实伦视频高清在线观看| 欧美日韩亚洲高清精品| 亚洲熟女精品中文字幕| 免费少妇av软件| 精品国产三级普通话版| a级一级毛片免费在线观看| 在线观看三级黄色| 亚洲国产精品999| 激情 狠狠 欧美| 免费av不卡在线播放| 国产高清有码在线观看视频| 人妻夜夜爽99麻豆av| 2021少妇久久久久久久久久久| 国产黄片美女视频| 久久ye,这里只有精品| 国产精品99久久99久久久不卡 | 直男gayav资源| 亚洲欧洲国产日韩| 国产av一区二区精品久久 | 狠狠精品人妻久久久久久综合| 一区在线观看完整版| 免费少妇av软件| 各种免费的搞黄视频| 高清视频免费观看一区二区| 在线观看国产h片| 在线观看人妻少妇| 只有这里有精品99| 国产黄色视频一区二区在线观看| 亚洲av成人精品一二三区| 99热这里只有精品一区| 国国产精品蜜臀av免费| 成人一区二区视频在线观看| 欧美zozozo另类| 国产成人freesex在线| 18禁裸乳无遮挡免费网站照片| 久久久久久久久久人人人人人人| 国产中年淑女户外野战色| 国产黄色免费在线视频| 国产免费一区二区三区四区乱码| 午夜福利在线在线| 麻豆精品久久久久久蜜桃| 色吧在线观看| videos熟女内射| 各种免费的搞黄视频| 成年免费大片在线观看| 亚洲自偷自拍三级| 国产成人a区在线观看| 久久99热这里只频精品6学生| 91狼人影院| 国产精品熟女久久久久浪| 欧美精品一区二区大全| 亚洲美女搞黄在线观看| 日本wwww免费看| 国产午夜精品一二区理论片| 18禁在线播放成人免费| 日本黄大片高清| 亚洲四区av| 亚洲天堂av无毛| 人妻系列 视频| 身体一侧抽搐| 国产人妻一区二区三区在| 成年美女黄网站色视频大全免费 | 亚洲精品视频女| 久久久久久人妻| 少妇丰满av| 18+在线观看网站| 久久ye,这里只有精品| 美女脱内裤让男人舔精品视频| 毛片女人毛片| 99热国产这里只有精品6| 亚洲精品乱码久久久v下载方式| 国产亚洲最大av| 尾随美女入室| 免费观看的影片在线观看| 亚洲第一区二区三区不卡| 亚洲精品中文字幕在线视频 | 色网站视频免费| 亚洲国产最新在线播放| 在线观看免费日韩欧美大片 | 国产免费福利视频在线观看| 久久精品久久久久久久性| 日本av免费视频播放| 亚洲丝袜综合中文字幕| 亚洲精品国产成人久久av| 欧美日韩一区二区视频在线观看视频在线| 国内精品宾馆在线| 最新中文字幕久久久久| 亚洲国产精品国产精品| 色视频在线一区二区三区| 免费大片黄手机在线观看| 国产成人aa在线观看| 蜜臀久久99精品久久宅男| 夜夜爽夜夜爽视频| 十分钟在线观看高清视频www | 亚洲av日韩在线播放| 国产中年淑女户外野战色| 国产欧美另类精品又又久久亚洲欧美| 夜夜爽夜夜爽视频| 18禁在线无遮挡免费观看视频| 亚洲久久久国产精品| 下体分泌物呈黄色| 日产精品乱码卡一卡2卡三| 国产色婷婷99|