范立維，梁忠偉，劉曉初，吳俊，吳子軒，耿晨，謝鑫成

表面強(qiáng)化技術(shù)

基于正態(tài)分布的GCr15軸承鋼強(qiáng)化研磨殘余應(yīng)力場(chǎng)數(shù)值模擬

范立維a,b,c，梁忠偉a,b,c，劉曉初a,b,c，吳俊a,b,c，吳子軒a,b,c，耿晨a,b,c，謝鑫成a,b,c

（廣州大學(xué) a.機(jī)械與電氣工程學(xué)院 b.廣州市金屬材料強(qiáng)化研磨高性能加工重點(diǎn)實(shí)驗(yàn)室 c.廣東省強(qiáng)化研磨高性能微納加工工程技術(shù)研究中心，廣州 510006）

探索強(qiáng)化研磨不同工藝參數(shù)下定點(diǎn)噴射對(duì)GCr15軸承鋼殘余應(yīng)力場(chǎng)的影響規(guī)律。采用圖像處理技術(shù)分析了不同工藝參數(shù)下強(qiáng)化研磨定點(diǎn)噴射表面覆蓋率的分布特征。采用二維正態(tài)分布函數(shù)描述強(qiáng)化研磨定點(diǎn)噴射下鋼珠的分布特征，運(yùn)用Python/Opencv確定了在不同工藝參數(shù)下有限元模型所需的鋼珠數(shù)量，基于Abaqus/Python構(gòu)建出強(qiáng)化研磨正態(tài)分布有限元模型。運(yùn)用所建立的正態(tài)分布模型分析不同噴射速度、鋼珠直徑及覆蓋率對(duì)殘余應(yīng)力場(chǎng)的影響。當(dāng)噴射速度從45 m/s增加到70 m/s時(shí)，表面殘余壓應(yīng)力從?683.5 MPa增加到?902.4 MPa，最大殘余壓應(yīng)力從?981.6 MPa增加到?1330.6 MPa，殘余壓力層厚度從89 μm增加到151 μm，最大殘余壓應(yīng)力深度從30 μm移動(dòng)到70 μm。當(dāng)鋼珠直徑從0.4 mm增加到1.0 mm時(shí)，表面殘余壓應(yīng)力先增大后減小，最大殘余壓應(yīng)力從?1063.5 MPa增加到?1240.7 MPa，最大殘余壓應(yīng)力深度從30 μm增加到60 μm，殘余壓應(yīng)力層厚度從103 μm增加到147 μm，其中鋼珠直徑從0.8 mm增加到1.0 mm，最大殘余壓應(yīng)力保持不變。當(dāng)噴射覆蓋率從100%到300%時(shí)，表面殘余壓應(yīng)力、最大殘余壓應(yīng)力及最大殘余壓應(yīng)力深度略有增加，殘余壓應(yīng)力層厚度幾乎保持不變。將正態(tài)分布模型、隨機(jī)分布模型仿真值與實(shí)驗(yàn)值進(jìn)行比較，發(fā)現(xiàn)三者的表面殘余壓應(yīng)力、最大殘余壓應(yīng)力深度及殘余壓應(yīng)力厚度幾乎一致，最大殘余壓應(yīng)力隨機(jī)分布模型的仿真值比實(shí)驗(yàn)值高32.1%，正態(tài)分布模型的仿真值比實(shí)驗(yàn)值高18.9%。強(qiáng)化研磨正態(tài)分布有限元模型能夠較為準(zhǔn)確地預(yù)測(cè)殘余應(yīng)力變化過(guò)程，能夠?yàn)閺?qiáng)化研磨工藝參數(shù)優(yōu)化提供一定的指導(dǎo)。

強(qiáng)化研磨；殘余應(yīng)力場(chǎng)；正態(tài)分布；數(shù)值模擬；GCr15軸承鋼

強(qiáng)化研磨是一種金屬材料精密加工技術(shù)，其加工原理為連續(xù)不斷的鋼珠、研磨粉和研磨液混合組成的研磨料在高壓氣體的驅(qū)動(dòng)下對(duì)材料表面進(jìn)行撞擊，從而使材料表面形成有利于抗疲勞、抗腐蝕、抗磨損的表面類織構(gòu)層（油囊、紋理）[1-2]。在強(qiáng)化研磨過(guò)程中，由于鋼珠的高速?zèng)_擊產(chǎn)生強(qiáng)化作用，會(huì)使材料表面引入殘余壓應(yīng)力，殘余壓應(yīng)力對(duì)裂紋擴(kuò)展起抑制作用，能顯著提高材料的表面性能。在課題組的前期研究中，通過(guò)實(shí)驗(yàn)方法探索了不同工藝參數(shù)下工件的殘余應(yīng)力分布[3]，但在強(qiáng)化研磨殘余應(yīng)力場(chǎng)數(shù)值模擬方面鮮有學(xué)者進(jìn)行研究，目前缺少一種能可靠預(yù)測(cè)強(qiáng)化研磨殘余應(yīng)力場(chǎng)的有限元模型，導(dǎo)致強(qiáng)化研磨實(shí)驗(yàn)加工之后仍要進(jìn)行大量實(shí)驗(yàn)，這嚴(yán)重制約了強(qiáng)化研磨技術(shù)的研究效率。因此，建立一種可靠的強(qiáng)化研磨殘余應(yīng)力場(chǎng)的數(shù)值模擬模型至關(guān)重要。

近年來(lái)，賀占蜀等[4]、李雁淮等[5]建立了單丸粒和雙丸粒噴丸模型，研究了不同噴丸參數(shù)對(duì)殘余應(yīng)力場(chǎng)的影響。董星等[6]、王永軍等[7]、MAJZOOBI等[8]、HAN K等[9]分別建立了周向等距覆蓋、陣列覆蓋、對(duì)稱覆蓋、多層覆蓋等多種多彈丸規(guī)律覆蓋噴丸模型，分別采用不同簡(jiǎn)化方法模擬了噴丸動(dòng)態(tài)碰撞過(guò)程，為噴丸有限元模型探索提供了多種思路。何嘉禧等[10]、LIN Qin-jie等[11]建立了噴丸三維隨機(jī)覆蓋仿真模型，探究了不同噴丸工藝參數(shù)對(duì)其表面完整性的影響。AHMAD等[12]、洪滔等[13]運(yùn)用離散元的方法對(duì)彈丸碰撞過(guò)程建立了FEM-DEM耦合仿真模型，研究了不同覆蓋率對(duì)其殘余應(yīng)力場(chǎng)的影響。然而，在實(shí)驗(yàn)中，為避免噴嘴堵塞，會(huì)使噴嘴隔靶材一段距離，因此強(qiáng)化研磨定點(diǎn)噴射表面覆蓋率的分布并不是隨機(jī)均勻分布的，而是越靠近中心區(qū)域覆蓋率越高。眾多研究人員通常定點(diǎn)噴射一定時(shí)間后，取靶材受噴中心覆蓋率達(dá)到98%及以上區(qū)域作為樣品，再測(cè)量樣品的殘余應(yīng)力等實(shí)驗(yàn)值。在實(shí)際中，通常也有需要進(jìn)行定點(diǎn)噴射的零件。在定點(diǎn)噴射過(guò)程中，測(cè)試區(qū)除了本區(qū)域內(nèi)鋼珠會(huì)對(duì)此區(qū)域產(chǎn)生碰撞強(qiáng)化，測(cè)試區(qū)之外的鋼珠與靶材碰撞產(chǎn)生應(yīng)力波的傳遞和碰撞變形，也會(huì)對(duì)觀測(cè)區(qū)的殘余應(yīng)力場(chǎng)、晶粒尺寸、位錯(cuò)密度等強(qiáng)化效果評(píng)價(jià)指標(biāo)值產(chǎn)生影響。以上研究中，單丸粒、多丸粒碰撞模型均不能體現(xiàn)強(qiáng)化研磨加工碰撞中鋼珠在空間分布的隨機(jī)性，而三維隨機(jī)碰撞模型、FEM-DEM耦合仿真模型只研究了鋼珠在空間上的隨機(jī)分布，不能體現(xiàn)強(qiáng)化研磨定點(diǎn)噴射表面覆蓋率逐漸遞減的隨機(jī)分布。

