The sensitivity of endpoint forces produced by the extrinsic muscles of the thumb to posture

This study utilizes a biomechanical model of the thumb to estimate the force produced at the thumb-tip by each of the four extrinsic muscles. We used the principle of virtual work to relate joint torques produced by a given muscle force to the resulting endpoint force and compared the results to two separate cadaveric studies. When we calculated thumb-tip forces using the muscle forces and thumb postures described in the experimental studies, we observed large errors. When relatively small deviations from experimentally reported thumb joint angles were allowed, errors in force direction decreased substantially. For example, when thumb posture was constrained to fall within ±15° of reported joint angles, simulated force directions fell within experimental variability in the proximal–palmar plane for all four muscles. Increasing the solution space from ±1° to an unbounded space produced a sigmoidal decrease in error in force direction. Changes in thumb posture remained consistent with a lateral pinch posture, and were generally consistent with each muscle’s function. Altering thumb posture alters both the components of the Jacobian and muscle moment arms in a nonlinear fashion, yielding a nonlinear change in thumb-tip force relative to muscle force. These results explain experimental data that suggest endpoint force is a nonlinear function of muscle force for the thumb, support the continued use of methods that implement linear transformations between muscle force and thumb-tip force for a specific posture, and suggest the feasibility of accurate prediction of lateral pinch force in situations where joint angles can be measured accurately.     

Higher medially-directed joint reaction forces are a characteristic of dysplastic hips: A comparative study using subject-specific musculoskeletal models

Acetabular dysplasia is a known cause of hip osteoarthritis. In addition to abnormal anatomy, changes in kinematics, joint reaction forces (JRFs), and muscle forces could cause tissue damage to the cartilage and labrum, and may contribute to pain and fatigue. The objective of this study was to compare lower extremity joint angles, moments, hip JRFs and muscle forces during gait between patients with symptomatic acetabular dysplasia and healthy controls. Marker trajectories and ground reaction forces were measured in 10 dysplasia patients and 10 typically developing control subjects. A musculoskeletal model was scaled in OpenSim to each subject and subject-specific hip joint centers were determined using reconstructions from CT images. Joint kinematics and moments were calculated using inverse kinematics and inverse dynamics, respectively. Muscle forces and hip JRFs were estimated with static optimization. Inter-group differences were tested for statistical significance (p ≤ 0.05) and large effect sizes (d ≥ 0.8). Results demonstrated that dysplasia patients had higher medially directed JRFs. Joint angles and moments were mostly similar between the groups, but large inter-group effect sizes suggested some restriction in range of motion by patients at the hip and ankle. Higher medially-directed JRFs and inter-group differences in hip muscle forces likely stem from lateralization of the hip joint center in dysplastic patients. Joint force differences, combined with reductions in range of motion at the hip and ankle may also indicate compensatory strategies by patients with dysplasia to maintain joint stability.     

A model-based study of passive joint properties on muscle effort during static stance

This study examined the impact of lower extremity joint stiffnesses and simulated joint contractures on the muscle effort required to maintain static standing postures after a spinal cord injury (SCI). Static inverse computer simulations were performed with a three-dimensional 15 degree of freedom musculoskeletal model placed in 1600 different standing postures. The required lower extremity muscle forces were calculated through an optimization routine that minimized the sum of the muscle stresses squared, which was used as an index of the muscle effort required for each standing posture. Joint stiffnesses were increased and decreased by 100 percent of their nominal values, and contractures were simulated to determine their effects on the muscle effort for each posture. Nominal muscle and passive properties for an individual with a SCI determined the baseline muscle effort for comparisons. Stiffness changes for the ankle plantar flexion/dorsiflexion, hip flexion/extension, and hip abduction/adduction directions had the largest effect on reducing muscle effort by more than 5 percent, while changes in ankle inversion/eversion and knee flexion/extension had the least effect. For erect standing, muscle effort was reduced by more than 5 percent when stiffness was decreased at the ankle plantar flexion/dorsiflexion joint or hip flexion/extension joint. With simulated joint contractures, the postural workspace area decreased and muscle effort was not reduced by more than 5 percent for any posture. Using this knowledge, methods can be developed through the use of orthoses, physical therapy, surgery or other means to appropriately augment or diminish these passive moments during standing with a neuroprosthesis.     

Normal hip joint contact pressure distribution in single-leg standing--effect of gender and anatomic parameters.

