劉綱 顧嘉偉 李孟珠 賀成華

摘 要：針對(duì)傳統(tǒng)損傷識(shí)別需布設(shè)接觸式傳感器而不利于古建筑木結(jié)構(gòu)保護(hù)的問(wèn)題，利用古建筑自身彩繪圖案探索基于數(shù)字圖像相關(guān)原理（DIC）的非接觸式損傷識(shí)別方法。通過(guò)數(shù)碼相機(jī)拍攝古建筑彩繪木梁振動(dòng)時(shí)序圖片，采用整、亞像素位移搜索算法解算木梁上各點(diǎn)的位移響應(yīng)時(shí)程，利用傅里葉變換得到幅頻、相頻曲線，進(jìn)而通過(guò)幅頻幅值和相頻相位角提取木梁的工作模態(tài);利用DIC方法可獲取沿梁長(zhǎng)方向所有像素模態(tài)信息的優(yōu)勢(shì)，引入振型轉(zhuǎn)角模態(tài)、曲率模態(tài)作為損傷特征指標(biāo)進(jìn)行損傷識(shí)別。數(shù)值算例和彩繪簡(jiǎn)支梁實(shí)驗(yàn)結(jié)果表明，采用DIC方法獲得的工作模態(tài)能準(zhǔn)確識(shí)別木梁的損傷部位和損傷程度的相對(duì)大小，為古建筑木結(jié)構(gòu)預(yù)防性保護(hù)提供了新技術(shù)手段。

關(guān)鍵詞：古建筑木結(jié)構(gòu);損傷識(shí)別;彩繪梁;數(shù)字圖像相關(guān);工作模態(tài)

中圖分類(lèi)號(hào)：TU366.2 文獻(xiàn)標(biāo)志碼：A 文章編號(hào)：2096-6717（2022）02-0148-10

基金項(xiàng)目：國(guó)家自然科學(xué)基金（52078084、51578095）;重慶市技術(shù)創(chuàng)新與應(yīng)用發(fā)展專(zhuān)項(xiàng)（cstc2020jscx-msxm0584）

作者簡(jiǎn)介：劉綱（1977- ），教授，博士生導(dǎo)師，主要從事結(jié)構(gòu)健康監(jiān)測(cè)、振動(dòng)控制研究，E-mail：gliu@cqu.edu.cn。

Abstract： Considering that traditional damage detection needs to install contact sensors， which is not conducive for protection of ancient timber structures， a non-contact damage identification method based on digital image correlation （DIC） by using painted patterns of ancient buildings itself was proposed. Firstly， the time series images of the painted wooden beam during vibration are captured by the digital camera. Then， use the integer and sub-pixel displacement search algorithm to calculate the displacement response time history of each point on the wooden beam as well as the Fourier transform to obtain the amplitude-frequency and phase-frequency curves. Finally， the operational mode of the wooden beam is extracted by the amplitude-frequency amplitude and phase-frequency phase angle.This study then adopts the rotation mode and the curvature mode as damage characteristic indicators for damage identification by the advantage of using the DIC method to obtain modal information of all pixels along the beam length. Results from simulation and experiment of painted simply supported beam show that the operational mode obtained by the DIC method can accurately identify the damage location and the relative damage degree of the wooden beam， which provides a new technique for preventive protection of ancient timber structures.

Keywords：ancient timber structures;damage identification;painted beams;digital image correlation;operational mode

木結(jié)構(gòu)是中國(guó)古建筑的主要結(jié)構(gòu)形式，是傳承中華燦爛文明的重要載體。受木材老化、環(huán)境及生物侵蝕的長(zhǎng)期作用和突發(fā)災(zāi)害的過(guò)載作用，古建筑木結(jié)構(gòu)普遍存在結(jié)構(gòu)性能退化的問(wèn)題，并伴有不同程度的殘損[1]。但作為中國(guó)特有的歷史文化遺產(chǎn)和寶貴的旅游資源，不少古建筑仍需對(duì)游客開(kāi)放，長(zhǎng)期“帶病”工作可能加劇木結(jié)構(gòu)的損傷。因此，利用先進(jìn)技術(shù)手段實(shí)時(shí)掌握木結(jié)構(gòu)的安全狀態(tài)，是實(shí)現(xiàn)古建筑木結(jié)構(gòu)預(yù)防性保護(hù)的根本前提[2-3]。

