Distributed Body Weight Over the Whole Spine for Improved Inference in Spine Modelling

Numerous models used to investigate the causes of low back pain are based upon the concept of the lever model that considers equilibrium of forces and moments about a single intervertebral joint. Consideration of forces and moments at each intervertebral joint is essential if a more realistic idea of the loading on the spine is to be obtained. This will also allow the role of the curvature of the spine to be investigated. A distributed loading pattern for forces due to body weight for the whole spine has not been investigated before. In this paper, a distributed loading pattern for the whole spine for various postures is investigated and the potential impact on the calculations is discussed.     

Diurnal variations in intervertebral disc height affect spine flexibility,intradiscal pressure and contact compressive forces in the facet joints

Diurnal changes of intervertebral disc height are caused by high compressive loading during the day, which expulses fluid from the disc, and by osmotic pressure, which imbibes fluid into the disc at low loading. The aim of the present study was to determine the magnitude of diurnal changes in spine flexibility, intradiscal pressures and contact forces in the facet joints. A validated osseoligamentous finite element model of the lumbar spine was used to determine these quantities for morning and evening situations. Disc height varied by 10% for these two situations. Spine flexibility and facet joint forces were markedly higher in the evening than in the morning. Intradiscal pressures were higher in the morning than in the evening. The different spine flexibilities in the morning and evening should be taken into account during kinematical measurements. Predicted facet joint forces may be used for the designing and pre-clinical testing of artificial facet joint replacements.     

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.     

Effect of a pedicle-screw-based motion preservation system on lumbar spine biomechanics: A probabilistic finite element study with subsequent sensitivity analysis

Pedicle-screw-based motion preservation systems are often used to support a slightly degenerated disc. Such implants are intended to reduce intradiscal pressure and facet joints forces, while having a minimal effect on the motion patterns.In a probabilistic finite element study with subsequent sensitivity analysis, the effects of 10 input parameters, such as elastic modulus and diameter of the elastic rod, distraction of the segment, level of bridged segments, etc. on the output parameters intervertebral rotations, intradiscal pressures, and facet joint forces were determined. A validated finite element model of the lumbar spine was employed. Probabilistic studies were performed for seven loading cases: upright standing, flexion, extension, left and right lateral bending and left and right axial rotation.The simulations show that intervertebral rotation angles, intradiscal pressures and facet joint forces are in most cases reduced by a motion preservation system. The influence on intradiscal pressure is small, except in extension. For many input parameter combinations, the values for intervertebral rotations and facet joint forces are very low, which indicates that the implant is too stiff in these cases. The output parameters are affected most by the following input parameters: loading case, elastic modulus and diameter of the elastic rod, distraction of the segment, and angular rigidity of the connection between screws and rod.The designated functions of a motion preservation system can best be achieved when the longitudinal rod has a low stiffness, and when the connection between rod and pedicle screws is rigid.     

Constraining spine stability levels in an optimization model leads to the prediction of trunk muscle cocontraction and improved spine compression force estimates

A major limitation of optimization models of the spine has been the inability to accurately predict trunk muscle co-activity. Antagonist muscle activity is thought to be necessary to maintain adequate levels of spine stability but, in turn, creates increased loading on the spine. It is thus hypothesized that the CNS attempts to optimize the relationship between spine loading and spine stability in determining muscular activation patterns. This study presents an optimization model of the spine in which stability was constrained to target levels predicted from regression equations of independent loading variables. Objective functions were set to either minimize the sum of the cubed muscle forces or minimize the sum of the squared intervertebral forces at the L4-L5 disc level. Results demonstrate that the inclusion of stability constraints in optimization simulations produced realistic predictions of antagonist muscle activity and predictions of spine compression levels that agree more closely with EMG-based estimates, compared to simulations in which stability was unconstrained. It was concluded that spinal stability is a vital consideration for the CNS when dictating trunk muscle recruitment patterns.     

Function of the spine

