Analytical Study on the Cornering Behavior of an Articulated Tracked Vehicle †

: Articulated tracked vehicles have been traditionally studied and appreciated for the extreme maneuverability and mobility ﬂexibility in terms of grade and side slope capabilities. The articulation joint represents an attractive and advantageous solution, if compared to the traditional skid steering operation, by avoiding any trust adjustment between the outside and inside tracks. This paper focuses on the analysis and control of an articulated tracked vehicle characterized by two units connected through a mechanical multiaxial joint that is hydraulically actuated to allow the articulated steering operation. A realistic eight degrees of freedom mathematical model is introduced to include the main nonlinearities involved in the articulated steering behavior. A linearized vehicle model is further proposed to analytically characterize the cornering steady-state and transient behaviors for small lateral accelerations. Finally, a hitch angle controller is designed by proposing a torque-based and a speed-based Proportional Integral Derivative (PID) logics. The controller is also veriﬁed by simulating maneuvers typically adopted for handling analysis.


Introduction
The cornering behavior of tracked vehicles has peculiar characteristics that are different from wheeled vehicles. The steering operation of these vehicles can be accomplished through different mechanisms: skid-steering and steering by articulation. In skid steering, a turning yaw moment is obtained by applying a different longitudinal thrust force between left and right track sides [1,2]. For tracked vehicles with two or more units, the steering operation may be achieved by a relative yaw rotation between units through a specific mechanism on the connecting joint [3,4]. This solution is also preferred to the skid steering since it does not require a thrusts adjustment between the outside and inside tracks so that the resultant forward thrust can be maintained during a turning maneuver, as shown in [5]. Nowadays, articulated tracked vehicles have been widely adopted in several engineering fields, e.g., planetary exploration, military [6], agriculture [7] and construction [8].
The articulated steering is usually accomplished by installing a hydraulic steering system able to provide the power required to overcome vehicle lateral resistances arising during the relative rotation between the two units [8]. The kineto-dynamic behavior of the hydraulic steering mechanism on articulated vehicles was analyzed in [9,10], thus showing a correlation between the hydraulic steering dynamics with the oscillatory articulation response and weaving motion of the vehicle. Moreover, the increasing research focus on autonomous driving has encouraged the development of an automatic steering controller for path tracking or path following applications. Conventional hydraulic steering systems can be replaced or adapted to allow an automatic control: in [11][12][13] a pressure following control is obtained by exploiting the current onboard vehicle hydraulic systems, meanwhile an automatic steering was proposed by [14] through the design of a multifunctional hydraulic steering circuit; an autonomous guidance strategy was also presented in [15] for a load haul dump vehicle where the relationship between vehicle stability and speed is obtained to improve the vehicle dynamic response.
However, these research activities are limited to articulated wheeled vehicles and they cannot be directly extended to articulated tracked vehicles (ATVs) due to the continuous track-terrain contact distribution. Indeed, only few similar studies are available in literature about tracked articulated vehicles. A mathematical model for plane motion to predict the steerability and mobility of an ATV was developed in [5] through numerical simulations. More recently, an improved ATV dynamics model was proposed in [16] by considering the shear stress-shear displacement relation of the soil at the track-terrain interface. The soil deformation on track-soil interaction was also developed in [17], where the side bulldozing effect was also included for improving the model accuracy in simulating the ATV steering operation. A non-linear mathematical model of an ATV was also described by [18,19] where the effect of a hitch angle, i.e., the relative yaw angle between the ATV units, feedback controller was analyzed through steady-state cornering maneuvers. The control of the ATV hitch angle was also adapted for autonomous applications as described in [20], where a path tracking control was designed based on the distance deviation and the heading angle deviation between the ATV and the desired path.
To the best knowledge of the authors, there is a general lack of comprehensive analytical investigation of the ATV transient and steady-state behaviors. Most of the results available in literature are numerically obtained through detailed mathematical models but without any analytical correlation to ATV parameters and operative conditions, e.g., vehicle speed. The activity presented in this paper aims at bridging the gap by providing an analytical approach for analyzing and controlling the ATV cornering response. The paper content represents an extension of the preliminary work introduced in [19] by providing further details related to the non-linear ATV model and by introducing a simplified linearized model, valid for small lateral accelerations. Moreover, the concept of the understeer characteristics, typical of wheeled passenger cars, is also extended to ATV categories. Although the main paper contribution focuses on the methodology for an analytical evaluation of the ATV lateral dynamics behavior, a hitch angle controller is presented to explore the influence of a hydraulic actuation system for the steering operation. The hitch angle controller is designed and calibrated on the non-linear model through an optimization procedure.
The paper structure is organized as follows: a non-linear mathematical model of the ATV is presented in Section 2 by including the effect of rigid bodies planar dynamics, hydraulic steering dynamics, and angular tracks dynamics. The ATV model is then linearized and an analytical solution for the transient and steady-state responses is provided for small lateral accelerations in Section 3. A hitch angle controller is then designed in Section 4 to achieve the desired performance. The controller is then validated through simulations with the non-linear model. Finally, some conclusions are drawn in Section 5.

ATV Non-Linear Mathematical Model
A simplified schematic of the ATV considered within the rest of the paper is reported in Figure 1.
The vehicle architecture is characterized by two units connected through a multi degrees of freedom steering joint. The front unit hosts the driveline components so that the second unit may fulfill the payload functionalities of goods and personnel transportation. The weight of each ATV unit is supported by four or five couples of road wheels. Each road wheel is connected to the hull through a torsion bar suspension system that, together with the set of bushing, represents the elastic-damping components of the ATV. The road wheels are driven by the two sprocket wheels of each unit through a track mechanism. The engine torque is assumed to be equally distributed between the two units and further split in equal parts between the two sprocket wheels of each unit. The ATV steering operation Machines 2021, 9,38 3 of 23 is entrusted to a hydraulic system mounted on the steering joint that provides a steering torque to each unit through the regulation of a flow rate proportional valve.
Machines 2021, 9, x FOR PEER REVIEW 3 of 23 wheels are driven by the two sprocket wheels of each unit through a track mechanism. The engine torque is assumed to be equally distributed between the two units and further split in equal parts between the two sprocket wheels of each unit. The ATV steering operation is entrusted to a hydraulic system mounted on the steering joint that provides a steering torque to each unit through the regulation of a flow rate proportional valve. The ATV dynamics was analyzed through a non-linear mathematical model with eight degrees of freedom (8-DOF) by including the longitudinal, lateral and yaw motion of the front unit, the yaw motion of the rear unit and the angular motion of the four sprocket wheels. The main assumptions for the ATV non-linear model are: • the two units are considered as rigid bodies with mass and and moment of inertia around the vertical axis through the center of gravity (CoG) and ; • the road is considered rigid and flat (the effect of terrain slopes or sinkage is neglected); • the continuous track-terrain contact force distribution is discretized with four contact points for each track; • the steering joint is considered as an ideal yaw rotational hinge placed at a distance of from the front CoG and from the rear CoG; • the front and rear tracks width is equal to ; • the number of road wheels is = 4 for each track; • flexible deformation of the tracks is not considered.

