Impact of Functional Electrical Stimulation of Lower Limbs during Sitting Pivot Transfer Motion for Paraplegic People

Abstract

Individuals with Spinal Cord Injury (SCI), perform Sitting Pivot Transfer (SPT) motion around fifteen times a day using upper extremities. It can lead to upper limbs pain and often shoulder complications. In this paper, we investigate the influence of Functional Electrical Stimulation (FES) on SPT motion of a paraplegic person. First, we proposed to develop a dynamic optimization method in order to predict SPT motion of an able-bodied subject. This approach has been validated by comparing the computed SPT trajectories with the ones measured during the experiment with an able-bodied subject. Then, we used the optimization tool to analyze the influence of FES on the SPT maneuver of paraplegic persons. Our results suggest that FES can decrease arm participation during the transfer motion of a paraplegic person.

Keywords

Multi-contact Motion optimization FES Paraplegia Sitting Pivot Transfer

1. Introduction

Among all transfer activities, SPT is reported to be one of the most commonly used by individuals with SCI Bromley (1998). The performance of SPT may rank among the most demanding functional mobility activities for the upper limbs Gagnon, Nadeau, Desjardins & Noreau (2008). Transferring from a wheelchair to a treatment table, bed, tub/shower bench, toilet seat, car seat and vice versa represents typical SPT performed by individuals with SCI. When initiating SPT, paraplegic individuals move their wheelchair as close as possible to the target seat. They move the buttocks forward close to the front edge of the wheelchair seat and, with the help of their arms, place firmly their feet on the floor. Then, they place one hand, called the trailing hand, on a stable position on the wheelchair and the other hand, known as the leading hand, on the target surface far enough to leave sufficient space for the buttocks. From this starting position, subjects bend their trunk forward and sideways, meanwhile they lift up their body and sustain their weight with the arms. After that, with a very rapid twisting motion they place the buttocks on the target seat. The transfer is concluded when the subjects reach again a seated postural stability. The performance of SPT varies from one person to the other depending on the subject's features, like lower limb spasticity, and environmental factors, like wheelchair design Gagnon, Nadeau, Noreau, Eng & Gravel (2008); Koontz, Kankipati, Lin, Cooper & Boninger (2011). SCI persons perform around fifteen SPT motions a day Gagnon, Nadeau, Noreau, Dehail & Gravel (2008); Gagnon, Nadeau, Noreau, Eng & Gravel (2008).

Usually SPT is divided into three distinct phases:

Pre-lift is the preparatory phase that ends when the trunk bends forward,

Lift is the phase when most of the body weight is supported by the upper limbs,

Post-lift is the rebalancing phase, when the subject is sitting on the new seat.

The lift phase lasts forty percent of the entire transfer duration Gagnon, Nadeau, Noreau, Eng& Gravel (2008). A lot of scientific works have been focused on experimental studies in which the kinetics and the kinematics of the SPT motion has been analyzed Gagnon, Nadeau, Desjardins &Noreau (2008); Gagnon, Nadeau, Noreau, Dehail & Gravel (2008); Koontz, Kankipati, Tsai & Lin (2011); Koontz, Kankipati, Lin, Cooper & Boninger (2011); Koontz, Lin, Kankipati, Boninger & Cooper (2011); Tanimoto et al. (2008); van Drongelen et al. (2011). Though Functional Electrical Stimulation enables restoration of the motion in paralyzed limbs Popovic & Sinkjær (2000), to the best of our knowledge, the influence of FES on SPT motion has not been investigated yet. It has been documented that paraplegic patients have upper limbs pain and often shoulder complications Koontz, Kankipati, Lin, Cooper & Boninger (2011); van Drongelen et al. (2011). Hence, the aim of the present study is to enhance the scientific research concerning paraplegics' SPT with assistance of functional electrical stimulation using dynamic optimization approach.

