Ee-Chon Teo

Evaluation of the role of ligaments, facets and disc nucleus in lower cervical spine under compression and sagittal moments using finite element method.

Cervical spinal instability due to ligamentous injury, degenerated disc and facetectomy is a subject of great controversy. There is no analytical investigation reported on the biomechanical response of cervical spine in these respects. Parametric study on the roles of ligaments, facets, and disc nucleus of human lower cervical spine (C4-C6) was conducted for the very first time using noninvasive finite element method.A three-dimensional (3D) finite element (FE) model of the human lower cervical spine, consisted of 11,187 nodes and 7730 elements modeling the bony vertebrae, articulating facets, intervertebral disc, and associated ligaments, was developed and validated against the published data under three load configurations: axial compression; flexion; and extension. The FE model was further modified accordingly to investigate the role of disc, facets and ligaments in preserving cervical spine motion segment stability in these load configurations. The passive FE model predicted the nonlinear force displacement response of the human cervical spine, with increasing stiffness at higher loads. It also predicted that ligaments, facets or disc nucleus are crucial to maintain the cervical spine stability, in terms of sagittal rotational movement or redistribution of load. FE method of analysis is an invaluable application that can supplement experimental research in understanding the clinical biomechanics of the human cervical spine.

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 in vitro 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.

First cervical vertebra (atlas) fracture mechanism studies using finite element method.

Injury mechanisms and stress distribution patterns are important in the clinical evaluation of spinal injuries. Recognition and interpretation of the failure patterns help to determine spinal instability and consequently the choice of treatment. Although, the biomechanics responses of the atlas have received much attention, it has not been investigated using theoretical modeling. Mathematical techniques such as finite element model will provide further understanding to the injury mechanisms of the atlas, which is important for the prevention, diagnosis, and treatment of spinal injuries. In the present study, a detailed three-dimensional finite element model of the human atlas (C1) was constructed, with the geometrical data obtained using a three-dimensional digitizer. Anterior arch, superior/inferior articular processes, transverse processes, posterior arch and posterior tubercule were modeled using eight-noded brick elements. Using the material properties from literature, the 7808-finite element model was exercised under three simulated axial compressive mode of pressure loading and boundary conditions to investigate the sites of failure reported in vivo and in vitro. This report demonstrates high concentration of localized stress at the anterior and posterior archs of the atlas, which agrees well with those reported in the literature. Furthermore, under simulated hyperextension, our results agreed well with the experimental findings, which show that the groove of the posterior arch is subjected to enormous bending moment. The close agreement of the failure location provided confidence to perform further analysis and in vitro experiments. These results may be potentially used to supplement experimental research in understanding the clinical biomechanics of the atlas.

Development and Validation of A C0–C7 FE Complex for Biomechanical Study

In this study, the digitized geometrical data of the embalmed skull and vertebrae (C0-C7) of a 68-year old male cadaver were processed to develop a comprehensive, geometrically accurate, nonlinear C0-C7 FE model. The biomechanical response of human neck under physiological static loadings, near vertex drop impact and rear-end impact (whiplash) conditions were investigated and compared with published experimental results. Under static loading conditions, the predicted moment-rotation relationships of each motion segment under moments in midsagittal plane and horizontal plane agreed well with experimental data. In addition, the respective predicted head impact force history and the S-shaped kinematics responses of head-neck complex under near-vertex drop impact and rear-end conditions were close to those observed in reported experiments. Although the predicted responses of the head-neck complex under any specific condition cannot perfectly match the experimental observations, the model reasonably reflected the rotation distributions among the motion segments under static moments and basic responses of head and neck under dynamic loadings. The current model may offer potentials to effectively reflect the behavior of human cervical spine suitable for further biomechanics and traumatic studies.

Finite element modeling of a 3D coupled foot–boot model

Increasingly, musculoskeletal models of the human body are used as powerful tools to study biological structures. The lower limb, and in particular the foot, is of interest because it is the primary physical interaction between the body and the environment during locomotion. The goal of this paper is to adopt the finite element (FE) modeling and analysis approaches to create a state-of-the-art 3D coupled foot-boot model for future studies on biomechanical investigation of stress injury mechanism, foot wear design and parachute landing fall simulation. In the modeling process, the foot-ankle model with lower leg was developed based on Computed Tomography (CT) images using ScanIP, Surfacer and ANSYS. Then, the boot was represented by assembling the FE models of upper, insole, midsole and outsole built based on the FE model of the foot-ankle, and finally the coupled foot-boot model was generated by putting together the models of the lower limb and boot. In this study, the FE model of foot and ankle was validated during balance standing. There was a good agreement in the overall patterns of predicted and measured plantar pressure distribution published in literature. The coupled foot-boot model will be fully validated in the subsequent works under both static and dynamic loading conditions for further studies on injuries investigation in military and sports, foot wear design and characteristics of parachute landing impact in military.

Finite element application in implant research for treatment of lumbar degenerative disc disease

Surgical treatment for disc degeneration can be roughly grouped as fusion, disc replacement and dynamic stabilization. The clinical efficacy and biomechanical features of the implants used for disc degenerations can be evaluated through short- or long-term follow up observation, in vitro and in vivo experiments and computational simulations. Finite element models are already making an important contribution to our understanding of the spine and its components. Models are being used to reveal the biomechanical function of the spine and its behavior when healthy, diseased or damaged. They are also providing support in the design and application of spinal instrumentation. The article reviewed the most recent studies in the application of FE models that address the issue of implant research for treatment of low back pain. The published studies were grouped and reviewed thoroughly based on the function of implants investigated. The considerations of the finite element analysis in these studies were further discussed.

The biomechanics of lumbar graded facetectomy under anterior-shear load

In this paper, an anatomically accurate three-dimensional finite-element (FE) model of the human lumbar spine (L2-L3) was used to study the biomechanical effects of graded bilateral and unilateral facetectomies of L3 under anterior shear. The intact L2-L3 FE model was validated under compression, tension, and shear loading and the predicted responses matched well with experimental data. The gross external (translational and coupled) responses, flexibilities, and facet load were delineated for these iatrogenic changes. Results indicted that unilateral facetectomy of greater than 75% and bilateral facetectomy of 75% or more resection markedly alter the translational displacement and flexibilities of the motion segment. This study suggests that fixation or fusion to restore strength and stability of the lumbar spine may be required for surgical intervention of greater than 75% facetectomy.

Finite-element analysis for lumbar interbody fusion under axial loading

A parametric study was conducted to evaluate axial stiffness of the interbody fusion, compressive stress, and bulging in the endplate due to changes in the spacer position with/without fusion bone using an anatomically accurate and validated L2-L3 finite-element model exercised under physiological axial compression. The results show that the spacer plays an important role in initial stability for fusion, and high compressive force is predicted at the ventral endplate for the models with the spacer and fusion bone together. By varying the positioning of the spacer anteriorly along anteroposterior axis, no significant change in terms of axial stiffness, compressive stress, and bulging of the endplate are predicted for the implant model. The findings suggest that varying the spacer position in surgical situations does not affect the mechanical behavior of the lumbar spine after interbody fusion.

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.

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

Methods: The study was 1 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.