Prediction of biomechanical parameters in the lumbar spine during static sagittal plane lifting.

A combined approach involving optimization and the finite element technique was used to predict biomechanical parameters in the lumbar spine during static lifting in the sagittal plane. Forces in muscle fascicles of the lumbar region were first predicted using an optimization-based force model including the entire lumbar spine. These muscle forces as well as the distributed upper body weight and the lifted load were then applied to a three-dimensional finite element model of the thoracolumbar spine and rib cage to predict deformation, the intradiskal pressure, strains, stresses, and load transfer paths in the spine. The predicted intradiskal pressures in the L3-4 disk at the most deviated from the in vivo measurements by 8.2 percent for the four lifting cases analyzed. The lumbosacral joint flexed, while the other lumbar joints extended for all of the four lifting cases studied (rotation of a joint is the relative rotation between its two vertebral bodies). High stresses were predicted in the posterolateral regions of the endplates and at the junctions of the pedicles and vertebral bodies. High interlaminar shear stresses were found in the posterolateral regions of the lumbar disks. While the facet joints of the upper two lumbar segments did not transmit any load, the facet joints of the lower two lumbar segments experienced significant loads. The ligaments of all lumbar motion segments except the lumbosacral junction provided only marginal moments. The limitations of the current model and possible improvements are discussed.     

Effects of ankle position during static stretching for the hamstrings on the decrease in passive stiffness

Static stretching is frequently performed to improve flexibility of the hamstrings, although the ankle position during hamstring stretching has not been fully investigated. We investigated the effects of ankle position during hamstring stretching on the decrease in passive stiffness. Fourteen healthy men performed static stretching for the hamstrings with the ankle dorsiflexed and plantar-flexed in a randomized order on different days. The hip was passively flexed to the maximum angle which could be tolerated without stretch pain with the knee fully extended; this was maintained for 5 min, with 1-min stretching performed in 5 sessions. Final angles and passive stiffness were measured before and after stretching. The final angle was defined as that formed by the tibia and horizontal plane when the knee was passively extended from hip and knee angles at 90° flexion to the maximum extension angle which could be tolerated without stretch pain. Passive stiffness was determined by the slope of torque–angle curve during the measurement of the final angle. The final angle significantly increased after stretching with the ankle dorsiflexed and plantar-flexed, whereas passive stiffness significantly decreased only after stretching with the ankle planter-flexed. The results suggest that passive stiffness decreases after stretching with the ankle planter-flexed but not after stretching with the ankle dorsiflexed, although the range of joint motion increases regardless of the ankle position during 5-min stretching for the hamstrings. These results indicate that static stretching should be performed with the ankle plantar-flexed when aiming to decrease passive stiffness of the hamstrings.     

Biarticular hip extensor and knee flexor muscle moment arms of the feline hindlimb

Moment arms are important for understanding muscular behavior and for calculating internal muscle forces in musculoskeletal simulations. Biarticular muscles cross two joints and have moment arms that depend on the angle of both joints the muscles cross. The tendon excursion method was used to measure the joint angle-dependence of hamstring (biceps femoris, semimembranosus and semitendinosus) moment arm magnitudes of the feline hindlimb at the knee and hip joints. Knee angle influenced hamstring moment arm magnitudes at the hip joint; compared to a flexed knee joint, the moment arm for semimembranosus posterior at the hip was at most 7.4 mm (25%) larger when the knee was extended. On average, hamstring moment arms at the hip increased by 4.9 mm when the knee was more extended. In contrast, moment arm magnitudes at the knee varied by less than 2.8 mm (mean=1.6 mm) for all hamstring muscles at the two hip joint angles tested. Thus, hamstring moment arms at the hip were dependent on knee position, while hamstring moment arms at the knee were not as strongly associated with relative hip position. Additionally, the feline hamstring muscle group had a larger mechanical advantage at the hip than at the knee joint.     

An improved inverse dynamics formulation for estimation of external and internal loads during human sagittal plane movements

