Design
A three-dimensional finite element analysis.
Time and settings
The experiment was conducted in 2019 at the Digital Medical Center of Inner Mongolia Medical University.
Subjects
One patient with lumbar sacral vertebrae was selected, a 16-year-old male of Han nationality, with a body mass of 60 kg, height of 1.75 m and body mass index (19.59 kg/m2). The volunteer had no previous history of lumbar trauma and related diseases. We obtain written informed consent from patient and his guardian. The research was approved by the Ethics Committee of Inner Mongolia Medical University (approval number YKD2018015) on March 5, 2018.
Methods
Acquisition of two-dimensional data
In this study, the Lightspeed dual-source 64-slice spiral CT (GE, USA) was used to scan the upper edge of the L1 vertebra to the lower edge of the L5 vertebra using the volume scan mode. Lie on your back with your head forward. Scanning parameters: slice thickness 5 mm, layer spacing 5 mm, scanning field 20 cm2, matrix 512 × 512, standard algorithm reconstruction. At the end of the scan, the data is divided into layer thickness 0.625 mm, layer spacing 0.625 mm and saved in DICOM format (Figure 1).
3D reconstruction of 2D data
Import the original DICOM data into Mimics16.0 (Materialise Belgium), through the steps of threshold setting, image segmentation, region growth, image filling and so on, the L4-sacrum is segmented and smoothed, and then the Mask, of each vertebra is obtained through the Calculate3D module to get the 3D image of the vertebra. Then use the 3-matic software in Mimics16.0 software to optimize the model and export it in STL format (Figure 2).
Reconstruction of 3D solid model
After Mimics image extraction and 3-matic model repair, the complete vertebral models of L5 and sacrum were obtained. At present, the appearance of each model is close to the real vertebra, with regular vertebrae, transverse process, spinous process, upper and lower articular process and other main structures. After introducing the model into Geomagicstudio (Geomagic, USA) the model unit is first selected, which is millimeter by default, and the system automatically prompts for grid doctor examination. After this processing, the polygon model of feature and smoothness balance is obtained. Simplify the sacrum model, use the plane clipping function, retain only the part of the sacrum in contact with S1, and carefully close the clipping plane, using the grid doctor to check to avoid causing the model to be unclosed (Figure 3).
Making of lumbar intervertebral disc model
The real lumbar vertebra is a complex motion system composed of vertebrae, intervertebral discs, ligaments, muscles and so on. If you want to carry out finite element analysis, the most simplified model should at least include vertebrae and intervertebral discs. It is very easy to make intervertebral disc models with Geomagic Studio software. The shape of the intervertebral disc model made in the experiment is lifelike, showing a disk shape similar to the solid, the leading edge is thick, the trailing edge is flat, and the edge is smooth without edges and corners. The nucleus pulposus is located slightly behind the center of the intervertebral disc on the cross section and in the middle of the intervertebral disc on the coronal plane. After measurement, the relevant parameters of intervertebral disc and nucleus pulposus are consistent with the real human body. After assembling, the intervertebral disc and adjacent vertebrae can be closely attached, the relative position of intervertebral disc and vertebra is correct, and there are no pathological manifestations such as intervertebral disc herniation, which can meet the needs of further experiments (Table 1 and Figure 4).
Assembly and volume meshing of lumbosacral transitional vertebrae
Through the previously described work, a high-quality triangular model of the lumbar spine and intervertebral disc was obtained, but it is not a CAD model and it needs to be further transformed into a NURBS surface solid model and exported in IGES, sat, X_T, and other formats in order for the data to be smoothly imported into the WorkBench finite element analysis software (ANSYS, USA). The specific workflow is as follows, using the L4/5 disc as an example:
(1) After importing the model into Geomagic Studio, the Exact Surfaces Phase is clicked and Exact Surface is clicked.
(2) Click the Detect Curvature method and select Auto Estimate and Simplify Contour.
(3) Using the Construct Patch tool, the area was re-divided. There were occasional errors such as intersecting paths, different sizes of surface patches in the structure, or mismatch between the distribution pattern and the actual shape of the model, and further manual editing of the surface patches was required.
(4) Use the Move Panel function under the Move option to edit the patch.
(5) Construct the grid, click on the Construct Grid button, click on Repair Intersection Area, and check the geometry option, and use a resolution of 20 to get the model image that covers the grid.
(6) Click on the fitted surface to transform the model formed by the grid into a NURBS surface and optimize the smoothness option to further smooth the model. Then, save the NURBS surface in IGES or STEP format for the next step.
Import of finite element model into WorkBench and calculations
The assembled model was imported into WorkBench as described previously[11-13]. The corresponding material properties were assigned to the vertebrae, the annulus fibrosus, and the nucleus pulposus. Next, the assembled model and properties were assigned, the mesh was defined, and the contact surface, boundary conditions, and load were determined. For loading, a 2 N·mm force was applied for the analysis and calculation under the conditions of flexion, extension, lateral bending to the left, and lateral bending to the right. The stress and overall deformation of the model as a whole and the L5 and S1 finite element models were calculated when the flexion, extension, and bending moments were 2 N·mm (Table 2 and Figure 5).
Main outcome measures
To analyze the mechanical properties and stress-strain rules of lumbar vertebrae under simulated normal exercise conditions.