Biomechanical modeling and sensitivity analysis of bipedal running ability. I. Extant taxa

I used a simple mathematical model of the inverse dynamics of locomotion to estimate the minimum muscle masses required to maintain quasi-static equilibrium about the four main limb joints at mid-stance of fast running. Models of 10 extant taxa (a human, a kangaroo, two lizards, an alligator, and five birds) were analyzed in various bipedal poses to examine how anatomy, size, limb orientation, and other model parameters influence running ability. I examined how the muscle masses required for fast running compare to the muscle masses that are actually able to exert moments about the hip, knee, ankle, and toe joints, to see how support ability varies across the limb. I discuss the assumptions and limitations of the models, using sensitivity analysis to see how widely the results differed with feasible parameter input values. Even with a wide range of input values, the models validated the analysis procedure. Animals that are known to run bipedally were calculated as able to preserve quasi-static equilibrium about their hindlimb joints at mid-stance, whereas non-bipedal runners (iguanas and alligators) were recognized as having too little muscle mass to run quickly in bipedal poses. Thus, this modeling approach should be reliable for reconstructing running ability in extinct bipeds such as nonavian dinosaurs. The models also elucidated how key features are important for bipedal running capacity, such as limb orientation, muscle moment arms, muscle fascicle lengths, and body size. None of the animals modeled had extensor muscle masses acting about any one joint that were 7% or more of their body mass, which provides a reasonable limit for how much muscle mass is normally apportioned within a limb to act about a particular joint. The models consistently showed that a key biomechanical limit on running ability is the capacity of ankle extensors to generate sufficiently large joint moments. Additionally, the analysis reveals how large ratite birds remain excellent runners despite their larger size; they have apomorphically large extensor muscles with relatively high effective mechanical advantage. Finally, I reconstructed the evolution of running ability in the clade Reptilia, showing that the ancestors of extant birds likely were quite capable runners, even though they had already reduced key hip extensors such as M. caudofemoralis longus.     

Biarticular hip extensor and knee flexor muscle moment arms of the feline hindlimb

Moment arms are important for understanding muscular behavior and for calculating internal muscle forces in musculoskeletal simulations. Biarticular muscles cross two joints and have moment arms that depend on the angle of both joints the muscles cross. The tendon excursion method was used to measure the joint angle-dependence of hamstring (biceps femoris, semimembranosus and semitendinosus) moment arm magnitudes of the feline hindlimb at the knee and hip joints. Knee angle influenced hamstring moment arm magnitudes at the hip joint; compared to a flexed knee joint, the moment arm for semimembranosus posterior at the hip was at most 7.4 mm (25%) larger when the knee was extended. On average, hamstring moment arms at the hip increased by 4.9 mm when the knee was more extended. In contrast, moment arm magnitudes at the knee varied by less than 2.8 mm (mean=1.6 mm) for all hamstring muscles at the two hip joint angles tested. Thus, hamstring moment arms at the hip were dependent on knee position, while hamstring moment arms at the knee were not as strongly associated with relative hip position. Additionally, the feline hamstring muscle group had a larger mechanical advantage at the hip than at the knee joint.     

Biomechanical modeling and sensitivity analysis of bipedal running ability. II. Extinct taxa

Using an inverse dynamics biomechanical analysis that was previously validated for extant bipeds, I calculated the minimum amount of actively contracting hindlimb extensor muscle that would have been needed for rapid bipedal running in several extinct dinosaur taxa. I analyzed models of nine theropod dinosaurs (including birds) covering over five orders of magnitude in size. My results uphold previous findings that large theropods such as Tyrannosaurus could not run very quickly, whereas smaller theropods (including some extinct birds) were adept runners. Furthermore, my results strengthen the contention that many nonavian theropods, especially larger individuals, used fairly upright limb orientations, which would have reduced required muscular force, and hence muscle mass. Additional sensitivity analysis of muscle fascicle lengths, moment arms, and limb orientation supports these conclusions and points out directions for future research on the musculoskeletal limits on running ability. Although ankle extensor muscle support is shown to have been important for all taxa, the ability of hip extensor muscles to support the body appears to be a crucial limit for running capacity in larger taxa. I discuss what speeds were possible for different theropod dinosaurs, and how running ability evolved in an inverse relationship to body size in archosaurs.     

Quantitative functional anatomy of the lower limb with application to human gait

