Configuration and Kinematics of a 3-DOF Generalized Spherical Parallel Mechanism for Ankle Rehabilitation

The kinematic equivalent model of an existing ankle-rehabilitation robot is inconsistent with the anatomical structure of the human ankle, which influences the rehabilitation effect. Therefore, this study equates the human ankle to the UR model and proposes a novel three degrees of freedom (3-DOF) generalized spherical parallel mechanism for ankle rehabilitation. The parallel mechanism has two spherical centers corresponding to the rotation centers of tibiotalar and subtalar joints. Using screw theory, the mobility of the parallel mechanism, which meets the requirements of the human ankle, is analyzed. The inverse kinematics are presented, and singularities are identified based on the Jacobian matrix. The workspaces of the parallel mechanism are obtained through the search method and compared with the motion range of the human ankle, which shows that the parallel mechanism can meet the motion demand of ankle rehabilitation. Additionally, based on the motion-force transmissibility, the performance atlases are plotted in the parameter optimal design space, and the optimum parameter is obtained according to the demands of practical applications. The results show that the parallel mechanism can meet the motion requirements of ankle rehabilitation and has excellent kinematic performance in its rehabilitation range, which provides a theoretical basis for the prototype design and experimental verification.


Introduction
The human ankle, which is one of the three joints of the lower extremity, is of fundamental importance for balance, support, and propulsion.Nevertheless, they are particularly susceptible to musculoskeletal and neurological injuries [1].The traditional treatment of ankle injury highly depends on a one-by-one physical therapy by a doctor, which requires many human resources and heavy labor intensity.Moreover, the rehabilitation physician cannot objectively evaluate the rehabilitation status.
Evidence suggests that without sufficient rehabilitation, 44% of people will experience problems in the future, and approximately 38% of people will have recurrent activity limitations that affect their functioning [2].However, a machine-assisted treatment can collect real-time data during the entire rehabilitation process to provide targeted rehabilitation training for the patient, and it can also ensure the accuracy and repeatability of rehabilitation training.
Extensive research has been conducted on ankle-rehabilitation robots.Girone [3] developed a 6-DOF ankle-rehabilitation robot based on the Stewart platform; however, its structure was complex and difficult to control.Saglia et al. [4] proposed a 3-UPS/U (U, P, and S stand for universal, prismatic, and spherical pairs, respectively) parallel mechanism for ankle rehabilitation; however, the degrees of freedom of the mechanism was insufficient.Malosio et al. [5,6] designed an ankle-rehabilitation robot based on a 3-RRR (R stands for revolute pair) spherical mechanism.The robot had a remote rotation center, and the ankle center was located at the center of movement of the mechanism during rehabilitation.Li et al. [7][8][9] presented 3-RRS, 2-UPS/ RRR, and 3-UPS/RRR parallel mechanisms for ankle rehabilitation.These three mechanisms equated an ankle to an S-pair, and all had a remote rotation center.Fang [10] proposed a rope-driven parallel mechanism for ankle rehabilitation.The mechanism also equated the ankle to an S-pair, which solved the problem of inertial impact caused by rigid rods.Bian et al. [11,12] developed a biological fusion ankle-rehabilitation mechanism that uses the motion characteristics of the ankle to constrain the degrees of freedom of movement of the mechanism.Chen et al. [13] equated the ankle with a spatial RR model and proposed an anklerehabilitation mechanism based on the 3-UPU.
In summary, although the current ankle-rehabilitation robot meets basic rehabilitation training requirements, it still has certain limitations.Most ankle-rehabilitation robots equate the ankle and an S-pair [14] or a spatial RR model [15].The ankle is one of the most complex joints of the human body.Its structural complexity and motion particularity prevent its axes of motion from intersecting at one point; however, there are multiple instantaneous rotation axes [16,17].The equivalence of the ankle to the S-pair or spatial RR model generates a human-machine interaction force, which makes the rehabilitation training effect unsatisfactory, and even causes secondary damage to patients.To improve the fitting accuracy of humanmachine motion, this study equates the ankle with the UR model based on ankle anatomy and motion characteristics and proposes a novel 2-UPU/[RR][RRR]/PRPS generalized spherical parallel mechanism for ankle rehabilitation.
The remainder of this paper is organized as follows: In Section 2, the UR model is described with respect to the anatomical structure of the ankle.In Section 3, a parallel mechanism based on the UR model is proposed, the structure of the parallel mechanism is described in detail, and a systematic analysis of the parallel mechanism is performed, including the inverse kinematics, Jacobian matrix, singularity, and workspace.In Section 4, based on motion-force transmissibility, the parameters of the parallel mechanism are optimized.Finally, the conclusions are presented in Section 5.

