犁耕作业大马力拖拉机驱动轮滑转率控制方法

张 硕，武仲斌，陈 军，李 臻，朱忠祥，宋正河，毛恩荣

·农业装备工程与机械化·

犁耕作业大马力拖拉机驱动轮滑转率控制方法

张 硕1，武仲斌2，陈 军1，李 臻2※，朱忠祥2，宋正河2，毛恩荣2

（1. 西北农林科技大学机械与电子工程学院，杨凌 712100；2. 中国农业大学工学院现代农业装备优化设计北京市重点实验室，北京 100083）

针对犁耕作业时大马力拖拉机驱动轮易产生过度滑转的问题，该研究以大马力拖拉机电液悬挂机组为研究对象，考虑“拖拉机-农具-土壤”系统的强非线性特征，在建立大马力拖拉机犁耕作业机组非线性系统动力学模型的基础上，提出基于滑模变结构控制的大马力拖拉机驱动轮滑转非线性控制方法；并以模糊PID控制为对比，采用Matlab/Simulink验证本文动力学模型的正确性和控制算法的有效性；以Lovol-TG1254型大马力拖拉机为载体，搭建犁耕作业大马力拖拉机驱动轮滑转控制平台，开展田间对比试验，并分析不同控制方法下的滑转控制效果，验证滑模变结构控制算法的控制精度和稳定性。试验结果表明：在2.17 m/s的犁耕作业工况下，与模糊PID控制算法相比，滑模变结构控制算法将拖拉机驱动轮滑转率有效控制在最优值0.2，平均绝对值偏差为0.008，减小了约27%，最大偏差为0.028，减小了约49%；耕深、液压缸位移和水平牵引力调节变化量分别减小了27%、36%、42%。该研究提出的基于滑模变结构的大马力拖拉机驱动轮滑转控制方法可实现犁耕作业驱动轮滑转最优目标控制。

农业机械；试验；拖拉机；犁耕作业；滑转控制；滑模变结构控制

0 引 言

大马力拖拉机电液悬挂机组的犁耕作业是北方大田作物田间生产的重要农艺环节，拖拉机驱动轮的滑转程度极大地影响作业机组的作业效率和耕地质量[1-2]。随着农艺农机融合对土壤保护要求和大中型拖拉机智能化程度的不断提高[3-4]，大马力拖拉机田间作业的驱动防滑逐渐成为电液悬挂作业机组自动控制领域的研究热点[5]。

目前，国内外针对大马力拖拉机电液悬挂自动控制系统的研究主要涉及耕深或牵引力-耕深控制等[6-8]，少部分研究引入滑转率作为控制参量，采用PID[9]、模糊控制等方法[10-11]，其本质仍是依靠耕深调节实现将滑转率限定在某一阈值或区间的等效控制。不同于汽车等道路行驶车辆，犁耕作业时拖拉机驱动轮的过度滑转主要来源于“拖拉机-农具-土壤”系统的耦合作用，具有非常明显的非线性特征[12-13]。铧式犁等农具在与土壤接触、摩擦并发生剪切的过程中，受耕深变化和田间地面附着条件影响负载较大波动，加之地面附着条件不能产生足够支撑驱动轮前进的附着力，驱动轮就会发生过度滑转。在现有研究中，大多以线性化方法近似描述拖拉机、电液悬挂系统或作业机组的非线性特征，部分涉及非线性建模的研究也做出了较多假设和简化，普遍适用性不够[13-14]。也有部分研究建立作业机组数学模型以及应用非线性控制方法开展更为精确的驱动轮滑转控制，但是并没有充分考虑大马力拖拉机作业机组系统的强非线性特征，控制精度较低。

滑模变结构控制（Sliding Mode Variable Structure Control，SMVSC）是现代控制理论中一种较为成熟的非线性控制方法[15-16]，在农业机械自动控制领域已有效解决拖拉机主动减振[17]、路径跟踪控制[18]和自动导航[19]等非线性控制问题，提高了系统抵抗不确定参数扰动和外界激励的鲁棒性；也在电动车辆[20]、月球车[21]和特种装备[22]的驱动防滑中得到了较为成熟的应用。本文以“拖拉机-农具-土壤”系统为研究对象，在建立作业机组非线性系统动力学模型的基础上，提出了一种基于滑模变结构的大马力拖拉机犁耕作业驱动轮滑转控制方法实现面向犁耕作业工况的拖拉机驱动轮防滑控制，以提高大马力拖拉机电液悬挂机组的作业效率和耕地质量。

1 大马力拖拉机电液悬挂作业机组动力学建模

