發布時間：2022-04-13所屬分類：工程師職稱論文瀏覽：1次
摘 要： 摘要: 利用 DEFORM 有限元軟件對某型號雙列調心滾子軸承外圈感應淬火進行模擬，采用響應曲面法對工藝參數進行優化，并根據優化后工藝參數得到深冷處理前后滾道中心和密封槽底淬硬層深度以及滾道中心次表面殘余壓應力。結果表明: 二次響應曲面能夠很好地擬合外圈感應淬
摘要: 利用 DEFORM 有限元軟件對某型號雙列調心滾子軸承外圈感應淬火進行模擬，采用響應曲面法對工藝參數進行優化，并根據優化后工藝參數得到深冷處理前后滾道中心和密封槽底淬硬層深度以及滾道中心次表面殘余壓應力。結果表明: 二次響應曲面能夠很好地擬合外圈感應淬火模擬結果，淬硬層深度各影響因素顯著程度由高到低依次為電流頻率、電流密度和線圈轉速，殘余壓應力各影響因素顯著程度由高到低依次為線圈轉速、電流密度和電流頻率; 最優工藝參數為線圈轉速 0. 049 r/min、電流頻率 246. 667 kHz 和電流密度 6. 5 × 107 A/m2 ; 深冷處理后，滾道中心淬硬層深度、密封槽底淬硬層深度、滾道中心次表面殘余壓應力以及滾道表面硬度分別為 1. 642 mm，0. 608 mm，120 MPa 和 62. 5 HRC; 深冷處理能夠有效增加淬硬層深度，降低殘余壓應力并提高淬火硬度; 數值模擬結果與試驗結果相比，誤差在 13% 以內，具有良好的一致性。
關鍵詞: 滾動軸承; 雙列調心滾子軸承; 高氮鋼; 感應淬火; 數值模擬; 淬硬層; 殘余應力; 深冷處理
0 引言
滾動軸承是飛機機體機械設備支承和維持轉動的關鍵部件，近年來，飛行速度的提升使得對飛機安全性與可靠性的要求越來越高，意味著對滾動軸承性能要求更加嚴格。感應淬火工藝是提高零件表面硬度、強度以及耐磨性的重要手段，是一個涉及電磁、傳熱、組織轉變和應力應變的多場相互作用過程，特別是受熱材料的物理性質隨其狀態和外加場的不同而劇烈變化，很難采用解析的方式精確描述。
文獻[1 - 2]針對多頻感應淬火，通過受溫度影響的電導率和焦耳加熱項對電磁和熱進行耦合，建立了包含電磁及熱力學效應的鋼感應淬火的數學模型，研究了相變引起的相變應變和相變塑性等，結果表明該淬火方式可以很好地控制復雜工件的淬硬層輪廓。文獻[3]采用有限元分析方法預測了感應淬火過程中的組織演變，并根據電磁 - 熱耦合有限元模型開發了用于計算感應淬火加熱和冷卻過程中組織轉變的專用程序，據此分析了二維軸對稱幾何模型的感應淬火過程，仿真與試驗結果有較好的一致性。文獻[4]采用響應曲面法優化了連續感應淬火后的相變區域，研究了進給速度、輸入功率、間隙、曲率和流速對溫度、顯微組織、顯微硬度和幾何形狀的影響，結果表明進給速度對連續感應淬火過程的組織均勻性、顯微硬度、氧化帶寬的影響較大，而工件曲率對其中心溫度的影響最大。文獻[5 - 7]首先基于電磁 - 熱耦合的數值模擬方法對圓柱體感應淬火過程進行了模擬，結果表明淬火前自由冷卻能夠降低零件表層溫度梯度，據此對淬火過程進行優化，避免奧氏體分解成非馬氏體微觀組織; 之后對直齒圓柱齒輪異步雙頻感應淬火進行數值模擬，結果表明該淬火方式能夠獲得均勻的硬化層輪廓; 最后考慮材料蠕變特性，對軸類零件感應淬火和回火過程進行數值模擬，并基于所提出的數值模型分析了相變和蠕變對回火應力松弛的影響。文獻[8]建立了 55CrMo 鋼精密滾珠絲杠單感應線圈加熱的有限元數值模擬模型，據此繪制了滾珠絲杠溝道感應加熱過程中的溫度變化曲線，并針對單感應線圈在加熱效率與深度上的缺陷提出雙感應線圈更適合滾珠絲杠感應淬火。文獻[9]利用 MSC. Marc 軟件，采用逐步循環加載方法建立了 55CrMo 鋼絲杠感應加熱過程的有限元數值模擬模型，研究了 5010 型絲杠感應加熱過程中工藝參數對材料奧氏體化情況的影響，并基于響應曲面法建立了工藝參數與淬火性能之間關系的回歸方程，獲得了最優工藝參數。文獻[10]建立了雙列球軸承內圈感應淬火過程耦合的數值模擬模型，并據此研究了淬火過程中各部位溫度、組織及應力隨時間的變化以及工藝參數對淬火性能的影響，結果表明，感應加熱的升溫速度隨著電流頻率和密度的增大而增大，淬硬層深度隨著電流頻率的增大逐漸減小，在殘余奧氏體較多的部位殘余應力較大，感應加熱時工件各部位的升溫速度與感應器形狀有關。文獻[11]基于電磁場、溫度場、組織場、應力應變場耦合模型，利用 DEFORM 有限元軟件模擬了雙列調心滾子軸承內圈感應淬火過程，提出了分段電流密度的淬火工藝，研究了內圈感應淬火過程中溫度變化、組織演變、表面和次表面硬度、殘余應力和殘余奧氏體等，結果表明感應淬火電流密度分段設置適用于小型復雜零件的表面淬火，能夠得到均勻的淬火組織。
綜上所述，國內外學者雖然從耦合數學模型和有限元模擬等方面，并結合試驗在感應淬火過程各物理場建模和淬火后工件物理性能變化等方面開展了大量研究，但到目前為止對軸承等特殊結構零件特別是高氮鋼材料的感應淬火工藝研究卻少見報道。鑒于此，本文利用 DEFORM 有限元分析軟件，建立感應淬火電磁場、溫度場、組織場及應力應變場耦合的有限元數值模擬模型，研究某型號雙列調心滾子軸承外圈感應淬火工藝參數對淬火性能的影響，并據此對參數進行優化。
1 建立感應淬火過程數值模擬模型
1. 1 感應淬火過程耦合物理場
感應淬火過程包括電磁、溫度、組織及應力應變 4 個物理場，其耦合相互作用關系如圖 1 所示，圖中: 1 為材料電磁特性，2 為焦耳熱，3 為溫度變化引起的相變，4 為相變潛熱，5 為相變引起的相變應力和相變塑性，6 為應力應變誘導相變，7 為機械功生熱，8 為熱應力，9 和 10 為熱交換。試驗材料的溫度分布受焦耳熱、相變潛熱以及變形生熱共同作用的影響，材料自身性能隨溫度的變化反過來又影響電磁場能量的變化，從而影響加熱過程，且溫度變化引起材料的相變，隨之產生相變應力，同時工件的變形也影響溫度，從而引起溫度場的變化。
1. 2 計算模型
以某型號雙列調心滾子軸承外圈為例，其結構參數見表 1，由于需對外圈密封槽底淬硬層深度進行研究，故不能去除滾道部位倒角。由于外圈主要淬火區域為滾道，且該外圈滾道尺寸較大，結構簡單，故采用仿形感應淬火法。因此，外圈感應線圈設計如圖 2 所示，在感應加熱過程中線圈繞外圈軸線轉動，線圈徑向截面的長度為 4 mm、寬度為 3 mm 和 45°倒角為 1 mm。由于滾道兩側存在密封槽底，為達到淬火性能的要求，感應線圈在滾道處與淬火表面的距離為 1. 0 mm，滾道兩側與淬火表面距離為 0. 5 mm。
采用 DEFORM 有限元分析軟件中的邊界元法對雙列調心滾子軸承外圈滾道表面進行感應淬火分析。考慮外圈為回轉體，對外圈進行簡化以提高計算效率，取外圈 1 /12 為研究對象建立三維有限元計算模型，對其進行四面體網格劃分，如圖 3a 所示�？紤]感應淬火集膚效應，對外圈滾道表面感應淬火區網格進行細化，總單元數為 96 568，節點數為 21 462，感應線圈單元數為 15 694，節點數為 3 704，并對外圈截面設置為對稱約束，端面設置為固定約束。為方便后續分析，在外圈密封槽底的徑向方向及滾道中心的徑向方向標記 Ps1， Pe1，Ps2，Pe2特征節點 ( 圖 3b) ，分別表征外圈密封槽底和滾道中心位置感應淬火特征的梯度變化及淬硬層深度( 馬氏體體積分數為 50% ) 的測量點。
1. 3 感應淬火工藝參數
軸承外圈材料為 X30CrMoN15 - 1 高氮鋼，感應淬火前外圈的初始組織為珠光體。針對外圈仿形感應淬火特征，線圈繞外圈軸線轉動。根據多次重復模擬，線圈轉速 n、電流頻率 f 和電流密度 J 取值見表 2。另外，假設環境溫度為 20 ℃，感應加熱溫度為 1 000 ℃，淬火液跟隨線圈轉動將外圈冷卻至室溫。
2 響應曲面法模擬結果與分析
2. 1 響應曲面法模擬方案及結果
根據表 2 工藝參數的取值，進行組合模擬，并選用響應曲面法對模擬結果進行分析。根據本文研究對象及響應曲面法試驗設計要求，利用 Design Expert 軟件進行 BBD 模擬[12]，模擬方案及結果見表 3。
2. 3 感應淬火工藝參數優化
以感應淬火后外圈滾道中心位置和密封槽底的淬硬層深度及滾道中心次表面殘余壓應力為優化目標，對分析結果進行多目標有約束優化。為保證感應淬火后淬硬層深度合理，要求滾道中心位置淬硬層深度取值為 0. 950 ～ 1. 600 mm，密封槽底淬硬層深度取值為 0. 300 ～ 0. 650 mm，建立如( 5) 式所示的優化目標函數和約束函數。
實際生產中，在合理的感應淬火范圍內，線圈轉速直接影響生產效率，轉速越高，生產效率越高，且電流密度越小，越節約能源，所以為降低加工成本，應選用較高的線圈轉速與較小的電流密度。綜合考慮，選取表 7 中第 4 組工藝參數為最終優化結果，在 2. 4 節中采用該組參數進行感應淬火數值模擬，研究淬火后深冷處理( - 196 ℃ × 4 h) 對外圈性能的影響。
2. 4 深冷處理對外圈淬火性能的影響
深冷處理是在鋼試樣冷卻到室溫后，繼續在 0 ℃以下的介質中冷卻的熱處理工藝，能夠最大限度地降低淬火后鋼試樣的殘余奧氏體，進一步提高熱處理性能。
滾道中心和密封槽底處深冷處理前后殘余奧氏體在深度方向的梯度分布如圖 4 所示，殘余奧氏體含量( 體積分數) 隨著深度的增加逐漸下降，淬火后表面殘余奧氏體含量分別為 8. 53% 和 8. 52%，深冷處理后分別降低為 1. 43% 和 1. 42% 。
滾道中心和密封槽底處深冷處理前后馬氏體在深度方向的梯度分布如圖 5 所示，淬火后表面馬氏體含量( 體積分數) 分別為 89. 61% 和 89. 59% ，深冷處理后分別提高為 95. 91% 和95. 89% 。將馬氏體含量超過 50% 的區域視為淬硬層，則滾道中心和密封槽底處淬硬層深度分別為 1. 609，0. 595 mm，深冷處理后分別為 1. 642，0. 608 mm。這是由于在深冷處理過程中，一部分殘余奧氏體將繼續轉化為馬氏體，使外圈馬氏體含量進一步提高，馬氏體含量超過 50% 的區域隨之增大，導致淬硬層深度增加。
本文來源于：《軸承》(月刊)創刊于1958年，是中國唯一國內外發行的軸承科技及其應用期刊，涵蓋了軸承設計、試驗及應用等部門�？d論文堅持技術超前和實用相結合。設有：產品設計與應用、工藝與裝備、試驗與分析、試驗分析、專題綜述、中國軸承論壇、國外軸承科技等欄目。
滾道中心和密封槽底處深冷處理前后硬度在深度方向的梯度分布圖 6 所示，淬火后表面硬度分別為 60. 4，60. 3 HRC，深冷處理后分別為 62. 5， 62. 3 HRC。由于深冷處理后馬氏體含量有所提高，而馬氏體含量是決定材料硬度的重要因素，所以外圈表面硬度在深冷處理后也略微增大。
滾道中心和密封槽底處深冷處理前后殘余壓應力在深度方向的分布如圖 7 所示，淬火后表面殘余壓應力分別為 78，76 MPa，深冷處理后分別為 54，59 MPa。最大殘余壓應力出現在次表面，淬火后最大殘余壓應力分別為 146，199 MPa，深冷處理后最大殘余壓應力分別為 120，172 MPa。深冷過程中外圈內少量殘余奧氏體繼續轉化為馬氏體，組織的體積增大產生組織應力，會使其內部殘余壓應力略微上升; 但深冷處理還會使材料組織細化，位錯密度增加，內部缺陷減少，這些又會降低殘余壓應力; 綜合作用下，深冷處理后外圈殘余壓應力呈降低趨勢。在密封槽底處殘余壓應力在次表面出現 2 個峰值( 圖 7b 和圖 8) ，這是由于該位置厚度較大，而滾道部位厚度較小，冷卻過程中受滾道部位外表面影響較大，故該位置第二次表面出現小應力峰值，深冷處理前后分別為151，134 MPa。
2. 5 小結
深冷處理前，根據 2. 3 節優化的工藝參數，使用二次響應曲面擬合公式預測的滾道中心淬硬層深度、密封槽底淬硬層深度和滾道中心次表面殘余壓應力分別為 1. 547 mm，0. 614 mm，138 MPa，模擬值分別為 1. 609 mm，0. 595 mm，146 MPa，預測值與模擬值的誤差在 6% 以內。
3 試驗驗證
使用第 2. 3 節優化后的工藝參數，采用圖 9 所示感應淬火設備對 X30CrMoN15 -1 鋼制雙列調心滾子軸承外圈進行感應淬火，設備包括控制系統、工件裝夾移動裝置、感應線圈及噴液冷卻裝置。為驗證軸承外圈感應淬火數值模擬的正確性，對感應淬火試驗后滾道表面殘余應力、殘余奧氏體含量、硬度及淬硬層深度進行測量，試驗用軸承外圈數量為 3 個。在深冷處理后的外圈滾道表面中心圓周上( 即 Ps2所在滾道位置) 取 3 個相距 120°的測量點。使用 IXRD 殘余應力 X 射線衍射儀，測量其殘余應力和殘余奧氏體含量，試驗值與模擬值對比結果分別見表8 和表9，誤差在13%以內。
通過線切割將軸承外圈沿軸向方向分割為 3 個試樣，使用 HR - 150DT 型電動洛氏硬度計對熱處理后外圈試樣的滾道位置的硬度進行測量，取 3 個試樣的平均硬度值作為該外圈滾道位置的硬度，試驗值與模擬值對比結果見表 10，誤差在 8% 以內。
取第 2. 3 節所得滾道 50% 馬氏體處對應的硬度值作為測量參考，試驗測量硬度值在參考硬度值以上區域為淬硬層。在硬度測量用的 3 個外圈上繼續進行淬硬層深度測量，測量位置為外圈滾道中心，取 3 個外圈淬硬層深度平均值作為該滾道測量位置的淬硬層深度值，試驗值與模擬值對比結果見表 11，誤差在 6% 以內。
4 結論
針對雙列調心滾子軸承外圈的仿形感應淬火過程進行了數值模擬，通過響應曲面法獲得了最優工藝參數，據此分析了深冷處理前后外圈性能的變化，得到如下結論:
1) 采用二次響應曲面對外圈感應淬火模擬結果進行回歸分析，擬合度較好。淬硬層深度各影響因素顯著程度由高到低依次為電流頻率、電流密度、線圈轉速，殘余壓應力各影響因素顯著程度由高到低依次為線圈轉速、電流密度、電流頻率。最優工藝參數為線圈轉速 0. 049 r/min、電流頻率 246. 667 kHz、電流密度 6. 5 × 107 A/m2 。
2) 深冷處理前，使用二次響應曲面擬合的公式預測滾道中心淬硬層深度、密封槽底淬硬層深度和滾道中心次表面殘余壓應力值分別為 1. 547 mm，0. 614 mm，138 MPa，模擬值分別為 1. 609 mm， 0. 595 mm，146 MPa，預測值與模擬值誤差在 6% 以內。采用最優工藝參數進行模擬，得到深冷處理前滾道表面殘余奧氏體含量約為 8. 53% ，表面硬度約為 60. 4 HRC。
3) 深冷處理后，滾道中心淬硬層深度、密封槽底淬硬層深度和滾道中心次表面殘余壓應力值分別為 1. 642 mm，0. 608 mm，120 MPa，滾道表面殘余奧氏體含量約為 1. 43% ，表面硬度約為 62. 5 HRC。深冷處理能夠有效降低殘余壓應力，增加淬硬層深度，且能夠降低殘余奧氏體含量，提高馬氏體含量，提高淬火硬度。
模擬結果與試驗結果誤差在 13% 以內，證明本文研究成果的可靠性與正確性，能夠為該型號雙列調心滾子軸承外圈感應淬火工藝提供理論依據。——論文作者：楊虎1，2，3 ，張文虎4，5 ，謝鵬飛1，2，3 ，鄭金濤4 ，張松5
參考文獻:
[1] HMBERG D. A mathematical model for induction hardening including mechanical effects[J]. Nonlinear Analysis: Real World Applications，2004，5( 1) : 55 - 90.
[2] HMBERG D，LIU Q Z，MONTALVO - URQUIZO J，et al. Simulation of multi-frequency-induction-hardening including phase transitions and mechanical effects[J]. Finite Elements in Analysis and Design，2016，121: 86 - 100.
[3] SPEZZAPRIA M，FORZAN M，DUGHIERO F. Numerical simulation of solid-solid phase transformations during induction hardening process[J]. IEEE Transactions on Magnetics，2016，52( 3) : 1 - 4.
[4] QIN X P，GAO K，ZHU Z H，et al. Prediction and optimization of phase transformation region after spot continual induction hardening process using response surface method[J]. Journal of Materials Engineering and Performance，2017，26( 9) : 4578 - 4594.
[5] TONG D M，GU J F，TOTTEN G E. Numerical simulation of induction hardening of a cylindrical part based on multi-physics coupling[J]. Modelling and Simulation in Materials Science and Engineering，2017，25( 3) : 1 - 20.
[6] TONG D M，GU J F，TOTTEN G E. Numerical investigation of asynchronous dual - frequency induction hardening of spur gear[J]. International Journal of Mechanical Sciences，2018，142 /143: 1 - 9.
[7] TONG D M，GU J F，YANG F. Numerical simulation on induction heat treatment process of a shaft part: involving induction hardening and tempering[J]. Journal of Materials Processing Technology，2018，262: 277 - 289.
[8] 賀連芳，李輝平，蓋康，等. 55CrMo 鋼感應淬火工藝的數值模擬及工藝優化[J]. 材料熱處理學報，2015， 36( 1) : 199 - 204.
[9] 蓋康，賀連芳，張春芝，等. 基于 RSM 的絲杠感應淬火工藝數值模擬及參數優化[J]. 材料熱處理學報， 2016，37( 增刊 1) : 146 - 152.
[10] 臧樂航，鄧四二，張文虎. 雙列球軸承中頻感應淬火數值模擬[J]. 材料熱處理學報，2018，39( 1) : 137 - 144.
SCISSCIAHCI
