A combined numerical and experimental technique for estimation of the forces and moments in the lumbar intervertebral disc

Evaluation of the loads on lumbar intervertebral discs (IVD) is critically important since it is closely related to spine biomechanics, pathology and prosthesis design. Non-invasive estimation of the loads in the discs remains a challenge. In this study, we proposed a new technique to estimate in vivo loads in the IVD using a subject-specific finite element (FE) model of the disc and the kinematics of the disc endplates as input boundary conditions. The technique was validated by comparing the forces and moments in the discs calculated from the FE analyses to the in vitro experiment measurements of three corresponding lumbar discs. The results showed that the forces and moments could be estimated within an average error of 20%. Therefore, this technique can be a promising tool for non-invasive estimation of the loads in the discs and may be extended to be used on living subjects.     

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.     

Trunk muscle activation and associated lumbar spine joint shear forces under different levels of external forward force applied to the trunk.

High anterior intervertebral shear loads could cause low back injuries and therefore the neuromuscular system may actively counteract these forces. This study investigated whether, under constant moment loading relative to L3L4, an increased externally applied forward force on the trunk results in a shift in muscle activation towards the use of muscles with more backward directed lines of action, thereby reducing the increase in total joint shear force. Twelve participants isometrically resisted forward forces, applied at several locations on the trunk, while moments were held constant relative to L3L4. Surface EMG and lumbar curvature were measured, and an EMG-driven muscle model was used to calculate compression and shear forces at all lumbar intervertebral joints. Larger externally applied forward forces resulted in a flattening of the lumbar lordosis and a slightly more backward directed muscle force. Furthermore, the overall muscle activation increased. At the T12L1 to L3L4 joint, resulting joint shear forces remained small (less than 200N) because the average muscle force pulled backward relative to those joints. However, at the L5S1 joint the average muscle force pulled the trunk forward so that the increase in muscle force with increasing externally applied forward force caused a further rise in shear force (by 102.1N, SD=104.0N), resulting in a joint shear force of 1080.1N (SD=150.4N) at 50Nm moment loading. It is concluded that the response of the neuromuscular system to shear force challenges tends to increase rather than reduce the shear loading at the lumbar joint that is subjected to the highest shear forces.     

Trunk strength,muscle activity and spinal loads in maximum isometric flexion and extension exertions: A combined in vivo-computational study

Determination of the trunk maximum voluntary exertion moment capacity and associated internal spinal forces could serve in proper selection of workers for specific occupational task requirements, injury prevention and treatment outcome evaluations. Maximum isometric trunk exertion moments in flexion and extension along with surface EMG of select trunk muscles are measured in 12 asymptomatic subjects. Subsequently and under individualized measured harness-subject forces, kinematics and upper trunk gravity, an iterative kinematics-driven finite element model is used to compute muscle forces and spinal loads in 4 of these subjects. Different co-activity and intra-abdominal pressure levels are simulated. Results indicate significantly larger maximal resistant moments and spinal compression/shear forces in extension exertions than flexion exertions. The agonist trunk muscles reach their maximum force generation (saturation) to greater extent in extension exertions compared to flexion exertions. Local lumbar extensor muscles are highly active in extension exertions and generate most of the internal spinal forces. The maximum exertion attempts produce large spinal compression and shear loads that increase with the antagonist co-activity level but decrease with the intra-abdominal pressure. Intra-abdominal pressure decreases agonist muscle forces in extension exertions but generally increase them in flexion exertions.     

Recovering from Laboratory-Induced slips and trips causes high levels of lumbar muscle activity and spine loading

