Spine Model Research Articles

Implications of range of motion requirements for the laxity of ligaments in a lumbar finite element model

Finite element models of the lumbar spine often adopt ligament properties from tensile tests without accounting for possible differences between testing and in situ initial ligament length. Such differences could result in laxities or preloads at the beginning of a simulation that would affect the ligament forces, tangent stiffness, and the posture at which they fail. In vivo and in vitro human experimental data reported laxities or preloads. However, laxities or preloads, which could also result from postural differences, are often neglected in simulation studies. This study proposes a numerical methodology to identify ranges of ligament laxities or preloads compatible with the selected tensile ligament properties, the model, and the range of motion (RoM) the model aims to simulate. The approach assumes that ligaments should remain in a safe elongation range for the complete RoM, and that each ligament should play a significant mechanical role in at least one load case. The methodology was applied to the functional spinal unit (FSU) models using the RoM from healthy subjects and ligament properties from the literature. Without laxity, some ligaments reached their elongation at failure within the RoM. Laxity ranges varied considerably (from −9.2 mm preload to 10.7 mm laxity) and flexion was the most critical load case to determine them. Their effect on the mobility response was also assessed. The effect on the mobility of a FSU was also assessed. While the proposed method cannot determine an exact laxity value, it is simple and it can be applied to any model to identify a plausible range of ligament initial length.

Trunk muscle forces and spinal loads during heavy deadlift: Effects of personalization, muscle wrapping, muscle lever arm, and lumbopelvic rhythm.

Heavy deadlift is used as a physical fitness screening tool in the U.S. Army. Despite the relevance of such a screening tool to military tasks performed by Service Members, the biomechanical impact of heavy deadlift and its risk of low-back injury remain unknown. A kinematics-driven musculoskeletal model of spine was implemented to investigate biomechanics of the lower back in a volunteer (23 years old, height of 1.82 m, and body mass of 98.8kg) during a 68 kg deadlift. In search of protective mechanisms, effects of model personalization and variations in trunk musculature and lumbopelvic rhythm were also investigated. The net moment, compression and shear forces at the L5-S1 reached peaks of 684 Nm, 17.2 and 4.2kN, respectively. Geometrical personalization and changes in lumbopelvic rhythm had the least effects on predictions while increases in muscle moment arms (40%) had the largest effects that caused, respectively, 32% and 36% decrease in the maximum compressive and shearing forces. Initiating wrapping of back muscles at farther distances from the spine had opposing effects on spinal loads; peak compression at the L5-S1 decreased by 12% whereas shear increased by 19%. Despite mechanisms considered, spinal loads during heavy deadlift exceed the existing evidence concerning the threshold of injury for spinal segments, suggesting the vulnerability to injury. Chronic exposure to such high-spinal loads may lead to (micro) fractures, degeneration, pathoanatomical changes and finally low-back pain.

Finite element modelling and biodynamic response prediction of the seated human body exposed to whole-body vibration

Biodynamic modelling of seat-occupant systems can assist in seat comfort design. A finite element (FE) model of the seated human body, including detailed modelling of the lumbar spine, was established to reflect the human response to vibration and biodynamic response of the lumbar spine under whole-body vibration (WBV). The lumbar spine model was established and validated against the in-vitro results and calculated data. The posture of the lumbar spine was adjusted according to the radiological research results, and the adjusted model was combined to establish a FE model of the seated human body. The present seated human model with backrest inclination angles of 10, 20, and 30°, validated by comparing the measured apparent mass and seat-to-lumbar spine transmissibility, was used to calculate the biodynamic response of the lumbar spine with three inclined backrests under WBV. The results showed that the model could characterise the apparent mass, seat-to-lumbar spine transmissibility, and the biodynamic response of the lumbar spine. Practitioner summary: Biodynamic models can represent dynamic characteristics of the human body exposed to vibration and assist in seat comfort design. The three-dimensional FE model of the human body can be used to explore the human response to vibration and the biodynamic response of the lumbar spine under WBV.

