Abaqus喷丸强化：Python二次开发构建随机弹丸和粒子生成器，实现结果数据自动提取

导读：航空发动机关键零部件的疲劳失效一直是亟待解决的问题，提高材料疲劳强度可以减少部件维护次数从而降低成本，同时可以减少由于疲劳而导致的发动机故障。

其中喷丸强化、激光冲击强化等表面强化技术通过表面改性可以显著提高材料疲劳寿命，具体表现为在材料表面引入残余压应力抑制裂纹的萌生和扩展；提高表面硬度从而提升材料的抗微动疲劳性能；细化晶粒使表面强度提升。然而，喷丸强化工艺参数繁多，不同材料对应的最优喷丸工艺不同，实验条件下确定材料喷丸工艺需要消耗大量材料（待强化材料和ALMEN试片），同时测量残余应力和硬度等表征结果也会消耗大量人力、物力。

本文介绍通过python二次开发构建弹丸和通过粒子生成器生成弹丸是喷丸强化仿真常用的仿真手段，以及这两种常用建模手段进行单弹丸的喷丸对比。诚邀朋友们关注我在仿真秀官网原创首发机械喷丸/随机弹丸/粒子生成器/二次开发/后处理视频教程《ABAQUS机械喷丸参数化建模仿真12讲：Python二次开发构建随机弹丸和粒子生成器》，详情见后文。

一、数值模拟在喷丸强化中的应用

数值模拟手段通过大量的迭代发展，已经在航空航天、汽车和医学等领域得到广泛应用。近年来涌现出大量文献研究仿真手段在喷丸强化领域的应用。喷丸强化时间短，应变率高，变形较大的瞬时非线性特点导致其计算只能使用显式求解。ABAQUS和LS-DYNA在显式计算领域具有独特优势，是喷丸强化仿真计算的首选。激光冲击强化在调整步长较小的情况下可以采用隐式求解，求解结果相对显式求解精度较高，但求解效率较低。

复合曲率试件的有限元模型（图3-来自迈特莱斯官网）

根据有限元分析指导加工的试样（图4-来自迈特莱斯官网）

二、ABAQUS在喷丸强化仿真中的优势

ABAQUS软件相比LS-DYNA而言灵活度更高，其可以通过python进行前后处理二次开发，Fortran进行求解二次开发以及修改关键字进行输入文件的二次开发。例如喷丸强化首先可以通过编写python脚本进行随机弹丸的建模和后处理残余应力和粗糙度的自动化提取，其次可以通过修改inp文件的关键字进行粒子生成器的定义，通过python和VDLOAD子程序的协同作用进行激光冲击前处理的参数化建模、光斑的路径定义、激光工艺参数的参数化定义以及脉冲曲线的参数化定义。

发动机曲轴的横截面图，显示了由于偏心引起的不平衡残余应力，及由于不平衡残余应力而预测的变形（图5-来自迈特莱斯官网）

三、单弹丸的喷丸对比

通过python二次开发构建弹丸和通过粒子生成器生成弹丸是喷丸强化仿真常用的仿真手段。本文就这两种常用建模手段进行单弹丸的喷丸对比，单位制为mm、t、s。

1、构建靶材模型

首先构建靶材模型，在ABAQUS软件构建1mm×1mm×0.5mm的靶材，网格大小为0.01，材料选用工业中常用的中国牌号为TC4的钛合金，其材料参数如下。

通过CAE软件构建半径为0.15mm的铸钢丸，并对其进行刚体约束。其密度为7.85e-9，杨氏模量和泊松比分别为2.1e5和0.3，网格大小为0.01。冲击速度设定为80000mm/s。通过GUI界面构建的模型如下。

2、粒子生成器单弹丸模型

构建粒子生成器的单弹丸模型有很多办法，控制粒子生成器的流率可以生成单个弹丸，其流率协调公式如下。

其流率为粒子速度、密度和喷嘴阈值的函数关系。然而，通过控制流率的方式去生成单个弹丸的过程较为繁琐，本文直接通inp文件的节点坐标去构建粒子生成器生成的弹丸。粒子生成器生成弹丸的方式本质上属于生成随机节点，弹丸表示生成颗粒的球心位置，边界条件施加在生成的节点上。

因此，首先在inp文件中构建part以及相应的节点号以及节点所在空间位置。

下面对节点的单元类型和集 合进行定义。粒子的单元类型为PD3D，粒子集 合定义为DEM。

对粒子进行相关属性的定义。对上述集 合DEM内的粒子定义其密度为7.85e-9，alpha阻尼先设定为0。0.15表示粒子半径为0.15。

对粒子进行杨氏模量和泊松比的定义，对于刚体而言，其并无太大价值，但其对粒子的法向接触行为的定义具有重要意义。ABAQUS帮助文档详细解释了Hertz接触与Johnson-Kendall-Roberts adhesive接触，具体公式计算见帮助文档。

下面对粒子表面进行定义，对其表面命名为DG_Surf，对应的集 合为DEM，结束颗粒的part建模。

