韓茂盛，陶 歡

（中國船舶重工集團公司 第七二五研究所，河南 洛陽 471000）

B10銅合金高溫流變行為及BP神經(jīng)網(wǎng)絡本構模型

韓茂盛，陶 歡

（中國船舶重工集團公司 第七二五研究所，河南 洛陽 471000）

利用Gleeble-1500D熱模擬機對B10銅合金進行熱壓縮實驗，研究了該合金在高溫塑性變形過程中的流變應力行為。實驗溫度為800～950℃，應變速率為0.1～15s-1。研究結果表明，B10銅合金的流變應力隨著變形溫度的增加而減小，隨著應變速率的增大而增大?；贐P神經(jīng)網(wǎng)絡建立了該合金的本構關系模型，預測值與實驗值對比表明BP神經(jīng)網(wǎng)絡具有很高的預測精度，所建立的本構模型平均相對誤差在1%以內。該模型能夠客觀真實地描述B10銅合金的高溫塑性變形行為，為該合金熱變形分析提供基礎。

B10銅合金；流變行為；本構關系；BP神經(jīng)網(wǎng)絡

B10合金是以鎳、鐵為主要合金元素的銅合金，具有優(yōu)良的耐海水腐蝕性能和抗海生物污損性能，因此廣泛應用于艦船、海上石油工程的海水管路系統(tǒng)[1,2]。目前，B10合金主要是以管材及管系附件（如凸緣、法蘭等）的形式應用，而其成形方式需要經(jīng)過高溫塑性加工成形[3]。金屬熱變形流變應力是材料在高溫下的基本性能之一，不僅受變形溫度、變形程度、應變速率和合金化學成分的影響，也是變形體內部顯微組織演變的綜合反映，因此，研究材料的高溫流變行為是制定熱加工工藝的關鍵。

材料的本構關系是材料的流動應力與變形熱力參數(shù)之間的函數(shù)關系，是聯(lián)系材料塑性變形行為與各種成形參數(shù)的橋梁，在以有限元為代表的現(xiàn)代塑性加工力學中，金屬的流變應力值是提高理論計算精度的關鍵[4,5]。目前，建立材料本構關系時多根據(jù)實驗數(shù)據(jù)，借助于回歸模型來獲得一些經(jīng)驗公式以滿足工程技術要求。但是材料的變形過程是一個非常復雜的非線性問題，在實驗過程中容易受到各種隨機因素的干擾，變形條件與流變應力的關系很難用簡單的數(shù)學模型準確描述。

人工神經(jīng)網(wǎng)絡具有簡單的非線性處理單元的復合映射即可獲得復雜的非線性處理能力，不需要給出數(shù)學模型，通過對給定的樣本數(shù)據(jù)進行學習，從大量的數(shù)據(jù)中提取規(guī)則，以一組權重形式形成一種網(wǎng)絡的穩(wěn)定狀態(tài)，然后通過聯(lián)想記憶和推廣能力來獲取所需數(shù)據(jù)。BP算法是普遍使用的一種神經(jīng)網(wǎng)絡學習算法，具有一定的聯(lián)系容錯能力。因此，利用BP算法來建立材料本構關系獲得了較為廣泛的應用[6-8]。

本文利用Gleeble-1500D熱模擬機對B10合金進行高溫等溫壓縮實驗，研究其在熱變形條件下的流變應力行為，并以不同變形溫度、不同應變速率和應變下的等溫恒應變速率壓縮實驗數(shù)據(jù)作為學習樣本，利用BP神經(jīng)網(wǎng)絡建立B10合金的本構關系模型，并對其進行誤差分析，為合理制定B10合金熱加工工藝提供指導依據(jù)。

1 實驗

