數(shù)字圖像相關(guān)的不銹鋼真實應(yīng)力應(yīng)變關(guān)系測量

白鵬翔, 倪英薦, 雷 冬

(河海大學(xué)力學(xué)與材料學(xué)院，南京 210098)

金屬材料在工業(yè)鄰域和科研方面一向扮演著極其重要的角色，尤其是近些年隨著建筑工程行業(yè)的蓬勃發(fā)展，各種鋼結(jié)構(gòu)建筑以及各種建筑結(jié)構(gòu)里的鋼構(gòu)件中不銹鋼的使用量大幅上升，對不銹鋼材料的性能和力學(xué)行為的研究也越來越深入[1]。

應(yīng)力應(yīng)變的本構(gòu)關(guān)系是不銹鋼材料的基本力學(xué)性能之一，決定著材料受載后變形的力學(xué)行為。工程上廣泛應(yīng)用的應(yīng)力應(yīng)變曲線實際上是一種近似的名義應(yīng)力應(yīng)變曲線[2-6]，可以近似反映材料在彈性階段的應(yīng)力應(yīng)變關(guān)系，但不能真實反映材料在塑性階段的本構(gòu)關(guān)系。當(dāng)材料在拉伸過程中，對應(yīng)變的實時測量通常采用應(yīng)變片、引伸計、非接觸光學(xué)測量方法等技術(shù)。隨著載荷的增大應(yīng)變也在不斷提高，一方面由于泊松效應(yīng)試件的橫截面會產(chǎn)生變化，使名義應(yīng)力應(yīng)變曲線與真實的應(yīng)力應(yīng)變曲線有所區(qū)別，如王元清等[7]使用引伸計的方式對鋼材料的真實應(yīng)力應(yīng)變曲線進(jìn)行測量；另一方面由于明顯的局部頸縮現(xiàn)象，截面的變化會更加劇烈，此外局部頸縮現(xiàn)象還會使應(yīng)力和應(yīng)變不再沿著試件均勻分布，所以不能再使用引伸計測量試件的伸長量來估計在頸縮區(qū)的真實應(yīng)變。引伸計與應(yīng)變片均為接觸式測量技術(shù)，在安裝和使用過程中不可避免地會對試件產(chǎn)生一定的損傷或限制變形，產(chǎn)生額外的誤差。在各種非接觸測量技術(shù)中，只有光彈法可以直接對應(yīng)力進(jìn)行測量，王煦等[8]卻由于缺乏相應(yīng)的變形數(shù)據(jù)而無法獲得完整的真實應(yīng)力應(yīng)變關(guān)系曲線[9-10]。因而最好使用一種無損、非接觸測量技術(shù)對試件變形過程進(jìn)行檢測以及對變形場進(jìn)行計算。

針對上述情況，在選取不銹鋼標(biāo)準(zhǔn)板試件作為試驗材料的基礎(chǔ)上，選擇了一種無損、高精度的非接觸變形測量技術(shù)——數(shù)字圖像相關(guān)技術(shù)來對試件加載過程中變形場數(shù)據(jù)進(jìn)行測量。在變形場中選取目標(biāo)點組成光學(xué)引伸計對試件的橫向和縱向應(yīng)變進(jìn)行測量，并根據(jù)相應(yīng)的實驗結(jié)果計算得到真實應(yīng)力應(yīng)變曲線。最后通過擬合真實應(yīng)力應(yīng)變曲線，并結(jié)合Johnson-Cook材料模型，獲取材料的彈性及塑性本構(gòu)關(guān)系。

1 數(shù)字圖像相關(guān)測量技術(shù)

數(shù)字圖像相關(guān)(digital image correlation, DIC)技術(shù)[11-15]是結(jié)合了計算機(jī)視覺技術(shù)的一種圖像測量方法，主要是采用現(xiàn)代光學(xué)手段對被測樣品表面全場進(jìn)行測量。

該方法是在20世紀(jì)80年代由日本和美國學(xué)者提出的一種光學(xué)測量技術(shù)，是利用被測試件表面隨機(jī)分布或者人工制作的無規(guī)則特征點作為定位標(biāo)記，通過圖像相關(guān)性函數(shù)來匹配試件表面變形前后的散斑特征，跟蹤表面幾何點的運動來得到位移信息，并計算不同時刻特征點的位置以得到全場各點的位移場。相應(yīng)的全場應(yīng)變則是對位移場數(shù)據(jù)分析而得。理論上應(yīng)變精度可以到達(dá)大約100個微應(yīng)變量級。

