CN114741643B - 基于应变滑动空间相关模型的桥梁结构状态交叉诊断方法 - Google Patents
基于应变滑动空间相关模型的桥梁结构状态交叉诊断方法 Download PDFInfo
- Publication number
- CN114741643B CN114741643B CN202210377068.6A CN202210377068A CN114741643B CN 114741643 B CN114741643 B CN 114741643B CN 202210377068 A CN202210377068 A CN 202210377068A CN 114741643 B CN114741643 B CN 114741643B
- Authority
- CN
- China
- Prior art keywords
- window
- bridge structure
- space window
- diagnosed
- space
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/16—Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge
- G01B11/18—Measuring arrangements characterised by the use of optical techniques for measuring the deformation in a solid, e.g. optical strain gauge using photoelastic elements
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M99/00—Subject matter not provided for in other groups of this subclass
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
- Bridges Or Land Bridges (AREA)
Abstract
基于应变滑动空间相关模型的桥梁结构状态交叉诊断方法,属于土木工程结构运营安全智能监测领域。以主梁结构密集测点的应变监测数据为基础,将适当长度的主梁区段及该区段关于桥梁纵桥向对称轴对称的区域划分为一个空间窗;构建所划分空间窗内分布式应变的空间相关模型,利用该模型的应变残差,计算该空间窗对应主梁区域的桥梁结构状态诊断因子;采用交叉诊断方法实现该空间窗对应桥梁结构对称区域的状态诊断;从桥梁端部开始不断滑动空间窗,使空间窗遍历主梁全域,实现主梁所有区域的结构状态诊断。本发明解决了利用分布式应变实现桥梁结构状态的有效诊断和损伤区域的识别问题。
Description
技术领域
本发明属于土木工程结构运营安全智能监测领域,特别涉及一种基于应变滑动空间相关模型的桥梁结构状态交叉诊断方法。
背景技术
随着我国经济的快速发展,我国桥梁的设计与建设水平不断提升,已经创造了多项世界第一,当前全世界前10座最大跨度的悬索桥中国即占据6座,这些建设成果标志着我国从桥梁第一大国逐步迈向了第一强国。然而,随着服役年限的上升,我国桥梁的运营状况也将面临严峻的局面,特别是以往基于无损探伤和人工视察为主的检测和评估手段已经远远不能满足桥梁结构运营安全的要求。因此,如何采用有效手段实现桥梁结构状态的诊断已经成为亟待解决的关键问题。
随着桥梁健康监测技术的日渐成熟,桥梁结构健康监测系统已经成为桥梁结构状态诊断的有力工具。但是,目前桥梁状态诊断方法几乎都是利用点式传感设备获取的有限监测信息,而仅利用几处关键断面的监测数据通常难以有效识别结构发生的潜在损伤,因此,如何利用先进监测技术准确诊断复杂环境下桥梁结构的状态亟需被进一步的研究。
发明内容
本发明的目的是为了解决目前的桥梁状态诊断方法难以有效识别结构发生的潜在损伤,提供一种基于应变滑动空间相关模型的桥梁结构状态交叉诊断方法,利用分布式应变监测数据实现桥梁结构状态的有效诊断,并对损伤区域进行识别。
为实现上述目的,本发明采取的技术方案如下:
基于应变滑动空间相关模型的桥梁结构状态交叉诊断方法,所述方法包括以下步骤:
步骤一:将分布式传感光纤布设在桥梁主梁上,获取主梁结构密集测点的应变监测数据,将适当长度的主梁区段及该区段关于桥梁纵桥向对称轴对称的区域划分为一个空间窗;
步骤二:构建步骤一所划分空间窗内分布式应变的空间相关模型,以该模型的应变残差为基础,计算该空间窗对应主梁区域的桥梁结构状态诊断因子;
步骤三:以步骤一所划分空间窗内对称区域的桥梁结构状态诊断因子为基础,采用交叉诊断方法实现该桥梁结构对称区域的状态诊断;
步骤四:从桥梁端部开始不断滑动空间窗,并重复步骤一~步骤三,最终使空间窗遍历主梁全域,实现主梁所有区域的结构状态诊断。
本发明相对于现有技术的有益效果是:首先,利用主梁结构密集测点的应变监测数据,将适当长度的主梁区段及该区段关于桥梁纵桥向对称轴对称的区域划分为一个空间窗;其次,构建所划分空间窗内分布式应变的空间相关模型,利用该模型的应变残差,计算该空间窗对应主梁区域的桥梁结构状态诊断因子;在此基础上,采用交叉诊断方法实现该空间窗对应桥梁结构对称区域的状态诊断;最后,从桥梁端部开始不断滑动空间窗,使空间窗遍历主梁全域,实现主梁所有区域的结构状态诊断。与现有技术相比,解决了利用分布式应变实现桥梁结构状态的有效诊断和损伤区域的识别问题。
附图说明
图1为本发明的基于应变滑动空间相关模型的桥梁结构状态交叉诊断方法的流程图。
图2为实施例中桥梁结构关键点示意图。
图3为实施例中桥梁结构有限元模型示意图。
图4为实施例中空间窗的大小和位置示意图。
图5为实施例中桥梁结构左侧区域B点处的诊断结果图。
图6为实施例中桥梁结构左侧区域C点处的诊断结果图。
图7为实施例中桥梁结构右侧区域B′点处的诊断结果图。
图8为实施例中桥梁结构右侧区域C′点处的诊断结果图。
图9为实施例中桥梁结构左侧区域的状态诊断结果汇总图。
图10为实施例中桥梁结构右侧区域的状态诊断结果汇总图。
具体实施方式
下面结合实施例和附图对本发明的技术方案作进一步的说明,但并不局限于此,凡是对本发明技术方案进行修正或等同替换,而不脱离本发明技术方案的精神范围,均应涵盖在本发明的保护范围之中。
具体实施方式一:本实施方式记载的是基于应变滑动空间相关模型的桥梁结构状态交叉诊断方法,如图1所示,所述方法包括以下步骤:
步骤一:将分布式传感光纤布设在桥梁主梁上,获取主梁结构密集测点的应变监测数据,将适当长度的主梁区段及该区段关于桥梁纵桥向对称轴对称的区域划分为一个空间窗;
步骤二:构建步骤一所划分空间窗内分布式应变的空间相关模型,以该模型的应变残差为基础,计算该空间窗对应主梁区域的桥梁结构状态诊断因子;
步骤三:以步骤一所划分空间窗内对称区域的桥梁结构状态诊断因子为基础,采用交叉诊断方法实现该桥梁结构对称区域的状态诊断;
步骤四:从桥梁端部开始不断滑动空间窗,并重复步骤一~步骤三,最终使空间窗遍历主梁全域,实现主梁所有区域的结构状态诊断。
具体实施方式二:具体实施方式一所述的基于应变滑动空间相关模型的桥梁结构状态交叉诊断方法,步骤一中,所述空间窗的划分方法为:
(1)将靠近桥梁端部的空间窗边界称为外边界,远离桥梁端部的空间窗边界称为内边界;将空间窗内的两个相互对称的主梁区段分别称为左侧区域和右侧区域;如果空间窗还未开始滑动,则将外边界设置为桥梁端部;如果空间窗已经开始滑动,则将外边界设置为上一次空间窗所在位置的内边界;
(2)按下式计算空间窗内边界所有可能位置对应的空间窗长度判别指标,
式中,Cr为空间窗长度判别指标;为分布式应变空间相关模型参数;cov(·)为矩阵的协方差计算函数;/>为空间窗左侧区域第j个测点在i时刻的应变测量值;/>为i时刻空间窗右侧区域第j个测点的应变测量值;n为空间窗左侧区域的结构分布式应变测点总数;
(3)Cr最小时的内边界为最优内边界,由外边界和最优内边界完成该空间窗的划分。
具体实施方式三:具体实施方式一或二所述的基于应变滑动空间相关模型的桥梁结构状态交叉诊断方法,步骤二中,所述桥梁结构状态诊断因子的计算方法为:
(1)按下式建立窗内分布式应变的空间相关模型,
(2)按下式计算窗内分布式应变空间相关模型的残差ePU和eUP,
(3)将空间窗均匀划分为多个区域,各个区域的模型残差如下式所示,
式中,为空间窗左侧区域第v个待诊断区域的应变残差向量;/>为空间窗左侧区域内第n个测点在第s时刻的应变残差值,该测点属于第v个待诊断区域;/>为空间窗右侧区域第v个待诊断区域的应变残差向量;/>为空间窗右侧区域内第n个测点在第s时刻的应变残差值,该测点属于第v个待诊断区域;
(4)按下式计算左侧区域和右侧区域的桥梁结构状态诊断因子,
具体实施方式四:具体实施方式三所述的基于应变滑动空间相关模型的桥梁结构状态交叉诊断方法,步骤三中,所述的桥梁结构状态交叉诊断方法为:
(1)采集参考状态下桥梁结构分布式应变监测数据,计算相应的桥梁结构状态诊断因子,并将桥梁结构状态诊断因子从小到大排列,如下式所示,
(2)按下式计算空间窗左侧区域和右侧区域的桥梁结构状态诊断阈值,
式中,TYv,PU、TYv,UP分别为空间窗左侧区域第v个待诊断区域和空间窗右侧区域第v个待诊断区域的桥梁结构状态诊断阈值;分别为空间窗左、右侧区域第v个待诊断区域由小到大排列的桥梁结构状态诊断因子序列中取统计数据95%置信概率的分位数;
(3)按下式实现桥梁结构状态的诊断,
y′v,PU>TYv,PU (8)
y′v,UP>TYv,UP (9)
式中,y′v,PU、y′v,UP分别为待诊断状态下空间窗左侧区域第v个待诊断区域和空间窗右侧区域第v个待诊断区域的桥梁结构状态诊断因子;
若公式(8)成立,则记录空间窗左侧区域第v个待诊断区域异常,若公式(9)成立,则记录空间窗右侧区域第v个待诊断区域异常,
(4)按公式(8)和(9)可对所有待诊断区域进行诊断,如果空间窗左侧各个待诊断区域均有异常,而右侧仅有部分待诊断区域异常,则可判断桥梁结构损伤发生在空间窗右侧出现异常的待诊断区域;如果空间窗右侧各个待诊断区域均有异常,而左侧仅有部分待诊断区域异常,则可判断桥梁结构损伤发生在空间窗左侧出现异常的待诊断区域。
实施例1:
本实施例选取某实际桥梁结构的有限元模型作为算例,验证本发明提出的基于应变滑动空间相关模型的桥梁结构状态交叉诊断方法的有效性。该桥梁的结构关键点示意如图2所示,有限元模型如图3所示。采用本发明所提方法,从桥端逐步移动空间窗,每次移动时空间窗的大小和位置如图4所示。设置损伤发生在桥梁结构左侧区域B点,采用本发明所提方法诊断该桥梁结构状态,桥梁结构左侧区域B点和C点处的诊断结果分别如图5和图6所示,桥梁结构右侧区域B′点和C′点处的诊断结果分别如图7和图8所示,桥梁结构左侧区域的状态诊断结果汇总如图9所示,桥梁结构右侧区域的状态诊断结果汇总如图10所示。
由图5~10可以看出,桥梁结构左侧区域只有B点处的桥梁结构状态诊断因子超过阈值,而对桥梁结构右侧区域B点所在的空间窗,其绝大部分测点处的桥梁结构状态诊断因子都超过阈值。因此,根据本发明所提交叉诊断方法可以判断桥梁结构出现损伤,且损伤位于左侧区域B点,这一结果与实际状况一致,验证了本发明所提方法的有效性。
Claims (1)
1.基于应变滑动空间相关模型的桥梁结构状态交叉诊断方法,其特征在于:所述方法包括以下步骤:
步骤一:将分布式传感光纤布设在桥梁主梁上,获取主梁结构密集测点的应变监测数据,将适当长度的主梁区段及该区段关于桥梁纵桥向对称轴对称的区域划分为一个空间窗;所述空间窗的划分方法为:
(1)将靠近桥梁端部的空间窗边界称为外边界,远离桥梁端部的空间窗边界称为内边界;将空间窗内的两个相互对称的主梁区段分别称为左侧区域和右侧区域;如果空间窗还未开始滑动,则将外边界设置为桥梁端部;如果空间窗已经开始滑动,则将外边界设置为上一次空间窗所在位置的内边界;
(2)按下式计算空间窗内边界所有可能位置对应的空间窗长度判别指标,
式中,Cr为空间窗长度判别指标;为分布式应变空间相关模型参数;cov(·)为矩阵的协方差计算函数;为空间窗左侧区域第j个测点在i时刻的应变测量值;为i时刻空间窗右侧区域第j个测点的应变测量值;n为空间窗左侧区域的结构分布式应变测点总数;
(3)Cr最小时的内边界为最优内边界,由外边界和最优内边界完成该空间窗的划分;
步骤二:构建步骤一所划分空间窗内分布式应变的空间相关模型,以该模型的应变残差为基础,计算该空间窗对应主梁区域的桥梁结构状态诊断因子;所述桥梁结构状态诊断因子的计算方法为:
(1)按下式建立窗内分布式应变的空间相关模型,
(2)按下式计算窗内分布式应变空间相关模型的残差ePU和eUP,
(3)将空间窗均匀划分为多个区域,各个区域的模型残差如下式所示,
式中,为空间窗左侧区域第v个待诊断区域的应变残差向量;为空间窗左侧区域内第n个测点在第s时刻的应变残差值,该测点属于第v个待诊断区域;为空间窗右侧区域第v个待诊断区域的应变残差向量;为空间窗右侧区域内第n个测点在第s时刻的应变残差值,该测点属于第v个待诊断区域;
(4)按下式计算左侧区域和右侧区域的桥梁结构状态诊断因子,
步骤三:以步骤一所划分空间窗内对称区域的桥梁结构状态诊断因子为基础,采用交叉诊断方法实现该桥梁结构对称区域的状态诊断;所述的桥梁结构状态交叉诊断方法为:
(1)采集参考状态下桥梁结构分布式应变监测数据,计算相应的桥梁结构状态诊断因子,并将桥梁结构状态诊断因子从小到大排列,如下式所示,
(2)按下式计算空间窗左侧区域和右侧区域的桥梁结构状态诊断阈值,
式中,TYv,PU、TYv,UP分别为空间窗左侧区域第v个待诊断区域和空间窗右侧区域第v个待诊断区域的桥梁结构状态诊断阈值;分别为空间窗左、右侧区域第v个待诊断区域由小到大排列的桥梁结构状态诊断因子序列中取统计数据95%置信概率的分位数;
(3)按下式实现桥梁结构状态的诊断,
y′v,PU>TYv,PU (8)
y′v,UP>TYv,UP (9)
式中,y′v,PU、y′v,UP分别为待诊断状态下空间窗左侧区域第v个待诊断区域和空间窗右侧区域第v个待诊断区域的桥梁结构状态诊断因子;
若公式(8)成立,则记录空间窗左侧区域第v个待诊断区域异常,若公式(9)成立,则记录空间窗右侧区域第v个待诊断区域异常,
(4)按公式(8)和(9)可对所有待诊断区域进行诊断,如果空间窗左侧各个待诊断区域均有异常,而右侧仅有部分待诊断区域异常,则可判断桥梁结构损伤发生在空间窗右侧出现异常的待诊断区域;如果空间窗右侧各个待诊断区域均有异常,而左侧仅有部分待诊断区域异常,则可判断桥梁结构损伤发生在空间窗左侧出现异常的待诊断区域;
步骤四:从桥梁端部开始不断滑动空间窗,并重复步骤一~步骤三,最终使空间窗遍历主梁全域,实现主梁所有区域的结构状态诊断。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210377068.6A CN114741643B (zh) | 2022-04-11 | 2022-04-11 | 基于应变滑动空间相关模型的桥梁结构状态交叉诊断方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210377068.6A CN114741643B (zh) | 2022-04-11 | 2022-04-11 | 基于应变滑动空间相关模型的桥梁结构状态交叉诊断方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114741643A CN114741643A (zh) | 2022-07-12 |
CN114741643B true CN114741643B (zh) | 2023-03-24 |
Family
ID=82281352
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210377068.6A Active CN114741643B (zh) | 2022-04-11 | 2022-04-11 | 基于应变滑动空间相关模型的桥梁结构状态交叉诊断方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114741643B (zh) |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108009566B (zh) * | 2017-11-09 | 2020-04-28 | 华南理工大学 | 一种时空窗口下的改进型pca损伤检测方法 |
EP3644009A1 (de) * | 2018-10-24 | 2020-04-29 | Hahn-Schickard-Gesellschaft für angewandte Forschung e.V. | Dehnungsmessstruktur mit einem strukturierten träger |
CN110706213B (zh) * | 2019-09-23 | 2020-05-29 | 济南城建集团有限公司 | 基于应变响应累积分布函数差的桥梁集群结构损伤判别方法 |
CN114048776B (zh) * | 2021-11-10 | 2022-07-05 | 哈尔滨工业大学 | 基于损伤诊断指标非线性窄域特征的桥梁损伤识别方法 |
-
2022
- 2022-04-11 CN CN202210377068.6A patent/CN114741643B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN114741643A (zh) | 2022-07-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111737909B (zh) | 基于时空图卷积网络的结构健康监测数据异常识别方法 | |
CN106897543B (zh) | 模态柔度曲率矩阵范数的梁结构损伤识别方法 | |
CN108399277B (zh) | 一种基于温度与应变关联性的桥梁损伤识别方法 | |
CN105579922A (zh) | 信息处理装置以及分析方法 | |
US8386204B2 (en) | Method for the diagnosis of the EGR cooler efficiency in a diesel engine | |
US9969507B2 (en) | Method for performing diagnostics of a structure subject to loads and system for implementing said method | |
CN111354496B (zh) | 核电厂事故在线诊断及状态跟踪预测方法 | |
CN115950609B (zh) | 一种结合相关性分析和神经网络的桥梁挠度异常检测方法 | |
CN108022058A (zh) | 一种风力机状态可靠性评估方法 | |
CN111814870B (zh) | 一种基于卷积神经网络的cps模糊测试方法 | |
CN106932337A (zh) | 多光纤光栅钢筋腐蚀传感器系统的数据处理方法 | |
CN114741643B (zh) | 基于应变滑动空间相关模型的桥梁结构状态交叉诊断方法 | |
CN108009566A (zh) | 一种时空窗口下的改进型pca损伤检测方法 | |
CN110243598B (zh) | 列车轴承温度的处理方法、装置及存储介质 | |
CN115597901A (zh) | 一种桥梁伸缩缝损伤的监测方法 | |
CN111709944B (zh) | 一种基于数字化图像异常识别的桥隧结构裂缝监测方法 | |
CN109752383B (zh) | 一种基于多重交叉验证的桥梁损伤识别方法 | |
CN117665592A (zh) | 一种锂离子电池热失控故障诊断与预警方法 | |
CN114048776A (zh) | 基于损伤诊断指标非线性窄域特征的桥梁损伤识别方法 | |
KR102674766B1 (ko) | 인공지능을 이용한 밸브 고장진단 시스템 및 이를 이용한 밸브 고장진단 방법 | |
CN114692465B (zh) | 桥梁损伤位置的无损识别方法、存储介质及设备 | |
CN115456104A (zh) | 一种对传感器故障不敏感的桥梁拉索异常识别定位方法 | |
CN114693151A (zh) | 基于分布式应变空间相关性的桥梁桩基础状态诊断方法 | |
CN116226758A (zh) | 一种基于改进粒子滤波和卷积网络的雕刻机异常检测方法 | |
CN118013808B (zh) | 一种海上风力发电结构螺栓松动健康状态鉴别方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |