2. 瞬态动力学控制方程与求解原理

瞬态结构动力学要回答的问题，本质上是：一个结构在随时间变化的载荷作用下，每一个时刻的位移、速度、加速度分别是多少。与静力学只需求解 K·u = F 不同，瞬态分析必须同时考虑惯性效应和阻尼耗散，因此控制方程是一个关于时间的二阶常微分方程组：

M·ü(t) + C·u̇(t) + K·u(t) = F(t)

其中 M 为质量矩阵，C 为阻尼矩阵，K 为刚度矩阵，ü、u̇、u 分别是节点加速度、速度、位移向量，F(t) 是随时间变化的外载荷向量。理解这三个矩阵的物理意义，是理解整个求解过程的前提。

三个矩阵分别代表什么

质量矩阵 M 描述结构的惯性。它决定了结构抵抗速度变化的程度，也决定了系统的固有频率。在 ANSYS 中，质量矩阵通常由单元密度积分得到，可采用一致质量矩阵（consistent mass）或集中质量矩阵（lumped mass）。一致质量矩阵精度更高但计算量更大；集中质量矩阵为对角阵，求逆快，常用于显式方法。质量矩阵直接决定惯性项 M·ü 的大小，当载荷变化极快时，惯性项不可忽略，这正是瞬态分析与静力分析的根本区别。

阻尼矩阵 C 描述能量耗散。真实结构中的阻尼来源复杂（材料内摩擦、连接面滑移、空气阻力等），难以精确建模。工程上常用 Rayleigh 阻尼，即把阻尼矩阵表示为质量矩阵与刚度矩阵的线性组合：

C = α·M + β·K

其中 α 为质量比例阻尼系数，β 为刚度比例阻尼系数。α 主要抑制低频响应，β 主要抑制高频响应。在 ANSYS 中通过 ALPHAD 和 BETAD 两个实常数输入。若给定两个目标频率 ω₁、ω₂ 及对应阻尼比 ξ₁、ξ₂，可反解出 α、β。例如取 ω₁=2π·10 rad/s、ω₂=2π·100 rad/s、ξ₁=ξ₂=0.02，则 α≈2.28、β≈5.8e-5（示例计算值，实际需按结构频率确定）。

刚度矩阵 K 描述结构抵抗变形的能力，与静力分析中的刚度矩阵一致，由单元刚度组装而成。它决定了弹性恢复力 K·u。当结构发生大变形或接触状态变化时，K 可能随位移变化，此时进入非线性瞬态分析。

时间积分：把连续时间离散成步

控制方程无法解析求解，必须用数值时间积分逐步推进。核心思路是：已知 t 时刻的位移、速度、加速度，通过某种假设推算 t+Δt 时刻的响应。不同算法对加速度在时间步内的变化方式做了不同假设，由此带来稳定性与精度的差异。

中心差分法（显式）：用 t 时刻的平衡方程求解 t+Δt 的位移，加速度用中心差分近似：

ü(t) = [u(t+Δt) - 2u(t) + u(t-Δt)] / Δt²

代入运动方程后可显式解出 u(t+Δt)，无需迭代求解方程组，每步计算量小。但它是有条件稳定的，时间步必须小于临界步长 Δt_cr = 2/ω_max（ω_max 为结构最高阶固有频率）。因此显式方法适合冲击、碰撞、爆炸等短时高频问题，步长通常取 1e-6～1e-7 s 量级。

Newmark 法（隐式）：通过两个参数 γ、β 控制加速度假设：

u̇(t+Δt) = u̇(t) + [(1-γ)ü(t) + γü(t+Δt)]·Δt
u(t+Δt) = u(t) + u̇(t)·Δt + [(0.5-β)ü(t) + βü(t+Δt)]·Δt²

当 γ=1/2、β=1/4 时，即平均加速度法，无条件稳定，是 ANSYS 瞬态分析默认设置。当 γ=1/2、β=1/6 时为线性加速度法，有条件稳定。Newmark 法每步需要求解隐式方程组，计算量大，但步长可以取得较大，适合结构响应持续时间较长的常规瞬态问题。

HHT 算法（α 法）：HHT 在 Newmark 基础上引入参数 α，对刚度项和阻尼项在 t 与 t+Δt 之间做加权，同时保持对惯性项的 t+Δt 时刻取值。其稳定性条件为 α∈[-1/3, 0]，当 α=0 时退化为 Newmark 平均加速度法。HHT 的优势在于可引入数值阻尼，有效抑制高频虚假振荡，同时不显著影响低频精度。在 ANSYS 中通过 TINTP 命令或 TRNOPT 设置，α 通常取 -0.05～-0.1。

稳定性与精度的取舍

稳定性指误差不会随时间步推进而无限放大。显式中心差分法必须满足 CFL 条件，一旦步长超过临界值，结果会迅速发散。隐式 Newmark（β=1/4）和 HHT 无条件稳定，步长不受稳定性限制，但步长过大会损失高频精度。因此实际选择步长时，稳定性只是下限，精度才是决定因素。工程经验是：若关心结构的第 n 阶模态，时间步应小于该阶周期 T_n 的 1/20，即 Δt ≤ T_n/20。例如关心 100 Hz 响应，T=0.01 s，则 Δt 应不大于 5e-4 s。

精度方面，Newmark 平均加速度法在低频段精度良好，但会引入周期性数值振荡；HHT 通过 α 参数引入可控数值阻尼，能有效衰减这些非物理高频振荡，因此在 ANSYS 中常作为默认推荐。中心差分法在满足步长条件时精度也较高，但步长限制使其更适合短时问题。

总结来说，瞬态动力学的求解过程就是：用 M、C、K 三个矩阵描述结构动力学特性，用 F(t) 描述外载，再选择合适的时间积分算法逐步推进。显式中心差分适合短时高频冲击，隐式 Newmark 适合常规结构响应，HHT 则在隐式框架下提供数值阻尼以抑制高频噪声。理解三者的适用条件、稳定性和精度差异，才能在 ANSYS 中做出正确的求解设置。