隨著傳感及信息技術(shù)的快速發(fā)展，結(jié)構(gòu)監(jiān)測(cè)及損傷識(shí)別技術(shù)已在航空航天、機(jī)械及土木工程結(jié)構(gòu)中得到廣泛應(yīng)用，其通過(guò)傳感設(shè)備實(shí)時(shí)獲取結(jié)構(gòu)響應(yīng)和環(huán)境信息，利用損傷識(shí)別和安全評(píng)估及時(shí)掌握結(jié)構(gòu)的安全狀態(tài)，從而可有效預(yù)防結(jié)構(gòu)安全事故[4-5]。近年來(lái)，該技術(shù)被逐步引入古建筑木結(jié)構(gòu)領(lǐng)域，Jiang等[6]利用光纖光柵傳感技術(shù)對(duì)古建筑木結(jié)構(gòu)的梁偏轉(zhuǎn)、柱傾斜、關(guān)節(jié)錯(cuò)位等進(jìn)行監(jiān)測(cè)，并開(kāi)發(fā)了適用于穿斗式木結(jié)構(gòu)的監(jiān)測(cè)系統(tǒng)。楊娜等[7]基于某藏式古建木結(jié)構(gòu)監(jiān)測(cè)系統(tǒng)，采用統(tǒng)計(jì)過(guò)程控制理論對(duì)監(jiān)測(cè)數(shù)據(jù)進(jìn)行異常診斷，結(jié)果表明，游客數(shù)量對(duì)木結(jié)構(gòu)的穩(wěn)定存在較大影響。薛建陽(yáng)等[8]通過(guò)殿堂式古建筑縮尺模型振動(dòng)臺(tái)試驗(yàn)與西安鐘樓簡(jiǎn)化力學(xué)模型，提出了柱腳滑移狀態(tài)下的結(jié)構(gòu)等效剛度識(shí)別方法，利用速度、加速度傳感器的響應(yīng)識(shí)別出結(jié)構(gòu)層間抗側(cè)剛度值，結(jié)果表明，該方法可在地震倒塌前發(fā)出預(yù)警。Li等[9]在兩層木框架模型的6個(gè)梁柱節(jié)點(diǎn)處分別布設(shè)加速度傳感器，在隨機(jī)激勵(lì)下利用小波單譜熵指標(biāo)識(shí)別出了木結(jié)構(gòu)的損傷。王鑫等[10]建立了木結(jié)構(gòu)框架模型，對(duì)環(huán)境激勵(lì)下梁上各節(jié)點(diǎn)的加速度響應(yīng)信號(hào)進(jìn)行小波包分解，采用小波包能量曲率差指標(biāo)實(shí)現(xiàn)了木梁的損傷定位。

但以上監(jiān)測(cè)方法均需在古建筑木結(jié)構(gòu)中安裝傳感器，對(duì)木結(jié)構(gòu)自身有一定損傷，不利于古建筑木結(jié)構(gòu)的保護(hù)。近十年來(lái)，隨著攝像技術(shù)的快速發(fā)展，基于數(shù)字圖像相關(guān)（digital image correlation，DIC）的監(jiān)測(cè)技術(shù)得到長(zhǎng)足發(fā)展。該技術(shù)采用數(shù)字相機(jī)非接觸拍攝被測(cè)對(duì)象表面圖像，通過(guò)圖像的像素匹配獲得被測(cè)對(duì)象的全場(chǎng)變形信息，進(jìn)而提取應(yīng)變、模態(tài)等測(cè)試參量，具有非接觸、全場(chǎng)測(cè)量及便于實(shí)施等優(yōu)勢(shì)，并采用人工噴涂散斑等方式增密被測(cè)對(duì)象的表面特征，以大幅提高測(cè)量精度，目前已在現(xiàn)代結(jié)構(gòu)變形、振動(dòng)監(jiān)測(cè)中展示了巨大潛力[11]。Yu等[12]在懸臂鋁板上噴涂散斑點(diǎn)，利用高速相機(jī)進(jìn)行全場(chǎng)振動(dòng)測(cè)量，準(zhǔn)確獲取了被測(cè)物體的自振頻率、振型等模態(tài)參數(shù);Khadka等[13]在風(fēng)力渦輪機(jī)葉片表面制作散斑點(diǎn)，利用DIC方法辨識(shí)了葉片在旋轉(zhuǎn)狀態(tài)下的模態(tài)參數(shù);Molina-Viedma等[14]以一榀框架結(jié)構(gòu)為研究對(duì)象，利用DIC方法提取了隨機(jī)振動(dòng)激勵(lì)下結(jié)構(gòu)的工作模態(tài)，采用局部模態(tài)濾波方法識(shí)別和定位了框架梁的損傷。

DIC測(cè)量方法可借助被測(cè)物體表面特定的灰度特征實(shí)現(xiàn)變形、模態(tài)參數(shù)等測(cè)試，而古建筑木結(jié)構(gòu)表面往往繪制了彩色圖案。為滿(mǎn)足古建筑木結(jié)構(gòu)保護(hù)要求，筆者直接利用木結(jié)構(gòu)表面自身彩繪圖案，探索采用DIC非接觸式測(cè)量方法辨識(shí)木梁工作模態(tài)的可行性。在此基礎(chǔ)上，選取轉(zhuǎn)角模態(tài)和曲率模態(tài)作為損傷指標(biāo)，利用DIC測(cè)試的工作模態(tài)沿梁長(zhǎng)方向連續(xù)的特征，實(shí)現(xiàn)木結(jié)構(gòu)早期損傷的準(zhǔn)確定位。數(shù)值模擬和實(shí)驗(yàn)結(jié)果表明，DIC方法可辨識(shí)木梁工作模態(tài)并實(shí)現(xiàn)微小損傷定位，從而為古建筑木結(jié)構(gòu)預(yù)防性保護(hù)提供新的監(jiān)測(cè)手段。

1 基于DIC原理的模態(tài)參數(shù)辨識(shí)方法

1.1 DIC位移測(cè)試原理

數(shù)字圖像相關(guān)法（DIC）是近年發(fā)展起來(lái)的基于視覺(jué)的新測(cè)試技術(shù)。以二維圖像位移解算為例，DIC方法先建立參考圖像（未變形）與目標(biāo)圖像（變形后）統(tǒng)一的x-y坐標(biāo)系，再進(jìn)行整、亞像素位移解算。

