摘要
目前,绝大多数基于深度学习的结构损伤识别方法依靠深度神经网络自动提取结构的损伤敏感特征,并通过损伤状态之间特征的差异实现模式分类识别.然而,这些方法面临着损伤量化难度大的挑战,并且需要大量的模型训练数据.本文提出基于模型嵌入循环神经网络(Model-Embedding Recurrent Neural Network,MERNN)的损伤识别方法.首先,通过数据驱动的卷积神经网络(Convolutional Neural Network,CNN)建立荷载-响应之间的映射关系,然后,利用龙格库塔法改进传统的循环神经网络,建立基于循环神经网络架构的数值计算单元.最后,基于结构响应计算值与实测响应残差构成的损失函数与神经网络的自动微分机制来实现结构刚度参数的更新,进而实现结构损伤识别.数值模拟框架与实验室的3层剪切型框架的损伤识别结果表明,本文提出的方法能基于少量响应数据准确量化结构损伤.
结构健康监测系统在结构服役过程中积累了大量监测数据,基于监测数据的结构损伤识别和状态评估方法得到越来越多的关注. 深度学习方法作为一种有效的数据特征提取方法,在结构损伤识别领域取得了丰富的成果. 2017年,Lin
以有限元模型修正为代表的损伤识别方法具有有效量化损伤的能力. Zhu
现有基于时域响应的损伤量化方法大多需要外荷载的先验知识如荷载形式等,且完成荷载识别过程时间长,无法直接从测得的结构响应中捕捉损伤-荷载-响应之间的关联性,存在计算效率低、参数多的问题. 深度学习方法具有强大的非线性拟合和高维特征提取能力,可以直接从监测数据中捕捉荷载时间序列与响应时间序列的联系,实现不基于结构参数的荷载识别. 卷积神经网络的参数共享机制,对提取多维特征向量中的相似特征并减少操作中的参数非常有利,被广泛用于响应重构、荷载识别等时间序列建模中. Li
本文提出了一种基于模型嵌入循环神经网络MERNN的结构损伤识别方法. 利用卷积神经网络建立多种结构响应与荷载之间的联系,然后建立基于循环神经网络架构的数值计算单元. 构建基于计算响应与实测响应残差的损失函数,通过神经网络的误差反向传播和自动微分机制来进行基准模型参数的更新,通过更新前后0参数的变化来实现结构损伤识别. 数值模拟框架与实验室的3层剪切型框架的损伤识别结果表明,本文提出的方法能基于少量响应数据准确量化结构损伤.
1 基于MERNN的损伤识别方法
线性动力学系统的动力响应方程可表示为如下形式:
(1) |
(2) |
式中:、、分别为结构的质量矩阵、阻尼矩阵和刚度矩阵,瑞利阻尼通常可以表示为,其中和为瑞利阻尼系数;分别为结构的位移、速度和加速度响应;为结构上作用的荷载.
实际中,作用在结构上的外荷载往往为未知状态. 因此要按
1.1 基于卷积神经网络的动力荷载识别
为了识别未知的动力荷载的时程,本文采取传统的包含多层2D卷积层、池化层和全连接层的CNN架构来构建测得的结构响应与外荷载之间的映射关系,根据估计的荷载和荷载的真实值之间的差异计算荷载估计误差,利用误差更新网络的权重参数. 本文通过估计荷载和荷载真实值之间的均方误差来构建CNN荷载识别网络的损伤函数.
(3) |
1.2 基于龙格库塔法改进的循环神经网络单元
通过卷积神经网络获得结构上作用的未知荷载参数后,使用常见的显式积分求解方法如Newmark-法或四阶龙格库塔法对数值微分方程[
(4) |
式中:、、分别代表第n个、第n+个和第n+1个时间点;、分别表示第n个时间点与第n+1个时间点的结构响应;~为龙格库塔法的数值计算点间的斜率;h为步长.
在给定初始状态后,显式迭代运算可以得到结构在每个时刻的位移等响应.
为了嵌入结构响应计算方法,本文中选择使用四阶龙格库塔法代替循环神经网络的内部计算流程. LSTM单元是经典的循环神经网络单元,通过遗忘门、输入门和输出门控制状态向量和输出向量在网络中的流动,解决了梯度消失问题.
MERNN单元首先基于每个时刻的初位移和初速度向量按
(5) |
(6) |
式中:;分别为第i个时间点与第i+1个时间点上第k个子步的位移、速度和加速度;分别为第i-1个时间点的位移和速度;;为时间间隔.
计算出两个时间点内四个子步的位移、速度和加速度向量后,如
(7) |
式中:可表示为结构的位移、速度或加速度向量,;B为一系数矩阵,.
传统的长短时记忆网络或循环神经网络可以表示一种显式迭代计算过程,因此可以使用显式的四阶龙格库塔法即式(5)~

图1 基于循环神经网络架构的响应计算单元
Fig.1 Response calculation cell based on RNN architecture
其中MERNN层中的MERNN单元将CNN模型识别到的荷载作为输入,将结构响应作为输出. 在数值计算单元中,此时荷载项等效于LSTM单元中的单元输入,由结构的位移和速度组成的状态向量等效于LSTM单元中的记忆向量,改进后的循环神经网络可以用于结构响应的计算.
1.3 结构损伤识别方法
改进的MERNN层和常规的神经网络层一样,可以在内部设置可训练的参数,故首先将未损伤结构的质量矩阵、刚度矩阵和阻尼矩阵作为基准模型参数引入网络,在模型训练的过程中基于卷积神经网络预测出的荷载与结构参数计算出结构响应,与真实的结构响应进行对比,根据时域响应之间的差异对结构的刚度参数进行更新. 基于计算响应与真实响应之间的误差构建MERNN网络的损失函数,如
(8) |
(9) |
(10) |
(11) |
式中:和分别为实测响应和计算响应;分别为参数更新前刚度、参数更新后刚度;表示第i个参数的梯度;代表使用的学习率.
通过上述的方式即可在神经网络中建立可训练的结构局部刚度参数,利用误差反向传播算法对由于损伤引起的结构刚度变化进行量化,实现损伤识别,整体的损伤识别流程如

图2 整体架构图
Fig.2 Overall framework of proposed method
首先,基于结构实测响应结合卷积神经网络进行荷载识别,之后基于无损结构的刚度、质量、阻尼矩阵对MERNN进行初始化,其中部分结构刚度参数作为网络中的可训练参数,网络结合荷载估计值可以计算出结构响应,利用计算结构响应和实测结构响应的差异对可训练参数进行更新,通过更新前后刚度参数的对比即可实现结构损伤量化.
2 数值模拟
采用一个3层剪切型钢框架来验证提出的方法,每一层柱的高度为600 mm,柱的几何尺寸为4 mm×50 mm,横梁的长度为800 mm,横梁的几何尺寸为40 mm×50 mm,框架的钢材的弹性模量为190 GPa,密度为7 850 kg/
为了验证提出的方法对不同类型的荷载均有良好的效果,按照荷载形式的不同设置了两种不同的工况,具体的工况形式如下:
1)工况一(周期荷载工况):对底层柱的侧向刚度折减40%来模拟损伤,即有,在第3层梁柱节点处施加周期荷载,采样频率为100 Hz,采样时间取10 s,为生成网络训练用的数据集,通过改变周期荷载的幅值和频率来生成数据,荷载相关参数满足,. 使用第2层和第3层框架节点的位移响应和外荷载时程作为训练数据,共有450个训练样本.
2)工况二(冲击荷载工况):结构损伤状态与工况一相同,荷载施加位置相同,荷载的大小满足,冲击荷载的作用时间为第a秒到第a+0.05秒,采样频率为100 Hz,采样时间为10 s,通过改变冲击荷载的大小和作用时间生成数据,荷载相关参数满足[11,12,…,100],[0,0.1,0.2,0.3,04]. 共有500个训练样本,每一条时间序列包含1 000个数据点.
通过卷积神经网络可以充分捕捉荷载和结构响应之间的映射关系,对周期荷载和冲击荷载均具有良好的识别效果. 周期荷载均方估计误差为0.009,冲击荷载的均方估计误差为0.005.
将卷积神经网络识别的荷载和实测结构响应一起作为后续MERNN损伤识别网络的输入, 3层框架模型的柱的侧向刚度作为MERNN中可训练的参数,以无损情况下的3层框架基准模型的刚度、质量、阻尼矩阵作为网络的初始化参数,结合荷载估计值,损伤识别网络通过对比基准模型计算的响应和实测响应,对神经网络中可训练的柱侧向刚度参数进行更新.
图

