Contributions of the individual muscles of the shoulder to glenohumeral joint stability during abduction

The aim of this study was to determine the relative contributions of the deltoid and rotator cuff muscles to glenohumeral joint stability during arm abduction. A three-dimensional model of the upper limb was used to calculate the muscle and joint-contact forces at the shoulder for abduction in the scapular plane. The joints of the shoulder girdle-sternoclavicular joint, acromioclavicular joint, and glenohumeral joint-were each represented as an ideal three degree-of-freedom ball-and-socket joint. The articulation between the scapula and thorax was modeled using two kinematic constraints. Eighteen muscle bundles were used to represent the lines of action of 11 muscle groups spanning the glenohumeral joint. The three-dimensional positions of the clavicle, scapula, and humerus during abduction were measured using intracortical bone pins implanted into one subject. The measured bone positions were inputted into the model, and an optimization problem was solved to calculate the forces developed by the shoulder muscles for abduction in the scapular plane. The model calculations showed that the rotator cuff muscles (specifically, supraspinatus, subscapularis, and infraspinatus) by virtue of their lines of action are perfectly positioned to apply compressive load across the glenohumeral joint, and that these muscles contribute most significantly to shoulder joint stability during abduction. The middle deltoid provides most of the compressive force acting between the humeral head and the glenoid, but this muscle also creates most of the shear, and so its contribution to joint stability is less than that of any of the rotator cuff muscles.     

Glenohumeral stability in simulated rotator cuff tears

Rotator cuff tears disrupt the force balance in the shoulder and the glenohumeral joint in particular, resulting in compromised arm elevation torques. The trade-off between glenohumeral torque and glenohumeral stability is not yet understood. We hypothesize that compensation of lost abduction torque will lead to a superior redirection of the reaction force vector onto the glenoid surface, which will require additional muscle forces to maintain glenohumeral stability. Muscle forces in a single arm position for five combinations of simulated cuff tears were estimated by inverse dynamic simulation (Delft Shoulder and Elbow Model) and compared with muscle forces in the non-injured condition. Each cuff tear condition was simulated both without and with an active modeling constraint for glenohumeral stability, which was defined as the condition in which the glenohumeral reaction force intersects the glenoid surface. For the simulated position an isolated tear of the supraspinatus only increased the effort of the other muscles with 8%, and did not introduce instability. For massive cuff tears beyond the supraspinatus, instability became a prominent factor: the deltoids were not able to fully compensate lost net abduction torque without introducing destabilizing forces; unfavorable abductor muscles (i.e. in the simulated position the subscapularis and the biceps longum) remain to compensate the necessary abduction torque; the teres minor appeared to be of vital importance to maintain glenohumeral stability. Adverse adductor muscle co-contraction is essential to preserve glenohumeral stability.     

Comparison of an EMG-based and a stress-based method to predict shoulder muscle forces

The estimation of muscle forces in musculoskeletal shoulder models is still controversial. Two different methods are widely used to solve the indeterminacy of the system: electromyography (EMG)-based methods and stress-based methods. The goal of this work was to evaluate the influence of these two methods on the prediction of muscle forces, glenohumeral load and joint stability after total shoulder arthroplasty. An EMG-based and a stress-based method were implemented into the same musculoskeletal shoulder model. The model replicated the glenohumeral joint after total shoulder arthroplasty. It contained the scapula, the humerus, the joint prosthesis, the rotator cuff muscles supraspinatus, subscapularis and infraspinatus and the middle, anterior and posterior deltoid muscles. A movement of abduction was simulated in the plane of the scapula. The EMG-based method replicated muscular activity of experimentally measured EMG. The stress-based method minimised a cost function based on muscle stresses. We compared muscle forces, joint reaction force, articular contact pressure and translation of the humeral head. The stress-based method predicted a lower force of the rotator cuff muscles. This was partly counter-balanced by a higher force of the middle part of the deltoid muscle. As a consequence, the stress-based method predicted a lower joint load (16% reduced) and a higher superior–inferior translation of the humeral head (increased by 1.2 mm). The EMG-based method has the advantage of replicating the observed cocontraction of stabilising muscles of the rotator cuff. This method is, however, limited to available EMG measurements. The stress-based method has thus an advantage of flexibility, but may overestimate glenohumeral subluxation.     

