吳永強(qiáng)，安紅萍

(太原科技大學(xué)材料科學(xué)與工程學(xué)院，太原030024)

疏松缺陷是在鋼錠凝固過程中不可避免的現(xiàn)象，在熱成形過程中若不能有效控制，將會(huì)嚴(yán)重降低鍛件的使用壽命甚至使鍛件報(bào)廢。長期以來，國內(nèi)外學(xué)者對鍛造過程中存在疏松缺陷的的演變及控制做了大量的研究，總結(jié)出疏松缺陷閉合的規(guī)律和提出疏松缺陷閉合的變形臨界條件，為制定合理的鍛造工藝提供可靠依據(jù)[1]。有限元數(shù)值模擬方法是研究疏松缺陷閉合的常用方法，它即能夠獲得塑性加工制品內(nèi)部疏松缺陷的閉合過程，又能夠獲得影響其閉合過程相應(yīng)的變形條件(溫度、應(yīng)變速率等參數(shù))。事實(shí)上，大鍛件與其內(nèi)部存在的疏松缺陷尺寸相差懸殊，利用有限元分析軟件在相對較大的工件中很難作出尺寸很小的疏松缺陷，即疏松缺陷尺寸大小與整體工件尺寸大小比例條件達(dá)不到實(shí)際生產(chǎn)條件，數(shù)值模擬分析的精確性受到質(zhì)疑。很多研究者利用有限元數(shù)值模擬方法對鍛件內(nèi)部存在的孔洞閉合現(xiàn)象進(jìn)行了模擬計(jì)算，例如:清華大學(xué)王祖唐利用數(shù)值模擬的方法對大型鍛件孔洞閉合過程進(jìn)行分析，并利用鋼坯內(nèi)人工造孔及劃刻網(wǎng)格的方法進(jìn)行數(shù)據(jù)分析得出孔洞沿最大壓應(yīng)變方向閉合的結(jié)果[2-3];崔振山、任廣升等用剛塑性有限元的數(shù)值模擬和塑泥試樣的物理模擬分析了在鐓粗過程中工件內(nèi)部孔洞缺陷演化過程，得到孔洞周圍存在一個(gè)應(yīng)力應(yīng)變集中區(qū)域，并且應(yīng)變集中程度隨著孔洞尺寸的變小而越劇烈[4-5];西北工業(yè)大學(xué)呂毅利用有限元計(jì)算細(xì)觀力學(xué)法(有限元法、細(xì)觀力學(xué)和材料科學(xué)相結(jié)合的方法)做了《基于有限元計(jì)算細(xì)觀力學(xué)的復(fù)合材料宏觀性能一體化預(yù)測》，該課題的方法關(guān)鍵是以先進(jìn)的前后處理軟件MSC.PATRAN-R2為平臺，在該平臺上用 PCL(patran command language)為開發(fā)工具建立合適的代表性特征體積元(RVE)，實(shí)現(xiàn)了代表性特征體積元(RVE)參數(shù)化自動(dòng)建模，建立了代表性特征體積元(RVE)庫[6]。

本文針對代表性特征體積元(RVE)施加適當(dāng)?shù)倪吔鐥l件進(jìn)行有限元數(shù)值模擬分析，利用代表性特征體積元(RVE)應(yīng)力應(yīng)變場對宏觀鍛件內(nèi)部疏松缺陷的閉合行為進(jìn)行預(yù)測，為以后根據(jù)溫度、應(yīng)變速率、疏松缺陷大小合理制定在鍛造過程中相應(yīng)的工藝參數(shù)做準(zhǔn)備，以及為制造質(zhì)量優(yōu)良的塑性加工制品提供幫助。

1 宏、細(xì)觀模型的建立

