盧毓江

(蒙西華中鐵路股份有限公司 運(yùn)輸管理部，北京 100073)

模型參數(shù)對(duì)軌道列車吸能結(jié)構(gòu)特性的影響

盧毓江

(蒙西華中鐵路股份有限公司 運(yùn)輸管理部，北京 100073)

為深入地研究材料本構(gòu)模型參數(shù)、材料、單元尺寸參數(shù)和網(wǎng)格密度分布等對(duì)車輛碰撞仿真結(jié)果的影響，基于碰撞界面的非線性接觸算法，建立了典型吸能裝置的動(dòng)態(tài)沖擊動(dòng)力學(xué)模型，在進(jìn)行材料動(dòng)態(tài)沖擊試驗(yàn)的前提下，考察了虛擬數(shù)字計(jì)算中涉及到的變形體本構(gòu)模型、應(yīng)變率效應(yīng)等關(guān)鍵技術(shù)對(duì)吸能部件的響應(yīng)影響形式.研究表明，變形體的材料本構(gòu)模型的簡(jiǎn)化差異對(duì)動(dòng)態(tài)沖擊響應(yīng)具有顯著的影響；在結(jié)構(gòu)設(shè)計(jì)相同的的情況下，低碳鋼材料的吸能容量要遠(yuǎn)遠(yuǎn)超過鋁合金材料；而過大的網(wǎng)格尺寸不能準(zhǔn)確描述結(jié)構(gòu)細(xì)微的塑性變形.

軌道列車；吸能裝置；耐碰撞性能

0 引言

耐撞性軌道列車的吸能區(qū)可分為車鉤緩沖器吸能區(qū)、特殊的吸能裝置吸能區(qū)、車端結(jié)構(gòu)吸能區(qū)，低速碰撞或高速碰撞初始階段，鉤緩裝置和吸能裝置成為主要吸收列車的碰撞動(dòng)能的部件，而吸能裝置是車輛專用于車輛碰撞吸能的保障裝置.

車輛碰撞過程中，虛擬仿真計(jì)算涉及到幾個(gè)關(guān)鍵技術(shù)：首先是與靜態(tài)加載作用下的強(qiáng)度分析相比，材料本構(gòu)模型發(fā)生了變化，動(dòng)載作用下的本構(gòu)關(guān)系具有一系列不同的力學(xué)性能，特別地，在沖擊加載條件下，大多數(shù)金屬材料的屈服極限呈現(xiàn)增大的趨勢(shì)，而屈服出現(xiàn)滯后現(xiàn)象，這種現(xiàn)象被稱為材料的應(yīng)變率效應(yīng).目前應(yīng)用在軌道列車車體上的常用材料很多，在相同的結(jié)構(gòu)下，吸能裝置所用材料的不同對(duì)該裝置的吸能特性影響也很大[1- 3]；其次，在車輛碰撞仿真中，由于其大變形，大位移以及大轉(zhuǎn)動(dòng)的特性，數(shù)字模型的單元尺寸和網(wǎng)格密度分布等參數(shù)對(duì)模型的仿真精度和計(jì)算效率有重要影響，因此深入地研究材料本構(gòu)模型參數(shù)、材料、單元尺寸參數(shù)和網(wǎng)格密度分布等對(duì)車輛碰撞結(jié)果的影響顯得十分重要.

鑒于此，本文基于碰撞界面的非線性接觸算法，建立了典型吸能裝置的動(dòng)態(tài)沖擊動(dòng)力學(xué)模型，在進(jìn)行材料動(dòng)態(tài)沖擊試驗(yàn)的前提下，考察虛擬數(shù)字計(jì)算中涉及到的變形體本構(gòu)模型、應(yīng)變率效應(yīng)等關(guān)鍵技術(shù)對(duì)軌道列車典型吸能部件的響應(yīng)影響形式.

1 材料本構(gòu)模型參數(shù)化分析

1.1 變形體本構(gòu)模型的影響分析

