李郁冯云杰

武汉理工大学武汉

摘要：为了分析研究螺旋卸船机水平螺旋至溜筒过程的设备性能，选取50mm煤颗粒为作业对象，建立水平螺旋与溜筒的离散元仿真模型，完成了水平螺旋在上倾20°、水平、下倾20°等3个位置上，转速为rpm与前方供料t/h条件下的输送过程仿真。通过对比水平螺旋输送质量流量，轴向输送速度、流动状态、以及溜筒内流动轨迹与溜筒内壁切向累积接触能量，分析水平螺旋与溜筒间的转接通过性、水平螺旋转速及溜筒内壁磨损情况。得出结论：水平螺旋转速满足生产率要求，水平螺旋与溜筒转接处通过性良好，并给出溜筒内壁需要进行耐磨加工的具体位置。

关键词：螺旋卸船机；水平螺旋；溜筒；离散元法；仿真

中图分类号：TH文献标识码：A文章编号：-（）21--05

0引言

螺旋卸船机是一种高效环保的散货码头作业设备，其输送方式主要借助螺旋输送，即利用无挠性牵引构件实现散货的连续输送，可作业各种粉状、粒状及小块状物料[1]。对于螺旋输送的工作性能及颗粒运动，国内外学者已经通过离散元仿真方法进行了深入研究。Shimizu等[2]首先通过离散元法模拟螺旋输送机，分析各种螺旋输送机模型的力矩、功耗等，并与理论公式进行了对比分析。Cleary等[3]以螺杆转速，填充水平，倾斜度为性能参数，用离散元法预测了螺旋输送机的性能，分析了颗粒速度，质量流量，能耗，功耗方面的变化。LatoPezoa等[4]通过离散元仿真结合聚类分析等数值分析方法评估了3种设计变量对螺旋输送机颗粒混合性能的影响。余书豪等[5]研究了螺距、供料方式、螺旋叶片等性能参数对螺旋输送机性能的影响，并建立响应面模型进行优化设计，但对于诸如具体卸船设备的研究还较为缺乏。

本文在前期工作基础上[7-15]，通过离散元仿真实验，分析螺旋卸船机水平输送段及转接溜筒过程中的颗粒流运动状态，进而提出对设备结构及生产操作条件的改进意见，研究方法与结果可对相关领域的设备性能优化提供一定借鉴。

1离散元仿真建模

1.1离散元颗粒接触模型

以某电厂常用煤物料为研究对象，将其简化为直径50mm的球颗粒，颗粒间接触模型使用Hertz-Mindlin无滑移接触模型，如图1所示。

对颗粒进行受力分析，采用牛顿第二定律建立颗粒运动方程，对运动方程积分，不断更新颗粒的位置[6]，即有

式中：mi、Ii为颗粒i的质量和转动惯量，vi、ωi为颗粒i的平动速度和转动速度，ni为与颗粒i接触的颗粒数，Fijc为颗粒i与j之间的接触力，Fig为颗粒i所受重力，Fif为颗粒受流体的作用力，Mit为切向碰撞力产生的力矩，Min为法向碰撞产生的力矩，Mir为滚动摩擦力矩。

图1HM无滑移接触模型

1.2仿真参数

通过文献[8]及试验标定，确定仿真材料本征参数与接触参数，如表1、表2所示。

1.3水平螺旋输送与溜筒模型

仿真研究对象选取某额定生产率为t/h的螺旋卸船机，其水平螺旋转速为rpm，螺旋叶片采用双头螺旋，螺距为mm，设计仿真实体结构如图2所示。螺旋卸船机在正常工作中，水平螺旋会根据料堆位置进行一定幅度变化（变幅范围为-20°～+20°）。

为分析水平螺旋位于不同位置时的物料输送状态及与转接过程中的物料通过性，选取水平螺旋倾角分别为-20°、0°、20°下的3种特征位置，定义为A、B、C等3种输送模式。其中，输送模式A为水平螺旋倾角-20°状态，输送模式B为水平螺旋倾角0°状态，输送模式C为水平螺旋倾角20°状态。为了观察物料流运动状态，在水平螺旋的中间输送段与溜筒承接处分别设置质量流量传感器Q1和Q2，如图3所示。

1.进料口2.水平输料管3.水平螺旋4.转料口5.溜筒6.导料板7.出料口

图2水平螺旋与溜筒结构简图

图3水平螺旋与溜筒仿真模型

2结果与分析

2.1水平螺旋与溜筒转接通过性分析

观察质量流量传感器，通过提取稳定输送状态下质量流量数据，分别得到水平输送段和溜筒承接处两个位置处的平均质量流量1Q与2Q。3种输送模式下稳定输送状态下的质量流量如图4～图6所示。

从表3仿真结果分析可得，3种模式下水平输送段和溜筒承接处两个位置处质量流量的相对偏差分别为0.48%、0.28%、0.08%，说明在3种模式下的水平螺旋转溜筒处均未发生明显颗粒堵塞，颗粒通过性良好。此外，输送状态C即水平螺旋上倾20°时，通过性达到最佳。

图4模式A稳定时段流量

图5模式B稳定时段流量

图6模式C稳定时段流量

2.2水平螺旋输送状态分析

3种模式下水平螺旋不同倾角下的稳定输送状态如图7所示，其平均轴向速度见图8。从图中可知，输送模式A即水平螺旋向下倾斜20°时，物料颗粒在输送管底部密集，个别颗粒有径向跳跃，颗粒间速度无明显差别；输送模式B即水平螺旋处于水平时，部分颗粒出现明显径向跳跃；输送模式C即水平螺旋向上倾斜20°时，物料分布于输送管整个断面，物料流整体有较大径向速度。