A biomechanical comparison of posterior fixation approaches in lumbar fusion using computed tomography based lumbosacral spine modelling

Extreme lateral interbody fusion (XLIF) may be performed with a standalone interbody cage, or with the addition of unilateral or bilateral pedicle screws; however, decisions regarding supplemental fixation are predominantly based on clinical indicators. This study examines the impact of posterior supplemental fixation on facet micromotions, cage loads and load-patterns at adjacent levels in a L4-L5 XLIF at early and late fusion stages. CT data from an asymptomatic subject were segmented into anatomical regions and digitally stitched into a surface mesh of the lumbosacral spine (L1-S1). The interbody cage and posterior instrumentation (unilateral and bilateral) were inserted at L4-L5. The volumetric mesh was imported into finite element software for pre-processing, running nonlinear static solves and post-processing. Loads and micromotions at the index-level facets reduced commensurately with the extent of posterior fixation accompanying the XLIF, while load-pattern changes observed at adjacent facets may be anatomically dependent. In flexion at partial fusion, compressive stress on the cage reduced by 54% and 72% in unilateral and bilateral models respectively; in extension the reductions were 58% and 75% compared to standalone XLIF. A similar pattern was observed at full fusion. Unilateral fixation provided similar stability compared to bilateral, however there was a reduction in cage stress-risers with the bilateral instrumentation. No changes were found at adjacent discs. Posterior supplemental fixation alters biomechanics at the index and adjacent levels in a manner that warrants consideration alongside clinical information. Unilateral instrumentation is a more efficient option where the stability requirements and subsidence risk are not excessive.

Biomechanical evaluation of two-level oblique lumbar interbody fusion combined with posterior four-screw fixation：A finite element analysis

ObjectiveBy constructing the three-dimensional finite element model of two-level OLIF lumbar spine, the aim of this study was to demonstrate the feasibility and effectiveness of posterior four-screw fixation for treatment of two-level lumbar degenerative diseases from the perspective of biomechanics. MethodsAn intact L3-S1 segment nonlinear lumbar finite element model (M0) was constructed from the CT scanning data of a healthy adult. After verification, two-level OLIF procedure were simulated, and three patterns of finite element analysis models were constructed: two-level stand-alone OLIF group (M1), two-level OLIF + four-screw fixation group (M2) and two-level OLIF + six-screw fixation group (M3). Range of motion, stress of the cage, and stress of fixation were evaluated in the different models. ResultsUnder various motion modes，the ROM of M2 and M3 were significantly lower than those of M1. The ROM reduction of M2 relative to M1 was much greater than that of M3 relative to M2. Moreover, the peak von Mises stresses of endplates in M2 were almost the same as those in M3. In terms of the maximum stresses of cages, M2 and M3 were essentially identical. Besides, the maximum stresses of posterior instrumentation in M2 and M3 were similar, which were mainly concentrated at the root of pedicle screws. ConclusionThere were no significant differences between M2 and M3 from the biomechanical analysis. In two-level OLIF, posterior four-screw fixation can replace six-screw fixation, which reduces surgical trauma and decreases economic burden of patients, and will be a cost-effective alternative.

Effects of cervical rotatory manipulation on the cervical spinal cord complex with ossification of the posterior longitudinal ligament in the vertebral canal: A finite element study.

Background: There are few studies focusing on biomechanism of spinal cord injury according to the ossification of the posterior longitudinal ligament (OPLL) during cervical rotatory manipulation (CRM). This study aimed to explore the biomechanical effects of CRM on the spinal cord, dura matter and nerve roots with OPLL in the cervical vertebral canal. Methods: Three validated FE models of the craniocervical spine and spinal cord complex were constructed by adding mild, moderate, and severe OPLL to the healthy FE model, respectively. We simulated the static compression of the spinal cord by OPLL and the dynamic compression during CRM in the flexion position. The stress distribution of the spinal cord complex was investigated. Results: The cervical spinal cord experienced higher von Mises stress under static compression by the severe OPLL. A higher von Mises stress was observed on the spinal cord in the moderate and severe OPLL models during CRM. The dura matter and nerve roots had a higher von Mises stress in all three models during CRM. Conclusion: The results show a high risk in performing CRM in the flexion position on patients with OPLL, in that different occupying ratios in the vertebral canal due to OPLL could significantly increase the stress on the spinal cord complex.

On the Finite Element Modeling of the Lumbar Spine: A Schematic Review

Finite element modelling of the lumbar spine is a challenging problem. Lower back pain is among the most common pathologies in the global populations, owing to which the patient may need to undergo surgery. The latter may differ in nature and complexity because of spinal disease and patient contraindications (i.e., aging). Today, the understanding of spinal column biomechanics may lead to better comprehension of the disease progression as well as to the development of innovative therapeutic strategies. Better insight into the spine’s biomechanics would certainly guarantee an evolution of current device-based treatments. In this setting, the computational approach appears to be a remarkable tool for simulating physiological and pathological spinal conditions, as well as for various aspects of surgery. Patient-specific computational simulations are constantly evolving, and require a number of validation and verification challenges to be overcome before they can achieve true and accurate results. The aim of the present schematic review is to provide an overview of the evolution and recent advances involved in computational finite element modelling (FEM) of spinal biomechanics and of the fundamental knowledge necessary to develop the best modeling approach in terms of trustworthiness and reliability.

Assessment of thoracic spinal curvatures in static postures using spatially tracked 3D ultrasound volumes: a proof-of-concept study

The assessment of spinal posture is a difficult endeavour given the lack of identifiable bony landmarks for placement of skin markers. Moreover, potentially significant soft tissue artefacts along the spine further affect the accuracy of marker-based approaches. The objective of this proof-of-concept study was to develop an experimental framework to assess spinal postures by using three-dimensional (3D) ultrasound (US) imaging. A phantom spine model immersed in water was scanned using 3D US in a neutral and two curved postures mimicking a forward flexion in the sagittal plane while the US probe was localised by three electromagnetic tracking sensors attached to the probe head. The obtained anatomical ‘coarse’ registrations were further refined using an automatic registration algorithm and validated by an experienced sonographer. Spinal landmarks were selected in the US images and validated against magnetic resonance imaging data of the same phantom through image registration. Their position was then related to the location of the tracking sensors identified in the acquired US volumes, enabling the localisation of landmarks in the global coordinate system of the tracking device. Results of this study show that localised 3D US enables US-based anatomical reconstructions comparable to clinical standards and the identification of spinal landmarks in different postures of the spine. The accuracy in sensor identification was 0.49 mm on average while the intra- and inter-observer reliability in sensor identification was strongly correlated with a maximum deviation of 0.8 mm. Mapping of landmarks had a small relative distance error of 0.21 mm (SD = ± 0.16) on average. This study implies that localised 3D US holds the potential for the assessment of full spinal posture by accurately and non-invasively localising vertebrae in space.

Cervical Spine Finite Element Models for Healthy Subjects: Development and Validation

Abstract Finite element modeling is a popular method for predicting kinematics and kinetics in spine biomechanics. With the advancement of powerful computational equipment, more detailed finite element models have been developed for the various spine segments. In this study, five detailed finite element models of the cervical spine are developed and validated. The geometric boundaries of the vertebrae are determined from computed tomography (CT) scans of five female subjects. The models include the C2–C7 vertebrae, intervertebral discs, nuclei, endplates, and five major ligaments (anterior longitudinal ligament (ALL), posterior longitudinal ligament (PLL), ligamentum flavum (LF), interspinous ligament (ISL), and capsular ligament (CL)). The ligaments follow nonlinear stress–strain curves whereas all other parts adopt linear material properties. All the material properties are taken from existing literature. The mesh convergence test is performed under flexion/extension. For flexion/extension motion, a pure moment is applied at the top surface of the odontoid process of the C2 vertebra while nodes at the bottom surface of the C7 vertebra are fixed in all directions. The models are extensively validated in flexion/extension, lateral bending, and axial rotation against experimental and finite element studies in the literature. Intervertebral rotation and range of motion are studied under different loading conditions found in the literature. This research also investigates intersubject variability for the cervical spine among five finite element models from five different subjects. Predicted angular displacements and ranges of motion of the current models are consistent with the literature. The validated models are expected to be applicable to simulate neck-related trauma like whiplash and high-g acceleration, among other scenarios.

