2. 机械系统动力学建模基础

要对一台机械系统做模态与频响分析，第一步不是打开软件，而是把真实结构抽象成可计算的动力学模型。这一步决定了后续所有频率、振型与响应结果是否可信。本节围绕单自由度（SDOF）与多自由度（MDOF）系统的运动方程、质量/刚度/阻尼矩阵以及模态叠加原理展开，为模态分析和频响计算建立统一的数学描述。

一、单自由度系统：最小可用的动力学模型

把结构简化为一个集中质量 m、一根等效弹簧 k 和一个等效阻尼器 c，外力 f(t) 作用于质量块，位移 x(t) 从静平衡位置量起，向下为正。由牛顿第二定律可得单自由度运动方程：

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

其中 m 的单位为 kg，k 为 N/m，c 为 N·s/m。无阻尼自由振动时令 f=0、c=0，得到固有圆频率 ω_n=√(k/m)，固有频率 f_n=ω_n/(2π)。例如 m=2 kg、k=8000 N/m，则 ω_n=63.2 rad/s，f_n≈10.1 Hz。阻尼比 ζ=c/(2√(mk))，当 c=8 N·s/m 时 ζ≈0.063，属于小阻尼系统，共振峰尖锐。

单自由度模型的价值在于：它给出了“质量—刚度—阻尼”三要素与固有频率、共振放大之间的定量关系。工程中常把复杂结构在某阶模态附近等效为单自由度系统，用于快速估算共振转速或隔振效果。

二、多自由度系统：矩阵形式的运动方程

当结构存在多个明显振型或需要描述空间分布时，必须用多自由度模型。设系统有 n 个自由度，位移向量 x=[x₁, x₂, …, x_n]^T，外力向量 f(t)，则运动方程为：

M·ẍ(t) + C·ẋ(t) + K·x(t) = f(t)

其中 M、C、K 分别为 n×n 的质量、阻尼和刚度矩阵。以两自由度弹簧—质量系统为例，质量块 m₁、m₂ 由弹簧 k₁、k₂、k₃ 串联连接，则：

M = [m₁  0 ]      K = [k₁+k₂   -k₂  ]
    [0   m₂]          [-k₂    k₂+k₃]

矩阵元素有明确物理意义：M 对角元为各自由度集中质量，非对角元表示惯性耦合；K 对角元为某自由度位移引起的恢复力系数，非对角元为自由度之间的弹性耦合，负号表示相邻质量反向运动时弹簧被拉伸。阻尼矩阵 C 通常难以直接写出，工程上常用比例阻尼假设 C=αM+βK，其中 α、β 由实验模态阻尼比拟合得到。

三、模态叠加原理：把耦合方程解耦

无阻尼自由振动时，令 f=0、C=0，设解为 x=φe^jωt，代入得广义特征值问题：

(K − ω²M)·φ = 0

求解得到 n 个特征值 ω_r² 和对应特征向量 φ_r，即第 r 阶固有频率与振型。振型具有正交性：φ_r^TMφ_s=0、φ_r^TKφ_s=0（r≠s）。利用这一性质，把物理坐标 x 展开为各阶模态的线性组合：

x(t) = Σ φ_r · q_r(t) = Φ · q(t)

其中 Φ 为模态矩阵，q 为模态坐标。将展开式代入运动方程并左乘 Φ^T，在比例阻尼下可得到一组相互独立的单自由度方程：

m_r·q̈_r + c_r·q̇_r + k_r·q_r = F_r(t)

这里 m_r=φ_r^TMφ_r 为模态质量，k_r=φ_r^TKφ_r 为模态刚度，F_r=φ_r^Tf 为模态力。每个模态坐标 q_r 可独立求解，再叠加回物理坐标，这就是模态叠加法的核心。

四、从方程到频响：建模流程与判断依据

实际建模时，可按以下顺序推进：确定边界条件与自由度 → 建立 M、K 矩阵 → 求解特征值得到 ω_r、φ_r → 拟合阻尼 → 用模态叠加计算频响函数。频响函数 H_ij(ω) 表示在 j 点激励、i 点响应的传递关系，其模态展开形式为：

H_ij(ω) = Σ [ φ_ir · φ_jr / (m_r · (ω_r² − ω² + j·2ζ_r·ω_r·ω)) ]

观察频响曲线时，峰值位置对应固有频率，峰值高度受模态阻尼比控制，峰的数量反映关注频段内的模态阶数。若某阶模态在激励点或响应点振型值接近零，则该阶在频响中可能“消失”，这正是模态实验中需要多点激励、多点测量的原因。

建模阶段还需注意：质量矩阵应包含结构质量与附加质量（如配重、液体），刚度矩阵应反映实际连接刚度（螺栓连接、焊缝、轴承支承），边界条件要与实际安装状态一致。若模型固有频率与实测偏差超过 5%～10%，应优先检查连接刚度与边界约束，而非直接调整阻尼。

五、建模中的关键约定与常见误区

坐标系与符号约定必须在建模开始前明确并全文统一。本节约定：位移 x 以静平衡位置为零点，向下为正；弹簧力以拉伸为正；阻尼力与速度方向相反。若改用向上为正，则重力项与弹簧预压缩量需重新处理，否则固有频率计算会出现偏差。

常见误区有三类。其一，把集中质量模型直接用于高频段，当激励频率接近或超过结构局部模态时，单自由度或少量自由度模型失效，需要细化网格或增加自由度。其二，忽略连接刚度，把螺栓连接简化为刚性连接，导致计算频率偏高；实际连接刚度可通过锤击试验识别。其三，阻尼矩阵随意赋值，比例阻尼系数 α、β 应由实测模态阻尼比反算，而非凭经验填写。

掌握上述方程、矩阵与叠加原理后，后续的模态分析（求解特征值）与频响分析（计算 H(ω)）就有了统一的数学基础。建模阶段多花时间核对质量、刚度与边界，远比在结果异常时反复调整参数更有效。