环境振动模态阻尼辨识中时域、频域与子空间估计器的基准比对

《Measurement and Control》：Benchmarking time-domain, frequency-domain, and subspace estimators for ambient-vibration modal damping identification

　　

阻尼比是结构健康监测(SHM)与运行模态分析(OMA)中最难可靠辨识且离散度最大的模态参数；可信基线依赖于多估计器交叉验证与不确定性量化。针对高层钢筋混凝土(RC)框架-剪力墙建筑在低幅值环境振动下的实测阻尼数据稀缺问题，研究人员以一栋26层RC框架-剪力墙住宅建筑为对象，基于八层实测三分量环境振动记录与顶层连续参考点，构建了仅依赖开源Python库的运行模态分析流程。研究人员应用三种相互独立的估计器——随机递减技术(RDT)、协方差驱动随机子空间辨识(SSI-COV)与半功率带宽法(HPBW)——进行基准比对，估算水平与竖向模态阻尼比并量化其不确定性。测得环境模态阻尼比水平向约0.26%–0.6%，竖向约1.0%–1.6%，竖向阻尼显著高于水平向。对频谱孤立模态，RDT与SSI-COV一致性高，单模差异≤0.2%；对基本水平模态（其结构峰值邻接另一谱峰），采用参数化子空间估计，而HPBW系统性偏高仅作粗参考。由此为该类中高层建筑提供了带量化不确定性的可信环境阻尼基线及基于开源库（无商业软件）的辨识工作流，支撑后续SHM基线建立与长期状态跟踪。

研究背景方面，全球城市化推进使高层RC框架-剪力墙结构成为大中城市居住与公共建筑主流形式，设计使用年限50–70年，服役期内经历多次中小地震、长期风载、混凝土收缩徐变与材料老化，动力特性缓慢演化；SHM通过周期模态参数辨识早期探测性能退化已成全生命周期安全工程必需。OMA基于环境振动无需人工激励、不干扰正常使用且可重复，是土木结构动力特性辨识基本工具。三大模态参数中，频率与振型由质量刚度分布主导可从峰位与传递关系稳定获取，阻尼比反映材料内摩擦、连接界面摩擦、非结构构件与土-基础相互作用综合效应，对辨识方法、激励水平与噪声高度敏感，估计离散度最大。已有研究多源于大风强震较大幅值响应且对象多为钢或混合超高层，高层RC框架-剪力墙建筑低幅环境振动实测模态阻尼数据仍相对稀缺。阻尼比常在刚度可观测变化前因微裂缝扩展、连接松动、非结构损伤与材料老化而改变耗能机制，是潜在早损敏感指标，但其用于退化预警须以健康态实测基线与普通波动范围为前提。因此研究人员开展本研究，以填补该类建筑低幅环境阻尼实测空白并提供多估计器交叉验证的不确定性量化基线。

研究人员以一栋新建26层RC框架-剪力墙住宅（总高78 m，平面40 m×22 m，桩基，C30–C50混凝土，设计使用年限70年，7度0.10g设防，II类场地）为对象，采用IGU-BD3C-5三分量宽带智能地震仪在F0、F4、F7、F10、F14、F17、F20、F23八层布点及F26顶层连续参考，每点约40 min记录，共九测点沿高分布，日常环境振动为激励。原始4000 Hz经8阶Butterworth抗混低通后抽稀至100 Hz，Welch法（60 s Hann窗，50%重叠，Δf≈0.0167 Hz）估功率谱密度，各session均79段，模态信噪比12.6–19.3 dB。

