{'fraction_non_alphanumeric': 0.030955425376906837, 'fraction_numerical': 0.007344851170121026, 'mean_word_length': 4.694209278699628, 'pattern_counts': {'":': 0, '<': 2, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 6, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
An Analysis of Power Consumption of Fluid-Driven Robotic Arms Using Isotropy Index: A Proof-of-Concept Simulation-Based Study

The manipulability of a robotic arm may be defined based on ease of motion in different directions or ease of applying force/torque. In this study, we use manipulability measures to investigate how the kinematics of robots can be employed to calculate the optimal power required to drive the actuation systems of their arms. We hypothesize that the isotropy measure is related to the power consumption of the robotic arm. In addition to theoretical aspects, we consider practical applications that can minimize power consumption in robotic systems. Since the method is simple to implement and has no assumption on the robot’s work environment or dependence on information on the main power supply, manipulability measures can be used as a tool to predict the power consumption of robotic manipulators.


Introduction
Measuring power consumption in robotic systems creates advantages such as increasing the working time of robots without changing the power supplied [1]. Power consumption measures can inform designers of how robotic arms' weight, size and configuration impact maneuverability and energy efficiency. This is particularly important for robotic systems situated remotely. More efficient energy use has positive environmental impacts, such as decreasing greenhouse gas emissions and slowing the growing demand for energy, especially in emerging economies [2].
To optimize power consumption in robotic systems, it is of the utmost importance to identify quantifiable tools which control and minimize the amount of power required. For instance, if the configuration of a stationary robot is correlated to the power consumed, adding a mobile base to the stationary manipulator would help in altering the configuration. In practice, relocating the base in stationary manipulators is not practical. Redundant stationary robots, however, can be positioned so their posture is optimal for a given task using the path-planning algorithm [3].
In this study, we define quantitative tools related to the energy required to run a system. These tools, the manipulability ellipsoids, help us to visualize how the manipulator configuration of a robot can contribute to its task execution. They visualize the directions and particular properties of arm end effectors [4]. A task is executable if the vectors stay within their ellipsoids.
Manipulability velocity ellipsoids (MVE) and manipulability force ellipsoids (MFE) are very well-known and are valuable tools to analyze, design and control robotic manipulators [5]. The shape of MVE and MFE reflect the constraints to which the robot will be subjected, such as moving or applying force in a specific direction [6]. The ellipsoids do not provide exact numerical information about maximum velocity or force values but suggest directions in which the robot can access or apply force more efficiently [7].
Manipulability measures play significant roles in quantitatively evaluating the efficiency of robotic systems. The manipulability measures indicate the ability to position and orient the end effectors or to apply force or torque. We use Yoshikawa's measure of manipulability [5], condition number [5], isotropy index [8] and the eccentricity measure of manipulability [9]. In this study, we normalize the manipulability indices for ease of comparison.
The manipulability measures are highly related to the structure and configuration of the robot [10]. For example, when the manipulator is in a singular configuration, the manipulability measures of MVE in some directions become extremely poor, while in others they are very good. Conversely, the robot may be placed in a configuration in which the measures of MVE are identical in all directions. Dexterity analysis is key when designing and evaluating the performance of robotic manipulators. Researchers have used the concept of manipulability to analyze the motion of robotic manipulators for different purposes, often to avoid singular configurations [11][12][13][14].
Although the manipulability ellipsoids are helpful tools in motion analysis of robotic manipulators, they do not offer sufficient information on how fast or slow the manipulators may move along the arbitrary direction or if they suffer from a physical inconsistency when position and orientation information are combined into a single scalar performance parameter [15]. The power manipulability ellipsoid was based on a new parameter that is independent of the selection of mechanisms' physical unit to overcome the drawback of the physical inconsistency [16]. Specifically, Mansouri et al. [16] presented the concept of a power manipulability ellipsoid: a homogenous tensor defined in six dimensions that includes both translational and rotational components (only motion components of the robotic manipulator). This concept uses a hybrid presentation of the manipulability velocity and force ellipsoids to introduce a vectorial presentation of the power consumed by a robotic manipulator.
The power manipulability ellipsoid had the same drawback as the manipulability ellipsoid and required further modifications to enable us to extract quantitative information on the power and/or kinematic measures. The concept was improved by the authors [17] by forming a quaternion formulation of the power manipulability, previously introduced as a new homogeneous performance index of robot manipulators [14].
In this paper, we utilize the concept of MVE and MFE to define the manipulability power ellipsoid (MPE) to show whether the power consumption and manipulability of a fluid-driven robotic manipulator could be theoretically correlated. We hypothesized that when the isotropy measure of an MPE increases, the robot's consumed power reduces.
We further conduct a comparative simulation study between the properties of previously used ellipsoids (MVE and MFE) and the newly defined ellipsoid (MPE). The study is conducted using numerical techniques based on assumptions, such as considering quasistatic properties for the hydraulic system. Compared to Mansouri et al. [16,17], our proposed power ellipsoid does not require complex evaluations and mainly focuses on the specification of the fluid-driven actuation system.

Materials and Methods
The following subsections explain the concept of manipulability of robotic manipulators and the manipulability ellipsoids and describe the manipulability measures. We then address the power measurements required to drive a fluid-actuated robotic arm. Finally, we define the novel power ellipsoid and describe the test rig and robot equations. Let [x e y e z e θ x θ y θ z ] T denote the generalized coordinates of the robot's end effector with respect to the reference frame, x R y R z R , which is composed of parameters describing Robotics 2022, 11, 32 3 of 15 the robot geometry and joint variables. The reference frame is shown in Figure 1 in the side view (x R z R plane) of a typical hydraulic manipulator with two degrees of freedom (DOFs): Let [ ] denote the generalized coordinates of the r with respect to the reference frame, , which is composed of par the robot geometry and joint variables. The reference frame is shown side view ( plane) of a typical hydraulic manipulator with two d (DOFs): , and , . In general modelling, the robot's end effector poss angular velocity components [ , , , , , ] with respect to th { }. Configuration and velocity/torque limitations are limiting facto bot's motion. Increasing the velocities or accelerations of the end effecto trajectory intensifies the effects of the manipulator kinematics and dy violate the workability of actuators. Locating the manipulator in an im tion can also increase these effects.
An improper configuration is defined when the robot's configur singular situation, e.g., the isometric index is close to zero. Enabling m proper trajectories by measuring and improving quantitative tools com bot's disability. Since the formulations used in this study are based equations, the following presents the manipulator's kinematics concep The manipulator DOF is denoted by n. Meanwhile, the end effecto space is shown by s. Assuming the robot is nonredundant, which mean matrix of the manipulator, ( ), is given by: In (1), ⃗ = [ … ] ∈ ℛ , ⃗ = [ … ] ∈ ℛ represents the v lar or linear displacement and velocity for revolute or prismatic join deed, ⃗ describes the relative displacement between two adjacent lin 1).
The end effector velocity vector, ⃗ ∈ ℛ , is generally defined as when s = 6. Note that if (det ( ) ≠ 0), i.e., when the Jacobian matrix has Configuration and velocity/torque limitations are limiting factors affecting the robot's motion. Increasing the velocities or accelerations of the end effector along the desired trajectory intensifies the effects of the manipulator kinematics and dynamics, which may violate the workability of actuators. Locating the manipulator in an improper configuration can also increase these effects.
An improper configuration is defined when the robot's configuration is close to the singular situation, e.g., the isometric index is close to zero. Enabling movement along improper trajectories by measuring and improving quantitative tools compensate for the robot's disability. Since the formulations used in this study are based on the kinematics equations, the following presents the manipulator's kinematics concept.
The manipulator DOF is denoted by n. Meanwhile, the end effector DOF in Cartesian space is shown by s. Assuming the robot is nonredundant, which means n ≤ s, the Jacobian matrix of the manipulator, J(q), is given by: In (1), q n ] T ∈ R n represents the vector of joint angular or linear displacement and velocity for revolute or prismatic joints, respectively. Indeed, → q describes the relative displacement between two adjacent links (see q j in Figure 1). θ z T when s = 6. Note that if (detJ(q) = 0), i.e., when the Jacobian matrix has full rank, the robot would be in a nonsingular condition. In this case, the system can be kinematically evaluated. Hereafter, the dependence on q will be omitted for notational compactness.

Manipulability Ellipsoids
In order to correlate the consumed energy of a fluid-power-driven manipulator to the manipulability concept, we define the existing manipulability ellipsoids: the manipulability velocity ellipsoid (MVE) and manipulability force ellipsoid (MFE). Having these ellipsoids, we can visualize the directions of velocity (MVE) or force (MFE) at the end effector [18]. The definitions of MVE and MFE will later be used to correlate the consumed power to the manipulability of the robot. This helps us observe the variation of the power required to run the robotic manipulator in terms of variations in the robot's kinematic parameters.
In order to analyze the properties of each manipulability ellipsoid, we use quantifiable tools such as determinants, eigenvalues and eigenvectors. This generalizes the singular value decomposition (SVD) of the mapping matrix (e.g., K and K L ) in the manipulability ellipsoid [19]. Using SVD, some manipulability measures are defined to quantitatively observe the variation of the consumed power.
The manipulability velocity ellipsoid is obtained using kinematics equations and helps us visualize the velocity of the end effector. In (1), if we map the unit circle . → X e , it will result in the manipulability velocity ellipsoid, which is expressed by: This is illustrated in Figure 2. As shown, applying the norm, || . q|| = 1 , results in a circle in the joint space. The circle is mapped through the Jacobian matrix into an ellipsoid in the Cartesian space, called the manipulability velocity ellipsoid. By mapping the normalized value of torque by J −T , we obtain the MFE.
Robotics 2021, 10, x FOR PEER REVIEW effector [18]. The definitions of MVE and MFE will later be used to correlate the power to the manipulability of the robot. This helps us observe the variation of required to run the robotic manipulator in terms of variations in the robot's parameters.
In order to analyze the properties of each manipulability ellipsoid, we use ble tools such as determinants, eigenvalues and eigenvectors. This generalizes lar value decomposition (SVD) of the mapping matrix (e.g., and ) in the bility ellipsoid [19]. Using SVD, some manipulability measures are defined t tively observe the variation of the consumed power.
The manipulability velocity ellipsoid is obtained using kinematics equ helps us visualize the velocity of the end effector. In (1), if we map the ⃗ ∈ ℛ |‖ ‖ = ∑ = 1 by to the space of end effector velocities, ⃗ , it w the manipulability velocity ellipsoid, which is expressed by: This is illustrated in Figure 2. As shown, applying the norm, ‖ ‖ = 1, r circle in the joint space. The circle is mapped through the Jacobian matrix into a in the Cartesian space, called the manipulability velocity ellipsoid. By mappin malized value of torque by , we obtain the MFE. As concluded from (2) and in Figure 2, in MVE, a Jacobian matrix is sh mapping matrix of velocity. Therefore, three submatrices of are derived as fo thogonal matrices and ∈ ℛ × is a diagonal matrix with non-negative rea on the diagonal and can be written as: As concluded from (2) and in Figure 2, in MVE, a Jacobian matrix is shown as the mapping matrix of velocity. Therefore, three submatrices of J are derived as follows: T ∈ R s×s are orthogonal matrices and Σ MVE ∈ R n×s is a diagonal matrix with non-negative real numbers on the diagonal and can be written as: Assuming that the manipulator is in a nonsingular condition, the following relationship can be expressed as: Moreover, the force ellipsoid can be defined by the eigenvalues of JJ T −1 , which are equal to 1/σ n (n = 1, 2, . . . , n). Thus, the principal axes of the force ellipsoid are → u 1 /σ 1 , . . . , → u n /σ n . While the MVE reflects the uniformity of the velocity of the end effector, the MFE reflects the force/torque applied from/to the end effector.

Manipulability Measures
The manipulability measure is used as a quantitative index to evaluate the performance of the robotic manipulator. For each ellipsoid, we can define measures based on the matrices using SVD. The measures we used in evaluating the dexterity of robotic manipulators are in Table 1. Table 1. Manipulability measures used in the simulation study.

Measure Equation Description
Yoshikawa's measure Yoshikawa [18] is a measure of manipulability which is in relation to the volume of the ellipsoid (Σ Y = c n ω) where ω = σ 1 . . . σ n and c n .
Condition number The condition number of a mapping matrix measuring the directional uniformity of the ellipsoid.

Isotropy index
The ratio of the length of minor semiaxis to the length of major semiaxis of the manipulability velocity ellipsoid.
The eccentricity of the ellipsoid and the ability of the end effector to move in a desired direction.
1 σ 1 and σ n are the largest and smallest eigenvalues of JJ T .

Power Consumption in Fluid-Power-Driven Robotic Arms
In the code developed, the hydraulic system is considered to have quasistatic properties. The actuator dynamics are slowly changing, and the capacitances of the cylinder Robotics 2022, 11, 32 6 of 15 chambers and hydraulic lines are neglected. Moreover, it is assumed that the mechanical power that enables the actuator to move back and forth is measured at the inlet of the actuator. Since the loss in the values of piping and connections cannot be correlated to the kinematics of the robot, they are ignored.
The 4th order Runge-Kutta method is used to generate the numerical simulation model described in this paper. System parameters used in simulations, formulation of the actuators and the method of modelling are described in detail in Banthia et al. [20].
For the hydraulic actuator shown in Figure 3, the mechanical power supplied by the hydraulic pump is the product of the input pressure and flow of the hydraulic fluid. The equation used to calculate the required power depends on the direction of the link motion when the piston either moves back ( .
where P O is the input pressure, supplied by the pump, and Q O is the flow of the hydraulic fluid. By substituting Q O = .
x s A O into (2), the following equation results: where A O denotes the annulus area of the piston. As shown in (7), the mechanical power required to run the hydraulic link is calculated by measuring the pressure in the inlet chamber, P O and calculating the linear velocity of the rod, .
x s . According to Figure 3, when .
x s,i < 0, the consumed power is obtained using P I .
x s A I .
Robotics 2021, 10, x FOR PEER REVIEW Eccentricity measure The eccentricity of the ellipsoid ability of the end effector to mo desired direction. 1 and are the largest and smallest eigenvalues of .
where is the input pressure, supplied by the pump, and is the flow of the hy fluid. By substituting = into (2), the following equation results: = where denotes the annulus area of the piston. As shown in (7), the mechanica required to run the hydraulic link is calculated by measuring the pressure in t chamber, and calculating the linear velocity of the rod, . According to Figure   , < 0, the consumed power is obtained using . Bernoulli's equation is used to express the power in terms of robot kinematic bles. To apply Bernoulli's equation, it is assumed that the flow is inviscid and stea locity pattern constant). The fluid is incompressible, and the friction losses are neg In this study, it is also assumed that the energy required for fluid elevation is ne with respect to the total energy of the system. By applying these assumptions, Ber equation for the system shown in Figure 3 is written as in (8). 1 where denotes the hydraulic fluid density. and are the linear velocity fluid at the outlet of the pump and the inlet chamber of the actuator. The pump pressure and the inlet chamber pressure are denoted by and . Depending on Bernoulli's equation is used to express the power in terms of robot kinematics variables. To apply Bernoulli's equation, it is assumed that the flow is inviscid and steady (velocity pattern constant). The fluid is incompressible, and the friction losses are negligible. In this study, it is also assumed that the energy required for fluid elevation is negligible with respect to the total energy of the system. By applying these assumptions, Bernoulli's equation for the system shown in Figure 3 is written as in (8). x s are the linear velocity of the fluid at the outlet of the pump and the inlet chamber of the actuator. The pump output pressure and the inlet chamber pressure are denoted by P p and P s . Depending on the direction of the link motion, P s can be either P O or P I . The pump pressure, P p , is set to be constant.
By applying the continuity equation to the system, given that there is no leakage in the system, the following equation can be written: Combining (8) and (9) gives the pressure at the inlet chamber: where Substituting (6) into (2) results in: The angular velocity of the joint, . q, and the linear velocity in the hydraulic actuator, .
x s , are related using the geometrical correlation mechanism, shown in Figure 1, as follows: For a robot with n DOFs, the total power needed to run all n hydraulic actuators is calculated by summing the power needed for each actuator, i.e., W = n ∑ i=1 W i . Equation (13), for the n-DOF robot, can be rewritten as: . → where . → q n ] T ∈ R n . The matrix of velocity, D ∈ R n×n , maps the joints angular velocities to the linear velocity of hydraulic actuators and is given by: Combining (12) and (14) results in the symbolic form of the power-velocity equation as below: where → W = [W 1 . . . W n ] T ∈ R n is the vector of power consumed by each actuator. The mapping matrix of energy, K, is defined as follows: where In accordance with the properties of the real hydraulic system, it is observed that P p . Therefore, the linearized form of the mapping matrix of energy is expressed as: The concept presented in (17) and (18) is then used to define the PE.

