CN110319990B - 基于倾角仪优化布置的桥梁动挠度监测方法 - Google Patents

基于倾角仪优化布置的桥梁动挠度监测方法 Download PDF

Info

Publication number
CN110319990B
CN110319990B CN201910461900.9A CN201910461900A CN110319990B CN 110319990 B CN110319990 B CN 110319990B CN 201910461900 A CN201910461900 A CN 201910461900A CN 110319990 B CN110319990 B CN 110319990B
Authority
CN
China
Prior art keywords
bridge
information entropy
displacement
optimal
sensors
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
Application number
CN201910461900.9A
Other languages
English (en)
Other versions
CN110319990A (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
China Railway Design Corp
Original Assignee
Harbin Institute of Technology
China Railway Design Corp
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, China Railway Design Corp filed Critical Harbin Institute of Technology
Priority to CN201910461900.9A priority Critical patent/CN110319990B/zh
Publication of CN110319990A publication Critical patent/CN110319990A/zh
Application granted granted Critical
Publication of CN110319990B publication Critical patent/CN110319990B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M5/00Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings
    • G01M5/0008Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings of bridges
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M5/00Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings
    • G01M5/0041Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings by determining deflection or stress
    • G01M5/005Investigating the elasticity of structures, e.g. deflection of bridges or air-craft wings by determining deflection or stress by means of external apparatus, e.g. test benches or portable test systems
    • 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
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/06Power analysis or power optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Architecture (AREA)
  • Structural Engineering (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Civil Engineering (AREA)
  • Testing Of Devices, Machine Parts, Or Other Structures Thereof (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明提出了一种基于倾角仪优化布置的桥梁动挠度监测方法,包括:建立倾角‑动挠度转化方程,在桥梁有限元模型中模拟桥面过车,提取出倾角仪布设候选节点的倾角时程数据与桥梁关键断面的位移时程数据作为遗传算法优化训练的数据集,优化过程中采用信息熵作为适应度函数,从而得到固定数量下的倾角仪最优布设位置以及相对应的信息熵。以预测的有限元模型中桥梁关键断面时程位移的相对误差等于5%为界限确定出临界信息熵,在倾角仪最优布置下的各传感器数量所对应的信息熵中找出小于等于临界信息熵的传感器数量的最小值即为最优传感器数量,最优传感器数量对应的最优布设位置即为倾角仪的布设位置。本方法能够对桥梁的动挠度进行准确的预测。

Description

基于倾角仪优化布置的桥梁动挠度监测方法
技术领域
本发明涉及桥梁工程健康监测领域,具体涉及一种基于倾角仪优化布置的桥梁动挠度监测方法。
背景技术
桥梁在现代交通运输中有着极其重要的战略地位。近些年来,越来越多的高速公路、铁路桥梁建成并服役,极大地提高了我国的交通运输效率、促进了经济的高速增长。
桥梁动挠度监测是桥梁健康监测系统的重要组成部分,是评定桥梁承载能力、车辆行驶安全性的一项重要指标,但是一直没有通用的监测方法,如基于拉线式位移计测量动挠度的方法需要把吊锤固定于桥下,对于桥下有水的桥梁或者桥面远离桥下地面的桥梁的测试受到了限制;激光测量桥梁动挠度则因为需要一个静态的参考点而只能测量离桥头比较近的桥梁;GPS测量桥梁动挠度的测试精度最多只能达到厘米级别,不能够满足桥梁动挠度的测试精度要求。因此,如何测试横跨大江河等复杂地质条件下的桥梁动挠度是个亟待解决的难题。开展复杂地质条件下及复杂结构的桥梁动挠度监测方法的研究具有重要的现实意义。
近年来随着桥梁健康监测技术的发展,倾角仪被越来越广泛地应用在监测桥梁动挠度的桥梁健康监测系统中。然而,其倾角-动挠度转化原理大多是基于桥梁倾角的曲线积分,此方法对于大型复杂桥梁结构往往不能够适用。此外倾角仪的布设位置与数量往往由工程师凭借其工程经验来完成布设,这种方法存在着一些不可避免的缺点:一是可能造成传感器布设的过度冗余,从而造成一定程度的浪费;二是对于大型复杂结构桥梁,即使布置较多的传感器所预测的桥梁动挠度也不一定能够达到精度要求。因此若仅只依靠工程师的经验来决定倾角仪的布设位置与数量不仅可能造成不必要的人力、物力、财力的浪费,而且在预测桥梁动挠度、判断桥梁是否安全问题中存在着巨大的安全隐患。
发明内容
为了解决现有的倾角-动挠度监测方法不能够很好地适用于大型复杂结构桥梁的动挠度监测、因而在预测桥梁动挠度、判断桥梁是否安全问题中存在着巨大的安全隐患,而倾角仪的数量和位置往往按照工程师的经验来选择和布置,可能造成传感器布设的过度冗余,从而造成一定程度的浪费的问题,本发明提出的一种基于倾角仪优化布置的桥梁动挠度监测方法。
本发明所采用的技术如下:一种基于倾角仪优化布置的桥梁动挠度监测方法,步骤如下:
步骤一、建立倾角仪所在桥梁节点的倾角来预测桥梁关键断面位移的结构动力学方程;
步骤二、在桥梁有限元模型中创建节点动力荷载来模拟桥面过车,提取出有限元模型中倾角仪布设候选节点的倾角时程数据与桥梁关键断面的位移时程数据分别作为训练集和验证集;
步骤三、采用遗传算法进行倾角仪布设位置的优化,优化过程中采用桥梁关键断面预测位移误差的信息熵函数作为遗传算法中的适应度函数,从而得到固定数量下的倾角仪最优布设位置以及相对应的信息熵。
步骤四、以预测的有限元模型中桥梁关键断面时程位移的相对误差等于5%为界限确定出临界信息熵,在倾角仪最优布置下的各传感器数量所对应的信息熵中找出小于等于临界信息熵的传感器数量的最小值即为最优传感器数量,最优传感器数量对应的最优布设位置即为倾角仪的布设位置。
本发明还具有如下技术特征:
1.如上所述的步骤一中结构动力学方程建立步骤如下:
步骤一一、对于梁式结构,根据结构动力学中的振型叠加法建立运动方程:
vN=ΦNqm (1)
其中vN∈RN×1代表竖向位移矩阵,ΦN∈RN×m代表竖向位移振型矩阵,qm∈Rm×1代表竖向位移矩阵的广义坐标向量矩阵,N代表布设倾角仪的个数,m代表所用阵型向量的阶数,在竖直平面内,ΦN是关于桥梁纵桥向长度坐标x的函数矩阵,qm是关于时间t的函数矩阵。
步骤一二、将公式(1)两边同时对x求偏导得:
Figure BDA0002078309700000021
其中θN∈RN×1是纵桥向竖直平面内的倾角矩阵,ΨN∈RN×m代表纵桥向竖直平面内的倾角振型矩阵,采用穆尔-彭罗斯广义逆来求解ΨN的伪逆。
当ΨN为列满秩矩阵时:
qm=(ΨN TΨN)-1ΨN TθN (3)
当ΨN为行满秩矩阵时:
qm=ΨN TNΨN T)-1θN (4)
由于qm仅仅为时间t的函数,因此可以用来预测同一时刻桥梁其他P个断面的位移即:
vP=ΦPqm (5)
其中,其中ΦP∈RP×m,P代表沿桥梁关键断面测点的P个位置。
2.如上所述的步骤二,步骤二具体包括:
步骤二一、在桥梁有限元模型中创建节点动力荷载来模拟桥面过车,节点动力荷载的大小与载重车的重量相一致,在桥面上的移动速度与载重车的过车速度相一致。
步骤二二、将有限元模型运行分析,提取出有限元模型中倾角仪布设候选节点的倾角时程数据与桥梁关键断面的位移时程数据分别作为训练集和验证集。
3.如上所述的步骤三,步骤三具体包括:
步骤三一、遗传算法训练过程中采用信息熵作为适应度函数,其理论依据在于:对于服从均值为零,方差为σ2的随机变量s,其概率密度函数为p(s),则其信息熵为:
Figure BDA0002078309700000031
设预测桥梁关键断面位移的误差矩阵ε=[ε1,···,εmP]∈RH×P(m=1,2,···,P)中每个误差向量εm服从均值为零,方差为
Figure BDA0002078309700000032
的高斯分布,其中P代表桥梁关键断面的个数,此时模型的信息熵为:
Figure BDA0002078309700000033
其中
Figure BDA0002078309700000041
为误差向量εm方差估计值;
步骤三二、遗传算法训练过程中对某一特定数量的倾角仪寻找出最优布设位置的理论依据:通过引入信息熵,将特定数量下倾角仪的最优布设位置问题转化为使预测的P个关键断面位移的误差的信息熵最小的问题,优化目标函数表示为:
Figure BDA0002078309700000042
步骤三三、运用遗传算法对倾角传感器进行优化布设的过程如下:
(1)初始化:随机生成N个传感器布设方案作为N个个体,作为第一代群体。
(2)适应度计算:根据每个个体的传感器布设方案,基于有限元模型中倾角仪布设候选节点的时程倾角数据,通过步骤一中建立的倾角转换位移的结构动力学方程来预测桥梁关键断面的时程位移,然后基于预测的时程位移的误差来计算信息熵,以信息熵的倒数作为每个个体的适应度,即信息熵越小,个体的适应度越大。
(3)选择:基于每个个体的适应度,以概率的方式决定选择能保留基因至下一代的个体,适应度越大的个体,基因遗传到下一代的概率就越大。
(4)交叉:在每一代个体之间随机进行基因交换,即个体之间的传感器位置存在一定程度的概率互换。
(5)变异:每个个体的一个或几个传感器布设位置有一定概率的随机变动。
经过上述一轮操作后进化得到下一代群体。
(6)终止条件判断:当进化代数达到最大进化代数时,以该代群体中的最大适应度个体作为最优解输出。
4.如上所述的步骤四,步骤四具体包括:
步骤四一、以预测的有限元模型中桥梁关键断面时程位移的相对误差等于5%为界限确定出临界信息熵的步骤:
εH×P=0.05×vH×P (9)
其中vH×P为P个桥梁关键断面的实际时程位移矩阵,它的每一列代表一个关键断面的时程位移数据,其中H为每个节点时程位移数据的长度;
Figure BDA0002078309700000043
模型的临界信息熵为:
Figure BDA0002078309700000051
其中ΘT代表模型的临界信息熵;
步骤四二、确定出临界信息熵后,在倾角仪最优布置下的各传感器数量所对应的信息熵中找出小于等于临界信息熵的传感器数量的最小值即为最优传感器数量,最优传感器数量对应的最优布设位置即为倾角仪的布设位置。
本发明的有益效果为:本发明提出了能够适用于所有桥梁结构的倾角-动挠度转化方程,以及在满足桥梁动挠度预测精度的前提下,在桥上布设最佳数量与最佳位置的倾角仪传感器,从而保证能够准确地预测出桥梁关键断面的动挠度。极大地避免了不必要的人力、物力、财力的浪费,极大地消除了在预测桥梁动挠度、判断桥梁是否安全问题中所存在的巨大的安全隐患。本发明还能满足桥梁健康监测实时预警子系统的需求,直接实时地对桥梁关键断面的预测挠度是否超限进行判断识别。本发明提高了桥梁健康监测实时预警子系统智能识别的自动化、智能化、准确性和鲁棒性,为桥梁健康监测动挠度在线实时预警子系统的建立提供了解决方案。
附图说明
图1为三跨连续刚构桥的有限元模型及其关键断面示意图;
图2为三跨连续钢构桥的累积有效竖向模态质量参与率图;
图3为不同数量的倾角仪最优布设位置下的信息熵变化曲线图;
图4为本发明算法辨识有限元模型中关键断面时程位移与理论时程位移对比图:(a)左边跨L/4断面;(b)左边跨L/2断面;(c)左边跨3L/4断面;(d)中跨L/4断面;(e)中跨L/2断面;(f)中跨3L/4断面;(g)右边跨L/4断面;(h)右边跨L/2断面;(i)右边跨3L/4断面;
图5为遗传算法优化出的7个倾角仪的最佳布设位置示意图;
图6为倾角仪传感器所测的原始信号图;
图7为小波去噪处理后的倾角信号图;
图8为本发明算法辨识实验室某三跨连续钢构模型桥关键断面时程位移与实测时程位移对比图:(a)左边跨L/4断面;(b)左边跨L/2断面;(c)左边跨3L/4断面;(d)中跨L/4断面;(e)中跨L/2断面;(f)中跨3L/4断面;(g)右边跨L/4断面;(h)右边跨L/2断面;(i)右边跨3L/4断面;
具体实施方式
下面结合附图举例对本发明做进一步的说明:
实施例1
一种基于倾角仪优化布置的桥梁动挠度监测方法,步骤如下:
步骤一、建立倾角仪所在桥梁节点的倾角来预测桥梁关键断面位移的结构动力学方程;
结构动力学方程建立步骤如下:
步骤一一、对于梁式结构,根据结构动力学中的振型叠加法建立运动方程:
vN=ΦNqm (1)
其中vN∈RN×1代表竖向位移矩阵,ΦN∈RN×m代表竖向位移振型矩阵,qm∈Rm×1代表竖向位移矩阵的广义坐标向量矩阵,N代表布设倾角仪的个数,m代表所用阵型向量的阶数,在竖直平面内,ΦN是关于桥梁纵桥向长度坐标x的函数矩阵,qm是关于时间t的函数矩阵。
步骤一二、将公式(1)两边同时对x求偏导得:
Figure BDA0002078309700000061
其中θN∈RN×1是纵桥向竖直平面内的倾角矩阵,ΨN∈RN×m代表纵桥向竖直平面内的倾角振型矩阵,采用穆尔-彭罗斯广义逆来求解ΨN的伪逆。
当ΨN为列满秩矩阵时:
qm=(ΨN TΨN)-1ΨN TθN (3)
当ΨN为行满秩矩阵时:
qm=ΨN TNΨN T)-1θN (4)
由于qm仅仅为时间t的函数,因此可以用来预测同一时刻桥梁其他P个断面的位移即:
vP=ΦPqm (5)
其中ΦP∈RP×m,P代表沿桥梁关键断面测点的P个位置。
步骤二、在桥梁有限元模型中创建节点动力荷载来模拟桥面过车,节点动力荷载的大小与载重车的重量相一致,在桥面上的移动速度与载重车的过车速度相一致,将有限元模型运行分析,提取出有限元模型中倾角仪布设候选节点的倾角时程数据与桥梁关键断面的位移时程数据分别作为训练集和验证集。
步骤三、采用遗传算法进行倾角仪布设位置的优化,优化过程中采用桥梁关键断面预测位移误差的信息熵函数作为遗传算法中的适应度函数,从而得到固定数量下的倾角仪最优布设位置以及相对应的信息熵,具体步骤如下:
步骤三一、遗传算法训练过程中采用信息熵作为适应度函数,其理论依据在于:对于服从均值为零,方差为σ2的随机变量s,其概率密度函数为p(s),则其信息熵为:
Figure BDA0002078309700000071
假设预测桥梁关键断面位移的误差矩阵ε=[ε1,···,εmP]∈RH×P(m=1,2,···,P)中每个误差向量εm服从均值为零,方差为
Figure BDA0002078309700000072
的高斯分布,其中P代表桥梁关键断面的个数,此时模型的信息熵为:
Figure BDA0002078309700000073
其中
Figure BDA0002078309700000074
为误差向量εm方差估计值;
步骤三二、遗传算法训练过程中对某一特定数量的倾角仪寻找出最优布设位置的理论依据:通过引入信息熵,将特定数量下倾角仪的最优布设位置问题转化为使预测的P个关键断面位移的误差的信息熵最小的问题,优化目标函数表示为:
Figure BDA0002078309700000075
步骤三三、运用遗传算法对倾角传感器进行优化布设的过程如下:
(1)初始化:随机生成N个传感器布设方案作为N个个体,作为第一代群体。
(2)适应度计算:根据每个个体的传感器布设方案,基于有限元模型中倾角仪布设候选节点的时程倾角数据,通过步骤一中建立的倾角转换位移的结构动力学方程来预测桥梁关键断面的时程位移,然后基于预测的时程位移的误差来计算信息熵,以信息熵的倒数作为每个个体的适应度,即信息熵越小,个体的适应度越大。
(3)选择:基于每个个体的适应度,以概率的方式决定选择能保留基因至下一代的个体,适应度越大的个体,基因遗传到下一代的概率就越大。
(4)交叉:在每一代个体之间随机进行基因交换,即个体之间的传感器位置存在一定程度的概率互换。
(5)变异:每个个体的一个或几个传感器布设位置有一定概率的随机变动。
经过上述一轮操作后进化得到下一代群体。
(6)终止条件判断:当进化代数达到最大进化代数时,以该代群体中的最大适应度个体作为最优解输出。
4.如上所述的步骤四,步骤四具体包括:
步骤四一、以预测的有限元模型中桥梁关键断面时程位移的相对误差等于5%为界限确定出临界信息熵的步骤:
εH×P=0.05×vH×P (9)
其中vH×P为P个桥梁关键断面的实际时程位移矩阵,它的每一列代表一个关键断面的时程位移数据,其中H为每个节点时程位移数据的长度;
Figure BDA0002078309700000081
于是模型的临界信息熵为:
Figure BDA0002078309700000082
其中ΘT代表模型的临界信息熵;
步骤四二、确定出临界信息熵后,在倾角仪最优布置下的各传感器数量所对应的信息熵中找出小于等于临界信息熵的传感器数量的最小值即为最优传感器数量,最优传感器数量对应的最优布设位置即为倾角仪的布设位置。
实施例2
实验室某三跨连续刚构桥,其长度为8米,宽度为0.35米。用Midas Civil 2017建立了三跨连续刚构桥的有限元模型,三跨连续刚构桥的有限元模型由200个梁单元和201个节点组成,桥面上共有127个节点作为安装倾角仪的候选位置,共有9个关键断面作为桥梁动挠度的监测断面,关键断面的位置分别取每一跨的L/4、L/2、3L/4处的断面,三跨连续刚构桥的有限元模型及其关键断面示意图如图1所示。
步骤一、建立倾角仪所在桥梁节点的倾角来预测桥梁关键断面位移的结构动力学方程;
本实例中建立实验室中三跨连续钢构桥布设的M(M=3,4,...,15)个倾角仪传感器所在的节点倾角来预测其9个关键断面动挠度的结构动力学方程,在确定倾角仪的最佳布设位置之前,选择合适的阵型向量数量是非常重要的。采用累积有效竖向模态质量参与率用于确定阵型向量数量。图2显示了三跨连续钢构桥前10阶竖向里兹振型向量的累积有效竖向模态质量参与率。发现前10阶竖向里兹振型向量的累积有效竖向模态质量参与率达到了91%。而当累积有效竖向模态质量参与率大于90%时,可以认为已经获得了足够数量的竖向里兹振型向量来反映结构的模态信息。因此,在此三跨连续钢构桥中使用了前10阶竖向里兹振型向量。竖向位移振型矩阵ΦN与竖向转角振型矩阵ΨN取Midas有限元模型中的前10阶竖向位移里兹振型向量与纵桥向竖直平面内转角里兹振型向量。
步骤二、在有限元模型中创建节点动力荷载来模拟桥面过车。节点动力荷载的大小为40kg,移动速度为0.12m/s,节点动力荷载从Midas模型桥的一侧加载到另一侧。令采样频率为50Hz,从Midas模型桥中提取出127个候选布设节点的竖向时程倾角数据作为有限元模型的训练集,提取出9个关键断面所在节点的竖向时程位移数据作为有限元的验证集。
步骤三、采用遗传算法进行倾角仪布设位置的优化,从127个布设倾角仪传感器的候选位置中,优化出了3个传感器到15个传感器的最佳布设位置,并计算出了3个传感器到15个传感器的最佳布设位置下所对应的信息熵如图3所示;
步骤四、以预测的有限元模型中桥梁关键断面时程位移的相对误差等于5%为界限确定出临界信息熵ΘT=-11.66,从图3中可以看出,当处于最佳布设位置的倾角仪传感器的数量增加到7时,信息熵开始小于临界信息熵且其变化趋于稳定。在有限元模型的数值模拟下,遗传算法优化出的7个最佳布设位置下的倾角仪所预测的桥梁9个关键断面所在节点的竖向时程位移曲线与理论竖向时程位移曲线之间的关系如图4。从图4可以看出,所预测的有限元模型关键断面的竖向时程位移曲线与理论竖向时程位移曲线很好地吻合在一起。因此,从理论上可认为安装7个最佳布设位置的倾角仪传感器足以满足预测动挠度的精度要求。
将遗传算法优化出的7个最佳布设位置的倾角仪安装在实验室的三跨连续钢构桥上(示意图如图5)。载重小车用于模拟经过桥梁的车辆,砝码与模型车的质量均为20kg,载重小车的总质量为40kg,模型载重小车在电机与滑轮装置的共同牵引下在连续钢构桥上以恒定的速度来回移动,载重车的移动速度为0.12m/s。接触式位移计安装在连续钢构桥的9个关键断面下,以测量载重小车在运动过程中桥梁关键断面的实际竖向变形。在实验中,倾角仪和接触式位移计的采样频率均为50Hz。
在实验室中,当载重小车以恒定速度在桥上做往返运动时,倾角仪传感器所测的原始信号如图6,对倾角的原始信号做了小波去噪处理后的信号如图7。实验室中三跨连续刚构桥的9个关键断面测点的预测竖向时程位移与实测竖向时程位移的对比图如图8。
由图8可知,桥梁关键断面的预测竖向时程位移曲线与实测竖向时程位移曲线很好地吻合在一起,验证了本发明所提算法的准确性。

Claims (5)

1.一种基于倾角仪优化布置的桥梁动挠度监测方法,其特征在于,步骤如下:
步骤一、建立倾角仪所在桥梁节点的倾角来预测桥梁关键断面位移的结构动力学方程;
步骤二、在桥梁有限元模型中创建节点动力荷载来模拟桥面过车,提取出有限元模型中倾角仪布设候选节点的倾角时程数据与桥梁关键断面的位移时程数据分别作为训练集和验证集;
步骤三、采用遗传算法进行倾角仪布设位置的优化,优化过程中采用桥梁关键断面预测位移误差的信息熵函数作为遗传算法中的适应度函数,从而得到固定数量下的倾角仪最优布设位置以及相对应的信息熵;
步骤四、以预测的有限元模型中桥梁关键断面时程位移的相对误差等于5%为界限确定出临界信息熵,在倾角仪最优布置下的各传感器数量所对应的信息熵中找出小于等于临界信息熵的传感器数量的最小值即为最优传感器数量,最优传感器数量对应的最优布设位置即为倾角仪的布设位置。
2.根据权利要求1所述的一种基于倾角仪优化布置的桥梁动挠度监测方法,其特征在于,步骤一中结构动力学方程建立步骤如下:
步骤一一、对于梁式结构,根据结构动力学中的振型叠加法建立运动方程:
vN=ΦNqm (1)
其中vN∈RN×1代表竖向位移矩阵,ΦN∈RN×m代表竖向位移振型矩阵,qm∈Rm×1代表竖向位移矩阵的广义坐标向量矩阵,N代表布设倾角仪的个数,m代表所用阵型向量的阶数,在竖直平面内,ΦN是关于桥梁纵桥向长度坐标x的函数矩阵,qm是关于时间t的函数矩阵;
步骤一二、将公式(1)两边同时对x求偏导得:
Figure FDA0002769993910000011
其中θN∈RN×1是纵桥向竖直平面内的倾角矩阵,ΨN∈RN×m代表纵桥向竖直平面内的倾角振型矩阵,采用穆尔-彭罗斯广义逆来求解ΨN的伪逆;
当ΨN为列满秩矩阵时:
qm=(ΨN TΨN)-1ΨN TθN (3)
当ΨN为行满秩矩阵时:
qm=ΨN TNΨN T)-1θN (4)
用来预测同一时刻桥梁其他P个断面的位移:
vP=ΦPqm (5)
其中ΦP∈RP×m,P代表沿桥梁关键断面测点的P个位置。
3.根据权利要求1所述的一种基于倾角仪优化布置的桥梁动挠度监测方法,其特征在于,步骤二包括如下步骤:
步骤二一、在桥梁有限元模型中创建节点动力荷载来模拟桥面过车,节点动力荷载的大小与载重车的重量相一致,在桥面上的移动速度与载重车的过车速度相一致;
步骤二二、将有限元模型运行分析,提取出有限元模型中倾角仪布设候选节点的倾角时程数据与桥梁关键断面的位移时程数据分别作为训练集和验证集。
4.根据权利要求1所述的一种基于倾角仪优化布置的桥梁动挠度监测方法,其特征在于,步骤三包括如下步骤:
步骤三一、遗传算法训练过程中采用信息熵作为适应度函数,其理论依据在于:对于服从均值为零,方差为σ2的随机变量s,其概率密度函数为p(s),则其信息熵为:
Figure FDA0002769993910000021
设预测桥梁关键断面位移的误差矩阵ε=[ε1,···,εmP]∈RH×P(m=1,2,···,P)中每个误差向量εm服从均值为零,方差为
Figure FDA0002769993910000022
的高斯分布,其中P代表桥梁关键断面的个数,此时模型的信息熵为:
Figure FDA0002769993910000023
其中
Figure FDA0002769993910000024
为误差向量εm方差估计值;
步骤三二、遗传算法训练过程中对某一特定数量的倾角仪寻找出最优布设位置的方法:通过引入信息熵,将特定数量下倾角仪的最优布设位置问题转化为使预测的P个关键断面位移的误差的信息熵最小的问题,优化目标函数表示为:
Figure FDA0002769993910000031
步骤三三、运用遗传算法对倾角传感器进行优化布设的过程如下:
步骤三三一、初始化:随机生成N个传感器布设方案作为N个个体,作为第一代群体;
步骤三三二、适应度计算:根据每个个体的传感器布设方案,基于有限元模型中倾角仪布设候选节点的时程倾角数据,通过步骤一中建立的倾角转换位移的结构动力学方程来预测桥梁关键断面的时程位移,然后基于预测的时程位移的误差来计算信息熵,以信息熵的倒数作为每个个体的适应度,即信息熵越小,个体的适应度越大;
步骤三三三、选择:基于每个个体的适应度,以概率的方式决定选择能保留基因至下一代的个体,适应度越大的个体,基因遗传到下一代的概率就越大;
步骤三三四、交叉:在每一代个体之间随机进行基因交换,即个体之间的传感器位置存在一定程度的概率互换;
步骤三三五、变异:每个个体的一个或几个传感器布设位置有一定概率的随机变动;经过上述一轮操作后进化得到下一代群体;
步骤三三六、终止条件判断:当进化代数达到最大进化代数时,以该代群体中的最大适应度个体作为最优解输出。
5.根据权利要求4所述的一种基于倾角仪优化布置的桥梁动挠度监测方法,其特征在于,步骤四具体包括:
步骤四一、以预测的有限元模型中桥梁关键断面时程位移的相对误差等于5%为界限确定出临界信息熵的步骤:
εH×P=0.05×vH×P (9)
其中vH×P为P个桥梁关键断面的实际时程位移矩阵,它的每一列代表一个关键断面的时程位移数据,其中H为每个节点时程位移数据的长度;
Figure FDA0002769993910000032
模型的临界信息熵为:
Figure FDA0002769993910000041
其中ΘT代表模型的临界信息熵;
步骤四二、确定出临界信息熵后,在倾角仪最优布置下的各传感器数量所对应的信息熵中找出小于等于临界信息熵的传感器数量的最小值即为最优传感器数量,最优传感器数量对应的最优布设位置即为倾角仪的布设位置。
CN201910461900.9A 2019-05-30 2019-05-30 基于倾角仪优化布置的桥梁动挠度监测方法 Active CN110319990B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910461900.9A CN110319990B (zh) 2019-05-30 2019-05-30 基于倾角仪优化布置的桥梁动挠度监测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910461900.9A CN110319990B (zh) 2019-05-30 2019-05-30 基于倾角仪优化布置的桥梁动挠度监测方法

