Journal of Healthcare Engineering
Volume 2018 (2018), Article ID 8647591, 11 pages
https://doi.org/10.1155/2018/8647591

Research Article

Dynamic Modeling and Interactive Performance of PARM: A Parallel Upper-Limb Rehabilitation Robot Using Impedance Control for Patients after Stroke

Hui Guang,¹ Linhong Ji,¹ Yingying Shi,² and Berno J. E. Misgeld³

¹Department of Mechanical Engineering, Tsinghua University, Beijing, China
²Department of Mechanical Engineering, Beihang University, Beijing, China
³Helmholtz-Institute for Biomedical Engineering, RWTH Aachen University, Aachen, Germany

Correspondence should be addressed to Linhong Ji; jilh@tsinghua.edu.cn

Received 1 June 2017; Revised 15 September 2017; Accepted 13 February 2018; Published 5 April 2018

Academic Editor: Carlo Ferraresi

Copyright © 2018 Hui Guang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

The robot-assisted therapy has been demonstrated to be effective in the improvements of limb function and even activities of daily living for patients after stroke. This paper presents an interactive upper-limb rehabilitation robot with a parallel mechanism and an isometric screen embedded in the platform to display trajectories. In the dynamic modeling for impedance control, the effects of friction and inertia are reduced by introducing the principle of virtual work and derivative of Jacobian matrix. To achieve the assist-as-needed impedance control for arbitrary trajectories, the strategy based on orthogonal deviations is proposed. Simulations and experiments were performed to validate the dynamic modeling and impedance control. Besides, to investigate the influence of the impedance in practice, a subject participated in experiments and performed two types of movements with the robot, that is, rectilinear and circular movements, under four conditions, that is, with/without resistance or impedance, respectively. The results showed that the impedance and resistance affected both mean absolute error and standard deviation of movements and also demonstrated the significant differences between movements with/without impedance and resistance (). Furthermore, the error patterns were discussed, which suggested that the impedance environment was capable of alleviating movement deviations by compensating the synergetic inadequacy between the shoulder and elbow joints.

1. Introduction

Stroke is caused by cerebrovascular accident and is one of the leading diseases of disability, motor disorder, and deterioration of activities of daily living (ADL). The incidences in the European Union and the United States are approximately one million and 0.8 million per year, respectively [1, 2], and thirty percent of patients suffer recurrent attacks, which results in increasing demand for rehabilitation services.

For patients after stroke, the task-repetitive training has been demonstrated to be effective in improving their upper and lower extremity functions and ADL [3]. To meet the requirement for repetitive training, various upper-limb rehabilitation robots have been developed over the past twenty years, which are generally classified into two categories [4]: end-effector robots, such as DIAGNOBOT [5], CARR [6], MIT-MANUS [7], MIME [8], GENTLE/s [9], and exoskeleton robots, such as CADEN-7 [10], RUPERT [11], BONES [12], and ARMin [13]. Since the robotic rehabilitation exhibits the advantages in terms of high-dosage, high-intensity, and task-specific training [14], randomized controlled trials comparing the robot-assisted and conventional therapy have yielded significant effects of robots on the improvements of limb function [15, 16] and even ADL [17].

Although many robots for the upper-limb rehabilitation have been developed, mechanical design, control, and training methods remain an area of interest. As pointed out by Belda-Lois et al. [18], robot-assisted rehabilitation could be enhanced by means of precisely controllable assistance or resistance, enhanced training motivation through interactive feedback, and quantifiable and objective measures of subject performance. Besides, cost should also be considered [19].

Generally, the exoskeleton robots take individual joint motions into account to minimize abnormal postures and joint motions. Nevertheless, due to the complexity of the human upper-limb anatomy, the instantaneous centers of rotation of the upper-limb joints are changed with movement [20], which causes the inconvenience of joint axis alignments and raises interactive force between human and robots [21], thereby obstructing the development and application of the exoskeleton robots. In contrast to exoskeleton robots, end-effector robots are simple and cost-effective and can adapt to patients with diverse somatotypes [4]. Despite the disadvantage of end-effector robots in joint training, extensive research has also demonstrated their effectiveness and superiority for improving upper-limb function and ADL in comparison with conventional therapies [14].