A scheme was developed to classify muscles according to their primary, secondary and tertiary functions, e.g. a muscle which produces primarily a flexion moment may also produce secondary abduction and tertiary internal rotation moments. The functions of muscles crossing the hip and knee joints were computed based upon the changing relative positions of joint centers and muscle origins and insertions during one gait cycle. The function of several of the major muscles crossing the hip and knee joints is reported for the different limb positions corresponding to normal gait. It was found that the amount of force necessary to produce a given moment about a joint was dependent upon the limb position. In addition, the muscle functions changed significantly with limb position. Electrical stimulation of muscles of a paralyzed subject gave qualitative support to the results.     

More than one way to be a giant: Convergence and disparity in the hip joints of saurischian dinosaurs

Saurischian dinosaurs evolved seven orders of magnitude in body mass, as well as a wide diversity of hip joint morphology and locomotor postures. The very largest saurischians possess incongruent bony hip joints, suggesting that large volumes of soft tissues mediated hip articulation. To understand the evolutionary trends and functional relationships between body size and hip anatomy of saurischians, we tested the relationships among discrete and continuous morphological characters using phylogenetically corrected regression. Giant theropods and sauropods convergently evolved highly cartilaginous hip joints by reducing supraacetabular ossifications, a condition unlike that in early dinosauromorphs. However, transitions in femoral and acetabular soft tissues indicate that large sauropods and theropods built their hip joints in fundamentally different ways. In sauropods, the femoral head possesses irregularly rugose subchondral surfaces for thick hyaline cartilage. Hip articulation was achieved primarily using the highly cartilaginous femoral head and the supraacetabular labrum on the acetabular ceiling. In contrast, theropods covered their femoral head and neck with thinner hyaline cartilage and maintained extensive articulation between the fibrocartilaginous femoral neck and the antitrochanter. These findings suggest that the hip joints of giant sauropods were built to sustain large compressive loads, whereas those of giant theropods experienced compression and shear forces.     

Power and work produced in different leg muscle groups when rising from a chair

The aim of this study was to determine the power output and work done by different muscle groups at the hip and knee joints during a rising movement, to be able to tell the degree of activation of the muscle groups and the relationship between concentric and eccentric work. Nine healthy male subjects rose from a chair with the seat at knee level. The moments of force about the hip and knee joints were calculated semidynamically. The power output (P) and work in the different muscle groups surrounding the joints was calculated as moment of force times joint angular velocity. Work was calculated as: work = f Pdt. The mean peak concentric power output was for the hip extensors 49.9 W, hip flexors 7.9 W and knee extensor 89.5 W. This power output corresponded to a net concentric work of 20.7 J, 1.0 J and 55.6 J, respectively. There was no concentric power output from the knee flexor muscles. Energy absorption through eccentric muscle action was produced by the hip extensors and hip flexors with a mean peak power output of 4.8 W and 7.4 W, respectively. It was concluded that during rising, the hip and knee muscles mainly worked concentrically and that the greatest power output and work were produced during concentric contraction of the knee and hip extensor muscles. There was however also a demand for eccentric work by the hip extensors as well as both concentric and eccentric work by the hip flexors. The knee flexor muscles were unloaded.     

The role of intersegmental dynamics during rapid limb oscillations

The interactive dynamic effects of muscular, inertial and gravitational moments on rapid, multi-segmented limb oscillations were studied. Using three-segment, rigid-body equations of motion, hip, knee and ankle intersegmental dynamics were calculated for the steady-state cycles of the paw-shake response in adult spinal cats. Hindlimb trajectories were filmed to obtain segmental kinematics, and myopotentials of flexors and extensors at each of the three joints were recorded synchronously with the ciné film. The segmental oscillations that emerged during the paw-shake response were a consequence of an interplay between active and passive musculotendinous forces, inertial forces, and gravity. During steady-state oscillations, the amplitudes of joint excursions, peak angular velocities, and peak angular accelerations increased monotonically and significantly in magnitude from the proximal joint (hip) to the most distal joint (ankle). In contrast to these kinematic relationships, the maximal values of net moments at the hip and knee were equal in magnitude, but of significantly lower magnitude than the large net moment at the ankle joint. At both the ankle and the knee, the flexor and extensor muscle moments were equal, but at the hip the magnitude of the peak flexor muscle moment was significantly greater than the extensor muscle moment. Muscle moments at the hip not only acted to counterbalance accelerations of the more distal segments, but also acted to maintain the postural orientation of the hindlimb. Large muscle moments at the knee functioned to counterbalance the large inertial moments generated by the large angular accelerations of the paw. At the ankle, the muscle moments dominated the generation of the paw accelerations. At the ankle and the knee, muscle moments controlled limb dynamics by slowing and reversing joint motions, and the active muscle forces contributing to ankle and knee moments were derived from lengthening of active musculotendinous units. In contrast to the more distal joints, the active muscles crossing the hip predominantly shortened as a result of the interplay among inertial forces and gravitational moments. The muscle function and kinetic data explain key features of the complex interactions that occur between central control mechanisms and multi-segmented, oscillating limb segments during the paw-shake response.     