實驗用原材料為B10合金鑄錠，其化學成分見表1。將鑄錠在950℃進行均勻化處理后，加工成?10mm×15mm的圓柱試樣。實驗設備采用Gleeble-1500D熱模擬實驗機，合金的變形溫度T分別為800℃、850℃、800℃、950℃，應變速率ε˙分別為0.1、1、10、15s-1，變形量均為60%。以10℃/s的加熱速率將試樣加熱到變形溫度，保溫3min，變形結束后水淬。

表1 實驗合金化學成分（Wt%）

2 B10銅合金高溫流變行為

圖1為B10合金熱壓縮變形時的真應力應變曲線。從圖中可以看出，當變形溫度和應變速率一定時，真應力總體趨勢隨著應變的增加先迅速增大，隨后增大速率減緩，最后基本趨于穩(wěn)定，這是因為熱壓縮開始階段，由于位錯的交滑移和攀爬引起的軟化不足以補償由位錯密度增加帶來的硬化，因此真應力急劇增大；隨著變形量的增加，晶粒內參與滑移的可動位錯數(shù)量增加，軟化作用增強，與加工硬化趨于平衡，真應力亦趨于平穩(wěn)。由圖中還能發(fā)現(xiàn)，峰值應力及穩(wěn)態(tài)應力隨著變形溫度的降低和應變速率的提高呈增加趨勢，且溫度越低，增加的幅度越明顯，表明該合金為正應變敏感材料。當變形溫度為800℃～850℃時，隨著真應變的增加，真應力一直在增加，表明變形過程中動態(tài)硬化占主導地位，塑性變形較為困難，因而難以進行變形加工；當變形溫度達到900℃以上時，真應力達到峰值后不隨著真應變的增加而增大，甚至出現(xiàn)動態(tài)軟化，特別是應變速率較低、變形溫度較高時，應力下降的趨勢尤為明顯，這一方面是熱壓縮變形時發(fā)生了動態(tài)回復和動態(tài)再結晶的作用，另一方面由于熱壓縮變形中的絕熱溫升現(xiàn)象導致實際溫度隨著壓縮過程的進行有所升高，一定程度上導致了變形抗力的下降。

3 B10銅合金B(yǎng)P神經(jīng)網(wǎng)絡本構模型

3.1 模型建立

BP神經(jīng)網(wǎng)絡即采用BP算法的多層感知器神經(jīng)網(wǎng)絡模型，由輸入層、隱含層和輸出層組成，其訓練過程由信息正向傳播和誤差反向傳播兩部分組成，正向傳播時，輸入信息由輸入層經(jīng)隱含層單元處理后傳向輸出層，每一層神經(jīng)元的狀態(tài)只影響下一層神經(jīng)元狀態(tài)。若在輸出層得不到希望輸出，則轉入誤差反向傳播，將誤差信號沿原神經(jīng)元連接通路返回。返回過程中，逐一修改各層神經(jīng)元連接的權值，通過上述過程的不斷迭代，最終使得信號誤差達到允許范圍內，網(wǎng)絡結束學習訓練過程[9-10]。本文采用三層BP神經(jīng)網(wǎng)絡模型，輸入層包含變形溫度T、真應變ε和應變速率ε˙三個節(jié)點，輸出層只有一個輸出量，為真應力σ，如圖2所示。

3.2 數(shù)據(jù)預處理

圖1 B10合金在不同應變速率下熱壓縮變形的真應力-應變曲線

從真應力-應變曲線中選取624組數(shù)據(jù)點，其中488組作為訓練樣本，176組作為檢驗樣本。由于變形溫度、變形速率和真應變的數(shù)值相差較大，如果同時用于網(wǎng)絡的輸入節(jié)點，勢必會湮沒小數(shù)據(jù)對徑向函數(shù)的作用，造成隱含層與輸入層之間權值調節(jié)上的困難，從而影響網(wǎng)絡的收斂速度和精度，因此根據(jù)BP神經(jīng)網(wǎng)絡輸入輸出規(guī)范要求，在建立網(wǎng)絡之前必須對所有數(shù)據(jù)進行歸一化處理，使各特征變量在0～1范圍內取值，避免網(wǎng)絡模型迭代計算時數(shù)值過于集中在某些神經(jīng)元和權值上而降低模型的解算精度，計算公式如下：

