This paper presents a generalized, multibody dynamics model for a tracked vehicle equipped with a towing winch and control strategies that enhance vehicle mobility by regulating track slip based on real-time terrain characterization and automating winch use. The vehicle model is validated under conditions where no action is taken by the winch. Thereafter, two mobility enhancing control strategies are outlined. The first strategy regulates track slip to a real-time estimated value that generates maximum net traction. This is done by computing state-force estimates from a Kalman filter that are compared to terrain traction models using a Bayesian hypothesis selection approach. If the vehicle is traction limited and the first strategy fails, a second strategy that automates winch use is activated. Simulation results are shown for both scenarios.

## Introduction

Tracked vehicles have good performance in off-road locomotion over soft terrain when compared to wheeled vehicles. This is because the vehicle weight is distributed over a large area, which reduces sinkage and increases traction [1,2]. The demand for reliable autonomy for heavy class tracked vehicles is growing for operations involving agriculture, construction, and logistics. These tasks require towing heavy loads where immobilization can lead to long down times.

To maintain mobility autonomously, two approaches are proposed in this paper. The first is a traction control architecture outlined in Fig. 1. Here, a discrete time Kalman filter (DTKF) computes state-force estimates and compares these to terrain traction models stored in memory prior to operation. The mostly likely models are selected in a Bayesian hypothesis selection method. Based on the models selected, a “peak slip” point is approximated where the vehicle can operate with a track slip ratio that generates maximum net traction. This is the reference input to a traction control system that uses a proportional-integral-derivative (PID) controller to regulate track slip. Individual pieces of such an architecture have been proposed and implemented via simulation, but none implement traction control using a validated model that incorporates vehicle-terrain dynamics and limitations of the vehicle powertrain as presented in this paper. Related work in terrain characterization has focused on lighter weight, wheeled platforms [3–7] where the underlying terrain models assumed do not hold well compared with heavy-duty vehicles considered here. The approach of computing state-force estimates and comparing them to traction models stored offline to characterize terrain is similar to Ref. [6]. However, the formulation here uses a DTKF instead of an extended Kalman filter (EKF) so that the filter can run with fixed gain to reduce the real-time computation cost. Furthermore, we evaluate 320 terrain traction models in the Bayesian hypothesis selection where only 21 were evaluated in Ref. [6]. Previous work in traction controllers for tracked vehicles is limited. Fan derives a sliding mode controller and proves convergence via Lyaponuv stability theory [8,9]. However, it is implied that the track slip ratio can be measured perfectly for state feedback and does not take into account realistic limitations of the vehicle powertrain as an actuator.

The second approach employed to maintain mobility autonomously uses a winch control architecture. A picture of a CASE Quadtrac equipped with a winch is shown in Fig. 2. Hydraulic winches allow heavy-duty vehicles to reduce loading at the tractor's hitch or drawbar temporarily in soft terrain by releasing cable. Once firmer terrain is reached, the winch reels cable in to recover the payload. The winch control approach here uses a low bandwidth logic controller based on slip estimates from the DTKF. This gives the traction control mode time to respond when the two control modes are active simultaneously so that it can respond to varying drawbar loading. Previous work in cable towing systems has been investigated for marine vehicles where one vehicle tows another [10–13]. In these studies, the path of the towed vehicle is the control variable of interest. Therefore, knowing the cable position and trajectory is essential for the desired control performance. Modeling of the cable includes finite-segment, rigid link lumped parameter elements. The control objective of interest motivating the modeling presented in this paper is mobility of the towing vehicle. Therefore, the simplified model of the cable here treats it as a rigid body that can only support tensile loads.

The work is motivated by the National Science Foundation South Pole Overland Traverse (SPoT) and Greenland Inland Traverse (GrIT) which use AGCO MT865 and CASE Quadtrac tractors, respectively. The traverses resupply research stations in Antarctica and Greenland, respectively, and optimizing the mobility of the tractor fleets has the potential to increase the efficiency of operations. The model and control development is for an MT865 tractor but the same approaches can be used for other vehicles of the same class.

