Energy-E ﬃ cient Torque Distribution Optimization for an Omnidirectional Mobile Robot with Powered Caster Wheels

: A mobile robot with no less than two powered caster wheels (PCWs) has the ability to perform omnidirectional motions and belongs to a redundantly actuated system. Redundant actuation will bring the issue of non-uniqueness of actuating torque distribution, and inappropriate choices of torque distribution schemes will lead to unexpected large required actuating torques and extra energy consumption. This paper proposes a new torque distribution optimization approach based on a gradient projection method (GPM) for the omnidirectional mobile robot (OMR) with direct drive PCWs. It can signiﬁcantly reduce the maximal required actuating torque and the energy consumption of the system. The modular kinematic and dynamic modeling method is presented ﬁrst, which is suitable for an arbitrary number of employed PCWs, as well as their install positions in the chassis. The detailed energy consumption model of the OMR, including output energy consumption and electrical energy loss, is formulated through experimental testing. The e ﬀ ectiveness of the proposed algorithms is validated by simulation examples. Lastly, the computational e ﬃ ciency of the method is veriﬁed.


Introduction
Omnidirectional mobile robots (OMRs) are widely employed to perform tasks in narrow and congested space for their ability to instantaneously move in any direction regardless of the current poses. Among a variety of categories of OMRs, the mobile robot with powered castor wheels (PCWs) has a simple and efficient wheel design to achieve omnidirectional motions [1]. It can carry a heavy payload and is less sensitive to the ground conditions owing to the continuous contact between the wheels and the ground [2]. For fully utilizing their maneuverability and agility, OMRs are powered by onboard batteries. Energy-efficient becomes a vital performance index for OMRs as they are usually restricted by the heavy and expensive batteries [3].
It is of great significance to study the energy consumption of the mobile robots, and energy-saving strategies have been conducted by many researchers via different aspects, for example, trajectory planning [4][5][6], motion control [7,8], and mechanical design [9]. In order to study the energy-efficient strategies, a detailed and accurate energy consumption model needs to first be established for the robot system. Hou et al. [10] propose an energy consumption model that incorporates three major

Design and Modelling of the OMR
The PCWs for the OMR are designed as a direct driving type. Two direct driving motors are integrated in a PCW, as shown in Figure 1. The upper motor is responsible for the steering motion, and the hub wheel motor is responsible for the rolling motion. All the PCWs are installed in the chassis of the OMR, the control system and other devices, such as the battery and sensors, are installed inside the mobile platform. By virtue of the direct driving PCWs, there is no transmission mechanism in the system, and the efficiency and stability of motion of the OMR are improved. Moreover, the PCW is integrated with less components, and the reliability of the direct driving PCW is also better than the conventional ones. The prototype of the OMR with two direct driving PCWs is shown in Figure 2.

Kinematic Model of the OMR
The kinematic analysis of the OMR with PCWs will be presented in this section. Because the OMR moves on a horizontal surface, its motion can be described in a 2D frame. The OMR is designed in a modular pattern, which means we can choose the appropriate number of PCWs as well as their installing positions according to specific task requirements. A schematic diagram that illustrates the kinematic model of that modular designed OMR is presented in Figure 3. We assume that there are n-PCWs installed in the chassis of the OMR, and the points

Design and Modelling of the OMR
The PCWs for the OMR are designed as a direct driving type. Two direct driving motors are integrated in a PCW, as shown in Figure 1. The upper motor is responsible for the steering motion, and the hub wheel motor is responsible for the rolling motion. All the PCWs are installed in the chassis of the OMR, the control system and other devices, such as the battery and sensors, are installed inside the mobile platform. By virtue of the direct driving PCWs, there is no transmission mechanism in the system, and the efficiency and stability of motion of the OMR are improved. Moreover, the PCW is integrated with less components, and the reliability of the direct driving PCW is also better than the conventional ones. The prototype of the OMR with two direct driving PCWs is shown in Figure 2.

Kinematic Model of the OMR
The kinematic analysis of the OMR with PCWs will be presented in this section. Because the OMR moves on a horizontal surface, its motion can be described in a 2D frame. The OMR is designed in a modular pattern, which means we can choose the appropriate number of PCWs as well as their installing positions according to specific task requirements. A schematic diagram that illustrates the kinematic model of that modular designed OMR is presented in Figure 3. We assume that there are n-PCWs installed in the chassis of the OMR, and the points ( ) represent the install positions for the PCWs, which are also called support points. The support points can also be regarded as the distal ends of the PCWs. A PCW is actuated by two active joints, which are responsible for the rolling motion and steering motion, respectively. It can produce a driving wrench to the support point in the platform. The resultant forces by all PCWs will drive the OMR to perform the desired trajectory.

