得到结果

S,U,RF,LE这四个value的成员名称竟然完全一模一样，里面竟然都有mises(应该只有应力才有mises应力的说法)。这实在是说不通，唯一说得通的可能是，这些value的成员都是应力。但经过实验发现他们各自不相等，不可能都代表应力。
但根据官网这个说法，那么实际上确实都是应力，只不过有的不合理

进行个数统计的结果：

得到

就个人感觉，官方文档似乎也没有很详细的把abaqus script每个涉及到的内置的结构体讲的很清楚。
这里还有一本书

是讲解这方面知识的。

那么，提取位移，应变等正确的办法是什么呢？

v.data。就可以提取出想要的3个自由度的位移了
应变出了一点问题:
field Output request里面明明有E，但是最后

# E_output = o1.steps['Step-1'].frames[-1].fieldOutputs['E'] # 报错——根本没有E，即使在fieldoutput中增加E request，最后字段里也没有,但是好消息是，有EE。我大概明白这个逻辑了，就是在abaqus里面E选上之后，最后就会同时有PE，EE.但似乎也不太对。。。

原来在这里
http://130.149.89.49:2080/v2016/books/ker/default.htm?startat=book01.html#ker
有详尽的介绍

FieldValue object
http://130.149.89.49:2080/v2016/books/ker/default.htm?startat=book01.html#ker

根据这个说法，那么实际上确实都是应力。

abaqus里面的 conjugate data 是指
the imaginary portion of a complex result

最终答案:
通过查文档以及
https://www.cnpython.com/qa/171659
得到

S11_output = o1.steps['Step-1'].frames[-1].fieldOutputs['S'].getScalarField(componentLabel='S11') # 这一行仍然报错。

经过这三行的实验，成功确定了 elementfieldvalues.data
第一个元素（共四个元素）是S11，
第四个元素是S12,
第二个元素是S22，
第三个元素是S33。
S13,S23两个直接就不存在

print S_output.getSubset(region=ElementSet[i]).values # 这一行返回2个 <type 'FieldValue'> array

我推测为什么fieldValueArray有2个元素，是因为一个是全局的，一个是局部的。经过实验确认第一个是全局的。
通过自己搭模型也可以确认很多事情——例如values[0]是全局的，values[1]不知道是什么东西，等等.
现在明确了 values的成员是不同积分点的输出。
另外，abaqus很特殊的一点是：其位移输出是全局坐标系下的量，但是对复合材料，**它的应力和应变输出都是主方向坐标轴的量。**对于普通材料，师兄给我演示的例子显示似乎输出的是全局坐标系的量

哈哈，这个不是和式分解吗？（abaqus theory guide）

S4R单元在我本次实验3方向的应力和应变似乎都被强制设置为0.

在撤去BC-3（z=0），以及增加横向载荷（magnitude=10）后，静力分析的结果不再收敛。调整magnitude=1后，结果收敛，而且εz\varepsilon_zεz即ε33\varepsilon_{33}ε33仍然为0.如下图
通过下图这段话，结合上面的应力、应变均为0的实验结果，可以做这样一个判断：abaqus认为我们使用的S4R单元是thin-only shell elements.

http://130.149.89.49:2080/v2016/books/usb/default.htm?startat=pt06ch29s06alm19.html#usb-elm-eusingshellsection
这是这个页面有用的地方——如何定义泊松比：
OdbFile可以用来生成模型
odb的显示问题
原来odb的信息还可以直接在窗口查看（如下,查看了所有的fieldoutput，果然没有E）
控制台的输出：

