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

    懸浮控制參數(shù)對(duì)磁浮車輛與道岔梁耦合振動(dòng)的影響分析

    2022-03-11 03:32:02楊志南馮洋劉東生趙春發(fā)
    機(jī)械 2022年2期
    關(guān)鍵詞:控制參數(shù)電磁鐵磁懸浮

    楊志南,馮洋,劉東生,趙春發(fā)*,

    懸浮控制參數(shù)對(duì)磁浮車輛與道岔梁耦合振動(dòng)的影響分析

    楊志南1,馮洋2,劉東生2,趙春發(fā)*,2

    (1.中鐵第一勘察設(shè)計(jì)院集團(tuán)有限公司,陜西 西安 710043;2.西南交通大學(xué) 牽引動(dòng)力國(guó)家重點(diǎn)實(shí)驗(yàn)室,四川 成都 610031)

    建立了包含PID懸浮控制系統(tǒng)的兩節(jié)編組磁浮列車動(dòng)力學(xué)模型和細(xì)致的道岔梁有限元模型,仿真分析了不同懸浮控制參數(shù)條件下磁浮列車以30 km/h通過(guò)道岔梁時(shí)的動(dòng)力學(xué)響應(yīng)。仿真結(jié)果表明,間隙反饋系數(shù)取值偏小時(shí),電磁鐵懸浮間隙的波動(dòng)幅值較大,不利于行車安全;增大間隙反饋系數(shù),電磁懸浮系統(tǒng)的特征頻率逐漸接近道岔梁的一階垂彎模態(tài)頻率,導(dǎo)致車岔耦合振動(dòng)更強(qiáng)烈。間隙速度反饋系數(shù)取值較小時(shí),電磁懸浮系統(tǒng)的阻尼偏小,車體振動(dòng)加速度較大;間隙速度反饋系數(shù)取值偏大時(shí),微分環(huán)節(jié)超前調(diào)節(jié)作用過(guò)強(qiáng),車岔耦合振動(dòng)劇烈。為緩解磁浮車輛與道岔梁的耦合振動(dòng),建議間隙反饋系數(shù)在7000~8000之間選取,速度反饋系數(shù)在45~60之間選取。

    磁浮列車;道岔;反饋控制;耦合振動(dòng);數(shù)值仿真

    常導(dǎo)電磁懸浮型(EMS,Electromagnetic Suspension)磁浮列車采用電磁吸力支承和導(dǎo)向車輛,以直線感應(yīng)電機(jī)驅(qū)動(dòng)列車前進(jìn),具有振動(dòng)小、噪聲低、乘坐舒適及爬坡能力強(qiáng)等優(yōu)點(diǎn),近年來(lái)在國(guó)內(nèi)外得到快速發(fā)展與應(yīng)用。目前,我國(guó)已開(kāi)通運(yùn)營(yíng)上海高速磁浮線、長(zhǎng)沙磁浮快線和北京地鐵S1磁浮線,還有正在建設(shè)的清遠(yuǎn)磁浮旅游線和鳳凰磁浮旅游線。EMS磁浮交通的技術(shù)優(yōu)勢(shì)在上述工程應(yīng)用中已得到證實(shí),但同時(shí)也有一些工程問(wèn)題沒(méi)有得到很好的解決。例如,國(guó)內(nèi)外EMS磁浮列車研發(fā)和試運(yùn)營(yíng)期間,都曾出現(xiàn)過(guò)車輛與鋼梁強(qiáng)烈耦合振動(dòng)的問(wèn)題[1-9],嚴(yán)重時(shí)甚至導(dǎo)致懸浮失穩(wěn),出現(xiàn)電磁鐵砸軌現(xiàn)象。

    為減緩磁浮車輛與鋼梁的耦合振動(dòng),工程上采用了提高鋼梁剛度及結(jié)構(gòu)阻尼或安裝外部吸能裝置等措施,如上海高速磁浮線對(duì)維修基地鋼框架梁進(jìn)行加固[3],在鋼道岔梁上加裝調(diào)諧質(zhì)量減振器(TMD,Tuned Mass Damper)[4]。長(zhǎng)沙磁浮快線在道岔主動(dòng)梁中間增設(shè)支承臺(tái)車,并且在梁上安裝液體質(zhì)量雙調(diào)諧減振器(TLMD,Tuned Liquid Mass Damper)[6-9]。這些被動(dòng)減振措施增加了磁浮道岔的制造安裝成本,不利于EMS磁浮交通的推廣應(yīng)用,因此,很有必要從優(yōu)化懸浮控制的角度探索減緩車岔耦合振動(dòng)的方法。事實(shí)上,已有研究表明全狀態(tài)反饋?zhàn)顑?yōu)控制、瞬時(shí)最優(yōu)控制、自適應(yīng)控制等能有效抑制磁浮車軌耦合振動(dòng)[10-13],但這些控制算法相對(duì)復(fù)雜,工程實(shí)施困難,且參數(shù)調(diào)節(jié)不方便,故工程實(shí)際仍采用經(jīng)典的PID控制算法。

    對(duì)于采用PID懸浮控制的EMS磁浮列車,已有試驗(yàn)研究對(duì)不同工況下車輛和道岔梁的振動(dòng)響應(yīng)進(jìn)行了細(xì)致分析[6-9],但極少研究控制參數(shù)對(duì)系統(tǒng)動(dòng)力學(xué)響應(yīng)的影響,相關(guān)的數(shù)值模擬研究也很少見(jiàn),因此,如何調(diào)整懸浮控制參數(shù)來(lái)減緩車岔耦合振動(dòng)還缺乏成熟的理論指導(dǎo)。為此,本文結(jié)合國(guó)內(nèi)某中低速磁浮線工程,建立兩節(jié)編組磁浮列車多體動(dòng)力學(xué)模型、PID懸浮控制模型和三段鉸接式磁浮道岔有限元模型,仿真分析低速運(yùn)行條件下磁浮車輛與道岔梁的耦合振動(dòng)響應(yīng),重點(diǎn)研究懸浮控制器比例系數(shù)、微分系數(shù)對(duì)車岔耦合振動(dòng)的影響,從避免磁浮車岔耦合振動(dòng)過(guò)大的角度提出懸浮控制參數(shù)的取值建議。

    1 車岔耦合振動(dòng)分析模型

    1.1 車輛動(dòng)力學(xué)模型

    中低速磁浮車輛由車體、二系懸掛系統(tǒng)、迫導(dǎo)向機(jī)構(gòu)、懸浮架、懸浮電磁鐵、直線電機(jī)定子件以及各種輔助裝備組成。通常每車配備五個(gè)懸浮架,懸浮架由側(cè)梁、托臂、抗側(cè)滾裝置等組成。懸浮電磁鐵與懸浮架托臂、托臂與側(cè)梁之間通過(guò)螺栓連接,直線電機(jī)定子件安裝于側(cè)梁底面,組成具有懸浮、導(dǎo)向和牽引功能的懸浮模塊??箓?cè)滾裝置由左右兩片抗側(cè)滾梁和吊桿組成,在抑制懸浮模塊側(cè)滾的同時(shí),允許左右模塊之間發(fā)生小范圍的相對(duì)運(yùn)動(dòng)。車輛左右側(cè)第2、5位空簧上部的承臺(tái)與車體底部固結(jié),其他位空簧上部通過(guò)可橫向移動(dòng)的滑臺(tái)與車體相連。在第1、3位滑臺(tái)和第4、6位滑臺(tái)之間還設(shè)置了迫導(dǎo)向機(jī)構(gòu)。迫導(dǎo)向機(jī)構(gòu)由T型臂、橫向拉桿和鋼纜組成,用以提高磁浮車輛曲線通過(guò)能力。

    在磁浮車輛結(jié)構(gòu)分析的基礎(chǔ)上,考慮到單節(jié)磁浮車輛的長(zhǎng)度約為15 m,而磁浮道岔主動(dòng)梁的長(zhǎng)度約為19 m,因此,在滿足道岔主動(dòng)梁可滿跨受力的條件下,為提高模型計(jì)算效率,使用SIMPACK軟件建立了兩節(jié)編組磁浮列車動(dòng)力學(xué)模型。單節(jié)車輛的動(dòng)力學(xué)模型如圖1所示,其中將車體、懸浮模塊、抗側(cè)滾梁、可移動(dòng)滑臺(tái)、T型臂、橫向拉桿、牽引拉桿等視為剛體,將二系空氣彈簧簡(jiǎn)化為線性彈簧阻尼器。車體考慮除縱向運(yùn)動(dòng)之外的5個(gè)自由度,懸浮模塊考慮6個(gè)運(yùn)動(dòng)自由度,抗側(cè)滾梁具有獨(dú)立的搖頭自由度,移動(dòng)滑臺(tái)相對(duì)車體僅具有橫向平移自由度,每車共有141個(gè)剛體運(yùn)動(dòng)自由度。中低速磁浮車輛模型的主要參數(shù)如表1所示。

    圖1 中低速磁浮車輛模型

    1.2 懸浮控制模型

    中低速磁浮車輛每個(gè)懸浮電磁鐵包含四個(gè)線圈繞組,分為前后兩組,分別受控于前、后懸浮控制器,即每個(gè)懸浮電磁鐵包含兩個(gè)相對(duì)獨(dú)立的懸浮控制點(diǎn),因此,可針對(duì)單點(diǎn)懸浮系統(tǒng)(圖2)進(jìn)行懸浮控制系統(tǒng)的設(shè)計(jì)與建模分析。

    表1 中低速磁浮車輛模型主要參數(shù)

    u為電磁鐵線圈電壓,V;i為電磁鐵線圈電流,A;δ為懸浮間隙,mm。

    運(yùn)用牛頓第二定律,可推導(dǎo)出單點(diǎn)懸浮系統(tǒng)的垂向動(dòng)力學(xué)方程為:

    根據(jù)電磁感應(yīng)定律,可得到電磁鐵線圈的電壓和電流關(guān)系式為:

    對(duì)于式(1)、式(2)組成的三階電磁懸浮系統(tǒng),一般采用如圖3所示的雙環(huán)控制器進(jìn)行懸浮控制[14-15]。

    δ0為額定間隙,mm;kc為電流環(huán)比例系數(shù);i0為初始電流,A;s為拉氏變換因子。

    圖3中,電流環(huán)采用比例或比例-積分反饋控制,目的是使線圈電流以最快速度跟蹤控制電壓,避免電流延時(shí)引起懸浮失穩(wěn);位置環(huán)可設(shè)計(jì)為懸浮間隙、間隙速度和電磁鐵加速度的反饋控制,主要目的是實(shí)現(xiàn)恒定的懸浮間隙。

    對(duì)于式(3)所示被控系統(tǒng),引入電流負(fù)反饋和比例控制,可得電流環(huán)的閉環(huán)傳遞函數(shù)為:

    對(duì)比式(3)、式(4)可知,選取合適的比例系數(shù)后,電流環(huán)的時(shí)間常數(shù)將遠(yuǎn)小于開(kāi)環(huán)電壓-電流系統(tǒng)的時(shí)間常數(shù),線圈電流可快速跟蹤控制電壓,使得懸浮系統(tǒng)演變?yōu)槎A系統(tǒng),從而可以將電流視為位置環(huán)的控制量。

    將式(1)中的電磁力計(jì)算項(xiàng)在平衡點(diǎn)(0,0)附近線性化,得到單點(diǎn)懸浮系統(tǒng)的線性化模型為:

    對(duì)線性化系統(tǒng)采用PID控制算法,同時(shí)附加電磁鐵加速度反饋以壓縮系統(tǒng)頻帶,則電磁鐵線圈電流的控制律為:

    將式(6)代入式(5),忽略消除系統(tǒng)穩(wěn)態(tài)誤差的間隙積分項(xiàng)和壓縮系統(tǒng)頻帶的加速度反饋?lái)?xiàng),可得到平衡點(diǎn)附近閉環(huán)懸浮系統(tǒng)的動(dòng)力學(xué)方程為:

    由式(7)可得到:

    式中:為主動(dòng)控制磁懸浮系統(tǒng)的特征頻率,Hz;ξ為阻尼比。

    由式(8)可知,如果增大間隙反饋系數(shù),但間隙微分反饋系數(shù)保持不變,則增大,ξ減小;若增大間隙微分反饋系數(shù),但間隙反饋系數(shù)保持不變,則不變,ξ增大。

    按照上述雙環(huán)控制器設(shè)計(jì)原理,使用MATLAB軟件建立懸浮控制模型,并利用SIMPACK軟件的MatSIM接口導(dǎo)入建立的磁浮車輛動(dòng)力學(xué)模型。

    1.3 道岔梁有限元模型

    中低速磁浮交通三段鉸接式道岔主要由主動(dòng)梁、第一從動(dòng)梁、第二從動(dòng)梁、垛梁、驅(qū)動(dòng)裝置、鎖定裝置以及電氣系統(tǒng)等組成[16]。道岔主體結(jié)構(gòu)圖4所示。

    1.活動(dòng)端垛梁;2.主動(dòng)梁;3.第二從動(dòng)梁;4.第一從動(dòng)梁;5.固定端垛梁。

    圖4中,主動(dòng)梁、第一和第二從動(dòng)梁的長(zhǎng)度分別為19.8 m、4.8 m和5 m,主動(dòng)梁中間處梁高1.4 m,隔板間距1.2 m,腹板厚度24 mm。使用ANSYS軟件建立道岔有限元模型如圖5所示,其中支座處采用彈性約束,梁體選用SHELL181殼單元建模,網(wǎng)格尺寸為100 mm,單元數(shù)量約4.1萬(wàn),詳細(xì)建模過(guò)程參見(jiàn)文獻(xiàn)[17]。

    圖5 磁浮道岔梁有限元模型

    有限元模態(tài)分析表明,磁浮道岔的前六階模態(tài)頻率分別為8.60 Hz、15.77 Hz、18.99 Hz、20.90 Hz、25.47 Hz和31.99 Hz。文獻(xiàn)[18]中實(shí)測(cè)獲得的道岔主動(dòng)梁第一階橫向彎曲自振頻率為8.5 Hz,第一階垂向自振頻率均值為15.7 Hz,與有限元模型前兩模態(tài)頻率的誤差分別為1%和0.4%,多個(gè)實(shí)測(cè)高階特征頻率與計(jì)算結(jié)果也較接近,說(shuō)明建立的道岔有限元模型符合實(shí)際。

    在后文的仿真計(jì)算中,為了模擬磁浮列車進(jìn)出道岔的動(dòng)力學(xué)響應(yīng),在圖5模型的前后端均建立了2跨跨度25 m的混凝土軌道梁模型。此外,基于德國(guó)高速低干擾軌道譜生成磁浮軌道不平順樣本[18],將其作為動(dòng)力學(xué)仿真計(jì)算時(shí)的線路不平順激擾。

    1.4 耦合振動(dòng)計(jì)算方法

    利用作者所在課題組開(kāi)發(fā)的磁浮交通系統(tǒng)分布式協(xié)同仿真平臺(tái)實(shí)現(xiàn)SIMPACK車輛模型和ANSYS有限元模型的協(xié)同仿真,仿真平臺(tái)包括兩臺(tái)工作站、SIMPACK軟件、ANSYS軟件以及自主開(kāi)發(fā)的接口程序和耦合計(jì)算器[19]。車輛模塊和道岔模塊通過(guò)接口程序與耦合計(jì)算器進(jìn)行數(shù)據(jù)交互和耦合計(jì)算,耦合器控制不同計(jì)算機(jī)上的任務(wù)模塊輸入輸出數(shù)據(jù)的收發(fā)和仿真時(shí)序,并通過(guò)網(wǎng)絡(luò)通知相應(yīng)的執(zhí)行機(jī)(SIMPACK計(jì)算器和ANSYS計(jì)算器)。

    2 車岔耦合振動(dòng)響應(yīng)特性

    磁浮車輛與道岔梁強(qiáng)烈耦合振動(dòng)多發(fā)生在車輛低速通過(guò)時(shí),速度較高時(shí)耦合振動(dòng)并不明顯,因此,本節(jié)計(jì)算了車輛以30 km/h速度通過(guò)道岔時(shí)的動(dòng)力學(xué)響應(yīng)。選取了兩組懸浮控制參數(shù),A組參數(shù)為:=9000,=70,=500,對(duì)應(yīng)的=10.8 Hz,ξ=0.04。B組參數(shù)為:=13000,取值與工況A相同,變?yōu)?4.0 Hz,ξ=0.03。

    2.1 振動(dòng)響應(yīng)時(shí)域分析

    圖6分別給出了采用A、B組控制參數(shù)時(shí)頭車車體質(zhì)心垂向加速度、頭車右側(cè)第5位電磁鐵前端的懸浮間隙和懸浮力的時(shí)程曲線。

    圖6 以30 km/h速度過(guò)岔時(shí)車輛動(dòng)力學(xué)響應(yīng)

    由圖6(a)可知,A組和B組控制參數(shù)對(duì)應(yīng)的車體質(zhì)心加速度均較小,最大值分別為0.04 m/s2和0.14 m/s2,后者較前者增加250%。由圖6(b)可知,列車進(jìn)入道岔主動(dòng)梁以后,兩組控制參數(shù)條件下懸浮間隙的調(diào)整頻率和幅度均明顯增大,間隙波動(dòng)幅值分別為0.51 mm和0.81 mm,小于其允許波動(dòng)限值4 mm。由圖6(c)可知,列車進(jìn)出道岔時(shí)電磁鐵懸浮力有明顯的突變,這是因?yàn)榛炷淋壍懒汉偷啦砹褐g存在較大的豎向剛度差;同樣,車輛進(jìn)出道岔主動(dòng)梁時(shí),懸浮力的波動(dòng)較為劇烈,這是因?yàn)閺膭?dòng)梁的豎向剛度明顯大于主動(dòng)梁,主動(dòng)梁滿跨受載時(shí)其梁端轉(zhuǎn)角也較大。其中,采用B組控制參數(shù)時(shí),主動(dòng)梁上單個(gè)線圈最大懸浮力達(dá)到11.32 kN,較其額定懸浮力8 kN增大41.5%。

    圖7是磁浮車輛以30 km/h速度通過(guò)時(shí)主動(dòng)梁跨中的垂向位移和垂向加速度時(shí)程響應(yīng)。計(jì)算結(jié)果表明,采用A組控制參數(shù)時(shí),主動(dòng)梁跨中垂向位移最大值為1.92 mm,垂向加速度最大值為6.5 m/s2;采用B組參數(shù)時(shí),跨中垂向位移最大值為2.7 mm,垂向加速度最大值達(dá)到13.6 m/s2;采用B組參數(shù)時(shí),主動(dòng)梁跨中垂向位移和垂向加速度幅值較A組參數(shù)分別增加40.6%和109.2%。由圖7還可知,當(dāng)間隙反饋系數(shù)由9000增大到13000以后,列車駛出主動(dòng)梁時(shí)的振動(dòng)響應(yīng)更加劇烈,這與圖6(c)中懸浮力的響應(yīng)規(guī)律一致。

    綜合分析圖6和圖7可知,兩組控制參數(shù)條件下,磁浮車輛和道岔梁均發(fā)生了明顯的耦合振動(dòng)。特別當(dāng)懸浮間隙反饋系數(shù)增大到13000以后,由于電磁懸浮系統(tǒng)的特征頻率更接近道岔梁的豎向基頻,磁軌動(dòng)力作用更為劇烈,故車輛系統(tǒng)和道岔主動(dòng)梁振動(dòng)響應(yīng)更強(qiáng)烈。

    2.2 振動(dòng)響應(yīng)頻譜分析

    圖8是頭車右側(cè)第5位電磁鐵前端懸浮間隙響應(yīng)的頻譜圖??梢钥闯觯瑑山M控制參數(shù)條件下懸浮間隙響應(yīng)的優(yōu)勢(shì)頻率均在5 Hz以內(nèi),它們與車輛二系懸掛頻率和懸浮架固有頻率有關(guān);由于懸浮間隙響應(yīng)包含了道岔梁的垂向位移響應(yīng),主動(dòng)梁的第一、二階垂彎模態(tài)頻率在頻譜圖中也有明顯表現(xiàn)。

    圖7 以30 km/h速度過(guò)岔時(shí)主動(dòng)梁動(dòng)力學(xué)響應(yīng)

    圖8 第5位電磁鐵前端懸浮間隙響應(yīng)頻譜

    圖9給出了頭車右側(cè)第5位電磁鐵前端的懸浮力響應(yīng)頻譜圖。采用A組控制參數(shù)時(shí),電磁懸浮力的第一、二階主頻分別為31.23 Hz和15.19 Hz;采用B組控制參數(shù)時(shí),懸浮力的第一、二階主頻分別為15.67 Hz和31.23 Hz??梢?jiàn),懸浮力的響應(yīng)主頻接近于道岔主動(dòng)梁的前兩階垂彎模態(tài)頻率。進(jìn)一步對(duì)主動(dòng)梁跨中的垂向加速度響應(yīng)進(jìn)行頻譜分析,發(fā)現(xiàn)其主頻分布與圖9相似,僅幅值略有差別,不再給出。

    圖9 第5位電磁鐵前端懸浮力響應(yīng)頻譜

    上述頻譜分析表明,A組和B組控制參數(shù)條件下,電磁懸浮力的調(diào)整頻率總是與主動(dòng)梁的低階垂彎自振頻率接近。對(duì)于B組控制參數(shù),由于平衡點(diǎn)附近電磁懸浮系統(tǒng)的特征頻率更為接近主動(dòng)梁的第一階垂彎模態(tài)頻率,導(dǎo)致磁浮車輛與道岔梁的耦合振動(dòng)更為劇烈。

    3 控制參數(shù)對(duì)耦合振動(dòng)的影響

    從改變電磁懸浮系統(tǒng)特征頻率(剛度)和阻尼的角度出發(fā),開(kāi)展間隙反饋系數(shù)、間隙微分反饋系數(shù)對(duì)磁浮車岔耦合振動(dòng)的影響分析,其中k取值范圍為5000~15000,k取值范圍為20~80[20]。

    3.1 間隙反饋系數(shù)影響分析

    為探究k對(duì)磁浮車岔耦合振動(dòng)的影響,將k設(shè)定為70、k設(shè)定為500,k分別取5000、7000、9000、11000和13000。相應(yīng)的,平衡點(diǎn)附近,在4.7~14 Hz之間變化,ξ在0.40~0.79之間變化。

    圖10顯示了頭車和尾車車體質(zhì)心垂向加速度隨k的變化曲線??芍^、尾車車體質(zhì)心的垂向加速度均隨k的增加而增大;當(dāng)k>9000后,車體質(zhì)心加速度幅值隨k的增長(zhǎng)速度變大,說(shuō)明增大電磁懸浮的剛度會(huì)導(dǎo)致車體垂向振動(dòng)加速度增大,從而降低列車運(yùn)行平穩(wěn)性。

    圖11給出了頭車和尾車右側(cè)第5位電磁鐵前端的懸浮間隙波動(dòng)量隨k的變化曲線??梢园l(fā)現(xiàn),懸浮間隙波動(dòng)幅值總體上隨k的增大先快速減小后趨于穩(wěn)定;當(dāng)k>7000以后,懸浮間隙波動(dòng)幅值差別不大??梢?jiàn),k取值較大時(shí),電磁懸浮系統(tǒng)的剛度較大,從而降低了懸浮間隙的波動(dòng)幅值,有利于預(yù)防磁軌發(fā)生碰撞。

    圖10 車體垂向加速度與kp的關(guān)系

    圖11 懸浮間隙波動(dòng)幅值與kp的關(guān)系

    圖12給出了頭、尾車右側(cè)第5位電磁鐵前端的懸浮力方差隨k的變化曲線。結(jié)果表明,k越大,電磁懸浮力的方差越大,這與圖10中車體質(zhì)心加速度的變化規(guī)律一致,說(shuō)明增大k將加劇車輛與道岔梁之間的動(dòng)力作用水平。

    圖12 懸浮力方差與kp的關(guān)系

    圖13是主動(dòng)梁1/4跨、跨中和3/4跨處梁體的垂向加速度幅值與k的關(guān)系曲線。可以發(fā)現(xiàn),不同位置處梁體的垂向加速度總體上均隨k的增大而增大;當(dāng)k>7000后,1/4跨和3/4跨處主動(dòng)梁的振動(dòng)加速度明顯大于跨中處,說(shuō)明列車進(jìn)出主動(dòng)梁時(shí)的車岔耦合振動(dòng)更劇烈。

    圖13 主動(dòng)梁垂向加速度與kp的關(guān)系

    綜合分析圖10~13可知,單獨(dú)增大懸浮控制器的間隙反饋系數(shù),相當(dāng)于增大了電磁懸浮的剛度,電磁懸浮系統(tǒng)的特征頻率隨之增大,并逐漸接近道岔主動(dòng)梁的一階垂彎模態(tài)頻率,此時(shí),盡管電磁鐵懸浮間隙波動(dòng)量不會(huì)明顯增大,但車輛與道岔梁的振動(dòng)加速度、電磁懸浮力方差均明顯增大。

    3.2 間隙微分反饋系數(shù)影響分析

    k設(shè)定為9000,k設(shè)定為500,k分別取30、45、60和70。相應(yīng)的,平衡點(diǎn)附近,=10.8 Hz,ξ在0.17~0.4之間變化。

    圖14給出了頭、尾車車體質(zhì)心垂向加速度與k的關(guān)系曲線??梢钥闯觯?i>k增大后,ξ增大,車體質(zhì)心加速度隨k的增大呈現(xiàn)先減小后緩慢增大的趨勢(shì),但車體質(zhì)心垂向加速度幅值均很小,總體上中等大小的k有利于降低磁浮車輛的振動(dòng)響應(yīng)。

    圖14 車體質(zhì)心垂向加速度與kd的關(guān)系

    圖15是頭、尾車右側(cè)第5位電磁鐵前端的懸浮間隙波動(dòng)量與k的關(guān)系曲線,可以看出,電磁鐵懸浮間隙波動(dòng)量總體上隨k的增大而變化不大。圖16給出了頭、尾車右側(cè)第5位電磁鐵前端的懸浮力方差與k的關(guān)系曲線,可以看出,懸浮力方差隨k的增大先減小后增大,k在45~60范圍內(nèi)取值時(shí)磁浮車輛和道岔梁的動(dòng)力作用水平較低。

    圖15 懸浮間隙波動(dòng)幅值與kd的關(guān)系

    圖16 懸浮力方差與kd的關(guān)系

    圖17給出了主動(dòng)梁1/4跨、跨中和3/4跨處垂向加速度與k的關(guān)系曲線??梢园l(fā)現(xiàn),不同位置處梁體的垂向加速度均隨k的增大先緩慢減小后快速增大;當(dāng)k>60后,盡管電磁懸浮系統(tǒng)的阻尼更大,但較大的k加劇了系統(tǒng)的振蕩,使得車岔耦合振動(dòng)更為劇烈。

    圖17 主動(dòng)梁垂向加速度與kd的關(guān)系

    由圖14~17可知,隨著k的逐步增大,由于電磁懸浮系統(tǒng)的阻尼隨之增大,車岔耦合振動(dòng)響應(yīng)先呈現(xiàn)為減小趨勢(shì);但是,當(dāng)進(jìn)一步增大k時(shí),超前調(diào)節(jié)作用過(guò)強(qiáng),使得系統(tǒng)動(dòng)態(tài)性能惡化;因此,過(guò)小或過(guò)大的k均有可能導(dǎo)致車岔耦合振動(dòng)強(qiáng)烈。

    4 結(jié)論

    結(jié)合國(guó)內(nèi)某中低速磁浮線工程,仿真分析了磁浮列車以30 km/h速度通過(guò)道岔梁時(shí)的動(dòng)力學(xué)響應(yīng),得到以下結(jié)論:

    (1)對(duì)比分析了間隙反饋系數(shù)為9000和13000時(shí)的車岔耦合振動(dòng)響應(yīng),由于后者的懸浮剛度更大,電磁懸浮系統(tǒng)的特征頻率與道岔主動(dòng)梁的一階垂彎模態(tài)頻率接近,故磁浮車岔耦合振動(dòng)更為劇烈。道岔主動(dòng)梁梁端轉(zhuǎn)角、主動(dòng)梁與混凝土梁之間的豎向剛度差導(dǎo)致磁浮列車出岔時(shí)車岔耦合振動(dòng)明顯增大,其中主動(dòng)梁振動(dòng)加速度增長(zhǎng)最為顯著。電磁懸浮力、道岔主動(dòng)梁振動(dòng)加速度的響應(yīng)頻率以主動(dòng)梁的前兩階豎彎模態(tài)為主,而懸浮間隙的響應(yīng)頻率則以5 Hz以下的低頻為主。

    (2)當(dāng)間隙微分反饋系數(shù)不變、間隙反饋系數(shù)由5000逐漸增加到13000時(shí),車體豎向加速度、電磁力方差、主動(dòng)梁豎向加速度的響應(yīng)幅值均隨間隙反饋系數(shù)的增大而增大,但懸浮間隙波動(dòng)量先快速減小后平緩變化。間隙反饋系數(shù)在7000~8000范圍內(nèi)取值時(shí),磁浮車岔耦合振動(dòng)響應(yīng)較小。

    (3)當(dāng)間隙反饋系數(shù)不變、間隙微分反饋系數(shù)由30逐漸增加到70時(shí),車體垂向加速度和懸浮間隙波動(dòng)量隨間隙反饋系數(shù)的增大先明顯減小后平緩變化,懸浮力方差和主動(dòng)梁加速度先減小后增大。間隙微分反饋系數(shù)在45~60之間取值時(shí),磁浮車岔耦合振動(dòng)響應(yīng)較小。

    [1]翟婉明,趙春發(fā). 現(xiàn)代軌道交通工程科技前沿與挑戰(zhàn)[J]. 西南交通大學(xué)學(xué)報(bào),2016,51(2):209-226.

    [2]Fichtner K,Pichlmeier F. The Transrapid guideway switch test and verification[C]. Shanghai:Proceedings of the 18th International Conference on Magnetically Levitated Systems and Linear Drives,2004:624-631.

    [3]李莉,孟光. 慢起慢落時(shí)磁浮車輛與鋼軌道框架耦合共振分析[J]. 振動(dòng)與沖擊,2006,25(6):46-48.

    [4]殷月俊,羅漢中,黃醒春. 高速磁浮道岔振動(dòng)響應(yīng)的原位實(shí)測(cè)[J]. 上海交通大學(xué)學(xué)報(bào),2007,41(4):658-663.

    [5]Zhao Chunfa,Gu Xingtao,Xiao Zhou. Dynamic analysis of maglev train and switch beam coupled system[C]. Chengdu:Advance in Environmental Vibration,the 5th International Symposium on Environmental Vibration,2011:528-533.

    [6]靖仕元. 多重調(diào)諧質(zhì)量調(diào)諧阻尼器的磁浮道岔減振方案[J]. 鐵道工程學(xué)報(bào),2019,36(4):80-83,89.

    [7]柴小鵬,汪正興,王波,等. 磁浮工程道岔梁的TLMD減振技術(shù)研究[J]. 世界橋梁,2017,45(2):60-65.

    [8]劉大玲. 中低速磁浮道岔動(dòng)載試驗(yàn)方法研究[J]. 城市軌道交通研究,2017,20(5):89-92.

    [9]羅華軍,吳志會(huì),佟來(lái)生,等. 中低速磁浮交通車岔耦合振動(dòng)研究[J]. 電力機(jī)車與城軌車輛,2018,41(1):5-8.

    [10]王輝,鐘曉波,沈鋼. 彈性軌道梁上磁懸浮控制方法[J]. 交通運(yùn)輸工程學(xué)報(bào),2013,13(5):33-38,46.

    [11]黎松奇,張昆侖,劉國(guó)清,等. 基于逆系統(tǒng)方法的磁浮列車非線性控制[J]. 控制工程,2017,24(8):1542-1546.

    [12]李曉龍,翟明達(dá),郝阿明. 基于輸出飽和條件的磁浮列車懸浮控制參數(shù)優(yōu)化[J]. 國(guó)防科技大學(xué)學(xué)報(bào),2017,39(4):149-153.

    [13]周丹峰,李杰,余佩倡,等. 磁浮交通軌排耦合自激振動(dòng)分析及自適應(yīng)控制方法[J]. 自動(dòng)化學(xué)報(bào),2019,45(12):2328-2343.

    [14]李云鋼,常文森. 磁浮列車懸浮系統(tǒng)的串級(jí)控制[J]. 自動(dòng)化學(xué)報(bào),1999,25(2):107-111.

    [15]劉恒坤,常文森. 磁懸浮列車的雙環(huán)控制[J]. 控制工程,2007,14(2):198-2000.

    [16]曾國(guó)鋒,袁亦竑,吉文,等. 長(zhǎng)沙中低速磁浮工程的道岔設(shè)計(jì)與調(diào)試[J]. 城市軌道交通研究,2016,19(5):44-48.

    [17]張宇生,趙春發(fā),周文,等. 不同約束條件下中低速磁浮道岔主動(dòng)梁自振特性[J]. 鐵道標(biāo)準(zhǔn)設(shè)計(jì),2020,64(4):22-27.

    [18]任曉博,趙春發(fā),馮洋,等. 中低速磁浮車輛-軌道-橋梁垂向耦合振動(dòng)仿真分析[J]. 鐵道標(biāo)準(zhǔn)設(shè)計(jì),2019,63(2):70-76.

    [19]李倩,黃海于,馮洋,等. 磁浮交通系統(tǒng)動(dòng)力學(xué)分布式協(xié)同仿真接口的設(shè)計(jì)與實(shí)現(xiàn)[J]. 計(jì)算機(jī)應(yīng)用,2019,39(S1):164-167.

    [20]趙春發(fā),翟婉明. 常導(dǎo)電磁懸浮動(dòng)態(tài)特性研究[J]. 西南交通大學(xué)學(xué)報(bào),2004,39(4):464-468.

    Influence of Levitation Control Parameters on Coupled Vibration Between Maglev Vehicles and Switch Girder

    YANG Zhinan1,F(xiàn)ENG Yang2,LIU Dongsheng2,ZHAO Chunfa2

    (1.China Railway First Survey and Design Institute Group Co., Ltd., Xi’an 710043, China;2.State Key Laboratory of Traction Power, Southwest Jiaotong University, Chengdu 610031, China)

    In this paper, a dynamic model of two-car maglev trainset is built considering the PID levitation control system, meanwhile, a detailed finite element model of the switch girder is also established. Then, they are used to calculate dynamic response of maglev train with different levitation control parameters whenit passes through the switch girder at a speed of 30km/h. The simulated results show that the levitation airgap fluctuation become larger when the airgap feedback coefficient is small, which is unfavorable for train ride safety. With the increase of airgap feedback coefficient, the characteristic frequency of electromagnetic levitation system gradually approaches the first-order vertical bending modal frequency of the switch girder, which leads to stronger coupled vibration. When the airgap velocity feedback coefficient is small, the damping of electromagnetic levitation system is small, so the vibration acceleration of carbody become larger. If the airgap velocity feedback coefficient is quite large, the coupled vibration is violent due totoo strong lead compensation effect. In order to alleviate the coupled vibration between maglev vehicles and the switch girder, it is suggested to select the airgap feedback coefficient between 7000 and 8000 and the airgap velocity feedback coefficient between 45 and 60.

    maglev train;switch;feedback control;coupled vibration;numerical simulation

    U237;U441+.3

    A

    10.3969/j.issn.1006-0316.2022.02.006

    1006-0316 (2022) 02-0038-09

    2021-07-09

    中鐵第一勘察設(shè)計(jì)院集團(tuán)有限公司課題:磁浮道岔優(yōu)化創(chuàng)新及選型研究(院科18-05);國(guó)家重點(diǎn)研發(fā)計(jì)劃課題子任務(wù):高速磁浮車輛動(dòng)力學(xué)性能匹配設(shè)計(jì)理論與技術(shù)(2016YFB1200602-15)

    楊志南(1991-),男,天津薊州人,工程師,主要從事鐵路站場(chǎng)設(shè)計(jì)工作,E-mail:1569156195@qq.com。*通訊作者:趙春發(fā)(1973-),男,湖北仙桃人,博士,研究員,主要研究方向?yàn)檐壍澜煌üこ虅?dòng)力學(xué),E-mail:cfzhao@swjtu.edu.cn。

    猜你喜歡
    控制參數(shù)電磁鐵磁懸浮
    自制簡(jiǎn)易電磁鐵
    奇妙的電磁鐵
    有趣的磁懸浮陀螺
    高超聲速飛行器滑??刂茀?shù)整定方法設(shè)計(jì)*
    JB/T 10159—2019《交流本整濕式閥用電磁鐵》解讀
    第二章 磁懸浮裝置
    一種電磁鐵吸取減速頂調(diào)整墊的工具
    Birkhoff系統(tǒng)穩(wěn)定性的動(dòng)力學(xué)控制1)
    神奇的磁懸浮
    神奇的磁懸浮
    高潮久久久久久久久久久不卡| 三级国产精品欧美在线观看 | 亚洲av第一区精品v没综合| 五月玫瑰六月丁香| 在线视频色国产色| 两个人的视频大全免费| 三级男女做爰猛烈吃奶摸视频| 国内精品一区二区在线观看| 最近最新中文字幕大全电影3| 国产欧美日韩一区二区精品| 中国美女看黄片| bbb黄色大片| 亚洲成a人片在线一区二区| 毛片女人毛片| 天天添夜夜摸| 少妇裸体淫交视频免费看高清 | 99在线视频只有这里精品首页| 禁无遮挡网站| 国产成人aa在线观看| 国产区一区二久久| 久久草成人影院| 成人三级做爰电影| 免费在线观看黄色视频的| www.www免费av| 大型av网站在线播放| 亚洲精品色激情综合| 欧美日韩亚洲综合一区二区三区_| 国产亚洲精品一区二区www| 久久久久久久久免费视频了| 91在线观看av| 欧美日韩精品网址| 亚洲av成人精品一区久久| 日韩欧美在线乱码| 欧美成人午夜精品| 国产精品影院久久| 在线免费观看的www视频| 天堂√8在线中文| 久久中文字幕一级| 久久精品成人免费网站| 制服诱惑二区| av福利片在线| 久9热在线精品视频| 欧美3d第一页| 国产激情久久老熟女| 亚洲第一电影网av| 免费在线观看日本一区| 9191精品国产免费久久| 亚洲一区二区三区不卡视频| 不卡av一区二区三区| 国产一级毛片七仙女欲春2| 一本大道久久a久久精品| 九九热线精品视视频播放| 亚洲一区二区三区不卡视频| 男女那种视频在线观看| 在线观看美女被高潮喷水网站 | 国产一区二区在线观看日韩 | 国产精品,欧美在线| 久久国产乱子伦精品免费另类| 麻豆久久精品国产亚洲av| 91字幕亚洲| 在线观看www视频免费| 国产亚洲精品久久久久久毛片| 不卡一级毛片| 欧美精品啪啪一区二区三区| 精品高清国产在线一区| 亚洲欧美精品综合一区二区三区| 一级毛片女人18水好多| 成人精品一区二区免费| 不卡av一区二区三区| 日本成人三级电影网站| 好男人电影高清在线观看| 欧美在线黄色| 88av欧美| 亚洲av日韩精品久久久久久密| 日日干狠狠操夜夜爽| 久久亚洲真实| 操出白浆在线播放| 午夜老司机福利片| 亚洲av电影不卡..在线观看| 看免费av毛片| 国产免费男女视频| 午夜视频精品福利| 日韩成人在线观看一区二区三区| 亚洲片人在线观看| 变态另类丝袜制服| 亚洲精品中文字幕一二三四区| 久久亚洲真实| 毛片女人毛片| 麻豆成人午夜福利视频| 久久精品国产99精品国产亚洲性色| 亚洲成人久久爱视频| 国产成人aa在线观看| 国产激情久久老熟女| 50天的宝宝边吃奶边哭怎么回事| 女人爽到高潮嗷嗷叫在线视频| 亚洲成人免费电影在线观看| 亚洲成人精品中文字幕电影| 两个人免费观看高清视频| av有码第一页| 日本a在线网址| 国产精华一区二区三区| 99久久精品国产亚洲精品| 精品久久久久久久久久久久久| 国产激情欧美一区二区| 曰老女人黄片| 国产一区二区三区视频了| 国产精品av视频在线免费观看| 黄色视频不卡| 91av网站免费观看| 国产亚洲av嫩草精品影院| 两人在一起打扑克的视频| 正在播放国产对白刺激| 女人爽到高潮嗷嗷叫在线视频| 啦啦啦韩国在线观看视频| 黑人欧美特级aaaaaa片| 精品久久久久久久人妻蜜臀av| 日日夜夜操网爽| 特级一级黄色大片| 悠悠久久av| 午夜免费观看网址| 日日干狠狠操夜夜爽| 一级作爱视频免费观看| 校园春色视频在线观看| 日日夜夜操网爽| 亚洲 欧美 日韩 在线 免费| 国产三级在线视频| 五月玫瑰六月丁香| 久久亚洲真实| 精品欧美国产一区二区三| 欧美3d第一页| 亚洲精品国产一区二区精华液| 欧美大码av| 国产激情偷乱视频一区二区| 少妇熟女aⅴ在线视频| 天天躁狠狠躁夜夜躁狠狠躁| 1024香蕉在线观看| 久久精品综合一区二区三区| 午夜亚洲福利在线播放| 日韩高清综合在线| 国产麻豆成人av免费视频| 级片在线观看| 国产伦一二天堂av在线观看| 手机成人av网站| 久久久国产欧美日韩av| 搡老岳熟女国产| 国产私拍福利视频在线观看| 免费高清视频大片| 欧美日韩亚洲综合一区二区三区_| 精品免费久久久久久久清纯| 中文字幕精品亚洲无线码一区| 一级毛片高清免费大全| 麻豆成人av在线观看| 国产精品久久久久久久电影 | 精品一区二区三区四区五区乱码| 久久久精品大字幕| 国产伦一二天堂av在线观看| 老熟妇仑乱视频hdxx| 国产成人系列免费观看| 五月玫瑰六月丁香| 国产av又大| 精品第一国产精品| 全区人妻精品视频| 精品乱码久久久久久99久播| 国产精品久久视频播放| 色av中文字幕| 岛国在线观看网站| www.自偷自拍.com| 色综合亚洲欧美另类图片| 亚洲七黄色美女视频| 久久精品国产99精品国产亚洲性色| 一本大道久久a久久精品| 日韩 欧美 亚洲 中文字幕| 久久久国产精品麻豆| 久久久国产成人免费| 91老司机精品| 欧美日韩黄片免| 免费高清视频大片| 国产成人aa在线观看| 久久精品国产综合久久久| 久久欧美精品欧美久久欧美| 亚洲精品国产精品久久久不卡| 听说在线观看完整版免费高清| а√天堂www在线а√下载| av超薄肉色丝袜交足视频| 午夜免费激情av| 日本撒尿小便嘘嘘汇集6| 观看免费一级毛片| 亚洲色图av天堂| 亚洲一区二区三区色噜噜| 亚洲av成人不卡在线观看播放网| 男女下面进入的视频免费午夜| 午夜福利免费观看在线| 俺也久久电影网| 中文在线观看免费www的网站 | 欧美黑人欧美精品刺激| 一边摸一边抽搐一进一小说| 欧美日韩亚洲国产一区二区在线观看| 国产亚洲精品综合一区在线观看 | 国产av又大| 国产精品一区二区免费欧美| 亚洲电影在线观看av| 国产精品,欧美在线| 男插女下体视频免费在线播放| 国内久久婷婷六月综合欲色啪| 老司机在亚洲福利影院| 黄色女人牲交| 1024视频免费在线观看| 亚洲第一电影网av| 久久香蕉精品热| 亚洲男人天堂网一区| a在线观看视频网站| 日韩欧美国产在线观看| 国产精品一区二区精品视频观看| 亚洲精品一区av在线观看| 91成年电影在线观看| bbb黄色大片| 窝窝影院91人妻| ponron亚洲| 美女黄网站色视频| 又大又爽又粗| 亚洲国产中文字幕在线视频| 精品不卡国产一区二区三区| 欧美黑人欧美精品刺激| 成年人黄色毛片网站| 国产精品永久免费网站| 小说图片视频综合网站| 久久久久久久久久黄片| 国产成人影院久久av| 夜夜躁狠狠躁天天躁| 悠悠久久av| 长腿黑丝高跟| 狂野欧美白嫩少妇大欣赏| www.精华液| 老司机在亚洲福利影院| 草草在线视频免费看| 精品久久久久久久毛片微露脸| 久99久视频精品免费| 国产亚洲精品综合一区在线观看 | 最近在线观看免费完整版| 日韩欧美国产在线观看| 国产亚洲精品久久久久久毛片| 国产精品一及| 久久精品91无色码中文字幕| 精品欧美一区二区三区在线| 黑人操中国人逼视频| av天堂在线播放| 国产av不卡久久| 久久精品成人免费网站| 精品久久久久久成人av| 日本成人三级电影网站| 欧美日韩中文字幕国产精品一区二区三区| 制服诱惑二区| netflix在线观看网站| 日本成人三级电影网站| 亚洲男人天堂网一区| 香蕉国产在线看| 亚洲国产中文字幕在线视频| 手机成人av网站| 美女午夜性视频免费| 亚洲精品国产精品久久久不卡| 又大又爽又粗| 天堂影院成人在线观看| 老司机在亚洲福利影院| 两个人免费观看高清视频| 国产一区二区在线av高清观看| 啪啪无遮挡十八禁网站| 国产av一区在线观看免费| 欧美日韩精品网址| 国产精品一及| 别揉我奶头~嗯~啊~动态视频| 欧美日韩亚洲国产一区二区在线观看| av免费在线观看网站| 亚洲全国av大片| 天天躁夜夜躁狠狠躁躁| www日本黄色视频网| 亚洲av熟女| 又爽又黄无遮挡网站| 国产精品98久久久久久宅男小说| 黄片小视频在线播放| 18禁黄网站禁片午夜丰满| 色综合婷婷激情| 丁香欧美五月| 1024视频免费在线观看| 亚洲精品av麻豆狂野| 亚洲自拍偷在线| 精品久久久久久久末码| 国产精品永久免费网站| 婷婷亚洲欧美| 久久婷婷人人爽人人干人人爱| 久久精品91无色码中文字幕| 男女那种视频在线观看| 国产爱豆传媒在线观看 | 一本精品99久久精品77| www.自偷自拍.com| 黄频高清免费视频| 久久香蕉激情| 99久久99久久久精品蜜桃| 美女免费视频网站| 免费看十八禁软件| 国内久久婷婷六月综合欲色啪| 黑人巨大精品欧美一区二区mp4| 国产成年人精品一区二区| 非洲黑人性xxxx精品又粗又长| 在线十欧美十亚洲十日本专区| 宅男免费午夜| 别揉我奶头~嗯~啊~动态视频| 韩国av一区二区三区四区| 成人国产一区最新在线观看| 亚洲av日韩精品久久久久久密| 无人区码免费观看不卡| 国产伦一二天堂av在线观看| 亚洲av成人不卡在线观看播放网| 亚洲aⅴ乱码一区二区在线播放 | 国产不卡一卡二| www.www免费av| 韩国av一区二区三区四区| 国产一区二区在线观看日韩 | 中文资源天堂在线| 精华霜和精华液先用哪个| 一级片免费观看大全| 国产激情欧美一区二区| 少妇粗大呻吟视频| 国产成人精品久久二区二区91| 亚洲欧洲精品一区二区精品久久久| 人妻夜夜爽99麻豆av| 19禁男女啪啪无遮挡网站| 亚洲欧美精品综合久久99| 啦啦啦韩国在线观看视频| 一区二区三区高清视频在线| а√天堂www在线а√下载| 特大巨黑吊av在线直播| 日本黄色视频三级网站网址| 天天添夜夜摸| 免费在线观看黄色视频的| 日日爽夜夜爽网站| 亚洲一区二区三区色噜噜| 18禁国产床啪视频网站| 18禁美女被吸乳视频| 欧美色视频一区免费| 日本一本二区三区精品| 成人av一区二区三区在线看| 国产精品日韩av在线免费观看| 欧美人与性动交α欧美精品济南到| 窝窝影院91人妻| 一进一出抽搐gif免费好疼| 美女高潮喷水抽搐中文字幕| 免费看十八禁软件| 国产精品久久电影中文字幕| 亚洲av五月六月丁香网| 草草在线视频免费看| 黄片大片在线免费观看| 中文字幕久久专区| 999精品在线视频| 黄色女人牲交| 两性夫妻黄色片| 真人做人爱边吃奶动态| 神马国产精品三级电影在线观看 | av视频在线观看入口| 国产aⅴ精品一区二区三区波| 久久久久久九九精品二区国产 | 日韩大尺度精品在线看网址| 国产一区二区三区视频了| 在线观看美女被高潮喷水网站 | 最新在线观看一区二区三区| 成人高潮视频无遮挡免费网站| 欧美成人午夜精品| 香蕉丝袜av| 老司机午夜福利在线观看视频| 国产99白浆流出| 狠狠狠狠99中文字幕| 亚洲在线自拍视频| 成人三级黄色视频| 国产亚洲av嫩草精品影院| 亚洲国产日韩欧美精品在线观看 | 亚洲成人中文字幕在线播放| 久久九九热精品免费| 欧美日韩国产亚洲二区| 免费电影在线观看免费观看| 欧美丝袜亚洲另类 | 久久久久久久久免费视频了| 男女午夜视频在线观看| 国产精品久久视频播放| 国内久久婷婷六月综合欲色啪| 久久亚洲精品不卡| 国产高清videossex| 国产私拍福利视频在线观看| 精品第一国产精品| 国产黄a三级三级三级人| 一级毛片女人18水好多| 九九热线精品视视频播放| 久久久久久久久中文| 亚洲av电影在线进入| 亚洲国产日韩欧美精品在线观看 | 哪里可以看免费的av片| 一级毛片高清免费大全| 香蕉丝袜av| 99国产精品一区二区三区| 岛国在线免费视频观看| 麻豆成人av在线观看| 国产成人精品久久二区二区91| 欧美高清成人免费视频www| 岛国视频午夜一区免费看| 桃色一区二区三区在线观看| 在线a可以看的网站| 国产三级黄色录像| 精品久久久久久久久久免费视频| 中文字幕人成人乱码亚洲影| 高清毛片免费观看视频网站| 9191精品国产免费久久| 日韩精品免费视频一区二区三区| 99久久无色码亚洲精品果冻| 国产av一区在线观看免费| 不卡一级毛片| 国产一区二区三区在线臀色熟女| 成年女人毛片免费观看观看9| 国产精品1区2区在线观看.| 亚洲男人天堂网一区| 三级毛片av免费| 岛国在线观看网站| 1024香蕉在线观看| 精品熟女少妇八av免费久了| 国产成+人综合+亚洲专区| 国产精品久久久人人做人人爽| 国产真实乱freesex| 欧美性长视频在线观看| 久久久精品欧美日韩精品| 日韩大尺度精品在线看网址| 欧美黑人精品巨大| 亚洲一区高清亚洲精品| 日韩欧美 国产精品| 两个人视频免费观看高清| 欧美高清成人免费视频www| 男女床上黄色一级片免费看| 操出白浆在线播放| 国产精品国产高清国产av| 国产精品爽爽va在线观看网站| 精品第一国产精品| 国产主播在线观看一区二区| 国产区一区二久久| 变态另类成人亚洲欧美熟女| 欧美成人性av电影在线观看| 欧美高清成人免费视频www| 色综合亚洲欧美另类图片| 又爽又黄无遮挡网站| 男人舔奶头视频| 毛片女人毛片| 老熟妇乱子伦视频在线观看| 叶爱在线成人免费视频播放| 50天的宝宝边吃奶边哭怎么回事| 国产单亲对白刺激| 日韩欧美国产一区二区入口| 麻豆国产av国片精品| 黑人巨大精品欧美一区二区mp4| 免费搜索国产男女视频| 特级一级黄色大片| 中文字幕熟女人妻在线| 精品久久久久久久人妻蜜臀av| 99在线人妻在线中文字幕| 亚洲国产欧美网| 两个人视频免费观看高清| 小说图片视频综合网站| 女警被强在线播放| 中文字幕人成人乱码亚洲影| 亚洲av成人精品一区久久| 国产高清视频在线观看网站| 9191精品国产免费久久| 欧美黑人巨大hd| 国产一区二区三区视频了| 亚洲欧美激情综合另类| 一区福利在线观看| 日韩大尺度精品在线看网址| 在线观看美女被高潮喷水网站 | 欧美性猛交╳xxx乱大交人| 757午夜福利合集在线观看| 十八禁人妻一区二区| 999久久久精品免费观看国产| av福利片在线| 五月伊人婷婷丁香| 国产又黄又爽又无遮挡在线| 久久久久久国产a免费观看| 女人爽到高潮嗷嗷叫在线视频| 夜夜爽天天搞| 人妻夜夜爽99麻豆av| 成年免费大片在线观看| 欧美午夜高清在线| 国产成人精品久久二区二区免费| 亚洲黑人精品在线| 一区二区三区高清视频在线| 精品欧美一区二区三区在线| 精品乱码久久久久久99久播| 欧美日本视频| 老鸭窝网址在线观看| 亚洲精品中文字幕在线视频| 女人被狂操c到高潮| 他把我摸到了高潮在线观看| 看黄色毛片网站| 免费人成视频x8x8入口观看| 999久久久精品免费观看国产| 99久久精品国产亚洲精品| 婷婷六月久久综合丁香| 国产99久久九九免费精品| 午夜亚洲福利在线播放| 三级男女做爰猛烈吃奶摸视频| 一级作爱视频免费观看| 啦啦啦免费观看视频1| 91大片在线观看| 国产一区二区三区在线臀色熟女| 黄色成人免费大全| 国产成人一区二区三区免费视频网站| 午夜福利欧美成人| 亚洲 欧美一区二区三区| 久久中文字幕一级| 日韩中文字幕欧美一区二区| 很黄的视频免费| 亚洲精品久久成人aⅴ小说| 99精品在免费线老司机午夜| 亚洲av熟女| 国产一区二区三区视频了| 性欧美人与动物交配| 丰满的人妻完整版| 国产亚洲精品av在线| 久久久久国内视频| 大型黄色视频在线免费观看| 亚洲av电影不卡..在线观看| 亚洲av成人av| 757午夜福利合集在线观看| 精品欧美国产一区二区三| 国产成人系列免费观看| 怎么达到女性高潮| 性色av乱码一区二区三区2| 久久久久久免费高清国产稀缺| 亚洲一区高清亚洲精品| 日本精品一区二区三区蜜桃| 亚洲国产日韩欧美精品在线观看 | 国产高清有码在线观看视频 | 99riav亚洲国产免费| 亚洲全国av大片| 国产久久久一区二区三区| АⅤ资源中文在线天堂| 曰老女人黄片| 无人区码免费观看不卡| 免费av毛片视频| 日韩欧美在线乱码| 久久精品成人免费网站| 2021天堂中文幕一二区在线观| 亚洲天堂国产精品一区在线| 两性午夜刺激爽爽歪歪视频在线观看 | www.熟女人妻精品国产| 国产黄a三级三级三级人| 婷婷亚洲欧美| 无人区码免费观看不卡| 久久这里只有精品中国| 国产1区2区3区精品| 老鸭窝网址在线观看| 国产三级中文精品| 亚洲国产精品久久男人天堂| 色综合站精品国产| av福利片在线观看| 国产麻豆成人av免费视频| 国产免费av片在线观看野外av| 99国产精品99久久久久| 亚洲五月天丁香| 伊人久久大香线蕉亚洲五| 国产片内射在线| 国产精品爽爽va在线观看网站| 黑人巨大精品欧美一区二区mp4| 两人在一起打扑克的视频| 精品国产美女av久久久久小说| 好看av亚洲va欧美ⅴa在| 午夜精品在线福利| 可以在线观看毛片的网站| 久久人妻福利社区极品人妻图片| 亚洲一区二区三区色噜噜| 欧美日韩亚洲国产一区二区在线观看| 亚洲精品国产精品久久久不卡| www日本在线高清视频| 亚洲成av人片免费观看| 一a级毛片在线观看| 成人欧美大片| 久久精品成人免费网站| 最好的美女福利视频网| xxx96com| 狠狠狠狠99中文字幕| 男人舔女人的私密视频| 午夜久久久久精精品| 一个人免费在线观看电影 | 欧美黄色片欧美黄色片| 在线十欧美十亚洲十日本专区| 女警被强在线播放| 他把我摸到了高潮在线观看| 亚洲自偷自拍图片 自拍| 欧美日韩福利视频一区二区| 国产欧美日韩一区二区三| 亚洲成av人片在线播放无| 欧美成狂野欧美在线观看| www日本黄色视频网| 中文字幕高清在线视频| 这个男人来自地球电影免费观看| 国产成人系列免费观看| 99国产精品99久久久久| 午夜福利视频1000在线观看| 99re在线观看精品视频| 欧美高清成人免费视频www| 香蕉国产在线看| 亚洲美女视频黄频| 一本综合久久免费| 成人av在线播放网站| 欧美午夜高清在线| 亚洲最大成人中文| 非洲黑人性xxxx精品又粗又长| 亚洲精品色激情综合| 午夜激情av网站| 欧美3d第一页| av有码第一页| 欧美成人午夜精品|