雷哓晖 吕晓兰 李雪

摘要：喷杆喷雾机是目前普遍使用的一种大田植保机械，为了增加农药有效利用率和减少雾滴飘失，喷杆多加装罩盖装置。针对喷杆喷雾机生产研发中出现的损坏问题，以山东绿田农业机械有限公司生产的罩盖导流式喷雾机拨禾板固定支架为研究对象，利用有限元理论找出固定支架应力集中部位，并在其原有结构基础上进行拓扑优化设计，以增强其结构强度。在受损情况最为严重的大豆田中，拓扑优化后固定支架最大应力由原来的652.4 MPa降至 123.8 MPa，降了81%，机械抗振性也有大幅度的增强。该结果为企业生产优质的拨禾板式喷杆喷雾机提供技术支持。

关键词：喷杆喷雾机;拨禾板;固定支架;拓扑优化;应力分析;数字模型;仿真;优化设计

中图分类号：S491   文献标志码： A  文章编号：1002-1302（2019）14-0243-05

田间施药是大田作物种植管理中最重要的环节，可有效地预防田间作物病虫草害的发生。目前国内大田施药机械主要有2种：一种是以无线遥控多旋翼飞行器为载体的航空植保机械，一种是高地隙喷杆喷雾机。航空植保机械多适用于山区丘陵地区，以作业速度快、地形适应性强为特点，缺点是农药雾滴易飘失、价格高昂且载药量少[1-2]。高地隙喷杆喷雾机以平原作业为主，高地隙底盘设计可有效减少机具作业过程中对作物的碾压伤害，为提高农药的有效利用率及减少农药雾滴的飘失，喷雾机喷杆多加装风幕装置[3]和罩盖装置[4]。

拨禾板导流式喷雾属于罩盖喷雾，是一种结构简单、较风幕式喷雾价格便宜、可减少雾滴飘失的喷雾方法。它对高秆和低矮作物都适用。当装载拨禾板的喷雾机作业于低矮作物时，拨禾板底端紧贴作物顶部，倾斜喷出的农药药液会穿透作物冠层;当机具作业于高秆作物时，拨禾板可以拨开作物冠层，减少上部冠层对农药药液的拦截，使药液直接喷施于冠层下部。文献[4]对拨禾板的安装倾角及宽度都有详细的定义，山东绿田农业机械有限公司参考文献[4]研发设计一款拨禾板式喷杆喷雾机，但拨禾板与喷杆之间的固定支架作业时强度效果欠佳。本研究以山东绿田农业机械有限公司设计的拨禾板固定支架为研究对象，运用计算机拓扑优化理论对其进行优化设计。

1 模型建立

1.1 工作原理

拨禾板导流式喷杆喷雾机工作原理如图1所示，喷雾机液压驱动喷杆升降使得拨禾板置于作物冠层当中，喷杆在图1中从左至右作业，喷杆带动拨禾板、拨禾板固定支架以及喷头整体移动，拨禾板在喷雾机的带动下将作物冠层拨开，使得喷头喷出的农药雾滴喷施到作物冠层内部及下部。

1.2 数字模型

大田喷杆喷雾机多使用平面扇形雾喷头进行作业，根据喷头流量及扇面角度的不同，喷头在喷杆上一般以0.5～1.0 m 为间隔等间距排列（本试验研究固定支架单体间隔为0.8 m）。拨禾板式喷杆如图2所示，拨禾板通过固定支架安装于喷杆下方，喷头安装于拨禾板上方并与固定支架固接。为了增加结构强度，固定支架设计为平行四边形结构，有时为了增加拨禾板作业灵活性，固定支架平行四边形结构角度设计成可调节模式。

运用机械三维设计软件PTC/Creo，根据实际尺寸对拨禾板固定支架进行建模。拨禾板式喷杆喷雾机固定支架（简称原模型）设计图如图3所示，单个支架质量为1.18 kg，宽度为30 mm。

1.3 作业工况分析

