基于EKF与UKF的土木工程结构参数辨识与损伤诊断系统

资 源 简 介

该项目专门针对土木工程结构在环境激励或强震作用下的动力响应特征，开发了一套基于非线性卡尔曼滤波理论的结构健康监测与损伤识别方案。其核心功能是通过实时处理传感器采集的加速度、速度或位移时程信号，在线估计结构系统的关键物理参数，如非线性刚度、阻尼系数以及本构模型参数。系统首先构建非线性状态空间模型，采用参数扩展技术将待识别的物理参数与位移、速度等状态量整合，构建增广状态向量。通过扩展卡尔曼滤波（EKF）或无迹卡尔曼滤波（UKF）算法，在存在测量噪声和过程噪声的环境下，利用预测-校正机制实现对结构损伤演化过程的实时跟踪。该项目能够有效捕获结构在极端载荷下的刚度退化轨迹，并能够对多自由度复杂结构的损伤位置和损伤程度进行定量评估。系统的应用领域涵盖超高层建筑的抗震监测、大跨度桥梁的健康诊断以及工程结构的精细化动力特性辨识，通过该算法可以显著提高结构损伤识别在噪声干扰下的精度与稳定性，为结构的维修加固和安全性评估提供核心算法支撑。

详 情 说 明

基于扩展与无迹卡尔曼滤波的土木工程结构非线性参数辨识与损伤诊断系统

项目介绍

本系统是一个基于非线性滤波理论的结构健康监测（SHM）仿真平台。它针对土木工程中典型的多自由度（MDOF）结构，通过实时处理含有噪声的动力响应数据，在线辨识结构的刚度和阻尼参数。系统能够敏锐地捕获由于极端荷载或长期老化引起的结构损伤（如刚度突降），并对比了扩展卡尔曼滤波（EKF）与无迹卡尔曼滤波（UKF）在计算精度、收敛速度和抗噪性方面的表现。

功能特性

• 双算法对比辨识：集成EKF与UKF两种主流非线性滤波算法，实现了对结构物理参数的高精度在线追踪。
• 损伤演化模拟：具备模拟结构刚度时变特性的能力，能够复现结构在运行过程中发生的突发性损伤过程。
• 多源观测融合：算法同时利用位移响应和加速度响应作为观测输入，增强了系统状态的可观测性和稳定性。
• 动态仿真引擎：内置基于Newmark-beta法的动力学仿真模块，用于生成逼真的结构动力响应基准数据。
• 定量损伤评估：系统根据辨识出的刚度退化轨迹，自动计算各个楼层的损伤程度百分比并输出分析报告。
• 可视化残差分析：提供直观的辨识结果对比图和状态估计误差分析，便于评估算法性能。

系统要求

• 软件环境：MATLAB R2016b 或更高版本。
• 核心需求：需要支持矩阵运算的基础环境，无需额外工具箱。

实现逻辑说明

系统的执行逻辑严格遵循以下流程：

1. 环境配置与结构定义：初始化采样频率、仿真时长，并定义一个二自由度（2-DOF）结构的质量、初始刚度和初始阻尼常数。
2. 损伤场景设置：人为设定刚度退化逻辑，例如在特定时间点减小层刚度，模拟结构受损。
3. 激励与动力仿真：施加正弦与随机噪声叠加的地面驱动力，利用Newmark-beta法求解二阶动力学方程，获得结构的真实位移、速度和加速度。
4. 观测模拟：对真实的响应信号添加高稳态高斯白噪声，模拟传感器在实际工程环境中的采集误差。
5. 增广状态构建：将结构位移、速度与待辨识的刚度、阻尼参数通过参数扩展技术整合为一个8维增广状态向量。
6. 迭代滤波优化：

* EKF逻辑：通过对状态转移方程和观测方程进行一阶泰勒展开，利用雅可比矩阵进行线性化近似，完成预测与校正循环。 * UKF逻辑：采用无迹变换（Unscented Transform），通过Sigma点采样捕捉非线性分布的均值和方差，避免了复杂的导数计算。

1. 结果产出：实时记录并绘制刚度与阻尼的辨识对比曲线，计算最终的损伤定量指标。

核心函数与算法细节分析

1. 动力学仿真模块 采用Newmark-beta数值积分算法。该算法通过控制参数（beta=0.25, gamma=0.5）确保计算的无条件稳定性，将连续的时间二阶微分方程离散化，逐时步求解多自由度结构在时变刚度条件下的精确动力响应。

2. 状态转移模型 构建了包含结构物理参数的非线性状态方程。状态更新考虑了层间相互作用力，将刚度和阻尼视为随机游走过程（Process Noise），使得滤波器能够灵活调整对参数变化的敏感度。

3. 扩展卡尔曼滤波 (EKF) 实现 核心在于雅可比矩阵（Jacobian）的实时求取。系统推导了状态转移矩阵Phi和观测矩阵H对位移、速度以及四个物理参数（k1, k2, c1, c2）的偏导数。通过一阶线性化，在每一个采样周期内进行时间预测和量测更新。

4. 无迹卡尔曼滤波 (UKF) 实现 基于对称采样策略生成2n+1个Sigma点。通过这些采样点在非线性映射后的加权平均来估计后验状态和协方差。相比EKF，该实现能更好地处理强非线性问题，且由于不需要手动计算复杂的雅可比矩阵，具有更高的鲁棒性。

5. 观测模型设计 观测模型将预测的位移状态和根据牛顿第二定律反算的加速度作为预测观测值。这种设计充分利用了加速度计在工程实践中精度高、安装简便的优势，显著提升了参数辨识的精度，尤其是在刚度发生突变的关键时刻。

6. 后处理与损伤定量评估 通过对比初始猜测值、环境激励下的稳定辨识值以及真实物理值，计算刚度损失比例。系统不仅展示了参数的收敛轨迹，还通过位移残差分析直观反映了模型估计的准确性。