王璐，易亞楠，劉廣彥，馬沁巍，馬少鵬

（1. 北京理工大學(xué) 宇航學(xué)院，北京 100081；2. 上海交通大學(xué) 船舶海洋與建筑工程學(xué)院，上海 200240）

高溫氣冷堆是第4 代先進核電技術(shù)的代表，具有安全性和發(fā)電效率高等優(yōu)點[1]. 核石墨作為慢化劑、結(jié)構(gòu)材料和反射層廣泛應(yīng)用于高溫氣冷堆的各種基礎(chǔ)建設(shè)中[2?3]，因此對核石墨力學(xué)性能的研究具有重要意義. 對于材料構(gòu)件受載變形的測量分析，目前普遍采用應(yīng)變片等接觸式方法進行測量，然后基于測量的離散點信息再進行統(tǒng)計分析. 過少的測量點往往會產(chǎn)生較大的測量噪聲，而且應(yīng)變片實施起來工序繁瑣，不利于工程大批量測試要求[4]. 為了更好地對核石墨材料變形狀態(tài)進行分析，研究非接觸、高性能、便捷的測量方法至關(guān)重要.

數(shù)字圖像相關(guān)（digital image correlation, DIC）[5]方法又稱為數(shù)字散斑相關(guān)方法(digital speckle correlation method，DSCM)，是一種基于圖像分析的非接觸變形場測量方法. DIC 利用數(shù)字相機記錄試件表面的散斑圖像，然后對散斑圖像進行相關(guān)分析獲得試件表面的變形場. DIC 方法以其對測試環(huán)境要求低、實施簡單、測量范圍易調(diào)節(jié)等優(yōu)勢[6]而被廣泛接受并應(yīng)用. 劉劉等[7]基于數(shù)字圖像相關(guān)技術(shù)實現(xiàn)了多孔氣凝膠基復(fù)合材料彈性力學(xué)常數(shù)的識別. 曾祥福等[8]基于三維數(shù)字圖像相關(guān)技術(shù)對鋁板及石墨圓柱試件動載下的材料力學(xué)行為進行了研究. DIC 方法對于簡單表面的簡單變形（如伸縮變形）的測量精度基本滿足工程測量的需求. 然而，在很多實際工況下，局部變形場（應(yīng)變高度集中區(qū)，如剪切帶[9]、試件缺口附近[10?11]和裂紋尖端[12]）的位移測量誤差可達到百分之幾十像素甚至整像素[13]. 這導(dǎo)致使用者很難判斷測量結(jié)果的可靠性. 核反應(yīng)堆中的構(gòu)件往往結(jié)構(gòu)復(fù)雜，存在大量溝槽結(jié)構(gòu)，局部具有應(yīng)力集中效應(yīng)，因此常規(guī)子區(qū)DIC 方法難以滿足測量要求. 大多復(fù)雜形狀試件或結(jié)構(gòu)的力學(xué)性能以及破壞過程與其邊界的應(yīng)變狀態(tài)有重要的關(guān)聯(lián). 如圓環(huán)受壓結(jié)構(gòu)，它的破壞往往是從圓環(huán)邊界開始的[14]. 因此復(fù)雜形狀試件或結(jié)構(gòu)邊界的變形成為很多實驗的觀測重點. 在多數(shù)實際工況下，局部變形場無法通過應(yīng)變片測量，因為應(yīng)變片尺寸內(nèi)高梯度的應(yīng)變會被平均. 此時，采用DIC 方法雖然能夠測量這些高應(yīng)變梯度區(qū)的局部變形場，但由于子區(qū)DIC 僅能計算獲得子區(qū)中心點的變形信息[15]，因此無法直接計算獲得邊界處的變形信息. 網(wǎng)格DIC[16?17]可以實現(xiàn)包含邊界的全場測量，但由于其計算過程復(fù)雜、運行效率低下，難以滿足工程上大規(guī)模測試的要求.

為了準確地獲得核石墨構(gòu)件表面變形場，本文發(fā)展了一種考慮力學(xué)約束的DIC 方法（mechanicsconstrained，DIC），以下簡稱為MC-DIC 方法. MC-DIC方法保留了子區(qū)DIC 計算便捷的優(yōu)勢，針對其無法實現(xiàn)邊界測量、測量結(jié)果不具有空間連續(xù)性等問題進行了改進. 文章通過模擬實驗驗證了算法的有效性，并將其成功地應(yīng)用于真實核石墨圓環(huán)對徑壓縮試件全場變形測量.

1 MC-DIC 基本原理

MC-DIC 的計算流程主要分為兩部分：子區(qū)DIC計算與考慮力學(xué)約束的位移場空間重構(gòu). 首先采用子區(qū)DIC 進行相關(guān)計算獲得試件表面位移場和相關(guān)系數(shù)場；其次通過相關(guān)系數(shù)[18]確定高精度位移點（以下簡稱信賴點）；之后采用8 節(jié)點等參單元在試件表面構(gòu)建網(wǎng)格，單元與單元之間通過節(jié)點連接形成力學(xué)約束；然后基于信賴點實現(xiàn)全部單元節(jié)點位移的反演；最后基于得到的單元節(jié)點位移實現(xiàn)整個變形場的重構(gòu)，從而實現(xiàn)試件表面變形場的求解. 流程圖如圖1 所示.

