列車撞擊荷載的有限元數(shù)值分析

晏啟祥， 李 彬， 張 蒙， 何 川， 楊文波， 耿 萍

（1．西南交通大學(xué)交通隧道工程教育部重點實驗室，四川成都610031；2．中鐵工程設(shè)計咨詢集團(tuán)有限公司，北京100020）

列車撞擊荷載的有限元數(shù)值分析

晏啟祥1， 李 彬1， 張 蒙2， 何 川1， 楊文波1， 耿 萍1

（1．西南交通大學(xué)交通隧道工程教育部重點實驗室，四川成都610031；2．中鐵工程設(shè)計咨詢集團(tuán)有限公司，北京100020）

為獲得列車脫軌撞擊荷載，分析列車撞擊時盾構(gòu)隧道的動力響應(yīng)，建立了列車編組的三維撞擊有限元模型，探討了不同列車編組、不同撞擊速度和撞擊角度下列車近似撞擊力時程曲線，分析了列車撞擊力最大值和撞擊時間與列車撞擊速度和角度的關(guān)系，并將典型撞擊荷載用于分析不同厚度二次襯砌管片襯砌的動力響應(yīng)．結(jié)果表明：列車編組數(shù)量一定時，列車斜向撞擊力最大值隨撞擊速度和撞擊角度增大而增大；當(dāng)撞擊角度增大到7．5°后，撞擊力作用時間隨撞擊速度增大而延長；根據(jù)列車撞擊力最大值出現(xiàn)時刻不同，可將撞擊力時程曲線劃分為2類特征曲線，其中第1類特征曲線（撞擊瞬間撞擊力達(dá)到最大）總體上符合高斯多峰擬合公式，可用10個參數(shù)近似擬合．二次襯砌厚度增大能有效減小管片襯砌應(yīng)力、速度、加速度等動力響應(yīng)以及拉、壓損傷區(qū)域．

列車撞擊荷載；特征曲線；雙層襯砌；動力響應(yīng)；拉壓損傷

Key words：train impact load；characteristic curve；double lining；dynamic response；tension and compression damages

目前，高速鐵路隧道的運行安全日益受到人們重視，尤其是在有較大彎道和坡道的盾構(gòu)隧道中，高速列車脫軌撞擊可能對盾構(gòu)隧道拼裝式柔性管片襯砌造成嚴(yán)重破壞，危及旅客人身安全．關(guān)于高速列車撞擊問題，國內(nèi)外已進(jìn)行了一些研究，主要集中在兩方面．一是列車自身的撞擊力學(xué)特性、列車耐撞性或吸能列車設(shè)計．如S W Kirkpatrick、M Schroeder等對車輛結(jié)構(gòu)的耐撞性能進(jìn)行了數(shù)值分析［1-2］；D Tyrell等對列車與列車撞擊進(jìn)行了模擬試驗［3］；A Sutton介紹了不同車輛的設(shè)計、制造以及車輛碰撞標(biāo)準(zhǔn)［4］；姚松、田紅旗、高廣軍等模擬了薄壁結(jié)構(gòu)車輛吸能部件碰撞問題，分析了列車多體耦合撞擊特性，并將研究成果應(yīng)用于吸能列車［5-6］；晏啟祥、李彬、張蒙等研究了時速200 km／h的高速列車的撞擊荷載，分析了盾構(gòu)隧道二次襯砌對管片襯砌的防護(hù)效果［7-8］．二是列車撞擊防撞結(jié)構(gòu)設(shè)計問題．如向俊、龔凱等通過對高速鐵路無砟軌道橋梁上高速列車脫軌全過程的數(shù)值模擬，得到了防撞墻受力的計算公式［9］；吳彪、王吉英對一座特大橋鋼筋混凝土墻式護(hù)欄翼緣板的橫向防撞強(qiáng)度進(jìn)行了驗算［10］．相關(guān)研究非常有限，且主要針對橋梁結(jié)構(gòu)．

