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

    動(dòng)車組閘片磨損分析及壽命預(yù)測(cè)

    2024-01-03 04:50:44張輝李斌史立尹劍張生芳
    關(guān)鍵詞:閘片動(dòng)車組摩擦

    張輝,李斌,史立,尹劍,張生芳

    (1.天津鐵路信號(hào)有限責(zé)任公司,天津 300300;2.大連交通大學(xué) 機(jī)械工程學(xué)院,遼寧 大連 116028)

    隨著制動(dòng)技術(shù)的發(fā)展,盤式制動(dòng)器以其摩擦穩(wěn)定性高、制動(dòng)噪聲小等眾多優(yōu)點(diǎn),逐漸取代了摩擦制動(dòng)方式中的閘瓦制動(dòng)器,成為動(dòng)車制動(dòng)系統(tǒng)的重要形式[1]。盤式制動(dòng)器在工作過(guò)程中需要在短時(shí)間內(nèi)實(shí)現(xiàn)能量轉(zhuǎn)換,產(chǎn)生的摩擦熱會(huì)使制動(dòng)盤/片處于高溫高壓環(huán)境中,在高溫高壓作用下,制動(dòng)盤/片內(nèi)會(huì)產(chǎn)生熱應(yīng)力、壓應(yīng)力以及微弱的變形。同時(shí),由于摩擦作用,盤式制動(dòng)摩擦副會(huì)產(chǎn)生不同程度的磨損。磨損的累積會(huì)改變摩擦副工作過(guò)程中的接觸面積,進(jìn)而影響制動(dòng)穩(wěn)定性,降低動(dòng)車運(yùn)行安全性。在動(dòng)車組運(yùn)維過(guò)程中,針對(duì)閘片的疲勞壽命監(jiān)測(cè)多為預(yù)防性檢修,需要經(jīng)常對(duì)動(dòng)車組閘片進(jìn)行拆卸、檢查、維護(hù),因此需要對(duì)動(dòng)車組閘片制動(dòng)過(guò)程中的熱、力、磨損及其壽命預(yù)測(cè)進(jìn)行相關(guān)研究。

    王磊等[2]利用Ansys建立了動(dòng)車組用閘片循環(huán)對(duì)稱模型,分析了閘片幾何尺寸對(duì)熱流密度分配、對(duì)流散熱和輻射散熱的影響。夏毅敏等[3]基于傳熱理論和有限元方法,考慮溫度對(duì)材料特性的影響,利用Ansys分析了閘片溫度場(chǎng)和應(yīng)力場(chǎng)在制動(dòng)過(guò)程的分布以及變化規(guī)律。張生芳等[4]用間接耦合法將熱分析結(jié)果代入結(jié)構(gòu)場(chǎng),分析閘片的溫度及應(yīng)力變化。

    摩擦必然伴隨著磨損,磨損的計(jì)算方法種類繁多。Rajesh等[5]建立了圓盤與銷軸接觸分析的有限元模型,用3D表面到表面的體-地連接和旋轉(zhuǎn)類型表示接觸面的滑動(dòng)行為。張方宇等[6]模擬了摩擦片和制動(dòng)盤的磨合過(guò)程,證明了應(yīng)力-磨損耦合分析方法的可行性。郝強(qiáng)[7]使用微元法改進(jìn)了Archard磨粒磨損計(jì)算公式,利用有限元法計(jì)算探討了宏觀尺度下的閘片表面磨損分布規(guī)律。尹家寶等[8]提出一種求解摩擦塊磨損的數(shù)值方法,利用Abaqus結(jié)合磨損子程序?qū)崿F(xiàn)了磨損量及磨損量累計(jì)過(guò)程的模擬,最終通過(guò)Matlab二次開(kāi)發(fā)直接提取了磨損數(shù)值。趙子衡等[9]對(duì)傳統(tǒng)的Archard磨損模型進(jìn)行改進(jìn),模擬了某鼓式制動(dòng)器制動(dòng)過(guò)程,結(jié)果表明完全耦合法能夠更真實(shí)地模擬制動(dòng)過(guò)程的熱-力-磨損耦合。

    基于統(tǒng)計(jì)分析的壽命預(yù)測(cè)主要通過(guò)大量采集初始數(shù)據(jù)建立統(tǒng)計(jì)模型,從而實(shí)現(xiàn)壽命的預(yù)測(cè)。曲線擬合、灰色預(yù)測(cè)和神經(jīng)網(wǎng)絡(luò)等方法能夠很好地實(shí)現(xiàn)磨損預(yù)測(cè)。Zhang等[10]提出一種基于深度學(xué)習(xí)的壽命預(yù)測(cè)方法,組建了混合學(xué)習(xí)結(jié)構(gòu)網(wǎng)絡(luò),通過(guò)多項(xiàng)式回歸模型計(jì)算預(yù)測(cè)。Chang等[11]利用磨損模型和磨損進(jìn)程的遞歸函數(shù),通過(guò)原始算法推理了磨損隨時(shí)間的變化。Fan等[12]利用遺傳算法優(yōu)化了長(zhǎng)短時(shí)記憶網(wǎng)絡(luò),利用訓(xùn)練好的網(wǎng)絡(luò)進(jìn)行電機(jī)軸承壽命的預(yù)測(cè)。

    本文基于熱-力耦合有限元模型,通過(guò)Archard磨損理論進(jìn)行二次開(kāi)發(fā),實(shí)現(xiàn)了閘片磨損的有限元計(jì)算,模擬分析了閘片的磨損演變分布規(guī)律。采用單一變量法研究了制動(dòng)參數(shù)對(duì)閘片磨損的影響,并根據(jù)閘片磨損仿真,建立磨損與制動(dòng)參數(shù)對(duì)應(yīng)關(guān)系的數(shù)據(jù)庫(kù),通過(guò)BP算法調(diào)取數(shù)據(jù)建立磨損預(yù)測(cè)模型,建立動(dòng)車組閘片磨損量與行駛里程的對(duì)應(yīng)關(guān)系,最后設(shè)計(jì)了動(dòng)車組閘片壽命預(yù)測(cè)系統(tǒng)。

    1 制動(dòng)盤熱-力耦合磨損計(jì)算模型

    1.1 熱傳導(dǎo)初始條件及邊界條件

    根據(jù)傳熱學(xué)理論,摩擦熱會(huì)以熱傳導(dǎo)、熱對(duì)流以及熱輻射的形式傳遞到外部環(huán)境。

    (1)熱傳導(dǎo)模型

    將閘片看作無(wú)內(nèi)熱源的各向同性材料,根據(jù)傳熱理論,其熱傳導(dǎo)方程為:

    (1)

    式中:T為溫度,K;t為時(shí)間,s;ρ為材料密度,kg/m3;C為材料比熱容,J/(kg·K);λ為材料導(dǎo)熱系數(shù),W/(m·K)。

    動(dòng)車組盤式制動(dòng)器主要有熱傳導(dǎo)、熱對(duì)流及熱輻射,模型初始邊界條件為:在制動(dòng)初始時(shí)刻t=0時(shí),T(x,y,z)=T0,初始溫度T0取20 ℃。

    在閘片的所有換熱界面上:

    (2)

    在摩擦副接觸摩擦界面上:

    (3)

    式中:T0為環(huán)境溫度;ni為各界面的法向單位向量;z為摩擦界面的法向單位向量;αi為各界面的對(duì)流換熱系數(shù);αz為摩擦界面的對(duì)流換熱系數(shù);εi、εz為輻射換熱系數(shù),W/(m2·K);σ為斯蒂芬-波爾茲曼常數(shù);q為閘片吸收的熱流密度,W/m2。

    (2)對(duì)流換熱模型

    設(shè)強(qiáng)制對(duì)流換熱系數(shù)為αc,對(duì)于普朗特?cái)?shù)Pr>0.6的流體,當(dāng)雷諾數(shù)Rel

    (4)

    當(dāng)雷諾數(shù)Rel>Rexc時(shí),對(duì)流換熱系數(shù)為:

    (5)

    式中:Re為雷諾數(shù);Pr為空氣的普朗特?cái)?shù),取0.7;λ為空氣的熱導(dǎo)率,取λ=6.14×10-6W/(m·K);l為閘片的特征長(zhǎng)度。

    (3)輻射換熱系數(shù)

    制動(dòng)過(guò)程中存在輻射換熱現(xiàn)象。根據(jù)斯蒂芬-波爾茲曼定律,閘片與空氣之間的輻射換熱關(guān)系式為:

    (6)

    式中:ε為物體發(fā)射率;σ0為斯蒂芬-玻爾茲曼系數(shù),σ0=5.67×10-8W/(m2·k4)。

    1.2 動(dòng)車組盤式制動(dòng)器有限元模型的建立

    動(dòng)車組盤式制動(dòng)系統(tǒng)在工作過(guò)程中發(fā)生了大量的能量轉(zhuǎn)換,能量的轉(zhuǎn)換又包含多種形式,根據(jù)研究條件,提出如下假設(shè):

    a.閘片及制動(dòng)盤為各向同性材料,摩擦系數(shù)恒定為0.37。

    b.模擬制動(dòng)過(guò)程是勻減速。

    c.忽略摩擦過(guò)程中磨屑帶走的能量。

    d.制動(dòng)壓力在閘片背面均布分布。

    e.忽略制動(dòng)過(guò)程中環(huán)境溫度的變化,初始值設(shè)置為20 ℃。

    f.忽略振動(dòng)對(duì)結(jié)果的影響。

    制動(dòng)盤外半徑為320 mm,內(nèi)半徑為180 mm,制動(dòng)盤原型厚度為80 mm,制動(dòng)盤厚度為40 mm;閘片外半徑為310 mm,內(nèi)半徑為188 mm,厚度為20 mm,摩擦面積為28 500 mm2。

    (1)摩擦副材料參數(shù)

    摩擦副材料初始參數(shù)見(jiàn)表1,制動(dòng)盤材料參數(shù)隨溫度變化見(jiàn)表2。

    表1 摩擦副材料初始參數(shù)[4]

    (2)有限元模型前處理

    制動(dòng)盤和閘片采用六面體網(wǎng)格,制動(dòng)盤單元數(shù)量為2 400個(gè),閘片單元數(shù)量為540個(gè),采用減縮積分單元C3D8T,見(jiàn)圖1。

    圖1 制動(dòng)盤和閘片網(wǎng)格模型

    表面熱交換條件膜層散熱系數(shù)為0.001 W/(m·K),熱輻射發(fā)射率設(shè)置為0.2。圖2為有限元模型邊界條件示意圖,閘片四周表面受x、y方向的轉(zhuǎn)動(dòng)和位移約束,RP點(diǎn)約束在x、y方向的位移,施加轉(zhuǎn)速。制動(dòng)盤非接觸表面施加z方向的對(duì)稱約束。

    圖2 有限元模型邊界條件示意圖

    1.3 磨損計(jì)算方法

    Archard的磨損模型[13]形式如下:

    (7)

    式中:V為磨損體積;s為磨損位移;k為無(wú)量綱的磨損系數(shù);FN為法向載荷;H為磨損表面的硬度。

    將閘片接觸面劃分為N個(gè)無(wú)窮小的接觸面積,在單一無(wú)窮小的接觸面積ΔA上,時(shí)間增量為dt,則磨損模型演變的微分形式為:

    (8)

    式中:dV為磨損體積增量;ds為滑動(dòng)位移增量;FN為法向接觸力。

    假設(shè)ΔA在dt下的磨損厚度增量為dh,則dV=ΔAdh,代入式(7)得:

    (9)

    式中:FN/ΔA為無(wú)窮小接觸面積處的接觸應(yīng)力,令FN/ΔA為pc;k/H為有量綱磨損系數(shù),用kd表示,則式(9)表示為:

    dh=kdpcΔs

    (10)

    取時(shí)間增量為Δt,磨損厚度增量為Δh,相對(duì)滑動(dòng)位移增量為Δs,得:

    Δh=kdpcΔs

    (11)

    首先計(jì)算接觸應(yīng)力pc。當(dāng)Δs較小時(shí),認(rèn)為pc為不變量,求得對(duì)應(yīng)的Δh,根據(jù)h變化之后的模型,獲得接觸應(yīng)力pc。

    設(shè)閘片和制動(dòng)盤間接觸點(diǎn)i=1,…,N,歐拉積分步數(shù)j=1,…,M,第i個(gè)接觸點(diǎn)第j個(gè)積分步的接觸應(yīng)力為pc(i,j),磨損深度為Δh(i,j),j次磨損深度為h(i,j),制動(dòng)盤轉(zhuǎn)動(dòng)位移為Δs(i,j)。則:

    Δh(i,j)=kdpc(i,j)Δs(i,j)

    (12)

    h(i,j+1)=h(i,j)+Δh(i,j)

    (13)

    在ALE自適應(yīng)網(wǎng)格選項(xiàng)中設(shè)置磨損邊界為速度約束,將式(14)左右兩側(cè)同時(shí)除以第j個(gè)時(shí)間步長(zhǎng)Δt(j),得:

    (14)

    (15)

    2 閘片熱-力耦合及磨損結(jié)果分析

    2.1 閘片瞬態(tài)溫度場(chǎng)應(yīng)力場(chǎng)仿真結(jié)果分析

    本節(jié)模擬計(jì)算制動(dòng)初速度為120 km/h、制動(dòng)壓力為23 000 N時(shí)的制動(dòng)工況。圖3為制動(dòng)過(guò)程中閘片溫度場(chǎng)分布云圖。制動(dòng)開(kāi)始后,摩擦產(chǎn)熱遠(yuǎn)大于摩擦散熱,閘片表面的溫度不斷上升,隨著摩擦的持續(xù),摩擦表面溫度分布有明顯的變化。在1 s時(shí)摩擦面靠近閘片外徑處開(kāi)始出現(xiàn)溫度梯度;在13 s左右,由于轉(zhuǎn)速不斷降低,摩擦熱流的輸入與熱量的耗散達(dá)到平衡,達(dá)到了制動(dòng)過(guò)程中的最高溫度;之后摩擦產(chǎn)熱小于熱量的耗散,溫度開(kāi)始下降。整個(gè)制動(dòng)過(guò)程中,閘片表面最高溫度達(dá)到了209.4 ℃,制動(dòng)結(jié)束時(shí)最高溫度為122.6 ℃,閘片的最高溫度出現(xiàn)在摩擦半徑273 mm處。

    圖3 制動(dòng)過(guò)程中閘片溫度分布云圖

    圖4為制動(dòng)過(guò)程中閘片應(yīng)力分布云圖。在制動(dòng)過(guò)程中,由于摩擦作用,閘片摩擦表面累積了大量的熱量,在摩擦表面存在相應(yīng)的溫度梯度,高溫使得閘片的變形向徑向擴(kuò)展,閘片摩擦區(qū)域同時(shí)產(chǎn)生很大的熱應(yīng)力。制動(dòng)初期,由于摩擦入口先發(fā)生滑動(dòng),閘片摩擦入口區(qū)域產(chǎn)生局部熱應(yīng)力;摩擦制動(dòng)持續(xù)進(jìn)行之后,摩擦出口區(qū)域應(yīng)力逐漸大于其余位置,閘片徑向的溫度梯度同時(shí)影響閘片應(yīng)力分布,在摩擦半徑273 mm附近出現(xiàn)了最大應(yīng)力。從制動(dòng)開(kāi)始13 s后閘片摩擦表面的徑向應(yīng)力分布呈明顯的梯度分布,由于應(yīng)力集中的存在,周向應(yīng)力的分布存在較大的差異。

    圖4 制動(dòng)過(guò)程中閘片應(yīng)力分布云圖

    2.2 閘片磨損仿真結(jié)果分析

    本節(jié)將基于上述有限元模型,設(shè)置制動(dòng)初速度為160 km/h、制動(dòng)壓力為14 000 N、仿真時(shí)長(zhǎng)為10 s,不同時(shí)刻閘片磨損分布云圖見(jiàn)圖5。閘片左側(cè)為摩擦入口區(qū)域,右側(cè)為摩擦出口區(qū)域。在制動(dòng)初期,壓力大、速度快,磨損值增加迅速,閘片中間部分磨損最小,隨后在溫度以及應(yīng)力共同作用下,閘片的上下邊緣部分磨損程度最小。制動(dòng)初期摩擦入口區(qū)域磨損大,由于摩擦方向向右,隨著相對(duì)運(yùn)動(dòng)的持續(xù),逐漸演變成摩擦出口區(qū)域的磨損大。

    圖5 不同時(shí)刻閘片磨損分布云圖

    選擇閘片不同周向以及徑向節(jié)點(diǎn),周向節(jié)點(diǎn)從左到右編號(hào)1~19,徑向節(jié)點(diǎn)從下到上編號(hào)1~11,見(jiàn)圖6。閘片表面磨損規(guī)律見(jiàn)圖7, 閘片徑向磨損呈

    圖6 閘片表面磨損分析節(jié)點(diǎn)圖

    (a) 周向磨損

    (b) 徑向磨損圖7 閘片表面磨損規(guī)律

    現(xiàn)先增大后減小的趨勢(shì),在制動(dòng)過(guò)程中閘片內(nèi)徑的線速度相對(duì)較小,因此節(jié)點(diǎn)1磨損值最小,隨著各節(jié)點(diǎn)摩擦半徑的增大,磨損程度相比于內(nèi)側(cè)均有所增大。制動(dòng)過(guò)程中溫度及應(yīng)力最大值區(qū)域位于摩擦半徑273 mm(節(jié)點(diǎn)8)附近,溫度與應(yīng)力大小與磨損深度成正比,因此節(jié)點(diǎn)8處磨損最大,節(jié)點(diǎn)8向兩側(cè)呈現(xiàn)磨損梯度值。

    2.3 制動(dòng)初速度對(duì)閘片磨損的影響

    將初始?jí)毫υO(shè)定為23 kN,隨著制動(dòng)初速度的增大,閘片表面磨損深度不斷增加,但閘片摩擦表面磨損分布趨勢(shì)無(wú)明顯變化。速度從80到180 km/h,最大磨損深度為1.40×10-2、1.69×10-2、2.01×10-2、2.27×10-2、2.53×10-2、2.80×10-2mm。制動(dòng)過(guò)程中,初始速度直接決定摩擦副之間的相對(duì)滑動(dòng)速度,在相同時(shí)間、相同制動(dòng)壓力下,初速度越大,摩擦副表面相對(duì)滑動(dòng)距離越大,導(dǎo)致閘片表面磨損加劇,不同制動(dòng)初速度下閘片磨損曲線見(jiàn)圖8。

    圖8 不同制動(dòng)初速度下閘片磨損曲線

    從圖8可以看出,在制動(dòng)初期,不同轉(zhuǎn)速下最大磨損量差值較小,隨著制動(dòng)時(shí)間的增加,其差值逐漸增大,且差值存在一定的差異。以80 km/h為起點(diǎn),20 km/h為一個(gè)跨度,速度每增加20 km/h,磨損深度分別增加2.91×10-3、3.19×10-3、2.66×10-3、2.52×10-3、2.78×10-3mm。如圖9所示,隨著速度的增加,磨損深度差值的變化呈現(xiàn)先上升后下降再上升的趨勢(shì),但磨損深度的變化率逐漸減小,后期趨于平穩(wěn)。速度增大,磨損深度變化較大,速度達(dá)到120 km/h以后,摩擦副溫度較大。同時(shí),由于閘片節(jié)點(diǎn)磨損,摩擦副接觸面積發(fā)生變化,摩擦力大小發(fā)生一定的變化,導(dǎo)致磨損深度增加量減小,因此磨損變化率逐漸減小。由于速度不斷上升,摩擦副相對(duì)滑動(dòng)速度逐漸變大,滑動(dòng)速度的增大將會(huì)加劇閘片接觸表面的磨損,磨損深度進(jìn)一步增大。

    (a) 磨損深度差值

    (b) 磨損變化率圖9 磨損變化曲線

    2.4 制動(dòng)壓力對(duì)閘片磨損的影響

    將初始速度設(shè)定為120 km/h,閘片摩擦表面磨損分布趨勢(shì)隨著制動(dòng)壓力的增大無(wú)明顯變化。圖10為不同制動(dòng)壓力下閘片磨損曲線,不同壓力下,磨損深度的變化趨勢(shì)相同,摩擦面摩擦力隨著制動(dòng)壓力的增大而增大,進(jìn)一步增大摩擦表面接觸應(yīng)力,迫使閘片表面發(fā)生變形。隨著制動(dòng)壓力的增大而增大,擠壓摩擦力和剪切力也增大,進(jìn)一步加劇了磨損程度。

    圖10 不同制動(dòng)壓力下閘片磨損曲線

    在10 s時(shí)間內(nèi),制動(dòng)壓力為17 kN時(shí)最大磨損深度為1.52×10-2mm,隨后壓力每增加3 kN,最大磨損深度分別增加2.313 8×10-3、2.580 6×10-3、2.231 4×10-3、2.332 9×10-3、2.253 7×10-3mm,利用差值計(jì)算磨損變化率繪圖,見(jiàn)圖11。隨著制動(dòng)壓力的增大,摩擦力線性增大,摩擦制動(dòng)使得閘片溫升變大,導(dǎo)致閘片發(fā)生熱變形,在摩擦副相對(duì)滑動(dòng)的同時(shí),由于熱變形的影響,摩擦接觸面積一直變化,因此磨損差值會(huì)呈現(xiàn)出上下波動(dòng),但其磨損變化率呈現(xiàn)逐漸降低的趨勢(shì)。

    (a) 磨損深度差值

    (b) 磨損變化率圖11 磨損變化曲線

    3 動(dòng)車組閘片壽命預(yù)測(cè)系統(tǒng)設(shè)計(jì)

    3.1 動(dòng)車組閘片磨損預(yù)測(cè)模型的建立

    數(shù)據(jù)預(yù)處理保證數(shù)據(jù)庫(kù)完整性的同時(shí)還可以減小數(shù)據(jù)庫(kù)的量,數(shù)據(jù)預(yù)處理包含數(shù)據(jù)去噪、特征選取以及特征數(shù)據(jù)處理3個(gè)過(guò)程。

    (1)輸入輸出層結(jié)構(gòu)

    以閘片各節(jié)點(diǎn)的磨損值為輸出參數(shù),輸入層為制動(dòng)初速度、制動(dòng)壓力、制動(dòng)時(shí)間。輸出層神經(jīng)元包括閘片209個(gè)節(jié)點(diǎn)的磨損仿真數(shù)值。

    (2)隱藏層結(jié)構(gòu)

    本文采用單隱層、3層輸入、多輸出的BP神經(jīng)網(wǎng)絡(luò)建立預(yù)測(cè)模型。目前,隱層神經(jīng)元數(shù)量有幾種相關(guān)的經(jīng)驗(yàn)公式可供參考,本文參考公式(16):

    (16)

    式中:m、n為輸入層、輸出層神經(jīng)元數(shù)量;a∈[1,10]。

    確定該網(wǎng)絡(luò)結(jié)構(gòu)的隱藏層神經(jīng)元數(shù)量個(gè)數(shù)為15~24,該模型選擇18個(gè)隱層神經(jīng)元。

    (3)激勵(lì)函數(shù)

    選擇S型正切函數(shù)tansig作為隱藏層神經(jīng)元的激勵(lì)函數(shù),選取S型對(duì)數(shù)函數(shù)logsig作為輸出層神經(jīng)元的激勵(lì)函數(shù)。

    (4)學(xué)習(xí)率

    學(xué)習(xí)率大小直接影響學(xué)習(xí)精度、學(xué)習(xí)速度。通常情況下,學(xué)習(xí)率的范圍大概取0.01~0.70即可,本模型學(xué)習(xí)率為0.01。

    選擇3種制動(dòng)壓力,分別為14、23、32 kN;初速度選擇范圍為60~200 km/h,以20 km/h為一個(gè)跨度;仿真時(shí)長(zhǎng)20 s,以2 s為一個(gè)時(shí)間跨度選擇數(shù)據(jù)進(jìn)行訓(xùn)練,共提取出230組數(shù)據(jù)進(jìn)行模型的訓(xùn)練、驗(yàn)證以及測(cè)試。

    表3 為Progress區(qū)域訓(xùn)練結(jié)果統(tǒng)計(jì)表,Epoch表示實(shí)際訓(xùn)練次數(shù)(19次);Time表示訓(xùn)練總時(shí)間(31 min 9 s);Performance為性能指標(biāo),該模型訓(xùn)練均方誤差為1.37×10-7;Gradient為梯度,進(jìn)度條中顯示當(dāng)前的梯度值為8.37×10-8,右側(cè)設(shè)定的梯度值為1×10-7,梯度值達(dá)到設(shè)定值則停止訓(xùn)練。

    表3 網(wǎng)絡(luò)訓(xùn)練結(jié)果數(shù)據(jù)統(tǒng)計(jì)表

    3.2 磨損量-行駛里程對(duì)應(yīng)關(guān)系的建立

    目前復(fù)興號(hào)動(dòng)車組閘片厚度原型為20 mm,其余動(dòng)車組閘片厚度原型為17 mm。其中毫米磨耗行駛里程=壽命萬(wàn)公里/(原型厚度-磨耗到限)。根據(jù)最大磨損深度計(jì)算失效磨耗量,以CRH380BL上的博科瑞B(yǎng)KR380FA為例,動(dòng)車組閘片最大壽命為79.3 萬(wàn)km(取79 萬(wàn)km),拖車為44.2 萬(wàn)km(取44 萬(wàn)km)。閘片原型厚度為17 mm,允許磨耗到限厚度最薄為5 mm。對(duì)上述數(shù)據(jù)進(jìn)行計(jì)算,閘片每磨耗1 mm,可行駛里程動(dòng)車取n1=658 33 km,拖車取n2=36 666 km。根據(jù)動(dòng)車與拖車行駛里程和單位磨損量的關(guān)系,當(dāng)動(dòng)車組閘片磨損x時(shí),計(jì)算得到拖車閘片磨損x2=xn1/n2≈1.8x。綜合以上選擇數(shù)據(jù),其壽命計(jì)算為:

    y1=(12-x)×65 833

    (17)

    y2=(12-1.8x)×36 666

    (18)

    式中:y1、y2分別為動(dòng)車組閘片的剩余行駛里程和拖車閘片的剩余行駛里程;x為動(dòng)車組閘片磨損厚度。

    3.3 系統(tǒng)軟件界面設(shè)計(jì)結(jié)果及應(yīng)用

    壽命預(yù)測(cè)計(jì)算App基于Matlab中的圖形用戶界面設(shè)計(jì),在用戶與計(jì)算機(jī)之間建立交互操作。使用面板、圖標(biāo)按鈕等功能框進(jìn)行設(shè)計(jì)界面布局,通過(guò)編輯回調(diào)函數(shù),邏輯選擇的switch分支結(jié)構(gòu)見(jiàn)圖12。

    圖12 switch分支結(jié)構(gòu)

    圖13為預(yù)測(cè)系統(tǒng)使用案例??梢钥闯?當(dāng)輸入制動(dòng)初速度為160 km/h、制動(dòng)壓力為32 kN、制動(dòng)時(shí)間為20 s時(shí),磨損預(yù)測(cè)模型磨損深度預(yù)測(cè)結(jié)果最大值為0.030 76 mm,期望輸出最大目標(biāo)值為0.030 741 mm,誤差值為1.9×10-5mm;磨損深度預(yù)測(cè)結(jié)果平均值為0.024 52 mm,期望輸出目標(biāo)的平均值為0.024 485 mm,誤差值為3.5×10-5mm,可以看出磨損深度預(yù)測(cè)精度較高。隨后針對(duì)壽命計(jì)算模塊進(jìn)行型號(hào)選擇,閘片型號(hào)選擇DC20,動(dòng)車型號(hào)選擇CRH380AL,選擇型號(hào)之后動(dòng)車、拖車磨損極限彈出初始值為12 mm,此時(shí)計(jì)算按鈕的回調(diào)函數(shù)根據(jù)switch語(yǔ)句選擇出CRH380AL車型上DC20閘片對(duì)應(yīng)的壽命計(jì)算方式,隨后根據(jù)輸入磨損量計(jì)算剩余行駛里程以及剩余磨損厚度,并自動(dòng)更新迭代至動(dòng)車、拖車磨損極限。

    圖13 預(yù)測(cè)系統(tǒng)使用案例

    4 結(jié)論

    (1) 本文基于熱-力耦合仿真模型,結(jié)合自適應(yīng)網(wǎng)格技術(shù)進(jìn)行二次開(kāi)發(fā),模擬計(jì)算了閘片磨損過(guò)程,分析了制動(dòng)初始參數(shù)對(duì)于閘片磨損的影響。結(jié)果表明:該模型磨損仿真結(jié)果與試驗(yàn)閘片磨損分布趨勢(shì)相同,閘片的摩擦出口區(qū)域出現(xiàn)了較嚴(yán)重的磨損,磨損在徑向呈現(xiàn)梯度分布。制動(dòng)初速度從80 km/h增大到180 km/h時(shí),最大磨損深度增加了1.40×10-2mm。制動(dòng)壓力從17 kN增大到32 kN時(shí),最大磨損深度增加了0.94×10-2mm。

    (2)本文主要分析動(dòng)車組閘片壽命預(yù)測(cè)的實(shí)現(xiàn)途徑及實(shí)現(xiàn)過(guò)程,根據(jù)閘片磨損仿真結(jié)果,建立磨損與制動(dòng)參數(shù)對(duì)應(yīng)關(guān)系的數(shù)據(jù)庫(kù),基于BP神經(jīng)網(wǎng)絡(luò)模型,運(yùn)用數(shù)據(jù)庫(kù)數(shù)據(jù)對(duì)預(yù)測(cè)模型進(jìn)行了訓(xùn)練、驗(yàn)證和測(cè)試,建立了閘片磨損預(yù)測(cè)模型。根據(jù)動(dòng)車組運(yùn)行數(shù)據(jù)建立了閘片磨損量-行駛里程的對(duì)應(yīng)關(guān)系。通過(guò)編程語(yǔ)言調(diào)用訓(xùn)練好的模型,編寫了各車型閘片的壽命計(jì)算邏輯,建立了集磨損預(yù)測(cè)、壽命計(jì)算為一體的動(dòng)車組閘片壽命預(yù)測(cè)系統(tǒng)。

    猜你喜歡
    閘片動(dòng)車組摩擦
    動(dòng)車組閘片壽命預(yù)測(cè)研究
    干摩擦和濕摩擦的區(qū)別
    廣州地鐵三號(hào)線制動(dòng)閘片斷裂問(wèn)題分析及優(yōu)化
    “95后”動(dòng)車組女司機(jī)的首個(gè)春運(yùn)
    海峽姐妹(2020年2期)2020-03-03 13:36:34
    神奇的摩擦起電
    條分縷析 摩擦真相
    動(dòng)車組BTM帶內(nèi)干擾的排查與整治
    降低動(dòng)車組閘片脫落風(fēng)險(xiǎn)的有效措施及其運(yùn)用效果
    解讀摩擦起電
    CRH3型動(dòng)車組輪對(duì)壓裝曲線研究
    中文字幕人妻熟人妻熟丝袜美 | 亚洲精品乱码久久久v下载方式 | 99精品欧美一区二区三区四区| 成人特级av手机在线观看| 搡女人真爽免费视频火全软件 | 超碰av人人做人人爽久久 | 午夜福利成人在线免费观看| 欧美丝袜亚洲另类 | 人人妻人人看人人澡| 两人在一起打扑克的视频| 一本一本综合久久| 国内久久婷婷六月综合欲色啪| 91九色精品人成在线观看| 国产极品精品免费视频能看的| 精品国产三级普通话版| 国产精品久久久久久人妻精品电影| 在线观看日韩欧美| 啦啦啦韩国在线观看视频| 国产精品国产高清国产av| 国产男靠女视频免费网站| 国产欧美日韩精品亚洲av| 国产精品久久视频播放| 日本在线视频免费播放| 国语自产精品视频在线第100页| 日本撒尿小便嘘嘘汇集6| 热99re8久久精品国产| 精品久久久久久久久久免费视频| 99视频精品全部免费 在线| 亚洲精品日韩av片在线观看 | 色尼玛亚洲综合影院| 内射极品少妇av片p| 欧美+亚洲+日韩+国产| www.999成人在线观看| 国模一区二区三区四区视频| 桃色一区二区三区在线观看| 午夜福利成人在线免费观看| 久久久精品大字幕| 伊人久久精品亚洲午夜| 又黄又粗又硬又大视频| 99久国产av精品| 久久草成人影院| 久久久国产成人精品二区| 成年女人毛片免费观看观看9| 国产精品久久电影中文字幕| 欧美激情久久久久久爽电影| 亚洲美女视频黄频| 久久久久国内视频| 中文在线观看免费www的网站| 亚洲av五月六月丁香网| 搡女人真爽免费视频火全软件 | 欧美乱码精品一区二区三区| 村上凉子中文字幕在线| 偷拍熟女少妇极品色| 免费在线观看亚洲国产| 欧美乱色亚洲激情| 国产不卡一卡二| 岛国在线免费视频观看| 日韩欧美在线二视频| 少妇高潮的动态图| 久久亚洲精品不卡| 亚洲五月天丁香| 内地一区二区视频在线| 国产毛片a区久久久久| 亚洲国产精品999在线| 日韩欧美免费精品| 亚洲精品456在线播放app | 成人18禁在线播放| 国产在视频线在精品| 美女被艹到高潮喷水动态| av福利片在线观看| 午夜福利18| 国产成人福利小说| 18美女黄网站色大片免费观看| 非洲黑人性xxxx精品又粗又长| 中出人妻视频一区二区| 国产熟女xx| 国产精品亚洲一级av第二区| 精华霜和精华液先用哪个| 欧美性猛交黑人性爽| 青草久久国产| 成年免费大片在线观看| 免费看十八禁软件| 国产成人aa在线观看| 亚洲久久久久久中文字幕| 国产色婷婷99| 久久精品国产亚洲av涩爱 | 成熟少妇高潮喷水视频| 国内精品久久久久精免费| 少妇的逼水好多| 欧美最黄视频在线播放免费| 国产毛片a区久久久久| 18禁黄网站禁片免费观看直播| 90打野战视频偷拍视频| 久久精品91无色码中文字幕| 少妇的丰满在线观看| 在线国产一区二区在线| 波野结衣二区三区在线 | 国内毛片毛片毛片毛片毛片| 欧美日韩福利视频一区二区| 中文在线观看免费www的网站| 波野结衣二区三区在线 | 51午夜福利影视在线观看| 午夜激情福利司机影院| 日本三级黄在线观看| 国产成+人综合+亚洲专区| 黄色成人免费大全| 日本黄色视频三级网站网址| 亚洲专区中文字幕在线| 午夜福利欧美成人| 亚洲,欧美精品.| 欧美激情在线99| 国产麻豆成人av免费视频| 特大巨黑吊av在线直播| 国产高清videossex| 久久久久久久亚洲中文字幕 | 女同久久另类99精品国产91| 日本免费一区二区三区高清不卡| 国产伦精品一区二区三区视频9 | 日韩亚洲欧美综合| 日本在线视频免费播放| 一个人看视频在线观看www免费 | 久久久久久久精品吃奶| 无人区码免费观看不卡| 亚洲精品乱码久久久v下载方式 | 国产精品一区二区三区四区久久| 国产三级黄色录像| 男插女下体视频免费在线播放| 国产中年淑女户外野战色| 两个人看的免费小视频| 淫秽高清视频在线观看| 国产精品99久久99久久久不卡| 久久精品91蜜桃| 国产精品香港三级国产av潘金莲| 俄罗斯特黄特色一大片| 国产亚洲精品久久久com| 国产精品美女特级片免费视频播放器| 淫妇啪啪啪对白视频| 两个人的视频大全免费| 听说在线观看完整版免费高清| 51国产日韩欧美| 国内精品久久久久精免费| av福利片在线观看| 欧美日本亚洲视频在线播放| 97超视频在线观看视频| 狠狠狠狠99中文字幕| 国产成人aa在线观看| 夜夜躁狠狠躁天天躁| 国产高清videossex| 老司机深夜福利视频在线观看| 国产乱人视频| 美女黄网站色视频| 91久久精品国产一区二区成人 | 亚洲av成人精品一区久久| 久久久久久人人人人人| 亚洲av第一区精品v没综合| 国产久久久一区二区三区| 全区人妻精品视频| 亚洲av电影在线进入| 熟女少妇亚洲综合色aaa.| 婷婷精品国产亚洲av| 久久99热这里只有精品18| 日本一二三区视频观看| 国产在视频线在精品| netflix在线观看网站| 亚洲18禁久久av| 日本免费一区二区三区高清不卡| 搞女人的毛片| 亚洲精品一卡2卡三卡4卡5卡| 女同久久另类99精品国产91| 国产高清视频在线播放一区| 国产精品av视频在线免费观看| 男女之事视频高清在线观看| 免费在线观看成人毛片| 久久人妻av系列| 在线天堂最新版资源| 天堂动漫精品| 亚洲最大成人手机在线| 亚洲av电影不卡..在线观看| 亚洲国产高清在线一区二区三| 午夜福利免费观看在线| 精品日产1卡2卡| 中文亚洲av片在线观看爽| 国产成人福利小说| 久久久久精品国产欧美久久久| 久久精品人妻少妇| 美女cb高潮喷水在线观看| 91在线观看av| 精品乱码久久久久久99久播| 欧美色欧美亚洲另类二区| 三级男女做爰猛烈吃奶摸视频| 免费高清视频大片| 高清在线国产一区| 不卡一级毛片| 亚洲av熟女| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国产探花极品一区二区| 国产高清视频在线观看网站| 欧美极品一区二区三区四区| 亚洲 国产 在线| 男女视频在线观看网站免费| 亚洲av电影不卡..在线观看| 国产乱人伦免费视频| 最好的美女福利视频网| 好男人电影高清在线观看| 国产av在哪里看| 色综合站精品国产| 午夜激情福利司机影院| 中文字幕av在线有码专区| 成人国产一区最新在线观看| or卡值多少钱| 午夜久久久久精精品| 伊人久久大香线蕉亚洲五| 嫩草影院入口| 看黄色毛片网站| 一区福利在线观看| 18禁黄网站禁片午夜丰满| 亚洲av二区三区四区| 色综合站精品国产| 国产在线精品亚洲第一网站| 男女做爰动态图高潮gif福利片| 亚洲人与动物交配视频| 在线天堂最新版资源| 亚洲av一区综合| 国产av一区在线观看免费| 久久精品人妻少妇| 一本一本综合久久| 欧美日韩综合久久久久久 | 悠悠久久av| 日韩av在线大香蕉| 国产成+人综合+亚洲专区| 五月玫瑰六月丁香| 中文字幕人妻熟人妻熟丝袜美 | 日本三级黄在线观看| 波多野结衣高清无吗| 国产精品久久久久久人妻精品电影| av专区在线播放| 国内揄拍国产精品人妻在线| 国产视频内射| 国产高清视频在线播放一区| 男插女下体视频免费在线播放| 国产午夜精品论理片| 国产91精品成人一区二区三区| 别揉我奶头~嗯~啊~动态视频| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲精品美女久久久久99蜜臀| 亚洲久久久久久中文字幕| 日韩欧美国产一区二区入口| 久久午夜亚洲精品久久| 99视频精品全部免费 在线| 免费观看人在逋| 欧美最黄视频在线播放免费| 深夜精品福利| 在线播放国产精品三级| 精华霜和精华液先用哪个| 精品熟女少妇八av免费久了| 动漫黄色视频在线观看| 最近最新免费中文字幕在线| 天天一区二区日本电影三级| 免费大片18禁| 久久久久久久久中文| 丰满的人妻完整版| 色av中文字幕| 村上凉子中文字幕在线| 很黄的视频免费| 午夜两性在线视频| 午夜福利在线观看吧| 欧美乱码精品一区二区三区| 欧美成人一区二区免费高清观看| 亚洲va日本ⅴa欧美va伊人久久| 亚洲内射少妇av| 国产精品 欧美亚洲| 亚洲真实伦在线观看| 禁无遮挡网站| 亚洲中文字幕日韩| 欧美成人一区二区免费高清观看| 亚洲国产日韩欧美精品在线观看 | 9191精品国产免费久久| 久久久久久久午夜电影| 中文字幕久久专区| 97超视频在线观看视频| 色尼玛亚洲综合影院| 精品久久久久久,| 别揉我奶头~嗯~啊~动态视频| 麻豆成人av在线观看| 俄罗斯特黄特色一大片| 久久精品国产综合久久久| 精品乱码久久久久久99久播| 熟女少妇亚洲综合色aaa.| 最好的美女福利视频网| 日本撒尿小便嘘嘘汇集6| 国产高清视频在线播放一区| 99在线视频只有这里精品首页| 成人精品一区二区免费| 在线免费观看的www视频| 99精品久久久久人妻精品| 午夜两性在线视频| 在线观看美女被高潮喷水网站 | 久久国产精品影院| 黄片小视频在线播放| 女警被强在线播放| 又黄又粗又硬又大视频| 久久中文看片网| 亚洲国产色片| 国产亚洲精品综合一区在线观看| 日本五十路高清| 69av精品久久久久久| 夜夜躁狠狠躁天天躁| 国产男靠女视频免费网站| 色精品久久人妻99蜜桃| 国产激情偷乱视频一区二区| 少妇人妻一区二区三区视频| 亚洲成a人片在线一区二区| 12—13女人毛片做爰片一| 国产伦人伦偷精品视频| 国产av一区在线观看免费| 伊人久久大香线蕉亚洲五| 午夜a级毛片| 少妇丰满av| 日本 欧美在线| 国内久久婷婷六月综合欲色啪| 久久亚洲真实| 蜜桃久久精品国产亚洲av| 久久精品国产综合久久久| 久久精品夜夜夜夜夜久久蜜豆| 久久性视频一级片| 午夜视频国产福利| 日韩欧美国产一区二区入口| 午夜激情欧美在线| 99国产极品粉嫩在线观看| 亚洲成人久久爱视频| 欧美三级亚洲精品| 男人的好看免费观看在线视频| 黄色成人免费大全| 在线观看午夜福利视频| 久久国产精品影院| 日本 欧美在线| 99久久九九国产精品国产免费| 搡老妇女老女人老熟妇| 熟妇人妻久久中文字幕3abv| 免费大片18禁| 亚洲va日本ⅴa欧美va伊人久久| 国产亚洲欧美在线一区二区| 99热这里只有精品一区| 99视频精品全部免费 在线| 久久香蕉精品热| 国产精品乱码一区二三区的特点| 三级男女做爰猛烈吃奶摸视频| 国产欧美日韩精品一区二区| 国产精品98久久久久久宅男小说| 哪里可以看免费的av片| 女人被狂操c到高潮| 午夜两性在线视频| 淫秽高清视频在线观看| 少妇的丰满在线观看| 亚洲欧美日韩高清在线视频| 18禁裸乳无遮挡免费网站照片| 久久草成人影院| 欧美成人免费av一区二区三区| 精品国产亚洲在线| 天天添夜夜摸| 九色国产91popny在线| 中文字幕av成人在线电影| 久久亚洲真实| 在线国产一区二区在线| 国产欧美日韩精品亚洲av| 欧美一级a爱片免费观看看| 宅男免费午夜| 一级毛片女人18水好多| 久久久久久久午夜电影| 在线播放无遮挡| 麻豆国产97在线/欧美| 啦啦啦韩国在线观看视频| 亚洲精品粉嫩美女一区| 波多野结衣巨乳人妻| 久久久国产精品麻豆| 国产精品野战在线观看| 黄色视频,在线免费观看| 亚洲无线观看免费| 成人午夜高清在线视频| 成人国产综合亚洲| а√天堂www在线а√下载| 内射极品少妇av片p| 噜噜噜噜噜久久久久久91| 久久精品国产综合久久久| 18禁在线播放成人免费| 久久久久国产精品人妻aⅴ院| 一夜夜www| 欧美日韩国产亚洲二区| 丝袜美腿在线中文| 99视频精品全部免费 在线| 男人的好看免费观看在线视频| 精品国产超薄肉色丝袜足j| 亚洲av免费高清在线观看| av在线蜜桃| 禁无遮挡网站| 日韩免费av在线播放| 午夜精品在线福利| 亚洲精品亚洲一区二区| 国产单亲对白刺激| 一a级毛片在线观看| 亚洲av电影不卡..在线观看| www.色视频.com| 久久久久精品国产欧美久久久| 色吧在线观看| 国产色婷婷99| 久久精品综合一区二区三区| 欧美乱色亚洲激情| 嫩草影院入口| 999久久久精品免费观看国产| 麻豆成人午夜福利视频| 男人舔奶头视频| 国产69精品久久久久777片| 久久国产乱子伦精品免费另类| 精品午夜福利视频在线观看一区| 久久久久久久久大av| 国产在视频线在精品| 88av欧美| 日韩欧美在线乱码| 最近最新中文字幕大全免费视频| 99视频精品全部免费 在线| 桃红色精品国产亚洲av| 精品人妻偷拍中文字幕| 97超级碰碰碰精品色视频在线观看| 少妇人妻一区二区三区视频| svipshipincom国产片| 丰满人妻一区二区三区视频av | 国产爱豆传媒在线观看| 欧美丝袜亚洲另类 | 97碰自拍视频| 中出人妻视频一区二区| 他把我摸到了高潮在线观看| 日日夜夜操网爽| or卡值多少钱| 成人国产综合亚洲| 亚洲最大成人手机在线| 亚洲第一欧美日韩一区二区三区| 丰满的人妻完整版| 窝窝影院91人妻| 黄色视频,在线免费观看| 亚洲欧美日韩东京热| www.www免费av| 日本成人三级电影网站| 久久久色成人| 成人国产一区最新在线观看| www国产在线视频色| 男人舔女人下体高潮全视频| 国产色爽女视频免费观看| 亚洲精品国产精品久久久不卡| 在线国产一区二区在线| 制服丝袜大香蕉在线| 国语自产精品视频在线第100页| 午夜福利视频1000在线观看| 国产精品综合久久久久久久免费| 一级黄片播放器| 欧美日本亚洲视频在线播放| 免费在线观看亚洲国产| 亚洲成人免费电影在线观看| 人妻夜夜爽99麻豆av| 欧美成人性av电影在线观看| 9191精品国产免费久久| 丰满人妻一区二区三区视频av | 99久国产av精品| 国产99白浆流出| 黄片大片在线免费观看| 白带黄色成豆腐渣| 日日干狠狠操夜夜爽| 成人精品一区二区免费| 久久婷婷人人爽人人干人人爱| 亚洲成a人片在线一区二区| 国内精品美女久久久久久| 中国美女看黄片| 婷婷亚洲欧美| 狠狠狠狠99中文字幕| 久久久久免费精品人妻一区二区| 床上黄色一级片| 夜夜躁狠狠躁天天躁| 国产精品,欧美在线| 免费看日本二区| 他把我摸到了高潮在线观看| 国产精品野战在线观看| 亚洲欧美一区二区三区黑人| 全区人妻精品视频| x7x7x7水蜜桃| 热99在线观看视频| 天堂网av新在线| 欧美一区二区亚洲| 黄色日韩在线| 欧美乱妇无乱码| 免费看a级黄色片| 成人18禁在线播放| 国产一区二区亚洲精品在线观看| www日本在线高清视频| www.www免费av| 国产免费一级a男人的天堂| 国产亚洲欧美98| 国产精品久久久久久精品电影| 可以在线观看毛片的网站| 18禁美女被吸乳视频| 夜夜看夜夜爽夜夜摸| 国内精品久久久久久久电影| 国产免费av片在线观看野外av| av专区在线播放| 日日干狠狠操夜夜爽| 国产99白浆流出| 国产69精品久久久久777片| 欧美午夜高清在线| a级一级毛片免费在线观看| 观看美女的网站| 夜夜看夜夜爽夜夜摸| 国产日本99.免费观看| 久久99热这里只有精品18| 美女被艹到高潮喷水动态| 国产高清有码在线观看视频| 色哟哟哟哟哟哟| 国产单亲对白刺激| 欧美黄色片欧美黄色片| 成人永久免费在线观看视频| 黄片大片在线免费观看| 欧美日本亚洲视频在线播放| 久久国产精品人妻蜜桃| 舔av片在线| 麻豆成人av在线观看| 在线免费观看不下载黄p国产 | 嫩草影视91久久| 欧美中文综合在线视频| 九色成人免费人妻av| 国产探花极品一区二区| 一本一本综合久久| av国产免费在线观看| 天堂影院成人在线观看| 美女黄网站色视频| 人妻丰满熟妇av一区二区三区| 一级毛片高清免费大全| 国内少妇人妻偷人精品xxx网站| 日本一二三区视频观看| 午夜两性在线视频| 国产成人影院久久av| 白带黄色成豆腐渣| 亚洲五月天丁香| 欧美区成人在线视频| 免费观看的影片在线观看| 日韩中文字幕欧美一区二区| 日韩欧美精品免费久久 | 老司机在亚洲福利影院| 岛国在线免费视频观看| 网址你懂的国产日韩在线| 日本a在线网址| 99国产精品一区二区蜜桃av| 国产精品99久久久久久久久| 欧美成人性av电影在线观看| 国产精品99久久99久久久不卡| 中亚洲国语对白在线视频| 国产伦一二天堂av在线观看| 午夜亚洲福利在线播放| 手机成人av网站| 丝袜美腿在线中文| 亚洲欧美日韩东京热| 午夜福利在线观看吧| 动漫黄色视频在线观看| 久99久视频精品免费| 色哟哟哟哟哟哟| 精品无人区乱码1区二区| 五月伊人婷婷丁香| 亚洲不卡免费看| 免费电影在线观看免费观看| 黑人欧美特级aaaaaa片| 亚洲av美国av| 怎么达到女性高潮| 90打野战视频偷拍视频| 久久精品国产亚洲av涩爱 | 中文字幕av成人在线电影| 日本三级黄在线观看| 乱人视频在线观看| 国产99白浆流出| 欧美黄色淫秽网站| 国内揄拍国产精品人妻在线| 最近最新中文字幕大全电影3| 午夜福利高清视频| 国产亚洲欧美在线一区二区| 一区二区三区国产精品乱码| 黄色视频,在线免费观看| 欧美高清成人免费视频www| 两个人视频免费观看高清| 国产精品亚洲一级av第二区| 欧美日韩福利视频一区二区| 国产男靠女视频免费网站| 亚洲精品乱码久久久v下载方式 | 内地一区二区视频在线| 丰满的人妻完整版| 欧美最黄视频在线播放免费| 观看美女的网站| 中文亚洲av片在线观看爽| 国产精品美女特级片免费视频播放器| 老汉色∧v一级毛片| 亚洲成a人片在线一区二区| 99久久精品热视频| 可以在线观看毛片的网站| 国产av不卡久久| 欧美+日韩+精品| 免费观看的影片在线观看| 又粗又爽又猛毛片免费看| 亚洲人成网站在线播| 欧美黄色片欧美黄色片| 中文字幕久久专区| 国产伦人伦偷精品视频| 国产视频一区二区在线看| 黄色片一级片一级黄色片| 国产三级在线视频| 日韩有码中文字幕| 一区二区三区高清视频在线| 亚洲第一欧美日韩一区二区三区| 中国美女看黄片| 日本黄色视频三级网站网址| 狂野欧美白嫩少妇大欣赏| 色噜噜av男人的天堂激情|