基于激光条带追踪的桥梁振动非接触测量方法

唐亮 吴桐 杜鹏 刘年 周志祥 戚小磊

唐亮, 吴桐, 杜鹏, 等. 基于激光条带追踪的桥梁振动非接触测量方法 [J]. 工程科学与技术, 2022, 54(6): 135-145. doi: 10.15961/j.jsuese.202100936
引用本文: 唐亮, 吴桐, 杜鹏, 等. 基于激光条带追踪的桥梁振动非接触测量方法 [J]. 工程科学与技术, 2022, 54(6): 135-145. doi: 10.15961/j.jsuese.202100936
TANG Liang, WU Tong, DU Peng, et al. Non-contact Measurement Method of Bridge Vibration Based on Laser Strip Tracking [J]. Advanced Engineering Sciences, 2022, 54(6): 135-145. doi: 10.15961/j.jsuese.202100936
Citation: TANG Liang, WU Tong, DU Peng, et al. Non-contact Measurement Method of Bridge Vibration Based on Laser Strip Tracking [J]. Advanced Engineering Sciences, 2022, 54(6): 135-145. doi: 10.15961/j.jsuese.202100936

基于激光条带追踪的桥梁振动非接触测量方法

基金项目: 国家自然科学基金项目(51708068);重庆交通大学研究生科创基金项目(2021B0004)
详细信息
    • 收稿日期:  2021-09-16
    • 网络出版时间:  2022-10-09 11:23:00
  • 作者简介:

    唐亮(1983—),男,副教授,博士. 研究方向:桥梁健康监测与损伤识别. E-mail:tangliang@cqjtu.edu.cn

    通信作者:

    吴桐, E-mail: wutong@mails.cqjtu.edu.cn

  • 中图分类号: U446.3

