Zheng-cheng Zhong - Academia.edu (original) (raw)

Papers by Zheng-cheng Zhong

Research paper thumbnail of Biomechanical Comparison of Two New Stand-alone Anterior Lumbar Interbody Fusion Cages with Established Fixation Techniques-A Three-dimensional Finite Element Analysis

Journal of Orthopaedic Surgery Taiwan, 2009

Purpose: Anterior lumbar interbody fusion (ALIF) device avoids posterior fixation by enhancing st... more Purpose: Anterior lumbar interbody fusion (ALIF) device avoids posterior fixation by enhancing stability via the locking screws or threaded cylinder into the vertebral bodies. Two new stand-alone cages were compared with established fixation methods based on a finite element (FE) analysis to investigate the stabilization effect among all cage designs. Methods: A validated FE model of the L1-L5 was used to implant three types of ALIF cages (SynCage-Open, SynFix, and Stabilis) at the L3-4 level; the SynCage-Open model was supplemented with pedicle screw implant, as established fixation. The ends of the L5 vertebral bodies of each model were fixed and a 400N follower load plus a 10N-m moment were applied to the top of the L1 vertebra to produce flexion-extension, lateral bending, and axial rotation, respectively. The FE software (ANSYS 9.0) was used to calculate parameters, including range of motion (ROM), facet contact force, and stresses of the annulus and implants, for comparison wi...

Research paper thumbnail of Biomechanical Comparison between Isobar and Dynamic-Transitional Optima (DTO) Hybrid Lumbar Fixators: A Lumbosacral Finite Element and Intersegmental Motion Analysis

BioMed Research International

Biomechanical performance of longitudinal component in dynamic hybrid devices was evaluated to di... more Biomechanical performance of longitudinal component in dynamic hybrid devices was evaluated to display the load-transfer effects of Dynesys cord spacer or Isobar damper-joint dynamic stabilizer on junctional problem based on various disc degenerations. The dynamic component was adapted at the mildly degenerative L3–L4 segment, and the static component was fixed at the moderately degenerative L4–L5 segment under a displacement-controlled mode for the finite element study. Furthermore, an intersegmental motion behavior was analyzed experimentally on the synthetic model under a load-controlled mode. Isobar or DTO hybrid fixator could reduce stress/motion at transition segment, but compensation was affected at the cephalic adjacent segment more than the caudal one. Within the trade-off region (as a motion-preserving balance between the transition and adjacent segments), the stiffness-related problem was reduced mostly in flexion by a flexible Dynesys cord. In contrast, Isobar damper aff...

Research paper thumbnail of Biomechanical Analysis of the Lumbar Spine on the Facet Joint Force and Intradiscal Pressure— a Realistic Finite Element Study

The lumbar spine often leads to a high incidence of disc diseases, such as herniated disc, sciati... more The lumbar spine often leads to a high incidence of disc diseases, such as herniated disc, sciatica and low back pain. Biomechanical analyses of lumbar spine have been used to explore the related problems using in vitro studies or finite element simulations. Most of the previous FE studies have used simplified models such as a quarter of the vertebrae and discs, a half of the vertebrae, or a simple whole spine model. This study developed and validated an FE model of the lumbosacral spine with a realistic geometric shape to simulate the lumbar spine subjected to several loading conditions. The effects of symmetric postures such as left and right axial rotations on the facet joint forces at various levels of the lumbar spine were compared. In addition, we investigated the effect of postures on the intradiscal pressures in the nuclei pulposi.

Research paper thumbnail of Biomechanical evaluation of a novel pedicle screw-based interspinous spacer: A finite element analysis

Medical engineering & physics, Aug 13, 2017

Interspinous spacers have been designed to provide a minimally invasive surgical technique for pa... more Interspinous spacers have been designed to provide a minimally invasive surgical technique for patients with lumbar spinal stenosis or foraminal stenosis. A novel pedicle screw-based interspinous spacer has been developed in this study, and the aim of this finite element experiment was to investigate the biomechanical differences between the pedicle screw-based interspinous spacer (M-rod system) and the typical interspinous spacer (Coflex-F™). A validated finite element model of an intact lumbar spine was used to analyze the insertions of the Coflex-F™, titanium alloy M-rod (M-Ti), and polyetheretherketone M-rod (M-PEEK), independently. The range of motion (ROM) between each vertebrae, stiffness of the implanted level, the peak stress at the intervertebral discs, and the contact forces on spinous process were analyzed. Of all three devices, the Coflex-F™ provided the largest restrictions in extension, flexion and lateral bending. For intervertebral disc, the peak stress at the impla...

Research paper thumbnail of Biomechanical Behavior Between Prodisc-L and Physio-L in Total Disc Replacement

Journal of Bone Joint Surgery British Volume, Jun 1, 2012

