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

    基于擴張狀態(tài)觀測器的磁懸浮球的反步控制

    2021-05-23 12:21:20李冰林曾勵張鵬銘竺志大
    河北科技大學學報 2021年2期

    李冰林 曾勵 張鵬銘 竺志大

    摘 要:針對單自由度磁懸浮系統(tǒng)的非線性、不確定性和易受擾動等特點,提出一種基于擴張狀態(tài)觀測器的反步控制方法以提高系統(tǒng)的控制性能。在系統(tǒng)受到不確定性擾動的情況下,運用擴張狀態(tài)觀測器實時估計懸浮球的位置、速度和擾動信息,并將這種估計值與控制器設計相結(jié)合,然后采用反步法設計磁懸浮球的懸浮位置跟蹤控制律,以Lyapunov方法證明系統(tǒng)的跟蹤誤差最終有界收斂。仿真結(jié)果表明,ESO-BS控制與PID控制相比,前者系統(tǒng)調(diào)節(jié)時間為0.01 s,后者調(diào)節(jié)時間為0.08 s,明顯偏長,因此,ESO-BS控制的動態(tài)特性要優(yōu)于PID控制。在系統(tǒng)存在不確定性的條件下,所設計的控制律能實現(xiàn)磁懸浮球的穩(wěn)定懸浮,并能根據(jù)要求的懸浮高度位置實現(xiàn)位置跟蹤。

    關鍵詞:控制系統(tǒng)仿真技術(shù);磁懸浮球;擴張狀態(tài)觀測器;反步控制;位置跟蹤

    中圖分類號:TH117.1 文獻標識碼:A

    doi:10.7535/hbkd.2021yx02007

    Back-stepping control based on extended state observer for magnetic levitation ball

    LI Binglin1,2,3,ZENG Li2,ZHANG Pengming3,ZHU Zhida2

    (1.College of Automobile and Traffic Engineering,Nanjing Forestry University,Nanjing,Jiangsu 210037,China;2.College of Mechanical Engineering,Yangzhou University,Yangzhou,Jiangsu 225009,China; 3.Zhejiang Jinggong Science & Technology Company Limited,Shaoxing,Zhejiang 312030,China)

    Abstract:Aiming at the characteristics of non-linearity,uncertainty and susceptibility to disturbance of single-degree-of-freedom magnetic levitation ball system,a back-stepping control method based on extended state observer (ESO-BS) was proposed to improve the control performance of the system.The extended state observer was used to estimate the position,velocity and disturbance information of the levitation ball in real time when the system was disturbed by uncertainty.The estimation was combined with the controller design,then the back-stepping method was used to design the levitation position tracking control law of the magnetic levitation ball,and the Lyapunov method was used to prove the ultimate bounded convergence of the tracking error of the system.The simulation results show that,compared with PID control,the adjusting time of ESO-BS control system is 0.01 s,while that of PID control is 0.08 s,which is obviously longer,so the dynamic characteristic of ESO-BS control is better than that of PID control.When there is uncertainty in the system,the designed control law can realize the stable suspension of the ball,and can realize the position tracking according to the required suspension height position.

    Keywords:

    control system simulation technology; magnetic levitation ball; extended state observer(ESO); back-stepping control; position tracking

    磁懸浮能使機械實現(xiàn)無摩擦、無需潤滑高速運轉(zhuǎn),在交通、國防、工業(yè)等領域有較大的應用前景[1-2]。單自由度的懸浮系統(tǒng)具有非線性、開環(huán)不穩(wěn)定的特性[3],通常方法都是將系統(tǒng)在平衡點處線性化,然后再設計控制算法,如文獻[4—5]采用的PID控制和模糊PID控制,盡管這些控制有一定的魯棒性能,但容易使系統(tǒng)性能變得超調(diào)??紤]模型擾動的魯棒控制雖然對模型的攝動有一定的容許度,但是該控制方法也需要合理選擇靈敏度加權(quán)函數(shù)[6-7]?;?刂祈憫?、對參數(shù)變化不敏感,但需要解決變結(jié)構(gòu)時的抖振問題[8]。神經(jīng)網(wǎng)絡控制需對參數(shù)進行大量訓練,對實時性和硬件要求較高[9-10]。以上控制都建立在懸浮球處在平衡位置小擾動下,一旦外擾強度增大,系統(tǒng)偏離平衡點的幅值過大,再由于系統(tǒng)非線性的因素,會導致系統(tǒng)失去穩(wěn)定,進而懸浮失效。擴張狀態(tài)觀測器(ESO)將系統(tǒng)的內(nèi)外干擾信號視為一種總擾動,并能通過觀測器將這一總擾動進行估計,但在參數(shù)整定時比較復雜[11]。反步控制能夠有效消除參數(shù)變化和外界擾動對系統(tǒng)的影響,可以處理不確定性以及帶有未知參數(shù)的非線性系統(tǒng),在四旋翼飛行器等控制領域有較為廣泛的應用[12-13]。

    本研究針對單自由度的懸浮系統(tǒng),在外界存在擾動的情況下,采用擴張狀態(tài)觀測器取代對這種擾動進行估計,然后基于反步法設計磁懸浮球的位置控制器和移動速度控制器,抑制外界擾動,實現(xiàn)對磁懸浮球高度位置的精確懸浮控制,使得控制的位置誤差最終有界。該方案不需要獲取準確的線性化參數(shù),同時不需要估計系統(tǒng)不確定項的上界,簡化了系統(tǒng)控制,最后通過仿真驗證該控制算法的可行性。

    1 單自由度懸浮系統(tǒng)模型及問題描述

    1.1 系統(tǒng)模型

    1)懸浮球與電磁鐵數(shù)學模型

    磁懸浮控制系統(tǒng)的工作原理如圖1所示, F(x0,I0)是電磁鐵在平衡點(I0,x0)處的電磁力,I0是線圈中的偏置電流。當系統(tǒng)工作時,懸浮球受到擾動,并從平衡位置x0向下偏移微小位移,經(jīng)傳感器檢測,送入到控制器中。控制器緊接著向功率放大器輸入控制信號,然后功率放大器將信號放大處理并向磁極線圈中輸入電流為(I0+i),產(chǎn)生的電磁力將浮子拉回到平衡位置。

    根據(jù)電磁學理論,忽略鐵心的影響,懸浮球上的電磁力為[14]

    F(x,I)=μ0S0(NI2x)2=14μ0n2S0I2x2=KI2x2,

    式中:K=14μ0n2S0,電流I增加,電磁力增加;間隙x增加,電磁力減小。F(x,I)在平衡點(x0,I0)進行泰勒展開,電磁力等于懸浮球重力,即F(x0,I0)=mg,有:

    F(x,I)=mg+μ0n2S0I0I-I0 2x20-μ0n2S0I20x-x0 2x30 。(1)

    設工作過程中系統(tǒng)受到外力fd擾動,被懸浮對象的運動學方程為

    md2xdt2=mg-F(x,I)+fd,

    代入式(1),進一步得到:

    md2xdt2+kxx+kII=fd,(2)

    式中:k,kI分別為磁懸浮的位移剛度系數(shù)和電流剛度系數(shù);kx=-μ0N2I20A/(2x30);kI=μ0N2I20A/(2x20)。

    令x1=x,x·1=x2,u=I,單自由度模型式(2)可表示為

    x·1=x2,x·2=fd-kxxmx1-kIImu,(3)

    式中:fd為外擾動;-kxxmx1視為內(nèi)部擾動;f=fd-kxxmx1視為內(nèi)外擾動之和。

    2)傳感器數(shù)學模型[15]

    光電位置傳感器具有線性特性,信號轉(zhuǎn)換環(huán)節(jié)為電壓放大器,屬于比例環(huán)節(jié),所以其模型可以簡化為

    Gs(s)=xs(s)x(s)≈Ks 。

    3)功率放大器數(shù)學模型

    采用電壓-電流型功率放大器,其為一階慣性環(huán)節(jié),由于時間常數(shù)Ta極小,其數(shù)學模型可以簡化為

    G0(s)=i(s)u(s)=1+TasKa≈1Ka 。

    由以上可知,模型中的傳感器數(shù)學模型的輸入、輸出成線性關系,而功率放大器模型的輸入、輸出成反比例關系。

    1.2 問題描述

    在實際工作過程中,電磁線圈存在漏磁現(xiàn)象,鐵心具有磁阻,存在磁飽和問題,這會增加系統(tǒng)的不確定性,同時外界的環(huán)境包括振動對系統(tǒng)控制都有一定的影響。這些問題都會導致系統(tǒng)建模存在偏差。為了降低這些因素對系統(tǒng)控制穩(wěn)定性的影響,設計了一種擴張狀態(tài)觀測器對系統(tǒng)狀態(tài)變量進行估計。

    2 擴張狀態(tài)觀測器(ESO)

    擴張狀態(tài)觀測器將系統(tǒng)未知的干擾信號擴充成新的狀態(tài)變量,本文將f作為狀態(tài)變量進行擴張。擴張狀態(tài)器的輸出動態(tài)補償干擾項的影響,最終可將含有未知干擾項的系統(tǒng)用非線性反饋化為“積分串聯(lián)”形式,然后用非線性狀態(tài)誤差反饋,設計理想的控制器[16-17]。

    對于二階非線性系統(tǒng),可以表示為

    x·=Ax+Bu+Ed,y=Cx,

    式中:A=010001000〗;B=0b00〗;C=100〗;E=001〗T;d為擾動信息。

    擴張狀態(tài)觀測器的增益為

    z·=Az+Bu+L(y-y^),y=Cx,

    式中:z=z1z2z3〗T,表示x=x1x2x3〗T的估計值;L=β01β02β03〗T,為觀測增益。選擇觀測器的極點在左半平面-w0處,則:

    sI-A=s3+β01s2+β02s+β03=(s+w0)3,(4)

    得β01=3w0,β02=3w20,β01=w30,其中w0為觀測器的帶寬,滿足多項式(4)是Hurwitz的,ESO能夠收斂。觀測器的輸出z3將跟蹤f。帶寬的選擇通常會考慮到信號的精度性和噪聲的敏感性,應進行合理選擇[18]。

    擾動反饋控制律為

    u=(u0-z3)/b0 ,

    此時系統(tǒng)變?yōu)殡p積分系統(tǒng)。

    傳統(tǒng)的非線性反饋控制的控制增益:

    u0=k1(v0-z1)+k1(v·0-z2) ,

    式中:k1=w2c;k2=2wc,為控制增益;wc為系統(tǒng)帶寬。

    3 控制器設計

    將擴張狀態(tài)觀測器的觀測值與控制器設計相結(jié)合,以實現(xiàn)磁懸浮系統(tǒng)對干擾的抑制能力。反步控制器主要是基于Lyapunov穩(wěn)定性理論[19],實現(xiàn)對懸浮位置控制和速度控制。

    Step1 根據(jù)式(3),希望找到一個x2d,使得t→SymboleB@,x1→x1d,根據(jù)擴展狀態(tài)觀測器的估計,即使z1→x1d,令:

    e=x1d-z1 。(5)

    對式(5)進行求導,則:

    e·=x·1d-x·1=x·1d-x2 。

    構(gòu)造Lyapunov函數(shù)如下:

    V1=12e2。(6)

    對式(6)進行求導,則:

    V·1=V1ededt=ee·=e(x·1d-x2)=-k1e2 k1>0,(7)

    所以V·1<0,負定,因此,位移誤差e逐漸趨近于0。

    Step2 由式(7)可知,需使x·1d-x2=-k1e,重新令:

    x2d=x·1d+k1e ,(8)

    使得t→SymboleB@,x2→x2d,再令:

    δ=x2d-x2,(9)

    將式(9)代入式(7)得

    V·1=e(x·1d-x2)=e·1d-(x2d-δ)〗,

    將其代入式(8)得:

    V·1=-k1e2+eδ 。

    對式(9)進行微分,可得:

    δ·=x·2d-x·2 。(10)

    將式(3)、式(8)代入式(10)可得:

    δ·=x··1d+k1e·-(fd-kxxmx1-kIImu)=x··1d+k1(x·1d-x2)-(fd-kxxmx1-kIImu) 。(11)

    再構(gòu)造Lyapunov函數(shù)如下:

    V2=V1+12δ2。(12)

    對式(12)進行微分,則:

    V·2=V·1+δδ·=-k1e2+δ(e+δ·) 。

    令:

    e+δ·=-k2δ, k2>0 ,(13)

    則:V·2=-(k1e2+k2δ2)<0,負定,因此,速度誤差δ逐漸趨近于0。

    Step3 求取控制律,將式(11)代入式(13)得:

    e+x··1d+k1(x·1d-x2)-(fd-kxxmx1-kiimu)=-k2δ 。(14)

    對式(14)進行整理,并同時根據(jù)擴張狀態(tài)觀測器對變量x1,x2的估計,可得控制律為

    u=mkII··1d-k1(x·1d-z2)+(fd-kxxmz1)〗 。(15)

    4 穩(wěn)定性證明

    將式(9)代入式(5)得:

    e·=x·1d-x2=x·1d-(x2d-δ)=x·1d-(x·1d+k1e-δ)=-k1e+δ 。(16)

    將式(15)代入式(10)得:

    δ·=x··1d+k1(x·1d-x2)-fd+kxxmx1+kIImmkII··1d-k1(x·1d-x2)+(fd-kxxmx1)〗=-e-k2δ 。(17)

    聯(lián)立式(16)、式(17)得

    e·δ·〗=-k1 1-1 -k2〗eδ〗 。(18)

    設λ1,λ2為式(18)系數(shù)矩陣的特征值,則有:

    λ1+λ2=-k1-k2,λ1λ2=k1k2+1。

    由于k1,k2>0,可知λ1,λ2<0,系統(tǒng)漸進穩(wěn)定。當t→SymboleB@時,e→0,δ→0。通過調(diào)節(jié)系統(tǒng)參數(shù)k1和k2可以改變所控制變量誤差的收斂速度。

    5 仿真試驗分析

    為了驗證所設計控制器的效果,使懸浮球穩(wěn)定地懸浮在離電磁線圈0.014 m的位置高度,在相同的初始條件下,仿真PID控制、Back-stepping(BS)控制以及基于擴張狀態(tài)觀測器的反步控制(ESO-BS)的控制效果,同時研究在懸浮位置的響應速度情況、抗擾能力以及跟蹤性能控制參數(shù)。ESO中wc=200,w0=300,k1=700;k2=700;PID控制器的參數(shù)為kp=-0.45,ki=0,kd=-0.005,磁懸浮的結(jié)構(gòu)和運行的具體參數(shù)值見表1[20]。

    1)位置響應

    圖2為ESO-BS和PID 2種控制方法的懸浮位置階躍響應曲線,懸浮球都能穩(wěn)定地懸浮在0.014 m的位置處。由圖2可知,ESO-BS的動態(tài)特性要優(yōu)于PID控制,系統(tǒng)調(diào)節(jié)時間為0.01 s,PID控制的調(diào)節(jié)時間明顯較長,為0.08 s,即對于單自由度的磁懸浮系統(tǒng),ESO-BS控制在快速響應上得到了顯著提高,遠優(yōu)于PID控制。

    圖3、圖4為擴張狀態(tài)觀測器的觀測值(z1,z2,f)與實際值的比較,可以看出,不論是觀測出的位置信息、速度信息(位置的微分)、內(nèi)外擾信息,兩者之間的差值均很小,擴張狀態(tài)觀測器很好地跟蹤到系統(tǒng)實際值、及微分信號和干擾信號。

    2)抗擾控制

    圖5為在0.05~0.06 s時間段,施加擾動值13 N(為其自重的多倍),基于狀態(tài)觀測器的反步控制的控制效果。由圖5可知,加入干擾后,磁懸浮球在加入擾動的時間段有微小波動,但最終仍維持在懸浮的0.014 m處。其他控制方式在加入同樣大的負載后系統(tǒng)無法穩(wěn)定。圖6為擾動條件下擴張狀態(tài)觀測器的觀測值(z1,z2,f)與實際值的比較,針對觀測器輸出的內(nèi)外擾信息,估計值基本能跟蹤到系統(tǒng)實際值。

    3)跟蹤控制

    為了使磁懸浮球高度可調(diào),使其跟蹤所需要的懸浮高度,并驗證控制算法的跟蹤性能,令所跟蹤的信號以0.014 m為基準高度,并對其施加振幅為0.01 m的正弦信號。

    圖7為擴張狀態(tài)觀測器輸出的觀測信號。由圖7可知,控制系統(tǒng)可以使懸浮球及時跟蹤輸入?yún)⒖剂康淖兓?,具有較好的跟蹤性能,同時各個觀測信號都表現(xiàn)出了較好的觀測性能。速度跟蹤方面,從初始階段開始,實際速度值和觀測速度值有很好的一致性,說明其有很好的跟蹤性能,如圖8所示。

    6 結(jié) 語

    針對懸浮球懸浮位置控制,設計了一種基于擴張狀態(tài)觀測器(ESO)的磁懸浮球的反步控制算法,通過觀測ESO內(nèi)外擾動信號,在Back-stepping控制器設計中以直接反饋的形式消除這種擾動,提高了系統(tǒng)的抗干擾能力。在位置響應控制中,該方法比PID控制調(diào)節(jié)時間短,響應快;在抗擾控制中,針對外界施加的擾動力矩,能對擾動進行估計,使懸浮球穩(wěn)定懸浮;跟蹤性能方面,能使磁懸浮球穩(wěn)定跟蹤給定位置信號。因此,該算法具有較好的控制性能。

    本文算法的效果主要體現(xiàn)在軟件仿真上,未來需在具體實物平臺上進行實驗驗證。

    參考文獻/References:

    [1] LEE H W,KIM K C,LEE J.Review of maglev train technologies[J].IEEE Transactions on Magnetics,2006,42(7):1917-1925.

    [2] 周天豪,陳磊,祝長生,等.基于自適應變步長最小均方算法的磁懸浮高速電機不平衡補償[J].電工技術(shù)學報,2020,35(9):1900-1911.

    ZHOU Tianhao,CHEN Lei,ZHU Changsheng,et al.Unbalance compensation for magnetically levitated high-speed motors based on adaptive variable step size least mean square algorithm[J].Transactions of China Electrotechnical Society,2020,35(9):1900-1911.

    [3] 呂治國,龍志強.磁懸浮球系統(tǒng)的非線性自適應控制方法[J].控制工程,2020,27(1):127-133.

    LU Zhiguo,LONG Zhiqiang.Method on nonlinear adaptive controller for maglev levitation ball system[J].Control Engineering of China,2020,27(1):127-133.

    [4] 邱洪,黃蘇丹,曹廣忠.基于DSP的磁懸浮球模糊PID數(shù)字控制器[J].武漢理工大學學報(信息與管理工程版),2009,31(6):933-936.

    QIU Hong,HUANG Sudan,CAO Guangzhong.Development of fuzzy PID controller for magnetic levitated ball system based on DSP[J].Journal of Wuhan University of Technology (Information & Management Engineering),2009,31(6):933-936.

    [5] SUMAN D,BHATT R.Robust I-fuzzy controller for magnetic ball levitation[C]//2017 International conference of Electronics,Communication and Aerospace Technology (ICECA).Coimbatore:[s.n.],2017:618-621.

    [6] 李群明,朱伶,徐震.磁懸浮球的魯棒控制器設計[J].中南大學學報(自然科學版),2007,38(5):922-927.

    LI Qunming,ZHU Ling,XU Zhen.Robust controller design of maglev ball system[J].Journal of Central South University (Science and Technology),2007,38(5):922-927.

    [7] 黎恒,肖伸平.磁懸浮球系統(tǒng)的滑模變結(jié)構(gòu)魯棒控制器設計[J].湖南工業(yè)大學學報,2014,28(1):58-61.

    LI Heng,XIAO Shenping.Design of sliding mode variable structure robust controller for magnetic levitation ball system[J].Journal of Hunan University of Technology,2014,28(1):58-61.

    [8] 張井崗,方線偉,趙志誠.磁懸浮球系統(tǒng)的分數(shù)階滑??刂芠J].南京理工大學學報(自然科學版),2014,38(1):72-77.

    ZHANG Jinggang,F(xiàn)ANG Xianwei,ZHAO Zhicheng.Fractional order sliding mode control of magnetic levitation ball system[J].Journal of Nanjing University of Science and Technology,2014,38(1):72-77.

    [9] CHIJINDU A T,ANGELA U E.Improving the intelligent control of magnetic levitiation ball using artificial neural network[J].International Journal of Engineering and Computer Science,2019,8(6):24679-24685.

    [10]沈昕璐,莫鑫,周亞南.基于模糊神經(jīng)網(wǎng)絡的磁懸浮球位置控制研究[J].農(nóng)業(yè)裝備與車輛工程,2018,56(5):59-64.

    SHEN Xinlu,MO Xin,ZHOU Yanan.Magnetic levitation ball position control research based on fuzzy neural network[J].Agricultural Equipment & Vehicle Engineering,2018,56(5):59-64.

    [11]張鋆豪,張文安.磁懸浮球系統(tǒng)的線性自抗擾控制與參數(shù)整定[J].系統(tǒng)科學與數(shù)學,2017,37(8):1741-1756.

    ZHANG Junhao,ZHANG Wenan.Linear active disturbance rejection control to magnetic levitation ball system and parameters tuning[J].Journal of Systems Science and Mathematical Sciences,2017,37(8):1741-1756.

    [12]姚達琛,堵威,唐漾.基于擴張狀態(tài)觀測器的四旋翼積分滑模反步控制[J].東華大學學報(自然科學版),2020,46(2):260-268.

    YAO Dachen,DU Wei,TANG Yang.Backstepping controller based on integral sliding mode controller for quadrotor with extended state observer[J].Journal of Donghua University(Natural Science),2020,46(2):260-268.

    [13]費愛玲,李檸,李少遠.固定翼無人機的自抗擾反步控制[J].控制理論與應用,2016,33(10):1296-1302.

    FEI Ailing,LI Ning,LI Shaoyuan.Active disturbance rejection back-stepping control of fixed-wing unmanned aerial vehicle[J].Control Theory and Applications,2016,33(10):1296-1302.

    [14]王曉樂.磁懸浮球系統(tǒng)控制算法的研究[D].成都:西南交通大學,2016.

    WANG Xiaole.Research of Control Algorithm of Magnetic Levitation Ball Control System[D].Chengdu:Southwest Jiaotong University,2016.

    [15]錢玉恒,楊亞非,張翔.離散滑模變結(jié)構(gòu)控制在磁懸浮球系統(tǒng)中的應用[J].中國慣性技術(shù)學報,2012,20(3):368-372.

    QIAN Yuheng,YANG Yafei,ZHANG Xiang.Application of discrete sliding mode variable structure control in magnetic levitation ball system[J].Journal of Chinese Inertial Technology,2012,20(3):368-372.

    [16]HUA X X,HUANG D R,GUO S H.Extended state observer based on ADRC of linear system with incipient fault[J].International Journal of Control,Automation and Systems,2020,18(6):1425-1434.

    [17]夏先齊,張葆,李賢濤,等.基于擴張狀態(tài)觀測器的永磁同步電機低速滑??刂芠J].光學精密工程,2019,27(12):2628-2638.

    XIA Xianqi,ZHANG Bao,LI Xiantao,et al.Low speed sliding mode control of permanent magnet synchronous motor based on extended state observer[J].Optics and Precision Engineering,2019,27(12):2628-2638.

    [18]TIAN G,GAO Z.Frequency response analysis of active disturbance rejection based control system[C]// IEEE International Conference on Control Applications.[S.l.]:IEEE,2007:1595-1599.

    [19]黃書童,江駒,李歐迅,等.近空間可變翼飛行器小翼切換滑模反步控制[J].飛行力學,2020,38(2):29-33.

    HUANG Shutong,JIANG Ju,LI Ouxun,et al.Sliding mode backstepping control for near space morphing vehicle winglet switching[J].Flight Dynamics,2020,38(2):29-33.

    [20]曹廣忠,潘劍飛,黃蘇丹,等.磁懸浮控制系統(tǒng)控制算法及實現(xiàn)[M].北京:清華大學出版社,2013.

    av天堂中文字幕网| 成人三级做爰电影| 国产黄片美女视频| 中文字幕人妻丝袜一区二区| 超碰成人久久| 欧美色欧美亚洲另类二区| 天天添夜夜摸| 亚洲精品一卡2卡三卡4卡5卡| 男女做爰动态图高潮gif福利片| 91麻豆av在线| 亚洲av成人av| 欧美中文综合在线视频| 亚洲 欧美 日韩 在线 免费| 免费观看的影片在线观看| 黄色 视频免费看| 男女那种视频在线观看| 男女视频在线观看网站免费| av天堂在线播放| 色尼玛亚洲综合影院| or卡值多少钱| 少妇人妻一区二区三区视频| 国产av麻豆久久久久久久| 最新在线观看一区二区三区| 国产亚洲av嫩草精品影院| 偷拍熟女少妇极品色| 午夜福利在线在线| 曰老女人黄片| 精品一区二区三区av网在线观看| 久久精品aⅴ一区二区三区四区| 亚洲av中文字字幕乱码综合| tocl精华| 97碰自拍视频| 在线观看免费午夜福利视频| 亚洲av美国av| 夜夜躁狠狠躁天天躁| 国产精品,欧美在线| 国内精品一区二区在线观看| 黄片大片在线免费观看| 国产精品久久久久久精品电影| 成人性生交大片免费视频hd| 国产爱豆传媒在线观看| 狂野欧美白嫩少妇大欣赏| 亚洲自偷自拍图片 自拍| 午夜福利在线观看吧| 欧美在线黄色| xxxwww97欧美| 97超视频在线观看视频| 757午夜福利合集在线观看| 88av欧美| 国产av一区在线观看免费| 久久久久久人人人人人| 一卡2卡三卡四卡精品乱码亚洲| 亚洲av片天天在线观看| 欧美黑人巨大hd| 五月玫瑰六月丁香| av黄色大香蕉| 日本精品一区二区三区蜜桃| 人妻久久中文字幕网| 国产伦人伦偷精品视频| 国产高清videossex| 床上黄色一级片| 免费看光身美女| 欧美另类亚洲清纯唯美| av视频在线观看入口| 国产一区二区三区视频了| 波多野结衣巨乳人妻| 亚洲七黄色美女视频| 999久久久精品免费观看国产| 给我免费播放毛片高清在线观看| 岛国在线观看网站| 国产高清有码在线观看视频| 在线观看免费视频日本深夜| 国产精品美女特级片免费视频播放器 | 午夜久久久久精精品| 亚洲国产欧美一区二区综合| 啦啦啦免费观看视频1| 欧美精品啪啪一区二区三区| 后天国语完整版免费观看| 身体一侧抽搐| 夜夜夜夜夜久久久久| 99国产极品粉嫩在线观看| 很黄的视频免费| 欧美黑人欧美精品刺激| 又黄又粗又硬又大视频| 精品一区二区三区四区五区乱码| 波多野结衣高清作品| 久久精品亚洲精品国产色婷小说| 无人区码免费观看不卡| 五月伊人婷婷丁香| 最新美女视频免费是黄的| 日本与韩国留学比较| 国产三级中文精品| 非洲黑人性xxxx精品又粗又长| 欧美黑人巨大hd| 欧美xxxx黑人xx丫x性爽| 午夜福利在线观看吧| 观看美女的网站| a级毛片在线看网站| 午夜福利高清视频| 在线免费观看不下载黄p国产 | 桃色一区二区三区在线观看| av中文乱码字幕在线| 午夜两性在线视频| 亚洲电影在线观看av| 亚洲中文字幕一区二区三区有码在线看 | 一级黄色大片毛片| 久久精品国产99精品国产亚洲性色| 亚洲成a人片在线一区二区| 91av网站免费观看| 久久久久亚洲av毛片大全| 国产亚洲精品久久久久久毛片| 亚洲 欧美一区二区三区| 国产麻豆成人av免费视频| 欧美激情在线99| 一二三四社区在线视频社区8| 88av欧美| 国产爱豆传媒在线观看| 婷婷精品国产亚洲av在线| 波多野结衣高清无吗| 国产三级中文精品| 久久人人精品亚洲av| 亚洲欧美日韩高清专用| 欧美性猛交╳xxx乱大交人| 免费在线观看视频国产中文字幕亚洲| 欧美成狂野欧美在线观看| 亚洲欧美精品综合一区二区三区| svipshipincom国产片| 午夜亚洲福利在线播放| 久久精品aⅴ一区二区三区四区| 亚洲国产精品sss在线观看| 久久久久久九九精品二区国产| 国产精品亚洲一级av第二区| 欧美日韩瑟瑟在线播放| 亚洲专区字幕在线| 欧美日韩国产亚洲二区| 久久久久久久久中文| 亚洲av中文字字幕乱码综合| 啦啦啦观看免费观看视频高清| 波多野结衣巨乳人妻| av片东京热男人的天堂| 国产美女午夜福利| 给我免费播放毛片高清在线观看| 亚洲精品中文字幕一二三四区| 窝窝影院91人妻| 一本久久中文字幕| 99精品久久久久人妻精品| 久久久久久久午夜电影| 免费av毛片视频| 亚洲av电影在线进入| 黄色视频,在线免费观看| 男女视频在线观看网站免费| 18禁黄网站禁片午夜丰满| www.熟女人妻精品国产| 久久久久久久久免费视频了| 日日摸夜夜添夜夜添小说| 欧美高清成人免费视频www| 国产精品99久久99久久久不卡| 人人妻,人人澡人人爽秒播| 国产精品香港三级国产av潘金莲| 久久久久久久精品吃奶| 色av中文字幕| 国产亚洲精品av在线| 成人欧美大片| 高潮久久久久久久久久久不卡| 亚洲人成网站高清观看| 亚洲国产中文字幕在线视频| 老司机福利观看| 法律面前人人平等表现在哪些方面| 国产成+人综合+亚洲专区| 国产视频一区二区在线看| 18禁黄网站禁片免费观看直播| 99在线人妻在线中文字幕| 别揉我奶头~嗯~啊~动态视频| 好男人在线观看高清免费视频| 给我免费播放毛片高清在线观看| 国产成人av教育| 亚洲av电影不卡..在线观看| 国产伦精品一区二区三区四那| 黑人操中国人逼视频| 日韩国内少妇激情av| 亚洲专区国产一区二区| 三级男女做爰猛烈吃奶摸视频| 色综合婷婷激情| 亚洲 欧美一区二区三区| 精品久久久久久久人妻蜜臀av| 日本三级黄在线观看| 一卡2卡三卡四卡精品乱码亚洲| 国产综合懂色| 757午夜福利合集在线观看| 丁香六月欧美| 国产伦一二天堂av在线观看| 99久国产av精品| 91麻豆av在线| 婷婷丁香在线五月| 丁香六月欧美| 婷婷精品国产亚洲av| 一个人看视频在线观看www免费 | 97超级碰碰碰精品色视频在线观看| 欧美中文日本在线观看视频| www国产在线视频色| 久久午夜综合久久蜜桃| 成人欧美大片| 又黄又爽又免费观看的视频| 久久久久久人人人人人| 久久国产乱子伦精品免费另类| 久久精品91蜜桃| 草草在线视频免费看| 丁香欧美五月| 黄频高清免费视频| 亚洲性夜色夜夜综合| 欧美日韩一级在线毛片| 美女cb高潮喷水在线观看 | 一二三四在线观看免费中文在| 午夜两性在线视频| 亚洲国产精品合色在线| 午夜影院日韩av| 两个人的视频大全免费| 欧美黑人欧美精品刺激| 亚洲激情在线av| av视频在线观看入口| 成人无遮挡网站| 国产精品自产拍在线观看55亚洲| 亚洲精品粉嫩美女一区| 亚洲精品美女久久久久99蜜臀| 欧美zozozo另类| 色综合站精品国产| 中文字幕熟女人妻在线| 美女高潮喷水抽搐中文字幕| 一区福利在线观看| 亚洲精品中文字幕一二三四区| av在线天堂中文字幕| 久久天躁狠狠躁夜夜2o2o| 亚洲精品中文字幕一二三四区| 91久久精品国产一区二区成人 | 日韩欧美 国产精品| 中文字幕人成人乱码亚洲影| 亚洲国产精品999在线| 精华霜和精华液先用哪个| 九九久久精品国产亚洲av麻豆 | 亚洲精华国产精华精| 久久久久久久精品吃奶| 欧美日韩福利视频一区二区| 精品欧美国产一区二区三| 在线播放国产精品三级| 日本三级黄在线观看| 搞女人的毛片| 国产熟女xx| 亚洲一区二区三区色噜噜| 日韩中文字幕欧美一区二区| 国产极品精品免费视频能看的| 国产高清视频在线播放一区| 我的老师免费观看完整版| 国产欧美日韩一区二区三| 国产精品免费一区二区三区在线| 欧美成人性av电影在线观看| 日韩欧美免费精品| 成人鲁丝片一二三区免费| 又爽又黄无遮挡网站| 欧美日本视频| 天天一区二区日本电影三级| 村上凉子中文字幕在线| 天天躁狠狠躁夜夜躁狠狠躁| 国产精品九九99| 国产精品九九99| 国产精品av久久久久免费| 国产伦人伦偷精品视频| 色播亚洲综合网| 99久久精品国产亚洲精品| 欧美色欧美亚洲另类二区| 免费观看精品视频网站| 欧美最黄视频在线播放免费| 国产精品98久久久久久宅男小说| 国产精品99久久久久久久久| 又黄又粗又硬又大视频| 一区二区三区激情视频| 亚洲电影在线观看av| 日本 av在线| 久久99热这里只有精品18| 两个人视频免费观看高清| 99久久综合精品五月天人人| 成人欧美大片| 麻豆一二三区av精品| 免费在线观看影片大全网站| 精品电影一区二区在线| 国产精品香港三级国产av潘金莲| 男女做爰动态图高潮gif福利片| 19禁男女啪啪无遮挡网站| 欧美性猛交╳xxx乱大交人| 69av精品久久久久久| 99精品久久久久人妻精品| 久久久久性生活片| 九色国产91popny在线| 一级毛片精品| 麻豆久久精品国产亚洲av| 精品人妻1区二区| 久久香蕉国产精品| 香蕉丝袜av| 欧美色欧美亚洲另类二区| 熟女少妇亚洲综合色aaa.| 免费在线观看视频国产中文字幕亚洲| 成人国产综合亚洲| 精品国产三级普通话版| 日本精品一区二区三区蜜桃| 亚洲成av人片免费观看| 国产日本99.免费观看| 一区二区三区国产精品乱码| 91在线观看av| 午夜福利18| 他把我摸到了高潮在线观看| 香蕉丝袜av| 久久午夜亚洲精品久久| 欧美丝袜亚洲另类 | 欧美zozozo另类| 精品欧美国产一区二区三| 精品午夜福利视频在线观看一区| 亚洲第一电影网av| 18禁观看日本| 日本五十路高清| 我要搜黄色片| 一进一出抽搐动态| 亚洲色图 男人天堂 中文字幕| 中出人妻视频一区二区| 国产综合懂色| 精品国产超薄肉色丝袜足j| 日本五十路高清| 99久久无色码亚洲精品果冻| 日韩欧美在线乱码| 国产99白浆流出| 久久人人精品亚洲av| 亚洲精品色激情综合| 最新在线观看一区二区三区| 一本精品99久久精品77| 国产视频一区二区在线看| 黄色 视频免费看| 噜噜噜噜噜久久久久久91| 国产高清有码在线观看视频| 成人鲁丝片一二三区免费| 黑人欧美特级aaaaaa片| 亚洲色图av天堂| 18禁观看日本| 深夜精品福利| 国产伦人伦偷精品视频| 亚洲一区高清亚洲精品| 淫妇啪啪啪对白视频| 国产精品电影一区二区三区| 女人高潮潮喷娇喘18禁视频| 听说在线观看完整版免费高清| 69av精品久久久久久| 女人被狂操c到高潮| 国产精品99久久久久久久久| 久久这里只有精品中国| 久久久国产精品麻豆| 一本一本综合久久| 欧洲精品卡2卡3卡4卡5卡区| 三级男女做爰猛烈吃奶摸视频| 91字幕亚洲| 哪里可以看免费的av片| 一本一本综合久久| 91九色精品人成在线观看| 18禁国产床啪视频网站| 熟女人妻精品中文字幕| 色老头精品视频在线观看| 欧美日韩中文字幕国产精品一区二区三区| 亚洲avbb在线观看| 午夜视频精品福利| 国产精品 国内视频| 日韩欧美在线乱码| 日韩大尺度精品在线看网址| 在线观看舔阴道视频| 欧美黄色片欧美黄色片| 国产激情久久老熟女| 午夜福利18| 日本撒尿小便嘘嘘汇集6| 亚洲国产高清在线一区二区三| 国产精品久久久久久人妻精品电影| 91在线观看av| 国产成人福利小说| 91老司机精品| 色综合亚洲欧美另类图片| 精品久久蜜臀av无| 国产高潮美女av| 亚洲黑人精品在线| 很黄的视频免费| 久久久久久久久免费视频了| 1024手机看黄色片| 特级一级黄色大片| 色噜噜av男人的天堂激情| 国产精品精品国产色婷婷| 国产成人欧美在线观看| 一本精品99久久精品77| 两个人的视频大全免费| 欧美成人性av电影在线观看| 可以在线观看毛片的网站| 亚洲欧美日韩无卡精品| 两性夫妻黄色片| 成人欧美大片| 悠悠久久av| 亚洲精品美女久久av网站| 亚洲国产色片| 欧美日韩瑟瑟在线播放| 波多野结衣高清作品| 99久久国产精品久久久| 中文字幕av在线有码专区| 男人舔女人的私密视频| АⅤ资源中文在线天堂| 亚洲成人久久爱视频| 在线十欧美十亚洲十日本专区| 老汉色av国产亚洲站长工具| 精品无人区乱码1区二区| 波多野结衣巨乳人妻| 午夜视频精品福利| 人人妻人人澡欧美一区二区| 亚洲av日韩精品久久久久久密| 成人鲁丝片一二三区免费| 黑人巨大精品欧美一区二区mp4| 亚洲精品色激情综合| www国产在线视频色| h日本视频在线播放| 国产成年人精品一区二区| 熟女人妻精品中文字幕| 国产精品一区二区三区四区免费观看 | 亚洲自拍偷在线| 国产精品九九99| 成人av在线播放网站| 欧美精品啪啪一区二区三区| 色噜噜av男人的天堂激情| 性色av乱码一区二区三区2| 亚洲欧美精品综合久久99| 欧美日本视频| 成人国产一区最新在线观看| 国产一区二区三区在线臀色熟女| 99久久国产精品久久久| 一区二区三区国产精品乱码| 天天躁狠狠躁夜夜躁狠狠躁| 欧美黑人巨大hd| 亚洲欧美日韩高清在线视频| 国产精品久久久久久亚洲av鲁大| 免费大片18禁| 国产高清videossex| 日韩欧美 国产精品| 少妇裸体淫交视频免费看高清| 午夜成年电影在线免费观看| 成人特级黄色片久久久久久久| 麻豆国产97在线/欧美| 少妇的丰满在线观看| 国产精品1区2区在线观看.| 丰满人妻一区二区三区视频av | 日本 av在线| 成人无遮挡网站| 可以在线观看毛片的网站| 亚洲精品一区av在线观看| 免费看美女性在线毛片视频| 免费观看精品视频网站| 51午夜福利影视在线观看| 真实男女啪啪啪动态图| 精华霜和精华液先用哪个| 欧美日韩福利视频一区二区| 欧美精品啪啪一区二区三区| 麻豆国产97在线/欧美| 欧洲精品卡2卡3卡4卡5卡区| 大型黄色视频在线免费观看| 亚洲五月天丁香| 中文字幕高清在线视频| 九九久久精品国产亚洲av麻豆 | 亚洲第一电影网av| 国产精品久久久久久亚洲av鲁大| 欧美xxxx黑人xx丫x性爽| 99久久精品一区二区三区| 亚洲av成人不卡在线观看播放网| 禁无遮挡网站| 看黄色毛片网站| 巨乳人妻的诱惑在线观看| 欧美成狂野欧美在线观看| 真人做人爱边吃奶动态| 男人舔女人下体高潮全视频| 亚洲五月天丁香| 人妻丰满熟妇av一区二区三区| 精品久久久久久,| 老鸭窝网址在线观看| 亚洲av成人一区二区三| 中文字幕人成人乱码亚洲影| 99riav亚洲国产免费| 天堂av国产一区二区熟女人妻| 免费看a级黄色片| 最近最新免费中文字幕在线| 亚洲天堂国产精品一区在线| 日韩成人在线观看一区二区三区| 在线观看免费午夜福利视频| 婷婷亚洲欧美| 嫩草影院精品99| 国产综合懂色| www日本黄色视频网| 国产真实乱freesex| 欧美日韩综合久久久久久 | 国产精品影院久久| 亚洲五月婷婷丁香| av天堂中文字幕网| 啪啪无遮挡十八禁网站| 日本熟妇午夜| 亚洲精品国产精品久久久不卡| 欧美中文综合在线视频| av天堂在线播放| 一边摸一边抽搐一进一小说| 少妇的逼水好多| 亚洲国产欧洲综合997久久,| 九色国产91popny在线| 亚洲美女黄片视频| 亚洲成av人片免费观看| 亚洲欧美日韩东京热| 日韩三级视频一区二区三区| 日本三级黄在线观看| 国产男靠女视频免费网站| 久久九九热精品免费| 国产精品一区二区三区四区免费观看 | 欧美黑人巨大hd| 免费在线观看视频国产中文字幕亚洲| 免费大片18禁| 99国产精品一区二区三区| 亚洲精品中文字幕一二三四区| 一本综合久久免费| 久久中文字幕一级| 欧美精品啪啪一区二区三区| 亚洲精华国产精华精| 免费看光身美女| 51午夜福利影视在线观看| 看片在线看免费视频| 久久99热这里只有精品18| 欧美中文综合在线视频| 午夜福利18| 欧美日韩一级在线毛片| 午夜福利在线观看免费完整高清在 | 黑人巨大精品欧美一区二区mp4| 久久久久久久久久黄片| 亚洲人成伊人成综合网2020| 老熟妇乱子伦视频在线观看| av在线天堂中文字幕| 国产欧美日韩精品一区二区| 变态另类丝袜制服| 91老司机精品| 免费电影在线观看免费观看| 日本在线视频免费播放| 午夜a级毛片| 亚洲第一欧美日韩一区二区三区| 亚洲av成人不卡在线观看播放网| 99精品久久久久人妻精品| 我的老师免费观看完整版| 久久久色成人| 日本成人三级电影网站| 亚洲精华国产精华精| 老司机在亚洲福利影院| 一夜夜www| 久久久久免费精品人妻一区二区| 亚洲精品美女久久av网站| 亚洲在线观看片| 国产精品,欧美在线| 亚洲第一电影网av| 亚洲美女视频黄频| 亚洲国产中文字幕在线视频| 99re在线观看精品视频| 在线免费观看不下载黄p国产 | av欧美777| 观看免费一级毛片| 1024香蕉在线观看| 国产激情欧美一区二区| 久久中文字幕一级| av天堂中文字幕网| 别揉我奶头~嗯~啊~动态视频| 日韩精品青青久久久久久| 中出人妻视频一区二区| 国产伦精品一区二区三区视频9 | 午夜精品在线福利| 久久久久久人人人人人| 又爽又黄无遮挡网站| 国产精品香港三级国产av潘金莲| 最新美女视频免费是黄的| 欧美日韩中文字幕国产精品一区二区三区| 无限看片的www在线观看| 91久久精品国产一区二区成人 | 99久久99久久久精品蜜桃| 国产亚洲欧美在线一区二区| 国产精品野战在线观看| 国产精品精品国产色婷婷| 国产黄片美女视频| 午夜精品久久久久久毛片777| 国产野战对白在线观看| 久久久久国产一级毛片高清牌| 午夜精品在线福利| 每晚都被弄得嗷嗷叫到高潮| 首页视频小说图片口味搜索| 久久天堂一区二区三区四区| 99久久精品一区二区三区| 最近最新中文字幕大全免费视频| 午夜两性在线视频| 久久这里只有精品19| 一a级毛片在线观看| 九色国产91popny在线| 欧美性猛交黑人性爽| 色老头精品视频在线观看| 禁无遮挡网站| 国产精品久久久久久亚洲av鲁大| 校园春色视频在线观看| 十八禁网站免费在线| 欧美日本视频| 国模一区二区三区四区视频 | 欧美乱妇无乱码| 久久久国产欧美日韩av| 91麻豆av在线| АⅤ资源中文在线天堂| 在线十欧美十亚洲十日本专区| 免费观看的影片在线观看| 又粗又爽又猛毛片免费看| 成人三级做爰电影| 国产av不卡久久| 亚洲欧美日韩东京热| 美女cb高潮喷水在线观看 |