吸能裝置在撞擊載荷作用下會(huì)發(fā)生大位移、大轉(zhuǎn)動(dòng)、大變形.這種變形可能是線性彈性變形、非線性彈性變形、塑性變形、粘塑性變形及其組合等.材料在進(jìn)入塑性狀態(tài)后，應(yīng)力-應(yīng)變關(guān)系是非線性的，這種應(yīng)力-應(yīng)變關(guān)系是與加載歷史和作用環(huán)境有關(guān)，致使材料真實(shí)的應(yīng)力-應(yīng)變行為模擬變得非常復(fù)雜，因此有必要抓住主要的響應(yīng)特性，對(duì)真實(shí)的應(yīng)力-應(yīng)變關(guān)系進(jìn)行簡(jiǎn)化，對(duì)簡(jiǎn)化的模型進(jìn)行數(shù)學(xué)分析，比較不同簡(jiǎn)化本構(gòu)模型對(duì)響應(yīng)特性的影響.

本文對(duì)某列車典型吸能元件建立了有限元模型，仿真分析場(chǎng)景為通過對(duì)剛性墻施加位移載荷實(shí)現(xiàn)對(duì)該吸能元件的動(dòng)態(tài)仿真，壓縮速度為1 m/s.選取材料密度為2.7×103kg/m3，泊松比為0.33.其有限元模型如圖1所示.

圖1 專有的吸能元件的有限元模型

本節(jié)主要對(duì)變形體各本構(gòu)模型進(jìn)行比較參數(shù)分析，現(xiàn)對(duì)材料參數(shù)分別作如下設(shè)定：

(1)理想彈塑性模型

材料理想的彈性-塑性模型的應(yīng)力-應(yīng)變曲線如圖2所示，其線段OA是彈性階段，AB為塑性階段.

圖2 理想的彈性-塑性模型的應(yīng)力-應(yīng)變曲線

將屈服點(diǎn)的屈服應(yīng)力設(shè)為125 MPa，OA是彈性階段的彈性模量E設(shè)為69 000 MPa.

(2)彈性-線性強(qiáng)化模型

彈性-線性強(qiáng)化模型的應(yīng)力-應(yīng)變曲線如圖3所示，其線段OA是彈性階段，AB為塑性強(qiáng)化階段.

圖3 彈性-線性強(qiáng)化模型的應(yīng)力-應(yīng)變曲線

將屈服點(diǎn)的屈服應(yīng)力設(shè)為125 MPa，OA是彈性階段的彈性模量設(shè)為69 000 MPa，AB是塑性強(qiáng)化階段的切線模型設(shè)為1 500 MPa.

(3)彈性-兩段線性強(qiáng)化模型

彈性-線性強(qiáng)化模型的應(yīng)力-應(yīng)變曲線如圖4所示，其主要分為三條折線階段，線段OA是彈性階段，AB、BC為塑性強(qiáng)化階段.

將屈服點(diǎn)的屈服應(yīng)力設(shè)為125 MPa，OA是彈性階段的彈性模量設(shè)為69 000 MPa，AB是塑性強(qiáng)化第一階段的切線模型 設(shè)為1 750 MPa，BC是塑性強(qiáng)化第二階段的切線模型設(shè)為500 MPa.

圖4 彈性-線性強(qiáng)化模型的應(yīng)力-應(yīng)變曲線

(4)彈性-多段線性強(qiáng)化模型

彈性-線性強(qiáng)化模型的應(yīng)力-應(yīng)變曲線如圖5所示，其主要分為三條折線階段，線段OA是彈性階段，AC為多段塑性強(qiáng)化階段.

圖5 彈性-線性強(qiáng)化模型的應(yīng)力-應(yīng)變曲線

將屈服點(diǎn)的屈服應(yīng)力設(shè)為125 MPa，OA是彈性階段的彈性模量設(shè)為69 000 MPa，AC是多段塑性強(qiáng)化階段，通過多段折線最接近的逼近與真實(shí)的應(yīng)力-應(yīng)變曲線.

