文章快速检索    
  同济大学学报(自然科学版)  2019, Vol. 47 Issue (8): 1115-1122, 1130.  DOI: 10.11908/j.issn.0253-374x.2019.08.005
0

引用本文  

樊学平, 屈广, 刘月飞. 桥梁时变可靠度指标的改进粒子滤波预测算法[J]. 同济大学学报(自然科学版), 2019, 47(8): 1115-1122, 1130. DOI: 10.11908/j.issn.0253-374x.2019.08.005.
FAN Xueping, QU Guang, LIU Yuefei. Improved Particle Filter Prediction Algorithm of Time-Variant Reliability Indices for Bridges[J]. Journal of Tongji University (Natural Science), 2019, 47(8): 1115-1122, 1130. DOI: 10.11908/j.issn.0253-374x.2019.08.005

基金项目

国家自然科学基金(51608243)

第一作者

樊学平(1983—), 男, 副教授, 硕士生导师, 工学博士, 主要研究方向为桥梁结构安全预后与损伤预后.E-mail:fxp_2004@163.com

文章历史

收稿日期:2018-04-29
桥梁时变可靠度指标的改进粒子滤波预测算法
樊学平 1,2, 屈广 1, 刘月飞 1,2     
1. 兰州大学 土木工程与力学学院, 甘肃 兰州 7300001;
2. 兰州大学西部灾害与环境力学教育部重点实验室, 甘肃 兰州 730000
摘要:基于健康监测时间序列数据, 提出了桥梁动态可靠度指标的改进粒子滤波预测方法.首先, 利用监测极值数据建立动态模型, 将其作为粒子滤波算法的状态方程和监测方程; 然后, 采用贝叶斯动态线性模型(BDLM)为粒子滤波器提供随时间更新的动态建议分布, 以解决传统粒子滤波算法的样本退化问题, 同时增加了粒子滤波算法的鲁棒性及自适应性; 进而利用改进的粒子滤波算法(IPF), 结合极值监测数据实现结构极值的动态预测, 并结合一次二阶矩(FOSM)可靠性方法, 实现桥梁结构可靠度指标的动态预测; 最后通过在役桥梁工程实例与设计试验对所提模型和方法的合理性与有效性进行验证.
关键词桥梁健康监测数据    动态模型    建议分布函数    改进粒子滤波算法    一次二阶矩方法    动态可靠度指标预测    
Improved Particle Filter Prediction Algorithm of Time-Variant Reliability Indices for Bridges
FAN Xueping 1,2, QU Guang 1, LIU Yuefei 1,2     
1. School of Civil Engineering and Mechanics, Lanzhou University, Lanzhou 730000, China;
2. Key Laboratory of Mechanics on Disaster and Environment in Western China of the Ministry of Education, Lanzhou University, Lanzhou 730000, China
Abstract: This paper proposes an improved particle filter (IPF) prediction approach of dynamic reliability indices for bridges based on monitoring time series data. First, the dynamic models, which can provide state equation and monitoring equation for the IPF, are built with the monitoring extreme data of bridges. Next, the Bayesian dynamic linear model (BDLM) is utilized to produce the real-time updated proposal distribution for IPF in order to solve the sample degradation problem and increase the robustness and adaptability of the traditional particle filter. After that, by using the IPF approach, the structural extreme information is dynamically predicted based on the monitoring extreme data, and dynamic reliability indices of bridges are predicted by using the first order second moment (FOSM) reliability method. Finally, three existing bridges and a designed experiment are provided to illustrate the feasibility and application of the proposed model and method.
Key words: bridge health monitoring data    dynamic models    proposal distribution function    improved particle filter(IPF) prediction method    first order second moment(FOSM) method    dynamic reliability indices prediction    

