岑升波

【摘? 要】ABAQUS是一款功能非常強(qiáng)大的有限元軟件，被廣泛地應(yīng)用在解決各種工程項(xiàng)目問(wèn)題當(dāng)中。論文以防撞梁的熱沖壓成形為例，詳細(xì)介紹了ABAQUS在處理復(fù)雜接觸問(wèn)題時(shí)的基本步驟和注意事項(xiàng)，展示了ABAQUS在處理高度非線(xiàn)性問(wèn)題時(shí)的優(yōu)越性，也為新手在處理相似問(wèn)題時(shí)提供參考和幫助。

【Abstract】ABAQUS is a finite element software with very powerful functions， which is widely used to solve various problems in engineering projects. This paper takes the hot stamping of the anti-collision beam as an example， explicitly introduces the basic steps and precautions of ABAQUS in dealing with complex contact problems， shows the advantages of ABAQUS in dealing with highly nonlinear problems， and also provides reference and help for novices in dealing with similar problems.

【關(guān)鍵詞】ABAQUS軟件;熱沖壓成形;接觸

1 引言

ABAQUS是一套功能非常強(qiáng)大的工程模擬仿真軟件，擁有各種類(lèi)型的材料模型庫(kù)，可以模擬各種工程領(lǐng)域的許多問(wèn)題，例如，熱傳導(dǎo)、熱電耦合分析、聲學(xué)分析、巖土力學(xué)分析以及壓電介質(zhì)分析[1，2]。ABAQUS擁有十分強(qiáng)大的處理高度復(fù)雜非線(xiàn)性問(wèn)題的能力，具有模擬復(fù)雜系統(tǒng)問(wèn)題的高度可靠性，被廣泛應(yīng)用于在各國(guó)工業(yè)和研究所中。但是在國(guó)內(nèi)ABAQUS的普及率遠(yuǎn)不如ANSYS，不僅中文的資料少，而且可供新手學(xué)習(xí)的例子少，尤其是對(duì)比較復(fù)雜的接觸問(wèn)題的介紹更少。

在現(xiàn)實(shí)世界中，有許許多多的工程項(xiàng)目問(wèn)題都涉及接觸問(wèn)題，如運(yùn)動(dòng)中火車(chē)輪與鋼軌的接觸、切削中車(chē)刀與工件的接觸、沖壓時(shí)模具與毛坯間的接觸，等等。在接觸過(guò)程中，兩個(gè)物體在接觸界面上的相互作用是復(fù)雜的力學(xué)現(xiàn)象，接觸問(wèn)題常常同時(shí)涉及三種非線(xiàn)性：大形變引起的材料非線(xiàn)性、幾何非線(xiàn)性和接觸界面的非線(xiàn)性，再加上接觸界面的事先未知性和接觸條件的不等式約束，決定了接觸分析過(guò)程屬于高度非線(xiàn)性的，不僅需要花費(fèi)較多的計(jì)算機(jī)資源，而且收斂難度較大。本文從汽車(chē)防撞梁的熱沖壓成形出發(fā)，基于A(yíng)BAQUS軟件對(duì)接觸問(wèn)題進(jìn)行分析研究。

2 ABAQUS接觸分析的關(guān)鍵問(wèn)題

