Lower extremity control during turns initiated with and without hip external rotation

The pirouette turn is often initiated in neutral and externally rotated hip positions by dancers. This provides an opportunity to investigate how dancers satisfy the same mechanical objectives at the whole-body level when using different leg kinematics. The purpose of this study was to compare lower extremity control strategies during the turn initiation phase of pirouettes performed with and without hip external rotation. Skilled dancers (n=5) performed pirouette turns with and without hip external rotation. Joint kinetics during turn initiation were determined for both legs using ground reaction forces (GRFs) and segment kinematics. Hip muscle activations were monitored using electromyography. Using probability-based statistical methods, variables were compared across turn conditions as a group and within-dancer. Despite differences in GRFs and impulse generation between turn conditions, at least 90% of each GRF was aligned with the respective leg plane. A majority of the net joint moments at the ankle, knee, and hip acted about an axis perpendicular to the leg plane. However, differences in shank alignment relative to the leg plane affected the distribution of the knee net joint moment when represented with respect to the shank versus the thigh. During the initiation of both turns, most participants used ankle plantar flexor moments, knee extensor moments, flexor and abductor moments at the push leg׳s hip, and extensor and abductor moments at the turn leg׳s hip. Representation of joint kinetics using multiple reference systems assisted in understanding control priorities.     

Changes in axial stiffness of the trunk as a function of walking speed

Research suggests that abnormal coordination patterns between the thorax and pelvis in the transverse plane observed in patients with Parkinson's disease and the elderly might be due to alteration in axial trunk stiffness. The purpose of this study was to develop a tool to estimate axial trunk stiffness during walking and to investigate its functional role. Fourteen healthy young subjects participated in this study. They were instructed to walk on the treadmill and kinematic data was collected by 3D motion analysis system. Axial trunk stiffness was estimated from the angular displacement between trunk segments and the amount of torque around vertical axis of rotation. The torque due to arm swing cancelled out the torque due to the axial trunk stiffness during walking and the thoracic rotation was of low amplitude independent of changes in walking speeds within the range used in this study (0.85-1.52 m/s). Estimated axial trunk stiffness increased with increasing walking speed. Functionally, the suppression of axial rotation of thorax may have a positive influence on head stability as well as allowing recoil between trunk segments. Furthermore, the increased stiffness at increased walking speed would facilitate the higher frequency rotation of the trunk in the transverse plane required at the higher walking speeds.     

Individuals with non-specific low back pain use a trunk stiffening strategy to maintain upright posture

There is increasing evidence that individuals with non-specific low back pain (LBP) have altered movement coordination. However, the relationship of this neuromotor impairment to recurrent pain episodes is unknown. To assess coordination while minimizing the confounding influences of pain we characterized automatic postural responses to multi-directional support surface translations in individuals with a history of LBP who were not in an active episode of their pain. Twenty subjects with and 21 subjects without non-specific LBP stood on a platform that was translated unexpectedly in 12 directions. Net joint torques of the ankles, knees, hips, and trunk in the frontal and sagittal planes as well as surface electromyographs of 12 lower leg and trunk muscles were compared across perturbation directions to determine if individuals with LBP responded using a trunk stiffening strategy. Individuals with LBP demonstrated reduced peak trunk torques, and enhanced activation of the trunk and ankle muscle responses following perturbations. These results suggest that individuals with LBP use a strategy of trunk stiffening achieved through co-activation of trunk musculature, aided by enhanced distal responses, to respond to unexpected support surface perturbations. Notably, these neuromotor alterations persisted between active pain periods and could represent either movement patterns that have developed in response to pain or could reflect underlying impairments that may contribute to recurrent episodes of LBP.     

The three-dimensional kinematics and flexibility characteristics of the human ankle and subtalar joints--Part I: Kinematics