整像素解算階段，在參考圖像中選取任意像素作為計(jì)算點(diǎn)P（x0， y0），并以計(jì)算點(diǎn)為中心選取（2M+1）×（2M+1）像素大小的矩形區(qū)域?yàn)閰⒖甲蛹?，其中M為參考子集的半徑，如圖1所示;然后，在目標(biāo)圖像中選取任意點(diǎn)為待測(cè)點(diǎn)P′（x0′， y0′），以該點(diǎn)為中心形成與參考子集相同尺寸的目標(biāo)子集，并計(jì)算相關(guān)系數(shù)CZNSSD。

式中：f（x0， y0）為參考圖像坐標(biāo)（x0， y0）處像素點(diǎn)灰度值;fm為參考子集灰度平均值;g（x0′， y0′）表示目標(biāo)圖像坐標(biāo)（x0′， y0′）處像素點(diǎn)灰度值;gm是對(duì)應(yīng)的目標(biāo)子集灰度平均值。通過(guò)整像素搜索算法，按既定搜索路徑依次篩選出相關(guān)系數(shù)極小值所對(duì)應(yīng)的待測(cè)點(diǎn)，即為參考圖像中P點(diǎn)在變形后的整像素匹配點(diǎn)，在x、y方向的整像素位移分別記為u0、v0。常用的整像素搜索算法有粗細(xì)搜索法、十字搜索法、逐點(diǎn)搜索法等，筆者采用計(jì)算效率和精度較穩(wěn)定的粗細(xì)搜索法[15]。

在亞像素計(jì)算階段，利用已獲得的整像位移值，通過(guò)插值和亞像素搜索方法進(jìn)行亞像素位移解算。對(duì)于二維圖像，需要在x、y方向分別進(jìn)行插值。其中常用的三次樣條插值函數(shù)為[16]

1.2 工作模態(tài)辨識(shí)方法

當(dāng)相機(jī)采集幀率滿(mǎn)足那奎斯特采樣定理時(shí)[19]，利用DIC方法解算出一系列圖片中目標(biāo)點(diǎn)的位移時(shí)程曲線;再通過(guò)傅里葉變換得到各像素點(diǎn)位移時(shí)程的幅頻、相頻曲線;然后根據(jù)各階固有頻率處幅頻曲線的幅值大小、相頻曲線的相位角判斷各像素點(diǎn)相對(duì)位置，從而獲得被測(cè)結(jié)構(gòu)工作模態(tài)。

多自由度體系在模態(tài)空間下的運(yùn)動(dòng)微分方程可表示為[20]

由式（13）可知，頻響函數(shù)矩陣中任意一行或一列包含所有的模態(tài)參數(shù)，例如，式（13）中第1列表示在1號(hào)測(cè)點(diǎn)y方向激勵(lì)，能夠獲得1～e號(hào)測(cè)點(diǎn)的y方向位移響應(yīng)。因此，可通過(guò)兩種方式進(jìn)行工作模態(tài)的辨識(shí)：其一，在結(jié)構(gòu)的某一固定點(diǎn)進(jìn)行激振，在其他各點(diǎn)進(jìn)行拾振，即可求得H的一列;其二，在結(jié)構(gòu)的某一固定點(diǎn)拾振，在其他各點(diǎn)激振，即可求得H的一行。

在固定點(diǎn)進(jìn)行激勵(lì)，利用DIC方法識(shí)別出各測(cè)點(diǎn)的位移響應(yīng)后，通過(guò)傅里葉變換獲取各測(cè)點(diǎn)在同一固有頻率下幅頻曲線的幅值大小并進(jìn)行歸一化處理，再利用各測(cè)點(diǎn)相頻曲線的相位角判斷其相對(duì)位置，即可實(shí)現(xiàn)工作模態(tài)的辨識(shí)。

2 基于模態(tài)參數(shù)的損傷識(shí)別方法

采用模態(tài)參數(shù)進(jìn)行損傷識(shí)別的原理在于，損傷將造成結(jié)構(gòu)質(zhì)量、剛度改變，進(jìn)而導(dǎo)致結(jié)構(gòu)模態(tài)參數(shù)發(fā)生變化，且在損傷處模態(tài)振型往往會(huì)出現(xiàn)奇異性峰值，從而判斷結(jié)構(gòu)損傷部位[21]。大量研究表明[22-24]，轉(zhuǎn)角模態(tài)和曲率模態(tài)對(duì)結(jié)構(gòu)早期損傷較為敏感，故選取這兩個(gè)指標(biāo)進(jìn)行損傷識(shí)別和定位。

2.1 轉(zhuǎn)角模態(tài)

轉(zhuǎn)角模態(tài)是位移模態(tài)的一階導(dǎo)數(shù)。因結(jié)構(gòu)損傷處剛度減少，導(dǎo)致剛度變化截面處左右側(cè)轉(zhuǎn)角不同，但結(jié)構(gòu)仍滿(mǎn)足整體變形協(xié)調(diào)條件[25]，即

