2. 结构动力学基础

在进入 ANSYS 模态与谐响应分析之前，必须先建立一套能描述“结构如何振动”的数学语言。有限元软件求解的并不是抽象的公式，而是这些方程在离散网格上的数值形式。理解单自由度（SDOF）与多自由度（MDOF）系统的运动方程、固有频率、振型、阻尼和频响函数，是判断仿真结果是否合理、设置是否正确的前提。

一、单自由度系统的运动方程

最简单的振动模型是质量—弹簧—阻尼器系统：质量块 m 只能沿一个方向平动，弹簧提供恢复力 kx，阻尼器提供与速度成正比的阻力 cẋ，外部激励为 F(t)。根据牛顿第二定律，运动方程为：

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

其中 x 为位移，ẋ 为速度，ẍ 为加速度。若忽略阻尼和外部激励，得到无阻尼自由振动方程 mẍ + kx = 0，其解为简谐形式，圆频率为：

ωn = √(k/m)   （rad/s）
fn = ωn / (2π) （Hz）

例如，一个质量 2 kg、刚度 800 N/m 的单自由度系统，其固有圆频率为 √(800/2)=20 rad/s，固有频率约为 3.18 Hz。这个数值只由质量和刚度决定，与初始条件无关，是系统本身的“身份特征”。

二、多自由度系统与矩阵方程

真实结构由无数质点组成，有限元将其离散为有限个自由度。此时运动方程扩展为矩阵形式：

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

其中 [M]、[C]、[K] 分别为质量矩阵、阻尼矩阵和刚度矩阵，{x} 为节点位移向量。忽略阻尼和载荷后，自由振动方程为 [M]{ẍ}+[K]{x}=0。假设解为 {x}={φ}·sin(ωt)，代入后得到广义特征值问题：

([K] − ω²[M])·{φ} = 0

该方程有非零解的条件是系数矩阵行列式为零，即 det([K]−ω²[M])=0。求解可得到 n 个特征值 ω₁²、ω₂²……ωₙ²，对应 n 个固有频率；每个特征值对应一个特征向量 {φᵢ}，即第 i 阶振型。在 ANSYS 中，这一步就是模态分析的核心求解过程。

三、固有频率与振型的物理含义

固有频率描述结构“倾向于以多快的节奏振动”，振型描述结构“以什么形状振动”。以一根两端固定的梁为例，一阶振型通常呈半个正弦波，中间振幅最大；二阶振型呈一个完整正弦波，中间存在一个节点（位移为零的位置）。振型不是绝对位移，而是各自由度位移的相对比例，因此可以任意缩放。ANSYS 输出的振型通常按质量矩阵归一化，使 {φ}ᵀ[M]{φ}=1。

振型之间具有正交性：{φᵢ}ᵀ[M]{φⱼ}=0，{φᵢ}ᵀ[K]{φⱼ}=0（i≠j）。这一性质使得多自由度系统可以通过模态叠加法解耦，也是谐响应分析采用模态叠加法的基础。

四、阻尼的类型与影响

阻尼耗散振动能量，使自由振动逐渐衰减，并在共振区限制响应幅值。结构动力学中常见的阻尼模型有三类：

• 粘性阻尼：阻尼力与速度成正比，F_d=cẋ。数学处理最方便，ANSYS 中常用。
• 结构阻尼（滞回阻尼）：阻尼力与位移成正比、与速度同相，用损耗因子 η 或结构阻尼比 g 描述，适合金属材料。
• 库仑阻尼：由干摩擦产生，阻尼力为常数，方向与运动相反，非线性特征明显。

工程中最常用的是粘性阻尼比 ζ=c/(2√(mk))。ζ 决定系统响应形态：ζ<1 为欠阻尼，自由振动呈衰减振荡；ζ=1 为临界阻尼，无振荡；ζ>1 为过阻尼，缓慢回到平衡位置。大多数金属结构 ζ 在 0.01～0.05 之间，焊接结构约 0.02，橡胶隔振器可达 0.1～0.3。在 ANSYS 模态分析中，阻尼比直接影响谐响应共振峰的幅值：阻尼越小，峰值越高、越尖锐。

五、频响函数（FRF）的定义

对运动方程两边做傅里叶变换，令 s=jω，可得频域关系：

(−ω²[M] + jω[C] + [K])·{X(ω)} = {F(ω)}

定义动刚度矩阵 [Z(ω)]=−ω²[M]+jω[C]+[K]，则 {X(ω)}=[H(ω)]·{F(ω)}，其中 [H(ω)]=[Z(ω)]⁻¹ 即为频响函数矩阵。对单自由度系统，位移频响函数为：

H(ω) = 1 / (k − mω² + jωc)

其幅值在 ω≈ωn 处达到峰值，峰值高度约为 1/(2ζk)。频响函数描述“单位简谐激励下系统的稳态响应”，是谐响应分析和实验模态分析的核心量。ANSYS 谐响应分析本质上就是求解 [H(ω)] 在给定频率点上的值，再乘以实际载荷谱得到响应。

六、从理论到 ANSYS 的对应关系

将上述概念映射到 ANSYS Workbench 操作中：材料密度和弹性模量决定 [M] 和 [K]；约束条件决定哪些自由度被保留；求解器提取的“频率”对应 ωᵢ/2π；振型云图对应 {φᵢ}；阻尼设置在“Analysis Settings”中通过恒定阻尼比或材料阻尼输入；谐响应分析的扫频范围应覆盖关注阶次的固有频率。理解这些对应关系，才能在软件中正确设置参数并判断结果是否可信。

掌握这些基础后，后续的 ANSYS 操作就不再是“点按钮”，而是有明确物理意义的求解过程：建模是在构造 [M] 和 [K]，网格是在离散连续体，求解是在解特征值问题，后处理是在读取振型和频响。只有把软件操作与运动方程对应起来，才能真正理解模态与谐响应分析的结果。