曲軸表面強(qiáng)化工藝仿真研究(北理工-廖日東)

1、曲軸表面強(qiáng)化工藝仿真研究曲軸表面強(qiáng)化工藝仿真研究Numerical Simulation of the Crankshaft Surface Strengthening 廖日東廖日東Liao Ri-dong北京理工大學(xué)北京理工大學(xué) 動力系統(tǒng)工程研究所動力系統(tǒng)工程研究所Research Center of Power Machinery, Beijing Institute of Technology 曲軸及其常見的失效形式曲軸及其常見的失效形式Crankshaft and Common Failure Forms 軸頸異常軸頸異常磨損磨損Abnormal Wear of Crankshaft

2、Neck扭轉(zhuǎn)扭轉(zhuǎn)疲勞斷裂疲勞斷裂Torsional Fatigue Fracture 彎曲彎曲疲勞斷裂疲勞斷裂Bending Fatigue Fracture2常見的曲軸表面強(qiáng)化工藝常見的曲軸表面強(qiáng)化工藝Common Crankshaft Surface Strengthening Methods感應(yīng)淬火感應(yīng)淬火Induction Quenching圓角滾壓圓角滾壓Fillet Rolling氣體滲氮氣體滲氮Gas Nitriding34不考慮強(qiáng)化的曲軸強(qiáng)度計算不考慮強(qiáng)化的曲軸強(qiáng)度計算Crankshaft Strength Calculation without considering Sur

3、face Strengthening5不考慮強(qiáng)化的曲軸強(qiáng)度計算不考慮強(qiáng)化的曲軸強(qiáng)度計算Crankshaft Strength Calculation without considering Surface Strengthening表面強(qiáng)化工藝仿真研究表面強(qiáng)化工藝仿真研究Studies on Surface Strengthening Simulations 從工藝參數(shù)設(shè)計出發(fā)，預(yù)測工件的彈塑性力學(xué)性從工藝參數(shù)設(shè)計出發(fā)，預(yù)測工件的彈塑性力學(xué)性能、斷裂韌性和疲勞壽命等宏觀力學(xué)量，也就是說能、斷裂韌性和疲勞壽命等宏觀力學(xué)量，也就是說是實現(xiàn)強(qiáng)化工藝的預(yù)測設(shè)計。是實現(xiàn)強(qiáng)化工藝的預(yù)測設(shè)計。涉及材料學(xué)、工

4、藝學(xué)、熱力學(xué)、固體物理、電磁涉及材料學(xué)、工藝學(xué)、熱力學(xué)、固體物理、電磁學(xué)、彈塑性力學(xué)、斷裂力學(xué)、流體力學(xué)等。學(xué)、彈塑性力學(xué)、斷裂力學(xué)、流體力學(xué)等。6圓角滾壓數(shù)值分析圓角滾壓數(shù)值分析Numerical Analysis on Fillet Rolling 78圓角滾壓數(shù)值分析圓角滾壓數(shù)值分析Numerical Analysis on Fillet Rolling 圓角滾壓數(shù)值分析圓角滾壓數(shù)值分析Numerical Analysis on Fillet Rolling 9感感應(yīng)應(yīng)加加熱熱淬淬火火仿仿真真研研究究感應(yīng)加熱過程仿真感應(yīng)加熱過程仿真強(qiáng)化效果分析與驗證強(qiáng)化效果分析與驗證淬火冷卻過程仿真淬火

