Computational Biomechanical Modeling of Scoliotic Spine: Challenges and Opportunities

Characterization of spine and torso stiffness via differentiable biomechanics

We present a methodology to personalize the stiffness response of a biomechanical model of the torso and the spine. In high contrast to previous work, the proposed methodology uses controlled force-deformation data that mimic the conditions of spinal bracing for scoliosis, which leads to personalized biomechanical models that are suitable for computational brace design. The novel methodology relies on several technical contributions. First, a prototype system that includes controlled force measurement and low-dose radiographs, with low-encumbrance for its implementation in the clinical protocol. Second, a model of differentiable biomechanics of the torso and the spine, which becomes the key building block for robust parameter estimation. And third, an optimization procedure for parameter estimation from force-deformation data, which relies on differentiability of the biomechanics and the image generation process. We demonstrate the application of the methodology to a cohort of 7 subjects who underwent scoliosis check-ups, and we show quantitative validation of the estimated personalized parameters and the improvement over default parameters from the bibliography.

Real-time prediction of postoperative spinal shape with machine learning models trained on finite element biomechanical simulations

PURPOSE

Adolescent idiopathic scoliosis is a chronic disease that may require correction surgery. The finite element method (FEM) is a popular option to plan the outcome of surgery on a patient-based model. However, it requires considerable computing power and time, which may discourage its use. Machine learning (ML) models can be a helpful surrogate to the FEM, providing accurate real-time responses. This work implements ML algorithms to estimate post-operative spinal shapes.

METHODS

The algorithms are trained using features from 6400 simulations generated using the FEM from spine geometries of 64 patients. The features are selected using an autoencoder and principal component analysis. The accuracy of the results is evaluated by calculating the root-mean-squared error and the angle between the reference and predicted position of each vertebra. The processing times are also reported.

RESULTS

A combination of principal component analysis for dimensionality reduction, followed by the linear regression model, generated accurate results in real-time, with an average position error of 3.75 mm and orientation angle error below 2.74 degrees in all main 3D axes, within 3 ms. The prediction time is considerably faster than simulations based on the FEM alone, which require seconds to minutes.

CONCLUSION

It is possible to predict post-operative spinal shapes of patients with AIS in real-time by using ML algorithms as a surrogate to the FEM. Clinicians can compare the response of the initial spine shape of a patient with AIS to various target shapes, which can be modified interactively. These benefits can encourage clinicians to use software tools for surgical planning of scoliosis.

A Review of Finite Element Modeling for Anterior Cervical Discectomy and Fusion

The cervical spine poses many complex challenges that require complex solutions. Anterior cervical discectomy and fusion (ACDF) has been one such technique often employed to address such issues. In order to address the problems with ACDF and assess the modifications that have been made to the technique over time, finite element analyses (FEA) have proven to be an effective tool. The variations of cervical spine FEA models that have been produced over the past couple of decades, particularly more recent representations of more complex geometries, have not yet been identified and characterized in any literature. Our objective was to present material property models and cervical spine models for various simulation purposes. The outlining and refinement of the FEA process will yield more reliable outcomes and provide a stable basis for the modeling protocols of the cervical spine.

The application of finite element analysis to determine the optimal UIV of growing-rod treatment in early-onset scoliosis

Objectives: To analyze the stress distribution in the proximal vertebral body and soft tissue of dual growing-rod (GR) with different upper instrumented vertebra (UIV) to determine the optimal UIV.

Methods: A ten-year-old male EOS case treated with GR was selected. Based on spiral computed tomography (CT) scanning performed in 0.6 mm thick slices, a finite element model (FEM) of the preoperative state (M0, the original spine state) of the patient was created. Subsequently, four models with different UIV fixations were numerically analyzed by FEM, including M1 (UIV = T1, i.e., the upper-end vertebrae (UEV) of the upper thoracic curve), M2 (UIV = T2), M3 (UIV = T3) and M4 (UIV = T4, i.e., the lower end vertebrae (LEV) of the upper thoracic curve). Displacement and maximum stress in the proximal vertebral body and soft tissue were measured and compared among the five models.