Rigid Bodies Dynamics
The free body diagram for the ATV planar motion modelling is shown in Figure 2 where a moving reference coordinate system ( , , ) is centered in the front unit CoG and a second moving reference coordinate system ( , , ) is placed in the rear unit CoG. Reprinted with permission from ref. [19]. Copyright 2020 Springer Nature Switzerland AG. The ATV dynamics was analyzed through a non-linear mathematical model with eight degrees of freedom (8-DOF) by including the longitudinal, lateral and yaw motion of the front unit, the yaw motion of the rear unit and the angular motion of the four sprocket wheels. The main assumptions for the ATV non-linear model are: • the two units are considered as rigid bodies with mass m 1 and m 2 and moment of inertia around the vertical axis through the center of gravity (CoG) I G1 and I G2 ; • the road is considered rigid and flat (the effect of terrain slopes or sinkage is neglected); • the continuous track-terrain contact force distribution is discretized with four contact points for each track; • the steering joint is considered as an ideal yaw rotational hinge placed at a distance of x G1 from the front CoG and x G2 from the rear CoG; • the front and rear tracks width is equal to T; • the number of road wheels is N w = 4 for each track; • flexible deformation of the tracks is not considered.

Rigid Bodies Dynamics
The free body diagram for the ATV planar motion modelling is shown in Figure 2 where a moving reference coordinate system R 1 (x 1 , y 1 , z 1 ) is centered in the front unit CoG and a second moving reference coordinate system R 2 (x 2 , y 2 , z 2 ) is placed in the rear unit CoG.
Machines 2021, 9, x FOR PEER REVIEW 3 of 23 wheels are driven by the two sprocket wheels of each unit through a track mechanism. The engine torque is assumed to be equally distributed between the two units and further split in equal parts between the two sprocket wheels of each unit. The ATV steering operation is entrusted to a hydraulic system mounted on the steering joint that provides a steering torque to each unit through the regulation of a flow rate proportional valve. The ATV dynamics was analyzed through a non-linear mathematical model with eight degrees of freedom (8-DOF) by including the longitudinal, lateral and yaw motion of the front unit, the yaw motion of the rear unit and the angular motion of the four sprocket wheels. The main assumptions for the ATV non-linear model are: • the two units are considered as rigid bodies with mass and and moment of inertia around the vertical axis through the center of gravity (CoG) and ; • the road is considered rigid and flat (the effect of terrain slopes or sinkage is neglected); • the continuous track-terrain contact force distribution is discretized with four contact points for each track; • the steering joint is considered as an ideal yaw rotational hinge placed at a distance of from the front CoG and from the rear CoG; • the front and rear tracks width is equal to ; • the number of road wheels is = 4 for each track; • flexible deformation of the tracks is not considered.

Rigid Bodies Dynamics
The free body diagram for the ATV planar motion modelling is shown in Figure 2 where a moving reference coordinate system ( , , ) is centered in the front unit CoG and a second moving reference coordinate system ( , , ) is placed in the rear unit CoG. Reprinted with permission from ref. [19]. Copyright 2020 Springer Nature Switzerland AG. Figure 2. Free body diagram of ATV (top view). Reprinted with permission from ref. [19]. Copyright 2020 Springer Nature Switzerland AG.
The ATV longitudinal and lateral equilibrium equations are described by (see also [18]): ..
where a x1 and a y1 are the front unit longitudinal and lateral acceleration components; a x2 and a y2 are the rear unit longitudinal and lateral acceleration components; u 1 and v 1 are the front unit longitudinal and lateral velocity components; u 2 and v 2 are the rear unit longitudinal and lateral velocity components. Front and rear velocity vectors, V 1 and V 2 , are inclined by an angle of β 1 and β 2 with respect to the corresponding longitudinal direction. ψ 1 and ψ 2 are the yaw angles of the front and rear units, respectively. α = ψ 1 − ψ 2 represents the hitch angle between the two ATV units. F I x,ji and F I I x,ji represent the longitudinal force components between the j th = L, R track and the terrain in the i th road wheel contact point for the front and rear units, respectively. F I y,ji and F I I y,ji represent the lateral force components between the j th = L, R track and the terrain in the i th road wheel contact point for the front and rear units, respectively; F aero is the longitudinal component of the aerodynamics force.
The front and rear yaw moment balance equations are described by: .. .. ..
where M I z and M I I z are the front and rear yaw moments generated by the track-terrain forces for the front and rear units, respectively. d I i and d I I i are the distances between the i th axle to the front and rear CoG, respectively, and they are assumed positive if placed in front of their respective CoG; C s is the steering torque applied by the hydraulic system.

Track Angular Dynamics
Each track rotational dynamics is described by the free body diagram in Figure 3.  The rotational equilibrium is expressed by: where is the angular speed of the sprocket wheel and = / is the road wheel angular speed on the left/right ( = / ) side of the = , unit. and are the sprocket wheels and road wheels radius, respectively.
is the engine torque, is the overall transmission efficiency, is the transmission ratio between engine shaft and each sprocket wheel. , is the dynamic vertical force on each road wheel.
is the equivalent moment of inertia around each sprocket wheel axis and includes all main inertial contributions of the whole powertrain system. Each track rolling resistance , is described with a hyperbolic tangent function by introducing the constant coefficient and the track speed ( = ) dependent coefficient . Finally, is a pre-defined threshold for the hyperbolic tangent shape factor.

