鐘軼寧，王惠敏，劉德貴，周文龍，付雪松，李志強(qiáng)

（1. 大連理工大學(xué) 材料科學(xué)與工程學(xué)院，遼寧 大連 116024；2. 北京航空制造工程研究所，北京 100024）

鈦合金具有比強(qiáng)度高、耐腐蝕和耐高溫性能好的特點(diǎn)，是航空飛行器及其發(fā)動(dòng)機(jī)的主要結(jié)構(gòu)材料之一[1—2]，但是鈦合金構(gòu)件在服役過(guò)程中容易發(fā)生疲勞斷裂，這是目前存在的主要問(wèn)題之一[3]。表面強(qiáng)化對(duì)于提高材料的抗疲勞能力效果顯著，在已有表面強(qiáng)化技術(shù)中，表面噴丸強(qiáng)化處理具有效果顯著、價(jià)格低廉、設(shè)備簡(jiǎn)單的優(yōu)點(diǎn)而備受關(guān)注[4—5]。

噴丸強(qiáng)化的主要工藝參數(shù)（彈丸尺寸、噴丸速度、噴丸覆蓋率等）之間存在復(fù)雜的交互作用，工藝參數(shù)設(shè)計(jì)不合理，不但達(dá)不到強(qiáng)化效果，還浪費(fèi)財(cái)力、物力，甚至破壞噴丸效果[6]。目前，我國(guó)對(duì)于鈦合金零件噴丸強(qiáng)化工藝參數(shù)（如覆蓋率等）對(duì)強(qiáng)化效果影響規(guī)律的研究尚不夠系統(tǒng)[7]，使噴丸強(qiáng)化的工程質(zhì)量難以保證，因此需要進(jìn)行大量的實(shí)驗(yàn)研究來(lái)優(yōu)化工藝參數(shù)，但是將耗費(fèi)大量的時(shí)間、人力和實(shí)驗(yàn)費(fèi)用，而采用數(shù)值模擬方法輔助實(shí)驗(yàn)研究，可以有效減少實(shí)驗(yàn)量，極大降低成本[8]。數(shù)值模擬噴丸過(guò)程的核心是如何建立模型。Meguid[9—10]最早提出了單彈丸與多彈丸錯(cuò)排的兩種模型，為噴丸強(qiáng)化的模擬研究奠定了基礎(chǔ)。H. Y. Miao[11]將Matlab與ANSYS結(jié)合，建立了彈丸在空間隨機(jī)分布的模型。洪滔[12]利用離散元研究大量彈丸間相互作用后的能量，Kovthaman Murugaratnam[13]將有限元與離散元結(jié)合，研究了大量彈丸的噴射過(guò)程，流體力學(xué)[14]與 SPH[15]方法等也被引入了到噴丸的模擬研究中，數(shù)值模擬在噴丸研究中發(fā)揮了越來(lái)越重要的作用。

近年來(lái)，陶瓷彈丸的濕噴丸技術(shù)因能減少摩擦、保護(hù)材料表面，以及其彈丸可回收、無(wú)粉塵污染[16]、彈丸硬度高、破碎率低[17]等特點(diǎn)，得到了關(guān)注。有關(guān)干/濕噴丸強(qiáng)化效果的對(duì)比[18]、濕噴丸殘余應(yīng)力[19]、濕噴丸后的組織和性能[16]等已有較多研究文獻(xiàn)。由于濕噴丸過(guò)程中涉及流固耦合問(wèn)題，并且接觸體間的相互作用較為復(fù)雜，目前針對(duì)濕噴丸強(qiáng)化過(guò)程數(shù)值模擬的研究報(bào)道較為罕見。文中利用ABAQUS軟件中的CEL方法，建立了一種濕噴丸的彈丸隨機(jī)分布模型，開展覆蓋率對(duì) TC4板材的殘余應(yīng)力以及表面粗糙度的模擬分析，并與實(shí)驗(yàn)結(jié)果進(jìn)行了對(duì)比。

1 濕噴丸模型的建立

