杜亞志， 王學(xué)濱,2， 董 偉， 馮威武， 侯文騰

(1. 遼寧工程技術(shù)大學(xué) 力學(xué)與工程學(xué)院， 遼寧 阜新 123000； 2. 遼寧工程技術(shù)大學(xué) 計(jì)算力學(xué)研究所， 遼寧 阜新 123000)

1 引 言

目前，CT技術(shù)、立體成像技術(shù)、數(shù)字圖像相關(guān)(digital image correlation, DIC)方法[1～5]等先進(jìn)的觀測手段均被用于測量剪切帶內(nèi)外微結(jié)構(gòu)特征、研究剪切帶的發(fā)展演化等。

文獻(xiàn)[6]針對(duì)DIC方法中圖像的亞像素插值誤差較大的問題，提出了一種基于B樣條插值核的插值濾波器;文獻(xiàn)[7]通過對(duì)散斑圖進(jìn)行平移，比較了梯度法、曲面擬合法及N-R(Newton-Raphson)法的精度;文獻(xiàn)[8]提出了一種散斑圖質(zhì)量評(píng)價(jià)參數(shù)——平均灰度梯度;根據(jù)文獻(xiàn)[8]提出的參數(shù)，文獻(xiàn)[9]提出了一種散斑圖質(zhì)量評(píng)價(jià)參數(shù)——平均一階及二階灰度梯度商;文獻(xiàn)[4]使用整像素精度的DIC方法測量了水平常應(yīng)變剪切帶的寬度;文獻(xiàn)[5]根據(jù)梯度塑性理論制作了含應(yīng)變梯度的水平剪切帶，比較了一階及二階DIC方法的位移場、應(yīng)變場及應(yīng)變梯度場的測量結(jié)果的差異。

為了解基于DIC法的剪切帶的測量精度，本文制作了簡單剪切條件下水平(傾角為0°)和傾斜(傾角為60°)常應(yīng)變剪切帶，并使用DIC法對(duì)剪切帶內(nèi)、外的位移場、應(yīng)變場和寬度進(jìn)行了測量；分析了形函數(shù)、子區(qū)尺寸及測點(diǎn)間距對(duì)測量結(jié)果的影響。

2 剪切帶的制作方法

常見的剪切帶模型為簡單剪切模型，在該模型中，剪切帶內(nèi)的變形是簡單剪切，剪切帶外僅有剛體位移。在水平常應(yīng)變剪切帶中，剪切位移為水平位移，若以剪切帶中心線上任意點(diǎn)為原點(diǎn)建立直角坐標(biāo)系x′O′y′,x′和y′分別為剪切帶的切向和法向，則有剪切帶水平位移計(jì)算公式：

s(y′)=γ0y′

(1)

式中：γ0為上述水平剪切帶的剪應(yīng)變；y′為剪切帶的法向坐標(biāo)，y′∈[-w/2,w/2]，w為剪切帶寬度。

在該模型中，剪切帶的傾角θ為零，而在通常情況下剪切帶是傾斜的。文獻(xiàn)[2]中僅簡單提及了傾斜剪切帶的制作方法，即“首先，將一張實(shí)驗(yàn)圖片旋轉(zhuǎn)60°；然后，在一定位置進(jìn)行常應(yīng)變簡單剪切；最后，再將圖像反向旋轉(zhuǎn)60°，得到包含傾斜剪切帶(傾角為60°)的圖像”。這里采用簡單剪切模型及文獻(xiàn)[2]提出的制作方法制作。具體步驟如下：

(1)定義剪切帶邊界。理論上，剪切位移是剪切帶法向坐標(biāo)的連續(xù)函數(shù)，而在剪切帶制作過程中，僅對(duì)圖像中離散的像素進(jìn)行操作。圖1給出了垂直于剪切帶的一列像素、剪切帶邊界的位置及該列像素位移的示意圖。在圖1中,像素用小正方形表示，像素的中心點(diǎn)用圓點(diǎn)表示，剪切帶邊界用兩條平行線段表示;圖1(a)為剪切帶邊界在兩個(gè)相鄰像素之間；圖1(b)為剪切帶邊界在像素的中心。于是產(chǎn)生了在位移曲線中剪切帶寬度的定義不同，寬度分別為w1和w2。在圖1(b)中，剪切帶內(nèi)靠近邊界的點(diǎn)的位移等于最大剪切位移，據(jù)此制作剪切帶，有利于與有關(guān)的理論結(jié)果作對(duì)比。

