數(shù)字圖像相關(guān)方法中數(shù)字散斑場的制作與應(yīng)用研究

錢 偉， 蔣 明

(1．蘇州科技大學(xué) 數(shù)理學(xué)院，江蘇 蘇州 215009;2．蘇州科技大學(xué) 土木學(xué)院，江蘇 蘇州 215009)

1 引 言

數(shù)字圖像相關(guān)(Digital image correlation, DIC)測量技術(shù)是通過試件表面的圖案或者紋理變化獲得其位移和應(yīng)變等信息。DIC實驗需要在試樣的表面生成具有隨機灰度分布的散斑圖案，即散斑場，在不影響試件本身力學(xué)行為的基礎(chǔ)上散斑場能夠隨著試件表面的變形而產(chǎn)生相應(yīng)的變化。散斑場通常采用人工噴漆方法進行制備，由于不同的人噴涂的手法和周圍環(huán)境的不同，從而導(dǎo)致噴漆散斑質(zhì)量參差不齊。

在DIC高溫環(huán)境的測量中，許多學(xué)者對相關(guān)問題進行了系統(tǒng)的研究和探索。時金娜[1]等人在混凝土表面噴涂油漆，利用DIC技術(shù)，對經(jīng)歷20，200，400，600 ℃后的混凝土在單軸受壓過程中的全場位移進行采集，分析了混凝土在不同溫度下、不同受力階段應(yīng)變云圖的變化特征和變形參數(shù)的變化規(guī)律；吳大方[2]等人將白色的Al2O3陶瓷板表面制備成黑色，之后將白色的Al2O3粉末加入高溫陶瓷粘接劑中攪拌均勻，通過濺射方式制備出散斑場，實現(xiàn)了1 200 ℃高溫氧化環(huán)境下部件受熱前表面的應(yīng)變測量；胡悅[3]等人利用不同顏色的高溫漆制斑，在不同的溫度節(jié)點外加不同光源及相應(yīng)的濾波片，采集并觀察圖像是否具有良好的對比度，實驗結(jié)果表明：采用紫外光照明可以實現(xiàn)DIC在高溫環(huán)境下的測量，同時利用黑色或者藍色散斑直接噴涂在試件上有著最佳的對比度。這些采用高溫粘接劑和高溫油漆等方法制作的散斑場，不僅成本高，散斑場質(zhì)量不穩(wěn)定，且耐受溫度有限；張平等人[4]利用激光打標(biāo)的方法在PBT材料上標(biāo)刻了不同的圖案，并且在750 ℃的高溫下標(biāo)刻圖案的質(zhì)量均表現(xiàn)為優(yōu)，說明了激光打標(biāo)圖案能夠適應(yīng)高溫環(huán)境。

運用軟件制作的數(shù)字散斑場(Digital Speckle Patterns, DSPs)能夠有效解決制斑差異大、質(zhì)量不穩(wěn)定等問題。采用光纖激光打標(biāo)技術(shù)進行數(shù)字散斑場的轉(zhuǎn)移。將數(shù)字散斑場運用于高溫焊接實驗，在焊接的瞬間，試件表面溫度迅速升高，散斑隨著試件一起變形，且散斑場質(zhì)量穩(wěn)定，沒有出現(xiàn)高溫融化現(xiàn)象。

2 數(shù)字散斑場的制作及轉(zhuǎn)移

2.1 數(shù)字散斑場的制作

DIC通過識別附著在物體表面隨機分布的散斑來獲取物體的變形信息，達到測量全場變形的目的。因此，散斑作為變形信息的載體，其質(zhì)量直接影響到DIC的變形測量精度和效率。

一般來說，形成散斑場的散斑點形狀可以是任意的，如三角形、正方形、多邊形和圓，但是考慮到散斑的各向同性，選擇的散斑點的形狀為圓形,采用圓點設(shè)計散斑場主要考慮3個要素：占空比ρ(一幅散斑場中散斑點的面積與整幅圖像面積的比值)、散斑直徑d和隨機度r[5-10]。本文設(shè)計的散斑場的散斑直徑相同，通過控制散斑點隨機出現(xiàn)的范圍來控制數(shù)字散斑場的隨機度。如圖1，在白色背景上規(guī)則地放置固定直徑的黑色圓點，形成規(guī)則分布的散斑點陣列(圖1(a))；將該散斑點陣列中的第i個點添加一個隨機的擾動，使得該點在指定范圍內(nèi)并且只能在該范圍內(nèi)移動，如圖1(a)中的紅色方框；遍歷所有散斑點，就可以得到隨機分布的數(shù)字散斑場(圖1(b))。

圖1 軟件生成的散斑場Fig.1 Software generated DSPs

