王晨羽,顏世晶,李金泉
(沈陽理工大學 機械工程學院,沈陽 110159)
材料的本構關系模型是反映材料物理力學行為的數(shù)學模型,描述材料應力與應變、應變率與溫度等變量之間的關系。目前,針對不同材料和不同實驗條件已經(jīng)得出了很多本構關系。Fields等提出了一種冪指數(shù)強化本構模型,求解過程比較簡便,但在高溫條件下其冪指數(shù)增長的特性不能準確描述材料的流變行為[1]。文獻[2]提出的Bodner-Partom本構模型,針對材料彈性和塑性變形特性,描述不同材料的力學行為,但參數(shù)較多,不宜應用。文獻[3]提出的Follansbee-Kocks本構模型,認為應變和初始加工硬化率線性相關,但在高應變率下難以獲得精確的模擬結果,同樣也存在材料參數(shù)過多,使用不方便等缺點。Zerilli-Armstrong本構模型是在材料微觀結構的位錯力學模型基礎上,增加了晶粒尺寸對本構的影響[4]。Johnson-Cook(J-C)模型[5]由于適合大多數(shù)金屬材料,且結構簡單便于應用,因此,獲得了廣泛的應用,而且針對不同材料產生的誤差分別修正了應變強化項[6]、應變率強化項[7]和溫度軟化項[8-10],得到了改進型的J-C本構模型,使J-C本構模型的應用更加廣泛,因此,本文根據(jù)霍普金森壓桿(SHPB)實驗研究黃銅的J-C本構關系模型,實驗表明:黃銅的擬合曲線與實驗曲線吻合較好,相對誤差小于9.38%。
試驗材料為黃銅,進行SHPB試驗,試樣為圓柱形,尺寸為Φ5mm×5mm,所用子彈長度為140mm,有效氣壓為0~ 0.60Mpa,溫度范圍0~ 800℃,應變率范圍2000 s-1~ 8000 s-1,將相同的條件進行多次重復試驗,取三組有效數(shù)據(jù),以此來確保試驗結果的可靠性。試樣成分見表1。
表1 黃銅的化學成分 wt/%
J-C熱塑性本構方程表達式為
(1)
σ(ε)=A+Bεn
(2)
根據(jù)常溫、準靜態(tài)下的σ-ε試驗數(shù)據(jù),當ε=0時的應力值即為參數(shù)A,用Matlab軟件對數(shù)據(jù)進行擬合,得到參數(shù)B、n的值。
常溫下(T=T0),式(1)簡化為
(3)
上述參數(shù)求得后,由式(1)可得
(4)
依據(jù)相應溫度下的一系列σ-ε數(shù)據(jù),設定一些應變和應變率為已知,可得相對應的應力值表達式,式(4)轉變?yōu)橐訲*和m為變量的函數(shù),改變一系列溫度,則可擬合出m值。
通過上述計算,得黃銅的本構關系為
(5)
圖1為黃銅在5500s-1應變率下,溫度分別為25℃、270℃、570℃和790℃時J-C本構模型的擬合曲線與真實應力應變曲線的比較,二者吻合良好。
圖1 黃銅的J-C模型曲線和真實應力-應變曲線
為便于比較,將圖1中不同溫度下的真實應力應變曲線繪在同一張圖中,如圖2所示。
由圖2中可以看出,在較低的溫度范圍內(25℃~ 270℃),應力隨應變的增大而增大,應力變化幅度類似,表現(xiàn)出相似的應變硬化效應,且在該溫度范圍內應力近似相等,熱軟化效應不明顯;當溫度增加到570℃時,應力降低幅度較大,熱軟化效應明顯,但應力仍隨應變的增大而增大,仍有應變硬化效應;溫度繼續(xù)增大到790℃時,應力繼續(xù)降低,且降低幅度較大,熱軟化效應明顯,但應力隨應變增大的變化很小,說明高溫時黃銅的應變硬化效應較小。圖3為不同應變下應力隨溫度變化的曲線。
圖2 不同溫度下真實應力-應變曲線
圖3 不同應變下溫度對黃銅H62應力的影響
由圖3可知在較低溫度范圍內熱軟化效應很小;當溫度增加到570℃和790℃時,應力降幅非常大,表現(xiàn)出明顯的熱軟化效應;而相同溫度下,應變對應力的影響很小,因此,在高溫環(huán)境下實際應用時,黃銅的工作溫度不宜超過270℃。
為更清楚的反映擬合本構模型的準確性,隨機選取每種實驗條件下的五組數(shù)據(jù),分析其相對誤差及誤差平均值,如表2 所示。
表2 本構模型預測值與實驗值的誤差
由表2可以看出,表中相對誤差均小于10%,同時,黃銅合金不同溫度下的平均誤差均在5.34%以下。通過誤差分析可知,本實驗擬合的本構關系準確度較高。
(1)利用SHPB裝置對黃銅(H62)進行了動態(tài)壓縮試驗,分析了黃銅的動態(tài)破壞特性;
(2)測試了黃銅的J-C模型曲線和真實應力-應變曲線,擬合得到了黃銅的J-C本構模型,黃銅的擬合曲線與實驗曲線吻合較好,相對誤差小于9.38%。
(3)黃銅(H62)在25℃~570℃范圍內,有明顯的應變硬化效應,25℃~270℃范圍內幾乎沒有熱軟化效應,而在570℃~790℃范圍內有較大的熱軟化效應。