Compared to serial mechanisms, parallel mechanisms exhibit inherent advantages of low inertia, high stiffness, and satisfactory payload capability [22, 23]. More importantly, as the end-effector is controlled in parallel, the errors of the joint control are not accumulated and amplified by serial counterparts, and thus the manipulator is less affected by joint clearance and has higher precision in aspects of position, stiffness, and interactive force control [12]. Therefore, parallel manipulators have been recently applied to rehabilitation robots, including shoulder [12, 24], wrist [25], hip [26], and upper-limb rehabilitation devices [27].

Another issue is that understanding sensorimotor physiology is more imperative prior to developing a rehabilitation robot. For instance, one aspect is how individual joints, as well as segments, are coordinated to achieve the task. In physiology, limb movements are perceived in an egocentric reference frame, in which targets are defined with respect to the trunk or head. In contrast, an allocentric reference frame represents the coordinate system external to the body [28]. However, for current training robots, target and actual trajectories are presented in a standing monitor, which is a virtual environment based on the allocentric reference frame for patients. Thus, patients are required to transform the targets and movements in the virtual environment to the egocentric reference frame to accomplish the task, causing difficulties in perception and sensorimotor control. Besides, it might weaken the effect of proprioceptive training since the actual positions do not directly correspond to virtual positions.

Based on the issues discussed above, a novel end-effector-based upper-limb rehabilitation robot, which is named PARM, is developed with a parallel mechanism and patient-frame-based interactive feedback to enhance training performance. Distinct from other rehabilitation robots, a monitor was embedded in the platform to show target and actual trajectories, providing isometric direct visual feedback for patients. The trajectories displayed on the platform screen were the same as the actual trajectories in movement space, particularly in the aspects of scale, position, and direction. Therefore, patients could perceive targets and movements in the egocentric reference frame, which should improve the motor recovery and proprioceptive training. As the precise control of position, stiffness, and force contributes to training effects [18], the impacts of friction and acceleration were incorporated to improve the control precision. Consistent with the robots such as MIT-MANUS [7], an assist-as-needed strategy was also introduced in PARM to improve interaction between patients and robots. In the assist-as-needed control, patients determine the manipulator in terms of position, velocity, and acceleration; thus, the reference positions are variable with the movement and associated with real-time deviations.

To summarize, the novelty of our work is the strategy based on the orthogonal deviation for assist-as-needed impedance control, which aims to obtain the equilibrium positions and calculate impedance force, and the hardware which adopts parallel mechanism and isometric visual feedback. Simulated and experimental results validated the dynamic modeling and impedance control. Since the mechanism of the impedance control contributing to the motor coordination is still less clear, the functional interaction between impedance control and movements was also discussed.

2. Apparatus and Specification

The rehabilitation robot PARM aims to improve the motor performance of stroke patients by enhancing movement interaction between the patients and the robot. This interactive robot incorporates multiple training modes for patients with diverse disability and recovery stages, which are summarized as patient-passive training and patient-active training (Figure 1). Arbitrary reference trajectories are predefined by therapists prior to training. In the patient-passive training, the movements are entirely actuated by the robot with position control, in which the robot is a mechanical admittance whereas the patient’s arm is regarded as an impedance. Contrastively, in the patient-active training, movements are initiated and actuated by patients with partially assistance or resistance. For instance, in the training with impedance and propulsion, the impedance force towards the target trajectory aims to rectify deviations, while the propulsive force towards the movement direction could reduce the active force of the patient, which decreases the task difficulty. Conversely, the resistance force against the movement direction increases the movement effort.

Figure 1: Training modes of PARM for stroke patients.

