潘為民，程東華，候森，張玉卓

(河南科技大學(xué)機(jī)電工程學(xué)院，河南洛陽 471003)

0 前言

熱處理工藝在零件的極限性能和極限服役性能中占有重要地位。對零件進(jìn)行熱處理，可使零件材料的內(nèi)部結(jié)構(gòu)得以改善，并增強(qiáng)零件的實(shí)用性。但在進(jìn)行熱處理工藝時，會或多或少地產(chǎn)生變形等質(zhì)量問題，進(jìn)而影響零件的使用性能[1]。大型空腔葉片在固溶熱處理過程中，由于葉片的空腔區(qū)域厚度較小且曲率變化相對較大，葉片實(shí)體部位與空腔連接處在加熱和冷卻過程中極易受熱不均，故在空腔區(qū)域和空腔與實(shí)體相連部位容易產(chǎn)生變形。

國外熱處理試件的數(shù)值模擬研究最早出現(xiàn)在20世紀(jì)70年代，瑞典科學(xué)家HILDENWALL、ERLCSSON[2]模擬驗(yàn)證了碳鋼的淬火應(yīng)力；FERNANDES對合金珠光體展開了相似的研究，以確認(rèn)離散性的可靠性。時至二十世紀(jì)后期，國外眾多學(xué)者已經(jīng)為非穩(wěn)態(tài)溫度場編制了許多計算機(jī)仿真模擬程序。RAMERSTORFER研究了熱處理工藝過程中的熱彈塑性，比較了熱應(yīng)力和相變、蠕變等對模擬仿真結(jié)果的影響[3-5]。自1990年以來，隨著計算機(jī)技術(shù)不斷創(chuàng)新發(fā)展，溫度場建模中最復(fù)雜的非線性問題逐漸得到解決，CIALKOWSI等研究了使用Trefez方法解決二維穩(wěn)態(tài)熱導(dǎo)傳反應(yīng)的逆向求解方法[6]。

本文作者基于大型空腔葉片在固溶熱處理工藝過程中的變形問題，通過對大型空腔葉片固溶熱處理過程進(jìn)行仿真模擬，使該過程“可視化”，得到大型空腔葉片鑄件冷卻后各部位的位移場分布情況，利用位移場的分布可以知道鑄件各方向收縮量的大小及鑄件的形狀畸變情況，分析并預(yù)測溫度場、應(yīng)力場對大型空腔葉片變形的影響機(jī)制，為探索熱處理工藝參數(shù)對大型空腔葉片成型的影響、優(yōu)化熱處理工藝、控制合理的加工余量提供了依據(jù)。

1 葉片結(jié)構(gòu)特征及熱處理工藝

1.1 葉片結(jié)構(gòu)特征

大型空腔葉片形似翼型，葉身中部是不規(guī)則的空腔結(jié)構(gòu)，且曲率較大，對流體的流動特性影響較大。為了減輕大型空腔葉片質(zhì)量，內(nèi)部型腔的曲率隨葉片外形曲率的改變而改變。由于葉片的型面復(fù)雜，且空腔曲面的曲率變化較大，又具有很大的體積和質(zhì)量，故在制造過程中，通常在葉片表面設(shè)計工藝預(yù)留孔和輔助立柱來減少葉片應(yīng)力集中和熱處理變形等現(xiàn)象。為了便于葉片在加工過程中的移動，在葉片上焊接有吊環(huán)和輔助墊塊，如圖1所示。

圖1 大型空腔葉片圖片

1.2 葉片熱處理工藝

該大型空腔葉片材料采用我國自主研制的607A鋼。607A不銹鋼在冶煉時適當(dāng)?shù)丶尤肓艘恍?qiáng)化元素如Mo、Al、Ti等，也可以把Ti改為Nb、Ta等元素形成馬氏體不銹鋼。607A鋼經(jīng)過固熔+時效處理后，葉片既能保持高強(qiáng)度，又具有高韌性。607A不銹鋼主要化學(xué)成分見表1。冶煉607A不銹鋼時為了提高冶煉成功率，降低了不銹鋼中硫磷的含量，把原規(guī)定的硫含量不大于0.030%改為硫含量不大于0.020%，通過對607A不銹鋼冶煉工藝的不斷改善，現(xiàn)在607A鋼硫含量甚至可以達(dá)到不大于0.015%，這大大提高了冶煉607A鋼的成功率，不僅降低了經(jīng)濟(jì)損失，還提高了葉片的使用壽命。