In motor control and biomechanics literature, a classical assumption is that human beings perform a motion according to certain optimal criteria, i.e. motion control can be related to a problem of minimizing a biomechanical cost function. Optimization processes usually provide a better understanding of human postural and locomotor system Flash & Hogan (1985); Kuzelicki et al. (2005); Lin et al. (1999); Martin et al. (2006); Pandy et al. (1995); Rosenbaum et al. (1995); Uno et al. (1989); Yamasaki et al. (2011). In human motion science domain, the proposed optimization algorithms generally minimize the time integral of the square of a quantitative function (jerk, acceleration, torque, torque-change). For example, Flash et al. suggested that the arm reaching motion minimizes the time integral of the hand position jerk Flash & Hogan (1985). Rosenbaum et al. computed motion with the joint angle jerk as the cost function Rosenbaum et al. (1995). Uno et al improved the optimization criterion by minimizing the torque-change, i.e. minimizing the sum of the joint torque derivatives Uno et al. (1989). Optimal sit-to-stand motion for able-bodied subjects was investigated by Pandy et al. Pandy et al. (1995). Authors used a three segment 2-D model which includes eight muscle groups. Optimal neural excitation signals were computed by minimizing the muscle forces and their derivatives. Kuzelicki et al. guarantee symmetric motion using a 3-D eleven-degrees of freedom-model by minimizing the sum of torque, torque-change and the difference between left and right ground reaction forces Kuzelicki et al. (2005). Yamasaki et al. suggested that a cost function which minimizes torque change could better describe rising motion than a cost function which minimizes joint jerk Yamasaki et al. (2011). As far as we know, there is no study which deals with optimal transfer trajectories in paraplegic population.

The goal of this work is to investigate, first the ability of optimization process to predict SPT trajectories in able-bodied subjects and then, the influence of FES on arm efforts during SPT motion of a paraplegic person. Hence, SPT motion of one able-bodied subject were measured. The trajectory of Center of Mass (CoM) position estimated using experimental data were compared with the ones obtained using optimization process. Once the results were validated, the influence of FES stimulation on the corresponding motion in paraplegic subjects was investigated.

2. Method

2.1 Kinematic model

In this study, we used the biomechanical model of a human shown on Figure 1. The model consists in nineteen segments and eighteen spherical joints, i.e. fifty four Degrees of Freedom (DoF), but we set a constant value to the six DoF of the neck. Despite the existence of accurate models for the knees and shoulders Gini et al. (2007); Murray & Larochelle (1998), we consider that spherical joints are good enough to describe the whole body motion. The spherical joints are considered as three rotations around X, Z and Y axis. The joint limits are the one presented in Table 1. The kinematic and dynamic properties of this model are scaled regarding the desired height and weight using the anthropometric Winter tables Winter (2009).

Table 1.

Joints limits for the biomechanical model (in radian).

Joint	θ _x	θ _y	θ _z
Pelvis	−0.87 : 0.61	−0.70 : 0.70	−0.09 : 0.09
Mid-spine	−0.61 : 0.47	−0.34 : 0.34	−0.63 : 0.63
Lower neck	−1.13 : 0.70	−0.61 : 0.61	−0.61 : 0.61
Upper neck	−1.13 : 0.70	−0.61 : 0.61	−0.61 : 0.61
left clacivula¹	−0.01 : 0.01	−0.15 : 0.15	−0.15 : 0.15
left shoulder¹	−2.26 : 3.14	−2.26 : 0.00	−1.04 : 0.52
left elbow¹	−0.001 : 2.79	−1.57 : 1.57	−0.01 : 0.01
left wrist¹	−1.22 : 0.87	−0.01 : 0.01	−0.52 : 0.35
left hip¹	−0.87 : 1.65	−0.35 : 1.13	−0.61 : 0.61
left knee¹	−2.26 : 0.001	−0.01 : 0.01	−0.01 : 0.01
left ankle¹	−0.69 : 0.52	−0.34 : 0.34	−0.34 : 0.34

