張文韜

摘要：根據(jù)汽車的正面碰撞特性，應(yīng)用有限元模擬技術(shù)，通過利用ANSA軟件建立某汽車前保險(xiǎn)杠與正面剛性墻的碰撞模型，并用顯式動(dòng)力學(xué)分析軟件LS-dyna求解器求解該模型，最后基于Hyperviews軟件進(jìn)行保險(xiǎn)杠強(qiáng)度的動(dòng)態(tài)分析，得到相應(yīng)部件應(yīng)變力的應(yīng)變?cè)茍D，對(duì)薄弱部分提出更改建議并驗(yàn)證其可行性，為整車碰撞分析提供參考。

關(guān)鍵詞：前保險(xiǎn)杠;碰撞;有限元分析

中圖分類號(hào)：U463.326? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ?文獻(xiàn)標(biāo)識(shí)碼：A? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 文章編號(hào)：1674-957X（2020）21-0055-02

0? 引言

汽車是伴隨著工業(yè)文明而出現(xiàn)的，是工業(yè)文明的標(biāo)志之一，繼蒸汽汽車之后，燃油汽車和電動(dòng)汽車先后誕生，汽車保有量日益增加，交通事故也不斷出現(xiàn)。由此可見，關(guān)于汽車的行駛安全正在被大家廣泛關(guān)注，因此對(duì)于如何將被動(dòng)安全性能提高，在標(biāo)準(zhǔn)法規(guī)內(nèi)，最大程度的減少駕駛?cè)藛T和行人的傷害成為汽車有限元分析中的重要課題。

據(jù)調(diào)查數(shù)據(jù)顯示正面碰撞發(fā)生的幾率是所有交通事故中最高的，并且前保險(xiǎn)作為在碰過程中承擔(dān)吸能和承載的主要部件，對(duì)于提高整車碰撞性能有非常重要的意義。

由于碰撞過程非常短，而表現(xiàn)的結(jié)構(gòu)非線性和材料非線性又非常復(fù)雜且始終伴隨動(dòng)態(tài)接觸[1]，因此本文采用有限元分析理論和顯式算法相結(jié)合來研究汽車前保險(xiǎn)杠的碰撞性能，在滿足法規(guī)的前提下對(duì)前保險(xiǎn)杠進(jìn)行優(yōu)化改進(jìn)，為以后的整車碰撞分析提供參考。

1? 接觸碰撞分析的基本過程

當(dāng)發(fā)生碰撞時(shí)，垂直于接觸界面的速度是瞬時(shí)不連續(xù)的，沿界面的切向速度也是不連續(xù)的[2]，這些特點(diǎn)會(huì)給基于離散方程進(jìn)行的時(shí)間積分帶來問題。針對(duì)這一問題，本文采用LS-dyna顯式動(dòng)力學(xué)分析軟件對(duì)保險(xiǎn)杠和剛性墻進(jìn)行受力狀況的求解。為降低計(jì)算機(jī)求解時(shí)間，導(dǎo)入的CATIA整車模型只保留由橫梁，吸能盒支架組成的保險(xiǎn)杠系統(tǒng)，然后將三維模型導(dǎo)入ANSA軟件中進(jìn)行離散化網(wǎng)格劃分，并對(duì)網(wǎng)格賦予材料、厚度、重量等屬性，不同部件之間采用焊點(diǎn)鏈接，建立剛性墻，生成K文件，遞交LS-dyna求解器進(jìn)行求解計(jì)算，最后將計(jì)算結(jié)果導(dǎo)入Hyperviews中進(jìn)行計(jì)算結(jié)果的后處理。

1.1 搭建保險(xiǎn)杠和剛性墻的有限元模型

將保險(xiǎn)杠系統(tǒng)相關(guān)部件的CATIA模型導(dǎo)入ANSA軟件中進(jìn)行網(wǎng)格劃分，網(wǎng)格尺寸選用10mm，由于保險(xiǎn)杠屬于鈑金件，所以使用BT殼單元?jiǎng)澐志W(wǎng)格;材料為低碳鋼板，對(duì)應(yīng)材料編號(hào)為MAT24，該材料的本構(gòu)關(guān)系如式（1）所示為：

（1）

式中：為屈服應(yīng)力;為有效塑性應(yīng)變;C、p為材料常數(shù)，通過設(shè)置C值和p值可以求取其動(dòng)態(tài)應(yīng)力應(yīng)變關(guān)系，針對(duì)本文的低碳鋼，可取C=40、p=5。