In spite of the considerable effort which has been invested in attempts to understand the mechanism of human spines, substantial controversy remains, particularly in connection with assumptions which have to be made by those engaged in biological modelling.The hypothesis presented here is that the living joint has stress sensors driving a feedback mechanism, an arrangement which could react to imposed loads by modifying muscular action in such a way as to minimize stress at the joints and therefore the risk of injury. A theory of this kind gives an image of the spine not in terms of a spatial picture, as would a CAT scan, but in terms of stresses, forces and moments acting at the intervertebral joints. Calculations show that the erectores spinae alone cannot support more than about 50 kg; there must be some other mechanism to explain man's ability substantially to exceed that load.It is suggested that the interaction between the erectores spinae and the abdominals are of fundamental importance in the function of the spine; how they are co-ordinated during the lifting of weights is examined in detail. The theory resulting from this hypothesis is used to relate spinal injury and an injured subject's posture and behaviour.A mathematical formulation permits an objective evaluation of the spine, and a procedure for determining an automatic diagnosis of lumbar spine disabilities is proposed.     

In vitro axial preload application during spine flexibility testing: towards reduced apparatus-related artefacts

Presently, there is little consensus about how, or even if, axial preload should be incorporated in spine flexibility tests in order to simulate the compressive loads naturally present in vivo. Some preload application methods are suspected of producing unwanted “artefact” forces as the specimen rotates and, in doing so, influencing the resulting kinematics. The objective of this study was to quantitatively compare four distinct types of preload which have roots in contemporary experimental practice. The specific quantities compared were the reaction moments and forces resulting at the intervertebral disc and specimen kinematics. The preload types incorporated increasing amounts of caudal constraint on the preload application vector ranging from an unconstrained dead-load arrangement to an apparatus that allowed the vector to follow rotations of the specimen. Six human cadaveric spine segments were tested (1-L1/L2, 3-L2/L3, 1-L3/L4 and 1-L4/L5). Pure moments were applied to the specimens with each of the four different types of compressive preload. Kinematic response was measured using an opto-electronic motion analysis system. A six-axis load cell was used to measure reaction forces and moments. Artefact reaction moments and shear forces were significantly affected by preload application method and magnitude. Unconstrained preload methods produced high artefact moments and low artefact shear forces while more constrained methods did the opposite. A mechanical trade-off is suggested by our results, whereby unwanted moment can only be prevented at the cost of shear force production. When comparing spine flexibility studies, caution should be exercised to ensure preload was applied in a similar manner for all studies. Unwanted moments or forces induced as a result of preload application method may render the comparison of two seemingly similar studies inappropriate.     

A multibody modelling approach to determine load sharing between passive elements of the lumbar spine

The human spinal segment is an inherently complex structure, a combination of flexible and semi-rigid articulating elements stabilised by seven principal ligaments. An understanding of how mechanical loading is shared among these passive elements of the segment is required to estimate tissue failure stresses. A 3D rigid body model of the complete lumbar spine has been developed to facilitate the prediction of load sharing across the passive elements. In contrast to previous multibody models, this model includes a non-linear, six degrees of freedom intervertebral disc, facet bony articulations and all spinal ligaments. Predictions of segmental kinematics and facet joint forces, in response to pure moment loading (flexion-extension), were compared to published in vitro data. On inclusion of detailed representation of the disc and facets, the multibody model fully captures the non-linear flexibility response of the spinal segment, i.e. coupled motions and a mobile instantaneous centre of rotation. Predicted facet joint forces corresponded well with reported values. For the loading case considered, the model predicted that the ligaments are the main stabilising elements within the physiological motion range; however, the disc resists a greater proportion of the applied load as the spine is fully flexed. In extension, the facets and capsular ligaments provide the principal resistance. Overall patterns of load distribution to the spinal ligaments are in agreement with previous predictions; however, the current model highlights the important role of the intraspinous ligament in flexion and the potentially high risk of failure. Several important refinements to the multibody modelling of the passive elements of the spine have been described, and such an enhanced passive model can be easily integrated into a full musculoskeletal model for the prediction of spinal loading for a variety of daily activities.     

A multibody modelling approach to determine load sharing between passive elements of the lumbar spine