式中：xl、xr分別為截面的左、右側(cè);v（x）為y方向的位移模態(tài);EI為梁的抗彎剛度;Kτ為剛度變化截面處的等效轉(zhuǎn)動(dòng)剛度，具體公式見(jiàn)文獻(xiàn)[26]。式（14）表明，在剛度變化截面處位移模態(tài)的一階導(dǎo)數(shù)不連續(xù)，即轉(zhuǎn)角模態(tài)將出現(xiàn)畸變，從而可通過(guò)該特性識(shí)別梁截面發(fā)生的損傷。

2.2 曲率模態(tài)

曲率模態(tài)為位移模態(tài)的二階導(dǎo)數(shù)，一般通過(guò)中心差分法近似計(jì)算[27]。

式中：ρ（x）為梁彎曲振動(dòng)曲線的曲率半徑;Δ為測(cè)點(diǎn)間距。當(dāng)結(jié)構(gòu)局部出現(xiàn)損傷或損傷程度增加時(shí)，梁損傷部位的抗彎剛度減少，故損傷位置處的曲率將增大，因而通過(guò)檢測(cè)曲率模態(tài)在損傷前后的變化即可判斷損傷程度與損傷位置。

3 數(shù)值模擬

3.1 數(shù)值模型

建立2 m單跨簡(jiǎn)支梁有限元模型，驗(yàn)證密集測(cè)點(diǎn)條件下?lián)p傷識(shí)別的有效性，矩形梁截面寬×高為5 cm×1 cm，如圖3所示。在梁底部設(shè)置寬1 mm的橫向通長(zhǎng)裂紋模擬損傷，以裂紋深度h與梁高H的比值表示損傷程度，以左端支座為零點(diǎn)表示損傷位置。損傷工況如表1所示。

為對(duì)比密集、稀疏測(cè)點(diǎn)下?lián)p傷識(shí)別效果，假設(shè)在簡(jiǎn)支梁上以0.25 m等間距布設(shè)加速度傳感器為稀疏測(cè)點(diǎn)情況，相當(dāng)于在被測(cè)結(jié)構(gòu)1/8等分點(diǎn)布設(shè)傳感器;而DIC方法能實(shí)現(xiàn)的測(cè)點(diǎn)最小間距由相機(jī)像素和拍攝物距決定，當(dāng)前單臺(tái)主流相機(jī)針對(duì)2 m長(zhǎng)木梁能夠識(shí)別出的物理距離（相當(dāng)于測(cè)點(diǎn)間距）可高達(dá)0.34 mm，但考慮到計(jì)算工作量及必要性，以0.05 m作為DIC測(cè)試的最大間距（密集測(cè)點(diǎn)）。若0.05 m間距無(wú)法實(shí)現(xiàn)預(yù)期的損傷識(shí)別，可進(jìn)一步選取更小的測(cè)點(diǎn)間距進(jìn)行試算。梁的損傷模型及傳感器測(cè)點(diǎn)布置如圖3所示。在數(shù)值算例中，直接采用有限元計(jì)算的模態(tài)振型進(jìn)行損傷識(shí)別，沒(méi)有采用DIC進(jìn)行模態(tài)振型識(shí)別。

3.2 損傷識(shí)別結(jié)果

考慮到古建筑木結(jié)構(gòu)梁具有截面大、跨度小的特點(diǎn)，實(shí)際結(jié)構(gòu)中往往無(wú)法有效對(duì)高階振型進(jìn)行激勵(lì)，故選取稀疏、密集測(cè)點(diǎn)下簡(jiǎn)支梁的第一階模態(tài)振型，其工況1～工況4下模態(tài)振型的變化如圖4所示，其中工況0表示無(wú)損傷狀態(tài)。稀疏測(cè)點(diǎn)下的轉(zhuǎn)角、曲率模態(tài)在損傷前后無(wú)明顯差別，而密集測(cè)點(diǎn)獲得的轉(zhuǎn)角模態(tài)與曲率模態(tài)在損傷位置處存在明顯的奇異性峰值，可以判斷損傷位置與損傷程度。

裂紋深度的增加會(huì)導(dǎo)致稀疏測(cè)點(diǎn)下的轉(zhuǎn)角模態(tài)和曲率模態(tài)幅值均產(chǎn)生微小變化，但在損傷位置并未存在奇異性峰值，故采用損傷前后模態(tài)差作為判斷指標(biāo)，以提高識(shí)別靈敏度，稀疏測(cè)點(diǎn)與密集測(cè)點(diǎn)的模態(tài)差如圖5所示，工況1/0.25表示工況1/測(cè)點(diǎn)間距為0.25 m。無(wú)噪聲情況下，稀疏測(cè)點(diǎn)的模態(tài)差在損傷位置處均未出現(xiàn)明顯的奇異性峰值，而密集測(cè)點(diǎn)的模態(tài)差可以準(zhǔn)確判斷出損傷位置;峰值的相對(duì)大小與損傷程度呈正比關(guān)系，由此可區(qū)分損傷程度，且曲率差指標(biāo)更為準(zhǔn)確。

