Computational and Mathematical Methods in Medicine

Volume 2013 (2013), Article ID 570878, 7 pages

http://dx.doi.org/10.1155/2013/570878

Research Article

Theoretical Hill-Type Muscle and Stability: Numerical Model and Application

S. Schmitt,1,2 M. Günther,1,3 T. Rupp,1,2 A. Bayer,¹ and D. Häufle¹

¹Department of Sports and Exercise Science, University of Stuttgart, Allmandring 28, 70569 Stuttgart, Germany

²Stuttgart Research Centre for Simulation Technology, University of Stuttgart, Pfaffenwaldring 5a, 70569 Stuttgart, Germany

³Institute of Sports Science, Science of Motion, University of Jena, Seidelstraß 20, 07749 Jena, Germany

Received 19 July 2013; Accepted 19 September 2013

Academic Editor: Zhonghua Sun

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

Abstract.

The construction of artificial muscles is one of the most challenging developments in today’s biomedical science. The application of artificial muscles is focused both on the construction of orthotics and prosthetics for rehabilitation and prevention purposes and on building humanoid walking machines for robotics research. Research in biomechanics tries to explain the functioning and design of real biological muscles and therefore lays the fundament for the development of functional artificial muscles. Recently, the hyperbolic Hill-type force-velocity relation was derived from simple mechanical components. In this contribution, this theoretical yet biomechanical model is transferred to a numerical model and applied for presenting a proof-of-concept of a functional artificial muscle. Additionally, this validated theoretical model is used to determine force-velocity relations of different animal species that are based on the literature data from biological experiments. Moreover, it is shown that an antagonistic muscle actuator can help in stabilising a single inverted pendulum model in favour of a control approach using a linear torque generator.

1. Introduction

Research in muscle biomechanics, a vital and broad field for over 80 years now (A.V. Hill 1922: Nobel prize in physiology and medicine for his discovery relating to the production of heat in the muscle), explains the function and design of real biological muscles and therefore lays the fundament for the development of functional artificial muscles. Nevertheless, structure and functioning of biological muscles are not (yet) fully understood.

In biology, microscopic muscle models are able to predict muscle characteristics and functioning of biological muscles quite well [1–9]. Unfortunately and as a tradeoff, they require a large number of parameters. In a bionics approach it is an enormous challenge to realise all these properties of biological muscle in one artificial muscle at once [10].

Macroscopic muscle models are commonly based on phenomenology. Macroscopic muscle models are indeed of (limited) predictive character but do not incorporate any structural knowledge. Recently, the nonlinear (hyperbolic-like) Hill-type force-velocity relation was derived from simple mechanical components [11]. It was shown that a contractile element (CE) consisting of a mechanical energy source (active element AE), a parallel damper element (PDE), and a serial element (SE) exhibits operating points with nonlinear (hyperbola-like) force-velocity dependency. In this concept, the force-velocity relation is no longer a phenomenological outcome of a black box (i.e., the CE) but rather a physical outcome of the interaction of the three elements AE, PDE, and SE. Based on this concept, it is now possible to describe in detail which structural arrangement is necessary to get a biology-like force-velocity relation on a macroscopic level. Therefore, this concept can be interpreted as a basic engineering design for the CE of a Hill-type artificial actuator [12–15]. In this paper, the meaning of the structural arrangement of the simple mechanical components already published will be revisited. Furthermore, by taking one first example of technical embodiment, it will be shown how this concept can help to construct more biologically-motivated artificial muscles. Most importantly, it can be shown that the application of a Hill-type muscle model could improve biomechanical stability. An antagonistic pair of our muscle model theoretically shows a first demonstration of how an artificial muscle could help in the stabilisation of a technical machine. As a result, the control of an inverted pendulum can be improved by the use of a muscle-like drive in favour of a linear torque generator.

2. Material and Methods

2.1. Hill’s Original Formulation of Muscle Dynamics

In his famous paper “The heat of shortening and the dynamic constants of muscle” [16], Hill firstly formulated the simple and convenient equation describing the muscle’s contraction dynamics:

where the symbol “” denotes the current muscle force, “” the respective contraction velocity of the muscle, “” the muscle’s maximum isometric force, and “” and “” Hill’s so-called dynamic constants of muscle, which we call the Hill parameters.

2.2. Derivation of the Hill Parameters