Effects of geometric individualisation of a human spine model on load sharing: neuro-musculoskeletal simulation reveals significant differences in ligament and muscle contribution

In spine research, two possibilities to generate models exist: generic (population-based) models representing the average human and subject-specific representations of individuals. Despite the increasing interest in subject specificity, individualisation of spine models remains challenging. Neuro-musculoskeletal (NMS) models enable the analysis and prediction of dynamic motions by incorporating active muscles attaching to bones that are connected using articulating joints under the assumption of rigid body dynamics. In this study, we used forward-dynamic simulations to compare a generic NMS multibody model of the thoracolumbar spine including fully articulated vertebrae, detailed musculature, passive ligaments and linear intervertebral disc (IVD) models with an individualised model to assess the contribution of individual biological structures. Individualisation was achieved by integrating skeletal geometry from computed tomography and custom-selected muscle and ligament paths. Both models underwent a gravitational settling process and a forward flexion-to-extension movement. The model-specific load distribution in an equilibrated upright position and local stiffness in the L4/5 functional spinal unit (FSU) is compared. Load sharing between occurring internal forces generated by individual biological structures and their contribution to the FSU stiffness was computed. The main finding of our simulations is an apparent shift in load sharing with individualisation from an equally distributed element contribution of IVD, ligaments and muscles in the generic spine model to a predominant muscle contribution in the individualised model depending on the analysed spine level.

The Role of Multifidus in the Biomechanics of Lumbar Spine: A Musculoskeletal Modeling Study.

The role of multifidus in the biomechanics of lumbar spine remained unclear. This study aimed to investigate the role of multifidus in the modeling of lumbar spine and the influence of asymmetric multifidus atrophy on the biomechanics of lumbar spine. This study considered five different multifidus conditions in the trunk musculoskeletal models: group 1 (with entire multifidus), group 2 (without multifidus), group 3 (multifidus with half of maximum isometric force), group 4 (asymmetric multifidus atrophy on L5/S1 level), and group 5 (asymmetric multifidus atrophy on L4/L5 level). In order to test how different multifidus situations would affect the lumbar spine, four trunk flexional angles (0°, 30°, 60°, and 90°) were simulated. The calculation of muscle activation and muscle force was done using static optimization function in OpenSim. Then, joint reaction forces of L5/S1 and L4/L5 levels were calculated and compared among the groups. The models without multifidus had the highest normalized compressive forces on the L4/L5 level in trunk flexion tasks. In extreme cases produced by group 2 models, the normalized compressive forces on L4/L5 level were 444% (30° flexion), 568% (60° flexion), and 576% (90° flexion) of upper body weight, which were 1.82 times, 1.63 times, and 1.13 times as large as the values computed by the corresponding models in group 1. In 90° flexion, the success rate of simulation in group 2 was 49.6%, followed by group 3 (84.4%), group 4 (89.6%), group 5 (92.8%), and group 1 (92.8%). The results demonstrate that incorporating multifidus in the musculoskeletal model is important for increasing the success rate of simulation and decreasing the incidence of overestimation of compressive load on the lumbar spine. Asymmetric multifidus atrophy has negligible effect on the lower lumbar spine in the trunk flexion posture. The results highlighted the fine-tuning ability of multifidus in equilibrating the loads on the lower back and the necessity of incorporating multifidus in trunk musculoskeletal modeling.

3D Spine Model Reconstruction Based on RGBD Images of Unclothed Back Surface.

