基于改进近似贝叶斯计算的损伤识别方法技术

本发明专利技术涉及一种基于改进近似贝叶斯计算的损伤识别方法。首先，利用近似贝叶斯计算使得参数后验概率分布的求解无需计算参数的似然函数，解决了实际工程应用中似然函数往往无法求解的难题；蒙特卡罗马尔科夫链抽样先建立与求解问题具有相似性的概率模型，再对该模型进行随机抽样，然后利用抽取的样本求出其统计特征估计值，并作为原问题的近似解。但对复杂工程问题来说，蒙特卡罗马尔科夫链抽样方法所需的计算量往往非常大，因此本发明专利技术在抽样过程中利用随机响应面快速计算参数样本所对应的结构响应统计特征值，避免了调用有限元模型进行数值求解，因而大幅提高了计算效率，解决了贝叶斯方法在多参数、大样本量下由于计算量过大而无法实现的问题。

【技术实现步骤摘要】

本专利技术涉及一种。
技术介绍
传统实际工程结构长期处于复杂的运营环境和外荷载作用中，不可避免地会发生 不同程度和类型的损伤。随着时间的增长，已有损伤会不断累积，导致结构性能不断退化。 若不能及时发现损伤并采取有效的加固措施，在极端条件下就可能导致结构发生灾难性事 故。损伤识别作为结构健康监测系统的核心内容 [1]，是近些年来相关领域的研究热点之一。 当前已有的损伤识别方法可以归为确定性和不确定性方法两大类[2_ 4]。确定性方法中参数 和响应都被看作是确定值，所构建的损伤识别指标也是确定性的，仅能有效地应用在实验 室试验模型或精确控制的特定结构上，目前还未有能够在实际工程中通用的可靠方法 [5]。 因此在现实中，理论上可行的方法在实际应用时的效果往往很不理想，最终导致工程人员 不得不重新依赖费时费力的常规检测技术。因而考虑到工程结构的构件本身、所处的工作 环境以及所受到的外荷载都不可避免地存在不确定性因素(如构件几何尺寸误差、材料离 散性、边界和连接条件的不确定、荷载随机性、环境噪声等)，在损伤识别过程中结合概率统 计分析方法是目前的发展趋势，其也可以看作是对确定性方法的拓展。 目前已有的概率损伤识别方法主要采用了经典概率统计分析理论[6]、贝叶斯理论 [7_ 12]或随机有限元方法[13_16]。经典概率统计方法基于现有的样本观测值，通过构建适当的 估计量和假设检验方法来求取未知参数的统计值。其在选取检验统计量时往往很困难，且 无法利用参数的先验知识，也不考虑后续样本所能提供的信息，因此在具体应用上带有一 定的局限性。随机有限元法主要通过对实测数据或模型参数进行摄动式的随机模拟来获取 参数的概率统计特征，其中模型误差和测量噪声对参数的影响可以通过对实测数据的统计 平均来降低。然而当参数摄动范围较大时，摄动法的精度会明显降低。同时对常用的一阶灵 敏度摄动而言，其分析结果往往是局部收敛的，并且参数初值的选取对结果有着很大影响， 因此在应用上具有较大的局限性。贝叶斯方法结合了参数的先验信息(具体表现为先验概 率分布)及当前的实测数据，再基于最优概率模型来确定参数的后验概率分布。其主要优点 是可以充分利用历史数据或专家经验（即先验信息)，并结合当前实测数据来不断更新参数 的概率分布，这一点非常符合在线健康监测和损伤识别的需求。但参数后验概率分布公式 中的正则化常数往往无法求解，需要采用马尔可夫链蒙特卡罗方法 [17]来求得后验分布的 近似解。当结构模型比较复杂且包含较多未知参数时，贝叶斯修正问题的求解也会变得十 分困难且计算量大幅增加。更常见的问题是无法获取似然函数的表达式，导致贝叶斯方法 的实用性大打折扣。为此，对传统贝叶斯方法进行改进，以求提出适用于工程实际问题的贝 叶斯损伤识别方法，是本专利技术的主要目的。 总体而言，贝叶斯方法由于能充分利用已有数据和新的测试数据，以此不断修正 参数的后验概率分布估计，这点对在线损伤识别是非常有利的。因此，若能对传统贝叶斯方 法进行改进，解决似然函数求解困难的问题，并提高样本响应的计算效率，对工程实际应用 来说具有重要的理论和实用意义。
技术实现思路
本专利技术的目的在于提供一种。 为实现上述目的，本专利技术的技术方案是：一种基于改进近似贝叶斯计算的损伤识 别方法，首先，根据专家经验或历史数据假设结构随机参数的先验概率分布；其次，基于概 率配点法和回归分析建立联系结构随机参数和响应的随机响应面；再次，基于参数先验概 率分布随机生成初始值，通过转换函数进行采样，并利用随机响应面快速计算结构响应的 统计特征值，根据目标函数及设置的阈值判断是否接受样本；然后，进一步计算样本的接受 概率，以确定是否最终接受该样本；随后不断重复前两步，实现循环抽样，直到形成一条稳 定的马尔科夫链，并根据包含在该链里的所有参数样本来计算参数的后验概率分布；最后， 根据估计的参数后验概率分布来构建损伤指标，判断结构是否发生损伤。 在本专利技术实施例中，该方法的具体实现步骤如下： 步骤Sl :根据历史数据或专家经验假设结构随机参数 @噙=H..的先验概率分布贫以建立参数的初始贝叶斯模型；其间若无法确 定先验分布的类型，可以先假设为均勻分布露(4)=U(O5I); 步骤S2 :将随机参数@用标准随机变量.表示，并基于概率配点法和回归分析建立随 机响应面，表现为联系随机参数和响应的显式表达式：本文档来自技高网...

【技术保护点】
一种基于改进近似贝叶斯计算的损伤识别方法，其特征在于：首先，根据专家经验或历史数据假设结构随机参数的先验概率分布；其次，基于概率配点法和回归分析建立联系结构随机参数和响应的随机响应面；再次，基于参数先验概率分布随机生成初始值，通过转换函数进行采样，并利用随机响应面快速计算结构响应的统计特征值，根据目标函数及设置的阈值判断是否接受样本；然后，进一步计算样本的接受概率，以确定是否最终接受该样本；随后不断重复前两步，实现循环抽样，直到形成一条稳定的马尔科夫链，并根据包含在该链里的所有参数样本来计算参数的后验概率分布；最后，根据估计的参数后验概率分布来构建损伤指标，判断结构是否发生损伤。

【技术特征摘要】
1. 一种基于改进近似贝叶斯计算的损伤识别方法，其特征在于：首先，根据专家经验 或历史数据假设结构随机参数的先验概率分布；其次，基于概率配点法和回归分析建立联 系结构随机参数和响应的随机响应面；再次，基于参数先验概率分布随机生成初始值，通过 转换函数进行采样，并利用随机响应面快速计算结构响应的统计特征值，根据目标函数及 设置的阈值判断是否接受样本；然后，进一步计算样本的接受概率，以确定是否最终接受该 样本；随后不断重复前两步，实现循环抽样，直到形成一条稳定的马尔科夫链，并根据包含 在该链里的所有参数样本来计算参数的后验概率分布；最后，根据估计的参数后验概率分 布来构建损伤指标，判断结构是否发生损伤。2. 根据权利要求1所述的基于改进近似贝叶斯计算的损伤识别方法，其特征在于：该 方法具体实现步骤如下： 步骤Sl :根据历史数据或专家经验假设结构随机参数 @噙=H..#)的先验概率分布贫以建立参数的初始贝叶斯模型；其间若无法确 定先验分布的类型，可以先假设为均勻分布露(4)=[/(o5i); 步骤S2 :将随机参数@用标准随机变量.表示，并基于概率配点法和回归分析建立随 机响应面，表现为联...

【专利技术属性】
技术研发人员：方圣恩，董照亮，姜绍飞，林友勤，
申请(专利权)人：福州大学，
类型：发明
国别省市：福建;35