Variation of rotation moment arms with hip flexion.

Excessive flexion and internal rotation of the hip is a common gait abnormality among individuals with cerebral palsy. The purpose of this study was to examine the influence of hip flexion on the rotational moment arms of the hip muscles. We hypothesized that flexion of the hip would increase internal rotation moment arms and decrease external rotation moment arms of the primary hip rotators. To test this hypothesis we measured rotational moment arms of the gluteus maximus (six compartments), gluteus medius (four compartments), gluteus minimus (three compartments) iliopsoas, piriformis, quadratus femoris, obturator internus, and obturator externus. Moment arms were measured at hip flexion angles of 0, 20, 45, 60, and 90 degrees in four cadavers. A three-dimensional computer model of the hip muscles was developed and compared to the experimental measurements. The experimental results and the computer model showed that the internal rotation moment arms of some muscles increase with flexion; the external rotation moment arms of other muscles decrease, and some muscles switch from external rotation to internal rotation as the hip is flexed. This trend toward internal rotation with hip flexion was apparent in 15 of the 18 muscle compartments we examined, suggesting that excessive hip flexion may exacerbate internal rotation of the hip. The gluteus maximus was found to have a large capacity for external rotation. Enhancing the activation of the gluteus maximus, a muscle that is frequently underactive in persons with cerebral palsy, may help correct excessive flexion and internal rotation of the hip.     

Modeling of the muscle/tendon excursions and moment arms in the thumb using the commercial software anybody

A biomechanical model of a thumb would be useful for exploring the mechanical loadings in the musculoskeletal system, which cannot be measured in vivo. The purpose of the current study is to develop a practical kinematic thumb model using the commercial software Anybody (Anybody Technology, Aalborg, Denmark), which includes real CT-scans of the bony sections and realistic tendon/muscle attachments on the bones. The thumb model consists of a trapezium, a metacarpal bone, a proximal and a distal phalanx. These four bony sections are linked via three joints, i.e., IP (interphalangeal), MP (metacarpophalangeal) and CMC (carpometacarpal) joints. Nine muscles were included in the proposed model. The theoretically calculated moment arms of the tendons are compared with the corresponding experimental data by Smutz et al. [1998. Mechanical advantage of the thumb muscles. J. Biomech. 31(6), 565–570]. The predicted muscle moment arms of the majority of the muscle/tendon units agree well with the experimental data in the entire range of motion. Close to the end of the motion range, the predicted moment arms of several muscles (i.e., ADPt and ADPo (transverse and oblique heads of the adductor pollicis, respectively) muscles for CMC abduction/adduction and ADPt and FPB (flexor pollicis brevis) muscle for MP extension/flexion) deviate from the experimental data. The predicted moment potentials for all muscles are consistent with the experimental data. The findings thus suggest that, in a biomechanical model of the thumb, the mechanical functions of muscle–tendon units with small physiological cross-sectional areas (PCSAs) can be well represented using single strings, while those with large PCSAs (flat-wide attachments, e.g., ADPt and ADPo) can be represented by the averaged excursions of two strings. Our results show that the tendons with large PCSAs can be well represented biomechanically using the proposed approach in the major range of motion.     

Gear ratios at the limb joints of jumping dogs