The human spinal segment is an inherently complex structure, a combination of flexible and semi-rigid articulating elements stabilised by seven principal ligaments. An understanding of how mechanical loading is shared among these passive elements of the segment is required to estimate tissue failure stresses. A 3D rigid body model of the complete lumbar spine has been developed to facilitate the prediction of load sharing across the passive elements. In contrast to previous multibody models, this model includes a non-linear, six degrees of freedom intervertebral disc, facet bony articulations and all spinal ligaments. Predictions of segmental kinematics and facet joint forces, in response to pure moment loading (flexion–extension), were compared to published in vitro data. On inclusion of detailed representation of the disc and facets, the multibody model fully captures the non-linear flexibility response of the spinal segment, i.e. coupled motions and a mobile instantaneous centre of rotation. Predicted facet joint forces corresponded well with reported values. For the loading case considered, the model predicted that the ligaments are the main stabilising elements within the physiological motion range; however, the disc resists a greater proportion of the applied load as the spine is fully flexed. In extension, the facets and capsular ligaments provide the principal resistance. Overall patterns of load distribution to the spinal ligaments are in agreement with previous predictions; however, the current model highlights the important role of the intraspinous ligament in flexion and the potentially high risk of failure. Several important refinements to the multibody modelling of the passive elements of the spine have been described, and such an enhanced passive model can be easily integrated into a full musculoskeletal model for the prediction of spinal loading for a variety of daily activities.     

Loading of the lumbar spine during backpack carriage

Backpack carriage is significantly associated with a higher prevalence of low back pain. Elevated compression and shear forces in the lumbar intervertebral discs are known risk factors. A novel method of calculating the loads in the lumbar spine during backpack carriage is presented by combining physical and numerical modelling. The results revealed that to predict realistic lumbar compression forces, subject-specific lumbar curvature data were not necessary for loads up to 40 kg. In contrast, regarding shear forces, using subject-specific lumbar curvature data from upright MRI measurements as input for the rigid body model significantly altered lumbar joint force estimates.     

Optimised loads for the simulation of axial rotation in the lumbar spine

Simplified loading modes (pure moment, compressive force) are usually applied in the in vitro studies to simulate flexion-extension, lateral bending and axial rotation of the spine. The load magnitudes for axial rotation vary strongly in the literature. Therefore, the results of current investigations, e.g. intervertebral rotations, are hardly comparable and may involve unrealistic values. Thus, the question 'which in vitro applicable loading mode is the most realistic' remains open. A validated finite element model of the lumbar spine was employed in two sensitivity studies to estimate the ranges of results due to published load assumptions and to determine the input parameters (e.g. torsional moment), which mostly affect the spinal load and kinematics during axial rotation. In a subsequent optimisation study, the in vitro applicable loading mode was determined, which delivers results that fit best with available in vivo measurements. The calculated results varied widely for loads used in the literature with potential high deviations from in vivo measured values. The intradiscal pressure is mainly affected by the magnitude of the compressive force, while the torsional moment influences mainly the intervertebral rotations and facet joint forces. The best agreement with results measured in vivo were found for a compressive follower force of 720N and a pure moment of 5.5Nm applied to the unconstrained vertebra L1. The results reveal that in many studies the assumed loads do not realistically simulate axial rotation. The in vitro applicable simplified loads cannot perfectly mimic the in vivo situation. However, the optimised values lead to the best agreement with in vivo measured values. Their consequent application would lead to a better comparability of different investigations.     

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.     

Joint contact loading in forefoot and rearfoot strike patterns during running

Research concerning forefoot strike pattern (FFS) versus rearfoot strike pattern (RFS) running has focused on the ground reaction force even though internal joint contact forces are a more direct measure of the loads responsible for injury. The main purpose of this study was to determine the internal loading of the joints for each strike pattern. A secondary purpose was to determine if converted FFS and RFS runners can adequately represent habitual runners with regards to the internal joint loading. Using inverse dynamics to calculate the net joint moments and reaction forces and optimization techniques to estimate muscle forces, we determined the axial compressive loading at the ankle, knee, and hip. Subjects consisted of 15 habitual FFS and 15 habitual RFS competitive runners. Each subject ran at a preferred running velocity with their habitual strike pattern and then converted to the opposite strike pattern. Plantar flexor muscle forces and net ankle joint moments were greater in the FFS running compared to the RFS running during the first half of the stance phase. The average contact forces during this period increased by 41.7% at the ankle and 14.4% at the knee joint during FFS running. Peak ankle joint contact force was 1.5 body weights greater during FFS running (p<0.05). There was no evidence to support a difference between habitual and converted running for joint contact forces. The increased loading at the ankle joint for FFS is an area of concern for individuals considering altering their foot strike pattern.     