Modulation of shoulder muscle and joint function using a powered upper-limb exoskeleton

Robotic-assistive exoskeletons can enable frequent repetitive movements without the presence of a full-time therapist; however, human-machine interaction and the capacity of powered exoskeletons to attenuate shoulder muscle and joint loading is poorly understood. This study aimed to quantify shoulder muscle and joint force during assisted activities of daily living using a powered robotic upper limb exoskeleton (ArmeoPower, Hocoma). Six healthy male subjects performed abduction, flexion, horizontal flexion, reaching and nose touching activities. These tasks were repeated under two conditions: (i) the exoskeleton compensating only for its own weight, and (ii) the exoskeleton providing full upper limb gravity compensation (i.e., weightlessness). Muscle EMG, joint kinematics and joint torques were simultaneously recorded, and shoulder muscle and joint forces calculated using personalized musculoskeletal models of each subject’s upper limb. The exoskeleton reduced peak joint torques, muscle forces and joint loading by up to 74.8% (0.113 Nm/kg), 88.8% (5.8%BW) and 68.4% (75.6%BW), respectively, with the degree of load attenuation strongly task dependent. The peak compressive, anterior and superior glenohumeral joint force during assisted nose touching was 36.4% (24.6%BW), 72.4% (13.1%BW) and 85.0% (17.2%BW) lower than that during unassisted nose touching, respectively. The present study showed that upper limb weight compensation using an assistive exoskeleton may increase glenohumeral joint stability, since deltoid muscle force, which is the primary contributor to superior glenohumeral joint shear, is attenuated; however, prominent exoskeleton interaction moments are required to position and control the upper limb in space, even under full gravity compensation conditions. The modeling framework and results may be useful in planning targeted upper limb robotic rehabilitation tasks.     

Recruitment,force and fatigue characteristics of quadriceps muscles of paraplegics isometrically activated by surface functional electrical stimulation

This study deals with the recruitment characteristics of unfatigued electrically stimulated quadriceps muscles of paraplegic subjects and with the time-dependent force output of these muscles under sustained stimulation conditions. Both these aspects of the performance of paralysed stimulated muscles were studied under isometric conditions and at different muscle lengths. The forces in the knee joint resulting from stimulation of the quadriceps were also calculated. Recruitment force curves due to a ramp-like stimulation function indicated a strong dependence on muscle length and demonstrated a sigmoid-shaped curve with three distinct regions: negligible force up to threshold stimulation intensity; rapid force increase; and levelling-off of the curve after which the force remains constant even though intensity is further increased. When normalized to the maximal force, recruitment was found to be independent of muscle length, generating a typical recruitment curve for every patient, under isometric stimulation. The peak forces were obtained at the same flexion angles previously published for normal subjects, but with much lower values. Muscle fatigue in tetanic isometric conditions, defined as the decrease in force due to sustained stimulation with fixed parameters, was found to be length dependent and to have a double exponential decay. The first is the acute force loss and is the more significant for functional purposes; the second is the more moderate and asymptotic region, in which partial force receovery in the form of bursts is observed.     

The relation between increased deltoid activation and adductor muscle activation due to glenohumeral cuff tears