盡管目前橋梁防撞墻設(shè)計模擬了列車撞擊行為，但并未給出列車撞擊荷載，而是對不同編組列車進(jìn)行整車建模，然后實施撞擊．這種方式數(shù)值建模復(fù)雜、運算工作量大，不易被工程設(shè)計人員掌握，可普及程度低．目前，無論是研究列車撞擊作用下列車車輛本身還是橋梁結(jié)構(gòu)的防撞問題，國內(nèi)外都未明確給出列車脫軌撞擊荷載的時程曲線．這一現(xiàn)狀嚴(yán)重制約了列車脫軌撞擊周邊結(jié)構(gòu)的動力學(xué)效應(yīng)研究．因此，通過有限元數(shù)值分析，獲取近似的列車撞擊荷載及其特征曲線，對于分析列車撞擊隧道或橋梁結(jié)構(gòu)具有實用價值和工程意義．

本文利用ABAQUS有限元軟件，建立了簡化的列車-剛性墻撞擊三維模型［11］，探討了不同列車編組、不同撞擊速度和不同撞擊角度下的列車撞擊力時程曲線，以揭示列車撞擊的荷載特性，并將獲得的撞擊力作為列車撞擊荷載作用于雙層襯砌盾構(gòu)隧道，對二次襯砌厚度對外部管片襯砌的防護(hù)效果進(jìn)行了分析．

1 列車撞擊數(shù)值模型

為獲取列車撞擊力時程曲線，用CATIA軟件建立列車的三維有限元模型，并用ABAQUS軟件建立列車-剛性墻撞擊動力分析模型，分析了列車不同編組、不同撞擊速度、不同撞擊角度下的撞擊力．采用Hilber-Hughes-Taylor時間積分法［12］（H．H．T法）求解，其動力平衡方程及其位移和速度分別為：

式中：M、C和K分別為質(zhì)量矩陣、阻尼矩陣和剛度矩陣；P為列車撞擊力；u、˙u和¨u分別為位移、速度和加速度；β和γ為由權(quán)重因子α決定的參數(shù)．

為提高計算效率，對列車模型進(jìn)行適當(dāng)簡化，不計小倒角和對撞擊無主導(dǎo)影響的部件，用鋁合金型材薄殼單元模擬車體，采用嵌入泡沫芯材和玻璃鋼等型材實現(xiàn)車體剛度近似等效和質(zhì)量分布相似；車廂之間的車鉤用非線性彈簧模擬，彈簧剛度和阻尼參數(shù)根據(jù)國產(chǎn)密接式車鉤緩沖裝置力學(xué)性能曲線確定［13］，分別取2 000 kN／m和40 kN／（s／m）．

圖1為典型列車編組模型，由1節(jié)機(jī)車和8節(jié)車輛組成．機(jī)車長25 700 mm，寬3 400 mm，高3 700 mm，車鉤中心線間距25 000 mm；中間車輛長24 500mm，寬3 400mm，高3 700mm，車鉤中心間距25 400 mm．列車編組的材料參數(shù)見表1．

圖1 列車編組三維模型Fig．1 A three-dimensional model of train formation

表1 列車車體材料參數(shù)Tab．1 Material parameters of train body

2 撞擊力時程曲線

為獲得不同列車編組、不同撞擊速度、不同撞擊角度下的列車撞擊力時程曲線，取7、9、11、13、15和17節(jié)共6種列車編組，50、100、150、200、250和300 km／h共6種撞擊速度，5．0°、7．5°、10．0°、12．5°、15．0°和17．5°共6種撞擊角度，總計216種組合工況進(jìn)行分析．本文中，撞擊角度是指列車行駛撞擊方向與被撞平面的夾角．

2．1 不同列車編組時的撞擊力時程曲線

圖2為不同編組列車以300 km／h的速度、以17．5°角撞擊剛性墻的撞擊力時程曲線．可見，6種編組的撞擊力時程曲線趨勢相同，數(shù)值接近．相同撞擊速度和撞擊角度下6種編組的撞擊力時程曲線具有相同特點，說明撞擊力時程曲線主要與前幾節(jié)車輛有關(guān)．原因是撞擊過程時間極短，后續(xù)車輛的慣性力通過車鉤傳遞到機(jī)車撞擊部位時已經(jīng)是撞擊過程的后段，難以形成疊加效應(yīng)．

圖2 6種編組的列車撞擊力時程曲線Fig．2 Time-history curves of impact force for six kinds of train formations

2．2 不同撞擊速度和角度時的撞擊力時程曲線

圖3為9節(jié)編組列車分別在5．0°、7．5°、 10．0°、12．5°、15．0°和17．5°撞擊角度下，不同撞擊速度時的撞擊力時程曲線．

圖3 撞擊速度和角度不同時的列車撞擊力時程曲線Fig．3 Time-history curves of impact force under different impact velocities and impact angles

可見，撞擊力最大值與撞擊速度緊密相關(guān)，撞擊速度越大，撞擊力越大；同一撞擊速度下，撞擊角度越大，撞擊力越大，以300 km／h、17．5°角撞擊時的撞擊力最大值可達(dá)40 MN；撞擊力作用時間與撞擊角度密切相關(guān)，撞擊角度在5．0°左右時，撞擊作用時間幾乎不隨撞擊速度改變，都在25 ms左右；但當(dāng)撞擊角度達(dá)到7．5°后，撞擊力作用時間隨撞擊速度增大而增大，以300 km／h的速度撞擊時，撞擊力作用時間可達(dá)37．5 ms；撞擊速度大于等于200 km／h時，撞擊力曲線震蕩明顯，呈現(xiàn)出明顯的多峰特性．

2．3 撞擊力特征曲線及其擬合公式

從圖3可見，撞擊力時程曲線呈現(xiàn)2類特征曲線的特點．第1類特征曲線：撞擊力在撞擊發(fā)生瞬間急劇增大并達(dá)到最大值，隨后急劇下降，然后上升并在若干峰值附近持續(xù)震蕩，最后逐漸減小直至消失，其中，較明顯的峰值有3個，見圖4．圖3中的大多數(shù)曲線都呈現(xiàn)這一特征．第2類特征曲線：撞擊發(fā)生瞬間撞擊力緩慢增大，然后略有下降，隨之再緩慢增大，撞擊后若干毫秒才達(dá)到最大值，最后緩慢降低直至消失，如圖5．這種情況一般在50 km／h等低速撞擊條件下出現(xiàn)．

圖4 第1類撞擊力特征曲線Fig．4 The first kind of characteristic curve of impact force

圖4 的撞擊力特征曲線中，撞擊發(fā)生瞬間撞擊力最大，達(dá)17．30 MN，緊接著開始減小，并在低位震蕩，約11 ms后重新上升到一個高位震蕩，震蕩時間約5．5 ms．原因是后續(xù)多節(jié)車廂的慣性撞擊力相繼通過自動車鉤連接傳向車頭的緣故．撞擊力在撞擊發(fā)生后約32 ms消失．

圖5的撞擊力特征曲線中，撞擊瞬間后緩慢增大，3ms后增大到3．23MN，然后略有下降，隨后再次緩慢增大到最大值4．56 MN，最后緩慢降低直至消失．

圖5 第2類撞擊力特征曲線Fig．5 The second kind of characteristic curve of impact force

由第1類特征曲線可知，在列車撞擊過程中，將會出現(xiàn)多個峰值，但較為明顯的峰值有3個（圖4）．采用高斯多峰擬合方法［14］，對撞擊力曲線的3個顯著峰值進(jìn)行分峰擬合，可得到3個積分強(qiáng)度不同的擬合峰，見圖6．撞擊力擬合公式（4）中共涉及10個擬合參數(shù)：

式中：P為擬合的撞擊力；P0為多峰擬合基線；A1、A2和A3分別為3個擬合峰的單峰面積即積分強(qiáng)度；w1、w2和w3分別為3個擬合峰的半峰寬；tc1、tc2和tc3分別為3個擬合峰的峰頂位置．

圖6 第1類撞擊力特征曲線的高斯多峰擬合Fig．6 Gauss multiple peak fitting of the first kind of characteristic curve

為便于應(yīng)用，針對圖3中撞擊力時程曲線，以目前國內(nèi)城際客運專線最低速度200 km／h作為基準(zhǔn)，考慮200、250和300 km／h三種速度，并選擇10．0°、12．5°和15．0°三種撞擊角度，9種情況下擬合公式中的參數(shù)見表2．采用擬合公式可避免通過復(fù)雜的列車撞擊模型獲取列車撞擊荷載．

表2 高斯多峰擬合參數(shù)Tab．2 Parameters of Gauss multiple peak fitting

3 列車撞擊荷載在二次襯砌動力響應(yīng)分析中的應(yīng)用

