王果，黃嘉敏，陳振明，趙海東*

低壓鑄造鋁合金輪轂?zāi)＞邿嶙冃蔚臄?shù)值模擬

王果1，黃嘉敏1，陳振明2，趙海東1*

（1.華南理工大學(xué) 機械與汽車工程學(xué)院，廣州 510641； 2.佛山南海奔達模具有限公司，廣東 佛山 528200）

針對模具高溫變形影響低壓鑄造鋁合金輪轂尺寸精度的問題，通過數(shù)值模擬和實驗相結(jié)合的方式，揭示低壓鑄造過程中模具的變形規(guī)律。建立了低壓鑄造過程中描述模具熱力學(xué)行為的數(shù)學(xué)模型以及鑄造多循環(huán)計算方法，進行了輪轂低壓鑄造實時藍光掃描驗證，通過模擬與實驗結(jié)果的對比驗證，研究了輪轂低壓鑄造過程中的模具溫度、應(yīng)力及變形規(guī)律。模擬與實測溫度對比結(jié)果表明，兩者變化趨勢基本一致，側(cè)模和下模在開模時刻的最高溫度分別為486 ℃和512 ℃。側(cè)模的模擬與實測溫度峰值的相對誤差為2%，下模模擬與實測溫度最大相差23 ℃。變形模擬與實時掃描結(jié)果均表明，側(cè)模熱變形呈曲率減小趨勢，按變形大小劃分，側(cè)模型腔表面存在小變形（?0.04～0.16 mm）、過渡（0.16～0.48 mm）和大變形（0.48～ 0.89 mm）3個變形區(qū)域，側(cè)模變形的模擬精度約80%；下模平面沿輪轂軸向下沉0.4 mm，變形模擬精度約75%。所建立的熱變形模擬方法可以較好地計算鋁合金鑄造模具熱變形，闡明鋁輪轂低壓鑄造模具變形規(guī)律，為后續(xù)模具反變形設(shè)計、提高輪轂尺寸精度提供研究基礎(chǔ)。

熱變形；模具；低壓鑄造；數(shù)值模擬；藍光實時掃描

低壓鑄造（LPDC）是鋁合金輪轂的主要成形方式[1-3]。隨著輪轂尺寸的增大，模具尺寸增大，同時，為實現(xiàn)快速冷卻和細化組織，通常采用冷卻水道的方法，使模具溫度分布不均勻性加大，這會加劇模具變形，影響鋁合金輪轂的尺寸精度，該變形和尺寸偏差增加了后續(xù)加工量并浪費了材料[4-6]。因此，研究輪轂低壓鑄造模具的熱變形規(guī)律，對獲得尺寸精確的輪轂具有重要意義。

鑄造是涉及流動、凝固、傳熱、應(yīng)力-應(yīng)變等多物理場的復(fù)雜過程[7-11]，鑄件/模具的傳熱和應(yīng)力-應(yīng)變模擬引起了研究者們的高度關(guān)注。Yoon等[12]完成了AC7A和AC4C鋁合金輪轂鑄造過程的數(shù)值模擬，結(jié)果表明，AC7A和AC4C鋁合金輪轂鑄件具有相似的溫度分布，但AC7A輪轂具有更小的變形和應(yīng)力。Song等[13]結(jié)合殼體變形、模具流膜、氣隙動態(tài)分布等特點，建立了寬厚板坯連鑄殼體凝固的二維分切-移動瞬態(tài)熱機械耦合模型，得到了氣隙和殼體變形的分布規(guī)律。Jayakrishna等[14]利用非線性運動學(xué)和各向異性塑性模型，結(jié)合蠕變模型，建立了一個三維熱機模型，用于計算漏斗形銅模和冷卻水腔組合中的瞬態(tài)循環(huán)變形、殘余應(yīng)力和塑性破壞應(yīng)變。Anglada等[15]提出了一種基于熱機械模擬對壓鑄件最終尺寸進行預(yù)測的方法，并使用3D掃描以點云的形式獲取了室溫下實際鑄件的尺寸，完成了對鑄件變形預(yù)測的驗證。綜上所述，盡管研究者們在鑄造模具和鑄件變形的數(shù)值計算方面開展了很多工作，但是因為對模具高溫變形的實時測量存在困難，關(guān)于鑄造中高溫模具變形預(yù)測的驗證研究較少。

