二、振动基础理论：单自由度系统、多自由度系统、连续系统的振动方程与求解方法

振动控制工程的核心任务，是建立对结构或机械系统动态行为的数学描述，并据此预测、评估和抑制振动。无论面对的是简单的单自由度弹簧-质量系统，还是复杂的连续弹性体，其分析路径都遵循“建模→列方程→求解→解读”的范式。本章将系统梳理三类典型系统的振动方程建立方法、求解策略及其工程适用边界，为后续的隔振、减振与主动控制技术奠定坚实的理论基础。

2.1 单自由度系统：振动分析的基石

单自由度（SDOF）系统是振动理论中最基本、也是工程应用最广泛的模型。它用一个质量块、一个弹簧和一个阻尼器来表征系统的惯性、恢复力和能量耗散特性。尽管实际结构都是连续体，但在许多工程场景下——例如设备安装在隔振器上、仪表盘固定在支架上——我们关心的是其主导模态，此时将其简化为SDOF模型不仅计算简便，而且物理意义清晰，能够快速给出振动响应的量级估计。

2.1.1 核心原理：运动方程的建立

根据牛顿第二定律或达朗贝尔原理，SDOF系统的运动方程可写为：

m·ẍ(t) + c·ẋ(t) + k·x(t) = F(t)

其中，m为等效质量（kg），c为粘性阻尼系数（N·s/m），k为刚度（N/m），x(t)为位移响应，F(t)为外部激励力。方程左侧三项分别代表惯性力、阻尼力和弹性力，右侧为激振力。该方程是二阶常系数线性常微分方程，其通解由自由振动响应（齐次解）和强迫振动响应（特解）叠加而成。

定义无阻尼固有频率ωn = √(k/m) 和阻尼比ζ = c/(2√(km))，方程可改写为标准形式：

ẍ + 2ζωn·ẋ + ωn²·x = F(t)/m

这两个参数是描述SDOF系统动态特性的核心指标。ωn决定了系统共振频率的位置，ζ则决定了共振峰的尖锐程度和衰减速率。工程中，ζ通常小于0.1（欠阻尼），此时系统在受到冲击后会以ωd = ωn·√(1-ζ²)的频率衰减振荡。

2.1.2 求解方法：从时域到频域

针对不同的激励类型，SDOF系统的求解策略有所不同：

• 自由振动（F(t)=0）：解为 x(t) = e^(-ζωn·t)·[A·cos(ωd·t) + B·sin(ωd·t)]，其中A、B由初始位移和初始速度决定。该解用于评估系统在受到初始扰动后的衰减特性。
• 简谐激励（F(t)=F₀·sin(ωt)）：稳态响应为 x(t) = X·sin(ωt - φ)，其中幅值X = F₀/k / √[(1-r²)² + (2ζr)²]，相位φ = arctan[2ζr/(1-r²)]，r = ω/ωn为频率比。该解是振动控制中最常用的公式，用于计算隔振传递率。
• 任意周期激励：将F(t)展开为傅里叶级数，对每个谐波分量分别求解后叠加。
• 瞬态激励（冲击、阶跃）：使用杜哈梅积分（卷积积分）x(t) = (1/(m·ωd))·∫₀ᵗ F(τ)·e^(-ζωn(t-τ))·sin[ωd(t-τ)]dτ，或采用拉普拉斯变换法。

在工程软件中，数值积分法（如Newmark-β法、龙格-库塔法）被广泛用于处理非线性或时变参数系统，但理解解析解有助于建立直观的物理图像。

2.1.3 案例拆解：精密仪器隔振设计

某半导体光刻机工作台需要隔离来自地面的振动，已知工作台总质量m=800kg，四个支撑隔振器总刚度k=1.2×10⁶ N/m，阻尼比ζ=0.05。地面振动主频为25Hz（由附近空调机组引起），振幅为0.1mm。要求评估工作台的振动传递率。

分析过程：

