Structural Optimisation and Design of a Cable-Driven Hyper-Redundant Manipulator for Confined Semi-Structured Environments

Structural optimisation of robotic manipulators is critical for any manipulator used in confined semi-structured environments, such as in agriculture. Many robotic manipulators utilised in semi-structured environments retain the same characteristics and dimensions as those used in fully-structured industrial environments, which have been proven to experience low dexterity and singularity issues in challenging environments due to their structural limitations. When implemented in environments other than fully-structured industrial environments, conventional manipulators are liable to singularity, joint limits and workspace obstacles. This makes them inapplicable in confined semi-structured environments, as they lack the flexibility to operate dexterously in such challenging environments. In this paper, structural optimisation of a hyper-redundant cable-driven manipulator is proposed to improve its performance in semi-structured and challenging confined spaces, such as in agricultural settings. The optimisation of the manipulator design is performed in terms of its manipulability and kinematics. The lengths of the links and the joint angles are optimised to minimise any error between the actual and desired position/orientation of the end-effector in a confined semi-structured task space, as well as to provide optimal flexibility for the manipulators to generate different joint configurations for obstacle avoidance in confined environments. The results of the optimisation suggest that the use of a redundant manipulator with rigid short links can result in performance with higher dexterity in confined, semi-structured environments, such as agricultural greenhouses.


Introduction
The introduction of robotics in the middle of the twentieth century promoted the use of robots in many industrial and commercial automated contexts [1][2][3][4][5]. However, they are mostly used in fully-structured industrial environments [6,7], where the robot workspaces are clearly defined and known a priori. Only recently have robots been further utilised in semi-structured environments [8][9][10][11], such as in agriculture, in which the surroundings are dynamic and can be uncertain. Some work has presented strategies for introducing robotic manipulators into such environments [8,10,12,13]. Nevertheless, all the manipulators used are similar to conventional manipulators, with the 5-6 degrees of freedom (DOFs) that are found in a fully-structured industrial environment. Such robots have been proven to be inefficient, as the manipulators suffer from a lack of dexterity appropriate for the given workspace [10,14,15].
A semi-structured environment is defined as an environment with known landmarks that has no static topology, where the surroundings cannot be classified into, for example, rooms and walls [16]. It is instead a dynamic environment and subject to change during the operation of the manipulators [17]. This means any interactions with the environment Before implementing the robotic manipulators in semi-structured environments, their dexterity and manipulability must be considered in the optimisation process to ensure appropriacy and applicability. In the earlier work on semi-structured environments, the manipulability criteria, despite their critical role, have not been considered. The research to date only considered conventional industrial manipulators, which are not applicable in confined semi-structured environments due to the long distance between joints. The distance between every two consecutive joints should be optimised for better manipulability in confined semi-structured environments. Therefore, this article aims at extending the aforementioned works by optimising the structure of hyper-redundant manipulators to enable them to work flexibly with high manipulability within their given workspace in confined semi-structured environments. The introduced optimisation methods can result in an optimal design of robotic manipulators as the manipulability measures are considered in the optimisation process. This makes the optimised manipulator highly applicable in such confined environments as the manipulability measures are critical in terms of assessing the performance of robotic manipulators in this type of environment.
In this study, the robotic arm's joint angles and link lengths are optimised using three manipulability measures and inverse kinematics optimisation to ensure the manipulator can move flexibly within its workspace with different configurations. Two algorithms are considered for the optimisation: Artificial Bees Colony (ABC) [38] and Fmincon [39], which is a built-in optimisation algorithm in MATLAB. Then, the optimisation results are validated using a hyper-redundant manipulator, which was created based on the optimisation results. The performance of the developed manipulator is tested in terms of manipulability by generating a collision-free trajectory in a confined semi-structured environment and analysing its manipulability at certain points within the environment. The results show that the optimised manipulator can achieve high manipulability even though the environment used for the validation is different from those used for optimisation. Based on the optimisation results, this implies that the manipulator can be used in any given confined semi-structured environment as it has high manipulability measures.
The remaining sections of this article are organised as follows: Section 2 presents the structural optimisation of a hyper-redundant manipulator to find the optimal length of links and joint angles that enable the manipulator to move with high manipulability. Section 3 shows the results of the proposed optimisation methods, along with the validation of the results. The conclusion is drawn in Section 4.