數(shù)字散斑場是通過控制斑點數(shù)量、圓心坐標(biāo)和圓的半徑來設(shè)計制作。數(shù)字散斑場由下面4個公式生成：

(1)

(2)

(3)

n=ρA/(0.25·πd2)，

(4)

2.2 數(shù)字散斑場的轉(zhuǎn)移

圖2 光纖激光打標(biāo)機Fig.2 Fiber optic laser marking machine

目前常用的散斑場轉(zhuǎn)移方法(如水轉(zhuǎn)印[11]、絲網(wǎng)印刷等)制作出的散斑場不能滿足超高溫環(huán)境的測量需求。采用光纖激光打標(biāo)機轉(zhuǎn)移散斑場，選擇的光纖激光打標(biāo)機的激光波長為1 064 nm，功率為100 W，標(biāo)刻深度為0.01～1.6 mm，重復(fù)定位精度為±0.002 mm，圖2為實驗中用到的光纖激光打標(biāo)機的裝置圖。激光打標(biāo)機轉(zhuǎn)移散斑場是通過高能量密度的激光燒蝕試件表面，在試件表面形成圖案，并且不改變試件的原有材質(zhì)。為了不改變試件的力學(xué)性質(zhì),在標(biāo)刻時選擇的標(biāo)刻深度為0.01 mm，遠小于試件的厚度2 mm，所以將散斑場標(biāo)刻在試件表面不會改變試件本身的力學(xué)性質(zhì)，耐受溫度與試件一致。圖3為標(biāo)刻完成的散斑場。

圖3 轉(zhuǎn)移的散斑場圖Fig.3 Enlarged image of transferred DSPs

3 數(shù)字散斑場精度分析

3.1 數(shù)字散斑場打印精度分析

激光打標(biāo)技術(shù)轉(zhuǎn)移數(shù)字散斑場就是將制作好的散斑圖導(dǎo)入打標(biāo)軟件中，然后調(diào)整好激光的功率和標(biāo)刻的視場，通過激光燒蝕試件表面，獲得散斑場。激光打標(biāo)技術(shù)理論的最小打標(biāo)點為原點，但是在實際應(yīng)用中，由于溫度和周圍環(huán)境的影響，可能會對散斑的大小和形狀造成影響。激光打標(biāo)的精度可以從散斑的圓度和面積比兩個方面來進行分析。

3.1.1 圓度

圓度體現(xiàn)了散斑點的再現(xiàn)狀態(tài)，也體現(xiàn)了散斑點的保真度，圓度的表示公式如下所示：

(5)

式中，A為散斑點面積，P為散斑點周長。當(dāng)散斑點的圓度為1時，散斑點處于理想狀態(tài)。散斑點圓度越接近于1，其毛細現(xiàn)象、羽化現(xiàn)象、滲透現(xiàn)象就越少，圖像質(zhì)量就越高。首先通過Matlab軟件讀入拍攝的散斑圖(圖3)并進行二值化處理，將灰度值高于閾值大小的變?yōu)?，低于閾值大小的變?yōu)?，求解所有散斑點的周長和面積并求取圓度值，然后求出圓度的平均值。

導(dǎo)入的散斑場中共有1 845個散斑點，軟件識別到1 789個散斑點，識別率為96.96%，通過分析得到制作的散斑圖的平均圓度為0.978 2。

3.1.2 面積比

在實際應(yīng)用中，散斑點的打印面積由激光打標(biāo)機的精度決定，散斑點的指定面積是指一個方形像素在實際打印中的理論面積，它的邊長為分辨率的倒數(shù)。散斑點的理想面積指的是像素矩形外接圓的面積。散斑點的面積比如式(6)所示：

(6)

選取上述散斑圖，統(tǒng)計二值化之后散斑圖的所有白色點的面積，求取白色面積和整幅圖像面積的比值。最后分析得到黑色部分面積/全部面積=0.526 7，面積比=0.5/0.626 7=0.949 3。

通過圓度和面積比的分析，激光打標(biāo)方法轉(zhuǎn)移數(shù)字散斑場的精度高，滿足實驗要求，而且轉(zhuǎn)移方法安全環(huán)保。

3.2 平均灰度梯度法

根據(jù)潘兵[12]等人的理論與實驗研究，平均灰度梯度可以有效評價一副散斑圖的質(zhì)量，而且平均灰度梯度越大，散斑圖的質(zhì)量越高，該參數(shù)的定義如下：

(7)

