基于应变响应累积分布函数差的桥梁集群结构损伤判别方法与流程

本发明属于时变环境下桥梁结构损伤判别领域，特别涉及一种基于应变响应累积分布函数差的桥梁集群结构损伤判别方法。

背景技术：

由于结构健康监测技术可以持续可靠地提供桥梁结构响应及环境荷载信息、识别设计与施工中的结构缺陷、诊断与评估运营中的结构损伤，所以该技术已在桥梁工程中得以广泛应用。随着桥梁结构运营安全监测系统的集群化发展，区域路网或城市环线内的多座桥梁均建立了结构运营安全监测系统，即建立了桥梁集群结构运营安全监测系统。但是，现有方法均是利用单座桥梁结构响应监测数据与其历史参考状态进行对比来诊断损伤。而利用集群内多座桥梁结构响应的监测数据，在不使用历史参考状态的情况下，准确判别集群内全部桥梁结构损伤的方法还没有这方面的直接研究成果。

技术实现要素：

本发明是为了解决集群内有监测系统的相似桥梁结构损伤判别的问题，现提供基于应变响应累积分布函数差的桥梁集群结构损伤判别方法。

基于应变响应累积分布函数差的桥梁集群结构损伤判别方法，该方法包括以下步骤：

步骤一：采集健康状态下桥梁集群内有监测系统的全部相似桥梁所有测点的结构应变监测数据；

步骤二：采用基于数据相似度的监测应变响应聚类分析方法对结构应变监测数据的每一列进行聚类分析，获得所有测点按照数据相关度的聚类划分结果；

步骤三：根据步骤二的聚类划分结果，采集待诊断状态下，某一类别内所有测点的结构应变监测数据；

步骤四：根据步骤三获得的某一类别内所有测点的结构应变监测数据，利用基于桥梁结构应变监测数据集合重构的降噪方法，获得该类别减弱噪声干扰后的结构应变监测数据；

步骤五：在步骤四的基础上，在该类别内随机选择参考点，构建待诊断状态下该类别内所有测点的损伤判别因子及阈值；

步骤六：将该类别内所有测点的损伤判别因子与损伤判别阈值进行比较，当测点的损伤判别因子大于损伤判别阈值时，则该测点处发生了损伤，获得该类别内全部测点的损伤判别结果；

步骤七：重复步骤三至步骤六，获得集群内所有类别下的测点的损伤判别结果。

本发明所述的基于应变响应累积分布函数差的桥梁集群结构损伤判别方法，采用健康状态下的桥梁结构应变监测数据，进行聚类划分，找到所受环境因素影响相似的测点；利用待诊断状态下的桥梁结构应变监测数据，建立同一类别下的桥梁结构损伤判别因子及阈值；通过对不同测点的结构损伤判别因子与阈值进行比较，诊断桥梁在测点处是否产生结构损伤。

本发明所述方法利用相同监测时段下多座桥梁环境温度相似的特点，能够剔除环境温度对桥梁结构监测响应的影响；同时，本发明所述方法通过基于桥梁结构应变监测数据集合重构的降噪方法，能够减弱测试噪声对结构应变监测数据的干扰。

本发明适用于解决实际运营的桥梁集群结构损伤判别问题。本发明还能够直接应用于桥梁集群结构健康监测系统，实现对桥梁集群结构状态的在线实时诊断。

附图说明

图1为桥梁集群结构有限元模型图。

图2为环境温度与混凝土弹性模量的关系图。

图3为底板与顶板二个区域的温度时程图。

图4为减弱噪声干扰后结构应变监测数据相似度聚类分析结果。

图5为损伤状态下所有测点的损伤判别特征。

图6为损伤状态下所有测点的损伤判别因子。

图7为具体实施方式一所述的基于应变响应累积分布函数差的桥梁集群结构损伤判别方法的流程图。

具体实施方式

具体实施方式一：本实施方式所述的基于应变响应累积分布函数差的桥梁集群结构损伤判别方法，该方法包括以下步骤：