在大马力拖拉机电液悬挂作业机组数学建模过程中，考虑拖拉机在田间作业多为直线运动，且北方平原地区田间地势相对平坦，将其简化为沿拖拉机前进方向上的直线运动，忽略拖拉机前进过程中的侧滑、侧倾等侧向运动。如图1所示，以拖拉机后轮轴心点为坐标原点，建立拖拉机电液悬挂作业机组运动坐标系O-O。将大马力拖拉机电液悬挂作业机组的运动简化分解为2部分：一部分为拖拉机随基点沿前进方向上的平移运动，另一部分为拖拉机悬挂机构绕基点的转动运动，并以拖拉机在水平地面行驶时的相对位置关系来表达拖拉机各主要点的相对位置。由刚体平面运动学的基本理论可知，对于拖拉机及其悬挂机构杆件的各主要节点，其速度可表达为该点与基点在连线方向上的速度投影分量，其加速度可表达为基点的加速度与该点绕点的转动加速度的矢量和。

分别对各主要节点进行运动学分析，可得到拖拉机各主要节点的运动学关系。其中，后轮驱动拖拉机的驱动轮滑转率K如式（1）所示。

式中Ot、O分别为驱动轮沿水平前进方向的理论速度和瞬时速度，m/s；O、K分别为驱动轮的几何半径和动力半径，m；K为驱动轮角速度，rad/s。

注：O-O为拖拉机运动坐标系；为后轮轴心；1为前轮轴心；为提升液压缸与拖拉机的铰接点，、分别为上、下拉杆与拖拉机的铰接点，为提升液压缸活塞杆末端，、为提升杆的铰接点，为提升液压缸的提升轴轴心，、为农具的上下挂接点，为五铧犁的铧尖，T为拖拉机的质心，为农具质心；T为拖拉机俯仰角，rad；β为连线与垂直方向的夹角，rad；O1、O为前、后轮的几何半径，m；K1、K为前、后轮的动力半径，m；O1、O为垂直载荷引起的前、后轮着地点处土壤沉陷量，m；K1、K分别为前、后轮轮胎的弹簧刚度，N·mm-1；K1、K分别为前、后轮轮胎的阻尼系数，N·(m·s-1)-1；T为拖拉机相对其质心的转动惯量，N·m²；J为农具相对其质心的转动惯量，N·m²；T为拖拉机的重力，N；mg为农具的重力，N；K、K、K1分别为后轮所受的垂直地面反作用力、水平滚动阻力和驱动力，N；K1、K1、K1分别为前轮所受的垂直地面反作用力、水平滚动阻力和驱动力，N；H、V分别为农具在水平和竖直方向上所受的土壤阻力，N。

Note:O-Ois the tractor motion coordinate system;isthe axis of the rear wheel;1isthe axis of the front wheel;is the hinge point between the lifting hydraulic cylinder and the tractor;,is the hinge point between the upper and lower pull rods and the tractor;is the end of the piston rod of the lifting hydraulic cylinder;,is the hinge point of the lift lever;is the lifting shaft axis of the lifting hydraulic cylinder;,is the upper and lower contact of the plough;is the tip of the five-furrow plough;Tis the centroid of the tractor;is the center of mass of agricultural tools;Tis the pitching angle of tractor, rad;βis the angle betweenand the vertical direction, rad;O1,Oare the geometric radius of the front and rear wheels, m;K1,Kare the power radius of the front and rear wheels, m;O1,Oare the soil subsidencecaused by vertical loading at front and rear wheel locations，m；K1,Kare the spring stiffness of the front and rear tires respectively, N·mm-1;K1,Kare the damping coefficients of front and rear tires, N·(m·s-1)-1;Tis the moment of inertia of the tractor around the center of mass, N·m²;Jis the moment of inertia of the plough around the center of mass，N·m²;Tis the force of gravity on the tractor, N;is the force of gravity on the farm tools, N;K,K,Kare the vertical ground reaction force, horizontal rolling resistance and driving force of the rear wheel respectively, N;K1,K1,K1are the vertical ground reaction force, horizontal rolling resistance and driving force of the front wheel respectively, N;H,Vare soil resistance of farm tools in horizontal and vertical directions respectively, N.

图1 拖拉机电液悬挂作业机组运动学及动力学分析简图

Fig.1 Simple diagram of kinematic and dynamic analysis for high-power tractor

