底肥定点深施机地轮组件优化
宋朱军1, 吴雪梅1,*, 李国昌1, 王芳1, 刘红芸1, 张富贵1, 张培坤2
1.贵州大学 机械工程学院,贵州 贵阳 550025
2.福建永顺机械有限公司,福建 南平 353500
*通信作者,吴雪梅,E-mail: xm_wu@163.com

作者简介:宋朱军(1992—),男,河北邯郸人,硕士研究生,主要从事机械电子工程方面的研究。E-mail: 815645779@qq.com

摘要

针对底肥定点深施肥机地轮组件驱动力不足,造成施肥点与标记点误差较大的问题,采用一次回归正交设计进行参数优化试验。选择影响施肥点与标记点误差的主要因素(机具前进速度、叶片深度、轮毂半径、叶片面积)作为试验因素,以地轮组件所受的驱动力矩为评价指标,进行4因素4水平的正交试验。仿真试验表明,各因素对施肥点与标记点误差的影响从大到小依次为叶片面积>轮毂半径>叶片深度>机具前进速度。在试验范围内,最优参数组合为:机具前进速度2.3 m·s-1,叶片深度8 cm,轮毂半径17.50 cm,叶片面积52 cm2。在此条件下,地轮组件的转动力矩为101.96 N·m,较优化前增大11.76%。田间验证试验表明,优化前后,标记点与施肥点的位置误差由4.23 cm缩小至1.99 cm,实现了施肥点的准确标记功能。研究结果可为底肥定点深施机的结构及作业参数优化提供参考依据。

关键词: 底肥定点深施肥机; 地轮; 优化; 正交试验
中图分类号:S224.2 文献标志码:A 文章编号:1004-1524(2019)05-0829-08
Optimization of ground wheel components of base fertilizer fixed-point deep applicator
SONG Zhujun1, WU Xuemei1,*, LI Guochang1, WANG Fang1, LIU Hongyun1, ZHANG Fugui1, ZHANG Peikun2
1. College of Mechanical Engineering, Guizhou University, Guiyang 550025, China
2. Fujian Yongshun Machinery Co., Ltd., Nanping 353500, China
Abstract

In order to solve the large deviation of fertilization point and marking point caused by the lack of driving force of the ground wheel assembly of base fertilizer fixed-point deep fertilization machine, an optimal parameter experiment is carried out by using a regression orthogonal design. The main factors affecting the deviation between the fertilization point and the marking point were selected as the advance speed of the apparatus, the blade depth, the hub radius and the blade area. The four-factor four-level orthogonal experiment was carried out with the driving torque of the ground wheel assembly as the evaluation index. Simulation experiments showed that the influence of factors on deviation between the fertilization point and the marking point decreased as the blade area>hub radius>blade depth>advance speed of the apparatus. The optimal combination of parameters were as follows: the advance speed of the apparatus was 2.3 m·s-1, the blade depth was 8 cm, the hub radius was 17.50 cm, and the blade area was 52 cm2. Under this condition, the driving torque provided by the soil for the ground wheel assembly was 101.96 N·m, which was increased by 11.76%. Field experiment results indicated that the position deviation between the marking point and the fertilization point was reduced from 4.23 cm to 1.99 cm after optimization, which realized the accurate marking function. The research could provide reference for the optimization of the structure and operation parameters of the base fertilizer fixed-point deep applicator.

Keyword: base fertilizer fixed-point deep applicator; ground wheel; optimization; orthogonal experiment

合理的底肥深施技术对促进作物生长有积极作用[1]。对于有株距要求的作物, 如烤烟、玉米、辣椒等, 移栽前进行底肥定点深施, 可促进作物根系生成, 进而提高作物品质和产量[2, 3, 4]。目前, 底肥深施机多为条施, 不能实现定点深施和标记施肥点[5, 6, 7]。针对此现象, 研究人员曾开发出一款自走式底肥定点深施机[8, 9]。该施肥机的地轮组件采用空心轮毂与叶片螺纹连接, 工作时通过链条传动将动力传至外槽轮式排肥器驱动其排肥, 但由于地轮组件与垄顶接触时产生的驱动力不足, 会出现打滑现象, 造成施肥点与标记点间误差较大。