The outline of this paper is as follows: Section 2 summarizes the vehicle-terrain model and validation of the model using experimental data from an AGCO MT865 tractor in Antarctica. Section 3 describes the tractor-winch-sled model, which incorporates hydraulic dynamics and hybrid dynamics of the winch system along with rigid-body sled dynamics. Section 4 describes an algorithm for real-time terrain characterization. Section 5 explains the design of the feedback traction controller. Section 6 describes criteria for entering and exiting the traction control modes based on the operating terrain parameter space. Section 7 shows simulation results of the traction control mode for five different terrain types to show the robustness of the control approach. Section 8 describes the design and justification for the rule set used in the logic-based winch controller. Sections 9 shows simulation results using both the traction and winch control modes, where the use of modes provides a robust mobility solution for heavy-tracked vehicles. Section 10 concludes the paper.

## Vehicle-Terrain Model Validation

*, normalized driver input throttle Π*

_{E}*, the transmission gear number or ratio*

_{I}*g*

_{GR}, the driver speeds $\phi \u02d9$, tractor ground speed

*v*, and drawbar load at the hitch DB. These values allow longitudinal motion validation of the model presented in Ref. [1]. The equations of motion reduce to

_{T}*F*is the gross traction generated by both tracks,

*R*is the resistance from the terrain imparted on the track, DB is the drawbar load,

*m*= 25,000 kg is the tractor's nominal mass,

_{T}*τ*is the torque of the driver sprocket in the track,

*ζ*= 10 (kg/m

^{2})/s a damping coefficient, and

*J*= 300 kg/m

_{S}^{2}is the track's moment of inertia. The equations for gross traction and terrain resistance come from Wong–Bekker terramechanics theory [1,14] and are given by

where *A *=* *2.97 m^{2} is the nominal contact area of one track, *c* is the terrain cohesion in kPa, *W* is the vehicle nominal weight in kN, Φ is the terrain friction angle, *l *=* *3.25 m is the nominal contact length a track, *i* is the track slip ratio, *K* is the terrain shear deformation modulus in cm, *k _{c}* is the cohesion modulus in kN/m

^{n}^{+1},

*k*

_{Φ}is the friction angle modulus in kN/m

^{n}^{+2},

*n*is the unitless exponent modulus,

*z*is the total track sinkage in m,

*S*is the unitless slip-sinkage coefficient, and

*z*is the static track sinkage in m. Equation (5) extends the original Wong–Bekker theory to account for slip-sinkage effects.

_{S}*z*is the static sinkage and the

_{S}*iSz*term predicts additional sinkage due to track slippage and associated terrain excavation.

_{S}*S*values used in this study are based on experimental data in Ref. [15].

*g*

_{GR}is the gear number,

*g*

_{FD}is the final drive gear ratio and

*ω*is the engine speed in rad/s [16]. The torque delivered to the tracks is a function of the engine speed and is calculated by

_{E}*τ*

_{E}_{,max}is the maximum available torque at the engine's current operating speed that is calculated by interpolation of tabulated torque-speed engine data provided in Ref. [16], and

*τ*

_{load}is the load for the winch hydraulics. An overriding engine controller, Π

*is modeled to prevent stalling. Π and Π*

_{C}*are given by MathWorks [17] and*

_{C}For a more detailed discussion of the coupling of the powertrain and vehicle terrain dynamics, refer to Refs. [18] and [19].

The model is validated by simulation where a tractor starts from rest or zero initial conditions at *t *=* *0 using the inputs of the driver during data collection of Π* _{I}* and

*g*

_{GR}and measured DB. Other inputs to the model for simulation are the estimated terrain parameters for compacted firn on which experiments were conducted:

*K*=

*2 cm,*

*c*=

*10 kPa, and Φ = 19 deg and the total resistance*

*R*

_{total}=

*R*+

*DB*. In Fig. 3, model predictions or simulation results are plotted against experimental data as blue and black lines, respectively, for 4 min of data during which the tractors starts from rest and ramps up to a steady velocity while shifting through gears. Measured and simulated responses for

*v*, $\phi \u02d9$,

_{T}*i*, Ω

*, and engine power*

