孟志軍，胡茂良，牛連杰，王曄

（1. 西安諾博爾稀貴金屬材料有限公司，西安 710201；2. 哈爾濱理工大學(xué) 材料科學(xué)與工程學(xué)院，哈爾濱 150008；3. 河北建筑工程學(xué)院 數(shù)理系，河北 張家口 075000）

Cu-Ag合金是典型的二元共晶體系合金，銅銀合金擁有優(yōu)異的機(jī)械性能、較高的硬度、非常優(yōu)良的導(dǎo)電和導(dǎo)熱性能以及較好的耐磨性。未來(lái)高速列車用接觸線材、強(qiáng)磁場(chǎng)磁體系統(tǒng)、脈沖強(qiáng)磁場(chǎng)線圈以及大規(guī)模集成電路引線框架等方面將會(huì)優(yōu)先選擇潛在導(dǎo)體材料[1—3]，同時(shí)該材料也是電機(jī)換向器和接觸線的理想材料，也常常用作高壓器件和大電流繼電器的接點(diǎn)材料[4—6]。銅銀基合金制備的電接觸材料在電力系統(tǒng)、電器工業(yè)中發(fā)揮著比較重要的作用[7—9]，在焊接過(guò)程中也展現(xiàn)出了優(yōu)良的流動(dòng)性和浸潤(rùn)性，此外銅銀合金也廣泛應(yīng)用于航天工業(yè)中高溫合金和陶瓷之間的焊接釬料[10]，由于其應(yīng)用前景非常廣闊，近年來(lái)已經(jīng)得到了研究者們的廣泛關(guān)注。顯然 Cu-Ag合金的鑄態(tài)顯微組織將決定材料的各方面性能，因此制備出組織成分均勻的鑄件是獲得材料高強(qiáng)高導(dǎo)性能的關(guān)鍵[11—12]，因此，如何快速而又準(zhǔn)確地預(yù)測(cè)Cu-Ag合金中常出現(xiàn)的成分偏析缺陷，已成為鑄態(tài) Cu-Ag合金生產(chǎn)中的重要環(huán)節(jié)。

隨著計(jì)算機(jī)技術(shù)的發(fā)展，鑄件缺陷的預(yù)測(cè)和相關(guān)數(shù)據(jù)的分析已經(jīng)在科學(xué)研究和工業(yè)生產(chǎn)中扮演著越來(lái)越重要的角色。根據(jù)大量數(shù)據(jù)對(duì)鑄件質(zhì)量進(jìn)行準(zhǔn)確和有效預(yù)測(cè)，能夠?yàn)楣に噮?shù)的調(diào)整和控制提供高效、可靠的依據(jù)[13]。此外，作為一種應(yīng)用比較廣泛的統(tǒng)計(jì)模型，多元回歸分析對(duì)于解決鑄造生產(chǎn)過(guò)程中的缺陷預(yù)測(cè)問題來(lái)說(shuō)，是最佳的選擇[14]。與傳統(tǒng)方法相比，該方法具有簡(jiǎn)便、快捷、可靠和準(zhǔn)確的特點(diǎn)，同時(shí)回歸分析能夠建立出表達(dá)各類工藝參數(shù)與最終目標(biāo)相互關(guān)系的函數(shù)關(guān)系式，從而可以直接通過(guò)計(jì)算獲得與真實(shí)值接近的預(yù)測(cè)值，而且不受給定試驗(yàn)樣本的范圍限制[15]，更重要的是，該方法還能夠自動(dòng)略去試驗(yàn)樣本數(shù)據(jù)中那些不穩(wěn)定的數(shù)值[16]。近十年來(lái)，許多學(xué)者采用多元回歸分析方法對(duì)鑄件質(zhì)量和性能進(jìn)行預(yù)測(cè)并優(yōu)化鑄造過(guò)程中的參數(shù)。V. D. Tsoukalas根據(jù)多元線性回歸分析和遺傳運(yùn)算方法提出了一種有效方法，并通過(guò)這種方法針對(duì)壓力鑄造AlSi9Cu3合金中的縮孔問題進(jìn)行了工藝參數(shù)的優(yōu)化[17]。M. Abdel Aziz通過(guò)實(shí)驗(yàn)設(shè)計(jì)方法研究了 AA6063-(TiC/Al2O3)合金的腐蝕特性，并根據(jù)線性回歸分析建立了回歸方程，其實(shí)驗(yàn)結(jié)果與回歸方程的預(yù)測(cè)結(jié)果一致[18]。M. S. Huang和T. Y. Lin基于回歸模型提出了一種穩(wěn)健優(yōu)化方法，該方法能夠有效減少干擾因素對(duì)噴射成形過(guò)程的影響[19]。舒大禹等采用正交試驗(yàn)對(duì)Ti8LC合金的熱處理制度進(jìn)行了設(shè)計(jì)，并對(duì)Ti8LC合金的熱處理工藝與抗拉強(qiáng)度進(jìn)行了回歸分析，建立了熱處理制度與抗拉強(qiáng)度之間的回歸方程[20]。David Hudaka和 Murat Tiryakioglub運(yùn)用線性回歸分析方法和概率估計(jì)量，比較了兩種類型形狀參數(shù)（威布爾模數(shù)）的無(wú)偏估計(jì)[21]。D. P. Mondal等通過(guò)線性回歸分析建立了Al-Zn合金鑄件中外加載荷、磨損尺寸、Zn元素濃度以及各參數(shù)之間的交互作用與磨損速率之間的函數(shù)關(guān)系式[22]。凌成智采用多元回歸擬合配合改進(jìn)的模擬退火算法，進(jìn)行工藝參數(shù)的優(yōu)化，獲得了較好的優(yōu)化效果[23]。

