Design, Implementation, and Kinematics of a Twisting Robot Continuum Arm Inspired by Human Forearm Movements

In this article, a soft robot arm that has the ability to twist in two directions is designed. This continuum arm is inspired by the twisting movements of the human upper limb. In this novel continuum arm, two contractor pneumatic muscle actuators (PMA) are used in parallel, and a selfbending contraction actuator (SBCA) is laid between them to establish the twisting movement. The proposed soft robot arm has additional features, such as the ability to contract and bend in multiple directions. The kinematics for the proposed arm is presented to describe the position of the distal end centre according to the dimensions and positions of the actuators and the bending angle of the SBCA in different pressurized conditions. Then, the rotation behaviour is controlled by a high precision controller system.


Introduction
Biological inspiration leads to an extensive number of inventions around the world. A rigid robot is just one example of an invention that has benefited us, either in industrial applications or in our social lives. However, using these machines in close proximity to humans has dramatically increased during the last decade, which has led to an increased risk of injury. Soft robots are a good alternative to rigid robots when it comes to health and safety due to their impressive advantages.
The softness of the material that is typically used to build soft robotics provides the term 'soft' and is normally used to cover the entire robot body without any rigid joints. One widely-used soft actuator is the pneumatic muscle actuator (PMA), which offers various benefits, such as being lightweight, having the variety and flexibility to be built into different structures [1][2][3][4][5], having high degrees of freedom (DoF) [6][7][8], safe to individuals and can be implemented by various dimensions [2,5].
Despite the general similarity in soft PMAs, each design can show a unique performance, one of these prototypes is the OctArm by [14], which shows adaptability, flexibility and efficient performance. Nevertheless, it is not easy to build and control as it is constructed from many actuators. This design shows mechanical difficulties. The Air-Octor soft arm in [15] is noticeably less complicated to build and control because the body of the arm is a single soft PMA and it is actuated by tendons. Though, the friction of the cables decreases its performance, which produces unpredicted changes in cable binding and moves the trunk in indefinite performance. A soft arm of three sections of six extensor PMA in the first two sections and three extensor actuators in the Section 3 is proposed considered a fixed bone, while the radius rotates longitudinally about its shaft from the elbow and in semi-circular rotation from the wrist joint.
Twisting toward a clockwise direction is called pronation and the opposite rotation is called supination. The diagonal intersection between radius and ulna is occurring during the pronation as shown in Figure 1. A small fibrous connects the two bones at halfway between the wrist and the elbow and it is studied by academics as an extra joint [44]. The twisting behaviour of the two bones is illustrated in Figure 2. The bones are also linked by an interosseous layer to enlarge the region of the muscular connection along with maintaining a specific correlation between the radius and the ulna.
Robotics 2022, 11, x FOR PEER REVIEW 3 of 23 considered a fixed bone, while the radius rotates longitudinally about its shaft from the elbow and in semi-circular rotation from the wrist joint. Twisting toward a clockwise direction is called pronation and the opposite rotation is called supination. The diagonal intersection between radius and ulna is occurring during the pronation as shown in Figure 1. A small fibrous connects the two bones at halfway between the wrist and the elbow and it is studied by academics as an extra joint [44]. The twisting behaviour of the two bones is illustrated in Figure 2. The bones are also linked by an interosseous layer to enlarge the region of the muscular connection along with maintaining a specific correlation between the radius and the ulna.

The Design of a Twisting Soft Robot Arm
The human arm-twisting behaviour due to an elbow joint that is described in the previous section has inspired us to design a soft robot arm that can be twisted. The structure of the proposed arm is given in Figure 3. A self-bending contraction actuator (SBCA)

The Design of a Twisting Soft Robot Arm
The human arm-twisting behaviour due to an elbow joint that is described in the previous section has inspired us to design a soft robot arm that can be twisted. The structure of the proposed arm is given in Figure 3. A self-bending contraction actuator (SBCA) is implemented to establish the proposed soft arm's rotating behaviour. The bending actuator is constructed by adding a reinforcement flexible rod of 20 cm long, 0.6 cm width, and 0.2 cm thickness to the 20 cm contraction PMA. The rod is sewn to the braided sleeve to prevent the contraction from that side, as shown in Figure 3a.
When the actuator is pressurised, the addition of the reinforcing rod means that one side of the muscle will be prevented from contracting and remains at its maximum length. This difference in contraction on the two sides of the actuator results in the actuator bending. Further details of the operation of the SBCA can be found in [20].