_{E}*P*are provided in Fig. 3. The blue line on the Π plot is the output of the overriding engine controller, Π

_{E}*.*

_{C}The overriding engine controller overrides the driver input for the first 25–30 s in third gear and at several other points during gear shifts. This highlights the model's simplifying assumptions as presented, which assume constant track resistance forces. It is plausible that different terrain was encountered during the beginning of data collection such that this assumption does not hold for all time. In addition, it is postulated that the discrepancies during gear shifts can be accounted for by the use of oversimplified dynamics that do not account for control mechanisms. However, predictions at all other points prove to be accurate and show that the model sufficiently couples the powertrain with the vehicle-terrain dynamics.

## Tractor–Winch–Sled Model

Two important model elements are presented in this section to simulate the tractor–winch–sled system: the multibody vehicle dynamics that incorporate the winch and the hydraulics components that power it.

### Rigid Body Dynamics.

**u**= [

*τ*,

*τ*]

_{W}*, where*

^{T}*ψ*is the winch angular position, $\psi \u02d9$ is the winch angular velocity,

*v*

_{SD}is sled's velocity,

*τ*is the winch torque, and the remaining variables have been defined previously. The governing differential equations and constraint are given by

_{W}**M**

_{1}

where *J _{W}* = 0.82 kg

^{2}is the winch rotational inertia,

*r*= 0.19 m is the winch radius,

_{W}*ζ*= 3 (kg/m

_{W}^{2})/s is the winch damping coefficient,

*L*is the total cable available, and

_{C}*ψr*indicates the amount of cable let out by the winch.

_{W}**u**=

*τ*. The governing differential equations are

**M**

_{2}

### Winch Hydraulics.

*g*= 1, a fixed displacement piston motor connected to the winch through the gear ratio

_{P}*g*= 18, a four-way two-position valve in the center, a pressure relief valve B1, and a braking valve B2 [20–22]. The schematic for the proposed system is shown in Fig. 5 in two modes of operation: reeling in the winch and braking. In both modes, the kinematic constraints for the pump and motor speeds are

_{M}*P*and

_{H}*P*and activates the braking valve denoted as B2. The equations describing the pressure-rise rate within the section of hose

_{O}*P*in this mode of operation are given by [23]

_{O}The model used for the braking valve B2 is the same as for a pressure relief valve B1 to regulate hydraulic pressure. The valve set relief pressure, *P*_{set} is an input that modulates how much braking torque to use for the hydraulic motor in the winch.

*I*) isolates the sections of hose,

*P*and

_{H}*P*to create a pressure differential driven by the piston pump connected to the engine. The equations in this mode of operation are given by

_{O}This mode of operation for the hydraulics allows the winch to pull in the towed load. The value *P*_{set} in Eq. (30) is an input to the variable displacement, pressure-compensated pump to maintain a pressure set point for *P _{H}*. This modulates the amount of torque the motor uses to pull in the towed load.

## Real-Time Terrain Characterization

Characterizing terrain mobility conditions in real-time requires knowledge of track slip and terrain-vehicle forces. However, neither of these can be directly measured, and therefore, force and slip are estimated from predicted and measured trajectory motion using a Kalman filter. These estimates are then compared with terrain slip–force curve models stored prior to operation where the most likely terrain and its associated peak slip are identified recursively using Bayes rule.

*v*, $\phi \u02d9$, and

_{T}*DB*. By defining the net traction as

*F*

_{net}=

*F*–

*R*and the resistance torque at the driver sprocket as $\tau res=F+\zeta \phi \u02d9$, Eqs. (1) and (2) can be rewritten as

**K**

*. Therefore, the filter reduces to*

_{ss}The parameter estimate $p\u0302k$ along with the slip mesh previously defined are used to recompute Eqs. (3)–(6) where the slip value that corresponds to the maximum *F*_{net} value is referred to as the estimated peak slip $i\u0302pk$ and is the reference input for the traction controller outlined in Sec. 5.