綜上所述，目前回歸分析方法的研究已在生產(chǎn)成本和質(zhì)量控制中進(jìn)行了廣泛的應(yīng)用，并取得了良好的效果，因此，文中將直接采用回歸分析方法以及數(shù)值模擬技術(shù)對(duì) Cu-45Ag合金金屬型鑄件的偏析缺陷進(jìn)行研究，獲得主要凝固參數(shù)與 Cu-45Ag合金鑄件偏析缺陷之間的定量關(guān)系，為合理優(yōu)化鑄造工藝參數(shù)提供理論支撐。

1 試驗(yàn)方法

1.1 材料及物理性能

所采用的合金為 Cu-45Ag二元過(guò)共晶合金，二元合金的相圖見圖1，可以看出，該二元合金的共晶點(diǎn)在Ag的原子數(shù)分?jǐn)?shù)為39.9%時(shí)發(fā)生共晶反應(yīng)，共晶的反應(yīng)溫度為779.3 ℃，則可以明確本實(shí)驗(yàn)所采用Ag-55%Cu（質(zhì)量分?jǐn)?shù)）二元過(guò)共晶合金的固相線溫度為779.3 ℃。過(guò)橫坐標(biāo)Cu的質(zhì)量分?jǐn)?shù)為55%的位置作垂線，與相圖中的液相線交于一點(diǎn)，過(guò)該點(diǎn)做水平線與溫度軸的交點(diǎn)為885.8 ℃，該溫度即為過(guò)共晶合金的液相線溫度。對(duì)材料的熱物性參數(shù)進(jìn)行查找和實(shí)際測(cè)定，獲得的鑄件的熱物性參數(shù)見表1。鑄型采用鑄鐵，其相應(yīng)的熱物性參數(shù)可直接通過(guò)模擬軟件的數(shù)據(jù)庫(kù)來(lái)選擇。

圖1 Ag-Cu二元合金相圖Fig.1 Phase diagram of Ag-Cu binary alloy

表1 Cu-45Ag二元合金的熱物性參數(shù)Tab.1 Thermophysical parameters of Cu-45Ag binary alloy

1.2 模型的建立

鑄件和鑄型三維造型以及網(wǎng)格剖分結(jié)果的示意圖見圖2，采用Solid Works三維造型軟件對(duì)擠壓鑄造的鑄件/鑄型系統(tǒng)部分進(jìn)行三維造型，造型方法為實(shí)體造型，所獲得的三維實(shí)體包括鑄件、涂層以及鑄型部分，然后將各個(gè)部分在軟件中進(jìn)行裝配，并將裝配好的文件和各部分分別保存為擴(kuò)展名為 step的文件。實(shí)際鑄件將在金屬型中進(jìn)行澆注，鑄錠的尺寸為厚34 mm，寬125 mm，高260 mm。澆注過(guò)程采用頂注式澆注，不再單獨(dú)設(shè)置澆冒口，如圖 2a所示。有限元網(wǎng)格的剖分結(jié)果見圖 2b，其中鑄型的網(wǎng)格長(zhǎng)度為8 mm，鑄件的網(wǎng)格長(zhǎng)度為2 mm，總網(wǎng)格數(shù)約為170萬(wàn)個(gè)。

