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

    基于輪軌時(shí)變接觸剛度的車(chē)輛?軌道系統(tǒng)垂向振動(dòng)分析

    2017-11-01 14:18:45賀天龍陳憲麥徐磊
    關(guān)鍵詞:密貼平順輪軌

    賀天龍,陳憲麥,徐磊

    ?

    基于輪軌時(shí)變接觸剛度的車(chē)輛?軌道系統(tǒng)垂向振動(dòng)分析

    賀天龍1, 2,陳憲麥1, 2,徐磊3

    (1. 中南大學(xué)土木工程學(xué)院,湖南長(zhǎng)沙,410075;2. 高速鐵路建造技術(shù)國(guó)家工程實(shí)驗(yàn)室,湖南長(zhǎng)沙,410075;3. 西南交通大學(xué)牽引動(dòng)力國(guó)家重點(diǎn)實(shí)驗(yàn)室,四川成都,610031)

    采用輪軌時(shí)變接觸剛度代替輪軌非線性接觸力的動(dòng)力效應(yīng),模擬輪軌之間的自然接觸狀態(tài),將有限元理論和能量變分法用于建立車(chē)輛?軌道系統(tǒng)垂向統(tǒng)一方程。通過(guò)輸入脈沖型短波及中、長(zhǎng)波不平順激勵(lì),計(jì)算系統(tǒng)的動(dòng)力響應(yīng),驗(yàn)證本文模型的正確性,并與國(guó)內(nèi)常用的輪軌密貼模型進(jìn)行比較分析。研究結(jié)果表明:輪軌密貼模型在分析中、長(zhǎng)波不平順激勵(lì)工況下的系統(tǒng)動(dòng)力響應(yīng)是可行的;但在短波不平順激勵(lì)工況下,輪軌之間的剛性處理方法將放大短波激勵(lì)效應(yīng),車(chē)輪“跳軌”處理方法尚需改進(jìn)。

    車(chē)輛?軌道耦合系統(tǒng);時(shí)變接觸剛度;垂向振動(dòng);能量變分法;統(tǒng)一方程

    車(chē)輛?軌道系統(tǒng)的耦合振動(dòng)分析是鐵路動(dòng)力學(xué)領(lǐng)域非常重要的研究?jī)?nèi)容,國(guó)內(nèi)的車(chē)/軌系統(tǒng)動(dòng)力學(xué)建模工作以車(chē)輛?軌道耦合動(dòng)力學(xué)[1]為基本理論指導(dǎo),已日趨完善。針對(duì)車(chē)輛?軌道系統(tǒng)的橫向相互作用,應(yīng)該考慮非線性蠕滑力及輪軌橫向相對(duì)位移。然而,對(duì)于車(chē)輛?軌道系統(tǒng)的垂向相互作用,國(guó)內(nèi)主要采用輪軌密貼[2?4]和輪軌彈性接觸(基于Hertz非線性接觸理論)[1]這2種建模方法。由于鐵路鋼軌存在的磨耗、扣件失效、軌道板離縫及路基沉降等一系列鐵路軌道基礎(chǔ)結(jié)構(gòu)動(dòng)態(tài)性能演變問(wèn)題都與車(chē)/軌系統(tǒng)的垂向動(dòng)力作用密切相關(guān),所以,有必要對(duì)這2種垂向輪軌關(guān)系作用下的系統(tǒng)動(dòng)力響應(yīng)進(jìn)行比較分析。為此,本文作者基于有限元方法對(duì)不同鐵路線路工況的強(qiáng)模擬性能及求解軌道結(jié)構(gòu)任意位置的動(dòng)力響應(yīng),采用有限元方法,同時(shí)考慮輪軌之間的彈性接觸剛度(即考慮輪軌相對(duì)位移),以板式軌道結(jié)構(gòu)為例,對(duì)這2種垂向輪軌關(guān)系進(jìn)行計(jì)算。實(shí)際上,人們對(duì)軌道的有限元?jiǎng)恿M(jìn)行了大量研究,如:羅震等[5]基于多剛體系統(tǒng)動(dòng)力學(xué)原理建立了車(chē)輛動(dòng)力學(xué)計(jì)算模型,以梁?板?板有限元模擬板式軌道,通過(guò)輪軌界面形成車(chē)/軌系統(tǒng)的動(dòng)態(tài)耦合,分別采用新型快速顯示積分方法[6]和Newmark-β法[7]求解了車(chē)輛系統(tǒng)和軌道系統(tǒng)的動(dòng)力響應(yīng);向俊等[8]采用橫向有限條和板段單元模型分析了車(chē)/軌垂向動(dòng)力響應(yīng),在正常行車(chē)時(shí)為輪軌密貼,而在車(chē)輪跳軌時(shí),考慮輪軌相對(duì)位移。此外,人們?cè)诳紤]輪軌接觸彈性時(shí),常將輪軌非線性接觸彈簧時(shí)不變等效成線性化處理[9?12]。這些方法均沒(méi)有在輪軌接觸非線性、有限元方法及動(dòng)力系統(tǒng)方程的統(tǒng)一求解上達(dá)到統(tǒng)一。本文作者采用輪/軌時(shí)變接觸剛度的方法模擬輪軌之間的自然接觸狀態(tài),通過(guò)能量變分和有限元建立車(chē)輛?軌道系統(tǒng)的統(tǒng)一方程,驗(yàn)證模型的正確性。建立能考慮車(chē)輪“跳軌”的輪軌密貼模型,比較本文模型與輪軌密貼模型在脈沖型短波及中、長(zhǎng)波不平順激勵(lì)下的系統(tǒng)動(dòng)力響應(yīng)。

    1 垂向耦合動(dòng)力學(xué)模型

    將車(chē)輛?板式軌道作為一個(gè)整體,建立如圖1所示的車(chē)輛?板式軌道垂向耦合動(dòng)力學(xué)模型。高速車(chē)輛模型采用具有二系懸掛的車(chē)輛模型。板式軌道因其結(jié)構(gòu)左右對(duì)稱,可取一股軌道進(jìn)行研究,以鋼軌模擬無(wú)限長(zhǎng)且離散黏彈性點(diǎn)支承的Bernoulli?Euler梁,彈性支承在鋼軌扣件處的軌下襯墊上;將軌道板模擬成連續(xù)黏彈性支承的有限長(zhǎng)Bernoulli-Euler梁,CA砂漿等效為連續(xù)分布線性彈簧和阻尼器。

    圖1 車(chē)輛?板式軌道垂向耦合動(dòng)力學(xué)模型

    2 垂向耦合振動(dòng)方程的建立

    對(duì)于車(chē)/軌系統(tǒng)的垂向或橫向振動(dòng),輪軌相互作用力確是車(chē)輛?軌道系統(tǒng)的內(nèi)力,翟婉明[1]將每一積分步的輪軌力“顯示”表達(dá),分別將其作用于車(chē)輛系統(tǒng)與軌道系統(tǒng),獲得了與實(shí)測(cè)值接近的結(jié)果,得到廣泛應(yīng)用。曾慶元[2]認(rèn)為輪軌相互作用力是“隱式”的,由于輪軌作用力與輪軌相對(duì)位移存在強(qiáng)非線性關(guān)系,對(duì)輪軌作用力在相對(duì)位移上的做功無(wú)法有效地進(jìn)行一階變分,只能采用輪軌密貼及線性蠕滑的假定。

    在車(chē)?軌(橋)系統(tǒng)的數(shù)值積分步長(zhǎng)內(nèi),所有的力及動(dòng)力參數(shù)是基本不變的。本文假定輪軌之間存在1個(gè)輪/軌接觸彈簧,隨著積分步變化,其接觸彈簧剛度是時(shí)變的,但在1個(gè)積分步長(zhǎng)內(nèi),認(rèn)為此彈簧接觸剛度為定值,這符合輪軌之間的真實(shí)接觸狀態(tài)。

    2.1 輪軌接觸關(guān)系的處理方法

    2.1.1 輪軌時(shí)變接觸剛度模型

    根據(jù)Hertz非線性彈性接觸理論,可以確定輪軌之間的相互作用力()為

    式中:為輪軌接觸常數(shù)(m/N2/3);δ()為輪軌間的彈性壓縮量(m)。采用磨耗型踏面車(chē)輪,則

    圖2 輪軌時(shí)變接觸剛度與輪軌相對(duì)位移kwr(t)的關(guān)系

    輪軌時(shí)變接觸剛度的計(jì)算基于Hertz非線性接觸理論,這種方法完全將車(chē)輛系統(tǒng)與軌道系統(tǒng)耦合成一個(gè)整體,輪軌之間的彈性變形能可由輪軌接觸剛度和輪軌相對(duì)位移計(jì)算獲得。這是一種輪軌力“隱式”表達(dá)的方法,輪軌力作為輪軌之間的相互作用力是車(chē)輛?軌道系統(tǒng)的1個(gè)內(nèi)力。當(dāng)采用能量原理建立系統(tǒng)的平衡方程時(shí),只需將輪軌變形能計(jì)入系統(tǒng)總能量即可,這兼顧了輪軌非線性接觸本質(zhì)和振動(dòng)方程解適定性的要求[13]。

    2.1.2 輪軌密貼模型

    曾慶元等[2?3]采用的模型均為輪軌密貼模型。曾慶元[2]認(rèn)為軌道不平順不是時(shí)間的函數(shù),軌道不平順附加速度及加速度均為0 m/s2;婁平等[3]采用較復(fù)雜的輪軌約束方程,不僅考慮了軌道不平順的附加速度及加速度項(xiàng),而且加入了形函數(shù)矩陣的時(shí)間求導(dǎo)項(xiàng)。本文在文獻(xiàn)[3]的基礎(chǔ)上,采用簡(jiǎn)化的輪軌約束方程:

    輪軌密貼模型在時(shí)域數(shù)值積分過(guò)程中,輪軌力是不計(jì)算和代入系統(tǒng)動(dòng)力方程荷載列陣的。其輪軌力可根據(jù)每一步的計(jì)算結(jié)果,基于達(dá)朗貝爾原理進(jìn)行計(jì)算:

    式中:w1,w2,w3和w4分別為第1~4輪/軌接觸力;1,2,3和4分別為第1~4輪對(duì)的靜輪載;其他變量含義見(jiàn)文獻(xiàn)[1]。

    2.2 車(chē)輛和軌道系統(tǒng)矩陣的建立

    2.2.1 Hamilton原理

    Hamilton原理是力學(xué)中應(yīng)用最廣泛和最重要的積分形式的變分原理,提供了從所有可能運(yùn)動(dòng)中找出真實(shí)運(yùn)動(dòng)的1個(gè)準(zhǔn)則。對(duì)于完整的保守系統(tǒng),Hamilton最小作用量原理定義為[14]

    式中:為拉格朗日泛函。實(shí)際上,從Hamilton最小作用量原理成立的基本條件來(lái)說(shuō),由于車(chē)/軌系統(tǒng)存在黏滯阻尼力、庫(kù)侖摩擦力及干擾力等非有勢(shì)力,不再符合上述Hamilton完整保守系統(tǒng)的定義。但當(dāng)積分時(shí)間足夠小(1→0)時(shí),進(jìn)行瞬時(shí)處理。在某一瞬時(shí),所有的內(nèi)外力及應(yīng)力是不變的,只產(chǎn)生滿足勢(shì)能駐值原理[15]的位移及應(yīng)變,此時(shí),非有勢(shì)力(黏滯阻尼力、摩擦力及干擾力)可近似作有勢(shì)力處理。曾慶元[16]基于這種瞬時(shí)狀態(tài)的理想化處理,建立了彈性系統(tǒng)動(dòng)力學(xué)總勢(shì)能不變值原理,在車(chē)輛?軌道(橋梁)振動(dòng)分析領(lǐng)域獲得了廣泛應(yīng)用。這里的拉格朗日泛函可為

    式中:M為系統(tǒng)的動(dòng)能;I為系統(tǒng)的彈性勢(shì)能;C為系統(tǒng)的阻尼力勢(shì)能;F為外力勢(shì)能。

    由式(3)可得

    根據(jù)式(9)的基本表達(dá)形式及“對(duì)號(hào)入座”法則[2],可以建立車(chē)輛和軌道的系統(tǒng)矩陣。

    2.2.2 車(chē)輛系統(tǒng)矩陣

    車(chē)輛系統(tǒng)的動(dòng)力矩陣可參考文獻(xiàn)[17]中的方法獲得。這里以車(chē)輛系統(tǒng)剛度矩陣為例,簡(jiǎn)單說(shuō)明其形成過(guò)程。在車(chē)輛-軌道系統(tǒng)的時(shí)域數(shù)值積分過(guò)程中,將輪軌之間的接觸彈性變形能在車(chē)輛系統(tǒng)矩陣的分析過(guò)程中予以考慮。

    車(chē)輛系統(tǒng)的初始振動(dòng)位移參數(shù)列陣如下:

    式中:c,gq和gh分別為車(chē)體、前構(gòu)架、后構(gòu)架的沉浮自由度;c,gq和gh分別為車(chē)體、前構(gòu)架、后構(gòu)架的點(diǎn)頭自由度;w1,w2,w3和w4為第1~4輪對(duì)的沉浮自由度;r1,r2,r3和r4為第1~4輪對(duì)對(duì)應(yīng)的鋼軌振動(dòng)位移;irr1,irr2,irr3和irr4為第1~4輪對(duì)對(duì)應(yīng)的軌道垂向不平順。

    需指出的是:鋼軌振動(dòng)位移及軌道不平順在原則上是屬于軌道系統(tǒng)的振動(dòng)位移參數(shù),本文為編制程序方便,將輪軌接觸點(diǎn)處對(duì)應(yīng)輪對(duì)位移、鋼軌位移及鋼軌不平順均納入車(chē)輛系統(tǒng)中。

    對(duì)于多剛體的車(chē)輛系統(tǒng),其剛度矩陣的形成主要通過(guò)彈簧單元變形能的位移變分獲得,故本文的車(chē)輛系統(tǒng)共10彈簧變形單元,如車(chē)體與前構(gòu)架、前構(gòu)架與第1輪對(duì)、第1輪/軌接觸彈簧的彈性變形能(以彈簧受壓為正)分別為

    式中:s1,gq1和wr1分別為車(chē)體與前構(gòu)架、前構(gòu)架與第1輪對(duì)、第1輪/軌接觸點(diǎn)之間的彈簧剛度;c和t分別為車(chē)輛二構(gòu)架中心距之一半、同一構(gòu)架所屬二輪對(duì)軸距之一半。對(duì)式(18)進(jìn)行位移變分,得

    前構(gòu)架與第1輪對(duì)之間的彈簧剛度矩陣和第1輪軌接觸對(duì)的彈簧剛度矩陣可按類似方法獲得。

    仿照上述方法,根據(jù)位移參數(shù)在矩陣中的位置,可以組裝車(chē)輛系統(tǒng)的總剛度矩陣。質(zhì)量矩陣與阻尼矩陣的形成與此類似。

    2.2.3 軌道系統(tǒng)矩陣

    取2個(gè)扣件之間的軌段為軌道系統(tǒng)的1個(gè)有限單元,若不考慮軌道板與軌道板之間的縫隙,則無(wú)論多長(zhǎng)的軌道系統(tǒng)矩陣,均可以基于此有限單元的質(zhì)量、剛度、阻尼矩陣快速組裝形成整個(gè)軌道系統(tǒng)的動(dòng)力矩陣。

    由于考慮輪軌之間的時(shí)變接觸剛度,在整個(gè)數(shù)值積分計(jì)算過(guò)程中,車(chē)/軌系統(tǒng)矩陣是非線性變化的,此時(shí),有限單元?jiǎng)澐珠L(zhǎng)度需要滿足一定的精細(xì)化要求,否則將得不到正確的計(jì)算結(jié)果。

    梁?jiǎn)卧獜澢魏瘮?shù)采用Hermitian三次插值函數(shù),梁?jiǎn)卧獌?nèi)任意一點(diǎn)的豎向位移都可以通過(guò)節(jié)點(diǎn)的豎向位移和轉(zhuǎn)角位移表示:

    軌道系統(tǒng)矩陣的形成見(jiàn)文獻(xiàn)[18?19]。

    2.2.4 車(chē)輛系統(tǒng)和軌道系統(tǒng)的耦合

    與文獻(xiàn)[1]中每一步顯示求解輪軌力的方法不同,本文的車(chē)輛系統(tǒng)與軌道系統(tǒng)的垂向耦合通過(guò)輪/軌之間的時(shí)變剛度系數(shù)分配的方法實(shí)現(xiàn)。通過(guò)對(duì)第1~4輪軌接觸彈簧變形能的位移變分,可獲得輪軌之間彈簧剛度矩陣。如第1輪軌接觸對(duì)的彈性變形能變分為

    經(jīng)過(guò)重分配后的第1輪軌接觸剛度矩陣為

    通過(guò)式(23)~(25)實(shí)現(xiàn)了輪對(duì)垂向振動(dòng)位移、輪對(duì)所在的鋼軌接觸單元節(jié)點(diǎn)位移及輪/軌接觸點(diǎn)的軌道不平順等位移參數(shù)的銜接。其余3個(gè)輪/軌接觸對(duì)的剛度矩陣可按類似方法獲得。這樣,便解決了車(chē)輛系統(tǒng)和軌道系統(tǒng)的垂向耦合問(wèn)題。

    需要注意的是:式(26)并非最終的車(chē)輛?軌道耦合振動(dòng)方程,因?yàn)轳詈舷到y(tǒng)的剛度矩陣中含有非獨(dú)立位移參數(shù)的彈性變形能,如式(17)中的后8個(gè)振動(dòng)位移參數(shù)。應(yīng)將軌道不平順實(shí)測(cè)或反演值代入式(26),形成車(chē)輛?軌道系統(tǒng)的內(nèi)部激勵(lì),然后劃去非獨(dú)立矩陣方程,形成最終的耦合振動(dòng)方程,其具體原理見(jiàn)文獻(xiàn)[2]。本文采用wilson?法求解此振動(dòng)方程,積分步長(zhǎng)為0.05 ms。

    采用較小的積分步長(zhǎng)及輪軌密貼的方法,相當(dāng)于給定了輪軌之間的約束條件,在較大的時(shí)間步長(zhǎng)(如0.001~0.010 s)[2?3]時(shí)也能得到收斂解。而采用輪軌時(shí)變接觸剛度完全模擬輪軌之間的自然接觸狀態(tài),在時(shí)域數(shù)值積分時(shí),每一積分步內(nèi)所有的力素是不變的,而實(shí)際上此步長(zhǎng)內(nèi)的輪軌力是非線性變化的,時(shí)間步長(zhǎng)越長(zhǎng),與實(shí)際的輪軌接觸狀態(tài)差別越大。通過(guò)仿真計(jì)算,認(rèn)為采用輪軌時(shí)變接觸剛度時(shí),積分時(shí)間0.05 ms較合理。

    3 模型驗(yàn)證結(jié)果與計(jì)算結(jié)果比較

    3.1 脈沖型短波不平順激勵(lì)

    為驗(yàn)證基于輪軌接觸時(shí)變剛度系數(shù)的車(chē)輛?軌道耦合振動(dòng)分析模型的正確性,采用與文獻(xiàn)[18]中相同的計(jì)算條件,以長(zhǎng)鋼軌焊接區(qū)軌面短波不平順(短波波長(zhǎng)為0.1 m,長(zhǎng)波和短波波幅分別為0.3 mm和0.2 mm)為軌道不平順激勵(lì)波形,見(jiàn)圖3。跟蹤計(jì)算第1輪對(duì)及其對(duì)應(yīng)軌道結(jié)構(gòu)的動(dòng)力響應(yīng)時(shí)程(如輪軌垂向作用力、鋼軌加速度、軌道板加速度和鋼軌位移等),計(jì)算結(jié)果見(jiàn)圖4。

    圖3 鋼軌焊接區(qū)軌面短波不平順激勵(lì)波形

    從圖4及文獻(xiàn)[21]中計(jì)算結(jié)果可知:本文計(jì)算所獲得的輪軌作用力、鋼軌加速度響應(yīng)波幅及鋼軌位移響應(yīng)波形幾乎一致,但軌道板振動(dòng)加速度及鋼軌位移響應(yīng)幅值與文獻(xiàn)[21]中結(jié)果存在較大差別,這與本文有限元模型劃分精度及積分步長(zhǎng)有一定的關(guān)系??偟膩?lái)說(shuō),本文的計(jì)算模型是正確的。

    同時(shí),比較本文模型與輪軌密貼模型的計(jì)算結(jié)果(圖5,積分步長(zhǎng)0.5 ms)可知輪軌密貼模型的計(jì)算結(jié)果與文獻(xiàn)[21]中結(jié)果和本文模型計(jì)算結(jié)果均存在較大差異,其輪軌作用力產(chǎn)生的振動(dòng)響應(yīng)遠(yuǎn)比輪軌彈性接觸模型的的振動(dòng)響應(yīng)大。其原因在于輪軌密貼是一種輪軌剛性接觸的處理方法。同時(shí),由式(3)~(5)所示的輪軌約束方程可知,脈沖型短波不平順的附加速度及加速度(見(jiàn)圖6)對(duì)輪軌之間的相互作用產(chǎn)生很大的影響,而實(shí)際的輪軌彈性接觸狀態(tài)能極大地弱化這種作用。雖然輪軌密貼模型計(jì)算結(jié)果表明存在跳軌現(xiàn)象,但其跳軌持續(xù)時(shí)間及響應(yīng)幅值沒(méi)有得到正確反映。

    3.2 軌道隨機(jī)不平順中波和長(zhǎng)波激勵(lì)

    這里仍然采用與3.1節(jié)中相同的系統(tǒng)動(dòng)力參數(shù)。通過(guò)統(tǒng)計(jì)、歸納武廣高速鐵路2013年近9月的軌道高低不平順譜,采用逆傅里葉變換方法[22],反演獲得軌道高低隨機(jī)不平順,波長(zhǎng)范圍為1~60 m,可作為中波和長(zhǎng)波不平順激勵(lì)源,見(jiàn)圖7。2種模型的積分時(shí)間均取為0.05 ms。

    (a) 輪軌作用力響應(yīng);(b) 鋼軌加速度響應(yīng);(c) 軌道板加速度響應(yīng);(d) 鋼軌位移響應(yīng)

    (a) 輪軌作用力響應(yīng);(b) 鋼軌加速度響應(yīng);(c) 軌道板加速度響應(yīng);(d) 鋼軌位移響應(yīng)

    圖6 脈沖短波不平順的附加加速度

    圖7 軌道高低不平順波形

    3.2.1 車(chē)輛系統(tǒng)振動(dòng)響應(yīng)

    車(chē)體和前構(gòu)架的振動(dòng)位移及振動(dòng)加速度計(jì)算波形及幅值見(jiàn)圖8。由于輪軌時(shí)變接觸剛度系數(shù)模型和輪軌密貼模型對(duì)輪對(duì)位移參數(shù)的考慮完全不同,這里不進(jìn)行不同車(chē)輛振動(dòng)系統(tǒng)之間的比較。

    從圖8可知:2種模型在中、長(zhǎng)波不平順激勵(lì)下產(chǎn)生的車(chē)輛系統(tǒng)振動(dòng)響應(yīng)基本一致。這說(shuō)明在鐵路線路狀態(tài)較緩和、無(wú)短波高頻激勵(lì)時(shí),輪軌的彈性或剛性接觸對(duì)車(chē)輛系統(tǒng)的振動(dòng)響應(yīng)影響不大。

    3.2.2 軌道系統(tǒng)振動(dòng)響應(yīng)

    某里程處鋼軌及軌道板單元中點(diǎn)的振動(dòng)位移及加速度響應(yīng)波形與幅值比較見(jiàn)圖9。

    從圖9可知:2種模型獲得的軌道系統(tǒng)計(jì)算結(jié)果相差不大。但作為輪軌系統(tǒng)的直接接觸單元,輪軌密貼模型的鋼軌位移及加速度響應(yīng)比輪軌彈性接觸模型的稍大。

    3.2.3 輪軌力響應(yīng)

    輪軌力是車(chē)輛?軌道系統(tǒng)振動(dòng)特征最顯著和最直接的表達(dá)參數(shù),這里對(duì)2種模型第1輪對(duì)的輪軌力時(shí)域及頻域特征進(jìn)行比較分析,見(jiàn)圖10。

    從圖10可知:2種模型的頻率特征分布基本一致;在輪軌力的時(shí)域分布上,輪軌時(shí)變剛度模型的輪軌力基本在輪對(duì)軸重附近波動(dòng);但由式(6)~(9)計(jì)算獲得的輪軌密貼模型輪軌力的整體波動(dòng)均值比輪軌時(shí)變剛度系數(shù)模型的大(這與鋼軌振動(dòng)響應(yīng)較大的計(jì)算結(jié)果一致),改變積分步長(zhǎng)及運(yùn)行速度,仍得到類似結(jié)果,并且運(yùn)行速度越快,積分步長(zhǎng)越短,這種差別越明顯。

    (a) 車(chē)體振動(dòng)位移;(b) 前構(gòu)架振動(dòng)位移;(c) 車(chē)體振動(dòng)加速度;(d) 前構(gòu)架振動(dòng)加速度

    (a) 鋼軌中點(diǎn)振動(dòng)位移;(b) 軌道板振動(dòng)位移;(c) 鋼軌中點(diǎn)振動(dòng)加速度;(d) 軌道板振動(dòng)加速度

    (a) 輪軌力; (b) 輪軌力頻譜幅值

    需要指出的是:基于達(dá)朗貝爾原理,采用式(6)~(9)計(jì)算輪軌力是正確的。本文采用輪軌時(shí)變接觸剛度模型,根據(jù)每一積分步獲得的車(chē)輛系統(tǒng)振動(dòng)響應(yīng)用式(6)~(9)計(jì)算相應(yīng)的輪軌力,并與按赫茲接觸理論計(jì)算結(jié)果進(jìn)行比較,見(jiàn)圖11。

    從圖11可知:基于赫茲接觸理論和基于達(dá)朗貝爾原理的輪軌力計(jì)算值基本一致,由此表明式(6)~(9)用于輪軌力計(jì)算的正確性。另一方面,圖11所示計(jì)算結(jié)果表明了力學(xué)原理的統(tǒng)一性,無(wú)論輪軌之間的相互作用多么復(fù)雜,若單獨(dú)隔離車(chē)輛系統(tǒng)或軌道系統(tǒng),則在某一瞬時(shí),它必定處在動(dòng)力平衡狀態(tài)。文獻(xiàn)[2]采用達(dá)朗貝爾原理計(jì)算的輪軸橫向搖擺力與實(shí)測(cè)結(jié)果接近,說(shuō)明車(chē)輛?軌道系統(tǒng)的垂向、橫向振動(dòng)均滿足這一基本原理。

    (a) 第1輪對(duì)輪軌力時(shí)程;(b) 第3輪軌輪軌力時(shí)程

    4 結(jié)論

    1) 對(duì)于輪軌非線性接觸問(wèn)題,采用“局部線性?整體非線性”的方法是可行的,解決了輪軌接觸位移函數(shù)的選取問(wèn)題,便于用能量變分法建立統(tǒng)一的車(chē)輛?軌道耦合振動(dòng)方程。

    2) 取輪軌相互作用力為內(nèi)力,軌道(橋梁)的邊界條件為車(chē)?軌(橋)系統(tǒng)的邊界條件。將輪軌力“隱式”處理為輪軌時(shí)變接觸剛度系數(shù)的動(dòng)力分析方法是可行的。

    3) 輪軌時(shí)變剛度模型與輪軌密貼模型在中、長(zhǎng)波不平順激勵(lì)下,產(chǎn)生的系統(tǒng)動(dòng)力響應(yīng)基本一致,但當(dāng)輪軌密貼模型在積分步長(zhǎng)較小時(shí),其剛性接觸的假設(shè)將逐漸脫離輪軌實(shí)際的彈性接觸狀態(tài),并且動(dòng)力響應(yīng)較大。

    4) 輪軌密貼模型中所采用輪軌剛性接觸方法不適合線路短波不平順激勵(lì)工況。

    [1] 翟婉明. 車(chē)輛-軌道耦合動(dòng)力學(xué)[M]. 北京: 科學(xué)出版社, 2007: 32?55. ZHAI Wanming. Vehicle and track coupling dynamics[M]. Beijing: Science Press, 2007: 32?55.

    [2] 曾慶元. 列車(chē)橋梁時(shí)變系統(tǒng)振動(dòng)分析理論與應(yīng)用[M]. 北京: 中國(guó)鐵道出版社, 1999: 49?56. ZENG Qingyuan. Theory of vibration analysis of train-bridge time- variant system and its applications[M]. Beijing: China Railway Press, 1999: 49?56.

    [3] 婁平, 曾慶元. 車(chē)輛?軌道?橋梁系統(tǒng)豎向運(yùn)動(dòng)方程的建立[J]. 鐵道學(xué)報(bào), 2004, 26(5): 71?80. LOU Ping, ZENG Qingyuan. Formulation of equations of vertical motion for vehicle-track-bridge system[J]. Journal of the China Railway Society, 2004, 26(5): 71?80.

    [4] 張敏, 張楠, 夏禾. 大跨度鐵路懸索橋風(fēng)?車(chē)?橋耦合動(dòng)力分析[J]. 中國(guó)鐵道科學(xué), 2013, 34(4): 14?21. ZHANG Min, ZHANG Nan, XIA He. Analysis on wind-vehicle-bridge dynamic interaction for long-span railway suspension bridge[J]. China Railway Science, 2013, 34(4): 14?21.

    [5] 羅震, 翟婉明, 蔡成標(biāo), 等. 車(chē)輛作用下板式軌道動(dòng)力分析模型及驗(yàn)證[J]. 鐵道學(xué)報(bào), 2010, 32(4): 70?73. LUO Zhen, ZHAI Wanming, CAI Chengbiao, et al. Dynamics model of slab track subjected to moving vehicle and its verification[J]. Journal of the China Railway Society, 2010, 32(4): 70?73.

    [6] ZHAI Wanming. Two simple fast integration methods for large-scale dynamic problems in engineering[J]. International Journal for Numerical Methods in Engineering, 1996, 39(24): 4199?4214.

    [7] NEWMARK N M. A method of computation for structural dynamics[J]. Journal of the Engineering Mechanical Division, ASCE, 1959, 85(2): 67?94.

    [8] 向俊, 赫丹, 曾慶元. 橫向有限條與無(wú)砟軌道板單元的車(chē)軌系統(tǒng)豎向振動(dòng)分析法[J]. 鐵道學(xué)報(bào), 2007, 29(4): 64?69. XIANG Jun, HE Dan, ZENG Qingyuan. Analysis method of vertical vibration of train and ballastless track system with the lateral finite strip and slab segment element[J]. Journal of the China Railway Society, 2007, 29(4): 64?69.

    [9] 雷曉燕, 張斌, 劉慶杰. 軌道過(guò)渡段動(dòng)力特性的有限元分析[J]. 中國(guó)鐵道科學(xué), 2009, 30(5): 15?21. LEI Xiaoyan, ZHANG Bin, LIU Qingjie. Finite element analysis on the dynamic characteristics of the track transition[J]. China Railway Science, 2009, 30(5): 15?21.

    [10] 楊新文, 宮全美, 周順華, 等. 高速列車(chē)作用下雙塊式無(wú)砟軌道與路基垂向耦合振動(dòng)分析[J]. 鐵道學(xué)報(bào), 2014, 36(8): 75?83. YANG Xinwen, GONG Quanmei, ZHOU Shunhua, et al. Analysis on vertical vibration of coupled double-block ballastless track and subgrade system under high-speed train running[J]. Journal of China Railway Society, 2014, 36(8): 75?83.

    [11] 楊宏印, 張海龍, 陳志軍, 等. 考慮徐變影響的車(chē)-線-橋耦合有限元分析[J]. 中國(guó)鐵道科學(xué), 2013, 34(6): 13?19. YANG Hongyin, ZHANG Hailong, CHEN Zhijun, et al. Finite element analysis of train-track-bridge interaction system considering the influence of creep[J]. China Railway Science, 2013, 34(6): 13?19.

    [12] CECíLIA VALE, RUI CAL?ADA. A dynamic vehicle-track interaction model for predicting the track degradation process[J]. Journal of Infrastructure Systems, 2014, 20: 1?13.

    [13] 周智輝, 曾慶元. 列車(chē)-橋梁(軌道)系統(tǒng)振動(dòng)方程解的適定性分析[J]. 力學(xué)與實(shí)踐, 2005, 27(5): 46?49. ZHOU Zhihui, ZENG Qingyuan. Well-posedness of vibration equation of train-bridge (track) system[J]. Mechanics in Engineering, 2005, 27(5): 46?49.

    [14] 錢(qián)偉長(zhǎng). 變分法及有限元[M]. 北京: 科學(xué)出版社, 1980: 26?37. QIAN Weichang. Variational methodand finite element method[M]. Beijing: Science Press, 1980: 26?37.

    [15] 普齊米尼斯基. 矩陣結(jié)構(gòu)分析理論[M]. 王德榮,譯. 北京: 國(guó)防工業(yè)出版社, 1974: 25?36. PRZEMIENIECKI J S. Theory of matrix structural analysis[M]. WANG Derong, trans. Beijing: National Defense Industry Press, 1974: 25?36.

    [16] 曾慶元. 彈性系統(tǒng)動(dòng)力學(xué)總勢(shì)能不變值原理[J]. 華中理工大學(xué)學(xué)報(bào), 2000, 28(1): 1?3. ZENG Qingyuan. The principle of total potential energy with stationary value in elastic system dynamics[J]. Journal of Huangzhong University of Science and Technology, 2000, 28(1): 1?3.

    [17] 李東平. 車(chē)輛多體系統(tǒng)動(dòng)力學(xué)方程的有限元法[J]. 中國(guó)鐵道科學(xué), 2004, 25(5): 33?38. LI Dongping. The finite element method for establishing dynamical equations of vehicle multi-body system[J]. China Railway Science, 2004, 25(5): 33?38.

    [18] 婁平, 曾慶元. 移動(dòng)荷載作用下板式軌道的有限元分析[J]. 交通運(yùn)輸工程學(xué)報(bào), 2004, 4(1): 29?33. LOU Ping, ZENG Qingyuan. Finite element analysis of slab track subjected to moving load[J]. Journal of Traffic and Transportation Engineering, 2004, 4(1): 29?33.

    [19] 趙懷耘, 劉建新, 翟婉明. 板式軌道動(dòng)力響應(yīng)分析方法[J]. 交通運(yùn)輸工程學(xué)報(bào), 2007, 7(5): 19?23. ZHAO Huaiyun, LIU Jianxin, ZHAI Wanming. Analysis method of dynamic responses on slab track[J]. Journal of Traffic and Transportation, 2007, 7(5): 19?23.

    [20] 曾慶元, 楊平. 形成矩陣的“對(duì)號(hào)入座”法則與桁梁空間分析的桁段單元法[J]. 鐵道學(xué)報(bào), 1986, 8(2): 48?59. ZENG Qingyuan, YANG Ping. The “set-in-right-position” rule for forming structural matrices and the finite truss-element method for space analysis of truss bridges[J]. Journal of the China Railway Society, 1986, 8(2): 48?59.

    [21] 翟婉明, 韓衛(wèi)軍, 蔡成標(biāo), 等. 高速鐵路板式軌道動(dòng)力特性研究[J]. 鐵道學(xué)報(bào), 1999, 21(6): 65?69. ZHAI Wanming, HAN Weijun, CAI Chengbiao, et al. Dynamic properties of high-speed railway slab tracks[J]. Journal of the China Railway Science, 1999, 21(6): 65?69.

    [22] 陳果, 翟婉明. 鐵路軌道不平順隨機(jī)過(guò)程的數(shù)值模擬[J]. 西南交通大學(xué)學(xué)報(bào), 1999, 34(2): 138?142. CHEN Guo, ZHAI Wanming. Numerical simulation of the stochastic process of railway track irregularities[J]. Journal of Southwest Jiaotong University, 1999, 34(2): 138?142.

    (編輯 陳燦華)

    Vehicle-track system vertical vibration analysis based on wheel-railtime-varying contact stiffness

    HE Tianlong1, 2, CHEN Xianmai1, 2, XU Lei3

    (1. School of Civil Engineering, Central South University, Changsha 410075, China;2. National Engineering Laboratory for High Speed Railway Construction, Changsha 410075, China;3. Track Power State Key Laboratory, Southwest Jiaotong University, Chengdu 610031, China)

    A vertical and unified equation of the vehicle-track system was established with finite element theory and energy variation method, and wheel-rail time-varying contact stiffness replaced wheel-track nonlinear contact force on dynamic effects to simulate the natural contact between wheel and rail. The validity of the model was verified by inputting the long and medium-wave irregularity of track and pulsed shortwave in this model and the dynamic response was calculated. The tight contact model between wheel and rail, which was popular in China was compared and analyzed.The results show that analyzing dynamic response of the vehicle-track system at the tight contact model for the long and medium-wave irregularity of track is feasible. However, the rigid method between wheel and track will amplify the incentive effect of shortwave for the shortwave irregularity of track, the treatment method about jumping rail should be improved.

    vehicle-track coupling system; time-varying contact stiffness; vertical vibration; energy variation method; unified equation

    10.11817/j.issn.1672?7207.2017.08.023

    U213.2;U211.5

    A

    1672?7207(2017)08?2142?10

    2016?08?23;

    2016?10?22

    國(guó)家自然科學(xué)基金資助項(xiàng)目(51478482);中國(guó)鐵路總公司科技研究開(kāi)發(fā)計(jì)劃項(xiàng)目(Z2013-G006);武漢鐵路局科技發(fā)展計(jì)劃項(xiàng)目(13G04)(Project(51478482) supported by the National Natural Science Foundation of China; Project(Z2013-G006) supported by the Science and Technology Research and Development of China Railway Corporation; Project(13G04) supported by the Science and Technology Development Plan of Wuhan Railway Administration)

    陳憲麥,副教授,從事軌道動(dòng)力學(xué)、線路評(píng)估及養(yǎng)護(hù)維修等研究;E-mail:xianmaichen@aliyun.com

    猜你喜歡
    密貼平順輪軌
    基于Simulink的汽車(chē)行駛平順性研究
    中低速磁浮道岔與輪軌道岔的差異
    淺析密貼檢查器及防潮不良引起的道岔故障
    平順植保站:開(kāi)展粟灰螟防治
    ZYJ7型道岔維修與機(jī)械故障處理
    中低速磁浮與輪軌交通信號(hào)系統(tǒng)的差異
    廣佛線密貼檢查器電路改造方案
    非線性穩(wěn)態(tài)曲線通過(guò)時(shí)輪軌滾動(dòng)接觸的數(shù)值求解方法
    不同輪軌沖角下高速輪軌穩(wěn)態(tài)滾動(dòng)接觸的蠕滑特性
    小號(hào)碼道岔扳動(dòng)力隨密貼段剛度變化規(guī)律研究
    欧美日韩av久久| av在线观看视频网站免费| 国产成人午夜福利电影在线观看| 久久久a久久爽久久v久久| 日本黄色片子视频| av不卡在线播放| 18禁在线播放成人免费| 国产亚洲欧美精品永久| 99久久中文字幕三级久久日本| 免费观看a级毛片全部| 一级毛片黄色毛片免费观看视频| 精品一品国产午夜福利视频| 在线亚洲精品国产二区图片欧美 | h视频一区二区三区| 啦啦啦啦在线视频资源| 一区二区三区免费毛片| 久久影院123| 亚洲欧洲日产国产| 免费人妻精品一区二区三区视频| 国产探花极品一区二区| 免费黄频网站在线观看国产| 麻豆乱淫一区二区| 最近手机中文字幕大全| 久久精品久久精品一区二区三区| 在线观看免费日韩欧美大片 | 亚洲伊人久久精品综合| 两个人的视频大全免费| 亚洲,欧美,日韩| 久久狼人影院| 亚洲人成网站在线播| 少妇熟女欧美另类| 一本久久精品| 欧美精品国产亚洲| 这个男人来自地球电影免费观看 | 久久精品人人爽人人爽视色| 男人添女人高潮全过程视频| 亚洲精品456在线播放app| 免费看av在线观看网站| 欧美日韩亚洲高清精品| 另类精品久久| av电影中文网址| √禁漫天堂资源中文www| 亚洲精品日韩av片在线观看| 国产成人精品无人区| 99视频精品全部免费 在线| 国产成人精品在线电影| 丝袜喷水一区| 亚洲精品中文字幕在线视频| 一区二区三区乱码不卡18| 在线观看美女被高潮喷水网站| 亚洲av电影在线观看一区二区三区| 亚洲欧洲精品一区二区精品久久久 | 一二三四中文在线观看免费高清| 日日撸夜夜添| 国产女主播在线喷水免费视频网站| tube8黄色片| 999精品在线视频| 简卡轻食公司| 久久综合国产亚洲精品| 丰满迷人的少妇在线观看| 国产日韩一区二区三区精品不卡 | 精品国产一区二区久久| 一级,二级,三级黄色视频| 免费观看在线日韩| 国产欧美亚洲国产| 国产精品人妻久久久久久| 久久久久久久亚洲中文字幕| 香蕉精品网在线| 成年人午夜在线观看视频| 国产 精品1| 99九九在线精品视频| 国产精品一区二区在线不卡| 日韩精品有码人妻一区| 男女边摸边吃奶| 人体艺术视频欧美日本| 亚洲性久久影院| 久久婷婷青草| 极品少妇高潮喷水抽搐| 欧美成人精品欧美一级黄| 欧美xxxx性猛交bbbb| 国产在视频线精品| 在线观看三级黄色| 国产黄色免费在线视频| 少妇精品久久久久久久| 如何舔出高潮| 亚洲精品第二区| 水蜜桃什么品种好| 最近手机中文字幕大全| 日日摸夜夜添夜夜爱| 国产精品久久久久久av不卡| 精品卡一卡二卡四卡免费| 国产精品久久久久久av不卡| 丰满乱子伦码专区| freevideosex欧美| 亚洲精品亚洲一区二区| 熟女电影av网| 国产精品无大码| 国产一区二区在线观看av| 一级a做视频免费观看| 少妇的逼好多水| 啦啦啦在线观看免费高清www| 亚洲国产精品一区三区| 精品国产国语对白av| 国产精品99久久99久久久不卡 | 午夜久久久在线观看| 亚州av有码| 亚洲精品,欧美精品| 少妇丰满av| 我的女老师完整版在线观看| 一本—道久久a久久精品蜜桃钙片| 欧美精品一区二区大全| 九九在线视频观看精品| 99精国产麻豆久久婷婷| 久久久久久久久久人人人人人人| 国产高清不卡午夜福利| 欧美三级亚洲精品| av在线观看视频网站免费| 国产爽快片一区二区三区| 一区二区三区免费毛片| 一本久久精品| 日韩免费高清中文字幕av| 97精品久久久久久久久久精品| 少妇人妻精品综合一区二区| 欧美3d第一页| 国产成人aa在线观看| 久久久久精品性色| 久久精品国产亚洲网站| 亚洲国产欧美日韩在线播放| av黄色大香蕉| 亚洲图色成人| 乱码一卡2卡4卡精品| 免费黄频网站在线观看国产| 免费黄频网站在线观看国产| 一级爰片在线观看| tube8黄色片| 美女中出高潮动态图| 人妻 亚洲 视频| 国产毛片在线视频| 欧美激情极品国产一区二区三区 | 亚洲精品中文字幕在线视频| 青春草国产在线视频| 人妻系列 视频| 麻豆精品久久久久久蜜桃| 亚洲精品第二区| 亚洲欧美日韩另类电影网站| 国产精品久久久久久精品电影小说| 免费av中文字幕在线| 国产欧美另类精品又又久久亚洲欧美| 免费久久久久久久精品成人欧美视频 | 欧美 亚洲 国产 日韩一| 亚洲国产最新在线播放| 精品国产国语对白av| 99热6这里只有精品| 91成人精品电影| 国产探花极品一区二区| 人人妻人人澡人人爽人人夜夜| 大陆偷拍与自拍| 欧美三级亚洲精品| 欧美日韩在线观看h| 国产成人91sexporn| 九色成人免费人妻av| 国产精品人妻久久久影院| 又粗又硬又长又爽又黄的视频| av天堂久久9| 99久久综合免费| 免费黄色在线免费观看| 亚洲国产精品国产精品| 久久久a久久爽久久v久久| 乱码一卡2卡4卡精品| 免费少妇av软件| 在线观看免费视频网站a站| 久久久欧美国产精品| 国产成人精品在线电影| 18在线观看网站| 最近最新中文字幕免费大全7| 亚洲欧洲日产国产| 校园人妻丝袜中文字幕| 中文乱码字字幕精品一区二区三区| 少妇 在线观看| 三上悠亚av全集在线观看| 99久久精品一区二区三区| 大陆偷拍与自拍| av免费观看日本| 日本猛色少妇xxxxx猛交久久| 国产日韩一区二区三区精品不卡 | 午夜福利影视在线免费观看| 欧美3d第一页| 欧美变态另类bdsm刘玥| 寂寞人妻少妇视频99o| 欧美日韩在线观看h| 中文精品一卡2卡3卡4更新| 免费观看性生交大片5| 国产免费视频播放在线视频| 丰满少妇做爰视频| 男人爽女人下面视频在线观看| 在线免费观看不下载黄p国产| 如何舔出高潮| 中文字幕人妻熟人妻熟丝袜美| 欧美97在线视频| 久久久久久久久久成人| 亚洲第一区二区三区不卡| 成人黄色视频免费在线看| 国产精品秋霞免费鲁丝片| 99久久中文字幕三级久久日本| 亚洲综合色网址| 亚洲国产精品一区三区| 免费观看无遮挡的男女| 一本久久精品| 午夜激情av网站| 久久99热6这里只有精品| 久久久久久久大尺度免费视频| 涩涩av久久男人的天堂| 久久精品熟女亚洲av麻豆精品| 热99久久久久精品小说推荐| 亚洲国产日韩一区二区| 欧美日韩综合久久久久久| av在线观看视频网站免费| 亚洲国产av新网站| 一级二级三级毛片免费看| 成人无遮挡网站| 麻豆乱淫一区二区| 日韩成人av中文字幕在线观看| av国产久精品久网站免费入址| 大香蕉97超碰在线| 亚洲不卡免费看| 色5月婷婷丁香| 成人毛片60女人毛片免费| 日本午夜av视频| 免费av中文字幕在线| 大片免费播放器 马上看| 亚洲精品第二区| 一级毛片aaaaaa免费看小| 国产欧美亚洲国产| 日本黄大片高清| 亚洲精品乱码久久久久久按摩| 视频区图区小说| 久久影院123| 亚洲精品成人av观看孕妇| 一级爰片在线观看| www.av在线官网国产| 亚洲av二区三区四区| 亚洲人成77777在线视频| 久久99一区二区三区| 国产精品成人在线| 97超视频在线观看视频| 91精品国产国语对白视频| 国产亚洲一区二区精品| 国产精品人妻久久久久久| 我的老师免费观看完整版| 美女国产视频在线观看| 亚洲情色 制服丝袜| 国产精品99久久99久久久不卡 | 夫妻午夜视频| 丝袜美足系列| 亚洲精品乱码久久久v下载方式| kizo精华| 人成视频在线观看免费观看| 精品久久久久久电影网| 国产亚洲欧美精品永久| videossex国产| 久热久热在线精品观看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | av在线播放精品| 丝袜脚勾引网站| 亚洲av在线观看美女高潮| 国产毛片在线视频| 成人毛片a级毛片在线播放| 精品国产一区二区久久| 久久99热这里只频精品6学生| 亚洲人与动物交配视频| 伊人久久国产一区二区| 国产成人午夜福利电影在线观看| 国产男女超爽视频在线观看| 久久99精品国语久久久| 日韩三级伦理在线观看| 一区二区三区四区激情视频| 简卡轻食公司| av线在线观看网站| 亚洲无线观看免费| 午夜福利在线观看免费完整高清在| 在线观看免费日韩欧美大片 | av线在线观看网站| 午夜免费男女啪啪视频观看| 纵有疾风起免费观看全集完整版| 啦啦啦在线观看免费高清www| 99久久精品国产国产毛片| 乱码一卡2卡4卡精品| 日韩在线高清观看一区二区三区| 色网站视频免费| 少妇被粗大猛烈的视频| 欧美少妇被猛烈插入视频| 美女国产高潮福利片在线看| 青春草视频在线免费观看| 亚洲av日韩在线播放| 男女边吃奶边做爰视频| a 毛片基地| 麻豆成人av视频| 人妻夜夜爽99麻豆av| 各种免费的搞黄视频| 春色校园在线视频观看| 午夜影院在线不卡| 亚洲av日韩在线播放| 免费久久久久久久精品成人欧美视频 | 亚洲av综合色区一区| 日韩视频在线欧美| 精品熟女少妇av免费看| 亚洲av男天堂| 国产精品久久久久成人av| 丝袜喷水一区| 欧美三级亚洲精品| 人人妻人人澡人人看| 搡老乐熟女国产| 欧美xxxx性猛交bbbb| 日韩精品免费视频一区二区三区 | 99热这里只有精品一区| 久久女婷五月综合色啪小说| 亚洲国产精品专区欧美| 一本大道久久a久久精品| 亚洲少妇的诱惑av| 国产成人精品在线电影| 国产欧美日韩综合在线一区二区| 好男人视频免费观看在线| 性色av一级| 国产伦理片在线播放av一区| 日本与韩国留学比较| 综合色丁香网| 啦啦啦啦在线视频资源| 岛国毛片在线播放| 美女主播在线视频| 国产综合精华液| 秋霞在线观看毛片| 国产av一区二区精品久久| 亚洲av国产av综合av卡| 久久久久人妻精品一区果冻| 亚洲成人一二三区av| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 婷婷色麻豆天堂久久| 一级a做视频免费观看| a级毛片黄视频| 国产淫语在线视频| 久久热精品热| 久久久久人妻精品一区果冻| 亚洲国产欧美日韩在线播放| 国产高清不卡午夜福利| 精品久久久久久电影网| 成人无遮挡网站| 高清毛片免费看| 少妇高潮的动态图| 精品久久久久久久久亚洲| 欧美日韩视频高清一区二区三区二| 一级黄片播放器| 午夜视频国产福利| 成人毛片a级毛片在线播放| 自拍欧美九色日韩亚洲蝌蚪91| 春色校园在线视频观看| 国产亚洲精品第一综合不卡 | 午夜福利,免费看| 国产精品偷伦视频观看了| a级毛片免费高清观看在线播放| 五月天丁香电影| 国产色婷婷99| 人妻夜夜爽99麻豆av| 中文精品一卡2卡3卡4更新| 免费日韩欧美在线观看| 91成人精品电影| 一本大道久久a久久精品| 大又大粗又爽又黄少妇毛片口| 亚洲激情五月婷婷啪啪| 中文字幕最新亚洲高清| 国产一级毛片在线| 人妻制服诱惑在线中文字幕| 国产 一区精品| 最近中文字幕高清免费大全6| a级毛色黄片| 99re6热这里在线精品视频| 菩萨蛮人人尽说江南好唐韦庄| 下体分泌物呈黄色| 日日摸夜夜添夜夜添av毛片| 日本午夜av视频| 国产免费视频播放在线视频| av线在线观看网站| 亚洲欧美一区二区三区黑人 | 国产免费又黄又爽又色| 91久久精品电影网| 91久久精品国产一区二区三区| a 毛片基地| 亚洲三级黄色毛片| 制服丝袜香蕉在线| 国产日韩一区二区三区精品不卡 | 欧美激情国产日韩精品一区| 久久久久精品久久久久真实原创| 一区二区三区免费毛片| 伊人久久国产一区二区| 亚洲国产精品一区二区三区在线| 一级毛片电影观看| 嘟嘟电影网在线观看| 大香蕉久久网| 亚洲国产欧美在线一区| 欧美3d第一页| 插逼视频在线观看| 久久久精品94久久精品| 三上悠亚av全集在线观看| 老司机影院毛片| 午夜久久久在线观看| 久热久热在线精品观看| 亚洲成人av在线免费| 日韩三级伦理在线观看| 国产黄色视频一区二区在线观看| 国产成人freesex在线| 18禁裸乳无遮挡动漫免费视频| 久久精品国产亚洲网站| 最近中文字幕2019免费版| 亚洲精品国产色婷婷电影| 久久久久久久大尺度免费视频| 夜夜看夜夜爽夜夜摸| 啦啦啦在线观看免费高清www| 在线 av 中文字幕| 欧美人与善性xxx| 免费不卡的大黄色大毛片视频在线观看| 国产 精品1| 日韩熟女老妇一区二区性免费视频| 制服人妻中文乱码| 日本欧美视频一区| 婷婷色麻豆天堂久久| 能在线免费看毛片的网站| 国产精品久久久久久久久免| 国产男女超爽视频在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 国产精品国产三级国产专区5o| 久久久久久久久久久免费av| 啦啦啦在线观看免费高清www| 热re99久久精品国产66热6| 国产男女超爽视频在线观看| 精品久久久久久久久av| 韩国高清视频一区二区三区| 精品国产一区二区三区久久久樱花| 九色成人免费人妻av| 亚洲精品456在线播放app| 国产女主播在线喷水免费视频网站| 国产高清有码在线观看视频| av黄色大香蕉| 久久免费观看电影| 久久97久久精品| 一级爰片在线观看| 亚洲精品,欧美精品| 亚洲婷婷狠狠爱综合网| 成人综合一区亚洲| 黄片无遮挡物在线观看| 精品少妇内射三级| 精品人妻在线不人妻| 日本色播在线视频| 国内精品宾馆在线| 亚洲久久久国产精品| 亚州av有码| 国产精品久久久久久精品古装| 自拍欧美九色日韩亚洲蝌蚪91| 最新的欧美精品一区二区| 99九九在线精品视频| 夜夜爽夜夜爽视频| 成人黄色视频免费在线看| 日本欧美国产在线视频| 亚洲性久久影院| 女人精品久久久久毛片| 天堂俺去俺来也www色官网| 一级毛片黄色毛片免费观看视频| 日本爱情动作片www.在线观看| 99久久精品一区二区三区| 国产伦精品一区二区三区视频9| 三级国产精品片| 肉色欧美久久久久久久蜜桃| 在线观看美女被高潮喷水网站| 一个人免费看片子| 亚洲熟女精品中文字幕| 亚洲国产色片| 岛国毛片在线播放| 18禁在线播放成人免费| 日本欧美视频一区| 亚洲欧美成人综合另类久久久| 建设人人有责人人尽责人人享有的| 亚洲欧美一区二区三区黑人 | 大又大粗又爽又黄少妇毛片口| 精品久久久久久久久亚洲| 九色亚洲精品在线播放| 中文字幕人妻丝袜制服| 久久韩国三级中文字幕| 久久久精品免费免费高清| 日韩不卡一区二区三区视频在线| 亚洲人与动物交配视频| 久久人妻熟女aⅴ| 一区二区三区四区激情视频| 熟女av电影| 一级a做视频免费观看| 极品少妇高潮喷水抽搐| 国产精品成人在线| 曰老女人黄片| 亚洲经典国产精华液单| 国产片内射在线| 欧美一级a爱片免费观看看| 人人妻人人澡人人看| 国产无遮挡羞羞视频在线观看| 日韩 亚洲 欧美在线| 精品熟女少妇av免费看| 久久 成人 亚洲| 久久毛片免费看一区二区三区| 国产精品国产三级国产av玫瑰| 亚洲av不卡在线观看| 免费高清在线观看视频在线观看| 99久国产av精品国产电影| 亚洲中文av在线| 久久综合国产亚洲精品| 亚洲精品久久成人aⅴ小说 | 久久99热这里只频精品6学生| 亚洲人成网站在线观看播放| 日韩中文字幕视频在线看片| 国产爽快片一区二区三区| 欧美日韩在线观看h| 最近的中文字幕免费完整| 制服丝袜香蕉在线| 午夜福利影视在线免费观看| 国产在线视频一区二区| a级毛片免费高清观看在线播放| 极品少妇高潮喷水抽搐| 69精品国产乱码久久久| 午夜激情福利司机影院| 精品熟女少妇av免费看| 黄色毛片三级朝国网站| 欧美老熟妇乱子伦牲交| 又粗又硬又长又爽又黄的视频| 高清毛片免费看| 国产精品一国产av| 少妇 在线观看| 国产免费福利视频在线观看| 好男人视频免费观看在线| 91精品三级在线观看| 一级毛片电影观看| 国产日韩欧美在线精品| 亚洲怡红院男人天堂| 国产精品久久久久久久久免| 视频在线观看一区二区三区| 婷婷色综合www| 久久女婷五月综合色啪小说| 最近中文字幕高清免费大全6| 18+在线观看网站| 大片电影免费在线观看免费| 高清在线视频一区二区三区| 美女主播在线视频| 在线观看www视频免费| 亚洲精品视频女| 亚洲国产精品成人久久小说| 夜夜看夜夜爽夜夜摸| 王馨瑶露胸无遮挡在线观看| 免费大片黄手机在线观看| 亚洲国产精品一区三区| 亚洲人成网站在线观看播放| 免费久久久久久久精品成人欧美视频 | 美女大奶头黄色视频| 在现免费观看毛片| 成年美女黄网站色视频大全免费 | 性高湖久久久久久久久免费观看| 久久综合国产亚洲精品| 午夜免费男女啪啪视频观看| 午夜久久久在线观看| 免费人成在线观看视频色| 精品久久国产蜜桃| 欧美变态另类bdsm刘玥| 男女国产视频网站| 国产午夜精品一二区理论片| 亚洲av电影在线观看一区二区三区| 亚洲性久久影院| 亚洲精品久久午夜乱码| 黄色一级大片看看| 男男h啪啪无遮挡| 国产一区二区三区综合在线观看 | 国产视频首页在线观看| 两个人免费观看高清视频| 亚洲欧洲精品一区二区精品久久久 | 综合色丁香网| 久久久久久久久久成人| 久久人人爽av亚洲精品天堂| 国产伦理片在线播放av一区| 成人18禁高潮啪啪吃奶动态图 | 久久人妻熟女aⅴ| 国产亚洲精品久久久com| 中文字幕人妻丝袜制服| 2021少妇久久久久久久久久久| 国产免费福利视频在线观看| 亚洲国产色片| 亚洲激情五月婷婷啪啪| 国产日韩欧美视频二区| 超色免费av| 超碰97精品在线观看| 美女主播在线视频| 内地一区二区视频在线| 极品人妻少妇av视频| 丝袜喷水一区| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲三级黄色毛片| 日韩强制内射视频| 亚洲国产精品999| av国产精品久久久久影院| 日韩熟女老妇一区二区性免费视频| av国产久精品久网站免费入址| 肉色欧美久久久久久久蜜桃| 美女国产高潮福利片在线看| 91午夜精品亚洲一区二区三区| 久久99热6这里只有精品| 哪个播放器可以免费观看大片| 午夜精品国产一区二区电影| 亚洲欧美成人综合另类久久久| 黑人巨大精品欧美一区二区蜜桃 | 国产熟女午夜一区二区三区 | 汤姆久久久久久久影院中文字幕| 国产视频内射| 国产一区二区在线观看日韩| 亚洲国产精品国产精品|