This approach requires a priori knowledge about the range of terrain parameters expected since only a finite number of models can be considered. The a priori knowledge assumed here are bounds on the six-dimensional terrain parameter space defined by the terrain parameter vector $p=[c\Phi nkeqKS]T$. Discretization of the terrain parameter space is as follows: $c\u2208{1,2,\u2026,8},\u2009\Phi \u226120\u2009deg,\u2009n\u22611,\u2009(kc/b)+k\Phi \u2261keq\u2208{100,200,\u2026,500},\u2009K\u2208{0.5,1.5,\u2026,7.5},\u2009S\u2261(60%/33%)z$.

The justification for the values chosen here comes from Refs. [1], [2], [15], and [27]. Estimated values for Antarctic snow terrain cohesion are given in Ref. [27]. Data from Wong [1,2] show values slightly above and below 20 deg for Φ, and Lever estimates the range to be 15 deg ≤ Φ ≤ 25 deg [27]. Furthermore, *c* and Φ cannot be uniquely estimated or defined to produce the same traction forces due to the parameterization in Eq. (3), so there is no added benefit to discretizing both parameters and determining the likelihood of different combinations. The Bekker sinkage exponent *n* is defined to be constant as [27] estimates *n *=* *1 or close to this value for all Antarctic snow. *k*_{eq} is discretized to match the variety of different sinkage conditions experienced by SPoT defined by the Bekker relationship in Eq. (4). The discretization of *K* is larger due to its effect on the slip point at which the maximum net traction occurs referred to as the peak slip point, *i _{pk}*. In order to allow for a range of

*i*,

_{pk}*K*between 0.5 cm and 7.5 cm is considered. The last parameter

*S*is defined as a constant due to the large amount of experimental data in Ref. [15] that shows consistent dynamic slip-sinkage behavior.

## Gain-Scheduled Traction Control

A block diagram of the closed-loop longitudinal tractor dynamics under traction control is shown in Fig. 6. Here, the throttle Π in the input to the engine and gear ratio selection *g*_{GR} is the input to the transmission. The tractor speed *v _{T}* and driver speed $\phi \u02d9$ are the states of the rigid body dynamics. The engine is considered as a time-varying gain depending on the engine operating speed, and the transmission is a piecewise constant gain value that is the product of

*g*

_{GR}and

*g*

_{FD}. The approximate continuous range for the engine gain is 1750 ≤

*K*≤ 2400 and the product of

_{E}*g*

_{GR}

*g*

_{FD}takes on discrete values dependent on the gear number (1–16).

The traction controller takes the form of a PIDF controller. The limited range of continuous engine gains can be handled by a single PIDF controller with fixed gains. However, the range of the 16 different transmission gains *g*_{GR}*g*_{FD} is too large for a fixed-gain controller, and the controller is a gain scheduled based on the gear ratio *g*_{GR}.

This feed-forward term is the throttle input that is in effect just prior to entering the traction control mode and provides a stable transition between the leader–follower and traction control modes.

*K*

_{p}_{,}

*= 0.28,*

_{t}*K*

_{i}_{,}

*= 0.17,*

_{t}*K*

_{d}_{,}

*= 0.03, and the filter time constant is*

_{t}*t*= 0.3 s. These gains are divided by the gear ratio

_{f}*g*

_{GR}in each gear number to have a PIDF controller gain schedule for each gear. This way, the controller uses smaller gain values in gear numbers that have higher gear ratios and larger gain values in higher gear numbers with lower gear ratios. The aim here is that controller performance should be similar across different gear selections. This controller is discretized to operate the traction control loop in Fig. 6 at 20 Hz using the difference equation given by

The final result is 16 sets of controller coefficients *a*_{0}, *a*_{1}, and *a*_{2} and *b*_{0}, *b*_{1}, and *b*_{2} for all 16 gear numbers, which are provided in Ref. [28].

Gear selections in the traction control mode force the engine to operate within a bounded engine speed range of Ω* _{E}* ∈ [1350, 1950] so that the engine cannot stall or operate against the upper rpm limit (2100 rpm in gears 1–15 and 2300 rpm in gear 16). If the engine operating speed drifts below the lower bound of 1350 rpm, the transmission shifts down in gear number and if the engine speed goes above the upper bound of 1950 rpm, it shifts up in gear number. Gear shifts are limited to one change in gear number every 2 s, while detection of needed gear shifts occurs every 0.05 s so as to avoid limiting gear changes to every 2 s. For example, if a gear shift occurs at a time of

