Is it possible to estimate the compressive force in the lumbar spine from intradiscal pressure measurements? A finite element evaluation

Spinal loads during dynamic full flexion and return to standing posture in different age and sex groups: A musculoskeletal model study

During forward flexion, spine motion varies due to age and sex differences. Previous studies showed that lumbar/pelvis range of flexion (RoF) and lumbo-pelvic ratio (L/P) are age/sex dependent. How variation of these parameters affects lumbar loading in a normal population requires further assessment. We aimed to estimate lumbar loads during dynamic flexion-return cycle and the differences in peak loads (compression) and corresponding trunk inclinations due to variation in lumbar/pelvis RoF and L/P. Based on in vivo L/P (0.11-3.44), temporal phases of flexion (early, middle, and later), the lumbar (45-55°) and hip (60-79°) RoF; full flexion-return cycles of six seconds were reconstructed for three age groups (20-35, 36-50 and 50+ yrs.) in both sexes. Six inverse dynamic analyses were performed with a 50th percentile model, and differences in peak loads and corresponding trunk inclinations were calculated. Peak loads at L4-L5 were 179 N higher in younger males versus females, but 228 N and 210 N lower in middle-aged and older males, respectively, compared to females. Females exhibited higher trunk inclinations (6°-20°) than males across all age groups. Age related differences in L4-L5 peak loads and corresponding trunk inclinations were found up to 415 N and 19° in males and 152 N and 13° in females. With aging, peak loads were reduced in males but were found non-monotonic in females, whereas trunk inclinations at peak loads were reduced in both sexes from young to middle/old age groups. In conclusion, lumbar loading and corresponding trunk inclinations varied notably due to age/sex differences. Such data may help distinguishing normal or pathological condition of the lumbar spine.

The role of the nucleus pulposus in intervertebral disc recovery: Towards improved specifications for nucleus replacement devices

Nucleus replacement devices (NRDs) have potential to treat degenerated or herniated intervertebral discs (IVDs). However, IVD height loss is a post-treatment complication. IVD height recovery involves the nucleus pulposus (NP), but the mechanism of this in response to physiological loads is not fully elucidated. This study aimed to characterise the non-linear recovery behaviour of the IVD in intact, post-nuclectomy, and post-NRD treatment states, under physiological loading. 36 bovine IVDs (12 intact, 12 post-nuclectomy, 12 post-treatment) underwent creep-recovery protocols simulating Sitting, Walking or Running, followed by 12 h of recovery. A rheological model decoupled the fluid-independent (elastic, fast) and fluid-dependent (slow) recovery phases. In post-nuclectomy and post-treatment groups, nuclectomy efficiency (ratio of NP removed to remaining NP) was quantified following post-test sectioning. Relative to intact, post-nuclectomy recovery significantly decreased in Sitting (-0.3 ± 0.4 mm, p < 0.05) and Walking (-0.6 ± 0.3 mm, p < 0.001) coupled with significant decreases to the slow response (p < 0.05). Post-nuclectomy, the fast and slow responses negatively correlated with nuclectomy efficiency (p < 0.05). In all protocols, the post-treatment group performed significantly worse in recovery (-0.5 ± 0.3 mm, p < 0.01) and the slow response (p < 0.05). Results suggest the NP mainly facilitates slow-phase recovery, linearly dependent on the amount of NP present. Failure of this NRD to recover is attributed to poor fluid imbibition. Additionally, unconfined NRD performance cannot be extrapolated to the in vitro response. This knowledge informs NRD design criteria to provide high osmotic pressure, and encourages testing standards to incorporate long-term recovery protocols.

Dynamic assessment for low back-support exoskeletons during manual handling tasks

Exoskeletons can protect users' lumbar spine and reduce the risk of low back injury during manual lifting tasks. Although many exoskeletons have been developed, their adoptability is limited by their task- and movement-specific effects on reducing burden. Many studies have evaluated the safety and effectiveness of an exoskeleton using the peak/mean values of biomechanical variables, whereas the performance of the exoskeleton at other time points of the movement has not been investigated in detail. A functional analysis, which presents discrete time-series data as continuous functions, makes it possible to highlight the features of the movement waveform and determine the difference in each variable at each time point. This study investigated an assessment method for exoskeletons based on functional ANOVA, which made it possible to quantify the differences in the biomechanical variables throughout the movement when using an exoskeleton. Additionally, we developed a method based on the interpolation technique to estimate the assistive torque of an exoskeleton. Ten men lifted a 10-kg box under symmetric and asymmetric conditions five times each. Lumbar load was significantly reduced during all phases (flexion, lifting, and laying) under both conditions. Additionally, reductions in kinematic variables were observed, indicating the exoskeleton's impact on motion restrictions. Moreover, the overlap F-ratio curves of the lumbar load and kinematic variables imply that exoskeletons reduce the lumbar load by restricting the kinematic variables. The results suggested that at smaller trunk angles (<25°), an exoskeleton neither significantly reduces the lumbar load nor restricts trunk movement. Our findings will help increasing exoskeleton safety and designing effective products for reducing lumbar injury risks.

Computational lumbar spine models: A literature review

BACKGROUND

Computational spine models of various types have been employed to understand spine function, assess the risk that different activities pose to the spine, and evaluate techniques to prevent injury. The areas in which these models are applied has expanded greatly, potentially beyond the appropriate scope of each, given their capabilities. A comprehensive understanding of the components of these models provides insight into their current capabilities and limitations.

