CN111967185A - 基于索力和位移分布相关性建模的斜拉桥状态评估方法 - Google Patents

基于索力和位移分布相关性建模的斜拉桥状态评估方法 Download PDF

Info

Publication number
CN111967185A
CN111967185A CN202010794104.XA CN202010794104A CN111967185A CN 111967185 A CN111967185 A CN 111967185A CN 202010794104 A CN202010794104 A CN 202010794104A CN 111967185 A CN111967185 A CN 111967185A
Authority
CN
China
Prior art keywords
cable
cable force
distribution
vertical displacement
main beam
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.)
Granted
Application number
CN202010794104.XA
Other languages
English (en)
Other versions
CN111967185B (zh
Inventor
李惠
徐阳
田亚迪
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Harbin Institute of Technology
Original Assignee
Harbin Institute of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Harbin Institute of Technology filed Critical Harbin Institute of Technology
Priority to CN202010794104.XA priority Critical patent/CN111967185B/zh
Publication of CN111967185A publication Critical patent/CN111967185A/zh
Application granted granted Critical
Publication of CN111967185B publication Critical patent/CN111967185B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/27Design optimisation, verification or simulation using machine learning, e.g. artificial intelligence, neural networks, support vector machines [SVM] or training a model
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/13Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Software Systems (AREA)
  • Medical Informatics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Artificial Intelligence (AREA)
  • Architecture (AREA)
  • Civil Engineering (AREA)
  • Structural Engineering (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)

Abstract

本发明涉及一种基于索力和位移分布相关性建模的斜拉桥状态评估方法,对拉索索力及主梁位移监测进行去趋势项处理,然后进行归一化,去除标准差以内的数据点,获得车致效应概率分布;搭建无监督图像变换模型,组装成分布矩阵作为网络输入;选取训练集数据输入至无监督图像变换模型进行训练,使用损失函数为根据网络特点所构造的综合函数,优化算法为ADAM算法;将测试集中的主梁竖向位移分布输入至经过训练的无监督图像变换模型中,计算预测与真实拉索索力分布之间的EMD距离,作为桥梁状态改变的评价指标,评估桥梁状态是否发生改变。本发明可以在桥梁结构***及荷载时空信息均未知的情况下,实现基于分布相关性建模的斜拉桥状态智能评估。

Description

基于索力和位移分布相关性建模的斜拉桥状态评估方法
技术领域
本发明涉及桥梁工程领域,具体涉及一种基于索力和位移分布相关性建模的斜拉桥状态评估方法。
背景技术
随着国民经济的快速发展,我国迎来了大跨桥梁的建设高潮。斜拉桥由于其跨度大、造型美观,格外受到设计选型青睐。斜拉桥的主要关键受力构件是拉索,由于拉索质量小、阻尼小、柔度大,在风雨及车辆交变载荷的作用下极易发生往复振动,引起疲劳损伤,进而有可能发生断裂。另外,斜拉桥在全寿命服役期内会不可避免地遭受环境侵蚀、材料老化等复杂因素的耦合作用,发生渐变损伤的萌生、发展和累积,导致结构产生抗力衰减,甚至在极端情况下发生灾难性事故。
为了保证斜拉桥的安全服役,桥梁管理部门会投入大量的财力在斜拉桥上安装结构健康监测***,对各种不同类型桥梁结构响应进行长期监测,并根据监测数据进行桥梁状态评估。结构健康监测***可以对不同桥梁结构响应(如拉索索力、主梁竖向位移等)进行监测,并在运营期内积累海量结构健康监测大数据。目前已有一些基于数据驱动的桥梁结构健康监测和状态评估方法。然而,这些方法一部分是针对有限元仿真或实验室条件下的模拟或测试数据,在实际大跨度斜拉桥上无法得到真正有效的应用,原因在于:一方面,实际大跨度斜拉桥结构十分复杂,数值模拟或实验室模型试验中都做了大量简化;另一方面,大跨度斜拉桥所处的现场环境非常复杂,结构响应数据不仅受温度作用影响明显,而且运营期间的实际车辆及风等动力荷载的大小很难准确获得。再者,基于传统力学分析的结构响应关系建模方法需要已知桥梁脉冲响应函数或频响函数,以及准确的荷载大小及位置信息。所以,传统力学分析方法在实际桥梁结构中也很难应用。近年来也有学者提出了一些基于单一类型监测数据的桥梁状态评估方法,但是此类方法的评估效果大都会高度依赖于监测数据质量。考虑到桥梁结构健康监测***的实际情况,传感器故障和数据异常是实际监测数据中必然要面对的情况。因此,基于单一类型监测数据的桥梁状态评估方法极易受到数据质量的约束:当某些通道的传感器数据丢包或异常时,评估结果往往不准确。以上这些因素均为基于实测结构健康监测数据的斜拉桥状态评估带来了巨大困难。
拉索作为斜拉桥的关键构件,将主梁自重及车辆荷载等传递到桥塔上,起到了极其重要的连接作用;主梁竖向位移则是在桥梁施工及运营期间直接衡量桥梁状态的重要指标。目前国内外学者所进行的桥梁健康监测数据挖掘研究通常是针对单一变量进行的,并且大都是从连续时程数据分析角度出发,对两种及以上变量簇之间的关系建模研究甚少,特别是从数据分布角度出发的研究几乎没有:即使在部分数据或传感器异常状态下,监测数据的部分时程信息虽然丢失,但是其整体分布信息仍然可以保持基本不变。因此,对拉索索力与主梁竖向位移这两种重要结构响应监测变量的分布相关性进行建模,更能反映斜拉桥的整体及局部状态。若这两种重要结构响应监测变量的分布相关性发生改变,则说明桥梁状态可能发生改变,并且对应不同的模式。因此,如何在桥梁结构***及荷载时空信息均未知的情况下,利用拉索索力与主梁竖向位移的监测数据,建立基于分布相关性建模的斜拉桥状态评估方法,提出斜拉桥状态改变的评价指标并挖掘其内在诱因,为大跨度斜拉桥由于状态改变是否需要进行维修提供桥梁管理和交通管制的决策依据,包括但不限于诸如桥梁健康监测***升级、传感器异常、拉索断丝、地震、强风超限振动、堵车、超载、船撞等异常事件。
发明内容
基于以上不足之处,本发明的目的是公开一种基于索力和位移分布相关性建模的斜拉桥状态评估方法,在桥梁结构***及荷载时空信息均未知的情况下,利用拉索索力与主梁竖向位移的监测数据,建立基于分布相关性建模的斜拉桥状态评估方法,提出斜拉桥状态改变的评价指标并挖掘其内在诱因,解决了实际大跨度斜拉桥面对异常事件情况下的维护管理决策问题,为实际大跨度斜拉桥由于状态改变而需要进行维修决策和交通管制提供依据,可以在桥梁结构***及荷载时空信息均未知的情况下,实现基于分布相关性建模的斜拉桥状态智能评估。
一种基于索力和位移分布相关性建模的斜拉桥状态评估方法,包括:
步骤一、对拉索索力及主梁位移进行预处理:首先进行去趋势项处理,根据数据特点去除温度影响,获得车致拉索索力与主梁竖向位移响应;然后对数据进行归一化处理,去除0.3倍标准差以内的数据点,获得车致拉索索力与主梁竖向位移的概率分布;
步骤二、搭建无监督图像变换模型,选取所有可用通道的拉索索力与主梁竖向位移数据分布,按行组装成拉索索力与主梁竖向位移分布矩阵,分别作为网络的两个输入:主梁竖向位移分布矩阵与拉索索力分布矩阵;
步骤三、从经过步骤一处理的数据中选取一部分构成训练集,输入至无监督图像变换模型中进行训练,训练过程中使用的损失函数为根据网络特点所构造的综合函数,优化算法为ADAM算法,并使用默认的超参数设置;
步骤四、将测试集中的主梁竖向位移分布输入至经过训练的无监督图像变换模型中,得到预测的拉索索力分布作为输出,并计算预测与真实拉索索力分布之间的EMD距离(Earth Mover’s Distance),并以预测与真实拉索索力分布之间的EMD距离为桥梁状态改变的评价指标,评估桥梁状态是否发生改变及相应的作用模式。
本发明还具有如下特征:
1、如上的步骤一具体包括:
步骤一一、对拉索索力和主梁竖向位移时程数据采用相同的预处理流程,以30分钟为时间窗,将该时间窗内的所有数据点线性归一化到[0,1]范围,然后画出该段数据的频率柱状图,并以频率峰值位置处为中心,左右各取0.03,并挑选出落在此段频率范围内的原始数据极值点;
步骤一二、对步骤一一选取出的极值点进行中值滤波、线性插值和光滑处理,光滑窗口取3600,得到温度趋势项,再从原始数据中减去温度趋势项,作为经过预处理的车致竖向位移响应和拉索索力;
步骤一三、对各个通道的车致竖向位移响应和拉索索力数据分别进行归一化处理,使各通道数值保持在[-1,1]范围内,并去除0.3倍标准差以内的噪声数据点,获得车致拉索索力与主梁竖向位移的概率分布。
2、如上的步骤二在无监督图像变换模型时,选取所有可用的27个主梁竖向位移与139个拉索索力通道,将各通道数据分布按行整合成矩阵后作为网络的两个输入:主梁竖向位移分布矩阵与拉索索力分布矩阵。
搭建无监督图像变换模型主要分为两大模块,步骤二具体包括:
步骤二一、搭建变分自编码器模块:对于主梁竖向位移,以主梁竖向位移概率分布矩阵为输入,编码器由4个卷积层、4个残差连接层和2个共享卷积层组成,输入经过编码器后成为隐式特征向量然后再进入解码器,解码器由2个共享卷积层、4个残差连接层和4个卷积层组成;类似的,对于拉索索力,以拉索索力概率分布矩阵为输入,编码器由4个卷积层、4个残差连接层和2个共享卷积层组成,输入经过编码器后成为隐式特征向量然后再进入解码器,解码器由2个共享卷积层、4个残差连接层和4个卷积层组成。变分自编码器中的所有共享卷积层和隐式特征向量,对于主梁竖向位移和拉索索力共享参数。
步骤二二、搭建生成对抗网络模块:对于主梁竖向位移,以变分自编码器的解码器作为生成器,判别器由8个卷积层组成;类似的,对于拉索索力,以变分自编码器的解码器作为生成器,判别器由8个卷积层组成。同步骤二一中所述,主梁竖向位移和拉索索力的生成器(即变分自编码器中的解码器)存在参数共享。
3、如上的步骤三中,主梁竖向位移和拉索索力的原始监测年限为2006-2015年,采样频率为3Hz。选取2006-2007年的数据作为训练集,其余数据用作测试,来验证所建立模型的普适性。主梁竖向位移和拉索索力均采用每天凌晨3点至8点的监测数据,时间窗口为3小时,滑动步长为10min,因此每天5个小时的监测数据可以产生13个训练样本。最终,训练集和测试集分别包括2986和6236个样本。
步骤三中,模型训练所采用的误差函数具体为:
Figure BDA0002624845010000041
式中,L为损失函数,下标1代表主梁竖向位移,下标2代表拉索索力;VAE代表变分自编器模块(由编码器E和解码器G组成),GAN代表生成对抗网络模块(由生成器G和判别器D组成);CC代表循环约束项,由于存在参数和隐式特征向量共享,会存在以下情况:由主梁竖向位移的编码器E1和拉索索力的解码器G2生成拉索索力,以及由拉索索力的编码器E2和主梁竖向位移的解码器G1生成主梁竖向位移;min和max分别代表取最小值和最大值运算。
优化算法为ADAM算法,模型更新公式具体为:
Figure BDA0002624845010000042
Figure BDA0002624845010000043
式中,t代表第t个时间步,ωt代表网络模型参数,
Figure BDA0002624845010000044
Figure BDA0002624845010000045
分别代表第t个时间步的误差一阶梯度和二阶梯度,mt和vt分别为考虑梯度一阶、二阶矩的中间更新变量,β1和β2分别为对应的一阶、二阶衰减因子,η0为初始学习率,ε为分母非零控制参数。超参数设置默认值为β1=0.9,β2=0.99,ε=10-60=10-3
4、如上的步骤四中,在模型训练完毕后,以车致拉索索力与主梁竖向位移的概率分布作为输入,可以得到预测拉索索力的概率分布,然后计算预测与真实拉索索力概率分布之间的EMD距离(Earth Mover’s Distance),具体计算方式为:
Figure BDA0002624845010000046
式中,EMD表示EMD距离,P和Q分别代表模型预测与真实监测数据得到的拉索索力概率分布,m和n分别代表分布P和Q的分段区间数量,dij和fij分别代表由第i个P分布区段到第j个Q分布区段的位移变化值和概率数值变化值,min代表取最小值运算。
本发明的优点及有益效果如下:
(1)本发明针对大跨度斜拉桥斜拉索索力与主梁竖向位移的力学相关性建模难题,实现了在桥梁结构***及荷载时空信息均未知的情况下,建立了由主梁竖向位移和斜拉索索力概率分布的无监督图像变换模型;
(2)本发明提出预测与真实拉索索力概率分布之间的EMD距离作为桥梁状态改变的指标,能够体现结构响应监测变量概率分布相关性的变化,即说明桥梁状态的改变;
(3)本发明对监测数据的温度趋势项做了去趋势处理,因此对温度变化不敏感,普遍适用于大跨桥梁结构主梁竖向位移和拉索索力的全年监测数据;
(4)本发明的整个关系建模过程为数据驱动,显著降低了传统力学分析中由于模型假定不合理、不准确而引起的误差;
(5)本发明满足大跨度斜拉桥索力的在线预测和实时数据处理需求,即在模型训练完毕后,不必进行训练集更新,直接将采集到的主梁竖向位移和拉索索力监测数据作为输入,即可得到预测索力分布;
(6)本发明提高了大跨度斜拉桥拉索索力与主梁竖向位移力学相关性建模的自动化和智能化程度以及准确性和鲁棒性,为大跨度斜拉桥服役状态的自主智能评估提供了技术支撑。
(7)本发明可以诊断出由于桥梁健康监测***升级、传感器异常、拉索断丝、地震、强风超限振动、堵车、超载、船撞等异常事件而造成的大跨度斜拉桥结构状态改变,从而为桥梁交通管理和巡检维修决策提供科学依据,大大提高桥梁管理的智能化程度。
附图说明
图1为本发明的一个实施例的流程图;
图2为本发明步骤一中数据预处理去趋势项的示意图;
图3为本发明步骤二中无监督图像变换模型的示意图;
图4为本发明基于无监督图像变换模型获得的部分索力分布预测效果图;
图5为本发明基于无监督图像变换模型获得的部分索力分布EMD距离变化图。
具体实施方式
实施例1
本实施方式为基于索力和位移分布相关性建模的斜拉桥状态评估方法,如图1所示,包括:
步骤一、对拉索索力及主梁位移进行预处理:首先进行去趋势项处理,根据数据特点去除温度影响,获得车致拉索索力与主梁竖向位移响应,去趋势项过程如图2所示;然后对数据进行归一化处理,去除0.3倍标准差以内的数据点,获得车致拉索索力与主梁竖向位移的概率分布;
步骤二、搭建无监督图像变换模型,模型示意图如图3所示,选取所有可用通道的拉索索力与主梁竖向位移数据分布,按行组装成拉索索力与主梁竖向位移分布矩阵,分别作为网络的两个输入(主梁竖向位移分布矩阵与拉索索力分布矩阵);
步骤三、从经过步骤一处理的数据中选取一部分构成训练集,输入至无监督图像变换模型中进行训练,训练过程中使用的损失函数为根据网络特点所构造的综合函数,优化算法为ADAM算法,并使用默认的超参数设置;
步骤四、将测试集中的主梁竖向位移分布输入至经过训练的无监督图像变换模型中,得到预测的拉索索力分布作为输出,部分索力分布预测效果如图4所示,并计算预测与真实拉索索力分布之间的EMD距离(Earth Mover’s Distance),并以预测与真实拉索索力分布之间的EMD距离为桥梁状态改变的评价指标,评估桥梁状态是否发生改变及相应的作用模式,部分索力分布EMD距离变化如图5所示。
实施例2
本实施方式与具体实施方式一不同的是:
步骤一具体包括:
步骤一一、对拉索索力和主梁竖向位移时程数据采用相同的预处理流程,以30分钟为时间窗,将该时间窗内的所有数据点线性归一化到[0,1]范围,然后画出该段数据的频率柱状图,并以频率峰值位置处为中心,左右各取0.03,并挑选出落在此段频率范围内的原始数据极值点;
步骤一二、对步骤一一选取出的极值点进行中值滤波、线性插值和光滑处理,光滑窗口取3600,得到温度趋势项,再从原始数据中减去温度趋势项,作为经过预处理的车致竖向位移响应和拉索索力;
步骤一三、对各个通道的车致竖向位移响应和拉索索力数据分别进行归一化处理,使各通道数值保持在[-1,1]范围内,并去除0.3倍标准差以内的噪声数据点,获得车致拉索索力与主梁竖向位移的概率分布。
其它步骤与具体实施方式一相同。
实施例3
本实施方式与具体实施方式一或二不同的是:
搭建无监督图像变换模型主要分为两大模块,步骤二具体包括:
步骤二一、搭建变分自编码器模块:对于主梁竖向位移,以主梁竖向位移概率分布矩阵为输入,编码器由4个卷积层、4个残差连接层和2个共享卷积层组成,输入经过编码器后成为隐式特征向量然后再进入解码器,解码器由2个共享卷积层、4个残差连接层和4个卷积层组成;类似的,对于拉索索力,以拉索索力概率分布矩阵为输入,编码器由4个卷积层、4个残差连接层和2个共享卷积层组成,输入经过编码器后成为隐式特征向量然后再进入解码器,解码器由2个共享卷积层、4个残差连接层和4个卷积层组成。变分自编码器中的所有共享卷积层和隐式特征向量,对于主梁竖向位移和拉索索力共享参数。
步骤二二、搭建生成对抗网络模块:对于主梁竖向位移,以变分自编码器的解码器作为生成器,判别器由8个卷积层组成;类似的,对于拉索索力,以变分自编码器的解码器作为生成器,判别器由8个卷积层组成。同步骤二一中所述,主梁竖向位移和拉索索力的生成器(即变分自编码器中的解码器)存在参数共享。
其它步骤与具体实施方式一至二相同。
实施例4
本实施方式与具体实施方式一至三之一不同的是:
步骤三中,主梁竖向位移和拉索索力的原始监测年限为2006-2015年,采样频率为3Hz。选取2006-2007年的数据作为训练集,其余数据用作测试,来验证所建立模型的普适性。主梁竖向位移和拉索索力均采用每天凌晨3点至8点的监测数据,时间窗口为3小时,滑动步长为10min,因此每天5个小时的监测数据可以产生13个训练样本。最终,训练集和测试集分别包括2986和6236个样本。
步骤三中,模型训练所采用的误差函数具体为:
Figure BDA0002624845010000071
式中,L为损失函数,下标1代表主梁竖向位移,下标2代表拉索索力;VAE代表变分自编器模块(由编码器E和解码器G组成),GAN代表生成对抗网络模块(由生成器G和判别器D组成);CC代表循环约束项,由于存在参数和隐式特征向量共享,会存在以下情况:由主梁竖向位移的编码器E1和拉索索力的解码器G2生成拉索索力,以及由拉索索力的编码器E2和主梁竖向位移的解码器G1生成主梁竖向位移;min和max分别代表取最小值和最大值运算。
优化算法为ADAM算法,模型更新公式具体为:
Figure BDA0002624845010000081
Figure BDA0002624845010000082
式中,t代表第t个时间步,ωt代表网络模型参数,
Figure BDA0002624845010000083
Figure BDA0002624845010000084
分别代表第t个时间步的误差一阶梯度和二阶梯度,mt和vt分别为考虑梯度一阶、二阶矩的中间更新变量,β1和β2分别为对应的一阶、二阶衰减因子,η0为初始学习率,ε为分母非零控制参数。超参数设置默认值为β1=0.9,β2=0.99,ε=10-60=10-3
其它步骤与具体实施方式一至三相同。
具体实施方式五:
本实施方式与具体实施方式一至四之一不同的是:
步骤四中,在模型训练完毕后,以车致拉索索力与主梁竖向位移的概率分布作为输入,可以得到预测拉索索力的概率分布,然后计算预测与真实拉索索力概率分布之间的EMD距离(Earth Mover’s Distance),具体计算方式为:
Figure BDA0002624845010000085
式中,EMD表示EMD距离,P和Q分别代表模型预测与真实监测数据得到的拉索索力概率分布,m和n分别代表分布P和Q的分段区间数量,dij和fij分别代表由第i个P分布区段到第j个Q分布区段的位移变化值和概率数值变化值,min代表取最小值运算。
其它步骤与具体实施方式一至四相同。
本发明针对实际大跨度斜拉桥面对异常事件情况下的维护管理决策问题,通过对斜拉桥关键构件响应(即主梁竖向位移和拉索索力)进行数据挖掘,建立主梁竖向位移和拉索索力之间力学相关性的深度学习代理模型,并通过观测代理模型评价指标的变化来评估大跨度斜拉桥的状态变化,为实际大跨度斜拉桥由于状态改变而需要进行维修决策和交通管制提供依据,包括但不限于诸如桥梁健康监测***升级、传感器异常、拉索断丝、地震、强风超限振动、堵车、超载、船撞等异常事件。

Claims (7)

1.一种基于索力和位移分布相关性建模的斜拉桥状态评估方法,其特征在于,步骤如下:
步骤一、对拉索索力及主梁位移进行预处理:首先进行去趋势项处理,根据数据特点去除温度影响,获得车致拉索索力与主梁竖向位移响应;然后对数据进行归一化处理,去除0.3倍标准差以内的数据点,获得车致拉索索力与主梁竖向位移的概率分布;
步骤二、搭建无监督图像变换模型,选取所有可用通道的拉索索力与主梁竖向位移数据分布,按行组装成拉索索力与主梁竖向位移分布矩阵,分别作为网络的两个输入:主梁竖向位移分布矩阵、拉索索力分布矩阵;
步骤三、从经过步骤一处理的数据中选取一部分构成训练集,输入至无监督图像变换模型中进行训练,训练过程中使用的损失函数为根据网络特点所构造的综合函数,优化算法为ADAM算法,并使用默认的超参数设置;
步骤四、将测试集中的主梁竖向位移分布输入至经过训练的无监督图像变换模型中,得到预测的拉索索力分布作为输出,并计算预测与真实拉索索力分布之间的EMD距离,并以预测与真实拉索索力分布之间的EMD距离为桥梁状态改变的评价指标,评估桥梁状态是否发生改变及相应的作用模式。
2.根据权利要求1所述的基于索力和位移分布相关性建模的斜拉桥状态评估方法,其特征在于,步骤一具体包括:
步骤一一、对拉索索力和主梁竖向位移时程数据采用相同的预处理流程,以30分钟为时间窗,将该时间窗内的所有数据点线性归一化到[0,1]范围,然后画出该段数据的频率柱状图,并以频率峰值位置处为中心,左右各取0.03,并挑选出落在此段频率范围内的原始数据极值点;
步骤一二、对选取出的原始数据极值点进行中值滤波、线性插值和光滑处理,光滑窗口取3600,得到温度趋势项,再从原始数据中减去温度趋势项,作为经过预处理的车致竖向位移响应和拉索索力;
步骤一三、对各个通道的车致竖向位移响应和拉索索力数据分别进行归一化处理,使各通道数值保持在[-1,1]范围内,并去除0.3倍标准差以内的噪声数据点,获得车致拉索索力与主梁竖向位移的概率分布。
3.根据权利要求1或2所述的基于索力和位移分布相关性建模的斜拉桥状态评估方法,其特征在于,步骤二在无监督图像变换模型时,选取所有可用的27个主梁竖向位移与139个拉索索力通道,将各通道数据分布按行整合成矩阵后作为网络的两个输入:主梁竖向位移分布矩阵与拉索索力分布矩阵。
4.根据权利要求3所述的基于索力和位移分布相关性建模的斜拉桥状态评估方法,其特征在于,步骤二具体包括:
步骤二一、搭建变分自编码器模块:对于主梁竖向位移,以主梁竖向位移概率分布矩阵为输入,编码器由4个卷积层、4个残差连接层和2个共享卷积层组成,输入经过编码器后成为隐式特征向量然后再进入解码器,解码器由2个共享卷积层、4个残差连接层和4个卷积层组成;类似的,对于拉索索力,以拉索索力概率分布矩阵为输入,编码器由4个卷积层、4个残差连接层和2个共享卷积层组成,输入经过编码器后成为隐式特征向量然后再进入解码器,解码器由2个共享卷积层、4个残差连接层和4个卷积层组成,变分自编码器中的所有共享卷积层和隐式特征向量,对于主梁竖向位移和拉索索力共享参数;
步骤二二、搭建生成对抗网络模块:对于主梁竖向位移,以变分自编码器的解码器作为生成器,判别器由8个卷积层组成;对于拉索索力,以变分自编码器的解码器作为生成器,判别器由8个卷积层组成,主梁竖向位移和拉索索力的生成器参数共享。
5.根据权利要求1或2所述的基于索力和位移分布相关性建模的斜拉桥状态评估方法,其特征在于,步骤三中,主梁竖向位移和拉索索力的原始监测年限为2006-2015年,采样频率为3Hz,选取2006-2007年的数据作为训练集,其余数据用作测试,来验证所建立模型的普适性,主梁竖向位移和拉索索力均采用每天凌晨3点至8点的监测数据,时间窗口为3小时,滑动步长为10分钟,每天5个小时的监测数据产生13个训练样本,最终,训练集和测试集分别包括2986和6236个样本,模型训练所采用的误差函数具体为:
Figure FDA0002624843000000021
式中,L为损失函数,下标1代表主梁竖向位移,下标2代表拉索索力;
VAE代表变分自编器模块,由编码器E和解码器G组成;
GAN代表生成对抗网络模块,由生成器G和判别器D组成;
CC代表循环约束项,由于存在参数和隐式特征向量共享,会存在以下情况:由主梁竖向位移的编码器E1和拉索索力的解码器G2生成拉索索力,以及由拉索索力的编码器E2和主梁竖向位移的解码器G1生成主梁竖向位移;
min和max分别代表取最小值和最大值运算。
6.根据权利要求1或2所述的基于索力和位移分布相关性建模的斜拉桥状态评估方法,其特征在于,步骤三中,优化算法为ADAM算法,模型更新公式具体为:
Figure FDA0002624843000000031
式中,t代表第t个时间步,ωt代表网络模型参数,
Figure FDA0002624843000000034
Figure FDA0002624843000000032
分别代表第t个时间步的误差一阶梯度和二阶梯度,mt和vt分别为考虑梯度一阶、二阶矩的中间更新变量,β1和β2分别为对应的一阶、二阶衰减因子,η0为初始学习率,ε为分母非零控制参数,超参数设置默认值为β1=0.9,β2=0.99,ε=10-60=10-3
7.根据权利要求1或2所述的基于索力和位移分布相关性建模的斜拉桥状态评估方法,其特征在于,步骤四中,在模型训练完毕后,以车致拉索索力与主梁竖向位移的概率分布作为输入,得到预测拉索索力的概率分布,然后计算预测与真实拉索索力概率分布之间的EMD距离,具体计算方式为:
Figure FDA0002624843000000033
式中,
EMD表示EMD距离,
P和Q分别代表模型预测与真实监测数据得到的拉索索力概率分布,
m和n分别代表分布P和Q的分段区间数量,
dij和fij分别代表由第i个P分布区段到第j个Q分布区段的位移变化值和概率数值变化值,min代表取最小值运算。
CN202010794104.XA 2020-08-10 2020-08-10 基于索力和位移分布相关性建模的斜拉桥状态评估方法 Active CN111967185B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010794104.XA CN111967185B (zh) 2020-08-10 2020-08-10 基于索力和位移分布相关性建模的斜拉桥状态评估方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010794104.XA CN111967185B (zh) 2020-08-10 2020-08-10 基于索力和位移分布相关性建模的斜拉桥状态评估方法

Publications (2)

Publication Number Publication Date
CN111967185A true CN111967185A (zh) 2020-11-20
CN111967185B CN111967185B (zh) 2022-08-02

Family

ID=73364975

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010794104.XA Active CN111967185B (zh) 2020-08-10 2020-08-10 基于索力和位移分布相关性建模的斜拉桥状态评估方法

Country Status (1)

Country Link
CN (1) CN111967185B (zh)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113570121A (zh) * 2021-07-09 2021-10-29 桂林理工大学 一种基于卷积神经网络的斜拉索索力趋势预测方法
CN113865683A (zh) * 2021-12-01 2021-12-31 苏州博宇鑫交通科技有限公司 基于机器视觉的城市高架桥梁超重超载动态预警方法
CN114004467A (zh) * 2021-10-15 2022-02-01 河南工业大学 一种基于监测数据的预制装配式桥梁结构性能分析方法
CN114280061A (zh) * 2021-12-27 2022-04-05 交通运输部公路科学研究所 一种斜拉桥索梁锚固区技术状况的观测方法及监视窗
CN114896844A (zh) * 2022-05-11 2022-08-12 昆明理工大学 一种拱桥扣背索索力数据处理方法、***和存储介质
CN115048998A (zh) * 2022-06-13 2022-09-13 大连理工大学 一种基于监测数据的斜拉桥群索索力异常识别定位方法
CN116881819A (zh) * 2023-09-07 2023-10-13 成都理工大学 一种基于孤立森林的斜拉索工作状态监测方法
CN117521188A (zh) * 2023-07-31 2024-02-06 中铁建大桥工程局集团第四工程有限公司 一种三主桁钢桁拱桥高程线形调控方法以及***

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2444787A1 (en) * 2010-10-22 2012-04-25 Vilnius Gediminas Technical University Method and device for bridge state evaluation using dynamic method
CN105303462A (zh) * 2015-11-03 2016-02-03 国家电网公司 基于emd距离的燃煤发电机组班组煤耗对比分析方法
CN108225811A (zh) * 2018-01-30 2018-06-29 四川天佑慧桥科技有限公司 一种基于车辆荷载的桥梁结构安全评估***
CN110570329A (zh) * 2019-09-04 2019-12-13 四川智通路桥工程技术有限责任公司 一种桥梁施工全过程控制及可视化预警***
CN111062080A (zh) * 2019-12-20 2020-04-24 哈尔滨工业大学 斜拉索索力与主梁竖向位移时空相关深度学习建模方法
CN111310273A (zh) * 2020-03-16 2020-06-19 河北省交通规划设计院 一种基于多源数据的全桥结构安全状态监测方法及***

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2444787A1 (en) * 2010-10-22 2012-04-25 Vilnius Gediminas Technical University Method and device for bridge state evaluation using dynamic method
CN105303462A (zh) * 2015-11-03 2016-02-03 国家电网公司 基于emd距离的燃煤发电机组班组煤耗对比分析方法
CN108225811A (zh) * 2018-01-30 2018-06-29 四川天佑慧桥科技有限公司 一种基于车辆荷载的桥梁结构安全评估***
CN110570329A (zh) * 2019-09-04 2019-12-13 四川智通路桥工程技术有限责任公司 一种桥梁施工全过程控制及可视化预警***
CN111062080A (zh) * 2019-12-20 2020-04-24 哈尔滨工业大学 斜拉索索力与主梁竖向位移时空相关深度学习建模方法
CN111310273A (zh) * 2020-03-16 2020-06-19 河北省交通规划设计院 一种基于多源数据的全桥结构安全状态监测方法及***

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
YONG SHI 等: "Large-Power ZVZCS Full-Bridge Three-Level DC–DC Converter With Wide Operation Range and Its Application in Sapphire Crystal Furnace Power Supply", 《IEEE JOURNAL OF EMERGING AND SELECTED TOPICS IN POWER ELECTRONICS》 *
施洲 等: "大型公铁两用斜拉桥健康监测数据的时域峰值统计及频域融合处理", 《中国铁道科学》 *

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113570121A (zh) * 2021-07-09 2021-10-29 桂林理工大学 一种基于卷积神经网络的斜拉索索力趋势预测方法
CN114004467A (zh) * 2021-10-15 2022-02-01 河南工业大学 一种基于监测数据的预制装配式桥梁结构性能分析方法
CN113865683A (zh) * 2021-12-01 2021-12-31 苏州博宇鑫交通科技有限公司 基于机器视觉的城市高架桥梁超重超载动态预警方法
CN114280061A (zh) * 2021-12-27 2022-04-05 交通运输部公路科学研究所 一种斜拉桥索梁锚固区技术状况的观测方法及监视窗
CN114280061B (zh) * 2021-12-27 2022-10-14 交通运输部公路科学研究所 一种斜拉桥索梁锚固区技术状况的观测方法及监视窗
CN114896844A (zh) * 2022-05-11 2022-08-12 昆明理工大学 一种拱桥扣背索索力数据处理方法、***和存储介质
CN114896844B (zh) * 2022-05-11 2023-06-23 昆明理工大学 一种拱桥扣背索索力数据处理方法、***和存储介质
CN115048998A (zh) * 2022-06-13 2022-09-13 大连理工大学 一种基于监测数据的斜拉桥群索索力异常识别定位方法
CN117521188A (zh) * 2023-07-31 2024-02-06 中铁建大桥工程局集团第四工程有限公司 一种三主桁钢桁拱桥高程线形调控方法以及***
CN116881819A (zh) * 2023-09-07 2023-10-13 成都理工大学 一种基于孤立森林的斜拉索工作状态监测方法
CN116881819B (zh) * 2023-09-07 2023-11-14 成都理工大学 一种基于孤立森林的斜拉索工作状态监测方法

Also Published As

Publication number Publication date
CN111967185B (zh) 2022-08-02

Similar Documents

Publication Publication Date Title
CN111967185B (zh) 基于索力和位移分布相关性建模的斜拉桥状态评估方法
CN111274737A (zh) 一种机械设备剩余使用寿命预测方法及***
CN108873859B (zh) 基于改进关联规则的桥式抓斗卸船机故障预测模型方法
CN108444662B (zh) 一种基于日温度效应的桥梁损伤在线监测方法
CN110457786B (zh) 基于深度置信网络的卸船机关联规则故障预测模型方法
CN111691679B (zh) 基于数字孪生的预应力钢结构智能张拉方法
Yue et al. Mechanics-Guided optimization of an LSTM network for Real-Time modeling of Temperature-Induced deflection of a Cable-Stayed bridge
CN113900381B (zh) 一种基于物联网的钢结构远程健康监测平台及应用方法
CN114282398A (zh) 一种基于大数据的桥梁健康监测***及方法
CN116150897A (zh) 一种基于数字孪生的机床主轴性能测评方法及***
CN111062080B (zh) 斜拉索索力与主梁竖向位移时空相关深度学习建模方法
CN114186475B (zh) 一种基于Attention-MGRU的回转支承寿命预测方法
CN114548375B (zh) 基于双向长短期记忆神经网络的斜拉桥主梁动挠度监测方法
CN116341272A (zh) 一种用于配网工程数字化的施工安全风险管控***
CN116451514A (zh) 桥梁结构状态评估方法、装置以及电子设备
Wang et al. Digital prediction model of temperature-induced deflection for cable-stayed bridges based on learning of response-only data
CN116933152B (zh) 一种基于多维emd-pso-lstm神经网络的海浪信息预测方法及***
CN116663126A (zh) 一种基于通道注意力BiLSTM模型的桥梁温度效应预测方法
CN115217152B (zh) 一种沉管隧道管节接头张合变形的预测方法及装置
CN116049958A (zh) 历史建筑结构监测数据异常诊断与修复***
CN114880767B (zh) 一种基于注意力机制Dense-GRU网络的航空发动机剩余使用寿命预测方法
CN114399024B (zh) 油气浓度大数据智能检测***
Desjardins Structural health assessment of bridges by long-term vibration monitoring and automated operational modal analysis
Chun et al. An application of data fusion technology in structural health monitoring and damage identification
CN113673010A (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