疏松缺陷的變形過程是由疏松周圍基體材料的運(yùn)動(dòng)引起的，在塑性加工過程中由于外力加載，基體材料產(chǎn)生流動(dòng)并填充到疏松缺陷中，疏松缺陷閉合。為分析疏松閉合過程建立宏觀與細(xì)觀兩種尺度下模型，具體做法是首先通過有限元分析軟件Deform-2D建立宏觀模擬模型，而細(xì)觀模型中代表性特征體積元(RVE)的邊界條件從宏觀模型中獲得，宏細(xì)觀兩種尺度下的模型關(guān)系如圖1所示。

圖1 宏細(xì)觀模型關(guān)系Fig.1 The relationship between the macro model and the mesoscopic one

圖2表示細(xì)觀模型中代表性特征體積單元(RVE)由基體和疏松缺陷組成。這個(gè)代表性特征體積單元(RVE)處于無限大的基體之中，周圍材料對它產(chǎn)生約束作用。提取宏觀模型中特定點(diǎn)處的加載歷史作為細(xì)觀模型中代表性體積元(RVE)的邊界條件，研究與宏觀模型一致的變形過程中疏松閉合情況。通過一些簡單的本構(gòu)關(guān)系描述代表性特征體積單元(RVE)在變形過程中的行為(如:力學(xué)行為等)，再用無數(shù)具有簡單本構(gòu)關(guān)系的小單元疊加就可以得到整個(gè)塑性變形工件的整體行為。

圖2 含疏松的代表性特征體積單元Fig.2 The representative volume element with loose defects

2 邊界施加方法

本課題研究的是利用細(xì)觀模型再現(xiàn)宏觀模型的變化過程，它的具體要求是在代表性特征體積元(RVE)上加載邊界條件，使它在與宏觀模型相同的時(shí)間內(nèi)有相同的變形規(guī)律，即具有相同的等效應(yīng)變、等效應(yīng)變速率、平均應(yīng)力等表征塑性加工制件在變形過程中的參數(shù)?；诤暧^基體材料的細(xì)觀結(jié)構(gòu)周期性假設(shè)，選取適當(dāng)?shù)拇硇蕴卣黧w積元(RVE)，在保證宏、細(xì)觀模型之間能量等效前提下，采取弱化的連續(xù)性邊界條件以及周期性邊界條件，建立和求解細(xì)觀力學(xué)方程。借助均勻化理論將大量代表性特征體積元(RVE)的表征參數(shù)平均化，建立宏觀與細(xì)觀兩種尺度坐標(biāo)范圍內(nèi)的一系列方程，得到宏觀和細(xì)觀兩種某表征參數(shù)下的細(xì)觀場表達(dá)式。在有限元分析軟件Deform-2D中，代表性特征體積元(RVE)的邊界條件從宏觀模型中獲得，為保證邊界條件數(shù)據(jù)傳遞的正確性，本文做了大量的邊界條件加載方案，最后得到了可行性的方案[7-9]。下面列舉本文做過的三種邊界條件加載方案:

(1)在Deform-2D軟件中建立宏觀模型，選擇特征單元，提取表達(dá)其變形歷史的應(yīng)力-時(shí)間及速度—時(shí)間數(shù)據(jù)作為細(xì)觀模型的邊界條件，建立如圖2所示的含有疏松缺陷的細(xì)觀模型(為簡化取1/4分析)，如圖3所示。

圖3 方案一:代表性特征體積元的四分之一Fig.3 Scheme one:One quarter of the representative volume element

在代表性特征體積元(RVE)的1/4的上邊和右邊同時(shí)加載應(yīng)力邊界條件或者同時(shí)加載速度邊界條件;或者在代表性特征體積元(RVE)的1/4的上邊加載應(yīng)力邊界條件右邊加載速度邊界條件;或者在其上邊加載速度邊界條件右邊加載應(yīng)力邊界條件。而在以上所有加載方法下代表性特征體積元(RVE)的1/4模型的下邊和左邊總是加載應(yīng)力邊界條件使細(xì)觀模型固定。實(shí)驗(yàn)結(jié)果是代表性特征體積元(RVE)變形的有關(guān)參數(shù)(等效應(yīng)變、等效應(yīng)力、平均應(yīng)力等)是宏觀模型的對應(yīng)參數(shù)的幾倍甚至是十幾倍，即從宏觀模型提取的邊界條件加載到細(xì)觀模型的數(shù)據(jù)傳遞不準(zhǔn)確，細(xì)觀模擬過程不能完全再現(xiàn)宏觀模擬過程。宏觀工件中右上角一點(diǎn)與在用該方案加載邊界條件的代表性特征體積元(RVE)中等效應(yīng)變在同一時(shí)刻的相對誤差值如表1所示:

表1 相對誤差百分比Tab.1 The relative error percentage

(2)與方案一相同，選擇宏觀模型中的單元提出其變形歷史構(gòu)建細(xì)觀模型。不同之處為在代表性特征體積元(RVE)的四邊同時(shí)加載應(yīng)力邊界條件或者同時(shí)加載速度邊界條件;或者采用上邊和右邊加載應(yīng)力邊界條件，左邊和下邊加載速度邊界條件的混合加載方式。有限元模擬結(jié)果是代表性特征體積元(RVE)變形的有關(guān)參數(shù)(等效應(yīng)變、等效應(yīng)力、平均應(yīng)力等)逐漸接近宏觀模型相對應(yīng)參數(shù)，即從宏觀模型提取的邊界條件到代表性特征體積元(RVE)的數(shù)據(jù)傳遞準(zhǔn)確性提高如表2所示，但細(xì)觀模擬過程還不能完全準(zhǔn)確再現(xiàn)宏觀模擬過程如圖4所示。

表2 相對誤差百分比Tab.2 The relative error percentage

圖4 方案二:對整個(gè)單元進(jìn)行加載Fig.4 Scheme two:the loading of the whole element

(3)建立宏觀模型，選擇特征節(jié)點(diǎn)，提取表達(dá)其變形歷史的應(yīng)力-時(shí)間及速度—時(shí)間數(shù)據(jù)作為細(xì)觀模型的邊界條件，由于宏觀模型比代表性特征體積元(RVE)的尺寸大很多，宏觀模型中這個(gè)點(diǎn)可視為一個(gè)微小的含缺陷區(qū)域。在建立細(xì)觀模型時(shí)可不再拘泥于宏觀模型中單元的尺寸大小如圖5所示:

為了方便有限元模擬分析把細(xì)觀模型尺寸定位1mm×1mm.在進(jìn)行邊界條件加載時(shí)如圖6所示黑色箭頭表示該代表性特征體積元上(RVE)加載的是速度邊界條件，紅色橫線表示該代表性特征體積元(RVE)加載的邊界條件是應(yīng)力邊界條件，中間位置紫色短道橫線表示疏松缺陷鍛合完畢后狀態(tài)。

圖5 方案三:宏觀模型中取其中一點(diǎn)Fig.5 Scheme three:one point in the macroscopic model

圖6 加載方案Fig.6 The loading project

在有限元分析軟件Deform-2D中，在加載邊界條件時(shí)，邊界條件的加載間隔步數(shù)過大容易使剛開始幾步的表征變形參數(shù)(等效應(yīng)變、等效應(yīng)力、平均應(yīng)力等)發(fā)生突變，與宏觀模型相應(yīng)參數(shù)相差很大。故不采用等間隔步數(shù)加載邊界條件，而使前幾步加密處理，讓邊界條件加載的間隔步數(shù)減少。宏細(xì)觀邊界條件數(shù)據(jù)傳遞準(zhǔn)確性的證明如圖7所示，其中縱坐標(biāo)誤差百分比即相對誤差百分比。

圖7 傳遞中的誤差Fig.7 The error in transfer process