利用ABAQUS軟件中的CEL方法，建立了TC4鈦合金板濕噴丸過(guò)程的模型。此方法兼具歐拉網(wǎng)格與拉格朗日網(wǎng)格的優(yōu)點(diǎn)，能簡(jiǎn)單而有效地實(shí)現(xiàn)流固耦合。所采用的磨液比（水和陶瓷彈丸的質(zhì)量比）為20%，陶瓷彈丸半徑R=0.15 mm，根據(jù)水和陶瓷彈丸的密度，得出特定區(qū)域0.8 mm×0.8 mm×4 mm內(nèi)陶瓷彈丸的數(shù)量，并利用rand隨機(jī)函數(shù)（見式(1)），確定空間中隨機(jī)分布彈丸的中心坐標(biāo)。式(1)保證了所有的陶瓷彈丸都完全處于參考體內(nèi)，去除彈丸在空間中發(fā)生重疊的坐標(biāo)點(diǎn)。為了避免陶瓷彈丸在碰撞板材之前即脫離水域，且考慮到水對(duì)板材產(chǎn)生的撞擊作用有限，將參考體的尺寸擴(kuò)大，設(shè)為1.6 mm×1.6 mm×4 mm（見圖1）。TC4板材的密度為4.51 g/cm3，彈性模量為121.4 GPa，泊松比為0.3，陶瓷彈丸的密度為3.8 g/cm3，彈性模量為350 GPa，泊松比為0.26。由于陶瓷彈丸的硬度較高，且其變形問(wèn)題并不是文中研究的重點(diǎn)，將其設(shè)為剛體，TC4板材為變形體。

圖1 彈丸隨機(jī)分布模型Fig.1 Random arrangement model

考慮到板材在濕噴丸過(guò)程中發(fā)生了塑性變形，需要采用真實(shí)應(yīng)力與應(yīng)變對(duì)其塑性進(jìn)行定義，真實(shí)應(yīng)力-應(yīng)變與名義應(yīng)力-應(yīng)變間的轉(zhuǎn)換公式見式(2)和式(3)。由式(2)和式(3)、拉伸實(shí)驗(yàn)及彈性模量，可得TC4的塑性參數(shù)，真實(shí)應(yīng)力分別為 973, 1014, 1047, 1108,1156, 1198, 1234, 1264, 1281 MPa的塑性應(yīng)變分別為0, 0.012, 0.022, 0.042, 0.061, 0.081, 0.101, 0.121,0.130。

式中：σ為真實(shí)應(yīng)力；ε為真實(shí)應(yīng)變；σnom為名義應(yīng)力；εnom為名義應(yīng)變。

材料的體積響應(yīng)是由狀態(tài)方程控制的，假設(shè)水是不可壓縮粘性流體，采用ABAQUS/Explicit中的線性Us-Up狀態(tài)方程（Us為沖擊波波速，Up為波后粒子速度，Us=c0+sUp），描述流場(chǎng)從而實(shí)現(xiàn)水的建模。線性Us-Up形式的 Mie-Grüneisen狀態(tài)方程見式(4)[20]。水的密度為1000 kg/m3，粘度為0.001 Pa……s，水中的聲速為1450 m/s，斜率系數(shù)和材料參數(shù)均為0。

式中：p為水的壓力；ρ0為水的初始密度；c0和s為Us-Up關(guān)系的系數(shù)；c0為水中的聲速；s為斜率系數(shù)；Г0為材料參數(shù)；Em為比內(nèi)能；η為名義體積壓縮應(yīng)變，其表達(dá)式為η=1-ρ0?ρ，ρ為擾動(dòng)后水的密度。

分別采用C3D8R, C3D4, EC3D8單元對(duì)TC4板材、陶瓷彈丸、歐拉域進(jìn)行網(wǎng)格劃分。為了解決減縮積分單元的沙漏問(wèn)題，對(duì)板材以及歐拉域的中心區(qū)域進(jìn)行網(wǎng)格細(xì)化。利用面接觸定義陶瓷彈丸與板材間的接觸，摩擦因數(shù)為 0.2，利用通用接觸定義水與其他部件間的接觸，忽略摩擦作用。水和陶瓷彈丸的初始速度均為80 m/s。通過(guò)對(duì)預(yù)定場(chǎng)中初始狀態(tài)的定義，實(shí)現(xiàn)計(jì)算結(jié)果的導(dǎo)入，以達(dá)到高覆蓋率。