A practical and easy-to-use analysis technique that can study the patient's hip joint contact force/pressure distribution would be useful to assess the effect of abnormal biomechanical conditions and anatomical deformities on joint contact stress for treatment planning purpose. This technique can also help to establish the normative database on hip joint contact pressure distribution in men and women in different age groups. Twelve anatomic parameters and seven biomechanical parameters of the hip joint in a normal population (41 females, 15 males) were calculated. The inter-parameter correlations were investigated. The pressure distribution in the hip joint was calculated using a three-dimensional discrete element analysis (DEA) technique. The 3D contact geometry of the hip joint was estimated from a 2D radiograph by assuming that the femoral head and the acetabular surface were spherical in shape. The head-trochanter ratio (HT), femoral head radius, pelvic height, the joint contact area, the normalized peak contact pressure, abductor force, and the joint contact force were significantly different between men and women. The normalized peak contact pressure was correlated both with acetabular coverage and head-trochanter ratio. Change of abductor force direction within normal variation did not affect the joint peak contact pressure. However, in simulated dysplastic conditions when the CE angle is small or negative, abductor muscle direction becomes very sensitive in joint contact pressure estimation. The models and the results presented can be used as the reference base in computer simulation for preoperative planning in pelvic or femoral osteotomy.     

Mechanical interaction of center of pressure and force direction in the upright human

Humans maintain upright bipedal posture by producing appropriate force against the environment through the interaction of neural controlled muscle force with the mechanics of the skeletal system. Characterizing these mechanics facilitates understanding of the neural control. We used a mechanical model of an upright human to analyze how the mechanical linkage aspects of the human body affect the force between the feet and the ground (F). Key parameters of F that directly regulate upright body posture are the direction of F (θ(F)) and its point of application (x(CP), anterior-posterior position of the center of pressure). Instantaneous analysis of the equations of motion demonstrated that θ(F) varied systematically with x(CP) such that the F vectors intersected at a point called the Posture-specific force Intersection point or PI (Π). The Π was located above the center of mass when the hip and knee joints were modeled as rigid and was located near the knee when the hip and knee torques were held constant. Limb posture and the knee torque affected the location of Π. This Π behavior quantifies the purely mechanical effect of anterior-posterior center of pressure shifts on the direction of F, which has consequences for the control of whole body posture.     

Crouched postures reduce the capacity of muscles to extend the hip and knee during the single-limb stance phase of gait

Many children with cerebral palsy walk in a crouch gait that progressively worsens over time, decreasing walking efficiency and leading to joint degeneration. This study examined the effect of crouched postures on the capacity of muscles to extend the hip and knee joints and the joint flexions induced by gravity during the single-limb stance phase of gait. We first characterized representative mild, moderate, and severe crouch gait kinematics based on a large group of subjects with cerebral palsy (N=316). We then used a three-dimensional model of the musculoskeletal system and its associated equations of motion to determine the effect of these crouched gait postures on (1) the capacity of individual muscles to extend the hip and knee joints, which we defined as the angular accelerations of the joints, towards extension, that resulted from applying a 1N muscle force to the model, and (2) the angular acceleration of the joints induced by gravity. Our analysis showed that the capacities of almost all the major hip and knee extensors were markedly reduced in a crouched gait posture, with the exception of the hamstrings muscle group, whose extension capacity was maintained in a crouched posture. Crouch gait also increased the flexion accelerations induced by gravity at the hip and knee throughout single-limb stance. These findings help explain the increased energy requirements and progressive nature of crouch gait in patients with cerebral palsy.     

Effect of hip joint angle at seat-off on hip joint contact force during sit-to-stand movement: a computer simulation study

Background

Sit-to-stand movements are a necessary part of daily life, and excessive mechanical stress on the articular cartilage has been reported to encourage the progression of osteoarthritis. Although a change in hip joint angle at seat-off may affect hip joint contact force during a sit-to-stand movement, the effect is unclear. This study aimed to examine the effect of the hip joint angle at seat-off on the hip joint contact force during a sit-to-stand movement by using a computer simulation.

Methods

A musculoskeletal model was created for the computer simulation, and eight muscles were attached to each lower limb. Various sit-to-stand movements were generated using parameters (e.g., seat height and time from seat-off to standing posture) reported by previous studies. The hip joint contact force for each sit-to-stand movement was calculated. Furthermore, the effect of the hip joint angle at seat-off on the hip joint contact force during the sit-to-stand movement was examined. In this study, as the changes to the musculoskeletal model parameters affect the hip joint contact force, a sensitivity analysis was conducted.

