軌道客車能量分配快速分析方法及應用

李本懷

(中車長春軌道客車股份有限公司 工程實驗室，吉林 長春 130062)*

?

軌道客車能量分配快速分析方法及應用

李本懷

(中車長春軌道客車股份有限公司 工程實驗室，吉林 長春 130062)*

為實現方案設計階段軌道客車吸能系統的快速設計，基于LS-DYNA的碰撞接觸分析技術，提出一種應用剛體和非線性彈簧相結合的客車碰撞性能分析方法.該方法根據列車重量、吸能元件的壓潰力、可用壓潰空間等輸入參數，快速分析軌道客車的碰撞能量.采用該方法研究兩列6輛編組地鐵列車以25 km/h相對速度碰撞的能量分配問題，結果表明該方法對列車吸能系統能量分配高效可行.

軌道客車；吸能系統；碰撞；壓潰行程；能量分配

0 引言

隨著國內城市化進程的推進，各大中型城市因軌道交通大容量、速度快和準時的優勢，相繼將軌道交通作為公共交通優先發展[1].與此同時，車輛運行安全性也備受關注[2].近年來，歐美、日本等發達國家列車被動安全設計技術已經成熟，根據相應的事故研究和分析，制定了各自的碰撞法規[3].列車各界面的緩沖裝置的設計對列車整體碰撞有較大影響，英國OLEO公司和日本東京都開發了自己的一維碰撞仿真分析系統，用于車鉤緩沖、吸能裝置的設計.國內各鐵路院校及主機廠采用三維多車編組模型模擬了列車碰撞過程，對列車碰撞安全性能進行研究[4].但在列車吸能方案的設計階段，三維分析模型規模大，計算時間長，無法滿足主機廠列車吸能系統的設計要求.

基于上述原因，有必要對列車碰撞能量分配方法進行深入研究.基于LS-DYNA的碰撞接觸分析技術[5]，本文提出一種應用剛體和非線性彈簧進行列車碰撞性能分析的方法，可在方案設計階段對列車各斷面的吸能量快速進行預測，為列車車鉤、防爬器及吸能區所需的壓潰力和壓潰行程設計提供依據.通過研究兩列6輛編組地鐵列車相對碰撞的能量分配，驗證了該方法的實用性.

1 能量快速分配方法

能量快速分配方法的基本方程為：

式中:γ為相對體積；ρ為當前質量密度；ρ0初始質量密度.

由于上述偏微分方程的非線性性質和碰撞接觸問題邊界條件的復雜性，無法用解析法來求出精確解.有限元法是適用于求解上述問題的有效數值解法，其特點是首先建立和原問題基本方程及相應定解條件相等效的積分提法，然后據此建立近似解法.如果原問題的方程具有某些特定的性質，則它的等效積分提法可以歸結為某個泛函的變分，相應的近似解法實際上是求解泛函的駐值問題.當問題的微分方程和邊界條件已知，但是變分的泛函尚未找到或者根本不存在的情況，用伽遼金(Galerkin)法(即利用近似解的試探函數序列作為權函數)來確定單元特性和建立求解方程.

基于碰撞接觸分析技術，利用剛體和非線性彈簧進行列車碰撞能量的分配方法，能夠快速預測不同速度的軌道客車碰撞時，車輛端部吸能裝置的吸能狀態，指導吸能結構的優化設計，同時也能驗證軌道客車編組有限元模型碰撞仿真分析結果[6].能量快速分配方法的主要流程：

(1)模型建立：利用HYPERMESH軟件，建立車剛性模型、車鉤緩沖裝置模型(通過對彈簧單元的選取、力-行程曲線編輯和單元組合實現車鉤緩沖器和壓潰管的力學特性)、輪軌模型，主要包括車體重量信息、車體材料自由度狀態、模型單元信息、停放制動信息、運行輪軌信息等；

