Abaqus自动网格划分插件：Python二次开发实战

简介：当前CAE前处理中，复杂几何模型的网格划分往往依赖手动反复操作，工作量大且结果受经验影响。这套基于Python开发的ABAQUS网格自动划分插件，正是面向这类批量前处理与自动化建模需求，适用于具备Python基础的仿真工程师和希望扩展ABAQUS二次开发能力的学习者。资源包共4个文件，内含2个Python脚本（插件注册主程序与网格工具模块）与2个GIF演示动图，压缩包仅73KB，轻量小巧便于快速查看效果。动图展示了切割网格与轮廓提取等关键功能，脚本则覆盖了从加载几何模型、设定网格尺寸到生成并检查网格的完整逻辑；借助脚本内封装的函数，可针对不同几何特征生成六面体或四面体网格，并执行质量检查与局部细化。已有3313人学习下载，说明该方案在同类需求中具有一定参考价值。通过研读meshTools_plugin.py及配套代码，可以了解插件在ABAQUS中的集成方式与常用二次开发思路，同时，清晰封装的接口也便于直接复用或改造到实际项目中，有助于缩短建模时间、提高网格划分的规范性与可重复性。

1. Abaqus 网格自动划分插件：把三小时手工分网压到八分钟

某次结构强度评审前，一块带孔安装板的网格一直过不了关。手工分区、调种子、重画，前后试了三个多小时，网格还是歪的。后来我把分区规则写成一个基于 Python 的 Abaqus 插件，从导入几何到网格完成只需要八分钟，雅可比稳定大于 0.7。这个插件做的事不复杂：脚本识别孔洞、台阶和薄壁特征，自动补充分区面，再按特征密度布种。适合做强度分析和反复改几何的工程师，新手能直接跑通，熟手可以改分区和播种参数来适配自家模型。

2. 运行机制：分区策略、几何识别与 Python API 三件套

2.1 插件不是一个大脚本，而是三个文件的协作

写过 Abaqus 脚本的人都知道，一个分网流程往往就是一个几百行的 .py 文件，用abaqus cae noGUI=xxx.py从头跑到尾。但做成交互插件后，事情就不一样了。一个能住进 Plug-ins 菜单的自动分网插件，通常由三块拼出来：内核脚本（kernel）、GUI 对话框（RSG 生成）、注册文件（plugin）。内核脚本负责所有几何操作和网格指令，GUI 负责收集参数，注册文件告诉 Abaqus 菜单里放哪个按钮。

文件	职责	改动频率
my_auto_mesh.py	识别特征、分区、播种、生成网格	每次调算法都改这里
my_auto_mesh_gui.py	RSG 生成的参数面板	很少改
my_auto_mesh_plugin.py	菜单注册、参数传递	增加参数时才动

my_auto_mesh/ ├── my_auto_mesh.py # 内核脚本：识别、分区、布种、生成 ├── my_auto_mesh_plugin.py # 注册文件：菜单入口与参数传递 ├── my_auto_mesh_gui.py # GUI 对话框：参数面板 └── my_auto_mesh_plugin.png # 菜单图标，可选

我一般会把内核脚本写成「去掉 GUI 也能独立跑」的纯逻辑。原因很实际：交互调参的时候用 GUI，批量改模型的时候用 noGUI 脚本，两者共用同一个内核函数，维护成本最低。如果直接把参数面板的逻辑写进内核脚本，后面做批量计算时会非常痛苦，因为每次都得去点一遍对话框。

2.2 自动分区的核心问题：特征先于网格

手工分网时，人脑会在几秒内完成一个判断：哪里是孔、哪里是倒角、哪里是薄壁，然后决定分区边界和种子密度。脚本要把这个「人脑判断」编码成规则，也是整个插件最有争议的部分——处理得太糙，分区失败；处理得太细，碰上奇怪几何就崩。

常见的做法是遍历 edges 和 faces，做两层特征识别：

• 圆形边（edge.geometryType == CIRCLE）视为孔特征，取其圆心、半径和轴向。
• 半径突变或夹角突变处视为台阶/圆角边界，后续播种时单独加密。
• 壳体部件没有 cells，识别到part_obj.cells为空时，要走 face 流程，不能按实体硬切。

