Skip to main content

METHODS article

Front. Surg., 25 January 2021
Sec. Orthopedic Surgery

Development of a Computer-Aided Design and Finite Element Analysis Combined Method for Affordable Spine Surgical Navigation With 3D-Printed Customized Template

  • 1National Center for Spinal Disorders, Buda Health Center, Budapest, Hungary
  • 2In Silico Biomechanics Laboratory, National Center for Spinal Disorders, Buda Health Center, Budapest, Hungary
  • 3School of Ph.D. Studies, Semmelweis University, Budapest, Hungary
  • 4Do3D Innovations Ltd., Budapest, Hungary
  • 5Department of Mechanical Engineering, INSIGNEO Institute for In Silico Medicine, The University of Sheffield, Sheffield, United Kingdom
  • 6Department of Spinal Surgery, Semmelweis University, Budapest, Hungary

Introduction: Revision surgery of a previous lumbosacral non-union is highly challenging, especially in case of complications, such as a broken screw at the first sacral level (S1). Here, we propose the implementation of a new method based on the CT scan of a clinical case using 3D reconstruction, combined with finite element analysis (FEA), computer-assisted design (CAD), and 3D-printing technology to provide accurate surgical navigation to aid the surgeon in performing the optimal surgical technique by inserting a pedicle screw at the S1 level.

Materials and Methods: A step-by-step approach was developed and performed as follows: (1) Quantitative CT based patient-specific FE model of the sacrum was created. (2) The CAD model of the pedicle screw was inserted into the sacrum model in a bicortical convergent and a monocortical divergent position, by overcoming the geometrical difficulty caused by the broken screw. (3) Static FEAs (Abaqus, Dassault Systemes) were performed using 500 N tensile load applied to the screw head. (4) A template with two screw guiding structures for the sacrum was designed and manufactured using CAD design and 3D-printing technologies, and investment casting. (5) The proposed surgical technique was performed on the patient-specific physical model created with the FDM printing technology. The patient-specific model was CT scanned and a comparison with the virtual plan was performed to evaluate the template accuracy

Results: FEA results proved that the modified bicortical convergent insertion is stiffer (6,617.23 N/mm) compared to monocortical divergent placement (2,989.07 N/mm). The final template was created via investment casting from cobalt-chrome. The template design concept was shown to be accurate (grade A, Gertzbein-Robbins scale) based on the comparison of the simulated surgery using the patient-specific physical model and the 3D virtual surgical plan.

Conclusion: Compared to the conventional surgical navigation techniques, the presented method allows the consideration of the patient-specific biomechanical parameters; is more affordable, and the intraoperative X-ray exposure can be reduced. This new patient- and condition-specific approach may be widely used in revision spine surgeries or in challenging primary cases after its further clinical validation.


Spinal fixation is a routine procedure for the treatment of unstable spine due to trauma, congenital malformations, degenerative diseases, and tumors (1). The accurate placement of screws in the spine is challenging, given the risk of damage to neighboring anatomical structures (spinal cord, nerve roots, arteries, and veins) (2, 3). Computer-assisted surgery (CAS) has been adopted as a safe and accurate guiding system for the placement of pedicle and lateral mass screws in the spine (4). CAS navigation systems use optical tracking via infrared cameras incorporating 3D geometries from pre-operative CT scans or in combination with fluoroscopy-based imaging (5, 6) or intraoperative CT scans (7). Optimal registration of the spine geometry to the navigational instruments is crucial for precise screw insertion. During surgery, it is often required to perform intraoperative CT scans or use fluoroscopy to re-register the system (57). Surgical manipulation after obtaining the intraoperative CT or fluoroscopy images may cause CAS registration errors, which can result in screw malposition. This phenomenon cannot be completely excluded even with a state of the art intraoperative CT technology (7). First concept of individual templates was first introduced by Radermacher et al. (8) in the early 90's by using computer controlled milling device for the manufacturing process. Currently, the 3D-printed patient-specific surgical navigation templates are accurate (9, 10), decrease surgical time, reduce intraoperative X-ray exposure (11), and can be more accessible compared to traditional CT or fluoroscopy-based systems (12, 13). The decline in the costs of 3D-printing technology is expected to continue due to its continuous and fast development (1416). The MySpine (Medacta International SA, Castel San Pietro, CH) patient- matched pedicle targeting guide for pedicle screw placement (17) is an already clinically available device for the large international spine surgical community. However, in less developed areas of the world, where complex spinal deformity is relatively common and advanced CAS technology is not available (11, 18) 3D-printed templates are still not as widely implemented in the clinical practice, as it would be desirable.

