李波，楊祖建

(1.貴州師范大學機械與電氣工程學院，貴州貴陽550000;2.貴州航宇科技發(fā)展股份有限公司，貴州貴陽550000)

?

基于Simufact的GH3044異形截面環(huán)件軋制過程數值模擬*

李波1，楊祖建2

(1.貴州師范大學機械與電氣工程學院，貴州貴陽550000;2.貴州航宇科技發(fā)展股份有限公司，貴州貴陽550000)

摘要：利用有限元數值仿真分析軟件Simufact.forming對GH3044合金異形截面環(huán)坯的塑性成形過程進行動態(tài)模擬仿真，分析軋制過程中環(huán)件的溫度場、應變場和軋制力的分布情況及隨時間的變化規(guī)律。實時觀測環(huán)件在軋制過程中成形和缺陷生成的動態(tài)過程，獲得在各個軋制時刻各個節(jié)點位置的位移、應力、應變、速度的值和可視化云圖，其模擬結果對實際塑性加工過程有預見性和前瞻性。對產品的工藝參數的設計和優(yōu)化起到積極的指導作用。

關鍵詞:Simufact環(huán)件軋制數值模擬軋制工藝

0引言

環(huán)件軋制是一種常用的環(huán)形件塑性成形加工工藝，其工作原理如圖1所示。驅動輥在給定的轉速下繞軸線旋轉，通過摩擦力帶動環(huán)坯旋轉，芯輥在給定的進給速度下，平移擠壓環(huán)坯，同時在環(huán)坯摩擦力的作用下，被動旋轉，環(huán)坯在主動輥、芯輥、導向輥、錐輥的同時擠壓下，產生連續(xù)局部塑性變形，使環(huán)件壁厚減小、直徑擴大,直到厚度和直徑達到設計值后，軋制結束。其中錐輥的作用是在軸向進行擠壓，端面軋輥作旋轉端面軋制運動和軸向進給運動，保證環(huán)件的指定高度。導向輥為可自由轉動的從動輥，其作用是在軋制開始階段，對環(huán)坯進行定心，保證軋制過程平穩(wěn)進行，在軋制后期，讓環(huán)件歸圓，保證環(huán)件的精度。

圖1　徑軸向軋制設備圖

GH3044是Ni—Cr基固溶強化型變形高溫合金，使用溫度在900℃以下。合金中含有大量鎢元素和少量鋁和鈦元素。合金具有高的塑性和中等的熱強性，并具有優(yōu)良的抗氧化性能和良好的沖壓、焊接工藝性能。適宜制作在900℃以下長期工作的航空發(fā)動機主燃燒室和加力燃燒室零部件，以及隔熱屏、導向葉片等[1]。

本文采用Simufact有限元分析軟件，模擬了GH3044合金異形截面環(huán)件在一個生產周期內的軋制過程，并對軋制過程的溫度場、應變場及芯輥軋制力進行了分析研究。模擬仿真分析主要包括以下步驟[2]：

1)生成環(huán)件軋制的分析項目；

2)導入工件及模具幾何模型；

3)工件及模具位置定義；

4)定義工件及模具材料；

5)定義模具運動方式；

6)定義模具和工件的熱參數；

7)定義摩擦系數；

8)工件網格劃分；

9)工藝分析條件設置；

10)求解分析；

11)分析結果后處理。

1環(huán)件軋制幾何模型的建立

圖2　軋制模擬模型

軋制模擬幾何模型的建立如圖2所示，模型由主輥、芯輥、抱輥、錐輥和坯料組成。模型中各軋輥及毛坯均采用等分六面體自動掃描劃分網格，網格單元數為3 956。環(huán)件軋制溫度為1 140℃，模具預熱溫度為300℃。

2Simufact.forming數值仿真分析邊界條件

圖3　軋制零件網格模型

工件材料為GH3044，模具為模具鋼。

驅動輥轉速：2 rad/s，繞Z軸旋轉。

芯輥進給速度為1 mm/s,沿Y向進給。

