鎢銅合金熱變形行為及本構(gòu)方程

摘 要：【目的】鎢銅合金是航空航天和軍工領(lǐng)域用關(guān)鍵材料，但鎢銅合金流動(dòng)應(yīng)力應(yīng)變的研究采用的本構(gòu)方程的準(zhǔn)確性較差，需要對(duì)其加以研究?！痉椒ā坎捎谜婵绽煸囼?yàn)機(jī)對(duì)鎢銅合金進(jìn)行了不同溫度、不同應(yīng)變速率下的高溫拉伸實(shí)驗(yàn)，研究其熱變形行為，并建立了鎢銅合金的本構(gòu)方程?！窘Y(jié)果】結(jié)果表明，鎢銅合金存在著明顯的溫度軟化效應(yīng)和加工硬化效應(yīng)，其流變應(yīng)力隨溫度升高和應(yīng)變速率降低而降低。利用相關(guān)系數(shù)和絕對(duì)平均誤差可衡量本構(gòu)方程的準(zhǔn)確性?！窘Y(jié)論】應(yīng)變補(bǔ)償程可大幅提高鎢銅合金本構(gòu)的準(zhǔn)確度。

關(guān)鍵詞：鎢銅合金；熱變形行為；本構(gòu)方程

中圖分類號(hào)：TG146.1" " "文獻(xiàn)標(biāo)志碼：A" " 文章編號(hào)：1003-5168（2024）24-0069-10

DOI：10.19968/j.cnki.hnkj.1003-5168.2024.24.014

Thermal Deformation Behavior and Constitutive Equation of Tungsten-Copper Alloy

Abstract： [Purposes] Tungsten-copper alloy is a key material used in the aerospace and military fields. However， the accuracy of the constitutive equations adopted in the research on the flow stress and strain of tungsten-copper alloy is rather poor， which needs to be further studied. [Methods] High-temperature tensile experiments on tungsten-copper alloys were carried out under different temperatures and strain rates using a vacuum tensile testing machine to study the thermal deformation behavior， and the constitutive equation of the tungsten-copper alloy was established. [Findings] The results show that the tungsten-copper alloy has obvious temperature softening effect and work hardening effect， and its flow stress decreases with the increase of temperature and the decrease of strain rate. The accuracy of the constitutive equation can be measured by using the correlation coefficient and the absolute average error. [Conclusions] Strain compensation can significantly improve the accuracy of the constitutive equation of tungsten-copper alloy.

Keywords： tungsten-copper alloy; thermal deformation behavior; constitutive equation

0 引言

鎢銅合金是由互不相容的金屬鎢和銅組成的，其既具有鎢的高強(qiáng)度、高導(dǎo)熱性能，又具有銅的高塑性、高導(dǎo)電性等優(yōu)異性能，被廣泛應(yīng)用于航空航天、軍工、電工電子等領(lǐng)域［1］。鎢銅合金多采用粉末冶金的方法制備，存在致密度低、內(nèi)部孔隙等缺陷，一般采用變形強(qiáng)化的方法消除這些缺陷［2］。采用試驗(yàn)的方法對(duì)鎢銅合金在不同變形條件下的變形行為進(jìn)行研究，既耗材又耗費(fèi)時(shí)間，因此，科研人員多采用ABAQUS、DEFORM等商業(yè)有限元軟件對(duì)塑性成形過(guò)程進(jìn)行仿真分析［3］。對(duì)于有限元仿真來(lái)說(shuō)，建立合適的材料模型無(wú)疑是重要的一步［4］。一般來(lái)說(shuō)，可以通過(guò)構(gòu)建材料的本構(gòu)模型來(lái)描述材料在不同溫度和應(yīng)變速率下的流動(dòng)應(yīng)力應(yīng)變狀態(tài)［5］。用于預(yù)測(cè)不同溫度和應(yīng)變速率下的流變應(yīng)力的本構(gòu)方程主要有基于現(xiàn)象學(xué)的本構(gòu)模型、物理學(xué)本構(gòu)模型等［6］?；诂F(xiàn)象學(xué)的本構(gòu)模型是基于材料的宏觀變形建立的，主要基于經(jīng)驗(yàn)觀察和數(shù)學(xué)函數(shù)模型，不需要考慮熱變形中的物理意義，在很多材料上得到了應(yīng)用。常見(jiàn)的現(xiàn)象學(xué)的本構(gòu)模型有：Johnson-Cook模型（JC模型）［7］、修正后的Johnson-Cook模型［8］、和Arrhenius模型［9］。物理學(xué)本構(gòu)模型考慮了材料熱變形過(guò)程中的實(shí)際物理意義，能夠有效模擬宏觀變形與微觀組織變化之間的聯(lián)系。常見(jiàn)的物理學(xué)本構(gòu)模型有Zerilli-Armstrong模型［10］（ZA模型）、PB模型［11］和MTS模型［12］等。這些本構(gòu)模型都可以通過(guò)試驗(yàn)和數(shù)據(jù)分析求得相關(guān)參數(shù)，進(jìn)而得到完整的本構(gòu)模型。