本文針對鋁合金輪轂低壓鑄造過程，建立了描述其熱力學(xué)行為的數(shù)學(xué)模型及鑄造多循環(huán)過程的計算方法，完成了低壓鑄造過程中模具變形的實時藍光掃描，并對模擬結(jié)果進行了驗證分析，討論了鋁輪轂低壓鑄造模具的變形規(guī)律。

1 實驗

1.1 低壓鑄造

輪轂材料為A356.2（AlSi7Mg）鑄造鋁合金，其化學(xué)成分如表1所示。LPDC輪轂?zāi)＞呤疽鈭D如圖1a所示。下模和上模采用H13鋼制造，側(cè)模采用42CrMo鋼制造。側(cè)模由4個模具組成，這4個模具在鑄造時裝配成完整的輪轂圓周結(jié)構(gòu)。使用K型熱電偶（直徑0.5 mm）和CoMo Injection高精度數(shù)據(jù)采集系統(tǒng)（瑞士Kistler公司）測量下模和側(cè)模采樣點的溫度變化，其位置如圖1b所示。在每個采樣點位置都設(shè)置一組相鄰的3個熱電偶，熱電偶通過螺栓固定在模具上。熱電偶測溫位置距每個界面（鑄件/下模、鑄件/側(cè)模、鑄件/上模）的距離為2 mm。采用ATOS Ⅲ Triple Scan 8M（德國GOM GmbH）藍光掃描設(shè)備對鑄造開模后的模具進行實時掃描檢測，檢測分辨率為0.01 mm，并利用Free Gom inspect軟件對高溫下的掃描尺寸進行分析，以獲取低壓鑄造過程中模具的實時變形。

1.2 42CrMo鋼高溫拉伸

鋁合金輪轂?zāi)＞叩墓ぷ鳒囟确秶鸀?00～550 ℃，不同溫度下的應(yīng)力-應(yīng)變關(guān)系對模具變形至關(guān)重要。針對實驗輪轂?zāi)＞卟牧?，采用Qayyum等[16]研究中列出的數(shù)據(jù)確定H13鋼在不同溫度下的本構(gòu)關(guān)系，通過不同溫度下的拉伸實驗獲得側(cè)模42CrMo鋼的本構(gòu)關(guān)系。使用Gleeble-1500熱/力模擬試驗機在75、175、275、375、475、575 ℃溫度下進行拉伸實驗。拉伸試樣直徑為6 mm，標距段長度為42 mm。拉伸試樣兩端帶有連接螺紋，其形狀和尺寸如圖1d和圖1e所示，符合GB/T 4338—2006[17]。拉伸應(yīng)變率為2.5×10?3s?1，總變形量控制在5%左右，在上述的每個溫度下重復(fù)進行3次以上的拉伸實驗以保證拉伸數(shù)據(jù)的準確性和可重復(fù)性。

表1 A356.2鋁合金輪轂鑄件的化學(xué)成分

2 數(shù)值模擬

2.1 溫度場模擬

多周期LPDC鑄造包括充型、凝固和開模等基本循環(huán)過程。采用直接有限差分法（FDM）對充型過程進行數(shù)值模擬。充型過程中的動量守恒、質(zhì)量守恒和能量守恒公式如式（1）～（3）所示[18]。

式中：為溫度，℃；為時間，s；()為隨溫度變化的材料密度，kg/m3；()為隨溫度變化的比熱容，J/(kg·℃)；()為隨溫度變化的導(dǎo)熱系數(shù)，W/(m·℃)；為凝固潛熱，W/m3；s為固相分數(shù)，%。

在解決模具開模過程中的傳熱問題時，涉及模腔與空氣的對流換熱以及模腔與環(huán)境的輻射換熱。為了便于描述開模過程中模腔與環(huán)境之間的傳熱行為，引入一個綜合的界面換熱系數(shù)（IHTC）來解決傳熱問題，如式（5）所示。

式中：為熱流密度，W/m2；e為環(huán)境溫度，℃；mold為模具型腔溫度，℃。