圖1 MC-DIC 計算流程圖Fig. 1 Flow chart of MC-DIC

相關(guān)系數(shù)計算如下所示

式中：f(x,y) 與g(x+u,y+v)分別為變形前后各像素點的灰度值；與為變形前后計算子區(qū)內(nèi)的灰度平均值；u，v為子區(qū)中心點的位移；2N為參考圖像子區(qū)的寬度.

圖1 中黑色方形代表非信賴點，即相關(guān)系數(shù)較低的像素點；灰色圓形代表信賴點，即相關(guān)系數(shù)較高的像素點；黑色圓點代表網(wǎng)格節(jié)點. 在傳統(tǒng)子區(qū)DIC計算過程中，由于各種因素的影響，許多散斑點往往無法實現(xiàn)很好地匹配，尤其是試件邊界附近，使得這些計算點的相關(guān)系數(shù)較低. 考慮到試件表面的變形場必然服從力學(xué)約束，因此通過信賴點再加上力學(xué)先驗約束即可求解出高質(zhì)量的全場變形. 相關(guān)系數(shù)（式（1））表征了該點位移計算的置信度，相關(guān)系數(shù)越大（越接近1），表明該位移值計算得越準確. 為了獲得高精度的表面變形場，本文選取相關(guān)系數(shù)C>0.95的計算點作為信賴點. 考慮到已有大量關(guān)于子區(qū)DIC 原理介紹的文獻，此處僅重點介紹位移場空間重構(gòu)的基本原理.

MC-DIC 空間位移場重構(gòu)借鑒了有限單元法的思想，將試件表面劃分為若干個單元，單元與單元之間通過節(jié)點連接，保證了計算區(qū)域位移的連續(xù)性，符合試件表面真實位移的空間分布[19]. 由于節(jié)點位于單元的邊界，因此將單元節(jié)點布置于試件邊界處時即可計算試件邊界處的變形. MC-DIC 以每個單元節(jié)點的位移作為參數(shù)空間，將單元內(nèi)的信賴點通過形函數(shù)與各自單元的8 個節(jié)點建立位移約束關(guān)系，基于已知信賴點位移通過優(yōu)化的方式反演出所有節(jié)點的位移，之后再基于節(jié)點位移重構(gòu)出單元中所有像素點的位移，如圖2. 考慮到實際構(gòu)件中單元形狀的不規(guī)則性，此處采用了有限單元法中等參元的思想，將每個單元從全局坐標系xoy中映射到局部坐標系ξoη中，2 個坐標系的數(shù)學(xué)關(guān)系如下

圖2 8 節(jié)點等參單元中全局坐標系和局部坐標系的關(guān)系Fig. 2 Relationship of global and local coordinates in the 8-node quadrilateral elements

式中：ξ和η為局部坐標系的橫縱坐標軸；x和y為全局坐標系坐標軸，Ni(ξ,η)為單元內(nèi)第i個節(jié)點的形函數(shù)，與有限元法的形函數(shù)相同.

在全局坐標系xoy中，像素點(x,y)的位移可以用局部坐標系ξoη中相對應(yīng)的點(ξ,η)來表示，即

式中：u和v分別為節(jié)點沿橫軸和縱軸方向的位移；Ω為計算區(qū)域；m為單元編號；i為單元內(nèi)節(jié)點編號.

目標函數(shù)可表達為

式中：t和h分別為子區(qū)DIC 計算得到的信賴點和與之對應(yīng)的重構(gòu)位移點；M為目標區(qū)域信賴點的數(shù)量.因此，位移場的求解最終可以轉(zhuǎn)化為求解Jloss全局最小值問題，此處可以采用NR 算法[20]、LM 算法[21]和遺傳算法[22]等，完成求解后可一次性獲得全場變形.

2 模擬實驗

為驗證上述MC-DIC 方法的有效性，首先進行了開孔方板壓縮模擬實驗，模型示意圖及網(wǎng)格劃分如圖3 所示，線載荷密度為1.4×107N /pixel，楊氏模量E=1010Pa，泊松比v=0.33. 實驗以式（6）作為目標函數(shù)，采用LM算法進行優(yōu)化計算，當(dāng)Jloss的改變量小于10?6或待優(yōu)化參數(shù)變化量小于10?3時候迭代終止.實驗所用的參考圖像來源于真實的散斑圖像，圖像采用8 bit 編碼，灰度介于0～255 之間. 變形圖像通過參考圖像使用立方插值生成，每個整像素點的位移由開孔壓縮理論解（式（7）～（9））生成

圖3 開孔方板模型示意圖及MC-DIC 網(wǎng)格布置圖Fig. 3 The geometry of open-hole specimen and meshing scheme in MC-DIC

式中：

