楊 志,邢宗義,王 露,龍 靜

(1.南京理工大學(xué)自動化學(xué)院,南京 210094； 2.廣州地鐵集團(tuán)有限公司運營總部 ,廣州 510380)

基于磨耗數(shù)據(jù)統(tǒng)計模型的鏇輪決策優(yōu)化

楊 志1,邢宗義1,王 露1,龍 靜2

(1.南京理工大學(xué)自動化學(xué)院,南京 210094； 2.廣州地鐵集團(tuán)有限公司運營總部 ,廣州 510380)

針對目前地鐵車輛車輪等級鏇修而導(dǎo)致的鏇修不合理問題,提出一種基于磨耗數(shù)據(jù)統(tǒng)計模型鏇輪決策優(yōu)化方法?；趶V州地鐵8號線的輪對磨耗數(shù)據(jù),首先利用SPSS分析輪徑磨耗和輪緣磨耗與輪緣厚相關(guān)性；其次建立基于狀態(tài)轉(zhuǎn)移過程的輪緣磨耗模型和基于數(shù)理統(tǒng)計模型的輪徑磨耗模型；然后研究單個車輪的單級和多級鏇修控制限策略,并通過蒙特卡羅仿真分析比較兩種控制限策略的優(yōu)劣。結(jié)果表明,多級鏇修控制限策略大大提高車輪使用壽命,節(jié)約列車運營成本。

輪對磨耗；優(yōu)化；單級鏇修策略；多級鏇修策略

1 概述

隨著軌道交通的不斷發(fā)展，車輪的磨耗一直都是國內(nèi)外學(xué)者關(guān)注的問題。Ansari等[1]建立基于多體動力學(xué)理論的輪對磨損模型，通過分析德黑蘭地鐵車輛輪對磨耗數(shù)據(jù)驗證了模型的可行性。Ward等[2]提出基于輪軌動力學(xué)的車輪磨耗數(shù)據(jù)算法，實現(xiàn)車輪磨耗仿真分析。Arizon等[3]采用解析模型對車輪磨耗進(jìn)行分析，實現(xiàn)車輪磨損故障建模和壽命預(yù)測。黃永達(dá)[4]和刁曉明[5]通過統(tǒng)計武廣高速動車組踏面磨耗數(shù)據(jù)，分析得出踏面磨耗較大時輪緣厚度的磨耗也較大。上述文獻(xiàn)主要從動力學(xué)模型和實際磨耗數(shù)據(jù)等方面進(jìn)行了大量研究，為車輪鏇修優(yōu)化奠定了基礎(chǔ)。

TAO等[6]通過分析D20E型機車車輪的磨耗數(shù)據(jù)，建立基于數(shù)理統(tǒng)計磨耗模型，實現(xiàn)不同類型機車車輪鏇修周期預(yù)測。Braghin等[7]提出了一種輪對磨耗數(shù)值模擬方法，研究表明列車行駛里程達(dá)20萬km時進(jìn)行鏇修，將有效提高車輪使用壽命。員華[8]、朱士友[9]和孫小康[10]均對等級鏇修可行性進(jìn)行了分析，認(rèn)為應(yīng)根據(jù)輪緣厚度選擇不同的鏇修等級。上述方法從車輪磨耗角度進(jìn)行了統(tǒng)計分析并給出鏇修建議，但沒有觸及車輪鏇修決策優(yōu)化的實質(zhì)。

由于列車自身車型的不同和運營線路的差異，僅僅基于動力學(xué)性能分析來制定車輪的鏇修方案無法滿足實際應(yīng)用，而通過統(tǒng)計分析具體線路實際運營中的列車車輪磨耗數(shù)據(jù)來制定或調(diào)整鏇輪維修周期是有效可行的。因此，基于廣州地鐵8號線的磨耗數(shù)據(jù)進(jìn)行分析，提出一種車輪多級鏇修控制限策略模型，以滿足鏇輪維修的實際需求。

2 車輪磨耗與輪緣厚相關(guān)性分析

許宏[11]、王凌[12]和Ansari M[13]研究表明軌道交通的車輪磨耗速率與輪緣厚度有關(guān)，但磨耗率與輪緣厚度的趨勢和相關(guān)性卻不盡相同，車輪磨耗情況與列車型號、運營線路環(huán)境、車輪材質(zhì)等多種因素相關(guān)。因此，建立車輪磨耗模型之前需要對輪徑磨耗和輪緣磨耗相關(guān)性進(jìn)行分析。

將廣州地鐵8號線列車的輪對尺寸數(shù)據(jù)作為樣本數(shù)據(jù)進(jìn)行分析，由式(1)和式(2)計算每間隔30 d輪徑磨耗率和輪緣磨耗率

式中，Di+1和Di分別為不包含鏇輪情況下ti+1和ti時刻輪徑測量值；Sd，i+1和Sd，i分別為不包含鏇輪情況下ti+1和ti時刻輪緣厚度測量值；vd，i和vsd，i分別表示各個時間段輪徑磨耗率和輪緣磨耗率。

為了方便分析車輪磨耗與輪緣厚度的相關(guān)性，將vd，i近似成輪徑為(Di+1-Di)/2的輪徑磨耗率，同理，vsd，i近似成輪緣厚度為(Sd，i+1-Sd，i)/2的輪緣磨耗率。

本文利用SPSS軟件分別對輪徑磨耗率和輪緣磨耗率與輪緣厚度的相關(guān)性進(jìn)行分析，通過對式(1)和式(2)近似得到的輪徑磨耗率和輪緣磨耗率進(jìn)行統(tǒng)計計算，得到輪徑磨耗率與輪緣厚度皮爾森相關(guān)性為-0.05，而輪緣磨耗率與輪緣厚度皮爾森相關(guān)性為-0.286，因此可以認(rèn)為輪緣磨耗率與輪緣厚度是相關(guān)的，而輪徑磨耗率與輪緣厚度是獨立的。

3 輪緣磨耗模型和輪徑磨耗模型

根據(jù)輪緣磨耗率與輪緣厚度相關(guān)，輪徑磨耗率與輪緣厚度相互獨立的分析結(jié)論，分別對車輪輪緣和輪徑的磨耗規(guī)律進(jìn)行建模。

3.1 基于狀態(tài)轉(zhuǎn)移過程的輪緣磨耗模型

根據(jù)地鐵車輛輪對尺寸安全限度要求：輪緣厚度Sd應(yīng)滿足26 mm≤Sd≤32 mm，將輪緣厚安全域分為20個以0.3 mm為長度的狀態(tài)區(qū)間，分別為：(31.7，32]，(31.4，31.7]，(30.1，31.4]，…，(26.3，26.6]，[26，26.3]，各狀態(tài)分別記為S1，S2，…S20。

將上述均勻時間間隔的輪緣厚度樣本數(shù)據(jù)進(jìn)行多狀態(tài)區(qū)間統(tǒng)計，在車輪經(jīng)過30 d磨耗后，輪緣厚度會發(fā)生一定的改變。由于車輪磨耗存在大小差異，因此，經(jīng)磨耗的車輪輪緣厚度值可能會停留在原狀態(tài)區(qū)間，記為S(i，i)，也有可能轉(zhuǎn)移至另外一個狀態(tài)區(qū)間，記為S(i，i+m)，狀態(tài)轉(zhuǎn)移示意如圖1所示。

由于每個狀態(tài)區(qū)間長度均為0.3 mm，若車輪在30 d磨耗后為狀態(tài)停留，對應(yīng)磨耗量為0～0.3 mm，若為狀態(tài)一級轉(zhuǎn)移，則對應(yīng)磨耗量為0.3～0.6 mm，若為狀態(tài)二級轉(zhuǎn)移，則對應(yīng)磨耗量為0.6～0.9 mm，以此類推。將樣本數(shù)據(jù)根據(jù)狀態(tài)轉(zhuǎn)移方式進(jìn)行統(tǒng)計，得到各狀態(tài)轉(zhuǎn)移的頻數(shù)統(tǒng)計情況，進(jìn)而得到每個初始狀態(tài)對應(yīng)的狀態(tài)轉(zhuǎn)移概率。由于狀態(tài)變換情況主要集中在狀態(tài)停留、狀態(tài)一級轉(zhuǎn)移、狀態(tài)二級轉(zhuǎn)移，為方便分析，舍棄三級轉(zhuǎn)移及三級以上的小概率狀態(tài)轉(zhuǎn)移情況。將狀態(tài)停留、狀態(tài)一級轉(zhuǎn)移、狀態(tài)二級轉(zhuǎn)移根據(jù)統(tǒng)計權(quán)重進(jìn)行重新分配，得到調(diào)整后的狀態(tài)轉(zhuǎn)移概率表，如表1所示。