The revision surgery of a lumbosacral non-union can be complicated by an implant related failure, with a broken pedicle screw. In the S1 segment, the convergent bicortical screw trajectory provides superior anchoring compared to any other directions, but the proper insertion of the new screws in a revision surgery due to the broken screw is extremely difficult without surgical navigation. Here, we present a complex clinical case in which the accurate surgical technique required the development of a computer-aided design (CAD) and finite-element analysis (FEA) combined method for affordable spine surgical navigation with a 3D-printed customized navigation template.


Clinical Case

The study was approved by the National Ethics Committee of Hungary and the National Institute of Pharmacy and Nutrition (reference number: OGYÉI/163-4/2019). Informed consent was obtained from the patient. A 38-years-old patient underwent multiple spine surgeries at the L5–S1 level over a 5-years period with transforaminal interbody fusion (TLIF). During the latest surgery, implant removal and S1 left side nerve root decompression were performed. Six months later, the patient was referred to our institution due to manifestation of mechanical low back pain, with no sign of sensorimotor deficit. Medical imaging at admission (Figure 1) demonstrated a broken S1 left side pedicle screw deep in the sacral bone, and a non-union in the L5–S1 intervertebral space. A revision surgery aiming at the re-fusion of the LV/SI segment was decided.


Figure 1. Clinical case of a 38-years-old male patient suffering from low back pain. The patient previously underwent multiple surgeries at the L5–S1 level. A broken left sacral screw can be identified on the standing X-ray images of the patient (A. Coronal, B. Sagittal plane). Signs of non-union are identifiable in the intervertebral space on the CT scan images of the L5 vertebra and the sacrum (C. Coronal, D. Sagittal plane).

Patient-Specific 3D Geometry Definition

For the study Quantitative Computed Tomography (QCT) scans were used, performed with a Hitachi Presto CT machine (Hitachi Presto, Hitachi Medical Corporation, Tokyo, Japan) using an in-line calibration phantom with five cylindrical insertions of known mean equivalent bone mineral density (BMD) values (0, 0.5, 0.1, 0.15, and 0.2 g/cm3) with an intensity of 225 mA and voltage of 120 kV. The imaging protocol was previously defined in the MySpine project (ICT-2009.5.3 VPH, Project ID: 269909) (19, 20), and the images were reconstructed with a voxel size of 0.6 × 0.6 × 0.6 mm3. The data were extracted from the hospital PACS in DICOM file format. To comply with the ethical approval of the patient data protection, deidentification of the DICOM data was performed using the freely available Clinical Trial Processor software (Radiological Society of North America, (21). The thresholding algorithm and manual segmentation tools (erase, paint, fill, etc.) in Mimics image analysis software (Mimics Research, Mimics Innovation Suite v21.0, Materialize, Leuven, Belgium) were used (Figure 2) to define the geometry of the sacrum and the broken screw.


Figure 2. Patient-specific geometry and FE model definition. QCT based segmentation was used to define the sacrum geometry. Hounsfield Unit (HU) values of the QCT images were converted into bone mineral density (BMD) equivalent values. Elastic properties of the sacral bone were estimated using a set of density to elasticity relationships from the literature to convert the BMD equivalent value at each element of the FE mesh to Apparent Density (ρapp) (22, 23) and then to the Elastic Modulus (E).

The resulting masks (group of voxels) were homogenously filled by preserving the outer contour of the geometrical border in 2D. From the mask, a triangulated surface mesh was automatically generated. On the 3D geometries surface smoothing (iteration: 6, smooth factor: 0.7, with shrinkage compensation) and uniform remeshing was applied (target triangle edge length 0.6 mm, sharp edge preservation, sharp edge angle 60°).

Surgical Planning and FE Model Generation

A CD Horizon Legacy (Medtronic) polyaxial pedicle screw, 45 mm long and 6.5 mm in diameter, was scanned with the ScanBox 3D scanner (Smart Optics Sensortechnik GmbH, Bochum, Germany). The model of the screw was reconstructed and modified (from polyaxial to monoaxial head) in 3-matic (Mimics Research, Mimics Innovation Suite v21.0, Materialize, Leuven, Belgium) software. The triangulated surface mesh of the screw model was uniformly re-meshed (target triangle edge length: 0.6 mm, sharp edge preservation, sharp edge angle: 60°) (Figure 3A). The screw model was virtually inserted into the 3D model of the patient's sacrum in two positions (convergent: S1, divergent: ALA), using the Mimics software's STL import tool (Figure 3B) with the consideration of the broken screw. Two non-manifold assemblies were created in the Mimics software containing the broken screw, implanted screw, and sacrum for the convergent (S1) and divergent (ALA) positions.


Figure 3. Virtual pedicle screw insertion into the patient-specific sacrum model. (A) Modified (monoaxial) virtual model of the pedicle screw. (B) Pedicle screw insertion in the convergent position (S1) and divergent position (ALA), the geometrical difficulty caused by the broken screw was overcome in both insertions. (C) Boundary condition of the FEA, the sacrum was fixed on the S1 endplate and the caudal 1/3 of the sacrum, 500 N tensile load was applied on the screw head.

The assembly was exported to the 3-matic software where nine FE meshes were generated for each of the implantation scenarios (S1, ALA). The broken screw, inserted implant, and sacrum-implant interface had a triangle set with an edge length of 0.6 mm. The outer surface of the sacrum mesh was changed in the nine models by defining the uniform triangle mesh edge length as 2.0, 2.5, 3.0, 3.5, 4.0, 4.5, 5.0, 5.5, and 6.0 mm. Adaptive meshing protocol was used for the volume mesh creation with 10-node tetrahedral elements. The maximum edge length of the meshing process corresponded with the initial edge length of the sacrum surface mesh (Supplementary Figure 1), for the screw and the broken screw the same FE mesh parameters was used in all models.

The material property assignment for the volumetric elements representing the sacral bone tissue was performed in two steps (Figure 2): first, conversion of the HU (Hounsfield Unit) values to BMD values based on the in-line phantom was performed, the conversion curve was assumed to be linear according to studies (22, 24). The obtained relationship between the HU and the apparent bone density for each element was ρapp = −0.0829 + 0.0026 HUapp [g/cm3]). Then, the bone tissue was assumed to be isotropic and linearly elastic with a Poisson's ratio of 0.3 (25). Conversion curves between the density and the elastic modulus of the bone were based on the correlation established by Kopperdahl et al. (23), E = −34.7 + 3,230·ρapp, (bone elastic modulus = E [MPa]). The FE models were exported to the Abaqus/CAEv11 (Dassault Systemes, Simulia Corp, Providence, RI, USA). For the broken and the inserted pedicle screws the material properties were defined as follows: Poisson's ratio of 0.3 (26), elastic modulus of 114,000 MPa (26). Between the screws and the sacrum tie connections were used. The finite element model was subjected to a static of 500 N tensile load applied to the screw head and it was fixed at the S1 endplate and lower third of the sacrum (Figure 3C).

Navigation Template Design, Manufacturing, and Accuracy Evaluation

The template design was based on the axis of the virtually inserted screw, individual geometry, and surface of the cranial/dorsal part of the sacrum. In the 3-matic software the two axes and surface for the template/sacrum contact were defined based on the STL assembly (broken screw, inserted implant, sacrum). The contact surface and the axes were exported to the Autodesk Fusion 360 (Autodesk Inc., California, U.S.A.) CAD software which was used for the finalization of the design (Figure 4A). The virtual model of the template was printed with masked stereolithography (MSLA) technology based 3D-printing machine (VOXEL L 3D-Printer; Parameters: building size: 125 × 65 × 65 mm, layer thickness: 0.05 mm; Material: Voxeltek Cast Resin; Do3D, Hungary) (Figure 4B). The used photopolymer resin can be used as a pattern for investment casting. Finally, the model was produced in a dental laboratory via investment casting (Hexacast induction centrifugal casting machine; Parameters: start torque: 0–21 Nm, maximum melting mass: 100 g, max heating: 1,750°C, dimensions (width × height × depth): 660 × 390 × 645 mm; Material: CoCr; PiDental, Hungary) from cobalt-chrome (Figures 4C.I,II). The accuracy of the casted part was tested via 3D scanning ScanBox 3D scanner (Smart Optics Sensortechnik GmbH, Bochum, Germany) and compared to the 3D-printed model. The point clouds resulting from the scanning were aligned and compared in the 3-matic software with the part comparison module (Figure 5).


