RESEARCH ARTICLE
F.E.M. Stress-Investigation of Scolios Apex
A. Daghighi1, *, H. Tropp2, #, N. Dahlström3, A. Klarbring4
Article Information
Identifiers and Pagination:
Year: 2018Volume: 12
First Page: 51
Last Page: 71
Publisher ID: TOBEJ-12-51
DOI: 10.2174/1874120701812010051
Article History:
Received Date: 24/04/2018Revision Received Date: 15/07/2018
Acceptance Date: 16/07/2018
Electronic publication date: 31/07/2018
Collection year: 2018
open-access license: This is an open access article distributed under the terms of the Creative Commons Attribution 4.0 International Public License (CC-BY 4.0), a copy of which is available at: https://creativecommons.org/licenses/by/4.0/legalcode. This license permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Abstract
Background:
In scoliosis, kypholordos and wedge properties of the vertebrae should be involved in determining how stress is distributed in the vertebral column. The impact is logically expected to be maximal at the apex.
Aim:
To introduce an algorithm for constructing artificial geometric models of the vertebral column from DICOM stacks, with the ultimate aim to obtain a formalized way to create simplistic models, which enhance and focus on wedge properties and relative tilting.
Material/Methods:
Our procedure requires parameter extraction from DICOM image-stacks (with PACS,IDS-7), mechanical FEM-modelling (with Matlab and Comsol). As a test implementation, models were constructed for five patients with thoracal idiopathic scoliosis with varying apex rotation. For a selection of load states, we calculated a response variable which is based upon distortion energy.
Results:
For the test implementation, pairwise t-tests show that our response variable is non-trivial and that it is chiefly sensitive to the transversal stresses (transversal stresses where of main interest to us, as opposed to the case of additional shear stresses, due to the lack of explicit surrounding tissue and ligaments in our model). Also, a pairwise t-test did not show a difference (n = 25, p-value≈0.084) between the cases of isotropic and orthotropic material modeling.
Conclusion:
A step-by-step description is given for a procedure of constructing artificial geometric models from chest CT DICOM-stacks, such that the models are appropriate for semi-global stress-analysis, where the focus is on the wedge properties and relative tilting. The method is inappropriate for analyses where the local roughness and irregularities of surfaces are wanted features. A test application hints that one particular load state possibly has a high correlation to a certain response variable (based upon distortion energy distribution on a surface of the apex), however, the number of patients is too small to draw any statistical conclusions.