*t*=

*12*

*s*in the simulation, and the engine operates outside desired bounds at

*t*=

*14.5*

*s*but not

*t*=

*14*

*s*, a gear shift can occur before

*t*=

*16*

*s*. These changes in gear number are required for operation in the traction control mode, if there are changes in the tractor's ground speed so that the controller can track updated $\phi \u02d9ref$ values.

## Traction Control Mode Activation

The activation and deactivation of the traction control mode requires criteria for detecting soft terrain. There are a wide variety of possible approaches to this problem. Here, slip thresholds or limits are used to enter and exit the control mode. These thresholds are based on the bounds and discretization of the terrain parameter space in Sec. 4 and simulation experiments. Figure 7 shows a histogram of the peak slip *i _{pk}* for the maximum values of

*F*

_{net}for each model or hypothesis considered for estimation. The histogram shows that the upper and lower bounds for peak slip values are approximately

*i*∈ [5,40]. Since most of the terrain models have

_{pk}*i*≤ 25% traction control is initiated once $i\u0302>25%$. Higher and lower threshold values were evaluated. Higher values led to increased failure rates of the control mode since the controller turns on too late to prevent immobilization. Lower values risk the traction controller turning on when there is sufficient traction to pull the payload, and the controller may not be able to converge to $i\u0302pk$. The slip threshold used for exiting the traction control mode is $i\u0302<5%$. This is chosen since data collected in Antarctica in Fig. 3 show slip values below this threshold on firm terrain, and it is below almost all

_{pk}*i*values in Fig. 7.

_{pk}## Traction Control Mode Simulation

*t*=

*19.45 s, respectively, from open-loop throttle and gear commands. These are Π = [0.24, 0.26, 0.3, 0.4, 0.5, 0.6]*

^{T}and

*g*

_{GR}= [2, 3, 4, 5, 6, 7]

^{T}for the time vector

*t*= [0, 4, 8, 12, 16, 18]

^{T}. All tractors start on a firm terrain and travel into lower traction terrains stretching from 30 m to 130 m whose terrain parameters are in Table 1. Terrain parameters do not match any of the hypotheses of the 320 models from Sec. 4. Plots of slip

*i*versus

*F*

_{net,DB}for these lower traction terrains are provided in Fig. 8 where

*F*

_{net,DB}is defined as

This is done by rearranging the equations from Sec. 3 for the dynamics in phase 2. The zoomed in plot in Fig. 8(b) highlights the narrow domain of slip values where tractors must operate with positive *F*_{net,DB} values in order to maintain mobility.

A two-dimensional plot of lateral and longitudinal tractor positions taken at snapshots in time is shown in Fig. 9 for the five tractors colored magenta, cyan, green, blue, and red. All tractors are able to maintain tractor mobility using the traction control mode architecture except the red tractor, which operates on a terrain with the lowest cohesion *c* and shear deformation modulus *K* and becomes immobilized ∼53 m into the 100 m soft terrain stretch. Figure 10 shows the trajectories of all five tractors for the east or lateral position *X*, tractor speed *v _{T}*, engine speed in rpm Ω

*, selected gear ratio*

_{E}*g*

_{GR}, slip ratio

*i*, and throttle input Π. The cyan and green tractors maintain mobility across the soft terrain stretch and continuously increase their ground speeds and shift up in gear. The magenta and blue tractors also maintain tractor mobility but are only able to maintain and not increase their ground speeds. The red tractor, however, continuously shifts down in its gear selection to try and regulate to its $i\u0302pk$ slip ratio.

The performance of the traction control mode for each tractor depends on the accuracy of the peak slip point estimate $i\u0302pk$ from the Bayesian hypothesis selection. Figure 11 plots the estimated peak slip point $i\u0302pk$ and maximum net traction estimate $F\u0302net,max$ for each tractor. Solid colored lines correspond to true values and dots to estimated values. Since none of the terrain parameter vectors given in Table 1 match any of the hypotheses, the peak slip point *i _{pk}* cannot be perfectly estimated. In fact, each of the five terrain conditions differs in all of its terrain parameters from the 320 hypotheses except in the parameter