Validation of lumbar spine loading from a musculoskeletal model including the lower limbs and lumbar spine

Low back mechanics are important to quantify to study injury, pain and disability. As in vivo forces are difficult to measure directly, modeling approaches are commonly used to estimate these forces. Validation of model estimates is critical to gain confidence in modeling results across populations of interest, such as people with lower-limb amputation. Motion capture, ground reaction force and electromyographic data were collected from ten participants without an amputation (five male/five female) and five participants with a unilateral transtibial amputation (four male/one female) during trunk-pelvis range of motion trials in flexion/extension, lateral bending and axial rotation. A musculoskeletal model with a detailed lumbar spine and the legs including 294 muscles was used to predict L4-L5 loading and muscle activations using static optimization. Model estimates of L4-L5 intervertebral joint loading were compared to measured intradiscal pressures from the literature and muscle activations were compared to electromyographic signals. Model loading estimates were only significantly different from experimental measurements during trunk extension for males without an amputation and for people with an amputation, which may suggest a greater portion of L4-L5 axial load transfer through the facet joints, as facet loads are not captured by intradiscal pressure transducers. Pressure estimates between the model and previous work were not significantly different for flexion, lateral bending or axial rotation. Timing of model-estimated muscle activations compared well with electromyographic activity of the lumbar paraspinals and upper erector spinae. Validated estimates of low back loading can increase the applicability of musculoskeletal models to clinical diagnosis and treatment.     

Role of trunk muscles in generating follower load in the lumbar spine of neutral standing posture

Recently, experimental results have demonstrated that the load carrying capacity of the human spine substantially increases under the follower load condition. Thus, it is essential to prove that a follower load can be generated in vivo by activating the appropriate muscles in order to demonstrate the possibility that the stability of the spinal column could be maintained through a follower load mechanism. The aim of this study was to analyze the coordination of the trunk muscles in order to understand the role of the muscles in generating the follower load. A three-dimensional finite element model of the lumbar spine was developed from T12 to S1 and 117 pairs of trunk muscles (58 pairs of superficial muscles and 59 pairs of deep muscles) were considered. The follower load concept was mathematically represented as an optimization problem. The muscle forces required to generate the follower load were predicted by solving the optimization problem. The corresponding displacements and rotations at all nodes were estimated along with the follower forces, shear forces, and joint moments acting on those nodes. In addition, the muscle forces and the corresponding responses were investigated when the activations of the deep muscles or the superficial muscles were restricted to 75% of the maximum activation, respectively. Significantly larger numbers of deep muscles were involved in the generation of the follower load than the number of superficial muscles, regardless of the restriction on muscle activation. The shear force and the resultant joint moment are more influenced by the change in muscle activation in the superficial muscles. A larger number of deep trunk muscles were activated in order to maintain the spinal posture in the lumbar spine. In addition, the deep muscles have a larger capability to reduce the shear force and the resultant joint moment with respect to the perturbation of the external load or muscle fatigue compared to the superficial muscles.     

The effect of control strategies for an active back-support exoskeleton on spine loading and kinematics during lifting

With mechanical loading as the main risk factor for LBP, exoskeletons (EXO) are designed to reduce the load on the back by taking over part of the moment normally generated by back muscles. The present study investigated the effect of an active exoskeleton, controlled using three different control modes (INCLINATION, EMG & HYBRID), on spinal compression forces during lifting with various techniques.Ten healthy male subjects lifted a 15 kg box, with three lifting techniques (free, squat & stoop), each of which was performed four times, once without EXO and once each with the three different control modes. Using inverse dynamics, we calculated L5/S1 joint moments. Subsequently, we estimated spine forces using an EMG-assisted trunk model.Peak compression forces substantially decreased by 17.8% when wearing the EXO compared to NO EXO. However, this reduction was partly, by about one third, attributable to a reduction of 25% in peak lifting speed when wearing the EXO. While subtle differences in back load patterns were seen between the three control modes, no differences in peak compression forces were found. In part, this may be related to limitations in the torque generating capacity of the EXO. Therefore, with the current limitations of the motors it was impossible to determine which of the control modes was best. Despite these limitations, the EXO still reduced both peak and cumulative compression forces by about 18%.     

