Musculoskeletal Model of the Upper Limb Based on the Visible Human Male Dataset

A mathematical model of Ihe human upper limb was developed based on high-resolution medical images of the muscles and bones obtained from the Visible Human Male ( HM) project. Three-dimensional surfaces of the muscles and bones were reconstructed from Computed Tomography (CT) images and Color Cryosection images obtained from the VHM cadaver. Thirteen degrees of freedom were used to describe the orientations of seven bones in the model: clavicle, scapula, humerus, radius, ulna, carpal bones, and hand. All of the major articulations from the shoulder girdle down to the wrist were included in the model. The model was actuated by 42 muscle bundles, which represented the actions of 26 muscle groups in the upper limb. The paths of the muscles were modeled using a new approach called the Obstacle-set Method (33) The calculated paths of the muscles were verified by comparing the muscle moment arms computed in the model with the results of anatomical studies reported in the literature, In-vivo measurements of maximum isometric muscle torques developed at the shoulder, elbow, and wrist were also used to estimate the architectural properties of each musculotendon actuator in the model. The entire musculoskeletal model can be reconstructed using the data given in this paper, along with information presented in a companion paper which defines the kinematic structure of the model (26)     

Effects of postural muscle fatigue on the relation between segmental posture and movement.

The purpose of this study was to examine whether fatigue of postural muscles might influence the coordination between segmental posture and movement. Seven healthy adults performed series of fifteen fast wrist flexions and extensions while being instructed to keep a dominant upper limb posture as constant as possible. These series of voluntary movements were performed before and after a fatiguing submaximal isometric elbow flexion, and also with or without the help of an elbow support. Surface EMG from muscles Delto?deus anterior, Biceps brachii, Triceps brachii, Flexor carpi ulnaris, Extensor carpi radialis were recorded simultaneously with wrist, elbow and shoulder accelerations and wrist and elbow displacements. Fatigue was evidenced by a shift of the elbow and shoulder muscles EMG spectra towards low frequencies. Kinematics of wrist movements and corresponding activations of wrist prime-movers, as well as the background of postural muscle activation before wrist movement were not modified. There were only slight changes in timing of postural muscle activations. These data indicate that postural fatigue induced by a low-level isometric contraction has no effect on voluntary movement and requires no dramatic adaptation in postural control.     

Upper limb muscle volumes in adult subjects

Muscle force-generating properties are often derived from cadaveric studies of muscle architecture. While the relative sizes of muscles at a single upper limb joint have been established in cadaveric specimens, the relative sizes of muscles across upper limb joints in living subjects remain unclear. We used magnetic resonance imaging to measure the volumes of the 32 upper limb muscles crossing the glenohumeral joint, elbow, forearm, and wrist in 10 young, healthy subjects, ranging from a 20th percentile female to a 97th percentile male, based on height. We measured the volume and volume fraction of these muscles. Muscles crossing the shoulder, elbow, and wrist comprised 52.5, 31.4, and 16.0% of the total muscle volume, respectively. The deltoid had the largest volume fraction (15.2%+/-1%) and the extensor indicis propius had the smallest (0.2%+/-0.05%). We determined that the distribution of muscle volume in the upper limb is highly conserved across these subjects with a three-fold variation in total muscle volumes (1427-4426cm(3)). When we predicted the volume of an individual muscle from the mean volume fraction, on average 85% of the variation among subjects was accounted for (average p=0.0008). This study provides normative data that forms the basis for investigating muscle volumes in other populations, and for scaling computer models to more accurately represent the muscle volume of a specific individual.     

Moment-generating capacity of upper limb muscles in healthy adults

Muscle strength and volume vary greatly among individuals. Maximum isometric joint moment, a standard measurement of strength, has typically been assessed in young, healthy subjects, whereas muscle volumes have generally been measured in cadavers. This has made it difficult to characterize the relationship between isometric strength and muscle size in humans. We measured maximum isometric moments about the shoulder, elbow, and wrist in 10 young, healthy subjects, ranging in size from a 20th percentile female to a 97th percentile male. The volumes of 32 upper limb muscles were determined from magnetic resonance images of these same subjects, and grouped according to their primary function. The maximum moments produced using the shoulder adductors (67.9+/-28.4 Nm) were largest, and were approximately 6.5(+/-1.2) times greater than those produced using the wrist extensors (10.2+/-4.6 Nm), which were smallest. While there were substantial differences in moment-generating capacity among these 10 subjects, moment significantly covaried with muscle volume of the appropriate functional group, explaining between 95% (p<0.0001; shoulder adductors) and 68% (p=0.004; wrist flexors) of the variation in the maximum isometric joint moments among subjects. While other factors, such as muscle moment arms or neural activation and coordination, can contribute to variation in strength among subjects, they either were relatively constant across these subjects compared to large differences in muscle volumes or they covaried with muscle volume. We conclude that differences in strength among healthy young adults are primarily a consequence of variation in muscle volume, as opposed to other factors.     

