Material property sensitivity analysis on resonant frequency characteristics of the human spine.

The aim of this study is to investigate the effect of material property changes in the spinal components on the resonant frequency characteristics of the human spine. Several investigations have reported the material property sensitivity of human spine under static loading conditions, but less research has been devoted to the material property sensitivity of spinal biomechanical characteristics under a vibration environment. A detailed three-dimensional finite element model of the human spine, T12-pelvis, was built and used to predict the influence of material property variation on the resonant frequencies of the human spine. The simulation results reveal that material properties of spinal components have obvious influences on the dynamic characteristics of the spine. The annulus ground substance is the dominant component affecting the vertical resonant frequencies of the spine. The percentage change of the resonant frequency relative to the basic condition was more than 20% if Young's modulus of disc annulus is less than 1.5 MPa. The vertical resonant frequency may also decrease if Poisson's ratio of nucleus pulposus of intervertebral disc decreases.

Vibration characteristics of the human spine under axial cyclic loads: effect of frequency and damping.

STUDY DESIGN: A nonlinear finite element model of lumbar spine segment L3-L5 was developed. The effects of upper body mass, nucleus injury, damping, and different vibration frequency loads were analyzed for the whole body vibration. OBJECTIVES: To analyze the influence of whole body vibration on facets of lumbar spine and to analyze the influence of nucleus injury, upper body mass, and damping on the dynamic characteristics of lumbar spine. SUMMARY OF BACKGROUND DATA: Many studies have investigated whole body vibration for lumbar spine. However, very few investigations analyzed the influence of whole body vibration on facets and vibration characteristics of the injured spine. METHODS: The nonlinear finite element model of the L3-L5 segment was constructed based on the embalmed vertebra geometry and validated. Besides static and modal analyses, transient dynamic analyses were also conducted on the model with an upper body mass under damping and different frequency cyclic loads. RESULTS: In the period of human spine vibration, the vibration effects of different regions of the lumbar spine are not the same. Anterior regions of the L3-L5 segment show small vibration amplitudes, but posterior regions show large amplitudes. The vibration amplitude of facet contact force is more than 2.0-fold as large as that of displacement and stress on vertebrae or discs. To decrease the weight of the upper body will increase the resonant frequency. To remove the nucleus will decrease the resonant frequencies. The vibration displacement, stress, and facet contact force will reduce generally by 50% using damping ratio 0.08. CONCLUSIONS: The posterior regions of intervertebral discs of the lumbar spine are easy to injure during long-term whole body vibration compared to anterior regions. The vibration of human spine is more dangerous to facets, especially during whole body vibration approximating a sympathetic vibration, which may lead to abnormal remodeling and disorder of the lumbar spine.

Effect of facetectomy on lumbar spinal stability under sagittal plane loadings.

STUDY DESIGN: A study using an anatomically accurate finite-element model of a L2-L3 motion segment to investigate the biomechanical effects of graded bilateral and unilateral facetectomies of L3 under flexion and extension loadings. OBJECTIVE: To predict the amount of facetectomy on lumbar motion segment that would cause segmental instability, therefore enhancing the understanding concerning the role of the facet under sagittal loadings. SUMMARY OF BACKGROUND DATA: This study provides a quantitative study on the role of facets in preserving segmental lumbar stability. Previous analytical models lack of three-dimensional structural characterization and insufficient element representation for facet joints. METHODS: A validated finite-element L2-L3 model was subjected to sagittal loadings at 7.5 Nm. Effects of ligaments and facets were examined to establish their relative importance on segment response. The effect of iatrogenic changes (graded unilateral and bilateral facetectomy) was then investigated under these loadings to predict the alterations in terms of gross external (angular and coupled) responses, flexibilities, and facet load. RESULTS: This study shows the importance of preserving ligaments to prevent rotational instabilities for motion segment under flexion. The effect of the facetectomy on the motion segment is insignificant under flexion. In extension, unilateral facetectomy and resection on contralateral facet markedly alters the rotational motion and flexibilities as well as coupled motions. Also, unilateral complete facetectomy with resection of less than 100% on contralateral facet generates high facet load. CONCLUSIONS: Clinically, this study suggests that it may be appropriate to incorporate additional stabilization procedure in restoring the spinal strength and stability for surgical intervention of unilateral complete facetectomy and resection on contralateral facet. The exploitation of the finite-element method to simulate clinically related situations permits an improved understanding of lumbar spinal stability to assist in defining clinical expectation for various forms of surgical intervention of the operative procedures.