接觸分析中需要注意的常見(jiàn)問(wèn)題包括：?jiǎn)卧?lèi)型的選擇和網(wǎng)格質(zhì)量的控制、接觸關(guān)系的建立、相對(duì)滑動(dòng)、主從面的定義等[3]。ABAQUS軟件提供了非常豐富的單元種類(lèi)，在給用戶(hù)提供更多選擇的同時(shí)，也增加了用戶(hù)選擇適合自己模型的單元類(lèi)型的難度，人們要清楚自己計(jì)算的目的是什么，然后選擇合適的單元類(lèi)型。例如，熱沖壓成形模型，涉及溫度變化和形變，就這需要選擇溫度—位移耦合類(lèi)型。網(wǎng)格質(zhì)量的控制，ABAQUS軟件是通過(guò)Verify Mesh檢查模型網(wǎng)格質(zhì)量，可以檢查出分析過(guò)程中會(huì)導(dǎo)致錯(cuò)誤和警告信息的單元，想要獲得精確的計(jì)算結(jié)果，完美的單元質(zhì)量是必不可少的。接觸關(guān)系的建立，人們?cè)诮⒔佑|對(duì)時(shí)，首先要搞清楚主從面的定義，一般選擇剛度較大的一面作為主面，柔性材料作為從面，而解析面和剛體面必須作為主面。定義接觸屬性時(shí)，一定要清楚接觸里面包含了哪些內(nèi)容，如熱沖壓成形模型中模具和毛坯之間的互相接觸，既有摩擦又會(huì)發(fā)生熱傳遞，在定義屬性時(shí)不能遺漏。相對(duì)滑動(dòng)問(wèn)題，ABAQUS包括有限滑動(dòng)和小滑動(dòng)兩種選擇，其中，有限滑動(dòng)是最通用的選擇，其在接觸面間允許任意大的滑動(dòng)和旋轉(zhuǎn)，而小滑動(dòng)是在接觸面間允許小的相對(duì)滑動(dòng)。在接觸分析中需要注意的問(wèn)題還有很多，稍微不注意都可能導(dǎo)致計(jì)算不收斂，會(huì)在下文防撞梁的熱沖壓成形例子中提及。

3 熱沖壓成形案例

本案例基于某轎車(chē)前防撞梁熱沖壓成形的過(guò)程，詳細(xì)介紹了利用ABAQUS軟件解決接觸類(lèi)問(wèn)題的基礎(chǔ)步驟和注意事項(xiàng)。

3.1 案例的背景

汽車(chē)防撞梁是車(chē)輛受到碰撞時(shí)吸收碰撞能量的一種裝置，能夠減少駕乘人員受到的沖擊，是提高汽車(chē)安全可靠性的重要部件。為了實(shí)現(xiàn)車(chē)輛的輕量化且安全性能不降低，防撞梁一般都使用高強(qiáng)鋼22MnB5，采用熱沖壓成形的工藝沖壓成形。防撞梁的成形工藝為：將加熱至850℃左右完全奧氏體化的22MnB5毛坯置于25℃的凹模之上，然后壓邊圈以15000N的力向下運(yùn)動(dòng)直至最低點(diǎn)，毛坯發(fā)生變形且四周被壓邊圈壓緊，接著凸模以40m/s的速度向下運(yùn)動(dòng)，毛坯在凸模的作用下發(fā)生塑性變形，最后凸模、凹模完全閉合，進(jìn)行保壓和淬火。淬火階段是在凸模和凹模的冷卻管道中通入冷卻水，工件上大量的熱量通過(guò)熱傳遞被循環(huán)冷卻水帶走實(shí)現(xiàn)淬火，本文不做討論。在本案例中，毛坯和模具之間的接觸界面既存在力的相互作用又存在熱傳遞作用，在計(jì)算過(guò)程中，既有材料大變形帶來(lái)的非線(xiàn)性問(wèn)題又有溫度帶來(lái)的材料屬性的非線(xiàn)性問(wèn)題，這是屬于比較復(fù)雜的熱力耦合問(wèn)題。

3.2 模型的建立

防撞梁的三維模型如圖1所示，長(zhǎng)度方向的尺寸為1068.5mm，寬度方向的尺寸為135mm，弧度為22°，毛坯厚度為1mm。本文采用SolidWorks軟件建立防撞梁熱沖壓成形的三維模型，分別為凸模、凹模、壓邊圈和毛坯，如圖2所示。由于毛坯的長(zhǎng)寬都遠(yuǎn)遠(yuǎn)大于厚度，為了節(jié)約計(jì)算時(shí)間，將模型轉(zhuǎn)變?yōu)闅んw導(dǎo)入ABAQUS中。在導(dǎo)入部件的過(guò)程中應(yīng)注意，如果僅僅關(guān)注的是防撞梁成形的過(guò)程，不考慮溫度場(chǎng)的話(huà)，可以將凸模、凹模和壓邊圈都設(shè)置為離散剛體，可大大提高計(jì)算速度;但是如果還要關(guān)注防撞梁成形過(guò)程中毛坯和模具之間熱傳遞的溫度場(chǎng)的話(huà)，所有部件都要設(shè)為可變體，這對(duì)后面網(wǎng)格單元類(lèi)型的選擇很重要。