綜上，采用方案三即采用左邊和下邊加載速度邊界條件，上邊和右邊加載應(yīng)力邊界條件，該方案有限元模擬出來的代表性特征體積元(RVE)的表征變形參數(shù)(等效應(yīng)變、等效應(yīng)力、平均應(yīng)力等)與宏觀模型的相應(yīng)表征變形參數(shù)非常接近，最大的誤差為10%，當(dāng)變形穩(wěn)定后，大多數(shù)誤差范圍所在數(shù)量級為10-2，有的宏細(xì)觀參數(shù)間誤差甚至達(dá)到千分之幾，即該細(xì)觀模型的模擬過程可以再現(xiàn)宏觀模型的變形歷史，而且能夠獲得缺陷閉合過程。在建立正確的有限元數(shù)值模型后，研究者就可以根據(jù)此方案在不同的溫度、應(yīng)變速率、疏松缺陷尺寸大小條件下進(jìn)行模擬分析，為消除鍛件內(nèi)部疏松缺陷和制定合理鍛壓工藝提供依據(jù)[10]。疏松缺陷閉合時(shí)形狀變化過程如圖8所示。

圖8 疏松缺陷的閉合過程Fig.8 The closing process of the loose defects

4 結(jié)論

本文基于有限元分析軟件Deform-2D在宏細(xì)觀兩種尺度下建立模型，細(xì)觀模型為含有疏松缺陷的代表性特征體積元(RVE)，通過在其上施加適當(dāng)?shù)倪吔缂s束條件和載荷，使從宏觀模擬提取的邊界條件數(shù)據(jù)恰當(dāng)?shù)厥┘拥酱硇蕴卣黧w積單元(RVE)上，從而保證從宏觀模型到細(xì)觀模型數(shù)據(jù)傳遞的準(zhǔn)確性即變形表征參數(shù)(等效應(yīng)變、等效應(yīng)變速率、等效應(yīng)力、平均應(yīng)力等)誤差極小，這樣就說明在鍛造過程中細(xì)觀模型的變形過程可以再現(xiàn)宏觀模型中的某一單元變形過程。研究者在消除大型鍛件內(nèi)部存在疏松缺陷時(shí)，可以對工件中存在嚴(yán)重缺陷的某一小范圍進(jìn)行模擬分析以便制定鍛造工藝，這樣既節(jié)省了設(shè)計(jì)時(shí)間也減少了在對整個(gè)工件全部分析過程中帶來的不準(zhǔn)確性。

[1]郭會(huì)光，曲宗實(shí).我國大鍛件制造業(yè)的發(fā)展[J].大型鑄鍛件，2003(1):1-2.

[2]王祖唐，劉莊，任猛.大型鋼錠內(nèi)部空洞性缺陷鍛合過程的實(shí)驗(yàn)研究[J].機(jī)械工程學(xué)報(bào)，1989，25(4):47-53.

[3]王祖唐，劉莊，任猛.大鋼錠內(nèi)部空洞鍛合過程的數(shù)值模擬[J].機(jī)械工程學(xué)報(bào)，1989，25(3):51-55.

[4]任廣升，譚紅，李運(yùn)興.鐓粗過程中鍛件內(nèi)孔洞缺陷鄰域應(yīng)變分布的模擬研究[J].塑性工程學(xué)報(bào)，1994，1(3):14-19.

[5]崔振山，任廣升，徐秉業(yè)，等.圓柱體內(nèi)部空洞的熱鍛閉合條件[J].清華大學(xué)學(xué)報(bào)，2003，43(2)，227-229，233.

[6]呂毅，呂國志，孫龍生.基于有限元計(jì)算細(xì)觀力學(xué)的RVE庫的建立與應(yīng)用[J].機(jī)械強(qiáng)度，2009.31(5):851-856.

[7]LUO CHUNHUI.Evolution of voids close to an inclusion in hot deformation of metals[J].Computational Materials Science，2001，21:360-374.

[8]LUO CHUNHUI，ULF STAHLBERG.Deformation of inclusion during hot rolling of steels[J].Journal of Materials Processing Technology，2001，114:87-97.

[9]LUO CHUNHUI.Modeling the behavior of inclusions in plastic deformation of steels[D].The doctor’s degree in KTH，2001.

[10]任運(yùn)來.大型鍛件內(nèi)部缺陷修復(fù)條件和修復(fù)方法研究[D].河北秦皇島:燕山大學(xué)，2003.