The in-vitro, three dimensional kinematic characteristics of the human ankle and subtalar joint were investigated in this study. The main goals of this investigation were: 1) To determine the range of motion of the foot-shank complex and the associated range of motion of the ankle and subtalar joints; 2) To determine the kinematic coupling characteristics of the foot-shank complex, and 3) To identify the relationship between movements at the ankle and subtalar joints and the resulting motion produced between the foot and the shank. The tests were conducted on fifteen fresh amputated lower limbs and consisted of incrementally displacing the foot with respect to the shank while the motion of the articulating bones was measured through a three dimensional position data acquisition system. The kinematic analysis was based on the helical axis parameters describing the incremental displacements between any two of the three articulating bones and on a joint coordinate system used to describe the relative position between the bones. From the results of this investigation it was concluded that: 1) The range of motion of the foot-shank complex in any direction (dorsiflexion/plantarflexion, inversion/eversion and internal rotation/external rotation) is larger than that of either the ankle joint or the subtalar joint.; 2) Large kinematic coupling values are present at the foot-shank complex in inversion/eversion and in internal rotation/external rotation. However, only a slight amount of coupling was observed to occur in dorsiflexion/plantarflexion.; 3) Neither the ankle joint nor the subtalar joint are acting as ideal hinge joints with a fixed axis of rotation.; 4) Motion of the foot-shank complex in any direction is the result of rotations at both the ankle and the subtalar joints. However, the contribution of the ankle joint to dorsiflexion/plantarflexion of the foot-shank complex is larger than that of the subtalar joint and the contribution of the subtalar joint to inversion/eversion is larger than that of the ankle joint.; 5) The ankle and the subtalar joints have an approximately equal contribution to internal rotation/external rotation movements of the foot-shank complex.     

Segment and joint angles of hind limb during bipedal and quadrupedal walking of the bonobo (Pan paniscus)

We describe segment angles (trunk, thigh, shank, and foot) and joint angles (hip, knee, and ankle) for the hind limbs of bonobos walking bipedally ("bent-hip bent-knee walking," 17 sequences) and quadrupedally (33 sequences). Data were based on video recordings (50 Hz) of nine subjects in a lateral view, walking at voluntary speed. The major differences between bipedal and quadrupedal walking are found in the trunk, thigh, and hip angles. During bipedal walking, the trunk is approximately 33-41 degrees more erect than during quadrupedal locomotion, although it is considerably more bent forward than in normal human locomotion. Moreover, during bipedal walking, the hip has a smaller range of motion (by 12 degrees ) and is more extended (by 20-35 degrees ) than during quadrupedal walking. In general, angle profiles in bonobos are much more variable than in humans. Intralimb phase relationships of subsequent joint angles show that hip-knee coordination is similar for bipedal and quadrupedal walking, and resembles the human pattern. The coordination between knee and ankle differs much more from the human pattern. Based on joint angles observed throughout stance phase and on the estimation of functional leg length, an efficient inverted pendulum mechanism is not expected in bonobos.     

Mobility of the subtalar joint in the intact ankle complex

A previous study by these authors showed that the calcaneus follows a unique path of unresisted coupled motion relative to the tibia during passive flexion and that most of this motion occurred at the ankle level. Subtalar motion in the intact ankle complex was observed only when perturbations from this path were induced by the application of force to the calcaneus. Relative motion of the bones of the ankle complex was tracked by stereophotogrammetry in seven specimens. Anatomical landmarks, reference frames and joint angles were defined by standard techniques. Sequential moments were applied to the calcaneus about the long axis of the tibia. Measured movements at subtalar level demonstrated plantarflexion coupled to supination and internal rotation (inversion) and dorsiflexion coupled to pronation and external rotation (eversion). These movements were resisted and were fully recovered when the external load was removed. Subtalar motion diminished as the ankle approached maximal dorsi- and plantarflexion. Two clearly distinguished mean axes of rotation were observed for inversion and eversion runs. The axes of inversion and eversion of the subtalar complex changed orientation along a preferred and repeatable path. The subtalar joint complex occupied only a single stable position in the unloaded state and with no range of unresisted motion. It is inferred that mobility was possible only by the stretching and lengthening of the ligaments and the indentation of the articular surfaces, requiring the application of loads. The subtalar joint complex behaves like a flexible structure.     

Reduced knee joint moment in ACL deficient patients at a cost of dynamic stability during landing

The current study aimed to examine the effect of anterior cruciate ligament deficiency (ACLd) on joint kinetics and dynamic stability control after a single leg hop test (SLHT). Twelve unilateral ACLd patients and a control subject group (n=13) performed a SLHT over a given distance with both legs. The calculation of joint kinetics was done by means of a soft-tissue artifact optimized rigid full-body model. Margin of stability (MoS) was quantified by the difference between the base of support and the extrapolated center of mass. During landing, the ACLd leg showed lower external knee flexion moments but demonstrated higher moments at the ankle and hip compared to controls (p<0.05). The main reason for the joint moment redistribution in the ACLd leg was a more anterior position of the ground reaction force (GRF) vector, which affected the moment arms of the GRF acting about the joints (p<0.05). For the ACLd leg, trunk angle was more flexed over the entire landing phase compared to controls (p<0.05) and we found a significant correlation between moment arms at the knee joint and trunk angle (r2 = 0.48;p<0.01). The consequence of this altered landing strategy in ACLd legs was a more anterior position of the center of mass reducing the MoS (p<0.05). The results illustrate the interaction between trunk angle, joint kinetics and dynamic stability during landing maneuvers and provide evidence of a feedforward adaptive adjustment in ACLd patients (i.e. more flexed trunk angle) aimed at reducing knee joint moments at the cost of dynamic stability control.     

