3、有限元方法入门：显式动力学分析原理、时间步长与稳定性
好，咱们进入正题。这一章聊的是显式动力学分析，说白了就是「怎么算冲击」。
我刚开始接触这个的时候，也觉得很玄乎。什么时间步长、稳定性，听着就头大。但后来做了几个项目，踩过坑，才慢慢摸到门道。今天我把这些经验掰开揉碎了讲给你听。
3.1 显式动力学到底在算什么？
先问一个问题：一个子弹打在钢板上，0.1毫秒内发生了什么？
传统静力学算不了，因为时间太短，惯性力、波传播这些效应都跑出来了。这时候就得用显式动力学。
它的核心思想其实很简单：把时间切成极小的片段，在每个片段里，用上一时刻的状态推算下一时刻。
公式长这样：
a(t) = M⁻¹ · [F_ext(t) - F_int(t)]
v(t+Δt) = v(t) + a(t) · Δt
u(t+Δt) = u(t) + v(t+Δt) · Δt你看，没有迭代，没有求解方程组。一步到位。这就是「显式」的含义——当前步的结果只依赖上一步的已知量。
关键点：显式动力学不需要组装全局刚度矩阵，也不需要求解大型线性方程组。它是一步一步「推」出来的。
我在做装甲板抗冲击仿真时，模型有几十万个单元。如果用隐式算法，一个时间步可能要算几个小时。换成显式，虽然步数多，但每步极快，整体反而高效得多。
3.2 时间步长——为什么不能随便设？
你可能会想：时间步长设小一点，算得准；设大一点，算得快。那我取个中间值不就行了？
嗯，这里要注意。显式动力学的时间步长不是你想设多少就设多少的。
它有一个上限，叫「临界时间步长」。超过这个值，计算就会发散，结果直接飞掉。
临界时间步长的公式：
Δt_critical = L_min / c其中：
- L_min：网格中最小的单元特征尺寸
- c：材料中的波速（纵波速度）
波速c的计算：
c = √(E / ρ)E是弹性模量，ρ是密度。材料越硬、越轻，波速越快，临界时间步长就越小。
警告：实际计算时，一般取Δt = 0.5~0.8倍的Δt_critical。别卡着上限算，容易翻车。我曾经有一次为了省时间，设了0.95倍，结果算了两个小时，最后一步直接炸了——所有节点速度变成NaN。白干了。
3.3 稳定性——显式动力学的命根子
显式动力学是条件稳定的。什么意思？就是时间步长必须小于某个值，否则误差会指数级增长。
为什么会这样？
你可以把显式算法想象成一个「滚雪球」的过程。每步都有微小的截断误差。如果时间步长合适，这些误差会被「消化」掉。但如果步长太大，误差就会像滚雪球一样越滚越大，最后整个系统崩溃。
影响稳定性的因素：
- 网格质量：细长单元、畸形单元会大幅降低临界步长
- 材料刚度：越硬的材料，波速越快，步长越小
- 密度：密度越小，波速越快，步长越小
技巧：如果模型中有非常小的单元导致步长过小，可以考虑使用「质量缩放」。就是人为增加这些单元的密度，降低波速，从而增大步长。但注意——增加的质量不能超过总质量的5%，否则惯性效应失真。
3.4 知识体系框架
下面这张图，是我自己总结的显式动力学核心逻辑。你看一遍，基本就通了。
3.5 实战中的避坑指南
讲几个我踩过的坑，你记一下：
- 网格尺寸突变：大网格和小网格挨在一起，小网格决定了全局步长。我曾经有个模型，99%的单元尺寸是5mm，就几个倒角处是0.1mm。结果步长被拖慢了50倍。后来我把那几个小网格重新划分，效率直接起飞。
- 材料参数错误：显式动力学对材料参数很敏感。尤其是密度和弹性模量，差一个数量级，步长就差一个数量级。我建议你每次建模前，先手动算一下波速，心里有个底。
- 初始穿透：接触定义时，如果初始状态就有穿透，显式算法会在第一步产生巨大的接触力，直接把模型「弹飞」。检查初始接触状态，这是基本功。
总结一句话：显式动力学，时间步长是命，网格质量是根。步长选不对，算再多也是白费。
好了，这一章就到这里。你先把显式动力学的原理和步长概念消化掉，后面我们讲具体怎么在软件里操作。