这样做能显著控制单元数量。一个直径 100mm 的板，全局种子 2mm，孔边种子 0.5mm，孔周围的单元密度是远区的 16 倍，但只在孔附近产生，整体单元数反而比均匀加密少了 30% 以上。先把特征找出来，再决定怎么分区，这是自动分网的底层逻辑。

3. 核心代码解读：从几何识别到播种完成的三个函数

3.1 几何识别函数 find_holes

自动化第一步是让脚本「看见」孔。Abaqus 的边对象自带几何类型属性，圆孔的边是CIRCLE类型，可以通过edge.getRadius()拿到半径、圆心和法向。注意不是所有圆边都是真的孔——倒角、圆角、圆柱凸台也会有圆形边，所以要用半径下限过滤掉微特征。

# -*- coding: utf-8 -*- from abaqusConstants import CIRCLE def find_holes(part_obj, radius_min=0.1): """识别部件里的圆孔特征，返回 [(center, radius, axis), ...]""" holes = [] for edge in part_obj.edges: if edge.geometryType == CIRCLE: try: radius, center, normal = edge.getRadius() except AttributeError: # 旧版本接口可能不含 getRadius，跳过而不是中断 continue if radius >= radius_min: holes.append((center, radius, normal)) return holes

这里radius_min的取值很关键。默认 0.1mm，会把螺栓孔、定位孔都识别出来，但也会把加工圆角算进去。如果模型里有大量 R0.5 以下的圆角，我一般会把它调到 1.0，用孔边种子处理真正需要关心的孔。getRadius()返回的顺序是 (radius, center, normal)，不同小版本可能有差异，保险做法是先打印出来确认一下。

3.2 自动分区函数 partition_by_holes

识别到孔之后，下一步是把孔周围的几何切成独立块。这里用的思路是：过孔心、沿孔轴线方向做一个平面切分。切分后孔变成「独立块」的一部分，局部加密时可以直接针对块区域操作，而不是靠边选择一个个去碰运气。

def partition_by_holes(part_obj, holes, axis_map=(0.0, 0.0, 1.0)): cur_cells = [c for c in part_obj.cells] for center, radius, normal in holes: # 过孔心、沿孔轴向做平面切分 part_obj.PartitionCellByPlaneOriginNormal( origin=center, normal=axis_map, cells=cur_cells) cur_cells = [c for c in part_obj.cells] # 注意：每次切分后部件特征对象会刷新，必须重新取 cells return cur_cells

代码里的坑在于cur_cells的刷新。Abaqus 每次执行分区后，原来拿到的 cell 对象引用就会失效，如果继续用旧列表，第二次切分时大概率报cells have been deleted。axis_map是孔轴方向，薄板模型通常 Z 向垂直板面，直接传(0, 0, 1)即可；斜孔多的模型，最好用hole里返回的normal代替，避免切分平面和孔轴不垂直导致切偏。

3.3 播种与网格生成函数 auto_seed_and_mesh

分区完成后进入播种环节。Abaqus 的播种机制是全局种子加局部覆盖：先对整个部件布一个较粗的全局种子，再用seedEdgeBySize对孔附近的边做局部加密。这比直接全部件细播效率高，也能保证孔边捕捉几何。

def auto_seed_and_mesh(part_obj, holes, global_size=2.0, hole_size=0.5, deviation=0.05): # 第一步：全局粗种子 part_obj.seedPartBySize(size=global_size, deviationFactor=deviation, minSizeFactor=0.1) # 第二步：对每个孔周围的边做局部加密 for center, radius, normal in holes: bbox_min = (center[0] - radius * 2, center[1] - radius * 2, center[2] - radius * 2) bbox_max = (center[0] + radius * 2, center[1] + radius * 2, center[2] + radius * 2) edge_list = part_obj.edges.getByBoundingBox( bbox_min[0], bbox_min[1], bbox_min[2], bbox_max[0], bbox_max[1], bbox_max[2]) if len(edge_list): part_obj.seedEdgeBySize(edges=edge_list, size=hole_size, deviationFactor=deviation) # 第三步：统一用四面体兜底，确保自动流程不死在单元类型上 elem_type = ElemType(elemCode=C3D10, elemset=(part_obj.cells,)) part_obj.setMesh(regions=part_obj.cells, elemTypes=(elem_type,)) part_obj.generateMesh(regions=part_obj.cells)

