A new PMHS model for lumbar spine injuries during vertical acceleration

Ejection from military aircraft exerts substantial loads on the lumbar spine. Fractures remain common, although the overall survivability of the event has considerably increased over recent decades. The present study was performed to develop and validate a biomechanically accurate experimental model for the high vertical acceleration loading to the lumbar spine that occurs during the catapult phase of aircraft ejection. The model consisted of a vertical drop tower with two horizontal platforms attached to a monorail using low friction linear bearings. A total of four human cadaveric spine specimens (T12-L5) were tested. Each lumbar column was attached to the lower platform through a load cell. Weights were added to the upper platform to match the thorax, head-neck, and upper extremity mass of a 50th percentile male. Both platforms were raised to the drop height and released in unison. Deceleration characteristics of the lower platform were modulated by foam at the bottom of the drop tower. The upper platform applied compressive inertial loads to the top of the specimen during deceleration. All specimens demonstrated complex bending during ejection simulations, with the pattern dependent upon the anterior-posterior location of load application. The model demonstrated adequate inter-specimen kinematic repeatability on a spinal level-by-level basis under different subfailure loading scenarios. One specimen was then exposed to additional tests of increasing acceleration to induce identifiable injury and validate the model as an injury-producing system. Multiple noncontiguous vertebral fractures were obtained at an acceleration of 21 g with 488 g/s rate of onset. This clinically relevant trauma consisted of burst fracture at L1 and wedge fracture at L4. Compression of the vertebral body approached 60% during the failure test, with -6,106 N axial force and 168 Nm flexion moment. Future applications of this model include developing a better understanding of the vertebral injury mechanism during pilot ejection and developing tolerance limits for injuries sustained under a variety of different vertical acceleration scenarios.     

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.     

Twente spine model: A complete and coherent dataset for musculo-skeletal modeling of the lumbar region of the human spine

Musculo-skeletal modeling can greatly help in understanding normal and pathological functioning of the spine. For such models to produce reliable muscle and joint force estimations, an adequate set of musculo-skeletal data is necessary. In this study, we present a complete and coherent dataset for the lumbar spine, based on medical images and dissection measurements from one embalmed human cadaver. We divided muscles into muscle-tendon elements, digitized their attachments at the bones and measured morphological parameters. In total, we measured 11 muscles from one body side, using 96 elements. For every muscle element, we measured three-dimensional coordinates of its attachments, fiber length, tendon length, sarcomere length, optimal fiber length, pennation angle, mass, and physiological cross-sectional area together with the geometry of the lumbar spine. Results were consistent with other anatomical studies and included new data for the serratus posterior inferior muscle. The dataset presented in this paper enables a complete and coherent musculo-skeletal model for the lumbar spine and will improve the current state-of-the art in predicting spinal loading.     

Validation of an OpenSim full-body model with detailed lumbar spine for estimating lower lumbar spine loads during symmetric and asymmetric lifting tasks

There is currently no validated full-body lifting model publicly available on the OpenSim modelling platform to estimate spinal loads during lifting. In this study, the existing full-body-lumbar-spine model was adapted and validated for lifting motions to produce the lifting full-body model. Back muscle activations predicted by the model closely matched the measured erector spinae activation patterns. Model estimates of intradiscal pressures and in vivo measurements were strongly correlated. The same spine loading trends were observed for model estimates and reported vertebral body implant measurements. These results demonstrate the suitability of this model to evaluate changes in lumbar loading during lifting.     

Intra-abdominal pressure mechanism for stabilizing the lumbar spine

