ANSYS的絲錐淬火冷卻過程計(jì)算機(jī)模擬
2013-06-10 by:廣州有限元分析、培訓(xùn)中心-1CAE.COM 來源:仿真在線
作者: 韓斌慧*閻獻(xiàn)國(guó) 來源: 萬方數(shù)據(jù)
關(guān)鍵字: 絲錐 淬火 溫度場(chǎng) 應(yīng)力 計(jì)算機(jī)模擬
以現(xiàn)有的絲錐淬火工藝規(guī)程為依據(jù),使用ANSYS軟件對(duì)絲錐淬火冷卻過程進(jìn)行了模擬,得出了溫度場(chǎng)以及應(yīng)力場(chǎng)的等值線分布圖,為正確預(yù)測(cè)刀具溫度場(chǎng)、應(yīng)力場(chǎng),解決高硬度刀具溫度場(chǎng)難于在線測(cè)量問題提供了比較可行、可靠而且簡(jiǎn)便的辦法。該方法對(duì)于同類型刀具具有普遍意義。
熱處理是保證制造業(yè)產(chǎn)品的內(nèi)在質(zhì)量,提高其使用性能和可靠性的重要環(huán)節(jié)。傳統(tǒng)的熱處理是一種依靠經(jīng)驗(yàn)型技術(shù)實(shí)現(xiàn)的熱處理,經(jīng)驗(yàn)型的技術(shù)很難溶人以信息化為主導(dǎo)的先進(jìn)制造技術(shù)之中。因此迫切需要將熱處理改造為基于知識(shí)的熱處理技術(shù),從而減弱熱處理信息化在制造業(yè)中的瓶頸問題。
對(duì)熱處理過程的計(jì)算機(jī)模擬是保證熱處理信息化的關(guān)鍵。金屬切削刀具是制造業(yè)中最基本的生產(chǎn)要素,絲錐則屬于金切刀具中復(fù)雜刀具的一種,對(duì)其熱處理過程的計(jì)算機(jī)模擬一旦有所突破,將會(huì)最終帶動(dòng)整個(gè)刀具制造行業(yè),從而影響到整個(gè)制造業(yè)的熱處理信息化進(jìn)程。由于淬火過程直接影響到刀具的最終使用性能,而且屬于熱處理過程中溫度和應(yīng)力變化最大的一個(gè)環(huán)節(jié),組織變化復(fù)雜,所以以此作為研究的依據(jù)。
熱處理過程還具有“瞬息萬變”的特性,欲實(shí)時(shí)獲取某瞬時(shí)、某特定冷卻條件下的溫度場(chǎng)、應(yīng)力場(chǎng)分布情況,在現(xiàn)有的測(cè)試手段、試驗(yàn)方法以及測(cè)試儀器下還顯得難度較大。
隨著計(jì)算機(jī)軟、硬件水平的提高,伴隨著數(shù)值計(jì)算方法、計(jì)算傳熱學(xué)、熱彈塑性理論、相變動(dòng)力學(xué)、計(jì)算流體力學(xué)等學(xué)科的同步跟進(jìn),對(duì)于絲錐的淬火過程計(jì)算機(jī)模擬已成為可行、可靠的技術(shù)手段。
1 淬火過程計(jì)算機(jī)模擬的客觀依據(jù)
計(jì)算機(jī)模擬只是一種手段,其模擬結(jié)果的可靠程度要通過試驗(yàn)或現(xiàn)行熱處理工藝進(jìn)行檢驗(yàn)。此處嚴(yán)格按照某工具廠M24絲錐淬火工藝規(guī)程進(jìn)行:對(duì)于W9Mo3Cr4V材料制成的絲錐,加熱過程采用在400℃、800℃、1200℃三級(jí)加熱,每個(gè)溫度保溫20min,然后以90S的時(shí)間用鹽浴冷卻的方式降至室溫。
2 淬火過程計(jì)算機(jī)模擬的理論依據(jù)
該問題屬于瞬態(tài)熱分析問題。與絲錐相接觸的硝鹽水溶液的溫度和換熱系數(shù)是已知的,所以該問題屬于第三類邊界條件的問題。第三類邊界條件是指與物體相接觸的流體介質(zhì)的溫度和換熱系數(shù)a均為已知,這樣物體與流體介質(zhì)之間的換熱量就為已知。用公式表示為
其中:和α可以是常數(shù),也可以是隨時(shí)間和空間變化的函數(shù)。
又由于材料溫度將在一個(gè)較大的溫度區(qū)間變化,材料熱性能隨溫度變化的非線性因素不能忽略,因此該問題屬于比較復(fù)雜的瞬態(tài)非線性熱分析問題。
3 絲錐模型的創(chuàng)建以及有限元模型的導(dǎo)入
3.1 原始模型的建立
精確建模是正確分析的前提。ANSYS軟件在創(chuàng)建絲錐表面的螺紋、容屑槽等結(jié)構(gòu)要素時(shí)非常麻煩,甚至不可能。為此采用當(dāng)前流行的三維實(shí)體造型軟件Solidworks2006,按照GB967-83的標(biāo)準(zhǔn)精確建模,由此產(chǎn)生的實(shí)體能夠保證足夠的模擬精度。創(chuàng)建出的模型如圖1所示。
圖2 M24絲錐有限元分析模型
3.2 ANSYS分析模型的建立
在絲錐實(shí)際淬火過程中,出現(xiàn)的廢品多緣于切削刃尖角處應(yīng)力過大造成開裂。所以模擬時(shí),關(guān)注的是含有螺紋部分的切削刃。為了抓主要矛盾,合理利用計(jì)算機(jī)資源,提高計(jì)算速度,所以選取切削刃部分1/2長(zhǎng)度,同時(shí)截取徑向的1/4進(jìn)行分析,截取的時(shí)候保證容屑槽的對(duì)稱性,這樣,可以將分析部位的網(wǎng)格劃分的更細(xì),既可以提高計(jì)算精度,又不至于使計(jì)算負(fù)荷過重,同時(shí)還能方便地利用ANSYS的后處理功能,從軸向、徑向?qū)τ诟鞣N場(chǎng)量的分布進(jìn)行直觀的觀察判斷。ANSYS模型如圖2所示。
4 模擬過程中相關(guān)參數(shù)的設(shè)置
在對(duì)溫度場(chǎng)求解的過程中,必要的熱物性參數(shù)是導(dǎo)熱系數(shù)、比熱、對(duì)流換熱系數(shù)、以及材料的密度。其中,因?yàn)槊芏入S溫度變化較小,模擬時(shí)采用恒定值來近似,其余參數(shù)則根據(jù)熱處理手冊(cè)獲取。根據(jù)需要,對(duì)于起始和終止溫度點(diǎn)處的熱物性參數(shù),使用Origin7.5軟件進(jìn)行必要的插值后獲取。
在進(jìn)行應(yīng)力場(chǎng)分析的過程中,材料的彈性模量、泊松比以及線膨脹系數(shù)是必要參數(shù)。泊松比隨溫度變化的影響較小,采用恒定值來處理。彈性模量和線膨脹系數(shù)則在已有值的基礎(chǔ)上作插值補(bǔ)充。
5 ANSYS模擬實(shí)現(xiàn)過程
5.1 定義單元類型
模擬時(shí)一次求解完成兩種場(chǎng)的分析,所以必須采用耦合單元。這里采用了ScalarTet98耦合單元,該單元可以很好地用于曲面模型的模擬。
5.2 定義材料性能參數(shù)
依次選擇MainMenu→Processor→MaterialProps→MeterialModels命令,出現(xiàn)DefineMaterialModelBe.havior對(duì)話框。在MaterialModelAvailable列表框中依次輸入淬火溫度區(qū)間需要的導(dǎo)熱系數(shù)、彈性模量、泊松比、線膨脹系數(shù)、比熱容、密度等所有參數(shù)。如圖3、4。
圖3 淬火過程熱導(dǎo)事值輸入
圖4 淬火過程彈性模量和泊松比輸入
模擬過程中,高溫部分比熱變化不大,取值480J/(kg·K),密度取值,線膨脹系數(shù)在已有資料的基礎(chǔ)上,使用Origin7.5軟件外插值后取值。
5.3 加載求解
(1)設(shè)置分析類型
圖5 淬火過程時(shí)闊步長(zhǎng)設(shè)置
分析類型設(shè)置為Transient,為了完成瞬態(tài)分析,打開自動(dòng)時(shí)間積分項(xiàng)。按圖5進(jìn)行時(shí)間步長(zhǎng)的設(shè)置。
(2)定義初始條件
點(diǎn)擊Main Menu→Solution→Define Loads→Apply→Structural→Temperature→Uniform Temp菜單,在[TUNIF]Uniform temperature輸人框中輸入1200,表示冷卻過程開始時(shí)絲錐具有各個(gè)部位一致的均勻初始溫度。
(3)定義載荷
淬火過程中所受到的熱載荷是處于室溫下的淬火油(工藝規(guī)程規(guī)定是硝鹽水溶液)和高溫下的絲錐之間的對(duì)流換熱載荷,因此分析時(shí)要將對(duì)流系數(shù)和冷卻介質(zhì)初始溫度作為對(duì)流載荷施加給絲錐。此處所用到的硝鹽水溶液的對(duì)流換熱系數(shù)值可以從熱處理手冊(cè)獲得,并按照定義數(shù)組參數(shù)的形式輸入到ANSYS中,如圖6所示。載荷設(shè)置完畢選擇所有表面,加載求解。
圖6 定義對(duì)流邊界條件
6 模擬結(jié)果及分析
6.1 查看溫度場(chǎng)結(jié)果
求解結(jié)束后,打開通用后處理器,讀入最后一步計(jì)算結(jié)果,可以得到淬火冷卻過程結(jié)束溫度場(chǎng)分布等值線圖,如圖7所示。
圖7 淬火結(jié)束芯部溫度場(chǎng)等值線圖
從圖中可看到,經(jīng)過90s冷卻后,刀尖點(diǎn)溫度最低為20℃,芯部的溫度相對(duì)比較高。此處假設(shè)介質(zhì)的溫度為20℃左右,所以模擬結(jié)果和實(shí)際比較吻合。打開時(shí)間歷程后處理器,選取圖2中的四個(gè)特征點(diǎn)為變量,得到溫度隨時(shí)間變化的溫度曲線,如圖8。
可以發(fā)現(xiàn),絲錐冷卻過程中,選定各點(diǎn)的冷卻規(guī)律和趨勢(shì)基本是一致的,差別僅僅在于一個(gè)時(shí)間上的推移。
圖8 淬火過程選定點(diǎn)溫度變化歷程
6.2 查看應(yīng)力場(chǎng)結(jié)果
打開通用后處理器,選擇查看應(yīng)力項(xiàng),得到等效應(yīng)力等值線圖,如圖9所示。
圖9 淬火結(jié)束等效應(yīng)力場(chǎng)等值線圖
圖10 淬火過程牙底最大等效應(yīng)力點(diǎn)應(yīng)力變化曲線
從圖9中可以發(fā)現(xiàn),等效應(yīng)力最大的部位出現(xiàn)在切削刃槽底部位和中心部位,選取淬火后等效應(yīng)力最大的一個(gè)牙底處作為考察點(diǎn)(模型中的1115單元以及該單元包含的178節(jié)點(diǎn)區(qū)域),得到圖10所示的牙底淬火過程中的應(yīng)力變化曲線。
從圖10中可以看出,隨著淬火時(shí)間的延長(zhǎng),牙底部位的應(yīng)力值在開始階段出現(xiàn)波動(dòng),40S以后開始出現(xiàn)快速上升,直到淬火結(jié)束時(shí)等效應(yīng)力值達(dá)到最大。這一結(jié)果和實(shí)際生產(chǎn)中的結(jié)論也是比較吻合的。
7 結(jié)論
(1)通過上述模擬過程可以看出:使用ANSYS軟件對(duì)絲錐淬火冷卻過程的模擬是可行的,無論從模擬所得的溫度場(chǎng)分布云圖,還是應(yīng)力場(chǎng)分布云圖來看,基本上和實(shí)際(該課題有相應(yīng)的試驗(yàn)測(cè)定了溫度場(chǎng))比較吻合,由此可以推廣到別的熱處理過程甚至別的切削工具的熱處理模擬中去,只要能夠測(cè)定出材料的熱物性參數(shù),即可獲得比較滿意的結(jié)果;
(2)ANSYS軟件可以比較直觀地把淬火冷卻過程中的溫度、應(yīng)力變化情況顯示出來,可以為熱處理信息化實(shí)現(xiàn)過程中合理選擇熱處理參數(shù)提供參考,從而為推動(dòng)熱處理信息化的發(fā)展創(chuàng)造條件;
(3)為了驗(yàn)證模擬的可靠性,先前做過相應(yīng)的試驗(yàn),但是由于應(yīng)力場(chǎng)測(cè)量難度較大,存在一定的誤差,而且模擬過程中做了一定的簡(jiǎn)化,比如:沒有考慮組織轉(zhuǎn)變引起的相變潛熱的影響,也會(huì)產(chǎn)生一定的誤差。隨著后續(xù)的進(jìn)一步努力,必將在熱處理計(jì)算機(jī)模擬方面有所突破。
相關(guān)標(biāo)簽搜索:ANSYS的絲錐淬火冷卻過程計(jì)算機(jī)模擬 Fluent、CFX流體分析 HFSS電磁分析 Ansys培訓(xùn) Abaqus培訓(xùn) Autoform培訓(xùn) 有限元培訓(xùn) Solidworks培訓(xùn) UG模具培訓(xùn) PROE培訓(xùn) 運(yùn)動(dòng)仿真