Ti6Al4V合金粉末高溫高壓成形過程中粉固界面及其耦合變形研究

(北京航空航天大學(xué) 機(jī)械工程及自動(dòng)化學(xué)院，北京，100083)

鈦合金具有較高的比強(qiáng)度，同時(shí)鈦合金的耐熱性、耐蝕性、抗彈性良好[1]，廣泛應(yīng)用于航空航天、軍事和船舶等領(lǐng)域。傳統(tǒng)加工方式對(duì)鈦合金的加工存在一些問題，如鑄造會(huì)產(chǎn)生縮孔、縮松，鍛造后續(xù)機(jī)加工周期長等。熱等靜壓作為粉末近凈成形的一種有效方法，克服了鈦合金傳統(tǒng)成形工藝上的不足[2]。熱等靜壓技術(shù)是一種在高溫、高壓條件下通過包套對(duì)粉末材料進(jìn)行直接固結(jié)成形的技術(shù)。熱等靜壓產(chǎn)品致密度高、力學(xué)性能優(yōu)越，極大提高了材料的利用率，成為鈦合金、鋁合金等零件成形的最有效方法之一[3]。熱等靜壓過程中，粉末顆粒在壓力與溫度的共同作用下轉(zhuǎn)變?yōu)楦叨戎旅芑w，正是由于這種轉(zhuǎn)變，粉固界面的接觸狀態(tài)在熱等靜壓過程中的每一步都是不明確的，而且粉末與包套之間的摩擦阻礙了最優(yōu)致密度的形成。摩擦對(duì)于粉末成形過程的理解和分析具有重要的作用。相關(guān)研究報(bào)告顯示，許多粉末制件的最終性能強(qiáng)烈依賴于粉末與壁面的摩擦[4]，而且熱等靜壓件的一些缺陷例如裂紋、包套破裂等的產(chǎn)生與粉固界面的接觸狀態(tài)密切相關(guān)。粉末熱等靜壓成形過程是一個(gè)高度復(fù)雜的非線性化過程[5]，難以通過數(shù)值解析方法控形，而采用反復(fù)試驗(yàn)的方法不僅費(fèi)時(shí)，而且極大地增加了成本。數(shù)值模擬成為研究粉末熱等靜壓成形過程的一種快速、有效的手段。而在數(shù)值模擬過程中，粉固界面的接觸處理是關(guān)鍵一步，粉固界面接觸對(duì)嚴(yán)重影響著數(shù)值模擬結(jié)果，不當(dāng)?shù)奶幚韲?yán)重降低數(shù)值模擬的可靠性，失去其指導(dǎo)意義。目前，針對(duì)熱等靜壓過程中粉固界面的研究鮮有報(bào)道，熱等靜壓過程中包套與粉末的接觸受到多方面的影響，這種接觸的具體模型也很難從實(shí)驗(yàn)中直接得出，王剛[6]認(rèn)為熱等靜壓過程中2A12合金粉末與包套之間沒有相對(duì)位移，不必考慮包套與粉末的摩擦；劉國承等[7]采用庫侖摩擦模型模擬316L粉末與包套的接觸，并得到了近似尺寸；瞿宗宏等[8]采用粘滑摩擦模型處理FGH4097合金粉末與包套之間的接觸并得到了準(zhǔn)確的成形尺寸及變形趨勢，而關(guān)于Ti6Al4V合金粉固界面的研究，還沒有相關(guān)報(bào)道。因此，有必要展開對(duì)Ti6Al4V粉末與包套之間接觸以及摩擦力的研究。本文作者針對(duì)Ti6Al4V合金成形開展數(shù)值模擬研究工作，建立熱等靜壓力學(xué)模型，分析粉固界面的接觸狀態(tài)，研究粉固界面接觸狀態(tài)對(duì)熱等靜壓成形過程的影響，對(duì)Ti6Al4V合金熱等靜壓包套設(shè)計(jì)及工藝參數(shù)優(yōu)化具有重要意義。

1 接觸模型及材料本構(gòu)方程