此外，眾多學(xué)者研究了噴丸表面覆蓋率對(duì)強(qiáng)化效果各評(píng)價(jià)指標(biāo)的影響。研究結(jié)果表明，隨表面覆蓋率的增加，其殘余應(yīng)力場(chǎng)各特征參數(shù)、位錯(cuò)密度等噴丸強(qiáng)化效果評(píng)價(jià)指標(biāo)值均緩慢增大[14-15]。但是，當(dāng)表面覆蓋率增加到一定程度時(shí)會(huì)使材料表面循環(huán)軟化，部分位錯(cuò)纏結(jié)消失，殘余應(yīng)力更容易松馳，從而降低噴丸強(qiáng)化效果[16-17]?？梢钥闯觯瑥?qiáng)化研磨工藝的強(qiáng)化效果和效率與表面覆蓋率是否合理有很大關(guān)系，而強(qiáng)化研磨定點(diǎn)噴射試驗(yàn)過(guò)程中表面覆蓋率從噴射中心逐漸遞減，各區(qū)域表面覆蓋率分布不一致，當(dāng)多點(diǎn)強(qiáng)化研磨時(shí)會(huì)使表面覆蓋率的均勻程度難以控制，容易造成某些區(qū)域表面覆蓋率過(guò)高或者過(guò)低，從而影響工件的抗疲勞性能。

綜上所述，目前國(guó)內(nèi)外學(xué)者為研究噴丸動(dòng)態(tài)碰撞過(guò)程探索了多種有限元仿真模型，且取得了相應(yīng)的研究成果，但尚未見到體現(xiàn)強(qiáng)化研磨定點(diǎn)噴射表面覆蓋率逐漸遞減的鋼珠隨機(jī)分布的數(shù)值仿真模型。本文采用漸進(jìn)式研究方法，首先研究強(qiáng)化研磨定點(diǎn)噴射表面覆蓋率的分布特征，進(jìn)而建立一種新的符合強(qiáng)化研磨定點(diǎn)噴射表面覆蓋率逐漸遞減分布的隨機(jī)碰撞仿真模型，利用新的仿真模型探究不同強(qiáng)化研磨工藝參數(shù)對(duì)工件殘余應(yīng)力場(chǎng)變化規(guī)律的影響，最后通過(guò)實(shí)驗(yàn)驗(yàn)證仿真模型的合理性。旨在探索一種新的強(qiáng)化研磨定點(diǎn)噴射殘余應(yīng)力場(chǎng)數(shù)值模擬方法，更準(zhǔn)確地預(yù)測(cè)定點(diǎn)噴射中殘余應(yīng)力場(chǎng)隨不同強(qiáng)化研磨工藝參數(shù)的變化過(guò)程，以提高強(qiáng)化研磨加工技術(shù)的研究效率。

1 強(qiáng)化研磨定點(diǎn)噴射表面覆蓋率的分布特征

為了建立強(qiáng)化研磨鋼珠的隨機(jī)分布模型，首先需獲取強(qiáng)化研磨定點(diǎn)噴射靶材表面覆蓋率的分布特征。實(shí)驗(yàn)設(shè)備為廣東省強(qiáng)化研磨高性能微納加工工程技術(shù)研究中心自主研制的強(qiáng)化研磨機(jī)（如圖1所示），實(shí)驗(yàn)材料為100 mm×75 mm×10 mm熱處理后的GCr15軸承鋼板，GCr15軸承鋼板的硬度為55HRC。鋼珠在高壓氣體的驅(qū)動(dòng)下經(jīng)文丘里管噴嘴以一定速度撞擊靶材表面，在撞擊區(qū)域內(nèi)，凹坑密度從撞擊區(qū)域中心逐漸遞減，且撞擊區(qū)域呈圓形。強(qiáng)化研磨工藝參數(shù)為：文丘里管管口直徑8 mm，噴嘴到受噴靶材表面的距離80 mm，鋼珠直徑0.8 mm，噴射角度90°，噴射壓力0.5 MPa，噴射流量2.4 kg/min，加工時(shí)間8 s。經(jīng)強(qiáng)化研磨加工后，對(duì)樣品表面進(jìn)行二值化處理，如圖2a所示。參照文獻(xiàn)[18]的單次噴丸表面覆蓋率分布特征的確定方法，在靶材的圓形受噴表面選取4個(gè)區(qū)域，分別記為Ⅰ、Ⅱ、Ⅲ和Ⅳ，每個(gè)區(qū)域內(nèi)有若干個(gè)采樣點(diǎn)，如圖2b所示，其中區(qū)域Ⅰ的采樣點(diǎn)直徑為6 mm，其他區(qū)域采樣點(diǎn)的直徑均為3 mm。

圖1 強(qiáng)化研磨加工設(shè)備

圖2 沿寬度方向上強(qiáng)化研磨定點(diǎn)噴射表面覆蓋率的測(cè)試方法

圖3 強(qiáng)化研磨定點(diǎn)噴射表面覆蓋率沿寬度方向上的分布規(guī)律

一維表面上覆蓋率正態(tài)分布概率密度函數(shù)為：

保持其他強(qiáng)化研磨加工參數(shù)不變，利用上述表面覆蓋率確定方法測(cè)得不同鋼珠直徑、不同噴射壓力下強(qiáng)化研磨定點(diǎn)噴射表面覆蓋率分布特征，如圖4所示。其中P0.5D0.8代表壓力0.5 MPa、鋼珠直徑0.8 mm。以下加工參數(shù)都將采用此種縮寫。

2 仿真建模

強(qiáng)化研磨加工是一個(gè)復(fù)雜的循環(huán)動(dòng)態(tài)碰撞接觸過(guò)程，涉及到材料非線性、幾何非線性、邊界非線性等問(wèn)題，建立完全符合實(shí)際情況的仿真模型極為困難。因此，為了簡(jiǎn)化仿真計(jì)算過(guò)程，在不影響研究結(jié)果的前提下，眾多學(xué)者做出了如下假設(shè)[2,6]：（1）強(qiáng)化研磨噴射過(guò)程中起強(qiáng)化作用的主要是鋼珠對(duì)靶材的沖擊碰撞，研磨粉與研磨液的影響可忽略不計(jì)；（2）強(qiáng)化研磨噴射過(guò)程中所有鋼球沖擊靶材速度的大小和方向均相同，無(wú)速度差；（3）強(qiáng)化研磨噴射過(guò)程中忽略鋼球間的相互作用；（4）靶面為理想光滑平面；（5）靶體為均質(zhì)的各向同性材料。強(qiáng)化研磨定點(diǎn)噴射仿真模型建立過(guò)程如圖5所示。

圖5 強(qiáng)化研磨定點(diǎn)噴射仿真模型建立過(guò)程

2.1 幾何模型

2.1.1 靶材模型

