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

    鋼軌波磨區(qū)段高速輪軌瞬態(tài)滾動(dòng)接觸高頻動(dòng)態(tài)特性

    2018-10-13 03:22:04王衛(wèi)東劉金朝
    中國鐵道科學(xué) 2018年5期
    關(guān)鍵詞:波磨軌枕輪軌

    于 淼,王衛(wèi)東,劉金朝

    (1.中國鐵道科學(xué)研究院,北京 100081;2.中國鐵道科學(xué)研究院 基礎(chǔ)設(shè)施檢測研究所,北京 100081)

    鋼軌波磨是指新鋪設(shè)的鋼軌在使用一段時(shí)間后,在接觸表面沿其縱向面出現(xiàn)的波浪形不均勻磨損,具有明顯的準(zhǔn)周期形態(tài)[1]。高速鐵路鋼軌波磨的波長一般為50~150 mm[2]。當(dāng)列車運(yùn)行速度較高時(shí),鋼軌波磨易引起軌道—車輛系統(tǒng)的劇烈振動(dòng),縮短車輛及軌道部件的使用壽命,增加鐵路養(yǎng)護(hù)維修費(fèi)用,嚴(yán)重時(shí)對(duì)行車安全也是潛在的危害。因此,探究鋼軌波磨機(jī)理以及明確鋼軌波磨引起的軌道—車輛動(dòng)態(tài)響應(yīng)規(guī)律是解決鐵路現(xiàn)場波磨問題的關(guān)鍵。

    從19世紀(jì)末期,國內(nèi)外開始對(duì)鋼軌波磨展開研究。波磨成因復(fù)雜,與多個(gè)因素相關(guān),表現(xiàn)形式也有較大差異。2009年,Grassie S L[3]按照波長固定機(jī)理和損傷機(jī)理將波磨成因分成了6類,其中,由于Pinned-Pinned共振引起的“響軌波磨”,波長一般為25~80 mm。該類型波磨在直線軌道以及高速鐵路曲線區(qū)段時(shí)有發(fā)生,在高速鐵路波磨中占主要部分。目前,對(duì)波磨成因的研究方法包括現(xiàn)場觀測、波磨試驗(yàn)、理論分析、數(shù)值仿真等。其中,仿真方法主要分為多剛體動(dòng)力學(xué)仿真、剛?cè)狁詈戏抡嬉约坝邢拊抡?。將輪軌考慮為柔性體時(shí),多用于輪軌自激振動(dòng)、接觸共振及瞬態(tài)滾動(dòng)接觸分析。B?hmer[4]利用ABAQUS有限元軟件從頻域角度研究了塑性變形對(duì)鋼軌波磨發(fā)展的影響。Gómez[5]利用線性模型和“有限條法”對(duì)波磨發(fā)展趨勢進(jìn)行預(yù)測,計(jì)算了在各個(gè)輪軌固有頻率下波磨的預(yù)測結(jié)果。Correa[6]等利用有限元方法分析了同一轉(zhuǎn)向架不同車輪的位移導(dǎo)納和模態(tài)特性與鋼軌波磨的關(guān)聯(lián)關(guān)系。陳光雄[7]建立了輪對(duì)通過小半徑曲線的穩(wěn)態(tài)有限元模型,研究發(fā)現(xiàn)蠕滑力飽和時(shí)輪軌系統(tǒng)摩擦自激振動(dòng)將引起曲線鋼軌波磨。李霞[8]利用有限元模型研究了軌道振型和共振頻率與鋼軌波磨的內(nèi)在關(guān)聯(lián)。以上利用有限元法建模主要從頻域角度分析鋼軌波磨問題,或利用線性模型對(duì)其進(jìn)行時(shí)域分析。Kalker和Groβ-Thebing[9]指出當(dāng)輪對(duì)運(yùn)動(dòng)波長與輪軌接觸斑長度的比值小于10時(shí)需要考慮輪軌非穩(wěn)態(tài)滾動(dòng)接觸問題。鋼軌短波波磨是典型的輪軌非穩(wěn)態(tài)滾動(dòng)接觸問題之一。因此,建立輪軌非穩(wěn)態(tài)滾動(dòng)接觸模型是研究高速鐵路鋼軌短波波磨問題的關(guān)鍵。近幾年,趙鑫[10-11]利用三維高速輪軌有限元模型從輪軌瞬態(tài)滾動(dòng)接觸角度解釋了波磨在出現(xiàn)后進(jìn)入穩(wěn)定的現(xiàn)象。

    本文針對(duì)高速鐵路鋼軌波磨問題建立輪軌瞬態(tài)滾動(dòng)接觸三維有限元模型,考慮軌道—車輛系統(tǒng)耦合作用以及輪軌真實(shí)幾何,利用隱式和顯式相結(jié)合方法進(jìn)行輪軌動(dòng)態(tài)滾動(dòng)接觸計(jì)算。進(jìn)而分析列車高速通過鋼軌波磨時(shí)的高頻動(dòng)態(tài)特性以及列車通過頻率對(duì)波磨發(fā)展趨勢的影響,為探究鋼軌波磨發(fā)展趨勢以及防治措施提供理論依據(jù)。

    1 高速輪軌瞬態(tài)滾動(dòng)接觸有限元模型

    1.1 模型建立

    針對(duì)高速鐵路軌道—車輛系統(tǒng),采用了LMA型踏面車輪以及帶有波磨區(qū)段的無砟軌道。車輪直徑為860 mm,鋼軌廓形為CN60,設(shè)1/40軌底坡。為了避免鋼軌兩端引起的應(yīng)力波效應(yīng),且考慮軌道扣件整體剛度及阻尼作用,建立的軌道全長為15.21 m,軌枕間距為0.65 m,共24組軌枕,波磨區(qū)段長為1 m。設(shè)置輪對(duì)以第9個(gè)軌枕為起點(diǎn),運(yùn)行里程約為3.7 m,當(dāng)輪軌接觸由靜態(tài)到滾動(dòng)狀態(tài)時(shí)會(huì)出現(xiàn)初始激擾,車輛運(yùn)行速度越高,動(dòng)態(tài)效應(yīng)越強(qiáng),所以在波磨區(qū)段與車輪初始位置間設(shè)置了2.4 m的動(dòng)態(tài)松弛區(qū)。為了將本文模型結(jié)果與文獻(xiàn)[10]相比對(duì),將鋼軌波磨波長設(shè)為0.08 m,波深0.14 mm,縱向及截面幾何參照文獻(xiàn)[10],模型選取的系統(tǒng)結(jié)構(gòu)和材料屬性參數(shù)基本與其一致,見表1,且均不考慮輪軌橫向自由度。文獻(xiàn)[10]利用ANSYS/LS-DYNA有限元軟件建立了單輪單軌有限元模型,為了更接近實(shí)際情況以及下一步開展考慮輪對(duì)橫向自由度的研究,本文利用ABAQUS有限元軟件建立了考慮軌道—車輛系統(tǒng)耦合作用的軌道—輪對(duì)模型如圖1所示。

    利用Hypermesh軟件劃分軌道—輪對(duì)系統(tǒng)有限單元,主要為8節(jié)點(diǎn)縮減積分線性實(shí)體單元,輪軌初始接觸以及鋼軌波磨區(qū)段采用最小尺寸單元為1.2 mm,整個(gè)有限元模型節(jié)點(diǎn)約為2.0×106個(gè)。由于鋼軌波磨主要引起軌道—車輛耦合系統(tǒng)高頻振動(dòng),車輛一系懸掛以上車輛系統(tǒng)振動(dòng)頻率較低,因此,將一系懸掛以上結(jié)構(gòu)設(shè)為施加載荷的質(zhì)點(diǎn)且能夠在垂向自由振動(dòng),考慮慣性作用。在車輛系統(tǒng)結(jié)構(gòu)中,軸箱與車軸通過軸承連接,受力較均勻。在有限元模型中一系懸掛由5根彈簧并聯(lián)模擬,每根彈簧懸掛質(zhì)點(diǎn)重量為1.6 t,輪對(duì)軸重為8 t。鋼軌扣件結(jié)構(gòu)由20根彈簧(4排×5列)并聯(lián)模擬,僅考慮垂向剛度和阻尼,仿真模型中共建立48組扣件系統(tǒng)(一側(cè)鋼軌24組),通過編譯腳本程序自動(dòng)生成。固定CA砂漿底面,軌道兩端施加對(duì)稱邊界條件。

    表1 模型參數(shù)

    圖1 軌道—車輛系統(tǒng)有限元模型

    1.2 計(jì)算方法

    利用隱式和顯式相結(jié)合方法模擬高速輪軌瞬態(tài)滾動(dòng)接觸過程。隱式算法利用Newton迭代且滿足動(dòng)力學(xué)平衡方程,基于面—面接觸算法計(jì)算輪對(duì)在初始位置的輪軌接觸狀態(tài)。初始增量步設(shè)為0.1 s,采用自動(dòng)增量步調(diào)整。由于隱式算法是無條件穩(wěn)定的,所以時(shí)間增量相比于顯式算法要大一些[12]。但是在解決非線性問題時(shí),由于每1個(gè)增量步需要經(jīng)過多次迭代才能滿足給定的容許誤差要求,所以計(jì)算耗時(shí)較長。

    輪對(duì)動(dòng)態(tài)滾動(dòng)階段采用顯式積分方法。ABAQUS/Explicit采用中心差分法計(jì)算接觸問題。對(duì)一般的非線性問題,增量開始時(shí)需要搜索接觸面并判斷接觸狀態(tài)。內(nèi)力和接觸力根據(jù)上一步的狀態(tài)進(jìn)行遞推,因此不需要隱式算法的迭代過程。利用該方法進(jìn)行顯式時(shí)程積分,計(jì)算波磨區(qū)段輪軌受力狀態(tài)。計(jì)算過程主要分為4個(gè)步驟,具體如下。

    步驟1:在t時(shí)刻,求出第i個(gè)增量步的節(jié)點(diǎn)加速度[12],為

    ü(i)=M-1(FP(i)-FI(i))

    (1)

    式中:u為節(jié)點(diǎn)位移;M為節(jié)點(diǎn)質(zhì)量矩陣;FP為外力;FI為內(nèi)力。

    步驟2:對(duì)加速度在時(shí)域積分。在1個(gè)增量步內(nèi),加速度假定為恒值,Δt為時(shí)間增量,則節(jié)點(diǎn)速度為

    (2)

    步驟3:對(duì)速度在時(shí)域積分。在1個(gè)增量步內(nèi),速度假定為恒值,則節(jié)點(diǎn)位移為

    (3)

    步驟4:通過應(yīng)變變化率計(jì)算應(yīng)變增量dε,從而計(jì)算應(yīng)力σ,然后將由式(5)計(jì)算所得的節(jié)點(diǎn)內(nèi)力作為下1個(gè)增量步的初始內(nèi)力,進(jìn)入下1個(gè)增量步計(jì)算。

    σ(i+1)=f(σ(i),dε)

    (4)

    FI(i+1)=BTσ(i+1)

    (5)

    式中:B為應(yīng)變矩陣。

    1.3 模型驗(yàn)證及對(duì)比

    隱式與顯式算法在計(jì)算輪軌接觸時(shí)有所不同,將靜態(tài)結(jié)果作為初始條件導(dǎo)入動(dòng)態(tài)顯式計(jì)算時(shí)會(huì)引入初始激擾。若此時(shí)加載輪對(duì)初始速度會(huì)使初始激擾更加劇烈,延長所需的動(dòng)態(tài)松弛區(qū)長度,增加計(jì)算時(shí)長。為此,在采用隱式算法計(jì)算輪軌靜平衡后,增加采用顯式算法計(jì)算輪軌平衡的過程。

    將利用隱式算法獲得的軌道—車輛系統(tǒng)靜態(tài)平衡時(shí)的結(jié)果作為初始條件,由于彈簧狀態(tài)無法作為初始條件導(dǎo)入,需基于軸重計(jì)算設(shè)置彈簧預(yù)壓力。利用顯式方法計(jì)算輪軌接觸達(dá)到平衡時(shí)的時(shí)域過程如圖2所示。從圖2可以看出:輪軌接觸存在初始激擾,當(dāng)計(jì)算時(shí)長為0.03 s時(shí),輪對(duì)垂向振動(dòng)加速度已經(jīng)趨于0,且輪軌垂向接觸力基本平穩(wěn),約為83.5 kN。因此,將顯式方法計(jì)算輪軌接觸平衡的時(shí)長設(shè)為0.03 s即可滿足計(jì)算輪軌平衡的要求。

    圖2 基于顯式算法的輪軌接觸平衡過程

    設(shè)置輪對(duì)初始運(yùn)行速度為300 km·h-1,包括沿縱向的平動(dòng)速度及繞軸心的轉(zhuǎn)動(dòng)角速度。為避免輪對(duì)在運(yùn)行中一系懸掛作用點(diǎn)出現(xiàn)偏心情況,需添加一系懸掛載荷質(zhì)點(diǎn)與輪對(duì)作用點(diǎn)在縱向及橫向運(yùn)動(dòng)方向上的位移等式約束。在不施加牽引扭矩情況下,輪對(duì)沿鋼軌縱向自由滾動(dòng),輪軌瞬態(tài)滾動(dòng)接觸力如圖3所示。

    圖3 輪軌接觸力

    從圖3可以看出:由于輪對(duì)設(shè)置了初始運(yùn)行速度,所以輪軌接觸力存在初始高頻振動(dòng);輪軌垂向接觸力和縱向接觸力在運(yùn)行里程為1.2 m左右時(shí)收斂趨于平穩(wěn);輪軌縱向接觸力基本為0,車輪與鋼軌接觸為純滾動(dòng),扣件對(duì)輪軌縱向接觸力影響不大;由于軌枕處扣件系統(tǒng)剛度與阻尼作用,離散支撐軌道系統(tǒng)存在固有剛度不平順,在扣件附近輪軌垂向接觸力出現(xiàn)波谷,幅值約為4 kN,且相對(duì)于扣件位置存在約0.2π的相位滯后。

    達(dá)到輪軌接觸平穩(wěn)后的輪對(duì)繼續(xù)沿鋼軌向前滾動(dòng),在2.4 m處進(jìn)入鋼軌波磨區(qū)段,得到光滑鋼軌表面(鋼軌無波磨)和波磨條件下的輪軌接觸力,并與文獻(xiàn)[10]結(jié)果對(duì)比如圖4所示。

    圖4 輪軌垂向接觸力對(duì)比結(jié)果

    從圖4可以看出:本文輪軌滾動(dòng)接觸仿真模型與文獻(xiàn)[10]中輪軌垂向接觸力波動(dòng)的趨勢基本一致。從光滑鋼軌表面處輪軌垂向接觸力可以看出,文獻(xiàn)[10]中結(jié)果在輪軌接觸力平衡時(shí)仍存在峰峰值約為5 kN的高頻振動(dòng),而本文輪軌垂向接觸力相對(duì)平穩(wěn),這是造成文獻(xiàn)[10]的鋼軌波磨區(qū)段輪軌接觸力波動(dòng)大于本文的主要原因之一。由于離散支撐軌道系統(tǒng)存在固有剛度不平順,在扣件附近光滑鋼軌表面輪軌垂向接觸力出現(xiàn)波動(dòng)。文獻(xiàn)[10]中考慮了輪對(duì)附屬部件質(zhì)量為3.3 kN,因此輪軌垂向接觸力均值大于本文結(jié)果。同時(shí),由圖4還可以看出:對(duì)應(yīng)波磨波峰附近輪軌垂向接觸力出現(xiàn)波峰,對(duì)應(yīng)波磨波谷附近輪軌垂向接觸力出現(xiàn)波谷;輪軌垂向接觸力波峰相對(duì)于波磨峰值存在約0.25π的相位超前,進(jìn)而促使鋼軌波磨向車輪滾動(dòng)方向發(fā)展;輪軌垂向接觸力在軌枕附近整體振動(dòng)較大,最大峰值約為126 kN,動(dòng)靜比約為1.5;在軌枕跨間整體振動(dòng)較小,最小峰值約為109 kN,動(dòng)靜比約為1.3。波磨區(qū)段輪軌垂向接觸力整體振型與“拍”振類似,可能存在波磨引起的激振力與軌道系統(tǒng)某固有頻率相接近的情況,即可能發(fā)生了系統(tǒng)共振,形成了時(shí)強(qiáng)時(shí)弱的信號(hào)。為此,下一部分將進(jìn)一步研究軌道—車輛系統(tǒng)可能存在的共振現(xiàn)象。

    2 輪軌瞬態(tài)響應(yīng)分析

    2.1 波磨區(qū)段列車通過頻率分析

    我國高速鐵路運(yùn)營速度為200~350 km·h-1,因此,為了研究波磨區(qū)段列車通過頻率對(duì)軌道—車輛系統(tǒng)動(dòng)態(tài)響應(yīng)的影響,選取列車在100~500 km·h-1運(yùn)行速度下波磨區(qū)段輪軌滾動(dòng)接觸狀態(tài)進(jìn)行對(duì)比分析。當(dāng)列車運(yùn)行速度為v(km·h-1)時(shí),通過波長為λ(m)的波磨區(qū)段所對(duì)應(yīng)的列車通過頻率fw(Hz)為

    (6)

    利用式(6)計(jì)算可得出列車運(yùn)行速度為100~500 km·h-1時(shí),波長為0.08 m波磨區(qū)段的列車通過頻率如圖5所示。從圖5可以看出:列車的運(yùn)行速度與通過頻率成正比,波長為0.08 m的短波波磨引起軌盜—車輛系統(tǒng)高頻振動(dòng)。

    圖5 不同速度下鋼軌波磨區(qū)段列車的通過頻率

    2.2 鋼軌Pinned-Pinned振動(dòng)分析

    鋼軌Pinned-Pinned振動(dòng)是指振型的節(jié)點(diǎn)位于軌枕處,鋼軌似乎被軌枕處的節(jié)點(diǎn)固定,即鋼軌在軌枕作用下的周期性固定模態(tài)振動(dòng)形式。在2個(gè)軌枕跨中處施加激勵(lì),易激發(fā)鋼軌Pinned-Pinned振型,且頻率響應(yīng)存在1個(gè)顯著峰值;在軌枕上方施加激勵(lì)會(huì)出現(xiàn)反諧振振型[11]。文獻(xiàn)[3]中提出鋼軌Pinned-Pinned頻率的計(jì)算方法,APD Man對(duì)其進(jìn)行了簡化,并添加了Pinned-Pinned頻率的模態(tài)階數(shù)系數(shù)n,鋼軌n階Pinned-Pinned頻率fp[13]為

    (7)

    式中:m和EI分別為單位長度鋼軌的質(zhì)量及鋼軌的抗彎曲剛度;L為軌枕或扣件間距。

    通過現(xiàn)場測試和仿真計(jì)算,APD Man進(jìn)一步提出了更接近實(shí)際的Pinned-Pinned頻率估算公式,為

    fp=10.2L-1.61(EI)0.33m-0.33

    (8)

    從式(8)可以看出,除去鋼軌固有屬性外,鋼軌Pinned-Pinned振動(dòng)頻率主要與軌枕間距有關(guān),因此,可以考慮通過縮短軌枕間距來提高鋼軌Pinned-Pinned振動(dòng)頻率。在本文模型中,L=0.65 m,EI=6.62×106N·m2 [14],m=60.643 kg·m-1。則由式(7)可得1階Pinned-Pinned頻率的計(jì)算結(jié)果約為1 228 Hz,由式(8)可得的計(jì)算結(jié)果約為938 Hz。

    對(duì)軌道系統(tǒng)有限元仿真模型進(jìn)行模態(tài)計(jì)算,獲得軌道結(jié)構(gòu)中鋼軌垂向Pinned-Pinned振動(dòng)模態(tài)如圖6所示,圖中放大變形系數(shù)為50。

    圖6 鋼軌垂向Pinned-Pinned振動(dòng)模態(tài)

    通過模態(tài)計(jì)算可得鋼軌垂向Pinned-Pinned頻率為949.43 Hz,波長為扣件間距的2倍,且在2個(gè)扣件之間振幅最大,在扣件處振幅為0。仿真所得模態(tài)結(jié)果即鋼軌垂向Pinned-Pinned頻率與式(8)的計(jì)算結(jié)果相接近,因此式(8)可用于估算鋼軌Pinned-Pinned頻率。無砟軌道扣件間距一般為0.65 m,最小不宜小于0.60 m,由此利用式(8)計(jì)算可得高速鐵路無砟軌道Pinned-Pinned頻率一般約為938~1 067 Hz。結(jié)合式(6)可得當(dāng)列車以270~307 km·h-1速度通過0.08 m左右波長的波磨鋼軌時(shí),鋼軌會(huì)出現(xiàn)Pinned-Pinned共振模態(tài)。由于系統(tǒng)存在較大的阻尼,所以即使在共振區(qū)域,共振頻率的動(dòng)力放大因數(shù)也會(huì)較小,軌道—車輛系統(tǒng)動(dòng)態(tài)響應(yīng)振幅變化不會(huì)過于劇烈。

    2.3 鋼軌Pinned-Pinned振動(dòng)對(duì)輪軌接觸力的影響

    利用三維輪軌瞬態(tài)滾動(dòng)接觸仿真模型,研究鋼軌Pinned-Pinned振動(dòng)對(duì)軌道—車輛系統(tǒng)動(dòng)態(tài)響應(yīng)的影響,包括輪軌垂向接觸力和縱向接觸力。當(dāng)不施加牽引扭矩時(shí),輪對(duì)為純滾運(yùn)動(dòng),輪軌縱向接觸力為0,因而無法研究鋼軌Pinned-Pinned振動(dòng)對(duì)輪軌縱向接觸力的影響。所以,當(dāng)計(jì)算動(dòng)態(tài)輪軌滾動(dòng)接觸時(shí),在輪對(duì)車軸處施加牽引扭矩,模擬列車牽引過程。模型中輪軌摩擦系數(shù)設(shè)為0.5,為避免車輪出現(xiàn)打滑現(xiàn)象,則取牽引系數(shù)為0.3。利用顯式算法計(jì)算輪軌接觸平衡時(shí)的輪軌接觸力可求出所要施加的牽引扭矩大小。為了減緩施加牽引扭矩時(shí)引起的初始輪軌接觸波動(dòng),牽引扭矩在輪對(duì)初始位置時(shí)設(shè)為0,然后成線性增加,在輪對(duì)滾動(dòng)0.005 s時(shí)達(dá)到最大值[10],即牽引系數(shù)為0.3時(shí)對(duì)應(yīng)的牽引扭矩值,約為21 500 N·m。之后,牽引扭矩保持不變。在牽引扭矩作用下,列車運(yùn)行速度如圖7所示。

    圖7 牽引扭矩作用下列車運(yùn)行速度

    從圖7可以看出:由于初始階段牽引扭矩較小,在輪軌摩擦力作用下列車運(yùn)行速度會(huì)略有下降;當(dāng)牽引扭矩保持最大值不變時(shí),列車運(yùn)行速度會(huì)逐漸增加,且在波磨區(qū)段隨波磨幾何不平順呈現(xiàn)輕微周期性波動(dòng)。由于模型中波磨區(qū)段長度為1 m,該區(qū)段列車運(yùn)行速度增加不超過0.2 km·h-1,速度變化相對(duì)較小,因此,在后文的分析中可忽略牽引扭矩引起的速度改變所帶來的影響。

    由第1.3部分可知,輪對(duì)滾動(dòng)基本進(jìn)入穩(wěn)態(tài)后,在通過軌枕上方鋼軌時(shí)由軌道離散支撐作用引起的剛度不平順會(huì)導(dǎo)致輪軌接觸力出現(xiàn)波動(dòng)。當(dāng)施加牽引扭矩時(shí),列車以不同速度運(yùn)行,輪對(duì)在未進(jìn)入鋼軌波磨區(qū)段前,不同速度下的輪軌垂向接觸力如圖8所示。

    從圖8可以看出:低速時(shí)輪軌垂向接觸力振動(dòng)收斂較迅速,隨著列車運(yùn)行速度的提高,輪軌垂向接觸力初始高頻振動(dòng)將愈發(fā)劇烈;在不同運(yùn)行速度下,列車運(yùn)行約1.5 m后基本進(jìn)入平穩(wěn)狀態(tài);當(dāng)輪對(duì)通過2 m附近處軌枕上方鋼軌時(shí),輪軌垂向接觸力再次出現(xiàn)波動(dòng)。當(dāng)運(yùn)行速度為100 km·h-1時(shí),輪軌垂向接觸力波動(dòng)不明顯。隨著運(yùn)行速度的提高,由于軌道—車輛系統(tǒng)剛度和阻尼的作用,扣件附近輪軌垂向接觸力出現(xiàn)波動(dòng)越明顯,且輪軌垂向接觸力逐漸減小,相對(duì)于扣件位置存在相位滯后。

    圖8 鋼軌非波磨區(qū)段不同速度下的輪軌垂向接觸力

    達(dá)到輪軌接觸平穩(wěn)后的輪對(duì)繼續(xù)沿鋼軌向前滾動(dòng),在2.4 m處進(jìn)入鋼軌波磨區(qū)段。在牽引扭矩作用下,列車以不同速度運(yùn)行時(shí)波磨區(qū)段輪軌垂向接觸力如圖9所示。

    圖9 鋼軌波磨區(qū)段不同速度下的輪軌垂向接觸力

    從圖9可以看出:當(dāng)列車運(yùn)行速度為100 km·h-1時(shí),鋼軌波磨區(qū)段輪軌垂向接觸力隨波磨幾何不平順基本成周期等幅振動(dòng),受軌枕處扣件支撐作用影響較小。對(duì)應(yīng)波磨波峰附近的輪軌垂向接觸力出現(xiàn)波峰;對(duì)應(yīng)波磨波谷附近的輪軌垂向接觸力出現(xiàn)波谷。輪軌垂向接觸力波峰相對(duì)于波磨峰值存在相位超前,進(jìn)而促使鋼軌波磨向車輪滾動(dòng)方向發(fā)展。隨著列車運(yùn)行速度的提高,輪軌垂向接觸力振幅整體呈上升趨勢。當(dāng)速度為200~300 km·h-1時(shí),輪軌垂向接觸力最大值出現(xiàn)在軌枕上方附近,2個(gè)軌枕跨中處迅速減小;當(dāng)速度為400~500 km·h-1時(shí),輪軌垂向接觸力最大值出現(xiàn)在2個(gè)軌枕之間中后部分,在軌枕上方較小,與前面提到的鋼軌離散支撐規(guī)律一致,是由于軌枕處扣件系統(tǒng)剛度及阻尼作用。當(dāng)速度為300 km·h-1時(shí),對(duì)應(yīng)列車通過頻率為1 042 Hz,相對(duì)接近于鋼軌1階垂向Pinned-Pinned振動(dòng)頻率,輪軌垂向接觸力出現(xiàn)“拍”振形式,軌道—車輛系統(tǒng)出現(xiàn)了Pinned-Pinned共振。由于軌道系統(tǒng)存在較大的阻尼,所以即使在共振區(qū),共振頻率的動(dòng)力放大因數(shù)也會(huì)相對(duì)較小,振幅變化不會(huì)過于劇烈。

    當(dāng)列車運(yùn)行速度分別為300和400 km·h-1時(shí),利用Wigner-Ville分布對(duì)輪軌垂向接觸力進(jìn)行時(shí)頻分析,并進(jìn)行500 Hz高通濾波后的結(jié)果如圖10所示。

    圖10 輪軌垂向接觸力時(shí)頻圖

    從圖10可以看出:當(dāng)列車運(yùn)行速度為300 km·h-1時(shí),主頻約為1 042 Hz的輪軌垂向接觸力振動(dòng)能量主要集中在2.6 m處以及3.25 m處的軌枕上方附近,Pinned-Pinned共振將造成輪軌垂向接觸力較大波動(dòng),會(huì)加速扣件系統(tǒng)傷損或疲勞斷裂;當(dāng)列車運(yùn)行速度為400 km·h-1時(shí),主頻約為1 389 Hz的輪軌垂向接觸力振動(dòng)能量主要集中在2個(gè)軌枕之間中后部,約3 m附近。

    列車以不同速度運(yùn)行時(shí)輪軌縱向接觸力如圖11所示。

    圖11 不同速度下鋼軌波磨區(qū)段的輪軌縱向接觸力

    從圖11可以看出:列車以不同速度通過鋼軌波磨區(qū)段時(shí),輪軌縱向接觸力隨波磨幾何不平順呈周期性波動(dòng),對(duì)應(yīng)波磨波峰附近的輪軌縱向接觸力出現(xiàn)波峰,對(duì)應(yīng)波磨波谷附近的輪軌縱向接觸力出現(xiàn)波谷,且相位略有差異。當(dāng)速度為300 km·h-1時(shí),軌枕上方附近的輪軌縱向接觸力明顯大于其他速度結(jié)果,且“拍”振形式較明顯,即與鋼軌Pinned-Pinned振動(dòng)發(fā)生了共振,輪軌縱向接觸力最大峰值主要出現(xiàn)在軌枕上方后0.2 m之內(nèi),約為2個(gè)軌枕跨間最小峰值的1.3倍。由于軌道系統(tǒng)阻尼較大,共振現(xiàn)象不會(huì)過于劇烈。對(duì)比圖9和圖11可以看出,Pinned-Pinned共振對(duì)輪軌縱向接觸力影響相比于垂向接觸力更加明顯。

    在鋼軌波磨區(qū)段,車軸處施加的牽引扭矩載荷為定值,是由牽引系數(shù)0.3計(jì)算獲得。而當(dāng)列車通過鋼軌波磨區(qū)段時(shí),由于輪軌非穩(wěn)態(tài)滾動(dòng)接觸作用,當(dāng)牽引扭矩為定值時(shí),牽引比(輪軌縱向接觸力與垂向接觸力的比值)[3]出現(xiàn)了波動(dòng)。計(jì)算速度為300 km·h-1時(shí)波磨區(qū)段的牽引比如圖12所示。

    圖12 波磨區(qū)段牽引比

    從圖12可以看出:波磨區(qū)段牽引比出現(xiàn)明顯波動(dòng),波動(dòng)范圍在0.20~0.48之間,波峰均未超過輪軌摩擦系數(shù)0.5,說明輪軌接觸既沒有發(fā)生純滾也沒有出現(xiàn)全滑動(dòng)現(xiàn)象。輪軌接觸出現(xiàn)了周期性黏滑振動(dòng),進(jìn)而導(dǎo)致鋼軌表面的不均勻磨耗[1]。與輪軌垂向接觸力和縱向接觸力分布不同,牽引比與波磨幾何不平順基本上為反相位,對(duì)應(yīng)波磨波谷附近的牽引比出現(xiàn)波峰;而對(duì)應(yīng)波磨波峰附近的牽引比出現(xiàn)波谷。波磨區(qū)段的牽引比在軌枕附近整體波動(dòng)較大,最大峰值約為0.47,接近于輪軌摩擦系數(shù)0.5,則波磨波谷處輪軌接觸斑內(nèi)滑動(dòng)區(qū)與黏著區(qū)的比值較大,將會(huì)加劇波磨波谷處滑動(dòng)磨損,而波磨波峰處牽引比較小,作用不明顯。牽引比在2個(gè)軌枕跨間整體波動(dòng)較小,最小波峰約為0.33。因此,列車通過波磨區(qū)段時(shí),在牽引扭矩作用下,輪軌接觸存在周期性黏滑振動(dòng),當(dāng)軌道—車輛系統(tǒng)出現(xiàn)Pinned-Pinned共振時(shí),會(huì)促使軌枕附近鋼軌波磨波谷處的滑動(dòng)磨損。

    3 結(jié)論及展望

    (1)列車以不同速度通過鋼軌波磨區(qū)段時(shí),輪軌垂向接觸力和縱向接觸力隨波磨幾何不平順呈周期性波動(dòng),相位略有差異,對(duì)應(yīng)波磨波峰附近的輪軌接觸力出現(xiàn)波峰,對(duì)應(yīng)波磨波谷附近的輪軌接觸力出現(xiàn)波谷;而牽引比波動(dòng)與波磨幾何不平順基本呈現(xiàn)反相位。

    (2)當(dāng)列車通過波磨區(qū)段所對(duì)應(yīng)的通過頻率與軌道Pinned-Pinned頻率相近時(shí)引起系統(tǒng)共振,輪軌接觸力出現(xiàn)“拍”振特性;輪軌接觸力在軌枕附近整體波動(dòng)較大,在2個(gè)軌枕跨間整體波動(dòng)較小,會(huì)加速鋼軌扣件傷損。

    (3)列車通過波磨區(qū)段時(shí),在牽引扭矩作用下,輪軌接觸存在周期性黏滑振動(dòng)。當(dāng)軌道—車輛系統(tǒng)出現(xiàn)Pinned-Pinned共振時(shí),波磨區(qū)段的牽引比在軌枕附近整體波動(dòng)較大,波磨波谷處牽引比接近于輪軌摩擦系數(shù),將會(huì)加劇波磨波谷處滑動(dòng)磨損,加快波磨的發(fā)展。

    (4)由于本文主要針對(duì)直線軌道鋼軌波磨問題,且兩側(cè)鋼軌采用相同波磨幾何不平順,因此輪軌橫向蠕滑作用較小,同時(shí),為了更好地與文獻(xiàn)[10]結(jié)果相比對(duì),本文在建立高速輪軌非穩(wěn)態(tài)滾動(dòng)接觸模型時(shí)約束了輪軌橫向自由度。在后續(xù)工作中,解決曲線波磨問題以及兩側(cè)鋼軌存在不同波磨幾何不平順時(shí)將進(jìn)一步研究橫向蠕滑對(duì)鋼軌波磨形成及發(fā)展的影響。

    猜你喜歡
    波磨軌枕輪軌
    城軌小半徑曲線鋼軌波磨的成因與整治研究綜述
    裝配式可調(diào)間距的X形抗滑動(dòng)軌枕研究
    北京地鐵鋼軌波磨現(xiàn)狀及防治措施
    中低速磁浮道岔與輪軌道岔的差異
    軌枕外觀及修補(bǔ)工藝研究
    中低速磁浮與輪軌交通信號(hào)系統(tǒng)的差異
    非線性穩(wěn)態(tài)曲線通過時(shí)輪軌滾動(dòng)接觸的數(shù)值求解方法
    高速鐵路鋼軌波磨對(duì)車輛—軌道動(dòng)態(tài)響應(yīng)的影響
    一種新的混凝土軌枕裝卸機(jī)設(shè)計(jì)
    減緩高速鐵路鋼軌波磨的仿真分析
    国产成+人综合+亚洲专区| 国内精品一区二区在线观看| 黄色配什么色好看| 亚洲自拍偷在线| 深夜精品福利| 色精品久久人妻99蜜桃| aaaaa片日本免费| 国产av一区在线观看免费| 日日摸夜夜添夜夜添av毛片 | 国产精品久久久久久久电影| 精品一区二区免费观看| 日韩av在线大香蕉| 女人被狂操c到高潮| 99热精品在线国产| 欧美日韩综合久久久久久 | 欧美午夜高清在线| 嫁个100分男人电影在线观看| 成人亚洲精品av一区二区| 国产精品久久电影中文字幕| 亚洲男人的天堂狠狠| 尤物成人国产欧美一区二区三区| 亚洲精品一卡2卡三卡4卡5卡| 亚洲成人中文字幕在线播放| 91在线精品国自产拍蜜月| aaaaa片日本免费| 国产激情偷乱视频一区二区| 嫩草影视91久久| 欧美精品国产亚洲| 国产人妻一区二区三区在| 国产 一区 欧美 日韩| 在线十欧美十亚洲十日本专区| 欧美日韩黄片免| 欧美午夜高清在线| 一本久久中文字幕| 男女那种视频在线观看| 波野结衣二区三区在线| 欧美bdsm另类| 亚洲五月天丁香| 午夜免费激情av| 亚洲不卡免费看| 看黄色毛片网站| 亚洲一区二区三区色噜噜| 亚洲国产欧洲综合997久久,| 亚洲最大成人av| av女优亚洲男人天堂| 老女人水多毛片| 在线a可以看的网站| 最新中文字幕久久久久| 在线观看午夜福利视频| 国产色爽女视频免费观看| 亚洲精品一卡2卡三卡4卡5卡| 久久精品久久久久久噜噜老黄 | 很黄的视频免费| 怎么达到女性高潮| 精品一区二区三区人妻视频| 久久精品综合一区二区三区| 久9热在线精品视频| 老女人水多毛片| 亚洲人成网站在线播放欧美日韩| 国产私拍福利视频在线观看| 欧美激情在线99| 97热精品久久久久久| 国产精品嫩草影院av在线观看 | a级毛片免费高清观看在线播放| 欧美性猛交黑人性爽| 黄色视频,在线免费观看| 午夜老司机福利剧场| 老司机午夜福利在线观看视频| 我的老师免费观看完整版| 亚洲美女黄片视频| 2021天堂中文幕一二区在线观| 日韩中字成人| 哪里可以看免费的av片| 久久午夜福利片| 精品欧美国产一区二区三| 男插女下体视频免费在线播放| 99热6这里只有精品| 亚洲在线自拍视频| 亚洲最大成人中文| 成人特级黄色片久久久久久久| av在线蜜桃| 我的老师免费观看完整版| 91av网一区二区| 午夜影院日韩av| 国产免费一级a男人的天堂| www日本黄色视频网| 美女黄网站色视频| a级毛片a级免费在线| 91av网一区二区| 老女人水多毛片| 午夜免费男女啪啪视频观看 | 搡老岳熟女国产| 久久久成人免费电影| www.999成人在线观看| 乱人视频在线观看| 久久亚洲真实| 日韩欧美精品免费久久 | 欧美不卡视频在线免费观看| 免费av不卡在线播放| 永久网站在线| 午夜免费成人在线视频| 国产一区二区在线av高清观看| 听说在线观看完整版免费高清| 欧美黑人巨大hd| 午夜精品在线福利| 国产中年淑女户外野战色| 欧美另类亚洲清纯唯美| 午夜精品一区二区三区免费看| 18+在线观看网站| 国语自产精品视频在线第100页| 搡老岳熟女国产| 亚洲国产精品成人综合色| 亚洲人成网站在线播| 国产精品98久久久久久宅男小说| 天堂动漫精品| 老司机午夜十八禁免费视频| 看十八女毛片水多多多| 国产精品一区二区三区四区久久| 国内精品久久久久久久电影| 黄色一级大片看看| 亚洲激情在线av| 高清毛片免费观看视频网站| 中文字幕高清在线视频| 一a级毛片在线观看| 免费搜索国产男女视频| 国产在视频线在精品| 国语自产精品视频在线第100页| 男人舔女人下体高潮全视频| 午夜免费激情av| av黄色大香蕉| 不卡一级毛片| 久久久久久大精品| 成年女人毛片免费观看观看9| 国产探花极品一区二区| 九九久久精品国产亚洲av麻豆| 69av精品久久久久久| 99久国产av精品| 亚洲最大成人手机在线| 亚洲一区高清亚洲精品| 国产精品国产高清国产av| 搡老岳熟女国产| 日日干狠狠操夜夜爽| 国产v大片淫在线免费观看| 内射极品少妇av片p| 亚洲美女黄片视频| 中亚洲国语对白在线视频| 色尼玛亚洲综合影院| 亚洲av二区三区四区| 国产高清激情床上av| 午夜福利在线观看吧| 黄色丝袜av网址大全| 欧美最新免费一区二区三区 | 久久国产精品人妻蜜桃| 欧美日本视频| 美女大奶头视频| 最近中文字幕高清免费大全6 | 嫁个100分男人电影在线观看| 国产又黄又爽又无遮挡在线| 日韩欧美国产一区二区入口| 91在线精品国自产拍蜜月| 久久伊人香网站| 永久网站在线| 国产黄色小视频在线观看| 久久精品国产亚洲av天美| 九九在线视频观看精品| av在线天堂中文字幕| 午夜福利成人在线免费观看| av在线蜜桃| 极品教师在线免费播放| 国产色爽女视频免费观看| 成人毛片a级毛片在线播放| 极品教师在线免费播放| 久久国产精品人妻蜜桃| 精品久久久久久,| 91在线精品国自产拍蜜月| 亚洲天堂国产精品一区在线| 久久久久久九九精品二区国产| 欧美性猛交╳xxx乱大交人| 精品午夜福利视频在线观看一区| 一个人看视频在线观看www免费| 久久精品人妻少妇| 国产精品一及| а√天堂www在线а√下载| 桃红色精品国产亚洲av| 日本黄色片子视频| 欧美潮喷喷水| 色哟哟哟哟哟哟| 国内揄拍国产精品人妻在线| 亚洲最大成人av| 久9热在线精品视频| 久久伊人香网站| 国产探花极品一区二区| 亚洲精品在线观看二区| 国产欧美日韩精品亚洲av| 最近最新免费中文字幕在线| 男人狂女人下面高潮的视频| 国产激情偷乱视频一区二区| 午夜日韩欧美国产| 欧美一区二区精品小视频在线| 一a级毛片在线观看| 麻豆国产97在线/欧美| 国内毛片毛片毛片毛片毛片| 一本久久中文字幕| 又爽又黄无遮挡网站| 亚洲av中文字字幕乱码综合| 老司机深夜福利视频在线观看| 美女高潮的动态| 国产在线精品亚洲第一网站| 又紧又爽又黄一区二区| 亚洲自拍偷在线| 舔av片在线| 日韩欧美精品免费久久 | 国产精品人妻久久久久久| 亚洲一区二区三区色噜噜| 国产欧美日韩精品亚洲av| 久久久色成人| 丁香六月欧美| 久久欧美精品欧美久久欧美| 国产精品一及| 色播亚洲综合网| 可以在线观看的亚洲视频| 真实男女啪啪啪动态图| 中文亚洲av片在线观看爽| 日韩欧美一区二区三区在线观看| 最近最新免费中文字幕在线| 精品久久久久久成人av| 成人国产综合亚洲| 一级a爱片免费观看的视频| 欧美一区二区亚洲| 99热这里只有是精品在线观看 | 久久精品综合一区二区三区| 精品熟女少妇八av免费久了| 亚洲三级黄色毛片| 欧美在线黄色| 国产亚洲欧美98| 亚洲 欧美 日韩 在线 免费| 免费搜索国产男女视频| 免费av观看视频| 美女黄网站色视频| 免费黄网站久久成人精品 | 人妻久久中文字幕网| 成人av一区二区三区在线看| 欧美xxxx黑人xx丫x性爽| 久久久久久久久中文| 国产精品久久久久久精品电影| 免费在线观看亚洲国产| 熟女电影av网| 69人妻影院| a级一级毛片免费在线观看| 婷婷六月久久综合丁香| 国产亚洲欧美98| 身体一侧抽搐| www.www免费av| 国产精品亚洲美女久久久| 十八禁网站免费在线| 一卡2卡三卡四卡精品乱码亚洲| 91字幕亚洲| www.www免费av| 欧美绝顶高潮抽搐喷水| 精品久久久久久久久av| 亚洲在线观看片| eeuss影院久久| 国产成人aa在线观看| 又紧又爽又黄一区二区| 啪啪无遮挡十八禁网站| 亚洲,欧美,日韩| www.999成人在线观看| 免费黄网站久久成人精品 | 丰满人妻熟妇乱又伦精品不卡| 午夜免费成人在线视频| 亚洲成av人片免费观看| 午夜两性在线视频| 久久6这里有精品| 小蜜桃在线观看免费完整版高清| 少妇的逼好多水| 在线播放国产精品三级| 久久精品久久久久久噜噜老黄 | 国产亚洲精品av在线| 亚洲,欧美精品.| 欧美另类亚洲清纯唯美| 3wmmmm亚洲av在线观看| 97超级碰碰碰精品色视频在线观看| 亚洲精品在线美女| 免费人成视频x8x8入口观看| netflix在线观看网站| 97超视频在线观看视频| 中文在线观看免费www的网站| 亚洲欧美日韩无卡精品| 欧美色欧美亚洲另类二区| 色哟哟哟哟哟哟| 一区二区三区免费毛片| 黄色日韩在线| 男女视频在线观看网站免费| 黄色一级大片看看| 岛国在线免费视频观看| 久久精品国产亚洲av涩爱 | 伦理电影大哥的女人| 一个人免费在线观看的高清视频| 精品久久久久久久末码| 99久久九九国产精品国产免费| a级毛片a级免费在线| 国内久久婷婷六月综合欲色啪| 亚洲成人久久性| 国产91精品成人一区二区三区| 女人十人毛片免费观看3o分钟| 51国产日韩欧美| 成人美女网站在线观看视频| 欧美日韩黄片免| 久久久久久国产a免费观看| 国产成人欧美在线观看| 国产精品98久久久久久宅男小说| 亚洲av成人av| 久久99热6这里只有精品| 欧美不卡视频在线免费观看| 国产精品一及| 51国产日韩欧美| 成人国产综合亚洲| 五月玫瑰六月丁香| 午夜激情福利司机影院| 国产毛片a区久久久久| 亚洲av免费在线观看| 日本撒尿小便嘘嘘汇集6| h日本视频在线播放| 亚洲精品粉嫩美女一区| 如何舔出高潮| 12—13女人毛片做爰片一| 欧美日韩乱码在线| 天堂网av新在线| 亚洲精品成人久久久久久| 欧美成人一区二区免费高清观看| 少妇熟女aⅴ在线视频| 中文字幕熟女人妻在线| 欧美日本视频| 超碰av人人做人人爽久久| 中文字幕高清在线视频| 露出奶头的视频| 精品国内亚洲2022精品成人| 国产av不卡久久| 中文在线观看免费www的网站| 热99在线观看视频| 亚洲美女搞黄在线观看 | avwww免费| 国产精品嫩草影院av在线观看 | 欧美日韩综合久久久久久 | 三级男女做爰猛烈吃奶摸视频| 亚洲av美国av| 欧美色欧美亚洲另类二区| 小说图片视频综合网站| 久久久久久久久大av| 国产精品久久久久久人妻精品电影| 老女人水多毛片| 国产老妇女一区| 免费av毛片视频| 国产日本99.免费观看| 国产视频一区二区在线看| 成人特级av手机在线观看| 久久伊人香网站| 日本五十路高清| 99riav亚洲国产免费| 91在线观看av| 丰满人妻一区二区三区视频av| 久久久久精品国产欧美久久久| 草草在线视频免费看| 亚洲久久久久久中文字幕| 最新中文字幕久久久久| 高清在线国产一区| 97超视频在线观看视频| 欧美最黄视频在线播放免费| 成年女人永久免费观看视频| 国产在线男女| 噜噜噜噜噜久久久久久91| 亚洲av中文字字幕乱码综合| 久久久精品欧美日韩精品| 国产精品人妻久久久久久| 淫秽高清视频在线观看| 一本综合久久免费| 亚洲av熟女| 黄色一级大片看看| 一区二区三区免费毛片| 精品国产亚洲在线| 男人的好看免费观看在线视频| 精品不卡国产一区二区三区| 成年人黄色毛片网站| 亚洲电影在线观看av| 久久久久九九精品影院| 欧美日韩黄片免| 最近在线观看免费完整版| 亚洲成人久久性| 欧美日韩福利视频一区二区| 国产精品久久久久久亚洲av鲁大| 99久久精品国产亚洲精品| av中文乱码字幕在线| aaaaa片日本免费| 2021天堂中文幕一二区在线观| 最新在线观看一区二区三区| 人人妻,人人澡人人爽秒播| 老司机午夜十八禁免费视频| 少妇人妻一区二区三区视频| 蜜桃久久精品国产亚洲av| 久久婷婷人人爽人人干人人爱| 99久久成人亚洲精品观看| 我的女老师完整版在线观看| 精品久久久久久久久久免费视频| 麻豆国产97在线/欧美| 亚洲中文日韩欧美视频| 天美传媒精品一区二区| 欧美激情久久久久久爽电影| 最近最新免费中文字幕在线| 亚洲 欧美 日韩 在线 免费| 久99久视频精品免费| 成人国产一区最新在线观看| 两人在一起打扑克的视频| 成人三级黄色视频| 亚洲国产精品久久男人天堂| 88av欧美| 国产激情偷乱视频一区二区| 亚洲av电影在线进入| 中文资源天堂在线| 91久久精品国产一区二区成人| 在线观看66精品国产| 日韩高清综合在线| 悠悠久久av| 热99re8久久精品国产| 免费看a级黄色片| 亚洲第一欧美日韩一区二区三区| 国内揄拍国产精品人妻在线| 国产高清三级在线| 琪琪午夜伦伦电影理论片6080| 中文字幕av在线有码专区| 国内毛片毛片毛片毛片毛片| 国产乱人伦免费视频| 俺也久久电影网| 亚洲成人久久性| 国产中年淑女户外野战色| 午夜福利18| 精品久久久久久成人av| 亚洲av免费在线观看| 99精品在免费线老司机午夜| a级毛片a级免费在线| 国产免费男女视频| 欧美xxxx黑人xx丫x性爽| 999久久久精品免费观看国产| 别揉我奶头~嗯~啊~动态视频| 久久久久亚洲av毛片大全| 免费看日本二区| 日本 av在线| 天堂影院成人在线观看| 天堂√8在线中文| 赤兔流量卡办理| 偷拍熟女少妇极品色| 人妻制服诱惑在线中文字幕| 日韩高清综合在线| 日韩人妻高清精品专区| 亚洲第一欧美日韩一区二区三区| 91麻豆精品激情在线观看国产| 99热这里只有是精品50| 淫秽高清视频在线观看| 欧美激情久久久久久爽电影| 国产亚洲精品久久久com| 又粗又爽又猛毛片免费看| 美女高潮的动态| 热99在线观看视频| 国产探花在线观看一区二区| 一级毛片久久久久久久久女| 三级男女做爰猛烈吃奶摸视频| 日韩亚洲欧美综合| 亚洲乱码一区二区免费版| 中文字幕免费在线视频6| 99热只有精品国产| 人人妻人人看人人澡| 日日干狠狠操夜夜爽| 久久香蕉精品热| 亚洲成av人片在线播放无| 中文字幕熟女人妻在线| 亚洲色图av天堂| 嫁个100分男人电影在线观看| 国产乱人视频| 色综合欧美亚洲国产小说| 国产精品一及| 9191精品国产免费久久| 亚洲精华国产精华精| 1000部很黄的大片| 亚洲精品456在线播放app | 午夜免费激情av| 亚洲 国产 在线| 欧美成人免费av一区二区三区| 我的女老师完整版在线观看| 日韩欧美在线乱码| 嫩草影院入口| 日韩人妻高清精品专区| 校园春色视频在线观看| 亚洲熟妇中文字幕五十中出| 熟妇人妻久久中文字幕3abv| 嫁个100分男人电影在线观看| 国产三级中文精品| 亚洲av不卡在线观看| 亚洲av免费在线观看| 亚洲成av人片在线播放无| 国产精品久久久久久人妻精品电影| 亚洲中文字幕一区二区三区有码在线看| 免费在线观看成人毛片| 男女做爰动态图高潮gif福利片| 亚洲欧美日韩高清在线视频| 90打野战视频偷拍视频| 亚洲经典国产精华液单 | 亚洲精品成人久久久久久| 真实男女啪啪啪动态图| 深爱激情五月婷婷| 中文亚洲av片在线观看爽| 午夜激情福利司机影院| 黄色女人牲交| 99久久精品一区二区三区| 岛国在线免费视频观看| av在线蜜桃| 成年女人毛片免费观看观看9| 黄片小视频在线播放| 欧美日本视频| 99riav亚洲国产免费| 婷婷精品国产亚洲av在线| 亚洲国产精品久久男人天堂| 色综合站精品国产| 老熟妇乱子伦视频在线观看| 久99久视频精品免费| 男人舔女人下体高潮全视频| 亚洲人成电影免费在线| 搡老岳熟女国产| 国产成人a区在线观看| 国内毛片毛片毛片毛片毛片| 夜夜看夜夜爽夜夜摸| 免费在线观看成人毛片| 国产三级在线视频| 色视频www国产| 中文字幕人成人乱码亚洲影| 91麻豆av在线| av国产免费在线观看| 特大巨黑吊av在线直播| 3wmmmm亚洲av在线观看| 国产精品自产拍在线观看55亚洲| 又爽又黄无遮挡网站| 色噜噜av男人的天堂激情| 女人被狂操c到高潮| 舔av片在线| 一本综合久久免费| 欧美一区二区亚洲| 成人国产综合亚洲| 亚洲精品一区av在线观看| 最后的刺客免费高清国语| 中文亚洲av片在线观看爽| 亚洲精华国产精华精| 我的老师免费观看完整版| 桃红色精品国产亚洲av| 国产爱豆传媒在线观看| 桃色一区二区三区在线观看| 国产一区二区在线观看日韩| 乱码一卡2卡4卡精品| 久久久久久久久久黄片| 亚洲在线观看片| 国产精品日韩av在线免费观看| netflix在线观看网站| 日日干狠狠操夜夜爽| 色av中文字幕| 久久中文看片网| 波多野结衣巨乳人妻| 国产高清有码在线观看视频| 好看av亚洲va欧美ⅴa在| 夜夜爽天天搞| 亚洲成人免费电影在线观看| 日韩亚洲欧美综合| 中文字幕av成人在线电影| 日本撒尿小便嘘嘘汇集6| 91久久精品电影网| 免费大片18禁| 在线观看66精品国产| 成人av在线播放网站| 婷婷精品国产亚洲av在线| 首页视频小说图片口味搜索| 日韩欧美一区二区三区在线观看| 看免费av毛片| 久久精品综合一区二区三区| 免费看光身美女| 日本一二三区视频观看| 欧美乱色亚洲激情| 最近中文字幕高清免费大全6 | 男人的好看免费观看在线视频| 国产高清视频在线播放一区| 日本黄色片子视频| 日本 av在线| 欧美最黄视频在线播放免费| 观看美女的网站| 午夜两性在线视频| 美女 人体艺术 gogo| 午夜a级毛片| 免费观看精品视频网站| 亚洲av一区综合| .国产精品久久| 麻豆国产av国片精品| 久久国产乱子伦精品免费另类| 亚洲国产精品久久男人天堂| 亚洲18禁久久av| 在线国产一区二区在线| 欧美性感艳星| 欧美日韩综合久久久久久 | 亚洲成av人片免费观看| av天堂在线播放| 成人亚洲精品av一区二区| 精品午夜福利在线看| 国产午夜精品论理片| 性插视频无遮挡在线免费观看| 三级男女做爰猛烈吃奶摸视频| 中国美女看黄片| 悠悠久久av| 两人在一起打扑克的视频| eeuss影院久久| 又爽又黄a免费视频| 长腿黑丝高跟|