目前，鎢銅合金高溫下的流動(dòng)應(yīng)力應(yīng)變研究較多。趙艷平等［13］采用熱壓縮試驗(yàn)探究了溫度和應(yīng)變速率對(duì)鎢銅合金變形的影響，結(jié)果表明：抗拉強(qiáng)度隨應(yīng)變速率的減小或溫度的增加而減小。劉勇等［14］進(jìn)一步研究了不同鎢含量的鎢銅合金熱變形行為，發(fā)現(xiàn)在所有鎢含量下鎢銅合金均對(duì)溫度和應(yīng)變速率敏感，并結(jié)合應(yīng)力應(yīng)變曲線，建立了Arrhennius本構(gòu)方程。王凌浩等［15］對(duì)鎢銅合金進(jìn)行熱模擬壓縮試驗(yàn)，指出W-20Cu、W-35Cu的熱加工溫度在700～900 ℃比較合適。并建立了W-35Cu的Johnson-Cook本構(gòu)方程及W-20Cu的Arrhennius本構(gòu)方程，并采用了四次多項(xiàng)式擬合了Arrhennius本構(gòu)方程中參數(shù)與應(yīng)變的關(guān)系，對(duì)Arrhennius本構(gòu)方程進(jìn)行修正。彭付申等［16］對(duì)溶滲法制備的W-20Cu進(jìn)行熱壓縮試驗(yàn)，建立了Arrhennius本構(gòu)方程，并采用了五次多項(xiàng)式擬合了材料參數(shù)與應(yīng)變的關(guān)系，對(duì)Arrhennius本構(gòu)進(jìn)行修正。Lin［17］等通過(guò)高溫下的應(yīng)變-應(yīng)力關(guān)系，確定了W-10%Cu合金的熱變形方程（Arrhenius本構(gòu)方程），其中應(yīng)變速率與靜態(tài)應(yīng)力、變形溫度和激活能有關(guān)。劉輝明等［18］通過(guò)對(duì)W-40Cu和W-50Cu的動(dòng)態(tài)壓縮試驗(yàn)，探究了應(yīng)變速率和溫度對(duì)鎢銅合金屈服強(qiáng)度的影響，并建立了JC本構(gòu)方程。馬竇琴［19］通過(guò)熱模擬機(jī)測(cè)定了W-25wt %Cu燒結(jié)坯的應(yīng)力應(yīng)變曲線，建立了基于現(xiàn)象學(xué)的Arrhenius本構(gòu)方程，并對(duì)鎢銅合金包套擠壓過(guò)程進(jìn)行了仿真模擬。

上述鎢銅合金流動(dòng)應(yīng)力應(yīng)變的研究所采用的W含量一般在80 wt%以下且一般構(gòu)建Arrhennius本構(gòu)方程，其準(zhǔn)確性有待定量計(jì)算。基于此，本研究對(duì)鎢銅合金開(kāi)展了不同溫度和不同速率的高溫拉伸試驗(yàn)，以此建立鎢銅合金的不同本構(gòu)方程，并對(duì)其準(zhǔn)確性進(jìn)行定量計(jì)算。

1 試驗(yàn)材料與方法

本研究采用浙江五洲新春提供的鎢銅合金。拉伸試樣的尺寸如圖1所示。采用高溫電子真空拉伸試驗(yàn)機(jī)（UZDL50，中機(jī)設(shè)備裝備股份有限公司，中國(guó)，長(zhǎng)春）進(jìn)行拉伸試驗(yàn)。拉伸試驗(yàn)按照《GB/T 228.2—2015，金屬材料" 拉伸試驗(yàn)" 第2部分：高溫試驗(yàn)方法》［20］。對(duì)鎢銅合金進(jìn)行不同溫度及不同應(yīng)變速率的高溫拉伸試驗(yàn)時(shí)，先將制取的試樣利用砂紙打磨至2 000目，確保表面光滑，無(wú)缺陷等引起的應(yīng)力集中。在進(jìn)行高溫拉伸試驗(yàn)時(shí)，在真空中進(jìn)行，升溫速率選擇8 ℃/min，升溫至指定溫度后保溫20 min，隨后進(jìn)行高溫拉伸試驗(yàn)。

2 試驗(yàn)結(jié)果

2.1 高溫拉伸行為

鎢銅合金高溫拉伸測(cè)試在高溫電子真空拉伸試驗(yàn)機(jī)UZDL50上進(jìn)行，高溫拉伸試驗(yàn)溫度分別為400 ℃、500 ℃、600 ℃、800 ℃，應(yīng)變速率為0.006 7 s-1、0.067 s-1、0.67 s-1，共10組試驗(yàn)。

鎢銅合金在不同溫度及應(yīng)變速率下的高溫拉伸應(yīng)力應(yīng)變曲線如圖2所示。由圖2（a）可以看出，在應(yīng)變速率保持不變時(shí)，峰值應(yīng)力隨試驗(yàn)溫度的降低而增大。在試驗(yàn)溫度保持不變時(shí)，峰值應(yīng)力隨應(yīng)變速率的減小而減小。這是由于在相同應(yīng)變速率下，試驗(yàn)溫度越高，軟化速率越快，應(yīng)力值越小。此外，在試驗(yàn)溫度不變時(shí)，應(yīng)變速率引起的加工硬化隨著應(yīng)變速率的提高不斷增大，從而引起峰值應(yīng)力的提高。這也可能與變形過(guò)程中銅的動(dòng)態(tài)回復(fù)和動(dòng)態(tài)再結(jié)晶存在一定的關(guān)聯(lián)。相關(guān)研究顯示銅在熱變形過(guò)程中的動(dòng)態(tài)再結(jié)晶現(xiàn)象對(duì)銅的軟化作用不可忽視［21］。而試驗(yàn)溫度和應(yīng)變速率的提高會(huì)導(dǎo)致銅相的動(dòng)態(tài)再結(jié)晶增多，從而導(dǎo)致鎢銅合金在熱變形過(guò)程中軟化速率增大。

鎢銅合金的屈服強(qiáng)度、抗拉強(qiáng)度在高溫條件下對(duì)溫度較為敏感，在整體上是隨著溫度升高而降低，并且在400～600 ℃的溫度區(qū)間內(nèi)下降最快，表明材料有溫度軟化效應(yīng)。在400 ℃時(shí)，較高應(yīng)變速率下鎢銅合金出現(xiàn)了明顯的屈服階段，表明材料擁有很高的應(yīng)變速率強(qiáng)化效應(yīng)，如圖2（b）和圖2（c）所示。隨著鎢銅合金的應(yīng)變?cè)黾?，出現(xiàn)加工硬化，應(yīng)力先是快速增大，這個(gè)階段鎢銅合金積累了大量的位錯(cuò)，從而導(dǎo)致位錯(cuò)運(yùn)動(dòng)受到較大的阻力，鎢銅合金此時(shí)應(yīng)變的增加需要更大的應(yīng)力。而當(dāng)應(yīng)變到達(dá)一定的程度時(shí)，加工硬化導(dǎo)致鎢銅合金的能量不斷增加而進(jìn)入一種不穩(wěn)定的狀態(tài)，為位錯(cuò)的運(yùn)動(dòng)提供了一定的驅(qū)動(dòng)力，鎢銅合金產(chǎn)生軟化變形所需的應(yīng)力也趨于穩(wěn)定，此時(shí)應(yīng)力的增加較為緩慢，直至斷裂［19］。

2.2 Johnson-Cook本構(gòu)

Johnson-Cook本構(gòu)模型參數(shù)相對(duì)較少且表達(dá)方式比較簡(jiǎn)單，是目前應(yīng)用最廣泛的高溫本構(gòu)模型之一。Johnson-Cook本構(gòu)模型在1983年被首次提出，主要用于材料在高溫及不同應(yīng)變速率下的流動(dòng)行為，其考慮了應(yīng)變硬化、應(yīng)變速率強(qiáng)化和熱軟化的影響，在許多材料及仿真軟件中都得到了實(shí)現(xiàn)。Johnson-Cook本構(gòu)的表達(dá)式為式（1）。

