封偉強，王偉，孟松鶴，許承海，解維華

（哈爾濱工業(yè)大學 復合材料與結(jié)構(gòu)研究所，哈爾濱 150080）

基于數(shù)字圖像相關(guān)方法的石墨高溫楊氏模量測試

封偉強，王偉，孟松鶴，許承海，解維華

（哈爾濱工業(yè)大學 復合材料與結(jié)構(gòu)研究所，哈爾濱 150080）

目的 為了得到石墨在高溫環(huán)境下的楊氏模量，為飛行器熱防護系統(tǒng)和高溫熱結(jié)構(gòu)的設(shè)計提供可靠的技術(shù)保障。方法 基于數(shù)字圖像相關(guān)（DIC）方法和通電電阻加熱技術(shù)，建立一套測量材料在高溫環(huán)境下力學性能的測試系統(tǒng)。利用該系統(tǒng)，測量高溫下石墨試樣的表面應變場和應力-應變曲線，計算相應溫度下的楊氏模量。在1400 ℃實驗環(huán)境下，采用基于DIC與高溫引伸計的方法同時測量超高溫陶瓷試件的拉伸應變數(shù)據(jù)，并進行比較驗證。結(jié)果 在 1400 ℃實驗環(huán)境下，采用基于 DIC與高溫引伸計方法測得超高溫陶瓷的應變-時間曲線吻合良好，方差為 1.3×10-7。1200～1900 ℃高溫環(huán)境下，石墨的楊氏模量隨溫度的升高呈線性增長趨勢。結(jié)論 采用基于DIC方法準確有效，該方法可方便快速地實現(xiàn)對石墨材料在高溫環(huán)境下楊氏模量的測量。另外，該方法也可應用于其他導電材料楊氏模量的測試。

石墨；高溫應變場；數(shù)字圖像相關(guān)方法；楊氏模量；應力-應變曲線

近年來，隨著我國航空航天工業(yè)的迅速發(fā)展，各種再入飛行器和高超聲速飛行器的飛行速度越來越快，其所處的服役環(huán)境也越來越復雜[1]。高超聲速飛行器以高馬赫數(shù)飛行時，其表面溫度可達上千度，因此對熱防護系統(tǒng)的設(shè)計和可靠性評估非常重要[2]。

數(shù)字圖像相關(guān)（Digital Image Correlation）測量方法是20世紀80年代初提出的一種非接觸的變形測量方法。它使用相機分別記錄下被測試件表面變形前后的灰度信息，然后通過軟件處理圖像信息，以得到所需的應變場[3]。DIC方法具有實驗設(shè)備和過程簡單、能實現(xiàn)全場測量、對測量環(huán)境和隔振要求低等優(yōu)點[4]，被廣泛用于各個學科的研究[5]。如 S.Mallon等[6]運用此方法對正交編織玻璃纖維增強復合材料的動態(tài)斷裂行為進行了研究，N.S.Ha等人[7]用該方法對人造仿生甲殼蟲尾翼進行了模態(tài)分析，其他國內(nèi)外學者對該方法在高溫實驗中的應用也進行了大量的研究工作[8—11]。另外，研究人員對誤差影響因素分析、散斑制作改進算法和提高測量精度等方面進行了許多嘗試和改進[12—17]。

文中針對等靜壓石墨材料，采用數(shù)字圖像相關(guān)方法測量了其在 1200～1900 ℃高溫環(huán)境下的單軸拉伸變形信息，結(jié)合載荷數(shù)據(jù)可得到石墨材料在高溫環(huán)境下的應力-應變曲線，進而計算得到不同溫度下的楊氏模量。

1 實驗

1.1 高溫實驗系統(tǒng)

為了測量待測材料的高溫力學性能，文中將數(shù)字圖像相關(guān)（Digital Image Correlation）方法和電阻通電加熱技術(shù)相結(jié)合，建立了一套高溫力學性能測試系統(tǒng)，如圖1所示。該系統(tǒng)由硬件和軟件兩個子系統(tǒng)組成，硬件子系統(tǒng)主要由高溫力學實驗機、Point Gray公司的CCD相機、圖形采集卡、光源和分析計算機組成；軟件子系統(tǒng)主要包括Vic-Snap圖像采集軟件和Vic-2D分析軟件。

