关键词EDEM粉末铺展刮刀振动PBF仿真偏度峰度一、文章简要介绍本文基于Journal of Manufacturing and Materials Processing期刊2024年CC BY 4.0开放获取发表的研究论文介绍利用离散元法DEM对粉末床熔融PBF增材制造中刮刀铺粉过程进行仿真优化的方法。研究采用Altair EDEM软件结合田口实验设计Taguchi DoE与方差分析ANOVA系统研究了铺粉速度、刮刀振动频率/振幅和后角对氧化铝粉末层质量的影响规律提出了以表面偏度Ssk和峰度Sku作为等效铺层质量评价指标的新思路。二、仿真步骤步骤一粉末建模与参数标定。选用SA-ZL-20球形氧化铝粉末D108.2µm、D5021µm、D9047.5µm通过对数正态分布精确描述粒径分布。颗粒间及颗粒-壁面的静摩擦、滚动摩擦、恢复系数均通过实验标定休止角仿真偏差仅1.6%验证了粉末模型的真实性。选用Hertz-Mindlin with JKR粘附接触模型基于Tabor参数分析确定接触模型适用性。步骤二仿真几何搭建。模拟区域为1mm×1mm方形基板理论层厚100µm。刮刀为类立方体形状含100µm水平底边和可调后角0°、5°、10°。刮刀叠加水平平移0.01/0.05/0.1 m/s与垂直振动500/1000/2000 Hz1/2.5/5 µm的复合运动。初始时刻在虚拟粉末工厂内按粒径分布生成氧化铝颗粒自由落体形成铺粉斜坡。步骤三L27田口正交实验设计。四因素三水平的全因子设计通过L27正交阵列缩减为27组仿真。每次仿真完成后导出粉末层的颗粒坐标与粒径数据通过C语言脚本进行射线投射计算层顶面点云1001×1001网格节点据此量化四个层质量指标层厚偏差LTD、表面覆盖率SCR、均方根粗糙度Sq-RMS和真实堆积密度PDtr。步骤四表面形貌统计分析。在Sq-RMS基础上进一步计算算术平均高度Sa、表面偏度Ssk和峰度Sku。通过最小二乘法建立三种归一化质量指标的加权平均值与Ssk/Sku之间的线性关系以决定系数R²优化权重组合验证偏度与峰度作为等效铺层质量评价指标的有效性。三、结果解读图1 模拟几何模型复现结果图1为EDEM粉末铺展仿真几何模型包含基板、刮刀、四周边界板及虚拟粉末工厂半透明立方体刮刀底部设有100µm水平底边以增强振动压实效果。图2 粉末自由落体与斜坡形成复现结果图2展示粉末生成后的自由落体过程初始时刻左粉末在虚拟立方体内按粒径分布生成随后在重力作用下自由落体至基板上中最终在刮刀前方形成待铺展的粉末斜坡右。图3 铺粉层俯视图按高度着色复现结果图3为铺粉完成后粉末层的俯视图颗粒按中心Z坐标着色红色高位、蓝色低位直观展示了铺粉层的厚度均匀性。颗粒高度分布可用于量化层厚偏差和表面粗糙度。图4 不同铺粉速度的层质量对比复现结果图4对比了不同铺粉平移速度下的层表面质量高速0.08 m/s上导致明显沟槽和覆盖缺陷低速0.01 m/s下可获得均匀致密的粉末层直观印证了铺粉速度对层质量的显著影响ANOVA贡献率87.5%。图5 三类质量指标与表面偏度的关系复现结果图5给出了三类归一化质量指标层厚偏差、覆盖率不足率、RMS粗糙度与表面偏度Ssk的关系。Ssk越负即深谷多、表面平整质量指标的加权平均值越小铺层质量越高线性拟合R²0.887。图6 三类质量指标与表面峰度的关系复现结果图6给出了三类质量指标与表面峰度Sku的关系。Sku越大加权平均值越小铺层质量越高。线性拟合R²0.950表明峰度是比偏度更可靠的等效质量指标。图7 振动辅助铺粉与无振动层质量对比复现结果图7直观对比了无振动铺粉左与最优振动参数铺粉右2000Hz/5µm的层质量差异。振动辅助使表面缺陷大幅减少覆盖率从98.6%提升至99.3%RMS粗糙度从20.9µm降低至15.7µm。四、我们提供的服务我们精通EDEM、Rocky等离散元软件及其与CFD、FEA的多物理场耦合在粉末铺展、颗粒流动、磨损破碎、混合搅拌等场景拥有丰富的仿真复现经验。无论是增材制造粉末床质量优化、制药粉末混合均匀性分析还是矿山破碎机衬板磨损预测我们均能提供从模型搭建、参数标定到DOE优化、后处理分析的完整仿真服务。