表1 簡化狀態(tài)轉(zhuǎn)移概率

從表1可以得出，狀態(tài)二級轉(zhuǎn)移的概率基本保持在一個穩(wěn)定范圍內(nèi)，說明異常磨耗情況是相對穩(wěn)定的。而狀態(tài)停留、狀態(tài)一級轉(zhuǎn)移會隨輪緣厚度的變化有明顯的變化趨勢，且輪緣磨耗率與輪緣厚度之間相關(guān)性較大，因此二者存在一定的變化規(guī)律。根據(jù)表1作圖2觀察狀態(tài)停留和狀態(tài)一級轉(zhuǎn)移下輪緣磨耗隨輪緣厚度的變化規(guī)律。

圖2 輪緣磨耗與輪緣厚度趨勢變化曲線

圖2為磨耗規(guī)律擬合曲線，能大致反映輪緣磨耗隨輪緣厚度變化過程。狀態(tài)停留概率隨輪緣厚度的減小呈現(xiàn)先增大再減小趨勢，在初始狀態(tài)S6～S9范圍內(nèi)時，其狀態(tài)停留概率較大，說明其磨耗率越低。從總體趨勢上看，狀態(tài)停留輪緣磨耗隨輪緣厚度的減小呈先減小再增大的趨勢。而由于狀態(tài)二級轉(zhuǎn)移概率基本穩(wěn)定，因此狀態(tài)一級轉(zhuǎn)移的變化規(guī)律必然與狀態(tài)停留變化規(guī)律呈相反趨勢。

通過統(tǒng)計各初始狀態(tài)下經(jīng)過磨耗的狀態(tài)轉(zhuǎn)移概率，可以得到車輪輪緣磨耗過程中每個初始狀態(tài)分別對應(yīng)3種轉(zhuǎn)移的可能性，每種轉(zhuǎn)移的可能性由統(tǒng)計的概率決定，滿足均勻分布。

3.2 基于數(shù)理統(tǒng)計理論的輪徑磨耗模型

根據(jù)廣州地鐵8號線近5年輪徑磨耗率數(shù)據(jù)進(jìn)行統(tǒng)計分析，得出輪徑磨耗率基本在-3～3 mm/30 d，為方便統(tǒng)計分析，將區(qū)間[-3,3]平均分為60個小區(qū)間，每個區(qū)間長度為0.1 mm。對每個區(qū)間內(nèi)的輪徑磨耗率的頻數(shù)進(jìn)行統(tǒng)計，得出如圖3所示輪徑磨耗率統(tǒng)計直方圖。其中，輪徑磨耗率在-0.2～-0.1 mm/30 d的頻數(shù)最大，主要集中在-0.3～0.1 mm/30 d，輪徑磨耗率向兩側(cè)延伸，磨耗率頻數(shù)逐漸減小。

為得到輪徑磨耗率統(tǒng)計規(guī)律，分別采用正態(tài)分布、對數(shù)正態(tài)分布、威爾布分布和Γ分布對輪徑磨耗進(jìn)行建模。將輪徑磨耗率頻數(shù)統(tǒng)計結(jié)果導(dǎo)入Matlab中，可以擬合計算得到4種分布概率密度函數(shù)具體參數(shù)。

圖3 輪徑磨耗率頻數(shù)統(tǒng)計直方圖

正態(tài)分布概率密度函數(shù)

式中，μ=-0.173 3，σ=0.543 8。

對數(shù)正態(tài)分布概率密度函數(shù)

式中，μ=1.0188，σ=0.2134。

威布爾分布概率密度函數(shù)

式中，α=4.8187，β=0.0047。

Γ分布概率密度函數(shù)

式中，α=24.8086，β=0.1139。

表2 擬合優(yōu)度檢驗結(jié)果

根據(jù)表2可以得到，通過χ2擬合優(yōu)度檢驗，輪徑磨耗分布服從正態(tài)分布和Γ分布，而對數(shù)正態(tài)分布和威布爾分布有顯著差異。因此舍棄對數(shù)正態(tài)分布和威布爾分布模型，對正態(tài)分布和Γ分布進(jìn)行比較分析。