公式（1）為普朗克（M. Planck）輻射定律描述的光譜輻射力的表達式，黑體的光譜輻射力隨溫度的升高急劇增長。另外，高的光譜輻射力大部分處于紅光以及紅外波段，這樣會影響測量的準確度，為了消除這些干擾，在相機鏡頭上加設(shè)藍光帶通濾鏡，只允許波段在425～475 nm內(nèi)的藍光通過，同時使用兩個藍色 LED燈補強照明亮度，如圖 1所示。

圖1 高溫力學性能測試系統(tǒng)Fig.1 High-temperature mechanical performance test system

1.2 數(shù)字圖像相關(guān)方法原理

數(shù)字圖像相關(guān)（DIC）測量法的基本原理是以初始圖像為基準，以感興趣的像素點為中心選取一個正方形圖像子區(qū)，并通過追蹤對比該圖像子區(qū)在變形后圖像（或目標圖像）中的位置得到圖像子區(qū)中心點的位移矢量[18]，按同樣的方法對參考圖像中感興趣區(qū)域內(nèi)的多個像素點進行同樣的相關(guān)計算即可獲得該計算區(qū)域的位移場。在實際分析計算時，為獲得精確可靠的匹配結(jié)果，所選擇的圖像子區(qū)應該足夠大以包含充分多的灰度變化信息，從而確保該圖像子區(qū)在變形后的圖像中能被唯一識別。這就要求被測物體表面必須覆蓋有隨機分布的散斑圖，在實驗過程中散斑隨試件表面協(xié)同變形，承載了試件的變形信息。實踐證明，由于散斑的質(zhì)量對數(shù)字圖像相關(guān)法的位移測量精度有很大的影響，因此制作高質(zhì)量的散斑圖在整個實驗過程中尤為重要。試件散斑如圖2所示。

為了表征圖像子區(qū)在變形前后兩幅散斑圖像的相似程度，定義了如式（2）所示的相關(guān)函數(shù)[19]：

圖2 試件散斑Fig.2 Speckles on specimens

1.3 散斑制備技術(shù)

高中電磁學反映的是電場與磁場的相互作用關(guān)系，揭示的是物理宏觀世界的本質(zhì)規(guī)律，抽象的概念與規(guī)律缺乏表象基礎(chǔ)，注重外部表面特征和非本質(zhì)聯(lián)系的高中學生對于可見度不高的電磁現(xiàn)象及電磁規(guī)律顯得很茫然。只有認識公式、理解公式，準確應用公式表述問題，才能避免亂套公式、混淆公式的現(xiàn)象。

為了獲得較為可靠的散斑，實驗前用高標號砂紙將試件標距區(qū)表面進行適度打磨，并清洗干凈，采用高溫耐火涂料手工制作散斑。具體做法是將氧化鋁粉末和膠體按照質(zhì)量比3∶1混合均勻制成耐火涂料，用細毛刷蘸取涂料并對準試件，保持一定的距離，用手指撥動毛刷使涂料飛濺到試件表面，操作時注意控制毛刷和試件的距離以及撥動的力度，即可形成大小適中、分布均勻且非重復的散斑。上述制作過程可重復，直至散斑效果理想為止。一般散斑的大小范圍較為寬泛，但如果散斑太大，某一分析子區(qū)可能為純黑或者純白區(qū)域，不利于VIC-2D軟件分析；如果散斑太小，則相機的分辨率不足以呈現(xiàn)試件表面形貌。因此理想的散斑應在保證采集到的圖像中有至少 3～4個像素點的前提下盡可能小[20]。已制作好散斑的石墨材料高溫拉伸板狀試件如圖3所示。

圖3 石墨材料高溫拉伸板狀試件Fig.3 Graphite material sheet-like tensile specimens at high temperature

