張 順，劉小剛

（南京航空航天大學(xué)江蘇省航空動(dòng)力系統(tǒng)重點(diǎn)實(shí)驗(yàn)室，南京210016）

0 引言

超塑性成形（SPF）具有塑性高、變形抗力小、可一次精密成形等優(yōu)點(diǎn)，已在航空航天等工程領(lǐng)域得到廣泛應(yīng)用。其中，TC4鈦合金以其優(yōu)異的超塑性性能成為航空工業(yè)重要的超塑性成形件制造材料[1]。例如，超塑性成形/擴(kuò)散連接（SPF/DB）組合工藝成形的TC4合金寬弦空心風(fēng)扇葉片已成為國(guó)外先進(jìn)航空發(fā)動(dòng)機(jī)關(guān)鍵技術(shù)之一。鈦合金材料在超塑性成形時(shí)，需要在一定溫度、一定壓力下保持一段時(shí)間，成形工藝較為復(fù)雜。合理的工藝參數(shù)設(shè)置對(duì)獲得變形均勻的高質(zhì)量件有重要影響。采用數(shù)值模擬方法，研究材料超塑性貼膜過(guò)程及成形后的變形規(guī)律，可為零件質(zhì)量控制提供依據(jù)；此外，還可以為工藝設(shè)置提供參考，大大降低試制成本。

國(guó)內(nèi)外研究者針對(duì)材料超塑性成形過(guò)程開(kāi)展了一些數(shù)值模擬工作。Bonet J等[2]介紹了幾種模擬3維厚板零件超塑性成形所需的數(shù)值模擬技術(shù)，給出了超塑性行為的黏塑性特征方程；施曉琦等[3]結(jié)合有限元分析軟件MARC對(duì)單層板加強(qiáng)結(jié)構(gòu)和2層板盒形件結(jié)構(gòu)的超塑性成形過(guò)程進(jìn)行模擬,對(duì)構(gòu)件的形狀和尺寸進(jìn)行了預(yù)測(cè),獲得了最佳應(yīng)變速率下加載的壓力-時(shí)間曲線；趙毅等[4]對(duì)SPF/DB成形原理和過(guò)程進(jìn)行討論，利用ANSYS有限元軟件的建模和數(shù)值模擬,提取夾層板作為研究對(duì)象,建立模型，定義材料性能參數(shù)以及接觸和加載方式；黃鋼華[5]利用MARC對(duì)TC4鈦合金超塑性成形與超塑性成形/擴(kuò)散連接進(jìn)行數(shù)值模擬,預(yù)測(cè)成形件的壁厚分布情況,并獲得了最大等效應(yīng)變速率恒定條件下的優(yōu)化加載曲線；崔元杰[6]借助MARC分析了面板與芯板初始厚度比值r1以及擴(kuò)散連接寬度與面板初始厚度比值r2對(duì)TC4合金3層板“溝槽”缺陷的影響規(guī)律；Jiang[7]采用MARC軟件設(shè)計(jì)不同形狀的預(yù)成形模具，利用有限元模擬方法研究預(yù)成形模具和終成形模具的表面摩擦因數(shù)對(duì)筒形件超塑性成形側(cè)壁厚度分布的影響規(guī)律；K Kibble[8]、DC Chen[9]分別利用ABAQUS和DEFORM[TM]3D軟件，針對(duì)TC4合金圓柱模型以及橢圓柱模型進(jìn)行超塑性成形數(shù)值模擬，同時(shí)對(duì)成形之后的結(jié)構(gòu)件的厚度分布進(jìn)行預(yù)測(cè)。

綜上，研究者對(duì)一些杯狀、薄壁圓筒狀等簡(jiǎn)單形體的超塑性成形過(guò)程進(jìn)行了數(shù)值模擬，并成功預(yù)測(cè)了成形件的壁厚分布。本文在此基礎(chǔ)上，采用MARC軟件通過(guò)控制最大應(yīng)變速率的方法對(duì)TC4合金空心夾芯結(jié)構(gòu)SPF/DB過(guò)程進(jìn)行數(shù)值模擬，通過(guò)合理控制參數(shù)以獲得減薄率均勻的超塑性成形件；研究不同工藝參數(shù)對(duì)成形過(guò)程的影響規(guī)律，從而為寬弦空心葉片結(jié)構(gòu)的成形提供參考。

1 超塑性成形基本原理

