王 蕊， 張曉云， 馮 浩, 陳建國

(1.上海交通大學(xué) 機(jī)械與動力學(xué)院,上海 200240； 2. 司法部 司法鑒定科學(xué)技術(shù)研究所,上海 200063)

交通事故不僅給生活、生產(chǎn)帶來較大危害，且對社會穩(wěn)定、發(fā)展亦造成較大負(fù)面影響。因此，對交通事故過程分析及人體傷害防范的數(shù)值仿真技術(shù)廣受關(guān)注。目前，交通事故數(shù)值仿真技術(shù)研究主要集中于小客車碰撞事故，忽視了大客車交通事故及安全設(shè)計(jì)需求。實(shí)際上我國大、中型客車所致事故數(shù)量及死亡人數(shù)占有較大比例。調(diào)查表明，造成“群死群傷”的惡性重大交通事故多與大客車碰撞、翻滾等事故形態(tài)有關(guān)。如，2011年9月14日上海發(fā)生的大客車側(cè)翻事故，死傷24人；2011年10月7日天津發(fā)生相同事故，死亡35人。兩起事故中大客車均在側(cè)翻后與路邊護(hù)欄發(fā)生碰撞，故而傷亡慘重。因此，針對大客車碰撞事故，結(jié)合車身及車內(nèi)變形等重要信息，研究適用事故仿真及防范方法具有重要的理論意義與實(shí)用價值。

本文類比此兩起交通事故，以某12 m長大客車車身骨架為對象，在HYPERMESH中模擬事故發(fā)生環(huán)境用LS-DYNA進(jìn)行運(yùn)算，后處理結(jié)果在HYPERVIEW中觀察仿真大客車側(cè)翻碰撞護(hù)欄并向前滑行過程。通過改進(jìn)材料及設(shè)計(jì)，探索減小事故損失方法。

1 側(cè)翻事故有限元模型建立

1.1 側(cè)翻事故客車模型前處理

大客車側(cè)翻事故研究對象主要包括大客車及碰撞護(hù)欄兩大部分，整體碰撞事故環(huán)境為：

(1)以大客車車身骨架及車架為板殼單元，底盤車輪、軸、發(fā)動機(jī)為實(shí)體單元，前后6個空氣彈簧懸架及整體橋結(jié)構(gòu)的導(dǎo)向桿系分別為剛性單元、梁單元，并忽略蒙皮、玻璃、內(nèi)飾板等非承載件，見圖1。板殼單元基準(zhǔn)面為實(shí)際結(jié)構(gòu)中面，結(jié)構(gòu)間用共節(jié)點(diǎn)、剛性連接、點(diǎn)焊等模擬連接關(guān)系[1-7]。施加質(zhì)量單元模擬大客車重力，使模型總質(zhì)量、重心位置與實(shí)車一致。

(2)據(jù)文獻(xiàn)[8]規(guī)定，選單坡型路側(cè)混凝土護(hù)欄(防撞等級為SB級)，尺寸參數(shù)按標(biāo)準(zhǔn)確定，建立橫截面為等腰梯形的護(hù)欄幾何模型，并劃分為合適大小六面體實(shí)體單元，置于大客車車輪同一水平面處?；炷磷o(hù)欄為剛性護(hù)欄，用有限元建模護(hù)欄可對大客車碰撞護(hù)欄事故充分評價[9]。

(3)為方便大客車旋轉(zhuǎn)定義，在大客車輪胎下方建立翻滾平臺。翻滾平臺側(cè)邊與相鄰護(hù)欄底邊間距2597 mm，此為高度仿真大客車側(cè)窗處與護(hù)欄碰撞事故，可在一定程度內(nèi)變動。

(4)文獻(xiàn)[10]對大客車內(nèi)生存空間位置與尺寸有嚴(yán)格規(guī)定。而生存空間為可保證乘員生命安全、衡量車體合格與否的重要標(biāo)準(zhǔn)。據(jù)生存空間形狀的不變性，為簡化計(jì)算，略去建立此部分模型，而改用在仿真結(jié)果圖中進(jìn)行繪制示意方式(圖中尺寸、位置嚴(yán)格按建模比例縮放獲得)。

(5)側(cè)翻事故大客車有限元模型見圖1、圖2。共有單元683907個，節(jié)點(diǎn)695298個。其中三角形單元數(shù)量9719，占1.42%，滿足小于5%要求。

圖1 大客車有限元模型

1.2 材料參數(shù)設(shè)置

車體模型中主要采用兩種材料：①M(fèi)AT20剛體材料即大客車輪胎、發(fā)動機(jī)及增加的護(hù)欄等不變形部件；②MAT24分段線性塑性材料即車身各變形部件。底盤部分用20#碳素鋼，車身上部用Q235[6]，改進(jìn)材料Q345，其屬性見表1。參照實(shí)際事故中混凝土護(hù)欄基本不變形事實(shí)(圖3)，為提高運(yùn)算效率，護(hù)欄模型材料選MAT20剛體材料。