The Kinematics Models of the Twisting Soft Robot Arm
The key behaviour of the proposed soft robot continuum arm is the ability to twist anticlockwise and return clockwise by about 91°, while the human elbow produces about 70° in both directions [45]. However, unlike the human arm, where the radius and ulna are rigid, in this system, they are both formed from flexible actuators, which can be adjusted in length. This means that, in addition to the supination-pronation motion, the arm can also bend as the relative lengths of the radius and ulna are adjusted. The proposed continuum arm provides a wide range of movements according to the quantity of air pressure in every one of the three muscles. These movements represent extra characteristics of the suggested robot arm. Pressurising both contraction actuators at a similar amount leads to shortness of the continuum arm. The maximum safe work pressure is 500 kPa and the maximum contraction ratio for both contraction PMAs is 30% (11 cm). While applied Z Y X Two contraction actuators of 37 cm in length are used in parallel and one SBCA actuator of 20 cm is positioned between them to establish the twisting. The contraction actuators are made similarly by using a 42 cm inner rubber tube and 42 cm braided mesh and two solid ends of 2.5 cm each. After construction, the unpressurised length is 37 cm (see Figure 3b). Like the ulna bone of the person's forearm, the right contraction PMA is fixed at the elbow and wrist joints as shown in Figure 3c. the other air muscle is laid in parallel but has been connected to the encoder from the top and the rotating joint from the bottom to twist along its upright axis from the upper, and over the right actuator from the base. This operation behaves similar to the radius of a human arm. To test the twisting angle of the proposed continuum arm, air pressure is applied via the solenoid valve on the SBCA, and the twisting angle is observed by the encoder. The experiment gives that pressurising both contraction actuators simultaneously at any value shows no effect on the twisting angle. The range of the rotating angle is 0 • to 91 • . As a result, the supination behaviour is occurring in a pressurised condition and the pronation occurs during the venting process.

The Kinematics Models of the Twisting Soft Robot Arm
The key behaviour of the proposed soft robot continuum arm is the ability to twist anticlockwise and return clockwise by about 91 • , while the human elbow produces about 70 • in both directions [45]. However, unlike the human arm, where the radius and ulna are rigid, in this system, they are both formed from flexible actuators, which can be adjusted in length. This means that, in addition to the supination-pronation motion, the arm can also bend as the relative lengths of the radius and ulna are adjusted. The proposed continuum arm provides a wide range of movements according to the quantity of air pressure in every one of the three muscles. These movements represent extra characteristics of the suggested robot arm. Pressurising both contraction actuators at a similar amount leads to shortness of the continuum arm. The maximum safe work pressure is 500 kPa and the maximum contraction ratio for both contraction PMAs is 30% (11 cm). While applied different air pressure on the contraction actuators shows bending behaviour in the right or left direction, the arm bends toward the actuator of high pressure. The SBCA provides the rotation behaviour and shows multi-degree of freedom (DoF) movements with a combination of different air pressure in the contraction PMAs.
The position and orientation for the centre of the free end change under the impact of the different pressures in the three actuators. Figure 4 shows the proposed soft robot arm under different pressurising conditions. The position and orientation for the centre of the free end change under the impact of the different pressures in the three actuators. Figure 4 shows the proposed soft robot arm under different pressurising conditions.
To analyse the position and the orientation of the centre point of the distal end, the kinematics for the proposed arm is discussed from three views: Y-Z, X-Y and X-Z planes. Figure 5a shows the front view of the proposed arm in non-pressurised conditions. Where the length of the left and right contraction actuators at no pressure is defined as L01 and L02 respectively, r1 is the length that lay between the left contraction PMA and the centre of the soft arm, and r2 is the distance between the right contraction PMA and the rotating centre. L1 and L2 are lengths under pressurised conditions. Figure 5b-d illustrates the Y-Z, X-Y and X-Z views respectively, at pressurised conditions of the SBCA. l1 and l2 are the linear fronts and side view of L1 and L2 respectively, λ1 is the resultant side distance between the initial point of the first contraction actuator at relaxed conditions of the SBCA and the new location of the first contraction air muscle at pressurised conditions of the SBCA, γ is the deviation angle of the left contraction PMA due to the bending of the SBCA, r1 and r2 are the length between the spinning centre and the centre To analyse the position and the orientation of the centre point of the distal end, the kinematics for the proposed arm is discussed from three views: Y-Z, X-Y and X-Z planes. Figure 5a shows the front view of the proposed arm in non-pressurised conditions. Where the length of the left and right contraction actuators at no pressure is defined as L 01 and L 02 respectively, r 1 is the length that lay between the left contraction PMA and the centre of the soft arm, and r 2 is the distance between the right contraction PMA and the rotating centre. L 1 and L 2 are lengths under pressurised conditions. Figure 5b-d illustrates the Y-Z, X-Y and X-Z views respectively, at pressurised conditions of the SBCA. l 1 and l 2 are the linear fronts and side view of L 1 and L 2 respectively, λ 1 is the resultant side distance between the initial point of the first contraction actuator at relaxed conditions of the SBCA and the new location of the first contraction air muscle at pressurised conditions of the SBCA, γ is the deviation angle of the left contraction PMA due to the bending of the SBCA, r 1 and r 2 are the length between the spinning centre and the centre of the contraction air muscles, and β represents the bending position of the bending actuator. The contraction ratio ε of the contraction actuators is defined as: Assuming both contraction actuators are identical, then: From Figure 5b: And: From Figure 5c: Then: And from Figure 5e: φ is the angle between the second contraction muscle and the vertical layout at pressurised conditions. Applying air pressure to the contraction actuators causes a bending angle of ϕ degrees in the direction of the shortest contraction PMA, as illustrated in Figure 6.
∆Z and ∆Y are the displacements in both Z and Y directions. from the set of equations above, it is easy to find the new position and angle of the rotating centre.
The effect of the deformation of the three actuators on the kinematics of the proposed twisting arm can be ignored due to the rotation ability of the arm and the stiffness of the actuators at pressurised conditions. Δ = − cos (10) ΔZ and ΔY are the displacements in both Z and Y directions. from the set of equations above, it is easy to find the new position and angle of the rotating centre.
The effect of the deformation of the three actuators on the kinematics of the proposed twisting arm can be ignored due to the rotation ability of the arm and the stiffness of the actuators at pressurised conditions.