2.2 Dynamic modeling and balance

To study the contact forces of the hands, we start from the the dynamic model as presented:

Figure 1.

Representation of the kinematic model and joint frames for zero joint positions.

[\begin{array}{l} Γ \\ O \end{array}] = [\begin{array}{l} D_{1} (q, \dot{q}, \overset{..}{q}) \\ D_{2} (q, \dot{q}, \overset{..}{q}) \end{array}] + [\begin{array}{l} J_{1}^{T} (q) \\ J_{2}^{T} (q) \end{array}] F

(1)

Where q ∈ ℝⁿ is a vector containing the joint positions, Γ ∈ ℝⁿ the vector of the joint torque D₁ ∈Rⁿ and D₂ ∈Rⁿ the dynamics effects (sum of inertia, Coriolis, centrifugal and gravity) due to the joint trajectories, J₁ ∈R^n×3Nf and J₂ ∈R^6×3Nf the components of the Jacobian matrix and F={F_1,F₂,…} the vector of the contact forces. Note that the generalized formulation of the inverse dynamic model usually contains the position and orientation of the waist that are computed from the joint angles and the position of the right foot knowing its constant position and orientation.

2.3 Computation of the contact forces

Equation (1) emphasizes the link between the joint trajectories q(t) , the contact forces F and the joint torquesr. From the joint trajectories, there is an infinity of couple forces-torques. Here, we present how to find a set of contact forces that compensates the dynamic effects, ensures desired torques and encourages the balance as much as possible.

Due to the non-planar contact points, we cannot use classical method such as the Zero Moment Point Vukobratovic & Borovac (2004) to characterize balance. In order to ensure the balance and desired torque values, the contact forces must counterpart the dynamic effects and fulfill the friction and unilateral constraint to avoid undesired sliding or taking off, such as:

D_{2} + J_{2}^{T} F = 0

(2)

\forall e [Γ_{e}] = D_{1, e} + J_{1, e}^{T} F

(3)

\forall i \in {1, ..., N_{f}} {\begin{cases} F_{i}^{n} > 0 \\ F_{i}^{t^{2}} \leq μ_{i}^{2} F_{i}^{n^{2}} \end{cases}}

(4)

with F_iⁿ and F_i^t are the normal and tangential components of the contact forces F_i, μ_i the friction coefficient and D_ie and J₁^T,_e are the e^th line of D₁ and J₁ ^T. Equations (2) and (4) refer to the balance criteria presented in Harada et al. (2007) which states that the contact wrench sum must remain in the contact wrench cone to ensure balance.

To solve Equations (2) and (3), one can use the Moore-Penrose pseudo-inverse matrix of the Jacobian matrix.

However, the pseudo-inverse matrix minimizes the instantaneous norm of the contact forces, without any effect on the friction or sliding constraint. We rather consider the contact forces that are as close as possible to the normal direction to the contact surface, i.e. that are the solution to the following problem:

\begin{array}{l} \min \frac{1}{2} \sum_{i} β_{i} (α_{i} | | F_{i}^{t} | |^{2} + F_{i}^{n^{2}}) \\ Σ_{i} ([\begin{array}{c} {\overset{\land}{P}}_{i} A_{i} \\ A_{i} \end{array}] [F_{i}]) + [D_{2}] = 0 \\ Σ_{i} (η_{e, i} [F_{i}]) + D_{1, e} - [Γ_{e}] = 0 \end{array}

(5)

Where P̂_iA_i∈R^3×3 and η_e,i∈R³ appear in the decomposition of the Jacobian matrix J^T₂ and J^T₁. β_i is a weight value to modify the repartition of the different contact forces and α_i gives more importance to the tangential components regarding to the normal one for each contact forces. The solution to problem (5) counterparts the dynamics effects and produces the desired joint torque Γ_e while fulfilling as best as possible the friction and unilateral constraint (by setting a_i >> 1). The solution of problem (5) is:

