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

    軸盤制動(dòng)對(duì)高速列車車輪多邊形磨耗的影響

    2022-03-31 05:48:24康熙陳光雄朱琪董丙杰
    表面技術(shù) 2022年3期
    關(guān)鍵詞:有限元振動(dòng)模型

    康熙,陳光雄,朱琪,董丙杰

    摩擦磨損與潤(rùn)滑

    軸盤制動(dòng)對(duì)高速列車車輪多邊形磨耗的影響

    康熙1,2,陳光雄1,2,朱琪1,2,董丙杰1,2

    (1.西南交通大學(xué) 機(jī)械工程學(xué)院,成都 610031;2.摩擦學(xué)研究所,成都 610031)

    研究高速列車軸盤制動(dòng)引起車輪多邊形磨耗的形成機(jī)理,并提出相應(yīng)的抑制措施?;谀Σ磷约ふ駝?dòng)引起車輪多邊形磨耗的觀點(diǎn),建立高速列車拖車輪對(duì)-軸盤制動(dòng)-軌道系統(tǒng)的有限元模型。采用復(fù)特征值法,分析制動(dòng)工況下制動(dòng)盤和制動(dòng)片摩擦激勵(lì)的振動(dòng)。根據(jù)等效阻尼比判斷摩擦自激振動(dòng)的不穩(wěn)定性,等效阻尼比越小,則不穩(wěn)定振動(dòng)發(fā)生趨勢(shì)越強(qiáng)。當(dāng)?shù)刃ё枘岜刃∮讪C0.001時(shí),不穩(wěn)定振動(dòng)的振幅會(huì)克服系統(tǒng)阻尼逐漸增大。為了考慮模型中非線性因素的影響,采用瞬時(shí)動(dòng)態(tài)仿真,獲得制動(dòng)時(shí)輪軌間的法向接觸力,通過(guò)功率譜密度分析,獲得輪軌振動(dòng)主頻。此外,分析軸盤制動(dòng)系統(tǒng)安裝位置和3種類型的制動(dòng)片對(duì)車輪多邊形磨耗的影響。軸盤制動(dòng)系統(tǒng)摩擦制動(dòng)容易激勵(lì)出637 Hz左右的不穩(wěn)定振動(dòng),由于復(fù)特征值分析與瞬時(shí)動(dòng)態(tài)分析求解方法不同,因此該不穩(wěn)定振動(dòng)頻率的計(jì)算結(jié)果存在6%左右的相對(duì)誤差。軸盤制動(dòng)系統(tǒng)的安裝位置對(duì)于不穩(wěn)定振動(dòng)的發(fā)生趨勢(shì)具有重要影響,考慮到軸盤制動(dòng)系統(tǒng)實(shí)際安裝空間,當(dāng)制動(dòng)壓力角為–10°~10°時(shí),637 Hz左右的振動(dòng)對(duì)應(yīng)的等效阻尼比隨壓力角的增大而減小。采用多個(gè)蜂窩狀制動(dòng)單元組成的制動(dòng)片,在制動(dòng)時(shí)可引起602 Hz左右的不穩(wěn)定振動(dòng)。當(dāng)制動(dòng)片表面存在復(fù)合溝槽結(jié)構(gòu)時(shí),在550~650 Hz內(nèi),沒(méi)有等效阻尼比小于–0.001的不穩(wěn)定振動(dòng)。當(dāng)高速列車運(yùn)行速度為300 km/h時(shí),軸盤摩擦制動(dòng)引起的637 Hz左右的不穩(wěn)定振動(dòng)可通過(guò)輪對(duì)傳導(dǎo)至輪軌系統(tǒng)中,引起輪軌摩擦功周期性波動(dòng),從而導(dǎo)致拖車車輪發(fā)生22~23階多邊形磨耗。在滿足制動(dòng)系統(tǒng)安裝要求的條件下,適當(dāng)增大壓力角,可減輕由軸盤制動(dòng)引起的車輪多邊形磨耗。采用多個(gè)蜂窩狀制動(dòng)單元組成的制動(dòng)片,容易導(dǎo)致拖車車輪發(fā)生20~21階多邊形磨耗。在制動(dòng)片表面添加復(fù)合溝槽結(jié)構(gòu),可抑制由軸盤制動(dòng)引起的19~23階車輪多邊形磨耗。

    軸盤制動(dòng)系統(tǒng);車輪多邊形磨耗;摩擦自激振動(dòng);數(shù)值仿真;高速列車;表面溝槽

    掃碼查看文章講解

    近年來(lái),我國(guó)高速鐵路建設(shè)蓬勃發(fā)展,但與此同時(shí),多樣、復(fù)雜化的列車運(yùn)行環(huán)境導(dǎo)致高鐵列車車輪踏面損傷[1]現(xiàn)象越來(lái)越嚴(yán)重。車輪多邊形磨耗,也被稱為車輪非圓磨耗或橢圓形磨耗,指新造或鏇修后的車輪使用一段時(shí)間后,車輪與鋼軌滾動(dòng)接觸的踏面周向出現(xiàn)不均勻磨耗的現(xiàn)象。高鐵車輪多邊形磨耗引起的強(qiáng)迫振動(dòng)不僅會(huì)影響乘坐舒適性,而且當(dāng)振動(dòng)頻率接近列車零部件的固有頻率時(shí),會(huì)產(chǎn)生共振,降低車軸、軸承等零部件的使用壽命,這對(duì)于高速行駛的列車而言,具有很大的安全隱患[2]。在車輪多邊形治理方面,運(yùn)營(yíng)公司普遍采用鏇修的方法來(lái)消除車輪多邊形磨耗,但該方法不僅會(huì)縮短車輪壽命,增加運(yùn)營(yíng)成本,而且鏇修后的車輪在投入使用一段時(shí)間后,仍會(huì)出現(xiàn)多邊形磨耗現(xiàn)象。只有充分認(rèn)識(shí)車輪多邊形磨耗的形成機(jī)理及其影響因素,才能從根本上解決問(wèn)題。

    關(guān)于車輪多邊形磨耗的形成機(jī)理,國(guó)內(nèi)外學(xué)者進(jìn)行了詳細(xì)的研究。由于車輪多邊形問(wèn)題涉及了車輛系統(tǒng)動(dòng)力學(xué)、輪軌接觸力學(xué)、振動(dòng)力學(xué)、材料科學(xué)、摩擦學(xué)等多學(xué)科的內(nèi)容,目前國(guó)內(nèi)外學(xué)者對(duì)于車輪多邊形磨耗的形成機(jī)理仍未達(dá)成共識(shí)[3]。Bommundt[4]基于車輪初始不圓順與車輪轉(zhuǎn)動(dòng)慣量共同作用導(dǎo)致車輪多邊形磨耗的觀點(diǎn),采用攝動(dòng)技術(shù)進(jìn)行了研究。結(jié)果表明,車輛速度越快,車輪非圓浪涌低次諧波形成越快。Meywerk[5]采用輪對(duì)-鋼軌彈性模型,對(duì)車輪廓面不圓順的發(fā)展進(jìn)行了研究,發(fā)現(xiàn)輪對(duì)左右車輪廓面的多邊形相位差越大,車輪多邊形形成越快,輪對(duì)的第一、二階彎曲振動(dòng)模態(tài)對(duì)車輪多邊形磨耗有重要影響。Meinders等[6]建立了輪對(duì)-鋼軌彈性系統(tǒng)接觸模型,使用車輪磨損反饋圈,將初始模型推廣為長(zhǎng)期磨損模型,研究了高速車輛車輪初始不圓順和輪對(duì)失衡對(duì)車輪多邊形磨耗的影響,發(fā)現(xiàn)車輪初始不平順對(duì)車輪多邊形磨耗的發(fā)展具有重要影響。初始二階和二階以上的高階車輪多邊形相對(duì)穩(wěn)定,靜態(tài)、動(dòng)態(tài)不平衡對(duì)車輪多邊形磨耗的發(fā)展影響較小。Johansson[7]對(duì)瑞典鐵路所使用的99種運(yùn)營(yíng)里程超過(guò)100 000 km的車輪進(jìn)行了車輪不圓度實(shí)測(cè),發(fā)現(xiàn)客車、貨車、通勤車、地鐵車輛均出現(xiàn)了不同程度的車輪多邊形磨耗,其中高速客車最為嚴(yán)重。在加工過(guò)程中,三角卡盤固定不當(dāng)會(huì)導(dǎo)致部分地鐵車輪形成初始多邊形。Jin等[8]通過(guò)試驗(yàn)與理論計(jì)算結(jié)合的方法,對(duì)地鐵車輛車輪多邊形磨耗的形成機(jī)理進(jìn)行了研究,發(fā)現(xiàn)72 Hz左右的輪對(duì)一階彎曲共振導(dǎo)致了地鐵車輪9階多邊形磨耗。Ma等[9]通過(guò)理論分析和動(dòng)力學(xué)仿真對(duì)地鐵車輪多邊形磨耗的產(chǎn)生機(jī)理進(jìn)行了研究,提出車輪滾動(dòng)多周的振動(dòng)導(dǎo)致地鐵車輪9階多邊形磨耗的觀點(diǎn)。Tao等[10]通過(guò)試驗(yàn)與數(shù)值模型分析對(duì)地鐵車輛車輪多邊形磨耗機(jī)理進(jìn)行了研究,研究發(fā)現(xiàn),P2共振是導(dǎo)致車輪形成5~8階多邊形的原因,通過(guò)改善制動(dòng)閘片與車輪踏面的匹配關(guān)系、降低牽引力、改善鋼軌焊接處不平順、降低扣件剛度等可抑制車輪多邊形磨耗。趙曉男等[11-12]建立了輪對(duì)-軌道系統(tǒng)有限元模型,對(duì)多邊形車輪形成機(jī)理進(jìn)行了研究,認(rèn)為當(dāng)高速線路發(fā)生制動(dòng)滑動(dòng)時(shí),輪軌間的蠕滑力趨于飽和,可能會(huì)導(dǎo)致輪軌系統(tǒng)發(fā)生摩擦自激振動(dòng),進(jìn)而引起車輪多邊形磨耗。Wu等[13]采用有限元方法研究了制動(dòng)引起的不穩(wěn)定振動(dòng)對(duì)輪軌磨耗的影響,研究表明,列車制動(dòng)可能引起車輪踏面產(chǎn)生波浪型磨耗。Ma等[14]采用車輛-軌道多體動(dòng)力學(xué),仿真研究了高速列車高階車輪多邊形磨耗的形成和發(fā)展,并分析了季節(jié)和車輪在車輛縱向分布位置的影響。研究發(fā)現(xiàn),同一轉(zhuǎn)向架輪對(duì)間鋼軌的三階彎曲模態(tài)可導(dǎo)致高階車輪多邊形磨耗的形成,相鄰轉(zhuǎn)向架輪對(duì)間鋼軌的彎曲共振會(huì)加速該磨耗的進(jìn)一步發(fā)展。此外,多雨和高溫導(dǎo)致了車輪多邊形磨耗在夏季比其他季節(jié)更嚴(yán)重。列車在包含曲線區(qū)段的高速線路上往返運(yùn)營(yíng)導(dǎo)致了第一位和第四位的車輪多邊形磨耗更為嚴(yán)重。

    在先前研究中,軸盤制動(dòng)系統(tǒng)有限元模型被用于預(yù)測(cè)可導(dǎo)致制動(dòng)尖叫噪聲的摩擦自激振動(dòng)。通過(guò)對(duì)比仿真結(jié)果與現(xiàn)場(chǎng)測(cè)試[15],分析、驗(yàn)證了該模型的準(zhǔn)確性和可靠性。本文以哈爾濱-大連客運(yùn)專線上運(yùn)行的CRH3型動(dòng)車組拖車為原型,建立了直線區(qū)段高速鐵路拖車輪對(duì)-軸盤制動(dòng)-軌道系統(tǒng)的有限元模型,對(duì)制動(dòng)工況下軸盤制動(dòng)系統(tǒng)中因摩擦激勵(lì)的振動(dòng)與車輪多邊形磨耗之間的關(guān)系進(jìn)行了研究。研究發(fā)現(xiàn),制動(dòng)盤-制動(dòng)片摩擦自激振動(dòng)可通過(guò)輪對(duì)傳遞至輪軌系統(tǒng)中,從而引起輪軌摩擦功周期性變化,導(dǎo)致車輪多邊形磨耗。

    1 理論模型和仿真模型

    1.1 拖車輪對(duì)-軸盤制動(dòng)-軌道系統(tǒng)的接觸模型

    直線區(qū)段高速列車拖車輪對(duì)-軸盤制動(dòng)-軌道系統(tǒng)的接觸模型如圖1所示。圖1中,SVL、SVR分別為輪對(duì)左右兩端的垂向懸掛力;L、R和L、R分別為左右輪軌間的法向接觸力和蠕滑力;L、R分別為左右輪軌間的接觸角;RL、RL和RV、RV分別為扣件對(duì)鋼軌橫向和垂向的支撐剛度、阻尼;F、F分別為地基對(duì)軌道板的垂向支撐剛度、阻尼,部分參數(shù)的具體數(shù)值[13]見(jiàn)表1。

    圖1 拖車輪對(duì)-軌道系統(tǒng)的接觸模型

    表1 拖車輪對(duì)-軸盤制動(dòng)-軌道系統(tǒng)具體參數(shù)

    Tab.1 Parameters of finite element model of a trailer wheelset-track system with an axle-mounted disc brake system

    1.2 拖車輪對(duì)-軸盤制動(dòng)-軌道系統(tǒng)的有限元模型

    在Abaqus中建立的直線線路上高速列車拖車輪對(duì)-軸盤制動(dòng)-軌道系統(tǒng)的有限元模型如圖2所示。該模型主要由拖車輪對(duì)、制動(dòng)盤、制動(dòng)片、閘片托、制動(dòng)支架、鋼軌和軌道板組成,網(wǎng)格單元類型為C3D8I。在制動(dòng)系統(tǒng)中,起連接作用的銷部件采用鉸鏈單元模擬。制動(dòng)盤與制動(dòng)片之間相互接觸,制動(dòng)片與閘片托之間采用綁定約束。閘片托與制動(dòng)支架前端通過(guò)鉸鏈單元連接,僅保留兩部件之間的相對(duì)轉(zhuǎn)動(dòng)自由度。制動(dòng)支架中間孔與機(jī)架之間通過(guò)鉸鏈單元連接,除轉(zhuǎn)動(dòng)以外的其他自由度均被約束,制動(dòng)力作用在支架末端的銷孔上,制動(dòng)時(shí),通過(guò)杠桿原理,將力放大后傳遞到制動(dòng)片上。車輪和鋼軌、制動(dòng)盤和制動(dòng)片之間的接觸關(guān)系均為摩擦關(guān)系??奂偷来驳闹蝿偠?、阻尼由彈簧、阻尼單元模擬,有限元模型中的材料參數(shù)[13]見(jiàn)表2。

    圖2 輪對(duì)-軸盤制動(dòng)-軌道系統(tǒng)的有限元模型

    表2 有限元模型的材料參數(shù)

    Tab.2 Material parameters of finite element model

    1.3 摩擦自激振動(dòng)的復(fù)特征值分析法

    采用Abaqus軟件對(duì)制動(dòng)工況下高速列車拖車輪對(duì)-軸盤制動(dòng)-軌道系統(tǒng)的穩(wěn)定性進(jìn)行復(fù)特征值分析[16-17],在摩擦耦合的作用下,系統(tǒng)的運(yùn)動(dòng)方程為[18-20]:

    式中:和分別為運(yùn)動(dòng)方程式(1)的特征值和特征向量。利用QZ法對(duì)式(2)進(jìn)行求解,可得通解為:

    2 結(jié)果與討論

    2.1 拖車輪對(duì)-軸盤制動(dòng)-軌道系統(tǒng)的不穩(wěn)定振動(dòng)分布及模態(tài)

    高速列車在正常運(yùn)行時(shí),車輪與鋼軌間存在較小的滾動(dòng)摩擦力。列車制動(dòng)時(shí),一般情況下車輪與鋼軌之間不發(fā)生滑動(dòng),輪軌間仍為滾動(dòng)摩擦,而制動(dòng)片與制動(dòng)盤之間存在較大摩擦力。為了排除車輛在制動(dòng)工況下輪軌間摩擦力的影響,在有限元模型中設(shè)置拖車輪對(duì)和軸盤制動(dòng)系統(tǒng)的平移速度為300 km/h,車輪半徑為0.46 m,輪對(duì)的轉(zhuǎn)動(dòng)速度為181.159 rad/s,使輪軌間縱向蠕滑率為0。通過(guò)復(fù)特征值分析,計(jì)算了制動(dòng)時(shí)拖車輪對(duì)-軸盤制動(dòng)-軌道系統(tǒng)可能發(fā)生的摩擦自激振動(dòng)在頻域上的分布情況,如圖3所示。637 Hz左右的不穩(wěn)定振動(dòng)對(duì)應(yīng)的負(fù)等效阻尼比最小為–0.051 55,故可認(rèn)為列車在制動(dòng)時(shí)這個(gè)不穩(wěn)定振動(dòng)模態(tài)最容易被激發(fā),即制動(dòng)系統(tǒng)最容易產(chǎn)生頻率為637 Hz左右的不穩(wěn)定振動(dòng)。該不穩(wěn)定振動(dòng)模態(tài)如圖4所示,鋼軌無(wú)明顯變化,變形主要發(fā)生在輪對(duì)上,其中軸盤沿車軸方向的變形最為嚴(yán)重,說(shuō)明制動(dòng)引起的不穩(wěn)定振動(dòng)主要發(fā)生在軸盤制動(dòng)系統(tǒng)中。該振動(dòng)可通過(guò)安裝在車軸上的制動(dòng)盤傳導(dǎo)至車輪上,從而引起輪軌系統(tǒng)發(fā)生不穩(wěn)定振動(dòng),導(dǎo)致車輪多邊形磨耗。

    圖3 不穩(wěn)定振動(dòng)的分布

    圖4 不穩(wěn)定振動(dòng)模態(tài)

    輪軌表面的損傷問(wèn)題通常被認(rèn)為是由輪軌間摩擦功的變化導(dǎo)致的,文獻(xiàn)[22]提出的磨損公式為:

    式中:為單位時(shí)間內(nèi)的磨損量;為磨損常數(shù);為摩擦功;為長(zhǎng)時(shí)摩擦功,是恒定常數(shù)。根據(jù)式(5)可知,單位時(shí)間內(nèi)輪軌的磨損量主要由摩擦功決定,它的計(jì)算公式為:

    式中:為輪軌間法向接觸力;Δ為輪軌間相對(duì)滑移速度;、可通過(guò)牽引-滑移率試驗(yàn)曲線獲得,均為常數(shù)。由式(6)可知,在恒定縱向蠕滑力和相對(duì)滑移速度的假設(shè)下,輪軌間法向接觸力的周期性波動(dòng)會(huì)導(dǎo)致摩擦功以相同的頻率發(fā)生周期性變化,進(jìn)而導(dǎo)致車輪踏面發(fā)生非均勻磨耗。因此,為了進(jìn)一步分析制動(dòng)時(shí)制動(dòng)片與制動(dòng)盤間的摩擦力誘發(fā)的自激振動(dòng)在輪軌系統(tǒng)中的傳導(dǎo)情況,對(duì)拖車輪對(duì)在直線線路上的軸盤制動(dòng)過(guò)程進(jìn)行瞬時(shí)動(dòng)態(tài)分析,得到0.2 s內(nèi)左右輪軌間法向接觸力的變化情況,如圖5所示。通過(guò)功率譜密度分析,得出左右輪軌間法向接觸力的振動(dòng)主頻在600 Hz左右,如圖6所示。將瞬時(shí)動(dòng)態(tài)分析獲得的輪軌振動(dòng)主頻與復(fù)特征值分析獲得的最易發(fā)生的不穩(wěn)定振動(dòng)頻率進(jìn)行對(duì)比,發(fā)現(xiàn)這兩種方法預(yù)測(cè)的結(jié)果比較接近。由于復(fù)特征值分析未考慮非線性因素,而瞬時(shí)動(dòng)態(tài)分析是通過(guò)顯式時(shí)間積分求解系統(tǒng)的動(dòng)態(tài)響應(yīng),因此這兩種方法的計(jì)算結(jié)果存在6%左右的相對(duì)誤差。

    圖5 輪軌間法向接觸力

    圖6 輪軌間法向接觸力PSD分析結(jié)果

    一般情況下,直線區(qū)段上高速列車的運(yùn)行速度取300 km/h,拖車車輪的名義滾動(dòng)圓半徑為0.46 m,則車輪多邊形磨耗階數(shù)的計(jì)算公式為:

    當(dāng)頻率為637 Hz左右時(shí),通過(guò)計(jì)算可以得到對(duì)應(yīng)的車輪多邊形磨耗階數(shù),即該不穩(wěn)定振動(dòng)頻率可能引起車輪22~23階多邊形磨耗。統(tǒng)計(jì)結(jié)果[23]顯示,當(dāng)高速列車的運(yùn)行速度為300 km/h時(shí),車輪多邊形引起的振動(dòng)的主要頻率為550~650 Hz,對(duì)應(yīng)車輪踏面上19~23階多邊形磨耗,與仿真結(jié)果中預(yù)測(cè)的不穩(wěn)定振動(dòng)頻率近乎一致。故可認(rèn)為本文的車輪多邊形磨耗預(yù)測(cè)結(jié)果與實(shí)際情況較為接近,即預(yù)測(cè)模型具有較高的可靠性。

    2.2 制動(dòng)系統(tǒng)安裝位置對(duì)高速列車拖車車輪多邊形磨耗的影響

    文獻(xiàn)[24]指出,制動(dòng)壓力角對(duì)制動(dòng)系統(tǒng)的不穩(wěn)定振動(dòng)發(fā)生趨勢(shì)有重要影響。壓力角是指制動(dòng)片壓力線與水平線之間的夾角,規(guī)定順時(shí)針為正。為了研究制動(dòng)系統(tǒng)安裝位置對(duì)高速列車拖車車輪多邊形磨耗的影響,根據(jù)制動(dòng)系統(tǒng)實(shí)際安裝空間,分別取壓力角為10°、5°、–5°、–10°。不同制動(dòng)壓力角對(duì)應(yīng)系統(tǒng)的不穩(wěn)定振動(dòng)模態(tài)分布如圖7所示。從圖7中可以看出,最可能發(fā)生的不穩(wěn)定振動(dòng)頻率均在637 Hz左右,對(duì)應(yīng)的負(fù)等效阻尼比分別為–0.029 19、–0.042 89、–0.056 62、–0.062 23,即負(fù)等效阻尼比隨著壓力角的增大而減小。當(dāng)=10°時(shí),637 Hz左右的不穩(wěn)定振動(dòng)發(fā)生趨勢(shì)最弱,在–10°~10°的范圍內(nèi)。壓力角為10°時(shí),制動(dòng)引起車輪22~24階多邊形磨耗的可能性最低。

    圖7 制動(dòng)系統(tǒng)安裝位置對(duì)不穩(wěn)定振動(dòng)分布的影響

    2.3 不同類型的制動(dòng)片對(duì)高速列車拖車車輪多邊形磨耗的影響

    為了研究不同類型的制動(dòng)片對(duì)高速列車拖車車輪多邊形磨耗的影響,分別研究了采用A、B、C 3種類型制動(dòng)片的軸盤制動(dòng)系統(tǒng)在制動(dòng)時(shí)的不穩(wěn)定振動(dòng)分布情況,制動(dòng)片和閘片托的有限元模型如圖8所示。其中,B型制動(dòng)片由多個(gè)蜂窩狀制動(dòng)單元組成,C型制動(dòng)片由A型制動(dòng)片改造(在表面添加復(fù)合溝槽結(jié)構(gòu))而成。實(shí)驗(yàn)研究[25]表明,溝槽結(jié)構(gòu)會(huì)擾亂摩擦系統(tǒng)連續(xù)的自激振動(dòng),對(duì)制動(dòng)尖叫噪聲有良好的抑制作用。

    當(dāng)車輛以300 km/h的速度運(yùn)行時(shí),采用3種不同制動(dòng)片對(duì)應(yīng)的系統(tǒng)在制動(dòng)工況下的不穩(wěn)定振動(dòng)頻率分布如圖9所示。從圖9中可以看出,采用B型制動(dòng)片系統(tǒng)的不穩(wěn)定振動(dòng)頻率(602 Hz)對(duì)應(yīng)的負(fù)等效阻尼比最小,發(fā)生的可能性最高。通過(guò)公式(7)可以計(jì)算出,該不穩(wěn)定振動(dòng)頻率對(duì)應(yīng)的車輪多邊形磨耗階數(shù)為20~21階,即采用B型制動(dòng)片的軸盤制動(dòng)系統(tǒng)在制動(dòng)時(shí)可導(dǎo)致拖車車輪20~21階多邊形磨耗。采用C型制動(dòng)片的系統(tǒng)制動(dòng)時(shí),在高速列車車輪多邊形主要振動(dòng)頻率范圍(550~650 Hz)內(nèi),沒(méi)有等效阻尼比小于–0.001的不穩(wěn)定振動(dòng),即制動(dòng)系統(tǒng)采用C型制動(dòng)片可抑制由軸盤制動(dòng)引起的高速列車車輪19~23階多邊形磨耗。

    圖8 3種類型制動(dòng)片和閘片托的有限元模型

    圖9 不同類型的制動(dòng)片對(duì)不穩(wěn)定振動(dòng)分布的影響

    3 結(jié)論

    1)軸盤制動(dòng)容易引起頻率為637 Hz左右的不穩(wěn)定振動(dòng)。當(dāng)高速列車運(yùn)行速度為300 km/h時(shí),可能導(dǎo)致拖車車輪發(fā)生22~23階多邊形磨耗。

    2)軸盤制動(dòng)系統(tǒng)的安裝位置對(duì)高速列車拖車車輪多邊形磨耗有重要影響。在–10°~10°內(nèi),適當(dāng)增大制動(dòng)壓力角,可減緩由軸盤制動(dòng)引起的車輪多邊形磨耗的發(fā)生。

    3)由多個(gè)蜂窩狀制動(dòng)單元組成的B型制動(dòng)片在制動(dòng)時(shí)容易引起602 Hz左右的不穩(wěn)定振動(dòng),導(dǎo)致拖車車輪發(fā)生20~21階多邊形磨耗。制動(dòng)系統(tǒng)采用帶有復(fù)合溝槽結(jié)構(gòu)的C型制動(dòng)片,可抑制由摩擦自激振動(dòng)引起的19~23階車輪多邊形磨耗。

    本文主要采有限元仿真研究了軸盤制動(dòng)對(duì)車輪多邊形磨耗的影響。在下一階段的工作中,將采用縮尺輪軌實(shí)驗(yàn)對(duì)結(jié)論進(jìn)一步驗(yàn)證。

    [1] 王延朋, 丁昊昊, 鄒強(qiáng), 等. 列車車輪踏面滾動(dòng)接觸疲勞研究進(jìn)展[J]. 表面技術(shù), 2020, 49(5): 120-128.

    WANG Yan-peng, DING Hao-hao, ZOU Qiang, et al. Research Progress on Rolling Contact Fatigue of Railway Wheel Treads[J]. Surface Technology, 2020, 49(5): 120- 128.

    [2] 金學(xué)松, 吳越, 梁樹林, 等. 高速列車車輪多邊形磨耗、機(jī)理、影響和對(duì)策分析[J]. 機(jī)械工程學(xué)報(bào), 2020, 56(16): 118-136.

    JIN Xue-song, WU Yue, LIANG Shu-lin, et al. Character-istics, Mechanism, Influences and Countermeasures of Polygonal Wear of High-Speed Train Wheels[J]. Journal of Mechanical Engineering, 2020, 56(16): 118-136.

    [3] 朱海燕, 胡華濤, 尹必超, 等. 軌道車輛車輪多邊形研究進(jìn)展[J]. 交通運(yùn)輸工程學(xué)報(bào), 2020, 20(1): 102-119.

    ZHU Hai-yan, HU Hua-tao, YIN Bi-chao, et al. Research Progress on Wheel Polygons of Rail Vehicles[J]. Journal of Traffic and Transportation Engineering, 2020, 20(1): 102-119.

    [4] BROMMUNDT E. A Simple Mechanism for the Poly-gonalization of Railway Wheels by Wear[J]. Mechanics Research Communications, 1997, 24(4): 435-442.

    [5] MEYWERK M. Polygonalization of Railway Wheels[J]. Archive of Applied Mechanics, 1999, 69(2): 105-120.

    [6] MEINDERS T, MEINKE P. Rotor Dynamics and Irregular Wear of Elastic Wheelsets[M]. Berlin: Springer, 2003.

    [7] JOHANSSON A. Out-of-Round Railway Wheels—Asses-s-ment of Wheel Tread Irregularities in Train Traffic[J]. Journal of Sound and Vibration, 2006, 293(3-5): 795-806.

    [8] JIN Xue-song, WU Lei, FANG Jian-ying, et al. An Inv-e-stigation into the Mechanism of the Polygonal Wear of Metro Train Wheels and Its Effect on the Dynamic Beh-aviour of a Wheel/Rail System[J]. Vehicle System Dyna-mics, 2012, 50(12): 1817-1834.

    [9] MA Wei-hua, SONG Rong-rong, LUO Shi-hui. Study on the Mechanism of the Formation of Polygon-Shaped Wh-eels on Subway Vehicles[J]. Proceedings of the Institution of Mechanical Engineers, Part F: Journal of Rail and Ra-pid Transit, 2016, 230(1): 129-137.

    [10] TAO Gong-quan, WEN Ze-feng, LIANG Xi-ren, et al. An Investigation into the Mechanism of the Out-of-Round Wheels of Metro Train and Its Mitigation Measures[J]. Vehicle System Dynamics, 2019, 57(1): 1-16.

    [11] 趙曉男, 陳光雄, 崔曉璐, 等. 高速列車車輪多邊形磨耗的形成機(jī)理及影響因素探究[J]. 表面技術(shù), 2018, 47(8): 8-13.

    ZHAO Xiao-nan, CHEN Guang-xiong, CUI Xiao-lu, et al. Formation Mechanism and Influencing Factors of the Polygonal Wear of High-Speed Train Wheels[J]. Surface Technology, 2018, 47(8): 8-13.

    [12] ZHAO X N, CHEN G X, LV J Z, et al. Study on the Mechanism for the Wheel Polygonal Wear of High-Speed Trains in Terms of the Frictional Self-Excited Vibration Theory[J]. Wear, 2019, 426-427: 1820-1827.

    [13] WU B W, QIAO Q F, CHEN G X, et al. Effect of the Un-s-table Vibration of the Disc Brake System of High-Speed Trains on Wheel Polygonalization[J]. Proceedings of the Institution of Mechanical Engineers, Part F: Journal of Rail and Rapid Transit, 2020, 234(1): 80-95.

    [14] MA Chao-zhi, GAO Liang, CUI Ri-xin, et al. The Initi-a-tion Mechanism and Distribution Rule of Wheel High- Order Polygonal Wear on High-Speed Railway[J]. Engin-eering Failure Analysis, 2021, 119: 104937.

    [15] CHEN G X, LV J Z, ZHU Q, et al. Effect of the Braking Pressure Variation on Disc Brake Squeal of a Railway Ve-hicle: Test Measurement and Finite Element Analysis [J]. Wear, 2019, 426-427: 1788-1796.

    [16] YUAN Y. An Eigenvalue Analysis Approach to Brake Squeal Problem[C]// Proc of the 29th ISATA Conf on Automotive Braking Systems. Florence: [s. n.], 1996: 3-9.

    [17] 康熙, 陳光雄, 呂金洲, 等. 縮尺輪軌模型中鋼軌波磨的相似性[J]. 西南交通大學(xué)學(xué)報(bào), 2020, 55(6): 1320- 1327.

    KANG Xi, CHEN Guang-xiong, Lü Jin-zhou, et al. Sim-ilarity of Small-Scale Wheelset-Track Model for Investig-ation of Rail Corrugation[J]. Journal of Southwest Jiao-tong University, 2020, 55(6): 1320-1327.

    [18] LOU G, WU T W, BAI Z. Disk Brake Squeal Prediction Using the ABLE Algorithm[J]. Journal of Sound and Vib-ration, 2004, 272(3-5): 731-748.

    [19] CHEN G X, ZHANG S, WU B W, et al. Field Measure-ment and Model Prediction of Rail Corrugation[J]. Pro-ceedings of the Institution of Mechanical Engineers, Part F: Journal of Rail and Rapid Transit, 2020, 234(4): 381- 392.

    [20] OUYANG Hua-jiang, NACK W, YUAN Yong-bin, et al. Numerical Analysis of Automotive Disc Brake Squeal: A Review[J]. International Journal of Vehicle Noise and Vibration, 2005, 1(3/4): 207.

    [21] 夏晨光, 陳光雄, 朱旻昊, 等. 車輪輻板阻尼涂層對(duì)鋼軌波磨的影響[J]. 表面技術(shù), 2020, 49(3): 134-140.

    XIA Chen-guang, CHEN Guang-xiong, ZHU Min-hao, et al. Effect of Damped Coating of Wheel Web on Rail Corrugation[J]. Surface Technology, 2020, 49(3): 134-140.

    [22] BROCKLEY C A, KO P L. An Investigation of Rail Corrugation Using Friction-Induced Vibration Theory[J]. Wear, 1988, 128(1): 99-106.

    [23] 李大地, 戴煥云. 基于鋼軌模態(tài)振動(dòng)的車輪高階多邊形研究[J]. 機(jī)車電傳動(dòng), 2018(2): 61-66.

    LI Da-di, DAI Huan-yun. Research on Wheel Polygoniz-ation Based on Rail Modal Oscillation Analysis[J]. Elec-tric Drive for Locomotives, 2018(2): 61-66.

    [24] 陳光雄, 戴煥云, 曾京, 等. 踏面制動(dòng)尖叫噪聲的有限元分析[J]. 機(jī)械工程學(xué)報(bào), 2009, 45(3): 203-207.

    CHEN Guang-xiong, DAI Huan-yun, ZENG Jing, et al. Finite Element Analysis of the Squeal Propensity for a Pad-Wheel Tread Brake System[J]. Journal of Mechanical Engineering, 2009, 45(3): 203-207.

    [25] 王安宇, 王東偉, 范志勇, 等. 表面溝槽與潤(rùn)滑劑協(xié)同作用對(duì)摩擦振動(dòng)和噪聲特性的影響[J]. 表面技術(shù), 2019, 48(8): 16-22.

    WANG An-yu, WANG Dong-wei, FAN Zhi-yong, et al. Synergistic Effect of Surface Groove and Lubricant on Friction-Induced Vibration and Noise Characteristics[J]. Surface Technology, 2019, 48(8): 16-22.

    Influence of Axle-Mounted Disc Brake on Polygonal Wear of High-Speed Train Wheels

    1,2,1,2,1,2,1,2

    (1. School of Mechanical Engineering, Southwest Jiaotong University, Chengdu 610031, China; 2. Tribology Research Institute, Chengdu 610031, China)

    In order to study the formation mechanism of wheel polygonal wear on high-speed trains and propose corres-ponding inhibition measures. Based on the viewpoint that the friction-induced self-excited vibration causes wheel polygonal wear, a finite element model including high-speed train trailer wheelset, axle-mounted disc brake system and track was established. The vibrations excited by the friction between discs and pads under brake conditions were analyzed by the complex eigenvalue method. According to the effective damping ratio to judge the instability of excited vibrations, the smaller the effective damping ratio, the stronger the occurrence trend of unstable vibration. When the effective damping ratio was less than –0.001, the amplitude of unstable vibrations increased gradually overcoming the system damping. The wheel-rail normal contact forces during braking were obtained by the transient dynamic simulation, which took into account the influence of nonlinear factors. And then, the main frequencies of wheel-rail vibrations were calculated by the power spectral density (PSD) analysis of the wheel-rail normal contact forces. In addition, it was investigated that the effect of the position of the axle-mounted disc brake system and three types of brake pads on wheel polygonal wear. Friction braking of axle-mounted disc brake system was easy to excite unstable vibration about 637 Hz. Because of the difference between complex eigenvalue analysis and instantaneous dynamic analysis, there was a relative error of about 6% in the calculation results of the unstable vibration frequency. The installation position of the axle-mounted disc brake system had an important influence on the occurrence trend of unstable vibration. Considering the actual installation space of the axle-mounted disc brake system, when the brake pressure angle was in the range of –10° to 10°, the equivalent damping ratio corresponding to the vibration about 637 Hz decreased with the increase of the pressure angle. Brake pads composed of several honeycomb brake units could cause unstable vibration of about 602 Hz during braking. When there was a composite groove structure on the surface of the brake pad, there was no unstable vibration with equivalent damping ratio less than –0.001 in the frequency range of 550~650 Hz. Results show that when the high-speed train runs at 300 km/h, the unstable vibration of 637 Hz caused by axle-disc friction braking can be transmitted to the wheel-rail system through the wheelset, which causes the periodic fluctuation of wheel-rail friction work and leads to 22~23 order polygonal wear of trailer wheels. Under the condition of meeting the installation requirements of braking system, increasing the pressure angle properly can reduce the wheel polygon wear caused by axle-disc braking. Using the brake pads consisted of multiple honeycomb brake units can easily cause trailer wheels to occur 20~21 order polygonal wear, and 19~23 order wheel polygonal wear caused by axle-mounted disc brake is greatly suppressed by using the brake pads with groove-textured surface.

    axle-mounted disc brake system; wheel polygonal wear; friction-induced self-excited vibration; numerical simulation; high-speed trains; surface groove

    2021-03-19;

    2021-08-30

    KANG Xi (1996—), Male, Doctoral candidate, Research focus: frictional vibration and noise.

    陳光雄(1962—),男,博士,教授,主要研究方向?yàn)槟Σ琳駝?dòng)與噪聲。

    CHEN Guang-xiong(1962—), Male, Doctor, Professor, Research focus: frictional vibration and noise.

    康熙, 陳光雄, 朱琪, 等.軸盤制動(dòng)對(duì)高速列車車輪多邊形磨耗的影響[J]. 表面技術(shù), 2022, 51(3): 43-50.

    U271.91;TH17

    A

    1001-3660(2022)03-0043-08

    10.16490/j.cnki.issn.1001-3660.2022.03.003

    2021-03-19;

    2021-08-30

    國(guó)家自然科學(xué)基金(51775461)

    Fund:The National Natural Science Foundation of China (51775461)

    康熙(1996—),男,博士研究生,主要研究方向?yàn)槟Σ琳駝?dòng)與噪聲。

    KANG Xi, CHEN Guang-xiong, ZHU Qi, et al. Influence of Axle-Mounted Disc Brake on Polygonal Wear of High-Speed Train Wheels[J]. Surface Technology, 2022, 51(3): 43-50.

    猜你喜歡
    有限元振動(dòng)模型
    一半模型
    振動(dòng)的思考
    重要模型『一線三等角』
    振動(dòng)與頻率
    重尾非線性自回歸模型自加權(quán)M-估計(jì)的漸近分布
    中立型Emden-Fowler微分方程的振動(dòng)性
    3D打印中的模型分割與打包
    磨削淬硬殘余應(yīng)力的有限元分析
    UF6振動(dòng)激發(fā)態(tài)分子的振動(dòng)-振動(dòng)馳豫
    基于SolidWorks的吸嘴支撐臂有限元分析
    国产极品天堂在线| 成人美女网站在线观看视频| 天堂俺去俺来也www色官网| 热re99久久精品国产66热6| 一区二区三区乱码不卡18| 在线观看国产h片| 激情五月婷婷亚洲| 我要看日韩黄色一级片| 久久影院123| 欧美国产精品一级二级三级 | 国产日韩欧美视频二区| 最近2019中文字幕mv第一页| 国产亚洲一区二区精品| 亚洲人成网站在线播| 超碰97精品在线观看| 久久精品国产a三级三级三级| 国产免费又黄又爽又色| 国产伦在线观看视频一区| 日本黄大片高清| 亚洲精品久久久久久婷婷小说| 久久精品国产自在天天线| 精品午夜福利在线看| 十分钟在线观看高清视频www | 桃花免费在线播放| 日韩av免费高清视频| 9色porny在线观看| 午夜福利影视在线免费观看| 亚洲欧美中文字幕日韩二区| 国产精品久久久久久精品电影小说| 一本久久精品| 一本一本综合久久| 亚洲欧美日韩卡通动漫| 久久国内精品自在自线图片| 99久国产av精品国产电影| 国产淫片久久久久久久久| 色视频www国产| 在线观看人妻少妇| 午夜福利,免费看| 熟女电影av网| 日本午夜av视频| 一级a做视频免费观看| 亚洲,一卡二卡三卡| 女性被躁到高潮视频| 成人二区视频| 国产高清三级在线| 国产伦在线观看视频一区| 你懂的网址亚洲精品在线观看| 午夜日本视频在线| 日韩伦理黄色片| 好男人视频免费观看在线| 久久国产精品大桥未久av | 夜夜看夜夜爽夜夜摸| 亚洲欧美日韩东京热| .国产精品久久| 黑丝袜美女国产一区| 两个人的视频大全免费| 69精品国产乱码久久久| 午夜日本视频在线| videos熟女内射| 美女中出高潮动态图| 久热这里只有精品99| 人人妻人人爽人人添夜夜欢视频 | 亚洲精品中文字幕在线视频 | 久久人人爽人人片av| 国产精品熟女久久久久浪| 另类精品久久| 乱人伦中国视频| 亚洲久久久国产精品| 国产av码专区亚洲av| 黄片无遮挡物在线观看| 在线观看三级黄色| 国产精品久久久久久av不卡| 亚洲国产精品成人久久小说| 999久久久精品免费观看国产| 在线观看人妻少妇| 亚洲久久久国产精品| 在线看a的网站| 最黄视频免费看| 中文精品一卡2卡3卡4更新| 亚洲精品久久午夜乱码| 丝袜脚勾引网站| www.自偷自拍.com| 国产不卡av网站在线观看| 黄色视频,在线免费观看| 国产成人影院久久av| 久久久久视频综合| 丁香六月欧美| 亚洲精品成人av观看孕妇| 精品亚洲成a人片在线观看| 人人妻人人爽人人添夜夜欢视频| 性色av乱码一区二区三区2| 国产黄频视频在线观看| 中文精品一卡2卡3卡4更新| 国产1区2区3区精品| 97精品久久久久久久久久精品| 涩涩av久久男人的天堂| 欧美黑人欧美精品刺激| 天天操日日干夜夜撸| bbb黄色大片| 青春草亚洲视频在线观看| 亚洲专区字幕在线| 12—13女人毛片做爰片一| 久久精品aⅴ一区二区三区四区| 国产精品自产拍在线观看55亚洲 | 日韩视频在线欧美| 90打野战视频偷拍视频| 午夜福利在线观看吧| 桃花免费在线播放| 中文字幕av电影在线播放| 侵犯人妻中文字幕一二三四区| 亚洲成人免费av在线播放| 十八禁网站免费在线| 丝袜人妻中文字幕| 国产成人啪精品午夜网站| 国产成人啪精品午夜网站| xxxhd国产人妻xxx| 美女国产高潮福利片在线看| 久久久久久免费高清国产稀缺| 精品一区二区三卡| 美女高潮到喷水免费观看| 成人手机av| 丰满人妻熟妇乱又伦精品不卡| 老司机深夜福利视频在线观看 | 69av精品久久久久久 | 高清黄色对白视频在线免费看| 国产精品久久久久久人妻精品电影 | 午夜福利视频精品| 一级黄色大片毛片| 午夜影院在线不卡| 大陆偷拍与自拍| 日韩有码中文字幕| 50天的宝宝边吃奶边哭怎么回事| 97精品久久久久久久久久精品| 91成年电影在线观看| a级毛片在线看网站| 亚洲国产欧美在线一区| 丰满迷人的少妇在线观看| 亚洲精品国产av成人精品| 好男人电影高清在线观看| 99国产综合亚洲精品| 国产男人的电影天堂91| 国产成人欧美| 精品久久久久久久毛片微露脸 | 91字幕亚洲| 狠狠婷婷综合久久久久久88av| 国产亚洲午夜精品一区二区久久| 欧美 日韩 精品 国产| 日韩视频一区二区在线观看| 亚洲国产中文字幕在线视频| 久久人人97超碰香蕉20202| 性色av乱码一区二区三区2| av一本久久久久| 悠悠久久av| 久久中文字幕一级| 欧美少妇被猛烈插入视频| 蜜桃国产av成人99| 大型av网站在线播放| 91大片在线观看| 亚洲 国产 在线| 狠狠狠狠99中文字幕| 久久ye,这里只有精品| 老熟妇仑乱视频hdxx| 肉色欧美久久久久久久蜜桃| 伦理电影免费视频| 无遮挡黄片免费观看| 中文字幕另类日韩欧美亚洲嫩草| 又大又爽又粗| 热99re8久久精品国产| 亚洲国产中文字幕在线视频| 热re99久久国产66热| 国产精品国产三级国产专区5o| 国产精品香港三级国产av潘金莲| 性色av乱码一区二区三区2| 亚洲国产av影院在线观看| 亚洲 欧美一区二区三区| 性色av一级| 久久毛片免费看一区二区三区| 大型av网站在线播放| 可以免费在线观看a视频的电影网站| 午夜91福利影院| av天堂在线播放| 午夜福利在线免费观看网站| 热99re8久久精品国产| 精品少妇内射三级| 午夜福利乱码中文字幕| 极品人妻少妇av视频| 亚洲国产精品成人久久小说| 国产一区二区三区在线臀色熟女 | 中文字幕人妻丝袜一区二区| 久久久久视频综合| 亚洲全国av大片| 国产成人精品无人区| 天堂8中文在线网| 精品亚洲成a人片在线观看| 丝袜人妻中文字幕| 亚洲精品国产av成人精品| 久久中文看片网| 精品久久久久久电影网| 国产亚洲一区二区精品| 最近最新中文字幕大全免费视频| 脱女人内裤的视频| kizo精华| 欧美另类亚洲清纯唯美| 亚洲伊人色综图| 男女边摸边吃奶| 亚洲中文日韩欧美视频| 久久久精品免费免费高清| 在线观看免费午夜福利视频| 久热爱精品视频在线9| 五月开心婷婷网| 国产精品免费大片| 天天躁夜夜躁狠狠躁躁| 两性夫妻黄色片| 超碰成人久久| 久久久国产精品麻豆| 丰满人妻熟妇乱又伦精品不卡| 国产深夜福利视频在线观看| 亚洲色图 男人天堂 中文字幕| 夜夜骑夜夜射夜夜干| 午夜福利,免费看| 亚洲成国产人片在线观看| 热re99久久精品国产66热6| 国产伦人伦偷精品视频| 国产成人免费无遮挡视频| 久久精品aⅴ一区二区三区四区| 一进一出抽搐动态| 久久人人爽人人片av| 精品国产一区二区三区久久久樱花| 91大片在线观看| 一个人免费看片子| 午夜老司机福利片| 每晚都被弄得嗷嗷叫到高潮| 免费人妻精品一区二区三区视频| 老司机在亚洲福利影院| 狠狠婷婷综合久久久久久88av| 50天的宝宝边吃奶边哭怎么回事| 久久毛片免费看一区二区三区| 日本vs欧美在线观看视频| 亚洲五月色婷婷综合| 久久国产精品大桥未久av| 久久热在线av| 国产av一区二区精品久久| 极品人妻少妇av视频| 久久久精品区二区三区| 亚洲一区二区三区欧美精品| 国产成+人综合+亚洲专区| 99re6热这里在线精品视频| av网站在线播放免费| 嫩草影视91久久| 汤姆久久久久久久影院中文字幕| 欧美人与性动交α欧美软件| 亚洲人成77777在线视频| 亚洲va日本ⅴa欧美va伊人久久 | 久久久久网色| 成人国产一区最新在线观看| 国产亚洲精品久久久久5区| 精品一区在线观看国产| 国产成人欧美| 国产三级黄色录像| 一本色道久久久久久精品综合| 色老头精品视频在线观看| 精品乱码久久久久久99久播| 欧美日韩视频精品一区| 91大片在线观看| 不卡一级毛片| 一个人免费在线观看的高清视频 | 亚洲精品国产av成人精品| 国产精品香港三级国产av潘金莲| 麻豆av在线久日| 免费一级毛片在线播放高清视频 | 亚洲第一av免费看| 男女之事视频高清在线观看| 大片电影免费在线观看免费| 97人妻天天添夜夜摸| 国产xxxxx性猛交| 99re6热这里在线精品视频| 美女主播在线视频| 午夜久久久在线观看| 亚洲自偷自拍图片 自拍| 一级片'在线观看视频| 午夜两性在线视频| 欧美精品人与动牲交sv欧美| 又黄又粗又硬又大视频| 成人免费观看视频高清| 午夜免费鲁丝| 欧美黄色片欧美黄色片| 久久影院123| 美女午夜性视频免费| 日韩一卡2卡3卡4卡2021年| 美女视频免费永久观看网站| tocl精华| 亚洲成人免费电影在线观看| cao死你这个sao货| 无限看片的www在线观看| av有码第一页| 久久国产亚洲av麻豆专区| 久久亚洲国产成人精品v| 中文精品一卡2卡3卡4更新| 婷婷成人精品国产| 制服诱惑二区| 999精品在线视频| 色精品久久人妻99蜜桃| 日韩 亚洲 欧美在线| 视频区欧美日本亚洲| 性少妇av在线| 视频区图区小说| 免费在线观看黄色视频的| 久久国产亚洲av麻豆专区| av线在线观看网站| 男女免费视频国产| 国产一区二区 视频在线| 国产欧美日韩一区二区三区在线| 国产精品.久久久| 成年女人毛片免费观看观看9 | 极品少妇高潮喷水抽搐| 欧美精品一区二区大全| 国产成人av激情在线播放| 午夜精品国产一区二区电影| 俄罗斯特黄特色一大片| 黑人猛操日本美女一级片| 久久久久久人人人人人| 国产精品久久久久久人妻精品电影 | 国产日韩一区二区三区精品不卡| 他把我摸到了高潮在线观看 | 777久久人妻少妇嫩草av网站| 免费观看av网站的网址| 亚洲自偷自拍图片 自拍| 精品国产一区二区三区久久久樱花| 国产1区2区3区精品| 国产一区二区三区综合在线观看| 看免费av毛片| 成年人免费黄色播放视频| 国产成人av激情在线播放| a级片在线免费高清观看视频| 午夜成年电影在线免费观看| 在线永久观看黄色视频| 久久午夜综合久久蜜桃| 午夜福利,免费看| 欧美大码av| 国产精品久久久人人做人人爽| 老司机福利观看| 丝袜喷水一区| 老司机午夜福利在线观看视频 | 亚洲第一青青草原| 亚洲情色 制服丝袜| 国产又色又爽无遮挡免| 久久国产精品人妻蜜桃| a级片在线免费高清观看视频| 男女国产视频网站| 精品乱码久久久久久99久播| 搡老乐熟女国产| 久久毛片免费看一区二区三区| 午夜成年电影在线免费观看| 少妇被粗大的猛进出69影院| 国产人伦9x9x在线观看| 亚洲av美国av| 精品国内亚洲2022精品成人 | 自线自在国产av| 视频区图区小说| 嫩草影视91久久| 一区二区av电影网| videos熟女内射| 国产区一区二久久| 欧美在线黄色| www.熟女人妻精品国产| 欧美黄色片欧美黄色片| 国产av国产精品国产| 岛国毛片在线播放| 大片免费播放器 马上看| 国产成人精品在线电影| 夜夜骑夜夜射夜夜干| 精品国产乱码久久久久久男人| 一区二区三区乱码不卡18| 欧美激情高清一区二区三区| www日本在线高清视频| 亚洲欧美精品自产自拍| 水蜜桃什么品种好| 丝袜人妻中文字幕| 欧美亚洲日本最大视频资源| 日本一区二区免费在线视频| 午夜免费成人在线视频| 777米奇影视久久| 国产xxxxx性猛交| 国产男女内射视频| 国产又色又爽无遮挡免| 50天的宝宝边吃奶边哭怎么回事| 女人久久www免费人成看片| 免费一级毛片在线播放高清视频 | 十分钟在线观看高清视频www| 亚洲成人免费av在线播放| 久久久久视频综合| 精品人妻在线不人妻| 国产高清videossex| 久久中文看片网| 精品一区二区三区av网在线观看 | 久久99一区二区三区| 欧美激情久久久久久爽电影 | 久久国产精品人妻蜜桃| 日韩三级视频一区二区三区| 王馨瑶露胸无遮挡在线观看| 日本a在线网址| 色老头精品视频在线观看| 电影成人av| 日韩一区二区三区影片| 中文字幕人妻熟女乱码| 国产深夜福利视频在线观看| 成年女人毛片免费观看观看9 | 亚洲精品第二区| 免费人妻精品一区二区三区视频| 在线亚洲精品国产二区图片欧美| 亚洲av男天堂| 亚洲精品久久成人aⅴ小说| 欧美另类一区| 国产欧美日韩精品亚洲av| 丁香六月欧美| 欧美激情高清一区二区三区| 欧美日韩精品网址| 91国产中文字幕| 成人av一区二区三区在线看 | 在线观看免费视频网站a站| 国产高清videossex| 在线精品无人区一区二区三| 99久久99久久久精品蜜桃| 久久人妻福利社区极品人妻图片| 高清在线国产一区| 欧美日本中文国产一区发布| 99久久99久久久精品蜜桃| 性色av乱码一区二区三区2| 9热在线视频观看99| 国产成人影院久久av| 亚洲精品中文字幕在线视频| 日韩一卡2卡3卡4卡2021年| 午夜影院在线不卡| 国产精品一区二区在线不卡| 国产三级黄色录像| av不卡在线播放| a级毛片在线看网站| 啦啦啦啦在线视频资源| 菩萨蛮人人尽说江南好唐韦庄| 丝袜喷水一区| 正在播放国产对白刺激| 大香蕉久久网| 我要看黄色一级片免费的| 母亲3免费完整高清在线观看| 亚洲欧美激情在线| 国产成人av教育| 午夜福利视频在线观看免费| 欧美国产精品一级二级三级| 日本黄色日本黄色录像| 欧美乱码精品一区二区三区| 国产精品 国内视频| 另类精品久久| 日本撒尿小便嘘嘘汇集6| 青春草视频在线免费观看| 亚洲av电影在线观看一区二区三区| 啦啦啦 在线观看视频| 人妻人人澡人人爽人人| 亚洲精品国产一区二区精华液| 久久久久网色| 久久女婷五月综合色啪小说| 黄色怎么调成土黄色| kizo精华| 一级毛片电影观看| 国产淫语在线视频| 18禁裸乳无遮挡动漫免费视频| 丝瓜视频免费看黄片| 最黄视频免费看| 丝袜在线中文字幕| 成人手机av| 国产亚洲一区二区精品| 丝袜人妻中文字幕| 一二三四在线观看免费中文在| 亚洲国产精品一区二区三区在线| 日韩免费高清中文字幕av| 国产精品自产拍在线观看55亚洲 | 亚洲精品久久久久久婷婷小说| 亚洲欧美精品综合一区二区三区| 老鸭窝网址在线观看| 久久精品亚洲熟妇少妇任你| 久久人妻福利社区极品人妻图片| 一本久久精品| 少妇裸体淫交视频免费看高清 | 亚洲欧美色中文字幕在线| 高清av免费在线| 黄片小视频在线播放| 午夜福利视频精品| 母亲3免费完整高清在线观看| 国产一级毛片在线| 爱豆传媒免费全集在线观看| 国产伦理片在线播放av一区| 咕卡用的链子| 亚洲一区二区三区欧美精品| 黑人巨大精品欧美一区二区蜜桃| 高清av免费在线| 精品少妇黑人巨大在线播放| 免费在线观看视频国产中文字幕亚洲 | 丝袜美足系列| 久久久久久久久久久久大奶| 12—13女人毛片做爰片一| av欧美777| 亚洲精品久久午夜乱码| 成年美女黄网站色视频大全免费| 热re99久久精品国产66热6| 操美女的视频在线观看| 久久毛片免费看一区二区三区| 久久国产亚洲av麻豆专区| 在线看a的网站| 婷婷成人精品国产| 亚洲第一欧美日韩一区二区三区 | 国产在线视频一区二区| 国产精品熟女久久久久浪| 中文欧美无线码| 亚洲成人免费av在线播放| 成人18禁高潮啪啪吃奶动态图| 天天添夜夜摸| 一边摸一边做爽爽视频免费| 日本wwww免费看| 考比视频在线观看| 大片免费播放器 马上看| av超薄肉色丝袜交足视频| 久久久久久久大尺度免费视频| 各种免费的搞黄视频| 一本一本久久a久久精品综合妖精| 天天躁日日躁夜夜躁夜夜| 最新的欧美精品一区二区| 人人妻,人人澡人人爽秒播| 国产片内射在线| 色94色欧美一区二区| 一级a爱视频在线免费观看| 亚洲欧美成人综合另类久久久| 午夜福利乱码中文字幕| 日日夜夜操网爽| 大型av网站在线播放| 欧美日韩亚洲高清精品| 老司机福利观看| 日韩三级视频一区二区三区| 一边摸一边抽搐一进一出视频| 免费观看a级毛片全部| 天天影视国产精品| 亚洲 欧美一区二区三区| 大型av网站在线播放| 老熟妇仑乱视频hdxx| 欧美另类亚洲清纯唯美| 国产精品.久久久| 国产精品免费视频内射| 老司机午夜十八禁免费视频| 久久毛片免费看一区二区三区| 亚洲国产欧美网| 久久免费观看电影| 亚洲精品一区蜜桃| 丰满迷人的少妇在线观看| 肉色欧美久久久久久久蜜桃| 蜜桃国产av成人99| 大片电影免费在线观看免费| 女人久久www免费人成看片| 国产精品免费视频内射| 99热网站在线观看| 人妻一区二区av| 热99国产精品久久久久久7| 麻豆国产av国片精品| 日韩精品免费视频一区二区三区| 午夜福利视频精品| 久久久久国内视频| tube8黄色片| 我的亚洲天堂| 交换朋友夫妻互换小说| 老熟女久久久| 欧美久久黑人一区二区| 免费女性裸体啪啪无遮挡网站| 国产一级毛片在线| 日本猛色少妇xxxxx猛交久久| videosex国产| 麻豆乱淫一区二区| 老司机午夜福利在线观看视频 | 亚洲七黄色美女视频| 老熟妇仑乱视频hdxx| 久久久精品94久久精品| 最近中文字幕2019免费版| 亚洲精品中文字幕在线视频| 美女扒开内裤让男人捅视频| 亚洲av电影在线进入| 韩国高清视频一区二区三区| 色综合欧美亚洲国产小说| 亚洲av日韩精品久久久久久密| 国产视频一区二区在线看| 18禁裸乳无遮挡动漫免费视频| 日韩,欧美,国产一区二区三区| 美女国产高潮福利片在线看| 午夜老司机福利片| 亚洲av美国av| 下体分泌物呈黄色| 操美女的视频在线观看| 国产欧美日韩精品亚洲av| 日韩 亚洲 欧美在线| 久久久久久亚洲精品国产蜜桃av| 亚洲精品美女久久av网站| 亚洲欧美精品综合一区二区三区| 亚洲免费av在线视频| 久久久水蜜桃国产精品网| 午夜激情av网站| 国产97色在线日韩免费| 国产欧美亚洲国产| 侵犯人妻中文字幕一二三四区| 各种免费的搞黄视频| 女性生殖器流出的白浆| 一区二区三区激情视频| 人妻久久中文字幕网| 成人亚洲精品一区在线观看| 在线观看www视频免费| 午夜久久久在线观看| 欧美日韩亚洲综合一区二区三区_| 12—13女人毛片做爰片一| 亚洲欧洲精品一区二区精品久久久| 51午夜福利影视在线观看| 久久ye,这里只有精品|