如圖6，對(duì)比分析四種材料模型的計(jì)算結(jié)果，可以看出理想彈塑性模型內(nèi)能吸收遠(yuǎn)遠(yuǎn)小于其他三種方案，其吸收的能量為88.50 kJ，該方案在能量吸收方面產(chǎn)生的誤差較大，而其他三種方案其能量吸收相當(dāng)，其中彈性-線性強(qiáng)化模型吸收的能量為170.41 kJ，彈性-兩段線性強(qiáng)化模型吸收的能量為158.20 kJ，彈性-多段線性強(qiáng)化模型吸收的能量為166.20 kJ，通過誤差分析可以得知，以彈性-多段線性強(qiáng)化模為基準(zhǔn)，彈性-線性強(qiáng)化模型的誤差為-2.50%，彈性-兩段線性強(qiáng)化模型的誤差為4.81%，可見其吸收能量的誤差率都低于5%.

圖6 能量吸收能力對(duì)比

由圖7可以看出對(duì)于不同的本構(gòu)模型方案，其界面力變化也有較大的差異，其中最為明顯的是理想彈塑性模型，在整體上都低于其他三個(gè)方案，同時(shí)與彈性-多段線性強(qiáng)化模型相比，彈性-兩段線性強(qiáng)化模型較為接近實(shí)際情況，而彈性-線性強(qiáng)化模型相差的較大.

圖7 界面力情況對(duì)比

根據(jù)以上對(duì)比可知，變形體本構(gòu)模型不同，其變形模式、能量吸收能力、界面力大小都有差異，與彈性-多段線性強(qiáng)化模型相比，彈性-兩段線性強(qiáng)化模型和彈性-線性強(qiáng)化模型較為接近實(shí)際情況，其能量吸收誤差率在5%以內(nèi)，從整體響應(yīng)方面來看，彈性-兩段線性強(qiáng)化模型更為真實(shí)的反映實(shí)際情況；特別地，對(duì)于理想彈塑性模型雖然其造成的誤差較大，但是從這次實(shí)例仿真中可以看出，如果實(shí)際材料的塑性屈服面是平穩(wěn)的，結(jié)構(gòu)的吸能和變形模式響應(yīng)是較為理想的.

1.2 基于實(shí)驗(yàn)數(shù)據(jù)的應(yīng)變率影響分析

車輛在碰撞過程中，結(jié)構(gòu)要承受高速的沖擊影響，而在動(dòng)態(tài)沖擊時(shí)隨著沖擊速度的變化，材料的應(yīng)變率也隨之改變，實(shí)踐表明，應(yīng)變率在10-6～10-2s-1內(nèi)，金屬材料的力學(xué)性能沒有發(fā)生明顯的改變，當(dāng)應(yīng)變率大于10-2s-1時(shí)，金屬的力學(xué)性能將發(fā)生顯著地改變，這就必須考慮由于應(yīng)變率變化導(dǎo)致材料力學(xué)性能改變以及帶來的一系列結(jié)構(gòu)性能的變化.文獻(xiàn)[5]指出，在有限元法應(yīng)用于車輛碰撞過程當(dāng)中，當(dāng)在高速撞擊時(shí)，必須考慮由于應(yīng)變率效應(yīng)所帶來的影響.

為比較材料應(yīng)變率對(duì)結(jié)構(gòu)碰撞響應(yīng)的影響，本文首先對(duì)鋁合金材料5083-H111進(jìn)行了材料動(dòng)態(tài)沖擊試驗(yàn)，通過沖擊加載，獲取鋁合金沖擊載荷下的應(yīng)力-應(yīng)變曲線，揭示不同應(yīng)變率加載下的鋁合金材料的沖擊力學(xué)性能，為以后的應(yīng)用研究打下堅(jiān)實(shí)基礎(chǔ).