Structural Optimisation of Hyper-Redundant Manipulators
Optimisation is the process of determining the optimal solution among diverse alternatives; it is necessary for robotic manipulator design so that the optimal structure of joint parameters can be obtained. Several optimisation algorithms are available, including Genetic Algorithms (GA), Particle Swarm Optimisation (PSO) and Artificial Bees Colony (ABC). No optimisation algorithm can be classified as the unqualified best, as the performance of each specific algorithm depends on the problem structure, which means some objective functions can perform faster than others, depending on the problem. However, ABC generally showed successful and stable results when compared with GA and PSO when it was tested to minimise the makespan for job scheduling, and protein production problems [40,41], as it showed advantages in stability and not being trapped in local minima. Moreover, it does not employ external parameters such as crossover and mutation rates, as in GA [42]. Another widely used optimisation solver is the Fmincon solver; it is the main optimisation solver in the MATLAB Optimisation Toolbox [39], which performs better for nonlinearly constrained problems. Fmincon is based on interior point techniques for nonlinear minimisation. The flowcharts of both optimisation algorithms are shown in Figures 1 and 2.       Any optimisation problem can be defined as finding the vector of parameters that minimises or maximises an objective function: min or max f (z) = (z 1 , z 2 , . . . , z n ) (1) and constrained by: where z represents the vector consisting of parameters to be optimized. In this study, z i is defined as x i = [q, a] T , where q = (q 1 , q 2 , . . . , q n ) is the joint angles in degrees, and a is the link length of a hyper-redundant manipulator in mm. In this case, all the links are assumed to have an equal length. The objective function used in this study to optimise the joint parameters for n DOF is given by: where f (.) is the manipulability measure, which is used as the objective function, and q, q and a, a are the lower bound and upper bound for joint angles and the length of the link, respectively. When choosing lower and upper bound values, the joint structure and limitations were considered, as well as the designs in previous studies that showed a lack of flexibility to work in their given workspace [14,24,43]. Thus, the lower bound was chosen to avoid being too small to provide sufficient space for other components, such as pulleys, wires and joint sensors, and the upper bound was sufficiently large that its size would not affect its manipulability. That is, if the optimized length of the link were too small, pulleys, wires and sensors cannot be added to the joint due to space limitations; while, if it is too long, it will result in a manipulator that is not flexible enough to handle its workspace as with those in the aforementioned studies.
For the proposed redundant manipulator, in order to achieve high manipulability and dexterity, the optimised objective functions are based on manipulability measures and inverse kinematics. Based on the definitions of dexterity and manipulability, the objective functions are dependent on the joint parameters. The minimum and maximum values used in both optimisation algorithms are ±90 degrees for joint angles and 80 mm and 250 mm for the length of the link, respectively. The minimum value is chosen, as mentioned above, to provide sufficient space for the joint components, whilst the dexterity performance is used to define the maximum value as the link length, for which greater than 250 mm was proven to lack dexterity in any workspace [14].
The manipulator includes universal joints, as described further in Section 3. Thus, every joint is driven by two motors-one for each axis. The joint is attached to two pulleys that are connected to the motors by cables to enable each motor to drive one axis. The last joint is considered a general end-effector so that it can be incorporated into any type of bespoke gripper. Another aspect of the design is the base, which houses the motors, electronics and power equipment. The cables that are used in the design are of the Bowden cable type and are assumed in this work to have no effect on the kinematics and dynamics.