Track-Terrain Contact Model
The track-terrain contact distribution is modeled by an equivalent finite number of contact patches equal to the number of road wheels [21]. Both contact patch force components , and , are modeled as hyperbolic tangent functions, depending on their longitudinal slip ratios and slip angles : where the coefficients , and , are introduced to consider the saturation of the available adhesion to the vertical load. The peak longitudinal and lateral forces occur at and at , respectively. and coefficients enable the longitudinal/lateral combined slips influence.
Each road wheel slip ratio and slip angle is defined as: The rotational equilibrium is expressed by: where ω k sj is the angular speed of the sprocket wheel and ω k W j = ω k sj R s /R W is the road wheel angular speed on the left/right (j = L/R) side of the k th = I, I I unit. R s and R W are the sprocket wheels and road wheels radius, respectively. C M is the engine torque, η is the overall transmission efficiency, τ is the transmission ratio between engine shaft and each sprocket wheel. F k z,ji is the dynamic vertical force on each road wheel. I s is the equivalent moment of inertia around each sprocket wheel axis and includes all main inertial contributions of the whole powertrain system. Each track rolling resistance C k Roll,j is described with a hyperbolic tangent function by introducing the constant coefficient f 0 and the track speed (V k j = ω k sj R s ) dependent coefficient f 2 . Finally, ω th is a pre-defined threshold for the hyperbolic tangent shape factor.

Track-Terrain Contact Model
The track-terrain contact distribution is modeled by an equivalent finite number of contact patches equal to the number of road wheels [21]. Both contact patch force components F k x,ji and F k y,ji are modeled as hyperbolic tangent functions, depending on their longitudinal slip ratios σ k ji and slip angles β k ji : where the coefficients a 1x , a 2x and a 1y , a 2y are introduced to consider the saturation of the available adhesion to the vertical load. The peak longitudinal and lateral forces occur at σ max and at β max , respectively. k β and k σ coefficients enable the longitudinal/lateral combined slips influence. Each road wheel slip ratio and slip angle is defined as: where u k ji and v k ji are, respectively, the longitudinal and lateral components of the i th road wheel in the j th = L/R side of the k th unit.

Dynamic Vertical Forces
The vertical forces on each road wheel are calculated as the sum of the static and the dynamic load contributions. The static load contribution F k zST,ji is evaluated by considering the mass distribution among the four axles of each unit. Due to the statically indeterminate nature of the problem, i.e., cantilever beam supported by a redundant number of constraints, the mass distribution is correlated to the suspension system between the sprung (ATV hull) and the unsprung (track mechanisms) masses. Indeed, each road wheel is connected to the ATV hull by a torsion bar, as shown in Figure 4.
where and are, respectively, the longitudinal and lateral components of the road wheel in the = / side of the unit.

Dynamic Vertical Forces
The vertical forces on each road wheel are calculated as the sum of the static and the dynamic load contributions. The static load contribution , is evaluated by considering the mass distribution among the four axles of each unit. Due to the statically indeterminate nature of the problem, i.e., cantilever beam supported by a redundant number of constraints, the mass distribution is correlated to the suspension system between the sprung (ATV hull) and the unsprung (track mechanisms) masses. Indeed, each road wheel is connected to the ATV hull by a torsion bar, as shown in Figure 4. Each torsion bar reacts to the relative motion between the road wheel and the ATV hull, thus behaving as an equivalent vertical spring. The following assumptions are considered for the ATV suspension system:

•
Each suspension torsion bar is modeled with a linear characteristic of the equivalent spring deflection.

•
Each equivalent spring has no static preload.

•
The static deflection is equal for each spring of the same unit (no static pitch and/or roll).

•
The corresponding right and left springs of each unit's axle share the same stiffness value.
Under these assumptions, the static load contribution is then calculated as: where and are the equivalent vertical spring stiffness for the front and rear units, respectively; and represent the mass distribution on the axle for the Each torsion bar reacts to the relative motion between the road wheel and the ATV hull, thus behaving as an equivalent vertical spring. The following assumptions are considered for the ATV suspension system:

•
Each suspension torsion bar is modeled with a linear characteristic of the equivalent spring deflection.

•
Each equivalent spring has no static preload.

•
The static deflection is equal for each spring of the same unit (no static pitch and/or roll).

•
The corresponding right and left springs of each unit's axle share the same stiffness value.
Under these assumptions, the static load contribution is then calculated as: where K I Si and K I I Si are the i th equivalent vertical spring stiffness for the front and rear units, respectively; q I i and q I I i represent the mass distribution on the i th axle for the front and rear units, respectively. It is remarkable to note that the assumptions made for the ATV suspension system imply that the static mass distribution corresponds to the spring stiffness distribution among each unit's axle. During cornering maneuvers, the lateral acceleration implies a load transfer between the inside and outside tracks for the front unit, ∆F I z = m 1 a y1 h G1 /T, and for the rear unit, ∆F I I z = m 2 a y2 h G2 /T. Even if the roll dynamics is not included in the model, the steadystate roll angle represents an internal variable required to evaluate the spring deflection on each unit's track. The steady-state roll angle is calculated through the moment balance equation of each unit sprung mass around the correspondent roll axis, assumed in a fixed position from the ground.
By assuming small roll angles, i.e., neglecting the distance between CoG and CoG' in Figure 5, the moment balance equation of each unit sprung mass around its respective roll axis is expressed by: where φ 1 and φ 2 are the roll angles for the front and rear units, respectively; h roll1 and h roll2 are the roll axis height for the front and rear units, respectively; K I roll i = K I Si T 2 /2 and K I I roll i = K I I Si T 2 /2 represent the equivalent roll stiffness produced by the two side springs of each i th axle for the front and rear units, respectively.
Machines 2021, 9, x FOR PEER REVIEW 7 of 23 front and rear units, respectively. It is remarkable to note that the assumptions made for the ATV suspension system imply that the static mass distribution corresponds to the spring stiffness distribution among each unit's axle. During cornering maneuvers, the lateral acceleration implies a load transfer between the inside and outside tracks for the front unit, Δ = ℎ / , and for the rear unit, Δ = ℎ / . Even if the roll dynamics is not included in the model, the steadystate roll angle represents an internal variable required to evaluate the spring deflection on each unit's track. The steady-state roll angle is calculated through the moment balance equation of each unit sprung mass around the correspondent roll axis, assumed in a fixed position from the ground. By assuming small roll angles, i.e., neglecting the distance between CoG and CoG' in Figure 5, the moment balance equation of each unit sprung mass around its respective roll axis is expressed by: where and are the roll angles for the front and rear units, respectively; ℎ and ℎ are the roll axis height for the front and rear units, respectively; = /2 and = /2 represent the equivalent roll stiffness produced by the two side springs of each axle for the front and rear units, respectively. The moment balance equation of each unsprung mass around the roll axis provides the lateral load transfer on each axle for the front (Δ ) and rear (Δ ) units, as shown in the free body diagram of Figure 5: By assuming ℎ ≪ ℎ and ℎ ≪ ℎ the load transfer on each axle for the front and rear units is finally defined as: The moment balance equation of each unsprung mass around the roll axis provides the lateral load transfer on each axle for the front (∆F I zi ) and rear (∆F I I zi ) units, as shown in the free body diagram of Figure 5: By assuming h roll1 h G1 and h roll2 h G2 the load transfer on each axle for the front and rear units is finally defined as: Equation (11) states the load transfer due to lateral acceleration is distributed as the mass distribution among each unit's axles. Finally, the total vertical force on each track road wheel for the front and rear units is expressed by: where the sign of the lateral load transfers ∆F I zi and ∆F I I zi is referred to the left (j = L) or to the right (j = R) side of the front and rear units, respectively, as shown in Figure 5.

