卜春陽(yáng)，袁戰(zhàn)偉，李 晶，何 凱

(1.金堆城鉬業(yè)股份有限公司技術(shù)中心，陜西 西安 710077)(2.長(zhǎng)安大學(xué)材料科學(xué)與工程學(xué)院，陜西 西安 710061)

0 引 言

金屬材料的熱加工是材料加工的重要方式，在成型過程中，變形溫度、加載速率、變形程度等均有顯著的影響[1，2]，合理的變形條件最終決定了材料成型的質(zhì)量，材料熱加工工藝是材料熱成型的基礎(chǔ)，一直以來(lái)也是材料研究熱點(diǎn)[3-6]。此外，在變形過程中，摩擦作用對(duì)材料的變形也有著重要的影響，一方面影響變形過程中加載載荷，而另一方面也對(duì)成型質(zhì)量有影響[7-9]。

對(duì)于摩擦作用對(duì)材料變形的影響，已有不少學(xué)者做了研究。李達(dá)人等[7]通過數(shù)值模擬方法確定了W-40%Cu粉末燒結(jié)材料在熱加工數(shù)值模擬過程中的摩擦因子。鄧華紅等[8]通過數(shù)值模擬研究了葉片精鍛過程中摩擦的作用，發(fā)現(xiàn)摩擦對(duì)溫度場(chǎng)和載荷形成曲線均有較大影響。馬勇等[9]采用有限元軟件分析了不同摩擦條件對(duì)7075鋁合金等通道角擠壓過程的影響，發(fā)現(xiàn)隨著摩擦因數(shù)增大，載荷峰值明顯增大甚至成倍增長(zhǎng)，且載荷值波動(dòng)加劇，等效應(yīng)力應(yīng)變分布不均勻。本文結(jié)合Deform 3D有限元軟件，研究了摩擦系數(shù)對(duì)2024鋁合金的熱模擬壓縮過程的影響。

1 有限元模型

本研究根據(jù)實(shí)際熱模擬壓縮過程采用Deform 3D有限元軟件建立模型，在模型中分別建立了上砧板、下砧板和試樣幾何模型，其中試樣大小為φ10×15 mm，材料為2 024鋁合金，試樣上的網(wǎng)格采用四面體單元，單元數(shù)量為11 200，畫網(wǎng)格后的模型如圖1所示。在加載過程中，上砧向下加載，下砧固定不動(dòng)，加載位移為9 mm。

圖1 有限元分析模型

在實(shí)際熱模擬壓縮試驗(yàn)過程中，采用石墨片作為潤(rùn)滑劑置于兩個(gè)砧子和試樣之間以減少摩擦作用，然而由于摩擦作用的存在，坯料變形后不可避免出現(xiàn)鼓肚，因此有必要對(duì)摩擦的作用進(jìn)行進(jìn)一步研究。在本研究中，為了研究摩擦作用對(duì)熱變形的影響，將摩擦系數(shù)f分別設(shè)為0.2、0.3、0.4和0.5，其他設(shè)置保持一致。

材料中的損傷用Deform 3D軟件中所帶的Cockcroft-Latham準(zhǔn)則損傷參數(shù)進(jìn)行表征，其表達(dá)式為[10]：

(1)

(2)

式中σ1，σ2，σ3三主應(yīng)力。

2 結(jié)果與分析

2.1 模擬結(jié)果云圖

整個(gè)模擬過程91步完成，下面以摩擦系數(shù)為0.2的模型結(jié)果為例，分別對(duì)第30、60和91步的結(jié)果進(jìn)行分析。圖2為模擬后不同壓下量下的等效應(yīng)力分布云圖，整個(gè)變形過程中等效應(yīng)力分布在66.8～83.0 MPa之間。從圖2(a)可以發(fā)現(xiàn)，在第30步時(shí)，試樣出現(xiàn)了一定程度的鼓肚，上下端面位置應(yīng)力較大，心部最小，在試樣鼓肚位置應(yīng)力均勻分布。圖2(b)中的應(yīng)力分布和圖2(a)中的分布相似，只是最小和最大應(yīng)力值有所提高。在圖2(c)中，試樣中的最小等效應(yīng)力分布在試樣鼓肚位置，與圖2(a)、(b)存在明顯差別，端面處的等效應(yīng)力有所提高，心部的應(yīng)力分布相對(duì)較高。

圖3為模擬后不同壓下量下的等效應(yīng)變分布云圖，從圖3(a)可以發(fā)現(xiàn)，在第30步時(shí)，上下端面位置應(yīng)變較大，心部最小，在試樣鼓肚位置應(yīng)力均勻分布。圖3(b)中的應(yīng)變分布有所增加，其最小和最大應(yīng)變值均有所提高。在圖3(c)中，試樣中應(yīng)變分布與圖3(b)相似，應(yīng)變值進(jìn)一步增加。

圖2 不同分析步下的等效應(yīng)力分布云圖(a)第30步；(b)第60步；(c)第91步

圖3 不同分析步下的等效應(yīng)變分布云圖(a)第30步；(b)第60步；(c)第91步

圖4 不同分析步下的損傷分布云圖(a)第30步；(b)第60步；(c)第91步

圖5 不同摩擦作用下的結(jié)果(a)鼓肚；(b)等效應(yīng)力；(c)等效應(yīng)變；(d)損傷最大值

圖4為不同壓下量的Cockcroft-Latham損傷結(jié)果，對(duì)比3個(gè)不同壓下量下結(jié)果可以發(fā)現(xiàn)，在試樣鼓肚及附近位置出現(xiàn)明顯損傷，特別是鼓肚位置的損傷值較大，而坯料上下端未見明顯損傷。隨著壓下量的增加，鼓肚位置的損傷值不斷增加。分析其原因，是由于在鼓肚位置由于沒有約束，處于自由變形狀態(tài)，在正的應(yīng)力三軸度的作用下，此處的損傷將不斷增加。當(dāng)達(dá)到材料的臨界損傷值時(shí)，試樣中將出現(xiàn)明顯的裂紋。

2.2 結(jié)果分析

通過后處理，將變形過程不同階段下的鼓肚(試樣上下端面直徑平均值d1和中間鼓肚位置的最大直徑d2之間的差值)、等效應(yīng)力、等效應(yīng)變及損傷結(jié)果進(jìn)行分析，具體結(jié)果如圖5所示。從圖5(a)所示結(jié)果可以發(fā)現(xiàn)，隨著加載的進(jìn)行，不同摩擦作用下的鼓肚變化趨勢(shì)基本一致，均是先增大后趨于平衡；而不同摩擦系數(shù)之間的差異在于，摩擦系數(shù)越高，所引起的鼓肚幅值也越大，分析原因，主要是鼓肚位置為徑向自由變形，而由于摩擦作用的存在，使得坯料和砧板結(jié)合處材料徑向流動(dòng)受阻，摩擦作用越大，材料受阻越嚴(yán)重，與鼓肚位置材料流動(dòng)差異也就越大，鼓肚就越明顯。圖5(b)中為試樣中平均值變化，誤差棒表示試樣中的最大值和最小值，可以發(fā)現(xiàn)不同摩擦作用下，平均應(yīng)力的變化不大，而摩擦作用越大，最大值和最小值之間的差異越大，試樣中的應(yīng)力均勻性相對(duì)降低。在圖5(c)中，隨著加載的進(jìn)行，試樣中的等效應(yīng)變逐漸增加，等效應(yīng)變最大值和最小值(誤差棒中的最大值和最小值)的差異也逐漸增大，說(shuō)明試樣中應(yīng)變均勻性也逐漸降低，且摩擦系數(shù)越大，應(yīng)變均勻性越差。圖5(d)中給出了基于Cockcroft-Latham準(zhǔn)則的損傷最大值結(jié)果，隨著變形的進(jìn)行，試樣中損傷最大值不斷增加，且摩擦的作用影響非常明顯，摩擦系數(shù)越大，損傷最大值越大。因此，從以上結(jié)果分析，可以發(fā)現(xiàn)，由于摩擦作用的存在，對(duì)材料變形、等效應(yīng)力、等效應(yīng)變及損傷均不同程度有影響，為了使得材料變形均勻及避免損傷，應(yīng)當(dāng)減少摩擦作用的影響。

3 結(jié) 論

為了分析摩擦作用對(duì)材料在熱變形過程中的影響，本文結(jié)合Deform 3D有限元軟件對(duì)其進(jìn)行研究，分別設(shè)置摩擦系數(shù)f為0.2、0.3、0.4和0.5，對(duì)變形過程中的變形鼓度、應(yīng)力、應(yīng)變及Cockcroft-Latham損傷結(jié)果進(jìn)行分析，發(fā)現(xiàn)隨著加載壓下量的增加，不同摩擦系數(shù)條件下試樣中的鼓度、應(yīng)力、應(yīng)變和損傷值均逐漸升高；而隨著摩擦系數(shù)的增加，所引起的鼓肚幅值也越大，應(yīng)力和應(yīng)變的最大值和最小值之間的差異越大，應(yīng)力應(yīng)變均勻性相對(duì)降低，且損傷最大值越大。