Inverse dynamic investigation of voluntary leg lateral movements in weightlessness: a new microgravity-specific strategy

This study deals with the quantitative assessment of exchanged forces and torques at the restraint point during whole body posture perturbation movements in long-term microgravity. The work was based on the results of a previous study focused on trunk bending protocol, which suggested that the minimization of the torques exchanged at the restraint point could be a strategy for movement planning in microgravity (J. Biomech. 36(11) (2003) 1691). Torques minimization would lead to the optimization of muscles activity, to the minimization of energy expenditure and, ultimately, to higher movement control capabilities. Here, we focus on leg lateral abduction from anchored stance. The analysis was based on inverse dynamic modelling, leading to the estimation of the total angular momentum at the supporting ankle joint. Results agree with those obtained for trunk bending movements and point out a consistent minimization of the torques exchanged at the restraint point in weightlessness. Given the kinematic features of the examined motor task, this strategy was interpreted as a way to master the rotational dynamic effects on the frontal plane produced by leg lateral abduction. This postural stabilizing effects was the result of a multi-segmental compensation strategy, consisting of the counter rotation of the supporting limb and trunk accompanying the leg raising. The observed consistency of movement-posture co-ordination patterns among lateral leg raising and trunk bending is put forward as a novel interpretative issue of the adaptation mechanisms of the motor system to sustained microgravity, especially if one considers the completely different kinematics of the centre of mass, which was observed in weightlessness for these two motor tasks.     

Articulated external fixation of the ankle: minimizing motion resistance by accurate axis alignment

This study describes how an optimal single hinge axis position can be established for the application of articulated external fixation to the ankle joint. By deliberately introducing various amounts of relative mal-alignment between the optimal talocrural joint axis and the actual fixator hinge axis, it was possible to measure the corresponding amounts of additional resistance to joint motion. In a cadaveric study of six ankle specimens, we determined the instant axis of rotation of the talocrural joint from 3-D kinematic data. acquired by an electromagnetic motion tracking system. For each specimen, an optimal fixator hinge position was calculated from these motion data. Compared to the intact natural joint, aligning the fixator along the optimized axis position caused a moderate increase in energy (0.14 J) needed to rotate the ankle through a prescribed plantar/dorsiflexion range. However, malpositioning the hinge by 10 mm caused more than five times that amount of increase in motion resistance. While articulated external fixation with limited internal fixation can establish a favorable environment for the repair of severe injuries such as tibial pilon fractures, the large additional resistance to motion accompanying a malpositioned fixator axis suggests the development of untoward intra-articular forces that could act to disturb fragment alignment.     

Intralimb compensation strategy depends on the nature of joint perturbation in human hopping

Due to the well-described spring-mass dynamics of bouncing gaits, human hopping is a tractable model for elucidating basic neuromuscular compensation principles. We tested whether subjects would employ a multi-joint or single-joint response to stabilize leg stiffness while wearing a spring-loaded ankle-foot orthosis (AFO) that applied localized resistive and assistive torques to the ankle. We analyzed kinematics and kinetics data from nine subjects hopping in place on one leg, at three frequencies (2.2, 2.4, and 2.8Hz) and three orthosis conditions (freely articulating AFO, AFO with plantarflexion resistance, and AFO with plantarflexion assistance). Leg stiffness was invariant across AFO conditions, however, compensation strategy depended upon the nature of the applied load. Biological ankle stiffness increased in response to a resistive load at twice the rate that it decreased with an assitive load. Ankle adjustments alone fully compensated for an assistive load with no net change in combined (biological plus applied) total ankle stiffness (p > or =0.133). In contrast, a resistive load resulted in a 7.4-9.0% increase in total ankle stiffness across frequencies and a concomitant 10-15% increase in knee joint stiffness at each frequency (p< or =0.037). The increased knee joint stiffness in response to resistive ankle load allowed subjects to maintain a more flexed knee at mid-stance, which attenuated the effect of the increased total ankle joint stiffness to preserve leg stiffness and whole limb biomechanical performance. Our findings suggest humans maintain invariant leg stiffness in bouncing gaits through different intralimb compensation strategies that are specific to the nature of the joint loading.     

