Finite element biomechanical analysis of 3D printed intervertebral fusion cage in osteoporotic population.

OBJECTIVE: To study the biomechanical characteristics of each tissue structure when using different 3D printing Cage in osteoporotic patients undergoing interbody fusion. METHODS: A finite element model of the lumbar spine was reconstructed and validated with regarding a range of motion and intervertebral disc pressure from previous in vitro studies. Cage and pedicle screws were implanted and part of the lamina, spinous process, and facet joints were removed in the L4/5 segment of the validated mode to simulate interbody fusion. A 280 N follower load and 7.5 N·m moment were applied to different postoperative models and intact osteoporotic model to simulate lumbar motion. The biomechanical characteristics of different models were evaluated by calculating and analyzing the range of motion of the fixed and cephalic adjacent segment, the stress of the screw-rod system, the stress at the interface between cage and L5 endplate, and intervertebral disc pressure of the adjacent segment. RESULTS: After rigid fixation, the range of motion of the fixed segment of model A-C decreased significantly, which was much smaller than that of the osteoporotic model. And with the increase of the axial area of the interbody fusion cages, the fixed segment of model A-C tended to be more stable. The range of motion and intradiscal pressure of the spinal models with different interbody fusion cages were higher than those of the complete osteoporosis model, but there was no significant difference between the postoperative models. On the other hand, the L5 upper endplate stress and screw-rod system stress of model A-C show a decreasing trend in different directions of motion. The stress of the endplate is the highest during flexion, which can reach 40.5 MPa (model A). The difference in endplate stress between models A-C was the largest during lateral bending. The endplate stress of models A and B was 150.5% and 140.9% of that of model C, respectively. The stress of the screw-rod system was the highest during lateral bending (model A, 102.0 MPa), which was 108.4%, 102.4%, 110.4%, 114.2% of model B and 158.5%, 110.1%, 115.8%, 125.4% of model C in flexion, extension, lateral bending, and rotation, respectively. CONCLUSIONS: For people with osteoporosis, no matter what type of cage is used, good immediate stability can be achieved after surgery. Larger cage sizes provide better fixation without significantly increasing ROM and IDP in adjacent segments, which may contribute to the development of ASD. In addition, larger cage sizes can disperse endplate stress and reduce stress concentration, which is of positive significance in preventing cage subsidence after operation. The cage and screw rod system establish a stress conduction pathway on the spine, and a larger cage greatly enhances the stress-bearing capacity of the front column, which can better distribute the stress of the posterior spine structure and the stress borne by the posterior screw rod system, reduce the stress concentration phenomenon of the nail rod system, and avoid exceeding the yield strength of the material, resulting in the risk of future instrument failure.

Biomechanical evaluation of different sizes of 3D printed cage in lumbar interbody fusion-a finite element analysis.

OBJECTIVE: To study the biomechanical characteristics of various tissue structures of different sizes of 3D printed Cage in lumbar interbody fusion. METHODS: A finite element model of normal spine was reconstructed and verified. Pedicle screws and Cage of different sizes were implanted in the L4/5 segment to simulate lumbar interbody fusion. The range of motion of the fixed and cephalic adjacent segment, the stress of the screw-rod system, the stress at the interface between cage and L5 endplate, and intervertebral disc pressure of the adjacent segment were calculated and analyzed. RESULTS: The range of motion and intervertebral disc pressure of the adjacent segment of each postoperative model were larger than those of the intact model, but there was not much difference between them. The stress of cage-endplate interface was also larger than that of the intact model. However, the difference is that the stress of the endplate and the screw-rod system has a tendency to decrease with the increase of the axial area of cage. CONCLUSIONS: Cage with larger axial area in lumbar interbody fusion can reduce the stress of internal fixation system and endplate, but will not increase the range of motion and intervertebral disc pressure of adjacent segment. It has a certain effect in preventing the cage subsidence, internal fixation system failure and screw rod fracture.