Optimisation Algorithms
The first optimisation algorithm used in our work to determine the optimal joint parameters is ABC as it can avoid being trapped in local minima wherein different joint parameters are tested and compared. ABC is a metaheuristic method inspired by foraging honeybees. There are three groups of bees in ABC: employed bees, onlooker bees and scout bees. Onlooker bees and scout bees are also called unemployed bees. Initially, the scout bees look for food until they find its source. Then, the employed bees and onlooker bees exploit the abundance (fitness) of the food source. The employed bees are linked to a certain food source and share information about the food by means of dancing. Their dancing is observed by onlooker bees to determine the best source using probability. The location of the food source reflects a solution in the ABC algorithm, and the abundance of food represents the fitness of the solution.
For the proposed structural optimisation, the initialisation stage of the algorithm starts by assigning initial values to the joint angles and length of the link. Then, it initialises the population and the limits of the variables. Next, it starts by assigning the initial solution, which is, in this study, zero for manipulability optimisation, with infinity as the distance between the current and desired end-effector coordinates. The corresponding position for every employed bee represents the values of the joint parameters. Abundance is shown by the manipulability index and the distance between the current and desired position of the end-effector.
After the initialisation process, the employed bees start searching for a new food source v m , which is more abundant than x m , which is already in their memory, and evaluate its i-th fitness using the following equation [38]: where i is a randomly selected parameter index, x k is a randomly chosen food source, and ϕ mi is chosen randomly between −1 and 1. Then, after the new food source is identified using Equation (3), the fitness of the new food source is calculated, and a greedy selection is performed between v mi and x mi , as follows: where f it m (x m ) is the fitness of the new food source, and f m (x m ) is the objective function value for the solution x m . Then, the information is shared by the employed bees with the onlooker bees, which determine the best solution among the fitness values using probability, p m . The best solution is determined using the roulette wheel fitness selection method [44], as in Equation (5): where SP is the size of the population, which was chosen, based on several trials, to be 500, divided equally into employed and onlooker bees so the employed bees share the information about their food source to recruit the corresponding onlooker bees. Based on the shared information, the onlooker bees choose the food source. Then, the employed bees, whose solution cannot be improved after certain trials, abandon their solution using abandonment criteria. As a result, they become scout bees, looking for a new source of food. The abandonment criteria define the control limit for jumping from local minima to global minima; that is, if the solution of the employed bees does not change after a number of iterations, they are converted to scout bees. The overall flowchart of the algorithm is shown in Figure 1.
The number of iterations in ABC is 300 for manipulability measures and 1000 for inverse kinematics. The objective function values for manipulability measures did not increase from around 250 iterations, which means the best value is reached at around 250 iterations. This was observed after several trials when testing the algorithm. Therefore, the number of iterations was reduced to 300 to save computational time. In Fmincon, on the other hand, the initial number of iterations was set to 5000, and the exit criteria depend on the tolerance, which is 1e-16. This means that although the number of iterations was set too high, the optimisation can stop after a few iterations as long as the objective function value is not increasing, as the results show in the next section. If, for example, the algorithm starts by assigning the number of joint angles, based on the given length of the manipulator workspace, the length of the links is identified. The next step is identifying the distance between every two consecutive joints at each iteration. Then, the algorithm starts by assigning the minimum number of joints and checking the manipulability and dexterity with iteratively assigned joint parameter values and storing the value of the corresponding objective function. Next, it assigns a different number of joint parameters and evaluates their objective function and compares the value of the current objective function with the previous objective function value. As the manipulability depends on the value of the specific joint angles, the subroutine has to iteratively assign different angles and check the manipulability measures.
Finally, after considering the optimisation results, a computer-aided design (CAD) model of the manipulator will be created. Every link should have two axes (pitch and yaw) so the manipulator can move freely and flexibly and cover the area within the range of the manipulator. As highlighted previously, the main reason for the low level of dexterity is the long length of the links. In other words, every consecutive similar motion (for example, yaw-yaw) should be close to one another, based on the optimised values of each joint's parameters, so the manipulator will be sufficiently dexterous for the specific task.

