球墨鑄鐵曲軸多維振動(dòng)鑄造充型過程的數(shù)值模擬

張子鵬 崔驕建

摘要：為研究多維振動(dòng)條件下金屬液在型腔中的充型過程，以汽車發(fā)動(dòng)機(jī)用曲軸為研究對(duì)象，運(yùn)用離散元軟件對(duì)曲軸充型過程進(jìn)行數(shù)值模擬分析，并采用正交試驗(yàn)法分析了各振動(dòng)參數(shù)對(duì)充型過程的影響。結(jié)果表明，各振動(dòng)參數(shù)對(duì)充型過程影響的大小順序是振動(dòng)幅度、振動(dòng)自由度、振動(dòng)頻率。當(dāng)振動(dòng)頻率為50Hz、振動(dòng)自由度為3、振動(dòng)幅度為0.5mm時(shí)，曲軸鑄件的成型質(zhì)量較高，為離散元數(shù)值模擬技術(shù)以及多維振動(dòng)在鑄件充型過程中的應(yīng)用提供參考。

關(guān)鍵詞：球墨鑄鐵曲軸;多維振動(dòng);離散單元法

中圖分類號(hào)：TG255? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 文獻(xiàn)標(biāo)識(shí)碼：A? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? ? 文章編號(hào)：1674-957X（2020）21-0004-02

0? 引言

鑄造是獲得機(jī)械零部件的一種快速成型方法，它被廣泛應(yīng)用于工業(yè)生產(chǎn)中。傳統(tǒng)鑄造生產(chǎn)中，通常采用“試錯(cuò)+試制”法來確定合理的工藝，但這種方法制造周期長、生產(chǎn)效率低、生產(chǎn)成本高[1]。研究表明，在金屬液充型過程中施加一定的振動(dòng)可以提高鑄件質(zhì)量，降低鑄件廢品率[2-5]。然而，在現(xiàn)有的機(jī)械振動(dòng)及鑄件充型技術(shù)中，存在振動(dòng)自由度單一、不能實(shí)時(shí)觀察金屬液在型腔中的流動(dòng)狀況和充型過程等問題，使得中大型復(fù)雜鑄件的成型質(zhì)量受到影響。實(shí)驗(yàn)研究固然會(huì)使實(shí)驗(yàn)更具說服力，然而大量實(shí)驗(yàn)勢必會(huì)使得實(shí)驗(yàn)成本增加、結(jié)果論證周期延長以及實(shí)驗(yàn)可觀性變差等后果。為此，安徽理工大學(xué)王成軍等提出了一種基于離散元軟件模擬機(jī)械振動(dòng)對(duì)鑄件充型過程影響的方法，通過設(shè)置振動(dòng)自由度、振動(dòng)頻率和振動(dòng)幅度等參數(shù)，并運(yùn)用細(xì)小顆粒模擬金屬液在鑄件型腔中的流動(dòng)，有效解決傳統(tǒng)鑄造過程以及現(xiàn)有模擬軟件中存在的缺點(diǎn)[6]。本文運(yùn)用離散元數(shù)值模擬軟件，以汽車用曲軸為研究對(duì)象，采用控制變量法，分析了曲軸在多維振動(dòng)條件下的鑄造充型過程，并通過正交試驗(yàn)法分析振動(dòng)自由度、振動(dòng)頻率和振動(dòng)幅度等振動(dòng)參數(shù)對(duì)曲軸充型過程影響程度的大小。

1? 離散單元法的概述及數(shù)值模擬

離散單元法（Discrete Element Method，簡稱：DEM）是由美國Cundall P. A教授在1971年首次提出的一種基于顆粒離散體物料分析的方法。自研究巖土工程和地質(zhì)力學(xué)兩大傳統(tǒng)學(xué)科領(lǐng)域開始，經(jīng)過50多年的研究與深入，現(xiàn)已被專家和學(xué)者用于在礦物加工、土木工程、化學(xué)工程等領(lǐng)域模擬顆粒系統(tǒng)單元粒子流動(dòng)，同時(shí)分析顆粒剪切和填充特性[7]。離散單元法的基本原理是將研究對(duì)象離散成一個(gè)個(gè)剛性元素，然后根據(jù)牛頓第二運(yùn)動(dòng)定律，采用動(dòng)態(tài)、靜態(tài)松弛法等時(shí)間步迭代法求解各元素的運(yùn)動(dòng)學(xué)方程[8-9]。