The effect of lifting load on the kinematic characteristics of lumbar spinous process in vivo.

BACKGROUND: There are limited data on the in vivo natural kinematics of the lumbar spinous process. This paper intends to explore the effect of lifting load on the in vivo movement mode of the lumbar spinous process and its biomechanical changes. METHODS: Ten asymptomatic subjects between the ages of 25 and 39 underwent CT scans of the lumbar spine in the supine position, and 3D models of L3-L5 were constructed. Using a Dual Fluoroscopy Imaging System (DFIS), instantaneous orthogonal fluoroscopic images of each subject's flexion-extension, left-right bending, and left-right rotational movements were taken under different loads (0 kg, 5 kg, 10 kg). The supine CT model was matched, using computer software, to the bony contours of the images from the two orthogonal views, so that the instantaneous 3D vertebral position at each location could be quantified. A Cartesian coordinate system was ultimately constructed at the tip of the spinous process to obtain the 6DOF kinematic data of the spinous process. RESULTS: In different postural movements of the trunk, there was no significant difference in the rotation angle and translation range of the lumbar spinous process under different loads (P > 0.05). In flexion to extension motion, spinous processes mainly rotate < 4° along the medial and lateral axes and translate < 4 mm along the craniocaudal direction. In the left-right bending motion, spinous processes mainly rotate < 5° along the anterior and posterior axes, and the translation is mainly coupling < 2 mm. In the rotational motion, the spinous process is mainly coupled motion, the rotation range is less than 3°, and the translation range is less than 2 mm. The distance between spinous processes measured in the supine position was 6.66 ± 2.29 mm at L3/4 and 5.08 ± 1.57 mm at L4/5. CONCLUSION: The in vivo kinematics of the lumbar spinous process will not change significantly with increasing low load. In complex motion, the spinous process is dominated by coupling motion.

Investigation of geometric deformations of the lumbar disc during axial body rotations.

BACKGROUND: Quantitative data on in vivo vertebral disc deformations are critical for enhancing our understanding of spinal pathology and improving the design of surgical materials. This study investigated in vivo lumbar intervertebral disc deformations during axial rotations under different load-bearing conditions. METHODS: Twelve healthy subjects (7 males and 5 females) between the ages of 25 and 39 were recruited. Using a combination of a dual fluoroscopic imaging system (DFIS) and CT, the images of L3-5 segments scanned by CT were transformed into three-dimensional models, which matched the instantaneous images of the lumbar spine taken by a double fluorescent X-ray system during axial rotations to reproduce motions. Then, the kinematic data of the compression and shear deformations of the lumbar disc and the coupled bending of the vertebral body were obtained. RESULTS: Relative to the supine position, the average compression deformation caused by rotation is between + 10% and - 40%, and the shear deformation is between 17 and 50%. Under physiological weightbearing loads, different levels of lumbar discs exhibit similar deformation patterns, and the deformation patterns of left and right rotations are approximately symmetrical. The deformation patterns change significantly under a 10 kg load, with the exception of the L3-4 disc during the right rotation. CONCLUSION: The deformation of the lumbar disc was direction-specific and level-specific during axial rotations and was affected by extra weight. These data can provide new insights into the biomechanics of the lumbar spine and optimize the parameters of artificial lumbar spine devices.

Biomechanical behaviour of tension-band-reconstruction titanium plate in open-door laminoplasty: a study based on finite element analysis.