根據(jù)圖3知輪徑磨耗率頻數(shù)在-0.2～-0.1 mm/月所占比例最大，而由各分布的概率密度函數(shù)可得正態(tài)分布和Γ分布分別在輪徑磨耗率為-0.173 3 mm/月和-0.301 6 mm/月時所占比例最大，可知正態(tài)分布更符合輪徑磨耗分布。且通過計算實際輪徑磨耗率在-3～1 mm/30 d概率與正態(tài)分布和Γ分布進(jìn)行比較，得出實際輪徑所占總體概率為0.974 6，而正態(tài)分布和Γ分布分別為0.984 5和0.954 7。綜上計算結(jié)果，確定正態(tài)分布為輪徑的磨耗模型。

4 單輪等級鏇修控制限策略模型

通過輪緣磨耗模型和輪徑磨耗模型模擬仿真車輪磨耗過程，建立單個車輪的單級和多級鏇修控制限策略模型，并通過蒙特卡羅仿真對比分析兩種策略方案優(yōu)劣得出最優(yōu)的控制限策略方案。

4.1 單級鏇修控制限策略模型

單個車輪的鏇修策略模型主要考慮車輪的輪緣厚度的控制限范圍，即車輪何時鏇修，鏇修恢復(fù)至多少，兩者分別對應(yīng)輪緣厚度控制限的下限和上限。在建立輪緣磨耗模型時提出了狀態(tài)轉(zhuǎn)移的概念，并對車輪輪緣厚度的安全閾進(jìn)行了狀態(tài)劃分，現(xiàn)做出如下描述：當(dāng)車輪的輪緣厚度磨耗至某一狀態(tài)Sj，再將其恢復(fù)至某一狀態(tài)Si。其中，Si對應(yīng)控制限上限狀態(tài)，記為Sh；Sj對應(yīng)控制限下限狀態(tài)，記為Sl，且控制限上下限狀態(tài)分別對應(yīng)各自的區(qū)間范圍，其區(qū)間分別表示為[Shr，Shl]和[Slr，Sll]。

單級控制限僅對應(yīng)一個上限和一個下限，即車輪輪緣厚度磨耗至一個統(tǒng)一的下限，恢復(fù)至同一個上限。當(dāng)車輪的輪緣厚度經(jīng)過一個磨耗周期(30 d)磨耗至小于等于控制限下限狀態(tài)區(qū)間左值Sll時，對車輪進(jìn)行鏇修，將輪緣厚度恢復(fù)至控制限上限狀態(tài)區(qū)間Shl處，單個車輪單級鏇修控制限策略示意如圖4所示。輪緣在磨耗過程中，遵循輪緣磨耗狀態(tài)轉(zhuǎn)移概率表。

圖4 單級鏇修控制限策略示意

根據(jù)車輪單級鏇修控制限原理計算得到單級鏇修有210種策略方案。本文鏇修控制限策略的制定原則是保證輪緣厚度維持在26～32 mm安全閾范圍內(nèi)，同時得到較優(yōu)的控制限策略，以延長車輪的使用壽命并保證較小的鏇修次數(shù)。因而將車輪使用壽命和鏇修次數(shù)作為優(yōu)化目標(biāo)，經(jīng)過蒙特卡羅仿真得到每種鏇修策略方案的期望使用壽命和期望鏇修次數(shù)，從而選出較優(yōu)的策略方案，仿真流程如圖5所示。

圖5 單級鏇修控制限蒙特卡羅仿真流程

單級鏇修蒙特卡羅仿真步驟如下。

①參數(shù)初始化：設(shè)置未經(jīng)磨耗的車輪輪徑和輪緣厚度初始值分別為D(0)=840 mm，Sd(0)=32 mm，壽命累計周期T=0，鏇修次數(shù)N=0，兩次鏇修間隔時間t=0；控制限上下限狀態(tài)區(qū)間為Sh和Sl，對應(yīng)鏇修恢復(fù)上下限值為Shl和Sll，此處將分別導(dǎo)入210種策略對應(yīng)的值。

②本次磨耗時間周期累計t=t+1。

