The goal of this paper is to study computationally how blood vessels adapt when they are exposed to a mechanobiological insult, namely, a sudden change of their biomechanical conditions such as proteolytic injuries or implantation.Adaptation occurs through growth and remodeling (G&R), consisting of mass production or removal of structural proteins, such as collagen, until restoring the initial homeostatic biomechanical conditions. In some circumstances, the initial conditions can never be recovered, and arteries evolve towards unstable pathological conditions, such as aneurysms, which are responsible for significant morbidity and mortality. Therefore, computational predictions of G&R under different circumstances can be helpful in understanding fundamentally how arterial pathologies progress. For that, we have developed a low-cost open-source finite-element 2D axisymmetric shell model (FEM) of the arterial wall. The constitutive equations for static equilibrium used to model the stress-strain behavior and the G&R response are expressed within the homogenized constrained mixture theory. The originality is to integrate the layer-specific behavior of both arterial layers (media and adventitia) into the model. Considering different mechanobiological insults, our results show that the resulting arterial dilatation is strongly correlated with the media thickness. The adaptation to stent implantation is particularly interesting. For large stent oversizing ratios, the artery cannot recover from the mechanobiological insult and dilates forever, whereas dilatation stabilizes after a transient period for more moderate oversizing ratios. We also show that stent implantation induces a different response in an aneurysm or in a healthy artery, the latter yielding more unstable G&R. Finally, our G&R model can efficiently predict, with very low computational cost, fundamental aspects of arterial adaptation induced by clinical procedures. KEYWORDSarterial growth and remodeling, axisymmetric shell model, homogenized constrained mixture theory, layer-specific behavior, stent implantation Int J Numer Meth Biomed Engng. 2020;36:e3282.wileyonlinelibrary.com/journal/cnm

Evolution of mechanical and structural properties in the Ascending Thoracic Aorta (ATA) is the results of complex mechanobiological processes. In this work, we address some numerical challenges in order to elaborate computational models of these processes. For that, we extend the state of the art of homogenized constrained mixture (hCM) models. In these models, prestretches are assigned to the mixed constituents in order to ensure local mechanical equilibrium macroscopically, and to maintain a homeostatic level of tension in collagen fibers microscopically. Although the initial prestretches were assumed as homogeneous in idealized straight tubes, more elaborate prestretch distributions need to be considered for curved geometrical models such as patient-specific ATA. Therefore, we introduce prestretches having a three-dimensional gradient across the ATA geometry in the homeostatic reference state. We test different schemes with the objective to ensure stable growth and remodeling (G&R) simulations on patient-specific curved vessels. In these simulations, aneurysm progression is triggered by tissue changes in the constituents such as mass degradation of intramural elastin. The results show that the initial prestretches are not only critical for the stability of numerical simulations, but they also affect the G&R response. Eventually, we submit that initial conditions required for G&R simulations need to be identified regionally for ensuring realistic patient-specific predictions of aneurysm progression.

Progressive aneurysmal dilatation is a well-recognized complication in patients with chronic type B aortic dissection (cTBAD), which may lead to a delayed rupture and create a life-threatening condition. However, our understanding of such aortic expansion in cTBAD remains weak. In the present paper, we propose to use numerical simulations to study the role of growth and remodeling (G&R) in aneurysmal dilatation after cTBAD. We set up a 3D finite-element model of G&R for aortic dissection within an open-source code. Constitutive equations, momentum balance equations, and equations related to the mechanobiology of the artery were formulated based on the homogenized constrained mixture theory. The model was first applied to idealized aortic geometries with cylindrical and toric shapes to demonstrate its feasibility and efficiency. The model was then applied to a patient-specific aortic segment to show its potential in more relevant and complex patient-specific clinical applications. It was found that the G&R tends to naturally trigger the aneurysmal dilatation after dissection, in order to restore its tensional equilibrium. Our results indicated that the value of the gain parameter, related to collagen G&R, plays an important role in the stability of aortic expansion after cTBAD. A small gain parameter will induce an excessive aneurysmal degeneration whilst a large gain parameter helps to recover a stabilized state of the artery after dissection. Finally, it was found that other mechanobiology-related parameters, such as the circumferential length of the dissection, as well as the pressure in the false lumen, may also be determinant for the stability of aneurysmal dilatation after cTBAD. Both a wide tear and an elevated false lumen pressure favor an unstable development of aortic expansion after cTBAD. As future work, the present model will be validated through predictions of aneurysmal dilatation in patient-specific clinical cases, in comparison with datasets followed over a significant period of time.

Aneurysm shrinkage is clinically observed after successful endovascular aortic aneurysm repair (EVAR). However, global understanding of post‐operative aneurysm evolutions remains weak. In this work, we propose to study these effects using numerical simulation. We set up a 3D finite‐element model of post‐EVAR vascular adaptation within an open‐source finite‐element code, which was initially developed for growth and remodeling (G&R). We modeled the endograft with a set of uniaxial prestrained springs that apply radial forces on the inner surface of the artery. Constitutive equations, momentum balance equations, and equations related to the mechanobiology of the artery were formulated based on the homogenized constrained mixture theory. We performed a sensitivity analysis by varying different selected parameters, namely oversizing and compliance of the stent‐graft, gain parameters related to collagen G&R, and the residual pressure in the aneurysm sac. This permitted us to evaluate how each factor influences post‐EVAR vascular adaptation. It was found that oversizing, compliance or gain parameters have a limited influence compared to that of the residual pressure in the aneurysm sac, which was found to play a critical role in the stability of aneurysm after stent‐graft implantation. An excessive residual pressure larger than 50 mmHg can induce a continuous expansion of the aneurysm while a moderate residual pressure below this critical threshold yields continuous shrinkage of the aneurysm. Moreover, it was found that elderly patients, with relatively lower amounts of remnant elastin in the arterial wall, are more sensitive to the effect of residual pressure. Therefore, these results show that elderly patients may present a higher potential risk of aortic sac expansion due to intra‐aneurysm sac pressure after EVAR than younger patients.

No abstract