OBJECTIVE: To investigate and evaluate the biomechanical behaviour of tension-band-reconstruction (TBR) and ordinary titanium plates in open-door laminoplasty by finite element (FE) analysis. METHODS: TBR titanium plate and ordinary titanium plate were implanted into a validated finite element model of healthy adult cervical vertebrae. Among them, 5 ordinary titanium plate were used in model A, 2 TBR titanium plates and 3 ordinary titanium plates were used in model B, and 5 TBR titanium plates were used in model C. The same loading conditions was applied identically to all models. Range of motion (ROM) of the vertebral body, stress distribution of the titanium plate and intradiscal pressure (IDP) were compared in flexion, extension, lateral bending and rotation. RESULTS: The ROM of model B and C was similar in flexion and extension, and both were smaller than that of model A. The highest von Mises stress in the titanium plate appears is in model C. The IDP in C2/3 was significantly higher than that in other segments in flexion. There was no significant difference in IDP among three models in left lateral bending and left axial rotation. CONCLUSION: Application of TBR titanium plate in open-door laminoplasty can reduced ROM in flexion, extension and axial rotation of the cervical vertebrae. But the increase of stress in TBR titanium plate could lead to higher risk of adverse events such as titanium plate deformation. Moreover, compared with complete TBR titanium plate, the combination of TBR titanium plate for C3 and C7 with ordinary titanium plate for the other vertebrae largely reduce the stress of the titanium plates by ensuring stability. The proposed FE model (C2-T1) exhibits a great potential in evaluating biomechanical behaviour of TBR titanium plate for open-door laminoplasty.

Biomechanical comparison of spinal column shortening - a finite element study.

BACKGROUND: At present, research on spinal shortening is mainly focused on the safe distance of spinal shortening and the mechanism of spinal cord injury, but there is no research on the biomechanical characteristics of different shortening distances. The purpose of this study was to study the biomechanical characteristics of spine and internal fixation instruments at different shortening distances by the finite element (FE) method. METHODS: An FE model of lumbar L1-S was established and referred to the previous in vitro experiments to verify the rationality of the model by verifying the Intradiscal pressure (IDP) and the range of motion (ROM) of the motion segment. Five element models of spinal shortening were designed under the safe distance of spinal shortening, and the entire L3 vertebra and both the upper and lower intervertebral discs were resected. Model A was not shortened, while models B-E were shortened by 10%, 20%, 30% and 50% of the vertebral body, respectively. Constraining the ROM of the sacrum in all directions, a 7.5 N ·m moment and 280 N follower load were applied on the L1 vertebra to simulate the motion of the lumbar vertebrae in three planes. The ROM of the operated segments, the Von Mises stress (VMS) of the screw-rod system, the VMS of the upper endplate at the interface between the titanium cage and the L4 vertebral body, and the ROM and the IDP of the adjacent segment (L5/S) were recorded and analysed. RESULTS: All surgical models showed good stability at the operated segments (L1-5), with the greatest constraint in posterior extension (99.3-99.7%), followed by left-right bending (97.9-98.7%), and the least constraint in left-right rotation (84.9-86.3%) compared with the intact model. The VMS of the screw-rod system and the ROM and IDP of the distal adjacent segments of models A-E showed an increasing trend, in which the VMS of the screw-rod system of model E was the highest under flexion (172.5 MPa). The VMS of the endplate at the interface between the cage and L4 upper endplate of models A-E decreased gradually, and these trend were the most obvious in flexion, which were 3.03, 2.95, 2.83, 2.78, and 2.61 times that of the intact model, respectively. CONCLUSION: When performing total vertebrae resection and correcting the spinal deformity, if the corrected spine has met our needs, the distance of spinal shortening should be minimized to prevent spinal cord injury, fracture of internal fixations and adjacent segment disease (ASD).

Optimization of in-brace corrective force in adolescents with Lenke type 5 curve using finite element model.