③根據(jù)T+t-1時刻的車輪輪徑值D(T+t-1)和輪緣厚度值Sd(T+t-1)，由輪徑磨耗正態(tài)分布模型隨機生成輪徑磨耗概率pD以及狀態(tài)轉(zhuǎn)移概率均勻分布隨機生成狀態(tài)轉(zhuǎn)移概率pSd，并根據(jù)磨耗概率得到輪緣磨耗率vd和輪徑磨耗率vSd，通過D(T+t)=D(T+t-1)-vd和Sd(T+t)=Sd(T+t-1)-vSd計算當(dāng)前時刻輪徑值和輪緣厚度值。

④判斷當(dāng)前時刻輪徑值D(T+t)是否≥770 mm，若是，則轉(zhuǎn)入第⑤步，否則結(jié)束此次仿真，得到最后的壽命累積周期Tend和鏇修次數(shù)Nend。

⑤判斷當(dāng)前時刻輪緣厚度值Sd(T+t)是否小于Sll，若是，轉(zhuǎn)入第⑥步，否則轉(zhuǎn)入第②步。

⑥對車輪進(jìn)行鏇修，恢復(fù)輪緣厚度至Shl，并對參數(shù)進(jìn)行初始化：T=T+t，N=N+1，t=0，D(T)=D(T)-k×(Shl-Sd(T))，Sd(T)=Shl，初始化完成后轉(zhuǎn)入第②步。其中，k為鏇修比例系數(shù)，根據(jù)龐松林[14]和趙文杰[15]研究成果，取鏇修比例系數(shù)為4.2。

4.2 多級鏇修控制限策略模型

單級控制限策略在鏇修過程中以損失較大車輪直徑為代價，縮短車輪使用壽命，且單級控制限策略僅針對單個車輪，若因為單個車輪磨耗超過控制限下限就進(jìn)行鏇修工作，則會造成列車頻繁下線，給實際鏇修工作和運營計劃帶來困難。因此，為了延長列車使用壽命以及滿足實際鏇修需求，提出多級鏇修控制限策略模型。

對應(yīng)單級鏇修控制限，多級鏇修控制限有多個鏇修上限和多個鏇修下限，并且上限和下限成對出現(xiàn)，多級控制限策略示意如圖6所示。

圖6 多級鏇輪控制限示意

在多級鏇修控制限模型中，控制策略種類繁多，兩級控制鏇修控制限策略有1 140種，三級控制鏇修控制限策略有3 876種。為簡化工作以及檢驗方法的有效性，以兩級鏇修控制限策略為例進(jìn)行蒙特卡羅仿真。與單級鏇修控制限策略仿真過程不同，兩級鏇修控制限策略初始化時需要對兩級鏇修上下限進(jìn)行設(shè)置。在判斷鏇修條件時會進(jìn)行兩級輪緣厚度閾值范圍判斷，即當(dāng)前輪緣厚度如果滿足一級鏇修條件則進(jìn)行一級鏇修，如果當(dāng)前輪緣厚度值不滿足一級鏇修條件，則轉(zhuǎn)入二級鏇修條件進(jìn)行判斷，并進(jìn)行相應(yīng)的操作。通過上述過程實現(xiàn)單個車輪的多級鏇修控制限策略仿真。

5 仿真與結(jié)果分析

對單級和兩級鏇修控制限進(jìn)行蒙特卡羅仿真分析，每種策略方案設(shè)置1 000次的仿真次數(shù)，取均值得到每種策略方案的期望使用壽命E(T)和期望鏇修次數(shù)E(N)。

單級鏇修控制限策略仿真如圖7所示。其中圖7(a)為210種策略方案編號示意，且紅點代表控制限上限左值，藍(lán)點代表控制限下限左值。圖7(b)和圖7(c)分別為每種策略方案所對應(yīng)的車輪期望使用壽命周期和車輪期望鏇修次數(shù)。

圖7 單級鏇修控制限策略仿真結(jié)果

根據(jù)單級鏇修控制限策略仿真圖可以得出以下結(jié)論。

①車輪最大使用壽命策略方案編號為137，最小鏇修次數(shù)策略方案編號為20，可知使用壽命最大方案鏇修次數(shù)卻不是最小，同理鏇修次數(shù)最小方案使用壽命卻不是最大。因此在選擇鏇修方案應(yīng)綜合考慮使用壽命和鏇修次數(shù)。

