西北旱区残膜干法制粒装置机理分析与关键结构参数优化

《AgriEngineering》：Mechanism Analysis and Parameter Optimization of a Residual-Film Dry Granulation Device

　　

中国西北地区干旱农田回收的残留塑料薄膜因干法清理后仍残留较多土壤与秸秆，无法采用水洗熔融制粒工艺进行处理。本研究对一款自主研发的干法制粒装置进行优化，该装置通过摩擦诱导生热将此类薄膜转化为颗粒。首先，研究人员分析了喂入、切割与搅拌过程，确定动刀转速、定刀刃角和搅拌叶片上端长度为关键结构参数。其次，利用Fluent、LS-DYNA和EDEM（离散元法，Discrete Element Method）仿真确定各参数的合理范围。第三，以残膜利用率（Y₁）和制粒率（Y₂）为响应值开展单因素试验与Box-Behnken试验，并拟合回归模型（R²=0.984和0.989）。最优参数组合为665.86 r/min、45.58°和97.24 mm，模型预测Y₁=92.77%、Y₂=95.58%。田间验证结果分别为91.71%和94.08%，相对误差为1.14%和1.57%。研究结果可为残膜干制粒设备的结构设计与参数选取提供依据。

研究背景与问题
地膜覆盖（plastic-film mulching）是西北干旱半干旱区增产的有效措施（保墒、增温、抑草），但我国地膜年用量已超200万吨、约占全球总量的75%，大量回收残膜被堆存、填埋或倾倒，经风力和水力重新分布后造成跨界污染，降低土壤容重、孔隙度与透气性并污染秸秆基牲畜饲料；将回收残膜转化为可运输产品由此成为闭合残膜回收循环的关键。当前主导的转化技术为水洗熔融制粒（melt granulation），要求原料杂质含量低于约0.2%，必须清洗后制粒，耗水量大并产生污水；滚筒筛、圆盘筛、风选机等干法除杂（dry impurity removal）设备虽可替代水洗，却均无法完全清除细土和粉尘，干法清理后残膜杂质仍远超水洗制粒工艺的容忍限度。这一“干法清理高效、水洗制粒工艺”的错配，导致现有设备喂入性能差、制粒效果不理想，参数调优仍依赖试错法。

研究开展、结论与意义
针对上述问题，研究人员以自主研发的残膜干法制粒装置为对象，按照理论受力分析、数值仿真、响应面试验与样机田间验证四个步骤开展研究：解析喂入、切割、搅拌三个过程的受力以识别控制性结构参数，通过仿真再现气流场、切割变形与颗粒运动以缩小参数范围，进而明确动刀转速、定刀刃角和搅拌叶片上端长度对残膜利用率（Y₁）和制粒率（Y₂）的影响并确定最优组合。研究建立了切割阻力与刀具刃角、膜厚之间的定量关系（F_f=μF_n2+0.5kl₂²sinθ），提出了基于仿真的参数范围确定流程和结构参数—性能指标回归模型，并给出经验证的最优参数组合（665.86 r/min、45.58°、97.24 mm）。相比前期仅报道结构设计的方案，本研究揭示了干法制粒装置的工作机理，建立了定量化的参数—性能关系并提供了经验证的优化流程，可为残膜干制粒设备的结构设计与参数选择提供依据。

关键技术与方法
主要方法包括：理论受力分析，即喂入过程气流推力与压差阻力的平衡模型、基于胡克定律（Hooke's law）的切割“压缩—切割”两阶段受力模型以及搅拌推料的几何体积模型；Fluent（计算流体动力学（CFD），标准k-ε湍流模型）模拟喂入口气流场，LS-DYNA显式动力学（等效厚度叠加简化模型、最大等效塑性应变失效准则）模拟定刀切割，EDEM（离散元法（DEM），Hertz-Mindlin接触模型）模拟搅拌颗粒运动，三者相互独立、无数据传递；基于响应面法（RSM）的单因素与Box-Behnken试验，以残膜利用率（Y₁）和制粒率（Y₂）为响应值拟合二次回归模型并多目标寻优。试验样本为联合残膜回收—秸秆粉碎机械回收并经干法除杂的残膜捆。

3.1 单因素试验与结果分析（Single-Factor Experiment and Result Analysis）
以动刀转速、定刀刃角和搅拌叶片上端长度为因素开展单因素试验（每组三次重复，相对偏差均小于2%）。当定刀刃角为45°、搅拌叶片上端长度为100 mm时，残膜利用率随动刀转速升高而下降，制粒率先升后降，综合确定动刀最优转速范围为550~850 r/min；动刀转速为700 r/min时，定刀刃角优选范围为42.5~47.5°、搅拌叶片上端长度优选范围为87.5~112.5 mm。

