Inertial characteristics of adolescent male body segments

In response to the presently limited information on body segment inertial characteristics of children and adolescents this investigation estimated the mass, centre of mass and principal moments of inertia of adolescent male body segments. Significant prediction equations based on anthropometric measurements were then sought. Thirteen subjects were measured at 6-monthly intervals for 2.5 yr to provide inertial characteristics for the leg, thigh, lower trunk and upper trunk segments. These characteristics were derived using an elliptical zone modelling technique. Following a correlation analysis, significant prediction equations of segment inertial parameters were derived from five, or fewer, anthropometric measurements. For all cases, more than 84% of the variance in the dependent variable was accounted for with a maximum R2 value of 94% being recorded for the prediction of thigh segment mass. The use of these prediction equations offered accurate and convenient estimates of body segment inertial characteristics within the limitations applicable to all modelling approaches. In contrast to recent studies, these equations accommodated the current morphological status of the subject.     

Measuring morphological parameters of the pelvic floor for finite element modelling purposes

The goal of this study was to obtain a complete data set needed for studying the complex biomechanical behaviour of the pelvic floor muscles using a computer model based on the finite element (FE) theory. The model should be able to predict the effect of surgical interventions and give insight into the function of pelvic floor muscles. Because there was a lack of any information concerning morphological parameters of the pelvic floor muscle structures, we performed an experimental measurement to uncover those morphological parameters. Geometric parameters as well as muscle parameters of the pelvic floor muscles were measured on an embalmed female cadaver. A three-dimensional (3D) geometric data set of the pelvic floor including muscle fibre directions was obtained using a palpator device. A 3D surface model based on the experimental data, needed for mathematical modelling of the pelvic floor, was created. For all parts of the diaphragma pelvis, the optimal muscle fibre length was determined by laser diffraction measurements of the sarcomere length. In addition, other muscle parameters such as physiological cross-sectional area and total muscle fibre length were determined. Apart from these measurements we obtained a data set of the pelvic floor structures based on nuclear magnetic resonance imaging (MRI) on the same cadaver specimen. The purpose of this experiment was to discover the relationship between the MRI morphology and geometrical parameters obtained from the previous measurements. The produced data set is not only important for biomechanical modelling of the pelvic floor muscles, but it also describes the geometry of muscle fibres and is useful for functional analysis of the pelvic floor in general. By the use of many reference landmarks all these morphologic data concerning fibre directions and optimal fibre length can be morphed to the geometrical data based on segmentation from MRI scans.These data can be directly used as an input for building a mathematical model based on FE theory.     

Short-term repeatability of parameters extracted from radial displacement of muscle belly.

The aim of this study is to analyse the short-term repeatability of the parameters extracted from radial muscle belly response on a stimulation pulse. The method uses a prefixed tension of a displacement sensor tip to the muscle and is being developed for noninvasive and selective evaluation of skeletal muscle contraction properties. Five parameters were extracted and statistically evaluated from the measured displacement: maximal displacement, delay time, contraction time, sustain time and half relaxation time. Care has been taken to leave sufficient time between stimulation pulses in order to reduce the effect of muscle fatigue and a constant pre-tension was assured by controlled step motor in consecutive measurement by withdrawal and anew placement of the sensor to the muscle belly after each measurement. Intra-class correlation coefficient (ICC) and normalized squared error (NSEM) were used as measures of short-term repeatability and accuracy (precision) of the measurements while fatigue rate was evaluated using area ratio fatigue index and normalized slope. All five measured parameters have been found highly repeatable (ICC from 0.86 to 0.98) and can be measured with high precision (NSEM from 0.43 to 1.93). Maximal displacement and half relaxation time show largest influence to muscle fatigue rate and are also expected to be the best measure of the fatigue rate. This investigation should serve as an initial study of repeatability of the presented method that should help in subsequent investigations and use of the method.     

Inter-rater reliability of muscle contractile property measurements using non-invasive tensiomyography