Most 3D spine reconstruction methods require X-ray images as input, which usually leads to high cost and radiation damage. Therefore, these methods are hardly ever applied to large scale scoliosis screening or spine pose monitoring during treatment. We propose a novel, low-cost, easy-to-operate and none-radioactive 3D spine model reconstruction method, which is based on human back surface information without requiring X-ray images as input. Our method fits a pre-built Spine Priors Model (SPrM) to human back surface information and reconstructs the main part of spine with 17 vertebrae: lumbar vertebrae L1-L5 and thoracic vertebrae T1-T12. The Spine Priors Model is constructed according to human spine priors, including Statistical Spine Shape Model (SSSM), Spine Pose Model (SPM) and Spine Biomechanical Simplified Model (SBSM). The spine-related information on back surface, including back surface spinous curve and local symmetry nearby spinous curve is extracted from the RGBD images of human back surface. We formulate the spine optimization constraints from spine-related feature on back surface and spine priors, then optimize the spine model by gradient descent to get the optimal personalized shape parameters and pose parameters of the Spine Priors Model (SPrM). We assess our reconstruction by scoliosis Cobb angle error, and the result is comparable to current X-ray based methods.

Biomechanical Analysis of Two-Level Novel Cage-Type Implant for Anterior Cervical Discectomy and Fusion: A Finite Element Analysis.

One of the standard treatments for spinal diseases is anterior cervical discectomy and fusion (ACDF). ACDF is a secure and successful operation that prevents patients to improve their pain and function. The mechanical goal of the ACDF is to prevent motion between adjoining vertebrae by a novel cage-screw implant. The objective of this study is to analyze the biomechanical flexibility in terms of the range of motion (ROM) of two-level ACDF fixation using the finite element method (FEM). A CT scan-based FEM model of the cervical spine (C2-C7) is used and two-level cage is implanted at C4-C6 segments. A 50-N compressive force and 1-Nm moment are applied on C2 vertebrae and C7 is fixed in all directions. The ROM at two-level fixation (C4-C5-C6) is reduced by 55 to 88% compared with intact spine during all physiological movement. The ROM slightly increase (3-9%) at the adjacent segment. The maximum von Mises stress variations are 25-65 MPa during flexion-extension, lateral bending, and axial rotations under given loading. The maximum von Mises stress found in cage and screw is below the yield stress during all physiological movement.

Biomechanical analysis of customized cage conforming to the endplate morphology in anterior cervical discectomy fusion: A finite element analysis

In anterior cervical discectomy and fusion (ACDF), an interbody fusion device is an essential implant. An unsuitable interbody fusion device can cause postoperative complications, including subsidence and nonunion. We designed a customized intervertebral fusion device to reduce postoperative complications and validated it by finite element analysis. Herein, we built a non-homogeneous model of the C3-7 cervical spine. Three implant models (customized cage, commercial cage, and bone graft cage) were constructed and placed in the C45 cervical segment after ACDF surgery. The simulated range of motion (ROM), stress at the cage-bone interface, and stress on the cage and implants were compared under different conditions. The commercial cage showed maximum stress peaks at 40.3 MPa and 43.2 MPa in the inferior endplate of C4 and superior endplate of C5 under rotational conditions, higher compared to 29.7 MPa and 26.4 MPa, respectively, in the customized cage. The ROM was not significantly different between the three cages placed after ACDF. The stresses on the commercial cage were higher compared to the other two cages under all conditions. The bone graft in the customized cage was subject to higher stress than the commercial cage under all conditions, particularly lateral bending, wherein the maximum stress was 5.5 MPa. These results showed that a customized cage that better conformed to the vertebral anatomy was promising for reducing the risk of stress shielding and the occurrence of subsidence.

Methodology for development of single cell dendritic spine (SCDS) synaptic tagging and capture model using Virtual Cell (VCell)

