关键词：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球形氧化铝粉末（D10=8.2µm、D50=21µm、D90=47.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 Hz，1/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优化、后处理分析的完整仿真服务。