為了考察這種散斑的溫度適用范圍，事先分別采集了散斑在不同溫度下的圖像。首先將散斑已干燥完成的石墨試件裝載進實驗艙，將夾具裝載完畢，預加200 N拉伸載荷，使試件與電極之間緊密接觸，避免通電時因接觸不良導致打火；然后使用真空機械泵抽取實驗艙內(nèi)空氣，使其內(nèi)部氣壓降至5 Pa以下，盡量減少高溫時石墨被氧化對圖像采集的影響，使用紅外高溫計通過實驗艙窗口對石墨試件標距段的表面溫度進行實時監(jiān)測，使用溫控設(shè)備進行實時控溫；接通電源開始進行升溫，使用CCD相機采集下不同溫度下散斑的分布圖像，如圖4所示。

圖4 不同溫度下散斑的分布圖像Fig.4 Images of speckle distribution at different temperatures

通過觀察對比發(fā)現(xiàn)，在1700 ℃以下散斑性能良好，基本不會脫落；在1800 ℃左右試件表面的散斑開始剝離，個別附著不牢的散斑首先脫落；若溫度繼續(xù)升高，散斑迅速剝離，大約到1900 ℃時大尺寸散斑基本已剝離，僅剩小尺寸散斑可辨識。因此，此種方法制備的散斑可用于1900 ℃以下的實驗環(huán)境。

經(jīng)過分析，散斑在高溫下剝離的可能原因為：升溫引起的熱膨脹導致試件和散斑的熱變形不匹配；材料高溫下發(fā)生氧化、熱解；散斑點厚度及尺寸不均勻造成附著力略有偏差。

結(jié)合以上原因，針對于不同材料制作高溫散斑，應充分考慮到材料和散斑在高溫下的特性，采取可靠的工藝加以改進，得到質(zhì)量高、性能可靠的散斑，這需要在今后的研究工作中加以改進。

高溫拉伸實驗過程主要分為三段，1200 ℃實驗過程如圖5所示。首先是預加載階段（圖5階段I），室溫下以10 N/s的加載速度將載荷加至200 N，保證在接下來的通電加熱和保溫過程中試件與上下主電極之間的良好接觸，防止出現(xiàn)打火現(xiàn)象；然后是升溫階段（圖5階段II），保持200 N預載的同時，將試件溫度升至目標溫度，紅外比色計與控溫軟件配合可以實現(xiàn)20 ℃/s的升溫速率；最后是高溫加載階段（圖5階段III），試件溫度達到目標溫度后，即開始以 1 mm/min的速率對試件進行加載，同時使用CCD相機進行試件表面散斑圖像的采集。

圖5 高溫拉伸實驗過程Fig.5 Schematic diagram of high-temperature tensile testing

2 結(jié)果與討論

2.1 準確性驗證

該實驗選用的CCD相機為Point Gray公司的500萬像素相機，最大分辨率為2448(H)×2048(V)，像素大小為4.40 μm×4.40 μm，最大FPS為30 FPS，快門速度 0.02 ms～10 s，鏡頭為 Schneider Kreuznach Xeonplan f/2.8-50mm。

DIC方法為非接觸式變形測量方法，為驗證此種方法測量高溫應變的準確性，采用高溫引伸計這種接觸式變形測量系統(tǒng)進行驗證，在1400 ℃實驗環(huán)境下使用 DIC和引伸計同時測量了超高溫陶瓷試件的拉伸應變數(shù)據(jù)，如圖6所示。兩種方式測得1400 ℃實驗環(huán)境下超高溫陶瓷的應變-時間曲線如圖7所示，可以看出，兩種測量手段得到的應變-時間曲線吻合良好，方差為，足以表明DIC方法測量高溫應變的準確性。

圖6 引伸計測量超高溫陶瓷應變Fig.6 Strain measurement of UHTC by Extensometers

圖7 1400 ℃時兩種方法測得的應變數(shù)據(jù)Fig.7 Strain measured by two methods at 1400 ℃

2.2 石墨楊氏模量測試

選取的石墨材料為等靜壓制備，密度為 1.8 g/cm3。標距區(qū)尺寸為30 mm×6 mm×5 mm，按照上述方法在試件表面制作高溫散斑。

