近年来，云南省蔬菜产业快速发展，云南已成“南菜北运”和粤港澳大湾区的重要“菜篮子”基地^[1]。2021年云南省蔬菜产量超过2753.51万t，同比增长9.8%^[2]。然而蔬菜育苗是一项劳动力密集、劳动强度大的工作^[3]。究其原因，与缺乏专业化小型机械作业，基础设施不足等因素相关^[4]。工厂化育苗播种能够提升蔬菜育苗效率和经济效益^[5]。同时采取新型的管理技术与新型穴盘育苗技术结合，可以大大缩减人力物力^[6]。但对于精量播种的蔬菜播种机研发，必须紧密结合农机农艺相融合的特点^[7]。我们在用蔬菜穴盘育苗播种机实施播种时发现，该机因充种过程中会产生漏充及多充现象导致充种合格率低，尤其多充现象明显，但多充可采用适宜的清种装置避免这一情况的发生。因此设计一款合适的清种装置对于提高蔬菜穴盘育苗播种机播种精度和效率显得尤为重要。

精密播种的清种装置主要有机械式、气吸式和电控式3种。王计新^[8]进行了斜面重力的清种型排种器研究，运用正交试验得出了最优参数组合；许建等^[9]通过对圆盘式排种器建立清种过程数学模型，得出排种盘半径、垂直倾角和排种盘转速是影响排种器清种能力的主要因素，为机械式排种器关键部件的设计提供了理论指导；Braunack等^[10]研究得出内充种式排种器具有结构简单、造价低、使用维护方便等优势，同时配合所需清种装置可适用于中耕作物的精播作业；刘云强等^[11]为解决蔬菜穴盘育苗精密播种机播种非球形种子时重播率较高的问题，设计了一种清种装置，并通过对清种装置内流场仿真模拟优化了清种装置结构参数，在吸种气压为15.7 kPa，清种气压为3.3 kPa，振动频率为50 Hz时效果最好；康施为等^[12]考虑棉花排种器易堵塞等问题，设计了排种器及清种结构，并借助气流进行清种，保证单穴单粒的播种合格率；Yazgi等^[13]利用响应曲面法对气吹式排种器清种性能进行研究，提供了一种优化精密播种机性能的方法，所建立的方程可用于估计播种机工作性能相关变量，通过试验验证并优化了参数；贾洪雷等^[14]设计了气吸机械复合式大豆精密排种器，并得出当导种槽倾角为 45°、深度为 2 mm时，可最大程度提高种群上升区的种子运动速度，增加种群离散程度，进而减小种群内摩擦力，提高充种效率；杨铭^[15]通过图像处理及电控装置对排种器充种过程实现振动式清种，同时也可有效清理孔眼堵塞现象。

通过对国内外研究现状分析可知，以上设计的排种装置对清种过程的理论分析较少，而且大多是针对单级清种装置，存在清种不彻底、清种效果还不够好等问题。为了优选出提高播种合格率、降低种植成本的相应参数，本文针对蔬菜穴盘育苗播种机，设计机械式、气吹式、电控式三级调节清种装置，并开展四因素三水平二次正交旋转组合台架试验，以得出播种的最佳工作参数，从而为实际生产提供理论指导。

1.   三级清种装置结构和工作原理

因机械式清种装置便于根据种子三轴尺寸调节，且考虑到团簇形种子在清落过程中极易出现前后及左右2粒种子同时存在的情况，设计了机械式清种装置用于第一级清种；因左右2粒种子在吸附过程中，易因气流场变化而清落其中受吸附力较小的1粒种子，故设计了气吹式清种装置；因前后2粒种子在运动过程中有圆周切向的种子位置偏差，无法由机械式及气吹式清种装置实现清落，但可通过步进电机经传动装置控制滑轴，清落前后2粒种子中的后排种子，故设计了电控式清种装置；且电控式清种会对左右2粒种子因水平交错偏差产生错误清种的现象，故第二级为气吹式清种装置，第三级为电控式清种装置，结构如图1所示。

图  1  清种装置结构示意图

1：机架；2：滤轴；3：滑轨；4：气嘴；5：凸台；6：凸轮；7：齿轮；8：步进电机

Figure  1.  Schematic diagram of seed cleaning device

1: Frame; 2: Filter shaft; 3: Slide; 4: Air nozzle; 5: Boss; 6: Cam; 7: Gear; 8: Stepper motor

下载: 全尺寸图片 幻灯片