Hydraulic Steering System
An example of the hydraulic steering circuit is represented in Figure 6.
Equation (11) states the load transfer due to lateral acceleration is distributed as the mass distribution among each unit's axles. Finally, the total vertical force on each track road wheel for the front and rear units is expressed by: where the sign of the lateral load transfers Δ and Δ is referred to the left ( = ) or to the right ( = ) side of the front and rear units, respectively, as shown in Figure 5.

Hydraulic Steering System
An example of the hydraulic steering circuit is represented in Figure 6. The steering torque is generated by the two hydraulic piston forces as results of pressure and by: where and are, respectively, the minimum distance of the right and the left piston rods from the hinge (see Appendix A for their analytical expression). and represent the right and left forces, respectively, applied by the corresponding piston rod to ATV bodies: Two hydraulic cylinders are linked to the front unit through C and D hinges and their piston rods are linked to the rear unit through A and B hinges. The fluid pressure p 1 and p 2 between the cylinder chambers depends on flow rates Q 1 and Q 2 controlled by a flow rate proportional valve VP 2 through the regulation of the X v command. A constant displacement pump provides the flow rate Q p to cope with the VP 2 valve request, Q VP2,in . The excessive flow rate, Q VP1 , is absorbed by a second flow rate proportional valve VP 1 that receives the complimentary command 1 − X v . This solution also prevents the pump delivery pressure p P from growing excessively. For safety reasons, all pressures are saturated by pressure relief valves.
The steering torque C s is generated by the two hydraulic piston forces as results of pressure p 1 and p 2 by: where f R and f L are, respectively, the minimum distance of the right and the left piston rods from the hinge O (see Appendix A for their analytical expression). F R and F L represent the right and left forces, respectively, applied by the corresponding piston rod to ATV bodies: where A p is the cylinder bore area and A s is the piston rod area. c oil is the damping coefficient due to the hydraulic fluid. v R and v L are the right and left piston rod extension/compression speeds, respectively (see Appendix A for their analytical expression). Pressure p 1 , p 2 and p P dynamics are defined by the following expressions: where β oil is the fluid bulk modulus. V 1 and V 2 are the equivalent cylinder chambers corresponding to p 1 , p 2 , respectively. Their variation is kinematically imposed on the ATV bodies' dynamics (see Appendix A for their analytical expression). V pipes is the equivalent pipes volume used for the pump pressure, p P , dynamics. Q 1 and Q 2 are the flow rates delivered to chambers V 1 and V 2 .
in and Q VP2.out are the flow rates through the valve VP 2 pump-connected and tank-connected ports, respectively. Q VP1 is the flow rate delivered through the VP 1 valve. Q VP2.in , Q VP2.out and Q VP1 are calculated by: where c q is the coefficient of discharge; A max is VP 1 and VP 2 maximum opening valve area. VP 2 pressure drops are ∆p VP2,in = p P − p 1 and ∆p VP2,out = p 2 − p atm when X v < 0, otherwise ∆p VP2,in = p P − p 2 and ∆p VP2,out = p 1 − p atm . VP 1 pressure drop is ∆p VP1 = p P − p atm .

Stable Equilibrium Points
The ATV non-linear model was implemented in MATLAB® and Simulink® environments. The ATV equilibrium conditions were then obtained and represented as characteristic plots between the desired quantities versus the hitch angle, as shown in Figure 7.
Each subplot marker in Figure 7 (square markers refer to the front unit quantities and circle markers to the rear unit) represents an ATV equilibrium condition. Sideslip angles, yaw rates, curvatures, lateral accelerations, and yaw moments characteristics show a linear dependence at least for a limited range of the hitch angle. On the other hand, the relation between the steering torque C s and the hitch angle α is strongly non-linear. The presence of a peak, occurring at 10 deg for u 1 = 20 km/h, represents a critical point for the relation between the hitch angle and the steering torque: before the peak, the steering torque monotonically increases with the hitch angle, meanwhile, after that peak, an increment of the hitch angle will cause a reduction of the steering torque. This trend leads to an inversion of the steering torque sign that can cause an ATV jackknifing effect for higher hitch angles (see also [18]). Based on these simulation results, a linearized model of the ATV in different equilibrium points may simplify the analysis of the vehicle lateral dynamics.

Linearized ATV Lateral Dynamics Model
This section aims to present a linearized version of the ATV non-linear model for lateral dynamics analysis. Since the paper focus is on lateral dynamics, the following hypothesis are considered within the rest of the paper: The effect of the hydraulic steering dynamics, introduced for the non-linear model in Equation (15), was not included in the linearized ATV model: an ideal actuator is assumed for the steering operation so that the actuator applies a desired steering torque with an infinite torque bandwidth that is not influenced by the internal feedback between the actuator and system dynamics. This represents a strong hypothesis since the actuation system may influence the overall vehicle dynamics, but the main purpose of the ATV linearized model is to analyze only the ATV rigid bodies dynamics. The effect of the hydraulic steering dynamics is partially assessed through numerical simulations with the non-linear model in Section 4, where the hydraulic pressure dynamics influences the design of the hitch angle controller.
Equations (2)- (4) show that it is convenient to introduce the axle forces instead of single road wheel contact forces. Given the ATV geometrical symmetry, the left and right road wheels of each axle have almost the same slip angle value ( = ≅ ) so that the ℎ axle lateral force for each ℎ unit is calculated as follows: , ( , ) = , + , = − [ 1 ( , 2 + , 2 ) + 2 ( , + , )] tanh 3 where is the equivalent slip angle of the ℎ axle for the ℎ = , unit:

Linearized ATV Lateral Dynamics Model
This section aims to present a linearized version of the ATV non-linear model for lateral dynamics analysis. Since the paper focus is on lateral dynamics, the following hypothesis are considered within the rest of the paper:
Neglected longitudinal acceleration (a x1 = 0) and rear unit longitudinal force The longitudinal speed of the front unit u 1 is considered a time-independent quantity: u 1 is kept constant during each simulation/maneuver.
The effect of the hydraulic steering dynamics, introduced for the non-linear model in Equation (15), was not included in the linearized ATV model: an ideal actuator is assumed for the steering operation so that the actuator applies a desired steering torque with an infinite torque bandwidth that is not influenced by the internal feedback between the actuator and system dynamics. This represents a strong hypothesis since the actuation system may influence the overall vehicle dynamics, but the main purpose of the ATV linearized model is to analyze only the ATV rigid bodies dynamics. The effect of the hydraulic steering dynamics is partially assessed through numerical simulations with the non-linear model in Section 4, where the hydraulic pressure dynamics influences the design of the hitch angle controller.
Equations (2)- (4) show that it is convenient to introduce the axle forces instead of single road wheel contact forces. Given the ATV geometrical symmetry, the left and right road wheels of each axle have almost the same slip angle value (β k Li = β k Ri ∼ = β k i ) so that the i th axle lateral force for each k th unit is calculated as follows: where β k i is the equivalent slip angle of the i th axle for the k th = I, I I unit: Each axle characteristics F k y,i a y , β k i , together with their cornering stiffness, are reported in Figure 8 for the front unit (left subplot) and the rear unit (right subplot). Each axle characteristics , ( , ), together with their cornering stiffness, are reported in Figure 8 for the front unit (left subplot) and the rear unit (right subplot). The axle force characteristics shows a linear behavior for small slip angles, which tends to saturate for higher values. Based on that, each axle lateral force in Equation (17) is linearized around the linearization point identified by 0 : where , 0 , 0 and represent the value of the ℎ axle lateral force, slip angle and cornering stiffness, respectively, calculated at the linearization point. The cornering stiffness is expressed as: where , 0 and , 0 represent the left and right vertical loads, respectively, on each ℎ road wheel calculated at the linearization point. By introducing the track-terrain lateral forces linearization and by considering that = 1 − 2 , the ATV lateral and yaw dynamics is given by: cos( ) + 1 2 2 2̈2 cos( ) + 1 2 2 1̈1 + 1 2 1 ,   The axle force characteristics shows a linear behavior for small slip angles, which tends to saturate for higher values. Based on that, each axle lateral force in Equation (17) is linearized around the linearization point identified by β k i0 : where F k y,i0 , β k i0 and C k i represent the value of the i th axle lateral force, slip angle and cornering stiffness, respectively, calculated at the linearization point. The cornering stiffness C k i is expressed as: where F k z,Li 0 and F k z,Ri 0 represent the left and right vertical loads, respectively, on each i th road wheel calculated at the linearization point. By introducing the track-terrain lateral forces linearization and by considering that α = ψ 1 − ψ 2 , the ATV lateral and yaw dynamics is given by: .. .. ..
where the vector g F k y,i0 , β k i0 , C k i = [g 1 , g 2 , g 3 ] collects all the terms related to the linearization point in Equation (20). The rear unit vehicle speed components, u 2 and v 2 , are kinematically related to the front unit vehicle speed by the following relations: By substituting Equation (25) into Equations (21)-(23), the following non-linear system of four equations is obtained: ψ 2 , α is the states vector and w = C s is the input from the hydraulic steering system.
The ATV dynamics with linear track-terrain lateral forces described in Equation (26) is still non-linear due to the presence of trigonometric functions and states/input products. These non-linearities are then simplified by introducing the first order Taylor expansion around the point P 0 = (x 0 , w 0 ): The matrices A, B and E represent the state, the input and the linearization matrices and contain vehicle geometrical and inertial parameters, axles cornering stiffness and the front unit longitudinal speed u 1 .
The first step for analyzing the system dynamics is the identification of the steadystate equilibrium condition so that . x = 0. A stable ATV steady-state cornering behavior is characterized by a constant hitch angle, i.e., . α = . ψ 1 − . ψ 2 = 0, meaning that the ATV is behaving as a single solid unit with a steering joint locked at a constant hitch angle α. The ATV steady-state equilibrium is then obtained from Equation (27) by: where v 1,ss , . ψ 1,ss = . ψ 2,ss . ψ ss and α ss represent the steady-state front unit lateral speed, vehicle yaw rate and hitch angle, respectively.