目前, 针对地轮性能的研究比较广泛[10, 11, 12], 但多停留在理论分析上, 少有研究者利用动力学仿真来分析地轮各部件的受力情况。为此, 本研究针对前述自走式底肥定点深施机地轮组件打滑, 造成施肥点与标记点间误差较大的问题, 借助工程离散元仿真软件EDEM 6.0建立仿真土壤模型, 对该施肥机的地轮组件进行颗粒离散元法的细观力学分析。类似方法已广泛应用于农业机械的各个方面。研究证明, 通过观察土壤颗粒的运动情况和位移情况, 可以了解地轮组件的受力情况[13]。此外, 本研究还采用4因素4水平的正交试验方法, 应用Design-Expert软件对试验数据进行曲面响应分析, 以优化地轮组件的结构参数设计, 实现精准标记施肥点的目的。

1 地轮组件结构与工作原理
1.1 地轮组件结构

地轮组件主要由4部分组成(图1), 主要部分为空心轮毂和15个叶片。叶片通过螺纹连接固定在轮毂的对称面上, 在设计上呈上窄下宽。这种形状的设计可以保证叶片有更大的面积接触土壤, 从而增大其工作中的驱动力。叶片上端为有螺纹的连接杆, 使叶片与轮毂连接。轮毂设计为空心, 存放标记所用的石灰粉。从石灰进料孔填装石灰粉, 石灰粉通过石灰标记孔对施肥位置进行标记。

图1 地轮组件结构三维示意图
1, 叶片; 2, 石灰标记孔; 3, 空心轮毂; 4, 石灰进料孔。
Fig.1 Schematic of three dimensional structure of ground wheel components
1, Blade; 2, Lime marking hole; 3, Hollow wheel hub; 4, Lime feed hole.

1.2 工作原理

底肥定点深施机的地轮组件是施肥装置的动力来源装置, 其依靠叶片与垄顶土壤的相互作用使地轮转动, 并通过链轮链条驱动排肥器排肥, 实现被动施肥的目的。

工作时, 地轮组件靠自身重力压附在垄的顶部, 在整机动力源的牵引力作用下, 叶片在土壤中接触产生相互作用, 从而驱动地轮组件转动。地轮组件通过链条传动将产生的动力传送至外槽轮式排肥器, 驱动其排肥。通过旋转叶片可以调节叶片与轮毂的距离, 改变轮毂中心与垄顶土壤的距离, 实现株距调节。地轮转动过程中, 空心轮毂中的石灰粉通过标记孔落在深施肥土壤上方, 以实现对施肥点进行标记的功能。

2 地轮组件设计及结构参数分析
2.1 空心轮毂

为了保证地轮组件在工作中能够满足精准标记施肥点的目的, 地轮采用空心设计, 并在轮毂表面开设3个对称石灰进料孔, 端盖一端固定于轮毂表面。在轮毂的中心设有用来固定链轮的伸出轴。地轮未转动前, 在垄顶上的受力分析[14]如图2所示。

图2 地轮静止在垄顶时的受力分析Fig.2 Force analysis of ground wheel components at ridge top in a static state

整机开始工作时, 轮毂并未转动, 处于临界状态, 其受力分析下:

G-Fy-Ry=0; (1)

Fx-Rx=0。 (2)

式(1)、(2)中:G, 地轮重力; FxFy, 分别表示牵引力F在水平和竖直方向的分力; RxRy, 分别表示土壤对地轮的作用力R在水平和竖直方向的分力。

μ 为土壤与轮毂摩擦系数。当地轮开始转动时,

Fx> Rx+μ (G-Fy-Ry)。(3)

此时, 轮毂与垄顶土壤由静摩擦转为滚动摩擦。轮毂工作力矩为

M1=(Rx-Fx)Rd, (4)

式(4)中Rd为轮毂半径。

2.2 叶片

叶片作为地轮组件的重要受力部件, 15个叶片分别间隔24° 呈对称分布。叶片插入垄顶土壤后, 受到土壤对其表面的挤压, 在工作时其与土壤的作用为驱动地轮转动提供主要驱动力。其作用力取决于垄顶土壤的紧实度及土壤与叶片的接触面积。在地轮工作过程中, 叶片与土壤的相互作用[14]如图3所示。

图3 工作中地轮与叶片的力矩分析Fig.3 Torque analysis of ground wheel components and blade under working conditions

叶片工作力矩为

M2=q(Rd+h), (5)

式(5)中:q, 土壤作用力; h, 叶片深度, 即地轮组件中通过螺纹与轮毂连接后叶片中心到轮毂面的距离。

由式(5)可知, 影响叶片工作力矩的主要因素为qh。通过改变旋入轮毂深度, 可调节叶片深度。

2.3 参数分析

在式(1)、(2)中, 牵引力F由发动机功率P确定, 即

F= Pv, (6)

式(6)中v为整机前进速度, 则式(4)可以改写为

M1=(Rx-sinθ · Pv)Rd, (7)

式(5)中θ 如图2所示, 为牵引力方向与水平方向的夹角。

土壤作用力q由土壤紧实度Q和叶片面积S决定, 即

q=Q· S, (8)

故式(5)可写为

M2=QS(Rd+h), (9)

则地轮组件工作总力矩可写为

M=M1+M2=(Rx-sinθ · Pv)Rd+QS(Rd+h)。 (10)

针对底肥定点深施机存在的地轮驱动力不足的问题, 结合地轮组件结构参数(图4), 分析式(10)可和, 地轮驱动力矩由整机前进速度v、轮毂半径Rd、叶片面积S和叶片深度h确定。

图4 地轮组件结构参数
1, 机具前进速度; 2, 叶片深度; 3, 轮毂半径; 4, 叶片面积。
Fig.4 Schematic of structural parameters of ground wheel assembly
1, Advance speed of apparatus; 2, Blade depth; 3, Hub radius; 4, Blade area.

机具前进速度是影响整机工作效率和稳定性的重要参数。根据式(10)可知, 要增加转动力矩, 可通过增加机具前进速度实现。研究发现, 机具前进速度过快容易导致地轮打滑[9]。由式(10)可知:轮毂半径越大, 转动力矩越大, 但轮毂的尺寸会同时影响地轮组件的整体规格与质量; 叶片面积决定着机具与土壤接触和作用的面积, 直接影响叶片转动力矩; 叶片深度增加, 同样会增大地轮组件转动力矩。

3 地轮组件参数优化仿真试验

采用EDEM 6.0软件建立地轮组件的仿真模型, 分析各结构参数变化对地轮驱动力矩的影响, 优化地轮组件的结构参数。

3.1 仿真参数确定

3.1.1 地轮组件和土壤颗粒模型

为了真实反映田间土壤, 采用离散元法建立起垄后的土壤模型。起垄后的土壤颗粒表面黏附力较小, 具有散粒体物料特性和压缩性, 因此确定以Hertz-Mindlin with bonding作为土壤与地轮组件的接触模型。离散元仿真参数包括颗粒参数的本征参数和土壤颗粒接触参数, 确定模型参数设置详见表1

表1 离散元仿真微观参数 Table 1 Microscopic parameter of discrete element simulation

3.1.2 接触力学模型及其参数选取

通过文献参考和颗粒参数标定等方法确定离散元仿真参数[15, 16, 17, 18, 19, 20, 21]。通过确定仿真接触模型及参数, 同时为了减少仿真运行时间和存储空间, 建立了地轮组件作业的离散元虚拟土壤模型(1 800 mm × 500 mm × 600 mm), 如图5所示。应用Pro/E软件分别创建16种参数匹配下的地轮组件模型, 并保存为.igs格式导入EDEM 6.0软件中进行仿真, 图5为其中一种参数匹配下, 地轮组件与土壤相互作用时的仿真过程。

图5 土壤-地轮组件相互作用模型Fig.5 Interaction model of soil-ground wheel components

3.2 地轮组件参数优化正交试验及仿真分析

3.2.1 试验方案设计

根据前期分析, 采用4因素4水平的正交试验, 以机具前进速度、叶片深度、轮毂半径、叶片面积为影响因素, 以地轮组件所受转动力矩为评价指标, 选用L16(45)正交表, 开展地轮组件所受工作力矩的仿真。试验因素水平设计如表2所示。

表2 试验因素水平 Table 2 Experiment factors and levels

选取16种参数下的地轮组件, 开展仿真试验。试验分为2个步骤:先通过颗粒工厂生成仿真模型, 然后将地轮组件以设计的参数开始仿真运动。地轮组件仿真时只进行直线运行, 以模拟实际工作中当地轮组件发生打滑或卡住时, 土壤对地轮组件产生的最大转动力矩。

