软件
产品
你跑完模态分析，得到一堆频率和振型数据，怎么快速整理成报告？手动在Visualization界面里翻找节点位移？太慢了。2026年，用Abaqus自带的Python库（Abapy）写几行脚本，就能自动提取结果，输出格式整洁的文本文件。下面我把完整代码拆开讲，附带一个真实工程案例。
模态分析后，你需要知道每阶频率（Hz）和对应的振型向量（节点位移）。在Abaqus/CAE可视化模块里，可以逐个查看振型动画，但想批量导出节点数据就麻烦了——你得创建XY Data、保存、整理格式……一个10阶模态的模型，手动折腾半小时。
一个真实场景：某汽车排气管支架模态分析，模型有2万节点，需要提取支架安装点处的振型参与系数。手动复制粘贴了3个小时，还漏了一阶。后来用脚本自动处理，3秒跑完，准确率100%。
下面这段Python脚本，可以在Abaqus自带的CAE界面里运行（File→Run Script），或者用abaqus python命令在命令行执行。
from abaqus import *from abaqusConstants import *from odbAccess import *odb = openOdb(path='myODB.odb')abaqusConstants提供了枚举值（如SCALAR），odbAccess用来读写结果文件。path替换成你自己的ODB文件路径。
stepName = odb.steps.keys()[-1] # 取最后一个分析步（通常模态分析在最末）step = odb.steps[stepName]frame = step.frames[-1] # 最后一帧包含所有模态frequency = frame.frequency # 提取各阶频率（单位Hz）注意：frame.frequency是一个数组，存放所有提取的模态频率，顺序与Eigenvalue输出顺序一致。
nodeSet = odb.rootAssembly.instances['myInst'].nodeSets['myNodeSet']eigenResult = frame.fieldOutputs['EIGEN']eigenResults = eigenResult.getSubset(region=nodeSet)eigenVectors = eigenResults.valuesfilename = 'modeShapes.txt'with open(filename, 'w') as file: file.write('Modal frequencies (Hz):\n') for i, freq in enumerate(frequency): file.write('{:>3d}: {:>10.4f}\n'.format(i+1, freq)) file.write('\n') file.write('Modal vectors:\n') for i, vector in enumerate(eigenVectors): file.write('Vector {:>2d}:\n'.format(i+1)) for j, displacement in enumerate(vector.data): file.write('{:>8.4f}'.format(displacement)) if (j+1) % 6 == 0: file.write('\n') file.write('\n')每个vector.data包含该节点在X、Y、Z方向的位移分量（实部、虚部或幅值，取决于求解器设置）。这里每行输出6个数据后换行，让文件更易读。
输出示例：
Modal frequencies (Hz): 1: 45.2332 2: 127.8956 3: 301.4471Modal vectors:Vector 1: 0.0123 0.0456 0.0789 0.0001 0.0002 0.0003 0.0112 0.0441 0.0765 0.0001 0.0002 0.0003...from abaqus import *from abaqusConstants import *from odbAccess import *# 1. 打开ODB文件odb = openOdb(path='your_model.odb')# 2. 获取最后一阶模态分析步的最后一帧stepName = odb.steps.keys()[-1]step = odb.steps[stepName]frame = step.frames[-1]# 3. 提取频率freqs = frame.frequency# 4. 指定节点集（根据你的模型修改）myAssembly = odb.rootAssemblyinstanceName = 'PART-1-1' # 实例名nodeSetName = 'SET-1' # 节点集名nodeSet = myAssembly.instances[instanceName].nodeSets[nodeSetName]# 5. 提取振型eigenField = frame.fieldOutputs['EIGEN']eigenSubset = eigenField.getSubset(region=nodeSet)vectors = eigenSubset.values# 6. 输出到文本outputFile = 'modal_results.txt'with open(outputFile, 'w') as f: f.write('=== Modal Analysis Results ===\n') f.write('Frequencies (Hz):\n') for i, fval in enumerate(freqs): f.write(f'{i+1:3d}: {fval:10.4f}\n') f.write('\nNode Displacements (mode shapes):\n') for vecIdx, vec in enumerate(vectors): f.write(f'\nMode {vecIdx+1}:\n') data = vec.data for dIdx, val in enumerate(data): f.write(f'{val:12.6f}') if (dIdx+1) % 6 == 0: f.write('\n') if len(data) % 6 != 0: f.write('\n')print(f'Results saved to {outputFile}')问题1：找不到实例或节点集 报错KeyError。检查实例名和节点集名是否与CAE中一致。可以在CAE命令行里运行mdb.models['Model-1'].rootAssembly.instances.keys()查看所有实例名。
问题2：frame.frequency为空 说明你读的不是模态分析步，或者模态提取没有成功。检查ODB里是否存在EIGEN场输出。可以在Visualization模块里，Create Field Output→EIGEN，看是否有数据。
问题3：振型数据过多导致文本文件巨大 如果节点集包含上千个节点，每个节点有6个分量（3方向实部+3方向虚部），文件会非常大。建议只提取关键节点（如响应点、约束点）的振型。或者在脚本里只写入最大位移节点的数据。
一个工程案例 某风力机叶片模态分析，叶片有15万节点。全提取的话文本文件超过200MB。改为提取叶尖和叶根处两个节点集的振型（共12个节点），文件只有8KB，但足以用于后续的强迫响应计算。
你可以把数据写入CSV，方便导入Excel或MATLAB：
import csvwith open('modes.csv', 'w', newline='') as csvfile: writer = csv.writer(csvfile) writer.writerow(['Mode', 'Freq_Hz'] + [f'Dof_{i}' for i in range(6)]) for i, vec in enumerate(vectors): row = [i+1, freqs[i]] + list(vec.data) writer.writerow(row)然后用Excel画振型曲线，或者用MATLAB的plot函数可视化。
用Abaqus的Python接口（Abapy）自动后处理模态分析，三步走：打开ODB、定位节点集和特征向量、写入文本。相比手动导出，脚本可以把数小时的工作压缩到几秒，而且杜绝人为遗漏。2026年的仿真工程师，掌握这种自动化技能，会让你从重复劳动中解放出来。把上面的代码保存成.py文件，下次模态分析完直接运行，振型和频率一键到手。
武汉格发信息技术有限公司，格发许可优化管理系统可以帮你评估贵公司软件许可的真实需求，再低成本合规性管理软件许可,帮助贵司提高软件投资回报率，为软件采购、使用提供科学决策依据。支持的软件有: CAD,CAE,PDM,PLM,Catia,Ugnx, AutoCAD, Pro/E, Solidworks 等。
技术文档
推荐好文
155-2731-8020