下文中若無(wú)特殊說(shuō)明，符號(hào)及計(jì)量單位均與式（1）中相同。

本部分以溫度為400 ℃為參考溫度計(jì)算Johnson-Cook本構(gòu)，其參數(shù)求解過(guò)程如下：

①求解A、B、n。

在溫度為400 ℃，應(yīng)變速率為0.006 7 s-1時(shí)，此時(shí)式（1）簡(jiǎn)化為式（2）。

[σ=A+Bεn] （2）

即忽略應(yīng)變速率強(qiáng)化和熱軟化效應(yīng)。整理式（2）并將兩邊取自然對(duì)數(shù)可得式（3）。

[lnσ?A=lnB+nlnε] （3）

式中：A為此溫度及應(yīng)變速率下的鎢銅合金的屈服強(qiáng)度，由圖2可知A為250.63 MPa。將應(yīng)變速率為0.006 7 s-1、溫度為400 ℃的試驗(yàn)數(shù)據(jù)及A的值代入式（3），并對(duì)其進(jìn)行擬合，擬合曲線如圖3所示。n的值即是擬合曲線的斜率，B的值可以由擬合曲線的截距獲得，擬合曲線的截距即為lnB。求得A為250.63，n為0.907，B為1 152.858。

②求解C。

參考溫度為400 ℃時(shí)，式（1）可以整理變形為式（4）。

在3種不同應(yīng)變速率（0.006 7 s-1、0.067 s-1和0.67 s-1）下，自0.03至斷裂時(shí)的應(yīng)變，間隔0.01取應(yīng)變值和應(yīng)力值，將這些數(shù)據(jù)代入式（4），并進(jìn)行線性擬合，得到直線的斜率就是所求的C。求得C為0.075 75。

③求解m。

當(dāng)應(yīng)變速率為0.006 7 s-1，式（1）變形為式（5）。

[σ=（A+Bεn）（1?T?m）] （5）

對(duì)上式進(jìn)行整理，并兩邊取對(duì)數(shù)可得式（6）。

將在此應(yīng)變速率下的500 ℃、600 ℃、800 ℃的試驗(yàn)數(shù)據(jù)按照求解C的取值方法，取得應(yīng)力值和應(yīng)變值，并將其帶入式（6），進(jìn)行線性并置擬合，直線的斜率即為所求的m值，求得m值為0.798。

最終求得的鎢銅合金的Johnson-Cook本構(gòu)方程的參數(shù)值見(jiàn)表1。

根據(jù)Johnson-Cook模型，建立應(yīng)力、應(yīng)變、變形速率與變形溫度之間的關(guān)系為式（7）。

Johnson-Cook模型流動(dòng)應(yīng)力試驗(yàn)值與仿真值的對(duì)比如圖4所示。由圖4可知，除了在參考變形條件下或接近參考變形條件下的仿真值可以接受外，其余仿真值與試驗(yàn)值偏差較大.隨著應(yīng)變速率或溫度的增加，這種偏差會(huì)變大。造成這種現(xiàn)象的原因是應(yīng)變速率和溫度對(duì)應(yīng)力的影響可能是相互關(guān)聯(lián)的。隨著溫度或應(yīng)變速率的升高，這種相關(guān)性將變得更加顯著。Johnson-Cook模型假設(shè)應(yīng)變速率和溫度對(duì)應(yīng)力的影響是沒(méi)有關(guān)聯(lián)的，因此Johnson-Cook本構(gòu)模型在高溫或高應(yīng)變速率下并不適用。

Johnson-Cook模型雖然形式簡(jiǎn)單，物理解釋明顯，但很難對(duì)流動(dòng)應(yīng)力進(jìn)行精確預(yù)測(cè)。這是因?yàn)镴ohnson-Cook模型假設(shè)應(yīng)變硬化、應(yīng)變速率硬化和熱軟化是3個(gè)互不相關(guān)的效應(yīng)，忽略了變形溫度、應(yīng)變速率和應(yīng)變對(duì)材料流動(dòng)行為的耦合效應(yīng)，而變形溫度、應(yīng)變速率和應(yīng)變可能同時(shí)影響動(dòng)態(tài)再結(jié)晶過(guò)程［22］，這是本應(yīng)考慮的。因此，相關(guān)學(xué)者嘗試對(duì)Johnson-Cook模型進(jìn)行改進(jìn)，改進(jìn)的Johnson-Cook模型一般引入新的指標(biāo)變量，可以更加準(zhǔn)確地對(duì)流動(dòng)應(yīng)力進(jìn)行預(yù)測(cè)。王偉等［23］在熱壓燒結(jié)Zr-6Al-0.1B合金Johnson-Cook本構(gòu)模型的建立中引入了變形溫度敏感系數(shù)D，此時(shí)的Johnson-Cook本構(gòu)效率如式（8）。

