利用熱應(yīng)力場(chǎng)模擬分析對(duì)壓鑄模具局部區(qū)域壽命進(jìn)行估算
孫永明，陳振宇，陳國(guó)恩，李四娣，汪學(xué)陽(yáng)
（廣東鴻圖科技股份有限公司，廣東高要 526108）
摘要：壓鑄模具表面在壓鑄生產(chǎn)過(guò)程中周期性的重復(fù)受到拉應(yīng)力、壓應(yīng)力影響從而導(dǎo)致出現(xiàn)熱疲勞。本文利用MAGMASOFT軟件對(duì)某大型動(dòng)力總成部件模具熱應(yīng)力場(chǎng)進(jìn)行數(shù)值模擬分析，通過(guò)熱疲勞壽命計(jì)算可以估算出模具距離進(jìn)澆口不同距離的部位實(shí)際使用壽命，并針對(duì)使用壽命較低區(qū)域提出設(shè)計(jì)優(yōu)化措施，從而提高整套模具使用壽命。
關(guān)鍵詞：熱應(yīng)力分析；模具壽命；鑲件；
中圖分類號(hào)： 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：
中圖分類號(hào)： 文獻(xiàn)標(biāo)識(shí)碼：A 文章編號(hào)：
The simulation analysis of thermalstress field which is used to estimate the local life of die - casting dies
Yongmingsun,Zhenyu chen,Guoen chen,Sidi li,Xueyangwang
(GuangdongHongtu Technology(Holdings) Co.,Ltd., Guangdong 526108)
Keywords:thermal stress analysis; Die life; mold Insert;
Abstract:The surface of die - casting mould is affected by tensile stress andcompressive stress periodically in the process of die-casting which will leadto thermal fatigue. In thispaper, MAGMASOFT software is used to conduct numerical simulation analysis ofthe thermal stress field of a large powertrain component mold. The actualservice life of the parts with different distance from the mold to the gate canbe estimated by calculating the thermal fatigue life. Finally, designoptimization measures for areas with low servicelife are put forward to improve the service lifeof the whole mold.
壓鑄模具的主要工作環(huán)境是高溫、高壓，因此對(duì)壓鑄模具提出了熱疲勞強(qiáng)度高、耐腐蝕性好等要求。由于壓鑄模具的制作成本高、周期長(zhǎng)，例如類似大型復(fù)雜的動(dòng)力總成模具制作周期往往需要超過(guò)120天，模具材料選用H13情況下，成本超過(guò)100萬(wàn)元，因此模具的使用壽命非常的重要。模具的失效形式主要有：塑性變形、熱疲勞、腐蝕、沖蝕及損傷等。實(shí)際生產(chǎn)過(guò)程中，在多種因素共同作用下影響了模具實(shí)際使用壽命，例如材料的選型、材料熱處理工藝、壓鑄工藝條件、表面強(qiáng)化處理、結(jié)構(gòu)設(shè)計(jì)、模具加工工藝、過(guò)程維護(hù)保養(yǎng)等。
熱疲勞是壓鑄模具早期主要的失效模式。對(duì)我國(guó)目前壓鑄模具失效原因進(jìn)行分析統(tǒng)計(jì)，發(fā)現(xiàn)因熱疲勞導(dǎo)致模具失效占所有模具失效模式近70%。因此，本文主要針對(duì)壓鑄模具生產(chǎn)過(guò)程中受到熱應(yīng)力堆積所造成的熱疲勞失效進(jìn)行研究分析。模具在壓鑄生產(chǎn)過(guò)程中，受到金屬液循環(huán)的高溫高速?zèng)_刷，以致模具表面受到持續(xù)的熱沖擊，模具表面在反復(fù)地加熱、散熱的過(guò)程急熱、急冷的熱量傳遞過(guò)程，模具表面溫度與內(nèi)部溫度變化不同步存在相對(duì)膨脹約束作用，導(dǎo)致模具表面溫度上升時(shí)候受到壓縮應(yīng)力，在溫度下降時(shí)候受到拉伸應(yīng)力。所產(chǎn)生的應(yīng)力可用表達(dá)式（4-1）來(lái)表示：
（1）
式中：-- 模具受到的應(yīng)力，Pa；
-- 溫度梯度，°C/mm；
-- 彈性模量，Pa；
-- 模具的熱膨脹系數(shù)，1/°C；
-- 模具平均溫度，°C。
1實(shí)驗(yàn)方法
本文研究中使用了MAGMASOFT軟件對(duì)模具在壓鑄生產(chǎn)過(guò)程中的溫度場(chǎng)進(jìn)行了數(shù)值模擬分析，隨后將溫度場(chǎng)模擬結(jié)果轉(zhuǎn)化為溫度載荷作為輸入，并利用有限元法對(duì)模具熱應(yīng)力場(chǎng)進(jìn)行數(shù)值模擬分析。
2 模型建立及參數(shù)設(shè)置
2.1分析模型的選擇
根據(jù)分析對(duì)象的不同熱應(yīng)力模擬模型有以下幾種：熱彈性模型、熱粘性模型、熱彈塑性模型及熱彈粘塑性模型，其中運(yùn)用范圍最廣的是熱彈塑性模型。熱彈塑性模型不考慮材料粘性效果，并假設(shè)材料在屈服前變形為彈性變形，屈服后的變形為塑性變形。材料的彈性模量及屈服應(yīng)力只考慮其與溫度計(jì)算過(guò)程中的函數(shù)關(guān)系，且在材料臨近液相溫度時(shí)，對(duì)應(yīng)彈性模量及屈服應(yīng)力均為0。本文研究對(duì)象為H13的壓鑄模具，根據(jù)H13鋼熔點(diǎn)及工作溫度等相關(guān)條件分析，本實(shí)驗(yàn)所進(jìn)行的熱應(yīng)力分析計(jì)算方法應(yīng)選用熱彈性模型[1]。
2.2模擬初始條件
初始條件的假設(shè)：
（1）將熱彈塑性模型中材料應(yīng)力和應(yīng)變非線性曲線簡(jiǎn)化為雙線性，材料在彈性和塑性階段均為線性；
（2）材料屈服滿足Von Mises屈服條件，即金屬通過(guò)累積達(dá)到單位體積內(nèi)的應(yīng)變能，過(guò)渡到塑性狀態(tài)，塑性狀態(tài)下金屬材料的等效應(yīng)力為定值，其大小不受應(yīng)力的狀態(tài)所影響；
（3）材料滿足各向同性硬化法則，材料在進(jìn)入塑性狀態(tài)后，加載屈服面相對(duì)于初始屈服面的中心及形狀均維持不變；其表達(dá)式為：
（2）
式中：k為硬化參數(shù)。
（4）材料滿足塑性流動(dòng)法則，流動(dòng)法則規(guī)定了塑性狀態(tài)下塑性應(yīng)變?cè)隽颗c屈服曲面法向方向的關(guān)系，從而預(yù)測(cè)材料塑性應(yīng)變的大小及擴(kuò)展方向；其表達(dá)式為：
（3）
式中：為系數(shù)；為塑性勢(shì)，σ是應(yīng)力的函數(shù)。
2.3模擬計(jì)算的本構(gòu)方程
（4）
變分方程：
（5）
2.3模擬材料的熱物性設(shè)置
本文研究的壓鑄模具材料為H13，通過(guò)相關(guān)文獻(xiàn)資料查閱，鑄件與模具換熱系數(shù)設(shè)置為1000W/(m2·K)；定義環(huán)境及噴涂劑溫度為20°C條件下，模具與空氣換熱系數(shù)設(shè)置為20W/(m2·K)，模具與噴涂劑換熱系數(shù)設(shè)置為600W/(m2·K)，澆注溫度為660°C，模具預(yù)熱溫度為180°C，冷卻水道參數(shù)按模具3D數(shù)模結(jié)構(gòu)[2]；
3 數(shù)值模擬結(jié)果分析
由式（1）可知，模具表面的應(yīng)力大小與對(duì)應(yīng)的溫度梯度成正比。因此，根據(jù)模具表面溫度場(chǎng)模擬結(jié)果分別在模具定模模面的充填前段、中段及尾段取點(diǎn)進(jìn)行熱應(yīng)力模擬分析，并對(duì)比3個(gè)區(qū)域的模具表面熱應(yīng)力變化曲線進(jìn)行分析，溫度場(chǎng)模擬結(jié)構(gòu)如圖2所示。具體取點(diǎn)位置如圖1所示，A區(qū)域中模具表面分別選取A點(diǎn)與間隔3mm的B點(diǎn)；C區(qū)域中模具表面分別選取C點(diǎn)與間隔3mm的D點(diǎn)；E區(qū)域中模具表面分別選取E點(diǎn)與間隔3mm的F點(diǎn)。
圖1 按距離進(jìn)澆口不同位置取A點(diǎn)、C點(diǎn)、E點(diǎn)
圖2 A點(diǎn)、C點(diǎn)、E點(diǎn)單循環(huán)模具表面溫度變化
圖3為在應(yīng)力數(shù)值模擬結(jié)果基礎(chǔ)上輸出的A區(qū)域、C區(qū)域及E區(qū)域在多個(gè)壓鑄循環(huán)周期內(nèi)的熱應(yīng)力變化曲線圖。從圖1可以看出，A-B區(qū)域、C-D區(qū)域及E-F區(qū)域在壓鑄循環(huán)周期內(nèi)所受到的最大壓縮應(yīng)力出現(xiàn)在合模階段金屬液剛開始充填進(jìn)型腔模具表面急熱的時(shí)候，最大壓縮應(yīng)力分別為126.5MPa、120.1MPa及119.2MPa，最大拉應(yīng)力出現(xiàn)在模具開模后模具表面急冷的時(shí)候，最大拉應(yīng)力分別為26.5MPa、20.1MPa及19.2MPa。由模擬結(jié)果可以看到，模具表面所受到的最大壓應(yīng)力與拉應(yīng)力相差較大。這主要是因?yàn)?，合模期間，模具定模面不僅受到熱膨脹內(nèi)部的約束力還受到了動(dòng)模的鎖模約束力；而在模具打開以后，定模模面只受到收縮內(nèi)部的約束力。圖4為A區(qū)域、C區(qū)域及E區(qū)域在1個(gè)壓鑄循環(huán)周期內(nèi)的模具表面應(yīng)變變化曲線圖。由圖可知，模具3個(gè)區(qū)域內(nèi)最大應(yīng)變分別為0.034mm、0.025mm及0.022mm。同時(shí)發(fā)現(xiàn)模具表面的應(yīng)變最大值出現(xiàn)在模具剛打開的時(shí)候。因?yàn)槎Ｔ谀＞叽蜷_后不再受到動(dòng)模鎖模力的約束，模具表面從壓應(yīng)力轉(zhuǎn)換為拉應(yīng)力過(guò)程是應(yīng)力大量釋放的過(guò)程，過(guò)程中模具變形在屈服區(qū)外時(shí)會(huì)伴隨出現(xiàn)彈性變形，因此循環(huán)周期后期應(yīng)變量會(huì)出現(xiàn)回落[3]。
圖3 模具A、C、E區(qū)域多個(gè)周期內(nèi)表面所受應(yīng)力變化曲線
圖4 模具A、C、E區(qū)域表面應(yīng)變值變化
