動態斷裂韌性實驗中DIC技術應用研究

曲　嘉，李東昌，黃　超

（1.哈爾濱工程大學，黑龍江 哈爾濱 150001；2.哈爾濱理工大學，黑龍江 哈爾濱 150080）

動態斷裂韌性實驗中DIC技術應用研究

曲嘉1，李東昌1，黃超2

（1.哈爾濱工程大學，黑龍江 哈爾濱 150001；2.哈爾濱理工大學，黑龍江 哈爾濱 150080）

高速沖擊動態斷裂韌性的加載和測試技術一直是近年關注的熱點，隨著計算機和光學傳感器的發展，采用數字圖像相關方法測量材料的動態斷裂韌性已成為重要選擇。該文基于分離式Hopkinson壓桿原理的加載技術，通過高速攝影機拍攝高速沖擊下三點彎曲試樣裂紋的起裂和擴展，運用數字圖像相關技術分析裂尖場的散斑圖像，計算得到相應的應變場變化，試樣外表面處于平面應力狀態，其裂尖應變場呈現“0”型，而非試樣對稱面所處于的平面應力狀態下呈“8”型。結果表明，DIC技術可以應用于動態斷裂韌性實驗，也證明裂尖場粘貼的應變片測量試樣起裂的有效性。

數字圖像相關；動態斷裂韌性；分離式霍普金森壓桿；不銹鋼

0　引　言

斷裂韌性是材料抵抗破壞的重要參數，隨著斷裂力學工程應用的逐步深入，已成為研究者關注的熱點［1］。三點彎曲試樣作為測定準靜態載荷作用下的標準試件已經得到了廣泛應用，但由于動態斷裂題的復雜性，目前還沒有標準的測試方法。盡管如此，由于三點彎曲試樣比較簡單，加載也很方便，因此采用三點彎曲試樣，并通過分離式霍普金森壓桿技術進行動態斷裂韌性實驗研究。

數字圖像相關測量方法是由美國 Peters和Ranson教授及日本的Yamaguchi在20世紀80年代初期共同獨立提出的［2-3］。隨著數字圖像相關識別技術和計算機數字圖像處理技術的不斷進步和發展，數字圖像相關技術成為一種運用計算機對采集圖像進行數字化分析的技術［4］。本文采用高速攝影與數字圖像相結合，監測試樣的起裂和裂紋擴展，獲得在裂紋尖端所引起的應變場變化。結果分析驗證了裂尖場粘貼的應變片測量試樣起裂的有效性。數字圖像相關方法具有試驗環境易于實現、測量范圍易于控制和試驗結果圖形化等諸多優點，能滿足實驗力學及其工程測量的實際需要。

1　動態斷裂韌性實驗

1.1試樣疲勞預制裂紋和噴涂散斑

本次實驗設備為分離式Hopkinson壓桿，用不銹鋼材料制備含有預制裂紋的三點彎曲試樣，試樣的厚度為10 mm，高度為20 mm，跨距為80 mm。使用Instron-8801電液伺服試驗機對試樣以正弦波的方式進行加載，預制后的裂紋總長度為10mm。在預制裂紋的試件一面噴涂隨機散斑，處理完成的試件如圖1所示。

圖1　預制裂紋和噴涂散斑試件圖

1.2霍普金森壓桿技術和實驗裝置

分離式霍普金森壓桿（SHPB）原型是由Hopkinson于1914年提出，用于測量載荷的脈沖波形。1949年，Kolsky提出將壓桿分成兩段，試樣置于兩桿中間，從而使這一裝置可用于測量材料在沖擊載荷作用下的應力-應變關系［5-6］。SHPB試驗裝置示意圖如圖2所示，通過炮筒內發射的短桿子彈以一定速度撞擊入射桿，再由入射桿對三點彎曲試樣進行加載，短桿子彈的速度由激光分束測速裝置測得。在試樣裝夾處用高速攝影機拍攝試樣整個斷裂過程，以便使用數字圖形相關技術分析起裂和裂紋擴展過程，圖3為SHPB試驗裝置圖和三點彎曲試樣裝夾圖。