剛性墻由于定義為不可壓縮，不可變形，該結(jié)構(gòu)部件的相對(duì)位置保持不變。因此，使用MAT20材料模型。

1.2 沙漏控制

考慮到計(jì)算時(shí)間，為了提升計(jì)算的效率，絕大部分積分計(jì)算一般采用部分積分類型來進(jìn)行計(jì)算。因此，在采用部分積分過程中有可能出現(xiàn)沙漏，當(dāng)出現(xiàn)的沙漏問題較嚴(yán)重時(shí)則會(huì)導(dǎo)致模型的能量不守恒，影響了模型的計(jì)算精度甚至導(dǎo)致其計(jì)算結(jié)果的可信度降低，所以必須對(duì)計(jì)算沙漏進(jìn)行相關(guān)控制。其中，沙漏控制分為兩種一個(gè)是對(duì)單個(gè)件的局部沙漏控制，另外一個(gè)是對(duì)整體的總體沙漏控制兩種[3]。

1.3 時(shí)間步長(zhǎng)控制

時(shí)間步長(zhǎng)為每一步有限元積分的時(shí)間長(zhǎng)度[4]。顯式中心差分法的穩(wěn)定性決定于時(shí)間長(zhǎng)度，如果積分步長(zhǎng)大于則會(huì)發(fā)生計(jì)算不穩(wěn)定，導(dǎo)致計(jì)算精度迅速下降。其中為系統(tǒng)阻尼比，為系統(tǒng)最大固有頻率。

其中板殼單元的最大穩(wěn)定時(shí)間步長(zhǎng)為：

（2）

Ls為單元的特征長(zhǎng)度，ρ為質(zhì)量密度，E為彈性模量。

1.4 接觸定義

通過罰函數(shù)因子來計(jì)算得到接觸反力的大小，實(shí)現(xiàn)阻止相關(guān)節(jié)點(diǎn)對(duì)面的穿透。但存在著罰剛度過大可導(dǎo)致迭代次數(shù)增加，甚至引起剛度矩陣病態(tài)的問題，所以在選擇罰剛度時(shí)不能過大，應(yīng)該是逐漸增長(zhǎng)的，取結(jié)果穩(wěn)定后的最小值[5]。

由于保險(xiǎn)杠各部件之間的鏈接均為點(diǎn)焊連接，在ANSA軟件中通過關(guān)鍵字SPOTWELD模擬點(diǎn)焊連接;在此碰撞過程中，不考慮車身位移，因而為了節(jié)省計(jì)算時(shí)間，將車身簡(jiǎn)化為一個(gè)剛性面，對(duì)車身剛性面和保險(xiǎn)杠支架支撐處的六個(gè)自由度進(jìn)行約束，以模擬保險(xiǎn)杠和車體的鏈接。將保險(xiǎn)杠橫梁和吸能盒厚度均設(shè)置為1.8mm，行駛初速度設(shè)置為4m/s，碰撞時(shí)長(zhǎng)設(shè)置為100ms，完成模型搭建工作，提交LS-dyna進(jìn)行求解計(jì)算。最終搭建的保險(xiǎn)杠與剛性墻模型如圖1所示。

2? 保險(xiǎn)杠撞擊剛性墻的分析結(jié)果

汽車在猛烈的撞擊中，沖擊力非常大，響應(yīng)時(shí)間也非常短，保險(xiǎn)杠各個(gè)部件的位移量、速度、加速度都發(fā)生了較大的變化。本次分析計(jì)算最終求解時(shí)間為2.5h，碰撞時(shí)長(zhǎng)為100ms，利用Hyperviews后處理軟件得到圖2在不同時(shí)間段內(nèi)，保險(xiǎn)杠橫梁和吸能盒的變形狀況。

由應(yīng)力應(yīng)變?cè)茍D可知，在發(fā)生碰撞的瞬間，保險(xiǎn)杠橫梁產(chǎn)生了較大的塑性形變，同時(shí)，該材料的屈服應(yīng)力極限不足以支撐吸能盒變形的局部應(yīng)力，也消耗了碰撞所產(chǎn)生的巨大能量，降低了整車動(dòng)能。這種情況下，對(duì)汽車翼子板、車燈、前機(jī)罩、白車身以及乘員安全都有較大的影響。所以，為了充分發(fā)揮保險(xiǎn)杠系統(tǒng)的吸能效果，有必要對(duì)其部件的相關(guān)參數(shù)進(jìn)行調(diào)整，提高保險(xiǎn)杠系統(tǒng)的吸能效果和碰撞特性。

3? 優(yōu)化及驗(yàn)證

