我入行前两年，一直觉得结构仿真就是“把力加上去，看应力云图”。直到有一次做设备底座的跌落分析，静力算出来的最大应力连屈服强度一半都不到，按静力标准绰绰有余，可样机实测摔了几次就出现裂纹。回头看才意识到，这类问题根本不能当静力算——它必须用 ANSYS Workbench 的瞬态结构分析，把载荷随时间的变化过程完整复现出来。这也是今天这篇要聊的时域分析。
无论你是刚接触 Workbench 的新手，还是平时做静力分析比较多、想往动力学方向再走一步的工程师，这篇内容都能给你一条走得通的路。我的计划是：先用物理直觉说清楚“什么时候必须做瞬态分析”，再讲透时域求解的底层逻辑，然后完整走一遍 Workbench 瞬态分析流程，最后把时间步长这个最关键的命门掰开揉碎，顺带附上我实际踩过的几个坑和排查思路。
1. 为什么静力分析会“算错”：时域分析的应用边界
1.1 静力分析先天的两个假设
静力分析能成立，靠的是两个隐含前提。第一个前提是载荷施加得足够慢，慢到结构在任意瞬间都处于平衡状态，加速度带来的惯性力小到可以忽略；第二个前提是结构各点的响应可以看成同步发生，不用考虑应力波在结构里跑来跑去的时间。
在这两个前提下，一个弹簧压到底，载荷多大位移就多大，不存在“过头”和“滞后”。设备底座螺栓预紧、压力容器受内压、建筑梁柱承受恒载，这些工况用静力分析是又快又准。可一旦对象变成跌落、碰撞、冲击振动、电机启停、地震响应，这两个前提就崩了。
举一个特别直观的例子。你用两只手慢慢压一个弹簧，弹簧的变形完全跟得上手的位移，每个时刻都是平衡的；但你用拳头快速砸上去，弹簧还没反应过来，力已经撤了，位移却还在振荡。静力分析处理不了第二种情况，因为“还没反应过来”这件事，本质上是质量和加速度在起作用。
1.2 一个简单的判据：加载时间与固有周期
那怎么判断一个工况到底能不能用静力分析？我的经验是拿载荷持续时间 (t_p) 和结构固有周期 (T) 比一比。
固有周期可以从模态分析里拿到，一阶频率的倒数就是基频周期。如果 (t_p > 10T)，结构有足够时间“跟上”载荷的变化，可以按准静态处理，用静力分析足够；如果 (t_p) 和 (T) 在同一个量级，甚至 (t_p < T)，那就必须走瞬态，否则会严重低估结构内部的动态放大效应。
拿很多同行做过的机箱跌落来算：跌落冲击的接触脉冲往往只有几毫秒，而机箱壳体的一阶频率大概在几十到几百赫兹，对应周期也是几毫秒到几十毫秒。两者完全重叠，所以跌落分析从来没人敢用静力去做，这是有物理原因在里面的。
1.3 时域和频域怎么分工
有人会问：那既然要解决动态问题，谐响应分析也是频域的，和时域分析有什么区别？
谐响应分析只能算稳态正弦激励，而且假设结构始终按固定频率做简谐振动，线性范围内才能用。它适合对标准正弦振动台试验、旋转机械的简谐激励做响应评估。但很多真实工况不是稳态正弦，而是冲击、突变、启停、非周期载荷，这些只能放到时域里一点点积分。更关键的是，接触、塑性、大变形这些非线性行为在频域里根本无法直接体现，只有时域分析能处理。
所以这里可以划一条很实用的界线：问题里若存在明显的“先后顺序”“冲击过程”“非线性接触演变”，就直接选时域分析；问题只是求稳态简谐响应幅值，才考虑频域工具。
提示：不要因为学了瞬态分析，就什么工况都往上套。单调加载的强度校核，静力分析仍然是成本最低、结果最稳的选择。
2. 动力学方程背后的物理直觉：质量、阻尼和时间积分
2.1 质量、阻尼、刚度到底在方程里干什么
瞬态结构分析求解的核心方程是：
[ [M]{\ddot{u}} + [C]{\dot{u}} + [K]{u} = {F(t)} ]
很多人一看到矩阵和微分符号就头疼，其实每一项的物理含义都很直白。([M]{\ddot{u}}) 是惯性项，加速度越大、质量越大，产生的惯性力越大；它解释了为什么重物砸到地面上的冲击力远大于它本身的重力。([C]{\dot{u}}) 是阻尼项，对应结构内部材料摩擦、连接面滑移、空气阻力等耗能机制；阻尼决定了冲击过后结构会不会一直来回振荡、振荡多久才衰减。([K]{u}) 是刚度项，和静力分析里一样，位移越大弹性回复力越大。
在静力分析里，前两项被直接扔掉了，只剩 ([K]{u}={F})。而时域分析把三项全部拉回来参与力平衡，这就是“动态平衡”的本质。
我比较喜欢用推秋千来记。秋千摆起来之后你不推它，它也能自己来回荡很久，这就是惯性；荡着荡着幅度越来越小，最后停下来，就是阻尼；而你把它拉到旁边，手一松，它想弹回来的那股劲，就是刚度。三者同时存在，才构成一个完整的动态系统。
2.2 瑞利阻尼参数的计算方法
瞬态分析里最难拿到准确数值的参数就是阻尼。很多工程问题的材料阻尼数据一片空白，这时候最常见的处理是用瑞利阻尼，即把阻尼矩阵近似成质量矩阵和刚度矩阵的线性组合：
[ [C] = \alpha[M] + \beta[K] ]
(\alpha) 主要压低频响应，(\beta) 主要压高频响应。知道了目标频带内两个频率对应的模态阻尼比 (\xi)，就能反推：
[ \alpha = \frac{2\xi\omega_1\omega_2}{\omega_1+\omega_2}, \quad \beta = \frac{2\xi}{\omega_1+\omega_2} ]
举个例子。假设某结构关心 10~200Hz 这个频带，取阻尼比 (\xi=0.03)，那么 (\omega_1=62.8\text{ rad/s})、(\omega_2=1256.6\text{ rad/s})，算下来：
[ \alpha = \frac{2 \times 0.03 \times 62.8 \times 1256.6}{62.8 + 1256.6} \approx 3.55, \quad \beta = \frac{0.06}{1319.4} \approx 4.5\times 10^{-5} ]
按这个参数填到 Workbench 的 Damping 里，就能让结构在这段频带内的响应和真实结构大致贴合。金属结构阻尼比一般取 0.01~0.05，不同行业有自己的习惯值；橡胶等阻尼材料则要单独处理，不能简单套用瑞利阻尼。
注意：瑞利阻尼在频带之外的抑制效果可能跟实际差别很大。如果你只取一个频率点来反推参数，低频段或者高频段可能会被过度衰减，所以条件允许时至少要取两个分布合理的频率点。
2.3 直接积分法是如何“一步一个脚印”算出来的
时域分析用的数值方法是直接时间积分。你可以把它理解成把连续时间切成很多离散的时间步，每一步都基于当前时刻的位移、速度、加速度，推算下一个时刻的状态。
Workbench 默认的程序控制求解器采用的就是 Newmark 族算法（HHT），属于隐式积分。隐式积分的显著优点是在一定条件下是无条件稳定的，但代价是每一时间步都要解一次大型方程组。所以时间步长减小一半，计算量大约翻一倍，这个数量级概念一定要有，后面选步长时会反复用到。
有了这个动态系统视角，你再看 Workbench 里那些设置项就不会觉得神秘了。接下来我们进软件里实际操作一遍，把流程走通。
3. Workbench瞬态分析从建模到结果动画：一步步怎么走
3.1 分析类型与材料参数：密度缺失是新手第一大坑
在 Workbench 左侧的 Toolbox 里找到 Transient Structural，拖到项目原理图中，然后关联几何模型。这个地方有个常见误区：很多人沿用静力分析的习惯，材料只填弹性模量和泊松比，把密度给漏了。
瞬态分析里质量矩阵全靠密度生成，密度不填，惯性项全部为零，整个动力学方程就退化成了静力方程，算出来的结果自然毫无意义。所以每次新建瞬态分析，我第一件事就是到 Engineering Data 里确认密度已经填好，同时检查弹性模量、泊松比，非线性问题还要带上屈服强度和切线模量。
单位问题也是重灾区。几何用毫米时，Workbench 内部常见配套单位是力用 N、应力用 MPa、时间是 s，密度则需要换算成 tonne/mm³，也就是把 (\text{kg/m}^3) 的数值除以 (10^9)。钢材密度 7850 (\text{kg/m}^3)，换算后要填 (7.85\times 10^{-9})。Workbench 界面有时会按材料库帮你处理，但一旦你自己手输数据，就很容易在这个数量级上翻车，导致结果差出好几个数量级。
3.2 几何与接触：冲击类问题要格外注意的设置
几何来源可以是从 DesignModeler、SpaceClaim 或外部 CAD 导入。复杂装配体我的建议是先把倒角、小孔、螺纹这类小特征简化掉再进仿真，因为这些特征到了网格划分阶段会逼着你把网格局部化得非常细，而瞬态分析一步一个方程，网格变细之后时间步长往往还得跟着降，计算量直接爆炸。
接触设置是瞬态分析最容易埋雷的地方。绑定接触适合模拟螺栓连接、焊接这类没有相对滑移的场景；摩擦接触和无摩擦接触适合冲击过程中可能发生开缝、滑移的结构。绑定接触会把结构算得过刚，但作为第一轮初步分析完全够用；如果重点就是看接触面的撞击力和应力波传递，那就要按实际接触行为设置。
我遇到过多次因为接触初始穿透导致求解器报错的案例。几何建模时两个面差了一丁点，默认接触探测模式下求解器会认为初始就有穿透，一开算就弹错误。碰到这种情况，先到 Contact 里看 Gap/Penetration 的初始状态，有穿透就调整接触偏移或修正几何，不要硬着头皮求解。
3.3 网格尺寸与单元类型选择的实际考量
网格策略在瞬态分析里比静力分析更重要，因为它和时间步长是联动的。我的习惯是：整体网格可以适当粗一些，但在冲击点、应力集中区域、连接根部这些关键部位必须局部加密；单元类型优先六面体，复杂几何用四面体也能跑，只是计算效率和精度略逊一筹。
这里要引入一个波长概念。结构响应中你想捕捉高频成分，网格就必须能分辨对应波长的空间分布，规则是单元尺寸要远小于波长。钢材中的纵波波速大约 5000 m/s，如果关心 1000 Hz 的振动成分，对应波长 5 m，那 0.1 m 的网格就绰绰有余；但如果关心 100 kHz 的应力波，波长只剩 50 mm，网格就要细到 5 mm 以下。
这也是很多新手做跌落仿真时结果跳变、振荡得离谱的原因：网格太粗，高频应力波被数值方法直接过滤掉或者严重畸变。
3.4 载荷表格与初始条件：让时间概念进入模型
在静力分析里，载荷就是“把力加上去”；在瞬态分析里，载荷是一张时间表。Workbench 的载荷定义窗口里有 Tabular Data，你可以直接输入时间-载荷曲线。比如 0 到 0.001 s 线性上升到 1000 N，0.001 到 0.002 s 保持 1000 N，0.002 s 后卸载为 0，那就需要在这几个转折点都定义数据行。
如果载荷历程来自试验实测的 CSV 文件，也可以直接把两列数据复制进表格，Workbench 会自动按线性插值处理。需要特别注意的是，载荷突变的地方一定要手动加密时间点，否则插值会把这个突变拉平，冲击效果就完全没了。
初始条件容易被忽略。默认情况下结构从静止零应力开始，这符合大多数工况。但跌落问题通常要给定初始速度：在 Environment 上右键插入 Initial Condition，给零件赋一个速度初值，模拟它在接触瞬间已经带有落地速度的状态；或者借助后续的时间步把重力、接触一点点加载进去。这块内容展开比较多，我计划在后面的系列文章里单独写，本篇先把概念建立起来。
3.5 求解设置与后处理：动画、探针与时间历程
Analysis Settings 里最重要的几个参数是 Number of Steps、Step End Time 和 Time Step。简单场景一个载荷步从头到尾算完，复杂载荷历程则要拆成多个阶段，每个阶段对应一个 Step，设定不同的结束时间和子步策略。这个我下一章详细讲。
求解完成后，后处理不要只盯着最大应力云图。瞬态分析的黄金判断标准是看响应随时间的变化过程：插入 Total Deformation、Equivalent Stress，然后播放动画，看结构是否像实际物体那样有波动的变形状态；再插入一个 Probe 探针，放在你关心的节点或面上，绘制这个位置的时间历程曲线，看响应峰值在哪个时刻出现、衰减趋势对不对。
如果动画里看到应力在相邻单元之间呈棋盘状跳变，或者曲线在高频区间疯狂抖动，多半是时间步长或者网格的问题，不是真实物理现象。这类细节靠照本宣科学不到，得自己多算几个案例才有感觉。
4. 时间步长取多少：时域分析真正的分水岭
4.1 时间步长的四条约束，按场景选
我把确定时间步长的方法总结成一张表，实战中照着选就行：
| 约束类型 | 适用场景 | 经验规则 |
|---|---|---|
| 载荷频率约束 | 所有瞬态分析 | 载荷中最高频率成分，每周期至少 10~20 步 |
| 响应频率约束 | 关心结构固有响应 | 关注到第几阶模态，该阶周期内至少 10~20 步 |
| 波传播约束 | 高速冲击、应力波传播明显 | 时间步长小于最小网格尺寸除以波速 |
| 数值稳定性约束 | 隐式求解器 | 程序自动控制，一般不用手动干预 |
载荷频率约束怎么理解？如果冲击脉冲宽度是 2 ms，它的主要频谱大体覆盖到几百赫兹。等效最高频率按脉冲宽度的倒数按量级估计，再乘以 20 作为采样率要求，就能反推出时间步长。脉冲越窄，频谱越宽，需要的步长就越小，计算越昂贵。
响应频率约束是同一个道理，只是看结构本身。关心到第几阶模态，就要保证那阶模态周期内有足够多的时间点。与其事后发现曲线毛刺太多，不如事前先估算。
4.2 一个具体算例，把步长算给你看
假设一块 500 mm × 200 mm 的钢制薄板，撞击点处网格细化到 5 mm，第一阶固有频率大约 50 Hz，冲击脉冲宽度 2 ms。
先按冲击谱来：2 ms 脉冲对应的频率量级约 500 Hz，取 500 Hz 最大频率，则时间步长至少要 (\frac{1}{20 \times 500} = 0.1\text{ ms})。再按结构响应来：如果关心到 300 Hz，对应周期 3.33 ms，则步长需要小于 (\frac{3.33}{20}=0.17\text{ ms})。两个里面取更小者，也就是 0.1 ms。
这个案例不需要波传播约束，因为 5 mm 网格对应的时间步长约束是 (\frac{0.005}{5000} = 1\times 10^{-6}\text{ s})，也就是 0.001 ms，比 0.1 ms 小了 100 倍。这就是两种分析模式的本质区别——如果你真的在关心应力波局部传播，那么所有网格单元都成了传感器，计算规模会急剧上升，这种问题通常已经不是 Workbench 隐式瞬态的舒适区了，而应该考虑显式动力学工具。
总时长如果按 20 ms 算，0.1 ms 的步长对应 200 个子步，这对瞬态分析来说是一个相当轻松的计算量。所以时域分析并不一定都贵，贵不贵取决于你对时间分辨率的要求。
4.3 载荷步、子步、自动时间步的工作关系
Workbench 里 Step 是一段固定的求解时间区间，每个 Step 内可以进一步划分成多个子步，每个子步就是求解器真正执行数值积分的一个时间跳跃。多个 Step 的典型用途是模拟“先预紧、再冲击、后观测”的过程：
- Step 1：0~1 s，施加螺栓预紧力，步长可以放得比较大；
- Step 2：1~1.1 s，施加冲击载荷，这是整个过程的核心，步长要按上文方法取小；
- Step 3：1.1~5 s，载荷撤除，观察自由振动衰减，步长可以适当放大。
每一步在 Analysis Settings 里单独设置 End Time 和 Time Step。打开自动时间步（Auto Time Stepping）后，还要设置 Initial Time Step、Minimum Time Step、Maximum Time Step 三个值。我的经验是 Initial 给目标步长的十分之一，让求解器起步时容易收敛；Maximum 按目标步长填，防止它自作主张跳太大；Minimum 作为兜底，通常设为目标步长的千分之一即可。
“Program Controlled”大多数情况下能跑，但它不了解你的物理问题，给出的步长一般偏保守或者偏大。知道目标步长时，强烈建议手动填，这一小步能让你的求解效率提升不少。
5. 我踩过的时域分析坑：异常关闭、不收敛与算不动
5.1 几何结构编辑器异常关闭：一次系统的排查顺序
很多人在 Workbench 里建模建得好好的，突然 DesignModeler 闪退了。遇到这个问题先别急着重装，我的排查顺序是这样：
第一步，重启 Workbench 再双击 Geometry 单元格，如果是偶发性的进程崩溃，这一步就恢复了。第二步，检查项目路径。Workbench 老版本对中文路径和特殊字符的兼容性不太好，把整个工程拷到D:\CAE_Projects\drop_test这类纯英文路径下再打开，这个操作解决过不少“疑难杂症”。
如果还不行，就考虑显卡兼容问题。在 DM 的 Tools、Options、Graphics 里关掉硬件加速，再重新打开几何，市面上很多工作站的显卡驱动版本跟 ANSYS 不兼容，关闭硬件加速后稳定很多。再不行就换建模入口：用 SpaceClaim 或者外部 CAD 把模型另存成 x_t、step 格式，再关联到瞬态分析的 Geometry 里，没必要死磕 DesignModeler。
最后才考虑重装。为了不打断分析流程，我一般会在工程里保留一份备用的外部模型，DM 挂了就直接换路径重新关联，项目接着跑，损失也就是十几分钟。
5.2 求解不收敛的正确排查链路
瞬态分析不收敛的错误弹窗信息往往很笼统，按链路排查比瞎试高效得多。
先查约束。如果模型存在刚体位移，求解器会疯狂迭代然后报错。判断方法很简单：看变形动画是不是整体在平移或旋转，或者看求解信息里有没有弱弹簧警告。约束不足就补约束，该加固定支撑的加支撑，该用远端位移约束的用远端位移。
再查接触。接触初始穿透、接触刚度不足、主从面设置不当都会导致不收敛。求解前先检查 Contact 里的 Gap/Penetration 值，穿透量过大的话调整几何或者接触偏移。
然后看子步。错误信息里经常提示增加子步数量，这时候把自动时间步打开，或者把步长调小一点。很多时候不是模型错了，而是时间步大到求解器根本跨不过接触突变。
最后检查大变形开关。如果结构发生了明显的大转动、大挠度，默认关闭的大变形效应会带来问题，需要手动开启。
我的总原则是：先从最简单的模型开始验证。接触和非线性全开之前，边界条件和载荷本身至少要能在一个线性瞬态分析里跑通，然后再一级一级往上加复杂度，这样就算出问题，也能定位在刚刚加的某一层上。
5.3 计算时间爆炸与结果文件庞大怎么办
瞬态分析的结果文件经常大得吓人，因为每子步都存一份完整结果。Workbench 的 Output Controls 里可以把瞬态结果设置成等间隔保存，或者只保存固定数量的时间点。这样动画可能会跳帧，但该看的峰值和趋势一个不少，文件体积却能缩到原来的几分之一。
计算时间超出预期时，先看网格。动态分析里网格数量和时间步长双重制约，减少过渡区网格、把非关注区的网格放大到合理限度，往往比单纯加 CPU 更有效。再检查接触对数量，不必要的接触对用绑定替代摩擦，能省下大量非线性迭代时间。然后看时间范围，别把不需要的尾段算进去，及时截断分析时长。最后才是启用多核并行，Workbench 在多核下对大型瞬态模型的加速是很明显的，但并行核数增加后收益会递减，不必盲目追高。
还有一个很实用的做法：如果结构满足平面应力或者轴对称条件，优先降成二维模型做参数扫描，结果成熟后再回到三维做最终确认。我做过的一些薄板结构跌落分析，二维模型算一版只要十几分钟，三维模型要跑一晚上，效率和趋势判断的便利性差距非常大。
我后来养成了一个习惯，凡是做瞬态分析，先顺手把模型的模态分析做了。模态频率直接告诉你时间步长该按哪个基准选，振型告诉你关注点应该放在哪个位置，还能提前判断激励频谱会不会撞上共振峰。这个步骤花不了几分钟，但对后续参数设置帮助太大了。
另外一个小经验是给结构加一个很低的瑞利阻尼，比如阻尼比 0.5%~2%，否则纯无阻尼系统算出来的响应曲线全是尖锐的高频毛刺，乍一看还以为是计算错误。阻尼参数可以在整个项目里统一保存成自定义材料库，后面几个类似项目直接调用，省得反复输。时域分析这条路，第一次走通是门槛，走通之后你会发现它的思路其实比静力分析更接近真实世界。
