Finding the Quickest Straight-Line Trajectory for a Three-Wheeled Omnidirectional Robot under Input Voltage Constraints

1 1st R&D Institute-3, Agency for Defense Development, 111 Sunam-dong, Yuseong-gu, Daejeon 305-600, Republic of Korea 2Department of Computer Science and Engineering, Kwangwoon University, 20 Kwangwoon-ro, Nowon-gu, Seoul 139-701, Republic of Korea 3Department of Computer Engineering, College of Information Technology, Gachon University, 1342 Seongnamdaero, Sujeong-gu, Seongnam-si, Gyeonggi-do 461-701, Republic of Korea


Introduction
Currently used in many applications, mobile robots can be characterized by the motions of which they are capable.Unlike nonholonomic differential drive mobile robots [1][2][3], holonomic mobile robots such as three-wheeled omnidirectional mobile robots (TOMRs) can perform independent translational and rotational motion from any starting configuration [4].
The dynamics of omnidirectional mobile robots have been widely studied, as a prerequisite for effective controller design.Dynamic models of the orthogonal type of omnidirectional mobile robot were developed by [5,6].A model in which simplified geometric relationships determine the robot's position has since been introduced [7].The effect of slipping is taken into account in a recent model [8], a simple model of the friction between the robot's wheels and the floor.The general dynamic characteristics of -wheeled omnidirectional mobile robots with control redundancy have also been studied [9].A fuller model of TOMR dynamics, including motor parameters has recently been introduced in the global coordinate frame [10].
In this paper, we address the time-optimality of TOMR motions.However, it is worth remembering that other characteristics of a motion, such as energy consumption, safety, simplicity of planning, and accuracy, may also be of importance.Nevertheless Balkcom et al. [11] note that timeoptimality is a fundamental characteristic of a robot.
Most research on time-optimality of robot has focused on velocity and acceleration constraints [12,13].The rapidity of a motion may also be limited by torque or curvature constraints [14,15].However, since most mobile robots are battery powered, constraints on voltage may be the most crucial [16].
Many researchers have tried to minimize the motion times of two-wheeled differential-drive mobile robots [12][13][14][15][16][17][18][19].Recently, a few researchers have made progress with the same problem for TOMRs.The near-optimal strategy described by Kalmár-Nagy et al. [7] is intended to be computationally efficient, making it suitable for real-time applications; however, the input voltage constraint is not fully considered, and the underlying dynamic model omits the Coriolis terms.A different approach to optimizing TOMR times, using a genetic algorithm [20] and nonlinear programming, was suggested by Fu et al. [21]; however, it is computationally intensive, and its numerical nature does not capture the general properties of TOMR motion.Balkcom et al. [11] classified trajectories designed to minimize motion times in terms of spin, circular arc, and tangential motion.They proposed fast trajectories based on "spin, " "roll, " "shuffle, " and "tangent" motions.They did not determine which trajectory is optimal for specific start and end configurations.
The quickest trajectory for a TOMR is not usually a straight line.However, in practice, a path is often expressed as a combination of simple paths, such as straight lines or clothoid arcs.Thus, it is worthwhile to look at straight-line trajectories.
In this paper, using the full dynamics of a TOMR under input voltage constraints in global coordinates, we show how to find the quickest straight-line trajectory.We also consider motions with and without rotation of the robot.When rotation is allowed, we do not specify the final heading of the TOMR.If we were to specify its final heading, we would be faced with a multiobjective optimization: we leave that for future work.We address the problem, with or without rotation, using conventional optimization theory.
The remainder of this paper is organized as follows: in Section 2, we review the dynamics of a TOMR based on the analysis of [10].In Sections 3 and 4, respectively, we formulate and analyze the problem of minimizing the time required for a straight-line motion without and with rotation.Simulation results are presented in Section 5 and we draw some conclusions in Section 6.

