动力分析（Dynamic Analysis） — CAE术语解说
动力分析概述
动力分析和静力分析有什么区别？静力分析也可以施加荷载看变形，对吗？
静力分析假设荷载缓慢施加，加速度为零，因此忽略质量×加速度的惯性力。但实际结构会承受振动、冲击等，加速度不可忽视。这时惯性力变得重要，需要用动力分析来处理。
具体在什么场景下需要动力分析？
比如汽车碰撞试验、建筑地震响应、电子设备坠落冲击、发动机旋转振动——这些都涉及"随时间变化的荷载"问题。实务中当有人问"这个荷载真的是静态的吗？"时，就需要考虑动力分析。
支配方程
动力分析的基本方程是什么？
结构动力分析的根本是以下运动方程：
$$\mathbf{M}\ddot{\mathbf{u}} + \mathbf{C}\dot{\mathbf{u}} + \mathbf{K}\mathbf{u} = \mathbf{F}(t)$$
$\mathbf{M}$、$\mathbf{C}$、$\mathbf{K}$ 分别代表什么？
$\mathbf{M}$ 是质量矩阵（惯性大小），$\mathbf{C}$ 是阻尼矩阵（能量耗散），$\mathbf{K}$ 是刚性矩阵（弹簧硬度）。右边的 $\mathbf{F}(t)$ 是随时间变化的外力。静力分析可以看作是将 $\mathbf{M}\ddot{\mathbf{u}}$ 和 $\mathbf{C}\dot{\mathbf{u}}$ 置零的特殊情况。
阻尼 $\mathbf{C}$ 怎样确定？材料数据表里没有啊。
实务中常用Rayleigh阻尼 $\mathbf{C} = \alpha\mathbf{M} + \beta\mathbf{K}$ 来近似。$\alpha$ 和 $\beta$ 通过目标频率段的阻尼比（钢铁通常1~2%，橡胶5~10%）反推计算。没有实验数据就很难精确确定，这是个棘手的问题。
模态分析（固有值分析）
模态分析是求什么的？
在无外力、无阻尼的条件下，求解 $\mathbf{M}\ddot{\mathbf{u}} + \mathbf{K}\mathbf{u} = \mathbf{0}$，得到固有振动数（以何种频率共振）和振型（该频率下的振动形状）。这就像结构物的"振动指纹"。
固有值问题表述为：
$$(\mathbf{K} - \omega_i^2\mathbf{M})\boldsymbol{\phi}_i = \mathbf{0}$$
其中 $\omega_i$ 是第 $i$ 阶模态的固有圆频率，$\boldsymbol{\phi}_i$ 是对应的振型向量。固有振动数为 $f_i = \omega_i / (2\pi)$ [Hz]。
知道固有振动数有什么用？
比如发动机转速3000转/分对应激励频率50Hz。如果结构一阶固有频率也是50Hz，就会发生共振导致损坏。所以模态分析用来"这个频率段有固有频率吗？"来评估安全性。汽车、飞机、家电——各行业的首个动力分析通常都是模态分析。
振型有什么用处？
它告诉你哪些部位变形最大，对于确定加强位置很有帮助。比如汽车仪表板抖动，看振型就能发现"这个位置变形最大"，然后在那里加肋板来增加刚度。
频率响应分析（谐波响应）
频率响应分析和模态分析怎样不同？
模态分析是"无外力"求固有频率。频率响应分析（Harmonic Response）是"施加正弦波外力 $\mathbf{F}(t) = \mathbf{F}_0 e^{i\omega t}$，求各频率下的响应"。也就是计算输入输出关系，或称传递函数。
听说有模态叠加法和直接法两种？
模态叠加法（Modal Superposition）先用模态分析的结果解耦方程，计算速度快，但有模态截断误差。直接法（Direct Method）每个频率步都直接求解原方程组，精度高但计算量大。宽频带用模态法，局部细节用直接法，根据需求选择。
实际应用的例子？
旋转机械不平衡响应是典型例子。通过改变转速来寻找"哪个转速时振动最大"和"轴承部位振幅是否超限"。汽车NVH（Noise, Vibration, Harshness）分析中也用频率响应来看车内声压在不同频率的分布。
瞬态响应分析（时间积分）
瞬态响应分析是干什么用的？
在时间域逐步积分运动方程，直接求解"时刻 $t$ 时的位移、速度、加速度、应力"。冲击荷载、脉冲荷载、地震波这种无法用正弦波描述的任意时刻历荷载都能处理。
时间步长 $\Delta t$ 怎样确定？太小了计算时间会爆炸。
一般经验是，对目标最高频率 $f_{max}$，时间步应满足 $\Delta t \le T/20$，其中 $T = 1/f_{max}$。比如要追踪1000Hz的信号，则 $\Delta t \le 0.05\,\text{ms}$。显式法还受Courant条件约束。实务中通常边看结果边调，确保荷载波形没有被"切掉"。
比如手机坠落冲撃，用瞬态响应分析吗？
完全是。手机坠落试验、汽车碰撞模拟都属于瞬态响应分析。毫秒级现象用显式法，地震这样几十秒的现象用隐式法。这是下一个话题。
随机振动分析
随机振动分析听起来很复杂，是什么意思？
当荷载无法确定性描述、只能用统计来表示时使用。输入定义为功率谱密度（PSD），输出为RMS值（均方根）。比如运输途中的卡车床随机振动、飞机乱流响应——无法预测"何时会有什么波形"的环境评估。
输出只有RMS值，最大值怎样评估？
假设高斯分布，用3σ值（RMS的3倍）作为最大值估计。这样能覆盖99.7%的概率。MIL规格和JAXA环境试验标准都采用这个方法。也可以用Miles方程做单自由度系统的快速估算，但多自由度系统必须用有限元正确求解。
显式法和隐式法
显式法（Explicit）和隐式法（Implicit）怎样选择？
简单来说，差别如下：
| 项目 | 显式法 (Explicit) | 隐式法 (Implicit) |
|---|---|---|
| 代表求解器 | LS-DYNA, Abaqus/Explicit, PAM-CRASH | Nastran, Abaqus/Standard, ANSYS Mechanical |
| 时间步长 | 极小（Courant条件 $\Delta t \le L_{min}/c$） | 相对较大（无条件稳定） |
| 单步计算成本 | 小（无需求解线性方程组） | 大（每步求解线性方程组） |
| 适用现象 | 碰撞、爆炸、贯穿（~毫秒） | 振动、地震、热过渡（~秒及以上） |
| 稳定性 | 条件稳定 | 无条件稳定（Newmark-β等） |
Courant条件是什么？
显式法的稳定条件。时间步必须不超过最小单元尺寸 $L_{min}$ 除以声速 $c$ 的值：
$$\Delta t \le \frac{L_{min}}{c}, \qquad c = \sqrt{\frac{E}{\rho}}$$
钢的 $c \approx 5000\,\text{m/s}$，单元1毫米的话 $\Delta t \le 0.2\,\mu\text{s}$？太小了！
正是。所以显式法解长时间现象时，时间步数会天文数字，不现实。反之，10毫秒碰撞用显式法只需5万步左右，完全可行。关键是"现象时间尺度"和"需要的时间步"的平衡。
Newmark-β法经常听说，是什么算法？
隐式法的代表时间积分格式。通过参数 $\beta$ 和 $\gamma$ 的选择改变特性。$\beta = 1/4$, $\gamma = 1/2$ 是平均加速度法（无条件稳定、二阶精度），$\beta = 1/6$, $\gamma = 1/2$ 是线性加速度法。Nastran内部也用Newmark类算法。
实务应用总结
到底怎样选择用哪种动力分析方法？能总结一下吗？
这样判断比较好：
- 「要确认是否会共振」 → 模态分析
- 「想知道某频率激励下的响应」 → 频率响应分析
- 「要输入地震波或冲击时刻历」 → 瞬态响应分析
- 「要评估统计振动环境」 → 随机振动分析
- 「碰撞、爆炸、毫秒级超高速现象」 → 显式法瞬态响应
- 「振动、地震、秒级及以上现象」 → 隐式法瞬态响应
先从模态分析开始，再逐步进展到其他分析，这样吗？
完全同意。模态分析是动力分析的基础，计算成本也小。掌握固有频率和振型后再进行频率响应或随机振动，这是标准做法。不要直接跳到瞬态响应分析，否则结果很难解释。
准确理解CAE术语是团队沟通的基础。——Project NovaSolver也致力于支持实务工程师的学习。
动力分析实务中的课题，请告诉我们
Project NovaSolver致力于解决CAE工程师日常面临的课题——建模复杂性、计算成本、结果解释——的解决方案。您的实务经验是开发更好工具的源动力。
联系我们（筹备中）相关话题
价值
详细
错误