Computational method for muscle-path representation in musculoskeletal models

 This paper presents a new and efficient method to calculate the line-of-action of a muscle as it wraps over bones and other tissues on its way from origin to insertion. The muscle is assumed to be a one-dimensional, massless, taut string, and the surfaces of bones that the muscle may wrap around are approximated by cross-sectional boundaries obtained by slicing geometrical models of bones. Each cross-sectional boundary is approximated by a series of connected line segments. Thus, the muscle path to be calculated is piecewise linear with vertices being the contact points on the cross-sectional boundaries of the bones. Any level of geometric accuracy can be obtained by increasing the number of cross sections and the number of line segments in each cross section. The algorithm is computationally efficient even for large numbers of cross sections. Received: 9 July 2001 / Accepted in revised form: 11 March 2002     

The application of muscle wrapping to voxel-based finite element models of skeletal structures

Finite elements analysis (FEA) is now used routinely to interpret skeletal form in terms of function in both medical and biological applications. To produce accurate predictions from FEA models, it is essential that the loading due to muscle action is applied in a physiologically reasonable manner. However, it is common for muscle forces to be represented as simple force vectors applied at a few nodes on the model’s surface. It is certainly rare for any wrapping of the muscles to be considered, and yet wrapping not only alters the directions of muscle forces but also applies an additional compressive load from the muscle belly directly to the underlying bone surface. This paper presents a method of applying muscle wrapping to high-resolution voxel-based finite element (FE) models. Such voxel-based models have a number of advantages over standard (geometry-based) FE models, but the increased resolution with which the load can be distributed over a model’s surface is particularly advantageous, reflecting more closely how muscle fibre attachments are distributed. In this paper, the development, application and validation of a muscle wrapping method is illustrated using a simple cylinder. The algorithm: (1) calculates the shortest path over the surface of a bone given the points of origin and ultimate attachment of the muscle fibres; (2) fits a Non-Uniform Rational B-Spline (NURBS) curve from the shortest path and calculates its tangent, normal vectors and curvatures so that normal and tangential components of the muscle force can be calculated and applied along the fibre; and (3) automatically distributes the loads between adjacent fibres to cover the bone surface with a fully distributed muscle force, as is observed in vivo. Finally, we present a practical application of this approach to the wrapping of the temporalis muscle around the cranium of a macaque skull.     

The influence of limb alignment and transfemoral amputation technique on muscle capacity during gait

Many factors influence successful outcomes following transfemoral amputation. One factor is surgical technique. In this study, the influence of limb alignment and surgical technique on a muscle’s capacity to generate force was examined using musculoskeletal modeling. Non-amputee and transfemoral amputee models were analyzed while hip adduction, femur length, and reattached muscle wrap position, tension and stabilization technique were systematically varied. With muscle tension preserved, wrap position and femur length had little influence on muscle capacity. However, limb alignment, muscle tension and stabilization technique notably influenced muscle capacity. Overall, myodesis stabilization provided greater muscle balance and function than myoplasty stabilization.     

Muscle wrapping on arbitrary meshes with the heat method

Muscle paths play an important role in musculoskeletal simulations by determining a muscle’s length and how its force is distributed to joints. Most previous approaches estimate the way in which muscles ‘wrap’ around bones and other structures with smooth analytical wrapping surfaces. In this paper, we employ Newton’s method with discrete differential geometry to permit muscle wrapping over arbitrary polygonal mesh surfaces that represent underlying bones and structures. Precomputing distance fields allows us to speed up computations for the common situation where many paths cross the same wrapping surfaces. We found positive results for the accuracy, robustness, and efficiency of the method. However the method did not exhibit continuous changes in path length for dynamic simulations. Nonetheless this approach provides a valuable step toward fast muscle wrapping on arbitrary meshes.     

Optimal design parameters of the bicycle-rider system for maximal muscle power output

The purpose of this study was to find the optimal values of design parameters for a bicycle-rider system (crank length, pelvic inclination, seat height, and rate of crank rotation) which maximize the power output from muscles of the human lower limb during bicycling. The human lower limb was modelled as a planar system of five rigid bodies connected by four smooth pin joints and driven by seven functional muscle groups. The muscles were assumed to behave according to an adapted form of Hill's equation. The dependence of the average power on the design parameters was examined. The instantaneous power of each muscle group was studied and simultaneous activity of two seemingly antagonistic muscle groups was analyzed. Average peak power for one full pedal revolution was found to be around 1100 W. The upper body position corresponding to this peak power output was slightly reclined, and the pedalling rate was 155 rpm for a nominal crank length of 170 mm.     

Using a Spread Sheet Program to Model the Interaction of Low-Frequency Electric Fields with Inhomogeneous,Irregularly Shaped Objects

This paper describes the use of an ordinary business spread sheet program to calculate, using a two-stage finite difference method, the electric fields and current densities produced inside irregularly shaped models of the upper arm and forearm. The limb interiors are inhomogeneous, being represented as realistic cross-sections of bone and muscle.

A spread sheet forms a two-dimensional array. Each cell of the sheet can correspond to a physical element of space. Equations entered into the cells then represent the relationships among the potentials of the corresponding spatial elements.

The model is validated for a two-layer, lossy dielectric cylinder. Good agreement is obtained between the numerical and analytical solutions in this case except near the boundaries of the outer layer. The electric field within the limb depends on the shape and orientation of the limb relative to the applied field. The flow of current around the less conductive bones in the forearm can be observed.     

A Computational Approach for Simulating Muscle Morphologic Changes in Musculoskeletal Modeling

Digital data characterizing the geometry of anatomic structures (e.g., bones, muscles and tendons) are becoming readily available from Magnetic Resonance Imaging and Computerized Tomography technology. These data can be useful in forward simulations of limb dynamics to study the interaction between tissue morphology and limb dynamics, but only if computational tools are available to manipulate these data to simulate various structural changes. The objective of this project was to develop a computational approach to simulate physiological changes in muscle volume and limb cross-section. A previously reported method for calculating the area and area centroid location of complex shapes was combined with a newly derived algorithm that simulates muscle hypertrophy or atrophy. The new algorithm modifies the cross-sectional areas of specified muscles and the spatial orientation of other muscles, as appropriate, to simulate desired muscle volume changes. An approach of this type is needed to facilitate musculoskeletal modeling and computer simulations of movement designed to address questions related to the interactions between muscle morphology, limb inertial properties and limb dynamics.     

Mathematical model that predicts lower leg motion in response to electrical stimulation

Electrical stimulation of skeletal muscles of patients with upper motor neuron lesions can be used to restore functional movements such as standing or walking. Mathematical muscle models can assist in designing stimulation patterns that will enable patients to perform particular tasks more efficiently. In this study we extend our previous model to allow us to predict changes in knee joint angle in response to electrical stimulation of the human quadriceps femoris muscle. The model was tested both with and without inertial loads placed around the ankle joints of healthy subjects. Results showed that the model predicted the knee extensions with a RMS angle error that was generally 相似文献   

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.     

Motion Analytics of Trapezius Muscle Activity in an 18-Year-Old Female with Extended Upper Brachial Plexus Birth Palsy

Background  The trapezius muscle is often utilized as a muscle or nerve donor for repairing shoulder function in those with brachial plexus birth palsy (BPBP). To evaluate the native role of the trapezius in the affected limb, we demonstrate use of the Motion Browser, a novel visual analytics system to assess an adolescent with BPBP. Method  An 18-year-old female with extended upper trunk (C5–6–7) BPBP underwent bilateral upper extremity three-dimensional motion analysis with Motion Browser. Surface electromyography (EMG) from eight muscles in each limb which was recorded during six upper extremity movements, distinguishing between upper trapezius (UT) and lower trapezius (LT). The Motion Browser calculated active range of motion (AROM), compiled the EMG data into measures of muscle activity, and displayed the results in charts. Results  All movements, excluding shoulder abduction, had similar AROM in affected and unaffected limbs. In the unaffected limb, LT was more active in proximal movements of shoulder abduction, and shoulder external and internal rotations. In the affected limb, LT was more active in distal movements of forearm pronation and supination; UT was more active in shoulder abduction. Conclusion  In this female with BPBP, Motion Browser demonstrated that the native LT in the affected limb contributed to distal movements. Her results suggest that sacrificing her trapezius as a muscle or nerve donor may affect her distal functionality. Clinicians should exercise caution when considering nerve transfers in children with BPBP and consider individualized assessment of functionality before pursuing surgery.     