Anatomical Structure and Equivalent Model of the Human Ankle
Research of the anatomical structure and motion forms of the ankle is the premise and key to the design of an ankle-rehabilitation mechanism, and it is also an important basis for evaluating the safety, man-machine coordinates, and comfort of ankle-rehabilitation mechanisms.

Anatomical Structure of the Human Ankle
As shown in Figure 1, the ankle is a highly adaptive uniaxial joint that is mainly composed of the tibiotalar and subtalar joints [18].The tibiotalar joint forms the junction between the distal tibia and fibula of the lower leg and talus, and the subtalar joint is between the talus and calcaneus.Both the tibiotalar and subtalar joints collectively bear the weight and movement of the lower extremity.The main motion forms of the ankle include dorsiflexion/plantarflexion, occurring in the sagittal plane; inversion/eversion, occurring in the transverse plane; and adduction/abduction, occurring in the frontal plane [19].The axes of motion of the ankle are not fixed and have no confluence but change constantly with the ankle movement [18].

Equivalent Model of the Human Ankle
Currently, scholars in the field of ankle rehabilitation consider the spatial RR model or S-pair as the equivalent model of the ankle.The two rotation axes of the spatial RR model correspond to the rotation axes of the tibiotalar and subtalar joints of the ankle; however, they ignore the motion forms of the ankle.The S-pair considers the motion forms of the ankle and ignores the real structure of the ankle.Liu et al. [20] proposed a series of equivalent ankle models with respect to the physiological structure and motion characteristics of the ankle.According to the anatomical structure of the human ankle, the direction and position of the rotation axis of the tibiotalar joint change with its activities, whereas the rotation axis of the subtalar joint is relatively stable [21].Therefore, this study adopts the UR model as the kinematic equivalent model of the human ankle, in which the tibiotalar joint is equivalent to the U-pair and the subtalar joint is equivalent to the R-pair.The distance from the geometric center of the U-pair to the R-pair is the size of the human talus, as shown in Figure 2. Owing to the limitation of the size of the human talus, it is difficult for the serial equivalent model to meet rehabilitation requirements.Therefore, a parallel mechanism was designed with respect to the motion forms of the serial equivalent model to improve the rehabilitation accuracy of the ankle.With point O as the origin of the coordinate system, the X-axis coincides with OA 3 , the Z-axis is perpendicular to the fixed platform upward, and the Y-axis satisfies the right-hand rule; thus, the fixed coordinate system O-XYZ is established.Point O 1 is taken as the origin of the coordinates, the x-axis is parallel to the direction of B 2 B 1 , the z-axis is perpendicular to the platform moving upward, and the y-axis satisfies the right-hand rule; thus, the moving coordinate system O 1 -xyz is established.