1.1 接觸模型

1.1.1 Bilinear接觸模型

Bilinear模型是庫侖摩擦模型在有限元中的表現(xiàn)形式，庫侖摩擦模型是目前塑性成形過程中常用的模型之一[9]，由 COULOMB 在 AMONTONS 和LEONARDO 的基礎(chǔ)上，提出“摩擦力與作用在接觸面上的正壓力成正比，且與相對(duì)滑移速度無關(guān)”，進(jìn)一步發(fā)展成庫侖摩擦模型[10]。

式中：σn為接觸點(diǎn)法向應(yīng)力；σfr為切向應(yīng)力；μ為摩擦因數(shù)；t為相對(duì)滑動(dòng)速度方向上的切向單位矢量。

當(dāng)法向力給定后，摩擦力出現(xiàn)階梯函數(shù)狀變化，數(shù)值模擬計(jì)算中引入這種不連續(xù)性，往往會(huì)導(dǎo)致計(jì)算困難，因此，采用1個(gè)滑動(dòng)面來消除這種這種不連續(xù)性，Bilinear模型圖如圖1所示。圖1中：ft為剪切力；為切向位移增量；δ為黏性極限距離。

圖1 Bilinear模型示意圖Fig.1 Sketch map of Bilinear friction model

1.1.2 Stick-slip接觸模型

Stick-slip接觸模型是基于庫侖摩擦模型的針對(duì)粘滑摩擦狀態(tài)的修正模型，假設(shè)初始狀態(tài)時(shí)處于黏結(jié)狀態(tài)，經(jīng)過一定切向位移轉(zhuǎn)變?yōu)榛瑒?dòng)摩擦狀態(tài)，能夠?qū)崿F(xiàn)黏性摩擦到滑動(dòng)摩擦的突變，如圖2所示。圖2中：Ft為摩擦力；Fn為法向力；α為乘子，一般取1.05；β為滑動(dòng)摩擦到黏性摩擦的相對(duì)位移過度區(qū)域，一般取10-6。

圖2 Stick-slip模型示意圖Fig.2 Sketch map of Stick-slip friction model

1.1.3 Glue模型

Glue接觸模型將接觸處理為始終黏結(jié)狀態(tài)，在數(shù)值模型中表現(xiàn)為接觸節(jié)點(diǎn)之間固連，沒有相對(duì)位移產(chǎn)生。

1.2 本構(gòu)模型

粉末致密化過程實(shí)質(zhì)上是由非連續(xù)體轉(zhuǎn)變?yōu)檫B續(xù)體的過程，在這個(gè)過程中體積不斷發(fā)生變化。非連續(xù)體的變形非常復(fù)雜，需要考慮所有顆粒的應(yīng)力與應(yīng)變[6]，因此，將粉末看作“連續(xù)的多孔材料”，應(yīng)用連續(xù)介質(zhì)理論描述粉末體。

在經(jīng)典Von Mises理論基礎(chǔ)上，SHIMA等[11]考慮多孔材料相對(duì)密度和靜水壓力的影響，提出Shima屈服模型。

式中：σ1，σ2和σ3為3個(gè)主應(yīng)力；f為靜水壓力對(duì)屈服的影響程度；σm為靜水壓力；F為屈服應(yīng)力。本文采用MSC.marc軟件對(duì)成形過程進(jìn)行模擬，軟件中嵌入的是修正后的Shima模型：

式中：σy為單軸屈服應(yīng)力；σd為應(yīng)力偏張量；p為靜水壓力；β和γ為相對(duì)密度的函數(shù)，

單軸壓縮試驗(yàn)可獲得β和γ與ρ之前的關(guān)系數(shù)據(jù)，進(jìn)而獲得q1～q4和b1～b4共8個(gè)系數(shù)的數(shù)值。

2 結(jié)果與分析

試驗(yàn)采用的Ti6Al4V預(yù)合金粉末由陜西海寶特種金屬材料有限公司提供，粉末粒徑在60～300 μm之內(nèi)，包套材料為不銹鋼。Ti6Al4V的化學(xué)成分如表1所示。

