轻量化三七移栽机车架有限元分析及拓扑优化设计
1.
2.
3.
Finite element analysis and topology optimization design of lightweight Panax Notoginseng transplanting machine frame
1.
2.
3.
通讯作者: 苏微,女,1979年生,博士,副教授;研究方向为农业机械化生产和管理。E-mail:laisubo@163.com
收稿日期: 2024-03-01 修回日期: 2024-04-10
基金资助: |
|
Received: 2024-03-01 Revised: 2024-04-10
作者简介 About authors
马遥,男,1999年生,硕士研究生;研究方向为农业机械装备与计算机测控。E-mail:1348124520@qq.com
三七是我国种植规模最大,开发利用较多的大宗药材之一,而云南省作为我国三七主要的产区,由于地形和农艺要求的影响,传统农业机具难以进入移栽区域进行作业,导致目前三七移栽仍采用传统人工作业方式,为此,研制一款三七移栽机对于促进三七产业化发展尤为重要。车架作为移栽机车身重要的承载结构,对整车性能的优劣有着至关重要的影响。为此,针对移栽机车架进行结构分析,旨在提升车架结构及整车的性能,为车架结构设计提供理论依据。利用SolidWorks软件建立车架三维模型,并导入ANSYS软件中进行静态有限元分析,通过静态电测试验确定有限元分析的应力值和试验真实应力值的相对误差范围后,对不同工况下车架动载性能进行分析以及前8阶模态振动变形分析,基于分析结果可知车架具有良好强度性能,但座椅位置存在较大变形,结构的刚度不足的问题。采用拓扑优化对车架进行优化设计,在保证车架应力分布合理的情况下减少车架变形情况通过改变斜撑支架结构布置实现减小变形的目的。根据优化设计结果,使车架总质量增加8.739%,变形量减小88.268%,最大应力减小11.693%。改进后的车架最大应力有所降低,并且车架变形有何明显改善,研究表明通过有限元法以及拓扑优化技术对于指导移栽机车架的结构设计具有较好的适用性。
关键词:
Panax Notoginseng is one of the most extensively cultivated and utilized bulk medicinal materials in China, and Yunnan Province is the main producing area for Panax Notoginseng in the country. Due to the influence of terrain and agronomic requirements, traditional agricultural machinery faces difficulties in entering the planting areas for operation, resulting in the current manual operation of Panax Notoginseng transplantation. Therefore, the development of a Panax Notoginseng transplanting machine is crucial for promoting the industrialization of Panax Notoginseng. As the crucial load-bearing structure of the transplanting machine, the frame significantly affects the overall performance of the vehicle. This paper focuses on the structural analysis of the transplanting machine frame, aiming to improve the performance of the frame and the entire vehicle, and provide theoretical basis for frame structural design. A three-dimensional model of the frame is established using SolidWorks software and imported into ANSYS software for static finite element analysis. After determining the relative error range between the stress values obtained from static electrical tests and experimental stress values, the dynamic loading performance of the frame under different working conditions and the modal vibration deformation analysis of the first eight modes are conducted. The analysis results indicate that the frame exhibits good strength performance, but significant deformation occurs at the seat position, indicating insufficient stiffness of the structure. Topology optimization is employed to optimize the frame design, aiming to reduce frame deformation while ensuring reasonable stress distribution. By altering the arrangement of diagonal brace brackets, the goal of reducing deformation is achieved. According to the optimization design results, the total mass of the frame increases by 8.739%, while the deformation is reduced by 88.268%, and the maximum stress is decreased by 11.693%. The improved frame exhibits reduced maximum stress and significantly improved deformation, demonstrating the applicability of finite element method and topology optimization technology in guiding the structural design of transplanting machine frames.
Keywords:
本文引用格式
马遥, 苏微, 赖庆辉, 张弦, 于庆旭, 宛渊.
MA Yao, SU Wei, LAI Qinghui, ZHANG Xian, YU Qingxu, WAN Yuan.
0 引言
目前三七移栽作业面临着无机可用的环境,因此研制一种适用于丘陵山地的自走式轻量化三七移栽机具有十分重要的意义。当前国内针对丘陵山地移栽机结构件的设计和研究主要集中在栽植机构,大多数移栽机以牵引式为主,对于自走式移栽机车架的研究较少。而车架是整机的重要支撑结构,是保障移栽作业质量的关键部件,它不仅承载移栽机各部件,还要承受外界对移栽机的各种载荷冲击[5],更是决定了三七移栽作业的工作效率与安全。
针对这一现状,笔者通过有限元法结合拓扑优化方法对移栽机车架进行分析及优化,通过应力试验验证仿真结果准确性,并根据优化结果提出车架结构改进方案。
1 车架有限元模型建立
1.1 基于SolidWorks软件的车架三维模型建立
本文研究对象为昆明理工大学自主研发的丘陵山地轻量化三七移栽机,如图1所示。
图1
图1
三七轻量化移栽机
1.电池总成
2.传动装置 3.电机 4.后轮 5.座椅位置
6.移栽装置 7.前轮 8.机架
Figure 1
Panax Notoginseng lightweight transplanting
该移栽机采用直流减速电机作为动力源,在作业时通过传动装置驱动移栽机前进。操作人员将三七种苗放置在移栽装置输送带上,输送带的速度与移栽机的行进速度相匹配,从而实现随速控制。以确保三七移栽株的行距符合规范,机具质量为85 kg,换行作业时可由人工抬起转向。
该移栽机质量较小,但其车架强度与刚度没有足够的数据支撑,因此,通过对车架结构进行有限元分析,探究车架结构性能薄弱点,从而进行优化改进。车架长×宽×高为:1 600 mm×830 mm×420 mm,主要由45Mn方钢(20方钢厚2 mm)和钢板(厚5 mm)焊接而成。利用SolidWorks软件按照实际测量的车架尺寸绘制整体车架三维图。
1.2 基于SolidWorks软件的车架有限元模型建立
图2
图2
移栽机车架几何模型
1.前轮支撑横梁
2.移栽装置总成支撑横梁 3.踏板
4.踏板支撑纵梁
5.座椅总成支撑横梁 6.电池总成支撑纵梁
7.电机支撑横梁
8.传动总成支撑纵梁 9.竖直支撑梁
10.斜支撑梁
11.纵梁 12.前轮支撑板
Figure 2
Geometric model of transplanting frame
表1 45Mn碳素钢材料参数
Table 1
材料参数 | 数值 |
---|---|
拉伸强度/MPa | 600 |
屈服强度/MPa | 355 |
弹性模量/MPa | 2.1×105 |
泊松比 | 0.269 |
密度/(kg·m-3) | 7.85×103 |
考虑到移栽机车架的结构形状和受力情况,在ANSYS软件中采用壳单元对车架进行网格划分。选取网格划分尺寸为5 mm,车架模型共划分415 512个节点,207 072个单元。
2 移栽机车架结构静态与模态分析
2.1 载荷及工况约束处理
2.1.1 载荷处理
表2 有限元分析中各载荷数值
Table 2
序号 | 载荷名称 | 载荷数值/N | 施加位置 |
---|---|---|---|
1 | 直流减速电机 | 137.2 | 电机支撑横梁 |
2 | 第一级传动装置总成 | 147 | 传动总成支撑纵梁 |
3 | 驾驶员与座椅总成 | 686 | 座椅总成支撑横梁 |
4 | 移栽装置总成 | 392 | 移栽装置总成支撑横梁 |
5 | 铅酸电池 | 137.4 | 电池总成支撑纵梁 |
6 | 种苗箱及种苗 | 98 | 踏板支撑纵梁 |
2.1.2 约束与工况分析
三七移栽机属于农用车辆,作业环境为丘陵山地等复杂地形,移栽机主要在畦床上行驶,移栽前对畦床进行平整,土质松软。本文所提及三七移栽机整机质量为85 kg,换行作业时由人工抬起进行转向,不考虑转向行驶。因此确定以下两种工况对其进行应力分析[10]。
图3
图3
移栽机车架动载荷边界条件处理
Figure 3
Processing of dynamic load boundary conditions for transplanting frames
2.2 车架结构有限元分析
在两种不同的工况条件下,对车架模型进行了有限元分析,并得到了以下结果。
图4
图4
满载弯曲工况下应力与变形分布图
Figure 4
Stress and deformation distribution diagram of full load bending conditions
2.3 移栽机车架强度静态电测试验
采用试验验证有限元模型载荷施加和边界条件设置的正确性。根据电阻应变原理[14],选择箔式电阻应变片、DCS-100A数据采集软件和PCD-300B数据分析仪对数据进行收集和处理。
在有限元分析的基础上,为了更精确地评估车架模型与实际结构关于应力分布和变化准确性,选取了一批应力变化最为显著的关键位置作为静态电测试验的测量点,共计选取12个测量点,车架应力测点位置分布如图6所示。根据试验结果,整理出表3,涵盖了车架上12个测试点的试验值和有限元分析值数据。根据对比结果显示,除了测点10和12的相对误差超过20%外,其余测点的相对误差均在20%以内。这进一步证明了所采用的有限元模型的简化、约束的确定、载荷的简化以及处理方法基本正确[15]。个别测点的测试值与应力值存在偏差的原因为:(1)边界条件未能准确地反映实际情况;(2)载荷施加的加载条件可能与实际测试中的加载条件存在差异。
图5
图5
满载扭转工况下应力与变形分布图
Figure 5
Stress and deformation distribution diagram under full load torsion conditions
图6
图6
车架应力测点位置分布示意图
Figure 6
Stress distribution of observation points
表3 静态应变试验结果与有限元分析结果对比
Table 3
测点号 | 试验值/MPa | 计算值/MPa | 相对误差/% |
---|---|---|---|
1 | 0.998 53 | 1.137 9 | -13.958 |
2 | 49.628 | 44.741 | 9.847 3 |
3 | 22.887 | 24.294 | -6.147 6 |
4 | 97.427 | 91.561 | 6.020 9 |
5 | 14.869 | 13.731 | 7.653 5 |
6 | 2.989 2 | 3.412 3 | -14.15 4 |
7 | 2.482 | 2.361 1 | 4.871 1 |
8 | 38.992 | 40.778 | -4.580 4 |
9 | 103.35 | 93.456 | 9.573 3 |
10 | 25.237 | 36.308 | -43.868 |
11 | 23.181 | 24.181 | -4.313 9 |
12 | 43.877 | 53.888 | -22.816 |
2.4 移栽机车架模态分析
在保持线性不变的前提下,对于一个连续的且非线性的物理结构,经过离散化处理后,转化为一个具有无限维度(当N趋于无穷时)的线性弹性系统。该系统的动力学特性可以由微分方程准确描述,如
式中:[M]——质量矩阵;
[C]——阻尼矩阵;
[K]——刚度矩阵;
x(t)——位移响应向量;
f(t)——激励力向量。
本文所探讨的是一个基于多自由度、无阻尼的自由振动系统模型。这种无阻尼的自由振动系统可更清晰地分析系统内部的动态响应,进而为理解更复杂系统的振动特性提供基础。因此,上述微分方程可表达为
在有限元法中,模态分析的核心理论在于利用坐标变换方式,将描述物理系统动态的原始物理坐标空间转化为模态坐标。并对该模型解耦,使耦合方程组变成n个独立的微分方程,再对它们逐一进行求解[20]。将物理坐标{x(t)}用模态坐标[u]Ф表示,如
对无阻尼系统和比例阻尼系统来讲,振型矩阵u具有模态正交性[20],通过将振型矩阵u进行正交化处理,可以将系统的运动方程分解为一系列独立的振动模态,如
解上述联立方程式,可解得:
式中:A,φ——代表任意常数;
{u}——模态向量。
将
其中:
对模态向量进行正规化处理,将模态向量的长度标准化为1,则
将n个正规化的模态向量,按顺序排列成一个方阵,构成模态矩阵[u],由上述论述的变化过程,可知方程(7)可以写成:
在经过上述一系列步骤后,将无阻尼的自由振动系统微分方程转换成了解耦方程组。解耦方程组的形式简化了问题的复杂性,便于进一步分析和求解。运输机车架的前8阶振形图如图7所示。
图7
图7
移栽机车架的前8阶模态振形图
Figure 7
First 8 order model vibration shape of transport plane frame
表4 车架结构的各阶模态频率表
Table 4
阶数 | 频率f/Hz | 振型特征 |
---|---|---|
1 | 16.085 | 二阶绕z轴弯曲振动 |
2 | 17.597 | 扭转振动 |
3 | 29.554 | 绕x弯曲振动 |
4 | 40.022 | 二阶扭转振动 |
5 | 42.237 | 绕z轴弯曲振动 |
6 | 55.044 | 扭转振动 |
7 | 70.228 | 绕x弯曲振动 |
8 | 77.446 | 弯扭组合振动 |
传统的汽车发动机是汽车振动主要的激励[22],本研究将传统发动机替换为直流减速电机作为移栽机的动力源。由于电机平衡性较好,因此有效减少了整车的振动激励。移栽机振动激励主要有3部分:
1) 道路表面的不平整造成的位移激励;
2) 主减速器及其后续的传动系统的激励;
3) 行驶系统的转动部件的转动不平衡量引起的额外激励[23]。
对于道路激励频率一般小于3 Hz,车轮不平衡所引起的激励一般小于5 Hz[24]。
根据车架的分析结果可知,车架的前三阶的振动频率为16.085~29.554 Hz之间,与移栽机在正常行驶过程中可能遭受的主要激励频率进行了对比。有效地避开了那些可能导致车架与整车发生共振的频率区间,从而显著提升了车架在实际应用中的稳定性和安全性。
3 移栽机车架拓扑优化设计
结构优化是指在满足设计要求和约束条件的前提下,通过调整结构的形状、尺寸、材料等,使得结构在性能、重量、成本等方面达到最优。结构优化主要包括以下几种类型:
1) 尺寸优化:尺寸优化是通过调整结构的尺寸参数(如梁的截面尺寸、板的厚度等)来达到最优设计。其目标通常是材料使用量的最小化。
2) 形状优化:形状优化是通过改变结构的几何形状,使得结构在力学性能(如应力分布、刚度等)和其他性能指标(如气动性能、流体阻力等)上达到最优。
3) 材料优化:材料优化是通过选择和组合不同材料,使得复合材料结构在性能和成本上达到最优。
4) 拓扑优化:拓扑优化是在一定空间区域(骨架结构或连续体)内寻求材料最合理分布的一种优化方法,它是一个迭代的过程。其目的是寻求结构的刚度在设计空间内最佳的分布形式或结构的最佳传力路线,以优化结构的某些性能或减轻结构的重量[25]。
图8
根据拓扑优化结果,考虑各总成的布置和实际作业要求,改变车架斜支撑梁布置方式,通过延长斜支撑梁至前车轮处,并添加两根竖直支撑梁,建立如图10所示的新车架结构优化方案。
图9
图10
将优化后的设计方案导入到ANSYS中,分别进行满载弯曲工况和扭转工况的有限元分析。载荷和约束方式与原车架相同。改进后车架结构在两种工况下的等效应力云图和变形云图见图11,车架在弯曲工况下的最大等效应力为134.87 MPa,最大变形量为1.264 mm,满足了汽车设计的要求。而在扭转工况下,车架的最大等效应力为209.79 MPa,最大变形量为4.287 5 mm。
图11
图11
结构改进后车架的等效应力云图和变形云图
Figure 11
Equivalent stress and displacement of improved frame
表明车架的强度和刚度相比优化前均有所提升,这对于增强车辆的稳定性和安全性至关重要。得出结果与原车架分析结果对比如表5所示。
表5 车架模型优化前后各主要性能参数对比值
Table 5
变量名称 | 优化前 | 优化后 | 变化情况/% |
---|---|---|---|
车架总质量/kg | 19.910 | 21.650 | 增加8.739 |
弯曲最大等效应力/MPa | 154.11 | 134.87 | 减少12.484 |
弯曲最大变形量/mm | 10.774 | 1.264 | 减少88.268 |
扭转最大等效应力/MPa | 237.57 | 209.79 | 减少11.693 |
扭转最大变形量/mm | 11.599 | 4.287 5 | 减少63.035 |
一阶扭转频率/Hz | 17.597 | 29.131 | 增加65.545 |
一阶弯曲频率/Hz | 29.554 | 41.433 | 增加40.194 |
从表5可知,优化后的车架质量增加了8.739%,通过对运输机上车架的关键部位进行优化,扭转和弯曲工况下最大等效应力有一定程度上的减小,说明优化后车架的整体最大变形量有着显著的变化,弯曲、扭转最大变形量分别为1.264 mm,4.287 5 mm,均满足汽车设计理论要求。
此外,优化后低阶模态频率为29.131 Hz以上,更好避免了由振动激励引起的车架共振的发生。车架各梁材料为45Mn,屈服极限为355 MPa,安全系数较高,说明优化后车架的整体性能得到很大的提升。整体而言,采用这个优化方案可以更好地改善移栽机车架的刚度不足的情况。在保证车架安全性能的情况下,实现车架轻量化的目的。也意味着车架在面对各种复杂的工作环境时,将表现出更高的稳定性和可靠性。
4 结论
1) 采用了成熟的有限元理论,并利用ANSYS软件,对移栽机车架的力学特性进行了深入分析。建立了车架的有限元模型,并对其在弯曲和扭转工况下的静态和动态特性进行了详尽研究。这项研究为车架结构试验和优化设计提供了可靠的理论基础。在验证了模型准确性后,我们对车架进行了拓扑优化,并提出了优化方案。根据分析结果显示,优化后的车架整体性能得到了显著提升。
2) 根据有限元分析结果,本文重点关注原车架存在的严重变形问题,并采用拓扑优化方法对其进行了设计优化,优化后,移栽机车架的变形量大幅减少了88.268%,为1.264 mm;最大应力减小了11.693%,为134.87 MPa;车架质量增加了8.739%,为21.65 kg。在保障车架整体性能的情况下,将车架质量控制在可接受范围,实现了优化目标。
3) 在丘陵山地等复杂地形条件下,移栽机需要承受更大的载荷和更频繁的振动。通过轻量化设计,可以有效降低车架的自重,减轻整机负担,从而提高移栽机的稳定性和操控性。同时,轻量化车架的研究还可以推动农机行业的技术进步和产业升级,促进农机行业的可持续发展。
参考文献
导苗管式三七种苗栽植机构参数优化试验
[J]. ,
Parameter optimization experiment of seedling guiding tube transplanting machine of Panax notoginseng seedling
[J]. ,
构建农机农艺融合的工程技术体系, 促进云南三七产业可持续发展
[C]//
旱地移栽机械取苗机构发展现状与趋势
[J]. ,
Development status and trends of the development of seedling pick-up mechanism of transplanter in dry land
[J]. ,
浅山丘陵作业车辆底盘关键部件的设计与分析研究
[D].
Design and analysis of key components of chassis for light-duty vehicles in shallow mountainous terrain
[D].
Vegetable transplanters for use in developing countries—A review
[J]. ,
机动式喷杆喷雾机机架的轻量化设计
[J]. ,
Lightweight design of chassis frame for motor boom sprayer
[J]. ,
山地果园轮式运输机车架结构分析与优化
[J]. ,
Analysis and optimization of frame structure for wheeled transporter in hill orchard
[J]. ,
车辆与路面相互作用下路面结构动力分析
[D].
Research on dynamic response of pavement structure due to vehicle and pavement interaction
[D].
车辆对路面作用的动载荷研究
[D].
Study on dynamic road of vehicles on pavement
[D].
标准四驱水田自走底盘转向驱动桥设计与工程结构分析
[J]. ,
Design and engineering structure analysis on paddy-field chassis steering drive axle with standard all-wheel-drive
[J]. ,
重型载货汽车车架有限元分析及拓扑优化
[D].
Finite elements analysis and topological optimization of heavy-duty truck's frame
[D].
基于有限元法的某型大学生方程式赛车车架优化设计
[J]. ,
Optimization of a FSAE car frame based on the finite element method
[J]. ,
基于振动特性的联合收获机车架的结构优化
[D].
Structural optimization of combine harvester frame based on vibration characteristics
[D].
客车骨架结构分析及优化研究
[D].
Research on the analysis and optimization of the coach-body frame
[D].
基于ANSYS的除草机车架轻量化研究
[J]. ,
Lightweight research of mower frame based on ANSYS
[J]. ,
基于发动机安装位置的货车车架减振优化
[J]. ,
Optimization of vibration of truck's frame based on engine mounting location
[J]. ,
多功能电动车车架结构有限元分析及优化
[D].
The finite element analysis and optimization of a multi-function electric vehicle frame
[D].
高速水稻插秧机车架的轻量化设计
[J]. ,
Light design of frame for self-propelled chassis rice transplanter
[J]. ,
基于Hyperworks和Ansys的汽车车架有限元分析
[D].
FEM analysis of truck frame model based on Hyperworks and Ansys
[D].
/
〈 | 〉 |