• <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线在线观看网站| 亚洲成色77777| 婷婷色麻豆天堂久久| 啦啦啦在线观看免费高清www| 国产av国产精品国产| 秋霞在线观看毛片| 成人综合一区亚洲| 久久99一区二区三区| 七月丁香在线播放| 午夜福利视频精品| 欧美成人午夜免费资源| 亚洲精华国产精华液的使用体验| 十八禁网站网址无遮挡| 久久久久精品人妻al黑| 日韩,欧美,国产一区二区三区| a级毛片在线看网站| 中国三级夫妇交换| 中文字幕精品免费在线观看视频 | 国产成人精品久久久久久| 精品一区二区三区四区五区乱码 | 18禁国产床啪视频网站| 午夜日本视频在线| 国产片内射在线| 下体分泌物呈黄色| 国产淫语在线视频| 天天躁夜夜躁狠狠久久av| 中国三级夫妇交换| 校园人妻丝袜中文字幕| 亚洲精品日本国产第一区| 麻豆精品久久久久久蜜桃| 午夜免费观看性视频| 国产成人精品久久久久久| 亚洲精品乱久久久久久| a级毛片黄视频| 国产高清三级在线| 亚洲国产精品一区二区三区在线| 有码 亚洲区| 国产片特级美女逼逼视频| 又黄又爽又刺激的免费视频.| 欧美日韩综合久久久久久| 国产麻豆69| 美女内射精品一级片tv| tube8黄色片| 国产精品久久久久久久电影| 中文精品一卡2卡3卡4更新| 亚洲av免费高清在线观看| 男女午夜视频在线观看 | 下体分泌物呈黄色| 国产精品人妻久久久久久| 日韩av不卡免费在线播放| 国产男女内射视频| 一二三四在线观看免费中文在 | 亚洲av.av天堂| 精品国产国语对白av| 99久久人妻综合| 亚洲经典国产精华液单| 99热网站在线观看| 在线免费观看不下载黄p国产| 男女午夜视频在线观看 | 人妻系列 视频| 少妇的逼水好多| 男女边摸边吃奶| 日本vs欧美在线观看视频| 国产一区二区三区av在线| 亚洲成人一二三区av| 人妻少妇偷人精品九色| 热99国产精品久久久久久7| 9热在线视频观看99| 精品一区二区三卡| 亚洲精品乱码久久久久久按摩| 国产男人的电影天堂91| 亚洲精品成人av观看孕妇| 久久久欧美国产精品| 亚洲激情五月婷婷啪啪| 伦理电影免费视频| 天堂俺去俺来也www色官网| 午夜视频国产福利| 国产男女超爽视频在线观看| 欧美精品人与动牲交sv欧美| 日产精品乱码卡一卡2卡三| 我的女老师完整版在线观看| 国产日韩一区二区三区精品不卡| 国产免费视频播放在线视频| 日韩 亚洲 欧美在线| 国产深夜福利视频在线观看| 精品视频人人做人人爽| 不卡视频在线观看欧美| 边亲边吃奶的免费视频| av免费观看日本| 2018国产大陆天天弄谢| 久久精品久久精品一区二区三区| av不卡在线播放| 韩国av在线不卡| av在线老鸭窝| 少妇被粗大的猛进出69影院 | 国产福利在线免费观看视频| 男人操女人黄网站| 女人精品久久久久毛片| 国产综合精华液| 亚洲国产色片| tube8黄色片| 看十八女毛片水多多多| 国产精品人妻久久久影院| 国产片特级美女逼逼视频| 欧美变态另类bdsm刘玥| 国产 一区精品| 卡戴珊不雅视频在线播放| 亚洲av在线观看美女高潮| 一级毛片黄色毛片免费观看视频| 一级毛片电影观看| 在线亚洲精品国产二区图片欧美| 亚洲欧美一区二区三区黑人 | 中文乱码字字幕精品一区二区三区| 成年女人在线观看亚洲视频| 成人手机av| 在线观看美女被高潮喷水网站| 欧美 日韩 精品 国产| 亚洲经典国产精华液单| 中文精品一卡2卡3卡4更新| 国产 精品1| 精品一品国产午夜福利视频| 国产精品.久久久| 亚洲美女视频黄频| 午夜免费鲁丝| 女人被躁到高潮嗷嗷叫费观| 国产免费福利视频在线观看| 日韩三级伦理在线观看| 久久人人爽av亚洲精品天堂| 亚洲成色77777| 纵有疾风起免费观看全集完整版| 亚洲激情五月婷婷啪啪| 亚洲熟女精品中文字幕| 亚洲国产欧美在线一区| www日本在线高清视频| 精品亚洲成国产av| 成人午夜精彩视频在线观看| 看免费成人av毛片| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 91精品国产国语对白视频| 日日摸夜夜添夜夜爱| 一本大道久久a久久精品| 免费黄色在线免费观看| 纯流量卡能插随身wifi吗| 亚洲精品一二三| 边亲边吃奶的免费视频| 免费黄频网站在线观看国产| 久久国产精品男人的天堂亚洲 | 精品少妇久久久久久888优播| 男人舔女人的私密视频| 22中文网久久字幕| 在线观看美女被高潮喷水网站| 黑人巨大精品欧美一区二区蜜桃 | 成人二区视频| 大话2 男鬼变身卡| 日日爽夜夜爽网站| 男女无遮挡免费网站观看| 国产黄频视频在线观看| 毛片一级片免费看久久久久| av视频免费观看在线观看| 成人手机av| 18禁在线无遮挡免费观看视频| 在线观看一区二区三区激情| 这个男人来自地球电影免费观看 | 美女xxoo啪啪120秒动态图| 久久人人97超碰香蕉20202| 肉色欧美久久久久久久蜜桃| 精品一区在线观看国产| 亚洲伊人色综图| 高清黄色对白视频在线免费看| 制服人妻中文乱码| 精品亚洲成a人片在线观看| 美女大奶头黄色视频| 亚洲,一卡二卡三卡| 极品少妇高潮喷水抽搐| 欧美日韩成人在线一区二区| 丝袜在线中文字幕| 精品人妻熟女毛片av久久网站| 女人精品久久久久毛片| 97精品久久久久久久久久精品| 精品酒店卫生间| 五月玫瑰六月丁香| 久久久久久久大尺度免费视频| 国产一区二区三区综合在线观看 | 夜夜骑夜夜射夜夜干| 日韩免费高清中文字幕av| 国产精品一二三区在线看| 精品国产一区二区三区四区第35| 18+在线观看网站| 亚洲第一区二区三区不卡| av片东京热男人的天堂| 美女内射精品一级片tv| 自线自在国产av| 成人免费观看视频高清| 美女主播在线视频| 中文字幕av电影在线播放| 亚洲国产欧美日韩在线播放| 亚洲久久久国产精品| 国产精品成人在线| 国产精品一区二区在线不卡| 欧美成人午夜精品| 人人妻人人爽人人添夜夜欢视频| 制服诱惑二区| 免费观看无遮挡的男女| 国产在线一区二区三区精| 亚洲av欧美aⅴ国产| 久久 成人 亚洲| 亚洲国产看品久久| 亚洲美女视频黄频| 狂野欧美激情性bbbbbb| av.在线天堂| 9色porny在线观看| 午夜福利乱码中文字幕| 久久精品久久久久久噜噜老黄| 久久久久精品性色| 大陆偷拍与自拍| 制服人妻中文乱码| 亚洲国产毛片av蜜桃av| 日韩不卡一区二区三区视频在线| 国产乱来视频区| 国产av码专区亚洲av| 成人漫画全彩无遮挡| 精品亚洲成a人片在线观看| 国产精品成人在线| 久久精品夜色国产| 一级毛片黄色毛片免费观看视频| 视频中文字幕在线观看| 久久久精品免费免费高清| 国产麻豆69| 日韩一区二区三区影片| av网站免费在线观看视频| 中文字幕av电影在线播放| 考比视频在线观看| 亚洲综合色网址| 狠狠精品人妻久久久久久综合| 精品人妻偷拍中文字幕| 欧美人与性动交α欧美软件 | 丝袜美足系列| 欧美人与善性xxx| 亚洲综合色网址| 内地一区二区视频在线| 两性夫妻黄色片 | 99精国产麻豆久久婷婷| 黄色毛片三级朝国网站| 亚洲,欧美,日韩| 亚洲国产欧美在线一区| 国产亚洲精品第一综合不卡 | 少妇精品久久久久久久| 草草在线视频免费看| 日韩一区二区三区影片| 日韩av不卡免费在线播放| 精品久久久精品久久久| 久久久久久久久久久免费av| av在线播放精品| 80岁老熟妇乱子伦牲交| 国产亚洲欧美精品永久| 久久久a久久爽久久v久久| 午夜av观看不卡| 久久久久精品性色| 亚洲欧美精品自产自拍| 美国免费a级毛片| 不卡视频在线观看欧美| 黄片播放在线免费| 日韩一区二区视频免费看| 精品一区二区三区四区五区乱码 | 日本黄大片高清| 日日摸夜夜添夜夜爱| h视频一区二区三区| 亚洲精品自拍成人| 在线观看美女被高潮喷水网站| 国产精品一二三区在线看| 少妇被粗大的猛进出69影院 | 狂野欧美激情性xxxx在线观看| 黄网站色视频无遮挡免费观看| 亚洲精品美女久久av网站| 日本wwww免费看| 99久国产av精品国产电影| 成年动漫av网址| 国产精品无大码| 国产成人免费无遮挡视频| 精品亚洲成国产av| 欧美精品一区二区大全| 亚洲av电影在线观看一区二区三区| 晚上一个人看的免费电影| 亚洲人成77777在线视频| 自拍欧美九色日韩亚洲蝌蚪91| 99热全是精品| 国产成人精品在线电影| 国产精品三级大全| 999精品在线视频| 一级毛片黄色毛片免费观看视频| 国产精品无大码| 日本wwww免费看| 国产精品一区www在线观看| 性高湖久久久久久久久免费观看| 人人妻人人澡人人爽人人夜夜| 一级a做视频免费观看| 亚洲欧美一区二区三区国产| 看免费成人av毛片| 爱豆传媒免费全集在线观看| 三级国产精品片| 亚洲,欧美,日韩| 国产精品国产三级国产av玫瑰| 婷婷成人精品国产| 大话2 男鬼变身卡| 成人影院久久| 国产精品一区www在线观看| 国产精品一区www在线观看| 国产精品.久久久| 丰满乱子伦码专区| 超碰97精品在线观看| 久久久久国产网址| 一本色道久久久久久精品综合| 丝袜脚勾引网站| 久久精品国产亚洲av涩爱| 女人精品久久久久毛片| 侵犯人妻中文字幕一二三四区| 精品人妻熟女毛片av久久网站| 99久久中文字幕三级久久日本| 日韩 亚洲 欧美在线| 九草在线视频观看| 中文字幕另类日韩欧美亚洲嫩草| 一本色道久久久久久精品综合| 99热6这里只有精品| 国产精品一区二区在线观看99| 亚洲综合色惰| 国产 精品1| 国产成人精品在线电影| 18禁在线无遮挡免费观看视频| 大香蕉久久网| 啦啦啦视频在线资源免费观看| 高清视频免费观看一区二区| 美女xxoo啪啪120秒动态图| 免费看不卡的av| 黑人猛操日本美女一级片| 国产淫语在线视频| 免费av中文字幕在线| 精品一品国产午夜福利视频| 免费av中文字幕在线| 成年动漫av网址| 久久人人爽人人爽人人片va| 青春草视频在线免费观看| 久久久久视频综合| 夜夜爽夜夜爽视频| 国产欧美另类精品又又久久亚洲欧美| 欧美精品高潮呻吟av久久| 最近最新中文字幕免费大全7| 韩国精品一区二区三区 | 久久99热6这里只有精品| 久久久亚洲精品成人影院| 9热在线视频观看99| 永久免费av网站大全| 老司机影院成人| 满18在线观看网站| 欧美精品亚洲一区二区| 国产黄色免费在线视频| 啦啦啦啦在线视频资源| 波多野结衣一区麻豆| 婷婷色av中文字幕| 少妇被粗大猛烈的视频| 内地一区二区视频在线| 七月丁香在线播放| 欧美97在线视频| 大片电影免费在线观看免费| 少妇的逼水好多| 狂野欧美激情性xxxx在线观看| kizo精华| 日本av手机在线免费观看| 深夜精品福利| 卡戴珊不雅视频在线播放| videosex国产| av电影中文网址| 免费观看av网站的网址| 国产精品久久久久久精品电影小说| 成人毛片a级毛片在线播放| 欧美日韩av久久| 国产精品 国内视频| 97精品久久久久久久久久精品| 国产黄色免费在线视频| 亚洲av综合色区一区| 国产 一区精品| 女的被弄到高潮叫床怎么办| 菩萨蛮人人尽说江南好唐韦庄| 高清在线视频一区二区三区| 黑人高潮一二区| 日本欧美视频一区| 国产男人的电影天堂91| 久久毛片免费看一区二区三区| 亚洲精品国产色婷婷电影| 国产精品不卡视频一区二区| xxxhd国产人妻xxx| 国产男女超爽视频在线观看| 国内精品宾馆在线| 母亲3免费完整高清在线观看 | 久久久久久人妻| 亚洲美女视频黄频| 一级爰片在线观看| 国产精品欧美亚洲77777| 少妇被粗大的猛进出69影院 | 亚洲精品久久成人aⅴ小说| 精品久久国产蜜桃| 人妻一区二区av| 观看美女的网站| www.av在线官网国产| 一二三四中文在线观看免费高清| 久久久久久久久久人人人人人人| www日本在线高清视频| 伊人亚洲综合成人网| 秋霞在线观看毛片| 少妇被粗大猛烈的视频| 乱码一卡2卡4卡精品| 久久久久视频综合| 午夜免费男女啪啪视频观看| 黄色怎么调成土黄色| 国产成人精品福利久久| 99九九在线精品视频| 狠狠婷婷综合久久久久久88av| 2021少妇久久久久久久久久久| 午夜福利网站1000一区二区三区| 亚洲五月色婷婷综合| 欧美精品一区二区免费开放| 99久久人妻综合| 黑人高潮一二区| 日韩一区二区视频免费看| 成人18禁高潮啪啪吃奶动态图| 久久久久精品人妻al黑| 日韩欧美精品免费久久| 成人毛片60女人毛片免费| 熟女人妻精品中文字幕| 日本av免费视频播放| 久久人人97超碰香蕉20202| 最近最新中文字幕免费大全7| 男人添女人高潮全过程视频| 两个人看的免费小视频| 国产精品一区二区在线观看99| 国产亚洲精品久久久com| 国产永久视频网站| 国产精品久久久久久精品古装| av在线观看视频网站免费| 老司机影院成人| 国产福利在线免费观看视频| 99热这里只有是精品在线观看| 精品人妻偷拍中文字幕| 人妻少妇偷人精品九色| 国产精品国产三级专区第一集| 欧美 日韩 精品 国产| 亚洲综合色惰| 欧美变态另类bdsm刘玥| 高清不卡的av网站| 亚洲欧美一区二区三区黑人 | 深夜精品福利| 久久精品国产综合久久久 | 美女视频免费永久观看网站| 精品一区二区三区四区五区乱码 | 成人国语在线视频| 免费看不卡的av| 国精品久久久久久国模美| 国产熟女欧美一区二区| 最近的中文字幕免费完整| 国产精品无大码| 伊人久久国产一区二区| 有码 亚洲区| 国产精品.久久久| 国产男女内射视频| av在线老鸭窝| 极品人妻少妇av视频| a级毛色黄片| 亚洲av欧美aⅴ国产| 最后的刺客免费高清国语| 性高湖久久久久久久久免费观看| 男女国产视频网站| 亚洲第一av免费看| 极品人妻少妇av视频| 在线亚洲精品国产二区图片欧美| 91aial.com中文字幕在线观看| 国产精品国产三级国产专区5o| 永久免费av网站大全| 大香蕉97超碰在线| 欧美日韩成人在线一区二区| av线在线观看网站| 黑人猛操日本美女一级片| 成年人免费黄色播放视频| 一区二区三区乱码不卡18| 国产精品 国内视频| 飞空精品影院首页| xxx大片免费视频| 尾随美女入室| 一区二区三区乱码不卡18| 久久久亚洲精品成人影院| 久久久久久人人人人人| 中文精品一卡2卡3卡4更新| 国产毛片在线视频| 午夜福利影视在线免费观看| 国产精品国产三级专区第一集| 熟妇人妻不卡中文字幕| 色哟哟·www| 十八禁网站网址无遮挡| av不卡在线播放| 久久99蜜桃精品久久| 亚洲图色成人| 色哟哟·www| www.av在线官网国产| 色婷婷久久久亚洲欧美| 亚洲熟女精品中文字幕| 天天影视国产精品| 免费大片18禁| 精品熟女少妇av免费看| 丰满饥渴人妻一区二区三| 亚洲综合色惰| 人妻系列 视频| 国产一区二区在线观看av| 日韩欧美一区视频在线观看| 日日摸夜夜添夜夜爱| 国产成人精品无人区| 免费看av在线观看网站| av又黄又爽大尺度在线免费看| 熟女av电影| 好男人视频免费观看在线| 欧美精品人与动牲交sv欧美| 狠狠精品人妻久久久久久综合| 高清av免费在线| 不卡视频在线观看欧美| 在线 av 中文字幕| 99香蕉大伊视频| 精品国产一区二区三区四区第35| 亚洲 欧美一区二区三区| 寂寞人妻少妇视频99o| 国产精品久久久久久精品古装| 黑人高潮一二区| 91精品三级在线观看| 欧美国产精品va在线观看不卡| 日本av免费视频播放| 午夜影院在线不卡| 最近中文字幕高清免费大全6| 满18在线观看网站| 国产欧美日韩一区二区三区在线| 亚洲内射少妇av| 伦理电影大哥的女人| 国产精品99久久99久久久不卡 | av在线观看视频网站免费| 亚洲少妇的诱惑av| 国产精品一区二区在线不卡| 你懂的网址亚洲精品在线观看| 免费看av在线观看网站| 自拍欧美九色日韩亚洲蝌蚪91| 日韩,欧美,国产一区二区三区| 曰老女人黄片| 亚洲av成人精品一二三区| 久久毛片免费看一区二区三区| 成人无遮挡网站| 欧美亚洲 丝袜 人妻 在线| 午夜福利视频在线观看免费| 久久久久久久精品精品| 成人毛片a级毛片在线播放| 内地一区二区视频在线| 中文字幕最新亚洲高清| 肉色欧美久久久久久久蜜桃| 国产成人欧美| 色哟哟·www| 欧美最新免费一区二区三区| 99热这里只有是精品在线观看| 日韩三级伦理在线观看| 欧美日本中文国产一区发布| 欧美日韩av久久| 国产一区二区在线观看av| 亚洲,欧美精品.| 一本大道久久a久久精品| av免费观看日本| 久久这里有精品视频免费| videossex国产| 黄片无遮挡物在线观看| 亚洲三级黄色毛片| 精品卡一卡二卡四卡免费| 2022亚洲国产成人精品| 一个人免费看片子| 国产成人av激情在线播放| 天天操日日干夜夜撸| 老熟女久久久| 亚洲国产精品一区三区| 满18在线观看网站| 老熟女久久久| 国产成人av激情在线播放| 交换朋友夫妻互换小说| 亚洲伊人色综图| 视频中文字幕在线观看| 22中文网久久字幕| 国产亚洲av片在线观看秒播厂| 精品少妇久久久久久888优播| 国产精品 国内视频| 精品亚洲成a人片在线观看| 精品少妇久久久久久888优播| 天天躁夜夜躁狠狠久久av| 乱码一卡2卡4卡精品| av在线播放精品| 午夜日本视频在线| 亚洲第一av免费看| 嫩草影院入口| 午夜av观看不卡| 80岁老熟妇乱子伦牲交| 中文字幕精品免费在线观看视频 | 日韩电影二区| 精品国产国语对白av| 亚洲熟女精品中文字幕| 两个人免费观看高清视频| 一本大道久久a久久精品| 久久久久精品久久久久真实原创| 久久亚洲国产成人精品v| 国内精品宾馆在线| 日本wwww免费看| 十八禁网站网址无遮挡| 久久久久国产精品人妻一区二区| 色吧在线观看| 精品酒店卫生间| 国产亚洲一区二区精品| 国产成人精品婷婷| 美女脱内裤让男人舔精品视频|