在獲得不同列車速度和撞擊角度下列車近似撞擊力的基礎(chǔ)上，選取9節(jié)編組、12．5°撞擊角度、200 km／h撞擊速度的列車撞擊力時程曲線，針對國內(nèi)某水下雙層襯砌盾構(gòu)隧道，用ABAQUS軟件分析管片襯砌在5種二次襯砌厚度下的應(yīng)力、速度、加速度、拉壓損傷動力響應(yīng)．隧道內(nèi)凈空半徑4．37 m，管片襯砌厚度0．48 m，二次襯砌厚度分別為200、250、300、350和400 mm．圖7為二次襯砌厚0．30 m時的典型襯砌構(gòu)造．

隧道圍巖為第四系上更新統(tǒng)Q3和全新統(tǒng)Q4沉積物，主要由砂性土、風(fēng)化泥質(zhì)粉砂等組成，用莫爾-庫倫彈塑性模型模擬．管片襯砌和二次襯砌用塑性損傷模型模擬［15］．?dāng)?shù)值分析模型長、寬、高分別為90、60和60 m，邊界采用彈簧-阻尼器人工邊界．管片襯砌環(huán)間接頭通過弱化該處混凝土襯砌的方式近似模擬．砂性土、風(fēng)化泥質(zhì)粉砂、管片襯砌和二次襯砌的物理力學(xué)參數(shù)見表3．

將以上選取的撞擊力時程曲線以近似平均面力的形式施加在二次襯砌上．為分析方便，選取管片襯砌內(nèi)側(cè)一系列點作為數(shù)值分析監(jiān)測點．以二次襯砌撞擊區(qū)域中心橫向水平對應(yīng)的管片襯砌點為基點，水平縱向前后兩側(cè)各等距離取7個數(shù)值分析監(jiān)測點，間距為0．5 m，15個監(jiān)測點編號依次為1～15，撞擊中心點對應(yīng)8號監(jiān)測點（圖8）．列車從左至右行駛撞擊在隧道上．

3．1 外部管片襯砌的動力響應(yīng)

圖7 典型襯砌構(gòu)造Fig．7 Structure diagram of typical lining

表3 材料參數(shù)Tab．3 Material parameters

圖8 數(shù)值分析監(jiān)測點Fig．8 Monitoring points for numerical analysis

圖9為不同二次襯砌厚度下管片襯砌各監(jiān)測點的Mises應(yīng)力峰值、速度峰值和加速度峰值（峰值指整個撞擊時程中監(jiān)測點的最大值）．可見，Mises應(yīng)力峰值出現(xiàn)在撞擊位置（8號監(jiān)測點），而速度、加速度響應(yīng)峰值出現(xiàn)在10號監(jiān)測點，這與列車撞擊力作用方向有關(guān)；二次襯砌厚度對管片襯砌動力響應(yīng)的影響較大，管片襯砌的Mises應(yīng)力隨厚度增大明顯下降，表明二次襯砌對分散撞擊荷載、改善管片襯砌受力、避免應(yīng)力集中具有重要作用．

各監(jiān)測點的速度峰值和加速度峰值同樣隨二次襯砌厚度增大而減小，撞擊點附近減小效果最明顯．400 mm厚二次襯砌相較200 mm厚二次襯砌，外部管片襯砌的Mises應(yīng)力、速度和加速度峰值分別減小35．73%、18．68%和14．52%．

3．2 外部管片襯砌的損傷特性

圖10為二次襯砌厚度為300 mm時管片襯砌的受壓、受拉損傷云圖．可見，管片襯砌受拉損傷主要分布在撞擊點右側(cè)，撞擊區(qū)域上下兩側(cè)損傷值最大，但面積較?。渌我r砌厚度下管片襯砌的撞擊損傷具有相近的特點，不再贅述．

圖10 二次襯砌厚300 mm時管片襯砌的拉、壓損傷云圖Fig．10 Tension and compression damage images of segmental lining with 300 mm-thick secondary lining

表4為管片襯砌拉、壓損傷與二次襯砌厚度的關(guān)系．表4表明：管片襯砌受拉損傷顯著大于受壓損傷；隨二次襯砌厚度增大，拉、壓損傷面積減小．與厚度200 mm的二次襯砌相比，厚度400 mm時管片襯砌受拉損傷面積減小35．73%，受壓損傷面積減小23．21%．管片襯砌受拉和受壓損傷最大值隨二次襯砌厚度變化的幅度較?。M管增加二次襯砌厚度對降低拉、壓損傷的效果不明顯，但可明顯減小管片襯砌的拉、壓損傷區(qū)域，因而增大二次襯砌厚度能提高管片襯砌抵抗損傷破壞的能力．