BACKGROUND: Pelvic parameters have been taken into consideration for the evaluation of the outcomes of bracing in AIS. To discuss the stress required to correct the pelvic deformity related to Lenke5 adolescent idiopathic scoliosis (AIS) by finite element analysis, and provide a reference for the shaping of the pelvic region of the brace. METHODS: An three-dimensional (3D) corrective force on the pelvic area was defined. Computed tomography images were used to reconstruct a 3D model of Lenke5 AIS. Computer-aided engineering software Abaqus was used to implement finite element analysis. By adjusting the magnitude and position of corrective forces, coronal pelvic coronal plane rotation (PCPR) and Cobb angle (CA) of lumbar curve in the coronal plane, horizontal pelvic axial plane rotation, and apical vertebra rotation (AVR) were minimized to achieve the best effect on the spine and pelvic deformity correction. The proposed corrective conditions were divided into three groups: (1) forces applied on X-axis; (2) forces applied both in the X- and Y-axis; and (3) forces applied along the X-, Y-, and Z-axis at the same time. RESULTS: In three groups, CA correction reduced by 31.5%, 42.5%, and 59.8%, and the PCPR changed to 12°, 13°, and 1° from 6.5°, respectively. The best groups of correction forces should simultaneously locate on the sagittal, transverse, and coronal planes of the pelvis. CONCLUSIONS: For Lenke5 AIS, 3D correction forces can sufficiently reduce scoliosis and pelvic asymmetrical state. Force applied along the Z-axis is vital to correct the pelvic coronal pelvic tilt associated with Lenke5 AIS.

The effect of various weight-bearing activities on the motion of lumbar facet joints in vivo.

BACKGROUND: Lumbar facet joints (LFJs) are usually related to the pathogenesis of the spine. The purpose of this paper is to study the effects of lifting load on the motion of lower lumbar facet joints in vivo. METHODS: Ten healthy volunteers aged 25 ≤ 39 years, 5 males and 5 females, were recruited. Using a dual fluoroscopy imaging system (DFIS) combined with CT, firstly, the L3-S1 segment image scanned by CT was converted into a three-dimensional model. Then, the lumbar motion images of L3-S1 vertebrae taken by the DFIS under different loads (0 kg, 5 kg, 10 kg) and different body postures (maximum flexion and extension, maximum left and right bending, and maximum left and right torsion) were captured. Finally, in the Rhino software, the instantaneous motion state of the lumbar spine is reproduced by translation and rotation according to the anatomical structure of the lumbar spine and the previous images. With the help of computer software, a Cartesian coordinate system was placed in the center of each articular surface to measure the kinematics of the articular process and to obtain 6DOF data under different loads (0 kg, 5 kg, 10 kg) in the lumbar facet joints. RESULTS: In the flexion and extension of the trunk, weight bearing reduced the translational range in the mid-lateral direction. In the L3/4 segment, the lateral translational range of the left and right facet joints gradually decreased with increasing load, and the translational range at 0 kg was significantly greater than that at 10 kg (left side: 0 kg, 0.86° ± 0.57°, 10 kg, 0.24° ± 0.26°, p = 0.01; right side: 0 kg, 0.86° ± 0.59°, 10 kg, 0.26° ± 0.27°, p = 0.01). In the L5/S1 segment, the translation range of the LFJ at 0 kg was significantly greater than that at 10 kg (p = 0.02). Other bending and rotation movements were not found to cause differential changes in the 6DOF of the LFJ. In bending, the rotation range was the largest in the L3/4 segment (p < 0.05) and gradually decreased from top to bottom. At the same level, there were significant differences in the translation range of the left and right facets in the anterior posterior and craniocaudal directions (p < 0.05). CONCLUSION: Increasing the load has a significant impact on the coupled translational movement of lumbar facet joints. The asymmetric translational movement of the left and right facet joints may be a factor that accelerates the degeneration of facet joints.

Application of dual-trajectory screws in revision surgery for lumbar adjacent segment disease: a finite element study.