此Johnson-Cook本構(gòu)效率中其他參數(shù)的求法不變，在求解熱軟化指數(shù)m時(shí)，可將式（8）轉(zhuǎn)化為式（9）。

將求解m的試驗(yàn)數(shù)據(jù)帶入式9，并進(jìn)行線性并置擬合，可得D=0.60。此時(shí)的JC本構(gòu)方程為式（10）所示。

加入變形溫度敏感系數(shù)的Johnson-Cook模型流動(dòng)應(yīng)力試驗(yàn)值與仿真值的對(duì)比如圖5所示。此時(shí)的應(yīng)變速率為0.006 7 s-1，由圖5可以看出，相較于圖4（a），在不同溫度下的仿真值更接近試驗(yàn)值，在高溫下誤差也更小。

2.3 基于應(yīng)變補(bǔ)償?shù)腁rrhenius本構(gòu)

Arrhenius方程是描述金屬熱變形的最通用的本構(gòu)方程。在Arrhenius模型中，固定應(yīng)變或峰值應(yīng)力下應(yīng)變率與應(yīng)力的關(guān)系可由下式確定，其通用的表達(dá)形式為式（11）。

當(dāng)流動(dòng)應(yīng)力較低時(shí)（aslt;0.8），式（11）可以簡(jiǎn)化為式（12）。

當(dāng)流動(dòng)應(yīng)力較高時(shí)（asgt;1.2），式（11）可以簡(jiǎn)化為式（13）。

以上式中：a、A、n、A1、A2、n1、n2為材料常數(shù)；Q 為材料的熱變形激活能/J·mol-1；R為氣體常數(shù)；取值為8.314 J·mol-1·K-1；

Zener等［24］指出應(yīng)變速率與試驗(yàn)溫度對(duì)材料的流動(dòng)應(yīng)力的耦合影響可以由溫度補(bǔ)償變形速率因子表示，即Z參數(shù)計(jì)算公式為式（14）。

對(duì)式（11）至式（13）分別兩端取對(duì)數(shù)，可以得到式（15）至式（17）。

在確定的溫度下，由式（15）至式（17）得出n1、n2和n的值可由式（18）至式（20）得到。

選取參考應(yīng)變?yōu)閑=0.11，將應(yīng)變?yōu)?.11時(shí)的數(shù)據(jù)代入式（18）至式（20），然后對(duì)[ln?ln.]、[?ln.]、[ln[sinh（）]?ln]分別繪制散點(diǎn)圖并進(jìn)行線性擬合，結(jié)果如圖6（a）至圖6（c）所示。由式（18）至式（20）可知，材料常數(shù)n1、n2及n分別為圖6（a）、6（b）和6（c）中擬合曲線斜率的倒數(shù)。求得的n1、n2及n的值列于表2中，由此得出n1、n2及n分別為：21.699 82、0.055 90、16.397 92，由[α=]n2/n1可以得出a的值為0.002 576。

在確定的應(yīng)變速率下，式（15）可以轉(zhuǎn)化為式（21）。

熱變形激活能Q是能夠反映材料在某些確定的變形條件下塑性變形難易程度的參數(shù)，對(duì)式（21）求偏導(dǎo)可得式（22）。

將求得的n與a的數(shù)值代入上式（22），對(duì)[ln[sinh （）]?1 000/T]進(jìn)行擬合，擬合結(jié)果如圖6（d）所示?？汕蟮肣的值，Q=155 877.4 J·mol-1。

將所得Q的值代入式（14）并對(duì)兩邊取對(duì)數(shù)可得式（23）。

假設(shè)溫度和速率均恒定，由式（14）與式（23）可得式（24）。

[lnZ=lnA+nlnsin?ασ] （24）

由式（24）可知[ln[sinh（）]?lnZ]呈線性關(guān)系。以[ln[sin?（）]?lnZ]繪制散點(diǎn)圖，并進(jìn)行線性擬合，如圖7所示。lnA為擬合曲線的截距，由此可得出lnA=17.648 23，則A= 4.618 8×107。

