阻尼自由振动(阻尼固有振动数)模拟器 返回
振动工程

阻尼自由振动(阻尼固有振动数)模拟器

由弹簧、质量和阻尼器组成的一自由度系统被拨动后释放时"阻尼自由振动"的体验工具。改变质量、弹簧常数、阻尼系数,可实时了解非阻尼固有振动数、阻尼比 ζ、阻尼固有振动数、对数阻尼率 δ,观看振幅指数衰减的振动波形。

参数设置
质量 m
kg
振动物体的等效质量
弹簧常数 k
N/m
系统的恢复刚度。越大振动越快
阻尼系数 c
N·s/m
阻尼器等装置吸收振动能量的强度
评估周期数 n
评估多少个周期后的振幅比
计算结果
非阻尼固有振动数 (Hz)
阻尼比 ζ
阻尼固有振动数 (Hz)
对数阻尼率 δ
振幅半减周期数 (次)
n 个周期后的振幅比 (%)
弹簧-质量-阻尼器系统与阻尼振动 — 动画

轻敲质量后释放,它会以阻尼固有振动数振动,同时振幅按指数衰减。显示衰减包络线(虚线)、连续峰值和对数阻尼率 δ。

阻尼自由振动波形 — 位移 vs 时间
阻尼固有振动数 vs 阻尼比 ζ
理论·主要公式

$$\omega_n=\sqrt{\frac{k}{m}},\qquad \zeta=\frac{c}{2\sqrt{m\,k}}$$

非阻尼固有圆频率 ωₙ 和阻尼比 ζ。m:质量,k:弹簧常数,c:阻尼系数。临界阻尼系数为 c_c = 2√(mk),阻尼比为 ζ = c/c_c。

$$\omega_d=\omega_n\sqrt{1-\zeta^{2}},\qquad \delta=\frac{2\pi\zeta}{\sqrt{1-\zeta^{2}}}$$

阻尼固有圆频率 ω_d(在不足阻尼 ζ < 1 时有效)和对数阻尼率 δ。阻尼略微降低振动频率,δ 表示振幅每个周期的衰减速率。

$$x(t)=x_0\,e^{-\zeta\omega_n t}\cos(\omega_d t),\qquad \frac{x_n}{x_0}=e^{-n\delta}$$

阻尼自由振动的位移 x(t) 和 n 个周期后的振幅比。包络线 e^(−ζωₙt) 包围振动并按指数衰减。

阻尼固有振动数是什么

🙋
用手拉动弹簧上的重物然后释放,它会摇晃一阵子然后停下来。那个"逐渐停下"就是阻尼吧?
🎓
完全正确。在理想的无摩擦世界里,一旦拨动重物,它会以自身的固有振动数永远摇晃。但现实中总是存在"阻尼"。空气阻力、接触面的摩擦、材料内部能量损失,或者故意装上的阻尼器,这些都会逐渐吸走振动的能量。所以现实的自由振动呈现"振动但越来越弱,最后消失"的阻尼振动形态。
🙋
那"阻尼固有振动数"是和固有振动数完全不同的东西吗?
🎓
好问题。阻尼实际上同时做两件事,区分这两件事很重要。首先,阻尼会"略微减缓"振动。无阻尼的理想固有振动数是 fₙ,但有阻尼的系统自由振动时是稍低的"阻尼固有振动数" f_d。公式是 f_d = fₙ·√(1−ζ²),其中 ζ 是阻尼比。阻尼越大振动越慢,当 ζ = 1(临界阻尼)时 f_d 变成零,系统停止振动。但对大多数实际结构的轻阻尼,这个差异微乎其微。阻尼比 5% 时,振动数只下降 0.1% 左右,实工程常直接用简单的非阻尼公式 fₙ。
🙋
差别这么小的话,阻尼固有振动数就不用管了?
🎓
轻阻尼的话可以。但阻尼增大时差别就不能忽视了。特别是到了临界阻尼这个特殊点——阻尼比正好等于 1——阻尼固有振动数会掉到零,系统再也不振动,只是平滑地回到静止位置。左边的阻尼系数拉到最大,你就会看到阻尼比超过 1,阻尼固有振动数变成 0 Hz,波形从"振动"切换到"不振动"的景象。
🙋
阻尼的"另一个功能"是什么?你说有两个。
🎓
第二个是更能看见的。阻尼"减少振幅",衡量这个速度的就是对数阻尼率 δ。连续两个峰值振幅比的自然对数,定义为 δ = ln(xₙ/xₙ₊₁)。每摇一次周期,振幅就变成原来的 e^(−δ) 倍。默认参数(m=10 kg、k=10000 N/m、c=40 N·s/m)下 δ≈0.398,10 个周期后振幅只剩最初的约 1.9%。上面的图中,指数衰减的包络线就是由 δ 决定的。
🙋
对数阻尼率是只能算出来的东西,还是在实际中有用处?
🎓
这就是它大显身手的地方——"反推"。敲结构一下,记录自由振动消退的波形。从峰值振幅衰减的速度可以测出 δ,然后反推阻尼比 ζ。这叫"对数阻尼率法",是振动试验中最简便实用的阻尼测量手段之一。敲一下、录波形、数几个峰值,就能把"看不见的"阻尼这种物性拎出来。