F_{i} = W_{i}^{- 1} [{\hat{P}}_{i} A_{i} A_{i} η_{e, i}] Ω^{- 1} [\begin{array}{l} D_{2} \\ D_{1, e} - [Γ_{e}] \end{array}]

(6)

with:

Ω = \underset{i}{Σ} ([\begin{array}{c} {\hat{P}}_{i} A_{i} \\ A_{i} \\ η_{e, i} \end{array}] W_{i}^{- 1} [{\hat{P}}_{i} A_{i} A_{i} η_{e, i}])

(7)

with W_i = diag(β_i, α_i, β_i, α_ii, β_i) (we assume that the z-axis is the normal direction of the contact forces) and $Ω \in R^{(6 + N_{e}) \times (6 + N_{e})}$ is a are matrix that can be easil inverted, by using, for instance, the Gauss-Jordan algorithm. From the computation of the contact forces, we can compute the joints torques using Equation (1).

2.4 Optimization process

In mathematics or computational science, optimization refers to the selection of the best element from a set of available alternatives. In our case, the optimization problem consists in minimizing a real function while ensuring some constraints, such as respecting joint limits, ensuring balance during the motion, etc. As mentioned in Introduction section, this approach was already used in the case of motion generation in the robotics field and field of human motion science Lee et al. (2005); Lengagne, Kheddar & Yoshida (2011); Mombaur et al. (2005).

In this paper, we refer to the optimization process in order to generate SPT motion for the biomechanical model as presented in Section 2. This optimal motion, described through a parameter vector X, minimizes a cost function C and ensures a set of discrete z_k, and continuous inequality g_i and equality h_j constraints. Thus the optimal motion is the solution of:

\begin{array}{c} a r g m i n C (X) \\ \forall i, \forall t \in [0, T_{f}] g_{i} (X, t) \leq 0 \\ \forall j, \forall t \in [0, T_{f}] h_{j} (X, t) = 0 \\ \forall k z_{k} (X, t_{k}) \leq 0 \end{array}

(8)

The parameter vector $X \in R^{n} = {p_{1, 1}, \dots, p_{1, m}, \dots, β, T_{f}}$ contains the parameters of the joint trajectories p_j,i, β is a set of parameters as presented in Section 2.3 and T_f is the motion duration. As already used in Lee et al. (2005), we represent the joint trajectories q (t) using B-Spline basis functions De Boor (1978):

q_{j} (t) = \sum_{i = 1}^{m} b_{i}^{K} (t) p_{j, i}

(9)

As presented in Lengagne, Kheddar & Yoshida (2011), we used a time-interval discretization based on a polynomial approximation to consider a finite set of constraints that properly represents the continuous constraints. We used the optimization software IPOPT Wächter & Biegler (2006) that is a free and effective solution with C++ interface. We consider cost function C as the weighted sum of joint torques and joint jerks, as presented:

C (q) = a \int_{0}^{T} \sum_{i} Γ_{i}^{2} d t + b \int_{0}^{T} \sum_{i} \overset{... 2}{q i} d t

(10)

with a=1e-2and b=1e-5. The values a and b were set heuristically as explained in Lengagne, Kheddar, Druon & Yoshida (2011).

2.5 Experimental Validation

As mentioned above, the SPT can be divided into three phases: pre-lift, lift and post-lift phases. For each phase an initial and final position can be defined:

Pre-lift phase: the subject starts sitting on the initial chair with the feet on the floor, the right hand on the handle on the right side of the initial chair and left hand on the handle on the left side of the target chair;

Lift phase: the buttocks are moved from the initial to the target chair, while the positions of the feet and hands do not change;

Post-lift phase: the subject regains a sitting posture on the target chair, position of the feet and hands unchanged.