目前,结构健康监测(structural health monitoring, SHM)领域的研究主要集中在利用传感器采集数据和监测数据的应用两个方面.前者主要集中在数据压缩、数据恢复、数据获得技术和系统组装技术等方面[1-5],目前在硬件和技术上已经趋于成熟,对于后者而言,国内外学者的研究主要集中在结构损伤识别、结构模态参数识别以及结构模型修正等方面[6-8].基于SHM数据的桥梁可靠性预测及评定,国内外学者多是基于离线监测信息展开研究[9-13],而基于实时监测信息的结构性能动态预估相对较少,且离线的结构性能预测及评定模型在实际工程应用中存在一定局限性,如:基于一次回归函数与常值函数的桥梁可靠性预测分析[9-10],均未考虑监测变量的随机性和非平稳性;基于ARMA(autoregressive moving average model)模型的桥梁构件荷载效应预测[11],以及基于单一或组合贝叶斯动态模型的桥梁可靠性动态预测[12-13],多是基于平稳监测信息进行可靠性分析,而考虑到桥梁动力响应的复杂性,基于非线性及非平稳监测信息的结构性能预测方法还需进一步研究.近年来,粒子滤波器因在非线性问题中出众的性能而倍受关注[14].在SHM领域中,粒子滤波器在结构应力预测、结构可靠性动态预测、结构系统识别以及结构损伤识别等方面已取得一些研究成果,如:基于高斯混合粒子滤波器、改进高斯混合粒子滤波器以及折扣高斯粒子滤波器的桥梁可靠性预测[12, 15-17].但粒子滤波器假定重要性采样是指能够从一个合理的后验建议密度函数中得到一组覆盖真实状态的样本集合.因此,找到一个最优建议密度函数来指导采样过程,就能大大提高粒子的利用率,进而可以有效降低粒子退化的影响,提高系统的识别精度.桥梁长期运营中,由于其构件受到多种荷载耦合作用,SHM系统所获取的结构动力响应数据十分复杂,并且受诸如气候条件、桥下水位等外在环境因素影响较大,因此桥梁应力(或挠度)状况时常会出现变化,这就要求所选取的粒子滤波建议分布具有长期稳定性、受外界干扰小、精度高等特点,并且具有能随桥梁状况变化而不断更新的良好自适应性.因而,如何选取具有以上特性的建议分布就成为利用粒子滤波方法预测桥梁时变可靠性的难点所在.

为合理应用监测信息对桥梁动态可靠度指标进行实时在线预测分析,基于桥梁监测极值应力(或挠度)信息,利用贝叶斯动态线性模型和粒子滤波器相融合的改进粒子滤波预测算法,结合FOSM(first order second moment)方法对桥梁可靠度指标进行了动态预测,并进行了实例和设计试验验证分析.改进粒子滤波方法利用贝叶斯动态线性预测模型为动态粒子提供建议分布,并充分利用动态健康监测数据对建议分布进行递推更新,有效提高了粒子的利用率,并且引入贝叶斯动态模型折扣因子[12],使得粒子滤波在递推过程中具有较好的自适应性.

1 监测数据的动态模型

SHM系统在长期运营过程中积累了大量监测应力(或挠度)信息,定义每小时(或每天)监测应力(或挠度)的极大值为监测极值应力(或挠度)信息.定义{Xt, t=1, 2, …}为极值应力(或挠度)状态时间序列,Xt表示t时刻的极值应力(或挠度)状态.所建动态方程基于以下两点假设:

(1) {Xt, t=1, 2, …}具有马尔科夫性.

(2) 监测信息{Zt, t=1, 2, …}之间相互独立,且Zt仅和状态变量Xt相关,ZtXt呈线性关系.

桥梁极值应力(或挠度)的动态模型如下:

状态方程为

$ X_{t}=f\left(X_{t-1}\right)+w_{t}, w_{t} \sim N\left(0, W_{t}\right) $ (1)

监测方程为

$ Z_{t}=X_{t}+v_{t}, v_{t} \sim N\left(0, V_{t}\right) $ (2)

初始状态信息为

$ p\left(X_{t-1} | D_{t-1}\right)=N\left(m_{t-1}, C_{t-1}\right) $ (3)

式(1)~(3)中:Xt为状态变量;f(·)为状态转移函数,可以通过极值应力(或挠度)监测数据利用五点三次平滑处理之后的时间序列近似回归拟合得到[12]wt为状态白噪声; Wt为状态误差的方差;Ztt时刻的极值应力(挠度)监测值;vt为监测白噪声;Vt为监测误差的方差;初始状态变量的概率密度函数(probability density function, PDF) p(Xt-1Dt-1)可根据既有极值应力(或挠度)信息平滑得到的时间序列样本统计得到;N(·)为正态PDF.

基于先验信息集Dt和当前监测数据Zt+1估计状态Xt+1的后验PDF, 即p(Xt+1Dt+1),求解过程包含以下两个阶段:

(1) 预测.已知t时刻状态的后验PDF p(XtDt),利用式(1)可得t+1时刻状态Xt+1的先验PDF为

$ p\left(X_{t+1} | D_{t}\right)=\int p\left(X_{t+1} | X_{t}\right) p\left(X_{t} | D_{t}\right) \mathrm{d} X_{t+1} $ (4)

(2) 测量更新.在t+1时刻,利用当前监测极值Zt+1,由Bayes公式更新式(4),可得后验PDF为

$ p\left(X_{t+1} | D_{t+1}\right)=\frac{p\left(X_{t+1} | D_{t}\right) p\left(Z_{t+1} | X_{t+1}\right)}{p\left(Z_{t+1} | D_{t}\right)} $ (5)

其中

$ \begin{aligned} p\left(Z_{t+1} | D_{t}\right)=& \int p\left(X_{t+1} | D_{t}\right)\cdot \\ & p\left(Z_{t_{t+1}} | X_{t+1}\right) \mathrm{d} X_{t+1} \end{aligned} $ (6)