Power Ellipsoid (PE)
As presented in (17) and (18), the power supplied by the hydraulic pump is directly correlated to the angular velocity of the robot's joints. In this equation, if we map the unit Figure 4 depicts the concept of MPE. As seen, the mapping matrix of energy is K in the nonlinear form and K L in the linear form. As presented in (17) and (18)  Here, similar measures for MPE are used to correlate the power consumption and manipulability of the robot. According to the definitions used for the manipulability of robotic manipulators, it can be concluded that by decreasing the isotropy measure of MPE, , , the consumed power by the robot reduces. Table 2 summarizes the qualitative relationship between the manipulability measures of MPE and the power required to run the hydraulic robot. The signs shown in this table indicate the increasing (↑) and decreasing (↓) relationship between the measures and energy.

Simulation Setup
The 3-DOF hydraulic robot was modelled and simulated to see whether the ellipsoids had any correlation with the power consumed. Figure 5 depicts the manipulator Here, similar measures for MPE are used to correlate the power consumption and manipulability of the robot. According to the definitions used for the manipulability of robotic manipulators, it can be concluded that by decreasing the isotropy measure of MPE, Σ I,PE , the consumed power by the robot reduces. Table 2 summarizes the qualitative relationship between the manipulability measures of MPE and the power required to run the hydraulic robot. The signs shown in this table indicate the increasing (↑) and decreasing (↓) relationship between the measures and energy.

