后处理 — CAE用语解释

这是一个很大的误解。后处理的本质是从大量数值数据中提取"有意义的信息"，并将其可视化为"支持决策所需的依据"。例如，不仅要显示应力分布，还要定量表示在降伏应力为245MPa的条件下，最大Mises应力在哪里达到320MPa，该区域的体积占整体的百分之几。这才是支持设计判断的关键。

这因分析领域而异。结构分析看位移、应变、安全系数、接触压力。热分析看温度分布、热流量、传热系数。流体分析（CFD）看速度矢量、压力分布、流线、涡度、湍流能等。以汽车轮毂空气动力分析为例，使用Q准则或λ₂准则等涡识别方法可视化旋转轮毂周围的复杂涡流结构，评估其对空气阻力（Cd值）的影响。

计算公式因材料和失效准则而异。对延性材料来说，最简单的是基于Mises应力的安全系数。设屈服应力为σ_y，等效应力为σ_eq，则安全系数n为

这叫做"插值"或"平滑化"处理。最常见的方法是使用单元内的形状函数进行插值。对于一阶四面体单元（tetra），单元内任意点的物理量φ可用节点值φ_i和形状函数N_i计算为

使用"裁切"或"探针"功能。例如，Ansys的"Section Plane"可计算用户定义的平面（如x=0.1m）与网格单元的交线。交线上的点通过所在单元的形状函数对物理量进行插值。对于沿路径的图表（"曲线"），先将路径离散化为细密的点群，然后对每个点进行相同的插值计算，最后绘制成折线图。

基本上是这样。将位移矢量u乘以缩放系数α来计算显示位移u_disp：

不行。首先要从"模型是否按意图工作"的基本检查开始。具体来说：1) 反力确认：例如悬臂梁在端部施加10N荷载时，固定端的反力总和（矢量和）应该是10N；2) 质量确认：从材料密度和体积计算的FEM总质量是否与CAD数据相差不大（例如CAD是1.5kg但FEM是0.01kg，说明网格有缺陷）。如果跳过这些检查，后续的所有结果都将毫无意义。

大多数软件都有"符号显示"功能。Ansys中，固定支持显示为三角形标记，强制位移显示为箭头，分布压力显示为覆盖表面的均匀箭头群。必须检查这些符号的方向和大小（缩放）。特别是旋转约束和力矩荷载难以看清，所以要通过位移结果（即使变形很小也要放大显示）来确认行为。例如施加力矩后，应该出现扭转或弯曲变形。

用不同网格尺寸（粗、标准、细）计算的关键最大应力值或最大位移值要汇总成表格或折线图。以单元数或节点数为横轴，应力值为纵轴。如果细化网格后结果变化在2%以内，则认为已收敛。未收敛则需在应力集中部进一步细化网格。Abaqus/CAE的"XY Data"功能可制作此图。

至少要遵守这些规则：1) 色条的最大/最小值要设置为可与材料许可值（如屈服应力）或规范基准值比较的范围。例如S45C的屈服应力是365MPa，就设最大值为400MPa。2) 必须标注单位（MPa、mm、°C）。3) 在图中明确标注变形缩放系数（Deformation Scale Factor: 200等）。忽略这些会导致变形程度被误解。可参考JIS Z 8316（技术制图—显示的一般原则）。

理念差异很大。Ansys Mechanical面向"工程师"，提供安全系数、疲劳寿命（与nCode联动）、自动报告生成等与设计决策直接相关的功能。Abaqus/CAE面向"研究人员和开发者"，强项是用户定义场输出（UVARM）和Python脚本的高度定制化后处理。COMSOL的卖点是"多物理场集成"，可在同一画面上叠加显示不同物理场（结构、电磁、流体）的结果，直观评估相互影响。

根本区别是"计算功能的有无"。查看器（Ensight、FieldView、ParaView）纯粹用于可视化，擅长读取已计算的数据并漂亮地显示，但无法计算新的物理量。而有偿软件的后处理模块（Ansys的"Mechanical"或"CFD-Post"）有强大的"派生结果"功能，可从结果数据进一步计算。例如，可从应力张量现场计算主应力或Mises应力，或从时间历程数据进行傅里叶变换求频率响应。

ParaView是开源软件，定制性高，可用Python脚本自动化批处理和复杂可视化（如粒子追踪的复杂过滤），适合研究和新手法开发。STAR-CCM+的查看器针对汽车等行业的标准评估指标（如下压力系数Cl的时间序列图、偏航角下的空力特性曲线）有丰富的模板，工作效率压倒性地高。定常业务报告建议用商用查看器。

首先怀疑"数值奇点"。尖锐角承受集中荷载，或肋的根部等尖锐边缘网格粗疏的情况下，理论上应力趋于无穷大。这种表观高应力网格相关性强，细化网格时值只会不断上升。对策：1) 改为现实的面压荷载；2) 角部加圆角；3) 观察结果时采用"排除尖峰的平均值"或"距离奇点数mm处的值"。Ansys的"Probe"功能可在距奇点数mm处取应力。

很可能是"接触设置"问题。初级错误是没定义接触对，或定义了但"接触算法"不当。Ansys中，本应设为"Bonded"（不分离）的地方误设为"Frictionless"（无摩擦滑动）就会穿透。对策：1) 重新确认接触区，正确设置接触和目标面；2) 用"Adjust to Touch"调整初始接触防止尖角；3) 非线性接触中，用小步长（子步）施加荷载，每步检查接触状态是否正确。

通常不是计算误差，而是"网格粗"造成的。壁面附近速度梯度急峻，网格粗时速度矢量本身就是离散的，流线积分就会产生不自然的轨迹。对策有二。其一，改变后处理设置："流线起点数增加""积分步长细化"。根本解决是其二：正确设置边界层网格。例如用乱流模型时，无量纲壁面距y+设为约1，第一层厚度设为0.01mm左右，层数设为15～20层。

这是显示缩放问题。位移或应力的等高线显示中，色条最大/最小值基于"整个模型"自动设置。大零件变形10mm、小螺栓变形0.1mm时，螺栓几乎全是最小值颜色，看不到细节。对策：1) 仅选中该零件，将"显示对象"改为"选定的体"单独绘制；2) 手动设置色条范围（如0～0.5mm）。Abaqus中还有"移除上下限"的设置。