Analytical Solution for Small Lateral Accelaratons
The analysis of the linearized system described by Equation (27), and its steady-state solution in Equation (28), is carried out by considering a null linearization point P 0 , i.e., . ψ 2,0 , α 0 = 0 and w 0 = 0. This hypothesis is well representative for small ATV front unit lateral accelerations (a y1 < 1 m/s 2 ). The general expression of Equation (27) can be linked to the characteristic matrices of the dynamic system under investigation: where matrices L, M and N are defined by: Equation (28) provides the steady-state relation from the input steering torque C s to the front lateral speed v 1,ss (or, equivalently, the front sideslip angle β 1,ss ), the ATV yaw rate ψ ss ) and the hitch angle α ss . Nevertheless, the steady-state characteristics β 1,ss vs. α ss and . ψ ss vs. α ss represent a more attractive solution from the steering controller point of view. The explicit dependence of C s is removed by subtracting the third equation to the second one in Equation (29), thus obtaining the relation among the three states β 1,ss , . ψ ss and α ss : where a 11 , a 12 , a 21 , a 22 , b 1 and b 2 are the speed-dependent coefficients defined by: The steady-state response is then computed by solving the algebraic system of equations in Equation (33): where G β1 and G . ψ represent the steady-state sideslip angle and yaw rate gains, respectively. The steady-state gains of rear sideslip angle G β2 , front lateral acceleration G ay1 , front curvature G ρ , front axle slip angles G β I i and rear axle slip angles G β I I i are calculated from G β1 and G . ψ : The ATV linearized model for small lateral accelerations, through Equations (29), (35) and (36), is then analyzed by providing a realistic parametrization in Table 1. It is worth noting that each axle cornering stiffness is evaluated for a constant vertical force equal to the static contribution, which represents a realistic assumption if small lateral accelerations are considered.
Matrices A and B in Equation (29) are still influenced by the ATV front longitudinal speed u 1 . Despite the assumption of constant speed for the linearized mode, u 1 plays an important role for ATV lateral stability and dynamics. The linear model dynamics was then investigated by using the modal analysis, as shown in Figure 9, where the poles, the natural frequencies and the damping ratios are reported as functions of vehicle speed.  All damping ratios are positive, hence the ATV linear model is stable in the whole speed range. The lower natural frequency ω n,1 increases with ATV longitudinal speed up to 65 km/h, where a peak value of 1 rad/s is reached, after which it tends to slightly reduce for higher speed values. The first two poles are real and negative in the whole operative speed range; meanwhile, the third and the fourth eigenvalues are real and distinct for low vehicle speed, while they form a couple of complex conjugate eigenvalues for speeds greater than 33 km/h.
To evaluate the validation range of the linearized ATV model, a step steering torque was applied to both linear and non-linear models to compare their transient and steadystate responses. Simulation results are shown in Figure 10 for three different steering torque histories.  The steering input torques were selected to reach three stable operative points from the steady-state characteristics between steering torque and hitch angle in Figure 7: a step steering torque of 59 Nm (a), a step steering torque of 343 Nm (b), and a smoothed step steering torque of 343 Nm where an overshoot of 432 Nm is imposed during the initial phase. The perfect match between the non-linear and linear responses to the steering torque (a) shows that the assumptions behind the linearization procedure are verified for an ATV operative range characterized by low speed and hitch angle values. When a higher steering torque is applied, e.g., input (b), the starting phase of transient response is still well described by the linear model, but its steady-state response is not aligned with the non-linear model. This result is justified by the selection of the linearization point in the origin of the state space: the more the state deviation from the state space origin, the lower the accuracy level of the linear model in approximating the ATV non-linearities. The non-linear characteristics between the steering torque and the hitch angle is well emphasized by the different response between the linear and the non-linear models to the The steering input torques were selected to reach three stable operative points from the steady-state characteristics between steering torque and hitch angle in Figure 7: a step steering torque of 59 Nm (a), a step steering torque of 343 Nm (b), and a smoothed step steering torque of 343 Nm where an overshoot of 432 Nm is imposed during the initial phase. The perfect match between the non-linear and linear responses to the steering torque (a) shows that the assumptions behind the linearization procedure are verified for an ATV operative range characterized by low speed and hitch angle values. When a higher steering torque is applied, e.g., input (b), the starting phase of transient response is still well described by the linear model, but its steady-state response is not aligned with the non-linear model. This result is justified by the selection of the linearization point in the origin of the state space: the more the state deviation from the state space origin, the lower the accuracy level of the linear model in approximating the ATV non-linearities. The nonlinear characteristics between the steering torque and the hitch angle is well emphasized by the different response between the linear and the non-linear models to the steering input (c). Indeed, the steering input (c) is designed to reach the same final value of the steering step (b) but with a different time history that aims at moving from the left to the right side of the maximum peak in the C s vs. α characteristics shown in Figure 7. Even if the final value of the steering torque (c) is the same as (b), the steady-state hitch angle from the non-linear model is higher; meanwhile, the linear model approaches the same hitch angle.
The steady-state response of the linear model was further analyzed by comparing the yaw rate and sideslip angles characteristics with respect to the non-linear model, as reported in Figure 11: to the right side of the maximum peak in the vs. characteristics shown in Figure 7. Even if the final value of the steering torque (c) is the same as (b), the steady-state hitch angle from the non-linear model is higher; meanwhile, the linear model approaches the same hitch angle.
The steady-state response of the linear model was further analyzed by comparing the yaw rate and sideslip angles characteristics with respect to the non-linear model, as reported in Figure 11: Even if the yaw rate response of the non-linear model is well described by the ATV linearized model for a wider range, the sideslip angles of the front and rear units are only validated for hitch angles lower than 10 deg at a constant speed of 20 km/h. By considering the yaw rate values limited within this range, the liner approximation in Equation (29) is representative of the steady state ATV cornering behavior only for lateral accelerations lower than 1 m/s 2 .
By focusing the analysis within the linear ATV cornering behavior, the influence of the front vehicle speed on steady-state gains, defined in Equation (36), is shown in Figure  12.  Even if the yaw rate response of the non-linear model is well described by the ATV linearized model for a wider range, the sideslip angles of the front and rear units are only validated for hitch angles lower than 10 deg at a constant speed of 20 km/h. By considering the yaw rate values limited within this range, the liner approximation in Equation (29) is representative of the steady state ATV cornering behavior only for lateral accelerations lower than 1 m/s 2 .
By focusing the analysis within the linear ATV cornering behavior, the influence of the front vehicle speed on steady-state gains, defined in Equation (36), is shown in Figure 12.
Even if the final value of the steering torque (c) is the same as (b), the steady-state hitch angle from the non-linear model is higher; meanwhile, the linear model approaches the same hitch angle.
The steady-state response of the linear model was further analyzed by comparing the yaw rate and sideslip angles characteristics with respect to the non-linear model, as reported in Figure 11: Even if the yaw rate response of the non-linear model is well described by the ATV linearized model for a wider range, the sideslip angles of the front and rear units are only validated for hitch angles lower than 10 deg at a constant speed of 20 km/h. By considering the yaw rate values limited within this range, the liner approximation in Equation (29) is representative of the steady state ATV cornering behavior only for lateral accelerations lower than 1 m/s 2 .
By focusing the analysis within the linear ATV cornering behavior, the influence of the front vehicle speed on steady-state gains, defined in Equation (36), is shown in Figure  12.  The gains evaluated at zero longitudinal speeds represent the "kinematic" ATV steadystate behavior where the effect of inertial forces, due to lateral accelerations, are negligible. The non-null kinematic front and rear slip angle gains show that all the four axles slip on the terrain even at low or zero speed conditions. Moreover, from the curvature gain plot, the ATV shows an evident steady-state understeer behavior: the same value of hitch angle produces a larger curvature radius (inverse of the curvature) and higher front and rear sideslip angles (in module) for higher speeds.
The ATV cornering behavior is also well described by the equivalent understeer characteristics, typical of passenger vehicles. For the ATV, the understeer characteristics represents the dynamic hitch angle correction, with respect to the kinematic condition, required to keep a desired curvature radius for different lateral accelerations. The kinematic hitch angle, α kin = ρ G ρ,kin , that must be set to follow a desired trajectory at lower speeds (or lateral accelerations) is calculated by considering the kinematic curvature gain G ρ, kin : It must be noted that the single-track model assumptions commonly used for the handling analysis of passenger cars lead to a kinematic gain G ρ, kin equal to the inverse of the vehicle wheelbase. On the other hand, the ATV kinematic curvature gain depends on front and rear axle cornering stiffness as well as on vehicle geometrical parameters.
The dynamic hitch angle α dyn is then calculated as the difference between the current and the kinematic hitch angles: The understeer gradient, K US , is positive for understeering and negative for oversteering ATV behaviors. It can be analytically proved that the understeer gradient of Equation (38) is speed-independent for the linearized model, but the understeering behavior of the non-linear model is strongly influenced by the vehicle speed, especially at higher lateral accelerations, as shown in Figure 13. axles slip on the terrain even at low or zero speed conditions. Moreover, from the curvature gain plot, the ATV shows an evident steady-state understeer behavior: the same value of hitch angle produces a larger curvature radius (inverse of the curvature) and higher front and rear sideslip angles (in module) for higher speeds.
The ATV cornering behavior is also well described by the equivalent understeer characteristics, typical of passenger vehicles. For the ATV, the understeer characteristics represents the dynamic hitch angle correction, with respect to the kinematic condition, required to keep a desired curvature radius for different lateral accelerations. The kinematic hitch angle, = , , that must be set to follow a desired trajectory at lower speeds (or lateral accelerations) is calculated by considering the kinematic curvature gain , : It must be noted that the single-track model assumptions commonly used for the handling analysis of passenger cars lead to a kinematic gain , equal to the inverse of the vehicle wheelbase. On the other hand, the ATV kinematic curvature gain depends on front and rear axle cornering stiffness as well as on vehicle geometrical parameters.
The dynamic hitch angle is then calculated as the difference between the current and the kinematic hitch angles: The understeer gradient, , is positive for understeering and negative for oversteering ATV behaviors. It can be analytically proved that the understeer gradient of Equation (38) is speed-independent for the linearized model, but the understeering behavior of the non-linear model is strongly influenced by the vehicle speed, especially at higher lateral accelerations, as shown in Figure 13. In the range of low lateral accelerations ( 1 ≤ 1 m/s 2 ), the understeer gradient of the non-linear model is also speed-independent and aligned with the linear model. For higher In the range of low lateral accelerations (a y1 ≤ 1 m/s 2 ), the understeer gradient of the non-linear model is also speed-independent and aligned with the linear model. For higher lateral accelerations (a y1 > 1 m/s 2 ), the non-linear model shows a more understeering behavior at 40 km/h and an oversteering behavior at 20 km/h if compared to the linear model. However, the magnitude of the dynamic hitch angle obtained with the non-linear model at high lateral accelerations is not so significant if compared to the absolute hitch angle α ss . For example, at a y1 = 3.5 m/s 2 the absolute hitch angle obtained with the non-linear model (at u 1 = 20 km/h) is α NL ss = 30 deg.