本文選取汽車發(fā)動(dòng)機(jī)用曲軸為實(shí)驗(yàn)研究對(duì)象，運(yùn)用Solidworks建立曲軸零件三維模型。曲軸質(zhì)量為44kg，外形尺寸為Φ656×190mm，主軸頸零件尺寸為Φ85mm，連桿軸頸尺寸為Φ70mm，平衡塊壁厚為25mm。

1.1 模擬實(shí)驗(yàn)參數(shù)的設(shè)置

為真實(shí)反映金屬液在鑄件型腔中的流動(dòng)，本文結(jié)合離散元數(shù)值模擬的特點(diǎn)，設(shè)置重力加速度大小9.81m/s2，方向+Y;接觸模型為顆粒與鑄件;顆粒半徑0.2mm，其生產(chǎn)方式為固定數(shù)量（顆粒數(shù)量：20000）;顆粒進(jìn)入鑄件型腔內(nèi)的速度為5m/s，方向?yàn)?X;仿真時(shí)間為1s，網(wǎng)格尺寸0.8mm，來研究顆粒在鑄件型腔內(nèi)的充型距離，進(jìn)而間接反映金屬液在鑄件型腔中的充型能力，設(shè)置鑄件型腔和顆粒的材料，具體包括材料的泊松比、密度、剪切模量等參數(shù)的設(shè)置。

1.2 模擬實(shí)驗(yàn)的過程

為減少離散元軟件的仿真時(shí)間，仿真前本文適當(dāng)?shù)缺壤s小仿真模型，并將曲軸模型另存為IGS格式導(dǎo)入離散元軟件中。充型開始時(shí)，模擬金屬液的細(xì)小顆粒從與冒口相連的法蘭端以一定的速度進(jìn)入曲軸鑄件的型腔，并在自身重力作用下下降。隨著顆粒的進(jìn)一步生成，顆粒在型腔內(nèi)流動(dòng)并與型腔內(nèi)壁發(fā)生碰撞，使得顆粒流動(dòng)速度降低并逐漸產(chǎn)生堆積。當(dāng)顆粒運(yùn)動(dòng)到曲軸另一端時(shí)，如果振動(dòng)參數(shù)設(shè)置合理、計(jì)算時(shí)間充足，顆粒將逐漸填滿鑄件整個(gè)型腔，如圖1所示。

1.3 不同振動(dòng)參數(shù)對(duì)曲軸充型能力的影響

研究不同振動(dòng)參數(shù)對(duì)鑄件充型能力的影響，采用控制變量法，分析該因素對(duì)顆粒在型腔中充型距離的影響程度，從而間接表征該因素對(duì)金屬液充型能力的影響。本文將選擇在相同澆注速度和相同模擬總步數(shù)的條件下，分析顆粒在型腔中的充型距離，從而進(jìn)一步用來判斷金屬液在型腔中充型能力的大小。

1.3.1 振動(dòng)自由度的影響

振動(dòng)自由度的類型有單自由度（X、Y、Z），兩自由度（XY、XZ、YZ）和三自由度（XYZ）三種。本文選取f=50Hz，A=0.5mm來研究振動(dòng)自由度對(duì)曲軸充型過程的影響，結(jié)果如表2所示。當(dāng)振動(dòng)自由度DOF=3時(shí)，顆粒充型距離取得最大值，說明顆粒充型距離隨振動(dòng)自由度的增加而增加。

1.3.2 振動(dòng)頻率的影響

由表2可知，當(dāng)自由度DOF=3時(shí)，充型距離取得最大值，故本實(shí)驗(yàn)取DOF=3，A=0.5mm來研究振動(dòng)頻率對(duì)曲軸充型能力的影響，結(jié)果如表3所示。當(dāng)振動(dòng)頻率f=50Hz時(shí)，顆粒充型距離取得最大值。若進(jìn)一步增加振動(dòng)頻率，金屬液產(chǎn)生紊流的可能性增大，反而會(huì)使得金屬液的充型能力降低。