1.3 初始條件及接觸設(shè)置

將翻滾平臺靠近護(hù)欄側(cè)及護(hù)欄全約束。為計(jì)算效率及穩(wěn)定，將翻滾平臺預(yù)先自豎直位置旋轉(zhuǎn)45°作為初始位置；在0～0.1 s使翻滾平臺以5°/s角速度旋轉(zhuǎn)，以保證大客車能翻過臨界位置而不影響計(jì)算結(jié)果。整車模型有沿x方向的初速度10 m/s(36 km/h)。

表1 材料屬性表

圖3 實(shí)際交通事故中護(hù)欄

單面接觸用于單個物體表面的自接觸或與另一物體表面的相互接觸，無需定義接觸面與目標(biāo)面，適用于求解預(yù)先未知接觸表面大變形問題，計(jì)算時精確、快速[5]。因此在汽車碰撞仿真分析中，為防止各部件間發(fā)生相互穿透或部件產(chǎn)生大變形時自身產(chǎn)生穿透現(xiàn)象，采用單面接觸以提高數(shù)值計(jì)算的穩(wěn)定性及效率。碰撞分析時，采用自動接觸方式，可減少工作量[7]。側(cè)翻過程中，車體與地面接觸定義為剛性墻接觸，并定義車身與地面間摩擦數(shù)為0.3。車體與護(hù)欄間采用自動搜索類型接觸[11]。

2 大客車側(cè)翻碰撞護(hù)欄事故仿真分析

2.1 仿真結(jié)果驗(yàn)證

大客車由初始位置受重力作用開始側(cè)翻，0.6 s后與護(hù)欄發(fā)生碰撞。由圖4看出，側(cè)窗立柱首先因拉伸強(qiáng)度不足發(fā)生斷裂，護(hù)欄侵入車內(nèi)乘員的生存空間；側(cè)窗立柱與車頂連接處發(fā)生斷裂，原因可能為兩處間焊接強(qiáng)度不足。而車頂與側(cè)窗立柱構(gòu)成的閉環(huán)結(jié)構(gòu)可抵抗變形，車頂坍塌進(jìn)一步擠壓車內(nèi)乘員的生存空間；大客車因慣性向前滑行，更增加事故的嚴(yán)重性。

圖4 大客車側(cè)翻事故仿真序列

按規(guī)定，生存空間以外其他部分不得侵入大客車的生存空間[10]，但模擬最終損毀情況顯示側(cè)窗立柱及護(hù)欄均侵入生存空間，見圖5。與圖6相比，此類交通事故對乘員的生命安全威脅極大。

圖5 大客車側(cè)翻事故損壞細(xì)節(jié)圖

圖6 實(shí)際事故中大客車生存空間壓縮圖

2.2 能量分析

為檢驗(yàn)仿真精度，需觀察能量是否守恒及非物理能量如沙漏能等是否超過此總能量5%[10]。側(cè)翻過程的能量變化曲線見圖7。由圖7看出，仿真過程能量守恒，沙漏能在側(cè)翻變形中占總能量比未超1.65%，在許可范圍5%以內(nèi)，對分析結(jié)果影響微弱。因此，該仿真計(jì)算結(jié)果可信。

圖7 大客車側(cè)翻事故能量分析圖

2.3 模型材料改進(jìn)仿真分析

將大客車側(cè)窗立柱材料改為Q345，其他同初始設(shè)置，仿真結(jié)果見圖8。與圖5對比看出，在側(cè)翻碰撞護(hù)欄過程中，大客車側(cè)窗護(hù)欄仍出現(xiàn)斷裂進(jìn)而擠壓車內(nèi)乘員生存空間情況，護(hù)欄侵入車內(nèi)高度343.965 mm，較Q235側(cè)窗立柱侵入量小，可減輕事故損失。

圖8 材料改進(jìn)后仿真圖

2.4 模型厚度改進(jìn)仿真分析

將大客車模型側(cè)窗立柱材料由1.5 mm加厚至3.0 mm，其他條件不變，仿真結(jié)果見圖9。由圖9看出，護(hù)欄侵入高度294.961 mm，護(hù)欄、側(cè)窗立柱未侵入生存空間，滿足設(shè)計(jì)規(guī)定，與圖8相比，車頂塌陷較少，損毀程度進(jìn)一步減小。

圖9 側(cè)窗立柱厚度改進(jìn)仿真圖

3 改進(jìn)建議