为极值的一步预测PDF.

2 改进的粒子滤波算法

结合重要性采样与重采样方法,改进的粒子滤波器(IPF)通过蒙特卡罗方法以参数化的分布逼近式(4)和式(5)中的预测PDF与后验状态PDF,是式(4)和(5)近似求解的一种策略.

2.1 基于BDLM(Bayesian dynamic linear model)的近似重要性采样分布

粒子滤波的有效性依赖于重要性采样分布的合理性.较为合理的重要性采样分布能够有效减少具有偏差的样本,提升粒子滤波性能.因此,引入BDLM[12]在重要性采样阶段近似产生重要性采样PDF.假设采样过程中抽取了n个粒子,则更新粒子集合的步骤如下:

(1) t时刻每个粒子的状态后验PDF为

$ p\left(X_{t}^{(i)} | D_{t}\right)=N\left(m_{t}^{(i)}, C_{t}^{(i)}\right), i=1, 2, \cdots, n $ (7)

式中:考虑到粒子的随机性,mt(i)Ct(i)分别为t时刻第i个粒子服从高斯分布的均值和方差.

(2) t+1时刻的状态先验PDF为

$ p\left(X_{t-1}^{(i)} | D_{t}\right)=N\left(a_{t}^{(i)}, R_{t}^{(i)}\right) $ (8)

式中: $a_{t + 1}^{(i)} = {{\boldsymbol{G}}_{t + 1}}m_t^{(i)}{\boldsymbol{G}}_{t + 1}^{\rm{T}};R_{t + 1}^{(i)} = {\boldsymbol{G}}_{t + 1}^{(i)}{\boldsymbol{C}}_t^i{\boldsymbol{G}}_{t + 1}^{(i){\rm{T}}} + W_{t + 1}^{(i)};{\boldsymbol{G}}_{t + 1}^{(i)} $为已知状态转移矩阵;状态误差方差Wt+1(i)Wt+1(i)=Ct(i)(1/δ-1), δ为折扣因子,在非线性系统中,取δ=0.8能够大大增加算法的自适应性.

(3) t+1时刻的一步预测PDF为

$ p\left(Z_{t+1} | D_{t}\right)=N\left(f_{t+1}^{(i)}, Q_{t+1}^{(i)}\right) $ (9)

式中: $ f_{t+1}^{(i)}=\boldsymbol{G}_{t+1} a_{t+1}^{(i)} ; Q_{t+1}^{(i)}=R_{t+1}^{(i)}+V, V$为已知监测误差方差.

(4) t+1时刻的状态后验PDF为

$ p\left(X_{t+1}^{(i)} | D_{t+1}\right) \approx N\left(m_{t+1}^{(i)}, C_{t+1}^{(i)}\right) $ (10)

式中:$m_{t+1}^{(i)}=a_{t+1}^{(i)}+\boldsymbol{A}_{t+1}^{(i)} e_{t+1}^{(i)} ; C_{t+1}^{(i)}=R_{t+1}^{(i)}-\boldsymbol{A}_{t+1}^{(i)} \boldsymbol{A}_{t+1}^{(i)} $·$Q_{t+1}^{(i)} ; \boldsymbol{A}_{t+1}^{(i)}=R_{t+1}^{(i)} / Q_{t+1}^{(i)} ; e_{t+1}^{(i)}=Z_{t+1}-f_{t+1}^{(i)} $为一步预测误差.

式(10)用来近似模拟重要性采样分布.

2.2 改进的粒子滤波算法

粒子滤波器能够为监测值与状态值的后验概率分布进行高精度表达[18],在非线性系统中具有出色的表现.然而,传统粒子滤波器以给定的高斯分布作为建议分布,并没有充分利用最新监测信息对建议分布进行动态更新[19],因此通常需要大量粒子样本来近似接近系统的后验PDF,并且在状态更新时计算量较大,具有一定盲目性.特别是当桥梁的应力(或挠度)环境发生变化时,由粒子滤波产生的样本就会出现较大偏差.因而有效减少样本数量的自适应采样策略是该算法的重点所在.因此,将自适应性及鲁棒性良好的BDLM应用到粒子滤波中,为其提供重要性采样,不仅可以获取更为准确的建议分布,而且能够提升粒子滤波算法的鲁棒性与自适应性,同时大大减少了粒子滤波算法的样本需求量,降低了算法复杂度.

改进后的算法步骤如下:

(1) 初始化,t=0,从先验PDF p(X0)中抽取n个初始化状态粒子,得到初始采样集合{Xt(i)}i=1n.

(2) 由式(7)~(10)递推得到每个粒子的建议分布,亦即在t时刻用BDLM以及最新监测信息Zt+1更新粒子,则第i个粒子更新为