模具溫度：300℃。

工件溫度：1 140℃。

環(huán)境參數：20℃。

網格類型：六面體網格。

網格單元數量：3 956。

3環(huán)件軋制模擬仿真求解

有限元模型準備好后，施加邊界條件，然后調用Simufact.forming求解器進行求解。

圖4　求解器設置圖

圖5　求解器參數設置圖

4軋件軋制模擬仿真求解后處理

Simufact模擬仿真過程中，可以通過云圖、圖表、曲線、動畫等方式把環(huán)件塑性變形的過程實時的動態(tài)顯示出來，以便觀察者可以直觀的觀測每一個時刻每一位置的溫度場、應變場、應力場、模具軋制力等的值。以下分別把溫度場、應變場和軋制力的分布情況及隨時間的變化過程進行分項分析。

4.1溫度場變化及其分布規(guī)律

取幾個具有代表性時刻的溫度場分布情況進行研究，即開始軋制1 s時刻、中間時刻5 s時刻、尾聲時刻8 s時刻和結束軋制10 s時刻的溫度場分布，如圖6-圖9所示。

通過圖6到圖9的溫度場分布云圖結果可以看出，隨著軋制過程的進行，環(huán)件的整體溫度逐漸上升，到軋制結束時刻，溫度達到了最高。這與實際的物理實驗得到的溫度分布情況吻合。導致溫度上升的原因是環(huán)件在塑性變形過程中，一方面變形能轉化為熱能，另外一方面是由于環(huán)件與模具接觸的地方，存在摩擦，摩擦產生的機械能轉換為熱量，也導致了環(huán)件的整體溫度上升。

圖6　t=1 s時刻的環(huán)件溫度場分布云圖

圖7　t=5 s時刻的環(huán)件溫度場分布云圖

圖8　t=8 s時刻的環(huán)件溫度場分布云圖

圖9　t=10 s時刻(軋制結束)的環(huán)件溫度場分布云圖

對環(huán)件取7個跟蹤點，觀測每個跟蹤點隨時間的變化規(guī)律及溫度的分布情況，如圖10。

圖10　采集的7個跟蹤點位置示意圖

圖11　跟蹤點的溫度場隨時間的變化曲線圖

從圖11可以看出，環(huán)件壁厚方向上，中間位置的跟蹤點P2、P4、P6的溫度最高，并且隨著軋制過程的繼續(xù)，在環(huán)件與模具咬入前，中部溫度平緩，咬入開始剛剛進入軋制時，溫度突然上升，中心點轉動繞過主動輥和芯輥的擠壓位置后，溫度平緩略下降，直到下一個咬入周期后，溫度又陡然上升，如此循環(huán)，直到軋制結束。而環(huán)件內壁和外壁的跟蹤點 P1、P3、P5、P7從空行程到咬入之間，溫度急劇下降，當咬入時刻，由于高溫環(huán)件與低溫模具接觸，產生熱傳導，溫度從高溫區(qū)環(huán)件外壁向低溫區(qū)主動輥傳遞，內部向低溫區(qū)芯輥傳遞，導致跟蹤點位置的溫度更加急劇下降，如此循環(huán)，直到軋制結束。在軸向方向上，上下端面的溫度比中部的溫度低，壁厚方向，內部溫度比外部高，分析結果與實際的物理實驗結果吻合。這種溫度分布現(xiàn)象是由于環(huán)件置放于空氣中，表面與空氣直接接觸，而環(huán)件的溫度比空氣溫度高，由于對流及輻射的作用，環(huán)件自身的熱量及軋制過程中變形能和摩擦產生的熱量都將與周圍的低溫空氣進行熱交換，從而帶走環(huán)件表面的熱量，導致環(huán)件表面的溫度能快速的降低[3]，而環(huán)件壁厚內部直接是金屬原子接觸，高度密集的固體原子不能快速的與外界空間進行熱交換，溫度梯度小，較差的散熱條件導致熱量不能被快速帶走，因此內部溫度較高[4]。模擬結果對實際塑性加工過程有預見性和前瞻性，對加工過程的工況有明顯的指導性。