Publications (2)

Publication Number Publication Date
CN110319990A CN110319990A (zh) 2019-10-11
CN110319990B true CN110319990B (zh) 2021-01-22

Family

ID=68119258

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910461900.9A Active CN110319990B (zh) 2019-05-30 2019-05-30 基于倾角仪优化布置的桥梁动挠度监测方法

Country Status (1)

Country Link
CN (1) CN110319990B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
IT202100025799A1 (it) * 2021-10-08 2023-04-08 Aisico S R L Metodo per il controllo delle deformazioni delle attrezzature tramite inclinometro

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110991108B (zh) * 2019-11-22 2023-04-25 湖南城市学院 一种机械臂关节扭矩传感器结构设计方法
CN111964638A (zh) * 2020-06-30 2020-11-20 同恩(上海)工程技术有限公司 一种基于拟合反演方式的挠度监测方法及系统
CN112229587A (zh) * 2020-10-10 2021-01-15 哈尔滨工业大学 一种基于倾角仪间接测量高铁桥梁动挠度的方法
CN112945195B (zh) * 2021-01-26 2022-11-04 北京讯腾智慧科技股份有限公司 一种轨道桥梁的列车过桥时的坡度测量方法及装置
CN113283130B (zh) * 2021-04-21 2022-04-15 中国铁路设计集团有限公司 基于应变模态的高速铁路32米标准梁动挠度监测方法

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
KR20060102581A (ko) * 2005-03-24 2006-09-28 한국유지관리 주식회사 교량 및 노반의 상태평가를 위한 하중재하/측정 자동화시스템
KR100850292B1 (ko) * 2007-07-31 2008-08-04 (주)대한이엔씨 구조물 안전진단용 처짐각 측정장치
JP4866434B2 (ja) * 2009-02-24 2012-02-01 ヤマトスチール株式会社 山形鋼の曲がり矯正装置及び山形鋼の製造方法
CN105910778A (zh) * 2016-04-25 2016-08-31 上海同济建设工程质量检测站 一种桥梁挠度的检测方法
CN106338239B (zh) * 2016-11-02 2019-07-23 广州建设工程质量安全检测中心有限公司 一种用于建设工程中横向挠度曲线测量的装置及方法
CN106629404B (zh) * 2016-11-23 2019-03-22 江西乔田重工有限公司 一种起重机械安全监控系统及架桥机主梁的挠度测量方法
CN106644326A (zh) * 2017-02-17 2017-05-10 中国地震局工程力学研究所 基于动挠度的桥梁限载快速监测系统
CN106769157A (zh) * 2017-02-23 2017-05-31 上海喆之信息科技有限公司 一种基于无线网络的桥梁结构可靠性评定系统
CN109060279B (zh) * 2018-06-05 2020-02-18 广东华交科工程科技有限公司 一种倾角传感器测量桥梁挠度的误差分析方法
CN109029882B (zh) * 2018-07-13 2020-05-08 中设设计集团股份有限公司 一种提高基于倾角仪的桥梁挠度测试精度的方法

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
IT202100025799A1 (it) * 2021-10-08 2023-04-08 Aisico S R L Metodo per il controllo delle deformazioni delle attrezzature tramite inclinometro