為得到石墨在不同溫度下的楊氏模量，分別進行了1200，1500，1700，1800，1900 ℃的拉伸實驗。利用DIC方法分別得到的上述5種溫度時試件表面的豎向應變場如圖 8所示。石墨材料是一種脆性材料，在其斷裂前變形都很小，采集到的變形場圖像信息會受到噪聲干擾，另外高溫下散斑也會有一定的剝離現(xiàn)象，都可能導致應變云圖的起伏甚至缺失。

圖8 不同溫度下試件表面豎向應變場Fig.8 The surface vertical strain field at different temperatures

由高溫力學試驗機和 DIC方法分別輸出試件標距段的載荷和應變數(shù)據(jù)，然后根據(jù)時間軸的對應關(guān)系繪制出應力-應變曲線。不同實驗溫度下的石墨試件的應力-應變曲線如圖 9所示，經(jīng)過線性擬合后得到石墨在不同實驗溫度下的楊氏模量。使用該方法測得石墨在不同溫度下的楊氏模量如圖 10所示，可見石墨的楊氏模量隨溫度的升高呈線性增長趨勢。

圖9 不同溫度下石墨材料的應力-應變曲線Fig.9 Stress-strain curves of graphite at different temperatures

圖10 不同溫度下石墨材料楊氏模量Fig.10 Elastic modulus of graphite at different temperatures

2.3 誤差分析和改進措施

由于該實驗采用VIC-2D商業(yè)軟件進行分析，故不考慮算法本身的誤差，僅考慮實驗過程中設(shè)備、環(huán)境和操作等因素引起的測量誤差，其可能的影響因素為：高溫散斑的質(zhì)量直接影響軟件的運算結(jié)果；溫度升高導致試件本身的熱膨脹引起的變形測量誤差；試件升溫產(chǎn)生的紅光對圖像采集有一定影響；實驗艙內(nèi)空氣在高溫下導致試件表面氧化，對圖像采集有影響；各種環(huán)境噪聲對CCD相機的數(shù)據(jù)采集引起隨機誤差。

針對上述影響因素，為降低測量誤差，對應提出以下改進措施：制作分布均勻、對比度高、顆粒大小適中的高質(zhì)量散斑；在試件升溫至目標溫度后保持一段時間，待熱膨脹變形穩(wěn)定后再開始正式加載，消除熱膨脹引起的誤差；在鏡頭前加裝藍色帶通濾鏡，消除高溫產(chǎn)生的紅外干擾，同時為彌補光照不足，增加兩臺藍光LED燈進行補光，調(diào)整鏡頭光圈大小，使其既有充足的進光量，又不會曝光過度；用機械泵抽取實驗艙內(nèi)空氣使壓強降至5 Pa以下，提高實驗艙真空度，降低空氣的影響；選取較小的加載速率，將相機三腳架的各個節(jié)點固定牢固，降低試驗機和其他設(shè)備的振動對CCD的不利影響。

3 結(jié)語

1）將DIC方法與通電加熱技術(shù)結(jié)合搭建了一套高溫測試系統(tǒng)，可實現(xiàn)高溫下導電材料楊氏模量的測量。

2）對比DIC方法和高溫引伸計測得的應變數(shù)據(jù)，說明了文中所提出的測量方法有效、準確。

3）使用文中方法分別測量了等靜壓石墨不同溫度下的應變場，繪制了相應的應力-應變曲線，進一步計算得到了石墨在 1200～1900 ℃溫度范圍內(nèi)的楊氏模量隨溫度升高的線性增長規(guī)律。

4）由于通電加熱的特性，該系統(tǒng)還可適用于其他導電材料如C/C復合材料在1900 ℃以下的高溫力學性能測試，為飛行器熱防護系統(tǒng)及高溫熱結(jié)構(gòu)設(shè)計提供了數(shù)據(jù)支持。

[1] 王偉. 數(shù)字圖像相關(guān)方法在熱結(jié)構(gòu)材料高溫變形測試中的應用[D]. 哈爾濱: 哈爾濱工業(yè)大學, 2014. WANG Wei. The Application of Digital Image Correlation Method on High Temperature Deformation Test of Thermal Structure Material[D]. Harbin: Harbin Institute of Technology, 2014.