The emerging of non-fusion surgery is aimed to solve the long-term complication of fusion surgery... more The emerging of non-fusion surgery is aimed to solve the long-term complication of fusion surgery that may bring the adjacent disc degeneration. Among several kinds of artificial discs developed in these years, the majority in the market is Prodisc-L (Synthes Inc.) which is designed with the purpose to restore the motions including anteroposterior translation, lateral bending, and axial rotation. These is also one artificial disc called Physio-L (Nexgen Spine) which were hyper-elastic material (Polycarbonate Polyurethanes) and is designed to restore the motions maintioned above plus axial loading. The concept of using hyper-elastic material as disc is to mimic the material properties of intervetebral discs so that this disc both absorb the axial loading and also restore the physiological range of motion. Few studies focused on the biomechanical behavior of hyper-elastic artificial discs have yet been reported. Therefore, the purpose of this study is to compare the biomechanical behavior between Prodisc-L and Physio-L. A validated three-dimensional finite element model of the L1-L5 lumbar intact spine was used in this study with ANSYS software [Fig.1]. Total disc replacement surgery, partial discectomy, total nuclectomy and removal of the anterior longitudinal ligament were performed at the L3/L4 segment of this intact model, and the Prodisc-L and Physio-L was implanted into L3/L4 segment, respectively. In addition, hyper-elastic materials adopted by Physio-L are usually categorized by their hardness into soft and hard [Fig.2]. Therefore, two kinds of Physio-L were studied. A 400 N follower load and a 10 N-m moment were applied to the intact model to obtain four physiological motions as comparison baseline. The implanted models were subjected to 400 N follower load and specific moments in accordance with the hybrid test method. For the Prodisc-L model in the surgical segment, the range of motion (ROM) varied by -26%, +17%, -0.01%, and -0.04% in flexion, extension, lateral bending, and axial rotation, respectively, as compared to intact model [Fig.3]. For the Physio-L (soft) model, ROM varied by +10%, +8%, +3%, and +19% in four physiological motions, respectively. For the physio-L (hard) model, ROM varied by +1%, +8%, +1%, and +11% in four physiological motions, respectively. For the Prodisc-L model in the adjacent segments, ROM varied by +4% ∼ +10%, -2% ∼ -5%, -1% ∼ -4%, and +1% ∼ -2% in four physiological motions, respectively. For the Physio-L (soft) model, ROM varied by 0% ∼ -5%, -2% ∼ -5%, -0% ∼ -5%, and -9% ∼ -11% in four physiological motions, respectively. For the physio-L (hard) model, ROM varied by +4% ∼ -2%, +8% ∼ -5%, +1 ∼ -5%, and +11% ∼ -6% in four physiological motions, respectively. As seemed in the simulation, the behavior of Physio-L (both soft and hard) is similar to that of intact model under flexion and extension, but not in axial rotation. In addition, Physio-L (hard) model is more similar to intact model as compared to Physio-L (soft) model.

Research paper thumbnail of Finite element analysis of the lumbar spine with a new cage using a topology optimization method

Medical Engineering & Physics, 2006

Research paper thumbnail of Preserving Posterior Complex Can Prevent Adjacent Segment Disease following Posterior Lumbar Interbody Fusion Surgeries: A Finite Element Analysis

PloS one, 2016

To investigate the biomechanical effects of the lumbar posterior complex on the adjacent segments... more To investigate the biomechanical effects of the lumbar posterior complex on the adjacent segments after posterior lumbar interbody fusion (PLIF) surgeries. A finite element model of the L1-S1 segment was modified to simulate PLIF with total laminectomy (PLIF-LAM) and PLIF with hemilaminectomy (PLIF-HEMI) procedures. The models were subjected to a 400N follower load with a 7.5-N.m moment of flexion, extension, torsion, and lateral bending. The range of motion (ROM), intradiscal pressure (IDP), and ligament force were compared. In Flexion, the ROM, IDP and ligament force of posterior longitudinal ligament, intertransverse ligament, and capsular ligament remarkably increased at the proximal adjacent segment in the PLIF-LAM model, and slightly increased in the PLIF-HEMI model. There was almost no difference for the ROM, IDP and ligament force at L5-S1 level between the two PLIF models although the ligament forces of ligamenta flava remarkably increased compared with the intact lumbar sp...

Research paper thumbnail of Finite element analysis of the spondylolysis in lumbar spine

Bio-medical materials and engineering, 2006

Spondylolysis is a fracture of the bone lamina in the pars interarticularis and has a high risk o... more Spondylolysis is a fracture of the bone lamina in the pars interarticularis and has a high risk of developing spondylolisthesis, as well as traction on the spinal cord and nerve root, leading to spinal disorders or low back pain when the lumbar spine is subjected to high external forces. Previous studies mostly investigated the mechanical changes of the endplate in spondylolysis. However, little attention has been focused on the entire structural changes that occur in spondylolysis. Therefore, the purpose of this study was to evaluate the biomechanical changes in posterior ligaments, disc, endplate, and pars interarticularis between the intact lumbar spine and spondylolysis. A total of three finite element models, namely the intact L2-L4 lumbar spine, lumbar spine with unilateral pars defect and with bilateral pars defect were established using a software ANSYS 6.0. A loading of 10 N.m in flexion, extension, left torsion, right torsion, left lateral bending, and right lateral bendin...

Research paper thumbnail of articleBiomechanical analysis of the lumbar spine on facet joint force and intradiscal pressure-a finite element study

Background: Finite element analysis results will show significant differences if the model used i... more Background: Finite element analysis results will show significant differences if the model used is performed under various material properties, geometries, loading modes or other conditions. This study adopted an FE model, taking into account the possible asymmetry inherently existing in the spine with respect to the sagittal plane, with a more geometrically realistic outline to analyze and compare the biomechanical behaviour of the lumbar spine with regard to the facet force and intradiscal pressure, which are associated with low back pain symptoms and other spinal disorders. Dealing carefully with the contact surfaces of the facet joints at various levels of the lumbar spine can potentially help us further ascertain physiological behaviour concerning the frictional effects of facet joints under separate loadings or the responses to the compressive loads in the discs. Methods: A lumbar spine model was constructed from processes including smoothing the bony outline of each scan imag...

Research paper thumbnail of Finite Element Analysis Using a Follower Load on Fusion and Non-Fusion Lumbar Spinal Implants

Purpose: The follower load tangent to the spinal curvature mimics the physiologic compressive loa... more Purpose: The follower load tangent to the spinal curvature mimics the physiologic compressive loads for in vivo tests. In vitro tests for spinal implants have advanced from the traditional load-control method (LCM) to the displacement-control method (DCM). This study used a follower load in finite element (FE) models of fusion and non-fusion spinal implants to evaluate the differences of the LCM and DCM at each motion segment. Methods: A FE model of the intact L1-L5 (INT) was validated and used to implant artificial disc replacement (ADR), or anterior lumbar interbody fusion (ALIF) cage supplement with pedicle screws fixation at the L3-4 level. The follower load was applied by using thermo-isotropic link elements passing through the instantaneous center of rotation at each motion segment with the L5 bottom constrained. The LCM imposed 10 Nm moments of flexion, extension, bending and rotation with a 400 N follower preload. The DCM applied motion that matched the angular displacement ...

Research paper thumbnail of Correction: Preserving Posterior Complex Can Prevent Adjacent Segment Disease following Posterior Lumbar Interbody Fusion Surgeries: A Finite Element Analysis

PloS one, 2017

The seventh author's name is spelled incorrectly. The correct name is: Zhe-Chen Li.