圖6為稀疏、密集測(cè)點(diǎn)對(duì)簡(jiǎn)支梁多處損傷識(shí)別的對(duì)比結(jié)果。對(duì)于圖6（a）中的轉(zhuǎn)角模態(tài)，稀疏測(cè)點(diǎn)無(wú)法實(shí)現(xiàn)損傷識(shí)別;當(dāng)間距為0.01 m時(shí)，密集測(cè)點(diǎn)出現(xiàn)奇異性峰值，但對(duì)于靠近跨中處的損傷，轉(zhuǎn)角模態(tài)幅值出現(xiàn)減小的情況，并無(wú)奇異性峰值，對(duì)損傷不夠敏感。對(duì)于圖6（b）中的曲率模態(tài)，當(dāng)損傷發(fā)生在稀疏測(cè)點(diǎn)位置時(shí)，不能識(shí)別損傷，而當(dāng)損傷發(fā)生在測(cè)點(diǎn)之間時(shí)，會(huì)導(dǎo)致?lián)p傷的錯(cuò)誤定位以及損傷程度的辨識(shí)錯(cuò)誤;當(dāng)間距為0.05 m時(shí)，根據(jù)密集測(cè)點(diǎn)曲率模態(tài)的峰值位置和大小，能夠準(zhǔn)確區(qū)分不同位置和不同程度的損傷。

4 實(shí)驗(yàn)驗(yàn)證

選取長(zhǎng)×寬×高為1.5 m×0.05 m×0.01 m的簡(jiǎn)支木梁為研究對(duì)象，在梁側(cè)面從左至右分別粘貼旋子、蘇式、和璽3種常見(jiàn)的古建筑彩繪圖，共同作為DIC拍攝區(qū)域，3種彩繪圖的長(zhǎng)度均為0.25 m。采用1 920×1 080像素的SONY HDR-CX670攝像機(jī)進(jìn)行拍攝，幀率取為50 Hz。先以0.05 m等間距布置的測(cè)點(diǎn)作為密集測(cè)點(diǎn)，同時(shí)在梁上以0.25 m等間距布置5個(gè)加速度傳感器，設(shè)置其采樣頻率為64 Hz。在梁跨中施加豎向初位移，以木梁自由振動(dòng)方式激振。實(shí)驗(yàn)布置如圖7所示。

以數(shù)值模擬算例定義的損傷方式在木梁底部設(shè)置1 mm深度的橫向通長(zhǎng)裂縫，損傷工況如表2所示。

攝像機(jī)光軸垂直于木梁表面進(jìn)行拍攝時(shí)，表示物理位移與圖像位移之間關(guān)系的尺度因子λ=dknown/Iknown[28]，其中，dknown為物體表面已知長(zhǎng)度，mm;Iknown為在圖像上對(duì)應(yīng)的像素長(zhǎng)度，pixel。實(shí)驗(yàn)?zāi)玖旱暮穸葹?0 mm，在圖像上對(duì)應(yīng)像素長(zhǎng)度為16 pixel，故尺度因子λ=0.625 mm/pixel，像素間距為0.625 mm，因此，DIC方法能夠精確實(shí)現(xiàn)等間距密集測(cè)量。

在木梁跨中布置激光位移傳感器（LDS）以驗(yàn)證DIC方法測(cè)量位移的準(zhǔn)確性。LDS的型號(hào)為L(zhǎng)R-TB2000，采樣頻率為64 Hz。根據(jù)計(jì)算的尺度因子，將圖像位移轉(zhuǎn)換為物理位移，跨中處自由衰減的位移時(shí)程曲線如圖8（a）所示。由圖8（a）可知，激光位移傳感器與DIC方法獲取的跨中位置處位移時(shí)程曲線基本一致，從而說(shuō)明DIC測(cè)量結(jié)構(gòu)位移具有較好的準(zhǔn)確性。

采用DIC方法識(shí)別出的位移時(shí)程曲線如圖8（b）所示，圖例中數(shù)值代表距左端支座的距離。從圖8（b）可知，0.25與1.25 m處、0.5與1 m處的測(cè)點(diǎn)振動(dòng)幅值基本吻合，即對(duì)稱(chēng)位置處的位移時(shí)程曲線一致，與理論相符，表明不同類(lèi)型的古建筑彩繪圖對(duì)DIC的識(shí)別精度影響不大。

通過(guò)識(shí)別梁上密集測(cè)點(diǎn)的位移時(shí)程曲線，進(jìn)行FFT變換，獲取各像素點(diǎn)的幅頻曲線與相頻曲線。根據(jù)相位角正負(fù)與幅值大小進(jìn)行歸一化處理，獲得梁的一階工作模態(tài)，無(wú)損狀態(tài)下的工作模態(tài)如圖9所示。由該圖可知，通過(guò)DIC密集測(cè)點(diǎn)獲得的工作模態(tài)曲線更光滑，包含更多的局部信息。

根據(jù)數(shù)值模擬實(shí)驗(yàn)，考慮轉(zhuǎn)角模態(tài)與曲率模態(tài)對(duì)損傷的相對(duì)敏感程度，以及實(shí)驗(yàn)中環(huán)境噪聲與光照條件的影響，采用較為敏感的曲率差損傷指標(biāo)判斷損傷識(shí)別效果，基于第一階振型不同損傷程度的單處損傷識(shí)別結(jié)果如圖10所示。圖10（a）中，根據(jù)不同工況下、不同測(cè)點(diǎn)位置處的曲率差峰值并不能確定損傷位置在0.5 m處，相反，在其他測(cè)點(diǎn)出現(xiàn)了錯(cuò)誤的損傷定位。原因在于傳感器與木梁的接觸面較大，安裝傳感器時(shí)不能準(zhǔn)確布置在目標(biāo)采集點(diǎn)位置上，導(dǎo)致利用傳感器獲取的位移模態(tài)進(jìn)行中心差分時(shí)，測(cè)點(diǎn)間距的較大誤差會(huì)造成曲率模態(tài)峰值發(fā)生明顯改變，從而導(dǎo)致?lián)p傷位置的定位錯(cuò)誤，并且根據(jù)峰值的相對(duì)大小也難以判斷損傷程度，因此，利用傳統(tǒng)加速度傳感器難以識(shí)別木梁的損傷。