*n*. However, the cyan and green tractors $i\u0302pk$ estimates converge close to their true

*i*values of 24.5% and 21.5% so the traction control mode is able to maintain mobility and increase tractor ground speeds across the soft terrain patches. The $i\u0302pk$ estimates for the magenta and blue tractors converge to ∼24% and ∼22.5% when the true values for their terrains are 29% and 17%. Even though both of these estimates are off by 5%, the tractors to maintain mobility. The red tractor, however, converges to $i\u0302pk$ value of ∼21% when the true value is 8.5%. This error is much larger and leads to tractor immobilization.

_{pk}The value of *i _{pk}* for the red tractor's soft terrain is significantly different than the other soft terrains due to the very small shear deformation modulus

*K*=

*0.7 cm. This small*

*K*value makes it difficult to accurately estimate the correct

*i*value for two reasons. First, once the tractor operates under closed-loop traction control at a higher slip, determining whether or not the underlying terrain has a high or low

_{pk}*K*value is difficult since this parameter implies a steep slope of force-slip curves at lower slip values.

## Winch Control Mode

Section 7 presented the architecture of the traction control mode and demonstrated its effectiveness in simulation results. However, if there is no domain of slip ratios that provide positive *F*_{net,DB} values, this mode will not prevent the tractor from being immobilized and the drawbar load on the vehicle must be reduced, which requires a towing winch to release the payload. An approach to automate winch use is proposed in this section and is referred to as the winch control mode. The hybrid rigid body dynamics and hydraulic actuator modeling for simulating this mode was discussed in Sec. 3.

The winch control mode uses a low bandwidth 4 Hz logic-based controller. This allows winch torques to be varied iteratively based on the vehicle's response to drawbar loading conditions. The low bandwidth of the controller provides the traction control mode time to respond to new loading conditions when the two are active simultaneously.

There are two stages during winch use. The first stage starts by activating the winch control mode. At this point, the tractor detects that the traction control mode alone cannot maintain vehicle mobility when $i\u0302\u226545%$ and lets out the payload. Subsequently, the hydraulic winch slowly increases its hydraulic brake torque so as to bring the sled back into motion and reduce the rate at which the cable is being let out until it reaches zero. This is done by slowly increasing the hydraulic pressure of the brake valve. However, if the brake torque is increased so as to overload the drawbar, the hydraulic pressure is reduced, and then, the pressure is slowly increased again at a reduced rate. If the tractor reaches firm enough terrain, this will eventually result in the sled payload being brought back into motion and the winch will stop once the tractor and sleds speed match. The second stage of winch use pulls the payload back in. If stage 1 reaches completion, the terrain has sufficient shear strength to fully recover the payload and stage 2 will recover the payload more aggressively. The specifics of each stage are discussed in this section.

### Stage 1: Winch Activation and Initializing Payload Recovery.

Since recovering lost cable length can be difficult and time consuming to recover, there should be a strong indication that the traction control architecture alone cannot maintain tractor mobility before using the winch. At the same time, waiting too long to activate the winch could result in an immobilized tractor that has excavated itself into the terrain. To strike a balance between both concerns, the bounds on the terrain parameter space and their resulting slip–force curves are used as a guide. Figure 7 shows a histogram of the peak slip points *i _{pk}* for all terrain hypothesis models. The largest value of

*i*in Fig. 7 is

_{pk}*i*≃ 41%. Therefore, once the estimates of $i\u0302\u226545%$, this indicates that traction control alone cannot maintain tractor mobility, and the winch is activated to let out payload and reduce drawbar load.

_{pk}The starting position of the four-way two-position valve in the winch hydraulics is shown in Fig. 5(b) before winch activation. This position is maintained once the control mode has been activated and the *P*_{set} value transitions from *P*_{set} = 2700 psi to *P*_{set} = 0 psi in Eq. (22) for the braking valve B2. This allows more hydraulic fluid to flow out of the *P _{O}* node to reduce fluid pressure and subsequently the braking torque