Estimates of Achilles tendon moment arm differ when axis of ankle rotation is derived from ankle motion

The plantarflexor moment arm of the Achilles tendon determines the mechanical advantage of the triceps surae and also indirectly affects muscle force generation by setting the amount of muscle-tendon shortening per unit of ankle joint rotation. The Achilles tendon moment arm may be determined geometrically from an axis (or center) of joint rotation and the line of action of the tendon force, but such moment arms may be sensitive to the location of the joint axis. Using motion analysis to track an ultrasound probe overlying the Achilles tendon along with markers on the shank and foot, we measured Achilles tendon moment arm during loaded and unloaded dynamic plantarflexion motions in 15 healthy subjects. Three representations of the axis or center of rotation of the ankle were considered: (1) a functional axis, defined by motions of the foot and shank; (2) a transmalleolar axis; and (3) a transmalleolar midpoint. Moment arms about the functional axis were larger than those found using the transmalleolar axis and transmalleolar midpoint (all p < 0.001). Moment arms computed with the functional axis increased with plantarflexion angle (all p < 0.001), and increased with loading in the most plantarflexed position (p < 0.001) but these patterns were not observed when either using a transmalleolar axis or transmalleolar midpoint. Functional axis moment arms were similar to those estimated previously using magnetic resonance imaging, suggesting that using a functional axis for ultrasound-based geometric estimates of Achilles tendon moment arm is an improvement over landmark-based methods.     

Activation of walking by electrical stimulation in humans under the conditions of muscle unloading and its variations under the effect of afferent influences

The possibility of initiating an involuntary walking rhythm in a suspended human leg by electrical stimulation was studied. The subjects lay on the side with one leg suspended in an exoskeleton allowing horizontal rotation in three joints: the hip, knee, and ankle ones. To evoke involuntary walking of the suspended leg, two methods were used: continuous vibration of the quadriceps muscle of the hip and electrical stimulation of the cutaneous nerves innervating the foot of the immobile leg. The hip and ankle were involved in the involuntary movements, with reciprocal bursts of electromyographic activity being also observed in the antagonistic muscles of the hip. The application of an external load (4 N or 8 N) to the foot caused a perceptible intensification of its movements. An additional weight (0.5 kg) or a rubber band wrapped around the foot caused no substantial change in the pattern of stimulated walking. Electrical stimulation is an effective means of activating walking movements, and their characteristics confirm the assumption that the walking rhythm is of central origin. Additional afferentation from the sole’s receptors plays an important role in the modulation of the induced movements and the modification of the general walking pattern under the conditions of muscle unloading.     

Quantitative measurement of spastic ankle joint stiffness using a manual device: A preliminary study

Quantitative measurement of ankle joint stiffness following stroke could prove useful in monitoring the progress of a rehabilitation programme. The objective of this study was to design a manual device for use in the clinical setting. Manual measurement of spastic ankle joint stiffness has historically been conducted using hand-held dynamometers or alternative devices, but some difficulties have been reported in controlling the velocity applied to the ankle during the measurement. In this study, a manually operated device was constructed with a footplate, a torquemeter and a potentiometer. It was mechanically designed to rotate around an approximated axis of the ankle joint and to measure ankle joint angular position and its corresponding resistive torque. Two stroke hemiplegic subjects pariticapted in a pilot study. The results suggested that difficulty in controlling the applied velocity might be complemented by presenting torque data as a function of peak angular velocity in each stretching cycle. Moreover, the results demonstrated that the device could potentially apply a wide range of angular velocities and provide potentially useful clinical information. Quantitative data successfully acquired using this method included the approximate ankle angular position, where the velocity-dependent characteristics of stiffness was notably initiated and its corresponding torque and velocity.     

A geometric model of the human ankle joint.

