- Open Access
Control of thumb force using surface functional electrical stimulation and muscle load sharing
© Westerveld et al.; licensee BioMed Central Ltd. 2013
- Received: 16 May 2013
- Accepted: 3 October 2013
- Published: 9 October 2013
Stroke survivors often have difficulties in manipulating objects with their affected hand. Thumb control plays an important role in object manipulation. Surface functional electrical stimulation (FES) can assist movement. We aim to control the 2D thumb force by predicting the sum of individual muscle forces, described by a sigmoidal muscle recruitment curve and a single force direction.
Five able bodied subjects and five stroke subjects were strapped in a custom built setup. The forces perpendicular to the thumb in response to FES applied to three thumb muscles were measured. We evaluated the feasibility of using recruitment curve based force vector maps in predicting output forces. In addition, we developed a closed loop force controller. Load sharing between the three muscles was used to solve the redundancy problem having three actuators to control forces in two dimensions. The thumb force was controlled towards target forces of 0.5 N and 1.0 N in multiple directions within the individual’s thumb work space. Hereby, the possibilities to use these force vector maps and the load sharing approach in feed forward and feedback force control were explored.
The force vector prediction of the obtained model had small RMS errors with respect to the actual measured force vectors (0.22±0.17 N for the healthy subjects; 0.17±0.13 N for the stroke subjects). The stroke subjects showed a limited work range due to limited force production of the individual muscles. Performance of feed forward control without feedback, was better in healthy subjects than in stroke subjects. However, when feedback control was added performances were similar between the two groups. Feedback force control lead, especially for the stroke subjects, to a reduction in stationary errors, which improved performance.
Thumb muscle responses to FES can be described by a single force direction and a sigmoidal recruitment curve. Force in desired direction can be generated through load sharing among redundant muscles. The force vector maps are subject specific and also suitable in feedforward and feedback control taking the individual’s available workspace into account. With feedback, more accurate control of muscle force can be achieved.
- Load sharing
- Muscle recruitment
- Force control
Stroke has become a major cause of morbidity and mortality in the western world. Incidence of stroke also increases in less developed countries as a result of changing life-styles . Graying of society and improved health-care are likely to result in an increase of stroke survivors. Functional independence of stroke survivors is highly influenced by their ability to perform a successful grasp. In many activities of daily living, like drinking or opening a door, grasp and release is an essential part of the required movement.
Functional electrical stimulation (FES) of hand muscles can be helpful to train grasp and release in stroke subjects [2–4]. Depending on the ability of the individual patient, the assistance may be (selectively) increased or decreased in order to maximize the voluntary activity which is important in relearning movements [5, 6].
Grasping comprises coordinated finger and thumb motion and controlled force exertion on the object to be held. As muscles initiate human movement, accurate control of muscle force is a prerequisite for movement control. For grasping tasks the fingers can be regarded as single degree of freedom (DoF) joints, since movement of the individual phalanges is coupled because of the under actuation of the finger. Furthermore, rotation along the flexion-extension axis of the finger is by far the most important movement for grasping and releasing objects. The thumb, however, requires a different approach as it moves along multiple axes. Controlling force and movement of the thumb will be most challenging and may serve as a model, which may be generalized/reduced to the single DoF case for the other fingers.
A healthy thumb is actuated in several directions by nine muscles in total [7, 8]. However, not all nine muscles can be targeted properly with surface FES. Mainly, because of overlying muscles and nearby sensory nerves making stimulation uncomfortable. Therefore, only a small subset of thumb muscles is available for FES with surface electrodes. This limits the movements which can be controlled with FES. However, thumb movements relevant for grasping (mainly opposition) are feasible with surface electrodes.
Force distribution over multiple muscles is commonly applied in biomechanical modelling, solving actuator redundancy problems for a given task [9, 10]. This load sharing approach might also be useful for activating a redundant muskuloskeletal system. In addition, by sharing the load over all available muscles we maximize the available range of force. However, to our knowledge, load sharing has not been applied to external activation of muscles with surface electrical stimulation. We will evaluate this possibility and expect this approach to result in accurate force control with a force distribution over the individual muscles optimized by minimizing the sum of squared recruitment over all muscles.
Recently, Lujan et al.  measured thumb forces evoked by three thumb muscles in healthy subjects and one spinal cord injured patient. Using the measured forces they trained an artificial neural network (ANN) for feed forward force control. They showed good control of the isometric thumb force in 2D. With the current study we aim at a more transparent approach: using linear combinations of estimated muscle force vectors instead of using a black-box ANN. This approach gives us the benefit of learning more of the underlying physiological system, by comparing combined muscle responses with individual muscle responses. In addition, it might allow for a more generally applicable approach, without the need of training an ANN.
The goal of the current study is twofold: 1) Is it possible to describe thumb muscle responses to FES by a sigmoidal muscle recruitment curve and a single direction of force? And if so, are these so called muscle force maps subject specific, suitable for stroke subjects and time-invariant? And 2) Are muscle force maps suitable for use in 2D thumb force control with FES applying load sharing? And if so, is feed forward control only sufficient and is the approach also suitable for stroke subjects?
We will introduce the proposed generalized muscle force model for thumb force control and muscle load sharing first. Thereafter we will describe the experimental evaluation of this model in both healthy subjects and stroke subjects.
Generalized muscle force model
Feedforward thumb force model
In Eq. 2, the predicted thumb force vector , is the vector sum of the individual muscle forces (n=3), modelled as a recruitment fraction, x i , of the maximal muscle force magnitudes, .
The combination of obtained stimulation amplitudes, A i , is the combination which theoretically would produce a force equal to the reference force, , or at least the force which is minimizing Eq. 3 when the system has reached its boundaries of operation. The constant C represents the offset term as introduced in Eq. 1.
Characteristics of included stroke subjects
Either the dominant arm (healthy subjects) or the affected arm (stroke subjects) was strapped in a custom built device. This setup was used to fixate the wrist and the hand in neutral pronosupination, and to measure the isometric thumb force in two directions perpendicular to the axis of the thumb. Forces were measured by two 45.3 N load cells (Futek, Irvine) preloaded with springs. See Figure 1.
A special built 3 channel asynchronous biphasic electrical stimulator (TIC Medizin, Dorsten, Germany) was used to apply the electrical stimulation pattern. Stimulation was applied at a constant frequency (30 Hz) and pulse width (150 μ s). The amplitude could be controlled via custom built controllers within the stimulator’s range [0-30m A] in steps of 0.125m A. A single 50×50 m m anode was used together with 16×19 m m cathodes for each channel. Electrodes with similar size showed good results on both selectivity and comfort in a simulation study .
An EtherCAT I/O system (Beckhoff Automation GmbH, Verl, Germany) using Matlab/xPC (The Mathworks, Nattick, USA) as EtherCAT master device was used to control the stimulator parameters and to capture analog data from the force sensors.
Force vector map determination
The subject specific force map (see Figure 1 for an example) was determined in the isometric setup, with the thumb visually positioned at 30 degrees of abduction and 30 degrees of extension. The threshold and maximal stimulation amplitude for each muscle were determined first: we stimulated (30 H z;150 μ s) each muscle individually for 1 second, followed by 0.5 second without stimulation. Every 1.5 second the amplitude was increased by 1m A. When either a saturation in the force response was observed or the subject reported unpleasant discomfort, the stimulation was stopped.
The range between the threshold minus 1m A and the maximal amplitude was divided in ten equidistant stimulation levels for each muscle. We applied these 30 stimulations (10 amplitudes per muscle) randomly and measured the exerted thumb forces.
From this initialization measurement, we determined the force direction of each individual muscle and the recruitment curve relating muscle stimulation to exerted force. The recruitment curves were described with a sigmoidal function having three parameters, using Eq. 1. Parameter values were obtained with a least-squares fit, using the Levenberg-Marquardt algorithm . See Figure 1 for an example of muscle recruitment curves and force directions. This force vector map indicates the ability to control the thumb force in different directions for a specific subject.
Individual muscle controllers
In Eq. 7 the open loop gain, K o , is calculated from the normalized input magnitude, X0, the measured steady state output magnitude, M u , the time until the output responds, τ dead and the time between the first response and the output reaching the steady state, τ. As suggested in , the proportional gain, K c , for each muscle was calculated as 90% of the open loop gain and the integration time for the PI-controller, T i , was set as 3.3 times τ dead . .
For every muscle and subject the inverse of the recruitment curve compensates the non-linear and subject and muscle specific recruitment. In this way the non-linear elements and maximal force levels are compensated within the control loop leading to a linear feedback controller between observed force error and reference force. Furthermore it is expected that range of control gains between the different muscles and different subjects is relatively small, since the muscle and subject specific recruitment curve transforms the outputs of the PI controllers (forces) into the required stimulation amplitudes.
After determining the initial gains for each muscle, in total four single muscle tests were done for each muscle to be able to analyze performances of the individual muscle controllers: 1) step response reference pattern with feedback control, 2) 0.5 Hz sinusoidal reference pattern with feedback controller, 3) step response reference pattern with a combination of feedforward and feedback control, and 4) 0.5 Hz sinusoidal reference pattern with a combination of feedforward and feedback control.
When oscillatory behaviour was observed during the first test, the proportional gain was lowered systematically and the test was repeated until good tracking of the reference was observed without severe oscillations. In some cases the integration time T i was increased slightly for further fine tuning.
2D thumb force targets
For evaluation of the 2D controllers, 5 second constant reference force targets were used. The targets were set at 0.5 N and 1.0 N in different directions within the workspace of the subject. Initially, directions were chosen at –90°, –60°, –30°, 0°, 30° and 60°. Angles outside the theoretical workspace of the subject were not measured. When less than four target directions were theoretically feasible, intermediate angles (15° step size) were also evaluated.
Feedforward thumb force control
The applicability of the thumb force model was evaluated first in an experiment based on feed forward control of the three muscles. In this experiment control was based on the measured muscle parameters and the thumb model described in Eq. 2. Based on the previously determined force map, target angles greater than the angle of the long abductor muscle or smaller. The experiment was repeated three times to explore the reproducibility of the methods. The target sequence was the same in each repetition. The sharing of the load was calculated by implementing Eq. 5 in a real-time quadratic programming (QP) problem solver using the online active set strategy .
Feedforward and Feedback thumb force control
RMS errors were calculated from the magnitude of the error vector between measured muscle force during the initialisation procedure and muscle force estimate based on the obtained parameters. In addition, the area of the theoretical work range resulting from the muscle force vectors obtained during the first initialization procedure was calculated and compared between subjects.
In which F pre and F post , are the observed forces during the procedures at the beginning and the end of the session, respectively. The forces are summed over all applied input amplitudes during the initialization procedure. The slope, m, is an estimate of the ratio between initial force generation and final force generation for a given muscle.
Single muscle control performances were evaluated based on the sine tracking tasks. RMS errors between the actual and reference forces were calculated. The 2D controller performances were evaluated based on the stationary error of the responses. This stationary error was defined as the average magnitude of the force error vector during the last 10 percent of the in total 5 seconds lasting step response.
Due to the relatively small sample size, non-parametric statistics was applied. We used Mann Whitney U tests to statistically evaluate improvement with feedback control over feedforward control only and also to evaluate performance in stroke subjects with respect to healthy subjects.
Force vector maps
F x correlation
F y correlation
F x correlation
F y correlation
Force controller evaluation
Single muscle controllers
The averaged proportional gain over all healthy subjects was 0.22±0.28. For the stroke subjects the average proportional gain was 1.04±1.16, note that these values are dimensionless as the feedback controller has a force both as input and as output, since the inverse recruitment is placed after the controller. The average integral times were 0.56 ± 0.12s and 0.62 ± 0.45s for healthy subjects and stroke subjects respectively.
During the single muscle control experiments, some saturation effects (stimulation reaching predetermined maximal amplitude) were observed, leading to a non-linear feedback system. Disregarding the cases were this saturation occurred, the estimated controller gains were 0.17±0.12 and 0.57±0.12s on average for all subjects for proportional gain and integral time respectively.
Combined muscle controllers
We showed the possibility to describe responses to electrical stimulation of individual thumb muscles as a force vector map with a single activation direction and a sigmoidal recruitment curve. As expected the variability between subjects is relatively large (Figure 4) due to anatomical differences. As a result, force maps always need to be determined for each individual subject. Within subject the results are repeatable, demonstrating the feasibility of our approach (Figure 4 and Table 2). Note that for subsequent sessions it is required to redo the initialization, since the response is highly dependent on exact electrode position . However, in stroke subjects the AbPL muscle was difficult to target. In the subjects in which we were able to target the muscle initially, the responses during the second initialization procedure differed greatly from the initial procedure as indicated by the low correlation coefficients in Table 2 and in Figure 4. Therefore the AbPL muscle seems less reliable for use in 2D force control tasks compared to the other muscles.
The load sharing approach resulted in the muscle being pulled nicely towards the target force by the feedback controller. Since the error vector was used as input for the feedback load sharing, the appropriate ratio of muscle activations was calculated to generate force in the right direction. To our knowledge this load sharing approach is a novel application in electrically stimulated muscle. In our opinion this could be an appropriate solution to solve redundancy problems in activation of multi-dimensional muskuloskeletal systems with FES and simultaneously take the boundaries of the individual force sources into account. The variation of controller gains over different muscles and different subjects was low, which gives the possibility to use fixed values for these parameters when applying the methodology presented here. Either as a true fixed value of as a starting point for further fine tuning instead of the Ziegler-Nichols methods  which were currently used. Thereby further reducing the tunable parameters and setup time.
Performance of the 2D feedforward force controller was worse for the stroke subjects compared to the healthy subjects. For the stroke subjects, adding feedback terms reduced stationary errors. For the healthy subjects the differences between feedforward control only and combined with feedback control were small, see Figure 9. However, depending on the model accuracy of the individual muscle’s input-output relation, the feedback controller also reduced the control performance in certain cases. An example of this can be observed from Figure 7 where the feedback controller negatively influences the force direction for the 0.5 N targets. This is likely a result of a mismatch in the FPB model, causing the thumb being pulled in a more negative direction than needed. Therefore we recommend estimating model accuracy before starting the control trials, and redo the initialization if necessary.
We measured forces in two directions in a plane perpendicular to the thumb. Therefore we neglected the forces perpendicular to this plane. Due to this fact we might have made some errors in absolute force recordings. However, since we are using the same setup in both model identification and control, we expect that the influence of these non-measured forces on our performance observations are minimal.
Forces in unmeasured direction could have led to the relatively low observed forces compared to other studies . However, we expect that these unmeasured forces were small. The stimulated muscles are responsible for thumb movement Therefore the force component in line width the thumb will be small compared to the perpendicular force components. A more likely cause is the fact that we aimed at selective activation with small electrodes leading to relatively low current densities and low muscle activation. Even though the observed forces and the evaluated targets of 0.5 N and 1.0 N are relatively low, they are sufficient for positioning the thumb for functional grasping of objects compared to the evaluated force levels during grasping in [16, 17]. Recently, we have shown applicability of a similar approach during grasp and release of objects .
In all subjects, the FPB muscle showed a steep recruitment curve: when the stimulation came above threshold force increase was high for an increase in stimulation amplitude. This will have resulted in a bigger influence of FPB modelling errors on the output force errors. The steeper recruitment compared to other muscles is likely a result from differences in neural innervation. The FPB muscle is innervated from the recurrent branch of the median nerve which is very superficial before entering the FPB muscle. The OpP muscle is innervated by the same nerve branch, but laterally the branch runs less superficial . The AbPL muscle is innervated by the posterior interosseus nerve which is also less superficial.
We reduced the experiment length by only testing specific points along the recruitment curve during the initialization phase. We did not specifically optimize this method of recruitment curve sampling. However, the results in pilot measurements where we compared our current approach with more dense sampling of the muscle recruitment resulted in only minor differences between the obtained recruitment curves. Recently, Schearer and colleagues  compared different methods of recruitment curve sampling extensively. Application of methods described there might further improve the accuracy of the obtained recruitment curves of individual muscles, which then could also improve the accuracy of the controllers.
The stroke subjects showed smaller workspaces compared to the healthy subjects (Figure 5). This is likely a result of non-use after stroke, which could have been overcome partially by additional muscle training prior to the experiment. However, since we only analyzed performance from the trials where the target force vector was within the theoretical workspace, this has not affected our current findings.
The ARAT scores of the stroke subjects had a broad range. Therefore the subjects cannot be considered as a homogeneous group. However, the emphasizes of the current approach lies on modelling subject specific recruitment relations. Therefore we did not observe lower stimulation responses related to lower ARAT scores. Furthermore, this is supported by the fact that the subjects with the best ARAT scores showed the smallest theoretical work range for the selected muscles.
We expect the remainder of the variation to have a physiological cause. The most likely one is a non-linear additive relation between the individual muscle directions and recruitments. We expect that the linear addition of two individual force magnitudes to produce the desired combined force magnitude had the largest contribution to the remainder of the observed variability.
Lujan and Crago  were able to control the thumb forces in two directions by using an artificial neural network. They also observed differences between the measured force of combined muscle activation and the sum of the individual components, which suggested a non linear additive relation. Lujan and Crago stimulated different muscles (Extensor Pollicis Longus, Abductor Pollicis Brevis and Adductor Pollicis). The evoked forces in that study are about five times higher than the forces which we found, possibly caused by higher stimulation frequencies (50 Hz compared to 30 Hz in our study) and the different set of stimulated muscles. This makes a good comparison between results difficult. Lujan and Crago only report 2D control RMS errors of one healthy subject and one spinal cord injured (SCI) patient, having implanted electrodes. The RMS error of the SCI patient was 0.89 N, which is very low compared to our results in stroke subjects when relating to the achieved force range. However implanted electrodes are known to produce higher muscle selectivity and more direct muscle activation, which makes this comparison unfair. The healthy subject they presented showed an RMS error of 2.65 N, which is (taking the factor 5 into account) within the same range as the stationary RMS errors we observed. However, we were able to obtain that similar performance without training and optimizing an artificial neural network but with a more transparent model consisting of only four parameters per muscle.
Schearer and colleges  recently published a single case study on controlling multiple degrees of freedom (in the shoulder) in a SCI subject with implanted electrodes using a feedforward controller. They also solved for redundancy by using a quadratic program and showed initial RMS errors of 5.29 N. As shoulder muscles are much stronger than thumb muscles, this value is again difficult to compare with our results. Given the range of their target forces (-18 N to 4.5 N in the x direction, -18 N to 22.5 N in the y direction and -9 N to 0 N in the z direction) one could say that the performance of their controller was slightly better than ours, which seems logical given the fact that the electrodes used by Schearer and colleagues were implanted. Therefore their stimulation was likely to result in more selective and accurate activation of individual muscles. In addition, Schearer et al suggest to improve the performance by adding a feedback controller, which is exactly what we did in the current study. We showed that adding the feedback path can indeed improve performance when the feedforward model is not accurate enough.
This study is a framework for evaluating multi-dimensional control of joints with electrical stimulation. To be clinically applicable in post-stroke rehabilitation, the method needs several extensions. First of all, we currently addressed only thumb muscles. For functional grasp and release training the finger muscles are of course equally important. However, compared to the thumb, those joints do not have the redundancy in actuation: mainly one extensor muscle and one superficial flexor muscle. Therefore the current method could easily be extended to the fingers, which we also evaluated recently .
When using additional electrodes for (selective) finger flexion and extension, the number of electrodes will increase quickly. Since, electrode placement is subject dependent and can be time consuming, the time required for setup will also increase rapidly. From a practical point of view, time can be gained with the application of electrode arrays and an approach to automatically search for proper electrode locations .
Finally, the relations between stimulation and movement and control of movement for grasp and release are also important. However proper force control is a fundamental prerequisite for proper control of movement. Therefore the current study can be seen as an intermediate step towards an approach for assisting grasp and release movements and next steps in our research will focus on directly mapping muscle activation to evoked movements.
Stroke subjects showed a limited workspace in our study. Since they did not have severe spasticity, it is likely that their muscle force have decreased dramatically due to long time non-use after their stroke. Therefore, we expect that results in more acute stroke subjects lie closer to those of the healthy subjects in the current experiment. However, this needs further evaluation and likely a subject specific approach will lead to the best results.
The aim of this study was to evaluate the possibility to predict thumb muscle force responses to FES and to control thumb muscle force in 2D in both healthy and stroke subjects. For a single muscle, the static relation between muscle force and activation was described by a sigmoidal muscle recruitment curve and a single direction of force. Subsequently, load sharing was used to combine the activation of individual muscles to actively control thumb force in 2D.
From our results we can conclude that it is possible to describe the thumb muscle responses to FES by a single force direction and a sigmoidal recruitment curve. The large variations between subjects indicate that these force maps are highly subject specific, likely due to anatomical differences, requiring an individual approach. The relatively small variation within subjects demonstrates the feasibility and time-invariance of our approach. Effects of muscle fatigue were observed, especially in stroke patients, so the approach presented here is applicable mainly for short sessions (up to 30 minutes).
To our knowledge this is the first study applying a load sharing paradigm in controlling multiple muscles with surface FES in a multidimensional biomechanical system. The load sharing approach controlled the thumb towards the target forces in the 2D control experiments. With feedforward force control only, errors were larger in stroke subjects, compared to healthy subjects. However, with added feedback control, significant differences in control performance had disappeared. Therefore the methodology for multi-dimensional feedback force control presented here has potential applicability as part of post stroke rehabilitation techniques. Especially when applied earlier after stroke and muscles are stronger.
This study was part of the Interreg IV MIAS-ATD project, part of the European regional development fund. We would like to acknowledge the support of all project partners, which is greatly appreciated.
- Ovbiagele B, Nguyen-Huynh M: Stroke epidemiology: advancing our understanding of disease mechanism and therapy. Neurotherapeutics 2011, 8: 319-329. .[10.1007/s13311-011-0053-1] [http://0-dx.doi.org.brum.beds.ac.uk/10.1007/s13311-011-0053-1] .[10.1007/s13311-011-0053-1] 10.1007/s13311-011-0053-1PubMed CentralView ArticlePubMedGoogle Scholar
- Crago PE, Nakai RJ, Chizeck H: Feedback regulation of hand grasp opening and contact force during stimulation of paralyzed muscle. IEEE Trans Biomed Eng 1991, 38: 17-28. 10.1109/10.68205View ArticlePubMedGoogle Scholar
- Popović D, Sinkjær T, Popović M: Electrical stimulation as a means for achieving recovery of function in stroke patients. NeuroRehabilitation 2009, 25: 45-58.PubMedGoogle Scholar
- Micera S, Keller T, Lawrence M, Morari M, Popovic D: Wearable neural prostheses. Eng Med Biol Mag 2010,29(3):64-69.View ArticleGoogle Scholar
- Wolbrecht ET, Chan V, Reinkensmeyer DJ, Bobrow JE: Optimizing compliant, model-based robotic assistance to promote neurorehabilitation. IEEE Trans Neural Syst Rehabil Eng 2008,16(3):286-297.View ArticlePubMedGoogle Scholar
- Westerveld AJ, Schouten AC, Veltink PH, van der Kooij H: Selectivity and resolution of surface electrical stimulation for grasp and release. IEEE Trans Neural Syst Rehabil Eng 2012, 20: 94-101.View ArticlePubMedGoogle Scholar
- Kaufman K, An K, Litchy W, Cooney W, Chao E, et al.: In-vivo function of the thumb muscles. Clin Biomech 1999,14(2):141-150. 10.1016/S0268-0033(98)00058-8View ArticleGoogle Scholar
- Pearlman J, Roach S, Valero-Cuevas F: The fundamental thumb-tip force vectors produced by the muscles of the thumb. J Orthop Res 2004,22(2):306-312. 10.1016/j.orthres.2003.08.001View ArticlePubMedGoogle Scholar
- Happee R: Inverse dynamic optimization including muscular dynamics, a new simulation method applied to goal directed movements. J Biomech 1994,27(7):953-960. 10.1016/0021-9290(94)90267-4View ArticlePubMedGoogle Scholar
- Prilutsky BI, Zatsiorsky VM: Optimization-based models of muscle coordination. Exerc Sport Sci Rev 2002, 30: 32. 10.1097/00003677-200201000-00007PubMed CentralView ArticlePubMedGoogle Scholar
- Lujan J, Crago P: Automated optimal coordination of multiple-DOF neuromuscular actions in feedforward neuroprostheses. Biomed Eng IEEE Trans 2009, 56: 179-187.View ArticleGoogle Scholar
- Kuhn A, Keller T, Lawrence M, Morari M: The influence of electrode size on selectivity and comfort in transcutaneous electrical stimulation of the forearm. Neural Syst Rehabil Eng IEEE Trans 2010,18(3):255-262.View ArticleGoogle Scholar
- Seber GAF, Wild CJ: Nonlinear Regression. Hoboken, NJ: Wiley-Interscience; 2003.Google Scholar
- Ziegler JG, Nichols N: Optimum settings for automatic controllers. Trans ASME 1942, 64: 759-768.Google Scholar
- Ferreau H, Bock H, Diehl M: An online active set strategy to overcome the limitations of explicit MPC. Int J Robust Nonlinear Control 2008,18(8):816-830. 10.1002/rnc.1251View ArticleGoogle Scholar
- Flanagan JR, Burstedt MK, Johansson RS: Control of fingertip forces in multidigit manipulation. J Neurophysiol 1999,81(4):1706-1717.PubMedGoogle Scholar
- Singh G, Boddu S, Chakravorty I, Bairy G, Ganesh M: An instrumented glove for monitoring forces during object manipulation. In Point-of-Care Healthc Technol (PHT) IEEE. Bangalore, India: IEEE; 2013:212-215.View ArticleGoogle Scholar
- Westerveld AJ, Kuck A, Schouten A, Veltink P, van der Kooij H: Grasp and release with surface functional electrical stimulation using a Model Predictive Control approach. In Engineering in Medicine and Biology Society (EMBC), 2012 Annual International Conference of the IEEE. San Diego, California, USA: IEEE; 2012:333-336.View ArticleGoogle Scholar
- Kozin SH: The anatomy of the recurrent branch of the median nerve. J Hand Surg 1998,23(5):852-858. [http://0-www.sciencedirect.com.brum.beds.ac.uk/science/article/pii/S0363502398801627]  10.1016/S0363-5023(98)80162-7View ArticleGoogle Scholar
- Schearer E, Liao Y, Perreault E, Tresch M, KM L: Optimal sampling of recruitment curves for functional electrical stimulation control. In 34th Annual International Conference of the IEEE EMBS. San Diego, California, USA: IEEE; 2012:329-332.Google Scholar
- Schearer E, Liao YW, Perreault E, Tresch M, Memberg W, Kirsch R, Lynch K: System identification for 3D force control of a human arm neuroprosthesis using functional electrical stimulation. In Robotics and Automation (ICRA), 2012 IEEE International Conference on. St. Paul, Minnesota, USA: IEEE; 2012:3698-3705.View ArticleGoogle Scholar
- Malesevic N, Maneski LZ, Ilic V, Jorgovanovic N, Bijelic G, Keller T, Popovic D: A multi-pad electrode based functional electrical stimulation system for restoration of grasp. J Neuro Eng Rehabil 2012, 9: 66. [http://www.jneuroengrehab.com/content/9/1/66]  10.1186/1743-0003-9-66View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.