BACKGROUND: Advancements in medicine and the popularity of lumbar fusion surgery have made lumbar adjacent segment disease (ASDz) increasingly common, but there is no mature plan for guiding its surgical treatment. Therefore, in this study, four different finite element (FE) ASDz models were designed and their biomechanical characteristics were analysed to provide a theoretical basis for clinical workers to choose the most appropriate revision scheme for ASDz. METHODS: According to whether internal fixation was retained, different FE models were created to simulate ASDz revision surgery, and flexion, extension, axial rotation and lateral bending were simulated by loading. The biomechanical characteristics of the adjacent segments of the intervertebral disc and the internal fixation system and the range of motion (ROM) of the lumbar vertebrae were analysed. RESULTS: The difference in the ROM of the fixed segment between FE models that did or did not retain the original internal fixation was less than 0.1°, and the difference was not significant. However, the stress of the screw-rod system when the original internal fixation was retained and prolonged fixation was performed with dual-trajectory screws was less than that when the original internal fixation was removed and prolonged fixation was performed with a long bar. Especially in axial rotation, the difference between models A and B is the largest, and the difference in peak stress reached 30 MPa. However, for the ASDz revision surgery segment, the endplate stress between the two models was the lowest, and the intradiscal pressure (IDP) of the adjacent segment was not significantly different between different models. CONCLUSION: Although ASDz revision surgery by retaining the original internal fixation and prolonging fixation with dual-trajectory screws led to an increase in stress in the fusion segment endplate, it provides stability similar to ASDz revision surgery by removing the original internal fixation and prolonging fixation with a long bar and does not lead to a significant change in the IDP of the adjacent segment while avoiding a greater risk of rod fracture.

Kinematic Characteristics and Biomechanical Changes of Lower Lumbar Facet Joints Under Different Loads.

OBJECTIVE: To explore the kinematic biomechanical changes and symmetry in the left and right sides of the facet joints of lumbar spine segments under different functional loads. METHODS: Participants (n = 10) performing standing flexion and extension movements were scanned using computed tomography (CT) and dual fluoroscopy imagine system. Instantaneous images of the L3 -S1 vertebrae were captured, and by matching a three-dimensional CT model with contours from dual fluoroscopy images, in vivo facet joint movements were reproduced and analyzed. Translations and rotations of lumbar vertebral (L3 and L4 ) facet joints of data were compared for different loads (0, 5, 10 kg). The participants performed flexion and extension movements in different weight-bearing states, the translations and angles changes were calculated respectively. RESULTS: From standing to extension, there were no statistical differences in rotation angles for the facet joint processes of different vertebral segment levels under different weight loads (P > 0.05). Mediolateral axis and cranio-caudal translations under different weight loads were not statistically different for vertebral segment levels (P > 0.05). Anteroposterior translations for L3 (1.4 ± 0.1 mm) were greater than those for L4 (1.0 ± 0.1 mm) under the different load conditions (P = 0.04). Bilaterally, mediolateral, anteroposterior, and cranio-caudal translations of the facet joints under different weights (0, 10 kg) for each segment level (L3 and L4 ) were symmetric (P > 0.05). From flexion to standing, there were no statistical differences in rotation angles for different weights (0, 5, 10 kg) for each level (L3 and L4 ) (P > 0.05). There were no statistical differences between mediolateral, anteroposterior, and cranio-caudal translations at each segment level (L3 and L4 ) under different loads (P > 0.05). Under the condition of no weight (0 kg), L3 mediolateral translations on the left side (1.7 ± 1.6 mm) were significantly greater (P = 0.03) than those on the right side (1.6 ± 1.6 mm). Left side (1.0 ± 0.7 mm) L4 mediolateral translations were significantly smaller (P = 0.03) than those on the right side (1.1 ± 0.7 mm). There were no statistical differences between different weights for either anteroposterior and cranio-caudal translations (P > 0.05). There were no statistical differences for mediolateral, anteroposterior, and cranio-caudal translations for 10 kg (P > 0.05). CONCLUSION: Lumbar spine facet joint kinematics did not change significantly with increased loads. Anteroposterior translations for L3 were greater than those for L4 of the vertebral segments are related to the coronal facet joint surface. Changes in facet surface symmetry indicates that the biomechanical pattern between facet joints may change.