Simulation Setup
The 3-DOF hydraulic robot was modelled and simulated to see whether the ellipsoids had any correlation with the power consumed. Figure 5 depicts the manipulator with coordinate frames assigned to the links. Using the Denavit-Hartenberg notation, the forward kinematics is calculated by multiplying the transformation matrices. The generalized coordinate frames are attached to the robot, which starts from the base frame, {x o y o z o }, and ends at the end effector, {x e y e z e }. The transformation matrix, which correlates the position of the end effector to the base frame, is given by Maddahi et al. [21]: , 10, x FOR PEER REVIEW To evaluate the manipulability of the robot and to measures, the Jacobian matrix is used and is expressed as in The specification of the analyzed part of the robot is list are adopted from CAD drawings. The required parameters draulic system are shown in Table 4. All actuators have th exclusion of different strokes [22].  In (21), l 1 , l 2 and l 3 are the lengths of robot links, which are shown in Figure 5. c i and s i denote cos(θ i ) and sin(θ i ), and c ij and s ij represent cos θ i + θ j and sin θ i + θ j , respectively.
To evaluate the manipulability of the robot and to calculate the manipulability measures, the Jacobian matrix is used and is expressed as in (22).
The specification of the analyzed part of the robot is listed in Table 3. The dimensions are adopted from CAD drawings. The required parameters for each actuator and the hydraulic system are shown in Table 4. All actuators have the same parameters with the exclusion of different strokes [22].