Results and conclusions

The hip joint contact force during the sit-to-stand movement increased approximately linearly as the hip flexion angle at the seat-off increased. Moreover, the normal sit-to-stand movement and the sit-to-stand movement yielding a minimum hip joint contact force were approximately equivalent. The effect of the changes to the musculoskeletal model parameters on the main findings of this study was minimal. Thus, the main findings are robust and may help prevent the progression of hip osteoarthritis by decreasing mechanical stress, which will be explored in future studies.

     

Sensitivity of temporomandibular joint force calculations to errors in muscle force measurements

A two-dimensional, five-muscle model was used to determine the degree of precision required for accurate calculation of temporomandibular joint force magnitude and direction. The sensitivity of the calculations to each variable were assessed by incrementing each variable through its presumed biological range and were expressed as rate of change in the joint force per unit change in each variable. Sensitivity of the calculations to variables depends upon both bite force direction and bite position. The bite force direction with maximum precision for joint force magnitude produced minimal precision for joint force direction. The accuracy needed for each muscle force varied greatly. The effect of error for each muscle parameter depended upon the magnitude, direction, and moment arm length of the muscle force relative to those of the resultant muscle force. If each of the five muscle forces was known to the nearest 1% of total muscle force magnitude, 1 degree of muscle force direction, and 1 mm of moment arm length, temporomandibular joint force magnitude could be calculated to the nearest 4 kg and joint force direction to the nearest 7 degrees. It is not known whether this precision for the muscle forces is possible.     

Anterior hip joint force increases with hip extension, decreased gluteal force, or decreased iliopsoas force

Abnormal or excessive force on the anterior hip joint may cause anterior hip pain, subtle hip instability and a tear of the acetabular labrum. We propose that both the pattern of muscle force and hip joint position can affect the magnitude of anterior joint force and thus possibly lead to excessive force and injury. The purpose of this study was to determine the effect of hip joint position and of weakness of the gluteal and iliopsoas muscles on anterior hip joint force. We used a musculoskeletal model to estimate hip joint forces during simulated prone hip extension and supine hip flexion under four different muscle force conditions and across a range of hip extension and flexion positions. Weakness of specified muscles was simulated by decreasing the modeled maximum force value for the gluteal muscles during hip extension and the iliopsoas muscle during hip flexion. We found that decreased force contribution from the gluteal muscles during hip extension and the iliopsoas muscle during hip flexion resulted in an increase in the anterior hip joint force. The anterior hip joint force was greater when the hip was in extension than when the hip was in flexion. Further studies are warranted to determine if increased utilization of the gluteal muscles during hip extension and of the iliopsoas muscle during hip flexion, and avoidance of hip extension beyond neutral would be beneficial for people with anterior hip pain, subtle hip instability, or an anterior acetabular labral tear.     

A model of semitendinosus muscle sarcomere length, knee and hip joint interaction in the frog hindlimb

The interaction between the semitendinosus muscle and both hip and knee joint angles was examined in the frog (Rana pipiens) hindlimb. Sarcomere length was measured by laser diffraction in passive muscle during hip and knee rotation. A model was then developed to predict semitendinosus sarcomere length as a function of both hip and knee flexion angle. Based on published frog muscle fiber length-tension [Gordon, A. M. et al., J. Physiol. 184, 170-192 (1966)] and force-velocity [Edman, K. A. P., J. Physiol. 291, 143-159 (1979)] properties, and published joint angles during hopping [Calow, L. J. and Alexander, R. McN., J. Zool. (Lond.) 171, 293-321 (1973)], muscle sarcomere length, force and hip and knee torque during a hop were predicted. The semitendinosus muscle generally operated on the descending limb of the length-tension curve at normal joint angle combinations. The model predicted that, during a single coordinated movement, a period of sarcomere shortening (concentric) was followed by a period of sarcomere lengthening (eccentric). Based on calculated torque profiles at the hip and knee joints, this study suggested that the semitendinosus muscle probably functions more as a hip extensor than a knee flexor. In addition, based on the nature of the shortening-lengthening cycle, the semitendinosus may act to mechanically link the force of knee extension to hip extension.     

Maximum voluntary joint torque as a function of joint angle and angular velocity: model development and application to the lower limb