In a recent paper [11] it was demonstrated that the phenomenologically found hyperbolic force-velocity relation of a concentrically contracting assembly of activated muscle fibres [16] can be derived from the simple mechanical arrangement (Figure 1(a)) of an arbitrary force generating (active) element (AE) to which a damper (PDE) is connected in parallel and a serial element (SE) in series fulfilling the force equilibrium:

where the symbol “” denotes a force produced by the element denoted by a corresponding index and the kinematic relation for the lengths (symbols “”) of the elements AE, PDE, and SE:

with representing the contractile element length. Note that a dot symbol “” denotes the first time derivative of a length , that is, an element’s contraction velocity. Please refer to the appendix or [13, 15] for a more detailed description.

(a)

(b)

Figure 1: (a) Theoretical construct of the CE [11]. The CE consists of three elements: active element AE, parallel damping element PDE, and serial element SE. is the origin of the CE, the length of the AE/PDE, and the length of the whole CE. By choosing in theory, we can turn the SE off in order to represent a contractile element without any compliance. (b) Hardware design. AE and PDE were realised with an electric motor, SE with a mechanical spring. A second electric motor was used to exert a defined external force on the CE. The electromagnet held the muscle at constant length until release at time .

2.3. Technical Embodiment

The hardware implementation (Figure 1(b)) was done analogously. Both AE and PDE were realised each with an electric motor (Maxon ECmax40) [14]. The motor torque () was controlled by Maxon digital EC-motor control units (DEC 70/10). Both motors were mounted from opposite sides to the same disc with the radius . The disc was used to coil up a steel rope and exert a force

on the rope. The force characteristics of the PDE and AE and were implemented in MATLAB Simulink through Real-Time Workshop and Real-Time Windows Target. Thus, the motors could exert the specified force on the steel rope as required by the theoretical construct. For the SE, a spring () was tied to a steel rope. Another motor could exert a defined external force on the CE construct. As sensor signals, the motor shaft positions were recorded by optical encoders (Scancon 2RMHF 5000 pulses/revolution), representing the internal degree of freedom and the total CE length . A load cell (Transducer Techniques MLP 25 with amplifier TM0-1-24) was used to calibrate motor torques and exerted forces. All sensor data were recorded with MATLAB Simulink via a Sensoray 626 AD I/O at 1 kHz.

To investigate the force-velocity characteristics of the artificial CE, two types of experiments had to be performed. The first experiment was an isometric contraction (contraction with constant CE length: .). For this purpose, the CE end was fixed to the electromagnet guaranteeing a constant CE length. Then the AE activation was set to (maximum activation) and the shortening of the AE (rotation of the motors) was recorded. The time from the beginning of the activation until the end of AE shortening and the maximum isometric force were evaluated.

Isotonic quick release experiments were performed to guarantee a defined (velocity ratio; please see appendix for a definition of ). Each isotonic quick release contraction experiment started like an isometric contraction, only that the CE was released at ( s) by releasing the electromagnet. CE contraction velocity and force were evaluated shortly after at  s. The values and were extracted. The experiment was performed with different external forces, ten repetitions each. The curve versus for all external forces represents the force-velocity characteristics of the artificial CE (Figure 3(a), crosses).

2.4. Representing the Variety of Biological Muscles

In a further evaluation of our theoretical approach we scaled the model parameters to represent various biological muscles of different animals (Figure 3(b)). The model parameters and were calculated from and values determined in experiments for different biological muscles (Table 1).

Table 1: Muscle parameters (, , , ) determined in experiments (see reference) and muscle model parameters (, ), respectively, calculated (see the appendix).

Muscle
[N]
[m/s]
[N]
[m]

Reference

Piglet gastrocnemius
3.0
0.015
30.0
0.015
0.003
2200
[17]

Cat soleus
4.8
0.042
21.0
0.033
0.011
620
[18]

Cat tenuissimus
0.05
0.057
0.18
0.032
0.600
4
[19]

Rat gastrocnemius
2.68
0.042
13.4
0.013
0.167
386
[20]

Rat tibialis anterior
4.3
0.053
4.3
0.027
0.076
162
[19]

Frog sartorius
0.18
0.012
0.67
0.031
0.287
72
[19]

2.5. Control of the Inverted Pendulum

