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

    軋機(jī)液壓AGC-LP耦合系統(tǒng)的H∞控制器設(shè)計(jì)

    2017-07-07 13:28:01魯凌云李維剛黃衛(wèi)華嚴(yán)???/span>
    關(guān)鍵詞:張應(yīng)力活套機(jī)架

    魯凌云,李維剛,2,黃衛(wèi)華,嚴(yán)???/p>

    (1.武漢科技大學(xué)信息科學(xué)與工程學(xué)院,湖北 武漢,430081;2.武漢科技大學(xué)冶金工業(yè)過程系統(tǒng)科學(xué)湖北省重點(diǎn)實(shí)驗(yàn)室,湖北 武漢,430065)

    ?

    軋機(jī)液壓AGC-LP耦合系統(tǒng)的H∞控制器設(shè)計(jì)

    魯凌云1,李維剛1,2,黃衛(wèi)華1,嚴(yán)???

    (1.武漢科技大學(xué)信息科學(xué)與工程學(xué)院,湖北 武漢,430081;2.武漢科技大學(xué)冶金工業(yè)過程系統(tǒng)科學(xué)湖北省重點(diǎn)實(shí)驗(yàn)室,湖北 武漢,430065)

    根據(jù)熱連軋機(jī)自動厚度控制(AGC)系統(tǒng)和液壓活套控制(LP)系統(tǒng)各自的特點(diǎn),結(jié)合這兩個子系統(tǒng)的模型,建立液壓AGC-LP耦合系統(tǒng)模型。參照廣義系統(tǒng)的標(biāo)準(zhǔn)選取了合適的控制目標(biāo),基于H∞控制理論對組合系統(tǒng)進(jìn)行分析,采用MATLAB魯棒控制工具箱中線性矩陣不等式(LMI)求解方法設(shè)計(jì)了H∞狀態(tài)反饋控制器,使得系統(tǒng)達(dá)到性能指標(biāo)最優(yōu)。將使用H∞控制器與傳統(tǒng)PID控制器的AGC-LP耦合系統(tǒng)進(jìn)行仿真比較,結(jié)果表明,使用H∞控制器的系統(tǒng)有更好的動態(tài)和穩(wěn)態(tài)性能,兩個子系統(tǒng)相關(guān)變量的耦合也有所改善。

    熱連軋機(jī);AGC;活套系統(tǒng);H∞控制器;LMI;解耦

    為了保證軋制過程的順利進(jìn)行,現(xiàn)代主流的熱連軋機(jī)大都通過液壓活套裝置(LP)的緩沖作用來吸收咬鋼過程中形成的套量,并保持恒定的小張應(yīng)力控制,同時通過自動厚度控制系統(tǒng)(AGC)來提高帶鋼厚度控制的精度。由于張應(yīng)力的變化是引起帶鋼厚度變化的主要原因之一,因此AGC-LP系統(tǒng)的控制難點(diǎn)在于LP系統(tǒng)的高度控制和張應(yīng)力控制之間的耦合以及AGC系統(tǒng)與LP系統(tǒng)之間的相互影響?,F(xiàn)有對AGC-LP系統(tǒng)耦合的研究大都考慮了組合系統(tǒng)的線性成分和非線性成分,并改善了系統(tǒng)的魯棒性,但組合模型強(qiáng)耦合非線性的特性為系統(tǒng)的解耦帶來了困難。

    由于H∞控制理論考慮了實(shí)際系統(tǒng)與數(shù)學(xué)模型之間的不確定性,降低了系統(tǒng)建模的復(fù)雜程度,且H∞方法將復(fù)雜的系統(tǒng)解耦問題轉(zhuǎn)換成了線性矩陣不等式求解問題,并且利用MATLAB的LMI工具箱可以精確地求解,為此,本文基于H∞控制理論對AGC-LP耦合模型進(jìn)行控制器設(shè)計(jì)。首先,根據(jù)AGC和LP系統(tǒng)子模型建立了AGC-LP耦合系統(tǒng)模型,然后通過H∞控制方法為AGC-LP控制系統(tǒng)設(shè)計(jì)了狀態(tài)反饋控制器,最后對給定輸入信號下組合系統(tǒng)的輸出響應(yīng)進(jìn)行仿真,并通過與傳統(tǒng)PID控制器相關(guān)性能指標(biāo)的比較來驗(yàn)證所設(shè)計(jì)控制器的有效性。

    1 AGC-LP耦合系統(tǒng)模型的建立

    以熱連軋機(jī)相鄰機(jī)架作為分析對象,涉及到的變量參數(shù)如圖1所示。圖1中,l1和l2分別為活套支撐輥接觸面到Fi機(jī)架出口和Fi+1機(jī)架入口的距離;σ為Fi機(jī)架出口帶鋼張應(yīng)力;H和h分別為Fi機(jī)架入口和出口的帶鋼厚度;vi和vi+1分別為帶鋼在Fi機(jī)架的出口速度和Fi+1機(jī)架的入口速度;R為軋輥的半徑;θ1和θ2分別為Fi機(jī)架出口和Fi+1機(jī)架入口的連線與軋線之間的夾角;l為活套臂長;θ為活套高度角;a為活套器中心到相鄰機(jī)架的直線距離;L=2a;ω為軋輥的角速度;d為活套支撐輥與軋線的距離。

    圖1 活套機(jī)構(gòu)參數(shù)示意圖

    1.1 液壓閥控缸系統(tǒng)的建立

    LP系統(tǒng)和AGC系統(tǒng)均由液壓閥控缸驅(qū)動,兩者的液壓閥控缸模型相同,只是負(fù)載和對象不同。關(guān)于閥控缸的建模過程詳見文獻(xiàn)[1],其中幾個關(guān)鍵方程如下:伺服閥的動態(tài)特性可以用一個二階振蕩環(huán)節(jié)來描述:

    (1)

    式中:Gv為傳遞函數(shù);s為拉普拉斯算子;ωv為伺服閥的截止頻率;δv為伺服閥閥芯的阻尼系數(shù)。

    液壓缸的流量方程為

    QL=KqXv-KcΔPL

    (2)

    式中:Kq為伺服閥的流量增益系數(shù);Kc為伺服閥的流量壓力系數(shù);Xv為伺服閥的位移量;ΔPL為系統(tǒng)壓力與回油壓力之差。

    根據(jù)牛頓第二定律,液壓缸的運(yùn)動方程為

    (3)

    式中:Apst、Arod分別為液壓缸無桿腔和有桿腔的橫截面積;Ppst、Prod分別為液壓缸無桿腔和有桿腔的壓力;M為負(fù)載的等效質(zhì)量;B為運(yùn)動黏滯系數(shù);G為彈性阻力系數(shù);y為活塞位移;t為時間;fLd為負(fù)載作用力。

    1.2 AGC系統(tǒng)模型

    Fi機(jī)架出口帶鋼厚度可用彈跳方程來表示,彈跳方程給出了活套張應(yīng)力與活套高度hi的關(guān)系[2]:

    (4)

    式中:S為輥縫設(shè)定值;Bw為帶鋼寬度;Y為變形抗力系數(shù);Q為應(yīng)力狀態(tài)系數(shù);τ為入口帶鋼張應(yīng)力;C為軋機(jī)剛度系數(shù)。

    對于AGC系統(tǒng)來說,式(3)所示力平衡方程中ΔF即為系統(tǒng)提供的軋制力。如果只考慮軋件變形抗力而忽略其他因素,則負(fù)載作用力fLd可以表示為

    (5)

    若忽略軋機(jī)的軋輥偏心量Δx,則軋件的出口厚度h便等于工作輥的負(fù)載輥縫SL,而SL為

    (6)

    1.3 LP系統(tǒng)模型的關(guān)鍵方程

    設(shè)第Fi機(jī)架的出口速度為Vout(i),第Fi+1機(jī)架的入口速度為Vin(i+1),那么活套的套量Lloop可以表示為

    (7)

    張應(yīng)力σ可以表示為

    (8)

    式中:E為材料的彈性模數(shù)。

    熱連軋自動厚度控制過程中,對活套的補(bǔ)償有本架補(bǔ)償和后架補(bǔ)償兩種方式,采取后架補(bǔ)償即在第Fi+1機(jī)架投入AGC能獲得更高的速度補(bǔ)償精度,因此有如下的速度補(bǔ)償公式[3]:

    (9)

    式中:vi為Fi機(jī)架帶鋼的出口速度;vRi為Fi機(jī)架的軋輥速度;vi+1為Fi+1機(jī)架帶鋼的入口速度;vR(i+1)為Fi+1機(jī)架的軋輥速度;fi為Fi機(jī)架的前滑值;βi+1為Fi+1機(jī)架的后滑值。

    由式(7)~式(9)可得活套套量和帶鋼速度以及張應(yīng)力的關(guān)系為

    (10)

    式中:

    活套的角加速度與綜合力矩M密切相關(guān),對于活套有如下的力矩平衡關(guān)系:

    (11)

    (12)

    1.4 AGC-LP模型的建立

    綜合式(1)~式(12),已經(jīng)找出了關(guān)聯(lián)變量,建立起了自動厚度控制系統(tǒng)與活套系統(tǒng)的聯(lián)系。將AGC系統(tǒng)寫成狀態(tài)空間表達(dá)式,選取輥縫SL、帶鋼速度V、活套輸出力矩M作為輸入控制量,即u=[u1u2u3]T=[SLVM]T,令張應(yīng)力σ、活套高度角θ、軋輥角速度ω作為狀態(tài)變量,即x=[x1x2x3]T=[σθω]T,以帶鋼出口厚度h、張應(yīng)力σ和活套高度角θ作為輸出量,即y=[y1y2y3]T=[hσθ]T,則可得AGC-LP的狀態(tài)方程為

    (13)

    輸出方程為

    (14)

    2 H狀態(tài)反饋控制器設(shè)計(jì)

    將AGC-LP系統(tǒng)視作如圖2所示的廣義系統(tǒng),其中P(s)為系統(tǒng)狀態(tài)空間矩陣,K為反饋控制器矩陣,K(s)為控制器的傳遞函數(shù)。整個閉環(huán)系統(tǒng)由以下的狀態(tài)空間描述:

    (15)

    式中:x∈n為狀態(tài)向量;u∈n為控制輸入;y∈p為測量輸出,A、C1、C2為系統(tǒng)狀態(tài)變量的系數(shù)矩陣;B1、D11、D21為系統(tǒng)擾動的系數(shù)矩陣;B2、D12、D22為系統(tǒng)輸入的系數(shù)矩陣。

    圖2 廣義系統(tǒng)

    由于AGC-LP系統(tǒng)是可觀測的,因此y較容易得到;在H∞控制器的設(shè)計(jì)過程中,測量輸出都是作為系統(tǒng)的反饋輸入,z∈r是有限大的被調(diào)輸出,也是在系統(tǒng)的仿真過程中所得到的輸出響應(yīng);ω∈q為具有有限能量的外部擾動,這里將以有限能量的白噪聲作為外部擾動。

    H∞控制的目的就是設(shè)計(jì)一個控制器u(s)=K(s)y(s),使得閉環(huán)系統(tǒng)滿足以下性質(zhì):

    (1)閉環(huán)系統(tǒng)是內(nèi)部穩(wěn)定的,即閉環(huán)系統(tǒng)狀態(tài)矩陣的所有特征值均在左半開復(fù)平面中。

    (2)從擾動輸入ω到被調(diào)輸出z的閉環(huán)傳遞函數(shù)Twz(s)的H∞范數(shù)小于1,即

    (16)

    具有這樣性質(zhì)的控制器u(s)=K(s)y(s)稱為式(15)所示系統(tǒng)的一個H∞控制器[2]。

    傳統(tǒng)的活套系統(tǒng)都采用PID控制器,沒有考慮活套高度和張應(yīng)力之間的相互作用,系統(tǒng)的抗擾動性能有所欠缺。相比之下,H∞狀態(tài)反饋設(shè)計(jì)是在Ricatti方程的求解條件上增加由Bode積分定理所規(guī)定的魯棒性約束,即達(dá)到式(16)所要求的性能指標(biāo),以此為目標(biāo)來設(shè)計(jì)H∞反饋控制以消除建模誤差和對象零極點(diǎn)漂移的影響。

    由于系統(tǒng)的狀態(tài)是可以直接測量到的,因此采用H∞狀態(tài)反饋控制器。首先,要滿足以下的反饋控制率:

    u=Kx

    使得相應(yīng)的閉環(huán)系統(tǒng)

    是漸近穩(wěn)定的。

    閉環(huán)傳遞函數(shù)Twz(s)的H范數(shù)滿足:

    (17)

    當(dāng)且僅當(dāng)存在一個對稱正定矩陣X和矩陣W,使得矩陣不等式

    (18)

    成立,則u=W*(X*)-1x是式(15)所示系統(tǒng)的一個狀態(tài)反饋H∞控制器[4]。

    在加入了H∞狀態(tài)反饋控制器后AGC-LP系統(tǒng)的框圖如圖3所示。

    要確定矩陣不等式(18)中的各子矩陣,由式(13)可得:

    圖3 帶H∞反饋控制器的AGC-LP系統(tǒng)框圖

    (19)

    取活套高度角為10°,彈性模量E=2.1×105MPa,GD2=1.5×105kN·m;σ=25.9 MPa,L、l1、l2、R分別為6、2.2、1.8、0.8 m,J=78.5 N·m,B=8.5,Bw=0.082,C=8.25,H=0.3m,h=0.156 m,Ks=50,則可求得系統(tǒng)矩陣A、B2、C1,而外部擾動的系數(shù)矩陣B1、D11、D12未知;系統(tǒng)的測量輸出對應(yīng)的系數(shù)矩陣C2、D21、D22均可通過觀測得到。在熱連軋系統(tǒng)中,擾動的因素很多,例如溫度、軋輥形變等,這些擾動綜合起來影響軋制帶鋼的厚度,基于這一點(diǎn),用近似的白噪聲和正弦擾動來模擬帶鋼在軋制過程中的厚度波動,則可求得系數(shù)矩陣B1、D11。

    得到了式(19)所需的各個矩陣后,再利用MATLAB的LMI工具箱相關(guān)函數(shù)進(jìn)行求解。因?yàn)锳GC-LP系統(tǒng)的優(yōu)化問題是一個具有線性矩陣不等式約束的線性目標(biāo)函數(shù)最小化問題,故使用工具箱中的求解器mincx,將AGC-LP系統(tǒng)的狀態(tài)空間矩陣代入到式(18)的矩陣不等式進(jìn)行求解。經(jīng)過反復(fù)的迭代運(yùn)算,得到系統(tǒng)的最優(yōu)H∞性能指標(biāo)γ1的值如表1所示。

    表1 線性目標(biāo)在LMI約束下的最小化

    Table 1 Minimization of linear target under LMI constraints

    迭代計(jì)算次數(shù)當(dāng)前γ1的最小值16642662122384.73102091.8487463.29??479448.017279488165.4510004981.797114??880.404939890.410422

    最終,γ1達(dá)到最小值0.410 422,得到了滿足H∞性能指標(biāo)的最優(yōu)解,可以求得優(yōu)化系統(tǒng)性能的閉環(huán)反饋控制矩陣,為了方便仿真,再通過MATLAB指令:

    [gopt,k]=hinflmi(P,r)

    (20)

    得到系統(tǒng)的最優(yōu)H∞控制器K(s),對應(yīng)的系統(tǒng)矩陣可以由函數(shù)[ak,bk,ck,dk]=ltiss(K(s))求得,這樣就能得到一個使H∞性能指標(biāo)γ1最優(yōu)的控制器矩陣。

    3 仿真結(jié)果及分析

    3.1 H∞控制器與PID控制器對AGC-LP系統(tǒng)響應(yīng)的影響

    采用Simulink軟件進(jìn)行仿真,對使用H∞控制器與PID控制器時AGC-LP系統(tǒng)的輸出進(jìn)行比較。按照圖3所示的系統(tǒng)結(jié)構(gòu)搭建閉環(huán)仿真框圖,主控制器K(s)的系統(tǒng)矩陣和液壓AGC-LP的系統(tǒng)矩陣用兩個狀態(tài)空間模塊表示,其中主控制器的系統(tǒng)矩陣分別與式(20)所求得的ak、bk、ck、dk相對應(yīng)。

    以鞍鋼F3機(jī)架出口為例,設(shè)定帶鋼出口厚度為7.52 mm,張應(yīng)力為25.9 MPa,活套高度角為9°。經(jīng)過仿真,在持續(xù)擾動下,AGC活套系統(tǒng)的張應(yīng)力、活套高度角、軋件厚度3個變量在傳統(tǒng)PID控制和H∞控制下的響應(yīng)曲線如圖4所示。由圖4可以看出,對于活套高度角和軋件厚度的響應(yīng),使用H∞控制器不但減少了震蕩和調(diào)節(jié)時間,同時活套角的超調(diào)量減小,軋件厚度的穩(wěn)態(tài)精度也得到提高,張應(yīng)力的調(diào)節(jié)時間由0.9 s縮短到0.6 s??傮w上看,H∞控制較PID控制調(diào)節(jié)時間縮短,過渡過程比較平穩(wěn),過渡過程中的震蕩次數(shù)減少,震蕩幅度也有所降低,有效避免了張應(yīng)力的瞬間突變,超調(diào)及響應(yīng)時間也同時縮短,可以控制熱連軋過程盡快進(jìn)入穩(wěn)態(tài)。

    (a)張應(yīng)力 (b)活套高度角 (c)軋件厚度

    圖4 系統(tǒng)在H∞和PID控制方式下的輸出比較

    Fig.4 Comparison of output of system under the H infinity and PID control

    3.2 H∞控制器對AGC-LP系統(tǒng)的解耦效果

    對活套張應(yīng)力和活套高度角分別附加1 MPa和3°的階躍擾動,仿真對比使用H∞控制器和傳統(tǒng)PID控制器時AGC-LP系統(tǒng)的響應(yīng),以驗(yàn)證H∞控制器對AGC-LP系統(tǒng)的解耦效果,結(jié)果如圖5所示。由圖5中可以看出,與PID控制器相比,本文所設(shè)計(jì)H∞控制器對AGC-LP系統(tǒng)的解耦有較大的幫助。

    (a)高度角受擾時高度角和張應(yīng)力響應(yīng)曲線

    (b)張應(yīng)力受擾時高度角和張應(yīng)力響應(yīng)曲線

    3.3 H∞控制器控制下AGC系統(tǒng)的輸出響應(yīng)

    為了驗(yàn)證加入H∞控制器后液壓AGC系統(tǒng)的厚度控制效果,按照實(shí)際情況,仿真中在帶鋼進(jìn)入精軋前施加擾動,主要考慮中間坯頭尾厚度波動、溫度以及高頻隨機(jī)擾動等干擾信號模擬帶鋼厚度變化情況。設(shè)帶鋼原始厚度為H0,則帶鋼厚度的波動可近似表示為

    模擬鞍鋼集團(tuán)朝陽鞍凌鋼鐵有限公司熱連軋機(jī)組F1~F5機(jī)架的出口厚度給定量[5],每一個機(jī)架都按照各自的設(shè)定給定相應(yīng)的輥縫、活套高度角以及厚度值,仿真結(jié)果如圖6所示。由圖6中可知,雖然存在來料厚度不均勻以及壓力和彈跳的波動等影響因素,但通過H∞控制器控制的AGC系統(tǒng)的調(diào)節(jié),機(jī)架出口的帶鋼厚度基本穩(wěn)定,由此證明了系統(tǒng)模型的準(zhǔn)確性與控制器的有效性。

    圖6 F1~F5機(jī)架出口厚度仿真結(jié)果

    Fig.6 Simulation effort of F1~F5 mills’ outgoing thickness

    4 結(jié)語

    本文研究了基于LMI的狀態(tài)反饋H∞魯棒控制方法在AGC-LP控制系統(tǒng)中的應(yīng)用,通過建立系統(tǒng)的狀態(tài)空間模型,使用LMI工具箱進(jìn)行優(yōu)化求解,讓系統(tǒng)性能指標(biāo)優(yōu)化的過程更加簡便。仿真結(jié)果表明采用H∞控制理論設(shè)計(jì)的控制器不僅可以讓系統(tǒng)具有良好的性能指標(biāo),還對耦合系統(tǒng)有較好的解耦效果,這對于AGC-LP系統(tǒng)性能的改進(jìn)有一定的實(shí)際意義。

    [1] 童朝南,武延坤,宗勝悅,等. 熱連軋中液壓活套系統(tǒng)數(shù)學(xué)模型的研究[J]. 系統(tǒng)仿真學(xué)報(bào), 2008, 20(6):1381-1385,1389.

    [2] 曲蕾,王京. 多變量非線性厚度-活套系統(tǒng)的魯棒逆控制[J]. 控制理論與應(yīng)用, 2009, 26(5):562-566.

    [3] 童朝南,孫一康,陳百紅. 熱連軋AGC控制中活套補(bǔ)償?shù)膬煞N觀點(diǎn)[J]. 軋鋼, 2002, 19(4):47-48.

    [4] 俞立. 魯棒控制-線性矩陣不等式處理方法[M].北京:清華大學(xué)出版社, 2002.

    [5] 莊野,張輝,劉新婷. 鞍凌熱連軋自動厚度控制系統(tǒng)實(shí)時仿真[J]. 冶金自動化, 2013, 37(6):54-58.

    [責(zé)任編輯 鄭淑芳]

    Design of H∞controller for hydraulic AGC-looper coupling system for hot strip mills

    LuLingyun1,LiWeigang1,2,HuangWeihua1,YanBaokang1

    (1. College of Information Science and Engineering, Wuhan University of Science and Technology, Wuhan 430081, China; 2. Hubei Province Key Laboratory of Systems Science in Metallurgical Process, Wuhan University of Science and Technology,Wuhan 430065, China )

    According to the characteristics of the automatic gauge control(AGC) system and the looper control system of hot strip mills, a model of the hydraulic AGC-LP system was constructed by the models of two subsystems. The appropriate control target was selected by the criterion of the generalized system to establish the generalized controlled object, and then the H∞state feedback controller was designed to achieve the optimal performance, with the linear matrix inequality(LMI) of robust control toolbox in Matlab. Comparing the PID control system with H∞controll system in simulation, the results show that the system with the designed H∞controller has a better dynamic and static performance, in which the relevant decoupling performance among variables of two subsystem has been improved.

    hot strip mill; AGC; looper system; H∞controller; LMI; decouple

    2017-03-25

    湖北省教育廳科學(xué)技術(shù)研究計(jì)劃重點(diǎn)項(xiàng)目(D20161103);武漢市青年科技晨光計(jì)劃資助項(xiàng)目(2016070204010099);武漢科技大學(xué)冶金工業(yè)過程系統(tǒng)科學(xué)湖北省重點(diǎn)實(shí)驗(yàn)室開放基金資助項(xiàng)目(Z201501).

    魯凌云(1993-),男,武漢科技大學(xué)碩士生.E-mail:121063261@qq.com

    李維剛(1977-),男,武漢科技大學(xué)教授,博士.E-mail:liweigang.luck@foxmail.com

    10.3969/j.issn.1674-3644.2017.04.008

    TP273

    A

    1674-3644(2017)04-0279-06

    猜你喜歡
    張應(yīng)力活套機(jī)架
    熱連軋機(jī)精軋活套輥斷裂分析與改進(jìn)
    礦冶(2020年5期)2020-12-29 02:55:38
    無張力軋制在棒線材上的應(yīng)用
    什么是不銹鋼的應(yīng)力腐蝕開裂?
    別忽略它的存在!“意大利新一代架皇”BAS Accordeon(雅歌頓)XL4 2.0發(fā)燒機(jī)架
    冷軋立式活套控制原理
    山西冶金(2019年1期)2019-04-26 08:25:32
    內(nèi)聚力-張力學(xué)說中關(guān)于負(fù)壓的幾點(diǎn)疑問
    異常血流動力對TLR4/NF—κB信號傳導(dǎo)通路及其下游炎癥因子的影響
    熱軋拉矯機(jī)機(jī)架加工討論
    雙機(jī)架平整機(jī)板形控制算法及其應(yīng)用
    上海金屬(2013年6期)2013-12-20 07:58:02
    全連續(xù)冷連軋機(jī)組活套設(shè)計(jì)
    欧美+日韩+精品| 内地一区二区视频在线| 欧美成人a在线观看| 一个人观看的视频www高清免费观看| 中国国产av一级| 国产高清视频在线观看网站| 91久久精品电影网| 国产一级毛片在线| 亚洲欧美日韩东京热| 欧美又色又爽又黄视频| 丰满乱子伦码专区| 美女 人体艺术 gogo| 欧美一级a爱片免费观看看| 99久久九九国产精品国产免费| 国产精品福利在线免费观看| 日本免费a在线| 国产中年淑女户外野战色| 欧美最黄视频在线播放免费| 日本黄色视频三级网站网址| 久久精品影院6| 最后的刺客免费高清国语| 能在线免费观看的黄片| 婷婷色av中文字幕| 国产精品美女特级片免费视频播放器| 亚洲国产精品成人综合色| 99热只有精品国产| 男女下面进入的视频免费午夜| 日日啪夜夜撸| 91午夜精品亚洲一区二区三区| av国产免费在线观看| 乱人视频在线观看| 久久草成人影院| 伦精品一区二区三区| 国产乱人视频| 在线观看一区二区三区| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 午夜福利高清视频| 色5月婷婷丁香| 日韩在线高清观看一区二区三区| 日产精品乱码卡一卡2卡三| 亚洲国产欧美人成| 免费在线观看成人毛片| 国产精品人妻久久久久久| 亚洲成人久久性| 国产三级在线视频| 高清午夜精品一区二区三区 | 亚洲性久久影院| 久久99热6这里只有精品| 乱系列少妇在线播放| 嫩草影院精品99| 中文精品一卡2卡3卡4更新| 岛国在线免费视频观看| 国产黄片美女视频| 国产精品不卡视频一区二区| 在线免费观看不下载黄p国产| 午夜免费激情av| 日本欧美国产在线视频| 人妻系列 视频| av专区在线播放| 边亲边吃奶的免费视频| 淫秽高清视频在线观看| 亚洲一区高清亚洲精品| 在线播放无遮挡| 国产乱人视频| 最近视频中文字幕2019在线8| 最近2019中文字幕mv第一页| 一本精品99久久精品77| 高清日韩中文字幕在线| 欧美最新免费一区二区三区| 99久久九九国产精品国产免费| av在线亚洲专区| 小蜜桃在线观看免费完整版高清| 十八禁国产超污无遮挡网站| 亚洲国产精品合色在线| 欧美bdsm另类| 综合色av麻豆| 性色avwww在线观看| 久久99热6这里只有精品| av在线老鸭窝| 午夜视频国产福利| 少妇高潮的动态图| 国产高潮美女av| 长腿黑丝高跟| 啦啦啦啦在线视频资源| 亚洲性久久影院| 12—13女人毛片做爰片一| 网址你懂的国产日韩在线| 国产成人一区二区在线| 国产伦精品一区二区三区四那| 天堂网av新在线| 少妇的逼水好多| 欧美变态另类bdsm刘玥| 在线a可以看的网站| 久久99热这里只有精品18| 精品久久久久久久末码| 青春草视频在线免费观看| 国产高潮美女av| 久久久精品94久久精品| 午夜福利成人在线免费观看| 伦理电影大哥的女人| 亚洲精品日韩av片在线观看| 中国美女看黄片| 真实男女啪啪啪动态图| 久久精品人妻少妇| 日日摸夜夜添夜夜爱| 精品久久久久久成人av| 1000部很黄的大片| 亚洲精品国产成人久久av| 亚洲成a人片在线一区二区| 国产精品99久久久久久久久| 久久婷婷人人爽人人干人人爱| 国产精品久久视频播放| 久久这里有精品视频免费| 午夜福利高清视频| 搡女人真爽免费视频火全软件| 国产一区二区在线观看日韩| 免费看a级黄色片| 国产精品.久久久| 少妇高潮的动态图| 亚洲经典国产精华液单| 久久久久久九九精品二区国产| 精品国产三级普通话版| 亚洲国产精品成人久久小说 | 亚洲av男天堂| 亚洲aⅴ乱码一区二区在线播放| 最近手机中文字幕大全| 亚洲av第一区精品v没综合| 欧美日韩乱码在线| 综合色丁香网| 嫩草影院入口| 97热精品久久久久久| 免费黄网站久久成人精品| 99视频精品全部免费 在线| 久久人人爽人人爽人人片va| 国模一区二区三区四区视频| 欧美色欧美亚洲另类二区| 中文字幕制服av| a级毛色黄片| 成年女人永久免费观看视频| 亚洲国产欧洲综合997久久,| 麻豆精品久久久久久蜜桃| 免费看光身美女| 国产伦精品一区二区三区视频9| 最近手机中文字幕大全| 国产伦理片在线播放av一区 | 亚洲欧美清纯卡通| 中文资源天堂在线| 久久久久久伊人网av| av在线播放精品| 国产中年淑女户外野战色| 亚洲人成网站在线播| 成人二区视频| 久久草成人影院| 中文字幕制服av| 99热这里只有精品一区| 可以在线观看毛片的网站| 成人亚洲欧美一区二区av| 在线观看一区二区三区| 深爱激情五月婷婷| 国产成人a∨麻豆精品| 亚洲五月天丁香| 1024手机看黄色片| 亚洲欧美中文字幕日韩二区| av在线老鸭窝| 99久久九九国产精品国产免费| 成人毛片60女人毛片免费| 日日撸夜夜添| 国产日本99.免费观看| 国内精品一区二区在线观看| 国产黄a三级三级三级人| 久久99热6这里只有精品| 内射极品少妇av片p| av天堂在线播放| 美女大奶头视频| 男女啪啪激烈高潮av片| 精品久久久久久久久亚洲| 日韩三级伦理在线观看| 免费大片18禁| 国产精品久久久久久精品电影小说 | 偷拍熟女少妇极品色| АⅤ资源中文在线天堂| 亚洲熟妇中文字幕五十中出| 欧美一区二区亚洲| 日韩 亚洲 欧美在线| 国产成人精品久久久久久| 国产精品三级大全| 国产成人精品婷婷| 国产精品一区二区三区四区久久| 午夜福利在线在线| 亚洲av二区三区四区| 波多野结衣巨乳人妻| 又黄又爽又刺激的免费视频.| 免费看日本二区| 中出人妻视频一区二区| av黄色大香蕉| 国产精品一区二区在线观看99 | 日韩高清综合在线| 国产欧美日韩精品一区二区| 国产蜜桃级精品一区二区三区| 国产麻豆成人av免费视频| www.av在线官网国产| 一边亲一边摸免费视频| 中出人妻视频一区二区| 此物有八面人人有两片| 午夜福利在线观看吧| 又黄又爽又刺激的免费视频.| 亚洲经典国产精华液单| 只有这里有精品99| 边亲边吃奶的免费视频| 欧美bdsm另类| 99九九线精品视频在线观看视频| 国产精品久久久久久精品电影小说 | 久久久成人免费电影| 嫩草影院新地址| 好男人视频免费观看在线| 99久国产av精品国产电影| 日韩 亚洲 欧美在线| 校园春色视频在线观看| www.av在线官网国产| 五月伊人婷婷丁香| 国产精品一及| 一级毛片我不卡| 波多野结衣巨乳人妻| 此物有八面人人有两片| 淫秽高清视频在线观看| 一区二区三区免费毛片| videossex国产| 亚洲va在线va天堂va国产| 亚洲一区高清亚洲精品| 亚洲成av人片在线播放无| 亚洲国产高清在线一区二区三| 午夜久久久久精精品| 青青草视频在线视频观看| 国产片特级美女逼逼视频| 亚洲第一电影网av| 欧美高清成人免费视频www| 亚洲成av人片在线播放无| 身体一侧抽搐| 床上黄色一级片| 1024手机看黄色片| 国产成人a区在线观看| 你懂的网址亚洲精品在线观看 | 好男人视频免费观看在线| 97超碰精品成人国产| 亚洲一级一片aⅴ在线观看| 国产视频内射| 人人妻人人澡人人爽人人夜夜 | 女人十人毛片免费观看3o分钟| 日本欧美国产在线视频| 欧美成人免费av一区二区三区| 午夜免费男女啪啪视频观看| 大型黄色视频在线免费观看| 亚洲欧美成人综合另类久久久 | 国产亚洲精品久久久com| 国产伦在线观看视频一区| 精品久久久久久久久亚洲| 国产一区二区在线av高清观看| a级毛色黄片| 韩国av在线不卡| 国产精品久久久久久久久免| 美女被艹到高潮喷水动态| 成人高潮视频无遮挡免费网站| 久久精品人妻少妇| 免费电影在线观看免费观看| 国产亚洲av嫩草精品影院| 欧美最黄视频在线播放免费| 国产私拍福利视频在线观看| 欧美高清性xxxxhd video| 久久欧美精品欧美久久欧美| 久久久国产成人免费| 国产精品1区2区在线观看.| 黄色视频,在线免费观看| 久久久久久久久久久丰满| 身体一侧抽搐| 国产精品日韩av在线免费观看| 国产高清视频在线观看网站| 观看美女的网站| 久久精品夜色国产| 天堂av国产一区二区熟女人妻| 人体艺术视频欧美日本| 国产真实伦视频高清在线观看| 91久久精品国产一区二区成人| 不卡视频在线观看欧美| 国产精品乱码一区二三区的特点| or卡值多少钱| 国产精品野战在线观看| 亚洲内射少妇av| 人妻久久中文字幕网| 欧美区成人在线视频| 美女国产视频在线观看| 一级毛片我不卡| 成熟少妇高潮喷水视频| 国产欧美日韩精品一区二区| 欧美日韩综合久久久久久| 九九久久精品国产亚洲av麻豆| 男女视频在线观看网站免费| 国产精品蜜桃在线观看 | 国产av麻豆久久久久久久| 欧美区成人在线视频| 国产精品人妻久久久影院| 精品日产1卡2卡| 神马国产精品三级电影在线观看| 熟女人妻精品中文字幕| 99久久久亚洲精品蜜臀av| 99riav亚洲国产免费| 国产精品一区二区三区四区免费观看| 日日干狠狠操夜夜爽| 在线播放无遮挡| 亚洲欧洲日产国产| 国内精品久久久久精免费| 国产私拍福利视频在线观看| 乱码一卡2卡4卡精品| 26uuu在线亚洲综合色| 99热这里只有是精品在线观看| ponron亚洲| 亚洲av免费高清在线观看| 国产精品一区二区三区四区免费观看| 精品午夜福利在线看| 麻豆成人午夜福利视频| 中文字幕久久专区| 亚洲国产高清在线一区二区三| av免费在线看不卡| 高清毛片免费观看视频网站| 亚洲av不卡在线观看| .国产精品久久| 99热只有精品国产| 欧美色视频一区免费| 日韩欧美在线乱码| 亚洲av.av天堂| 99热这里只有精品一区| 热99在线观看视频| 亚洲美女搞黄在线观看| 搡女人真爽免费视频火全软件| 日本黄色片子视频| 欧美日韩综合久久久久久| 成人漫画全彩无遮挡| 午夜精品国产一区二区电影| 中国三级夫妇交换| 插逼视频在线观看| 各种免费的搞黄视频| 丰满少妇做爰视频| tube8黄色片| 中文字幕人妻丝袜制服| 亚洲一级一片aⅴ在线观看| 曰老女人黄片| 色5月婷婷丁香| 亚洲av成人精品一二三区| 男人爽女人下面视频在线观看| 日本免费在线观看一区| 久久国内精品自在自线图片| 欧美日韩一区二区视频在线观看视频在线| 国产在视频线精品| 日本黄大片高清| 久久国内精品自在自线图片| 满18在线观看网站| 欧美日韩av久久| 三级国产精品欧美在线观看| 国产精品一区www在线观看| 成人毛片60女人毛片免费| av不卡在线播放| 插逼视频在线观看| 亚洲av免费高清在线观看| 精品国产露脸久久av麻豆| 91在线精品国自产拍蜜月| 人妻一区二区av| 日韩不卡一区二区三区视频在线| 久久国内精品自在自线图片| 两个人免费观看高清视频| 国产永久视频网站| 国产一区二区三区av在线| 亚洲欧洲精品一区二区精品久久久 | 在线看a的网站| 丰满少妇做爰视频| 亚洲欧美一区二区三区黑人 | 亚洲成人手机| 国产成人精品在线电影| av电影中文网址| 日本免费在线观看一区| 国产老妇伦熟女老妇高清| 男女免费视频国产| 狂野欧美白嫩少妇大欣赏| 视频在线观看一区二区三区| 国产精品国产三级国产专区5o| 日产精品乱码卡一卡2卡三| 国产精品一二三区在线看| 亚洲美女黄色视频免费看| 国产精品久久久久久精品古装| 国产一区二区三区av在线| 欧美亚洲日本最大视频资源| 亚洲欧美成人综合另类久久久| 亚洲图色成人| 制服诱惑二区| 人妻少妇偷人精品九色| 国产精品.久久久| 亚洲av福利一区| 久久久精品区二区三区| 另类亚洲欧美激情| av免费观看日本| 夜夜骑夜夜射夜夜干| 少妇的逼水好多| 少妇人妻 视频| 亚洲精品日本国产第一区| 免费大片黄手机在线观看| 多毛熟女@视频| 亚洲少妇的诱惑av| 伊人亚洲综合成人网| 能在线免费看毛片的网站| 日韩,欧美,国产一区二区三区| 亚洲精品av麻豆狂野| 人人妻人人澡人人爽人人夜夜| 一级毛片我不卡| 欧美xxⅹ黑人| 少妇人妻精品综合一区二区| 制服诱惑二区| 亚洲国产精品国产精品| av又黄又爽大尺度在线免费看| 精品一品国产午夜福利视频| 国产片内射在线| 麻豆成人av视频| av国产精品久久久久影院| 久久精品人人爽人人爽视色| 日本欧美国产在线视频| 国产黄色免费在线视频| 中文字幕亚洲精品专区| 国产精品国产三级国产av玫瑰| 精品久久久噜噜| 黑人猛操日本美女一级片| 日日摸夜夜添夜夜添av毛片| 久久国产亚洲av麻豆专区| 最近中文字幕2019免费版| av不卡在线播放| 日韩av不卡免费在线播放| av电影中文网址| 99视频精品全部免费 在线| xxxhd国产人妻xxx| 赤兔流量卡办理| 亚洲综合色惰| 久久人人爽人人片av| 免费av不卡在线播放| 日韩视频在线欧美| 亚洲欧美成人综合另类久久久| 麻豆精品久久久久久蜜桃| 亚洲精品aⅴ在线观看| 亚洲精品国产av成人精品| 国产成人精品在线电影| 黄色毛片三级朝国网站| 不卡视频在线观看欧美| 久久人人爽av亚洲精品天堂| 国产成人午夜福利电影在线观看| 九色成人免费人妻av| 夜夜爽夜夜爽视频| 免费不卡的大黄色大毛片视频在线观看| 久久精品国产鲁丝片午夜精品| 国产国语露脸激情在线看| 国产欧美日韩一区二区三区在线 | 日本欧美国产在线视频| 夫妻性生交免费视频一级片| 亚洲色图综合在线观看| 五月玫瑰六月丁香| 免费日韩欧美在线观看| 色94色欧美一区二区| 日韩成人av中文字幕在线观看| 新久久久久国产一级毛片| 最黄视频免费看| 国产成人aa在线观看| 欧美最新免费一区二区三区| 日韩成人av中文字幕在线观看| 18+在线观看网站| 夫妻午夜视频| 99re6热这里在线精品视频| 亚洲成人手机| 男的添女的下面高潮视频| 国产精品久久久久久精品电影小说| 免费观看av网站的网址| 亚洲高清免费不卡视频| 日韩熟女老妇一区二区性免费视频| 在线观看国产h片| 国产成人午夜福利电影在线观看| 国产片内射在线| 久久久精品免费免费高清| √禁漫天堂资源中文www| 亚洲性久久影院| 久久97久久精品| 中国美白少妇内射xxxbb| 亚洲国产精品一区三区| 久久久精品免费免费高清| 各种免费的搞黄视频| 亚洲av在线观看美女高潮| 午夜av观看不卡| 久久人人爽人人片av| 欧美xxⅹ黑人| 日韩成人av中文字幕在线观看| 久久狼人影院| 亚洲美女视频黄频| 观看av在线不卡| 在线免费观看不下载黄p国产| 午夜福利网站1000一区二区三区| 热re99久久精品国产66热6| 国产成人91sexporn| 成人午夜精彩视频在线观看| av在线播放精品| 国产男女超爽视频在线观看| 精品国产国语对白av| 国产免费现黄频在线看| 国产黄片视频在线免费观看| 欧美变态另类bdsm刘玥| 一级毛片电影观看| 亚洲人成网站在线观看播放| 日韩三级伦理在线观看| 亚洲精品乱码久久久v下载方式| 国产 一区精品| 建设人人有责人人尽责人人享有的| 在线看a的网站| 人妻人人澡人人爽人人| 国产av精品麻豆| 一区二区三区免费毛片| 久久国产精品大桥未久av| 亚洲欧美日韩卡通动漫| 尾随美女入室| 日韩成人av中文字幕在线观看| 97超碰精品成人国产| 国产精品一区二区在线不卡| 亚洲,欧美,日韩| 亚洲欧美清纯卡通| 欧美精品一区二区免费开放| 午夜影院在线不卡| xxxhd国产人妻xxx| 我的老师免费观看完整版| 飞空精品影院首页| 亚洲内射少妇av| 精品酒店卫生间| 男人爽女人下面视频在线观看| 欧美日韩精品成人综合77777| 中文字幕久久专区| 黄色毛片三级朝国网站| av.在线天堂| 国产精品一区www在线观看| 亚洲精品久久久久久婷婷小说| 久久99热6这里只有精品| 汤姆久久久久久久影院中文字幕| a 毛片基地| 久久免费观看电影| 好男人视频免费观看在线| 久久久久久久久久成人| 国产精品一区二区在线不卡| 少妇被粗大的猛进出69影院 | 国产又色又爽无遮挡免| 搡老乐熟女国产| 久久久欧美国产精品| 青春草亚洲视频在线观看| 欧美激情国产日韩精品一区| 99久国产av精品国产电影| 国产熟女午夜一区二区三区 | 视频中文字幕在线观看| 满18在线观看网站| 18禁观看日本| 国产亚洲最大av| 熟女电影av网| 欧美丝袜亚洲另类| 亚洲av福利一区| 黄色配什么色好看| 狂野欧美激情性bbbbbb| 色婷婷av一区二区三区视频| 亚洲无线观看免费| 五月天丁香电影| 国产国语露脸激情在线看| 欧美变态另类bdsm刘玥| 国语对白做爰xxxⅹ性视频网站| 国产午夜精品久久久久久一区二区三区| 国产成人精品一,二区| 久久精品久久精品一区二区三区| 亚洲av欧美aⅴ国产| 久久精品国产a三级三级三级| 国国产精品蜜臀av免费| 国产午夜精品一二区理论片| 最近中文字幕高清免费大全6| 岛国毛片在线播放| av国产久精品久网站免费入址| 亚洲中文av在线| 午夜福利视频精品| 色5月婷婷丁香| 日韩强制内射视频| 赤兔流量卡办理| 人人妻人人添人人爽欧美一区卜| 欧美激情国产日韩精品一区| 日日摸夜夜添夜夜爱| 国产男女超爽视频在线观看| 日韩免费高清中文字幕av| av.在线天堂| 97精品久久久久久久久久精品| 大香蕉久久成人网| 18禁动态无遮挡网站| 我的女老师完整版在线观看| 高清视频免费观看一区二区| 精品一区二区三区视频在线| 亚洲国产av影院在线观看| 伦精品一区二区三区| 美女国产高潮福利片在线看| 亚洲美女视频黄频| 高清视频免费观看一区二区| 人妻少妇偷人精品九色| 国产成人免费观看mmmm| 亚洲国产日韩一区二区| 水蜜桃什么品种好| 一级毛片 在线播放| 国产欧美亚洲国产| 一边摸一边做爽爽视频免费| 爱豆传媒免费全集在线观看| 国产精品一区二区三区四区免费观看| 久久精品国产亚洲网站| 精品久久蜜臀av无| 国产成人精品在线电影| 婷婷成人精品国产| 男女边吃奶边做爰视频| 欧美日韩一区二区视频在线观看视频在线|