免耕播种机抛送装置数值模拟与试验研究.docVIP

免耕播种机抛送装置数值模拟与试验研究 　　摘要：为了了解全秸秆覆盖地免耕播种机叶片式抛送装置气流流场的分布情况以及抛料管道气流速度的影响因素，采用FLUENT对抛送装置进行数值模拟并对其进行了试验研究。结果表明，气流速度沿叶轮径向方向由内向外逐渐升高；圆形壳体出料口外侧的气流速度较内侧处的高；模拟所得结果和试验值的分布趋势基本一致，验证了数值模拟的可靠性；抛送叶轮的转速越高、叶片直径越大，各测点的气流速度就越大，越有利于秸秆物料的抛送；当叶片数为4、叶片倾角为后倾10°时，更有利于秸秆物料的抛送。 　　关键词：全秸秆覆盖地免耕播种机；抛送装置；数值模拟；试验研究 　　中图分类号： S223.2文献标志码： 　　文章编号：1002-1302（2016）08-0410-04 　　近些年来，我国面临着较为严重的空气污染问题，每到秋末冬初季节，全国多个省市都会被持续的雾霾天气笼罩着，这对人们的生产生活以及身体健康造成了极大的影响。农田废弃秸秆的焚烧是造成这种问题的主要原因之一。由于前茬农作物收获后，会有大量的碎秸、残茬留在地表，对现行的免耕播种机而言，大量的碎秸、残茬会引发机具入土部件挂草壅堵，进而导致架种、晾种等，这不仅影响作业质量而且影响作业顺畅性[1]，所以就地焚烧秸秆成了农民“最省事”的选择。国家和农民都迫切需要能够满足工况的机具。农业部南京农业机械化研究所正在试验研制的全秸秆覆盖地免耕播种机，能够很好地解决现有的问题。该机具主要作业工艺：秸秆经过秸秆粉碎装置粉碎后，在粉碎装置高速旋转产生的气流带动下进入横向输送搅龙，然后被推送至叶片式抛送装置，越过种箱、肥箱，均匀地抛撒于播后区域[2]。抛送装置是该机的一个重要环节，其高速旋转的抛送叶轮产生的气流是影响碎秸能否顺利抛撒以及功率消耗的重要因素。因此，本研究对抛送装置进行数值模拟以及试验研究，了解其基本特性。 　　1理论计算方法 　　1.1数学模型 　　由于叶片式抛送装置工作?r，气体在抛送叶轮内部的流动是非定常流，因此，采用MRF（多参考系模型）作为定常流动近似模拟此旋转流场[3]。当采用标准的k-ε模型求解时，控制方程包括连续性方程（质量守恒方程）、动量守恒方程以及k方程和ε方程，方程如 　　采用有限体积法离散控制方程，对该抛送装置的流场进行稳态隐式、非耦合求解。考虑到壁面附近区域受到分子黏性的影响，运用标准壁面函数法模拟近壁区域的流动、动量、耗散率和湍动能的离散格式均设置为精度较高的二阶迎风格式，对压力―速度耦合运用SIMPLE算法求解[7～9]。 　　2抛送装置气流流场的试验研究 　　为了了解叶片式抛送装置抛料管道处气流速度的大小与抛送叶轮转速、叶轮外径、叶片数目以及叶片倾角之间的关系，同时也为了验证数值模拟结果的可靠性，对抛送装置的气流流场进行试验研究。 　　2.1试验设备与仪器 　　试验是在自制的试验台上进行（图2）。测量时没有物料的喂入。气流速度测试采用的仪器为JX-2000皮托管风速仪（上海金枭电子）。 　　2.2试验方案的选择 　　抛送叶轮的转速以及叶片的结构形式是影响抛送装置气 　　2.3试验方法 　　选取抛料管道处的A截面为测试位置，如图1-a所示。共选择12个点作为测量点，如图1-b所示，在A截面中D上取4个点，E线上取4个点，F线上取4个点。测量风速时，需在A截面侧面打孔进行测量。每个测点重复测试3次，取3次测试结果的平均值为最终结果。 　　3计算结果与分析 　　3.1模拟结果与分析 　　由于FLUENT中的可视化信息基本是以平面为基础，所以，当数值模拟计算区域后，为了显示以及输出计算结果，需生成各种类型的平面。图3-a、图3-b、图3-c分别为Z=-0.07、0、0.07 m的速度矢量图。由图3可知，气流速度沿叶轮径向方向由内向外逐渐升高，这主要是由高速旋转的抛送叶轮使气体分子在离心力的作用下获得动能所造成的。另外，圆形壳体出料口处外侧的气流速度较内侧处的高，这主要是因为高速旋转的叶片末端靠近外侧。 　　从图3-c即Z=0.07 m平面处的速度矢量图中可以观察到较为明显的涡流，这是由于其靠近进料口处以及叶片数有限。进口处的气流沿着轴向方向进入叶轮后，在抛送叶轮的作用下改变90°方向沿着叶轮径向方向流动。从 Z=-0.07 m 和Z=0 m平面的速度矢量图也可以看出少量的涡流，造成这种现象的主要原因是叶片数目有限。所以，可以通过增加叶片数量来减少叶轮区的涡流，进而减少能量损失，得到更有利的秸秆物料抛送的气流流场。 　　[FK（W41][TPHZD3.tif][FK）] 　　3.2模拟结果与试验值的比较分析 　　当转速为1 600、2 000 r/min时，抛送装置中A截面上D、E、F线的速度分布散点图如图