步骤一：采集健康状态下桥梁集群内有监测系统的全部相似桥梁所有测点的结构应变监测数据；

步骤二：采用基于数据相似度的监测应变响应聚类分析方法对结构应变监测数据的每一列进行聚类分析，获得所有测点按照数据相关度的聚类划分结果；

步骤三：根据步骤二的聚类划分结果，采集待诊断状态下，某一类别内所有测点的结构应变监测数据；

步骤四：根据步骤三获得的某一类别内所有测点的结构应变监测数据，利用基于桥梁结构应变监测数据集合重构的降噪方法，获得该类别减弱噪声干扰后的结构应变监测数据；

步骤五：在步骤四的基础上，在该类别内随机选择参考点，构建待诊断状态下该类别内所有测点的损伤判别因子及阈值；

步骤六：将该类别内所有测点的损伤判别因子与损伤判别阈值进行比较，当测点的损伤判别因子大于损伤判别阈值时，则该测点处发生了损伤，获得该类别内全部测点的损伤判别结果；

步骤七：重复步骤三至步骤六，获得集群内所有类别下的测点的损伤判别结果。

基于结构应变响应累积分布函数差的桥梁集群结构损伤判别方法，利用聚类分析，找到所受环境温度一致的测点，在健康状态下，即线弹性系统下，结构应变监测数据的概率分布可由温度的概率分布线性变换得到，因此聚类后可以建立多个具有相同概率分布的结构应变监测数据，在每一类别下，经过减弱噪声处理后的结构应变监测数据，以累计分布函数残差为损伤判别特征构建损伤判别因子，并与损伤判别阈值进行比较，从而判别桥梁结构损伤与否。

具体实施方式二：本实施方式是对具体实施方式一所述的基于应变响应累积分布函数差的桥梁集群结构损伤判别方法作进一步说明，本实施方式中，步骤二所述的获得所有测点按照数据相关度的聚类划分结果的方法为：

步骤二一：根据步骤一获得的健康状态下的结构应变监测数据集合x，其矩阵中每一列代表同一个测点，每一行代表同一时刻，矩阵x可表达式为：

x＝[ω1,ω2,…,ωi,…,ωn]m×n(1)

式中，ωi为第i个测点健康状态下的结构应变监测数据；1≤i≤n；n为测点个数；m为健康状态采样点个数。

步骤二二：根据健康状态下的结构应变监测数据集合x的任意列向量为ωi，即任意两个测点ωi和ωj的相关系数为ρij：

式中，σi²和σj²分别为ωi和ωj的方差；σij为ωi和ωj的协方差；1≤i≤n；1≤j≤n；i≠j；

步骤二三：根据不同测点的应变响应监测数据的相关系数ρij作为测点聚类划分标准，对健康状态下结构应变监测数据集合x进行聚类划分。

具体实施方式三：本实施方式是对具体实施方式一所述的基于应变响应累积分布函数差的桥梁集群结构损伤判别方法作进一步说明，本实施方式中，步骤四所述的获得该类别减弱噪声干扰后的结构应变监测数据的方法为：

步骤四一：假设待诊断状态下某一类别的结构应变监测数据集合为矩阵y，则矩阵y的构造矩阵定义为：

a＝y^ty(3)

式中，(·)^t表示对矩阵转置；y为待诊断状态下结构应变监测数据集合，u为待诊断状态采样点个数；v为某一类别的测点数；

步骤四二：对构造矩阵a进行奇异值分解：

a＝usvv^t(4)

式中，u和v分别为构造矩阵a的左奇异向量矩阵和右奇异向量矩阵；sv为矩阵a的奇异值按由大到小排列所构成的对角阵；

步骤四三：保留sv中前r个奇异值，得到重构矩阵a′：

式中，a′为构造矩阵a降噪后的重构矩阵；sr为重构矩阵a′保留的前r个奇异值构成的对角阵；