表1 607A不銹鋼成分(質(zhì)量分?jǐn)?shù)) 單位：%

為生產(chǎn)出滿足高性能要求的葉片鑄件，通過實(shí)驗(yàn)分析大型空腔葉片的生產(chǎn)工藝過程，對607A鋼采用固溶處理的熱處理工藝。607A鋼在熱處理爐加熱至1 100 ℃，在1 100 ℃下保溫240 min然后進(jìn)行空冷。為了改善機(jī)加工的性能，一般在固溶熱處理后進(jìn)行時效處理，以達(dá)到最大程度的軟化，為了得到最佳的強(qiáng)化效果，需要進(jìn)行二次時效[7]。鑄件時效：鑄件在(720±10)℃持續(xù)加熱8 h，再降至(620±10)℃保溫8 h后空冷。

2 葉片熱處理過程數(shù)值模擬

2.1 有限元模型的建立

采用有限元軟件SolidWorks建立大型空腔葉片有限元模型，另存為STL格式導(dǎo)入DEFORM軟件中，選用四面體單元類型[8]，劃分網(wǎng)格數(shù)80 000個，節(jié)點(diǎn)數(shù)5 902個，元素19 974，多面體42 820個，端點(diǎn)21 412個。應(yīng)變場采用直接耦合法，大型空腔葉片材料選擇彈塑性材料，圖2所示為大型空腔葉片三維模型。計算中，在葉片進(jìn)水邊施加x、y、z三個方向的約束將其固定。圖3所示為葉片網(wǎng)格劃分模型。

圖2 大型空腔葉片三維模型

2.2 初始條件的設(shè)定

文中熱處理模擬工藝采用固溶熱處理，607A鋼在熱處理爐加熱至1 100 ℃，在1 100 ℃下保溫240 min然后進(jìn)行空冷。工件表面初始溫度為 20 ℃，葉片加熱時與外界存在輻射，降溫過程中，與空氣之間存在輻射和自然對流。

對流可表示為

q=hf(TS-T0)

(1)

其中：hf為工件表面的對流換熱系數(shù)(W/(m2·℃))；TS為工件表面溫度(℃)；T0為周圍流體溫度(℃)。

工件在加熱過程中，材料的輻射換熱系數(shù)和對流換熱系數(shù)均隨溫度改變而改變。在進(jìn)行模擬時需要分析輻射和對流的影響，并設(shè)置換熱系數(shù)為H[9]。因此，大型空腔葉片外表面與外界之間的換熱系數(shù)可表示為

q=Hf(TS-T0)

(2)

熱處理過程中葉片表面溫度變化較大，內(nèi)部會因?yàn)榻M織應(yīng)力和熱應(yīng)力的影響產(chǎn)生較大的熱應(yīng)力。當(dāng)熱應(yīng)力較大時，材料易產(chǎn)生彈性變形，熱應(yīng)力達(dá)到一定程度時會發(fā)生塑性變形。此時，需要利用平衡微分方程，分析應(yīng)力與應(yīng)變及應(yīng)變與位移的關(guān)系[10]，模擬時hf設(shè)置為0.1，Hf設(shè)置為0.02。

2.3 溫度場數(shù)值模擬及分析

模擬大型空腔葉片固溶處理過程，圖4所示為葉片表面固溶處理升溫過程中的溫度模擬圖。不難看出：加熱初期，葉片厚度較大的進(jìn)水邊升溫慢，厚度較小的曲面空腔升溫快；隨著爐內(nèi)溫度不斷升高，葉片整體溫度基本上都達(dá)到1 100 ℃。葉片出爐空冷時，環(huán)境溫度為室溫30 ℃，因此葉片表面迅速降溫，厚度較小的曲面空腔冷卻速度較快，進(jìn)水邊冷卻速度較慢。隨著冷卻時間的延長，整個葉片的溫度均快速下降，且空腔和實(shí)體部位的溫度差逐漸縮小，溫度分布也逐漸趨于平緩，最終葉片整體達(dá)到室溫30 ℃左右，如圖5所示。

圖4 升溫至1 100 ℃

圖5 冷卻至室溫

2.4 應(yīng)力場數(shù)值模擬及分析