hahaha
({'baseElementTypes': ('S4R',), 'bulkDataBlocks': ' object', 'componentLabels': (), 'description': 'Active yield flag', 'isComplex': OFF, 'locations': 'FieldLocationArray object', 'name': 'AC YIELD', 'type': SCALAR, 'validInvariants': (), 'values': 'FieldValueArray object'})
({'baseElementTypes': (), 'bulkDataBlocks': ' object', 'componentLabels': ('CF1', 'CF2', 'CF3'), 'description': 'Point loads', 'isComplex': OFF, 'locations': 'FieldLocationArray object', 'name': 'CF', 'type': VECTOR, 'validInvariants': (MAGNITUDE,), 'values': 'FieldValueArray object'})
({'baseElementTypes': (), 'bulkDataBlocks': ' object', 'componentLabels': ('CM1', 'CM2', 'CM3'), 'description': 'Point moments', 'isComplex': OFF, 'locations': 'FieldLocationArray object', 'name': 'CM', 'type': VECTOR, 'validInvariants': (MAGNITUDE,), 'values': 'FieldValueArray object'})
({'baseElementTypes': ('S4R',), 'bulkDataBlocks': ' object', 'componentLabels': ('LE11', 'LE22', 'LE33', 'LE12'), 'description': 'Logarithmic strain components', 'isComplex': OFF, 'locations': 'FieldLocationArray object', 'name': 'LE', 'type': TENSOR_3D_PLANAR, 'validInvariants': (MAX_INPLANE_PRINCIPAL, MIN_INPLANE_PRINCIPAL, OUTOFPLANE_PRINCIPAL, MAX_PRINCIPAL, MID_PRINCIPAL, MIN_PRINCIPAL), 'values': 'FieldValueArray object'})
({'baseElementTypes': ('S4R',), 'bulkDataBlocks': ' object', 'componentLabels': ('PE11', 'PE22', 'PE33', 'PE12'), 'description': 'Plastic strain components', 'isComplex': OFF, 'locations': 'FieldLocationArray object', 'name': 'PE', 'type': TENSOR_3D_PLANAR, 'validInvariants': (MAX_INPLANE_PRINCIPAL, MIN_INPLANE_PRINCIPAL, OUTOFPLANE_PRINCIPAL, MAX_PRINCIPAL, MID_PRINCIPAL, MIN_PRINCIPAL), 'values': 'FieldValueArray object'})
({'baseElementTypes': ('S4R',), 'bulkDataBlocks': ' object', 'componentLabels': (), 'description': 'Equivalent plastic strain', 'isComplex': OFF, 'locations': 'FieldLocationArray object', 'name': 'PEEQ', 'type': SCALAR, 'validInvariants': (), 'values': 'FieldValueArray object'})
({'baseElementTypes': ('S4R',), 'bulkDataBlocks': ' object', 'componentLabels': (), 'description': 'Magnitude of plastic strain', 'isComplex': OFF, 'locations': 'FieldLocationArray object', 'name': 'PEMAG', 'type': SCALAR, 'validInvariants': (), 'values': 'FieldValueArray object'})
({'baseElementTypes': (), 'bulkDataBlocks': ' object', 'componentLabels': ('RF1', 'RF2', 'RF3'), 'description': 'Reaction force', 'isComplex': OFF, 'locations': 'FieldLocationArray object', 'name': 'RF', 'type': VECTOR, 'validInvariants': (MAGNITUDE,), 'values': 'FieldValueArray object'})
({'baseElementTypes': (), 'bulkDataBlocks': ' object', 'componentLabels': ('RM1', 'RM2', 'RM3'), 'description': 'Reaction moment', 'isComplex': OFF, 'locations': 'FieldLocationArray object', 'name': 'RM', 'type': VECTOR, 'validInvariants': (MAGNITUDE,), 'values': 'FieldValueArray object'})
({'baseElementTypes': ('S4R',), 'bulkDataBlocks': ' object', 'componentLabels': ('S11', 'S22', 'S33', 'S12'), 'description': 'Stress components', 'isComplex': OFF, 'locations': 'FieldLocationArray object', 'name': 'S', 'type': TENSOR_3D_PLANAR, 'validInvariants': (MISES, MAX_INPLANE_PRINCIPAL, MIN_INPLANE_PRINCIPAL, OUTOFPLANE_PRINCIPAL, MAX_PRINCIPAL, MID_PRINCIPAL, MIN_PRINCIPAL, TRESCA, PRESS, INV3), 'values': 'FieldValueArray object'})
({'baseElementTypes': (), 'bulkDataBlocks': ' object', 'componentLabels': ('U1', 'U2', 'U3'), 'description': 'Spatial displacement', 'isComplex': OFF, 'locations': 'FieldLocationArray object', 'name': 'U', 'type': VECTOR, 'validInvariants': (MAGNITUDE,), 'values': 'FieldValueArray object'})
({'baseElementTypes': (), 'bulkDataBlocks': ' object', 'componentLabels': ('UR1', 'UR2', 'UR3'), 'description': 'Rotational displacement', 'isComplex': OFF, 'locations': 'FieldLocationArray object', 'name': 'UR', 'type': VECTOR, 'validInvariants': (MAGNITUDE,), 'values': 'FieldValueArray object'})
['baseElementTypes', 'bulkDataBlocks', 'componentLabels', 'description', 'isComplex', 'locations', 'name', 'type', 'validInvariants', 'values']
<type 'FieldValue'>
['baseElementType', 'conjugateData', 'conjugateDataDouble', 'data', 'dataDouble', 'elementLabel', 'face', 'instance', 'integrationPoint', 'inv3', 'localCoordSystem', 'localCoordSystemDouble', 'magnitude', 'maxInPlanePrincipal', 'maxPrincipal', 'midPrincipal', 'minInPlanePrincipal', 'minPrincipal', 'mises', 'nodeLabel', 'outOfPlanePrincipal', 'position', 'precision', 'press', 'sectionPoint', 'tresca', 'type']
['baseElementType', 'conjugateData', 'conjugateDataDouble', 'data', 'dataDouble', 'elementLabel', 'face', 'instance', 'integrationPoint', 'inv3', 'localCoordSystem', 'localCoordSystemDouble', 'magnitude', 'maxInPlanePrincipal', 'maxPrincipal', 'midPrincipal', 'minInPlanePrincipal', 'minPrincipal', 'mises', 'nodeLabel', 'outOfPlanePrincipal', 'position', 'precision', 'press', 'sectionPoint', 'tresca', 'type']
['baseElementType', 'conjugateData', 'conjugateDataDouble', 'data', 'dataDouble', 'elementLabel', 'face', 'instance', 'integrationPoint', 'inv3', 'localCoordSystem', 'localCoordSystemDouble', 'magnitude', 'maxInPlanePrincipal', 'maxPrincipal', 'midPrincipal', 'minInPlanePrincipal', 'minPrincipal', 'mises', 'nodeLabel', 'outOfPlanePrincipal', 'position', 'precision', 'press', 'sectionPoint', 'tresca', 'type']
['baseElementType', 'conjugateData', 'conjugateDataDouble', 'data', 'dataDouble', 'elementLabel', 'face', 'instance', 'integrationPoint', 'inv3', 'localCoordSystem', 'localCoordSystemDouble', 'magnitude', 'maxInPlanePrincipal', 'maxPrincipal', 'midPrincipal', 'minInPlanePrincipal', 'minPrincipal', 'mises', 'nodeLabel', 'outOfPlanePrincipal', 'position', 'precision', 'press', 'sectionPoint', 'tresca', 'type']
3200
1681
3200
1681
0
({'coordinates': array([254.0, 247.649993896484, 0.0], 'f'), 'instanceName': 'MENGPI_REAL-1', 'label': 1640})