Research paper thumbnail of Using Finite Element Method to Determine Pad Positions in a Boston Brace for Enhancing Corrective Effect on Scoliotic Spine: A Preliminary Analysis

Scoliosis, the three-dimensional (3D) deformation of the spinal column, is likely to worsen witho... more Scoliosis, the three-dimensional (3D) deformation of the spinal column, is likely to worsen without early intervention. The Boston brace has been widely used to correct abnormal spinal curvature and prevent deformity progression. However, two-dimensional (2D) X-ray images and an experienced therapist are required for a proper Boston brace to be fabricated. There has been little research on how to determine the pad positions of a Boston brace to yield a better corrective effect on a scoliotic spine. This study uses the finite element (FE) method to determine the optimal pad positions for a Boston brace. The software ANSYS 11.0 is employed to establish an FE model of the trunk and brace from S-shaped scoliosis. After the FE model is validated, fifteen FE models with pad positions rotated in a range of 0 to 30 degrees and translated in a range of 0 to 30 mm are created to investigate the effects on correcting a scoliotic spine. A 3D FE model of a scoliotic spine is established to calcu...

Research paper thumbnail of Finite element analysis of the spondylolysis in lumbar spine

Bio-medical materials and engineering, 2006

Spondylolysis is a fracture of the bone lamina in the pars interarticularis and has a high risk o... more Spondylolysis is a fracture of the bone lamina in the pars interarticularis and has a high risk of developing spondylolisthesis, as well as traction on the spinal cord and nerve root, leading to spinal disorders or low back pain when the lumbar spine is subjected to high external forces. Previous studies mostly investigated the mechanical changes of the endplate in spondylolysis. However, little attention has been focused on the entire structural changes that occur in spondylolysis. Therefore, the purpose of this study was to evaluate the biomechanical changes in posterior ligaments, disc, endplate, and pars interarticularis between the intact lumbar spine and spondylolysis. A total of three finite element models, namely the intact L2-L4 lumbar spine, lumbar spine with unilateral pars defect and with bilateral pars defect were established using a software ANSYS 6.0. A loading of 10 N.m in flexion, extension, left torsion, right torsion, left lateral bending, and right lateral bendin...

Research paper thumbnail of Biomechanical Comparison of Instrumented Posterior Lumbar Interbody Fusion With One or Two Cages by Finite Element Analysis

Spine, 2006

Study Design. Using finite element models to study the biomechanics of lumbar instrumented poster... more Study Design. Using finite element models to study the biomechanics of lumbar instrumented posterior lumbar interbody fusion (PLIF) with one or two cages.

Research paper thumbnail of Biomechanical comparison between lumbar disc arthroplasty and fusion

Medical Engineering & Physics, 2009

The artificial disc is a mobile implant for degenerative disc replacement that attempts to lessen... more The artificial disc is a mobile implant for degenerative disc replacement that attempts to lessen the degeneration of the adjacent elements. However, inconsistent biomechanical results for the neighboring elements have been reported in a number of studies. The present study used finite element (FE) analysis to explore the biomechanical differences at the surgical and both adjacent levels following artificial disc replacement and interbody fusion procedures. First, a three-dimensional FE model of a five-level lumbar spine was established by the commercially available medical imaging software Amira 3.1.1, and FE software ANSYS 9.0. After validating the five-level intact (INT) model with previous in vitro studies, the L3/L4 level of the INT model was modified to either insert an artificial disc (ProDisc II; ADR) or incorporate bilateral posterior lumbar interbody fusion (PLIF) cages with a pedicle screw fixation system. All models were constrained at the bottom of the L5 vertebra and subjected to 150N preload and 10Nm moments under four physiological motions. The ADR model demonstrated higher range of motion (ROM), annulus stress, and facet contact pressure at the surgical level compared to the non-modified INT model. At both adjacent levels, ROM and annulus stress were similar to that of the INT model and varied less than 7%. In addition, the greatest displacement of posterior annulus occurred at the superior-lateral region. Conversely, the PLIF model showed less ROM, less annulus stress, and no facet contact pressure at the surgical level compared to the INT model. The adjacent levels had obviously high ROM, annulus stress, and facet contact pressure, especially at the adjacent L2/3 level. In conclusion, the artificial disc replacement revealed no adjacent-level instability. However, instability was found at the surgical level, which might accelerate degeneration at the highly stressed annulus and facet joint. In contrast to disc replacement results, the posterior interbody fusion procedure revealed possibly accelerative degeneration of the annulus and facet joint at both adjacent levels.

Research paper thumbnail of Influence of Dynesys System Screw Profile on Adjacent Segment and Screw

Journal of Spinal Disorders & Techniques, 2010

Displacement-controlled finite element analysis was used to evaluate the mechanical behavior of t... more Displacement-controlled finite element analysis was used to evaluate the mechanical behavior of the lumbar spine after insertion of the Dynesys dynamic stabilization system. This study aimed to investigate whether different depths of screw placement of Dynesys would affect load sharing of screw, range of motion (ROM), annulus stress, and facet contact force. In clinical follow-up, a high rate of screw complications and adjacent segment disease were found after using Dynesys. The pedicle screw in the Dynesys system is not so easy to implant into the standard position and causes the screw to protrude more prominently from the pedicle. Little is known about how the biomechanical effects are influenced by the Dynesys screw profile. The Dynesys was implanted in a 3-dimensional, nonlinear, finite element model of the L1 to L5 lumbar spine. Different depths of screw position were modified in this model by 5 and 10 mm out of the pedicle. The model was loaded to 150 N preload and controlled the same ROMs by 20, 15, 8, and 20 degrees in flexion, extension, torsion, and lateral bending, respectively. Resultant ROM, annulus stress, and facet contact force were analyzed at the surgical and adjacent level. Under flexion, extension, and lateral bending, the Dynesys provided sufficient stability at the surgical level, but increased the ROM at the adjacent level. Under flexion and lateral bending, the Dynesys alleviated annulus stress at the surgical level, but increased annulus stress at the adjacent level. Under extension, the Dynesys decreased facet loading at the surgical level but increased facet loading at the adjacent level. This study found that the Dynesys system was able to restore spinal stability and alleviate loading on disc and facet at the surgical level, but greater ROM, annulus stress, and facet loading were found at the adjacent level. In addition, profile of the screw placement caused only a minor influence on the ROM, annulus stress, and facet loading, but the screw stress was noticeably increased.