(a) 工况一周期荷载识别结果

(b) 工况一参数更新前后的第3层节点位移响应对比

(c) 工况一结构损伤识别结果

(d) 工况二冲击荷载识别结果

(e) 工况二参数更新前后的第3层节点位移响应对比

(f) 工况二损伤识别结果

(g) 工况一中损失函数的变化情况
图3 3层剪切型框架荷载和损伤识别结果
Fig.3 Identification results of load and damage detection of three-story frame
由
3 3层框架试验
为了进一步验证提出方法的有效性,利用实验室的3层钢框架进行研究.如

图4 试验结构与损伤位置
Fig.4 Experiment structure and damage location
由于初始模型与试验模型存在一定程度的误差,直接进行基于时域数据的模型修正和损伤识别困难较大,故此处先基于采集到的加速度响应进行基于频域的模型修正,通过敲击框架第3层右侧梁柱节点,得到的加速度响应如

图5 试验中测得的加速度响应
Fig.5 Acceleration response sampled from test
模态阶数 | 实测频率/Hz | 修正前频率/Hz | 修正后频率/Hz | 误差/% |
---|---|---|---|---|
1 | 1.94 | 2.19 | 1.97 | 1.55 |
2 | 5.68 | 5.88 | 5.65 | 0.53 |
3 | 8.62 | 8.50 | 8.62 | 0.00 |
修正前各层柱的侧向刚度均为1.287 5 kN/m,修正后从上到下的侧向刚度分别为0.907 6 kN/m、1.281 7 kN/m和1.409 0 kN/m. 如
损伤工况设置在底层柱处,如
(12) |
式中:为等效宽度;为矩阵的特征值;与分别为等效前后的单元刚度矩阵.