在图1中，对拖拉机电液悬挂作业机组进行受力分析可知，拖拉机（包括悬挂机构）所受外力主要有：拖拉机的重力T（为重力加速度，m/s2），农具的重力mg，前轮所受的垂直地面反作用力K1、水平滚动阻力K1和驱动力K1，后轮所受的垂直地面反作用力K、水平滚动阻力K和驱动力K，农具所受的土壤阻力H和V，农具对悬挂机构在上、下悬挂点处的作用力F、Fx和Fy。因此，以后驱动轮几何中心为中心，可得到水平和垂直方向的力平衡方程和力矩平衡方程，如式（2）所示。

式中Tx、Ty分别为拖拉机俯仰角在水平和垂直方向上的分量，rad；Wx、Wy分别为农具质心加速度在水平和垂直方向上的分量，m·s-2；P、P分别为点到后轮轴心的水平和垂直距离，m；sV、sH为犁体受力点到点处的水平和竖直距离，m；O1、O1前轮轴心到后轮轴心的水平和垂直距离，m；T为拖拉机质心到后轮轴心的水平距离，m；W、W为农具质心到后轮轴心在水平和垂直方向上的距离，m。

根据以上分析，并参考前期研究成果和建模方法[17,23]，将拖拉机车体、车轮、悬挂机构和农具等各子系统方程带入式（2），可得到大马力拖拉机电液悬挂作业机组动力学微分方程组，如式（3）所示。

2 滑模变结构控制算法

2.1 总体方案

在建立的大马力拖拉机电液悬挂作业机组驱动轮滑转非线性系统动力学模型中，存在部分难以准确获得的时变运动参数和不确定项，且系统方程存在强非线性。当以最优滑转率为控制目标，实现大马力拖拉机犁耕作业驱动轮滑转自动控制时，需采用一种适应强非线性特点、对外界扰动不敏感的非线性控制方法。在运动控制系统的非线性控制方法中，滑模变结构控制方法对模型的不确定性和外界随机具较高的可靠性和鲁棒性和[15-16]。如图2所示，采用滑模变结构控制方法设计了最优目标滑转率控制系统。

图2 基于滑模变结构控制的驱动轮滑转控制原理简图

2.2 滑模控制器设计

在北方平原地区小麦、玉米等大田作物田间犁耕作业时，在麦茬、玉米茬轮耕地块土壤条件下，当驱动轮滑转率维持在0.2附近时，拖拉机作业机组工作效率和耕地质量均能较好地满足农艺要求[24-25]。因此，本文以最优滑转率opt=0.2作为控制目标，定义控制系统输出误差如式（4）所示。

滑模变结构控制方法是控制输出误差趋近0，由式（3）可知，滑转控制系统动力学模型为一阶非线性系统，将滑模变结构控制关于时间的切换函数定义为

式中为滑模系数，且＞0，通过多次仿真确定=1。

在由滑转率及其1阶导数构成的滑模控制相平面中，切换线是以−为斜率，并且过（opt，0）的直线。滑模变结构控制可通过选取等效控制量，使得拖拉机在作业切换线滑动趋近控制目标（opt，0）。

对式（5）求导，根据广义滑模条件可知：

将式（3）中驱动轮滑转率的1阶微分方程表达式代入式（6）可得：

控制系统从任意初始状态不断趋近滑模面时，采用指数趋近律来改善趋近运动的动态效果。同时，在()持续接近0时，引入等速趋近律，确定趋近速度为某一非0等速度，以保证控制系统以某一速度尽快到达滑模面。所采用的指数趋近律如式（8）所示。

式中表示系统的运动点的等速趋近律，且＞0；表示指数趋近律，且＞0；sgn(m())表示符号函数。

联立式（7）和式（8），并选取提升液压缸所需油腔内活塞压力LL作为等效控制量，可得到：

式中1、2、3、4、5、6、7为拖拉机结构参数相关的系数，其计算如式（11）所示。

根据指数趋近律即（8），采用正定函数作为Lyapunov函数，具体可表达为

=2/2 （12）

式中表示滑模控制器的滑模面。

对式（12）求时间导数有：

式中为任意小的正整数，>0。

根据Lyapunov稳定性理论可知，在平衡点的邻域内，Lyapunov函数为正定，且其时间导数为负定，则系统局部渐进稳定，所设计的控制系统是稳定的。

当系统进入滑模运动并沿着滑模面不断切换时，为减少系统“抖动”的影响，采用连续饱和函数sat(/)代替传统等速趋近率不连续的符号函数sgn(m())，其定义如式（14）所示。

式中表示边界层厚度。

2.3 控制量的等效转换