One able-bodied subject (male, age 29, height 1.75m, weight 72kg) participated in this study after signing an informed consent. The angle between the two chairs was 20 degrees. Kinematic variables were measured using a stereophotogrammetric system (8 Mx cameras, VICON©) at 100Hz sampling rate. In order to estimate accurately small postural modifications, 35 reflective markers have been located on the subject's body using Plug-In-Gait template http://www.vicon.com/ (n.d.). The experimental set-up is presented in Figure 2.

Figure 2.

Experimental setup and initial position of the subject.

The subject was instructed to assume comfortable feet placement, to place both hands on the handles and to keep the same position of the feet and hands during the experiment. The initial position of the subject is shown in Figure 2. The subject performed the SPT at his own speed and way five times.

The Vicon Nexus© and LIFEMOD© softwares have been used to estimate CoM positions from experimental data. Each trial was time-normalized to 100 samples. Root Mean Square error (RMS) and Pearson's correlation coefficients between experimental and optimized CoM positions have been calculated in order to compare those variables.

3. Results

3.1 The ability of optimization process to predict SPT trajectories

Optimal SPT trajectories were calculated using described cost function. In Fig. 3, trajectories of CoM position in Anterior-Posterior (A/P), vertical, and Medio-Lateral (M/L) directions are presented. It can be observed that differences between the optimized CoM positions and the ones calculated from measured data in anterior-posterior and vertical direction are not bigger than the variability among the different trials of the same subject. The difference between the optimized CoM position in M/L direction and the corresponding ones estimated from experimental data is slightly bigger than in the two other directions.

Figure 3.

Computed (red line) CoM position in A/P, vertical and M/L direction in able-bodied subject. The mean CoM trajectory, in A/P, vertical and M/L direction, calculated from experimental data (black line) and gray lines marking the plus/minus standard deviation interval are shown.

Table 2 presents Root Mean Square error and Person's correlation coefficients calculated between CoM position computed using our optimization process and CoM position calculated using experimental data for all five trials. The RMS error is smaller than 10cm in A/P direction (Trial 5) and M/L direction (Trial 1). Observing Person's correlation coefficients, it can be noticed that our optimization process well reproduces the behavior of our subject in M/L and vertical direction.

Table 2.

RMS error and Person's correlation coefficients calculated between CoM position computed using optimization process and CoM position calculated using experimental data.

	Number of Trial	1	2	3	4	5
CoM A/P	RMS Correlation Coefficients	0.1 0.6732	0.0873 0.2951	0.0741 0.1069	0.0507 0.7178	0.081 0.6402
CoM Vertical	RMS Correlation Coefficients	0.0581 0.6552	0.0511 0.6732	0.0409 0.6706	0.0.056 0.7395	0.0452 0.6908
CoM M/L	RMS Correlation Coefficients	0.0863 0.9017	0.0716 0.9722	0.0481 0.9939	0.0.0687 0.9191	0.1078 0.8643

3.2 Influence of FES assistance on hand forces during SPT motion

The optimal STP trajectories are calculated for nine different scenarios (Sc). The first scenario, Sc 1, represents the behavior of an able bodied subject, allowing the variations of knee joint torque. In the other scenarios, the value of the knee joint torque is constant during the motion. In the second scenario, Sc 2, the knee torque is 0 Nm and represents the behavior of a paraplegic subject performing the motion without FES assistance. Other scenarios, from the third (Sc 3) to the ninth (Sc 9), represent a paraplegic subject performing FES-assisted SPT task while the stimulation parameters were increased. The torque values are following. Sc 3: 5 Nm, Sc 4: 10 Nm, Sc 5: 15 Nm, Sc 6: 20 Nm, Sc 7: 30 Nm, Sc 8: 40 Nm and Sc 9: 50 Nm. The virtually stimulated muscles are quadriceps and biceps femoris.

The following assumptions have been made:

for Sc 2 -Sc 9, the voluntary control of the knee joints is not possible,

the virtually stimulated bi-articular muscles produce torque control only at the knee joints,

knee joint torques do not change during the lift phase of the transfer (it has zero value during pre and post lift phase) and no stimulation leads to a null knee joint torque.