圖2 鑄件的三維模型Fig.2 3D-model of casting

1.3 R.G.L開放模塊

ProCAST模擬軟件在其后處理模塊中提供了一個(gè)可供用戶自行設(shè)計(jì)判據(jù)的求解器，即R.G.L開放模塊，模塊中的求解器主要針對(duì)主要凝固參數(shù)進(jìn)行設(shè)計(jì)和計(jì)算，即該模塊可直接計(jì)算出冷卻速率、溫度梯度、凝固速率這3個(gè)主要凝固參數(shù)，并根據(jù)用戶需求，通過(guò)其提供的開放公式進(jìn)行計(jì)算和預(yù)測(cè)模型的開發(fā)，開放公式見式(1)。

式中：M為所建立預(yù)測(cè)模型的評(píng)價(jià)值，即因變量；U為凝固速率；G為溫度梯度；R為冷卻速率；a,c,d為常數(shù)。

顯然對(duì)于偏析缺陷而言，其主要形成于合金凝固的糊狀區(qū)內(nèi)，因此判據(jù)公式中所需要的各凝固參數(shù)的計(jì)算就應(yīng)當(dāng)確定在這個(gè)溫度范圍內(nèi)，溫度梯度和凝固速率的計(jì)算則需要在特定的溫度點(diǎn)下進(jìn)行。根據(jù)圖1給出的 Cu-Ag的二元相圖，各凝固參數(shù)的選取溫度范圍和取值分別由式(2)和式(3)給出。

式中：TL為液相線溫度；TS為固相線溫度；Tupper和Tlower則分別表示計(jì)算冷卻速率溫度范圍的上限值和下限值；TR,G為計(jì)算凝固速率和溫度梯度所給定的溫度值。文中將基于這一開放模塊，對(duì)相應(yīng)缺陷的預(yù)測(cè)模型進(jìn)行相應(yīng)的研究和建立。

1.4 多元非線性回歸分析

作為一種常用的數(shù)據(jù)分析統(tǒng)計(jì)處理模型，多元非線性回歸分析在工程中的應(yīng)用非常廣泛，是研究解決零件生產(chǎn)中制造工藝參數(shù)與產(chǎn)品質(zhì)量相關(guān)關(guān)系的一種強(qiáng)有力的數(shù)理統(tǒng)計(jì)方法。與其他方法相比，該方法不僅同樣具有簡(jiǎn)捷方便、精確度和可靠性高等特點(diǎn)，而且還能建立合適的函數(shù)關(guān)系式來(lái)表明各個(gè)工藝參數(shù)與產(chǎn)品質(zhì)量之間的相互聯(lián)系[13—17]，該方程可以突破試驗(yàn)設(shè)計(jì)對(duì)工藝參數(shù)范圍的限制，從而得到理想的優(yōu)化工藝參數(shù)，因此將通過(guò)這種工程統(tǒng)計(jì)學(xué)的方法，分析研究Cu-Ag合金金屬型鑄造過(guò)程的凝固參數(shù)與宏觀偏析缺陷之間的定量關(guān)系，最終獲得快速預(yù)測(cè)該缺陷且具有工程統(tǒng)計(jì)學(xué)意義的數(shù)學(xué)模型，從而達(dá)到指導(dǎo)實(shí)際生產(chǎn)過(guò)程、提高鑄件質(zhì)量、減少試制時(shí)間的目的。

采用專業(yè)數(shù)據(jù)處理軟件 SAS軟件對(duì)模擬計(jì)算所獲得的數(shù)據(jù)進(jìn)行多元非線性回歸分析，分析過(guò)程分為以下幾個(gè)步驟： ①模型假設(shè)，首先對(duì)回歸方程的形式進(jìn)行假設(shè)，這部分將主要參考一些宏觀偏析缺陷預(yù)測(cè)的模型； ②參數(shù)估計(jì)，在SAS軟件中進(jìn)行編程，給出假設(shè)的回歸方程，計(jì)算出相應(yīng)參數(shù)的估計(jì)值，多元非線性回歸采用最小二乘法進(jìn)行估計(jì)； ③顯著性檢驗(yàn)，軟件運(yùn)行后，將在結(jié)果中直接輸出顯著性的檢驗(yàn)結(jié)果，根據(jù)結(jié)果對(duì)所假設(shè)的模型進(jìn)行評(píng)價(jià)； ④模型修正，如果步驟①中假設(shè)的模型沒有通過(guò)檢驗(yàn)，則需要對(duì)回歸方程的形式進(jìn)行修正，并重復(fù)上述步驟，直到回歸方程通過(guò)檢驗(yàn)為止； ⑤試驗(yàn)驗(yàn)證，通過(guò)實(shí)際試驗(yàn)結(jié)果對(duì)回歸方程的適用性進(jìn)行驗(yàn)證。