To increase the benefit of robot-aided therapy, control schemes should be customized for individuals and adopted to patients’ deficits in upper-limb motor function, based on their poststroke stages and clinical assessments. For subacute patients, since they are generally unable to perform voluntary arm-reaching tasks due to dystonia, training is mainly executed in patient-passive modality. For chronic patients (more than 6 months poststroke), robot-aided therapies are generally performed in patient-active modalities to enhance patient engagement. Specifically, when patients could perform inaccurate arm-reaching tasks, impedance control is applied for this training stage to rectify deviations and improve their abilities in motor control. Besides, for patients with low strength, propulsion in movement direction is included to reduce the movement effort. However, for the patients having coordinated motor control, impedance control is removed. Instead, resistance in movement direction might be involved to match their motor function and improve training outcomes.

PARM mainly consisted of a lifting platform, two monitors, a five-bar parallel mechanism with two motors and actuators, and a three-axis force sensor (Figure 2). A horizontal monitor was embedded in the platform to display the reference trajectory and actual trajectories, providing direct visual feedback for patients, while a standing monitor was used to display the configuration of training parameters and quantitative assessments. Additionally, the inclination and height of the platform could be adjusted to make the trajectories conveniently observed for patients.

(a)

(b)

(c)

Figure 2: Mechanical construction of PARM. (a) Robot overview. (b) Robot parallel mechanism. (c) Motor and joint connections.

The five-bar parallel mechanism RPRPR (revolute-prismatic-revolute-prismatic-revolute) is shown in Figures 2(b) and 2(c). Linear rails were adopted to increase movement range of end-effector and improve kinematic precision, and linear bearing blocks constituted the prismatic joints to reduce friction. Each side consisted of three prismatic joints and two linear rails, and thereby the minimal length of each side was the length of a rail, while the maximal length was the sum of two rails and a link. During movement, the length of two sides was accordingly adjusted to the two revolute joints controlled by the motors.

The parallel mechanism was actuated by two Maxon RE50 DC motors with shaft keys, and connectors concatenated the linear rails and motors by screws and shaft keys, respectively. The motors were fixed on the platform and in serial with angle encoders, and the nominal voltage, maximum torque, and torque constant of the motor were 36 V, 418 mNm, and 60.4 mNm/A, respectively. The Maxon gearboxes EP52C, whose gear ratios were 43 : 1, modulated the motor outputs. The motors were actuated by Maxon EPOS2 70/10, and control programs were coded in LabVIEW (NI, USA).

In the patient-active training, an assisted-as-needed strategy was introduced by employing impedance control. To improve the control precision of manipulator impedance, a three-dimensional force sensor was mounted on the end-effector, which could additionally record the interactive force between patients and the robot for quantitative assessments. The sensor signals were collected by Arduino board (Mega 2560) and subsequently transmitted to PC through USB serial communication. The end-effector components are shown in Figure 3. Since the angle between the upper-limb and the links changed with movements when patients hold the handle, the handle had one independently rotational degree of freedom (DoF) with respect to the links, and the revolute friction was reduced by thrust bearings. The cone below the end-effector was used to indicate movement positions, and Teflon was adopted to reduce the friction with screen. In addition, the revolute joint of two links was constituted by an axis, and friction was also reduced by thrust bearings.

Figure 3: End-effector components.

3. Kinematic and Dynamic Modeling

3.1. Kinematics of PARM

PARM has two DoF actuated by two servo motors, and the kinematic diagram is shown in Figure 4. The end-effector position was determined by the joints , which is given by

where means the distance between two joints.

Figure 4: Kinematic diagram of PARM.

Patient-passive control is based on the inverse kinematics of the robot arm. For continuous predefined trajectories , controlled joints is calculated as

The calculated joint angles are implemented with position control of servo motors.

Deriving (1), Jacobian matrix denotes the relationship between the end-effector and joint velocity, which is given by

Besides, link lengths and were adjusted automatically to the joint angles and as

Similarly, the relation between the elongation velocity of the two links and the joint angular velocity is expressed by Jacobian matrix as

3.2. Dynamic Modeling of PARM

For dexterous and accurate control of a manipulator, inertia and friction should be considered. In this study, it is hypothesized that three components constituted the motor torques. Namely, the first component counteracted friction; another component compensated the inertia of the end-effector, links, and motor-gear system; the last component generated the manipulator impedance and achieved the flexibility. The dynamic diagrams of the mechanism and motor-gear system are shown in Figures 5 and 6, respectively, where the arrows indicate the positive references.