Dynamics of Three-Wheeled Omnidirectional Mobile Robots
The TOMR considered in this paper consists of a cylindrical body mounted on three omnidirectional wheels, each driven by its own DC electric motor, located at equal intervals.We assume that the mass of the robot around the circumference of the robot's bone is uniformly distributed about its center and that the wheels are of the (orthogonal) type [4], meaning that each wheel can roll freely perpendicular to the direction in which it is driven by its axle.This idealized TOMR is depicted schematically in Figure 1, where  is the mass of the robot,  is the inertia of robot rotation,  is the radius of wheel,  is the distance between the center of robot and one wheel, and   is the force acting on the tangential direction in each wheel [10].
The location of the robot in the global coordinates is given by x = [,,]  , and x m = [  ,   , ]  represents its location in the moving robot coordinate frame, where   is perpendicular to the rotational direction of the first wheel.The relationship between x and x m is where We use the TOMR dynamics as suggested in [10], where three wheels are actuated by control input to motors u = [ 1 ,  2 ,  3 ]  .Here each component of u is normalized by input voltage.The dynamic equations for TOMR including the Coriolis terms in global coordinate frame are represented as follows: The coordinated inputs to the dynamic equations [  ,   ,   ]  are represented as Here  and  are, respectively, the approximate decay constants of linear and angular TOMR velocity.The constant parameter ℎ is specified by motor specification, supply voltage, wheel radius, and floor.When there is no change in the heading of the TOMR, these become first-order differential equations, and the Coriolis terms disappear.The dynamics (3) with Coriolis terms are more accurate than those suggested by Kalmár-Nagy et al. [7].

Quickest Straight-Line Trajectory for a TOMR with Input Voltage Constraints, without Rotation
We now find the shortest time for a straight-line motion by a TOMR which is not allowed to rotate.For simplicity, we consider the motion to take place along a straight line in the positive  direction, and therefore we can set ẏ = 0 and φ = 0.
For convenience, let the state vector be z = [, ẋ ]  .From the corresponding costate vector  = [  ,   ]  , we can construct the Hamiltonian Since − λ = /z, the costate terms are Unfortunately, the structure of the Hamiltonian does not allow it to satisfy the stationary condition, so we must pursue the minimum-time solution directly.We recognize Problem 1 as a typical one-dimensional bang-bang control problem [22], and, thus, using (7) and (13), the Hamiltonian can be rewritten as From ( 4)-( 6) and ( 9)- (11), we can find a shortest-time solution to (14), provided that the optimal input vector is at least 1-extreme.That means that the input to a system is in the form of bang-bang control signals with more than one extreme value.Table 1 shows which component of the optimal battery input vector is extreme, depending on the heading  0 .
The magnitude of the optimal vector  *  depending on the initial heading of the TOMR.
We can now determine the general form of the optimal input  *  that will produce the quickest transit.Looking first at  0 ∈ (0, /3), we can determine the nonextreme components of the battery input vector when  3 = 1:  1 = − sin( 0 )/sin( 0 + /3) and  2 = sin( 0 − /3)/sin( 0 + /3).As a result,  *  becomes 1.5/sin( 0 + /3).If we similarly examine the other ranges of  0 in Table 1, we can construct a generalized expression for  *  for any  0 , as follows: where ⌊⌋ is the floor function.The translational magnitude of  *  is plotted with respect to  0 in Figure 2. We see that the optimal translational magnitude of  *  varies with the heading: for |  | ≤ 1, we can assert 1.5 ≤  *  ( 0 ) ≤ √ 3. We can now say that, when a TOMR moves in a straight line without rotating, the time required for the quickest transit depends on the initial heading.This time is the shortest when the initial heading is 0, ±/3, ±2/3, or .
The input constraints (9) constitute a cube in   space, which we will call the u cube, as shown in Figure 3.Here we define terminology for later use as follows.
Definition 1.In a straight-line motion of a TOMR without rotation, six angles relate to the direction of the force on the robot to the maximum translational velocity.These angles are called the maximum translational angles (MTAs).We call the input control vectors corresponding to the MTAs the maximum translational input points (MTIPs).
The MTAs correspond to the maximum values of  *  :  0 = /3,  = 0, 1, . . ., 5. In particular, the optimal input voltage vector u * at an MTA is one of the 2-extreme MTIPs: 3 shows the six MTIPs corresponding to MTAs in the u cube.
Finally, if we define the switching time of the optimal input vector to be   , we can obtain optimal input profile u * () and the velocity profile ẋ * () from standard optimization theory applied to bang-bang systems: Furthermore, using the boundary conditions of (8), we can derive explicit expressions for the minimum motion time   and the switching time   : where The derivation of ( 17) and ( 18) is explained in the appendix.