METHODS

The objective of this review was to provide a critical assessment of the different characteristics of model elements employed across the spectrum of lumbar spine modeling and in newer combined methodologies to help better evaluate existing studies and delineate areas for future research and refinement.

FINDINGS

A total of 155 studies met selection criteria and were included in this review. Most current studies use either highly detailed Finite Element models or simpler Musculoskeletal models driven with in vivo data. Many models feature significant geometric or loading simplifications that limit their realism and validity. Frequently, studies only create a single model and thus can't account for the impact of subject variability. The lack of model representation for certain subject cohorts leaves significant gaps in spine knowledge. Combining features from both types of modeling could result in more accurate and predictive models.

INTERPRETATION

Development of integrated models combining elements from different model types in a framework that enables the evaluation of larger populations of subjects could address existing voids and enable more realistic representation of the biomechanics of the lumbar spine.

Biomechanical Properties of Paraspinal Muscles Influence Spinal Loading—A Musculoskeletal Simulation Study

Paraspinal muscles are vital to the functioning of the spine. Changes in muscle physiological cross-sectional area significantly affect spinal loading, but the importance of other muscle biomechanical properties remains unclear. This study explored the changes in spinal loading due to variation in five muscle biomechanical properties: passive stiffness, slack sarcomere length (SSL), in situ sarcomere length, specific tension, and pennation angle. An enhanced version of a musculoskeletal simulation model of the thoracolumbar spine with 210 muscle fascicles was used for this study and its predictions were validated for several tasks and multiple postures. Ranges of physiologically realistic values were selected for all five muscle parameters and their influence on L4-L5 intradiscal pressure (IDP) was investigated in standing and 36° flexion. We observed large changes in IDP due to changes in passive stiffness, SSL, in situ sarcomere length, and specific tension, often with interesting interplays between the parameters. For example, for upright standing, a change in stiffness value from one tenth to 10 times the baseline value increased the IDP only by 91% for the baseline model but by 945% when SSL was 0.4 μm shorter. Shorter SSL values and higher stiffnesses led to the largest increases in IDP. More changes were evident in flexion, as sarcomere lengths were longer in that posture and thus the passive curve is more influential. Our results highlight the importance of the muscle force-length curve and the parameters associated with it and motivate further experimental studies on in vivo measurement of those properties.

Sensitivity of the Cervical Disc Loads, Translations, Intradiscal Pressure, and Muscle Activity Due to Segmental Mass, Disc Stiffness, and Muscle Strength in an Upright Neutral Posture

Musculoskeletal disorders of the cervical spine have increased considerably in recent times. To understand the effects of various biomechanical factors, quantifying the differences in disc loads, motion, and muscle force/activity is necessary. The kinematic, kinetic, or muscle response may vary in a neutral posture due to interindividual differences in segmental mass, cervical disc stiffness, and muscle strength. Therefore, our study aimed to develop an inverse dynamic model of the cervical spine, estimate the differences in disc loads, translations, intradiscal pressure, and muscle force/activity in a neutral posture and compare these results with data available in the literature. A head–neck complex with nine segments (head, C1–T1) was developed with joints having three rotational and three translational degrees of freedom, 517 nonlinear ligament fibers, and 258 muscle fascicles. A sensitivity analysis was performed to calculate the effect of segmental mass (5th to 95th percentile), translational disc stiffness (0.5–1.5), and muscle strength (0.5–1.5) on the cervical disc loads (C2–C3 to C7–T1), disc translations, intradiscal pressure, and muscle force/activity in a neutral posture. In addition, two axial external load conditions (0 and 40 N) were also considered on the head. The estimated intradiscal pressures (0.2–0.56 MPa) at 0 N axial load were comparable to in vivo measurements found in the literature, whereas at 40 N, the values were 0.39–0.93 MPa. With increased segmental mass (5th to 95th), the disc loads, translations, and muscle forces/activities increased to 69% at 0 N and 34% at 40 N axial load. With increased disc stiffness (0.5–1.5), the maximum differences in axial (&lt;1%) and shear loads (4%) were trivial; however, the translations were reduced by 67%, whereas the differences in individual muscle group forces/activities varied largely. With increased muscle strength (0.5–1.5), the muscle activity decreased by 200%. For 40 vs. 0 N, the differences in disc loads, translations, and muscle forces/activities were in the range of 52–129%. Significant differences were estimated in disc loads, translations, and muscle force/activity in the normal population, which could help distinguish between normal and pathological cervical spine conditions.

Comparison of In Vivo Intradiscal Pressure between Sitting and Standing in Human Lumbar Spine: A Systematic Review and Meta-Analysis

Background: Non-specific low back pain (LBP) is highly prevalent today. Disc degeneration could be one of the causes of non-specific LBP, and increased intradiscal pressure (IDP) can potentially induce disc degeneration. The differences in vivo IDP in sitting and standing postures have been studied, but inconsistent results have been reported. The primary objective of this systematic review is to compare the differences in vivo IDP between sitting and standing postures. The secondary objective of this review is to compare effect size estimates between (1) dated and more recent studies and (2) healthy and degenerated intervertebral discs. Methods: An exhaustive search of six electronic databases for studies published before November 2021 was conducted. Articles measuring in vivo IDP in sitting and standing postures were included. Two independent researchers conducted the screening and data extraction. Results: Ten studies that met the inclusion criteria were included in the systematic review, and seven studies with nine independent groups were included in meta-analyses. The sitting posture induces a significantly higher IDP on the lumbar spine (SMD: 0.87; 95% CI = [0.33, 1.41]) than the standing posture. In studies published after 1990 and subjects with degenerated discs, there are no differences in vivo IDP between both postures. Conclusions: Sitting causes higher loads on the lumbar spine than standing in the normal discs, but recent studies do not support this conclusion. Furthermore, the degenerated discs showed no difference in IDP in both postures.