圖2　SHPB試驗裝置圖

圖3　試件裝夾圖

1.3試驗過程及結果

實驗過程中，先對桿系進行固定，利用桿系的調節裝置對桿系進行調節，使入射桿與透射桿保持同軸并且在同一水平面上。在入射桿和透射桿上分別粘貼兩組應變片，粘貼時保證每根桿上的兩個應變片在桿的同一截面相對稱的位置，這樣可以消除由于桿的彎曲變形對應變信號測量所產生的影響。通過這兩組應變片測量入射桿和透射桿的應變信號。在對引線以及屏蔽線的焊接過程中，一定要保證連接點要光滑，最大程度的消除人為因素對試驗結果產生的影響。使用高速攝影機進行圖像采集，采集到的數字圖像如圖4所示，并存入計算機。

圖4是使用高速攝影機采集到的數字圖像，（a）～（d）分別為激發采集后不同時間下三點彎曲試樣的表面圖像，可以看出在激發后458.302 ms時表面首先開始開裂，在此之前沒有開裂，本文使用的高速攝影機前一幀為458.227ms，選取此時的試樣圖像為原始圖像，開裂后再經過0.025 ms，裂紋擴展已經可以從表面明顯看出，在激發后458.401 ms時裂紋擴展已經特別明顯，取這幾張有代表性的圖像作為動態斷裂韌性試驗的結果。

2　數字圖像相關方法分析裂紋擴展

2.1數字圖像相關基本原理

數字圖像相關方法的基本原理［10-14］是通過比較試件變形前后物體表面分布的散斑圖的變化，來獲得位移和應變等變形信息。數字圖像相關實現過程是在變形前的圖像中，以所要計算的點為中心選取子區，利用子區中的灰度場，在變形后的圖像中，通過相關運算找出最值的目標圖像的對應區域，變形前后對應以確定該點的位移值。再通過亞像素插值方法獲得連續位移值，計算位移梯度獲得應變值［7-9］。

圖4　動態三點彎曲試樣起裂和裂紋擴展圖

圖5　數字圖像相關原理圖

相關運算的核心是相關公式的選擇，本文采用的相關系數為

式中：r——參考圖像灰度大??；

d——變形后圖的灰度大??；

（x′，y′）——原始圖像中心點的坐標；

（u′，ν′）——理論位移，計算最大相關系數對

應的位移是（u，ν）。

采用Matlab編制的數字圖像相關程序對采集到的圖像進行數據處理［10］。

2.2數字圖像技術分析

圖6為試樣噴涂散斑區域和試樣計算區域，該區域也是實驗中貼片的位置，通過分析數字圖像，采用有限元法對位移插值獲得應變場變化。

圖6　圖像采集區域和DIC計算區域

圖7中（a）為激發后的時間458.302 ms時的起裂時刻產生的應變，其中最大應變為1.4%，隨著裂紋擴展，激發后的時間458.327ms時產生的最大應變為1.6%，最后458.401ms時裂紋擴展到很大時相較于原始圖像產生30%應變。從裂紋擴展過程可以看出金屬材料動態斷裂過程中，表面上裂紋是沿著裂紋線開裂，裂紋尖端前移導致應變區域向裂紋擴展方向移動。試樣外表面處于平面應力狀態，其裂尖應變場的呈現“0”型，而非試樣對稱面所處于的平面應力狀態下呈“8”型，由于本文測試的圖像為試樣表面，因此，計算顯示應變場是呈現“0”型擴展的。

圖7　動態三點彎曲試樣應變場變化圖

3　結束語

本文基于分離式Hopkinson壓桿原理的加載技術，通過高速攝影機拍攝高速沖擊下裂紋的起裂和擴展，運用數字圖像相關技術分析裂紋擴展圖像，獲得在裂紋尖端所引起的應變場變化。結果表明，DIC技術可用于動態斷裂韌性實驗，測得處于平面應力狀態的試樣外表面裂尖應變場呈現“0”型，而非試樣對稱面所處于的平面應力狀態下呈“8”型，也證明了裂尖場粘貼的應變片測量試樣起裂的有效性。同時，從裂紋擴展過程可以看出金屬材料動態斷裂的過程是沿著裂紋線開裂，由于處于平面應力狀態，應變場是呈現環狀擴展的。據此，DIC技術應用于動態斷裂韌性可為工程中測量材料的斷裂特性提供一種實用可靠的方法。