三级清种装置的工作原理：采用负压力充种后，进入清种阶段，清种阶段又分为3个过程。首先是机械式清种阶段，滤轴通过螺栓安装在机架上并可旋转，通过滤轴中的滤孔对簇形种子在清落过程中极易形成的前后或左右2粒种子进行过滤清种；然后是气吸式清种阶段，由两端圆形凸台固接在机架，通过气流和排种器吸附左右2粒种子的溢出气流形成对冲，以清除左右2粒种子；最后是电控式清种阶段，通过步进电机带动齿轮发生转动，使得偏心凸轮带动滑轴在滑轨上根据排种器周向的2排种子的时间差实现往复运动，从而清落前后2粒种子的其中1粒。

2.   关键部件设计与结构参数确定

2.1   蔬菜种子的机械物理特性测定

种子物料特性是研究的重要依据^[16]，根据蔬菜穴盘育苗播种样机对蔬菜种子的适应性，选取似圆度较高的甘蓝、花椰菜、娃娃菜等包衣蔬菜种子为研究对象，因不同包衣种子大小不同，测量各种子的三轴尺寸及千粒质量，记录数据取均值，测算包衣种子平均长度为2.12 mm、平均宽度1.93 mm、平均高度1.84 mm、平均千粒质量165.86 g，经式(1)得出包衣种子球度为92.45%。可近似将上述包衣种子看作球体来计算^[17]。

式中，S_p为种子平均球度，L为种子平均长度，D为种子平均宽度，H为种子平均高度。

2.2   机械式清种装置

滤轴是清除团簇形种子的关键部件，其与排种器间隙配合，故滤轴的长度可由排种器结构尺寸确定；为了保证单粒种子通过并阻隔多粒种群，因此滤孔的大小需满足农艺要求；滤轴直径与滤轴承力有关，为保证滤轴在与排种器接触过程中有足够的疲劳强度，滤孔的边界不能超过滤轴中心线位置，则滤轴尺寸参数应满足：

式中，L_p为排种器轴向长度，L_h为配气盘外壁厚度，L_c为固接齿轮厚度，L_J为固定补偿，L_L为滤轴总长，L_min为不同种子三轴尺寸最小值，L_max为不同种子三轴尺寸最大值，L_k为滤孔直径大小，D为滤轴直径，L_{k max}为滤孔直径最大值。

根据式(2)，测量排种器参数可得滤轴长度为290 mm；由种子三轴尺寸可得滤孔直径范围为2.13~3.38 mm，考虑排种器孔径大小及为提高清种装置的适应性，确定滤孔直径范围为1.50~5.00 mm；因滤轴直径过大将会造成安装位移偏差过大，故由上式可取滤轴直径为10 mm。最终滤轴结构模型如图2所示。

图  2  滤轴结构模型图

L_p：排种器轴向长度；L_h：配气盘外壁厚度；L_c：固接齿轮厚度；L_J：固定补偿；L_L：滤轴总长

Figure  2.  Structural model of filter shaft

L_p: Axial length of seeding device; L_h: Outer wall thickness of valve disc; L_c: Thickness of solid gear; L_J: Fixed compensation; L_L: Total length of filter shaft

下载: 全尺寸图片 幻灯片

2.3   气吹式清种装置

2.3.1   起点处运动学分析

因气吹式清种装置出气口气流受排种器的溢出气流影响，故对排种器吸孔处的气流进行分析，因排种器吸孔接触种箱准备吸附时的一瞬间，阻抗力是最大的，故对此时的受力状态进行平衡分析(图3)。

图  3  种子刚吸附时受力平衡图

F_b1：种箱对种子的支持力；F_b2：排种器对种子的支持力；F_X：负压吸附力；F_n：种群力；α：休止角；G：单粒种子重力

Figure  3.  Force balance map of the seed just adsorbed

F_b1: Supporting force of the seed box on the seed; F_b2: Supporting force of the seeding device on the seed; F_X: Negative pressure adsorption force; F_n: Seed group force; α: Angle of repose; G: Gravity of a single seed

下载: 全尺寸图片 幻灯片

式中，F_b2为排种器对种子的支持力，F_X为负压吸附力，F_n为种群力，α为休止角，F_b1为种箱对种子的支持力，G为单粒种子的重力。

由上式可知，若种子附着于排种器转动，则：

种子吸附时，会存在机器振动及大气压强，为增加吸附可靠性，故考虑0.2倍气压压力损耗，又因排种器负压过大会造成种子大量破损，故：

因实际种子受力接触流通截面为半球面，考虑不同种子的接触流通截面相差不大，简算为孔径流通截面(A)^[18-19]，则：

式中，A为实际流通截面积；r为吸孔半径；β为气流场中心线夹角；V_R为相对于吸孔外气流场的气体流速，m/s；Q_v为吸孔内部体积流量，m³/s。