Kinematics of Hyper Redundant Manipulators
The optimisation in this study is based on inverse kinematics and manipulability. The different manipulability measures used for optimisation are based on the Jacobian matrix. The kinematics of the proposed hyper-redundant manipulator are solved using Denavit-Hartenberg (DH) parameters, as shown in Table 1. The model consists of n revolute joints along the z-axis of frames { 1} to { n} . As the manipulator is a universal joint-type manipulator, n is an even number. The forward kinematics function, h(.), is: where y = x, y, z, θ x , θ y , θ z T is the pose vector (the translational pose in the x, y and z axes, and the rotational angles with respect to the x, y and z axes) in the task space, and q = (q 1 , q 2 , . . . , q n−1 , q n ) T is the joint angle vector in the configuration space. Joint q n d n a n (mm) α n (deg) The transformation between the base frame and any other frame is given by the transformation matrix as: where d n is the offset along the Z n−1 axis; a n is the length of the common normal, which in our case is the length of the link connecting the previous yaw movement axis to the next pitch movement axis; and α n is the angle about the common normal, from the Z n−1 axis to the Z n axis. The inverse kinematics are used to find the function h −1 , which maps the task space to the configuration space: The configuration space of this type of manipulator is larger than the task space, which results in infinite inverse kinematics solutions with v number of equations and u unknowns, where v > u, with u = 6, which corresponds to 3 translation and 3 angular motions in a 3D space. Therefore, iterative methods using the Jacobian matrix are used, as analytical solutions do not exist. Deriving Equation (6) with respect to time gives: where J is the Jacobian matrix that maps the translational and the angular velocities of the joints to the pose velocities.
As the inverse of the Jacobian matrix can only exist for non-singular square matrices, the Moore-Penrose pseudoinverse is taken to give: y (11) where J + is the pseudoinverse of the Jacobian matrix, which is equal to: where J T is the transpose of the Jacobian matrix.