Kinematic Model of the OMR
The kinematic analysis of the OMR with PCWs will be presented in this section. Because the OMR moves on a horizontal surface, its motion can be described in a 2D frame. The OMR is designed in a modular pattern, which means we can choose the appropriate number of PCWs as well as their installing positions according to specific task requirements. A schematic diagram that illustrates the kinematic model of that modular designed OMR is presented in Figure 3. We assume that there are n-PCWs installed in the chassis of the OMR, and the points A i (i = 1, 2, · · · , n) represent the install positions for the PCWs, which are also called support points. The support points can also be regarded as the distal ends of the PCWs. A PCW is actuated by two active joints, which are responsible for the rolling motion and steering motion, respectively. It can produce a driving wrench  The global reference frame XYZ is used to describe the configuration of the OMR, and , , i j k    represent the three orthogonal axes of the local frame attached to the platform. As long as the number of PCWs is no less than two for singularity avoidance [20], the proposed kinematic analysis for the modular OMR is not restricted by the number of PCWs. Conventionally, the PCWs are arranged symmetrically for good dynamic performance and increasing the number of PCWs can improve the load capacity of the OMR. It should be noted that the kinematic model in this paper refers to the instantaneous kinematics, which maps the joint velocity q to the operational space velocity x . As shown in Figure  respect to the X-axis. r denotes the wheels' radius and b denotes the offset of the PCW, l is the vector for the center of the OMR to the support point, and α is the angle between l and i  . Then for the i th PCW, the linear speed at the center of the platform can be derived by Note that ω , representing the angular velocity of the OMR, can be computed by Combining (1) and (2), the kinematics of a PCW can be expressed in a single matrix form written as -sin sin cos sin cos cos sin cos 1 0 1 Equation (3) presents the forward kinematic model of wheel i and the kinematic models of the rest of the wheels can be derived in the same manner. The single castor wheel Jacobian matrix i J is always full rank as its determinant remains a nonzero constant. Thus, the inverse matrix of i J always exists. In practice, it would be difficult to measure or control i σ as it represents a passive rotational motion. This variable needs to be eliminated from the inverse kinematics of a single PCW. The The global reference frame XYZ is used to describe the configuration of the OMR, and → i , → j , → k represent the three orthogonal axes of the local frame attached to the platform. As long as the number of PCWs is no less than two for singularity avoidance [20], the proposed kinematic analysis for the modular OMR is not restricted by the number of PCWs. Conventionally, the PCWs are arranged symmetrically for good dynamic performance and increasing the number of PCWs can improve the load capacity of the OMR. It should be noted that the kinematic model in this paper refers to the instantaneous kinematics, which maps the joint velocity . q to the operational space velocity .
x. As shown in Figure 3, ρ i (i = 1, 2, · · · , n) are defined as the rolling angles of the wheels and ϕ i are defined as the steering angles; σ i are defined as the angular displacements of the ith wheels with respect to the X-axis. r denotes the wheels' radius and b denotes the offset of the PCW, l is the vector for the center of the OMR to the support point, and α is the angle between l and → i . Then for the ith PCW, the linear speed at the center of the platform can be derived by Note that ω, representing the angular velocity of the OMR, can be computed by Combining (1) and (2), the kinematics of a PCW can be expressed in a single matrix form written as Equation (3) presents the forward kinematic model of wheel i and the kinematic models of the rest of the wheels can be derived in the same manner. The single castor wheel Jacobian matrix J i is always full rank as its determinant remains a nonzero constant. Thus, the inverse matrix of J i always exists. In practice, it would be difficult to measure or control σ i as it represents a passive rotational motion. This variable needs to be eliminated from the inverse kinematics of a single PCW. The controllable joint velocity denoted by q i , including steering and rolling variables only, can be expressed as where C i is called the wheel constraint matrix [16] because it describes the constraints at each contact point between the wheel and the ground. All wheels' equations are combined to get the unified inverse kinematics, as follows: x.