In patients with rotator cuff tears lost elevation moments are compensated for by increased deltoid activation. Concomitant proximal directed destabilizing forces at the glenohumeral joint are suggested to be compensated for by ‘out-of-phase’ adductor activation, preserving glenohumeral stability. Aim of this study was to demonstrate causality between moment compensating deltoid activation and stability compensating ‘out-of-phase’ adductor muscle activation.A differential arm loading with the same magnitude of forces applied at small and large moment arms relative to the glenohumeral joint was employed to excite deltoid activation, without externally affecting the force balance. Musculoskeletal modeling was applied to analyze the protocol in terms of muscle forces and glenohumeral (in)stability. The protocol was applied experimentally using electromyography (EMG) to assess muscle activation of healthy controls and cuff tear patients.Both modeling and experiments demonstrated increased deltoid activation with increased moment loading, which was higher in patients compared to controls. Model simulation of cuff tears demonstrated glenohumeral instability and related ‘out-of-phase’ adductor muscle activation which was also found experimentally in patients when compared to controls.Through differential moment loading, the assumed causal relation between increased deltoid activation and compensatory adductor muscle activation in cuff tear patients could be demonstrated. ‘Out-of-phase’ adductor activation in patients was attributed to glenohumeral instability. The moment loading protocol discerned patients with cuff tears from controls based on muscle activation.     

Relation between the ankle joint angle and the maximum isometric force of the toe flexor muscles

The purpose of this study was to investigate the relationships between the ankle joint angle and maximum isometric force of the toe flexor muscles. Toe flexor strength and electromyography activity of the foot muscles were measured in 12 healthy men at 6 different ankle joint angles with the knee joint at 90 deg in the sitting position. To measure the maximum isometric force of the toe flexor muscles, subjects exerted maximum force on a toe grip dynamometer while the activity levels of the intrinsic and extrinsic plantar muscles were measured. The relation between ankle joint angle and maximum isometric force of the toe flexor muscles was determined, and the isometric force exhibited a peak when the ankle joint was at 70–90 deg on average. From this optimal neutral position, the isometric force gradually decreased and reached its nadir in the plantar flexion position (i.e., 120 deg). The EMG activity of the abductor hallucis (intrinsic plantar muscle) and peroneus longus (extrinsic plantar muscle) did not differ at any ankle joint angles. The results of this study suggest that the force generation of toe flexor muscles is regulated at the ankle joint and that changes in the length-tension relations of the extrinsic plantar muscle could be a reason for the force-generating capacity at the metatarsophalangeal joint when the ankle joint angle is changed.     

Towards a model for force predictions in the human shoulder.

In this paper the concept of a three-dimensional biomechanical model of the human shoulder is introduced. This model is used to analyze static load sharing between the muscles, the bones and the ligaments. The model consists of all shoulder structures, which means that different positions and different load situations may be analyzed using the same model. Solutions can be found for the complete range of shoulder motion. However, this article focuses only on elevation in the scapular plane and on forces in structures attached to the humerus. The intention is to expand the model in future studies to also involve the forces acting on the other shoulder bones: the scapula and the clavicle. The musculoskeletal forces in the shoulder complex are predicted utilizing the optimization technique with the sum of squared muscle stresses as an objective function. Numerical results predict that among the muscles crossing the glenohumeral joint parts of the deltoideus, the infraspinatus, the supraspinatus, the subscapularis, the pectoralis major, the coracobrachialis and the biceps are the muscles most activated during this sort of abduction. Muscle-force levels reached values of 150 N when the hand load was 1 kg. The results from the model seem to be qualitatively accurate, but it is concluded that in the future development of the model the direction of the contact force in the glenohumeral joint must be constrained.     

Static optimization underestimates antagonist muscle activity at the glenohumeral joint: A musculoskeletal modeling study