因吸种孔孔径很小，根据伯努利定理及气体流动状态可推断出孔径内部气体流速V_C：

式中，V_C为孔径内部气体流速，m/s；K为绝热指数；R_T为气流常数，J/(kg·K)；T_e为空气隔热温度，K；P_i为排种器气流压力，Pa；P_e为大气压强，Pa。

流经单个吸孔的体积流量为单个吸孔的截面面积与吸孔气体流速的乘积。流经单个吸孔的体积流量则为：

由上述理论条件下计算可得，单个内孔压强≥5.92 Pa，故排种器整体气压需在959.04 Pa以上。

因蔬菜种子根据形状可分为圆形种子和非圆形种子，本文主要针对包衣后的似圆形种子进行动力学分析，定义沿排种器轴向靠左种子为C₁，靠右种子为C_2。

因左右2粒种子同时被吸附，其间存在孔眼间隙，必定会减小充种时的理论气压，故增加${\xi _1}$，由图4a受力平衡可得：

图  4  气吹式清种区左(C₁)右(C₂)种子受力图

f_a2：C₁种子抵抗C₂种子作用力后与排种器间的摩擦力；F_n3：C₁种子对C₂种子的作用力；γ₂：气吹式清种区起点与水平面的夹角；δ₂：C₁种子中心与排种器轴线的垂直连线和排种器气孔径间的夹角；F_b4：C₁种子的理论支持力；f_a1：C₂种子抵抗C₁种子作用力后与排种器间的摩擦力；F_n2：C₂种子对C₁种子的作用力；δ₁：C₂种子中心与排种器轴线的垂直连线和排种器气孔径间的夹角；F_b5：C₂种子的理论支持力；F_X：负压吸附力

Figure  4.  Stress map of left (C₁) and right (C₂) seeds in air blowing seed cleaning area

f_a2: Friction force between C₁ seed and seeding device after resisting the force of C₂ seed; F_n3: Force of C₁ seed on C₂ seed; γ₂: Angle between the starting point of the gas-blown seed cleaning area and the horizontal plane; δ₂: Angle between the vertical line connecting the C₁ seed center and the axis of the seeding device and the gas aperture of the seeding device; F_b4: Theoretical support of C₁ seed; f_a1: Friction force between C₂ seed and seeding device after resisting C₁ seed force; F_n2: Force of C₂ seed on C₁ seed; δ₁: Angle between the vertical line connecting C₂ seed center and the axis of the seeding device and the gas aperture of the seeding device; F_b5: Theoretical support of C₂ seed; F_X: Negative pressure adsorption force

下载: 全尺寸图片 幻灯片

式中，f_a2为C₁种子抵抗C₂种子作用力后与排种器间的摩擦力，F_n3为C₁种子对C₂种子的作用力，${\gamma _2}$为气吹式清种区起点与水平面的夹角，${\delta _2}$为C₁种子中心与排种器轴线的垂直连线和排种器气孔径间的夹角，因转换成平面夹角计算，故添加关系系数${\xi _2}$，F_b4为C₁种子的理论支持力。

由图4b受力平衡可得：

式中，f_a1为C₂种子抵抗C₁种子作用力后与排种器间的摩擦力，F_n2为C₂种子对C₁种子的作用力，${\delta _1}$为C₂种子中心与排种器轴线的垂直连线和排种器气孔径间的夹角，因转换成平面夹角计算，故添加关系系数${\xi _3}$，F_b5为C₂种子的理论支持力。

通过考虑了排种器负压系数${\xi _1}$，结合分析F_n2和F_n3，排除了同时将双粒种子清落的情况发生。

2.3.2   气吹式清种装置的结构参数确定

气流在喷嘴出口处会产生射流现象，在大气阻力下，气体流速会随着距离的增大而变小^[8]。因此种子越靠近喷嘴气流出口处，种子受到的气流越集中，故种子与喷嘴气流出口之间的距离(L₁)应尽可能小，取L₁=10 mm, 为了达到稳压的目的，需要在气压口处设置圆弧式喷嘴，但却增加了气流路径及气流管道内的损耗，导致气流场减弱，而气流场的强弱可由流入气压腔的气流体积流量大小决定，由此推断需补偿圆弧式喷嘴出气口气压，使气流场增强，因此需要进气口气流体积流量的增大，故取进气口直径D₁=8 mm。

根据气流场流量分布，可得：

式中，H为气流体积流量比例系数，Q_v1为清种装置进气口流入的气流体积流量，${\beta _1}$为进气口处气流场锥角。

因流经单个吸孔的体积流量为单个吸孔的截面面积与吸孔气体流速的乘积，故可得：

因此气流方程可表示为：

