基于物理和数据双驱动的桥梁内部响应重构方法技术

本发明专利技术公开了一种基于物理和数据双驱动的桥梁内部响应重构方法，首先利用加速度传感器获得桥梁未知荷载下作用的加速度响应；其次基于频响函数，将物理逻辑嵌入到神经网络中；把物理方程以及相应的边界条件和初始条件作为惩罚项放入损失函数，并以此来限制可行性解的空间；然后训练神经网络模型，在将传感器得到的已知点加速度响应集输入到网络中预测未知点的加速度。本发明专利技术通过将直接求解控制方程的问题转换为损失函数的优化问题来找到方程的解，能有效解决桥梁內部响应难以测量、过分依赖测量数据等问题，并能提高桥梁内部响应预测的准确性和鲁棒性。测的准确性和鲁棒性。测的准确性和鲁棒性。

【技术实现步骤摘要】
基于物理和数据双驱动的桥梁内部响应重构方法


[0001]本专利技术涉及桥梁结构健康监测领域，具体地说是一种桥梁响应重构方法。

技术介绍

[0002]近十年来，结构内部响应重构引起了人们的广泛关注。对于工程结构的综合诊断和预测，需要准确估计结构内部关键部位的健康状态，以确定维修计划。由于结构内部点的响应不可达，即传感器无法测量到结构内部的响应。因此，发展一种在工作状态下对难以直接测得的结构内部响应进行预示的方法便尤为重要。在建筑结构中，虚拟传感器方法通过采集部分容易测量的物理量进行处理和计算，得到需要监测的其他物理量的值。根据虚拟传感器的思想，使用部分可以测得的结构位置的响应，来预示难测得位置的响应，从而解决目标位置响应无法在工作状态下直接测量的问题。
[0003]针对以上问题，主要有两种建模方法：基于物理/方程的建模和基于数据驱动的机器学习。基于模型的方法假设一个精确的系统模型的可用性，而数据驱动的方法基于机器学习。这两种方法都有一定的局限性，基于物理/方程的建模，基于预先已知的机理的知识对于需要监测的过程进行建模。然而，在实际应用中，即使有了对整个模型的知识，还是很难通过这些已有的知识建立一个完备的机理模型。其中的原因，一是因为通常情况下整个模型结构十分复杂，很难建立一个非常准确的机理模型，二是因为在工作过程中，会产生很多的不确定性，也使得机理模型很难用于解决实际的工程问题中。数据驱动的虚拟传感器是通过整理与处理收集的历史数据来进行分析与回归。由于其不需要对原系统结构具有完全的了解，而是需要大量的数据对模型进行训练，因此其效率较高，然而，由于系统的输入参数和输出变量的离散性，数据永远不可能是完整的，不管数据获取是否昂贵，数据永远无法填充状态空间。因此，系统和它产生的数据之间的鸿沟在数据驱动模型中根深蒂固。它们往往会在训练数据之外失败。因此缺乏推断力。导致实际响应预测的结果并不准确。

技术实现思路