數(shù)字圖像相關(guān)技術(shù)中最基礎(chǔ)也是最核心的是相關(guān)函數(shù)，相關(guān)函數(shù)是評價圖像之間相關(guān)性的最直接的手段，學(xué)者對各種不同的相關(guān)函數(shù)進(jìn)行了研究和發(fā)展[16]。使用式(1)中的歸一化最小平方距離相關(guān)函數(shù):

(1)

數(shù)字圖像相關(guān)技術(shù)的優(yōu)點有非接觸、無損檢測、全場表面變形測量、精度高、自動化程度高等，現(xiàn)在被廣泛運用在材料科學(xué)、生物醫(yī)學(xué)、土木交通等各個領(lǐng)域中[17-21]。

2 實驗過程

對試件進(jìn)行拉伸實驗來測試該材料的力學(xué)性能。在試驗中，拉伸試驗機(jī)為Instron Bluehill 3367型號，采用的圖像采集設(shè)備為3D-DIC非接觸測量系統(tǒng)，使用美國CSI公司vic-3d測量系統(tǒng)，裝置包含2個point grey的CCD相機(jī)及高分辨率鏡頭，2個照明光源，計算機(jī)控制系統(tǒng)，以及Vic-3D圖像相關(guān)后處理軟件，現(xiàn)場裝置如圖1所示。

圖1 實驗裝置

試件采用鋼制標(biāo)準(zhǔn)拉伸試件，如圖2所示，其長度為190 mm，拉伸段長度為100 mm，寬度為20 mm，夾持段寬度26 mm，厚度為1.5 mm。

圖2 試件尺寸

實驗過程中采取準(zhǔn)靜態(tài)、勻速加載的方式，加載模式為按照位移加載，拉伸速率設(shè)置為5 mm·min-1，圖像采集速率為3 s·幀-1。通過軟件記錄時間與荷載，其變化關(guān)系如圖3所示，其中最大荷載分別為21 401.7、21 365.8 N。

圖3 載荷時間曲線

使用CCD相機(jī)拍攝試件的拉伸變形過程，將收集到的圖像通過Vic-3D軟件進(jìn)行處理，顯示出其拉伸過程中的變形狀態(tài)。圖4、圖5分別為兩組實驗試件沿著縱向拉伸方向的應(yīng)變場。

結(jié)合上述的應(yīng)變場與荷載-時間曲線，可以很明顯地看出以下結(jié)果。

(1)在彈性階段，由于材料服從胡克定律，當(dāng)荷載迅速地增加時，試件的形變與荷載成一定的線性關(guān)系。

(2)在屈服階段，斜面上的最大切應(yīng)力超過了材料的極限值，此時隨著變形迅速發(fā)展，即試件的變形很快地進(jìn)入到塑性變形階段。

(3)在塑性變形階段，隨著荷載的增加，試件的變形發(fā)展緩慢，在這個過程中，試件的形變量是顯著增加的。

(4)當(dāng)荷載施加到最大值附近時，試件的橫截面顯著縮小，出現(xiàn)“頸縮”現(xiàn)象，試件繼續(xù)變形所需的荷載迅速減小，直到試件在最小橫截面處斷裂。

(5)在斷裂前，從圖5(e)中可以看出，最大應(yīng)變區(qū)在頸縮處成斜向分布，試件的斷裂面并非整齊的切面，并且與橫向成一定角度，這表明導(dǎo)致試件斷裂破壞的不僅僅有拉應(yīng)力，還有剪切應(yīng)力的影響。

圖4 試件1的縱向應(yīng)變場

圖5 試件2的縱向應(yīng)變場

3 真實應(yīng)力應(yīng)變關(guān)系

3.1 應(yīng)力應(yīng)變關(guān)系的獲取

使用數(shù)字圖像相關(guān)(DIC)獲得試件拉伸過程的全場應(yīng)變，通過相關(guān)公式計算，可以獲得試件指定區(qū)域上變形的平均應(yīng)變。采用DIC虛擬引伸計方法，選取圖像上的兩點連線為計算區(qū)域，虛擬引伸計的方向包括橫向與縱向方向(圖6)，計算在該區(qū)域上的平均應(yīng)變，結(jié)果如圖7所示。

圖7 虛擬引伸計的應(yīng)變曲線

在真實變形情況下，試件在斷面處的面積明顯縮小，使得其真實強(qiáng)度往往大于理論上的強(qiáng)度，材料表現(xiàn)出更強(qiáng)的性能。因此對于彈塑性材料來說，需要建立其真實應(yīng)力應(yīng)變關(guān)系來更準(zhǔn)確地描述試件的變形情況及材料性能。