Figure 5: Robotic dynamics of the mechanism. In this diagram, and denote the motor torques; and are the friction in the prismatic joint; and denote the friction between the end-effector and platform; and denote the interactive force between patients and the robot; , , and denote the inertia of the link, end-effector, and the moment of inertia of the motor shaft, respectively; and represent the acceleration of the joints and end-effector, respectively.

Figure 6: Dynamics of the motor-gear system. and denote the moments of inertia of the motor and gear shafts, respectively; , and , are the friction torques and angular accelerations of the two shafts, respectively; and represent the motor and load torques, respectively; is the gear ratio.

3.2.1. Friction Component

In this study, the principle of virtual work was utilized to deduce the equilibrium relations. Specifically, and were virtual displacement of the motor joints and end-effector, respectively, and was the corresponding virtual change of link length. In the patient-active training, the end-effector was mechanical impedance, and conversely, the human arm was regarded as mechanical admittance [29]. Thus, the equation can be written as

where denotes the joint torques counteracting friction and external force; is the external force acting on the end-effector; represents the friction in the prismatic joint; means the friction between the end-effector and platform, and . Based on (3) and (5), (6) can be deduced as

To calculate the component of the joint torques which only counteracted the friction, the external force should be excluded. Eliminating the term , (7) can be written as

where the joint-parameter matrixes and .

According to the motor-gear system shown in Figure 6, motor torque counteracting the friction is derived as

where and denote the friction torques of motor and gear shafts, respectively, and is the gear ratio.

3.2.2. Inertia Component

The joint torque was assumed to compensate the inertia force generated by joint angular acceleration and end-effector acceleration , which can be calculated as

where denotes the mass of the end-effector, and and are the moments of inertia of left and right links, respectively.

The end-effector acceleration could be obtained by the derivative of (3), which is given by

In summary, the motor torque counteracting the inertia force could be calculated as

where and are the moments of inertia of the motor and gear shafts, respectively.

4. Impedance Control for Assist-as-Needed Training

In impedance control, the end-effector behaves as a damped spring-mass system, which is represented in a single DoF system as

where denotes the external force; parameters , , and are the dynamic parameters of the end-effector corresponding to mass, damping, and spring, respectively; represents the desired equilibrium position, while denotes the actual end-effector position.

In robot-aided training, the predefined target trajectory meant the movement that patients were expected to track, which was, however, supposed to be different with the actual trajectory due to movement error. Actual trajectories were obtained by joint sensors and forward kinematics. Assume denoted the desired position on the predefined target trajectory, when was the actual end-effector position. Since the reference position determined the direction and magnitude of the impedance force, it was significant to search the appropriate reference position. In assist-as-needed training, patients determined the manipulator in terms of position, velocity, and acceleration; thus, the reference positions were variable with movements and associated with the real-time deviations. In this study, the strategy based on the orthogonal deviation was proposed to define the desired equilibrium positions for arbitrary predefined trajectories. As illustrated in Figure 7, the curve represents the predefined target trajectory, and force and are the interactive force detected by the three-axis force sensor in and directions, respectively. In this strategy, the direction of the actual position relative to the reference point was orthogonal to the tangent of the predefined trajectory at the equilibrium position, which indicated that the reference position was the point on the predefined trajectory closest to the current position . Therefore, the tangent component of the external force with respect to the equilibrium point provided the propulsive force along the trajectory, while the normal component force was supposed to be the impedance force () shown in (13), and denotes the deviation input to the impedance control.

Figure 7: Diagram of the strategy based on the orthogonal deviation.

The dynamic modeling was used to calculate the motor torques to generate the required impedance according to deviations and impedance parameters. However, even though the dynamic modeling incorporated friction and inertia, control errors inevitably occurred in experiments; thus, the force sensor was utilized to obtain actual interactive force as feedback to reduce the errors. Since the acceleration and friction components have been discussed in Section 3, let be the motor torque-generating manipulator impedance, which is given by

