[发明专利]基于QMR算法的桥面多轴移动荷载的识别方法

说明书

技术领域

本发明属于桥面移动荷载识别技术领域，尤其涉及一种由桥梁位移识别桥面多轴移动荷载的方法。

背景技术

我国桥梁现状是“重建轻养”，从1999年到2013年，国内媒体公开报道我国因各种原因垮塌的桥梁多达110余座，其中尚不包括汶川地震引起的桥梁垮塌。引起桥梁损伤与破坏原因可归纳为外部因素和内部因素，其中外部因素中由于汽车超载导致桥梁疲劳损伤和耐久性降低占据主导地位，内部因素则主要是桥梁自身承载力降低和材料强度退化。

随着我国公路交通的爆发式增长，许多桥梁实际承受的车流量较早期设计值增加很多，车速和车重的增加均会对桥梁产生不利影响，而大型多轴车辆尤其是超载多轴车辆的出现明显加剧了桥梁破坏的风险。

我国公路超限站在控制车辆超重方法做出许多工作，但目前测量方法多是采用地磅技术，即通过停车称重来实现车辆总重的测量。在发展快速交通的趋势下，如何在车辆行驶过程中精确车辆荷载具有重要的工程实际意义，尤其是对多轴货车各轴荷载的精确测量对保护桥梁的安全性和耐久性都有很大帮助。

现有的移动荷载识别技术多针对常规两轴车辆进行识别，不能对多轴车辆荷载进行识别，因此急需一种能够对桥面多轴移动车辆荷载进行识别的方法。

发明内容

本发明的目的是提供一种仅需测量桥梁位移响应即可快速高效的识别桥面多轴移动车辆荷载，识别精度高且不影响桥面车辆正常通行。

为达到上述目的，本发明采用的技术方案是：一种基于QMR算法的桥面多轴

移动荷载的识别方法，包括以下步骤：

1)、在桥梁底面对应位置x₁,x₂,…x_m处分别粘贴m个位移传感器，测得桥面多轴移动车辆荷载fk(t)在x位置处t时刻的位移为v(x,t)，k＝1,2,3…，为车辆轴数；

2)、建立车桥系统振动微分方程：取桥梁长度为L，抗弯刚度为EI，桥梁单位长度质量为ρ，考虑粘性阻尼并取阻尼系数为C，忽略桥梁的剪切变形和转动惯量，桥面多轴移动车辆荷载f_k(t)以速度c自梁左端支承处向右移动，则车桥系统的振动微分方程为：

其中δ(x-ct)是狄拉克函数；

方程(1)的边界条件为：

v(0,t)＝0，v(L,t)＝0，v(x,0)＝0，

3)、对方程(1)求解；

4)、建立桥梁在k轴车辆荷载作用下，由位移响应识别多轴移动荷载系统方程：v_(m×1)＝S_(m×k)·f_(k×1)(2)

v_(m×1)为移动荷载f_k(t)在x₁,x₂,…x_m处的实际位移，且m≥k；S_(m×k)为已知的系统矩阵；f_(k×1)为所求的k轴移动荷载；

式(2)的离散形式表示为：

其中

5)、采用QMR算法求得多轴移动荷载的精确值；

对方程(2)的系统矩阵S和位移响应v采用基于Lanczos双A-正交的QRM算法求解：

首先给定f₀，f₀可以由最小二乘法得到，计算r₀＝v-Sf₀，依次可求得r₁＝v-Sf₁，r₂＝v-Sf₂等，r₀、r₁、r₂分别是残差第一步，第二步，第三步迭代产生的残差。令β＝||r₀||₂，令ε₁＝r₀/β，选取使得标准内积<ω₁，Sε₁>＝1，通过如下迭代产生向量ε_j和ω_j，以及标量δ_j，β_j，j＝1,2,…,b

δ_j+1ε_j+1＝Sε_j-β_jε_j-1-α_jε_j