Single cell dendritic spine modelling methodology has been adopted to explain structural plasticity and respective change in the neuronal volume previously. However, the single cell dendrite methodology has not been employed previously to explain one of the important aspects of memory allocation i.e., Synaptic tagging and Capture (STC) hypothesis. It is difficult to relate the physical properties of STC pathways to structural changes and synaptic strength. We create a mathematical model based on earlier reported synaptic tagging networks. We built the model using Virtual Cell (VCell) software and used it to interpret experimental data and investigate the behavior and characteristics of known Synaptic tagging candidates.•We investigate processes associated with synaptic tagging candidates and compare them to the assumptions based on the STC hypothesis.•We assess the behavior of several reported synaptic tagging candidates against the requirements outlined in the synaptic tagging hypothesis.

Fabrication of Rigid Polyurethane Foam Lumbar Spine Model for Surgical Training using Indirect Additive Manufacturing

Lumbar model is an artificial bone that is commonly used in surgical training to simulate working with the human-like bone for the trainer. The common lumbar model is made of rigid polyurethane (PU) foam and is produced using casting. However, the current lumbar model is expensive and has limitations in representing the real human lumbar, especially in geometry, visuals, and haptics. Therefore, an alternative method of fabricating lumbar models made of rigid polyurethane for surgical training using indirect additive manufacturing will be investigated in this paper. The proposed indirect additive manufacturing is a combination of 3D printing and casting methods. The main process of this method is started by fabricating a mold made of polyvinyl alcohol (PVA) using fused deposition modeling (FDM) 3D printing and subsequently casting PU foam material into the 3D printed PVA mold. Accordingly, the aim of this study is to find the optimized casting process parameters, especially for injecting the material into the mold, to achieve a better quality of lumbar model. The study was conducted using a Design of Experiment (DoE) Taguchi Orthogonal Array to optimize the casting process. The geometrical measurements of middle end-plate depth, upper end-plate width, spinal canal width, spinal canal depth, and lower pedicle length show the error ranged from 0.14% to 0.85%. The average porosity, measured from the body, lamina, and spinous, was found to be non-uniform. It is ranged from 19.58% to 21.94% on the middle part and 39.78% to 45.41% on the subsurface of lumbar model. The density was increased by 64.89% compared to the reference open molded PU foam.

The Effect of Anterior-Only, Posterior-Only, and Combined Anterior Posterior Fixation for Cervical Spine Injury with Soft Tissue Injury: A Finite Element Analysis

This finite element analysis aimed to investigate the effects of surgical procedures for cervical spine injury. A three-dimensional finite element model of the cervical spine (C2-C7) was created from computed tomography. This model contained vertebrae, intervertebral discs, anterior longitudinal ligament, and posterior ligament complex. To create the cervical spine injury model, posterior ligament complex and anterior longitudinal ligament at C3-C4 were resected and the center of the intervertebral disc was resected. We created posterior-only fixation (PF), anterior-only fixation (AF), and combined anterior-posterior fixation (APF) models. A pure moment with a compressive follower load was applied, and range of motion, annular/nucleus stress, instrument stress, and facet forces were analyzed. In all motion except for flexion, range of motion of PF, AF, and APF models decreased by 80%-95%, 85%-93%, and 97%-99% compared with the intact model. C3-C4 annulus stress of PF, AF, and APF models decreased by 28%-72%, 96%-100%, and 99%-100% compared with the intact model. Facet contact forces of PF, AF, and APF models decreased by 77%-79%, 97%-99%, and 77%-86% at C3-C4 compared with the intact model. Screw stress in the PF model was higher than in the APF model, and plate stress in the AF model was lower than in the APF model, but bone graft stress in the AF model was higher than in the APF model. Cervical stabilization was preserved by the APF model. Regarding range of motion, the PF model had an advantage compared with the AF model except for flexion. An understanding of biomechanics provides useful information for the clinician.

Morphable models of the lumbar spine to vary geometry based on pathology, demographics, and anatomical measurements