式中，R为气流场任意点到气流场中心的距离。

由式(14)可知，气流场中任意质点的速度与气流场中心距离的平方成反比，因此可通过缩短气管的长度来保证气压的稳定。为使气压腔内从每一个圆弧式喷嘴流出的气压尽可能保持一致，由伯努利定理得：

式中：P为流体中某点的压强，v为流体该点的流速，ρ为流体密度，g为重力加速度，h为该点所在高度，C是一个常量。由式(15)可知，在同一等高线上，气体流速越大，气压越小。

因测量的种子平均千粒质量为165.86 g，排种器单个内孔最低压强为5.92 Pa，气吹式清种装置中种子与喷嘴气流出口之间的距离为10 mm，进气口直径为8 mm，同时假定出气口直径为2 mm，取${\xi _1}$为0.8，${\xi _2}$、${\xi _3}$为0.95，得进气口压强P_j为0.3 kPa。吸附的种子在排种器运动过程中会减小与气吹式清种装置出气口的距离，使得种子随排种器转动角度更容易被清落，为避免清落双粒种子的可能性，在式(10)、(11)中已考虑了排种器负压系数${\xi _1}$，故排除了同时将双粒种子清落的情况发生。

2.4   电控式清种装置

2.4.1   电控式清种装置的结构参数确定

电控式清种装置的结构如图5所示，因选用步进电机转轴为8 mm，故加工时其中主动锥齿轮的孔径需少于50丝与转轴进行过盈配合，从动锥齿轮的转轴一端则需要连接偏心凸轮。凸轮的偏心距需配合滑轴的运动轨迹位移来实现，而滑轴的运动轨迹位移则需根据排种器轴向相邻两孔径中心轴间距及种子最大三轴尺寸直径来确定，为防止出现滑轴弯轴偏差及似圆形种子半径偏差，增加种子三轴尺寸最大半径补偿。因此可知满足要求的偏心凸轮行程为3.2 mm，故滑轴的移动距离为3.2 mm。滑轴上装有拨片，根据实际考察及结合农机农艺相融合的要求，确定排种器转速为0.11~0.17 r/s，那么排种器转1圈最长需9.09 s，因此滑轴的速度须达到0.29 mm/s以上。为能够安装滑轴，抵住弹簧压缩，且便捷拨片运动，因此滑轨孔径设置为7 mm，面向排种器的平面需开2 mm的拨片沟槽。

图  5  电控式清种装置结构示意图

Figure  5.  Schematic diagram of the structure of the electronically controlled seed cleaning device

下载: 全尺寸图片 幻灯片

2.4.2   电控式清种装置的运动设计

偏心凸轮旋转180°的清种过程需在0.091 s内完成，考虑偏心凸轮的恢复过程会受到下一排种子存在吸附前后双粒种子状态的影响，根据排种器的结构设计，以及为简算程序设计，则固定回程为0.091 s，且在4个回程后可对后排种子再次进行清种。本次设计选择三菱57步进电机，采用GX Works2软件进行程序读写，在机器与电控装置同时启动的情况下，实现步进电机经传动部分对滑轴无间隔清种。

3.   气腔结构Ansys仿真分析

3.1   Ansys仿真分析

气吹式清种装置气腔进气口与出气口结构对气流场稳定性有极大影响，故通过软件Ansys对气腔内气流场的气流分布、气体流速及气压大小进行仿真分析，以确定气吹式清种装置气腔整体结构。

借助Solidworks软件将4组不同形式的气吹式清种装置三维模型保存为x_t文件，进而导入Ansys软件的ICEM CFD模块中，对气流入口处及出口处设置四面体网格为0.6 mm，其余位置设置最大网格尺寸为2 mm的网格划分，同时使用修复手段将网格质量低于0.4的网格删除，以提高划分网格的质量。进而求解网格模型，导出Mesh网格文件。把已生成的网格文件导入Fluent软件中，设置气流进口处采用Pressure-inlet边界，气流大小由进口流量和结构尺寸确定为0.3 kPa，气流出口处采用Pressure-outlet边界，气腔表面采用Wall边界。计算分析时，为保证气腔管道内部流动是紊流状态，因此设置计算模型为标准$\kappa - \varepsilon $湍流模型^[20]，湍流的强度设置为0.05。初始化后，为计算便捷有效，默认实际大气压强，设置仿真步长0.01 s，求解步数定为100步。本研究截取不同鸭嘴气体流经气腔管道的中心截面进行仿真，比较不同进气口与出气口的压力分布图。