(2)對(duì)剪切帶外上下兩盤進(jìn)行剛體平移，平移量為γ0w2。

(3)使用仿射變換方法對(duì)圖像進(jìn)行變形，使用光滑性較好的雙三次樣條函數(shù)對(duì)圖像進(jìn)行插值。

圖1 剪切帶的兩種邊界的定義

圖2給出了源圖像及含剪切帶的圖像。圖2(a)為一幅散斑圖[10]，據(jù)此制作了水平剪切帶(圖2(b))和傾斜剪切帶(圖2(c))，其中，γ0=-0.1745，w=30像素，θ=60°。

圖2 源圖像及含剪切帶的圖像

3 結(jié)果分析

在圖像坐標(biāo)系xOy下，從圖2(a)變形至圖2(b)，根據(jù)式(1)可計(jì)算出最大水平位移為2.62像素，垂直位移、水平線應(yīng)變?chǔ)舩及垂直線應(yīng)變?chǔ)舮的理論值均為零，剪應(yīng)變的理論值γxy=γ0=-0.174 5。 使用基于粒子群優(yōu)化算法及N-R迭代的粗-細(xì)DIC方法[3]分別對(duì)從圖2(a)變形至圖2(b)、從圖2(a)變形至圖2(c)的位移和應(yīng)變進(jìn)行了計(jì)算，計(jì)算區(qū)域如圖2(a)所示(水平剪切帶，使用計(jì)算區(qū)域1；傾斜剪切帶，使用計(jì)算區(qū)域2)。

3.1 水平剪切帶的結(jié)果分析

3.1.1 一階DIC方法的結(jié)果分析

圖3給出了基于一階DIC方法的從圖2(a)變形至圖2(b)的結(jié)果，其中兩條虛線代表剪切帶的邊界。子區(qū)尺寸為31像素×31像素(簡稱為31像素)，測點(diǎn)數(shù)量為23×19，測點(diǎn)間距為5像素。

圖3 一階DIC方法的水平剪切帶的位移場、應(yīng)變場及相關(guān)系數(shù)

從圖3(a)和3(b)中可以發(fā)現(xiàn):圖像上、下兩部分的位移矢量場有較大差別,上部分位移向左，下部分位移向右，且兩部分位移的大小基本相等，這與帶外區(qū)域平動(dòng)相符；在上、下兩部分之間，存在一個(gè)明顯的位移梯度帶，帶內(nèi)的水平位移在法向變化較大，這與剪切帶的位置(圖3(b)中的兩條平行線段之間的區(qū)域)對(duì)應(yīng)。該位移梯度帶的寬度稍大于剪切帶寬度，這表明，在剪切帶邊界附近，位移場存在一定的誤差。

從圖3(c)中可以發(fā)現(xiàn):在遠(yuǎn)離剪切帶及剪切帶中心線附近的區(qū)域，垂直位移均在零附近，這與理論值接近；在剪切帶邊界附近，垂直位移在-0.15和0.2像素之間隨機(jī)變化，表明在剪切帶邊界垂直位移的誤差限較大。

從圖3(d)和3(e)中可以發(fā)現(xiàn):在遠(yuǎn)離剪切帶及剪切帶中心線附近的區(qū)域，兩種線應(yīng)變?cè)诹愀浇?/p>

在剪切帶邊界附近，兩種線應(yīng)變?cè)?0.02～0.02之間隨機(jī)變化。從圖3(f)中發(fā)現(xiàn):在剪切帶內(nèi)，γxy在-0.16～-0.08之間，大于理論值-0.174 5；在剪切帶邊界附近，γxy在-0.08～0之間，表明γxy存在一定的誤差。

從圖3(g)中可以發(fā)現(xiàn):在遠(yuǎn)離剪切帶及剪切帶中心線附近的區(qū)域，相關(guān)系數(shù)均大于0.999，表明位移誤差限較小；在剪切帶邊界附近，相關(guān)系數(shù)在0.993～0.999之間，表明位移誤差限較大。

3.1.2 形函數(shù)對(duì)位移和應(yīng)變的影響

為了減小隨機(jī)誤差的影響，對(duì)圖3中相同y坐標(biāo)的19個(gè)點(diǎn)的位移和應(yīng)變求均值。圖4為不同DIC方法的位移和應(yīng)變。