Static optimization is commonly employed in musculoskeletal modeling to estimate muscle and joint loading; however, the ability of this approach to predict antagonist muscle activity at the shoulder is poorly understood. Antagonist muscles, which contribute negatively to a net joint moment, are known to be important for maintaining glenohumeral joint stability. This study aimed to compare muscle and joint force predictions from a subject-specific neuromusculoskeletal model of the shoulder driven entirely by measured muscle electromyography (EMG) data with those from a musculoskeletal model employing static optimization. Four healthy adults performed six sub-maximal upper-limb contractions including shoulder abduction, adduction, flexion, extension, internal rotation and external rotation. EMG data were simultaneously measured from 16 shoulder muscles using surface and intramuscular electrodes, and joint motion evaluated using video motion analysis. Muscle and joint forces were calculated using both a calibrated EMG-driven neuromusculoskeletal modeling framework, and musculoskeletal model simulations that employed static optimization. The EMG-driven model predicted antagonistic muscle function for pectoralis major, latissimus dorsi and teres major during abduction and flexion; supraspinatus during adduction; middle deltoid during extension; and subscapularis, pectoralis major and latissimus dorsi during external rotation. In contrast, static optimization neural solutions showed little or no recruitment of these muscles, and preferentially activated agonistic prime movers with large moment arms. As a consequence, glenohumeral joint force calculations varied substantially between models. The findings suggest that static optimization may under-estimate the activity of muscle antagonists, and therefore, their contribution to glenohumeral joint stability.     

Effect of abducting and adducting muscle activity on glenohumeral translation, scapular kinematics and subacromial space width in vivo

It is currently unknown in which ways activity of the ab- and adductor shoulder muscles affects shoulder biomechanics (scapular kinematics and glenohumeral translation), and whether these changes are relevant for alterations of the subacromial space width. The objective of this experimental in vivo study was thus to test the hypotheses that potential changes of the subacromial space width (during antagonistic muscle activity) are caused by alterations of scapular kinematics and/or glenohumeral translation. The shoulders of 12 healthy subjects were investigated with an open MRI-system at 30 degrees, 60 degrees, 90 degrees, 120 degrees and 150 degrees of arm elevation. A force of 15N was applied to the distal humerus, once causing isometric contraction of the abductors and once contraction of the adductors. The scapulo-humeral rhythm, scapular tilting and glenohumeral translation were calculated from the MR image data for both abducting and adducting muscle activity. Adducting muscle activity led to significant increase of the subacromial space width in all arm positions. The scapulo-humeral rhythm (2.2-2.5) and scapular tilting (2-4 degrees) remained relatively constant during elevation, no significant difference was found between abducting and adducting muscle activity. The position of the humerus relative to the glenoid was, however, significantly (p < 0.05) different (inferior and anterior) for adducting versus abducting muscle activity in midrange elevation (60-120 degrees). These data show that the subacromial space can be effectively widened by adducting muscle activity, by affecting the position of the humerus relative to the glenoid. This effect may be employed for conservative treatment of the impingement syndrome.     

Velocity-dependent cost function for the prediction of force sharing among synergistic muscles in a one degree of freedom model

Prediction of accurate and meaningful force sharing among synergistic muscles is a major problem in biomechanics research. Given a resultant joint moment, a unique set of muscle forces can be obtained from this mathematically redundant system using nonlinear optimization. The classical cost functions for optimization involve a normalization of the muscle forces to the absolute force capacity of the target muscles, usually by the cross-sectional area or the maximal isometric force. In a one degree of freedom model this leads to a functional relationship between moment arms and the predicted muscle forces, such that for constant moment arms, or constant ratios of moment arms, agonistic muscle forces increase or decrease in unison. Experimental studies have shown however that the relationship between muscle forces is highly task-dependent often causing forces to increase in one muscle while decreasing in a functional agonist, likely because of the contractile conditions and contractile properties of the involved muscles. We therefore, suggest a modified cost function that accounts for the instantaneous contraction velocity of the muscles and its effect on the instantaneous maximal force. With this novel objective function, a task-dependent prediction of muscle force distribution is obtained that allows, even in a one degree of freedom system, the prediction of force sharing loops, and simultaneously increasing and decreasing forces for agonist pairs of muscles.     

Force during stretches of rat skeletal muscles after hypertonia at short and long lengths