Measurements of human strength can be important during analyses of physical activities. Such measurements have often taken the form of the maximum voluntary torque at a single joint angle and angular velocity. However, the available strength varies substantially with joint position and velocity. When examining dynamic activities, strength measurements should account for these variations. A model is presented of maximum voluntary joint torque as a function of joint angle and angular velocity. The model is based on well-known physiological relationships between muscle force and length and between muscle force and velocity and was tested by fitting it to maximum voluntary joint torque data from six different exertions in the lower limb. Isometric, concentric and eccentric maximum voluntary contractions were collected during hip extension, hip flexion, knee extension, knee flexion, ankle plantar flexion and dorsiflexion. Model parameters are reported for each of these exertion directions by gender and age group. This model provides an efficient method by which strength variations with joint angle and angular velocity may be incorporated into comparisons between joint torques calculated by inverse dynamics and the maximum available joint torques.     

A new parametric approach for modeling hip forces during gait

An analytical parametric model was developed to estimate the natural biological variations in muscle forces and their effect on the hip forces subject only to physiological constraints and not predefined optimization criterion. Force predictions are based on the joint kinematics and kinetics of each subject, a previously published muscle model, and physiological constraints on the muscle force distributions. The model was used to determine the hip contact forces throughout the stance phase of gait of a subject with a total hip replacement (THR). The parametrically modeled peak hip force without antagonistic muscle activity varied from 2.7 to 3.2 Body Weights (mean 2.9 Body Weights), which agreed well with published in vivo measurements from instrumented THRs in other subjects. For every 10% increase in antagonistic activity, the mean peak hip force increased by 0.2 Body Weights. The parametric model allows one to examine the effect of specific muscle weaknesses or increased antagonistic muscle activity on the hip forces. The model also provides a tool for studying the effect of gait adaptations on hip forces, as predictions are based on each individual's gait data. Differences in peak forces between subjects can then be evaluated relative to the uncertainty in not knowing the precise muscle force distributions.     

Contributions of individual muscles to hip joint contact force in normal walking

The human hip joint withstands high contact forces during daily activity and is therefore susceptible to injury and structural deterioration over time. Knowledge of muscle-force contributions to hip joint loading may assist in the development of strategies to prevent and manage conditions such as osteoarthritis, femoro-acetabular impingement and fracture. The main aim of this study was to determine the contributions of individual muscles to hip contact force in normal walking. Muscle contributions to hip contact force were calculated based on a previously published dynamic optimization solution for normal walking, which provided the time histories of joint motion, ground reaction forces, and muscle forces during the stance and swing phases of gait. The force developed by each muscle plus its contribution to the ground reaction force were used to determine the muscle’s contribution to hip contact force. Muscles were the major contributors to hip contact force, with gravitational and centrifugal forces combined contributing less than 5% of the total contact force. Four muscles that span the hip – gluteus medius, gluteus maximus, iliopsoas, and hamstrings – contributed most significantly to the three components of the hip contact force and hip contact impulse (integral of hip contact force over time). Three muscles that do not span the hip – vasti, soleus, and gastrocnemius – also contributed substantially to hip joint loading. These results provide additional insight into lower-limb muscle function during walking and may also be relevant to studies of cartilage degeneration and bone remodelling at the hip.     

Computer determination of contact stress distribution and size of weight bearing area in the human hip joint

The mathematical models and the corresponding computer program for determination of the hip joint contact force, the contact stress distribution, and the size of the weight bearing area from a standard anteroposterior radiograph are described. The described method can be applied in clinical practice to predict an optimal stress distribution after different operative interventions in the hip joint and to analyze the short and long term outcome of the treatment of various pathological conditions in the hip. A group of dysplastic hips and a group of normal hips were examined, with respect to the peak contact stress normalized by the body weight, and with respect to the functional angle of the weight bearing area. It is shown that both these parameters can be used in the assessment of hip dysplasia.     

Influence of altered gait patterns on the hip joint contact forces

Children who exhibit gait deviations often present a range of bone deformities, particularly at the proximal femur. Altered gait may affect bone growth and lead to deformities by exerting abnormal stresses on the developing bones. The objective of this study was to calculate variations in the hip joint contact forces with different gait patterns. Muscle and hip joint contact forces of four children with different walking characteristics were calculated using an inverse dynamic analysis and a static optimisation algorithm. Kinematic and kinetic analyses were based on a generic musculoskeletal model scaled down to accommodate the dimensions of each child. Results showed that for all the children with altered gaits both the orientation and magnitude of the hip joint contact force deviated from normal. The child with the most severe gait deviations had hip joint contact forces 30% greater than normal, most likely due to the increase in muscle forces required to sustain his crouched stance. Determining how altered gait affects joint loading may help in planning treatment strategies to preserve correct loading on the bone from a young age.     

