基于数字图像相关的古建筑彩绘梁损伤识别方法

刘纲 顾嘉伟 李孟珠 贺成华

摘 要：针对传统损伤识别需布设接触式传感器而不利于古建筑木结构保护的问题，利用古建筑自身彩绘图案探索基于数字图像相关原理（DIC）的非接触式损伤识别方法。通过数码相机拍摄古建筑彩绘木梁振动时序图片，采用整、亚像素位移搜索算法解算木梁上各点的位移响应时程，利用傅里叶变换得到幅频、相频曲线，进而通过幅频幅值和相频相位角提取木梁的工作模态;利用DIC方法可获取沿梁长方向所有像素模态信息的优势，引入振型转角模态、曲率模态作为损伤特征指标进行损伤识别。数值算例和彩绘简支梁实验结果表明，采用DIC方法获得的工作模态能准确识别木梁的损伤部位和损伤程度的相对大小，为古建筑木结构预防性保护提供了新技术手段。

关键词：古建筑木结构;损伤识别;彩绘梁;数字图像相关;工作模态

中图分类号：TU366.2 文献标志码：A 文章编号：2096-6717（2022）02-0148-10

基金项目：国家自然科学基金（52078084、51578095）;重庆市技术创新与应用发展专项（cstc2020jscx-msxm0584）

作者简介：刘纲（1977- ），教授，博士生导师，主要从事结构健康监测、振动控制研究，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

木結构是中国古建筑的主要结构形式，是传承中华灿烂文明的重要载体。受木材老化、环境及生物侵蚀的长期作用和突发灾害的过载作用，古建筑木结构普遍存在结构性能退化的问题，并伴有不同程度的残损[1]。但作为中国特有的历史文化遗产和宝贵的旅游资源，不少古建筑仍需对游客开放，长期“带病”工作可能加剧木结构的损伤。因此，利用先进技术手段实时掌握木结构的安全状态，是实现古建筑木结构预防性保护的根本前提[2-3]。

随着传感及信息技术的快速发展，结构监测及损伤识别技术已在航空航天、机械及土木工程结构中得到广泛应用，其通过传感设备实时获取结构响应和环境信息，利用损伤识别和安全评估及时掌握结构的安全状态，从而可有效预防结构安全事故[4-5]。近年来，该技术被逐步引入古建筑木结构领域，Jiang等[6]利用光纤光栅传感技术对古建筑木结构的梁偏转、柱倾斜、关节错位等进行监测，并开发了适用于穿斗式木结构的监测系统。杨娜等[7]基于某藏式古建木结构监测系统，采用统计过程控制理论对监测数据进行异常诊断，结果表明，游客数量对木结构的稳定存在较大影响。薛建阳等[8]通过殿堂式古建筑缩尺模型振动台试验与西安钟楼简化力学模型，提出了柱脚滑移状态下的结构等效刚度识别方法，利用速度、加速度传感器的响应识别出结构层间抗侧刚度值，结果表明，该方法可在地震倒塌前发出预警。Li等[9]在两层木框架模型的6个梁柱节点处分别布设加速度传感器，在随机激励下利用小波单谱熵指标识别出了木结构的损伤。王鑫等[10]建立了木结构框架模型，对环境激励下梁上各节点的加速度响应信号进行小波包分解，采用小波包能量曲率差指标实现了木梁的损伤定位。

但以上监测方法均需在古建筑木结构中安装传感器，对木结构自身有一定损伤，不利于古建筑木结构的保护。近十年来，随着摄像技术的快速发展，基于数字图像相关（digital image correlation，DIC）的监测技术得到长足发展。该技术采用数字相机非接触拍摄被测对象表面图像，通过图像的像素匹配获得被测对象的全场变形信息，进而提取应变、模态等测试参量，具有非接触、全场测量及便于实施等优势，并采用人工喷涂散斑等方式增密被测对象的表面特征，以大幅提高测量精度，目前已在现代结构变形、振动监测中展示了巨大潜力[11]。Yu等[12]在悬臂铝板上喷涂散斑点，利用高速相机进行全场振动测量，准确获取了被测物体的自振频率、振型等模态参数;Khadka等[13]在风力涡轮机叶片表面制作散斑点，利用DIC方法辨识了叶片在旋转状态下的模态参数;Molina-Viedma等[14]以一榀框架结构为研究对象，利用DIC方法提取了随机振动激励下结构的工作模态，采用局部模态滤波方法识别和定位了框架梁的损伤。

