👨🎓个人主页
💥💥💞💞欢迎来到本博客❤️❤️💥💥
🏆博主优势：🌞🌞🌞博客内容尽量做到思维缜密，逻辑清晰，为了方便读者。
⛳️座右铭：行百里者，半于九十。
💥1 概述
如果单自由度（SDOF）系统的自由衰减响应（FDR）不能直接获得，则可以使用环境振动数据来估计模态阻尼比。这里使用了随机递减技术（RDT）[1]以及自然激发技术（NExT）[2]。首先，使用[3]在时域中模拟SDOF对白噪声的响应。然后使用 RDT 或 NExT 计算 IRF。最后，将指数衰减拟合到IRF的包络上，得到模态阻尼比。
环境振动估算阻尼比（SDOF）研究
一、研究背景与意义
阻尼比是描述系统能量耗散特性的关键参数，直接影响结构的振动响应、共振振幅及动力稳定性。在工程结构动力学与振动控制领域，准确估算阻尼比对于理解结构物的动力响应、评估其抵抗振动的能力以及预测其长期服役性能至关重要。传统的阻尼比测量方法通常需要施加人为激励，如激振器激振或冲击锤敲击，但在大型结构物或难以接近的结构物上，施加人为激励往往存在困难，且可能对结构物本身造成损伤。因此，利用环境振动信号（如自然风荷载、地面脉动、设备运行扰动等随机激励产生的振动）估算阻尼比，具有无需人工施加激励、不干扰结构正常运行、适用于大型或现役结构等显著优势，已成为结构动态特性识别的重要手段。
二、SDOF系统振动特性与阻尼比物理意义
单自由度（SDOF）系统是指仅需一个独立坐标即可描述运动状态的简化模型，如简支梁的一阶振动、悬臂结构的横向振动等。环境振动作为一种随机激励，其信号具有宽频带、非平稳、低幅值等特点，对应的SDOF系统响应同样表现为随机过程。与脉冲或正弦激励下的响应不同，环境振动响应不含明显的瞬态衰减过程，需通过功率谱分析、自相关函数等统计方法提取振动特性。
阻尼比表征了结构物能量耗散的能力，其物理意义在于：
- 表征结构物抵抗振动的能力：阻尼比越大，结构物的能量耗散能力越强，受到外部激励时产生的振动响应越小。
- 直接影响结构的共振峰值和振动衰减速度：阻尼比的变化可以反映结构物的损伤程度，例如材料疲劳、连接松动等。
三、环境振动信号特性与数据预处理
- 信号采集：环境振动信号通常通过加速度传感器采集，因为加速度响应信噪比更高，便于后续分析。采样频率需满足奈奎斯特准则，通常为系统最高固有频率的3-5倍，对于土木结构，采样频率一般取10-100 Hz。
- 预处理步骤：
- 去趋势：消除信号中的线性或多项式趋势项，如温度变化导致的结构缓慢变形，可通过最小二乘法拟合趋势并减去实现。
- 滤波：采用带通滤波器去除低频干扰（如传感器漂移）和高频噪声（如电子噪声），保留目标频段（围绕系统固有频率）的信号。
- 去噪：对于含突发噪声的信号，可采用小波阈值去噪或奇异值分解（SVD）方法，在保留有效振动信息的同时抑制噪声。
- 平稳性检验：环境振动信号可能因激励源变化（如阵风、车流突变）呈现非平稳特性，而多数阻尼比估算方法基于平稳随机过程假设。常用的平稳性检验方法包括视觉检查（观察信号时域波形，若振幅和频率成分无明显趋势性变化，则近似为平稳）、统计检验（计算不同时间段的均值和方差，若差异在允许范围内，则认为平稳）以及功率谱估计（采用短时傅里叶变换（STFT）分析功率谱随时间的变化，若谱峰位置和幅值稳定，则满足平稳性）。
四、基于环境振动的阻尼比估算方法
- 峰值法：
- 原理：基于结构传递函数的峰值与阻尼比之间的关系。通过传感器采集结构物在环境激励下的振动响应数据，对采集到的时域信号进行快速傅里叶变换（FFT）等频谱分析，得到结构的功率谱密度（PSD）函数或传递函数。在频谱图上识别结构的固有频率，即峰值对应的频率。利用峰值与背景噪声的比值来估算阻尼比。
- 优点：方法简单，计算量小，易于实现。
- 缺点：精度较低，受噪声影响较大。背景噪声的准确估计是关键，在噪声较大的情况下，峰值法可能失效。另外，该方法假设频谱峰值形状接近理想的洛伦兹曲线，在实际应用中，这种假设可能不成立。
- 适用场景：适用于环境激励强度较强、噪声水平较低、结构响应清晰的简单结构。
- 半功率带宽法：
- 原理：基于结构传递函数峰值附近频率的带宽与阻尼比之间的关系。通过传感器采集结构物在环境激励下的振动响应数据，对采集到的时域信号进行快速傅里叶变换（FFT）等频谱分析，得到结构的功率谱密度（PSD）函数或传递函数。在频谱图上识别结构的固有频率，即峰值对应的频率。在峰值两侧找到功率谱密度下降到峰值一半的两个频率点，分别记为f1和f2。计算阻尼比：ζ≈ (f2 - f1) / (2 * fn)，其中fn为结构的固有频率。
- 优点：精度相对较高，比峰值法更可靠。对噪声的敏感性低于峰值法。
- 缺点：对频谱分辨率要求较高，需要准确识别半功率点。
- 适用场景：适用于频谱分辨率较高、峰值形状较好的结构。
- 频域分解法（FDD）：
- 原理：基于奇异值分解（SVD）的模态参数识别方法。其基本思想是将结构在环境激励下的响应信号进行频谱分析，然后对每个频率点的功率谱密度矩阵进行奇异值分解，从而识别结构的固有频率、振型和阻尼比。通过多个传感器采集结构物在环境激励下的振动响应数据，利用采集到的多通道时域信号，计算结构的功率谱密度矩阵。对每个频率点的功率谱密度矩阵进行奇异值分解，得到奇异值和奇异向量。通过观察奇异值谱，识别结构的固有频率，即奇异值谱中的峰值对应的频率。与固有频率对应的奇异向量可以近似作为结构的振型。采用曲线拟合的方法，将奇异值谱在固有频率附近的曲线拟合为单自由度体系的频率响应函数，然后利用半功率带宽法或更复杂的曲线拟合方法来估算阻尼比。
- 优点：可以同时识别多个模态参数，包括固有频率、振型和阻尼比。对噪声具有一定的鲁棒性。
- 缺点：计算量较大，需要多个传感器进行数据采集。阻尼比的估算精度受曲线拟合方法的影响较大。在模态密集的情况下，FDD方法可能难以准确分离各个模态。
- 适用场景：适用于需要识别多个模态参数、数据采集较为方便的结构。
- 随机子空间辨识法（SSI）：
- 原理：基于状态空间模型的模态参数识别方法。该方法首先将结构的动力学方程转化为状态空间模型，然后利用状态空间模型的观测数据来辨识系统的模态参数，包括固有频率、振型和阻尼比。通过传感器采集结构物在环境激励下的振动响应数据，利用采集到的时域信号，构建观测矩阵。利用随机子空间辨识算法（例如Unweighted Principal Component, UPCA; Canonical Variate Analysis, CVA; or Balanced Realization, BR）辨识状态空间模型。从辨识得到的状态空间模型中提取结构的模态参数，包括固有频率、振型和阻尼比。
- 优点：对噪声具有较强的鲁棒性，可以在噪声较大的环境下准确识别模态参数。能够处理非平稳环境激励下的结构响应。
- 缺点：计算量较大，需要较多的数据样本。算法的复杂度较高，需要一定的理论基础。
- 适用场景：适用于噪声较大、环境激励非平稳的复杂结构。
五、研究挑战与未来展望
- 非平稳环境振动下的阻尼比识别：实际环境振动信号往往具有非平稳特性，如何从非平稳信号中准确提取阻尼比是当前研究的一个难点。
- 多传感器数据融合方法：利用多个传感器采集的数据进行融合处理，可以提高阻尼比估算的精度和鲁棒性。
- 极端环境下的阻尼特性变化规律：研究结构在强风、地震等极端环境下的阻尼特性变化规律，对于结构的安全评估和抗震设计具有重要意义。
📚2 运行结果
clf;close all;
figure
subplot(211)
plot(t,y)
xlabel('time (s)')
ylabel(' displ (m)')
subplot(212)
pwelch(y,[],[],[],1/dt)
set(gcf,'color','w')
编辑
clf;close all;
figure
hold on; box on;
plot(newT,IRF,'b',newT,envelop,'k');
xlabel('time (s)')
ylabel('normalized displacement')
xlim([0,Ts])
set(gcf,'color','w')
编辑
clf;close all;
figure
hold on; box on;
plot(newT,IRF,'b',newT,envelop,'k');
xlabel('time (s)')
ylabel('normalized displacement')
set(gcf,'color','w')
编辑
clf;close all;
figure
hold on; box on;
plot(newT,IRF,'b',newT,envelop,'k');
xlabel('time (s)')
ylabel('normalized displacement')
set(gcf,'color','w')
编辑
🎉3 参考文献
部分理论来源于网络，如有侵权请联系删除。
[1] Ibrahim, S. R. (1977). Random decrement technique for modal identification of structures. Journal of Spacecraft and Rockets, 14(11), 696-700.
[2] James III, O. H., & Came, T. G. (1995). The natural excitation technique (next) for modal parameter extraction from operating structures.
🌈4 Matlab代码实现
资料获取，更多粉丝福利，MATLAB|Simulink|Python资源获取【请看主页然后私信】