3.2 多因素试验与结果分析（Multi-Factor Experiment and Result Analysis）
开展Box-Behnken响应面试验：残膜利用率范围为90.22%~93.66%（平均92.12%），制粒率范围为90.59%~96.20%（平均93.87%）；利用率模型R²=0.9842（调整R²=0.9639），制粒率模型R²=0.9890（调整R²=0.9749）。方差分析（ANOVA）表明，影响利用率的因素显著性排序为动刀转速、搅拌叶片上端长度，影响制粒率的排序为动刀转速、定刀刃角、搅拌叶片上端长度。响应面分析显示：定刀主要影响切割效率、对造粒筒内部流场影响小，故刃角与转速的交互作用较弱；动刀转速升高会增强出口气流扰动、吹走残膜，使利用率下降；制粒率先升后降，转速超过700 r/min后摩擦热过剩而搅拌叶片传热能力有限，导致造粒筒底部颗粒软化、粘附和团聚。

3.3 验证试验与结果分析（Verification Experiment and Result Analysis）
采用复合满意度函数D=√(d₁·d₂)（两响应等权重）进行多目标优化，最优参数组合为动刀转速665.86 r/min、定刀刃角45.58°、搅拌叶片上端长度97.24 mm，模型预测利用率92.77%、制粒率95.58%；调整为666 r/min、46°、97 mm后，在闭环残膜处理示范基地（新疆生产建设兵团第六师芳草湖团场）重复5次验证，平均利用率91.71%、制粒率94.08%，相对误差分别为1.14%和1.57%。相比原设计上限1000 r/min工况（利用率约87.4%、制粒率约90.3%），优化组合在转速降低三分之一的同时将两项指标分别提高4.3和3.8个百分点。

3.4 讨论（Discussion）
受力分析表明定刀刃角是切割过程的关键几何变量：切割阻力随sinθ增大，小刃角虽减小受力但削弱撕裂，大刃角改善撕裂却降低切割效率；切割仿真再现了这一权衡，刃角由30°增至60°时切割时间延长、平均变形增大，裂纹尖端变形因局部断裂释放应力而非单调变化。气流仿真显示喂入口气流速度随动刀转速近似线性增加，高转速虽改善切割与生热却将残膜反吹而降低利用率，与单因素结果一致。与前期设计研究相比，本研究将三个结构参数与两项性能指标定量关联，动刀转速为主导因素（F=337.04和58.14），“转速×刃角”“转速×叶片长度”交互作用显著；最优参数并非单独最大化任一指标，而是在限制反吹与能耗的同时保持制粒率高于94%。研究同时指出局限：EDEM仅追踪5个示踪颗粒，搅拌分析为说明性而非统计性；试验仅采用单一地区、单一地膜类型，设备为间歇运行需批间冷却；回归模型为局部二次近似，不可在测试范围外推。

结论（Conclusions）
（1）通过分析喂入、切割与搅拌过程，识别出动刀转速、定刀刃角和搅拌叶片上端长度三个控制性结构参数；推导的切割阻力模型F_f=μF_n2+0.5kl₂²sinθ表明切割阻力随刃角与膜厚平方增大，解释了刃角须保持适度范围的原因。
（2）仿真确定三个参数的范围：喂入临界转速约1000 r/min（范围400~1000 r/min），40°~50°刃角兼顾切割效率与撕裂性能，搅拌叶片上端长约100 mm时颗粒向上运动速度最持续；三次仿真相互独立、代码间无数据传递。
（3）单因素试验将参数范围缩小至550~850 r/min、42.5°~47.5°和87.5~112.5 mm；Box-Behnken模型显著（R²=0.984和0.989）且失拟（lack of fit）不显著，动刀转速对两项指标均为主导因素。
（4）满意度最优组合（665.86 r/min、45.58°、97.24 mm）预测Y₁=92.77%、Y₂=95.58%，田间验证结果为91.71%和94.08%（相对误差1.14%和1.57%）；优化设定以三分之一的较低转速将两项指标分别提高4.3和3.8个百分点。主要局限为EDEM分析仅采用5个示踪颗粒，且试验仅使用单一地膜类型与单一地区的样本，全颗粒统计分析与多地区试验留待未来工作。