Quickest Straight-Line Trajectory for a TOMR with Input Voltage Constraints, Allowing Rotation
We now find the quickest straight-line trajectory of a TOMR that can rotate during the motion.We consider an arbitrary initial heading, and there is no constraint on the final heading.The Coriolis terms become more important if the TOMR rotates.Again we consider a straight-line motion in the positive  direction, and so ẏ = 0. From the second equation in (3), we can formulate the constraint which is required to compensate for the Coriolis term φ ẋ in order to maintain () ≡ 0. Hence we can define the quickest trajectory problem with rotation as follows.
Problem 2. Given   (> 0), find the motion that takes the least time (  ), subject to the following constraints: dynamic constraint: configuration constraint: input constraints: (2 The notation "−" in (22) indicates that we do not care about the heading of the TOMR or its derivative at the end of the trajectory.
As we did in the previous section, we use the maximum principle in our analysis.If the corresponding state vector is z = [, , ẋ , φ ]  and the costate vector is  = [  ,   ,   ,   ]  , then the Hamiltonian is formulated as follows: From − λ = /z and ( 21), the costate functions can be expressed as follows: Here  is a Lagrangian constant which maintains the equality constraint in (24).Note that     in (27),  φ in (28), and  ẋ in (29) all stem from the Coriolis terms.Using ( 21), (25), and ( 26)-(29), we can rewrite the Hamiltonian as follows: where To minimize the Hamiltonian at any time in the motion, we can choose an input vector that minimizes     or     .In other words, all optimal input vectors must maximize the absolute value of   or   .We will refer to these two classes of input vector as max   and max   .From (30), we can conclude that max   or max   should be selected as an input vector for time-optimality of the TOMR.When rotation was not permitted, the optimum input vector was always of class max   but, for Problem 2, vectors of class max   and max   are both candidates.

Geometric Analysis of the Quickest Trajectory with an
Input Vector of Class max   .From the geometric properties of the input constraints, we find that all optimal candidates must be 2-extreme.The constraint imposed by   in ( 5) and ( 24) forces solutions to lie on a plane through the u cube, with the equation which depends on the heading .When the heading angle is between −60 ∘ and 0 ∘ , the   = 0 planes corresponding to several different headings are shown in Figure 4(a).Because  1 ,  2 , and  3 are orthogonal, any input vector to correspond to the quickest trajectory must be located where an edge of the u cube intersects the cutting plane; thus such a vertex must have at least two extreme components as white starts in Figure 4(a).
We will now select one of the two input vectors that are specified by (30) as candidates to correspond to the quickest trajectory.We first consider the vector of class max   .For any heading angle , an input vector that maximizes (4) on a cutting plane generated by the   constraint of   determines, by (32), a unique point in the u cube corresponding to that input vector (see Table 2).
At this stage we observe an important property of the optimal input vector.The sign of   in an optimal input vector of class max   depends on the initial heading of the TOMR.For example, suppose the  plane (shown in light gray in Figure 4(b)) is   = 0 (equivalent to  1 +  2 +  3 = 0).This contains six MTIPs.If the initial heading of the TOMR is between 0 and /3, then the sign of   of the optimal input vector of class max   is negative, because the extremal point is in the region  1 + 2 + 3 < 0 (the volume enclosed by darkgray planes in Figure 4(b)).Conversely, if the initial heading is between −/3 and 0, then the extremal point lies in the region  1 +  2 +  3 > 0, so the sign of   is positive.This geometric analysis leads to the following.Property 1.Consider a TOMR moving in a straight line, with allowed rotation.If we select a candidate optimal input vector of class max   , that vector is 2-extreme and the sign of its rotational element   is fixed throughout the motion.
Additionally, when the optimal candidate is of class max   , we can expect the heading of the TOMR to tend to 0, 2/3, or 4/3.For example, when the initial TOMR heading is between −/3 and 0, the sign of the optimal input vector's rotational element is always positive.Thus the heading of the TOMR approaches 0 ∘ , and then the extremal input vector moves along the edge of the u cube driving the motion, as shown by the white arrow in Figure 4(a).We introduce Definition 2 and Theorem 3 to explain this phenomenon.
Definition 2. Three of the six MTAs, 0, 2/3, and 4/3, are called stable maximum translational angles (SMTAs).The other three MTAs are called unstable maximum translational angles (UMTAs).Then the input vectors corresponding to SMTAs and UMTAs are called SMTIPs and UMTIPs, respectively.These vectors are shown in Figure 4(b).
Theorem 3. Suppose a TOMR moves in a straight line, and rotation is allowed.For all initial headings other than the MTAs, if the optimal input vector is of class max   , then the TOMR heading ultimately converges to the nearest SMTA.Furthermore, the sign of the first derivative of the heading remains the same throughout the motion.
Proof.Considering the symmetry of a TOMR, it is only necessary to consider initial headings in the range −/3 <  0 < /3.We also continue to assume that the TOMR moves in the positive  direction.From Property 1, the optimal input vector of class max   converges to [ 1 , −1, 1]  .
From ( 4)-( 6), (24), and the dynamics of φ , we can formulate the following two equations for the   and   elements of [ 1 , −1, 1]  : If we now introduce a positive definite function (, φ ), of the form then we can use (33) to find the derivative of (, φ ) with respect to time, which is Since ẋ is positive over the entire motion, V(, φ ) < − φ 2 .Near SMTA Notation:   , φ  , and ẋ  are, respectively, the heading, its first derivative, and the linear velocity at a time .
Thus, for any initial heading and any first derivative of that heading, (, φ ) converges to an invariant set [23]; that is, {(, φ ) | φ () ≡ 0}.From the dynamics of φ , if φ = 0 then   = 0, and thus  = 0. Hence the heading of the TOMR converges to 0 ∘ .Furthermore, from the properties of an invariant set, the sign of φ is invariant.
To summarize, if we select a candidate optimal input vector of class max   , then the wheel of the TOMR closest to the direction of travel gradually becomes aligned to the straight-line path, and the rotational component of the control vector has a fixed sign.This effect is illustrated in Figure 5.For example, if the TOMR moves to the right upper direction of the direction of travel, the sign of the rotational element of the input vector is positive because the first wheel is nearest to the direction of travel.During the motion, the heading of the TOMR converges monotonically to the direction of travel, as shown by the broad arrow in Figure 5.
Finally, we need to compare whether or not using a 2extreme input vector corresponding to max   is better than disallowing rotation.Since the optimal input vector without rotation is only 1-extreme, we would expect the max   route to produce a quicker transit.Unfortunately, the Coriolis terms in TOMR dynamics make it difficult to determine which method produces a quicker transit.Nevertheless, we can compare each method based on the kinematics (or the dynamics without the Coriolis terms) of a TOMR; after that, we will compare the approaches through simulations.Property 2. Suppose a TOMR moves in a straight line with no constraints on its final heading.Then, a trajectory with allowed rotation using an optimal input vector of class max   is not slower than a trajectory without rotation, if we consider the kinematics, or the dynamics without the Coriolis terms.
Proof.We only need to show that the magnitude of   in a 2-extreme input vector is no smaller than that in a 1-extreme vector over any time interval.Without loss of generality, we can assume the range of initial headings to 0 ≤  0 ≤ /3 and assume that the TOMR moves in the positive  direction.
From (15), the translational magnitude  *  of the optimal input vector without rotation can be written as follows: When rotation is allowed and the input vector is of class max   , we can express   based on kinematics (or dynamics without the Coriolis terms) as follows: From ( 36) and (37), we see that the magnitude of   in the 2-extreme control vector is never smaller than that in the 1-extreme vector.Thus we can state that a quicker trajectory is obtained by allowing rotation and using an optimal input vector of class max   than by disallowing rotation.

Combination of Candidate Optimal Input Vectors of
Classes max   and max   .So far we have not considered an input vector of class max   as the voltage input to the TOMR.This is because the Hamiltonian provides no explicit indication whether to switch candidates, if the dynamics being represented are complicated.It is impossible to use only max   to produce the quickest trajectory, so we now consider combining input vectors of classes max   and max   to produce an optimal input vector.We pursue an insight into switching candidates by looking at the kinematics of a TOMR.Theorem 4. Suppose, as before, that a TOMR moves in a straight line, with allowed rotation, and without any constraints on its final heading.The best input vector obtainable from a purely kinematic model is of class max   , for all initial headings except the UMTAs.
Proof.Again we only need to consider initial headings in the range 0 <  0 < /3, and we again assume that the TOMR moves in the positive  direction.Then, by Theorem 3, the heading of the TOMR converges to the nearest SMTA if we use an optimal input vector of class max   .If the TOMR heading decreases, then the motion takes longer.Thus, we only consider input vector of class max   , so that the TOMR heading increases.
Let Δ be the increase in the heading during the time interval  Δ .To compare a loss associated with max   with a new gain associated with max   , suppose that the optimal candidate input vector of class max   is used until the heading of the TOMR reaches  0 + Δ, and then the TOMR uses the optimal candidate corresponding to max   .
Initially, suppose the gain associated with max   and the loss associated with max   are Δ + (Δ) and Δ − (Δ), respectively.From the kinematics of a TOMR, we can obtain expressions for the elements   and   in an input vector of class max   as follows: Then, by a simple time-slicing method, we can easily show that the maximum value of Δ + (Δ) is achieved when  0 +Δ is /6.Now if we formulate expressions for Δ + (Δ) and Δ − (Δ) using (38), we obtain From these inequalities, we deduce that Δ + (Δ)/Δ − (Δ) < 1.Hence, based on TOMR kinematics, the optimal candidate input vector of class max   is never better than that of class max   .
Some investigators have doubted whether a straight-line motion of a TOMR is ever quicker if rotation is allowed.Our claim to the contrary is supported by the following corollary, the proof of which follows from Theorem 4.

Mathematical Problems in Engineering
Corollary 5. Assume that a TOMR moves in a straight line and rotation is allowed, with no constraints on its final heading.The quickest trajectory derived from TOMR kinematics occurs when the orientation of a straight-line path does not cross any UMTA.
Finally, we consider the case when we obtain a control vector corresponding to max   , when the initial heading is one of the UMTAs.If the heading of the TOMR is slightly perturbed from a UMTA, a quicker trajectory is available by increasing   .The following property describes how the quickest trajectory can be obtained when the initial heading is a UMTA.It can easily be proved from Property 2 and Theorem 4.
Property 3. Suppose that a TOMR, which is allowed to rotate, moves in a straight line without any constraint on its final heading.When the initial heading is one of the UMTAs, the quickest trajectory is achieved if the optimal input vector of class max   continues to be used after the TOMR heading is perturbed by max   .
In fact, the straight-line movement of a TOMR that maintains its heading along an MTA can be regarded as a singular motion because one component of the control vector is identically zero during a nontrivial interval [24].Thus, Property 3 suggests that a TOMR should avoid singular motion at a UMTA to obtain the best solution to Problem 2 using a straight-line path.Nevertheless, a singular motion can form a useful element of a combined trajectory [11].
Thus far we have investigated several properties of the quickest straight-line motion of a TOMR when rotation is allowed.We cannot find an analytic minimum-time solution to Problem 2. However, combining the results obtained in Section 4, we can construct the optimal input vector u * at any time, as shown in Table 2.For example, if  0 ∈ (−/3, 0), the optimal vectors for Problem 2 are shown as the stars on a u cube in Figure 4(a), and the trace of the optimal vectors is shown as the white-hollow arrow in Figure 4(a).

Simulations
We confirmed our theoretical results by simulations.In particular, we show that Property 2 and Theorem 4, while based only on the kinematics of a TOMR, show close agreement with the simulated dynamics.In all the simulations the TOMR travels over a straight path of 5 meters.The resolution of the simulations is 1 millisecond, and the simulation terminates when the error falls below 0.1%.Values of the TOMR parameters are listed in Table 3 [10].
The first simulation compared cases with and without rotation.Exploiting symmetry as before, the initial headings of the TOMR are chosen in the range −60 ∘ to 60 ∘ .Figure 6 shows that allowing rotation always produces a faster transit, for any initial heading except the SMTAs, which confirms Property 2. When the initial heading is ±30 ∘ , the travelling time without rotation is 14.4% longer than it is if rotation is permitted.We also see that the minimum transit time without   rotation is achieved when the TOMR heading is aligned to an MTA, and the optimal input vector is 2-extreme, as predicted by our analysis in Section 3.
When rotation is permitted, the quickest transit is also achieved when the initial TOMR heading is close to a UMTA.Hence we see a deviation from the results obtained from kinematics and Theorem 4. Theorem 4 suggests that the transit time will decrease as the initial heading deviates from the nearest SMTA to the direction of travel.But the simulation of dynamics shows that transit time does not always decrease proportionally to the difference between the initial heading and the nearest SMTA (see Case 2 in Figure 6).We attribute this difference to the Coriolis terms in the TOMR dynamics.
The second simulation was intended to verify the optimal input and velocity profile obtained in Section 4. Figure 7 shows the profile of the optimal voltages when rotation is allowed and the initial heading is −50 ∘ .As the identity of the extreme component changes, the input vector is instantaneously 3-extreme and then the nonextreme component gradually converges to zero.During the transit, the TOMR heading converges to the nearest SMTA (i.e.,   ≅ 0).The nonextreme component of the input vector produces rotation of the TOMR, and so this result conforms to Property 1 and Theorem 3 in Section 4.
In contradiction of Property 2, based on kinematics, we find that the magnitude of the -coordinated input can be smaller than √ 3, due to the influence of the Coriolis Linear velocity without rotation (case 1) Linear velocity with allowed rotation (case 2) Angular velocity with allowed rotation (case 2) terms in the dynamics.However, we know that the Coriolis terms of TOMR dynamics do not contradict the overall properties obtained from kinematics. Figure 8 shows that the translational (linear) velocity profiles of TOMR are directly compared in the cases of with and without rotation.In addition, the angular velocity profile of TOMR with rotation is presented as expected in Theorem 3. Figures 7 and 8 show that allowing rotation produces more extreme components, resulting in a quicker transit.The third simulation is intended to evaluate the validity of Theorem 4, based on TOMR dynamics.Figure 9 shows that  when initial heading of TOMR is not a UMTA, the fastest transit for TOMR is achieved by using max   only.The offset angle of the horizontal axis in Figure 9 means the angle changed by max   prior to the translational movement of the TOMR.The sign of the offset angle signifies a change in the initial heading towards (−) or away from (+) the nearest UMTA/SMTA.The quickest transit is not achieved by using an input vector of class max   for any initial heading other than a UMTA.Furthermore, the circle and the decreasing edge of the 60 ∘ line in Figure 9 correspond, respectively, to the fastest transit when there is no perturbation (singular motion at a UMTA) and perturbation; these results support Property 3.

Conclusions
We studied the problem of generating the quickest straightline trajectory for a TOMR with input voltage constraints.We believe this is the first attempt to solve this problem analytically.
We used dynamic equations with Coriolis terms for the motion of a TOMR, which has practical limits to motor input voltage.We then formulated the problem of finding the quickest straight-line trajectory for a TOMR which is not allowed to rotate.We found that the quickest transit occurs when the input vector to the dynamics is at least 1extreme.We show how to construct voltage input and velocity profiles for this optimal trajectory.We also discovered that the minimum transit time depends on the TOMR's initial heading and verified this observation through simulations.
A TOMR can rotate while it moves.We therefore considered the problem of finding the quickest straight-line trajectory, if the TOMR is allowed to rotate.However, the complexity of the dynamics with the Coriolis terms prevents us to explicit solution to this problem.Nevertheless, we were able to determine that the shortest transit occurs when the input vector of the TOMR is at least 2-extreme.In effect, input vectors of class max   or max   are required to achieve conditions for the quickest trajectory.From the conditions which attach to this trajectory, we derive several important properties of a TOMR.
(i) When the optimal vector is of class max   , the nearest wheel to the direction of travel gradually becomes aligned to a straight-line path and the rotational element of the input vector is of fixed sign.
(ii) When a vector of class max   is 2-extreme, the quickest transit is not slower than the quickest transit without rotation.
(iii) For any initial heading other than the UMTAs, the quickest transit is achieved by an input vector of class max   .
(iv) If the initial heading is a UMTA, the fastest transit is achieved by an input vector of class max   after a vector of class max   has been used as a perturbation.
We verified and illustrated these properties by means of several simulations.We do not think that a straight-line trajectory is usually the fastest way to get a TOMR from one position to another, but we see straight-line motions as necessary components of a path-planning scheme.We also recognize that there may be constraints on a TOMR's final heading after a motion.We expect that the solution to this problem will require switching between input vectors of classes max   and max   .Hence, in future work, we intend to study the quickest straight-line trajectory of the TOMR with two objectives.Our research can be further extended to nonstraight paths, such as a clothoid arc or the combination of clothoid arcs and straight lines.

Figure 4 :
Figure 4: Variation of the optimal input vector with heading,   realized geometrically in terms of the u cube: (a) five cutting planes corresponding to variants   constraints, and the trace of a 2-extreme input vector (white arrow), when the heading angle is between −60 ∘ and 0 ∘ .(b)  plane corresponding to   = 0;   > 0 in the region containing u = [1, 1, 1]  ; and   < 0 in the region containing u = [−1, −1, −1]  .
of the rotational element of the input vector

Figure 5 :
Figure 5: TOMR motion in a straight line related to the sign of the rotational element of the candidate control vector max   .

2 :
with rotation (free of final heading)

Figure 6 :
Figure 6: Variation in transit time with respect to the initial heading of TOMR.

Figure 7 :
Figure 7: Component voltages of the control vector when  0 = −50 ∘ and rotation is allowed.
Minimum-time trend of TOMR by the combination of optimal candidates Transit time (s) Offset angle by max  (deg) Initial heading: 25 ∘ Initial heading: 30 ∘ Initial heading: 35 ∘ Initial heading: 60 ∘

Figure 9 :
Figure 9: Shortest transit times obtained by combining input vectors of class max   and max   .When the offset angle is 0 ∘and the initial heading is not a UMTA, max   alone is used; when the offset angle is 0 ∘ and the initial heading angle is a UMTA, perturbation is used.The circle corresponds to no perturbation at a UMTA.

Table 1 :
Extreme input related to the heading  0 .

Table 2 :
Optimal control vectors for Problem 2.

Table 3 :
The parameters of TOMR.