After validating our approach as presented previously, we computed the contact forces on the hands using the biomechanical model presented in Section 2, and the cost function descibed in Section 2.4. The results of the computation for the nine scenarios are presented in Figure 4. As expected, the contact forces under the hands are lower in the case of an able-bodied person (Sc 1) than in the case of a paraplegic person performing the motion without assistance of FES (Sc 2). In addition, analyzing other scenarios, it appears that the stimulation of the knee joints has an impact during the lift phase during which the contact forces of the hands decrease when the equivalent torque produced by the FES of lower limbs muscles increases.

Figure 4.

Representation of the contact forces for the several cases.

4. Conclusion

The results of this study indicate that it is possible to describe sitting pivot transfer of an able-bodied subject within the frame of optimization theory. Taking into account the body dynamics and kinetics, the proposed method was able to reproduce motion with the CoM trajectories quite closed to the one measured during experiments or estimated from the optimization process. Using it, we have shown that functional electrical stimulation applied on paralyzed lower limbs has an impact on arm efforts during the SPT motion. Our study indicates that using FES decreases arm efforts during the SPT motion. The ability to achieve those transfers with minimal participation of the upper limbs would help to preserve long-term shoulder integrity and improve the life of paraplegic patients.

Besides modeling of dynamics and kinematics of human body, Future studies will include modelling of muscle properties and finding the optimal FES pattern in order to minimize the fatigue of stimulated lower limbs muscles during SPT motion.

Footnotes

5. Acknowledgments

The authors thank the subject for participating to this study and making it possible. This work is partially supported by grants from the German Research Foundation (DFG: Deutsche Forschungsgemeinschaft).

1

for the right part of the model the values of y joint limits are opposite.

References

Bromley

(1998). Tetraplegia and paraplegia: a guide for physiotherapists, 5th ededn, Churchill Livingstone, Edinburgh.

De Boor

(1978). A Pratical Guide to Splines, Vol. 27, Springer-Verlag, New York.

Flash

& Hogan

(1985). The coordination of arm movements: An experimentally confirmed mathematical model. J Neurosci 5 (7): 1688–703.

Gagnon

Nadeau

Desjardins

& Noreau

(2008). Biomechanical assessment of sitting pivot transfer tasks using a newly developed instrumented transfer system among long-term wheelchair users, J Biomech 41(5): 1104–10.

Gagnon

Nadeau

Noreau

Dehail

& Gravel

(2008). Quantification of reaction forces during sitting pivot transfers performed by individuals with spinal cord injury, J Rehabil Med 40(6): 468–76.

Gagnon

Nadeau

Noreau

Eng

J. J.

& Gravel

(2008). Trunk and upper extremity kinematics during sitting pivot transfers performed by individuals with spinal cord injury, ClinBiomech (Bristol, Avon) 23(3): 279–90.

Gini

Scarfogliero

& Folghraiter

(2007). Human-oriented biped robot design:

Insights into the development of a truly antropomorphic leg, IEEE International Conference on Robotics and Automation, pp. 2910–2915.

Harada

Hattori

Hirukawa

Morisawa

Kajita

& Yoshida

(2007). Motion planning for walking pattern generation of humanoid, Intelligent Robots and Systems, 2007. IROS 2007. IEEE/RSJ International Conference on, pp. 4227–4233. http://www.vicon.com/ (n.d.).

10.

Koontz

Kankipati

Tsai

C.-Y.

& Lin

Y.-S.

(2011). Biomechanical analysis of sitting pivot wheelchair transfers, University of Southern California, Long Beach, CA.

11.

Koontz

A. M.

Kankipati

Lin

Y.-S.

Cooper

R. A.

& Boninger

M. L.

(2011). Upper limb kinetic analysis of three sitting pivot wheelchair transfer techniques, ClinBiomech (Bristol, Avon) 26 (9): 923–9.

12.

Koontz

A. M.

Lin

