logo资料库

iLeg-下肢康复机器人:概念验证.pdf

第1页 / 共8页
第2页 / 共8页
第3页 / 共8页
第4页 / 共8页
第5页 / 共8页
第6页 / 共8页
第7页 / 共8页
第8页 / 共8页
资料共8页,全文预览结束
IEEE TRANSACTIONS ON HUMAN-MACHINE SYSTEMS, VOL. 46, NO. 5, OCTOBER 2016 761 iLeg—A Lower Limb Rehabilitation Robot: A Proof of Concept Feng Zhang, Zeng-Guang Hou, Long Cheng, Weiqun Wang, Yixiong Chen, Jin Hu, Liang Peng, and Hongbo Wang Abstract—In this paper, a robot, namely iLeg, is designed for the pur- pose of rehabilitation of patients with hemiplegia or paraplegia. The iLeg is composed of one reclining seat and two leg orthoses, and each leg ortho- sis has three degrees of freedom, which correspond to the hip, knee, and ankle. Based on this robotic system, two controllers, i.e., passive training controller and active training controller, are proposed. The former takes ad- vantage of the proportional-integral control method to solve the trajectory tracking problem, and the latter employs the surface electromyography sig- nals to achieve active training. Two simplified impedance controllers, i.e., damping-type velocity controller and spring-type position controller, are designed for active training. A perceptron neural network detects move- ment intentions. The performance of the controllers was investigated with one able-bodied male. The results showed that the leg orthosis tracked the predefined trajectory based on the passive training controller, with the er- ror rates of 0.45%, 0.44%, and 0.27%, respectively, for the hip, knee, and ankle. The active training controller whose loop rate is 6.67 Hz can move the leg orthosis smoothly, and the average recognition error of the perceptron neural network is less than 5%. Index electromyography (EMG), Terms—Active training, rehabilitation robot, spinal cord injury (SCI). I. INTRODUCTION S TROKE and spinal cord injury (SCI) often lead to long- term limb dysfunctions, especially hemiplegia and para- plegia. To improve neurorehabilitation and motor recovery, and to avoid disuse atrophy of the lower limbs, repetitive and inten- sive rehabilitation exercises with the disabled limbs are indis- pensable. Traditional rehabilitation exercises that are conducted manually with the help of therapists are not only labor intensive, but also very costly. Rehabilitation bikes are often used in China to help stroke or SCI patients do treadmill training, which has a positive role on activities of daily living in convalescent stroke patients [1]. The lower limb rehabilitation robot (LLRR) could provide a more flexible alternative to rehabilitation bikes. Manuscript received December 17, 2014; revised April 25, 2015, October 21, 2015, and January 10, 2016; accepted April 8, 2016. Date of publica- tion July 7, 2016; date of current version September 14, 2016. This work was supported in part by the National Natural Science Foundation of China under Grant 61225017, Grant 61403378, Grant 61422310, Grant 61421004, and Grant 61533016, and by the Strategic Priority Research Program of the Chinese Academy of Sciences under Grant XDB02080000. This paper was rec- ommended by Associate Editor J. L. Contreras-Vidal. (Corresponding author: Zeng-Guang Hou.) F. Zhang, Z.-G. Hou, L. Cheng, W. Wang, Y. Chen, J. Hu, and L. Peng are with the State Key Laboratory of Management and Control for Complex Sys- tems, Institute of Automation, Chinese Academy of Sciences, Beijing 100190, China (e-mail: feng.zhang@ia.ac.cn; zengguang.hou@ia.ac.cn; long.cheng@ ia.ac.cn; weiqun.wang@ia.ac.cn; yixiong.chen@ia.ac.cn; jin.hu@ia.ac.cn; liang.peng@ia.ac.cn). H. Wang is with the School of Mechanical Engineering, Yanshan University, Qinhuangdao 066004, China (e-mail: hongbo_w@ysu.edu.cn). Color versions of one or more of the figures in this paper are available online at http://ieeexplore.ieee.org. Digital Object Identifier 10.1109/THMS.2016.2562510 The LLRRs can be categorized into three types according to their driving mechanisms and locomotor training styles. The first one is the sitting/lying type that features a reclining seat and two leg orthoses, e.g., MotionMaker [2]. It can be used for single- and multiple-joint rehabilitation. The second one is the stand- ing/walking type that usually comprises a step-posture control system and a body weight support system, e.g., Lokomat [3], LOPES [4], LokoHelp [5], ReoAmbulator [6], Gait Trainer [7], HapticWalker [8], and WalkTrainer [9]. This type of LLRR as- sists walking movements of gait-impaired patients and combines intensive functional locomotion therapy with patient assessment and feedback tools. The third one is the wearable exoskeleton type that can be worn by the patient and assists him/her to stand or walk, e.g., Ekso [10], Rex [11], Robot Suit HAL [12], and ReWalk [13]. The effectiveness of LLRR is still an open topic in the current literature. While some studies show their effec- tiveness [3], [15], [16], other studies suggest that robotics are not superior to traditional physical therapy [14]. The sitting/lying-type LLRR could provide a useful alterna- tive to rehabilitation bikes, as it can help patients do not only passive treadmill training, but also other types of training, such as active training or functional electrical stimulation (FES). For example, the MotionMakerTM takes advantage of FES to mimic natural exercise during the rehabilitation [16]. Voluntary or active motor training has been proven more beneficial than passive motor training in eliciting performance improvements and cortical reorganization [17]; thus, an LLRR that can help patients accomplish not only passive training, but also active training would be more helpful. Electromyography (EMG) that contains rich information of muscle activity and voluntary intention can be used for this purpose. For example, Leonard et al. developed a novel EMG-driven hand exoskeleton for bilateral rehabilitation of grasping in stroke, and the EMG signals from the nonparetic hand were used to estimate the grasping force and then replicated as robotic assistance for the paretic hand by means of the hand-exoskeleton [18]. Studies also suggest that substantial motor control information can be extracted from paretic muscles of stroke survivors by EMG signals [19]. In this study, we introduce a novel sitting/lying-type LLRR, namely iLeg. iLeg is intended for handicapped and hemiplegic patients, and its main functions are training and rehabilitation of articular mobility and movement coordination, as well as muscular strength. The purpose of developing iLeg is to re- place rehabilitation bikes and provide more rehabilitation train- ing methods. Two rehabilitation training control methods, i.e., passive training and active training, are proposed. Two sim- plified impedance control methods, i.e., spring-type position control and damping-type velocity control based on EMG sig- nals, are implemented to control the robot’s motion. Finally, 2168-2291 © 2016 IEEE. Personal use is permitted, but republication/redistribution requires IEEE permission. See http://www.ieee.org/publications standards/publications/rights/index.html for more information.
762 IEEE TRANSACTIONS ON HUMAN-MACHINE SYSTEMS, VOL. 46, NO. 5, OCTOBER 2016 Fig. 1. Mechanism of iLeg. (a) Virtual prototype of iLeg. (b) Interior view of left leg orthosis mechanism. (c) Exterior view of left leg orthosis mechanism. (d) Crus length adjusting mechanism. TABLE I ADJUSTMENT RANGES OF iLeg JOINT ANGLE RANGES AND SPECIFICATIONS OF THE MOTORS OF iLeg TABLE II Item Width Thigh link Crus link Ankle link ◦ Item Minimal ( ◦ ) Maximal ( · ) Torque (N m) Speed (r/min) Power (W) Reduction Minimum (mm) Maximum (mm) 400 540 315 395 325 420 90 90 Hip Knee Ankle 0 −120 −15 80 0 35 1.333 0.55 0.09 3200 4860 8600 447 280 79.4 120 180 190 the system performance is investigated in a proof-of-concept evaluation. II. SYSTEM DESIGN A. Mechanical System Design The virtual prototype of iLeg, which consists of one reclining seat, two leg orthoses, two support mechanisms, one control box, and a base, is shown in Fig. 1(a). The modular-design supports assembling, debugging, and maintenance of iLeg. The seat back ◦ and 160 ; can be adjusted to an optional angle between 90 thus, the patient can sit or recline in the seat. Each leg orthosis is mounted on an associated support mechanism, which can move bidirectionally on the base; thus, the distance between two leg orthoses can be adjusted. The distance is defined as the width of iLeg (see Table I). It is adjusted to be consistent with the hip width of the patient during training. The maximum adjustment facilitates getting ON and OFF. ◦ B. Leg Orthosis The two leg orthoses’ mechanical structures are symmetrical. The virtual prototype of the left leg orthosis is shown in Fig. 1(b) and (c). Each leg orthosis has three degrees of freedom, which correspond to the hip, knee, and ankle. The leg orthosis can simulate leg movements such as hip flexion/extension, knee flexion/extension, ankle flexion/extension, and the coordinated movements of the three joints. By considering the physiological movement range of each joint of human leg [20], the movement ranges of iLeg are given in Table II. Each joint of the orthosis is driven by a DC servo motor, and drive belts are used for energy transmission. By using the drive belt, the motors can be located far away from the associ- ated joints; thus, the mechanical structures of the thigh and crus links can be smaller. Mechanical interference between the drive units and robot links is eliminated. The ankle motor (32SYK60 of Saegmotory Co., Ltd., Shanghai, China) is located in the pedal, and its output power is delivered to the ankle through a driving belt which is hidden in the ankle link [see Fig. 1(b)]. This addresses the limited space in the ankle. The hip and knee joints are also driven this way. The drive unit of the knee is located behind the hip joint. This novel design not only reduces the size and complexity of the joint, but also offsets some weight of the leg orthosis. The knee motor (48SYK91 of Saegmotory Co., Ltd., Shanghai, China) speed of the orthosis is reduced through the harmonic transmission and multistage gear mechanism, in- creasing output torque. The drive unit of the hip is placed on the associated support mechanism to simplify the mechanical structure of the orthosis, as shown in Fig. 1(a). The hip motor (ID33004 of MCG, Inc., MN, USA) speed is reduced through a harmonic speed reducer, and then, the power is delivered to the hip through a driving belt. The specifications of the motors of iLeg are given in Table II. C. Adjusting Mechanism The link length of the leg orthosis should be adjustable to fit patients of different height (150–190 cm). iLeg can adjust lengths of the thigh and crus by the associated electric length
IEEE TRANSACTIONS ON HUMAN-MACHINE SYSTEMS, VOL. 46, NO. 5, OCTOBER 2016 763 where x and y denote the coordinates of point P , lthigh and lcrus are the lengths of thigh and crus links, and θhip and θknee are the joint angles of hip and knee, respectively. According to (1), we obtain the inverse kinematics of leg orthosis as θhip = arctan y x θknee = −arccos l2 thigh + arccos x2 + y2 − l2 thigh 2lthigh lcrus − l2 crus + x2 + y2 2lthigh − l2 x2 + y2 crus . (2) Fig. 2. Kinematic model of the leg orthosis. B. Passive Training Controller By (2), we can obtain the joint trajectories once the trajectory of point P is known. adjusting mechanism; thus, the leg orthosis somewhat matches the human leg. According to the respective proportion of thigh and crus lengths to body height [21], the designed adjustment ranges of each link are given in Table I. The length adjusting mechanisms are implemented through telescopic mechanisms, and each link can be controlled in a sliding groove to change the length. The lead screw that is driven by a stepping motor and associated screw nut is used to drive the mechanisms. The prototype of the drive unit of the crus length adjusting mecha- nism is shown in Fig. 1(d). The encoder in Fig. 1(d) is used to measure the link length, and it moves together with the screw nut along the rack. Limit switches are used to detect the limit positions. One end of the crus length adjusting mechanism is fixed to the knee joint, and the other end is connected to the crus link through the screw nut. The thigh length adjusting mecha- nism uses a similar idea, although it is implemented differently. By using this driving method, the length adjusting mechanisms move stably and smoothly. The ankle link shown in Fig. 1(b) is the segment between the ankle joint and the pedal, and it corresponds to the patient’s foot height. As there is no adjust- ing mechanism, layers on the pedal can support patients with different foot heights. III. CONTROL METHODS A. Kinematics The leg orthosis has three serial links: thigh, crus, and ankle [see Fig. 1(b)]. The kinematic model of the leg orthosis is shown ◦ when the thigh is parallel to in Fig. 2. The hip joint angle is 0 the seat plane. Trajectory planning of the leg orthosis in the Cartesian plane usually refers to the trajectory of point P which is at the end of the crus link. The movement of the ankle joint is often planned according to specific applications. For example, for passive cycling training, the trajectory of point P is a circle in the Cartesian plane, while the ankle joint can keep still during the cycling movement. Thus, a simplified two-link model can be used for the kinematic analysis of the leg orthosis. We obtain the forward kinematics of leg orthosis as x = lthigh cos(θhip) + lcrus cos(θhip + θknee) y = lthigh sin(θhip) + lcrus sin(θhip + θknee) (1) The passive training controller solves the trajectory tracking problem of the leg orthosis. A proportional-integral (PI) con- troller is designed in the motor driver layer for this purpose. The PI controller is composed of the current loop, velocity loop, and position loop (see Fig. 3). The position loop is in the outer layer; it receives position command θd from the trajectory generator, and calculates the instantaneous profile velocity ˙θd and acceleration ¨θd. These signals along with the actual position feedback θ are processed by the position loop to generate a velocity command vc. From Fig. 3, we obtain the velocity command vc as ˙θd + Af vc = Pp(θd − θ) + Vf ¨θd (3) where Pp, Vf , and Af are the proportional gain, velocity feed forward gain, and acceleration feed forward gain, respectively. The primary effect of Pp is reducing the rise time, Vf is reduc- ing tracking error during constant velocity, and Af is reducing tracking error during acceleration and deceleration. The velocity loop in the middle layer accepts the velocity command which is generated by the position loop, subtracts the actual velocity v and produces a velocity error signal Δv. The error signal is processed by using the integral and proportional gains to produce a current command ic: ic = Vp(vc − v) + Vi t (vc − v)dt (4) 0 where Vp and Vi are the velocity loop proportional and integral gains, respectively. The current loop is in the inner layer, and it is similar to the velocity loop in structure. This loop tracks the current command ic by adjusting the PWM command. The current i generated by this loop can be described by i = Cp(ic − i) + Ci 0 t (ic − i)dt (5) where Cp and Ci are the current loop proportional and integral gains, respectively. The parameters in the PI controller are tuned manually in a specific order: current loop, velocity loop, position loop (see Table III).
764 IEEE TRANSACTIONS ON HUMAN-MACHINE SYSTEMS, VOL. 46, NO. 5, OCTOBER 2016 Fig. 3. Structure of the PI controller. TABLE III PARAMETERS OF THE PI CONTROLLER Joint Hip Knee Ankle P P V f A f 1800 1500 1300 6000 5200 4000 16384 15124 12571 V p 400 500 300 V i 10 20 15 C p C i 68 100 200 138 120 158 C. Active Training Controller 1) Controller Design: For active training, the robot trajec- tories are controlled by the patient. The robot, iLeg, detects the patient’s motion intention and assists him/her to achieve it. Impedance control is a popular method for active control of a manipulator’s interactive behavior [22], [24]. The basic idea of impedance control is to establish a mass–damper–spring rela- tionship between the Cartesian position Δx and the Cartesian torque/force τ as τ = MΔ¨x + DΔ ˙x + KΔx (6) with the positive-definite matrices M, D, K representing the virtual inertia, damping, and stiffness of the interactive system, respectively. These matrices can be chosen by the control system designer according to the goals to be achieved by the robot; thus, impedance control does not attempt to track motion and force trajectories but to regulate the mechanical impedance specified by a target model [24]. Force/torque signals are important but not essential for impedance control. Due to the virtuality of M, D, and K, the Cartesian torque/force τ can also be replaced with other signals. Surface electromyography (sEMG) is a physiological signal, which represents the muscle strength directly; thus, it is a suitable alternative for obtaining the Cartesian torque/force. The human joint flexion and extension movements are con- trolled by two muscle groups, i.e., agonist and antagonist, whose roles are exchanged when the movement changes. Thus, sEMG signals from the two muscle groups should be acquired simulta- neously, and the role of each muscle group must be recognized firstly. In this paper, the normalized amplitudes of sEMG are directly used as the Cartesian forces. Finally, we obtain the sEMG-based impedance controller for active training as CA = MΔ¨θ + DΔ ˙θ + KΔθ ⎡ ⎢ ⎣ c1 −c2 0 0 0 0 C = 0 0 c3 −c4 0 0 0 0 0 0 c5 −c6 ⎤ ⎥ ⎦ Fig. 4. Damping-type velocity controller. a1 θ1 A = θ = a2 θ2 a3 θ3 T T a4 a5 a6 (7) where C is a matrix classifying the agonist and antagonist of each joint; ci, i = 1, ..., 6, are logical variables (0 or 1), and they correspond to the hip flexor, hip extensor, knee extensor, knee flexor, ankle flexor, and ankle extensor, respectively. If ci = 1, it means that the corresponding muscle group is the agonist; otherwise, it is the antagonist. The flexor and extensor of each joint cannot be agonists at the same time. A is a vector, and ai, i = 1, ..., 6, represent the normalized sEMG amplitudes of the muscles. θ is a vector, and θi, i = 1, 2, 3, represent joint angles of the hip, knee, and ankle. The matrices M, D, and K affect the maximum acceler- ation, maximum speed, and maximum position, respectively. Two simplified impedance controllers have been designed: one is the damping-type velocity controller, and the other is the spring-type position controller. The first controller is obtained by setting the matrices M and K in (7) to zero matrices, and the sEMG amplitudes are converted into joint angular velocities. Thus, the patient can drive the leg orthosis to the specific posi- tions at different speeds by contracting the associated muscles. The control law of this method can be written as CA = D( ˙θd − ˙θ0) (8) where D is a positive-definite diagonal matrix, ˙θd is the desired joint angular speed, and ˙θ0 is the reference joint angular speed. Traditionally, ˙θ0 is set to a zero vector. The damping parameters in D can be adjusted to fit patients with different levels of muscle strength. When the matrix D is larger, the robot will be less sensitive to changes in sEMG. It means that more sEMG changes are necessary to achieve the same speed. The control block diagram of the damping-type velocity controller is shown in Fig. 4, from which we obtain the input of the PI controller θd as t t θd = = 0 0 ˙θd dt = 0 t ( ˙θf + ˙θ0)dt −1CA + ˙θ0)dt (D (9)
IEEE TRANSACTIONS ON HUMAN-MACHINE SYSTEMS, VOL. 46, NO. 5, OCTOBER 2016 765 Fig. 5. Spring-type position controller. gi = g(yi) = g ⎝ Each perceptron has three inputs: sEMG amplitudes of the corresponding flexor and extensor, and the threshold value bi, i = 1, 2, 3. The output of each perceptron can be written as ⎛ 2i j=2i−1 ⎞ wj aj + bi ⎠ ⎛ ⎝sgn ⎛ ⎝ 2i j=2i−1 ⎞ ⎠ + 1 ⎞ ⎠ wj aj + bi = 1 2 1, 0, yi > 0 yi ≤ 0, = (12) The following two linear functions f1(·) and f2(·) are used to obtain the outputs of the network: i = 1, 2, 3. Fig. 6. Structure of the neural network. where ˙θf is the speed command generated by the sEMG signals. Similarly, the second controller is obtained by setting the matrices M and D in (7) to zero matrices, and the sEMG ampli- tudes are converted into joint angles. There is a reference angle for each joint. The orthosis deviates from the reference angle when the subject contracts the relevant muscles, and the orthosis moves back to the reference angle when the relevant muscles are relaxed. The leg orthosis works like a spring, and the control law of this method can be written as CA = K(θd − θ0) (10) where K is a positive-definite diagonal matrix, θd is the desired joint angle, and θ0 is the reference joint angle. The stiffness parameters in K can be adjusted to fit patients with different levels of muscle strength. Similarly, when the matrix K is larger, the robot will be less sensitive to changes in sEMG. It means that more sEMG changes are necessary to achieve the same difference between θd and θ0. The control block diagram of the spring-type position controller is shown in Fig. 5. According to Fig. 5 and (10), we obtain the input of the PI controller θd as θd = θf + θ0 = K −1CA + θ0 (11) where θf is the error between θd and θ0, which is generated by the sEMG signals. The patient is in loop of the controllers, and he/she can control the position or velocity of the leg orthosis voluntarily by con- tracting the agonists (see Figs. 4 and 5). Both controllers can motivate the patient to actively participate in the rehabilitation by setting appropriate training goals. 2) Agonist Recognition: Since the movements of the leg or- thosis are controlled by sEMG from the agonists, the agonist recognition, namely computation of the matrix C in (7)–(11) must be completed first. Considering the independence of the joint movements, a neural network (see Fig. 6) that consists of three perceptrons with the same structure is designed for the agonist recognition. f1(x) = x, f2(x) = 1 − x, x = 0, 1 x = 0, 1. (13) IV. PROOF-OF-CONCEPT PROTOTYPE EVALUATION A proof-of-conceptual prototype evaluation was conducted with an able-bodied 30-year-old 175-cm-tall male. The study was approved by the Institutional Review Board of China Re- habilitation Center (Beijing, China). The experiments include passive training and active training, and all were conducted un- der the guidance of physician. A. Passive Training The purpose of passive training evaluation is to test the trajec- tory tracking performance of the PI controller. Because cycling training can dramatically improve aerobic capacity and func- tional performance of stroke patients [25], [26], the leg orthosis is controlled to follow a predefined cycling movement in this experiment. For cycling movement, the trajectory of point P (x, y), shown in Fig. 2, is a circle, which can be described as x = xc + r cos (ωt) y = yc + r sin (ωt) (14) with (xc, yc) the center’s coordinates, r the radius, and ω the ◦ between angular velocity. There is a phase difference of 180 two leg orthoses. By substituting (14) into (2), we can obtain trajectories of the hip and knee joints. The trajectory of ankle joint is planned as that of the hip, except that the amplitude is reduced by half. Such planning can result in the ankle flexion when the ankle approaches the hip, and ankle extension when the ankle leaves the hip. The thigh and crus links lengths were adjusted to the subject’s height, in this case, lthigh = 48 cm, lcrus = 42 cm. In this ex- periment, (xc , yc) = (65, 0) cm, r = 15 cm, and ω = −0.5 π/s. Trajectory tracking performance is summarized in Fig. 7. The predefined and actual trajectories of each joint are shown in Fig. 7(a); each joint moves smoothly and follows the predefined
766 IEEE TRANSACTIONS ON HUMAN-MACHINE SYSTEMS, VOL. 46, NO. 5, OCTOBER 2016 Fig. 7. Trajectory tracking performance of the passive training controller. (a) Actual (dashed lines) and predefined (solid lines) trajectories of each joint. (b) Tracking errors of each joint. (c) Tracking performance in the Cartesian plane. TRACKING ERRORS OF THE CYCLING MOVEMENT DATA USED TO TRAIN AND TEST THE NEURAL NETWORK TABLE IV TABLE V Joint Hip Knee Ankle ◦ A m i n ( ) ◦ A m a x ( ) ◦ E m a x ( ) E rate (%) 21.68 −112.89 0.84 58.29 −54.66 19.14 0.16 0.26 0.05 0.45 0.44 0.27 Item TF/GM RF/BF TA/GA ◦ ) θh i p ( ◦ θk n e e ( ◦ θankle ( Times ) ) 0 −90 0 20 30 −90 0 20 60 −90 0 20 30 −30 0 20 30 −60 0 20 30 −90 0 20 30 −90 −10 20 30 −90 0 20 30 −90 20 20 trajectory. The curves of tracking errors are shown in Fig. 7(b). Fig. 7(c) shows that the actual trajectory is an approximate el- lipse which is very close to the desired circle in the Cartesian plane. The specific tracking errors are summarized in Table IV, where Amin, Amax, Emax, and Erate denote the minimum angle, maximum angle, maximum error, and error rate, respectively. The error rate Erate is defined by Erate = Emax Amax − Amin × 100%. (15) Although the knee joint obtains the largest maximum error of 0.26◦ , its tracking performance is still better than that of the hip joint due to a smaller error rate (see Table IV). Besides the parameters of the PI controller, the tracking error is also influenced by the joint moment of inertia. That is why, the ankle joint has the smallest error rate. B. Active Training The purpose of the active training evaluation is to test the agonist recognition rate of the proposed neural network, and to test the movement characteristics of the leg orthoses con- trolled by sEMG signals. sEMG signals from the tensor fasciae latae (TF), gluteus maximus (GM), rectus femoris (RF), biceps femoris (BF), tibialis anterior (TA), and gastrocnemius (GA) were acquired for this training evaluation. 1) Neural Network Training: Only the right leg was used for testing the agonist recognition. The subject sat in iLeg with the seat back adjusted to 150◦ . Both legs were attached to the leg orthoses with strips of velcro, and the left leg kept a relaxed pos- ture. The subject was asked to contract the selected muscles at three different postures (see Table V). Each muscle contraction lasted for 1–2 s and was repeated 60 times. The subject rested SPECIFIC PARAMETERS OF THE NEURAL NETWORK TABLE VI w 1 0.15 w 2 −0.11 w 3 0.31 w 4 −0.20 w 5 0.25 w 6 −0.30 b1 0.06 b2 −0.05 b3 0.03 for 30 s every ten repetitions. Due to the leg orthosis keeping its posture during each repetition, the subject could only contract the leg muscles without moving the leg. All sEMG signals were preprocessed before use (filtering, full-wave rectification, and smoothing, as in [27]). After pre- processing, sEMG signals were normalized by dividing the maximum amplitudes, which were acquired under the muscles’ maximum voluntary contractile force. As there were 60 samples for flexion or extension of each joint, 120 samples were used for each perceptron. They were divided in half randomly. One part was used for neural network training, and the other was used for neural network validation. The delta rule was used to update the weights of the proposed neural network [23]. The parameters of the neural network after training are given in Table VI. The trained neural network was validated by using the samples for validation, and the agonist recognition errors are shown in Table VII. The average recognition errors Emean were 5.71% ± 1.30%, 4.05% ± 1.07%, and 3.06% ± 1.05%, respec- tively, for the hip, knee, and ankle. The errors were mainly caused by cocontraction of the antagonists, which happened when the maximum sEMG amplitudes were less than 30% of MVC. 2) Spring-Type Position Controller: Only the experiment on the knee joint is presented. According to (10), the reference joint angle θ0 should be set firstly. To ensure the range of motion of
IEEE TRANSACTIONS ON HUMAN-MACHINE SYSTEMS, VOL. 46, NO. 5, OCTOBER 2016 767 TABLE VII AGONIST RECOGNITION ERROR Joint Posture 1 (%) Posture 2 (%) Posture 3 (%) Hip Knee Ankle 7.15 5.26 3.52 5.36 3.65 3.80 4.62 3.23 1.85 E mean (%) 5.71 ± 1.30 4.05 ± 1.07 3.06 ± 1.05 Fig. 9. Results of the damping-type velocity controller. (a) sEMG amplitudes of RF. (b) sEMG amplitudes of BF. (c) Knee joint angles of the leg orthosis. to the matrix K and sEMG changes. The response time is always less than 1/fc. If the current movement is not finished in 1/fc, the controller will stop it and start the next loop. Fig. 8. Results of the spring-type position controller. (a) sEMG amplitudes of RF. (b) sEMG amplitudes of BF. (c) Knee joint angles of the leg orthosis. the knee joint in both directions, it is better to set θ0 near the middle of the knee joint. In this experiment, θ0 = −46◦ , which means that the knee joint will remain in this position unless the agonists are activated. Due to normalization of the sEMG amplitudes, each element of vector A in (10) is between 0 and 1; thus, the maximum angle from the reference joint angle is determined by the matrix K. When tuning K, we need to make sure that θd in (11) is in the range of motion of the leg orthosis. In this experiment, K = 0.04, which means that the maximum angle deviation from the reference joint angle is K −1 = 25◦ . The results are shown in Fig. 8. At t0, the RF was recognized as the agonist. The knee began to extend, and the extension movement led to the increase of the joint angle. At t1, the sEMG amplitude of RF reached the maximum value. The extension movement did not stop until t2 because the desired joint angle was always larger than the current joint angle between t1 and t2. After t2, the knee began to move toward the reference angle due to the decrease of the sEMG amplitude of RF. Similarly, at t3, the BF was recognized as the agonist, and the knee began to flex. Although the sEMG amplitude of BF began to decrease after t4, the flexion movement did not stop until t5. After t5, the knee joint moved toward the reference angle when the sEMG amplitude of BF continued to decrease. Fig. 8 shows that there is a delay of about 0.3 s between the control signals and the desired joint angle of the leg orthosis. The delay is mainly caused by the control loop rate fc and the response time of the leg orthosis. The control loop rate is also the rate that the controller processes the sEMG signals. Due to the instability of sEMG signals, the leg orthosis tends to oscillate if the loop rate is too high. When tuning the loop rate, we should consider the tradeoff between the delay and oscillation. In this experiment, we get a feasible fixed loop rate of 6.67 Hz after repeated tests. The response time is also the time that the leg orthosis would take to finish the movement. Thus, it is relevant 3) Damping-Type Velocity Controller: Only data from the knee joint are presented. According to (8) and (9), the sEMG amplitudes are used to control the angular velocity of the leg orthosis, and the maximum angular velocity is determined by −1 the matrix D. When tuning D, we need to make sure that D does not exceed the maximum angular speed of the leg orthosis. Traditionally, we set the values of D according to the subject’s subjective feelings. D = 0.0625; thus, we obtain the maximum angular velocity as D −1 = 16 ◦ /s. The results are shown in Fig. 9. Different from the spring- type position controller, this controller moved the leg orthosis continuously in the same direction until the agonists were ex- changed or until the agonists were relaxed. The control loop rate of this controller was also 6.67 Hz. As shown in Fig. 9, the RF was recognized as the agonist at t0, and the knee was controlled to move toward the direction of extension. The joint angular velocity changed with the sEMG amplitudes of RF dur- ing the movement, and the movement did not stop until the RF was relaxed (at t1). Similarly, the BF was recognized as the agonist at t2, and the knee was controlled to move toward the direction of flexion. This movement stopped at t3 when the BF was relaxed. V. DISCUSSION In this paper, we have introduced a novel sitting/lying-type LLRR, namely iLeg. iLeg has a similar appearance to Motion- Maker, but the mechanical structures and drive mechanisms of the leg orthosis are quite different. iLeg uses drive belts for en- ergy transmission, whereas MotionMaker uses lead screw for energy transmission. By using the drive belts, the moment arm of each joint will remain unchanged during the joint movement. Besides, iLeg can adjust lengths of the thigh and crus by the as- sociated electric length adjusting mechanism, which makes the leg orthosis more flexible. The purpose of developing iLeg is to replace rehabilitation bikes, which are widely used in China for training and rehabilitation of articular mobility and movement coordination, and provide more rehabilitation training methods for handicapped and hemiplegic patients. There are two main
768 IEEE TRANSACTIONS ON HUMAN-MACHINE SYSTEMS, VOL. 46, NO. 5, OCTOBER 2016 differences between iLeg and the current popular LLRRs. One is the novel mechanism design, and the other is the passive and active training control methods. Thus, iLeg has not only all the functions of a rehabilitation bike, but also can provide sEMG- based active rehabilitation training methods. Because China has the largest population in the world, iLeg is expected to be widely used in China. Compared with the standing/walking-type LLRRs, such as the Lokomat or LOPES, the training trajectories of iLeg are more diverse, such as cycling movement. This novel system may help patients perform uniarticular training, which provides convenience for targeted joint rehabilitation, such as foot drop. The sEMG-based impedance controllers detect not only the subject’s movement intention, but also voluntary participation, which may encourage a patient to move his/her legs. Compared to existing techniques, the active training controllers proposed in this paper do not need the user to remember the relationship between the muscle contraction and the motion. Moreover, since the position or velocity of the leg orthosis is determined by the sEMG amplitudes, the user may try his/her best to contract the muscles to generate strong sEMG amplitudes during each rehabilitation session. Such voluntary participation could make the rehabilitation training more effective. Muscle spasm should be considered seriously during active rehabilitation. The sEMG should be treated as abnormal signal when muscle spasm happens, and the controllers should stop iLeg immediately to avoid secondary damage to the subject. Although studies have suggested that substantial motor control information can be extracted from paretic muscles of stroke survivors by EMG signals [19], there are still many challenges when impaired subjects try and use iLeg. People with stroke and SCI typically are weak and have relatively small EMG amplitudes and abnormal muscle coordination patterns that may make it difficult for them to control iLeg. To deal with these challenges, we should focus on two aspects. One is choosing the appropriate EMG signals that can describe the motion intention of patients, and the other is improving the motion intention recognition algorithm. In our future work, we will focus on the clinical validation of this system. ACKNOWLEDGMENT The authors are grateful to Dr. Y. Hong, J. Zhang, and Z. Lu with the China Rehabilitation Research Center, Beijing, China, for providing suggestions. REFERENCES [1] G. Yan, H. Shen, X. Zhao, Q. Wei, Y. Kang, Z. Jia, L. Song, and M. Huang, “Effects of treadmill training on ADL of convalescent stroke patients,” Chin. J. Rehabil., vol. 22, no. 3, pp. 163–164, 2007. [2] C. Schmitt, P. M´etrailler, A. Al-Khodairy, R. Brodard, J. Fournier, M. Bouri, and R. Clavel, “The Motion MakerTM: A rehabilitation system combining an orthosis with closed-loop electrical muscle stimulation,” in Proc. 8th Vienna Int. Workshop Funct. Elect. Stimul., Vienna, Austria, 2004, pp. 117–120. [3] M. A. Maestro, A. E. Ruz, R. M. Casado-Lopez, A. M. Gonzalez, G. P. Mateos, E. G. Valdizan, and J. L. Martin, “Lokomat robotic-assisted versus overground training within 3 to 6 months of incomplete spinal cord lesion: Randomized controlled trial,” Neuralrehabil. Neural Repair, vol. 26, no. 9, pp. 1058–1063, 2012. [4] J. F. Veneman, R. Kruidhof, E. E. G. Hekman, R. Ekkelenkamp, E. H. F. van Asseldonk, and H. van der Kooij, “Design and evaluation of the LOPES exoskeleton robot for interactive gait rehabilitation,” IEEE Trans. Neural Syst. Rehabil. Eng., vol. 15, no. 3, pp. 379–386, Sep. 2007. [5] S. Freivogel, J. Mehrholz, T. Husak-Sotomayor, and D. Schmalohr, “Gait training with the newly developed ‘LokoHelp’-system is feasible for non- ambulatory patients after stroke, spinal cord and brain injury: A feasibility study, Brain Injury, vol. 22, nos. 7/8, pp. 625–632, 2008. [6] G. R. West, “Powered gait orthosis and method of utilizing same,” Patent 6 689 075, 2004. [7] S. Hesse, C. Werner, D. Uhlenbrock, S. Frankenberg, A. Bardeleben, and B. Brandl-Hesse, “An electromechanical gait trainer for restoration of gait in hemiparetic stroke patients: Preliminary results,” Neurorehabil. Neural Repair, vol. 15, no. 1, pp. 39–50, 2001. [8] H. Schmidt, “HapticWalker—A novel haptic device for walking simula- tion,” in Proc. EuroHaptics, Munich Germany, 2004, pp. 66–70. [9] Y. Stauffer, Y. Allemand, M. Bouri, J. Fournier, R. Clavel, P. Metrailler, R. Brodard, and F. Reynard, “The WalkTrainer—A new generation of walking reeducation device combining orthoses and muscle stimulation,” IEEE Trans. Neural Syst. Rehabil. Eng., vol. 17, no. 1, pp. 38–45, Feb. 2009. [10] E. Strickland, “Good-bye, wheelchair,” IEEE Spectr., vol. 49, no. 1, pp. 30–32, Jan. 2012. [11] Rex Bionics REX web site, 2014. [Online]. Available: http://www. rexbionics.com [12] A. Tsukahara, R. Kawanishi, Y. Hasegawa, and Y. Sankai, “Sit-to-stand and stand-to-sit transfer support for complete paraplegic patients with robot suit HAL,” Adv. Robot., vol. 24, no. 11, pp. 1615–1638, 2010. [13] Argo Medical Technologies ReWalk web site, 2014. [Online]. Available: http://www.argomedtec.com [14] B. H. Dobkin and P. W. Duncan, “Should body weight-supported treadmill training and robotic-assistive steppers for locomotor training trot back to the starting gate?” Neurorehabil. Neural Repair, vol. 26, no. 4, pp. 308– 317, 2012. [15] S. Freivoqel, D. Schmalohr, and J. Mehrholz, “Improved walking ability and reduced therapeutic stress with an electromechanical gait device,” J. Rehabil. Med., vol. 41, no. 9, pp. 734–739, 2009. [16] P. Metrailler et al., “Improvement of rehabilitation possibilities with the MotionMakerTM,” in Proc. 1st IEEE/RAS-EMBS Int. Conf. Biomed. Robot. Biomechatron., 2006, pp. 359–364. [17] M. Lotze, C. Braun, N. Birbaumer, S. Anders, and L. G. Cohen, “Mo- tor learning elicited by voluntary drive,” Brain, vol. 126, pp. 866–872, 2003. [18] D. Leonardis et al., “An EMG-controlled robotic hand exoskeleton for bilateral rehabilitation,” IEEE Trans. Haptics, vol. 8, no. 2, pp. 140–151, Apr.–Jun. 2015. [19] X. Zhang and P. Zhou, “High-density myoelectric pattern recogni- tion toward improved stroke rehabilitation,” IEEE Trans. Biomed. Eng., vol. 59, no. 6, pp. 1649–1657, Jun. 2012. [20] A. Roaas and G. B. J. Andersson, “Normal range of motion of the hip, knee and ankle joints in male subjects, 30–40 years of age,” Acta Orthopaedica Scandinavica, vol. 53, pp. 205–208, 1982 [21] Y-C. Lin, M-J. J. Wang, and E. M. Wang, “The comparisons of anthropo- metric characteristics among four peoples in East Asia,” Appl. Ergonom., vol. 35, no. 2, pp. 173–178, 2004. [22] N. Hogan, “Impedance control: An approach to manipulation: Parts I-III,” J. Dyn. Syst. Meas. Control, vol. 107, no. 1, pp. 1–24, 1985. [23] B. Widrow and M. E. Hoff, “Adaptive switching circuits,” in Proc. IRE WESCON Convention Rec., 1960, pp. 96–104. [24] C. C. Cheah and D. Wang, “Learning impedance control for robotic ma- nipulators,” IEEE Trans. Robot. Autom., vol. 14, no. 3, pp. 452–465, Jun. 1998. [25] K. L. Michal, “The influence of early cycling training on balance in stroke patients at the subacute stage: Results of a preliminary trial,” Clin. Rehabil., vol. 20, no. 5, pp. 398–405, 2006. [26] W. J. Thomas, J. M. Beltman, P. Elich, P. A. Koppe, H. Konijnenbelt, A. de Haan, and K. H. Gerrits, “Effects of electric stimulation-assisted cycling training in people with chronic stroke,” Arch. Phys. Med. Rehabil., vol. 89, no. 3, pp. 463–469, 2008. [27] F. Zhang, P. Li, Z.-G. Hou, Z. Lu, Y. Chen, Q. Li, and M. Tan, “sEMG- based continuous estimation of joint angles of human legs by using BP neural network,” Neurocomputing, vol. 78, no. 1, pp. 139–148, 2012.
分享到:
收藏