[2] 潘兵, 吳大方, 高鎮(zhèn)同, 等. 1200 ℃高溫熱環(huán)境下全場變形的非接觸光學測量方法研究[J]. 強度與環(huán)境, 2010, 38(1): 52—59. PAN Bing, WU Da-fang, GAO Zhen-tong, et al. Study of Non-contact Optical Metrology for Full-field Deformation Measurement at 1200 ℃[J]. Structure & Environment Engineering, 2010, 38(1): 52—59.

[3] 席濤, 熊宸, 孔繁羽, 等. 基于時序變形預測的數(shù)字圖像相關(guān)加速方法[J]. 實驗力學, 2014, 29(6): 711—718. XI Tao, XIONG Chen, KONG Fan-yu, et al. Accelerating Digital Image Correlation Method Based on the Timing of Deformation Prediction[J]. Journal of Experimental Mechanics, 2014, 29(6): 711—718.

[4] 唐翃翔. 基于數(shù)字圖像的高溫變形場與溫度場的測量[D]. 長春: 吉林農(nóng)業(yè)大學, 2014. TANG Hong-xiang. The Measurement of Deformation Field and Temperature Field under High-temperature Based on Digital Image Process[D]. Changchun: Jilin Agricultural University, 2014.

[5] 周曉峰. 基于散斑數(shù)字圖像相關(guān)的平面全場應變測量方法及應用[D]. 廣州: 華南理工大學, 2012. ZHOU Xiao-feng. Method and Application of In-plane Strain Measurement Based on Speckle Digital Image Correlation[D]. Guangzhou: South China University of Technology, 2012.

[6] MALLON S, KOOHBOR B, KIDANE A, et al. Fracture Behavior of Prestressed Composites Subjected to Shock Loading: A DIC-based Study[J]. Experimental Mechanics, 2015, 55(1): 211—225.

[7] HA N S, VANG H M, GOO N S. Modal Analysis Using Digital Image Correlation Technique: An Application to Artificial Wing Mimicking Beetle’s Hind Wing[J]. Experimental Mechanics, 2015, 55(5): 989—998.

[8] BLABER J, ADAIR B S, ANTONIOU A. A Methodology for High Resolution Digital Image Correlation in High Temperature Experiments[J]. Review of Scientific Instruments, 2015, 86(3): 035111(6pp).

[9] HOLZWEISSING M J, KANAGARAJAH P, MAIER H J. Digital Image Correlation at High Temperatures for Fatigue and Phase Transformation Studies[J]. Journal of Strain Analysis, 2014, 49(4): 204—211.

[10] LYONS J S, LIU J, SUTTON M A. High-temperature Deformation Measurements Using Digital-image Correlation[J]. Experimental Mechanics, 1996, 36(1): 64—70.

[11] PAN Bing, WU Da-fang, WANG Zhao-yang, et al. High-temperature Digital Image Correlation Method for Full-field Deformation Measurement at 1200 ℃[J]. Measurement Science and Technology, 2011, 22(1): 015701(11pp).

[12] GUO Xiang, LIANG Jin, TANG Zheng-zong, et al. High-temperature Digital Image Correlation Method for Full-field Deformation Measurement Captured with Filters at 2600 ℃ Using Spraying to Form Speckle Patterns[J]. Optical Engineering, 2014,53(6): 063101(12pp).

[13] GRANT B M B, STONE H J, WITHERS P J, et al. High-temperature Strain Field Measurement Using Digital Image Correlation[J]. Journal of Strain Analysis, 2009, 44(4): 263—271.

[14] PAN Bing, WU Da-fang, GAO Jian-xin. High- temperature Strain Measurement Using Active Imaging Digital Image Correlation and Infrared Radiation Heating[J]. Journal of Strain Analysis, 2014, 49(4): 224—232.

[15] SU Yun-quan, YAO Xue-feng, WANG Shen, et al. Improvement on Measurement Accuracy of High- temperature DIC by Grayscale-average Technique[J]. Optics and Lasers in Engineering, 2015,75: 10—16.

[16] CHEN Guang, DING Ke-qin, FENG Qi-bo. Research on High-accuracy Two-dimensional Digital Image Correlation Hardware Measurement Systems Used in the Engineering Practice[J]. Optical Metrology and Inspection for Industrial Applications III, 2014, 9726: 92760L(6pp).

