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

    受迫振動非線性特性的教學拓展

    2021-01-25 03:58:44丁紅勝王佳曦張師平董軍軍
    大學物理 2021年2期
    關鍵詞:擺輪受迫振動振幅

    丁紅勝,王佳曦,張師平,董軍軍

    (1. 北京科技大學 數理學院應用物理系,北京 100083;2. 北京科技大學 自然科學基礎實驗中心,北京 100083)

    振動是一種普遍存在于客觀世界的物質運動形式,系統(tǒng)在周期性外力的持續(xù)作用下的振動即為受迫振動[1].波耳共振儀是高校物理實驗教學中研究受迫振動運動規(guī)律較為普遍的儀器,也是學者探討較多的實驗項目[2-4].實驗測量中,儀器的軸承摩擦和彈簧非線性效應的影響是造成實驗數據偏離預期的主要因素.如果在受迫振動的分析中考慮這兩種因素的影響,則非常有助于學生理解實驗原理和結果偏差.本文基于受迫振動常規(guī)實驗,延伸探討其非線性特性,通過Matlab數值分析的方法來研究杜芬方程的混沌特性,以此提高學生對受迫振動的認知,并加深對其非線性特性的理解.

    1 受迫振動

    1.1 受迫振動方程

    波耳共振儀[5](如圖1所示)的擺輪可在彈性力矩作用下自由擺動,若同時加上阻尼力矩和驅動力矩,擺輪可做受迫振動.

    1. 光電門A;2. 長凹槽;3. 短凹槽;4. 銅制擺輪;5. 搖桿; 6. 蝸卷彈簧;7. 機架;8. 阻尼線圈;9. 連桿;10. 搖桿調節(jié)螺絲; 11. 光電門B;12. 角度盤;13. 有機玻璃轉盤;14. 底座; 15. 彈簧夾持螺絲;16閃光燈圖1 波耳共振儀的結構示意圖[5]

    (1)

    式中,J為擺輪的轉動慣量,-kθ為彈性力矩,k為彈簧的勁度系數,M0為驅動力矩的幅值,ω為驅動力的角頻率.

    (2)

    當mcosωt=0時,即在無周期性外力矩作用時,式(2)為阻尼振動方程;當阻尼系數為零,即無阻尼時,式(2)為簡諧振動方程,ω0即為振動系統(tǒng)的固有頻率.

    方程(2)的通解為

    θ=θ1e-βtcos(ω1t+α)+θ2cos(ωt+φ)

    (3)

    通解含有兩個部分:θ1e-βtcos(ω1t+α)表示阻尼振動,隨著時間的演化,振動會逐漸衰減,直至忽略不計;θ2cos(ωt+φ)表示簡諧振動,驅動力矩持續(xù)對擺輪做功,振動系統(tǒng)接收驅動力所傳送的能量,使振動系統(tǒng)最終達到一個穩(wěn)定的狀態(tài).

    振幅為定值,

    (4)

    系統(tǒng)相位差為

    (5)

    由式(4)和式(5)可看出,穩(wěn)定振動時,驅動力的幅值M0、阻尼系數β、驅動力的頻率ω(或驅動力周期T)和系統(tǒng)固有頻率ω0(或固有周期T0)4個參量共同決定系統(tǒng)振幅和相位差的大小,與振動的起始狀態(tài)無關.

    (6)

    可以看出,阻尼系數越小,驅動力的角頻率越接近系統(tǒng)的固有角頻率,振動幅值也越大.

    1.2 阻尼系數、固有周期、固有頻率的測量

    實驗中通常選擇阻尼開關位置為2,將波耳共振儀有機玻璃盤上零度標志線對準0°,用手將擺輪轉動至140°~150°之間,放手后控制箱會自動連續(xù)記錄擺輪完成阻尼振動10次的振幅θ0~θ9.數據如表1所示.

    表1 阻尼系數測量數據

    選擇“自由振蕩”進行測量,將有機玻璃盤上零度標志線保持在0°,用手將擺輪轉動至140°~150°之間,放手后記錄每次振動振幅值及其相應周期,并計算相應頻率值.數據如表2所示.

    表2 振幅、周期、頻率關系數據

    2 受迫振動方程非線性修正

    2.1 彈簧非線性效應

    一般情況下,可以認為彈簧的勁度系數k為常數,但是在上述彈簧自由振動的實驗中,通常隨著振幅的減小,振動系統(tǒng)固有周期會逐漸增大,這意味著振動系統(tǒng)的彈性回復力隨著振幅的減小而減弱[2].因此不妨假設彈簧的勁度系數為振幅的函數k=k(θ),則回復力矩為

    Fk=k(θ)θ

    (7)

    將Fk=k(θ)θ進行泰勒展開至一階項,

    (8)

    (9)

    2.2 軸承摩擦效應

    軸承摩擦效應普遍存在于定軸轉動的運動中,當軸承轉動速度較小時軸承靜摩擦力較大,而當轉速達到一定臨界值后,軸承摩擦力將保持穩(wěn)定[6].

    實驗發(fā)現,軸承摩擦力隨軸承角速度的變化率較緩,因此可采用二次函數的形式來描述軸承摩擦力矩隨角速度的關系

    (10)

    其中,固定系數B>0,C>0,在此情況下,即表明摩擦力矩隨角速度的增加而緩慢降低.

    2.3 阻尼振動方程非線性修正

    綜合彈簧的非線性效應與軸承摩擦效應對振動方程的修正,可以得到修正后的阻尼振動方程為

    (11)

    (12)

    式(12)即為非線性項修正后的阻尼振動方程,是一個待定系數的二階微分方程.

    在此情況下,基于Matlab軟件,用數值方法計算其數值解,通過多次給定a、b、c的數值進行模擬計算,擬合阻尼振動實驗所測得的數據.

    (13)

    利用Matlab軟件內置的ode45函數,根據四階龍格-庫塔法,多次變步長數值積分,可以得到,當待定系數滿足a=0.014,b=0.002,c=0.005時,數值方法解得的阻尼振動方程與實驗數據擬合的最為接近,特別是在振幅為90°~140°時,彈簧形變較大,彈簧非線性效應較為明顯的區(qū)域.

    圖2中指數曲線為阻尼振動的振幅-時間實測數據的擬合曲線,周期性曲線為非線性微分方程的數值解.

    圖2 阻尼振動數值積分與擬合曲線

    2.4 受迫振動方程非線性修正

    在修正后的阻尼振動方程(12)的右邊加上一個周期性力矩,即為受迫振動方程,

    (14)

    該方程中m為未知參量,需要通過實驗儀器進行數據測量加以確定.

    選取“阻尼2”檔位,打開電機,設定在某一轉速,當觀察到擺輪的受迫振動周期與電機的周期趨于一致時,即受迫振動達到穩(wěn)定時,開始測量.讀出每次擺輪受迫振動的幅值θ和驅動力矩每10次振動周期10T,利用頻閃法測量擺輪受迫振動的位移與驅動力之間的相位差φ,并記錄此時電動機的轉速值,結果如表3所示.

    表3 幅頻特性與相頻特性測量數據

    至此獲得了在一定情況下,受迫振動動力學方程的全部參量:ω0=3.775 rad·s-1,β=0.081,m=1.02,a=0.014,b=-0.002,c=0.005.

    (15)

    把上述一階微分方程組輸入程序,利用數值積分求解其動力學特性,結果如圖3所示.

    圖3 系統(tǒng)時間響應與相平面曲線

    圖3中左側為受迫振動系統(tǒng)的時間響應曲線,右側為系統(tǒng)的相平面曲線.受迫振動系統(tǒng)的時間響應曲線表示了振動振幅隨時間變化的特性,可以看出,隨著振動時間的增加,振動頻率逐漸穩(wěn)定,振幅逐漸增加直到最大值并且保持穩(wěn)定.

    受迫振動系統(tǒng)的相平面曲線圖表示了速度隨振幅變化的關系.可以看出隨著時間的增加,速度與振幅也會趨近一個穩(wěn)定狀態(tài).

    如果改變其中的一個參量,例如將與彈簧非線性效應相關的參數a從0.014調整為10,采用同樣的方法對系統(tǒng)進行數值分析,得到的結果如圖4所示.原本穩(wěn)定的受迫振動系統(tǒng)在振動后期開始出現一定的隨機運動,即混沌現象.

    圖4 系統(tǒng)時間響應與相平面曲線

    基于波耳共振儀的受迫振動,由于其振動方程中涉及的非線性項偏多,不利于定性和定量分析受迫振動系統(tǒng)中的混沌現象,因此可以選擇典型的非線性受迫振動方程——杜芬方程進行探討.

    3 受迫振動方程的混沌

    3.1 杜芬方程

    振動問題中的很多數學模型都可以通過轉化為杜芬方程來進行研究和分析,杜芬方程是研究和分析某些非線性動力學系統(tǒng)的基礎.

    杜芬方程[7]的標準形式為

    (16)

    其中,δ>0,g(x)是含x3項的函數,f(x,t)是周期性函數.

    標準化后的杜芬方程稱之為杜芬系統(tǒng).杜芬系統(tǒng)是一類從簡單受迫振動的物理模型中歸納出的非線性振動系統(tǒng)模型,形式簡潔,具有代表性,其應用范圍十分廣泛.許多動力學過程都與杜芬系統(tǒng)的模型極其相似,比如化學鍵的斷裂、建筑結構的震顫和輪船的橫搖運動、電路周期振蕩、電路信號檢測、電路系統(tǒng)模擬分析和電路的故障檢測等等.

    3.2 硬彈簧型杜芬方程

    杜芬方程一般可以分成四種基本類型,可以選取其中的硬彈簧型杜芬方程[8]進行探討,

    (17)

    將式(17)轉化成一階微分方程組:

    (18)

    其中ω、f、a、β、c為可變參量.

    給定不同的初始條件,利用數值計算的方法來分析硬彈簧型杜芬系統(tǒng)的受迫振動.

    當振動時間為300 s時,結果如圖5所示,從上至下依次是不同f取值下的系統(tǒng)時間響應圖、系統(tǒng)相平面軌跡圖、系統(tǒng)功率密度譜.

    圖5 系統(tǒng)時間響應、相平面軌跡、功率密度譜隨f值的演化

    由圖5可以看出,隨著系統(tǒng)參量f的變化,系統(tǒng)逐漸由周期性運動進入混沌狀態(tài).

    3.3 杜芬方程的混沌

    從圖5的相平面軌跡圖可以看出,系統(tǒng)從只有一個相對穩(wěn)定的平衡位置(原點),逐漸變?yōu)榫哂袃蓚€相對平衡的位置(焦點),這體現了無序混沌狀態(tài)中存在一定程度上的有序特性.相平面軌跡圖中相對平衡位置發(fā)生改變的同時,相軌跡也從圓形逐漸轉化為類似于莫比烏斯環(huán)的橢圓形狀.進一步推測,系統(tǒng)會擁有兩個相對平衡的焦點,相軌跡圍繞這兩個焦點運動,而且可將這兩個焦點類比為三維非線性系統(tǒng)中奇異吸引子[9],如圖6所示.

    圖6 奇異吸引子

    奇異吸引子是混沌系統(tǒng)中獨有的一種吸引子,也被稱之為混沌吸引子,它代表著混沌狀態(tài)中的一種定態(tài).具體表現在系統(tǒng)從任一初始狀態(tài)開始,最終都會被奇異吸引子吸引到相空間中的某一特定區(qū)域.奇異吸引子與一般吸引子有所區(qū)別,當混沌狀態(tài)的軌跡線進入奇異吸引子之后,軌跡線并不會像一般吸引子一樣,吸引軌跡線繼續(xù)圍繞其運動,最終形成一個相對封閉的圖形,而是會讓兩條接近的軌跡線發(fā)生指數分離.從系統(tǒng)外部看,奇異吸引子是在聚集軌跡線;從系統(tǒng)內部看,是軌跡線發(fā)散的過程.奇異吸引子常常成對出現,兩個奇異吸引子同時對軌跡線的作用,共同構筑了混沌系統(tǒng)一定程度上的有序性.

    圖7 從左至右分別為f=1, f=3, f=5時的相空間軌跡線

    從圖7可以看出,杜芬方程中的兩個焦點在含有時間項的坐標下,展開成了兩條軸線,體現了其時間平移不變性,參量f的變化引起相空間軌跡明顯改變,表明混沌系統(tǒng)對參量具有很強的敏感性,但在速度-振幅的相平面軌跡圖中,兩個焦點是穩(wěn)定的,不隨時間發(fā)生變化,這與奇異吸引子的性質非常類似.

    4 結論

    基于波耳共振儀受迫振動實驗,探討其中的非線性影響因素,對阻尼振動方程進行修正,考慮軸承摩擦和彈簧非線性效應,結合Matlab數值計算結果與波耳共振儀實測數據,確定波耳共振儀系統(tǒng)的非線性項參數,分析其非線性效應.通過引入非線性理論的思想和數值分析的計算方法,有助于學生準確認知實驗現象,加深對實驗原理的理解.

    為了豐富實驗教學內容,引入彈簧型杜芬方程,分析其非線性特性,利用數值計算方法探討受迫振動模型進入混沌狀態(tài)的條件及其運動特性.受迫振動是否會進入混沌狀態(tài),可以根據系統(tǒng)響應時間曲線(振幅-時間)、系統(tǒng)相平面曲線(速度-振幅)和功率密度譜曲線(功率-頻率)綜合判斷.如果受迫振動做周期運動,其響應時間曲線會呈現規(guī)則的周期性,其相軌跡圖具有極限環(huán)的特性,其功率圖譜峰值呈現離散的狀態(tài);如果系統(tǒng)進入混沌狀態(tài),其時間曲線會具有無規(guī)則的、類似于噪聲的特性,其相平面曲線在有限區(qū)域內是無周期性的,隨著振動時間的延長,其軌跡線會逐漸鋪滿整個區(qū)域,功率圖譜含有多個峰值;如果系統(tǒng)的響應時間曲線和相平面曲線在上述兩種情況之間,則功率圖譜呈現離散狀態(tài),系統(tǒng)做的是擬周期運動,如果改變它的系統(tǒng)參量或者初始條件,則有可能轉變?yōu)橹芷谶\動或混沌運動.一般情況下,受迫振動受非線性效應影響越大,最終進入混沌狀態(tài)的可能性越高.

    猜你喜歡
    擺輪受迫振動振幅
    基于模糊PID的爬樓輪椅后腿擺輪機構同步控制策略*
    擺輪分揀機技術發(fā)展與分揀控制技術應用研究
    論手表的快擺與慢擺
    科學與技術(2019年5期)2019-04-23 06:49:52
    機械手表擺輪游絲的設計計算
    利用自制教具探究受迫振動與共振
    十大漲跌幅、換手、振幅、資金流向
    十大漲跌幅、換手、振幅、資金流向
    滬市十大振幅
    永磁行星齒輪傳動系統(tǒng)受迫振動分析
    有限差分法模擬電梯懸掛系統(tǒng)橫向受迫振動
    毛片一级片免费看久久久久| 日日摸夜夜添夜夜爱| 久久精品国产自在天天线| 日本黄大片高清| 国产男女超爽视频在线观看| 欧美不卡视频在线免费观看| 亚洲性久久影院| 国产单亲对白刺激| 亚洲av电影不卡..在线观看| 欧美极品一区二区三区四区| 国产精品久久久久久久久免| 国产精品久久久久久久久免| 麻豆av噜噜一区二区三区| 国产黄色视频一区二区在线观看| 2021天堂中文幕一二区在线观| 午夜免费激情av| 少妇裸体淫交视频免费看高清| 亚洲av免费在线观看| 久久热精品热| 成人亚洲欧美一区二区av| 亚洲在线观看片| 婷婷色av中文字幕| 51国产日韩欧美| 中文欧美无线码| 大话2 男鬼变身卡| 肉色欧美久久久久久久蜜桃 | av国产免费在线观看| 男人舔奶头视频| 老司机影院毛片| 内射极品少妇av片p| 观看免费一级毛片| 中文字幕制服av| 国产精品国产三级国产专区5o| 日本一二三区视频观看| 国产精品福利在线免费观看| 三级国产精品欧美在线观看| 97超碰精品成人国产| videos熟女内射| 男人和女人高潮做爰伦理| 欧美另类一区| 男女下面进入的视频免费午夜| av免费观看日本| 亚洲av日韩在线播放| 亚洲色图av天堂| 亚洲在线自拍视频| 久99久视频精品免费| 偷拍熟女少妇极品色| 69人妻影院| 99久久精品热视频| 一边亲一边摸免费视频| 国产高清国产精品国产三级 | 亚洲精品亚洲一区二区| 国产毛片a区久久久久| av福利片在线观看| 搞女人的毛片| 可以在线观看毛片的网站| 国产大屁股一区二区在线视频| 91精品国产九色| 青春草国产在线视频| 春色校园在线视频观看| 天天躁日日操中文字幕| 少妇人妻一区二区三区视频| 欧美精品一区二区大全| 99视频精品全部免费 在线| 伊人久久精品亚洲午夜| 国产精品99久久久久久久久| 可以在线观看毛片的网站| 十八禁国产超污无遮挡网站| 亚洲av不卡在线观看| 最近的中文字幕免费完整| 国产精品久久久久久av不卡| 一级毛片久久久久久久久女| 国产视频首页在线观看| 亚洲av中文av极速乱| 我的老师免费观看完整版| 女人久久www免费人成看片| 日日摸夜夜添夜夜添av毛片| 黄色欧美视频在线观看| 亚洲人成网站高清观看| 久久久久久久国产电影| 成年人午夜在线观看视频 | av一本久久久久| 激情 狠狠 欧美| 3wmmmm亚洲av在线观看| 在线免费观看的www视频| 精品一区在线观看国产| 国产成人免费观看mmmm| 免费播放大片免费观看视频在线观看| 精品一区二区三区人妻视频| .国产精品久久| 国产午夜精品论理片| 国产在线男女| 成年女人在线观看亚洲视频 | 国产综合精华液| 99九九线精品视频在线观看视频| 国产高清国产精品国产三级 | 内地一区二区视频在线| 国产乱来视频区| 日韩 亚洲 欧美在线| 亚洲av电影不卡..在线观看| 天堂av国产一区二区熟女人妻| 精品久久久噜噜| 免费大片黄手机在线观看| 国产精品一区二区三区四区免费观看| 美女脱内裤让男人舔精品视频| 22中文网久久字幕| 一区二区三区高清视频在线| 国产女主播在线喷水免费视频网站 | 99久国产av精品国产电影| 久久精品国产亚洲网站| 欧美xxxx性猛交bbbb| 国产片特级美女逼逼视频| 午夜福利视频精品| 超碰av人人做人人爽久久| 精品一区在线观看国产| 成人漫画全彩无遮挡| 国产黄片视频在线免费观看| 91精品伊人久久大香线蕉| 日韩欧美精品免费久久| 国产麻豆成人av免费视频| 午夜福利视频精品| 久久精品久久精品一区二区三区| 国产精品久久久久久av不卡| 激情五月婷婷亚洲| 亚洲最大成人av| 搡老乐熟女国产| 久久久久久久午夜电影| 两个人视频免费观看高清| 精品久久久精品久久久| 建设人人有责人人尽责人人享有的 | 日本一本二区三区精品| 久久久久久久亚洲中文字幕| 伊人久久国产一区二区| 久热久热在线精品观看| 蜜桃久久精品国产亚洲av| 嘟嘟电影网在线观看| 欧美3d第一页| 成年免费大片在线观看| 国产单亲对白刺激| 乱系列少妇在线播放| a级毛色黄片| 秋霞在线观看毛片| 国产视频内射| 久久99热6这里只有精品| 亚洲一区高清亚洲精品| 国产黄色视频一区二区在线观看| 国产亚洲av嫩草精品影院| 久久久色成人| 午夜激情福利司机影院| 能在线免费观看的黄片| 男女视频在线观看网站免费| 亚洲国产精品成人综合色| 六月丁香七月| 在线观看免费高清a一片| 午夜老司机福利剧场| 97热精品久久久久久| 精品久久久精品久久久| 一个人观看的视频www高清免费观看| 午夜福利视频精品| 91在线精品国自产拍蜜月| 少妇猛男粗大的猛烈进出视频 | 欧美区成人在线视频| 久久久久久久久久久免费av| 亚洲经典国产精华液单| 国产av码专区亚洲av| 天天躁夜夜躁狠狠久久av| 只有这里有精品99| 亚洲伊人久久精品综合| 亚洲综合色惰| 国产黄a三级三级三级人| 免费不卡的大黄色大毛片视频在线观看 | 精品99又大又爽又粗少妇毛片| 天堂俺去俺来也www色官网 | 久久久久久国产a免费观看| 国产伦精品一区二区三区视频9| 日韩不卡一区二区三区视频在线| 看非洲黑人一级黄片| 免费电影在线观看免费观看| 亚洲18禁久久av| 激情五月婷婷亚洲| 性插视频无遮挡在线免费观看| 亚洲aⅴ乱码一区二区在线播放| 青青草视频在线视频观看| 久久久久久久久久久丰满| 欧美日韩视频高清一区二区三区二| 麻豆成人午夜福利视频| 免费黄网站久久成人精品| 久久久a久久爽久久v久久| 熟女电影av网| 一级毛片 在线播放| av卡一久久| av线在线观看网站| av专区在线播放| 国产一级毛片七仙女欲春2| 亚洲综合色惰| 欧美成人一区二区免费高清观看| 欧美极品一区二区三区四区| 亚洲性久久影院| 精品一区二区三区人妻视频| 热99在线观看视频| 免费看美女性在线毛片视频| 极品教师在线视频| 日本熟妇午夜| 亚洲欧美精品专区久久| 日韩av免费高清视频| 插阴视频在线观看视频| 日本三级黄在线观看| 国产精品99久久久久久久久| 日韩人妻高清精品专区| 国产男人的电影天堂91| 在线观看美女被高潮喷水网站| 能在线免费看毛片的网站| 亚洲精品乱码久久久v下载方式| 国产黄片美女视频| 国产精品.久久久| av网站免费在线观看视频 | 亚洲精品国产成人久久av| 少妇人妻精品综合一区二区| 中文字幕久久专区| 国产成人精品久久久久久| 欧美精品国产亚洲| 好男人在线观看高清免费视频| 免费黄频网站在线观看国产| 久久久久国产网址| 中文天堂在线官网| 亚洲av男天堂| 亚洲精品国产av成人精品| 成人鲁丝片一二三区免费| 亚洲欧美成人精品一区二区| 永久免费av网站大全| 一级片'在线观看视频| 国产亚洲精品久久久com| 国产一级毛片在线| 色5月婷婷丁香| 日韩人妻高清精品专区| 男女国产视频网站| 一本一本综合久久| 欧美97在线视频| 日产精品乱码卡一卡2卡三| 亚洲av电影在线观看一区二区三区 | 亚洲内射少妇av| 亚洲av免费高清在线观看| 极品少妇高潮喷水抽搐| 久久亚洲国产成人精品v| 日本免费在线观看一区| 国内揄拍国产精品人妻在线| 在线观看免费高清a一片| 亚洲久久久久久中文字幕| 91狼人影院| 日韩一本色道免费dvd| 日韩制服骚丝袜av| 欧美成人a在线观看| 熟妇人妻不卡中文字幕| 一个人看的www免费观看视频| 亚洲综合色惰| 好男人视频免费观看在线| av线在线观看网站| 三级经典国产精品| av免费观看日本| 在线观看美女被高潮喷水网站| 嘟嘟电影网在线观看| 欧美最新免费一区二区三区| 免费播放大片免费观看视频在线观看| 18禁裸乳无遮挡免费网站照片| 亚洲久久久久久中文字幕| 熟妇人妻久久中文字幕3abv| 精品一区二区三区人妻视频| 国产国拍精品亚洲av在线观看| 五月天丁香电影| 高清在线视频一区二区三区| 国产激情偷乱视频一区二区| 麻豆久久精品国产亚洲av| 有码 亚洲区| 99久国产av精品国产电影| 成人美女网站在线观看视频| 边亲边吃奶的免费视频| 亚洲av成人精品一二三区| 免费大片18禁| 成人二区视频| 久久鲁丝午夜福利片| 日韩欧美精品v在线| 欧美性猛交╳xxx乱大交人| 韩国av在线不卡| 天堂av国产一区二区熟女人妻| 汤姆久久久久久久影院中文字幕 | 成年人午夜在线观看视频 | 身体一侧抽搐| 高清日韩中文字幕在线| 99久国产av精品国产电影| 成年女人在线观看亚洲视频 | 久久久久久久久大av| 午夜福利视频1000在线观看| 日韩视频在线欧美| 身体一侧抽搐| 亚洲怡红院男人天堂| 好男人视频免费观看在线| 蜜桃久久精品国产亚洲av| 少妇高潮的动态图| av播播在线观看一区| 欧美日韩视频高清一区二区三区二| 国产精品蜜桃在线观看| 精品国产一区二区三区久久久樱花 | 青春草亚洲视频在线观看| 亚洲人成网站在线播| 国产午夜福利久久久久久| 激情五月婷婷亚洲| 久热久热在线精品观看| 日本免费在线观看一区| 国产成人福利小说| 免费在线观看成人毛片| 欧美高清性xxxxhd video| 欧美不卡视频在线免费观看| 久久99热这里只有精品18| 美女xxoo啪啪120秒动态图| 特大巨黑吊av在线直播| 女人久久www免费人成看片| 神马国产精品三级电影在线观看| 国产黄a三级三级三级人| 亚洲成人久久爱视频| 午夜视频国产福利| av在线蜜桃| 精品久久久久久久久av| 国产乱来视频区| 国产精品.久久久| 亚洲伊人久久精品综合| 欧美成人一区二区免费高清观看| 日韩欧美国产在线观看| 纵有疾风起免费观看全集完整版 | 久久6这里有精品| 久久久久久久国产电影| 国产av码专区亚洲av| 精品午夜福利在线看| 国产麻豆成人av免费视频| 欧美潮喷喷水| 国产黄色免费在线视频| 一个人看视频在线观看www免费| 国产高潮美女av| 亚州av有码| 能在线免费观看的黄片| 99re6热这里在线精品视频| 国产 一区精品| 成人午夜高清在线视频| 精品久久久久久久久久久久久| 日本-黄色视频高清免费观看| 国产成年人精品一区二区| 久久久久精品性色| 国产又色又爽无遮挡免| 久久久久精品久久久久真实原创| 69人妻影院| 成人综合一区亚洲| 久久99精品国语久久久| av卡一久久| 国产大屁股一区二区在线视频| 午夜日本视频在线| 大片免费播放器 马上看| 精品久久久久久久人妻蜜臀av| 国产一区二区在线观看日韩| 国产片特级美女逼逼视频| 人妻制服诱惑在线中文字幕| 高清av免费在线| 成人亚洲精品av一区二区| 欧美三级亚洲精品| 一级毛片电影观看| 97精品久久久久久久久久精品| 伦理电影大哥的女人| 国产精品一区二区三区四区久久| 乱系列少妇在线播放| 少妇熟女aⅴ在线视频| 亚洲内射少妇av| 啦啦啦啦在线视频资源| 日本免费在线观看一区| 国产亚洲午夜精品一区二区久久 | 日韩欧美精品免费久久| 亚洲丝袜综合中文字幕| 91午夜精品亚洲一区二区三区| 免费大片18禁| 直男gayav资源| av网站免费在线观看视频 | 中文乱码字字幕精品一区二区三区 | 国产男女超爽视频在线观看| 高清av免费在线| 天堂俺去俺来也www色官网 | 中文欧美无线码| 精品久久久久久久末码| 免费人成在线观看视频色| 国产黄频视频在线观看| 一区二区三区免费毛片| 亚洲性久久影院| 777米奇影视久久| 国产单亲对白刺激| 我的老师免费观看完整版| 丝袜美腿在线中文| 青春草视频在线免费观看| 国产一区二区三区av在线| 久久久久性生活片| 免费观看a级毛片全部| 国产一级毛片在线| 欧美日韩在线观看h| 又爽又黄无遮挡网站| 亚洲欧美精品专区久久| 亚洲成人一二三区av| 亚洲成人中文字幕在线播放| 亚洲精品国产av蜜桃| av免费观看日本| 亚洲精品成人av观看孕妇| 国产av码专区亚洲av| 亚洲精品国产成人久久av| av黄色大香蕉| 国产精品无大码| 久久久久精品性色| 国产老妇女一区| 国产欧美另类精品又又久久亚洲欧美| 国产精品国产三级国产专区5o| 三级经典国产精品| 精品一区在线观看国产| 99热这里只有是精品50| 亚洲熟女精品中文字幕| 精品人妻熟女av久视频| 亚洲经典国产精华液单| 校园人妻丝袜中文字幕| 成人综合一区亚洲| 免费av毛片视频| 亚洲精品中文字幕在线视频 | 亚洲怡红院男人天堂| 在线观看人妻少妇| 国产成人福利小说| 欧美人与善性xxx| 亚洲aⅴ乱码一区二区在线播放| 熟妇人妻久久中文字幕3abv| 乱人视频在线观看| 99久久精品一区二区三区| 亚洲国产精品sss在线观看| 国产麻豆成人av免费视频| 两个人的视频大全免费| 精品久久久精品久久久| 成人高潮视频无遮挡免费网站| 欧美精品国产亚洲| 精品久久久久久成人av| 国产永久视频网站| 欧美日韩国产mv在线观看视频 | 在线免费观看的www视频| a级毛色黄片| 男的添女的下面高潮视频| 成人无遮挡网站| 久久久久久久久久成人| 一级a做视频免费观看| 99热这里只有是精品50| 久久这里只有精品中国| 免费黄频网站在线观看国产| 久久久精品94久久精品| 免费大片黄手机在线观看| 午夜老司机福利剧场| 日韩欧美 国产精品| 国产乱来视频区| 黄片wwwwww| 亚洲美女搞黄在线观看| 天堂√8在线中文| 91精品一卡2卡3卡4卡| 2018国产大陆天天弄谢| 少妇熟女欧美另类| 边亲边吃奶的免费视频| 性色avwww在线观看| 黄色配什么色好看| 免费大片18禁| 日韩在线高清观看一区二区三区| 床上黄色一级片| 成人综合一区亚洲| 国产精品一区二区三区四区久久| 免费看不卡的av| 成人午夜高清在线视频| 午夜福利在线观看免费完整高清在| 乱系列少妇在线播放| 精品人妻视频免费看| 别揉我奶头 嗯啊视频| 亚洲av成人av| 国产精品99久久久久久久久| 一级毛片我不卡| 在线观看av片永久免费下载| 久热久热在线精品观看| 免费不卡的大黄色大毛片视频在线观看 | 汤姆久久久久久久影院中文字幕 | 免费看a级黄色片| 熟妇人妻不卡中文字幕| 2018国产大陆天天弄谢| 男人狂女人下面高潮的视频| 亚洲一级一片aⅴ在线观看| 日韩一区二区视频免费看| 久久综合国产亚洲精品| 亚洲成人av在线免费| 亚洲欧美中文字幕日韩二区| 国产一区二区三区av在线| 如何舔出高潮| 国产精品一区二区三区四区久久| 高清视频免费观看一区二区 | 91午夜精品亚洲一区二区三区| 久热久热在线精品观看| av又黄又爽大尺度在线免费看| 亚洲国产日韩欧美精品在线观看| 精品久久久久久久末码| 男女啪啪激烈高潮av片| 亚洲人成网站在线播| 国产精品国产三级国产av玫瑰| 亚洲丝袜综合中文字幕| 成人午夜高清在线视频| 欧美区成人在线视频| 高清欧美精品videossex| 精品熟女少妇av免费看| 久久6这里有精品| 看免费成人av毛片| 免费观看性生交大片5| 91精品一卡2卡3卡4卡| 777米奇影视久久| 国产高清三级在线| 热99在线观看视频| 三级国产精品片| 纵有疾风起免费观看全集完整版 | 亚洲图色成人| 日韩精品有码人妻一区| 别揉我奶头 嗯啊视频| 久久久成人免费电影| 久久久久久久久久久丰满| 亚洲国产精品专区欧美| 国产精品久久久久久久久免| 亚洲欧美精品专区久久| 国产精品伦人一区二区| 日韩伦理黄色片| 欧美成人午夜免费资源| 久久久久久久久久黄片| 内射极品少妇av片p| 国产亚洲精品久久久com| av播播在线观看一区| 观看免费一级毛片| 国产精品精品国产色婷婷| 成年版毛片免费区| 久久午夜福利片| av.在线天堂| 午夜精品在线福利| 全区人妻精品视频| 国产在视频线精品| 色网站视频免费| 国模一区二区三区四区视频| 色网站视频免费| 午夜亚洲福利在线播放| 久久久久久久久久成人| 日韩成人av中文字幕在线观看| 可以在线观看毛片的网站| 床上黄色一级片| 国产免费又黄又爽又色| 国产国拍精品亚洲av在线观看| 久久久久久久久中文| 视频中文字幕在线观看| 色综合色国产| 一本久久精品| 男的添女的下面高潮视频| 国产精品国产三级国产专区5o| 97人妻精品一区二区三区麻豆| 中文字幕亚洲精品专区| 亚洲欧美成人综合另类久久久| videos熟女内射| 国产男人的电影天堂91| 免费大片黄手机在线观看| 成人性生交大片免费视频hd| 婷婷六月久久综合丁香| 91精品一卡2卡3卡4卡| 亚洲欧美清纯卡通| 国产一区二区三区av在线| 日韩视频在线欧美| 在线观看免费高清a一片| 日韩精品青青久久久久久| 亚洲天堂国产精品一区在线| 一个人免费在线观看电影| 午夜福利网站1000一区二区三区| 欧美成人精品欧美一级黄| 国产精品三级大全| 成人亚洲欧美一区二区av| 麻豆av噜噜一区二区三区| 国产精品一区二区三区四区久久| 国产精品一区二区三区四区免费观看| 最近的中文字幕免费完整| 亚洲国产av新网站| kizo精华| 国内精品宾馆在线| 欧美日韩精品成人综合77777| 午夜免费观看性视频| 成人无遮挡网站| 深夜a级毛片| 美女被艹到高潮喷水动态| 性插视频无遮挡在线免费观看| av网站免费在线观看视频 | 中国美白少妇内射xxxbb| 男女视频在线观看网站免费| 日本黄大片高清| 99热这里只有是精品在线观看| 国产成人免费观看mmmm| 不卡视频在线观看欧美| 亚洲国产色片| 亚洲国产精品专区欧美| 亚洲在久久综合| 熟女人妻精品中文字幕| 99热网站在线观看| 日韩欧美一区视频在线观看 | 两个人视频免费观看高清| 色综合色国产| 午夜免费激情av| 神马国产精品三级电影在线观看| 麻豆av噜噜一区二区三区| 亚洲成人中文字幕在线播放| 日本三级黄在线观看| 亚洲精品色激情综合| 国产激情偷乱视频一区二区| 免费看美女性在线毛片视频| 99re6热这里在线精品视频| 亚洲aⅴ乱码一区二区在线播放| 在线a可以看的网站| 日韩欧美精品免费久久|