首先對采集到的散斑圖進行亞像素平移，每一幅散斑圖連續(xù)精確平移20次，相鄰兩幅散斑圖之間的平移量為0.05個像素，最大平移量為1個像素。圖4為不同的散斑圖以及相對應(yīng)的平均灰度梯度曲線，其中圖4(a)選用的是激光打標(biāo)在試件上的散斑圖，圖中的激光打標(biāo)散斑場的散斑直徑為0.3 mm，占空比為0.5，隨機度為0.35。圖4(b)選取的是平時實驗中質(zhì)量較好的一張噴漆散斑圖。表1為不同散斑進行數(shù)值模擬之后的平均灰度梯度值。由表格數(shù)據(jù)和平均灰度曲線圖可以看出，在相同的計算模板下數(shù)字散斑場的質(zhì)量高于噴漆散斑。

(a)激光打標(biāo)散斑圖及相對應(yīng)的平均灰度梯度曲線(a)Laser marking speckle patterns and the corresponding average gray gradient curve

(b)人工散斑圖及相對應(yīng)的平均灰度梯度曲線(b)Speckle patterns drawn by spray painting and the corresponding average gray gradient curve圖4 不同散斑圖以及對應(yīng)的平均灰度梯度曲線Fig.4 Different speckle patterns and corresponding average gray gradient curves

表1 兩幅不同散斑圖的平均灰度梯度Tab.1 Average gray gradient of two different DSPs

4 數(shù)字散斑的可靠性實驗

4.1 拉伸實驗

圖5 試件實際尺寸Fig.5 Actual size of test piece

拉伸實驗選擇的試件為6061T6鋁板，尺寸如圖5所示，設(shè)計板厚2 mm。數(shù)字散斑場散斑物理直徑為0.3 mm，隨機度為0.35，占空比為0.5。實驗采用3種散斑場形式：①試件全部是激光打標(biāo)的數(shù)字散斑場；②試件一半采用激光打標(biāo)數(shù)字散斑場，一半采用人工噴涂散斑場；③試件全部為人工噴涂散斑場，具體散斑場如圖6所示。

(a)激光打標(biāo)散斑圖(a) Laser marking DSPs

(b)一半激光打標(biāo)，一半噴漆散斑圖。(b) Half laser marking, half spray paint speckle patterns.

(c) 噴漆散斑圖(c) Spray paint speckle patterns圖6 試件及表面散斑制備Fig.6 Specimen and surface speckle preparation

圖7為實驗環(huán)境圖，實驗設(shè)備主要包括拉伸實驗機、ARAMIS數(shù)字圖像相關(guān)系統(tǒng)、兩個工業(yè)CCD(2 048 pixels×2 048 pixels)、兩個50 mm鏡頭等。實驗中先將兩個相機進行雙目立體視覺標(biāo)定[14]，然后采用PCIE8620數(shù)據(jù)采集卡同步控制實驗機力傳感器，實現(xiàn)CCD圖像與力值數(shù)據(jù)的同步采集[15]。加載方案參照《GB/T228.1-2010金屬材料拉伸實驗第一部分：室溫實驗方法》，拉伸時控制應(yīng)變速率εe=0.00 025/s。CCD的采樣頻率為1幀/s，每采集一組圖像記為一個stage，記錄的第一組圖像記為stage 0，每個stage對應(yīng)一個力值數(shù)據(jù)[16-17]。記錄的荷載時間曲線見圖8。開始拉伸的同時進行散斑圖的采集，當(dāng)試件被拉斷之后，結(jié)束圖像采集。拉斷之后的斷裂處的局部放大圖像如圖9所示。可以看到，斷裂的過成中激光打標(biāo)的數(shù)字散斑場隨著試件一起變形，由原來的圓形變成了橢圓形，而在噴漆的試件的斷口處出現(xiàn)了噴漆散斑局部剝落以及變模糊的現(xiàn)象。后處理的相關(guān)計算由德國GOM公司的3D DIC商業(yè)軟件ARAMIS完成。

圖7 實驗環(huán)境Fig.7 Lab environment

圖8 荷載時間(stage)曲線Fig.8 Load time (stage) curve

圖9 斷裂處散斑圖Fig.9 Speckle patterns at break

圖10為計算軟件計算的左右相機拍攝的拉斷處的應(yīng)力云圖，通過應(yīng)力云圖不難發(fā)現(xiàn)，由于斷裂處噴漆的掉落導(dǎo)致相機無法識別，應(yīng)力云圖有殘缺，但是激光打標(biāo)的數(shù)字散斑卻表現(xiàn)良好。為了避免由于試件、實驗、操作等差異造成的顯著性差異，對測試的3個試件選取了在局部頸縮階段時中間部位的豎直線，對于第二個試件則選取了兩條平行線；圖11比較了激光打標(biāo)數(shù)字散斑場和傳統(tǒng)噴漆散斑場的差異性。在頸縮階段，由于受到人工散斑噴漆剝落的影響，第二個試件的右半部分和第三個試件的頸縮部位變形數(shù)據(jù)缺失，第一個試件和第二個試件的左半部分激光打標(biāo)數(shù)字散斑場表現(xiàn)較好，由此可以看出激光打標(biāo)數(shù)字散斑場更加適合用作金屬試件的大變形測試。