deviationFactor控制网格对曲线的贴合程度，默认 0.05 表示允许边偏离几何 5% 的种子尺寸。minSizeFactor=0.1是限制最小单元尺寸为全局种子的 10%，防止小边、微缝把单元尺寸拖到极小。最后用 C3D10 是故意设的兜底——自动分区的结果不保证每个 cell 都能扫掠成六面体，硬上 C3D8R 会导致大片失败，退回四面体反而稳定。如果你的模型拓扑简单且确认可扫掠，把elemCode改成 C3D8R 即可。

4. 安装部署与参数调整：从 GUI 调用到命令行传参

4.1 插件安装路径与菜单注册

插件装错了路径是新手最常见的翻车点。Abaqus 的插件搜索目录有两个：一个是 CAE 安装目录下的site-packages/abaqus_plugins，全局生效；另一个是用户主目录下的~/abaqus_plugins，只对当前系统用户生效。我建议用用户目录，避免权限问题和污染安装目录。

cp -r my_auto_mesh ~/abaqus_plugins/

重启 CAE 后，Plug-ins 菜单下就会出现入口。如果没出现，先确认插件目录里有没有__init__.py，没有的话 Abaqus 不会把它识别成插件。另一个常见问题是注册文件名和内核文件名不一致，注册文件里写死的模块名一旦对不上，菜单项会直接消失，而且不报错。

4.2 参数调整与典型场景预设

GUI 面板上暴露的参数不需要太多，五个就够：全局种子尺寸、孔边种子尺寸、偏差系数、最小尺寸系数、单元类型。参数之间是强耦合的，调整时要看模型的量级而不是盲目抄数值。

场景	全局种子	孔边种子	偏差系数	说明
静力学分析	5~8mm	1~2mm	0.05	单元数控制在 10 万以内
疲劳/裂纹扩展	2~4mm	0.5mm	0.02	孔边必须加密，偏差越小越好
模态分析	8~12mm	2~3mm	0.1	追求低单元数，高频模态才加密

minSizeFactor一般固定 0.1，调小会让网格在小特征处爆量，调大又可能丢失几何细节。如果模型里有很多 R0.2 以下的微小圆角，建议先做几何清理，而不是靠调这个参数硬扛。

4.3 命令行方式：批量改几何时的正确姿势

交互点菜单适合调试，真正干活时还是命令行靠谱。插件内核脚本设计成了可独立运行，于是批量处理可以用 noGUI 模式，把模型名和参数通过命令行传进去。

abaqus cae noGUI=auto_mesh_batch.py -- plate.cae Part-1 2.0 0.5

对应的脚本入口只需要十几行：

# auto_mesh_batch.py import sys from abaqus import mdb model_name = sys.argv[-4] part_name = sys.argv[-3] global_size = float(sys.argv[-2]) hole_size = float(sys.argv[-1]) p = mdb.models[model_name].parts[part_name] holes = find_holes(p) cells = partition_by_holes(p, holes) auto_seed_and_mesh(p, holes, global_size, hole_size) mdb.save()

注意sys.argv的前几位是 Abaqus 内部参数，从--之后开始取才是自定义参数，所以这里从倒数第四个开始读。批量跑多个模型时，外面套一个 for 循环改模型名即可，全程不需要打开 CAE 界面。调试时建议先跑一个模型，确认分区和播种结果没问题再上批量，否则翻车了只能看日志猜原因。

5. 避坑记录：五个反复踩中的网格划分问题

5.1 特征掩码导致脚本回放失败

现象：手工在 GUI 里点选面、边后录制的宏，直接拿到脚本里跑就报错，或者换个模型就选中了错误对象。

原因：Abaqus 在 GUI 交互时生成的特征选择指令是getSequenceFromMask，这个掩码和当前 CAE 会话的数据库状态强绑定，换机器、换模型、重启 CAE 后掩码就可能失配。

解决：脚本里一律改用findAt(coordinates)或getByBoundingBox(...)选取对象。用坐标选对象虽然慢一点，但可重复执行，不依赖界面状态。这是自动分网插件能稳定跑起来的先决条件。

5.2 中面模型被当成实体划分

现象：从三维软件导出的钣金件抽中面后导入 Abaqus，插件直接报cells are empty，或者把面壳当成实体去分区，越切越乱。

原因：中面模型、壳模型没有cells，只有faces和edges。插件没做类型判断，一上来就取part_obj.cells，自然拿不到任何有效对象。