1.3.3 振動(dòng)幅度的影響

由表3可知，當(dāng)DOF=3、f=50Hz時(shí)充型距離取得最大值，故本次實(shí)驗(yàn)取DOF=3、f=50Hz來研究振動(dòng)幅度對(duì)曲軸充型能力的影響。由表4可知，當(dāng)振動(dòng)幅度A=0.5mm時(shí)，顆粒充型距離達(dá)到最大值。

1.4 模擬實(shí)驗(yàn)研究的結(jié)果

為減少實(shí)驗(yàn)次數(shù)，提高實(shí)驗(yàn)數(shù)據(jù)的可信度，本文采用正交試驗(yàn)法來確定振動(dòng)自由度、振動(dòng)頻率、振動(dòng)幅度等因素對(duì)鑄件充型能力的影響程度，令B為振動(dòng)自由度;令C為振動(dòng)頻率;令D為振動(dòng)幅度，實(shí)驗(yàn)方案采用L9（34）的正交試驗(yàn)表，如表5所示，采用正交試驗(yàn)得出最佳的實(shí)驗(yàn)方案為B3C3D1，各參數(shù)對(duì)顆粒充型距離影響大小的順序依次為振動(dòng)頻率、振動(dòng)自由度、振動(dòng)幅度。

2? 結(jié)論

本文采用正交試驗(yàn)法分析各振動(dòng)參數(shù)對(duì)顆粒充型距離的影響，結(jié)果表明，各振動(dòng)參數(shù)對(duì)顆粒在鑄件型腔中充型能力的影響程度從大到小依次為：振動(dòng)頻率>振動(dòng)自由度>振動(dòng)幅度。當(dāng)振動(dòng)頻率為50Hz、振動(dòng)自由度為3、振動(dòng)幅度為0.5mm時(shí)，顆粒在鑄件型腔中的充型距離取得最大值。與傳統(tǒng)鑄造相比，離散單元法不僅真實(shí)模擬金屬液在鑄件型腔中的充型過程，實(shí)時(shí)觀察金屬液在鑄件型腔中的流動(dòng)情況，而且可以分析多維振動(dòng)對(duì)金屬液充型能力的影響。因此，本文研究結(jié)果可為離散單元法在多維振動(dòng)鑄件充型過程中的進(jìn)一步應(yīng)用提供理論依據(jù)。

參考文獻(xiàn)：

[1]傳海軍，劉高峰，楊新強(qiáng)，等.操縱殼體壓鑄充型和凝固過程的數(shù)值模擬[J].鑄造，2015，64（02）：130-134.

[2]吳士平，王汝佳，陳偉，等.振動(dòng)過程的數(shù)值模擬在金屬凝固中應(yīng)用的研究進(jìn)展[J].金屬學(xué)報(bào)，2018，54（02）：247-264.

[3]王成軍，崔驕建.振動(dòng)技術(shù)在金屬鑄造成形過程中的應(yīng)用與研究[J].鑄造技術(shù)，2018，39（10）：2240-2243.

[4]孫勇，朱留憲，吳代建.機(jī)械振動(dòng)對(duì)AZ31合金充型和組織性能的影響[J].鑄造技術(shù)，2017，38（12）：2936-2939.

[5]戴貴鑫.振動(dòng)傳輸對(duì)高強(qiáng)鋁合金低壓鑄造縮松缺陷的影響[D].哈爾濱工業(yè)大學(xué)，2018.

[6]崔驕建.曲軸鑄件鑄造成型過程模擬及工藝研究[D].安徽理工大學(xué)，2019.

[7]胡國明.顆粒系統(tǒng)的離散元素法分析仿真[M].武漢：武漢理工大學(xué)出版社，2010.

[8]王國強(qiáng)，郝萬軍，王繼新.離散元法及其在EDEM上的實(shí)踐[M].西安：西北工業(yè)大學(xué)出版社，2010.

[9]孫其誠，王光謙.顆粒物質(zhì)力學(xué)導(dǎo)論[M].北京：科學(xué)出版社， 2009.