滚筒杀青机杀青过程的温度场分析.docVIP

滚筒杀 青机杀 青过程 的温度场分析 茶叶杀青要做到“杀透、杀匀”，通过模拟茶叶杀青的过程.可以得到 杀青过程中温度场的变化情况，对于改进杀青机结构和运动参数，提高杀青 效果至关重要?本章将对杀青过程中的温度场及其影响参数进行分析。 4.1滾筒杀青机模型的构建 茶叶杀靑数值模拟主要包括以下步骤：建立杀青机和茶叶几何模型，定 义杀青机和茶叶颗粒材料属性，定义约束条件和力学边界条件，加载载荷, 设置求解条件，求解和数据分析等。其流程如图4?1所示。下面将分别介绍 主要步骤的具体内容和PFC3D数值模拟的相关设置。 图41离散元数值模拟流程图 研究以用于名茶生产的三段电加热6CST50型滚筒杀青机为原型展开， 梵基本参数如表4」所示，杀青过程中杀青叶与周围环境之间发生着复杂的 物理和化学反应，因此在进行数值模拟分析时，忽略一些对温度场、应力场 彫响较小的因素。对离散元模型进行适当简化处理，并进行如下假设： 1、由于所使用的PFC3D版木不支持并行计算，且能处理的模型的基本 单元为球体，为研究方便和缩短导热计算时间，采用木构参数与茶叶相同的 当量球形颗粒取代杀青叶，使其导热过程近似一致； 2、忽略茶叶颗粒之间的热对流: 3、受软件自身限制，不能分别设置筒壁和颗粒的热传导率，因此采用 平均热传导率进行研究。 本文采用的滚筒杀青机参数如表4?1° 表4-1 6CST50烈滚筒杀青机基本参数 名称 基本参数 筒长（mm） 2200 樂厚（mm） 2 滾筒外径（mm） 500 转速（r/min） 33 滚筒倾角（°） 2 表中：为了缩短杀青仿真时间，滾筒倾角取2。 通过对滚筒中正在杀青的茶叶进行动力学计算分析⑹】【刑，得到与茶叶本 构参数相同的当量球形颗粒来替代杀青叶，在杀青过程数值模拟中，茶叶颗 粒当量直径为0?5?0.8cm,呈正态分布，如图4-2所示。利用PROE建立滚 筒杀青机模型，如图4?3所示，导入PFC3D,再由程序生成不同大小的茶叶 颗粒，杀青机杀青仿真模型如图4?4所示。 0.0080 0.0075 0.0070 0.0065 0.0055 0.0050 0 200 400 600 800 1000 1200 颗粒分布 图4?2茶叶额粒直徨分布图 !%| 4-3滚简杀冇机仿真模空 图4?4滚筒杀青机仿真模空 4.1.2材料属性定义 材料的力学性能参数和热物理性能参数能在很大程度上影响杀青数值模拟 的结果和计算过程，本文采用的茶叶颗粒本构参数如表4?2。 表4?2杀青叶模型本构參数 名称 参数 颗粒密度 532.4 颗粒间摩擦系数 1.0 颗粒切向刚度＜N/m) 1.0e4 颗粒法向刚度(Min) 1.5e4 颗粒与璧面摩擦系数 0.8 颗粒比热(KJ/kgQ 3.44 颗粒线性热膨胀系数(1/C) 8e-6 热传导率(W/m t) 0.52064 4.1.3边界条件定义 在PFC3D计算域内的加热过程由边界条件驱动，解决实际问题的过程 就是将边界数据推至计算域内部的过程。因此边界条件决定收敛速度与结果 的准确性。 物理模型采用温度和速度结合的混合边界。数值模拟中，设置滚筒筒 壁切向和法向刚度分别为5.958xlO6N/m和6?799乂10咖加，三段筒壁加热温度 分别为265C、260匸、150-C,倾斜角为2°，同时滚筒以匀速转动? 4.1.4载荷加载和求解条件设? 茶叶杀青过程数值模拟时步包括用于热计算的温度时步和用于机械计 算的机械时步，足够小的计算时步，可以保证数值模拟过程的稳定性，但时 步太小，将使得计算工作量大大增加。 由于茶叶杀青过程也是茶叶加热的过程，经时间山后茶叶温度7；⑷大于 0,由式(2-26)可知: (4-1) 其中：旳为杀青叶的质量，C,为茶叶的比热，-为连接热源的管道长度，Z为 管道单位长度的热阻。 由于热惯性远大于机械惯性，根据PFC3D要求，热计算时步和机械计算 时步应满足： 其中：4縉为热计算的时步，人匚祐为机械计算的时步。 在数值模拟过程中，取热计算最小时步7.0s3s,机械计算最小时步7e-6so 一般茶叶在杀青前含水率为72-75%左右，杀青后含水率为5840%,叶 温变化和失水率的关系可近似表示为线性关系【斶?】 2/=(-7.88x10) AT (4-3) 其中：AM为失水率，AT为叶温变化 4.2滾筒杀青机杀青过程温度场分析 为了在仿真过程中追踪颗粒的温度变化情况，给不同温度的颗粒标匕颜色, 每计算一个时步，颗粒的颜色更新一次，本论文用这种方法荻得颗粒的温度场分 布情况。颗粒温度犬于7(TC为红色，68C至70C为橙色，66匸至68C为黄色, 631至66￡为绿色，60*C至63￡为浅蓝色，低于60￡为蓝色卩Y 为了分析滚筒在不同结构和运动参数时是否均匀杀青，统计不同参数F茶叶