Tensiomyography (TMG) is a relatively novel technique to assess muscle mechanical response based on radial muscle belly displacement consecutive to a single electrical stimulus. Although intra-session reliability has been found to be good, inter-rater reliability and the influence of sensor repositioning and electrodes placement on TMG measurements is unknown. The purpose of this study was to analyze the inter-rater reliability of vastus medialis muscle contractile property measurements obtained with TMG as well as the effect of inter-electrode distance (IED). Five contractile parameters were analyzed from vastus medialis muscle belly displacement–time curves: maximal displacement (Dm), contraction time (Tc), sustain time (Ts), delay time (Td), and half-relaxation time (Tr). The inter-rater reliability and IED effect on these measurements were evaluated in 18 subjects. Intra-class correlation coefficients, standard errors of measurement, Bland and Altman systematic bias and random error as well as coefficient of variations were used as measures of reliability. Overall, a good to excellent inter-rater reliability was found for all contractile parameters, except Tr, which showed insufficient reliability. Alterations in IED significantly affected Dm with a trend for all the other parameters. The present results legitimate the use of TMG for the assessment of vastus medialis muscle contractile properties, particularly for Dm and Tc. It is recommended to avoid Tr quantification and IED modifications during multiple TMG measurements.     

Increased cardiovascular response to static contraction of larger muscle groups

The purpose of this study was to investigate the influence of the size of the active muscle mass on the cardiovascular response to static contraction. Twelve male subjects performed one-arm handgrip (HG), two-leg extension (LE), and a "dead-lift" maneuver (DL) in a randomly assigned order for 3 min at 30% of maximal voluntary contraction. O2 uptake (VO2), heart rate (HR), and mean intra-arterial blood pressure (MABP) were measured at rest and, in addition to absolute tension exerted, throughout contraction. There was a direct relationship between the size of the active muscle mass and the magnitude of the increases in VO2, HR, and MABP, even though all contractions were performed at the same relative intensity. Tension, VO2, HR, and MABP increased progressively from HG to LE to DL. It was concluded that at the same percentage of maximal voluntary contraction, the magnitude of the cardiovascular response to isometric exercise is directly influenced by the size of the contracting muscle mass.     

High-frequency oscillations as a consequence of neglected serial damping in Hill-type muscle models

High-frequency vibrations e.g., induced by legs impacting with the ground during terrestrial locomotion can provoke damage within tendons even leading to ruptures. So far, macroscopic Hill-type muscle models do not account for the observed high-frequency damping at low-amplitudes. Therefore, former studies proposed that protective damping might be explained by modelling the contractile machinery of the muscles in more detail, i.e., taking the microscopic processes of the actin–myosin coupling into account. In contrast, this study formulates an alternative hypothesis: low but significant damping of the passive material in series to the contractile machinery—e.g., tendons, aponeuroses, titin—may well suffice to damp these hazardous vibrations. Thereto, we measured the contraction dynamics of a piglet muscle–tendon complex (MTC) in three contraction modes at varying loads and muscle–tendon lengths. We simulated all three respective load situations on a computer: a Hill-type muscle model including a contractile element (CE) and each an elastic element in parallel (PEE) and in series (SEE) to the CE pulled on a loading mass. By comparing the model to the measured output of the MTC, we extracted a consistent set of muscle parameters. We varied the model by introducing either linear damping in parallel or in series to the CE leading to accordant re-formulations of the contraction dynamics of the CE. The comparison of the three cases (no additional damping, parallel damping, serial damping) revealed that serial damping at a physiological magnitude suffices to explain damping of high-frequency vibrations of low amplitudes. The simulation demonstrates that any undamped serial structure within the MTC enforces SEE-load eigenoscillations. Consequently, damping must be spread all over the MTC, i.e., rather has to be de-localised than localised within just the active muscle material. Additionally, due to suppressed eigenoscillations Hill-type muscle models taking into account serial damping are numerically more efficient when used in macroscopic biomechanical neuro-musculo-skeletal models.     

The simulation of aerial movement--II. A mathematical inertia model of the human body