A model of an inverted pendulum was used to investigate what effect muscle-like actuator characteristics could have on the control of robotic stance. For quiet stance, the task was to keep an upright posture while deflecting the ground to which the pendulum was suspended with a hinge joint. The model consisted of two rigid segments connected with a hinge joint (Figure 2). S1 had a mass of  kg, a moment of inertia of  kgm² (calculated around hinge joint axis), and the center of gravity was at  m. The initial orientation of the leg/trunk segment S1 was vertical and horizontal for the foot segment S2. The pendulum could be perturbed by rotating S2 about the joint by the angle . Three perturbations were considered: (a) a linear ramp increase of (for , where is the time) and (), (b) a sinusoidal oscillation , and (c) a sinusoidal oscillation .

(a)

(b)

Figure 2: Model of the inverted pendulum. S1 represents the leg-trunk segment; S2 represents the foot. COG indicates the center of gravity location of S1. is the angle of the foot (perturbation) and the deviation from the upright position of S2. (a) The joint is actuated by a direct torque generator with linear characteristics. (b) The joint is actuated by two antagonistic muscles.

(a)

(b)

Figure 3: (a) Ten and plots for quick-release contraction experiments using 19 different external forces were recorded. Based on those and plots the force-velocity curve (crosses) was calculated. In direct comparison with the biological experiments (rat gastrocnemius muscle [20]) and the predictions from theory, the artificial muscle data shows a good match for both, with . A hyperbola fit of the artificial muscle data results in [13, 15]. (b) The strength of the presented approach is shown by a comparison of curves calculated for different biological muscle parameters. The respective - curve can be plotted by just taking and values from experiments and calculating the parameters and . The latter two parameters are necessary to build a technical muscle of that respective type.

The hinge joint could be actuated either by a direct torque generator or by an antagonistic pair of muscles (Figures 2(a), and 2(b)). The muscles were represented by two macroscopic muscle models of the same type. These muscle models incorporate the contraction dynamics, as well as a serial and a parallel elastic element representing the tendon and the passive elastic properties of soft muscle tissue. The muscle model was already described in detail in [11]. The parameters used for the muscle models are listed in Table 2. Both muscles were connected to a simple geometry as depicted in Figure 2.

Table 2: Parameters for the muscle model were based on a human tibialis anterior muscle (see [21] for a detailed description). The muscle model used for the study was described in detail in [17].

0.1 m
10000 N
0.57
4.0
0.25
2.25 
0.23 m
0.1825
0.073
10000 N
0.3
0.01

Muscles and direct torque generator were controlled based on a feedback signal measuring the deviation of segment S1 from the vertical orientation. A physiological delay of  s was considered. Three different controllers were applied: no feedback is provided, a simple proportional feedback (P controller), and a PID controller. MATLAB Simulink embedded ODE5 (Dormand-Prince) solver with 1 ms step size was used to solve the differential equations.

3. Results

The relation between muscle output force and its contraction velocity is the common criterion for the comparison of macroscopic muscle models. Therefore, we calculated the - curve (Figure 3(a)) first. The - curve of the presented functional artificial muscle shows a very good match with both the prediction from theory and biological experiments.

By comparing our artificial muscle prototype’s force-velocity relation as shown above, we consider our approach as quite successful. The functional artificial muscle prototype exhibits contraction dynamics similar to Hill’s model characteristics (Figure 3(a) [13, 15]).

In a model of the inverted pendulum, muscle-like nonlinear actuator characteristics were compared against a direct torque generator (linear characteristics). The muscle-driven model did not topple, not even without feedback (first row, Figure 4). Also, when using the simple P controller (middle row) the muscle-driven model performed better during all perturbations and was able to cope better with the feedback delay of  s. Using the PID controller solutions were only found, where the direct torque controller performed better (bottom row, Figure 4). Here, a solution with high gain for the integral part of the PID controller is presented. Therefore, slow perturbations could be compensated very effectively.

(a)

(b)

(c)

Figure 4: Model reactions to perturbations in foot orientation . Control target is the upright posture (). Left column shows the reaction to a ramp perturbation, middle column to a 1 Hz, and right column to a 0.1 Hz sinusoidal perturbation. Top row is without feedback, middle row with a simple P controller (direct torque controller gain: P 500; muscle controller gain: P 1), and bottom row with a PID controller (direct torque controller gains: P 500, I 50, D 500; muscle controller gains: P 1, I 0.3, D 0.3). The feedback signal is delayed by  s.

