基于随机车流下宏应变能的桥梁损伤快速识别方法与流程

[0001]
本发明涉及桥梁损伤识别技术领域，具体是一种基于随机车流下宏应变能的桥梁损 伤快速识别方法。


背景技术：

[0002]
因为长期往复的上部车辆荷载作用会引起桥梁性能的逐步下降，为了避免潜在事故 的发生，保障上部通行人员和车辆安全，需要对桥梁进行实时检测。但目前常用的方法 是通过荷载试验来完成，但实践下来发现荷载试验开展效率低，作业时间长，需要中断 交通，会严重影响正常车辆通行。因此基于健康监测的桥梁快速损伤识别算法尤为必要。 目前，许多学者已经开发出很多损伤识别算法，但这其中很大一部分的使用工况受限， 对桥梁上部车流状况有严格的限制，只能应用于类似于荷载试验中跑车试验的简单工 况，而无法直接在随机车流工况下工作，因此应用场景有限。


技术实现要素：

[0003]
发明目的：为解决上述技术问题，本发明提供一种基于随机车流下宏应变能的桥梁 损伤快速识别方法。
[0004]
技术方案：一种基于随机车流下宏应变能的桥梁损伤快速识别方法，包括以下步骤：
[0005]
步骤1：在桥梁下方沿桥跨方向连续安装n个宏应变传感器:s1,s2,
…
sn；
[0006]
步骤2：采集一段随机车流经过时桥梁的宏应变响应mr1,mr2,
…
mrn；
[0007]
步骤3：对各个传感器采集到的宏应变响应在时域内进行积分，得到宏应变能向量 mev；
[0008]
步骤4：利用宏应变能向量的无穷范数，对宏应变能向量进行正则化，得到正则化 的宏应变能向量rmev；
[0009]
步骤5：根据识别得到的正则化的宏应变能向量rmev和标准值对比，得到桥梁损 伤位置和程度。
[0010]
优选地，步骤1中的宏应变指结构表面一段区域内各点应变的积分，宏应变与点应 变的关系如下：
[0011][0012]
其中，为宏应变，ε(x)为点应变。
[0013]
优选地，步骤2中，一段随机车流从正常车流中截取，满足车流的完整性，在车流 的起始和结束后桥上保持空车状态。
[0014]
优选地，步骤3中，宏应变能向量由各个传感器测得的宏应变能me组成，宏应变 能me计算公式为：
[0015]
[0016]
其中为第i个传感器采集到的宏应变时程，t1，t2为单车上下桥时间。
[0017]
优选地，步骤4中，正则化的宏应变能rmev向量计算公式为：
[0018][0019]
优选地，步骤5中，正则化的宏应变能向量rmev的标准值为1，当rmev大于1 时则说明相应位置存在损伤，损伤程度β计算公式为：
[0020]
损伤程度：
[0021]
有益效果：本发明可以在随机车流工况下有效准确快速判断结构损伤位置和程度， 无需中断交通或限制交通状况。
附图说明
[0022]
图1为本发明的基于随机车流下宏应变能的桥梁损伤快速识别方法的步骤流程示意 图；
[0023]
图2为当车辆通过完好桥梁方法应用示意图；
[0024]
图3为当车辆通过带损伤桥梁方法应用示意图；
[0025]
图4为实际车辆产生的静态宏应变响应组成示意图；
[0026]
图5为桥梁空间受力示意图；
[0027]
图6为桥梁空间荷载横向分布示意图；
[0028]
图7为桥梁空间受力解耦示意图。
具体实施方式
[0029]
下面结合附图对本发明作进一步说明。
[0030]
如图1-7所示，本发明的一种基于随机车流下宏应变能的桥梁损伤快速识别方法， 包括以下步骤：
[0031]
步骤1：在桥梁下方沿桥跨方向连续安装n个宏应变传感器:s1,s2,
…
sn，如图2 所示；其中宏应变传感器可以选用大标距的电阻式应变片、光纤光栅宏应变传感器或其 它可用于测量宏应变的传感设备。n个宏应变传感器在布置时最佳方案为首尾相连布满 桥梁整个跨径，如图2和图3所示，当成本有限时也可以仅在结构关键区域布置，如跨 中、四分之一跨等。
[0032]
步骤2：采集一段随机车流经过时桥梁的宏应变响应macro-strain response(mr): mr1,mr2,
…
mrn，如图2和图3所示；其中一段随机车流需要从正常车流中截取，需 要满足车流的完整性，在车流的起始和结束后桥上保持空车状态。
[0033]
步骤3：对各个传感器采集到的宏应变响应在时域内进行积分，得到宏应变能向量 macro-strain energy vector(mev)；宏应变能向量由各个传感器测得的宏应变能 (macro-strain energy,me)组成，如图2、图3，宏应变能计算公式为：
[0034]
[0035]
其中为第i个传感器采集到的宏应变时程，t1，t2为单车上下桥时间。
[0036]
步骤4：利用mev的无穷范数，对mev进行正则化，得到正则化的宏应变能向量 regularized mev(rmev)，如图2，图3；
[0037]
正则化的宏应变能向量rmev计算公式为：
[0038][0039]
步骤5：根据识别得到的rmev和标准值对比，得到桥梁损伤位置和程度。其中 rmev的标准值为1，当rmev大于1时则说明相应位置存在损伤，损伤程度β计算公 式为：
[0040]
损伤程度：
[0041]
本发明中，首先根据结构力学中应变影响线的定义推导宏应变影响线方程。根据结 构力学，在简支梁任意截面xi处的应变影响线方程为：
[0042][0043]
其中(ei)
i
为截面xi处的截面刚度，h代表中和轴高度。根据宏应变的定义，第m 个传感器标距范围内的宏应变与标距范围内结构各点的应变值间的转化关系如下式：
[0044][0045]
其中lg为宏应变传感器的标距长度，m
·
lg和(m+1)
·
lg分别为第m个宏应变传感器 左右两端的坐标。因此可得到第m个标距范围内的宏应变影响线表达式：
[0046][0047]
其中代表传感器标距内结构等效抗弯刚度。
[0048]
当将车辆单个车轴简化为一集中力时，车轴移动过桥这一过程正好可以用影响线来 描述。此时需要将移动力的位置x用时间t来表示。如何表示完全取决于车辆的行进速 度，当车辆匀速行驶时，移动力位置函数则可以表示为x(t)＝v
·
t。对于大多数公路桥梁， 因为它们以中小跨桥梁占主导，跨度有限，车辆通过桥梁时用时很短，在这期间速度基 本上难以发生较大变化，因此当桥梁上部有车辆通过时，基本可将其视为匀速通过。故 当一匀速移动车轴作用在桥梁上，在第m个传感器标段内产生的静态宏应变时程可表示 为：
表示为：
[0058][0059]
其中表示第i辆车第k轴轴重分配的荷载，如图6所示，和分别是被监测 梁体与第i辆车左、右侧车轮间的相对距离。
[0060]
对于具体的随机交通车流，其响应即可通过这两种基本工况经随机组合得到。通过 对随机车流作用下宏应变响应沿时间进行积分，推导出静态宏应变响应时程与时间轴间 覆盖的总面积，即宏应变能me为：
[0061][0062]
其中t
s
和t
e
分别为整条随机车流开始和结束的时间，和分别是第i辆车第k 个车轴开始过桥和完全通过桥梁的时间。单位静态应变响应所包围的时程面积 等于：
[0063][0064]
其中f(m,lg,h,l)是一与几何参数m，lg，h和l相关的函数。
[0065]
可以得到随机车流产生的宏应变能me等于：
[0066][0067]
该式表明由第m个传感器测得的宏应变能与传感器覆盖范围内结构的平均刚度成 反比。基于这一关系可以发现，该值可用来作为损伤指标来识别结构损伤状况。对一段 随机车流而言，各个车辆的横向位置、车速、轴重都是恒定的，而对一座桥梁而言，当 传感器安装好后，函数fm＝f(m,lg,h,l)的值就已固定下来。但此时由n个宏应变传感 器阵列测得的宏应变能向量mev在不同的随机车流下会因荷载不同而变化，难以进行 比较。为便于比较，实现在不同随机车流状况下评估桥梁状况，这里进一步通过正则化 的宏应变能向量rmev用于桥梁损伤识别。
[0068]
经过正则化，rmev会始终大于等于1，在不同随机车流状况下也可以相互比较， 不受车流状况影响，只与桥梁纵向刚度分布情况相关，亦即只与结构损伤状况相关。由 于结构刚度和指标rmev之间的反比例关系，当在某一传感器标距内的结构出现损伤， 结构刚度降低，相对应的rmev值将会升高表征出结构损伤位置。
[0069]
对于结构损伤程度，假设第m个标距内刚度退化程度为β，根据指标rmev的定 义，结构损伤程度可以由下式计算得：
[0070][0071]
其中代表了第m个传感器标距内的结构刚度，为结构完好时的 刚度。
[0072]
以上所述仅是本发明的优选实施方式，应当指出：对于本领域的普通技术人员来说， 在不脱离本发明原理的前提下，还可以做出若干改进和润饰，这些改进和润饰也应视为 本发明的保护范围。