A mathematical inertia model which permits the determination of personalized segmental inertia parameter values from anthropometric measurements is described. The human body is modelled using 40 geometric solids which are specified by 95 anthropometric measurements. A 'stadium' solid is introduced for modelling the torso segments using perimeter and width measurements. This procedure is more accurate than the use of elliptical discs of given width and depth and permits a smaller number of such solids to be used. Inertia parameter values may be obtained for body models of up to 20 segments. Errors in total body mass estimates from this and other models are discussed with reference to the unknown lung volumes.     

MRI measures of perfusion-related changes in human skeletal muscle during progressive contractions.

Although skeletal muscle perfusion is fundamental to proper muscle function, in vivo measurements are typically limited to those of limb or arterial blood flow, rather than flow within the muscle bed itself. We present a noninvasive functional MRI (fMRI) technique for measuring perfusion-related signal intensity (SI) changes in human skeletal muscle during and after contractions and demonstrate its application to the question of occlusion during a range of contraction intensities. Eight healthy men (aged 20-31 yr) performed a series of isometric ankle dorsiflexor contractions from 10 to 100% maximal voluntary contraction. Axial gradient-echo echo-planar images (repetition time = 500 ms, echo time = 18.6 ms) were acquired continuously before, during, and following each 10-s contraction, with 4.5-min rest between contractions. Average SI in the dorsiflexor muscles was calculated for all 240 images in each contraction series. Postcontraction hyperemia for each force level was determined as peak change in SI after contraction, which was then scaled to that obtained following a 5-min cuff occlusion of the thigh (i.e., maximal hyperemia). A subset of subjects (n = 4) performed parallel studies using venous occlusion plethysmography to measure limb blood flow. Hyperemia measured by fMRI and plethysmography demonstrated good agreement. Postcontraction hyperemia measured by fMRI scaled with contraction intensity up to approximately 60% maximal voluntary contraction. fMRI provides a noninvasive means of quantifying perfusion-related changes during and following skeletal muscle contractions in humans. Temporal changes in perfusion can be observed, as can the heterogeneity of perfusion across the muscle bed.     

Are fixed limb inertial models valid for dynamic simulations of human movement?

During human movement, muscle activation and limb movement result in subtle changes in muscle mass distribution. Muscle mass redistribution can affect limb inertial properties and limb dynamics, but it is not currently known to what extent. The objectives of this study were to investigate: (1) how physiological alterations of muscle and tendon length affect limb inertial characteristics, and (2) how such changes affect dynamic simulations of human movement. To achieve these objectives, a digital model of a human leg, custom software, and Software for interactive musculoskeletal modeling were used to simulate mass redistribution of muscle–tendon structures within a limb segment during muscle activation and joint movement. Thigh and shank center of mass and moments of inertia for different muscle activation and joint configurations were determined and compared. Limb inertial parameters representing relaxed muscles and fully active muscles were input into a simulated straight-leg movement to evaluate the effect inertial parameter variations could have on movement simulation results. Muscle activation and limb movement altered limb segment center of mass and moments of inertia by less than 0.04 cm and 1.2%, respectively. These variations in limb inertial properties resulted in less than 0.01% change in maximum angular velocity for a simulated straight-leg hip flexion task. These data demonstrate that, for the digital human leg model considered, assuming static quantities for segment center of masses and moments of inertia in movement simulations appear reasonable and induce minimal errors in simulated movement dynamics.     

Myoelectric signal conduction velocity and spectral parameters: influence of force and time

Reports on measurement of muscle fiber conduction velocity in humans are scarce. Inferences on the behavior of conduction velocity have been drawn from the behavior of myoelectric spectral parameters. The present report contains information on conduction velocity and spectral parameters studied at various muscle contraction levels and during and after sustained contractions. The following results have been obtained from measurements on the tibialis anterior muscle. Conduction velocity demonstrated a positive correlation with limb circumference and with muscle force output. Thus we conclude that the diameters of the muscle fibers of high-threshold motor units are, on an average, larger than those of low-threshold motor units. The study of a sustained contraction and on the recovery after such a contraction revealed that conduction velocity consistently decreased during a strong contraction as did various myoelectric spectral parameters. However, the spectral parameters decreased approximately twice as much as did the conduction velocity, and we conclude that factors other than the conduction velocity along the muscle fibers affect the myoelectric signal during a high-level contraction. These other factors appertain to changes in the firing statistics of individual motor units as well as the correlation between the firings of different motor units.     