主要关键技术方法包括：信号预处理（抗混滤波、抽稀、Welch功率谱）；模态频带定位（在多楼层多session谱上锚定突出峰，基本水平模态搜索带收窄以避开邻接非结构共模线）；三类阻尼估计器并行——频域HPBW按(f₂?f₁)/(2f₀)取?3 dB半功率点；时域RDT经零相位4阶Butterworth带通（半带宽≈0.25f₀）、1σ正斜率穿越触发、≥40次触发、Hilbert包络对数线性拟合（峰至15%）得ζ_RDT=?s/(2πf₀)；SSI-COV将顶层参考与楼层测点六通道再抽稀至25 Hz，构50块行Toeplitz协方差阵，SVD取前n阶可观矩阵，模型阶20–80步长2，稳定图判据Δf<1%、Δζ<5%，连续时间转换ζ_i=?Re(λ_i)/|λ_i|；最后以八session为独立单元聚合并以session间标准差为统一不确定性度量，全流程仅用numpy、scipy、obspy开源库。

研究结果部分，模态频带定位小节指出，各方向分量工程频带内突出谱峰作为阻尼估算锚点，水平E–W与N–S合并报，峰频仅作内部锚点不报结论。模态身份校验小节通过顶层连续参考与逐层同步记录，用跨谱频域分解(FDD)非参振型与SSI-COV参振型比对（MAC基本水平1.00、二水平0.86），各层与顶层幅值平方相干近1确认全局性；排除一条三分量同频、竖向仅低18 dB、无沿高一致振型与低相干的固定频率非结构共模线（判为设备源）；二竖向模态沿高相干约0.46偏低，作为频带等效阻尼而非纯轴模报告。RDT阻尼估计小节得各模session均值±std：一水平0.68±0.16%、二水平0.59±0.19%、一竖向1.19±0.23%、二竖向1.58±0.13%，竖向明显高于水平，一水平因邻峰窄带未能全分离致偏高散大。SSI-COV阻尼估计小节聚八session全阶稳定极中位值得：一水平0.26%（session n=7，右偏std 0.76%但图内MAD≈0.06%）、二水平0.54±0.12%、一竖向1.06±0.25%、二竖向1.39±0.47%；孤立三模与RDT差≤0.2%，一水平因参数化分离结构极显著低于RDT。HPBW交叉检验小节系统性偏高（一水平1.58±0.07%、二水平1.38±0.89%、一竖向1.35±0.59%、二竖向2.95±0.76%），偏差源于有限频率分辨率离散拓宽(≈Δf/(2f₀))与邻峰/噪声底非理想峰形两类，仅作粗参考。实测阻尼基线小节综合取RDT与SSI-COV为可靠基线、一水平采SSI-COV中位0.26%为值、RDT 0.68%为上界，整体水平0.26%–0.6%、竖向1.0%–1.6%，沿高各模散度无单调趋势。参数鲁棒性小节以一竖向为例扫RDT触发0.5σ–2σ（1.23%–1.28%）、签名10–25周期（0.94%–1.56%）、带宽0.20–0.35f₀（1.15%–1.25%）及SSI块行40–70（1.43%–1.59%），证多数控参不敏，签名长度为中等敏感固有权衡，全四模扫掠结论一致。

讨论与结论翻译部分，研究人员指出本工作贡献不在新算法而在为稀缺数据类建筑提供带量化不确定性的实测交叉验证环境阻尼基线：三估计器框架中RDT与SSI-COV原理互异结果互证，分歧作诊断（如一水平邻接非结构线暴露），HPBW有限分辨率下系统偏高被定量刻画；给出该类高层RC框架-剪力墙低幅环境振动水平0.26%–0.6%、竖向1.0%–1.6%可信基线及散度；全流程开源库实现无商业软件。测得竖向阻尼显著高于水平向，孤立模态RDT与SSI-COV差≤0.2%，一水平因邻峰采SSI-COV 0.26%为基、RDT 0.68%为上界；HPBW仅粗参考。该基线支撑后续SHM阻尼状态变量参考与长期条件跟踪，更紧约束一水平需更长更净记录，二竖向轴模/楼板参与分离需更密竖向阵列，此为局限。论文发表于《Measurement and Control》。