在滑模变结构控制算法的设计过程中，推导了等效控制量LL的表达式。电液悬挂自动控制系统的控制执行机构采用自行研发的螺纹插装式比例控制阀（包括比例提升阀和比例下降阀，最高工作压力25 MPa，额定流量为80 L/min），通过比例放大器（驱动器）的驱动电压，实现阀芯开度和流量调节，最终达到等效控制量即提升液压缸所需油腔活塞压力的理想值。在室内试验中，根据比例控制阀的动静态特性，调节连接比例控制阀出油口的比例溢流阀开启压力来模拟负载变化，改变比例放大器（驱动器）的输入电压，得到不同驱动电压下比例控制的稳态流量。根据试验结果进行回归分析，可得到控制电压与输出流量的拟合关系如式（15）~（16）所示。

式中up为比例提升阀的驱动电压，V；down为比例下降阀的驱动电压，V；up为液压缸提升过程中，从比例提升阀流进液压缸的液压油流量，m3/s；down为液压缸下降过程中，从液压缸流入比例下降控制阀的流量，m3/s；Δ为比例换向阀两端压差，取值1.5 MPa；L为液压缸负载压力，Pa；0为回油压力，取值0；为液压油液的密度，kg/m3。

提升液压缸的流量连续性方程如式（17）所示。

式中为油液的体积弹性模量，Pa；L为提升液压缸油腔活塞的有效作用面积，m2；L为提升液压缸油腔及油路中的油液容积，m3；tL为提升液压缸的外泄漏系数，m3/(Pa·s)。

联立式（10）和式（16）求得L。提升阀开启时L=up()；下降阀开启时L=−down()。综合式（15）~（16），最终得到控制系统输入信号即电液比例控制阀的控制电压。

3 控制算法仿真分析

为了验证所建立的非线性动力学模型的准确性和所设计的控制算法的有效性，在Matlab/Simulink中建立犁耕作业大马力拖拉机驱动轮滑转控制系统仿真模型，如图3所示。在仿真模型中，根据北方平原地区小麦、玉米等大田土壤条件[24-25]，采用滤波白噪声法并根据拖拉机运动过程中前后轮的时间延迟关系，建立处于E级与F级标准路面之间的田间随机路面激励。

为了验证控制算法有效性，以最优滑转率opt=0.2为控制目标，分别采用SMVSC和模糊PID控制算法进行仿真，分析滑转控制系统对外界扰动即土壤比阻输入变化的消扰特性和动态响应特性。经过多次调整，确定SMVSC相关控制参数为=1.2、=0.01、=0.01，模糊PID控制初始参数为P=12、I=0.5、D=1.5。仿真初始条件设置为拖拉机车速2.2 m/s，液压缸活塞杆的初始位移量8.83 cm，初始耕深20 cm，仿真时间50 s。将土壤比阻的变化作为外界输入扰动，设定土壤比阻稳态值为30 000 N/m2[9,24]。仿真开始时，加入为幅值8 000 N/m2、周期10 s、占空比50%的土壤比阻扰动脉冲信号，仿真结果如图4所示。

当仿真开始时，土壤比阻从30 000 N/m2阶跃变化到38 000 N/m2，随后以10 s周期方波持续变化。由图4a可知，采用模糊PID控制时，当土壤比阻每隔5 s阶跃变化时，系统响应时间为1.5 s左右，虽然能将滑转率控制在0.2附近，但滑转率波动较大，最大误差达到0.02；SMVSC控制下的系统响应时间约为0.5 s，在控制初期超调量几乎为0，在0.5 s后能持续稳定追踪控制目标0.2，能较好地抵抗土壤比阻扰动变化的影响。另外，在土壤比阻处在30 000 N周期内时，由于加入幅值为8 000 N/m2的随机路面土壤比阻激励扰动，模糊PID控制根据偏差的大小时刻进行调整，驱动轮滑转率的稳态误差达到0.015，最大超调量0.018；而SMVSC控制下的驱动轮滑转率始终在滑模面附近，稳态误差0.002，对外界扰动的鲁棒性较好，能持续稳定在0.2，最大超调量0.008。

由图4b、图4c、图4d可知，在SMVSC控制系统中，土壤比阻阶跃变化时，液压缸活塞杆位移量由8.99 cm伸出至9.56 cm，耕深由21 cm减小至16.21 cm，通过减小耕深来抑制滑转率增大的趋势，水平牵引力稳定在7 056.54 N左右；模糊PID控制系统的液压缸活塞杆位移瞬间增大，直接导致耕深变化和滑转率控制的超调量均较大，响应时间较长。

图3 大马力拖拉机犁耕作业驱动轮滑转控制系统仿真模型

图4 大马力拖拉机犁耕作业驱动轮滑转控制仿真结果