在鑄造軟件ICAST Beta版（開發(fā)版）中進行溫度場計算。實際的LPDC是連續(xù)的多模次過程，模具的溫度經(jīng)過若干模次后才能穩(wěn)定，因此需要進行多模次循環(huán)計算，以獲得穩(wěn)定的模具循環(huán)溫度場。由于模擬充型過程所需的計算時間較長，但實際充型過程所需時間較短，在循環(huán)過程中模具溫度變化對鋁液充型后溫度分布的影響不大，因此，提出了如圖2所示的鋁合金LPDC循環(huán)過程模擬計算方法。在第一模次，計算“充型-凝固-開?！?個階段完整的溫度場，從第二模次開始，僅計算凝固和開模2個階段的溫度場直至達到設(shè)定的目標模次。在第一模次計算時，記錄充型結(jié)束后鑄件單元的溫度，在后續(xù)循環(huán)凝固計算時更新模具溫度數(shù)據(jù)，將鑄件單元的溫度設(shè)置為第一模次充型結(jié)束后的溫度。

圖1 輪轂LPDC實驗和42CrMo高溫拉伸

圖2 低壓鑄造多模次循環(huán)計算示意圖

在預(yù)處理模塊中設(shè)置網(wǎng)格劃分、邊界條件和計算參數(shù)。網(wǎng)格劃分后模型的六面體網(wǎng)格總數(shù)為892 800，側(cè)模節(jié)點數(shù)為47 899，下模節(jié)點數(shù)為57 468。溫度場模擬中使用的物理參數(shù)如表2所示。在溫度場計算過程中，每5 s輸出側(cè)模和下模的溫度場數(shù)據(jù)。

2.2 應(yīng)力-應(yīng)變模擬

模具熱變形模擬是基于熱彈塑模型的。根據(jù)塑性增量理論，假定材料應(yīng)變在時間步長中的增量如式（6）所示[15,21-22]。

式中：d為總應(yīng)變增量；del為符合廣義胡克定律的彈性應(yīng)變增量；dpl為符合流動準則的塑性應(yīng)變增量，如式（7）所示，其對應(yīng)的應(yīng)力符合Von Mises屈服準則；dth為由溫度載荷引起的熱應(yīng)變增量，如式（8）所示。

式中：pl為比例系數(shù)；為塑性勢；為應(yīng)力；為溫度的等效熱膨脹系數(shù)。

基于有限元軟件ANSYS，使用溫度場結(jié)果作為載荷計算模具的應(yīng)力和變形。使用反距離加權(quán)（IDW）插值算法[23-25]實現(xiàn)FDM到FEM的溫度場轉(zhuǎn)換。IDW插值算法的基本思想是假設(shè)未采樣節(jié)點的屬性值是鄰域已知值的加權(quán)平均，詳細算法請見文獻[26]。

42CrMo鋼不同溫度下的真應(yīng)力-應(yīng)變曲線如圖3a所示。真應(yīng)力-應(yīng)變曲線被擬合并簡化為2個部分，如圖3b所示，即彈性部分和塑性部分，用于模具熱變形模擬。下模H13鋼的應(yīng)力-應(yīng)變關(guān)系采用文獻[16]中的數(shù)據(jù)，該曲線也被擬合并簡化為彈性和塑性兩部分。熱變形計算中使用的材料性能參數(shù)如表3所示。

表2 材料的熱物性參數(shù),

Tab.2 Physical properties of materials

圖3 42CrMo鋼的真應(yīng)力-應(yīng)變曲線（a）及其擬合簡化后的彈塑性曲線（b）

為了適應(yīng)模具復(fù)雜的幾何結(jié)構(gòu)，保證計算精度，選擇二階四面體元素（SOLID187）進行網(wǎng)格劃分，劃分后側(cè)模節(jié)點數(shù)為86 833，下模節(jié)點數(shù)為184 494。側(cè)模和下模的底面在軸方向的位移被設(shè)置為0；在常溫下側(cè)模兩側(cè)的45°平面留有空隙，所以這個面沒有設(shè)置邊界約束條件。由于鋁合金在凝固過程中的彈性模量比模具鋼的彈性模量小得多，所以模具與鑄件的接觸面被設(shè)定為可以自由變形。

表3 熱變形模擬中模具的材料物理參數(shù)。

Tab.3 Material physical parameters of mold used in thermal deformation simulation

3 結(jié)果與討論

3.1 溫度場

LPDC輪轂鑄造過程中鑄件和模具的溫度場計算結(jié)果如圖4所示。隨著A356.2鋁液充填型腔，鑄件的溫度梯度逐漸增大，當充填結(jié)束（13.2 s）時，鑄件的最低溫度為652.5 ℃。在凝固過程中，輪輞上緣遠離底部中央澆口，降溫速度較快，整體上鑄件呈現(xiàn)出自上而下的凝固順序。在開模時刻（265 s），輪轂鑄件上最高溫度為521 ℃，所有節(jié)點溫度均低于固相線溫度，鑄件已完全凝固。如圖4d和圖4e所示，側(cè)模和下模在開模時刻的最高溫度分別為486 ℃和512 ℃。

