Minimizing errors during in vitro testing of multisegmental spine specimens: considerations for component selection and kinematic measurement

Apparatus-induced artifacts may invalidate standard spine testing protocols. Kinematic measurements may be compromised by the configuration of motion capture equipment. This study has determined: (1) the influence of machine design (component friction) on in vitro spinal kinetics; (2) the sensitivity of kinematic measurements to variations in the placement of motion capture markers. A spinal loading simulator has been developed to dynamically apply pure bending moments (three axes) with or without a simultaneous compressive preload. Two linear slider types with different friction coefficients, one with caged ball bearings and one with high-precision roller bearings on rails, were mounted and specimen response compared in sequential tests. Three different optoelectronic marker cluster configurations were mounted on the specimen and motion data was captured simultaneously from all clusters during testing. A polymer tube with a uniform bending stiffness approximately equivalent to a polysegmental lumbar spine specimen was selected to allow reproducible behavior over multiple tests. The selection of sliders for linear degrees of freedom had a marked influence on parasitic shear forces. Higher shear forces were recorded with the caged-bearing design than with the high-precision rollers and consequently a higher moment was required to achieve a given rotation. Kinematic accuracy varied with each marker configuration, but in general higher accuracy was achieved with larger marker spacings and situations where markers moved predominantly parallel to the camera's imaging plane. Relatively common alternatives in the mechanical components used in an apparatus for in vitro spine testing can have a significant influence on the measured kinematic and kinetics. Low-magnitude parasitic shear forces due to friction in sliders induces a linearly increasing moment along the length of the specimen, precluding the ideal of pure moment application. This effect is compounded in polysegmental specimens. Kinematic measurements are highly sensitive to marker design and placement, despite equivalent absolute precision of individual marker measurements, however marker configurations can be designed to minimize errors related to spatial distribution and system bias.     

Analysis of large compression loads on lumbar spine in flexion and in torsion using a novel wrapping element

Axial compression on the spine could reach large values especially in lifting tasks which also involve large rotations. Experimental and numerical investigations on the spinal multi motion segments in presence of physiological compression loads cannot adequately be carried out due to the structural instability and artefact loads. To circumvent these problems, a novel wrapping cable element is used in a nonlinear finite element model of the lumbosacral spine (L1-S1) to investigate the role of moderate to large compression loads on the lumbar stiffness in flexion and axial moments/rotations. The compression loads up to 2,700 N was applied with no instability or artefact loads. The lumbar stiffness substantially increased under compression force, flexion moment, and axial torque when applied alone. The presence of compression preloads significantly stiffened the load-displacement response under flexion and axial moments/rotations. This stiffening effect was much more pronounced under larger preloads and smaller moments/rotations. Compression preloads also increased intradiscal pressure, facet contact forces, and maximum disc fibre strain at different levels. Forces in posterior ligaments were, however, diminished with compression preload. The significant increase in spinal stiffness, hence, should be considered in biomechanical studies for accurate investigation of the load partitioning, system stability, and fixation systems/disc prostheses.     

Structural behavior of human lumbar spinal motion segments

The objectives of this study were to obtain linearized stiffness matrices, and assess the linearity and hysteresis of the motion segments of the human lumbar spine under physiological conditions of axial preload and fluid environment. Also, the stiffness matrices were expressed in the form of an 'equivalent' structure that would give insights into the structural behavior of the spine. Mechanical properties of human cadaveric lumbar L2-3 and L4-5 spinal motion segments were measured in six degrees of freedom by recording forces when each of six principal displacements was applied. Each specimen was tested with axial compressive preloads of 0, 250 and 500 N. The displacements were four slow cycles of +/-0.5mm in anterior-posterior and lateral displacements, +/-0.35 mm axial displacement, +/-1.5 degrees lateral rotation and +/-1 degrees flexion-extension and torsional rotations. There were significant increases with magnitude of preload in the stiffness, hysteresis area (but not loss coefficient) and the linearity of the load-displacement relationship. The mean values of the diagonal and primary off-diagonal stiffness terms for intact motion segments increased significantly relative to values with no preload by an average factor of 1.71 and 2.11 with 250 and 500 N preload, respectively (all eight tests p<0.01). Half of the stiffness terms were greater at L4-5 than L2-3 at higher preloads. The linearized stiffness matrices at each preload magnitude were expressed as an equivalent structure consisting of a truss and a beam with a rigid posterior offset, whose geometrical properties varied with preload. These stiffness properties can be used in structural analyses of the lumbar spine.     

Pure moment testing for spinal biomechanics applications: Fixed versus sliding ring cable-driven test designs

In vitro multi-axial bending testing using pure moment loading conditions has become the standard in evaluating the effects of different types of surgical intervention on spinal kinematics. Simple, cable-driven experimental set-ups have been widely adopted because they require little infrastructure. Traditionally, “fixed ring” cable-driven experimental designs have been used; however, there have been concerns with the validity of this set-up in applying pure moment loading. This study involved directly comparing the loading state induced by a traditional “fixed ring” apparatus versus a novel “sliding ring” approach. Flexion-extension bending was performed on an artificial spine model and a single cadaveric test specimen, and the applied loading conditions to the specimen were measured with an in-line multiaxial load cell. The results showed that the fixed ring system applies flexion-extension moments that are 50–60% less than the intended values. This design also imposes non-trivial anterior–posterior shear forces, and non-uniform loading conditions were induced along the length of the specimen. The results of this study indicate that fixed ring systems have the potential to deviate from a pure moment loading state and that our novel sliding ring modification corrects this error in the original test design. This suggests that the proposed sliding ring design should be used for future in vitro spine biomechanics studies involving a cable-driven pure moment apparatus.     

Search for critical loading condition of the spine--a meta analysis of a nonlinear viscoelastic finite element model

The relative vulnerability of spinal motion segments to different loading combinations remains unknown. The meta-analysis described here using the results of a validated L2-L3 nonlinear viscoelastic finite element model was designed to investigate the critical loading and its effect on the internal mechanics of the human lumbar spine. A Box-Behnken experimental design was used to design the magnitude of seven independent variables associated with loads, rotations and velocity of motion. Subsequently, an optimization method was used to find the primary and secondary variables that influence spine mechanical output related to facet forces, disc pressure, ligament forces, annulus matrix compressive/shear stresses and anulus fibers strain. The mechanical responses with respect to the two most-relevant variables were then regressed linearly using the response surface quadratic model. Axial force and sagittal rotation were identified as the most-relevant variables for mechanical responses. The procedure developed can be used to find the critical loading for finite element models with multi input variables. The derived meta-models can be used to predict the risk associated with various loading parameters and in setting safer load limits.     

A frontal plane model of the lumbar spine subjected to a follower load: implications for the role of muscles

Compression on the lumbar spine is 1000 N for standing and walking and is higher during lifting. Ex vivo experiments show it buckles under a vertical load of 80-100 N. Conversely, the whole lumbar spine can support physiologic compressive loads without large displacements when the load is applied along a follower path that approximates the tangent to the curve of the lumbar spine. This study utilized a two-dimensional beam-column model of the lumbar spine in the frontal plane under gravitational and active muscle loads to address the following question: Can trunk muscle activation cause the path of the internal force resultant to approximate the tangent to the spinal curve and allow the lumbar spine to support compressive loads of physiologic magnitudes? The study identified muscle activation patterns that maintained the lumbar spine model under compressive follower load, resulting in the minimization of internal shear forces and bending moments simultaneously at all lumbar levels. The internal force resultant was compressive, and the lumbar spine model, loaded in compression along the follower load path, supported compressive loads of physiologic magnitudes with minimal change in curvature in the frontal plane. Trunk muscles may coactivate to generate a follower load path and allow the ligamentous lumbar spine to support physiologic compressive loads.     

Biomechanical analysis of manual lifting tasks

This paper presents the findings of a study conducted to determine peak forces generated in the human spine while the individual is engaged in lifting maximum acceptable weight. Calculations of forces and moments, acting on each body segment, were based on film data collected on four individuals for twelve variations of the manual lifting task. The variations were defined by: box-size (three different boxes were used), presence or absence of handles, and symmetry and asymmetry of the lifting task (sagittal and nonsagittal lifting). In general, lower loads were accepted for lift when lifting asymmetrically or when lifting boxes without handles or when lifting bigger boxes. However, peak forces (compressive and shear forces in the spine and ground reaction forces) for these situations were not always lower than those generated when handling either compact boxes or boxes with handles or when lifting boxes symmetrically in the sagittal plane. On the basis of these results, it was concluded that lifting loads asymmetrically or in boxes without handles or in bulky boxes is relatively much more stressful than lifting the same load symmetrically or in boxes with handles or in compact boxes.     

Ex vivo measurement of lumbar intervertebral disc pressure using fibre-Bragg gratings

Methods were developed to measure intervertebral disc pressure using optical fibre-Bragg gratings (FBGs). The FBG sensor was calibrated for hydrostatic pressure in a purpose-built apparatus and the average sensitivity was determined to be -5.7 +/- 0.085 pm/MPa (mean +/- SD). The average coefficient of determination (r(2)) for the calibration data was 0.99, and the average hysteresis of the sensor was 2.13% of full scale. The FBG was used to measure intradiscal pressure response to compressive load in five lumbar functional spine units. The pressure measured by the FBG sensor varied linearly with applied compressive load with coefficients of determination ranging from 0.84 to 0.97. The FBG sensor's sensitivity to compressive load ranged from 0.702 +/- 0.043 kPa/N (mean +/- SD) in a L1-L2 specimen, to 1.07 +/- 0.069 kPa/N in a L4-L5 specimen. These measurements agree with those of previous studies in lumbar spines. Two strain gauge pressure sensors were also used to measure intradiscal pressure response to compressive load. The measured pressure sensitivity to load ranged from 0.251 kPa/N (L4-L5) to 0.850 kPa/N (L2-L3). The average difference in pressure sensitivity to load between Sensors 1 and 2 was 12.9% of the value for Sensor 1, with a range from 1.1% to 20.4%, which suggests that disc pressure was not purely hydrostatic. This may have contributed to the difference between the responses of the FBG and strain gauge sensors.     

A dynamic biomechanical evaluation of lifting maximum acceptable loads

A biomechanical evaluation of the job-related stresses imposed upon a worker is a potential means of reducing the high incidence rates of manual material handling injuries in industry. A biomechanical model consisting of seven rigid links joined at six articulations has been developed for this purpose. Using data from cinematographic analysis of lifting motions the model calculates: (1) body position from articulation angles, (2) angular velocities and accelerations, (3) inertial moments and forces, and (4) reactive moments and forces at each articulation, including the L5/S1 joint. Results indicated effects of the common task variables. Larger load and box sizes increased the rise times and peak values of both vertical ground reaction forces and predicted L5/S1 compressive forces. However, boxes with handles resulted in higher L5/S1 compressive forces than for boxes without handles. Also, in lifting the larger boxes the subjects did not sufficiently compensate with reduced box weights in order to maintain uniform L5/S1 compressive forces. Smoothed and rectified EMG of erector spinae muscles correlated significantly with L5/S1 compressive forces, while predicted and measured vertical ground reaction forces also correlated significantly, indicating the validity of the model as a tool for predicting job physical stresses.     

The use of surface strain data and a neural networks solution method to determine lumbar facet joint loads during in vitro spine testing

A new method for determining facet loads during in vitro spine loading using strain gauges and a neural networks solution method was investigated. A test showed that the new solution method was more robust than and as accurate as a previously presented graphical solution method for computing facet loads using surface strain. The technique was subsequently utilized to assess facet loads at L1-L2 during flexibility testing [7.5Nm pure moments in flexion (FL), extension (EX), right and left axial rotation (AR), and right and left lateral bending (LB)], and stiffness testing (FL-EX with 400N compressive follower load) of six human lumbar spine segments (T12-L2). In contrast to other techniques, such as thin film sensors or pressure-sensitive film, the strain-gauge method leaves the facet joint capsule intact during data collection, presumably allowing more natural load transmission. During flexibility tests, the mean (+/-standard deviation) calculated facet loads (in N) were 46.1+/-41.3 (FL), 51.5+/-39.0 (EX), 70.3+/-43.2 (AR-contralateral side), 31.3+/-33.4 (AR-ipsilateral side), 30.6+/-29.1 (LB-contralateral side), and 32.0+/-44.4 (LB-ipsilateral side). During stiffness tests, the calculated facet loads were 45.5+/-40.4 (upright), 46.6+/-41.9 (full FL), and 75.4+/-39.0 (full EX), corresponding to an equivalent of 11.4%, 11.6%, and 18.8% of the compressive follower load (upright, full FL and EX, respectively). The error associated with this technique, which was below 11N for loads up to 125N, is comparable to that reported with other techniques. The new method shows promise for assessing facet load during in vitro spine testing, an important parameter when evaluating new implant systems and surgical techniques.     

Sensitivity of muscle and intervertebral disc force computations to variations in muscle attachment sites

Abstract

The current paper aims at assessing the sensitivity of muscle and intervertebral disc force computations against potential errors in modeling muscle attachment sites. We perturbed each attachment location in a complete and coherent musculoskeletal model of the human spine and quantified the changes in muscle and disc forces during standing upright, flexion, lateral bending, and axial rotation of the trunk. Although the majority of the muscles caused minor changes (less than 5%) in the disc forces, certain muscle groups, for example, quadratus lumborum, altered the shear and compressive forces as high as 353% and 17%, respectively. Furthermore, percent changes were higher in the shear forces than in the compressive forces. Our analyses identified certain muscles in the rib cage (intercostales interni and intercostales externi) and lumbar spine (quadratus lumborum and longissimus thoracis) as being more influential for computing muscle and disc forces. Furthermore, the disc forces at the L4/L5 joint were the most sensitive against muscle attachment sites, followed by T6/T7 and T12/L1 joints. Presented findings suggest that modeling muscle attachment sites based on solely anatomical illustrations might lead to erroneous evaluation of internal forces and promote using anatomical datasets where these locations were accurately measured. When developing a personalized model of the spine, certain care should also be paid especially for the muscles indicated in this work.     

The dynamic flexion/extension properties of the lumbar spine in vitro using a novel pendulum system

The biomechanical properties of the ligamentous cadaver spine have been previously examined using a variety of experimental testing protocols. Ongoing technical challenges in the biomechanical testing of the spine include the application of physiologic compressive loads and the application of dynamic bending moments while allowing unconstrained three-dimensional motion. The purpose of this study was to report the development of a novel pendulum apparatus that addressed these challenges and to determine the effects of various axial compressive loads on the dynamic biomechanical properties of the lumbar functional spinal unit (FSU). Lumbar FSUs were tested in flexion and extension under five axial compressive loads chosen to represent physiologic loading conditions. After an initial rotation, the FSUs behaved as a dynamic, underdamped vibrating elastic system. Bending stiffness and coefficient of damping increased significantly as the compressive pendulum load increased. The apparatus described herein is a relatively simple approach to determining the dynamic bending properties of the FSU, and potentially disc arthroplasty devices. It is capable of applying physiologic compressive loads at dynamic rates without constraining the kinematics of the joints, crucial requirements for testing FSUs in vitro.     

Estimating the L5S1 flexion/extension moment in symmetrical lifting using a simplified ambulatory measurement system

Mechanical loading of the spine has been shown to be an important risk factor for the development of low-back pain. Inertial motion capture (IMC) systems might allow measuring lumbar moments in realistic working conditions, and thus support evaluation of measures to reduce mechanical loading. As the number of sensors limits applicability, the objective of this study was to investigate the effect of the number of sensors on estimates of L5S1 moments.Hand forces, ground reaction forces (GRF) and full-body kinematics were measured using a gold standard (GS) laboratory setup. In the ambulatory setup, hand forces were estimated based on the force plates measured GRF and body kinematics that were measured using (subsets of) an IMC system. Using top-down inverse dynamics, L5S1 flexion/extension moments were calculated.RMSerrors (Nm) were lowest (16.6) with the full set of 17 sensors and increased to 20.5, 22 and 30.6, for 8, 6 and 4 sensors. Absolute errors in peak moments (Nm) ranged from 17.7 to 16.4, 16.9 and 49.3 Nm, for IMC setup’s with 17, 8, 6 and 4 sensors, respectively. When horizontal GRF were neglected for 6 sensors, RMSerrors and peak moment errors decreased from 22 to 17.3 and from 16.9 to 13 Nm, respectively.In conclusion, while reasonable moment estimates can be obtained with 6 sensors, omitting the forearm sensors led to unacceptable errors. Furthermore, vertical GRF information is sufficient to estimate L5S1 moments in lifting.     

Simulating the response of a standing operator to vibration stress by means of a biomechanical model

Under vibration stress the compressive forces transmitted in the joints of a standing operator are composed of nearly static and oscillating force parts. Because these forces can hardly be measured they were assessed by means of a biomechanical model. In the model, 27 rigid bodies with 103 degrees of freedom represent the segments of the human body. 106 force elements imitate the muscles of the trunk and the legs. At first, the model parameter were varied so that for the simulated sitting posture the model fits the seat-to-head transmissibility given in the literature and in ISO/CD 5982. For the standing posture, the transfer functions between the ground acceleration and the oscillating forces in the ankle, the knee, the hip, and the motion segment L3-L4 were computed. According to the moduli of these functions the forces in the ankles are higher than those in the knees or the hips and they nearly come up to the forces in the lumbar spine. Further the results of the simulation indicate that under equal vibration stress in the standing and the sitting posture the differences between the compressive forces in the lumbar spine are small.     

Pure moment testing for spinal biomechanics applications: fixed versus 3D floating ring cable-driven test designs

Pure moment testing has become a standard protocol for in vitro assessment of the effect of surgical techniques or devices on the bending rigidity of the spine. Of the methods used for pure moment testing, cable-driven set-ups are popular due to their low requirements and simple design. Fixed loading rings are traditionally used in conjunction with these cable-driven systems. However, the accuracy and validity of the loading conditions applied with fixed ring designs have raised some concern, and discrepancies have been found between intended and prescribed loading conditions for flexion-extension. This study extends this prior work to include lateral bending and axial torsion, and compares this fixed ring design with a novel "3D floating ring" design. A complete battery of multi-axial bending tests was conducted with both rings in multiple different configurations using an artificial lumbar spine. Applied moments were monitored and recorded by a multi-axial load cell at the base of the specimen. Results indicate that the fixed ring design deviates as much as 77% from intended moments and induces non-trivial shear forces (up to 18 N) when loaded to a non-destructive maximum of 4.5 Nm. The novel 3D floating ring design largely corrects the inherent errors in the fixed ring design by allowing additional directions of unconstrained motion and producing uniform loading conditions along the length of the specimen. In light of the results, it is suggested that the 3D floating ring set-up be used for future pure moment spine biomechanics applications using a cable-driven apparatus.     

The effect of lateral eccentricity on failure loads,kinematics, and canal occlusions of the cervical spine in axial loading

Current neck injury criteria do not include limits for lateral bending combined with axial compression and this has been observed as a clinically relevant mechanism, particularly for rollover motor vehicle crashes. The primary objectives of this study were to evaluate the effects of lateral eccentricity (the perpendicular distance from the axial force to the centre of the spine) on peak loads, kinematics, and spinal canal occlusions of subaxial cervical spine specimens tested in dynamic axial compression (0.5 m/s). Twelve 3-vertebra human cadaver cervical spine specimens were tested in two groups: low and high eccentricity with initial eccentricities of 1 and 150% of the lateral diameter of the vertebral body. Six-axis loads inferior to the specimen, kinematics of the superior-most vertebra, and spinal canal occlusions were measured. High speed video was collected and acoustic emission (AE) sensors were used to define the time of injury. The effects of eccentricity on peak loads, kinematics, and canal occlusions were evaluated using unpaired Student t-tests. The high eccentricity group had lower peak axial forces (1544±629 vs. 4296±1693 N), inferior displacements (0.2±1.0 vs. 6.6±2.0 mm), and canal occlusions (27±5 vs. 53±15%) and higher peak ipsilateral bending moments (53±17 vs. 3±18 Nm), ipsilateral bending rotations (22±3 vs. 1±2°), and ipsilateral displacements (4.5±1.4 vs. −1.0±1.3 mm, p<0.05 for all comparisons). These results provide new insights to develop prevention, recognition, and treatment strategies for compressive cervical spine injuries with lateral eccentricities.     

Lifting over an obstacle: effects of one-handed lifting and hand support on trunk kinematics and low back loading

Mechanical loading of the low back during lifting is a common cause of low back pain. In this study two-handed lifting is compared to one-handed lifting (with and without supporting the upper body with the free hand) while lifting over an obstacle. A 3-D linked segment model was combined with an EMG-assisted trunk muscle model to quantify kinematics and joint loads at the L5S1 joint. Peak total net moments (i.e., the net moment effect of all muscles and soft tissue spanning the joint) were found to be 10+/-3% lower in unsupported one-handed lifting compared to two-handed lifting, and 30+/-8% lower in supported compared to unsupported one-handed lifting. L5S1 joint forces also showed reductions, but not of the same magnitude (18+/-8% and 15+/-10%, respectively, for compression forces, and 15+/-17% and 11+/-14% respectively, for shear forces). Those reductions of low back load were mainly caused by a reduction of trunk and load moment arms relative to the L5S1 joint during peak loading, and, in the case of hand support, by a support force of about 250 N. Stretching one leg backward did not further reduce low back load estimates. Furthermore, one-handed lifting caused an 6+/-8 degrees increase in lateral flexion, a 9+/-5 degrees increase in twist and a 6+/-6 degrees decrease in flexion. Support with the free hand caused a small further increase in lumbar twisting. It is concluded that one-handed lifting, especially with hand support, reduces L5S1 loading but increases asymmetry in movements and moments about the lumbar spine.     

Investigation of optimal follower load path generated by trunk muscle coordination

It has been reported that the center of rotation of each vertebral body is located posterior to the vertebral body center. Moreover, it has been suggested that an optimized follower load (FL) acts posterior to the vertebral body center. However, the optimal position of the FL with respect to typical biomechanical characteristics regarding spinal stabilization, such as joint compressive force, shear force, joint moment, and muscle stress, has not been studied. A variation in the center of rotation of each vertebra was formulated in a three-dimensional finite element model of the lumbar spine with 117 pairs of trunk muscles. Then, the optimal translation of the FL path connecting the centers of rotations was estimated by solving the optimization problem that was to simultaneously minimize the compressive forces, the shear forces, and the joint moments or to minimize the cubic muscle stresses. An upright neutral standing position and a standing position with 200N in both hands were considered. The FL path moved posterior, regardless of the optimization criteria and loading conditions. The FL path moved 5.0 and 7.8mm posterior in upright standing and 4.1mm and 7.0mm posterior in standing with 200N in hands for each optimization scheme. In addition, it was presented that the optimal FL path may have advantages in comparison to the body center FL path. The present techniques may be important in understanding the spine stabilization function of the trunk muscles.     

Spinal loads and trunk muscles forces during level walking – A combined in vivo and in silico study on six subjects

During level walking, lumbar spine is subjected to cyclic movements and intricate loading of the spinal discs and trunk musculature. This study aimed to estimate the spinal loads (T12–S1) and trunk muscles forces during a complete gait cycle.Six men, 24–33 years walk barefoot at self-selected speed (4–5 km/h). 3D kinematics and ground reaction forces were recorded using a motion capturing system and two force plates, implemented in an inverse dynamic musculoskeletal model to predict the spinal loads and trunk muscles forces. Additionally, the sensitivity of the intra-abdominal pressure and lumbar segment rotational stiffness was investigated.Peak spinal loads and trunk muscle forces were between the gait instances of heel strike and toe off. In L4–L5 segment, sensitivity analysis showed that average peak compressive, antero-posterior and medio-lateral shear forces were 130–179%, 2–15% and 1–6%, with max standard deviation (±STD) of 40%, 6% and 3% of the body weight. Average peak global muscles forces were 24–55% (longissimus thoracis), 11–23% (iliocostalis thoracis), 12–16% (external oblique), 17–25% (internal oblique) and 0–8% (rectus abdominus) of body weight whereas, the average peak local muscles forces were 11–19% (longissimus lumborum), 14–31% (iliocostalis lumborum) and 12–17% (multifidus). Maximum ± STD of the global and local muscles forces were 13% and 8% of the body weight.Large inter-individual differences were found in peak compressive and trunk muscles forces whereas the sensitivity analysis also showed a substantial variation.     

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.