The Validation and the Experiment Results
In this section, air pressure is applied to the contraction actuators and the SBCA via a (3/3 Matrix) solenoid valve and the pressure in each actuator is measured by an analogue pressure sensor. Since the right contraction PMA is fixed from both sides, the MPU 6050 sensor is fixed to the end of the left contraction PMA to measure the γ angle. The bending angle of the SBCA β is measured by the encoder, which is fixed instead to the upper joint (see Figure 4). Several experiments are done to illustrate the behaviour of the presented twisting soft robot arm. The pressure is increased in the SBCA from 0 to 400 kPa by steps of 50 kPa, and the angle β is measured experimentally. Its data is listed in Table 1. The angles γ, ϕ, and φ are measured experimentally by the MPU sensors and calculated from Equations (6), (7), and (8), respectively.   Table 1 at zero air pressure on both contraction actuators. The model results show significant parallels with the experimental data, with small errors, high correlation factors, and similar performances. Despite what the experiment shows when both the contraction actuators are in the relaxed condition, the rotating centre shows a displacement varying in Y and Z directions according to (9) and (10), and are listed in Figure 8.

The Validation and the Experiment Results
In this section, air pressure is applied to the contraction actuators and the SBCA via a (3/3 Matrix) solenoid valve and the pressure in each actuator is measured by an analogue pressure sensor. Since the right contraction PMA is fixed from both sides, the MPU 6050 sensor is fixed to the end of the left contraction PMA to measure the γ angle. The bending angle of the SBCA β is measured by the encoder, which is fixed instead to the upper joint (see Figure 4). Several experiments are done to illustrate the behaviour of the presented twisting soft robot arm. The pressure is increased in the SBCA from 0 to 400 kPa by steps of 50 kPa, and the angle β is measured experimentally. Its data is listed in Table 1. The angles γ, φ, and ϕ are measured experimentally by the MPU sensors and calculated from Equations (6), (7), and (8), respectively.  Figure 7a-c shows the orientation angles of the rotating centre as a function of the bending angle of the SBCA according to the data in Table 1 at zero air pressure on both contraction actuators. The model results show significant parallels with the experimental data, with small errors, high correlation factors, and similar performances. Despite what the experiment shows when both the contraction actuators are in the relaxed condition, the rotating centre shows a displacement varying in Y and Z directions according to (9) and (10), and are listed in Figure 8.  Applying air pressure to both the contraction PMAs affects all of the parameters (see Equations (6)-(10)), and to demonstrate this impact, air pressure is applied to the left and the right contraction actuators individually. The length variations for both the contraction actuators are listed in Table 2. The length variation between the two contractor actuators is ignored due to the similarity in design and implementation.   Applying air pressure to both the contraction PMAs affects all of the parameters (see Equations (6)-(10)), and to demonstrate this impact, air pressure is applied to the left and the right contraction actuators individually. The length variations for both the contraction actuators are listed in Table 2. The length variation between the two contractor actuators is ignored due to the similarity in design and implementation.  Applying air pressure to both the contraction PMAs affects all of the parameters (see Equations (6)-(10)), and to demonstrate this impact, air pressure is applied to the left and the right contraction actuators individually. The length variations for both the contraction actuators are listed in Table 2. The length variation between the two contractor actuators is ignored due to the similarity in design and implementation. Three different β values are chosen to demonstrate the effects of L 1 on the presented continuum arm parameters at relaxed conditions for the right contraction actuator. Figure 9a-d illustrates these effects at β equal to (20 • , 40 • , and 60 • ) for the calculated and measured data. The angles except for β in these figures have been measured manually three times and the average has been drawn.
The length change of the contraction actuators is defined as: where: s is the number of actuators. Since both the contraction actuators have a similar initial length, then: And Figure 9a shows the behaviour of the γ angle with the raised change in length of the right contraction PMA at three different values of β. This figure illustrates that the increase in ∆L 1 and/or β leads to an increase in γ, and that is a validation of (6) and Figure 7a. Similarly, the increment in ϕ is due to the increase in ∆L 1 and/or β (see (8) and Figure 7c). The variations of the rotating centre are similarly increased in the Y and Z directions.
The experiment is repeated but at the relaxed pressurising conditions of L 1 , and the air pressure is increased in the right contraction actuator according to Table 2. The experiment results are used to calculate φ from (7) at three different values of β, then, the angle ϕ and the changes in the Y and Z positions of the rotating centre are found for (L 1 = 37 cm). The calculating results are found according to (6)-(10) and the measured data is found manually three times, then the average is recorded. Figure 10a-d demonstrates the parameters of the twisting soft arm at different pressurised conditions of the right contraction actuator. Figure 10a shows that the increase in ∆L 2 and/or β leads to a decrease in the φ angle at a fixed length of the left actuator, according to (7). Otherwise, the increment in one or both of ∆L 2 and β leads to a decrease in the value of ϕ and makes the orientation around X of the rotating centre negative. The comparison between Figures 9b and 10b illustrates that the ϕ values in a positive direction are more than the values in the negative direction because the rise in β leads to the rise in ϕ in the positive direction (see Figure 7c). On the other hand, the variation of Y and Z is also raised with the increase in ∆L 2 and/or β. Nevertheless, the variation in Z, due to the change in L 2 , is more than the variation due to the change in L 1 , because of the effects of the increment of β (see Figure 8).
The experiment results due to the variety in both L 1 and L 2 are validated by comparing the results in Figures 9 and 10 and listed in Table 3 in the terms of coefficient of determination (R 2 ).
On the other hand, an experiment to observe the load effect on the contraction and twisting angle is done by attaching a load to the end of the proposed arm. The load is varied from zero to 2 kg. As shown in Table 4, the length of the contraction actuators is slightly affected by increasing the attached load, while the variations in twisting angle can be ignored.
In a comparison with the literature, the proposed arm mimics the performances of the human arm at high force. The twisting arm of [21] provides a high twisting angle but at 0.4 N, 2.5 N for the robot arm in [23], and 5 N in [24]. While the contraction PMA provides about 400 N [46] Robotics 2022, 11, x FOR PEER REVIEW 13 of 23 In a comparison with the literature, the proposed arm mimics the performances of the human arm at high force. The twisting arm of [21] provides a high twisting angle but at 0.4 N, 2.5 N for the robot arm in [23], and 5 N in [24]. While the contraction PMA provides about 400 N [46].