Manipulability Optimisation
Manipulability optimisation has been undertaken to ensure a high manipulability measure. Three different manipulability measures were used as objective functions to optimise the best joint parameter values. The first measure used is the Yoshikawa manipulability measure [32], which is given by: where det J J T is the determinant of the multiplication of the Jacobian matrix and its transpose. The Yoshikawa measure is equal to the product of the singular values of the Jacobian matrix. As this index is proportional to the velocity ellipsoid, the high singular values mean a high end-effector velocity for the unit norm: a manipulator is considered to have a high manipulability if its velocity ellipsoid is isotropic, i.e., has a sphere-shape, which means it can change its end-effector position in all three axes equally. The objective function values and the optimised joint parameters for the optimisation based on the Yoshikawa manipulability measure using Fmincon and ABC are shown in Figures 3 and 4. However, the Yoshikawa measure is not an accurate manipulability measure, as the velocity ellipsoid cannot depict the exact status of every singular value. For example, if one singular value is very small, the velocity ellipsoid can still be large as the other singular values are high. Furthermore, this measure is unbounded: therefore, another measure was used to compensate for the downside of the Yoshikawa manipulability measure: the Local Conditioning Index (LCI). The LCI is the ratio of the minimum and maximum singular values [45,46], Equation (14):   However, the Yoshikawa measure is not an accurate manipulability measure, as the velocity ellipsoid cannot depict the exact status of every singular value. For example, if one singular value is very small, the velocity ellipsoid can still be large as the other singular values are high. Furthermore, this measure is unbounded: therefore, another measure was used to compensate for the downside of the Yoshikawa manipulability measure: the Local Conditioning Index (LCI). The LCI is the ratio of the minimum and maximum singular values [45,46], Equation (14): where and are the minimum and maximum singular values, respectively. The result obtained from the LCI manipulability measure as an objective function for parameter optimisation is shown in Figure 5 and Figure 6:   However, the Yoshikawa measure is not an accurate manipulability measure, as the velocity ellipsoid cannot depict the exact status of every singular value. For example, if one singular value is very small, the velocity ellipsoid can still be large as the other singular values are high. Furthermore, this measure is unbounded: therefore, another measure was used to compensate for the downside of the Yoshikawa manipulability measure: the Local Conditioning Index (LCI). The LCI is the ratio of the minimum and maximum singular values [45,46], Equation (14): where and are the minimum and maximum singular values, respectively. The result obtained from the LCI manipulability measure as an objective function for parameter optimisation is shown in Figure 5 and Figure 6:   However, the Yoshikawa measure is not an accurate manipulability measure, as the velocity ellipsoid cannot depict the exact status of every singular value. For example, if one singular value is very small, the velocity ellipsoid can still be large as the other singular values are high. Furthermore, this measure is unbounded: therefore, another measure was used to compensate for the downside of the Yoshikawa manipulability measure: the Local Conditioning Index (LCI). The LCI is the ratio of the minimum and maximum singular values [45,46], Equation (14): where and are the minimum and maximum singular values, respectively. The result obtained from the LCI manipulability measure as an objective function for parameter optimisation is shown in Figure 5 and Figure 6:  The third manipulability measure that was used to optimise the parameters is a new measure introduced by Mariass et al. [47]. This measure is bounded and can be written as a function of the joint coordinates. The measure is given by: where and are the harmonic and arithmetic mean of eigenvalues, respectively, and is the trace of the matrix. Figures 7 and 8 show the result for the measure when used as an objective function for parameter optimisation.  The third manipulability measure that was used to optimise the parameters is a new measure introduced by Mariass et al. [47]. This measure is bounded and can be written as a function of the joint coordinates. The measure is given by: where H and A are the harmonic and arithmetic mean of eigenvalues, respectively, and tr is the trace of the matrix. Figures 7 and 8 show the result for the ∆H measure when used as an objective function for parameter optimisation.  The third manipulability measure that was used to optimise the parameters is a new measure introduced by Mariass et al. [47]. This measure is bounded and can be written as a function of the joint coordinates. The measure is given by: where and are the harmonic and arithmetic mean of eigenvalues, respectively, and is the trace of the matrix. Figures 7 and 8 show the result for the measure when used as an objective function for parameter optimisation.  It can be seen from the results of the three manipulability measures that the optimisation algorithms continued to reduce the links' length until it reached the lower bound limit of length between every two joints, resulting in a 16-DOF manipulator. This means the shorter the length of links, the better the manipulability that can be achieved. However, in practice, the length of the link should be sufficient to avoid design limitations by providing space for joint sensors, cables and pulleys. Hence, the lower bound of the link lengths is restricted to 80 mm. It can be seen from the results of the three manipulability measures that the optimisation algorithms continued to reduce the links' length until it reached the lower bound limit of length between every two joints, resulting in a 16-DOF manipulator. This means the shorter the length of links, the better the manipulability that can be achieved. However, in practice, the length of the link should be sufficient to avoid design limitations by providing space for joint sensors, cables and pulleys. Hence, the lower bound of the link lengths is restricted to 80 mm.

Inverse Kinematics (IK) Optimisation
Inverse kinematics optimisation was conducted to verify the performance of the manipulator, considering certain task space coordinates. It is based on the results from the manipulability optimisation in Section 2.2 in terms of the optimal number of joints. The objective function and constraints for joint parameters in this part are similar to equation (2), with (. , . , . ), Equation (18), as the objective function, which minimises the error between the current and desired positions of the end-effector. Additional workspace constraints were added to ensure the manipulator can perform flexibly in confined and enclosed workspaces. The workspace was constrained to resemble a box with dimensions of 1000 mm, 700 mm and 900 mm in the x, y and z coordinates, respectively. This means the end-effector should not go beyond this workspace. Another constraint was added to limit the middle joint to spaces no larger than 500 mm, 350 mm and 450 mm in the x, y and z coordinates, respectively.
where and are the current and desired end-effector positions in the x-axis, respectively. Similarly, and , and and are the current and desired end-effector coordinates in the y-axis and z-axis, respectively. The desired coordinates were 500 mm, 300 mm and 200 mm in the x, y and z coordinates, respectively. The result obtained by Fmicon for inverse kinematics optimisation was more accurate than the result obtained by ABC, as the error was nearly 0.002 mm using the Fmincon solver, while it was about 0.4 mm in ABC, as shown in Figures 9 and 10:

Inverse Kinematics (IK) Optimisation
Inverse kinematics optimisation was conducted to verify the performance of the manipulator, considering certain task space coordinates. It is based on the results from the manipulability optimisation in Section 2.2 in terms of the optimal number of joints. The objective function and constraints for joint parameters in this part are similar to Equation (2), with G(., ., .), Equation (18), as the objective function, which minimises the error between the current and desired positions of the end-effector. Additional workspace constraints were added to ensure the manipulator can perform flexibly in confined and enclosed workspaces. The workspace was constrained to resemble a box with dimensions of 1000 mm, 700 mm and 900 mm in the x, y and z coordinates, respectively. This means the end-effector should not go beyond this workspace. Another constraint was added to limit the middle joint to spaces no larger than 500 mm, 350 mm and 450 mm in the x, y and z coordinates, respectively.
where x and x d are the current and desired end-effector positions in the x-axis, respectively. Similarly, y and y d , and z and z d are the current and desired end-effector coordinates in the y-axis and z-axis, respectively. The desired coordinates were 500 mm, 300 mm and 200 mm in the x, y and z coordinates, respectively. The result obtained by Fmicon for inverse kinematics optimisation was more accurate than the result obtained by ABC, as the error was nearly 0.002 mm using the Fmincon solver, while it was about 0.4 mm in ABC, as shown in Figures 9 and 10:    All the optimised values of parameters obtained by the different objective functions used in this study are summarised in Table 2.  All the optimised values of parameters obtained by the different objective functions used in this study are summarised in Table 2. Based on the results of the different objective functions, either from manipulability measures or inverse kinematics, it is observable that the best performance of the manipulator can be obtained within a range between 90 mm and 150 mm for the length of the links. The joint angles also reached the lower and upper bounds of ±90 degrees. Therefore, considering the obtained results, the manipulator designed in this work has a length of 100 mm between every two consecutive similar joints, which satisfies the high dexterity needs of the manipulator.

Results Validation
Based on the optimisation results, a model for the optimised manipulator was developed using Autodesk Inventor, as shown in Figure 11. It consists of the housing for the actuators and electronic equipment, along with the manipulator. To account for dexterity, the universal joint is considered in the design, as shown in Figure 11 (right). The joint consists of two identical yokes joined by an octagon in the centre. Every yoke is connected to a pulley so it can perform two-axis (pitch and yaw) movements. The distance between every two consecutive joints is 100 mm, and the distance from the last joint to the tip is also 100 mm, as shown in Figure 12; the inner diameter of the joint is 60 mm, while its outer diameter is 90 mm. The joints' movements are animated, and they can move ±90 degrees in both the pitch and the yaw axes. Based on the results of the different objective functions, either from manipulability measures or inverse kinematics, it is observable that the best performance of the manipulator can be obtained within a range between 90 mm and 150 mm for the length of the links. The joint angles also reached the lower and upper bounds of ±90 degrees. Therefore, considering the obtained results, the manipulator designed in this work has a length of 100 mm between every two consecutive similar joints, which satisfies the high dexterity needs of the manipulator.