Degrees of Freedom
The DOF of the 2-UPU/[RR][RRR]/PRPS parallel mechanism for ankle rehabilitation was analyzed using screw theory [22].In screw theory, a unit screw $ is defined as where s is the unit vector along the direction of the screw axis, r is the position vector of any point on the screw axis, and h is the pitch.
The motion screws of Limb 1 are given in the fixed coordinate system by where s OA1 is a unit vector along the direction of the first rotation axis of the U-pair in Limb 1 connected to the fixed platform; s 12 is a unit vector along the direction of the second rotation axis of the U-pair in Limb 1 connected to the fixed platform; r A1 is the position vector of point A 1 ; s A1B1 is a unit vector along the direction of the P-pair in Limb 1; r B1 is the position vector of point B 1 ; s O1B1 is a unit vector along the direction of the first rotation axis of the U-pair in Limb 1 connected to the moving platform; and r O1 is the position vector of point O 1 .
Then, by employing the reciprocal screw theory, the constraint screw of Limb 1 $ C1 can be obtained as where G 1 is the intersection of the rotation axes in Limb 1 connected to the fixed and moving platforms, and r G1 is the position vector of point G 1 .
The motion screws of Limb 2 can be expressed as where s OA2 is a unit vector along the direction of the first rotation axis of the U-pair in Limb 2 connected to the fixed platform; s 22 is a unit vector along the direction of the second rotation axis of the U-pair in Limb 2 connected to the fixed platform; r A2 is the position vector of point A 2 ; s A2B2 is a unit vector along the direction of the P-pair in Limb 2; r B2 is the position vector of point B 2 ; and s O2B2 is a unit vector along the direction of the first rotation axis of the U-pair in Limb 2 connected to the moving platform.
In a similar approach, the constraint screw of Limb 2 $ C2 can be obtained as where G 2 is the intersection of the rotation axes in Limb 2 connected to the fixed and moving platforms, and r G2 is the position vector of point G 2 .
The motion screws of Limb 3 can be written as where s OA3 is a unit vector along the direction of the rotation axis of the R-pair in Limb 3 connected to the fixed platform; s OB3 is a unit vector along the direction of the rotation axis of the second R-pair in Limb 3; s O1C3 is a unit vector along the direction of the rotation axis of the third R-pair in Limb 3; s O1D3 is a unit vector along the direction of the rotation axis of the fourth R-pair in Limb 3; and s O1E3 is a unit vector along the direction of the rotation axis of the fifth R-pair in Limb 3.
The constraint screw of Limb 3 $ C3 can be obtained as where s OO1 is a unit vector along the direction of the 'double-centered line.Limb 4 is a PRPS limb with 6-DOF and no constraint screw.
The motion and constraint screws of the UPU limb and the [RR][RRR] limb are shown in Figure 5.According to Eqs. (3), (5), and (7), the parallel mechanism has three linearly independent constraint forces that are arbitrarily distributed in space.Under the restriction of three constraint forces, the mechanism can realize 2-DOF motions of the moving spherical center O 1 on the spherical surface with the fixed spherical center O as the center and OO 1 as the radius, which can rotate around the moving spherical center O 1 with 1-DOF, which is consistent with the motion characteristics of the UR equivalent model.The DOF of the parallel mechanism can be obtained using the modified G-K formula [23]: where F denotes the number of DOF, d denotes the rank of the parallel mechanism, n denotes the number of ( 5) components, g denotes the number of joints, f i denotes the number of DOF for the ith joint, v denotes the number of redundant constraints, and ζ denotes the number of isolated DOF.

Inverse Kinematic Analysis
Let the distance from point O to A 1 be l OA1 .Then, distance from point O to A 2 is denoted by l OA2 ; the distance from point O 1 to B 1 is denoted by l O1B1 ; the distance from point O 1 to B 2 is denoted by l O1B2 ; the length of the 'double-centered line' is denoted by l; the distance from point O 1 to E 3 is denoted by l O1E3 ; and the distance from point E 3 to H is denoted by l EH .The angle between the axes of the U-pair connected to the moving platform and the moving platform is φ.l OA1 =l OA2 , l O1B1 =l O1B2 The RPY angles α, β, and γ are used to represent the orientation of the moving platform relative to the fixed platform.The initial attitude of the moving coordinate system O 1 -xyz is like that of the fixed coordinate system O-XYZ.First, O 1 -xyz is rotated by α about the X-axis, then by β about the Y-axis, and finally by γ about the Z-axis.Thus, the rotation transformation matrix can be expressed as follows: where C represents a cosine function and S represents a sine function.
Let the normal vectors of planes OA 1 B 1 O 1 and OA 2 B 2 O 1 be n 1 and n 2 , respectively.Then, The intersection line of planes OA 1 B 1 O 1 and OA 2 B 2 O 1 is a straight line where the 'double-centered line' OO 1 is located with respect to the structural characteristics of the parallel mechanism, as shown in Figure 6.Thus, the 'double-centered line' OO 1 can be obtained as The position vector of the origin (point O 1 ) of the moving coordinate system in the fixed coordinate system can be obtained as follows: The position vectors of points B 1 , B 2 , and H in the moving coordinate system can be written as Through the transformation of the RPY angles, the position vectors of points B 1 , B 2 , and H in the fixed coordinate system can be expressed as: Then, the lengths of Limbs 1 and 2, l 1 and l 2 , are obtained as The angle between vector MH and the guide rail can be expressed as Then, the slider movement distance l 3 can be obtained as

