频率响应函数(FRF)计算 返回
振动分析

频率响应函数(FRF)计算

输入固有频率与阻尼比,实时绘制FRF波德图(幅值与相位)。自动计算共振频率、Q值、半功率带宽及动力放大系数。可选添加第二阶模态进行多自由度分析。

第一阶模态参数
固有频率 fn₁
Hz
阻尼比 ζ₁
%
质量 m₁
kg
添加第二阶模态
计算结果(第一阶)
计算结果
共振频率 fr (Hz)
Q值(品质因数)
半功率带宽 Δf (Hz)
动力放大系数 DAF
最大导纳 (μm/N)
刚度 k (kN/m)
波德图 — 幅值(导纳 |H|)
波德图 — 相位 (°)
理论与主要公式

$H(\omega) = \dfrac{1/k}{1 - r^2 + 2j\zeta r}$


$r = \omega/\omega_n,\quad Q = \dfrac{1}{2\zeta}$

什么是频率响应函数(FRF)

🙋
“频率响应函数”听起来好复杂啊,它到底是什么?
🎓
简单来说,FRF就是描述一个结构“怕什么频率”的函数。想象一下你推一个秋千,如果推的节奏正好对上它晃悠的节奏,它就会荡得特别高。FRF就是告诉你,在不同的推搡频率下,这个秋千(或者任何结构)会晃得多厉害。在实际工程中,我们用它来预测和分析振动问题。
🙋
诶,真的吗?那这个模拟器里那些“固有频率”、“阻尼比”是干嘛用的?
🎓
问得好!“固有频率”就是结构自己最喜欢晃悠的那个频率,比如汽车发动机的某个部件可能有个200Hz的固有频率。“阻尼比”可以理解为晃悠时遇到的阻力大小,就像秋千在空气里晃和在水里晃的区别。你试着拖动上面“阻尼比ζ₁”的滑块,把它从0.01调到0.5,你会看到共振峰从又高又尖变得又矮又平,这就是阻尼在“吃掉”振动能量。
🙋
那“双自由度系统”和“质量比”又是怎么回事?一个秋千还能有两个晃法?
🎓
哈哈,可以这么理解!一个秋千是单自由度,只能前后晃。但现实中的结构复杂多了,比如汽车底盘,它既有上下跳动的“点头”模式,也有前后俯仰的“仰头”模式,这就是两个自由度。在模拟器里,你切换到“双自由度”模式,设置第二个固有频率fn₂和质量比m₂/m₁,你会看到图上出现两个共振峰!这就像给结构做了个“振动体检”,能看出它有几个“怕”的频率点。

物理模型与关键公式

对于最基本的单自由度弹簧-质量-阻尼系统,其位移导纳(柔度)形式的频率响应函数如下,它描述了在频率为ω的正弦力激励下,产生的位移响应幅值与相位。

$$H(\omega) = \frac{X(\omega)}{F(\omega)}= \dfrac{1/k}{1 - r^2 + 2j\zeta r}$$

其中,$k$是刚度(N/m),$r = \omega/\omega_n$是频率比,$\omega_n = 2\pi f_n$是固有角频率,$\zeta$是阻尼比。分母中的虚数单位$j$决定了响应的相位滞后。

共振时的动力放大系数和衡量共振峰尖锐程度的品质因数Q,是工程中评估振动特性的关键指标。

$$DAF \approx \frac{1}{2\zeta} = Q, \quad \Delta f = f_n / Q$$

$DAF$是动力放大系数,表示共振位移是静变形的多少倍。$Q$是品质因数,$Q$值越大,共振峰越尖锐。$\Delta f$是半功率带宽(-3dB带宽),可以直接从实测的FRF曲线中量取,并用于反算阻尼比:$\zeta \approx \Delta f / (2 f_n)$。

现实世界中的应用

汽车NVH(噪声、振动与声振粗糙度)分析:工程师通过计算车身和底盘部件的FRF,预测在发动机怠速或路面激励下,哪些频率可能引起方向盘抖动或座椅发麻,从而在设计阶段就优化结构或增加阻尼材料。

航空航天结构模态试验:在飞机机翼或火箭整流罩上粘贴激振器和传感器,实测其FRF,识别出结构的固有频率和阻尼。这是验证CAE仿真模型、确保飞行中不会因气动载荷引发共振(颤振)的关键步骤。

精密仪器与机床隔振:制造芯片的光刻机或高精度机床对环境振动极其敏感。通过分析其FRF,可以设计出专门的隔振平台,让平台的固有频率远低于机器敏感的频率,从而“过滤”掉地面传来的有害振动。

桥梁与建筑健康监测:通过在大型桥梁上布置传感器阵列,监测其在风或车辆载荷下的振动响应,提取FRF。长期对比FRF的变化(如固有频率下降、阻尼增加),可以早期发现结构刚度退化或损伤,实现预警。

常见误解与注意事项

首先,切勿认为“频率响应函数(FRF)只有一种”。实际上,根据输入与输出的组合,可以定义多种FRF,例如位移/力(柔度)、速度/力(导纳)、加速度/力(惯性)等。仿真中通常显示的是基于位移的“柔度”。例如在振动试验中使用加速度计时,处理的便是惯性FRF,其在低频段噪声显著等特性会有所不同。其次,阻尼比ζ与共振峰值高度并非简单的比例关系。即使ζ从0.1倍增到0.2,峰值的衰减幅度会略大于一半。具体而言,ζ=0.01时的峰值振幅约为50(1/(2ζ)),而ζ=0.1时约为5。最后,设置二自由度系统的“质量比”和“耦合弹簧”时需格外谨慎。例如,若将质量比极端放大(如m2/m1=100),较低的共振模态会近似表现为仅小质量(m1)振动的“局部振动”状态,而较高模态中两个质量会呈反相位振动。在实际工程中设定参数时,若未准确估算预期物理构件的质量与刚度,可能会陷入分析结果与实际完全不符的误区。

使用指南

  1. 在第一模态输入框输入固有频率(Hz)和阻尼比(0-1范围),例如钢梁第一阶共振频率fn1=45Hz、阻尼比ζ=0.08
  2. 输入质量参数m1(kg),系统自动计算FRF幅值和相位角,实时生成波德图
  3. 若分析多自由度系统,勾选"启用第二模态"并输入fn2、ζ2参数,计算耦合响应
  4. 点击"计算"按钮获得各频率点的幅值响应(mm/N)、相位延迟(°)、共振峰值Q值

具体计算示例

某混凝土框架结构第一层,固有频率fn1=32Hz、阻尼比ζ1=0.06、质量m1=50000kg。在共振频率处,FRF幅值为|H(jω)|=0.512mm/N;Q值=1/(2ζ)≈8.33;半功率带宽BW=ζ×fn×2=3.84Hz。若激振力F=5kN,则共振响应幅度约=0.512×5000=2560mm(实际需考虑非线性),相位角在共振点处φ=-90°。

实务注意事项

  1. 阻尼比超过0.707时系统过阻尼,无明显共振峰,适用于隔振器设计验证(例如橡胶隔震支座ζ=0.15)
  2. 多自由度模态间距过小(Δf<5Hz)时易产生驻波干涉,需验证传感器布置和激励方向
  3. 工程测量时,实测FRF应考虑器件幅频特性,加速度传感器通常工作带宽100-5000Hz内
  4. 对于大跨度桥梁(L>500m),环境风振可激发低阶模态,应采用随机减量法(RDT)滤波低频噪声