國民經(jīng)濟(jì)和現(xiàn)代科學(xué)技術(shù)的快速發(fā)展,對我國機(jī)械制造業(yè)提出了精密、高效、經(jīng)濟(jì)、清潔、產(chǎn)業(yè)化等方面的要求,對熱加工零件的形狀、尺寸精度和力學(xué)性能等方面的要求也越來越高。熱處理工藝作為主要的熱加工工藝,在控制熱加工零件的形狀、尺寸精度和力學(xué)性能等方面發(fā)揮著重要作用。對于大型、精密、復(fù)雜的熱加工零件,僅依靠經(jīng)驗(yàn)和傳統(tǒng)的理論知識完成其熱處理工藝設(shè)計(jì)有較大難度。熱處理工藝是改善材料性能的重要手段,但熱處理變形和殘余應(yīng)力一直困擾著工程技術(shù)人員,是制約機(jī)械零件、工裝模具制造精度及使用壽命的難題之一。
熱處理技術(shù)的發(fā)展更傾向于利用虛擬的設(shè)計(jì)-制造-驗(yàn)證一體化環(huán)境, 將真實(shí)的設(shè)計(jì)、制造、材料、驗(yàn)證、應(yīng)用乃至維修和全生命周期管理等諸多環(huán)節(jié)統(tǒng)一起來,從而最大限度地縮短新產(chǎn)品研發(fā)周期,降低研發(fā)成本, 提高產(chǎn)品的市場競爭力。在這個過程中,計(jì)算機(jī)輔助工程(CAE)技術(shù)已成為創(chuàng)新設(shè)計(jì)、數(shù)字化設(shè)計(jì)和材料制造技術(shù)的核心之一。CAE已被廣泛應(yīng)用至鍛造、擠壓、熱沖壓、軋制、熱處理等熱加工工藝的設(shè)計(jì),并取得了較好的效果。美國在2010年發(fā)布的新版熱處理技術(shù)路線圖中將虛擬熱處理作為重點(diǎn)發(fā)展方向。中國工程院在2013年年底制定的中國熱處理技術(shù)與表層改性路線圖,也將虛擬熱處理作為我國熱處理領(lǐng)域的十二個重點(diǎn)研究內(nèi)容之一。目前,高等學(xué)校的師生和企業(yè)的技術(shù)人員越來越重視材料熱加工過程中的溫度、組織、相變、應(yīng)力等物理量的數(shù)值模擬。
編寫一本包含一定的理論知識和工程背景的熱處理數(shù)值模擬方面的專著, 是作者多年的夢想。本書的內(nèi)容多為作者近年來發(fā)表的一些研究結(jié)果、學(xué)習(xí)心得以及指導(dǎo)研究生的成果,并吸收了國內(nèi)外同行的研究成果。本書包括溫度、組織轉(zhuǎn)變、應(yīng)力/應(yīng)變、濃度場、晶粒長大、多物理場耦合分析等方面的數(shù)值模擬技術(shù),可作為高等學(xué)校的教材,也可以作為企業(yè)技術(shù)人員的指導(dǎo)讀物。在本書的撰寫過程中,山東大學(xué)博士生導(dǎo)師趙國群教授給了我們悉心的指導(dǎo), 山東大學(xué)和山東科技大學(xué)的其他老師也給了我們不少幫助,在此向他們表示深深的感謝。
由于計(jì)算機(jī)在熱處理領(lǐng)域中的應(yīng)用非常廣泛,且計(jì)算機(jī)軟硬件技術(shù)、數(shù)值模擬技術(shù)的發(fā)展日新月異,熱處理數(shù)值模擬的新技術(shù)和新方法不斷出現(xiàn), 加之編者學(xué)識有限,書中難免有不當(dāng)之處,敬請讀者批評指正。
著者
2017年5月
第1章 概述 1
1.1 引言 1
1.2 虛擬熱處理的基本概念 2
1.3 淬火工藝模擬技術(shù)的國內(nèi)外研究現(xiàn)狀 3
1.3.1 國外淬火工藝模擬研究現(xiàn)狀 3
1.3.2 國內(nèi)淬火工藝模擬研究現(xiàn)狀 3
1.3.3 國內(nèi)外的熱處理軟件包 4
1.3.4 熱處理過程的優(yōu)化 6
1.4 淬火過程數(shù)值模擬的難點(diǎn)及存在的問題 7
第2章 淬火工藝溫度場模擬技術(shù) 9
2.1 引言 9
2.2 淬火過程導(dǎo)熱偏微分方程 10
2.2.1 溫度場控制方程 10
2.2.2 初始條件 10
2.2.3 邊界條件 11
2.3 瞬態(tài)溫度場的變分 12
2.3.1 平面瞬態(tài)溫度場的變分 12
2.3.2 軸對稱瞬態(tài)溫度場的變分 18
2.4 瞬態(tài)溫度場的求解 19
2.4.1 差分方法 19
2.4.2 系數(shù)矩陣的存儲方法 20
2.4.3 溫度場數(shù)值振蕩問題 21
2.5 熱物性參數(shù)的選擇 28
2.6 溫度場計(jì)算流程框圖 28
2.7 溫度場有限元模擬程序驗(yàn)證 30
2.7.1 變熱導(dǎo)率定常內(nèi)熱的一維穩(wěn)態(tài)熱傳導(dǎo)問題 30
2.7.2 內(nèi)熱二維瞬態(tài)熱傳導(dǎo)問題 31
本章小結(jié) 33
第3章 淬火工藝相變過程模擬技術(shù) 34
3.1 引言 34
3.2 TTT 曲線 35
3.3 相變過程的數(shù)學(xué)模型 35
3.3.1 擴(kuò)散型轉(zhuǎn)變 35
3.3.2 非擴(kuò)散型轉(zhuǎn)變 36
3.3.3 馬氏體相變溫度的計(jì)算 36
3.3.4 貝氏體相變溫度的計(jì)算 37
3.3.5 相變潛熱的計(jì)算與處理 37
3.4 Scheil疊加法則 37
3.5 杠桿定律 39
3.6 淬火過程的相變塑性 40
3.7 淬火力學(xué)性能計(jì)算 41
3.8 組織場模擬流程框圖 41
3.9 P20端淬工藝模擬與實(shí)驗(yàn)研究 43
3.9.1 端淬工藝模擬 43
3.9.2 端淬實(shí)驗(yàn)研究 44
3.9.3 相變潛熱對溫度場和組織場的影響 52
本章小結(jié) 53
第4章 淬火過程冷卻曲線的采集及介質(zhì)傳熱系數(shù)的計(jì)算 55
4.1 引言 55
4.2 計(jì)算模型及計(jì)算方法 56
4.2.1 計(jì)算模型的建立 56
4.2.2 傳熱系數(shù)優(yōu)化區(qū)間的確定 56
4.2.3 傳熱系數(shù)最佳值的確定 58
4.2.4 黃金分割法迭代次數(shù)的分析 60
4.3 傳熱系數(shù)的求解 61
4.4 實(shí)驗(yàn)裝置 65
4.4.1 實(shí)驗(yàn)工裝 65
4.4.2 熱電偶 66
4.4.3 熱電偶調(diào)理板 68
4.4.4 數(shù)據(jù)采集卡 68
4.5 冷卻曲線的采集及傳熱系數(shù)計(jì)算 71
本章小結(jié) 76
第5章 淬火過程應(yīng)力/應(yīng)變場的模擬技術(shù) 78
5.1 引言 78
5.2 淬火過程力學(xué)基本方程 79
5.3 熱彈塑性本構(gòu)關(guān)系 80
5.3.1 彈性區(qū)的應(yīng)力應(yīng)變關(guān)系 80
5.3.2 塑性區(qū)的應(yīng)力應(yīng)變關(guān)系 82
5.3.3 過渡區(qū)的彈塑性比例系數(shù)的計(jì)算 84
5.4 應(yīng)力/應(yīng)變場有限元基本理論與技術(shù) 86
5.4.1 單元和形函數(shù) 86
5.4.2 單元應(yīng)變速率矩陣 87
5.4.3 等效應(yīng)變速率矩陣 89
5.4.4 邊界條件 89
5.5 熱彈塑性問題求解 89
5.5.1 變分方程及剛度矩陣 89
5.5.2 增量變剛陣方法 90
5.5.3 迭代收斂準(zhǔn)則 92
5.6 預(yù)應(yīng)力淬火過程的應(yīng)力、應(yīng)變計(jì)算 92
5.7 應(yīng)力、應(yīng)變計(jì)算流程圖 94
5.8 應(yīng)力/應(yīng)變計(jì)算程序檢驗(yàn) 95
本章小結(jié) 98
第6章 淬火過程溫度、相變和應(yīng)力的耦合分析 100
6.1 引言 100
6.2 耦合分析程序流程框圖 102
6.3 耦合分析有限元模型 103
6.4 溫度、相變及應(yīng)力應(yīng)變耦合分析 104
6.4.1 溫度場的模擬 104
6.4.2 組織場的模擬 106
6.4.3 應(yīng)力/應(yīng)變場模擬 108
6.5 彈塑性區(qū)域的演變 115
6.6 淬火零件的變形 117
本章小結(jié) 119
第7章 滲碳工藝有限元模擬關(guān)鍵技術(shù)研究 121
7.1 引言 121
7.2 滲碳工藝有限元模擬 123
7.2.1 基本條件 123
7.2.2 瞬態(tài)濃度場的變分 123
7.2.3 有限差分法 130
7.2.4 濃度場的數(shù)值振蕩問題 131
7.3 有限元模擬程序的實(shí)驗(yàn)驗(yàn)證 132
7.3.1 圓柱體的實(shí)驗(yàn)與模擬 132
7.3.2 齒輪的實(shí)驗(yàn)與模擬 133
第8章 基于MC方法的組織模擬關(guān)鍵技術(shù)研究 139
8.1 晶粒長大MC Potts模型模擬關(guān)鍵技術(shù) 139
8.1.1 傳統(tǒng)晶粒長大模型關(guān)鍵技術(shù) 139
8.1.2 對傳統(tǒng)晶粒長大Exxon MC Potts模型的改進(jìn) 141
8.1.3 新模型模擬計(jì)算機(jī)算法流程 149
8.2 再結(jié)晶MC Potts新模型模擬關(guān)鍵技術(shù) 150
8.2.1 傳統(tǒng)再結(jié)晶模型關(guān)鍵技術(shù) 150
8.2.2 再結(jié)晶新模型模擬流程 151
本章小結(jié) 152
第9章 氣體淬火過程工藝參數(shù)的優(yōu)化 153
9.1 引言 153
9.2 曲面響應(yīng)模型 153
9.3 回歸模型的方差分析 154
9.4 逐步回歸分析 156
9.5 氣體淬火工藝及工藝參數(shù)評估 157
9.5.1 氣體淬火技術(shù) 157
9.5.2 有限元模型 158
9.5.3 目標(biāo)函數(shù)的建立 158
9.5.4 工藝參數(shù)評估 159
9.6 階段性傳熱系數(shù)模型 161
9.6.1 設(shè)計(jì)變量的確定 161
9.6.2 Box-Behnken實(shí)驗(yàn)設(shè)計(jì) 162
9.6.3 響應(yīng)曲面的擬合 164
9.6.4 優(yōu)化目標(biāo)函數(shù)的建立 166
9.6.5 工藝參數(shù)的優(yōu)化結(jié)果 166
9.7 區(qū)域性傳熱系數(shù)模型 170
9.7.1 設(shè)計(jì)變量的確定 170
9.7.2 中心復(fù)合實(shí)驗(yàn)設(shè)計(jì) 171
9.7.3 響應(yīng)曲面的擬合 173
9.7.4 優(yōu)化目標(biāo)函數(shù)的建立 175
9.7.5 工藝參數(shù)的優(yōu)化結(jié)果 176
本章小結(jié) 179
參考文獻(xiàn) 182