Simulation Results
This simulation aimed to measure the power consumed by the hydraulic actuators to run the robot, to calculate the manipulability measures and validate the relationship between the required power and manipulability measures. A program was developed in C++ to emulate the characteristics of this robotic manipulator. The end effector of the robot was programmed to move along a circular path. The actual trajectory of the robot's end effector is illustrated in Figure 6. The angular displacement of each joint is shown in Figure 7.

Simulation results
This simulation aimed to measure the power consumed by to run the robot, to calculate the manipulability measures and v between the required power and manipulability measures. A pro C++ to emulate the characteristics of this robotic manipulator. The was programmed to move along a circular path. The actual traje effector is illustrated in Figure 6. The angular displacement of each 7.

Simulation results
This simulation aimed to measure the power consumed by to run the robot, to calculate the manipulability measures and v between the required power and manipulability measures. A pro C++ to emulate the characteristics of this robotic manipulator. The was programmed to move along a circular path. The actual traje effector is illustrated in Figure 6. The angular displacement of each 7.        By using forces and linear velocities, the power consumed in each actuator is calculated. The angular velocity of each joint is read by the encoders. By using (23), the linear velocity of each actuator is obtained. The variation of pressure is depicted in Figure 9. The total power is the summation of W 1 , W 2 , and W 3 . velocity of each actuator is obtained. The variation of pressure is depicted in Figure 9. The total power is the summation of , , and . Figure 10 depicts the values of manipulability measures for the developed force ellipsoid, MFE, the manipulability velocity ellipsoid, MVE, and the manipulability energy ellipsoid, MPE. The measures include Yoshikawa's measure ( , ), condition number ( ), isotropy index ( ) and eccentricity measure ( ). As illustrated, the relationship between the energy measures and velocity or force measures can, respectively, be found.    Figure 10 depicts the values of manipulability measures for the developed force ellipsoid, MFE, the manipulability velocity ellipsoid, MVE, and the manipulability energy ellipsoid, MPE. The measures include Yoshikawa's measure (Σ Y, ), condition number (Σ C ), isotropy index (Σ I ) and eccentricity measure (Σ E ). As illustrated, the relationship between the energy measures and velocity or force measures can, respectively, be found.
The relationship between the total mechanical power required to run the robot and the manipulability measures is illustrated in Figure 11. The experimental results confirm the analytical conclusions, which are presented in Table 1.
By increasing the isotropy index of PE, the required power of the pump increases as shown in the circled area A in Figure 11. Comparing Figures 9 and 10, we find the energy directly changes by altering the Yoshikawa's measure and isotropy index and has an adverse relationship with the condition number and eccentricity measure. By comparing |J| and |K L |, we conclude that in addition to less energy consumption, which results in minimizing the torque/force of actuators, changing the manipulability can avoid singularity configurations.  The relationship between the total mechanical power required to run the robot and the manipulability measures is illustrated in Figure 11. The experimental results confirm the analytical conclusions, which are presented in Table 1.
By increasing the isotropy index of PE, the required power of the pump increases as shown in the circled area A in Figure 11. Comparing Figures 9 and 10, we find the energy directly changes by altering the Yoshikawa's measure and isotropy index and has an adverse relationship with the condition number and eccentricity measure. By comparing | | and | |, we conclude that in addition to less energy consumption, which results in minimizing the torque/force of actuators, changing the manipulability can avoid singularity configurations. Robotics 2021, 10, x FOR PEER REVIEW Figure 11. Total mechanical power required to run the robot and corresponding determin mapping matrices of velocity and power.