本試驗(yàn)所選材料為5083-H111，其化學(xué)成分為，硅：0.4%、鐵：0.4%、銅：0.10%、錳：0.3～1.0%、鎂：4.4%～4.9%、鉻：0.05%～0.25%、鋅：0.25%、鈦：0.15%，其余為A1.試驗(yàn)臺(tái)如圖8所示，試驗(yàn)儀要采用分離式霍普金森壓桿系統(tǒng)(SHPB)進(jìn)行，其組成主要有：炮管、子彈、入射桿、透射桿、吸收桿、超動(dòng)態(tài)應(yīng)變儀和計(jì)算機(jī)等，通過對(duì)氣缸氣壓的控制，來達(dá)到控制子彈不同沖擊速度進(jìn)而控制應(yīng)變率的目的，由一維應(yīng)力波理論，以及超動(dòng)態(tài)應(yīng)變儀獲得輸入桿和輸出桿上的測(cè)試數(shù)據(jù)，根據(jù)這些波形數(shù)據(jù)，即可獲得試樣的平均應(yīng)變率、平均應(yīng)變和試件兩端面應(yīng)力的平均值.

圖8 SHPB實(shí)驗(yàn)系統(tǒng)

通過動(dòng)態(tài)沖擊試驗(yàn)獲得了5083-H111 試樣在加載應(yīng)變率分別為1 943、2 215 和3 317 s-1的應(yīng)力-應(yīng)變曲線，如圖9所示.從圖中可看出，隨著應(yīng)變率增大流動(dòng)應(yīng)力水平也相應(yīng)變大， 即應(yīng)變率對(duì)流動(dòng)應(yīng)力遞增的影響，同時(shí)，最大塑性應(yīng)變值和最大應(yīng)力也同應(yīng)變率成正比規(guī)律變化.

圖9 不同應(yīng)變率下的應(yīng)力-應(yīng)變曲線

將該試驗(yàn)數(shù)據(jù)導(dǎo)入到LS-DYNA材料模型，對(duì)考慮應(yīng)變率情況下的結(jié)構(gòu)碰撞響應(yīng)與不考慮應(yīng)變率情況下的結(jié)構(gòu)碰撞響應(yīng)進(jìn)行比較，計(jì)算工況為吸能裝置以20 m/s的速度撞擊剛性墻，其仿真結(jié)果如圖10所示，可以看出，考慮應(yīng)變率和未考慮應(yīng)變率的吸能裝置能量吸收和界面力值有差異，在吸能能力方面，考慮應(yīng)變率的吸收311.76 kJ，未考慮應(yīng)變率的吸收243.75 kJ，誤差率在27.90%；在界面力方面，初始峰值力相差較小，但是從界面力-時(shí)間特性曲線來看，其曲線的變化規(guī)律誤差也較大.

(a)能量吸收能力對(duì)比

(b)界面力對(duì)比

為了后續(xù)計(jì)算的正確有效性，根據(jù)文獻(xiàn)[5- 6]數(shù)據(jù)，本文增加了一個(gè)低碳鋼材料的虛擬仿真試驗(yàn)，以考察低碳鋼材料的應(yīng)變率效應(yīng)對(duì)結(jié)構(gòu)的碰撞響應(yīng)特性的影響.根據(jù)文獻(xiàn)數(shù)據(jù)，本文計(jì)算應(yīng)變率模型應(yīng)用Cowper-Symonds模型，模型參數(shù)取C=40.4、P=5，碰撞場(chǎng)景采用吸能裝置以20 m/s的速度撞擊剛性墻，其能量吸收及界面力值計(jì)算結(jié)果如圖11所示.

(a)能量吸收能力對(duì)比

(b)界面力對(duì)比

從圖11中可以看出，對(duì)于低碳鋼材料，考慮應(yīng)變率和未考慮應(yīng)變率的吸能裝置能量吸收和界面力值有較大差異，在吸能能力方面，考慮應(yīng)變率的吸收816.07 kJ，未考慮應(yīng)變率的吸收400.18 kJ，誤差率在103.92%；在界面力方面，初始峰值力和界面力-時(shí)間特性曲線相差都較大，其曲線的變化規(guī)律誤差也較大.

因此，方案試驗(yàn)和計(jì)算仿真方面表明，材料應(yīng)變率效應(yīng)對(duì)結(jié)構(gòu)的碰撞響應(yīng)特性有顯著的影響.

