BoLei space

08 Visualization：后处理、取值口径与结果验证

Properties2
description读取 ODB，区分位移分量与模量、积分点与节点应力，并用反力和理论解验证结果。
tagsabaqus, 基础笔记, 后处理

8 min read

08 Visualization：后处理、取值口径与结果验证

后处理既是出图，也是检查模型是否回答了原问题。开始前先确认 目标步骤已成功完成，再明确要提取的区域、变量、分量、帧和输出位置。

本页使用 已完成的悬臂梁案例，不涉及新的求解。

1. 打开正确的结果

  1. File → Open 选择对应 .odb，核对路径、作业名和实例。
  2. 选择正确的 Step 与 Frame；本例为 Load 最后一帧，步时间 1.0。
  3. 先显示未变形和变形形状，检查整体运动与支承位置。
  4. 再打开 Contours on Deformed Shape，选择变量和分量。
  5. 指定变形放大倍数、视角和图例范围，再导出图片。

ODB 可能只包含失败前的部分帧；也可能属于同目录中的旧作业。帧时间必须与当前分析目标对应。仅看图中“有颜色”不能证明满载结果存在。

2. 位移分量与模量

变量含义适合回答
U1、U2、U3指定坐标方向的位移分量向哪个方向移动多少
U, Magnitude位移向量模长，非负位移总量最大在哪里
UR有相应自由度的节点转角梁、壳、RP 等对象的转动

本例沿 −Y 加载，因此端部下挠查看 U2。位移模量不能保留向上/向下符号。若使用局部坐标变换，应明确图表中的方向约定。

悬臂梁 U2 云图，单位 mm，变形放大 30 倍

本图变形放大 30 倍，图例数值仍是实际计算结果。报告截图时应注明单位、变量、帧和放大倍数；跨方案比较尽量统一视角与色标范围。

3. Mises、分量和主应力

  • S, Mises：等效应力，常用于金属屈服相关的结果讨论，不显示拉压符号。
  • S11 等分量：与指定坐标系有关，可用于检查梁的轴向弯曲应力。
  • 主应力：用于观察主方向上的拉压状态，不能与 Mises 混为同一指标。

查看壳或梁时，还要选择截面点或层位置；“最大应力”如果不说明位置与分量，信息是不完整的。

悬臂梁 Mises 云图，单位 MPa，变形放大 30 倍

4. 积分点、外推和平均

本例实体应力在积分点计算。为了形成连续云图，CAE 会按相应显示设置将应力外推到节点，并可能对相邻单元值做平均。外推和平均后的峰值不必等于原始积分点峰值。官方后处理示例

取值方式本例峰值
原始积分点最大 Mises56.4186 MPa
当前节点外推与平均云图图例约 62.38 MPa

图中 Avg: 75% 表示平均阈值设置，不表示将应力乘以 75%。取消平均也不一定就是原始积分点值：还要确认输出位置是否仍为外推的单元节点值。

因此比较两个结果前，统一变量、坐标系、位置、截面点、外推及平均方式。对于塑性问题，节点外推应力可能高于材料当前屈服应力，不能仅凭这一显示现象断言本构积分错误。

5. 如何提取端部位移和反力

本例使用装配节点集 TIPNODES 和 FIXEDNODES。以所有 36 个自由端节点的 U2 算术平均值作为端部下挠指标；固定端 RF2 则求和而不是平均。

在 Report Field Output 中选定变量、步、帧、区域和输出位置，可导出文本。对固定格式的任务，脚本更容易保持取值一致：

# Abaqus 2022 / Python 2.7；odb 已打开
# 显式使用内置函数，避免 from abaqus import * 后的同名覆盖
import __builtin__ as bi
frame = odb.steps['Load'].frames[-1]
tip = frame.fieldOutputs['U'].getSubset(
    region=odb.rootAssembly.nodeSets['TIPNODES']).values
fixed = frame.fieldOutputs['RF'].getSubset(
    region=odb.rootAssembly.nodeSets['FIXEDNODES']).values
if not tip or not fixed:
    raise RuntimeError('Required node set is empty')
mean_u2 = bi.sum(v.data[1] for v in tip) / len(tip)
reaction_y = bi.sum(v.data[1] for v in fixed)
max_mises = bi.max(v.mises for v in frame.fieldOutputs['S'].values)

该 max_mises 片段针对本例 ODB 中保存的积分点 S。其他 ODB 应先检查实际位置。新版本若使用 Python 3，内置函数模块名改为 builtins。

节点平均是这里明确选择的指标，并非面积加权平均。非均匀网格或其他物理问题应重新确定适当的平均方法。

6. 结果必须做哪些核对

核对本例证据能说明什么
变形方向自由端 U2 为负与 −Y 载荷方向一致
力平衡固定端 RF2 约 +100 N与 −100 N 外力平衡
刚度基准平均下挠 0.1902596 mm与梁理论接近
应力趋势根部外侧较大、沿梁长降低与弯矩趋势一致

理论截面惯性矩及端部挠度：

平均下挠相对理论偏差约 −0.1137%。最大位移模量为 0.1907993 mm，它与平均 U2 是不同指标。

本次没有完成网格收敛和强度评定。位移吻合不保证固定端局部应力已收敛；没有屈服数据和评价准则，也不能直接给出安全系数。网格检查方法见 06 Mesh。

7. 发布一份可复核的结果

至少保留：模型假设与单位、边界和载荷、网格与单元、求解完成证据、取值区域与口径、数值表、带说明的云图、理论或实验对比及局限。

用小型 JSON/CSV 保存关键结果，正文表格从数据生成；图片使用笔记库内相对路径；大型 CAE、ODB 单独归档。后处理脚本应检查目标帧和集合是否存在，缺失时停止，不要静默生成空表。

完整数据、复现资料与工作流讨论：abaqus-cantilever-static。返回系列入口：01 Part。

Graph View

Backlinks