根据流体力学流体运动学知识，出口气流分析(图6)可知，只有中间端进气圆弧式喷嘴产生细微的气压差，同时右边的两侧气流场与左边基本相同，外侧气流场对种子的作用就会更明显，内侧气压并不会影响气吹式清种装置对种子的清落效果。

图  6  不同结构喷嘴气流分布图

Figure  6.  Airflow distribution of nozzles with different structures

下载: 全尺寸图片 幻灯片

3.2   Ansys气流仿真数据分析

对4组不同的气吹式清种装置结构在Ansys软件中气流仿真进行数据分析，定义两端进气圆弧式喷嘴结构为Z₁，两端进气直孔式喷嘴结构为Z₂，中间端进气直孔式喷嘴结构为Z₃，中间端进气圆弧式喷嘴结构为Z₄，进而分析气吹式清种装置气流出口绝对气压、入口处气体流速及出口处气体流速的变化情况。借助Ansys软件从Fluent模块中导出上述指标在各步数下的数值，运用软件Origin进行数据分析，如图7所示。

图  7  仿真100步数下的散点图

Figure  7.  Simulation of scatter plot under 100 steps

下载: 全尺寸图片 幻灯片

由图7a可知，Z₁绝对气压范围为101 517~101 556 Pa，Z₂数值范围为101 525~101 601 Pa，Z₃数值范围为101 501~101 632 Pa，Z₄数值范围为101 515~101 534 Pa。因Z₄比Z₃的数值波动范围小，Z₁比Z₂的数值波动范围小，说明圆弧式喷嘴出口处绝对气压较直孔式喷嘴稳定；由图7b可知，4组不同结构的气体流速幅度无明显差异，但入口处气体流速有较大差别，除开气体流速为0 m/s的情况下，Z₁数值范围为4.323 90~4.681 77 m/s，Z₂数值范围为4.638 45~5.150 55 m/s，Z₃数值范围为9.625 35~9.902 55 m/s，Z₄数值范围为8.779 15~9.065 80 m/s。因Z₃比Z₂的入口处气体流速要快约一倍，Z₄比Z₁的要快约一倍，说明中间端进气相对于两端进气有着减弱损耗的优势；由图7c可知，气流经气压腔在到达出口处前，经过阶梯气腔，会造成气流堵塞，因对流作用此时入口处气体流速为0 m/s，出口处只存在大气压，故气体流速在进气压作用下会加快，此时流速将是最大值，进而慢慢下降，恢复稳定，但在气压腔内经对流作用气压释放完的一瞬间，出口处气体流速将为0 m/s，因无气体可释放，这会大大影响清落种子的可靠性，故出口处气体流速为0 m/s的次数出现越少越好。若出口处气体流速与入口处相比越小，则说明该结构的气流损耗越小；气体流速变化幅度越小，说明流阻系数越小。除开气体流速为0 m/s的情况，Z₁数值范围为12.758 5~16.796 7 m/s，Z₂数值范围为10.779 3~17.759 8 m/s，Z₃数值范围为11.687 7~17.216 1 m/s，Z₄数值范围为8.336 2~15.848 3 m/s，且0 m/s数值在Z₁、Z₂、Z₃和Z₄中分别有22、55、50和23个，因此Z₄数值较为稳定。

经上述仿真分析可得：Z₄在阶梯气腔内气流相对稳定，气压损耗最少，对于气吹式清种装置的稳压作用是最可靠的。

4.   台架试验

4.1   试验台架搭建与测试

所设计的三级清种装置的清种效果需通过试验检验，本研究按原样机尺寸搭建播种试验台进行测试。对于甘蓝、花椰菜、娃娃菜等种子，为保证有利于种苗生长且适应现有市面上移栽机械工作需求，蔬菜育苗一般使用50~288孔穴盘^[21]，根据原样机参数选取162穴。试验台的机架部分通过机械加工完成，滤轴、气吹式清种装置、滑轴及滑轨、排种器、配气盘稳压盘、除尘风机接口等复杂零件通过3D打印机打印，再购买其余标准零部件，组装成试验台架，进行单粒充种台架试验。

所搭建的试验台的可靠性需要进行测试，本研究先开展无清种装置下的测试试验，对应原样机的工作参数，设置排种器转速为0.17 r/s，排种器负压为0.6 kPa，种子直径为2.5 mm的单粒充种试验，重复试验10次。单圈充种穴数均值及单圈单粒充种穴数均值为每组试验下以轴向的1圈为研究对象，转1圈计算1次，转10圈求得平均值。经过高速视频目标跟踪测量软件V1.0进行测算，计算得无清种装置下的平均单粒充种合格率为(81.78±2.2)%，与实地试播测算的原样机的充种合格率(80%)相比基本一致。因此，该试验台架可用于后续的四因素三水平二次回归正交台架试验。