Joint Space Dynamic Model of the OMR
In order to formulate the dynamic model of the entire system, the dynamic model of a single PCW should first be obtained. In Figure 4, the schematic model of the PCW is presented. A PCW consists of three main components: a wheel, a bracket, and a "virtual" link from the upper side of the bracket to the mass center of the OMR. We call the link "virtual" because it only exists conceptually. The PCWs are installed directly in the chassis of the robot. The virtual link just provides the position information of the support points from the mass center of the robot.
where i C is called the wheel constraint matrix [16] because it describes the constraints at each contact point between the wheel and the ground. All wheels' equations are combined to get the unified inverse kinematics, as follows:

Joint Space Dynamic Model of the OMR
In order to formulate the dynamic model of the entire system, the dynamic model of a single PCW should first be obtained. In Figure 4, the schematic model of the PCW is presented. A PCW consists of three main components: a wheel, a bracket, and a "virtual" link from the upper side of the bracket to the mass center of the OMR. We call the link "virtual" because it only exists conceptually. The PCWs are installed directly in the chassis of the robot. The virtual link just provides the position information of the support points from the mass center of the robot. is always orthogonal to the steering motion i ϕ  , thus they can be decoupled and calculated by where r I is the inertial moment of the wheel about its rolling axis, s I is the inertial moment of the PCW about the steering axis, and w m is the mass of the PCW. The dynamic model of a PCW in the joint space can be derived as Because i K is only in regard to i q  , / 0 The dynamic model of a single PCW can be formulated as The equation of motion for a PCW can be obtained through the Lagrangian method. The OMR can only achieve planner motion, hence there only exists kinetic energy K i . The rolling motion . ρ i is always orthogonal to the steering motion . ϕ i , thus they can be decoupled and calculated by where I r is the inertial moment of the wheel about its rolling axis, I s is the inertial moment of the PCW about the steering axis, and m w is the mass of the PCW. The dynamic model of a PCW in the joint space can be derived as d dt Because K i is only in regard to . q i , ∂K i /∂q i = 0 holds. The dynamic model of a single PCW can be formulated as M i .. where In (8), τ i represents the actuating torques and τ f i denotes the friction torques, which consist of Coulomb friction and viscous friction. P i is the matrix that maps the ith joint velocity to the support point velocity, and −F i is the external wrench applied by the platform.
The platform of the OMR can then be regarded as a manipulated object by the PCWs. The force resulted from the distal end of ith PCW is F w i . The dynamic model of the platform only can be expressed as Λ p ..
where Λ p = diag m o m o I o is the inertial term of the platform, and µ p (x, . x) is the Coriolis and centrifugal term. F w o is the resultant force applied by all the PCWs.

Operational Space Dynamic Model of the OMR
The operational space dynamic model of a PCW is expressed in the form as where C i is the operational space Coriolis and centripetal terms, and F i is the operational space wrench. The OMR is assumed to move in a plane surface, so the gravity component is neglected here.
After deriving the dynamic models of the PCWs, the augmented object model (AOM) [21] can be employed to obtain the dynamic model of the entire system. The AOM declares that the systematic dynamic model of the entire system in operational space can be obtained by adding all the operational space dynamics of the PCWs and the platform only at the operational space frame. Thus, the augmented mass matrix Λ ⊕ and Coriolis and centripetal terms µ ⊕ of the OMR are given as where Λ p and µ p are the mass matrix and the Coriolis and centripetal terms of the platform only expressed in the operational space, respectively. F is the resultant wrench applied to the OMR, expressed in the operational space frame. Finally, the operational dynamics of the OMR can be written as