常见问题

非阻尼固有振动数 fₙ 是完全无阻尼理想系统自由振动时的频率,由 fₙ = (1/2π)·√(k/m) 决定。阻尼固有振动数 f_d 是实际阻尼系统自由振动时的频率,由 f_d = fₙ·√(1−ζ²) 表示,略低一些。ζ 是阻尼比。阻尼越大,振动越慢,当 ζ = 1(临界阻尼)时,f_d 变为零,系统不再振动。对于轻阻尼(ζ 为几个百分点),f_d 与 fₙ 的差别极小,工程中常用非阻尼公式就足够了。
对数阻尼率 δ 是阻尼自由振动连续两个峰值振幅比的自然对数,定义为 δ = ln(xₙ/xₙ₊₁) = 2πζ/√(1−ζ²)。它表示每个周期振幅衰减的快度,是个无量纲量。δ 越大,振动衰减越快。实际应用中最有用的是"反演":轻敲结构一下,记录自由振动波形,测量峰值振幅的衰减速率,就能得到 δ,进而求出阻尼比 ζ。这种方法称为对数阻尼率法,是最基本的实验阻尼测量方法之一。
当阻尼比 ζ 达到 1 时为临界阻尼,系统不再振动,只是以最快速度回到静止位置。此时阻尼固有振动数 f_d 为零。ζ 超过 1 时为过阻尼,系统同样不振动但回位更慢。临界阻尼系数为 c_c = 2·√(m·k),阻尼比为 ζ = c/c_c。门自动关闭装置、汽车悬挂、仪表指针等为避免超调,通常设计在临界阻尼附近(ζ ≈ 0.6~0.8)。
大多数实际结构的阻尼比 ζ = 0.01~0.05 很小,此时修正项 √(1−ζ²) 接近 1。例如 ζ = 0.05 时,振动数只降低约 0.1%,小于测量误差。因此计算固有振动数时,常用更简单、无需知道阻尼就能求解的非阻尼公式 fₙ = (1/2π)·√(k/m)。但对于阻尼较大的隔振系统或制振设计,以及 ζ 接近 1 的情况,阻尼固有振动数的降低就不能忽略,需要正确处理 f_d。

现实世界中的应用

实验模态分析和打击试验:研究桥梁、建筑、机械部件的动特性时,用冲击锤敲打结构,记录自由振动,从振幅衰减速率测量对数阻尼率 δ,再反推阻尼比 ζ。这就是"对数阻尼率法",是时间域标准的阻尼测量方法,与频率域的半值宽度法互补。对单一清晰模态的结构,只需一次敲打和几个峰值就能求出阻尼比,是现场最便利的手段。