[0004]本专利技术是为了解决上述现有技术存在的不足之处，提供一种基于物理和数据双驱动的桥梁内部响应重构方法，以期能有效解决桥梁內部响应难以测量、过分依赖测量数据等问题，并能提高桥梁内部响应预测的准确性和鲁棒性。
[0005]本专利技术为达到上述专利技术目的，采用如下技术方案：
[0006]本专利技术一种基于物理和数据双驱动的桥梁内部响应重构方法的特点是，包括以下步骤：
[0007]步骤1：在桥梁的m个节点上分别布置传感器，用于记录桥梁的加速度响应；
[0008]步骤2：利用有限元模型对桥梁进行建模，得到桥梁有限元模型，并对所述桥梁有限元模型施加随机荷载，得到m个外部节点的加速度响应矩阵及其频响函数矩阵H
k
＝[H1,H2,
…
H
i
,
…
H
m
]T
、n个桥梁内部的加速度响应矩阵
及其频响函数矩阵H
u
＝[H
′1,H
′2,
…
H
′
j
,
…
H
′
n
]T
，其中，u
i
为第i个外部节点处长度为w的加速度响应，u
′
j
为第j个内部节点处长度为w的加速度响应，H
i
为第m个外部节点处长度为w的频率响应，H
j
为第n个内部节点处长度为w的频率响应，令总加速度响应集令总频响函数集H＝{H
k
,H
u
}；
[0009]步骤3：构建卷积神经网络模型，依次包括：输入层、a层卷积层、a
‑
1层最大池化层、e层全连接层、输出层；且相邻两个卷积层之间添加有一层所述最大池化层；各全连接层之间采用relu激活函数；
[0010]将作为输入集输入所述卷积神经网络中，并依次经过交叉连接的a层卷积层和a
‑
1层最大池化层的卷积、池化处理后，由第a层卷积层输出维度为b
×
c
×
d的特征向量；其中b表示输入信号的通道数，c表示卷积产生的通道数，c表示卷积核的大小；
[0011]维度为b
×
c
×
d的特征向量集通过扁平处理后，得到维度为b
×
c的特征向量集；
[0012]维度为b
×
c的特征向量集依次经过e层全连接层的处理后，得到维度为w
×
n的内部节点响应预测向量记为
[0013]步骤4：构建差值函数；
[0014]步骤4.1：对中的每个加速度响应进行降维处理后，再进行傅里叶转换，并得到维度为m
×
w的频域上m个外部节点的加速度响应矩阵其中，u
i
(w)为第i个节点处长度为w的频域上的加速度响应；
[0015]步骤4.2：对进行傅里叶变换，得到维度为n
×
w的频域上n个内部响节点应预测向量其中，u
′
j
(w)为网络预测的第j个内部节点处长度为w的频域上的加速度响应；
[0016]步骤4.3：利用式(1)构建结构动力学中外部节点响应与内部节点响应的关系方程：
[0017][0018]式(1)中，T
aku
(ω)为传递矩阵，将内部点的加速度频响函数矩阵即H
u
与外部节点的加速度频响函数矩阵H
k
的伪逆相乘后，得到真实的加速度响应传递矩阵T
aku
(w)，从而利用式(2)计算差值f：
[0019][0020]步骤4.4：将与的伪逆相乘后，得到一个预测的加速度响应传递矩阵T
net
(w)＝[a1(w),a2(w),
…
,a
j
(w),
…
a
m
(w)]；其中，a
j
(w)为第j个外部节点处长度为n的加速度响应传递函数；
[0021]利用式(3)得到差值T：
[0022]T＝T
net
(w)
‑
T
aku
(ω)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(3)
[0023]步骤5：利用式(4)生成卷积神经网络的总损失L
all
：
[0024]L
all
＝L1+L2+L3+L4ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(4)
[0025]式(4)中，L1表示第一损失，并由式(5)得到，L2表示第二损失，并由式(6)得到，L3表示第三损失，并由式(7)得到，L4表示第四损失，并由式(8)得到；
[0026]L1＝(f
‑
f0)
T
(f本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种基于物理和数据双驱动的桥梁内部响应重构方法，其特征是，包括以下步骤：步骤1：在桥梁的m个节点上分别布置传感器，用于记录桥梁的加速度响应；步骤2：利用有限元模型对桥梁进行建模，得到桥梁有限元模型，并对所述桥梁有限元模型施加随机荷载，得到m个外部节点的加速度响应矩阵及其频响函数矩阵H
k
＝[H1,H2,
…
H
i
,
…
H
m
]
T
、n个桥梁内部的加速度响应矩阵及其频响函数矩阵H
u
＝[H
′1,H
′2,
…
H
′
j
,
…
H
′
n
]
T
，其中，u
i
为第i个外部节点处长度为w的加速度响应，u
′
j
为第j个内部节点处长度为w的加速度响应，H
i
为第m个外部节点处长度为w的频率响应，H
j
为第n个内部节点处长度为w的频率响应，令总加速度响应集令总频响函数集H＝{H
k
,H
u
}；步骤3：构建卷积神经网络模型，依次包括：输入层、a层卷积层、a
‑
1层最大池化层、e层全连接层、输出层；且相邻两个卷积层之间添加有一层所述最大池化层；各全连接层之间采用relu激活函数；将作为输入集输入所述卷积神经网络中，并依次经过交叉连接的a层卷积层和a
‑
1层最大池化层的卷积、池化处理后，由第a层卷积层输出维度为b
×
c
×
d的特征向量；其中b表示输入信号的通道数，c表示卷积产生的通道数，c表示卷积核的大小；维度为b
×
c
×
d的特征向量集通过扁平处理后，得到维度为b
×
c的特征向量集；维度为b
×
c的特征向量集依次经过e层全连接层的处理后，得到维度为w
×
n的内部节点响应预测向量记为步骤4：构建差值函数；步骤4.1：对中的每个加速度响应进行降维处理后，再进行傅里叶转换，并得到维度为m
×
w的频域上m个外部节点的加速度响应矩阵其中，u
i
(w)为第i个节点处长度为w的频域上的加速度响应；步骤4.2：对进行傅里叶变换，得到维度为n
×
w的频域上n个内部响节点应预测向量其中，u
′
j
(w)为网络预测的第j个内部节点处长度为w的频域上的加速度响应；步骤4.3：利用式(1)构建结构动力学中外部节点响应与内部节点响应的关系方程：式(1)中，T
aku
(ω)为传递矩阵，将内部点的加速度频响函数矩阵即H

【专利技术属性】
技术研发人员：贺文宇，陈文瀚，户东阳，李怡帆，李志东，王佐才，张静，
申请(专利权)人：合肥工业大学，
类型：发明
国别省市：