Controller Design
While the previous section analyzed the ATV dynamic behavior through the linearized model, the present section evaluates the influence of the hydraulic actuation system by introducing a hitch angle controller. The hitch angle controller was designed based on the non-linear model, due to the presence of strong nonlinearities that characterize the ATV and the hydraulic pressure dynamics. Two alternatives logics are proposed for controlling the ATV hitch angle, as shown in Figure 14. linear model (at = 20 km/h) is = 30 deg.

Controller Design
While the previous section analyzed the ATV dynamic behavior through the linearized model, the present section evaluates the influence of the hydraulic actuation system by introducing a hitch angle controller. The hitch angle controller was designed based on the non-linear model, due to the presence of strong nonlinearities that characterize the ATV and the hydraulic pressure dynamics. Two alternatives logics are proposed for controlling the ATV hitch angle, as shown in Figure 14. Integral Derivative (PID) control logic, where represents the desired hitch angle. One PID logic is designed to elaborate a steering torque , named "Torque-Based Logic" (TBL), and the second PID calculates a hydraulic command valve , named as "Speed-Based Logic" (SBL). In the first case, the hydraulic steering dynamics is excluded by considering an ideal actuator that directly provides the desired steering torque, as assumed for the ATV linearized model and in [22], where an electric motor is directly mounted on the steering joint. For the SBL, the steering torque applied to the ATV depends on the pressure dynamics within the hydraulic circuit of the steering system. The reason behind the selection of these two logics aims at analyzing the impact of a direct steering torque on the ATV lateral dynamics if compared to the steering action applied by the hydraulic circuit, which is influenced by the internal feedback between the actuator and system dynamics. A speed controller is also included in Figure 14 to reduce the speed error = − by regulating the gas pedal position (not described in this paper). Furthermore, all the feedback signals shown in Figure 14 are assumed to be fully available but feasible estimators are required for practical implementations, as also described in [23]. The two alternative logics are expressed by: Both hitch angle controllers use the error e α = α re f − α as input for a Proportional Integral Derivative (PID) control logic, where α re f represents the desired hitch angle. One PID logic is designed to elaborate a steering torque C s , named "Torque-Based Logic" (TBL), and the second PID calculates a hydraulic command valve X v , named as "Speed-Based Logic" (SBL). In the first case, the hydraulic steering dynamics is excluded by considering an ideal actuator that directly provides the desired steering torque, as assumed for the ATV linearized model and in [22], where an electric motor is directly mounted on the steering joint. For the SBL, the steering torque applied to the ATV depends on the pressure dynamics within the hydraulic circuit of the steering system. The reason behind the selection of these two logics aims at analyzing the impact of a direct steering torque on the ATV lateral dynamics if compared to the steering action applied by the hydraulic circuit, which is influenced by the internal feedback between the actuator and system dynamics. A speed controller is also included in Figure 14 to reduce the speed error e u = u re f − u 1 by regulating the gas pedal position X g (not described in this paper). Furthermore, all the feedback signals shown in Figure 14 are assumed to be fully available but feasible estimators are required for practical implementations, as also described in [23]. The two alternative logics are expressed by: The PID's gains K P,TBL , K D,TBL , K I,TBL , K P,SBL , and K I,SBL were obtained through the Simulink® "Response Optimization" toolbox. The optimizations process is designed to achieve the desired response to a step input on α re f at a desired ATV speed u 1 , kept constant by the speed controller shown in Figure 14, based on the requirements reported in Table 2. Optimal solution for the set of PID gains depends on the boundary conditions for the optimization problem, such as the ATV speed u 1 and the final hitch angle value α re f , f inal . An example for optimal PID gains selection is provided in Table 3. The comparison between the two logics shows that the TBL requires an extreme PID adaptation, in terms of gains magnitude and sign variation, to the maneuver operative conditions. This is mainly due to the strong non-linear characteristics between the steering torque and the hitch angle. On the other hand, the SBL does not require such a drastic PID gains variation since the hydraulic steering circuit automatically adapts the fluid pressure to provide the needed steering torque. This result implies that the SBL structure is more stable than the TBL. A further verification of this concept is represented by the ATV non-linear model response to a step input when a fixed set of PID gains, e.g., the first row in Table 3, is applied for each final hitch angle values. Results are shown in Figure 15.

