山地履带自走式辣椒收获机底盘动力学仿真与试验∗吉旭1,唐勇2,林蜀云2,张德俊1,张太华1,徐卫平1(1.贵州师范大学机械与电气工程学院,贵阳市,550025;2.贵州省山地农业机械研究所,贵阳市,550000)摘要:针对贵州山地土质松软,田块面积小㊁坡度大,现有辣椒收获机在山地上行驶困难等问题,设计一种适用于丘陵山地的履带自走式辣椒收获机,并以该收获机的底盘为对象,研究收获机在山地行驶过程中的通过性和稳定性㊂利用R e c u r D yn 对底盘在横坡行驶㊁纵坡行驶㊁翻越垂直壁和跨越壕沟等过程进行仿真㊂仿真结果表明,收获机在黏土路面上满载行驶时,横坡行驶最大坡度角为22ʎ,纵坡上坡最大坡度角为30ʎ,纵坡下坡最大坡度角为21ʎ,翻越垂直壁最大高度为510m m ,跨越壕沟最大宽度为1020m m ㊂田间试验结果表明,收获机纵坡上坡㊁翻越垂直壁和跨越壕沟的极限值与仿真结果的相对误差分别为10%㊁1.96%和3.92%,吻合度较高㊂试验验证了收获机在行驶过程中具有较好的稳定性和通过性,能够满足现阶段贵州山地辣椒采收要求㊂关键词:辣椒收获机;履带底盘;R e c u r D yn ;稳定性;通过性中图分类号:S 225.99㊀㊀文献标识码:A㊀㊀文章编号:20955553(2022)07007209吉旭,唐勇,林蜀云,张德俊,张太华,徐卫平.山地履带自走式辣椒收获机底盘动力学仿真与试验[J ].中国农机化学报,2022,43(7):72-80J iX u ,T a n g Y o n g ,L i nS h u y u n ,Z h a n g D e j u n ,Z h a n g T a i h u a ,X u W e i p i n g .D y n a m i cs i m u l a t i o na n dt e s t o f t h e t r a c k e d c h a s s i s o f t h em o u n t a i n o u s s e l f -p r o p e l l e d p e p p e rh a r v e s t e r [J ].J o u r n a l o fC h i n e s eA g r i c u l t u r a lM e c h a n i z a t i o n ,2022,43(7):72-80收稿日期:2021年7月31日㊀㊀修回日期:2021年9月14日∗基金项目:国家自然科学基金资助项目(72061006);黔科合平台人才([2017]5708);黔教合(Y J S C X J H 2020 105);黔科合支撑([2021]一般265)第一作者:吉旭,女,1996年生,白族,贵州毕节人,硕士研究生;研究方向为农业机械装备及自动化㊂E -m a i l :1909132458@q q .c o m 通讯作者:徐卫平,男,1967年生,贵州毕节人,教授,硕导;研究方向为农业机械装备及自动化㊂E -m a i l :2609550082@q q.c o m D yn a m i c s i m u l a t i o na n d t e s t o f t h e t r a c k e d c h a s s i s o f t h e m o u n t a i n o u s s e l f -p r o p e l l e d p e p pe r h a r v e s t e r J i X u 1,T a n g Y o n g 2,L i nS h u y u n 2,Z h a n g D e j u n 1,Z h a n g T a i h u a 1,X uW e i p i n g1(1.C o l l e g e o f M e c h a n i c a l a n dE l e c t r i c a lE n g i n e e r i n g ,G u i z h o uN o r m a lU n i v e r s i t y ,G u i y a n g ,550025,C h i n a ;2.G u i z h o uM o u n t a i nA g r i c u l t u r a lM a c h i n e r y R e s e a r c hI n s t i t u t e ,G u i y a n g ,550000,C h i n a )A b s t r a c t :A i m i n g a t t h e p r o b l e m so f s o f t s o i l i nG u i z h o u m o u n t a i n s ,s m a l l l a n da r e a ,l a r g es l o p e ,a n dt h e i n a b i l i t y o f t h ee x i s t i n gp e p p e r h a r v e s t e r t o t r a v e l o nt h em o u n t a i n ,ac r a w l e rs e l f -p r o p e l l e d p e p p e rh a r v e s t e rs u i t a b l e f o rh i l l y an d m o u n t a i n o u sa r e a sw a s d e s i g n e d .T h e c h a s s i s o f t h e h a r v e s t e r i s u s e d a s t h e o b j e c t ,a n d t h e p a s s a b i l i t y a n d s t a b i l i t y o f t h e h a r v e s t e r d u r i n g m o u n t a i nd r i v i n ga r e s t u d i e d .R e c u r D y n i su s e d t o s i m u l a t e t h e p r o c e s s o f c h a s s i s d r i v i n g o nh o r i z o n t a l s l o p e s ,v e r t i c a l s l o pe s ,o v e r v e r t i c a lw a l l s ,a n d c r o s s i n g t r e n c h e s .T h e s i m u l a t i o n r e s u l t s s h o wt h a tw h e n t h e h a r v e s t e r i s r u n n i n g o n a c l a y r o a dw i t h af u l l l o a d ,t h em a x i m u ms l o p e a ng l e o f c r o s s -s l o p e d r i v i n g i s 22ʎ,th em a xi m u ms l o p ea n g l eo f v e r t i c a l s l o p eu p h i l l i s30ʎ,t h em a x i m u ms l o p ea n g l eo f d o w n h i l l i s 21ʎ,t h em a x i m u mh e i g h t o v e r a v e r t i c a l w a l l i s 510m m ,a n d t h em a x i m u m w i d t h a c r o s s t h e t r e n c h i s 1020m m.T h e f i e l d t e s t r e s u l t s o f t h eh a r v e s t e r s h o wt h a t t h er e l a t i v ee r r o rb e t w e e nt h em a x i m u m v a l u eo f t h eh a r v e s t e ru p h i l l ,v e r t i c a lw a l l c r o s s i n g,a n dd i t c h c r o s s i n g a n d t h e s i m u l a t i o nr e s u l t s i s10%,1.96%,a n d3.92%.T h et e s t r e s u l t sa r es i m i l a r t ot h es i m u l a t i o nr e s u l t s .T h et e s t v e r i f i e s t h a t t h eh a r v e s t e rh a s g o o ds t a b i l i t y a n d p a s s a b i l i t y d u r i n g t h ed r i v i n g a n dc a n m e e t t h ec u r r e n tr e q u i r e m e n t so fG u i z h o u m o u n t a i n p e p p e r h a r v e s t i n g.K e yw o r d s :p e p p e r h a r v e s t e r ;t r a c k e d c h a s s i s ;R e c u r D y n ;s t a b i l i t y ;p a s s a b i l i t y Copyright ©博看网. All Rights Reserved.0㊀引言近年来,辣椒作为贵州省 十四五 规划中十二个特色优势扶贫产业之一,种植面积居全国第一[12]㊂辣椒采摘期短,刺激皮肤,人工收获效率低㊁成本高,机械化收获是实现贵州省辣椒全程机械化的关键[3]㊂目前制约贵州山地辣椒机械化收获的原因主要有:一是喀斯特地形阻碍了农业机械化发展;二是辣椒种植分散㊁杂乱,规模小;三是现有辣椒收获设备和技术不适用于贵州山地和土壤[4]㊂国内外广泛使用的辣椒收获机不能完全符合贵州山地地形和土壤类型㊂美国十方和P i k-R i t e两个公司生产的辣椒收获机功能齐全,收获的辣椒品质高㊁含杂少㊁破损率低,但体型庞大,不适用于山地丘陵作业[56]㊂国内辣椒收获机发展相对滞后,以4J Z-3600A自走式辣椒收获机和雷神4J Z-2300型辣椒收获机为代表,机器结构紧凑㊁配置合理,可一次性完成辣椒的采摘㊁输送和收集[79]㊂陈长林等[10]设计的刷辊式辣椒收获机㊁黄晓鹏等[11]设计的滚筒弹齿式辣椒收获机以及张祥军等[12]设计的梳指式辣椒收获机均采用轮式底盘㊂轮式辣椒收获机虽在新疆㊁河北和河南等平地取得良好的推广应用,但轮式底盘与地面为线接触,接触面积较小,在丘陵山地这种非结构化的黏土地面上行走时存在稳定性差㊁容易发生倾翻等问题[13]㊂履带底盘具有接地面积大㊁接地比压小㊁附着性能好㊁爬坡能力强㊁转弯半径小等特点,在农业机械上得到广泛应用,而基于履带底盘的动力学分析也受到越来越多的关注㊂欧阳益斌[14]针对丘陵山地果林深松作业效率低㊁通过性差等问题,分析了深松机具及其履带底盘的受力情况,并基于R e c u r D y n软件进行了动力学仿真试验㊂王峰[15]利用R e c u r D y n对履带底盘进行动力学仿真分析,进一步论证了履带底盘在丘陵山地果园中的良好应用效果㊂基于贵州山地地形的复杂性及土质的特殊性,本文以山地履带自走式辣椒收获机为研究对象,利用R e c u r D y n对收获机底盘在横坡行驶㊁纵坡行驶㊁翻越垂直壁和跨越壕沟过程进行仿真,并通过田间试验来验证仿真分析结果,以期为山地辣椒实现机械化收获提供理论和实践基础㊂1㊀整机结构根据贵州辣椒种植农艺㊁地形和土壤等实际情况,确定收获机行走装置为具有更高稳定性和通过性的倒梯形履带底盘㊂整机由采摘装置㊁输送装置㊁清选装置㊁收集装置和行走装置组成,可一次性完成辣椒的采摘㊁输送㊁清选和收集,收获机尺寸为5570m mˑ1852m mˑ3372m m,结构如图1所示㊂收获机作业时,发动机动力传送到行走系统,驱动底盘行驶,液压缸调整采摘装置位置,液压马达驱动采摘滚筒转动,滚筒上的弹齿拉动辣椒并克服辣椒与茎干间的连接力,使辣椒从茎干上脱离[16]㊂采摘后的辣椒通过滚筒转动被抛送至输送带,经过清选后输送到料箱,完成采摘㊂当辣椒收装至集料箱体积三分之二时,液压缸举升集料箱进行卸料,完成收获全过程㊂图1㊀收获机结构示意图F i g.1㊀S c h e m a t i c d i a g r a mo f t h e s t r u c t u r e o f t h eh a r v e s t e r1.机架㊀2.集料箱㊀3.二级输送带㊀4.清选风机㊀5.驾驶室6.一级输送带㊀7.收获装置㊀8.发动机㊀9.履带底盘收获机履带底盘主要由支架㊁驱动轮㊁导向轮㊁拖带轮㊁支重轮㊁游离三角和张紧装置等部件组成,整体尺寸为1800m mˑ350m mˑ580m m,结构如图2所示,结构参数如表1所示㊂图2㊀履带底盘结构示意图F i g.2㊀S c h e m a t i c d i a g r a mo fc r a w l e r c h a s s i s s t r u c t u r e1.车架㊀2.驱动轮㊀3.车体部分㊀4.拖带轮㊀5.导向轮6.张紧装置㊀7.支重轮㊀8.悬架㊀9.游离三角㊀10.履带表1㊀履带底盘结构参数㊀㊀辣椒收获机底盘运动时,驱动轮将液压马达的能量转化为底盘行驶的动能,拖带轮起到支撑履带的作用,导向轮和张紧装置保证了履带的卷绕方向并承受行驶时履带受到的冲击力㊂支重轮通过游离三角与车架铰接,其位置可根据地形的变化进行调节,避免了支Copyright©博看网. All Rights Reserved.重轮悬空时对底盘产生的载荷分布不均,从而提高了底盘的稳定性[17]㊂2㊀底盘稳定性及通过性分析山地履带自走式辣椒收获机底盘需要具有良好的通过性和行驶平顺性[1819]㊂为简化分析,将机体受到的作用力分解在沿横向和纵向两个方向,即从横坡和纵坡两个方向分别讨论其不同工况下的行驶稳定性[20],从翻越垂直壁和跨越壕沟两个方面分别讨论其通过性㊂2.1㊀稳定性分析收获机在横坡上行驶时受自身的重力㊁坡面的支撑力和地面的摩擦力,受力分析如图3所示㊂图3㊀辣椒收获机横坡行驶受力简图F i g .3㊀S i m p l i f i e dd i a gr a mo f t h e f o r c e o n t h e c r o s s s l o p e o f t h e p e p pe r h a r v e s t e r ㊀注:O 为质心;G 为重力;α为横向坡度;S 为轨距;e 为质心偏移量;h 为质心高度;N 1,N 2为地面对左右支重轮沿横向方向的支撑力;f 1,f 2为地面对履带的摩擦力㊂在横坡上匀速行驶,当坡度超过极限时,底盘侧向附着力会小于重力的下滑分量,收获机将发生侧滑或侧翻[21]㊂由经验公式得最大滑移坡度临界值a n 需要满足关系式φG c o s αn =G s i n αn (1)式中:φ附着系数;a n最大坡度㊂当收获机在横坡上稳定行驶时,车体受力平衡,对A 点取力矩可得N 2S -G c o s α(0.5S +e )+Gh s i n α=0(2)式(2)化简后得N 2=G c o s α(0.5S +e )-G h s i n αS(3)可知车体发生侧翻的临界条件是N 2=0,即稳定行驶条件为G c o s α(0.5S +e )-G h s i n α>0,设αm a x 为车体在横坡上稳定行驶的极限坡度,且e =0,则分析可得αm a x =ar c t a n 0.5S +e h =a r c t a n S 2h=23ʎ(4)收获机纵向上坡和下坡时受自身的重力㊁坡面的支撑力和地面的摩擦力,受力分析如图4所示㊂(a)上坡(b)下坡图4㊀辣椒收获机纵坡上坡㊁下坡行驶受力图F i g .4㊀S i m p l i f i e dd i a g r a mo f t h e v e r t i c a l s l o pe d r i v i n gf o r c e o f t h e p e p pe r h a r v e s t e r ㊀注:β为纵向坡度;m 为支重轮C 与支持力N 3间的距离;L 为前后支重轮间距;a ,b 为支重轮与质心位置的横向距离;N 3,N 4为地面对履带底盘的支持力;F 为地面对履带的摩擦力㊂收获机在纵坡上稳定上坡时,整车受力平衡,对支重轮C 列力矩方程可得N 3m -G a c o s β+G h s i n β=0(5)由于整车受力平衡,因此在垂直于斜坡方向上的合力为零,即N 3=G c o s β,代入式(5)得m =a c o s β-h s i n βc o s β(6)当坡度β超过一定的限度时,车体将发生倾翻,因此车体发生倾翻的临界条件是m =0,即a c o s β-h s i n β=0(7)则纵向上坡最大坡度角βm a x 上=a r c t a n (a /h )=28ʎ(8)同理,纵向下坡最大坡度角βm a x 下=a r c t a n (b /h )=21.5ʎ(9)2.2㊀翻越垂直壁通过性分析底盘在翻越垂直障碍物时,翻越过程包括底盘前端越障㊁支持段履带越障和机体质心越过障碍物[22]㊂底盘前端越障是指底盘导向轮与障碍物边缘接触时,履带在牵引力的作用下转动,底盘持续行驶到底盘前支重轮爬上障碍物的过程㊂支持段履带越障阶段是指前进过程中底盘与地面的夹角逐渐增大,机体质心不断上升且缓慢前移,质心与障碍物的垂直面重合时的过程㊂机体质心越过障碍物阶段是指底盘落到障碍物的水平表面到越障完成过程㊂质心与障碍物垂直面重Copyright ©博看网. All Rights Reserved.合时的临界状态如图5所示㊂图5㊀底盘翻越垂直壁临界状态F i g .5㊀As i m p l i f i e dd i a g r a mo f t h e f o r c e a p pl i e d t o t h e p e p p e r h a r v e s t e r i n t h e p r o c e s s o f c r o s s i n g ob s t ac l e s ㊀注:θ1为底盘极限偏转角度;H 为底盘极限越障高度;r 1为驱动轮半径,r 2为支重轮半径,h 1为导向轮高度;L 1为驱动轮到导向轮的距离;N 5为障碍物对底盘的支撑力㊂底盘翻越最大高度的临界状态时,整车受力平衡,对支重轮E 列力矩方程N 412L -h t a n θæèöø+(h -r 1)G s i n θ-12L G c o s θ=0(10)底盘达最大临界越障高度的临界条件为N 4=0,则(h -r 1)G s i n θ1-12L G c o s θ1=0(11)其中,履带极限偏转角度θ1=ar c t a n L 2(h 1-r 2)(12)越障高度表达式h (θ)=L 2+r 2t a n θ2-h 1t a n θæèöø㊃s i n θ(13)θ=30ʎ时,实际越障高度H 实=m i n (h 1,h (θ)m a x )=530m m (14)0<θ<θ1,最大高度不得超过底盘前轮中心离地高度㊂2.3㊀跨越壕沟通过性分析底盘跨越壕沟宽度与其质心位置㊁支持段履带长度㊁支重轮大小和行驶速度有关[23]㊂机体水平跨越壕沟,当其质心与壕沟近侧边界线垂直面重合时,履带发生偏转,底盘前支重轮与壕沟的远侧边界线接触,此时壕沟宽度为机体跨越壕沟临界宽度,如图6所示㊂图6㊀底盘跨越壕沟临界状态F i g .6㊀As c h e m a t i c d i a g r a mo f t h e f o r c e a p pl i e d t o t h e p e p p e r h a r v e s t e r d u r i n g t h e p r o c e s s o f c r o s s i n g th e d i t c h 机体质心越过壕沟远侧边界线后,底盘落在壕沟远侧的水平路面上,跨沟完成㊂收获机履带底盘跨越壕沟最大宽度M =L 12æèöø2+h 12+r 1=1059m m (15)3㊀底盘动力学模型的建立底盘动力学建模时,由于本身的零件和自由度数量过多,使用传统的动力学方程无法对其进行推导㊂目前常见的动力学建模研究方法分别有牛顿 欧拉(N e w t o n -E u l e r )法㊁拉格朗日(L a g r a n ge )法㊁凯恩(K a n e )法以及罗伯逊 维登伯格(R o b e r s o n -W i t t e n b e r g)法等,其中前两种方法最为常用㊂利用牛顿 欧拉法建立动力学方程时推导过程较为简单,但需要隔离机体中的每一个刚体,且每个刚体相互间都存在着约束力,建立的方程较多㊂利用拉格朗日法建立动力学方程时,基于系统能量的概念,以简单的形式表达系统的动力学状态,方程较少,但推导过程较为复杂㊂综合,选用拉格朗日法来建立底盘动力学方程,并利用R e c u r D y n 中的T r a c k (L M )模块对整机进行参数化建模并完成装配㊂3.1㊀约束力的设置将收获机三维模型从S o l i d W o r k s 中简化后导入R e c u r D yn ,得到收获机动力学仿真模型㊂定义各部件之间的约束,底盘和车体各部件间的约束关系如表2所示㊂表2㊀各部件间约束关系T a b .2㊀C o n s t r a i n t r e l a t i o n s h i p b e t w e e n c o m p o n e n t s 基体活动体接触类型约束个数车体驱动轮R e v o l u t e J o i n t 2车体履带架F i x e d J o i n t8张紧装置导向轮R e v o l u t e J o i n t 2履带架拖带轮R e v o l u t e J o i n t 2履带架支重轮1R e v o l u t e J o i n t 2履带架支重轮2R e v o l u t e J o i n t 2履带架支重轮3R e v o l u t e J o i n t 2履带架支重轮4R e v o l u t e J o i n t 2履带架支重轮5R e v o l u t e J o i n t 2履带架支重轮6R e v o l u t e J o i n t 2履带架支重轮7R e v o l u t e J o i n t 2履带架支重轮8R e v o l u t e J o i n t 23.2㊀驱动函数的设置在R e c u r D yn 中建立三维模型和设置路面环境后,根据收获机实际行驶的不同档位速度对履带两侧的驱动轮添加驱动函数[24]㊂Ⅰ挡速度为1.8k m /h ,设置的驱动函数为S T E P (T I M E ,0,0D ,0,-400D );Copyright ©博看网. All Rights Reserved.Ⅱ挡速度为3.9k m /h ,设置的驱动函数为S T E P(T I M E ,0,0D ,0,-867D );Ⅲ挡的速度为5.3k m /h ,设置的驱动函数为S T E P (T I M E ,0,0D ,0,-1178D )㊂3.3㊀路面模型的建立根据经典土壤力学理论[25]得到黏土路面物理性能参数如表3㊂土壤的承压和剪切力学特性参数直接影响履带底盘的通过性和牵引性,为了给收获机提供更准确的地面力学模型,采用M.G.B e k k e r 压力 沉陷关系模型进行模拟[2627],即P =k cB +k φæèöøz n (16)式中:P 接地比压,k P a;Z 土壤沉陷量,m ;n土壤变形指数;B 承载面的短边,m ;k c粘聚变形模量,k N /(m n +1);k φ摩擦变形模量,k N /(m n +2)㊂表3㊀黏土路面物理性能参数T a b .3㊀P h y s i c a l p e r f o r m a n c e p a r a m e t e r s o f c l a ypa v e m e n t 参数数值土壤变形指数n 0.5剪切变形模量/m 5剪切内阻力角φ/(ʎ)20粘聚变形模量k c /(k N ㊃m -n -1)4.127摩擦变形模量k φ/(k N ㊃m -n -2)0.433内聚力P /P a0.0694㊀动力学仿真分析底盘的稳定性体现在横坡和纵坡的行驶状况,通过性体现在跨越壕沟和翻越垂直障碍物的能力[2829]㊂通过理论分析计算和整机满载动力学仿真试验,得出不同路况下的理论计算结果和仿真结果如表4所示㊂分析仿真中收获机侧向偏移量㊁俯仰角和质心垂向速度的变化,得出机体在不同路况下的行驶状况㊂表4㊀不同路况理论及仿真结果T a b .4㊀T h e o r y an d s i m u l a t i o n r e s u l t s o f d i f f e r e n t r o a d c o n d i t i o n s指标理论值仿真值横坡行驶最大坡度角αm a x/(ʎ)2222纵坡上坡最大坡度角βm a x 上/(ʎ)2830纵坡下坡最大坡度角βm a x 下/(ʎ)21.521翻越垂直壁最大高度H /m m 530510跨越壕沟最大宽度M /m m105910204.1㊀横坡行驶收获机在长度为100m ,坡度为15ʎ㊁20ʎ和25ʎ的黏土横坡上仿真时的机体侧向偏移量如图7所示㊂图7㊀不同坡度下车体侧向偏移量F i g .7㊀V e r t i c a l d i s pl a c e m e n t o f t h e c e n t e r o f m a s su n d e r d i f f e r e n t s l o pe s 根据G B /T15370.4 2012,收获机偏移量不得超过6%㊂分析图8可知,收获机在15ʎ和20ʎ横坡上行驶时,滑移量分别为1.27%和4.84%;收获机在25ʎ横坡上行驶时,行驶一段时间后出现急剧滑移并发生侧向倾翻,滑移量高达7.38%㊂试验结果表明,收获机在坡度为25ʎ的横坡上行驶时存在横向附着力较小,容易发生倾翻等问题,出于安全性能考虑,收获机应在小于20ʎ的横坡上行驶㊂4.2㊀纵坡行驶在20ʎ的黏土纵坡上,对机体不同档位下的纵向上坡过程进行仿真,仿真时长为25s ,步长为1000步,得到不同档位下俯仰角和质心垂向速度的变化如图8所示㊂分析图8(a )可知,机体以Ⅰ挡上坡时,0~6s 内在平地上加速启动,俯仰角在0ʎ上下波动;6~18s 机体处于上坡阶段,先从平地过渡到斜坡,俯仰角随之上升到22ʎ,爬上斜坡后,俯仰角在20ʎ左右;18~25s 机体上坡完成并在高处水平面上行驶,俯仰角回到0ʎ附近㊂由图8(b )可知,机体在行驶过程中发生抖动导致质心垂向速度出现了正负值变化,当机体以Ⅰ挡上坡时,速度最大值为749.4m m /s ;以Ⅱ挡上坡时,速度最大值为1381.4m m /s;以Ⅲ挡上坡时,机体发生倾翻㊂由此可见,上坡时速度越大则时间越短,但速度越大质心垂向速度波动也越大,底盘稳定性越差,机体极易发生倾翻㊂因此,机体纵向上坡时应以低速挡行驶㊂在15ʎ的黏土纵坡上,对机体不同档位的下坡过程进行仿真,仿真时长为20s ,步长为1000步,得到不同档位下俯仰角和质心垂向速度的变化如图9所示㊂Copyright ©博看网. All Rights Reserved.(a)俯仰角变化(b)质心垂向速度变化图8㊀收获机爬坡20ʎ时动力学仿真F i g .8㊀D y n a m i c s i m u l a t i o no f h a r v e s te r c l i m b i n g 20ʎ(a)俯仰角变化(b)质心垂向速度变化图9㊀收获机下坡15ʎ时动力学仿真F i g .9㊀D yn a m i c s i m u l a t i o nw h e n t h e h a r v e s t e r g o e s d o w n h i l l a t 15ʎ分析图9可知,机体以Ⅰ挡下坡时,0~6s 内在平地上加速启动,此时俯仰角在0ʎ上下波动;6~10s 从平地过渡到斜坡,俯仰角最小值达-28ʎ;10~16s 在斜坡上匀速下坡,俯仰角在-16ʎ左右;16~20s 完成下坡动作,随后过渡到低处水平面上继续行驶,俯仰角为0ʎ,质心垂向速度最大值为260m m /s ,最小值为-1536m m /s ㊂机体以Ⅱ挡下坡时,俯仰角波动较大,最大值达到20.2ʎ,最小值为33.3ʎ,质心垂向速度最大值为1237.7m m /s ,最小值为-2696.4m m /s ㊂机体以Ⅲ挡下坡时发生倾翻㊂由仿真得,机体纵向下坡时速度越大时间越短,质心垂向速度波动也越大,底盘稳定性越差,机体极易发生倾翻㊂因此,机体纵向下坡时应以低速挡行驶㊂4.3㊀翻越垂直壁在450m m 高的垂直壁黏土路面上,机体以不同档位进行翻越垂直壁仿真,仿真时长为10s,步长为1000步,得到不同档位下俯仰角和质心垂向速度的变化如图10所示㊂(a)俯仰角变化(b)质心垂向速度变化图10㊀收获机翻越垂直壁450m m 时动力学仿真F i g .10㊀D yn a m i c s i m u l a t i o nw h e n t h eh a r v e s t e r c r o s s e s t h e v e r t i c a lw a l l 450m m收获机在翻越垂直壁时,要保证翻过垂直壁后机体还能恢复到在平地上行驶㊂分析图10(a )可知,机体以Ⅰ挡翻越垂直壁时最大俯仰角为14.8ʎ,以Ⅱ挡翻越时的最大俯仰角16.6ʎ,可见机体在Ⅰ㊁Ⅱ挡位下能Copyright ©博看网. All Rights Reserved.够顺利通过450m m的垂直壁障碍并且恢复平地行驶状态㊂机体在Ⅲ挡下越障时由于俯仰角超过了底盘的纵向倾翻角,机体发生倾翻㊂分析图10(b)可知,机体以Ⅰ挡越障时质心垂向速度最大值为947.9m m/s,最小值为-305.4m m/s;以Ⅱ挡越障时质心垂向速度最大值为815.6m m/s,最小值为-481m m/s,可见机体以Ⅰ挡越障时更为平缓,低速挡越障更为安全㊂4.4㊀跨越壕沟机体以不同档位速度通过700m m宽的壕沟时,仿真结果如图11所示,仿真时长为10s,步长为1000步㊂(a)俯仰角变化(b)质心垂向速度变化图11㊀收获机跨越壕沟700m m时动力学仿真F i g.11㊀D y n a m i c s i m u l a t i o nw h e n t h eh a r v e s t e rc r o s s e s t h e t r e n c h700m m收获机在跨越壕沟时,要保证跨过壕沟后还能恢复到在平地上行驶㊂分析图11(a)可知,机体以Ⅰ挡跨越壕沟时最大俯仰角为-9.37ʎ,以Ⅱ挡跨越壕沟时最大俯仰角为-7.99ʎ,机体以Ⅰ㊁Ⅱ挡位能够顺利跨越700m m宽的壕沟㊂当以Ⅲ挡跨越壕沟时,由于速度过大,机体前支重轮在垂直方向上的位移远大于水平方向,在惯性力的作用下,机体前部落入壕沟,跨沟失败㊂分析图11(b)可知,机体以Ⅰ挡跨越壕沟时质心垂向速度最大值为532.8m m/s,最小值为-421.5m m/s;以Ⅱ挡跨越壕沟时质心垂向速度最大值为591.2m m/s,最小值为-662m m/s㊂机体在以Ⅱ挡跨越壕沟时的质心垂向速度波动比Ⅰ挡时大,表明在机体低速挡下跨越壕沟更为平缓,应控制跨越壕沟速度在Ⅰ到Ⅱ挡之间㊂5㊀田间试验试验目的:对贵州山地履带自走式辣椒收获机在纵向上坡㊁翻越垂直壁和跨越壕沟等方面的性能进行试验,验证其结构设计和技术参数的合理性,分析收获机在行驶稳定性和通过性方面存在的问题和不足,并比较与理论计算和仿真试验间的差异㊂试验条件:2021年4月20日于遵义市播州区辣椒种植基地对收获机进行纵向上坡㊁翻越垂直壁和跨越壕沟试验㊂试验基地土质为黏土,测得土壤含水率约24%,土壤硬度值约48N/c m3㊂5.1㊀纵向爬坡试验本试验对收获机的纵向上坡通过性进行测试,在试验场地选择不同坡度的纵坡,收获机以不同档位通过斜坡,记录上坡过程消耗的时间,上坡试验结果如图12所示㊂图12㊀上坡时间与坡度角关系F i g.12㊀R e l a t i o n s h i p b e t w e e nu p h i l l t i m e a n d s l o p e a n g l e分析图12可知,收获机以Ⅰ挡进行上坡时,在20ʎ以内随着坡度的缓慢增大上坡时间无明显变化,坡度在20ʎ之后上坡时间明显增大,上坡度增加到27ʎ时达到最大上坡角度㊂以Ⅱ挡和Ⅲ挡进行上坡时,机体本身牵引力较小,在上坡度为20ʎ以后就不能再上坡成功㊂试验得出的极限上坡度与前文理论计算得到的极限上坡度28ʎ相比误差为3.57%,与前文仿真得出的极限上坡度30ʎ相比误差为10%,误差产生的原因是理论计算和仿真时得到的极限上坡度是在理想环境下得出的,而实际的田间土壤湿度大㊁黏度重,在进行上坡时路面湿滑导致履带的地面附着力较小㊂5.2㊀翻越垂直壁试验本试验对收获机翻越不同高度垂直壁进行试验,以收获机底盘前端接触垂直壁到底盘后端离开垂直壁的时间作为评价指标,试验结果如表5所示㊂翻越垂直壁田间试验中,当以Ⅲ挡越障时速度过快会导致对Copyright©博看网. All Rights Reserved.底盘的冲击力过大,在越障过程中容易发生倾翻㊂因此选择Ⅰ挡和Ⅱ挡进行试验㊂表5㊀不同档位翻越垂直壁结果T a b.5㊀R e s u l t s o f c r o s s i n g t h e v e r t i c a lw a l lw i t hd i f f e r e n t g e a r s档位垂直壁高度/m m 300350400450470500510Ⅰ档速度(1.8k m/h)通过通过通过通过通过缓慢通过无法通过通过时间/s5.405.636.557.958.9010.220Ⅱ档速度(3.9k m/h)通过通过通过缓慢通过无法通过通过时间/s5.556.918.5310.040㊀㊀分析表5可知,当通过时间超过10s为缓慢翻越垂直壁㊂收获机以Ⅰ挡翻越高度为500m m垂直壁时机体越障困难,当以Ⅱ挡翻越垂直壁时,垂直壁高度为450m m时机体越障困难㊂翻越垂直壁田间试验实测值500m m与理论值530m m相比误差为5.66%,与仿真值510m m相比误差为1.96%㊂误差产生的原因是在田间试验中履带的地面附着力较小,当达到一定越障高度后,越障过程中机体的前端履带会被垂直壁的棱顶回,此时履带出现打滑现象,导致越障失败㊂5.3㊀跨越壕沟试验对收获机进行跨越壕沟田间试验时,以收获机底盘前端接触壕沟近侧到底盘后端离开壕沟远侧的时间作为评价指标,试验结果如表6所示㊂使用Ⅲ挡跨越壕沟时速度过快会使机体受到猛烈撞击,悬挂压力猛增,造成悬挂变形,而且Ⅲ挡的牵引力有限,会导致机体有陷入壕沟的危险㊂因此选择Ⅰ挡和Ⅱ挡进行试验㊂表6㊀不同档位跨越壕沟结果T a b.6㊀R e s u l t s o f c r o s s i n g t r e n c h e s i nd i f f e r e n t g e a r s档位壕沟宽度/m m 5006007008009009509801000Ⅰ档速度(1.8k m/h)通过通过通过通过通过通过缓慢通过无法通过通过时间/s5.916.806.937.819.349.8910.220Ⅱ档速度(3.9k m/h)通过通过通过通过通过缓慢通过无法通过通过时间/s6.006.707.698.679.6110.110㊀㊀分析表6可知,当通过时间超过10s为缓慢跨越壕沟㊂当不断加大壕沟宽度至接近理论值时,收获机以Ⅰ挡跨越宽度为980m m的壕沟时机体跨沟困难,收获机以Ⅱ挡跨越宽度为950m m的壕沟时机体跨沟困难㊂跨越壕沟田间试验实测值980m m与理论值1059m m相比,误差为7.46%,与仿真值1020m m 相比,误差为3.92%㊂误差原因是收获机实际跨越壕沟过程中后支重轮掉入壕沟内时前轮摩擦推力不足导致跨沟失败㊂6㊀结论本文利用R e c u r D y n对山地履带自走式辣椒收获机的横坡行驶㊁纵坡上坡㊁纵坡下坡㊁翻越垂直壁和跨越壕沟等进行动力学仿真,并对收获机进行田间试验,得出以下结论㊂1)对辣椒收获机行驶稳定性及通过性进行理论分析,得到横坡行驶最大坡度角为22ʎ,纵坡上坡最大坡度角为28ʎ,下坡最大坡度角为21.5ʎ,翻越垂直壁最大高度为530m m,跨越壕沟最大宽度为1059m m㊂2)对辣椒收获机进行动力学仿真,得到横坡行驶最大坡度角为22ʎ,纵坡上坡最大坡度角为30ʎ,下坡最大坡度角为21ʎ,翻越垂直壁最大高度为510m m,跨越壕沟最大宽度为1020m m㊂辣椒收获机驱动速度从Ⅰ挡(1.8k m/h)上升到Ⅱ挡(3.9k m/h)时,在纵向上坡㊁下坡㊁翻越垂直壁和跨越壕沟时的质心垂向加速度的最大值分别增大了632m m/s㊁977.7m m/s㊁132.3m m/s和58.4m m/s㊂3)对辣椒收获机进行田间试验,得到收获机在实际作业过程中纵向上坡极限坡度角为27ʎ,翻越垂直壁最大高度为500m m,跨越壕沟最大宽度为980m m㊂试验结果与理论计算及仿真分析结果相吻合,能够满足贵州丘陵山地的辣椒收获的要求㊂参㊀考㊀文㊀献[1]王洪亮,吴康云.对贵州辣椒产业技术体系建设的思考[J].辣椒杂志,2020,18(2):14-17.W a n g H o n g l i a n g W uK a n g y u n.T h o u g h t so n t h e c o n s t r u c t i o n o f g u i z h o u p e p p e r i n d u s t r i a l t e c h n o l o g y s y s t e m[J].J o u r n a l o f C h i n aC a p s i c u m,2020,18(2):14-17.Copyright©博看网. All Rights Reserved.[2]贵州:从辣椒大省迈向辣椒强省[J].中国蔬菜,2020 (9):105.[3]金学欢.贵州省辣椒产业发展研究[D].贵阳:贵州大学,2018.[4]周丕东,刘春波,张佩.贵州山地农业机械化发展现状及对策建议[J].中国农机化学报,2020,41(7):231-236. Z h o uP i d o n g L i u C h u n b o Z h a n g P e i.P r e s e n ts i t u a t i o na n d c o u n t e r m e a s u r e so fa g r i c u l t u r a lm e c h a n i z a t i o ni n m o u n t a i n o u s a r e a s o fG u i z h o uP r o v i n c e[J].J o u r n a l o fC h i n e s eA g r i c u l t u r a l M e c h a n i z a t i o n,2020,41(7):231-236.[5]H e r c kL V,K u r t s e rP,W i t t e m a n sL,e t a l.C r o p d e s i g nf o ri m p r o v e dr o b o t i ch a r v e s t i n g:A c a s es t u d y o fs w e e t p e p p e rh a r v e s t i n g[J].B i o s y s t e m s E ng i n e e r i n g,2020, 192:294-308.[6]K i m T H,K i m DC,C h oY.P e r f o r m a n c e c o m p a r i s o na n de v a l u a t i o no ft w os m a l lc h i l i p e p p e rh a r v e s t e r p r o t o t y p e s t h a t a t t a c ht o w a l k i n g c u l t i v a t o r s[J].A p p l i e dS c i e n c e s, 2020,10(7):2570.[7]胡爽吉.梳齿式辣椒采摘装置的设计研究[D].石河子:石河子大学,2012.[8]张俊三,阿力木㊃买买提吐尔逊,李谦绪,等.辣椒收获机(4J Z-3600A自走式)[P].中国专利:C N306106326S, 2020-10-16.[9]赵飞龙.辣椒收获机液压系统设计研究[D].石河子:石河子大学,2020.[10]陈长林,石磊,孙勇飞,等.一种刷辊式辣椒收获机[P].中国专利:C N209994889U,2020-01-31. [11]黄晓鹏,万芳新,戴立勋,等.一种线辣椒收获机[P].中国专利:C N211129013U,2020-07-31.[12]张祥军,张胜楠.一种辣椒收获机[P].中国专利:C N109287269A,2019-02-01.[13]解福祥,徐萌萌,张训林,等.果园运输机器人的结构设计与试验[J].中国农机化学报,2020,41(5):170-175.X i e F u x i a n g X u M e n g m e n g Z h a n g X u n l i n,e t a l. S t r u c t u r a l d e s i g n a n d s i m u l a t i o n a n a l y s i s o f s h r i m p o r c h a r d t r a n s p o r t r o b o t[J].J o u r n a l o f C h i n e s eA g r i c u l t u r a lM e c h a n i z a t i o n,2020,41(5):170-175.[14]欧阳益斌.果园深松机底盘典型工况动力学分析与试验[D].长沙:中南林业科技大学,2018. [15]王锋.丘陵山地果园动力底盘的坡地通过性研究[D].重庆:西南大学,2020.[16]宋韩韩.基于A D AM S的刚柔耦合整车模型平顺性仿真研究[D].锦州:辽宁工业大学,2015.[17]E d w i nP,S h a n k a rK,K a n n a nK.S o f t s o i l t r a c k i n t e r a c t i o nm o d e l i n g i ns i n g l er i g i d b o d y t r a c k e dv e h i c l e m o d e l s[J]. J o u r n a l o fT e r r a m e c h a n i c s,2018,77:1-14. [18]鲍向东,刘锡勇.贵州丘陵山区小型多功能田间作业机研究[J].农业机械,2011(28):120-121.[19]刘平义,彭凤娟,李海涛,等.丘陵山区农用自适应调平底盘设计与试验[J].农业机械学报,2017,48(12): 42-47.L i u P i n g y i,P e n g F e n g j u a n,L i H a i t a o.D e s i g n a n de x p e r i m e n t of a d a p t i v e l e v e l i ng ch a s si s f o rh i l l y a r e a[J].T r a n s a c t i o n s o f t h e C h i n e s e S o c i e t y f o r A g r i c u l t u r a l M a c h i n e r y,2017,48(12):42-47.[20]孙术发,任春龙,李涛,等.基于履带式底盘的改进型森林消防车通过性[J].农业工程学报,2018,34(17): 61-67.S u nS h u f a,R e nC h u n l o n g,L iT a o,e ta l.T r a f f i c a b i l i t y a n a l y s i so f i m p r o v e df o r e s tf i r ee n g i n eb a s e do nc r a w l e r c h a s s i s[J].T r a n s a c t i o n s o f t h e C h i n e s e S o c i e t y o fA g r i c u l t u r a l E n g i n e e r i n g,2018,34(17):61-67.[21]张拓.小型山地履带底盘行驶性能仿真分析及试验研究[D].重庆:重庆理工大学,2019.[22]黄晗,李建桥,吴宝广,等.轻载荷条件下轻型车辆车轮牵引通过性模型的建立与验证[J].农业工程学报,2015, 31(12):64-70.H u a n g H a n,L i J i a n q i a o,W u B a o g u a n g,e t a l.C o n s t r u c t i o na n dv e r i f i c a t i o no f l i g h t w e i g h tv e h i c l ew h e e l t r a c t i v e t r a f f i c a b i l i t y m o d e l u n d e r l i g h t l o a d[J].T r a n s a c t i o n s o f t h e C h i n e s e S o c i e t y o f A g r i c u l t u r a lE n g i n e e r i n g,2015,31(12):64-70.[23]吴哲,马岩,杨春梅,等.简易履带式行走机构的设计及通过性分析[J].林业机械与木工设备,2013,41(12): 14-16.W uZ h e,M aY a n,Y a n g C h u n m e i,e t a l.D e s i g na n d p a s s a g e a n a l y s i s o f as i m p l e c r a w l e rw a l k i n g m e c h a n i s m[J].F o r e s t r yM a c h i n e r y&W o o d w o r k i n g E q u i p m e n t,2013,41(12): 14-16.[24]李浩.一种轮履复合救援机器人底盘的动力学分析与优化[D].天津:天津理工大学,2015.[25]王颖,李建桥,张广权,等.仿生步行足沙地力学特性研究[J].农业机械学报,2016,47(2):384-389.W a n g Y i n g,L i J i a n q i a o,Z h a n g G u a n g q u a n,e t a l.M e c h a n i c a l c h a r a c t e r i s t i c so fb i o n i c w a l k i n g f o o to ns o f t s a n d[J].T r a n s a c t i o n s o f t h e C h i n e s e S o c i e t y f o rA g r i c u l t u r a lM a c h i n e r y,2016,47(2):384-389.[26]芮强,王红岩,王钦龙,等.履带车辆转向性能参数分析与试验研究[J].机械工程学报,2015,51(12):127-136.R u i Q i a n g,W a n g H o n g y a n,W a n g Q i n l o n g,e t a l.R e s e a r c h o n t h e a c q u i s i t i o n o f s t e e r i n g p e r f o r m a n c e p a r a m e t e r s o f a r m o r e dv e h i c l eb a s e do ne x p e r i m e n t s[J].T r a n s a c t i o n s o f t h e C h i n e s e S o c i e t y o f A g r i c u l t u r a lE n g i n e e r i n g,2015,51(12):127-136.[27]谢铌.小型山地履带底盘坡地行驶性能分析与试验研究[D].重庆:重庆理工大学,2020.[28]崔雪斌,张宏,石涛.基于链环不均匀系数的履带车辆行驶平顺性分析[J].工程设计学报,2018,25(1):71-78.C u i X u e b i n,Z h a n g H o n g,S h i T a o.R i d e c o m f o r t a n a l y s i s o f t r a c k e d v e h i c l eb a s e do nn o n u n i f o r mc o e f f i c i e n to f l i n k[J].C h i n e s e J o u r n a l o f E n g i n e e r i n gD e s i g n,2018,25(1): 71-78.[29]姚禹.基于履带 地面耦合系统的低速履带车辆通过性研究[D].长春:吉林大学,2016.Copyright©博看网. All Rights Reserved.。
