时间步长 — CAE术语解释
时间步长的理论基础
时间步长的物理意义
「时间步长」只是计算的细化划分吗？它在物理上有什么意义呢？
不仅仅是计算上的划分。它是正确捕捉「物理现象传播速度」最重要的离散化参数。例如在冲击波传播问题中，时间步长
CFL条件如何具体确定？当网格单元尺寸差异很大时，以哪一个为基准？
最严格（最小的）条件控制整体时间步长。公式为：
如此小的时间步长，要仿真1秒的现象需要几十万步，计算永远都不会结束。实际工程中如何平衡？
这正是工程实践的关键。主要有两个对策。首先，网格划分时避免物理上不必要的「微小单元」。例如铆钉的圆角，如果对现象无影响，应予以简化。其次，利用求解器的「质量缩放」功能。在Ansys LS-DYNA等中，可以人为提高最小单元的刚性，扩大允许的时间步长，从而缩短计算时间。但必须评估引入的人工质量对结果的影响。
时间步长的数值计算方法
显式方法与隐式方法中的处理差异
听说「显式方法」和「隐式方法」中时间步长的重要性完全不同。具体有什么区别？
本质区别在于稳定性。显式方法（如中心差分法）是「条件稳定」的，严格需要满足CFL条件等时间步长上限。隐式方法（如Newmark法、后向欧拉法）号称「无条件稳定」，理论上可以使用任意大的时间步长而不会发散。但隐式方法在步长过大时会产生过度数值阻尼，「压制」物理振荡，引发另一类问题。
隐式方法时间步长过大时振荡被压制，这是什么机制？能用公式说明吗？
以后向欧拉法为例，离散化后的「振幅衰减率」取决于时间步长。设单振动的圆频率为
那么隐式方法应该如何设置时间步长？「不发散就行」显然不够。
完全同意。隐式方法需要根据「待解问题的最高频率成分」来设置。经验法则：解振荡现象时，至少用20～30步划分一个周期。即
时间步长的实际应用
合适的时间步长确定工作流程
实际分析中，如何初始化时间步长，又如何逐步调整？请给出具体步骤。
实务流程如下：1) 网格划分后，用前处理器（如LS-PrePost或HyperMesh）确认最小单元尺寸；2) 查询材料的音速。例如钢材纵波速度约5900 m/s；3) 用CFL条件估算
经常看到「自动时间步长控制」功能。这个能信任吗？原理是什么？
可以信任，但要理解其原理。Ansys Mechanical的「自动时间步」监视前一步的收敛速度（迭代次数），收敛快则增大下一步步长，收敛慢则减小。COMSOL的「自动步长」由用户设定「允许局部截断误差」，求解器动态调整步长使其落在误差范围内。但在初接触或破裂等急剧变化时，自动控制可能跟不上，导致发散。最好的做法是在这些关键时段手动锁定较小的步长。
接触问题需要特别细化时间步长，为什么？
接触发生时，边界条件瞬间改变，产生应力波。这种波含有极高频率成分，粗步长无法捕捉。另外，显式方法中接触算法（如惩罚法）自身引入数值「弹簧」，其固有振动频率可能要求更短的步长。以汽车门闩啮合分析为例，接触前后的1毫秒常需用平常的1/10甚至更短的步长计算。
时间步长的软件比较
主要求解器中的时间步长设置差异
Ansys LS-DYNA和Abaqus/Explicit都是显式方法代表，时间步长设置和默认行为有区别吗？
差别很大。LS-DYNA默认进行某种程度的「质量缩放」，试图达到用户指定的「目标时间步长」。通过关键字`*CONTROL_TIMESTEP`中的`DT2MS`参数，可以锁定所需的步长，求解器会相应调整最小单元的质量。而Abaqus/Explicit默认不进行质量缩放，计算稳定限界步长，再乘以安全系数0.9作为实际步长。用户若要锁定步长，需指定`*FIXED TIME INCREMENTATION`。这体现了两者的哲学差异：LS-DYNA重视「固定计算时间」，Abaqus重视「物理精度优先」。
隐式方法求解器呢？Abaqus/Standard和Ansys Mechanical的时间增量控制比较如何？
Abaqus/Standard的自动时间增量控制很成熟。在`*STEP`中设置`INC`（最大增分数）和`INITIAL`（初始增分），求解器根据收敛情况自动在`MIN`和`MAX`之间调整步长。Ansys Mechanical也有类似功能，但参数通常以「弱非线性」「强非线性」等预设形式提供，不如Abaqus直观可控。对于包含复杂接触的非线性动力学分析，Abaqus的控制更加鲁棒。
开源求解器（如CalculiX、Code_Aster）相比商业软件如何？
基本概念相同，但自动化程度和用户界面差异大。CalculiX (CCX)的动态显式分析中，关键字`*DYNAMIC`需手动设置`ALPHA`、`BETA`（阻尼参数）和`DELTMX`（最大增分），基本无自动控制。Code_Aster更成熟，可定义`PAS_APRIORI`（先验步长）和`PAS_MINI`（最小步长），设置`ADAPTATION='OUI'`启用误差自适应。但缺乏图形化引导，需深入理解理论和手册，学习曲线陡峭。
时间步长的故障排除
由时间步长引起的典型错误和对策
显式法分析出现「Negative Volume」或「Hourglass」错误，是不是时间步长太大？
「Negative Volume」典型症状就是步长过大，一步内单元剧烈变形、翻转（雅可比行列式为负）。先尝试将步长减半。「Hourglass」成因不同，常见于降积分单元。这类单元存在「零能量模式」（不产生应变能的变形），步长过大时此模式失控。对策是缩小步长，同时加强Hourglass控制（如LS-DYNA中提高`IHQ`参数）。
隐式非线性静力分析报「不收敛」，缩小荷载增分就一定能解决吗？
多数情况有效，但不绝对。减小增分使解变化更接近线性，牛顿-拉夫逊法易收敛。但在接触或材料不连续性极强的问题中，可能存在「硬极限」，再小的增分也不收敛。此时需1) 调整接触刚性；2) 切换到动态隐式法（Abaqus/Standard的`*DYNAMIC`）加入阻尼；3) 启用高级收敛辅助选项如「放松」或「线搜索」。Ansys中打开`SOLCONTROL=ON`可自动化这些设置。
计算全程完成，但结果应力和变形充满明显「噪声」，物理上不合理。这也是时间步长问题？
是的，特别在动力分析中。显式法步长过粗会无法正确传播高频物理模式，噪声由此产生。隐式法则可能过度数值阻尼压制了低频，导致结果「钝化」。对策：首先细化时间步长重新计算，检查噪声是否减少；如仍不行，可在后处理阶段对结果应用低通滤波。例如，汽车碰撞评估标准SAE J211/ISO 6487规定使用通道滤波等级（CFC）处理数据，这本质上是为清除数值噪声。
相关主题
更详细
错误