Evaluation of disuse atrophy of rat skeletal muscle based on muscle energy metabolism assessed by 31P-MRS

The purpose of this study was to evaluate disuse atrophy of skeletal muscle using a hind-limb suspension model, with special reference to energy metabolism. Twenty-four Sprague-Dawley rats were divided into four groups: control group (C), hind-limb suspended for 3 days (HS-3), for 7 days (HS-7) and for 14 days (HS-14). The gastrocnemius-plantaris-soleus (GPS) muscles in each group were subjected to the following measurements. After a 2-min rest, contraction of the GPS muscles was induced by electrical stimulation of the sciatic nerve at 0.25 Hz for 10 min, then the frequency was increased to 0.5 and 1.0 Hz every 10 min. During the stimulation, twitch forces were recorded by a strain gauge, and 31P-MRS was performed simultaneously. Maximum tension was measured at the muscle contraction induced at 0.25 Hz; the wet weight of the whole and each muscle in the GPS muscles was also measured. From the 31P-MR spectra during muscle contraction, the oxidative capacity was calculated and compared among the groups. The weights of the whole GPS muscles in C, HS-3, HS-7 and HS-14, were 2.66 +/- 0.09, 2.39 +/- 0.21, 2.34 +/- 0.21 and 2.18 +/- 0.14 (g) respectively. Thus, the muscle mass significantly decreased with time (p < 0.05). Among the GPS muscles, the decrease in weight of the soleus muscle was especially remarkable; in the HS-14 group its weight decreased to 60% of that in the C group. We evaluated maximum tension and oxidative capacity as the muscle function. The maximum tensions in C, HS-3, HS-7 and HS-14 were 519 +/- 43, 446 +/- 66, 450 +/- 23 and 465 +/- 29 (g), respectively. This was significantly greater in the C group than in any other groups, however there were no significant differences among the three HS groups. The oxidative capacity during muscle contraction in the C group was higher than in any HS group and it did not further decrease even if the suspension of the limbs was prolonged beyond 3 days. The present study showed that in disuse atrophy, muscle mass and muscle function did not change simultaneously. Thus, it is necessary to develop countermeasures to prevent muscle atrophy and muscle function deterioration independently.     

A novel method of studying fascicle architecture in relaxed and contracted muscles

A muscle’s architecture, described by geometric variables such as fascicle pennation angles or lengths, plays a crucial role in its functionality. Usually, single parameters are used to estimate force vectors or lengthening rates, thereby assuming that they represent the architecture properly and are constant during contraction. To describe muscle architecture in more detail and compare relaxed and contracted states, we developed and validated a new approach. The m. soleus of the laboratory rat was shock-frozen while relaxed and under isometric contraction, reconstructed three-dimensionally from histological sections, and fascicle lengths, curvatures and pennation angles, as well as the shape of the aponeuroses were analysed. Remarkable differences in volume distribution and the shapes of the aponeuroses as well as locally varying changes in the fascicle architecture were observed. While the mean pennation angle increased by only 2° due to contraction, local changes of up to 4° were observed. Fascicle curvature increased in the distal but remained unchanged in the proximal parts. Our approach may help to identify functional subunits within the muscle, i.e., regions with homogeneous architectural properties. Our results are discussed regarding the input parameters essential for realistic muscle modelling and challenge maximum isometric force estimations that are based on the physiological cross-sectional area or the Hill-model.     

Muscle-contraction properties in overarm throwing movements

