Abstract
Recently, the specific problems arising to mobile robots are the following: determining the position and orientation of the robot on the field, planning an optimal path of motion and stationary or moving objects avoidance.
This paper realizes a computational analysis of the Cauchy problem associated to a mobile robot kinematics. The ‘phase portrait’ graphical tool of the mathematical soft MAPLE11 points out the influence of the initial conditions: the initial velocities of the driving (left and right) wheels of the robot on the robot trajectory. Considering a pair of two simulation cases for the initial conditions brings a good reliability of the analysis.
MSC: 68T40, 70B15, 70B10, 37M05.
Keywords:
mobile robot; differential drive; kinematic model; phaseportrait1 Introduction
The basic goal of this work is to get a consistent computational analysis of the trajectory of a twowheeled differential robot. This implies a few levels: analytical, computational and experimental. For the moment, we shall focus on the computational level.
The mobile robot placing in a workspace represents the definition of the position and orientation possibilities that can be achieved in that space. Controllability of a mobile robot defines the possible trajectories from its workspace. In order to move in the workspace, the mobile robots need certain mechanical components to allow precise movement (or movement would be chaotic). Currently, there are several components, including wheel, track and legs. The most popular designs are the wheeled mobile robots. The wheeled mobile robots are also better controlled than other types of robots. Disadvantages of wheeled robots are that they cannot navigate well over obstacles such as rocky terrain, sharp slopes, or areas with low friction. Usually, robots of this type are used during competitions (e.g., a soccer game, a sumo fight etc.) and are most popular in the consumer market. Robots can have any number of wheels, but three wheels are sufficient for static and dynamic balance. Additional wheels can add to balance; however, additional mechanisms will be required to keep all the wheels on the ground when the terrain is not flat.
Mobile robot kinematics is defined by the number, type and arrangement of the wheel mobile platform. Combination and arrangement of these different types of wheels determine the type of mobile robot kinematic model. The most popular design among the mobile robots is a twowheeled robot with differential drive. This differential motion mechanism has the following advantages:
– ensures smooth motion;
– achieved by propelling and directing by the angular speed difference between the wheels.
Besides the advantages mentioned above, the mechanism of differential motion has some disadvantages:
– the transverse displacements impossibility;
– the risk of slipping.
A differential twowheeled robot is a mobile robot whose movement is based on two separately driven wheels placed on either side of the robot body. It can thus change its direction by varying the relative rate of rotation of its wheels and hence does not require an additional steering motion. One or more free wheels (or ‘castor’ wheels) assure the robot equilibrium. Each castor wheel is independently mounted on a vertical nondrive axis of the body and it is automatically and freely aligned on the route as a result of the forces developed by the two ‘drive wheels’. If both wheels are driven in the same direction and at the same speed, the robot will go in a straight line. Otherwise, depending on the speed of rotation and its direction, the centre of rotation may fall anywhere in the line that joins the two wheels. Since the direction of the robot is dependent on the rate and direction of rotation of the two driven wheels, these quantities should be sensed and controlled precisely. This usually creates some problem. If both wheels are turned with equal speed in opposite directions, the robot will rotate about the central point of the axis.
In what follows, a few mechanisms found in the structure of the mobile robots are presented.
(a) A widespread mobile robot is the twowheeled mobile robot. Its kinematic scheme is presented in Figure 1.
Figure 1. Kinematical scheme of the twowheeled mobile robot.
An important component present in this structure type is ICRinstantaneous center of rotation, with Rradius from ICR.
The velocity of each wheel can be varied. The robot must rotate about a point that lies along the common axis of left and right wheels in order to perform rolling motion. This point that the robot rotates around is known as the ICR (Instantaneous Center of Rotation).
There are the following three interesting cases concerning the left and right velocities and :
1. If , then and the robot moves in a straight line.
2. If , then and the robot rotates about the midpoint of the wheel axis.
3. If , then the robot rotates about the left wheel and if , then the robot rotates about the right wheel. In this case .
The synchro drive configuration is a popular arrangement of wheels in indoor mobile robot applications. In a synchronous drive robot, each wheel is capable of being driven and steered, as depicted in Figure 2. Typical configurations have three steered wheels arranged as vertices of an equilateral triangle often mounted by a cylindrical platform. All the wheels turn and drive in unison [1].
Figure 2. Synchronous drive robot.
A typical tricycle drive robot has three wheels with odometers on the two rear wheels with steering and power provided through the front wheel (as depicted in Figure 3).
Figure 3. Tricycle drive robot.
The omnidirectional robots are able to move in any direction at any time and they are also holonomic. They can be realized by either using spherical, castor, or Swedish wheels. The omnidirectional robot depicted in Figure 4 is based on three spherical wheels, each activated by one motor. In this design, the spherical wheels are suspended by three contact points, two given by spherical bearings and one by a wheel connected to the motor axle. This concept provides excellent maneuverability and is simple in design. However, it is limited to flat surfaces and small loads, and it is quite difficult to find round wheels with high friction coefficients.
Figure 4. Omnidirectional drive robot.
2 The kinematicmathematical model of the twowheeled differential drive mobile
In what follows some considerations regarding the mathematical models for a class of mobile robots, namely twowheeled differential drive mobile robots, some of the most utilized mechanical structures in mobile robotics practice, are presented.
In this section the basic mathematical model of the nonslipping and purerolling twowheeled mobile robot is concisely exposed.
The position and direction of a mobile robot is given by which consists of the coordinates x, y of the robot and the angle θ between the orientation of the robot and the Xaxis. According to the motion principle of rigid body kinematics, the motion of a mobile robot can be described using the formulas (1) and (2) below [1,2]:
where is the angular velocity of the left wheel, is the angular velocity of the right wheel, ω is the angular velocity of the robot center, is the velocity of the left wheel, is the velocity of the right wheel, v is the velocity of the robot, r is the radius of the wheel, L is the distance between the two wheels.
From equations (1) and (2), we obtain
The dynamic function of the robot is defined as [1,2]
The velocity vector U is defined as . The posture vector P is defined as , which consists of the coordinate of the robot and the angle between the orientation of the robot and the Xaxis.
By combining (3), (4) and (5), we obtain the relation of the velocity vector U and the posture vector P associated with the robot kinematics:
Equations (5) and (6) describe the kinematical model of the twowheeled mobile robot.
3 Methods
In the computational analysis of the mobile robot trajectory, some specific and fast tools of MAPLE11 soft are used. These tools are generally based on a wide spread numeric method  the Fehlberg fourthfifth order RungeKutta method  the socalled ‘rkf45’ method  with degree four interpolant [3].
A widespread and fast graphical tool is the ‘phaseportrait’ tool. The ‘phaseportrait’ is a plot builder which realizes the phaseportrait for a system of differential equations. It is a fast procedure based on specific numeric methods for approximating the solution of the studied differential system. For the present computational aim, the classical method of Euler has been chosen, in fact, the ‘forward Euler’ method.
The calling sequence has the following form:
The parameters are the following:
deqns  a list or a set of firstorder ordinary differential equations, or a single differential equation of any order;
vars  dependent variable, or a list or a set of dependent variables;
trange  range of the independent variable;
inits  a set or a list of lists; initial conditions for solution curves;
options  (optional) equations of the form keyword = value.
The default method of integration is . Other methods can be specified in the optional equations. Note that because numerical methods are used to generate plots, the output is subject to the characteristics of the numerical method in use. In particular, unusual output may occur when dealing with asymptotes of solution curves. This also means that the initial conditions of the problem must be given in a standard form, that is, the function values and all derivatives up to the differential order of the differential equation at the same point minus one.
By default, plots are produced with boxed axes. In contrast with DEplot tool, this tool does not produce a direction field constituted by field arrows. This is not possible in 3D case, but its role is taken by the ‘scene’ parameter. Below some of the basic parameters of the procedure are detailed.
The inits parameter must take the form of
inits is a list (or a set) of lists, each sublist specifying one group of initial conditions (for one solution curve).
The ‘scene’ parameter has the following form:
scene specifies the plot to be viewed. For example, indicates that the plot of x versus y versus z is to be plotted, with t implicit, while plots t versus y versus z (t explicit). This option can also be used to change the order in which to plot the variables. If vars is entered as a set, there is no default ordering; if entered as a list, the given ordering is used.
The parameter ‘stepsize’ has the following form:
and specifies the distance between mesh points to be used in generating the graph. For , the default stepsize value is . If the stepsize specified is too large, the default is used.
In this paper the ‘phaseportrait’ graphical tool is tested in order to get comparative graphical analysis with other MAPLE tools. This plot tool is appropriate for the proposed model as it produces an appropriate representation of the trajectory of the studied mobile.
4 Computational simulations of the kinematic model
4.1 The Cauchy problem associated to the model
From mathematical standpoint, equations (3) and (6) represent a system of firstorder differential equations together with its basic parameters.
The Cauchy problem associated to equation (3) is considered:
where the associated initial condition is deduced from the initial velocity of the robot, that means, from the left and right velocity pairs given by the relation (2) above.
For the simplicity of the calculus, we take in the following analysis the notation
and thus the Cauchy problem becomes
where
4.2 Computational modeling in case studies
In this section the results of computational simulations of the mobile robot trajectory in some specific conditions are presented. These conditions are in agreement with the experimental simulations recently realized in the University of Craiova in order to analyze the trajectory of this robot type.
In our experiment, we use a YSRA mobile robot which is small in size () and a black wooden rectangular playground. A host computer controls the robot by commanding the robot velocities like a radiocontrolled car. There are three different programs included in the YSRA robot; a program receives commands of the host computer, the other program counts the numbers of Encoder, and another program controls the motor’s movement by speed commands from the host computer.
The robot system consists of four parts: microcontroller, motor driver, communication, and power unit. The functions of the components are the following:
1. the communication part: receives communication commands from the host computer;
2. the power part: makes power to operate the motor and microcontroller;
3. the microcontroller part: translates the communication commands from the host computer and makes the control commands of the motor. The microcontroller is Intel 80C296SA;
4. motor part: impresses electricity to the motor as much as wanted by commands of the motor.
The experimental results obtained with the YSRA robot are presented in snapshots in Figures 1627. This experiment is filmed by a SAMSUNG Digital Color CCD camera (SDC410ND) and Figures 1627 are captured from a clip. Following the pictures in their successive order, it comes easy that the robot trajectory is indeed circular.
The mobile robot trajectory is simulated in an interval of time units, divided into five parts, representing the stages of observation and simulation cases too.
The initial parameters are as follows: the distance between the two wheels , the initial position and orientation for the first stage being , . The ‘phaseportrait’ procedure is used, where the initial conditions at the next stage depend on the velocities calculated at the precedent stage. For the velocities pair, a few trigonometric variation functions were tested [4]. For the present aim, the left and right velocities were varied in two separate situations:
We have chosen these combinations, taking into account that the basic trigonometric functions sin and cos are periodic and their variation is formed by loops. Figures 5 and 6 show the functions for the two cases.
Figure 5. Case 1 for the variation functions of the left and right velocities of the mobile.
Figure 6. Case 2 for the variation functions of the left and right velocities of the mobile.
For each of the above cases, the ‘phaseportrait’ procedure was applied step by step for the Cauchy problem (9) on the intervals: , , , , . At each step, the initial conditions, the values of the velocities pair , were calculated based on the solution of the differential system (9) at the preceding step. The ‘scene’ parameter was chosen in order to give for the mobile trajectory a ‘seen from above’ image.
First there are presented the plots for the first simulation case, in Figures 711. Each case and specific observations are labeled on the figure.
Figure 7. ,.
Figure 8. ,. The trajectory records negative values and seems to be routed several times.
Figure 9. ,. Important modifications of the trajectory.
Figure 10. ,.
Figure 11. ,. The trajectory seems to be formed by arcs.
For the second simulation case, the results are presented in Figures 1215.
Figure 12. ,.
Figure 13. ,. Important modification of the trajectory.
Figure 14. ,. Special modification in the allure and values of the trajectory.
Figure 15. ,. Significant modification comparing to the preceding cases.
In Figures 1627, the trajectory evolution for the YSRA robot is presented in snapshots. As specified above, the figures exhibit, in their successive order, a circular trajectory.
Figure 16. Frame 1 captured from clip using the variation functions of the left and right wheel velocities proposed in Case 1.
Figure 17. Frame 2 captured from clip using the variation functions of the left and right wheel velocities proposed in Case 1.
Figure 18. Frame 3 captured from clip using the variation functions of the left and right wheel velocities proposed in Case 1.
Figure 19. Frame 4 captured from clip using the variation functions of the left and right wheel velocities proposed in Case 1.
Figure 20. Frame 5 captured from clip using the variation functions of the left and right wheel velocities proposed in Case 1.
Figure 21. Frame 6 captured from clip using the variation functions of the left and right wheel velocities proposed in Case 1.
Figure 22. Frame 7 captured from clip using the variation functions of the left and right wheel velocities proposed in Case 1.
Figure 23. Frame 8 captured from clip using the variation functions of the left and right wheel velocities proposed in Case 1.
Figure 24. Frame 9 captured from clip using the variation functions of the left and right wheel velocities proposed in Case 1.
Figure 25. Frame 10 captured from clip using the variation functions of the left and right wheel velocities proposed in Case 1.
Figure 26. Frame 11 captured from clip using the variation functions of the left and right wheel velocities robot proposed in Case 1.
Figure 27. Frame 12 captured from clip using the variation functions of the left and right wheel velocities of the YSRA robot proposed in Case 1.
5 Discussions
At first sight, it must be observed that at the beginning of the analysis, when starting with the initial velocities , the mobile robot trajectory seems to be circular, and at the further stages it changes, recording a periodic trend or a ‘loop’ trajectory. The trajectory modification records itself significant modifications. Thus an important influence of the initial conditions of the Cauchy problem on the mobile trajectory is noticed.
In the first simulation case, for , the trajectory becomes thick, recording a periodicity trend. This is very interesting since there is no other case with this situation in the analysis. After that, the mobile tends to have not a complete circle. Moreover, for , a trajectory formed by line segments is recorded.
In Case 2 of simulation, it is very important to notice the special trajectory form for : a half loop with positive values, and a complete different allure for with negative values. Again the trajectory allure has a strong modification from a simulation stage to the next one. And all these in no special conditions for the time unit simulation parameter: in MAPLE, the time unit can have from the smallest values to bigger ones.
In fact, it is no rule for this important change of the trajectory, all these special events issue in a random way. Thus we can say that the Cauchy problem associated to the mobile kinematics is in fact a sensitivecase model.
6 Conclusion
The basic conclusion of the above analysis is that the Cauchy problem analyzed here brings a new approach of the studied mobile robot kinematics. The ‘phaseportrait’ graphical tool is extremely helpful because it is fast and has very flexible parameters. Note that here only a situation for the ‘scene’ parameter of the procedure, namely ‘’, was used. Other possibilities for ‘scene’ would give another perspective for the mobile trajectory analysis. We must take into account that we can assimilate the mobile in his move with a material point lying in the middle of the axis that joins the centers of the left and right driving wheels, and rotating around its ICR.
Thus, from qualitative reasoning, some basic ideas for next target have to be outlined:
• To increase the number of simulations, for more variation functions for the pair .
• To use another analytical or numerical MAPLE tools for getting more analysis data on the Cauchy problem behavior. Testing step by step the Cauchy problem was possible due to the flexible structure of the MAPLE11 graphical/computational tools [5]. At the same time, this shows that these repetitive simulations are relatively easy to perform.
It is important to notice the issues of repetitive phenomena: in the above analysis, a specific allure of the trajectory is repeated. The further changes are due to the initial conditions variations. Repetitive phenomena can be collected in order to realize a global panel of random distributed events in the kinematics of a twowheeled differential drive mobile robot and to use the statistical observations in the further analysis. Finally, the feasibility and effectiveness of the proposed computational analysis for the mobile robot trajectory is demonstrated by the experimental results.
Competing interests
The authors declare they have no competing interests.
Authors’ contributions
DC realized the basic conception of the model analysis, established the analysis steps and contributed to the interpretation of the results. Also, she realized the experimental simulations with the specified mobile robot in order to match the mathematical model. AI realized the computational simulations in MAPLE11, contributed to the final analysis of results and final manuscript arrangement.
Acknowledgements
Dedicated to Professor Hari M Srivastava.
References

Fahimi, F: The Autonomous Robots: Modeling, Path Planning, and Control, Springer, Berlin (2008)

Kim, JH, Kim, DH, Seow, KT: Soccer Robotic, Springer, Berlin (2004)

Abell, ML, Braselton, JP: Maple by Example, Elsevier, San Diego (2005)

Coman, D, Ionescu, A: Motion control of the soccer robot based on fuzzy logic. In: Simos T, Maroulis G (eds.) Topics in Fuzzy Logic, 6th International Conference on Computational Methods in Science and Engineering, Crete, Greece, September 2008. Lecture Notes in Advances Computational Science 1148. pp. 57–60. Melville, New York (2009)

Hindmarsh AC, Stepleman RS (eds.): Odepack, a Systemized Collection of ODE Solvers, NorthHolland, Amsterdam (1983)