圖2 BP神經(jīng)網(wǎng)絡模型

式中，X為實驗所獲得的原始數(shù)據(jù)，Xmin和Xmax為原始相應數(shù)據(jù)的最小值和最大值，X′為參數(shù)X歸一化處理后的數(shù)值。

公式（1）用來歸一化溫度和流變應力，由于應變已經(jīng)在0～1之間，不需要歸一化處理。但是，由于應變速率變化很大，而且歸一化后數(shù)值很小，因此，采用公式（2）進行歸一化。

3.3 模型訓練

采用Matlab將歸一化后的訓練樣本對BP神經(jīng)網(wǎng)絡進行訓練，轉移函數(shù)為“tan sigmoid”和“pure linear”，訓練函數(shù)為"Trainlm"，訓練目標誤差為10-4。在網(wǎng)絡模型訓練時，隱含層神經(jīng)元數(shù)目的選擇是一個至關重要的問題，數(shù)目太少會導致訓練出來的網(wǎng)絡容錯性差，不能識別參與未參與訓練的樣本；數(shù)目過多會導致學習時間過長，誤差不一定最佳，因此隱含層單元數(shù)存在一個最佳值[11]。為了確定隱藏層的神經(jīng)元的數(shù)量，從兩個神經(jīng)元開始，采用試錯法進行一一訓練，圖3為不同隱藏神經(jīng)元的網(wǎng)絡性能。從圖中可以看出，隨著隱含層神經(jīng)元個數(shù)的增加，網(wǎng)絡性能開始呈現(xiàn)顯著提高的趨勢，當隱藏層存在16個神經(jīng)元時，具有最小的平均方根誤差，因此，一個具有16個神經(jīng)元的隱藏層用于預測B10銅合金的流變應力。BP網(wǎng)絡模型經(jīng)過1200個循環(huán)達到穩(wěn)定狀態(tài)。

圖3 不同隱含層神經(jīng)元個數(shù)的網(wǎng)絡性能

3.4 結果分析

由于BP神經(jīng)網(wǎng)絡預測輸出的真應力值在0～1之間，因此需要反歸一化處理后才能與實驗值進行比較。將公式（1）進行移項處理，并替換其中的相應項，便可得到反歸一化運算公式：

式中，σ0為預測真應力值，σmax、σmin為檢驗樣本真應力的最大值和最小值，σn為網(wǎng)絡輸出的歸一化處理后的真應力值。

圖4為采用BP神經(jīng)網(wǎng)絡模型預測的B10銅合金流變應力與實驗值的對比圖，從圖中可以看出，采用BP神經(jīng)網(wǎng)絡模型預測的流變應力值與實驗值吻合良好，且具有較高的精度。

圖5為BP神經(jīng)網(wǎng)絡模型預測的流變應力值與實驗值的關聯(lián)性曲線，從圖中可以看出，BP神經(jīng)網(wǎng)絡預測的流變應力與實驗數(shù)據(jù)的相關性較好，擬合相關系數(shù)達到0.9999。

對于模型精度的評價，可用平均相對誤差表示：

式中，Ei為流變應力的實驗值，Pi為通過BP模型計算的預測數(shù)據(jù)，N為實驗分析數(shù)據(jù)點的個數(shù)。經(jīng)計算，得到平均相對誤差為0.56%，進一步證明BP神經(jīng)網(wǎng)絡模型具有較高的精度，可滿足工程計算要求。

4 結論

圖4 B10銅合金流變應力值與實驗值的對比

圖5 BP神經(jīng)網(wǎng)絡模型預測流變應力值與實驗值的關聯(lián)性曲線