2 材料型號(hào)變化的影響分析

地鐵車輛用材料一般包括低碳鋼材料、不銹鋼材料和鋁合金材料，不同材料的力學(xué)特性也相應(yīng)有不同因此在車輛發(fā)生碰撞過程中，其結(jié)構(gòu)的碰撞響應(yīng)也會(huì)有不同的變化規(guī)律，本節(jié)主要對(duì)低碳鋼材料和鋁合金材料進(jìn)行了對(duì)比分析，在結(jié)構(gòu)尺寸不變的情況下，探討不同材料特性對(duì)結(jié)構(gòu)碰撞響應(yīng)特性的影響.

本節(jié)利用上節(jié)的仿真模型，通過對(duì)剛性墻施加位移載荷實(shí)現(xiàn)對(duì)該吸能元件的動(dòng)態(tài)仿真，壓縮速度為20 m/s.選取低碳鋼材料密度為7.85×103kg/m3，泊松比為0.3；鋁合金材料密度為2.7×103kg/m3，泊松比為0.33.

低碳鋼材料和鋁合金材料的吸收能量情況和界面力特性曲線如圖12所示，從圖可以看出，低碳鋼材料和鋁合金材料，在相同結(jié)構(gòu)下，其吸能能力和界面力值具有較大的差異，在相同的結(jié)構(gòu)條件下，低碳鋼材料的吸能容量要遠(yuǎn)遠(yuǎn)超過鋁合金材料，這與文獻(xiàn)[7]結(jié)果是一致的.

(a)能量吸收能力對(duì)比

(b)界面力況對(duì)比

3 單元尺寸的影響分析

運(yùn)用LS-DYNA進(jìn)行有限計(jì)算仿真時(shí)，從各種參數(shù)在計(jì)算時(shí)間所占的比例中可以看出，其所消耗的計(jì)算時(shí)長(zhǎng)主要與參與計(jì)算單元的數(shù)據(jù)傳遞和接觸計(jì)算有關(guān)系，而單元尺寸的大小不僅與模型所占有的單元數(shù)量有聯(lián)系，同時(shí)在計(jì)算時(shí)的穩(wěn)定性方面也有較大影響.例如，四邊形殼單元的極限時(shí)間步長(zhǎng)的計(jì)算公式為：

(1)

其中,le是單元e的特征長(zhǎng)度，ce是單元e的當(dāng)前波速.

因此，為保障車輛碰撞計(jì)算的準(zhǔn)確性，應(yīng)該考慮單元尺寸對(duì)碰撞響應(yīng)的影響.圖13中所示為設(shè)計(jì)的吸能裝置與剛性墻的碰撞模擬仿真圖，參照模型的結(jié)構(gòu)總體尺寸，分別選取了3、5、7、12 mm特征尺寸的4節(jié)點(diǎn)殼單元對(duì)仿真模型進(jìn)行了網(wǎng)格劃分，其有限元單元數(shù)分別為138 978、52 389、26 082、10 498，四種不同單元尺寸劃分的吸能裝置在碰撞最后時(shí)刻的仿真結(jié)果如圖14所示.

圖13 不同尺寸影響因素的仿真模型

(a) 3 mm(b) 5 mm

(c) 7 mm(d)12 mm

圖14 不同尺寸仿真變形模式比較

分析和比較圖14的仿真結(jié)果可以看出：特征尺寸為12 mm的吸能裝置變形模式與其他三種有較大的差異，吸能裝置在局部失穩(wěn)后，未有產(chǎn)生有效的壓塌失效過程，從變形模式方面來看，該尺寸下也為捕捉到吸能裝置在碰撞中產(chǎn)生的小“皺褶”，單元的變化類似“鋸齒形”；特征尺寸為3 mm的變形模式類似于“皺褶”，能較好地反映實(shí)際屈曲變形過程，這與文獻(xiàn)的實(shí)驗(yàn)結(jié)果也吻合一致[8- 9].