大型空腔葉片試件模擬結(jié)果如圖6—7所示，大型空腔葉片經(jīng)固溶處理冷卻后，等效應(yīng)力最大值主要集中在葉片上冠的尖角A、B、C、D四個區(qū)域，最大值約為152 MPa；較大值主要分布在空腔與實(shí)體連接部位的進(jìn)水邊和曲面曲率變化較大的部位，值在77.1～115 MPa之間；最小值主要集中在曲率變化較小的葉片曲面中心，應(yīng)力值為1.91～39.5 MPa。應(yīng)力較大的區(qū)域容易發(fā)生變形，因此，葉片上冠的尖角區(qū)域和空腔與實(shí)體連接的部位容易發(fā)生變形。

圖6 等效應(yīng)力分布前視圖

圖7 等效應(yīng)力分布后視圖

2.5 應(yīng)變場數(shù)值模擬及分析

由葉片固溶熱處理的變形分布圖可知：熱處理工藝在葉片鑄造成型時影響較大，變形較大的區(qū)域主要在葉身曲率較大的部位以及葉片空腔與實(shí)體連接的部位。這表明大型空腔葉片在進(jìn)行熱處理工藝時產(chǎn)生較大形變，主要是由于葉片內(nèi)腔壁厚較小且曲率變化相對較大，而葉片進(jìn)水端和上冠又是實(shí)體，葉片實(shí)體部位與空腔連接處在加熱和冷卻過程中極易受熱不均，故在空腔區(qū)域和空腔與實(shí)體相連部位容易產(chǎn)生變形。而且葉片空腔部位和實(shí)體部位存在不同的散熱率，也會導(dǎo)致葉片在不同區(qū)域發(fā)生不同的應(yīng)力應(yīng)變，這也是葉片在此易產(chǎn)生較大變形的原因之一。大型空腔葉片經(jīng)固溶熱處理后，其中葉片上冠的尖角A、B、C、D四個區(qū)域在Y軸正方向變形較大，值約為5.57 mm；葉片的進(jìn)水端E、F、G、H四個區(qū)域易在Y軸負(fù)方向變形較大，值約為6.05 mm。圖8為葉片經(jīng)固溶熱處理后Y軸方向變形分布的前視圖，圖9為葉片經(jīng)固溶熱處理后Y軸方向變形分布的后視圖。

圖8 葉片Y軸方向變形前視圖

圖9 葉片Y軸方向變形后視圖

3 實(shí)驗(yàn)驗(yàn)證

通過測量大型空腔葉片熱處理前后的點(diǎn)云數(shù)據(jù)，根據(jù)測得的點(diǎn)云數(shù)據(jù)重構(gòu)葉片的三維模型，對比葉片在熱處理前和熱處理后的三維模型，進(jìn)而得到熱處理前后葉片的變形情況。

根據(jù)大型空腔葉片的外形特征，葉片表面曲率驟變處變化較大難以測量，并且掃描儀測不到激光照射不到的區(qū)域，故采用多個基準(zhǔn)點(diǎn)從不同角度轉(zhuǎn)換測量大型空腔葉片，獲取葉片熱處理前和熱處理后的點(diǎn)云數(shù)據(jù)。三維激光掃描測量示意如圖10所示。

圖10 三維激光掃描測量

測量之前，利用墊塊調(diào)整葉片的測量角度，并確定葉片的放置姿態(tài)和測量方法。由于大型空腔葉片體積較大且葉身表面曲率變化大又存在測量死角，故利用兩側(cè)轉(zhuǎn)站測量的方法[11]對大型空腔葉片進(jìn)行測量。轉(zhuǎn)站測量時，為了提高葉片的測量精度，根據(jù)被測葉片的形狀特征取5個固定點(diǎn)為基準(zhǔn)點(diǎn)。由于測量獲得的原始數(shù)據(jù)數(shù)量極大，存在許多噪聲點(diǎn)和冗余點(diǎn)，這些點(diǎn)將不利于后續(xù)葉片三維模型的重建，因此利用K-鄰域搜索算法(如圖11所示)對這些噪聲點(diǎn)和冗余點(diǎn)進(jìn)行去除[12]，從而有效精簡大型空腔葉片熱處理前后的的點(diǎn)云數(shù)據(jù)。

圖11 K-鄰域搜索算法