where denotes the error feedback coefficient.

As the tangent force illustrated in Figure 7, for propulsion/resistance control, motor torque is implemented as

where denotes corresponding assistive/resistive tangent force along the predefined trajectory.

Summarizing (9), (12), (14), and (15), as the control scheme shown in Figure 8, the motor torque for impedance control is calculated as

(a)

(b)

Figure 8: Control schemes of impedance and propulsion/resistance control. (a) Overview. (b) Controller diagram. , , , and are calculated by the strategy based on orthogonal deviation (Figure 7), and denotes propulsion/resistance force.

5. Simulations and Experiments

5.1. Impedance Parameter Determination

Impedance parameters , , and , which determined the dynamic behavior of the manipulator, were optimized by simulations. Specifically, it was supposed that the end-effector was released from the initial coordinates m, while the equilibrium position was m. was set to 0.8 kg, which was the approximately actual mass of the end-effector, whereas the damping and stiffness coefficients and ranged from 10 to 50 Ns/m and 200 to 600 N/m, respectively. The dynamic responses of the manipulator in the absence of external interaction are shown in Figure 9. The result showed that the oscillation deteriorated with larger and smaller and the response time and overshoot were the least when  N/m and  Ns/m. Therefore, and were set to 200 N/m and 30 Ns/m for experiments, respectively.

Figure 9: Manipulator responses with different impedance parameters.

5.2. Comparison between Experimental and Desired Responses

To validate the dynamic modeling and impedance control, an experiment of the deviation-regression response, which was the same as the simulation introduced in Section 5.1, was conducted. The comparisons between the experimental and desired results indicated by simulations are shown in Figure 10. As shown in Figures 10(a) and 10(b), the experimental responses of the end-effector and joints are approximately the same as the desired response, and the steady-state errors are approximately zero, indicating the accuracy and validity of the dynamic modeling and impedance control. Figure 10(c) shows the simulated and experimental motor currents, where and denote the motor currents of motors 1 and 2, which actuate the and , respectively. Consistent with Figures 10(a) and 10(b), the current responses also demonstrate the consistency between the actual dynamic performance and the modeling.

(a)

(b)

(c)

Figure 10: Experimental and desired results during the response from the deviated position to the equilibrium position with impedance control. (a) End-effector response. (b) Joint response. (c) Motor current response.

5.3. Movement Experiments

In order to discuss how impedance interacted with participants and the error pattern during movements, movement experiments were performed. A healthy male subject, who was 23 years old and left-handed, participated the experiments. The subject performed the movements of two representative types, that is, rectilinear repetitions and clockwise circular repetitions. Furthermore, each moment was performed under four conditions, that is, without impedance or resistance, with impedance only, with resistance only, and with impedance and resistance simultaneously. The resistance force was set to 8 N, which was implemented with (15), and the reference line in rectilinear repetition was set as  mm, while the radius of circular movements was 125 mm with respect to the center at  mm. Since the movement speed could affect the accuracy, the repetition frequencies of rectilinear and circular movements were set to 0.5 Hz and 0.35 Hz, respectively. Each experiment lasted 100 seconds, and the interval time between two experiments was 1 hour to eliminate experimental interactions. The experimental protocol was approved by the Ethics Committee of Tsinghua University, Beijing, China.

The trajectories, errors, and interactive force of rectilinear and circular movements are represented in Figures 11 and 12, respectively. The repetitions in rectilinear movements were indicated by the alternation of (Figure 11), while the force alternations in and directions both suggested the repetitions of circular movements (Figure 12). In circular movements, the signals of and were both sinusoidal with time, and the phrase of lagged behind that of by . The maximal in rectilinear movements and the maximal and in circular moments with resistance were larger than those without resistance by 8 N in average, which also validated the force control. Additionally, the results indicated that fewer errors were observed in the presence of impedance, whereas the performance deteriorated in the presence of resistance.

(a)

(b)

(c)

(d)