目前通常用2種定義方式來(lái)描述超塑性：從拉伸試驗(yàn)的延伸率角度來(lái)定義，認(rèn)為δ＞200%即為超塑性；從應(yīng)變速率感應(yīng)系數(shù)m角度來(lái)定義，只要滿足m＞0.3，就是超塑性過(guò)程。當(dāng)然還有學(xué)者認(rèn)為抗頸縮能力強(qiáng)，就可以認(rèn)為是超塑性[1]。

超塑性成形的產(chǎn)生要求有一定的變形溫度和較小的應(yīng)變速率，同時(shí)合金本身還需要具有直徑小于5 μm的等軸晶粒，這種超塑性稱為超細(xì)晶粒超塑性；有些材料在特定溫度下，相結(jié)構(gòu)會(huì)發(fā)生轉(zhuǎn)變，此時(shí)發(fā)生的超塑性稱為相變超塑性[10]。本文中TC4合金的超塑性變形屬于前者。

在超塑性成形過(guò)程中變形板材作為1個(gè)非牛頓黏性的材料[11]，其基體的流變應(yīng)力為

式中：K為材料常數(shù)；m為應(yīng)變率敏感性指數(shù)；n為應(yīng)變硬化指數(shù)；˙和分別為有效應(yīng)變和基質(zhì)的有效應(yīng)變速率。

本文中3層板材料為T(mén)C4細(xì)晶鈦合金，蔡云[12]系統(tǒng)總結(jié)了TC4合金最佳超塑溫度在900℃附近，最佳變形速率在9.8×10-4附近，采用恒應(yīng)變速率法和速度突變法對(duì)m值進(jìn)行求解，求得能夠準(zhǔn)確描述TC4合金的m值分別為0.54和0.55。所以，在此基礎(chǔ)上本構(gòu)關(guān)系中：

在SPF/DB過(guò)程中，摩擦系數(shù)也是1個(gè)需要考慮的參數(shù)，摩擦系數(shù)是指2個(gè)表面間的摩擦力和作用在其一表面上的垂直力之比值。與表面的粗糙度有關(guān)，而與接觸面積的大小無(wú)關(guān)。本文采用COULOMB模型[13]，即反正切摩擦模型，定義為

式中：σt、ft為摩擦應(yīng)力；摩擦力 σn、fn為正應(yīng)力、正力；μ為摩擦系數(shù)；為相對(duì)滑動(dòng)速度方向矢量，vr為相對(duì)滑移速度；vrc為工件黏附開(kāi)始時(shí)的臨界相對(duì)速度，該值一般取vr的1%～10%。取值太大，會(huì)減小有效摩擦值；如取值太小，則模擬較難收斂，本文取值為2×105mm/s。

2 TC4合金3層板有限元數(shù)值模擬

TC4合金空心夾芯結(jié)構(gòu)的超塑性成形過(guò)程如圖1所示。該夾芯結(jié)構(gòu)由上、下面板加1層芯板構(gòu)成，成形時(shí)先將發(fā)生超塑性變形的部位涂上止焊劑，然后將3層板置于成形模具中，加熱至一定溫度，未涂止焊劑的部位發(fā)生擴(kuò)散連接（模擬計(jì)算時(shí)可在相應(yīng)位置施加約束來(lái)實(shí)現(xiàn)），其余部位面板和芯板在氣體壓力的作用下發(fā)生超塑性變形。

2.1 3層板有限元模型建立

模擬計(jì)算在有限元軟件MARC中進(jìn)行，首先對(duì)幾何模型進(jìn)行有限元網(wǎng)格劃分，如圖2所示。模型上、下面板厚度為2 mm，芯板厚度為1.3 mm，擴(kuò)散連接寬度為4 mm，上、下模具高度為8 mm。在超塑性過(guò)程中，上、下面板與芯板之間包括擴(kuò)散連接部位，擴(kuò)散帶寬度為4 mm；充分考慮到3層板材的對(duì)稱性，如圖3所示，有限元模型一共劃分為37600個(gè)單元，46512個(gè)節(jié)點(diǎn)。單元格類型選取SOLID單元（8節(jié)點(diǎn)）。

2.2 有限元參數(shù)設(shè)定

接觸條件：設(shè)定上、下模具為剛體，3層板為可變形體，把芯板劃分為未擴(kuò)散連接、上面板擴(kuò)散連接和下面板擴(kuò)散連接3部分，將其分別定義為變形體。設(shè)置接觸表面，定義上、下面板和模具相接觸，摩擦系數(shù)為0.2。上、下面板與芯板擴(kuò)散連接部分為黏著接觸。

邊界條件：如圖4所示，（對(duì)稱邊界）在板面對(duì)稱面上，沿X軸方向，Y=0；（固定邊界）邊緣節(jié)點(diǎn)上X，Y，Z方向位移為0；（面載荷）在未焊接部位選取單元面，施加面壓力。