利用精簡后的點(diǎn)云數(shù)據(jù)得到擬合曲線[13]，再由曲面生成曲面片，最后經(jīng)過拼接得到完整的大型空腔葉片的曲面模型。利用最小二乘法[14]同時使用初始匹配和精確匹配算法[15]分析重構(gòu)的葉片模型精度，可以得到重構(gòu)模型各部位的精確程度及變形趨勢。通過對比數(shù)據(jù)發(fā)現(xiàn)，數(shù)值模擬變形最大值為5.57 mm，最小值為-6.05 mm，在進(jìn)行固溶熱處理時，葉片的空腔區(qū)域剛性不足，空腔壁厚較小且曲率變化相對較大，熱處理后殘余的應(yīng)力和應(yīng)變無法完全去除，導(dǎo)致葉片的尺寸發(fā)生改變，有85.41%的點(diǎn)集中分布在空腔區(qū)域且變形量在-2.66～2.29 mm之間。葉片自身的表面粗糙度也會在一定程度上影響重構(gòu)模型的評定精度[16]。綜合考慮，重構(gòu)的葉片三維模型可以達(dá)到該葉片的設(shè)計精度要求，即可以用重構(gòu)后的葉片三維模型來計算分析該葉片熱處理前后的變形程度。

通過測量大型空腔葉片熱處理前和熱處理后的點(diǎn)云數(shù)據(jù)，根據(jù)測得的點(diǎn)云數(shù)據(jù)重構(gòu)葉片的三維模型，對比葉片在熱處理前和熱處理后的三維模型，進(jìn)而得到熱處理前后葉片的變形情況。大型空腔葉片熱處理變形的測量結(jié)果如圖12—13所示?？梢钥闯觯鹤冃屋^大的點(diǎn)主要集中在葉片上冠尖角區(qū)域及空腔與實(shí)體連接的部位，最大值為4.76 mm，最小值為-5.13 mm，有85.41%的點(diǎn)變形程度在-2.66～2.29 mm之間，這些點(diǎn)主要集中在曲率變化較小的空腔部位。

圖12 熱處理變形評定結(jié)果前視圖

圖13 熱處理變形評定結(jié)果后視圖

通過分析大型空腔葉片固溶熱處理加工后呈現(xiàn)出的變形趨勢，發(fā)現(xiàn)該葉片經(jīng)固溶熱處理后變形較大，主要是由于葉片內(nèi)腔壁厚較小且曲率變化相對較大，而葉片進(jìn)水端和上冠又是實(shí)體，葉片實(shí)體部位與空腔連接處在加熱和冷卻過程中極易受熱不均，故在空腔區(qū)域和空腔與實(shí)體相連部位容易產(chǎn)生變形。此外，在加熱和冷卻過程中，葉片受熱不均也會使不同區(qū)域產(chǎn)生不同的應(yīng)力應(yīng)變，這也是葉片產(chǎn)生應(yīng)力和應(yīng)變的主要原因，熱處理后殘余的應(yīng)力和應(yīng)變無法完全去除，必定導(dǎo)致葉片的尺寸變化。故大型空腔葉片在經(jīng)過固溶熱處理后會在葉片尖角區(qū)域和空腔與實(shí)體相連的部位產(chǎn)生較大的變形量。

4 結(jié)論

(1)通過對大型空腔葉片固溶熱處理工藝進(jìn)行數(shù)值模擬，獲得了固溶熱處理后的應(yīng)力、應(yīng)變場。與用三維激光掃描儀測量的葉片數(shù)據(jù)進(jìn)行了對比，模擬與實(shí)驗(yàn)吻合度較好，表明了葉片數(shù)值模擬的有效性。

(2)大型空腔葉片經(jīng)固溶熱處理后，變形較大的區(qū)域主要在葉片上冠的尖角區(qū)域以及葉片空腔與實(shí)體連接的部位。數(shù)值模擬變形最大值為5.57 mm，最小值-6.05 mm，實(shí)驗(yàn)測量最大值為4.76 mm，最小值-5.13 mm，變形程度基本一致，變形分布趨勢一致。

(3)在進(jìn)行固溶熱處理時，葉片的空腔區(qū)域剛性不足，空腔壁厚較小且曲率變化相對較大，熱處理后殘余的應(yīng)力和應(yīng)變無法完全去除，導(dǎo)致葉片的尺寸發(fā)生改變，有85.41%的點(diǎn)集中分布在空腔區(qū)域且變形量在-2.66～2.29 mm之間。