為了提高計(jì)算效率，將仿真模型等比例縮小，由于實(shí)驗(yàn)受噴區(qū)域呈正態(tài)分布，所以等比例縮小的仿真模型也呈對(duì)應(yīng)的正態(tài)分布?；谏逃糜邢拊浖嗀baqus/CAE，采用Python程序建立鋼珠隨機(jī)正態(tài)分布有限元模型。靶材模型尺寸為=3.5 mm、=1.5 mm的圓柱，分為觀測(cè)區(qū)（區(qū)域1）、沖擊區(qū)（區(qū)域2）、過(guò)渡區(qū)（區(qū)域3）、無(wú)限元區(qū)（區(qū)域4）等4個(gè)區(qū)域。觀測(cè)區(qū)為殘余應(yīng)力場(chǎng)結(jié)果的提取區(qū)域，覆蓋率在98%及以上。在P0.5D0.8的強(qiáng)化研磨噴射條件下，觀測(cè)區(qū)大小為=0.43 mm、=0.6 mm。不同強(qiáng)化研磨加工參數(shù)下觀測(cè)區(qū)的不一樣，為不同加工參數(shù)的正態(tài)分布曲線下覆蓋率達(dá)到98%及以上的強(qiáng)化研磨定點(diǎn)噴射寬度。沖擊區(qū)為鋼珠束沖擊區(qū)域，即大部分鋼珠的質(zhì)心被限制在該區(qū)域中。過(guò)渡區(qū)為包裹著沖擊區(qū)的過(guò)渡區(qū)域，可能會(huì)有極少部分鋼珠分布在此區(qū)域內(nèi)。最外圍區(qū)域?yàn)闊o(wú)限元區(qū)，由于目標(biāo)的幾何尺寸遠(yuǎn)小于實(shí)際構(gòu)件尺寸，會(huì)產(chǎn)生應(yīng)力波反射重新進(jìn)入模型而影響仿真結(jié)果的準(zhǔn)確性，所以在最外圍區(qū)域使用了無(wú)限元單元作為反射邊界。圖6為三維仿真模型。

圖6 三維強(qiáng)化研磨正態(tài)分布隨機(jī)碰撞模型

2.1.2 鋼珠分布模型

2.1.2.1 鋼珠分布函數(shù)的確定

式中：為鋼珠數(shù)；為單顆鋼珠噴射引起的凹坑直徑；為靶材料對(duì)應(yīng)的沖擊區(qū)域直徑。

鋼珠任意點(diǎn)（,）二維正態(tài)分布概率密度函數(shù)表示為：

式中：=(,)t，和均為一維正態(tài)隨機(jī)變量；為的期望，=(μ, μ)t；是的協(xié)方差矩陣，二維正態(tài)分布函數(shù)協(xié)方差矩陣計(jì)算公式為：

在噴射角度為90°的強(qiáng)化研磨定點(diǎn)噴射中，μ=μ==0，(,)=(,)=22=0.682，(,) =(,)=0。

圖7 鋼珠坐標(biāo)點(diǎn)（x,y）服從二維正態(tài)分布

2.1.2.2 鋼珠數(shù)量及三維坐標(biāo)的確定

強(qiáng)化研磨覆蓋率是一個(gè)重要的工藝參數(shù)，它是根據(jù)強(qiáng)化研磨噴射后表面所有凹坑面積與表面受噴總面積的比值來(lái)計(jì)算的。在工程實(shí)踐中，認(rèn)為當(dāng)覆蓋率達(dá)到98%時(shí)，表面已被完全覆蓋，只需將強(qiáng)化研磨時(shí)間調(diào)整為完全覆蓋所需的2倍或3倍，即可獲得200%或300%的覆蓋率[22]。在這項(xiàng)研究中，假設(shè)鋼珠是一個(gè)球體，并且在擊中目標(biāo)后，表面會(huì)形成一個(gè)半徑為的規(guī)則圓形凹坑。利用Python/Opencv圖像處理，如圖8a所示，在圓心（0,0）半徑內(nèi)整個(gè)區(qū)域?yàn)榧t色像素點(diǎn)，然后遵循二維正態(tài)分布隨機(jī)生成圓心為（[][],[][]）、半徑為的黑色圓。Python程序中鋼珠（[][],[][]）坐標(biāo)二維正態(tài)分布函數(shù)如式（6）所示，每生成一個(gè)黑色圓，利用ROI提取區(qū)域內(nèi)的圖像信息生成圖8b，計(jì)算1次半徑內(nèi)的Coverage值，Coverage值計(jì)算如式（7）所示。當(dāng)Coverage值小于時(shí)，則在生成了個(gè)黑色圓的圖像基礎(chǔ)上繼續(xù)生成第+1個(gè)黑色圓（遵循二維正態(tài)分布函數(shù)隨機(jī)生成一個(gè)），然后循環(huán)上述命令去計(jì)算Coverage值，當(dāng)Coverage值大于或等于（=98%）時(shí)，即退出此次循環(huán)，記錄單次循環(huán)內(nèi)鋼珠的全部坐標(biāo)（[][],[][]），及整個(gè)模型所需要的鋼珠總數(shù)

圖8 強(qiáng)化研磨定點(diǎn)噴射表面覆蓋率的計(jì)算方法

和影響Coverage值的鋼珠數(shù)，鋼珠數(shù)量確定的整個(gè)過(guò)程如圖9所示。

式中：np.random.multivariate_normal為python語(yǔ)言中的二維正態(tài)分布函數(shù)；1為黑色像素點(diǎn)數(shù)；2為紅色像素點(diǎn)數(shù)。

圖9 鋼珠數(shù)量確定流程

圖10描繪了6種不同強(qiáng)化研磨加工參數(shù)下遵循二維正態(tài)分布在區(qū)域內(nèi)覆蓋率達(dá)到98%、循環(huán)1000次下所需的鋼珠數(shù)的分布，發(fā)現(xiàn)在固定覆蓋條件下所需鋼珠總數(shù)量的概率分布近似呈正態(tài)分布。LIN Qin-jie等[11]在使用另外一種統(tǒng)計(jì)法計(jì)算隨機(jī)分布下鋼珠的數(shù)量時(shí)也有類似的分布特征。Avrami方程[23]是一種計(jì)算隨機(jī)分布所需鋼珠數(shù)量常用的經(jīng)驗(yàn)方法，用公式（2）—（3）表示。圖11將所提出的Python/ Opencv統(tǒng)計(jì)方法在特定區(qū)域?qū)崿F(xiàn)全覆蓋過(guò)程循環(huán)1000次的鋼珠數(shù)量和的平均數(shù)、中位數(shù)與Avrami方程進(jìn)行了比較，可以清楚地看到，循環(huán)1000次的鋼珠數(shù)量和的平均數(shù)、中位數(shù)均近似相等；在覆蓋面積為98%的區(qū)域內(nèi)，python統(tǒng)計(jì)法記錄的略高于公式法，這是因?yàn)檎龖B(tài)分布模型越靠近中心區(qū)域，鋼珠生成的概率越大，生成的概率從圓心處逐漸遞減，而Avrami方程是用于隨機(jī)分布的計(jì)算，每個(gè)區(qū)域生成的鋼珠概率相等；其次是因?yàn)閜ython統(tǒng)計(jì)法中計(jì)算的為小于+內(nèi)的圓心點(diǎn)，所以會(huì)略高于公式法。

圖10 循環(huán)1000次所需鋼珠總數(shù)的分布特征

圖11 統(tǒng)計(jì)方法與Avrami方程達(dá)到特定覆蓋率所需的鋼珠總數(shù)及其比較

