Fibre operating lengths of human lower limb muscles during walking

Muscles actuate movement by generating forces. The forces generated by muscles are highly dependent on their fibre lengths, yet it is difficult to measure the lengths over which muscle fibres operate during movement. We combined experimental measurements of joint angles and muscle activation patterns during walking with a musculoskeletal model that captures the relationships between muscle fibre lengths, joint angles and muscle activations for muscles of the lower limb. We used this musculoskeletal model to produce a simulation of muscle-tendon dynamics during walking and calculated fibre operating lengths (i.e. the length of muscle fibres relative to their optimal fibre length) for 17 lower limb muscles. Our results indicate that when musculotendon compliance is low, the muscle fibre operating length is determined predominantly by the joint angles and muscle moment arms. If musculotendon compliance is high, muscle fibre operating length is more dependent on activation level and force-length-velocity effects. We found that muscles operate on multiple limbs of the force-length curve (i.e. ascending, plateau and descending limbs) during the gait cycle, but are active within a smaller portion of their total operating range.     

A computationally efficient strategy to estimate muscle forces in a finite element musculoskeletal model of the lower limb

Concurrent multiscale simulation strategies are required in computational biomechanics to study the interdependence between body scales. However, detailed finite element models rarely include muscle recruitment due to the computational burden of both the finite element method and the optimization strategies widely used to estimate muscle forces. The aim of this study was twofold: first, to develop a computationally efficient muscle force prediction strategy based on proportional-integral-derivative (PID) controllers to track gait and chair rise experimental joint motion with a finite element musculoskeletal model of the lower limb, including a deformable knee representation with 12 degrees of freedom; and, second, to demonstrate that the inclusion of joint-level deformability affects muscle force estimation by using two different knee models and comparing muscle forces between the two solutions. The PID control strategy tracked experimental hip, knee, and ankle flexion/extension with root mean square errors below 1°, and estimated muscle, contact and ligament forces in good agreement with previous results and electromyography signals. Differences up to 11% and 20% in the vasti and biceps femoris forces, respectively, were observed between the two knee models, which might be attributed to a combination of differing joint contact geometry, ligament behavior, joint kinematics, and muscle moment arms. The tracking strategy developed in this study addressed the inevitable tradeoff between computational cost and model detail in musculoskeletal simulations and can be used with finite element musculoskeletal models to efficiently estimate the interdependence between muscle forces and tissue deformation.     

Preparatory co-activation of the ankle muscles may prevent ankle inversion injuries

Ankle inversion sprains are the most frequent acute musculoskeletal injuries occurring in physical activity. Interventions that retrain muscle coordination have helped rehabilitate injured ankles, but it is unclear which muscle coordination strategies, if any, can prevent ankle sprains. The purpose of this study was to determine whether coordinated activity of the ankle muscles could prevent excessive ankle inversion during a simulated landing on a 30° incline. We used a set of musculoskeletal simulations to evaluate the efficacy of two strategies for coordinating the ankle evertor and invertor muscles during simulated landing scenarios: planned co-activation and stretch reflex activation with physiologic latency (60-ms delay). A full-body musculoskeletal model of landing was used to generate simulations of a subject dropping onto an inclined surface with each coordination condition. Within each condition, the intensity of evertor and invertor co-activity or stretch reflexes were varied systematically. The simulations revealed that strong preparatory co-activation of the ankle evertors and invertors prior to ground contact prevented ankle inversion from exceeding injury thresholds by rapidly generating eversion moments after initial contact. Conversely, stretch reflexes were too slow to generate eversion moments before the simulations reached the threshold for inversion injury. These results suggest that training interventions to protect the ankle should focus on stiffening the ankle with muscle co-activation prior to landing. The musculoskeletal models, controllers, software, and simulation results are freely available online at http://simtk.org/home/ankle-sprains, enabling others to reproduce the results and explore new injury scenarios and interventions.     

To what extent is joint and muscle mechanics predicted by musculoskeletal models sensitive to soft tissue artefacts?

