转载点物料颗粒运动的数值模拟

本资料包含pdf文件1个，下载需要1积分

Numerical Simulation of the M ovement of M ateriel Particles on theTransfer StationPiao Xianglan， Guo Yue(Institute ofEngineering ofYah bian University，Vanji Jilin，133002，China)Abstract：To analyze the movement of materiel particles on the fansfer station，theinteraction between particles.interaction between particles and boundaries and the energY terms inthe movement process，the procedure of particle flow code(PFC3D)based on discrete elementmethod(DEM)is used to build the discrete element model of the tram，；fer station.Numericalsimulation of the materiel Particles is directed towards the different velocities of the particle flOWand whether there is a ba用e in front of the roller or not.BY means of simulation.it iS found tohave some advantages：(1)The solid particlesmotion can be folowed beRer(the position，velocitY，average contact force and average unbalance force of each particle).The maximum ofthe average contact force increases along with the increase of particlevelocity，but the averageunbalance force has litle change；(2)Folowing the course of the energy，the transformation ofvarious kinds of energies can be looked into；and in all kinds of energies，kinetic energy variesmost；the bafl e will carry weight to kinetic energy，friction energy and strain energy 。

Key words： transfer station；bulk solid materiel；discrete element method；numericalsimulation收稿日期：2010-05-04；定稿日期：2010-07-02基金项目：延边大学博士启动基金资助项目 (012800-950201003)作者简介：朴香兰 (1964)，女，吉林延边人，教授，博士，主要研究方向为现代设计理论与方法。E。mail：piaoxl###ybu·edu·。

第 3期 朴香兰等：转载点物料颗粒运动的数值模拟转送站固体颗粒物料的实际运动十分复杂，除受到输送带的速度、滚筒半径、输送机布置形式、物料与输送带的摩擦系数等影响外，还受到物料内部运动的影响。Rovertst A.W Ll WensrichC.M .L2J M aton A.E.j 等分别从不同角度研究了物料颗粒的运动及边界的设计。周家海L4J、宋伟刚5J分别对带式输送机卸料时的物料运动轨迹做了研究，并用不同方法绘制了轨迹线。采用连续介质的力学方法，只能把散装物料颗粒群体作为- 个整体来考虑，无法分析散粒群体中每个颗粒的运动过程、相互作用及其对整体运动的影响。

本文利用基于离散元理论的颗粒流代码程序，建立了转送站物料运动的离散元模型，在此基础上研究了散料颗粒的运动过程、相互作用及其挡板对颗粒流能量损耗的影响，并模拟了卸料的整个过程，最后得出结论。

1 离散单元法 (DEM)简介离散单元法 (Discrete Element Method，简称DEM)是把不连续体分离为刚性元素的集合，使各个刚性元素满足运动方程，用时步迭代的方法求解各刚性元素的运动方程，继而求得不连续体的整体运动形态。

离散元法与连续介质力学方法相比有本质上的不同。它根据离散体的离散特性来建立数值模型，因此更符合离散物质本身的性质，在分析具有离散性质的物质材料方面表现出了极大的优越性。使用离散元法对颗粒材料进行模拟分析，可以直接获得大量离散物质的复杂行为信息，从而可以进-步假定和分析离散物质的细观结构，为准确预测和分析现有连续介质理论无法解释和分析的物质力学行为提供了可行的途径。

离散元法不需要过多的假设，使用简单的方程就可以对高度复杂系统的准静态和动态行为进行模拟。离散单元法的原理较为简单，只需要满足两组基本方程：第-是运动方程 (牛顿第二运动定律)；第二是物理方程 (力-位移关系)。

目前有二维、三维粘弹性接触模型、Hertz理论接触模型、填隙流体的湿颗粒作用模型等[ 。

2 数值模拟及结果分析2.1 转送站离散元模型的建立[ ]以 PFC3D 为平台，采用基本单元为球形颗粒，在限定范围内随机生成 1000个物料颗粒，颗粒半径在 2535mm之间。滚筒的模型是采用- 个圆柱面和两个有向面来构成，如图 1所示。

生成的每个颗粒在重力加速度下降落到带子上方，达到初始平衡状态，如图2所示。

图2 初始平衡状态2.2 参数设定设定仿真所需的-些参数，如几何特性、材料属性、摩擦系数、阻尼系数等，如表 1所示。

表 1 模拟参数2-3 数值模拟与分析2.3.1 卸料轨迹的模拟采用带速为 3m/s，当运行-定时步后，卸料过程仿真状态如图3所示。图中黄颜色表示颗粒间、颗粒与边界间的相互作用力常当滚筒前设有-与料流方向成 90。的挡板时，模拟物料的运动轨迹如图4所示。通过理论计算绘制的轨迹线如图 5所示l9]。由图3可见，在滚筒圆的水平轴线方向上，即z0的方向上，物料中心的轨迹离滚筒中心的距离约等于 1m，这与图 5物料中心图 学 学 报 201 3丘轨迹线基本相符。这说明仿真时设定的参数基本正确、模拟可信。