喷杆喷雾机作业过程中作物冠层主要对拨禾板施加水平向左的推力，当作物茎秆折断时，阻力消失。以大田常见的秸秆种植作物（玉米、高粱、小麦、水稻、大豆、油菜）为研究对象，分析其种植密度和抗折力，为后续拨禾板支架的受力提供数据支持。有研究表明，玉米适合机械化作业且可增高其产量的种植密度范围为7.5万～9.0万株/hm2[5-6]，主茎秆机械抗折力因施氮量不同而有所改变，最高可达 26.88 N[7]。高粱适宜的最大种植密度为14.66万株/hm2[8]，主茎秆抗折力因品种的不同而有所不同，最大约 30 N[9]。小麦高产栽培的种植密度为240万株/hm2[10]，主茎秆抗折力最大值为 6.2 N[11]。高密度栽植水稻的种植密度可达33万株/hm2[12]，主茎秆抗折力最大为18.9 N[13]。大豆最大种植密度为 27万株/hm2[14]，主茎秆抗折力最大为 87 N[15]。油菜最大种植密度为60万株/hm2[16]，主茎秆抗折力最大为33.66 N[17]。经单位统一后，6种秸秆作物单位种植面积内最大抗折力总和如表1所示。由表1可知，大豆单位面积抗折力明显比其他秆作物大。

2 仿真与优化

2.1 理论模型

2.2 原模型仿真分析

使用有限元商业分析软件Hyperworks对原模型进行应力分析，以模型侧面为基准，通过壳网格拉伸的方法绘制边长为4 mm的六面体实体网格：网格单元总数2 198个，结点数 5 000个，网格质量良好。因喷雾机拨禾板固定支架材质为Q235，故所加载的材料属性为杨氏模量2.08×105 MPa、泊松比0.277、密度7.86×103 kg/m3、屈服强度235 MPa。图1中喷雾机水平向右行驶，若拨禾板完全深入到作物冠层中，则拨禾板受到作物茎秆水平向左的阻力最大，进而传递到各个固定支架底端。因拨禾板与地面成40°夹角，故单个固定支架受到宽度为B、长度为图2中L（0.35×sin50°=0.27 m）的面积内作物阻力的作用。分析中以作物最大抗折力进行计算，固定支架所受作用力及支架产生的应力如表2所示。由表2可知，B=0.8 m 时，因应力超过屈服强度，固定支架在小麦、大豆、油菜田中不可作业;B=0.5 m时，固定支架在小麦、大豆、油菜田中不可作业。

固定支架在6种受力工况下的应力集中区皆出现在平行四边形结构右上角喷杆连接处，图4为支架在玉米田中最大受力时的应力云图和位移云图。固定支架玉米田中受力在6種作物中最小，但此时位移云图中最大处已达17.98 mm，可见支架作业过程中位移变形量较大，需后续进行模态分析。

2.3 优化设计

考虑到原模型应力集中区出现在平行四边形结构右上角喷杆连接处，设计如图5-a所示毛坯件来进行优化设计，以弥补原模型中的不足。图5-a在原模型基础上增添了A区域和C区域，厚度与原模型相同，为30 mm。在Hyperworks软件中对待优化毛坯件进行应力分析，材料属性及固定约束设置与原模型相同，加载B=0.8 m时大豆田中作业的受力工况。毛坯件受力云图及位移云图如图5-b、 图5-c所示，毛

坯件最大应力为50.22 MPa，最大位移为0.23 mm，最大应力处未达到235 MPa，可进行拓扑优化设计。

在Hyperworks软件Topology模块中以最小化零件材料体积分数为优化目标，应力值不超过235 MPa为约束条件，对待优化毛坯件进行拓扑优化设计。经80步迭代计算后，得出如图6-a所示的优化密度图。从图6可以看出，优化件相比于原模型，在平行四边形结构内部增加了若干加强筋，而且在与喷杆连接处上方增加了1条加强筋。在Hyperworks软件Post选项板中的OSSmooth模块将图6拓扑优化模型以*.iges曲面模型格式导出。在PTC/Creo中将上述优化曲面模型导入，参照导入曲面模型轮廓进行三维绘图设计，优化模型设计完成后删除原曲面模型并导出二维机械图纸。考虑到加工难易程度，加工后的零件结构如图6-b所示，质量为1.87 kg。

进一步对优化设计模型进行应力分析，材料属性、固定约束设置与原模型相同，继续加载B=0.8 m时大豆田中作业的受力工况，应力云图和位移云图如图6-c、图6-d所示。图6显示，最大应力为123.8 MPa，最大位移为3.9 mm，满足设计要求。

2.4 模态分析

考虑到拨禾板固定支架原模型在作业过程中位移变形量较大，对原模型和优化模型进行前5阶模态分析，以对比支架优化前后模态频率变化情况，如表3所示。从表3可以看出，优化后支架模态频率皆有大幅度提升，即支架的抗振性增强。

3 讨论与结论