Jacobian Matrix
The Jacobian matrix of a parallel mechanism, which represents the mapping between the joint input rates and moving platform output velocity, is an important tool for analyzing the kinematic performance and singularity of the parallel mechanism.Screw theory is used to establish the complete Jacobian matrix of the parallel mechanism.
The instantaneous screw of the moving platform [24] can be expressed as $ p = [w T v T ] T , which should be a linear combination of screws in Limb i (i=1-4).
where θij is the rotational angular velocity of the jth (j=1- 6) joint in the ith limb; ḋij is the linear velocity of the pris- matic in the ith limb; and $ ij is the unit screw of the jth joint in the ith limb.
As mentioned in Section 3.2, the parallel mechanism has three constraint screws.The constraint Jacobian matrix is obtained based on $ Ci ° $ p = 0. ( 14) If we lock the actuated prismatic joint of the ith limb, a new unit-transmission screw $ Ti can be obtained, which represents the intermediate medium that transmits the power from the input to the output terminal [25].
Taking the reciprocity product on both sides of Eqs. ( 18) and (20) with $ Ti yields: According to Eqs. ( 21) and ( 23), the velocity model of the 2-UPU/[RR][RRR]/PRPS parallel mechanism can be expressed as follows: where J p denotes the forward Jacobian matrix, J q denotes the inverse Jacobian matrix, q denotes the velocity of the input joint.
When the parallel mechanism is away from singularities, we have

Singularity Analysis
When the parallel mechanism approaches the singularity configuration, the DOF is reduced or increased, which makes the parallel mechanism uncontrollable in certain directions.Based on the determinant of the Jacobian matrix, the singularity configuration of the parallel (21) mechanism is divided into three types [26,27].The first type, inverse singularity, occurs when det(J q ) = 0; the second type, forward kinematic singularity, occurs when det(J p ) = 0; and the third type, called combined singularity, occurs when det(J q ) = 0 and det(J p ) = 0.Because the human subtalar joint is always below the tibiotalar joint, this study considers the singularity of the workspace when the moving spherical center is lower than the fixed spherical center.

Inverse Singularity
Because matrix J q is a diagonal matrix, inverse singularity occurs whenever any of the diagonal elements becomes zero.Moreover, all actuators of the parallel mechanism are prismatic pairs, and the direction of the real unit of $ Ti is consistent with the direction of the ith actuator.

Forward Kinematic Singularity
Matrix J p consists of three zero-pitch constraint screws and three zero-pitch transmission force screws.A zeropitch screw can be represented as a line; therefore, the linear dependency among the screws becomes equivalent to the dependency between the lines they represent.Therefore, the forward kinematic singularity can be identified using the Grassmann line theory [28].
The three constraint forces provided by the four limbs to the moving platform are not always in the same plane as the change in the parallel mechanism configuration or position of the moving platform.According to the Grassmann line theory, the rank of the matrix composed of three constraint screws is always equal to three.Therefore, the condition for the forward kinematic singularity is that the rank of matrix J p should be 3-5.
(1) The rank of matrix J p is three, and there are three possible cases.
Case 1: The union of two flat pencils with a line in common but lying in distinct planes and with distinct centers.Because the three constraint screws are not always in the same plane, $ C1 is perpendicular to $ T1 , $ C2 is perpendicular to $ T2 , and the 'double-centered line' OO 1 with A 1 B 1 and A 2 B 2 can never constitute a plane.Therefore, this case does not occur in the parallel mechanism.
Case 2: All lines pass through a point but are not coplanar.This case does not exist because the three transmission-force screws are consistent with the direction of the actuators, and the three constraint screws do not meet at one point.Case 3: All lines are in a plane but do not constitute a planar pencil of lines.Because the three constraint screws are not always in the same plane, this case does not exist.
(2) Matrix J p is ranked forth, and there are three possible cases.
Case 1: All lines in a plane pass through one point on that plane.According to the structure of the parallel mechanism, only three transmission-force screws may be in the same plane, but the three constraint screws never meet at one point, so this case does not exist.
Case 2: A one-parameter family of flat pencils, with one line in common and forming a variety.According to the structure of the parallel mechanism, $ C1 is perpendicular to $ T1 and is a non-coplanar straight line with $ C3 , $ T2 , and $ T3 .In the same way, $ C2 is perpendicular to $ T2 , and is a non-coplanar straight line with $ C3 , $ T1 and $ T3 .Moreover, $ C1 and $ C2 can only form a plane at the initial position.At this position, $ T1 and $ T2 can form a plane, but $ T3 and $ C3 cannot form a plane.Therefore, this case does not exist.
Case 3: All lines are concurrent with two skewed lines.According to the analysis of Case 2, this case does not exist.
(3) The rank of the matrix J p is five, and there is one possible case: all the lines intersect one line.Because any one of these six lines has at least one set of noncoplanar lines, this case does not exist.
Through the above analysis, the parallel mechanism has no forward kinematic singularity.

Combined Singularity
A combined singularity occurs when forward and inverse kinematic singularities occur simultaneously.Because there is no inverse singularity or forward kinematic singularity in the parallel mechanism, it has no combined singularity.

Workspace
The workspace is a necessary condition for measuring whether the parallel mechanism can meet the rehabilitation requirements of the human ankle.In this study, the orientational workspace, moving spherical center, and moving platform workspace of the parallel mechanism for ankle rehabilitation were solved.The structural parameters of the parallel mechanism are presented in Table 1.
Under the condition that the moving distances of the P-pairs in Limbs 1 and 2 are limited to ±80 mm and the rotation range of the U-pair is ±45º, the workspaces of the parallel mechanism are solved using the search method with the help of MATLAB.The results are shown in Figure 7.According to Ref. [29], the motion ranges of the human ankle and the parallel mechanism are listed in Table 2.
From the data in Table 2, the workspace of the parallel mechanism fully meets the motion range requirements of the human ankle.Therefore, a parallel mechanism for ankle rehabilitation can complete the rehabilitation exercises of the ankle.

Performance Analysis and Optimization
Currently, the performance evaluation indices applied to parallel mechanisms mainly include the condition number of the Jacobian matrix [30], dexterity [31], and motion/force transmissibility [32].Based on screw theory, the motion-force transmissibility reflects the transmission efficiency of the motion and force of the parallel mechanism from the input to output.In this study, the motion/force transmissibility was used as the evaluation standard of the kinematic performance of the ankle-rehabilitation parallel mechanism, and the parameters were optimized based on this index.

