Conﬁguration Design of an Upper Limb Rehabilitation Robot with a Generalized Shoulder Joint

: For stroke patients with upper limb motor dysfunction, rehabilitation training with the help of rehabilitation robots is a social development trend. Existing upper limb rehabilitation robots have difﬁculty fully ﬁtting the complex motion of the human shoulder joint and have poor human– robot compatibility. In this paper, based on the anatomical structure of the human upper limb, an equivalent mechanism model of the human upper limb is established. The conﬁguration synthesis of the upper limb rehabilitation mechanism was carried out, a variety of shoulder joint man–machine closed-chain Θ s and shoulder elbow human–machine closed-chain Θ se conﬁguration combinations were synthesized, and the conﬁguration model with compatibility and reduced moment conduction attenuation was selected from them. Two conﬁgurations, 2 P a 1 P 3 R a and 5 R a 1 P , are proposed for the generalized shoulder joint mechanism of the robot. The closed-chain kinematic models of the two conﬁgurations are established, and the velocity Jacobian matrix is obtained. Motion performance analysis, condition reciprocal analysis and operability ellipsoid analysis of different conﬁguration design schemes were carried out in different operation planes. The results show that in the normal upper limb posture of the human body, the 5 R a 1 P conﬁguration of the shoulder joint has better kinematic performance. Finally, on this basis, an upper limb rehabilitation robot prototype with good human–computer compatibility is developed, and its moving space was veriﬁed.


Introduction
Neuromuscular injury can lead to disability or inconvenient movements, such as stroke and spinal cord injury, which have become important problems worldwide [1]. Limb motor dysfunction is caused by a wide range of nervous system diseases, and both patients and communities bear a huge economic burden [2]. Currently, there are more than 33 million stroke patients in the world [3]. Robots play a key role in the development of rehabilitation technology [4].
The human arm can complete complex movements, such as grasping, lifting, touching the eyes and touching the back of the head, in a complex environment. This ability is due to the complex bones and abundant muscle groups of the arm [5]. Since the 21st century, some studies on the anatomy of the human shoulder joint have shown that the rotation axis of the human glenohumeral joint drifts when it moves [6,7]. The joint movements of the scapula, clavicle and humerus are called the shoulder humeral rhythm [8,9].
Because the axis drift of the glenohumeral joint only occurs in the range of motion of the shoulder joint, if the motion of the shoulder blade belt is considered, it increases the complexity of the shoulder joint of the rehabilitation robot, so at present, most upper limb rehabilitation robots widely use six degrees of freedom. The representative product is Armeo Power developed by HOCOMA [10,11]. On this basis, HOCOMA developed the Armeo Spring training evaluation device [12], which is currently the most commonly used upper limb rehabilitation robot in clinics. Additionally, the 7-DOF (Degree-of-Freedom) upper limb rehabilitation robot, CADEN-7, was developed by the University of Washington [13]. The team of Pro. Pignolo in Italy developed the ARAMIS dual arm rehabilitation robot system [14], and the representative product in China is the SFRobot developed by Anyang Co. Shen Fang.
However, 6-DOF upper limb rehabilitation robots ignore the movement of the scapula, and the joint axis does not match the movement of the upper limbs of the human body, which may lead to a joint compensation function of the patient and affect the reconstruction of the normal movement function of the upper limb of the patient.
In research results involving the movement of the scapula, various types of mechanism configurations have been proposed to realize the drift of the axis of the shoulder joint movement. Canada Queen's University designed MEDARM [15], which has a 5-DOF shoulder joint composite structure. In MEDARM, the 2-DOF sternoclavicular joint and the 3-DOF shoulder joint can track the movement of the robot's shoulder joint rotation to the human body's generalized shoulder joint rotation. The dual-arm rehabilitation robot Harmony [16] developed by the University of Texas adds 2-DOF of shoulder girdle raising/lowering and extension/retraction in the shoulder joint and can adjust the swing centre of the boom on a spherical surface. LIMPACT [17], developed by the University of Twente in the Netherlands, adds two passive DOFs at the shoulder joint of the mechanism, which can realize movement of the centre of the ball and socket joint at any position on the horizontal plane. The Intelli-Arm [18], developed by Northwestern University in the United States, adds three mutually perpendicular mobile joints to the shoulder joints of the organization; of these joints, two are passive and one is active, and they can adapt to changes in the scapulohumeral rhythm. Others, including Arizona State University [19], Columbia University [20], Wonkwang University [21], Southeast University [22], North University of China [23], Beijing University of Technology [24], Shanghai JiaoTong University [25,26], Hubei University of Technology [27], etc., have conducted research in this area.
In existing research, the mechanical configuration of shoulder rehabilitation has been more in-depth, and most of the upper limb rehabilitation robots with bionic shoulder strap structure are relatively large, which is due to the introduction of the driving scapula structure. From the perspective of kinetic energy and energy consumption, there is a great challenge for the development of a wearable robot.
In conclusion, the rehabilitation robot considering the axis of the human glenohumeral joint as the floating point is a future development direction, and it can avoid the uncertain binding force/torque at the human-machine constraint caused by the axis floating of the human glenohumeral joint in the training process, which can cause secondary injury to patients.
In this paper, the anatomical and physiological structure of the shoulder joint is analysed, and the equivalent configuration of the upper limb is established. Different numbers and forms of passive joints are introduced to form a man-machine closed chain between the robot and upper limb of the human body. On this basis, the configuration scheme of the exoskeleton mechanism in line with the generalized motion law of the human shoulder joint is established. The kinematic performance of different design schemes is analysed and compared, and the optimal man-machine design scheme is proposed, which provides a theoretical reference and practical application value for the development of upper limb rehabilitation robot prototypes in the future.