［1］許澤建，李玉龍.高強鋼40Cr和30CrMnSiNi2A I型和II型動態斷裂韌性的研究［J］.金屬學報，2006，42（9）：965-970.

［2］YAMAGUCHI I.Speckle displacement and deformation in the diffraction and image fieldsfor small object deformation［J］.Acta Optica Sinica，1981，28（10）：1359-1376.

［3］SUTTON M A，WOLTERS W J，PETERS W H，et al. Determination of displacementsusing an improved digital correlation method［J］.Image and Vision Computing，1983，1（3）：133-139.

［4］高建新.數字散斑相關方法及其在力學測量中的應用［D］.北京：清華大學，1989.

［5］CHRISTENSEN R J，SWANSON S R，BROWN W S. Split-Hopkinson-Bar test on rocks under confining pressure［J］.Experimental Mechanics，1972（11）：508-513.

［6］宋博，姜錫權，陳為農.霍普金森壓桿實驗中的脈沖整形技術［C］∥第三屆全國爆炸力學實驗技術交流會，2004.

［7］芮嘉白，金觀昌，徐秉業.一種新的數字散斑相關方法及其應用［J］.力學學報，1994，26（5）：599-607.

［8］潘兵，謝惠民.數字圖像相關中基于位移場局部最小二乘擬合的全場應變測量［J］.光學學報，2007，27（11）：1980-1986.

［9］RAND J L，GRANT D A．Optical measurement of biaxialstraininehinfilmpolymers［J］.ASTMSpecial Technical Publication，1997，1318（8）：123-137.

［10］周曉峰.基于散斑數字相關的平面全場應變測量方法及應用［D］.廣州：華南理工大學，2012.

［11］張少軍，艾矯健，李忠富，等.利用數字圖像處理技術測量幾何尺寸［J］.北京科技大學學報，2002，24（3）：284-287.

［12］孫偉，何小元，胥明.數字圖像相關方法在膜材拉伸試驗中的應用［J］.工程力學，2007，24（2）：34-38.

［13］熊磊.基于數字圖像相關方法的材料拉伸變形測量技術［D］.長春：吉林大學，2012.

［14］PANB，QIANKM，XIEHM，etal.Twodimensionaldigitalimagecorrelationforin-plane displacementandstrainmeasurement：Areview［J］. Measurement Science&Technology，2009，20（6）：152-154.

（編輯：徐柳）

Application of digital image correlation method in dynamic fracture toughness testing

QU Jia1，LI Dongchang1，HUANG Chao2
（1.Harbin Engineering University，Harbin 150001，China；2.Harbin University of Science and Technology，Harbin 150080，China）

The high speed dynamic fracture toughness loading and testing technology has become the focus of the public.With the development of computer and optical sensors，measuring the dynamic fracture toughness of materials with digital image correlation method has become an importantchoice.BasedontheprincipleofsplitHopkinsonpressurebar（SHPB）loading technique，a high-speed camera was used to photograph crack initiation and propagation process of three-point bending specimen under high-speed shocking，and analyze the speckle image of crack tip by digital image correlation method.The developing processing of strain field was obtained.In conclusion，the outside surface of the specimen is in plane stress state and the strain field of the crack tip is in“0”shape instead of the"8"shape.The results showed that the DIC technique could be applied to the dynamic fracture toughness test，and it also proves the effectiveness of the strain gage measuring the crack initiation.

digital image correlation；dynamic fracture toughness；SHPB；stainless steel

A

1674-5124（2016）10-0045-04

10.11857/j.issn.1674-5124.2016.10.009

2016-05-03；

2016-06-19

國家自然科學基金（11302055）哈爾濱市科技創新人才專項基金（2013RFQXJ009）

曲嘉（1979-），男，黑龍江哈爾濱市人，副教授，博士，研究方向為沖擊動力學及強度理論。