圖4 不同DIC方法的位移及應(yīng)變

從圖4(a)中可以發(fā)現(xiàn):當(dāng)測點(diǎn)遠(yuǎn)離剪切帶(y>225像素或y<160像素)時(shí)，兩種DIC方法的水平位移均等于理論值；在剪切帶內(nèi)及附近(y=160～225像素)，兩種DIC方法的水平位移均有一定的誤差;在剪切帶中心，位移誤差均基本為零；從剪切帶中心到剪切帶邊界，位移誤差限逐漸增大，且一階DIC方法比二階DIC方法的水平位移誤差限大;從剪切帶邊界到遠(yuǎn)離剪切帶的測點(diǎn)，一階DIC方法的水平位移誤差限逐漸減小，趨于零，最大誤差限約為0.6像素,而二階DIC方法的水平位移最大誤差限約為0.3像素。整體上，在剪切帶內(nèi)及附近，二階DIC方法比一階DIC方法的位移誤差限小。從圖4(a)、4(b)中可以發(fā)現(xiàn):與水平位移誤差限相比，垂直位移誤差限較小;圖4(b)中一階DIC方法及二階DIC方法的垂直位移的范圍分別為-0.05～0.21像素和-0.05～0.05像素,這也表明，在剪切帶內(nèi)及附近，二階DIC方法比一階DIC方法的位移誤差限小。

由圖4(c)可以發(fā)現(xiàn):在剪切帶內(nèi)，兩種DIC方法的剪應(yīng)變均是變化的，在剪切帶中心剪應(yīng)變最小,一階DIC方法的剪應(yīng)變均大于理論值，而二階DIC方法的剪應(yīng)變小于理論值;在剪切帶邊界，兩種DIC方法的剪應(yīng)變基本相等，約為理論值的一半;從剪切帶邊界至帶外約20像素范圍內(nèi)，一階DIC方法的剪應(yīng)變逐漸增加，均小于零(理論值),而二階DIC方法的剪應(yīng)變較為復(fù)雜。

3.1.3 子區(qū)尺寸對(duì)位移和應(yīng)變的影響

圖5給出了從圖2(a)變形至圖2(b)的剪切帶的位移和應(yīng)變的分布規(guī)律，子區(qū)尺寸分別為21、41和61像素，其余計(jì)算參數(shù)同圖3。與圖4類似，已對(duì)相同y坐標(biāo)的19個(gè)點(diǎn)的位移和應(yīng)變求平均。

從圖5(a)、5(b)中可以發(fā)現(xiàn):隨著子區(qū)尺寸的增加，剪切帶內(nèi)及附近區(qū)域(y=150～230像素)的絕大部分測點(diǎn)的位移誤差限逐漸增加;當(dāng)子區(qū)尺寸不同時(shí)，位移及位移誤差均中心對(duì)稱，對(duì)稱中心為剪切帶中心(y=192像素)。下面以剪切帶中心上部(y<192像素)為例進(jìn)行分析。剪切帶中心附近測點(diǎn)的位移誤差基本為零;當(dāng)測點(diǎn)逐漸遠(yuǎn)離剪切帶中心時(shí)，位移誤差先增加，后減小，最后趨于零。在剪切帶邊界(y=177像素)，位移誤差達(dá)到峰值；在位移誤差峰值附近(y=147～192像素)，位移誤差基本對(duì)稱分布。對(duì)這種現(xiàn)象進(jìn)行解釋如下：

圖5 子區(qū)尺寸對(duì)位移和剪應(yīng)變的影響