4.2   试验实施与结果分析

在实际工作中，排种器负压对种子吸附状态有重大影响，排种器转速变化会影响种子在排种器运动过程中的受力情况，种子大小所带来的质量变化会对种子在气孔处的受力平稳性产生影响，气吹式清种装置正压会使得种子在不同吸附状态下与排种器吸力产生气流对冲现象，故该试验的参数选取为排种器转速(A)、种子直径(B)、排种器负压(C)及清种装置正压(D)。结合原样机工作参数及种子测量数据、前期吸附试验确定，排种器转速为0.11~0.17 r/s，种子直径为1.5~2.5 mm，排种器负压为0.2~1.0 kPa，清种装置正压为0.1~0.5 kPa。

为探寻本研究的清种装置的清种效果，以清种装置正压、排种器负压、排种器转速、种子大小作为试验因素，以单粒充种合格率、重播率作为试验指标，开展四因素三水平二次回归正交试验，试验因素水平如表1所示。采用Box-Behnken design(RSM-BBD)设计了试验方案，试验共进行29组，每组试验重复3次，试验方案如表2所示。

表  1  清种试验因素水平

Table  1.  Factor and level of seed cleaning test

水平 Level	排种器转速/(r·s−1) Speed of seeding device (A)	种子直径/mm Seed diameter (B)	排种器负压/kPa Negative pressure of seeding device (C)	清种装置正压/kPa Positive pressure of seed cleaning device (D)
−1	0.11	1.5	0.2	0.1
0	0.14	2.0	0.6	0.3
1	0.17	2.5	1.0	0.5

下载: 导出CSV 

| 显示表格

表  2  清种试验方案与结果

Table  2.  Seed cleaning test design and result

试验方案 Testing program	排种器转速/(r·s−1) Speed of seeding device (A)	种子直径/mm Seed diameter (B)	排种器负压/kPa Negative pressure of seeding device (C)	清种装置正压/kPa Positive pressure of seed cleaning device (D)	单粒充种合格率/% Qualified rate of single seed filling (Y1)	重播率/% Reseeding rate (Y2)
1	0.11	1.5	0.6	0.3	87.04	0.10
2	0.17	1.5	0.6	0.3	95.31	0.20
3	0.11	2.5	0.6	0.3	86.89	0.13
4	0.17	2.5	0.6	0.3	89.70	0.30
5	0.14	2.0	0.2	0.1	87.84	0.20
6	0.14	2.0	1.0	0.1	95.14	0.77
7	0.14	2.0	0.2	0.5	87.70	0.10
8	0.14	2.0	1.0	0.5	87.37	0.20
9	0.11	2.0	0.6	0.1	86.62	0.10
10	0.17	2.0	0.6	0.1	90.76	0.60
11	0.11	2.0	0.6	0.5	86.41	0.10
12	0.17	2.0	0.6	0.5	85.36	0.20
13	0.14	1.5	0.2	0.3	87.74	0.07
14	0.14	2.5	0.2	0.3	94.64	0.70
15	0.14	1.5	1.0	0.3	93.73	0.60
16	0.14	2.5	1.0	0.3	90.27	0.37
17	0.11	2.0	0.2	0.3	86.75	0.07
18	0.17	2.0	0.2	0.3	90.59	0.40
19	0.11	2.0	1.0	0.3	93.38	0.57
20	0.17	2.0	1.0	0.3	88.19	0.20
21	0.14	1.5	0.6	0.1	88.39	0.17
22	0.14	2.5	0.6	0.1	92.21	0.57
23	0.14	1.5	0.6	0.5	88.57	0.30
24	0.14	2.5	0.6	0.5	86.72	0.10
25	0.14	2.0	0.6	0.3	93.76	0.60
26	0.14	2.0	0.6	0.3	92.73	0.60
27	0.14	2.0	0.6	0.3	93.56	0.53
28	0.14	2.0	0.6	0.3	92.97	0.60
29	0.14	2.0	0.6	0.3	94.70	0.83

下载: 导出CSV 

| 显示表格

根据二次正交台架试验，得到单圈充种穴数及单圈单粒充种穴数，对29组试验中的每组单圈单粒充种合格率，重播率取平均值，如表2所示。

运用Design-Expert 10.0.3软件进行数据分析，得到单粒充种合格率(Y₁)和重播率(Y₂)的方差表(表3)。Y₁的回归模型达到极显著(P<0.01)，说明该回归模型与自变量具有显著的函数关系；因素A、C均对Y₁影响显著(P<0.05)，D对Y₁影响极显著，BC交互作用极显著，AC、CD交互作用显著，A²和D²影响极显著，其他因素影响不显著。将表3中不显著因素剔除后，得到回归方程如下：