$ \hat{X}_{t+1}^{(i)} \sim q\left(X_{t+1}^{(i)} | D_{t+1}\right) \approx N\left(m_{t+1}^{(i)}, C_{t+1}^{(i)}\right) $ (11)

利用式(2)、(8)和(11)可得粒子权重为

$ \omega_{t+1}^{(i)}=\frac{p\left(Z_{t+1} | \hat{X}_{t+1}^{(i)}\right) N\left(\hat{X}_{t-1}^{(i)} | a_{t-1}^{(i)}, R_{t-1}^{(i)}\right)}{q\left(\hat{X}_{t+1}^{(i)} | D_{t+1}\right)} $ (12)

式中:$p\left(Z_{t+1} | \hat{X}_{t+1}^{(i)}\right)=N\left(Z_{t+1}, \hat{X}_{t+1}^{(i)}, V\right) $.

将权重进行归一化,可得

$ \hat \omega _{t + 1}^i = \frac{{\omega _{t + 1}^i}}{{\sum\limits_{i = 1}^n {\omega _{t + 1}^i} }} $ (13)

(3) 利用重采样算法[20],根据归一化权值的大小对粒子集合进行复制和淘汰.

(4) 粒子集合均值Xt+1t+1时刻后验状态值.

(5) 利用式(1)~(2)及(7)~(9)式可得t+2时刻的一步预测PDF $ N\left(f_{t+2}^{(i)}, Q_{t+2}^{(i)}\right)$,则t+2时刻的预测值为

$ {\hat Z_{t + 2}} = \frac{1}{n}\sum\limits_{i = 1}^n {f_{t + 2}^{(i)}} $ (14)

(6) 重复(2)~(5)步,可实现监测变量的动态预测.

一步预测精度可由均方误差(EMSE)来衡量,即

$ {E_{{\rm{MSE}}}} = \frac{1}{T}\sum\limits_{t = 1}^T {{{\left( {{Z_t} - {{\hat Z}_t}} \right)}^2}} $ (15)

显然,均方误差越小,预测精度越高.

改进算法流程如图 1所示.

图 1 改进的粒子算法流程 Fig.1 Flowchart of IPF
3 可靠度指标预测计算 3.1 工程实例1预测公式

基于天津富民桥主缆的应力监测数据,结合FOSM方法对主缆可靠性进行动态预测.主缆钢丝的抗拉强度按照[σ]=1 670 MPa进行计算,主缆截面失效模式对应的功能函数为

$ g\left(\sigma_{\mathrm{p}, t}\right)=[\sigma]-\sigma_{\mathrm{p}, t} $ (16)

式中:σp, tt时刻预测的极值应力,且随机变量之间相互独立.

主缆钢丝抗拉强度的平均值为μσ=1 670 MPa,变异系数δσ=0.15,则利用FOSM方法可得主缆的动态监测可靠度指标βt

$ {\beta _t} = \frac{{{\mu _\sigma } - {\mu _{{\sigma _t}}}}}{{\sqrt {{V_\sigma }} }} $ (17)

式中:μσVσ分别为主缆钢丝强度的均值与方差,Vσ=μσ2δσ2μσtt时刻的极值应力监测值.

考虑到预测过程的不确定性与随机性,则动态可靠度指标预测公式为

$ {\beta _t} = \frac{{{\mu _\sigma } - {\mu _{{\sigma _{{\rm{p}}, t}}}}}}{{\sqrt {{V_\sigma } + {V_{{\sigma _{{\rm{p}}, t}}}}} }} $ (18)

式中:μσp, tVσp, t分别为t时刻主缆极值应力基于改进粒子滤波器的一步预测均值和方差.

本算例中的极值应力为每小时监测应力的极大值.

3.2 工程实例2预测公式

参考已有文献[12],可得美国I-39北桥第二跨横梁可靠度指标预测公式为

$ {\beta _{\rm{p}}} = \frac{{{\mu _{\rm{R}}} - {\mu _{\rm{S}}} - {\mu _{\rm{C}}} - {\gamma _{\rm{M}}}{\mu _{\rm{M}}}}}{{\sqrt {\sigma _{\rm{R}}^2 + \sigma _{\rm{S}}^2 + \sigma _{\rm{C}}^2 + {{\left( {{\gamma _{\rm{M}}}{\sigma _{\rm{M}}}} \right)}^2}} }} $ (19)

式中:μMσM分别为基于改进粒子滤波器预测得到的极值应力平均值和标准差;μRσR分别为按照规范计算的抗力的平均值和标准差;μSσS分别为由钢板恒载所引起的应力的平均值和标准差;μCσC分别为由混凝土恒载所引起的应力的平均值和标准差;γM=1.15是传感器的修正系数.