On the basis of dynamic and kinematic data, this study identifies the type of muscle contraction in unloaded overarm throwing movements. An unloaded throw or nearly unloaded throw is defined as the throw in which the external resistance is too small (e.g., the team handball, baseball, and water polo throws as well as the tennis and badminton smashes). A special arm-force-measuring apparatus was constructed to imitate an overarm throw. Forty-two subjects were placed into 3 groups: untrained subjects, weight-trained athletes, and team handball players. The measured parameters included the velocity of the initial movement, the release velocity, the velocity of the first 50 milliseconds of the concentric phase, the force value at the moment of deceleration of the initial movement, and the impulse values during the eccentric and concentric phases of the test movement. Statistically significant higher values of the above parameters (p < 0.05) were determined in that test at which the initial speed of movement was higher. Also, the correlation coefficients of the parameters of the initial phase of the throw movement were very high (p < 0.001), especially the parameters related with the movement's first 50 milliseconds. The results support the thesis that the stretch-shortening cycle is the type of muscle contraction in unloaded overarm throws. Furthermore, it is possible to increase the throw velocity by increasing the velocity of the initial movement (i.e., by provoking higher inertia forces).     

Local longitudinal muscle shortening of the human esophagus from high-frequency ultrasonography

We analyzed local longitudinal shortening by combining concurrent ultrasonography and manometry with basic principles of mechanics. We applied the law of mass conservation to quantify local axial shortening of the esophageal wall from ultrasonically measured cross-sectional area concurrently with measured intraluminal pressure, from which correlations between local contraction of longitudinal and circular muscle are inferred. Two clear phases of local longitudinal shortening were observed during bolus transport. During luminal filling by bolus fluid, the muscle layer distends and the muscle thickness decreases in the absence of circular or longitudinal muscle contraction. This is followed by local contraction, first in longitudinal muscle, then in circular muscle. Maximal longitudinal shortening occurs nearly coincidently with peak intraluminal pressure. Longitudinal muscle contraction begins before and ends after circular muscle contraction. Larger longitudinal shortening is correlated with higher pressure amplitude, suggesting that circumferential contractile forces are enhanced by longitudinal muscle shortening. We conclude that a peristaltic wave of longitudinal muscle contraction envelops the wave of circular muscle contraction as it passes through the middle esophagus, with peak longitudinal contraction aligned with peak circular muscular contraction. Our results suggest that the coordination of the two waves may be a physiological response to the mechanical influence of longitudinal shortening, which increases contractile force while reducing average muscle fiber tension by increasing circular muscle fiber density locally near the bolus tail.     

Muscle power output properties using the stretch-shortening cycle of the upper limb and their relationships with a one-repetition maximum bench press

The purpose of this study was to examine the output properties of muscle power by the dominant upper limb using SSC, and the relationships between the power output by SSC and a one-repetition maximum bench press (1 RM BP) used as a strength indicator of the upper body. Sixteen male athletes (21.4+/-0.9 yr) participated in this study. They pulled a load of 40% of maximum voluntary contraction (MVC) at a stretch by elbow flexion of the dominant upper limb in the following three preliminary conditions: static relaxed muscle state (SR condition), isometric muscle contraction state (ISO condition), and using SSC (SSC condition). The velocity with a wire load via a pulley during elbow flexion was measured accurately using a power instrument with a rotary encoder, and the muscle power curve was drawn from the product of the velocity and load. Significant differences were found among all evaluation parameters of muscle power exerted from the above three conditions and the parameters regarding early power output during concentric contraction were larger in the SSC condition than the SR and ISO conditions. The parameters on initial muscle contraction velocity when only using SSC significantly correlated with 1 RM BP (r=0.60-0.62). The use of SSC before powerful elbow flexion may contribute largely to early explosive power output during concentric contraction. Bench press capacity relates to a development of the above early power output when using SSC.     

A myocybernetic control model of skeletal muscle

A mathematical model of skeletal muscle is presented which contains the two physiological control parameters stimulation rate and motor unit recruitment. The model is complete in the sense that it adequately describes all possible contractive states normally occurring in living muscle. The modelling procedure relies entirely on established myo-physiological facts and each assumption made is substantiated by experimental data. Extensive simulation studies reveal that the model is capable of correctly predicting practically all known phenomena of the muscular force-output. A simplified version of the model is also presented, particularly suitable for inclusion as the driving structure in complex musculoskeletal link systems. This version was successfully tested in the prediction of an optimal human motion. The present control model is believed to fill a gap in the literature on models of muscle, and may be expected to provide a sound basis for research into the optimal control aspects of muscular contraction, and to stimulate such research.     