An increase in gear ratio of the limb extensor muscles during joint extension has been suggested to be a mechanism that facilitates optimal power production by skeletal muscles. The objectives of this study were to: (1) measure gear ratios at the wrist, elbow, shoulder, ankle, knee, and hip joints of jumping dogs, (2) compute the work performed by each of these joints, and (3) measure muscle shortening velocity for a joint exhibiting an increasing gear ratio during joint extension. The gear ratio out-lever was computed by dividing the ground reaction force (GRF) moment by the GRF, whereas the in-lever was directly measured as the perpendicular distance from the joint center to the line of action of the extensor muscle. In addition, changes in fascicle length were measured from the vastus lateralis muscle using sonomicrometry. Of the joints examined, only the gear ratios at the shoulder and knee joints increased during jumping in a manner that could facilitate peak power production of actively shortening muscles. The vastus lateralis was found to shorten at an average velocity of 3.20 muscle lengths per second. This is similar to estimates of shortening velocity that produce peak muscular power in mammals the size of dogs. Additionally, the knee extensors were found to produce a large proportion (26.6%) of the positive external work of the limbs. These observations suggest that dynamic gearing in jumping dogs may allow the extensor muscles of the knee joint to shorten in a way that maximizes their power production.     

Muscular contributions to hip and knee extension during the single limb stance phase of normal gait: a framework for investigating the causes of crouch gait

Crouch gait, a troublesome movement abnormality among persons with cerebral palsy, is characterized by excessive flexion of the hips and knees during stance. Treatment of crouch gait is challenging, at present, because the factors that contribute to hip and knee extension during normal gait are not well understood, and because the potential of individual muscles to produce flexion or extension of the joints during stance is unknown. This study analyzed a three-dimensional, muscle-actuated dynamic simulation of walking to quantify the angular accelerations of the hip and knee induced by muscles during normal gait, and to rank the potential of the muscles to alter motions of these joints. Examination of the muscle actions during single limb stance showed that the gluteus maximus, vasti, and soleus make substantial contributions to hip and knee extension during normal gait. Per unit force, the gluteus maximus had greater potential than the vasti to accelerate the knee toward extension. These data suggest that weak hip extensors, knee extensors, or ankle plantar flexors may contribute to crouch gait, and strengthening these muscles--particularly gluteus maximus--may improve hip and knee extension. Abnormal forces generated by the iliopsoas or adductors may also contribute to crouch gait, as our analysis showed that these muscles have the potential to accelerate the hip and knee toward flexion. This work emphasizes the need to consider how muscular forces contribute to multijoint movements when attempting to identify the causes of abnormal gait.     

Estimation of the effective static moment arms of the tendons in the index finger extensor mechanism

A novel technique to estimate the contribution of finger extensor tendons to joint moment generation was proposed. Effective static moment arms (ESMAs), which represent the net effects of the tendon force on joint moments in static finger postures, were estimated for the 4 degrees of freedom (DOFs) in the index finger. Specifically, the ESMAs for the five tendons contributing to the finger extensor apparatus were estimated by directly correlating the applied tendon force to the measured resultant joint moments in cadaveric hand specimens. Repeated measures analysis of variance revealed that the finger posture, specifically interphalangeal joint angles, had significant effects on the measured ESMA values in 7 out of 20 conditions (four DOFs for each of the five muscles). Extensor digitorum communis and extensor indicis proprius tendons were found to have greater MCP ESMA values when IP joints are flexed, whereas abduction ESMAs of all muscles except extensor digitorum profundus were mainly affected by MCP flexion. The ESMAs were generally smaller than the moment arms estimated in previous studies that employed kinematic measurement techniques. Tendon force distribution within the extensor hood and dissipation into adjacent structures are believed to contribute to the joint moment reductions, which result in smaller ESMA values.     

Contribution of mono- and biarticular muscles to extending knee joint moments in runners and cyclists.

Motor actions are governed by coordinated activation of mono- and biarticular muscles. This study considered differences in mono- and biarticular knee extensors between runners and cyclists in the context of adaptations to task-specific movement requirements. Two hypotheses were tested: 1) the length-at-use hypothesis, which is that muscle adapts to have it operate around optimal length; and 2) the contraction-mode hypothesis, which is that eccentrically active muscles prefer to operate on the ascending limb of the length-force curve. Ten runners and ten cyclists performed maximal, isometric knee extensions on a dynamometer at five knee and four hip joint angles. This approach allowed the separation of the contribution of mono- and biarticular extensors. Three major differences occurred: 1) compared with runners, monoarticular extensors of cyclists reach optimal length at larger muscle length; 2) in runners, optimal length of the biarticular extensor is shifted to larger lengths; and 3) the moment generated by monoarticular extensor was larger in cyclists. Mono- and biarticular extensors respond to different adaptation triggers in runners and cyclists. Monoarticular muscles seem to adapt to the length-at-use, whereas biarticular muscles were found to be sensitive to the contraction-mode hypothesis.     