載荷工況：設(shè)定最大氣壓為2.5 MPa；設(shè)定目標(biāo)應(yīng)變速率為1×10-3s-1，總工況時(shí)間為2000 s，步長(zhǎng)200步。

3 結(jié)果分析

3.1 成形過(guò)程

3層板結(jié)構(gòu)件在0～2000 s的超塑性變形過(guò)程的應(yīng)力如圖5所示。從圖中可見(jiàn)，TC4合金3層板超塑性成形過(guò)程分為自由膨脹階段、貼膜階段和充填階段。

在自由脹形過(guò)程中，隨著壓力的增大，模具對(duì)工件約束最少的地方將首先變形貼模。如圖5(a)所示，在t=0～130 s時(shí)間段，上面板發(fā)生超塑性變形，最先貼模；如圖5(b)所示，當(dāng)t=150 s時(shí)，下面板開(kāi)始碰觸剛性模具。由應(yīng)力云圖可知，在該過(guò)程中，芯板受力變形較大，出現(xiàn)頸縮現(xiàn)象；在t=1000 s時(shí)，工件貼模完成，進(jìn)入第3階段即充填階段。此時(shí)，結(jié)構(gòu)件與模具倒角區(qū)域出現(xiàn)較大的接觸應(yīng)力。

3.2 不同工藝參數(shù)對(duì)成形結(jié)果的影響

超塑性成形過(guò)程的影響因素包括工藝因素和設(shè)計(jì)因素。主要有摩擦系數(shù)、應(yīng)變速率敏感指數(shù)m、目標(biāo)應(yīng)變速率以及擴(kuò)散連接寬度和芯板厚度等。本文將著重從應(yīng)變速率敏感指數(shù)、目標(biāo)應(yīng)變速率、擴(kuò)散連接寬度角度，分析各因素對(duì)SPF過(guò)程的影響。

3.2.1 應(yīng)變速率敏感指數(shù)m

應(yīng)變敏感性指數(shù)對(duì)SPF超塑性成形過(guò)程有著極為重要的影響，其定義為塑性變形時(shí)材料的流變應(yīng)力對(duì)于應(yīng)變速率的敏感性參數(shù)，也就是隨著應(yīng)變速率增加板材所具有的硬化傾向系數(shù)[14]。對(duì)于普通金屬材料，m=0.02～0.2；而對(duì)于許多超塑性金屬材料，m=0.3～0.9。[15]

不同敏感指數(shù)下各面板壁厚分布如圖6所示，圖中橫軸為結(jié)構(gòu)件水平方向長(zhǎng)度；超塑性過(guò)程壓力時(shí)間分布如圖7所示；m=0.77時(shí)，結(jié)構(gòu)件成形后的不同視角如圖8所示。從圖6、7中可見(jiàn)，中間帶壁厚隨著m值的增大而增加，均勻程度有所降低。隨著m值的增大，擴(kuò)散連接處的減薄率顯著增大，頸縮現(xiàn)象嚴(yán)重，壁厚均勻程度明顯降低；與此同時(shí)，壓力達(dá)到2.5 MPa的時(shí)間隨著m的增大而增大；下面板中間帶，在m=0.77時(shí)，最大減薄率達(dá)到17.5%，面板均勻程度隨著m值的增大而降低；芯板隨m的增大，擴(kuò)散帶減薄率增加，當(dāng)m=0.77時(shí)，減薄率達(dá)到42%，而當(dāng)m=0.37時(shí)，芯板厚度較為均勻。但此時(shí)，構(gòu)件還并未完全貼模；當(dāng)m=0.57時(shí)，構(gòu)件貼模完成，且壁厚分布較為均勻。

從圖8中可見(jiàn)，在m=0.77時(shí)，構(gòu)件雖然已經(jīng)貼模完畢，但在Y軸方向上出現(xiàn)較為嚴(yán)重的縱向延伸。因此在進(jìn)行TC4合金超塑性成形模擬及試驗(yàn)時(shí)，選擇適當(dāng)?shù)膍值至關(guān)重要。在m=0.57時(shí)得到的模擬結(jié)果較為理想。

3.2.2 目標(biāo)應(yīng)變速率

在不同目標(biāo)應(yīng)變速率下各面板壁厚分布如圖9所示；超塑性過(guò)程壓力時(shí)間分布如圖10所示。在超塑性過(guò)程中，目標(biāo)應(yīng)變速率值一般選取為1×10-4～1×10-2s-1[16]。