表  3  清种试验的方差分析¹⁾

Table  3.  Analysis of variance of seed filling test

方差来源 Source	单粒充种合格率 Qualified rate of single seed filling					重播率 Reseeding rate
	平方和 SS	均方 MS	F	P		平方和 SS	均方 MS	F	P
模型 Model	247.92	17.71	6.30	0.0007**		1.45	0.10	7.46	0.0003**
A	13.70	13.70	4.88	0.0444*		0.06	0.06	4.14	0.0612
B	0.01	0.01	3.63×10−3	0.9528		0.04	0.04	3.20	0.0951
C	13.70	13.70	4.88	0.0444*		0.11	0.11	8.23	0.0124*
D	29.55	29.55	10.52	0.0059**		0.17	0.17	11.95	0.0038**
A×B	7.45	7.45	2.65	0.1256		1.23×10−3	1.23×10−3	0.09	0.7706
A×C	20.39	20.39	7.26	0.0175*		0.12	0.12	8.84	0.0101*
A×D	6.73	6.73	2.40	0.1439		0.04	0.04	2.89	0.1115
B×C	26.83	26.83	9.55	0.0080**		0.18	0.18	13.34	0.0026**
B×D	8.04	8.04	2.86	0.1129		0.09	0.09	6.49	0.0232*
C×D	14.55	14.55	5.18	0.0391*		0.06	0.06	3.98	0.0658
A2	53.38	53.38	19.00	0.0007**		0.38	0.38	27.28	0.0001**
B2	7.74	7.74	2.76	0.1191		0.17	0.17	12.19	0.0036**
C2	4.43	4.43	1.58	0.2299		0.04	0.04	3.10	0.1000
D2	73.12	73.12	26.03	0.0002**		0.23	0.23	16.26	0.0012**
残差 Residual	39.33	2.81				0.19	0.01
失拟项 Lack of fit	36.95	3.70	6.22	0.0465*		0.14	0.01	1.07	0.5162
误差 Error	2.38	0.59				0.05	0.01
总和 Total	287.24					1.64
1) A：排种器转速；B：种子直径；C：排种器负压；D：清种装置正压；“*”和“**”分别表示0.05和0.01的显著水平 　1) A: Rotation speed of seed feeder; B: Seed diameter; C: Negative pressure of seed feeder; D: Positive pressure of seed cleaning device; “*” and “**” indicate 0.05 and 0.01 significance levels respectively

下载: 导出CSV 

| 显示表格

Y₂的回归模型达到极显著，说明该回归模型与自变量具有显著的函数关系；因素C对Y₂影响显著，D对Y₂影响极显著，BC交互作用极显著，AC、BD交互作用显著， A²、B²和D²影响极显著，其他因素影响不显著。通过F检验，得出影响重播率的主次顺序为清种装置正压、排种器负压、排种器转速、种子直径。将表中不显著因素剔除后，得到回归方程如下：

应用Design-Expert10.0.3软件分析排种器转速、种子直径、排种器负压、清种装置正压4个因素分别对单粒充种合格率和重播率的交互影响，得到单粒充种合格率响应曲面如图8所示。

图  8  各因素交互作用于单粒充种合格率的响应曲面

Figure  8.  Response surface of the interactive effects of different factors on qualified rate of single seed filling

下载: 全尺寸图片 幻灯片

由图8a可知，单粒充种合格率随着播种器转速和排种器负压增加而先增加后减小，排种器转速在0.14~0.17 r/s、排种器负压在0.6~1.0 kPa时单粒充种合格率较高；由图8b可知，单粒充种合格率随排种器负压和种子直径增加而先增加后减少，种子直径在1.5~2.0 mm、排种器负压在0.6~1.0 kPa时单粒充种合格率较高；由图8c可知，单粒充种合格率随着清种装置正压和排种器负压的增加而先增加后减少，清种装置正压在0.1~0.3 kPa、排种器负压在0.6~1.0 kPa时单粒充种合格率较高。

重播率的响应曲面如图9所示。由图9a可知，重播率随着播种器转速和排种器负压增加而先增加后减小，排种器转速在0.11~0.14 r/s、排种器负压在0.2~0.6 kPa时重播率较低；由图9b可知，重播率随排种器负压和种子直径增加而先增加后减少，种子直径在1.5~2.0 mm、排种器负压在0.2~0.6 kPa时重播率较低；由图9c可知，重播率随着清种装置正压和种子直径的增加而先增加后减少，清种装置正压在0.1~0.3 kPa、排种器负压在1.5~2.0 kPa时重播率较低。

