李 飛,吳興文,羅 贇,劉開成,謝晨希
(西南交通大學 牽引動力國家重點實驗室,成都 610031)
地鐵作為城市軌道交通的重要交通工具,具有運量大、準時、安全等一系列優點[1],近年來隨著地鐵運營公司行車密度及地鐵車輛運量的提高,加速了地鐵線路的惡化,更容易出現鋼軌波磨、車輪多邊形磨耗等輪軌短波不平順。這些輪軌短波不平順極易導致輪軌高頻激勵,從而導致結構高頻共振疲勞問題。
國內學者針對輪軌短波不平順引發結構振動疲勞問題做了大量理論和試驗研究工作。連青林等[2]針對客車轉向架安全吊座孔附近產生的疲勞裂紋問題進行了有限元分析和線路試驗,結果表明輪軌沖擊頻率與安全吊桿振動頻率相近,從而導致共振疲勞失效問題。石懷龍等[3]針對地鐵車輛軸箱吊耳斷裂問題采用線路試驗開展斷裂機理研究;通過分析振動激擾源和結構響應特性,發現波長為61.5 mm 鋼軌波磨是導致車輛振動水平激增的主要原因,同時鋼軌波磨通過頻率在吊耳固有模態頻帶內,進一步導致結構共振從而引發疲勞破壞;作者同時提出鋼軌打磨、優化吊耳結構設計和使用管理條件等解決措施。何偉等[4]針對某地鐵線路車輛在運行中出現一系鋼彈簧疲勞斷裂的現象,建立了一系鋼彈簧柔性的剛柔耦合動力學模型,結果表明地鐵車輛車輪多邊形的階次、波深以及列車的運行速度對一系鋼彈簧的疲勞壽命都有很大的影響,且當車輪多邊形磨耗通過頻率與一系鋼彈簧固有頻率接近時,其壽命顯著降低。覃才等[5]針對某型高速動車組吊耳裝置異常振動且個別出現明顯疲勞裂紋問題,采用有限元法建立吊耳裝置有限元分析模型,分析吊耳裝置的振動特性、應力特性和疲勞壽命,研究吊耳裝置的厚度和倒圓角半徑對其振動特性、應力水平和疲勞壽命的影響,為吊耳裝置的設計提供參數范圍。
文中針對某地鐵車輛轉向架構架端部開裂問題,采用線路試驗的方法,從時域和頻域的角度分析振動激擾源和軸箱、構架端部的響應特性,確定了構架端部開裂原因;同時進一步通過試驗手段驗證了改進方案的效果。
研究對象為某B 型地鐵轉向架構架端部,其在運營過程中出現了多起疲勞失效問題,如圖1 所示。該轉向架構架端部為懸臂梁結構,主要用于安裝輪噴裝置和天線梁等附屬部件;對于沒有安裝附屬部件的構架端部,則為明顯的懸臂梁結構。通過對失效部位的分析,所有端部的失效位置均為懸臂梁與構架下蓋板的T 型焊接接頭處。在實際運營過程中,該位置可能承受從一系懸掛(鋼彈簧和一系液壓垂向減振器)傳遞的輪軌高頻激勵,從而激發構架端部懸臂結構的模態振動。

