Numerical Methods and Applications in Biomechanical ModelingView this Special Issue
Theoretical Hill-Type Muscle and Stability: Numerical Model and Application
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.
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 .
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 . 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” , 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  it was demonstrated that the phenomenologically found hyperbolic force-velocity relation of a concentrically contracting assembly of activated muscle fibres  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.
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) . 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).
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 kgm2 (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 .
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 . 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.
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.
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.
4.1. Element Representation
A brushless dc electric drive was used for the active element (AE) which was formulated in theory . 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  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, ?
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 . It is to see how the artificial muscle prototype behaves when including a serial element like that observed in biology and postulated in theory .
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, . Constantly improving lab techniques allowed to observe muscle phenomena (even) in greater detail, for example, . Microscopic muscle models deduced from basic assumptions of muscle structure and/or functional relationships of single variables come into play shortly after, for example, . However, the benchmark of muscle dynamics used for those microscopic models is still the phenomenological Hill relation . 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 . 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 , 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  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.
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  with respect to only one feature. In contrast to in , 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  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 .
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.
A. F. Huxley, “Muscle structure and theories of contraction,” Progress in Biophysics and Biophysical Chemistry, vol. 7, pp. 255–318, 1957.View at: Google Scholar
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.View at: Google Scholar
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.View at: Google Scholar
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.View at: Google Scholar
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.View at: Google Scholar
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.View at: Publisher Site | Google Scholar
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.View at: Google Scholar
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.View at: Google Scholar
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.View at: Google Scholar
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.View at: Google Scholar
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.View at: Google Scholar
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.View at: Google Scholar
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.View at: Google Scholar
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.View at: Google Scholar
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.View at: Publisher Site | Google Scholar