Human Upper Limb Anatomical Structure and Equivalent Mechanism Model
Broadly speaking, the shoulder joint is complex and involves the glenohumeral joint, acromioclavicular joint, sternoclavicular joint and scapular-chest wall joints [28] as shown in Figure 1. These joints coordinate and cooperate to complete a complex and flexible shoulder movement. The glenohumeral joint (GH) is the main joint in the shoulder joint. The GH is composed of the ball and socket joint between the humeral head and the scapular glenoid [29]. The acromioclavicular joint (AC) is composed of the acromion end surface of the scapula and the end surface of the clavicle acromion joint [29]. The sternoclavicular joint (SC) is composed of the clavicle end face of the sternum and the clavicle notch of the sternum stem [30]. The scapulothoracic joint (ST) is a junction point between the front of the scapula and the back of the rib cage. The motion form of each subdivision joint of the scapula belt is shown in Figure 3. proposed, which provides a theoretical reference and practical application value for the development of upper limb rehabilitation robot prototypes in the future.

Human Upper Limb Anatomical Structure and Equivalent Mechanism Model
Broadly speaking, the shoulder joint is complex and involves the glenohumeral joint, acromioclavicular joint, sternoclavicular joint and scapular-chest wall joints [28] as shown in Figure 1. These joints coordinate and cooperate to complete a complex and flexible shoulder movement. The glenohumeral joint (GH) is the main joint in the shoulder joint. The GH is composed of the ball and socket joint between the humeral head and the scapular glenoid [29]. The acromioclavicular joint (AC) is composed of the acromion end surface of the scapula and the end surface of the clavicle acromion joint [29]. The sternoclavicular joint (SC) is composed of the clavicle end face of the sternum and the clavicle notch of the sternum stem [30]. The scapulothoracic joint (ST) is a junction point between the front of the scapula and the back of the rib cage. The motion form of each subdivision joint of the scapula belt is shown in Figure 2. Each joint in the shoulder joint is equivalent to the kinematic pair in the mechanism, in which the AC has 3-DOF and is equivalent to the ball pair S4. The SC has 3-DOF; however, its backward rotation is small, and it is less involved in the overall joint movement. Therefore, the SC is simplified to form Hooke's hinge with 2 DOF, R5, R6. The ST has 3-DOF, which can be simplified to two translation pairs, P1, P2 and a rotating pair, R3. The GH is a typical ball and socket joint, which is equivalent to a ball joint, S7. Therefore, the equivalent mechanism model of the anatomical structure of the scapula belt is shown in Figure 3. According to the Grubler-Kutzbach formula, the DOF of the parallel mechanism can be calculated as shown below.
where F is the DOF of the mechanism; λ is the number of DOFs the mechanism rods should have in space, generally 6; n is the number of mechanism rods; g is the number of motion pairs in the mechanism; fi is the number of DOFs of the i-th motion pair. Since the translation trajectory of P1, P2 is an arc trajectory centred on the axis of rotation of R5 and R6, R5 and R6 can be regarded as a virtual constraint of the translation pair of P1, P2. Since the R3 rotation axis is very close to the centre of the GH ball pair S7, it is usually regarded as the same position, so the R3 rotation effect overlaps the S7 motion range. Biologically, because GH S7 has a certain range of motion, scapula uprotation/downrotation R3 can be used to compensate for the range of motion of the GH. Therefore, the equivalent mechanism model of the scapula belt can be simplified to a 5-DOF tandem mechanism, as shown in Figure 4.  Each joint in the shoulder joint is equivalent to the kinematic pair in the mechanism, in which the AC has 3-DOF and is equivalent to the ball pair S 4 . The SC has 3-DOF; however, its backward rotation is small, and it is less involved in the overall joint movement. Therefore, the SC is simplified to form Hooke's hinge with 2 DOF, R 5 , R 6 . The ST has 3-DOF, which can be simplified to two translation pairs, P 1 , P 2 and a rotating pair, R 3 . The GH is a typical ball and socket joint, which is equivalent to a ball joint, S 7 . Therefore, the equivalent mechanism model of the anatomical structure of the scapula belt is shown in Figure 4. According to the Grubler-Kutzbach formula, the DOF of the parallel mechanism can be calculated as shown below.
where F is the DOF of the mechanism; λ is the number of DOFs the mechanism rods should have in space, generally 6; n is the number of mechanism rods; g is the number of motion pairs in the mechanism; f i is the number of DOFs of the i-th motion pair. Since the translation trajectory of P 1 , P 2 is an arc trajectory centred on the axis of rotation of R 5 and R 6 , R 5 and R 6 can be regarded as a virtual constraint of the translation pair of P 1 , P 2 . Since the R 3 rotation axis is very close to the centre of the GH ball pair S 7 , it is usually regarded as the same position, so the R 3 rotation effect overlaps the S 7 motion range. Biologically, because GH S 7 has a certain range of motion, scapula uprotation/downrotation R 3 can be used to compensate for the range of motion of the GH. Therefore, the equivalent mechanism model of the scapula belt can be simplified to a 5-DOF tandem mechanism, as shown in Figure 5.  In addition, upper limb exercises also include elbow joint flexion/extension, forearm pronation/supination, wrist joint dorsiflexion/palm flexion and radial flexion/ulnar flexion. Since the rotation axes of the 2-DOF of the wrist joint intersect, to reduce the design difficulty and production cost, the wrist joint of the rehabilitation robot adopts a single DOF, and the wrist joints are separately trained with different DOFs. Thus far, the equivalent mechanism model of the human upper limb can be established as a 9-DOF series motion model, as shown in Figure 6. To facilitate the measurement of upper limb kinematic information and unified clinical application communication, the International Society of Biomechanics (ISB) developed a description method for upper limb movement [31]. The lifting surface angle α, uplift angle within the uplift plane β and internal rotation/external rotation angle γ are used to describe the current position of the upper arm relative to the GH joint. Klopcar [32] et al. used experiments and synthetically obtained the kinematic coupling relationship between the SC and the GH rotation angle. The functional relationship between the uplift angle β and the forward/retracted angle ϕ fb and the uplift/sink angle ϕ ud of the SC is: Appl. Sci. 2021, 11, 2080 6 of 20