結(jié)合式（14）、（15），根據(jù)雙曲正弦函數(shù)的定義，可將流變應(yīng)力σ定義為含Z參數(shù)的函數(shù)，具體表達(dá)式為式（25）。

為得到各應(yīng)變下的本構(gòu)方程，計(jì)算了各應(yīng)變（0.09-0.14）下的參數(shù)a、n、A、Q與e的關(guān)系。

本文根據(jù)實(shí)際的數(shù)據(jù)情況采用不同的多項(xiàng)式擬合各應(yīng)變下的參數(shù)與應(yīng)變的關(guān)系。擬合所用的多項(xiàng)式如式（26）至式（29）所示，擬合得到的曲線如圖8所示，最終得到的多項(xiàng)式的系數(shù)見(jiàn)表3。

[α=C1ε+C2ε2+C3]" （26）

[n=D1ε+D2ε2?D3ε3+D4ε4+D5] （27）

[Q=E1ε+E2ε2?E3ε3+E4]" （28）

[lnA=F1ε+F2ε2?F3ε3+F4]" （29）

綜上所述，建立的基于應(yīng)變補(bǔ)償?shù)腁rrhenius本構(gòu)模型為：

其中：

基于應(yīng)變補(bǔ)償?shù)腁rrhenius模型流動(dòng)應(yīng)力試驗(yàn)值與預(yù)測(cè)值的對(duì)比如圖9所示。由圖4可知，在不同溫度及應(yīng)變速率下，基于應(yīng)變補(bǔ)償?shù)腁rrhenius模型均具有很好的擬合準(zhǔn)確度，這是因?yàn)榛趹?yīng)變補(bǔ)償?shù)腁rrhenius模型綜合考慮了溫度及應(yīng)變速率和應(yīng)變的耦合效應(yīng)。

2.4 本構(gòu)模型精度評(píng)估分析

本研究分別建立了鎢銅合金的Johnson-Cook本構(gòu)模型、基于應(yīng)變補(bǔ)償修正的Arrhenius模型，并通過(guò)對(duì)比試驗(yàn)及仿真的流變應(yīng)力曲線，分析本構(gòu)方程對(duì)鎢銅合金變形時(shí)的應(yīng)力應(yīng)變預(yù)測(cè)能力的差異。為進(jìn)一步了解不同本構(gòu)方程的預(yù)測(cè)準(zhǔn)確性，可以采取定量計(jì)算的方法，即計(jì)算不同本構(gòu)方程的平均絕對(duì)誤差（AARE）和相關(guān)系數(shù)（R）來(lái)衡量準(zhǔn)確性［25］。平均相對(duì)誤差和相關(guān)系數(shù)的計(jì)算公式分別如式（36）及式（37）。

經(jīng)過(guò)計(jì)算和擬合，本構(gòu)方程相關(guān)系數(shù)和平均相對(duì)誤差見(jiàn)表4，本構(gòu)方程仿真值與試驗(yàn)值的相關(guān)性如圖10所示。由表4可知，相關(guān)系數(shù)最低的為JC本構(gòu)，其值為0.939 14，最高為Arrhenius本構(gòu)模型，其值為0.971 2。這表明兩種本構(gòu)均能夠?qū)︽u銅合金的熱應(yīng)變行為進(jìn)行預(yù)測(cè)。但是Johnson-Cook本構(gòu)模型的平均絕對(duì)誤差較大，為15.682 5，這是因?yàn)樵诟邷丶案邞?yīng)變速率下，Johnson-Cook本構(gòu)模型的仿真值與試驗(yàn)值相差較大。加入溫度敏感系數(shù)可有效提高Johnson-Cook本構(gòu)模型仿真值的準(zhǔn)確性。

從表4中本構(gòu)的平均絕對(duì)誤差和相關(guān)系數(shù)的比較可知，基于應(yīng)變補(bǔ)償?shù)腁rrhenius本構(gòu)模型的相關(guān)系數(shù)最大，平均絕對(duì)誤差更小，擁有更好的精度。如前所述，基于應(yīng)變補(bǔ)償?shù)腁rrhenius本構(gòu)模型綜合考慮了溫度、應(yīng)變速率和應(yīng)變對(duì)流變應(yīng)力的影響。而Johnson-Cook本構(gòu)模型認(rèn)為3者之間的影響是獨(dú)立的，因此對(duì)于仿真參考溫度及應(yīng)變速率差異較大的部分，難以有效預(yù)測(cè)。