Planar musculoskeletal models are common in the inverse dynamics analysis of human movements such as walking, running and jumping. The continued interest in such models is justified by their simplicity and computational efficiency. Related to a human planar model, a unified formulation for both the flying and support phases of the sagittal plane movements is developed. The actuation involves muscle forces in the lower limbs and the resultant muscle torques in the other body joints. The dynamic equations, introduced in absolute coordinates of the segments, are converted into useful compact forms using the projective technique. The solution to a determinate inverse dynamics problem allows for the explicit determination of the external reactions (presumed to vanish during the flying phases) and the resultant muscle torques in all the model joints. The indeterminate inverse dynamics problem is then focused on the assessment of muscle forces and joint reaction forces selectively in the supporting lower limb. Numerical results of the inverse dynamics simulation of sample sagittal plane movements are reported to illustrate the validity and effectiveness of the improved formulation.     

Brief communication: Cineradiographic analysis of the chimpanzee (Pan troglodytes) talonavicular and calcaneocuboid joints

During terrestrial locomotion, chimpanzees exhibit dorsiflexion of the midfoot between midstance and toe‐off of stance phase, a phenomenon that has been called the “midtarsal break.” This motion is generally absent during human bipedalism, and in chimpanzees is associated with more mobile foot joints than in humans. However, the contribution of individual foot joints to overall foot mobility in chimpanzees is poorly understood, particularly on the medial side of the foot. The talonavicular (TN) and calcaneocuboid (CC) joints have both been suggested to contribute significantly to midfoot mobility and to the midtarsal break in chimpanzees. To evaluate the relative magnitude of motion that can occur at these joints, we tracked skeletal motion of the hindfoot and midfoot during passive plantarflexion and dorsiflexion manipulations using cineradiography. The sagittal plane range of motion was 38 ± 10° at the TN joint and 14 ± 8° at the CC joint. This finding indicates that the TN joint is more mobile than the CC joint during ankle plantarflexion–dorsiflexion. We suggest that the larger range of motion at the TN joint during dorsiflexion is associated with a rotation (inversion–eversion) across the transverse tarsal joint, which may occur in addition to sagittal plane motion. Am J Phys Anthropol 154:604–608, 2014. © 2014 Wiley Periodicals, Inc.     

An improved multi-joint EMG-assisted optimization approach to estimate joint and muscle forces in a musculoskeletal model of the lumbar spine

Muscle force partitioning methods and musculoskeletal system simplifications are key modeling issues that can alter outcomes, and thus change conclusions and recommendations addressed to health and safety professionals. A critical modeling concern is the use of single-joint equilibrium to estimate muscle forces and joint loads in a multi-joint system, an unjustified simplification made by most lumbar spine biomechanical models. In the context of common occupational tasks, an EMG-assisted optimization method (EMGAO) is modified in this study to simultaneously account for the equilibrium at all lumbar joints (M-EMGAO). The results of this improved approach were compared to those of its conventional single-joint equivalent (S-EMGAO) counterpart, the latter method being applied to the same lumbar joints but one at a time. Despite identical geometrical configurations and passive contributions used in both models, computed outcomes clearly differed between single- and multi-joint methods, especially at larger trunk flexed postures and during asymmetric lifting. Moreover, muscle forces predicted by L5-S1 single-joint analyses do not maintain mechanical equilibrium at other spine joints crossed by the same muscles. Assuming that the central nervous system does not attempt to balance the external moments one joint at a time and that a given muscle cannot exert different forces at different joints, the proposed multi-joint method represents a substantial improvement over its single-joint counterpart. This improved approach, hence, resolves trunk muscle forces with biological integrity but without compromising mechanical equilibrium at the lumbar joints.     

Relationship between hip and knee strength and knee valgus during a single leg squat