3.2.2 试验结果与分析

试验方案及结果详见表3

表3 试验方案与结果 Table 3 Experimental design and results

地轮组件开始进入及离开土壤模型时, 工作阻力仿真数值变化较大, 不符合运动规律; 所以, 在数据处理时, 舍弃开始与结束时的几个点, 提取中间有效的数据点并取平均值。

利用Design-Expert 8.0.6软件进行数据处理和统计, 结果如表4所示, 一次回归模型为

y=61.41+7.61x1+9.51x2+10.23x3+11.71x4, (11)

式(11)中:y, 力矩; x1, 行进速度编码值; x2, 叶片深度编码值; x3, 轮毂半径编码值; x4, 叶片面积编码值。

表4 各因素对转动力矩影响的方差分析 Table 4 Variance analysis of influence of various factors on driving torque

模型的显著性检验结果为F=13.42, P< 0.01, 决定系数R2=0.829 9, 表明回归模型达到极显著水平(P< 0.01), 拟合程度较好。

分析表4可知:x2x3x4对转动力矩的影响达到极显著水平(P< 0.01), x1对转动力矩的影响达到显著水平(P< 0.05), 影响程度从大到小依次为x4> x3> x2> x1。由此可知, 在对地轮组件进行改进时, 适当增加轮毂半径和叶片面积对提高其转动力矩具有明显效果。

3.3 地轮组件最优工作参数试验

根据地轮组件仿真试验结果对各因素进行优化, 以地轮转动力矩为指标, 以机具前进速度、叶片深度、轮毂半径、叶片面积为影响因素, 应用Design-Expert 8.0.6的曲面优化功能做进一步优化, 最优值为:整机前进速度2.3 m· s-1, 叶片深度8 cm, 轮毂半径17.50 cm, 叶片面积52 cm2。此时, 土壤为地轮组件提供的最大力矩为101.96 N· m。优化过后的结果如图6所示, 优化后的最大转动力矩较优化前增大了11.76%。

图6 优化后的试验结果Fig.6 Optimized results

4 田间验证试验

为了验证优化后参数组合的合理性, 进行田间施肥点与标记点试验。于2017年7月在贵州大学农场选取一块长30 m、宽10 m的烤烟田, 使用优化改进后的机具(如图7-a所示)开展试验, 去除前后5 m准备区, 对20 m的有效区域进行测量。测量之前, 在每个石灰标记的中心插入标记杆, 标记各个标记点的位置(如图7-b所示), 再除去标记点附近的土壤, 以显现出肥料集中点。

图7 田间试验
a, 改进后机具的田间试验; b, 施肥点与标记点的测量。
Fig.7 Field experiment with improved ground wheel components
a, Field experiment with improved fertilizer applicator; b, Measurement of fertilization points and marking points.

选择施肥点与标记点的位置误差、施肥长度、施肥点间距3个参数进行测量, 对比优化前后样机的试验结果, 如表5所示。

表5 田间试验结果对比 Table 5 Comparison of field experiment results with fertilizer applicator before and after optimization

表5可以看出, 改进后的底肥定点深施机施肥点与标记点的位置平均误差、标准差均小于改进前, 说明本试验提出的改进方法对于减小底肥定点深施机施肥点与标记点的误差是有效的。改进前后施肥长度及施肥点间距基本保持不变。

5 结论

利用离散元分析软件EDEM 6.0建立仿真土壤模型, 进行地轮组件仿真。采用4因素4水平的正交试验方法, 获得优化目标与4个影响因素之间的回归方程及因素之间交互作用的响应曲面。分析各因素对转动力矩影响的方差发现:叶片深度、轮毂半径和叶片面积对转动力矩的影响达极显著水平。在对地轮组件改进时, 适当增加轮毂半径和叶片面积对提高其转动力矩具有明显效果。

运用Design-Expert 8.0.6软件对回归模型进行优化及验证。当整机行进速度为2.3 m· s-1、叶片深度为8 cm、轮毂半径为17.50 cm、叶片面积为52 cm2时, 土壤为地轮组件提供的最大转动力矩为101.96 N· m, 较优化前增大了11.76%。