（1）通過熱壓縮實驗得到了B10銅合金在應變速率為0.1～15s-1和變形溫度為800～900℃的真應力-應變曲線，發(fā)現(xiàn)B10銅合金的流變應力隨著變形溫度的增加而降低，隨著應變速率的增大而增大。

（2）基于實驗獲得的樣本數(shù)據(jù)，采用BP神經(jīng)網(wǎng)絡模型對樣本數(shù)據(jù)進行訓練和預測并進行驗證，結果顯示計算值和預測值的擬合相關系數(shù)達0.9999，得到的平均相對誤差在1%以內，表明BP神經(jīng)網(wǎng)絡模型具有較高的精度，能較好地滿足工程計算要求。

[1]林樂耘，徐 杰，趙月紅.國產(chǎn)B10銅合金海水腐蝕行為研究[J].中國腐蝕與防護學報，2000，20（6）：361-364.

[2]鄭俊濤，張三平，周學杰，等.白銅B10和B30沖刷腐蝕對比研究[J].裝備環(huán)境工程，2010，27(1)：40-42.

[3]王祝堂，田容璋.銅合金及其加工手冊[M].長沙：中南大學出版社，2001.

[4]Sha W,Edwrads K L.The use of artificial neural networks in materials science based research[J].Materials and Design,2007, 28:1747-1752.

[5]Mordike B L,Ebert T.Magnesium:properties-applicationspotential[J].Materials Science and Technology,1997,13:210.

[6]田希杰，高 軍.基于人工神經(jīng)網(wǎng)絡的冷擠壓精度預測[J].鍛壓裝備與制造技術，2007，6.

[7]李慧麗，吳朋越，謝水生.利用BP神經(jīng)網(wǎng)絡預測銅合金擠壓力[J].鍛壓裝備與制造技術，2005，2（22）.

[8]Bariani P F,Bruschi S,Negro T.Prediction of nickel-base superalloys'rheologicalbehavior under hot forging conditions using artificial neural networks[J].Journal of Materials Processing Technology,2004,152:395-400.

[9]張興全，彭穎紅.Ti-17合金本構關系的人工神經(jīng)網(wǎng)絡模型[J].中國有色金屬學報，1999，9（3）：590-595.

[10]劉學峰，馬勝軍，劉錦平.Cu-12%Al合金高溫壓縮變形過程中本構關系的 BP神經(jīng)網(wǎng)絡模型[J].材料工程，2009，1：10-13.

[11]樓順天，施 陽.基于MATLAB的系統(tǒng)分析與設計-神經(jīng)網(wǎng)絡[M].西安：西安電子科技大學出版社，1998.

High temperature flow stress behavior of B10 copper alloy and BP neural network constitutive model

HAN Maosheng,TAO Huan
（Luoyang Ship Material Research Institute,Luoyang 471000,Henan China）

The hot compression test has been conducted to B10 copper alloy by use of Gleeble-1500D simulator under temperature range from 800℃ to 950℃ andstrain rate range from 0.1s-1 to 15s-1.The flow stress behavior of this alloyduring plastic deformation at elevated temperature has been studied.The experimental results show that the flow stress of B10 copper alloy decreases with the increase of deformation temperature but enhances with the increase of the stain rate.The constitutive relationship model has been built to this alloy on the basis of BP neural network.By comparison between the predicted value with the experiment value,it is shown that BP neural network has high prediction accuracy.The average relative error of the constitutive model is less than 1%.This model can objectively and truly describe the plastic deformation behavior of the B10 copper alloy at high temperature.It provides basis for hot deformation analysisof this alloy.

B10 copper alloy;Flow stress behavior;Constitutive relationship;BP neural network

TG146.1

A

10.16316/j.issn.1672-0121.2016.06.028

1672-0121（2016）06-0112-04

2016-04-27；

2016-06-15

韓茂盛（1987-），男，碩士，工程師，從事材料制備、加工方面研究。E-mail：cfdeng2000@163.com