由文献[12]可知:μR=380 MPa, σR=26.6 MPa,结合由结构自重产生的应力分布参数, μS=116.3 MPa, σS=116.3 MPa×0.04=4.65 MPa, μC=108.8 MPa以及σC=108.8 MPa×0.04=4.35 MPa,可得第二跨横梁的可靠度指标预测公式为

$ \beta_{\mathrm{p}}=\frac{155-1.15 \mu_{\mathrm{M}}}{\sqrt{27.351^{2}+\left(1.15 \sigma_{\mathrm{M}}\right)^{2}}} $ (20)

参考式(17),可得第二跨横梁的可靠度指标监测值βt

$ \beta_{t}=\frac{155-1.15 M_{t}}{\sqrt{27.351^{2}}} $ (21)

式中:Mtt时刻的极值应力监测值.

本算例中的极值应力为每天监测应力的极大值.

3.3 工程实例3预测公式

参考已有文献[12],可得伊通河桥上游主梁动态可靠度指标预测公式为

$ \beta_{\mathrm{p}}=\frac{\mu_{\omega}-\mu_{f(t)}}{\sqrt{\sigma_{\omega}^{2}+\sigma_{f(t)}^{2}}}=\frac{162.5-\mu_{f(t)}}{\sqrt{18.6875^{2}+\sigma_{f(t)}^{2}}} $ (22)

式中:μf(t)σf(t)2分别为基于改进粒子滤波器的极值挠度一步预测平均值与方差;μωσω2分别为按照规范计算的挠度极限平均值和方差.

参考式(17)和(21),可得主梁动态可靠度指标监测值βt

$ \beta_{t}=\frac{162.5-f(t)}{\sqrt{18.6875^{2}}} $ (23)

式中:f(t)为t时刻的极值挠度监测值.

本算例中的极值挠度为每小时监测挠度的极大值.

4 算例分析 4.1 工程实例1计算分析

天津富民桥的结构及其示意图分别如图 2图 3所示.富民桥总长340.3 m,桥宽40 m.桥梁动力响应复杂,结构受温度荷载影响较大.基于此悬索桥主缆的监测信息,对主缆可靠性进行动态预测.

图 2 天津富民桥 Fig.2 Fumin Bridge in Tianjin
图 3 富民桥结构示意图 Fig.3 Structure of Fumin Bridge

此主缆共安装了2个传感器,取其中历史监测极值应力绝对值较大的FBG01192传感器作为实时监控对象.定义每小时监测应力的最大值为监测极值应力.从2009年8月1日到2009年9月3日,对此桥进行了810 h的健康监测.其中前300 h的监测极值应力如图 4所示.对其采用五点三次平滑方法[9]处理后的数据如图 5所示,处理后的数据近似作为初始状态信息,此初始状态信息用正态PDF模拟.

图 4 极值应力监测值(工程实例1) Fig.4 Monitoring extreme stress data
图 5 极值应力初始值与监测值(工程实例1) Fig.5 Initial data and monitoring data

采用前300 h的极值应力信息,并结合式(1)~(3),建立的动态模型如下:

状态方程为

$ X_{t}=X_{t-1}-0.0066+w_{t}, w_{t} \sim N\left(0, W_{t}\right) $ (24)

监测方程为

$ Z_{t}=X_{t}+v_{t}, v_{t} \sim N(0, 0.1032) $ (25)

初始状态信息为

$ X_{t-1} | D_{t-1} \sim N(147.6341, 6.2617) $ (26)

利用改进的粒子滤波算法,结合式(7)~(15)与式(24)~(26),基于第300小时到第809小时的监测极值应力数据,对第301小时到第810小时的极值应力进行实时动态一步预测,结果如图 6所示.一步预测均方误差可由式(15)计算,连续50次滤波结果的EMSE图 7.可以看出, 一步预测值与实测值几乎重合,并且均方误差较小(0.068左右波动)且平稳可观,验证了模型对于监测应力预测的有效性.将预测结果与基本粒子滤波(PF)[21]及差分自回归移动平均模型(ARIMA)[22]进行比较,预测结果对比如图 8所示.经计算得到粒子滤波的均方误差为0.146 8(平均50次滤波EMSE的均值),ARIMA模型的均方误差为0.475 8.由此可见改进粒子滤波器的预测精度十分可观,进一步验证了所提模型的准确性.

图 6 基于改进粒子滤波器的极值应力预测值与实测值(工程实例1) Fig.6 Monitoring and predicted extreme stresses based on IPF
图 7 连续50次滤波的预测均方误差(工程实例1) Fig.7 MSEs of 50 consecutive predictions
图 8 3种算法极值应力预测结果比较(工程实例1) Fig.8 Comparison of three algorithms

结合式(17)和式(18)对主缆可靠性进行实时预测分析,得到结果如图 9所示,可以看出可靠度指标预测值与实测值几乎重合,很好地预测了可靠度指标的变化趋势,而考虑了不确定性与随机性的可靠度指标预测值略低于实测值,但由于考虑了不确定性,预测结果更符合工程实际.

图 9 可靠度指标的实测值与一步预测值(工程实例1) Fig.9 Monitored and predicted reliability indices
4.2 工程实例2计算分析

采用美国I-39北桥第二跨横梁83 d的极值应力监测数据对所提算法进行验证.极值应力数据如图 10所示.

图 10 极值应力监测值(工程实例2) Fig.10 Monitored extreme stress data

首先基于前50 d的监测极值应力数据,结合式(1)~(3)建立动态模型如下:

状态方程为

$ X_{t}=X_{t-1}-0.0563+w_{t}, w_{t} \sim N\left(0, W_{t}\right) $ (27)

监测方程为

$ Z_{t}=X_{t}+v_{t}, v_{t} \sim N(0, 8.814) $ (28)

初始状态信息为

$ X_{t-1} | D_{t-1} \sim N(24.3351, 7.1387) $ (29)

然后采用所提IPF预测算法,结合式(7)~(15)与式(27)~(29),基于第50天到第82天的监测极值应力信息,对第51天到第83天的极值应力进行动态预测,预测结果如图 11所示.并与PF[21]及ARIMA[22]进行比较,预测结果如图 12所示.利用式(15)得到IPF、PF以及ARIMA(0, 0, 0)的EMSE分别为18.897 7、37.612 7、27.238 3,由此可见IPF准确性较好.

图 11 基于改进粒子滤波器的极值应力预测值与实测值(工程实例2) Fig.11 Monitored and predicted extreme stresses
图 12 3种算法极值应力预测结果比较(工程实例2) Fig.12 Comparison of three algorithms

结合式(20)~(21),可得可靠度指标预测值和监测值如图 13所示.可以看出, 可靠度指标预测值较好反映了监测可靠度指标的变化趋势,而考虑了不确定性与随机性的可靠度指标预测值略低于实测值,但由于考虑了不确定性,预测值更符合工程实际.

图 13 可靠度指标的实测值与一步预测值(工程实例2) Fig.13 Monitored and predicted reliability indices
4.3 工程实例3计算分析

采用吉林长春伊通河桥上游主梁70 h的极值挠度监测数据对所提算法进行验证.极值挠度数据如图 14所示.

图 14 极值挠度监测值(工程实例3) Fig.14 Monitored extreme deflection data

首先基于前50 h的监测极值挠度数据,结合式(1)~(3)建立动态模型如下:

状态方程为

$ X_{t}=X_{t-1}-0.0495+w_{t}, w_{t} \sim N\left(0, W_{t}\right) $ (30)

监测方程为

$ Z_{t}=X_{t}+v_{t}, v_{t} \sim N(0, 1.4213) $ (31)

初始状态信息为

$ X_{t-1} | D_{t-1} \sim N(45.8692, 0.9539) $ (32)

然后采用所提IPF预测算法,结合式(7)~(15)与式(30)~(32),基于第50小时到第69小时的监测极值挠度信息,对第51小时到第70小时的极值挠度进行动态预测,预测结果如图 15所示.并与PF[21]及ARIMA[22]进行比较,预测结果对比如图 16所示.利用式(15)得到IPF、PF以及ARIMA(1, 0, 2)的EMSE分别为0.907 6、1.628 3、1.427 7,由此可见IPF预测性能较好.

图 15 基于改进粒子滤波器的极值挠度预测值与实测值(工程实例3) Fig.15 Monitored and predicted extreme deflections based on IPF
图 16 3种算法极值挠度预测结果比较(工程实例3) Fig.16 Comparison of three algorithms

结合式(22)和式(23),可得可靠度指标分析结果如图 17所示.可以看出可靠度指标预测值很好地反映了监测可靠度指标的变化趋势,而考虑了不确定性与随机性的可靠度指标预测值略低于实测值,但由于考虑了不确定性,预测结果更符合工程实际.

图 17 可靠度指标的实测值与一步预测值(工程实例3) Fig.17 Monitored and predicted reliability indices
4.4 桥梁环境突变的假设试验

由于桥梁的应力环境受外界因素影响较大,当在役桥梁气候条件、桥下水位等外在环境因素发生突变时,桥梁的应力状况会也会随之出现变化.假设4.1中富民桥在监测到第600小时候桥梁环境发生突变,例如受到洪水影响,现假定其实际监测信息受到如下变化:

$ \widetilde{Z}_{t}=Z_{t}+28-\frac{(t-700)^{2}}{1000}+\varepsilon_{t}, \varepsilon_{t} \sim N(0, 1) $ (33)

式中: $\widetilde{Z}_{t}, t \geqslant 600 $表示从第600小时之后监测到的极值应力信息;用$ 28-\frac{(t-700)^{2}}{1000}$近似表示受到一个在100 h以内逐渐增长而在100 h后逐渐降低的水流冲击效应影响;εt为其他随机扰动.

现利用改进的粒子滤波算法,结合式(7)~(15)、式(24)~(26)以及式(33),基于第300小时到第809小时的监测极值应力数据,对第301小时到第810小时的监测极值应力进行实时动态一步预测,同时与PF方法结果进行对比,预测结果如图 18所示.可以明显地看出,由于IPF的建议分布在更新过程中借助了突变后的应力信息,因此仍能产生较为合理的建议分布,有效地对突变后的监测信息进行预测,而PF算法由于粒子集合的分布不能很好地覆盖真实值,因此当监测信息突变时出现了难以避免的的滤波发散.

图 18 桥梁环境突变下的应力预测结果 Fig.18 Bridge stress prediction under catastrophic conditions
5 结论

提出了一种用于桥梁动态可靠度指标在线预测的改进粒子滤波算法.基于天津富民桥主缆、美国I-39北桥第二跨横桥以及伊通河桥上游主梁的监测数据对所提算法进行了验证分析,为证实算法的适用性和鲁棒性,进行了假设试验验证,通过假设试验可以模拟多种情况下桥梁极值应力(或挠度)的预测分析.

基于监测极值应力(或挠度)信息对极值应力(或挠度)进行了预测分析,预测结果与实测值变化趋势基本保持一致,预测结果的均方误差也较为理想,且通过对比发现改进的粒子滤波器预测精度十分可观;此外,在桥梁极端情况下,IPF的预测性能良好,具有良好的预测稳定性以及对突变数据的灵敏性与自适应性.而对桥梁结构可靠性进行预测分析时,可靠度指标的一步预测值与实测值标化趋势十分吻合,而考虑了随机性影响所得到的预测可靠度指标从整体上略微低于实测值,但更符合实际情况.