同时，图8反映出3种模式下的平均轴向速度存在一定差别。其中，输送模式C下的水平螺旋物料轴向速度明显偏小，这是由于该模式下螺旋转速过快，明显大于标准转速，导致物料颗粒产生较大的垂直于输送方向的径向速度分量。螺旋叶片在进行水平输送的同时，对物料产生了搅拌干扰，以致输送效率的明显降低。

要达到理想的输送状态，保证颗粒流平稳移动，采用输送模式C，在卸船机额定生产率下工作时，建议适当降低水平螺旋转速。

图7水平螺旋输送颗粒流动状态

图8水平螺旋输送轴向输送速度

2.3溜筒内颗粒运动

1）出口处颗粒运动速度

表4列出了3种输送模式下溜筒出口处测得的颗粒平均速度，由此得知3种输送状态下溜管出口处颗粒运动速度和方向无明显变化，稳定在10m/s左右。

2）溜筒内颗粒运动轨迹

3种输送模式下溜筒内颗粒流运动轨迹如图9所示。图中标识为颗粒流与溜筒内壁发生集中接触的区域，输送模式A下颗粒流运动平滑，较为集中，无明显撞击分散；输送模式B及C下，颗粒流首次冲击位置改向圆弧之下，颗粒流与内壁发生剧烈冲撞，导致溜筒中段区域物料颗粒较为分散，颗粒运动方向紊乱，加剧下部的碰撞摩擦，因此建议适当增大溜筒入口处圆弧尺寸。

图9溜筒内颗粒运动轨迹

3）溜筒内壁磨损分析

在长期作业下，高速颗粒流会对溜筒内壁造成一定磨损。为分析得出不同模式下溜筒发生磨损的关键位置，本次研究中提取了溜筒内壁切向累计接触能量，用以表征各部位的相对磨损量，即通过计算几何模型三角形网格上累积的力和接触来计算得出几何体某部位的相对磨损量。

图10所示为3种输送模式下颗粒稳定流动至14s时溜筒各部位的相对磨损量。可以观察到溜筒磨损区域与颗粒流冲击区域一致；3种输送模式下溜筒中下部的磨损区域基本相同，溜筒入口处的磨损区域略有不同，但因该区域颗粒速度有限，其总体磨损量较小。

分析得出：对于几处易磨损位置，可在溜筒中间直管全内壁，下部斜管左侧壁，左导管右侧壁，右导管左侧壁增加耐磨层，从而抵抗泻料过程中的磨损。

图10溜筒内壁切向累积接触能量图

3结论

1）水平螺旋各变幅位置处，水平螺旋转溜筒处物料通过性良好，不会发生积料堵塞现象。

2）水平螺旋转速满足生产率要求，但在C模式下，螺旋转速过快，造成显著径向跳跃，可适当降低水平螺旋的转动速度以提高输送效率。

3）增大溜筒入口位置的圆弧尺寸可在B、C输送模式下使溜筒中的颗粒流更加平滑，减少碰撞；3种模式下溜筒出口处的颗粒速度大小与方向基本一致。

4）建议在溜筒中间直管全内壁，下部斜管左侧壁，左导管右侧壁，右导管左侧壁分别进行耐磨处理，以减缓溜筒内壁的磨损。

参考文献

[1]蒋琼珠.连续运输机[M].北京：人民交通出版社，.

[2]ShimizuY，CundallPA.Three-dimensionalDEMsimulationofbulkhandlingscrewconveyors[J].JournalofEngineeringMechanics，，(9)：-.

[3]OwenPJ，ClearyPW.PredictionofscrewconveyorperformanceusingtheDiscreteElementMethod(DEM)[J].PowderTechnology，()：–.

[4]MiladaPezo.Discreteelementmodelofparticletransportandpremixingactioninmodifiedscrewconveyors[J].PowderTechnology，()：-.

[5]余书豪.垂直螺旋输送机输送特性及螺旋体的优化研究[D].绵阳：西南科技大学，.

[6]LiJ，KuipersJAM.Gas-particleinteractionsindensegasfluidizedbeds[J].ChemicalEngineeringScience，，58(3)：-.

[7]李郁，李勇智.螺旋气力取料机的工作原理及试验研究[J].起重运输机械，(8)：25-27.

[8]夏鹏，李郁，杨公波.散粒物料堆积角离散元仿真研究[J].起重运输机械，(2)：-.

[9]杨公波，李郁，陈定方.基于EDEM仿真的斗轮堆取料机取料机理研究[J].武汉理工大学学报，，38(3)：-.

[10]叶方平，李郁，胡吉全，等.颗粒凝聚力对筒仓装卸过程影响的研究[J].武汉理工大学学报，，38(7)：91-96.

[11]叶方平，李郁，胡吉全，等.基于颗粒动力学理论的气力输送特性研究[J].武汉理工大学学报，(5)：10.

[12]FangpingY，YuL，JiquanH，etal.InvestigationofParticleDynamicsinaDiscRotatingDevicebyMeansofExperimentsandNumericalSimulationsUsingDEM[J].JournalofChemicalEngineeringOfJapan，，51(8)：-.

[13]LiY，YangG，HuJ，etal.DischargingTrajectoryoftheBucketWheelStacker-ReclaimerBasedonDEM[C]//InternationalConferenceonDiscreteElementMethods.SpringerSingapore，.

[14]邓聪，李郁.考虑气流影响的斗式提升机卸料过程仿真研究[J].起重运输机械，(15)：53-58.

[15]邱俊霖，胡吉全，李郁.伸缩溜筒下料仿真分析及其结构优化[J].起重运输机械，(11)：53－58.