2. 单自由度系统：无阻尼自由振动、固有频率、阻尼自由振动、对数衰减率

各位同学，咱们今天聊聊单自由度系统。这是振动分析的入门砖，也是我这些年做项目最常用的基础模型。说白了，很多复杂结构在特定条件下，都可以简化成一个质量块、一根弹簧和一个阻尼器。你想想看，是不是这个理？

2.1 无阻尼自由振动

先看最简单的情况——没有阻尼。一个质量块挂在弹簧上，你把它拉开然后松手，它就会来回晃悠。这就是无阻尼自由振动。

运动方程很简单：

m·ẍ + k·x = 0

其中 m 是质量，k 是弹簧刚度，x 是位移。这个方程的解是：

x(t) = A·sin(ωₙ·t + φ)

A 是振幅，φ 是初相位，ωₙ 就是我们要找的固有角频率。

我个人习惯把 ωₙ 记作：

ωₙ = √(k/m)

对应的固有频率 fₙ = ωₙ / (2π)，单位是赫兹（Hz）。

2.2 固有频率——系统的"灵魂"

固有频率是什么？说白了，就是系统自己最"喜欢"的振动频率。你给系统一个冲击，它就会以这个频率振动起来。

我在项目中遇到过一件事：某设备在测试时总在某个转速下振动特别大。一开始大家以为是装配问题，后来我算了一下固有频率，发现正好和那个转速对上。嗯，这就是共振。

固有频率只取决于系统的质量和刚度：

• 质量越大，固有频率越低（越"笨重"）
• 刚度越大，固有频率越高（越"硬朗"）

举个例子：

参数	数值	固有频率 fₙ
m = 1 kg, k = 100 N/m	ωₙ = 10 rad/s	fₙ ≈ 1.59 Hz
m = 4 kg, k = 100 N/m	ωₙ = 5 rad/s	fₙ ≈ 0.80 Hz
m = 1 kg, k = 400 N/m	ωₙ = 20 rad/s	fₙ ≈ 3.18 Hz

2.3 阻尼自由振动

现实世界哪有没阻尼的？空气阻力、材料内摩擦、连接处的摩擦……这些都会让振动慢慢停下来。

加上阻尼后的运动方程：

m·ẍ + c·ẋ + k·x = 0

c 是阻尼系数。这个方程的解取决于阻尼比 ζ：

ζ = c / (2·√(m·k))

根据 ζ 的大小，有三种情况：

• ζ < 1（欠阻尼）：振动逐渐衰减，这是最常见的情况
• ζ = 1（临界阻尼）：最快回到平衡位置，没有振荡
• ζ > 1（过阻尼）：慢慢回到平衡位置，也不振荡

欠阻尼情况下，振动频率会变慢：

ω_d = ωₙ · √(1 - ζ²)

ω_d 叫阻尼固有频率。你会发现，阻尼越大，振动频率越低。

2.4 对数衰减率——量化阻尼的工具

怎么测量阻尼比？最实用的方法就是对数衰减率。

你记录下振动波形，取两个相邻的峰值 x₁ 和 x₂，然后算：

δ = ln(x₁ / x₂)

δ 就是对数衰减率。它和阻尼比的关系是：

δ = 2π·ζ / √(1 - ζ²)

如果阻尼很小（ζ < 0.2），可以近似为：

δ ≈ 2π·ζ

反过来，知道了 δ 就能算出 ζ：

ζ = δ / √(4π² + δ²)

举个例子：测得 x₁ = 10 mm，x₂ = 6 mm，那么：

δ = ln(10/6) ≈ 0.511
ζ = 0.511 / √(4π² + 0.511²) ≈ 0.081

阻尼比约 8.1%，属于中等阻尼。

2.5 知识体系总览

下面这张图把本章的核心逻辑串起来了：

这张图把本章的知识点串起来了。从无阻尼到有阻尼，从固有频率到对数衰减率，每一步都是环环相扣的。

2.6 实用建议

最后，给各位几个我在项目中总结的经验：

1. 先算后测：做实验前，先用手算或仿真估算固有频率，这样测试时心里有数
2. 多取几个峰值：算对数衰减率时，别只取两个点，用多个周期拟合更准
3. 注意阻尼的影响：如果阻尼比较大，测到的频率和固有频率差不少，记得修正
4. 别忽视非线性：实际系统多少都有非线性，大振幅和小振幅下的固有频率可能不一样