熱等靜壓設(shè)備為ABB公司的QIH-5熱等靜壓機(jī)，最高溫度可達(dá)2 000 ℃，最高壓力為200 MPa。包套材質(zhì)為301不銹鋼，厚度為2 mm。為保證粉末充填密實(shí)，充填過程在VSR-200時(shí)效振動(dòng)系統(tǒng)上完成，振動(dòng)頻率為30 Hz，振動(dòng)時(shí)間1 h，振動(dòng)后粉末振實(shí)相對(duì)密度為0.65。為充分排出實(shí)驗(yàn)件內(nèi)空氣，在500 ℃下進(jìn)行高溫脫氣8 h，使得包套內(nèi)真空度達(dá)到10-5Pa，隨后，采用電子束焊密封包套。將已經(jīng)脫氣的包套進(jìn)行熱等靜壓壓制，熱等靜壓溫度為 920 ℃，壓力為120 MPa，升溫升壓、保溫保壓、降溫降壓時(shí)間分別為2，3和2 h。

表1 Ti6Al4V合金粉末的化學(xué)成分(質(zhì)量分?jǐn)?shù))Table1 Chemical composition of Ti6Al4V aluminum alloy powders %

Ti6Al4V合金粉末試件為軸對(duì)稱圓柱試件，所以選取模型的 1/8為計(jì)算單元，利用大型有限元軟件MSC.marc進(jìn)行模擬計(jì)算，并將粉末試件相對(duì)密度的模擬計(jì)算結(jié)果與試驗(yàn)結(jié)果進(jìn)行對(duì)比，結(jié)果如圖3所示。

試驗(yàn)結(jié)果如圖3(a)所示。在熱等靜壓過程中，盡管包套外表面都處在等靜壓載荷作用下，但由于包套對(duì)壓力具有屏蔽效應(yīng)[12]，包套內(nèi)部粉末則處在不等載荷的作用下，這也是試驗(yàn)件內(nèi)部產(chǎn)生相對(duì)密度梯度的原因。這種屏蔽效應(yīng)在包套和頂部直壁所組成的“L”狀結(jié)構(gòu)處尤為明顯，因此，該區(qū)域相對(duì)密度較低。其他區(qū)域粉末收縮比較明顯，致密化程度較高[13]。

由圖3可以看出：Bilinear模型與Stick-slip模型得到的坯料變形規(guī)律與HIP試驗(yàn)結(jié)果高度相符，Glue模型對(duì)坯料邊角處變形的預(yù)測結(jié)果與試驗(yàn)結(jié)果存在一定差距。在熱等靜壓過程中，隨著壓力與溫度上升，粉末體收縮，包套產(chǎn)生相應(yīng)變形。采用Glue模型時(shí)，粉末與包套之間處于一種黏結(jié)狀態(tài)，在伴隨包套發(fā)生變形的同時(shí)，由于這種黏結(jié)效應(yīng)，粉末產(chǎn)生一定程度的塑性變形；而對(duì)于Blinear模型與Stick-slip模型，粉末與包套處于摩擦接觸狀態(tài)，粉末與包套產(chǎn)生一定程度的相對(duì)滑動(dòng)。Glue模型的這種粉末塑性變形對(duì)于外部等靜壓的阻礙作用遠(yuǎn)遠(yuǎn)大于摩擦力的阻礙作用，因此，這種阻礙作用是與包套變形程度直接相關(guān)的。在“L”狀結(jié)構(gòu)處，包套的變形程度較大，這種阻礙作用尤為明顯。在其他區(qū)域，包套變形程度較低，阻礙作用較小。這也是 3種模型在“L”狀區(qū)域變形程度差異較大，而軸向與徑向尺寸差異較小的主要原因。表2所示為采用3種不同接觸模型的熱等靜壓數(shù)值模擬結(jié)果與試驗(yàn)結(jié)果的對(duì)比。

圖3 相對(duì)密度試驗(yàn)及模擬結(jié)果圖Fig.3 Test and simulation results of relative density