参考文献
[1]
BAO Y Q, BECK J L, LI H. Compressive sampling for accelerometer signals in structural health monitoring[J]. Structural Health Monitoring: An International Journal, 2011, 10(3): 235 DOI:10.1177/1475921710373287
[2]
ZHANG Y, LI J. Wavelet-based vibration sensor data compression technique for civil infrastructure condition monitoring[J]. Journal of Computing in Civil Engineering, 2006, 20(6): 390 DOI:10.1061/(ASCE)0887-3801(2006)20:6(390)
[3]
TROPP J A, GILBERT A C. Signal recovery from random measurements via orthogonal matching pursuit[J]. IEEE Transactions on Information Theory, 2007, 53(12): 4655 DOI:10.1109/TIT.2007.909108
[4]
李惠, 周文松, 欧进萍, 等. 大型桥梁结构智能健康监测系统集成技术研究[J]. 土木工程学报, 2006, 39(2): 46
LI Hui, ZHOU Wensong, OU Jinping, et al. A study on system integration technique of intelligent monitoring systems for soundness of long-span bridges[J]. China Civil Engineering Journal, 2006, 39(2): 46 DOI:10.3321/j.issn:1000-131X.2006.02.010
[5]
李爱群, 缪长青, 李兆霞, 等. 润扬长江大桥结构健康监测系统研究[J]. 东南大学学报, 2003, 33(5): 544
LI Aiqun, MIU Changqing, LI Zhaoxia, et al. Health monitoring system for the Runyang Yangtse River bridge[J]. Journal of Southeast University, 2003, 33(5): 544 DOI:10.3321/j.issn:1001-0505.2003.05.003
[6]
宗周红, 褚福鹏, 牛杰. 基于响应面模型修正的桥梁结构损伤识别方法[J]. 土木工程学报, 2013, 46(2): 115
ZONG Zhouhong, CHU Fupeng, NIU Jie. Damage identification methods of bridge structures using response surface based on finite element model updating[J]. China Civil Engineering Journal, 2013, 46(2): 115
[7]
宗周红, 周儒勉, 郑沛娟. 基于健康监测的桥梁结构损伤预后和安全预后研究进展及挑战[J]. 中国公路学报, 2014, 27(12): 46
ZONG Zhouhong, ZHOU Rumian, ZHENG Peijuan. Damage and safety prognosis of bridge structures based on structural health monitoring: progress and challenges[J]. China Journal of Highway and Transport, 2014, 27(12): 46 DOI:10.3969/j.issn.1001-7372.2014.12.006
[8]
魏锦辉, 任伟新. 基于响应面方法的桥梁静动力有限元模型修正[J]. 公路交通科技, 2015, 32(2): 68
WEI Jinhui, REN Weixin. Static and dynamic bridge finite element model updating based on response surface method[J]. Journal of Highway and Transportation Research and Development, 2015, 32(2): 68 DOI:10.3969/j.issn.1002-0268.2015.02.011
[9]
FRANGOPOL D M, STRAUSS A, KIM S Y. Use of monitoring extreme data for the performance prediction of structures: general approach[J]. Engineering Structures, 2008, 30: 3644 DOI:10.1016/j.engstruct.2008.06.010
[10]
STRAUSS A, FRANGOPOL D M, KIM S Y. Use of monitoring extreme data for the performance prediction of structures: Bayesian updating[J]. Engineering Structures, 2008, 30: 3654 DOI:10.1016/j.engstruct.2008.06.009
[11]
赵卓.基于ARMA模型的伊通河桥监测数据建模与可靠度分析[D].哈尔滨: 哈尔滨工业大学, 2012.
ZHAO Zhuo. Health monitoring data modeling and reliability analysis for Yitong River Bridge based on ARMA model [D]. Harbin: Harbin Institute of Technology, 2012. http://cdmd.cnki.com.cn/Article/CDMD-10213-1013037996.htm
[12]
樊学平.基于验证荷载和监测数据的桥梁可靠性修正与贝叶斯预测[D].哈尔滨: 哈尔滨工业大学, 2014.
FAN Xueping. Bridge reliability updating and Bayesian prediction based on proof loads and monitored data [D]. Harbin: Harbin Institute of Technology, 2014. http://cdmd.cnki.com.cn/Article/CDMD-10213-1014085198.htm
[13]
FAN X P. Bridge extreme stress prediction based on Bayesian dynamic linear models and non-uniform sampling[J]. Structural Health Monitoring: An International Journal, 2017, 16(3): 253 DOI:10.1177/1475921716688166
[14]
李桂芳, 孙勇成, 黄圣国. 粒子滤波理论综述[J]. 控制工程, 2010, 17(S2): 92
LI Guifang, SUN Yongcheng, HUANG Shengguo. A survey of particle filter theory[J]. Control Engineering of China, 2010, 17(S2): 92
[15]
樊学平, 刘月飞, 吕大刚. 应用高斯粒子滤波器的桥梁可靠性在线预测[J]. 哈尔滨工业大学学报, 2016, 48(6): 164
FAN Xueping, LIU Yuefei, LÜ Dagang. On-line reliability prediction of bridges based on gaussian particle filter[J]. Journal of Harbin Institute of Technology, 2016, 48(6): 164
[16]
樊学平, 刘月飞, 吕大刚. 桥梁极值应力的改进高斯混合粒子滤波器动态预测[J]. 同济大学学报(自然科学版), 2016, 44(11): 1660
FAN Xueping, LIU Yuefei, LÜ Dagang. Improved gaussian mixed particle filter dynamic prediction of bridge monitored extreme stress[J]. Journal of Tongji University (Natural Science), 2016, 44(11): 1660 DOI:10.11908/j.issn.0253-374x.2016.11.003
[17]
FAN X P, LIU Y F. Dynamic extreme stress prediction of bridges based on nonlinear mixed gaussian particle filtering algorithm and structural health monitoring data[J]. Advances in Mechanical Engineering, 2016, 8(6): 1
[18]
XIONG X, NAVON I M, UZUNOGLU B. A note on the particle filter with posterior gaussian resampling[J]. Dynamic Meteorology and Oceanography, 2006, 58(4): 456 DOI:10.1111/j.1600-0870.2006.00185.x
[19]
陈莹, 郑金荣. 基于观测迭代的插值粒子滤波算法[J]. 信号处理, 2010, 26(1): 146
CHEN Ying, ZHENG Jinrong. Improved divided difference particle filter based on observation interation[J]. Signal Processing, 2010, 26(1): 146 DOI:10.3969/j.issn.1003-0530.2010.01.028
[20]
于金霞, 刘文静, 汤永利. 粒子滤波重采样算法研究[J]. 微计算机信息, 2010, 26(16): 44
YU Jinxia, LIU Wenjing, TANG Yongli. Study on the resampling algorithms of particle filter[J]. Microcomputer Information, 2010, 26(16): 44 DOI:10.3969/j.issn.2095-6835.2010.16.017
[21]
郭晓松, 李奕芃, 郭君斌. 粒子滤波算法及其应用研究[J]. 计算机工程与设计, 2009, 30(9): 2264
GUO Xiaosong, LI Yipeng, GUO Junbin. Research on particle filter algorithm and its applications[J]. Computer Engineering & Design, 2009, 30(9): 2264
[22]
DHOTE V, SATANANDMISHRA, SHUKLA J P, et al. Runoff prediction using big data analytics based on ARIMA model[J]. Indian Journal of Geo-marine Sciences, 2018, 47(11): 2163