Stress analysis of intervertebral disc during occupational activities

BACKGROUND AND OBJECTIVE

Manual material handling activities cause large compression of the intervertebral disc of the lumbar spine. Intradiscal pressure (IDP) has generally been employed to predict the risk of low back injury. As an alternative to in vivo measurements, either motion analysis or finite element (FE) analysis has been used to estimate IDP. The purpose of this study is to propose a new biomechanical method that integrates FE analysis with motion analysis, in order to estimate the stresses and deformations of the intervertebral disc of the lumbar spine during occupational activities.

METHODS

In the proposed method, motion analysis is performed first by using motion capture data, and the results are employed as input data to FE analysis at specific times of interest during motion. In this method, an in-house interface program is used to scale an initial reference FE model to the subject of experiment, and transformed to the corresponding posture at a specific time during motion. The muscle forces and GRF obtained from motion analysis are applied to FE analysis as boundary and loading conditions. For a total of eighteen occupational activities, the IDP, shear stress, and strain of the L4-L5 segment are estimated.

RESULTS

Under each in vivo activity, the predicted IDP was in overall agreement with the available in vivo data. For lifting activities according to lift origin position, the maximum IDP occurred in the far-knee position immediately after lifting. As the lift origin position moved away from the spine, the stresses and strains in the disc increased.

CONCLUSIONS

This new proposed method is expected to allow the estimation of the stresses and deformations in the intervertebral disc during various occupational activities.

Calibration and validation of a novel hybrid model of the lumbosacral spine in ArtiSynth-The passive structures

In computational biomechanics, two separate types of models have been used predominantly to enhance the understanding of the mechanisms of action of the lumbosacral spine (LSS): Finite element (FE) and musculoskeletal multibody (MB) models. To combine advantages of both models, hybrid FE-MB models are an increasingly used alternative. The aim of this paper is to develop, calibrate, and validate a novel passive hybrid FE-MB open-access simulation model of a ligamentous LSS using ArtiSynth. Based on anatomical data from the Male Visible Human Project, the LSS model is constructed from the L1-S1 rigid vertebrae interconnected with hyperelastic fiber-reinforced FE intervertebral discs, ligaments, and facet joints. A mesh convergence study, sensitivity analyses, and systematic calibration were conducted with the hybrid functional spinal unit (FSU) L4/5. The predicted mechanical responses of the FSU L4/5, the lumbar spine (L1-L5), and the LSS were validated against literature data from in vivo and in vitro measurements and in silico models. Spinal mechanical responses considered when loaded with pure moments and combined loading modes were total and intervertebral range of motions, instantaneous axes and centers of rotation, facet joint contact forces, intradiscal pressures, disc bulges, and stiffnesses. Undesirable correlations with the FE mesh were minimized, the number of crisscrossed collagen fiber rings was reduced to five, and the individual influences of specific anatomical structures were adjusted to in vitro range of motions. Including intervertebral motion couplings for axial rotation and nonlinear stiffening under increasing axial compression, the predicted kinematic and structural mechanics responses were consistent with the comparative data. The results demonstrate that the hybrid simulation model is robust and efficient in reproducing valid mechanical responses to provide a starting point for upcoming optimizations and extensions, such as with active skeletal muscles.

An open-source musculoskeletal model of the lumbar spine and lower limbs: a validation for movements of the lumbar spine

Musculoskeletal models of the lumbar spine have been developed with varying levels of detail for a wide range of clinical applications. Providing consistency is ensured throughout the modelling approach, these models can be combined with other computational models and be used in predictive modelling studies to investigate bone health deterioration and the associated fracture risk. To provide precise physiological loading conditions for such predictive modelling studies, a new full-body musculoskeletal model including a detailed and consistent representation of the lower limbs and the lumbar spine was developed. The model was assessed against in vivo measurements from the literature for a range of spine movements representative of daily living activities. Comparison between model estimations and electromyography recordings was also made for a range of lifting tasks. This new musculoskeletal model will provide a comprehensive physiological mechanical environment for future predictive finite element modelling studies on bone structural adaptation. It is freely available on https://simtk.org/projects/llsm/.

Development of a multiscale model of the human lumbar spine for investigation of tissue loads in people with and without a transtibial amputation during sit-to-stand

Quantification of lumbar spine load transfer is important for understanding low back pain, especially among persons with a lower limb amputation. Computational modeling provides a helpful solution for obtaining estimates of in vivo loads. A multiscale model was constructed by combining musculoskeletal and finite element (FE) models of the lumbar spine to determine tissue loading during daily activities. Three-dimensional kinematic and ground reaction force data were collected from participants with ([Formula: see text]) and without ([Formula: see text]) a unilateral transtibial amputation (TTA) during 5 sit-to-stand trials. We estimated tissue-level load transfer from the multiscale model by controlling the FE model with intervertebral kinematics and muscle forces predicted by the musculoskeletal model. Annulus fibrosis stress, intradiscal pressure (IDP), and facet contact forces were calculated using the FE model. Differences in whole-body kinematics, muscle forces, and tissue-level loads were found between participant groups. Notably, participants with TTA had greater axial rotation toward their intact limb ([Formula: see text]), greater abdominal muscle activity ([Formula: see text]), and greater overall tissue loading throughout sit-to-stand ([Formula: see text]) compared to able-bodied participants. Both normalized (to upright standing) and absolute estimates of L4-L5 IDP were close to in vivo values reported in the literature. The multiscale model can be used to estimate the distribution of loads within different lumbar spine tissue structures and can be adapted for use with different activities, populations, and spinal geometries.