表2 參數(shù)模擬結(jié)果與試驗(yàn)結(jié)果的對(duì)比Table 2 Comparison between simulation results and test result of parameters

由表2可知：參數(shù)數(shù)值模擬結(jié)果和試驗(yàn)數(shù)據(jù)高度一致，最大相對(duì)誤差不超過 5%，在工程允許得誤差范圍之內(nèi)，說明數(shù)值模擬在預(yù)測Ti6Al4V合金粉末熱等靜壓致密化過程方面具有一定的指導(dǎo)意義。

2.1 粉固界面的接觸理論分析

目前，還沒有針對(duì)熱等靜壓過程中粉固界面的接觸理論研究，而針對(duì)粉末壓實(shí)過程中粉固界面的接觸研究較多，可采用壓制過程的理論分析模型計(jì)算摩擦力[14]，而對(duì)于圓柱粉末熱等靜壓件，由于包套外部受到等靜壓作用，而且在熱等靜壓后，包套徑向尺寸收縮僅為10%，因此，在粉末壓制圓柱模型側(cè)曲面施加等同于頂部平面的壓力，這種修改后的力學(xué)狀態(tài)等同于熱等靜壓圓柱模型1/2結(jié)構(gòu)的力學(xué)狀態(tài)。由于熱等靜壓圓柱件側(cè)曲面整體上收縮程度基本保持一致，而平面將會(huì)出現(xiàn)中心凹陷等不均勻收縮，因此，這種理論應(yīng)用僅局限于簡單圓柱模型?？蓪A柱粉末單向壓實(shí)理論引入熱等靜壓并針對(duì)熱等靜壓受力、變形的特點(diǎn)進(jìn)行修正補(bǔ)充。

圓柱粉末件的簡化模型及受力狀態(tài)如圖4所示。圖4中：Q為粉末體橫截面積；h為1/2粉末體高度；p為熱等靜壓的壓強(qiáng)。從距離上表面z處取一厚度為dz薄層進(jìn)行受力分析。由于圓柱體表面受到相等壓力，圓柱體粉末與包套的接觸狀態(tài)關(guān)于中心面對(duì)稱，因此，取1/2圓柱體進(jìn)行研究。

圖4 粉末體受力示意圖Fig.4 Force diagram of powder body

在熱等靜壓過程中，粉體受力平衡。上表面所受到的壓應(yīng)力pu=p，壓力fu=Qp，因此，薄層上表面的軸向?qū)嶋H致密化力為

式中：fμ為包套與粉末之間的摩擦力。

軸向?qū)嶋H致密化力的增量為

摩擦力正比于徑向壓力，

式中：M為圓柱體的包絡(luò)面面積。

考慮徑向應(yīng)力與軸向應(yīng)力的比值，令pr/pa=k，代入式(8)得

式中：C為圓柱體在z處的周長，由f(0)=fu，對(duì)式(9)進(jìn)行積分得：

k為徑向壓力與軸向壓力的比值，它是材料和幾何的相關(guān)參數(shù)。對(duì)于純彈性圓柱體，進(jìn)行軸向加載時(shí)，徑向會(huì)產(chǎn)生膨脹，根據(jù)胡克定律有

但是，彈性變形不是熱等靜壓過程中粉末致密化的主要機(jī)制，因此，由式(12)得到的固定值k不適用于熱等靜壓。

粉末致密化過程是一個(gè)高度非線性化的過程，由于塑形變形是粉末致密化的主要機(jī)制，因此，可將塑形數(shù)學(xué)理論應(yīng)用于粉末致密化。粉末壓制已經(jīng)采用該理論取得較好的結(jié)果[4]。針對(duì)粉末壓制過程中所得k進(jìn)行修正，可應(yīng)用于熱等靜壓過程中。