DIC测量方法可借助被测物体表面特定的灰度特征实现变形、模态参数等测试，而古建筑木结构表面往往绘制了彩色图案。为满足古建筑木结构保护要求，笔者直接利用木结构表面自身彩绘图案，探索采用DIC非接触式测量方法辨识木梁工作模态的可行性。在此基础上，选取转角模态和曲率模态作为损伤指标，利用DIC测试的工作模态沿梁长方向连续的特征，实现木结构早期损伤的准确定位。数值模拟和实验结果表明，DIC方法可辨识木梁工作模态并实现微小损伤定位，从而为古建筑木结构预防性保护提供新的监测手段。

1 基于DIC原理的模態参数辨识方法

1.1 DIC位移测试原理

数字图像相关法（DIC）是近年发展起来的基于视觉的新测试技术。以二维图像位移解算为例，DIC方法先建立参考图像（未变形）与目标图像（变形后）统一的x-y坐标系，再进行整、亚像素位移解算。

整像素解算阶段，在参考图像中选取任意像素作为计算点P（x0， y0），并以计算点为中心选取（2M+1）×（2M+1）像素大小的矩形区域为参考子集，其中M为参考子集的半径，如图1所示;然后，在目标图像中选取任意点为待测点P′（x0′， y0′），以该点为中心形成与参考子集相同尺寸的目标子集，并计算相关系数CZNSSD。

式中：f（x0， y0）为参考图像坐标（x0， y0）处像素点灰度值;fm为参考子集灰度平均值;g（x0′， y0′）表示目标图像坐标（x0′， y0′）处像素点灰度值;gm是对应的目标子集灰度平均值。通过整像素搜索算法，按既定搜索路径依次筛选出相关系数极小值所对应的待测点，即为参考图像中P点在变形后的整像素匹配点，在x、y方向的整像素位移分别记为u0、v0。常用的整像素搜索算法有粗细搜索法、十字搜索法、逐点搜索法等，笔者采用计算效率和精度较稳定的粗细搜索法[15]。

在亚像素计算阶段，利用已获得的整像位移值，通过插值和亚像素搜索方法进行亚像素位移解算。对于二维图像，需要在x、y方向分别进行插值。其中常用的三次样条插值函数为[16]

1.2 工作模态辨识方法

当相机采集帧率满足那奎斯特采样定理时[19]，利用DIC方法解算出一系列图片中目标点的位移时程曲线;再通过傅里叶变换得到各像素点位移时程的幅频、相频曲线;然后根据各阶固有频率处幅频曲线的幅值大小、相频曲线的相位角判断各像素点相对位置，从而获得被测结构工作模态。

多自由度体系在模态空间下的运动微分方程可表示为[20]

由式（13）可知，频响函数矩阵中任意一行或一列包含所有的模态参数，例如，式（13）中第1列表示在1号测点y方向激励，能够获得1～e号测点的y方向位移响应。因此，可通过两种方式进行工作模态的辨识：其一，在结构的某一固定点进行激振，在其他各点进行拾振，即可求得H的一列;其二，在结构的某一固定点拾振，在其他各点激振，即可求得H的一行。

在固定点进行激励，利用DIC方法识别出各测点的位移响应后，通过傅里叶变换获取各测点在同一固有频率下幅频曲线的幅值大小并进行归一化处理，再利用各测点相频曲线的相位角判断其相对位置，即可实现工作模态的辨识。

2 基于模态参数的损伤识别方法

采用模态参数进行损伤识别的原理在于，损伤将造成结构质量、刚度改变，进而导致结构模态参数发生变化，且在损伤处模态振型往往会出现奇异性峰值，从而判断结构损伤部位[21]。大量研究表明[22-24]，转角模态和曲率模态对结构早期损伤较为敏感，故选取这两个指标进行损伤识别和定位。

2.1 转角模态