Also Published As

Publication number Publication date
CN110319990A (zh) 2019-10-11

Similar Documents

Publication Publication Date Title
CN110319990B (zh) 基于倾角仪优化布置的桥梁动挠度监测方法
Zhu et al. Recent developments in inverse problems of vehicle–bridge interaction dynamics
Brown et al. Monitoring of structures using the global positioning system.
Deng et al. Identification of dynamic vehicular axle loads: Demonstration by a field study
CN106289691A (zh) 一种基于微波雷达装置的桥梁分块冲击振动检测方法及检测装置
CN108920766B (zh) 一种基于基函数表示和稀疏正则化的桥梁影响线识别方法
Jiang et al. Identification of vehicles moving on continuous bridges with rough surface
CN105975702B (zh) 斜拉桥健康监测传感器优化布置方法和系统
Li et al. A convolutional neural network‐based full‐field response reconstruction framework with multitype inputs and outputs
Yang et al. Parameters identification of moving load using ANN and dynamic strain
CN110321593B (zh) 基于累积模态质量参与率的桥梁动挠度振型矩阵构造方法
Ntotsios et al. Structural identification of Egnatia Odos bridges based on ambient and earthquake induced vibrations
Li et al. Deep learning models for time-history prediction of vehicle-induced bridge responses: A comparative study
CN117077272A (zh) 一种车桥耦合数值解预测方法
CN113935384B (zh) 一种信号自适应分解和识别的桥梁健康监测方法及系统
Zhao et al. Safety analysis of high-speed trains on bridges under earthquakes using a LSTM-RNN-based surrogate model
CN109406076A (zh) 一种利用位移传感器列阵输出的移动主成分进行梁桥结构损伤定位的方法
CN116822024B (zh) 一种铁路桥上多线列车最不利交会位置的确定方法
Akintunde et al. Singular value decomposition and unsupervised machine learning for virtual strain sensing: Application to an operational railway bridge
CN110344327B (zh) 一种斜拉桥上轨道控制网cpiii点实时高程计算方法
CN116186833B (zh) 索承桥主梁竖向变形相关模式识别与健康诊断机器学习方法
Wan et al. Binary segmentation for structural condition classification using structural health monitoring data
Chen et al. Recent Research and Applications of Numerical Simulation for Dynamic Response of Long‐Span Bridges Subjected to Multiple Loads
Li et al. A machine learning based approach for efficient safety evaluation of the high speed train and short span bridge system
CN114216634A (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
CB03 Change of inventor or designer information

Inventor after: Su Wei

Inventor after: Wang Fei

Inventor after: Meng Fanzeng

Inventor after: Yang Xuelin

Inventor after: Li Zhonglong

Inventor after: Li Hui

Inventor after: Li Shunlong

Inventor after: Wang Xin

Inventor after: Yi Yi

Inventor after: Liu Hongzhan

Inventor after: Wang Xu

Inventor after: Gu Jinshen

Inventor after: Di Hao

Inventor after: Wei Jianfeng

Inventor before: Su Wei

Inventor before: Wang Fei

Inventor before: Meng Fanzeng

Inventor before: Yang Xuelin

Inventor before: Li Zhonglong

Inventor before: Li Hui

Inventor before: Li Shunlong

Inventor before: Wang Xin

Inventor before: Yi Yi

Inventor before: Liu Hongzhan

Inventor before: Wang Xu

Inventor before: Gu Jinshen

Inventor before: Di Hao

Inventor before: Wei Jianfeng

CB03 Change of inventor or designer information
GR01 Patent grant
GR01 Patent grant