Motion/Force Transmissibility
Motion-force transmissibility can be divided into the input and output transmissibility with respect to the transfer of objects between the input and output of the parallel mechanism.The input and output transmissibility of a single limb of the parallel mechanism can be expressed as: where λ i denotes the input transmissibility, η i denotes the output transmissibility, $ Ii denotes the input screw of the ith driving limb, $ Ti denotes the transmission-force screw (26) of the ith driving limb, and $ Oi denotes the output screw of the ith actuated limb.Both λ i and η i were within the range of [0, 1].
To ensure that each limb has excellent input and output transmissibility, the local transmission index (LTI) [33] of the parallel mechanism is defined as where Γ is a dimensionless index, independent of the coordinate system.The closer Γ is to one, the better the motion-force transmissibility of the parallel mechanism.The closer Γ is to zero, the closer the parallel mechanism is to singularity.
Furthermore, if any two actuated prismatic joints are locked and only the actuated prismatic joint of the ith limb is retained, then only the motion from the ith actuated prismatic joint can be transmitted to the moving platform under the action of the ith transmission force.In this situation, the parallel mechanism becomes a 1-DOF mechanism and the unit instantaneous screw of the moving platform can be expressed by the unit output screw $ Oi .Therefore, according to the reciprocity of the motion and force screws, the unit output screw $ Oi can be obtained as follows: The LTI of the parallel mechanism can be obtained by substituting the output motion screws into Eqs.( 26)- (28).The distribution diagrams of the LTI of the parallel mechanisms in pure dorsiflexion/plantarflexion (α=0º), inversion/eversion (β=0º), and adduction/abduction (γ=0º) are illustrated in Figure 8(a)-(c), respectively.It can be observed from the figures that the LTI of the parallel mechanism in the central area of the workspace is greater than 0.7, and its value decreases progressively as it leans on the boundary of the workspace.There is no Γ=0 position in the target workspace, indicating that the parallel mechanism does not have singularity in the workspace.

Parameter Optimization
However, the LTI can only reflect the motion-force transmissibility of the parallel mechanism in a particular configuration but not in the entire workspace.Therefore, the global motion-force transmissibility, ζ, was introduced.When Г ≥ 0.7, the configuration sets of the parallel mechanism are defined as the high-quality transmission workspace [34].The ratio of the high-quality transmission workspace to the entire workspace is the global motion-force transmissibility:   closer ζ is to one, the better is the kinematic performance of the parallel mechanism.By substituting the structural parameters in the previous section into Eq.( 30), the global motion-force transmissibility of the parallel mechanism was 0.33.This shows that the motion-force transmissibility of the parallel mechanism is not excellent under such structural parameters.Therefore, further optimization is required.
In this study, the optimization method proposed in Ref. [35] is used to optimize three structural parameters: the distance from point O to point A i , l OA , distance from point O 1 to B i , ; and distance from point H to E 3 , l EH .First, the structural parameters are dimensionless and are treated as follows: where D is the normalized factor and r i is the dimensionless parameter of the three optimization parameters.
Considering the interference between the various components of the parallel mechanism and stroke of the push link, the dimensionless parameters should satisfy the following conditions: According to Eq. ( 32), the parameter-optimization region of the parallel mechanism can be obtained, as shown in Figure 9(a).It is projected onto the s-t coordinates to obtain Figure 9(b); therefore, the three-dimensional space is reduced to a two-dimensional space to reduce the optimization parameters.The mapping relationship is obtained as follows: (31) MATLAB was used to conduct an iterative search over the entire s-t region, as shown in Figure 9(b).The global motion-force transmissibility of the corresponding size was calculated, and performance atlases were plotted, as shown in Figure 10.The figure shows that the global motion-force transmissibility of the parallel mechanism in the range of 0.6 ≤ s ≤ 1.22 and 0 ≤ t ≤0.2 is greater than 0.5.This shows that the parameters (33) Table 2 Motion range of the human ankle and the parallel mechanism Considering the rationality of the layout between the links of the parallel mechanism and workspace, s=0.8, t=0.2, and D=200 mm were selected.The corresponding mechanism parameters are l OA =D•r 1 ≈255 mm, l O1B =D•r 2 ≈215 mm, and l EH =D•r 3 ≈131 mm.Then, the parameters are substituted into Eqs.( 26)-( 28), and the distribution diagrams of the LTI are shown in Figure 11.
By comparing the distribution diagrams of the LTI of the parallel mechanism before optimization, the area of the high-quality transfer workspace after optimization was significantly increased, and the global motion-force transmissibility of the parallel mechanism after optimization reached 0.59.This indicates that the kinematic performance of the optimized parallel mechanism was significantly improved.