Muscle force estimation in clinical gait analysis using AnyBody and OpenSim

A variety of musculoskeletal models are applied in different modelling environments for estimating muscle forces during gait. Influence of different modelling assumptions and approaches on model outputs are still not fully understood, while direct comparisons of standard approaches have been rarely undertaken. This study seeks to compare joint kinematics, joint kinetics and estimated muscle forces of two standard approaches offered in two different modelling environments (AnyBody, OpenSim). It is hypothesised that distinctive differences exist for individual muscles, while summing up synergists show general agreement. Experimental data of 10 healthy participants (28 ± 5 years, 1.72 ± 0.08 m, 69 ± 12 kg) was used for a standard static optimisation muscle force estimation routine in AnyBody and OpenSim while using two gait-specific musculoskeletal models. Statistical parameter mapping paired t-test was used to compare joint angle, moment and muscle force waveforms in Matlab. Results showed differences especially in sagittal ankle and hip angles as well as sagittal knee moments. Differences were also found for some of the muscles, especially of the triceps surae group and the biceps femoris short head, which occur as a result of different anthropometric and anatomical definitions (mass and inertia of segments, muscle properties) and scaling procedures (static vs. dynamic). Understanding these differences and their cause is crucial to operate such modelling environments in a clinical setting. Future research should focus on alternatives to classical generic musculoskeletal models (e.g. implementation of functional calibration tasks), while using experimental data reflecting normal and pathological gait to gain a better understanding of variations and divergent behaviour between approaches.     

Component inertia modeling of segmental wobbling and rigid masses

A modification to an existing mathematical model is described, which permits the determination of subject-specific inertia parameters for wobbling and rigid masses of female body segments. The model comprises segment-specific soft tissue, bone, and lung components. A total of 59 geometric solids (40 soft tissue, 17 bone, 2 lung) were used to represent the body components. Ninety-five anthropometric measurements were collected from 7 female participants and were used to develop and evaluate the model. The success of the model is evaluated using predicted mass and mass distribution. The overall absolute accuracy in predicted whole body mass was better than 3.0%, with a maximum error of 4.9%. The appropriateness of the cadaver-based density values used in the model is addressed and the accuracy of the component inertia model in relation to uniform density models is discussed. The model offers a novel approach for determining component inertia parameters, which have been used successfully in a wobbling mass model to produce realistic kinetic analyses of drop-landings.     

Prehistoric mongoloid dispersals. Edited by Takeru Akazawa and Emöke J. E. Szathmáry. New York: Oxford University Press. 1996. 387 pp. ISBN 0-19-852318-1 $150.00 (cloth)

Studies of the dynamics of locomotor performances depend on knowledge of the distribution of body mass within and between limb segments. However, these data are difficult to derive. Segment mass properties have generally been estimated by modelling limbs as truncated cones, but this approach fails to take into account that some segments are of elliptical, not circular, cross section; and further, the profiles of real segments are generally curved. Thus, they are more appropriately modelled as solids of revolution, described by the rotation in space of convex or concave curves, and the possibility of an elliptical cross section needs to be taken into account. In this project we have set out to develop a general geometric model which can take these factors into account, and permit segment inertial properties to be derived from cadavers by segmentation, and from living individuals using linear external measurements. We present a model which may be described by up to four parameters, depending o the profile and serial cross section (circular or ellipsoidal) of the individual segments. The parameters are obtained from cadavers using a simplified complex-pendulum technique, and from intact specimens by calculation from measurements of segment diameters and lengths. From the parameters, the center of mass, moments of inertia, and radii of gyration may be derived, using simultaneous equations. Inertial properties of the body segments of four Pan troglodytes and a single Pongo were determined, and contrasted to comparable findings for humans. Using our approach, the mass distribution characteristics of any individual or species may be represented by a rigid-link segment model or “android.” If this is made to move according to motion functions derived from a real performance of the individual represented, we show that recordings of resulting ground reaction forces may be quite closely simulated by predictive dynamic modelling. © 1996 Wiley-Liss, Inc.