2. 单自由度系统：弹簧-质量-阻尼模型

各位同学，今天我们来聊聊振动分析里最基础、也最核心的一个模型——单自由度系统。说白了，就是把一个复杂的无人机结构，简化成一个弹簧、一个质量块、一个阻尼器。

你可能会问：无人机那么复杂，一个弹簧能代表啥？

嗯，我刚开始做结构设计时也有这个疑问。后来在一次机臂振动测试中，我发现某个模态的响应，用这个简单模型算出来的结果，跟实测误差不到5%。从那以后，我对这个模型就再也不敢小看了。

2.1 模型建立

先看这个系统的组成：

• 质量块 m：代表无人机某个部件的等效质量
• 弹簧 k：代表结构的刚度
• 阻尼器 c：代表结构的内摩擦、空气阻力等耗能因素

我个人习惯把阻尼器想象成一个「阻力器」——你推它越快，它反抗越厉害。这个比喻在项目中帮了不少忙。

2.2 运动方程建立

根据牛顿第二定律，我们写出运动方程：

m·x''(t) + c·x'(t) + k·x(t) = F(t)

其中：

• x(t) —— 位移
• x'(t) —— 速度
• x''(t) —— 加速度
• F(t) —— 外部激励力

这个方程看着简单，但它的物理意义很深刻。左边三项分别对应：惯性力、阻尼力、弹性力。这三股力量在振动过程中不断博弈，决定了系统的行为。

2.3 固有频率与阻尼比

我们先看无阻尼自由振动的情况（c=0, F=0）。方程简化为：

m·x'' + k·x = 0

解这个方程，得到固有角频率：

ωn = √(k/m)

固有频率（Hz）则为：

fn = ωn / (2π) = (1/2π)·√(k/m)

这里有个经验：我在设计无人机机臂时，发现如果机臂的固有频率跟旋翼转速频率接近，共振会非常明显。有一次试飞，机臂振幅大到肉眼可见，吓得我赶紧改了结构。

接下来引入阻尼比 ζ：

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

阻尼比是个无量纲参数，它决定了系统振动的「性格」。

2.4 三种阻尼状态

根据阻尼比的大小，系统呈现三种完全不同的行为：

状态	条件	行为特征
欠阻尼	ζ < 1	振动衰减，有振荡
临界阻尼	ζ = 1	最快回到平衡，无振荡
过阻尼	ζ > 1	缓慢回到平衡，无振荡

2.4.1 欠阻尼（ζ < 1）

这是最常见的情况。系统受到扰动后，会以固有频率附近振荡，振幅按指数规律衰减。

位移解为：

x(t) = A·e^(-ζ·ωn·t)·sin(ωd·t + φ)

其中 ωd = ωn·√(1-ζ²) 是阻尼固有频率。

我记得有一次做无人机悬停测试，发现机身有轻微的抖动。用加速度计一测，发现阻尼比只有0.02左右。这就是典型的欠阻尼——振动衰减很慢，需要加阻尼材料来改善。

2.4.2 临界阻尼（ζ = 1）

临界阻尼是个分界线。此时系统刚好不振荡，以最快速度回到平衡位置。

位移解为：

x(t) = (A + B·t)·e^(-ωn·t)

临界阻尼在工程中很难精确实现，但它是设计减震器时的理想目标。我曾经设计过一款无人机云台减震器，目标就是让阻尼比尽量接近1。虽然最后只做到0.85左右，但效果已经非常好了。

2.4.3 过阻尼（ζ > 1）

过阻尼系统不会振荡，但回到平衡位置的速度比临界阻尼慢。阻尼越大，恢复越慢。

位移解为：

x(t) = A·e^(-λ1·t) + B·e^(-λ2·t)

其中 λ1, λ2 是两个不同的实数衰减系数。

2.5 工程应用要点

在实际项目中，我总结了几条经验：

1. 先测固有频率：用锤击法或扫频法，确定结构的固有频率
2. 再算阻尼比：用半功率带宽法，从频响曲线上提取阻尼比
3. 最后做调整：根据目标频率和阻尼比，修改结构或添加阻尼材料

举个例子，某次我处理一款六旋翼的机臂振动问题。实测发现机臂一阶弯曲频率是45Hz，而电机转速对应的激振频率是50Hz。两者只差5Hz，共振风险很高。

解决方案很简单：在机臂根部加了一圈橡胶阻尼垫，把阻尼比从0.015提升到0.06。虽然固有频率只降到了43Hz，但共振峰值被大幅抑制了。

这就是单自由度模型的魅力——用最简单的数学，解决最实际的工程问题。

下一章我们会讨论多自由度系统，到时候你会发现，单自由度系统的理解是基础中的基础。打好这个基础，后面的内容就水到渠成了。