3.3 材料屬性的定義

毛坯22MnB5鋼和模具的密度為7850kg/m3，線(xiàn)膨脹系數(shù)為1.3×10-5。由于防撞梁的熱沖壓過(guò)程是采用熱力耦合分析的方法，溫度的變化對(duì)材料的屬性影響非常大，根據(jù)查詢(xún)文獻(xiàn)資料[4]，22MnB5鋼在不同溫度下的真實(shí)塑性應(yīng)力應(yīng)變數(shù)據(jù)如表1所示。

22MnB5鋼的楊氏模量、泊松比以及熱傳導(dǎo)率、熱熔的數(shù)據(jù)[5，6]，如表2所示。模具的熱導(dǎo)率為52W·（m·K）-1，熱熔為450J·（kg·K）-1。在設(shè)定材料屬性時(shí)，最容易出問(wèn)題的是各參數(shù)單位的統(tǒng)一，很容易出現(xiàn)數(shù)量級(jí)的錯(cuò)誤。

3.4 分析步的設(shè)定

分析步設(shè)為兩步，分析步的類(lèi)型選擇動(dòng)力，溫度—位移，顯式：第一步為壓邊圈向下運(yùn)動(dòng)，時(shí)間為0.001s;第二步為凸模向下運(yùn)動(dòng)，時(shí)間為0.0175s。在編輯分析步時(shí)，把幾何非線(xiàn)性打開(kāi)。如果想節(jié)約計(jì)算時(shí)間，可以將質(zhì)量縮放打開(kāi)，縮放比例的大小以計(jì)算結(jié)果的動(dòng)能占總能量的5%以?xún)?nèi)合適。

3.5 相互作用的設(shè)定

進(jìn)入相互作用模塊，對(duì)模具和毛坯之間設(shè)定相互作用關(guān)系。首先，相互作用屬性管理器定義接觸關(guān)系的屬性，其切向行為中的摩擦公式選擇罰公式，摩擦系數(shù)為0.05，其熱傳導(dǎo)系數(shù)當(dāng)距離為0時(shí)系數(shù)為0.75;其次，定義毛坯和模具與周?chē)h(huán)境對(duì)流的熱交換系數(shù)，毛坯與空氣的自然對(duì)流系數(shù)為29W/m2/K，模具與冷卻水的熱交換系數(shù)為1200W/m2/K;再次，定義毛坯和模具之間的接觸面，在這里要強(qiáng)調(diào)的是模具面一定是主面，毛坯面為從面;最后，創(chuàng)建三個(gè)參考點(diǎn)，創(chuàng)建約束類(lèi)型，選擇剛體，分別將凸模、凹模和壓邊圈綁定起來(lái)。

3.6 施加載荷

進(jìn)入載荷模塊，對(duì)模具和毛坯施加邊界條件，在凹模的參考點(diǎn)上施加位移/轉(zhuǎn)角邊界條件，將所有方向固定。在壓邊圈的參考點(diǎn)上也施加位移/轉(zhuǎn)角邊界條件，分析步選擇第一步，設(shè)定壓邊圈向下運(yùn)動(dòng)24mm壓住毛坯表面，其他方向固定。同時(shí)，設(shè)定幅值，幅值的類(lèi)型注意選擇平滑分析步，否則容易出現(xiàn)不收斂問(wèn)題。在凸模的參考點(diǎn)上施加位移/轉(zhuǎn)角邊界條件，分析步選擇第二步，設(shè)定凸模向下運(yùn)動(dòng)70mm進(jìn)行沖壓動(dòng)作，其他方向固定，幅值的設(shè)定時(shí)間與分析步一致。進(jìn)入預(yù)定義場(chǎng)管理器，毛坯的溫度設(shè)定為850℃，模具溫度為25℃。

3.7 網(wǎng)格劃分