2 數(shù)值模擬與分析

2.1 覆蓋率對(duì)殘余應(yīng)力的影響

不同覆蓋率下板材表面的殘余應(yīng)力分布見圖2。當(dāng)覆蓋率由50%增加到100%時(shí)，板材表面節(jié)點(diǎn)的最大壓應(yīng)力由968 MPa增加至1148 MPa。當(dāng)覆蓋率繼續(xù)增加至 200%時(shí)，表面的最大殘余壓應(yīng)力減小，變?yōu)?047 MPa。

由于彈丸位置呈隨機(jī)狀態(tài)、且變形區(qū)之間會(huì)相互影響，所以板材表面各點(diǎn)的殘余應(yīng)力沿厚度方向的分布狀態(tài)以及隨覆蓋率的變化規(guī)律可能會(huì)不同。圖3a和圖3b分別為板材上表面的(0, 0)點(diǎn)以及(0.08,-0.08)點(diǎn)在覆蓋率為 200%時(shí)殘余應(yīng)力的分布。點(diǎn)(0.08, 0.08)厚度方向的最大殘余壓應(yīng)力出現(xiàn)在次表層，并且明顯大于點(diǎn)(0, 0)中厚度方向的最大殘余壓應(yīng)力，而點(diǎn)(0, 0)厚度方向的最大殘余壓應(yīng)力出現(xiàn)在表層。兩點(diǎn)在不同覆蓋率時(shí)，厚度方向殘余應(yīng)力分布曲線見圖4?？梢姡?dāng)取點(diǎn)位置不同時(shí)，殘余應(yīng)力隨著覆蓋率的變化規(guī)律存在差異。為了避免取點(diǎn)不同的影響，文中對(duì)同一厚度方向各節(jié)點(diǎn)的殘余應(yīng)力進(jìn)行了取平均值處理。

對(duì)同一厚度內(nèi)所有節(jié)點(diǎn)取平均值處理后的板材厚度方向的殘余應(yīng)力分布曲線見圖5?？梢姡?dāng)覆蓋率由50%增加至100%時(shí)，厚度方向的最大殘余壓應(yīng)力增大，由737 MPa增至850 MPa，位置由距離上表面20 μm左右處延伸至距離40 μm處。當(dāng)覆蓋率繼續(xù)增大到200%時(shí)，最大殘余壓應(yīng)力為829 MPa，所處深度并無(wú)明顯增加。隨著濕噴丸覆蓋率的增加，板材厚度方向的殘余壓應(yīng)力層的深度在 80～100 μm 的范圍內(nèi)略有變化。此外，3種覆蓋率下表面殘余壓應(yīng)力未見明顯差異。

圖2 不同覆蓋率下板材表面殘余應(yīng)力分布Fig.2 Residual stress contours on plate surface under different coverage

圖3 板材表面兩點(diǎn)厚度方向應(yīng)力Fig.3 Residual stress of two different points on plate surface

圖4 兩點(diǎn)厚度方向殘余應(yīng)力變化曲線Fig.4 Residual stress along the thickness direction of two different points on plate surface

圖5 不同覆蓋率下板材厚度方向殘余應(yīng)力分布曲線Fig.5 Residual stress along the plate thickness under different coverage

2.2 覆蓋率對(duì)粗糙度的影響

表面粗糙度的評(píng)價(jià)指標(biāo)主要有兩種，分別為輪廓的算術(shù)偏差和輪廓的最大高度。文中采用輪廓的最大高度，即PV粗糙度來(lái)對(duì)表面粗糙度進(jìn)行評(píng)價(jià)。不同覆蓋率時(shí)，板材表面各節(jié)點(diǎn)沿z軸方向的位移見圖6?？梢?，當(dāng)濕噴丸覆蓋率由50%增加至100%時(shí)，變形區(qū)域擴(kuò)大，波峰高度增加明顯，由4.24 μm增加至5.87 μm。當(dāng)覆蓋率由 100%增加至 200%時(shí)，板材表面波峰高度和波谷深度均發(fā)生明顯增加，這與覆蓋率增加至 200%時(shí)，在采用的模型中，板材表面的節(jié)點(diǎn)受到了與之前相同的撞擊作用有關(guān)。PV粗糙度變化曲線見圖7，當(dāng)覆蓋率由50%增加至200%時(shí)，表面的PV粗糙度分別為10.93, 12.29, 16.32 μm，逐漸增加。

