基于离散元法的旋耕刀三向工作阻力仿真分析与试验.PDF
《基于离散元法的旋耕刀三向工作阻力仿真分析与试验.PDF》由会员分享,可在线阅读,更多相关《基于离散元法的旋耕刀三向工作阻力仿真分析与试验.PDF(9页珍藏版)》请在得力文库 - 分享文档赚钱的网站上搜索。
1、第34卷 第18期 农 业 工 程 学 报 Vol.34 No.18 2018年 9月 Transactions of the Chinese Society of Agricultural Engineering Sep. 2018 113 基于离散元法的旋耕刀三向工作阻力仿真分析与试验 熊平原1,2,杨 洲1,孙志全1,张倩倩1,黄杨清1,张卓伟1 (1. 华南农业大学工程学院,广州 510642;2. 仲恺农业工程学院机电工程学院,广州 510225) 摘 要:为分析旋耕刀所受三向工作阻力及其变化规律,该文通过实测南方果园土壤颗粒参数,逆向重构旋耕刀三维实体,基于离散元颗粒接触理论,构建
2、了适应南方土质环境的旋耕刀-土壤相互作用仿真模型。土槽扭矩对比试验表明,仿真值与试验值的变化趋势相同,扭矩均随转速增加而变大,最大相对误差10%;扭矩先从0增加到某个最大值,接着逐步减小到一个低值,随后又快速增加到一个高值,最后回落,该变化过程同旋耕刀与土壤之间的接触状态相关。单刀受力仿真分析表明,水平阻力方向与前进方向相同,侧向阻力方向为由刀具弯折区内侧面指向刀体,垂直阻力方向为先垂直土面向上后转为向下;水平阻力和侧向阻力在最大耕深处出现最大值,而垂直阻力在入土后转动约30时出现最大值;水平阻力和垂直阻力的仿真波形与理论计算、土槽试验结果比照表明,对应曲线的变化趋势基本一致,且仿真结果与土槽
3、试验结果更为接近,水平阻力相对误差为11.3%,垂直阻力相对误差为16.8%;水平阻力最大值大于侧向和垂直方向阻力最大值,水平阻力是功率消耗的主要因素;随着转速的增加,3个方向阻力最大值均增大,当转速高于250 r/min时,增速加快;侧向阻力和垂直阻力随前进速度增加而平稳增大,水平阻力却出现下滑趋势;耕深对三向阻力的影响比较显著,增加耕深会急剧增大三向阻力值。相关试验数据可为旋耕机能耗分析、机体作业振动及刀片磨损等研究提供参考。 关键词:农业机械;土壤;计算机仿真;旋耕刀;三向阻力;离散元法 doi:10.11975/j.issn.1002-6819.2018.18.014 中图分类号:S2
4、22.3 文献标志码:A 文章编号:1002-6819(2018)-18-0113-09 熊平原,杨 洲,孙志全,张倩倩,黄杨清,张卓伟. 基于离散元法的旋耕刀三向工作阻力仿真分析与试验J. 农业工程学报,2018,34(18):113121. doi:10.11975/j.issn.1002-6819.2018.18.014 http:/www.tcsae.org Xiong Pingyuan, Yang Zhou, Sun Zhiquan, Zhang Qianqian, Huang Yangqing, Zhang Zhuowei. Simulation analysis and expe
5、riment for three-axis working resistances of rotary blade based on discrete element methodJ. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(18): 113121. (in Chinese with English abstract) doi:10.11975/j.issn.1002-6819.2018.18.014 http:/www.tcsae.
6、org 0 引 言 手扶式旋耕机由于体型轻巧、操作方便,在南方丘陵山区被广泛使用,主要用于果园、菜地和设施农业等小田块的耕整作业1-4。旋耕刀片是旋耕机的关键零部件,一方面绕主轴定速旋转,另一方面随机具匀速前进。刀体与土壤直接接触,周期性地完成入土、切土和抛土动作,土壤对刀体具有强烈的反作用力,表现出机具总功耗大、振动强和刀面磨损快5-9。中国普遍采用C型旋耕弯刀,正反刀面均为空间曲面,土壤对刀体的作用力系组成复杂,且随着刀体在土壤中的位置变化而改变。为便于研究,可将该复合作用反力分解成水平阻力、垂直阻力和侧向阻力,水平阻力即牵引阻力,是机具功耗的主要组成部分,垂直阻力和侧向阻力影响着机体振动
7、和刀片磨损状况,分析旋耕刀所受3个方向工作阻力值及其变化规律,是进行旋耕机能耗分析、机体作业振动研究及刀片减摩优化的前期基础10-12。 收稿日期:2018-04-26 修订日期:2018-07-31 基金项目:广东省现代农业产业技术体系创新团队(2017LM2153);国家现代农业产业技术体系建设专项资金(CARS-31) 作者简介:熊平原,男,湖北大悟人,博士生,主要从事农业机械化研究。Email: 通信作者:杨 洲,男,山西襄汾人,教授,博士生导师,主要从事农业机械化研究。Email: 目前国内外学者对旋耕刀工作阻力研究主要集中在 牵引阻力和功耗分析,对刀体所受三向阻力的变化规律研究较少
8、。Roul 等13设计出一种基于反向传播学习的神经网络模型,只需输入刀具宽度、耕深、运行速度、土壤含水率和紧实度等参数,就能够准确预测出不同类型的土壤耕作机具在砂质黏土中所需牵引力,试验表明该模型对铧式犁、中耕机及圆盘耙的牵引阻力预测误差均低于6.5%,比美国农业与生物协会给出的预测结果有很大提高14。Ahmadi 等15运用经典力学理论,从动力学视角推导出旋耕刀片动态转矩计算模型,转矩计算结果与Chertkiattipol 等16完成的试验结果相比,误差仅为5%。 Matin 等17分析了C型弯刀、半宽刀及直刀在四种转速下的扭矩、功率和能耗特性,并借助高速摄像机记录旋耕刀从入土到出土的全过程
9、。王增辉等18采用正交试验法,确定了影响旋耕碎茬单刀扭矩的主要试验因素,构建了功耗与刀片转速、耕深之间的函数关系。蒋建东等19提出了一种对旋耕机施加振动载荷以减小切削阻力的方法,并研究了外加激励的振型、频率和振幅对水平阻力的影响规律。方会敏等20构建了秸秆-土壤-旋耕刀相互作用模型,从土壤和秸秆运动及旋耕刀受力角度对秸秆、土壤和刀具的相互作用进行了分析。 已有研究仅给出了旋耕刀水平阻力受试验参数的影响规律及水平阻力值的预测模型,而对垂直阻力和侧向阻力的变化规律及3个方向阻力值的对比分析研究不足。本课题组前期通过改变旋耕刀在土壤中的相位角,研究农业工程学报(http:/www.tcsae.org
10、) 2018年 114 了三向工作阻力与各试验因素之间的变化规律,得到试验因素对单位幅宽阻力的影响等级21,但由于测量阻力值时,旋耕刀只有前进运动而不能随刀轴旋转,测量结果与实际受力情况存在偏差。因此,准确描述旋耕刀与土壤之间的相互作用状态,并实时获取刀片所受三向阻力值,成为急需解决的问题。旋耕刀工作时,刀体深入土层内部运动,直接用三维传感器实时检测刀具所受三向阻力值显得十分困难22-24。文献25-27研究表明,离散元法在土壤散体颗粒领域具有很好的应用前景。本文首先试验测量南方果园土壤颗粒特性参数,采用逆向设计法进行旋耕刀三维造型,运用离散元软件EDEM构建旋耕刀-土壤相互作用仿真模型;接着
11、进行旋耕扭矩仿真试验和土槽试验,对比试验结果以确定仿真模型的可靠性;最后仿真分析各试验因素对旋耕刀所受三向阻力的影响规律,并通过土槽试验辅助验证仿真分析结果的正确性。本研究构建了适应南方土质环境的旋耕刀-土壤相互作用模型,通过设定土壤颗粒黏结半径以提高模型表达精度。旋耕刀三向工作阻力仿真分析结果可为旋耕机能耗分析、机体振动和刀片磨损等研究提供数据参考。 1 模型构建 1.1 土壤颗粒参数 土壤物理和力学性质影响着旋耕刀-土壤力学模型相互作用效果28,为准确模拟南方丘陵土质环境,保证仿真结果可靠性,选取广东省农业科学院果树研究所试验地黏性红土进行参数测量。2016年3月8日实测0300 mm深度
12、土层平均含水为16.8%,平均紧实度为355 kPa;通过环刀法试验,测得土壤颗粒湿密度为1.78103 kg/m3,干密度为1.55103 kg/m3;采用筛析法进行颗粒度分析,先将采集土样烘干至恒质量,研碎团结体,使用孔径为2.0、1.0、0.5、0.25、0.075 mm的细筛进行筛分,结果表明直径大于2 mm的颗粒占比9.1%,在0.52 mm区间的占比68%,在0.10.5 mm区间的占20.8%,其余占比2.1%,故颗粒直径可取区间平均值1.25 mm;文献29中定义了土壤颗粒的9种基本状, 本研究在电子显微镜下观察采集土样发现,土壤颗粒浑圆少孔,呈粒状,可统一简化为球体;通过土壤
13、直剪试验,测得土壤颗粒内摩擦角为18,内聚力为0.9105 Pa;通过土壤单轴压缩试验,测得法向临界应力为32 kPa,切向临界应力为16 kPa;依据文献30设定土壤剪切模量为1106 Pa,泊松比为0.38。 1.2 旋耕刀建模 旋耕刀是旋耕作业的主体受力部件,其外形尺寸影响着刀体与土壤之间的相互作用效果。C型旋耕弯刀刃口曲线和刀身曲面相对复杂,如采用三维软件直接建模,通常需要简化部分曲线或曲面,降低了模型精度。本研究采用逆向工程建模法,首先将IT245旋耕刀片(天津市津旋农机具有限公司生产)用砂纸打磨光亮,喷涂反光材料,放入LPX-600型三维激光扫描仪(日本罗兰公司)进行外形构造,导出
14、成*.stl格式文件,然后使用Geomagic Studio软件对残缺部位进行修复,最后导入到UG NX10实例化。由于文献31推荐旋耕刀侧切刃采用阿基米德曲线,为校验刀体模型准确性,选取模型侧切刃曲线与标准阿基米德曲线进行对比,计算出相关系数R为0.99,误差较小,表明采用逆向建模法构建的旋耕刀体三维模型精度较高。 1.3 颗粒黏结模型 南方土壤含水率高、黏性强,为描述该固有特性,宜采用Hertz-Mindlin with Bonding黏结模型对土壤颗粒进行约束32-33。这种约束能够阻止颗粒间法向和切向的相对运动,当达到最大法向和切向应力时,约束被破坏,两颗粒才被视为硬球对彼此产生作用。
15、设t0时刻颗粒被黏结,法向黏结力Fn、法向力矩Tn、切向黏结力Ft和切向力矩Tt均随时步从0开始按式(1)进行调整。 24212n n nt t tn n tt t nBBF S A tF S A tT SJ tJT S tA RJ R (1) 式中vn、vt 为土壤颗粒的法向和切向速度,m/s;Sn、St为土壤颗粒的法向和切向刚度,N/m;n、t为土壤颗粒的法向和切向角速度,rad/s;A为土壤颗粒间接触区域面积,m2;J 为土壤颗粒的惯性矩,m4;RB为土壤颗粒间黏结半径,m;t为时间步长,s。Tn、Tt为土壤颗粒受到的法向和切向力矩,Nm。 当某一时刻土壤颗粒间黏结所承受的法向和切向应力
16、超过最大应力值,黏结就被破坏,条件为: maxmax2n tBt n BF T RA JF TRA J (2) 式中max、max为土壤颗粒间黏结所能承受的最大法向应力和切向应力,Pa;Fn、Ft为土壤颗粒间法向和切向黏结力,N。该土壤黏结模型中,颗粒之间由黏结带形成约束,黏结半径应比颗粒的实际半径大。图1为颗粒O1,O2的黏结状态,中间充满黏结带,在约束力作用下,2个土壤颗粒的实际接触距离(R1B+R2B)大于理论接触距离(R1+ R2)。 注:R1B、R2B为土壤颗粒间形成的黏结半径,m;R1、R2为土壤颗粒的实际半径,m。 Note: R1B, R2B represent the bon
17、ding radiuses of soil particles, m. R1, R2 represent the actual radiuses of soil particles, m. 图1 土壤颗粒黏结状态图 Fig.1 Bonding state of soil particles 在一定体积的土壤中,土壤颗粒形成的黏结带可转化为由于水分的存在而产生的黏结力。由前述土壤颗粒第18期 熊平原等:基于离散元法的旋耕刀三向工作阻力仿真分析与试验 115 参数测定结果可知,约70%的粒径分布在某一特定缩小区间,且土壤颗粒外形浑圆,故所有土壤颗粒可简化为相同直径的球体,则黏结半径可由土壤含水率式
18、(3)推出。 3 3223 3 31 2 1 24 4 3 34 4 4 3 3 3BBR Rmm m R R R (3) 式中为土壤含水率,%;m1为土壤颗粒质量,kg;m2为水分质量,kg;1为土壤颗粒干密度,kg/m3;2为水分密度,kg/m3;R为土壤颗粒半径,m;RB为土壤颗粒间黏结半径,m。 1.4 EDEM仿真模型 在UG NX10中装配土槽和旋耕刀,导入到离散元软件EDEM2.6,仿真模型见图2。旋耕刀型号为IT245,材料为65Mn;仿真土槽尺寸依据实际的自制土槽试验台缩小为300 mm300 mm600 mm,材料为45钢;参照土壤颗粒度分析,土壤颗粒半径可取为区间平均值0
19、.625 mm,结合EDEM软件功能特点,适当放大颗粒半径,减小仿真时间,对仿真结果影响不大34,故将仿真土壤颗粒半径放大1倍,即R=1.25 mm,计算得仿真土壤颗粒黏结半径为RB=1.38 mm,共生成颗粒150 000个;材料之间的动摩擦系数、静摩擦系数及恢复系数依次设定为35:0.04、0.5、0.28(土壤颗粒-65Mn),0.25、0.4、0.2(土壤颗粒-土壤颗粒),0.04、0.5、0.28(土壤颗粒-45钢);由于生成土壤颗粒数目较多,设定07.5s为颗粒生成时间,7.58 s为黏结生成时间,88.5 s为仿真时间;通过土壤压缩仿真预试验,为降低孔隙率,使颗粒填充更加紧密,生
20、成颗粒时,在土槽长度方向施加幅度为1.5 mm,频率为2 Hz的振动。 1. 土壤颗粒 2. 土槽 3. 旋耕刀 4. 刀轴 1. Soil particles 2. Soil bin 3. Rotary blade 4. Cutter shaft 图2 EDEM仿真模型 Fig.2 Simulation model of EDEM 2 仿真模型验证 2.1 土槽试验 为检验EDEM仿真模型的正确性,在华南农业大学工程学院自制小型土槽试验台上进行了旋耕扭矩对比试验,试验台结构详见文献36。土槽长8 000 mm,宽650 mm,高600 mm;土壤分3层铺放,每一层均调湿、夯实,调整后的土壤平
21、均含水率为17%,平均紧实度为318 kPa;试验台车前进速度由M440变频器(7.5 kW,德国西门子公司)控制,旋耕刀转速由RF300A变频器(15 kW,上海红旗泰电子科技有限公司)调节,旋耕作业时刀轴所受扭矩由WDH型转矩转速传感器(300 Nm,北京沃德行世纪科技有限公司)测量;采用津旋245旋耕刀,单排安装,1把左弯刀,2把右弯刀,交错分布,见图3;调节耕深为100 mm,台车前进速度为0.5 m/s,分别测量转速为150、200、250、300和350 r/min时刀轴所受扭矩。 1. 土槽 2. 减速器 3. 传动轴 4. 左弯刀 5. 右弯刀1 6. 刀轴 7.右弯刀2 1.
22、Soil bin 2.Speed reducer 3.Transmission shaft 4.Left-hand blade 5.Right-hand blade 1 6.Cutter shaft 7.Right-hand blade 2 图3 旋耕刀安装图 Fig.3 Installation of rotary blades 2.2 扭矩对比分析 仿真模型设定旋耕刀最大入土深度为100 mm,刀轴直线移动速度为0.5 m/s,刀轴转速与土槽试验相同,数据保存间隔为0.001 s,做5组仿真试验,在数据分析模块获取试验结果,截取88.3 s内的数据,得刀轴所受扭矩变化规律如图4所示。 由图
23、4可知,5组试验中刀轴所受扭矩的变化趋势基本一致,均是先从0增加到最大值,后逐步减小到一个低值,接着又快速增加到一个高值,后回落,这与旋耕刀在切土过程中与土壤的接触状态相吻合。旋耕切土过程大致分4个阶段:第1阶段,刀片未切土前,扭矩几乎为0,在刀轴带动下,侧切刃纵向切入土层,切土体积逐步增加,扭矩变大,接着主切刃进入横向切土状态,主、侧切刃同时与土壤相互作用,扭矩出现骤然变大,随着刀片旋转,当旋耕刀到达最大耕深时,扭矩出现最大值;第2阶段,刀片越过最大耕深时,从侧切刃开始慢慢退出土层,切土体积变少,扭矩相应减小;第3阶段,扭矩降低到某一低值时,下一把刀片又开始进入切土状态,前、后2把刀同时切土
24、,扭矩又快速回升;第4阶段,当前一把刀片主切刃完全退出土层时,只剩下后一把刀片工作,扭矩又逐步回落。刀轴上对称分布3把刀片,此过程周而复始进行。 从图4a4e中扭矩最大值可知,刀具转速越高,所需的旋耕作业功耗就越大,与文献37-38的研究结论相一致。同时,随着转速变大,单位时间内的切土次数增加,出现峰值的时间间隔缩短,与实际情况相符。将仿真结果与土槽试验结果进行比,如图5所示,2条曲线相关系数为0.994,仿真值与试验值的变化趋势相同,扭矩均随转速增加而变大,最大相对误差为10%。仿真结果能较好的反映旋耕功耗变化规律,较准确的表达扭矩值,故用本文构建的EDEM仿真模型来研究旋耕刀所受三向工作阻
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 基于 离散 旋耕刀三 工作 阻力 仿真 分析 试验
限制150内