A two-dimensional four-bar linkage model of the ankle joint is formulated to describe dorsi/plantarflexion in unloaded conditions as observed in passive tests on ankle complex specimens. The experiments demonstrated that the human ankle joint complex behaves as a single-degree-of-freedom system during passive motion, with a moving axis of rotation. The bulk of the movement occurred at the level of the ankle. Fibres within the calcaneofibular and tibiocalcaneal ligaments remained approximately isometric. The experiments showed that passive kinematics of the ankle complex is governed only by the articular surfaces and the ligaments. It was deduced that the ankle is a single-degree-of-freedom mechanism where mobility is allowed by the sliding of the articular surfaces upon each other and the isometric rotation of two ligaments about their origins and insertions, without tissue deformation. The linkage model is formed by the tibia/fibula and talus/calcaneus bone segments and by the calcaneofibular and tibiocalcaneal ligament segments. The model predicts the path of calcaneus motion, ligament orientations, instantaneous axis of rotation, and conjugate talus surface profile as observed in the experiments. Many features of ankle kinematics such as rolling and multiaxial rotation are elucidated. The geometrical model is a necessary preliminary step to the study of ankle joint stability in response to applied loads and can be used to predict the effects of changes to the original geometry of the intact joint. Careful reconstruction of the original geometry of the ligaments is necessary after injury or during total ankle replacement.     

In vivo determination of the anatomical axes of the ankle joint complex: An optimization approach

This study investigates the feasibility of a subject-specific three-dimensional model of the ankle joint complex for kinematic and dynamic analysis of movement. The ankle joint complex was modelled as a three-segment system, connected by two ideal highe joints: the talocrural and the subtalar joint. A mathematical formulation was developed to express the three-dimensional translation and rotation between the foot and shank segments as a function of the two joint angles, and 12 model parameters describing the locations of the joint axes. An optimization method was used to fit the model parameters to three-dimensional kinematic data of foot and shank markers, obtained during test movements throughout the entire physiological range of motion of the ankle joint. The movement of the talus segment, which cannot be measured non-invasively, is not necessary for the analysis.

This optimization method was used to determine the position and orientation of the joint axes in 14 normal subjects. After optimization, the discrepancy between the best fitting model and actual marker kinematics was between 1 and 3 mm for all subjects. The predicted inclination of the subtalar joint axis from the horizontal plane was 37.4±2.7°, and the medial deviation was 18.0±16.2°. The lateral side of the talucrural axis was directed slightly posteriorly (6.8±8.1°), and inclined downward by 7.0±5.4°. These results are similar to previously reported typical results from anatomical, in vitro, studies. Reproducibility was evaluated by repeated testing of one subject, which resulted in variations of about one-fifth of the standard deviation within the group, the inclination of the subtalar joint axis was significantly correlated to the arch height and a radiographic ‘tarsal index’. It is concluded that this optimization method provides the opportunity to incorporate inter-individual anatomical differences into kinematic and dynamic analysis of the ankle joint complex. This allows a more functional interpretation of kinematic data, and more realistic estimates of internal forces.     

Subject-specific axes of the ankle joint complex

The aim of this study was to use a two-axis ankle joint model and an optimisation process (van den Bogert et al., 1994) to calculate and compare the talocrural and subtalar hinge axes for non-weight-bearing ankle motion, weight-bearing ankle motion, and walking in normal, healthy adult subjects and to see which of the first two sets of axes better fit the walking data. Motion data for the foot and shank were collected on eight subjects whilst they performed the activities mentioned. After choosing the best marker sets for motion tracking, a two-hinge ankle joint model was fit to the motion data. Ankle joint ranges of motion were also calculated. It was found that the model fit the experimental data well, with non-weight-bearing motion achieving the best fit. Despite this, the calculated axis orientations were highly variable both between motion types and between subjects. No significant difference between the fit of the non-weight-bearing and weight-bearing models to the walking data was found, which implies that either set of functional axes is adequate for modeling walking; however, the subtalar deviation angle was significantly closer for the weight-bearing activity and walking than for the non-weight-bearing activity and walking, which suggests that it is marginally better to use the weight-bearing functional motions. The results lead to questions about the appropriateness of the two-hinge ankle model for use in applications in which the behaviour of the individual joints of the ankle complex, rather than simply the relative motion of the leg and foot, is important.     

A kinematic model of the human ankle

The spatial gross motion of the foot with respect to the shank is modelled as rotations about two fixed ankle axes: the upper ankle rotation axis (plantarflexion/dorsiflexion) and the subtalar rotation axis (inversion/eversion). The positions of the axes are determined by externally visible bony landmarks of the lower leg and are measured for a living subject. The model input data are the plantarflexion/dorsiflexion and inversion/eversion rotation angles; the model output is a 4 × 4 transformation matrix which quantitatively describes the relative position of a foot coordinate system with respect to a shank coordinate system.     

