3. 振动建模进阶：多自由度系统建模，模态分析基础，模态叠加原理

好，咱们接着聊振动建模。

上一节我们讲了单自由度系统，说白了就是一个弹簧上挂一个质量块。但实际伺服系统哪有这么简单？你想想看，一个工业机器人，光关节就有六七个，每个关节都有柔性，再加上手臂本身的弹性变形——这妥妥的是一个多自由度系统。

我个人习惯，遇到复杂问题先别慌。多自由度系统听着吓人，但核心思想跟单自由度是一样的。只不过从「一个方程」变成了「一组方程」。

3.1 多自由度系统建模

先看一个最简单的例子：两质量块系统。

两个质量块 m₁ 和 m₂，中间用弹簧 k₂ 连着，m₁ 左边还有个弹簧 k₁ 固定到墙上。每个质量块都有自己的位移 x₁ 和 x₂。

写运动方程，其实就是牛顿第二定律的矩阵形式：

M·ẍ + C·ẋ + K·x = F

其中：

• M 是质量矩阵，对角阵，m₁ 和 m₂ 在主对角线上
• K 是刚度矩阵，注意它不是对角阵——因为两个质量块之间有耦合
• C 是阻尼矩阵，实际工程中常用比例阻尼假设
• F 是外力向量

展开写就是：

m₁·ẍ₁ + (k₁+k₂)·x₁ - k₂·x₂ = F₁
m₂·ẍ₂ - k₂·x₁ + k₂·x₂ = F₂

看到了吗？x₁ 的方程里出现了 x₂，x₂ 的方程里出现了 x₁。这就是耦合。我在项目中遇到过，有些工程师一看到耦合就头大，想用解耦的方法硬拆。其实没必要——模态分析就是干这个的。

3.2 模态分析基础

为什么要做模态分析？说白了，就是找到系统「最自然的振动方式」。

你想想看，一个系统有 n 个自由度，理论上就有 n 个固有频率和 n 个振型。每个振型描述了系统在对应频率下「怎么振」——哪些地方振幅大，哪些地方振幅小，相位关系如何。

求解过程其实不复杂：忽略阻尼和外力，解特征值问题。

(K - ω²·M)·φ = 0

这里 ω² 是特征值，对应固有频率的平方。φ 是特征向量，对应振型。

我建议你记住一个重要的性质：振型关于质量矩阵和刚度矩阵是正交的。什么意思？

φᵢᵀ·M·φⱼ = 0  (i ≠ j)
φᵢᵀ·K·φⱼ = 0  (i ≠ j)

这个正交性，就是后面模态叠加法的数学基础。

3.3 模态叠加原理

好，现在到了最精彩的部分。

既然振型是正交的，那我们可以把系统的实际振动分解成各个模态的线性组合。就像傅里叶变换把时域信号分解成不同频率的正弦波一样——模态叠加把物理空间的振动分解成模态空间的独立振动。

数学表达：

x(t) = φ₁·q₁(t) + φ₂·q₂(t) + ... + φₙ·qₙ(t)

其中 qᵢ(t) 是模态坐标，每个 qᵢ 对应一个单自由度系统。

代入原方程，利用正交性，耦合的方程组就变成了 n 个独立的单自由度方程：

m̃ᵢ·q̈ᵢ + c̃ᵢ·q̇ᵢ + k̃ᵢ·qᵢ = f̃ᵢ

m̃ᵢ、c̃ᵢ、k̃ᵢ 分别是模态质量、模态阻尼、模态刚度。f̃ᵢ 是模态力。

为什么会这样？因为正交性把耦合项全部消掉了。每个模态独立振动，互不干扰。

3.4 模态截断与工程应用

实际工程中，我们不可能考虑所有模态。一般只保留前几阶低频模态，高频模态直接忽略。这叫「模态截断」。

我个人的经验法则：

• 伺服系统的带宽通常在 10-100Hz 范围内
• 保留带宽 2-3 倍以内的模态就够了
• 更高频率的模态，要么被滤波器衰减，要么对系统响应贡献极小

3.5 知识体系总览

下面这张图，是我自己总结的多自由度系统建模与模态分析的知识框架。你可以把它当作本章的「地图」：

嗯，这张图把整个流程串起来了。从物理建模开始，到数学方程，再到模态分析，最后落到工程应用。你跟着这个脉络走，就不会迷路。

3.6 一个小例子

最后给个简单的数值例子，帮你巩固一下。

假设两质量块系统：m₁=1, m₂=2, k₁=100, k₂=200。忽略阻尼。

质量矩阵和刚度矩阵：

M = [1  0]    K = [300  -200]
    [0  2]        [-200  200]

解特征值问题，得到：

ω₁² = 50.0  →  f₁ = 1.13 Hz
ω₂² = 350.0 →  f₂ = 2.98 Hz

振型：
φ₁ = [1.0, 0.5]ᵀ   → 两个质量同向运动
φ₂ = [1.0, -1.0]ᵀ  → 两个质量反向运动

你看，第一阶模态两个质量一起晃，第二阶模态它们对着晃。物理意义非常直观。

好了，多自由度建模和模态分析的基础就这些。记住三个关键词：耦合、正交、解耦。下一节我们会把这些理论用到实际的伺服减振算法中。

---

公众号:蓝海资料掘金营，微信deep3321