• <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ī)律研究
    99国产精品一区二区蜜桃av | 日日摸夜夜添夜夜添小说| 亚洲午夜精品一区,二区,三区| 美女视频免费永久观看网站| 波多野结衣av一区二区av| 亚洲精品一二三| 久久人妻熟女aⅴ| 中文欧美无线码| 国产黄频视频在线观看| 黄色视频不卡| 免费在线观看视频国产中文字幕亚洲| 多毛熟女@视频| av一本久久久久| 好男人电影高清在线观看| 在线观看舔阴道视频| 成人18禁在线播放| 美女主播在线视频| 满18在线观看网站| 欧美日韩中文字幕国产精品一区二区三区 | 国产欧美日韩精品亚洲av| 久久久欧美国产精品| 另类亚洲欧美激情| 男女午夜视频在线观看| 欧美+亚洲+日韩+国产| 亚洲国产成人一精品久久久| 亚洲第一青青草原| cao死你这个sao货| 少妇的丰满在线观看| 国产成人免费无遮挡视频| 国产视频一区二区在线看| 欧美在线黄色| 男女下面插进去视频免费观看| 性少妇av在线| 亚洲 欧美一区二区三区| 亚洲精品成人av观看孕妇| 美女午夜性视频免费| 日韩视频一区二区在线观看| 久久久精品国产亚洲av高清涩受| 国产成人精品久久二区二区免费| av线在线观看网站| 亚洲精品粉嫩美女一区| 日本一区二区免费在线视频| 精品少妇久久久久久888优播| 97人妻天天添夜夜摸| 成年版毛片免费区| 国产野战对白在线观看| 香蕉丝袜av| 青草久久国产| 肉色欧美久久久久久久蜜桃| 欧美一级毛片孕妇| 色婷婷av一区二区三区视频| 日本一区二区免费在线视频| 99久久国产精品久久久| 操出白浆在线播放| 免费人妻精品一区二区三区视频| 国产亚洲精品久久久久5区| videosex国产| 欧美日韩中文字幕国产精品一区二区三区 | 高清av免费在线| 久久久国产一区二区| 国产精品久久电影中文字幕 | 久久久久久久久久久久大奶| 久久久久国内视频| 久久久国产一区二区| 亚洲专区字幕在线| 日韩视频在线欧美| 精品福利永久在线观看| 高清在线国产一区| 国产精品亚洲一级av第二区| 蜜桃在线观看..| 亚洲av欧美aⅴ国产| 精品午夜福利视频在线观看一区 | 午夜福利,免费看| 色精品久久人妻99蜜桃| 久久人妻av系列| 一级a爱视频在线免费观看| 亚洲国产看品久久| 亚洲国产欧美一区二区综合| 法律面前人人平等表现在哪些方面| 免费日韩欧美在线观看| 日韩成人在线观看一区二区三区| 亚洲一码二码三码区别大吗| 久久久久视频综合| 激情在线观看视频在线高清 | 丝袜美腿诱惑在线| 亚洲国产av影院在线观看| 日日摸夜夜添夜夜添小说| 欧美日韩亚洲高清精品| 大香蕉久久网| 操美女的视频在线观看| 欧美 亚洲 国产 日韩一| 亚洲九九香蕉| 亚洲国产欧美在线一区| av免费在线观看网站| 国产高清激情床上av| 天天操日日干夜夜撸| 精品久久蜜臀av无| 欧美性长视频在线观看| 日韩欧美一区二区三区在线观看 | 男女边摸边吃奶| 国产亚洲午夜精品一区二区久久| 亚洲欧美激情在线| 亚洲视频免费观看视频| 免费看十八禁软件| 一级毛片女人18水好多| 日韩欧美国产一区二区入口| 国产伦人伦偷精品视频| 国产精品欧美亚洲77777| 我要看黄色一级片免费的| 一个人免费看片子| 亚洲七黄色美女视频| 色综合欧美亚洲国产小说| 欧美性长视频在线观看| 国产91精品成人一区二区三区 | 久久久久久久大尺度免费视频| 99re6热这里在线精品视频| 如日韩欧美国产精品一区二区三区| av天堂久久9| 一区二区三区激情视频| 欧美黑人欧美精品刺激| 美女主播在线视频| 老司机深夜福利视频在线观看| 国产精品久久久av美女十八| 久久精品成人免费网站| 国产99久久九九免费精品| 午夜视频精品福利| 国产精品自产拍在线观看55亚洲 | 91九色精品人成在线观看| 亚洲精品在线观看二区| 交换朋友夫妻互换小说| 国产一区二区 视频在线| 妹子高潮喷水视频| 久久国产精品大桥未久av| 热99re8久久精品国产| 在线观看舔阴道视频| 成人特级黄色片久久久久久久 | 黄色片一级片一级黄色片| 中文字幕高清在线视频| 午夜福利乱码中文字幕| 麻豆国产av国片精品| 久久国产亚洲av麻豆专区| 精品亚洲乱码少妇综合久久| 亚洲国产欧美一区二区综合| 黄色毛片三级朝国网站| 国产真人三级小视频在线观看| 国产成人精品无人区| 美女扒开内裤让男人捅视频| 欧美 日韩 精品 国产| 午夜福利免费观看在线| 老司机午夜福利在线观看视频 | 99精国产麻豆久久婷婷| 久久ye,这里只有精品| 黄网站色视频无遮挡免费观看| 捣出白浆h1v1| 中文字幕色久视频| 亚洲精品成人av观看孕妇| 亚洲男人天堂网一区| 最黄视频免费看| 黑人操中国人逼视频| 日韩成人在线观看一区二区三区| 大香蕉久久成人网| 午夜福利影视在线免费观看| 久久av网站| 日韩有码中文字幕| 天天躁日日躁夜夜躁夜夜| 国产成人系列免费观看| 后天国语完整版免费观看| 国产1区2区3区精品| 午夜福利视频在线观看免费| 一边摸一边抽搐一进一小说 | 12—13女人毛片做爰片一| 露出奶头的视频| 国产成人欧美在线观看 | 久久精品熟女亚洲av麻豆精品| 中亚洲国语对白在线视频| 天堂8中文在线网| 韩国精品一区二区三区| 亚洲久久久国产精品| a级毛片在线看网站| 日本vs欧美在线观看视频| 国产aⅴ精品一区二区三区波| 日韩免费高清中文字幕av| 成人国产一区最新在线观看| 91精品三级在线观看| 母亲3免费完整高清在线观看| 十八禁人妻一区二区| 一本色道久久久久久精品综合| 91大片在线观看| 老司机深夜福利视频在线观看| 新久久久久国产一级毛片| 国产精品一区二区在线观看99| 极品人妻少妇av视频| 午夜福利视频精品| 18禁观看日本| 成人亚洲精品一区在线观看| 99香蕉大伊视频| 真人做人爱边吃奶动态| 美女视频免费永久观看网站| 宅男免费午夜| 一边摸一边抽搐一进一小说 | 日本撒尿小便嘘嘘汇集6| 亚洲欧洲精品一区二区精品久久久| 亚洲成人手机| 国产一区二区三区视频了| 久久精品国产亚洲av香蕉五月 | 丝袜美腿诱惑在线| 国产欧美日韩综合在线一区二区| 亚洲成国产人片在线观看| 丁香六月欧美| 国产精品亚洲一级av第二区| 搡老岳熟女国产| 久久久国产精品麻豆| 久久性视频一级片| 午夜激情久久久久久久| 精品国产超薄肉色丝袜足j| 亚洲av电影在线进入| 大片电影免费在线观看免费| av天堂在线播放| 老鸭窝网址在线观看| 脱女人内裤的视频| 亚洲国产看品久久| 精品国产一区二区久久| 欧美成人午夜精品| 久久人妻熟女aⅴ| 国产在线精品亚洲第一网站| 免费少妇av软件| 国产精品成人在线| 免费日韩欧美在线观看| 人人妻人人澡人人爽人人夜夜| 亚洲第一欧美日韩一区二区三区 | 久久人妻福利社区极品人妻图片| 久久人人爽av亚洲精品天堂| 国产不卡一卡二| 亚洲熟女精品中文字幕| 欧美人与性动交α欧美软件| 精品一品国产午夜福利视频| 亚洲精华国产精华精| 久久99热这里只频精品6学生| 国产成人系列免费观看| 日本一区二区免费在线视频| 亚洲熟女毛片儿| 极品教师在线免费播放| 欧美久久黑人一区二区| 国产日韩欧美在线精品| 亚洲熟女精品中文字幕| 一级a爱视频在线免费观看| 久久精品人人爽人人爽视色| 搡老乐熟女国产| 精品亚洲乱码少妇综合久久| 久9热在线精品视频| 亚洲va日本ⅴa欧美va伊人久久| 亚洲成av片中文字幕在线观看| 在线av久久热| 午夜精品国产一区二区电影| 免费人妻精品一区二区三区视频| 激情在线观看视频在线高清 | 99精品久久久久人妻精品| 女性被躁到高潮视频| 老汉色av国产亚洲站长工具| 久久久精品94久久精品| 国产精品麻豆人妻色哟哟久久| 欧美成狂野欧美在线观看| 国产精品熟女久久久久浪| 亚洲精品美女久久久久99蜜臀| 婷婷丁香在线五月| 午夜免费成人在线视频| 久久久国产成人免费| 一边摸一边抽搐一进一出视频| 国产欧美日韩精品亚洲av| 国产亚洲欧美精品永久| 女人高潮潮喷娇喘18禁视频| 在线观看免费午夜福利视频| 一区二区三区激情视频| 久久久久久久久久久久大奶| 欧美 亚洲 国产 日韩一| 国产在线视频一区二区| 久久人妻av系列| 亚洲第一av免费看| 色婷婷av一区二区三区视频| 18禁观看日本| 一区在线观看完整版| 亚洲一区中文字幕在线| 肉色欧美久久久久久久蜜桃| 色综合婷婷激情| 亚洲,欧美精品.| 少妇的丰满在线观看| 中文字幕制服av| 亚洲精华国产精华精| 91九色精品人成在线观看| 亚洲成av片中文字幕在线观看| 国产亚洲欧美精品永久| 国产一区有黄有色的免费视频| 69精品国产乱码久久久| 亚洲成人免费电影在线观看| 久久免费观看电影| 久久久久网色| 精品免费久久久久久久清纯 | 欧美亚洲日本最大视频资源| 这个男人来自地球电影免费观看| 老司机在亚洲福利影院| 亚洲av日韩在线播放| 少妇猛男粗大的猛烈进出视频| 91麻豆精品激情在线观看国产 | 久久香蕉激情| 热re99久久精品国产66热6| 天天操日日干夜夜撸| 久久精品国产综合久久久| 亚洲,欧美精品.| 亚洲精品国产一区二区精华液| 乱人伦中国视频| 男人操女人黄网站| 热re99久久精品国产66热6| cao死你这个sao货| 五月天丁香电影| 久久国产精品大桥未久av| 一级毛片女人18水好多| 男女床上黄色一级片免费看| 久久国产精品人妻蜜桃| 久久天堂一区二区三区四区| √禁漫天堂资源中文www| 国产精品麻豆人妻色哟哟久久| 18禁黄网站禁片午夜丰满| 在线观看免费视频日本深夜| 一级片免费观看大全| 香蕉久久夜色| 露出奶头的视频| 日本av手机在线免费观看| 另类亚洲欧美激情| 丰满人妻熟妇乱又伦精品不卡| 久久久国产欧美日韩av| 亚洲国产av影院在线观看| 老司机靠b影院| 午夜福利在线观看吧| av不卡在线播放| 成人亚洲精品一区在线观看| 精品高清国产在线一区| 亚洲av日韩精品久久久久久密| 国产精品一区二区免费欧美| 久久av网站| 一级毛片电影观看| 午夜福利乱码中文字幕| 国产精品免费大片| 国产主播在线观看一区二区| 欧美中文综合在线视频| 国产主播在线观看一区二区| 国产成人啪精品午夜网站| 国产在线一区二区三区精| a级毛片在线看网站| 亚洲精品国产一区二区精华液| 别揉我奶头~嗯~啊~动态视频| 亚洲 欧美一区二区三区| 亚洲伊人色综图| 亚洲国产精品一区二区三区在线| 欧美成人免费av一区二区三区 | 日本vs欧美在线观看视频| 国产在视频线精品| 一区二区日韩欧美中文字幕| 久久午夜亚洲精品久久| 久久精品国产综合久久久| avwww免费| 久久精品亚洲av国产电影网| 日日摸夜夜添夜夜添小说| 在线观看一区二区三区激情| 又大又爽又粗| netflix在线观看网站| 久久亚洲精品不卡| tocl精华| 色播在线永久视频| 午夜福利乱码中文字幕| 亚洲五月婷婷丁香| 久久精品亚洲精品国产色婷小说| 精品国产乱码久久久久久男人| 亚洲中文日韩欧美视频| videosex国产| 99九九在线精品视频| 女性生殖器流出的白浆| 美女主播在线视频| 大型黄色视频在线免费观看| 国产一区二区在线观看av| 91字幕亚洲| 精品国产乱码久久久久久小说| 悠悠久久av| 欧美成狂野欧美在线观看| 岛国在线观看网站| 欧美在线一区亚洲| 国产精品98久久久久久宅男小说| 久久精品成人免费网站| 激情视频va一区二区三区| 又大又爽又粗| 久久久久久免费高清国产稀缺| 免费人妻精品一区二区三区视频| 国产精品久久久久久精品古装| 国产黄频视频在线观看| 国产一区二区三区视频了| 欧美精品一区二区免费开放| 欧美 日韩 精品 国产| 久久久国产成人免费| 悠悠久久av| 别揉我奶头~嗯~啊~动态视频| 12—13女人毛片做爰片一| 亚洲一卡2卡3卡4卡5卡精品中文| 午夜精品国产一区二区电影| 亚洲人成伊人成综合网2020| 一个人免费看片子| 美女扒开内裤让男人捅视频| 他把我摸到了高潮在线观看 | 国产在线免费精品| 亚洲午夜理论影院| 中文字幕人妻丝袜制服| 一边摸一边抽搐一进一出视频| 狠狠狠狠99中文字幕| 在线观看舔阴道视频| 叶爱在线成人免费视频播放| 国产精品熟女久久久久浪| 狠狠狠狠99中文字幕| 一边摸一边抽搐一进一出视频| 交换朋友夫妻互换小说| 日韩一区二区三区影片| 成年人黄色毛片网站| 老汉色∧v一级毛片| 久9热在线精品视频| 国产欧美日韩一区二区三区在线| 性少妇av在线| 免费少妇av软件| 人妻久久中文字幕网| 亚洲精品久久成人aⅴ小说| 欧美中文综合在线视频| 久久中文看片网| 精品第一国产精品| 美女高潮到喷水免费观看| 亚洲成人手机| 精品免费久久久久久久清纯 | 国产日韩欧美视频二区| 交换朋友夫妻互换小说| 中文亚洲av片在线观看爽 | 久久性视频一级片| 狠狠精品人妻久久久久久综合| 亚洲人成电影免费在线| 亚洲伊人色综图| 三级毛片av免费| 亚洲人成电影免费在线| 在线天堂中文资源库| 99久久精品国产亚洲精品| 久久久久久久国产电影| 一区二区三区精品91| 欧美精品亚洲一区二区| 国产精品av久久久久免费| 免费一级毛片在线播放高清视频 | 国产免费av片在线观看野外av| 国产成人一区二区三区免费视频网站| 国产日韩欧美视频二区| 亚洲欧美一区二区三区黑人| 黄色毛片三级朝国网站| a级片在线免费高清观看视频| 亚洲成国产人片在线观看| 天天躁夜夜躁狠狠躁躁| 高清毛片免费观看视频网站 | 天堂动漫精品| av一本久久久久| 天天影视国产精品| 亚洲国产av影院在线观看| 免费在线观看影片大全网站| 最新的欧美精品一区二区| 中文字幕精品免费在线观看视频| 男女免费视频国产| 亚洲三区欧美一区| 99久久国产精品久久久| 亚洲第一av免费看| 最近最新中文字幕大全电影3 | 人人妻人人爽人人添夜夜欢视频| 欧美国产精品一级二级三级| 搡老岳熟女国产| 久9热在线精品视频| 可以免费在线观看a视频的电影网站| 建设人人有责人人尽责人人享有的| 久久久久国内视频| 国产精品1区2区在线观看. | 久久久欧美国产精品| 国产男女内射视频| 国产又色又爽无遮挡免费看| av又黄又爽大尺度在线免费看| 亚洲欧美色中文字幕在线| 日本vs欧美在线观看视频| 亚洲伊人色综图| 日本撒尿小便嘘嘘汇集6| 国产精品影院久久| 日本撒尿小便嘘嘘汇集6| 涩涩av久久男人的天堂| 精品一区二区三区四区五区乱码| 女人被躁到高潮嗷嗷叫费观| 亚洲熟女精品中文字幕| 黄网站色视频无遮挡免费观看| 免费人妻精品一区二区三区视频| 国产亚洲一区二区精品| 国产免费福利视频在线观看| 99九九在线精品视频| 久久久水蜜桃国产精品网| 亚洲欧洲精品一区二区精品久久久| 99re在线观看精品视频| 首页视频小说图片口味搜索| 男人舔女人的私密视频| 成人黄色视频免费在线看| 精品国产亚洲在线| 男女边摸边吃奶| 曰老女人黄片| 久久人妻熟女aⅴ| 国产在视频线精品| 成在线人永久免费视频| 一级毛片精品| 国产欧美亚洲国产| av国产精品久久久久影院| 男女床上黄色一级片免费看| av天堂久久9| 热99久久久久精品小说推荐| 多毛熟女@视频| 天天影视国产精品| 亚洲第一欧美日韩一区二区三区 | 午夜福利在线免费观看网站| 热99国产精品久久久久久7| 免费人妻精品一区二区三区视频| 亚洲精品在线观看二区| www.熟女人妻精品国产| 女人久久www免费人成看片| 色婷婷久久久亚洲欧美| 成人黄色视频免费在线看| 免费黄频网站在线观看国产| 最新美女视频免费是黄的| 日韩熟女老妇一区二区性免费视频| 在线 av 中文字幕| 久9热在线精品视频| 国产精品久久久av美女十八| 久久人人爽av亚洲精品天堂| 日韩欧美免费精品| 成人特级黄色片久久久久久久 | 不卡一级毛片| 亚洲精品乱久久久久久| 亚洲男人天堂网一区| 欧美激情高清一区二区三区| 一边摸一边做爽爽视频免费| 性少妇av在线| 1024香蕉在线观看| 欧美日韩国产mv在线观看视频| 1024视频免费在线观看| 好男人电影高清在线观看| 免费久久久久久久精品成人欧美视频| 19禁男女啪啪无遮挡网站| 五月天丁香电影| 国产熟女午夜一区二区三区| 美女午夜性视频免费| 青青草视频在线视频观看| 国产在线视频一区二区| 在线观看免费高清a一片| 国产精品香港三级国产av潘金莲| 成人亚洲精品一区在线观看| 窝窝影院91人妻| 精品人妻熟女毛片av久久网站| 久久人妻av系列| 啦啦啦 在线观看视频| 精品国产国语对白av| 免费观看av网站的网址| 欧美日韩精品网址| a在线观看视频网站| 老汉色∧v一级毛片| 一区二区av电影网| 色在线成人网| 日韩中文字幕视频在线看片| 另类亚洲欧美激情| 亚洲 欧美一区二区三区| 国产精品av久久久久免费| 一边摸一边抽搐一进一出视频| 亚洲一区二区三区欧美精品| 一区二区三区激情视频| 69精品国产乱码久久久| 午夜免费成人在线视频| 国产亚洲欧美精品永久| 国产一卡二卡三卡精品| 国产99久久九九免费精品| av有码第一页| 女人爽到高潮嗷嗷叫在线视频| 亚洲欧美一区二区三区黑人| 亚洲一区中文字幕在线| 亚洲五月婷婷丁香| 亚洲国产欧美日韩在线播放| 麻豆国产av国片精品| 嫩草影视91久久| 一本—道久久a久久精品蜜桃钙片| 啦啦啦免费观看视频1| 久久午夜综合久久蜜桃| 日韩熟女老妇一区二区性免费视频| 高清欧美精品videossex| 久久天躁狠狠躁夜夜2o2o| 久久婷婷成人综合色麻豆| 在线观看免费午夜福利视频| 丝袜喷水一区| 免费久久久久久久精品成人欧美视频| 日韩欧美一区二区三区在线观看 | 国产精品免费视频内射| 亚洲九九香蕉| 中文字幕av电影在线播放| av视频免费观看在线观看| 国产单亲对白刺激| 亚洲国产成人一精品久久久| 男女边摸边吃奶| 欧美国产精品va在线观看不卡| 久久精品人人爽人人爽视色| xxxhd国产人妻xxx| 欧美午夜高清在线| 大片电影免费在线观看免费| 男女午夜视频在线观看| 18禁黄网站禁片午夜丰满| 国产精品久久久久久精品古装| 又大又爽又粗| 亚洲欧美精品综合一区二区三区| 狠狠婷婷综合久久久久久88av| 99九九在线精品视频|