Conclusions
(1) Based on the anatomical structure and motion characteristics of the human ankle, the UR equivalent model was selected, and a novel 2-UPU/[RR] [RRR]/PRPS generalized spherical parallel mechanism for ankle rehabilitation was presented.Screw theory analysis revealed that the parallel mechanism has three rotational DOFs that satisfy the demand for ankle rehabilitation.(2) The inverse kinematics of the parallel mechanism were analyzed using an analytical method.According to screw theory, a complete Jacobian matrix was established, and the singularity was investigated based on inverse singularity, forward kinematic singularity, and combined singularity, indicating that the parallel mechanism has no singularity.
Additionally, the workspaces were solved, and they showed that the parallel mechanism satisfied the motion range of the human ankle.(3) By considering the motion-force transmissibility, the global motion-force transmissibility was used as the performance evaluation criterion in this study, and the performance atlases were plotted in the parameter-optimal design space.Subsequently, according to the demands of practical applications, the optimum region was obtained.The results show that the ratio of the high-quality transmission workspace reached 0.59, which indicates that the parallel mechanism has excellent kinematic performance in the ankle-rehabilitation motion range.

Figure 1
Figure 1 Anatomical structure of the ankle

Figure 3
Figure 3 shows the structure of the 2-UPU/[RR][RRR]/ PRPS parallel mechanism for ankle rehabilitation based on the UR-equivalent model.The parallel mechanism has two rotation centers, namely a fixed spherical center O, which is equivalent to the rotation center of the tibiotalar joint, and a moving spherical center O 1 , which is equivalent to the rotation center of the subtalar joint.Point O is located at the intersection of the rotation axes on the fixed platform and point O 1 is located at the intersection point of the rotation axes on the moving platform.The line between points O and O 1 is called the 'double-centered line' , and its length is equal to the size of the talus.The parallel mechanism has four limbs between the fixed and moving platforms.Limbs 1 and 2 are UPU limbs, and the rotation axes connected to the fixed platform are perpendicular to each other and intersect at point O.The rotation axes connected to the moving platform intersect at point O 1 .The rotation axes connected to the P-pair are parallel to each other.Limb 3 is the [RR][RRR] limb, in which the first two rotation axes intersect at point O and the remaining rotation axes intersect at point O 1 .Limb 4 is an unconstrained PRPS limb in which the axis of the R-pair coincides with the direction of the guide rail.The P-pairs of Limbs 1 and 2 and the slider of Limb 4 are the actuators.As shown in Figure 4, point A i (i=1, 2) is the center of the U-pair connected to the fixed platform in the ith limb.Point B i (i=1, 2) is the center of the U-pair connected to the moving platform in the ith limb.Points A 3 and A 3 ′ are the centers of the R-pairs connected to the fixed platform in Limb 3. Point E 3 is the center of the R-pair connected to the moving platform in Limb 3. Point H is the center of the S-pair in Limb 4. Point M is the center of the guide rail.With point O as the origin of the coordinate system, the X-axis coincides with OA 3 , the Z-axis is perpendicular to the fixed platform upward, and the Y-axis satisfies the right-hand rule; thus, the fixed coordinate system O-XYZ is established.Point O 1 is taken as the origin of the coordinates, the x-axis is parallel to the direction of B 2 B 1 , the z-axis is perpendicular to the platform moving upward, and the y-axis satisfies the right-hand rule; thus, the moving coordinate system O 1 -xyz is established.

Figure 2 Figure 3 Figure 4
Figure 2 Equivalent model of the human ankle

Figure 5 Figure 6
Figure 5 Motion and constraint screws ) ζ = SG dW S dW , where W is the entire workspace, SG is the volume of the high-quality transfer workspace, S is the volume of the entire workspace, and ζ is within the range of [0, 1].The

Figure 10 Figure 11
Figure 10 Distribution of LTI in the parameters design

Table 1
The structural parameters of the parallel mechanism