The sit-up: complex kinematics and muscle activity in voluntary axial movement.

This paper describes the kinematics and muscle activity associated with the standard sit-up, as a first step in the investigation of complex motor coordination. Eight normal human subjects lay on a force table and performed at least 15 sit-ups, with the arms across the chest and the legs straight and unconstrained. Several subjects also performed sit-ups with an additional weight added to the head. Support surface forces were recorded to calculate the location of the center of pressure and center of gravity; conventional motion analysis was used to measure segmental positions; and surface EMG was recorded from eight muscles. While the sit-up consists of two serial components, 'trunk curling' and 'footward pelvic rotation', it can be further subdivided into five phases, based on the kinematics. Phases I and II comprise trunk curling. Phase I consists of neck and upper trunk flexion, and phase II consists of lumbar trunk lifting. Phase II corresponds to the point of peak muscle contraction and maximum postural instability, the 'critical point' of the sit-up. Phases III-V comprise footward pelvic rotation. Phase III begins with pelvic rotation towards the feet, phase IV with leg lowering, and phase V with contact between the legs and the support surface. The overall pattern of muscle activity was complex with times of EMG onset, peak activity, offset, and duration differing for different muscles. This complex pattern changed qualitatively from one phase to the next, suggesting that the roles of different muscles and, as a consequence, the overall form of coordination, change during the sit-up.     

In vivo determination of the Achilles tendon moment arm in three-dimensions

Two-dimensional methods have been applied to determine the Achilles tendon moment arm in previous studies, although the talocrural joint rotates in three-dimension. The purpose of this study was to develop a method for determining the Achilles tendon moment arm in three-dimensions (3DMA). A series of sagittal ankle images were obtained at ankle positions of -20°, -10° (dorsiflexed position), 0° (neutral position), +10°, +20°, and +30° (plantarflexed position). The talocrural joint axis was determined as the finite helical axis of the ankle joint over 20° of displacement, and the 3DMA was determined as the shortest distance from the talocrural joint axis to the line of action of the Achilles tendon force. The corresponding 2DMA was determined with the center of rotation method using the images captured on the sagittal plane passing through the mid-point of the medio-lateral width of the tibia. The 3DMA ranged from 35 to 41 mm across various ankle positions and was, on average, 11 mm smaller than 2DMA. The difference between the two measures was attributable primarily to the deviations of the talocrural joint axis from the anatomical medio-lateral direction. The deviations on the coronal plane (21.4±20.7°) and on the transverse planes (14.8±22.6°) accounted for the errors of 1.3 mm and 3.0 mm, respectively. In addition, selecting either a medially or laterally misaligned sagittal-plane image for determining the 2DMA gave rise to error by 3.5 mm. The remaining difference was accounted for by the random measurement error.     

Modeling the stability of the human glenohumeral joint during external rotation

An analytical model of the human glenohumeral joint was developed to predict glenohumeral kinematics and investigate how the glenohumeral capsule and articular contact between the humeral head and the glenoid stabilize the joint. This was performed during a simulation of an apprehension clinical exam or the cocked phase of throwing, when the humerus is susceptible to anterior instability or dislocation. Contact between the joint surfaces was modeled using a deformable articular contact method and the capsule was modeled as five elements with the ability to wrap around the surface of the humeral head. Experimental measurements (Novotny et al., Journal of Shoulder and Elbow surgery, 1998, 7, 629-639) provided geometric data from four in vitro specimens and kinematic results to validate model predictions. Material properties were taken from the literature. An equilibrium approach was used with the forces and moments produced by the ligaments and surface contact balanced against those applied externally to the humerus during external rotation of the abducted and extended humerus. The six equilibrium equations were solved for the position and orientation of the humerus. The center of the humeral head translated posteriorly and superiorly with external rotation. Model predictions for translational and rotational ranges of motion were not significantly different from experimental findings; however, at individual moment increments, the model underestimated the external rotation and overestimated the superior-inferior position of the humerus relative to the glenoid. The anterior band of the inferior glenohumeral ligament increased in tension with external rotation, while the axillary pouch and posterior band decreased in tension. Contact area, stress and force increased with external rotation and the contact area moved posteriorly and inferiorly around the rim of the glenoid. The model results provide information on how the relationship between the ligament element tensions and contact forces may act to avoid glenohumeral instability.