4. Discussion

4.1. Element Representation

A brushless dc electric drive was used for the active element (AE) which was formulated in theory [11]. The tradeoffs of these actuators are the torque-to-weight ratio and the necessity of a power supply, either over cable or by battery. Madden [22] gives an overview of the current state of the art of technical artificial muscles, their potential, and their tradeoffs. For our concept as of today, we are planning to use translational drives. Translational drives directly couple the driving forces to the movement direction and they are commercially available. However, for all electric drives one challenge remains: the storage of energy. Fortunately (or unfortunately), this is also a big issue in the automotive industry for the construction of electric vehicles. Therefore, we think that it is likely to see great improvement in the storage technology in the near future. This would also facilitate the use of electric drives as active element in functional artificial muscles.

The passive damping element (PDE) develops forces during the contraction, which may exceed over a very short period of time, the muscle output forces several times over. The question remains whether there are comparable forces internal to the artificial muscle in other technical embodiments. Unfortunately, this is commonly not reported in the literature. In our approach, we use an electric drive to produce the damping forces that are in fact a nonpassive damper. Are there any materials or other approaches possible instead of using an electric drive as presented in this approach? A magnetorheological damping element, for example, [23]?

Fortunately, the serial element (SE) seems to be the simplest challenge for a technical representation. This element should imply nonlinear force-displacement characteristics. Even a steel rod would show the dynamic characteristics similar to that of the serial element predicted in theory. However, as a must-have, this element needs to incorporate damping characteristics, yet very small [17]. It is to see how the artificial muscle prototype behaves when including a serial element like that observed in biology and postulated in theory [11].

4.2. What is Gained Using This Approach?

Understanding biological muscle characteristics and design is of great interest in biological science. Muscle models in general help to mathematically formulate muscle characteristics. The structure of our model is in essence purely mechanical. Therefore, it can serve as a functional starting point of bionic muscle design. Phenomenological models based on biological experiments were the first to define muscle characteristics, for example, [16]. Constantly improving lab techniques allowed to observe muscle phenomena (even) in greater detail, for example, [5]. Microscopic muscle models deduced from basic assumptions of muscle structure and/or functional relationships of single variables come into play shortly after, for example, [24]. However, the benchmark of muscle dynamics used for those microscopic models is still the phenomenological Hill relation [16]. One approach just recently succeeded in defining the macroscopic muscle characteristics without the need of any phenomenological information. In contrast, this approach was validated against the well-known experiments instead of being based on it [11]. Here, we used those findings to build a technical muscle and succeeded in the reproduction of crucial characteristics of biological muscle. With this approach, now, the macroscopic model can be iteratively improved accompanied by the technical muscle. In that, technical models can partly replace biological experiments.

4.3. Hill-Type Models for Control

Hill-type muscle models, as an alternative to joint torque generators, have been implemented in (multibody) computer models in order to generate movement. In this regard different control theories, that is, physiologically motivated ones, for example, equilibrium point hypothesis [25, 26], virtual model control [27], and others described above, come into operation. Hence, multibody models with Hill-type muscles as actuators allow for using control theories to generate movement. This way, control approaches [28, 29] are quantitatively tested and relevant control parameters [30] are determined. Furthermore, existing and/or newly developed control theories are compared.

In this study, different control approaches, that is, no feedback, P controller, and PID controller (see Section 2), were implemented and compared in two different inverted pendulums models, that is, one with muscles and the other one with direct torque generators. From this comparison of control and actuators, it can be concluded that the implementation of muscle-like characteristics changes the model’s inherent stability. Thus, it leads to a modification of successful control strategies to generate a similar movement. Furthermore, the presented arrangement of technical elements for the CE also allows for the investigation of structural changes in biological muscle used for movement control.

For further and more detailed conclusion, the presented approach will be implemented as muscle-like actuators in more complex (human) models to investigate (physiological motivated) control strategies and structural changes of muscle.

Appendix

A. Model Derivation

In order to end up with a hyperbolic relation, two further assumptions had to be made. First, the force of the PDE was assumed to be

where the damping coefficient of the PDE depends linearly on the current muscle force :

is the maximum (at ) and the normalised (to ) minimum (force independent) value of . Second, the gearing ratio

between internal (SE) and external (muscle) velocities was represented by an arbitrary parameter value .

The characteristics of the SE did not have to be specified. The AE is the source of mechanical energy. It may depend on length and on the macroscopic chemical state of the muscle, that is, the relative number of actively force-producing cross-bridges quantified by the normalised muscle activation .

In order to meet the conditions of our artificial muscle experiments presented in this paper, we had to modify the just reviewed model [11] with respect to only one feature. In contrast to   in [11], which related the isometric force (see (2) for ) as a linear function of contraction velocity to the AE force , we now assume the identity

Equation (A.4) is as consistent to the set of model equations (2), (3), (A.1), (A.2), and (A.3) as is in [11] to this very set.

Substituting (A.1), the explicit dependency of on force and model parameters (A.2), and (A.3) into (2) makes the latter force equilibrium (2) constitute a hyperbola

with the Hill parameters , and the isometric force being positive and consistently being negative in the shortening (concentric) case. The Hill parameters are

with the corresponding maximum shortening velocity

The unloaded contractile element () would contract concentrically with :

which is defined as the Hill parameter normalised to the current isometric force . Note that, for given , a concurrent parameter variation fulfilling meets the constraint , whereas solely the curvature is changed. In our model, this is equivalent to .

Acknowledgments

This work was supported by a Research Seed Capital (RisSC)—Tranche 2009 from the Ministry of Science, Research and Arts of Baden-Wuerttemberg, The University of Stuttgart (Kapitel 1403 Tit.Gr. 74), by the German Research Foundation (DFG) within the Cluster of Excellence in Simulation Technology (EXC310/1) at the University of Stuttgart, and by the German Research Foundation (DFG) within the funding programme Open Access Publishing.

References

A. F. Huxley, “Muscle structure and theories of contraction,” Progress in Biophysics and Biophysical Chemistry, vol. 7, pp. 255–318, 1957.

A. F. Huxley, “Mechanics and models of the myosin motor,” Philosophical Transactions of the Royal Society B, vol. 355, no. 1396, pp. 433–440, 2000.

H. E. Huxley, “Fifty years of muscle and the sliding filament hypothesis,” European Journal of Biochemistry, vol. 271, no. 8, pp. 1403–1415, 2004.

V. Lombardi, G. Piazzesi, M. A. Ferenczi, H. Thirlwell, I. Dobbie, and M. Irving, “Elastic distortion of myosin heads and repriming of the working stroke in muscle,” Nature, vol. 374, no. 6522, pp. 553–555, 1995.

V. Lombardi, G. Piazzesi, M. Reconditi et al., “X-ray diffraction studies of the contractile mechanism in single muscle fibres,” Philosophical Transactions of the Royal Society B, vol. 359, no. 1452, pp. 1883–1893, 2004.

G. Piazzesi and V. Lombardi, “A cross-bridge model that is able to explain mechanical and energetic properties of shortening muscle,” Biophysical Journal, vol. 68, no. 5, pp. 1966–1979, 1995.

M. Reconditi, M. Linari, L. Lucii et al., “The myosin motor in muscle generates a smaller and slower working stroke at higher load,” Nature, vol. 428, no. 6982, pp. 578–581, 2004.

I. A. Telley, J. Denoth, and K. W. Ranatunga, “Inter-sarcomere dynamics in muscle fibres: a neglected subject?” Advances in Experimental Medicine and Biology, vol. 538, pp. 481–500, 2004.

L. Tskhovrebova and J. Trinick, “Role of titin in vertebrate striated muscle,” Philosophical Transactions of the Royal Society B, vol. 357, no. 1418, pp. 199–206, 2002.

R. H. Baughman, “Materials science. Playing nature's game with artificial muscles,” Science, vol. 308, no. 5718, pp. 63–65, 2005.

M. Günther and S. Schmitt, “A macroscopic ansatz to deduce the Hill relation,” Journal of Theoretical Biology, vol. 263, no. 4, pp. 407–418, 2010.