本文針對(duì)汽車前保險(xiǎn)杠的優(yōu)化分析，目的在于保護(hù)其后面的汽車零部件如發(fā)動(dòng)機(jī)、翼子板、機(jī)罩等不受損害，保護(hù)駕駛?cè)藛T和乘員不受傷害。因此，在碰撞發(fā)生的過程中，應(yīng)使吸能盒的壓潰吸能盡可能多的吸收碰撞時(shí)所產(chǎn)生的能量，提高被動(dòng)安全性能。

此次優(yōu)化以保險(xiǎn)杠橫梁和吸能盒的壁厚作為優(yōu)化參數(shù)，取橫梁中下端第4154個(gè)節(jié)點(diǎn)相對(duì)于吸能盒末端即車體剛性面的相對(duì)位移為參考指標(biāo)，分析在三種不同壁厚的情況下，保險(xiǎn)杠的吸能狀況，三種方案，：case01橫梁壁厚1.5mm，吸能盒壁厚1.5mm，相對(duì)位移最大值171.329mm;case02：1.5mm，2.5mm，105.115mm;case03：2.5mm，2.5mm，102.945mm。

仿真結(jié)果在HyperGraph2D中可以得到相應(yīng)節(jié)點(diǎn)位移、時(shí)間和吸能、時(shí)間的歷程曲線，如圖3、圖4所示。

由三種方案可知，case01中的最大相對(duì)位移為171mm，case02和case03中的相對(duì)最大位移相差不多，遠(yuǎn)小于case01;從圖4可以看出case02和case03的吸能效果相差不大，而case03中的橫梁壁厚要比case02厚1mm，從車體輕量化和吸能效果的綜合考慮來看，選用case02中的壁厚參數(shù)最為合理。

4? 結(jié)語

①本文利用前處理軟件ANSA和顯式動(dòng)力學(xué)分析軟件LS-DYNA聯(lián)合建模，通過Hyperviews后處理軟件進(jìn)行優(yōu)化分析，得出最優(yōu)化保險(xiǎn)杠壁厚參數(shù)，既提高了吸能盒的吸能效果，又保證了整車體的輕量化，縮短了設(shè)計(jì)周期，減小了實(shí)驗(yàn)次數(shù)，為以后整車碰撞分析提供了重要的參考價(jià)值;②首次采用ANSA前處理軟件進(jìn)行建模工作，在網(wǎng)格劃分，建立接觸等方面的優(yōu)越性要強(qiáng)于傳統(tǒng)建模軟件，但在分析后的結(jié)果處理方面有明顯不足;③在本次論文中提到的保險(xiǎn)杠系統(tǒng)的改進(jìn)中，只做了對(duì)原本數(shù)模參數(shù)的更改，并沒有對(duì)數(shù)模本身進(jìn)行修改，在實(shí)際的企業(yè)級(jí)項(xiàng)目中，需要與相關(guān)設(shè)計(jì)人員進(jìn)行協(xié)調(diào)，達(dá)到更理想的效果。

參考文獻(xiàn)：

[1]李瑞強(qiáng).某乘用車后保險(xiǎn)杠有限元仿真與實(shí)驗(yàn)研究[D].吉林大學(xué)，2017.

[2]白金澤.LS-DYNA3D理論基礎(chǔ)與實(shí)例分析[M].北京：科學(xué)出版社，2005：50-51.

[3]胡遠(yuǎn)志.基于LS-DYNA和Hyperworks的汽車安全仿真與分析[M].北京：清華大學(xué)出版社，2011：15-16.

[4]田雨蒙.全鋁車身電動(dòng)汽車正面碰撞耐撞性分析及結(jié)構(gòu)優(yōu)化[D].東北林業(yè)大學(xué)，2017.

[5]樊建平，吳孟暢.基于界面接觸算法的車橋耦合動(dòng)力學(xué)分析[J].中國(guó)公路學(xué)報(bào)，2017，4：54-55.

[6]呂峰，陳方根.汽車前艙結(jié)構(gòu)耐撞性關(guān)鍵技術(shù)研究[J].交通信息與安全，2015，5：112-114.

[7]李兆凱.汽車前部結(jié)構(gòu)低速碰撞耐撞性與輕量化優(yōu)化[J].中國(guó)公路學(xué)報(bào)，2016，10：149-157.

[8]姚宙，郝玉敏，李紅建.基于能量管理與仿真的汽車前端結(jié)構(gòu)優(yōu)化設(shè)計(jì)[J].汽車安全與節(jié)能學(xué)報(bào)，2016，07：395-402.