Control the Twisting Angle
The encoder has been used instead of the upper joint in the rotating continuum arm. The encoder sends the rotating angle as feedback, which is proportional to the resistance value to the controller input as shown in Figure 11. The presented controller system, parallel neural network and proportional (PNNP), in the previous article [47] have been used to control the twisting angle. The training process is done by (trainlm) (MATLAB 2020a) for 100 iterations and at 10 −7 mean square error (MSE).

Control the Twisting Angle
The encoder has been used instead of the upper joint in the rotating continuum arm. The encoder sends the rotating angle as feedback, which is proportional to the resistance value to the controller input as shown in Figure 11. The presented controller system, parallel neural network and proportional (PNNP), in the previous article [47] have been used to control the twisting angle. The training process is done by (trainlm) (MATLAB 2020a) for 100 iterations and at 10 −7 mean square error (MSE). Figure 11. The block diagram of the controller system.
The controller system is built by a parallel structure of the NARMA-L2 Neural Network (NN) control system and proportional (P) controller. The NN ensures accuracy, while the P controller increases the system response [47].
The neural network control system has been trained under the approximation model. The twisting angle of the continuum arm is given by (14) as a function of the duty cycle.  Figure 11. The block diagram of the controller system.
The controller system is built by a parallel structure of the NARMA-L2 Neural Network (NN) control system and proportional (P) controller. The NN ensures accuracy, while the P controller increases the system response [47].
The neural network control system has been trained under the approximation model. The twisting angle of the continuum arm is given by (14) as a function of the duty cycle.
where y refers to the rotating angle and u * represents the duty cycle of the pulse width modulation (PWM) signal. The number (90) represents the maximum rotating angle, the 98% is set to be the maximum duty cycle instead of 100% to protect the solenoid valve from the continuous operation. (14) is found by applying different values to the duty cycle for one second and then recording the rotating angle individually by the potentiometer. The controller outputs for both the neural networks controller and the proportional controller are given below: The NARMA-L2 NN-controller output u1 can be defined as: where f ( ) and g( ) are approximated using neural networks. y( ) is the system output and: where m = 3 and n = 2 and they refer to the 3-delayed plant inputs, and 2-delayed plant outputs, respectively. On the other hand, the output of the proportional controller is given by: x is the maximum rotation angle (90 • in this case). And the overall controlled PWM signal is: Filling and venting the air muscle are two different processes. Therefore, two PNNP controllers are required, one controller for each of them. The control system calculates the error between the reference model (desired rotating angle) and the feedback. Positive error activates the fill controller system, while negative error initiates the vent system. Both branches work together in parallel to gain zero error. Applying sinusoidal validates the controller and rectangular pulse set signals at 0.25 and 0.5 Hz, as shown in Figure 12. Figure 12 shows the accuracy of the controller system to be used for various rotating angles. The sinusoidal response shows a low constant error because of the small changes in the signal over time.
On the other hand, the rectangular pulse response has a zero steady-state error because of its unchanged values at zero and 75 • . The pressure of the two contractors' PMAs does not affect the rotating feature. However, increasing the stiffness of the contraction actuators make them act like bones, similarly, the rectangular pulse response shows that the release time for the SBCA is more than the pressurised period for two reasons: the actuator hysteresis and the pressure difference between the PMA and outside environment. The static and dynamic characteristics of the solenoid valve and the pressure dynamics were ignored in the control process.