Figure 4. Design, manufacturing, and accuracy evaluation of the navigation template. (A) Template's virtual model created via CAD software. (B) 3D-printed (MSLA technology) template (red) fits exactly on the 3D-printed (FDM technology) patient-specific physical model. (C.I,II) Final navigation template created via investment casting from cobalt-chrome (C.I ventral surface polished, C.II dorsal surface). Evaluation of the drilling accuracy was performed on the physical model in the (D) convergent position (S1) and (E) divergent position (ALA).


Figure 5. 3D scanning based geometrical accuracy measurement. Cobalt-chrome investment casted navigation template's geometrical accuracy compared to the 3D-printed navigation template model created with MSLA technology. The color map (Scale; min = −1 mm, max = 1 mm) shows the geometrical difference, projected on the 3D-printed navigation template triangle based mesh model vertices (A ventral view, B dorsal view).

The accuracy of the template was tested on a patient-specific sacrum physical model, 3D-printed with a Fused Deposition Modeling (FDM) printer (Dimension 1200es 3D-Printer; Parameters: building size: 254 × 254 × 305 mm, layer thickness: 0.330–0.254 mm; Material: ABSplus/ivory; Stratasys, Israel). The drill template was placed on the FDM sacrum model; then, a cylinder inlet was connected to the template to support the drill bit, and the drilling of the model was performed according to the S1 and ALA positions (Figures 4D,E).

The template was removed and two CT scans were performed of the sacrum model with drill bits inserted in the S1 and ALA positions. The CT scan images were imported into the Mimics software where the segmentation (thresholding) and 3D reconstruction of the patient-specific FDM sacrum model geometry and drill bits were performed. The models were registered to the initial sacrum geometry derived from the QCT via point based rigid registration by selecting anatomical landmarks in the caudal part of the sacrum (Supplementary Figure 2). This step was followed by an automatic global registration inside the 3-matic software. The registration accuracy was measured with the part comparison module of the 3-matic software (Supplementary Figure 3). The centerline for the drill bit 3D geometry was defined and an analytical primitive (cylinder, 2.5 mm in diameter) was fitted to define the drilling axis. In 3-matic software 3D angle measurement tool, line to line module (World Coordinate System, XYZ coordinates) was used to quantify the accuracy of the screw insertion by defining the angels in the 3D space between the virtual screw centerline and the drill bit centerline.


Navigation Template Geometrical Accuracy and Performance

The investment casted cobalt-chrome drill template retains the geometrical properties of the pattern (3D-printed drill template model created with MSLA technology) based on the 3D scanning evaluation Figure 6. To evaluate the drill template's performance, we used a 3D-printed patient-specific physical model. The physical model with the two drilling positions was scanned with CT, segmented and aligned to the virtual surgical plan (Supplementary Figures 2, 3). The drill template allowed a highly accurate screw insertion in both investigated positions (Figure 6). The cylinders representing the drilling axes were not perfectly colinear and coincident with the screws in the virtual surgical plan, the 3D angle between the screw centerline and the drill bit centerline for the S1 was α = 4.42° and for the ALA was α = 2.4°.


Figure 6. Visualization of the navigation template compared to the virtual plan. The red cylinders represent the drill bits' axes in the (A) convergent position (S1) and (B) divergent position (ALA), based on the evaluation performed on the patient-specific physical model. The broken screw and implanted screw geometry are part of the virtual surgical plan based on the patient's QCT.

FEA Results

Nine models were created for each screw insertion scenario (N = 9, S1 and N = 9, ALA) with increasing element numbers based on the virtual surgical plan. The FE simulation results converged above 2*105 elements for both screw insertion scenarios at ~5 min solve times on two cores. The solve time at two cores for the S1 orientation was higher compared to the ALA (Figure 7A). The convergent bicortical screw insertion (S1) provided a stiffer (6,617.23 ± 1,106.24 N/mm) situation based on the nine FE model compared to the monocortical divergent screw position FE model values (2,989.07 ± 240.24 N/mm) (Figure 7B).


Figure 7. FE simulation results. (A) Convergence analysis for the average U, displacement magnitude (nodes of the middle 1/3 of the screw head) in convergent (US1) and divergent (UALA) screw positions at different mesh element numbers. Solve time distribution (right) at different mesh element numbers [convergent (TS1) and divergent (TALA) screw positions]. (B) The convergent screw insertion (S1) is stiffer compared to the divergent (ALA) insertion.