5、冷卻過程仿真工藝設(shè)計仿真分析試驗驗證感應(yīng)器結(jié)構(gòu)電工學(xué)參數(shù)冷卻液參數(shù)冷卻時間強(qiáng)化要求工件熱物性：電磁性能導(dǎo)/換熱性能力學(xué)性能電磁場分析溫度場分析相變分析應(yīng)力/位移場分析表面硬度測量淬裂缺陷檢測形位誤差檢測金相組織分析疲勞強(qiáng)度試驗殘余應(yīng)力測量材料性能分析相變應(yīng)力/塑性，與組織相關(guān)的力學(xué)性質(zhì)應(yīng)力/應(yīng)變誘導(dǎo)相變熱應(yīng)變/熱機(jī)械性能熱物性熱驅(qū)動力相變潛熱熱物性與相體積相關(guān)的電磁性質(zhì)塑性功生熱感應(yīng)渦流生熱熱流場分析溫度場分析相變分析應(yīng)力/位移場分析熱應(yīng)變/熱機(jī)械性能熱物性熱驅(qū)動力相變潛熱熱物性塑性功生熱沸騰/傳導(dǎo)換熱相變應(yīng)力/塑性，與組織相關(guān)的力學(xué)性質(zhì)應(yīng)力/應(yīng)變誘導(dǎo)相變10Simulation of Ind

6、uction Heating Quenching感應(yīng)感應(yīng)加熱加熱淬火淬火仿真仿真研究研究Simulation of Induction Heating Quenching11感應(yīng)加熱淬火數(shù)值分析感應(yīng)加熱淬火數(shù)值分析Numerical Analysis on Induction Heating Quenching12感應(yīng)加熱淬火數(shù)值分析感應(yīng)加熱淬火數(shù)值分析Numerical Analysis on Induction Heating Quenching13感應(yīng)加熱淬火數(shù)值分析感應(yīng)加熱淬火數(shù)值分析Numerical Analysis on Induction Heating Quenching14

7、感應(yīng)加熱淬火數(shù)值分析感應(yīng)加熱淬火數(shù)值分析Numerical Analysis on Induction Heating Quenching1542CrMo臨界冷卻速度感應(yīng)加熱淬火數(shù)值分析感應(yīng)加熱淬火數(shù)值分析Numerical Analysis on Induction Heating Quenching16典型位置各相淬火后硬度典型位置淬火后硬度感應(yīng)淬火感應(yīng)淬火仿真仿真研究研究Simulation of Induction Heating Quenching 材料的熱物性不可靠；沒有考慮加熱和冷卻速率（熱流密度變化率）的影響；在線的測量和驗證很困難，工藝的預(yù)測設(shè)計還沒有實現(xiàn)。17感應(yīng)加熱淬火感