N1～N3三個采樣點在10個模次內(nèi)隨時間變化的計算溫度曲線如圖5a所示。隨著模次的增加，模具溫度峰值持續(xù)上升，并在5個模次之后，循環(huán)開始和循環(huán)結(jié)束時的模具溫度基本相同，證明溫度開始穩(wěn)定。選取圖5a中溫度穩(wěn)定后第8模次的N1和N2的模擬溫度和實測溫度進行對比，結(jié)果如圖5b和圖5c所示?？梢钥吹?，模擬溫度與實測溫度（Exp）的變化趨勢基本一致。對于鑄件/側(cè)模界面附近（N1），模擬與實測溫度峰值的相對誤差為2%，由于鑄造過程影響因素眾多，模擬結(jié)果與實測誤差低于3%通常被認為是吻合較好[27]。但是，當單個鑄造過程結(jié)束時，下模N2的溫度與實測溫度相差23 ℃，其原因可能為下模安裝在模具底板上，在實際鑄造過程中，下模向底板傳熱，而底板向低壓鑄造機框架傳熱，而上述傳熱在模擬過程中未被考慮，所以，模擬結(jié)果高估了該點溫度。

圖4 溫度場模擬結(jié)果

3.2 模具變形計算結(jié)果與驗證

250 s（開模前15 s）、265 s（開模時間）和296 s（合模前）時刻的模具應(yīng)力分布情況如圖6所示，可以看到，三者應(yīng)力分布相似。開模后，模腔溫度降低，各部分的溫差也逐漸減小，模具逐漸均勻收縮，整體應(yīng)力變小。側(cè)模的最大應(yīng)力出現(xiàn)在水冷鑲件與側(cè)模接觸的尖角處，這是因為該處溫差大，結(jié)構(gòu)剛度低。下模凸臺和凸臺外沿位置的應(yīng)力較小，凸臺間的應(yīng)力較大。

側(cè)模在250 s（凝固過程）、265 s（開模時刻）和295 s（合模時刻）時的變形模擬結(jié)果分別如圖7a～c所示?？梢钥吹剑?個時刻的變形規(guī)律相似，均表現(xiàn)為型腔表面中部變形較小、表面邊緣變形較大，且兩側(cè)變形呈近似對稱分布的趨勢，三者變形量分別為0.098～1.28、0.15～1.20、0.096～1.12 mm。圖7c中白線是側(cè)模原始幾何模型的外輪廓線，其中在橢圓線框圈出的側(cè)模邊緣位置發(fā)生了明顯的向外側(cè)偏移變形。結(jié)合圖7d中點1～3的位置和圖7e中3個點隨時間變化的變形曲線可知，在整個鑄造過程中，點3、點2和點1的變形逐漸減小，即側(cè)模表現(xiàn)出越靠近邊緣部位變形越大、曲率變小的趨勢，成形的輪轂相應(yīng)部位的厚度將增大。

通過笛卡爾坐標向極坐標的轉(zhuǎn)換獲得了輪轂沿徑向各節(jié)點的模擬變形量，如圖8a所示。采用高精度藍光實時掃描獲得的模具實際變形結(jié)果如圖8b所示，其中變形方向為各表面的法線方向。根據(jù)圖8a中的徑向變形模擬結(jié)果，型腔表面被劃分為3個變形區(qū)：小變形區(qū)域I位于型腔中心；Ⅱ區(qū)域為過渡區(qū)；大變形區(qū)域Ⅲ位于型腔邊緣附近。一方面，型腔表面區(qū)域I的溫度最高（見圖4d），熱膨脹效應(yīng)導(dǎo)致區(qū)域I沿徑向偏離輪轂中心（向外），這使得型腔變大；另一方面，因整體變形呈曲率減小的規(guī)律，所以在變形協(xié)調(diào)作用下，輪轂區(qū)域Ⅲ的向外側(cè)變形會使區(qū)域I向內(nèi)變形，最終導(dǎo)致區(qū)域I有著最小的徑向變形，變形量為?0.04～0.16 mm，這是熱膨脹和變形協(xié)調(diào)共同作用的結(jié)果。區(qū)域Ⅲ表面的徑向變形最大，為0.48～ 0.89 mm（向外）。而區(qū)域Ⅱ為型腔表面區(qū)域I和區(qū)域Ⅲ之間的變形過渡區(qū)，變形量為0.16～0.48 mm（向外）。對應(yīng)圖8b中實時掃描結(jié)果，區(qū)域I的變形值為?0.02～ 0.23 mm，區(qū)域Ⅱ變形值為0.23～0.372 mm，邊緣區(qū)域Ⅲ變形值為0.372～0.714 mm。這表明側(cè)模型腔的實際變形也表現(xiàn)為偏離輪轂中心、曲率減小的趨勢，與模擬結(jié)果一致。兩者定量對比結(jié)果表明熱變形模擬精度約80%，與文獻[15]報道相當。因此，本文所建立的模型與計算方法可以對鋁合金鑄造熱變形進行有效計算。