转角模态是位移模态的一阶导数。因结构损伤处刚度减少，导致刚度变化截面处左右侧转角不同，但结构仍满足整体变形协调条件[25]，即

式中：xl、xr分别为截面的左、右侧;v（x）为y方向的位移模态;EI为梁的抗弯刚度;Kτ为刚度变化截面处的等效转动刚度，具体公式见文献[26]。式（14）表明，在刚度变化截面处位移模态的一阶导数不连续，即转角模态将出现畸变，从而可通过该特性识别梁截面发生的损伤。

2.2 曲率模态

曲率模态为位移模态的二阶导数，一般通过中心差分法近似计算[27]。

式中：ρ（x）为梁弯曲振动曲线的曲率半径;Δ为测点间距。当结构局部出现损伤或损伤程度增加时，梁损伤部位的抗弯刚度减少，故损伤位置处的曲率将增大，因而通过检测曲率模态在损伤前后的变化即可判断损伤程度与损伤位置。

3 数值模拟

3.1 数值模型

建立2 m单跨简支梁有限元模型，验证密集测点条件下损伤识别的有效性，矩形梁截面宽×高为5 cm×1 cm，如图3所示。在梁底部设置宽1 mm的横向通长裂纹模拟损伤，以裂纹深度h与梁高H的比值表示损伤程度，以左端支座为零点表示损伤位置。损伤工况如表1所示。

为对比密集、稀疏测点下损伤识别效果，假设在简支梁上以0.25 m等间距布设加速度传感器为稀疏测点情况，相当于在被测结构1/8等分点布设传感器;而DIC方法能实现的测点最小间距由相机像素和拍摄物距决定，当前单台主流相机针对2 m长木梁能够识别出的物理距离（相当于测点间距）可高达0.34 mm，但考虑到计算工作量及必要性，以0.05 m作为DIC测试的最大间距（密集测点）。若0.05 m间距无法实现预期的损伤识别，可进一步选取更小的测点间距进行试算。梁的损伤模型及传感器测点布置如图3所示。在数值算例中，直接采用有限元计算的模态振型进行损伤识别，没有采用DIC进行模态振型识别。

3.2 损伤识别结果

考虑到古建筑木结构梁具有截面大、跨度小的特点，实际结构中往往无法有效对高阶振型进行激励，故选取稀疏、密集测点下简支梁的第一阶模态振型，其工况1～工况4下模态振型的变化如图4所示，其中工况0表示无损伤状态。稀疏测点下的转角、曲率模态在损伤前后无明显差别，而密集测点获得的转角模态与曲率模态在损伤位置处存在明显的奇异性峰值，可以判断损伤位置与损伤程度。

裂纹深度的增加会导致稀疏测点下的转角模态和曲率模态幅值均产生微小变化，但在损伤位置并未存在奇异性峰值，故采用损伤前后模态差作为判断指标，以提高识别灵敏度，稀疏测点与密集测点的模态差如图5所示，工况1/0.25表示工况1/测点间距为0.25 m。无噪声情况下，稀疏测点的模态差在损伤位置处均未出现明显的奇异性峰值，而密集测点的模态差可以准确判断出损伤位置;峰值的相对大小与损伤程度呈正比关系，由此可区分损伤程度，且曲率差指标更为准确。

图6为稀疏、密集测点对简支梁多处损伤识别的对比结果。对于图6（a）中的转角模态，稀疏测点无法实现损伤识别;当间距为0.01 m时，密集测点出现奇异性峰值，但对于靠近跨中处的损伤，转角模态幅值出现减小的情况，并无奇异性峰值，对损伤不够敏感。对于图6（b）中的曲率模态，当损伤发生在稀疏测点位置时，不能识别损伤，而当损伤发生在测点之间时，会导致损伤的错误定位以及损伤程度的辨识错误;当间距为0.05 m时，根据密集测点曲率模态的峰值位置和大小，能够准确区分不同位置和不同程度的损伤。

4 实验验证