(2)施加列車自重：定義列車運動或停放制動狀態列車的運動性能，包括運動和靜止列車的重力加速度，通過加速度曲線，施加列車垂向重力加速度，保持列車在運動過程中始終有自重存在;

(3)建立接觸：建立列車碰撞過程中各項接觸關系，運動列車內部接觸、靜止列車內部接觸、運動與靜止列車碰撞過程中接觸、輪軌停放制動接觸和運行接觸；

(4)定義運動列車初始碰撞速度：定義運動列車的初始狀態，根據碰撞標準要求，施加運動列車初始速度；

(5)定義輸出信息：設置列車碰撞終止時間，碰撞輸出步長、打開各項輸出信息開關，通過列車節點集合信息定義速度、加速度、位移信息，通過單元集合定義列車各斷面力-時間曲線、力-行程曲線，通過PART組定義各項能量輸出；

(6)輸出計算文件提交LS-DYNA平臺計算：輸出模型K文件，通過命令提交LS-DYNA進行碰撞仿真計算；

(7)讀取計算結果信息：通過節點輸出文件讀取速度、減速度、位移信息，判定列車碰撞后是否達到相同速度(列車吸能結束)，列車各部位減速度是否滿足要求，兩列列車達到相同速度時行走距離等；通過單元輸出文件讀取列車碰撞過程中力-時間曲線，力-行程曲線，判定列車車鉤及緩沖裝置力和行程是否在許用范圍；通過能量輸出文件讀取列車動能、內能、摩擦耗能及各緩沖裝置吸能情況.依據計算結果進行分析判斷，如果滿足設計要求，分析結束，如果無法滿足設計要求返回(1)步修改設計參數重新計算.

2 工程應用

采用上述快速分配方法研究兩列6輛編組地鐵列車碰撞能量分配問題.列車采用4M2Tc 6輛編組，即：+Tc**M*M*M*M**Tc+，其中M重35 t，Tc重33 t，兩列碰撞編組模型等效形式為+Tc**M*M*M*M**Tc++Tc**M*M*M*M**Tc+，車輛與軌道間摩擦系數停放制動狀態為0.12，運動狀態為0.002，車鉤吸能參數見表1所示[6](“+”表示為全自動車鉤；“**”表示為雙壓潰管半永久牽引桿;“*”表示為單壓潰管半永久牽引桿).

表1 編組列車車鉤吸能參數

軌道客車碰撞吸能模型中，列車沿縱向自由度運動；軌道和車體采用剛體模型，車鉤通過彈簧單元等效模擬，通過列車間各斷面的能量吸收狀態確定壓潰力及行程[6].通過加速度曲線施加列車垂向重力加速度9 800 mm/s2，保持列車在運動過程中始終有自重存在.模型中的接觸包括：運動和靜止列車的內部自接觸、靜止列車與運動列車之間的接觸、列車與軌道間接觸.列車碰撞分析的終止時間為2 s，結果輸出步長為0.01 s.通過列車節點集合信息輸出速度、減速度、位移信息，通過單元集合輸出列車各斷面力-時間曲線、力-行程曲線，通過PART組輸出能量-時間曲線，輸出信息的節點和單元如表2和表3所示.

兩列車以25 km/h相對速度相撞時，車鉤緩沖器、壓潰管及防爬器全部參與吸能.依據表1可知，運動列車和靜止列車第一斷面行程為1 470 mm，即：兩個緩沖器110 mm，兩個壓潰管660 mm，兩個防爬器700 mm；運動列車和靜止列車第二、六斷面行程為510 mm，即兩個緩沖器110 mm，兩個壓潰管400 mm；運動列車和靜止列車第三、四、五斷面行程為310 mm，即兩個緩沖器110 mm，一個壓潰管200 mm.由圖1可知列車碰撞過程中運動列車和靜止列車第一斷面壓潰行程為1 432 mm，可知緩沖器、壓潰管壓潰完全，防爬器壓潰行程為331 mm，未超出最大壓縮行程350 mm；由圖2可知，運動列車和靜止列車第二斷面車鉤壓潰行程分別壓縮463 mm和457 mm，未超出車鉤最大壓縮行程510 mm；由圖3可知，運動列車和靜止列車第三段面車鉤壓潰行程分別壓縮148和146 mm，未超出車鉤最大壓縮行程310 mm，滿足設計要求.在0.78 s時，各列車速度曲線交匯成一條曲線，說明列車各車輛達到共同速度，吸能結束，如圖4所示.車輛吸能裝置完成吸能，剩余能量將由制動摩擦吸收.在0.78 s時，各車輛位移曲線斜率相同，兩輛列車以共同速度運行，列車達到共同速度，如圖5所示.