Currently, intra-abdominal pressure (IAP) is thought to provide stability to the lumbar spine but the exact principles have yet to be specified. A simplified physical model was constructed and theoretical calculations performed to illustrate a possible intra-abdominal pressure mechanism for stabilizing the spine. The model consisted of an inverted pendulum with linear springs representing abdominal and erector spinae muscle groups. The IAP force was simulated with a pneumatic piston activated with compressed air. The critical load of the model was calculated theoretically based on the minimum potential energy principle and obtained experimentally by increasing weight on the model until the point of buckling. Two distinct mechanisms were simulated separately and in combination. One was antagonistic flexor extensor muscle coactivation and the second was abdominal muscle activation along with generation of IAP. Both mechanisms were effective in stabilizing the model of a lumbar spine. The critical load and therefore the stability of the spine model increased with either increased antagonistic muscle coactivation forces or increased IAP along with increased abdominal spring force. Both mechanisms were also effective in providing mechanical stability to the spine model when activated simultaneously. Theoretical calculation of the critical load agreed very well with experimental results (95.5% average error). The IAP mechanism for stabilizing the lumbar spine appears preferable in tasks that demand trunk extensor moment such as lifting or jumping. This mechanism can increase spine stability without the additional coactivation of erector spinae muscles.     

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.     

A methodological approach for the biomechanical cause analysis of golf-related lumbar spine injuries

A new methodological approach employing mechanical work (MW) determination and relative portion of its elemental analysis was applied to investigate the biomechanical causes of golf-related lumbar spine injuries. Kinematic and kinetic parameters at the lumbar and lower limb joints were measured during downswing in 18 golfers. The MW at the lumbar joint (LJ) was smaller than at the right hip but larger than the MWs at other joints. The contribution of joint angular velocity (JAV) to MW was much greater than that of net muscle moment (NMM) at the LJ, whereas the contribution of NMM to MW was greater rather than or similar to that of JAV at other joints. Thus, the contribution of JAV to MW is likely more critical in terms of the probability of golf-related injury than that of NMM. The MW-based golf-related injury index (MWGII), proposed as the ratio of the contribution of JAV to MW to that of NMM, at the LJ (1.55) was significantly greater than those at other joints ( < 1.05). This generally corresponds to the most frequent occurrence of golf-related injuries around the lumbar spine. Therefore, both MW and MWGII should be considered when investigating the biomechanical causes of lumbar spine injuries.     

A technique for quantifying the bending moment acting on the lumbar spine in vivo

The bending properties of cadaveric lumbar spines were measured and used to convert in vivo measurements of lumbar flexion into bending moments ('stresses'). Forty-two lumbar motion segments were subjected to complex physiological loading and graphs were obtained of bending moment vs flexion angle. Variability was reduced by expressing both variables as a percentage of their values at the elastic limit. Data were averaged for each lumbar level, and a composite bending curve was compiled for the lumbar spine, L1-S1. A linear relationship was established between lumbar flexion measured in vitro and in vivo. This enabled values of 'per cent lumbar flexion' measured in vivo to be converted into 'per cent maximum bending moment' with a maximum likely error of about +/- 8%, which is equivalent to about +/- 5 Nm at L5-S1 for an average person. The technique was applied to 28 subjects, using dynamic measurements of lumbar flexion obtained with the '3-Space Isotrack' system. The bending moment at L5-S1 was 12 Nm on average when picking a pen up off the floor. Highly significant increases in bending moment were observed when heavier and bulkier objects were lifted.     

An EMG-assisted model of loads on the lumbar spine during asymmetric trunk extensions

An EMG-assisted, low-back, lifting model is presented which simulates spinal loading as a function of dynamic, asymmetric, lifting exertions. The purpose of this study has been to develop a model which overcomes the limitations of previous models including static or isokinetic mechanics, inaccurate predictions of muscle coactivity, static interpretation of myoelectric activity, and physiologically unrealistic or variable muscle force per unit area. The present model predicts individual muscle forces from processed EMG data, normalized as a function of trunk angle and asymmetry, and modified to account for muscle length and velocity artifacts. The normalized EMGs are combined with muscle cross-sectional area and intrinsic strength capacity as determined on a per subject basis, to represent tensile force amplitudes. Dynamic internal and external force vectors are employed to predict trunk moments, spinal compression, lateral and anterior shear forces. Data from 20 subjects performing a total of 2160 exertions showed good agreement between predicted and measured values under all trunk angle, asymmetry, velocity, and acceleration conditions. The design represents a significant step toward accurate, fully dynamic modeling of the low-back in multiple dimensions. The benefits of such a model are the insights provided into the effects of motion induced, muscle co-activity on spinal loading in multiple dimensions.     

Quantifying the centre of rotation pattern in a multi-body model of the lumbar spine

