RESUMEN
Temperature-based death time estimation is based either on simple phenomenological models of corpse cooling or on detailed physical heat transfer models. The latter are much more complex but allow a higher accuracy of death time estimation, as in principle, all relevant cooling mechanisms can be taken into account.Here, a complete workflow for finite element-based cooling simulation is presented. The following steps are demonstrated on a CT phantom: Computer tomography (CT) scan Segmentation of the CT images for thermodynamically relevant features of individual geometries and compilation in a geometric computer-aided design (CAD) model Conversion of the segmentation result into a finite element (FE) simulation model Computation of the model cooling curve (MOD) Calculation of the cooling time (CTE) For the first time in FE-based cooling time estimation, the steps from the CT image over segmentation to FE model generation are performed semi-automatically. The cooling time calculation results are compared to cooling measurements performed on the phantoms under controlled conditions. In this context, the method is validated using a CT phantom. Some of the phantoms' thermodynamic material parameters had to be determined via independent experiments.Moreover, the impact of geometry and material parameter uncertainties on the estimated cooling time is investigated by a sensitivity analysis.
Asunto(s)
Simulación por Computador , Análisis de Elementos Finitos , Cambios Post Mortem , Diseño Asistido por Computadora , Estudios de Factibilidad , Humanos , Modelos Biológicos , Tomografía Computarizada Multidetector , Fantasmas de ImagenRESUMEN
Computer Tomographic (CT) image data have become a standard basis for structural analyses of bony organs. In this context, regression functions between stiffness components and Hounsfields units (HU) from CT, related to X-ray attenuation coefficients, are widely used for the definition of the (actually inhomogeneous and anisotropic) material behavior inside the organ. Herein, we suggest to derive the functional dependence of the fully orthotropic stiffness tensors on the Hounsfield units from the physical information contained in the X-ray attenuation coefficients: (i) Based on voxel average rules for the X-ray attenuation coefficients, we assign to each voxel the volume fraction occupied by water (marrow) and that occupied by solid bone matrix. (ii) By means of a continuum micromechanics representation for bone, which is based on voxel-invariant (species and whole bone-specific) stiffness properties of solid bone matrix and of water, we convert the aforementioned volume fractions into voxel-specific orthotropic stiffness tensor components. The micromechanics model, in combination with the average rule for X-ray attenuation coefficients, predicts a quasi-linear relationship between axial Young's modulus and HU, and highly nonlinear relationships for both circumferential and radial Young's moduli as well as for the shear moduli in all principal material directions. Corresponding whole-organ Finite Element (FE) analyses of a partially edentulous human mandible characterized by atrophy of the alveolar ridge show that volumetric strain concentrations/peaks within the organ are decreased when considering material anisotropy, and increased when considering material inhomogeneity.