Research paper thumbnail of COMPARISON OF THE LOAD AND DISPLACEMENT CONTROLLED FINITE ELEMENT ANALYSES ON FUSION AND NON-FUSION SPINAL IMPLANTS

Journal of Biomechanics, 2007

ABSTRACT This study used finite element (FE) analysis with the load-controlled method (LCM) and t... more ABSTRACT This study used finite element (FE) analysis with the load-controlled method (LCM) and the displacement-controlled method (DCM) to examine motion differences at the implant level and adjacent levels between fusion and non-fusion implants. A validated three-dimensional intact (INT) L1-L5 FE model was used. At the L3-L4 level, the INT model was modified to surgery models, including the artificial disc replacement (ADR) of ProDisc II, and the anterior lumbar interbody fusion (ALIF) cage with pedicle screw fixation. The LCM imposed 10 Nm moments of four physiological motions and a 150 N preload at the top of L1. The DCM process was in accordance with the hybrid testing protocol. The average percentage changes in the range of motion (ROM) for whole non-operated levels were used to predict adjacent level effects (ALE%). At the implant level, the ALIF model showed similar stability with both control methods. The ADR model using the LCM had a higher ROM than the model using the DCM, especially in extension and torsion. At the adjacent levels, the ALIF model increased ALE% (at least 17 per cent) using the DCM compared with the LCM. The ADR model had an ALE% close to that of the INT model, using the LCM (average within 6 per cent), while the ALE% decreased when using the DCM. The study suggests that both control methods can be adopted to predict the fusion model at the implant level, and similar stabilization characteristics can be found. The LCM will emphasize the effects of the non-fusion implants. The DCM was more clinically relevant in evaluating the fusion model at the adjacent levels. In conclusion, both the LCM and the DCM should be considered in numerical simulations to obtain more realistic data in spinal implant biomechanics.

Research paper thumbnail of Effect of the cord pretension of the Dynesys dynamic stabilisation system on the biomechanics of the lumbar spine: a finite element analysis

European Spine Journal, 2011

The Dynesys dynamics stabilisation system was developed to maintain the mobility of motion segmen... more The Dynesys dynamics stabilisation system was developed to maintain the mobility of motion segment of the lumbar spine in order to reduce the incidence of negative effects at the adjacent segments. However, the magnitude of cord pretension may change the stiffness of the Dynesys system and result in a diverse clinical outcome, and the effects of Dynesys cord pretension remain unclear. Displacement-controlled finite element analysis was used to evaluate the biomechanical behaviour of the lumbar spine after insertion of Dynesys with three different cord pretensions. For the implanted level, increasing the cord pretension from 100 to 300 N resulted in an increase in flexion stiffness from 19.0 to 64.5 Nm/deg, a marked increase in facet contact force (FCF) of 35% in extension and 32% in torsion, a 40% increase of the annulus stress in torsion, and an increase in the high-stress region of the pedicle screw in flexion and lateral bending. For the adjacent levels, varying the cord pretension from 100 to 300 N only had a minor influence on range of motion (ROM), FCF, and annulus stress, with changes of 6, 12, and 9%, respectively. This study found that alteration of cord pretension affects the ROM and FCF, and annulus stress within the construct but not the adjacent segment. In addition, use of a 300 N cord pretension causes a much higher stiffness at the implanted level when compared with the intact lumbar spine.

Research paper thumbnail of The influence of different magnitudes and methods of applying preload on fusion and disc replacement constructs in the lumbar spine: a finite element analysis

Computer Methods in Biomechanics and Biomedical Engineering, 2013

In a finite element (FE) analysis of the lumbar spine, different preload application methods that... more In a finite element (FE) analysis of the lumbar spine, different preload application methods that are used in biomechanical studies may yield diverging results. To investigate how the biomechanical behaviour of a spinal implant is affected by the method of applying the preload, hybrid-controlled FE analysis was used to evaluate the biomechanical behaviour of the lumbar spine under different preload application methods. The FE models of anterior lumbar interbody fusion (ALIF) and artificial disc replacement (ADR) were tested under three different loading conditions: a 150 N pressure preload (PP) and 150 and 400 N follower loads (FLs). This study analysed the resulting range of motion (ROM), facet contact force (FCF), inlay contact pressure (ICP) and stress distribution of adjacent discs. The FE results indicated that the ROM of both surgical constructs was related to the preload application method and magnitude; differences in the ROM were within 7% for the ALIF model and 32% for the ADR model. Following the application of the FL and after increasing the FL magnitude, the FCF of the ADR model gradually increased, reaching 45% at the implanted level in torsion. The maximum ICP gradually decreased by 34.1% in torsion and 28.4% in lateral bending. This study concluded that the preload magnitude and application method affect the biomechanical behaviour of the lumbar spine. For the ADR, remarkable alteration was observed while increasing the FL magnitude, particularly in the ROM, FCF and ICP. However, for the ALIF, PP and FL methods had no remarkable alteration in terms of ROM and adjacent disc stress.

Research paper thumbnail of Biomechanical differences of Coflex-F and pedicle screw fixation combined with TLIF or ALIF – a finite element study

Computer Methods in Biomechanics and Biomedical Engineering, 2011

Lumbar interbody fusion is a common procedure for treating lower back pain related to degenerativ... more Lumbar interbody fusion is a common procedure for treating lower back pain related to degenerative disc diseases. The Coflex-F is a recently developed interspinous spacer, the makers of which claim that it can provide stabilisation similar to pedicle screw fixation. Therefore, this study compares the biomechanical behaviour of the Coflex-F device and pedicle screw fixation with transforaminal lumbar interbody fusion (TLIF) or anterior lumbar interbody fusion (ALIF) surgeries by using finite element analysis. The results show that the Coflex-F device combined with ALIF surgery can provide stability similar to the pedicle screw fixation combined with TLIF or ALIF surgery. Also, the posterior instrumentations (Coflex-F and pedicle screw fixation) combined with TLIF surgery had lower stability than when combined with ALIF surgery.