下面对接触进行定义。接触定义的关键字很常见，其关键字为：

本文采用的接触为，法向硬接触、切向摩擦系数为0.2。

粒子生成器生成的颗粒无需划分网格，因此，直接对粒子进行装配，其关键字为：

其含义为：降part的nozzle即粒子添加到装配体中。因为我们在part模块就已经对粒子的相对空间位置坐标进行了定义，在装配中无需定义其空间坐标。

下面就是对粒子的边界条件进行定义。对于刚体而言，在GUI界面我们通过对参考点施加相应的边界条件，对于喷丸工况而言，需要对参考点施加预定义场的初始速度。其逻辑是首先创建参考点的集 合，然后对集 合赋予相应的边界条件。因此，我们首先需要对粒子创建相应的集 合。由于粒子依附于节点，我们对节点进行set的集 合定义。

将其命名为nozzle-1_DEM。

下面进行初始速度的施加。3表面Z方向，-80000表示速度为80000mm/s，方向为Z轴负方向。

在计算时，粒子生成器经常会报错Zero element time increment estimate，这个时候我们讲增量步设定为固定增量步即可解决。粒子生成器弹丸模型构建结果如下。

3、两种建模方式的求解结果

下面对比一下两种方式的求解结果，左边为手动构建弹丸结果，右边为粒子生成器弹丸冲击结果。

（1）对于S11方向的残余应力而言，粒子生成器构建的弹丸冲击后的残余应力值较大。

（2）对于等效塑性应变而言，粒子粒子生成器冲击后的变形更大，最重要的一点是其弹坑形状非常怪异。

（3）对于U3位移结果而言，粒子生成器的弹丸其形状不规则，并且弹坑更深。

按道理而言，两种建模方法构建的弹丸其冲击结果应该相同，然而其结果却令人意外。上面我们提高阻尼，因此这可能是阻尼定义不当引起的怪异。本文对比了阻尼为0，0.1，0.5，1和7情况下，其塑性应变和残余应力。

提高阻尼，就残余应力而言，0.5的阻尼下残余应力接近手动建模下的弹丸冲击情况，但塑性应变仍无法满足。弹坑的不规则变形和残余应力显著影响后续的残余应力大小和粗糙度。

就残余应力深度而言，粒子生成器生成的弹丸也无法满足深度要求。

绘制了动能随时间变化的情况。可以发现动能在碰撞时两者的变化不同。这可能与阻尼有关，合适的阻尼值是粒子生成器应用的关键。作者目前看到的很多论文里面的粒子生成器并没有谈其阻尼的构建，然而所有的论文出乎意料的是，其残余应力与粗糙度均与实验相同。

4、粒子生成器在其他领域的应用

作者并未深入研究阻尼的定义，上述定义的阻尼下粒子生成器产生的粒子进行喷丸强化可能并不会产生理想的效果，但这并不影响粒子生成器在其他如岩土工程方面的应用。

然而，如果实际工况是大范围的如10mm×10mm的平板喷丸强化，这个时候所需弹丸可能高达几十万，通过python生成随机位置的颗粒并不现实，这对显卡有更高的要求，这个时候仍然需要采用粒子生成器作喷丸强化，这个时候有两种选择：查阅资料定义合适的颗粒阻尼或者在上述结果乘以相应的矫正系数，以贴合实验，但是第二种做法的合理性有待研究。

粒子生成器可以生成的粒子大小可以通过修改相应的函数去实现多种特定半径粒子的生成，一定半径范围内随机半径颗粒的生成等等。如果通过python二次开发，则还需定义手写几百行代码，而只需修改inp文件的关键字即可实现随机粒子的生成。

粒子生成器生成三种特定半径大小的弹丸。

粒子生成器生成一定半径区间的随机半径大小的弹丸。

此外，通过结合python和inp文件的二次开发，可以随机生成互不干涉的颗粒的球心坐标，进而写入inp文件，通过上述关键字的增改，完成颗粒的建模。下面是生成半径分别为0.04和0.06的空间中随机分布的粒子生成器颗粒。

本贴只是个人根据既有结果的初步探索研究，因水平有限，若内容有误，还请读者私信评论更正，以提高本贴质量。

四、我的ABAQUS喷丸强化应用教程

由笔者原创，且独家首发仿真的视频教程《ABAQUS机械喷丸参数化建模仿真12讲：Python二次开发构建随机弹丸和粒子生成器》本课程为ABAQUS机械喷丸参数化建模仿真精讲，对python构建随机弹丸和粒子生成器进行了step by step讲解。研究人员通过学习此课程，可以学会机械喷丸的建模和插件制作，同时通过脚本提取相应的残余应力和粗糙度，解放双手，提高效率。

以下是课程大纲（可加餐内容）

《ABAQUS机械喷丸参数化建模仿真12讲：Python二次开发构建随机弹丸和粒子生成器》

注意：课程中ABAQUS版本是2022。

1、您将得到

（6）可以根据VIP群用户需要，酌情加餐内容或者直播。