采用DIC方法識(shí)別結(jié)構(gòu)模態(tài)時(shí)，雖然圖像像素的最小間距為0.625 mm，但考慮到計(jì)算效率，先選取0.05 m作為密集測(cè)點(diǎn)，識(shí)別結(jié)果如圖10（b）所示。從圖10可知，通過(guò)DIC方法計(jì)算的曲率差在木梁距左端0.5 m處存在明顯的奇異性峰值，可實(shí)現(xiàn)損傷的準(zhǔn)確定位，且不同損傷情況下的峰值也存在差別，其大小與損傷程度呈正比關(guān)系。因此，可通過(guò)奇異性峰值大小判斷木梁損傷的相對(duì)大小。

在木梁存在兩處損傷的情況下，稀疏、密集測(cè)點(diǎn)的識(shí)別結(jié)果如圖11所示。圖11（a）中存在3處明顯峰值，其中0.75 m處的峰值最大，與實(shí)際損傷位置存在較大的誤差，原因在于根據(jù)曲率模態(tài)計(jì)算式（15），損傷造成的位移模態(tài)突變會(huì)導(dǎo)致其相鄰測(cè)點(diǎn)的曲率模態(tài)值發(fā)生改變，因此，0.5 m處的損傷導(dǎo)致位移模態(tài)在該處發(fā)生突變，0.75 m處測(cè)點(diǎn)的曲率模態(tài)由于該突變產(chǎn)生了峰值，同時(shí)0.8 m處的損傷也會(huì)對(duì)0.75 m處測(cè)點(diǎn)造成相同的影響，因此，稀疏測(cè)點(diǎn)易受多處損傷的影響而導(dǎo)致錯(cuò)誤的損傷定位，并且會(huì)造成損傷位置處的奇異性峰值疊加，從而導(dǎo)致?lián)p傷程度的誤判。

由圖11（b）可知，密集測(cè)點(diǎn)下，在損傷位置0.5、0.8 m處存在明顯奇異性峰值，可以判斷在這兩處位置存在損傷，并且對(duì)應(yīng)位置處的峰值大小與預(yù)設(shè)的裂紋深度成正比。因此，當(dāng)木梁存在多處損傷時(shí)，該方法不但可實(shí)現(xiàn)小尺度損傷的準(zhǔn)確識(shí)別，還能判斷不同損傷的相對(duì)大小。將測(cè)點(diǎn)的間距進(jìn)一步減小，能夠?qū)崿F(xiàn)兩處損傷引起的位移模態(tài)突變不會(huì)影響同一測(cè)點(diǎn)的曲率模態(tài)值，從而保證密集測(cè)點(diǎn)下的損傷準(zhǔn)確定位、損傷位置處奇異性峰值保持穩(wěn)定。

綜上，實(shí)驗(yàn)表明，使用傳統(tǒng)加速度傳感器進(jìn)行損傷識(shí)別受限于測(cè)點(diǎn)間距過(guò)大、布設(shè)位置偏離等因素，造成損傷的錯(cuò)誤定位;而DIC方法的密集測(cè)點(diǎn)能夠?qū)崿F(xiàn)多處損傷的定位以及損傷相對(duì)程度的辨別。

5 結(jié)論

針對(duì)接觸式傳感器進(jìn)行古建筑木結(jié)構(gòu)彩繪梁動(dòng)態(tài)損傷識(shí)別易造成木結(jié)構(gòu)損傷的問(wèn)題，利用木結(jié)構(gòu)梁表面自身彩繪圖案，選取不同模態(tài)參數(shù)作為損傷指標(biāo)，提出基于數(shù)字圖像相關(guān)（DIC）原理的非接觸式損傷識(shí)別方法，數(shù)值模擬和彩繪簡(jiǎn)支梁實(shí)驗(yàn)表明：

1）以古建筑彩繪圖案為對(duì)象，DIC方法能實(shí)現(xiàn)彩繪梁的無(wú)損全場(chǎng)動(dòng)態(tài)位移精準(zhǔn)測(cè)量，通過(guò)密集測(cè)點(diǎn)獲取損傷信息更多的工作模態(tài)，有效克服了傳統(tǒng)接觸式測(cè)量造成木梁表面損壞、測(cè)點(diǎn)不足的缺陷。

2）以規(guī)則排列的圖像像素為測(cè)點(diǎn)，能夠?qū)崿F(xiàn)絕對(duì)等間距測(cè)量，避免接觸式傳感器測(cè)點(diǎn)布置偏離導(dǎo)致的信息采集錯(cuò)誤。

3）DIC方法測(cè)點(diǎn)間距可達(dá)毫米級(jí)，有效避免了兩處損傷對(duì)同一測(cè)點(diǎn)的累積影響，能夠?qū)崿F(xiàn)木梁小尺度直裂縫的準(zhǔn)確定位和損傷相對(duì)大小的辨識(shí)。