3 結(jié)論

本研究對(duì)鎢銅合金的熱變形行為進(jìn)行了研究，通過(guò)應(yīng)力應(yīng)變曲線分析探討了鎢銅合金的高溫變形機(jī)制，建立了 Johnson-Cook本構(gòu)模型、基于應(yīng)變補(bǔ)償?shù)腁rrhenius本構(gòu)模型。最后通過(guò)仿真值和試驗(yàn)值的對(duì)比及定量計(jì)算分析了各種模型的準(zhǔn)確性。具體結(jié)論如下。

①鎢銅合金在高溫拉伸時(shí)有明顯的溫度軟化效應(yīng)和應(yīng)變速率強(qiáng)化效應(yīng)。當(dāng)應(yīng)變速率恒定時(shí)，隨著變形溫度的降低，流變應(yīng)力逐漸增大，相對(duì)應(yīng)的抗拉強(qiáng)度、屈服強(qiáng)度也逐漸增大。當(dāng)變形溫度恒定時(shí)，隨著應(yīng)變速率的降低，鎢銅合金的流變應(yīng)力逐漸降低。

②建立的塑性本構(gòu)方程模型均可對(duì)鎢銅合金的高溫塑性變形行為進(jìn)行預(yù)測(cè)。 Johnson-Cook本構(gòu)模型未考慮溫度和應(yīng)變速率的耦合效應(yīng)，在高溫及高變形速率下，仿真精度降低，相關(guān)系數(shù)為0.939 14，而平均相對(duì)誤差為15.682 5。加入溫度敏感系數(shù)，可顯著提高精確度，相關(guān)系數(shù)提高為0.942 6，平均誤差降低至11.563 6。

③ 基于應(yīng)變補(bǔ)償?shù)腁rrhenius模型的本構(gòu)方程平均絕對(duì)誤差（AARE）和相關(guān)系數(shù)（R）最小，這是因?yàn)樵跀M合時(shí)，考慮了溫度、應(yīng)變及應(yīng)變速率的耦合效應(yīng)，得到的結(jié)果更加接近實(shí)際值。進(jìn)行仿真計(jì)算時(shí)，可以采用基于應(yīng)變補(bǔ)償?shù)腁rrhenius模型的本構(gòu)方程。

參考文獻(xiàn)：

［1］林冰濤， 余小波， 張保紅， 等. 固體火箭發(fā)動(dòng)機(jī)用W8Cu鎢滲銅喉襯抗燒蝕性研究 ［J］. 兵器裝備工程學(xué)報(bào)， 2020， 41（12）： 214-219.

［2］ZHANG B G ， YANG K， HUANG Z F，" et al. Recent Advances in W–Cu Composites： A Review on the Fabrication， Application， Property， Densification， and Strengthening Mechanism ［J］. Advanced Engineering Materials， 2023， 26（1）.

［3］CHEN W G， FENG P， DONG L L， et al. Experimental and theoretical analysis of microstructural evolution and deformation behaviors of CuW composites during equal channel angular pressing ［J］. Materials amp; Design， 2018， 142： 166-176.

［4］SHOKRY A， GOWID S， MULKI H，" et al. On the Prediction of the Flow Behavior of Metals and Alloys at a Wide Range of Temperatures and Strain Rates Using Johnson-Cook and Modified Johnson-Cook-Based Models： A Review ［J］. Materials （Basel）， 2023， 16（4）.

［5］SAVAEDI Z， MOTALLEBI R， MIRZADEH H. A review of hot deformation behavior and constitutive models to predict flow stress of high-entropy alloys ［J］. Journal of Alloys and Compounds， 2022， 903： 163964.

［6］KAREEM S A， ANAELE J U， Aikulola Emmanuel Omosegunfunmi，" et al. Hot deformation behaviour， constitutive model description， and processing map analysis of superalloys： An overview of nascent developments ［J］. Journal of Materials Research and Technology， 2023， 26： 8624-8669.

［7］鄭春曉， 范景蓮， 龔星， 等. 細(xì)晶93W-4.9Ni-2.1Fe合金動(dòng)態(tài)本構(gòu)關(guān)系的研究 ［J］. 稀有金屬材料與工程， 2013， 42（10）： 2043-2047.