Energy Consumption Model of the OMR
The intact energy consumption model of the OMR consists of three parts: the sensors, the control system, and the motion system. The motion system consumes the overwhelming majority of energy in the OMR [10]. Hence, we mainly focus on the study of the influence of torque distributions to the energy consumption. The energy consumptions of the sensor system and the control system will not be discussed here. The formulation of energy consumption model is for predicting and evaluating the Energies 2019, 12, 4417 7 of 19 amount of energy consumed by the actuators in the OMR. The energy consumption can be calculated by integrating the power consumption about time. The power consumption model usually consist of output power consumption P output and the subordinate electrical power loss P loss [22,23].
where the output power consumption is the power consumed to attain and sustain the desired motion of the OMR, and refers to the mechanical power. P output can be computed by The definition of [a] * is that, if a is larger than zero, the value of [a] * is a. If a is smaller than zero, the outcome value of the function is zero. It means that the negative output power cannot be regenerated and is dissipated as heat [24].
The electrical power loss often contains three parts: coil loss, conduction loss, and switching loss. The coil loss results from the heat dissipation as the current flows through the resistances of a circuit. It is formulated as follows: where R i is the resistance of the coil in the ith actuator. I i is the current the of the ith actuator. The other two types of electrical power loss that happen in servo-amplifiers or servo-drives, including conduction loss and switching loss. Each loss happens when the current flows in an insulated-gate bipolar transistor of the servo-amplifier [25]. The conduction loss and switching loss are often formulated according to the voltage and current information. The literature of [25] provides a scheme for calculating the conduction and switching loss only dependent on the current information. The conduction loss and the switch loss can be expressed, respectively, by: where µ a i and µ b i are the coefficients for the conduction loss. The conduction loss is formulated in a quadratic polynomial manner with respect to the current. η i is the constant coefficient of switching with respect to the current. Hence, the electrical power loss can be computed by summing the coil loss, conduction loss, and switching loss.
All kinds of losses are formulated as the functions of current. By summing the coefficients according to the power of the current, the power loss can be expressed as where λ a i = µ a i + η i , and λ b i = µ b i + R i . λ a i and λ b i are the coefficients of the power loss of the ith actuator in terms of the current.
Because we mainly investigate the relationship between the torque distributions and the energy consumption in this paper, it is desired to express the power consumption in the manner of joint torques. Fortunately, for the torque motors employed in the PCWs, their outputs are linearly proportional to the current in the circuit. Therefore, the current term I i and I 2 i in (20) can be replaced by the torque term τ i and τ 2 i with their corresponding linear coefficients. Equation (20) can be reformulated as where κ a i and κ b i are the coefficients for the electrical power loss of the absolute value of torque and the square of torque in the ith actuator, respectively. The new electrical power loss model (21) in terms of the torque represents a combination of the electrical power loss model in terms of the current and the linear proportional relationship between the current and the torque in the actuator.
After obtaining the output power consumption model of the OMR, the coefficients of the model have to be specified via experimental measurement. For accuracy, the angular velocities of the actuators are set to be constant, and photoelectric encoders on the motors can be used to measure the velocities of motors and provide feedback to the control system in real time [10]. A strain-gauge sensor [24] is employed to measure the torque of each actuator. The coefficients in (21) are acquired through the power measurement while different static loads are applied to each actuator, and the experiment is executed as shown in Figure 5. The power measurement is implemented by a dynamometer (SUGAWARA PC-MCT6). The tested actuator is connected to the dynamometer by a shaft coupling. The industrial personal computer (IPC) is responsible for controlling and sending commands to the dynamometer. The actuator is commanded to rotate in constant angular velocity with several static loads induced by the dynamometer, and the real-time voltage and current of the actuator are shown in the oscilloscope (YOKOGAWA WT1800). The output power of the actuator can be calculated by multiplying its angular velocity and the output torque applied by the dynamometer. Figure 6a,b shows the experimental result of electrical power loss of the rolling actuator and the steering actuator with respect to different loads, respectively. The markers are measurement points of the experiment. The coefficients of the power consumption models are obtained through mathematical curve fitting by MATLAB by (21), and the results are shown in Table 1.
The energy consumption model for the OMR can be calculated by Each PCW has two motors for rolling and steering motions, respectively, so the number n is 2m for the OMR with m PCWs.

Torque Distribution Schemes for the OMR
The torque distribution schemes of robots can be derived from the forward kinematic model by where τ is the actuating torques applied by the joints and J is the Jacobian matrix of the robot. In the analysis of the torque distribution for the OMR, because the PCWs are direct driven, the impact of the transmission system in the PCWs and the dynamics of the motors are neglected. We assume that the torques produced by the actuators are the same torques applied to the joints. Owing to the parallel structure and redundant actuation of the OMR, the Jacobian matrix cannot be derived  From the above, the total power consumption of the motion system can be written as The energy consumption model for the OMR can be calculated by Each PCW has two motors for rolling and steering motions, respectively, so the number n is 2m for the OMR with m PCWs.