Musculoskeletal models are widely used to estimate joint kinematics, intersegmental loads, and muscle and joint contact forces during movement. These estimates can be heavily affected by the soft tissue artefact (STA) when input positional data are obtained using stereophotogrammetry, but this aspect has not yet been fully characterised for muscle and joint forces. This study aims to assess the sensitivity to the STA of three open-source musculoskeletal models, implemented in OpenSim.A baseline dataset of marker trajectories was created for each model from experimental data of one healthy volunteer. Five hundred STA realizations were then statistically generated using a marker-dependent model of the pelvis and lower limb artefact and added to the baseline data. The STA׳s impact on the musculoskeletal model estimates was finally quantified using a Monte Carlo analysis.The modelled STA distributions were in line with the literature. Observed output variations were comparable across the three models, and sensitivity to the STA was evident for most investigated quantities. Shape, magnitude and timing of the joint angle and moment time histories were not significantly affected throughout the entire gait cycle, whereas magnitude variations were observed for muscle and joint forces. Ranges of contact force variations differed between joints, with hip variations up to 1.8 times body weight observed. Variations of more than 30% were observed for some of the muscle forces.In conclusion, musculoskeletal simulations using stereophotogrammetry may be safely run when only interested in overall output patterns. Caution should be paid when more accurate estimated values are needed.     

The interaction of muscle moment arm,knee laxity,and torque in a multi-scale musculoskeletal model of the lower limb

IntroductionMusculoskeletal modeling allows insight into the interaction of muscle force and knee joint kinematics that cannot be measured in the laboratory. However, musculoskeletal models of the lower extremity commonly use simplified representations of the knee that may limit analyses of the interaction between muscle forces and joint kinematics. The goal of this research was to demonstrate how muscle forces alter knee kinematics and consequently muscle moment arms and joint torque in a musculoskeletal model of the lower limb that includes a deformable representation of the knee.MethodsTwo musculoskeletal models of the lower limb including specimen-specific articular geometries and ligament deformability at the knee were built in a finite element framework and calibrated to match mean isometric torque data collected from 12 healthy subjects. Muscle moment arms were compared between simulations of passive knee flexion and maximum isometric knee extension and flexion. In addition, isometric torque results were compared with predictions using simplified knee models in which the deformability of the knee was removed and the kinematics at the joint were prescribed for all degrees of freedom.ResultsPeak isometric torque estimated with a deformable knee representation occurred between 45° and 60° in extension, and 45° in flexion. The maximum isometric flexion torques generated by the models with deformable ligaments were 14.6% and 17.9% larger than those generated by the models with prescribed kinematics; by contrast, the maximum isometric extension torques generated by the models were similar. The change in hamstrings moment arms during isometric flexion was greater than that of the quadriceps during isometric extension (a mean RMS difference of 9.8 mm compared to 2.9 mm, respectively).DiscussionThe large changes in the moment arms of the hamstrings, when activated in a model with deformable ligaments, resulted in changes to flexion torque. When simulating human motion, the inclusion of a deformable joint in a multi-scale musculoskeletal finite element model of the lower limb may preserve the realistic interaction of muscle force with knee kinematics and torque.     

Force enhancement above the initial isometric force on the descending limb of the force-length relationship

Edman et al. (J. General Physiol. 80 (1982) 769) observed in single fibres of frog that the steady-state forces following active fibre stretch were greater than the purely isometric force obtained at the length from which the stretch was initiated. Operating on the descending limb of the force-length relationship, such a result can only be explained within the framework of the sarcomere length non-uniformity theory, if some fibre segments shortened during the fibre stretch. However, such a result was not found, leaving Edman's observation unexplained. Force enhancement above the initial isometric force has not been investigated systematically in whole muscle, and therefore it is not known whether this property is also part of whole muscle mechanics. The purpose of this study was to test if the steady-state forces following active stretch of cat semitendinosus were greater than the corresponding purely isometric forces at the muscle length from which the stretch was started. Cat semitendinosus was stretched by various amounts on the descending limb of the force-length relationship, and the steady-state forces following these stretches were compared to the corresponding isometric forces at the initial and final muscle lengths. In 109 of 131 tests, the steady-state forces following stretching were greater than the isometric forces at the initial muscle lengths. Force enhancement increased with increasing amounts of stretching, and force enhancement above the initial isometric force was more likely to occur following stretches of great compared to small amplitude. Passive forces following active muscle stretching were often significantly greater than the passive forces at the same muscle length following an isometric contraction or a passive stretching of the muscle. This observation was made consistently at the longest muscle lengths tested. It appears, therefore, that there is a passive force that accounts for part of the force enhancement above the isometric force at the initial muscle length, and that provides increased passive force when a muscle is actively, rather than passively, stretched at long muscle lengths. We conclude that cat semitendinosus demonstrates steady-state force enhancement above the corresponding purely isometric force at the initial muscle length on the descending limb of the force-length relationship for many contractile conditions, and that a unique, and so far undetected, passive, parallel element contributes to this force enhancement, particularly at long muscle lengths where muscle is assumed to be most vulnerable to injuries associated with sarcomere length instability.     