1. 计算固有频率：fₙ = (1/2π)·√(k/m) = (1/2π)·√(1.2×10⁶/800) ≈ 6.16 Hz
2. 计算频率比：r = f/fₙ = 25/6.16 ≈ 4.06
3. 计算传递率：T = √[(1+(2ζr)²)/((1-r²)²+(2ζr)²)] = √[(1+(2×0.05×4.06)²)/((1-4.06²)²+(2×0.05×4.06)²)] ≈ 0.066
4. 工作台振动幅值 = 0.1mm × 0.066 = 0.0066mm = 6.6μm

该结果满足光刻机对振动幅值小于10μm的要求。若地面振动频率降低至接近6Hz，则传递率将接近1甚至放大，此时需要重新设计隔振器刚度或增加质量。

2.1.4 常见误区与总结

• 误区一：忽略阻尼的影响。在共振区（r≈1），阻尼是限制振幅的唯一因素；在隔振区（r>√2），阻尼反而会降低隔振效果（因为传递率公式分子含2ζr项）。因此隔振器阻尼并非越大越好。
• 误区二：将实际结构直接视为SDOF。只有当系统的高阶模态频率远高于激励频率，且各阶模态耦合较弱时，SDOF近似才有效。否则应使用MDOF模型。
• 误区三：混淆固有频率与共振频率。有阻尼系统的共振频率为ωd = ωn·√(1-2ζ²)，而非ωn。对于小阻尼（ζ<0.1），两者差异小于1%，工程上可忽略。

总结：SDOF系统是振动控制的“乘法口诀”，所有复杂系统的分析最终都归结为对多个SDOF响应的叠加。掌握其频响特性、传递率公式和时域衰减规律，是后续学习的基础。

2.2 多自由度系统：模态叠加与耦合效应

当系统包含多个质量块，或需要同时关注多个位置的振动响应时，必须采用多自由度（MDOF）模型。MDOF系统与SDOF的本质区别在于：各自由度之间存在弹性耦合和阻尼耦合，导致运动方程不再是单个标量方程，而是一组联立微分方程。MDOF系统的振动特性由多个固有频率和对应的振型（模态）共同描述。

2.2.1 核心原理：矩阵形式的运动方程

对于n个自由度的线性系统，运动方程可写为矩阵形式：

[M]·{ẍ} + [C]·{ẋ} + [K]·{x} = {F(t)}

其中[M]、[C]、[K]分别为n×n阶质量矩阵、阻尼矩阵和刚度矩阵，{x}为位移列向量。质量矩阵通常为对角阵（集中质量模型）或一致质量矩阵（有限元模型）；刚度矩阵为对称正定矩阵（无刚体模态时）。

求解MDOF系统的核心思路是模态叠加法：通过求解无阻尼自由振动方程的特征值问题 [K]{φ} = ω²[M]{φ}，得到n个特征值ωᵢ²和对应的特征向量{φᵢ}（振型）。将物理坐标{x}变换为模态坐标{q}：{x} = [Φ]{q}，其中[Φ]为振型矩阵。利用振型关于[M]和[K]的正交性，可将耦合的方程组解耦为n个独立的SDOF方程：

Mᵢ·q̈ᵢ + Cᵢ·q̇ᵢ + Kᵢ·qᵢ = Pᵢ(t)   (i=1,2,...,n)

其中Mᵢ、Cᵢ、Kᵢ为模态质量、模态阻尼和模态刚度，Pᵢ(t)为模态力。每个模态方程均可按SDOF方法求解，最后叠加得到物理响应。

2.2.2 判断方法：模态截断与阻尼处理

工程中，MDOF系统的自由度可能成百上千（有限元模型），但实际参与振动的有效模态通常只有前几阶。模态截断准则：保留所有固有频率低于激励频率上限1.5~2倍的模态，或保留累计模态有效质量达到总质量90%以上的模态。