Controller Validation
The hitch angle controller with SBL design for , = 1 deg (first row in Table 3) is then validated through a more realistic multiple turn maneuver. Figure 16 reports the time histories of the main lateral dynamics' quantities during a multiple turn maneuver at a constant speed of 20 km/h. The simulated maneuver consists of a first swift Even if the PID gains are not adapted for the three operative conditions, the SBL provides a clearly better response if compared to the TBL, which does not respect all the desired performance requirements reported in Table 2, especially in terms of maximum settling time and overshoot. Although the stability of both logics was not assessed through an analytical formulation, the time domain responses described in Figure 15 represent a numerical verification of the controller performance and stability. The SBL shows a more stable response to the reference hitch angle variation with respect to the TBL, which tends to reduce the ATV stability when the amplitude of the reference signal rises, as can be seen by the increase of hitch angle oscillations.

Controller Validation
The hitch angle controller with SBL design for α re f , f inal = 1 deg (first row in Table 3) is then validated through a more realistic multiple turn maneuver. Figure 16 reports the time histories of the main lateral dynamics' quantities during a multiple turn maneuver at a constant speed of 20 km/h. The simulated maneuver consists of a first swift turn with a curvature radius of 40 m counterclockwise, a second one that brings the ATV on a straight line, a third one with a curvature radius of 20 m clockwise, and a final one that brings back the ATV on a straight line.

Figure 15.
Step response of the non-linear ATV model with TBL and SBL at 1 = 20 km/h.

Controller Validation
The hitch angle controller with SBL design for , = 1 deg (first row in Table 3) is then validated through a more realistic multiple turn maneuver. Figure 16 reports the time histories of the main lateral dynamics' quantities during a multiple turn maneuver at a constant speed of 20 km/h. The simulated maneuver consists of a first swift turn with a curvature radius of 40 m counterclockwise, a second one that brings the ATV on a straight line, a third one with a curvature radius of 20 m clockwise, and a final one that brings back the ATV on a straight line. Figure 16. Sideslip angles, yaw rates, hitch angle, lateral accelerations, curvatures, steering torque, and hydraulic command valve during a multiple turn maneuver, with the hitch angle controller with SBL and fixed PID gains applied to the non-linear model at constant ATV speed 1 = 20 km/h. Reprinted with permission from ref. [19]. Copyright 2020 Springer Nature Switzerland AG.
The ATV well follows the variable hitch angle required even if a fixed PID gains is adopted for the SBL. It is also remarkable to see that, given the strong non-linear relation between the steering torque and the hitch angle shown in Figure 7, the steering torque required to track the reference hitch angle during a multiple turn maneuver is extremely variable despite the maneuver not being considered as aggressive as proven by the low lateral acceleration condition. Figure 16. Sideslip angles, yaw rates, hitch angle, lateral accelerations, curvatures, steering torque, and hydraulic command valve during a multiple turn maneuver, with the hitch angle controller with SBL and fixed PID gains applied to the non-linear model at constant ATV speed u 1 = 20 km/h. Reprinted with permission from ref. [19]. Copyright 2020 Springer Nature Switzerland AG.
The ATV well follows the variable hitch angle required even if a fixed PID gains is adopted for the SBL. It is also remarkable to see that, given the strong non-linear relation between the steering torque and the hitch angle shown in Figure 7, the steering torque required to track the reference hitch angle during a multiple turn maneuver is extremely variable despite the maneuver not being considered as aggressive as proven by the low lateral acceleration condition.

Conclusions
The activity presented in this paper describes a non-linear model for a typical ATV and a methodology to linearize the model for the lateral dynamics analysis. The following conclusions are drawn from simulation results: • The linearized ATV model represents a promising tool for ATV dynamics analysis only in the range of small lateral accelerations (a y1 < 1 m/s 2 ).
• The linearized ATV model shows a stable behavior over the whole speed operative range, and it is characterized by a "slow" open-loop dynamics due to the presence of relatively small natural frequencies.

•
The steady-state ATV cornering behavior for small lateral accelerations is understeering with the proposed set of vehicle parameters. This is also proven by extending to the ATV the concept behind the understeer characteristics commonly adopted for passenger cars. • At higher lateral accelerations, the ATV steady-state understeering attitude obtained with the non-linear model is strongly influenced by the vehicle speed. An undesteering behavior is observed for 40 km/h and an oversteering attitude is obtained for 20 km/h.

•
The hitch angle controller with SBL represents a better solution if compared to the TBL, since the hydraulic steering circuit automatically adapts the pressure dynamics to the ATV load conditions without any direct steering torque regulation.
Future investigations are necessary to propose a similar methodology to design the hitch angle controller for higher lateral accelerations and in the presence of ATV parameters variation, e.g., the axles' cornering stiffness. The enhancement of the linearized model represents a fundamental requirement for extending the analytical approach to the hitch angle controller design and its stability and performance assessment as well as its robustness level against uncertain parameters and external disturbances. Moreover, other future developments of this work aim at improving the analytical methodology by including (a) the hydraulic actuation dynamics for enhancing the ATV linearized model accuracy in predicting the ATV lateral behavior, (b) the effect of the ATV track tension on slip angles generation and (c) the modification of the track-terrain contact forces when a soft-soil interaction is required with a non-negligible track sinkage.

Conflicts of Interest:
The authors declare no conflict of interest.

Appendix A
The aim of this appendix is to provide the non-linear kinematic relation imposed by the steering joint geometry as shown in Figure A1.
The right and left rod displacement, ∆X R = BC − BC and ∆X L = AD − AD are evaluated by considering the trianglesBOC andÂOD: V 2 are then calculated as: where V 0 is the total volume when the ATV is moving in a straight line (α = 0). Finally, the right and left piston force arms f R and f L with respect to the steering hinge are calculated by: