【正文】 =60s的最大應(yīng)力分別是 4555Mpa、 8547Mpa和 6345Mpa, ，分 別是 4686Mpa、 8791Mpa和 6529 Mpa。求解完畢后會彈出提示對話框。在主菜單中選擇【 Preprocessor】 【 Material Props】 【 Material Model】命令，選擇【 Structural】選項。 在主菜單中選擇【 General Postproc】 【 Read Results】 【 By Pick】 ,在彈出的 Results File 對話框，選擇時間為 30 的一項，再單擊 Read 按鈕。 (2) 設(shè)置結(jié)果輸出項。在主菜單中選擇【 Solution】 【 Define Loads】【 Apply】 【 Initial Condit ’n】 【 Define】命令，在彈出的節(jié)點拾取對話框中，單擊 Pick All 按鈕后會立即彈出初始條件定義對話框。選擇【 Thermal】選項，分別雙擊【 Conductivity】、【 Specific Heat】和【 Density】，輸入比熱容、導(dǎo)熱系數(shù)和密度值，最后結(jié)果如圖37。 （ 4） 設(shè)定分析模塊。℃ 1) 26 ：重軌冷噴前，表面溫度為 900℃ ，淬火冷空氣的溫度 5～ 10℃ ，冷氣的對流系數(shù) 100W/m2再運用熱 結(jié)構(gòu)耦合分析的方法，以溫度作為求解應(yīng)力的初始條件對有限元模型進行加載求解，可得到重軌的應(yīng)力場仿真結(jié)果。對于不同的單元，要作不同處理。 0)( ?? ? pdH ?? （ ） 4)塑 性區(qū)的行為服從流動法則。 2)重軌材料的屈服服從 Mises 屈服準(zhǔn)則。以往的研究結(jié)果表明，重軌及淬火時的換熱邊界條件是非線性的，主要是由淬火時相變的產(chǎn)生和淬火材料的物性參數(shù)隨溫度而變化等原因引起的。 （ 3）第三類邊界條件，又稱牛頓邊界條件，是指重軌與其相接觸的流體介質(zhì)間的熱對流系數(shù)以和介質(zhì)溫度乃為已知。 重軌淬火的邊界條件 重軌淬火時的邊界條件 【 13】 是指淬火工件外表面與周圍環(huán)境的熱交換情況。物理條件是指導(dǎo)熱重 軌的主要物理參數(shù)和物理特征。熱對流用牛頓冷卻方程來描述： )( BSn TThq ?? （ ） 式（ ）中， h 為對流換熱系數(shù)； ST 為固體表面溫度； BT 為周圍流體的溫度。 研究內(nèi)容 1． 依據(jù)傳熱原理，建立重軌的三維模型，用有限元方法研究重軌軌頭在淬火過程中的溫度場應(yīng)力場分布； 2. 研究不同的壓強以及淬火噴風(fēng)時間對重軌溫度場、應(yīng)力場的影響，找出最佳的淬火噴風(fēng)壓強； 3. 結(jié)合熱模擬實驗，為實際生產(chǎn)的正常運行提供理論依據(jù)。盡管從數(shù)學(xué)模型、物性參數(shù)、換熱系數(shù)等都進行了大量的簡化，但是在計算結(jié)果的總體趨勢還是非常有意義的。 國外對淬火試件淬火數(shù)值模擬計算比較早，數(shù)值模擬研究始于 70年代， 。這是因為，一方面是對冷卻介質(zhì)的應(yīng)用還存在著一定的片面性，另一方面是對冷卻介質(zhì)的應(yīng)用在水平上還沒有達到很深的層次上，甚至還停留在表面的層次上。淬火軌硬度可至合適并均勻一致。因此，淬火過程的深入研究對工程實際大有重要的指導(dǎo)意義，利 用計算機進行數(shù)值模擬 【 4】 有助于淬火工藝設(shè)計，便于選擇合適的淬火工藝調(diào)整方案，可以大大減少試驗量，具有一定的實用價值，已成為當(dāng)今熱處理領(lǐng)域的研究熱點之一。 重軌淬火是提高其韌性和耐磨性的主要途徑之一。計算了強制冷卻、空氣自然對流等淬火過程的溫度場和應(yīng)力場分布情況，分析淬火時間對溫度場和應(yīng)力場的影響。 在 數(shù)值模擬計算的過程中，輸入在不同的噴風(fēng)壓力下的對流換熱系數(shù)，得到相應(yīng)的溫度場和應(yīng)力場結(jié)果，并對結(jié)果進行了分析。 這 些淬火 過程出 現(xiàn) 的 問題 ，可能 會導(dǎo) 致重 軌軌頭 出 現(xiàn) 掉 塊 、裂 紋 等 現(xiàn) 象， 從 而影 響 列 車 的正常 運 行、降低 鐵 路的使用效率，甚至 對 列 車 的安全 運 行造成 隱 患。淬火過程是一個各種場相互耦合的復(fù)雜過程，要在理論上求解各場量的解析解是非常困難的，甚至是不可能的。當(dāng)工藝參數(shù)優(yōu)選后，只要重軌含碳量大于 0． 7％，使用壓縮空氣可以使奧氏體實現(xiàn)向索氏體的轉(zhuǎn)變，并且能夠保證熱處理后的材質(zhì)內(nèi)部結(jié)構(gòu)均勻。主要的問題仍然是淬火重軌的硬度偏低。它能對試件的溫度場、顯微組織場和內(nèi)應(yīng)力場進行耦合計算，給出每一瞬間的溫度場面、顯微組織場和內(nèi)應(yīng)力場，并能直接地觀察到各場量在淬火過程中的變化情況，這樣就可以在節(jié)省大量的人力、物力、財務(wù)和時問的情況下對試件進行全面的分析，預(yù)測試件淬火后的組織性能，從而可對淬火工藝方案進行優(yōu)化，使工藝更加高效合理。 20世紀 80年代中期，上海重機廠等單位開始采用有限差分法，對大型軋輥淬火過程溫度場進行了計算。 目前，對淬火試件內(nèi)應(yīng)力場的研究主要集中在淬火過程中試件的瞬態(tài)溫度、相應(yīng)、應(yīng)力之間的耦合作用，其目的是能準(zhǔn)確地反映出實際過程各場量的變化，較為準(zhǔn)確地預(yù)測出淬火后試件內(nèi)部組織、 性能、應(yīng)力和畸變。熱對流可以分為兩類：自然對流和強制對流。幾何條件是指參與導(dǎo)熱過程的重軌的幾何尺度、形狀。此時 )(00 zyxTT t 、?? （ ） 式（ ）中， )(0 zyxT 、 為已知溫度函數(shù)。q 7 的方向等同于邊界法線 n的方向，其表達式為 ）、（或者 zyxtqnTqnT ss ?? ?? ???????? （ 2. 7） 式（ ）中， ?q 為已知重軌表面的熱流密度，為定值 ( )/( 2mW )；）、（ zyxtq? 為己知重軌的熱流密度函數(shù)，隨位置、時間而變化。即使在線性換熱邊界條件中，當(dāng)已知的溫度、換熱系數(shù)、熱流密度等參數(shù)隨著材料的物性參數(shù)、導(dǎo)熱時間及材料 中溫度場的變化而發(fā)生變化時，線性換熱邊界條件就變?yōu)榉蔷€性換熱邊界條件。 根據(jù)體積不變定律，可導(dǎo)出塑性變形時的泊松系數(shù)為 ? =。在復(fù)雜應(yīng)力狀態(tài)下，設(shè)新的屈服應(yīng)力只與卸載前的等效塑性應(yīng)變總量有關(guān)，即只有當(dāng)應(yīng)力適合式時才會發(fā)生塑性變形。 用增量變剛度法求解彈塑性問題時，不是在每一次加載后所有單元都處于彈性或塑性狀態(tài)，實際上可能是下列三種狀態(tài)并存： (1)單元處于彈性狀態(tài)，稱為彈性單元； (2)單元處于塑性狀態(tài)，稱為塑性單元； (3)單元 在加載前處于彈性狀態(tài)，加載后處于塑性狀態(tài)，稱為過渡單元。利用 ANSYS 軟件，以已知的初始條件和邊界條件，得到重軌的溫度場仿真圖。 m1在通用菜單中選擇【 File】 【 Change Title】命令，彈出【 Change Title】對話框，輸入文件名“ zhonggui”（如圖 34），單擊【 OK】按鈕。在主菜單中選擇【 Preprocessor】 【 Material Props】 【 Material Model】命令，彈出材料屬性定義對話框。 （ 2） 設(shè)置重軌的初始溫度。單擊 OK 按鈕，完成時間和載 荷步的設(shè)置。單擊提示對話框中的 OK 按鈕，則求解結(jié)束。 線膨脹系數(shù) ，泊松比 ，彈性模量 220Gpa。 11. 求解，結(jié)構(gòu)分析 在主菜單中選擇【 Solution】 【 Solve】 【 Current LS】命令，彈出求解確認對話框，單擊 OK 按鈕，求解開始。 可以看出 冷卻能力比 ，在同樣加熱、保溫、冷卻時間相同條件下，前者溫度下降的快。計算模擬淬火過程中的溫度場和應(yīng)力場，并對仿真結(jié)果進行比較分析，了解了不同的噴風(fēng)壓力對溫度場 和應(yīng)力場的影響。 27 參考文獻 【 1】 王勖成．有限單元法【 M】．清華大學(xué)出版社． 20xx 【 2】 劉莊，吳肇基，吳景之等．熱處理過程的數(shù)值模擬【 M】．第一版．北京：科學(xué)出版社， 1996 【 3】 黃日升．鋼軌在整體加熱全長淬火熱處理過程中的應(yīng)力分布和變形的粗淺探討【 J】．包鋼科技， 1997 【 4】 謝建斌，程赫明，何天淳等． 1045 銅 ]淬火時溫度場的數(shù)值模擬【 J】．甘肅工業(yè) 大學(xué)學(xué)報， 20xx 【 5】 侯立軍． TIO 鋼高壓氣體淬火過程中的數(shù)值模擬【 D】．昆明理工大學(xué)，20xx 【 6】 Michael J． CIALKOWSKI, ndrzej FRACKOWIAK． Solution of the stationary 2D Inverse Heat Conduction Problem by Treffetz Method【 J】 ． J． of Thermal Science， 20xx 【 7】 薛守義．有限單元法【 M】．中國建材出版社． 20xx 【 8】 馮曉丹，張玉．基于 ANSYS軟件的 T8鋼淬火過程三維溫度場的模擬【 J】．模具制造， 20xx 【 9】 秦宇 . 基礎(chǔ)與實例教程 . 北京：化學(xué)工業(yè)出版社， 20xx 【 10】 楊世銘，陶文銓．傳熱學(xué)【 M】．第四版．北京高等教育出版社， 20xx 【 11】 梁木子，文杰，程先華等．蒸汽管道外接座溫度場和熱應(yīng)力場的有限元計算【 J】．上海交通大學(xué)學(xué)報， 20xx 【 12】 付斌，安濤，沈曉輝等．火車車輪淬火過程中的溫度場【 J】．鋼鐵研究學(xué)報， 20xx 【 13】 田東，胡明娟，潘劍生． T8鋼淬火過程三維溫度場計算及實驗【 J】 ．上海交通大學(xué)學(xué)報， 1998 【 14】 武斌斌，溫殿英，李連進．用有限元方法分析無縫鋼管的淬火殘余應(yīng)力【 J】．重型機械， 20xx