选取长×宽×高为1.5 m×0.05 m×0.01 m的简支木梁为研究对象，在梁侧面从左至右分别粘贴旋子、苏式、和玺3种常见的古建筑彩绘图，共同作为DIC拍摄区域，3种彩绘图的长度均为0.25 m。采用1 920×1 080像素的SONY HDR-CX670摄像机进行拍摄，帧率取为50 Hz。先以0.05 m等间距布置的测点作为密集测点，同时在梁上以0.25 m等间距布置5个加速度传感器，设置其采样频率为64 Hz。在梁跨中施加豎向初位移，以木梁自由振动方式激振。实验布置如图7所示。

以数值模拟算例定义的损伤方式在木梁底部设置1 mm深度的横向通长裂缝，损伤工况如表2所示。

摄像机光轴垂直于木梁表面进行拍摄时，表示物理位移与图像位移之间关系的尺度因子λ=dknown/Iknown[28]，其中，dknown为物体表面已知长度，mm;Iknown为在图像上对应的像素长度，pixel。实验木梁的厚度为10 mm，在图像上对应像素长度为16 pixel，故尺度因子λ=0.625 mm/pixel，像素间距为0.625 mm，因此，DIC方法能够精确实现等间距密集测量。

在木梁跨中布置激光位移传感器（LDS）以验证DIC方法测量位移的准确性。LDS的型号为LR-TB2000，采样频率为64 Hz。根据计算的尺度因子，将图像位移转换为物理位移，跨中处自由衰减的位移时程曲线如图8（a）所示。由图8（a）可知，激光位移传感器与DIC方法获取的跨中位置处位移时程曲线基本一致，从而说明DIC测量结构位移具有较好的准确性。

采用DIC方法识别出的位移时程曲线如图8（b）所示，图例中数值代表距左端支座的距离。从图8（b）可知，0.25与1.25 m处、0.5与1 m处的测点振动幅值基本吻合，即对称位置处的位移时程曲线一致，与理论相符，表明不同类型的古建筑彩绘图对DIC的识别精度影响不大。

通过识别梁上密集测点的位移时程曲线，进行FFT变换，获取各像素点的幅频曲线与相频曲线。根据相位角正负与幅值大小进行归一化处理，获得梁的一阶工作模态，无损状态下的工作模态如图9所示。由该图可知，通过DIC密集测点获得的工作模态曲线更光滑，包含更多的局部信息。

根据数值模拟实验，考虑转角模态与曲率模态对损伤的相对敏感程度，以及实验中环境噪声与光照条件的影响，采用较为敏感的曲率差损伤指标判断损伤识别效果，基于第一阶振型不同损伤程度的单处损伤识别结果如图10所示。图10（a）中，根据不同工况下、不同测点位置处的曲率差峰值并不能确定损伤位置在0.5 m处，相反，在其他测点出现了错误的损伤定位。原因在于传感器与木梁的接触面较大，安装传感器时不能准确布置在目标采集点位置上，导致利用传感器获取的位移模态进行中心差分时，测点间距的较大误差会造成曲率模态峰值发生明显改变，从而导致损伤位置的定位错误，并且根据峰值的相对大小也难以判断损伤程度，因此，利用传统加速度传感器难以识别木梁的损伤。

采用DIC方法识别结构模态时，虽然图像像素的最小间距为0.625 mm，但考虑到计算效率，先选取0.05 m作为密集测点，识别结果如图10（b）所示。从图10可知，通过DIC方法计算的曲率差在木梁距左端0.5 m处存在明显的奇异性峰值，可实现损伤的准确定位，且不同损伤情况下的峰值也存在差别，其大小与损伤程度呈正比关系。因此，可通过奇异性峰值大小判断木梁损伤的相对大小。

在木梁存在两处损伤的情况下，稀疏、密集测点的识别结果如图11所示。图11（a）中存在3处明显峰值，其中0.75 m处的峰值最大，与实际损伤位置存在较大的误差，原因在于根据曲率模态计算式（15），损伤造成的位移模态突变会导致其相邻测点的曲率模态值发生改变，因此，0.5 m处的损伤导致位移模态在该处发生突变，0.75 m处测点的曲率模态由于该突变产生了峰值，同时0.8 m处的损伤也会对0.75 m处测点造成相同的影响，因此，稀疏测点易受多处损伤的影响而导致错误的损伤定位，并且会造成损伤位置处的奇异性峰值叠加，从而导致损伤程度的误判。