本研究运用机械有限元仿真软件Hyperworks，针对生产中遇到的实际问题对拨禾板固定支架进行结构分析和拓扑优化设计。以大田作业受力工况入手，计算玉米、高粱、小麦、水稻、大豆、油菜6种作物单位面积抗折力，经过对固定支架原模型应力分析得出应力集中部位，并在此基础上建立优化模型毛坯件，通过计算机拓扑优化设计得出固定支架优化模型。原固定支架不能在密集的秸秆作物中作业，尤其是茎秆抗折力大的作物冠层中。以B=0.8 m时抗折力最大的大豆为例，拓扑优化后单个固定支架质量虽然增至1.87 kg，比原支架增加了58.5%，但最大应力由原来的652.4 MPa降至 123.8 MPa，降了81%，机械抗振性也有大幅度的增强。本研究将工业机械设计中常用的有限元理论应用在农业机械上，不仅增强了材料的有效利用率，而且缩短了开发设计周期，具有重要的现实意义。

参考文献：

[1]董云哲，李君兴，史云天，等. 植保无人机未来发展的思考[J]. 陕西农业科学，2016，62（3）：119-121.

[2]He X K，Bonds J，Herbst A，et al. Recent development of unmanned aerial vehicle for plant protection in East Asia[J]. International Journal of Agricultural and Biological Engineering，2017，10（3）：18-30.

[3]贾卫东，申 彬，周慧涛，等. 不同侧风和风幕风速对风幕式喷杆喷雾飘移的影响[J]. 农机化研究，2018（7）：62-67.

[4]张 京，李 伟，宋坚利，等. 挡板导流式喷雾机的防飘性能试验[J]. 农业工程学报，2008，24（5）：140-142.

[5]王 磊，李雅文，樊廷录，等. 适宜机械收获株行距对黄土旱塬春玉米产量及水分利用效率的影响[J]. 水土保持研究，2017，24（5）：363-370.

[6]张永科，黄文浩，何仲阳，等. 玉米密植栽培技术研究[J]. 西北农业学报，2004，13（4）：98-103.

[7]边大红，刘梦星，牛海峰，等. 施氮时期对黄淮海平原夏玉米茎秆发育及倒伏的影响[J]. 中国农业科学，2017，50（12）：2294-2304.

[8]沈 静. 甜高粱栽培技术要点[J]. 中国农业信息，2016（2）：79-80.

[9]崔凤娟，王振国，李 岩，等. 高粱茎秆性状及倒伏系数的研究[J]. 作物杂志，2014（2）：61-64.

[10]董海滨，张 煜，许为钢，等. 高产小麦品种蘖叶构型动态模式的探索[J]. 麦类作物学报，2017，37（12）：1555-1563.

[11]姚海坡，董志强，吕丽华，等. 不同植物生长调节剂对冬小麦茎秆特性和产量的影响[J]. 华北农学报，2015（增刊1）：152-156.

[12]秦 炎，秦亚平. 种植密度对水稻群体质量和产量形成的影响[J]. 湖北农业科学，2017，56（12）：2222-2225.

[13]肖 楠，董立强，丛 琳，等. 施氮量对旱直播粳稻抗倒伏性及产量的影响[J]. 沈阳农业大学学報，2017，48（6）：647-653.

[14]史 宏. 种植密度对机械化收获大豆底荚高度的影响[J]. 农学学报，2017，7（9）：62-67.

[15]黄中文，赵团结，喻德跃，等. 大豆抗倒伏性的评价指标及其QTL分析[J]. 作物学报，2008，34（4）：605-611.

[16]吴永成，徐亚丽，彭海浪，等. 播期及种植密度对直播油菜农艺性状和产量品质的影响[J]. 西南农业学报，2015，28（2）：534-538.

[17]陈新军，戚存扣，浦惠明，等. 甘蓝型油菜抗倒性评价及抗倒性与株型结构的关系[J]. 中国油料作物学报，2007，29（1）：54-57，62.

[18]毛 潭，冯志彬，尹宇鹤，等. 基于近似线性弹性理论的太阳能设备支架稳性分析方法[J]. 机械工程师，2015（2）：41-45.

[19]陈树勋，韦齐峰，黄锦成. 利用导重法进行结构拓扑优化[J]. 计算力学学报，2015，32（2）：160-166.

[20]杨 喜，王金丽，郭昌进，等. 基于Workbench的甘蔗叶粉碎机刀辊模态分析[J]. 中国农机化学报，2014，35（4）：36-39.