在粉末壓制過程中，粉末在軸向進(jìn)行壓縮時(shí)，徑向會(huì)產(chǎn)生膨脹的趨勢，因此，粉末壓制過程中的徑向力就是抑制這種趨勢而產(chǎn)生的力。而對(duì)于圓柱體熱等靜壓，其徑向尺寸在熱等靜壓前后變化僅為10%，可近似采用粉末壓制理論。包套軸向與徑向受到相同的壓力，粉末徑向就受到膨脹產(chǎn)生的力和外部壓力的雙重作用，因此，可得熱等靜壓過程中的k。

熱等靜壓過程中的k為

式中：σro為外部壓應(yīng)力；σre為徑向膨脹而產(chǎn)生的壓應(yīng)力。

由于σro≈σa，因此，

式中：kc為粉末壓制過程中的徑向力與軸向力的比值。

根據(jù)式(14)及粉末壓制過程中相關(guān)理論[4]，可獲得熱等靜壓過程中徑向力與軸向力比值隨相對(duì)密度的變化曲線，如圖5所示。

在數(shù)值模擬中所給定的摩擦因數(shù)往往是一個(gè)固定值，為了探究粉固界面的實(shí)際摩擦因數(shù)，針對(duì)Stick-slip與Bilinear模型(Glue模型假設(shè)為固聯(lián)狀態(tài)，不存在摩擦)的數(shù)值模擬結(jié)果進(jìn)行分析，將相關(guān)數(shù)據(jù)代入式(11)，得到摩擦因數(shù)隨時(shí)間變化關(guān)系。由于降溫降壓階段粉末已達(dá)到高度致密化狀態(tài)，此時(shí)的接觸為固體金屬之間的接觸，故只取升溫升壓階段及保溫保壓階段進(jìn)行分析，如圖6所示。

圖5 徑向力與軸向力比值Fig.5 Ratio of radial force to axial force

圖6 實(shí)際摩擦因數(shù)Fig.6 Actual friction coefficient

由圖6可以看出：采用2種接觸模型數(shù)值模擬結(jié)果所得實(shí)際摩擦因數(shù)差異較大；在升溫升壓階段，采用Stick-slip模型所得實(shí)際摩擦因數(shù)具有較大初始值，并且隨時(shí)間而降低，在4 000 s左右達(dá)到最小值，此時(shí)，粉末相對(duì)密度達(dá)0.9，粉末致密化過程基本完成；而采用 Bilinear模型所得實(shí)際摩擦因數(shù)，其初始值較小，隨著壓力與溫度升高而增加并在2 000 s達(dá)到最大值，然后出現(xiàn)下降趨勢，并在5 000 s后趨于穩(wěn)定；在保溫保壓階段，這2種模型所得的實(shí)際摩擦因數(shù)基本穩(wěn)定，采用 Stick-slip模型的實(shí)際摩擦因數(shù)為 0.15，采用Bilinear模型的實(shí)際摩擦因數(shù)為0.10。

粉固界面的摩擦狀態(tài)是由兩者之間的相互作用決定的，這種相互作用主要表現(xiàn)為粉末顆粒變形與滑移，在熱等靜壓過程中是不斷發(fā)生變化的。在熱等靜壓前，粉末顆粒與包套的實(shí)際接觸接近于點(diǎn)接觸，因此，粉固界面總的真實(shí)接觸面積非常小，當(dāng)施加壓力后，這些微小的接觸區(qū)域便承受巨大的應(yīng)力而發(fā)生黏附現(xiàn)象，使得熱等靜壓起始時(shí)具有較高的摩擦因數(shù)[15]；隨著壓力的不斷增加，粉末顆粒的接觸面積不斷增大，材料的硬化會(huì)產(chǎn)生一個(gè)更穩(wěn)定的接觸面[8]，進(jìn)而粉末與模具之間的接觸面上的孔隙逐漸被填平，使得接觸面逐漸平滑，從而摩擦因數(shù)降低[16]。伴隨著相對(duì)密度不斷提高，粉末由松散顆粒轉(zhuǎn)變?yōu)橹旅荏w，與包套的接觸狀態(tài)更接近于金屬實(shí)體的接觸，因而，在保溫保壓階段摩擦因數(shù)趨于1個(gè)定值。2個(gè)模型在升溫升壓階段的結(jié)果差異較大，升溫升壓階段粉末相對(duì)密度較低。而當(dāng)粉末相對(duì)密度較低時(shí)，粉末滑移可被看作粉末在靠近包套的邊界層中的剪切作用以及在粉固界面上的滑動(dòng)，這個(gè)過程有助于粉末在熱等靜壓早期階段的顆粒重排，并且對(duì)于復(fù)雜結(jié)構(gòu)的成形起著重要的作用。因此，綜合根據(jù)以上分析，Stick-slip接觸模型能更好地反映出摩擦因數(shù)的變化過程。