Following injection of tetanus toxin into rat gastrocnemius muscle to produce hypertonia, plantar flexor muscles were allowed to shorten (S, n=5) without restraint or held lengthened (L, n=3) by splinting. Saline injected rats served as control (n=5). One week after injection, peak forces during 3 stretches with passive muscles and acute isometric force deficits produced by 15 stretches of electrically stimulated muscles were examined under pentobarbital anesthesia. Isometric force and mass of plantar flexors were similar in S rats but 16% lower in L rats compared to control. Peak passive forces were highest in S rats but not different between L rats and control. At the end of the stretch protocol, isometric force deficits were 26% larger in S rats compared to L rats and 17% smaller in L rats compared to control. Acute isometric force deficits produced by stretches of active skeletal muscles were dependent on the muscle length maintained during hypertonia. Our animal model could be used to test rehabilitation interventions during hypertonia of skeletal muscles.     

Critical analysis of musculoskeletal modelling complexity in multibody biomechanical models of the upper limb

The inverse dynamics technique applied to musculoskeletal models, and supported by optimisation techniques, is used extensively to estimate muscle and joint reaction forces. However, the solutions of the redundant muscle force sharing problem are sensitive to the detail and modelling assumptions of the models used. This study presents four alternative biomechanical models of the upper limb with different levels of discretisation of muscles by bundles and muscle paths, and their consequences on the estimation of the muscle and joint reaction forces. The muscle force sharing problem is solved for the motions of abduction and anterior flexion, acquired using video imaging, through the minimisation of an objective function describing muscle metabolic energy consumption. While looking for the optimal solution, not only the equations of motion are satisfied but also the stability of the glenohumeral and scapulothoracic joints is preserved. The results show that a lower level of muscle discretisation provides worse estimations regarding the muscle forces. Moreover, the poor discretisation of muscles relevant to the joint in analysis limits the applicability of the biomechanical model. In this study, the biomechanical model of the upper limb describing the infraspinatus by a single bundle could not solve the complete motion of anterior flexion. Despite the small differences in the magnitude of the forces predicted by the biomechanical models with more complex muscular systems, in general, there are no significant variations in the muscular activity of equivalent muscles.     

Reduction of neuromuscular redundancy for postural force generation using an intrinsic stability criterion

Postural control requires the coordination of multiple muscles to achieve both endpoint force production and postural stability. Multiple muscle activation patterns can produce the required force for standing, but the mechanical stability associated with any given pattern may vary, and has implications for the degree of delayed neural feedback necessary for postural stability. We hypothesized that muscular redundancy is reduced when muscle activation patterns are chosen with respect to intrinsic musculoskeletal stability as well as endpoint force production. We used a three-dimensional musculoskeletal model of the cat hindlimb with 31 muscles to determine the possible contributions of intrinsic muscle properties to limb stability during isometric force generation. Using dynamic stability analysis we demonstrate that within the large set of activation patterns that satisfy the force requirement for posture, only a reduced subset produce a mechanically stable limb configuration. Greater stability in the frontal-plane suggests that neural control mechanisms are more highly active for sagittal-plane and for ankle joint control. Even when the limb was unstable, the time-constants of instability were sufficiently great to allow long-latency neural feedback mechanisms to intervene, which may be preferential for movements requiring maneuverability versus stability. Local joint stiffness of muscles was determined by the stabilizing or destabilizing effects of moment-arm versus joint angle relationships. By preferentially activating muscles with high local stiffness, muscle activation patterns with feedforward stabilizing properties could be selected. Such a strategy may increase intrinsic postural stability without co-contraction, and may be useful criteria in the force-sharing problem.     

Inter-joint coupling effects on muscle contributions to endpoint force and acceleration in a musculoskeletal model of the cat hindlimb