Torque Distribution Schemes for the OMR
The torque distribution schemes of robots can be derived from the forward kinematic model by where τ is the actuating torques applied by the joints and J is the Jacobian matrix of the robot. In the analysis of the torque distribution for the OMR, because the PCWs are direct driven, the impact of the transmission system in the PCWs and the dynamics of the motors are neglected. We assume that the torques produced by the actuators are the same torques applied to the joints. Owing to the parallel structure and redundant actuation of the OMR, the Jacobian matrix cannot be derived directly through kinematics analysis. Instead, it is obtained by computing the generalized inverse of the constrained matrix C of the OMR [15]. This means that if the task space velocity of the OMR is given, the joint velocity can be directly calculated by the constraint matrix. However, the task space velocity cannot be obtained from arbitrary joint velocity because of the kinematic constraints. The constraint matrix maps the operational space velocity .
x to the joint space velocity . q by .
The Jacobian matrix of the OMR can be obtained by computing the generalized inverse matrix of the constraint matrix: Because the OMR is redundantly actuated, the constraint matrix C is not square. Hence, there exists a non-unique Jacobian matrix for the OMR, and the different choices of the generalized inverse will lead to different torque distribution schemes.

Torque Distribution Schemes Resulting from the Generalized Inverse of the Jacobian
Two reasonable determinations of torque distribution schemes proposed by the published research [17] are first presented. The kinematics of the OMR with PCWs can be expressed in an equivalent way as A .
where B is a diagonal matrix and A is a non-square matrix. The first torque distribution scheme is expressed as where B −1 A LPI means computing the left Moore-Penrose pseudo-inverse of B −1 A. Using this model, the joint torque differences are minimized in a least-squares manner. The second torque distribution scheme is expressed as where A LPI means computing the left Moore-Penrose pseudo-inverse of A. Using this model, the joint torque as well as the interaction forces between different contact points are minimized, in a least-squares manner.

Gradient Projection-Based Torque Optimization
The gradient projection method is an efficient and widely-employed optimization method in the control of redundant robots [26,27]. It fully utilizes the null space of the Jacobian matrix to optimize varieties of performance criteria, such as manipulability, obstacle avoidance, and joint limit avoidance. It is approached by projecting the steepest descent direction of the performance criterion onto the null space and finding the best solution along the projected vector. The flowchart of the proposed torque distribution optimization algorithm is shown in Figure 7. For the torque distribution of the redundant actuated OMR, as mentioned before, the joint space torque is not directly projected from the operational space force by the transpose of the Jacobian matrix, but instead of the constraint matrix C, written by Energies 2019, 12, x FOR PEER REVIEW 12 of 21 The gradient projection method is an efficient and widely-employed optimization method in the control of redundant robots [26,27]. It fully utilizes the null space of the Jacobian matrix to optimize varieties of performance criteria, such as manipulability, obstacle avoidance, and joint limit avoidance. It is approached by projecting the steepest descent direction of the performance criterion onto the null space and finding the best solution along the projected vector. The flowchart of the proposed torque distribution optimization algorithm is shown in Figure 7. For the torque distribution of the redundant actuated OMR, as mentioned before, the joint space torque is not directly projected from the operational space force by the transpose of the Jacobian matrix, but instead of the constraint matrix C , written by