8、應(yīng)加熱淬火仿真仿真研究研究Simulation of Induction Heating Quenching qTtqttTttTc)(2022018感應(yīng)加熱淬火仿真研究感應(yīng)加熱淬火仿真研究Simulation of Induction Heating Quenching 含瞬變熱源的強(qiáng)瞬態(tài)熱傳導(dǎo)方程的求解；含瞬變熱源的強(qiáng)瞬態(tài)熱傳導(dǎo)方程的求解；電磁感應(yīng)淬火過程中電磁感應(yīng)淬火過程中“熱波熱波”的傳播機(jī)制；的傳播機(jī)制；非非Fourier效應(yīng)對電磁感應(yīng)淬火中多場耦合動力學(xué)響應(yīng)的影響。效應(yīng)對電磁感應(yīng)淬火中多場耦合動力學(xué)響應(yīng)的影響。 電磁感應(yīng)淬火過程的瞬態(tài)熱傳導(dǎo)瞬態(tài)熱傳導(dǎo)的非Fourier效應(yīng)研究 電磁

9、感應(yīng)淬火過程中材料相變材料相變的非Fourier效應(yīng)研究 電磁感應(yīng)淬火過程應(yīng)力和變形應(yīng)力和變形的非Fourier效應(yīng)研究 電磁感應(yīng)淬火過程仿真結(jié)果試驗驗證仿真結(jié)果試驗驗證研究19xnmxmeebJebJq2211222222222)(21)(21含與不含空氣區(qū)的電磁感應(yīng)加熱生熱率解析解含與不含空氣區(qū)的電磁感應(yīng)加熱生熱率解析解Analytical Solutions of the Heat Generation Rate of the Electromagnetic Induction Heating 20感應(yīng)加熱淬火感應(yīng)加熱淬火仿真仿真研究研究Simulation of Induction H

10、eating Quenching 溫度場控制方程 基于CV模型 格林函數(shù)212200220000( , )( , )000ttxx LTTTQ x tQ x tcttxtttTTTtTTtt222211440022222212214412-e18ecossin1cos11e114sinsin2sin2,tlmmmmmmmmllmlmlmlmla Q caea Qc Qc Qmacxxa 22200.54 10C/s2mpJ bck 63=7.841 10 kg/m90.47 10 J/kg. Cpc 111=1.3 10 m( )xpQ xkce1111()( )0Q xxxxQ xxx341

11、3190 J/ m.577sQ 1111.7 10 xm感應(yīng)加熱淬火瞬態(tài)傳熱非感應(yīng)加熱淬火瞬態(tài)傳熱非FourierFourier效應(yīng)效應(yīng)Non-fourier Effects on Temperature Field of Induction Heating Quenching 脈沖寬度不同時 熱流量不變內(nèi)熱源形式脈沖寬度不同時 不同位置處溫度隨時間的變化22隨著脈沖寬度的增加,溫度波的幅值減小,波動減弱。說明熱流傳播的波動性只在熱作用時間短到與材料熱弛豫時間相當(dāng)或更短的尺度才對傳熱過程起主導(dǎo)作用。0.1,0.5,1,5p0.1,0.5,1,5p感應(yīng)加熱淬火瞬態(tài)傳熱非感應(yīng)加熱淬火瞬態(tài)傳熱非Fo

12、urierFourier效應(yīng)效應(yīng)Non-Fourier Effects on Temperature Field of Induction Heating Quenching 氣體滲碳仿真研究氣體滲碳仿真研究Simulation of Gas CarburizingqTtTc202CDtC)(TD),(tzyxT碳的擴(kuò)散過程Carbon Diffusion Process(Fick第二定律)：(Fichs Second Law)工件加熱過程Workpiece Heating Process(瞬態(tài)熱傳導(dǎo)方程)：(Transient Heat Conduction Equation)23總結(jié)與展望

13、總結(jié)與展望Summaries and Prospects從Maxwell基本方程出發(fā)，推導(dǎo)了一維電磁感應(yīng)加熱過程生熱率?；陔姶鸥袘?yīng)加熱淬火過程由電磁感應(yīng)產(chǎn)生的指數(shù)型內(nèi)熱源，簡化為線性內(nèi)熱源，得到其非傅里葉導(dǎo)熱溫度場。建立了電磁感應(yīng)加熱有限元模型，利用梅尼爾模型，預(yù)測了水淬后的淬硬層和硬度。水淬后的硬度滿足設(shè)計要求，但是硬度高于設(shè)計值，需要采用比水冷卻性能弱的淬火液。24對于考慮材料非線性下，其內(nèi)熱源形式，得到相應(yīng)的非Fourier效應(yīng)的溫度場。對于電磁感應(yīng)加熱淬火熱傳導(dǎo)的非Fourier效應(yīng)，可以進(jìn)一步研究球體下的溫度場分布。研究不同非Fourier熱傳導(dǎo)數(shù)學(xué)模型對溫度場結(jié)果的影響。研究電磁

14、感應(yīng)加熱淬火過程中的非Fourier效應(yīng)25總結(jié)與展望總結(jié)與展望Summaries and Prospects氣體滲碳仿真研究氣體滲碳仿真研究Simulation of Gas Carburizingu滲碳問題有限元計算收斂性分析 由于滲碳層很薄，所以工件表層的網(wǎng)格一般需要加密處理,但是工件表層的網(wǎng)格需要加密到什么程度才能滿足精度要求？u大型工件碳濃度場模擬簡化計算 現(xiàn)有文獻(xiàn)都是針對非常簡單的結(jié)構(gòu)進(jìn)行碳濃度場的模擬，如何比較簡單的模擬大型工件的碳濃度場仍是一項尚待解決的問題。u溫度空間分布對碳濃度場的影響 為了準(zhǔn)確模擬碳濃度場，也為了對滲碳爐內(nèi)工件的布置進(jìn)行指導(dǎo)，考慮溫度場空間分布對碳濃度場的