The measurement and lateral comparison of the peak torque caused by the fast abduction exercise of stretched upper extremities in normal and trained adults

The maximal torque effect of the middle portion of action of the deltoid muscle while raising an outstretched upper limb was measured from left and right sides of normal untrained young adults and of the same age elite athletes. Seventeen strongly right-handed untrained males and females and 10 elite tennis players were tested. All participants were required to raise (abduct) one arm (right and then left, or vice versa) as fast as possible with maximal amplitude while standing on an electronic platform scale which measured to 0.001 kg. An assumed force at the centre of mass of the entire upper limb was considered. The force consisted of two components, namely static weight force of the upper limb and a dynamic force component created by upward acceleration of the limb. Using regression equations and scaling methods the static weight of the upper limb was derived and combined with the dynamic component to produce the total force, applied to the centre of mass of the limb. The total force multiplied by the distance from the centre of mass to point of rotation of the limb equated to the torque produced by deltoid muscle. Using video system analyses the angle of abduction was measured for each individual exercise. The additional anthropometrical tests identified proportionality and body mass indices for each participant.

There was no significant difference in dynamic force and torque between left and right limb from the three groups. Sportsmen demonstrated greater lateral abduction when performing the exercise from the dominant side of the body. Sportsmen also demonstrated greater range of abduction, bigger dynamic force and torque on both sides in comparison to untrained adults. Remarkably, the absolute and relative length of arms of athletes were shorter in comparison to untrained males, but the radius of gyration from the stretched upper limb (from its centre of gravity to the shoulder joint) were greater. This phenomenon may be due to distal shifting of the gravity center of the entire upper limb in elite athletes, perhaps, because greater investment of the distal portion of the limb with skeletal muscle tissue.     

湖南、湖北、贵州土家族成人的体成分比较

体成分指的是身体脂肪、蛋白质、肌肉、水等含量在体质量中所占的百分比.各成分之间的合理比例,对于维持机体的正常运行,十分重要.因此身体成分的研究一直是国内外研究的热点.体成分存在着人种和民族差异,因此它也是人类学研究的重要课题.目前的研究主要集中在中国青少年、大学生的体成分及体成分与其他身体指标的关系,而对于土家族的研究...     

A control theoretic model of the forearm

In this paper, a control theoretic model of the forearm is developed and analyzed, and a computational method for predicting muscle activations necessary to generate specified motions is described. A detailed geometric model of the forearm kinematics, including the carrying angle and models of how the biceps and the supinator tendons wrap around the bones, is used. Also, including a dynamics model, the final model is a system of differential equations where the muscle activations play the role of control signals. Due to the large number of muscles, the problem of finding muscle activations is redundant, and this problem is solved by an optimization procedure. The computed muscle activations for ballistic movements clearly recaptures the triphasic ABC (Activation-Braking-Clamping) pattern. It is also transparent, from the muscle activation patterns, how the muscles cooperate and counteract in order to accomplish desired motions. A comparison with previously reported experimental data is included and the model predictions can be seen to be partially in agreement with the experimental data.     

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.     

The Obstacle-Set Method for Representing Muscle Paths in Musculoskeletal Models

A computational method is introduced for modeling the paths of muscles in the human body. The method is based on the premise that the resultant muscle force acts along the locus of the transverse cross-sectional centroids of the muscle. The path of the muscle is calculated by idealizing its centroid path as a frictionless elastic band, which moves freely over neighboring anatomical constraints such as bones and other muscles. The anatomical constraints, referred to as obstacles, are represented in the model by regular-shaped, rigid bodies such as spheres and cylinders. The obstacles, together with the muscle path, define an obstacle set. It is proposed that the path of any muscle can be modeled using one or more of the following four obstacle sets: single sphere, single cylinder, double cylinder, and sphere-capped cylinder. Assuming that the locus of the muscle centroids is known for an arbitrary joint configuration, the obstacle-set method can be used to calculate the path of the muscle for all other joint configurations. The obstacle-set method accounts not only for the interaction between a muscle and a neighboring anatomical constraint, but also for the way in which this interaction changes with joint configuration. Consequently, it is the only feasible method for representing the paths of muscles which cross joints with multiple degrees of freedom such as the deltoid at the shoulder.     