在LIN Qin-jie等[11]針對(duì)隨機(jī)分布的鋼珠數(shù)量的研究中，取循環(huán)1000次所需鋼珠數(shù)量的平均數(shù)作為仿真模型的鋼珠數(shù)量，何嘉禧等[10]利用Avrami方程計(jì)算的鋼珠數(shù)量作為仿真模型所需的鋼珠數(shù)量，然后在仿真模型中利用隨機(jī)函數(shù)生成鋼球。但是從他們的研究結(jié)果及本文鋼珠數(shù)量的Python統(tǒng)計(jì)方法可知，達(dá)到特定覆蓋率所需的鋼珠數(shù)量并不固定，而是呈正態(tài)分布。如果將循環(huán)1000次的鋼珠數(shù)量平均數(shù)或者Avrami方程計(jì)算的鋼珠數(shù)量作為仿真模型所需要的鋼珠數(shù)量，在模型中利用隨機(jī)函數(shù)或者正態(tài)分布函數(shù)生成鋼珠，所生成的鋼珠不一定能達(dá)到特定覆蓋率，且每次仿真模型生成的鋼珠數(shù)量達(dá)到的覆蓋率都不一樣，會(huì)存在仿真誤差。為消除仿真誤差，本文取循環(huán)1000次去掉最高值，即循環(huán)999次的鋼珠總數(shù)量的中位數(shù)作為Abaqus模擬計(jì)算的鋼珠數(shù)。在Python/Opencv運(yùn)行結(jié)果中找到此次鋼珠數(shù)量下記錄的鋼珠的全部坐標(biāo)（[][],[][]），然后導(dǎo)入Abaqus中，方向遵循隨機(jī)分布，如式（8）所示，且2個(gè)鋼珠之間的距離＞2，基于Abaqus/Python建立了新的強(qiáng)化研磨隨機(jī)正態(tài)分布沖擊有限元模型。根據(jù) Klemenz[24]的經(jīng)驗(yàn)公式（9），可以建立實(shí)驗(yàn)參數(shù)與仿真參數(shù)的關(guān)系。

式中：random.uniform為Python語(yǔ)言中的隨機(jī)函數(shù)；為鋼珠的平均速度，m/s；為噴射壓力，MPa；為鋼珠的直徑，mm；為鋼球的進(jìn)給流量，kg/min。

本次強(qiáng)化研磨實(shí)驗(yàn)中，=2.4 kg/min。在D0.8P0.3、D0.8P0.5、D0.8P0.7強(qiáng)化研磨加工參數(shù)下，對(duì)應(yīng)的仿真參數(shù)分別約為D0.8V45、D0.8V60、D0.8V75，單顆鋼珠造成的凹坑直徑分別為0.24、0.26、0.28 mm。在D0.4P0.5、D0.6P0.5、D1.0P0.5強(qiáng)化研磨加工參數(shù)下，對(duì)應(yīng)的平均速度均約為60 m/s，對(duì)應(yīng)的凹坑直徑分別為0.14、0.2、0.3 mm。

2.2 材料模型

強(qiáng)化研磨加工過(guò)程中，鋼珠高速撞擊靶材表面，靶材受沖擊表面會(huì)在極短時(shí)間內(nèi)發(fā)生劇烈的塑性變形，在劇烈的塑性變形中材料發(fā)生應(yīng)變硬化、應(yīng)變率強(qiáng)化和溫度軟化效應(yīng)，而在Johnson-cook材料動(dòng)態(tài)本構(gòu)模型均有體現(xiàn)，能夠很好地用來(lái)描述材料在高應(yīng)變率條件下應(yīng)力-應(yīng)變的關(guān)系[25]，因此在強(qiáng)化研磨仿真模型中應(yīng)用該材料本構(gòu)模型，其具體表達(dá)式為：

表1 GCr15軸承鋼J-C本構(gòu)模型參數(shù)

Tab.1 J-C constitutive model parameters of GCr15 bearing steel

2.3 網(wǎng)格劃分與邊界條件

網(wǎng)格劃分上，在本模型中，鋼珠采用四節(jié)點(diǎn)線性四面體單元C3D4，靶材無(wú)限元區(qū)的網(wǎng)格為無(wú)限單元CIN3D8，靶材其余區(qū)域?yàn)榘斯?jié)點(diǎn)線性六面體減縮積分單元C3D8R，且進(jìn)行沙漏控制；對(duì)觀測(cè)區(qū)×內(nèi)進(jìn)行局部網(wǎng)格細(xì)化，網(wǎng)格最小尺寸為0.03 mm×0.03 mm× 0.01 mm；對(duì)靶材底面自由度進(jìn)行完全約束。鋼珠與靶材表面的相對(duì)切向運(yùn)動(dòng)使用庫(kù)倫摩擦模型描述，定義庫(kù)倫摩擦因數(shù)為0.25[27]。鋼珠與靶材的接觸為面-面接觸，力學(xué)約束公式化為罰接觸方法。

3 結(jié)果與分析

3.1 噴射速度對(duì)殘余應(yīng)力分布的影響

強(qiáng)化研磨之后，零件表層區(qū)域的殘余應(yīng)力場(chǎng)對(duì)材料表面性能的改善十分重要。因?yàn)閺?qiáng)化研磨主要是通過(guò)鋼珠對(duì)零件表面的沖擊來(lái)強(qiáng)化，所以強(qiáng)化研磨殘余應(yīng)力場(chǎng)的表征方式通常參考噴丸的描述方式，以表面殘余壓應(yīng)力、最大殘余壓應(yīng)力、最大殘余壓應(yīng)力深度和殘余壓應(yīng)力層厚度這4個(gè)特征參數(shù)的變化來(lái)描述殘余應(yīng)力場(chǎng)的分布特征[3]。殘余壓應(yīng)力對(duì)裂紋擴(kuò)展起抑制作用，能夠明顯地改善材料的表面性能。

在殘余壓力的測(cè)量過(guò)程中，對(duì)強(qiáng)化研磨后的試樣表面利用電腐蝕剝層處理。然后利用X射線衍射法測(cè)量不同深度的殘余應(yīng)力值，通常X射線的測(cè)量光斑直徑為1 mm，測(cè)量得到的殘余應(yīng)力值為同深度上X射線測(cè)量光斑范圍內(nèi)的應(yīng)力平均值。所以，在有限元仿真中，為了準(zhǔn)確地描述殘余應(yīng)力場(chǎng)在每個(gè)深度上的結(jié)果，利用Python程序提取觀測(cè)區(qū)在該深度水平上每層節(jié)點(diǎn)的平均應(yīng)力值，從而得到沿深度方向殘余應(yīng)力的分布曲線。

圖12顯示了利用建立的隨機(jī)正態(tài)分布仿真模型在=0.8 mm、=100%時(shí)不同噴射速度下殘余應(yīng)力的分布。結(jié)果表明，隨著噴射速度從45 m/s增加到75 m/s，殘余壓力層厚度從89 μm增加到151 μm。這表明強(qiáng)化研磨噴射速度的增加可以獲得更厚的殘余壓應(yīng)力層；同時(shí)，隨著噴射速度的增加，最大殘余壓應(yīng)力出現(xiàn)在較深的位置，從30 μm移動(dòng)到70 μm；最大殘余壓應(yīng)力深度和殘余應(yīng)力層厚度的增加是由于噴射速度的增加，使鋼珠攜帶的動(dòng)能增大，鋼珠沖擊材料表面，從而使材料表層塑性變形滲透層的深度增加[28]。另外，當(dāng)噴射速度從45 m/s增大到75 m/s時(shí)，表面殘余壓應(yīng)力從?683.5 MPa增加到?902.4 MPa，最大殘余壓應(yīng)力從?981.6 MPa增加到?1330.6 MPa。這是由于鋼珠動(dòng)能的增加使材料表層塑性變形程度增大，從而造成了表面殘余壓應(yīng)力和最大殘余壓應(yīng)力增加。從圖12中可以觀察到，噴射速度從60 m/s增大到75 m/s時(shí)，表面殘余壓應(yīng)力和最大殘余壓應(yīng)力的增加幅度較低。這是因?yàn)镚Cr15軸承鋼為高強(qiáng)度鋼，屈服強(qiáng)度比較高，不容易發(fā)生塑性變形，其值隨著速度的增加容易達(dá)到飽和狀態(tài)，這與鐘軼寧等人[28]的研究結(jié)果一致。另外有研究表明，表面殘余壓應(yīng)力值和次表層的最大殘余壓應(yīng)力值的最主要影響因素是材料的力學(xué)性能，而不是噴射強(qiáng)度[29]。

圖12 不同噴射速度下殘余應(yīng)力沿深度方向上的分布情況

3.2 鋼珠直徑對(duì)殘余應(yīng)力的影響