②在下限值相同的鏇修方案，隨著上限值增加，車輪使用壽命呈先增大后減小過程，與狀態(tài)轉(zhuǎn)移概率分布中狀態(tài)停留概率隨輪緣厚度先增大后減小的現(xiàn)象相符。

③在下限值相同的鏇修方案中，隨著上限值增加，車輪的期望鏇修次數(shù)均隨鏇修上限值的增大而減小，即鏇修策略上下限之間差值越小，期望鏇修次數(shù)反而越大。

綜合考慮車輪的期望使用壽命和鏇修次數(shù)，針對單個車輪的單級鏇修控制限策略選出4種較優(yōu)方案，如表3所示。

表3 單級鏇修控制限策略較優(yōu)方案

單個車輪兩級鏇修控制限策略多達(dá)1 140種，但由單級鏇修控制限策略仿真結(jié)果可知，車輪的期望壽命與狀態(tài)轉(zhuǎn)移概率分布中的狀態(tài)停留概率相關(guān)，控制限中狀態(tài)停留概率越大，則對應(yīng)的車輪期望使用壽命越大。根據(jù)圖2可知，狀態(tài)停留概率較大的狀態(tài)主要為S6～S9，則當(dāng)一級鏇修控制限范圍與二級鏇修控制限范圍過大或過小都不利于車輪使用壽命，因此選擇一級控制限下限值和二級控制限上限左值都比較集中的情況進(jìn)行分析。

選取控制限策略方案編號為749～832進(jìn)行分析，此時一級控制限下限左值為29.9，二級控制限上限左值為29.6，其兩級鏇修控制限策略仿真結(jié)果如圖8所示。其中圖8(a)為從749～832對應(yīng)的84種策略方案的編號示意，圖8(b)和圖8(c)分別為每種方案所對應(yīng)的期望使用壽命和期望鏇修次數(shù)。

圖8 兩級鏇修控制限策略仿真結(jié)果

根據(jù)兩級鏇修控制限策略仿真圖可以得出以下結(jié)論：

①與單級鏇修控制限策略仿真類似，使用壽命最大方案鏇修次數(shù)卻不是最小，鏇修次數(shù)最小方案使用壽命卻不是最大；

②在二級控制限上下限值相同時，隨著一級控制限上限值增大，車輪期望使用壽命先增大后減小，鏇修次數(shù)隨上限值的增大而減小；

③隨著二級控制限下限值的增大，期望使用壽命整體趨勢先增大后減小，而鏇修次數(shù)整體趨勢增大。

綜合考慮車輪的期望壽命和期望鏇修次數(shù)，針對兩級控制限策略優(yōu)選出4種較優(yōu)的方案，如表4所示。

對比單級與多級鏇修控制限策略，多級控制限策略總體優(yōu)于單級控制限策略，主要表現(xiàn)在以下兩個方面：

(1)從運營成本的角度，多級控制限策略的車輪期望使用壽命較單級控制限策略長，因此能更有效地節(jié)約列車的運營成本；

(2)從實際應(yīng)用角度，多級控制策略的期望鏇修次數(shù)雖然提高不大，并且隨著二級控制限下限值增大的同時，期望鏇修次數(shù)總體也呈增大的趨勢，但卻增加了鏇修時間的靈活性，更易于應(yīng)用于實際作業(yè)。

表4 兩級鏇修控制限策略較優(yōu)方案

6 結(jié)論

本文研究了一種基于磨耗數(shù)據(jù)統(tǒng)計模型的鏇輪決策優(yōu)化方法，建立了基于狀態(tài)轉(zhuǎn)移過程的輪緣磨耗模型和基于數(shù)理統(tǒng)計理論的輪徑磨耗模型，并在此基礎(chǔ)上建立單個車輪單級和多級控制限鏇修策略模型。通過蒙特卡羅仿真得到多級鏇修策略模型將有效提高車輪使用壽命，節(jié)約列車運營成本。在研究過程中，僅對兩級鏇修策略進(jìn)行研究，未對三級及三級以上鏇修進(jìn)行仿真分析，未檢驗是否有更優(yōu)的多級控制限策略。

[1] M Ansari， I A Hazrati， E Esmailzadeh， et al. Wear rate estimation of train wheels using dynamic simulations and field measurements[J]. Vehicle System Dynamics， 2008， 46(8):739-759.