Results: The spine model was fixed with the sacrum, and the gravity conditions were imposed on each vertebral body according to the research of Clin and Pearsall. The results are as follows:M4 model has the largest overall displacement, while M1 has the least displacement among the four models. Except M2, the maximum normalized stress of UIV increases with the downward movement of UIV. M1 has the lowerest annulus fibrosus stress and highest joint capsule stress, which is characterized by the vertebrae backward leaning, while M4 is the opposite. The supraspinous ligament stress of M3 and M4 is significantly higher than that of M1 and M2. This suggests that UIV downshift increases the tendency of the proximal vertebral bodies to bend forward, thereby increasing the tension of the posterior ligaments (PL).

Conclusion: The UIV of the GR is recommended to be close to the UEV of the upper thoracic curve, which can reduce the stress of the proximal PL, thereby reducing the occurrence of proximal junctional kyphosis (PJK).

Active Soft Brace for Scoliotic Spine: A Finite Element Study to Evaluate in-Brace Correction

Scoliosis is a spinal disorder that is conventionally treated using rigid or soft braces. Computational methods such as finite element-based models are used to investigate the mechanics of the spine and the effect of braces. Most spinal braces are either passive, static, or rigid and do not allow mobility to the spine, resulting in muscle atrophy, skin deterioration and other spine complexities. Lack of control over the amount of force being exerted by braces on the human spine could have adverse effects. Therefore, developing an active soft brace which allows mobility to the spine while applying controlled corrective forces could be a promising solution. This study presents finite element analysis (FEA) of an active soft brace that applies corrective forces using elastic bands. The pressure exerted by the brace on the spine can be controlled by varying the tensions in the elastic bands. The elastic band tensions are controlled using low-power, lightweight, and twisted string actuators (TSAs). This study aims to demonstrate the immediate corrections induced by the soft active brace using a scoliotic spine finite element (FE) model. A FE model of the patient’s trunk was created and validated with in vitro study. The brace model was installed on the simulated trunk to evaluate in-brace correction in both sagittal and coronal planes. The brace was evaluated under various load cases by simulating the actuator action.

Applications and limitations of current markerless motion capture methods for clinical gait biomechanics

BACKGROUND

Markerless motion capture has the potential to perform movement analysis with reduced data collection and processing time compared to marker-based methods. This technology is now starting to be applied for clinical and rehabilitation applications and therefore it is crucial that users of these systems understand both their potential and limitations. This literature review aims to provide a comprehensive overview of the current state of markerless motion capture for both single camera and multi-camera systems. Additionally, this review explores how practical applications of markerless technology are being used in clinical and rehabilitation settings, and examines the future challenges and directions markerless research must explore to facilitate full integration of this technology within clinical biomechanics.

METHODOLOGY

A scoping review is needed to examine this emerging broad body of literature and determine where gaps in knowledge exist, this is key to developing motion capture methods that are cost effective and practically relevant to clinicians, coaches and researchers around the world. Literature searches were performed to examine studies that report accuracy of markerless motion capture methods, explore current practical applications of markerless motion capture methods in clinical biomechanics and identify gaps in our knowledge that are relevant to future developments in this area.

RESULTS

Markerless methods increase motion capture data versatility, enabling datasets to be re-analyzed using updated pose estimation algorithms and may even provide clinicians with the capability to collect data while patients are wearing normal clothing. While markerless temporospatial measures generally appear to be equivalent to marker-based motion capture, joint center locations and joint angles are not yet sufficiently accurate for clinical applications. Pose estimation algorithms are approaching similar error rates of marker-based motion capture, however, without comparison to a gold standard, such as bi-planar videoradiography, the true accuracy of markerless systems remains unknown.

CONCLUSIONS

Current open-source pose estimation algorithms were never designed for biomechanical applications, therefore, datasets on which they have been trained are inconsistently and inaccurately labelled. Improvements to labelling of open-source training data, as well as assessment of markerless accuracy against gold standard methods will be vital next steps in the development of this technology.

Model-Based Radiation-Free Assessment of Scoliosis: A Principle Validation Study

Purpose