Hip contact forces and gait patterns from routine activities.

In vivo loads acting at the hip joint have so far only been measured in few patients and without detailed documentation of gait data. Such information is required to test and improve wear, strength and fixation stability of hip implants. Measurements of hip contact forces with instrumented implants and synchronous analyses of gait patterns and ground reaction forces were performed in four patients during the most frequent activities of daily living. From the individual data sets an average was calculated. The paper focuses on the loading of the femoral implant component but complete data are additionally stored on an associated compact disc. It contains complete gait and hip contact force data as well as calculated muscle activities during walking and stair climbing and the frequencies of daily activities observed in hip patients. The mechanical loading and function of the hip joint and proximal femur is thereby completely documented. The average patient loaded his hip joint with 238% BW (percent of body weight) when walking at about 4 km/h and with slightly less when standing on one leg. This is below the levels previously reported for two other patients (Bergmann et al., Clinical Biomechanics 26 (1993) 969-990). When climbing upstairs the joint contact force is 251% BW which is less than 260% BW when going downstairs. Inwards torsion of the implant is probably critical for the stem fixation. On average it is 23% larger when going upstairs than during normal level walking. The inter- and intra-individual variations during stair climbing are large and the highest torque values are 83% larger than during normal walking. Because the hip joint loading during all other common activities of most hip patients are comparably small (except during stumbling), implants should mainly be tested with loading conditions that mimic walking and stair climbing.     

Computer Determination of Contact Stress Distribution and Size of Weight Bearing Area in the Human Hip Joint

The mathematical models and the corresponding computer program for determination of the hip joint contact force, the contact stress distribution, and the size of the weight bearing area from a standard anteroposterior radiograph are described. The described method can be applied in clinical practice to predict an optimal stress distribution after different operative interventions in the hip joint and to analyze the short and long term outcome of the treatment of various pathological conditions in the hip. A group of dysplastic hips and a group of normal hips were examined, with respect to the peak contact stress normalized by the body weight, and with respect to the functional angle of the weight bearing area. It is shown that both these parameters can be used in the assessment of hip dysplasia.     

Sensitivity of intervertebral joint forces to center of rotation location and trends along its migration path

Translational vertebral motion during functional tasks manifests itself in dynamic loci for center of rotation (COR). A shift of COR affects moment arms of muscles and ligaments; consequently, muscle and joint forces are altered. Based on posture- and level-specific trends of COR migration revealed by in vivo dynamic radiography during functional activities, it was postulated that the instantaneous COR location for a particular joint is optimized in order to minimize the joint reaction forces. A musculoskeletal multi-body model was employed to investigate the hypotheses that (1) a posterior COR in upright standing and (2) an anterior COR in forward flexed posture leads to optimized lumbar joint loads. Moreover, it was hypothesized that (3) lower lumbar levels benefit from a more superiorly located COR.The COR in the model was varied from its initial position in posterior-anterior and inferior-superior direction up to ±6 mm in steps of 2 mm. Movement from upright standing to 45° forward bending and backwards was simulated for all configurations. Joint reaction forces were computed at levels L2L3 to L5S1. Results clearly confirmed hypotheses (1) and (2) and provided evidence for the validity of hypothesis (3), hence offering a biomechanical rationale behind the migration paths of CORs observed during functional flexion/extension movement. Average sensitivity of joint force magnitudes to an anterior shift of COR was +6 N/mm in upright and −21 N/mm in 30° forward flexed posture, while sensitivity to a superior shift in upright standing was +7 N/mm and −8 N/mm in 30° flexion. The relation between COR loci and joint loading in upright and flexed postures could be mainly attributed to altered muscle moment arms and consequences on muscle exertion. These findings are considered relevant for the interpretation of COR migration data, the development of numerical models, and could have an implication on clinical diagnosis and treatment or the development of spinal implants.     

Hip biomechanics in orthopaedic clinical practice

Radiographic and clinical studies, coupled with biomechanical assessment of the hip, are important tools for predicting the development of osteoarthitis of the hip. In order to better understand the treatment of hip dysplasia, it is necessary to determine the contact stress in the hip joint. In this study, a three-dimensional mathematical model was used to determine hip joint contact stress. Because of the discrepancy in the results of analyses of different radiographic indicators of hip dysplasia, the calculation of hip joint contact stress is proposed for a more accurate assessment of the severity of hip dysplasia.