Non-contact Measurement Method of Bridge Vibration Based on Laser Strip Tracking

  • 摘要: 针对计算机视觉的结构位移测量方法存在对靶点要求高及光照依赖性强的问题,提出一种通过竖直投射激光条带于桥梁结构表面,追踪光带中心以分析桥梁振动的新方法。通过大功率线激光器投射激光线至梁底面或顶面,倾斜拍摄激光线连续影像,利用改进灰度重心法亚像素追踪图像上激光条带中心线位置的改变,几何换算得到桥梁振动位移时程曲线,进而分析获得桥梁结构的动力参数。对简支梁模型桥进行锤击强迫振动室内试验,利用SONY 4K摄像机进行影像采集,利用LVDT位移传感器和加速度传感器进行同步数据采集。分析结果表明:激光条带横截面光强近似服从高斯分布,通过改进灰度重心法可得到激光条带中心亚像素坐标,追踪连续影像条带中心位置改变,进而得到其振动位移时程曲线;以LVDT采集数据为参考,SONY摄像机所得动位移测量结果相对误差极值为3.17%;以加速度传感器为基准对模态分析结果进行比较,摄像机对前两阶频率识别的相对误差极值为1.81%,表明动位移测量与低阶模态参数识别都取得了理想的效果;不同于加速度传感器所得模态振型空间精度受限于传感器数量,本文方法所得模态振型空间精度可达像素级,且几乎无环境光照改变干扰试验结果。该方法为非接触测量识别桥梁自振特性提供了一条经济高效的思路。

     

    Abstract: Vision-based structural displacement monitoring methods have the disadvantages of strong dependence on target and light conditions. A new method to analyze the bridge vibration by projecting a laser strip on the surface of the bridge and tracking the center of the light band was proposed to overcome the problems. Specifically, the laser line was projected to the bottom or top surface of the beam through the high-power line laser and the continuous images of the laser line were taken by the tilted camera. Through subpixel tracking of the change of the laser stripe center line position in the image, the bridge displacement time histories can be obtained by geometric conversion, and then the dynamic parameters can be obtained. Under indoor experimental conditions, the hammer impact vibration test of a simply supported beam model bridge was carried out. The Sony 4K camera was used as the image acquisition equipment, and the linear variable differential transformer (LVDT) displacement sensors and seven accelerometers were used as the traditional data acquisition equipment. The results demonstrate that the intensity of the laser stripe cross section approximately obeys the Gaussian distribution, and the sub-pixel coordinates of its center can be obtained by the improved gray center method. The time histories of the center can be acquired by tracking the change of its coordinates. Compared with the LVDT data, the relative error extreme value of dynamic displacement measurement results obtained by the Sony camera was 3.17%. The relative error extreme value of the first two frequency identification by the Sony camera was 1.81%. Both the dynamic displacement measurement and the low-order modal parameter identification had achieved ideal results. Different from the fact that the spatial accuracy of modal shapes obtained by accelerometers was limited by the number of sensors, the spatial accuracy of modal shapes obtained by this method could reach pixel level, and there was almost no disturbance caused by the interference of ambient light. The proposed method offered a high-efficiency and cost-effective way for non-contact displacement measurement and modal identification of bridges.

     

  • 桥梁在运营过程中会受到如车辆、风和地震等来自外界的各种荷载作用。结构健康监测(structural health monitoring,SHM)是通过对桥梁整体或局部构件的性能进行监测、评估和反馈,及时发现结构存在的安全隐患,并采取相应的维修加固处理,保证桥梁安全运营的根本举措。其中,在荷载作用下反映结构刚度的位移响应是SHM系统关注的重要参数[1-3]。通过对结构位移信号的测量并进行模态分析,可以得到结构的自振频率、振型、阻尼比等动力特性参数。

    传统的结构静/动力位移测量有两种方法:1)在结构表面安装百分表、LVDT等接触式位移传感器,这种方法测量结果准确,但是需要在桥下搭设支架以获得相对结构不动点,且需要布线为数据传输与仪器供电,实际操作较繁琐。2)利用全站仪、水准仪观测结构变形,但该方法只能获得结构的静力位移,无法测量在荷载作用下结构的实时位移响应[4-6]。为了克服接触式传感器的不足,研究人员对包括全球定位系统(GPS)[7-9]、激光多普勒测振仪(LDS)[10-11]、雷达干涉测量[12]和基于计算机视觉的非接触测量方法进行研究。其中,基于计算机视觉的位移测量方法以其大尺度测量和适用范围广等优势在工业领域得到了较广泛的应用。

    基于计算机视觉的位移测量方法根据目标类型可以分为点跟踪(PT)、数字图像相关方法(DIC)或无目标跟踪方法[13]。点跟踪是在结构上安装一系列光学目标(通常是高对比度圆形点或回复反射目标),采用椭圆定位算法跟踪光学目标的中心,测量目标中心的位置改变。Morlier等[14]利用2D点跟踪进行桥梁动位移响应测量,并进一步处理以提取固有频率。Busca等[15]在点跟踪方法中使用了方形光学目标跟踪其角点,测量桥上多点的动位移。Wu等[16]引入新的校准方法,提高摄影测量的精度,并对位于振动台上的结构进行2维测量,用点跟踪法测得结构响应的FFT谱图显示结构的前几阶频率。楚玺等[17]将欧拉运动放大算法应用于桥梁振动分析,通过对振动放大后的数字影像进行帧数分解,获得常规视频采样频率下的桥梁试件振动时间历程序列;最后,通过分析图像序列中标记点的边缘特征,得到试验桥梁标记点的动态位移测量结果。Dong等[18]提出通过追踪光学目标的多点同步测量结构动态位移的方法,利用实测多点动位移和FFT识别结构模态参数。Feng等[19]利用改进的模板匹配算法,通过追踪标志点实时测量一座铁路桥梁的动位移,并识别其频率。

    DIC基于追踪连续图案的灰度变化,对变形前图像中的感兴趣区域(ROI)进行网格划分,将每个子区域当作刚性运动;再针对每个子区域,通过一定的搜索方法按预先定义的相关函数进行相关计算,在变形后图像中寻找与该子区域的互相关系数为最大值的区域,即该子区域变形后的位置,进而获得该子区域的位移;以此对全部子区域进行计算,即可获得全场的变形信息。由于早期的DIC系统需要较大的变形值来识别响应,因此,早期与模态识别相关的研究大多集中在单一模态激励上。Warren等[20]利用振动台以共振频率激励结构,利用DIC系统测量激励过程中结构的变形形状,进而识别结构的模态。Wang等[21]利用2维DIC系统,测量复合材料板的振动模态,并与有限元模型进行比较。Tregidgo等[22]使用DIC测量飞机机翼在不同工况下的变形,并使用DIC测量机翼的响应。在这种情况下,振动的振幅通常小于谐波激励,模态提取较为困难。另外,许多DIC软件程序不具备对频域数据的信号处理功能,且由于DIC中的数据是点云的形式,商业模态分析软件程序无法进行数据处理。因此,研究人员使用定制的脚本对数据进行后处理以提取模态形状。Molina−Viedma等[23]利用DIC与运动放大算法,实现对结构微小振动的测量;随后,Molina−Viedma等[24-25]利用高速3维DIC系统和局部频域滤波算法,实现框架结构和悬臂梁的振动测量与损伤识别。Beberniss等[26]利用白噪声激励刚度较小的金属结构;为了提取振型,其利用多项式对DIC数据进行曲线拟合,并用这些多项式更新有限元模型。Hagara[27]和Trebuňa[28]等使用白噪声激励金属薄板,使用DIC测量结构的响应,并开发了一个能够从DIC中直接提取模型模态参数的软件程序。

    无目标方法是通过跟踪结构内部特征或边缘等自然特征来识别并跟踪目标对象[29]。Busca等[30]利用模板匹配、边缘检测和DIC 3种模板匹配算法,开发了一种基于视觉的位移传感器系统。利用视觉传感器跟踪固定在桥梁上的高对比度目标板,测量桥梁的竖向位移。Ribeiro等[31]利用随机抽样一致性算法(RANSAC),测量铁路桥梁的动态位移;基于鲁棒方向码匹配算法,通过跟踪结构表面的自然目标,开发一种用于实时位移测量的视觉传感器系统,避免了人工靠近结构物安装标靶。郑佳艳等[32]将欧拉运动放大算法与边缘检测算法相结合,测量了实验室一座自锚式悬索桥缩尺模型的动位移。Son等[33]提出一种摄像机对振动位移进行远程测量的边缘检测算法,并利用所提出的测量方法对圆柱结构振动位移的测量分辨率和测量误差进行估计。Yoon等[34]通过无人机搭载摄像机追踪结构上自然特征点,测量在荷载作用下铁路桁架桥的位移响应。Dong等[35]提出针对不同结构设置不同的虚拟标记点,利用光流追踪算法识别在人群跳跃工况下体育场的模态参数。

    然而,基于目标跟踪的方法通常需要人工靠近结构表面喷涂散斑或粘贴标识点等方式,以获得可跟踪的外部目标,在实际桥梁的布置上会花费大量时间,可测量的位移点数与标志点的数量有关。应用边缘检测、自然特征点跟踪及运动放大算法可实现无目标跟踪,但对光照和背景条件又有较高的要求;如果在测量过程中光照条件不理想或发生改变,该方法效果差强人意。

    针对以上存在的不足,本文提出一种通过在桥梁结构表面投射激光条带,追踪光带中心,测量桥梁位移并识别结构模态参数的方法。具体做法为:通过大功率线激光器投射激光线至桥梁底面或顶面,倾斜拍摄激光线连续影像;利用改进灰度重心法亚像素追踪拟合激光条带中心线位置在图像上的改变,几何换算得到桥梁振动位移时程曲线,进而分析获得桥梁结构的动力参数。激光线的远距离投射不需要靠近桥梁结构表面,且不依赖于环境光照和背景条件,可实现性好。在室内试验条件下,对简支梁模型桥进行锤击强迫振动试验,以SONY 4K摄像机作为影像采集设备,LVDT位移传感器和加速度传感器作为传统数据采集设备,进行同步数据采集试验,测量动位移响应,识别模型桥模态参数,并将测量结果与传统传感器进行对比。

    基于激光条带追踪的结构振动测量大致分为视频采集、图像预处理、激光条带中心提取、像素点位移计算和真实位移计算5个步骤,其流程图如图1所示,测量原理如图2所示。由图2(a)可见:线激光器在结构物正上方竖直向下投射,在结构表面形成一具有一定宽度的激光条带;当结构发生竖向振动时,激光条带会伴随结构一起振动;将摄像机倾斜拍摄,则激光条带的竖向振动在图像平面表现为斜向振动,其在结构表面的横向位置也会发生改变。图2(b)图2(a)的侧视图,假设A为投射在物体表面的激光线上一点,A'为图像中与A对应的点。当结构发生变形时,点A运动到点B,在图像平面表现为A'运动到B'。图像预处理包括图像裁剪、图像滤波和图像灰度化。采用灰度重心法提取激光条带中心线,将提取的中心线在时间尺度上进行串联,即可得到激光线各列像素中心点的像素位移时程曲线;再根据图像平面与真实物理平面的尺度转化因子,就能够计算得到结构表面激光条带中心线的真实位移时程曲线,也就是结构的位移时程曲线。

    图  1  基于激光条带中心追踪的结构振动测量流程图
    Fig.  1  Overall flowchart of vibration measurement method based on laser stripe tracking
    下载: 全尺寸图片
    图  2  基于激光条带追踪的结构振动测量原理图
    Fig.  2  Principle diagram of structural vibration measurement based on laser strip tracking
    下载: 全尺寸图片

    数字立体摄影测量是一种非接触测量方法,利用相机记录的一系列图像识别点、图案、图像中的特征或物体边界的坐标[29]。在测试期间,将相对位置固定在适当位置,同时记录图像(或视频)并保存。在图像中找到光学目标或特征,并跟踪以识别变形。

    图2(b)可见,假设相机针孔与物体表面距离为D,相机焦距为f,相机光轴与物体表面法线夹角为θ $ {d_{{\rm{pixel}}}} $ 为图像中单位像素对应的尺寸,若AB长度已知,A点与B点的真实坐标分别为zAzB,在图像平面对应点A'B'的像素坐标分别为IA'IB',则比例因子的计算公式为:

    $$ S {F_1} = \frac{{{{\textit{z}}_A} - {{\textit{z}}_B}}}{{{I_{A'}} - {I_{B'}}}}\qquad $$ (1)

    当倾斜角度 $\theta {\text{ = }}0$ 时,式(1)可简化为:

    $$ S {F_1} = \frac{D}{f}{d_{{\rm{pixel}}}} \qquad $$ (2)

    当倾斜角度 $\theta \ne 0$ 时,根据相似三角形原理,zAzB可以表示为:

    $$\left\{ \begin{aligned} &{{\textit{z}}_A} = \frac{{D{{\textit{z}}_{iA'}}}}{{f{{\cos }^2}\theta - {{\textit{z}}_{iA'}}\cos \;\theta \sin \;\theta }},\\ &{{\textit{z}}_B} = \frac{{D{{\textit{z}}_{iB'}}}}{{f{{\cos }^2}\theta - {{\textit{z}}_{iB'}}\cos \;\theta \sin \;\theta }} \end{aligned} \right.$$ (3)

    式中, $ {{\textit{z}}_{iA'}} $ $ {{\textit{z}}_{iB'}} $ 分别为图像平面的尺寸坐标, ${{\textit{z}}_{iA'}} = {I_{A'}}{d_{{\rm{pixel}}}}$ ${{\textit{z}}_{iB'}} = {I_{B'}}{d_{{\rm{pixel}}}}$ 。当θ很小(即 $\sin \;\theta \approx 0$ ),且有 $ {{\textit{z}}_{iA'}} \ll f $ $ {{\textit{z}}_{iB'}} \ll f $ 时,式(3)中的比例因子可以根据相机的内在参数和相机与结构之间的外在参数进一步简化:

    $$ \begin{aligned}[b] S {F_2} =& \frac{1}{{{I_{A'}} - {I_{B'}}}}\left( \frac{{D{{\textit{z}}_{iA'}}}}{{f{{\cos }^2}\theta - {{\textit{z}}_{iA'}}\cos \;\theta \sin\; \theta }} - \right. \\& \left.\frac{{D{{\textit{z}}_{iB'}}}}{{f{{\cos }^2}\theta - {{\textit{z}}_{iB'}}\cos \;\theta \sin\; \theta }} \right) \approx \frac{D}{{f{{\cos }^2}\theta }}{d_{{\rm{pixel}}}} \end{aligned} $$ (4)

    以Sony FDR AX–700摄像机为例,设图像分辨率为1 920×1 080,传感器高为8.8 mm,当 $\theta {\text{ = }}0$ 时,根据式(2)计算在不同焦距工况下比例因子SF随着摄像机与物体距离变化的关系如图3所示;当 $\theta \ne 0$ ,根据式(4)计算摄像机与物体距离为2.0 m时,不同焦距工况下比例因子随倾斜角度的变化关系如图4所示。

    图  3  不同焦距下比例因子与测量距离关系
    Fig.  3  Relationship between SF and measurement distance under different focal lengths
    下载: 全尺寸图片
    图  4  不同焦距下比例因子与倾斜角度关系
    Fig.  4  Relationship between SF and tilt angle under different focal lengths
    下载: 全尺寸图片

    图3可见:当焦距不变时,SF随测距增加线性增加;当测距不变时,SF随焦距增加而增加。由图4可见:随着倾斜角度的增大,比例因子SF表现为指数增长;但当 $\theta < {30^{{\circ }}}$ 时,SF增大不明显,可以大致认为相邻两点的比例因子不变。所以,根据式(2)和(4),当结构发生变形时,其真实变形值可以根据图像中的变形值×比例因子SF1SF2计算得到[19]。并且,实际变形值与图像中的变形值呈线性关系,即通过数字图像上点的位置变化可线性计算得到被测物表面点的变形值。因此,由图像上激光线的模态信息获知结构的模态信息是可行的。

    由于图像中包含了大量的背景及其他光噪声,为了提高光条中心定位的精度,首先需要对图像进行预处理,步骤如下:

    1)对图像感兴趣区域(ROI)进行裁剪。根据结构振动的先验知识可知,对于大部分桥梁而言,其振动幅度微小,且多为竖向振动,则激光条带中心位置的变化幅度范围较小,且仅在竖向发生变位。由于激光条带在图像长度方向所占尺寸远大于宽度方向尺寸,所以可以将图像裁剪长度设为激光条带全长,宽度为其宽度的5~10倍。

    2)采用中值滤波对图像进行降噪处理,其基本原理是把数字图像或数字序列中一点的值用该点的一个邻域中各点值的中值代替,让周围的像素值接近真实值,从而消除孤立的噪声点。

    3)对图像进行灰度变化。一张RGB图像是一个M×N×3的彩色像素数组,用矩阵可以表示为:

    $$ {\boldsymbol{g}}\left( {x,y} \right) = \left[ {\begin{array}{*{20}{c}} {R\left( {0,0} \right),G\left( {0,0} \right),B\left( {0,0} \right)}&{}& \cdots &{R\left( {0,N - 1} \right),G\left( {0,N - 1} \right),B\left( {0,N - 1} \right)} \\ {R\left( {1,0} \right),G\left( {1,0} \right),B\left( {1,0} \right)}&{}& \cdots &{R\left( {1,N - 1} \right),G\left( {1,N - 1} \right),B\left( {1,N - 1} \right)} \\ \vdots &{}& & \vdots \\ {R\left( {M - 1,0} \right),G\left( {M - 1,0} \right),B\left( {M - 1,0} \right)}&{}& \cdots &{R\left( {M - 1,N - 1} \right),G\left( {M - 1,N - 1} \right),B\left( {M - 1,N - 1} \right)} \end{array}} \right] $$ (5)

    对图像进行灰度转换是将3通道的彩色图像转换成2通道的灰度图像。通过灰度转换可以使图像中灰度值较大的区域变亮,灰度值较小的区域变暗,从而突出激光条带的灰度值比例,有助于进行下一步的计算和中心线的提取。目前常用的图像灰度处理方法主要包括平均值法、单一灰度法和加权平均值法。其中:平均值法是将3通道灰度的平均值作为转换后的灰度值;单一灰度法是仅保留3个通道中某一特定通道的灰度值;加权平均值法是将RGB 3分量以不同的权重进行加权平均,目前公认权重取值分别为0.299、0.587和0.114,则加权平均法的公式为:

    $$ g' = 0.299R + 0.587G + 0.114B $$ (6)

    式中,RGB分别为彩色图像中3通道的灰度值,g'为处理后的灰度。

    利用3种方法分别对试验的图像进行灰度化,如图5所示。从图5可以看出,由于激光线为红色,当仅保留R分量时,激光线与背景的对比度最高,便于后续计算。因此,在后续图像灰度化均采用R通道单一灰度值法。

    图  5  不同灰度转换法得到的灰度图像
    Fig.  5  Gray images obtained by different gray conversion methods
    下载: 全尺寸图片

    线激光器发出的光束投射到结构物表面时,形成一条一定宽度的光条,该光条中心线位置附近光强高,两侧边缘光强低。要想高精度追踪激光线的改变,需要高精度定位各时刻光条中心线位置。具体方法有基于光条灰度极大值附近区域灰度拟合曲线的极大值位置法(包括灰度极值法、抛物线拟合法和基于Hessian矩阵的Steger法)和基于光条区域灰度拟合曲线的灰度重心法(包括自适应阈值重心法、自适应重心迭代法、多项式拟合法、NURBS曲线插值法和方向模板法等)两大类[36],本文采用改进的灰度重心法提取激光条带中心坐标。

    g(x,y)表示经过预处理后的一幅灰度图像不同像素点的灰度值,对于一幅像素为M×N的图像,有:

    $$ {\boldsymbol{g}}\left( {x,y} \right) = \left[ {\begin{array}{*{20}{c}} {g\left( {0,0} \right)}&{g\left( {0,1} \right)}& \cdots &{g\left( {0,N - 1} \right)} \\ {g\left( {1,0} \right)}&{g\left( {1,1} \right)}& \cdots &{g\left( {1,N - 1} \right)} \\ \vdots & \vdots & & \vdots \\ {g\left( {M - 1,0} \right)}&{g\left( {M - 1,1} \right)}& \cdots &{g\left( {M - 1,N - 1} \right)} \end{array}} \right] $$ (7)

    传统灰度重心法是逐一计算激光条带各列的灰度重心,将每一列计算得到的激光条带灰度重心的纵坐标作为中心坐标。假设第k列的纵坐标分别为 ${y_i} $ $ (i = 0,1,2, \cdots ,M)$ ,对应灰度值分别为 ${g_i}$ $(i = 0,1,2, \cdots, M) $ ,则该列的灰度重心计算公式为:

    $$ y = \frac{{\displaystyle\sum\limits_{i = 0}^M ({{g_i} \cdot {y_i}}) }}{{\displaystyle\sum\limits_{i = 0}^M {{g_i}} }}\qquad $$ (8)

    采用传统灰度重心法提取激光条带中心线,将各列像素的灰度重心点沿长度直接连线。在实际操作过程中,由于激光线沿长度方向的截面宽度与亮度分布不均匀,使该方法极易受到噪声影响,且参与计算的像素点数量会发生改变。因此,为了提高提取激光条带中心线的精度,可通过设置动态灰度阈值 ${T_i}(i = 0,1,2, \cdots ,N)$ ,对图像进行自适应灰度矩阵重置,重新确定各列参与计算的光条像素点数量及边界像素:

    $$ {{{\boldsymbol{g}}}_{i}^{\prime }}\left(x,y\right)=\left\{ \begin{array}{l}{{\boldsymbol{g}}}_{i}\left(x,y\right),{\boldsymbol{g}}\left(x,y\right)\ge {T}_{i};\\ 0,{\boldsymbol{g}}\left(x,y\right) < {T}_{i}\end{array}\right. $$ (9)

    式中, $ {\boldsymbol{g}}\left( {x,y} \right) $ 为原始图像的灰度矩阵。

    激光条带截面光强分布近似于高斯分布[37]

    $$ I\left( y \right) = A{{\rm{e}}^{ - \frac{{{{\left( {y - {{\bar y}_0}} \right)}^2}}}{{{\sigma ^2}}}}} $$ (10)

    式中:y为图像平面上2维笛卡尔坐标系的y坐标;A为激光条带灰度值; $ {\bar y_0} $ 为强度分布中心的y坐标,即高斯曲线的中心坐标;σ为激光条带强度分布的宽度。当 $ y = \bar y $ 时, $ I\left( y \right) $ 达到最大值。

    因普通线结构光激光器的光强沿激光条带的纵向近似表现为高斯分布,即中间部分亮度大,两端亮度较低,则动态阈值的数学模型可以表示为:

    $$ T\left( x \right) = X{{\rm{e}}^{\frac{{{{\left( {x - b} \right)}^2}}}{{{c^2}}}}} $$ (11)

    式中,Xbc均为常数,且X>0。可以先计算部分离散阈值 $T$ ,根据高斯拟合确定参数Xbc,从而计算各列的阈值 ${T_i}$

    按照等间距选取y列像素,计算离散阈值 $T$ :首先,利用极值法得到光带中心的大致位置,设第i列激光宽度为s,在中心位置上下各取宽度s的像素点,使区域宽度为激光条带宽度的两倍,保证光带横向全部在该区域内;同时,再向左右分别扩展s列,可得到s×2s大小的矩形区域,利用最大类间方差法计算每一列的阈值 $T$ ;再根据式(11)进行高斯拟合,得到阈值函数T(x),确定各列像素的阈值,并确定各列参与计算的光条像素点数量及边界像素;然后,按照式(8)逐列计算光条的中心坐标。各列计算完毕,则可以得到光条的中心线提取结果。

    试验所用简支铝梁跨径为3.84 m;横截面为π形截面,顶板宽15 cm,厚3 mm;肋板高5 cm,厚3 mm;肋板间距为5 cm。材料弹性模量为7.2×104 MPa, 密度为2.73 g/cm3。采用索尼4K摄像机(FDR–AX700)追踪激光线位置的改变。摄像机拍摄视频为1 920×1 080像素、帧率为100 fps。同时,在简支梁跨中布置1个LVDT位移传感器,在梁下方等间距布置7个加速度传感器,采样频率均设置为200 Hz;配合江苏东华测试公司的动态采集仪DH 5902N采集位移响应作为验证。简支梁模型与试验布置如图6所示。

    图  6  简支梁试验现场照片
    Fig.  6  Photos of the simply supported beam test
    下载: 全尺寸图片

    俯视角度约30°时,分别进行拍摄距离为0.5和2.0 m自然光照下和2.0 m完全黑暗环境下共3种工况的锤击振动试验。试验过程中:用支架固定摄像机,调整好摄像机的倾斜角度之后,利用分辨率为0.01°的数字倾角仪测量摄像机倾角并记录;保持摄像机空间位置与角度不变,为防止手动拍摄带来的扰动,采用蓝牙遥控器控制拍摄。单次试验采用锤击法在模型梁跨径5/16处进行多次锤击激励,利用结构受迫振动后的自由衰减信号进行振动分析。

    将所录制的视频逐帧转换为灰度图像(图7)。激光线在简支梁上的投影为一灰度值明显大于其他区域的激光条带。选取第800列的激光条带分析截面光强,该截面光强分布如图8所示,利用高斯分布进行拟合,其相关系数为R=0.947 2。采用低通和中值滤波及图像分割后,得到激光条带,作为光条中心计算的基础。然后,对序列图像依次利用改进灰度重心法计算每列中的光强峰值点,可以得到各列光条的中心位置坐标(图9)。令 $ {h_x}\left( {y,t} \right) $ 为第x列中光强中心随时间的变化,可得到其振动时程曲线。

    图  7  激光线图像及灰度图像
    Fig.  7  Laser line image and its gray image
    下载: 全尺寸图片
    图  8  截面光强分布(x=800)
    Fig.  8  Light intensity distribution (x=800)
    下载: 全尺寸图片
    图  9  激光条带中心提取
    Fig.  9  Extraction of laser stripe center
    下载: 全尺寸图片

    将激光线中与LVDT位移传感器对应点的振动时程曲线与比例因子相乘,可得到该点位移时程曲线。将3种工况摄像机和LVDT的采集起点统一,各工况位移时程曲线如图10所示。由图10可以看出,由摄像机测量得到结构位移的总体趋势与位移传感器测量结果基本一致,正确反映了桥梁的振动幅值、振动衰减特性、起振点等信息。各次激振摄像机测量得到的位移与位移传感器结果吻合较好。

    图  10  各工况跨中测点位移时程曲线
    Fig.  10  Displacement time history curves of mid-span in each case
    下载: 全尺寸图片

    为了从整体上评估摄像机测量动位移的精度,对位移传感器测量数据进行重采样,与摄像机采样频率保持一致为100 Hz。使用归一化均方根误差( $ \delta _{\rm{NRMS}} $ )进行误差分析,计算公式为:

    $$ \delta_{\rm{NRMS}} = \dfrac{{\sqrt {\dfrac{1}{n}\displaystyle\sum\limits_{i = 1}^n {{{\left( {{x_i} - {y_i}} \right)}^2}} } }}{{{y_{\max }} - {y_{\min }}}} \times 100\% $$ (12)

    式中,n为测量数据的总数, ${x_i}$ ${y_i}$ 分别为由位移传感器和摄像机测量得到的第i个位移值, ${y_{\max }}$ ${y_{\min }}$ 分别为摄像机测量得到的最大位移值和最小位移值。各工况不同方法测得的位移值见表1。由表1可见:3种工况的 $\delta _{\rm{NRMS}}$ 分别为2.49%、2.58%和1.80%;通过分析各次锤击后的位移数据,可以得到工况1测量最大误差为1.82%,工况2最大误差为3.17%,工况3最大误差为2.49%;结构的绝对位移越大,其相对误差越小;在相同测距的情况下,黑暗环境的测量精度高于自然光条件。由此说明,利用本文算法进行动位移测量的精度较高,其结果与实测值吻合良好,满足工程精度要求。

    表  1  各工况位移测量值
    Table  1  Displacement measurements of each working condition
    试验
    工况
    激振
    顺序
    传感器位移
    测量值/mm
    摄像机测量值 $\delta _{\rm{NRMS}}$/%
    位移值/mm 误差/%
    工况1 1 1.702 1.671 1.82 2.49
    2 1.742 1.749 0.40
    3 1.664 1.669 0.30
    工况2 1 1.621 1.573 2.96 2.58
    2 1.484 1.437 3.17
    3 1.921 1.896 1.30
    工况3 1 2.329 2.271 2.49 1.80
    2 3.262 3.278 0.49
    3 2.874 2.869 0.17

    对所提取的振动时程曲线进行模态分析,得到试验梁在各工况的自振频率,如图11所示。由图11可以看出,幅–频曲线中主峰峰值明显突出,且各工况的峰值频率表现出高度的一致性。

    图  11  各工况跨中测点的功率谱密度和幅值谱曲线
    Fig.  11  Power spectral density and amplitude spectrum of mid−span point in each case
    下载: 全尺寸图片

    将识别的峰值频率列入表2表2中,3种工况下,摄像机识别出的1阶频率为6.832、6.902和6.908 Hz,2阶频率为27.520、27.660和27.664 Hz。将由位移传感器测试的数据所得结构自振频率作为参考值,可以看出:利用摄影机所得到两次试验3个工况第1阶自振频率的相对误差分别为1.81%、0.80%和0.71%;第2阶自振频率的相对误差均为0.25%左右;第3阶频率的识别效果不甚理想,因此并未列出。这与结构的刚度及摄像机的像素有关:由采样定理可知,由于摄像机采样频率fs为100 fps,最高分析频率为 ${f_{\rm{a}}} = {f_{\rm{s}}}/2 = 50\;{\rm{Hz}}$ ;而由加速度传感器所得数据及简支梁的频率计算公式可知,结构第3阶自振频率大于50 Hz。因此,用摄像机无法得到结构的第3阶频率,可通过高速摄像机替换。

    表  2  实测频率值及相对误差
    Table  2  Frequency values and relative errors
    试验
    工况
    频率
    阶数
    参考值
    /Hz
    摄像机频率
    /Hz
    误差
    /%
    工况1 1 6.958 6.832 1.81
    2 27.590 27.520 0.25
    工况2 1 6.958 6.902 0.80
    2 27.590 27.660 0.25
    工况3 1 6.958 6.908 0.71
    2 27.590 27.664 0.27

    通过对所有列进行模态分析,可以得到简支梁模型高密度的模态振型,图12为工况2由激光线和7个加速度传感器所得简支梁前两阶位移模态。

    图  12  工况2各阶模态振型
    Fig.  12  Mode shapes of case 2
    下载: 全尺寸图片

    图12可以看出:加速度传感所得振型的空间精度受限于传感器数量,而采用激光线所得振型空间精度可达到像素级;对于小尺度结构而言,传统接触式传感器会改变结构的质量分布,进一步改变结构的刚度分布及模态参数。所提出方法无需在结构表面布置传感器,因此,对结构进行模态参数识别不会改变结构的质量分布。

    利用激光条带中心追踪,进行桥梁振动测量与分析,得到如下结论:

    1)图像上点的变形与被测结构上对应点的实际变形值有严格几何意义上的线性关系;该线性关系与拍摄距离、角度及相机参数有关。通过数字图像上点的位置变化,可计算得到被测物表面点的变形值,由图像上激光线的模态信息获知结构的模态信息是可行的。

    2)光带的横截面光强近似服从高斯分布,利用改进灰度重心法得到激光条带中心,实现了所拍摄影像的亚像素计算和分析,通过追踪连续影像条带中心位置改变,可得到其振动位移时程曲线。

    3)摄像机测量所得动位移结果与LVDT位移传感器一致性良好,测量误差极值为3.17%,可满足工程测量精度的要求;通过对振动位移时程曲线分析,得到结构振动低阶频率,前两阶频率相对误差极值为1.81%。通过模态分析得到了试验梁的前两阶像素级精度模态振型,与加速度传感器结果吻合良好,但受限于设备拍摄帧率与像素,无法分析获得结构更高阶模态,可通过采用高速摄像机改善。

  • 图  1   基于激光条带中心追踪的结构振动测量流程图

    Fig.  1   Overall flowchart of vibration measurement method based on laser stripe tracking

    下载: 全尺寸图片

    图  2   基于激光条带追踪的结构振动测量原理图

    Fig.  2   Principle diagram of structural vibration measurement based on laser strip tracking

    下载: 全尺寸图片

    图  3   不同焦距下比例因子与测量距离关系

    Fig.  3   Relationship between SF and measurement distance under different focal lengths

    下载: 全尺寸图片

    图  4   不同焦距下比例因子与倾斜角度关系

    Fig.  4   Relationship between SF and tilt angle under different focal lengths

    下载: 全尺寸图片

    图  5   不同灰度转换法得到的灰度图像

    Fig.  5   Gray images obtained by different gray conversion methods

    下载: 全尺寸图片

    图  6   简支梁试验现场照片

    Fig.  6   Photos of the simply supported beam test

    下载: 全尺寸图片

    图  7   激光线图像及灰度图像

    Fig.  7   Laser line image and its gray image

    下载: 全尺寸图片

    图  8   截面光强分布(x=800)

    Fig.  8   Light intensity distribution (x=800)

    下载: 全尺寸图片

    图  9   激光条带中心提取

    Fig.  9   Extraction of laser stripe center

    下载: 全尺寸图片

    图  10   各工况跨中测点位移时程曲线

    Fig.  10   Displacement time history curves of mid-span in each case

    下载: 全尺寸图片

    图  11   各工况跨中测点的功率谱密度和幅值谱曲线

    Fig.  11   Power spectral density and amplitude spectrum of mid−span point in each case

    下载: 全尺寸图片

    图  12   工况2各阶模态振型

    Fig.  12   Mode shapes of case 2

    下载: 全尺寸图片

    表  1   各工况位移测量值

    Table  1   Displacement measurements of each working condition

    试验
    工况
    激振
    顺序
    传感器位移
    测量值/mm
    摄像机测量值 $\delta _{\rm{NRMS}}$/%
    位移值/mm 误差/%
    工况1 1 1.702 1.671 1.82 2.49
    2 1.742 1.749 0.40
    3 1.664 1.669 0.30
    工况2 1 1.621 1.573 2.96 2.58
    2 1.484 1.437 3.17
    3 1.921 1.896 1.30
    工况3 1 2.329 2.271 2.49 1.80
    2 3.262 3.278 0.49
    3 2.874 2.869 0.17

    表  2   实测频率值及相对误差

    Table  2   Frequency values and relative errors

    试验
    工况
    频率
    阶数
    参考值
    /Hz
    摄像机频率
    /Hz
    误差
    /%
    工况1 1 6.958 6.832 1.81
    2 27.590 27.520 0.25
    工况2 1 6.958 6.902 0.80
    2 27.590 27.660 0.25
    工况3 1 6.958 6.908 0.71
    2 27.590 27.664 0.27
  • [1] Wang Hao,Li Aiqun,Guo Tong,et al.Establishment and application of the wind and structural health monitoring system for the Runyang Yangtze River bridge[J].Shock and Vibration,2014,2014:421038. doi: 10.1155/2014/421038
    [2] Feng Dongming,Feng M Q.Model updating of railway bridge using in situ dynamic displacement measurement under trainloads[J].Journal of Bridge Engineering,2015,20(12):04015019.1–04015019.12. doi: 10.1061/(asce)be.1943-5592.0000765
    [3] 叶肖伟,董传智.基于计算机视觉的结构位移监测综述[J].中国公路学报,2019,32(11):21–39. doi: 10.19721/j.cnki.1001-7372.2019.11.002

    Ye Xiaowei,Dong Chuanzhi.Review of computer vision-based structural displacement monitoring[J].China Journal of Highway and Transport,2019,32(11):21–39 doi: 10.19721/j.cnki.1001-7372.2019.11.002
    [4] Psimoulis P A,Stiros S C.Measuring deflections of a short-span railway bridge using a robotic total station[J].Journal of Bridge Engineering,2013,18(2):182–185. doi: 10.1061/(asce)be.1943-5592.0000334
    [5] Fukuda Y,Feng M Q,Shinozuka M.Cost-effective vision-based system for monitoring dynamic response of civil engineering structures[J].Structural Control and Health Monitoring,2010,17(8):918–936. doi: 10.1002/stc.360
    [6] Kohut P,Holak K,Uhl T,et al.Monitoring of a civil structure’s state based on noncontact measurements[J].Structural Health Monitoring,2013,12(5/6):411–429. doi: 10.1177/1475921713487397
    [7] Gentile C,Bernardini G.An interferometric radar for non-contact measurement of deflections on civil engineering structures:Laboratory and full-scale tests[J].Structure and Infrastructure Engineering,2010,6(5):521–534. doi: 10.1080/15732470903068557
    [8] Brown C J,Karuma R,Ashkenazi V,et al.Monitoring of structures using the global positioning system[J].Proceedings of the Institution of Civil Engineers−Structures and Buildings,1999,134(1):97–105. doi: 10.1680/istbu.1999.31257
    [9] Meng X,Dodson A H,Roberts G W.Detecting bridge dynamics with GPS and triaxial accelerometers[J].Engineering Structures,2007,29(11):3178–3184. doi: 10.1016/j.engstruct.2007.03.012
    [10] Xia H,de Roeck G,Zhang N,et al.Experimental analysis of a high-speed railway bridge under Thalys trains[J].Journal of Sound and Vibration,2003,268(1):103–113. doi: 10.1016/S0022-460X(03)00202-5
    [11] Pieraccini M,Parrini F,Fratini M,et al.Static and dynamic testing of bridges through microwave interferometry[J].NDT & E International,2007,40(3):208–214. doi: 10.1016/j.ndteint.2006.10.007
    [12] Nassif H H,Gindy M,Davis J.Comparison of laser Doppler vibrometer with contact sensors for monitoring bridge deflection and vibration[J].NDT & E International,2005,38(3):213–218. doi: 10.1016/j.ndteint.2004.06.012
    [13] Helfrick M N,Niezrecki C,Avitabile P,et al.3D digital image correlation methods for full-field vibration measurement[J].Mechanical Systems and Signal Processing,2011,25(3):917–927. doi: 10.1016/j.ymssp.2010.08.013
    [14] Morlier J,Salom P,Bos F.New image processing tools for structural dynamic monitoring[J].Key Engineering Materials,2007,347:239–244. doi: 10.4028/www.scientific.net/kem.347.239
    [15] Busca G,Cigada A,Mazzoleni P,et al.Cameras as displacement sensors to get the dynamic motion of a bridge[C]//Bridge Maintenance,Safety,Management,Resilience and Sustainability,Boca Raton:CRC Press,2012:2835–2841.
    [16] Wu Lijun,Casciati F,Casciati S.Dynamic testing of a laboratory model via vision-based sensing[J].Engineering Structures,2014,60:113–125. doi: 10.1016/j.engstruct.2013.12.002
    [17] 楚玺,向小菊,周志祥,等.欧拉运动放大算法在桥梁振动分析中的试验研究[J].公路交通科技,2019,36(6):41–47.

    Chu Xi,Xiang Xiaoju,Zhou Zhixiang,et al.Experimental study on Euler motion amplification algorithm in bridge vibration analysis[J].Journal of Highway and Transportation Research and Development,2019,36(6):41–47
    [18] Dong C Z,Ye X W,Jin T.Identification of structural dynamic characteristics based on machine vision technology[J].Measurement,2018,126:405–416. doi: 10.1016/j.measurement.2017.09.043
    [19] Feng Dongming,Feng M Q,Ozer E,et al.A vision-based sensor for noncontact structural displacement measurement[J].Sensors,2015,15(7):16557–16575. doi: 10.3390/s150716557
    [20] Warren C,Niezrecki C,Avitabile P,et al.Comparison of FRF measurements and mode shapes determined using optically image based,laser,and accelerometer measurements[J].Mechanical Systems and Signal Processing,2011,25(6):2191–2202. doi: 10.1016/j.ymssp.2011.01.018
    [21] Wang Weizhuo,Mottershead J E,Ihle A,et al.Finite element model updating from full-field vibration measurement using digital image correlation[J].Journal of Sound and Vibration,2011,330(8):1599–1620. doi: 10.1016/j.jsv.2010.10.036
    [22] Tregidgo L,Wang Z,Gursul I.Unsteady fluid-structure interactions of a pitching membrane wing[J].Aerospace Science and Technology,2013,28(1):79–90. doi: 10.1016/j.ast.2012.10.006
    [23] Molina−Viedma A J,Felipe−Sesé L,López−Alba E,et al.High frequency mode shapes characterisation using Digital Image Correlation and phase-based motion magnification[J].Mechanical Systems and Signal Processing,2018,102:245–261. doi: 10.1016/j.ymssp.2017.09.019
    [24] Molina−Viedma Á J,López−Alba E,Felipe−Sesé L,et al.Full-field modal analysis during base motion excitation using high-speed 3D digital image correlation[J].Measurement Science and Technology,2017,28(10):105402. doi: 10.1088/1361-6501/aa7d87
    [25] 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. doi: 10.1002/stc.2586
    [26] Beberniss T,Eason T,Spottswood S.High-speed 3D digital image correlation measurements of long duration random vibration;recent advancements and noted limitations[C]//Proceedings,International Conference on Noise and Vibration Engineering(ISMA),Leuven:Katholieke Universiteit Leuven,2012.
    [27] Hagara M,Trebuňa F,Huňady R,et al.Experimental identification of modal parameters of thin metal sheets by using of DIC[J].Procedia Engineering,2012,48:180–188. doi: 10.1016/j.proeng.2012.09.503
    [28] Trebuňa F,Hagara M.Experimental modal analysis performed by high-speed digital image correlation system[J].Measurement,2014,50:78–85. doi: 10.1016/j.measurement.2013.12.038
    [29] Baqersad J,Poozesh P,Niezrecki C,et al.Photogrammetry and optical methods in structural dynamics—A review[J].Mechanical Systems and Signal Processing,2017,86:17–34. doi: 10.1016/j.ymssp.2016.02.011
    [30] Busca G,Cigada A,Mazzoleni P,et al.Vibration monitoring of multiple bridge points by means of a unique vision-based measuring system[J].Experimental Mechanics,2014,54(2):255–271. doi: 10.1007/s11340-013-9784-8
    [31] Ribeiro D,Calçada R,Ferreira J,et al.Non-contact measurement of the dynamic displacement of railway bridges using an advanced video-based system[J].Engineering Structures,2014,75:164–180. doi: 10.1016/j.engstruct.2014.04.051
    [32] 郑佳艳,毛若愚,吴桐,等.基于欧拉微动放大的桥梁全息变形形态获取方法[J].科学技术与工程,2019,19(34):347–354. doi: 10.3969/j.issn.1671-1815.2019.34.051

    Zheng Jiayan,Mao Ruoyu,Wu Tong,et al.Eulerian-based micro motion amplification method to measure holographic deformation form of bridge[J].Science Technology and Engineering,2019,19(34):347–354 doi: 10.3969/j.issn.1671-1815.2019.34.051
    [33] Son K S,Jeon H S,Park J H,et al.Vibration displacement measurement technology for cylindrical structures using camera images[J].Nuclear Engineering and Technology,2015,47(4):488–499. doi: 10.1016/j.net.2015.01.011
    [34] Yoon H,Shin J,Spencer B F Jr.Structural displacement measurement using an unmanned aerial system[J].Computer-Aided Civil and Infrastructure Engineering,2018,33(3):183–192. doi: 10.1111/mice.12338
    [35] Dong Chuanzhi,Celik O,Catbas F N.Marker-free monitoring of the grandstand structures and modal identification using computer vision methods[J].Structural Health Monitoring,2019,18(5/6):1491–1509. doi: 10.1177/1475921718806895
    [36] 吴家勇,王平江,陈吉红,等.基于梯度重心法的线结构光中心亚像素提取方法[J].中国图象图形学报,2009,14(7):1354–1360. doi: 10.11834/jig.20090719

    Wu Jiayong,Wang Pingjiang,Chen Jihong,et al.Method of linear structured light sub-pixel center position extracting based on gradient barycenter[J].Journal of Image and Graphics,2009,14(7):1354–1360 doi: 10.11834/jig.20090719
    [37] 高世一,杨凯珍.变边限高斯拟合提取激光条纹中心线方法的研究[J].仪器仪表学报,2011,32(5):1132–1137. doi: 10.19650/j.cnki.cjsi.2011.05.028

    Gao Shiyi,Yang Kaizhen.Research on central position extraction of laser strip based on varied-boundary Gaussian fitting[J].Chinese Journal of Scientific Instrument,2011,32(5):1132–1137 doi: 10.19650/j.cnki.cjsi.2011.05.028
图(12)  /  表(2)

本文结构

    /

    返回文章
    返回