4.2產生的應變大小及其分布規(guī)律

取幾個具有代表性時刻的應變場分布情況進行研究，即開始軋制1 s時刻、中間時刻5 s時刻、尾聲時刻8 s時刻和結束軋制10 s時刻的溫度場分布，詳見圖12。

圖12　t=1 s時刻環(huán)件等效應變場分布云圖

圖13　t=5 s時刻環(huán)件等效應變場分布云圖

圖14　t=8 s時刻環(huán)件等效應變場分布云圖

圖15　t=10 s(軋制結束)時刻環(huán)件等效應變場分布云圖

從圖12到圖15的等效應變的分布和變化云圖可知，開始軋制階段，應變主要發(fā)生在環(huán)件的外壁和內壁，因為外壁一開始就和主動輥、導向輥接觸，內壁一開始就和芯輥接觸，在主動輥和芯輥擠壓作用下，外壁和內壁與模具接觸區(qū)一開始就發(fā)生變形從而產生應變。由于環(huán)件內部與主動輥和芯輥的距離較遠，因此開始軋制階段，變形很小，隨著軋制的進行，變形越來越大，產生的應變也越來越大。在整個軋制過程中，外壁的應變一直大于內壁的應變，內壁的應變一直大于內部的應變。軋制結束，最大應變發(fā)生在環(huán)件小口的外壁處，達到0.641。

圖16　采集的7個跟蹤點位置示意圖

從圖17的跟蹤點的應變隨時間的變化圖可以看出，隨著軋制的進行，環(huán)件壁厚上面的應變越來越大。處于中部位置P2、P4、P6點的應變較小，處于外側的跟蹤點P1、P3、P5、P7的應變較大。應變在進入軋制區(qū)時刻，陡然上升，旋轉過了軋制區(qū)后，變得平緩，直到下一次轉到軋制區(qū)后，又陡然上升，如此循環(huán)，直到軋制結束。仿真的應變分布和隨時間的變化規(guī)律與實際的物理實驗分析結果吻合[5]。

圖17　跟蹤點的應變場隨時間的變化曲線圖

從圖中等效應變分布云圖中可以看出，變形初期只有環(huán)坯的內層和外層環(huán)面區(qū)域產生了塑性變形，隨著軋制過程的進行，塑性變形由環(huán)件內、外層向中層擴展，塑性變形區(qū)域逐漸擴大、變形程度逐漸增加，環(huán)件上等效應變也逐漸變大，但分布不均勻。軋制過程中環(huán)坯中間層由于遠離軋輥為小變形區(qū)，故等效應變也較小，而環(huán)件外層變形量始終大于內層，所以外層的等效應變大于內層[6]。

4.3產生的軋制力大小及與時間的變化規(guī)律

在軋制過程中，芯輥與環(huán)件直接接觸，并且進行平移進給，擠壓環(huán)件，導致芯輥上產生軋制力，其位置處于芯輥與環(huán)件的接觸區(qū)域，方向由芯輥指向主動輥。而驅動輥在旋轉過程中擠壓環(huán)件，也會產生軋制力矩，位置處于驅動輥的軸心位置，方向與主動輥旋轉方向相關。通過模擬仿真，可以實時的檢測到每一時刻芯輥的軋制力和驅動輥的軋制力矩的大小，以便指導生產過程中的軋制力和力矩的控制與調整[7]。

從圖18軋制過程中芯輥的軋制力隨時間的變化曲線圖和圖19的力矩圖可以看出，隨著軋制的進行，芯輥上面產生的軋制力越來越大，主動輥上面產生的軋制力矩也是越來越大。上升到一定的峰值后，又逐漸變小，并在一定范圍內波動。產生波動的主要原因是因為環(huán)坯的表面粗糙，與模具接觸的區(qū)域不均勻，呈動態(tài)的不均勻接觸，因此會導致軋制力與軋制力矩產生波動。由此可證明合金環(huán)件的軋制過程是一個非穩(wěn)態(tài)成形過程，由于軋制力對環(huán)件工藝參數的制定影響較大，因此軋制力的研究對企業(yè)在生產過程中驅動輥轉速、芯輥進給速度等主要工藝參數的制定具有很高的參考價值及一定的指導作用[8]。