15、影響是十分必要的。26滲碳問題有限元計算收斂性分析滲碳問題有限元計算收斂性分析Convergence Analysis of Finite Element Carburizing Model1.一維滲碳問題 Laplace變換法0.00.10.20.30.40.50.60.70.80.91.00.00.20.40.60.81.01.25000s4000s3000s2000s1000sC/%h/mm理論計算的缺點：1.只能求解一維問題2.對滲碳時間有限制理論求解：低碳鋼初始碳濃度一般為0.2%左右，經(jīng)過滲碳，工件表層碳含量發(fā)生變化，本文將工件表面至碳濃度變化量為0.01%（初始碳濃度的5%）處定

16、義為滲碳影響區(qū)域，該區(qū)域深度定義為DH 。t /s100015002000250030003500400045005000DH /mm0.490.610.710.790.870.941.011.071.13不同時間對應(yīng)的DH值DH值估算式27滲碳問題有限元計算收斂性分析滲碳問題有限元計算收斂性分析Convergence Analysis of Finite Element Carburizing Model數(shù)值模擬：表征滲碳區(qū)域內(nèi)節(jié)點碳濃度精確度網(wǎng)格尺寸/mm0.10.20.30.40.51000s0.006 0.009 0.017 0.028 0.025 1500s0.004 0.007 0

17、.012 0.021 0.028 2000s0.003 0.006 0.010 0.018 0.026 2500s0.003 0.005 0.008 0.015 0.019 3000s0.002 0.004 0.007 0.013 0.017 3500s0.002 0.004 0.007 0.012 0.015 4000s0.002 0.003 0.006 0.011 0.013 4500s0.002 0.003 0.006 0.010 0.012 5000s0.002 0.003 0.005 0.009 0.011 一維滲碳不同滲碳時間及網(wǎng)格尺寸條件下的值0.000 0.005 0.010

18、0.015 0.020 0.025 0.030 100020003000400050000.1mm0.2mm0.3mm0.4mm0.5mmt/s0.000 0.005 0.010 0.015 0.020 0.025 0.030 0.10.20.30.40.51000s1500s2000s2500s3000s3500s4000s4500s5000s網(wǎng)格大小/mm網(wǎng)格尺寸對節(jié)點碳濃度的影響28滲碳問題有限元計算收斂性分析滲碳問題有限元計算收斂性分析Convergence Analysis of Finite Element Carburizing Modelt /sDh(理論計算) /mmDh(模

19、擬計算) /mm0.10.20.30.40.510000.30 0.30 0.35 0.38 0.39 0.47 15000.37 0.38 0.39 0.46 0.47 0.48 20000.43 0.43 0.46 0.51 0.56 0.51 25000.48 0.49 0.52 0.54 0.61 0.61 30000.53 0.53 0.55 0.57 0.65 0.68 35000.57 0.58 0.59 0.60 0.69 0.73 40000.61 0.61 0.63 0.65 0.71 0.77 45000.65 0.66 0.67 0.70 0.74 0.80 50000

20、.69 0.69 0.71 0.74 0.76 0.83 一維滲碳Dh的理論值與模擬計算結(jié)果0.00.10.20.30.40.50.60.70.80.910002000300040005000理論值0.1mm0.2mm0.3mm0.4mm0.5mmDh/mmt/s網(wǎng)格尺寸對滲碳層深度的影響滲碳層深度Dh定義為碳濃度為0.3%處的深度 虛線表示Dh值與理論值相差0.1mm為得到較為精確的節(jié)點碳濃度值，網(wǎng)格尺寸只需小于Dh/329滲碳問題有限元計算收斂性分析滲碳問題有限元計算收斂性分析Convergence Analysis of Finite Element Carburizing Model數(shù)