從圖9中可見(jiàn)，當(dāng)目標(biāo)應(yīng)變速率為0.01和0.001時(shí)，TC4合金各層板的厚度分布情況基本一致。而在目標(biāo)應(yīng)變速率為0.0001時(shí)，發(fā)現(xiàn)上、下面板以及芯板厚度變化很小，模型在2000 s后幾乎沒(méi)有變形，也即在目標(biāo)應(yīng)變速率很小時(shí)，TC4合金的超塑性較差，無(wú)法滿足模擬過(guò)程的需要；反觀壓力時(shí)間曲線，從圖10中可見(jiàn)，目標(biāo)應(yīng)變速率為0.0001時(shí)的曲線并沒(méi)有達(dá)到設(shè)定的壓力最大值2.5 MPa。從而，在應(yīng)變速率為0.001時(shí)，板材超塑性性能達(dá)到峰值，也即在這一區(qū)間范圍內(nèi)，超塑性成形效果最優(yōu)。

3.2.3 擴(kuò)散連接寬度

為了有效地避免頸縮現(xiàn)象，應(yīng)當(dāng)選取適當(dāng)?shù)臄U(kuò)散連接寬度，這對(duì)設(shè)計(jì)工件尺寸提出了要求。不同擴(kuò)散連接寬度下各面板壁厚分布如圖11所示；超塑性過(guò)程的壓力時(shí)間分布如圖12所示。

從圖11中可見(jiàn)，上面板的中間連接段壁厚隨擴(kuò)散寬度的增大而加厚；隨著擴(kuò)散寬度增大，下面板最大減薄率增加，兩邊連接段的壁厚均勻程度隨擴(kuò)散寬度的增加而衰減，且中間帶兩側(cè)容易出現(xiàn)較大的波峰狀的壁厚分布帶；隨著芯板擴(kuò)散寬度的增大，縮頸現(xiàn)象逐漸突顯，最大減薄率分別為18.4%、26.9%、38.4%。此時(shí)，從圖12中可見(jiàn)，隨著擴(kuò)散連接寬度的增大，壓力分布趨勢(shì)不變，但是壓力達(dá)到2.5 MPa的時(shí)間則會(huì)相應(yīng)延長(zhǎng)。

4 結(jié)論

本文以SPF/DB組合工藝成形的TC4合金空心夾芯結(jié)構(gòu)為研究對(duì)象，采用MARC軟件對(duì)其超塑性成形過(guò)程進(jìn)行數(shù)值模擬，得到如下結(jié)論：

（1）TC4合金3層板超塑性成形過(guò)程分為自由膨脹階段、貼膜階段和充填階段。模擬發(fā)現(xiàn)芯板相對(duì)于上、下面板而言更容易發(fā)生縮頸現(xiàn)象，芯板與蒙皮擴(kuò)散連接位置為薄弱環(huán)節(jié)。

（2）3層板中間帶壁厚隨著m值的增大而增加，而且均勻程度有所降低；此外，擴(kuò)散連接處的減薄率顯著增大，縮頸現(xiàn)象嚴(yán)重，且壁厚均勻程度明顯降低；隨著m值的增大，SPF過(guò)程成形時(shí)間縮短，而材料延伸率普遍增加，從而導(dǎo)致材料在軸線方向上發(fā)生延伸變形。根據(jù)模擬結(jié)果可知，將應(yīng)變速率敏感指數(shù)控制在0.57左右時(shí)較為合適。

（3）應(yīng)變速率主要對(duì)板材的超塑性性能產(chǎn)生影響，應(yīng)當(dāng)控制應(yīng)變速率在10-3附近，此時(shí)超塑性成形效果最好。

（4）擴(kuò)散連接寬度主要對(duì)芯板壁厚分布有較大影響，即隨著擴(kuò)散連接寬度的增大，芯板縮頸現(xiàn)象愈加嚴(yán)重。

（5）在超塑性成形過(guò)程中，成形壓力控制的好壞決定了零件最終形狀與壁厚分布是否合理。合理的壓力-時(shí)間曲線可以保證得到壁厚分布均勻的構(gòu)件。利用有限元模擬技術(shù)可為得到最佳應(yīng)變速率下壓力-時(shí)間曲線，制定合理的氣壓控制參數(shù)提供必要的參考依據(jù)。

[1]許祥勝,趙振華,陳偉.外物損傷對(duì)TC4鈦合金的高周疲勞強(qiáng)度的影響[J].航空發(fā)動(dòng)機(jī),2017,43（3）:88-92.XU Xiangsheng,ZHAO Zhenhua,CHEN Wei.The influences of foreign damage on the high cycle fatigue behaviour of titanium alloy TC4[J].Aeroengine,2017,43（3）:88-92.(in Chinese)