Gradient Projection-Based Torque Optimization
As the 3-DOF OMR is actuated by ( ) 4 m m ≥ actuators, the joint torques can be expressed as where m v R ∈ is an arbitrary vector. On the basis of GPM, it can be replaced by the gradient of a performance criterion ( ) H τ , which yields As the 3-DOF OMR is actuated by m(m ≥ 4) actuators, the joint torques can be expressed as where v ∈ R m is an arbitrary vector. On the basis of GPM, it can be replaced by the gradient of a performance criterion H(τ), which yields where k is the searching step size and ∇H(τ) is the gradient of H(τ). Note that scalar k is positive if H(τ) needs to be maximized and is negative if H(τ) needs to be minimized. The detailed construction of H(τ) for the optimization will be presented later.
In order to implement GPM to the torque distribution of the OMR, the torque limits have to be determined first. The limit of the torque τ max and τ min is determined by two factors: the maximum and minimum produced torques of the actuator τ motormax and τ motormin ; and the maximum ground traction τ traction , which is expressed below: where N is the ground contact force of the wheel, µ c is the Coulomb friction coefficient, µ v is the viscous friction coefficient, and v c is the velocity of the wheel at the wheel-ground contact point. In practice, the maximum ground traction is smaller than the maximum produced torque, which means that, during the motions of the OMR, the actuating torques may exceed the ground traction and, consequently, slippage happens to the wheels. Employing the limit of the actuating torques and rewriting (22), it yields the following: Let τ = τ − (C T ) + F, where τ is a point in the null space of C T , then the optimization problem can be expressed as (35) Resolving the optimization problem of (35) will produce an optimal solution τ g in the null space of C T . After that, the optimal solution of τ g can be obtained by τ g = τ g + C T + F. After defining the optimization problem of the actuating torques of the OMR, the performance criterion H(τ ) need to be formulated. The boundary conditions in (35) declare that τ must be within the feasible zone, which is away from the actuating torque limit. It can be expressed as The performance criterion used to generate the feasible zone can be given as where i refers to the ith component of the corresponding torque vector. For such a performance criterion, the value of the function goes to infinity, while the actuating torque approaches either a maximal or minimal limit, and goes to minimal while the actuating torque reaches the midpoint of the feasible torque interval. This indicates that, while searching for the minimal point, the performance criterion will retain the actuating torque within the feasible torque zone and away from either the maximum or minimal limit. Hence, using this performance criterion, the maximal (or absolute value of minimal) required actuating torque of the OMR for a specific trajectory will be reduced.
In addition, it is our main purpose to take the energy consumption into account for the torque distribution of the OMR. Therefore, another term that resulted from (22) is introduced into the performance criterion, which is given as The modified performance criterion is used to optimize for a preference of better energy efficiency. It should be noted that the performance criterion function must be continuous while applying GPM, but the function in (38) is discontinuous at the point τ i = 0. This issue can be solved by using a curve fitting polynomial function to replace the original function at a small interval [−0.1, 0.1] around τ i = 0, with an acceptable sacrifice of accuracy. Now, with the objective function given in (38), the optimization problem expressed in (35) becomes a convex optimization problem because both the objective function and the inequality constraints are convex. In this situation, using the GPM can always find the global optimal solution. The GPM is executed by the following equation: where τ lg is the local optimal solution along the unit projective vector (P v ) unit , τ start is the starting point in the null space of C T , and k is the step size along (P v ) unit . On the basis of the definition of gradient, (P v ) unit is obtained by projecting the steepest descent direction of the objective function onto the null space of C T . In order to get the global optimal solution, set τ start = τ lg and repeat (39) again and again. The global optimal solution τ g will be obtained when the iteration stops at (P v ) unit = 0, indicating that the performance criterion has already come to its "lowest position". The last problem that needs to be coped with is the determination of the starting point τ start . The GPM will have a low computational cost and converge more rapidly if τ start is closer to the global optimal solution. A simple, but effective method is to project the torque distribution in (29) onto the null space of C T , and to set the result as the starting point for each motion in the entire trajectory, as the result in (29) has already been optimized in a least-squares manner. However, a more reasonable starting point in the current motion j of the planned trajectory can be obtained by projecting the optimal solution of last motion j − 1 onto the null space of the current motion j, expressed as As F and C T are continuous, τ start_ j will be very close to the new optimal solution τ g_ j if τ g_ j is the optimal solution in the last motion j − 1.