[17] BERKE B R, LAMBROS J. Ultraviolet Digital Image Correlation (UV-DIC) for High Temperature Applications[J]. Review of Scientific Instruments, 2014, 85(4): 045121(9pp).

[18] 潘兵, 吳大方, 夏勇.數(shù)字圖像相關(guān)方法中散斑圖的質(zhì)量評價研究[J]. 實驗力學, 2010, 25(2): 120—129. PAN Bing, WU Da-fang, XIA Yong. Evaluation of the Quality of Speckle Patterns with Digital Image Correlation Method[J]. Journal of Experienmental Mechanics, 2010, 25(2): 120—129.

[19] 潘兵，吳大方，高鎮(zhèn)同. 基于數(shù)字圖像相關(guān)方法的非接觸高溫熱變形測量系統(tǒng)[J]. 航空學報, 2010, 31(10): 1960—1967. PAN Bing, WU Da-fang, GAO Zhen-tong. A Non-contact High-temperature Deformation Measuring System Based on Digital Image Correlation Technique[J]. Acta Aeronautica Et Astronautica Sinica, 2010, 31(10): 1960—1967.

[20] 白曉虹. 數(shù)字圖像相關(guān)（DIC）測量方法在材料變形研究中的應用[D]. 沈陽: 東北大學, 2011. BAI Xiao-hong. Application of Digital Image Correlation Method in Study of Material Deformation[D]. Shenyang: Northeastern University, 2011.

Measurement of High-temperature Young's Modulus on Graphite Based on Digital Image Correlation

FENG Wei-qiang,WANG Wei,MENG Song-he,XU Cheng-hai,XIE Wei-hua
(Center for Composite Materials, Harbin Institute of Technology, Harbin 150080, China)

ObjectiveTo obtain the Young's modulus of graphite in high-temperature environment, and provide reliable technical guarantee for the design of aircraft thermal protection systems and high-temperature structure.MethodsCombining the digital image correlation (DIC) method and energized resistance heating technology, this paper established a test system for measuring mechanical properties in high-temperature environment. Using this system, this paper obtained the surface strain field and stress-strain curve of graphite at different temperatures, and calculated Young's modulus of graphite in high-temperature environment. The tensile strain data of Ultra High Temperature Ceramic specimens were measured using DIC and high tempera-ture extensometer method simultaneously at 1400 ℃ in the experimental environment, and comparison and verification were conducted.ResultsThe strain-time curves of UHTC measured by DIC and high-temperature extensometer method at 1400 ℃in the experimental environment were in good agreement, with a variance of 1.3×10-7. In high-temperature environments with temperatures of 1200 ～ 1900 ℃, the Young's modulus of graphite increased linearly with the increase of temperature.ConclusionsThe method based on DIC was accurate and effective, and this method could easily and quickly achieve the measurement of the Young's modulus of graphite material in high-temperature environment. Moreover, the method could also be applied to measure the Young's modulus of other conductive materials.

graphite; high-temperature strain field; digital image correlation method; Young's modulus; stress-strain curve

MENG Song-he (1969—), Male, from Inner Mongolia, Professor, Doctoral tutor, Research focus: testing and characterization of ultra-high temperature thermal protect materials.

10.7643/ issn.1672-9242.2016.03.003

TJ06

A

1672-9242(2016)03-0018-07

2016-01-25；

2016-02-22

Received：2016-01-25；Revised：2016-02-22

國家自然科學基金面上項目（11472092）；國家重點基礎(chǔ)研究發(fā)展計劃（2015CB655200）

Foundation：Supported by National Natural Science Foundation of China（11472092）and National Key Basic Research Development Program（2015CB655200）

封偉強（1990—），男，河北人，碩士研究生，主要研究方向為復合材料力學性能測試。

Biography：FENG Wei-qiang (1990—), Male, from Hebei, Master graduate student, Research focus: mechanical properties of composite materials.

孟松鶴（1969—），男，內(nèi)蒙古人，教授，博導，主要研究方向為超高溫防熱材料的測試與表征。