2 結(jié)果與討論

2.1 主要凝固參數(shù)

由于凝固過(guò)程中的凝固參數(shù)主要有凝固速率、溫度梯度和冷卻速率，而凝固速率可以通過(guò)溫度梯度和冷卻速率來(lái)進(jìn)行表示，因此，采用溫度梯度和冷卻速率這2個(gè)主要凝固參數(shù)來(lái)進(jìn)行宏觀偏析的預(yù)測(cè)。在仿真過(guò)程中在鑄件表面進(jìn)行選點(diǎn)，選點(diǎn)位置及編號(hào)見圖3。分別在鑄件正面、左側(cè)以及上下部 4個(gè)表面位置選取特征點(diǎn)中上表面的特征點(diǎn)命名為T1～T3，下表面的的特征點(diǎn)命名為 B1～B3，正面的特征點(diǎn)命名為Z11～Z34，左側(cè)特征點(diǎn)的命名為C1～C4，所有的特征點(diǎn)均勻分布。

圖3 各參考點(diǎn)的選取Fig.3 Position of each reference point

為了獲得實(shí)際鑄件的宏觀偏析缺陷的產(chǎn)生情況，采用澆注工藝參數(shù)：澆注溫度為1100 ℃、澆注時(shí)間為30 s、鑄型溫度300 ℃、涂料厚度1.0 mm進(jìn)行實(shí)際澆注，澆注完成后，利用直讀光譜分析儀對(duì)鑄件進(jìn)行Ag含量的分析，分析的位置與圖3給出的特征點(diǎn)位置相一致。利用 ProCAST采用相同的工藝參數(shù)進(jìn)行鑄件充型和凝固過(guò)程的模擬，計(jì)算完成后對(duì)計(jì)算結(jié)果中對(duì)各個(gè)面上的特征點(diǎn)位置進(jìn)行溫度梯度和冷卻速率數(shù)值提取，所獲取的各特征點(diǎn)的凝固信息與 Ag含量見表2，所給出的溫度梯度單位均為K/cm2，冷卻速率單位均為K/s。

表2 特征區(qū)域部位的凝固信息和Ag的質(zhì)量分?jǐn)?shù)Tab.2 Solidification information and Ag content of feature areas (mass fraction)

2.2 回歸分析

通過(guò)表2列出的實(shí)際澆注與數(shù)值模擬計(jì)算結(jié)果，可以看出對(duì)于上表面的凝固信息而言，其冷卻速度與側(cè)邊相近，但溫度梯度較小，而產(chǎn)生了較大的宏觀偏析缺陷。鑄件正面上部的各個(gè)特征點(diǎn)的冷卻速度和溫度梯度均較小，也同樣產(chǎn)生了較大的宏觀偏析，因此，可以認(rèn)為冷卻速度和溫度梯度越小，越容易產(chǎn)生宏觀偏析，且溫度梯度對(duì)宏觀偏析的形成影響更大。通過(guò)對(duì)數(shù)據(jù)進(jìn)行擬合，得到的回歸方程見式(4)。

式中：A=40.256 26,B=0.017 27。該表達(dá)式與文獻(xiàn)[24]所報(bào)道的在Al-7%Si合金鑄件中進(jìn)行宏觀偏析預(yù)測(cè)的判據(jù)形式相類似。此外，從模型中也可以看出，由于溫度梯度G的指數(shù)為2，因此溫度梯度對(duì)宏觀偏析形成的影響比冷卻速度要大。

回歸模型的顯著性檢驗(yàn)結(jié)果和參數(shù)估計(jì)的顯著性檢驗(yàn)結(jié)果分別見表3和表4?？梢钥闯?，不管是F檢驗(yàn)還是t檢驗(yàn)的p值均能夠滿足回歸方程(4)：p＜0.0001＜α。可見，提出的模型很好地符合表2中所列出的實(shí)測(cè)數(shù)據(jù)。