The biomechanical principles underlying the organization of muscle activation patterns during standing balance are poorly understood. The goal of this study was to understand the influence of biomechanical inter-joint coupling on endpoint forces and accelerations induced by the activation of individual muscles during postural tasks. We calculated induced endpoint forces and accelerations of 31 muscles in a 7 degree-of-freedom, three-dimensional model of the cat hindlimb. To test the effects of inter-joint coupling, we systematically immobilized the joints (excluded kinematic degrees of freedom) and evaluated how the endpoint force and acceleration directions changed for each muscle in 7 different conditions. We hypothesized that altered inter-joint coupling due to joint immobilization of remote joints would substantially change the induced directions of endpoint force and acceleration of individual muscles. Our results show that for most muscles crossing the knee or the hip, joint immobilization altered the endpoint force or acceleration direction by more than 90° in the dorsal and sagittal planes. Induced endpoint forces were typically consistent with behaviorally observed forces only when the ankle was immobilized. We then activated a proximal muscle simultaneous with an ankle torque of varying magnitude, which demonstrated that the resulting endpoint force or acceleration direction is modulated by the magnitude of the ankle torque. We argue that this simple manipulation can lend insight into the functional effects of co-activating muscles. We conclude that inter-joint coupling may be an essential biomechanical principle underlying the coordination of proximal and distal muscles to produce functional endpoint actions during motor tasks.     

Subject-specific thumb muscle activity during functional tasks of daily life

BackgroundThe trapeziometacarpal joint is subjected to high compressive forces during powerful pinch and grasp tasks due to muscle loading. In addition, muscle contraction is important for stability of the joint. The aim of the present study is to explore if different muscle activation patterns can be found between three functional tasks.MethodsIsometric forces and fine-wire electromyographic (fEMG) activity produced by three intrinsic and four extrinsic thumb muscles were measured in 10 healthy female volunteers. The participants performed isometric contractions in a lateral key pinch, a power grasp and a jar twist task. The tasks were executed with and without EMG recording to verify if electrode placement influenced force production.ResultsA subject-specific muscle recruitment was found which remained largely unchanged across tasks. Extrinsic thumb muscles were significantly more active than intrinsic muscles in all tasks. Insertion of the fEMG electrodes decreased force production significantly in all tasks.ConclusionThe thumb muscles display a high variability in muscle activity during functional tasks of daily life. The results of this study suggest that to produce a substantial amount of force, a well-integrated, but subject-specific, co-contraction between the intrinsic and extrinsic thumb muscles is necessary.     

A musculoskeletal shoulder simulation of moment arms and joint reaction forces after medialization of the supraspinatus footprint in rotator cuff repair

A non-anatomical reinsertion of the supraspinatus medially to the original footprint to avoid over-tensioning of the tendon in large and retracted tears is one surgical option in rotator cuff (RC) repair. The purpose of the study was to determine the biomechanical effects on the glenohumeral joint with regard to this surgical technique. A modified musculoskeletal computational shoulder model was used to evaluate the change in moment arms and muscle forces of the RC and the co-contracting muscles and the alteration of the joint reaction forces (compressive and shear forces) after reinsertion of the supraspinatus 5?mm, 10?mm, 15?mm and 20?mm medially to the original footprint. A medialization of the supraspinatus reduces its moment arm in glenohumeral abduction. In case of a medialization of the attachment of 15?mm and 20?mm, the supraspinatus restricts glenohumeral abduction at 54° and 68°. In glenohumeral forward flexion and in lower degrees of internal rotation the moment arm of the supraspinatus increases for a medialized tendon attachment and decreases in external rotation in relation to the anatomical condition. A medialization of the supraspinatus insertion point yields in an increase in muscle force for abduction, internal and external rotation. In the present model a medially non-anatomic reinsertion reduces significantly the compressive glenohumeral joint reaction and the glenohumeral stability. Moreover, the results show that a medialization of the supraspinatus leads to a reduction of the supraspinatus moment arm especially in abduction. This leads to an increase of a compensatory supraspinatus load for stabilization the humerus in space, which may potentially cause a postoperative overload of the tendon-bone-complex.     

Architectural properties of distal forelimb muscles in horses,Equus caballus