图6 切口柱的等效侧向刚度
Fig.6 Equivalent lateral stiffness of column with cut
首先利用冲击荷载下的数值模型的位移响应加上部分实测数据训练CNN荷载识别网络,训练样本生成方式与第2节类似,本例中采样频率取30 Hz,采样时间取为20 s,网络的学习率取为1×1

(a) 荷载识别结果

(b) 参数更新前后结构响应对比图

(c) 损伤程度识别结果
图7 荷载及损伤识别结果
Fig.7 Result of load and damage identification
由于识别到的荷载和真实荷载的差异、框架模型简化带来的误差,模型参数更新后位移响应依旧与实测响应存在一定的差异,但参数更新后,实测响应与更新后位移响应的相位贴合程度提高,计算位移响应与实测位移响应的相似度从27.32%增长到97.1%. 损伤识别结果表明底层左柱的侧向刚度相对于经过频域模型修正后的侧向刚度基准值折减了8%,其他柱的刚度与基准值的差别均小于4%,与底层柱损伤均摊后的结果近似. 由于CNN网络识别的荷载与实际荷载之间的误差以及学习率等超参数未精细调整等原因,参数更新后模型的响应与实测响应并未完全对齐,因此在损伤识别中引入了部分误差,但该方法仍能有效识别结构的损伤位置和程度.
通过上述算例可知,利用卷积神经网络获得结构的荷载信息后,基于MERNN的损伤识别网络构成了一个基于时域数据的模型嵌入-数据驱动的有限元模型修正模块. 基于神经网络的误差反向传播机制和Tensorflow中的自动微分机制,在模型训练的过程中,会根据模型输出值和目标值显式计算时域响应偏差关于模型可训练参数的梯度,根据
4 结 论
本文提出基于MERNN网络的结构损伤识别方法,通过传统的卷积神经网络构建响应与外荷载之间的联系,进而建立模型嵌入循环神经网络进行结构局部刚度参数的更新,通过局部刚度参数的变化来实现损伤量化. 3层剪切型框架损伤识别算例显示,本文提出的方法对周期荷载和冲击荷载均有良好的识别效果,面对结构损伤程度较大的工况,有着良好的损伤识别效果. 本文得到的主要结论如下:
1)本文提出的方法在荷载识别步中通过网络学习结构响应与荷载之间的映射关系,无须明确的结构参数,且有良好的抗噪性能,可以解决传统荷载识别与损伤识别相互耦合的问题.
2)基于MERNN的损伤识别网络本质上建立了一种基于时域响应数据的局部刚度参数修正方法,相对于传统基于优化算法和响应灵敏度结合的有限元模型修正方法,该方法通过神经网络的误差反向传播机制可以自动计算误差关于可训练参数的梯度,通过梯度直接指示结构参数的优化方向,该方法在一定程度上可以避免计算结构复杂的响应灵敏度.
参考文献
LIN Y Z, NIE Z H, MA H W. Structural damage detection with automatic feature-extraction through deep learning[J]. Computer-Aided Civil and Infrastructure Engineering, 2017,32(12):1025-1046. [百度学术]
MODARRES C, ASTORGA N, DROGUETT E L, et al.Convolutional neural networks for automated damage recognition and damage type identification[J]. Structural Control and Health Monitoring,2018,25(10):e2230. [百度学术]
SHU J P, ZHANG C G, CHEN X Y, et al. Model-informed deep learning strategy with vision measurement for damage identification of truss structures[J]. Mechanical Systems and Signal Processing, 2023, 196: 110327. [百度学术]
TENG S, CHEN X D, CHEN G F, et al. Structural damage detection based on transfer learning strategy using digital twins of bridges[J]. Mechanical Systems and Signal Processing, 2023,191:110160. [百度学术]
李佐强, 翁顺, 朱宏平, 等. 基于领域自适应支持向量机的跨域损伤识别[J]. 建筑钢结构进展, 2023, 25(9): 91-100. [百度学术]
LI Z Q, WENG S, ZHU H P,et al.Cross-domain damage identification based on domain adaption-support vector machine[J].Progress in Steel Building Structures,2023,25(9):91-100.(in Chinese) [百度学术]
ZHU H P,MAO L,WENG S.A sensitivity-based structural damage identification method with unknown input excitation using transmissibility concept[J].Journal of Sound and Vibration,2014,333(26):7135-7150. [百度学术]
CHEN Z D,WENG S,YU H,et al.Bayesian-based method for the simultaneous identification of structural damage and moving force[J].Mechanical Systems and Signal Processing,2023,185:109742. [百度学术]
GUO X Y,FANG S E.Structural parameter identification using physics-informed neural networks[J].Measurement,2023,220:113334. [百度学术]
LI Y X,NI P,SUN L M,et al.A convolutional neural network-based full-field response reconstruction framework with multitype inputs and outputs[J].Structural Control and Health Monitoring,2022,29(7):1-23. [百度学术]
LEI X M,SUN L M,XIA Y.Lost data reconstruction for structural health monitoring using deep convolutional generative adversarial networks[J].Structural Health Monitoring,2021, 20(4):2069-2087. [百度学术]
ROSAFALCO L,MANZONI A,MARIANI S,et al.An autoencoder-based deep learning approach for load identification in structural dynamics[J].Sensors,2021,21(12):4207. [百度学术]