【正文】 他們是我學(xué)習(xí)和生活上的良師益友，與他們的討論和交流都使我取得進(jìn)步 和提高。馬老師豐富的教學(xué)和實(shí)踐經(jīng)驗(yàn)，嚴(yán)謹(jǐn)?shù)闹螌W(xué)態(tài)度、為人處事的方法、忘我的工作熱情、一絲不茍的工作作風(fēng)和對(duì)學(xué)生多方面的關(guān)心深深地教育和影響我，讓我終身受益。重軌的淬火過(guò)程包括對(duì)重軌的加熱、保溫、噴風(fēng)冷卻和自然冷卻等方面，在未來(lái)若是對(duì)重軌的淬火過(guò)程進(jìn)行進(jìn)一步研究，則會(huì)對(duì)整個(gè)淬火過(guò)程的各個(gè)方面進(jìn)行更加全面的仿真，以更深的了解淬火過(guò)程的特點(diǎn)，并獲得更多的研究數(shù)據(jù)和結(jié)論，為實(shí)際生產(chǎn) 的正常運(yùn)行提供更有力的參數(shù)參考。本研究可 為實(shí) 25 際生產(chǎn)的正常運(yùn)行提供理論依據(jù)和數(shù)據(jù)參考。 (2) 在熱分析過(guò)程中充分考慮非線性的熱物性參數(shù)：比熱、導(dǎo)熱系數(shù)、對(duì)流化熱系數(shù)隨溫度的變化。主要研究結(jié)論如下： (1) 根據(jù)傳熱學(xué)及有限元理論，建立了重軌淬火過(guò)程中的三維傳熱模型，模擬真實(shí)工況，施加特定的邊界條件，并對(duì)多種噴風(fēng)壓力情況進(jìn)行有限元分析。 由 ANSYS仿真結(jié)果可知，在相同的加熱時(shí)間和保溫時(shí)間的條件下， ，具有冷熱快、殘余應(yīng)力小的優(yōu)點(diǎn)，所以在實(shí)際生產(chǎn)的淬火冷卻過(guò)程中更宜用 風(fēng)冷卻，這樣會(huì)提高生產(chǎn)中的淬火效率。 不同壓強(qiáng)下噴風(fēng)應(yīng)力場(chǎng)對(duì)比分析 （ 1） 圖 49 t=30s的 應(yīng) 力場(chǎng) 圖 410 t=45s的 應(yīng)力場(chǎng) 23 圖 411 t=60s的 應(yīng)力場(chǎng) （ 2） 圖 412 t=30s的 應(yīng)力場(chǎng) 圖 413 t=45s的 應(yīng)力場(chǎng) 2 24 圖 414 t=60s的 應(yīng)力場(chǎng) 由圖 49到 414可以得到 t=30s、 t=45s和 t =60s的最大應(yīng)力分別是 4555Mpa、 8547Mpa和 6345Mpa, ，分 別是 4686Mpa、 8791Mpa和 6529 Mpa。 圖 41 圖 42 不同壓強(qiáng)下噴風(fēng)溫度場(chǎng)對(duì)比分析 （ 1） 2 20 圖 43 t=30s的溫度場(chǎng) 圖 44 t=45s的溫度場(chǎng) 圖 45 t=60s的溫度場(chǎng) 21 （ 2） 圖 46 t=30s的溫度場(chǎng) 圖 47 t=45s的溫度場(chǎng) 圖 48 t=60s的 溫度場(chǎng) 2 22 從上面的這些 ANSYS通用后處理器得到的結(jié)果圖（ 43到 48），可以得到 t=30s、 t=45s和 t =60s的最小溫度分別是 874Mpa、 862Mpa和851Mpa, ，分別是 887Mpa、 881Mpa和 875 Mpa。 圖 318 重軌應(yīng)力場(chǎng)分布圖 19 第四章 重軌淬火過(guò)程的溫度場(chǎng)和應(yīng)力場(chǎng)分析 研究不同壓強(qiáng)下溫度場(chǎng)和應(yīng)力場(chǎng)的前提條件 重軌的表面換熱系數(shù)決 定重軌在冷淬過(guò)程中溫度的變化，圖 41 和 42 分別是重軌在 M pa 、 M pa 冷氣的條件下表面對(duì)流換熱系數(shù)隨重軌表面溫度的變化情況。 12. 查看應(yīng)力場(chǎng)分布 在主菜單中選擇【 General Postproc】 【 Plot Results】 【 Contour Plot】 【 Nodal Solution】命令，在彈出的 Contour Nodal Solution Data 對(duì)話框中，選擇【 Nodal Solution】 【 Stress】 【 von Mises stress】，如圖 318 所示。求解完畢后會(huì)彈出提示對(duì)話框。完成選擇后單擊【 OK】按鈕。 （ 2） 施加溫度載荷。再選擇【 Thermal Expansion】 【 Secant coefficient】 【 Isotropic】 ,如圖 317 所示輸入線膨脹系數(shù)。在主菜單中選擇【 Preprocessor】 【 Material Props】 【 Material Model】命令，選擇【 Structural】選項(xiàng)。 圖 315 單元轉(zhuǎn)換 （ 2） 為結(jié)構(gòu)分析材料屬性。在主菜單中選擇【 Preferences】命令，彈出設(shè)定分析模塊對(duì)話框，勾選【 Structural】復(fù)選框 .單擊【 OK】按鈕。單擊 OK，則得到第 30s 時(shí)重軌溫度場(chǎng)分布圖 (如圖 314)。 在主菜單中選擇【 General Postproc】 【 Read Results】 【 By Pick】 ,在彈出的 Results File 對(duì)話框，選擇時(shí)間為 30 的一項(xiàng)，再單擊 Read 按鈕。求解完畢后會(huì)彈出提示對(duì)話框。單擊 OK 按鈕，完成結(jié)果輸出控制的設(shè)置。在 Item 下拉列表中選擇 All Item 表明輸出所有計(jì)算結(jié)果。 (2) 設(shè)置結(jié)果輸出項(xiàng)。在主菜單中選擇【 Solution】 【 Load Step Opts】2 16 【 Time/Frequenc】 【 TimeTime Setup】命令，彈出時(shí)間設(shè)置對(duì)話框，在 Time一欄中輸入本分析的最終時(shí)間 60，在時(shí)間步長(zhǎng)一欄中輸入步長(zhǎng)時(shí)間 5；在加載方式單選按鈕中選中 Stepped 單選按鈕，如圖 313 所示。選中重軌表面節(jié)點(diǎn) ,單擊拾取對(duì)話框的 OK 按鈕，此時(shí)彈出對(duì)流邊界條件設(shè)置對(duì)話框，在其中 Film coefficient 一欄輸入 100，表明對(duì)流傳熱系數(shù)為 100；在 Bulk temperature 一欄輸入 5，表 明邊界溫度為 5 攝氏度，如圖 312。單擊【 OK】按鈕，完成初始溫度的設(shè)置 . （ 3）為重軌施加外界條件。在主菜單中選擇【 Solution】 【 Define Loads】【 Apply】 【 Initial Condit ’n】 【 Define】命令，在彈出的節(jié)點(diǎn)拾取對(duì)話框中，單擊 Pick All 按鈕后會(huì)立即彈出初始條件定義對(duì)話框。單擊【 OK】按鈕，完成分析類型的選擇。 圖 39 劃分網(wǎng)格 圖 310 網(wǎng)格圖 6. 施加邊界條件 15 （ 1） 設(shè)置分析類型。 圖 38 導(dǎo)入的模型 5. 劃分網(wǎng)格 在主菜單中選擇【 Preprocessor】 【 Meshing】 【 Mesh Tool】命令，彈出網(wǎng)格劃分選擇的對(duì)話框，如圖 39 勾選各項(xiàng)。選擇【 Thermal】選項(xiàng)，分別雙擊【 Conductivity】、【 Specific Heat】和【 Density】，輸入比熱容、導(dǎo)熱系數(shù)和密度值，最后結(jié)果如圖37?！?,密度為 7800kg/m3。 13 圖 36 選擇單元類型 3. 定義材料屬 性 定義材料的比熱容為 623J/kg在主菜單中選擇【 Preferences】命令，彈出設(shè)定分析模塊對(duì)話框，勾選【 Structural】和【 Thermal】復(fù)選框（如圖 35），然后單擊【 OK】按鈕完成分析模塊的選擇。 （ 4） 設(shè)定分析模塊。 2 12 （ 3） 設(shè)置工作標(biāo)題。 （ 2） 設(shè)置文件名。 ANSYS 仿真重軌溫度場(chǎng)和應(yīng)力場(chǎng)的基本步驟 1. 設(shè)置環(huán)境（設(shè)置工作目錄，設(shè)置文件名，設(shè)置工作標(biāo)題，設(shè)定分析模塊） （ 1） 設(shè)置工作目錄?！?1) 26 ：重軌冷噴前，表面溫度為 900℃ ，淬火冷空氣的溫度 5～ 10℃ ，冷氣的對(duì)流系數(shù) 100W/m2℃ 1) 478 487 505 533 579 650 772 760 623 K/(W 表 1 鋼軌的物性參數(shù) T/℃ 100 200 300 400 500 600 700 800 900 c/(J 圖 32 重軌的三維模型圖 確定重軌的各項(xiàng)材料參數(shù)及初始條件 50kg/m重軌密度為 7800kg /m3, 線膨脹系數(shù)。再運(yùn)用熱 結(jié)構(gòu)耦合分析的方法，以溫度作為求解應(yīng)力的初始條件對(duì)有限元模型進(jìn)行加載求解，可得到重軌的應(yīng)力場(chǎng)仿真結(jié)果。再將重軌的三維模型導(dǎo)入 ANSYS 軟件中。珠光體和貝氏體的等溫相變體積百分比按 MehlJohnson? 所提出的公式： )exp(1 knkk tbV ??? （ ） V為 k相的體積百分含量； kb 、 kn 仇為兩個(gè)與溫度有關(guān)的參數(shù)，它們可以從材料的等溫相變圖求得； k=1對(duì)應(yīng)珠光體， k=2對(duì)應(yīng)貝氏體。淬火的過(guò)程就是奧氏體向珠光體轉(zhuǎn)變的過(guò)程，使珠光體在較低的溫度下轉(zhuǎn)變，而得2 10 到強(qiáng)韌性和細(xì)片狀珠光體 (即索氏體 )組織，不得出現(xiàn)馬氏體、貝氏體有害組織，因此在淬火應(yīng)力模擬計(jì)算中僅考慮索氏體 (細(xì)珠光體 )的相變。對(duì)于不同的單元，要作不同處理。這種解法，每加載一次或迭代一次，剛度矩陣都要根據(jù)應(yīng)力水平重新計(jì)算一次，故稱這種方法為增量變剛度法。 增量變剛度法的特點(diǎn)是將載荷分段逐步增加，增加一次載荷，就會(huì)產(chǎn)生應(yīng)力和應(yīng)變?cè)隽????? 和 ???? ，只要增加載荷適當(dāng)?shù)男?，則可近似認(rèn)為在此計(jì)算步內(nèi)??? 保持不變，則 ? ? ? ??? ??? epD （ ） 成線性關(guān)系。具體表達(dá)式為： 變分方程： ? ?? ? ? ?hRK ??