The purpose of this study was to determine the relationship between hip and knee strength, and valgus knee motion during a single leg squat. Thirty healthy adults (15 men, 15 women) stood on their preferred foot, squatted to approximately 60 deg of knee flexion, and returned to the standing position. Frontal plane knee motion was evaluated using 3-D motion analysis. During Session 2, isokinetic (60 deg/sec) concentric and eccentric hip (abduction/adduction, flexion/extension, and internal/external rotation) and knee (flexion/extension) strength was evaluated. The results demonstrated that hip abduction (r2=0.13), knee flexion (r2=0.18), and knee extension (r2=0.14) peak torque were significant predictors of frontal plane knee motion. Significant negative correlations showed that individuals with greater hip abduction (r=-0.37), knee flexion (r=-0.43), and knee extension (r=-0.37) peak torque exhibited less motion toward the valgus direction. Men exhibited significantly greater absolute peak torque for all motions, excluding eccentric internal rotation. When normalized to body mass, men demonstrated significantly greater strength than women for concentric hip adduction and flexion, knee flexion and extension, and eccentric hip extension. The major findings demonstrate a significant role of hip muscle strength in the control of frontal plane knee motion.     

The effect of hamstring muscle compensation for anterior laxity in the ACL-deficient knee during gait

The hamstring muscles have been recognized as an important element in compensating for the loss of stability in the ACL-deficient knee, but it is still not clear whether the hamstring muscle force can completely compensate for the loss of ACL, and the consequences of increased hamstring muscle force. A two-dimensional anatomical knee model in the sagittal plane was developed to examine the effect of various levels of hamstring muscle activation on restraining anterior tibial translation in the ACL-deficient knee during level walking. The model included the tibiofemoral and patellofemoral joints, four major ligaments, the medial capsule, and five muscle units surrounding the knee. Simulations were conducted to determine anterior tibial translation and internal joint loading at a single selected position when the knee was under a peak external flexion moment during early stance phase of gait. Incremental hamstring muscle forces were applied to the modeled normal and the ACL-deficient knees. Results of simulations showed that the ACL injury increased the anterior tibial translation by 11.8mm, while 56% of the maximal hamstring muscle force could reduce the anterior translation of the tibia to a normal level during the stance phase of gait. The consequences of increased hamstring muscle force included increased quadriceps muscle force and joint contact force.     

Paired nonlinear behavior of active and passive joint torques associated with preparation for walk-to-run gait transition

This study investigated the lower extremity torque's active and passive features during the walk-to-run gait transition with continuously increased walking speed. Fourteen volunteers participated in the experiment. Kinematic and kinetic data were collected synchronously. Five strides leading up the gait transition were examined. Peaks of the passive (e.g., contact) and active (e.g., generalized muscle torques), along with net joint torque, and time to peak torques exhibited significant differences at the last stride before gait transition, compared to the first four strides, at the ankle, knee, and hip joints, respectively. Selected peak joint active and passive torques showed significant and opposite trends at critical events within a stride cycle: such ankle joint right after heel-contact, knee joint during weight acceptance, and both hip and knee joints right before toe-off. The magnitude and the corresponding time to active and passive peak torque changed in a nonlinear pattern before the transition from walk to run. The lower extremity segment-interaction during gait transition appeared to be an active reorganization exemplified by the interaction between the lower extremity's active and passive torque components.     

Simple and complex models for studying muscle function in walking

While simple models can be helpful in identifying basic features of muscle function, more complex models are needed to discern the functional roles of specific muscles in movement. In this paper, two very different models of walking, one simple and one complex, are used to study how muscle forces, gravitational forces and centrifugal forces (i.e. forces arising from motion of the joints) combine to produce the pattern of force exerted on the ground. Both the simple model and the complex one predict that muscles contribute significantly to the ground force pattern generated in walking; indeed, both models show that muscle action is responsible for the appearance of the two peaks in the vertical force. The simple model, an inverted double pendulum, suggests further that the first and second peaks are due to net extensor muscle moments exerted about the knee and ankle, respectively. Analyses based on a much more complex, muscle-actuated simulation of walking are in general agreement with these results; however, the more detailed model also reveals that both the hip extensor and hip abductor muscles contribute significantly to vertical motion of the centre of mass, and therefore to the appearance of the first peak in the vertical ground force, in early single-leg stance. This discrepancy in the model predictions is most probably explained by the difference in model complexity. First, movements of the upper body in the sagittal plane are not represented properly in the double-pendulum model, which may explain the anomalous result obtained for the contribution of a hip-extensor torque to the vertical ground force. Second, the double-pendulum model incorporates only three of the six major elements of walking, whereas the complex model is fully 3D and incorporates all six gait determinants. In particular, pelvic list occurs primarily in the frontal plane, so there is the potential for this mechanism to contribute significantly to the vertical ground force, especially during early single-leg stance when the hip abductors are activated with considerable force.     