Research paper thumbnail of Biomechanical Comparison of Two New Stand-alone Anterior Lumbar Interbody Fusion Cages with Established Fixation Techniques-A Three-dimensional Finite Element Analysis

Journal of Orthopaedic Surgery Taiwan, 2009

Purpose: Anterior lumbar interbody fusion (ALIF) device avoids posterior fixation by enhancing st... more Purpose: Anterior lumbar interbody fusion (ALIF) device avoids posterior fixation by enhancing stability via the locking screws or threaded cylinder into the vertebral bodies. Two new stand-alone cages were compared with established fixation methods based on a finite element (FE) analysis to investigate the stabilization effect among all cage designs. Methods: A validated FE model of the L1-L5 was used to implant three types of ALIF cages (SynCage-Open, SynFix, and Stabilis) at the L3-4 level; the SynCage-Open model was supplemented with pedicle screw implant, as established fixation. The ends of the L5 vertebral bodies of each model were fixed and a 400N follower load plus a 10N-m moment were applied to the top of the L1 vertebra to produce flexion-extension, lateral bending, and axial rotation, respectively. The FE software (ANSYS 9.0) was used to calculate parameters, including range of motion (ROM), facet contact force, and stresses of the annulus and implants, for comparison wi...

Research paper thumbnail of Biomechanical Comparison between Isobar and Dynamic-Transitional Optima (DTO) Hybrid Lumbar Fixators: A Lumbosacral Finite Element and Intersegmental Motion Analysis

BioMed Research International

Biomechanical performance of longitudinal component in dynamic hybrid devices was evaluated to di... more Biomechanical performance of longitudinal component in dynamic hybrid devices was evaluated to display the load-transfer effects of Dynesys cord spacer or Isobar damper-joint dynamic stabilizer on junctional problem based on various disc degenerations. The dynamic component was adapted at the mildly degenerative L3–L4 segment, and the static component was fixed at the moderately degenerative L4–L5 segment under a displacement-controlled mode for the finite element study. Furthermore, an intersegmental motion behavior was analyzed experimentally on the synthetic model under a load-controlled mode. Isobar or DTO hybrid fixator could reduce stress/motion at transition segment, but compensation was affected at the cephalic adjacent segment more than the caudal one. Within the trade-off region (as a motion-preserving balance between the transition and adjacent segments), the stiffness-related problem was reduced mostly in flexion by a flexible Dynesys cord. In contrast, Isobar damper aff...

Research paper thumbnail of Biomechanical Analysis of the Lumbar Spine on the Facet Joint Force and Intradiscal Pressure— a Realistic Finite Element Study

The lumbar spine often leads to a high incidence of disc diseases, such as herniated disc, sciati... more The lumbar spine often leads to a high incidence of disc diseases, such as herniated disc, sciatica and low back pain. Biomechanical analyses of lumbar spine have been used to explore the related problems using in vitro studies or finite element simulations. Most of the previous FE studies have used simplified models such as a quarter of the vertebrae and discs, a half of the vertebrae, or a simple whole spine model. This study developed and validated an FE model of the lumbosacral spine with a realistic geometric shape to simulate the lumbar spine subjected to several loading conditions. The effects of symmetric postures such as left and right axial rotations on the facet joint forces at various levels of the lumbar spine were compared. In addition, we investigated the effect of postures on the intradiscal pressures in the nuclei pulposi.

Research paper thumbnail of Biomechanical evaluation of a novel pedicle screw-based interspinous spacer: A finite element analysis

Medical engineering & physics, Aug 13, 2017

Interspinous spacers have been designed to provide a minimally invasive surgical technique for pa... more Interspinous spacers have been designed to provide a minimally invasive surgical technique for patients with lumbar spinal stenosis or foraminal stenosis. A novel pedicle screw-based interspinous spacer has been developed in this study, and the aim of this finite element experiment was to investigate the biomechanical differences between the pedicle screw-based interspinous spacer (M-rod system) and the typical interspinous spacer (Coflex-F™). A validated finite element model of an intact lumbar spine was used to analyze the insertions of the Coflex-F™, titanium alloy M-rod (M-Ti), and polyetheretherketone M-rod (M-PEEK), independently. The range of motion (ROM) between each vertebrae, stiffness of the implanted level, the peak stress at the intervertebral discs, and the contact forces on spinous process were analyzed. Of all three devices, the Coflex-F™ provided the largest restrictions in extension, flexion and lateral bending. For intervertebral disc, the peak stress at the impla...

Research paper thumbnail of Biomechanical Behavior Between Prodisc-L and Physio-L in Total Disc Replacement

Journal of Bone Joint Surgery British Volume, Jun 1, 2012

The emerging of non-fusion surgery is aimed to solve the long-term complication of fusion surgery... more The emerging of non-fusion surgery is aimed to solve the long-term complication of fusion surgery that may bring the adjacent disc degeneration. Among several kinds of artificial discs developed in these years, the majority in the market is Prodisc-L (Synthes Inc.) which is designed with the purpose to restore the motions including anteroposterior translation, lateral bending, and axial rotation. These is also one artificial disc called Physio-L (Nexgen Spine) which were hyper-elastic material (Polycarbonate Polyurethanes) and is designed to restore the motions maintioned above plus axial loading. The concept of using hyper-elastic material as disc is to mimic the material properties of intervetebral discs so that this disc both absorb the axial loading and also restore the physiological range of motion. Few studies focused on the biomechanical behavior of hyper-elastic artificial discs have yet been reported. Therefore, the purpose of this study is to compare the biomechanical behavior between Prodisc-L and Physio-L. A validated three-dimensional finite element model of the L1-L5 lumbar intact spine was used in this study with ANSYS software [Fig.1]. Total disc replacement surgery, partial discectomy, total nuclectomy and removal of the anterior longitudinal ligament were performed at the L3/L4 segment of this intact model, and the Prodisc-L and Physio-L was implanted into L3/L4 segment, respectively. In addition, hyper-elastic materials adopted by Physio-L are usually categorized by their hardness into soft and hard [Fig.2]. Therefore, two kinds of Physio-L were studied. A 400 N follower load and a 10 N-m moment were applied to the intact model to obtain four physiological motions as comparison baseline. The implanted models were subjected to 400 N follower load and specific moments in accordance with the hybrid test method. For the Prodisc-L model in the surgical segment, the range of motion (ROM) varied by -26%, +17%, -0.01%, and -0.04% in flexion, extension, lateral bending, and axial rotation, respectively, as compared to intact model [Fig.3]. For the Physio-L (soft) model, ROM varied by +10%, +8%, +3%, and +19% in four physiological motions, respectively. For the physio-L (hard) model, ROM varied by +1%, +8%, +1%, and +11% in four physiological motions, respectively. For the Prodisc-L model in the adjacent segments, ROM varied by +4% ∼ +10%, -2% ∼ -5%, -1% ∼ -4%, and +1% ∼ -2% in four physiological motions, respectively. For the Physio-L (soft) model, ROM varied by 0% ∼ -5%, -2% ∼ -5%, -0% ∼ -5%, and -9% ∼ -11% in four physiological motions, respectively. For the physio-L (hard) model, ROM varied by +4% ∼ -2%, +8% ∼ -5%, +1 ∼ -5%, and +11% ∼ -6% in four physiological motions, respectively. As seemed in the simulation, the behavior of Physio-L (both soft and hard) is similar to that of intact model under flexion and extension, but not in axial rotation. In addition, Physio-L (hard) model is more similar to intact model as compared to Physio-L (soft) model.