Comparative studies have been published in recent years (2729), demonstrating the reliably, efficacy, and advantages of 3D printed navigational templates compared to other navigational methods or free-hand technique. In this study, we present a technology development process in order to create a patient-specific drill template in a complex clinical case, in which a broken screw causes geometrical difficulty for new screw insertion. In order to safely insert the new screw, without compromising the local bone structure we developed a virtual surgical plan based on the QCT of the patient. This allowed us to test two different screw positions in the model and to design a drill template for safe screw insertion at the level of the first sacral vertebra with a geometrical difficulty caused by a broken screw from a previous surgery. The present study demonstrates the accuracy and applicability of a developed workflow which allows the creation of an affordable, metal, individualized navigational template by integration of FEA in the design and surgical planning process.

The integration of FEA in the pedicle screw intraoperative navigation was investigated by Van den Abbeele et al. (30), however, the application of FEA in the design process of a navigational template in spine surgery by integrating the patient bone mineral density related material properties is new. The results of the simulations showed that the convergent S1 insertion is significantly stiffer than the divergent ALA insertion. This finding is supported by cadaveric experimental studies (31, 32) and clinical experience as well (33). The biomechanical difference of the convergent and divergent insertions rely on the differences in the local bone mineral densities (34).

The combination of the 3D-printing technology and cobalt-chrome casting makes the manufacturing process more affordable. Investment casting of cobalt-chrome is a widely used technology in dental laboratories (35). 3D-printed patterns for casting is an accepted method in dentistry (35, 36); however, its application in spine surgery navigational templates is novel. The production of individualized metal navigational templates for screw insertion can be achieved via selective laser sintering 3D-printing technology of titanium-based alloys (37), but at a higher cost and lower accessibility compared to dental casting. Metal templates are robust, resistant to damage, and can also be easily autoclaved (37).

It is widely accepted in the literature to use cadavers for testing, evaluating the fitting accuracy of a navigational template (38). FDM technology can produce geometrically accurate spine physical models (39) and the different designs can be tested as well as the drilling accuracy can be evaluated. The use of FDM models for design process evaluation and development is advantageous due to the possibility to include retrospective patient imaging data with complex anatomical/geometrical variation (deformities, tumors, etc.) which is extremely difficult to control and integrate in the case of cadaveric specimen studies. Based on our FEA results, the S1 screw insertion's surgical plan and drill template position is recommended for surgical implementation. Despite the fact that the virtual screw axis and the drill bit centerline are not colinear and coincident (3D angle α > 0) according to the Gertzbein-Robbins scale (40) the template theoretically allows an accurate (grade A) screw insertion (Figure 6). We present the surgical technique for the screw insertion with the developed drill template (Figure 8 and Supplementary Video 1). The suggested screw insertion surgical technique uses the philosophy of the minimally invasive pedicle screw insertion techniques (MIS) by using a Kirschner wire, cannulated tap, and pedicle screw. This technique can easily be performed by any spine surgeon familiar with MIS pedicle screw insertion.


Figure 8. Proposed surgical technique for safe and accurate screw insertion in convergent position. (A) Transparent surface mesh of the patient's sacrum with the broken screw. (B) Section plane dimension and orientation, and drill template position on the sacrum. (C) Stainless steel cylinder inlet connected to the navigation template for the drill bit. (D) Stainless steel cylinder inlet connected to the template for the Kirschner wire. (E) The inlet cylinder and template are removed, the Kirschner wire's position is unchanged. (F) Cannulated tap is introduced along the Kirschner wire. (G) Cannulated pedicle screw is introduced in the sacrum along the Kirschner wire. (H) Final position of the screw. (I) Transparent surface mesh of the sacrum with the broken and convergently inserted pedicle screw geometry.

Limitations of this study include the fact that the developed template is presented using a single case, however the workflow can be applied for different parts of the spine with different geometrical difficulties/pathologies. The presented FEA models' loading conditions are simplified as well as the material property assignments; more complex FEA investigations would be desirable. In the future, a randomized study of specific subtypes of spinal pathologies (tumors, deformities, etc.) with a larger sample size would be preferred to demonstrate the clinical efficacy and cost-effectiveness of the developed methodology.


A patient-specific template for pedicle screw insertion allows the surgeon to insert the screw into its optimal position. The advantages of our technique compared to the conventional surgical navigation tools are the affordability, the potential to reduce intraoperative X-ray exposure, and the possibility for the consideration of patient-specific bone geometry and biomechanics. This new patient- and condition-specific approach can be widely used in revision spine surgeries or in challenging primary cases after its further clinical validations.