大客車側(cè)翻時車身上部結(jié)構(gòu)強(qiáng)度是決定事故損失程度的重要因素。GB 13094- 2007規(guī)定自2011年2月1日開始強(qiáng)制執(zhí)行對客車上部結(jié)構(gòu)強(qiáng)度要求，主要技術(shù)要求需滿足GB/T 17578《客車上部結(jié)構(gòu)強(qiáng)度的規(guī)定》，其修訂亦參照歐盟的ECE R663[12]?，F(xiàn)有標(biāo)準(zhǔn)中只對大客車側(cè)壁與地面碰撞進(jìn)行側(cè)翻試驗(yàn)，未考慮側(cè)翻碰撞護(hù)欄更慘重情況。建議大客車設(shè)計(jì)者在標(biāo)準(zhǔn)實(shí)驗(yàn)基礎(chǔ)上考慮此類事故，損失會進(jìn)一步減少。

4 結(jié) 論

(1)由仿真知，增加大客車側(cè)窗立柱材料強(qiáng)度及厚度均可減小損失。

(2)增加側(cè)窗立柱與客車頂部連接強(qiáng)度、側(cè)窗立柱數(shù)量、吸能元件等，可減少護(hù)欄浸入車廂擠壓生命空間的概率。通過面向真實(shí)交通事故數(shù)值仿真分析，可以為大客車結(jié)構(gòu)安全設(shè)計(jì)提供理論參考。

參 考 文 獻(xiàn)

[1]范子杰,馬 林,桂良進(jìn),等.大客車側(cè)翻安全性多目標(biāo)優(yōu)化[J].系統(tǒng)仿真學(xué)報(bào),2012,24(5):109-113.

FAN Zi-jie, MA Lin, GUI Liang-jin, et al. Multi-objective optimization for bus rollover crashworthiness [J]. Journal of System Simulation, 2012, 24(5):109-113.

[2]何漢橋.大客車車身結(jié)構(gòu)安全性仿真研究[D].長沙：湖南大學(xué)機(jī)械與汽車工程學(xué)院，2007.

[3]孫 信,王青春,范子杰.大客車翻滾安全性有限元分析[J].汽車技術(shù),2007,38(8):34-35.

SUN Xin,WANG Qing-chun,FAN Zi-jie, et al. Finite element analysis on rollover crashworthiness of a bus [J]. Automobile technology, 2007, 38(8):34-35.

[4]何漢橋,張維剛.高床大客車側(cè)翻結(jié)構(gòu)安全性仿真研究[J].機(jī)械科學(xué)與技術(shù),2007,26(7):922-925.

HE Han-qiao, ZHANG Wei-gang. Simulation of safety of high-bed bus under side overturning crash[J]. Mechanical Science and Technology for Aerospace Engineering, 2007, 26(7):922-925.

[5]余 躍,周鴻波,童水光,等.大客車側(cè)翻耐撞性的結(jié)構(gòu)設(shè)計(jì)方法[J].浙江大學(xué)學(xué)報(bào)(工學(xué)版), 2011,45(4):714-718.

YU Yue, ZHOU Hong-bo, TONG Shui-guang，et al. Structural design method for rollover crashworthiness of large bus[J]. Journal of Zhejiang University (Engineering Science), 2011, 45(4):714-718.

[6]盧琳兆,吳長風(fēng),丁守松,等.大客車側(cè)翻碰撞仿真分析及改進(jìn)[J].機(jī)電技術(shù),2012,35(1):26-28.

LU Lin-zhao, WU Chang-feng, DING Shou-song, et al. Simulation and improvement of bus rollover [J]. Mechanical and Electrical Technology, 2012,35(1):26-28.

[7]鄧兆祥，胡玉梅，王 攀,等.客車耐撞性結(jié)構(gòu)優(yōu)化設(shè)計(jì)[J]. 機(jī)械工程學(xué)報(bào),2005,41(11):217-220.

DENG Zhao-xiang, HU Yu-mei, WANG Pan, et al. Optimization design of mini car crashworthiness structure[J]. Chinese Journal of Mechanical Engineering,2005,41(11): 217- 220.

[8]JTG D81-2006,公路交通安全設(shè)施設(shè)計(jì)細(xì)則[S].

[9]閆書明.護(hù)欄碰撞中車體加速度敏感性分析[J].振動與沖擊,2011,30(7):1-4.

YAN Shu-ming. Sensitivity analysis of vehicle acceleration during impact to a barrier[J].Journal of Vibration and Shock,2011,30(7):1-4.

[10]Regulation No.66, Uniform technical prescriptions concerning the approval of large passenger vehicles with regard to the strength of their superstructure[S]. Geneva: United Nations Economic Commission for Europe, 2006.

[11]閆書明.有限元仿真方法評價護(hù)欄安全性能的可行性[J].振動與沖擊,2011,30(1): 152-156.

YAN Shu-ming. Feasibility analysis of barrier safety evaluation with finite element simulation method [J].Journal of Vibration and Shock,2011,30(1):152-156.

[12]李仕鋒,丁良旭,曹源文. 客車被動安全研究主要方法及內(nèi)容[J].客車技術(shù)與研究,2012,34(2):5-8.

LI Shi-feng, DING Liang-xu, CAO Yuan-wen. Main methods and contents of coach passive safety research [J]. Bus & Coach Technology and Research, 2012, 34(2):5-8.