表4 不同厚度管片襯砌損傷最大值與損傷面積Tab．4 Maximum damage and damage area of segment lining vs．secondary lining thickness

4 結(jié) 論

（1）列車斜向撞擊力主要與列車編組、撞擊速度和撞擊角度有關(guān)．但當(dāng)列車編組超過一定數(shù)量后，列車撞擊力與列車編組數(shù)量的相關(guān)性減?。?/p>

（2）列車編組數(shù)量一定時，列車斜向撞擊力最大值隨撞擊速度和撞擊角度增大而增大；撞擊力作用時間與撞擊角度有關(guān)，當(dāng)撞擊角度增大到一定值后，撞擊力作用時間隨撞擊速度增大而增大．

（3）列車斜向撞擊力時程曲線呈現(xiàn)2類特征，第1類特征曲線的撞擊力最大值出現(xiàn)在撞擊發(fā)生瞬間，可用高斯多峰擬合公式擬合；第2類特征曲線的撞擊力最大值出現(xiàn)在撞擊瞬間后數(shù)毫秒．

（4）在雙層襯砌盾構(gòu)隧道中，增加二次襯砌厚度可明顯降低列車撞擊荷載作用下管片襯砌應(yīng)力、速度和加速度等動力響應(yīng)，并可減小管片襯砌的拉、壓損傷區(qū)域，提高管片襯砌的抗損傷能力．

［1］ KIRKPATRICK SW，SCHROEDER M．Evaluation of passenger rail vehicle crashworthiness［J］．International Journal of Crashworthiness，2001，6（1）：95-106．

［2］ KIRKPATRICK SW，MACNEILL R A．Development of computer model for prediction of collision response of a railroad passengers car［C］∥Proceeding of Joint Rail Conference．Washington D．C．：IEEE，2002：9-16．

［3］ TYRELL D，JACOBSEN K，PARENT D．Preparations for a train-to-train impact test of crash-energy management passenger rail equipment［C］∥Proceedings of Joint Rail Conference．Pueblo，Colorado：IEEE，2005：107-116．

［4］ SUTTON A．The development of rail vehicle crash-worthiness［J］．Journal of Rail and Rapid Transit，2002，216（2）：97-108．

［5］ 姚松，田紅旗．車輛吸能部件的薄壁結(jié)構(gòu)碰撞研究［J］．中國鐵道科學(xué)，2001，22（2）：55-60．YAO Song，TIAN Hongqi．Crash research on thin-shelled structure as vehicle energy-absorbing components［J］．China Railway Science，2001，22（2）：55-60．

［6］ 高廣軍，田紅旗，姚松，等．列車多體耦合撞擊分析［J］．中國鐵道科學(xué)，2005，26（4）：93-97．GAO Guangjun，TIAN Hongqi，YAO Song，et al．Multi-car coupling collision analysis［J］．China Railway Science，2005，26（4）：93-97．

［7］ 晏啟祥，李彬，張蒙，等．200 km／h列車脫軌撞擊作用下盾構(gòu)隧道二次襯砌對管片襯砌的防護(hù)效果［J］．中國鐵道科學(xué)，2014，35（6）：70-78．YAN Qixiang，LI Bin，ZHANG Meng，et al．Protective effect of secondary lining of shield tunnel on segment lining under derailment impact at the speed of 200 km／h［J］．China Railway Science，2014，35（6）：70-78．

［8］ 張蒙．列車撞擊荷載下水下盾構(gòu)隧道的動力響應(yīng)特性研究［D］．成都：西南交通大學(xué)土木工程學(xué)院，2013．

［9］ 向俊，龔凱，毛建紅，等．高速列車運行安全性與橋梁防撞墻受力分析［J］．鐵道學(xué)報，2011，33（12）：83-87． XIANG Jun，GONG Kai，MIAO Jianhong，et al．Analysis on the running safety of high-speed train and the force of bridge collision-proof wall［J］．Journal of the China Railway Society，2011，33（12）：83-87．

