辅助充种带气吸式蚕豆精量排种器设计与试验.pdf

2023年月第5 4卷第7 期农报学业机械doi:10.6041/j.issn.1000-1298.2023.07.014辅助充种带气吸式蚕豆精量排种器设计与试验苏微赵庆辉赖庆辉谢观福田保宁王永杰（昆明理工大学现代农业工程学院，昆明6 5 0 5 0 0）摘要：针对蚕豆种子粒径大、三轴尺寸差异大，充种困难的问题，设计了一种带有平带辅助充种装置的气吸式蚕豆精量排种器。通过对充种过程中的动力学分析阐述了平带辅助充种装置及种子的运动机理；利用计算流体力学和离散元法双向耦合模拟的方法（Computationalfluiddynamics anddiscrete elementmethod,CFD-DEM），开展了单因素试验，确定了影响排种器充种性能的主要零部件参数并明晰了平带辅助充种机理；搭建试验台架，选取作业速度、平带输入轴转速和负压为试验因素，合格指数、重播指数、漏播指数为试验指标，进行了二次回归正交组合试验。试验结果表明，影响排种器合格指数的因素主次顺序为：作业速度、负压、平带输人轴转速。对试验结果进行多目标优化，得到最优参数组合为作业速度5.6 9km/h、平带输人轴转速3 95 r/min、负压3 8 45 Pa，对此结果进行排种器性能试验验证，此时合格指数为91.6%、重播指数为3.8%、漏播指数为4.6%，满足蚕豆播种要求关键词：蚕豆；排种器；耦合仿真；充种；气吸式；离散元中图分类号：S223.2文献标识码：A文章编号：10 0 0-12 98（2 0 2 3）0 7-0 144-12OSID：Design and Experiment of Air-suction Broad Bean Seed MeteringDevice with Flat Belt Auxiliary Seed-fillingSUWeiZHAOQinghuiLAI QinghuiXIE GuanfuTIAN BaoningWANG Yongjie(College of Modern Agricultural Engineering,Kunming University of Science and Technology,Kunming 650500,China)Abstract:An air-suction metering device for broad beans with a flat belt auxiliary seed-filling device wasdesigned to address the difficulties of sowing large-sized and highly different three-axis dimensional broadbean seeds.The motion mechanism of the flat belt auxiliary seed-filling device and seeds was elucidatedthrough a dynamic analysis of the seed-filling process.A one-factor test was conducted by using thecomputational fluid dynamics and discrete element method bidirectional coupling simulation method(CFD-DEM)to determine the main component parameters affecting the seed-filling performance andclarify the mechanism of the flat belt auxiliary seed-filling device.An experimental platform wasconstructed,and a quadratic regression orthogonal combination test was conducted with operating speed,flat belt input shaft speed,and negative pressure as the test factors and qualified index,reseeding index,and miss-seeding index as the test indicators.The test results showed that the primary and secondaryfactors affecting the qualified index of the air-suction metering device were operating speed,negativepressure,and flat belt input shaft speed.A multi-objective optimization of the test results yielded theoptimal parameter combination of the operating speed at 5.69 km/h,flat belt input shaft speed at395 r/min,and negative pressure at 3 845 Pa.The air-suction metering device performance was verifiedthrough sowing tests with a qualified index of 91.6%,reseeding index of 3.8%,and miss-seeding indexof 4.6%,which met the requirements for broad bean planting.Key words:broad bean;seed metering device;coupling simulation;seed-flling;air-suction;discreteelement method收稿日期：2 0 2 3-0 4-2 1修回日期：2 0 2 3-0 5-18基金项目：国家自然科学基金项目（5 2 16 5 0 3 1）作者简介：苏微（197 9一），女，副教授，主要从事农业机械化生产和管理研究，E-mail：l a i s u b o 16 3.c o m通信作者：赖庆辉（198 0），男，教授，主要从事农业机械装备与计算机测控研究，E-mail：l a i q i n g h u i 0 0 7 16 3.c o m145苏微等：辅助充种带气吸式蚕豆精量排种器设计与试验第7 期0引言蚕豆，可作为蔬菜、粮食、饲料、绿肥等，是一种重要的经济作物。我国蚕豆产量世界第一，但单产与英国等发达国家相比处于较低水平。其主要原因是我国蚕豆种植模式以人工点播为主，机械化程度低1-2 ，推进蚕豆播种机械发展是提高单产水平的重要手段之一。排种器是播种机的核心部件3-5 ，主要分为机械式和气力式两种6-7 。机械式排种器利用重力充种，对三轴尺寸差异大的种子适应性较差。气力式排种器具有通用性好、伤种率低、排种性能好及稳定等优点，是精量排种器的发展方向8 。蚕豆种子粒径大、三轴尺寸差异大，使用传统气吸圆盘式排种器结构存在充种困难、漏充率高的问题9。针对此类种子充种性能不佳的问题，国内外相关研究人员进行了一些研究：SATTI等【10-1 对影响气吸式玉米精量排种器各要素研究表明，排种器孔形、种盘转速以及真空压力是影响排种性能的主要因素；并在此基础上，进一步对排种器的重播、漏播等参数进行了试验、优化，对吸种、护种、投种等阶段进行动力学和运动学分析，最终得出排种器各项参数间的相互作用及影响关系。国内学者针对提高充填性能进行了大量的研究。赖庆辉等12 采用振动供种的方法，设计了一种振动供种装置，解决了微型薯的充种问题;LI等13 发现种子在吸附流中，吸附力对种子有显著影响的范围很小；史嵩等14 设计了一种带有曲线型导种槽的种盘，可以引导种子向型孔运动，改善了气吸式排种器的充种性能。陈进等15 利用电磁振动激振种群的方式降低种间作用力，达到了辅助充种的目的。谢东波等16 通过在种盘上增设扰种齿对大蒜种群进行扰动，提高了充种性能；刘瑞等17 通过设计扩容板增大充种区域并设计带有扰种台柱和中字型吸种孔的种盘，提高了充种性能；还有相关研究主要集中于对种群进行激励，降低种间作用力18-2 1;上述研究均可以一定程度提高排种器的充种性能，但大多是对种子进行无规则的激励，种群的运动不可控，存在种群对被吸附种子干扰的问题，针对上述问题，本文设计一种具有平带辅助充种装置的气吸式蚕豆精量排种器。通过对排种器充种过程进行理论分析，阐明平带辅助充种装置在充种过程中的作用机理；通过理论计算及CFD-DEM数值模拟确定并优化排种器主要结构参数；搭建台架进行二次回归正交组合试验，求得排种器的最佳工作参数组合，并进行验证。1排种器工作原理及平带搅种理论分析1.1排种器结构及工作原理辅助充种带气吸式蚕豆排种器结构如图1a所示。主要由种箱端壳体、种盘、护种板、主轴、种箱、清种装置、平带辅助充种装置及负压壳体构成，清种装置主要由锯齿形清种刀、调节柄构成。种盘根据工作状态分为充种区I、清种区、携种区与投种区IV。平带辅助充种装置位于充种区，由主动轮轴、平带、被动轮轴、张紧轴构成，通过一个步进电机带动同步带，如图1b所示。234569910111213148IV15(a)(b)图1平带辅助充种气吸式蚕豆排种器结构图Fig.1Structure diagrams of air-suction broad beanseed-metering device with flat belt auxiliary seed-filling1.种箱端壳体2.种盘3.护种板4.主轴由5.种箱6.清种装置7.平带辅助充种装置8.负压壳体9.同步带机构10.主动轮轴11.平带12.张紧轮轴13.平带辅助充种装置工作面14.被动轮轴15.步进电机工作时，种子由种箱落人充种区I的平带辅助充种装置上，被平带带动向种盘运动方向移动。众多种子在平带的带动下形成动态稳定的种子流。种子流上端保持在清种装置附近，种盘上的型孔靠近种子流组成的带状区域外侧，此处种子的运动方向与种盘转动方向一致，其中的1 2 颗种子会被型孔吸附，并随着种盘的吸附及种子流的托持下进入清种区。脱离种群后，由清种刀将多余的种子剔除，吸附最牢固的种子仍被吸附，被剔除的种子重新落回充种区中的种群，单粒种子在携种区由护种板护持至投种区IV，此时负压消失，种子在惯性和重力的作用下滑落，完成排种。1.2充种过程分析辅助充种带气吸式蚕豆精量排种器的充种过程尤为重要。如图2 所示，根据被吸附蚕豆种子的运动状态，可将充种过程分为3 个阶段：进种吸附阶段。托种跟随阶段。脱离种群阶段。平带辅助充种装置工作时根据种群的运动状态可分为：托种区、吸附区、回种区。进种吸附阶段：种子由种箱进人充种区初始端，在种群运动过程中，一些种子会靠近型孔；当种子达农机1462023年报学业械托种区吸附区回种区种盘转动方向种箱进种区托种区种子运动方向吸附区种子运动方向回种区种子运动方向主动轮轴转动方向图2充种过程分析Fig.2Analysis of seed-filling process到可以被气流吸附的临界位置时，种子会被种盘上的型孔吸附,最终被控制在型孔附近，并跟随种盘一起运动；在此阶段，需要对种群施加强制运动，使种群能够相对有序的运动，从而使种子更有效率地被吸附。从吸附流对种子的作用来看，种子自身向型孔附近运动是种子被更好地吸附的关键2 。托种跟随阶段：种子被吸附流稳定吸附开始跟随种盘运动，由于其运动与种群中其他种子的运动状态不同，开始出现种群对被吸附种子的阻力；此时托种区种子在平带的作用下向着种盘的运动方向移动并对被吸附的种子进行托持，回种区的种子向着种盘运动方向相反的方向移动，对被吸附的种子形成阻力。运动中回种区种子较为稀疏，托种区种子较为密集，因此托持力起到主要作用。在此阶段，托持力与种盘和平带的相对速度相关，此相对速度与被吸附种子的运动速度相适应是提高充种效率的关键因素。回种脱离阶段：被吸附的种子运动到平带末端，此时种子即将脱离种群；吸附区未被吸附的种子以及托持区即将脱离平带的种子在惯性、重力及种间作用力下横向移动进入回种区；被吸附种子会继续跟随种盘运动脱离种群进入清种区。在这3 个阶段中，被吸附的种子都在平带辅助充种装置的作用下降低了受到的种群阻力，同时将一部分种群阻力转换为有利于其继续跟随种盘运动的作用力。因此，为了解决蚕豆种子在气吸式排种器中充种困难的问题，本文将重点放在优化平带辅助充种装置的结构并确定其最佳工作参数。根据运动状态分析种群，托种区种子主要受到平带的摩擦力向着种盘旋转方向运动，并托持吸附区种子也向着此方向运动；当种子运动到平带的末端，托种区及吸附区未被吸附的种子将因为惯性及重力的作用进入回种区，最终回落到整个充种区的起始位置。平带辅助充种装置的工作状态主要与主动轮轴的转速及平带本身的材料属性相关。平带的运动一方面使种群处于运动状态，使种群活跃，增加了种子被吸附的概率；另一方面,托种区种子的运动可以托持在吸附区已经被吸附的种子向种盘运动方向移动，使得种群对被吸附种子的阻力作为有利于脱离种群的托持力。1.3种子受力分析通过以上分析可以看出，平带辅助充种装置对蚕豆种子的充种过程极为重要，因此通过受力分析、运动分析来确定平带辅助充种装置对种群运动的影响。可从充种区的3 个分区分别分析：设种子在角度上被吸附，如图3 所示，从被吸附的种子分析，观察经过型孔垂直于种盘半径的弦剖面可以得出种子能够被吸附并运移的临界条件为f=f,+mgsin+N,sing-N,sinaf=mR+mgsin+Nsing-Nsino(1)f=f.+f,-fx1f,=m(F,+N,coso+N,cosw)式中f吸附区种群对回种区种子的摩擦力m种子质量被吸附种子的总静摩擦力f.一托种区种群对种子的静摩擦力一f一一种盘对种子的静摩擦力f一一回种区种群对种子的静摩擦力N,一一托种区种群对种子的支持力N一回种区种群对种子的支持力一种盘角速度一R一一型孔与种盘中心距离种子与种盘的摩擦因数F气流对种子的作用力NNDraGV(a)吸附区被吸附种子(b)托种区、回种区种子受力情况受力情况图3 充种过程受力分析Fig.3Force analysis of seed-filling process种子受到的气流作用力包含电力、浮力、升力（M a g n u s 升力、Saffman升力）、Basset力、压力梯度力等多种力的作用2 3 。由于空气的密度、粘度较低且气吸式排种器常常工作在较高雷诺数下，种子受到的浮力、升力通常可以忽略，通常以曳力表征种子147苏微等：辅助充种带生气吸式蚕豆精量排种器设计与试验第7 期受到的吸附力2 4；颗粒受到的曳力分为摩擦阻力和压差阻力，现有的电力计算表达式通常是包含试验测得的阻力系数、迎风面积及流体与颗粒的相对速度的经验公式；Wenyu-Ergun曳力模型，综合了颗粒低体积分数和高体积分数2 5 ，适合排种器工况的计算。本文采用Wenyu-Ergun曳力模型描述种子受到的吸附力。电力计算式为F,=3PFC(00.850cp)FD20,ReTpPp(2）式中FD1Wenyu模型电力CD1Wenyu模型阻力系数Pf流体密度流体速度颗粒速度Tp颗粒等效半径Pp颗粒密度当前状态颗粒体积分数颗粒紧密堆积颗粒体积分数FD2Ergun模型电力Re-颗粒雷诺数对托种区及回种区的种子受力情况分析，托种区为其中回种区为2=N,+fui+Gh(5)IN,I=INhI式中f托种区种子受到的总静摩擦力f2托种区种子受到的种群摩擦力f种子与平带的最大静摩擦力平带线速度Na平带对托种区种子的支持力种子与平带的摩擦因数N2吸附区种群对托种区种子的支持力r平带工作面半径Jfa=fia+mgcose+N,cosy(3)Lf,=uaNd2Na=mgsino+N,siny+m(4)If.21=If,1LIN,I=IN,IGh回种区种子重力N一吸附区种群对回种区种子的支持力一从受力分析可以看出，力F，是被吸附种子能够稳定吸附在种盘上的关键因素，根据伯努利方程可知在相同的吸附姿态下负压直接决定来流速度，进而影响电力。此外，被吸附种子还受到托种区及回种区种子的种间作用；被吸附种子能够被种盘稳定运移的条件为：被吸附种子自身受到的吸附力而产生的对种盘的最大静摩擦力以及托种区种群对被吸附种子运动方向的托持力之和等于被吸附种子自身重力、运动前方种子的阻力以及回种区种群对被吸附种子的作用力之和。从式（3）（5）可以看出，平带的线速度越大,其对充种区种群的作用力越大。此外，种子存在滚动，可以描述为T;=-,F,R;W;(6)式中T:切向应力滚动摩擦因数R,一接触点到质心的距离接触点处的单位角速度F颗粒间法向作用力托种区种群对被吸附种子的作用力取决于二者的相对运动趋势，若托种区种群运动滞后于种盘的运动，则种群对被吸附种子的运移形成阻力，反之则形成有利于其运移的托持力。2关键零部件设计2.1型孔直径及数量气流流量一定时，型孔过大或型孔数量过多会导致负压降低，种子受到的吸附力不足。型孔过小单位面积上的负压较强，但气流作用面积较小，对粒径较大的种子吸附能力较弱。型孔数量过少在相同作业速度的情况下种盘转速过高，不利于充种。通常结合农业机械设计手册确定型孔直径及数量公式为2 6 d=(0.6 0.7)W(7)3npKzV=(8)50式中d-型孔直径W一蚕豆种子平均宽度V一作业速度np种盘转速K-一株距Z型孔数量通常蚕豆种子平均宽度为15.2 3 mm，株距为0.150.3m，设作业速度小于7 km/h，种盘转速小于5 0 r/min。根据以上分析确定型孔直径为10 mm，型孔数量为12 个。148农2023年机报学业械2.2种盘直径由于型孔数量较少，种盘直径主要考虑种群堆积数量及充种时长。蚕豆种子粒径较大，为了保证足够的种群堆积数量，提高吸附概率。综合农业机械设计手册，选择种盘直径为140 2 6 0 mm。考虑到蚕豆种子形状不规则且粒径较大，较大种盘直径有利于充种,选取种盘直径为D=240mm。2.3平带辅助充种结构平带辅助充种主要参数有平带材料、平带工作面弧长和平带工作面宽度。常用平带材料有PVC平带、尼龙编织带及橡胶黄绿片尼龙基带。对三者分别进行与蚕豆的摩擦试验，试验结果如表1所示。表1名各材料与种子的摩擦因数试验结果Tab.1Test results of friction coefficient betweenmaterials and seeds材料倾角/（）摩擦因数PVC130.231尼龙190.344橡胶290.554根据理论分析可知，种子与平带的摩擦因数越大,平带对种子的最大静摩擦力越大。从式（1）、（3）可知，平带对种子的最大静摩擦力越大，种子被吸附的临界条件中，被吸附种子受到的总摩擦力越易达到。因此选择橡胶黄绿片基带为平带材料。为了能充分与种群接触，平带底边为排种器圆周底端，顶边以理论计算及试验得到。由式（3）、（4）可知，平带工作面直接影响托种区种群运动，其在竖直面上对种群的影响范围取决于平带与种群的摩擦力，代人表1中平带材料的摩擦因数，得到f;=0.554mgsing+N,siny+m(9)7在平带工作面顶边，种群对平带的压力垂直作用于平带切向方向，即=90。设平带工作面转速等于实际作业最大速度，此时平带受到种群的作用最大。以单粒种子受到平带作用为准进行计算，当f为最大静摩擦力即fa=f,时，得到平带作用最大静摩擦力为0.50322nmgcos0=0.554mgsing+m(10)0.1200.5032T即Q=arcsin(11)0.120/2(ng+0.554g)4其中，n为该处种子数量。由于种群运动的不确定性，n难以确定，需进行数值模拟分析。3数值模拟利用计算流体力学和离散元法双向耦合仿真(Computational fluid dynamics and discrete elementmethod，C FD-D EM）进行耦合求解的方法进行数值模拟，该方法能够实现颗粒与流体的相互作用2 7 。求解器采用EDEM和Fluent,利用EDEM中的Bonding模型,将小颗粒粘结成一个实际种子，解决网格尺寸大于颗粒粒径的3 倍的计算要求。3.1仿真模型建立首先用粒径0.8 mm的球形颗粒将蚕豆种子填充，并导出颗粒坐标，颗粒建模如图4所示。利用SolidWorks建立排种器的简化模型，保存成ICES格式文件，导人EDEM中。流体域模型先由SolidWorks建立CAD模型，各部分按独立实体建立；利用ANSYS mesh工具进行网格划分，将划分的网格文件导入Fluent中，排种器简化模型及流体域网格模型如图5 所示；设置型孔部分为滑移网格，其转动中心为绕绝对坐标系原点指向Z轴正方向，转速根据试验水平分别设置；将型孔与种箱端及负压端的接触面设置为Interface，型孔部分接触面设置为与型孔相同的转速。(a）蚕豆种子(b）三维扫描(c）多球模型（d）元颗粒模型种子模型图4蚕豆种子颗粒建模Fig.4Broad bean seed particle model（a）排种器简化模型(b）流场网格划分图5气吸式排种器仿真模型Fig.5Simulation model of air-suction metering device现有研究中，通常先将粘接好的颗粒快速填充至充种区，再进行耦合仿真，这种方式不能很好地描述种子从种箱进入充种区的过程，不能维持种层高度，且初始填充的颗粒过多，极大地浪费了计算资源和存储空间。本文采用元颗粒（Meta一Particle）功能实现种子在仿真过程中的定量喂入，元颗粒可以与普通颗粒一样在颗粒工厂中生成,其本质是由一组坐标表示的颗粒团，并配合BondingV2模型，在颗粒团生成的瞬间同时生成Bonding键。通过预试验得出：设置2 个颗粒工厂，第1个设置为静态颗粒工厂，在充种区静态生成一定量的种子加快仿真进149苏微等：辅助充种带气吸式蚕豆精量排种器设计与试验第7 期程；第2 个设置为动态颗粒工厂，仿真过程中以一定速度生成种子,实现对种层高度的控制，使数值模拟更加符合实际工况，3.2仿真试验设计确定EDEM的时间步长为110-s,Fluent的时间步长为110-4s。仿真步数为8 0 0 0 0 步即仿真排种器实际工作时长8 s。利用Bonding模型描述种子时，种子的物理特征需要由Bonding模型中的切向刚度、法向刚度表征。仿真参数如表2 所示。表2 仿真参数Tab.2Parameters of simulation模型参数数值泊松比0.25蚕豆种子（颗粒）杨氏模量E/Pa1.059 108密度p/(kgm3)1645法向刚度/(Nm-3)5 1010颗粒-颗粒(bonding)切向刚度/(Nm=3)8109碰撞恢复系数0.402颗粒一颗粒静摩擦因数0.365滚动摩擦因数0.05碰撞恢复系数0.598颗粒一种盘静摩擦因数0.457滚动摩擦因数0.03碰撞恢复系数0.329颗粒一平带静摩擦因数0.554滚动摩擦因数0.08为保证工作面设计合理，仿真以工作面宽度为因素，设置平带工作面宽度分别为2 5、3 0、3 5、40、45mm进行单因素试验。在EDEM中对平带工作面添加ConveyorRotation运动，根据理论分析及预试验观察，分别设置作业速度为5.5 km/h、平带工作面转速为40 r/min，负压5 0 0 0 Pa。并选取其中被吸附种子的受力情况进行分析，以确定平带宽度的最佳参数。CFD-DEM仿真计算量大、耗时较长、占用存储空间较大，因此在保证试验结果可靠的情况下尽量减小仿真时间，通常选取10 0 2 0 0 颗种子作为统计数据。由于辅助充种带气吸式蚕豆排种器型孔数量较少，作业速度较慢，确定选取8 0 粒种子进行统计。计算完成后，通过EDEM后处理获取仿真结果，如图6 所示。3.3仿真结果分析从图6 可以看出，在平带工作面上端，种子的聚集量通常为3 4颗，即n最大为4，代人式（11）得0=2.452046，即平带辅助充种装置工作面弧长范围为2.45 2 0 46 0 90。过大的弧度会使得平带与排种器摩擦力增大，且不易加工安装，为了方便，图6EDEM后处理界面Fig.6EDEM post-processing取5 85。即工作面弧长l为16 7.47 mm。标记仿真过程中吸附区被吸附种子。选取吸附较好，稳定且姿态尽量一致的种子，同时标记平带工作面末端托种区种子，分别提取各水平种子的小颗粒受到的总力，如图7 所示。由式（1)可知，被吸附种子所受外力的量和模数越小，颗粒受到除吸附力以外其他外力影响越小。综合分析各水平总力随时间变化可知，种子在初始状态受到外力最大，此时种子处于种群运动中，受到种群影响大；被吸附后，种子进人吸附区，受到总力逐渐稳定；脱离种群瞬间，种子受到回流种子影响，总力出现跳变，接着进人清种区，出现二次跳变；种子进人携种区后总力恢复稳定，直至投种总力归零。由图7 可知，随着平带宽度的增加总力先增大后减小，总力波动先增大后减小。观察试验结果发现，随着平带宽度的增加，平带宽度方向上容纳种子数量变化为：容纳1颗种子并存有较大缝隙、恰好容纳2 颗种子、容纳2 颗种子并存在较大缝隙、恰好容纳3 颗种子。即当平带宽度略大于2 颗种子宽度之和时，托种区种子运动较为规律，对被吸附种子的影响较小，而在平带宽度方向上并排种子数量大于3颗或小于2 颗时，种群运动更加不规律，导致被吸附种子受到托种区种群外力冲击增大。综合考虑，设计平带宽度为3 5 mm。4台架试验4.1试验材料与设备根据理论计算及数值模拟试验研制平带辅助充种气吸式蚕豆排种器。采用云豆早7 蚕豆种子作为试验对象。在JPS-12型计算机视觉排种器性能试验台上进行试验，如图8 所示。通过合肥富煌君达高科信息技术有限公司生产的千眼狼5 F01型高速摄像机对排种器工作过程进行分析。根据CB/T6973一2 0 0 5 单粒（精密）播种机试验方法，在排种150农2023年机业报学械0.250.250.25,0.200.200.200.150.150.15N/YN/4N/0.100.100.100.050.050.0500.51.01.52.000.51.01.52.000.51.01.52.0仿真时间/s仿真时间/s仿真时间/s(a)b=25mm(b)b=30 mm(c)b=35 mm0.250.25,0.200.200.150.15N/9N/40.100.100.050.0500.51.01.52.000.51.01.52.0仿真时间/s仿真时间/s(d)b=40 mm(e)b=45mm图7仿真试验结果Fig.7Simulationtestresults3476图8 排种器性能试验装置Fig.8Seed-metering device performance test1.计算机2.排种器试验机操作台3.辅助充种带气吸式蚕豆排种器4.排种器主轴驱动电机5.高速摄像补光灯6.步进电机驱动器及单片机7.高速摄像机器稳定工作时，连续记录2 5 0 粒种子作为统计对象，每次试验重复5 次取平均值。4.2单因素试验4.2.1负压根据预试验确定负压分别为2 0 0 0、3 0 0 0、40005000、6 0 0 0 Pa。设置作业速度为5 km/h、平带输人轴转速为3 5 0 r/min。试验结果如图9所示。试验结果表明，随着负压增大，合格指数先增大后减小、重播指数逐渐增大、漏播指数逐渐减小。根据理论分析可得，随着负压增大，吸附力增大，种子吸附更牢固不易清种，因此重播指数上升，漏播指数下降。由于负压在3 0 0 0 5 0 0 0 Pa之间合格指数均大于8 5%，重播指数和漏播指数均为较低值，因此确定排种器负压范围在3 0 0 0 5 0 0 0 Pa之间。100-合格指数重播指数漏播指数80%/420020003000400050006000负压/Pa图9不同负压下试验结果Fig.9Test results under different negative pressures4.2.2作业速度根据预试验确定作业速度分别为4.3 2、5.0 4、5.76、6.48、7.2 0 k m/h。负压为40 0 0 Pa、平带输人轴转速为3 5 0 r/min。试验结果如图10 所示。100r一一合格指数重播指数漏播指数80%42004.325.045.766.487.20作业速度/(kmh-)图10不同作业速度下试验结果Fig.10Test results at different operating speeds试验结果表明，随着作业速度的增加，合格指数151苏微等：辅助充种带气吸式蚕豆精量排种器设计与试验第7 期先上升后下降，重播指数逐渐下降、漏播指数逐渐上升。当作业速度增大，种盘转速增大，离心力增大、充种时间减少，种子被吸附的概率降低；被吸附种子的运动速度增大导致其动量增大与种群碰撞时，受到的反作用力增大。因此随着作业速度的增加,重播指数下降、漏播指数上升。由于作业速度为5.0 46.48km/h时，合格指数均在8 5%以上，确定作业速度范围在5.0 4 6.48 km/h之间。4.2.3平带输入轴转速根据预试验确定平带输人轴转速分别为0、200、40 0、6 0 0、8 0 0 r/m in。作业速度为5.7 6 km/h、负压为40 0 0 Pa。讨试验结果如图11所示。100合格指数重播指数漏播指数%/802000200400600800平带输入轴转速/（rmin-）图11不同平带输人轴转速下试验结果Fig.11Test results of flat belt input shaft speed由图11可知，随着平带输入轴转速升高，合格指数先上升后下降、重播指数缓慢上升、漏播指数先下降后上升。当平带输人轴转速较低时，对种群的激励作用小，种间作用力大；随着平带输人轴转速的上升，种群“活跃”，种间作用力降低，且托种区种群对被吸附种子开始起到托持作用，漏播指数下降；随着平带输入轴转速的进一步提高，种群运动速度过快，托种区种群对被吸附种子的正向作用力过大，导致被吸附种子脱落，漏播指数上升。但随着平带输人轴转速的提高，种群“活跃性”始终上升，因此，重播指数一直缓慢上升。由于平带输入轴转速在252505r/min时，合格指数均在8 5%以上，确定平带输人轴转速范围在2 5 2 5 0 5 r/min之间。4.3二次回归正交组合试验4.3.1试验设计为了探究排种器的最佳工作数组合，进行二次回归正交组合试验。选取作业速度1、平带输入轴转速xz和负压x作为试验因素。合格指数Y、重播指数Y,和漏播指数Y，作为试验指标。试验因素和编码见表3,试验方案与结果见表4。X,、X2、X,为因素编码值。4.3.2回归方程建立及显著性检验利用Design-Expert12.0.1.0软件对试验结果表3试验因素编码Tab.3Experimental factors and codes因素编码作业速度平带输人轴转速负压x)/(kmh-1)2/(rmin-)x3/Pa-1.6825.042523 000-15.333033.40505.76378400016.1945445951.6826.485055000表4试验方案与结果Tab.4Text design and results因素序号Y/%Y2/%Y,/%XX2X31-1-1-167.610.122.321-1-146.96.946.23-11-174.013.512.5411-179.02.818.25-1-1176.322.51.261-1159.715.225.17-11153.638.28.2811158.427.414.291.6820078.518.62.9101.6820064.86.528.7110-1.682059.78.232.11201.682075.220.34.51300-1.68272.51.825.714001.68254.445.20.41500096.02.21.81600087.36.26.51700093.93.13.01800092.20.57.31900094.15.50.42000094.04.21.82100088.48.13.52200091.70.28.12300094.43.12.5进行回归分析，确定各因素的影响，得到合格指数Y、重播指数Y,和漏播指数Y,的回归方程。并进行回归方程的显著性分析，分析结果如表5 所示(1)合格指数Y,回归模型与显著性检验由表5 可知,该模型拟合度极显著（P0.01）。其中作业速度与负压交互项（X,X，)P大于0.0 5,对排种器工作时的合格指数不显著，其余项皆显著。说明相关试验因素对响应值的影响存在二次关系。失拟项P=0.2703不显著，说明不存在其他影响指标的主要因素。去除不显著的因素后进行多元拟合，建立合格指数Y的回归模型为Y,=92.48-3.7X,+2.97X,-3.66X,+5.89X,X2-7.81X,X,-7.68X,-9.16X2-10.58X3(12)农1522023年机报学业械表5 合格指数、重播指数与漏播指数的方差分析Tab.5Variance analysis of qualified index,reseeding index and miss-seeding index方差合格指数重播指数漏播指数来源平方和自由度FP平方和自由度FP平方和自由度FP模型5.251.25956.910.000 1*3085.37938.710.000 1*3311.65944.230.000 1*X187.04118.240.000 9*200.67122.660.000 4*775.17193.180.000 1*X2120.51111.750.004.5*165.56118.690.0008*568.56168.350.000 1*X3182.62117.810.001*1497.131169.060.000 1*633.98176.210.0001*X,X2277.3127.050.0002*15.1211.710.2139162.90119.580.0007*X,X31.910.18540.673 82.2010.2490.62610.011310.001 40.971 2X2X3488.28147.620.0001*102.24111.550.0048*143.65117.270.001 1*X936.43191.330.0001*146.47116.540.0013*342.21141.140.0001*X1 333.731130.080.000 1*210.20123.740.000 3*484.97158.30.0001*X1.777.261173.340.000 1*758.14185.610.000 1*213.84125.710.000 2*残差133.2913115.1313108.1413失拟66.1151.570.270 360.5751.780.223948.1951.290.3573误差67.18854.56859.968总和5384.54223200.5223419.822注：*表示影响显著（0.0 1P0.05），*表示影响极显著（P0.01）。根据式(12)可以得出，影响合格指数的因素主次顺序为作业速度、负压、平带输入轴转速。(2)重播指数Y,回归模型与显著性检验由表5 可知，该模型拟合度极显著（P0.01）。其中作业速度与平带输入轴转速交互项（X,X,）、作业速度与