In vivo visualization of cardiovascular structures is possible using medical images. However, one has to realize that the resulting 3D geometries correspond to in vivo conditions. This entails an internal stress state to be present in the in vivo measured geometry of e.g. a blood vessel due to the presence of the blood pressure. In order to correct for this in vivo stress, this paper presents an inverse method to restore the original zero-pressure geometry of a structure, and to recover the in vivo stress field of the final, loaded structure. The proposed backward displacement method is able to solve the inverse problem iteratively using fixed point iterations, but can be significantly accelerated by a quasi-Newton technique in which a least-squares model is used to approximate the inverse of the Jacobian. The here proposed backward displacement method allows for a straightforward implementation of the algorithm in combination with existing structural solvers, even if the structural solver is a black box, as only an update of the coordinates of the mesh needs to be performed.
@article{M2AN_2013__47_4_1059_0,
author = {Bols, Joris and Degroote, Joris and Trachet, Bram and Verhegghe, Benedict and Segers, Patrick and Vierendeels, Jan},
title = {Inverse modelling of image-based patient-specific blood vessels: zero-pressure geometry and in vivo stress incorporation},
journal = {ESAIM: Mathematical Modelling and Numerical Analysis - Mod\'elisation Math\'ematique et Analyse Num\'erique},
volume = {47},
year = {2013},
pages = {1059-1075},
doi = {10.1051/m2an/2012057},
language = {en},
url = {http://dml.mathdoc.fr/item/M2AN_2013__47_4_1059_0}
}
Bols, Joris; Degroote, Joris; Trachet, Bram; Verhegghe, Benedict; Segers, Patrick; Vierendeels, Jan. Inverse modelling of image-based patient-specific blood vessels: zero-pressure geometry and in vivo stress incorporation. ESAIM: Mathematical Modelling and Numerical Analysis - Modélisation Mathématique et Analyse Numérique, Tome 47 (2013) pp. 1059-1075. doi : 10.1051/m2an/2012057. http://gdmltest.u-ga.fr/item/M2AN_2013__47_4_1059_0/
[1] , , , , , , and , An integrated framework to quantitatively link mouse-specific hemodynamics to aneurysm formation in angiotensin II-infused ApoE -/- mice. Annal. Biomed. Eng. 39 (2011) 2430-2444.
[2] , , , , , and , On coupling a lumped parameter heart model and a three-dimensional finite element aorta model. Annal. Biomed. Eng. 37 (2009) 2153-2169.
[3] , , , and , Inverse modelling of an aneurysms stiffness using surrogate-based optimization and fluid-structure interaction simulations. Struct. Multidiscip. Optim. (2012) 1-13. | Zbl 1274.74102
[4] , and , Inverse elastostatic stress analysis in pre-deformed biological structures: Demonstration using abdominal aortic aneurysms. J. Biomech. 40 (2007) 693-6.
[5] , , , , and , Patient-specific initial wall stress in abdominal aortic aneurysms with a backward incremental method. J. Biomech. 40 (2007) 1081-1090.
[6] , , and , Prestressing in finite deformation abdominal aortic aneurysm simulation. J. Biomech. 42 (2009) 1732-1739.
[7] , , , , , and , Initial stress and nonlinear material behavior in patient-specific AAA wall stress analysis. J. Biomech. 42 (2009) 1713-1719.
[8] M.A.G. Merkx, M. van 't Veer, L. Speelman, M. Breeuwer, J. Buth and F.N. van de Vosse, Importance of initial stress for abdominal aortic aneurysm wall motion: Dynamic MRI validated finite element analysis. J. Biomech. 42 (2009) 2369-2373.
[9] and , Computational methods for inverse finite elastostatics. Comput. Methods Appl. Mech. Eng. 136 (1996) 47-57. | Zbl 0918.73117
[10] and , Computational methods for inverse deformations in quasi-incompressible finite elasticity. Inter. J. Numer. Methods Eng. 43 (1998) 821-838. | Zbl 0937.74064
[11] , and , Finite element modelling of inverse design problems in large deformations anisotropic hyperelasticity. Inter. J. Numer. Methods Eng. 74 (2008) 894-910. | MR 2389139 | Zbl 1158.74369
[12] , , and , The quasi-newton least squares method: A new and fast secant method analyzed for linear systems. SIAM J. Numer. Anal. 47 (2009) 2347-2368. | MR 2519606 | Zbl 1197.65041
[13] , and , Performance of a new partitioned procedure versus a monolithic procedure in fluid-structure interaction. Comput. Struct. 87 (2009) 793-801.
[14] , , and , On the similarities between the quasi-newton inverse least squares method and GMRes. SIAM J. Numer. Anal. 47 (2010) 4660-4679. | MR 2595053 | Zbl 1209.65048
[15] , , , and . Performance of partitioned procedures in fluid-structure interaction. Comput. Struct. 88 (2010) 446-457.
[16] , and , Non-invasive determination of zero-pressure geometry of arterial aneurysms. Annal. Biomedical Eng. 34 (2006) 1414-1419.
[17] , and , A computational strategy for prestressing patient-specific biomechanical problems under finite deformation. Inter. J. Numer. Methods Biomedical Eng. 26 (2010) 52-72. | Zbl 1180.92009
[18] , , , , and , Numerical framework for patient-specific computational modelling of vascular tissue. Inter. J. Numer. Methods Biomedical Eng. 26 (2010) 35-51. | Zbl 1180.92036
[19] , , , , , , , , , and , Initial stress in biomechanical models of atherosclerotic plaques. J. Biomech. 44 (2011) 2376-2382.
[20] , and , Computational modeling of mechanical anisotropy in the cornea and sclera. J. Cataract Refract. Surg. 31 (2005) 136-45.
[21] , , , , , and , Computational vascular fluid-structure interaction: methodology and application to cerebral aneurysms. Biomech. Model. Mechanobiol. 9 (2010) 481-498.
[22] and , Blood vessel tissue prestress modeling for vascular fluid-structure interaction simulation. Finite Elem. Anal. Des. 47 (2011) 593-599. | MR 2786927
[23] , , , , , and , Analysis of prolapse in cardiovascular stents: A constitutive equation for vascular tissue and finite-element modelling. J. Biomech. Eng. 125 (2003) 692-699.
[24] http://www.pyformex.org
[25] , , , , and , Full-hexahedral structured meshing for image-based computational vascular modeling. Medical Eng. Phys. 33 (2011) 1318-1325.