解决：在内核脚本最前面加一个分支判断，len(part_obj.cells) == 0时走面壳流程，用PartitionFaceBySketch或直接只播种不分区。壳单元的种子参数和实体也不一样，seedPartBySize在壳模型上作用的是面网格，elemCode要用S3、S4而不是C3D10。

5.3 布尔分区时提示 Feature failed

现象：PartitionCellByPlaneOriginNormal执行到某个孔时报错，提示无法用平面切分 cell，或者切出来的块几何形状怪异。

原因：原始几何存在微间隙、短边、零面积面，切分平面刚好穿过这些瑕疵区域。薄板模型尤其容易踩中，因为厚度方向只有一两毫米，平面稍微偏一点就和面重合了。

解决：先做几何修复。在 Abaqus CAE 里用 Geometry Edit 删掉多余的小面，或者把导入的 iges/step 模型先跑一遍 Repair 工具。插件侧的对策是：切分前检查center是否在 cell 内部，用cell.getSize()或cell.isOutOfPlane一类的接口做预判，不在内部就直接跳过这个孔，不强行分区。

5.4 曲率控制参数在低版本 Abaqus 里失效

现象：播种时勾选了曲率控制，设了曲率因子，但生成网格后孔边还是粗得离谱，仿佛参数没传进去。

原因：seedPartBySize的curvatureControl参数在 Abaqus 2019 之前的版本里默认关闭，而且部分旧版本即使传了参数也不生效。新版改了底层逻辑，但很多公司还在用旧版。

解决：不要依赖曲率控制实现局部加密，改用显式的seedEdgeBySize按边加密。这样无论版本高低，行为都一致。如果确实需要曲率适应，先确认版本支持，再显式传入：

part_obj.seedPartBySize(size=2.0, curvatureControl=(True, 0.5, 10.0))

5.5 单元类型设置导致回退到四面体

现象：插件明明设了六面体单元 C3D8R，生成网格时静默退化成 C3D10，或者直接报错无法生成。

原因：分区后每个 cell 并不都是可扫掠的六面体块。只要有一个 cell 拓扑不满足扫掠条件，整批 setMesh 就会失败或降级。这是自动分网和手工分网最大的差别——手工分网时人会控制分区形状，脚本不会。

解决：一种做法是全部用 C3D10 兜底，稳定但单元数偏多；另一种做法是先尝试C3D8R，生成失败再局部退回过C3D10。实现上可以用try/except包裹generateMesh，失败时把该 cell 单独重设一次单元类型。代价是代码复杂度和调试时间翻倍，我通常只在模型拓扑非常稳定时才启用六面体优先模式。

6. 进阶验证：网格质量批量检查与关键区域加密技巧

插件能出网格只是第一步，网格动不动得了、结果信不信得过，还得靠质量检查说话。Abaqus 自带的 Verify Mesh 面板交互好用，但批量场景下没法逐个点。我习惯在内核脚本里加一个质量统计函数，遍历所有单元，用节点坐标近似计算最长边和最短边的比值，作为纵横比的快速指标。

def quality_report(part_obj): coords = {n.label: n.coordinates for n in part_obj.nodes} aspects = [] for elem in part_obj.elements: labels = [n.label for n in elem.nodes] pts = [coords[l] for l in labels] dists = [] for i in range(len(pts)): for j in range(i + 1, len(pts)): d = sum((pts[i][k] - pts[j][k]) ** 2 for k in range(3)) ** 0.5 dists.append(d) aspects.append(max(dists) / min(dists)) return {'max_aspect': max(aspects), 'avg_aspect': sum(aspects) / len(aspects)}

这个函数是粗筛，判断标准比 Abaqus 内部算法宽松得多，但用来对比插件参数调整前后的趋势足够了。真正要出报告，还是要回到 CAE 里跑一遍 Verify Mesh，看雅可比和翘曲角。关键区域的加密也不止孔边一种：圆角根部、台阶过渡面这些位置，可以先elset或edgeset选出目标区域，再重复调用seedEdgeBySize，优先级天然覆盖全局种子。

我从那以后养成了一个习惯：每改一次分区或播种参数，就强制跑一遍quality_report，把max_aspect控制在 10 以内才往下做计算。网格不是越多越好，也不是越细越好，是每一处密度都花在该花的地方。这套插件用到现在，至少帮我省出了几十个手工分网的晚上，希望也能帮到你。

本文还有配套的精品资源，点击获取