21、值模擬：以0.1mm網(wǎng)格計算結(jié)果為標(biāo)準(zhǔn)二維滲碳不同滲碳時間及網(wǎng)格尺寸條件下的值網(wǎng)格尺寸對節(jié)點碳濃度的影響2.二維尖角滲碳網(wǎng)格尺寸 /mm0.20.30.40.50.62h0.002 0.007 0.011 0.013 0.011 3h0.002 0.005 0.009 0.013 0.016 4h0.001 0.003 0.007 0.011 0.015 5h0.001 0.003 0.005 0.009 0.013 6h0.001 0.002 0.005 0.007 0.011 7h0.001 0.002 0.004 0.006 0.009 8h0.001 0.002 0.003 0.005

22、 0.008 0.000 0.002 0.004 0.006 0.008 0.010 0.012 0.014 0.016 0.018 23456780.2mm0.3mm0.4mm0.5mm0.6mmt/h0.000 0.002 0.004 0.006 0.008 0.010 0.012 0.014 0.016 0.018 0.20.30.40.50.62h3h4h5h6h7h8h網(wǎng)格尺寸/mm30滲碳問題有限元計算收斂性分析滲碳問題有限元計算收斂性分析Convergence Analysis of Finite Element Carburizing Model網(wǎng)格尺寸對滲碳層深度的影響虛線表

23、示Dh值與0.1mm網(wǎng)格計算值相差0.1mm網(wǎng)格尺寸/mm0.10.20.30.40.50.6Dh /mm2h1.68 1.69 1.75 1.80 1.82 1.87 3h2.05 2.07 2.10 2.14 2.21 2.25 4h2.36 2.38 2.40 2.45 2.43 2.58 5h2.64 2.65 2.69 2.69 2.77 2.80 6h2.89 2.90 2.91 2.97 3.02 2.96 7h3.12 3.14 3.16 3.18 3.20 3.27 8h3.34 3.35 3.37 3.38 3.41 3.50 二維滲碳不同計算條件下的Dh值1.52.02.

24、53.03.523456780.1mm0.2mm0.3mm0.4mm0.5mm0.6mmDh/mmt/h為得到較為精確的節(jié)點碳濃度值，網(wǎng)格尺寸只需小于Dh/431數(shù)值模擬：仍以0.1mm網(wǎng)格計算結(jié)果為標(biāo)準(zhǔn)三維滲碳不同滲碳時間及網(wǎng)格尺寸條件下的值網(wǎng)格尺寸對節(jié)點碳濃度的影響3.三維尖角滲碳網(wǎng)格尺寸/mm0.20.30.40.50.62h0.003 0.008 0.014 0.018 0.016 3h0.002 0.005 0.010 0.016 0.021 4h0.002 0.004 0.008 0.012 0.018 5h0.001 0.003 0.006 0.010 0.015 6h0.001

25、 0.003 0.005 0.009 0.012 7h0.001 0.003 0.005 0.007 0.011 8h0.001 0.002 0.004 0.007 0.009 0.000 0.005 0.010 0.015 0.020 0.025 23456780.2mm0.3mm0.4mm0.5mm0.6mmt/h0.000 0.005 0.010 0.015 0.020 0.025 0.20.30.40.50.62h3h4h5h6h7h8ht/h32滲碳問題有限元計算收斂性分析滲碳問題有限元計算收斂性分析Convergence Analysis of Finite Element Car

26、burizing Model滲碳問題有限元計算收斂性分析滲碳問題有限元計算收斂性分析Convergence Analysis of Finite Element Carburizing Model網(wǎng)格尺寸對滲碳層深度的影響虛線表示Dh值與0.1mm網(wǎng)格計算值相差0.1mm三維滲碳不同計算條件下的Dh值網(wǎng)格尺寸/mm0.10.20.30.40.50.6Dh /mm2h1.94 1.98 2.01 2.04 2.18 2.12 3h2.37 2.38 2.44 2.51 2.51 2.67 4h2.72 2.74 2.79 2.79 2.90 2.95 5h3.04 3.06 3.08 3.16