圖1 構架端部開裂位置
為了弄清構架端部失效原因,應該首先弄清轉向架構架整體模態和構架端部局部模態之間的關系,轉向架一系懸掛系統高頻振動傳遞特性,以及構架端部振動加速度與端部關鍵位置動應力之間的關系;同時還需要揭示轉向架構架端部振動與輪軌表面狀態之間的關聯關系。為此,通過線路實測方法,測試了轉向架構架振動加速度和關鍵部位動應力數據,以及典型線路區段的鋼軌表面狀態。試驗數據表明,該轉向架車輪表面狀態較好,沒有明顯的車輪多邊形磨耗,因此需要重點關注鋼軌表面狀態—輪軌高頻激勵—構架端部失效的相互關聯。
轉向架構架工作模態為構架在外部激勵條件下表現出來的響應模態,其反映了構架在實際服役過程中主要表現的模態或者振動模式,對結構的失效分析具有重要意義。為此文中基于Poly-MAX 方法,利用加速度傳感器對轉向架構架工作模態進行了識別。
PolyMAX 方法是多參考點最小二乘復頻域法的擴展,利用輸出響應的功率譜密度近似代替頻響函數,使用極大似然估計使誤差最小化,實現全局模態的識別,有很好的抗干擾能力[6]。PolyMAX方法主要分為以下3 步[7-9]:
(1)建立線性化的直交矩陣模型。頻域系統中頻響函數H(ω)與輸入參數矩陣F(ω)和輸出參數矩陣X(ω) 之間的關系可以表示為式(1)~式(4)[7-8]:
式中:m為輸入通道總數;l為輸出通道總數;Δt為采樣間隔;ω為角頻率;Z為多項式基函數;P為數學模型的階次;αr,βr分別為分母和分子矩陣的多項式系數,可由最小二乘法求得。
(2)求解極點和模態參與因子。在上述求得分母矩陣系數αr的基礎上對其擴展的酉矩陣進行特征值分解,可得系統的極點及模態參與因子為式(5):
式中:V為特征向量矩陣;Λ為酉矩陣的特征值矩陣(對角矩陣),其對角元素λi(i=1,2,…,mP)與系統極點Pi和間存在的關系為式(6):
而模態阻尼比為式(7):
式中:ξi為阻尼比;σi為阻尼因子;ωi為無阻尼系統得固有頻率;ωdi為有阻尼系統固有頻率。
取特征向量矩陣[V]mP×mP最下面的m行,構成的子矩陣就是模態參與因子矩陣[L]mP×mP。
(3)獲取模態振型。采用最小二乘頻域法(FSFD)求解模態振型,其擬合函數為式(8):
式中:[H(ω) ]l×m為測量的頻響函數矩陣;ψi為待求的第i階模態振型列向量;<l>為模態參預因子行向量;LR和UR分別為分析頻帶外模態影響的下殘余項和上殘余項矩陣。
前文已求出極點Pi和模態參預因子向量<l>,故H(ω)可按不同采樣頻率列方程,最后通過線性最小二乘估計求得未知的模態振型ψi(i=1,2…,N),以及上下殘余項矩陣LR和UR。
為了獲得更加準確的構架端部運營模態,在構架主體粘貼29 個三向振動加速度傳感器,如圖2所示。截取車輛在典型區間內構架振動測點的加速度數據,獲取其頻響函數,然后基于PolyMAX 方法求解構架模態的穩態圖如圖3 所示。圖中o 表示未找到極點;v 表示頻率和模態參預因子穩定;f表示在給定精度內頻率穩定;d 表示頻率和阻尼穩定;s 表示3 種參數全部穩定[10]。

圖2 構架振動測點簡圖

圖3 穩態圖
綜合多個工況下的識別結果,識別出構架主要模態參數和振型見表1,其中構架端部存在頻率為222.9 Hz,阻尼比為0.72% 的縱向擺動固有模態。固有模態如圖4 所示。

表1 轉向架構架主要模態參數及振型

圖4 構架端部擺動模態
對該型地鐵車輛進行線路跟蹤測試,車輛在AW3 即重車狀態下在運營線上進行試驗,列車采用自動駕駛模式ATO 速度運行,常用速度為60 km/h,最大速度不超過120 km/h。由于軸箱振動加速度大小能夠很好地反映線路激擾水平,構架端部振動加速度和裂紋區域應力狀態能夠直接表明激擾對開裂區域的影響,故在軸箱和構架端部安裝加速度傳感器以獲取加速度時域數據;并在裂紋區域粘貼應變片,獲取關注部位動應力數據。加速度傳感器、應變片測點位置如圖5 所示,這些加速度和動應力信號最終通過信號線匯集在終端采集系統。

圖5 構架振動和動應力測點簡圖
選取軸箱振動加速度較大的某個區間作為典型區間進行分析。該區間軌道類型主要包括普通短軌枕整體道床和彈性短軌枕兩種。典型區間內軸箱垂向和構架端部縱向振動加速度時頻圖如圖6、圖7 所示,由圖可知18 mm 彈性短軌枕區段60~80 s 內振動加速度明顯大于其他區段,其中18 mm彈性短軌枕區段振動加速度較大區段頻帶范圍為190~230 Hz,普通短軌枕區段振動加速度較大區段頻帶范圍為80~100 Hz,其余區段振動能量較弱,對構架疲勞破壞影響較小。對彈性短軌枕區段60~75 s 內軸箱振動垂向和構架端部縱向加速度時域數據進行頻譜分析得其主頻分別為212.3 Hz和221.1 Hz。

圖6 軸箱垂向振動加速度特征

圖7 構架端部縱向振動加速度特征
構架端部裂紋區域兩個應力測點處的動應力時程曲線和高應力區段的頻譜如圖8、圖9 所示,其中18 mm 彈性短軌枕區段60~75 s 內動應力幅值明顯大于其他區段,應力范圍超過180 MPa,且頻次較多,表明該區域處于高應力水平,其他區段應力幅值較小,對高應力時段信號進行頻譜分析發現裂紋區域測點信號主頻均為221 Hz 左右。

圖8 構架端部D21 應力測點時頻圖

圖9 構架端部D22 應力測點時頻圖
對典型區間進行軌道調查,發現彈性短軌枕區段存在波長為80 mm 左右、波深為0.2 mm 左右的鋼軌波磨如圖10 所示,且波磨區段占彈性短軌枕區段總長的33%,其余彈性短軌枕區段無明顯鋼軌波磨。普通短軌枕區段存在波長為200 mm 左右的鋼軌波磨,但幅值較小,占比較低。已知該區段車速常用速度為60 km/h,由公式f=v/l可得波長為80 mm 的鋼軌波磨通過頻率為208.3 Hz,由此可見波磨通過頻率與軸箱垂向振動加速度主頻和構架端部縱向加速度主頻接近,說明來自線路波長為80 mm 的鋼軌波磨激擾通過軸箱傳遞到構架端部,并與構架端部固有頻率222.9 Hz 共振,引發構架端部振動疲勞問題。