2.2 接觸模型對(duì)相對(duì)密度的影響

由表2可以看出3種模型均能實(shí)現(xiàn)對(duì)整體幾何尺寸的預(yù)測。對(duì)于粉末制品來說，其最終整體成形密度與硬度、疲勞強(qiáng)度、韌性等性能密切相關(guān)，而密度分布的均勻性與缺陷、裂紋等因素密切相關(guān)，這些因素的存在甚至?xí)?dǎo)致熱等靜壓件因性能低劣而廢棄[16]。因此，在可靠預(yù)測幾何尺寸的前提下，預(yù)測最終密度的準(zhǔn)確度顯得尤為重要。在預(yù)測最高相對(duì)密度方面，Stick-slip模型與Bilinear模型均能實(shí)現(xiàn)精準(zhǔn)預(yù)測，而Glue模型預(yù)測效果較差。下面選取幾個(gè)典型位置的相對(duì)密度預(yù)測進(jìn)行分析，各節(jié)點(diǎn)位置如圖7所示。

熱等靜壓一般用于制造輪廓復(fù)雜的零件，節(jié)點(diǎn)A位于粉末體外沿，位于粉末體外廓急劇變化的部位，能夠代表復(fù)雜結(jié)構(gòu)部位粉末相對(duì)密度的變化趨勢；節(jié)點(diǎn)B位于粉末體中心位置，能夠代表內(nèi)部粉末相對(duì)密度的變化趨勢。

節(jié)點(diǎn)A和B處相對(duì)密度隨時(shí)間變化的曲線如圖8所示。從圖8可知：2個(gè)位置的相對(duì)密度都是呈現(xiàn)上升趨勢，但在熱等靜壓初期，相對(duì)密度有下降趨勢，這是由于隨著溫度上升，粉末體膨脹而產(chǎn)生的膨脹量大于壓力產(chǎn)生的壓縮量；隨著壓力的升高，壓縮量越來越大，相對(duì)密度逐漸上升。

圖7 FEM模型的節(jié)點(diǎn)位置Fig.7 Nodes positions in FEM model

圖8 熱等靜壓過程中相對(duì)密度變化Fig.8 Relative density changes in HIP

熱等靜壓完成后，這3種模型均能準(zhǔn)確預(yù)測節(jié)點(diǎn)B處的相對(duì)密度變化趨勢以及最終相對(duì)密度。而對(duì)于節(jié)點(diǎn)A處最終相對(duì)密度的預(yù)測結(jié)果，Glue模型僅為0.75，Stick-slip模型為0.80，Bilinear模型為0.83。在節(jié)點(diǎn)A區(qū)域進(jìn)行取樣，并采用阿基米德排水法[17]測得相對(duì)密度為0.88，因?yàn)橛晒?jié)點(diǎn)A附近區(qū)域的相對(duì)密度均比節(jié)點(diǎn)A的高，因此，測得值稍高于實(shí)際值。材料相對(duì)密度決定著材料的靜態(tài)強(qiáng)度、塑性、動(dòng)態(tài)強(qiáng)度、硬度等各項(xiàng)性能，而在相對(duì)密度為0.75時(shí)，對(duì)照相應(yīng)實(shí)驗(yàn)曲線[18]，平面應(yīng)力、抗拉強(qiáng)度、伸長率等參數(shù)均顯著低于正常值，此時(shí)，Ti6Al4V合金粉末未完全致密化，這與實(shí)際情況是不相符的。熱等靜壓多用于制備輪廓復(fù)雜零件，數(shù)值模擬對(duì)于以節(jié)點(diǎn)A為代表的結(jié)構(gòu)急劇變化位置的預(yù)測準(zhǔn)確性顯得尤為重要。由此可見，Glue模型所預(yù)測的最終相對(duì)密度與實(shí)際情況存在較大誤差。