D. F. B. Häufle, M. Günther, R. Blickhan, and S. Schmitt, “Proof of concept of an artificial muscle: theoretical model, numerical model, and hardware experiment,” in Proceedings of International Conference of Rehabilitation Robotics (ICORR '11), pp. 1–6, July 2011.

D. F. Häufle, M. Günther, R. Blickhan, and S. Schmitt, “Can quick release experiments reveal the muscle structure? A bionic approach,” Journal of Bionic Engineering, vol. 9, no. 2, pp. 211–223, 2012.

D. F. Häufle, M. Günther, R. Blickhan, and S. Schmitt, “Proof of concept: model based bionic muscle with hyperbolic force-velocity relation,” Applied Bionics and Biomechanics, pp. 1–7, 2012.

S. Schmitt, D. F. Häufle, R. Blickhan, and M. Günther, “Nature as an engineer: one simple concept of a bio-inspired functional artificial muscle,” Bioinspiration & Biomimetics, vol. 7, Article ID 036022, 2012.

A. V. Hill, “The heat of shortening and the dynamic constants of muscle,” Proceedings of the Royal Society of London. Series B, vol. 126, pp. 136–195, 1938.

M. Günther, S. Schmitt, and V. Wank, “High-frequency oscillations as a consequence of neglected serial damping in Hill-type muscle models,” Biological Cybernetics, vol. 97, no. 1, pp. 63–79, 2007.

G. K. Cole, T. van den Bogert, W. Herzog, and K. Gerritsen, “Modelling of force production in skeletal muscle undergoing stretch,” Journal of Biomechanics, vol. 29, pp. 1091–1104, 1996.

G. K. Klute, J. M. Czerniecki, and B. Hannaford, “Artificial muscles: actuators for biorobotic systems,” International Journal of Robotics Research, vol. 21, no. 4, pp. 295–309, 2002.

J. P. van Zandwijk, M. F. Bobbert, G. C. Baan, and P. A. Huijing, “From twitch to tetanus: performance of excitation dynamics optimized for a twitch in predicting tetanic muscle forces,” Biological Cybernetics, vol. 75, no. 5, pp. 409–417, 1996.

M. Günther and H. Ruder, “Synthesis of two-dimensional human walking: a test of the λ-model,” Biological Cybernetics, vol. 89, no. 2, pp. 89–106, 2003.

J. D. Madden, “Mobile robots: motor challenges and materials solutions,” Science, vol. 318, no. 5853, pp. 1094–1097, 2007.

E. Garcia, J. C. Arevalo, G. Muoz, and P. Gonzalez-De-Santos, “Combining series elastic actuation and magneto-rheological damping for the control of agile locomotion,” Robotics and Autonomous Systems, vol. 59, no. 10, pp. 827–839, 2011.

A. F. Huxley, “A note suggesting that the cross bridge attachment during muscle contraction may take place in two stages,” Proceedings of the Royal Society of London. Series B, vol. 183, no. 1070, pp. 83–86, 1973.

E. Bizzi, N. Hogan, F. A. Mussa-Ivaldi, and S. Giszter, “Does the nervous sytem use equilibriumpoint control to guide single and multiple joint movements?” Behavioral and Brain Sciences, vol. 15, pp. 603–613, 1992.

A. G. Fel'dman, “Functional tuning of the nervous system with control of movement or maintenance of a steady posture-II. Controllable parameters of the muscles,” Biophysics, vol. 11, no. 3, pp. 565–578, 1966.

J. Pratt, C.-M. Chew, A. Torres, P. Dilworth, and G. Pratt, “Virtual model control: an intuitive approach for bipedal locomotion,” International Journal of Robotics Research, vol. 20, no. 2, pp. 129–143, 2001.

M. F. Eilenberg, H. Geyer, and H. Herr, “Control of a powered ankle-foot prosthesis based on a neuromuscular model,” IEEE Transactions on Neural Systems and Rehabilitation Engineering, vol. 18, no. 2, pp. 164–173, 2010.

A. S. Ghafari, A. Meghdari, and G. Vossoughi, “Feedback control of the neuromusculoskeletal system in a forward dynamics simulation of stair locomotion,” Proceedings of the Institution of Mechanical Engineers, Part H: Journal of Engineering in Medicine, vol. 223, no. 6, pp. 663–675, 2009.

D. A. Kistemaker, A. J. Van Soest, and M. F. Bobbert, “A model of open-loop control of equilibrium position and stiffness of the human elbow joint,” Biological Cybernetics, vol. 96, no. 3, pp. 341–350, 2007.