Y.-S.

Kankipati

Boninger

M. L.

& Cooper

R. A.

(2011). Development of custom measurement system for biomechanical evaluation of independent wheelchair transfers, J Rehabil Res Dev 48 (8): 1015–28.

13.

Kuzelicki

Zefran

Burger

& Bajd

(2005). Synthesis of standing-up trajectories using dynamic optimization. Gait Posture 21 (1): 1–11.

14.

Lee

S.-H.

Kim

Park

Kim

& Bobrow

J. E.

(2005). Newton-type algorithms for dynamics-based robot movement optimization, IEEE Transactions on robotics, Vol. 21, pp. 657–667.

15.

Lengagne

Kheddar

Druon

& Yoshida

(2011). Emulating human leg impairments and disabilities in walking with humanoid robots, IEEE Int. Conf. on Robotics&BIOmimetics.

16.

Lengagne

Kheddar

& Yoshida

(2011). Considering floatting contact and un-modeled effects for multi-contact motion generation, Workshop on Humanoid Service Robot Navigation in Crowded and Dynamic Environments at the upcoming IEEE Humanoids Conference.

17.

Lin

C. J.

Ayoub

M. M.

& Bernard

T. M.

(1999). Computer motion simulation for sagittal plane lifting activities, International Journal of Industrial Ergonomics 24 (2): 141–155.

18.

Martin

Cahouët

Ferry

& Fouque

(2006). Optimization model predictions for postural coordination modes. J Biomech 39 (1): 170–6.

19.

Mombaur

K. D.

Bock

H. G.

Schloder

J. P.

& Longman

R. W.

(2005). Self-stabilizing somersaults, Robotics, IEEE Transactions on [see also Robotics and Automation, IEEE Transactions on] 21: 1148–1157.

20.

Murray

& Larochelle

(1998). A classification scheme for planar 4r, spherical 4 r, and spatial rccc linkages to facilitate computer animation, ASME Design engineering Technical Conferences.

21.

Pandy

M. G.

Garner

B. A.

& Anderson

F. C.

(1995). Optimal control of non-ballistic muscular movements: A constraint-based performance criterion for rising from a chair. J BiomechEng 117 (1): 15–26.

22.

Popovic

& Sinkjær

(2000). Control of movement for the physically disabled, Springer-Verlag London.

23.

Rosenbaum

D. A.

Loukopoulos

L. D.

Meulenbroek

R. G., J., V.

& Engelbrecht

S. E.

(1995). Planning reaches by evaluating stored postures. Psychol Rev 102 (1): 28–67.

24.

Tanimoto

Nanba

Tokuhiro

Ukida

& Yamamoto

(2008). Measurement system of transfer motion for patients with spinal cord injuries. IEEE Transaction on instrumentation and measurement. 57 (1): 213–219.

25.

Uno

Kawato

& Suzuki

(1989). Formation and control of optimal trajectory in human multijoint arm movement. minimum torque-change model. BiolCybern 61 (2): 89–101.

26.

vanDrongelen

van der Woude

L. H. V.

& Veeger

H. E. J.

(2011). Load on the shoulder complex during wheelchair propulsion and weight relief lifting, ClinBiomech(Bristol, Avon) 26 (5): 452–7.

27.

Vukobratovic

& Borovac

(2004). Zero-moment point: Thirty five years of its life, International Journal of Humanoid Robotics 1 (1): 157–173.

28.

Wächter

& Biegler

L. T.

(2006). On the implementation of a primal-dual interior point filter line search algorithm for large-scale nonlinear programming, Mathematical Programming 106: 22–57.

29.

Winter

D. A.

(2009). Anthropometry, in Biomechanics and Motor Control of Human Movement, Fourth Edition, John Wiley & Sons, Inc.

30.

Yamasaki

H. R.

Kambara

& Koike

(2011). Dynamic optimization of the sit-to-stand movement. J ApplBiomech 27 (4): 306–13.