[2]Bonet J,Bhargava P,Wood R D.Finite element analysis of the super plastic forming of thick sheet using the incremental flow formulation[J].International Journal for Numerical Methods in Engineering,1997,40（17）:3205-3228.

[3]施曉琦.鈦合金超塑成形/擴(kuò)散連接組合工藝研究[D].南京：南京航空航天大學(xué),2007.SHI Xiaoqi.Study on super plastic forming and diffusion bonding of titanium alloy[D].Nanjing：Nanjing University of Aeronautics Astronautics，2007.（in Chinese）

[4]趙毅,陶華.鈦合金板超塑性成形/擴(kuò)散連接成形因素影響的數(shù)值模擬[J].航空精密制造技術(shù),2007,43（1）:41-43.ZHAO Yi,TAO Hua.Numerical simulation of effect on SPF/DB of Ti alloy sheet[J].Aviation Precision Manufacturing Technology,2007,43（1）:41-43.（in Chinese）

[5]黃鋼華.鈦合金超塑性成形/擴(kuò)散連接的數(shù)值模擬及工藝研究[D].南京：南京航空航天大學(xué),2009.HUANG Ganghua.Numerical simulation and research on super plastic forming and diffusion bonding of titanium alloy[D].Nanjing:Nanjing University of Aeronautics Astronautics,2009.（in Chinese）

[6]崔元杰.TC4多層板結(jié)構(gòu)超塑成形/擴(kuò)散連接工藝數(shù)值模擬與試驗(yàn)研究[D].南京：南京航空航天大學(xué),2011.CUI Yuanjie.Numerical simulation and experimental study on super plastic forming and diffusion bonding process of TC4 multilayer structure[D].Nanjing：Nanjing University of Aeronautics Astronautics,2011.（in Chinese）

[7]JIANG Shaosong,LU Zhen,HE Xiaodong.Super plastic forming of Ti6Al4V titanium alloy cylinder with near uniform thickness distribution[J].Chinese Journal of Nonferrous Metals,2012,22（s2）:472-478.

[8]Chen Y,Kibble K,Hall R,et al.Numerical analysis of superplastic blow forming of Ti-6Al-4V alloys[J].Materials&Design,2001,22（8）:679-685.

[9]Chen D C.Finite element analysis of super-plastic blow-forming of Ti6Al4V sheet into closed ellip-cylindrical die[J].International Journal of Simulation Modelling,2010,9（1）:17-27.

[10]向毅斌,吳詩(shī)惇.超塑性脹形成形時(shí)間影響因素的數(shù)值分析[J].模具技術(shù),2000（2）:9-14.XIANG Yibin,WU Shichun.Numerical analysis of influencing factors on super plastic bulging forming time[J].Die Technology,2000（2）:9-14.（in Chinese）

[11]孫青平.空心風(fēng)扇葉片桁架式夾芯結(jié)構(gòu)件疲勞特性研究[D].南京：南京航空航天大學(xué),2014.SUN Qingping.Research on fatigue characteristics of hollow fan blade truss structure [D].Nanjing：Nanjing University of Aeronautics Astronautics,2014.（in Chinese）

[12]蔡云.TC4鈦合金高溫拉伸力學(xué)性能研究和組織演變[D].南京：南京航空航天大學(xué),2009.CAI Yun.Mechanical properties and microstructure evolution of TC4 titanium alloy under high temperature tensile test[D].Nanjing：Nanjing University of Aeronautics Astronautics,2009.（in Chinese）

[13]李靖誼,王衛(wèi)英,陳明和.控制超塑性充模脹形過(guò)程的計(jì)算機(jī)仿真技術(shù)[J].數(shù)據(jù)采集與處理,1998，13（4）：358-361.LI Jingyi,WANG Weiying,CHEN Minhe.Computer simulation of super plastic bulging process[J].Data Acquisition and Processing,1998，13（4）：358-361.（in Chinese）

[14]Vanderhasten M,Rabet L,Verlinden B.Ti-6Al-4V:Deformation map and modelisation of tensile behaviour[J].Materials&Design,2008,29（6）:1090-1098.

[15]Kuperfarb G,Germain Y,Abouaf M.A mechanical study of Ti-6Al-4V sheet[J].Journal of Mechanical Working Technology,1987,14（2）:159-172.

[16]Giuliano G.Constitutive equation for super-plastic Ti-6Al-4V alloy[J].Materials&Design，2008,29（7）:1330-1333.