Research paper thumbnail of Finite element analysis of the lumbar spine with a new cage using a topology optimization method

Medical Engineering & Physics, 2006

Research paper thumbnail of Preserving Posterior Complex Can Prevent Adjacent Segment Disease following Posterior Lumbar Interbody Fusion Surgeries: A Finite Element Analysis

PloS one, 2016

To investigate the biomechanical effects of the lumbar posterior complex on the adjacent segments... more To investigate the biomechanical effects of the lumbar posterior complex on the adjacent segments after posterior lumbar interbody fusion (PLIF) surgeries. A finite element model of the L1-S1 segment was modified to simulate PLIF with total laminectomy (PLIF-LAM) and PLIF with hemilaminectomy (PLIF-HEMI) procedures. The models were subjected to a 400N follower load with a 7.5-N.m moment of flexion, extension, torsion, and lateral bending. The range of motion (ROM), intradiscal pressure (IDP), and ligament force were compared. In Flexion, the ROM, IDP and ligament force of posterior longitudinal ligament, intertransverse ligament, and capsular ligament remarkably increased at the proximal adjacent segment in the PLIF-LAM model, and slightly increased in the PLIF-HEMI model. There was almost no difference for the ROM, IDP and ligament force at L5-S1 level between the two PLIF models although the ligament forces of ligamenta flava remarkably increased compared with the intact lumbar sp...

Research paper thumbnail of Finite element analysis of the spondylolysis in lumbar spine

Bio-medical materials and engineering, 2006

Spondylolysis is a fracture of the bone lamina in the pars interarticularis and has a high risk o... more Spondylolysis is a fracture of the bone lamina in the pars interarticularis and has a high risk of developing spondylolisthesis, as well as traction on the spinal cord and nerve root, leading to spinal disorders or low back pain when the lumbar spine is subjected to high external forces. Previous studies mostly investigated the mechanical changes of the endplate in spondylolysis. However, little attention has been focused on the entire structural changes that occur in spondylolysis. Therefore, the purpose of this study was to evaluate the biomechanical changes in posterior ligaments, disc, endplate, and pars interarticularis between the intact lumbar spine and spondylolysis. A total of three finite element models, namely the intact L2-L4 lumbar spine, lumbar spine with unilateral pars defect and with bilateral pars defect were established using a software ANSYS 6.0. A loading of 10 N.m in flexion, extension, left torsion, right torsion, left lateral bending, and right lateral bendin...

Research paper thumbnail of articleBiomechanical analysis of the lumbar spine on facet joint force and intradiscal pressure-a finite element study

Background: Finite element analysis results will show significant differences if the model used i... more Background: Finite element analysis results will show significant differences if the model used is performed under various material properties, geometries, loading modes or other conditions. This study adopted an FE model, taking into account the possible asymmetry inherently existing in the spine with respect to the sagittal plane, with a more geometrically realistic outline to analyze and compare the biomechanical behaviour of the lumbar spine with regard to the facet force and intradiscal pressure, which are associated with low back pain symptoms and other spinal disorders. Dealing carefully with the contact surfaces of the facet joints at various levels of the lumbar spine can potentially help us further ascertain physiological behaviour concerning the frictional effects of facet joints under separate loadings or the responses to the compressive loads in the discs. Methods: A lumbar spine model was constructed from processes including smoothing the bony outline of each scan imag...

Research paper thumbnail of Finite Element Analysis Using a Follower Load on Fusion and Non-Fusion Lumbar Spinal Implants

Purpose: The follower load tangent to the spinal curvature mimics the physiologic compressive loa... more Purpose: The follower load tangent to the spinal curvature mimics the physiologic compressive loads for in vivo tests. In vitro tests for spinal implants have advanced from the traditional load-control method (LCM) to the displacement-control method (DCM). This study used a follower load in finite element (FE) models of fusion and non-fusion spinal implants to evaluate the differences of the LCM and DCM at each motion segment. Methods: A FE model of the intact L1-L5 (INT) was validated and used to implant artificial disc replacement (ADR), or anterior lumbar interbody fusion (ALIF) cage supplement with pedicle screws fixation at the L3-4 level. The follower load was applied by using thermo-isotropic link elements passing through the instantaneous center of rotation at each motion segment with the L5 bottom constrained. The LCM imposed 10 Nm moments of flexion, extension, bending and rotation with a 400 N follower preload. The DCM applied motion that matched the angular displacement ...

Research paper thumbnail of Correction: Preserving Posterior Complex Can Prevent Adjacent Segment Disease following Posterior Lumbar Interbody Fusion Surgeries: A Finite Element Analysis

PloS one, 2017

The seventh author's name is spelled incorrectly. The correct name is: Zhe-Chen Li.