Investigation of lumbar spine biomechanics using global convergence optimization and constant loading path methods

Computational models and inverse dynamic optimization methods are used to predict in-vivo spinal loading. Spinal force is conventionally predicted using the constant loading path method, which is based on the concept that the physiological directions of the spine loads follow the same path of the spinal curve. However, the global convergence optimization method, in which the instantaneous center of rotation of the joint should be also predicted, is necessary for accurate prediction of joint forces of the human body. In this study, we investigate the joint forces, instantaneous centers of rotation, and muscle forces of the human lumbar spine using both global convergence optimization method and constant loading path method during flexion, upright standing, and extension postures. The joint forces predicted using the constant loading path method were 130%, 234%, and 253% greater than those predicted using the global convergence optimization method for the three postures. The instantaneous centers of rotation predicted using the global convergence optimization method were segment level-dependent and moved anteriorly in the flexion and posteriorly in the extension, whereas those predicted using the constant loading path method moved posteriorly in both the flexion and extension. The data indicated that compared to the global convergence optimization method, the constant loading path method introduces additional constraints to the spinal joint model, and thus, it results in greater joint and muscle forces.

Lumbar spinal ligament characteristics extracted from stepwise reduction experiments allow for preciser modeling than literature data

Lumbar ligaments play a key role in stabilizing the spine, particularly assisting muscles at wide-range movements. Hence, valid ligament force–strain data are required to generate physiological model predictions. These data have been obtained by experiments on single ligaments or functional units throughout the literature. However, contrary to detailed spine geometries, gained, for instance, from CT data, ligament characteristics are often inattentively transferred to multi-body system (MBS) or finite element models. In this paper, we use an elaborated MBS model of the lumbar spine to demonstrate how individualized ligament characteristics can be obtained by reversely reenacting stepwise reduction experiments, where the range of motion (ROM) was measured. We additionally validated the extracted characteristics with physiological experiments on intradiscal pressure (IDP). Our results on a total of in each case 160 ROM and 49 IDP simulations indicated superiority of our procedure (seven and eight outliers) toward the incorporation of classical literature data (on average 71 and 31 outliers).

Coupled artificial neural networks to estimate 3D whole-body posture, lumbosacral moments, and spinal loads during load-handling activities

Biomechanical modeling approaches require body posture to evaluate the risk of spine injury during manual material handling. The procedure to measure body posture via motion-analysis techniques as well as the subsequent calculations of lumbosacral moments and spine loads by, respectively, inverse-dynamic and musculoskeletal models are complex and time-consuming. We aim to develop easy-to-use yet accurate artificial neural networks (ANNs) that predict 3D whole-body posture (ANN_posture), segmental orientations (ANN_angle), and lumbosacral moments (ANN_moment) based on our measurements during load-handling activities. Fifteen individuals each performed 135 load-handling activities by reaching (0 kg) or handling (5 and 10 kg) weights located at nine different horizontal and five vertical (0, 30, 60, 90, and 120 cm from the floor) locations. Whole-body posture was measured via a motion capture system and lumbosacral moments were calculated via a 3D top-down eight link-segment inverse-dynamic model. ANN_posture, ANN_angle, and ANN_moment were trained (RMSEs = 6.7 cm, 29.8°, and 16.2 Nm, respectively) and their generalization capability was tested (RMSE = 7.0 cm and R² = 0.97, RMSE = 29.9° and R² = 0.85, and RMSE = 16.5 Nm and R² = 0.97, respectively). These ANNs were subsequently coupled to our previously-developed/validated ANN_load, which predicts spinal loads during 3D load-handling activities. The results showed outputs of the coupled ANNs for L4-L5 intradiscal pressure (IDPs) during a number of activities were in agreement with measured IDPs (RMSE = 0.37 MPa and R² = 0.89). Hence, coupled ANNs were found to be robust tools to evaluate posture, lumbosacral moments, spinal loads, and thus risk of injury during load-handling activities.

Musculoskeletal full-body models including a detailed thoracolumbar spine for children and adolescents aged 6-18 years

Currently available musculoskeletal inverse-dynamics thoracolumbar spine models are entirely based on data from adults and might therefore not be applicable for simulations in children and adolescents. In addition, these models lack lower extremities, which are required for comprehensive evaluations of functional activities or therapeutic exercises. We therefore created OpenSim-based musculoskeletal full-body models including a detailed thoracolumbar spine for children and adolescents aged 6-18 years and validated by comparing model predictions to in vivo data. After combining our recently developed adult thoracolumbar spine model with a lower extremity model, children and adolescent models were created for each year of age by adjusting segmental length and mass distribution, center of mass positions and moments of inertia of the major body segments as well as sagittal pelvis and spine alignment based on literature data. Similarly, muscle strength properties were adjusted based on CT-derived cross-sectional area measurements. Simulations were conducted from in vivo studies reported in the literature involving children and adolescents evaluating maximum trunk muscle strength (MTMS), lumbar disc compressibility (LDC), intradiscal pressure (IDP) and trunk muscle activity (MA). Model predictions correlated highly with in vivo data (MTMS: r ≥ 0.82, p ≤ 0.03; LDC: r = 0.77, p < 0.001; IDP: r ≥ 0.78, p < 0.001; MA: r ≥ 0.90, p < 0.001), indicating suitability for the reasonably accurate prediction of maximal trunk muscle strength, segmental loading and trunk muscle activity in children and adolescents. When aiming at investigating children or adolescents with pathologies such as idiopathic scoliosis, our models can serve as a basis for the creation of deformed spine models and for comparative purposes.