Figure 11: Trajectories, errors, and interactive force of the rectilinear movements (a) without impedance or resistance, (b) with impedance only, (c) with resistance only, and (d) with impedance and resistance simultaneously. and error are measured by mm.

Figure 12: Trajectories, errors, and interactive force of the circular movements (a) without impedance or resistance, (b) with impedance only, (c) with resistance only, and (d) with impedance and resistance simultaneously.

To assess the movement accuracy, the mean absolute error (MAE) and standard deviation (SD) were employed to evaluate the deviation of the movements. The brackets indicated the nonsignificant differences (), while the significances at were observed between other groups (Figure 13). The results showed that the rectilinear movement with impedance only had minimal MAE and SD. Specifically, the impedance could significantly decrease the MAEs, whereas the MAEs were significantly increased in the presence of resistance, for both rectilinear and circular movements. Significant differences were also noted between the two groups under the same condition. In addition, the results suggested that the SDs were larger in the absence of impedance and in the presence of resistance, demonstrating that the impedance and resistance mediated the movements by affecting the MAE and SD simultaneously.

Figure 13: Statistics of absolute errors, where the brackets indicate nonsignificant differences (), and the significances at are observed between other groups.

Since the impedance control could reduce the deviations significantly, it was essential to discuss the functional mechanism of the impedance for motor control, which incorporated two sides: theoretical and practical aspects. Theoretically, according to the impedance control proposed by Hogan [29], the main function of the impedance is to determine the interactive force given a deviation, in which case the manipulator is an impedance whereas the environment is an admittance. The generated impedance force is opposite to the deviation from the desired trajectory, which pulls the patient arm towards the desired trajectory as a “virtual damped spring” and alleviates movement deviations as compensative assistance. On the other hand, to interpret the functional mechanism of the impedance in practice, it is imperative to elucidate how the impedance interacts with participants, which is directly reflected by the error patterns in these movements. The kinematic analysis is shown in Figure 14, which indicates the statistic error patterns with respect to the positions. For rectilinear repetitions in direction, the mean errors with respect to positions presented the “arched deviations” within 300 to 550 mm, which reached the maximum at 425 mm (Figure 14(a)). The “arched effect” probably attributed to the inadequate synergy between the shoulder and elbow joints. Particularly, the contribution of the elbow joint motion was more than that of the shoulder joint; thus, the movements tended to present arc trajectories with respect to the elbow joint. As presented in Figure 14(a), the “arched effect” could be alleviated by impedance control. In circular repetitions, the maximal deviation occurred at polar angles of approximately 140 and 300 degrees (Figure 14(b)), which were close to the occasions when the elbow angles reached the maximum and minimum, respectively. Motor performance tended to decrease when close to the joint boundary, and the inadequacy of the elbow angles was supposed to be compensated by shoulder abduction and adduction, which might lead to movement errors but could be alleviated by impedance control. In summary, impedance environment was capable of alleviating movement deviations by compensating the synergetic inadequacy between shoulder and elbow joints, particularly when the movements were close to the joint boundary. Although the efficiency of impedance control regarding dynamic modeling and movement performance was demonstrated by experiments, nevertheless, the physical interaction between robot and human and the contribution to rehabilitation process still required further experiments performed on stroke patients.

(a)

(b)

Figure 14: Statistic error patterns of the rectilinear and circular movements. (a) Mean errors with respect to -position of rectilinear movements. (b) Mean absolute errors with respect to the polar angles of circular movements.

6. Conclusions

