★ FEATURED ARTICLE 我入行前两年一直觉得结构仿真就是“把力加上去看应力云图”。直到有一次做设备底座的跌落分析静力算出来的最大应力连屈服强度一半都不到按静力标准绰绰有余可样机实测摔了几次就出现裂纹。回头看才意识到这类问题根本不能当静力算——它必须用 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 这个频带取阻尼比 (\xi0.03)那么 (\omega_162.8\text{ rad/s})、(\omega_21256.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 N0.001 到 0.002 s 保持 1000 N0.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 10~1 s施加螺栓预紧力步长可以放得比较大Step 21~1.1 s施加冲击载荷这是整个过程的核心步长要按上文方法取小Step 31.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%否则纯无阻尼系统算出来的响应曲线全是尖锐的高频毛刺乍一看还以为是计算错误。阻尼参数可以在整个项目里统一保存成自定义材料库后面几个类似项目直接调用省得反复输。时域分析这条路第一次走通是门槛走通之后你会发现它的思路其实比静力分析更接近真实世界。
阅读完成 · 觉得有帮助？
