Research Article  Open Access
Robust Adaptive Fuzzy Control for Planetary Rovers While Climbing up Deformable Slopes with Longitudinal Slip
Abstract
Mobility control is one of the most essential parts of planetary roversâ€™ research and development. The goal of this research is to let the planetary rovers be able to achieve demand of motion from upper level with satisfied control performance under the rough and deformable planetary terrain that often lead to longitudinal slip. The longitudinal slip influences the mobility efficiency obviously, especially on the major deformable slopes. Compared with the past works on normal stiff terrains, properties of soil and interaction between wheels and soil should be considered additionally. Therefore, to achieve the final goal, in this paper, wheelsoil dynamic model for sixwheel planetary rovers while climbing up deformable slopes with longitudinal slip is first built and control based in order to account for slip phenomena. These latter effects are then taken into account within terramechanics theory, relying upon nonlinear control techniques; finally, a robust adaptive fuzzy control strategy with longitudinal slip compensation is developed to reduce the effects induced by slip phenomena and modeling error. Capabilities of this control scheme are demonstrated via full scale simulations carried out with a sixwheel robot moving on sloped deformable terrain, whose real time was computed relying uniquely upon RoSTDyn, a dynamic software.
1. Introduction
In the field of special mobile robots environment, including planetary exploration missions, caravan survey, polar expedition, and wild fire spreading, rovers may need to traverse on deformable terrains, and the interaction between rigid wheels and soft soil has become a meaningful research topic because of longitudinal slip influence mobility control obviously [1]. In the past works on normal stiff terrains, for example, Kanayama et al. [2] proposed a stable control scheme for an autonomous mobile robot under the assumption of perfect velocity tracking. Kim and Oh [3] proposed a modified inputoutput linearization method to solve the problem of a decoupling matrix using a generalized inverse that provided a leastsquares solution to the tracking control of twowheeled mobile robots. Raibert et al. [4] proposed a PID controller to solve the path tracking problem of a mobile robot using a simple linearized model of the mobile robot, which was composed of an integrator and a delay. Colombano et al. [5] proposed an outputfeedback controller that allowed a unicycle mobile robot to track a predefined path. However, all of these control methods based on normal stiff terrains hypothesis of nonholonomic mobile robot systems. Generally, there are two aspects concerning the control accuracy for rovers that traverse on deformable terrains. One is the modeling error. Built on terramechanics, the interaction dynamic model typically has a complicated integral and nonlinear form and must be simplified for utilization in control design. Modeling error is then unavoidable and should be kept in an acceptable range. The other is the uncertain parameters including soil parameters and mechanism parameters in the model. These soil parameters and mechanism parameters are timevariant as rover moves at various locations of deformable terrain. Possibility to deal with such problems is robust control and adaptive parameters estimation [1]. Among the system, many modern control theory techniques and stabilizing control algorithms have been proposed for the motion control of special mobile robots. The neural network control of mobile robots has been the subject of intense research in recent years [6, 7]. Bartsch et al. [6] investigated tracking control using an adaptive smart neural network for robot and they produced fine motion control based on partially unknown dynamics. Bekker [8] presented a neural networkbased model that combined the backstepping technique with a torque controller. For the posterior mobility control of rover, wheelterrain interaction model is addressed. In the relevant literature, McKyes provide wheelterrain mechanics with closed form [9]. The model describes a rigid wheel moving on deformable terrain with low cohesion as shown in Figure 4. Because the terrain is assumed as a plane, the dynamic model contains two force equations, normal stress and drawbar pull, and one driving torque equation. Such a model is decomposed by Terzaghi into three sequent formulations for the sake of convenience to identify the parameters [10]. Reece compared the simplified and original integral form by numerical analysis [11, 12]. Based on the wheelterrain mechanics, Brunskill and Lappas, propose online soil parameter estimation using the linear leastsquare method [13, 14]. Janosi and Hanamoto adopt function fitting method to approximate the integrals of the original model and then NewtonRaphson method to find the solution of parameters in the nonlinear equations numerically [15â€“18].
However, the abovementioned models are all built on flat terrain, although deformation of terrain is considered. In practical missions, slope terrain is hard to avoid. A slopebased wheelsoil dynamic model based on control for WMR is required to be built when the rover wheels are inclined [19, 20]. In this case, the model with two forces and one torque for flat terrain cannot fully describe the wheelterrain dynamics on a slope. On the other hand, noticing that the previous control method, such as robust control and realtime soil parameters estimation, is strictly taking the wheelterrain dynamics model into consenter, it may cost excessive time to find the satisfactory solution. This is unnecessary and hard to be implemented.
In this paper, we develop a slopebased wheelsoil dynamic model for WMR based on a simplified interaction model and uniwheel slopebased experiment, in which the drawbar pull can be denoted as a linear function of driving force, normal force, and slip; it is suitable for describing the rigid wheels on sloped deformable terrain. Then, the simplified slopebased wheelsoil interaction dynamic model is used for WMR slopebased control. Based on the simplified slopebased wheelsoil interaction dynamic model, control strategy can be designed simultaneously by adopting an adaptive fuzzy and slip compensation [8]. Robust control is used because the unmodeled error should be compensated timely by adaptive adjustment of the index function. The uncertain parameters include soil parameters and mechanism parameters, which make the state uncontrollable and are compensated by adaptive fuzzy strategy. The slip is the important state variable; in order to deal with this problem, we deemed as measurement model using the methods like the relevant literature. Capabilities of this control scheme is demonstrated via full scale simulations carried out with a sixwheel robot moving on sloped deformable terrain, whose real time computed relying uniquely upon RoSTDyn, a dynamic software.
This paper is organized as follows. The slopebased wheelsoil dynamic model for uniwheel is introduced in Section 2. In Section 3, the slopebased wheelsoil dynamic model for planetary rover is developed. The robust adaptive fuzzy control algorithm for planetary rover on sloped deformable terrain is established in Section 4, followed by the stability analysis in Section 5. The examples of simulations are detailed in Section 6.
2. SlopeBased WheelSoil Dynamic Model for Uniwheel
Assume that the wheels of a planetary rover are traversing perpendicularly on a sloped deformable terrain. The slope is positioned in the inertial frame . The local wheel frame is defined as , where is along the forward direction, along the lateral, and the normal. , , and represent the pitch, roll, and yaw angle of the wheel in the inertial coordinate system, as shown in Figure 1.
When roll angle Â°, the pitch angle makes the direction of the vertical load of the wheel deviate from that of the perpendicular stress to the slope; nevertheless, it does not violate the validity of terramechanics theory [2â€“4] to calculate the normal force , drawbar pull ; and wheel driving torque . Let and denote the entrance and exit angle which the wheel contact the soil; let and be the maximum and residual sinkage, respectively. Interaction between the wheel and terrain generates the normal stress toward the wheel center and shearing stress along the wheel tangential direction at every contact point. According to the terramechanics theory [2], these two can be calculated by where the soil sinkage exponent is a linear function of the slip ratio according to [3] and can be determined as in which and are real coefficients. and are the cohesive modulus and frictional modulus of sinkage. denote the wheel width. The angle of maximal stress is expressed by In (2), is the internal cohesion of the soil, and is the friction angle; denote the shear modulus. The slip ratio in (2)â€“(4) is defined as which explicates the relationship between angular rate and actual forward velocity . is the radial of the wheel. However, as one can see from Figure 2, different from the case of plane terrain, the entrance and exit angles are calculated by
For a moving wheel, there are three forces and three torques dynamic balance equations to completely describe the motion on slope. In Figure 2, it is obvious that the lateral force , the turning torque is uncontrolled, and roll torque . is passive and uncontrollable and should be monitored for the rollover avoidance. Thus, there are three equations that can be derived: normal stress , drawbar pull , and driving torque . From Figure 2, they have the following integral expression: in which is the moment of inertia of the wheel rotational axis.
3. SlopeBased WheelSoil Dynamic Model for Planetary Rover
When a planetary rover is traversing perpendicularly on a sloped deformable terrain, analysis force of slopebased wheelsoil is just like Figure 3. So the slopebased wheelsoil dynamic model was developed based on NewtonEuler method as (8)
Due to the force , are integrated model which is so complex that cannot be used for controller design directly; next, the integrated model should be simplified. Based on what our research had done, the force and model can be simplified as (9) and (10). Consider
denote soil thrust caused by electrical machine drive wheel leading to soil deformation; denote soil resistance, where is the coefficient of wheel resistance denoted by RC; RC reflects the counter of soil to wheel. The drowpull of the wheel is the thrust subtract the soil resistance. in formula (10) and (7) are accordant, so put the formula (10) into (8), the controlbased equation (11) is obtained as follows:
Analyzing the curve of RC, one can see from Figure 5 and Table 1 that RC can be donated as a linearity function, approximately, calculated by (12). The value of RC is 0.08â€“0.32, is 0.09â€“0.17, and is 0.20â€“0.26 approximately, from our previous experimentation. Consider