［8］SHOKRY A GOWID S， MULKI H，" et al. On the prediction of the flow behavior of metals and alloys at a wide range of temperatures and strain rates using johnson–cook and modified johnson-cook-based models： a review ［J］. 2023， 16（4）： 1574.

［9］SHARMA V NAMBURU S A S LALWANI P，" et al. Constitutive modelling and processing map analysis of tungsten heavy alloy （92.5W-5.25Ni-2.25Fe） at elevated temperatures ［J］. International Journal of Refractory Metals and Hard Materials， 2018， 76： 168-179.

［10］ZERILLI F J， ARMSTRONG R W. Dislocation-mechanics-based constitutive relations for material dynamics calculations ［J］. Journal of Applied Physics， 1987， 61（5）： 1816-1825.

［11］NEMAT-NASSER S， ISAACS J B. Direct measurement of isothermal flow stress of metals at elevated temperatures and high strain rates with application to Ta and TaW alloys ［J］. Acta Materialia， 1997， 45（3）： 907-919.

［12］FOLLANSBEE P S， KOCKS U F. A constitutive description of the deformation of copper based on the use of the mechanical threshold stress as an internal state variable ［J］. Acta Metallurgica， 1988， 36（1）： 81-93.

［13］趙艷平， 張立武， 杜曉斌. 熔滲法制備的W-Cu高溫流動(dòng)特性及其對(duì)塑性變形的影響 ［J］. 熱加工工藝， 2008（8）： 39-41，47.

［14］劉勇， 孫永偉， 田保紅， 等. 鎢含量對(duì)W-Cu復(fù)合材料高溫變形行為的影響 ［J］. 中國(guó)有色金屬學(xué)報(bào)， 2012， 22（9）： 2553-2558.

［15］王凌浩， 伍先明， 莫玉梅. W-20Cu復(fù)合材料熱變形行為及本構(gòu)關(guān)系 ［J］. 粉末冶金工業(yè)， 2021， 31（5）： 16-21.

［16］彭付申， 陳鑫， 袁戰(zhàn)偉， 等. W-20Cu復(fù)合材料熱變形行為及應(yīng)變補(bǔ)償本構(gòu)模型 ［J］. 兵器材料科學(xué)與工程， 2021， 44（6）： 41-46.

［17］LIN D G，HAN J S，KWON Y S，et al. High-temperature compression behavior of W-10wt.%Cu composite ［J］. International Journal of Refractory Metals and Hard Materials， 2015， 53： 87-91.

［18］劉輝明.細(xì)晶W-Cu合金高溫力學(xué)性能與動(dòng)態(tài)力學(xué)行為的研究［D］.長(zhǎng)沙：中南大學(xué)，2011.

［19］馬竇琴.細(xì)晶高致密鎢銅復(fù)合材料制備及電接觸性能研究［D］.鄭州：鄭州大學(xué)，2016.

［20］全國(guó)鋼標(biāo)準(zhǔn)化技術(shù)委員會(huì).金屬材料" " 拉伸試驗(yàn)" " 第2部分：高溫試驗(yàn)方法：GB/T 228.2—2015" ［S］.北京：中國(guó)標(biāo)準(zhǔn)出版社，2015

［21］ZHANG H M， WANG J， CHEN Q，" et al. Study of dynamic recrystallization behavior of T2 copper in hot working conditions by experiments and cellular automaton method ［J］. Journal of Alloys and Compounds， 2019， 784： 1071-1083.

［22］胡海.液相燒結(jié)90W-7Ni-3Fe高密度鎢合金組織性能研究［D］.哈爾濱：哈爾濱工業(yè)大學(xué)，2017.

［23］王偉， 杜帥帥， 閆華軍， 等. 熱壓燒結(jié)Zr-6Al-0.1B合金Johnson-Cook本構(gòu)模型的建立與精度分析 ［J］. 稀有金屬材料與工程， 2023， 52（6）： 2111-2117.

［24］ZENER C， HOLLOMON J. Herbert %J Journal of Applied Physics. Effect of Strain Rate Upon Plastic Flow of Steel ［J］. 1944， 15： 22-32.

［25］HARIKRISHNA K， BHOWMIK A， DAVIDSON M J，" et al. Evaluation of constitutive equations for modeling and characterization of microstructure during hot deformation of sintered Al-Zn-Mg alloy ［J］. Journal of Materials Research and Technology， 2024， 28： 1523-1537.