Understanding the kinematics of the spine provides paramount knowledge for many aspects of the clinical analysis of back pain. More specifically, visualisation of the instantaneous centre of rotation (ICR) enables clinicians to quantify joint laxity in the segments, avoiding a dependence on more inconclusive measurements based on the range of motion and excessive translations, which vary in every individual. Alternatively, it provides motion preserving designers with an insight into where a physiological ICR of a motion preserving prosthesis can be situated in order to restore proper load distribution across the passive and active elements of the lumbar region. Prior to the use of an unconstrained dynamic musculoskeletal model system, based on multi-body models capable of transient analysis, to estimate segmental loads, the model must be kinematically evaluated for all possible sensitivity due to ligament properties and the initial locus of intervertebral disc (IVD). A previously calibrated osseoligamentous model of lumbar spine was used to evaluate the changes in ICR under variation of the ligament stiffness and initial locus of IVD, when subjected to pure moments from 0 to 15 Nm. The ICR was quantified based on the closed solution of unit quaternion that improves accuracy and prevents coordinate singularities, which is often observed in Euler-based methods and least squares principles. The calculation of the ICR during flexion/extension revealed complexity and intrinsic nonlinearity between flexion and extension. This study revealed that, to accommodate a good agreement between in vitro data and the multi-body model predictions, in flexion more laxity is required than in extension. The results showed that the ICR location is concentrated in the posterior region of the disc, in agreement with previous experimental studies. However, the current multi-body model demonstrates a sensitivity to the initial definition of the ICR, which should be recognised as a limitation of the method. Nevertheless, the current simulations suggest that, due to the constantly evolving path of the ICR across the IVD during flexion–extension, a movable ICR is a necessary condition in multi-body modelling of the spine, in the context of whole body simulation, to accurately capture segmental kinematics and kinetics.     

Three-dimensional static modeling of the lumbar spine

This paper presents three-dimensional static modeling of the human lumbar spine to be used in the formation of anatomically-correct movement patterns for a fully cable-actuated robotic lumbar spine which can mimic in vivo human lumbar spine movements to provide better hands-on training for medical students. The mathematical model incorporates five lumbar vertebrae between the first lumbar vertebra and the sacrum, with dimensions of an average adult human spine. The vertebrae are connected to each other by elastic elements, torsional springs and a spherical joint located at the inferoposterior corner in the mid-sagittal plane of the vertebral body. Elastic elements represent the ligaments that surround the facet joints and the torsional springs represent the collective effect of intervertebral disc which plays a major role in balancing torsional load during upper body motion and the remaining ligaments that support the spinal column. The elastic elements and torsional springs are considered to be nonlinear. The nonlinear stiffness constants for six motion types were solved using a multiobjective optimization technique. The quantitative comparison between the angles of rotations predicted by the proposed model and in the experimental data confirmed that the model yields angles of rotation close to the experimental data. The main contribution is that the new model can be used for all motions while the experimental data was only obtained at discrete measurement points.     

Effect of lumbar fasciae on the stability of the lower lumbar spine

The biomechanical effect of tensioning the lumbar fasciae (LF) on the stability of the spine during sagittal plane motion was analysed using a validated finite element model of the normal lumbosacral spine (L4-S1). To apply the tension in the LF along the direction of the fibres, a local coordinate was allocated using dummy rigid beam elements that originated from the spinous process. Up to 10 Nm of flexion and 7.5 Nm of extension moment was applied with and without 20 N of lateral tension in the LF. A follower load of 400 N was additionally applied along the curvature of the spine. To identify how the magnitude of LF tension related to the stability of the spine, the tensioning on the fasciae was increased up to 40 N with an interval of 10 N under 7.5 Nm of flexion/extension moment. A fascial tension of 20 N produced a 59% decrease in angular motion at 2.5 Nm of flexion moment while there was a 12.3% decrease at 10 Nm in the L5-S1 segment. Its decrement was 53 and 9.6% at 2.5 Nm and 10 Nm, respectively, in the L4-L5 segment. Anterior translation was reduced by 12.1 and 39.0% at the L4-L5 and L5-S1 segments under 10 Nm of flexion moment, respectively. The flexion stiffness shows an almost linear increment with the increase in fascial tension. The results of this study showed that the effect of the LF on the stability of the spine is significant.     