4）受限于古建筑內(nèi)部場(chǎng)地大小，單個(gè)相機(jī)可能存在無(wú)法拍攝大跨度木梁的缺點(diǎn)，后期需進(jìn)一步研究使用多個(gè)相機(jī)進(jìn)行圖像拼接，從而實(shí)現(xiàn)大型古建筑結(jié)構(gòu)整體的模態(tài)測(cè)試。

參考文獻(xiàn)：

[1]薛建陽(yáng)， 吳占景， 張風(fēng)亮， 等. 碳纖維布加固古建筑木結(jié)構(gòu)基于結(jié)構(gòu)潛能和能量耗散地震破壞評(píng)估[J]. 土木建筑與環(huán)境工程， 2013， 35（6）： 103-111.

XUE J Y， WU Z J， ZHANG F L， et al. Potential and energy dissipation-based seismic damage evaluation of ancient timber structure strengthened with CFRP [J]. Journal of Civil， Architectural & Environmental Engineering， 2013， 35（6）： 103-111. （in Chinese）

[2]潘毅， 袁雙， 王慧琴， 等. 古建筑木結(jié)構(gòu)偷心造和計(jì)心造斗栱力學(xué)性能數(shù)值分析[J]. 土木建筑與環(huán)境工程， 2017， 39（5）： 9-15.

PAN Y， YUAN S， WANG H Q， et al. Numerical analysis of mechanical behavior of Tou-Xin-Zao and Ji-Xin-Zao Tou-Kung in Chinese ancient timber structures [J]. Journal of Civil， Architectural & Environmental Engineering， 2017， 39（5）： 9-15. （in Chinese）

[3]劉綱， 李孟珠， 蔣偉， 等. 古建筑彩繪梁整像素位移自適應(yīng)搜索解算方法[J]. 湖南大學(xué)學(xué)報(bào)（自然科學(xué)版）， 2020， 47（11）： 106-113.

LIU G， LI M Z， JIANG W， et al. Adaptive search algorithm method of whole-pixel deformation for ancient building painted beams [J]. Journal of Hunan University （Natural Sciences）， 2020， 47（11）： 106-113. （in Chinese）

[4]羅鈞， 劉綱， 黃宗明. 基于子結(jié)構(gòu)模型剪切型框架結(jié)構(gòu)損傷識(shí)別[J]. 振動(dòng). 測(cè)試與診斷， 2017， 37（2）： 294-300， 404.

LUO J， LIU G， HUANG Z M. Damage detection for shear frame structure based on the substructure ARMAX model [J]. Journal of Vibration，Measurement & Diagnosis， 2017， 37（2）： 294-300， 404. （in Chinese）

[5]單德山， 羅凌峰， 李喬. 橋梁健康監(jiān)測(cè)2019年度研究進(jìn)展[J]. 土木與環(huán)境工程學(xué)報(bào)（中英文）， 2020， 42（5）： 115-125.

SHAN D S， LUO L F， LI Q. State-of-the-art review of the bridge health monitoring in 2019 [J]. Journal of Civil and Environmental Engineering， 2020， 42（5）： 115-125. （in Chinese）

[6]JIANG S F， QIAO Z H， LI N L， et al. Structural health monitoring system based on FBG sensing technique for Chinese ancient timber buildings [J]. Sensors， 2019， 20（1）： 110.

[7]楊娜， 代丹陽(yáng)， 秦術(shù)杰. 古建筑木結(jié)構(gòu)監(jiān)測(cè)數(shù)據(jù)異常診斷[J]. 振動(dòng)工程學(xué)報(bào)， 2019， 32（1）： 64-71.

YANG N， DAI D Y， QIN S J. Monitoring data diagnosis of ancient wood structure [J]. Journal of Vibration Engineering， 2019， 32（1）： 64-71. （in Chinese）

[8]薛建陽(yáng)， 白福玉， 張錫成， 等. 古建筑木結(jié)構(gòu)地震損傷分析及抗側(cè)剛度識(shí)別[J]. 湖南大學(xué)學(xué)報(bào)（自然科學(xué)版）， 2019， 46（1）： 55-64.

XUE J Y， BAI F Y， ZHANG X C， et al. Seismic damage analysis and lateral stiffness identification for ancient wooden structures [J]. Journal of Hunan University （Natural Sciences）， 2019， 46（1）： 55-64. （in Chinese）

[9]LI S C， ZHANG L， TANG J， et al. Damage identification in timber structures based on wavelet singular spectrum entropy [J]. Mathematical Problems in Engineering， 2018， 2018： 1-9.

[10]王鑫， 胡衛(wèi)兵， 孟昭博. 基于小波包能量曲率差的古木結(jié)構(gòu)損傷識(shí)別[J]. 振動(dòng)與沖擊， 2014， 33（7）： 153-159， 186.

WANG X， HU W B， MENG Z B. Damage detection of an ancient wood structure based on wavelet packet energy curvature difference [J]. Journal of Vibration and Shock， 2014， 33（7）： 153-159， 186. （in Chinese）

[11]BING P. Digital image correlation for surface deformation measurement： Historical developments， recent advances and future goals [J]. Measurement Science and Technology， 2018， 29（8）： 082001.