图 3 物料的卸料过程仿真图4 设有挡板的料流模拟图5 物料不同层的轨迹线2.3.2 不同的运行速度对颗粒之间相互作用的影响颗粒之间的相互作用可通过平均接触力和平均不平衡力来描述。平均接触力是指颗粒与颗粒之间以及颗粒与边界之间的接触力在运行时步范围内的平均值，而平均不平衡力是指颗粒所受外加作用力未达到平衡时的平均值。

取带速分别为 3m/s、4m/s和 5rn/s时，平均接触力和平均不平衡力分别如图6所示。图中纵坐标表示平均接触力(红线)及平均不平衡力(黑线)。共同的时步范围是 6.O05x10 ～1×10 s。从模拟数据可知，当颗粒速度增加时，平均接触力的最小值基本-致，平均接触力的最大值随颗粒速度的增加而增加。但平均不平衡力无论是最大值、最小值随速度的变化不大。

233 不同运行速度对能量消耗的影响不同颗粒运行速度时的能量历程如图 7，纵坐标表示各种能量值。与图7对应的能量界限值如表 2。从图7及表 2可见：能量变化最明显的是动能 (黑线)，其次是实体能 (红线)、摩擦能(深蓝线)、应变能 (绿色)和边界能 (天蓝色)。

动能计入所有粒子的平移及旋转运动，实体能是由全部体力作用在集合粒子上而积聚的能量。全部体力包括重力荷载以及所有外加作用力和力矩。边界能是由所有边界(wal1)作用在集合粒子上而积聚的能量。随着颗粒速度的增加，其动能急剧增加，而边界能的变化很缓慢。这是因为动能与颗粒速度的平方成正比。但对边界能来讲，由于速度的增加，边界对颗粒的作用几率相对于颗粒之间的作用几率反而减少，因此边界能的变化很缓慢。

图6 平均接触力和平均不平衡力第 3期 朴香兰等：转载点物料颗粒运动的数值模拟(b)图7 能量历程表2 能量界限值2.3.4 任意颗粒在给定方向的速度及位置的变化 当颗粒速度分别是 3、4、5m/s时，对应的纵向速度 (黑线)及纵向位移 (红线)如图 8所示。图中纵坐标表示速度及位移(m)。任选质心坐标为 (0，0，0.5)的颗粒，对应 3种颗粒速度的颗粒 ID分别为 901、929、678。图8(a)纵向位移在给定的时步 (6.00510 ～1.00310 )内由0.5281m 降到 0.3018m；图 80)纵向位移在相同时步内由0.5252m降到 0.2736m；图 8(c1纵向位移由0.5274m降到 0.385lm。说明纵向位移随颗粒速度的改变不大。而速度在给定的时步内，虽然速度变化趋势大致相同，但随颗粒速度的变化，纵向速度的波动也越来越大。这是颗粒速度的增加，使得颗粒之间、颗粒与滚筒前挡板碰撞机会增加的缘故。

(a) Ball ID901 Co) Bal ID 929 (c) Ball ID678图 8 颗粒的速度、位移变化图 学 学 报2.3.5 挡板对能量消耗的影响在输送机头部滚筒前离轴线 lm处，设置与料流方向成 90。的挡板后对能量历程进行跟踪，其能量变化如表 3。与图 7及对应的能量界限值表 2进行对比可知，有挡板时动能的变化是随着速度的增加急剧上升，尤其是动能的最大值变化很大；其次是应变能和摩擦能，实体能和边界能几乎没有变化。

表 3 能量界限值3 结 论对转载点物料颗粒运动进行数值模拟可知，平均接触力最小值随速度的变化不大，平均接触力的最大值随着速度的增加而增加。但平均不平衡力的最大值、最小值随速度的增加变化很校任意颗粒的垂直速度随着速度的增加出现波动现象。通过能量历程的跟踪观察实体能、边界能、摩擦能、动能、应变能等的变化，发现颗粒速度的增加将导致颗粒运动动能的急剧增加，其次是实体能，摩擦能、应变能和边界能。加了挡板之后能耗的增加主要体现在动能、应变能和摩擦能的变化。

基于 PFC3D 的转载点物料颗粒运动的数值模拟，能直观地展现料流的宏观运动，可获得物理实验所不能取得的信息，如物料颗粒运动的力尝速度场，可以实时追踪所有变量并能存储起来，如能量历程的跟踪等，还能考察单个颗粒的运动形态。