Conclusions
The rotational performance of the individual forearm due to the elbow wrist

Conclusions
The rotational performance of the individual forearm due to the elbow wrist joints is described briefly in this article. This movement toward clockwise and counterclockwise has been used to design a soft continuum soft robot arm. The major element of the presented robot arm is the pronation and supination performances. In addition to this performance, the suggested continuum robot arm has the ability to bend and contract. These behaviours together with a twisting performance provide a high range of DoF. The novel arm has been constructed by two contraction PMAs in parallel and one SBCA, hence, the kinematics for the rotating behaviour are given under multiple actuation levels for the contraction actuators. To validate the performance of the presented robot arm, a control system is utilised to adjust the twisting angle.
AS illustrated through the article, the presented robot is easy to implement at a low cost, providing a high level of performance in terms of the degree of freedom, force to weight ratio, and robustness.
A potential future work, the presented continuum arm can be used in numerous industrial applications, and the twisting performance might be suitable for disabled people with an unhealthy forearm.

Data Availability Statement:
The datasets generated and/or analysed during the current study are available from the corresponding author on reasonable request. Acknowledgments: Many thanks to the University of Basrah and its computer engineering department to provide the lab facilities.

Conflicts of Interest:
The authors declare no conflict of interest.

Nomenclature ε
The contraction ratio L 0 The nominal length of the actuator [cm] L The length of the actuator [cm] r 1 , r 2 Distance between the rotating centre and contraction PMA [cm] ∆Z, ∆Y The displacement in both Z and Y directions [cm] y The output of the Plant system [Degree] u1 The controller output for the NN controller k p The proportional constant φ The angle between the second PMA and the vertical layout [Degree] u* The duty cycle of the controlled signal u2 The controller output for the P controller