*τ*.

_{W}Once the winch has been activated, attempts to recover the payload begin immediately. This is due to the fact that soft terrain patch lengths are unknown to the tractor. Furthermore, recovering all the cable let out once the winch has been activated can be a time-consuming process.

Algorithm 1 runs at 4 Hz and summarizes the logic used to initialize payload recovery once the winch control mode has been activated. Line 4 gathers the latest original, unsmoothed slip estimate $i\u0302$ from the DTKF. Once this has been obtained, the hydraulic braking valve regulates the pressure set point *P*_{set} value based on current and past slip values. When the tractor first enters the control mode, *P*_{set} = 0 psi and increments at levels of incPSI = 15 psi in line 6 until the track slip reaches 20%. This threshold is chosen so that if the current drawbar load is at its limit and slip increases, the traction control mode has a chance to maintain mobility without operating under continuous load increases. The limiting factor of performance here is that tractors must initialize payload recovery on terrains having positive *F*_{net} with *i *≤* *20% as is the case with all curves in Fig. 8 except the magenta curve. Otherwise the algorithm cannot recover the payload. It should be noted, however, that the magenta curve has a *K* on the edge of the bounded terrain parameter space. If the pressure reference *P*_{set} that has been maintained in line 8 causes the tractor to detect $i\u0302>35%$, the *P*_{set} value is significantly reduced until the drawbar load is reduced to a point where $i\u0302<35%$ in line 10. Each time *P*_{set} is significantly reduced, incPSI is reduced at a rate determined by *aCoeff* where incPSI ≥ 5 psi. Each time this occurs, the tractor is allowed to make another attempt to pull in the load without having to reset *P*_{set} = 0 psi and can increase *P*_{set} to reload the drawbar at a slower rate so that the tractor operates at the appropriate drawbar load for the terrain. While Algorithm 1 runs at 4 Hz, $i\u0302$ is calculated at 20 Hz to make sure the tractor is not operating above 45% track slip. If this occurs, the braking pressure is reset to *P*_{set} = 0

Stage 1 is considered complete once the winch has stopped rotary motion and is no longer releasing cable length. At this point, the tractor is ready to enter stage 2 of the winch control mode to make a full recovery of the payload.

### Stage 2: Full Payload Recovery.

Stage 1 uses the braking valve B2 to reduce the rate of cable loss $\psi \u02d9rW$ until $\psi \u02d9rW=0$. At this point, stage 2 is initiated to make a full recovery of the payload and cable. When stage 2 is initiated, the four-way two-position valve is put into position I. Figure 5(a) illustrates the valve position change from II to I. This allows the variable displacement, pressure compensated pump to push hydraulic fluid in the H node, increase the pressure, and apply torque to the fixed displacement hydraulic motor to pull in the payload. Once this occurs, the *P*_{set} value in Eq. 30 is set to *P*_{set} = 15 psi [20].

This small initial *P*_{set} value is used since *P _{H}* = 0 psi at the instant the valve position change occurs and the variable displacement, pressure compensated pump model puts the pump to its maximum displacement value

*D*

_{P}_{,Max}if

*P*≤

_{H}*P*

_{set}. This minimizes the amount of time the pump fully loads the engine to provide a stable transition and to avoid stalling out the tractor's engine. After the initial transition, the winch follows a similar procedure to Algorithm 1. This brief rule-based procedure is summarized in Algorithm 2 and is identical to the procedure used in stage 1 but omits the cautionary step of greatly reducing drawbar pull if the vehicle detects $i\u0302\u226535%$. There are several reasons for this. First, if the tractor has progressed from stages 1 to 2, the immediate terrain should provide enough traction to reel in the payload. Second, if the tractor is not able to recover the entire length of the cable due to another immediate soft terrain encounter, it is likely that the soft terrain patches that occur in the operating environment are too frequent and/or too long. Third, the drawbar requirements in the second stage are slightly different than in the first stage. In the first stage, the winch is attempting to bring the speed of the winch $\psi \u02d9rW$ to zero. To make progress, the second derivative $\psi \xa8rW$ must always be negative. This acceleration requires a drawbar load constantly in excess of what is required to pull the payload in order to complete the first stage. However, in stage 2, a constant winch speed results in progress toward full recovery of the payload. Therefore, additional drawbar load in excess of what is required to tow the payload or a nonzero value of $\psi \xa8rW$ is only needed briefly instead of constantly as in stage 1.