2.3 接觸模型對(duì)粉末體、包套之間相對(duì)位移的影響

包套與粉末體的相對(duì)位移對(duì)于外輪廓的形狀變化具有一定影響作用，在熱等靜壓中，存在著包套與坯料的相對(duì)位移過大而導(dǎo)致包套破裂的問題，通過對(duì)相對(duì)位移的預(yù)測能在一定程度上避免出現(xiàn)上述情況。包套與粉末體相對(duì)位移受到壓力、溫度、熱膨脹系數(shù)差異等多方面的影響，由節(jié)點(diǎn)C和D，E和F，G和H之間的徑向位移差得到相對(duì)位移，如圖9所示。從圖9可知：粉末體外沿的相對(duì)位移較大。因此，取節(jié)點(diǎn)C和D在3種模型下的數(shù)值模擬結(jié)果進(jìn)行分析，如圖10所示。

圖9 相對(duì)位移隨時(shí)間變化Fig.9 Variation of relative displacement with time

總體來看，Glue模型不能預(yù)測相對(duì)位移；Stickslip模型與 Bilinear模型的預(yù)測結(jié)果均呈現(xiàn)上升趨勢并最終趨于穩(wěn)定；在HIP升溫升壓階段，Bilinear模型與 Stick-slip模型預(yù)測的相對(duì)位移出現(xiàn)明顯上升趨勢，粉末的致密化主要發(fā)生在這個(gè)階段，相對(duì)位移的這種變化趨勢主要是粉末致密化的作用結(jié)果。在保溫保壓階段，相對(duì)位移趨于穩(wěn)定，主要是由于這個(gè)階段發(fā)生蠕變，粉末體相對(duì)密度變化較小[19]。在降溫降壓階段，Stick-slip模型預(yù)測的相對(duì)位移出現(xiàn)了上升趨勢，在降溫降壓階段，由于粉末材料與包套材料熱膨脹系數(shù)不同，當(dāng)溫度下降相同時(shí)產(chǎn)生的收縮量不同，從而導(dǎo)致一定相對(duì)位移的產(chǎn)生。因此，Glue模型不能預(yù)測相對(duì)位移，Stick-slip與Bilinear模型預(yù)測的最終相對(duì)位移基本一致，Stick-slip模型能更好地反映相對(duì)位移的具體變化趨勢。

圖10 節(jié)點(diǎn)C和D相對(duì)位移Fig.10 Relative displacement between nodeCandD

3 結(jié)論

1) 基于Shima模型的有限元技術(shù)在預(yù)測Ti6Al4V合金熱等靜壓零件變形以及相對(duì)密度方面具有較高的精度，徑向尺寸、軸向尺寸、相對(duì)密度的相對(duì)誤差均在 5%以內(nèi)。因此，數(shù)值模擬可以為熱等靜壓零件成形提供參考依據(jù)，實(shí)現(xiàn)精確控形。

2) 粉固界面之間的摩擦因數(shù)處于一個(gè)復(fù)雜變化的過程。初始階段摩擦因數(shù)最大，隨著粉末的致密化而逐漸減小，最終趨于1個(gè)穩(wěn)定值。Stick-slip模型能反映出這個(gè)變化趨勢，初始摩擦因數(shù)為0.27，穩(wěn)定摩擦因數(shù)為0.15。

3) Bilinear模型與Stick-slip模型能夠準(zhǔn)確預(yù)測相對(duì)密度變化趨勢及數(shù)值；Glue模型不能預(yù)測相對(duì)位移，Stick-slip模型能夠預(yù)測降溫降壓階段相對(duì)位移的具體變化趨勢。