Human-Machine Compatible Configuration Design of the Shoulder Joint
In an ideal state, the human body's affected limb and mechanical arm form a closed kinematic chain [33]. Based on the Hunt formula [34], the DOF of the human-machine closed chain is calculated by Equation (4).
where F is the number of DOFs, n is the number of joints in the closed-chain system, f i is the number of DOFs of the joints, d is the dimension of the motion space of the member (generally taken as 6), and l is the number of closed loops.
In the human generalized shoulder joint, the AC, SC, ST and GH are simplified into one joint; the motion of the joint is three rotations; the drift of the GH movement axis is in space. According to the above, the corresponding mechanical arm has 5 DOFs at the shoulder joint. Then, a closed chain of the man-machine can be formed, including the closed chain of the shoulder joint Θ s . By analogy, the elbow joint closed-chain Θ e and the wrist joint closed-chain Θ w are shown in Figure 7. For mechanisms with unknown DOFs, Formula (4) is expressed as Formula (5).
where f k and f uk are the DOFs of the known joints and the unknown joints of the humanmachine closed chain, respectively, and j is the number of known joints in the closed chain. According to Figure 5, a simplified, equivalent mechanism model can be divided into three human-machine closed chains, namely, Θ s , Θ se and Θ sew . It is known that there are five DOFs of the human shoulder joint and five DOFs of the shoulder joint mechanism in the closed shoulder joint chain, so f k,s = 10. The same is true for f k,se = 13 and f k,sew = 18. The institutional design requirements are F s = 5, F se = 6 and F sew = 8. Therefore, the unknown joint DOFs f uk,s , f uk,se and f uk,sew can be solved by Equation (5).
From Equations (6) and (7), it can be seen that the shoulder joint closed loop chain introduces at least one passive DOF and the shoulder-elbow closed loop chain must introduce five passive DOFs, and Equation (8) shows that the overall human-machine closed chain must introduce eight passive DOFs. Due to the introduction of five or more passive degrees of freedom in the independent closed-loop chain, the passive drive chain is longer. The efficiency of energy and force transmission will be lower [35], the overall robot control will become particularly complicated, the passive DOF introduced in the closed chain of the manipulator wrist joint is first eliminated, and the spatial flexibility of the handle is reduced to achieve better manipulator control accuracy. Second, the human-machine closed chain of the shoulder joint Θ s is required to be the just-constrained subsystem, and the passive DOF introduced by the human-machine closed chain should be: According to the generalized shoulder joint motion mechanism model, the main motion of the GH will be realized by the RRR mechanism with three rotation axes converging at one point, and there are many combinations to meet the configuration of the auxiliary motion branch chain. The combination methods are shown in Table 1. Table 1. Configuration combination of the auxiliary movement branch chain. (P-movement pair, R-rotation pair, U-Hooke's hinge, C-Cylindrical pair, S-ball pair). The ball pair is difficult to process in engineering and cannot be used as an active drive to realize the auxiliary movement of the shoulder joint, and the range of the angle of movement is small, which cannot meet the range of movement of the upper limbs; for the convenience of engineering, Hooke hinge U and cylinder pair S are usually replaced by two rotating pairs R or one rotating pair R and moving pair P. Therefore, only the R or P combination mechanism configuration in the table is selected.

Kinematics Analysis of the Equivalent Mechanism
Referring to the configuration and joint distribution of mature arm II and harmony robots at home and abroad, and considering the factors that the passive branch chain is as short as possible and the motion load of active joint is as small as possible, 2P a 1P3R a and 5R a 1P were selected as mechanical arm shoulder joint mechanism configurations with engineering application value, where P a represents the active mobile pair, P represents the passive mobile pair and R a represents the active mobile pair.

2P a 1P3R a Shoulder Joint Mechanism Configuration
The configuration of 2P a 1P3R a is shown in Figure 8. In the figure, R 1 a4 , R 1 a5 and R 1 a6 are the three active rotation pairs for the main motion. The three rotational axes are orthogonal to each other. P 1 a1 and P 1 a2 are the two active moving pairs, and the moving axes are orthogonal to each other, which can realize the upward/sinking movement of the GH motion axis in the coronal plane of the human body. P 1 3 is the passive moving pair, which is orthogonal to P 1 a1 , P 1 a2 and each other, and the passive sliding axis can realize forward/backward movement of the GH motion axis in the sagittal plane of the human body and can combine with P 1 a1 , P 1 a2 to realize the tracking of the GH motion axis at any position in space, thereby improving the human-machine compatibility.
is established in the SC U SC , a human shoulder joint equivalent mechanism connection system O 1 s − x 1 s y 1 s z 1 s is established in the GH S GH , and a connection system O 1 g − x 1 g y 1 g z 1 g is established at the same time here. The coordinate axes of the two coordinate systems coincide in the initial pose. A local coordinate system O 1 w − x 1 w y 1 w z 1 w is established on the frame of the robot, and a joint coordinate system O 1 i − x 1 i y 1 i z 1 i (i = 1, 2, . . . , 6) is established for each movement pair of the generalized shoulder joint mechanism, where the R 1 a4 ,R 1 a5 and R 1 a6 motion axes intersect at the axis of the shoulder joint. In the initial configuration of the mechanical arm, the coordinate origin O 1 i (i = 4, 5, 6) coincides with the origin O 1 s of the GH joint coordinate system. The origin O 1 w of the gantry fixed coordinate system coincides with the movement centre O 1 1 of the mobile pair P 1 a1 , as shown in Figure 8. The origin O 1 w of the frame local fixed coordinate system in the human body fixed coordinate system According to the pose transformation matrix, the pose transformation matrix b T 1 6 of the end coordinate system O 1 6 − x 1 6 y 1 6 z 1 6 of the mechanical arm branch chain relative to the human body fixed coordinate system where b T 1 w , w T 1 1 and i T 1 i+1 (i = 1, 2, 3, 4, 5 are the pose transformation matrices between adjacent coordinate systems; n 1 6 = n 1 6x , n 1 6y , n 1 6z T , s 1 6 = s 1 6x , s 1 6y , s 1 6z T and a 1 6 = a 1 6x , a 1 6y , a 1 6z T are the direction vectors of each axis end coordinate system O 1 6 − x 1 6 y 1 6 z 1 T is the position vector of the origin O 1 6 . For the branch chain of the human shoulder joint equivalent mechanism, the pose transformation matrix b T 1 g of the mechanism end coordinate system O 1 g − x 1 g y 1 g z 1 g relative to the human body fixed coordinate system can also be obtained as: Since the end coordinate systems of the two branch chains are closed and the pose is always the same, Equations (10) and (11) are simultaneously solved, so the kinematics constraint equation of the man-machine closed chain can be established as: Deriving both sides of Equation (12) with respect to time t and sorting yields: are the rotation speeds of the three rotation pairs R 1 a4 , R 1 a5 and R 1 a6 in 2P a 1P3R a , respectively, and J 1 ij (i, j = 1, 2, 3) is the polynomial of the human body rotation angle parameters α, β, γ.

5R a 1P Shoulder Joint Mechanism Configuration
The configuration of 5R a 1P is shown in Figure 9. In the figure, R 2 a3 , R 2 a4 and R 2 a5 are the three active rotation pairs of the main motion, and the three rotation axes are orthogonal to each other. R 2 a1 , R 2 a2 are two active rotation pairs, and the rotation axes are orthogonal to each other, where R 2 2 , R 2 2 , R 2 2 and R 2 a2 constitute a parallelogram structure, which translates a rotating pair in the human body U SC to the back of the human body, which is convenient for installation and control in engineering. P 2 6 is the passive moving pair, which is directly connected to the upper limb.
is established on the SC U SC , a local fixed coordinate system O 2 w − x 2 w y 2 w z 2 w is established on the frame, and a connected coordinate system O 2 i − x 2 i y 2 i z 2 i (i = 1, 2, . . . , 6) is established for each movement pair of the generalized shoulder joint mechanism, where the R 2 a4 , R 2 a5 and R 2 a6 motion axes intersect the machine at the axis of the arm shoulder joint. In the initial configuration of the mechanical arm, the coordinate origin O 2 i (i = 4, 5, 6) coincides with the origin O 2 s of the GH joint coordinate system. The origin O 2 w of the gantry coordinate system coincides with the movement centre O 2 1 of the rotating pair R 2 a1 , as shown in Figure 9. The origin O 2 w of the frame local fixed coordinate system is X 2 w , Y 2 w , Z 2 w in the human body fixed coordinate system For the branch chain of the human shoulder joint equivalent mechanism, the end coordinate system of the mechanism can also be obtained, and the pose transformation matrix b T 2 g of O 2 g − x 2 g y 2 g z 2 g relative to the fixed coordinate system of the human body where b T 2 s is the pose transformation matrix between adjacent coordinate systems, and b T 2 g is the pose transformation matrix of the new coordinate system O 2 g − x 2 g y 2 g z 2 g after the end coordinate system of the mechanism is rotated in space. n 2 g = n 2 gx , n 2 gy , n 2 gz T , s 2 g = s 2 gx , s 2 gy , s 2 gz T and a 2 g = a 2 gx , a 2 gy , a 2 gz T are the direction vectors of each axis end coordinate system a 2 g = a 2 gx , a 2 gy , a 2 is the position vector of the origin O 2 g . Simultaneous to Equations (14) and (15), the kinematics constraint equation of the man-machine closed chain can be established as: Deriving both sides of Equation (16) with respect to time t and sorting yields:  in 5R a 1P, respectively, and J 2 ij (i, j = 1, 2, 3) are the terms of the human body's rotational angle parameters α, β, γ.

Inverse Kinematic Solution and Flexibility Analysis of Two Configurations
Referring to the size of the Chinese adult human body (GB/T10000-1988) [36] and combining the physical drive module size of two different configurations, the dimensions to be used are shown in Table 2. At the same time, the patient's upper arm lift angle range is set during the rehabilitation robot-assisted movement process as 0 •~1 50 • . The internal rotation/external rotation angle of the upper arm is set to γ = 30 • , and the upper arm's natural sag position is set as the lift angle β = 0 • . In other lifting planes, this configuration is also selected as the initial lifting angle to determine the initial movement posture in each lifting plane.
Combining the size parameters listed in Table 2 and the inverse solution method of human-machine closed-chain kinematics given above, when the upper arm of the human body is lifted in the 0 • , 30 • , 60 • , 90 • and 120 • lifting planes and the relationship between the turning angle θ 1 i (i = 4, 5, 6) of the 2P a 1P3R a turning pair R 1 i (i = 4, 5, 6) and the lifting angle β are shown in Figure 10b-f. The relationship between the moving pair R 2 i (i = 1, 2) of the 5R a 1P configuration and the lift angle β is shown in Figure 11a. The relationship between the turning angle θ 2 i (i = 3, 4, 5) of the turning pair R 2 i (i = 3, 4, 5) and the lifting angle β is shown in Figure 11b-f. It can be seen from Figure 11 that when the upper arm performs a lifting motion in different lifting planes, the rotation angles of each joint of the manipulator arm move smoothly, and the corresponding rotation angles move similarly in different lifting surfaces. The changes in θ 1 and θ 2 follow the piecewise function of the scapula belt and the lift angle, respectively, showing the same change law, and the rotation angle changes smoothly. Among them, θ 3 and θ 4 decrease with increasing β, while θ 5 first increases and then decreases. In different uplift planes, the changing trend and range of motion of θ 3 remain unchanged, which is closely related to the range of β. θ 4 continuously increases the initial change amplitude as the lift angle increases.
The singularity σ of the mechanical arm Jacobi qualitatively describes the dexterity and performance of the mechanical arm [37]. At present, there are many quantitative indexes related to Jacobian singular value to reflect the dexterity and kinematic performance of the manipulator. Rojas uses the reciprocal of the condition number of the Jacobian matrix as an evaluation index to judge the dexterity of the mechanism, which is faster and clearer. The relationship between the condition number and the singular value of the Jacobian matrix is: where σ min and σ max are the smallest and largest singular values of J (q) , respectively. The relationship curve between the reciprocal of the condition number K −1 J1 and the lifting angle β are compared when the two different configurations of the mechanism move in the same lifting surface, as shown in Figure 12. It can be seen that K −1 J1 of 2P a 1P3R a is less affected by the change in the lifting surface angle and first increases and then decreases as β increases. As the internal rotation/external rotation angle γ increases, the overall movement flexibility of 2P a 1P3R a improves. The configuration reaches a maximum value at β = 90 • and γ = 90 • , K −1 J1 = 1, which means that the mechanism has reached isotropy at this time.  Figure 13 shows the relationship between K −1 J2 of the Jacobian matrix of 5R a 1P and the change in the lift angle β. It is observed that in the same uplift surface, as β increases, K −1 J2 first increases and then decreases and then continues to increase. When β = 90 • , K −1 J2 reaches its maximum value and then falls. The curve trend shows that the mechanism has the best movement flexibility when β = 90 • . In the same uplift surface, comparing different internal rotation/external rotation angles γ, the change trend of K −1 J2 with β is similar, and the larger the γ, the overall value of K −1 J2 increases; however, when 60 • < β < 120 • , K −1

J2
following γ, the range of change is small.  Figure 14 shows the relationship between K −1 J2 and the uplift angle β and the uplift surface angle α for 5R a 1P under a specific internal rotation/external rotation angle γ = 60 • ,90 • . The change trend of K −1 J2 with β is similar to that in Figure 13. The change range of K −1

J2
following α is larger than that following γ. In the case of the same γ, an increase in α will cause an increase in K −1 J2 , which improves the movement flexibility of the mechanism. Comparing Figure 12 to Figure 14, the lift angle α = 90 • and the internal rotation/external rotation angle γ = 60 • are taken as the comprehensive performance evaluation poses of 2P a 1P3R a and 5R a 1P. Figure 15 shows the comparison of the reciprocal change curves of the condition number k −1 of the two manipulator shoulder joint configurations under the comprehensive performance evaluation position. As seen in Figure 15, when 40 • < β < 70 • , the movement flexibility of 2P a 1P3R a is better than that of 5R a 1P. When β is at other angles, the movement flexibility of 5R a 1P is better than that of 2P a 1P3R a , and K −1 J1 of 2P a 1P3R a does not exceed 0.6. Compared with 5R a 1P, the movement flexibility is at a disadvantage.

Operational Ellipsoid Analysis
The reciprocal of the condition number of the Jacobian matrix can be used to judge whether the current posture of the mechanical arm is singular. If the condition number is large, then the current Jacobian matrix is ill conditioned. However, for some actions, the reciprocal of the condition number is too large and does not affect the normal operation of the mechanism, so the concept of operability is proposed. Yoshikawa [38] defines operability: where λ i (i = 1,2,···n) is the eigenvalue of matrix JJ T , and σ i (i = 1,2,···n) is the singular value of Jacobian matrix J.
To compare the isotropic motion capabilities of the two mechanisms in different poses, the operability ellipsoids of the two mechanisms are compared at different positions on the coronal, sagittal and horizontal planes of the human body. Figures 16 and 17 show the operability ellipsoids with 2P a 1P3R a and 5R a 1P in the lift plane α = 30 • of the human body at positions β = 30 • , 60 • , 90 • , 120 • . Through comparison, when β = 90 • , the operable ellipsoids of the two configurations are closest to the sphere. At this time, 5R a 1P has better isotropic movement ability in all directions than 2P a 1P3R a .  In the lifting plane α = 30 • , when the angles change as β = 30 • , 60 • , 90 • , 120 • , the operability ellipsoids of the two configurations have similar changes in shape, indicating that the two have the same movement flexibility in α = 30 • , and the same-sex athletic ability is the same in the angular position. Figures 18 and 19 show the operability ellipsoids with 2P a 1P3R a and 5R a 1P, respectively, in the lift plane α = 90 • (sagittal plane) of the human body at position β = 30 • , 60 • , 90 • , 120 • . Through comparison, it can be seen that when β = 90 • , the operable ellipsoid of the two configurations is closest to the sphere, and the operability ellipsoid of 5R a 1P is very close to the sphere compared to 2P a 1P3R a . In this lift plane, the isotropic movement ability in all directions is better.  In the sagittal plane, when the lift angle changes at β = 30 • , 90 • , 120 • , the change in the operability ellipsoid of 5R a 1P is smaller than that of 2P a 1P3R a . When β = 60 • , the change in the operability ellipsoid of 5R a 1P changes more than 2P a 1P3R a . However, the overall performance shows that 5R a 1P has slightly better movement flexibility in the α = 30 • lifting plane and has slightly better single-sex movement ability under the corresponding lifting angle posture.
It can be seen from Figures 20 and 21 that 2P a 1P3R a has the same operability ellipsoid in the horizontal plane, indicating that the motion performance of the mechanism configuration corresponding to different lift angles on the horizontal plane remains unchanged. In the horizontal plane, when the lifting surface angle changes at α = 30 • , 90 • , 120 • , the overall operability of the5R a 1P ellipsoid is closer to a spherical shape than 2P a 1P3R a . The overall performance shows that 5R a 1P has better movement flexibility in the horizontal plane, and the same-sex sports ability is stronger in the corresponding lift angle pose.  The changes in the operability ellipsoids of the 2P a 1P3R a and 5R a 1P in four different forms on three planes are summarized in Table 3. A comparison index S is proposed. S is the D-value between the sum of the changes in the operability ellipsoids of two mechanisms. It can be seen from Table 3 that on the 30 • coronal plane, 2P a 1P3R a is better than 5R a 1P; on the other two sides, 5R a 1P is better than 2P a 1P3R a .

Experiment
The kinematic analysis, inverse kinematics solution and kinematic performance analysis of the two configuration design schemes show that the kinematic performance of 5R a 1P is better, so a prototype of the 5R a 1P configuration was built, as shown in Figure 22. The prototype has a generalized shoulder joint, which can realize 5-DOF movement at the shoulder joint; the elbow joint, forearm joint and wrist joint each have one active DOF. Particularly, two passive DOFs are added to the forearm man-machine contact component, which are the freedom of rotation around the forearm axis and the freedom of movement in the axis direction, as shown in Figure 23. By setting four nylon rollers, the forearm support ring is fixed on the virtual circular track, the support ring is equipped with a linear guide and a slider, and the strap is internally fixed on the slider to realize the two passive DOFs of the forearm man-machine contact parts. A spring is added between the forearm support ring and the nylon roller support frame. Two identical compression springs are placed on both sides of the sliding block in the support ring. An elastic element is added to the mechanical structure of the passive pair to make the passive pair have a certain degree of flexibility during movement, relieving the uncertainty of exercise and the instability of exercise coordination. Friction was avoided between the support ring and the patient's skin, discomfort was avoided, and the patient was restrained. Due to the safety of patient movement, the prototype was designed with mechanical limits in structure. The rotational range of the shoulder 1 movement module is 0 •~1 10 • ; the rotational range of the shoulder 2 movement module is −90 •~3 0 • ; the rotational range of the shoulder 3 movement module is −60 •~9 0 • , in which the shoulder 1, shoulder 2 and shoulder 3 movement modules correspond to the three intersecting rotation axes of the shoulder joint in sequence so that the prototype can move along the maximum space boundary to realize several kinds of movement states of human upper limbs, as shown in Figure 24. Due to the safety of patient movement, the prototype was designed with mechanical limits in structure. The rotational range of the shoulder 1 movement module is 0°~110°; the rotational range of the shoulder 2 movement module is −90°~30°; the rotational range of the shoulder 3 movement module is −60°~90°, in which the shoulder 1, shoulder 2 and shoulder 3 movement modules correspond to the three intersecting rotation axes of the shoulder joint in sequence so that the prototype can move along the maximum space boundary to realize several kinds of movement states of human upper limbs, as shown in Figure 23.

Conclusions
Based on the anatomical structure of the upper limbs of the human body, the shoulder joint ergonomic closed-chain Θ s constrained and shoulder-elbow joint ergonomic closedchain Θ se under-constrained models are established. In the shoulder joint ergonomic closed chain, the mechanical arm has five active DOFs and one passive DOF, and the 2P a 1P3R a and 5R a 1P generalized shoulder joint mechanism configuration combinations were selected with engineering application value to achieve kinematic compatibility with human upper limbs.
According to the kinematics analysis, the inverse kinematics solution, the reciprocal of the condition number, operability ellipsoid analysis and a comparison of the two configuration schemes, the results show that in the normal upper limb posture of the human body, the kinematics performances, such as the mobility flexibility and isotropic movement ability, of the 5R a 1P configuration are better than those of 2P a 1P3R a , and this configuration was selected for prototype construction, which lays a theoretical foundation for the development of upper limb rehabilitation robots in the future. Informed Consent Statement: Written informed consent has been obtained from the patient(s) to publish this paper.