盧伍鋆

（福建省特種設備檢驗研究院龍巖分院，福建 龍巖 364000）

0 引言

拉伸強度是最基本的力學計算問題，其測定方法通常是在拉伸試驗機進行多次和小速度的拉壓實驗得到。拉伸試驗機測量拉伸強度雖然結(jié)果比較可靠，但是存在制樣過程耗時成本高和操作繁瑣問題。

隨著計算機硬件技術(shù)提高和有限元理論成熟，應用CAE軟件對此類實驗過程進行仿真模擬越來越受到人們的青睞[1]。這些CAE軟件的出現(xiàn)正逐步將研究學者們從復雜的實驗過程中解放出來，去研究材料力學性能的數(shù)值模擬，基于工藝模擬系統(tǒng)的Deform3D軟件是一種有效處理三維金屬流動問題的計算軟件。

有限元方法首先需要把研究對象離散化[2]，劃分為許多個細小的單元，在單元之間使用節(jié)點連接起來，通過數(shù)值迭代求解，從而得到整體物理量分布規(guī)律；有限元方法的應用大大地減少了研究時間，縮短了產(chǎn)品試樣的制造工時，降低了產(chǎn)品的生產(chǎn)成本。

為降低實驗試樣制作成本，文中采用對比方法首先進行低碳鋼試樣拉伸實驗和數(shù)值模擬，驗證數(shù)值模擬的正確性；接著進行15鋼拉伸試樣的數(shù)值模擬，研究其等效應力、拉伸力問題，為企業(yè)提供相應的力學參數(shù)。

1 拉伸實驗與數(shù)值模擬對比

圖1為某材料力學拉伸實驗，實驗棒如圖1（a）所示，拉伸試驗機如圖1（b）所示，實驗棒尺寸（Q235鋼）如圖1（c）所示，上下夾頭使用的拉伸速度為0.2 mm/min，從拉伸開始到強化階段整個實驗力與位移的變形過程曲線如圖1（d）所示。從圖中可以看出，拉伸時間約為6.22 s時，實驗力為37.10 kN。

圖1 材料力學拉伸實驗

如圖2所示，采用Deform3D數(shù)值軟件對圖1（b）試驗棒以實際拉伸速度0.2 mm/min進行數(shù)值分析。為提高分析精度，試驗棒中部采用加密網(wǎng)格（如圖2（a）），分析得到拉伸開始到強化階段整個實驗力與位移的變形過程曲線，如圖2（b）所示。從圖中可以看出，拉伸變形約為5.0 mm（約為6.22 s）時，實驗力最大約為34.6 kN，該數(shù)值比實際試驗37.10 kN少了2.5 kN，誤差約為6.7%，誤差存在原因是數(shù)值模擬采用了網(wǎng)格簡化模型和忽略夾頭摩擦因數(shù)等。

圖2 碳鋼的數(shù)值模擬

2 材料本構(gòu)關(guān)系

15鋼的彈性模量E=209 GPa，泊松比ν=0.28，密度ρ=7.85 g/cm3；由于拉伸試驗是在室溫下進行的，屬于準靜態(tài)過程，所以只考慮了應變率效應、應變硬化的影響，溫度忽略不計。

考慮應變、應變速率影響的Johnson-Cook（J-C）強度模型為

3 實驗試樣建模及相關(guān)參數(shù)設置

圖3為15鋼的標準拉伸試樣尺寸，根據(jù)該尺寸首先建立拉伸試樣和夾具三維模型并進行裝配；然后導入Deform軟件，進行多組拉伸速度的數(shù)值仿真。表1所示為給定的4組速度。

表1 拉伸速度mm/s

圖3 標準拉伸扭轉(zhuǎn)試樣零件圖

網(wǎng)格Size Ratio設置為1，網(wǎng)格最小尺寸設為1.5，定義接觸將上下模和試件固連到一起[6]，如圖4所示。

圖4 拉伸扭轉(zhuǎn)有限元模型

4 數(shù)值模擬結(jié)果分析

4.1 不同拉伸速度下的等效應力

在完成前處理相關(guān)參數(shù)設置后，提交FEM求解器計算，得到了拉伸試樣在不同速度下彈性、塑性變形過程。圖5（a）～圖5（d）分別拉伸速度為0.05、0.10、0.30、0.50 mm/s時的等效應力云圖。從圖中可以看出，隨著拉伸速度的增大，最大等效應力也隨之增大，且均出現(xiàn)在試樣中部；當拉伸速度為0.3 mm/s和0.5 mm/s時，材料發(fā)生了屈服。

圖5 拉伸時等效應力及最大等效應力隨速度變化關(guān)系

4.2 拉伸力變化分析

圖6（a）～圖6（d）為不同拉伸速度下拉伸力曲線。當拉伸速度為0.05 mm/s時，其彈性階段的拉伸力迅速上升，最大為24.05 kN，從彈性階段過渡到穩(wěn)定塑性變形階段，其拉力下降，如圖6（a）所示；從圖6（b）可知，當拉伸速度為0.1 mm/s時，其彈性階段的拉伸力快速上升，最大為84.57 kN，從彈性階段過渡到穩(wěn)定塑性變形階段，隨后拉力開始緩慢下降；從圖6（c）可知，當拉伸速度為0.3 mm/s時，其彈性階段的拉伸力緩慢上升，最大為52.3 kN，材料發(fā)生彈性變形、塑性變形、強化過程、頸縮過程；從圖6（d）可知，當拉伸速度為0.5 mm/s時，其彈性階段的拉伸力迅速上升，最大為53.6 kN，材料從彈性變形到塑性變形過程中拉力下降。

圖6 不同拉伸速度下的拉伸力曲線

5 結(jié)語

本文基于拉伸試驗機首先對低碳鋼拉伸試樣進行拉伸實驗，并借助Deform3D軟件對其進行數(shù)值模擬，然后對這兩種方法得到的拉伸位移與實驗力曲線結(jié)果進行對比，結(jié)果表明兩種方法的誤差約為6.7%，在工程許可范圍內(nèi)，驗證數(shù)值模擬的有效性。接著研究了某企業(yè)15鋼拉伸試樣在不同速度下拉伸過程，分析了拉伸時的等效應力、拉伸力等。結(jié)論為：隨著拉伸速度的增大，最大等效應力也隨之增大，且均出現(xiàn)在試樣中部；隨著拉伸速度的增加，其拉伸力也隨著增大。