Cross-correlation Analysis
Cross-correlation analysis was used to investigate the degree to which MPE s are correlated with the power consumption of the manipulator. Cross-correlation pares data of an MPE signal with a time-shifted (lagged) version of the power sign investigate whether a signal is stationary or non-stationary, values of four measures computed, as shown in Table 5.

Cross-Correlation Analysis
Cross-correlation analysis was used to investigate the degree to which MPE signals are correlated with the power consumption of the manipulator. Cross-correlation compares data of an MPE signal with a time-shifted (lagged) version of the power signal. To investigate whether a signal is stationary or non-stationary, values of four measures were computed, as shown in Table 5. Table 5. Measures used to investigate whether a signal (e.g., values of PE) is stationary or nonstationary.

Measure Equation Description
Mean value (first moment) Kurtosis is the measure of the peakedness of the probability distribution of a real-valued random variable [23]. 2 Skewness is the measure of the asymmetry of the probability distribution of a real-valued random variable [24].
The cross-correlation coefficient of the MPE signal (PE) and the power consumption signal (PC) for a given window of data is defined as in (24).
In (24), N is the size of the data set or window size of the time series PE and PC, and i refers to the i th element from two data sets. σ 1,PE and σ 1,PC are the mean values of PE and PC within the defined window. τ is an integer value for a time series analysis that determines how far apart the data points which are being compared are from each other (also known as the lag). As observed, the numerator of (24) is equivalent to the covariance of PE and PC with a delay τ, while the denominator is the product of the standard deviation of PE and PC. In this study, all signals were recorded with 2 ms sampling time or at a sampling frequency of 500 Hz.
Both PE and PC signals were examined to investigate whether they were stationary. Results show no upward or downward slope or jump in level in any of the four moments (σ 1 to σ 4 ). The observed trends were steep, exponential or approximately linear; therefore, the signals PE and PC were strongly stationary. It was found that the values of both pairs had a cross-correlation coefficient value of close to 1 (mean of 0.882), indicating that the two signals were similar.

