Contribution of non-extensor muscles of the leg to maximal-effort countermovement jumping
© Nagano et al; licensee BioMed Central Ltd. 2005
Received: 29 June 2005
Accepted: 06 September 2005
Published: 06 September 2005
The purpose of this study was to determine the effects of non-extensor muscles of the leg (i.e., muscles whose primary function is not leg extension) on the kinematics and kinetics of human maximal-effort countermovement jumping. Although it is difficult to address this type of question through experimental procedures, the methodology of computer simulation can be a powerful tool.
A skeletal model that has nine rigid body segments and twenty degrees of freedom was developed. Two sets of muscle models were attached to this skeletal model: all (most of) major muscles in the leg ("All Muscles" model) and major extensor muscles in the leg (i.e., muscles whose primary function is leg extension; "Extensors Only" model). Neural activation input signal was represented by a series of step functions with a step duration of 0.05 s. Simulations were started from an identical upright standing posture. The optimal pattern of the activation input signal was searched through extensive random-search numerical optimization with a goal of maximizing the height reached by the mass centre of the body after jumping up.
The simulated kinematics was almost two-dimensional, suggesting the validity of two-dimensional analyses when evaluating net mechanical outputs around the joints using inverse dynamics. A greater jumping height was obtained for the "All Muscles" model (0.386 m) than for the "Extensors Only" model (0.301 m). For the "All Muscles" model, flexor muscles developed force in the beginning of the countermovement. For the "All Muscles" model, the sum of the work outputs from non-extensor muscles was 47.0 J, which was 13% of the total amount (359.9 J). The quantitative distribution of the work outputs from individual muscles was markedly different between these two models.
It was suggested that the contribution of non-extensor muscles in maximal-effort countermovement jumping is substantial. The use of a computer simulation model that includes non-extensor muscles seems to be more desirable for the assessment of muscular outputs during jumping.
Jumping motions have been investigated by many researchers in the field of biomechanics in an effort to understand the coordination of the human body during explosive activities. A maximal-effort jumping is a suitable subject for this purpose, as the objective of a maximal-effort jumping can be defined in a very straightforward manner: "jump up as high as possible". Therefore less inter-subjects and intra-subject variability of body coordination is expected. In addition, jumping motions play important roles in many athletic activities such as track and field, basketball and volleyball. Therefore it is practically valuable to understand the biomechanics of the human body during jumping.
Researchers have reported many valuable insights regarding maximal-effort jumping motions using two-dimensional computer modelling and simulation [1–5]. Typically in these studies, leg muscles that have a primary function of leg extension (e.g., the m. gluteus maximus, m. rectus femoris, hamstrings, mm. vasti, m. gastrocnemius, m. soleus) were included in the model. In other words, other leg muscles that have a different primary function (e.g., joint flexion, abduction/adduction, rotation etc.; "non-extensor muscles") often were not explicitly implemented. Although it is true that the motion of the leg is mostly extension during jumping, there is a possibility that these non-extensor muscles do contribute to a jumping performance because of their three-dimensional anatomical configuration. Especially, when looking at the location of the origin, insertion and via-points of most muscles, it is observed that three-dimensional vectors instead of two-dimensional vectors better represent their line of action [6, 7]. Therefore it is valuable to investigate whether or not these non-extensor muscles of the leg make a substantial contribution to jumping performance. (Note that muscles whose primary function is not leg extension are called "non-extensor muscles" in this paper. This nomenclature does not imply that these muscles do not contribute to leg extension at all. In fact, as the human body is a linked segmental system, the activity of a muscle can affect the actions of multiple joints/degrees of freedom in the system . This paper utilized this nomenclature for the sake of simplicity.)
For that purpose, it is likely that the use of a three-dimensional neuromusculoskeletal model instead of a two-dimensional model is more straightforward. Anderson and Pandy  expanded their research on jumping motions using a three-dimensional model. Nagano et al.  also constructed a three-dimensional musculoskeletal model of a human ancestor's body that can be scaled up to represent the musculoskeletal system of modern humans . The purpose of this study was to evaluate the contribution of non-extensor muscles of the leg to maximal-effort countermovement jumping using a three-dimensional neuromusculoskeletal model.
The muscle parameter values used in this study. The values for each muscle are shown. Fmax: maximal isometric force of the contractile element. LCEopt: optimal length of the contractile element. αpen: pennation angle. Lslack: slack (unloaded) length of the series elastic element. ILIOP: m. iliopsoas. GMAXI: m. gluteus maximus. GMEDI: m. gluteus medius. GMIN: m. gluteus minimus. ADDLO: m. adductor longus. ADDMA: m. adductor magnus. ADDBR: m. adductor brevis. HEXRO: merged hip external rotator muscles. RECTF: m. rectus femoris. HAMST: merged hamstrings. VASTI: mm. vasti. BFESH: m. biceps femoris short head. GASTR: m. gastrocnemius. TIBAN: m. tibialis anterior. SOLEU: m. soleus. OPFLE: merged monoarticular planter flexor muscles other than m. soleus. All of these muscles were implemented in the "All Muscles" model. Muscles whose primary function is leg extension are noted as "Extensor". Only these muscles were implemented in the "Extensors Only" model.
Muscles investigated in  and  were considered for implementation. In order to perform computer simulation and numerical optimization within feasible computation time, it was necessary to limit the complexity of the model. Therefore muscles that have similar biomechanical function were merged to compose a single muscle group. For example, the m. vastus medialis, m. vastus intermedialis and m. vastus lateralis were merged as mm. vasti. Muscles or muscle groups whose maximal isometric force is greater than 500 N were selected. The m. biceps femoris short head, whose Fmax is smaller than 500 N, was also selected as the only mono-articular knee flexor muscle (Table 1). Coordinates of the origin, insertion and via-points of these muscles were derived from . Muscle parameter values, i.e., the optimal contractile element length (LCEopt), maximal isometric force of the contractile element (Fmax), pennation angle (αpen) and unloaded length of the series elastic element (Lslack), were derived from  and . A specific tension value of 31.5 N/cm2  was utilized. A bilateral symmetry was assumed between the right side and the left side of the body.
The interaction between the foot segments and the ground was modelled using the same form of equation as was reported in . Passive joint properties that function to limit the joint range of motion were adopted from .
Maximal-effort countermovement jump motions were generated through computer simulation with the "All Muscles" model and the "Extensors Only" model. A simulation was initiated from an upright posture with the hip, knee and ankle joints slightly flexed (5 degrees: dorsiflexed for the ankle joint) to facilitate the generation of countermovement. Simulations were performed from exactly the same initial posture for these two models. Muscle activation input profiles were modified through Bremermann's numerical optimization  in which the jumping height was maximized. The optimal combination of the activation input profiles for the muscles was searched. The optimization process was terminated when the objective function value had not improved for 10,000 successive iterations, which corresponds to approximately 60,000 function evaluations without any improvement .
The instantaneous power output value of the contractile element (PCE) was calculated as the product of the force development (FCE) and the shortening speed (VCE; positive value for shortening and negative value for lengthening) of the contractile element:
PCE = FCE·VCE (Eq. 1)
The work output of the contractile element (WCE) was calculated as the time integration of PCE from the start of simulation through the instant of take off:
The results of the numerical optimization obtained in this study. Hinit: initial height of the mass centre of the body. Hmax: maximal height reached by the mass centre of the body. Egain: energy gain of the mass centre of the body through the jumping motion calculated from the jump height, body mass and acceleration due to gravity. Δ: difference between the values for the "All Muscles" model and for the "Extensors Only" model.
Jump Height (m)
The amount of the mechanical work performed by the contractile element. The integrated values from the start of simulation through the instant of take-off. The added values for two legs (two contralateral muscles) are shown. Δ: difference between the values for the "All Muscles" model and for the "Extensors Only" model.
CE Work (two legs) (J)
The purpose of this study was to evaluate quantitatively the contribution of non-extensor muscles (muscles whose primary function is not leg extension) of the leg to maximal-effort countermovement jumping. Details of the simulation model utilized in this study have been described in . In that study, a countermovement jumping motion simulated with the "All Muscles" mode have been analyzed and compared with the experimental data reported in preceding studies, and the validity of the modelling and simulation has been discussed. The optimized jumping height was smaller for the "Extensors Only" model than for the "All Muscles" model by 0.084 m (Table 2; 28.0%). In this study, the decrement in performance was caused by the absence of non-extensor muscles. This implies that non-extensor muscles do have substantial contributions to a maximal-effort countermovement jumping performance.
In both cases, the general characteristics of the jumping kinematics obtained through the numerical optimization process (Figure 3) were similar to the ones obtained through experimental data collection of human countermovement jumping , although the motion of the body and joint excursions were greater for the "All Muscles" model than for the "Extensors Only" model (Figure 3). Only limited motions of the skeletal system were observed outside of the sagittal plane in this study. It should be noted that the computer simulation model utilized in this study has a capability to perform fully three-dimensional motions (e.g., hip joint abduction/adduction etc.). Nonetheless, the simulation model chose to perform almost two-dimensional motions. This finding supports that the two-dimensional inverse dynamic analyses on jumping performed in numerous preceding studies are mostly valid. Especially when calculating such mechanical variables as net joint reaction forces, net joint moments and power outputs of joints, reliable calculations can be assumed.
When performing computer simulation of jumping, it is assumed to be acceptable to construct a two-dimensional skeletal model of the human body for the same reason. However, when attaching muscle models to the skeletal model, it will be more appropriate to explicitly consider the contribution of non-extensor muscles of the leg. Implementing three-dimensional configuration of these muscles will be the most straightforward solution. Calculating the projection of the line of action of these muscles to the sagittal plane will be another option to accomplish this.
Regarding the profiles of ground reaction force (Figure 4), two peaks were observed during the push-off phase. The first peak was mostly caused through the interaction between the heel and the floor in the beginning of the push-off phase, whereas the second peak was mostly caused through the interaction between the toe and the floor in the last part of the push-off phase. This profile of ground reaction force with two peaks is often observed in ground reaction force data collected from human subjects during a maximal-effort countermovement jumping [1, 21]. The profile of ground reaction force was bumpy because each foot was modelled with only five contact points . The profile will become smoother with more contact points in a foot, although this modification will greatly increase the computation time.
In the beginning of countermovement, joint flexor muscles, such as the m. iliopsoas, m. biceps femoris short head and m. tibialis anterior, developed force (Figure 5). This resulted in a greater countermovement for the "All Muscles" model than for the "Extensors Only" model. This seems to suggest that the contribution of these flexor muscles in the beginning of countermovement should be considered when investigating the mechanism of maximal-effort countermovement jumping motion. For the "All Muscles" model, non-extensor muscles such as the m. gluteus medius, m. gluteus minimus, adductors and hip external rotators had relatively minor individual contributions in terms of mechanical work and power outputs (Table 3). However, when the work outputs of these muscles were added together, the amount was substantial (47.0 J in 359.9 J; 13%), suggesting that the contribution of these muscles in jumping motion is not negligible.
For the m. gluteus maximus, m. rectus femoris, hamstrings, m. gastrocnemius and other monoarticular plantar flexor muscles (OPFLE), the work output was greater for the "All Muscles" model than for the "Extensors Only" model (Table 3). This result is very reasonable considering that the jumping height was greater for the "All Muscles" model than for the "Extensors Only" model (Table 2). Generally speaking, to achieve a higher jumping performance in an optimally-coordinated movement, muscles need to perform more work. As the jumping height was greater for the "All Muscles" model, greater mechanical outputs of muscles are reasonably expected for this model than for the "Extensors Only" model.
However, there were two exceptions; for the mm. vasti and m. soleus, the work output was greater for the "Extensors Only" model than for the "All Muscles" model. This result came from the fact that the "Extensors Only" model underwent a smaller countermovement compared to the "All Muscles" model (Figure 3, Figure 4). In this study, the amount of work output was calculated as a net (positive and negative) value from the start of a motion through the instant of take-off (Eq. 1 and 2). As the magnitude of countermovement (negative phase) was smaller for the "Extensors Only" model, the net amount of work output of the mm. vasti and m. soleus was calculated to be greater for this model. This phenomenon can be observed in Figure 6, where a smaller negative power output of these muscles during the countermovement is exhibited for the "Extensors Only" model. As the mm. vasti is a major knee extensor and the m. soleus is a major ankle plantarflexor, these muscles had to function to brake the downward momentum generated during the countermovement. This had the effect of reducing the net work output of these muscles.
For the hamstrings, a positive work output (10.1 J) was calculated for the "All Muscles" model, whereas a negative (-36.1 J) value was calculated for the "Extensors Only" model. It is observed that the contractile element of the hamstrings was mostly stretched in an eccentric manner in the "Extensors Only" model in the latter phase of the countermovement, resulting in a negative power (Figure 6) and work (Table 3) outputs. This is because only a few muscles that can act as extensors were available to brake the countermovement of the trunk segment in this model. Specifically, the m. adductor longus, m. adductor brevis and m. adductor magnus had been removed from the model. Therefore the inertial load of the trunk segment (moving downwards) imposed on the hamstrings in the latter phase of the countermovement was so great as to stretch this muscle in an eccentric manner, although this muscle was vigorously activated during this period (Figure 5). In other words, the hamstrings was not strong enough to brake the downward momentum of the trunk segment in a concentric manner. This discussion is consistent with the muscle force development profile shown in Figure 5. The force development of the hamstrings was greater for the "Extensors Only" model than for the "All Muscles" model (Figure 5), which is reasonable considering that the action of this muscle was mostly eccentric for the "Extensors Only" model (eccentric part of the force-velocity relation; Figure 2). On the other hand, as the "All Muscles" model had more muscles to function to brake the countermovement, the hamstrings could shorten itself and produce positive (concentric) work and power outputs. For example, the m. adductor magnus did function to brake the countermovement in this model (Figure 6). This result is consistent with the observation that there was a greater countermovement for the "All Muscles" model than for the "Extensors Only" model. As there were not enough muscles to brake the downward momentum generated during the countermovement, the optimal magnitude of countermovement for the "Extensors Only" model was smaller than that for the "All Muscles" model.
As a result of this computer simulation study, it was found that the dynamics of the body motion is altered by the effects of non-extensor muscles. This finding is noteworthy considering that the overall kinematics of the body (Figure 3) and the ground reaction force profiles (Figure 4) were similar between the "All Muscles" model and the "Extensors Only" model. This result implies that it is desirable to consider explicitly the mechanical contribution of non-extensor muscles of the leg when investigating human jumping motions in terms of mechanical outputs of muscles. The use of a three-dimensional neuromusculoskeletal model seems to be more suitable for this purpose. However, results of this computer simulation study also supported that the nature of an optimally-coordinated countermovement jumping motion is mostly two-dimensional, which suggests the validity of the two-dimensional inverse dynamic analyses of net mechanical outputs around joints performed in many preceding studies.
AN would like to thank the special post-doctoral program of RIKEN. SY would like to thank the Junior Research Associate program of RIKEN. AN and SY would like to thank Professor Ryutaro Himeno at RIKEN for his supports. This study was partly supported by the Ministry of Education, Culture, Sports, Science and Technology in Japan (No: 16300205).
- Anderson FC, Pandy MG: Storage and Utilization of Elastic Strain-Energy during Jumping. J Biomech 1993, 26: 1413–1427. 10.1016/0021-9290(93)90092-SView ArticleGoogle Scholar
- Bobbert MF, Van Soest AJ: Effects of muscle strengthening on vertical jump height: a simulation study. Med Sci Sports Exerc 1994, 26: 1012–1020.View ArticleGoogle Scholar
- Pandy MG, Zajac FE: Optimal Muscular Coordination Strategies for Jumping. J Biomech 1991, 24: 1–10. 10.1016/0021-9290(91)90321-DView ArticleGoogle Scholar
- Pandy MG, Zajac FE, Sim E, Levine WS: An Optimal-Control Model for Maximum-Height Human Jumping. J Biomech 1990, 23: 1185–1198. 10.1016/0021-9290(90)90376-EView ArticleGoogle Scholar
- van Soest AJ, Schwab AL, Bobbert MF, van Ingen Schenau GJ: The influence of the biarticularity of the gastrocnemius muscle on vertical-jumping achievement. J Biomech 1993, 26: 1–8. 10.1016/0021-9290(93)90608-HView ArticleGoogle Scholar
- Brand RA, Crowninshield RD, Wittstock CE, Pedersen DR, Clark CR, van Krieken FM: A model of lower extremity muscular anatomy. J Biomech Eng 1982, 104: 304–310.View ArticleGoogle Scholar
- Delp SL, Loan JP, Hoy MG, Zajac FE, Topp EL, Rosen JM: An interactive graphics-based model of the lower extremity to study orthopaedic surgical procedures. IEEE Trans Biomed Eng 1990, 37: 757–767. 10.1109/10.102791View ArticleGoogle Scholar
- Yamaguchi GT: Dynamic modeling of musculoskeletal motion: a vectorized approach for biomechanical analysis in three dymensions. Edited by: Norwell MA. USA: Kluwer Academic Publishers; 2001.Google Scholar
- Anderson FC, Pandy MG: A Dynamic Optimization Solution for Vertical Jumping in Three Dimensions. Comput Methods Biomech Biomed Engin 1999, 2: 201–231.View ArticleGoogle Scholar
- Nagano A, Umberger BR, Marzke MW, Gerritsen KG: Neuromusculoskeletal computer modeling and simulation of upright, straight-legged, bipedal locomotion of Australopithecus afarensis (A.L. 288–1). Am J Phys Anthropol 2005, 126: 2–13. 10.1002/ajpa.10408View ArticleGoogle Scholar
- Nagano A, Komura T, Fukashiro S, Himeno R: Force, work and power output of lower limb muscles during human maximal-effort countermovement jumping. J Electromyogr Kinesiol 2005, 15: 367–376. 10.1016/j.jelekin.2004.12.006View ArticleGoogle Scholar
- de Leva P: Adjustments to Zatsiorsky-Seluyanov's segment inertia parameters. J Biomech 1996, 29: 1223–1230. 10.1016/0021-9290(95)00178-6View ArticleGoogle Scholar
- Inman VT: The joints of the ankle. Edited by: Baltimore MD. USA: Williams and Wilkins; 1976.Google Scholar
- Friederich JA, Brand RA: Muscle fiber architecture in the human lower limb. J Biomech 1990, 23: 91–95. 10.1016/0021-9290(90)90373-BView ArticleGoogle Scholar
- Brown IE, Satoda T, Richmond FJR, Loeb GE: Feline caudofemoralis muscle – Muscle fibre properties, architecture, and motor innervation. Exp Brain Res 1998, 121: 76–91. 10.1007/s002210050439View ArticleGoogle Scholar
- Nagano A, Gerritsen KGM: Effects of neuromuscular strength training on vertical jumping performance – A computer simulation study. J Appl Biomech 2001, 17: 113–128.Google Scholar
- Nagano A, Fukashiro S, Komura T: Contribution of series elasticity in human cyclic heel-raise exercise. J Appl Biomech 2003, 19: 340–352.Google Scholar
- He JP, Levine WS, Loeb GE: Feedback Gains for Correcting Small Perturbations to Standing Posture. Ieee Transactions on Automatic Control 1991, 36: 322–332. 10.1109/9.73565View ArticleGoogle Scholar
- Bremermann H: A method of unconstrained optimization. Math Biosci 1970, 9: 1–15. 10.1016/0025-5564(70)90087-8MathSciNetView ArticleGoogle Scholar
- Bobbert MF, Gerritsen KG, Litjens MC, Van Soest AJ: Why is countermovement jump height greater than squat jump height? Med Sci Sports Exerc 1996, 28: 1402–1412.View ArticleGoogle Scholar
- Fukashiro S, Komi PV, Jarvinen M, Miyashita M: In vivo Achilles tendon loading during jumping in humans. Eur J Appl Physiol Occup Physiol 1995, 71: 453–458. 10.1007/BF00635880View 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.