Musculoskeletal model of the upper limb based on the visible human male dataset

A mathematical model of the human upper limb was developed based on high-resolution medical images of the muscles and bones obtained from the Visible Human Male (VHM) 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].     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

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.     

A method for estimating torque-vector directions of shoulder muscles using surface EMGs

In this study, a new method is proposed to estimate the torque-vector directions of each shoulder muscle. The method is based on a multiple regression model that reconstructs shoulder torque, which is calculated from the hand force and posture, from the surface EMG of many muscles recorded simultaneously. The torque-vector directions of eleven shoulder muscles of four subjects were obtained at up to 30 different arm postures with this method. The mean confidence interval (p < 0.05) of the estimated torque-vector direction of each subject was 7.7-10.6 degrees. The correlation coefficient between the measured shoulder torque and reconstructed shoulder torque was between 0.76-0.84. The results for majority of the muscles were in accordance with previous studies, and reasonable from the viewpoint of anatomy. The torque-vector directions of a muscle, which are estimated with this method, have more of a functional meaning than a pure anatomical or mechanical one. These indicate the direction of the shoulder torque accompanying the muscle activation for a normal shoulder action that involves the cooperative contraction of many muscles.     

A method for estimating torque-vector directions of shoulder muscles using surface EMGs

In this study, a new method is proposed to estimate the torque-vector directions of each shoulder muscle. The method is based on a multiple regression model that reconstructs shoulder torque, which is calculated from the hand force and posture, from the surface EMG of many muscles recorded simultaneously. The torque-vector directions of eleven shoulder muscles of four subjects were obtained at up to 30 different arm postures with this method. The mean confidence interval ( p< 0.05) of the estimated torque-vector direction of each subject was 7.7-10.6 degrees. The correlation coefficient between the measured shoulder torque and reconstructed shoulder torque was between 0.76-0.84. The results for majority of the muscles were in accordance with previous studies, and reasonable from the viewpoint of anatomy. The torque-vector directions of a muscle, which are estimated with this method, have more of a functional meaning than a pure anatomical or mechanical one. These indicate the direction of the shoulder torque accompanying the muscle activation for a normal shoulder action that involves the cooperative contraction of many muscles.     

A musculoskeletal model of the hand and wrist: model definition and evaluation

Abstract

To improve our understanding on the neuromechanics of finger movements, a comprehensive musculoskeletal model is needed. The aim of this study was to build a musculoskeletal model of the hand and wrist, based on one consistent data set of the relevant anatomical parameters. We built and tested a model including the hand and wrist segments, as well as the muscles of the forearm and hand in OpenSim. In total, the model comprises 19 segments (with the carpal bones modeled as one segment) with 23 degrees of freedom and 43 muscles. All required anatomical input data, including bone masses and inertias, joint axis positions and orientations as well as muscle morphological parameters (i.e. PCSA, mass, optimal fiber length and tendon length) were obtained from one cadaver of which the data set was recently published. Model validity was investigated by first comparing computed muscle moment arms at the index finger metacarpophalangeal (MCP) joint and wrist joint to published reference values. Secondly, the muscle forces during pinching were computed using static optimization and compared to previously measured intraoperative reference values. Computed and measured moment arms of muscles at both index MCP and wrist showed high correlation coefficients (r?=?0.88 averaged across all muscles) and modest root mean square deviation (RMSD?=?23% averaged across all muscles). Computed extrinsic flexor forces of the index finger during index pinch task were within one standard deviation of previously measured in-vivo tendon forces. These results provide an indication of model validity for use in estimating muscle forces during static tasks.     

Upper limb H reflexes and somatosensory evoked potentials modulated by movement.

In the human lower limb, the magnitudes of both Hoffmann (H) reflexes and primary somatosensory evoked potentials (SEPs) from scalp electrodes, are reduced by active and/or passive movement. We surmised that similar effects occur for the upper limb and specifically hypothesised that amplitudes of median nerve induced flexor carpii radialis H reflexes and cortical SEPs are reduced with passive movement about the wrist or elbow. The results showed (P<0. 05) that either movement significantly attenuated mean magnitudes of SEPs elicited from stimulation at elbow or wrist and that reflex magnitudes attenuated with wrist movement. Thus, the upper limb shows similar movement-induced modulation to the lower limb. These attenuations of fast conducting sensory paths consequent to movement per se, may be a basic level of motor control, initiated from muscle mechanoreceptor discharge. Upon this basic level, more complex modulations then may be laid as appropriate for the particular characteristics of active motor tasks.     

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.     

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.