Research paper thumbnail of Using Finite Element Method to Determine Pad Positions in a Boston Brace for Enhancing Corrective Effect on Scoliotic Spine: A Preliminary Analysis

Scoliosis, the three-dimensional (3D) deformation of the spinal column, is likely to worsen witho... more Scoliosis, the three-dimensional (3D) deformation of the spinal column, is likely to worsen without early intervention. The Boston brace has been widely used to correct abnormal spinal curvature and prevent deformity progression. However, two-dimensional (2D) X-ray images and an experienced therapist are required for a proper Boston brace to be fabricated. There has been little research on how to determine the pad positions of a Boston brace to yield a better corrective effect on a scoliotic spine. This study uses the finite element (FE) method to determine the optimal pad positions for a Boston brace. The software ANSYS 11.0 is employed to establish an FE model of the trunk and brace from S-shaped scoliosis. After the FE model is validated, fifteen FE models with pad positions rotated in a range of 0 to 30 degrees and translated in a range of 0 to 30 mm are created to investigate the effects on correcting a scoliotic spine. A 3D FE model of a scoliotic spine is established to calcu...

Research paper thumbnail of Finite element analysis of the spondylolysis in lumbar spine

Bio-medical materials and engineering, 2006

Spondylolysis is a fracture of the bone lamina in the pars interarticularis and has a high risk o... more Spondylolysis is a fracture of the bone lamina in the pars interarticularis and has a high risk of developing spondylolisthesis, as well as traction on the spinal cord and nerve root, leading to spinal disorders or low back pain when the lumbar spine is subjected to high external forces. Previous studies mostly investigated the mechanical changes of the endplate in spondylolysis. However, little attention has been focused on the entire structural changes that occur in spondylolysis. Therefore, the purpose of this study was to evaluate the biomechanical changes in posterior ligaments, disc, endplate, and pars interarticularis between the intact lumbar spine and spondylolysis. A total of three finite element models, namely the intact L2-L4 lumbar spine, lumbar spine with unilateral pars defect and with bilateral pars defect were established using a software ANSYS 6.0. A loading of 10 N.m in flexion, extension, left torsion, right torsion, left lateral bending, and right lateral bendin...

Research paper thumbnail of Biomechanical Comparison of Instrumented Posterior Lumbar Interbody Fusion With One or Two Cages by Finite Element Analysis

Spine, 2006

Study Design. Using finite element models to study the biomechanics of lumbar instrumented poster... more Study Design. Using finite element models to study the biomechanics of lumbar instrumented posterior lumbar interbody fusion (PLIF) with one or two cages.

Research paper thumbnail of Biomechanical comparison between lumbar disc arthroplasty and fusion

Medical Engineering & Physics, 2009

The artificial disc is a mobile implant for degenerative disc replacement that attempts to lessen... more The artificial disc is a mobile implant for degenerative disc replacement that attempts to lessen the degeneration of the adjacent elements. However, inconsistent biomechanical results for the neighboring elements have been reported in a number of studies. The present study used finite element (FE) analysis to explore the biomechanical differences at the surgical and both adjacent levels following artificial disc replacement and interbody fusion procedures. First, a three-dimensional FE model of a five-level lumbar spine was established by the commercially available medical imaging software Amira 3.1.1, and FE software ANSYS 9.0. After validating the five-level intact (INT) model with previous in vitro studies, the L3/L4 level of the INT model was modified to either insert an artificial disc (ProDisc II; ADR) or incorporate bilateral posterior lumbar interbody fusion (PLIF) cages with a pedicle screw fixation system. All models were constrained at the bottom of the L5 vertebra and subjected to 150N preload and 10Nm moments under four physiological motions. The ADR model demonstrated higher range of motion (ROM), annulus stress, and facet contact pressure at the surgical level compared to the non-modified INT model. At both adjacent levels, ROM and annulus stress were similar to that of the INT model and varied less than 7%. In addition, the greatest displacement of posterior annulus occurred at the superior-lateral region. Conversely, the PLIF model showed less ROM, less annulus stress, and no facet contact pressure at the surgical level compared to the INT model. The adjacent levels had obviously high ROM, annulus stress, and facet contact pressure, especially at the adjacent L2/3 level. In conclusion, the artificial disc replacement revealed no adjacent-level instability. However, instability was found at the surgical level, which might accelerate degeneration at the highly stressed annulus and facet joint. In contrast to disc replacement results, the posterior interbody fusion procedure revealed possibly accelerative degeneration of the annulus and facet joint at both adjacent levels.

Research paper thumbnail of Influence of Dynesys System Screw Profile on Adjacent Segment and Screw

Journal of Spinal Disorders & Techniques, 2010

Displacement-controlled finite element analysis was used to evaluate the mechanical behavior of t... more Displacement-controlled finite element analysis was used to evaluate the mechanical behavior of the lumbar spine after insertion of the Dynesys dynamic stabilization system. This study aimed to investigate whether different depths of screw placement of Dynesys would affect load sharing of screw, range of motion (ROM), annulus stress, and facet contact force. In clinical follow-up, a high rate of screw complications and adjacent segment disease were found after using Dynesys. The pedicle screw in the Dynesys system is not so easy to implant into the standard position and causes the screw to protrude more prominently from the pedicle. Little is known about how the biomechanical effects are influenced by the Dynesys screw profile. The Dynesys was implanted in a 3-dimensional, nonlinear, finite element model of the L1 to L5 lumbar spine. Different depths of screw position were modified in this model by 5 and 10 mm out of the pedicle. The model was loaded to 150 N preload and controlled the same ROMs by 20, 15, 8, and 20 degrees in flexion, extension, torsion, and lateral bending, respectively. Resultant ROM, annulus stress, and facet contact force were analyzed at the surgical and adjacent level. Under flexion, extension, and lateral bending, the Dynesys provided sufficient stability at the surgical level, but increased the ROM at the adjacent level. Under flexion and lateral bending, the Dynesys alleviated annulus stress at the surgical level, but increased annulus stress at the adjacent level. Under extension, the Dynesys decreased facet loading at the surgical level but increased facet loading at the adjacent level. This study found that the Dynesys system was able to restore spinal stability and alleviate loading on disc and facet at the surgical level, but greater ROM, annulus stress, and facet loading were found at the adjacent level. In addition, profile of the screw placement caused only a minor influence on the ROM, annulus stress, and facet loading, but the screw stress was noticeably increased.