Slips, trips, and falls are some of the most substantial and prevalent causes of occupational injuries and fatalities, and these events may contribute to low-back problems. We quantified lumbar kinematics (i.e., lumbar angles relative to pelvis) and kinetics during unexpected slip and trip perturbations, and during normal walking, among 12 participants (6F, 6 M). Individual anthropometry, lumbar muscle geometry, and lumbar angles, along with electromyography from 14 lumbar muscles were used as input to a 3D, dynamic, EMG-based model of the lumbar spine. Results indicated that, in comparison with values during normal walking, lumbar range of motion, lumbosacral (L5/S1) loads, and lumbar muscle activations were all significantly higher during the slip and trip events. Maximum L5/S1 compression forces exceeded 2700 N during slip and trip events, compared with ∼ 1100 N during normal walking. Mean values of L5/S1 anteroposterior (930 N), and lateral (800 N) shear forces were also substantially larger than the shear force during the normal walking (230 N). These observed levels of L5/S1 reaction forces, along with high levels of bilateral lumbar muscle activities, suggest the potential for overexertion injuries and tissue damage during unexpected slip and trip events, which could contribute to low back injuries. Outcomes of this study may facilitate the identification and control of specific mechanisms involved with low back disorders consequent to slips or trips.     

Elevation and orientation of external loads influence trunk neuromuscular response and spinal forces despite identical moments at the L5–S1 level

A wide range of loading conditions involving external forces with varying magnitudes, orientations and locations are encountered in daily activities. Here we computed the effect on trunk biomechanics of changes in force location (two levels) and orientation (5 values) in 4 subjects in upright standing while maintaining identical external moment of 15 Nm, 30 N m or 45 Nm at the L5–S1. Driven by measured kinematics and gravity/external loads, the finite element models yielded substantially different trunk neuromuscular response with moderate alterations (up to 24% under 45 Nm moment) in spinal loads as the load orientation varied. Under identical moments, compression and shear forces at the L5–S1 as well as forces in extensor thoracic muscles progressively decreased as orientation of external forces varied from downward gravity (90°) all the way to upward (−25°) orientation. In contrast, forces in local lumbar muscles followed reverse trends. Under larger horizontal forces at a lower elevation, lumbar muscles were much more active whereas extensor thoracic muscle forces were greater under smaller forces at a higher elevation. Despite such differences in activity pattern, the spinal forces remained nearly identical (<6% under 45 Nm moment). The published recorded surface EMG data of extensor muscles trend-wise agreed with computed local muscle forces as horizontal load elevation varied but were overall different from results in both local and global muscles when load orientation altered. Predictions demonstrate the marked effect of external force orientation and elevation on the trunk neuromuscular response and spinal forces and questions attempts to estimate spinal loads based only on consideration of moments at a spinal level.     

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.     

Nonlinear gross response analysis of a lumbar motion segment in combined sagittal loadings

A 3-D nonlinear mathematical model is used to analyze the mechanical response of a lumbar L2-3 motion segment including the posterior elements when subjected to combined sagittal plane loads. The loadings consist of axial compression force, anterior and posterior shear forces, and flexion and extension moments. The facet articulation is modelled as a general moving contact problem and the ligaments as a collection of uniaxial elements. The disk nucleus is considered as an inviscid fluid and the annulus as a composite of collagenous fibers embedded in a matrix of ground substance. The presence of axial compression force reduces the segmental stiffness in flexion whereas a reverse trend is predicted in extension. In the presence of axial compression with and without sagittal shear force, flexion considerably increases the intradiscal pressure while extension reduces it. In other words, under an identical compression force, disk pressure is predicted to be noticeably larger in flexion than in extension. The segmental mechanical response in extension loadings is markedly influenced by the changes in the relative geometry of the articular surfaces at the lower regions. Finally, the deformation of the bony structures plays a significant role in the segmental mechanics under relatively large loads.     

Effect of the intra-abdominal pressure and the center of segmental body mass on the lumbar spine mechanics - a computational parametric study

Determination of physiological loads in human lumbar spine is critical for understanding the mechanisms of lumbar diseases and for designing surgical treatments. Computational models have been used widely to estimate the physiological loads of the spine during simulated functional activities. However, various assumptions on physiological factors such as the intra-abdominal pressure (IAP), centers of mass (COMs) of the upper body and lumbar segments, and vertebral centers of rotation (CORs) have been made in modeling techniques. Systematic knowledge of how these assumptions will affect the predicted spinal biomechanics is important for improving the simulation accuracy. In this paper, we developed a 3D subject-specific numerical model of the lumbosacral spine including T12 and 90 muscles. The effects of the IAP magnitude and COMs locations on the COR of each motion segment and on the joint/muscle forces were investigated using a global convergence optimization procedure when the subject was in a weight bearing standing position. The data indicated that the line connecting the CORs showed a smaller curvature than the lordosis of the lumbar spine in standing posture when the IAP was 0?kPa and the COMs were 10?mm anterior to the geometric center of the T12 vertebra. Increasing the IAP from 0 kPa to 10 kPa shifted the location of CORs toward the posterior direction (from 1.4?±?8.9 mm anterior to intervertebral disc (IVD) centers to 40.5?±?3.1 mm posterior to the IVD centers) and reduced the average joint force (from 0.78?±?0.11 Body weight (BW) to 0.31?±?0.07 BW) and overall muscle force (from 349.3?±?57.7 N to 221.5?±?84.2 N). Anterior movement of the COMs from -30 mm to 70 mm relative to the geometric center of T12 vertebra caused an anterior shift of the CORs (from 25.1?±?8.3 mm posterior to IVD centers to 7.8?±?6.2 mm anterior to IVD centers) and increases of average joint forces (from 0.78?±?0.1 BW to 0.93?±?0.1 BW) and muscle force (from 348.9?±?47.7 N to 452.9?±?58.6 N). Therefore, it is important to consider the IAP and correct COMs in order to accurately simulate human spine biomechanics. The method and results of this study could be useful for designing prevention strategies of spinal injuries and recurrences, and for enhancing rehabilitation efficiency.     

Strains and stresses in sub-dermal tissues of the buttocks are greater in paraplegics than in healthy during sitting

A pressure-related deep tissue injury (DTI) is a severe pressure ulcer, which initiates in muscle tissue overlying a bony prominence (e.g. the ischial tuberosities, IT) and progresses outwards through fat and skin, unnoticed by the paralyzed patient. We recently showed that internal strains and stresses in muscle and fat of individuals at anatomical sites susceptible to DTI can be evaluated by integrating Open-MRI scans with subject-specific finite element (FE) analyzes (Linder-Ganz et al., Journal of Biomechanics, 2007); however, sub-dermal soft tissue strains/stresses from paraplegics are still missing in literature. We hypothesize that the pathoanatomy of the buttocks in paraplegia increases the internal soft tissue loads under the IT, making these patients inherently susceptible to DTI. We hence compared the strain and stress peaks in the gluteus muscle and fat tissues under the IT of six healthy and six paraplegic patients, using the coupled MRI-FE method. Peak principal compression, principal tension, von Mises and shear strains in the gluteus were 1.2-, 3.1-, 1.4- and 1.4-fold higher in paraplegics than in healthy, respectively (p<0.02). Likewise, peak principal compression, principal tension, von Mises and shear stresses in the gluteus were 1.9-, 2.5-, 2.1- and 1.7-fold higher for the paraplegics (p<0.05). Peak gluteal compression and shear stresses decreased by as much as 70% when the paraplegic patients moved from a sitting to a lying posture, indicating on the effectiveness of recommending such patients to lie down after prolonged periods of sitting. This is the first attempt to compare internal soft tissue loads between paraplegic and healthy subjects, using an objective standardized bioengineering method of analysis. The findings support our hypothesis that internal tissue loads are significantly higher in paraplegics, and that postural changes significantly affect these loads. The method of analysis is useful for quantifying the effectiveness of various interventions to alleviate sub-dermal tissue loads at sites susceptible to pressure ulcers and DTI, including cushions, mattresses, recommendations for posture and postural changes, etc.     

Load-displacement properties of lower cervical spine motion segments

The load-displacement behavior of 35 fresh adult cervical spine motion segments was measured in compression, shear, flexion, extension, lateral bending and axial torsion tests. Motion segments were tested both intact and with posterior elements removed. Applied forces ranged to 73.6 N in compression and to 39 N in shear, while applied moments ranged to 2.16 Nm. For each mode of loading, principal and coupled motions were measured and stiffnesses were calculated. The effect of disc degeneration on motion segment stiffnesses and the moments required for motion segment failure were also measured. In compression, the stiffnesses of the cervical motion segments were similar to those of thoracic and lumbar motion segments. In other modes of loading, cervical stiffnesses were considerably smaller than thoracic or lumbar stiffnesses. Removal of the posterior elements decreased cervical motion segment stiffnesses by as much as 50%. Degenerated cervical discs were less stiff in compression and stiffer in shear than less degenerated discs, but in bending or axial torsion, no statistically significant differences were evident. Bending moments causing failure were an order of magnitude lower than those for lumbar segments.     

Predictive equations to estimate spinal loads in symmetric lifting tasks

Response surface methodology is used to establish robust and user-friendly predictive equations that relate responses of a complex detailed trunk finite element biomechanical model to its input variables during sagittal symmetric static lifting activities. Four input variables (thorax flexion angle, lumbar/pelvis ratio, load magnitude, and load position) and four model responses (L4–L5 and L5–S1 disc compression and anterior–posterior shear forces) are considered. Full factorial design of experiments accounting for all combinations of input levels is employed. Quadratic predictive equations for the spinal loads at the L4–S1 disc mid-heights are obtained by regression analysis with adequate goodness-of-fit (R²>98%, p<0.05, and low root-mean-squared-error values compared with the range of predicted spine loads). Results indicate that intradiscal pressure values at the L4–L5 disc estimated based on the predictive equations are in close agreement with available in vivo data measured under similar loadings and postures. Combinations of input (posture and loading) variable levels that yield spine loads beyond the tolerance compression limit of 3400 N are identified using contour plots. Ergonomists and bioengineers, faced with the dilemma of using either complex but more accurate models on one hand or less accurate but simple models on the other hand, have thereby easy-to-use predictive equations that quantifies spinal loads and risk of injury under different occupational tasks of interest.     

A comparison of lumbar spine and muscle loading between male and female workers during box transfers

There is a clear relationship between lumbar spine loading and back musculoskeletal disorders in manual materials handling. The incidence of back disorders is greater in women than men, and for similar work demands females are functioning closer to their physiological limit. It is crucial to study loading on the spine musculoskeletal system with actual handlers, including females, to better understand the risk of back disorders. Extrapolation from biomechanical studies conducted on unexperienced subjects (mainly males) might not be applicable to actual female workers. For male workers, expertise changes the lumbar spine flexion, passive spine resistance, and active/passive muscle forces. However, experienced females select similar postures to those of novices when spine loading is critical. This study proposes that the techniques adopted by male experts, male novices, and females (with considerable experience but not categorized as experts) impact their lumbar spine musculoskeletal systems differently. Spinal loads, muscle forces, and passive resistance (muscle and ligamentous spine) were predicted by a multi-joint EMG-assisted optimization musculoskeletal model of the lumbar spine. Expert males flexed their lumbar spine less (avg. 21.9° vs 30.3–31.7°) and showed decreased passive internal moments (muscle avg. 8.9% vs 15.9–16.0%; spine avg. 4.7% vs 7.1–7.8%) and increased active internal moments (avg. 72.9% vs 62.0–63.9%), thus producing a different impact on their lumbar spine musculoskeletal systems. Experienced females sustained the highest relative spine loads (compression avg. 7.3 N/BW vs 6.2–6.4 N/BW; shear avg. 2.3 N/BW vs 1.7–1.8 N/BW) in addition to passive muscle and ligamentous spine resistance similar to novices. Combined with smaller body size, less strength, and the sequential lifting technique used by females, this could potentially mean greater risk of back injury. Workers should be trained early to limit excessive and repetitive stretching of their lumbar spine passive tissues.     

Gait analysis parameters of healthy human subjects with asymmetric loads

This article focuses on the analysis of gait parameters, ground reaction forces (GRF), and motion signals, for the various asymmetric loads carried by healthy human subjects during walking. Eight asymptomatic human volunteers were enrolled in this study. They were asked to walk, at self-selected pace, with various weights ranging from 0 to 11.33 kg (25 lbs) in 2.26 kg (5 lbs) increments, in one hand on a wooden area equipped with a force platform. Moreover, motion data were recorded from lumbar L1 vertebrae at a frequency of 120 Hz. Three trials of data have been recorded for each subject. In order to quantify the effect of increasing loads on the GRF we define the compression area, restitution area, and coefficient of restitution (COR) for GRF curves. We observe an increase in the compression area with respect to the load and almost constant values for the COR. For motion signals analysis we employ wavelet theory. The signals obtained from the lumbar L1 sensor of the spine vertebrae show a decrease in the wavelet detail energy, for the levels 3, 4, and 5, with respect to increasing loads.     

Intervertebral reaction force prediction using an enhanced assembly of OpenSim models

OpenSim offers a valuable approach to investigating otherwise difficult to assess yet important biomechanical parameters such as joint reaction forces. Although the range of available models in the public repository is continually increasing, there currently exists no OpenSim model for the computation of intervertebral joint reactions during flexion and lifting tasks. The current work combines and improves elements of existing models to develop an enhanced model of the upper body and lumbar spine. Models of the upper body with extremities, neck and head were combined with an improved version of a lumbar spine from the model repository. Translational motion was enabled for each lumbar vertebrae with six controllable degrees of freedom. Motion segment stiffness was implemented at lumbar levels and mass properties were assigned throughout the model. Moreover, body coordinate frames of the spine were modified to allow straightforward variation of sagittal alignment and to simplify interpretation of results. Evaluation of model predictions for level L1–L2, L3–L4 and L4–L5 in various postures of forward flexion and moderate lifting (8 kg) revealed an agreement within 10% to experimental studies and model-based computational analyses. However, in an extended posture or during lifting of heavier loads (20 kg), computed joint reactions differed substantially from reported in vivo measures using instrumented implants. We conclude that agreement between the model and available experimental data was good in view of limitations of both the model and the validation datasets. The presented model is useful in that it permits computation of realistic lumbar spine joint reaction forces during flexion and moderate lifting tasks. The model and corresponding documentation are now available in the online OpenSim repository.     

Mathematical and empirical proof of principle for an on-body personal lift augmentation device (PLAD)

In our laboratory, we have developed a prototype of a personal lift augmentation device (PLAD) that can be worn by workers during manual handling tasks involving lifting or lowering or static holding in symmetric and asymmetric postures. Our concept was to develop a human-speed on-body assistive device that would reduce the required lumbar moment by 20-30% without negative consequences on other joints or lifting kinematics. This paper provides mathematical proof using simplified free body diagrams and two-dimensional moment balance equations. Empirical proof is also provided based on lifting trials with nine male subjects who executed sagittal plane lifts using three lifting styles (stoop, squat, free) and three different loads (5, 15, and 25kg) under two conditions (PLAD, No-PLAD). Nine Fastrak sensors and six in-line strap force sensors were used to estimate the reduction of compressive and shear forces on L4/L5 as well as estimate the forces transferred to the shoulders and knees. Depending on lifting technique, the PLAD applied an added 23-36Nm of torque to assist the back muscles during lifting tasks. The peak pelvic girdle contact forces were estimated and their magnitudes ranged from 221.3+/-11.2N for stoop lifting, 324.3+/-17.2N for freestyle lifts to 468.47+/-23.2N for squat lifting. The PLAD was able to reduce the compression and shear forces about 23-29% and 7.9-8.5%, respectively.     

Effects of an anatomically detailed erector spinae model on L4/L5 disc compression and shear

Biomechanical models utilized for analysis of tasks that load the lumbar spine often predict the resultant moment, disc compression and sometimes shear. Usually the extensor muscular and ligament forces of the lumbar spine are assumed to act 5 cm posterior to a disc centre of rotation. This study has re-examined the generation and pathways of muscular force transmission within the extensor musculature. The effects on L4/L5 disc compression and shear estimates of an anatomically and biomechanically justifiable range of tissue moment arms, lines of force and force generating capacity of muscle, input to a computer model, have been determined. Results indicated that L4/L5 compression estimates could be reduced by up to 35% when the output from a more realistic anatomical model of the erector spinae muscle group was compared with that from the frequently reported and simplified single muscle equivalent with a 5 cm moment arm. The shear force estimates could be altered from more than 500 N (L4 tending to shear anteriorly on L5) to less than 200 N with L4 tending to shear posteriorly on L5. Using the combination of input variables considered by the authors to be most feasible to estimate compression, a single 'equivalent' extensor soft tissue moment arm of 7.5 rather than 5 cm would be needed to equate the compression. This simplification of course, does not accommodate the shear force estimate problem.     

Computation of trunk equilibrium and stability in free flexion-extension movements at different velocities

Velocity of movement has been suggested as a risk factor for low-back disorders. The effect of changes in velocity during unconstrained flexion-extension movements on muscle activations, spinal loads, base reaction forces and system stability was computed. In vivo measurements of kinematics and ground reaction forces were initially carried out on young asymptomatic subjects. The collected kinematics of three subjects representing maximum, mean and minimum lumbar rotations were subsequently used in the kinematics-driven model to compute results during the entire movements at three different velocities. Estimated spinal loads and muscle forces were significantly larger in fastest pace as compared to slower ones indicating the effect of inertial forces. Spinal stability was improved in larger trunk flexion angles and fastest movement. Partial or full flexion relaxation of global extensor muscles occurred only in slower movements. Some local lumbar muscles, especially in subjects with larger lumbar flexion and at slower paces, also demonstrated flexion relaxation. Results confirmed the crucial role of movement velocity on spinal biomechanics. Predictions also demonstrated the important role on response of the magnitude of peak lumbar rotation and its temporal variation.     

Determination of trunk muscle forces for flexion and extension by using a validated finite element model of the lumbar spine and measured in vivo data

Muscle forces stabilize the spine and have a great influence on spinal loads. But little is known about their magnitude. In a former in vitro experiment, a good agreement with intradiscal pressure and fixator loads measured in vivo could be achieved for standing and extension of the lumbar spine. However, for flexion the agreement between in vitro and in vivo measurements was insufficient. In order to improve the determination of trunk muscle forces, a three-dimensional nonlinear finite element model of the lumbar spine with an internal fixation device was created and the same loads were applied as in a previous in vitro experiment. An extensive adaptation process of the model was performed for flexion and extension angles up to 20 degrees and -15 degrees, respectively. With this validated computer model intra-abdominal pressure, preload in the fixators, and a combination of hip- and lumbar flexion angle were varied until a good agreement between analytical and in vivo results was reached for both, intradiscal pressure and bending moments in the fixators. Finally, the fixators were removed and the muscle forces for the intact lumbar spine calculated. A good agreement with the in vivo results could only be achieved at a combination of hip- and lumbar flexion. For the intact spine, forces of 170, 100 and 600 N are predicted in the m. erector spinae for standing, 5 degrees extension and 30 degrees flexion, respectively. The force in the m. rectus abdominus for these body positions is less than 25 N. For more than 10 degrees extension the m. erector spinae is unloaded. The finite element method together with in vivo data allows the estimation of trunk muscle forces for different upper body positions in the sagittal plane. In our patients, flexion of the upper body was most likely a combination of hip- and lumbar spine bending.