Test-Retest Reliability of Innovated Strength Tests for Hip Muscles

The burden of hip muscles weakness and its relation to other impairments has been well documented. It is therefore a pre-requisite to have a reliable method for clinical assessment of hip muscles function allowing the design and implementation of a proper strengthening program. Motor-driven dynamometry has been widely accepted as the gold-standard for lower limb muscle strength assessment but is mainly related to the knee joint. Studies focusing on the hip joint are less exhaustive and somewhat discrepant with regard to optimal participants position, consequently influencing outcome measures. Thus, we aimed to develop a standardized test setup for the assessment of hip muscles strength, i.e. flexors/extensors and abductors/adductors, with improved participant stability and to define its psychometric characteristics. Eighteen participants performed unilateral isokinetic and isometric contractions of the hip muscles in the sagittal and coronal plane at two separate occasions. Peak torque and normalized peak torque were measured for each contraction. Relative and absolute measures of reliability were calculated using the intraclass correlation coefficient and standard error of measurement, respectively. Results from this study revealed higher levels of between-day reliability of isokinetic/isometric hip abduction/flexion peak torque compared to existing literature. The least reliable measures were found for hip extension and adduction, which could be explained by a less efficient stabilization technique. Our study additionally provided a first set of reference normalized data which can be used in future research.     

Effect of hip flexibility on optimal stalder performances on high bar

In the optimisation of sports movements using computer simulation models, the joint actuators must be constrained in order to obtain realistic results. In models of a gymnast, the main constraint used in previous studies was maximum voluntary active joint torque. In the stalder, gymnasts reach their maximal hip flexion under the bar. The purpose of this study was to introduce a model of passive torque to assess the effect of the gymnast's flexibility on the technique of the straddled stalder. A three-dimensional kinematics driven simulation model was developed. The kinematics of the shoulder flexion, hip flexion and hip abduction were optimised to minimise torques for four hip flexion flexibilities: 100°, 110°, 120° and 130°. With decreased flexibility, the piked posture period is shorter and occurs later. Moreover the peaks of shoulder and hip torques increase. Gymnasts with low hip flexibility need to be stronger to achieve a stalder; hip flexibility should be considered by coaches before teaching this skill.     

Role of heel lifting in standing balance recovery: A simulation study

Although lifting the heels has frequently been observed during balance recovery, the function of this movement has generally been overlooked. The present study aimed to investigate the functional role of heel lifting during regaining balance from a perturbed state. Computer simulation was employed to objectively examine the effect of allowing/constraining heel lifting on balance performance. The human model consisted of 3 rigid body segments connected by frictionless joints. Movements were driven by joint torques depending on current joint angle, angular velocity, and activation level. Starting from forward-inclined and static straight-body postures, the optimization goal was to recover balance effectively (so that ground projection of the mass center returned to the inside of the base of support) and efficiently by adjusting ankle and hip joint activation levels. Allowing/constraining heel lifting resulted in virtually identical movements when balance was mildly perturbed at the smallest lean angle (8°). At larger lean angles (8.5° and 9°), heel lifting assisted balance recovery more evidently with larger joint movements. Partial and altered timings of ankle/hip torque activation due to constraining heel lifting reduced linear and angular momentum generation for avoiding forward falling, and resulted in hindered balancing performance.     

Intervertebral reaction force prediction using an enhanced assembly of OpenSim models