这个图是shell element里面的thickness change，在membrane单元里也有类似的thickness change

20220219:终于被我找到藏在这里——

For plane stress, membrane, and shell elements, only the in-plane tensor components (11, 22, and
12 components) are stored by Abaqus/Standard. The out-of-plane direct component for stress (S33) is
reported as zero to the output database as expected, and the out-of-plane component of strain (E33) is
reported as zero even though it is not. This is because the thickness direction is computed based on
section properties rather than at the material level. The out-of-plane components can be requested for
field-type output and cannot be requested for history-type output. The out-of-plane stress components
are not reported to the data (.dat) file or to the results (.fil) file.

此外,Script User Guide 里搜索S11，或component或试试其他分量名称可以找到一张完整分量的表格。

addData方法：

向odb中增加数据。这个方法是在创建odb时使用的，像我只是读取job默认生成的odb文件（只是读取odb），用不到addData方法。

abaqus script提取应力应变位移 odb学习 addData相关推荐

ABAQUS中的应力应变描述
更多阅读:sppy.site 应力应变描述文档<Abaqus Analysis User's Guide>-><Introduction, Spatial Modeling, ...
提取abaqus中节点集合的应力应变
节点是没有应力应变的,只有积分点才有,所以直接在odb文件里面用FieldOutput是整不出来的(结果是个空列表),只能用session这种让abaqus自己帮你处理.具体的代码需要根据自己操作的时 ...
python提取abaqus数据_Python语言在ABAQUS数据提取中的简单应用
前面某期文章提到过ABAQUS中利用Python脚本进行后处理数据的提取,这个问题有必要做些讨论和总结. 在ABAQUS中,利用Python脚本做简单的处理是非常方便和易懂上手的,这里做简单介绍,参考 ...
利用python提取abaqus节点坐标的脚本_用于在Abaqus中提取结点力的Python程序
用于在Abaqus中提取结点力的Python程序 #coding=utf-8 from abaqus import * from abaqusConstants import * from odbAc ...
关于Abaqus+python提取s11残余应力（按路径提取）
目录写在前面代码过程与方法总结写在前面接上一篇文章,通过Abaqus+python提取残余应力S11的方法!今天又有新的方法,可能这是很常用的方法,但是我现在才发现. 先来看一张效果图: ...
[有限元] Ansys Workbench Mechanical 中的应力应变显示类型的文档翻译
应力和应变给定零件或整个装配体的模型和材料,以及特定结构载荷环境的情况,应力解决方案使您能够(结构的)预测安全系数.应力.应变和位移. 一般的三维应力状态,是根据与零件或装配体世界坐标系对齐的三个法 ...
应力应变基础理论分析
文章目录 1.基本假设与基本物理量 2.应变与位移的关系 3.应力与应变关系 5.变形连续方程(相容方程) 6.位移分量求解应力及边界条件在分析激光加工过程中的热应力问题之前,首先讨论一下弹性体的应 ...
comsol结构力学-应力应变仿真
今天聊聊comsol结构力学仿真部分,主要是后处理部分. 前面的几何.材料都比较简单,按照自己的需求绘制对应的图,选择相应的材料即可.物理场部分主要是施加约束条件,包括固定约束.边界载荷.指定位移等, ...
PCB应力应变测试分析结合IPC-9702和IPC-9704A标准
一,PCB应力应变测试背景应力应变测试的必要性,在现代社会中,电子设备充斥着我们的生活,包括汽车.计算机.智能手机.飞机等等.在这些产品中,都集成了印刷电路板(PCB),作为电子产品的核心部件,印刷 ...

abaqus script提取应力应变位移 odb学习 addData

addData方法：

abaqus script提取应力应变位移 odb学习 addData相关推荐

最新文章

热门文章