[12]YU L P， PAN B. Single-camera high-speed stereo-digital image correlation for full-field vibration measurement [J]. Mechanical Systems and Signal Processing， 2017， 94： 374-383.

[13]KHADKA A， FICK B， AFSHAR A， et al. Non-contact vibration monitoring of rotating wind turbines using a semi-autonomous UAV [J]. Mechanical Systems and Signal Processing， 2020， 138： 106446.

[14]MOLINA-VIEDMA J， PIECZONKA L， MENDROK K， et al. Damage identification in frame structures using high-speed digital image correlation and local modal filtration [J]. Structural Control and Health Monitoring， 2020， 27（9）： e2586.

[15]ZHANG Z F， KANG Y L， WANG H W， et al. A novel coarse-fine search scheme for digital image correlation method [J]. Measurement， 2006， 39（8）： 710-718.

[16]高越. 三維數(shù)字圖像相關(guān)法的關(guān)鍵技術(shù)及應(yīng)用研究[D]. 合肥： 中國(guó)科學(xué)技術(shù)大學(xué)， 2014.

GAO Y. Research on key technologies and applications of three-dimensional digital image correlation [D]. Hefei： University of Science and Technology of China， 2014. （in Chinese）

[17]CHEN W， JIANG Z， TANG L， et al. Equal noise resistance of two mainstream iterative sub-pixel registration algorithms in digital image correlation [J]. Experimental Mechanics， 2017， 57（6）： 979-996.

[18]PAN B. An evaluation of convergence criteria for digital image correlation using inverse compositional gauss-Newton algorithm [J]. Strain， 2014， 50（1）： 48-56.

[19]BARONE S， NERI P， PAOLI A， et al. Low-frame-rate single camera system for 3D full-field high-frequency vibration measurements [J]. Mechanical Systems and Signal Processing， 2019， 123： 143-152.

[20]李德葆， 陸秋海. 實(shí)驗(yàn)?zāi)B(tài)分析及其應(yīng)用[M]. 北京： 科學(xué)出版社， 2001.

LI D B， LU Q H.Experimental modal analysis and its application [M]. Beijing： Science Press， 2001. （in Chinese）

[21]PORCU M C， PATTERI D M， MELIS S， et al. Effectiveness of the FRF curvature technique for structural health monitoring [J]. Construction and Building Materials， 2019， 226： 173-187.

[22]徐華東， 王立海， 胡志棟. 運(yùn)用曲率模態(tài)技術(shù)的木梁損傷定量識(shí)別[J]. 振動(dòng). 測(cè)試與診斷， 2011， 31（1）： 110-114，132-133.

XU H D， WANG L H， HU Z D. Quantitative identification of wood beam damage using modal curvature theory [J]. Journal of Vibration， Measurement & Diagnosis， 2011， 31（1）： 110-114，132-133. （in Chinese）

[23]徐平. 局部裂紋損傷簡(jiǎn)支梁的曲率模態(tài)特性[J]. 固體力學(xué)學(xué)報(bào)， 2011， 32（Sup1）： 171-175.

XU P. Modal curvature characteristics of a simple beam with a local damaged crack [J]. Chinese Journal of Solid Mechanics， 2011， 32（Sup1）： 171-175. （in Chinese）

[24]繆炳榮， 楊樹(shù)旺， 王名月， 等. 利用振動(dòng)響應(yīng)的多種結(jié)構(gòu)損傷識(shí)別方法比較[J]. 振動(dòng)工程學(xué)報(bào)， 2020， 33（4）： 724-733.

MIAO B R， YANG S W， WANG M Y， et al. Comparison of various structural damage identification methods using vibration response [J]. Journal of Vibration Engineering， 2020， 33（4）： 724-733. （in Chinese）

[25]管德清， 黃燕. 基于轉(zhuǎn)角模態(tài)小波分析的彈性地基梁損傷識(shí)別研究[J]. 振動(dòng)與沖擊， 2008， 27（5）： 44-47，172-173.

GUAN D Q， HUANG Y. Damage identification of elastic foundation beams based on wavelet transform of rotating angle modes [J]. Journal of Vibration and Shock， 2008， 27（5）： 44-47，172-173. （in Chinese）

[26]金明凡， 趙玫. 基于響應(yīng)的梁損傷識(shí)別[J]. 振動(dòng)與沖擊， 2006， 25（1）： 86-89，103，168.

JIN M F， ZHAO M. Damage identification of beam structure based on vibration responses [J]. Journal of Vibration and Shock， 2006， 25（1）： 86-89，103，168. （in Chinese）

[27]項(xiàng)長(zhǎng)生， 李凌云， 周宇， 等. 基于模態(tài)曲率效用信息熵的梁結(jié)構(gòu)損傷識(shí)別[J]. 振動(dòng)與沖擊， 2020， 39（17）： 234-244.

XIANG C S， LI L Y， ZHOU Y， et al. Damage identification of beam structures based on modal curvature utility information entropy [J]. Journal of Vibration and Shock， 2020， 39（17）： 234-244. （in Chinese）

[28]FENG D M， FENG M Q. Experimental validation of cost-effective vision-based structural health monitoring [J]. Mechanical Systems and Signal Processing， 2017， 88： 199-211.

（編輯 黃廷）

3866500338216