Abaqus后处理输出任意参数平均值，自定义结果提取

0.出发点

数值仿真数据的合理处理对于结果非常重要，特别是细观模型，包括特殊结构、3D打印材料、颗粒、纤维填充等代表性单元，需要通过提取其平均值才能与试验进行对比。下面是个简单的算例，在y方向压缩极小曲面之Gyroid。

为了对比该极小曲面的应力水平，采用同样的材料做了单轴压缩，两种情况对比如图所示：

从图中可以看到，如果仅去极小曲面上的一个点作为其应力应变，其曲线甚至比实心立方体还高，显然这是不合理的。出现这种现象的主要原因是，此类细观结构或变形不均匀时，不能取一个点代表整个模型，不然很可能会出现较大的误差。对此问题，细观力学有些方法，比如作用力反作用力法，体积平均法等，但也有人认为对于细观结构，作用力反作用力法 不太合理。体积平均法的简单表达式如下：

即模型中每一个单元的应力（应变）对单元体积积分后，除以模型整个体积。上述应力应变曲线也证实，采用该方法能够得到较为真实的数据。

如果进提取某个方向的平均应力-应变，也许采用作用力反作用关系更快些；然而，有时候需要不同参量，如von mises应力等的平均值，此时需要通过python进行批量处理。

****基于以上需求，做了各abaqus后处理插件，可以提取任意参量的平均值。******

介绍

插件主界面如下图所示，功能主要包括：

• 定义任意step；
• 可以选择每隔n个帧提取，实际上这是为了加快提取速度；
• 可以提取标量：当所输出为标量时，比如用户自定义SDV1等，members可以不填写；同时支持多个变量，比如SDV1,SDV2等；
• 有份量的参量：比如应力S及真实应变LE等都是有方向的，此时不仅需要给出参量名称，也需要给定方向，比如要提取S11 则在Name填写S，members填写11; 如果相提取mises平均，则在Name先填写S，members 再填mises即可。同样支持多个量同时提取，比如 Name:S,E 而members：22同时提取S和E的2方向。
• 

• 可以自定义单元集合ElementSet或直接从屏幕选择：二选一；建议先定义好ElementSet，因为直接选取的话也会消耗时间（单元数量很多时）；如果前期没有定义set，后续将odb以可读写方式打开，如图所示，再定义set即可：
• 
• 可以对结果转换坐标输出：支持前处理或后处理自定义坐标后转换输出，如图所示。
• 
• dVi==1：当前处理没有输出EVOL或只想提取算数平均数时勾选，此时提取速度快；如果需要提取体积平均如式(1)，则需勾掉，并必须有EVOL。可以动态显示：
• 

下载地址：https://jmbwy.lanzouo.com/iZsHt1ymdlhe

密码:msib

免责声明：本文系网络转载或改编，未找到原创作者，版权归原作者所有。如涉及版权，请联系删