Abaqus Von Mises应力定义与S变量关联

Abaqus中Mises应力详解：从理论到实践

一、Mises应力的基本概念

Abaqus中的Mises应力，又称冯·米赛斯应力，是工程仿真中判断材料是否发生屈服的关键指标。它基于第四强度理论，将复杂的应力状态简化为一个标量值，便于不同位置的应力大小直接对比。

二、Mises应力的计算原理

Mises应力是基于剪切应变能理论计算的等效应力。其数值为正且不随坐标系变化，计算公式为：\[ \sigmav = \sqrt{\frac{(\sigma1 - \sigma2)^2 + (\sigma2 - \sigma3)^2 + (\sigma3 - \sigma_1)^2}{2}} \]\(\sigma1, \sigma2, \sigma_3\)为三个主应力。

三、Mises应力在工程中的应用

Mises应力用于判断延性好的金属材料（如钢材、铝合金等）是否发生屈服。当Mises等效应力达到材料的屈服强度时，材料开始进入塑性变形阶段。

四、在Abaqus中查看Mises应力云图

• 进入可视化模块：计算完成后，将模块切换到"Visualization"，视图区会显示模型轮廓。
• 绘制云图：单击工具箱中的"绘制云图在变形图上"按钮，或执行"绘图→云图→在变形图上"命令，视图区就会显示Mises应力云纹图。
• 切换显示方式：按住工具箱中的云图按钮展开，选择显示在变形前模型上或变形前后对比。

五、提取特定点的Mises应力值

• 使用查询工具：在Visualization模块选择工具栏的"Query"工具，点击"Probe values"，选"Node"或"Element"，在模型视图中点选你关心的节点，Abaqus会直接显示该点的Mises应力值。
• 使用Python脚本批量提取：编写Python脚本打开odb文件，fieldOutputs['S'].values[i].mises读取每个单元的Mises应力，再写入txt文件。

六、绘制应力随时间曲线

• 定义节点集合：将关心的节点定义为集合。
• 定义历史变量：在分析步中定义历史变量输出该节点集合的Mises应力。
• 查看曲线：计算完成后，在历史变量结果中查看曲线。

七、Mises应力值的意义

Mises应力是一个基于剪切应变能理论的标量值，它反映了材料因形状改变而积累的能量。与单轴应力不同，Mises应力在多向应力耦合时更为准确，不能直接用CPress或单轴应力分量替代。

八、常见问题及注意事项

• 管理员身份运行：在Abaqus中处理文件时，以管理员身份运行，以确保权限问题不会影响操作。
• 单元影响：一阶单元（如C3D8R）应力在积分点计算，外推至节点可能不平滑；二阶单元（如C3D10）精度更高但需注意接触面应力震荡问题。
• 自定义实现：若需修改屈服面或特殊硬化规则，需编写UMAT子程序手动实现Von Mises流动法则及一致切线刚度矩阵。