Soft-tissue flap coverage maximizes limb salvage after allograft bone extremity reconstruction

Limb salvage after extremity tumor ablation may include the use of allograft bone. The primary complication of this method is infection of the allograft, which can lead to limb loss in up to 50 percent of cases. The purpose of this study is to evaluate the efficacy of primary muscle flap coverage in the setting of allograft bone limb salvage surgery. This study is a prospective review of all patients with flap coverage of extremity allografts over the 10-year period 1991 to 2001. There were 20 patients (11 male and nine female patients) with an average age of 28 years (range, 6 to 72 years). Flap coverage was primary in 16 patients and delayed in four. Delayed coverage was performed for failed wounds that did not have a primary soft-tissue flap. Pathologic findings included osteosarcoma in nine patients, Ewing sarcoma in five patients, malignant fibrohistiocytoma in two patients, chondrosarcoma in two patients, synovial sarcoma in one patient, and leiomyosarcoma in one patient. Allograft reconstruction was performed for the upper extremity in 12 patients and for the lower extremity in eight patients. Flap reconstruction was accomplished with 20 pedicle flaps in 17 patients (latissimus dorsi, 12; gastrocnemius, four; soleus, three; and fasciocutaneous flap, one) and four free flaps (rectus abdominis, three; latissimus dorsi, one) in four patients. All pedicled flaps survived. There was one flap failure in the entire series, which was a free rectus abdominis flap. This case resulted in the only limb loss noted. The follow-up period ranged from 1 to 50 months (average, 12.35 months). At the time of final follow-up, three patients were dead of disease and 17 were alive with intact extremities. The overall limb salvage rate in the setting of bone allograft and soft-tissue flap coverage was 95 percent (19 of 20). Reoperation for bone-related complications was required in 50 percent (two of four) of cases receiving delayed flap coverage compared with 19 percent (three of 16) of patients with primary flap coverage (statistically not significant). The results of this study support the use of soft-tissue flap coverage for allograft limb reconstruction. In this series, no limb was lost in the setting of a viable flap. Reoperation was markedly reduced in the setting of primary flap coverage. Pedicled or microvascular transfer of well-vascularized muscle can be used to wrap the allograft and minimize devastating wound complications potentially leading to loss of allograft and limb.     

Transcranial Direct Current Stimulation Improves Ipsilateral Selective Muscle Activation in a Frequency Dependent Manner

Failure to suppress antagonist muscles can lead to movement dysfunction, such as the abnormal muscle synergies often seen in the upper limb after stroke. A neurophysiological surrogate of upper limb synergies, the selectivity ratio (SR), can be determined from the ratio of biceps brachii (BB) motor evoked potentials to transcranial magnetic stimulation prior to forearm pronation versus elbow flexion. Surprisingly, cathodal transcranial direct current stimulation (c-TDCS) over ipsilateral primary motor cortex (M1) reduces (i.e. improves) the SR in healthy adults, and chronic stroke patients. The ability to suppress antagonist muscles may be exacerbated at high movement rates. The aim of the present study was to investigate whether the selective muscle activation of the biceps brachii (BB) is dependent on altering frequency demands, and whether the c-tDCS improvement of SR is dependent on task frequency. Seventeen healthy participants performed repetitive isometric elbow flexion and forearm pronation at three rates, before and after c-tDCS or sham delivered to ipsilateral left M1. Ipsilateral c-tDCS improved the SR in a frequency dependent manner by selectively suppressing BB antagonist excitability. Our findings confirm that c-tDCS is an effective tool for improving selective muscle activation, and provide novel evidence for its efficacy at rates of movement where it is most likely to benefit task performance.     