Kinematics of the thoracic T10-T11 motion segment: locus of instantaneous axes of rotation in flexion and extension.

The purpose of this study was to determine the locations and loci of instantaneous axes of rotation (IARs) of the T10-T11 motion segment in flexion and extension. An anatomically accurate three-dimensional model of thoracic T10-T11 functional spinal unit (FSU) was developed and validated against published experimental data under flexion, extension, lateral bending, and axial rotation loading configurations. The validated model was exercised under six load configurations that produced motions only in the sagittal plane to characterize the loci of IARs for flexion and extension. The IARs for both flexion and extension under these six load types were directly below the geometric center of the moving vertebra, and all the loci of IARs were tracked superoanteriorly for flexion and inferoposteriorly for extension with rotation. These findings may offer an insight to better understanding of the kinematics of the human thoracic spine and provide clinically relevant information for the evaluation of spinal stability and implant device functionality.

Effects of cervical cages on load distribution of cancellous core: a finite element analysis.

METHODS: The study was designed to analyze the load distribution of the cancellous core after implantation of vertical ring cages made of titanium, cortical bone, and tantalum using the finite element (FE) method. The intact FE model of C5-C6 motion segment was validated with experimental results. RESULTS: The percentage of load distribution in cancellous core dropped by about one-third of the level for the intact model after the cage implantation. The difference among cages made of different materials (or different stiffnesses) was not very obvious. CONCLUSIONS: These results implied that the influence of the cage on the load transfer in the cancellous core is greatly related to the cage's dimensions and position within the intervertebral space. The dimension and position of the cage that least disturb the load distribution in cancellous core could be criteria in cage development.

Validation of T10-T11 finite element model and determination of instantaneous axes of rotations in three anatomical planes.

STUDY DESIGN: A finite element (FE) model of thoracic spine (T10-T11) was constructed and used to determine instantaneous axes of rotation (IARs). OBJECTIVES: To characterize the locations and loci of IARs in three anatomic planes. SUMMARY OF BACKGROUND DATA: The center of rotation is a part of a precise method of documenting the kinematics of a spinal segment for spinal stability and deformity assessments and for implant devices study. There is little information about loci of IARs in thoracic spine. METHODS: A FE model of thoracic spine (T10-T11) was developed and validated against published data. The validated model was then used to determine the locations and loci of IARs in three anatomic planes. RESULTS: Within the validated range, The IARs locations and loci were found to vary with the applied pure moments. Under flexion and extension pure moments, the loci of IARs were tracked anterosuperiorly for flexion and posteroinferiorly for extension with rotation between the superior endplate and the geometrical center of the inferior vertebra T11. Under left and right lateral bending pure moments, the loci were detected to diverge latero-inferiorly from the mid-height of the intervertebral disc, then converge medio-inferiorly toward the geometrical center of the inferior vertebra T11. For axial rotation, the IARs were located between anterior nucleus and anulus and found to diverge in opposite direction latero-posteriorly with increasing left and right axial torque. CONCLUSIONS: The results of IARs would provide further understanding to the kinematics and biomechanical responses of the human thoracic spine, which is important for the diagnosis of disc degeneration and implant study.

Finite element analysis of cervical spinal instability under physiologic loading.

The definition of cervical spinal instability has been a subject of considerable debate and has not been clearly established. Stability of the motion segment is provided by ligaments, facet joints, and disc, which restrict range of movement. Moreover, permanent damage to one of the stabilizing structures alters the roles of the other two. Although many studies have been conducted to investigate cervical injuries, to date there are only limited finite element investigations reported in the literature on the biomechanical response of the cervical spine in these respects. A comprehensive, geometric, nonlinear finite element model of the lower cervical spine has been successfully developed and validated under compression, anterior-posterior shear, and sagittal moments. Injury studies were done by varying each spinal component independently from the validated model. Seven analyses were conducted for each injury simulation (model without ligaments, model without facets, model without facets and ligaments, and model without disc nucleus). Results indicate that the role of the ligaments in resisting anterior and posterior shear and flexion and axial rotation moments is important. Under other physiologic loading (anterior-posterior shear, flexion-extension, lateral bending, and axial rotation), the disc nucleus is responsible for the initial stiffness of the cervical spine. The results also highlight the importance of facets in resisting compression at higher loads, anterior shear, extension, lateral bending, and torsion. The results provide new insight through injury simulation into the role of the various spinal components in providing cervical spinal stability. These findings seem to correlate well with experimental results as well as with common clinical experience.