Significance of spine stability criteria on trunk muscle forces following unilateral muscle weakening: A comparison between kinematics-driven and stability-based kinematics-driven musculoskeletal models

Two optimization-driven approaches were employed to develop kinematics-driven (KD) and stability-based kinematics-driven (SKD) musculoskeletal models of an adult thoracolumbar to ascertain the significance of spine stability in holding the upright-standing posture after muscular disuse atrophy. Both models were used to estimate muscle forces of the trunk with intact and unilaterally reduced longissimus thoracis pars thoracic (LGPT) and multifidus lumborum (MFL) muscles strength. A finite element model of the L5-S1 segment of the same kinematics was also developed to compare the joint stresses predicted by the KD and SKD models. Matching well with in vivo data, the SKD model predicted a 15% and 33% reduction in contralateral muscle forces to the 95% debilitated LGPT and MFL muscles, respectively. In contrast, the contralateral muscle force enhancement to the debilitated MFL muscle in the KD model was in contradiction with in vivo data, implying that the KD model is incapable of correctly predicting the muscular disorders. However, the similarity of both models' predictions of intradiscal pressures and intervertebral discs' stresses, which matched well with in vivo data, does indicate the feasibility of the KD model to investigate trunk muscle weakness effects on spinal loads, which could offer additional tools for research in ergonomics. Nonetheless, SKD models can be employed for assessment of contralateral muscle impotence in spinal neuromuscular disorders.

Twente Spine Model: A thorough investigation of the spinal loads in a complete and coherent musculoskeletal model of the human spine

Although in vivospinal loads have been previously measured, existing data are limited to certain lumbar and thoracic levels. A detailed investigation of spinal loads would assist with injury prevention and implant design but is unavailable. In this study, we developed a complete and coherent musculoskeletal model of the entire human spine and studied the intervertebral disc compression forces for physiological movements on three anatomical planes. This model incorporates the individual vertebrae at the cervical, thoracic, and lumbar regions, a flexible ribcage, and complete muscle anatomy. Intradiscal pressures were estimated from predicted compressive forces, and these were generally in close agreement with previously measured data. We found that compressive forces at the trunk discs increased during trunk lateral bending and axial rotation of the trunk. During flexion, compressive forces increased in the thoracolumbar and lumbar regions and slightly decreased at the middle thoracic discs. In extension, the forces generally decreased at the thoracolumbar and lumbar discs whereas they slightly increased at the upper and middle thoracic discs. Furthermore, similar to a previous biomechanical model of the cervical spine, our model predicted increased compression forces in neck flexion, lateral bending, and axial rotation, and decreased forces in neck extension.

Validation of Pre-operative Templating for Total Disc Replacement Surgery

Objectives

This was an analytic retrospective observational study. The aims were (1) to validate patient-specific templating process by comparing postoperative range of motion (ROM) with that predicted by the model, (2) to retrospectively determine the ideal implant size, height, configuration, and location to evaluate if the ROM achieved could have been improved, and (3) to correlate postoperative ROM and clinical outcome.

Background

Previous research revealed that after total disc replacement surgery, 34% of patients with less than 5° of postoperative ROM developed adjacent segment disease. The match between patient anatomy (size, facet orientation, disc height) and implant parameters are likely to affect postoperative ROM and clinical outcomes.

Methods

Seventeen consecutive patients were implanted with 22 ProDisc-L devices between 2008 and 2015. Three-dimensional finite element (FE) models of the implanted segment were constructed from preoperative computed tomography scans and virtually implanted with the ProDisc-L implant. ROM was determined with the endpoints of facet impingement in flexion and implant contact in extension. FE templating was used to determine the optimal implant size and location. ROM was then measured directly from flexion and extension radiographs and compared to predicted ROM. Pre and postoperative Oswestry Disability Index (ODI) data were used to correlate ROM with clinical outcomes.

Results

No significant difference was found between the actual and predicted ROM. The computational templating procedure identified an optimal ROM that was significantly greater than actual ROM. The ROM in our cohort could have been improved by an average of 1.2° or 12% had a different implant size or position been used.

Conclusions

FE analyses accurately predicted ROM in this cohort and can facilitate selection of the optimal implant size and location that we believe will increase the chance of achieving clinical success with the application of this technology.

Variation in human vertebral body strength for vertebral body samples from different locations in segments L1-L5

BACKGROUND

The human spine, in particular the lumbar spine, is subject to significant compressive and bending stresses, which affect the structure of the bone tissue of the vertebrae. The more heterogeneous the structure of the spongy bone tissue, the less resistant the whole vertebral body. It is therefore necessary to establish variations in bone strength parameters within one particular vertebral body.

METHODS