[2] A Ward， R Lewis， R S Dwyer. Incorporating a railway wheel wear model into multi-body simulations of wheelset dynamics[J]. Tribology， 2003，41(3):367-376.

[3] J De Arizon， O Verlinden， P Dehombreux. Prediction of wheel wear in urban railway transport: comparison of existing models[J]. Vehicle System Dynamics， 2007，45(9):849-866.

[4] TAO VAN CHIEN，李芾，戚壯，等.機車車輪磨耗統(tǒng)計數(shù)據(jù)處理方法與鏇修周期預(yù)測模型[J].鐵道學(xué)報，2015(12):14-19.

[5] 黃永達(dá).武廣客專輪對磨耗規(guī)律分析[J].甘肅科技，2012，28(18):48-49.

[6] 刁曉明，朱韶光，董孝卿.武廣客專動車組車輪磨耗及振動性能跟蹤研究[J].鐵道機車車輛，2013(2):1-6.

[7] F Braghin， R Lewis， R S Dwyer Joyce， et al. A mathematical model to predict railway wheel profile evolution due to wear[J]. Wear， 2006， 261(11-12):1253-1264.

[8] 員華，肖勝強，汪洋.基于磨耗量統(tǒng)計的輪對等級鏇修可行性分析[J].城市軌道交通研究，2006，9(1):43-45.

[9] 朱士友，潘麗莎，員華.輪對等級鏇修對車輛平穩(wěn)性的影響分析[J].城市軌道交通研究，2006，9(7):46-48.

[10] 孫小康.地鐵輪對踏面鏇修經(jīng)濟(jì)性分析[J].城市軌道交通研究，2016，19(5):91-92.

[11] 許宏，員華，王凌，等.基于高斯過程的地鐵車輛輪對磨耗建模及其鏇修策略優(yōu)化[J].機械工程學(xué)報，2010，46(24):88-95.

[12] 王凌，員華，那文波，等.基于磨耗數(shù)據(jù)驅(qū)動模型的輪對鏇修策略優(yōu)化和剩余壽命預(yù)報[J].系統(tǒng)工程理論與實踐，2011，31(6):1143-1152.

[13] M Ansari， I H Ashtiyani. Determination of the train wheel wear trend， comparing with field measurements[C]∥ Rail Conference， 2006. Proceedings of the 2006 IEEE/ASME Joint. IEEE， 2006:103-107.

[14] 龐松林.CRH2型高速動車組車輪鏇修策略的研究[D].北京：北京交通大學(xué)，2016.

[15] 趙文杰.地鐵車輛輪對磨耗故障預(yù)報和鏇修策略優(yōu)化[D].杭州：中國計量學(xué)院，2014.

Optimization of Wheel Re-profiling Strategy Based on Statistical Wear Model

YANG Zhi1, XING Zong-yi1, WANG Lu1, LONG Jing2

(1. College of Automation, Nanjing University of Science and Technology, Nanjing 210094, China; 2. Guangzhou City Underground Railway Corporation, Guangzhou 510380, China)

Aiming at improving the current unreasonable wheel-set re-profiling schedule due to class-lathing, a re-profiling strategy is proposed based on historical wear data. Firstly, the SPSS is used to analyze the interdependency between wheel diameter wear and wheel flange width based on the wheel wear data of Guangzhou Metro Line 8. Secondly, flange width wear model based on state transition and wheel diameter wear model based on mathematical statistics are built separately, and single-stage and multi-stage planned turns strategies of individual wheels are studied based on the two models. Finally, Monte Carlo simulation is conducted to compare the two strategies, and the result shows that multi-stage planned turns strategy can prolong expected life of wheels and save operating cost effectively.

Wheel-set wear; Optimization; Single-stage planned turns strategy; Multi-stage planned turns strategy

1004-2954(2018)01-0142-06

2017-03-17；

2017-06-22

國家重點研發(fā)計劃項目(2016YFB1200402)，廣州市產(chǎn)學(xué)研協(xié)同創(chuàng)新重大專項(201508010010)，國家工程實驗室建設(shè)項目(2016582)

楊 志(1992—)，男，江蘇連云港人，碩士研究生，從事軌道交通車輛檢測分析研究，E-mail:1145026112@qq.com。

U455.43

A

10.13238/j.issn.1004-2954.201703170006