圖10 應(yīng)力云圖Fig.10 Stress cloud

(a)試件1局部頸縮階段應(yīng)變曲線(a)Strain curve in the local necking stage for sample 1

(b)試件2局部頸縮階段應(yīng)變曲線(b)Strain curve in the local necking stage for sample 2

(c)試件3局部頸縮階段應(yīng)變曲線(c)Strain curve in the local necking stage for sample 3圖11 局部頸縮階段的應(yīng)變曲線Fig.11 Strain curve in the local necking stage

4.2 激光焊接實驗

激光焊接技術(shù)是采用偏光鏡反射激光產(chǎn)生的光束使其集中在聚焦裝置中產(chǎn)生巨大能量的光束，通過脈沖瞬間發(fā)出的激光焦點達到上千攝氏度，在幾毫秒內(nèi)將金屬材料融化。

圖12 焊接試件散斑圖Fig.12 Speckle patterns of welding specimen

將兩塊304不銹鋼板采用激光焊接的技術(shù)進行T形焊接，T型焊接筋板尺寸為80 mm×50 mm×2 mm，腹板尺寸為50 mm×50 mm×2 mm，預(yù)制數(shù)字散斑場范圍為50 mm×60 mm，其表面激光轉(zhuǎn)移散斑場的散斑物理直徑為0.3 mm，隨機度為0.35，占空比為0.5，如圖12所示，激光打標(biāo)之后試件表面依舊非常平滑，沒有出現(xiàn)凹凸不平的現(xiàn)象。焊接前使用丙酮除去表面油污及雜質(zhì)，并且在試件周圍布置了遮光板，排除了光強對實驗結(jié)果的影響。焊接時所選焊接功率為1 500 W，速度30 mm/s[18]，采用氬氣進行保護，保護氣流量為15 L/min, 焊接后自然冷卻至室溫。采用DIC系統(tǒng)進行測量，CCD靶面與試件表面平行，焊接方向從右向左焊，焊接開始時CCD也開始同步采集圖片，焊接過程中相機的采集頻率為40幀/s，焊接過程中共采集67組圖像；冷卻過程的采樣幀率為5幀/s，冷卻8 min共采集2 400組圖像。實驗環(huán)境如圖13所示。

圖13 實驗環(huán)境Fig.13 Lab environment

圖14 焊接之后應(yīng)變云圖Fig.14 Strain cloud after welding

圖15 不同焊接時刻的應(yīng)變曲線Fig.15 Strain curves at different welding moments

在試件正中間垂直于焊縫處畫一條2 cm的計算直線。圖14為焊接完冷卻8 min后CCD采集到的焊接圖片，中間淡藍色部分為應(yīng)變計算區(qū)域。圖15(a)(c)為激光焊接不同階段的應(yīng)變云圖以及計算直線處的應(yīng)變曲線。圖15(a)為剛開始焊接時，此時還沒有焊接到試件的正中間，所以畫線處的應(yīng)力趨于0；圖15(b)為焊接到紅色計算曲線時的應(yīng)變曲線，此時由于焊縫處溫度的上升，局部應(yīng)力增大，遠離焊縫區(qū)域未受影響，焊縫處應(yīng)力出現(xiàn)峰值，為208 MPa；圖15(c) 為冷卻8 min之后的應(yīng)變曲線，此時測出試件中間焊縫處的殘余應(yīng)力為103 MPa。

5 結(jié) 論

本文提出了一種應(yīng)用于高溫環(huán)境的數(shù)字散斑場，并且利用激光打標(biāo)機將其標(biāo)刻到試件的表面；然后通過圓度和面積比分析了轉(zhuǎn)移的數(shù)字散斑場，可以看出激光打標(biāo)轉(zhuǎn)移的數(shù)字散斑場圓度和面積比分別為0.978 2和0.949 3，標(biāo)刻精度高；通過分析數(shù)字散斑場和普通噴漆散斑場的平均灰度梯度值可以看出數(shù)字散斑場的質(zhì)量比較高。

將數(shù)字散斑場運用到拉伸和高溫焊接實驗，實驗結(jié)果表明，激光打標(biāo)機轉(zhuǎn)移的數(shù)字散斑場能夠隨著試件一起變形，并且在2 000 ℃的環(huán)境中沒有出現(xiàn)融化的現(xiàn)象，并且測出實驗中的不銹鋼試件的殘余應(yīng)力為103 MPa。