Adolescent scoliosis is one of the common pediatric spinal diseases which has a high risk of progression due to the rapid growth of the skeleton during the growing stage therefore needs regular clinical monitoring including X-rays. Because X-rays could lead to ionizing radiation-related health problems, an ionizing radiation-free, non-invasive method is presented here to estimate the degree of scoliosis and to potentially support the medical assessment.

Methods

The radiation-free body scanner provides a 3D surface scan of the torso. A basic 3D structure of the human ribcage and vertebral column was modeled and simulated with computer-aided design software and finite element method calculation. For comparison with X-rays, courses of vertebral columns derived from 3D torso images and 3D models were analyzed with respect to their apex positions and angles.

Results

The methods show good results in the estimation of the apex positions of scoliosis. Strong correlations (R = 0.8924) were found between the apex and Cobb angle from X-rays. Similar correlations (R = 0.8087) was obtained between the apex angles extracted from X-rays and the combination of torso scan images with 3D model simulations. Promising agreement was obtained between the spinal trajectories extracted from X-ray and 3D torso images.

Conclusions

Very strong correlations suggest that the apex angle could potentially be used for scoliosis assessment in follow-up examinations in complement to the Cobb angle. However, further improvements of the methods and tests on a larger number of data set are necessary before their introduction into the clinical application.

Computational modelling of the scoliotic spine: A literature review

Scoliosis is a deformity of the spine that in severe cases requires surgical treatment. There is still disagreement among clinicians as to what the aim of such treatment is as well as the optimal surgical technique. Numerical models can aid clinical decision-making by estimating the outcome of a given surgical intervention. This paper provided some background information on the modelling of the healthy spine and a review of the literature on scoliotic spine models, their validation, and their application. An overview of the methods and techniques used to construct scoliotic finite element and multibody models was given as well as the boundary conditions used in the simulations. The current limitations of the models were discussed as well as how such limitations are addressed in non-scoliotic spine models. Finally, future directions for the numerical modelling of scoliosis were addressed.

How Can Biomechanical Multibody Models of Scoliosis Be Accurate in Simulating Spine Movement Behavior While Neglecting the Changes of Spinal Length?

This paper studies how biomechanical multibody models of scoliosis can neglect the changes of spinal length and yet be accurate in reconstructing spinal columns. As these models with fixed length comprise rigid links interconnected by rotary joints, they resemble polygonal chains that approximate spine curves with a finite number of line segments. In mathematics, using more segments with shorter lengths can result in more accurate curve approximations. This raises the question of whether more accurate spine curve approximations by increasing the number of links/joints can yield more accurate spinal column reconstructions. For this, the accuracy of spine curve approximation was improved consistently by increasing the number of links/joints, and its effects on the accuracy of spinal column reconstruction were assessed. Positive correlation was found between the accuracy of spine reconstruction and curve approximation. It was shown that while increasing the accuracy of curve approximations, the representation of scoliosis concavity and its side-to-side deviations were improved. Moreover, reconstruction errors of the spine regions separated by the inflection vertebrae had minimal impacts on each other. Overall, multibody scoliosis models with fixed spinal lengths can benefit from the extra rotational joints that contribute toward the accuracy of spine curve approximation. The outcome of this study leads to concurrent accuracy improvement and simplification of multibody models; joint-link configurations can be independently defined for the regions separated by the inflection vertebrae, enabling local optimization of the models for higher accuracy without unnecessary added complexity to the whole model.

A one-dimensional collagen-based biomechanical model of passive soft tissue with viscoelasticity and failure

INTRODUCTION

Strains and sprains of soft tissues, including tendons and ligaments, are frequently occurring injuries. Musculoskeletal models show great promise in prediction and prevention of these injuries. However, these models rarely account for the viscoelastic properties of ligaments and tendons, much less their failure properties. The purpose of this project was to develop, simplify, and analyze a collagen-distribution model to address these limitations.

MODEL DEVELOPMENT

A distribution-moment approximation was applied to an existing partial differential equation model to reduce its computational complexity. The resulting model was equipped with a Voigt model in series, which endowed it with viscoelastic properties in addition to failure properties.

RESULTS

