張君媛,紀夢雪,王 楠,葉 威
(1.吉林大學(xué) 汽車仿真與控制國家重點實驗室,長春 130022;2.清華大學(xué) 汽車安全與節(jié)能國家重點實驗室,北京 100084;3.北京伊薩科技發(fā)展有限公司 銷售部,北京 100094)
在所有汽車正面碰撞死亡事故中,小偏置(25%重疊)正面碰撞約占正面碰撞總量的四分之一,但在各國安全性法規(guī)及評價規(guī)程中,僅美國公路安全保險協(xié)會(Insurance institute for highway safety,IIHS)從2012年開始對小偏置正面碰撞做評價[1]。
本文統(tǒng)計了2012~2014年間IIHS公布的30款車型正面全寬碰撞、正面40%偏置碰撞以及正面25%重疊率碰撞試驗的得分情況。其中正面25%重疊率碰撞工況的滿分率是最低的。在25%重疊率碰撞試驗中,車體結(jié)構(gòu)完整性的評價得分相對于約束系統(tǒng)與假人運動及假人傷害兩部分得分也最低。這種形式的碰撞對整車的考驗在于前端結(jié)構(gòu)(單側(cè))的吸能能力和側(cè)圍的縱向強度。國內(nèi)外對該工況的研究主要集中在車體結(jié)構(gòu)修改上,通常是根據(jù)某款車在該工況試驗中出現(xiàn)的吸能不足、乘員艙侵入情況,對結(jié)構(gòu)進行修改,再通過試驗或CAE方法進行驗證[2,3]。
本文建立了一種汽車正面25%重疊率碰撞工況下的簡化有限元分析模型。提出一種從乘員艙臨界侵入量角度出發(fā)的車體前端構(gòu)件設(shè)計方法,即先針對原車吸能不足的情況,以侵入量為約束條件,提出前端主要縱向薄壁梁構(gòu)件的目標吸能量,再將其分配到各子構(gòu)件。然后,根據(jù)吸能目標利用薄壁梁理論模型進行薄壁梁結(jié)構(gòu)斷面的快速設(shè)計,使之在結(jié)構(gòu)斷面設(shè)計之初即滿足總吸能量要求,以便最終滿足侵入量的要求。
如圖1所示,IIHS規(guī)定:該工況中正面重疊率為車身寬度的25%,碰撞速度為63.4~65.4 km/h,采用剛性壁障(其正面形狀為一寬1000 mm、高1524 mm的長方形;右端為一半徑150 mm、弧度115°的圓?。慌鲎裁驿摪搴穸葹?8.1 mm),且要求在駕駛員側(cè)放置一個50th百分位混Ⅲ男性假人[1]。
正面25%重疊率碰撞評價方法包括車體結(jié)構(gòu)完整性、約束系統(tǒng)與假人運動及假人傷害3部分。每一部分評價分為“優(yōu)秀”、“良好”、“及格”、“差”4個等級。其中,車體結(jié)構(gòu)完整性評價部分包括10個測量點(以下稱為侵入量監(jiān)測點)侵入量的評價:A柱上、下部,儀表板上、下部,歇腳板,左側(cè)地板,制動踏板,駐車踏板,門檻,轉(zhuǎn)向柱。
圖1 正面25%重疊率碰撞試驗工況Fig.1 25% overlap frontal collision test condition
本文以某乘用車為研究實例,在Pam-Crash/Safe軟件中建立了該車型有限元模型的IIHS正面25%重疊率碰撞工況仿真模擬環(huán)境。圖2為25%重疊率正面碰撞車輛-假人有限元基礎(chǔ)模型,該模型的坐標系與車身坐標系相同。采用車輛主動碰撞固定壁障的方式進行模擬,沖擊速度為64 km/h。圖3為有限元基礎(chǔ)模型在仿真模擬中的變形情況(為觀察乘員艙的侵入情況隱藏了左前輪):A柱上部發(fā)生嚴重的彎折變形,A柱下部(即車門鉸處)侵入量較大,左前門變形嚴重可能會導(dǎo)致車門不能順利開啟,同時前縱梁中段未充分變形,即車體前端構(gòu)件吸能不足導(dǎo)致了嚴重的乘員艙侵入。
圖2 正面25%重疊率碰撞車輛-假人有限元基礎(chǔ)模型Fig.2 Basic vehicle-dummy finite element modelfor 25% overlap frontal collision
圖3 基礎(chǔ)模型正面25%重疊率碰撞變形結(jié)果Fig.3 25% overlap frontal collisional deformationresults of basic model
圖4為有限元基礎(chǔ)模型在該工況下碰撞前、后乘員小腿和足部運動狀態(tài)對比圖。小腿和足部傷害與A柱下部、歇腳板、左側(cè)地板、制動踏板、駐車踏板、門檻(側(cè)向)、下儀表板這7個法規(guī)規(guī)定的侵入量相關(guān)。
圖4 碰撞前后乘員小腿和足部運動狀態(tài)對比Fig.4 Comparison of legs and feet movements ofoccupant before and after collision
分析相關(guān)構(gòu)件侵入量與小腿和足部傷害的關(guān)聯(lián),發(fā)現(xiàn)3條與小腿脛骨軸向力和脛骨橫向彎矩相關(guān)的力的傳遞路徑:第1條是保險杠與壁障接觸產(chǎn)生的力使翼子板和前縱梁向后擠壓,致使車門鉸和門檻處產(chǎn)生變形擠壓到小腿,使小腿受到脛骨橫向力矩;第2條是保險杠與壁障接觸產(chǎn)生的力傳至前縱梁、副車架、發(fā)動機等,使其向后擠壓侵入到前地板/前圍下部(包括左側(cè)地板、歇腳板、制動踏板、駐車踏板等法規(guī)規(guī)定的侵入量評價點),前地板/前圍下部變形將力傳至足部,使小腿受到脛骨軸向力;第3條是保險杠與壁障接觸后受力向后擠壓,導(dǎo)致發(fā)動機侵入乘員艙,使儀表板產(chǎn)生向后的位移,可能會擠壓乘員小腿產(chǎn)生脛骨軸向力。將基礎(chǔ)模型7個監(jiān)測點的侵入量與法規(guī)臨界侵入量[1](指法規(guī)規(guī)定的車體結(jié)構(gòu)評分達到優(yōu)秀的相關(guān)監(jiān)測點的侵入量最大值)進行對比,如圖5所示。該車駐車踏板、門檻及儀表板下部侵入量明顯超過了法規(guī)臨界侵入量。
圖5 乘員艙侵入量與法規(guī)臨界侵入量對比Fig.5 Comparison between intrusion of cabinand critical intrusion of regulations
基礎(chǔ)模型碰撞前總動能為203 kJ,碰撞中總的吸能量為151 kJ。本文將該工況下車身吸能較多的構(gòu)件分為3類:前端薄壁梁構(gòu)件、前端其他構(gòu)件和乘員艙構(gòu)件。這3部分的吸能情況見圖6。分析發(fā)現(xiàn)該工況下主要吸能構(gòu)件依次為:前縱梁、翼子板、輪罩、前車門、副車架、前指梁、保險杠總成等。乘員艙構(gòu)件吸能量之和約占總吸能量的40%,且該工況下乘員艙相當于單側(cè)吸能,再次驗證了車體前端構(gòu)件吸能不足是乘員艙侵入過大的原因。因此后續(xù)設(shè)計中應(yīng)降低乘員艙吸能量,將其分配給前端構(gòu)件。
圖6 三類構(gòu)件吸能量-時間曲線Fig.6 Absorbed energy-time curves ofthree kinds of components
為對主要吸能件進行快速優(yōu)化設(shè)計,本文對前縱梁、前指梁、副車架等部分進行簡化。參照文獻[4]的方法提取主要相關(guān)件的剛度特性(包括壓潰剛度特性和彎曲剛度特性)和慣性特性,并對剛度特性和慣性特性進行參數(shù)化,在Pam-Crash/Safe軟件中建立等效簡化模型。具體的構(gòu)件模擬及連接方式如下:主要吸能構(gòu)件中非板件采用梁單元模擬,材料類型選擇Beam214號材料,將參數(shù)化后的剛度特性曲線和慣性特性參數(shù)賦給相應(yīng)構(gòu)件的材料模型;板件則保留實際結(jié)構(gòu),與梁單元之間的連接方式采用Nodal Constraint節(jié)點連接。對吸能較小的構(gòu)件,采用Beam201號材料模擬,只輸入慣性特性參數(shù)。發(fā)生彎曲變形產(chǎn)生塑性鉸的部位采用Spring220彈簧單元模擬[5,6]。
以前縱梁為例說明等效簡化模型的建立過程。根據(jù)碰撞后變形特點,將前縱梁大致分成圖7所示的4段進行簡化:第1段主要為壓潰變形;第2段既有壓潰變形又有繞y軸和z軸的彎曲變形;第3段為繞y軸和z軸的彎曲變形;第4段為繞y軸和z軸的彎曲變形。等效簡化模型見圖8。采用同樣的方法對前指梁和副車架進行簡化。
圖7 根據(jù)前縱梁變形特點分段Fig.7 Segmentations of front characteristics of deformation
圖8 前縱梁簡化模型圖Fig.8 Simplified model side rail according toof front side rail
B柱之后的單元與前端侵入和乘員傷害的關(guān)聯(lián)性小,用質(zhì)量點代替。分別調(diào)整簡化模型車輛前端(A柱之前)、后端(B柱之后)及整車的質(zhì)量和質(zhì)心位置,使之與有限元基礎(chǔ)模型接近,得到最終的等效簡化模型。圖9為等效簡化后的模型,藍色的為被簡化的構(gòu)件。表1為有限元基礎(chǔ)模型與等效簡化模型質(zhì)量和質(zhì)心位置對比數(shù)據(jù)。
圖9 等效簡化模型Fig.9 Equivalent simplified model
圖10為等效簡化模型與有限元基礎(chǔ)模型的變形過程對比圖,兩種模型的變形形式比較一致。
表1 有限元基礎(chǔ)模型與等效簡化模型質(zhì)量和質(zhì)心位置對比Table 1 Comparison of masses and centroid locations between basic finite element model and equivalent simplified model
圖10 等效簡化模型與有限元基礎(chǔ)模型的變形過程對比Fig.10 Deformation Comparisons between basic finiteelement model and equivalent simplified model
等效簡化模型最終吸能量為158 kJ,誤差為4.6%。等效簡化模型7個相關(guān)監(jiān)測點的侵入量與有限元基礎(chǔ)模型相比最大誤差小于15%,如圖11所示。
圖11 等效簡化模型與有限元基礎(chǔ)模型監(jiān)測點侵入量對比Fig.11 Comparison of intrusion at monitoring sitesbetween equivalent simplified model andbasic finite element model
一般來說,正面25%重疊率碰撞工況不同于全寬碰撞,車輛脫離壁障時(本文將該時刻定義為碰撞結(jié)束時刻)仍有一定的速度v。車體前端總的目標吸能量按下式計算:
(1)
式中:E0為車體前端總的目標吸能量;m為等效簡化模型總質(zhì)量(見表1);v0為碰撞初速度;v為碰撞結(jié)束速度;η定義為車體前端目標吸能比例,即車體前端目標吸能量與整車碰撞前后動能差之比。
本文經(jīng)統(tǒng)計發(fā)現(xiàn),正面25%重疊率碰撞工況下80%以上的車輛在碰撞結(jié)束時刻的速度為20~30 km/h,本文取v為25 km/h。根據(jù)文獻[7],在正面全寬碰撞中為保證乘員艙不產(chǎn)生過大侵入,設(shè)計時一般令乘員艙吸能量占總吸能量的20%~30%。正面25%重疊率碰撞時,乘員艙相當于單側(cè)吸能,故本文將其吸能量限制在10%~20%,即車體前端目標吸能比例η為80%~90%,取η為0.85。由此得車體前端總目標吸能量為143 kJ。
原車輛-假人有限元基礎(chǔ)模型前端結(jié)構(gòu)除前指梁、前縱梁、副車架以外的構(gòu)件共吸能33 kJ,因此前指梁、前縱梁、副車架的總吸能量需盡可能達到143-33=110 kJ。
反復(fù)修改等效簡化模型中前指梁、前縱梁、副車架的特性參數(shù),使7個相關(guān)監(jiān)測點的侵入量不超過法規(guī)臨界侵入量(誤差控制在10%以內(nèi))。獲得此時前縱梁、前指梁及副車架的吸能量-時間曲線,如圖12所示。從圖12中得出前縱梁(含吸能盒)占這3個構(gòu)件總吸能量的目標吸能比例為55%,因此前縱梁(含吸能盒)的目標吸能量為60 kJ。
圖12 前端薄壁梁構(gòu)件吸能量-時間曲線Fig.12 Absorbed energy-time curves of frontalthin-walled beams
前縱梁是正面25%重疊率碰撞工況下最主要的吸能和傳力構(gòu)件,其前段和中段主要為薄壁直梁,通過壓潰變形吸收能量;后段與乘員艙地板縱梁相連,為曲梁,主要為彎曲變形。
考慮車身總布置情況,將前縱梁后段分成兩段來設(shè)計。將碰撞后等效簡化模型中前縱梁各段變形吸能量占整個前縱梁吸能量的比例作為前縱梁各段目標吸能比例,即第1段(含吸能盒)目標吸能比例為57%,第2段為36%,第3、4段為7%。
考慮該工況右側(cè)縱梁未充分變形只吸收很小的一部分能量,一般為左側(cè)前縱梁的30%,將左側(cè)前縱梁的目標吸能量(總目標吸能量減去右側(cè)前縱梁目標吸能量,此處為45 kJ)作為單側(cè)前縱梁的設(shè)計目標。為使設(shè)計的前端結(jié)構(gòu)也適應(yīng)正面全寬碰撞和40%偏置碰撞,左、右兩側(cè)前縱梁的設(shè)計保持一致。綜合以上,前縱梁第1段(含吸能盒)主要變形方式為壓潰變形,目標吸能量為26 kJ;第2段主要變形方式為壓潰變形,目標吸能量為16 kJ;第3、4段主要變形方式為彎曲變形,兩段總目標吸能量為3 kJ。采用相同方法算出前指梁和副車架各段目標吸能量,如表2所示。
前縱梁4段x向總長度為1188 mm。由發(fā)動機、車輪等布置要求確定前縱梁第1段(含吸能盒)、第2段、后兩段x向設(shè)計長度分別為403、330、455 mm。其中,吸能盒x向長度與有限元基礎(chǔ)模型相同,均為123 mm。各段吸能量、平均結(jié)構(gòu)力按下列公式計算:
(2)
Eqi=Fqi×Sqi
(3)
Sqi=kqi×Dqi
(4)
式中:Eq為前縱梁壓潰變形目標吸能量;Eqi為前縱梁第i段壓潰變形目標吸能量;Fqi為前縱梁第i段目標平均軸向結(jié)構(gòu)力;Sqi為前縱梁第i段壓潰長度;kqi為前縱梁第i段壓縮系數(shù);Dqi為前縱梁第i段原長度。
表2 前縱梁、前指梁和副車架各段變形方式及目標吸能量Table 2 Deformation modes and target energy absorption of each segment of front side rail,shotgun and sub-frame
前縱梁前兩段為壓潰變形,綜合考慮前兩段壓縮系數(shù)(一般為0.6~0.9[8])和逐級壓潰系數(shù)(為保障逐級壓潰一般為1.1~1.2[9])要求,由式(2)(3)(4)得到前縱梁第1段(含吸能盒)、第2段的目標平均軸向結(jié)構(gòu)力分別為92、101 kN。
為使吸能盒在碰撞過程中充分變形吸能,將其材料選為2.0T/DC01。本文采用低碳鋼矩形截面設(shè)計前端構(gòu)件。低碳鋼矩形截面薄壁梁平均軸向力計算公式如下[8,10]:
(5)
式中:F為平均軸向結(jié)構(gòu)力;σ0為平均流動應(yīng)力,低碳鋼取486MPa;b為矩形截面長寬均值;t為薄壁梁厚度;υ0為動態(tài)載荷;c為特征應(yīng)變率,低碳鋼取40s-1,p為材料敏感度度量值,低碳鋼取6。
取前縱梁1、2段厚度與有限元基礎(chǔ)模型相同,均為2.0 mm,由式(5)得到前縱梁1、2段斷面的長寬均值分別為161、117 mm,取第1段截面長162 mm、寬160 mm,第2段截面長140 mm、寬94 mm。前縱梁后兩段主要為彎曲變形,壓潰變形很小,且與地板縱梁相連,屬漸變面截面。設(shè)計這兩段斷面時綜合力的傳遞和彎曲變形吸能量要求,參照文獻[8]將厚度取為3.0 mm。最終前縱梁設(shè)計前后模型見圖13。
采用同樣的設(shè)計方法得到前指梁和副車架的設(shè)計方案,前指梁和副車架設(shè)計前后模型分別見圖14和圖15(設(shè)計后新模型均未進行工藝設(shè)計)。
圖13 改進設(shè)計前、后前縱梁模型Fig.13 Front side rail models before and afterimproved design
圖14 改進設(shè)計前、后前指梁模型Fig.14 Shotgun models before and after improveddesign
圖15 改進設(shè)計前、后副車架模型Fig.15 Sub-frame models before and after improved design
前縱梁、前指梁和副車架各段目標平均軸向結(jié)構(gòu)力、材料、厚度及截面尺寸等見表3。
為了驗證所設(shè)計車身前端結(jié)構(gòu)相關(guān)構(gòu)件的合理性,用設(shè)計的前縱梁、前指梁和副車架替換原有限元基礎(chǔ)模型中相應(yīng)的構(gòu)件得到新的有限元模型。
正面25%重疊率碰撞工況下,新模型7個監(jiān)測點侵入量除歇腳板和左側(cè)地板外其他監(jiān)測點實際侵入量均顯著減小,且未超過法規(guī)臨界侵入量(見圖16)。
新模型與原模型在正面全寬碰撞工況下的B柱加速度波形對比如圖17所示,加速度峰值略低于原模型,峰值點稍有提前,不會對約束系統(tǒng)匹配帶來新的問題[11]。
表3 前縱梁、前指梁和副車架各段目標平均軸向結(jié)構(gòu)力、材料、厚度及截面尺寸Table 3 Target average axial structural forces, materials, thickness and cross-sections sizes ofeach segment of front side rail, shotgun and sub-frame
圖16 改進結(jié)構(gòu)正面25%重疊率碰撞乘員艙侵入量Fig.16 Cabin intrusion of improved structure in25% overlap frontal collision
圖17 改進結(jié)構(gòu)與基礎(chǔ)模型正面全寬碰撞B柱加速度波形Fig.17 B-pillar acceleration waveforms of improvedstructure and basic model in frontal fullwidth collision
將新模型正面40%偏置碰撞乘員艙侵入量與法規(guī)臨界侵入量(該工況評價中達到滿分的侵入量)對比,結(jié)果表明,各監(jiān)測點侵入量均小于法規(guī)臨界值。
作者提出的等效簡化模型方法,一方面減少了優(yōu)化計算量;另一方面,該方法以乘員艙侵入量作為約束條件,反推主要吸能構(gòu)件的目標吸能量,再以此為目標,利用薄壁梁理論模型快速進行斷面設(shè)計,從一定程度上實現(xiàn)了“正向設(shè)計”。
[1] Insurance Institute for Highway Safety (IIHS). Small overlap frontal crashworthiness evaluation rating protocol (version I) [S].
[2] Munjurulimana D, Nagwanshi D, Marks M. Small overlap impact countermeasures for automobiles [C]∥SAE Paper, 2015-01-1491.
[3] Elliott E, Roche C, Reddy J. Small overlap impact countermeasure-front door hinge pillar dual box [C]∥SAE Paper, 2016-01-0402.
[4] 武櫟楠. 汽車正面抗撞性概念設(shè)計關(guān)鍵技術(shù)研究[D].長春:吉林大學(xué)汽車工程學(xué)院,2016.
Wu Li-nan. Research on key technology for conceptual design of vehicle frontal crashworthiness [D].Changchun: College Automotive Engineering, Jilin University,2016.
[5] Kecman D. Bending collapse of rectangular and square section tubes [J]. International Journal of Mechanical Sciences, 1983, 25(9/10): 623-636.
[6] Nikravesh P E, Chung I S, Benedict R L. Plastic hinge approach to vehicle crash simulation [J]. Computers & Structures, 1983, 16(1-4): 395-400.
[7] Malen D E. Fundamentals of Automobile Body Structure Design [M].Warrendale, Pennsylvania (USA):SAE International, 2011:224-226.
[8] 劉樂丹. 某乘用車正面碰撞波形優(yōu)化及其方法研究[D].長春:吉林大學(xué)汽車工程學(xué)院,2011.
Liu Le-dan. Research on optimization and methods for frontal crash waveforms of a passenger car [D].Changchun: College Automotive Engineering, Jilin University, 2011.
[9] 郝玉敏. A15EV電動車車體結(jié)構(gòu)正面抗撞性研究[D].長春:吉林大學(xué)汽車工程學(xué)院,2015.
Hao Yu-min. Research on frontal crashworthiness of A15EV electric vehicle body structure [D]. Changchun: College Automotive Engineering, Jilin University,2015.
[10] Abramowicz W. Thin-walled structures as impact energy absorbers [J]. Thin-Walled Structures, 2003, 41(2): 91-107.
[11] 王大志. 基于乘員保護的汽車正面碰撞結(jié)構(gòu)設(shè)計與變形控制研究[D].北京:清華大學(xué)機械工程學(xué)院,2006.
Wang Da-zhi. Research on structure design and deformation control for vehicle frontal collision based on occupants protection [D]. Beijing: College of Mechanical Engineering,Tsinghua University, 2006.