采用优化后的参数组合进行田间验证试验。结果显示, 优化后的机具有效地减小了施肥点与标记点的误差, 标记点与施肥点的位置误差由优化前的4.23 cm减小为1.99 cm, 实现了施肥点精准标记的功能。

The authors have declared that no competing interests exist.

参考文献
[1] 伟利国, 张小超, 苑严伟, . 2F-6-BP1型变量配肥施肥机的研制与试验[J]. 农业工程学报, 2012, 28(7): 14-18.
WEI L G, ZHANG X C, YUAN Y W, et al. Design and experiment of 2F-6-BP1 variable rate assorted fertilizer applicator[J]. Transactions of the Chinese Society of Agricultural Engineering, 2012, 28(7): 14-18. (in Chinese with English abstract) [本文引用:1]
[2] 罗原骏, 蒲玉琳, 龙高飞, . 施肥方式对土壤活性有机碳及碳库管理指数的影响[J]. 浙江农业学报, 2018, 30(8): 1389-1397.
LUO Y J, PU Y L, LONG G F, et al. Effects of fertilization on soil active organic carbon and carbon pool management index[J]. Acta Agriculturae Zhejiangensis, 2018, 30(8): 1389-1397. (in Chinese with English abstract) [本文引用:1]
[3] 袁文胜, 金梅, 吴崇友, . 国内种肥施肥机械化发展现状及思考[J]. 农机化研究, 2011, 33(12): 1-5.
YUAN W S, JIN M, WU C Y, et al. Development status and prospect of seed manure fertilizing mechanization in China[J]. Journal of Agricultural Mechanization Research, 2011, 33(12): 1-5. (in Chinese with English abstract) [本文引用:1]
[4] LORENCOWICZ E, UZIAK J. Repair cost of tractors and agricultural machines in family farms[J]. Agriculture and Agricultural Science Procedia, 2015, 7: 152-157. [本文引用:1]
[5] 金鑫, 李倩文, 苑严伟, . 2BFJ-24型小麦精量播种变量施肥机设计与试验[J]. 农业机械学报, 2018, 49(5): 84-92.
JIN X, LI Q W, YUAN Y W, et al. Design and test of 2BFJ-24 type variable fertilizer and wheat precision seed sowing machine[J]. Transactions of the Chinese Society for Agricultural Machinery, 2018, 49(5): 84-92. (in Chinese with English abstract) [本文引用:1]
[6] 李绪. 一垄双行深施肥机设计与试验[D]. 贵阳: 贵州大学, 2018.
LI X. Design and experiment of a fertilizer deep applicator for twin-row within one ridge[D]. Guiyang: Guizhou University, 2018. [本文引用:1]
[7] 韩云霞. 2BFJ-6型自动变量施肥机控制系统研究[D]. 长春: 吉林大学, 2004.
HAN Y X. Study on control system of 2BFJ-6 automatic variable rate fertilizer applicator[D]. Changchun: Jilin University, 2004. (in Chinese with English abstract) [本文引用:1]
[8] 夏志林, 李国昌, 吴雪梅, . 烟用定点施肥装置凸轮机构运动学分析与仿真[J]. 中国农机化学报, 2017, 38(4): 69-74.
XIA Z L, LI G C, WU X M, et al. Kinematics analysis and simulation of fertilizer placement unit's cam mechanism for tobacco[J]. Journal of Chinese Agricultural Mechanization, 2017, 38(4): 69-74. (in Chinese with English abstract) [本文引用:1]
[9] 李国昌, 吴雪梅, 殷坤勤, . 烟用自走式起垄定点施肥装置设计及试验[J]. 农机化研究, 2018, 40(4): 104-109, 115.
LI G C, WU X M, YIN K Q, et al. Design and experiment of self-propelled ridge and point fertilizing apparatus for tobacco[J]. Journal of Agricultural Mechanization Research, 2018, 40(4): 104-109, 115. (in Chinese with English abstract) [本文引用:2]
[10] 张晨光. 免耕播种机防滑地轮设计与试验研究[D]. 哈尔滨: 东北农业大学, 2017.
ZHANG C G. Design and experimental study on anti-slip ground of no-tillage seeder[D]. Harbin: Northeast Agricultural University, 2017. (in Chinese with English abstract) [本文引用:1]
[11] 夏连明, 耿端阳, 王相友, . 免耕播种机行星啮合式防滑地轮设计与试验[J]. 农业机械学报, 2012, 43(8): 56-60.
XIA L M, GENG D Y, WANG X Y, et al. Design and experiment of planetary geared anti-sliding ground wheel for no-tillage planter[J]. Transactions of the Chinese Society for Agricultural Machinery, 2012, 43(8): 56-60. (in Chinese with English abstract) [本文引用:1]
[12] 夏连明, 耿端阳, 王相友. 行星啮合式防滑地轮结构参数的优化及性能试验[J]. 农业工程学报, 2012, 28(10): 53-58.
XIA L M, GENG D Y, WANG X Y. Structure parameters optimization and performance experiment for planetary geared anti-sliding ground wheel[J]. Transactions of the Chinese Society of Agricultural Engineering, 2012, 28(10): 53-58. (in Chinese with English abstract) [本文引用:1]
[13] 杨军伟, 孙慧男, 张卓青. 离散元法及其在农业工程中的应用综述[J]. 现代食品, 2015(15): 28-33.
YANG J W, SUN H N, ZHANG Z Q. A review on fundamentals of distinct element method and its applications in agricultural engineering realm[J]. Modern Food, 2015(15): 28-33. (in Chinese with English abstract) [本文引用:1]
[14] 高玉璐. 免耕播种机地轮滑移现象的研究[D]. 北京: 中国农业大学, 2002.
GAO Y L. Study on ground wheel slide of no tillage seeding-machine[D]. Beijing: China Agricultural University, 2002. (in Chinese with English abstract) [本文引用:2]
[15] 夏连明. 玉米精量播种机关键部件研究[D]. 哈尔滨: 东北农业大学, 2011.
XIA L M. Research on key components for maize precision planter[D]. Harbin: Northeast Agricultural University, 2011. (in Chinese with English abstract) [本文引用:1]
[16] 方会敏, 姬长英, TAGAR Ahmed Ali, . 秸秆-土壤-旋耕刀系统中秸秆位移仿真分析[J]. 农业机械学报, 2016, 47(1): 60-67.
FANG H M, JI C Y, TAGAR A A, et al. Simulation analysis of straw movement in straw-soil-rotary blade system[J]. Transactions of the Chinese Society for Agricultural Machinery, 2016, 47(1): 60-67. (in Chinese with English abstract) [本文引用:1]
[17] 丁启朔, 任骏, BELAL Eisa Adam, . 湿粘水稻土深松过程离散元分析[J]. 农业机械学报, 2017, 48(3): 38-48.
DING Q S, REN J, BELAL E A, et al. DEM analysis of subsoiling process in wet clayey paddy soil[J]. Transactions of the Chinese Society for Agricultural Machinery, 2017, 48(3): 38-48. (in Chinese with English abstract) [本文引用:1]
[18] ZHANG X C, LI H W, DU R C, et al. Effects of key design parameters of tine furrow opener on soil seedbed properties[J]. International Journal of Agricultural and Biological Engineering, 2016, 9(3): 67-80. [本文引用:1]
[19] 赵淑红, 刘宏俊, 谭贺文, . 仿旗鱼头部曲线型开沟器设计与性能试验[J]. 农业工程学报, 2017, 33(5): 32-39.
ZHAO S H, LIU H J, TAN H W, et al. Design and performance experiment of opener based on bionic sailfish head curve[J]. Transactions of the Chinese Society of Agricultural Engineering, 2017, 33(5): 32-39. (in Chinese with English abstract) [本文引用:1]
[20] 黄玉祥, 杭程光, 苑梦婵, . 深松土壤扰动行为的离散元仿真与试验[J]. 农业机械学报, 2016, 47(7): 80-88.
HUANG Y X, HANG C G, YUAN M C, et al. Discrete element simulation and experiment on disturbance behavior of subsoiling[J]. Transactions of the Chinese Society for Agricultural Machinery, 2016, 47(7): 80-88. (in Chinese with English abstract) [本文引用:1]
[21] UCGUL M, FIELKE J M, SAUNDERS C. 3D DEM tillage simulation: validation of a hysteretic spring (plastic) contact model for a sweep tool operating in a cohesionless soil[J]. Soil and Tillage Research, 2014, 144: 220-227. [本文引用:1]