RESUMEN
OBJECTIVE: We propose a procedure for calibrating 4 parameters governing the mechanical boundary conditions (BCs) of a thoracic aorta (TA) model derived from one patient with ascending aortic aneurysm. The BCs reproduce the visco-elastic structural support provided by the soft tissue and the spine and allow for the inclusion of the heart motion effect. METHODS: We first segment the TA from magnetic resonance imaging (MRI) angiography and derive the heart motion by tracking the aortic annulus from cine-MRI. A rigid-wall fluid-dynamic simulation is performed to derive the time-varying wall pressure field. We build the finite element model considering patient-specific material properties and imposing the derived pressure field and the motion at the annulus boundary. The calibration, which involves the zero-pressure state computation, is based on purely structural simulations. After obtaining the vessel boundaries from the cine-MRI sequences, an iterative procedure is performed to minimize the distance between them and the corresponding boundaries derived from the deformed structural model. A strongly-coupled fluid-structure interaction (FSI) analysis is finally performed with the tuned parameters and compared to the purely structural simulation. RESULTS AND CONCLUSION: The calibration with structural simulations allows to reduce maximum and mean distances between image-derived and simulation-derived boundaries from 8.64 mm to 6.37 mm and from 2.24 mm to 1.83 mm, respectively. The maximum root mean square error between the deformed structural and FSI surface meshes is 0.19 mm. This procedure may prove crucial for increasing the model fidelity in replicating the real aortic root kinematics.
RESUMEN
Pedicle screw fixation is extensively performed to treat spine injuries or diseases and it is common for thoracolumbar fractures. Post-operative complications may arise from this surgery leading to back pain or revisions. Finite element (FE) models could be used to predict the outcomes of surgeries but should be verified when both simplified and realistic designs of screws are used. The aim of this study was to generate patient-specific Computed Tomography (CT)-based FE models of human vertebrae with two pedicle screws, verify the models, and use them to evaluate the effect of the screws' size and geometry on the mechanical properties of the screws-vertebra structure. FE models of the lumbar vertebra implanted with two pedicle screws were created from anonymized CT-scans of three patients. Compressive loads were applied to the head of the screws. The mesh size was optimized for realistic and simplified geometry of the screws with a mesh refinement study. Finally, the optimal mesh size was used to evaluate the sensitivity of the model to changes in screw's size (diameter and length) and geometry (realistic or simplified). For both simplified and realistic models, element sizes of 0.6 mm in the screw and 1.0 mm in the bone allowed to obtain relative differences of approximately 5% or lower. Changes in screw's length resulted in 4-10% differences in maximum deflection, 1-6% differences in peak stress in the screws, 10-22% differences in mean strain in the bone around the screw; changes in screw's diameter resulted in 28-36% differences in maximum deflection, 6-27% differences in peak stress in the screws, and 30-47% differences in mean strain in the bone around the screw. The maximum deflection predicted with realistic or simplified screws correlated very well (R 2 = 0.99). The peak stress in screws with realistic or simplified design correlated well (R 2 = 0.82) but simplified models underestimated the peak stress. In conclusion, the results showed that the diameter of the screw has a major role on the mechanics of the screw-vertebral structure for each patient. Simplified screws can be used to estimate the mechanical properties of the implanted vertebrae, but the systematic underestimation of the peak stress should be considered when interpreting the results from the FE analyses.