圖13描述了在=60 m/s、=100%時(shí)不同鋼珠直徑下殘余應(yīng)力的變化。從圖13中可以觀察到，當(dāng)鋼珠直徑從0.4 mm增加到1.0 mm時(shí)，最大殘余壓應(yīng)力從?1063.5 MPa增加到?1240.7 MPa，最大殘余壓應(yīng)力深度從30 μm增加到60 μm，殘余壓應(yīng)力層厚度從103 μm增加到147 μm。這是因?yàn)?，由鋼珠直徑與沖擊能量的關(guān)系式[30]可知，鋼珠的動(dòng)能與鋼珠直徑的立方成正比，當(dāng)鋼珠直徑增大時(shí)，鋼珠攜帶的動(dòng)能顯著增大，靶材表面因?yàn)槭艿戒撝闆_擊造成塑性變形的程度增大及塑性變形滲透層的影響深度顯著增大；其次，根據(jù)赫茲接觸理論可知，最大剪切應(yīng)力出現(xiàn)在次表層[31]，因此導(dǎo)致次表層的殘余壓應(yīng)力增大明顯。值得注意的是，從圖13中可以觀察到，鋼珠直徑從0.8 mm增加到1.0 mm時(shí)，表面殘余壓應(yīng)力反而減小，最大殘余壓應(yīng)力幾乎不變。這是因?yàn)樵趩挝幻娣e內(nèi)總能量相差不是很大的情況下，鋼珠直徑越大，則與靶材表面的接觸面積越大，會(huì)造成最表層區(qū)域塑性變形程度相對(duì)較小，次表面塑性變形總體程度幾乎不變，但又因?yàn)殡S直徑的增大，單顆鋼珠攜帶的動(dòng)能增大，所以塑性變形滲透層的影響深度增大，從而表現(xiàn)為鋼珠直徑太大表面殘余壓應(yīng)力反而降低，次表層殘余壓應(yīng)力高，最大殘余壓應(yīng)力深度向基體內(nèi)移動(dòng)。在張廣良[32]、ZHANG等[33]的研究中也表明，在總能量大致相同的情況下，鋼珠直徑越小，表面塑性變形程度越高，表面殘余壓應(yīng)力越大。

針對(duì)上述的研究結(jié)果，眾多學(xué)者研究了二次噴丸技術(shù)，即一次噴丸通過(guò)大尺寸鋼珠噴射，之后第二次噴丸采用小尺寸鋼珠噴射，這樣可以提高噴丸效果[10-11]，進(jìn)一步提高材料的表面性能[34-35]。

圖13 不同鋼珠直徑下殘余應(yīng)力沿深度方向上的分布情況

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

圖14描述了在=0.8 mm、=60 m/s時(shí)不同覆蓋率下殘余應(yīng)力場(chǎng)的變化。從圖14中可以觀察到，當(dāng)覆蓋率從100%增加到300%時(shí)，表面殘余壓應(yīng)力的增量很小，僅從?827 MPa增加到?898 MPa。這說(shuō)明表面覆蓋率的增加對(duì)表面殘余壓應(yīng)力的影響不大，這與LIN等[11]的研究結(jié)果一致。當(dāng)覆蓋率從100%增加到200%時(shí)，最大殘余壓應(yīng)力從?1236.4 MPa增加到?1305.6 MPa，最大殘余壓應(yīng)力深度略有增加，在60～110 μm深度范圍內(nèi)，殘余壓應(yīng)力均增幅較大，這可能是因?yàn)橄噜弲^(qū)域覆蓋率增加，使得相鄰區(qū)域內(nèi)的殘余拉應(yīng)力區(qū)域逐漸消失并且轉(zhuǎn)化為殘余壓應(yīng)力，所形成的殘余壓應(yīng)力層的影響深度與全覆蓋觀測(cè)區(qū)域的殘余壓應(yīng)力層的影響深度逐漸重疊；其次是因?yàn)殡S著覆蓋率的增加，材料表面循環(huán)硬化程度增強(qiáng)，用于塑性變形的能量減少，動(dòng)能轉(zhuǎn)化為材料內(nèi)部彈性能增多，赫茲壓力效應(yīng)更加顯著[36]，使得最大殘余壓應(yīng)力往內(nèi)移動(dòng)和殘余壓應(yīng)力增大。當(dāng)覆蓋率從200%增加至300%時(shí)，材料殘余壓應(yīng)力的增加量很小，逐漸趨向于穩(wěn)定，這是因?yàn)樵诖藝娚渌俣认職堄鄳?yīng)力場(chǎng)逐漸飽和所致。同時(shí)，隨著覆蓋率從100%增加到300%，殘余壓力層厚度保持不變。這是因?yàn)楦采w率的增加并沒有改變單顆鋼球的動(dòng)能，鋼珠動(dòng)能轉(zhuǎn)化為影響應(yīng)力波滲透層深度的能量不變，沒有足夠的能量去影響更深層區(qū)域的殘余應(yīng)力的分布。因此，覆蓋率的增加不能改變殘余壓應(yīng)力層厚度，這與WU Ji-zhan等[37]的發(fā)現(xiàn)相符。

圖14 不同強(qiáng)化研磨覆蓋率下殘余應(yīng)力沿深度方向上的分布情況

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

為了驗(yàn)證正態(tài)分布仿真模型的合理性，將隨機(jī)分布模型、正態(tài)分布模型的殘余應(yīng)力模擬結(jié)果與實(shí)驗(yàn)結(jié)果進(jìn)行對(duì)比。通過(guò)強(qiáng)化研磨機(jī)對(duì)GCr15軸承鋼進(jìn)行強(qiáng)化研磨處理，噴射壓力為0.5 MPa，鋼珠直徑為0.8 mm，加工時(shí)間為8 s，其他加工參數(shù)見第1節(jié)，然后取覆蓋率達(dá)到98%及以上的試樣中心區(qū)域進(jìn)行殘余應(yīng)力檢測(cè)實(shí)驗(yàn)的制樣。殘余應(yīng)力采用Proto I XRD射線應(yīng)力分析儀檢測(cè)。采用Cr靶Kα輻射，管電流為5 mA，管電壓為25 kV。通過(guò)控制電解腐蝕時(shí)間對(duì)試樣進(jìn)行剝層處理，以獲取試樣在不同深度下的殘余應(yīng)力值。

圖15為強(qiáng)化研磨后GCr15軸承鋼的殘余應(yīng)力仿真與實(shí)驗(yàn)結(jié)果。由圖15可知，正態(tài)分布模型、隨機(jī)分布模型與實(shí)驗(yàn)結(jié)果，在殘余應(yīng)力曲線的形狀與趨勢(shì)方面基本一致，三者表面殘余壓應(yīng)力、最大殘余壓應(yīng)力深度和殘余壓應(yīng)力層厚度基本一致，殘余應(yīng)力層厚度為120～140 μm，最大殘余壓應(yīng)力深度為40～50 μm，但正態(tài)分布模型與隨機(jī)分布模型的結(jié)果值都略高于實(shí)驗(yàn)值。這可能是因?yàn)?，在仿真模型中，為了?jì)算效率通常會(huì)把鋼珠約束為剛體，而實(shí)驗(yàn)中由于GCr15軸承鋼屬于高強(qiáng)度鋼，即使鋼珠與GCr15軸承鋼靶體有一定的硬度差，但在高速?zèng)_擊下鋼珠還是會(huì)發(fā)生塑性變形，鋼珠的塑性變形會(huì)消耗自身大量的動(dòng)能，所以會(huì)造成仿真模擬中靶體的塑性變形程度大于實(shí)驗(yàn)中強(qiáng)化研磨試樣的塑性變形程度，從而使仿真模型殘余壓應(yīng)力值偏大。WU Gang等[38]研究了多彈丸仿真模型下不同鋼珠硬度對(duì)殘余應(yīng)力分布的影響，發(fā)現(xiàn)不同硬度和約束為剛體狀態(tài)下的鋼珠對(duì)靶材沖擊時(shí)，殘余應(yīng)力的分布曲線相差較大。在仿真模擬中，鋼珠是以確定的速度和垂直的角度沖擊靶材，而噴射式強(qiáng)化研磨實(shí)驗(yàn)中使用文丘里管噴嘴且距離待加工工件有一定的距離，會(huì)使實(shí)驗(yàn)中鋼珠的噴射速度和噴射角度具有一定的隨機(jī)性，這2種誤差的共同作用下，造成實(shí)驗(yàn)值偏小。