Results Validation
Based on the optimisation results, a model for the optimised manipulator was developed using Autodesk Inventor, as shown in Figure 11. It consists of the housing for the actuators and electronic equipment, along with the manipulator. To account for dexterity, the universal joint is considered in the design, as shown in Figure 11 (right). The joint consists of two identical yokes joined by an octagon in the centre. Every yoke is connected to a pulley so it can perform two-axis (pitch and yaw) movements. The distance between every two consecutive joints is 100 mm, and the distance from the last joint to the tip is also 100 mm, as shown in Figure 12; the inner diameter of the joint is 60 mm, while its outer diameter is 90 mm. The joints' movements are animated, and they can move ±90 degrees in both the pitch and the yaw axes.  To evaluate the dexterity and manipulability of the proposed manipulator, a confined semi-structured environment was created, as shown in Figure 13. It is a different environment from the environment used in the optimisation process to show that the manipulator, based on the optimisation results, can perform effectively in any confined semistructured environment. The dimensions of the environment are 1.2 m in the x-axis, 0.6 m To evaluate the dexterity and manipulability of the proposed manipulator, a confined semi-structured environment was created, as shown in Figure 13. It is a different environment from the environment used in the optimisation process to show that the manipulator, based on the optimisation results, can perform effectively in any confined semi-structured environment. The dimensions of the environment are 1.  To evaluate the dexterity and manipulability of the proposed manipulator, a confined semi-structured environment was created, as shown in Figure 13. It is a different environment from the environment used in the optimisation process to show that the manipulator, based on the optimisation results, can perform effectively in any confined semistructured environment. The dimensions of the environment are 1.  Next, a collision-free trajectory was generated to drive the manipulator from the home configuration, Figure 13  Next, a collision-free trajectory was generated to drive the manipulator from the home configuration, Figure 13 Figure 15. Considering the restrictions of the confined workspace, the joint positions did not go to their maximum limit. The highest joint position was for joint one with 52.4 degrees, followed by joints nine, seven and eight with 40.7, 35.5 and −34.4 degrees, respectively, while joints four and eleven did not go above −28 degrees. Seven joints (two, six, ten, thirteen, fourteen, fifteen and sixteen) did not move at all, while joints three, five and twelve remained at almost zero throughout the trajectory. Considering the various obstacles in the given environment, the manipulator was able to move around obstacles and reach its goal without using the ±90-degrees full range of angles. Thus, implementing the manipulators resulting from the optimisation used in this study in environments with a different set of obstacles has been shown to result in a flexible performance given the high manipulability results.
(two, six, ten, thirteen, fourteen, fifteen and sixteen) did not move at all, while joints three, five and twelve remained at almost zero throughout the trajectory. Considering the various obstacles in the given environment, the manipulator was able to move around obstacles and reach its goal without using the ±90-degrees full range of angles. Thus, implementing the manipulators resulting from the optimisation used in this study in environments with a different set of obstacles has been shown to result in a flexible performance given the high manipulability results.   grees, respectively, while joints four and eleven did not go above -28 degrees. Seven joints (two, six, ten, thirteen, fourteen, fifteen and sixteen) did not move at all, while joints three, five and twelve remained at almost zero throughout the trajectory. Considering the various obstacles in the given environment, the manipulator was able to move around obstacles and reach its goal without using the ±90-degrees full range of angles. Thus, implementing the manipulators resulting from the optimisation used in this study in environments with a different set of obstacles has been shown to result in a flexible performance given the high manipulability results.    Figure 16 shows the velocity ellipsoids for the manipulator, with the design based on the optimisation results, at the middle and final pose configurations. The figure shows that the proposed manipulator has a high manipulability as the velocity ellipsoids are both isotropic (have a spherical shape). This implies that the manipulator can change its end-effector positions in all axes, using different configurations of joints, despite the restrictions imposed by the confined workspace. Therefore, the manipulator, with its design being based on the optimisation results, can move freely and generate high velocity in the x, y and z axes, regardless of the shape of the obstacles. the optimisation results, at the middle and final pose configurations. The figure shows that the proposed manipulator has a high manipulability as the velocity ellipsoids are both isotropic (have a spherical shape). This implies that the manipulator can change its end-effector positions in all axes, using different configurations of joints, despite the restrictions imposed by the confined workspace. Therefore, the manipulator, with its design being based on the optimisation results, can move freely and generate high velocity in the x, y and z axes, regardless of the shape of the obstacles.

