- Open Access
Contribution of non-extensor muscles of the leg to maximal-effort countermovement jumping
BioMedical Engineering OnLine volume 4, Article number: 52 (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.
A 3D neuromusculoskeletal model of the human body was constructed using DADS-3D (LMS CADSI, Coralville, Iowa, USA) with the FORTRAN-based USER.FORCE option. Detailed properties of this model have been reported in preceding studies [10, 11]. The musculoskeletal model consisted of nine rigid body segments (the head-arms-trunk (HAT), right and left upper legs, right and left lower legs, right and left feet and right and left toes) connected with frictionless joints (Figure 1). Body segmental parameter values were derived from  (body mass = 73.1 kg). The hip joints were modelled as ball-and-socket joints that have three degrees of freedom. The knee joints were modelled as hinge joints. The ankle joints were modelled as universal joints . The metatarsophalangeal joints were modelled as hinge joints with a tilted axis . The total number of degrees of freedom of the model was 20.
The body was actuated by two different sets of muscles to construct the "All Muscles" model and the "Extensors Only" model. Thirty-two Hill-type lower limb muscles (16 muscles in each leg) were implemented in the "All Muscles" model (Table 1). These include all (more precisely, most of) major muscles found in a human leg. Fourteen muscles (7 muscles in each leg) were implemented in the "Extensors Only" model (Table 1). These include major leg extensor muscles only (i.e., muscles whose primary function is leg extension). Note that such biarticular muscles as the hamstrings, m. rectus femoris and m. gastrocnemius were regarded as extensor muscles. Under the joint configurations assumed during jumping, these muscles do develop more joint extension moments than joint flexion moments.
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.
A muscle-tendon complex was composed of a contractile element (CE) and a series elastic element (SEE) serially connected with a pennation angle (αpen) (Figure 2). The mathematical model of the contractile element represented the force-length-velocity relations. Passive stress-strain property of the series elastic element was modelled with a quadratic function. A detailed mathematical representation of these components can be found in . Neural activation input to individual muscles was represented by a series of step functions with duration of 0.050 s . Excitation dynamics of the contractile element was modelled with a first-order ordinary differential equation as described in .
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 maximal height reached by the mass centre of the body measured from the floor was 1.317 m for the "All Muscles" model and 1.233 m for the "Extensors Only" model (Table 2). The jumping height measured from the starting posture was 0.386 m and 0.301 m, respectively. With the body mass of 73.1 kg and the gravitational acceleration of 9.81 m/s2, the energy gain of the mass centre of the body throughout the jumping motion was 277 J and 216 J, respectively.
Realistic kinematics of jumping was generated both for the "All Muscles" model and for the "Extensors Only" model. Sagittal views of the kinematics are presented as Figure 3. The motions of the segments/joints outside of the sagittal plane were small (~10 deg; not shown), suggesting that the motion of the skeletal system was mostly two-dimensional. Take-off occurred at 0.65 s and 0.61 s after the start of simulation for the "All Muscles" model and for the "Extensors Only" model, respectively. Ground reaction force profiles are shown in Figure 4.
In the "All Muscles" model, joint flexor muscles such as the m. iliopsoas, m. biceps femoris short head and m. tibialis anterior were activated in the beginning of the countermovement phase (Figure 5). For the hamstrings, mm. vasti and other plantar flexor muscles, force output was greater for the "Extensors Only" model than for the "All Muscles" model. For the m. rectus femoris, m. gastrocnemius and m. soleus, force output was greater for the "All Muscles" model than for the "Extensors Only" model.
For the "All Muscles" model, non-extensor muscles such as hip adductors and external rotators performed relatively little work (Table 3, Figure 6), although the sum of the work outputs was substantial (47.0 J). The behaviour of the hamstrings was markedly different between the "All Muscles" model and the "Extensors Only" model. Specifically, for the "All Muscles" model, the hamstrings exerted relatively small magnitude of positive work suggesting that the action of this muscle was mostly isometric. On the other hand, for the "Extensors Only" model, this muscle had relatively large negative work suggesting that the action of this muscle was mostly eccentric (Table 3, Figure 6).
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.
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-S
Bobbert MF, Van Soest AJ: Effects of muscle strengthening on vertical jump height: a simulation study. Med Sci Sports Exerc 1994, 26: 1012–1020.
Pandy MG, Zajac FE: Optimal Muscular Coordination Strategies for Jumping. J Biomech 1991, 24: 1–10. 10.1016/0021-9290(91)90321-D
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-E
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-H
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.
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.102791
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.
Anderson FC, Pandy MG: A Dynamic Optimization Solution for Vertical Jumping in Three Dimensions. Comput Methods Biomech Biomed Engin 1999, 2: 201–231.
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.10408
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.006
de Leva P: Adjustments to Zatsiorsky-Seluyanov's segment inertia parameters. J Biomech 1996, 29: 1223–1230. 10.1016/0021-9290(95)00178-6
Inman VT: The joints of the ankle. Edited by: Baltimore MD. USA: Williams and Wilkins; 1976.
Friederich JA, Brand RA: Muscle fiber architecture in the human lower limb. J Biomech 1990, 23: 91–95. 10.1016/0021-9290(90)90373-B
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/s002210050439
Nagano A, Gerritsen KGM: Effects of neuromuscular strength training on vertical jumping performance – A computer simulation study. J Appl Biomech 2001, 17: 113–128.
Nagano A, Fukashiro S, Komura T: Contribution of series elasticity in human cyclic heel-raise exercise. J Appl Biomech 2003, 19: 340–352.
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.73565
Bremermann H: A method of unconstrained optimization. Math Biosci 1970, 9: 1–15. 10.1016/0025-5564(70)90087-8
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.
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/BF00635880
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).
AN constructed the simulation model, performed computer simulation and other calculations, and drafted the manuscript. TK and SY participated in the process of model construction and numerical optimization. SF contributed valuable discussions and suggestions throughout this project, including the stage of manuscript writing. All authors read and approved the final manuscript.
Authors’ original submitted files for images
About this article
Cite this article
Nagano, A., Komura, T., Yoshioka, S. et al. Contribution of non-extensor muscles of the leg to maximal-effort countermovement jumping. BioMed Eng OnLine 4, 52 (2005). https://doi.org/10.1186/1475-925X-4-52
- Ground Reaction Force
- Flexor Muscle
- Work Output
- Trunk Segment
- Joint Reaction Force