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

    基于廣義預測控制的智能巖石試驗機研究

    2020-12-04 07:03:16高繼開姚志賓何本國
    液壓與氣動 2020年11期
    關(guān)鍵詞:墊塊試驗機控制算法

    高繼開,田 軍,姚志賓,何本國

    (東北大學深部金屬礦山安全開采教育部重點實驗室,遼寧沈陽 110819)

    引言

    巖石力學強度指標的測定是巖石力學與工程研究的基礎(chǔ),巖石力學試驗機是測量巖石力學性質(zhì)最常見的測試儀器。從系統(tǒng)原理來看,巖石壓力、拉力試驗機是典型的液壓系統(tǒng),試驗機已經(jīng)從原來的人工手動調(diào)節(jié)閥門流量控制油缸動作、輸出壓力,普遍發(fā)展為電液伺服數(shù)字控制系統(tǒng)[1]。國外數(shù)字控制巖石試驗機處于領(lǐng)先地位,以美國MTS公司、英國GCTS公司、法國TOP公司為代表占據(jù)了大量市場,近年來,國內(nèi)試驗機技術(shù)進展迅速,東北大學深部金屬礦山安全開采教育部重點實驗室聯(lián)合長春朝陽、科意、中機試驗等廠商研發(fā)了具有自主知識產(chǎn)權(quán)的國際領(lǐng)先系列真三軸設(shè)備[2-3]。目前國內(nèi)對控制精度要求高的試驗機大多采用德國DOLI公司的EDC控制器或美國DELTA公司的RMC控制器,控制器廠商設(shè)計的PID算法內(nèi)嵌在控制器中。國外控制器性能優(yōu)異,但價格昂貴,且控制技術(shù)封閉。國內(nèi)試驗機控制器需要軟件和硬件的雙重突破。

    試驗機控制器的PID參數(shù)通常由技術(shù)人員在出廠時設(shè)置好就不再更改,然而巖石力學實驗是一個復雜的過程,不同種類巖石強度差異甚遠,同一巖石試件在加載過程的不同階段也呈現(xiàn)出不同的性質(zhì)[4],另外還會出現(xiàn)設(shè)備維修老化等造成的實驗環(huán)境變化。巖石變形控制對控制精度要求極高,需要變換PID參數(shù)以滿足多變的應(yīng)用環(huán)境,然而試驗機PID參數(shù)調(diào)節(jié)過程復雜,即使有經(jīng)驗的技術(shù)人員也需要反復調(diào)節(jié)很長時間。當實驗環(huán)境發(fā)生變化,舊的參數(shù)不滿足控制需求時,只能耗費大量時間和財力通過儀器生產(chǎn)廠家再次調(diào)節(jié)。使用一種具有自適應(yīng)功能的控制算法及時調(diào)節(jié)控制策略,滿足不同條件下的控制需求是必要的。GPC廣義預測控制是在自適應(yīng)控制方法的研究中發(fā)展起來的一類預測控制算法。與傳統(tǒng)自校正控制算法相比,GPC采用實時辨識的建模策略,對控制對象的數(shù)學模型精度要求更低[5],GPC已經(jīng)在多種復雜工業(yè)過程獲得了應(yīng)用[6-7]。GPC實時辨識巖石試驗機加載過程模型參數(shù),根據(jù)模型預測值和未來期望值求算控制量,可以有效規(guī)避繁雜的PID參數(shù)調(diào)節(jié)過程并智能適應(yīng)實驗條件變化。

    巖石壓力試驗機是一個液壓控制系統(tǒng),包含液壓油源、控制器和作動器。液壓油源輸出動力油源,伺服閥通過流量的方向和大小調(diào)節(jié)作動器的位移速度,作動器通過壓頭加載在巖樣上,使巖樣產(chǎn)生變形。巖石試驗機不同于金屬萬能試驗機[8],其噸位大、試樣剛度變化大、控制精度要求更高,根據(jù)不同科研需要,加載過程復雜多變,尤其是測定巖石峰后承載能力時,是對非線性時變系統(tǒng)的高精度控制。

    圖1 巖石試驗機控制系統(tǒng)框圖

    圖1為試驗機伺服控制系統(tǒng)的框圖??刂破鹘邮掌谕冃沃祔sp和變形傳感器反饋的實際值y的誤差e,經(jīng)控制算法計算后輸出控制電壓u給伺服閥,進行方向和流量調(diào)節(jié)。

    圖2為巖樣加載初始狀態(tài)圖,變形傳感器為LVDT差動位移傳感器,作動器在試驗機最底部,加載時帶動底部墊塊、巖樣和球頭向上移動,球頭與上壓頭接觸,最終使巖樣產(chǎn)生變形。巖樣加載過程存在多個非線性環(huán)節(jié),一方面是由于試驗機本身造成的,如由電液伺服閥和系統(tǒng)摩擦產(chǎn)生的非線性環(huán)節(jié)[9-10],作動器與壓頭剛接觸時的受力狀態(tài)變化等;另一方面是由于巖樣本身性質(zhì)造成的,同種巖石從施加應(yīng)力到破壞完畢大多會呈現(xiàn)出彈性階段、塑性階段、脆性階段以及峰后階段等多種性質(zhì)變化[11],存在嚴重的時變特性,不僅如此,不同巖石加載過程中所呈現(xiàn)的力學性質(zhì)也存在巨大差異,且實驗過程常要求巖石恒速率變形,對控制精度要求極高,這都為試驗機控制算法帶來了巨大的挑戰(zhàn)。

    圖2 巖樣加載初始狀態(tài)圖

    1 巖石試驗機變形控制模型辨識

    系統(tǒng)存在諸多不確定非線性環(huán)節(jié),為簡化問題,將系統(tǒng)辨識為空載和加載時的2個線性模型,辨識方法采用最小二乘法。

    考慮式(1)單輸入輸出系統(tǒng):

    A(z-1)y(k)=B(z-1)u(t-nb)+e(k)

    式中,na—— 輸出階次

    nb—— 輸入階次

    y(k) —— 系統(tǒng)輸出

    u(k) —— 系統(tǒng)輸入

    e(k) —— 白噪聲輸入信號

    z-1—— 后移算子

    記:

    可得式(2):

    y(k)=hT(k)θ+e(k)

    e(k)=ψT(k)θ+v(k)

    (2)

    最小二乘的辨識的準則函數(shù)為:

    (4)

    式中,μ為遺忘因子,0<μ<1,一般取0.95~0.98之間;μ為1時是標準遞推最小二乘形式;P(k)為正定的協(xié)方差矩陣。

    分別在空載狀態(tài)和加載高剛度墊塊的情況下,使伺服閥控制電壓均勻變化,記錄變形傳感器的讀數(shù)。假定不同系統(tǒng)階次,選取部分輸入輸出數(shù)據(jù)進行參數(shù)辨識和誤差檢驗,另選取部分數(shù)據(jù)對辨識結(jié)果進行驗證,使用均方誤差值來表征模型辨識準確度。

    空載試驗下,辨識均方誤差隨階次變化如圖3所示;墊塊試驗下,辨識均方誤差隨階次變化如圖4所示。

    圖3 系統(tǒng)階次—均方誤差圖(空載)

    從空載實驗和墊塊加載實驗的辨識結(jié)果來看,系統(tǒng)取4階和5階時均方誤差較小,總是小于取更低階次時的誤差,加載試驗系統(tǒng)取5階時誤差比4階時更大??蛰d實驗的辨識誤差和檢驗誤差都小于墊塊加載實驗的辨識誤差,系統(tǒng)取4階時,墊塊加載實驗檢驗誤差約是空載實驗的15倍,說明墊塊加載的過程相對復雜,開環(huán)加載實驗使用的是質(zhì)地均勻的實心鐵塊,非均質(zhì)巖石加載過程模型會更加復雜,非線性和時變特征更強。最終,取均方差較小的階次,認為巖石驗機變形控制系統(tǒng)為4階系統(tǒng)。

    圖4 系統(tǒng)階次—均方誤差圖(加載)

    2 GPC仿真實驗

    2.1 GPC預測控制

    在GPC預測控制策略中,即使被控對象具有非線性或時變特性,依舊認為預測模型可以采用局部線性化和實時辨識的方式表征被控對象。被控對象當前時刻模型可通過式(4)遺忘因子最小二乘辨識。系統(tǒng)當前時刻的辨識結(jié)果形式符合式(1)差分方程形式。

    圖5 GPC預測控制結(jié)構(gòu)框圖

    1=Ej(z-1)A(z-1)Δ+z-jFj(z-1)

    (5)

    用EjΔzj乘以式(1)兩端,并忽略未來時刻噪聲,得到k時刻為起點的未來j步的最合適預測輸出值表達式:

    Fj(z-1)y(k)

    (6)

    式中,Gj(z-1)=Ej(z-1)B(z-1)。結(jié)合式(1)可得:

    =gj,0+gj,1z-1+…

    (7)

    因此,多項式Gj(z-1)中前j項的系數(shù)與對象階躍響應(yīng)前j項的采樣值一一對應(yīng),記做g1,…,gj,對于j=1,2,…,N,式(8)可以寫成矩陣方程:

    (8)

    i=1,2,…,N)

    其中多項式E,F,f可通過遞推求解丟番圖得到,參考文獻[14-15]。式(9)為GPC在k時刻優(yōu)化性能指標。

    式中,E(·) —— {·}的數(shù)學期望

    yr—— 系統(tǒng)的期望輸出

    N1—— 優(yōu)化時域起始時刻

    N2—— 優(yōu)化時域終止時刻

    M—— 控制時域,在M步后控制增量0

    λ(j) —— 控制指標的權(quán)系數(shù)

    式(9)中系統(tǒng)輸出值用預測值式(8)代替,則有當λI+GTG非奇異時,得到使性能指標式(9)最優(yōu)的解為:

    Δu(k)=(λI+GTG)-1GT(yr(k)-f(k)) (10)

    此時最優(yōu)控制量則可由式(11)給出:

    u(k)=u(k-1)+dT(yr(k)-f(k))

    (11)

    式中,dT為矩陣(λI+GTG)-1GT的第一行。

    2.2 仿真實驗

    考慮到GPC控制算法中存在大量矩陣運算,并且模型辨識環(huán)節(jié)實時進行,可以對控制對象局部線性化,因此在控制算法中使用2階模型在線辨識、預測輸出并計算控制量來降低計算量。

    巖石變形實驗多為不同速率下勻速加載過程,因此仿真實驗期望輸出采用斜坡函數(shù)。GPC仿真效果圖如圖6所示,系統(tǒng)初始變形為1 mm,在第2~10秒,10~14秒,14~22秒的各時段之間采用不同速率勻速加載,分別為0.005,0.003,0.001 mm/s,期望變形值在22 s后保持2.5 mm不變。

    圖6 GPC仿真過程曲線

    根據(jù)系統(tǒng)在第1章模型辨識結(jié)果,對2個4階被控對象模型仿真實驗。仿真開始階段,被控對象使用的是空載辨識模型,系統(tǒng)輸出跟隨性很好,變形速率與期望速率吻合。在第6秒,被控對象由空載模型突變?yōu)閴|塊加載模型,此時控制算法模型辨識參數(shù)處于調(diào)整階段,故在6~7秒時,雖然期望速率未發(fā)生變化,但原來已經(jīng)穩(wěn)定的跟蹤誤差發(fā)生了變化,這是控制器適應(yīng)被控對象模型變化的過程。

    對于仿真過程,空載模型和墊塊加載模型使用了完全相同的控制參數(shù),包括預測時域、控制時域、遺忘因子和性能指標均衡參數(shù),體現(xiàn)了控制算法的自適應(yīng)能力。GPC參數(shù)調(diào)節(jié)過程十分簡單,根據(jù)控制效果簡單調(diào)整后即可獲得較為理想的控制效果。從整體控制效果來看,基于GPC的變形控制系統(tǒng)具有很好的動態(tài)性能和穩(wěn)態(tài)性能,說明2階模型不僅當前時刻可以很好地表征被控對象,在預測時域的局部范圍內(nèi)對被控對象的預測輸出偏差也是滿足控制需求的。

    3 實驗結(jié)果

    通過前兩章的系統(tǒng)辨識與仿真,將GPC應(yīng)用到試驗機變形控制的可行性已經(jīng)得到驗證,并且控制結(jié)構(gòu)和相關(guān)參數(shù)已經(jīng)基本確定。使用C Sharp語言編寫控制算法,并開展巖石力學單軸實驗,驗證實際控制效果。軟件編程完畢后,采用高剛度墊塊的加載實驗進行設(shè)備調(diào)試,得到如圖7所示仿真結(jié)果。初始期望加載速率為0.01 mm/s,此時試驗機上壓頭還未與球頭接觸,直到第13.52秒時,墊塊開始受力,實驗力顯著增加,但是此過程參數(shù)變化平滑,且變化幅值較小,在控制量的及時調(diào)整下,系統(tǒng)輸出依舊保持了預定加載速率。隨后分別進行了0.004,0.001 mm/s和-0.001,-0.005 mm/s的加載、保載和卸載試驗。圖7b顯示跟蹤誤差在速率變換時較大,但也在保持在微米級,且控制器調(diào)節(jié)迅速,整體上取得了滿意的控制效果。

    圖7 GPC巖石試驗機實驗過程曲線(鐵塊)

    從圖8辨識模型參數(shù)變化過程來看,辨識模型參數(shù)實時變化,并與加載速率呈現(xiàn)出一定的相關(guān)性,變形速率為0時,模型參數(shù)幅值變化較小,趨于穩(wěn)定。

    圖8 GPC巖石試驗機模型參數(shù)變化曲線(鐵塊)

    設(shè)備調(diào)試完畢后,使用設(shè)備開展巖石力學實驗,測試紅砂巖和花崗巖的單軸抗壓強度。圖9為花崗巖強度測試過程曲線?;◢弾r恒速率變形加載實驗使用0.001 mm/s的變形速率進行加載,并在實驗過程中保載了117 s。在GPC控制策略作用下,巖石變形量與期望值基本重合,控制效果較好。巖石在1113秒時進入峰后破壞階段,強度變小,內(nèi)部開裂嚴重[16],導致在峰后階段出現(xiàn)了2處較大的跟蹤誤差,GPC及時調(diào)整了控制量,能夠在巖石峰后較弱的承載能力下繼續(xù)按期望加載而不使巖石發(fā)生立即完全破壞。

    圖9 GPC巖石試驗機實驗過程曲線(花崗巖)

    GPC辨識的模型參數(shù)變化如圖10所示,花崗巖加載過程中模型參數(shù)實時調(diào)整,在恒加載速率下,參數(shù)變

    圖10 GPC巖石試驗機模型參數(shù)變化曲線(花崗巖)

    化幅值較小,加載速率為0時,參數(shù)趨于穩(wěn)定。最終測得0.001 mm/s變形速率下,花崗巖試件實驗力峰值為114.76 kN,抗壓強度為73.76 MPa,同時得到了巖樣峰后壓力和變形值,取得了十分滿意的實驗效果。

    圖11為砂巖測試的實驗過程曲線,巖石變形按照設(shè)定的0.004 mm/s速率穩(wěn)定加載。由于砂巖樣品內(nèi)部存在裂隙,在58.52秒時引起了壓力值波動,出現(xiàn)了允許范圍內(nèi)的較大跟蹤誤差,控制器快速調(diào)整了控制策略,繼續(xù)按期望加載,完成了砂巖強度測試。圖12為砂巖實驗過程模型參數(shù)變化曲線,模型參數(shù)整體穩(wěn)定,在58.52秒時有小幅度突變。砂巖樣品試驗力峰值為127.25 kN,強度為 64.84 MPa。

    圖11 GPC巖石試驗機實驗過程曲線(砂巖)

    圖12 GPC巖石試驗機模型參數(shù)變化曲線(砂巖)

    圖13為巖石破壞后形態(tài),可以看出樣品中存在較大節(jié)理,出現(xiàn)提前脫落現(xiàn)象,極大降低了巖石強度[17]。2種巖石破碎效果差異明顯,花崗巖內(nèi)部存在更多節(jié)理裂隙,破碎后存在大量塊狀形態(tài),砂巖均質(zhì)性更好,多以粉末狀破碎。雖然巖石性質(zhì)差異巨大,但沒有影響控制效果,驗證了GPC對于不同種類巖石較好的適應(yīng)性。

    圖13 巖石破壞后形態(tài)

    4 結(jié)論

    通過對試驗機電液伺服控制系統(tǒng)的模型辨識和GPC廣義預測控制策略的研究,確認了系統(tǒng)階次,明確了算法可行性。巖石軸向變形控制系統(tǒng)是一個4階系統(tǒng),可以采用局部線性化的方式在GPC中使用2階模型實時辨識,減小控制算法計算量。

    用C Sharp語言編寫的GPC控制策略可以應(yīng)用于巖石試驗機電液伺服系統(tǒng),完成不同速率下的巖石力學恒速率加載實驗。GPC廣義預測控制,參數(shù)調(diào)節(jié)過程簡單,能夠有效規(guī)避PID參數(shù)調(diào)節(jié)的復雜過程,可針對不同巖石種類和巖石加載過程的非線性、時變性質(zhì)實時調(diào)整控制策略,自適應(yīng)和智能化程度高,具有較好的應(yīng)用推廣價值。

    猜你喜歡
    墊塊試驗機控制算法
    鐵路軸承試驗機
    哈爾濱軸承(2020年2期)2020-11-06 09:22:28
    接觸疲勞試驗機伺服電動缸加載系統(tǒng)設(shè)計
    哈爾濱軸承(2020年1期)2020-11-03 09:16:08
    軋鋼加熱爐墊塊的選型應(yīng)用及發(fā)展趨勢
    CL體系在施工過程中混凝土墊塊制作的施工工藝
    基于ARM+FPGA的模塊化同步控制算法研究
    多孔鉭杯及鉭金屬墊塊在Paprosky Ⅲ、Ⅳ型髖臼缺損翻修中的應(yīng)用研究
    橡膠整體式履帶動態(tài)嚙合試驗機的設(shè)計
    一種優(yōu)化的基于ARM Cortex-M3電池組均衡控制算法應(yīng)用
    一種非圓旋轉(zhuǎn)工件支撐裝置控制算法
    水泥電動抗折試驗機的應(yīng)用研究
    河南科技(2014年10期)2014-02-27 14:09:20
    内地一区二区视频在线| 欧美bdsm另类| 伦精品一区二区三区| 久久精品国产亚洲av涩爱| 99久久人妻综合| 最近最新中文字幕大全电影3| 青春草国产在线视频| av国产久精品久网站免费入址| 女人十人毛片免费观看3o分钟| 嫩草影院入口| av在线天堂中文字幕| 搡女人真爽免费视频火全软件| 欧美成人午夜免费资源| 亚洲av一区综合| 在线观看av片永久免费下载| 尾随美女入室| 亚洲经典国产精华液单| 亚洲精品视频女| 精华霜和精华液先用哪个| 最近最新中文字幕免费大全7| 18禁在线播放成人免费| 一本一本综合久久| 蜜桃亚洲精品一区二区三区| 国产免费又黄又爽又色| 久久亚洲国产成人精品v| 欧美日韩综合久久久久久| ponron亚洲| 婷婷六月久久综合丁香| 国产黄a三级三级三级人| 亚洲18禁久久av| xxx大片免费视频| 亚洲国产av新网站| 日韩电影二区| 特级一级黄色大片| 日本免费在线观看一区| 97人妻精品一区二区三区麻豆| 国产国拍精品亚洲av在线观看| 欧美性猛交╳xxx乱大交人| 久久久久久国产a免费观看| 嫩草影院新地址| 午夜免费激情av| 国产淫片久久久久久久久| 午夜爱爱视频在线播放| 女人被狂操c到高潮| 日本三级黄在线观看| 久久6这里有精品| 青春草亚洲视频在线观看| 最近中文字幕高清免费大全6| 97精品久久久久久久久久精品| 亚洲精品乱码久久久v下载方式| 日韩电影二区| 午夜福利在线观看免费完整高清在| 日韩一本色道免费dvd| 草草在线视频免费看| 有码 亚洲区| 亚洲国产最新在线播放| 久久久久性生活片| 免费人成在线观看视频色| 免费av不卡在线播放| av在线蜜桃| 91午夜精品亚洲一区二区三区| 只有这里有精品99| 久久久成人免费电影| 床上黄色一级片| 国产精品日韩av在线免费观看| 日产精品乱码卡一卡2卡三| 国产乱人视频| 日日啪夜夜撸| 少妇人妻一区二区三区视频| eeuss影院久久| 日产精品乱码卡一卡2卡三| 少妇高潮的动态图| 免费av观看视频| 波多野结衣巨乳人妻| 蜜桃久久精品国产亚洲av| 卡戴珊不雅视频在线播放| 免费不卡的大黄色大毛片视频在线观看 | 国产av不卡久久| 亚洲国产精品sss在线观看| av在线观看视频网站免费| 日本黄大片高清| 寂寞人妻少妇视频99o| 国产男女超爽视频在线观看| 久久久久久久久久久免费av| 亚洲色图av天堂| 国产伦理片在线播放av一区| 亚洲av中文字字幕乱码综合| 国产欧美日韩精品一区二区| 午夜福利在线观看免费完整高清在| 日韩中字成人| 少妇丰满av| 久久综合国产亚洲精品| 午夜激情欧美在线| 丝瓜视频免费看黄片| 亚洲精品一二三| 一区二区三区高清视频在线| 精品国产三级普通话版| 国产成人福利小说| 两个人视频免费观看高清| 亚洲国产色片| 99热6这里只有精品| 亚洲av在线观看美女高潮| 久久精品国产亚洲av天美| 在线观看免费高清a一片| 久久久久网色| 韩国av在线不卡| 国产免费又黄又爽又色| 免费av不卡在线播放| 久久国产乱子免费精品| 亚洲一级一片aⅴ在线观看| 视频中文字幕在线观看| 性色avwww在线观看| 成人鲁丝片一二三区免费| 最近最新中文字幕大全电影3| 久久久欧美国产精品| 波野结衣二区三区在线| freevideosex欧美| 97人妻精品一区二区三区麻豆| 色吧在线观看| 女人十人毛片免费观看3o分钟| 亚洲最大成人手机在线| 人妻夜夜爽99麻豆av| 久久精品人妻少妇| 亚洲真实伦在线观看| 国产成人精品婷婷| 日韩欧美 国产精品| 赤兔流量卡办理| 在线a可以看的网站| 午夜福利网站1000一区二区三区| 日韩制服骚丝袜av| 久久综合国产亚洲精品| 成年免费大片在线观看| 成年女人在线观看亚洲视频 | 在现免费观看毛片| 天堂影院成人在线观看| 国产老妇女一区| 亚洲美女视频黄频| 久久6这里有精品| 久久久久性生活片| 国产午夜精品一二区理论片| 人体艺术视频欧美日本| 午夜福利在线在线| 搡女人真爽免费视频火全软件| 2018国产大陆天天弄谢| 一个人观看的视频www高清免费观看| 一本一本综合久久| 欧美日韩亚洲高清精品| 99热这里只有是精品50| 国产亚洲精品av在线| 亚洲色图av天堂| 最近中文字幕2019免费版| 国产黄频视频在线观看| 少妇的逼好多水| 69人妻影院| 在线免费观看不下载黄p国产| 女人被狂操c到高潮| 免费观看精品视频网站| 一个人看的www免费观看视频| 日韩一区二区三区影片| 欧美成人a在线观看| 免费黄色在线免费观看| 久久国内精品自在自线图片| 女人被狂操c到高潮| 少妇高潮的动态图| 久久这里只有精品中国| 国精品久久久久久国模美| 亚洲欧美中文字幕日韩二区| 秋霞伦理黄片| 日韩精品青青久久久久久| 免费不卡的大黄色大毛片视频在线观看 | 精品一区二区免费观看| 九色成人免费人妻av| 99久国产av精品国产电影| 在线观看人妻少妇| 三级经典国产精品| 国内揄拍国产精品人妻在线| av卡一久久| 亚洲图色成人| 亚洲成人av在线免费| 精品久久久久久久久亚洲| 欧美一区二区亚洲| 亚洲精品成人av观看孕妇| 国产在线男女| 中文字幕av成人在线电影| 亚洲av福利一区| 国产亚洲av嫩草精品影院| 亚洲人成网站在线播| 九草在线视频观看| 91午夜精品亚洲一区二区三区| 欧美日韩精品成人综合77777| 国产精品福利在线免费观看| 日韩国内少妇激情av| 久久鲁丝午夜福利片| 女人被狂操c到高潮| 天堂影院成人在线观看| 国产白丝娇喘喷水9色精品| 男的添女的下面高潮视频| 亚洲真实伦在线观看| 中国国产av一级| 亚洲欧美清纯卡通| 成人性生交大片免费视频hd| 少妇人妻一区二区三区视频| 只有这里有精品99| 99热网站在线观看| 精品欧美国产一区二区三| 成人一区二区视频在线观看| 一级毛片黄色毛片免费观看视频| 国内揄拍国产精品人妻在线| av在线播放精品| 国产精品伦人一区二区| 日韩国内少妇激情av| 日本wwww免费看| 性色avwww在线观看| 日本一本二区三区精品| 久久久亚洲精品成人影院| 我要看日韩黄色一级片| 丝袜美腿在线中文| 日本猛色少妇xxxxx猛交久久| 在线天堂最新版资源| 日韩人妻高清精品专区| 亚洲精品国产成人久久av| 久久久久网色| 建设人人有责人人尽责人人享有的 | 一本一本综合久久| 秋霞伦理黄片| 亚洲av.av天堂| 欧美三级亚洲精品| 亚洲怡红院男人天堂| 午夜老司机福利剧场| 亚洲av免费高清在线观看| 免费少妇av软件| 2021天堂中文幕一二区在线观| 国产成人a∨麻豆精品| 夜夜爽夜夜爽视频| 亚洲精品亚洲一区二区| 两个人视频免费观看高清| 国产91av在线免费观看| h日本视频在线播放| 哪个播放器可以免费观看大片| 美女内射精品一级片tv| 午夜久久久久精精品| 国产熟女欧美一区二区| xxx大片免费视频| 午夜精品在线福利| 亚洲国产精品国产精品| 亚洲电影在线观看av| 亚洲av福利一区| 日韩亚洲欧美综合| 午夜免费男女啪啪视频观看| 日日啪夜夜撸| 男女边吃奶边做爰视频| 国产成人aa在线观看| 乱码一卡2卡4卡精品| 中文乱码字字幕精品一区二区三区 | 国产亚洲最大av| 国产亚洲一区二区精品| 久久精品国产亚洲av涩爱| 中文欧美无线码| 国产欧美另类精品又又久久亚洲欧美| 乱系列少妇在线播放| 夜夜爽夜夜爽视频| 国产男女超爽视频在线观看| 欧美日韩在线观看h| 亚洲内射少妇av| 极品教师在线视频| 在线a可以看的网站| 亚洲av成人精品一二三区| 性色avwww在线观看| 国产淫语在线视频| 中文欧美无线码| 亚洲欧美日韩卡通动漫| 女人被狂操c到高潮| 特级一级黄色大片| eeuss影院久久| 五月天丁香电影| 国产综合懂色| 寂寞人妻少妇视频99o| 美女国产视频在线观看| 一本久久精品| 尤物成人国产欧美一区二区三区| 亚洲成色77777| 国产成人91sexporn| 在线观看免费高清a一片| 三级国产精品片| 亚洲最大成人中文| 亚洲真实伦在线观看| 国产精品精品国产色婷婷| 99热6这里只有精品| 18+在线观看网站| 久久久亚洲精品成人影院| 亚洲国产欧美在线一区| av国产久精品久网站免费入址| 国产国拍精品亚洲av在线观看| 国产精品综合久久久久久久免费| 国产国拍精品亚洲av在线观看| 亚洲成人久久爱视频| 精品一区二区三卡| 亚洲欧美精品专区久久| 国语对白做爰xxxⅹ性视频网站| 青春草视频在线免费观看| 我的老师免费观看完整版| 色综合站精品国产| 老司机影院成人| 亚洲无线观看免费| 国产亚洲av片在线观看秒播厂 | 亚洲不卡免费看| 亚洲性久久影院| 中文字幕av成人在线电影| or卡值多少钱| 日韩一区二区三区影片| 人人妻人人看人人澡| 男女啪啪激烈高潮av片| 一级毛片久久久久久久久女| 国产91av在线免费观看| 国内少妇人妻偷人精品xxx网站| 国产黄片视频在线免费观看| 欧美一级a爱片免费观看看| 亚洲av中文av极速乱| 天美传媒精品一区二区| 亚洲av在线观看美女高潮| 在线观看av片永久免费下载| 三级男女做爰猛烈吃奶摸视频| 汤姆久久久久久久影院中文字幕 | 国产精品美女特级片免费视频播放器| 国产精品.久久久| 久久99热这里只有精品18| 在线 av 中文字幕| 免费大片黄手机在线观看| 又粗又硬又长又爽又黄的视频| 在线免费观看不下载黄p国产| 国产精品.久久久| 七月丁香在线播放| 久久久久网色| 日韩一区二区视频免费看| 亚洲精品乱码久久久久久按摩| 亚洲av免费在线观看| 在线观看免费高清a一片| 国产视频内射| 麻豆成人av视频| www.色视频.com| 亚洲国产精品成人久久小说| av又黄又爽大尺度在线免费看| 亚洲18禁久久av| 青春草国产在线视频| 亚洲色图av天堂| 最近最新中文字幕免费大全7| 最后的刺客免费高清国语| 亚洲av成人精品一区久久| 国产女主播在线喷水免费视频网站 | 七月丁香在线播放| 一本一本综合久久| 欧美人与善性xxx| 免费电影在线观看免费观看| 亚洲国产精品成人久久小说| 日韩在线高清观看一区二区三区| 亚洲精品成人久久久久久| 亚洲高清免费不卡视频| 人妻一区二区av| 亚洲av福利一区| 内地一区二区视频在线| 国产亚洲5aaaaa淫片| 性色avwww在线观看| 国产午夜精品论理片| 美女黄网站色视频| 麻豆乱淫一区二区| 国产成人一区二区在线| 美女黄网站色视频| 国产精品伦人一区二区| 国产一级毛片七仙女欲春2| 久久精品久久精品一区二区三区| 国产成人精品一,二区| 日韩欧美精品v在线| 在现免费观看毛片| 少妇高潮的动态图| 有码 亚洲区| 九九爱精品视频在线观看| 五月伊人婷婷丁香| 99久久人妻综合| 搡女人真爽免费视频火全软件| 国产91av在线免费观看| 亚洲欧洲国产日韩| 亚洲精品乱码久久久v下载方式| 在线观看一区二区三区| 尤物成人国产欧美一区二区三区| 亚州av有码| 久久久a久久爽久久v久久| av免费在线看不卡| 免费观看精品视频网站| 国内精品宾馆在线| 亚洲精华国产精华液的使用体验| 国产黄频视频在线观看| 免费观看a级毛片全部| 美女内射精品一级片tv| 国产伦一二天堂av在线观看| 国产一区二区亚洲精品在线观看| 亚洲精品日本国产第一区| 国内精品一区二区在线观看| 久久精品综合一区二区三区| 乱系列少妇在线播放| 熟女电影av网| 亚洲aⅴ乱码一区二区在线播放| 啦啦啦韩国在线观看视频| 日韩三级伦理在线观看| 麻豆精品久久久久久蜜桃| 国产一级毛片七仙女欲春2| 丝袜美腿在线中文| 国产又色又爽无遮挡免| 2021天堂中文幕一二区在线观| 久久精品综合一区二区三区| 国产在视频线在精品| 日韩中字成人| 五月天丁香电影| 不卡视频在线观看欧美| 精品午夜福利在线看| 久久久久久久国产电影| 汤姆久久久久久久影院中文字幕 | 久久久久性生活片| 国产av码专区亚洲av| 一区二区三区四区激情视频| 丝瓜视频免费看黄片| 精品一区二区免费观看| 黄色欧美视频在线观看| or卡值多少钱| 国内精品一区二区在线观看| 国产在视频线精品| 日韩人妻高清精品专区| 一级av片app| 精品国产三级普通话版| 在线观看人妻少妇| 国产精品日韩av在线免费观看| 尾随美女入室| 男女视频在线观看网站免费| 亚洲欧美精品自产自拍| 免费黄色在线免费观看| 五月伊人婷婷丁香| 国产 一区 欧美 日韩| 日韩av免费高清视频| 亚洲人成网站在线观看播放| 啦啦啦中文免费视频观看日本| 伦理电影大哥的女人| 亚洲va在线va天堂va国产| 大陆偷拍与自拍| 毛片一级片免费看久久久久| 精品少妇黑人巨大在线播放| 免费av不卡在线播放| 免费不卡的大黄色大毛片视频在线观看 | 97精品久久久久久久久久精品| 国产白丝娇喘喷水9色精品| 黑人高潮一二区| 欧美成人精品欧美一级黄| 久久久午夜欧美精品| 成年人午夜在线观看视频 | 亚洲国产av新网站| 国产单亲对白刺激| 午夜福利在线观看免费完整高清在| 亚洲婷婷狠狠爱综合网| 国产伦在线观看视频一区| 狂野欧美激情性xxxx在线观看| 中文字幕久久专区| 亚洲av福利一区| 中文字幕制服av| 国产黄频视频在线观看| 97精品久久久久久久久久精品| 亚洲国产色片| 欧美精品国产亚洲| 午夜久久久久精精品| 日本wwww免费看| 日本一本二区三区精品| or卡值多少钱| 欧美不卡视频在线免费观看| 春色校园在线视频观看| 久久精品人妻少妇| 精品99又大又爽又粗少妇毛片| 婷婷六月久久综合丁香| 在线观看人妻少妇| 激情五月婷婷亚洲| 欧美激情国产日韩精品一区| 欧美日韩国产mv在线观看视频 | 九草在线视频观看| 性色avwww在线观看| 国产午夜福利久久久久久| 亚洲国产精品成人久久小说| 亚洲精品日韩av片在线观看| 最近最新中文字幕大全电影3| 久久精品综合一区二区三区| 黄色一级大片看看| 99久久精品国产国产毛片| 日韩一本色道免费dvd| 女人被狂操c到高潮| 看非洲黑人一级黄片| 久久久a久久爽久久v久久| 综合色av麻豆| 亚洲av二区三区四区| 男人爽女人下面视频在线观看| 汤姆久久久久久久影院中文字幕 | 国产亚洲5aaaaa淫片| 亚洲不卡免费看| 久久久国产一区二区| 中国美白少妇内射xxxbb| 色播亚洲综合网| 成年版毛片免费区| 成人性生交大片免费视频hd| 美女主播在线视频| 国产成人a区在线观看| 噜噜噜噜噜久久久久久91| 伦理电影大哥的女人| 久久久久网色| 亚洲成人中文字幕在线播放| 伊人久久国产一区二区| 淫秽高清视频在线观看| 国产白丝娇喘喷水9色精品| 一区二区三区乱码不卡18| 人体艺术视频欧美日本| 毛片女人毛片| 日产精品乱码卡一卡2卡三| 午夜福利成人在线免费观看| 日韩不卡一区二区三区视频在线| 久久久久久久久久久丰满| 亚洲av成人av| 人妻制服诱惑在线中文字幕| 中文乱码字字幕精品一区二区三区 | 午夜免费激情av| 久久久午夜欧美精品| 亚洲av在线观看美女高潮| 色综合站精品国产| 免费av不卡在线播放| 久久精品夜色国产| 精品久久久久久久久亚洲| 国产免费又黄又爽又色| 视频中文字幕在线观看| .国产精品久久| 欧美xxxx性猛交bbbb| 搡女人真爽免费视频火全软件| 国产伦精品一区二区三区四那| 日韩av在线大香蕉| 国产免费一级a男人的天堂| 超碰97精品在线观看| 亚洲av免费在线观看| 国产伦在线观看视频一区| 亚洲欧美清纯卡通| 欧美日韩在线观看h| 国产日韩欧美在线精品| 免费观看性生交大片5| 国产一级毛片七仙女欲春2| 韩国高清视频一区二区三区| 免费不卡的大黄色大毛片视频在线观看 | 日韩大片免费观看网站| 欧美+日韩+精品| 最近手机中文字幕大全| 国产黄频视频在线观看| 久久久久久久亚洲中文字幕| 超碰av人人做人人爽久久| 国产熟女欧美一区二区| 免费看av在线观看网站| 国产一区亚洲一区在线观看| 国产一区二区三区综合在线观看 | 日本与韩国留学比较| 成年av动漫网址| 久久久成人免费电影| 天堂网av新在线| 中国国产av一级| 精品酒店卫生间| 婷婷色综合www| av天堂中文字幕网| 成人亚洲精品一区在线观看 | 亚洲高清免费不卡视频| 熟妇人妻久久中文字幕3abv| 国产黄色小视频在线观看| 天天一区二区日本电影三级| 一区二区三区四区激情视频| 亚洲精品第二区| 国产精品无大码| 麻豆乱淫一区二区| 97人妻精品一区二区三区麻豆| 蜜桃久久精品国产亚洲av| 我的女老师完整版在线观看| 真实男女啪啪啪动态图| 国产视频首页在线观看| 最近中文字幕2019免费版| 三级国产精品片| 在线免费观看的www视频| 天天一区二区日本电影三级| 日韩 亚洲 欧美在线| 国产成人免费观看mmmm| 人妻一区二区av| 国产淫片久久久久久久久| 国产伦精品一区二区三区视频9| 成人欧美大片| 丰满少妇做爰视频| 中文欧美无线码| 国产淫语在线视频| videossex国产| 男女啪啪激烈高潮av片| 国产大屁股一区二区在线视频| 亚洲精品久久午夜乱码| 日韩av不卡免费在线播放| 在线观看人妻少妇| 国产av国产精品国产| 国产乱人视频| 又爽又黄a免费视频| 亚洲精品久久午夜乱码| 男女啪啪激烈高潮av片| 国产色爽女视频免费观看| 女人十人毛片免费观看3o分钟| 国产免费一级a男人的天堂| 久久精品国产亚洲av天美| 青青草视频在线视频观看| 国产精品一区二区性色av| 蜜桃久久精品国产亚洲av| 深爱激情五月婷婷| 成人毛片a级毛片在线播放| 精品久久久久久久久久久久久| 一个人免费在线观看电影| 久久久亚洲精品成人影院| 又黄又爽又刺激的免费视频.|