Conclusions
This paper proposed a concept for structural optimisation of hyper-redundant manipulators that can be used in confined semi-structured environments, such as in agriculture. Two methods of optimisation were used, using the Fmincon and Artificial Bees Colony (ABC) algorithms, considering the manipulability measures and inverse kinematics, to ensure optimal joint parameters that enable the manipulator to work flexibly within its task space.
This optimisation method was conducted for manipulators with universal joints, as the universal joints provide a wide range of motion for the manipulator to cover two axes with a 360-degree workspace, so all points within the manipulator workspace are reachable. The optimised manipulator was assumed to have a 1 m length, and the optimisation method was conducted to ensure high performance within its workspace, so the reachability of the points within its workspace with the optimised optimal lengths of the links was shown to be successful. Therefore, for long manipulators, an additional number of links alone need to be added to the manipulator for it to cover its workspace dexterously.
A hyper-redundant manipulator was designed based on the results of optimisation to be used for validation in a confined semi-structured environment. Its performance was tested by assessing its manipulability in certain confined locations. It shows a high level of dexterity using different joint configurations, as it is depicted by manipulability ellipsoids. Unlike previous studies that only optimised conventional manipulators, which considered kinematic constraints alone, the optimisation in this study depends on the manipulability measures and inverse kinematics in a confined workspace. This shows the manipulator is dexterous and performs with high manipulability within any given workspace, as is shown as the testing environment was completely different from that in the optimisation scenario. This is demonstrated by the results of the validation, as the manipulator was able to reach its goal and move around the different obstacles with a range of

Conclusions
This paper proposed a concept for structural optimisation of hyper-redundant manipulators that can be used in confined semi-structured environments, such as in agriculture. Two methods of optimisation were used, using the Fmincon and Artificial Bees Colony (ABC) algorithms, considering the manipulability measures and inverse kinematics, to ensure optimal joint parameters that enable the manipulator to work flexibly within its task space.
This optimisation method was conducted for manipulators with universal joints, as the universal joints provide a wide range of motion for the manipulator to cover two axes with a 360-degree workspace, so all points within the manipulator workspace are reachable. The optimised manipulator was assumed to have a 1 m length, and the optimisation method was conducted to ensure high performance within its workspace, so the reachability of the points within its workspace with the optimised optimal lengths of the links was shown to be successful. Therefore, for long manipulators, an additional number of links alone need to be added to the manipulator for it to cover its workspace dexterously.
A hyper-redundant manipulator was designed based on the results of optimisation to be used for validation in a confined semi-structured environment. Its performance was tested by assessing its manipulability in certain confined locations. It shows a high level of dexterity using different joint configurations, as it is depicted by manipulability ellipsoids. Unlike previous studies that only optimised conventional manipulators, which considered kinematic constraints alone, the optimisation in this study depends on the manipulability measures and inverse kinematics in a confined workspace. This shows the manipulator is dexterous and performs with high manipulability within any given workspace, as is shown as the testing environment was completely different from that in the optimisation scenario. This is demonstrated by the results of the validation, as the manipulator was able to reach its goal and move around the different obstacles with a range of angles between ±60 degrees, which proves the high level of flexibility of the optimised manipulator in handling its workspace. Therefore, more complicated workspaces with different obstacles can be handled within the optimised manipulator as its range of angles is ±90 degrees. This means the manipulator can perform flexibly within any given confined semi-structured environment or workspace. To the best of our knowledge, this is the first work that optimises hyper-redundant manipulators in semi-structured environments. The previous works only consider conventional manipulators, and dexterity and manipulability measures have been neglected in their optimisation. As a future extension of this study, various control techniques will be investigated to ensure the accuracy and robustness of the manipulator's performance in semi-structured environments.