制振和隔振设计:加入制振材料、动态阻尼器或隔振垫,要验证阻尼是否增加了,用自由振动的衰减速度来确认。制振前后对比打击响应波形,看衰减变快了(对数阻尼率 δ 增大),说明制振有效;反之若波形还在摇晃,说明阻尼器位置或大小需调整。

瞄准临界阻尼的机构设计:门自动关闭装置、汽车减震器、模拟仪表指针、硬盘读头机构等,都要求无超调快速回位,所以设计在临界阻尼附近(阻尼比 ζ ≈ 0.6~0.8)。用这个工具拖动阻尼系数,能直观看到当 ζ 逼近 1 时,阻尼固有振动数降到零,波形从振荡变成平滑衰减的过程。

机械状态监测和故障诊断:回转机械和结构在螺栓松动、出现裂纹、支撑刚度变化时,阻尼特性会改变。定期记录打击响应,追踪对数阻尼率推算的阻尼比,阻尼的剧变往往是结构劣化的早期信号,比固有振动数偏移更敏感,可用于预知维护。

常见误区和注意

最常见的错误是"混淆阻尼固有振动数和非阻尼固有振动数"。两者由 f_d = fₙ·√(1−ζ²) 相关,但轻阻尼(ζ 为几个百分点)时补正项 √(1−ζ²) 接近 1,差别微乎其微,所以工程里常常不分,直接用非阻尼的 fₙ。但阻尼较大的隔振系统或 ζ 接近 1 时,这个差别就不能无视,必须用正确的阻尼固有振动数 f_d。还有个容易混淆的地方是,自由振动观测得到的就是"阻尼固有振动数 f_d"。打击试验从峰值间隔算出的频率是 f_d,不是 fₙ,当阻尼大时需要用 ζ 反推才能得到 fₙ。

其次是"只用相邻一对峰值测对数阻尼率"。δ 的严格定义是 δ = ln(xₙ/xₙ₊₁),即相邻一个周期的振幅比的对数,但实测时噪声干扰会使单对峰值的测量误差很大。更准确的方法是取相隔 n 个周期的两个峰值,用 δ = (1/n)·ln(x₀/xₙ)。阻尼越小,一周期的衰减越细微,越容易埋在噪声里,所以要用间距足够远的峰值对,或对多个峰值做对数轴直线拟合。只用一对峰的阻尼比常常差得离谱。

最后是"以为阻尼一定降低固有振动数"这个思想误区。阻尼固有振动数 f_d 确实低于非阻尼的 fₙ,但那个下降是由 √(1−ζ²) 造成的,轻阻尼时完全可以忽略。更需注意的是,强制振动下的"共振频率"是另一个量,会因阻尼发生不同的偏移。自由振动的阻尼固有振动数、强制振动的共振频率、非阻尼固有振动数这三个是严格不同的,阻尼大时混淆会导致设计失误。这个工具处理的是"自由振动的阻尼固有振动数",和强制激振的共振要区分清楚。

使用指南

  1. 用滑块设置质量(kg)、弹簧常数(N/m)、阻尼系数(N·s/m)。钢制框架结构的典型值是:质量500 kg、弹簧常数50000 N/m、阻尼系数2000 N·s/m。
  2. 非阻尼圆频率ωn和阻尼比ζ自动计算。ζ < 1时为阻尼自由振动,ζ ≥ 1时为过阻尼。
  3. 查看阻尼圆频率ωd、对数阻尼率δ、振幅半减周期数,设定任意周期数n,观察n个周期后的振幅残存率(%)。

具体计算例

建筑制震装置中,质量1000 kg、弹簧常数100000 N/m、阻尼系数5000 N·s/m的一自由度系统。非阻尼圆频率ωn≈1.59 Hz,阻尼比ζ≈0.79,阻尼圆频率ωd≈0.94 Hz,对数阻尼率δ≈0.436。振幅约1.6个周期后半减,10个周期后残存振幅为初值的约1.8%。

工程实际注意事项