A finite element model technique to determine the mechanical response of a lumbar spine segment under complex loads

This study presents a CT-based finite element model of the lumbar spine taking into account all function-related boundary conditions, such as anisotropy of mechanical properties, ligaments, contact elements, mesh size, etc. Through advanced mesh generation and employment of compound elements, the developed model is capable of assessing the mechanical response of the examined spine segment for complex loading conditions, thus providing valuable insight on stress development within the model and allowing the prediction of critical loading scenarios. The model was validated through a comparison of the calculated force-induced inclination/deformation and a correlation of these data to experimental values. The mechanical response of the examined functional spine segment was evaluated, and the effect of the loading scenario determined for both vertebral bodies as well as the connecting intervertebral disc.     

An upper bound computational model for investigation of fusion effects on adjacent segment biomechanics of the lumbar spine

Abstract

Prediction of the biomechanical effects of fusion surgery on adjacent segments is a challenge in computational biomechanics of the spine. In this study, a two-segment L3-L4-L5 computational model was developed to simulate the effects of spinal fusion on adjacent segment biomechanical responses under a follower load condition. The interaction between the degenerative segment (L4-5) and the adjacent segment (L3-4) was simulated using an equivalent follower spring. The spring stiffness was calibrated using a rigid fusion of a completely degenerated disc model at the L4-5 level, resulting in an upper bound response at the adjacent (L3-4) segment. The obtained upper bound equivalent follower spring was used to simulate the upper bound biomechanical responses of fusion of the disc with different degeneration grades. It was predicted that as the disc degeneration grade at the degenerative segment decreased, the effect on the adjacent segment responses decreased accordingly after fusion. The data indicated that the upper bound computational model can be a useful computational tool for evaluation of the interaction between segments and for investigation of the biomechanical mechanisms of adjacent segment degeneration after fusion.     

A myoelectrically based dynamic three-dimensional model to predict loads on lumbar spine tissues during lateral bending.

This work describes a dynamic model of the low back that incorporates extensive anatomical detail of the musculo-ligamentous-skeletal system to predict the load time histories of individual tissues. The dynamic reaction moment about L4/L5 was determined during lateral bending from a linked-segment model. This reaction moment was partitioned into restorative components provided by the disc, ligament strain, and active-muscle contraction using a second model of the spine that incorporated a detailed representation of the anatomy. Muscle contraction forces were estimated using both information from surface electromyograms, collected from 12 sites, and consideration of the modulating effects of muscle length, cross-sectional area and passive elasticity. This modelling technique is sensitive to the different ways in which individuals recruit their musculature to satisfy moment constraints. Time histories of muscle forces are provided. High muscle loads are consistent with the common clinical observation of muscle strain often produced by load handling. Furthermore, the coactivation measured in muscles on both sides of the trunk suggests that muscles are recruited to satisfy the lateral bending reaction torque in addition to performing other mechanical roles such as spine stabilization. If an estimate of the intervertebral joint compression is desired for assessment of lateral bends in industry, then a single equivalent lateral muscle with a moment arm of approximately 3.0-4.0 cm would conservatively capture the effects of muscle co-contraction quantified in this study.     

Parametric convergence sensitivity and validation of a finite element model of the human lumbar spine

The primary objective of this study was to generate a finite element model of the human lumbar spine (L1–L5), verify mesh convergence for each tissue constituent and perform an extensive validation using both kinematic/kinetic and stress/strain data. Mesh refinement was accomplished via convergence of strain energy density (SED) predictions for each spinal tissue. The converged model was validated based on range of motion, intradiscal pressure, facet force transmission, anterolateral cortical bone strain and anterior longitudinal ligament deformation predictions. Changes in mesh resolution had the biggest impact on SED predictions under axial rotation loading. Nonlinearity of the moment-rotation curves was accurately simulated and the model predictions on the aforementioned parameters were in good agreement with experimental data. The validated and converged model will be utilised to study the effects of degeneration on the lumbar spine biomechanics, as well as to investigate the mechanical underpinning of the contemporary treatment strategies.