27、3.21 3.19 6h3.32 3.35 3.38 3.40 3.42 3.56 7h3.59 3.61 3.62 3.68 3.73 3.81 8h3.83 3.85 3.89 3.93 3.98 4.00 1.5 2.0 2.5 3.0 3.5 4.0 23456780.1mm0.2mm0.3mm0.4mm0.5mm0.6mmDh/mmt/h為得到較為精確的節(jié)點碳濃度值，網(wǎng)格尺寸只需小于Dh/533大型工件碳濃度場模擬簡化計算大型工件碳濃度場模擬簡化計算Simplified Calculation of the Large Work Piece Carbon Concentration

28、Field大型工件結(jié)構(gòu)特征1結(jié)構(gòu)特征2一維滲碳對于大型工件滲碳模擬，在計算碳濃度場、金相組織場以及表面硬度時，可以將大型工件拆分成多個結(jié)構(gòu)特征分別進(jìn)行計算341.平面尖角滲碳ABFCED 0.00.20.40.60.81.01.2012345對角線一維C/%h/mm在遠(yuǎn)離尖角部位碳濃度分布呈現(xiàn)一維特性/608090100120140Lc/mm4.242.922.452.111.661.30平面外尖角不同尖角角度對應(yīng)的Lc值平面外尖角0.01.02.03.04.05.06080100120140Lc/mm/L定義為EF上某點到F點的距離35大型工件碳濃度場模擬簡化計算大型工件碳濃度場模擬簡化計算

29、Simplified Calculation of the Large Work piece Carbon Concentration Field0.00.20.40.60.81.01.201234對角線一維C/%h/mm/608090100120140Lc/mm1.061.021.000.980.940.84平面內(nèi)尖角不同尖角角度對應(yīng)的Lc值平面內(nèi)尖角1.平面尖角滲碳ABFCED在遠(yuǎn)離尖角部位碳濃度分布呈現(xiàn)一維特性0.60.70.80.91.01.16080100120140Lc/mm/L定義為CD上某點到點的距離36大型工件碳濃度場模擬簡化計算大型工件碳濃度場模擬簡化計算Simplifie

30、d Calculation of the Large Work piece Carbon Concentration Field平面外圓角2.平面圓角滲碳ABDECF2F1F 0.00.20.40.60.81.01.20123對角線一維C/%h/mm與平面外尖角相比，倒圓結(jié)構(gòu)有利于碳濃度分布更加均勻L定義為EF1上某點到F點的距離37大型工件碳濃度場模擬簡化計算大型工件碳濃度場模擬簡化計算Simplified Calculation of the Large Work piece Carbon Concentration Field/r/mm1234604.24 4.39 5.737.2480

31、2.933.164.075.14902.46 2.753.504.351002.14 2.453.003.691201.701.872.242.631401.31 1.391.561.79平面外圓角2.平面圓角滲碳平面外圓角不同尖角角度及倒圓半徑對應(yīng)的Lc值0.0 2.0 4.0 6.0 8.0 1234608090100120140Lc/mmr/mm0.0 2.0 4.0 6.0 8.0 6080100120140r=1r=2r=3r=4Lc/mm/38大型工件碳濃度場模擬簡化計算大型工件碳濃度場模擬簡化計算Simplified Calculation of the Large Work p

32、iece Carbon Concentration Field平面內(nèi)圓角2.平面圓角滲碳與平面內(nèi)尖角相比，倒圓結(jié)構(gòu)有利于碳濃度分布更加均勻 BACC1C2DEF0.00.20.40.60.81.01.20123對角線一維C/%h/mmL定義為C1B上某點到點的距離39大型工件碳濃度場模擬簡化計算大型工件碳濃度場模擬簡化計算Simplified Calculation of the Large Work piece Carbon Concentration Field平面內(nèi)圓角2.平面圓角滲碳平面內(nèi)圓角不同尖角角度及倒圓半徑對應(yīng)的Lc值/r/mm1234602.323.875.51 7.1880