图  9  各因素交互作用于重播率的响应曲面

Figure  9.  Response surface of the interactive effects of different factors on reseeding rate

下载: 全尺寸图片 幻灯片

为获得播种合格率最高和重播率最小的工作参数，以排种器转速等为自变量，以播种合格率最高和重播率最小为目标，通过Design-Expert 10软件对二次正交组合试验结果进行分析，得到播种机的最佳工作参数组合：当排种器转速在0.17 r/s、种子直径在1.5 mm、排种器负压在0.6 kPa、清种装置正压在0.3 kPa时，平均单粒充种合格率为95.31%，且平均重播率为0.2% (表2)。

4.3   对照试验

为探寻三级清种装置清种效果，开展单级清种装置试验，在最优工作参数：排种器转速在0.17 r/s、种子直径在1.5 mm、排种器负压在0.6 kPa、清种装置正压在0.3 kPa下，分别使用各单级清种装置进行3组重复试验，求得各装置单粒充种合格率、重播率 (表4)，形成对照数据，进一步验证试验效果。

表  4  单级清种装置试验数据表

Table  4.  Data sheet for single stage cleaning device test

清种装置类型 Type of seed cleaning device	单圈充种穴数 No. of seed filling holes per circle	单圈单粒充种穴数 No. of single seed filling holes per circle	单粒充种合格率/% Qualified rate of single seed filling	重播率/% Reseeding rate
机械式 Mechanical	51.1	48.3	89.50	5.19
气吹式 Air blown	51.3	50.2	92.96	2.10
电控式 Electrically controlled	50.8	49.3	91.23	2.53

下载: 导出CSV 

| 显示表格

对比表4与表2的最优参数结果可知，集机械式清种、气吹式清种、电控式清种为一体的三级清种装置相比于机械式清种装置，单粒充种合格率提高了5.81个百分点，重播率降低了4.99个百分点；相比于气吹式清种装置，单粒充种合格率提高了2.35个百分点，重播率降低了1.90个百分点；相比于电控式清种装置，单粒充种合格率提高了4.08个百分点，重播率降低了2.33个百分点。因此，本研究设计的三级清种装置具有提高单粒充种合格率、降低重播率的作用。

4.4   验证试验

为进一步验证最佳参数下清种装置的清种性能，选取排种器转速在0.17 r/s、种子直径在1.5 mm、排种器负压在0.6 kPa、清种装置正压在0.3 kPa时进行试验验证，进行5组试验，取平均值，得到对比试验结果如表5所示。

表  5  三级清种装置清种性能验证试验结果

Table  5.  Validation test result of seed cleaning performance of the three level seed cleaning device

项目 Itern	单粒充种合格率/% Qualified rate of single seed filling	重播率/% Reseeding rate
台架试验 Bench test	95.31	0.2
最优参数验证试验 Optimal parameter verification test	94.86	0.5
误差 Error	0.45	0.3

下载: 导出CSV 

| 显示表格

通过进行验证试验，结果表明最佳参数组合下蔬菜穴盘育苗播种机安装三级清种装置后平均单粒合格率为94.86%，平均重播率为0.5%，验证了台架试验的可靠性，且最佳工作参数组合下的试验指标能满足蔬菜播种的农艺要求。

5.   结论

1)通过理论分析得到三级清种装置的滤孔直径范围为1.50~5.00 mm，滤轴直径为10 mm；滑轴最低速度为0.29 mm/s，根据种子尺寸大小，偏心凸轮偏心距为3.2 mm，同时根据安装位置及机械疲劳强度，设置滑轴直径为6 mm，弹簧线径为8 mm，滑轨孔径为7 mm。运用流体动力学分析得出种子与喷嘴气流出口之间的距离为10 mm，进气口直径为8 mm，进气口压强为0.3 kPa。

2)针对气吹式清种装置，借助Ansys软件中的ICEM CFD模块及Fluent模块开展以气吹式清种装置进气口、出气口结构为试验因素，气腔内气流稳定性为试验指标的二因素二水平仿真试验，通过仿真效果对比得出，中间端进气圆弧式喷嘴气吹式清种装置最能满足理想清种效果。

3)通过开展正交台架试验，结果表明，当排种器转速在0.17 r/s、种子直径在1.5 mm、排种器负压在0.6 kPa、清种装置正压在0.3 kPa时，采用三级清种装置的平均单粒充种合格率为95.31%，平均重播率为0.2%，此时清种效果比较好。再通过单级清种对照试验、验证试验表明，安装集机械式清种、气吹式清种、电控式清种为一体的三级清种装置具有提高单粒充种合格率、降低重播率的作用。