OpenSim offers a valuable approach to investigating otherwise difficult to assess yet important biomechanical parameters such as joint reaction forces. Although the range of available models in the public repository is continually increasing, there currently exists no OpenSim model for the computation of intervertebral joint reactions during flexion and lifting tasks. The current work combines and improves elements of existing models to develop an enhanced model of the upper body and lumbar spine. Models of the upper body with extremities, neck and head were combined with an improved version of a lumbar spine from the model repository. Translational motion was enabled for each lumbar vertebrae with six controllable degrees of freedom. Motion segment stiffness was implemented at lumbar levels and mass properties were assigned throughout the model. Moreover, body coordinate frames of the spine were modified to allow straightforward variation of sagittal alignment and to simplify interpretation of results. Evaluation of model predictions for level L1–L2, L3–L4 and L4–L5 in various postures of forward flexion and moderate lifting (8 kg) revealed an agreement within 10% to experimental studies and model-based computational analyses. However, in an extended posture or during lifting of heavier loads (20 kg), computed joint reactions differed substantially from reported in vivo measures using instrumented implants. We conclude that agreement between the model and available experimental data was good in view of limitations of both the model and the validation datasets. The presented model is useful in that it permits computation of realistic lumbar spine joint reaction forces during flexion and moderate lifting tasks. The model and corresponding documentation are now available in the online OpenSim repository.     

Comparison of erector spinae and hamstring muscle activities and lumbar motion during standing knee flexion in subjects with and without lumbar extension rotation syndrome

The aim of this study was to compare the activity of the erector spinae (ES) and hamstring muscles and the amount and onset of lumbar motion during standing knee flexion between individuals with and without lumbar extension rotation syndrome. Sixteen subjects with lumbar extension rotation syndrome (10 males, 6 females) and 14 healthy subjects (8 males, 6 females) participated in this study. During the standing knee flexion, surface electromyography (EMG) was used to measure muscle activity, and surface EMG electrodes were attached to both the ES and hamstring (medial and lateral) muscles. A three-dimensional motion analysis system was used to measure kinematic data of the lumbar spine. An independent-t test was conducted for the statistical analysis. The group suffering from lumbar extension rotation syndrome exhibited asymmetric muscle activation of the ES and decreased hamstring activity. Additionally, the group with lumbar extension rotation syndrome showed greater and earlier lumbar extension and rotation during standing knee flexion compared to the control group. These data suggest that asymmetric ES muscle activation and a greater amount of and earlier lumbar motion in the sagittal and transverse plane during standing knee flexion may be an important factor contributing to low back pain.     

A computational model for dynamic analysis of the human gait

Biomechanical models are important tools in the study of human motion. This work proposes a computational model to analyse the dynamics of lower limb motion using a kinematic chain to represent the body segments and rotational joints linked by viscoelastic elements. The model uses anthropometric parameters, ground reaction forces and joint Cardan angles from subjects to analyse lower limb motion during the gait. The model allows evaluating these data in each body plane. Six healthy subjects walked on a treadmill to record the kinematic and kinetic data. In addition, anthropometric parameters were recorded to construct the model. The viscoelastic parameter values were fitted for the model joints (hip, knee and ankle). The proposed model demonstrated that manipulating the viscoelastic parameters between the body segments could fit the amplitudes and frequencies of motion. The data collected in this work have viscoelastic parameter values that follow a normal distribution, indicating that these values are directly related to the gait pattern. To validate the model, we used the values of the joint angles to perform a comparison between the model results and previously published data. The model results show a same pattern and range of values found in the literature for the human gait motion.     

Lumbar posterior ligament involvement during extremely heavy lifts estimated from fluoroscopic measurements.