Research paper thumbnail of COMPARISON OF THE LOAD AND DISPLACEMENT CONTROLLED FINITE ELEMENT ANALYSES ON FUSION AND NON-FUSION SPINAL IMPLANTS

Journal of Biomechanics, 2007

ABSTRACT This study used finite element (FE) analysis with the load-controlled method (LCM) and t... more ABSTRACT This study used finite element (FE) analysis with the load-controlled method (LCM) and the displacement-controlled method (DCM) to examine motion differences at the implant level and adjacent levels between fusion and non-fusion implants. A validated three-dimensional intact (INT) L1-L5 FE model was used. At the L3-L4 level, the INT model was modified to surgery models, including the artificial disc replacement (ADR) of ProDisc II, and the anterior lumbar interbody fusion (ALIF) cage with pedicle screw fixation. The LCM imposed 10 Nm moments of four physiological motions and a 150 N preload at the top of L1. The DCM process was in accordance with the hybrid testing protocol. The average percentage changes in the range of motion (ROM) for whole non-operated levels were used to predict adjacent level effects (ALE%). At the implant level, the ALIF model showed similar stability with both control methods. The ADR model using the LCM had a higher ROM than the model using the DCM, especially in extension and torsion. At the adjacent levels, the ALIF model increased ALE% (at least 17 per cent) using the DCM compared with the LCM. The ADR model had an ALE% close to that of the INT model, using the LCM (average within 6 per cent), while the ALE% decreased when using the DCM. The study suggests that both control methods can be adopted to predict the fusion model at the implant level, and similar stabilization characteristics can be found. The LCM will emphasize the effects of the non-fusion implants. The DCM was more clinically relevant in evaluating the fusion model at the adjacent levels. In conclusion, both the LCM and the DCM should be considered in numerical simulations to obtain more realistic data in spinal implant biomechanics.

Research paper thumbnail of Effect of the cord pretension of the Dynesys dynamic stabilisation system on the biomechanics of the lumbar spine: a finite element analysis

European Spine Journal, 2011

The Dynesys dynamics stabilisation system was developed to maintain the mobility of motion segmen... more The Dynesys dynamics stabilisation system was developed to maintain the mobility of motion segment of the lumbar spine in order to reduce the incidence of negative effects at the adjacent segments. However, the magnitude of cord pretension may change the stiffness of the Dynesys system and result in a diverse clinical outcome, and the effects of Dynesys cord pretension remain unclear. Displacement-controlled finite element analysis was used to evaluate the biomechanical behaviour of the lumbar spine after insertion of Dynesys with three different cord pretensions. For the implanted level, increasing the cord pretension from 100 to 300 N resulted in an increase in flexion stiffness from 19.0 to 64.5 Nm/deg, a marked increase in facet contact force (FCF) of 35% in extension and 32% in torsion, a 40% increase of the annulus stress in torsion, and an increase in the high-stress region of the pedicle screw in flexion and lateral bending. For the adjacent levels, varying the cord pretension from 100 to 300 N only had a minor influence on range of motion (ROM), FCF, and annulus stress, with changes of 6, 12, and 9%, respectively. This study found that alteration of cord pretension affects the ROM and FCF, and annulus stress within the construct but not the adjacent segment. In addition, use of a 300 N cord pretension causes a much higher stiffness at the implanted level when compared with the intact lumbar spine.

Research paper thumbnail of The influence of different magnitudes and methods of applying preload on fusion and disc replacement constructs in the lumbar spine: a finite element analysis

Computer Methods in Biomechanics and Biomedical Engineering, 2013

In a finite element (FE) analysis of the lumbar spine, different preload application methods that... more In a finite element (FE) analysis of the lumbar spine, different preload application methods that are used in biomechanical studies may yield diverging results. To investigate how the biomechanical behaviour of a spinal implant is affected by the method of applying the preload, hybrid-controlled FE analysis was used to evaluate the biomechanical behaviour of the lumbar spine under different preload application methods. The FE models of anterior lumbar interbody fusion (ALIF) and artificial disc replacement (ADR) were tested under three different loading conditions: a 150 N pressure preload (PP) and 150 and 400 N follower loads (FLs). This study analysed the resulting range of motion (ROM), facet contact force (FCF), inlay contact pressure (ICP) and stress distribution of adjacent discs. The FE results indicated that the ROM of both surgical constructs was related to the preload application method and magnitude; differences in the ROM were within 7% for the ALIF model and 32% for the ADR model. Following the application of the FL and after increasing the FL magnitude, the FCF of the ADR model gradually increased, reaching 45% at the implanted level in torsion. The maximum ICP gradually decreased by 34.1% in torsion and 28.4% in lateral bending. This study concluded that the preload magnitude and application method affect the biomechanical behaviour of the lumbar spine. For the ADR, remarkable alteration was observed while increasing the FL magnitude, particularly in the ROM, FCF and ICP. However, for the ALIF, PP and FL methods had no remarkable alteration in terms of ROM and adjacent disc stress.

Research paper thumbnail of Biomechanical differences of Coflex-F and pedicle screw fixation combined with TLIF or ALIF – a finite element study

Computer Methods in Biomechanics and Biomedical Engineering, 2011

Lumbar interbody fusion is a common procedure for treating lower back pain related to degenerativ... more Lumbar interbody fusion is a common procedure for treating lower back pain related to degenerative disc diseases. The Coflex-F is a recently developed interspinous spacer, the makers of which claim that it can provide stabilisation similar to pedicle screw fixation. Therefore, this study compares the biomechanical behaviour of the Coflex-F device and pedicle screw fixation with transforaminal lumbar interbody fusion (TLIF) or anterior lumbar interbody fusion (ALIF) surgeries by using finite element analysis. The results show that the Coflex-F device combined with ALIF surgery can provide stability similar to the pedicle screw fixation combined with TLIF or ALIF surgery. Also, the posterior instrumentations (Coflex-F and pedicle screw fixation) combined with TLIF surgery had lower stability than when combined with ALIF surgery.