The shape of the lumbar spine influences its function and dysfunction. Yet examining the influence of geometric differences associated with pathology or demographics on lumbar biomechanics is challenging in vivo where these effects cannot be isolated, and the use of simple anatomical measurements does not fully capture the complex three-dimensional geometry. The goal of this work was to develop and share morphable models of the lumbar spine that allow geometry to be varied according to pathology, demographics, or anatomical measurements. Partial least squares regression was used to generate statistical shape models that quantify geometric differences associated with pathology, demographics, and anatomical measurements from the lumbar spines of 87 patients. To determine if the morphable models detected meaningful geometric differences, the ability of the morphable models to classify spines was compared with models generated from random labels. The models for disc herniation (p < 0.04), spondylolisthesis (p < 0.001), and sex (p < 0.01) all performed significantly better than the random models. Age was predicted with a root mean square error of 14.1 years using the age-based model. The morphable models for anatomical measurements were able to produce instances with root mean square errors less than 0.8°, 0.3 cm2, and 0.7 mm between desired and resulting measurements. This method can be used to produce morphable models that enable further analysis of the relationship among shape, pathology, demographics, and function through computational simulations. The morphable models and code are available at https://github.com/aclouthier/morphable-lumbar-model.

Development of a Pedicle Screw Fixation Simulation Model for Surgical Training Using a 3-Dimensional Printer

Training surgeons in pedicle screw fixation (PSF) techniques during actual surgery is limited because of patient safety, complications, and surgical efficiency issues. Recent technical developments are leading the world to an era of personalized three-dimensional (3D) printing. This study aimed to evaluate the educational effect of using a 3D-printed spine model to train beginners in PSF techniques to improve screw accuracy and procedure time. Computed tomography (CT) scan data were used in a 3D printer to produce a life-size lumbar spine replica of L1-3 vertebrae. Four residents performed PSF thrice. Each resident performed 18 screw fixations on both sides (6 screws per trial). The time to complete the procedure and pedicle violation was recorded. The average time for the 3 procedures was 42.1±2.9 minutes, 38.8±3.3 minutes, and 32.1±2.5 minutes, respectively. Furthermore, the average pedicle screw score for the 3 procedures was 13.0±0.8, 14.5±0.6, and 16.0±0.8, respectively. As the trial was repeated, the procedure time decreased and the accuracy of screw fixation tended to be more accurate. It was possible to decrease the procedure time and increase accuracy through repeated training using the 3D-printed spine model. By implementing a 3Dprinted spine model based on the patient's actual CT data, surgeons can perform simulation surgery before the actual surgery. Therefore, this technology can be useful in educating residents to improve their surgical skills.

Head supported mass, moment of inertia, neck loads and stability: A simulation study

Occupations or activities where donning head-supported mass (HSM) is commonplace put operators at an elevated risk of chronic neck pain. Yet, there is no consensus about what features of HSM influence the relative contributions to neck loads. Therefore, we tested four hypotheses that could increase neck loads: (i) HSM increases gravitational moments; (ii) more muscle activation is required to stabilize the head with HSM; (iii) the position of the HSM centre of mass (COM) induces gravitational moments; and (iv) the added moment of inertia (MOI) from HSM increases neck loads during head repositioning tasks.We performed a sensitivity analysis on the C5-C6 compression evaluated from a 24-degree freedom cervical spine model in OpenSim for static and dynamic movement trials. For static trials, we varied the magnitude of HSM, the position of its COM, and developed a novel stability constraint for static optimization. In dynamic trials, we varied HSM and the three principle MOIs.HSM magnitude and compression were linearly related to one another for both static and dynamic trials, with amplification factors varying between 1.9 and 3.9. Similar relationships were found for the COM position, although the relationship between C5-C6 peak compression and MOI in dynamic trials was generally nonlinear.This sensitivity analysis uncovered evidence in favour of hypotheses (i), (ii) and (iii). However, the model’s prediction of C5-C6 compression was not overly sensitive to the magnitude of MOI. Therefore, the HSM mass properties may be more influential on neck compression than MOI properties, even during dynamic tasks.