由此可见，在土壤条件发生改变时，和模糊PID控制相比，滑模变结构控制对外界扰动的消扰特性较好，响应相对较快，验证了其控制的有效性和优越性。

4 田间试验

4.1 试验平台与试验方案

综合大马力驱动轮滑转控制系统的实际应用需求，以Lovol-TG1254型大马力拖拉机为载体，挂接东方红1LH-535型五铧犁搭建田间试验平台，如图5所示。采用Radar III型地面多普勒雷达测量拖拉机实际前进速度，脉冲计数式发动机转速传感器测量理论车速，并通过计算得到驱动轮实时滑转率；采用Wittower RE-38型绝对值型角度传感器测量提升臂转角，依据悬挂机构杆件关系计算得到耕深[9]；将应变片式牵引力传感器阻安装在两侧下拉杆与拖拉机的铰接点处，测量农具在水平和竖直方向上所受到的牵引阻力；采用文献[24-25]的经验公式，通过实时滑转率估算驱动力系数并计算驱动力。

2019年5月，在中国农业大学上庄试验站内典型北方大田轮种地块机型田间试验。试验地块大小60 m× 350 m，地表以下20 cm土壤坚实度平均约为130.5 kPa，拖拉机车速为B2挡、2.17 m/s。地块土壤硬度变化较大且地表杂草较多，在播种期浇水后的犁耕作业过程中，拖拉机驱动轮易产生过渡滑转。基于所搭建的田间试验平台，分别开展SMVSC和模糊PID控制算法下的犁耕作业大马力拖拉机驱动轮滑转控制对比试验。

1.多普勒雷达 2.发动机转速传感器 3.上位机控制界面 4.控制器 5.供电系统和驱动器 6.电液比例控制阀 7.耕深传感器 8.牵引力传感器

4.2 试验结果与分析

在试验过程中，首先由驾驶员操纵电液悬挂液压输出操纵杆，迫使犁具从地面以上以一定角度快速入土，在拖拉机作业机组开始犁耕作业以后，打开滑转自动控制系统。以0.2为最优控制目标，分别开展了基于SMVSC控制算法和模糊PID控制算法的驱动轮滑转控制对比试验，试验结果如图6所示。

由图6a可知，SMVSC控制下的滑转率平均值为0.2013，最大偏差为0.028，平均绝对值偏差为0.008，方差为0.000 1；模糊PID控制下的滑转率平均值为0.2045，最大偏差为0.055，平均绝对值偏差为0.011，方差为0.002。SMVSC控制下的滑转率平均绝对值偏差减小了27%，最大偏差减小约49%，滑转率方差也显著减小。虽然2种控制方法都能将滑转率控制在最优目标0.2附近，但是SMVSC控制控制下的滑转率波动幅度和控制偏差均相对较小，稳定性更好。

由图6b、图6c、图6d可知，SMVSC控制下的耕深调节的变化量为5.26 cm，远小于模糊PID控制的7.21 cm，减小了约27%；液压缸位移调节变化量为1.15 cm，小于模糊PID控制的1.8 cm，减小了约36%；水平牵引力的调节变化量为1 293.35N，小于模糊PID控制的2 217 N，减小了约42%。由此可见，基于滑模变结构的驱动轮滑转控制方法将滑转率稳定控制在0.2的同时，拖拉机作业过程中的调整量更小，作业状态更加平稳。

图6 大马力拖拉机犁耕作业驱动轮滑转控制试验结果

对比仿真和试验结果可知，在耕深20 cm的田间犁耕作业工况下，采用本文所设计的基于滑模变结构控制的大马力拖拉机驱动轮滑转控制方法，能够有效将拖拉机驱动轮滑转率控制在最优目标0.2。与模糊PID控制方法相比，控制超调量较小，达到稳定状态以后控制偏差较小，作业过程中的耕深调整量更小，控制精度和稳定性更高。

5 结 论

1）面向犁耕作业工况，以大马力拖拉机电液悬挂作业机组为研究对象，考虑“拖拉机-农具-土壤”复杂系统的强非线性特征，建立了适用于的大马力拖拉机驱动轮滑转的非线性系统动力学模型；在对动力学微分方程组和液压系统进行有效简化的基础上，选取提升液压缸油腔内活塞压力作为等效控制量，采用滑模变结构控制理论设计了驱动轮滑转非线性控制算法。

2）以最优滑转率0.2为控制目标，以模糊PID控制作为对比算法，应用Matlab/Simulink进行土壤比阻阶跃变化激励下的仿真分析；并搭建了滑转控制系统田间试验平台，开展大马力拖拉机犁耕作业驱动轮滑转控制田间对比试验。仿真和田间试验结果表明，基于滑模变结构控制的驱动轮滑转控制方法的控制误差更小，控制精度和稳定性更高。