表3 回歸模型的顯著性檢驗(yàn)Tab.3 Analysis of Variance for model

表4 回歸模型的參數(shù)估計(jì)的顯著性檢驗(yàn)Tab.4 Significance testing for parameter estimation of regression model

3 工藝改進(jìn)和試驗(yàn)驗(yàn)證

根據(jù)宏觀偏析預(yù)測(cè)模型，可以在實(shí)際澆注之前的前期模擬計(jì)算中對(duì)Ag元素的含量進(jìn)行大致的預(yù)測(cè)，從而判斷Cu-Ag合金中宏觀偏析的大致位置和程度。為了對(duì) Cu-Ag合金的澆注工藝進(jìn)行改進(jìn)，分別對(duì)影響冷卻速度和溫度梯度的主要澆注工藝參數(shù)，即澆注溫度、澆注時(shí)間（澆注速度）、鑄型溫度和涂料厚度分別進(jìn)行單因素調(diào)整模擬試驗(yàn)。計(jì)算結(jié)果表明，澆注時(shí)間（澆注速度）對(duì)主要凝固參數(shù)的影響較大，澆注溫度、鑄型溫度和涂料厚度的影響較小。此外，降低澆注速度可以有效提高冷卻速度和溫度梯度，從而有效減小宏觀偏析的傾向。

在通過(guò)確定最佳澆注工藝參數(shù)的同時(shí)，為了對(duì)模型的準(zhǔn)確性和實(shí)用性進(jìn)行驗(yàn)證，采用模擬后獲得的最佳澆注工藝方案進(jìn)行澆注，即澆注溫度為1100 ℃、澆注時(shí)間為120 s、鑄型溫度為200 ℃、涂料厚度為1.5 mm，并對(duì)實(shí)際澆注獲得的鑄件進(jìn)行檢測(cè)，結(jié)果表明，通過(guò)最優(yōu)澆注工藝參數(shù)獲得鑄件中的宏觀偏析缺陷已經(jīng)得到明顯的改善，測(cè)量獲得的數(shù)值與預(yù)測(cè)的數(shù)值基本接近，均位于預(yù)測(cè)優(yōu)化的數(shù)值范圍之內(nèi)。

4 結(jié)論

通過(guò)采用多元非線性回歸方法和數(shù)值模擬對(duì)Cu-45Ag合金金屬型鑄造過(guò)程產(chǎn)生的宏觀偏析缺陷的預(yù)測(cè)模型進(jìn)行研究，并根據(jù)獲得的回歸模型進(jìn)行了相應(yīng)的仿真研究和工藝改進(jìn)，通過(guò)對(duì)試驗(yàn)結(jié)果的分析研究，得到了以下結(jié)論。

1)通過(guò)回歸分析方法，獲得的Cu-45Ag合金金屬型鑄造宏觀偏析的數(shù)值模型為：

通過(guò)對(duì)回歸模型和參數(shù)估計(jì)的顯著性檢驗(yàn)結(jié)果進(jìn)行分析，表明該方程能夠很好地表達(dá)主要凝固參數(shù)與宏觀偏析之間的對(duì)應(yīng)關(guān)系。

2)針對(duì)優(yōu)化工藝參數(shù)進(jìn)行計(jì)算，其結(jié)果表明，澆注時(shí)間（澆注速度）對(duì)主要凝固參數(shù)的影響較大，澆注溫度、鑄型溫度和涂料厚度的影響較小。降低澆注速度可以有效提高冷卻速度和溫度梯度，從而有效減小宏觀偏析的傾向。

3)獲得的最佳澆注工藝方案為澆注溫度1100 ℃、澆注時(shí)間為120 s、鑄型溫度200 ℃、涂料厚度 1.5 mm。采用該方案進(jìn)行實(shí)際澆注，對(duì)所獲得鑄件的檢測(cè)結(jié)果表明，鑄件中的宏觀偏析缺陷得到了明顯的改善，驗(yàn)證試驗(yàn)的結(jié)果還同時(shí)表明，提出的宏觀偏析預(yù)測(cè)模型能夠有效對(duì) Cu-45Ag合金金屬型鑄件進(jìn)行偏析預(yù)測(cè)，為工藝設(shè)計(jì)人員進(jìn)行合理的工藝優(yōu)化提供幫助。