阻尼矩阵的处理是MDOF分析的难点。若阻尼为比例阻尼（[C] = α[M] + β[K]），则模态方程可解耦；否则需采用复模态分析或直接数值积分。工程中常用模态阻尼比ζᵢ来近似，即假设每个模态具有独立的阻尼比，通常由实验模态分析获得。

系统类型	方程形式	求解方法	典型应用
2自由度	2×2矩阵方程	解析法（求特征值）	动力吸振器、双级隔振
10~100自由度	稀疏矩阵方程	模态叠加法/子空间迭代	设备-基础耦合、多支撑管路
1000+自由度	大规模有限元	Lanczos法/Block Krylov	整机结构、复杂装配体

2.2.3 案例拆解：双级隔振系统设计

某船舶柴油发电机组产生宽频振动（20Hz~500Hz），单级隔振无法同时满足低频隔振效率和高频衰减要求。采用双级隔振：机组（m₁=2000kg）→一级隔振器（k₁=4×10⁶ N/m，ζ₁=0.05）→中间质量块（m₂=500kg）→二级隔振器（k₂=1×10⁶ N/m，ζ₂=0.05）→船体基础。

分析过程：

1. 建立2自由度方程，质量矩阵[M]=diag(2000, 500)，刚度矩阵[K]=[[4×10⁶, -4×10⁶],[-4×10⁶, 5×10⁶]]
2. 求解特征值问题，得到两阶固有频率：f₁≈8.2Hz，f₂≈18.7Hz
3. 计算频率比：对于20Hz激励，r₁=20/8.2=2.44，r₂=20/18.7=1.07（接近共振）
4. 计算传递率：T_total = T₁×T₂ ≈ 0.18×1.2 ≈ 0.22（20Hz处）
5. 对于100Hz激励：r₁=12.2，r₂=5.35，T_total ≈ 0.007×0.04 ≈ 0.00028

结果表明，双级隔振在100Hz以上具有极佳的衰减效果（衰减约71dB），但在20Hz附近由于第二阶模态接近激励频率，传递率反而略有放大。实际设计中需调整m₂或k₂，使第二阶固有频率避开主要激励频段。

2.2.4 常见误区与总结

• 误区一：认为模态越多越好。实际中，高阶模态的响应通常很小（模态质量参与因子低），且数值误差随模态阶数增加而增大。合理截断比盲目追求全模态更有效。
• 误区二：忽略振型的物理意义。振型描述了系统在各固有频率下的变形形状，是判断测点位置、激励方向是否有效的关键依据。例如，若激励力恰好作用在某一阶振型的节点上，则该阶模态不会被激发。
• 误区三：阻尼矩阵随意假设。比例阻尼只适用于阻尼分布均匀的系统。对于含有集中阻尼器（如减振器）的系统，必须使用非比例阻尼模型。

总结：MDOF系统的核心价值在于揭示系统的模态结构——哪些频率容易共振、哪些位置响应大、哪些模态对特定激励敏感。模态叠加法将复杂耦合问题分解为多个独立的SDOF问题，是工程分析最有力的工具。

2.3 连续系统：从偏微分方程到离散化

当结构的几何尺寸与振动波长可比时（如长轴、薄板、管路），必须将其视为连续系统。连续系统具有无限多个自由度，其运动由偏微分方程（PDE）描述。虽然解析求解PDE仅在少数简单边界条件下可行，但连续系统的模态分析为有限元法等数值方法提供了理论基础和验证基准。

2.3.1 核心原理：典型连续体的振动方程

以工程中最常见的欧拉-伯努利梁为例，其横向自由振动方程为：

EI·(∂⁴y/∂x⁴) + ρA·(∂²y/∂t²) = 0