因此，通過本節(jié)的分析，單元較小越容易捕捉到結(jié)構(gòu)碰撞中產(chǎn)生的小“皺褶”，而粗網(wǎng)格單元?jiǎng)傂暂^強(qiáng)，在對(duì)車輛進(jìn)行碰撞模擬時(shí)，應(yīng)充分考慮單元尺寸的細(xì)化對(duì)計(jì)算精度和仿真時(shí)間的合理分配.

4 結(jié)論

本文基于碰撞界面的非線性接觸算法，建立了典型吸能裝置的動(dòng)態(tài)沖擊動(dòng)力學(xué)模型，在進(jìn)行材料動(dòng)態(tài)沖擊試驗(yàn)的前提下，考察了虛擬數(shù)字計(jì)算中涉及到的變形體本構(gòu)模型、應(yīng)變率效應(yīng)等關(guān)鍵技術(shù)對(duì)吸能部件的響應(yīng)影響形式.研究表明，變形體的材料本構(gòu)模型的簡(jiǎn)化差異對(duì)動(dòng)態(tài)沖擊響應(yīng)具有顯著的影響，比較分析得出考慮應(yīng)變率的彈性-兩段線性強(qiáng)化模型更能反映實(shí)際；在結(jié)構(gòu)設(shè)計(jì)相同的的情況下，低碳鋼材料的吸能容量要遠(yuǎn)遠(yuǎn)超過鋁合金材料；而過大的網(wǎng)格尺寸不能準(zhǔn)確描述結(jié)構(gòu)細(xì)微的塑性變形，因此需要對(duì)計(jì)算的精度和速度進(jìn)行合理的分配.

[1]張志新．高速列車耐撞性結(jié)構(gòu)及安全性研究[D]．成都：西南交通大學(xué)，2012．

[2]劉艷文．軌道客車碰撞被動(dòng)安全性研究[D]．成都：西南交通大學(xué)，2013．

[3]王晉樂．電力機(jī)車車體結(jié)構(gòu)耐撞性設(shè)計(jì)[D]．成都：西南交通大學(xué)，2014．

[4]鐘志華．汽車耐撞性分析的有限元法[J]．汽車工程，1994，16(1)：1- 6．

[5]王成國(guó)，成艾國(guó)，高暉,等.材料應(yīng)變率對(duì)汽車碰撞性能影響的研究[J]．汽車工程， 2010，32(6)：482- 46．

[6]孟廣成，邰永剛．材料應(yīng)變率對(duì)護(hù)欄碰撞仿真結(jié)果影響研究[J]．交通工程，2012,6(6)：393- 395.

[7]柳艷杰．汽車低速碰撞吸能部件的抗撞性能研究[D]．哈爾濱：哈爾濱工程大學(xué)，2012．

[8]葉映臺(tái)．轎車前縱梁耐撞性仿真和結(jié)構(gòu)優(yōu)化的研究[D]．長(zhǎng)沙：湖南大學(xué)，2010.

[9]胡玉梅．汽車正面碰撞設(shè)計(jì)分析技術(shù)及應(yīng)用研究[D]．重慶：重慶大學(xué)，2002.

Numerical Simulation Study on Effects of Model Parameter on Railway Vehicle Energy-Absorbing Device

LU Yujiang

(Transportation Management Department，Mengxi-Huazhong Railway Co.，Ltd，Beijing 100073，China)

In order to study the effects of model parameters in crash simulation of energy-absorbing device, a typical FEM model based on nonlinear contact algorithm was built.By conducting impact test and considering the key parameter such as constitutive model and strain rate,it is found that material model has obvious influence on crash response.The structure of low-carbon steel can absorb more energy than the one of aluminium alloy in the same crash condition,and oversized mesh cannot simulate the plastic deformation accurately.

railway vehicle；energy-absorbing device；crash worthiness

1673- 9590(2017)04- 0045- 06

2017- 01- 09

盧毓江(1989-)，男，工程師，碩士，主要從事車輛工程相關(guān)研究E-mail:429076304@qq.com.

A