• <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)化
    減緩高速鐵路鋼軌波磨的仿真分析
    亚洲欧美精品综合久久99| 国产精品亚洲一级av第二区| 精华霜和精华液先用哪个| 国产一区二区激情短视频| 成人无遮挡网站| 久久草成人影院| 国产在线男女| 国产高清视频在线播放一区| 久久人人爽人人爽人人片va | 精品99又大又爽又粗少妇毛片 | 97碰自拍视频| 欧美性猛交黑人性爽| 国产精华一区二区三区| 天堂影院成人在线观看| 日本成人三级电影网站| 国内毛片毛片毛片毛片毛片| 久久久国产成人免费| 中文字幕免费在线视频6| 国产精品不卡视频一区二区 | 别揉我奶头~嗯~啊~动态视频| 日本三级黄在线观看| 中文资源天堂在线| 神马国产精品三级电影在线观看| 国产高清视频在线观看网站| 性色avwww在线观看| 麻豆成人午夜福利视频| 久久久精品大字幕| 国产精品一区二区性色av| 婷婷六月久久综合丁香| 欧美成人a在线观看| 极品教师在线视频| 99国产精品一区二区三区| 人人妻人人澡欧美一区二区| 国产伦人伦偷精品视频| 免费看日本二区| 国产伦一二天堂av在线观看| 日韩成人在线观看一区二区三区| а√天堂www在线а√下载| 日本黄色片子视频| 欧美日本视频| 村上凉子中文字幕在线| 性色av乱码一区二区三区2| 免费在线观看成人毛片| 国产伦在线观看视频一区| 好男人电影高清在线观看| 床上黄色一级片| 亚洲最大成人手机在线| 亚洲欧美日韩高清在线视频| 欧美日韩黄片免| 亚洲激情在线av| 91久久精品国产一区二区成人| 欧美3d第一页| 国产探花极品一区二区| 国产真实伦视频高清在线观看 | 国产人妻一区二区三区在| 又黄又爽又刺激的免费视频.| 又紧又爽又黄一区二区| avwww免费| 69人妻影院| 亚洲第一欧美日韩一区二区三区| 黄色一级大片看看| 精品国产三级普通话版| 性色avwww在线观看| 一本精品99久久精品77| 国内精品美女久久久久久| 国产高清视频在线观看网站| 国产成人av教育| 亚洲人成网站高清观看| 欧美日韩黄片免| 亚州av有码| 午夜福利在线观看吧| 欧美成人免费av一区二区三区| 欧美激情国产日韩精品一区| 亚洲精品成人久久久久久| 国产黄色小视频在线观看| 国产v大片淫在线免费观看| 两性午夜刺激爽爽歪歪视频在线观看| 欧美在线黄色| 亚洲av第一区精品v没综合| 美女高潮的动态| 99国产精品一区二区三区| av国产免费在线观看| 久久6这里有精品| 老司机福利观看| 国产精品99久久久久久久久| 亚洲成人精品中文字幕电影| av在线天堂中文字幕| 亚洲成a人片在线一区二区| 中文字幕高清在线视频| 人人妻人人澡欧美一区二区| 毛片一级片免费看久久久久 | 国产色婷婷99| 国产精品影院久久| 美女xxoo啪啪120秒动态图 | 美女黄网站色视频| 日韩有码中文字幕| 久久精品国产亚洲av涩爱 | 国产 一区 欧美 日韩| 动漫黄色视频在线观看| 精品久久久久久久久av| 亚洲成a人片在线一区二区| 久久精品国产亚洲av涩爱 | 日韩大尺度精品在线看网址| 神马国产精品三级电影在线观看| 欧美在线一区亚洲| 99久久成人亚洲精品观看| 久久香蕉精品热| 欧美黑人欧美精品刺激| 非洲黑人性xxxx精品又粗又长| 赤兔流量卡办理| 午夜福利在线在线| 色精品久久人妻99蜜桃| 天堂影院成人在线观看| 色哟哟哟哟哟哟| 国产黄a三级三级三级人| 久99久视频精品免费| 9191精品国产免费久久| 男女做爰动态图高潮gif福利片| 99riav亚洲国产免费| www.熟女人妻精品国产| 国产精品美女特级片免费视频播放器| 国产午夜精品论理片| 深夜a级毛片| 亚洲av熟女| 亚洲avbb在线观看| 悠悠久久av| 成人特级av手机在线观看| 欧美丝袜亚洲另类 | 神马国产精品三级电影在线观看| 亚洲aⅴ乱码一区二区在线播放| 日韩欧美在线二视频| 亚洲人成网站在线播放欧美日韩| 国产国拍精品亚洲av在线观看| 一级黄片播放器| 国产精品乱码一区二三区的特点| 久久亚洲精品不卡| 欧美xxxx性猛交bbbb| 91在线观看av| 国产免费av片在线观看野外av| 在现免费观看毛片| 久久久久久久久久成人| 亚洲国产精品999在线| 美女xxoo啪啪120秒动态图 | 亚洲自偷自拍三级| 97碰自拍视频| 看十八女毛片水多多多| 国产精品一区二区三区四区久久| 少妇人妻精品综合一区二区 | 国产精品一区二区三区四区久久| 18美女黄网站色大片免费观看| 国产亚洲av嫩草精品影院| 亚洲电影在线观看av| 动漫黄色视频在线观看| 不卡一级毛片| 国产爱豆传媒在线观看| 深夜a级毛片| 亚洲综合色惰| 美女高潮喷水抽搐中文字幕| www.色视频.com| 99久久精品热视频| 一个人免费在线观看的高清视频| 亚洲内射少妇av| а√天堂www在线а√下载| 日本免费a在线| 给我免费播放毛片高清在线观看| 成人鲁丝片一二三区免费| 国产熟女xx| 3wmmmm亚洲av在线观看| 亚洲五月天丁香| 亚洲在线观看片| 亚洲真实伦在线观看| 亚洲七黄色美女视频| 国产综合懂色| 此物有八面人人有两片| 中国美女看黄片| 亚洲三级黄色毛片| 97人妻精品一区二区三区麻豆| 最近最新免费中文字幕在线| 国产亚洲精品av在线| 热99在线观看视频| 非洲黑人性xxxx精品又粗又长| 日韩精品中文字幕看吧| 欧美性猛交黑人性爽| 欧美极品一区二区三区四区| 大型黄色视频在线免费观看| www.www免费av| 成人欧美大片| 无遮挡黄片免费观看| 久久久久国内视频| 99国产精品一区二区三区| 91麻豆精品激情在线观看国产| 三级毛片av免费| 欧美不卡视频在线免费观看| 久久久精品欧美日韩精品| 成人亚洲精品av一区二区| 欧美日韩黄片免| 在线观看美女被高潮喷水网站 | 久久精品久久久久久噜噜老黄 | 国产久久久一区二区三区| 亚洲第一区二区三区不卡| 美女xxoo啪啪120秒动态图 | 一夜夜www| 天堂影院成人在线观看| 久久人人精品亚洲av| 欧美激情在线99| 99热这里只有是精品50| 亚洲欧美日韩无卡精品| 熟女人妻精品中文字幕| 亚洲专区国产一区二区| 一夜夜www| 淫妇啪啪啪对白视频| 亚洲精品在线观看二区| 精品欧美国产一区二区三| 9191精品国产免费久久| 亚洲av.av天堂| 午夜激情欧美在线| 一个人看的www免费观看视频| 国产av不卡久久| 成熟少妇高潮喷水视频| 97碰自拍视频| 午夜福利在线观看吧| 亚洲中文日韩欧美视频| 久久中文看片网| 国产日本99.免费观看| 在线观看66精品国产| 国产亚洲欧美在线一区二区| 亚洲天堂国产精品一区在线| 国产三级在线视频| 熟女人妻精品中文字幕| 全区人妻精品视频| 超碰av人人做人人爽久久| 亚洲天堂国产精品一区在线| 亚洲乱码一区二区免费版| 国产精华一区二区三区| 黄色日韩在线| 欧美日韩乱码在线| 在线国产一区二区在线| 日韩欧美三级三区| 精品午夜福利视频在线观看一区| 麻豆久久精品国产亚洲av| 岛国在线免费视频观看| 午夜日韩欧美国产| 变态另类丝袜制服| 日本在线视频免费播放| 国产av在哪里看| 久久性视频一级片| 日韩国内少妇激情av| 变态另类丝袜制服| 99热这里只有是精品50| 3wmmmm亚洲av在线观看| 成年免费大片在线观看| 好男人电影高清在线观看| 亚洲无线观看免费| 九九在线视频观看精品| 国产精品亚洲av一区麻豆| 免费黄网站久久成人精品 | 亚洲精品色激情综合| 99国产精品一区二区蜜桃av| 免费在线观看成人毛片| 国内精品久久久久精免费| 精品一区二区三区视频在线观看免费| 精品久久久久久久久亚洲 | 国产成人福利小说| 色精品久久人妻99蜜桃| 成人精品一区二区免费| 亚洲 国产 在线| 淫秽高清视频在线观看| 精品无人区乱码1区二区| 欧美潮喷喷水| 国产爱豆传媒在线观看| 国产精品日韩av在线免费观看| 给我免费播放毛片高清在线观看| 色尼玛亚洲综合影院| 免费无遮挡裸体视频| 波多野结衣高清作品| 欧美成狂野欧美在线观看| 成人欧美大片| 午夜福利在线观看免费完整高清在 | 熟女电影av网| 国产91精品成人一区二区三区| 国产成人影院久久av| 九色国产91popny在线| 精品一区二区三区人妻视频| 亚洲精品粉嫩美女一区| 久久久色成人| 欧美日本视频| 久久亚洲精品不卡| 麻豆国产av国片精品| 免费黄网站久久成人精品 | 久久这里只有精品中国| 观看美女的网站| 色综合欧美亚洲国产小说| 女同久久另类99精品国产91| 免费高清视频大片| 99riav亚洲国产免费| 国产激情偷乱视频一区二区| www日本黄色视频网| www.www免费av| 国产伦精品一区二区三区四那| 欧美一区二区精品小视频在线| 国产欧美日韩精品亚洲av| 欧美xxxx性猛交bbbb| 在线十欧美十亚洲十日本专区| 久久99热这里只有精品18| 99在线视频只有这里精品首页| 少妇的逼水好多| 欧美日本亚洲视频在线播放| 国产伦人伦偷精品视频| 久久草成人影院| 最新在线观看一区二区三区| 久久久久国产精品人妻aⅴ院| 国产高清视频在线播放一区| 久久99热这里只有精品18| 精品人妻熟女av久视频| 在线播放无遮挡| 男女做爰动态图高潮gif福利片| 最近最新免费中文字幕在线| 久久精品国产亚洲av天美| 国产成+人综合+亚洲专区| 成人国产一区最新在线观看| 亚洲国产精品999在线| 国产三级黄色录像| 高清毛片免费观看视频网站| 日韩欧美在线二视频| 伦理电影大哥的女人| 好男人电影高清在线观看| 国内毛片毛片毛片毛片毛片| 亚洲美女搞黄在线观看 | 能在线免费观看的黄片| 精品久久久久久久久久免费视频| 在线观看av片永久免费下载| 国产黄色小视频在线观看| 99久久99久久久精品蜜桃| 精品久久国产蜜桃| 国产精品野战在线观看| 美女高潮喷水抽搐中文字幕| 成人美女网站在线观看视频| 亚州av有码| 天堂影院成人在线观看| 久久中文看片网| 99久久成人亚洲精品观看| 亚洲片人在线观看| 久久精品国产亚洲av天美| 最新中文字幕久久久久| 精品午夜福利视频在线观看一区| 亚洲精品成人久久久久久| 国产精品亚洲av一区麻豆| 人人妻人人澡欧美一区二区| xxxwww97欧美| 亚洲国产精品sss在线观看| 国产视频内射| 精品人妻视频免费看| 麻豆av噜噜一区二区三区| 日本五十路高清| 一进一出抽搐动态| 久久精品人妻少妇| 亚洲人成电影免费在线| 国产高清视频在线播放一区| 51国产日韩欧美| 免费在线观看影片大全网站| 精品一区二区三区视频在线| 午夜福利高清视频| 美女 人体艺术 gogo| 成人av一区二区三区在线看| 久99久视频精品免费| 99久久久亚洲精品蜜臀av| 好男人在线观看高清免费视频| 欧美日韩黄片免| 亚洲成av人片免费观看| 夜夜看夜夜爽夜夜摸| 美女cb高潮喷水在线观看| 成人特级黄色片久久久久久久| 亚洲五月天丁香| 一个人看的www免费观看视频| 91久久精品国产一区二区成人| 国产精品国产高清国产av| 女人十人毛片免费观看3o分钟| 亚洲,欧美,日韩| 免费观看人在逋| 久久精品国产亚洲av涩爱 | 全区人妻精品视频| 精品欧美国产一区二区三| 亚洲一区二区三区色噜噜| 真人一进一出gif抽搐免费| 动漫黄色视频在线观看| 嫁个100分男人电影在线观看| 久久国产精品影院| 亚洲成av人片免费观看| 成人毛片a级毛片在线播放| 久9热在线精品视频| 国产精品一区二区性色av| 国产欧美日韩一区二区精品| 91字幕亚洲| 99热6这里只有精品| 成人国产一区最新在线观看| 最新在线观看一区二区三区| 中文字幕人成人乱码亚洲影| 免费黄网站久久成人精品 | 中文在线观看免费www的网站| 国产蜜桃级精品一区二区三区| 久久久久精品国产欧美久久久| 我的女老师完整版在线观看| 国产色爽女视频免费观看| 精品乱码久久久久久99久播| 欧美一级a爱片免费观看看| 99在线人妻在线中文字幕| 99久久精品一区二区三区| 国产三级黄色录像| 欧美日韩国产亚洲二区| 日日干狠狠操夜夜爽| 成年版毛片免费区| 国产视频内射| 免费看a级黄色片| 国产视频内射| 国产不卡一卡二| 欧美国产日韩亚洲一区| 国产精品免费一区二区三区在线| 欧美最黄视频在线播放免费| 九九热线精品视视频播放| 在线观看舔阴道视频| 嫩草影院新地址| 国产探花极品一区二区| 午夜福利在线观看吧| 日本黄色视频三级网站网址| 真人一进一出gif抽搐免费| 成人av一区二区三区在线看| 国产视频内射| 国产真实伦视频高清在线观看 | 黄色配什么色好看| 18禁裸乳无遮挡免费网站照片| 欧美日韩黄片免| 一进一出抽搐动态| 久久国产乱子免费精品| 国产av一区在线观看免费| 亚洲av美国av| 午夜日韩欧美国产| 最近视频中文字幕2019在线8| 91九色精品人成在线观看| 观看美女的网站| 成年版毛片免费区| 欧美成狂野欧美在线观看| av在线蜜桃| 久99久视频精品免费| av中文乱码字幕在线| 特大巨黑吊av在线直播| 日韩欧美一区二区三区在线观看| 免费在线观看日本一区| 美女xxoo啪啪120秒动态图 | 一级作爱视频免费观看| 婷婷精品国产亚洲av在线| 亚洲一区高清亚洲精品| 99久久精品热视频| 欧美最黄视频在线播放免费| 超碰av人人做人人爽久久| 激情在线观看视频在线高清| 18禁裸乳无遮挡免费网站照片| 人妻夜夜爽99麻豆av| 很黄的视频免费| 国产aⅴ精品一区二区三区波| 宅男免费午夜| 看片在线看免费视频| 日韩欧美免费精品| 国语自产精品视频在线第100页| 日本一二三区视频观看| 国产精品久久久久久亚洲av鲁大| 欧美极品一区二区三区四区| 亚洲av免费在线观看| 波多野结衣巨乳人妻| 一夜夜www| 麻豆成人av在线观看| 一本一本综合久久| 日本与韩国留学比较| 精品日产1卡2卡| 宅男免费午夜| 国产精品人妻久久久久久| 亚洲欧美日韩卡通动漫| 999久久久精品免费观看国产| 中文字幕精品亚洲无线码一区| 精华霜和精华液先用哪个| 亚洲第一电影网av| 99在线人妻在线中文字幕| 国产伦在线观看视频一区| 天天躁日日操中文字幕| 1000部很黄的大片| 免费看日本二区| 我的女老师完整版在线观看| 日韩欧美在线乱码| 男女视频在线观看网站免费| 亚洲第一区二区三区不卡| 很黄的视频免费| 午夜a级毛片| 人妻丰满熟妇av一区二区三区| 极品教师在线免费播放| 两性午夜刺激爽爽歪歪视频在线观看| 一本久久中文字幕| 国内少妇人妻偷人精品xxx网站| 久久草成人影院| 久久中文看片网| 午夜日韩欧美国产| 久久久久久久精品吃奶| 国产一级毛片七仙女欲春2| 非洲黑人性xxxx精品又粗又长| 午夜福利欧美成人| 脱女人内裤的视频| 免费观看人在逋| 九九在线视频观看精品| 51国产日韩欧美| 日韩欧美精品v在线| 亚洲第一电影网av| 欧美一区二区国产精品久久精品| 老熟妇乱子伦视频在线观看| 久久久国产成人免费| 欧美又色又爽又黄视频| 丰满人妻熟妇乱又伦精品不卡| 国产一区二区激情短视频| 欧美激情久久久久久爽电影| 亚洲成av人片免费观看| 国产真实伦视频高清在线观看 | 精品不卡国产一区二区三区| 国产精品美女特级片免费视频播放器| 男女下面进入的视频免费午夜| 脱女人内裤的视频| 成人午夜高清在线视频| 不卡一级毛片| 国产v大片淫在线免费观看| 久久久久久久久久黄片| 午夜免费男女啪啪视频观看 | 88av欧美| 欧美高清成人免费视频www| 免费人成视频x8x8入口观看| 国产精品综合久久久久久久免费| 欧美日韩亚洲国产一区二区在线观看| 亚洲熟妇熟女久久| 欧美绝顶高潮抽搐喷水| 日韩亚洲欧美综合| 97人妻精品一区二区三区麻豆| 精品久久久久久久久久免费视频| 免费高清视频大片| 欧美最黄视频在线播放免费| 亚洲七黄色美女视频| 国内精品美女久久久久久| 久久草成人影院| 我的女老师完整版在线观看| netflix在线观看网站| 欧美成人一区二区免费高清观看| 好男人电影高清在线观看| 亚洲精品乱码久久久v下载方式| 国产一区二区亚洲精品在线观看| 国语自产精品视频在线第100页| 在现免费观看毛片| eeuss影院久久| 久久久久久久亚洲中文字幕 | h日本视频在线播放| 亚洲人与动物交配视频| 日韩欧美精品v在线| 十八禁国产超污无遮挡网站| 国产精品女同一区二区软件 | 久久中文看片网| 成人永久免费在线观看视频| 午夜免费激情av| 深夜a级毛片| 国产色婷婷99| 欧美一区二区亚洲| 在线天堂最新版资源| 久99久视频精品免费| 婷婷六月久久综合丁香| 亚洲国产日韩欧美精品在线观看| 国产精华一区二区三区| av国产免费在线观看| 成人欧美大片| 久久天躁狠狠躁夜夜2o2o| av中文乱码字幕在线| 国产精品一区二区免费欧美| 两个人视频免费观看高清| 久久久久久九九精品二区国产| 久久精品国产亚洲av天美| 最近视频中文字幕2019在线8| 成人特级黄色片久久久久久久| 久久精品国产99精品国产亚洲性色| 精品午夜福利在线看| 日本熟妇午夜| 国产淫片久久久久久久久 | 在线看三级毛片| 88av欧美| 老熟妇乱子伦视频在线观看| 亚洲人成伊人成综合网2020| 亚洲av.av天堂| 久久精品国产亚洲av涩爱 | 成人三级黄色视频| 可以在线观看毛片的网站| 久久天躁狠狠躁夜夜2o2o| 2021天堂中文幕一二区在线观| av在线老鸭窝| 久久6这里有精品| 男人舔奶头视频| 国产高清三级在线| 99国产精品一区二区蜜桃av| 此物有八面人人有两片| 日本一本二区三区精品| 91久久精品国产一区二区成人| 神马国产精品三级电影在线观看| 国产高清三级在线| 亚洲国产精品成人综合色| 91午夜精品亚洲一区二区三区 | 亚洲精品在线美女| 精品99又大又爽又粗少妇毛片 | 国产亚洲欧美98| 99在线人妻在线中文字幕| av天堂中文字幕网| 夜夜夜夜夜久久久久| 国产国拍精品亚洲av在线观看| av天堂中文字幕网| 国产精品伦人一区二区| 丰满乱子伦码专区| 国产午夜福利久久久久久|