其中E为弹性模量，I为截面惯性矩，ρ为密度，A为截面积，y(x,t)为横向位移。该方程是四阶空间导数和二阶时间导数的PDE。采用分离变量法，令y(x,t) = Y(x)·T(t)，代入后得到两个常微分方程：空间方程Y'''' - β⁴Y = 0（β⁴ = ρAω²/EI）和时间方程T̈ + ω²T = 0。空间方程的通解为：

Y(x) = C₁·sin(βx) + C₂·cos(βx) + C₃·sinh(βx) + C₄·cosh(βx)

代入边界条件（如简支、固支、自由）后，得到特征方程，其根βₙ对应固有频率ωₙ = βₙ²·√(EI/ρA)。例如，简支梁的固有频率为ωₙ = (nπ)²·√(EI/ρA·L⁴)，n=1,2,3,...

对于轴向振动（杆）和扭转振动（轴），方程分别为：

轴向: EA·(∂²u/∂x²) = ρA·(∂²u/∂t²)
扭转: GJ·(∂²θ/∂x²) = ρJ·(∂²θ/∂t²)

这些方程均为波动方程形式，其解为行波或驻波，取决于边界条件。

2.3.2 判断方法：何时必须使用连续模型

工程判断准则：当结构的最小特征尺寸（如梁的长度、板的跨度）小于振动波长的1/4时，连续效应显著，SDOF/MDOF模型误差超过10%。具体而言：

• 对于梁类结构，若激励频率超过一阶固有频率的3倍，则高阶模态贡献不可忽略，需使用连续模型或高精度有限元。
• 对于管路中的流体脉动，若声波波长与管长可比，必须采用分布参数模型。
• 对于旋转轴，当工作转速接近或超过一阶临界转速时，必须考虑轴的连续弹性。

连续系统的求解策略有两种：一是直接求解PDE（仅限简单边界），二是离散化——将连续体分割为有限个单元，用多项式插值近似位移场，得到近似的MDOF方程。后者即有限元法（FEM）的基本思想。

2.3.3 案例拆解：悬臂梁式传感器的振动分析

某压电式加速度传感器内部的核心敏感元件为悬臂梁结构，梁长L=15mm，宽b=5mm，厚h=0.3mm，材料为铍青铜（E=130GPa，ρ=8500kg/m³）。要求确定其基频，并评估在10kHz激励下的响应特性。

分析过程：

1. 计算截面惯性矩：I = b·h³/12 = 5×10⁻³×(0.3×10⁻³)³/12 = 1.125×10⁻¹⁴ m⁴
2. 计算单位长度质量：ρA = 8500×(5×10⁻³×0.3×10⁻³) = 0.01275 kg/m
3. 悬臂梁基频公式：f₁ = (1.875²/(2π))·√(EI/(ρA·L⁴)) = (3.516/(2π))·√(130×10⁹×1.125×10⁻¹⁴/(0.01275×(0.015)⁴))
4. 计算得：f₁ ≈ 8.2 kHz

该基频接近10kHz激励频率，意味着传感器在10kHz附近将产生明显共振，导致输出信号失真。改进方案：将梁厚增加至0.4mm（I增大至2.67×10⁻¹⁴ m⁴），则f₁ ≈ 12.6kHz，避开10kHz工作频段。

2.3.4 常见误区与总结

• 误区一：将连续系统简单等效为SDOF。等效质量法（如将梁的分布质量折算为端部集中质量）仅在基频附近有效，无法预测高阶模态响应。
• 误区二：忽略边界条件的影响。同一根梁，简支与固支的基频相差约2.47倍（简支f₁=π²/2π·√(EI/ρAL⁴)，固支f₁=3.516/2π·√(EI/ρAL⁴)）。边界条件错误将导致完全错误的预测。
• 误区三：认为有限元法可以完全替代解析解。有限元法需要网格收敛性验证，且对高阶模态的精度依赖于单元类型和网格密度。解析解是验证有限元模型正确性的“金标准”。

总结：连续系统理论揭示了振动在空间中的分布规律——模态形状、节点位置、波传播特性。虽然工程中大量使用有限元软件，但理解PDE的物理内涵和解析解的结构，有助于正确设置边界条件、评估离散误差、解读计算结果。