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

    地鐵小半徑曲線鋼彈簧浮置板軌道鋼軌波磨成因分析

    2024-01-11 15:41:00王志強雷震宇
    中國機械工程 2023年24期
    關(guān)鍵詞:波磨浮置輪軌

    王志強 雷震宇

    同濟大學(xué)鐵道與城市軌道交通研究院,上海,201804

    0 引言

    城市軌道交通的不斷發(fā)展使得人們對地鐵線路振動和噪聲水平的上限要求越來越嚴(yán)格。目前,國內(nèi)許多地鐵線路采用了各種各樣的新型軌道來實現(xiàn)減振降噪的目的,諸如套靴軌枕軌道、梯形軌枕軌道和浮置板軌道等。雖然改變軌道結(jié)構(gòu)形式在一定程度上緩解了振動噪聲問題,但現(xiàn)場調(diào)研卻發(fā)現(xiàn)減振軌道上出現(xiàn)了廣泛的鋼軌波磨現(xiàn)象[1-3]。鋼軌波磨發(fā)生在鋼軌表面,呈現(xiàn)為周期性的波狀磨耗。鋼軌波磨的存在容易誘發(fā)輪軌系統(tǒng)產(chǎn)生高頻振動和惱人噪聲,這有悖于減振軌道的最初目的,并且還會造成車輛軌道部件的疲勞失效[4]。如何有效預(yù)防鋼軌波磨一直都是鐵路行業(yè)的重要研究課題,而理解鋼軌波磨的成因機制則是實現(xiàn)控制鋼軌波磨的首要前提。

    針對鋼軌波磨形成原因,國內(nèi)外研究人員從輪軌系統(tǒng)的不同角度展開了探索。陳光雄等[5]提出了輪軌系統(tǒng)摩擦自激振動誘發(fā)鋼軌波磨的觀點,并結(jié)合線路測試驗證了上述觀點的有效性[6-7]。之后,CUI等[8]進一步拓展了摩擦自激振動理論,并研究了摩擦自激振動和軌面不平順反饋振動的聯(lián)合作用對鋼軌波磨的影響。結(jié)合現(xiàn)場測量和數(shù)值仿真,LIU等[9-10]發(fā)現(xiàn)曲線軌道上產(chǎn)生的波磨可歸因于摩擦誘發(fā)的扭轉(zhuǎn)振動引起的摩擦功率的周期性波動,并且受扣件影響的垂向動力學(xué)決定了扭轉(zhuǎn)振動的頻率,進而造成了實測線路上波磨波長隨扣件類型變化的現(xiàn)象。摩擦自激觀點能夠解釋現(xiàn)有地鐵系統(tǒng)中的大部分鋼軌波磨現(xiàn)象,并且能夠精確預(yù)測鋼軌波磨的波長和位置,但在表征波磨演化的輪軌微觀接觸方面還需進一步補充和完善。鋼軌波磨種類繁多,影響因素復(fù)雜,因而使用一種理論并不能完全描述其形成機制。早在1993年,根據(jù)波長固定機理和損傷機理,GRASSIE等[11]將當(dāng)時鐵路系統(tǒng)中的鋼軌波磨現(xiàn)象劃分為六類,分別為pinned-pinned共振、車轍、其他P2共振、重載、輕軌和特定軌道形式波磨,并分別給出了相應(yīng)的治理對策?;诖?大量針對性的案例研究在波磨領(lǐng)域涌現(xiàn)。通過建立輪軌非穩(wěn)態(tài)滾動接觸有限元模型,YU等[12]研究了輪軌接觸對短波波磨的影響以及列車通過頻率對系統(tǒng)動力學(xué)的影響,結(jié)果發(fā)現(xiàn)當(dāng)列車通過頻率接近鋼軌pinned-pinned頻率時,pinned-pinned共振會發(fā)生,并將促進扣件附近波谷處的滑動磨耗。關(guān)慶華等[13]詳細研究了車輛軌道系統(tǒng)P2共振的影響因素,并結(jié)合試驗提出和驗證了通過軌道振動屬性反推P2共振頻率的方法。金鋒等[14]對神朔重載鐵路波磨地段進行了長期跟蹤測試,發(fā)現(xiàn)其成因與有砟軌道P2共振頻率有關(guān),并將該線路上的波磨發(fā)展總結(jié)為初期萌生、中期快速發(fā)展和末期穩(wěn)定三個階段。此外,SUN等[15-16]、WANG等[17]、堯輝明等[18]從輪軌黏滑角度探究了波磨的形成機制,認(rèn)為輪軌微觀接觸界面的不穩(wěn)定黏滑振動促使了波磨的生成和發(fā)展。黏滑理論結(jié)合了微觀接觸和宏觀振動兩個尺度,能夠更為全面且直觀地表征鋼軌材料磨耗和系統(tǒng)振動特性,因而更有利于描述鋼軌波磨的發(fā)生和發(fā)展過程?,F(xiàn)有文獻中關(guān)于鋼軌波磨黏滑振動的研究主要基于半空間理論和穩(wěn)態(tài)假設(shè),且未考慮輪軌瞬態(tài)滾動接觸特性[19],鑒于此,本文采用有限元方法研究鋼軌波磨的形成原因,該方法能夠詳盡地考慮輪軌系統(tǒng)動力學(xué)、接觸力學(xué)與表面損傷之間的相互作用,這對解釋鋼軌波磨現(xiàn)象至關(guān)重要。

    目前地鐵系統(tǒng)中鋼彈簧浮置板軌道使用較為廣泛,且一般用于減振要求較高的特殊地鐵區(qū)段。鋼彈簧浮置板軌道主要通過浮置板受力后的慣性作用實現(xiàn)振動傳遞的衰減,從而減少浮置板傳遞到道床的振動能量。雖然鋼彈簧浮置板軌道的采用有效降低了輪軌系統(tǒng)周圍環(huán)境的振動水平,但它似乎對輪軌界面并不友好,大量的現(xiàn)場調(diào)研發(fā)現(xiàn)鋼彈簧浮置板軌道上出現(xiàn)了嚴(yán)重的鋼軌波磨現(xiàn)象,尤其是在小半徑曲線區(qū)間。為了闡明小半徑曲線鋼彈簧浮置板軌道上的波磨成因,本文采用有限元方法,從輪軌黏滑振動角度描述鋼軌波磨的演化過程。首先,根據(jù)現(xiàn)場線路情況,建立輪對-浮置板軌道三維有限元模型;然后,利用上述模型分析輪軌接觸黏滑特征和鋼軌縱向磨耗特性,從而量化波磨的發(fā)展趨勢;最后,借助復(fù)模態(tài)理論,研究輪軌系統(tǒng)的固有振動屬性與波磨發(fā)展的關(guān)系,以期理解浮置板軌道上的波磨機制。

    1 模型的建立

    1.1 波磨調(diào)研

    實測波磨區(qū)間位于天津某地鐵小半徑曲線線路,曲線半徑為350 m,外軌超高為95 mm,軌距為1435 mm,鋼軌類型為CHN60,軌下結(jié)構(gòu)為鋼彈簧浮置板,其上部通過ZX-2扣件與鋼軌相連,下部則直接與道床連接。線路運營車輛為地鐵B型車,運行速度為60 km/h?,F(xiàn)場調(diào)研發(fā)現(xiàn),小半徑曲線區(qū)間內(nèi)鋼軌波磨較為嚴(yán)重,且主要表現(xiàn)為短波波磨,波長約為25 mm。鋼軌波磨現(xiàn)場照片見圖1。

    (a)整體圖

    1.2 有限元模型

    參照波磨區(qū)間線路情況,利用有限元軟件ABAQUS建立了輪對-浮置板小半徑曲線軌道三維實體模型,如圖2所示,主要包括輪對、鋼軌、ZX-2扣件、浮置板、鋼彈簧隔振器和道床,其中扣件連接鋼軌和浮置板,鋼彈簧隔振器連接浮置板與道床。有限元模型中軌道線型為圓曲線,曲線半徑為350 m,外軌超高參考實際線路設(shè)為95 mm,軌距為1435 mm,軌道長度為25 m,其中鋼軌和道床沿模型全長布置,而浮置板共計4塊,每塊長度為6 m并沿軌道中心線連續(xù)布置。車輪型面為LM磨耗型,鋼軌型面為CHN60,外輪外軌接觸角設(shè)為18.6°,內(nèi)輪內(nèi)軌接觸角設(shè)為1.50°[20]。扣件間距為0.6 m,鋼彈簧隔振器間距為1.2 m,道床底部均勻布置接地彈簧與地基相連,以上連接部件均通過彈簧阻尼單元進行模擬,能夠考慮三個方向上的剛度和阻尼特性。模擬各部件的實體結(jié)構(gòu)均采用C3D8R單元進行離散,輪軌接觸區(qū)域網(wǎng)格尺寸設(shè)置為1 mm,其余輪軌區(qū)域網(wǎng)格尺寸逐漸過渡至30 mm,浮置板和道床結(jié)構(gòu)網(wǎng)格尺寸均勻設(shè)置為100 mm,整個模型共計279.4516萬個單元,325.3375萬個節(jié)點。有限元模型材料及結(jié)構(gòu)參數(shù)取值見表1[21-22],其中材料屬性均考慮為線彈性。

    表1 模型結(jié)構(gòu)參數(shù)取值Tab.1 Values of model structural parameters

    (a)等軸視圖

    有限元模型中輪軌接觸部分采用Coulomb摩擦模型進行表征,法向接觸為硬接觸,切向接觸采用罰函數(shù)法描述,輪軌摩擦因數(shù)設(shè)為0.35[23]。模型中車輪平移速度為60 km/h,轉(zhuǎn)動速度為39.68 rad/s(車輪名義滾動圓半徑為420 mm),且輪對兩端軸箱位置處承受幅值為69 kN的垂向力[24],以模擬一系懸掛力。模型邊界條件設(shè)置如下:軌道結(jié)構(gòu)縱向兩端面固定約束,浮置板和道床板橫向兩端面約束橫向位移,其他表面及結(jié)構(gòu)無任何約束。

    本文主要采用隱式積分方法分析輪軌瞬態(tài)滾動接觸行為,該方法在求解t+Δt時刻的控制方程時需要基于t和t+Δt兩個時刻的動態(tài)參量。因為上述時刻的動態(tài)參量隱含于控制方程,所以需要進行非線性計算。對于一般性的結(jié)構(gòu)問題,隱式積分方法在能夠給出可接受解時,通常需要更多的時間步,即需要更長的計算時間。此外,結(jié)構(gòu)響應(yīng)預(yù)測精度將隨著時間步長Δt相對典型響應(yīng)振型周期T的增大而減小,因此,在選擇最大允許時間步長時,應(yīng)著重考慮如下三個方面:①施加荷載的變化速率;②非線性阻尼和剛度屬性的復(fù)雜性;③結(jié)構(gòu)的典型振動周期。一般情況下,當(dāng)Δt/T<1/10時,隱式積分方法能夠提供可靠的計算結(jié)果[25]。在ABAQUS軟件中,隱式積分方法的時間步能夠在殘差半增量[25]的基礎(chǔ)上進行自動選擇,當(dāng)獲得t+Δt時刻的解后,通過對比平衡殘差在t+Δt/2時刻的值,便可以評估解的精確性并適當(dāng)?shù)卣{(diào)整時間步。

    隱式積分方法的動態(tài)求解過程主要由位移和速度積分的Newmark公式定義:

    u|t+Δt=

    (1)

    (2)

    基于以上計算方法,下面將著重探究輪軌接觸黏滑與鋼軌縱向磨耗特征及輪軌系統(tǒng)固有振動屬性,以期解釋地鐵小半徑曲線鋼彈簧浮置板軌道上的鋼軌波磨形成原因。

    2 接觸黏滑與磨耗特征分析

    2.1 輪軌接觸黏滑特性

    通過運用輪對-浮置板小曲線軌道有限元模型并計算,可以獲得輪對運行過程中各個參量的動態(tài)響應(yīng)。模型中輪對運行距離為6.25 m,即運行時間為0.375(6.25/60×3.6)s,其中包括1 m(5.25~6.25 m)的輪軌接觸求解區(qū)域;計算分析步長設(shè)為0.02 ms。上述輪對運行距離、輪軌接觸求解區(qū)域長度和計算分析步長參數(shù)的選擇主要依據(jù)文獻[26-27]。由于輪軌接觸求解區(qū)域的網(wǎng)格尺寸對模型數(shù)值解的精度具有重要影響,因此,在分析輪軌接觸黏滑特性之前,本節(jié)評估了網(wǎng)格尺寸大小的合理性。對于輪軌接觸求解區(qū)域,分別選擇大小為0.25 mm、0.50 mm和1.00 mm三種網(wǎng)格進行相同工況計算,相應(yīng)的結(jié)果(對應(yīng)接觸求解區(qū)域鋼軌中間位置)匯總于表2。由表2可得,三種網(wǎng)格尺寸的計算結(jié)果誤差較小,均在10%以內(nèi),但計算時間卻隨著網(wǎng)格尺寸減小而劇增,因此,考慮到計算效率,選擇1.00 mm的網(wǎng)格尺寸是合理的。以上計算所使用的處理器為Inter(R) Xeon(R) W-2145 CPU @ 3.70 GHz,機帶RAM為256.0 GB,處理器個數(shù)為8。

    表2 不同網(wǎng)格尺寸下的數(shù)值結(jié)果匯總Tab.2 Summary of numerical results for different mesh sizes

    輪軌接觸狀態(tài)是表征黏滑特性的重要參數(shù),因此,本節(jié)主要從輪軌接觸狀態(tài)變化出發(fā)分析內(nèi)外側(cè)輪軌界面的黏滑特性。輪對在接觸求解區(qū)域穩(wěn)態(tài)運行過程中5個代表性時刻(代表一個完整的接觸狀態(tài)變化過程),即0.350 28 s(時刻1)、0.350 70 s(時刻2)、0.351 12 s(時刻3)、0.351 54 s(時刻4)和0.351 96 s(時刻5)對應(yīng)的內(nèi)外側(cè)輪軌接觸狀態(tài)如圖3、圖4所示。需要補充的是,選擇上述時刻跨度的依據(jù)是其表征了內(nèi)外軌接觸狀態(tài)的最小變化周期。

    (a)時刻1 (b)時刻2 (c)時刻3

    由圖3、圖4可得,在輪對運行過程中,內(nèi)側(cè)輪軌接觸狀態(tài)呈現(xiàn)周期性變化,即輪軌黏滑分布呈現(xiàn)周期性變化(部分黏著-部分黏著-滑移-部分黏著-部分黏著),說明內(nèi)側(cè)輪軌系統(tǒng)發(fā)生了周期性的黏滑振動,進而可能導(dǎo)致內(nèi)軌表面出現(xiàn)周期性的波狀磨耗,最終形成鋼軌波磨;外側(cè)輪軌接觸狀態(tài)總體而言未有明顯變化,表明外側(cè)輪軌系統(tǒng)不會出現(xiàn)黏滑振動,從而不易形成鋼軌波磨。由于鋼軌波磨與微觀材料磨耗密切相關(guān),因此,為進一步量化鋼軌波磨的形成過程,以下將對鋼軌縱向磨耗特征進行分析。

    2.2 鋼軌縱向磨耗特征

    根據(jù)Archard磨耗模型[28]

    (3)

    可以得到磨耗體積Vw與法向載荷FN和滑移距離S成正比。其中,k為磨耗比例系數(shù),其大小與輪軌法向接觸應(yīng)力及相對滑移速度有關(guān);H為兩接觸物體中較軟材料的硬度。法向載荷決定了材料的磨耗程度,而滑移則決定了磨耗是否發(fā)生,即在黏著區(qū),輪軌間沒有相對滑移,故不會產(chǎn)生磨耗,而在滑移區(qū),由于節(jié)點之間存在相對滑移,故會造成材料的磨耗。然而實際上,黏著區(qū)會存在黏著磨損,且黏著磨損的存在也會對鋼軌波磨的形成產(chǎn)生影響,因此,為使計算結(jié)果更為精確,本文在建模過程中指定切向接觸罰函數(shù)法中的最大彈性滑移為表面特征尺寸的0.005%,以考慮黏著磨損的作用。鑒于此,本節(jié)主要根據(jù)鋼軌滑移云圖分析鋼軌縱向磨耗特征,內(nèi)外側(cè)鋼軌的滑移云圖見圖5、圖6。

    (a)鋼軌縱向

    對比圖5a和圖6a可發(fā)現(xiàn),在鋼軌縱向上,內(nèi)側(cè)鋼軌滑移量明顯大于外側(cè)鋼軌滑移量,表明內(nèi)側(cè)鋼軌縱向磨耗程度更大。同時,通過對比圖5b和圖6b可以發(fā)現(xiàn),在鋼軌橫向上,內(nèi)外側(cè)鋼軌滑移量相差不大,表明內(nèi)外側(cè)鋼軌橫向磨耗程度相近。由于鋼軌磨耗的發(fā)生是縱向和橫向滑移共同作用的結(jié)果,因此,通過提取軌面接觸區(qū)域節(jié)點縱橫向滑移量,借助數(shù)學(xué)軟件MATLAB,可繪制縱向和橫向滑移綜合作用下的滑移云圖,如圖7a所示。根據(jù)式(3)計算得到的鋼軌磨耗曲線如圖7b所示。其中滑移距離S包括黏著區(qū)中的彈性滑移(對應(yīng)黏著磨損),k和H的取值參考文獻[29]。由圖7可得,內(nèi)側(cè)鋼軌總滑移量稍大于外側(cè)鋼軌總滑移量,這是內(nèi)外側(cè)鋼軌較大且相近的橫向滑移量的補償作用所致。此外,由圖5~圖7容易看出,內(nèi)外側(cè)鋼軌表面縱橫向滑移云圖分布及磨耗曲線均表現(xiàn)出了一定程度的周期性,尤其是內(nèi)軌縱向滑移云圖,其波長28 mm(圖8)與實測波磨波長25 mm相近,這與內(nèi)側(cè)輪軌界面接觸黏滑狀態(tài)的周期特性相呼應(yīng),從而進一步說明了內(nèi)軌更容易發(fā)生嚴(yán)重的周期性磨耗,即鋼軌波磨。上述分析表明輪軌系統(tǒng)動態(tài)響應(yīng)始終存在周期特性,而周期特性如何形成還需進一步研究。

    (a)鋼軌總滑移云圖

    圖8 內(nèi)軌縱向滑移局部云圖Fig.8 Local nephogram of longitudinal slip of inner rail

    3 輪軌系統(tǒng)振動特性及參數(shù)敏感性分析

    3.1 振動特性分析

    通過運用復(fù)模態(tài)理論[5-6]并進行復(fù)特征值分析,提取了輪軌系統(tǒng)在0~1200 Hz范圍內(nèi)的不穩(wěn)定振型,如圖9、圖10所示。其中穩(wěn)定/不穩(wěn)定振型主要以復(fù)特征值實部的正負(fù)進行反映,如果復(fù)特征值實部大于0,則表明系統(tǒng)振型是不穩(wěn)定的,且數(shù)值越大越不穩(wěn)定,反之則相反。需要說明的是,輪軌接觸斑沿鋼軌縱向長度一般約為20 mm,由于接觸斑對波磨波長具有接觸過濾作用[30],所以在實測區(qū)間車輛運行速度為60 km/h的情況下,根據(jù)頻率計算公式

    圖9 664.1 Hz對應(yīng)的系統(tǒng)不穩(wěn)定振型Fig.9 Unstable mode of the system corresponding to 664.1 Hz

    圖10 665.8 Hz對應(yīng)的系統(tǒng)不穩(wěn)定振型Fig.10 Unstable mode of the system corresponding to 665.8 Hz

    (4)

    可得波磨波長對應(yīng)的最高頻率為833.3 Hz,由此可知復(fù)特征值分析選取的頻率范圍0~1200 Hz是足夠的。式(4)中,ω為波磨通過頻率,Hz;v為車輛運行速度,km/h;λ為波磨波長,mm。

    由圖9、圖10可知,在0~1200 Hz范圍內(nèi),輪軌系統(tǒng)共出現(xiàn)了兩個不穩(wěn)定振型,分別對應(yīng)頻率為664.1 Hz和665.8 Hz,振型均表現(xiàn)為內(nèi)輪相對于軌道的彎曲振動,且兩個不穩(wěn)定振型對應(yīng)的復(fù)特征值實部分別為15.9和16.2,這說明內(nèi)側(cè)輪軌系統(tǒng)在此頻率下發(fā)生了失穩(wěn)。同時,由1.1節(jié)可得實測波磨區(qū)間車輛運行速度為60 km/h,內(nèi)軌波磨波長約為25 mm,因此,由式(4)可得內(nèi)軌波磨通過頻率約為666.7 Hz。因為復(fù)特征值分析所得不穩(wěn)定振型對應(yīng)頻率與內(nèi)軌波磨通過頻率均相近,所以可以得出系統(tǒng)響應(yīng)周期特性的形成機制為輪軌系統(tǒng)664.1 Hz和665.8 Hz對應(yīng)的不穩(wěn)定振型所致,且該振型均表現(xiàn)為內(nèi)輪相對于軌道的彎曲振動。

    綜合輪軌接觸黏滑與磨耗特征及輪軌系統(tǒng)振動特性分析結(jié)果,可以得出實測地鐵小半徑曲線浮置板軌道上的內(nèi)軌波磨是由664.1 Hz和665.8 Hz對應(yīng)的輪軌系統(tǒng)不穩(wěn)定振型誘發(fā)的黏滑振動造成的,且該不穩(wěn)定振型均表現(xiàn)為內(nèi)輪相對于軌道的彎曲振動。

    3.2 參數(shù)敏感性分析

    由于鋼軌波磨是輪軌系統(tǒng)參數(shù)的一種集中體現(xiàn),因此為確定不穩(wěn)定黏滑振動的產(chǎn)生條件以實現(xiàn)鋼軌波磨的有效控制,本節(jié)分析了系統(tǒng)參數(shù)對不穩(wěn)定振動的影響。選取扣件和鋼彈簧三個方向上的剛度和阻尼特性作為代表性輸入?yún)?shù),各參數(shù)分布類型均為均勻分布,均值與表1中數(shù)值一致,最小/最大極限值如表3所示[31-33],其中Kfv、Kft、Kfl分別表示扣件垂向、橫向、縱向剛度;Dfv、Dft、Dfl分別表示扣件垂向、橫向、縱向阻尼;Ksv、Kst、Ksl分別表示鋼彈簧垂向、橫向、縱向剛度;Dsv、Dst、Dsl分別表示鋼彈簧垂向、橫向、縱向阻尼。基于3.1節(jié)中的數(shù)值模型并以復(fù)特征值實部作為分析目標(biāo),利用有限元隨機分析模塊并執(zhí)行概率設(shè)計分析,可以獲得反映鋼軌波磨發(fā)生趨勢的不穩(wěn)定振動的評價指標(biāo)-復(fù)特征值實部的統(tǒng)計特性,以作為參數(shù)敏感性分析的基礎(chǔ)。

    表3 輸入?yún)?shù)最小/最大極限值Tab.3 Minimum/maximum limit values of input parameters

    采用Spearman秩相關(guān)系數(shù)rsj[34]間接表示參數(shù)的敏感性,其計算公式如下:

    (5)

    j=1,2,…,12 -1≤rsj≤1

    式中,Pkj為第j個隨機參數(shù)的秩;Qkj為與Pkj對應(yīng)的響應(yīng)量的秩;n為總樣本數(shù),本節(jié)取為800。

    Spearman秩相關(guān)系數(shù)可用于度量變量之間聯(lián)系的強弱,如果計算所得系數(shù)接近1或-1,則認(rèn)為輸入變量對輸出變量影響顯著;如果計算所得系數(shù)趨近0,則認(rèn)為影響微弱。如果系數(shù)為正值,表明增大輸入變量的值,輸出變量的值也增大;如果系數(shù)為負(fù)值,表明增大輸入變量的值,輸出變量的值反而減小。

    將復(fù)特征值實部對各系統(tǒng)參數(shù)的敏感度匯總于表4,其中ξi表示復(fù)特征值實部。由表4可得,對復(fù)特征值實部影響較大的前兩個參數(shù)分別為扣件垂向剛度和鋼彈簧垂向剛度,且兩者對應(yīng)的敏感度為負(fù),表明它與復(fù)特征值實部成負(fù)相關(guān),這也可從圖11、圖12所示的復(fù)特征值實部隨扣件垂向剛度和鋼彈簧垂向剛度變化的散點圖中看出。由于減小復(fù)特征值實部可以有效抑制系統(tǒng)的不穩(wěn)定振動,從而控制鋼軌波磨的發(fā)生,因此,適當(dāng)?shù)卦龃罂奂弯搹椈傻拇瓜騽偠饶軌驅(qū)︿撥壊サ目刂破鸬椒e極作用,這也從系統(tǒng)參數(shù)角度說明鋼彈簧較低的垂向剛度特性誘發(fā)了實測線路上鋼軌波磨的形成(實測線路上ZX-2扣件為普通減振扣件,即非低剛度扣件)。其他參數(shù)對復(fù)特征值實部影響較小,說明這些參數(shù)對實測線路上的波磨形成影響不大。

    表4 復(fù)特征值實部ξi的敏感度Tab.4 Sensitivities of real part ξi of complex eigenvalues

    圖11 ξi隨Kfv變化散點圖Fig.11 Scatter diagrams of ξi varying with Kfv

    圖12 ξi隨Ksv變化散點圖Fig.12 Scatter diagrams of ξi varying with Ksv

    4 結(jié)論

    為解釋地鐵小半徑曲線鋼彈簧浮置板軌道上的波磨成因,本文從輪軌黏滑振動角度詳細地闡述了鋼軌波磨的形成過程。首先,依據(jù)現(xiàn)場線路情況完成了輪對-浮置板軌道三維有限元模型的構(gòu)建;然后利用上述模型分析了輪軌接觸黏滑特征和鋼軌縱向磨耗特性,從而量化了波磨的發(fā)生趨勢;最后,基于復(fù)模態(tài)理論研究了輪軌系統(tǒng)的固有振動屬性與波磨形成的關(guān)系。主要得到以下結(jié)論:

    (1)在輪對運行階段,內(nèi)側(cè)輪軌黏滑分布呈現(xiàn)周期性變化,說明內(nèi)側(cè)輪軌系統(tǒng)發(fā)生了周期性的黏滑振動,進而可能導(dǎo)致內(nèi)軌表面出現(xiàn)周期性的波狀磨耗,最終形成鋼軌波磨;外側(cè)輪軌黏滑分布總體上未有明顯變化,說明外側(cè)輪軌系統(tǒng)未出現(xiàn)黏滑振動,從而不易形成鋼軌波磨。

    (2)在鋼軌縱向上,內(nèi)軌表面滑移量明顯大于外軌表面滑移量,表明內(nèi)軌縱向磨耗程度更大;在鋼軌橫向上,內(nèi)外側(cè)鋼軌表面滑移量相差不大,表明內(nèi)外側(cè)鋼軌橫向磨耗程度相近。

    (3)內(nèi)外側(cè)鋼軌表面縱橫向滑移云圖分布均表現(xiàn)出了一定程度的周期性,尤其是內(nèi)軌縱向滑移云圖,其波磨波長為28 mm,與實測波磨波長25 mm相近,這與內(nèi)側(cè)輪軌界面接觸黏滑狀態(tài)的周期特性相呼應(yīng),從而進一步說明了內(nèi)軌更容易發(fā)生嚴(yán)重的周期性磨耗,即鋼軌波磨。

    (4)復(fù)模態(tài)分析顯示,內(nèi)側(cè)輪軌系統(tǒng)在與實測波磨通過頻率666.7 Hz相近頻率處發(fā)生了失穩(wěn),結(jié)合輪軌接觸黏滑與磨耗特征分析結(jié)果,可以得出實測地鐵小半徑曲線浮置板軌道上的內(nèi)軌波磨是由664.1 Hz和665.8 Hz對應(yīng)的輪軌系統(tǒng)不穩(wěn)定振型誘發(fā)的黏滑振動造成的,且該不穩(wěn)定振型均表現(xiàn)為內(nèi)輪相對于軌道的彎曲振動。

    (5)參數(shù)敏感性分析顯示扣件和鋼彈簧垂向剛度對復(fù)特征值實部影響較大,且與復(fù)特征值實部負(fù)相關(guān)。由于復(fù)特征值實部的減小能夠有效地抑制系統(tǒng)的不穩(wěn)定振動,從而控制波磨的發(fā)生,因此,適當(dāng)?shù)卦龃罂奂弯搹椈傻拇瓜騽偠饶軌驅(qū)Σタ刂破鸬椒e極作用。

    猜你喜歡
    波磨浮置輪軌
    城軌小半徑曲線鋼軌波磨的成因與整治研究綜述
    北京地鐵鋼軌波磨現(xiàn)狀及防治措施
    中低速磁浮道岔與輪軌道岔的差異
    移動荷載作用下離散支承浮置板軌道振動響應(yīng)研究
    中低速磁浮與輪軌交通信號系統(tǒng)的差異
    預(yù)制浮置板軌道隔振性能研究
    非線性穩(wěn)態(tài)曲線通過時輪軌滾動接觸的數(shù)值求解方法
    高速鐵路鋼軌波磨對車輛—軌道動態(tài)響應(yīng)的影響
    基于多島遺傳算法的橡膠浮置板彈性的優(yōu)化
    減緩高速鐵路鋼軌波磨的仿真分析
    www.999成人在线观看| 色综合站精品国产| 亚洲久久久国产精品| 国产一区二区激情短视频| 成人亚洲精品一区在线观看| 国产高清有码在线观看视频 | 日韩三级视频一区二区三区| 久久中文字幕人妻熟女| 一个人观看的视频www高清免费观看 | 黄色 视频免费看| 每晚都被弄得嗷嗷叫到高潮| 国产不卡一卡二| 宅男免费午夜| 俄罗斯特黄特色一大片| 亚洲av电影在线进入| 免费搜索国产男女视频| 日韩有码中文字幕| 国产av一区在线观看免费| 日日爽夜夜爽网站| 国产精品免费视频内射| 久久人人精品亚洲av| 一区二区日韩欧美中文字幕| 欧美乱色亚洲激情| 国产亚洲精品综合一区在线观看 | 色播在线永久视频| 人人妻人人澡欧美一区二区 | 久久久久国产精品人妻aⅴ院| 女人高潮潮喷娇喘18禁视频| 咕卡用的链子| 亚洲全国av大片| 午夜福利,免费看| 国产精品爽爽va在线观看网站 | 首页视频小说图片口味搜索| www.精华液| 亚洲 欧美一区二区三区| 欧美不卡视频在线免费观看 | 日本在线视频免费播放| 丝袜美足系列| 在线观看www视频免费| 精品国产一区二区三区四区第35| 一二三四社区在线视频社区8| 免费人成视频x8x8入口观看| 欧美人与性动交α欧美精品济南到| 狂野欧美激情性xxxx| 校园春色视频在线观看| 可以免费在线观看a视频的电影网站| 首页视频小说图片口味搜索| www.熟女人妻精品国产| 日日爽夜夜爽网站| 亚洲一区中文字幕在线| 欧美成人午夜精品| 国产欧美日韩一区二区三区在线| 91成人精品电影| 精品国产美女av久久久久小说| 99国产精品免费福利视频| 一个人免费在线观看的高清视频| 亚洲国产精品成人综合色| 欧美色视频一区免费| 狂野欧美激情性xxxx| 好看av亚洲va欧美ⅴa在| 亚洲av成人不卡在线观看播放网| 精品国产乱码久久久久久男人| aaaaa片日本免费| 在线国产一区二区在线| av中文乱码字幕在线| 一个人免费在线观看的高清视频| 亚洲成国产人片在线观看| 久久精品亚洲熟妇少妇任你| 日韩欧美一区二区三区在线观看| 成年人黄色毛片网站| 黄片播放在线免费| www.自偷自拍.com| 一夜夜www| 亚洲第一欧美日韩一区二区三区| 一边摸一边抽搐一进一出视频| 免费少妇av软件| 最近最新中文字幕大全免费视频| 精品一区二区三区四区五区乱码| 欧美日本视频| 亚洲欧美精品综合一区二区三区| 美女大奶头视频| 黄频高清免费视频| 久久久久国产一级毛片高清牌| 日韩国内少妇激情av| 国产成人系列免费观看| 麻豆av在线久日| 天堂影院成人在线观看| 999久久久国产精品视频| 国产亚洲欧美在线一区二区| 老司机午夜十八禁免费视频| 一边摸一边抽搐一进一出视频| 男女做爰动态图高潮gif福利片 | 咕卡用的链子| 999久久久国产精品视频| 国产欧美日韩一区二区三| 一个人免费在线观看的高清视频| 如日韩欧美国产精品一区二区三区| 精品一区二区三区av网在线观看| 老汉色∧v一级毛片| 国产黄a三级三级三级人| 免费在线观看黄色视频的| 久久国产亚洲av麻豆专区| 精品人妻1区二区| 俄罗斯特黄特色一大片| 亚洲欧洲精品一区二区精品久久久| 亚洲欧美激情在线| 欧美黑人欧美精品刺激| aaaaa片日本免费| 老司机深夜福利视频在线观看| 亚洲第一av免费看| 午夜免费激情av| 91老司机精品| av中文乱码字幕在线| а√天堂www在线а√下载| 日本在线视频免费播放| 午夜精品国产一区二区电影| 搡老妇女老女人老熟妇| 国产aⅴ精品一区二区三区波| 久久婷婷成人综合色麻豆| 国产一区二区激情短视频| 免费不卡黄色视频| 亚洲专区国产一区二区| 国产亚洲欧美在线一区二区| 亚洲国产中文字幕在线视频| 午夜两性在线视频| 久久久精品欧美日韩精品| 性欧美人与动物交配| 又黄又爽又免费观看的视频| 99国产精品免费福利视频| 午夜老司机福利片| 成年版毛片免费区| 国产精品一区二区三区四区久久 | 成熟少妇高潮喷水视频| 日本黄色视频三级网站网址| 成在线人永久免费视频| 级片在线观看| 亚洲全国av大片| 成人亚洲精品一区在线观看| www.熟女人妻精品国产| 亚洲少妇的诱惑av| 免费搜索国产男女视频| 亚洲中文日韩欧美视频| 久热这里只有精品99| 久久久久国内视频| 黄色成人免费大全| 国产一区二区三区综合在线观看| 精品久久蜜臀av无| 亚洲欧美激情在线| 性色av乱码一区二区三区2| 国产激情久久老熟女| 亚洲三区欧美一区| 精品免费久久久久久久清纯| 50天的宝宝边吃奶边哭怎么回事| 国产av在哪里看| 国产成人精品久久二区二区91| 亚洲精品国产色婷婷电影| 不卡一级毛片| cao死你这个sao货| 久久国产亚洲av麻豆专区| 美女大奶头视频| 久久久久九九精品影院| 久久久久久久久中文| 97人妻精品一区二区三区麻豆 | 欧美日韩福利视频一区二区| 欧美不卡视频在线免费观看 | 免费av毛片视频| 麻豆国产av国片精品| 亚洲国产日韩欧美精品在线观看 | 老鸭窝网址在线观看| 久久久久久久久免费视频了| 亚洲国产看品久久| 国产极品粉嫩免费观看在线| 午夜激情av网站| 变态另类丝袜制服| 国产成人啪精品午夜网站| 999久久久精品免费观看国产| 久久久国产精品麻豆| 亚洲avbb在线观看| www.自偷自拍.com| 一本久久中文字幕| 在线观看舔阴道视频| 日本 av在线| 波多野结衣一区麻豆| 欧美日本视频| 99国产精品免费福利视频| 12—13女人毛片做爰片一| 老汉色av国产亚洲站长工具| 校园春色视频在线观看| 午夜免费成人在线视频| 青草久久国产| 亚洲一区中文字幕在线| 亚洲人成77777在线视频| 欧美日韩亚洲国产一区二区在线观看| 好男人电影高清在线观看| 欧美成人性av电影在线观看| 国产成人免费无遮挡视频| 免费少妇av软件| 后天国语完整版免费观看| 搡老熟女国产l中国老女人| 国产av一区二区精品久久| 久久香蕉激情| 久久香蕉精品热| 91成年电影在线观看| 久久伊人香网站| 午夜精品国产一区二区电影| 高清毛片免费观看视频网站| 一区二区三区国产精品乱码| 少妇粗大呻吟视频| 国产欧美日韩一区二区三| 国产精品影院久久| 国产一区二区三区综合在线观看| a在线观看视频网站| 法律面前人人平等表现在哪些方面| 亚洲一区二区三区不卡视频| 亚洲 国产 在线| 美女高潮到喷水免费观看| 成年人黄色毛片网站| 黑人欧美特级aaaaaa片| 黑人巨大精品欧美一区二区mp4| 两个人免费观看高清视频| 亚洲aⅴ乱码一区二区在线播放 | 国产成年人精品一区二区| 亚洲五月婷婷丁香| 精品一区二区三区av网在线观看| 久久久久亚洲av毛片大全| 电影成人av| 色播在线永久视频| 国产欧美日韩一区二区三| 桃色一区二区三区在线观看| 国产高清视频在线播放一区| 禁无遮挡网站| 亚洲专区国产一区二区| 欧美成人性av电影在线观看| 国产熟女xx| 欧美激情 高清一区二区三区| 97超级碰碰碰精品色视频在线观看| 亚洲熟女毛片儿| 香蕉丝袜av| 91老司机精品| 男女下面进入的视频免费午夜 | 免费在线观看视频国产中文字幕亚洲| 在线免费观看的www视频| 欧美丝袜亚洲另类 | 两性夫妻黄色片| 超碰成人久久| 露出奶头的视频| 欧美国产日韩亚洲一区| 97碰自拍视频| 成人手机av| 午夜福利影视在线免费观看| 国产一区二区在线av高清观看| 国产亚洲欧美精品永久| 日本五十路高清| 桃红色精品国产亚洲av| 免费无遮挡裸体视频| 成人免费观看视频高清| 欧美最黄视频在线播放免费| 国产私拍福利视频在线观看| 欧美另类亚洲清纯唯美| 亚洲美女黄片视频| 日本vs欧美在线观看视频| 在线观看www视频免费| 亚洲aⅴ乱码一区二区在线播放 | 一区二区三区国产精品乱码| 亚洲精品美女久久av网站| 熟女少妇亚洲综合色aaa.| 国产精品影院久久| 国产麻豆成人av免费视频| 亚洲成人免费电影在线观看| 久久天躁狠狠躁夜夜2o2o| 国产亚洲精品综合一区在线观看 | 日韩欧美一区二区三区在线观看| 欧美日韩亚洲国产一区二区在线观看| 亚洲国产精品sss在线观看| 麻豆国产av国片精品| 久热爱精品视频在线9| 亚洲九九香蕉| 中文字幕高清在线视频| 老汉色av国产亚洲站长工具| 亚洲在线自拍视频| 午夜福利,免费看| 国产激情久久老熟女| 涩涩av久久男人的天堂| 无人区码免费观看不卡| 亚洲成人久久性| av超薄肉色丝袜交足视频| 欧美乱码精品一区二区三区| 天堂动漫精品| 最新美女视频免费是黄的| 午夜福利欧美成人| 国产单亲对白刺激| 在线视频色国产色| e午夜精品久久久久久久| 老司机靠b影院| 午夜福利在线观看吧| 无限看片的www在线观看| 色综合站精品国产| 91精品国产国语对白视频| 国产成年人精品一区二区| 丰满人妻熟妇乱又伦精品不卡| 午夜福利18| 国产精品秋霞免费鲁丝片| 精品一区二区三区视频在线观看免费| 又大又爽又粗| 97人妻天天添夜夜摸| 免费看美女性在线毛片视频| 咕卡用的链子| 久久中文字幕人妻熟女| 国产亚洲精品一区二区www| 日韩三级视频一区二区三区| 亚洲色图综合在线观看| 亚洲久久久国产精品| 天堂影院成人在线观看| 亚洲精品久久成人aⅴ小说| 成年女人毛片免费观看观看9| 嫩草影视91久久| 久久香蕉激情| 国产精品综合久久久久久久免费 | 国产精品一区二区精品视频观看| 亚洲av熟女| 搡老熟女国产l中国老女人| 女性生殖器流出的白浆| 日日爽夜夜爽网站| 色综合站精品国产| 色综合婷婷激情| 国产精品九九99| 久久久久久久久免费视频了| 99国产精品免费福利视频| 亚洲天堂国产精品一区在线| 免费av毛片视频| 欧洲精品卡2卡3卡4卡5卡区| av欧美777| 久久精品国产99精品国产亚洲性色 | 欧美黄色片欧美黄色片| 最近最新中文字幕大全电影3 | 久久久水蜜桃国产精品网| 国产日韩一区二区三区精品不卡| 国产一卡二卡三卡精品| 久久午夜综合久久蜜桃| 国产一卡二卡三卡精品| 久久人人爽av亚洲精品天堂| 欧美激情 高清一区二区三区| 在线播放国产精品三级| 国产精品久久久久久精品电影 | 精品少妇一区二区三区视频日本电影| 一个人观看的视频www高清免费观看 | 精品国产一区二区三区四区第35| 欧美另类亚洲清纯唯美| 久热这里只有精品99| 亚洲avbb在线观看| 国产免费av片在线观看野外av| 亚洲男人的天堂狠狠| 国产又色又爽无遮挡免费看| 欧美成人免费av一区二区三区| 亚洲国产精品sss在线观看| 十分钟在线观看高清视频www| 韩国精品一区二区三区| 精品一区二区三区四区五区乱码| www.自偷自拍.com| 两人在一起打扑克的视频| 成人国产一区最新在线观看| 精品国产一区二区久久| 99国产极品粉嫩在线观看| 18禁黄网站禁片午夜丰满| 日本免费一区二区三区高清不卡 | 啦啦啦免费观看视频1| 首页视频小说图片口味搜索| 美女午夜性视频免费| 久久精品国产99精品国产亚洲性色 | av在线天堂中文字幕| 中文字幕高清在线视频| 久久中文字幕一级| 超碰成人久久| 电影成人av| 一二三四在线观看免费中文在| 18禁美女被吸乳视频| 亚洲最大成人中文| a级毛片在线看网站| 757午夜福利合集在线观看| 国产一区二区三区综合在线观看| 久久国产精品人妻蜜桃| 国产精品二区激情视频| 丰满人妻熟妇乱又伦精品不卡| 乱人伦中国视频| 久久这里只有精品19| 成人国产一区最新在线观看| 女人精品久久久久毛片| 久久狼人影院| 久久午夜综合久久蜜桃| 亚洲熟妇熟女久久| 久久精品国产综合久久久| 一个人免费在线观看的高清视频| 欧美乱色亚洲激情| 啪啪无遮挡十八禁网站| 亚洲欧美日韩另类电影网站| 欧美成狂野欧美在线观看| 国产一级毛片七仙女欲春2 | 男男h啪啪无遮挡| 成人免费观看视频高清| 国产精品亚洲美女久久久| 性色av乱码一区二区三区2| 女性被躁到高潮视频| 一个人观看的视频www高清免费观看 | 婷婷六月久久综合丁香| 欧美久久黑人一区二区| 亚洲中文av在线| 一级毛片高清免费大全| 国产精品av久久久久免费| 久久精品亚洲精品国产色婷小说| 最近最新中文字幕大全免费视频| 97人妻天天添夜夜摸| 国产成人影院久久av| 亚洲精品中文字幕在线视频| 非洲黑人性xxxx精品又粗又长| 最好的美女福利视频网| 国产精品久久电影中文字幕| 亚洲,欧美精品.| 国产成人av激情在线播放| 老汉色av国产亚洲站长工具| 咕卡用的链子| 免费搜索国产男女视频| 亚洲精品在线美女| 国产欧美日韩综合在线一区二区| 国产av在哪里看| 大码成人一级视频| 午夜福利一区二区在线看| 黑人欧美特级aaaaaa片| 久久久国产精品麻豆| 久久精品亚洲精品国产色婷小说| 电影成人av| 亚洲 欧美 日韩 在线 免费| 在线观看免费午夜福利视频| 男女午夜视频在线观看| 人妻久久中文字幕网| 两个人免费观看高清视频| 国内久久婷婷六月综合欲色啪| 亚洲国产毛片av蜜桃av| 日本三级黄在线观看| 精品久久久久久久人妻蜜臀av | 首页视频小说图片口味搜索| 亚洲精品在线观看二区| 精品欧美一区二区三区在线| 国产免费av片在线观看野外av| 黄片播放在线免费| 久久久久久久午夜电影| 色婷婷久久久亚洲欧美| 大香蕉久久成人网| 男人舔女人的私密视频| 国产单亲对白刺激| 91麻豆精品激情在线观看国产| 成人18禁高潮啪啪吃奶动态图| 91精品三级在线观看| 国产成人免费无遮挡视频| 免费观看人在逋| 手机成人av网站| 18禁裸乳无遮挡免费网站照片 | 国产男靠女视频免费网站| 婷婷丁香在线五月| 午夜免费观看网址| 99国产精品99久久久久| 男女下面插进去视频免费观看| 成年版毛片免费区| 欧美色视频一区免费| 日本五十路高清| 熟女少妇亚洲综合色aaa.| 亚洲av美国av| 欧美乱码精品一区二区三区| 色播亚洲综合网| 国产一级毛片七仙女欲春2 | 亚洲一区二区三区不卡视频| 久久天躁狠狠躁夜夜2o2o| 悠悠久久av| 欧美色欧美亚洲另类二区 | 50天的宝宝边吃奶边哭怎么回事| 大陆偷拍与自拍| 精品福利观看| 国产欧美日韩一区二区精品| 精品久久久久久,| 日韩欧美在线二视频| 十八禁网站免费在线| 好看av亚洲va欧美ⅴa在| av超薄肉色丝袜交足视频| 悠悠久久av| 女人被狂操c到高潮| 一进一出好大好爽视频| 久久 成人 亚洲| 两个人视频免费观看高清| 日本黄色视频三级网站网址| 免费不卡黄色视频| 他把我摸到了高潮在线观看| 一级毛片精品| 亚洲七黄色美女视频| 国产亚洲av高清不卡| 国产成人一区二区三区免费视频网站| 日韩精品免费视频一区二区三区| 我的亚洲天堂| 一夜夜www| 一区二区三区激情视频| 免费在线观看亚洲国产| 国产av一区二区精品久久| 丁香六月欧美| 最近最新中文字幕大全免费视频| 亚洲精品中文字幕在线视频| 中文字幕色久视频| 男女床上黄色一级片免费看| 又黄又粗又硬又大视频| 成人手机av| av天堂在线播放| www.熟女人妻精品国产| 国产精品爽爽va在线观看网站 | 国产欧美日韩一区二区三区在线| 人妻久久中文字幕网| 亚洲精品美女久久久久99蜜臀| 国产伦一二天堂av在线观看| 久久久久九九精品影院| 涩涩av久久男人的天堂| 久久国产精品人妻蜜桃| 久久久精品欧美日韩精品| 在线观看66精品国产| 三级毛片av免费| 男人舔女人下体高潮全视频| 啦啦啦观看免费观看视频高清 | 少妇 在线观看| 日韩有码中文字幕| 国产主播在线观看一区二区| 精品久久久久久久人妻蜜臀av | 日韩国内少妇激情av| 国产精华一区二区三区| 男女午夜视频在线观看| 两个人看的免费小视频| 日本免费a在线| 亚洲三区欧美一区| 国产精品亚洲av一区麻豆| 亚洲男人天堂网一区| 欧美乱色亚洲激情| 精品欧美国产一区二区三| 美女免费视频网站| 午夜免费鲁丝| 男女之事视频高清在线观看| 老司机靠b影院| 老司机午夜福利在线观看视频| 国产精品一区二区三区四区久久 | 超碰成人久久| 夜夜爽天天搞| 亚洲五月色婷婷综合| 男人的好看免费观看在线视频 | 人人妻人人澡人人看| 亚洲电影在线观看av| 久9热在线精品视频| 日韩精品青青久久久久久| av片东京热男人的天堂| 一a级毛片在线观看| 在线播放国产精品三级| 久久精品影院6| 午夜福利影视在线免费观看| 国产精品久久久久久人妻精品电影| 黑人巨大精品欧美一区二区mp4| 99在线视频只有这里精品首页| 好男人电影高清在线观看| 国产单亲对白刺激| 欧美午夜高清在线| 可以在线观看毛片的网站| 视频区欧美日本亚洲| 亚洲最大成人中文| 日韩大码丰满熟妇| 老司机靠b影院| e午夜精品久久久久久久| 国产不卡一卡二| 国产av一区在线观看免费| 久久精品国产清高在天天线| 黑人欧美特级aaaaaa片| 十八禁网站免费在线| 亚洲专区字幕在线| 午夜精品在线福利| 久久国产精品男人的天堂亚洲| 免费观看精品视频网站| 老鸭窝网址在线观看| 日本五十路高清| 91老司机精品| 日韩av在线大香蕉| 人人妻,人人澡人人爽秒播| 9色porny在线观看| 男女下面插进去视频免费观看| 天天躁夜夜躁狠狠躁躁| 亚洲第一电影网av| 黄色视频,在线免费观看| 国产一区二区三区视频了| 美女扒开内裤让男人捅视频| 国产一区二区三区在线臀色熟女| 动漫黄色视频在线观看| 一卡2卡三卡四卡精品乱码亚洲| 久久久国产精品麻豆| 日本黄色视频三级网站网址| 精品日产1卡2卡| 琪琪午夜伦伦电影理论片6080| 黄色成人免费大全| 性欧美人与动物交配| 国产一区二区激情短视频| 国产视频一区二区在线看| 精品久久蜜臀av无| 久久久国产精品麻豆| 动漫黄色视频在线观看| 校园春色视频在线观看| 成在线人永久免费视频| 麻豆久久精品国产亚洲av| 午夜久久久在线观看| 亚洲狠狠婷婷综合久久图片| 亚洲精品中文字幕在线视频| svipshipincom国产片| 俄罗斯特黄特色一大片| 成人av一区二区三区在线看| 亚洲国产精品久久男人天堂| 91av网站免费观看| 亚洲 国产 在线| 黄片播放在线免费|