Simulation of extension,radial and ulnar deviation of the wrist with a rigid body spring model

A novel computational model of the wrist that predicts carpal bone motion was developed in order to investigate the complex kinematics of the human wrist.This rigid body spring model (RBSM) of the wrist was built using surface models of the eight carpal bones, the bases of the five metacarpal bones, and the distal parts of the ulna and radius, all obtained from computed tomography (CT) scans of a cadaver upper limb. Elastic contact conditions between the rigid bodies modeled the influence of the cartilage layers, and ligamentous structures were constructed using nonlinear, tension-only spring elements. Motion of the wrist was simulated by applying forces to the tendons of the five main wrist muscles modeled.Three wrist motions were simulated: extension, ulnar deviation and radial deviation. The model was tested and tuned by comparing the simulated displacement and orientation of the carpal bones with previously obtained CT-scans of the same cadaver arm in deviated (45°ulnar and 15°radial), and extended (57°) wrist positions. Simulation results for the scaphoid, lunate, capitate, hamate and triquetrum are presented here and provide credible prediction of carpal bone movement. These are the first reported results of such a model. They indicate promise that this model will assist in future wrist kinematics investigations. However, further optimization and validation are required to define and guarantee the validity of results.     

Long-term repeatability of force, endurance time and muscle activity during isometric contractions.

We determined the repeatability and correlations between force, endurance and muscle activity during isometric contractions over three years. Twenty-six subjects, with and without complaints of the shoulder and neck, performed standardized maximal and submaximal shoulder-abduction contractions and wrist extension-contractions at yearly intervals from 1997 to 1999. Peak forces developed during maximal contraction and the endurance times of submaximal contractions during shoulder abduction and wrist extension were measured. Electromyography (EMG) of muscle activity was recorded bilaterally from the upper trapezius, middle deltoid, and forearm extensor muscles. Root mean square EMG amplitudes were calculated. We found statistically significant associations between peak forces developed during wrist extension and shoulder abduction, and between endurance times of submaximal wrist extension and shoulder abduction. No statistically significant changes in peak force and EMG(peak) were found over the measurement years. The responses were not statistically significantly influenced by gender, or neck and shoulder pain. However, we observed considerable intra-individual variation in the inter-year measurements particularly for the responses to submaximal contraction. Such large variations represent a challenge when attempting to use the responses to interpret the effects of therapies.     

Investigation of the Peculiarities of Two-joint Muscles Using a 3 DOF Model of the Human Upper Limb in the Sagittal Plane: an Optimization Approach

Abstract