Propulsive adaptation to changing gait speed

Understanding propulsion and adaptation to speed requirements is important in determining appropriate therapies for gait disorders. We hypothesize that adaptations for changing speed requirements occur primarily at the hip. The slow, normal and fast gait of 24 healthy young subjects was analyzed. The linear power was analyzed at the hip joint. The anterior-posterior and vertical induced accelerations of the hip were also determined. Linear power and anterior-posterior-induced acceleration (IA) analyses of the hip reveal that the lower limb joint's moments contribute to body forward propulsion primarily during late swing and early stance. Propulsive adaptations to speed changes occur primarily at the hip and secondarily at the ankle. These analyses show that hip muscles, particularly the hip extensors, are critical to propulsion. They also show that ankle function is primarily for support, but is important to propulsion, especially at slow speeds.     

Locomotion in ornithischian dinosaurs: an assessment using three‐dimensional computational modelling

Ornithischian dinosaurs were primitively bipedal with forelimbs modified for grasping, but quadrupedalism evolved in the clade on at least three occasions independently. Outside of Ornithischia, quadrupedality from bipedal ancestors has only evolved on two other occasions, making this one of the rarest locomotory transitions in tetrapod evolutionary history. The osteological and myological changes associated with these transitions have only recently been documented, and the biomechanical consequences of these changes remain to be examined. Here, we review previous approaches to understanding locomotion in extinct animals, which can be broadly split into form–function approaches using analogy based on extant animals, limb‐bone scaling, and computational approaches. We then carry out the first systematic attempt to quantify changes in locomotor muscle function in bipedal and quadrupedal ornithischian dinosaurs. Using three‐dimensional computational modelling of the major pelvic locomotor muscle moment arms, we examine similarities and differences among individual taxa, between quadrupedal and bipedal taxa, and among taxa representing the three major ornithischian lineages (Thyreophora, Ornithopoda, Marginocephalia). Our results suggest that the ceratopsid Chasmosaurus and the ornithopod Hypsilophodon have relatively low moment arms for most muscles and most functions, perhaps suggesting poor locomotor performance in these taxa. Quadrupeds have higher abductor moment arms than bipeds, which we suggest is due to the overall wider bodies of the quadrupeds modelled. A peak in extensor moment arms at more extended hip angles and lower medial rotator moment arms in quadrupeds than in bipeds may be due to a more columnar hindlimb and loss of medial rotation as a form of lateral limb support in quadrupeds. We are not able to identify trends in moment arm evolution across Ornithischia as a whole, suggesting that the bipedal ancestry of ornithischians did not constrain the development of quadrupedal locomotion via a limited number of functional pathways. Functional anatomy appears to have had a greater effect on moment arms than phylogeny, and the differences identified between individual taxa and individual clades may relate to differences in locomotor performance required for living in different environments or for clade‐specific behaviours.     

A non-invasive protocol to determine the personalized moment arms of knee and ankle muscles

One difficulty that comes with predicting muscular forces is the accuracy of experimental data, particularly the assessment of muscle moment arms with respect to each joint rotation axis. This paper presents a non-invasive experimental protocol to obtain the personalized muscle moment arms with respect to the ankle and knee joints. A specific pointer is used by a specialist of lower limb anatomy in order to define the local portion of the line of action of the different muscles closed to the joint on the standing subject's lower limb. With this pointer, the three-dimensional coordinates of several points representing the line of action of 12 ankle and knee muscles are collected by a Motion Analysis system. The collection is done five times by the same operator and one time by two different operators. From this data, the intra and inter operator repeatabilities are tested. Relative (ICC) and absolute (SEM) reliabilities are determined in order to evaluate the intra operator repeatability of this non-invasive protocol. The ICC values obtained are higher than 0.91 for 10 among 12 muscles. The intra operator repeatability is thus confirmed. From the records realized by the two operators, the differences are negligible. Thus, the inter operator repeatability is also confirmed. The moments arms obtained using this non-invasive experimental protocol are compared with those calculated from origin and insertion points reported in the literature, according to the work of Whites, Pierrynowskis and Kepples, respectively. The estimations obtained using the non-invasive experimental protocol are found, for some muscles, more realistic than those calculated using the literature data and are always coherent with the role of the muscles described in anatomical books.     

