Abaqus用python进行.odb文件提取

目录

写在前面

Abaqus：

什么是有限元？

 .odb文件结构

Odb对象的成员信息如下：

rootAssembly对象的成员信息如下：

 使用：

what？S11？

方案分析：

 以下程序实现S11值的读取：

最近跟着老师做课题，过程中需要用到python对abaqus的.odb等各种文件进行读取，其他文件好说，Abaqus的.odb文件具有庞大数据结构，不能直接读取。幸运的是Abaqus和python之间有联动的方法。这篇文章主要记录了我在接触Abaqus和python联动的时遇到的问题及一些体会。

ABAQUS 是一套功能强大的工程模拟的有限元软件，其解决问题的范围从相对简单的线性分析到许多复杂的非线性问题。 ABAQUS 包括一个丰富的、可模拟任意几何形状的单元库。并拥有各种类型的材料模型库，可以模拟典型工程材料的性能，其中包括金属、橡胶、高分子材料、复合材料、钢筋混凝土、可压缩超弹性泡沫材料以及土壤和岩石等地质材料，作为通用的模拟工具， ABAQUS 除了能解决大量结构（应力 / 位移）问题，还可以模拟其他工程领域的许多问题，例如热传导、质量扩散、热电耦合分析、声学分析、岩土力学分析（流体渗透 / 应力耦合分析）及压电介质分析。

什么是有限元？https://baike.baidu.com/item/%E6%9C%89%E9%99%90%E5%85%83/1761759?fromModule=lemma_inlink

原计划的是提取每一个节点的坐标值，将他和对应的节点的S11值结合起来，但是结果无终。后来，学长给了另外两种思路：

1、创建一条路径，直接把这条路径上的S11值全部读取出来，再进行处理。

2、在前处理中在软件中得到需要读值的网格编号，在后处理中用前面得到的编号去对应S11的值。这样就可以得到网格的S11值？

疑问？S11的值不是应该分布在节点上吗？为什么通过网格编号能获取到S11的值呢？

答：获取的是积分点的S11值，当然也可以是其他点，或者全部节点。

在ABAQUS中对应力的部分理解 关于abaqus中 mises， s11 s22 s33 ，s12，tresca pressure， max principal， mid principal，min principal。简单地理解， 在ABAQUS中，一般是把X轴当做1轴，Y轴当做2轴，Z轴当做3轴；那么： S11就是X轴向的应力，正值为拉应力，负值为压应力; S22就是Y轴向的应力，正值为拉应力，负值为压应力; S33就是Z轴向的应力，正值为拉应力，负值为压应力; S12就是在YZ平面上，沿Y向的剪力; S13就是在YZ平面上，沿Z向的剪力; S23就是在XZ平面上，沿Z向的剪力; 由于剪力的对称性：S12=S21，S13=S31，S23=S32 Mises应力是即第四强度理论，根据能量守恒原理，用于判断材料是否屈服的应 力准则，即Mises准则，一般使用于判断延性比较好的材料，对于脆性材料，一 般采用第一强度理论。

原计划：提取每一个节点的坐标值，将他和对应的节点的S11值结合起来。

提取的结果是COORD的label和S11的label对应不上（如下）后来又查了很多资料，最后发现坐标完全可以从.inp文件中读取，.inp文件中还有C3D8类型，下面这篇文章就说了，C3D8和Nodes的关系，我想可以通过这个能够做到他们和应力坐标对应。具体的我就没再往下进行了，有兴趣的读者可以自己试一下！

手把手教你读abaqus的inp文件（以一个单元的inp文件为例） - 知乎手把手教你读abaqus的inp文件（以一个单元的inp文件为例）首先写出全部代码，一个单元的inp文件完整流程图（可先大致浏览后跳过，接下来分步骤进行详细说明）： *Heading ** Job name: 3pct2Mesh Model name: Mode…https://zhuanlan.zhihu.com/p/180455168

 方案1：创建一条路径，直接把这条路径上的S11值全部读取出来，再进行处理。

abaqus中好像可以通过path命令进行一条路径上的读取。具体的也没多接触。只因为我觉得方案2会更简单，更方便去实现！

方案2：在前处理中在软件中得到需要读值的网格编号，在后处理中用前面得到的编号去对应S11的值。

 在abaqus前处理中我们可以使用一下方法进行网格编号的查询：