The purpose of this paper is an investigation of the peculiarities of biarticular muscles by means of modelling and analytical solution of the indeterminate problem. The basic model includes 10 muscle elements performing flexio/extensio in the shoulder, elbow and wrist. Four of them are biarticular muscles. Two modifications of the model with only monoarticular muscles are developed. The indeterminate problem is solved analytically using the objective criterion σc_iF_i ² where F( is the module of the i-th muscle force and Cj is a weight coefficient. The predicted muscle forces, joint reactions and moments are compared in-between the basic model and its two modifications for different joint angles, external loading and weight coefficients. The main conclusions are: it is impossible to formulate strict advantages of the biarticular muscles under quasistatical conditions, their peculiarities depend on limb position, external loading and neural control; in general, monoarticular muscles are more powerful than biarticular ones; the biarticular muscles fine tune muscle coordination, their control is more precise and graceful; the presence of biarticular muscles leads to an increase of the joint reactions and moments, thus stabilizing the limb.     

Characterizing upper limb muscle volume and strength in older adults: a comparison with young adults

Aging is associated with the loss of muscle volume (MV) and force leading to difficulties with activities of daily living. However, the relationship between upper limb MV and joint strength has not been characterized for older adults. Quantifying this relationship may help our understanding of the functional declines of the upper limb that older adults experience. Our objective was to assess the relationship between upper limb MV and maximal isometric joint moment-generating capacity (IJM) in a single cohort of healthy older adults (age ≥ 65 years) for 6 major functional groups (32 muscles). MV was determined from MRI for 18 participants (75.1±4.3 years). IJM at the shoulder (abduction/adduction), elbow (flexion/extension), and wrist (flexion/extension) was measured. MV and IJM measurements were compared to previous reports for young adults (28.6±4.5 years). On average older adults had 16.5% less total upper limb MV compared to young adults. Additionally, older adult wrist extensors composed a significantly increased percentage of upper limb MV. Older adult IJM was reduced across all joints, with significant differences for shoulder abductors (p<0.0001), adductors (p=0.01), and wrist flexors (p<0.0001). Young adults were strongest at the shoulder, which was not the case for older adults. In older adults, 40.6% of the variation in IJM was accounted for by MV changes (p≤0.027), compared to 81.0% in young adults. We conclude that for older adults, MV and IJM are, on average, reduced but the significant linear relationship between MV and IJM is maintained. These results suggest that older adult MV and IJM cannot be simply scaled from young adults.     

3D finite element models of shoulder muscles for computing lines of actions and moment arms

Accurate representation of musculoskeletal geometry is needed to characterise the function of shoulder muscles. Previous models of shoulder muscles have represented muscle geometry as a collection of line segments, making it difficult to account for the large attachment areas, muscle–muscle interactions and complex muscle fibre trajectories typical of shoulder muscles. To better represent shoulder muscle geometry, we developed 3D finite element models of the deltoid and rotator cuff muscles and used the models to examine muscle function. Muscle fibre paths within the muscles were approximated, and moment arms were calculated for two motions: thoracohumeral abduction and internal/external rotation. We found that muscle fibre moment arms varied substantially across each muscle. For example, supraspinatus is considered a weak external rotator, but the 3D model of supraspinatus showed that the anterior fibres provide substantial internal rotation while the posterior fibres act as external rotators. Including the effects of large attachment regions and 3D mechanical interactions of muscle fibres constrains muscle motion, generates more realistic muscle paths and allows deeper analysis of shoulder muscle function.     

A musculoskeletal model of the upper extremity for use in the development of neuroprosthetic systems

Upper extremity neuroprostheses use functional electrical stimulation (FES) to restore arm motor function to individuals with cervical level spinal cord injury. For the design and testing of these systems, a biomechanical model of the shoulder and elbow has been developed, to be used as a substitute for the human arm. It can be used to design and evaluate specific implementations of FES systems, as well as FES controllers. The model can be customized to simulate a variety of pathological conditions. For example, by adjusting the maximum force the muscles can produce, the model can be used to simulate an individual with tetraplegia and to explore the effects of FES of different muscle sets. The model comprises six bones, five joints, nine degrees of freedom, and 29 shoulder and arm muscles. It was developed using commercial, graphics-based modeling and simulation packages that are easily accessible to other researchers and can be readily interfaced to other analysis packages. It can be used for both forward-dynamic (inputs: muscle activation and external load; outputs: motions) and inverse-dynamic (inputs: motions and external load; outputs: muscle activation) simulations. Our model was verified by comparing the model calculated muscle activations to electromyographic signals recorded from shoulder and arm muscles of five subjects. As an example of its application to neuroprosthesis design, the model was used to demonstrate the importance of rotator cuff muscle stimulation when aiming to restore humeral elevation. It is concluded that this model is a useful tool in the development and implementation of upper extremity neuroprosthetic systems.     

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.     

Long-latency reflexes of the human arm reflect an internal model of limb dynamics

A key feature of successful motor control is the ability to counter unexpected perturbations. This process is complicated in multijoint systems, like the human arm, by the fact that loads applied at one joint will create motion at other joints [1-3]. Here, we test whether our most rapid corrections, i.e., reflexes, address this complexity through an internal model of the limb's mechanical properties. By selectively applying torque perturbations to the subject's shoulder and/or elbow, we revealed a qualitative difference between the arm's short-latency/spinal reflexes and long-latency/cortical reflexes. Short-latency reflexes of shoulder muscles were linked exclusively to shoulder motion, whereas its long-latency reflexes were sensitive to both shoulder and elbow motion, i.e., matching the underlying shoulder torque. In fact, a long-latency reflex could be evoked without even stretching or lengthening the shoulder muscle but by displacing just the elbow joint. Further, the shoulder's long-latency reflexes were appropriately modified across the workspace to account for limb-geometry changes that affect the transformation between joint torque and joint motion. These results provide clear evidence that long-latency reflexes possess an internal model of limb dynamics, a degree of motor intelligence previously reserved for voluntary motor control [3-5]. The use of internal models for both voluntary and reflex control is consistent with substantial overlap in their neural substrates and current notions of intelligent feedback control [6-8].     

Moving muscle points provide accurate curved muscle paths in a model of the cervical spine

Muscle paths in musculoskeletal models have been modeled using several different methods; however, deformation of soft tissue with changes in posture is rarely accounted for, and often only the neutral posture is used to define a muscle path. The objective of this study was to model curved muscle paths in the cervical spine that take into consideration soft tissue deformation with changes in neck posture. Two subject-specific models were created from magnetic resonance images (MRI) in 5 different sagittal plane neck postures. Curved paths of flexor and extensor muscles were modeled using piecewise linear lines-of-action in two ways; (1) using fixed via points determined from muscle paths in the neutral posture and (2) using moving muscle points that moved relative to the bones determined from muscle paths in all 5 postures. Accuracy of each curved modeled muscle path was evaluated by an error metric, the distance from the anatomic (centroid) muscle path determined from the MRI. Error metric was compared among three modeled muscle path types (straight, fixed via and moving muscle point) using a repeated measures one-way ANOVA (α=0.05). Moving muscle point paths had 21% lower error metric than fixed via point paths over all 15 pairs of neck muscles examined over 5 postures (3.86 mm vs. 4.88 mm). This study highlights the importance of defining muscle paths in multiple postures in order to properly define the changing curvature of a muscle path due to soft tissue deformation with posture.     

Application of spherical and cylindrical wrapping algorithms in a musculoskeletal model of the upper limb

In the modelling of the upper limb, many muscles cannot be represented as a straight line from origin to insertion due to the complex morphology causing them to wrap around passive structures. The majority of bony contours that form these obstructions can be described adequately as simple geometric shapes such as spheres and cylinders.A novel technique for the parameterisation of muscle paths as they wrap around such shapes has been developed for use in an upper limb model. The new method involves the definition of moving co-ordinate systems in which the path of a wrapped muscle does not move, allowing simplified specification. In addition, an analytical calculation of the wrapping path around a cylinder is presented over previous approximate methods.Muscle moment arms were pre-calculated from vector considerations and within SIMM by tendon excursion. Close agreement between the two suggests that the proposed implementations accurately follow the theoretical relationship and can be used with confidence in musculoskeletal models.     

Muscle architecture of the upper limb in the orangutan

We dissected the left upper limb of a female orangutan and systematically recorded muscle mass, fascicle length, and physiological cross-sectional area (PCSA), in order to quantitatively clarify the unique muscle architecture of the upper limb of the orangutan. Comparisons of the musculature of the dissected orangutan with corresponding published chimpanzee data demonstrated that in the orangutan, the elbow flexors, notably M. brachioradialis, tend to exhibit greater PCSAs. Moreover, the digital II-V flexors in the forearm, such as M. flexor digitorum superficialis and M. flexor digitorum profundus, tend to have smaller PCSA as a result of their relatively longer fascicles. Thus, in the orangutan, the elbow flexors demonstrate a higher potential for force production, whereas the forearm muscles allow a greater range of wrist joint mobility. The differences in the force-generating capacity in the upper limb muscles of the two species might reflect functional specialization of muscle architecture in the upper limb of the orangutan for living in arboreal environments.     

Shoulder muscle function depends on elbow joint position: an illustration of dynamic coupling in the upper limb

Shoulder muscle function has been documented based on muscle moment arms, lines of action and muscle contributions to contact force at the glenohumeral joint. At present, however, the contributions of individual muscles to shoulder joint motion have not been investigated, and the effects of shoulder and elbow joint position on shoulder muscle function are not well understood. The aims of this study were to compute the contributions of individual muscles to motion of the glenohumeral joint during abduction, and to examine the effect of elbow flexion on shoulder muscle function. A three-dimensional musculoskeletal model of the upper limb was used to determine the contributions of 18 major muscles and muscle sub-regions of the shoulder to glenohumeral joint motion during abduction. Muscle function was found to depend strongly on both shoulder and elbow joint positions. When the elbow was extended, the middle and anterior deltoid and supraspinatus were the greatest contributors to angular acceleration of the shoulder in abduction. In contrast, when the elbow was flexed at 90°, the anterior deltoid and subscapularis were the greatest contributors to joint angular acceleration in abduction. This dependence of shoulder muscle function on elbow joint position is explained by the existence of dynamic coupling in multi-joint musculoskeletal systems. The extent to which dynamic coupling affects shoulder muscle function, and therefore movement control, is determined by the structure of the inverse mass matrix, which depends on the configuration of the joints. The data provided may assist in the diagnosis of abnormal shoulder function, for example, due to muscle paralysis or in the case of full-thickness rotator cuff tears.     

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.