33、1.81 2.783.874.98901.612.403.294.201001.442.122.793.541201.171.562.052.481400.95 1.131.351.640.02.04.06.08.01234608090100120140Lc/mmr/mm0.02.04.06.08.06080100120140r=1r=2r=3r=4Lc/mm/40大型工件碳濃度場模擬簡化計算大型工件碳濃度場模擬簡化計算Simplified Calculation of the Large Work piece Carbon Concentration Field3.軸對稱尖角滲碳ABFCED

34、軸對稱外凸尖角L定義為EF上某點到F點的距離0.00.20.40.60.81.01.2024對角線一維C/%h/mm軸對稱尖角結(jié)構(gòu)截面圖與平面尖角結(jié)構(gòu)相同，對稱軸位于EF線下方Rcm處41大型工件碳濃度場模擬簡化計算大型工件碳濃度場模擬簡化計算Simplified Calculation of the Large Work piece Carbon Concentration Field3.軸對稱尖角滲碳軸對稱外凸尖角/R/cm12345604.304.284.274.264.26802.942.932.932.932.93902.452.452.452.452.451002.112.112.

35、112.112.111201.661.661.661.661.671401.291.291.291.291.29軸對稱外凸尖角不同尖角角度及工件內(nèi)徑對應(yīng)的Lc值0.01.02.03.04.05.012345608090100120140Lc/mmR/cm0.01.02.03.04.05.06080100120140R=1R=2R=3R=4R=5Lc/mm/42大型工件碳濃度場模擬簡化計算大型工件碳濃度場模擬簡化計算Simplified Calculation of the Large Work piece Carbon Concentration Field3.軸對稱尖角滲碳軸對稱內(nèi)凹尖角L定

36、義為CD上某點到C點的距離0.00.20.40.60.81.01.20123對角線一維C/%h/mm/R/cm12345601.071.071.061.061.06801.041.031.031.031.03901.021.011.011.011.011001.000.990.990.990.991200.950.940.940.940.941400.850.850.850.840.84軸對稱內(nèi)凹尖角不同尖角角度及工件內(nèi)徑對應(yīng)的Lc值0.00.20.40.60.81.01.212345608090100120140Lc/mmR/cm0.00.20.40.60.81.01.26080100120

37、140R=1R=2R=3R=4R=5Lc/mm/43大型工件碳濃度場模擬簡化計算大型工件碳濃度場模擬簡化計算Simplified Calculation of the Large Work piece Carbon Concentration Field3.軸對稱圓角滲碳軸對稱外凸圓角L定義為EF上某點到F點的距離軸對稱外凸圓角結(jié)構(gòu)截面圖與平面外圓角結(jié)構(gòu)相同，對稱軸位于EF線下方Rcm處ABDECF2F1F0.00.20.40.60.81.01.20123對角線一維C/%h/mm44大型工件碳濃度場模擬簡化計算大型工件碳濃度場模擬簡化計算Simplified Calculation of th

38、e Large Work piece Carbon Concentration Field3.軸對稱圓角滲碳軸對稱外凸圓角軸對稱外凸圓角不同角度及工件內(nèi)徑對應(yīng)的Lc值/r/mm1234604.264.405.677.23802.943.144.035.07902.462.753.474.301002.142.422.983.651201.701.872.192.6411401.311.391.541.750123456781234608090100120140Lc/mmr/mm0123456786080100120140r=1r=2r=3r=4Lc/mm/45大型工件碳濃度場模擬簡化計算大型工

39、件碳濃度場模擬簡化計算Simplified Calculation of the Large Work piece Carbon Concentration Field3.軸對稱圓角滲碳軸對稱內(nèi)凹圓角L定義為CD上某點到C點的距離軸對稱內(nèi)凹圓角結(jié)構(gòu)截面圖與平面內(nèi)圓角結(jié)構(gòu)相同，對稱軸位于EF線下方Rcm處BACC1C2DEF0.00.20.40.60.81.01.20123對角線一維C/%h/mm46大型工件碳濃度場模擬簡化計算大型工件碳濃度場模擬簡化計算Simplified Calculation of the Large Work piece Carbon Concentration Fie