The research material comprised human L1-L5 lumbar vertebrae sampled from 15 donors aged 29-35. A total of 975 samples prepared from the collected material were subjected to compressive and bending strength tests. The samples for the tests were collected from carefully selected locations in order to discover the strength properties of various parts of the vertebral body.

FINDINGS

In the case of sample 2 (located in the posterior part of the vertebra, at mid-height) the stress values were the lowest and there were statistically significant differences compared to other samples. Moreover the value of compressive force in this case was lower for vertebrae with higher numbers. Top and bottom samples demonstrated statistically significant higher mean values of destructive stress. In terms of the bending strength test, the mean value of destructive stress in all lumbar vertebrae for all samples increased for vertebrae with higher numbers.

INTERPRETATION

The spongy tissue in healthy vertebral bodies has a very heterogeneous structure. This may be due to the presence of the nutrient canal and the arc structure allowing more springy movement and improved transfer of loads by the vertebral body.

Effect of disc degeneration on the mechanical behavior of the human lumbar spine: a probabilistic finite element study

BACKGROUND CONTEXT

Intervertebral disc degeneration has been subject to numerous in vivo and in vitro investigations and numerical studies during recent decades, reporting partially contradictory findings. However, most of the previous studies were limited in the number of specimens investigated and, therefore, could not consider the vast variety of the specimen geometries, which are likely to strongly influence the mechanical behavior of the spine.

PURPOSE

To complement the understanding of the mechanical consequences of disc degeneration, whereas considering natural variations in the major spinal geometrical parameters.

DESIGN/SETTING

A probabilistic finite element study.

METHODS

A parametric finite element model of a human L4-L5 motion segment considering 40 geometrical parameters was developed. One thousand individual geometries comprising four degeneration grades were generated in a probabilistic manner, and the influence of the severity of disc degeneration on the mechanical response of the motion segment to different loading conditions was statistically evaluated.

RESULTS

Variations in the individual structural parameters resulted in marked variations in all evaluated parameters within each degeneration grade. Nevertheless, the effect of degeneration in almost all evaluated response values was statistically significant. With degeneration, the intradiscal pressure progressively decreased. At the same time, the facet loads increased and the ligament tension was reduced. The initially nonlinear load-deformation relationships became linear whereas the segment stiffness increased.

CONCLUSIONS

Results indicate significant stiffening of the motion segment with progressing degeneration and gradually increasing loading of the facets from nondegenerated to moderately degenerated conditions along with a significant reduction of the ligament tension in flexion.

Prediction of in vivo lower cervical spinal loading using musculoskeletal multi-body dynamics model during the head flexion/extension, lateral bending and axial rotation

Cervical spine diseases lead to a heavy economic burden to the individuals and societies. Moreover, frequent post-operative complications mean a higher risk of neck pain and revision. At present, controversy still exists for the etiology of spinal diseases and their associated complications. Knowledge of in vivo cervical spinal loading pattern is proposed to be the key to answer these questions. However, direct acquisition of in vivo cervical spinal loading remains challenging. In this study, a previously developed cervical spine musculoskeletal multi-body dynamics model was utilized for spinal loading prediction. The in vivo dynamic segmental contributions to head motion and the out-of-plane coupled motion were both taken into account. First, model validation and sensitivity analysis of different segmental contributions to head motion were performed. For model validation, the predicted intervertebral disk compressive forces were converted into the intradiskal pressures and compared with the published experimental measurements. Significant correlations were found between the predicted values and the experimental results. Thus, the reliability and capability of the cervical spine model was ensured. Meanwhile, the sensitivity analysis indicated that cervical spinal loading is sensitive to different segmental contributions to head motion. Second, the compressive, shear and facet joint forces at C3–C6 disk levels were predicted, during the head flexion/extension, lateral bending and axial rotation. Under the head flexion/extension movement, asymmetric loading patterns of the intervertebral disk were obtained. In comparison, symmetrical typed loading patterns were found for the head lateral bending and axial rotation movements. However, the shear forces were dramatically increased during the head excessive extension and lateral bending. Besides, a nonlinear correlation was seen between the facet joint force and the angular displacement. In conclusion, dynamic cervical spinal loading was both intervertebral disk angle-dependent and level-dependent. Cervical spine musculoskeletal multi-body dynamics model provides an attempt to comprehend the in vivo biomechanical surrounding of the human head-neck system.

Prediction of Cervical Spinal Joint Loading and Secondary Motion Using a Musculoskeletal Multibody Dynamics Model Via Force-Dependent Kinematics Approach

STUDY DESIGN

A cervical spine biomechanical investigation using multibody dynamics.

OBJECTIVE

To develop a comprehensive cervical spine multibody dynamics model incorporated with the force-dependent kinematics (FDK) approach, and to study the influence of soft tissue deformation on the joint loading prediction.

SUMMARY OF BACKGROUND DATA

Musculoskeletal multibody dynamics models have been widely used to analyze joint loading. Current cervical spine musculoskeletal models, however, neglect the joint internal motion caused by soft tissue deformation. A novel FDK approach is introduced, which can predict joint internal motion and spinal joint loading simultaneously.

METHODS

A comprehensive cervical spine musculoskeletal model with the posterior facet joints and essential ligaments was developed. To quantify the influence of soft tissue structures on joint loading prediction, four different models with different features were created. These newly developed models were validated, under flexion-extension movement. The predicted intervertebral disc loads (from C3-C4 to C5-C6) were compared with the published cadaveric experimental results. Moreover, the predicted facet joint forces, ligament forces, and anterior-posterior translations of instantaneous centers of rotation were also studied.