Articular injuries in athletic horses are associated with large forces from ground impact and from muscular contraction. To accurately and noninvasively predict muscle and joint contact forces, a detailed model of musculoskeletal geometry and muscle architecture is required. Moreover, muscle architectural data can increase our understanding of the relationship between muscle structure and function in the equine distal forelimb. Muscle architectural data were collected from seven limbs obtained from five thoroughbred and thoroughbred-cross horses. Muscle belly rest length, tendon rest length, muscle volume, muscle fiber length, and pennation angle were measured for nine distal forelimb muscles. Physiological cross-sectional area (PCSA) was determined from muscle volume and muscle fiber length. The superficial and deep digital flexor muscles displayed markedly different muscle volumes (227 and 656 cm3, respectively), but their PCSAs were very similar due to a significant difference in muscle fiber length (i.e., the superficial digital flexor muscle had very short fibers, while those of the deep digital flexor muscle were relatively long). The ulnaris lateralis and flexor carpi ulnaris muscles had short fibers (17.4 and 18.3 mm, respectively). These actuators were strong (peak isometric force, Fmax=5,814 and 4,017 N, respectively) and stiff (tendon rest length to muscle fiber length, LT:LMF=5.3 and 2.1, respectively), and are probably well adapted to stabilizing the carpus during the stance phase of gait. In contrast, the flexor carpi radialis muscle displayed long fibers (89.7 mm), low peak isometric force (Fmax=555 N), and high stiffness (LT:LMF=1.6). Due to its long fibers and low Fmax, flexor carpi radialis appears to be better adapted to flexion and extension of the limb during the swing phase of gait than to stabilization of the carpus during stance. Including muscle architectural parameters in a musculoskeletal model of the equine distal forelimb may lead to more realistic estimates not only of the magnitudes of muscle forces, but also of the distribution of forces among the muscles crossing any given joint.     

Determination of muscle loading at the hip joint for use in pre-clinical testing

The stability of joint endoprostheses depends on the loading conditions to which the implant-bone complex is exposed. Due to a lack of appropriate muscle force data, less complex loading conditions tend to be considered in vitro. The goal of this study was to develop a load profile that better simulates the in vivo loading conditions of a "typical" total hip replacement patient and considers the interdependence of muscle and joint forces. The development of the load profile was based on a computer model of the lower extremities that has been validated against in vivo data. This model was simplified by grouping functionally similar hip muscles. Muscle and joint contact forces were computed for an average data set of up to four patients throughout walking and stair climbing. The calculated hip contact forces were compared to the average of the in vivo measured forces. The final derived load profile included the forces of up to four muscles at the instances of maximum in vivo hip joint loading during both walking and stair climbing. The hip contact forces differed by less than 10% from the peak in vivo value for a "typical" patient. The derived load profile presented here is the first that is based on validated musculoskeletal analyses and seems achievable in an in vitro test set-up. It should therefore form the basis for further standardisation of pre-clinical testing by providing a more realistic approximation of physiological loading conditions.     

Individual muscle force parameters and fiber operating ranges for elbow flexion-extension and forearm pronation-supination

We have quantified individual muscle force and moment contributions to net joint moments and estimated the operating ranges of the individual muscle fibers over the full range of motion for elbow flexion/extension and forearm pronation/supination. A three dimensional computer graphics model was developed in order to estimate individual muscle contributions in each degree of freedom over the full range of motion generated by 17 muscles crossing the elbow and forearm. Optimal fiber length, tendon slack length, and muscle specific tension values were adjusted within the literature range from cadaver studies such that the net isometric joint moments of the model approximated experimental joint moments within one standard deviation. Analysis of the model revealed that the muscles operate on varying portions of the ascending limb, plateau region, and descending limb of the force-length curve. This model can be used to further understand isometric force and moment contributions of individual muscles to net joint moments of the arm and forearm and can serve as a comprehensive reference for the forces and moments generated by 17 major muscles crossing the elbow and wrist.