你好！从你的代码和描述来看，核心问题出在对Abaqus应力场数据的理解和处理上，导致zip迭代时出现不匹配，同时也无法正确提取Mises应力值。下面我会详细分析问题原因，并给出修正后的完整代码。

问题根源分析 #

1. 错误的应力数据提取方式
   你当前直接调用的odb.steps['Belastung'].frames[-1].fieldOutputs['S']返回的是全应力张量场（每个节点对应一组应力分量，比如3D问题是6个分量：S₁₁、S₂₂、S₃₃、S₁₂、S₁₃、S₂₃），而非Mises应力的标量值。这意味着block.data的每个元素是一个数组，而非单个数值，直接用zip(block.nodeLabels, block.data)会导致迭代时的元素长度不匹配，进而引发错误。

2. 无效的数值判断逻辑
   你代码中的if mises_value > 300.0是在和一个数组做比较，这在Python中会引发ValueError（无法将数组与标量比较），这也是你感知到“zip不支持迭代”的间接原因。

修正后的完整代码 #

import csv
from odbAccess import openOdb
from abaqusConstants import *

#Extract Mises

def getMISES(output_path_mises, odb_path):
    # 建议添加readOnly=True避免意外修改ODB
    odb = openOdb(odb_path, readOnly=True)
    
    # 关键修改1：从应力张量场中提取Mises应力的标量场
    step_frame = odb.steps['Belastung'].frames[-1]
    stress_field = step_frame.fieldOutputs['S']
    mises_data = stress_field.getScalarField(invariant=MISES)
    
    # 目标装配体实例
    target_instance = odb.rootAssembly.instances['PRUEFKOERPER_4-1']
    node_coords = {node.label: node.coordinates for node in target_instance.nodes}
    # 预存目标实例节点标签集合，用于快速校验
    target_node_set = set(node_coords.keys())

    with open(output_path_mises, 'w') as csv_file:
        csv_writer = csv.writer(csv_file)
        csv_writer.writerow(['NodeLabel', 'x', 'y', 'z', 'MISES'])

        for block in mises_data.bulkDataBlocks:
            # 关键修改2：仅处理属于目标实例的应力数据块，提升效率
            if block.instance == target_instance:
                # 此时block.data每个元素是单个Mises数值，可与节点标签正确配对
                for node_label, mises_value in zip(block.nodeLabels, block.data):
                    # 关键修改3：浮点数比较逻辑正常生效
                    if mises_value > 300.0:
                        coords = node_coords.get(node_label, [None, None, None])
                        csv_writer.writerow([
                            node_label, 
                            coords[0], coords[1], coords[2], 
                            '{:.4f}'.format(float(mises_value))
                        ])

    odb.close()

关键修改点说明 #

1. 正确提取Mises标量场
   通过stress_field.getScalarField(invariant=MISES)从全应力张量中生成Mises应力的标量输出，此时block.data的每个元素都是单个浮点数，和block.nodeLabels的迭代项完全匹配，解决了zip迭代的核心问题。

2. 过滤目标实例的数据集
   增加if block.instance == target_instance:判断，只处理属于目标装配实例的应力数据块，避免遍历装配体中其他无关实例的数据，大幅提升脚本运行效率。

3. 修复数值过滤逻辑
   现在mises_value是单个浮点数，mises_value > 300.0的过滤条件可以正常工作，不会再出现类型不匹配的错误。

额外注意事项 #

• 你的Abaqus 2019完全支持getScalarField(invariant=MISES)方法，无需版本顾虑。
• 如果目标实例存在部分节点未被输出应力数据的情况，node_coords.get(node_label)会返回默认的[None, None, None]，你可以根据需求调整默认值或增加异常处理。

备注：内容来源于stack exchange，提问作者Jawo