RESULTS

The obtained intervertebral disc loads were varied among different models. Model 3 provided the closest prediction of joint loading to the experimental results. Moreover, the facet joint and ligament forces were in similar range of magnitude as literature findings. The predicted instantaneous centers of rotation translational changes were in accordance with the in vivo kinematics observation.

CONCLUSION

In the present study, a validated cervical spine musculoskeletal model was developed, using multibody dynamics and FDK approach. It can simulate the function of musculature and consider joint internal motion, and thus provides more reliable joint loading prediction. This newly developed cervical model can be used as an efficient tool to study the biomechanical behaviors of human cervical spine, and to understand the fundamental pathologies of spinal pains.

LEVEL OF EVIDENCE

N /A.

A Cervico-Thoraco-Lumbar Multibody Dynamic Model for the Estimation of Joint Loads and Muscle Forces

Computational musculoskeletal models have been developed to predict mechanical joint loads on the human spine, such as the forces and moments applied to vertebral and facet joints and the forces that act on ligaments and muscles because of difficulties in the direct measurement of joint loads. However, many whole-spine models lack certain elements. For example, the detailed facet joints in the cervical region or the whole spine region may not be implemented. In this study, a detailed cervico-thoraco-lumbar multibody musculoskeletal model with all major ligaments, separated structures of facet contact and intervertebral disk joints, and the rib cage was developed. The model was validated by comparing the intersegmental rotations, ligament tensile forces, facet joint contact forces, compressive and shear forces on disks, and muscle forces were to those reported in previous experimental and computational studies both by region (cervical, thoracic, or lumbar regions) and for the whole model. The comparisons demonstrated that our whole spine model is consistent with in vitro and in vivo experimental studies and with computational studies. The model developed in this study can be used in further studies to better understand spine structures and injury mechanisms of spinal disorders.

Development and Validation of a Musculoskeletal Model of the Fully Articulated Thoracolumbar Spine and Rib Cage

We developed and validated a fully articulated model of the thoracolumbar spine in opensim that includes the individual vertebrae, ribs, and sternum. To ensure trunk muscles in the model accurately represent muscles in vivo, we used a novel approach to adjust muscle cross-sectional area (CSA) and position using computed tomography (CT) scans of the trunk sampled from a community-based cohort. Model predictions of vertebral compressive loading and trunk muscle tension were highly correlated to previous in vivo measures of intradiscal pressure (IDP), vertebral loading from telemeterized implants and trunk muscle myoelectric activity recorded by electromyography (EMG).

On the Use of Biaxial Properties in Modeling Annulus as a Holzapfel-Gasser-Ogden Material

Besides the biology, stresses and strains within the tissue greatly influence the location of damage initiation and mode of failure in an intervertebral disk. Finite element models of a functional spinal unit (FSU) that incorporate reasonably accurate geometry and appropriate material properties are suitable to investigate such issues. Different material models and techniques have been used to model the anisotropic annulus fibrosus, but the abilities of these models to predict damage initiation in the annulus and to explain clinically observed phenomena are unclear. In this study, a hyperelastic anisotropic material model for the annulus with two different sets of material constants, experimentally determined using uniaxial and biaxial loading conditions, were incorporated in a 3D finite element model of a ligamentous FSU. The purpose of the study was to highlight the biomechanical differences (e.g., intradiscal pressure, motion, forces, stresses, strains, etc.) due to the dissimilarity between the two sets of material properties (uniaxial and biaxial). Based on the analyses, the biaxial constants simulations resulted in better agreements with the in vitro and in vivo data, and thus are more suitable for future damage analysis and failure prediction of the annulus under complex multiaxial loading conditions.

In vivo morphological features of human lumbar discs

Recent biomechanics studies have revealed distinct kinematic behavior of different lumbar segments. The mechanisms behind these segment-specific biomechanical features are unknown. This study investigated the in vivo geometric characteristics of human lumbar intervertebral discs. Magnetic resonance images of the lumbar spine of 41 young Chinese individuals were acquired. Disc geometry in the sagittal plane was measured for each subject, including the dimensions of the discs, nucleus pulposus (NP), and annulus fibrosus (AF). Segmental lordosis was also measured using the Cobb method.In general, the disc length increased from upper to lower lumbar levels, except that the L4/5 and L5/S1 discs had similar lengths. The L4/5 NP had a height of 8.6±1.3 mm, which was significantly higher than all other levels (P<0.05). The L5/S1 NP had a length of 21.6±3.1 mm, which was significantly longer than all other levels (P<0.05). At L4/5, the NP occupied 64.0% of the disc length, which was significantly less than the NP of the L5/S1 segment (72.4%) (P<0.05). The anterior AF occupied 20.5% of the L4/5 disc length, which was significantly greater than that of the posterior AF (15.6%) (P<0.05). At the L5/S1 segment, the anterior and posterior AFs were similar in length (14.1% and 13.6% of the disc, respectively). The height to length (H/L) ratio of the L4/5 NP was 0.45±0.06, which was significantly greater than all other segments (P<0.05). There was no correlation between the NP H/L ratio and lordosis. Although the lengths of the lower lumbar discs were similar, the geometry of the AF and NP showed segment-dependent properties. These data may provide insight into the understanding of segment-specific biomechanics in the lower lumbar spine. The data could also provide baseline knowledge for the development of segment-specific surgical treatments of lumbar diseases.

