烟苗移栽机机架轻量化研究_李世成.pdf

第 61 卷 第 3 期Vol.61 No.32023 年 3 月March 2023农业装备与车辆工程AGRICULTURAL EQUIPMENT&VEHICLE ENGINEERING0 引言目前的烟苗移栽机不能满足新的移栽技术要求，基本靠人工作业，导致劳动密集、劳动强度大、移栽质量不稳定、作业机械化程度低、效率低及劳动成本高等。为解决上述问题，课题组研制了一款可实现垄上打孔、投苗、施水等一体化作业的烟苗移栽机1。但该移栽机存在结构复杂、多种回转部件工作时振动较大、整机质量较大等一系列问题。机架是移栽机各总成装置的安装基体，在复杂工况载荷的激励作用下会导致机架发生塑性变形，从而影响移栽机工作性能。因此，机架轻量化设计对提升烟苗移栽机性能、节省能源有着重要作用2。以课题组自行研制的烟苗移栽机机架为研究对象，为实现移栽机正常工作前提下减少整体质量，基于 HyperMesh 建立烟苗移栽机机架的动静态力学有限元模型，采用响应面近似模型和多目标优化算法对机架各组件的厚度进行优化，在确保烟苗移栽机良好作业性能的条件下，避免机架产生共振，达到机架轻量化的目的。1 机架有限元模型的建立烟苗移栽机主要由机架、移栽装置、行走装置、控制装置、供能装置等构成，其中机架主要采用碳素结构钢Q235A的矩形管和钢板等型材焊接制成，图 1 为烟苗移栽机的三维结构模型。机架整体结构主要是由薄壁钢管型材焊接加工而成，对三维模型进行中面抽取得到二维中面组合而成的机架，再对中面赋予不同厚度尺寸代替原始三维机架模型；忽略一些非承载部件及机架倒角等对计算结果影响不大的因素，设置网格单元尺寸为 10 mm，网格单元类型为 QUAD4 和 CTRIA3，得到了 49 775 个单元和 49 432 个节点，网格平均doi:10.3969/j.issn.1673-3142.2023.03.008烟苗移栽机机架轻量化研究李世成（550025 贵州省 贵阳市 贵州大学 机械工程学院）摘要 为减轻机架质量和分析其动态特性，针对课题组研制的烟苗移栽机机架进行轻量化研究，通过建立机架有限元模型，展开机架动静态力学特性校核与分析，得到机架各项均满足设计要求，且材料冗余过大，存在很大的优化空间。基于 HyperStudy 搭建了机架静态工况与模态工况下的机架多目标优化模型，利用多目标遗传算法进行求解化设计。结果表明：优化后的机架经优化后的质量减轻了 40.56%，机架第 4 阶模态频率为 46.08 Hz，避开了共振频率，机架最大应力为 64.87 MPa，减小了 27.17%，表明优化后的机架具备了良好的力学性能。关键词 移栽机；机架；轻量化；动态特性；多目标优化 中图分类号 S223.9 文献标志码 A 文章编号 1673-3142(2023)03-0037-06引用格式：李世成.烟苗移栽机机架轻量化研究 J.农业装备与车辆工程,2023,61(3):37-42.Research on lightweight frame of tobacco seedling transplanterLI Shicheng(College of Mechanical Engineering,Guizhou University,Guiyang 550025,Guizhou,China)Abstract To reduce frame mass and analyze its dynamic characteristics,develop a lightweight research of the tobacco seedling transplanting machine frame.Through the establishment of finite element model of a frame,check and analyze the dynamic and static mechanical properties of the frame and find the chassis frame all meet the design requirements,and the material redundancy is too large,so there is much room for optimization.Based on HyperStudy,the multi-objective optimization model of the frame under static and modal conditions was built and the multi-objective genetic algorithm was used to solve the design.The results show that the mass of the optimized frame is reduced by 40.56%,the fourth-order modal frequency of the frame is 46.08 Hz,avoiding the resonance frequency,and the maximum stress of the frame is 64.87 MPa,which is reduced by 27.17%,indicating that the optimized frame has good mechanical properties.Key words transplanter;frame;lightweight;dynamic characteristics;multi-objective optimization基金项目：黔科合平台人才 20195616 号“贵州省山地智能农业装备研发与应用科技创新人才团队项目”收稿日期：2022-01-06 38农业装备与车辆工程 2023 年质量为 99.4%，表明网格质量较好，确保求解结果的可靠性。网格质量检测结果如图 2 所示3。2 机架的静力学分析移栽机机架为整个机体各关键部件的安装基体，为避免机架产生塑性变形，要求设计足够的强度。机架承受的主要载荷来自左侧边的操作平台、右侧的蓄电池、前部的汽油发电机、后部的水肥箱以及中部的移栽平台。查阅各驱动部件的参数并计算得到各关键部件的载荷，如表 1 所示。表 1 弯曲工况下机架外部载荷状况Tab.1 External load condition of the frame under bending condition载荷名称质量/kg汽油发电机20蓄电池25水肥箱25移栽平台50操作台80Y 向驱动装置12机架材料选用 Q235，其屈服强度 s=235 MPa。机架主要承受来自机架上安装部件的重力载荷，取动载荷系数 Kd=2.5，计算载荷 Fc=KdKg，在各关键承载部件的质量中心建立质量点，并赋予其相应的质量参数，通过建立的 rbe3 单元均布施加在机架的支承位置4。利用 OptiStruct 求解器对机架有限元模型进行求解，得到相应的位移云图如图 3 所示，应力云图如图 4 所示。由图3可知，机架变形量1.61 mm，小于2.5 mm，对机架整体性能影响很小，表明机架初始设计方案满足刚度要求。从图 4 可以看出，机架应力主要集中分布在底梁、中层纵梁，机架的最大应力为 89.21 MPa。在复杂的工作环境下，设计机架时按静强度准则应采用提高安全系数的方法，取n=1.5，则机架的许用应力=s/n=156.7 MPa。综上，机架正常工作所产生的最大应力值和最大变形值均满足设计要求，且存在较大的优化空间。3 机架模态分析激振源频率分析有助于评价机架结构的动态力学性能5。所有类型电机和回转部件工作时，可通过 f=n/60 计算机架上各驱动电机和主要回转部件的理论激振频率，汇总后如表 2 所示。1.机架 2.履带驱动电机 3.汽油发电机 4.打孔投苗装置5.控制柜 6.苗盘支撑架 7.水肥箱 8.Y向驱动装置9.履带 10.蓄电池图 1 烟苗移栽机的三维结构模型Fig.1 Three-dimensional structure model of tobacco seedling transplanter11029387456图 2 烟苗移栽机机架有限元网格模型Fig.2 Finite element mesh model of tobacco seedling transplanter frameElements Vlolating Thresholds:min size =0 (0.0%)max size =0 (0.0%)aspect ratio =0 (0.0%)warpage =0 (0.0%)skew =0 (0.0%)jacobian =0 (0.0%)max angle quad =12(0.0%)min angle quad =0 (0.0%)max angle trla =0 (0.0%)min angle tria =0 (0.0%)taper =0 (0.0%)YZXContour PlotDisplacement(Mag)Analysis systemMax=1.614E+00 Grids 2307Min=0.000E+00 Grids 233781.614E+001.435E+001.255E+001.076E+008.966E-017.173E-015.380E-013.587E-011.793E-010.000E+00No Result图 3 弯曲工况下机架位移云图Fig.3 Frame displacement cloud diagram under bending condition Contour Plot Element Stresses Analysis systemMax=8.921E+012D 22875Min=1.364E-022D 602053.921E+013.485E+013.050E+012.614E+012.179E+011.743E+011.308E+018.723E+004.368E+001.364E-02No Result图 4 弯曲工况下机架应力云图Fig.4 Frame stress cloud under bending condition39第 61 卷第 3 期表 2 激振源理论振动频率Tab.2 Theoretical vibration frequencies of excitation sources名称转速/(r/min)理论激振频率/HzY 向丝杆驱动电机3 00050Z 向丝杆驱动电机1 00016.7成穴钻头3005履带驱动轮911.5机架的环境决定了其约束条件比较复杂，应当采取自由模态计算方式求解机架的模态参数。通常高阶模态对结构强度影响较小，只提取其前 6 阶自由模态频率与最大变形量即可。利用 HyperMesh对建立好的机架进行模态分析、模型计算，得到机架前 6 阶自由模态振型云图如图 5 所示。汇总后的机架最大变形量与模态频率如表 3 所示。由表 3 可以看出，机架的前 6 阶自由模态频率介于 28.8881.63 Hz，第 4 阶模态频率下的机架操作台装置安装支架变形量较大，机架有发生共振几率，可能导致移栽机的工作性能，应进一步优化机架模态固有频率6。表 3 机架前 6 阶非刚体模态频率与变形量Tab.3 The first six orders of non-rigid body modal frequency and deformation of the frame模态阶次模态频率/Hz最大变形量/mm128.885.618243.784.330347.258.192448.1318.880564.308.110681.636.8394 基于 HyperStudy 的机架轻量化研究4.1 确定设计变量机架结构主要由底梁、中层纵梁、