由图11（b）可知，密集测点下，在损伤位置0.5、0.8 m处存在明显奇异性峰值，可以判断在这两处位置存在损伤，并且对应位置处的峰值大小与预设的裂纹深度成正比。因此，当木梁存在多处损伤时，该方法不但可实现小尺度损伤的准确识别，还能判断不同损伤的相对大小。将测点的间距进一步减小，能够实现两处损伤引起的位移模态突变不会影响同一测点的曲率模态值，从而保证密集测点下的损伤准确定位、损伤位置处奇異性峰值保持稳定。

综上，实验表明，使用传统加速度传感器进行损伤识别受限于测点间距过大、布设位置偏离等因素，造成损伤的错误定位;而DIC方法的密集测点能够实现多处损伤的定位以及损伤相对程度的辨别。

5 结论

针对接触式传感器进行古建筑木结构彩绘梁动态损伤识别易造成木结构损伤的问题，利用木结构梁表面自身彩绘图案，选取不同模态参数作为损伤指标，提出基于数字图像相关（DIC）原理的非接触式损伤识别方法，数值模拟和彩绘简支梁实验表明：

1）以古建筑彩绘图案为对象，DIC方法能实现彩绘梁的无损全场动态位移精准测量，通过密集测点获取损伤信息更多的工作模态，有效克服了传统接触式测量造成木梁表面损坏、测点不足的缺陷。

2）以规则排列的图像像素为测点，能够实现绝对等间距测量，避免接触式传感器测点布置偏离导致的信息采集错误。

3）DIC方法测点间距可达毫米级，有效避免了两处损伤对同一测点的累积影响，能够实现木梁小尺度直裂缝的准确定位和损伤相对大小的辨识。

4）受限于古建筑内部场地大小，单个相机可能存在无法拍摄大跨度木梁的缺点，后期需进一步研究使用多个相机进行图像拼接，从而实现大型古建筑结构整体的模态测试。

参考文献：

[1]薛建阳， 吴占景， 张风亮， 等. 碳纤维布加固古建筑木结构基于结构潜能和能量耗散地震破坏评估[J]. 土木建筑与环境工程， 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]潘毅， 袁双， 王慧琴， 等. 古建筑木结构偷心造和计心造斗栱力学性能数值分析[J]. 土木建筑与环境工程， 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]刘纲， 李孟珠， 蒋伟， 等. 古建筑彩绘梁整像素位移自适应搜索解算方法[J]. 湖南大学学报（自然科学版）， 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]罗钧， 刘纲， 黄宗明. 基于子结构模型剪切型框架结构损伤识别[J]. 振动. 测试与诊断， 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]单德山， 罗凌峰， 李乔. 桥梁健康监测2019年度研究进展[J]. 土木与环境工程学报（中英文）， 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]杨娜， 代丹阳， 秦术杰. 古建筑木结构监测数据异常诊断[J]. 振动工程学报， 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]薛建阳， 白福玉， 张锡成， 等. 古建筑木结构地震损伤分析及抗侧刚度识别[J]. 湖南大学学报（自然科学版）， 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]王鑫， 胡衛兵， 孟昭博. 基于小波包能量曲率差的古木结构损伤识别[J]. 振动与冲击， 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]高越. 三维数字图像相关法的关键技术及应用研究[D]. 合肥： 中国科学技术大学， 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]李德葆， 陆秋海. 实验模态分析及其应用[M]. 北京： 科学出版社， 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]徐华东， 王立海， 胡志栋. 运用曲率模态技术的木梁损伤定量识别[J]. 振动. 测试与诊断， 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]徐平. 局部裂纹损伤简支梁的曲率模态特性[J]. 固体力学学报， 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]缪炳荣， 杨树旺， 王名月， 等. 利用振动响应的多种结构损伤识别方法比较[J]. 振动工程学报， 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]管德清， 黄燕. 基于转角模态小波分析的弹性地基梁损伤识别研究[J]. 振动与冲击， 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]金明凡， 赵玫. 基于响应的梁损伤识别[J]. 振动与冲击， 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]项长生， 李凌云， 周宇， 等. 基于模态曲率效用信息熵的梁结构损伤识别[J]. 振动与冲击， 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