Static and dynamic bending responses of the human cervical spine.

The quasi-static and dynamic bending responses of the human mid-lower cervical spine were determined using cadaver intervertebral joints fixed at the base to a six-axis load cell. Flexion bending moment was applied to the superior end of the specimen using an electrohydraulic piston. Each specimen was tested under three cycles of quasi-static load-unload and one high-speed dynamic load. A total of five specimens were included in this study. The maximum intervertebral rotation ranged from 11.0 to 15.4 deg for quasi-static tests and from 22.9 to 34.4 deg for dynamic tests. The resulting peak moments at the center of the intervertebral joint ranged from 3.8 to 6.9 Nm for quasi-static tests and from 14.0 to 31.8 Nm for dynamic tests. The quasi-static stiffness ranged from 0.80 to 1.35 Nm/deg with a mean of 1.03 Nm/deg (+/- 0.11 Nm/deg). The dynamic stiffness ranged from 1.08 to 2.00 Nm/deg with a mean of 1.50 Nm/deg (+/- 0.17 Nm/deg). The differences between the two stiffnesses were statistically significant (p < 0.01). Exponential functions were derived to describe the quasi-static and dynamic moment-rotation responses. These results provide input data for lumped-parameter models and validation data for finite element models to better investigate the biomechanics of the human cervical spine.     

Prediction of biomechanical parameters in the lumbar spine during static sagittal plane lifting.

A combined approach involving optimization and the finite element technique was used to predict biomechanical parameters in the lumbar spine during static lifting in the sagittal plane. Forces in muscle fascicles of the lumbar region were first predicted using an optimization-based force model including the entire lumbar spine. These muscle forces as well as the distributed upper body weight and the lifted load were then applied to a three-dimensional finite element model of the thoracolumbar spine and rib cage to predict deformation, the intradiskal pressure, strains, stresses, and load transfer paths in the spine. The predicted intradiskal pressures in the L3-4 disk at the most deviated from the in vivo measurements by 8.2 percent for the four lifting cases analyzed. The lumbosacral joint flexed, while the other lumbar joints extended for all of the four lifting cases studied (rotation of a joint is the relative rotation between its two vertebral bodies). High stresses were predicted in the posterolateral regions of the endplates and at the junctions of the pedicles and vertebral bodies. High interlaminar shear stresses were found in the posterolateral regions of the lumbar disks. While the facet joints of the upper two lumbar segments did not transmit any load, the facet joints of the lower two lumbar segments experienced significant loads. The ligaments of all lumbar motion segments except the lumbosacral junction provided only marginal moments. The limitations of the current model and possible improvements are discussed.     

Realistic loading conditions for upper body bending

Different modes of load applications are used to simulate flexion and extension of the upper body. It is not clear which loading modes deliver realistic results and allow the comparison of different studies.In a numerical study, a validated finite element model of the lumbar spine, ranging from the vertebra L1 to the disc L5–S1 was employed. Each of six different loading modes was studied for simulating flexion and extension, including pure moments, an eccentric axial force, using a wedged fixture, and applying upper body weight plus follower load plus muscle forces. Intersegmental rotations, intradiscal pressures and facet joint contact forces were calculated. Where possible, results were compared to data measured in vivo.The results of the loading modes studied show a large variance for some values. Outcome measures such as flexion angle and intradiscal pressure differed at a segment by up to 44% and 88%, respectively, related to their maximum values. Intradiscal pressure is mainly determined by the magnitude of the applied compressive force. For flexion maximum contact forces between 0 and 69 N are predicted in each facet joint for different loading modes. For both flexion and extension, applying upper body weight plus follower load plus muscle forces as well as a follower load together with a bending moment delivers results which agreed well with in vivo data from the literature.Choosing an adequate loading mode is important in spine biomechanics when realistic results are required for intersegmental rotations, intradiscal pressure and facet joint contact forces. Only then will results of different studies be comparable.     

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.