Simulation Examples and Results Analysis
The simulation is presented in this section for validating the effectiveness of the proposed optimization method by MATLAB (r-2015b version, MathWorksS, Natick, US). The OMR used for simulation is shown in Figure 8. All PCWs in the OMR are identical. The platform of the OMR is designed in squared shape, and each PCW locates at one of the four corners of the platform. The parameters of the OMR are given as follows: the mass of the platform without any PCWs is 52 kg, the side length L of the platform is 0.5 m, the radius of each wheel is 0.05 m, the offset of the castor wheel is 0.046 m, and the mass of each PCW is 5 kg. The coefficients of the electrical power loss are referred from the experimental results in Figure 6. The lower and upper bounds of the actuating torque are ±6N · m for the rolling motors and ±8N · m for the steering motors. simulation is shown in Figure 8. All PCWs in the OMR are identical. The platform of the OMR is designed in squared shape, and each PCW locates at one of the four corners of the platform. The parameters of the OMR are given as follows: the mass of the platform without any PCWs is 52 kg, the side length L of the platform is 0.5 m, the radius of each wheel is 0.05 m, the offset of the castor wheel is 0.046 m, and the mass of each PCW is 5 kg. The coefficients of the electrical power loss are referred from the experimental results in Figure 6. The lower and upper bounds of the actuating torque are 6N m ± ⋅ for the rolling motors and 8N m ± ⋅ for the steering motors.  In order to further illustrate the results of the torque distribution optimization for the OMR with PCWs, a curvilinear trajectory is proposed next. The path is shown in Figure 9; the OMR moves from the original point to the endpoint, then turns around at the endpoint and moves back to the original point. The OMR is commanded to track the path with self-rotation for taking full advantage of holonomicity, and the red arrows in Figure 9 represent the head directions (vector i  in Figure 8) of the OMR. The function of the trajectory is expressed as follows: x=4 4cos(t/20) ; (meter) y=3 3cos(t/10) ; (meter) =3 3cos(t/15) ; (rad) (41) In order to further illustrate the results of the torque distribution optimization for the OMR with PCWs, a curvilinear trajectory is proposed next. The path is shown in Figure 9; the OMR moves from the original point to the endpoint, then turns around at the endpoint and moves back to the original point. The OMR is commanded to track the path with self-rotation for taking full advantage of holonomicity, and the red arrows in Figure 9 represent the head directions (vector → i in Figure 8) of the OMR. The function of the trajectory is expressed as follows: simulation is shown in Figure 8. All PCWs in the OMR are identical. The platform of the OMR is designed in squared shape, and each PCW locates at one of the four corners of the platform. The parameters of the OMR are given as follows: the mass of the platform without any PCWs is 52 kg, the side length L of the platform is 0.5 m, the radius of each wheel is 0.05 m, the offset of the castor wheel is 0.046 m, and the mass of each PCW is 5 kg. The coefficients of the electrical power loss are referred from the experimental results in Figure 6. The lower and upper bounds of the actuating torque are 6N m ± ⋅ for the rolling motors and 8N m ± ⋅ for the steering motors.  In order to further illustrate the results of the torque distribution optimization for the OMR with PCWs, a curvilinear trajectory is proposed next. The path is shown in Figure 9; the OMR moves from the original point to the endpoint, then turns around at the endpoint and moves back to the original point. The OMR is commanded to track the path with self-rotation for taking full advantage of holonomicity, and the red arrows in Figure 9 represent the head directions (vector i  in Figure 8) of the OMR. The function of the trajectory is expressed as follows: x=4 4cos(t/20) ; (meter) y=3 3cos(t/10) ; (meter) =3 3cos(t/15) ; (rad) (41) Compared with the series manipulators, the mobile robots usually have much larger tracking errors, and one of the primary reasons for the tracking errors is the slippage between the wheels and the ground. The slippage is often difficult to predict because of the uncertainty of the ground conditions. However, because the phenomenon often occurs when the actuating torques of the wheels exceed the traction limit of the ground, reducing the maximal required actuating torques of the wheels during the performance of a specific trajectory can decrease the tendency of slippage. According to the normal indoor ground condition and the limits of the actuators, the maximal torque limits are calculated by (33) and set to be ±5.3N · m and ±4.6N · m for the rolling and steering actuators, respectively. Figure 10 shows the maximal required actuating torque among all PCWs during the motion, with both forward and backward rotation. The maximal positive actuating torque is significantly reduced by the optimization, by approximately 11% at the highest point compared with the former torque distribution Schemes 1 and 2 expressed in (28) and (29), respectively. The maximal negative actuating torque is also reduced by 6% at the highest point. Moreover, the red dotted lines represent the maximal positive and negative torque limits of the motors in the PCW. The maximal torque resulted from the former two torque distribution schemes exceed the torque limit. As it is difficult to simultaneously consider the torque limits when using the pseudoinverse method. In contrast, the optimized torque is totally below the limit. For such a trajectory, the actuators are commanded to produce less torques for the backward rotation, so that all three torques' distributions have kept the maximal negative torques within the torque limits. For other joints with less actuating torques, the differences between the optimization outcomes and the results from pseudoinverse are not as significant as the differences shown in Figure 10. However, the torques of the other joints are still decreased to a smaller extent, and the data are not listed for concision. Figure 11 shows the power consumption of the OMR following the curvilinear trajectory adopting the three kinds of torque distribution schemes. As shown in Figure 11b,c, the power consumption is reduced by using the optimized torque distribution compared with the other two schemes, basically owing to the saving in electrical power loss part. The output power consumptions for the three torque distribution schemes are the same; because the OMR is commanded to follow the same trajectory by using the three different torque distribution schemes, the output power consumption that attains and sustains the motion of the OMR will be the same for the three cases. Besides being transferred into kinetic energy, the actuating torques are partly consumed to produce internal forces. The internal forces will be counteractive to each other and do not contribute to the motion of the OMR; the portion of the actuating torques corresponding to the internal forces still consumes energy in the actuators and is transformed to different kinds of losses expressed in (27). The total energy consumption as well as energy loss are presented in Table 2. By using the proposed optimization method, the total energy consumption can be reduced by 9.3% for the given trajectory. Furthermore, the torque distribution optimization has to be computationally efficient if it is applied in real-time dynamic control. Because the optimization algorithm is iteration-based, the computational time is determined by the number of iterations per step. The step size of the continuous motion is discretized in every 0.1 second, and the average numbers of iterations per motion step are shown in Table 3. In order to reveal the superiority of the proposed method for determining the starting point expressed in (40), the results of the number of iterations from two other choices of starting point are also presented as a contrast. The first choice is an arbitrary simple point, and the second choice is the pseudo-inverse solution resulting from (29). The computational efficiency of the proposed method for determining the starting point is apparently approved. Compared with the series manipulators, the mobile robots usually have much larger tracking errors, and one of the primary reasons for the tracking errors is the slippage between the wheels and the ground. The slippage is often difficult to predict because of the uncertainty of the ground conditions. However, because the phenomenon often occurs when the actuating torques of the wheels exceed the traction limit of the ground, reducing the maximal required actuating torques of the wheels during the performance of a specific trajectory can decrease the tendency of slippage. According to the normal indoor ground condition and the limits of the actuators, the maximal torque limits are calculated by (33) and set to be 5.3N m ± ⋅ and 4.6N m ± ⋅ for the rolling and steering actuators, respectively. Figure 10 shows the maximal required actuating torque among all PCWs during the motion, with both forward and backward rotation. The maximal positive actuating torque is significantly reduced by the optimization, by approximately 11% at the highest point compared with the former torque distribution schemes 1 and 2 expressed in (28) and (29), respectively. The maximal negative actuating torque is also reduced by 6% at the highest point. Moreover, the red dotted lines represent the maximal positive and negative torque limits of the motors in the PCW. The maximal torque resulted from the former two torque distribution schemes exceed the torque limit. As it is difficult to simultaneously consider the torque limits when using the pseudoinverse method. In contrast, the optimized torque is totally below the limit. For such a trajectory, the actuators are commanded to produce less torques for the backward rotation, so that all three torques' distributions have kept the maximal negative torques within the torque limits. For other joints with less actuating torques, the differences between the optimization outcomes and the results from pseudoinverse are not as significant as the differences shown in Figure 10. However, the torques of the other joints are still decreased to a smaller extent, and the data are not listed for concision. Figure 11 shows the power