Measures of internal lumbar load in professional drivers - the use of a whole-body finite-element model for the evaluation of adverse health effects of multi-axis vibration

The present study aimed to (1) employ the method for evaluation of vibration containing multiple shocks according to ISO/CD 2631-5:2014 (Model 1) and DIN SPEC 45697:2012 in a cohort of 537 professional drivers, (2) deliver the results for a re-analysis of epidemiological data obtained in the VIBRISKS study, (3) clarify the extent to which vibration acceleration and individual variables influence risk values, such as the daily compressive dose S(ed) and the risk factor R, and (4) compare the results with in vivo measurements and those obtained in previous studies with similar models. The risk factor R was influenced by the acceleration, lifetime exposure duration, sitting posture, age at the start of exposure and body mass/body mass index in order of decreasing effect. Age and annual and daily exposure duration had only a marginal effect. The daily compressive dose S(ed) and the risk factor R showed weak linear association with the daily vibration exposure A(8) and the vibration dose value VDV. The study revealed high shear forces in the lumbar spine.

PRACTITIONER SUMMARY

In a re-analysis of an epidemiological study of professional drivers, a software tool available with standards DIN SPEC 45697:2012 and ISO/CD 2631–5:2014 Model 1 was used to calculate the risk to the lumbar spine in terms of daily compressive dose S(ed) and risk factor R. The tool was found to be suitable for risk assessment in a large cohort.

Activities of everyday life with high spinal loads

Activities with high spinal loads should be avoided by patients with back problems. Awareness about these activities and knowledge of the associated loads are important for the proper design and pre-clinical testing of spinal implants. The loads on an instrumented vertebral body replacement have been telemetrically measured for approximately 1000 combinations of activities and parameters in 5 patients over a period up to 65 months postoperatively. A database containing, among others, extreme values for load components in more than 13,500 datasets was searched for 10 activities that cause the highest resultant force, bending moment, torsional moment, or shear force in an anatomical direction. The following activities caused high resultant forces: lifting a weight from the ground, forward elevation of straight arms with a weight in hands, moving a weight laterally in front of the body with hanging arms, changing the body position, staircase walking, tying shoes, and upper body flexion. All activities have in common that the center of mass of the upper body was moved anteriorly. Forces up to 1650 N were measured for these activities of daily life. However, there was a large intra- and inter-individual variation in the implant loads for the various activities depending on how exercises were performed. Measured shear forces were usually higher in the posterior direction than in the anterior direction. Activities with high resultant forces usually caused high values of other load components.

Comparison of eight published static finite element models of the intact lumbar spine: predictive power of models improves when combined together

Finite element (FE) model studies have made important contributions to our understanding of functional biomechanics of the lumbar spine. However, if a model is used to answer clinical and biomechanical questions over a certain population, their inherently large inter-subject variability has to be considered. Current FE model studies, however, generally account only for a single distinct spinal geometry with one set of material properties. This raises questions concerning their predictive power, their range of results and on their agreement with in vitro and in vivo values. Eight well-established FE models of the lumbar spine (L1-5) of different research centers around the globe were subjected to pure and combined loading modes and compared to in vitro and in vivo measurements for intervertebral rotations, disc pressures and facet joint forces. Under pure moment loading, the predicted L1-5 rotations of almost all models fell within the reported in vitro ranges, and their median values differed on average by only 2° for flexion-extension, 1° for lateral bending and 5° for axial rotation. Predicted median facet joint forces and disc pressures were also in good agreement with published median in vitro values. However, the ranges of predictions were larger and exceeded those reported in vitro, especially for the facet joint forces. For all combined loading modes, except for flexion, predicted median segmental intervertebral rotations and disc pressures were in good agreement with measured in vivo values. In light of high inter-subject variability, the generalization of results of a single model to a population remains a concern. This study demonstrated that the pooled median of individual model results, similar to a probabilistic approach, can be used as an improved predictive tool in order to estimate the response of the lumbar spine.

Discrepancies in anthropometric parameters between different models affect intervertebral rotations when loading finite element models with muscle forces from inverse static analyses

In only a few published finite element (FE) simulations have muscle forces been applied to the spine. Recently, muscle forces determined using an inverse static (IS) model of the spine were transferred to a spinal FE model, and the effect of methodical parameters was investigated. However, the sensitivity of anthropometric differences between FE and IS models, such as body height and spinal orientation, was not considered. The aim of this sensitivity study was to determine the influence of those differences on the intervertebral rotations (IVRs) following the transfer of muscle forces from an IS model to a FE model. Muscle forces were estimated for 20° flexion and 10° extension of the upper body using an inverse static musculoskeletal model. These forces were subsequently transferred to a nonlinear FE model of the spino-pelvic complex, which includes 243 muscle fascicles. Deviations of body height (±10 cm), spinal orientation in the sagittal plane (±10°), and body weight (±10 kg) between both models were intentionally generated, and their influences on IVRs were determined. The changes in each factor relative to their corresponding reference value of the IS model were calculated. Deviations in body height, spinal orientation, and body weight resulted in maximum changes in the IVR of 19.2%, 26% and 4.2%, respectively, relative to T12-S1 IVR. When transferring muscle forces from an IS to a FE model, it is crucial that both models have the same spinal orientation and height. Additionally, the body weight should be equal in both models.