圖15 強(qiáng)化研磨后GCr15軸承鋼殘余應(yīng)力仿真值與實(shí)驗(yàn)值

另外，正態(tài)分布模型仿真值的最大殘余壓應(yīng)力為?1236.4 MPa，實(shí)驗(yàn)值最大殘余壓應(yīng)力為?1040 MPa，模擬誤差為18.9%，隨機(jī)分布模型仿真值的最大殘余壓應(yīng)力為?1374.3 MPa，模擬誤差為32.1%?？梢?，正態(tài)分布模型仿真值與實(shí)驗(yàn)結(jié)果更加接近，而隨機(jī)分布模型與實(shí)驗(yàn)結(jié)果的誤差較大。這可能是因?yàn)?，如圖16所示不同仿真模型下軸平面的殘余應(yīng)力分布云圖，在鋼珠的高速?zèng)_擊下，材料表層發(fā)生了劇烈的塑性變形，為了抵抗塑性變形，材料表層會(huì)產(chǎn)生殘余應(yīng)力，材料表面表現(xiàn)為一定深度的殘余壓應(yīng)力場(chǎng)，同時(shí)材料內(nèi)部為了達(dá)到平衡，會(huì)在更深層區(qū)域形成一個(gè)殘余拉應(yīng)力場(chǎng)；在正態(tài)分布模型中，相鄰區(qū)域由于覆蓋率不高（低于98%），低于98%區(qū)域，在壓應(yīng)力評(píng)估區(qū)域也存在殘余拉應(yīng)力[28,38]，所形成的殘余壓應(yīng)力層厚度會(huì)低于全覆蓋觀測(cè)區(qū)域的殘余壓應(yīng)力層厚度，則在相鄰區(qū)域更深層區(qū)域形成的殘余拉應(yīng)力場(chǎng)所造成的拉應(yīng)力方向的波動(dòng)的影響深度會(huì)和觀測(cè)區(qū)所產(chǎn)生的殘余壓應(yīng)力場(chǎng)的影響深度部分重疊，并會(huì)抵消部分殘余壓應(yīng)力，董星等[6]、何嘉禧等[10]的仿真模型中也表明了存在此現(xiàn)象，所以正態(tài)分布模型的仿真值更接近實(shí)驗(yàn)值。

圖16 不同仿真模型下軸平面殘余應(yīng)力分布云圖

總的來(lái)說(shuō)，可認(rèn)為強(qiáng)化研磨正態(tài)分布模型在一定程度上能夠較為準(zhǔn)確地模擬強(qiáng)化研磨實(shí)驗(yàn)的殘余應(yīng)力變化過(guò)程，這對(duì)強(qiáng)化研磨工藝參數(shù)優(yōu)化具有一定的指導(dǎo)意義。

5 結(jié)論

1）仿真結(jié)果表明，強(qiáng)化研磨能夠使GCr15軸承鋼表層產(chǎn)生一定厚度的殘余壓應(yīng)力層，殘余壓應(yīng)力在次表面達(dá)到最大值，同時(shí)會(huì)在更深層區(qū)域產(chǎn)生殘余拉應(yīng)力層，以使材料內(nèi)部達(dá)到平衡。

2）隨著強(qiáng)化研磨噴射速度的增加，表面殘余壓應(yīng)力、最大殘余壓應(yīng)力深度及殘余壓應(yīng)力層厚度均有明顯提高。隨著鋼珠直徑的增大，表面殘余壓應(yīng)力先增大再減小，但能夠顯著提高最大殘余壓應(yīng)力、最大殘余壓應(yīng)力深度及殘余壓應(yīng)力層厚度。覆蓋率的增加對(duì)殘余應(yīng)力的分布影響最小。

3）強(qiáng)化研磨正態(tài)分布模型殘余應(yīng)力的仿真結(jié)果與實(shí)驗(yàn)結(jié)果有較好的一致性，能夠?qū)?qiáng)化研磨工藝參數(shù)優(yōu)化提供一定的指導(dǎo)。

[1] LIANG Zhong-wei, LIU Xiao-chu, XIAO Jin-rui, et al. Adaptive Prediction of Abrasive Impacting Pressure Effe-ctiveness in Strengthen Jet Grinding Using NSAE- ANFIS[J]. The International Journal of Advanced Manu-facturing Technology, 2020, 106(7-8): 2805-2828.

[2] 劉曉初, 劉鎮(zhèn), 梁忠偉, 等. 基于隨機(jī)碰撞的GCr15鋼強(qiáng)化研磨表面粗糙度數(shù)值模擬[J]. 表面技術(shù), 2021, 50(4): 385-392.

[3] LIU Xiao-chu, LIU Zhen, LIANG Zhong-wei, et al. Nu-merical Simulation of Surface Roughness of GCr15 Steel Strengthened Grinding Based on Random Impact[J]. Sur-face Technology, 2021, 50(4): 385-392.

[4] 黎達(dá)成. 強(qiáng)化研磨中噴射壓力對(duì)軸承鋼殘余應(yīng)力場(chǎng)的影響[D]. 廣州: 廣州大學(xué), 2019.

[5] LI Da-cheng. Influence of Injection Pressure on Residual Stress Field of Bearing Steel in Strengthen Grinding[D]. Guangzhou: Guangzhou University, 2019.

[6] 賀占蜀, 陳雷, 李超, 等. 單丸粒和雙丸粒噴丸模型的有限元模擬[J]. 塑性工程學(xué)報(bào), 2020, 27(4): 153-158.

[7] HE Zhan-shu, CHEN Lei, LI Chao, et al. Finite Element Simulation of Single and Double Shot Peening Model[J]. Journal of Plasticity Engineering, 2020, 27(4): 153-158.

[8] 李雁淮, 王飛, 呂堅(jiān), 等. 單丸粒噴丸模型和多丸粒噴丸模型的有限元模擬[J]. 西安交通大學(xué)學(xué)報(bào), 2007, 41(3): 348-352.

[9] LI Yan-huai, WANG Fei, Lyu Jian, et al. Finite Element Analysis for Shot-Peening Models of Single Shot Impact and Multiple Impingement[J]. Journal of Xi'an Jiaotong University, 2007, 41(3): 348-352.

[10] 董星, 劉雨慶, 段雄. 前混合水射流多彈丸噴丸模型及殘余應(yīng)力場(chǎng)的數(shù)值模擬[J]. 機(jī)械工程學(xué)報(bào), 2020, 56(4): 224-232.

[11] DONG Xing, LIU Yu-qing, DUAN Xiong. Numerical Simulation of Multi-Shot Shot Peening Models and Residual Stress Field of Premixed Water Jet[J]. Journal of Mechanical Engineering, 2020, 56(4): 224-232.

[12] 王永軍, 孫寶龍, 張煒, 等. 鋁鋰合金噴丸強(qiáng)化數(shù)值模擬及試驗(yàn)[J]. 航空動(dòng)力學(xué)報(bào), 2015, 30(3): 595-602.

[13] WANG Yong-jun, SUN Bao-long, ZHANG Wei, et al. Numerical Simulation and Experiment on Shot Peening of Al-Li Alloy[J]. Journal of Aerospace Power, 2015, 30(3): 595-602.