If one define and , from (11) and (12), the slopebased wheelsoil dynamic model for a planetary rover can be expressed as (13) where
4. Robust Adaptive Fuzzy Controller Design
The planetary rover shown in Figure 3 is a typical example of a nonholonomic mechanical system. It has six parallel driving wheels mounted on the same axis. The motion and orientation are controlled by independent actuators; that is, DC motors provide the necessary torque to the six wheels. The position of the robot in an inertial Cartesian frame is completely specified by the vector, where and are the coordinates of the center of the mass and is the local coordinate with an origin of with respect to the inertial basis. One can see from literature that kinematics equality constraints are independent with time and can be expressed as follows:
Let be a dimension and full rank matrix, formed by a set of smooth and linear independent vector fields spanning the null space of :
According to (15) and (16), it is possible to find an auxiliary vector time function such that for all time
Differentiating (17) yields
Then, the dynamic (13), which satisfies the nonholonomic constraint (15), can be rewritten in terms of the internal state variable as
By appropriate definitions one can rewrite equation as where
We assume that the external disturbances can be expressed as where is a known matrix.
Exploiting the structure of the dynamic equation (20), some properties are listed as follows.
Property 1. Matrix is symmetric and positivedefinite.
Property 2. Matrix is skewsymmetric.
Suppose that the desired trajectories in the task space are described by and . In order to track the desired trajectories, we define the sliding surface: where is the tracking error and is a constant positive definite diagonal matrix. Then define the reference state as
One now define a robust control input based on the sliding mode theory. Thus, we extend the sliding mode control law given in [15] and applied for planetary rover which is traversing perpendicularly on a sloped deformable terrain where and are constant positive definite diagonal matrices. The proposed control input consists of a continuous nominal control part and a discontinuous switching control part to make the system invariant to parametric uncertainty and external disturbances. The switching control part causes the chattering problem. To reduce chattering, the switching function sign() can be replaced by a continuous approximation in the neighbourhood of the sliding surface, such as sigmoidlike function. A complete stability analysis of the sigmoidlike function is given in [17].
However, the error of the unmodel influences the performance of the proposed sliding mode controller. In the literature of NN control, many different models of function approximations such as Gauss function have been established to approximate the actual condition, because of their inherent approximation capabilities. The radial basis function neural network is employed for function approximation, what is assumed as a smooth function where is the input vector, contains the ideal thresholds and weights ; and is the activation vector and is a vector of the function approximation errors. Commonly, the Gaussian function is chosen as the activation vector where and are the center and width of the th neuron.
To satisfy the universal approximation property for any continuous function , where is a compact set, there exists an estimate function such that where are optimal weights in the approximation, is the approximation error for the special case , and is a constant.
To eliminate the error influence on planetary rover traversing perpendicularly on a sloped deformable terrain, the following NN control law is considered: in which the continuous approximation law identifier is designed as where is a positive definite matrix, and the adaptive fuzzy controller is designed as with the individual fuzzy law as follows:â€‰IF is NB, THEN is NB;â€‰IF is NM, THEN is NM;â€‰IF is NS, THEN is NS;â€‰IF is ZE, THEN is ZE;â€‰IF is PS, THEN is PS;â€‰IF is PM, THEN is PM;â€‰IF is PB, THEN is PB;
where NB, NM, NS, ZE, PS, PM, and PB are fuzzy sets. BS stands for negative, P positive, B big, M medium, S small, and Z zero; they are all Gaussian membership functions defined as follows: can be written as in which is the amount of the rules.
The configuration of the proposed robust adaptive fuzzy system, which combines a kinematics controller proposed in , is summarized in Figure 3.
5. Stability Analysis
In this section, we do the system stability analysis for the closedloop behavior of the proposed control methodology as shown in Figure 6. Define the reference position , so that is the optimal compensation for ; according to Wangâ€™s theorem, there exists satisfying where can be as small as possible.
Define Therefore, (34) can be written as follows: Consider the following Lyapunov function candidate: Differentiating yields where Assume where ; then Putting (42) into (40), we have where , . , are guaranteed negative and this shows that and implies and as . Furthermore, (28) shows , if and only if . Therefore, global stability is guaranteed by the Lyapunov theorem. Both the system stability and velocity tracking convergence are guaranteed by control law (29), (31), and (32) driving the system (40) to closely track the desired motion trajectories.
6. Simulation Results
When a planetary rover is traversing perpendicularly on a sloped deformable terrain as shown in Figure 7, slippage occurs between the wheels and the ground, and most of the wheel velocities are influenced by the slip ratio. This section, the robust adaptive fuzzy control strategy with slip compensation designed upside on this paper, is demonstrated via full scale experiments carried out with the sixwheel robot moving on sloped deformable terrain. The feedback method of slip ratios is acquisition from [11]. Therefore, we assume that the reference trajectory is a line, where , , and and where the initial position is . The slope of the experimentation working condition is set as 0Â°, 8Â°, 15Â°, and 18Â°. The linear speed of planetary rover is set as â€‰m/min.
When a planetary rover is traversing perpendicularly on a 15Â° slope deformable terrain, the comparison of capabilities with/without wheelsoil dynamic model is denoted by Figure 8; form the capabilities one can see that the longitudinal slip influences the planetary rover mobility obviously.
(a)
(b)
(c)
(d)
(e)
The response of this controller is shown in Figure 9. We compensated for lost velocity due to wheel slippage. It is obvious that the tracking errors were almost the same as . The results show that slipcompensation plays a significant role. Thus, when slippage occurred between the wheels and the ground, the control method with slipcompensation and adaptive NN control laws performed well during tracking control.
(a)
(b)
(c)
(d)
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
This work was supported by the National Natural Science Foundation of China (50975059/61005080), China Postdoctoral Science Foundation (20100480994), and Special Foundation (201104431) and the â€ś111â€ť Project (B07018).
References
 NASA Mars Science Laboratory Factsheet, http://www.nasa.gov/centers/goddard/news/topstory/2008/radiometer_delivery.html.
 European Space Agency, ExoMars Science Management Plan: Doc. No: EXMMSPLESA00002, 12th February 2010.
 Google Lunar X Competition Rules, 2011, http://www.googlelunarxprize.org/les/downloads/lunar/GLXP Guidelines v3 Nov 20 2008.pdf.
 M. Raibert, K. Blankespoor, G. Nelson, and R. Playt, â€śBigDog, the roughterrain quaduped robot,â€ť in Proceedings of the17th World Congress: The International Federation of Automatic Control, Seoul, Korea, July 2008. View at: Google Scholar
 S. Colombano, F. Kirchner, and D. Spenneberg, â€śExploration of planetary terrains with a legged robot as a scout adjunct to a rover,â€ť in Proceedings of the Institute of Aeronautics and Astronautics, Space Conference, San Diego, Calif, USA, 2004. View at: Google Scholar
 S. Bartsch, T. Birnschein, F. Cordes et al., â€śSpaceClimber: development of a sixlegged climbing robot for space exploration,â€ť in Proceedings of the 41st International Symposium on Robotics, pp. 1265â€“1272, June 2010. View at: Google Scholar
 G. P. Scott and C. M. Saaj, â€śMeasuring and simulating the effect of variations in soil properties on microrover tracability,â€ť in Proceedings of the SPACE Conference, American Institute of Aeronautics and Astronautics, Pasadena, Calif, USA, 2009. View at: Google Scholar
 M. Bekker, Introduction to Terrain Vehicle Systems, Part 1—The Terrain & Part 2—The Vehicle, vol. 2, University of Michigan Press, Ann Arbor, Mich, USA, 1959.
 E. McKyes, Soil Cutting and Tillage, Elsevier Science, Amsterdam, The Netherlands, 1985.
 K. Terzaghi, Theoretical Soil Mechanics, John Wiley & Sons, London, UK, 3rd edition, 1943.
 A. R. Reece, â€śThe fundamental equation of earthmoving mechanics,â€ť Proceedings of the Institution of Mechanical Engineers, Conference Proceedings, vol. 179, no. 6, pp. 16â€“22, 1964. View at: Publisher Site  Google Scholar
 R. D. Grisso and J. V. Perumpral, â€śA soil tool interaction model for narrow tillage tools,â€ť American Society of Agricultural Engineers (ASAE), 1980. View at: Google Scholar
 C. Brunskill and V. Lappas, â€śThe effect of relative soil density on microrover tracability under low ground pressure conditions,â€ť in Prceedings of the 11th European Regional Conference of the International Society for TerrainVehicle Systems, Bremen, Germany, 2009. View at: Google Scholar
 J. Y. Wong, Theory of Ground Vehicles, John Wiley & Sons, New York, NY, USA, 4th edition, 2008.
 Z. Janosi and B. Hanamoto, â€śThe analytical determination of drawbar pull as a function of slip for tracked vehicles in deformable soils,â€ť in Proceedings of the ISTVS 1st International Conference on Mechanics of SoilVehicle Systems, pp. 707â€“736, Edizioni Minerva Tecnica, Torino, Italy, 1961. View at: Google Scholar
 R. Godbole and R. Alcock, â€śA device for the in situ determination of soil deformation modulus,â€ť Journal of Terramechanics, vol. 32, no. 4, pp. 199â€“204, 1995. View at: Publisher Site  Google Scholar
 A. R. Reece, Problems of SoilVehicle Mechanics, US Army Land Locomotion Lab, ATAC, Warren, Mich, USA, 1964.
 L. Richter, A. Ellery, Y. Gao, S. Michaud, N. Schmitz, and S. Weiss, A Predictive WheelSoil Interaction Model for Planetary Rovers Validated in Testbeds and Against MER Mars Rover Performance Data, Budapest, Hungary, 2006.
 M. Lyasko, â€śSlip sinkage effect in soilvehicle mechanics,â€ť Journal of Terramechanics, vol. 47, no. 1, pp. 21â€“31, 2010. View at: Publisher Site  Google Scholar
 A. Ellery, â€śEnvironmentrobot interaction—the basis for mobility in planetary microrovers,â€ť Robotics and Autonomous Systems, vol. 51, no. 1, pp. 29â€“39, 2005. View at: Publisher Site  Google Scholar
Copyright
Copyright © 2014 Li Zhengcai and Wang Yang. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.