基于高斯权值-混合粒子滤波的疲劳裂纹扩展预测方法技术

本发明专利技术公开了一种基于高斯权值‑混合粒子滤波的疲劳裂纹扩展预测方法，属于故障预测与健康管理技术领域。本发明专利技术提出了一种高斯权值‑混合粒子滤波算法，其根据获得的裂纹长度观测值定义一个观测概率密度，以该观测概率密度和先验转移概率密度的线性组合作为粒子滤波算法的重要性密度函数，并且通过与先验估计之间的高斯权值作为从观测概率密度中采样得到的粒子的权值，然后与基于主动导波的裂纹监测方法相结合，以实现疲劳裂纹扩展的准确预测。本发明专利技术所提出的高斯权值‑混合粒子滤波算法从重要性密度函数中采样时，同时考虑了先验信息和观测信息，降低了对准确定义疲劳裂纹扩展状态方程的依赖，在结构的疲劳裂纹扩展预测方面具有广泛的应用前景。

【技术实现步骤摘要】

本专利技术涉及一种基于高斯权值-混合粒子滤波的疲劳裂纹扩展预测方法，属于故障预测与健康管理

技术介绍
随着现代工程系统的复杂程度和综合化程度越来越高，传统的事后维修和定期维护策略已经无法满足其保障维护的要求。近年来，故障预测与健康管理(PrognosticsandHealthManagement，PHM)技术可以通过传感器在线获取系统的实际健康状态，并预测其退化情况，从而制定最佳的运行和维护策略，以实现视情维护。因此PHM技术受到越来越多的关注。工程结构作为系统的核心组成部分，其服役环境复杂并承受各种交变载荷，容易产生疲劳裂纹。疲劳裂纹的存在和扩展将严重削弱结构承载能力，甚至导致灾难性的事故的发生。因此对结构的疲劳裂纹进行在线准确地预测具有重要的理论意义和工程应用价值。然而疲劳裂纹扩展过程是一个包含各种不确定性的过程，比如材料性质的不确定性，服役环境的不确定性，以及载荷的不确定性等。同时基于结构健康监测方法的在线裂纹监测也受到各种不确定性因素的影响。这些不确定性因素严重影响了疲劳裂纹扩展预测的准确性。近年来，基于贝叶斯理论的概率方法通过结合裂纹扩展模型与实际的结构健康监测数据以消除这些不确定性的影响，受到越来越多的关注。在结构健康监测方法中，基于导波的结构健康监测方法具有小损伤敏感，监测范围广等优点，被认为是非常有前景的方法之一。在基于贝叶斯理论的方法中，粒子滤波方法由于不需要满足线性和高斯过程假设，被认为非常适用于解决非线性非高斯的疲劳裂纹扩展预测问题。但将粒子滤波方法应用于疲劳裂纹扩展预测时，常常难以准确地定义先验裂纹扩展状态方程，也就是描述当前结构疲劳裂纹扩展的模型，这会加剧粒子滤波算法的粒子退化和多样性匮乏问题，使得裂纹扩展预测结果存在较大的误差。
技术实现思路
本专利技术针对先验定义的疲劳裂纹扩展状态方程通常与结构实际的裂纹扩展存在较大偏差的问题，提出了一种基于高斯权值-混合粒子滤波的疲劳裂纹扩展预测方法，结合基于主动导波的裂纹监测方法实现结构疲劳裂纹扩展的准确预测。本专利技术为解决其技术问题采用如下技术方案：一种基于高斯权值-混合粒子滤波的疲劳裂纹扩展预测方法，包括如下步骤：(1)在结构处于服役条件下，采用基于主动导波的结构健康监测方法对结构的疲劳裂纹进行在线监测，通过获取的导波信号计算损伤因子，并以损伤因子作为裂纹长度的观测值；同时根据结构的实际服役情况定义结构的疲劳裂纹扩展状态空间模型；(2)在结构处于服役条件下，每当获得一个新的裂纹长度观测值，按照定义的混合系数λ，从均匀分布U(0,1)中随机采样得到随机数，如果该随机数大于λ，从观测概率密度中随机采样得到一个表征裂纹长度的粒子，否则从先验转移概率密度中随机采样得到粒子，重复该过程的得到确定数目粒子组成的粒子集；(3)对于从先验转移概率密度中采样得到的粒子，通过观测似然概率密度加权；而对于从观测概率密度中采样得到的粒子，根据先验估计定义一个高斯权值分布，通过粒子与先验估计之间的高斯权值作为粒子的权值，然后对所有粒子的权值进行归一化，得到粒子的归一化权值；(4)通过粒子集和相应的归一化权值表征裂纹长度的后验概率密度，并计算裂纹长度的后验估计；在获得的后验概率密度基础上，将粒子向未来时刻投影得到当前时刻的疲劳裂纹扩展预测结果；(5)根据每个粒子的归一化权值，使用系统重采样算法进行重采样；(6)重复上述步骤(2)，(3)，(4)，(5)，(6)实现在线的疲劳裂纹扩展预测。步骤(2)所述混合系数λ取值为0.5；此外，观测概率密度定义为p(xk|yk)，式中yk为k时刻获得的裂纹长度观测值，xk为k时刻的裂纹长度；观测概率密度的概率密度形式与表征观测不确定性的概率密度一致，以观测值yk对应的裂纹长度为均值，k时刻先验粒子集的标准差作为观测概率密度的标准差。步骤(3)中定义了一个高斯权值分布，对从观测概率密度中采样得到的粒子加权，权值计算方式如下，式中，是从观测概率密度p(xk|yk)中采样得到的第i个粒子，为第i个粒子在k-1时刻的权值；为第i个粒子在k时刻的权值；σv为表征观测不确定性的概率密度标准差；为定义的高斯权值分布，其均值为裂纹长度的先验估计标准差σp为经验设定值。本专利技术的有益效果如下：本专利技术提出了一种基于高斯权值-混合粒子滤波的疲劳裂纹扩展预测方法，以观测概率密度和先验转移概率密度的线性组合作为粒子滤波的重要性密度函数。从重要性密度函数中采样得到的粒子集同时考虑了观测信息和先验信息，降低了对准确定义疲劳裂纹扩展状态方程的依赖，同时结合了基于主动导波的裂纹监测方法，能有效地用于结构的疲劳裂纹扩展预测。附图说明图1为本专利技术提出的疲劳裂纹扩展预测方法流程图。图2为实施例中的结构尺寸图。图3为实施例中的结构疲劳裂纹扩展图。图4为实施例中结构健康监测方法获得的导波损伤因子图。图5为实施例中裂纹长度后验估计对比图。图6为实施例中预测的失效循环载荷数的相对误差对比图。具体实施方式下面结合附图对本专利技术创造的技术方案进行详细说明。图1所示，为一种基于高斯权值-混合粒子滤波的疲劳裂纹扩展预测方法的流程图，方法步骤如下，(1)在结构处于服役条件下，采用基于主动导波的结构健康监测方法对结构的疲劳裂纹进行在线监测。在结构的关键部位粘贴压电传感器，通过压电传感器激励和采集结构中的导波信号。利用结构健康监测系统实现导波信号的在线采集和调理，对获取的导波信号，计算其损伤因子作为裂纹长度的观测值；(2)根据结构的实际情况定义疲劳裂纹扩展状态空间模型，包括状态方程和观测方程两个部分，其中状态方程由结构的形状及其服役条件确定，观测方程由基于导波的结构健康监测方法确定，如下所示。状态方程：xk＝xk-1+Δxk-1(xk-1)·exp(ω)观测方程：yk＝g(xk)+v式中，k是离散的时间；xk是k时刻的裂纹长度；xk-1是k-1时刻的裂纹长度；Δxk-1(xk-1)是裂纹扩展增量，为xk-1的函数，由疲劳裂纹扩展模型定义，比如Paris模型，NASGRO模型；ω是一个随机变量，服从高斯分布表征疲劳裂纹扩展的不确定性，其均值选择为以保证exp(ω)的期望为1，指数项exp(ω)保证裂纹扩展增量非负；yk为k时刻获得的导波损伤因子；g(·)是观测映射，表征裂纹长度和损伤因子之间的关系；v是一个随机变量，表征观测的不确定性。(3)定义状态空间模型后，通过结构健康监测方法在线获取结构当前时刻k的裂纹长度观测值yk，每当获得一个新的裂纹长度观测值，从如下所示的混合建议分布中采样，式中，k是离散的时间；xk是k时刻的裂纹长度；xk-1是k-1时刻的裂纹长度，yk是k时刻的裂纹长度观测值；i为粒子序号，i＝1,...,Ns，Ns为粒子总数；为重要性密度函数；为先验转移概率密度；λ为混合系数取为λ＝0.5；p(xk|yk)为观测概率密度，其形式与表征观测不确定性的概率密度一致。对于观测概率密度，其均值为观测值yk对应的裂纹长度g-1(yk)；由上一时刻粒子集通过状态方程一步转移得到先验粒子集以的标准差作为观测概率密度的标准差。从混合建议分布中采样得到粒子集的方式如下：(a)令i＝1(b)从均匀分布U(0,1)中随机采样得到u(i)，如果u(i)＞λ，从p(xk|yk)中随机采本文档来自技高网...

【技术保护点】
一种基于高斯权值‑混合粒子滤波的疲劳裂纹扩展预测方法，其特征在于，包括如下步骤：（1）在结构处于服役条件下，采用基于主动导波的结构健康监测方法对结构的疲劳裂纹进行在线监测，通过获取的导波信号计算损伤因子，并以损伤因子作为裂纹长度的观测值；同时根据结构的实际服役情况定义结构的疲劳裂纹扩展状态空间模型；（2）在结构处于服役条件下，每当获得一个新的裂纹长度观测值，按照定义的混合系数，从均匀分布U(0,1)中随机采样得到随机数，如果该随机数大于，从观测概率密度中随机采样得到一个表征裂纹长度的粒子，否则从先验转移概率密度中随机采样得到粒子，重复该过程的得到确定数目粒子组成的粒子集；（3）对于从先验转移概率密度中采样得到的粒子，通过观测似然概率密度加权；而对于从观测概率密度中采样得到的粒子，根据先验估计定义一个高斯权值分布，通过粒子与先验估计之间的高斯权值作为粒子的权值，然后对所有粒子的权值进行归一化，得到粒子的归一化权值；（4）通过粒子集和相应的归一化权值表征裂纹长度的后验概率密度，并计算裂纹长度的后验估计；在获得的后验概率密度基础上，将粒子向未来时刻投影得到当前时刻的疲劳裂纹扩展预测结果；（5）根据每个粒子的归一化权值，使用系统重采样算法进行重采样；（6）重复上述步骤（2），（3），（4），（5），（6）实现在线的疲劳裂纹扩展预测。...

【技术特征摘要】
1.一种基于高斯权值-混合粒子滤波的疲劳裂纹扩展预测方法，其特征在于，包括如下步骤：（1）在结构处于服役条件下，采用基于主动导波的结构健康监测方法对结构的疲劳裂纹进行在线监测，通过获取的导波信号计算损伤因子，并以损伤因子作为裂纹长度的观测值；同时根据结构的实际服役情况定义结构的疲劳裂纹扩展状态空间模型；（2）在结构处于服役条件下，每当获得一个新的裂纹长度观测值，按照定义的混合系数，从均匀分布U(0,1)中随机采样得到随机数，如果该随机数大于，从观测概率密度中随机采样得到一个表征裂纹长度的粒子，否则从先验转移概率密度中随机采样得到粒子，重复该过程的得到确定数目粒子组成的粒子集；（3）对于从先验转移概率密度中采样得到的粒子，通过观测似然概率密度加权；而对于从观测概率密度中采样得到的粒子，根据先验估计定义一个高斯权值分布，通过粒子与先验估计之间的高斯权值作为粒子的权值，然后对所有粒子的权值进行归一化，得到粒子的归一化权值；（4）通过粒子集和相应的归一化权值表征裂纹长度的后验概率密度，并计算裂纹长度的后验估计；在获得的后验概率密度基础...

【专利技术属性】
技术研发人员：袁慎芳，陈健，邱雷，王卉，
申请(专利权)人：南京航空航天大学，
类型：发明
国别省市：江苏;32