進(jìn)入網(wǎng)格劃分模塊，先定義網(wǎng)格控制屬性，網(wǎng)格形狀盡量選擇四面體，劃分技術(shù)選擇結(jié)構(gòu)，有時(shí)需要拆分模型才能劃出想要的網(wǎng)格。進(jìn)入種子部件，給模型設(shè)置尺寸，然后劃分網(wǎng)格。一定要為網(wǎng)格指派單元類(lèi)型，本次計(jì)算為熱力耦合的問(wèn)題，且接觸里面定義了熱傳導(dǎo)屬性，所有模型的網(wǎng)格單元類(lèi)型一定要溫度—位移耦合，不然提交計(jì)算會(huì)報(bào)錯(cuò)“without temperature degree of freedom are not allowed in a heat transfer analysis”。

3.8 作業(yè)提交與結(jié)果分析

在工作管理器中提交作業(yè)，在監(jiān)控中可以實(shí)時(shí)觀(guān)察計(jì)算的進(jìn)展情況。當(dāng)計(jì)算結(jié)束以后進(jìn)入可視化模塊對(duì)結(jié)果進(jìn)行分析。當(dāng)熱沖壓結(jié)束后，防撞梁成形后凹槽底部溫度最低點(diǎn)的溫度—時(shí)間曲線(xiàn)圖如圖3所示，由于整個(gè)沖壓過(guò)程時(shí)間只有0.0185s，時(shí)間非常短暫，該節(jié)點(diǎn)的溫度從850℃下降到845.4℃，瞬間的冷卻速度達(dá)到248.6℃/s，遠(yuǎn)遠(yuǎn)大于馬氏體的轉(zhuǎn)變的冷卻速度，有利于提高22MnB5的馬氏體轉(zhuǎn)變率，提高防撞梁熱沖壓成形后的力學(xué)性能。

防撞梁熱沖壓成形后的應(yīng)力云圖如圖4所示，最大應(yīng)力出現(xiàn)在防撞梁凹槽底部圓角過(guò)渡兩側(cè)以及側(cè)壁上，最大應(yīng)力為169.9MPa，這是因?yàn)檫@些位置是毛坯與模具相互接觸、相互作用最為復(fù)雜的位置，也是材料形變最大的地方，所以其應(yīng)力也相對(duì)較大。

4 結(jié)語(yǔ)

ABAQUS軟件在處理高度非線(xiàn)性的復(fù)雜問(wèn)題上具有明顯的優(yōu)勢(shì)，本文以汽車(chē)防撞梁熱沖壓成形為例，詳細(xì)介紹了ABAQUS在處理復(fù)雜接觸問(wèn)題的過(guò)程和注意事項(xiàng)，為進(jìn)一步利用ABAQUS軟件研究汽車(chē)防撞梁的熱沖壓成形過(guò)程做好鋪墊，同時(shí)，也為新手處理類(lèi)似的問(wèn)題提供了參考和幫助。

【參考文獻(xiàn)】

【1】石亦平，周玉蓉.ABAQUS有限元分析實(shí)例詳解[M].北京：機(jī)械工業(yè)出版社，2006.

【2】莊茁，由小川，廖劍暉，等.基于A(yíng)BAQUS的有限元分析和應(yīng)用[M].北京：清華大學(xué)出版社，2009.

【3】劉展，朱精平，錢(qián)英莉，等.ABAQUS 6.6基礎(chǔ)教程與實(shí)例詳解[M].北京：中國(guó)水利水電出版社，2008.

【4】時(shí)立佳.B柱加強(qiáng)板熱沖壓成形工藝研究[D].秦皇島：燕山大學(xué)，2013.

【5】康永林，孫建林.高強(qiáng)度薄鋼研究的新進(jìn)展及其在汽車(chē)上的應(yīng)用[J].鋼鐵，2002，37（5）：65-70.

【6】谷諍巍，孟佳，李欣，等.超高強(qiáng)鋼熱成形奧氏體化加熱參數(shù)的優(yōu)化[J].吉林大學(xué)學(xué)報(bào)，2011（41）：194-197.