CN115017767A - 基于贝叶斯正则化的桥梁影响线识别与不确定性量化方法 - Google Patents
基于贝叶斯正则化的桥梁影响线识别与不确定性量化方法 Download PDFInfo
- Publication number
- CN115017767A CN115017767A CN202210621432.9A CN202210621432A CN115017767A CN 115017767 A CN115017767 A CN 115017767A CN 202210621432 A CN202210621432 A CN 202210621432A CN 115017767 A CN115017767 A CN 115017767A
- Authority
- CN
- China
- Prior art keywords
- bridge
- influence line
- bridge influence
- identification
- function
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 72
- 238000011002 quantification Methods 0.000 title claims abstract description 10
- 241001244708 Moroccan pepper virus Species 0.000 claims abstract description 50
- 238000009826 distribution Methods 0.000 claims abstract description 23
- 238000005457 optimization Methods 0.000 claims abstract description 11
- 230000006870 function Effects 0.000 claims description 64
- 230000004044 response Effects 0.000 claims description 31
- 239000011159 matrix material Substances 0.000 claims description 23
- 238000004422 calculation algorithm Methods 0.000 claims description 8
- 238000011156 evaluation Methods 0.000 claims description 7
- 238000005259 measurement Methods 0.000 claims description 7
- 230000008569 process Effects 0.000 claims description 7
- 239000013598 vector Substances 0.000 claims description 7
- 238000005315 distribution function Methods 0.000 claims description 6
- 238000004364 calculation method Methods 0.000 claims description 5
- 230000000694 effects Effects 0.000 claims description 4
- 239000000126 substance Substances 0.000 claims description 3
- 238000013139 quantization Methods 0.000 claims 5
- 238000010586 diagram Methods 0.000 description 12
- 238000012544 monitoring process Methods 0.000 description 8
- 238000011160 research Methods 0.000 description 5
- 230000036541 health Effects 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000000354 decomposition reaction Methods 0.000 description 3
- 238000013179 statistical model Methods 0.000 description 3
- 238000013398 bayesian method Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000013461 design Methods 0.000 description 2
- 238000013076 uncertainty analysis Methods 0.000 description 2
- 238000007476 Maximum Likelihood Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000010606 normalization Methods 0.000 description 1
- 230000010355 oscillation Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
- 238000010561 standard procedure Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/08—Probabilistic or stochastic CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/10—Numerical modelling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/10—Noise analysis or noise optimisation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T10/00—Road transport of goods or passengers
- Y02T10/10—Internal combustion engine [ICE] based vehicles
- Y02T10/40—Engine management systems
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)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Complex Calculations (AREA)
Abstract
本发明提出的基于贝叶斯正则化的桥梁影响线识别与不确定性量化方法,包括:基于桥梁影响线识别模型建立贝叶斯系统识别中的预测误差;利用预测误差的分布特性,构造似然函数和模型参数先验概率;根据似然函数和模型参数先验概率,基于贝叶斯层次模型建立后验概率密度函数,且通过后验概率密度函数,将桥梁影响线识别转化为目标函数优化;计算桥梁影响线的最优值(MPVs);量化桥梁影响线最优值的后验不确定性信息;计算标准化置信区间指标,评估不同类型桥梁影响线之间的识别质量,本发明方法可以快速、准确、自动地确定正则化系数进行桥梁影响线的识别,并且该方法能够不依赖基准影响线,仅通过影响线识别的后验不确定性信息定性评价识别质量。
Description
技术领域
本发明涉及工程结构健康监测领域,特别是指基于贝叶斯正则化的桥梁影响线识别与不确定性量化方法。
背景技术
桥梁影响线是桥梁的固有特性,基于桥梁这一特性在桥梁状态的评估和损伤识别等领域的应用具有非常好的前景。从桥梁监测数据中识别桥梁影响线一直是一个重要的研究领域。作为一个反问题,桥梁影响线识别的基本原理是从预先构建的车辆荷载矩阵和桥梁响应向量中重构桥梁影响线。有研究提出了桥梁影响线反演计算的思想,并使用最小二乘法,利用重型车辆通过简支梁桥时的应变响应,对其进行求解。然而,最小二乘法通常对测量噪声和误差非常敏感,这不可避免地引起桥梁影响线识别结果较大的振荡。因此,桥梁影响线识别通常是一个不适定的反问题,容易受到噪声的影响。为了解决该方法的病态性,通常需要正则化来求解病态反问题,以获得稳定的有界解。正则化可以是一种物理约束,比如模型参数必须是非负的,或者只是一种数学折衷。在过去几年中,几种正则化方法已被成功用于桥梁影响线识别。例如,通过引入惩罚项作为桥梁影响线的先验知识,提出了Tikhonov正则化模型。为了使其广泛适用于各种形状的桥梁影响线,提出了一种结合稀疏正则化和自适应B样条基函数的桥梁影响线识别模型,以及提出了一种基于正则化的最小二乘QR分解方法,通过整合经验模态分解来处理高速车辆引起的结构动力响应。其他先进技术也被用于改进影响线识别方法,例如结合计算机视觉、概率统计模型(最大似然估计)和频域方法。
但是,上述方法均属于确定性方法,存在自身的局限性:(1)不确定性的影响广泛地存在于影响线识别问题中,桥梁影响线多次识别结果存在不同程度的随机离散现象,确定性方法无法给出科学合理的解释。(2)为了能够客观评价影响线识别的准确性,现有的确定性方法通常需要对比桥梁基准影响线,事实上,桥梁基准影响线信息普遍缺失或者不够准确,例如成桥之初实测的基准影响线不一定能准确反映当前状况。(3)基于正则化的影响线识别方法,正则化系数的选择在一定程度决定着影响线识别的准确性。虽可通过引入辅助的数学模型优化正则化系数,但寻优的过程较耗时,大多采用系数矩阵的奇异值分解,计算效率会随着矩阵维数加大而显著降低。需要进行人工干预,结果好坏受模型使用者经验的影响。因此,需要建立影响线识别的概率方法,从而克服确定性方法的局限性。
概率方法通常分为两类:经典概率理论(即频率学派)与贝叶斯理论(即贝叶斯学派)。其中,贝叶斯学派和频率学派对于在应用中概率如何被赋值有着截然不同的看法,贝叶斯学派完全同意概率的公理化定义,但认为概率也可以用经验确定,可以根据逻辑学中某人对一个未知命题信任的程度来赋值概率,因而更符合实际情况,扩大了研究范围。这也是贝叶斯系统识别框架能够广泛应用于模型分析、损伤检测和结构模型更新等领域的原因。贝叶斯框架内的正则化技术能够使用MCMC技术或变分贝叶斯类方法自动确定最优正则化系数以及量化模型参数的后验不确定性信息。从不确定度的角度来看,统计度量可以作为评估解决方案质量的有用信息。近年来,贝叶斯正则化已成功应用于解决各个领域的关键问题,如荷载识别、图像处理和模型更新等。作为贝叶斯统计学理论的一个分支,贝叶斯正则化可以基于先验分布实现对不同正则化模型的概率解释。当先验分布为高斯分布时,贝叶斯正则化可以像Tikhonov正则化一样通过平衡原始目标函数和解的平滑条件来处理病态逆问题。近年来,基于贝叶斯统计学理论的应用越来越广泛,模型参数的识别已经非常成熟。例如,通过贝叶斯正则化实现了冲击荷载的识别,在识别的过程中应用序贯贝叶斯迭代算法自动确定正则化系数,大大提高了模型参数识别的准确性和鲁棒性。此外,不确定性分析方法的研究可以很大程度上提高结构动力分析方法的准确性和鲁棒性,当然,结构健康监测系统在应用不确定性分析时也可以得到非常可靠的分析依据。到目前为止,还鲜有研究提出基于贝叶斯框架的桥梁影响线识别方法。
发明内容
本发明的主要目的在于克服现有技术中的上述缺陷,提出一种基于贝叶斯正则化的桥梁影响线识别与不确定性量化方法,该方法可以快速、准确、自动地确定正则化系数进行桥梁影响线的识别,并且该方法能够不依赖基准影响线,仅通过影响线识别的后验不确定性信息定性评价识别质量。因此,该方法在结构健康监测领域具有广阔的应用前景。
本发明采用如下技术方案:
基于贝叶斯正则化的桥梁影响线识别与不确定性量化方法,包括:
基于桥梁影响线识别模型建立贝叶斯系统识别中的预测误差;
利用预测误差的分布特性,构造似然函数和模型参数先验概率;
根据似然函数和模型参数先验概率,基于贝叶斯层次模型建立后验概率密度函数,且通过后验概率密度函数,将桥梁影响线识别转化为目标函数优化;
计算桥梁影响线的最优值;
量化桥梁影响线最优值的后验不确定性信息;
计算标准化置信区间指标,评估不同类型桥梁影响线之间的识别质量。
具体地,所述基于桥梁影响线识别模型建立贝叶斯系统识别中的预测误差,具体为:
设车辆沿着给定的直线在桥梁上行驶,且给定车辆的各个轴独立作用于桥梁上,则桥梁响应近似为所有单轴荷载作用的总和:
Rs=AX
其中,Rs是特定位置的准静态响应向量;A是基于车辆信息构建的荷载矩阵;X表示桥梁影响线系数向量;
采用分段多项式插值函数拟合桥梁影响线:
X=Φf=Ωc
其中Ω为拟合矩阵;c为插值函数系数,Rs=AX改写为:
Rs=AΩc
将桥梁影响线识别模型嵌入贝叶斯框架,建立桥梁测量响应Rm与准静态响应Rs之间的预测误差τ,可如下所示:
τ=Rm-AΩc
具体地,利用预测误差的分布特性,构造似然函数和模型参数先验概率,具体为:
设定预测误差τ服从零均值的高斯分布,似然函数可以表示为:
其中是σ测量数据中包含的噪声和误差的标准差,Ng是响应采样点的数量,c为插值函数系数;
先验概率分布函数表示如下:
式中,μ2是表征插值函数系数c变化的尺度方差,Ns是插值函数系数的数量,参数σ2和μ2的概率密度函数由共轭先验原理确定,并由逆伽马分布表示:
其中(α0,β0)和(α1,β1)是逆伽马分布的超参数。
具体地,根据似然函数和模型参数先验概率,基于贝叶斯层次模型建立后验概率密度函数,且通过后验概率密度函数,将桥梁影响线识别转化为目标函数优化;具体为:
根据似然函数和先验分布,建立贝叶斯层次模型,构造系统的后验概率密度函数:
P(c,σ2,μ2|Rm)∝P(Rm|c,σ2)P(σ2)P(c|μ2)P(μ2)
进一步得到后验概率密度函数P(c,σ2,μ2|Rm)的表达式:
其中,||·||表示欧几里得范数;
基于后验概率分布函数,将桥梁影响线识别问题转化为目标函数优化;
P(c,σ2,μ2|Rm)的负对数如下所示:
具体地,计算桥梁影响线的最优值,具体为:
对如下公式进行最小化,
得到c,σ和μ的一阶导数,令c,σ和μ的一阶导数均为0,获得最优值方程组:
σ2=(2(α0+1)+Ng)-1(||AΩc-Rm||2+2β0)
μ2=(2(α1+1)+Ns)-1(||c||2+2β1)
具体地,所述量化桥梁影响线最优值的后验不确定性信息,具体为:
使用L对c,σ和μ的二阶导数获得Hessian矩阵:
取L的二阶导数,得到c的Hessian矩阵:
桥梁影响线MPVs的Hessian矩阵HXX和Hcc之间的关系可以表示为:
进一步代入,并对HXX求逆,进而量化桥梁影响线的后验不确定性信息:
具体地,所述计算标准化置信区间指标,评估不同类型桥梁影响线之间的识别质量,具体为:
识别的桥梁影响线的中心线(CL)、上限(UB)和下限(LB)可定义为:
CL=MPVs
UB,LB=CL±Zα2·σ'
式中,CL为桥梁影响线的中心线,UB为桥梁影响线的上限,LB为桥梁影响线的下限,MPVs为桥梁影响线的最优值,σ'是已识别桥梁影响线的标准偏差,Zα2是对应于可信水平α的可信系数,标准偏差σ'由等式(16)中的协方差矩阵确定;置信区间宽度WCI被定义为评估指标:
WCI=UB-LB=2Zα2·σ'
对WCI进行标准化,提出标准化置信区间SWCI指标:
其中
RMPVs=max(MPVs)-min(MPVs)
式中max(MPVS),min(MPVS)和RMPVs分别为已识别MPVs的最大值、最小值和范围。
由上述对本发明的描述可知,与现有技术相比,本发明具有如下有益效果:
本发明提出的基于贝叶斯正则化的桥梁影响线识别与不确定性量化方法,包括:基于桥梁影响线识别模型建立贝叶斯系统识别中的预测误差;利用预测误差的分布特性,构造似然函数和模型参数先验概率;根据似然函数和模型参数先验概率,基于贝叶斯层次模型建立后验概率密度函数,且通过后验概率密度函数,将桥梁影响线识别转化为目标函数优化;计算桥梁影响线的最优值;量化桥梁影响线最优值的后验不确定性信息;计算标准化置信区间指标,评估不同类型桥梁影响线之间的识别质量;本发明方法同确定性类方法相比,基于统计学模型识别桥梁影响线,能够快速、准确、自动地确定正则化系数进行桥梁影响线识别,并具有很强的鲁棒性和很高的准确性,相比现有的确定性类方法,该方法可以通过量化桥梁影响线MPVs的后验不确定性信息,不依赖基准影响线实现识别结果的定性评价,在结构健康监测领域中具有广阔的应用前景。
附图说明
图1为本发明实施例提供的三跨连续梁示意图;
图2为本发明实施例提供的考虑不同噪声干扰的挠度响应示意图,(a)SNR=30dB;(b)SNR=20dB;(c)SNR=10dB;
图3为本发明实施例提供的B点桥梁挠度影响线识别结果示意图;
图4为本发明实施例提供的A点挠度响应(SNR=20dB)示意图;
图5为本发明实施例提供的A点影响线MPVs置信区间示意图,(a)SNR=30dB;(b)SNR=20dB;(c)SNR=10dB;
图6为本发明实施例提供的指标变化趋势示意图,(a)为WCI指标变化趋势;(b)为SWCI指标变化趋势;
图7为本发明实施例提供的实桥图片;
图8为本发明实施例提供的传感器布置示意图,(a)为连续箱梁桥及传感器布置简图,(b)为传感器实际布置图;
图9为本发明实施例提供的快速公交车辆实车图片;
图10为本发明实施例提供的桥梁监测系统界面示意图;
图11为本发明实施例提供的桥梁应变响应示意图,(a)为1/4跨(L1),(b)为1/2跨(L2),(c)3/4跨(L3);
图12为本发明实施例提供的桥梁应变影响线识别结果示意图,(a)为1/4跨(L1),(b)为1/2跨(L2),(c)3/4跨(L3);
图13为本发明实施例提供的桥梁应变影响线MPVs的置信区间示意图,(a)为1/4跨(L1),(b)为1/2跨(L2),(c)3/4跨(L3);
图14为本发明实施例提供的MPVs的最大值、最小值和范围。
以下结合附图和具体实施例对本发明作进一步详述。
具体实施方式
本发明要解决的技术问题,就是提供一种基于贝叶斯正则化的桥梁影响线识别与不确定性量化方法。基于贝叶斯统计学理论将桥梁影响线识别模型嵌入贝叶斯框架,进而提出一种桥梁影响线识别的统计学方法。
解决上述技术问题,本发明采用的技术方案包含以下六个步骤:(1)基于桥梁影响线识别模型建立贝叶斯系统识别理论中的预测误差;(2)似然函数和模型参数先验概率的构造;(3)基于贝叶斯层次模型建立后验概率密度函数;(4)计算桥梁影响线的最优值(MPVs);(5)量化桥梁影响线MPVs的后验不确定性信息;(6)计算标准化置信区间(SWCI)指标。具体如下:
步骤S1,基于桥梁影响线识别模型建立贝叶斯系统识别中的预测误差;
假设车辆沿着给定的直线在桥梁上行驶,且给定车辆的各个轴独立作用于桥梁上,则桥梁响应可近似为所有单轴荷载作用的总和:
Rs=AX (1)
其中,Rs是特定位置的准静态响应向量;A是基于车辆信息(如轴重和轴距)构建的荷载矩阵;X表示桥梁影响线系数向量。
通过分析影响线的物理意义和理论推导可知,桥梁影响线具有光滑、连续、高阶可导的曲线性质。因此,拟采用分段多项式插值函数拟合影响线:
X=Φf=Ωc (2)
其中Ω为拟合矩阵;c为插值函数系数。式(1)可以改写为:
Rs=AΩc (3)
为了将桥梁影响线识别模型嵌入贝叶斯框架,建立桥梁测量响应Rm与准静态响应Rs之间的预测误差τ,可如下所示:
τ=Rm-AΩc (4)
步骤S2,根据似然函数和模型参数先验概率,基于贝叶斯层次模型建立后验概率密度函数,且通过后验概率密度函数,将桥梁影响线识别转化为目标函数优化;
预测误差τ的不确定性可以根据最大信息熵原理,用近似的概率分布的形式表示,一般地,假定预测误差τ服从零均值的高斯分布。因此,似然函数可以表示为:
其中是σ2测量数据中包含的噪声和误差的方差,Ng是响应采样点的数量。先验概率分布函数反映了研究者对待确定模型参数的先验知识。先验的一个普遍选择是多元正态分布。通常由标准技术根据共轭先验原理确定,表示如下:
式中,μ2是表征插值函数系数c变化的尺度方差,Ns是插值函数系数的数量。参数σ2和μ2的概率密度函数由共轭先验原理确定,并由逆伽马分布表示:
其中(α0,β0)和(α1,β1)是逆伽马分布的超参数。
步骤S3,基于贝叶斯层次模型建立后验概率密度函数
在贝叶斯统计学理论中,后验概率密度函数是由似然函数和先验分布通过贝叶斯定理得到的。建立贝叶斯层次模型,构造系统的后验概率密度函数:
P(c,σ2,μ2|Rm)∝P(Rm|c,σ2)P(σ2)P(c|μ2)P(μ2) (9)
通过将等式(5)到等式(8)代入等式(9),得到后验概率密度函数P(c,σ2,μ2|Rm)的表达式:
其中,||·||表示欧几里得范数。
基于后验概率分布函数,将桥梁影响线识别问题转化为目标函数优化问题。P(c,σ2,μ2|Rm)的负对数如下所示:
步骤S4,计算桥梁影响线的最优值
最小化式(11),可得式(11)关于c,σ和μ的一阶导数,让一阶导数为0,以获得最优值方程组:
σ2=(2(α0+1)+Ng)-1(||AΩc-Rm||2+2β0) (12-2)
μ2=(2(α1+1)+Ns)-1(||c||2+2β1) (12-3)
(1).给定基于先验知识的已知参数值,逆伽马分布中的超参数初始值进行下述假定:
步骤S5,量化桥梁影响线MPVs的后验不确定性信息
使用L对c,σ和μ的二阶导数获得Hessian矩阵:
取L的二阶导数,得到c的Hessian矩阵:
桥梁影响线MPVs的Hessian矩阵HXX和Hcc之间的关系可以表示为:
将式(12-2)和式(12-3)代入式(15),对HXX求逆,进而量化桥梁影响线的后验不确定性信息:
步骤S6,计算标准化置信区间指标,评估不同类型桥梁影响线之间的识别质量已识别桥梁影响线(如图1所示)的中心线(CL)、上限(UB)和下限(LB)可定义为:
CL=MPVs (17)
UB,LB=CL±Zα/2·σ' (18)
式中,σ'是已识别桥梁影响线的标准偏差,Zα/2是对应于可信水平α的可信系数。标准偏差σ'由等式(16)中的协方差矩阵确定。在本专利中,不确定性水平用于定性评估影响线识别质量,置信区间宽度WCI被定义为评估指标:
WCI=UB-LB=2Zα/2·σ' (19)
鉴于来自不同位置的不同类型的桥梁影响线在振幅方面可能存在实质性差异,需要进行标准化处理,进而可以实现不同桥梁影响线之间的比较。对WCI进行标准化,提出标准化置信区间(SWCI)指标:
其中
RMPVs=max(MPVs)-min(MPVs) (21)
式中max(MPVS),min(MPVS)和RMPVs分别为已识别MPVs的最大值、最小值和范围,如图14。值得一提的是,SWCI指标是无量纲的,也可用于评估不同类型桥梁影响线之间的识别质量,例如挠度和应变影响线。
实施算例1:三跨数值连续梁的影响线识别
如图1所示,采用三跨连续梁的数值算例用于验证本发明提出的影响线识别统计学模型。该连续梁由3个等跨梁组成,每跨10m,总长L=30m。横截面为矩形截面1m×0.5m。杨氏模量E=3.25×104MPa,密度ρ=2.25×103kg/m3。利用MATLAB编程建立平面梁单元的有限元模型,沿纵向均匀划分成300个梁单元,每个单元长0.1m。
用于数值模拟的四轴车辆的每轴轴重为2kN,相邻轴间距为0.5m。模拟竖向荷载从梁左边端缓慢移动至右边端的过程,可计算得到从A位置(左边跨跨中)、B位置(中跨跨中)、C位置(右边跨跨中)输出的移动荷载引起的挠度响应。模拟信噪比(SNR)分别为30dB、20dB和10dB的高斯白噪声信号,并加入到挠度响应中,可得位置B加噪后的挠度响应,如图2所示,其中(a)SNR=30dB;(b)SNR=20dB;(c)SNR=10dB。
在B点的响应中加入信噪比为20dB的高斯白噪声。根据车辆信息构造荷载矩阵A。根据序贯贝叶斯学习迭代算法计算桥梁影响线的MPVs,并绘制在图3中。为了进行比较,图中还绘制了使用确定性方法得到的桥梁影响线识别结果。可以看出有贝叶斯正则化方法得到的桥梁影响线识别结果更为准确。
为了验证所提出贝叶斯方法在正则化系数选择和桥梁影响线识别上的鲁棒性,将四种等级的高斯白噪声添加到位置B处的挠度响应中。应用L-curve法、GCV法和贝叶斯正则化方法选取正则化系数,并进行桥梁影响线识别。根据图3和表1中的影响线识别结果,对于5dB到30dB的SNR值,通过贝叶斯正则化获得的相对误差(RE)低于通过L-curve法和GCV法获得的RE。随着噪声水平的增加,基于贝叶斯正则化的RE值增长趋势明显低于其他方法。结果表明,与L-curve法和GCV法相比,本发明所提出的贝叶斯正则化方法对正则化系数的选择具有更好的鲁棒性。而且,在本案例中,传统方法平均需要约2min,而贝叶斯正则化平均在1s内获得最优的正则化系数,具有更高的计算效率。
随后,为了验证本发明提出的贝叶斯正则化方法的定性评价性能,在位置A处的响应中加入SNR值为10dB至30dB的高斯白噪声,SNR=20dB的响应在图4中已展示。根据序贯贝叶斯迭代算法计算桥梁影响线的MPVs,基于式(16)对影响线识别结果的后验不确定性信息进行量化。基于后验标准差,可以得到贝叶斯框架内MPVs的置信区间,如图5所示,其中(a)SNR=30dB;(b)SNR=20dB;(c)SNR=10dB。然后,计算MPVs的SWCI和WCI,如表2和图5所示。
图6(a)和表2中的结果表明,MPVs的WCI的变化趋势与噪声等级和RE值的变化趋势一致,且数值很小,这也表明了本发明所提方法可以高精度识别桥梁影响线。表2和图6(b)中SWCI的变化趋势与RE的变化趋势一致。当SWCI在阈值(0–0.1465)内时,影响线识别结果的质量较高,超过阈值后,质量逐渐下降。结果表明,基于贝叶斯框架内的不确定性信息,可以对影响线识别结果质量进行定性评价。
表1.正则化系数和相对误差(RE)值
表2.挠度影响线MPVs的WCI,SWCI和相对误差
在此验证中,连续梁数值算例的结果表明,该方法可以快速准确地确定正则化系数,进而高效地识别桥梁影响线,具有很强的鲁棒性和很高的准确性。而且,该方法通过量化桥梁影响线MPVs的后验不确定性信息,能够不依赖基准影响线实现识别结果的定性评价。
实施算例2:三跨连续箱梁实桥影响线识别
厦门快速公交2号线(同安线)项目起于滨海大道至同安枢纽,为快速公交车辆提供封闭的专用车道。真实桥梁如图7所示,该桥梁为三跨连续箱梁桥,各跨跨长均为30m。快速公交系统的桥梁、传感器布置简图以及实际布局如图8所示,在中跨的1/4处(L1)、1/2处(L2)以及3/4(L3)处分别布置了应变传感器,用于采集桥梁的应变响应数据,并上传至路网监测系统,路网监测系统界面如图9所示。快速公交车辆实车如图10所示。
通过现场调研获得的快速公交车辆信息如表3所示。为了消除乘客对荷载信息的影响,选择低峰时段的桥梁响应进行影响线路识别。应变响应如图11所示,其中(a)为1/4跨(L1),(b)为1/2跨(L2),(c)3/4跨(L3)。
为了验证本发明提出的贝叶斯正则化方法。首先,基于表3的车辆信息构造荷载矩阵A,并联合应变响应的监测数据作为模型输入,应用S4中的算法计算桥梁影响线MPVs,如图12所示,其中,(a)为1/4跨(L1),(b)为1/2跨(L2),(c)3/4跨(L3);然后,再基于式(16)量化桥梁影响线MPVs的后验不确定性信息,如图13所示,其中(a)为1/4跨(L1),(b)为1/2跨(L2),(c)3/4跨(L3);最后,基于式(20)计算SWCI指标,如表4所示。
表3.快速公交车辆信息
表4.应变影响线MPVs的WCI、SWCI和相对误差
图12中桥梁影响线MPVs与参考影响线的拟合效果表明本发明所提的方法具有很高的准确性,表4中三个位置MPVs的相对误差(RE)均低于5%,这也进一步验证了其准确性。
图13中参考影响线始终处于基于后验不确定性信息计算得到的WCI内,这也进一步验证了本发明所提出的贝叶斯正则化方法能够高精度识别桥梁影响线。表4的SWCI指标以RE为参考,其变化趋势与RE的变化趋势保持一致,因此,该点发现验证了本发明提出的贝叶斯正则化方法基于SWCI指标可以不依赖基准影响线定性评价识别结果的质量。
上述仅为本发明的具体实施方式,但本发明的设计构思并不局限于此,凡利用此构思对本发明进行非实质性的改动,均应属于侵犯本发明保护范围的行为。
Claims (7)
1.基于贝叶斯正则化的桥梁影响线识别与不确定性量化方法,其特征在于,包括:
基于桥梁影响线识别模型建立贝叶斯系统识别中的预测误差;
利用预测误差的分布特性,构造似然函数和模型参数先验概率;
根据似然函数和模型参数先验概率,基于贝叶斯层次模型建立后验概率密度函数,且通过后验概率密度函数,将桥梁影响线识别转化为目标函数优化;
计算桥梁影响线的最优值;
量化桥梁影响线最优值的后验不确定性信息;
计算标准化置信区间指标,评估不同类型桥梁影响线之间的识别质量。
2.根据权利要求1所述的基于贝叶斯正则化的桥梁影响线识别与不确定性量化方法,其特征在于,所述基于桥梁影响线识别模型建立贝叶斯系统识别中的预测误差,具体为:
设车辆沿着给定的直线在桥梁上行驶,且给定车辆的各个轴独立作用于桥梁上,则桥梁响应近似为所有单轴荷载作用的总和:
Rs=AX
其中,Rs是特定位置的准静态响应向量;A是基于车辆信息构建的荷载矩阵;X表示桥梁影响线系数向量;
采用分段多项式插值函数拟合桥梁影响线:
X=Φf=Ωc
其中Ω为拟合矩阵;c为插值函数系数,Rs=AX改写为:
Rs=AΩc
将桥梁影响线识别模型嵌入贝叶斯框架,建立桥梁测量响应Rm与准静态响应Rs之间的预测误差τ,可如下所示:
τ=Rm-AΩc
7.根据权利要求1所述的基于贝叶斯正则化的桥梁影响线识别与不确定性量化方法,其特征在于,所述计算标准化置信区间指标,评估不同类型桥梁影响线之间的识别质量,具体为:
识别的桥梁影响线的中心线(CL)、上限(UB)和下限(LB)可定义为:
CL=MPVs
UB,LB=CL±Zα2·σ'
式中,CL为桥梁影响线的中心线,UB为桥梁影响线的上限,LB为桥梁影响线的下限,MPVs为桥梁影响线的最优值,σ'是已识别桥梁影响线的标准偏差,Zα2是对应于可信水平α的可信系数,标准偏差σ'由等式(16)中的协方差矩阵确定;置信区间宽度WCI被定义为评估指标:
WCI=UB-LB=2Zα2·σ'
对WCI进行标准化,提出标准化置信区间SWCI指标:
其中
RMPVs=max(MPVs)-min(MPVs)
式中max(MPVS),min(MPVS)和RMPVs分别为已识别MPVs的最大值、最小值和范围。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210621432.9A CN115017767A (zh) | 2022-06-02 | 2022-06-02 | 基于贝叶斯正则化的桥梁影响线识别与不确定性量化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210621432.9A CN115017767A (zh) | 2022-06-02 | 2022-06-02 | 基于贝叶斯正则化的桥梁影响线识别与不确定性量化方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN115017767A true CN115017767A (zh) | 2022-09-06 |
Family
ID=83072088
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210621432.9A Pending CN115017767A (zh) | 2022-06-02 | 2022-06-02 | 基于贝叶斯正则化的桥梁影响线识别与不确定性量化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115017767A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116127277A (zh) * | 2023-04-12 | 2023-05-16 | 武汉工程大学 | 激波流场动态压力测量不确定度评定方法及系统 |
CN117852156A (zh) * | 2024-03-08 | 2024-04-09 | 江苏狄诺尼信息技术有限责任公司 | 一种基于贝叶斯的道路平面线位智能重构的方法及系统 |
Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105973619A (zh) * | 2016-04-27 | 2016-09-28 | 厦门大学 | 结构健康监测系统下基于影响线的桥梁局部损伤识别方法 |
CN107588915A (zh) * | 2017-10-18 | 2018-01-16 | 厦门大学 | 一种桥梁影响线识别方法及系统 |
WO2018013247A2 (en) * | 2016-06-02 | 2018-01-18 | Brown University | Physics informed learning machine |
CN108536971A (zh) * | 2018-04-13 | 2018-09-14 | 广州市建筑科学研究院有限公司 | 一种基于贝叶斯模型的结构损伤识别方法 |
CN108920766A (zh) * | 2018-06-05 | 2018-11-30 | 厦门大学 | 一种基于基函数表示和稀疏正则化的桥梁影响线识别方法 |
CN109540089A (zh) * | 2018-10-16 | 2019-03-29 | 华南理工大学 | 一种基于贝叶斯-克里金模型的桥面高程拟合方法 |
WO2019232737A1 (zh) * | 2018-06-07 | 2019-12-12 | 大连理工大学 | 一种基于迭代法的准静态桥梁影响线识别方法 |
CN111859250A (zh) * | 2020-06-19 | 2020-10-30 | 哈尔滨工业大学(深圳)(哈尔滨工业大学深圳科技创新研究院) | 一种环境及受迫激励下快速模态参数准确性评估方法、系统及存储介质 |
-
2022
- 2022-06-02 CN CN202210621432.9A patent/CN115017767A/zh active Pending
Patent Citations (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105973619A (zh) * | 2016-04-27 | 2016-09-28 | 厦门大学 | 结构健康监测系统下基于影响线的桥梁局部损伤识别方法 |
WO2018013247A2 (en) * | 2016-06-02 | 2018-01-18 | Brown University | Physics informed learning machine |
CN107588915A (zh) * | 2017-10-18 | 2018-01-16 | 厦门大学 | 一种桥梁影响线识别方法及系统 |
CN108536971A (zh) * | 2018-04-13 | 2018-09-14 | 广州市建筑科学研究院有限公司 | 一种基于贝叶斯模型的结构损伤识别方法 |
CN108920766A (zh) * | 2018-06-05 | 2018-11-30 | 厦门大学 | 一种基于基函数表示和稀疏正则化的桥梁影响线识别方法 |
WO2019232737A1 (zh) * | 2018-06-07 | 2019-12-12 | 大连理工大学 | 一种基于迭代法的准静态桥梁影响线识别方法 |
CN109540089A (zh) * | 2018-10-16 | 2019-03-29 | 华南理工大学 | 一种基于贝叶斯-克里金模型的桥面高程拟合方法 |
CN111859250A (zh) * | 2020-06-19 | 2020-10-30 | 哈尔滨工业大学(深圳)(哈尔滨工业大学深圳科技创新研究院) | 一种环境及受迫激励下快速模态参数准确性评估方法、系统及存储介质 |
Non-Patent Citations (1)
Title |
---|
施有志;林四新;叶代成;: "钢箱梁自行车桥的舒适度研究", 世界地震工程, no. 04, 15 December 2019 (2019-12-15) * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116127277A (zh) * | 2023-04-12 | 2023-05-16 | 武汉工程大学 | 激波流场动态压力测量不确定度评定方法及系统 |
CN117852156A (zh) * | 2024-03-08 | 2024-04-09 | 江苏狄诺尼信息技术有限责任公司 | 一种基于贝叶斯的道路平面线位智能重构的方法及系统 |
CN117852156B (zh) * | 2024-03-08 | 2024-05-17 | 江苏狄诺尼信息技术有限责任公司 | 一种基于贝叶斯的道路平面线位智能重构的方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN115017767A (zh) | 基于贝叶斯正则化的桥梁影响线识别与不确定性量化方法 | |
CN106650825B (zh) | 一种机动车尾气排放数据融合系统 | |
CN108960303B (zh) | 一种基于lstm的无人机飞行数据异常检测方法 | |
CN110895878B (zh) | 一种基于ge-gan的交通状态虚拟检测器的生成方法 | |
CN112926265A (zh) | 基于遗传算法优化神经网络的大气多孔探针测量校准方法 | |
Carraro et al. | Weight estimation on static B-WIM algorithms: A comparative study | |
CN111709549A (zh) | 一种基于svd-pso-lstm的短时交通流预测的导航提醒方法 | |
CN113009447B (zh) | 基于深度学习和探地雷达的道路地下空洞检测预警方法 | |
CN108846200B (zh) | 一种基于迭代法的准静态桥梁影响线识别方法 | |
CN116451514A (zh) | 桥梁结构状态评估方法、装置以及电子设备 | |
CN111967308A (zh) | 一种在线路面不平度辨识方法及系统 | |
CN116088453A (zh) | 生产质量预测模型训练方法及装置、生产质量监控方法 | |
CN114510778A (zh) | 基于混合智能优化lstm的轨道不平顺预测方法 | |
CN115545100A (zh) | 基于LSTM的GB-InSAR大气相位补偿方法 | |
CN114595721A (zh) | 一种基于路况识别的车辆零部件应变载荷预测方法 | |
CN109583741B (zh) | 行驶轨迹风险评估方法和系统 | |
CN114860715A (zh) | 一种流量实时预测的Lanczos时空网络方法 | |
Zheng et al. | A novel prediction model for car body vibration acceleration based on correlation analysis and neural networks | |
CN114692692B (zh) | 一种基于微波衰减信号融合核极限学习机的降雪识别方法 | |
CN116542139A (zh) | 一种液体喷射抛光表面粗糙度的预测方法及装置 | |
Li et al. | A Road Roughness Estimation Method based on PSO-LSTM Neural Network | |
Li et al. | Reliability optimization design method based on multi-level surrogate model | |
CN113705891A (zh) | 基于mra-bas-bp算法的城市商办综合体建筑停车需求预测方法 | |
Grashorn et al. | Transport map Bayesian parameter estimation for dynamical systems | |
WO2019232737A1 (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 |