步骤四四：由待诊断状态下某一类别的结构应变监测数据集合y和重构矩阵a′，获得减弱噪声干扰后的结构应变监测数据的重构矩阵y′：

y′＝(y^t)^-*·a′(6)

式中，y′为待诊断状态下某一类别的结构应变监测数据集合y降噪后的重构矩阵；(·)^-*为对矩阵求伪逆。

具体实施方式四：本实施方式是对具体实施方式一所述的基于应变响应累积分布函数差的桥梁集群结构损伤判别方法作进一步说明，本实施方式中，步骤五所述的构建待诊断状态下该类别内所有测点的损伤判别因子及阈值的方法为：

步骤五一：设待诊断状态下去除噪声后某一类别的结构应变监测数据集合为y′，定义如下：

式中，u为待诊断状态采样点个数；v为某一类别的测点数；p为待诊断状态的某一采样点，1≤p≤u；q为某一类别的某一测点，1≤q≤v；y′pq为待诊断状态下聚类降噪后第q个测点的结构应变监测数据的第p个采样点；

步骤五二：则第q个测点在第p个时刻结构应变监测数据的累计分布函数值定义为：

f(y′pq)＝∑p(≤y′pq)(8)

式中，p(y′pq)为待诊断状态下聚类降噪后第q个测点的结构应变监测数据的第p个采样点的概率；

步骤五三：设选择该类别内第c个测点作为参考点，则损伤判别特征εpq定义为：

εpq＝f(y′pq)-f(y′pc)(9)

式中，c为从v个测点中随机选取的参考测点；1≤c≤v；

步骤五四：根据损伤判别特征εku，获得损伤判别因子hq：

则待诊断状态下聚类组内所有测点的损伤判别因子集合η定义为：

η＝{h1,h2,…,hq,…,hv}(11)

步骤五五：根据待诊断状态下聚类组内所有测点的损伤判别因子集合，定义该类别下的阈值α为：

α＝β·{η}95％(12)

式中，β为保证系数根据经验选取，建议选择不小于2；{·}95％取95％的分位数，向下取整。

采用下述试验来验证本发明的效果：

本试验是以图1所示的三座结构形式相似的三跨混凝土连续梁桥作为桥梁集群为例，图1中为桥梁集群结构有限元模型图。为了便于模拟结构损伤，采用局部范围内抗弯刚度(弹性模量)的缩减来模拟结构可能出现的裂缝损伤或混凝土强度降低等损伤现象，损伤位置在#2桥梁的第26号测点，损伤程度由0％递增到1％。为了测试噪声的干扰，试验针对8％的噪声水平进行处理。假设环境温度与混凝土材料弹性模量线性相关，其相关关系如图2所示；同时，假设三座桥上温度划分为顶板温度和底板温度2个区域，底板和顶板的温度时程曲线如图3所示。

本试验具体如下：

采集桥梁集群结构健康状态下1～12月内所有测点的结构应变监测数据，建立健康状态下桥梁集群数据模型；

利用基于数据相似度的监测应变响应聚类分析方法对减弱噪声干扰后的结构应变监测数据的每一列进行聚类分析，获得所有测点按照数据相关度的聚类划分结果，具体结果如图4所示；

采集桥梁集群结构待诊断状态下13～24月内所用测点的结构应变监测数据，(本试验模拟24个月的监测结果，其中1～12月为集群内全部桥梁结构健康状态，13～24月为集群内有#2桥梁产生结构损伤的状态，该桥梁结构损伤程度由0％递增到1％，噪声水平为8％)，建立待诊断状态下桥梁集群数据模型；

利用基于监测应变响应数据重构的降噪算法处理结构应变监测数据，获得减弱噪声干扰后的结构应变监测数据；

在同一类别内确定参考点，确定损伤判别特征如图5所示，构建该类别内所有测点的损伤判别因子集合及阈值如图6所示，通过与阈值进行比较，实现桥梁集群结构该类别内所有测点的损伤判别。