40、ld3.軸對稱圓角滲碳軸對稱內(nèi)凹圓角軸對稱內(nèi)凹圓角不同角度及工件內(nèi)徑對應(yīng)的Lc值/r/mm1234602.323.855.447.07801.792.773.824.91901.592.383.254.141001.432.062.763.501201.161.541.982.451400.941.111.351.600123456781234608090100120140Lc/mmr/mm0123456786080100120140r=1r=2r=3r=4Lc/mm/47大型工件碳濃度場模擬簡化計算大型工件碳濃度場模擬簡化計算Simplified Calculation of the Lar

41、ge Work piece Carbon Concentration Field4.三維尖角滲碳在遠(yuǎn)離尖角部位碳濃度分布呈現(xiàn)一維特性三維尖角不同滲碳時間對應(yīng)的Lc值及DH值平面外尖角bacLt/h2345678Lc/mm1.691.992.272.512.722.923.11DH/mm1.371.681.942.172.382.572.740.00.51.01.52.02.53.03.52345678LcDHL/mmt/h0.00.51.01.52.02.53.03.51.4 1.6 1.8 2.0 2.2 2.4 2.6 2.8 3.0 LcDHL/mmt0.5/h0.5L值定義示意圖48大

42、型工件碳濃度場模擬簡化計算大型工件碳濃度場模擬簡化計算Simplified Calculation of the Large Work piece Carbon Concentration Field1.均勻溫度場溫度值對碳濃度場的影響-0.03-0.02-0.010.000.010.020.03-10 -8-6-4-202468101h2h3h4h5h6hC/%T/yx問題對碳擴(kuò)散的影響主要體現(xiàn)在擴(kuò)散系數(shù)上0.00.20.40.60.81.01.20123940930920910900C/%x/mm49溫度空間分布對碳濃度場的影響溫度空間分布對碳濃度場的影響Effects of Temper

43、ature Distributions on Carbon Concentration Field 溫度空間分布對碳濃度場的影響溫度空間分布對碳濃度場的影響Effects of Temperature Distributions on Carbon Concentration Field 2.沿滲碳方向溫度線性分布8208408608809009209409600246810T/x/mm滲碳表面保持940，溫度沿滲碳方向線性降低-0.025-0.020-0.015-0.010-0.0050.0000123456789 106h5h4h3h2h1hC/%-dT/dx(/mm)溫度分布-0.025

44、-0.020-0.015-0.010-0.0050.00001231/mm2/mm3/mm4/mm5/mm6/mm7/mm8/mm9/mm10/mmC/%x/mm最大碳濃度差值分布圖將不同溫度分布條件下的碳濃度分布與940均勻溫度場條件下的碳濃度分布進(jìn)行比較503.沿滲碳方向溫度拋物線形分布6407408409401040114012400246810T/x/mm沿x方向溫度分布呈x的一元二次函數(shù)形式：溫度分布將不同溫度分布條件下的碳濃度分布與940均勻溫度場條件下的碳濃度分布進(jìn)行比較-0.010-0.008-0.006-0.004-0.0020.0000.0020.0040.0060.008

45、0.0100123a=-2.5a=-2.0a=-1.5a=-1.0a=-0.5a=0a=0.5a=1.0a=1.5a=2.0a=2.5C/%x/mm-0.010-0.008-0.006-0.004-0.0020.0000.0020.0040.0060.0080.010-2.5-1.5-0.50.51.52.51h2h3h4h5h6hC/%a51溫度空間分布對碳濃度場的影響溫度空間分布對碳濃度場的影響Effects of Temperature Distributions on Carbon Concentration Field 4.垂直于滲碳方向溫度線性分布88090092094096098010000246810T/y/mm0.0000 0.0001 0.0002 0.0003 0.0004 0.0005 0.0006 2468101h2h3h4h5h6hC/%dT/dy(/mm)5.垂直于滲碳方向溫度拋物線形分布860880900920940960980100010