265 s開模時刻下模方向模擬的變形分布情況如圖9a和圖9b所示。由于265 s開模時，藍光掃描測量被側(cè)模阻擋（見圖1c），僅得到了大部分下模型腔表面的變形結(jié)果，如圖9c所示。模擬結(jié)果顯示，凸臺上表面相對于下模平面向方向的正變形超過0.56 mm。藍光掃描結(jié)果顯示，下模平面相對于凸臺上表面下沉了0.4 mm，兩者對比表明，變形模擬精度為75%。圖9b顯示，凸臺的斜側(cè)面朝向輪轂軸心（朝內(nèi)）的一側(cè)變形較小，朝外的一側(cè)變形大。圖9c顯示，凸臺側(cè)面的實際變形與模擬的變形趨勢相同。

圖7 側(cè)模變形的計算結(jié)果

圖8 側(cè)模的變形驗證

圖9 下模的變形計算結(jié)果與驗證

4 結(jié)論

建立了能在低壓鑄造過程中描述鋁合金模具熱力學(xué)行為的數(shù)學(xué)模型，提出了鑄造多循環(huán)溫度場計算方法，對實際低壓鑄造鋁合金輪轂進行了模擬應(yīng)用；完成了低壓鑄造過程中模具變形的實時藍光掃描，對溫度和變形模擬結(jié)果進行了驗證分析。結(jié)果表明，所建立的方法可以有效地計算鋁合金鑄造熱變形。低壓鑄造輪轂?zāi)＞咦冃文M和實驗結(jié)果表明，越靠近邊緣部位，側(cè)模變形越大，表現(xiàn)出曲率減小的趨勢，所成形的輪轂對應(yīng)側(cè)模邊緣的部位將增厚；下模平面下沉，成形的輪幅厚度增大。

[1] 樓華山, 王大紅, 樓江燕, 等. 高性能A356鋁合金輪轂開發(fā)技術(shù)研究進展之二——輪轂成形強化技術(shù)[J]. 特種鑄造及有色合金, 2020, 40(6): 619-622.LOU Hua-shan, WANG Da-hong, LOU Jiang-yan, et al. Progress in Research on Developing Automobile A356 Aluminum Alloy Wheel-enhanced Treatment of Wheel Forming[J]. Special Casting & Nonferrous Alloys, 2020, 40(6): 619-622.

[2] 秦鵬. 多輪輻鋁合金輪轂低壓鑄造有限元模擬及性能測試[D]. 秦皇島: 燕山大學(xué), 2018: 1-6. QIN Peng. Die Design and the Performance Test of Many Spokes Aluminum Alloy Hub Forming by Low Pressure Casting[D]. Qinhuangdao: Yanshan University, 2018: 1-6.

[3] 胡孟達. 鋁合金輪轂強力鑄造工藝研究[D]. 秦皇島: 燕山大學(xué), 2017: 1-5. HU Meng-da. Research on the Power Casting Technology of Aluminum Alloy Wheel[D]. Qinhuangdao: Yanshan University, 2018: 1-5.

[4] CAMPBELL J. Complete Casting Handbook: Metal Casting Processes, Metallurgy, Techniques and Design[M]. Oxford: Butterworth-Heinemann, 2015.

[5] KUMAR P, SINGH R, AHUJA I P S. Investigations on Dimensional Accuracy of the Components Prepared by Hybrid Investment Casting[J]. Journal of Manufacturing Processes, 2015, 20: 525-533.

