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

    基于AHP分析法的磁浮系統(tǒng)車軌耦合振動抑制方法

    2020-12-07 06:55:44汪科任羅世輝陳曉昊馬衛(wèi)華鄒瑞明
    鐵道學報 2020年11期
    關鍵詞:電磁鐵阻尼軌道

    汪科任,羅世輝,陳曉昊,馬衛(wèi)華,鄒瑞明

    (西南交通大學 牽引動力國家重點實驗室,四川 成都 610031)

    磁浮列車不同于有軌列車,是一種新型的軌道交通運行方式,具有運行平穩(wěn)、乘坐舒適和無軌道接觸等優(yōu)點,具有較好的市場運用前景[1]。EMS型磁浮列車依靠電磁懸浮力實現(xiàn)穩(wěn)定懸浮,懸浮力將懸浮間隙與懸浮電流耦合在一起[2]。由于懸浮電磁力同時作用在磁浮列車與軌道梁上,若懸浮控制器性能稍差,極易出現(xiàn)車軌耦合振動。美國AMT磁浮列車就曾出現(xiàn)過在剛度較強的水泥梁上能夠?qū)崿F(xiàn)穩(wěn)定懸浮,而在剛度較低的橋梁上卻失穩(wěn)的現(xiàn)象[3]。針對該問題大量學者進行了相關研究,結(jié)果表明:剛度較低的彈性軌道梁是導致系統(tǒng)懸浮失穩(wěn)的主要原因。為了防止懸浮系統(tǒng)在彈性軌道梁上懸浮失穩(wěn),工程應用中常常通過增加軌道梁剛度與質(zhì)量來抑制車軌耦合振動,但從經(jīng)濟性的角度來考慮,增加軌道梁剛度與質(zhì)量將導致軌道梁建設成本的大幅增加。現(xiàn)階段在磁浮列車的整個造價中,軌道梁的建設成本高達60%~80%[4]。所以設計出性能更好的懸浮控制器來抑制懸浮系統(tǒng)在柔性軌道梁上的耦合振動,促進磁浮列車更好的商業(yè)化發(fā)展是非常必要的[5]。

    車軌耦合振動常常發(fā)生在磁浮列車靜懸浮或以較低速度通過彈性橋梁時,為抑制這種車軌耦合振動,國內(nèi)外學者進行了大量的相關研究[6-16]。其大致可以分為以下2種:一種是從控制算法的角度對懸浮控制器進行優(yōu)化設計,研究者認為懸浮控制器本身的特性是導致磁浮系統(tǒng)在彈性橋梁上失穩(wěn)的主要原因;另一種是采用某種特定的算法對磁浮列車本身的系統(tǒng)結(jié)構進行優(yōu)化設計,從而提高磁浮系統(tǒng)的動力學特性。文獻[6]利用Simulink軟件建立了HSST型五懸浮架磁浮列車垂向動力學模型,對比分析了不同懸浮控制算法對磁浮系統(tǒng)懸浮穩(wěn)定性的影響,并指出不同的懸浮控制算法將對懸浮系統(tǒng)的穩(wěn)定性產(chǎn)生不同的影響。文獻[7-8]以單鐵懸浮系統(tǒng)為研究對象,采用雙環(huán)PID控制算法,從能量轉(zhuǎn)換的角度解釋了系統(tǒng)發(fā)生車軌耦合的原因。文獻[9]利用狀態(tài)觀測器對懸浮電磁鐵的速度進行了準確估計,通過動力學仿真得出:當懸浮控制器頻率與軌道梁主頻率接近時,懸浮系統(tǒng)將發(fā)生車軌耦合振動現(xiàn)象。文獻[10]以單鐵懸浮系統(tǒng)為研究對象,指出在一定程度上提高懸浮控制器的電氣懸浮剛度可在一定程度上抑制系統(tǒng)的車軌耦合振動。文獻[11-13]各自采用不同的控制算法,分別研究了車輛運行速度、軌道梁長度和軌道梁抗彎剛度等對磁浮列車動態(tài)響應特性的影響,從而優(yōu)化了磁浮系統(tǒng)的動力學特性。文獻[14]指出:采用越多的狀態(tài)反饋,系統(tǒng)具有更優(yōu)的動力學特性。文獻[15]以單鐵懸浮系統(tǒng)為研究對象,將車體對懸浮電磁鐵的作用力視為額定的外界擾動,采用卡爾曼濾波實現(xiàn)了系統(tǒng)在彈性軌道梁上的穩(wěn)定懸浮。文獻[16]在文獻[15]的基礎上考慮了車體的振動,采用魯棒控制器實現(xiàn)了系統(tǒng)在柔性梁上的穩(wěn)定懸浮。然而,通過以上文獻分析可知,在懸浮控制器設計過程中,鮮有文獻對二次型加權矩陣權重系數(shù)的選取以及不同控制算法對磁浮系統(tǒng)車軌耦合振動的影響進行很好的機理說明。文獻[17]通過對單懸浮架模型與單鐵模型的仿真對比得出,對磁浮系統(tǒng)車軌耦合振動機理的探究,單鐵系統(tǒng)是適用的。因此本文以彈性軌道梁-懸浮電磁鐵-車體垂向耦合振動模型為研究對象,采用AHP法(層次分析法)程序化的求取了控制器二次型性能指標中狀態(tài)加權矩陣的權重系數(shù),并以此為基礎分析了控制器1與控制器2對磁浮系統(tǒng)懸浮穩(wěn)定性的影響。將采用軌道梁、懸浮電磁鐵以及車體的振動狀態(tài)進行主動控制的反饋控制器定義為反饋控制器1,將僅對懸浮電磁鐵振動狀態(tài)進行主動控制的傳統(tǒng)控制方式構成的反饋控制器定義為控制器2,如文獻[10,18-20]等采用的控制器。

    針對磁浮系統(tǒng)中二次型性能指標中狀態(tài)加權矩陣選取困難的問題,采用AHP法來程序化的求取該值,在此基礎上利用二次最優(yōu)控制理論設計出2種不同的狀態(tài)反饋控制器,在一定程度上從能量激振的角度分析了2種不同狀態(tài)反饋方式對懸浮系統(tǒng)車軌耦合振動響應特性的影響。

    1 車軌耦合等效模型

    針對西南交通大學常州軌道交通研究院某等效單鐵懸浮系統(tǒng)(見圖1),取向下為系統(tǒng)正方向,根據(jù)牛頓第二定律可得軌道梁、懸浮電磁鐵和車體在額定懸浮位置處的動力學方程為

    (1)

    (2)

    (3)

    F=-Ps(x2-x1)+PiIM

    (4)

    根據(jù)Meisinger推導可得到系統(tǒng)電學方程為[2]

    (5)

    圖1 簡化的單鐵-軌道梁-車體模型

    懸浮電磁鐵線圈兩端電壓與電流的數(shù)學表達式為

    (6)

    (7)

    由式(6)、式(7)可得

    (8)

    將式(8)在額定位置處Taylor展開,得

    (9)

    式中:Um為實際加載在懸浮電磁鐵線圈兩端的電壓;φ為流過懸浮電磁鐵極面積的磁通量;Lm為電磁鐵瞬時電感;cm為實際懸浮間隙;UN為平衡位置處電磁鐵線圈兩端的電壓;c為相對額定懸浮位置氣隙的變化量。

    (10)

    (11)

    由于c=x2-x1,由式(10)、式(11)可得式(5)。

    (12)

    其中,

    2 狀態(tài)反饋控制器設計

    2.1 狀態(tài)反饋增益矩陣K的求取

    磁浮列車在運行過程中常常出現(xiàn)懸浮控制器性能惡化現(xiàn)象,尤其在溫度較高的夏季特別突出。根據(jù)相關分析,其主要是由懸浮控制器發(fā)熱量較大引起的。所以在確保懸浮系統(tǒng)具有較好的動態(tài)響應特性的同時,如何限制懸浮控制器的發(fā)熱量是非常必要的。幸運的是該問題可以通過線性二次最優(yōu)控制理論來處理,通過該理論設計出的控制器能夠?qū)崿F(xiàn)系統(tǒng)控制能量與系統(tǒng)動態(tài)特性之間的平衡。定義被控系統(tǒng)(式12)二次型性能指標J為

    (13)

    PA+ATP-PBRu-1BTP+Q=0

    (14)

    式中:矩陣P為正定對稱矩陣。狀態(tài)反饋增益矩陣K為

    K=Ru-1BTP

    (15)

    2.2 狀態(tài)加權矩陣Q的求取

    式(13)中,與控制輸入矩陣Ru相比,狀態(tài)加權矩陣Q的維數(shù)較高,導致其加權系數(shù)的選取通常比較困難,一般都是通過經(jīng)驗試湊,并沒有統(tǒng)一的方法。本文采用AHP法來程序化的解決狀態(tài)加權矩陣Q求取困難的問題。AHP法是一種科學計算各元素權重系數(shù)的方法,是在對復雜決策問題的本質(zhì)、影響因素及其內(nèi)在關系等進行深入分析的基礎上將問題劃分為目標層、準則層和方案層,同時對各層次的元素進行兩兩比較而最終確定各元素加權系數(shù)的一種數(shù)學分析方法,能有效地應用于難以用定量分析解決權重系數(shù)取值的研究課題[21]。其設計步驟為[22]

    Step1建立控制指標的層次結(jié)構。選取系統(tǒng)二次型性能指標即式(13)作為目標層;反饋的狀態(tài)變量作為準則層;以滿足總體合理的一致性檢驗準則作為方案層。

    Step2構造線性系統(tǒng)判斷矩陣,用于確定兩兩元素的重要性比較。在判斷矩陣的構造中主觀加權比例系數(shù)的確定是關鍵,令Hij(無單位)為指標i與j重要性的比較值,其遵循以下準則:

    (1)2個元素兩兩比較,若同等重要用標度1表示;稍重要用標度3表示;明顯重要用標度5表示;強烈重要用標度7表示;極端重要用標度9表示。

    (2)若2個指標的相對重要性位于2個比較值之間,則取值分別為2、4、6、8。

    Step3計算每一層的權重并進行權重一致性的檢驗。若一致性比例系數(shù)ηR小于0.1,則認為系統(tǒng)通過一致性檢驗;否則需要重新規(guī)劃各狀態(tài)量之間的合理權重,直到其通過一致性檢驗為止。

    令判斷矩陣H(判斷矩陣反映了設計者對各加權系數(shù)權重的分配)最大特征值相應的特征向量為WH,由于列車在實現(xiàn)穩(wěn)定懸浮后,乘坐舒適性是列車運行品質(zhì)中一個極其重要的指標,所以相對于其余狀態(tài)變量的權重系數(shù),與車體振動相關狀態(tài)變量的權重系數(shù)應相對較大。由于反饋控制器1將系統(tǒng)的所有狀態(tài)變量都進行了主動控制,其涉及的狀態(tài)反饋量最多,設計也最為復雜,所以本文首先對反饋控制器1進行求解。

    根據(jù)以上步驟可構造判斷矩陣H為

    對判斷矩陣H每一列歸一化后求行和再進行列歸一化處理后可得判斷矩陣H最大特征值對應的特征向量WH,為

    WH=[0.124 00.140 80.155 30.043 90.136 2

    0.199 9 0.199 9]T

    (16)

    根據(jù)式(17)可求取矩陣WH最大特征值λmax為7.488 1。根據(jù)式(18)求得一致性指標CI為0.081 35。其中n1為判斷矩陣H的維數(shù)。由于判斷矩陣H為7維矩陣,可知隨機一致性指標RI為1.32[22],帶入式(19)可得一致性比率CR為0.062,其值小于0.1,可得判斷矩陣的不一致性程度在允許的范圍內(nèi),即通過一致性檢驗,則根據(jù)矩陣WH所確定的懸浮系統(tǒng)的二次型性能指標的權重系數(shù)是合理的,可用于式(13)中加權矩陣Q值得確定。

    (17)

    (18)

    (19)

    (20)

    由式(20)可知,求取控制器2時,式(13)中狀態(tài)加權矩陣的值有:Q2=diag(0.155 3,0.043 9,0.136 2),最后結(jié)合式(14)和式(15)可求取控制器2為:u2=-K2X2。同理,由式(16)可知,求取控制器1時,式(13)中狀態(tài)加權矩陣的值有:Q1=diag(0.124 0,0.140 8,0.155 3,0.043 9,0.136 2,0.199 9,0.199 9),最后結(jié)合式(14)和式(15)可求取控制器1為:u1=-K1X1。系統(tǒng)結(jié)構間見圖2。

    圖2 系統(tǒng)結(jié)構圖

    3 動力學仿真分析

    3.1 系統(tǒng)閉環(huán)特征根分析

    考查閉環(huán)系統(tǒng)穩(wěn)定性有多種方法,如李亞普羅夫穩(wěn)定性判據(jù)、勞斯判據(jù)等,而通過驗證閉環(huán)系統(tǒng)矩陣特征根實部是否均為負數(shù)最直觀,所以本文選取該方法來定性判斷狀態(tài)反饋控制器1、2對系統(tǒng)懸浮穩(wěn)定性的影響。耦合系統(tǒng)相關參數(shù)見表1。若文中無特殊說明均選取軌道梁跨中位置,即x1=l/2作為靜態(tài)懸浮位置。

    基于AHP分析法,在反饋控制器1與反饋控制器2控制下閉環(huán)系統(tǒng)矩陣A特征根實部最大值變化規(guī)律,見圖3。由圖3可以看出,當軌道梁阻尼一定時,采用控制器1的反饋控制方式,當軌道梁剛度從500 ~1 500 kN/m遞增的過程中,系統(tǒng)矩陣A特征根實部最大值均為負值,系統(tǒng)是穩(wěn)定的。而采用控制器2的反饋控制方式系統(tǒng)矩陣A特征根實部最大值均為正值,系統(tǒng)是不穩(wěn)定的。但進一步分析發(fā)現(xiàn),采用反饋控制器2的控制方式,隨著軌道梁剛度的增加,系統(tǒng)特征根實部的最大值是逐漸減小的,所以工程施工中通常采用大幅增加軌道梁剛度和阻尼的措施來抑制車軌耦合振動,但這將在一定程度上增加建設成本[15]。

    表1 系統(tǒng)參數(shù)

    圖3 特征根實部最大值

    3.2 數(shù)值仿真分析

    3.1節(jié)分析中對2種不同狀態(tài)反饋方式進行了相應的定性分析,得出當軌道梁阻尼一定時,反饋控制器2不能使系統(tǒng)在較低軌道梁剛度下(如600 kN/m)實現(xiàn)穩(wěn)定懸浮,而控制器1能夠?qū)崿F(xiàn)懸浮系統(tǒng)的穩(wěn)定懸浮。所以這節(jié)將定量分析采用AHP法時控制器1在軌道梁剛度為600 kN/m時懸浮間隙、軌道梁位移,懸浮電磁鐵位移與車體振動加速度的動態(tài)響應特性,分別見圖4~圖7。假設懸浮電磁鐵初始位置距離軌道梁2 mm。根據(jù)式(14)~式(16)、式(20)可得控制器1與控制器2反饋增益矩陣K1、K2為

    (21)

    (22)

    圖4 懸浮間隙

    圖5 軌道梁位移

    圖6 電磁鐵位移

    圖7 車體加速度

    由圖4~圖7分析可知,對軌道梁、懸浮電磁鐵以及車體的振動狀態(tài)進行控制,將其振動狀態(tài)引入狀態(tài)反饋控制器的控制方式,基于AHP分析法能夠?qū)崿F(xiàn)懸浮系統(tǒng)在較低軌道梁剛度下的穩(wěn)定懸浮,僅需約0.5 s便能實現(xiàn)系統(tǒng)的穩(wěn)定懸浮,車體最大振動加速度僅0.18 m/s2。

    3.3 能量激振分析

    本節(jié)將利用文獻[7]的方法,在一定程度上從懸浮系統(tǒng)對軌道梁能量輸出的角度來分析采用反饋控制器1與反饋控制器2對系統(tǒng)懸浮穩(wěn)定性的影響。由于磁浮列車-軌道梁耦合系統(tǒng)僅包括懸浮系統(tǒng)與軌道梁2個部分。根據(jù)對唐山磁浮線路的實際測量,軌道梁的振動特性可以用一定頻率的正弦函數(shù)或余弦函數(shù)來表示[23]。由于能量的大小可以通過功率來反映,而功率的大小為作用力與被作用物體的速度大小的乘積來表示。為方便從能量的角度來探究不同軌道梁頻率段車軌耦合振動的本質(zhì),這里不妨以軌道梁的速度作為輸入,作用在軌道梁上的懸浮力作為輸出。該傳遞函數(shù)可以表示為

    HFV(s)=F(s)/V(s)

    (23)

    在某時刻,以振動頻率為ω,振動幅值為0.2 m/s2,無相位差的正弦函數(shù)來模擬軌道梁的振動特性,即v(t)=0.2sin(ωt)??芍饔迷谲壍懒号c懸浮系統(tǒng)之間的電磁力為

    F(t)=0.2HFV(jω)|sin[ωt+∠HFV(jω)]

    (24)

    進一步分析可知,在一個時間周期內(nèi)電磁懸浮力作用在彈性軌道梁上功率大小為

    (25)

    式中:T為積分時間。采用表1中的相關參數(shù),式(21)、式(22)對應的控制器1與控制器2時,懸浮系統(tǒng)作用在軌道梁上的平均功率見圖8。

    圖8 無軌道阻尼時懸浮系統(tǒng)輸出功率

    由圖8可以看出,采用傳統(tǒng)的控制方式即僅將懸浮電磁鐵的振動狀態(tài)進行主動控制,當軌道梁振動頻率約低于50 rad/s時,懸浮系統(tǒng)向軌道梁輸出能量,此時更容易激起車軌耦合振動。只有當軌道梁振動頻率約高于50 rad/s時,懸浮系統(tǒng)才吸收軌道梁振動的能量,但對能量的吸收也相當少,最大僅有10 W左右。而采用將軌道梁、懸浮電磁鐵以及車體的振動狀態(tài)同時引入反饋控制系統(tǒng)進行主動控制的方式,在整個軌道梁振動頻率范圍內(nèi),懸浮系統(tǒng)始終不向軌道梁輸出能量。同時在軌道梁振動頻率為48 rad/s附近,懸浮系統(tǒng)對軌道梁的振動能量最大可吸收約為110 W。說明采用控制器1的反饋控制方式能夠較好的抑制系統(tǒng)在低頻段的車軌耦合振動。

    由于實際的軌道梁還存在一定的阻尼,如果懸浮電磁鐵向軌道梁輸出的功率大于軌道梁阻尼消耗的功率,那么必然導致整個系統(tǒng)集聚的能量越來越大,從而導致整個系統(tǒng)懸浮失穩(wěn)。忽略模態(tài)阻尼的非線性特性時,經(jīng)計算軌道梁阻尼對橋梁的吸收功率為[23]

    Pb=0.01ξωm1

    (26)

    式中:ξ為軌道梁阻尼比。當考慮軌道梁阻尼消耗的能量后,懸浮系統(tǒng)作用在軌道梁上的平均功率見圖9。

    圖9 有軌道阻尼時懸浮系統(tǒng)輸出功率

    由圖9分析可知,由于軌道梁阻尼消耗的功率十分有限,采用控制器2的反饋孔子方式,即使考慮軌道梁阻尼對軌道振動能量的消耗,當軌道梁振動頻率低于50 rad/s左右時,懸浮系統(tǒng)仍然對軌道梁輸出能量,最終必然導致懸浮系統(tǒng)失穩(wěn)。而采用控制器1的反饋控制方式,懸浮系統(tǒng)在整個軌道梁振動頻率段內(nèi)都不向軌道梁輸出能量,能保持系統(tǒng)較好的穩(wěn)定懸浮。由表1可知,本文軌道梁振動圓頻率約為16 rad/s,由圖9可知在該頻率段內(nèi),采用控制器2的反饋控制方式,懸浮控制器將向軌道梁輸出能量從而導致懸浮失穩(wěn)。

    4 結(jié)論

    針對彈性軌道梁-懸浮電磁鐵-車體垂向耦合振動模型,通過本文的仿真分析可得出以下結(jié)論:

    (1)基于AHP法,在一定程度上可程序化的求取系統(tǒng)二次型性能指標中狀態(tài)加權矩陣權重系數(shù)求解困難的問題,且能實現(xiàn)系統(tǒng)的穩(wěn)定懸浮。

    (2)在軌道梁阻尼一定時,當軌道梁剛度較低時(圓頻率較低時),僅對懸浮電磁鐵的振動狀態(tài)進行主動控制,懸浮系統(tǒng)可能向軌道梁輸出能量,從而引入系統(tǒng)出現(xiàn)車軌耦合振動。但隨著軌道梁剛度的增加,懸浮系統(tǒng)將吸收軌道梁的振動能量,從而抑制系統(tǒng)的車軌耦合振動,可能這也是工程中常常采用增加軌道梁剛度來抑制車軌耦合振動的原因。

    (3)在軌道梁阻尼一定時,采用將軌道梁、懸浮電磁鐵以及車體的振動狀態(tài)引入反饋控制器進行主動控制的方式,懸浮系統(tǒng)能吸收軌道梁輸出的能量,從而抑制車軌耦合振動。同時能夠?qū)崿F(xiàn)系統(tǒng)在低軌道梁剛度上的穩(wěn)定懸浮。

    猜你喜歡
    電磁鐵阻尼軌道
    奇妙的電磁鐵
    JB/T 10159—2019《交流本整濕式閥用電磁鐵》解讀
    N維不可壓無阻尼Oldroyd-B模型的最優(yōu)衰減
    關于具有阻尼項的擴散方程
    具有非線性阻尼的Navier-Stokes-Voigt方程的拉回吸引子
    磁懸浮小車電磁鐵優(yōu)化設計
    防爆電機(2021年2期)2021-06-09 08:14:38
    一種電磁鐵吸取減速頂調(diào)整墊的工具
    基于單純形法的TLE軌道確定
    CryoSat提升軌道高度與ICESat-2同步運行
    朝美重回“相互羞辱軌道”?
    亚洲专区字幕在线| 国产在线一区二区三区精| 看免费av毛片| 久久久久国产一级毛片高清牌| 久久天躁狠狠躁夜夜2o2o| 夫妻午夜视频| 国产淫语在线视频| 一二三四在线观看免费中文在| 一级片免费观看大全| 日韩精品免费视频一区二区三区| 日本黄色视频三级网站网址 | 精品视频人人做人人爽| 亚洲专区国产一区二区| 久久久国产一区二区| 亚洲一区二区三区欧美精品| 欧美黄色片欧美黄色片| 12—13女人毛片做爰片一| 国产成人一区二区三区免费视频网站| 久久中文看片网| 悠悠久久av| 搡老熟女国产l中国老女人| 99热网站在线观看| 99国产精品一区二区蜜桃av | 亚洲人成电影免费在线| а√天堂www在线а√下载 | 99热国产这里只有精品6| 一夜夜www| 亚洲久久久国产精品| 亚洲成人免费av在线播放| 男男h啪啪无遮挡| 视频区图区小说| 日韩大码丰满熟妇| 亚洲第一青青草原| 一二三四在线观看免费中文在| 三级毛片av免费| 成人免费观看视频高清| 无遮挡黄片免费观看| 欧美 日韩 精品 国产| 好看av亚洲va欧美ⅴa在| 久久影院123| 免费久久久久久久精品成人欧美视频| 大香蕉久久网| 精品人妻1区二区| 欧美性长视频在线观看| 国产精品久久电影中文字幕 | 制服诱惑二区| 18禁美女被吸乳视频| 19禁男女啪啪无遮挡网站| 日韩 欧美 亚洲 中文字幕| 三上悠亚av全集在线观看| 欧美精品啪啪一区二区三区| 精品第一国产精品| 亚洲人成伊人成综合网2020| 久久性视频一级片| 满18在线观看网站| 日本黄色日本黄色录像| 天天躁狠狠躁夜夜躁狠狠躁| 国产一区二区三区视频了| 又紧又爽又黄一区二区| 国产成人精品久久二区二区免费| 黑人猛操日本美女一级片| 精品久久久久久电影网| 99国产精品一区二区蜜桃av | 国产精品久久电影中文字幕 | 欧美精品av麻豆av| 99香蕉大伊视频| 午夜福利在线观看吧| 大香蕉久久网| 久久精品国产清高在天天线| 老熟女久久久| 亚洲精品国产区一区二| 久久精品国产综合久久久| 久久热在线av| 大码成人一级视频| aaaaa片日本免费| 在线观看日韩欧美| 一区在线观看完整版| 久久99一区二区三区| 国产欧美亚洲国产| 国产精品免费一区二区三区在线 | 男女床上黄色一级片免费看| 国产精品国产av在线观看| 免费看a级黄色片| 欧美精品亚洲一区二区| 国产成人精品久久二区二区免费| 精品一区二区三区av网在线观看| 国产区一区二久久| 成人黄色视频免费在线看| 一级,二级,三级黄色视频| 久久人妻福利社区极品人妻图片| 成人免费观看视频高清| 老汉色av国产亚洲站长工具| 两个人免费观看高清视频| 十八禁高潮呻吟视频| 黄网站色视频无遮挡免费观看| 老熟妇乱子伦视频在线观看| 91精品国产国语对白视频| 最近最新中文字幕大全电影3 | 亚洲性夜色夜夜综合| 大型av网站在线播放| 久久久久久久午夜电影 | 丁香欧美五月| 一进一出抽搐gif免费好疼 | 9191精品国产免费久久| 久久国产精品人妻蜜桃| 黄色丝袜av网址大全| av片东京热男人的天堂| 操出白浆在线播放| 久久午夜亚洲精品久久| 69av精品久久久久久| 国内毛片毛片毛片毛片毛片| 久久久水蜜桃国产精品网| 18禁国产床啪视频网站| 精品乱码久久久久久99久播| 国产欧美日韩精品亚洲av| www.999成人在线观看| 久久亚洲精品不卡| 国产成人系列免费观看| 婷婷精品国产亚洲av在线| 国产老妇女一区| 午夜日韩欧美国产| 国产精品 欧美亚洲| 亚洲欧美激情综合另类| 国产成人影院久久av| 日韩高清综合在线| 女人十人毛片免费观看3o分钟| 噜噜噜噜噜久久久久久91| 少妇的逼好多水| 久久亚洲精品不卡| 亚洲成人久久性| 日本黄色视频三级网站网址| 中文资源天堂在线| 国产激情欧美一区二区| 制服人妻中文乱码| 欧美日韩亚洲国产一区二区在线观看| 国产精品久久久久久久久免 | 最好的美女福利视频网| ponron亚洲| 日韩精品青青久久久久久| 久久国产乱子伦精品免费另类| 大型黄色视频在线免费观看| 两个人的视频大全免费| 欧美不卡视频在线免费观看| 亚洲国产精品合色在线| 老熟妇乱子伦视频在线观看| 成人永久免费在线观看视频| 一级毛片高清免费大全| 两个人看的免费小视频| 国产色爽女视频免费观看| 99久久综合精品五月天人人| 国产真实乱freesex| 欧美性猛交╳xxx乱大交人| 亚洲熟妇熟女久久| 亚洲精品美女久久久久99蜜臀| 岛国视频午夜一区免费看| 亚洲自拍偷在线| 中文字幕人成人乱码亚洲影| 国产精品久久久人人做人人爽| 亚洲片人在线观看| 一进一出抽搐gif免费好疼| 国产毛片a区久久久久| 又黄又粗又硬又大视频| 精品人妻一区二区三区麻豆 | 国产伦精品一区二区三区视频9 | 国产精品一区二区三区四区久久| 午夜a级毛片| 午夜福利在线在线| 91字幕亚洲| ponron亚洲| 国产高清视频在线播放一区| 国产99白浆流出| 欧美不卡视频在线免费观看| 免费人成在线观看视频色| 中文字幕久久专区| 欧洲精品卡2卡3卡4卡5卡区| 婷婷丁香在线五月| 99久久综合精品五月天人人| 亚洲一区高清亚洲精品| 国产精品久久久久久久久免 | 有码 亚洲区| 国产三级黄色录像| 免费观看人在逋| 桃红色精品国产亚洲av| 久久久久九九精品影院| 国产欧美日韩一区二区精品| 国产97色在线日韩免费| 熟女少妇亚洲综合色aaa.| 男女那种视频在线观看| 欧美又色又爽又黄视频| 亚洲精品亚洲一区二区| 小说图片视频综合网站| 日韩中文字幕欧美一区二区| 99国产极品粉嫩在线观看| 日韩欧美国产一区二区入口| 男人舔奶头视频| 日韩中文字幕欧美一区二区| 香蕉久久夜色| 国产精品亚洲美女久久久| 搡老熟女国产l中国老女人| 亚洲专区国产一区二区| 欧美日韩黄片免| 亚洲 国产 在线| 国产不卡一卡二| 最近视频中文字幕2019在线8| 有码 亚洲区| 国产一区二区亚洲精品在线观看| 91麻豆精品激情在线观看国产| 国产三级黄色录像| 免费看a级黄色片| 一本综合久久免费| 精品无人区乱码1区二区| 久久久久九九精品影院| 欧美一区二区国产精品久久精品| 一进一出好大好爽视频| 他把我摸到了高潮在线观看| 国产主播在线观看一区二区| 国产伦精品一区二区三区视频9 | 国产av不卡久久| 欧美一区二区亚洲| 久久天躁狠狠躁夜夜2o2o| 久久久国产成人免费| 桃色一区二区三区在线观看| 国产色爽女视频免费观看| 人人妻人人看人人澡| 免费无遮挡裸体视频| 国产爱豆传媒在线观看| 国产精品综合久久久久久久免费| 神马国产精品三级电影在线观看| 91麻豆精品激情在线观看国产| 哪里可以看免费的av片| 欧美日韩乱码在线| 毛片女人毛片| 少妇的丰满在线观看| 国产综合懂色| 亚洲成av人片在线播放无| 国产黄片美女视频| ponron亚洲| 婷婷六月久久综合丁香| 久久精品影院6| 男女视频在线观看网站免费| 操出白浆在线播放| 亚洲不卡免费看| 啦啦啦韩国在线观看视频| 伊人久久精品亚洲午夜| 全区人妻精品视频| 亚洲美女黄片视频| 宅男免费午夜| 成人特级av手机在线观看| a在线观看视频网站| 亚洲人成电影免费在线| 欧美3d第一页| 99热这里只有精品一区| 在线播放国产精品三级| 国产在线精品亚洲第一网站| 亚洲av中文字字幕乱码综合| av天堂中文字幕网| 精品电影一区二区在线| av国产免费在线观看| 搡女人真爽免费视频火全软件 | 午夜精品一区二区三区免费看| 人人妻人人看人人澡| 国语自产精品视频在线第100页| 麻豆成人av在线观看| 国产乱人视频| 精品久久久久久,| 一区福利在线观看| 国产一区二区激情短视频| 国产一区二区在线av高清观看| 我要搜黄色片| 免费看十八禁软件| 99riav亚洲国产免费| 国产成人影院久久av| 欧美乱妇无乱码| 亚洲欧美精品综合久久99| 国产又黄又爽又无遮挡在线| 激情在线观看视频在线高清| 亚洲av二区三区四区| 99久久久亚洲精品蜜臀av| 88av欧美| 18禁美女被吸乳视频| 成年人黄色毛片网站| 午夜精品久久久久久毛片777| xxxwww97欧美| 色综合婷婷激情| 在线观看av片永久免费下载| 免费av毛片视频| 免费人成视频x8x8入口观看| 免费看日本二区| 91在线观看av| 成人高潮视频无遮挡免费网站| 无限看片的www在线观看| 怎么达到女性高潮| 中文字幕高清在线视频| 国产亚洲精品一区二区www| 国产精品久久久久久亚洲av鲁大| 精品一区二区三区视频在线 | 亚洲av一区综合| 可以在线观看毛片的网站| www.色视频.com| 狂野欧美激情性xxxx| 日韩欧美国产一区二区入口| 国产精品嫩草影院av在线观看 | 观看免费一级毛片| 国产一区在线观看成人免费| 国产单亲对白刺激| 欧美av亚洲av综合av国产av| 首页视频小说图片口味搜索| 色在线成人网| 女警被强在线播放| 不卡一级毛片| 久久久色成人| 国产免费一级a男人的天堂| 久久天躁狠狠躁夜夜2o2o| 欧美精品啪啪一区二区三区| 俄罗斯特黄特色一大片| 国产v大片淫在线免费观看| 人妻丰满熟妇av一区二区三区| 少妇熟女aⅴ在线视频| 五月伊人婷婷丁香| 丰满乱子伦码专区| 久久久色成人| 欧美色欧美亚洲另类二区| 天堂影院成人在线观看| 黄色成人免费大全| 精品久久久久久,| 欧美在线黄色| 少妇的逼好多水| 99国产综合亚洲精品| 久久性视频一级片| 欧美中文综合在线视频| 成人欧美大片| 波野结衣二区三区在线 | 精华霜和精华液先用哪个| 中文字幕熟女人妻在线| 欧美成狂野欧美在线观看| 欧美另类亚洲清纯唯美| 少妇熟女aⅴ在线视频| 美女免费视频网站| 女人十人毛片免费观看3o分钟| 精品午夜福利视频在线观看一区| 国产免费男女视频| 日韩成人在线观看一区二区三区| 一个人免费在线观看的高清视频| 可以在线观看的亚洲视频| 免费在线观看日本一区| 欧美日韩黄片免| 久久久久久久久久黄片| 午夜免费激情av| 制服人妻中文乱码| 女同久久另类99精品国产91| 亚洲美女视频黄频| 两个人视频免费观看高清| 成年女人毛片免费观看观看9| 免费av观看视频| 婷婷精品国产亚洲av在线| 中文字幕熟女人妻在线| 亚洲人成伊人成综合网2020| 国产精品电影一区二区三区| 国产99白浆流出| 丰满的人妻完整版| 男女视频在线观看网站免费| 日本三级黄在线观看| h日本视频在线播放| 男人舔女人下体高潮全视频| 91久久精品国产一区二区成人 | 亚洲在线自拍视频| 午夜福利在线观看免费完整高清在 | 又紧又爽又黄一区二区| av欧美777| 国产精品 国内视频| 熟女少妇亚洲综合色aaa.| 一级毛片高清免费大全| 又紧又爽又黄一区二区| 国产一区在线观看成人免费| 免费高清视频大片| 国产精品久久视频播放| 亚洲成人久久爱视频| 国产综合懂色| 在线观看免费视频日本深夜| 亚洲激情在线av| 日韩av在线大香蕉| 亚洲天堂国产精品一区在线| 免费av观看视频| 一进一出抽搐gif免费好疼| 美女高潮喷水抽搐中文字幕| 欧美成人一区二区免费高清观看| 在线看三级毛片| 波多野结衣巨乳人妻| 欧美一区二区精品小视频在线| 色综合亚洲欧美另类图片| 男女视频在线观看网站免费| 亚洲av不卡在线观看| 久久久久免费精品人妻一区二区| 亚洲国产精品sss在线观看| 欧美一区二区亚洲| 国产精品久久久久久久久免 | 亚洲欧美日韩高清专用| 天堂网av新在线| 久久久久久人人人人人| 露出奶头的视频| 日韩有码中文字幕| 免费人成在线观看视频色| 在线观看av片永久免费下载| 岛国视频午夜一区免费看| 久久九九热精品免费| 99久久精品一区二区三区| 村上凉子中文字幕在线| 性色av乱码一区二区三区2| 黄色片一级片一级黄色片| 亚洲久久久久久中文字幕| 一区二区三区国产精品乱码| 大型黄色视频在线免费观看| 19禁男女啪啪无遮挡网站| 中文字幕熟女人妻在线| 国产精品嫩草影院av在线观看 | 精品国内亚洲2022精品成人| 国产视频一区二区在线看| 婷婷六月久久综合丁香| 亚洲黑人精品在线| 一级毛片女人18水好多| 国产乱人视频| 国产精品自产拍在线观看55亚洲| www国产在线视频色| 国产亚洲精品久久久久久毛片| 99在线人妻在线中文字幕| 在线免费观看不下载黄p国产 | 露出奶头的视频| 蜜桃亚洲精品一区二区三区| 一级黄色大片毛片| 有码 亚洲区| 人妻久久中文字幕网| 久久精品人妻少妇| 偷拍熟女少妇极品色| 黄色丝袜av网址大全| 身体一侧抽搐| 性色av乱码一区二区三区2| 尤物成人国产欧美一区二区三区| 三级男女做爰猛烈吃奶摸视频| 搡女人真爽免费视频火全软件 | 18+在线观看网站| 亚洲精品在线观看二区| 最近最新中文字幕大全免费视频| 尤物成人国产欧美一区二区三区| 国产精品久久久久久久电影 | 99久久综合精品五月天人人| 欧美日韩综合久久久久久 | 国内毛片毛片毛片毛片毛片| 日韩欧美免费精品| 成人18禁在线播放| 久久久久久人人人人人| 91在线精品国自产拍蜜月 | 在线视频色国产色| 在线观看日韩欧美| 亚洲色图av天堂| 噜噜噜噜噜久久久久久91| 久久人妻av系列| 精品久久久久久,| 国产淫片久久久久久久久 | 成年女人毛片免费观看观看9| 好男人电影高清在线观看| 亚洲黑人精品在线| 免费人成在线观看视频色| 欧美bdsm另类| 51国产日韩欧美| 国产亚洲精品久久久久久毛片| 国产伦人伦偷精品视频| 午夜精品一区二区三区免费看| 亚洲国产精品999在线| 在线观看日韩欧美| 婷婷精品国产亚洲av| 一个人免费在线观看的高清视频| 国内精品久久久久精免费| 久久精品国产亚洲av香蕉五月| 成年女人毛片免费观看观看9| 欧美日本亚洲视频在线播放| x7x7x7水蜜桃| 午夜福利高清视频| 亚洲最大成人手机在线| 国产毛片a区久久久久| 噜噜噜噜噜久久久久久91| 淫妇啪啪啪对白视频| 亚洲片人在线观看| 亚洲国产色片| 国产精品一区二区三区四区免费观看 | 亚洲国产欧美网| 黄片小视频在线播放| 午夜福利成人在线免费观看| 午夜福利高清视频| 一边摸一边抽搐一进一小说| 国产精品亚洲美女久久久| 国产av不卡久久| 精品久久久久久久久久久久久| 成人三级黄色视频| 非洲黑人性xxxx精品又粗又长| 国产视频一区二区在线看| 岛国在线免费视频观看| 少妇的逼水好多| 丰满人妻一区二区三区视频av | 亚洲中文日韩欧美视频| 国产精品一区二区三区四区久久| 国产精品久久电影中文字幕| 日本五十路高清| 一级作爱视频免费观看| 欧美最新免费一区二区三区 | eeuss影院久久| 无人区码免费观看不卡| 每晚都被弄得嗷嗷叫到高潮| 日本一二三区视频观看| 久久精品国产亚洲av涩爱 | 成熟少妇高潮喷水视频| 精品电影一区二区在线| 高清毛片免费观看视频网站| 亚洲欧美日韩东京热| 欧美日韩中文字幕国产精品一区二区三区| 国产av一区在线观看免费| 精品久久久久久,| 午夜激情福利司机影院| 成年女人看的毛片在线观看| 91久久精品国产一区二区成人 | 人人妻,人人澡人人爽秒播| 午夜福利视频1000在线观看| 国产精品99久久99久久久不卡| 精品久久久久久久毛片微露脸| 久久久久久久亚洲中文字幕 | 欧美bdsm另类| 狂野欧美激情性xxxx| 此物有八面人人有两片| bbb黄色大片| 国内少妇人妻偷人精品xxx网站| 亚洲av免费高清在线观看| 婷婷精品国产亚洲av| 成年女人永久免费观看视频| 国产成人系列免费观看| 国产精品嫩草影院av在线观看 | 亚洲欧美日韩无卡精品| 99精品久久久久人妻精品| 国产99白浆流出| 一区二区三区免费毛片| 一本综合久久免费| 亚洲av日韩精品久久久久久密| 亚洲av成人精品一区久久| 久久这里只有精品中国| 亚洲第一欧美日韩一区二区三区| 男人舔奶头视频| 操出白浆在线播放| 悠悠久久av| 久9热在线精品视频| 日本黄色视频三级网站网址| 欧美成人a在线观看| 色老头精品视频在线观看| 狠狠狠狠99中文字幕| 天天一区二区日本电影三级| 嫩草影院精品99| 18+在线观看网站| 欧美一区二区国产精品久久精品| 淫妇啪啪啪对白视频| 国产精品一及| 成人18禁在线播放| 男人舔女人下体高潮全视频| 特大巨黑吊av在线直播| 国产av在哪里看| 97超级碰碰碰精品色视频在线观看| 好男人电影高清在线观看| 精品无人区乱码1区二区| 日韩免费av在线播放| 99在线视频只有这里精品首页| 每晚都被弄得嗷嗷叫到高潮| 午夜精品久久久久久毛片777| 久久久久久人人人人人| 麻豆成人午夜福利视频| 身体一侧抽搐| 给我免费播放毛片高清在线观看| 久久香蕉精品热| 嫩草影院入口| 成年版毛片免费区| 母亲3免费完整高清在线观看| 国产精品98久久久久久宅男小说| 欧美在线一区亚洲| 成人国产综合亚洲| 欧美在线一区亚洲| 国语自产精品视频在线第100页| 丝袜美腿在线中文| 欧美激情久久久久久爽电影| 欧美zozozo另类| 国产精品99久久99久久久不卡| 99久久无色码亚洲精品果冻| 网址你懂的国产日韩在线| 日本 欧美在线| 女同久久另类99精品国产91| 久久久色成人| 嫩草影视91久久| 在线看三级毛片| 国产真人三级小视频在线观看| 特大巨黑吊av在线直播| 国产精品香港三级国产av潘金莲| 蜜桃久久精品国产亚洲av| 亚洲片人在线观看| 亚洲第一电影网av| 97超视频在线观看视频| 亚洲在线观看片| 一二三四社区在线视频社区8| 日韩大尺度精品在线看网址| 99在线人妻在线中文字幕| 国产精品野战在线观看| 乱人视频在线观看| 伊人久久精品亚洲午夜| 久久6这里有精品| 国产黄片美女视频| 高清毛片免费观看视频网站| 亚洲 欧美 日韩 在线 免费| 久久久国产成人免费| 国产乱人视频| 亚洲一区高清亚洲精品| 欧美在线黄色| 精品久久久久久久久久久久久|