Development of a finite element musculoskeletal model with the ability to predict contractions of three-dimensional muscles

Representation of realistic muscle geometries is needed for systematic biomechanical simulation of musculoskeletal systems. Most of the previous musculoskeletal models are based on multibody dynamics simulation with muscles simplified as one-dimensional (1D) line-segments without accounting for the large muscle attachment areas, spatial fibre alignment within muscles and contact and wrapping between muscles and surrounding tissues. In previous musculoskeletal models with three-dimensional (3D) muscles, contractions of muscles were among the inputs rather than calculated, which hampers the predictive capability of these models. To address these issues, a finite element musculoskeletal model with the ability to predict contractions of 3D muscles was developed. Muscles with realistic 3D geometry, spatial muscle fibre alignment and muscle-muscle and muscle-bone interactions were accounted for. Active contractile stresses of the 3D muscles were determined through an efficient optimization approach based on the measured kinematics of the lower extremity and ground force during gait. This model also provided stresses and strains of muscles and contact mechanics of the muscle-muscle and muscle-bone interactions. The total contact force of the knee predicted by the model corresponded well to the in vivo measurement. Contact and wrapping between muscles and surrounding tissues were evident, demonstrating the need to consider 3D contact models of muscles. This modelling framework serves as the methodological basis for developing musculoskeletal modelling systems in finite element method incorporating 3D deformable contact models of muscles, joints, ligaments and bones.     

Subject-specific musculoskeletal modelling in patients before and after total hip arthroplasty*

The goal of this study was to define the effect on hip contact forces of including subject-specific moment generating capacity in the musculoskeletal model by scaling isometric muscle strength and by including geometrical information in control subjects, hip osteoarthritis and total hip arthroplasty patients. Scaling based on dynamometer measurements decreased the strength of all flexor and abductor muscles. This resulted in a model that lacked the capacity to generate joint moments required during functional activities. Scaling muscle forces based on functional activities and inclusion of MRI-based geometrical detail did not compromise the model strength and resulted in hip contact forces comparable to previously reported measured contact forces.     

Analysis of musculoskeletal loading in an index finger during tapping

Since musculoskeletal disorders of the upper extremities are believed to be associated with repetitive excessive muscle force production in the hands, understanding the time-dependent muscle forces during key tapping is essential for exploring the mechanisms of disease initiation and development. In the current study, we have simulated the time-dependent dynamic loading in the muscle/tendons in an index finger during tapping. The index finger model is developed using a commercial software package AnyBody, and it contains seven muscle/tendons that connect the three phalangeal finger sections. Our simulations indicate that the ratios of the maximal forces in flexor digitorum superficialis (FS) and flexor digitorum profundus (FP) tendons to the maximal force at the fingertip are 0.95 and 2.9, respectively, which agree well with recently published experimental data. The time sequence of the finger muscle activation predicted in the current study is consistent with the EMG data in the literature. The proposed model will be useful for bioengineers and ergonomic designers to improve keyboard design minimizing musculoskeletal loadings in the fingers.     

Computational modelling of muscle fibre operating ranges in the hindlimb of a small ground bird (Eudromia elegans), with implications for modelling locomotion in extinct species

The arrangement and physiology of muscle fibres can strongly influence musculoskeletal function and whole-organismal performance. However, experimental investigation of muscle function during in vivo activity is typically limited to relatively few muscles in a given system. Computational models and simulations of the musculoskeletal system can partly overcome these limitations, by exploring the dynamics of muscles, tendons and other tissues in a robust and quantitative fashion. Here, a high-fidelity, 26-degree-of-freedom musculoskeletal model was developed of the hindlimb of a small ground bird, the elegant-crested tinamou (Eudromia elegans, ~550 g), including all the major muscles of the limb (36 actuators per leg). The model was integrated with biplanar fluoroscopy (XROMM) and forceplate data for walking and running, where dynamic optimization was used to estimate muscle excitations and fibre length changes throughout both gaits. Following this, a series of static simulations over the total range of physiological limb postures were performed, to circumscribe the bounds of possible variation in fibre length. During gait, fibre lengths for all muscles remained between 0.5 to 1.21 times optimal fibre length, but operated mostly on the ascending limb and plateau of the active force-length curve, a result that parallels previous experimental findings for birds, humans and other species. However, the ranges of fibre length varied considerably among individual muscles, especially when considered across the total possible range of joint excursion. Net length change of muscle–tendon units was mostly less than optimal fibre length, sometimes markedly so, suggesting that approaches that use muscle–tendon length change to estimate optimal fibre length in extinct species are likely underestimating this important parameter for many muscles. The results of this study clarify and broaden understanding of muscle function in extant animals, and can help refine approaches used to study extinct species.     

Biophysical stimuli induced by passive movements compensate for lack of skeletal muscle during embryonic skeletogenesis

In genetically modified mice with abnormal skeletal muscle development, bones and joints are differentially affected by the lack of skeletal muscle. We hypothesise that unequal levels of biophysical stimuli in the developing humerus and femur can explain the differential effects on these rudiments when muscle is absent. We find that the expression patterns of four mechanosensitive genes important for endochondral ossification are differentially affected in muscleless limb mutants, with more extreme changes in the expression in the humerus than in the femur. Using finite element analysis, we show that the biophysical stimuli induced by muscle forces are similar in the humerus and femur, implying that the removal of muscle contractile forces should, in theory, affect the rudiments equally. However, simulations in which a displacement was applied to the end of the limb, such as could be caused in muscleless mice by movements of the mother or normal littermates, predicted higher biophysical stimuli in the femur than in the humerus. Stimuli induced by limb movement were much higher than those induced by the direct application of muscle forces, and we propose that movements of limbs caused by muscle contractions, rather than the direct application of muscle forces, provide the main mechanical stimuli for normal skeletal development. In muscleless mice, passive movement induces unequal biophysical stimuli in the humerus and femur, providing an explanation for the differential effects seen in these mice. The significance of these results is that forces originating external to the embryo may contribute to the initiation and progression of skeletal development when muscle development is abnormal.     

Knee muscle forces during walking and running in patellofemoral pain patients and pain-free controls

One proposed mechanism of patellofemoral pain, increased stress in the joint, is dependent on forces generated by the quadriceps muscles. Describing causal relationships between muscle forces, tissue stresses, and pain is difficult due to the inability to directly measure these variables in vivo. The purpose of this study was to estimate quadriceps forces during walking and running in a group of male and female patients with patellofemoral pain (n=27, 16 female; 11 male) and compare these to pain-free controls (n=16, 8 female; 8 male). Subjects walked and ran at self-selected speeds in a gait laboratory. Lower limb kinematics and electromyography (EMG) data were input to an EMG-driven musculoskeletal model of the knee, which was scaled and calibrated to each individual to estimate forces in 10 muscles surrounding the joint. Compared to controls, the patellofemoral pain group had greater co-contraction of quadriceps and hamstrings (p=0.025) and greater normalized muscle forces during walking, even though the net knee moment was similar between groups. Muscle forces during running were similar between groups, but the net knee extension moment was less in the patellofemoral pain group compared to controls. Females displayed 30–50% greater normalized hamstring and gastrocnemius muscle forces during both walking and running compared to males (p<0.05). These results suggest that some patellofemoral pain patients might experience greater joint contact forces and joint stresses than pain-free subjects. The muscle force data are available as supplementary material.     

A multibody knee model with discrete cartilage prediction of tibio-femoral contact mechanics

Combining musculoskeletal simulations with anatomical joint models capable of predicting cartilage contact mechanics would provide a valuable tool for studying the relationships between muscle force and cartilage loading. As a step towards producing multibody musculoskeletal models that include representation of cartilage tissue mechanics, this research developed a subject-specific multibody knee model that represented the tibia plateau cartilage as discrete rigid bodies that interacted with the femur through deformable contacts. Parameters for the compliant contact law were derived using three methods: (1) simplified Hertzian contact theory, (2) simplified elastic foundation contact theory and (3) parameter optimisation from a finite element (FE) solution. The contact parameters and contact friction were evaluated during a simulated walk in a virtual dynamic knee simulator, and the resulting kinematics were compared with measured in vitro kinematics. The effects on predicted contact pressures and cartilage–bone interface shear forces during the simulated walk were also evaluated. The compliant contact stiffness parameters had a statistically significant effect on predicted contact pressures as well as all tibio-femoral motions except flexion–extension. The contact friction was not statistically significant to contact pressures, but was statistically significant to medial–lateral translation and all rotations except flexion–extension. The magnitude of kinematic differences between model formulations was relatively small, but contact pressure predictions were sensitive to model formulation. The developed multibody knee model was computationally efficient and had a computation time 283 times faster than a FE simulation using the same geometries and boundary conditions.     

Modeling the effects of musculoskeletal geometry on scapulohumeral muscle moment arms and lines of action

Abstract

Biomechanical investigations examining shoulder function commonly observe a high degree of inter-individual variability in muscle activity and kinematic patterns during static and dynamic upper extremity exertions. Substantial differences in musculoskeletal geometry between individuals can alter muscle moment arms and lines of action that, theoretically, alter muscle activity and shoulder kinematics. The purposes of this research were to: (i) quantify model-predicted functional roles (moment arms, lines of action) of the scapulohumeral muscles, (ii) compare model predictions to experimental data in the literature, and (iii) evaluate sensitivity of muscle functional roles due to changes in muscle attachment locations using probabilistic modeling. Monte Carlo simulations were performed to iteratively adjust muscle attachment locations at the clavicle, scapula, and humerus of the Delft Shoulder and Elbow Model in OpenSim. Muscle moment arms and lines of action were quantified throughout arm elevation in the scapular plane. In general, model-predicted moment arms agreed well with the reviewed literature; however, notable inconsistencies were observed when comparing lines of action. Variability in moment arms and lines of action were muscle-specific, with 2 standard deviations in moment arm and line of actions as high as 25.8?mm and 18.8° for some muscles, respectively. Moment arms were particularly sensitive to changes in attachment site closest to the joint centre. Variations in muscle functional roles due to differences in musculoskeletal geometry are expected to require different muscle activity and movement patterns for upper extremity exertions.     

Modulation of passive force in single skeletal muscle fibres

In this study, we investigated the effects of activation and stretch on the passive force-sarcomere length relationship in skeletal muscle. Single fibres from the lumbrical muscle of frogs were placed at varying sarcomere lengths on the descending limb of the force-sarcomere length relationship, and tetanic contractions, active stretches and passive stretches (amplitudes of ca 10% of fibre length at a speed of 40% fibre length/s) were performed. The passive forces following stretch of an activated fibre were higher than the forces measured after isometric contractions or after stretches of a passive fibre at the corresponding sarcomere length. This effect was more pronounced at increased sarcomere lengths, and the passive force-sarcomere length relationship following active stretch was shifted upwards on the force axis compared with the corresponding relationship obtained following isometric contractions or passive stretches. These results provide strong evidence for an increase in passive force that is mediated by a length-dependent combination of stretch and activation, while activation or stretch alone does not produce this effect. Based on these results and recently published findings of the effects of Ca2+ on titin stiffness, we propose that the observed increase in passive force is caused by the molecular spring titin.     

Non-rigid deformation to include subject-specific detail in musculoskeletal models of CP children with proximal femoral deformity and its effect on muscle and contact forces during gait

To account for proximal femoral deformities in children with cerebral palsy (CP), subject-specific musculoskeletal models are needed. Non-rigid deformation (NRD) deforms generic onto personalized bone geometry and thereby transforms the muscle points. The goal of this study was to determine to what extent the models and simulation outcomes in CP patients differ when including subject-specific detail using NRD or Magnetic Resonance Imaging (MRI)-based models. The NRD models slightly overestimated hip contact forces compared to MRI models and differences in muscle point positions and moment arm lengths (MALs) remained, although differences were smaller than for the generic model.     

Transducer and base compliance alter the in situ 6 dof force measured from muscle during an isometric contraction in a multi-joint limb

Although musculoskeletal models are commonly used, validating the muscle actions predicted by such models is often difficult. In situ isometric measurements are a possible solution. The base of the skeleton is immobilized and the endpoint of the limb is rigidly attached to a 6-axis force transducer. Individual muscles are stimulated and the resulting forces and moments recorded. Such analyses generally assume idealized conditions. In this study we have developed an analysis taking into account the compliances due to imperfect fixation of the skeleton, imperfect attachment of the force transducer, and extra degrees of freedom (dof) in the joints that sometimes become necessary in fixed end contractions. We use simulations of the rat hindlimb to illustrate the consequences of such compliances. We show that when the limb is overconstrained, i.e., when there are fewer dof within the limb than are restrained by the skeletal fixation, the compliances of the skeletal fixation and of the transducer attachment can significantly affect measured forces and moments. When the limb dofs and restrained dofs are matched, however, the measured forces and moments are independent of these compliances. We also show that this framework can be used to model limb dofs, so that rather than simply omitting dofs in which a limb does not move (e.g., abduction at the knee), the limited motion of the limb in these dofs can be more realistically modeled as a very low compliance. Finally, we discuss the practical implications of these results to experimental measurements of muscle actions.     

A fibre reorientation model for orthotropic multiplicative growth

The main goal of this contribution consists in the development of a remodelling framework for orthotropic continua whereby the underlying symmetry group is incorporated via two fibre families. Special emphasis is placed on the modelling of biological tissues at finite deformations. Besides the incorporation of a referential mass source, anisotropic growth is addressed by means of a multiplicative decomposition of the overall deformation gradient into an elastic and a growth distortion. Projected quantities of a configurational growth stress tensor are advocated as driving forces for time-dependent saturation–type evolution of the principal values of the growth distortion. Moreover, the reorientation of both fibre families, which directly affects the strain energy as well as the growth distortion itself, is guided by analyzing critical energy points. In particular, a time-dependent formulation is developed which aligns the fibre directions according to the principal stretch directions. Finally, the proposed framework is embedded into a finite element context so that representative numerical examples, examining growth and resorption in volume and density together with fibre reorientation, close this study.     

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

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

Load-sharing in the lumbosacral spine in neutral standing & flexed postures – A combined finite element and inverse static study

Understanding load-sharing in the spine during in-vivo conditions is critical for better spinal implant design and testing. Previous studies of load-sharing that considered actual spinal geometry applied compressive follower load, with or without moment, to simulate muscle forces. Other studies used musculoskeletal models, which include muscle forces, but model the discs by simple beams or spherical joints and ignore the articular facet joints.This study investigated load-sharing in neutral standing and flexed postures using a detailed Finite Element (FE) model of the ligamentous lumbosacral spine, where muscle forces, gravity loads and intra-abdominal pressure, as predicted by a musculoskeletal model of the upper body, are input into the FE model. Flexion was simulated by applying vertebral rotations following spine rhythm measured in a previous in-vivo study, to the musculoskeletal model. The FE model predicted intradiscal pressure (IDP), strains in the annular fibers, contact forces in the facet joints, and forces in the ligaments. The disc forces and moments were determined using equilibrium equations, which considered the applied loads, including muscle forces and IDP, as well as forces in the ligaments and facet joints predicted by the FE model. Load-sharing was calculated as the portion of the total spinal load carried along the spine by each individual spinal structure. The results revealed that spinal loads which increased substantially from the upright to the flexed posture were mainly supported by the discs in the upright posture, whereas the ligaments’ contribution in resisting shear, compression, and moment was more significant in the flexed posture.