The model was able to reproduce the characteristic toe, linear, and failure regions ubiquitous throughout in-vitro tests on tissue specimens. In addition, it was able to reproduce a tri-phasic creep test consisting of an initial deformation, a steady-state, and failure. Stress-relaxation and hysteresis were also reproducible by the model.

DISCUSSION AND CONCLUSION

The ability to reproduce so many characteristics of biological tissues suggests more bio-fidelity was achieved by the reduced model was other currently available models. Future work to further improve its bio-fidelity is proposed for specific tendons and ligaments.

A mechanical analog thoracolumbar spine model for the evaluation of scoliosis bracing technology

Introduction

Thoracolumbar braces are used to treat Adolescent Idiopathic Scoliosis. The objective of this study was to design and validate a mechanical analog model of the spine to simulate a thoracolumbar, single-curve, scoliotic deformity in order to quantify brace structural properties and corrective force response on the spine.

Methods

The Scoliosis Analog Model used a linkage-based system to replicate 3D kinematics of spinal correction observed in the clinic. The Scoliosis Analog Model is used with a robotic testing platform and programmed to simulate Cobb angle and axial rotation correction while equipped with a brace. The 3D force and moment responses generated by the brace in reaction to the simulated deformity were measured by six-axis load cells.

Results

Validation of the model’s force transmission showed less than 6% loss in the force analysis due to assembly friction. During simulation of 10° Cobb angle and 5° axial rotation correction, the brace applied 101 N upwards and 67 N inwards to the apical connector of the model. Brace stiffness properties were 0.5–0.6 N/° (anteroposterior), 0.5–2.3 N/° (mediolateral), 23.3–26.5 N/° (superoinferior), and 0.6 Nm/° (axial rotational).

Conclusions

The Scoliosis Analog Model was developed to provide first time measures of the multidirectional forces applied to the spine by a thoracolumbar brace. This test assembly could be used as a future design and testing tool for scoliosis brace technology.

Development and validation of a geometrically personalized finite element model of the lower ligamentous cervical spine for clinical applications

Epidemiological and clinical studies show that the magnitude and scope of cervical disease are on the rise, along with the world's rising aging population. From a biomechanical perspective, the cervical spine presents a wide inter-individual variability, where its motion patterns and load sharing strongly depend on the anatomy. This study aimed to first develop and validate a geometrically patient-specific model of the lower cervical spine for clinical applications, and secondly to use the model to investigate the spinal biomechanics associated with typical cervical disorders. Based on measurements of 30 parameters from X-ray radiographs, the 3D geometry of the vertebrae and intervertebral discs (IVDs) were developed, and detailed finite element models (FEMs) of the lower ligamentous cervical spine for 6 subjects were constructed and simulated. The models were then used for the investigation of different grades of IVD alteration. The multi directional range of motion (ROM) results were in alignment with the in-vitro and in-Silico studies confirming the validity of the model. Severe disc alteration (Grade 3) presented a significant decrease in the ROM and intradiscal pressure (flexion, extension, and axial rotation) on the C5-C6 and slightly increase on the adjacent levels. Maximum stress in Annulus Fibrosus (AF) and facet joint forces increased for Grade 3 for both altered and adjacent levels. The novel validated geometrically-personalized FEM presented in this study potentially offers the clinical community a valuable quantitative tool for the noninvasive analyses of the biomechanical alterations associated with cervical spine disease towards improved surgical planning and enhanced clinical outcomes.

Validation of the AnyBody full body musculoskeletal model in computing lumbar spine loads at L4L5 level