该文提出的基于滑模变结构控制的大马力拖拉机犁耕作业驱动轮滑转控制方法，在实现滑转率最优目标稳定控制的同时，拖拉机电液悬挂系统调整幅度较小，作业过程更加平稳，对于提高田间犁耕作业质量和作业效率具有一定的实际生产指导意义。

[1] 张向前，杨文飞，徐云姬．中国主要耕作方式对旱地土壤结构及养分和微生态环境影响的研究综述[J]. 生态环境学报，2019，28(12)：2464-2472.

Zhang Xiangqian, Yang Wenfei, Xu Yunji. Effects of main tillage methods on soil structure, nutrients and micro-ecological environment of upland in China: A review[J]. Ecology and Environmental Sciences, 2019, 28(12): 2464-2472. (in Chinese with English abstract)

[2] 关劼兮，陈素英，邵立威，等．华北典型区域土壤耕作方式对土壤特性和作物产量的影响[J]. 中国生态农业学报，2019，27(11)：1663-1672.

Guan Jiexi, Chen Suying, Shao Liwei, et al. Soil tillage practices affecting the soil characteristics and yield of winter wheat and summer maize in North China[J]. Chinese Journal of Eco-Agriculture, 2019, 27(11): 1663-1672. (in Chinese with English abstract)

[3] 赵建军. 重型拖拉机电液提升器比例控制阀设计与特性研究[D]. 北京：中国农业大学，2015.

Zhao Jianjun. The Design and Characteristics Research of Proportional Control Valve for the Electro-hydraulic Hitch of Heavy Tractor[D]. Beijing: China Agricultural University, 2015. (in Chinese with English abstract)

[4] 丁肇，李耀明，唐忠. 轮式和履带式车辆行走对农田土壤的压实作用分析[J]. 农业工程学报，2020，36(5)：10-18.