Discussion
Manipulability measures and ellipsoids are useful tools to investigate the properties of robot end effectors such as velocity, acceleration and force. They are normally employed as indicators to quantitatively assess the robot's performance. One of the factors which can be considered to evaluate efficiency is the amount of power required to run the robot. In this study, the manipulability power ellipsoid, PE, was defined and correlated to the kinematic parameters of hydraulic manipulators.
To find the relationship between the MPE and consumed power, the coordination systems attached to the manipulator and the Jacobian matrix were first introduced. The manipulability of the robot was defined. We introduced the manipulability velocity, force and energy ellipsoids.
It was demonstrated that the increments of isotropy measure of MPE result in decreasing the consumed power. The concept was validated using a set of simulation studies performed on the main actuators of a 3-DOF robotic manipulator. We propose that the technique to correlate the power measurement to the robot kinematics is practical and provides useful insights for energy consumption. However, more mathematical modelling and validations are required before this correlation can be used as a tool to improve the energy efficiency in applications of hydraulic manipulators.
Quantification of straightforward energy consumption is an important step towards employing energy-efficient robotic systems in a sustainable manufacturing facility. An energy consumption model is normally case-dependent and complex to generate, while requiring continuous updates, as affecting parameters are changing, yet a key to analyze and improve the energy efficiency of robotic systems. Through a proof-of-concept study, we proposed that by implementing the proposed technique in this article, the energy consumption of fluid-driven robotic systems may simply be predicted by understanding and quantifying the characteristics of kinematics measures of the robotic arm, such as manipulability measures.
The proposed algorithm still needs further validations using robotic manipulators with different mechanisms. Therefore, future work will focus on validating the results of the simulation study presented in this paper, using the experimental platform that consists of a set of hydraulic actuators. We are currently conducting a comparative study to investigate the correlation among the MPE, its measures and other existing indices. We are working with machine learning techniques to design an optimal trajectory for a given task while minimizing power consumption.
In the future, this work can be used to find the optimal trajectory for a given task that is based on minimizing the power consumption. This can be done by combining the method and measures provided in this paper with an optimization method, such as genetic algorithms or artificial intelligence, and using the proposed manipulability power ellipsoid (MPE) as the cost function. Studies on the implementation of the same concept for other types of robotic systems, such as electric-driven robotcs, will be conducted in the future, as it is believed that the algorithms could be optimized and polished in order to be used for all types of robotoic manipulators.