[14] MAJZOOBI G H, AZIZI R, NIA A A. A Three-Dimen-sional Simulation of Shot Peening Process Using Multiple Shot Impacts[J]. Journal of Materials Processing Techno-logy, 2005, 164-165: 1226-1234.

[15] HAN K, OWEN D R J, PERIC D. Combined Finite/ Discrete Element and Explicit/Implicit Simulations of Peen Forming Process[J]. Engineering Computations, 2002, 19(1): 92-118.

[16] 何嘉禧, 汪舟, 甘進(jìn), 等. 二次噴丸42CrMo鋼表面完整性的數(shù)值模擬研究[J]. 表面技術(shù), 2020, 49(6): 216-223.

[17] HE Jia-xi, WANG Zhou, GAN Jin, et al. Numerical Simu-lation on Surface Integrity of 42CrMo Steel after Dual Shot Peening[J]. Surface Technology, 2020, 49(6): 216- 223.

[18] LIN Qin-jie, LIU Huai-ju, ZHU Cai-chao, et al. Effects of Different Shot Peening Parameters on Residual Stress, Surface Roughness and Cell Size[J]. Surface and Coatings Technology, 2020, 398: 126054.

[19] AHMAD A S, WU Yun-xin, GONG Hai. Coupled Finite and Discrete Element Shot Peening Simulation Based on Johnson-Cook Material Model[J]. Proceedings of the Institution of Mechanical Engineers, Part L: Journal of Materials: Design and Applications, 2020, 234(7): 974- 987.

[20] 洪滔, 王志偉, 袁巨龍. 噴丸強(qiáng)化過(guò)程的有限元和離散元模擬[J]. 中國(guó)機(jī)械工程, 2008, 19(11): 1321-1325.

[21] HONG Tao, WANG Zhi-wei, YUAN Ju-long. Numerical Simulation of Shot Peening Process Using Finite Element and Discrete Element Method[J]. China Mechanical En-gineering, 2008, 19(11): 1321-1325.

[22] WU Ji-zhan, LIU Huai-ju, WEI Pei-tang, et al. Effect of Shot Peening Coverage on Residual Stress and Surface Roughness of 18CrNiMo7-6 Steel[J]. International Jour-nal of Mechanical Sciences, 2020, 183: 105785.

[23] GHASEMI A, HASSANI-GANGARAJ S M, MAHMO-UDI A H, et al. Shot Peening Coverage Effect on Residual Stress Profile by FE Random Impact Analysis[J]. Surface Engineering, 2016, 32(11): 861-870.

[24] 王欣, 蔡建明, 王強(qiáng), 等. 噴丸表面覆蓋率對(duì)Ti60高溫鈦合金疲勞性能的影響[J]. 中國(guó)表面工程, 2011, 24(5): 58-63, 103.

[25] WANG Xin, CAI Jian-ming, WANG Qiang, et al. Effect of Shot Peening Surface Coverage on the Fatigue Pro-perty in Ti60 High-Temperature Titanium Alloy[J]. China Surface Engineering, 2011, 24(5): 58-63, 103.

[26] 李瑞鴻, 劉道新, 張煒, 等. 噴丸強(qiáng)化與表面完整性對(duì)300M鋼疲勞性能的影響[J]. 機(jī)械科學(xué)與技術(shù), 2011, 30(9): 1418-1423.

[27] Li Rui-hong, Liu Dao-xin, Zhang Wei, et al. Influence of Shot Peening and Surface Integrity on the Fatigue Properties of 300M Steel[J]. Mechanical Science and Technology, 2011, 30(9): 1418-1423.

[28] 盛湘飛, 李智, 周楠楠, 等. 基于正態(tài)分布的噴丸表面覆蓋均勻程度與強(qiáng)化效率研究[J]. 表面技術(shù), 2018, 47(5): 227-232.

[29] SHENG Xiang-fei, LI Zhi, ZHOU Nan-nan, et al. Unifor-mity of Surface Coverage and Strengthening Efficiency in Shot Peening Process Based on Normal Distribution[J]. Surface Technology, 2018, 47(5): 227-232.

[30] 謝鑫成. 基于強(qiáng)化研磨加工的工業(yè)機(jī)器人軸承內(nèi)圈滾道表面覆蓋率的研究[D]. 廣州: 廣州大學(xué), 2020.

[31] XIE Xin-cheng. Study on Raceway Surface Coverage of Industrial Robot Bearing Inner Ring Based on Strengthen Grinding Processing[D]. Guangzhou: Guangzhou Univer-sity, 2020.

[32] 吳文海, 藍(lán)天, 張?chǎng)? 等. 噴丸過(guò)程中噴嘴移速對(duì)顆粒分布的影響[J]. 表面技術(shù), 2020, 49(8): 342-349.

[33] WU Wen-hai, LAN Tian, ZHANG Ting, et al. Effect of No-zzle Moving Speed on Particle Distribution during Shot Peening[J]. Surface Technology, 2020, 49(8): 342- 349.

[34] OHTA T, INOUE A. Measurement and Functional Appro-ximation of Peening Intensity Distributions[J]. Materials Transactions, 2019, 60(6): 1004-1010.

[35] BAGHERIFARD S, GHELICHI R, GUAGLIANO M. On the Shot Peening Surface Coverage and Its Assessment by Means of Finite Element Simulation: A Critical Review and some Original Developments[J]. Applied Surface Science, 2012, 259: 186-194.

[36] KIRK D, ABYANEH M Y. Theoretical Basis of Shot Peening Coverage Control[J]. The Shot Peener, 1993, 2(9): 28-30.

[37] KLEMENZ M. Anwendung Der Simulation Der Rand-schichtausbildung Beim Kugelstrahlen Auf Die Absch?-t-zung Der Schwingfestigkeit Gekerbter Bauteile[D]. Kar-lsruhe: Karlsruhe University, 2009.

[38] SANJURJO P, RODRíGUEZ C, PE?UELAS I, et al. Influence of the Target Material Constitutive Model on the Numerical Simulation of a Shot Peening Process[J]. Surface and Coatings Technology, 2014, 258: 822-831.

[39] 張方圓. 硬切削GCr15軸承鋼表面變質(zhì)層微觀組織及抗磨損性能研究[D]. 大連: 大連理工大學(xué), 2019.

[40] ZHANG Fang-yuan. Study on Microstructure and Wear Resistance of the Metamorphic Layer of GCr15 Bearing Steel in Hard Cutting[D]. Dalian: Dalian University of Technology, 2019.

[41] 李源, 雷麗萍, 曾攀. 彈丸束噴丸有限元模型數(shù)值模擬及試驗(yàn)研究[J]. 機(jī)械工程學(xué)報(bào), 2011, 47(22): 43-48.

[42] LI Yuan, LEI Li-ping, ZENG Pan. Shot Stream Finite Element Model for Shot Peening Numerical Simulation and Its Experiment Study[J]. Journal of Mechanical Engi-neering, 2011, 47(22): 43-48.

[43] 鐘軼寧. 噴丸處理Ti-6Al-4V合金殘余應(yīng)力場(chǎng)特征與演變規(guī)律研究[D]. 大連: 大連理工大學(xué), 2018.

[44] ZHONG Yi-ning. Analysis of Characteristics and Evolu-tion of Residual Stress Field of Shot Peened Ti-6Al-4V Alloy[D]. Dalian: Dalian University of Technology, 2018.

[45] SEGURADO E, BELZUNCE F J, PARIENTE I F. Effects of Low Intensity Shot Peening Treatments Applied with Different Types of Shots on the Fatigue Performance of a High-Strength Steel[J]. Surface and Coatings Technology, 2018, 340: 25-35.

[46] JOHNSON W. Impact Strength of Materials[M]. UK: Ed-ward Arnold Publishers LTD, 1972.

[47] JOHN K L. Contact Mechanics[M]. Cambridge: Cam-bridge University Press, 1985.