Ding Zhao, Li Yaoming, Tang Zhong. Compaction effects of wheeled vehicles and tracked on farmland soil[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(5): 10-18. (in Chinese with English abstract)

[5] 倪向东，梅卫江. 国外大功率拖拉机新技术新结构和发展趋势的研究[J]. 农机化研究，倪向东，梅卫江. 国外大功率拖拉机新技术新结构和发展趋势的研究[J]. 农机化研究，2007，1：25-27.

Ni Xiangdong, Mei Weijing. Research on the new structure and technology of the development trend for the foreign large power tractor[J]. Journal of Agricultural Mechanization Research, 2007, 1: 25-27. (in Chinese with English abstract)

[6] 李立. 拖拉机后悬挂电液控制系统的研究[D]. 哈尔滨：哈尔滨工业大学，2012.

Li Li. Research on the Electro-hydraulic System in Tractor Rear Suspension[D]. Harbin: Harbin Institute of Technology, 2012. (in Chinese with English abstract)

[7] 王会明，侯加林，赵耀华，等. 拖拉机液压悬挂机构自动控制系统[J]. 农业机械学报，2006，37(10)：42-45.

Wang Huiming, Hou Jialin, Zhao Yaohua, et al. Study on automatic control system for hydraulic hitch equipment of tractor[J]. Transactions of the Chinese Society for Agricultural Machinery, 2006, 37(10): 42-45. (in Chinese with English abstract)

[8] Chanchal Gupta, Tewari V K, Ashok Kumar A, et al. Automatic tractor slip-draft embedded control system[J]. Computers and Electronics in Agriculture, 2019, 165: 1-11.

[9] 张硕，杜岳峰，朱忠祥，等. 后轮驱动大功率拖拉机牵引力-滑转率联合自动控制方法[J]. 农业工程学报，2016，32(12)：47-53.

Zhang Shuo, Du Yuefeng, Zhu Zhongxiang, et al. Integrated control method of traction & slip ratio for rear-driving high-power tractors[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2016, 32(12): 47-53. (in Chinese with English abstract)

[10] 白学峰，鲁植雄，张广庆，等. 基于滑转率的拖拉机驱动防滑模糊PID控制算法仿真分析[J]. 江西农业学报，2012，9：146-149，156.

Bai Xuefeng, Lu Zhixiong, Zhang Guangqing, et al. Simulative analysis of fuzzy PID control algorithm for antiskid driving system of tractor based on slip ratio[J]. Acta Agriculturae Jiangxi, 2012, 9: 146-149, 156. (in Chinese with English abstract)

[11] 殷新东，鲁植雄. 基于滑转率的四轮驱动拖拉机防滑模糊控制算法仿真研究[J]. 农业装备与车辆工程，2010，12：6-10.

Yin Xindong, Lu Zhixiong. Simulation research on acceleration slip regulation system for four-wheel drive tractor using fuzzy control method[J]. Agricultural Equipment & Vehicle Engineering, 2010, 12: 6-10. (in Chinese with English abstract)

[12] Mirko Simikic, Nebojsa Dedovic, Lazar Savin, et al. Power delivery efficiency of a wheeled tractor at oblique drawbar force[J]. Soil & Tillage Research, 2014, 141:32-43.

[13] 高大水，余群，喻谷源. 驱动轮滑转沉陷及其对牵引性能的影响[J]. 北京农业机械化学院学报，1984(3)：4-27.

Gao Dashui, Yu Qun, Yu Guyuan. Research of slip-sinkage of tractor driving wheels and its effect on tractor performance[J]. Journal of Beijing Agricultural Mechanization College, 1984(3): 4-27. (in Chinese with English abstract)

[14] Porteš P, Bauer F, Čupera J. Laboratory-experimental verification of calculation of force effects in tractor's three-point hitch acting on driving wheels[J]. Soil & Tillage Research, 2013, 128: 81-90.

[15] 刘金琨，孙富春. 滑模变结构控制理论及其算法研究与进展[J]. 控制理论与应用，2007，24(3)：407-418.

Liu Jinkun, Sun Fuchun. Research and development on theory and algorithms of sliding mode control[J]. Control Theory & Applications, 2007, 24(3): 407-418.

[16] 张丰. 滑模变结构控制理论在非线性系统中的应用[D]. 沈阳：沈阳理工大学，2013.

Zhang Feng. The Application of Sliding-Mode Variable Structure Control Theory in Nonlinear Systems[D]. Shenyang: Shenyang Ligong University, 2013. (in Chinese with English abstract)

[17] 承鉴，迟瑞娟，毛恩荣. 悬挂农具对电液悬挂系统拖拉机振动的影响[J]. 农业工程学报，2015，31(7)：24-32.

Cheng Jian, Chi Ruijuan, Mao Enrong. Influence of hanging farm implement on vibration of tractor with electro-hydraulic hitch system[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(7): 24-32. (in Chinese with English abstract)

[18] 牛雪梅，高国琴，鲍智达，等. 基于滑模变结构控制的温室喷药移动机器人路径跟踪[J]. 农业工程学报，2013，29(2)：9-16.

Niu Xuemei, Gao Guoqin, Bao Zhida, et al. Path tracking of mobile robots for greenhouse spraying controlled by sliding mode variable structure[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2013, 29(2): 9-16. (in Chinese with English abstract)

[19] 赵翾，杨珏，张文明，等. 农用轮式铰接车辆滑模轨迹跟踪控制算法[J]. 农业工程学报，2015，31(10)：198-203.

Zhao Xuan, Yang Jue, Zhang Wenming, et al. Sliding mode control algorithm for path tracking of articulated dump truck[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2015, 31(10): 198-203. (in Chinese with English abstract)

[20] 武仲斌，谢斌，迟瑞娟，等. 基于滑转率的双电机双轴驱动车辆转矩协调分配[J]. 农业工程学报，2018，34(15)：66-76.

Wu Zhongbin, Xie Bin, Chi Ruijuan, et al. Active modulation of torque distribution for dual-motor front and rear-axle drive type electric vehicle based on slip ratio[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2018, 34(15): 66-76. (in Chinese with English abstract)

[21] 黄郁馨，王彤宇，林琳，等. 基于滑模变结构控制的移动系统滑转率控制[J]. 兵工学报，2014，35(10)：1707-1715.

Huang Yuxin, Wang Tongyu, Lin Lin, et al. Slip ratio control of locomotion system based on sliding mode variable structure control[J]. Acta Armamentarii, 2014, 35(10): 1707-1715. (in Chinese with English abstract)

[22] 周兵，徐蒙，袁希文，等. 基于滑模极值搜索算法的车辆驱动防滑控制策略[J]. 农业机械学报，2015，46(2)：307-311.

Zhou Bing, Xu Meng, Yuan Xiwen, et al. Acceleration slip regulation based on extremum seeking control with sliding mode[J]. Transactions of the Chinese Society for Agricultural Machinery, 2015, 46(2): 307-311. (in Chinese with English abstract)

[23] 张硕. 基于滑模变结构的重型拖拉机犁耕作业滑转率控制方法研究[D]. 北京：中国农业大学，2018. Zhang Shuo. Study on Slip Rate Control of Heavy Tractor for Ploughing Based on Sliding Mode Variable Structure Control[D]. Beijing: China Agricultural University, 2018. (in Chinese with English abstract)

[24] 周志立，方在华. 拖拉机机组牵引动力学[M]. 北京：科学出版社，2010.

[25] 吴起亚，高行方. 拖拉机与农业机械牵引力学[M]. 北京：中国农业机械出版社，1985.

Control method of driving wheel slip rate of high-power tractor for ploughing operation

Zhang Shuo1, Wu Zhongbin2, Chen Jun1, Li Zhen2※, Zhu Zhongxiang2, Song Zhenghe2, Mao Enrong2

(1.,712100,; 2.,,100083,)

In north China, the ploughing operation of high-power tractor based on electro-hydraulic suspension system is the most common and important agricultural process in the field agricultural production. Due to the complex field working environment, the change of tillage depth and the fluctuation of soil specific resistance, the working load of tractor unit fluctuates greatly during ploughing operation, which is easy to cause excessive sliding of driving wheel， and seriously affects the traction efficiency and traction of the operation unit. In order to solve the problem of excessive driving wheel slip of high-power tractor for ploughing, taking the high-power tractor ploughing unit with electro-hydraulic hitch system as the research object, a sliding rate control method based on sliding mode variable structure control was proposed.this paper. Firstly, in view of the strong nonlinear characteristics of the complex system of “tractor-farm tools-soil”, the nonlinear dynamics model of tractor driving wheel sliding for ploughing was established based on tractor motion characteristics and the theories of vehicle dynamics. Then, the exponential reaching law was used to design the sliding rate nonlinear control algorithm based on sliding mode variable structure control theory. Especially, taking the hydraulic pressure of hydraulic cylinder as the equivalent control quantity, the control law of sliding mode variable structure control was derived after simplifying the hydraulic system. Through MATLAB/Simulink simulation, the reliability of the nonlinear dynamic model and control algorithm in the step change of soil specific resistance was verified. According to the soil data obtained from a specific site, the soil specific resistance was set as pulse signal, and the stable value was 30 000 N/m2. At the beginning of the simulation, a soil specific resistance with an amplitude of 8 000 N/m2, a period of 10 s and a duty cycle of 50% was added to block the impulse signal. The results showed that the slip rate of the driving wheel was always near the sliding surface, the steady-state error was 0.002, and the control overshoot was 0.008. Furthermore, the field test platform for automatic slip rate control system was built on Lovol-TG1254 tractor, and the field comparison tests between SMVSC control and Fuzzy PID control were carried out under the tractor speed of 2.17 m/s. The test results showed that compared with the fuzzy PID control, the mean absolute deviation of slip rate under the SMVSC control decreased by 27%, the maximum deviation decreased by about 49%, and the fluctuation range and control deviation of slip rate were smaller. The change of depth was 5.26 cm, which decreased by about 27% compared with that of the fuzzy PID control, the change of hydraulic cylinder displacement was 1.15 cm, which decreased by about 36%, and the adjustment change of traction was 1 293.35 N, which decreased by about 42%. The control deviation of the driving wheel sliding control method proposed in this paper was small, it can provide theoretical basis and technical support for improving the quality of ploughing operation.

agricultural machinery; experiments; tractor; ploughing operation; slip control; sliding mode variable structure control

张硕，武仲斌，陈军，等. 犁耕作业大马力拖拉机驱动轮滑转率控制方法[J]. 农业工程学报，2020，36(15)：47-55.doi：10.11975/j.issn.1002-6819.2020.15.006 http://www.tcsae.org

Zhang Shuo, Wu Zhongbin, Chen Jun, et al. Control method of driving wheel slip rate of high-power tractor for ploughing operation[J]. Transactions of the Chinese Society of Agricultural Engineering (Transactions of the CSAE), 2020, 36(15): 47-55. (in Chinese with English abstract) doi：10.11975/j.issn.1002-6819.2020.15.006 http://www.tcsae.org

2020-04-16

2020-08-13

中国博士后基金面上项目（2019M653764）；国家重点研发计划项目（2017YFD0700403）

张硕，博士，讲师，主要从事车辆智能控制、智能农机装备的相关研究。Email：zhangshuo@nwafu.edu.cn

李臻，博士，副教授，主要从事农业机械的设计、仿真，以及自动控制研究。Email：zhenli@cau.edu.cn

10.11975/j.issn.1002-6819.2020.15.006

S219.032.4; TP273+.2

A

1002-6819(2020)-15-0047-09