圖10 軌道不平順曲線圖
由上述數據分析可知,構架端部固有頻率與鋼軌波磨通過頻率共振引發構架端部疲勞破環。為了抑制構架端部縱向振動,在構架端部安裝加強筋并在同一線路相同工況下對改進后車輛進行線路測試,如圖11 所示。統計全線15 個區間構架端部縱向振動加速度水平,對每個站間加速度時間歷程進行0.2 s 滑移峰值計算,計算每個窗口內99.85%可靠度下的峰值,再取最大值,得到站間的最大值,構架端部改進前后加速度最值對比如圖12 所示。安裝加強筋前后典型區間內構架端部測點動應力和加速度時域數據對比如圖13~圖15 所錄,發現安裝加強筋后可以顯著降低端部縱向振動加速度和動應力。對改進方案60~75 s 內的實測數據進行頻譜分析,如圖16~圖18 所示,發現安裝加強筋后構架端部振動加速度和動應力主頻由原方案的221 Hz 左右下降到190 Hz 左右,能量降低90%以上,說明改進方案效果明顯。

圖11 構架端部安裝加強筋

圖12 構架端部縱向振動加速度區間最值對比

圖13 構架端部應力D21 測點對比

圖14 構架端部應力D22 測點對比

圖15 構架端部縱向振動加速度對比

圖16 改進方案D21 應力測點頻譜圖對比

圖17 改進方案D22 應力測點頻譜圖對比

圖18 改進方案構架端部縱向加速度頻譜圖對比
損傷是用來表征載荷作用對材料造成的傷害,通常認為結構在承受疲勞載荷的情況下損傷是逐漸累積的,文中基于實測數據對構架端部關鍵測點的動應力數據進行疲勞損傷分析。
首先對實測得到的動應力數據進行數據處理,去掉異常點和零飄,再基于雨流計數法對實測的動應力數據進行統計,得到對應測點的循環次數和應力范圍,然后以累積頻次為橫坐標,應力范圍為縱坐標,最終得到測點應力譜如圖19、圖20 所示,改進方案與原始方案總累計頻次值接近,但改進方案應力范圍大幅減小。最后基于Miner 線性損傷累計理論,參照結構焊接接頭的S-N曲線進行損傷計算。

圖19 D21 測點改進前后應力譜對比

圖20 D22 測點改進前后應力譜對比
由于構架端部的失效位置均為T 型焊接接頭,根據IIW 國際焊接學會提供的S-N曲線,選取T 型接頭對應的S-N曲線,2×106次循環對應的應力循環范圍為80 MPa;循環次數小于107次時,S-N曲線斜率m=3;同時考慮小載荷對高周疲勞的影響,對于循環次數大于107次時,S-N曲線斜率m=5,S-N曲線如圖21 所示。

圖21 S-N 曲線
Miner 累計損傷理論:在任何應力循環作用下結構都會產生疲勞損傷,疲勞損傷程度與應力幅作用的頻次有關,與載荷譜加載歷程、加載方式無關,應力加載所造成的損傷增量可線性疊加。對應的損傷可以表示為式(9):
結合IIW 標準給出的材料S-N曲線和Miner 線性累計損傷法則可得式(10):
式中:ni為載荷譜中應力范圍Δσi的循環次數;Ni為在該應力范圍時將導致損壞的循環總數;m、C為材料S-N曲線參數。當總的損傷D=1 時,結構發生破壞。
將Miner 線性累計損傷理論公式轉化為里程公式,可得式(11):
式中:L為安全運行公里數;L1為實測公里數。
構架端部動應力測點損傷對比如圖22 所示,結果表明:安裝加強筋后構架端部測點損傷值大幅降低,滿足轉向架的使用壽命要求。

圖22 構架端部動應力測點損傷對比
針對某B 型地鐵車輛轉向架構架端部開裂問題開展運營模態分析和線路試驗研究,分析了構架端部開裂原因及改進方案效果。結論表明:
(1)基于PolyMAX 方法識別出構架在實際運營的工作模態,其中構架端部存在頻率為222.9 Hz、阻尼比為0.72%的縱向局部擺動模態。其與彈性短軌枕區域波長80 mm 通過頻率為200 Hz 左右的鋼軌波磨激勵頻率會發生耦合共振,從而導致局部路段構架端部振動加劇,關鍵位置動應力增加,進一步導致構架端部振動疲勞問題。
(2)構架端部安裝加強筋后使構架端部局部剛度和阻尼比增加,在相同工況條件下改進后結構振動加速度和動應力均比原方案減小90% 左右,使得改進方案損傷顯著減小,滿足轉向架的使用壽命要求。