對(duì)于子區(qū)中心點(diǎn)在剪切帶邊界附近的情況，無論其位于剪切帶內(nèi)或帶外，一階DIC方法中的變形均勻性假設(shè)不適用，因此位移誤差限較大。當(dāng)子區(qū)中心點(diǎn)在剪切帶外時(shí)，子區(qū)覆蓋的一部分剪切帶區(qū)域是干擾區(qū)域;同樣，當(dāng)子區(qū)中心點(diǎn)在剪切帶內(nèi)時(shí)，子區(qū)覆蓋的一部分剪切帶外區(qū)域是干擾區(qū)域。當(dāng)上述子區(qū)的中心點(diǎn)與剪切帶邊界的距離相等時(shí)，子區(qū)中干擾區(qū)域尺寸占子區(qū)尺寸的比例相等，此時(shí)，位移誤差相等(圖5(b))。當(dāng)子區(qū)中心點(diǎn)在剪切帶邊界處時(shí)，該比例達(dá)到最大值0.5，此時(shí)位移誤差達(dá)到最大值或最小值。因此可以推測，位移誤差與干擾區(qū)域的尺寸占子區(qū)尺寸的比例有關(guān)。對(duì)于剪切帶中心附近的子區(qū)，當(dāng)其僅包含帶內(nèi)區(qū)域時(shí)，滿足變形均勻性假設(shè)，因此位移誤差限較??；隨著子區(qū)尺寸的增加，當(dāng)子區(qū)超出剪切帶時(shí)，子區(qū)覆蓋的兩塊剪切帶外區(qū)域是干擾區(qū)域。剪切帶中心附近測點(diǎn)的位移誤差限較小，這可能與干擾區(qū)域?qū)τ?jì)算結(jié)果的影響較小有關(guān)。

從圖5(d)、5(e)可以發(fā)現(xiàn):對(duì)于二階DIC方法的結(jié)果，隨著子區(qū)尺寸的增加，剪切帶內(nèi)及附近(y=147～237像素)測點(diǎn)的位移誤差限基本上逐漸增大。從圖5(b)、5(e)可以發(fā)現(xiàn):當(dāng)子區(qū)尺寸相同時(shí)，對(duì)于剪切帶內(nèi)及其附近測點(diǎn)的位移，二階DIC方法的位移誤差限通常比一階DIC方法的小，這是由于二階DIC方法中，子區(qū)的變形可以是不均勻的。

從圖5(c)、(f)可以發(fā)現(xiàn):在剪切帶中心，隨著子區(qū)尺寸的增加，一階及二階DIC方法的剪應(yīng)變的誤差基本上逐漸增加；在子區(qū)尺寸相同時(shí)，二階DIC方法的剪應(yīng)變誤差比一階DIC方法的小。

3.1.4 測點(diǎn)間距對(duì)應(yīng)變的影響

圖6給出了測點(diǎn)間距分別為2、5及10像素時(shí)一階及二階DIC方法的剪應(yīng)變。隨著測點(diǎn)間距的增加，在遠(yuǎn)離剪切帶的區(qū)域(y=135～150像素、235～250像素)，剪應(yīng)變接近理論值；在剪切帶內(nèi)，剪應(yīng)變-坐標(biāo)曲線由陡峭變得平緩。在剪切帶內(nèi)及附近，隨著測點(diǎn)間距的增加，一階DIC方法和二階DIC方法的剪應(yīng)變有所不同：一階DIC方法的剪應(yīng)變逐漸偏離理論值；在剪切帶邊界外各5像素，二階DIC方法的剪應(yīng)變的誤差隨著測點(diǎn)間距的增加而減小；在剪切帶中心附近(y=187～197像素)，二階DIC方法的剪應(yīng)變的誤差由負(fù)值變?yōu)檎怠?/p>

當(dāng)測點(diǎn)間距分別為2、5及10像素時(shí)，一階DIC方法的剪切帶寬度分別為70、75及90像素，為理論值的2.3～3倍；二階DIC方法的剪切帶寬度均為70像素，為理論值的2.3倍。兩種DIC方法的剪切帶寬度均大于理論值，這可解釋如下：在采用中心差分方法時(shí)，一點(diǎn)的應(yīng)變是根據(jù)周圍多個(gè)測點(diǎn)的位移得到的，剪切帶外測點(diǎn)的應(yīng)變可能被高估，由此導(dǎo)致剪切帶寬度被高估。

圖6 測點(diǎn)間距對(duì)剪應(yīng)變的影響

3.2 傾斜剪切帶的結(jié)果分析

圖7給出了二階DIC方法的從圖2(a)變形至圖2(c)(包含傾斜剪切帶)的位移矢量場、水平位移、垂直位移、εx、εy、γxy及相關(guān)系數(shù)的等值線圖，其中，兩條虛線代表剪切帶的兩個(gè)邊界。樣本子區(qū)尺寸=31像素×31像素，測點(diǎn)數(shù)量為63×19，測點(diǎn)間距為5像素?？梢园l(fā)現(xiàn)：(1)在遠(yuǎn)離剪切帶的位置，測點(diǎn)的位移比較均勻，應(yīng)變基本為零，相關(guān)系數(shù)較大(大于0.999 8)，表明這部分測量結(jié)果比較準(zhǔn)確；(2)在剪切帶內(nèi)及附近，應(yīng)變比較集中，應(yīng)變局部化帶的寬度的計(jì)算值均大于剪切帶寬度的理論值,相關(guān)系數(shù)的范圍為0.998 6～0.999 8，略小于遠(yuǎn)離剪切帶處的值。