## Winch Control Mode Simulation

In Sec. 7, results showed that all five tractors except the red one were able to maintain mobility across 100 m soft terrain stretch. However, the red tractor did traverse ∼53 of 100 m. Now, the red tractor will use its winch to traverse the remaining ∼47 m and beyond. Furthermore, it should be noted that the combined capability of the traction control mode and the winch control mode allows for terrain traversal without having to forfeit cable across the entire 100 m to minimize the amount of cable length required.

Figures 12–14 show trajectory data of the red tractor using both the traction control mode and winch control mode. Figure 12 shows a two-dimensional plot of the tractor position at snapshots in time. The terrain operating conditions for this red tractor are identical to those in Sec. 7 and can be referenced in Table 1 and Fig. 8. Here, the tractor operates as before until ∼60 s into the simulation when the winch control mode becomes activated in stage 1 as triggered by a slip estimate $i\u0302\u226545%$. This is evident in Figs. 13 and 14 where the slip *i*, pressure in the hydraulics *P _{O}*, and drawbar load DB all drop to values close to zero to let the payload out at ∼60 s. Subsequently,

*P*

_{set}is slowly incremented as outlined in Algorithm 1. This can be seen in the slow increase in

*P*(Fig. 13) and in the drawbar load DB (Fig. 14). At ∼83 s, the tractor has detected that the drawbar is overloaded and decreases

_{O}*P*

_{set}and incPSI in lines 10 and 11 of Algorithm 1. This rapidly decreases the load on the tractor and increases the drawbar load more slowly for another attempt to initialize the payload recovery. At ∼90 s, the tractor exits the soft terrain stretch and is able to bring the winch to a locked state where no more cable is being let out. This can be seen in Fig. 14 where

*ψr*is the amount of cable length let out and $\psi \u02d9rW$ is the speed at which cable is being let out. This process takes ∼15 s and stage 2 is entered ∼105 s into the simulation. The initial kickback from switching the discrete four-way two-position valve to position I creates a spike in the hydraulic pressure

_{W}*P*which highlights the reason for initializing the

_{H}*P*

_{set}value to a small value of 15 psi upon transition into stage 2. The

*P*

_{set}value is then slowly increased as outlined in Algorithm 2 to slowly increase the load on the engine

*τ*

_{load}and the speed at which the payload is pulled in $\psi \u02d9rW$.

## Conclusion

A comprehensive tracked vehicle dynamics model equipped with a towing winch has been developed to include the effects of the power-train, terrain–vehicle interaction, and payload on deformable terrain. This allows for an iterative, model-based design approach for a traction and winch control modes. The traction control mode implements a complete control architecture of terrain identification and slip setpoint control based on realistic sensor data, actuator limitations, and a validated model. Simulation results validate the controller in different terrain conditions. Since the traction control mode cannot perfectly estimate terrain conditions due to noisy data and the nonlinear nature of slip traction curves, a logic-based feedback controller is proposed for the winch so that vehicles can extend the domain under which mobility is retained. Simulation results showed the effectiveness of the controller in allowing a tractor to maintain mobility by reducing drawbar load, where it previously could not use only the traction control mode.

Future work will investigate more sophisticated approaches for detecting transitions between the activation and deactivation of control modes. In this work, threshold values were used based on the terrain parameter space. This would allow for a more general applicability of the control architectures presented here by eliminating the need for tuned, heuristic thresholds.

## Acknowledgment

The authors thank Jeffery Zimmerman, Darin Motz, Ravi Godbole, and others from AGCO corporation who provided intellectual support in accessing tractor MT865 CAN bus data in Antarctica. Furthermore, the collection of these data was also made possible through coordination efforts by South Pole Traverse manager Kory Hotel and tractor operators Brad Johnson and Robert Shaw who drove tractors during experiments.