Conclusions and Future Work
This paper presents the optimization of the torque distribution for the redundantly actuated OMR with PCWs based on GPM. The kinematic and dynamic models of the OMR are firstly formulated. In this paper, although the model used in simulation is an OMR with four PCWs, the modelling approach is in a modular way, which is suitable for the OMR with an arbitrary number of PCWs, as well as their mounting positions in the chassis. The detailed energy consumption model of the OMR including output energy consumption and electrical energy loss is then formulated, the the coefficients of the model are obtained by experimental test. Besides the energy consumption, the maximal required actuating torque to perform a specific trajectory is also non-negligible, as slippage often occurs when the actuating torques of the wheels exceed the ground traction limits. A performance criterion combined with torque limits and energy consumption is proposed in the GPM. Compared with two torque distribution schemes, which result from the Moore-Penrose pseudo-inverse method, the optimized torque distribution can significantly decrease the maximal required actuating torque as well as the total energy consumption, by 11% and 9.3%, respectively. Finally, the computational efficiency of the proposed optimization algorithm is validated. In further research, we intend to extend the torque distribution optimization method to an omnidirectional mobile manipulator, which means a serial manipulator is mounted on the OMR. While performing a task, the torque limit for each actuator of the PCW may change in real-time because the ground contact force of each PCW is related to the current configuration and motion of the manipulator. Moreover, the trajectory planning will become more complicated for that case, the motion of the OMR has to be decoupled from the motion of the entire system, and a multi-objective optimization method is expected to cope with the issues.