圖18　軋制過程中芯輥的軋制力隨時間的變化曲線圖

圖19　軋制過程中主動輥的軋制力矩隨時間的變化曲線圖

5結論

圖20　GH3044環(huán)件實物圖

通過利用Simufact有限元數值分析軟件對GH3044合金環(huán)件徑軸向軋制過程進行模擬仿真，實時的模擬了環(huán)件的動態(tài)軋制過程，對軋制成形過程中的溫度場、應力場和模具的軋制力的分布和變化規(guī)律進行了可視化的研究。模擬結果表明，環(huán)件在咬入之前的空行程階段、咬入后的平穩(wěn)軋制階段和軋制結束后的三個階段的溫度場、應力場和軋制力有明顯的區(qū)別，而且隨著軋制時間的推移，按照一定的變化規(guī)律和趨勢進行變化。通過觀測這些規(guī)律，可以指導工程師對軋制工藝參數、模具設計、環(huán)件的材料、散熱條件進行調整和優(yōu)化，通過企業(yè)設計制造，完成了GH3044環(huán)件的產品生產，其性能指標滿足標準要求。環(huán)件軋制后的實物如圖20所示。

參考文獻

[1]中國金屬學會高溫材料分會.中國高溫合金手冊：上卷[M].北京：中國標準出版社，2012：118-130.

[2]劉勁松，肖寒，段小亮.Simufact在材料成形與控制工程中的應用[M].中國水利水電出版社，2012：125-143.

[3]許思廣，曹起驥，連家創(chuàng).環(huán)件軋制中的溫度分析[J].鍛壓技術，1993(6)34-38.

[4]吳茂加，朱德芳，胡德政，等.鈦及其合金在熱軋過程中的溫降模型[J].稀有金屬材料與工程，1996,25(5):28-29.

[5]錢東升.異形截面環(huán)件冷軋力學原理和工藝理論研究[D].武漢：武漢理工大學，2009：39-48.

[6]李繞.金屬塑性成形原理[M].機械工業(yè)出版社，2010:56-82.

[7]王廷浦，齊克敏.金屬塑性加工學-軋制理論與工藝[M].冶金工業(yè)出版社，2012:50-85.

[8]華林.環(huán)件閉式軋制力和力矩上限計算[J].力學與實踐，1994,16(3)：39-43.

楊祖建(1984-)，男，貴州省貴陽市人，貴州航宇科技發(fā)展股份有限公司，工程師，主要從事金屬材料塑性成形研究。

Numerical simulation of the rolling process of GH3044 ring with special-shaped section based on Simufact

LI Bo, YANG Zujian

Abstract:In this paper, we carried out dynamic simulation of the rolling process of GH3044 ring with special-shaped section using finite element numerical simulation software Simufact, analyzed the distribution and changing rules of the temperature field, the strain field and the rolling force of the ring, observed the dynamic process of the forming of the ring and the appearing of flaws in real time, and obtained the values and visualized cloud pictures of the displacement, the stress, the strain and the speed of each node at each rolling moment. The results of the simulation provided foresight to the actual rolling process, and guidance on the design and optimization of the technical parameters of the products.

Keywords:Simufact; ring rolling; numerical simulation; rolling technique

收稿日期：2015-09-13

作者簡介：李波(1965-)，女，遼寧省丹東市人，就職于貴州師范大學機械與電氣工程學院，碩士研究生，副教授，主要從事金屬材料塑性成形研究。

基金項目：貴陽國家高新區(qū)(省級)科技合作計劃項目(GXKH2012-026)。

中圖分類號：TG337.1

文獻標識碼：A

文章編號：1002-6886(2016)02-0022-06