1. Introduction
Robotic systems have become an important part of modern industries and are expected to revolutionize the way household chores are done in the coming days. The ability of robotic systems to do repetitive and tedious tasks day in and day out makes them an asset. However, human safety around robotic systems is an important concern. Traditionally, in industrial settings, the robotic system and humans are separated by ensuring appropriate segregation measures in the workspace. However, for service robots involved in household chores, such segregation is difficult to achieve. Robotic systems may unwittingly harm humans around them due to abrupt contact or carelessness of human operators. To mitigate such safety risks, introducing compliance to the joints has been widely theorized and researched. One way to do this is by active compliance, which mimics mechanical compliance by using sensors and actuators [
1]. Shetty and Ang [
2] used a force feedback mechanism to achieve compliance for a rigid-joint robot. Zinn et al. [
3] on the other hand relied on distributed macro–mini actuation. They used two different actuators—a large and a small actuator. High-torque, low-frequency movement was ensured by the large actuator, while the small actuator was responsible for low-torque, high-frequency movements.
Active compliance brings a lot of flexibility to robotic systems. However, many times, due to failure of sensors or low sample frequency, active sensors may become unreliable, and inadvertently the safety of humans may be in question. Passive compliance mechanisms, in general, are more reliable from a safety perspective. Variable stiffness joints are relatively inexpensive and efficient in ensuring no abrupt contact with humans occurs. Upon impact with humans, variable stiffness joints rapidly lower the stiffness of joints to avoid any major damage to the human operator. Tonietti et al. [
4] demonstrated a robotic arm comprising linear actuator-aided variable stiffness. Rotary-type permanent magnets were used by Yun et al. [
5] to achieve variable stiffness. The variable stiffness joints should have appropriate stiffness during the movement of the manipulator so that the robotic system can carry the payload. According to Yoo et al. [
6], in a typical human living environment, variable stiffness joints should be capable of generating more than 10 Nm torque for supporting a 1 kg payload at a 1 m length of robot manipulator. However, using traditional electric motors may make the variable stiffness joints bulky. To mitigate this, Yoo et al. [
6] proposed variable stiffness joints made neodymium–iron–boron ring-type permanent magnets. Magnets are sometimes used in robotic systems along with direct current motors [
7].
In several works in the literature, to ensure the proper design of the variable stiffness joints, optimization theory is used. Hyun et al. [
8] used a response surface methodology (RSM) to derive empirical relations of torque and weight with respect to various design parameters of permanent magnets like inner stator width, outer stator width, and magnet height. Yoo et al. [
6] did a similar analysis and reported a significant improvement over the baseline model. Using a finite element modeling route, Choi and Yoo [
9] designed a Halbach magnet array. They used numerical optimization techniques to optimize their design. Very recently, Song et al. [
10] discussed a two-step optimization methodology to optimize the maximum speed and impact force reduction capability of variable stiffness robots.
Apart from variable stiffness joints, optimization theory has been applied to other facets of robotic system design as well. Hsiao et al. [
11] designed and optimized high-speed robotic arms to achieve weight reduction, moment of inertia reduction, and deformation reduction. They used a finite element analysis to develop the necessary sampling dataset to generate the empirical RSM models. They reported improvements of about 16, 23, and 20% in the weight reduction, moment of inertia reduction, and deformation reduction, respectively. Chau et al. [
12] designed and optimize a compliant planar spring using a Kriging model. They deployed a multi-objective genetic algorithm to achieve the optimization task. Zhao et al. [
13] used a non-dominated sorting genetic algorithm to design a wall-climbing robot for stability and weight reduction. They were able to reduce the vibration by 6.5% and the weight of the magnetic abortion unit by 9.7%.
From the above literature review, it is clear that recently researchers have started tapping into the potential of optimization methodologies to design and optimize robotic systems. However, the use of nature-inspired optimization techniques is very limited in this field [
14,
15]. Mostly, nature-inspired optimization techniques are used in path planning applications in robotics [
14,
15]. In this paper, a very recently developed nature-inspired metaheuristic algorithm called the Non-dominated Sorting Whale Optimization Algorithm (NSWOA) is applied. NSWOA is a previously developed algorithm by Jangir and Jangir [
16] in 2017 for continuous optimization problems. However, so far it has not been applied to robot design problems. Since NSWOA generates the optimized solutions in form of Pareto fronts, a multi-criteria decision-making method called MARCOS (Measurement of Alternatives and Ranking according to COmpromise Solution) is applied to extract solutions for predefined application scenarios. The proposed hybrid NSWOA-MARCOS method is then applied to a variable stiffness joint design and optimization problem. The rest of the paper is arranged as follows. The next section details the methodology used. The NSWOA and MARCOS methods are briefly discussed in this section. The description of the case study considered in this paper is discussed in
Section 3.
Section 4 begins with the understanding of the parametric effect of the design features on the responses. Next, the parametric optimization is carried out by NSWOA-MARCOS. Finally, the knowledge derived based on this study is summarized in the
Section 5.
2. Methodology
In this paper, a hybrid NSWOA-MARCOS method for the optimal design of robotic components is proposed. The method comprises two main parts: (1) a nature-inspired metaheuristic algorithm infused with a non-dominated solution sorting strategy to generate a Pareto front and (2) a multi-criteria decision-making method to extract the most viable solutions from the Pareto front depending on the application scenario.
2.1. Whale Optimization Algorithm
The whale optimization algorithm (WOA) was initially propounded by Mirjalili and Lewis [
17] in 2016 and is a nature-inspired metaheuristic algorithm that mimics the social behavior of humpback whales.
Humpback whales identify the prey’s location and encircle it [
18]. Similarly, the whale optimization algorithm assumes that the current best solution is the prey. The position of the other search agents or candidate solutions are updated with reference to the best solution. The following equation represents the mathematical equivalent of prey encirclement by humpback whales
where
and
are the location vector of the best whale and any other whale, respectively.
indicates the current iteration. At the
iteration, the location of the whale is given as
. The “.” between vectors represent an element-wise product. The “||” represents the absolute values. The coefficient vectors
and
can be calculated as
With the increase in iterations, linearly decreases from 2 to 0. is a random vector having a [0, 1] range. The various locations around the best solution so far are computed by adjusting and vectors.
The first phase of WOA is the exploitation phase based on the bubble-net feeding method of whales. Humpback whales make use of a shrinking circle and spiral-shaped path to swim around the prey [
18]. Two approaches are considered to model this feeding technique mathematically: shrinking encircling mechanism and spiral position updating [
19].
By decreasing the value of
in Equation (3), the shrinking encircling mechanism of humpback whales is realized [
17].
To update the spiral position [
17], the first step is to compute the distance between the whales and prey. The helix-shaped movement of the humpback whale is mathematically represented as
where
is a constant to define the logarithmic spiral’s shape,
is a random number within the range of [−1, 1].
At any given instant, the humpback whale may update its position either by the shrinking circle method or by the spiral path method. To model these two methods simultaneously, a 50% probability is assigned to choose between either method, which is mathematically represented as [
17].
where
is a random number within the range of [0, 1], the value of
is responsible for switching between a spiral or circular movement.
For the exploration phase, i.e., to search the prey, the same approach of variation of the
vector is used [
18]. The humpback whales randomly search for prey. If
, a random whale is chosen. However, if
, the best solution is selected when changing the location of the search agents. Mathematically this is achieved by [
17]
where
is a random location vector selected from the present population.
2.2. Non-Dominated Sorting Whale Optimization Algorithm
To model the multi-objective version of the WOA an archive is first incorporated into this algorithm. This archive is combined to store the best non-dominated solutions acquired so far. The solution procedure of the multi-objective whale optimization algorithm [
20,
21] is very analogous to that of the WOA. That is, solutions are established considering the bubble-net feeding method of the humpback whales. A leader selection mechanism is introduced to elect the solutions from the archive [
22]. A roulette wheel is also incorporated in the MOWOA to select solutions from the less occupied areas of the archive [
23]. The improvement of the distribution of solutions in the archive across all objectives is made using the following probability function [
24]
where
is a constant and should be greater than 1 and
is the number of solutions in the neighborhood of ith solution.
Meanwhile, the archive has a limit to storing non-dominating solutions, and it might become full with the progress of iterations. So, a mechanism is also implemented to eliminate the undesirable solutions from the archive. Undesirable solutions are those solutions that have many adjacent solutions. To discard the undesirable solutions from the archive of NSWOA a probability is employed to give a high value to the undesired solutions. The probability function is the inverse of the previous probability function (Equation (9)), which is used to select the optimal solution from the archive, and it is given as follows [
24]
2.3. MARCOS
The relation between any alternative and the ideal and anti-ideal solutions forms the backbone of the MARCOS method. Certain utility functions are determined based on these relationships and the best possible alternative is then determined by ranking all the alternatives.
The Pareto front obtained from the NSWOA serves as the input matrix to the MARCOS method. Any multi-criteria decision-making (MCDM) method has an input matrix or the initial decision matrix of the form
.
and
represent the number of alternatives and the criteria. Here, the solutions in the Pareto front are the alternatives, and the two responses, i.e., torque and weight, are the criteria. The initial decision matrix
for any MCDM can be written as,
here,
x11,
, etc. are the candidate solutions.
For MARCOS, an extended decision matrix
of the following form is used,
here,
,
, etc. are the anti-ideal (AAI) solutions and
,
, etc. are the ideal (AI) solutions.
The anti-ideal (AAI) solutions and the ideal (AI) solutions in the extended decision matrix are determined using the following relations,
here,
and
are the benefit and the cost criteria respectively. In terms of this paper, the torque is the benefit criteria and the weight is the cost criteria.
The extended decision matrix
is then normalized using the following relations to form the normalized decision matrix
.
The extended decision matrix
is of the following form,
The weighted decision matrix
is then determined by multiplying the weight vector with the normalized decision matrix
. The weight vector is expressed as,
here,
,
, etc. are the weights assigned to each of the criteria. Additionally, the sum of the weights must be equal to unity, i.e.,
The weighted decision matrix
can be written as,
Next, the utility degree of the alternatives
is determined using the following two equations,
here,
is the sum of the elements of the weighted decision matrix
.
Si is determined using the following equation,
The utility function of the alternatives
is then computed,
here,
and
are the utility functions in relation to anti-ideal (AAI) solutions and the ideal (AI) solutions.
Finally, the alternatives are ranked from rank 1 to in descending order of values.
3. Problem Description
The design and optimization of robotic components is an essential task due to the growing demand for robotic systems in modern world applications. However, given the intricate nature of associated systems, this task is a challenging one. In this work, a nature-inspired evolutionary algorithm, namely, the whale optimization algorithm, is employed for the design optimization of a variable stiffness joint in a robot manipulator. The case study considered here is adapted from Yoo et al. [
6].
Figure 1a shows the typical location of variable stiffness joints in the robot manipulator arm and
Figure 1b shows the architecture of the three-ring permanent magnet-type variable stiffness joint design considered by Yoo et al. [
6]. The variable stiffness joint in a robot manipulator is aimed to be used in a robot arm that can support a 1 kg payload. The objective of the optimization task is to minimize the overall weight of the variable stiffness joint. However, it is assumed that the rotational stiffness of the variable stiffness joint must be at least 10 Nm. In single-objective optimization, this can be achieved by modeling the torque requirement as a constraint to the problem and formulating the problem as a constrained optimization problem. However, in this paper, the optimization problem is formulated as a Pareto optimization problem where maximizing rotational stiffness is the second objective. Further, in the archiving step of the non-dominated solutions, along with Equation (10), an additional constraint of rejecting solutions of rotational stiffness values < 10 Nm is applied. This ensures the entry of only those solutions to the non-dominated solution archive that met the rotational stiffness constraint criteria.
Further, the variable stiffness joint system is assumed to be composed of the stator and the rotor connected with the joint. The inner stator width (
), outer stator width (
), and magnet height (
) are considered as the design parameters. For the optimization phase, the following constraints are applied to the design variables
Based on a set of experiments conducted as per central composite design, Yoo et al. [
6] proposed the following second-order empirical relationships for the stated design variables and maximum torque and weight. It should be noted that
,
, and
are in mm, whereas maximum torque and weight are in Nm and kg, respectively. However, Equations (26) and (27) are provided in coded form. Thus, the range of all the variables in Equations (26) and (27) is ±1.
The empirical relationships in Equations (26) and (27) are used as the objective functions for the NSWOA algorithm. Yoo et al. [
6] have reported excellent estimation power of the equations with
of 0.97 and 1 for Equations (26) and (27), respectively.
, i.e., the coefficient of determination, indicates the amount of variance in the data captured by the models.
5. Conclusions
In this article, a novel approach to design the optimization of robotic components is presented. Since robotic systems must consider multiple objectives, a Pareto optimal front generation capable, non-dominated sorting whale optimization algorithm (NSWOA) is used. The NSWOA generates a multitude of non-dominated solutions which are then passed through a MARCOS algorithm to select the optimal compromise solution for a pre-specified application. A case study of optimal variable stiffness joint by using permanent magnets is tackled in this paper. The inner stator width, outer stator width, and magnet height are the design parameters considered to minimize the weight and maximize the torque generated. It was found during the design process that weight and torque generated are conflicting objectives in nature i.e., in general, one cannot be improved without deteriorating the other. With respect to the previously reported best solution in the literature, an improvement of about 1.8% in weight reduction was observed in the new design. The improvement is about 41% with respect to the original design. The study is limited by consideration of only one nature-inspired algorithm. In the future, this study can be extended to test the hybridization of other metaheuristics like gray wolf optimizer, multiverse optimizer with MCDMs like MABAC, CoCoSo, etc. The concept of fuzzification can also be introduced to account for uncertainty in design.