[48] 張廣良. 不同強(qiáng)度鋼噴丸殘余應(yīng)力的有限元模擬[D]. 上海: 上海交通大學(xué), 2012.

[49] ZHANG Guang-liang. Study on Finite Element Simula-tion of Shot Peening of Different Yield Strength Steeland Its Residual Stress[D]. Shanghai: Shanghai Jiao Tong Uni-versity, 2012.

[50] ZHANG Ya-long, LAI Fu-qiang, QU Sheng-guan, et al. Effect of Shot Peening on Residual Stress Distribution and Tribological Behaviors of 17Cr2Ni2MoVNb Steel[J]. Surface and Coatings Technology, 2020, 386: 125497.

[51] 孫寒驍, 朱有利, 侯帥, 等. 疲勞載荷對(duì)17CrNiMo6噴丸強(qiáng)化層殘余應(yīng)力與組織的影響[J]. 中國(guó)表面工程, 2016, 29(4): 43-48.

[52] SUN Han-xiao, ZHU You-li, HOU Shuai, et al. Effects of Fatigue Load on Residual Stress and Microstructure of 17CrNiMo6 Shot Peening Strengthened Layer[J]. China Surface Engineering, 2016, 29(4): 43-48.

[53] WANG Xin, LI Si-qing, YANG Qing, et al. Effect of Dou-ble Shot Peening on Room-Temperature Notched Fatigue Property of TC4 Titanium Alloy[J]. Materials Science Forum, 2015, 817: 90-95.

[54] SHERAFATNIA K, FARRAHI G H, MAHMOUDI A H. Effect of Initial Surface Treatment on Shot Peening Resi-dual Stress Field: Analytical Approach with Experimental Verification[J]. International Journal of Mechanical Scien-ces, 2018, 137: 171-181.

[55] WU Ji-zhan, LIU Huai-ju, WEI Pei-tang, et al. Effect of Shot Peening Coverage on Hardness, Residual Stress and Surface Morphology of Carburized Rollers[J]. Surface and Coatings Technology, 2020, 384: 125273.

[56] WU Gang, WANG Zhou, GAN Jin, et al. FE Analysis of Shot-Peening-Induced Residual Stresses of AISI 304 Stainless Steel by Considering Mesh Density and Friction Coefficient[J]. Surface Engineering, 2019, 35(3): 242- 254.

Numerical Simulation of Residual Stress Field in Strengthened Grinding of GCr15 Bearing Steel Based on Normal Distribution

a,b,c,a,b,c,a,b,c,a,b,c,a,b,c,a,b,c,a,b,c

(a. School of Mechanical & Electric Engineering, b. Guangzhou Key Laboratory for Strengthened Grinding and High Performance Machining of Metal Material, c. Guangdong Engineering and Technology Research Centre for Strengthen Grinding and High Performance Micro-nanomachining, Guangzhou University, Guangzhou 510006, China)

This paper aims to explore the influence of fixed point injection on residual stress field of GCr15 bearing steel under different process parameters of strengthening grinding. Image processing technology was used to analyze the distribution characteristics of the surface coverage along the width direction of the strengthening grinding fixed point injection under different process parameters. The two-dimensional normal distribution function was used to describe the distribution characteristics of steel ball coordinates under the fixed point injection of strengthening grinding. Python/Opencv was used to determine the number of steel balls required by the finite element model under different process parameters. Based on Abaqus/ Python, the finite element model with normal distribution of strengthening grinding was constructed. Using the established normal distribution model, the influence of different injection velocity, ball diameter and coverage rate on the residual stress field was analyzed. When the injection velocity increased from 45 m/s to 70 m/s, the surface residual compressive stress increased from ?683.5 MPa to ?902.4 MPa, the maximum residual compressive stress increased from ?981.6 MPa to ?1330.6 MPa, and the thickness of residual pressure layer increased from 89 μm to 151 μm. The maximum residual compressive stress depth was moved from 30 μm to 70 μm. When the diameter of the steel ball increased from 0.4 mm to 1.0 mm, the surface residual compressive stress increased first and then decreased, the maximum residual compressive stress increased from ?1063.5 MPa to ?1240.7 MPa, and the maximum residual compressive stress depth increased from 30 μm to 60 μm. The thickness of residual compressive stress layer increased from 103 μm to 147 μm, and the diameter of steel ball increased from 0.8 mm to 1.0 mm, and the maximum residual compressive stress almost remained unchanged. When the injection coverage rate was from 100% to 300%, the surface residual compressive stress, the maximum residual compressive stress and the maximum residual compressive stress depth increased slightly, and the residual compressive stress layer thickness almost remained unchanged. The simulation values of the normal distribution model and the random distribution model were compared with the experimental values, and it was found that the surface residual compressive stress, maximum residual compressive stress depth and residual compressive stress thickness of the three models were almost the same. The simulation values of the random distribution model and the normal distribution model were 32.1% and 18.9% higher than the experimental values of the maximum residual compressive stress. The finite element model with normal distribution can accurately predict the change process of residual stress, which can provide some guidance for the optimization of process parameters of strengthening grinding.

strengthened grinding; residual stress field; normal distribution; numerical simulation; GCr15 bearing steel

2021-05-18；

2021-08-18

FAN Li-wei (1997—), Male, Postgraduate, Research focus: surface hardening technology of metal materials.

梁忠偉（1978—），男，博士，教授，主要研究方向?yàn)橹悄苤圃旒把b備技術(shù)。

LIANG Zhong-wei (1978—), Male, Doctor, Professor, Research focus: intelligent manufact uring and equipment technology.

范立維，梁忠偉，劉曉初，等. 基于正態(tài)分布的GCr15軸承鋼強(qiáng)化研磨殘余應(yīng)力場(chǎng)數(shù)值模擬[J]. 表面技術(shù), 2022, 51(3): 242-253.

Tg580.68

A

1001-3660(2022)03-0242-12

10.16490/j.cnki.issn.1001-3660.2022.03.026

2021-05-18；

2021-08-18

國(guó)家自然科學(xué)基金項(xiàng)目（51975136，51575116)；廣東省科技計(jì)劃項(xiàng)目（2017A010102014，2016A010102022）；國(guó)家重點(diǎn)研發(fā)計(jì)劃（2018YFB2000501）；廣州大學(xué)全日制研究生基礎(chǔ)創(chuàng)新項(xiàng)目資助（2020GDJC-M18）；廣東省高等學(xué)?？萍紕?chuàng)新團(tuán)隊(duì)項(xiàng)目（2017KCXTD025）；廣州市教育系統(tǒng)創(chuàng)新學(xué)術(shù)團(tuán)隊(duì)項(xiàng)目（1201610013）；廣東省高校重點(diǎn)領(lǐng)域?qū)ｍ?xiàng)（2019KZDZX1009）；廣州市科技計(jì)劃項(xiàng)目（201707010293）

Fund：National Natural Science Foundation of China (51975136, 51575116); Guangdong Science and Technology Program Project (2017A010102014, 2016A010102022); National Key Research and Development Project (2018YFB2000501); Innovation Research for the Postgraduates of Guangzhou University under Grant (2020GDJC-M18); Guangdong Province University Science and Technology Innovation Team Project (2017KCXTD025); Guangzhou Education System Innovation Academic Team Project (1201610013); Special Projects in Key Fields of Universities in Guangdong Province (2019KZDZX1009); Guangzhou Science and Technology Program Project (201707010293)

范立維（1997—），男，碩士研究生，主要研究方向?yàn)榻饘俨牧媳砻鎻?qiáng)化技術(shù)。

FAN Li-wei, LIANG Zhong-wei, LIU Xiao-chu, et al. Numerical Simulation of Residual Stress Field in Strengthened Grinding of GCr15 Bearing Steel Based on Normal Distribution[J]. Surface Technology, 2022, 51(3): 242-253.