In the panorama of available musculoskeletal modeling software, AnyBody software is a commercial tool that provides a full body musculoskeletal model which is increasingly exploited by numerous researchers worldwide. In this regard, model validation becomes essential to guarantee the suitability of the model in representing the simulated system. When focusing on lumbar spine, the previous works aimed at validating the AnyBody model in computing the intervertebral loads held several limitations, and a comprehensive validation is to be considered as lacking. The present study was aimed at extensively validating the suitability of the AnyBody model in computing lumbar spine loads at L4L5 level. The intersegmental loads were calculated during twelve specific exercise tasks designed to accurately replicate the conditions during which Wilke et al. (2001) measured in vivo the L4L5 intradiscal pressure. Motion capture data of one volunteer subject were acquired during the execution of the tasks and then imported into AnyBody to set model kinematics. Two different approaches in computing intradiscal pressure from the intersegmental load were evaluated. Lumbopelvic rhythm was compared with reference in vivo measurements to assess the accuracy of the lumbopelvic kinematics. Positive agreement was confirmed between the calculated pressures and the in vivo measurements, thus demonstrating the suitability of the AnyBody model. Specific caution needs to be taken only when considering postures characterized by large lateral displacements. Minor discrepancy was found assessing lumbopelvic rhythm. The present findings promote the AnyBody model as an appropriate tool to non-invasively evaluate the lumbar loads at L4L5 in physiological activities.

Semiautomated 3D Spine Reconstruction from Biplanar Radiographic Images: Prediction of Intervertebral Loading in Scoliotic Subjects

The present study proposes a semiautomatic software approach to reconstruct 3D subject-specific musculoskeletal model of thoracolumbar spine from radiographic digitized images acquired with EOS system. The approach is applied to evaluate the intervertebral loads in 38 standing adolescents with mild idiopathic scoliosis. For each vertebra, a set of landmarks was manually identified on radiographic images. The landmark coordinates were processed to calculate the following vertebral geometrical properties in the 3D space (i) location (ii) dimensions; and (iii) rotations. Spherical joints simulated disks, ligaments, and facet joints. Body weight distribution, muscles forces, and insertion points were placed according to physiological-anatomical values. Inverse static analysis, calculating joints' reactions in maintaining assigned spine configuration, was performed with AnyBody software. Reaction forces were computed to quantify intervertebral loads, and correlation with the patient anatomical parameters was then checked. Preliminary validation was performed comparing the model outcomes with that obtained from other authors in previous modeling works and from in vivo measurements. The comparison with previous modeling works and in vivo studies partially fulfilled the preliminary validation purpose. However, minor incongruities were pointed out that need further investigations. The subjects' intervertebral loads were found significantly correlated with the anatomical parameters in the sagittal and axial planes. Despite preliminary encouraging results that support model suitability, future investigations to consolidate the proposed approach are necessary. Nonetheless, the present method appears to be a promising tool that once fully validated could allow the subject-specific non-invasive evaluation of a deformed spine, providing supplementary information to the routine clinical examination and surgical intervention planning.

Evaluation of a Patient-Specific Finite-Element Model to Simulate Conservative Treatment in Adolescent Idiopathic Scoliosis

STUDY DESIGN

Retrospective validation study.

OBJECTIVES

To propose a method to evaluate, from a clinical standpoint, the ability of a finite-element model (FEM) of the trunk to simulate orthotic correction of spinal deformity and to apply it to validate a previously described FEM.

SUMMARY OF BACKGROUND DATA

Several FEMs of the scoliotic spine have been described in the literature. These models can prove useful in understanding the mechanisms of scoliosis progression and in optimizing its treatment, but their validation has often been lacking or incomplete.

METHODS

Three-dimensional (3D) geometries of 10 patients before and during conservative treatment were reconstructed from biplanar radiographs. The effect of bracing was simulated by modeling displacements induced by the brace pads. Simulated clinical indices (Cobb angle, T1-T12 and T4-T12 kyphosis, L1-L5 lordosis, apical vertebral rotation, torsion, rib hump) and vertebral orientations and positions were compared to those measured in the patients' 3D geometries.

RESULTS

Errors in clinical indices were of the same order of magnitude as the uncertainties due to 3D reconstruction; for instance, Cobb angle was simulated with a root mean square error of 5.7°, and rib hump error was 5.6°. Vertebral orientation was simulated with a root mean square error of 4.8° and vertebral position with an error of 2.5 mm.

CONCLUSIONS

The methodology proposed here allowed in-depth evaluation of subject-specific simulations, confirming that FEMs of the trunk have the potential to accurately simulate brace action. These promising results provide a basis for ongoing 3D model development, toward the design of more efficient orthoses.