Lower extremity biomechanics during a regular and counterbalanced squat

If the efficiency of human movement patterns could be improved using exercise, this could lead to more effective musculoskeletal disease-injury prevention and rehabilitation programs. It has been suggested that an efficient squat movement pattern emphasizes the use of the large hip extensors instead of the smaller knee extensors. The purpose of this study was to determine whether a counterbalanced squat (CBS) could produce a more hip-dominant and less knee-dominant squat movement pattern as compared with a regular squat (RS). There were 31 recreationally trained college-aged participants (15 male, 16 female) who performed 10 squats (5 CBS and 5 RS), while segment kinematics, ground reaction forces, and muscle (gluteus maximus [GM], quadriceps, hamstrings) electromyographic (EMG) activations were recorded. Peak sagittal plane net joint moments and joint ranges of motion at the hip, knee, and ankle joints along with peak and integrated EMG activation levels for all 3 muscles were compared using analysis of variance (squat type × sex). The results revealed that the CBS increased the hip joint moment and GM activation, while it decreased the knee joint moment and quadriceps activation as compared with the RS. Therefore, the CBS produces a more hip-dominant and less knee-dominant squat movement pattern and could be used in exercise programs aimed at producing more hip-dominant movement patterns.     

An experimental and analytic study of the dynamic properties of the human leg

Certain problems in the dynamics of the human body are characterized by large displacements of the parts of the body compared to the deformations of the tissues themselves. In such problems, it is convenient to think of the human body as a chain of rigid links, representing the anatomical segments, with joints between the rigid links representing the articulations of the human body.

Skeletal muscles are capable of creating torques at the joints of the body. The joint torques of the rigid link model should portray the static strength of skeletal muscle, the degradation of muscle strength with rate of shortening, the feedback control of the stretch reflex, and the viscoelastic properties of the muscles, tendons, and joint capsules.

In this study a sinusoidal test is performed upon the knee joints of nine subjects. The increment of knee moment required to flex and extend the knee slightly for various conditions of knee angle, knee angular velocity, and steady knee moment is measured. The hip angle is maintained constant. After the effects of leg inertia are removed, the resulting data are shown to obey a Maxwell fluid model in which the model coefficients depend upon the absolute value of the knee moment.     

Dimensions and moment arms of the hind- and forelimb muscles of common chimpanzees (Pan troglodytes).

This paper supplies quantitative data on the hind- and forelimb musculature of common chimpanzees (Pan troglodytes) and calculates maximum joint moments of force as a contribution to a better understanding of the differences between chimpanzee and human locomotion. We dissected three chimpanzees, and recorded muscle mass, fascicle length, and physiological cross-sectional area (PCSA). We also obtained flexion/extension moment arms of the major muscles about the limb joints. We find that in the hindlimb, chimpanzees possess longer fascicles in most muscles but smaller PCSAs than are predicted for humans of equal body mass, suggesting that the adaptive emphasis in chimpanzees is on joint mobility at the expense of tension production. In common chimpanzee bipedalism, both hips and knees are significantly more flexed than in humans, necessitating muscles capable of exerting larger moments at the joints for the same ground force. However, we find that when subject to the same stresses, chimpanzee hindlimb muscles provide far smaller moments at the joints than humans, particularly the quadriceps and plantar flexors. In contrast, all forelimb muscle masses, fascicle lengths, and PCSAs are smaller in humans than in chimpanzees, reflecting the use of the forelimbs in chimpanzee, but not human, locomotion. When subject to the same stresses, chimpanzee forelimb muscles provide larger moments at the joints than humans, presumably because of the demands on the forelimbs during locomotion. These differences in muscle architecture and function help to explain why chimpanzees are restricted in their ability to walk, and particularly to run bipedally.     

A musculoskeletal model customized for squatting task

Most musculoskeletal models (MSKM) are designed to evaluate gait and running, which have limited range of motion (ROM). The purpose of this study was to examine the effect of wrapping surfaces (WS) at the knee and hip joints in a MSKM, on the muscle moment arms (MA) and activations during squatting. The MSKM was then customized by changing parameters of the original WS and by implementing additional WS. The WS prevent muscles from crossing into the bones, providing realistic muscle MA for large ROM. The modified MSKM is suitable for analysis up to 138° hip and 145° knee flexions.