表2 輸出速度和位移信息的節點

表3 輸出車鉤吸能信息的單元

圖1 第一斷面行程-時間變化曲線圖

圖2 第二斷面車鉤行程-時間變化曲線

圖3 第三、四、五、六斷面車鉤行程-時間變化曲線

圖4 各車輛速度-時間變化曲線

圖5 各車輛位移-時間變化曲線

兩列車以25 km/h相對速度相撞時，前端車鉤在0.3 s時剪斷，防爬器接觸，防爬器壓潰管開始吸能，列車總動能由4.47 MJ變為2.25 MJ，共吸能約2.22 MJ，如圖6～圖8所示.

圖6 列車動能-時間變化曲線

圖7 列車內能-時間變化曲線

圖8 滑移能-時間變化曲線

3 結論

(1)基于碰撞接觸分析技術，提出一種適用于軌道客車吸能系統設計階段的能量分配分析方法.該方法的應用將縮短吸能系統方案的設計周

期，為后續編組列車有限元結構大變形碰撞接觸分析條件提供理論依據；

(2)兩列6輛編組地鐵以25 km/h相對速度相撞時，在0.3 s時前端車鉤剪斷，防爬器接觸，其壓潰管開始吸能，共吸能約2 200 kJ；在0.78 s吸能結束時，車鉤緩沖器、壓潰管及防爬器都沒有超出最大行程，各車輛以共同速度運行.

[1]王文斌，康康，趙洪倫．列車耐碰撞系統有限元和多體動力學聯合仿真[J]．同濟大學學報(自然科學版)，2011，39(10)：1552- 1556.

[2]ALEXANDER SCHOLES(英)．歐洲鐵路碰撞技術的開發[J]．國外鐵道車輛，1998(1)：23- 26.

[3]BRAUN JOHANNES，陳銘.有軌電車和輕軌車輛的防碰撞性[J].國外鐵道車輛，2005,42(5)：39- 41.

[4]李健，沈鋼．列車防碰撞裝置及動力學仿真[J]．鐵道車輛，2001，39(7):5- 8.

[5]趙海鷗.LS-DYNA動力分析指南[M].北京：兵工業出版社，2003.

[6]李本懷，王科飛，李成林．B型不銹鋼地鐵碰撞吸能研究 [C]．第八屆中國智能交通年會論文集，2013.

Fast Analysis Method of Impact Energy Distribution for Railway Vehicle

LI Benhuai

(CRRC Changchun Railway Vehicle Co.，Ltd,Changchun 130062,China)

Based on LS-DYNA collision analysis technology, an approach energy distribution analysis method is proposed for railway vehicle at design stage. It can quickly match energy distribution scheme for vehicles when their weight, crushing force and stroke of energy absorbers are known. A six-units train with velocity of 25 km/h impacting a six-units static train will be calculated for their vehicles energy distribution by the method. The analysis results show that the method is feasible and effective.

railway vehicle; energy absorber; crashworthiness; crush stroke; energy distribution

1673- 9590(2017)01- 0026- 04

2016- 05- 13

高速列車譜系化技術平臺及系列車型研制項目(2012AA112002)

李本懷(1978-)，男，高級工程師，碩士，主要從事車輛結構強度、疲勞、碰撞方面的研究

E-mail:lbh1008@163.com.

A