Innervation zones of the upper and lower limb muscles estimated by using multichannel surface EMG

The distribution of innervation zones was investigated in 3 subjects for 17 muscles and 8 muscle groups in the upper and lower limb, by detecting bi-directional propagation of motor unit action potentials (MUAPs) with the multichannel surface electrode array. Clarification of the distribution of innervation zones depended on the ease in detecting the propagation of MUAPs and the actual scattering of innervation zones, which were closely related with muscle morphology with respect to the arrangements of muscle fibers. In muscles having fibers running parallel to each other, such as the biceps brachii, intrinsic hand muscles, vastus lateralis and medialis, tensor fasciae latae, peronei, soleus, tibialis anterior, and hypothenar muscles in the foot, it was relatively easy to detect the propagating MUAPs, and the innervation zones were distributed in a relatively narrow band around muscle belly. On the other hand, in muscles with a complicated structure including pinnation of muscle fibers, in-series muscle fibers and aponeurotic tissues, such as the deltoid, flexors and extensors in the forearm, rectus femoris, sartorius, hamstrings and gastrocnemius, it was more difficult to detect the propagating MUAPs and to identify the innervation zones, which were widely scattered or distributed in complex configurations. The distribution of the innervation zones clarified in the present study can be used to find the optimal location of electrodes in surface EMG recordings and of stimulus electrodes in the functional and therapeutic electrical stimulations. It may also be useful in motor point biopsy for diagnosis of neuromuscular diseases as well as in the botulinum toxin injection for the treatment of spasticity.     

运动想象结合常规康复训练对脑卒中偏瘫患者上肢肌力恢复的影响

目的:探讨运动想象结合常规康复训练对脑卒中偏瘫患者上肢肌力恢复的影响.方法:选择31例病情稳定的脑卒中偏瘫患者,随机分为治疗组(n=16)和对照组(n=15)两组,两组患者均进行常规康复训练,其中治疗组采用常规康复训练联合运动想象治疗,对照组只进行常规康复训练.两组常规康复训练时间相同.分别在治疗前、治疗后4周、8周进行患侧上肢肌力评定,观察患侧上肢肌力恢复情况.结果:治疗前两组各项评分差异无统计学意义.治疗8周后,两组患者肌力评级与各组治疗前比较差异均有统计学意义(P＜0.05),且两组肌力比较差异也有统计学意义(P＜0.05).结论:运动想象结合常规康复训练有利于脑卒中偏瘫患者上肢肌力的恢复,与以往的研究结果一致.     

Subscapularis function in gibbons and chimpanzees: Implications for interpretation of humeral head torsion in hominoids

Chimpanzees and gibbons, along with other hominoids, share a number of features in the morphology of their shoulders that have generally been associated with use of the upper limb in overhead postural and locomotor activities. These include the position and shape of the scapula, as well as the morphology of the proximal end of the humerus. Results of an electromyographic (EMG) analysis of shoulder muscle activity patterns indicate that these two species of hominoids also share broad similarities in shoulder muscle function during locomotion and voluntary movements. Differences do exist, however, in the activity patterns of subscapularis, a medial rotator of the arm. These differences mainly involve a greater participation by the gibbon subscapularis in free arm movements. This greater participation is characterized by earlier onset of muscle activity, higher amplitude of recruitment, and involvement of more of the total mass of the muscle. These differences in muscle recruitment suggest that the shoulder of gibbons differs from that of chimpanzees in some manner that necessitates the greater contribution of a medial rotator to the production of motion in the upper limb. I suggest that the low degree of humeral head torsion in gibbons, compared to that of other hominoids, gives their elbow a “lateral set” that must be overcome by the action of subscapularis during free arm movements. I propose that this modest degree of humeral head torsion in gibbons reflects a compromise between necessary changes caused by the repositioning of the scapula onto the dorsum of the thorax and the demands for extreme positioning of the elbow during brachiation. In addition, I suggest that the greater amount of torsion in the chimpanzee humerus is an accommodation to quadrupedal habits, and finally, that the high degree of torsion in human humeri is an independently acquired trait related to use of the upper limb as a manipulatory organ.