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

    重載鐵路軌道剛度對貨物列車抗脫軌安全度的影響

    2017-02-09 06:06:02余翠英
    關(guān)鍵詞:道床軌枕扣件

    龔 凱 向 俊 余翠英,3

    (1中南大學(xué)土木工程學(xué)院, 長沙 410075)(2華東交通大學(xué)土木建筑學(xué)院, 南昌 330013)(3南昌理工學(xué)院建筑工程學(xué)院, 南昌 330044)

    重載鐵路軌道剛度對貨物列車抗脫軌安全度的影響

    龔 凱1,2向 俊1余翠英1,3

    (1中南大學(xué)土木工程學(xué)院, 長沙 410075)(2華東交通大學(xué)土木建筑學(xué)院, 南昌 330013)(3南昌理工學(xué)院建筑工程學(xué)院, 南昌 330044)

    為確定具有預(yù)防貨物列車脫軌功能的軌道剛度合理值,基于列車軌道系統(tǒng)空間振動計(jì)算模型及列車脫軌能量隨機(jī)分析方法,提出了貨物列車-軌道(FTT)系統(tǒng)橫向振動穩(wěn)定性分析方法,分析多組扣件及道床橫向剛度組合下FTT系統(tǒng)抗脫軌能力、FTT系統(tǒng)橫向振動穩(wěn)定性及其振動響應(yīng)的影響.結(jié)果表明:FTT系統(tǒng)抗脫軌能力、臨界車速及容許極限車速隨著扣件及道床橫向剛度的增大均有大幅度提高,但當(dāng)扣件和道床橫向剛度分別增大至90和10 MN/m時(shí),其提高幅度逐漸減小,且當(dāng)扣件和道床橫向剛度分別由120 MN/m增至150 MN/m、15 MN/m增至20 MN/m時(shí)FTT系統(tǒng)抗脫軌能力、臨界車速及容許極限車速僅提高了3.9%,1.8%和1.8%;另外,增大扣件和道床橫向剛度有助于減小軌道橫向位移.考慮日趨緊張的重載鐵路市場競爭,建議扣件橫向剛度取90~120 MN/m, 道床橫向剛度取10~15 MN/m.

    鐵道工程;重載鐵路;軌道剛度;FTT系統(tǒng)橫向振動穩(wěn)定性;貨物列車

    隨著重載鐵路行車密度、軸重及車速的增大,軌道結(jié)構(gòu)安全儲備日益下降,甚至危及行車安全.文獻(xiàn)[1]研究表明,橋梁橫向剛度不足易引起橋上列車脫軌.同理,軌道橫向剛度不夠亦會引起列車脫軌.為保證貨物列車安全、平穩(wěn)運(yùn)行,增大軌道橫向剛度十分必要.實(shí)踐表明,較大的軌道橫向剛度可控制列車脫軌,如高速鐵路無砟軌道的剛度非常大,直接保證了高速列車的運(yùn)行安全性.剛度過低會嚴(yán)重影響貨物列車的安全性及平穩(wěn)性,而剛度過高,雖然抗脫軌安全度很高,但會增加重載鐵路的造價(jià),而在日趨緊張的國際重載鐵路市場競爭中,價(jià)格優(yōu)勢至關(guān)重要.因此,亟需對具有預(yù)防貨物列車脫軌功能的重載鐵路軌道剛度合理值進(jìn)行研究.為此,本文通過探討現(xiàn)有軌道剛度合理值的評價(jià)方法及其不足,基于列車-軌道系統(tǒng)空間振動計(jì)算模型及列車脫軌能量隨機(jī)分析方法,提出貨物列車-軌道(freight train-track system, FTT)系統(tǒng)橫向振動穩(wěn)定性分析方法.采用此法分析了扣件及道床橫向剛度對FTT系統(tǒng)抗脫軌能力、臨界車速、容許極限車速以及FTT系統(tǒng)空間振動特性的影響,為具有預(yù)防貨物列車脫軌功能的重載鐵路軌道剛度分析方法及其合理值選取提供參考.

    1 FTT系統(tǒng)空間振動計(jì)算模型

    根據(jù)有砟軌道結(jié)構(gòu)特點(diǎn),將計(jì)算長度范圍內(nèi)軌道離散為n個(gè)等長的軌段單元,每個(gè)軌段單元自上而下均分為2層,一層為扣件層;另一層為軌枕與道床之間的連接層,該層假定為道床頂面的彈性層.采用線性彈簧及黏滯阻尼器模擬,相應(yīng)的豎向彈性系數(shù)和阻尼系數(shù)分別為K1,C1和K5,C5,橫向分別為K2,C2和K4,C4;K3,C3分別為軌枕與道床之間的縱向彈性系數(shù)和阻尼系數(shù).鋼軌視為彈性點(diǎn)支承Euler梁,軌枕視為彈性變形體,但不計(jì)軸向變形及扭轉(zhuǎn)變形,這樣可將每個(gè)軌段單元模擬為具有34個(gè)自由度的有限元模型.任一軌段單元的節(jié)點(diǎn)位移為

    (1)

    (2)

    (3)

    式中,下標(biāo)1, 2分別表示軌段單元的左端和右端節(jié)點(diǎn);上標(biāo)T,S分別表示鋼軌和軌枕;下標(biāo)R,L表示軌段單元的右側(cè)和左側(cè);U,V,W分別為沿X,Y,Z軸方向線位移;θ為轉(zhuǎn)角位移;下標(biāo)X,Y,Z分別為轉(zhuǎn)角位移所繞的坐標(biāo)軸;γ為鋼軌X軸方向的扭轉(zhuǎn)角變化率.軌道空間振動計(jì)算模型見圖1.

    (a) 三視圖

    (b) 端視圖

    根據(jù)模型(1)~(3)可導(dǎo)出每個(gè)軌段單元的空間振動勢能ΠTj及軌道空間振動總勢能ΠT[2].設(shè)t時(shí)刻,一列具有m節(jié)車廂的貨物列車運(yùn)行在長度為L的軌道上,此時(shí),軌道上的貨物列車空間振動勢能為ΠV[2].考慮輪軌相對位移銜接條件及輪軌“游間”影響,可得貨物列車-軌道系統(tǒng)在t時(shí)刻的空間振動總勢能為

    Π=ΠV+ΠT

    (4)

    根據(jù)式(4),按照彈性系統(tǒng)動力學(xué)總勢能不變值原理[3]及形成系統(tǒng)矩陣的“對號入座”法則[1],得出FTT系統(tǒng)在t時(shí)刻的矩陣方程為

    (5)

    2 現(xiàn)有軌道剛度合理值評價(jià)方法及存在問題分析

    2.1 現(xiàn)有軌道剛度合理值評價(jià)方法

    為確定軌道剛度合理值,國內(nèi)外學(xué)者針對軌道剛度問題進(jìn)行了許多研究,如Lopez-Pita等[4]研究了高速鐵路線路惡化后軌道豎向剛度的變化.Real等[5]、Dahlberg[6]及Giannakos[7]探討了軌道剛度對軌道豎向位移、各部件應(yīng)力的影響及消除軌道剛度突變的措施.Frohling等[8]研究了軌道剛度變化對車軌豎向動力響應(yīng)的影響.Wu等[9]分析了輪載重復(fù)加載下輪載附近軌道基礎(chǔ)剛度的變化規(guī)律.翟婉明等[10-11]探討了軌道豎向剛度對輪軌力、車體、構(gòu)架加速度等指標(biāo)的影響;同時(shí),基于極限環(huán)理論的臨界車速,研究了軌道結(jié)構(gòu)扣件橫向剛度和阻尼對高速客車臨界速度的影響及其合理取值.劉學(xué)毅[12]探討了軌道剛度對輪軌力、鋼軌動位移及加速度等指標(biāo)的影響.練松良[13]提出了軌下墊層、扣件和道床剛度的計(jì)算方法.向俊等[14-15]分別研究了有砟軌道、板式軌道等軌道剛度對車軌系統(tǒng)空間振動響應(yīng)的影響.張斌等[16]研究了CRTSⅡ型無砟軌道軌下墊板、CA砂漿層等結(jié)構(gòu)剛度對軌道結(jié)構(gòu)振動響應(yīng)的影響.王平等[17]研究了既有提速線軌道剛度對車輛軌道耦合系統(tǒng)低、中、高頻振動響應(yīng)的影響.張格明[18]和趙國堂[19]研究了軌道整體剛度對枕上壓力、輪軌力等動力指標(biāo)的影響.

    2.2 存在的主要問題

    軌道剛度合理值評價(jià)方法可歸納為如下3類:

    1) 動力學(xué)指標(biāo)評價(jià)法.通過分析軌道剛度對輪軌垂向力、輪軸橫向力、鋼軌位移等指標(biāo)的影響,進(jìn)而評價(jià)軌道剛度合理值.這些指標(biāo)能夠反映列車正常運(yùn)行時(shí)的通常值,但難以反映列車脫軌瞬間的極限值,并且未反映列車的脫軌信息.

    2) 列車安全性指標(biāo)評價(jià)法.采用脫軌系數(shù)、輪重減載率等安全性指標(biāo)對軌道剛度進(jìn)行評價(jià).但脫軌系數(shù)僅表達(dá)了輪軌一點(diǎn)接觸和開始懸浮的狀態(tài),其限值是基于單輪對假定荷載下的分析結(jié)果,與實(shí)際列車運(yùn)動趨勢存在差異,難以反映實(shí)際的脫軌信息;同時(shí),因脫軌試驗(yàn)難以實(shí)現(xiàn),脫軌瞬間的脫軌系數(shù)和輪重減載率未能夠測得.

    3) 基于極限環(huán)理論的臨界車速評價(jià)法.通過在軌道上施加軌道不平順來激發(fā)整個(gè)系統(tǒng)的振動, 然后使車輛在平直的、沒有軌道不平順的軌道上運(yùn)行,當(dāng)系統(tǒng)的振動響應(yīng)不再衰減到平衡位置而是趨于穩(wěn)定的極限環(huán)時(shí), 此時(shí)的速度即為車輛臨界速度.但系統(tǒng)振動響應(yīng)沒有衰減到平衡位置并不能反映列車脫軌,因?yàn)榱熊囓壍老到y(tǒng)是運(yùn)動系統(tǒng),運(yùn)動系統(tǒng)存在持續(xù)振動,運(yùn)動系統(tǒng)振動處于臨界失穩(wěn)狀態(tài)時(shí)所對應(yīng)的車速才是列車的臨界車速.

    3 重載鐵路FTT系統(tǒng)橫向振動穩(wěn)定性分析方法

    鑒于上述分析,提出重載鐵路FTT系統(tǒng)橫向振動穩(wěn)定性分析方法.文獻(xiàn)[20]針對高速列車運(yùn)行安全性問題提出了列車軌道系統(tǒng)橫向振動穩(wěn)定性分析方法,但重載鐵路有砟軌道不同于高速鐵路無砟軌道,并且車軌系統(tǒng)也不盡相同.

    3.1 FTT系統(tǒng)橫向振動穩(wěn)定性判別準(zhǔn)則

    文獻(xiàn)[1]表明,列車脫軌機(jī)理是FTT系統(tǒng)橫向振動狀態(tài)喪失穩(wěn)定.要判斷列車是否脫軌,必須判別FTT系統(tǒng)橫向振動狀態(tài)是否失穩(wěn).由物理概念及理論分析表明,當(dāng)干擾引起的系統(tǒng)抗力增量大于荷載增量時(shí),由干擾引起的位移及振動逐漸消失,系統(tǒng)恢復(fù)穩(wěn)定狀態(tài);而當(dāng)干擾引起的系統(tǒng)抗力增量小于荷載增量時(shí),由干擾引起的位移及振動不斷增長,系統(tǒng)穩(wěn)定狀態(tài)消亡.同理,因FTT系統(tǒng)具有自激性[1],當(dāng)列車在軌道上運(yùn)行時(shí)不斷有能量(即干擾)輸入FTT系統(tǒng),而此干擾使得FTT系統(tǒng)振動不斷增大,當(dāng)FTT系統(tǒng)橫向振動抗力做功等于輸入能量時(shí),FTT系統(tǒng)處于失穩(wěn)臨界狀態(tài).由此,FTT系統(tǒng)橫向振動狀態(tài)是否穩(wěn)定可通過FTT系統(tǒng)輸入能量及其抗力做功來判別.然而,FTT系統(tǒng)具有時(shí)變性且橫向振動方程數(shù)很多,難以建立輸入能量及抗力做功間的平衡準(zhǔn)則.設(shè)Δσpr為FTT系統(tǒng)橫向振動輸入能量增量,Δσcr為FTT系統(tǒng)極限抗力做功增量,則根據(jù)評判系統(tǒng)平衡狀態(tài)是否穩(wěn)定的能量增量思想[1],可建立FTT系統(tǒng)橫向振動穩(wěn)定性判別準(zhǔn)則:

    1) Δσcr>Δσpr,FTT系統(tǒng)橫向振動穩(wěn)定;

    2) Δσcr<Δσpr,FTT系統(tǒng)橫向振動不穩(wěn)定;

    3) Δσcr=Δσpr,FTT系統(tǒng)橫向振動處于失穩(wěn)臨界狀態(tài).

    3.2 FTT系統(tǒng)橫向振動輸入能量增量及抗力做功增量

    文獻(xiàn)[21]提出采用構(gòu)架蛇形波作為FTT系統(tǒng)橫向振動激振源,并以構(gòu)架蛇形波標(biāo)準(zhǔn)差σp作為FTT系統(tǒng)橫向振動的輸入能量,并繪制了構(gòu)架蛇形波標(biāo)準(zhǔn)差σp與車速V的關(guān)系曲線(σp-V曲線).隨著FTT系統(tǒng)輸入能量的積聚,FTT系統(tǒng)橫向振動逐漸增大,直至其橫向振動狀態(tài)失穩(wěn),列車脫軌,此時(shí),構(gòu)架蛇形波標(biāo)準(zhǔn)差為σpmax, 它是輸入FTT系統(tǒng)的最大輸入能量.然而,實(shí)踐表明有限次的試驗(yàn)難以實(shí)現(xiàn)列車脫軌,因而難以獲得相應(yīng)的σpmax.但由功能轉(zhuǎn)換原理可知,構(gòu)架蛇形波標(biāo)準(zhǔn)差σpmax必然大于σp,因而σpmax-V曲線必然位于σp-V曲線的上方,且2條曲線平行[1].因此,當(dāng)車速V提高至車速Vr時(shí),列車正常運(yùn)行時(shí)FTT系統(tǒng)橫向振動輸入能量增量Δσpr與列車脫軌時(shí)FTT系統(tǒng)橫向振動輸入能量增量Δσpr,max相等.

    另外,由靜力壓桿穩(wěn)定性分析中臨界荷載Pcr的確定過程可知, FTT系統(tǒng)抵抗輸入能量的極限能力(即FTT系統(tǒng)橫向振動極限抗力做功σc)也可通過不斷加載的方式確定.這樣,根據(jù)車輪脫軌幾何準(zhǔn)則[1],FTT系統(tǒng)橫向振動極限抗力做功σc可通過試算得出,試算過程見文獻(xiàn)[1].通過建立σc-V曲線,確定車速V~Vr時(shí)FTT系統(tǒng)橫向振動極限抗力做功σc0與σcr之差,即為FTT系統(tǒng)橫向振動極限抗力做功增量Δσcr.實(shí)際上,這里的σc反映了FTT系統(tǒng)抗脫軌的能力.

    3.3 FTT系統(tǒng)橫向振動失穩(wěn)臨界車速及容許極限車速

    按準(zhǔn)則1)~3),根據(jù)Δσpr及Δσcr,以不脫軌車速V為起點(diǎn),依次判別Vr車速下FTT系統(tǒng)橫向振動穩(wěn)定性.將Δσcr=Δσpr時(shí)對應(yīng)的車速稱為FTT系統(tǒng)橫向振動失穩(wěn)臨界車速Vcr,同時(shí),將Vcr除以安全系數(shù)K(本文采用脫軌系數(shù)進(jìn)行列車安全性分析時(shí)的系數(shù)為1.25)可得控制貨物列車安全運(yùn)行的容許極限車速VL.由此,建立貨物列車安全運(yùn)行控制指標(biāo)為臨界車速Vcr和容許極限車速VL.

    在運(yùn)用上述分析方法評價(jià)軌道剛度對貨物列車抗脫軌安全度的影響中,需分別計(jì)算軌道各部件剛度組合下的σc-V曲線.分析軌道剛度對FTT系統(tǒng)抗脫軌能力的影響,并計(jì)算相應(yīng)的Δσc及Δσp.通過判斷FTT系統(tǒng)橫向振動穩(wěn)定性,得出Δσc=Δσp時(shí)的Vcr和VL.同時(shí),以VL為基礎(chǔ),計(jì)算并分析軌道各部件剛度組合下貨物列車以小于VL的車速運(yùn)行時(shí)FTT系統(tǒng)空間振動響應(yīng),采用國際通用的Sperling平穩(wěn)性指標(biāo)評價(jià)軌道剛度對貨物列車運(yùn)行平穩(wěn)性的影響.將同時(shí)滿足貨物列車運(yùn)行安全性和平穩(wěn)性要求的軌道剛度作為能夠預(yù)防貨物列車脫軌的軌道剛度合理值.

    4 算例分析

    鑒于全列空載貨車發(fā)生脫軌概率最大[1],本文將列車計(jì)算編組考慮為1輛機(jī)車牽引1輛空載敞車.同時(shí),因貨物列車脫軌車速常大于50 km/h,故本文以50 km/h作為計(jì)算起點(diǎn)[1],以車速每增加10 km/h計(jì)算一次FTT系統(tǒng)橫向振動穩(wěn)定性.軌道計(jì)算長度為200 m,直線軌道,鋼軌單位質(zhì)量60 kg/m,Ⅱ型混凝土軌枕,軌枕間距為0.543 5 m,碎石道砟.K2分別取為30,60,90,120和150 MN/m,K4分別取為3.808,5,10,15和20 MN/m.FTT系統(tǒng)橫向振動輸入能量σp,當(dāng)車速小于等于90 km/h時(shí)σp參考文獻(xiàn)[1], 當(dāng)車速大于90 km/h時(shí),σp參考文獻(xiàn)[21],由此得到σp-V曲線如圖2所示.

    圖2 σp-V曲線

    因目前尚缺少軌道橫向剛度增大后FTT系統(tǒng)橫向振動輸入能量等相關(guān)資料,故在判別FTT系統(tǒng)橫向振動穩(wěn)定性時(shí), FTT系統(tǒng)橫向振動輸入能量σp按照圖2中的σp-V曲線取值.實(shí)際上,增大剛度會使輸入能量減小,若增大剛度后的FTT系統(tǒng)仍能承受剛度增大前的輸入能量,則實(shí)際軌道橫向剛度增大后形成的FTT系統(tǒng)將具備更大的安全裕量.

    4.1 扣件及道床橫向剛度對貨物列車運(yùn)行安全性的影響分析

    1)K2=30 MN/m.由圖3表明,當(dāng)K2=30 MN/m時(shí),隨K4及V的增大,σc隨之增大.當(dāng)K4為3.808,5,10,15和20 MN/m時(shí),σc最大分別為317,420,479,498和506 cm/s2.當(dāng)K4為5,10,15和20 MN/m時(shí),FTT系統(tǒng)抗脫軌能力分別較K4為3.808MN/m時(shí)提高了32.5%,51.1%,57.1%和59.6%.

    圖3 K2=30 MN/m時(shí)不同道床橫向剛度條件下σc隨V的變化

    由表1可知,當(dāng)K2為30 MN/m、K4為3.808 MN/m時(shí),在車速120~130 km/h范圍內(nèi)必然存在臨界車速使得Δσcr=Δσpr.按照內(nèi)插法,可得Vcr=122.72 km/h,考慮安全系數(shù)1.25,相應(yīng)的VL=98.18 km/h;同理,可得K2為30 MN/m、K4分別為5,10,15和20 MN/m時(shí),Δσcr=Δσpr對應(yīng)的Vcr分別為146.67,153.9,155.83及158.18 km/h,VL分別為117.34,123.1,124.67及126.56 km/h.

    表1 K2=30 MN/m時(shí)各道床橫向剛度條件下FTT系統(tǒng)橫向振動穩(wěn)定性計(jì)算結(jié)果 cm/s2

    注:Δσ=Δσc-Δσp.

    2)K2=60 MN/m.圖4表明,K2為60 MN/m、K4分別為3.808,5,10,15和20 MN/m時(shí),σc最大分別為425,467,505,515和520 cm/s2.相比K2,K4增大前,K2為60 MN/m、K4分別為3.808,5,10,15和20 MN/m時(shí)FTT系統(tǒng)抗脫軌能力分別提高了34.1%,47.3%,59.3%,62.5%及64.0%.

    由表2可知,當(dāng)K2為60 MN/m、K4為3.808 MN/m時(shí),在車速140~150 km/h范圍內(nèi)必然存在臨界車速使得Δσcr=Δσpr.通過計(jì)算可得Vcr為145.72 km/h,VL為116.58 km/h.當(dāng)K2為60 MN/m、K4分別為5,10,15及20 MN/m時(shí),Δσcr=Δσpr對應(yīng)的Vcr分別為153.08,155.29,157.5及159.17 km/h,VL分別為122.46,124.2,126及127.33 km/h.

    圖4 K2=60 MN/m時(shí)不同道床橫向剛度條件下σc隨V的變化

    3)K2=90 MN/m.5表明,K2為90 MN/m、K4分別為3.808,5,10,15及20 MN/m時(shí),FTT系統(tǒng)橫向振動極限抗力做功最大分別469,506,549,559和568 cm/s2.相比K2,K4增大前,K2為60 MN/m、K4分別為3.808,5,10,15及20 MN/m時(shí),FTT系統(tǒng)抗脫軌能力分別提高了47.9%,59.6%,73.1%,76.3%和79.2%.

    由表3可知,當(dāng)K2為90 MN/m、K4為3.808 MN/m時(shí),在車速150~160 km/h范圍內(nèi)必然存在臨界車速使得Δσcr=Δσpr.由計(jì)算可得Vcr為152.5 km/h、VL為122.0 km/h.當(dāng)K2為90 MN/m、K4分別為5,10,15及20 MN/m時(shí),Δσcr=Δσpr對應(yīng)的Vcr分別為161.82,163,164.45及166.0 km/h,VL分別為129.45,130.41,131.56及132.8 km/h.

    圖5 K2=90 MN/m時(shí)不同道床橫向剛度條件下σc隨V的變化

    表2 K2=60 MN/m時(shí)各道床橫向剛度條件下FTT系統(tǒng)橫向振動穩(wěn)定性計(jì)算結(jié)果 cm/s2

    表3 K2=90 MN/m時(shí)各道床橫向剛度條件下FTT系統(tǒng)橫向振動穩(wěn)定性計(jì)算結(jié)果 cm/s2

    4)K2=120 MN/m.由圖6表明,K2為120 MN/m時(shí)各道床橫向剛度對應(yīng)的σc最大分別480,521,554,567和581 cm/s2.相比K2,K4增大前,K2為120 MN/m、K4為3.808,5,10,15及20 MN/m條件下FTT系統(tǒng)抗脫軌能力分別提高了51.4%,64.4%,74.8%,78.9%和83.3%.

    圖6 K2=120 MN/m時(shí)不同道床橫向剛度條件下σc隨V的變化

    由表4可知,當(dāng)K2=120 MN/m、K4=3.808 MN/m、V為150~160 km/h時(shí),必然存在臨界車速使得Δσcr=Δσpr,此時(shí)Vcr為155.38 km/h、VL為124.31 km/h.同理,當(dāng)K2=120 MN/m、K4分別為5,10,15及20 MN時(shí),Δσcr=Δσpr對應(yīng)的Vcr分別為162.5,163.85,165.0及167.14 km/h,VL分別為130.0,131.08,132.0及133.71 km/h.

    5)K2=150 MN/m.由圖7表明,當(dāng)K2為150 MN/m、K4分別為3.808,5,10,15及20 MN/m時(shí),σc最大分別493,536,570,577和589 cm/s2.相比軌道結(jié)構(gòu)強(qiáng)化前,K2為150 MN/m時(shí)各道床橫向剛度對應(yīng)的FTT系統(tǒng)抗脫軌能力分別提高了55.5%,69.1%,79.8%,82.0%和85.7%.

    由表5可知,當(dāng)K2為150 MN/m、K4為3.808 MN/m時(shí),在V為150~160 km/h范圍內(nèi)必然存在臨界車速使得Δσcr=Δσpr,此時(shí)Vcr為157.14 km/h、VL為125.71 km/h.同理,當(dāng)K2為150 MN/m、K4分別為5,10,15及20 MN/m時(shí),Δσcr=Δσpr對應(yīng)的Vcr分別為163.63,164.29,166.25及168.0 km/h,VL分別130.91,131.43,133.0及134.4 km/h.

    表4 K2=120 MN/m時(shí)各道床橫向剛度條件下FTT系統(tǒng)橫向振動穩(wěn)定性計(jì)算結(jié)果 cm/s2

    表5 K2=150 MN/m時(shí)各道床橫向剛度條件下FTT系統(tǒng)橫向振動穩(wěn)定性計(jì)算結(jié)果 cm/s2

    綜上所述,F(xiàn)TT系統(tǒng)抗脫軌能力、臨界車速及容許極限車速隨著扣件及道床橫向剛度的增大均有大幅度提高,但當(dāng)扣件及道床橫向剛度增大至90和10 MN/m時(shí),其提高幅度逐漸減小,且當(dāng)扣件及道床橫向剛度分別由120 MN/m增至150 MN/m、15 MN/m增至20 MN/m時(shí)FTT系統(tǒng)抗脫軌能力、臨界車速及容許極限車速僅提高了3.9%,1.8%和1.8%.

    4.2 扣件及道床橫向剛度對FTT系統(tǒng)振動響應(yīng)的影響

    按照第3節(jié)的分析方法,計(jì)算V=80 km/h時(shí)多組K2,K4條件下貨物列車不脫軌(即V

    (a) 車體橫向Sperling平穩(wěn)性指標(biāo)

    (b) 鋼軌橫向位移

    (c) 軌枕橫向位移

    由圖8(a)可見,車體橫向Sperling指標(biāo)隨K2的增大而明顯減小,而隨K4的增大變化較小;當(dāng)K2為90~120 MN/m、K4為10~15 MN/m時(shí),車體橫向Sperling平穩(wěn)性指標(biāo)變化平緩.

    由圖8(b)可見,當(dāng)K4為10~20 MN/m時(shí),隨著K2的增大,鋼軌、軌枕橫向位移變化很小.K4增大到一定時(shí),K2增大對減小鋼軌、軌枕橫向位移的貢獻(xiàn)較小.同時(shí),當(dāng)K2為90~150 MN/m時(shí),K4>10 MN/m后,鋼軌、軌枕橫向位移減小趨勢平緩.

    5 結(jié)論

    1) 基于列車-軌道系統(tǒng)空間振動模型及列車脫軌能量隨機(jī)分析方法,提出了貨物列車-軌道FTT系統(tǒng)橫向振動穩(wěn)定性分析方法.

    2) 隨扣件及道床橫向剛度的增大,FTT系統(tǒng)抗脫軌能力、臨界車速及容許極限車速均有大幅提高,但隨扣件及道床橫向剛度繼續(xù)增大其增幅逐漸減小,且當(dāng)扣件及道床橫向剛度分別由120 MN/m增至150 MN/m、15 MN/m增至20 MN/m時(shí),FTT系統(tǒng)抗脫軌能力、臨界車速及容許極限車速僅提高了3.9%,1.8%和1.8%.

    3) 隨扣件及道床橫向剛度的增大,鋼軌、軌枕橫向位移均明顯減小;而道床橫向剛度對車體橫向Sperling平穩(wěn)性指標(biāo)的影響則較小.同時(shí),當(dāng)扣件橫向剛度大于90 MN/m及道床橫向剛度大于10 MN/m后,軌道結(jié)構(gòu)橫向位移減小幅度也逐漸變緩.

    綜上所述,考慮日趨緊張的國際重載鐵路市場競爭,建議扣件橫向剛度取90~120 MN/m、道床橫向剛度取10~15 MN/m.另外,本文提出的重載鐵路FTT系統(tǒng)橫向振動穩(wěn)定性分析方法可為重載鐵路軌道剛度合理取值提供一種新的評價(jià)方法.

    References)

    [1]曾慶元,向俊,周志輝,等. 列車脫軌分析理論與應(yīng)用[M]. 長沙:中南大學(xué)出版社,2006: 204-276.

    [2]Xiang Jun, Li Dejian, Zeng Qingyuan. Simulation of spatially coupling dynamic response of train track time variant system [J].JournalofCentralSouthUniversityofTechnology, 2003, 10(3): 226-230. DOI:10.1007/s11771-003-0014-x.

    [3]Zeng Qingyuan, Lou Ping, Xiang Jun. The principle of total potential energy with stationary value in elastic system dynamics and its application to the analysis of vibration and dynamic stability[J].JournalofHuangzhongUniversityofScience&Technology(UrbanScienceEdition), 2002, 19(1): 7-14.

    [4]Lopez-Pita A. The vertical stiffness of the track and the deterioration of high speed lines[J].RevistadeObrasPublicas, 2001, 148(3415): 7-26.

    [5]Real J I, Gomez L, Montalban L. Study of the influence of geometrical and mechanical parameters on ballasted railway tracks design[J].JournalofMechanicalScienceandTechnology, 2012, 26(9): 2837-2844. DOI:10.1007/s12206-012-0734-7.

    [6]Dahlberg T. Railway track stiffness variations—consequences and countermeasures[J].InternationalJournalofCivilEngineering, 2010, 8(1): 1-12.

    [7]Giannakos K. Track maintenance of heavy haul railways with resilient versus stiff fastenings[J].TransportationResearchRecord:JournaloftheTransportationResearchBoard, 2011, 2261: 155-162. DOI:10.3141/2261-18.

    [8]Frohling R D, Scheffel H, Ebersohn W. The vertical dynamic response of a rail vehicle caused by track stiffness variations along the track[J].VehicleSystemDynamics, 1996, 25(1): 175-187.DOI:10.1080/00423119608969194.

    [9]Wu T X, Thompson D J. The effects of local preload on the foundation stiffness and vertical vibration of railway track[J].JournalofSoundandVibration, 1999, 219(5): 881-904. DOI:10.1006/jsvi.1998.1939.

    [10]翟婉明,蔡成標(biāo),王開云. 軌道剛度對列車走行性能的影響[J]. 鐵道學(xué)報(bào),2000, 22(6): 80-83. DOI:10.3321/j.issn:1001-8360.2000.06.018. Zhai Wanming, Cai Chengbiao, Wang Kaiyun. Effect of track stiffness on train running behavior[J].JournaloftheChinaRailwaySociety, 2000, 22(6): 80-83. DOI:10.3321/j.issn:1001-8360.2000.06.018. (in Chinese)

    [11]王開云,翟婉明,蔡成標(biāo). 輪軌結(jié)構(gòu)參數(shù)對列車運(yùn)動穩(wěn)定性的影響[J]. 中國鐵道科學(xué), 2003, 24(1): 43-48. DOI:10.3321/j.issn:1001-4632.2003.01.009.. Wang Kaiyun, Zhai Wanming, Cai Chengbiao. Effect of wheel rail structure parameter on stability of train movement[J].ChinaRailwayScience, 2003, 24(1): 43-48. DOI:10.3321/j.issn:1001-4632.2003.01.009. (in Chinese)

    [12]劉學(xué)毅. 軌道剛度的影響分析及動力學(xué)優(yōu)化[J]. 西南交通大學(xué)學(xué)報(bào),2004, 39(1): 1-5. DOI:10.3969/j.issn.0258-2724.2004.01.001. Liu Xueyi. Effect analysis of track stiffness on dynamic characteristics of wheel-rail system and its dynamic optimization[J].JournalofSouthwestJiaotongUniversity, 2004, 39(1): 1-5. DOI:10.3969/j.issn.0258-2724.2004.01.001. (in Chinese)

    [13]練松良. 軌道結(jié)構(gòu)剛度的理論計(jì)算[J]. 中國鐵道科學(xué),2004, 25(1):67-71. DOI:10.3321/j.issn:1001-4632.2004.01.012. Lian Songliang. Theoretical calculation of the track structure stiffness[J].ChinaRailwayScience, 2004, 25(1):67-71. DOI:10.3321/j.issn:1001-4632.2004.01.012. (in Chinese)

    [14]向俊,曾慶元. 車輛軌道參數(shù)對列車軌道系統(tǒng)振動響應(yīng)的影響分析[J]. 鐵道車輛,2001,39(12):6-8,20. DOI:10.3969/j.issn.1002-7602.2001.12.002. Xiang Jun, Zeng Qingyuan. Analysis on the effect of car, rail parameters on the vibration response of the train-rail system[J].RollingStock, 2001,39(12):6-8,20. DOI:10.3969/j.issn.1002-7602.2001.12.002. (in Chinese)

    [15]向俊,赫丹.板式軌道在高速列車作用下的動力響應(yīng)分析[J]. 工程力學(xué),2008,25(12): 235-240. Xiang Jun, He Dan. Dynamic analysis of slab track under high speed train[J].EngineeringMechanics, 2008,25(12): 235-240. (in Chinese)

    [16]張斌,雷曉燕.基于車輛軌道單元的無砟軌道動力特性有限元分析[J]. 鐵道學(xué)報(bào),2011, 33(7): 78-85. DOI:10.3969/j.issn.1001-8360.2011.07.013. Zhang Bin, Lei Xiaoyan. Analysis on dynamic behavior of ballastless track based on vehicle and track elements with finite element method[J].JournaloftheChinaRailwaySociety, 2011, 33(7): 78-85. DOI:10.3969/j.issn.1001-8360.2011.07.013. (in Chinese)

    [17]王平,徐金輝,汪力,等. 軌道剛度對車輛軌道系統(tǒng)頻率響應(yīng)的影響[J].鐵道工程學(xué)報(bào),2014, 31(9):46-52. DOI:10.3969/j.issn.1006-2106.2014.09.010. Wang Ping, Xu Jinhui, Wang Li, et al. Effect of track stiffness on frequency response of vehicle-track coupling system[J].JournalofRailwayEngineeringSociety, 2014, 31(9):46-52. DOI:10.3969/j.issn.1006-2106.2014.09.010. (in Chinese)

    [18]張格明. 軌道剛度合理值評價(jià)指標(biāo)的研究[J]. 中國鐵道科學(xué),2002, 23(1): 51-57. DOI:10.3321/j.issn:1001-4632.2002.01.009. Zhang Geming. Research on right level of track structure stiffness and track-part stiffness[J].ChinaRailwayScience, 2002, 23(1): 51-57. DOI:10.3321/j.issn:1001-4632.2002.01.009. (in Chinese)

    [19]趙國堂. 鐵路軌道剛度的確定方法[J]. 中國鐵道科學(xué),2005, 26(1): 1-6. DOI:10.3321/j.issn:1001-4632.2005.01.001. Zhao Guotang. Method for determining the rigidity of railway track[J].ChinaRailwayScience, 2005, 26(1): 1-6. DOI:10.3321/j.issn:1001-4632.2005.01.001. (in Chinese)

    [20]向俊,曾慶元. 列車軌道(橋梁)時(shí)變系統(tǒng)橫向振動穩(wěn)定性與失穩(wěn)臨界車速分析方法[J]. 交通運(yùn)輸工程學(xué)報(bào),2011,11(1):19-24. Xiang Jun, Zeng Qingyuan. Analysis method of stability and critical speed for train-track (bridge) time-variable system under lateral vibration condition[J].JournalofTrafficandTransportationEngineering, 2011,11(1):19-24. (in Chinese)

    [21]曾慶元,郭向榮. 列車橋梁時(shí)變系統(tǒng)振動分析理論與應(yīng)用[M]. 北京:中國鐵道出版社,1999: 21-29.

    Effects of track stiffness on freight train anti-derailment safety degree in heavy-haul railway

    Gong Kai1,2Xiang Jun1Yu Cuiying1,3

    (1School of Civil Engineering, Central South University, Changsha 410075, China) (2School of Civil Engineering and Architecture, East China Jiaotong University, Nanchang 330013, China) (3School of Architectural Engineering, Nanchang Institute of Technology, Nanchang 330044, China)

    For the determination of the reasonable range of track stiffness to prevent the freight train derailment, the analysis method for the lateral vibration stability of freight train-track(FTT) system was proposed based on the train-track spatial vibration calculation model and the random energy analysis method for the train derailment. The anti-derailment capacity, the lateral vibration stability and the corresponding vibration response in FTT system on the multiple combination of the fastener and ballast bed lateral stiffnesses were analyzed. Results show that the anti-derailment capacity, critical speed and allowable limit speeds of the FTT system are significantly increased with the increase of the fastener and ballast bed lateral stiffnesses; while the fastener and ballast bed lateral stiffnesses, respectively, increase to 90 MN/m and 10 MN/m, the increased amplitudes decrease. Meanwhile, when the fastener and ballast bed lateral stiffnesses, respectively, increase from 120 MN/m to 150 MN/m, and from 15 MN/m to 20 MN/m, the anti-derailment capacity, critical and allowable limit speeds of the FTT system are increased by only 3.9% 1.8% and 1.8%, respectively. Moreover, increasing the lateral stiffnesses of the fastener and the ballast bed contributes to reducing the lateral displacement of the track structure. In conclusion, considering the rigorous market competition of heavy-haul railway, it is suggested that the fasteners lateral stiffness and the ballast bed lateral stiffness should be from 90 MN/m to 120 MN/m and from 10 MN/m to 15 MN/m, respectively.

    railway track engineering; heavy-haul railway; track stiffness; lateral vibration stability on freight train-track system; freight train

    第47卷第1期2017年1月 東南大學(xué)學(xué)報(bào)(自然科學(xué)版)JOURNALOFSOUTHEASTUNIVERSITY(NaturalScienceEdition) Vol.47No.1Jan.2017DOI:10.3969/j.issn.1001-0505.2017.01.031

    2016-06-14. 作者簡介: 龔凱(1986—),男,博士,講師;向俊(聯(lián)系人),男,博士,教授,博士生導(dǎo)師,jxiang@csu.edu.cn.

    國家自然科學(xué)基金委員會與神華集團(tuán)有限公司聯(lián)合資助項(xiàng)目(U1261113)、高等學(xué)校博士學(xué)科點(diǎn)專項(xiàng)科研基金資助項(xiàng)目(20100162110022)、牽引動力國家重點(diǎn)實(shí)驗(yàn)室開放課題資助項(xiàng)目(TPL0901,TPL1214)、華東交通大學(xué)科研啟動基金資助項(xiàng)目(2003416034)、江西省教育廳科技資助項(xiàng)目(GJJ151173,GJJ151175).

    龔凱,向俊,余翠英.重載鐵路軌道剛度對貨物列車抗脫軌安全度的影響[J].東南大學(xué)學(xué)報(bào)(自然科學(xué)版),2017,47(1):184-192.

    10.3969/j.issn.1001-0505.2017.01.031.

    U213.2

    A

    1001-0505(2017)01-0184-09

    猜你喜歡
    道床軌枕扣件
    裝配式可調(diào)間距的X形抗滑動軌枕研究
    科隆蛋扣件力學(xué)特性分析
    基于滲透性能的道床臟污評估標(biāo)準(zhǔn)研究
    一種改進(jìn)的LBP特征實(shí)現(xiàn)鐵路扣件識別
    CRTS—I型雙塊式無砟軌道到發(fā)線道床板施工方法的改進(jìn)
    科技資訊(2018年15期)2018-10-26 10:59:34
    高速鐵路膠粘道砟固化道床的動力學(xué)特性
    鐵道建筑(2018年2期)2018-03-16 09:36:33
    城市地鐵區(qū)間道床沉降處理施工技術(shù)
    居業(yè)(2016年9期)2016-12-26 21:56:44
    軌枕外觀及修補(bǔ)工藝研究
    非均布荷載下扣件式腳手架的隨機(jī)缺陷穩(wěn)定分析
    諧振式浮軌扣件的安裝工藝
    99国产极品粉嫩在线观看| 熟女电影av网| 国产精品一区二区精品视频观看| 黄色 视频免费看| 日韩精品免费视频一区二区三区| 国产精品av久久久久免费| 久久午夜亚洲精品久久| 中文字幕最新亚洲高清| 又黄又粗又硬又大视频| 亚洲国产精品sss在线观看| av有码第一页| 男插女下体视频免费在线播放| 亚洲欧美精品综合久久99| 在线国产一区二区在线| 老汉色∧v一级毛片| 国产精品久久久人人做人人爽| 窝窝影院91人妻| 无人区码免费观看不卡| 亚洲七黄色美女视频| 免费人成视频x8x8入口观看| 欧美zozozo另类| 男女视频在线观看网站免费 | 久久久国产精品麻豆| 国产亚洲精品第一综合不卡| 三级国产精品欧美在线观看 | 中文字幕精品亚洲无线码一区| 日本 欧美在线| 亚洲aⅴ乱码一区二区在线播放 | 国产激情偷乱视频一区二区| 亚洲天堂国产精品一区在线| 免费在线观看成人毛片| 一本一本综合久久| 一个人免费在线观看的高清视频| 久久99热这里只有精品18| 国产精品av视频在线免费观看| 久久婷婷人人爽人人干人人爱| 日韩欧美一区二区三区在线观看| 两个人免费观看高清视频| 精品久久久久久,| 欧美黑人精品巨大| 免费高清视频大片| 欧美黑人欧美精品刺激| 午夜免费观看网址| 一a级毛片在线观看| 久久久久久免费高清国产稀缺| 欧美一级a爱片免费观看看 | 欧美不卡视频在线免费观看 | www.999成人在线观看| 天堂动漫精品| 午夜福利免费观看在线| 丰满人妻一区二区三区视频av | 在线观看www视频免费| 日韩精品中文字幕看吧| 亚洲av中文字字幕乱码综合| 日韩大尺度精品在线看网址| 久久 成人 亚洲| 亚洲成av人片在线播放无| 老司机深夜福利视频在线观看| 日韩精品中文字幕看吧| 国产99白浆流出| 99久久99久久久精品蜜桃| 国产黄色小视频在线观看| 中亚洲国语对白在线视频| 国内毛片毛片毛片毛片毛片| 又大又爽又粗| 成年免费大片在线观看| 精品第一国产精品| 国产成人av激情在线播放| 777久久人妻少妇嫩草av网站| 又粗又爽又猛毛片免费看| 一区福利在线观看| 老汉色av国产亚洲站长工具| 亚洲九九香蕉| 久久久久国内视频| 一级作爱视频免费观看| 天堂影院成人在线观看| 搞女人的毛片| 成人特级黄色片久久久久久久| 最好的美女福利视频网| 久久久久久久久久黄片| 国内毛片毛片毛片毛片毛片| 午夜免费激情av| 国产麻豆成人av免费视频| 日韩欧美国产在线观看| 真人做人爱边吃奶动态| 日日摸夜夜添夜夜添小说| 91九色精品人成在线观看| 一级黄色大片毛片| 神马国产精品三级电影在线观看 | 可以在线观看的亚洲视频| 精品一区二区三区四区五区乱码| 成熟少妇高潮喷水视频| 国产在线精品亚洲第一网站| 91九色精品人成在线观看| 日韩中文字幕欧美一区二区| 亚洲精品一区av在线观看| 国产免费av片在线观看野外av| av在线播放免费不卡| 国产不卡一卡二| 亚洲成人精品中文字幕电影| 嫩草影视91久久| 国产欧美日韩精品亚洲av| 中文字幕熟女人妻在线| 麻豆av在线久日| 亚洲一码二码三码区别大吗| 2021天堂中文幕一二区在线观| 村上凉子中文字幕在线| 色噜噜av男人的天堂激情| 男女做爰动态图高潮gif福利片| www.999成人在线观看| 亚洲中文av在线| av福利片在线| 日本撒尿小便嘘嘘汇集6| 丰满人妻熟妇乱又伦精品不卡| 中文字幕久久专区| 又黄又粗又硬又大视频| 法律面前人人平等表现在哪些方面| 亚洲国产欧美一区二区综合| 1024手机看黄色片| 久久久国产成人精品二区| 日本一本二区三区精品| 在线十欧美十亚洲十日本专区| 在线观看www视频免费| 国产在线观看jvid| 国产午夜精品久久久久久| 69av精品久久久久久| 欧美日韩福利视频一区二区| xxx96com| 国产1区2区3区精品| 国产精品国产高清国产av| 国产精品久久久久久久电影 | 亚洲一码二码三码区别大吗| 久久性视频一级片| 操出白浆在线播放| 亚洲欧美精品综合久久99| 99精品在免费线老司机午夜| 99久久99久久久精品蜜桃| 曰老女人黄片| 老熟妇仑乱视频hdxx| 丰满人妻一区二区三区视频av | 九色国产91popny在线| 日韩欧美在线乱码| 一区二区三区激情视频| 精华霜和精华液先用哪个| 99在线人妻在线中文字幕| 最近在线观看免费完整版| 国产在线精品亚洲第一网站| 18禁观看日本| 国产成人精品久久二区二区91| 国产精品永久免费网站| 97碰自拍视频| 国产精品久久视频播放| 国产成人精品久久二区二区免费| 日韩精品免费视频一区二区三区| 日韩有码中文字幕| 成年版毛片免费区| 亚洲熟妇熟女久久| www国产在线视频色| av免费在线观看网站| 在线免费观看的www视频| 成人av在线播放网站| 欧美乱码精品一区二区三区| 亚洲激情在线av| 久9热在线精品视频| 亚洲av片天天在线观看| 亚洲国产精品sss在线观看| 色精品久久人妻99蜜桃| 色综合婷婷激情| 美女免费视频网站| 国产亚洲av嫩草精品影院| 中文字幕av在线有码专区| 国产成人系列免费观看| 色精品久久人妻99蜜桃| 黑人操中国人逼视频| 国产一区在线观看成人免费| www.www免费av| 国产1区2区3区精品| 精品久久久久久,| 免费看十八禁软件| 欧美日本视频| 美女高潮喷水抽搐中文字幕| 免费看日本二区| 国产亚洲av高清不卡| 麻豆国产97在线/欧美 | 村上凉子中文字幕在线| 老司机午夜十八禁免费视频| 成人永久免费在线观看视频| 欧美黑人欧美精品刺激| 人妻久久中文字幕网| 国产99白浆流出| 精品国产乱码久久久久久男人| √禁漫天堂资源中文www| 国产一级毛片七仙女欲春2| 亚洲人成77777在线视频| 亚洲精品国产一区二区精华液| 成人特级黄色片久久久久久久| 一个人免费在线观看的高清视频| 久久国产精品人妻蜜桃| 美女免费视频网站| 99精品久久久久人妻精品| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲av日韩精品久久久久久密| 久久精品国产99精品国产亚洲性色| 麻豆成人av在线观看| 亚洲av日韩精品久久久久久密| 色av中文字幕| 日韩欧美国产一区二区入口| 99riav亚洲国产免费| 午夜日韩欧美国产| 少妇被粗大的猛进出69影院| 欧美黄色淫秽网站| 精品国产超薄肉色丝袜足j| 舔av片在线| 在线播放国产精品三级| 无遮挡黄片免费观看| 久久久久九九精品影院| 亚洲国产精品合色在线| 亚洲av电影在线进入| 午夜老司机福利片| 丰满人妻熟妇乱又伦精品不卡| 午夜免费激情av| 成人三级做爰电影| 757午夜福利合集在线观看| 久久久久久人人人人人| 国产精品av视频在线免费观看| 舔av片在线| 美女 人体艺术 gogo| 久久热在线av| 曰老女人黄片| 伊人久久大香线蕉亚洲五| 女生性感内裤真人,穿戴方法视频| 丝袜人妻中文字幕| 国产精品久久电影中文字幕| 国产精品亚洲美女久久久| 国产成人精品无人区| 桃色一区二区三区在线观看| 天堂av国产一区二区熟女人妻 | 搡老岳熟女国产| 亚洲无线在线观看| a在线观看视频网站| 亚洲国产精品久久男人天堂| 亚洲专区字幕在线| 久久久久久人人人人人| 久久久久久久午夜电影| 很黄的视频免费| 国产熟女xx| 日韩成人在线观看一区二区三区| 91av网站免费观看| 香蕉国产在线看| 熟女少妇亚洲综合色aaa.| 窝窝影院91人妻| 亚洲人成网站在线播放欧美日韩| 国产精品久久久av美女十八| 51午夜福利影视在线观看| 欧美日韩亚洲国产一区二区在线观看| 看片在线看免费视频| 身体一侧抽搐| 五月玫瑰六月丁香| 久久精品成人免费网站| 99热这里只有精品一区 | 狂野欧美激情性xxxx| 国产亚洲精品一区二区www| 午夜免费成人在线视频| 十八禁网站免费在线| 后天国语完整版免费观看| 午夜久久久久精精品| 亚洲专区国产一区二区| 国产91精品成人一区二区三区| 禁无遮挡网站| 黄色毛片三级朝国网站| 午夜影院日韩av| 久久亚洲真实| 亚洲欧美精品综合久久99| 国产麻豆成人av免费视频| 日韩欧美在线乱码| 男人舔奶头视频| 嫩草影院精品99| 性欧美人与动物交配| 免费一级毛片在线播放高清视频| 欧美另类亚洲清纯唯美| 91九色精品人成在线观看| www.自偷自拍.com| 伊人久久大香线蕉亚洲五| 色综合站精品国产| 很黄的视频免费| 日日夜夜操网爽| 三级国产精品欧美在线观看 | 老熟妇乱子伦视频在线观看| 天天躁狠狠躁夜夜躁狠狠躁| 精品国产乱码久久久久久男人| 国产v大片淫在线免费观看| 国产免费男女视频| 在线观看一区二区三区| 一个人观看的视频www高清免费观看 | 国产三级中文精品| 在线观看美女被高潮喷水网站 | 一本一本综合久久| 成熟少妇高潮喷水视频| 国产免费男女视频| 久久久久亚洲av毛片大全| 黄色视频不卡| 可以在线观看的亚洲视频| 国产午夜精品论理片| 成人三级黄色视频| 国产伦人伦偷精品视频| 老司机福利观看| 又黄又爽又免费观看的视频| 午夜免费观看网址| 在线观看免费午夜福利视频| 搡老熟女国产l中国老女人| 国产又色又爽无遮挡免费看| 狠狠狠狠99中文字幕| 在线看三级毛片| 成人精品一区二区免费| 国产视频一区二区在线看| 丰满人妻一区二区三区视频av | 国产精品一区二区免费欧美| 欧美精品亚洲一区二区| 亚洲,欧美精品.| 99精品久久久久人妻精品| 香蕉久久夜色| 日本 av在线| 亚洲国产精品成人综合色| 91麻豆av在线| 日本熟妇午夜| av国产免费在线观看| 嫩草影院精品99| 狂野欧美激情性xxxx| 午夜激情av网站| 国内揄拍国产精品人妻在线| 男女视频在线观看网站免费 | 淫秽高清视频在线观看| 日日爽夜夜爽网站| 变态另类丝袜制服| 日日夜夜操网爽| 欧美黑人欧美精品刺激| 精品无人区乱码1区二区| 人人妻人人看人人澡| 国产精品1区2区在线观看.| 久久这里只有精品中国| 亚洲,欧美精品.| 国产av又大| 精品第一国产精品| 午夜视频精品福利| 国产午夜精品久久久久久| 国产视频内射| 日韩精品中文字幕看吧| 亚洲国产欧洲综合997久久,| 女人被狂操c到高潮| 国产午夜精品论理片| 亚洲美女黄片视频| 一夜夜www| 精品无人区乱码1区二区| 好看av亚洲va欧美ⅴa在| 亚洲美女视频黄频| 18禁裸乳无遮挡免费网站照片| 久久久久九九精品影院| 久久国产精品人妻蜜桃| 人妻夜夜爽99麻豆av| 久久久久久大精品| 日韩欧美在线二视频| 麻豆久久精品国产亚洲av| 嫩草影院精品99| 男女那种视频在线观看| 国产野战对白在线观看| 成年人黄色毛片网站| 欧美 亚洲 国产 日韩一| 在线永久观看黄色视频| 午夜福利视频1000在线观看| 观看免费一级毛片| 三级国产精品欧美在线观看 | av天堂在线播放| a级毛片在线看网站| 亚洲欧美日韩高清在线视频| 免费看a级黄色片| 久久九九热精品免费| 中文亚洲av片在线观看爽| 高清毛片免费观看视频网站| 久久婷婷成人综合色麻豆| 真人一进一出gif抽搐免费| 我的老师免费观看完整版| 亚洲一卡2卡3卡4卡5卡精品中文| 午夜两性在线视频| 51午夜福利影视在线观看| 精品国内亚洲2022精品成人| 午夜福利18| 日韩高清综合在线| 50天的宝宝边吃奶边哭怎么回事| 人人妻,人人澡人人爽秒播| 国产伦人伦偷精品视频| 久久热在线av| 亚洲欧美日韩无卡精品| 69av精品久久久久久| 免费观看精品视频网站| 一级a爱片免费观看的视频| 久久精品成人免费网站| 国产亚洲精品一区二区www| 国产一区二区在线观看日韩 | 手机成人av网站| 国产精品九九99| 日韩欧美在线二视频| 久久精品人妻少妇| 日韩欧美免费精品| 国产人伦9x9x在线观看| 亚洲成人久久性| 午夜亚洲福利在线播放| 丰满人妻熟妇乱又伦精品不卡| 精品国产超薄肉色丝袜足j| 香蕉丝袜av| 日韩大尺度精品在线看网址| 村上凉子中文字幕在线| 久久午夜亚洲精品久久| 老司机深夜福利视频在线观看| 91成年电影在线观看| 国产亚洲精品一区二区www| 久久香蕉国产精品| 级片在线观看| 色综合亚洲欧美另类图片| 国产亚洲欧美98| 黄色毛片三级朝国网站| 日韩精品免费视频一区二区三区| 18禁美女被吸乳视频| 成人18禁在线播放| 少妇熟女aⅴ在线视频| 美女大奶头视频| 亚洲国产欧美一区二区综合| av片东京热男人的天堂| 亚洲va日本ⅴa欧美va伊人久久| 国产精品一及| 成年人黄色毛片网站| 国产熟女xx| 亚洲精品美女久久久久99蜜臀| x7x7x7水蜜桃| 搞女人的毛片| 色综合站精品国产| 欧美丝袜亚洲另类 | 欧美高清成人免费视频www| 最好的美女福利视频网| 免费搜索国产男女视频| 麻豆一二三区av精品| 他把我摸到了高潮在线观看| 成人亚洲精品av一区二区| 亚洲人成网站在线播放欧美日韩| 国产黄色小视频在线观看| 国产精品久久久人人做人人爽| 中文亚洲av片在线观看爽| 88av欧美| 日本撒尿小便嘘嘘汇集6| 中文亚洲av片在线观看爽| 淫秽高清视频在线观看| 欧美极品一区二区三区四区| 日韩精品青青久久久久久| av在线天堂中文字幕| 欧美久久黑人一区二区| 无限看片的www在线观看| 悠悠久久av| 亚洲成人国产一区在线观看| 十八禁人妻一区二区| 一本精品99久久精品77| 亚洲真实伦在线观看| 1024手机看黄色片| 无限看片的www在线观看| 在线视频色国产色| 亚洲成人久久爱视频| 亚洲国产精品合色在线| 丝袜美腿诱惑在线| 精品久久久久久成人av| 麻豆久久精品国产亚洲av| 麻豆av在线久日| 精品高清国产在线一区| 日韩欧美一区二区三区在线观看| 禁无遮挡网站| 久久人妻av系列| 好看av亚洲va欧美ⅴa在| 2021天堂中文幕一二区在线观| 一个人观看的视频www高清免费观看 | 国产又黄又爽又无遮挡在线| 999久久久精品免费观看国产| 丝袜美腿诱惑在线| 脱女人内裤的视频| 精华霜和精华液先用哪个| 这个男人来自地球电影免费观看| 一级毛片精品| 黄色视频不卡| 精品国内亚洲2022精品成人| 一个人免费在线观看电影 | 动漫黄色视频在线观看| 午夜福利视频1000在线观看| 99国产极品粉嫩在线观看| 三级毛片av免费| 国产精品香港三级国产av潘金莲| 久久久久久亚洲精品国产蜜桃av| 1024香蕉在线观看| 亚洲中文日韩欧美视频| 国产精品精品国产色婷婷| 国产精品1区2区在线观看.| 免费人成视频x8x8入口观看| 国产高清视频在线播放一区| 午夜精品在线福利| 亚洲自拍偷在线| 成年人黄色毛片网站| 亚洲 欧美 日韩 在线 免费| 日本撒尿小便嘘嘘汇集6| 亚洲精品在线观看二区| 非洲黑人性xxxx精品又粗又长| 亚洲精品一卡2卡三卡4卡5卡| 亚洲五月天丁香| 999久久久国产精品视频| 色av中文字幕| 人人妻,人人澡人人爽秒播| 很黄的视频免费| 亚洲人成网站在线播放欧美日韩| 久久精品国产清高在天天线| 精品久久久久久成人av| 麻豆久久精品国产亚洲av| 亚洲中文av在线| 亚洲一区高清亚洲精品| 全区人妻精品视频| 久久精品国产亚洲av高清一级| 色av中文字幕| 国产精品久久久久久久电影 | 国产日本99.免费观看| 国产伦在线观看视频一区| 露出奶头的视频| 日本一本二区三区精品| 国产精品一及| 欧美一级毛片孕妇| 亚洲美女黄片视频| 99在线视频只有这里精品首页| 少妇粗大呻吟视频| 亚洲av熟女| 精品国内亚洲2022精品成人| 99在线人妻在线中文字幕| 午夜福利成人在线免费观看| 亚洲精品在线观看二区| 日韩av在线大香蕉| 精品一区二区三区av网在线观看| 欧美黑人精品巨大| 国产精品美女特级片免费视频播放器 | 国产成年人精品一区二区| 熟女少妇亚洲综合色aaa.| 亚洲国产看品久久| 中文字幕精品亚洲无线码一区| 久久精品人妻少妇| 嫁个100分男人电影在线观看| 久久国产精品影院| 免费观看人在逋| 午夜成年电影在线免费观看| 国产高清激情床上av| 18禁黄网站禁片午夜丰满| 一进一出抽搐gif免费好疼| 成人午夜高清在线视频| 亚洲中文日韩欧美视频| 欧美av亚洲av综合av国产av| 精品久久久久久久末码| 婷婷精品国产亚洲av| 最近视频中文字幕2019在线8| 欧美极品一区二区三区四区| 国产真人三级小视频在线观看| 麻豆成人午夜福利视频| 成人av一区二区三区在线看| 日韩大码丰满熟妇| 亚洲 国产 在线| 亚洲人成电影免费在线| 成人国产一区最新在线观看| 男人的好看免费观看在线视频 | 999久久久国产精品视频| 91av网站免费观看| 中文字幕最新亚洲高清| 亚洲成人中文字幕在线播放| 99精品久久久久人妻精品| 国产精品综合久久久久久久免费| 免费看十八禁软件| 国产在线精品亚洲第一网站| 男女床上黄色一级片免费看| www.精华液| 99久久精品热视频| 精品午夜福利视频在线观看一区| 国产一区二区在线av高清观看| 亚洲午夜理论影院| 久久九九热精品免费| 在线观看www视频免费| 欧美成人免费av一区二区三区| 欧美一区二区国产精品久久精品 | 这个男人来自地球电影免费观看| 欧美日本视频| 亚洲九九香蕉| 黑人欧美特级aaaaaa片| 亚洲真实伦在线观看| xxxwww97欧美| 999精品在线视频| 一级作爱视频免费观看| 国产视频一区二区在线看| 在线观看免费日韩欧美大片| 日韩精品中文字幕看吧| 精品乱码久久久久久99久播| 国产成人啪精品午夜网站| 日韩欧美一区二区三区在线观看| 欧美成人免费av一区二区三区| 国内揄拍国产精品人妻在线| 亚洲天堂国产精品一区在线| 久久精品国产清高在天天线| 精品日产1卡2卡| 精品久久久久久,| 操出白浆在线播放| 国产三级在线视频| 夜夜夜夜夜久久久久| 国产一区在线观看成人免费| aaaaa片日本免费| 欧美日本亚洲视频在线播放| 非洲黑人性xxxx精品又粗又长| 国产野战对白在线观看| 在线看三级毛片| 国产av一区二区精品久久| 欧美一级a爱片免费观看看 |