分別采用子區(qū)DIC 和MC-DIC 求解位移場，子區(qū)DIC 的相關(guān)窗設(shè)置為9 pixels，計算結(jié)果圖4 所示.圖4（a）為采用式（7）～（9）生成的理論解；圖4（b）為采用子區(qū)DIC 方法計算得到的位移場；圖4（c）為MC-DIC 計算得到的位移場，其選用相關(guān)系數(shù)0.95 的點作為信賴點. 圖4（d）為MC-DIC 計算得到的位移場，但其選用相關(guān)系數(shù)C>0 的點作為信賴點，即所有DIC 計算點都參與空間重構(gòu). 從圖中可以看出，圖4（b）所示位移場不光滑且不具有連續(xù)性，尤其是在模型邊界處變形更為明顯. 圖4（c）與圖4（a）最為接近，位移場光滑并具有連續(xù)性. 圖4（d）計算得到的位移場也光滑并具有連續(xù)性，但由于采用所有計算點進行重構(gòu)，大量低匹配度的位移點嚴重影響了重構(gòu)效果.為了更加直觀地展示出不同方法的計算效果，圖5 給出了圖3 黑線處計算得到的像素點位移（從左至右）. 從圖中可以看出DIC、MC-DIC(C>0.95)與理論解在非邊界區(qū)域計算結(jié)果較為吻合，考慮信賴點的MC-DIC(C>0.95)在邊界區(qū)域的計算結(jié)果與理論解更為接近. 這表明MC-DIC 方法在考慮信賴點的情況下可以很好地實現(xiàn)材料全場變形測量.

圖4 不同方法計算得到的位移場Fig. 4 Displacement fields obtained by different methods

圖5 不同方法計算得到的圖3 黑線處位移Fig. 5 Displacement obtained by different methods (blackline in Fig.3)

3 實 驗

為了獲得核石墨構(gòu)件表面高精度變形場，本節(jié)采用MC-DIC 方法對其表面變形場進行測量. 圖6 與表1 給出了核石墨圓環(huán)構(gòu)件的尺寸和材料參數(shù)，實驗采用WDW-50E 試驗機，利用位移加載方式，加載速度為0.05 mm/min，CCD 相機的分辨率為1 628×1 236. 試驗采用50 mm 長焦鏡頭來捕捉目標區(qū)域的散斑圖像，采集幀率為1 fps(圖7). 試件的表面噴涂不規(guī)則的散斑用于追蹤像素點的位移量. 獲得的圖片采用VIC-2D 進行處理，子區(qū)大小為47 pixels，計算步長為10.

圖6 核石墨圓環(huán)對徑壓縮實驗示意圖及MC-DIC 網(wǎng)格布置圖Fig. 6 Diametrical compression test diagram of the nuclear graphite ring and meshing scheme in MC-DIC

圖7 核石墨圓環(huán)對徑壓縮試驗布置Fig. 7 Experimental setup of nuclear graphite diametrical compression test

表1 核石墨圓環(huán)尺寸及材料參數(shù)Tab. 1 Ring specimen dimensions and mechanical properties of nuclear graphite

圖8 展示了子區(qū)DIC 與MC-DIC 計算得到核石墨圓環(huán)表面位移場. 從圖8 中可以看出MC-DIC 計算結(jié)果（圖8（c）8（d））更加平滑，其包含了圓環(huán)邊界的變形信息. 子區(qū)DIC 無法獲得邊界變形信息，由于計算位移點之間并沒有力學(xué)約束，因此計算得到的位移場并不連續(xù). 由于子區(qū)DIC 無法計算出核石墨邊界處的變形信息，因此無法通過子區(qū)DIC 計算結(jié)果準確獲得材料的力學(xué)參數(shù)，例如拉伸破壞應(yīng)變和抗拉強度等.

圖8 不同方法計算得到的核石墨圓環(huán)位移場Fig. 8 Displacement fields of nuclear graphite ring calculated by different methods

為了更加直觀地感受2 種方法計算的效果，圖9給出了圖6 黑色標記處的三維形貌圖. 從圖中可以直觀地看出，相比于子區(qū)DIC 計算結(jié)果凹凸不平，MC-DIC 計算結(jié)果更為平滑，更具有連續(xù)性.

圖9 不同方法計算得到的圖6 黑色標記處三維形貌圖(水平位移場U)Fig. 9 3D Morphology map (horizontal displacement field) obtained by different methods

4 結(jié) 論

本文在子區(qū)DIC 的基礎(chǔ)上發(fā)展了一種考慮力學(xué)約束的變形測量方法，通過模擬實驗對算法的準確性進行了驗證，并成功地應(yīng)用到真實核石墨試件表面變形場的測量中，得到如下主要結(jié)論：

①與傳統(tǒng)子區(qū)DIC 相比，考慮力學(xué)約束的MCDIC 計算結(jié)果具有空間連續(xù)性，符合材料真實變形行為，可以更好地實現(xiàn)材料表面高精度變形場測量.

②基于MC-DIC 方法可以實現(xiàn)材料邊界變形測量.

③MC-DIC 方法可以推廣到其他材料變形的測量中，為材料的大批量試驗及保證測試數(shù)據(jù)可靠性提供了方便和保障.