Data Availability Statement

The datasets generated for this study are available on request to the corresponding author.

Ethics Statement

The studies involving human participants were reviewed and approved by National Ethics Committee of Hungary and the National Institute of Pharmacy and Nutrition (reference number: OGYÉI/163-4/2019. The patients/participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.

Author Contributions

PEE, AL, PPV, MB, and DL: research design. PEE, AL, MB, BH, AJP, and LK: acquisition of data. PEE, AL, MB, BH, AJP, and DL: analysis and/or interpretation of data. All authors: drafting the paper or revising it critically and approval of the submitted and final versions.


The financial support from the Hungarian Scientific Research Fund (award number: OTKA FK123884) was gratefully acknowledged.

Conflict of Interest

MB was employed by the company DO3D Innovations Ltd.

The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.


The author would like to thank to Julia Szita, Ph.D. student, National Center for Spinal Disorders, Budapest for constructive criticism and the proofreading of the manuscript.

Supplementary Material

The Supplementary Material for this article can be found online at:


1. Gaines RW Jr. The use of pedicle-screw internal fixation for the operative treatment of spinal disorders. J Bone Joint Surg Am. (2000) 82:1458. doi: 10.2106/00004623-200010000-00013

CrossRef Full Text | Google Scholar

2. Kast E, Mohr K, Richter H-P, Börm W. Complications of transpedicular screw fixation in the cervical spine. Eur Spine J. (2006) 15:327–34. doi: 10.1007/s00586-004-0861-7

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Faraj AA, Webb JK. Early complications of spinal pedicle screw. Eur Spine J. (1997) 6:324–6. doi: 10.1007/BF01142678

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Merloz P, Tonetti J, Pittet L, Coulomb M, Lavallee S, Troccaz J, et al. Computer-assisted spine surgery. Comput Aided Surg Off J Int Soc Comput Aided Surg. (1998) 3:297–305. doi: 10.3109/10929089809148150

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Foley KT, Simon DA, Rampersaud YR. Virtual fluoroscopy: computer-assisted fluoroscopic navigation. Spine (Phila Pa 1976). (2001) 26:347–351. doi: 10.1097/00007632-200102150-00009

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Nolte L-P, Slomczykowski MA, Berlemann U, Strauss MJ, Hofstetter R, Schlenzka D, et al. A new approach to computer-aided spine surgery: fluoroscopy-based surgical navigation. Eur Spine J. (2000) 9:S078–88. doi: 10.1007/PL00010026

PubMed Abstract | CrossRef Full Text | Google Scholar

7. Ishikawa Y, Kanemura T, Yoshida G, Matsumoto A, Ito Z, Tauchi R, et al. Intraoperative, full-rotation, three-dimensional image (O-arm)–based navigation system for cervical pedicle screw insertion. J Neurosurg Spine. (2011) 15:472–8. doi: 10.3171/2011.6.SPINE10809

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Radermacher K, Portheine F, Anton M, Zimolong A, Kaspers G, Rau G, et al. Computer assisted orthopaedic surgery with image based individual templates. Clin Orthop Relat Res. (1998) 354:28–38. doi: 10.1097/00003086-199809000-00005

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Putzier M, Strube P, Cecchinato R, Lamartina C, Hoff EK. A new navigational tool for pedicle screw placement in patients with severe scoliosis: a pilot study to prove feasibility, accuracy, and identify operative challenges. Clin Spine Surg. (2017) 30:E430–9. doi: 10.1097/BSD.0000000000000220

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Jiang L, Dong L, Tan M, Qi Y, Yang F, Yi P, et al. A modified personalized image-based drill guide template for atlantoaxial pedicle screw placement: a clinical study. Med Sci Monit. (2017) 23:1325–33. doi: 10.12659/MSM.900066

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Garg B, Gupta M, Singh M, Kalyanasundaram D. Outcome and safety analysis of 3D-printed patient-specific pedicle screw jigs for complex spinal deformities: a comparative study. Spine J. (2019) 19:56–64. doi: 10.1016/j.spinee.2018.05.001

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Cristina Mancarella AL, Roberto Delfini FG, i, Mancarella C, Gregori F, Delfini R. Spinal neuronavigation and 3D-printed tubular guide for pedicle screw placement: a really new tool to improve safety and accuracy of the surgical technique? J Spine. (2015) 4:1–3. doi: 10.4172/2165-7939.1000e118

CrossRef Full Text | Google Scholar

13. Guo F, Dai J, Zhang J, Ma Y, Zhu G, Shen J, et al. Individualized 3D printing navigation template for pedicle screw fixation in upper cervical spine. PLoS ONE. (2017) 12:e0171509. doi: 10.1371/journal.pone.0171509

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Choonara YE, du Toit LC, Kumar P, Kondiah PPD, Pillay V. 3D-printing and the effect on medical costs: a new era? Expert Rev Pharmacoecon Outcomes Res. (2016) 16:23–32. doi: 10.1586/14737167.2016.1138860

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Marro A, Bandukwala T, Mak W. Three-dimensional printing and medical imaging: a review of the methods and applications. Curr Probl Diagn Radiol. (2016) 45:2–9. doi: 10.1067/j.cpradiol.2015.07.009

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Weller C, Kleer R, Piller FT. Economic implications of 3D printing: market structure models in light of additive manufacturing revisited. Int J Prod Econ. (2015) 164:43–56. doi: 10.1016/j.ijpe.2015.02.020

CrossRef Full Text | Google Scholar

17. Lamartina C, Cecchinato R, Fekete Z, Lipari A, Fiechter M, Berjano P. Pedicle screw placement accuracy in thoracic and lumbar spinal surgery with a patient-matched targeting guide: a cadaveric study. Eur Spine J. (2015) 24:937–41. doi: 10.1007/s00586-015-4261-y

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Sheha ED, Gandhi SD, Colman MW. 3D printing in spine surgery. Ann Transl Med. (2019) 7:S164. doi: 10.21037/atm.2019.08.88

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Castro-Mateos I, Pozo JM, Lazary A, Frangi A. 3D Vertebra segmentation by feature selection active shape model. In: Recent Advances in Computational Methods and Clinical Applications for Spine Imaging. Cham: Springer (2015). p. 241–5. doi: 10.1007/978-3-319-14148-0_22

CrossRef Full Text | Google Scholar

20. van Rijsbergen M, van Rietbergen B, Barthelemy V, Eltes P, Lazáry Á, Lacroix D, et al. Comparison of patient-specific computational models vs. clinical follow-up, for adjacent segment disc degeneration and bone remodelling after spinal fusion. PLoS ONE. (2018) 13:e0200899. doi: 10.1371/journal.pone.0200899

CrossRef Full Text | Google Scholar

21. Aryanto KYE, Oudkerk M, van Ooijen PMA. Free DICOM de-identification tools in clinical research: functioning and safety of patient privacy. Eur Radiol. (2015) 25:3685–95. doi: 10.1007/s00330-015-3794-0

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Rho JY, Hobatho MC, Ashman RB. Relations of mechanical properties to density and CT numbers in human bone. Med Eng Phys. (1995) 17:347–55. doi: 10.1016/1350-4533(95)97314-F

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Kopperdahl DL, Morgan EF, Keaveny TM. Quantitative computed tomography estimates of the mechanical properties of human vertebral trabecular bone. J Orthop Res. (2002) 20:801–5. doi: 10.1016/S0736-0266(01)00185-1

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Couteau B, Hobatho M-C, Darmana R, Brignola J-C, Arlaud J-Y. Finite element modelling of the vibrational behaviour of the human femur using CT-based individualized geometrical and material properties. J Biomech. (1998) 31:383–6. doi: 10.1016/S0021-9290(98)00018-9

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Charriere E, Sirey F, Zysset PK. A finite element model of the L5-S1 functional spinal unit: development and comparison with biomechanical tests in vitro. Comput Methods Biomech Biomed Eng. (2003) 6:249–61. doi: 10.1080/10255840310001606099

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Welsch G, Boyer R, Collings EW. Materials Properties Handbook: Titanium Alloys. Russell Township, OH: ASM International (1993).

Google Scholar

27. Fan Y, Du J, Zhang J, Liu S, Xue X, Huang Y, et al. Comparison of accuracy of pedicle screw insertion among 4 guided technologies in spine surgery. Med Sci Monit Int Med J Exp Clin Res. (2017) 23:5960. doi: 10.12659/MSM.905713

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Li Y, Lin J, Wang Y, Luo H, Wang J, Lu S, et al. Comparative study of 3D printed navigation template-assisted atlantoaxial pedicle screws versus free-hand screws for type II odontoid fractures. Eur Spine J. (2020). doi: 10.1007/s00586-020-06644-9. [Epub ahead of print].

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Wen ZJ, Gao ZC, Lu T, Wang YB, Liang H, He XJ. Comparison of the effect of navigation template assisted spinal pedicle fixation and traditional pedicle screw fixation: a meta-analysis. Zhongguo Gu Shang China J Orthop Traumatol. (2018) 31:1069–76. doi: 10.3969/j.issn.1003-0034.2018.11.017

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Van den Abbeele M, Valiadis JM, Lima LVPC, Khalifé P, Rouch P, Skalli W. Contribution to FE modeling for intraoperative pedicle screw strength prediction. Comput Methods Biomech Biomed Engin. (2018) 21:13–21. doi: 10.1080/10255842.2017.1414200

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Von Strempel A, Trenkmann S, Krönauer I, Kirsch L, Sukopp C. The stability of bone screws in the os sacrum. Eur Spine J. (1998) 7:313–20. doi: 10.1007/s005860050081

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Zhu Q, Lu WW, Holmes AD, Zheng Y, Zhong S, Leong JCY. The effects of cyclic loading on pull-out strength of sacral screw fixation: an in vitro biomechanical study. Spine (Phila Pa 1976). (2000) 25:1065–. doi: 10.1097/00007632-200005010-00005

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Kato M, Taneichi H, Suda K. Advantage of pedicle screw placement into the sacral promontory (tricortical purchase) on lumbosacral fixation. J Spinal Disord Tech. (2015) 28:E336–42. doi: 10.1097/BSD.0b013e31828ffc70

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Sabry FF, Xu R, Nadim Y, Ebraheim NA. Bone density of the first sacral vertebra in relation to sacral screw placement: a computed tomography study. Orthopedics. (2001) 24:475–7.

PubMed Abstract | Google Scholar

35. Misch CE. Dental Implant Prosthetics-E-Book. St. Louis, MO: Elsevier Health Sciences (2004).

Google Scholar

36. Bhaskaran E, Azhagarasan NS, Miglani S, Ilango T, Krishna GP, Gajapathi B. Comparative evaluation of marginal and internal gap of Co–Cr copings fabricated from conventional wax pattern, 3D printed resin pattern and DMLS tech: an in vitro study. J Indian Prosthodont Soc. (2013) 13:189–95. doi: 10.1007/s13191-013-0283-5

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Wang D, Wang Y, Wang J, Song C, Yang Y, Zhang Z, et al. Design and fabrication of a precision template for spine surgery using selective laser melting (SLM). Materials (Basel). (2016) 9:608. doi: 10.3390/ma9070608

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Sallent A, Ramírez M, Catalá J, Rodríguez-Baeza A, Bagó J, de Albert M, et al. Precision and safety of multilevel cervical transpedicular screw fixation with 3D patient-specific guides; a cadaveric study. Sci Rep. (2019) 9:15686. doi: 10.1038/s41598-019-51936-w

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Eltes PE, Kiss L, Bartos M, Gyorgy ZM, Csakany T, Bereczki F, et al. Geometrical accuracy evaluation of an affordable 3D printing technology for spine physical models. J Clin Neurosci. (2020) 72:438–46. doi: 10.1016/j.jocn.2019.12.027

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Gertzbein SD, Robbins SE. Accuracy of pedicular screw placement in vivo. Spine (Phila Pa 1976). (1990) 15:11–4. doi: 10.1097/00007632-199001000-00004

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: 3D printing, computed tomography, navigation, finite element simulation, spine surgery, surgical guidance/navigation

Citation: Eltes PE, Bartos M, Hajnal B, Pokorni AJ, Kiss L, Lacroix D, Varga PP and Lazary A (2021) Development of a Computer-Aided Design and Finite Element Analysis Combined Method for Affordable Spine Surgical Navigation With 3D-Printed Customized Template. Front. Surg. 7:583386. doi: 10.3389/fsurg.2020.583386

Received: 14 July 2020; Accepted: 18 December 2020;
Published: 25 January 2021.

Edited by:

Ralph Jasper Mobbs, University of New South Wales, Australia

Reviewed by:

Alexandre Terrier, École Polytechnique Fédérale de Lausanne, Switzerland
Konstantinos Markatos, Salamina Medical Center, Greece

Copyright © 2021 Eltes, Bartos, Hajnal, Pokorni, Kiss, Lacroix, Varga and Lazary. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Peter Endre Eltes,;