简介：这套Abaqus脚本压缩包面向需要频繁进行球体建模与装配的科研人员、工程师及高校学生，解决批量创建不同半径球体并精确控制装配位置的重复操作问题。脚本调用Abaqus的Python建模接口，使用者仅需修改少量参数即可生成任意半径的球体，并可按指定位置完成多球体装配；支持批量操作，能一次性产出多个模型，有效提升复杂结构力学分析和多球体装配模拟的前处理效率。压缩包共1个文件，为单个py脚本，大小仅1KB，轻量易用，适合下载后直接放到工作目录中调用或二次开发扩展。目前已有307人学习下载。借助脚本，读者可获得一套简洁的参数化建模思路，免去手动逐个建模与装配的繁琐流程，同时可作为Abaqus脚本化建模的入门范例，帮助理解Python与Abaqus交互建模的基本方法。
1. 把「任意半径球体 + 任意位置装配」写成 Abaqus 脚本：一次省一下午的思路
做细观模型或者颗粒类前处理时，经常要在 Abaqus 里生成几十个半径大小不一的球体，再按精确坐标摆到指定位置。如果靠 GUI 一个个创建 Part、切到 Assembly 再平移，半小时起步，坐标稍多一点还容易点错行。用 Abaqus 脚本建模，球体生成、半径赋值、装配坐标全变成数据，跑一遍脚本几秒钟收工，参数留在代码里随时可改。这个标题说的就是这件小事：把球体建模和装配从「手工操作」变成「参数驱动的脚本」，适合做骨料模型、滚珠简化、颗粒增强模拟或者批量球体前处理的工程师——哪怕你只写过几行 Python，照着把半径列表和坐标列表替换掉就能用。
2. 生成任意半径球体的最小脚本：从一步能跑通到真正参数化
2.1 一个能直接运行的球体建模骨架
Abaqus 的 Python 脚本入口是mdb这个模型数据库对象。做球体建模时，常规流程是：先建一个模型（Model），再在模型里建 Part，Part 用BaseSolidSphere特征直接生成实心球体。下面是去掉装配逻辑后的最简版本，直接在 CAE 的 File → Run Script 里跑就能在模型树里看到球体 Part。
# -*- coding: utf-8 -*- from abaqus import * from abaqusConstants import * # 模型名统一管理，方便后面删除旧模型重跑 model_name = 'SphereDemo' if model_name in mdb.models.keys(): del mdb.models[model_name] model = mdb.Model(name=model_name) # 创建 Part：三维、可变形实体 part_name = 'Ball_5' part = model.Part(name=part_name, dimensionality=THREE_D, type=DEFORMABLE_BODY) # 关键就这一行：生成半径 5.0 的实心球体 part.BaseSolidSphere(radius=5.0)这段代码里有三个点值得说明。del mdb.models[model_name]是为了让脚本可以反复运行而不报「模型已存在」错误，实际项目中脚本经常要调参数跑几十遍，没有这行会翻车。dimensionality=THREE_D表示创建三维部件，type=DEFORMABLE_BODY表示变形体，这是球体最常用的类型。BaseSolidSphere(radius=5.0)对应 GUI 里 Create Part 对话框的 Sphere 选项，它生成的球体基准点默认在球心，这一点后面装配时非常关键。
这个脚本跑完后，Abaqus 会创建一个名为Ball_5的可变形实体 Part，几何上就是一个完整球。如果只想验证建模，到这里就已经完成了标题里的前半件事。
2.2 半径参数化：把写死的 5.0 变成半径列表
「任意半径」不是让你每次改代码里的数字，而是把半径放到一个列表里，用循环批量生成。例如需要半径 5、8、3、10 四个球，常见做法是维护一个半径列表，循环建 Part 并打印结果。
# 半径列表：想生成多少种半径就加多少行 radii = [5.0, 8.0, 3.0, 10.0] for r in radii: part_name = 'Sphere_r%g' % r # 同一个模型里 Part 名不能重复 if part_name in model.parts.keys(): del model.parts[part_name] part = model.Part(name=part_name, dimensionality=THREE_D, type=DEFORMABLE_BODY) part.BaseSolidSphere(radius=r) pdb = model.parts[part_name] print('Part %s 创建完成，半径 %.3f' % (part_name, r))Sphere_r%g是一个字符串模板，%g会把浮点数转成不带多余零的字符串，比如 5.0 转成 5。这样每个半径对应一个独立 Part，名字一眼能看出大小。这里有个选择：如果多个球半径完全相同，可以只建一个 Part，装配时创建多个 Instance 指向同一个 Part，模型文件更小；如果每个球半径都不同，就每个半径建一个 Part。标题强调「任意半径」，所以上面代码按每个半径一个 Part 处理，简单直接。
这里要注意 Python 字符串格式化在 Abaqus 老版本（基于 Python 2）和新版本（基于 Python 3）里都兼容，%操作符两边都支持。不要把print(...)写成print ...这种 Python 2 老语法，新版会直接语法报错。
2.3 脚本的三种执行方式：GUI 运行、后台运行与批处理
写好的脚本怎么跑，直接影响使用习惯。日常调试推荐打开 CAE，用 File → Run Script 选择脚本文件，这样能立刻在视口里看到生成的 Part，出错时错误信息也完整。参数稳定后要批量跑，用命令行后台执行更省事：abaqus cae noGUI=make_spheres.py。这种模式下 Abaqus 不弹界面，跑完直接退出，适合晚上挂机批量生成模型。
第三种是配合abaqus cae script=make_spheres.py的方式，它会启动 CAE 但自动执行脚本，执行完不关闭界面，适合需要接着手工检查的场景。我一般按需切换：调试用 Run Script，正式生成用 noGUI，要检查装配位置就用 script 模式。无论哪种方式，脚本里print的内容都会输出到 Abaqus 的消息窗口或命令行终端，这就是后面排查问题的主要依据。
还有一点：脚本里不要依赖 GUI 才有的对象。比如caeModules模块里的某些交互对象，在 noGUI 模式不可用。纯建模装配的脚本，from abaqus import *和from abaqusConstants import *就够用了，不要多引入界面相关模块，否则换执行方式就报错。
3. 把球体放到任意位置：装配坐标的三种写法和批量装配脚本
3.1 装配坐标系与球心基准：位置语义比想象中重要
球体 Part 生成后，它的局部坐标系原点就在球心。装配到 Assembly 时，Abaqus 会把 Part 的局部坐标系对齐到装配的全局坐标系。也就是说，如果不做任何移动，球心就在全局坐标原点。要放到任意位置，本质就是把 Part 的基准点（球心）移动到目标坐标。
这里有一个最常见的认知偏差：有人以为 Part 的基准点在球体底部，于是目标坐标写的是球体底部位置，结果球体有一半埋进接触面。用脚本装配前，先记牢一句话：球体 Part 的基准点就是球心，平移向量就是球心的目标坐标。如果你希望球体底面落在某个平面，球心坐标要往上抬一个半径值。
Abaqus 里 Instance 的坐标位置有三种指定方式。第一种是创建 Instance 时直接传position参数；第二种是创建后调用translate平移；第三种是配合rotate调整姿态。球体是各向同性的几何体，旋转对最终模型没有影响，所以「任意位置」只需关心平移，这让装配代码比一般机械装配简单不少。
3.2 创建 Instance 并用 translate 移动到目标坐标
最常见的装配写法是先创建 Instance，再调translate移动。注意translate的vector是增量位移，不是绝对坐标。因为球体基准在球心且初始在原点，第一次平移的增量就等于目标坐标。
# 延续上一节的 model 对象，创建装配实例 asm = model.rootAssembly if not asm.datums: asm.DatumCsysByDefault(CARTESIAN) # 用已有的球体 Part 创建 Instance inst_name = 'Ball_5_Inst' asm.Instance(name=inst_name, part=model.parts['Sphere_r5'], dependent=ON) # 把球心移动到 (10.0, 20.0, 15.0) asm.translate(instanceList=(inst_name, ), vector=(10.0, 20.0, 15.0))model.rootAssembly是模型的根装配，所有部件都要挂到它下面。DatumCsysByDefault(CARTESIAN)是建立一个默认笛卡尔坐标系，没有它大多数情况下也能创建 Instance，但显式建立更稳妥。dependent=ON表示 Instance 依赖 Part 的几何，不复制一份几何数据，节省内存；后续如果对 Part 重新划分网格，Instance 会跟着更新。
translate的第三个参数是相对位移，这个机制很关键。脚本里为了不出错，我通常让每个 Instance 只创建一次、只移动一次，移动的目标值就是最终坐标。不要在同一个 Instance 上反复调translate去「修正」位置——那是增量，不是定位，调第二次等于在上一次基础上叠加位移，很容易跑偏。
3.3 批量装配：半径列表与坐标列表一一对应的完整脚本
把第 2 章和第 3 章的内容合并，就是一个完整可用的批量球体装配脚本。半径列表和坐标列表一一对应，循环里建 Part、建 Instance、移动到目标坐标，最后打印结果。
# -*- coding: utf-8 -*- from abaqus import * from abaqusConstants import * # ====== 用户只需要改这两个列表 ====== radii = [5.0, 8.0, 3.0, 10.0, 6.5] centers = [ (0.0, 0.0, 5.0), (15.0, 0.0, 8.0), (-6.0, 9.0, 3.0), (0.0, -12.0, 10.0), (3.0, 4.0, 6.5), ] # ================================== model_name = 'SphereModel' if model_name in mdb.models.keys(): del mdb.models[model_name] model = mdb.Model(name=model_name) asm = model.rootAssembly if not asm.datums: asm.DatumCsysByDefault(CARTESIAN) # 必须保证半径数量和坐标数量一致 if len(radii) != len(centers): raise ValueError('radii 与 centers 长度不一致') for i, (r, (x, y, z)) in enumerate(zip(radii, centers)): # 每个半径建一个 Part，名字带半径便于识别 part_name = 'Sphere_r%g' % r if part_name not in model.parts.keys(): part = model.Part(name=part_name, dimensionality=THREE_D, type=DEFORMABLE_BODY) part.BaseSolidSphere(radius=r) # 创建 Instance 并移动到目标球心坐标 inst_name = 'Ball_%02d' % i asm.Instance(name=inst_name, part=model.parts[part_name], dependent=ON) asm.translate(instanceList=(inst_name, ), vector=(x, y, z)) print('%s: r=%.3f 球心=(%.3f, %.3f, %.3f)' % ( inst_name, r, x, y, z))这段代码的逻辑：先检查半径列表和坐标列表长度是否一致，不一致直接抛异常，避免最后坐标错位。循环里，相同半径的 Part 只创建一次，后面同半径的球复用同一个 Part，既能多个球共享几何，又不至于每个球都重复建 Part。zip(radii, centers)把两个列表按位置配对，第 i 个球对应第 i 组坐标。
这里有一个实用细节：如果某个半径在后续模型里会成为单独的部件并单独划分网格，那就该为它建独立 Part；如果只是当作几何体放进装配不做精细网格，复用 Part 完全没问题。实际项目中，我一般先想清楚这些球后续要不要独立设置网格属性，再决定复用策略。
3.4 「球心坐标」与「落座坐标」：不同场景下的坐标换算
第 3.1 节留了一个尾巴：实际建模时，你手里的坐标未必是球心坐标。常见有两种来源：一种是从 CAD 或表格里拿到的球心位置，直接填入centers即可；另一种是要求球体「落座」在某个平面上，比如底面放在 z=0 的平台上，此时你手上的是接触点坐标，换算公式是：球心坐标 = 接触点坐标 + 径向偏移。
具体来说，球体底面要落在 z=0 平面，球心 z 值等于半径值。上面的示例里半径 5 的球，坐标 (0,0,5) 就是底面落座在 z=0 的写法；如果直接用 (0,0,0)，球就会有一半埋在平面下方。如果是球体侧面贴合 y=0 平面，球心 y 坐标等于半径，其他方向同理。
我建议在建列表之前先写一行注释标明坐标语义，比如# centers 全部为球心坐标或者# centers 全部为球体最低点坐标，代码里自动加半径。这个习惯能少踩很多坑，尤其是坐标一旦上百个，人很难逐个心算。
4. Abaqus 脚本球体装配的关键参数设置：半径、坐标、单位与浮点精度
4.1 半径和坐标的输入方式：脚本内列表、命令行参数、外部文件
参数从哪进来，决定了脚本的通用性。最简单的是脚本内的常量列表，打开脚本改数字就行，适合参数不频繁变的场景。参数要频繁试算时，从外部文件读入更舒服，既不碰脚本，又能用 Excel 编辑。
脚本内定义参数我一般放在文件顶部，用分隔线标注「用户只需要改这里」。这样做的好处是后续维护的人不用在几百行代码里找参数。外部文件读入的常见做法是读 CSV，每行一个球的半径和坐标，脚本循环读取：
import csv params = [] with open('spheres.csv', 'r') as f: reader = csv.reader(f) for row in reader: if not row or row[0].startswith('#'): continue r, x, y, z = [float(v) for v in row[:4]] params.append((r, x, y, z))关于命令行传参，Abaqus 的cae noGUI=script.py模式不同版本对sys.argv的透传行为不完全一致，我不建议把关键参数压在命令行上。真要传，可以走环境变量或临时文件，但都不如 CSV 直观。实际做项目时，我的习惯是「脚本只做执行，数据全在外部文件」，脚本里不出现任何硬编码的半径或坐标值，这样更换数据集不用动代码。
4.2 单位制的坑要从源头堵住：全模型必须统一长度单位
Abaqus 本身是无单位制的，你写 5.0，它不关心是 5 毫米、5 厘米还是 5 米，全模型数值一致即可。问题在于，球体半径来自设计图纸，装配坐标来自另一份表格，两份数据用了不同单位，装配出来就是一场灾难。
最常见的是毫米和米的混用：半径按毫米写，坐标按米写，或者反过来，球体会小到看不见或者大到盖住整个模型。解决思路是在脚本开头定义一个单位注释，统一写明这套模型采用的长度单位，并确保半径列表与坐标列表来自同一单位体系。我通常会加一行：# 单位制：mm / N / t / s / MPa，对应的导出应力单位是 MPa，质量单位是吨。如果你习惯用 kg，那就是 mm / N / kg / s / kPa？这个换算表容易记错，稳妥做法是在脚本里固定写LENGTH_UNIT = 'mm'，写入日志。
4.3 浮点精度与装配容差：容易被忽略的两个细节
浮点精度在球体装配脚本里翻车概率不高，但一旦翻车很难查。常见问题有两个：坐标用整数运算导致的小数截断，以及用==判断浮点数是否相等。前者出现在坐标由计算得出而非手动输入的场景，比如循环里x = i * 0.1，累积多次后可能变成 0.30000000000000004；后者出现在检查球体是否重叠时，用==比较两个坐标，永远不相等，判断失效。
处理方式很直接：坐标计算统一用浮点数，不要在中间环节强转 int；重叠判断不要比相等，比距离小于半径和。还有一点，装配完成后如果要检查球体是否贴合某个平面，不要依赖坐标完全相等，用差值绝对值小于一个很小的容差值判断，比如abs(z - target_z) < 1e-6。这个容差大小取决于模型尺寸，微米级模型用 1e-9 都不嫌小，米级模型 1e-4 也正常。
第 4 章讲的这些参数设置，本质上都是在减少「脚本能跑但结果不对」的概率。脚本能跑只是第一步，跑出来的模型对不对，更多取决于数据层面的一致性和精度控制。
5. 脚本建模球体装配避坑指南：我踩过的 5 个坑与排查套路
5.1 球体消失或大到盖住全模型：半径单位与直径/半径混用
现象：脚本跑完，视口里要么找不到球，要么球体巨大，整个模型只看到一个球的局部。
原因：一般是单位混用或半径/直径语义混用。有人从图纸上拿到的是直径 50，直接当半径写进BaseSolidSphere(radius=50)，球的实际尺寸比预期大一倍；或者半径按毫米设计，装配坐标却用了米，球在视口里缩成像素点。
解决：第一步，确认图纸给的是半径还是直径，BaseSolidSphere(radius=...)参数只认半径。第二步，统一单位，把脚本顶部单位注释当成硬规范，半径列表和坐标列表来自同一单位体系。排查时最快的方法是打印球心坐标和半径，用肉眼算一下球体边缘范围，看是否落在装配区域内部。
5.2 球心坐标明明对，球却埋进接触面一半
现象：球体位置看起来只差一点点，球的下半部分陷入边界平面，上半部分露在外面。
原因：这就是第 3.4 节说的坐标语义问题。你提供的坐标实际是球体接触点的位置，但脚本把它当球心用了。球体 Part 基准在球心，translate移动的是球心，不是底部或任意指定点。
解决：先确认每个坐标的语义。如果手头数据是接触点坐标，在代码里统一加上半径偏移，比如球体底面落在 z=0 平面，则z_center = z_contact + r。不要试图在个别球的坐标上手动改数值，一定在循环里统一处理，否则改错一个很难发现。
5.3 translate 第二次没反应：增量位移与固定约束的机制
现象：第一个球移动正常，想在同一 Instance 上再移动一次，脚本没报错，但球体纹丝不动，或位置完全和自己想的不一样。
原因：translate的向量是增量位移，不是绝对坐标。第二次调translate(..., vector=(10, 20, 30))是在当前位置上再叠加这个向量，而不是移动到 (10,20,30)。如果心里把它当绝对定位，结果就是位置偏到不知哪里。
解决：一个 Instance 只创建一次、只移动一次，目标值就是最终坐标。需要改变位置时，删除旧 Instance 重建，而不是重复移动。另外，如果 Instance 在装配里被加了固定约束，translate可能直接失效，这种情况也建议重建而不是尝试解除约束后移动。
5.4 脚本第二次运行报「模型已存在」或「Part 名重复」
现象：第一次运行正常，第二次运行报 KeyError，提示某个模型或 Part 已经存在。
原因：Abaqus 的mdb对象在进程内是持久的，同一个 CAE 会话里重复跑脚本，上次创建的模型和 Part 还在模型数据库里，名字冲突。
解决：脚本开头统一清理旧模型。核心三行代码：先判断模型名是否存在，存在就删除，再新建模型。Part 的清理可以依赖模型删除自动完成，不用逐个处理。只要每次跑脚本都从同一个模型名重建，就不会累积残留对象。
5.5 坐标写对了球却互相穿插：缺一个距离检查
现象：多个球装配后部分球体重叠穿插，网格生成时报几何干涉，或者后续接触仿真初始就穿透。
原因：坐标列表本身没做干涉校验。球心间距小于两球半径之和时，两个球必然穿透。如果这些球最终要参与接触分析或布尔运算，干涉是致命的；如果只是做简单的体积可视化，可能不明显。
解决：在批量装配脚本里增加距离检查逻辑，每放入一个球前，和已放球逐个计算球心距，小于半径和就报警。
placed = [] for i, (r, (x, y, z)) in enumerate(zip(radii, centers)): for (px, py, pz, pr) in placed: d = ((x-px)**2 + (y-py)**2 + (z-pz)**2) ** 0.5 if d < r + pr: print('警告: %s 与已有球重叠, 间距=%.3f' % (i, d)) placed.append((r, x, y, z))这段检查代码本身不复杂，复杂的是阈值怎么定。如果后续要做网格划分，建议判断条件从d < r + pr改成d < r + pr + gap，提前预留间隙，这个gap根据网格最小尺寸来定。否则即使不重叠，两球贴得太近也会影响网格质量。
6. 进阶玩法：把球体装配脚本变成参数化投放工具的 3 个技巧
6.1 用 CSV 管理参数，脚本只做执行不再写死数据
上一章提到的 CSV 输入，实际项目里我一般还会加一层：CSV 文件名通过脚本最上面的一个变量指定，方便切换数据集。列格式固定为radius,x,y,z，第一行可以是注释行，以#开头。这样每个工况一套 CSV，脚本完全不用改。几十个球的工况，用 Excel 维护坐标列表比在脚本里改数组直观得多。
6.2 随机不重叠投放：给「任意位置」加一个约束条件
有的需求不是给定坐标，而是给定一个区域，自动生成一组互不重叠的球。核心逻辑是随机生成候选球心和半径，检查与已有球不重叠才接受。
import random random.seed(42) def no_overlap(x, y, z, r, placed): for (px, py, pz, pr) in placed: if ((x-px)**2 + (y-py)**2 + (z-pz)**2) ** 0.5 < r + pr: return False return True placed = [] while len(placed) < 30: r = random.uniform(2.0, 6.0) x = random.uniform(0.0, 50.0) y = random.uniform(0.0, 50.0) z = random.uniform(r, 20.0) if no_overlap(x, y, z, r, placed): placed.append((r, x, y, z))注意random.seed(42)让随机序列可复现，同一个种子永远生成同一批位置，这在需要对比不同网格方案时非常重要。z = random.uniform(r, 20.0)保证了球体底面不会穿透 z=0 平面，这是落座语义的应用。
6.3 跑完先打印位置日志，再提交网格和计算
装配脚本跑完，我习惯把每个球的半径、球心坐标、编号打印成一份完整日志，并同步写一份 txt 文件存档。这样提交计算前可以快速核对，确认坐标和需求一致。其实我吃过一次亏：某天加班赶一个颗粒模型，坐标列表少了一行，最后一个球复用了上一个坐标，两个球叠在一起，接触计算跑了三天，结果完全不对。从那以后，我每次跑完脚本都会先看打印日志，确认球数、坐标数、半径数三者对得上，再往下一步走。这也是为什么示例代码里每创建完一个球就打一行 print，看着啰嗦，但关键时刻能救命。希望帮到你。
本文还有配套的精品资源，点击获取