圖7 二階DIC方法的傾斜剪切帶的位移場及應(yīng)變場

由圖3和圖7可以發(fā)現(xiàn)，針對(duì)水平和傾斜剪切帶的各種計(jì)算結(jié)果類似，不再贅述。

4 結(jié) 論

為了了解基于DIC法的剪切帶的測量精度，本文制作了水平和傾斜的常應(yīng)變剪切帶，并使用DIC法對(duì)剪切帶內(nèi)、外的位移場和應(yīng)變場進(jìn)行了計(jì)算，得到以下結(jié)論：

(1)在剪切帶邊界附近，一階或二階DIC方法的位移誤差限較大，這是因?yàn)檫@些區(qū)域的變形不均勻，子區(qū)覆蓋一部分剪切帶。子區(qū)尺寸越大，或測點(diǎn)距離剪切帶邊界越近，位移誤差限越大。一般地，對(duì)剪切帶邊界附近測點(diǎn)進(jìn)行計(jì)算時(shí)，二階DIC方法的位移誤差限比一階DIC方法的小，這應(yīng)與在二階DIC方法中允許子區(qū)的變形是不均勻的有關(guān)。

(2)隨著測點(diǎn)間距的增加，一階或二階DIC方法的剪切帶法向的剪應(yīng)變-坐標(biāo)曲線由陡峭變得平緩，剪切帶寬度逐漸增加。

[參考文獻(xiàn)]

[1] Bhandari A R, Inoue J. Experimental study of strain rates effects on strain localization characteristics of soft rocks[J].SoilsandFoundations, 2005, 45(1): 125-140.

[2] Rechenmacher A L, Abedi S, Chupin O,etal. Characterization of mesoscale instabilities in localized granular shear using digital image correlation[J].ActaGeotechnica, 2011, 6(4): 205-217.

[3] 王學(xué)濱, 杜亞志, 潘一山. 基于DIC粗-細(xì)搜索方法的單軸壓縮砂樣的應(yīng)變分布及應(yīng)變梯度的實(shí)驗(yàn)研究[J]. 巖土工程學(xué)報(bào), 2012, 34(11): 2050-2057.

[4] 王學(xué)濱, 杜亞志, 金石. 基于數(shù)字圖像相關(guān)方法的剪切帶寬度的影響因素分析[J].世界科技研究與發(fā)展, 2012, 34(3): 359-363.

[5] 王學(xué)濱, 杜亞志, 潘一山. 考慮一階和二階位移梯度的數(shù)字圖像相關(guān)方法在剪切帶測量中的比較[J]. 工程力學(xué), 2013, 30(7): 282-287.

[6] 任茂棟, 梁晉, 唐正宗, 等. 數(shù)字圖像相關(guān)法中的優(yōu)化插值濾波器[J]. 西安交通大學(xué)學(xué)報(bào), 2014, 48(7): 65-70.

[7] 資新運(yùn), 耿帥, 趙姝帆, 等. 3 種亞像素位移測量算法的比較研究[J]. 計(jì)量學(xué)報(bào), 2015, 36(3): 260-267.

[8] 潘兵, 吳大方, 夏勇. 數(shù)字圖像相關(guān)方法中散斑圖的質(zhì)量評(píng)價(jià)研究[J]. 實(shí)驗(yàn)力學(xué), 2010, 25(2): 120-129.

[9] 杜亞志, 王學(xué)濱, 馮威武, 等. 基于一階及二階灰度梯度的散斑圖質(zhì)量評(píng)價(jià)方法[J]. 光學(xué)技術(shù), 2017,43(2):169-175.

[10] Zhou P, Goodson K E. Subpixel displacement and deformation gradient measurement using digital image/speckle correlation[J].OpticalEngineering, 2001, 40(8): 1613-1620.

[11] Lu H, Cary P D. Deformation measurements by digital image correlation: Implementation of a second-order displacement gradient[J].ExperimentalMechanics, 2000, 40(4): 393-400.