[6] CAMPATELLI G, SCIPPA A. A Heuristic Approach to Meet Geometric Tolerance in High Pressure Die Casting[J]. Simulation Modelling Practice and Theory, 2012, 22: 109-122.

[7] KOPILA P, BARMAN M, BABU P M. Simulation of Liquid & Solid Metal Regions in a Continuous Casting Process Using COMSOL Multiphysics[J/OL]. Materials Today: Proceedings, 2023, [2023-01-14]. https://link-inghub.elsevier.com/retrieve/pii/S2214785323024914.

[8] JIANG J, NING G E, HUANG M, et al. Numerical Simulation of Squeeze Casting of Aluminum Alloy Flywheel Housing with Large Wall Thickness Difference and Complex Shape[J]. Transactions of Nonferrous Metals Society of China, 2023, 33(5): 1345-1360.

[9] DONG G, LI S, MA S, et al. Process Optimization of A356 Aluminum Alloy Wheel Hub Fabricated by Low-Pressure Die Casting with Simulation and Experimental Coupling Methods[J]. Journal of Materials Research and Technology, 2023, 24: 3118-3132.

[10] WANG J, ZUO Y, ZHU Q, et al. Numerical Simulation of the Flow Field and Temperature Field in Direct-Chill Casting of 2024 Aluminum Alloy under the Influence of Intensive Melt Shearing[J]. Journal of Materials Research and Technology, 2023, 22: 1491-1504.

[11] ?ABíK V, FUTá? P, PRIBULOVá A. Failure Analysis of a Clutch Wheel for Wind Turbines with the Use of Casting Process Simulation[J]. Engineering Failure Analysis, 2022, 135: 106159.

[12] YOON H S, CHOI J S, OH Y K. Thermal Deformation of Aluminum Alloy Casting Materials for Tire Mold by Numerical Analysis[J]. Transactions of Nonferrous Metals Society of China, 2012, 22: s859-s863.

[13] SONG J, CAI Z, PIAO F, et al. Heat Transfer and Deformation Behavior of Shell Solidification in Wide and Thick Slab Continuous Casting Mold[J]. Journal of Iron and Steel Research, International, 2014, 21: 1-9.

[14] JAYAKRISHNA P, CHAKRABORTY S, GANGULY S, et al. Computational Investigation of the Transient Cyclic Thermal Distortion of Funnel Shaped Mould in Thin Slab Continuous Casting Process[J]. Thermal Science and Engineering Progress, 2022, 36: 101508.

[15] ANGLADA E, MELéNDEZ A, VICARIO I, et al. Prediction and Validation of Shape Distortions in the Simulation of High Pressure Die Casting[J]. Journal of Manufacturing Processes, 2018, 33: 228-237.

[16] QAYYUM F, KAMRAN A, ALI A, et al. 3D Numerical Simulation of Thermal Fatigue Damage in Wedge Specimen of AISI H13 Tool Steel[J]. Engineering Fracture Mechanics, 2017, 180: 240-253.

[17] GB/T 4338—2006, 金屬材料高溫拉伸試驗方法[S].GB/T 4338-2006, High-temperature Tensile Test Method for Metallic Materials[S].

[18] ZHAO H D, WANG F, LI Y Y, et al. Experimental and Numerical Analysis of Gas Entrapment Defects in Plate ADC12 Die Castings[J]. Journal of Materials Processing Technology, 2009, 209(9): 4537-4542.

[19] KERMANPUR A, MAHMOUDi S, HAJIPOUR A. Numerical Simulation of Metal Flow and Solidification in the Multi-Cavity Casting Moulds of Automotive Components[J]. Journal of Materials Processing Technology, 2008, 206(1/2/3): 62-68.

[20] PRABHAKAR A, PAPANIKOLAOU M, SALONITIS K, et al. Sand Casting of Sheet Lead: Numerical Simulation of Metal Flow and Solidification[J]. The International Journal of Advanced Manufacturing Technology, 2020, 106(1): 177-189.

[21] LEWIS R W, HAN Z Q, GETHIN D T. Three- Dimension Finite Element Model for Metal Displacement and Heat Transfer in Squeeze Casting Processes[J]. Comptes Rendus Mecanique, 2017, 335(5): 287-294.

[22] POSTEK E W, LEWIS R W, GETHIN D T. Finite Element Modeling of the Squeeze Casting Process[J]. International Journal of Numerical Method for Heat & Fluid Flow, 2008, 18(3/4): 325-355.