The mechanical role of the lumbar posterior ligaments during lifting tasks remains controversial. This study was designed to assess the ligament and disc contribution in resisting trunk flexion moment during extremely heavy lifts performed by national class powerlifters. Direct measurements of lumbar vertebrae kinematics in sagittal plane were obtained from videofluoroscopy utilizing multiple digitizing, correction for optical distortions and digital filtering. Four experienced powerlifters executed three trials, resulting in about 72 mA s of total radiation exposure. In the first trial, joint angles were measured when subjects fully flexed their spines to a point where the passive tissues resisted the flexor moment creating myoelectric silence in the extensor musculature. Next, two conventional deadlift style lifts were executed with the barbell load ranging from 183.7 to 210.9 kg. Four vertebral corners were digitized at a sampling rate of 30 Hz. The relative intervertebral joint angles, distance between the ligament attachment points, shearing and compressive displacements were calculated from a rigid body motion approach. Analysis revealed that except for one trial of one subject, they accomplished their lifts with an amount of lumbar flexion between 1.5 and 13 degrees less than they demonstrated during full flexion. Resultant ligament lengths at the beginning of the lifts ranged from 56.1 to 99.8% of their lengths when the trunk was fully flexed. It was concluded that ligaments did not strain sufficiently to contribute substantial resistance to the trunk flexion moment, relegating this responsibility to the musculature.     

A musculoskeletal model of the human lower extremity: the effect of muscle, tendon, and moment arm on the moment-angle relationship of musculotendon actuators at the hip, knee, and ankle

We have developed a musculoskeletal model of the human lower extremity for computer simulation studies of musculotendon function and muscle coordination during movement. This model incorporates the salient features of muscle and tendon, specifies the musculoskeletal geometry and musculotendon parameters of 18 musculotendon actuators, and defines the active isometric moment of these actuators about the hip, knee, and ankle joints in the sagittal plane. We found that tendon slack length, optimal muscle-fiber length, and moment arm are different for each actuator, thus each actuator develops peak isometric moment at a different joint angle. The joint angle where an actuator produces peak moment does not necessarily coincide with the joint angle where: (1) muscle force peaks, (2) moment arm peaks, or (3) the in vivo moment developed by maximum voluntary contractions peaks. We conclude that when tendon is neglected in analyses of musculotendon force or moment about joints, erroneous predictions of human musculotendon function may be stated, not only in static situations as studied here, but during movement as well.     

A three-dimensional whole-body model to predict human walking on level ground

Predictive simulation of human walking has great potential in clinical motion analysis and rehabilitation engineering assessment, but large computational cost and reliance on measurement data to provide initial guess have limited its wide use. We developed a computationally efficient model combining optimization and inverse dynamics to predict three-dimensional whole-body motions and forces during human walking without relying on measurement data. Using the model, we explored two different optimization objectives, mechanical energy expenditure and the time integral of normalized joint torque. Of the two criteria, the sum of the time integrals of the normalized joint torques produced a more realistic walking gait. The reason for this difference is that most of the mechanical energy expenditure is in the sagittal plane (based on measurement data) and this leads to difficulty in prediction in the other two planes. We conclude that mechanical energy may only account for part of the complex performance criteria driving human walking in three dimensions.

     

Lower extremity biomechanics during a regular and counterbalanced squat

If the efficiency of human movement patterns could be improved using exercise, this could lead to more effective musculoskeletal disease-injury prevention and rehabilitation programs. It has been suggested that an efficient squat movement pattern emphasizes the use of the large hip extensors instead of the smaller knee extensors. The purpose of this study was to determine whether a counterbalanced squat (CBS) could produce a more hip-dominant and less knee-dominant squat movement pattern as compared with a regular squat (RS). There were 31 recreationally trained college-aged participants (15 male, 16 female) who performed 10 squats (5 CBS and 5 RS), while segment kinematics, ground reaction forces, and muscle (gluteus maximus [GM], quadriceps, hamstrings) electromyographic (EMG) activations were recorded. Peak sagittal plane net joint moments and joint ranges of motion at the hip, knee, and ankle joints along with peak and integrated EMG activation levels for all 3 muscles were compared using analysis of variance (squat type × sex). The results revealed that the CBS increased the hip joint moment and GM activation, while it decreased the knee joint moment and quadriceps activation as compared with the RS. Therefore, the CBS produces a more hip-dominant and less knee-dominant squat movement pattern and could be used in exercise programs aimed at producing more hip-dominant movement patterns.