This paper presented the design, dynamics, impedance control, and experiments of PARM: a parallel rehabilitation robot using impedance control to enhance interactive training. The parallel mechanism was introduced to reduce the inertia and improve the stiffness, capability, and precision. The motion perception and interaction could be improved by embedding an isometric screen. Apart from the mechanical design, the principle of virtual work and derivative of Jacobian matrix were incorporated to eliminate the frictional and inertial influence. Besides, the strategy based on orthogonal deviation was proposed to achieve the impedance control in assist-as-needed training. Comparisons between desired and experimental responses validated the dynamic modeling and impedance control. To investigate the influence of impedance for movements, movement experiments were also performed. The results showed that the errors of circular movements were mostly larger than those of rectilinear movements and demonstrated the significant differences between movements with/without impedance and resistance (), where the lowest and highest MAEs were noted in the presence of impedance and resistance, respectively. Furthermore, the “arched effect” was observed in rectilinear repetitions, and the deviation tended to occur when the motion was close to the joint boundary, but the impedance environment was capable of alleviating movement deviations by compensating the synergetic inadequacy between the shoulder and elbow joints. For the prospect of robot-assisted therapy, PARM could provide a reference for human-robot interaction in aspects of mechanical design, dynamic modeling, and assist-as-needed control.

Conflicts of Interest

The authors declare that they have no competing financial interest.

Acknowledgments

This research is supported by the National Natural Science Foundation of China (no. U1613207).

References

A. G. Thrift, T. Thayabaranathan, G. Howard et al., “Global stroke statistics,” International Journal of Stroke, vol. 12, no. 1, pp. 13–32, 2017.

D. Mozaffarian, E. J. Benjamin, A. S. Go et al., “Executive summary: heart disease and stroke statistics-2016 update: a report from the American Heart Association,” Circulation, vol. 133, no. 4, pp. 447–454, 2016.

L. H. Thomas, B. French, J. Coupe et al., “Repetitive task training for improving functional ability after stroke,” Stroke, vol. 48, no. 4, pp. e102–e103, 2017.

H. S. Lo and S. Q. Xie, “Exoskeleton robots for upper-limb rehabilitation: state of the art and future prospects,” Medical Engineering & Physics, vol. 34, no. 3, pp. 261–268, 2012.

M. E. Aktan and E. Akdoğan, “Design and control of a diagnosis and treatment aimed robotic platform for wrist and forearm rehabilitation: DIAGNOBOT,” Advances in Mechanical Engineering, vol. 10, no. 1, pp. 1–10, 2018.

M. Zhang, S. Q. Xie, X. Li et al., “Adaptive patient-cooperative control of a compliant ankle rehabilitation robot (CARR) with enhanced training safety,” IEEE Transactions on Industrial Electronics, vol. 65, no. 2, pp. 1398–1407, 2018.

H. I. Krebs, M. Ferraro, S. P. Buerger et al., “Rehabilitation robotics: pilot trial of a spatial extension for MIT-Manus,” Journal of Neuroengineering and Rehabilitation, vol. 1, no. 1, p. 5, 2004.

P. S. Lum, C. G. Burgar, M. V. der Loos, P. C. Shor, M. Majmundar, and R. Yap, “MIME robotic device for upper-limb neurorehabilitation in subacute stroke subjects: a follow-up study,” Journal of Rehabilitation Research and Development, vol. 43, no. 5, pp. 631–642, 2006.

R. Loureiro, F. Amirabdollahian, M. Topping, B. Driessen, and W. Harwin, “Upper limb robot mediated stroke therapy—GENTLE/s approach,” Autonomous Robots, vol. 15, no. 1, pp. 35–51, 2003.

J. C. Perry, J. Rosen, and S. Burns, “Upper-limb powered exoskeleton design,” IEEE/ASME Transactions on Mechatronics, vol. 12, no. 4, pp. 408–417, 2007.

S. Balasubramanian, R. Wei, M. Perez et al., “RUPERT: an exoskeleton robot for assisting rehabilitation of arm functions,” in 2008 Virtual Rehabilitation, pp. 163–167, Vancouver, BC, Canada, August 2008.

J. Hunt, H. Lee, and P. Artemiadis, “A novel shoulder exoskeleton robot using parallel actuation and a passive slip interface,” Journal of Mechanisms and Robotics, vol. 9, no. 1, article 011002, 2017.

T. Nef, M. Guidali, and R. Riener, “ARMin III–arm therapy exoskeleton with an ergonomic shoulder actuation,” Applied Bionics and Biomechanics, vol. 6, no. 2, 142 pages, 2009.