［10］ 吳彪，王吉英．從橋梁翼緣板的橫向防撞強(qiáng)度驗算談邊梁的橫向配筋［J］．中國市政工程，2004（2）：19-20．WU Biao，WANG Jiying．On transverse reinforcement based on computation checking of transverse crash strength of bridge flange slab［J］．China Municipal Engineering，2004（2）：19-20．

［11］ 王懷東，陳秉智，陳一萍，等．城軌列車虛擬碰撞研究［J］．大連交通大學(xué)學(xué)報，2011，32（3）：15-19．WANG Huaidong，CHEN Bingzhi，CHEN Yiping，et al．Study of urban railway vehicle virtual collision simulation［J］．Journal of Dalian Jiaotong University，2011，32（3）：15-19．

［12］ 黃勝，陳衛(wèi)忠，楊建平，等．地下工程地震動力響應(yīng)及抗震研究［J］．巖石力學(xué)與工程學(xué)報，2009，28（3）：483-490．HUANG Sheng，CHEN Weizhong，YANG Jianping，et al．Research on earthquake-induced dynamic responses and aseismic measures for underground engineering［J］．Chinese Journal of Rock Mechanics and Engineering，2009，28（3）：483-490．

［13］ 王娜娜，馬衛(wèi)華．自動車鉤的彈簧剛度和阻尼系數(shù)對列車縱向動力學(xué)的影響［J］．內(nèi)燃機(jī)車，2010（9）：1-3．WANG Nana，MA Weihua．Influence of spring stiffness and damping coefficient of automatic coupler on the longitudinal dynamics of train［J］．Diesel Locomotives，2010（9）：1-3．

［14］ LI Xiaoge，ZHANG Jianhua，LIU Zhaoqing，et al．Determination of critical micelle concentration of binary surfactant mixtures using UV-VIS spectrophotometry and Gaussian fitting method［J］．Global Journal of Physical Chemistry，2011，2（1）：34-38．

［15］ NEDJAR B．Elastoplastic-damage modelling including the gradient of damage：formulation and computational aspects［J］．International Journal of Solids and Structures，2001，38（30／31）：5412-5451．

（中、英文編輯：付國彬）

Numerical Analysis of Train Impact Load with Finite Element Method

YAN Qixiang1， LI Bin1， ZHANG Meng2， HE Chuan1， YANG Wenbo1， GENG Ping1
（1．Key Laboratory of Transportation Tunnel Engineering of Ministry of Education，Southwest Jiaotong University，Chengdu 610031，China；2．China Railway Engineering Consultants Group Co．，Ltd．，Beijing 100020，China）

In order to obtain the train impact load on shield tunnel to investigate the dynamic responses of shield tunnel segment under train impact load，a 3D numerical analysis model for train formation was established to obtain the time-history curves of train impact force at different impact velocities，train formations and oblique impact angles．The maximum train impact force and impact duration as functions of impact velocity and impact angle were researched．Based on the representative time-history curves，the dynamic responses of shield tunnel with different thickness double lining caused by train impact load were investigated．The research result shows that when train formation is certain，the maximum train impact force increases with the increases of impact velocity and impact angle；when impact angle is larger than 7．5°，impact duration extends with the increase of impact velocity．From different times at which the maximum train impact force appears，the time-history curves of train impact force can be divided into two kinds of characteristic curve，and the first kind of characteristic curve with the maximum impact force appearing at the impact moment is in accordance with the Gauss multiple peak fitting formula，and 10 parameters can be used to achieve its approximate fitting．To increase secondary lining thickness can effectively reduce the dynamic responses of external segment lining under train impact load such as stress，velocity and acceleration as well as tension and compression damage zones．

U260．11；U451．4

A

0258-2724（2016）01-0001-07

10．3969／j．issn．0258-2724．2016．01．001

2014-12-17

國家自然科學(xué)基金資助項目（51178400，51278425）；教育部新世紀(jì)人才資助項目（NCET-11-0713）

晏啟祥（1971—），男，教授，博士，博士生導(dǎo)師，研究方向為隧道工程，電話：13088020956，E-mail：764365015＠qq．com

晏啟祥，李彬，張蒙，等．列車撞擊荷載的有限元數(shù)值分析［J］．西南交通大學(xué)學(xué)報，2016，51（1）：1-7．