圖6 不同覆蓋率下板材表面沿z軸方向位移Fig.6 z displacement on plate surface under different coverage

3 實(shí)驗(yàn)結(jié)果與模擬結(jié)果的對(duì)比

為了驗(yàn)證模型的可靠性，利用JY-120WB型液體噴丸機(jī)進(jìn)行相關(guān)工藝參數(shù)的濕噴丸強(qiáng)化實(shí)驗(yàn)。利用EMPYREANX型X射線衍射儀，測(cè)定濕噴丸后板材厚度方向殘余應(yīng)力。使用NewView5022型ZYGO表面輪廓儀測(cè)定表面粗糙度，其垂直分辨率可達(dá) 0.1 nm。

不同覆蓋率下板材厚度方向殘余應(yīng)力的實(shí)驗(yàn)結(jié)果見圖8。在實(shí)驗(yàn)結(jié)果中，板材表面殘余應(yīng)力并未隨覆蓋率的增加而顯著變化，應(yīng)力值在700 MPa上下浮動(dòng)。覆蓋率超過(guò) 200%后，最大殘余壓應(yīng)力及其所處深度的變化亦不明顯。隨著覆蓋率的增加，殘余壓應(yīng)力層的厚度均有少許的增加。由此可見，實(shí)驗(yàn)和模擬的研究結(jié)果所反應(yīng)的規(guī)律較為一致。

圖7 不同覆蓋率下板材表面PV粗糙度Fig.7 PV roughness of plate surface under different coverage

圖8 噴丸覆蓋率對(duì)殘余應(yīng)力場(chǎng)影響Fig.8 Effect of wet peening coverage on residual stress

圖9 濕噴丸覆蓋率對(duì)粗糙度的影響Fig.9 Effect of wet shot peening coverage on roughness

不同覆蓋率下，受噴材料表面PV粗糙度的實(shí)驗(yàn)結(jié)果見圖9。比較實(shí)驗(yàn)和模擬的結(jié)果可以發(fā)現(xiàn)，雖然所得結(jié)果的具體數(shù)字存在差距，但得到的規(guī)律比較一致，即覆蓋率在一定范圍內(nèi)時(shí)，受噴材料表面粗糙度隨著覆蓋率增加而增大。在模擬分析結(jié)果中，覆蓋率的增加使粗糙度增大的速率較快，主要是因?yàn)閿?shù)值模擬模型中將覆蓋率設(shè)置為丸粒重復(fù)加載率，相當(dāng)于各節(jié)點(diǎn)上受到的再次撞擊率，導(dǎo)致受噴材料表面的峰高和谷深的變化更顯著，從而使PV粗糙度明顯增大。實(shí)際實(shí)驗(yàn)中丸粒撞擊表面所呈現(xiàn)的隨機(jī)狀態(tài)，其重復(fù)加載率比較低，存在較多的丸?；プ埠蛷椏訐崞降那闆r。為此，數(shù)學(xué)模型在這方面有待進(jìn)一步修正改進(jìn)。

4 結(jié)論

1) 利用ABAQUS中的CEL方法建立了TC4鈦合金板材的濕噴丸強(qiáng)化模型，通過(guò)函數(shù)方法獲得了陶瓷彈丸在空間隨機(jī)分布的位置坐標(biāo)，并利用導(dǎo)入預(yù)定義場(chǎng)的方式實(shí)現(xiàn)了模型的多次加載，獲得了不同覆蓋率下的濕噴丸強(qiáng)化結(jié)果。

2) 不同覆蓋率下的實(shí)驗(yàn)與模擬結(jié)果對(duì)比證實(shí)，濕噴丸過(guò)程模型可以較好地模擬濕噴丸強(qiáng)化后板材表層殘余壓應(yīng)力場(chǎng)和表面粗糙度的變化規(guī)律。

3) 隨著濕噴丸覆蓋率的增加，TC4板材表面殘余壓應(yīng)力值不隨覆蓋率發(fā)生明顯變化，殘余壓應(yīng)力層厚度略有增加，表面的PV粗糙度逐漸增加。