真實應(yīng)力是外加荷載F與某一時刻試件橫截面面積A的比值，真實應(yīng)變是試件某一時刻的伸長變形量ΔL與此時試件長度L的比值。根據(jù)在拉伸變形過程中的體積不變理論，取試件在頸縮前的數(shù)據(jù)進(jìn)行計算，公式如下：

(2)

式(2)中：t為當(dāng)前時刻；εt與σt為當(dāng)前時刻的應(yīng)變與應(yīng)力；L0與A0為初始長度與截面積；εe與σe為名義應(yīng)變與名義應(yīng)力。兩次試驗的應(yīng)力應(yīng)變曲線如圖8所示，取第二次試驗的數(shù)據(jù)進(jìn)行進(jìn)一步的分析。

圖8 兩試件的應(yīng)力應(yīng)變曲線

3.2 彈性模量

通過胡克定律，引入比例常數(shù)E，可得到：

(3)

式(3)中：ε為軸向線應(yīng)變；σ為橫截面上的正應(yīng)力；E為彈性模量。

彈性模量E是表征材料抵抗彈性變形的能力，在真實應(yīng)力應(yīng)變曲線中，取試件在彈性階段的數(shù)據(jù)進(jìn)行擬合，結(jié)果如圖9所示，可以得到彈性模量E=169 GPa。

圖9 彈性模量的獲得

從試件的拉伸情況來看，該材料并沒有明顯的屈服階段，故選取當(dāng)試件產(chǎn)生0.2%塑性應(yīng)變時的應(yīng)力為規(guī)定非比例伸長應(yīng)力，從圖9可得屈服應(yīng)力為236.8 MPa。

3.3 Johnson-Cook模型材料參數(shù)的測定

Johnson-Cook材料模型是一種能反映應(yīng)變率強(qiáng)化效應(yīng)和溫升軟化效應(yīng)的理想剛塑性強(qiáng)化模型，其表達(dá)式為

(4)

在靜載拉伸試驗中，不考慮應(yīng)變率及溫度效應(yīng)，上述模型可以簡化為

(5)

式(5)中：A、B、n為模型參數(shù)；σe為等效應(yīng)力，即材料的真實應(yīng)力；εp為材料的真實塑性應(yīng)變，當(dāng)塑性應(yīng)變εp=0時，A=σs即材料的屈服強(qiáng)度。

材料的變形由彈性變形及塑性變形兩部分組成，其中彈性變形滿足胡克定律。因此，試件的真實塑性應(yīng)變可由式(6)算出：

(6)

式(6)中:mE為應(yīng)力應(yīng)變曲線在彈性階段擬合的斜率。

將兩次試驗的數(shù)據(jù)進(jìn)行分析和擬合，結(jié)果如圖10所示。

圖10 塑性變形應(yīng)力應(yīng)變關(guān)系擬合

擬合可以得到相關(guān)參數(shù)為：A=236.8 MPa;B=1 825 MPa;n=0.89。

按照式(5)給出不銹鋼材料的真實塑性應(yīng)力應(yīng)變關(guān)系為

(7)

4 結(jié)論

使用數(shù)字圖像相關(guān)技術(shù)進(jìn)行標(biāo)準(zhǔn)試件的拉伸變形量測，得到其斷裂前的位移和應(yīng)變信息。對數(shù)據(jù)的處理與分析，獲取材料的力學(xué)特性，得到主要結(jié)論如下：

(1)得到了整個拉伸變形過程。在彈性階段，由于材料服從胡克定律，試件的形變與荷載成一定的線性關(guān)系；在屈服階段，斜面上的最大切應(yīng)力超過了材料的極限值，此時隨著變形迅速發(fā)展，即試件的變形很快地進(jìn)入到塑性變形階段；當(dāng)荷載施加到最大值附近時，試件的橫截面顯著縮小，出現(xiàn)“頸縮”現(xiàn)象，使得試件繼續(xù)變形所需的荷載迅速減小，直到試件在最小橫截面處斷裂。

(2)獲取真實應(yīng)力應(yīng)變關(guān)系。用虛擬引伸計計算平均應(yīng)變，再根據(jù)體積不變理論計算真實應(yīng)力應(yīng)變；用Johnson-Cook模型擬合材料的塑性本構(gòu)關(guān)系。

(3)由真實應(yīng)力應(yīng)變曲線提取出材料的力學(xué)特性參數(shù)，包括彈性模量，屈服應(yīng)力等。