W. H. Chang and Y. H. Kim, “Robot-assisted therapy in stroke rehabilitation,” Journal of Stroke, vol. 15, no. 3, pp. 174–181, 2013.

M. H. Milot, S. J. Spencer, V. Chan et al., “A crossover pilot study evaluating the functional outcomes of two different types of robotic movement training in chronic stroke survivors using the arm exoskeleton bones,” Journal of Neuroengineering and Rehabilitation, vol. 10, no. 1, p. 112, 2013.

R. Bertani, C. Melegari, M. C. De Cola, A. Bramanti, P. Bramanti, and R. S. Calabrò, “Effects of robot-assisted upper limb rehabilitation in stroke patients: a systematic review with meta-analysis,” Neurological Sciences, vol. 38, no. 9, pp. 1561–1569, 2017.

J. Mehrholz, A. Hädrich, T. Platz, J. Kugler, and M. Pohl, “Electromechanical and robot-assisted arm training for improving generic activities of daily living, arm function, and arm muscle strength after stroke,” Cochrane Database of Systematic Reviews, no. 6, article CD006876, 2012.

J. M. Belda-Lois, S. Mena-del Horno, I. Bermejo-Bosch et al., “Rehabilitation of gait after stroke: a review towards a top-down approach,” Journal of Neuroengineering and Rehabilitation, vol. 8, no. 1, p. 66, 2011.

Z. Qian and Z. Bi, “Recent development of rehabilitation robots,” Advances in Mechanical Engineering, vol. 7, no. 2, Article ID 563062, pp. 1–11, 2014.

R. A. R. C. Gopura, D. S. V. Bandara, K. Kiguchi, and G. K. I. Mann, “Developments in hardware systems of active upper-limb exoskeleton robots: a review,” Robotics and Autonomous Systems, vol. 75, pp. 203–220, 2016.

Y. Akiyama, Y. Yamada, and S. Okamoto, “Interaction forces beneath cuffs of physical assistant robots and their motion-based estimation,” Advanced Robotics, vol. 29, no. 20, pp. 1315–1329, 2015.

M. Luces, J. K. Mills, and B. Benhabib, “A review of redundant parallel kinematic mechanisms,” Journal of Intelligent & Robotic Systems, vol. 86, no. 2, pp. 175–198, 2017.

X. J. Liu and J. Wang, “A new methodology for optimal kinematic design of parallel mechanisms,” Mechanism and Machine Theory, vol. 42, no. 9, pp. 1210–1224, 2007.

J. Klein, S. Spencer, J. Allington, J. E. Bobrow, and D. J. Reinkensmeyer, “Optimization of a parallel shoulder mechanism to achieve a high-force, low-mass, robotic-arm exoskeleton,” IEEE Transactions on Robotics, vol. 26, no. 4, pp. 710–715, 2010.

A. Gupta, M. K. O'Malley, V. Patoglu, and C. Burgar, “Design, control and performance of RiceWrist: a force feedback wrist exoskeleton for rehabilitation and training,” The International Journal of Robotics Research, vol. 27, no. 2, pp. 233–251, 2008.

J. Li, S. Li, L. Zhang, C. Tao, and R. Ji, “Position solution and kinematic interference analysis of a novel parallel hip-assistive mechanism,” Mechanism and Machine Theory, vol. 120, pp. 265–287, 2018.

H. Taheri, S. A. Goodwin, J. A. Tigue, J. C. Perry, and E. T. Wolbrecht, “Design and optimization of partner: a parallel actuated robotic trainer for neurorehabilitation,” in 2016 38th Annual International Conference of the IEEE Engineering in Medicine and Biology Society (EMBC), pp. 2128–2132, Orlando, FL, USA, August 2016.

D. Colombo, S. Serino, C. Tuena et al., “Egocentric and allocentric spatial reference frames in aging: a systematic review,” Neuroscience & Biobehavioral Reviews, vol. 80, pp. 605–621, 2017.

N. Hogan, “Impedance control: an approach to manipulation: part I—theory,” Journal of Dynamic Systems, Measurement, and Control, vol. 107, no. 1, p. 1, 1985.