[23] LU G Y, WONG D W. An Adaptive Inverse-Distance Weighting Spatial Interpolation Technique[J]. Computers and Geosciences, 2008, 34(9): 1044-1055

[24] WITTEVEEN A S. Explicit and Robust Inverse Distance Weighting Mesh Deformation for CFD[J]. European Neuropsychopharmacology, 2013, 6(2): 98-98

[25] SI H M, CHO C, KWAHK S Y. A Hybrid Method for Casting Process Simulation by Combining FDM and FEM with an Efficient Data Conversion Algorithm[J]. Journal of Materials Processing Technology, 2003, 133(3): 311-321

[26] 王子超. 鋁合金壓鑄模具熱循環(huán)、應(yīng)力與疲勞的模擬研究[D]. 廣州: 華南理工大學(xué), 2021: 42-43. WANG Zi-chao. Simulation Study on Thermal Cycle, Stress and Fatigue of Aluminum Alloy Die-Casting Dies[D]. Guangzhou: South China University of Science and Technology, 2021: 42-43.

[27] LIN J H, ZHAO H D, HUANG J M. Spatial Interfacial Heat Transfer and Surface Characteristics during Gravity Casting of A356 Alloy[J]. Transactions of Nonferrous Metals Society of China, 2019, 29(1): 43-50.

Numerical Simulation of Thermal-deformation of Molds for Low Pressure Die-casting Aluminum Alloy Wheel Hubs

WANG Guo1,HUANG Jia-min1,CHEN Zhen-ming2,ZHAO Hai-dong1*

(1. School of Mechanical and Automotive Engineering, South China University of Technology, Guangzhou 510641, China; 2. Foshan Nanhai Superband Mould Co., Ltd., Guangdong Foshan 528200, China)

The work aims to reveal the rule of mold deformation during low pressure die casting process with numerical simulation and experiment to solve the effect of high-temperature deformation of the mold on the dimensional accuracy of low pressure die casting (LPDC) aluminum alloy wheel hubs. A numerical modelon thermodynamic behavior of the mold during low pressure die casting was established, and a simulation method for casting cycles was proposed. The real-time blue laser inspection verification of low pressure wheel hub casting was carried out. The mold temperature, stress and deformation rule during low-pressure wheel hub casting were studied by comparing the simulation and experimental results. The simulated and measured temperature trends were consistent with each other. The maximum temperatures of the lateral mold and the bottom mold at the moment of mold opening were 486 ℃ and 512 ℃, respectively. The relative error between the simulated and measured peak temperatures for the lateral mold was 2%, and the maximum difference between the simulated and measured temperatures for the bottom mold was 23 ℃. The simulation and inspection results showed that the thermal deformation of lateral mold showed a curvature decreasing trend. There were three deformation regions of small deformation (?0.04-0.16 mm), transition (0.16-0.48 mm) and large deformation (0.48-0.89 mm) on the surface of the cavity of the lateral mold according to the deformation size. The accuracy of deformation simulation for lateral mold was about 80%. The bottom mold plane was sunk 0.4 mm in the axial direction of the hub, and its simulation accuracy was about 75%. The established thermal-deformation simulation method can well describe deformation of aluminum alloy casting molds. The study shows the deformation rule of molds of LPDC aluminum wheel hubs, providing a research basis for the subsequent inverse design of mold deformation to improve dimensional accuracy of wheel hubs.

thermal deformation; mold; low pressure die casting; numerical simulation; real-time inspection of blue laser

10.3969/j.issn.1674-6457.2023.10.011

TG21+3

A

1674-6457(2023)10-0095-09

2023-08-04

2023-08-04

廣東省重點領(lǐng)域研發(fā)計劃（2020B010186002）；佛山市科技計劃（2016AG100932）

Key-Area Research and Development Program of Guangdong Province(2020B010186002); Foshan Science and Technology Program Project (2016AG100932)

王果, 黃嘉敏, 陳振明, 等. 低壓鑄造鋁合金輪轂?zāi)＞邿嶙冃蔚臄?shù)值模擬[J]. 精密成形工程, 2023, 15(10): 95-103.

WANG Guo, HUANG Jia-min, CHEN Zhen-ming, et al. Numerical Simulation of Thermal-deformation of Molds for Low Pressure Die-casting Aluminum Alloy Wheel Hubs[J]. Journal of Netshape Forming Engineering, 2023, 15(10): 95-103.

責任編輯：蔣紅晨