Computational modeling of progressive damage and rupture in fibrous biological tissues: application to aortic dissection

Fracture mechanics modeling of aortic dissection

Aortic dissection, a critical cardiovascular condition with life-threatening implications, is distinguished by the development of a tear and its propagation within the aortic wall. A thorough understanding of the initiation and progression of these tears, or cracks, is essential for accurate diagnosis and effective treatment. This paper undertakes a fracture mechanics approach to delve into the mechanics of tear propagation in aortic dissection. Our objective is to elucidate the impact of geometric and material parameters, providing valuable insights into the determinants of this pivotal cardiovascular event. Through our investigation, we have gained an understanding of how various parameters influence the energy release rate for tear propagation in both longitudinal and circumferential directions, aligning our findings with clinical data.

"What makes blood clots break off?" A Back-of-the-Envelope Computation Toward Explaining Clot Embolization

PURPOSE

One in four deaths worldwide is due to thromboembolic disease; that is, one in four people die from blood clots first forming and then breaking off or embolizing. Once broken off, clots travel downstream, where they occlude vital blood vessels such as those of the brain, heart, or lungs, leading to strokes, heart attacks, or pulmonary embolisms, respectively. Despite clots' obvious importance, much remains to be understood about clotting and clot embolization. In our work, we take a first step toward untangling the mystery behind clot embolization and try to answer the simple question: "What makes blood clots break off?"

METHODS

To this end, we conducted experimentally-informed, back-of-the-envelope computations combining fracture mechanics and phase-field modeling. We also focused on deep venous clots as our model problem.

RESULTS

Here, we show that of the three general forces that act on venous blood clots-shear stress, blood pressure, and wall stretch-induced interfacial forces-the latter may be a critical embolization force in occlusive and non-occlusive clots, while blood pressure appears to play a determinant role only for occlusive clots. Contrary to intuition and prior reports, shear stress, even when severely elevated, appears unlikely to cause embolization.

CONCLUSION

This first approach to understanding the source of blood clot bulk fracture may be a critical starting point for understanding blood clot embolization. We hope to inspire future work that will build on ours and overcome the limitations of these back-of-the-envelope computations.

Uncertainty quantification of the wall thickness and stiffness in an idealized dissected aorta

Personalized treatment informed by computational models has the potential to markedly improve the outcome for patients with a type B aortic dissection. However, existing computational models of dissected walls significantly simplify the characteristic false lumen, tears and/or material behavior. Moreover, the patient-specific wall thickness and stiffness cannot be accurately captured non-invasively in clinical practice, which inevitably leads to assumptions in these wall models. It is important to evaluate the impact of the corresponding uncertainty on the predicted wall deformations and stress, which are both key outcome indicators for treatment optimization. Therefore, a physiology-inspired finite element framework was proposed to model the wall deformation and stress of a type B aortic dissection at diastolic and systolic pressure. Based on this framework, 300 finite element analyses, sampled with a Latin hypercube, were performed to assess the global uncertainty, introduced by 4 uncertain wall thickness and stiffness input parameters, on 4 displacement and stress output parameters. The specific impact of each input parameter was estimated using Gaussian process regression, as surrogate model of the finite element framework, and a δ moment-independent analysis. The global uncertainty analysis indicated minor differences between the uncertainty at diastolic and systolic pressure. For all output parameters, the 4th quartile contained the major fraction of the uncertainty. The parameter-specific uncertainty analysis elucidated that the material stiffness and relative thickness of the dissected membrane were the respective main determinants of the wall deformation and stress. The uncertainty analysis provides insight into the effect of uncertain wall thickness and stiffness parameters on the predicted deformation and stress. Moreover, it emphasizes the need for probabilistic rather than deterministic predictions for clinical decision making in aortic dissections.

Mathematical modeling and numerical simulation of arterial dissection based on a novel surgeon's view

This paper presents a mathematical model for arterial dissection based on a novel hypothesis proposed by a surgeon, Axel Haverich, see Haverich (Circulation 135(3):205-207, 2017. https://doi.org/10.1161/circulationaha.116.025407 ). In an attempt and based on clinical observations, he explained how three different arterial diseases, namely atherosclerosis, aneurysm and dissection have the same root in malfunctioning Vasa Vasorums (VVs) which are micro capillaries responsible for artery wall nourishment. The authors already proposed a mathematical framework for the modeling of atherosclerosis which is the thickening of the artery walls due to an inflammatory response to VVs dysfunction. A multiphysics model based on a phase-field approach coupled with mechanical deformation was proposed for this purpose. The kinematics of mechanical deformation was described using finite strain theory. The entire model is three-dimensional and fully based on a macroscopic continuum description. The objective here is to extend that model by incorporating a damage mechanism in order to capture the tearing (rupture) in the artery wall as a result of micro-injuries in VV. Unlike the existing damage-based model of the dissection in the literature, here the damage is driven by the internal bleeding (hematoma) rather than purely mechanical external loading. The numerical implementation is carried out using finite element method (FEM).

A review on the biomechanical behaviour of the aorta

Large aortic aneurysm and acute and chronic aortic dissection are pathologies of the aorta requiring surgery. Recent advances in medical intervention have improved patient outcomes; however, a clear understanding of the mechanisms leading to aortic failure and, hence, a better understanding of failure risk, is still missing. Biomechanical analysis of the aorta could provide insights into the development and progression of aortic abnormalities, giving clinicians a powerful tool in risk stratification. The complexity of the aortic system presents significant challenges for a biomechanical study and requires various approaches to analyse the aorta. To address this, here we present a holistic review of the biomechanical studies of the aorta by categorising articles into four broad approaches, namely theoretical, in vivo, experimental and combined investigations. Experimental studies that focus on identifying mechanical properties of the aortic tissue are also included. By reviewing the literature and discussing drawbacks, limitations and future challenges in each area, we hope to present a more complete picture of the state-of-the-art of aortic biomechanics to stimulate research on critical topics. Combining experimental modalities and computational approaches could lead to more comprehensive results in risk prediction for the aortic system.

An image-based approach for the estimation of arterial local stiffness in vivo

The analysis of mechanobiology of arterial tissues remains an important topic of research for cardiovascular pathologies evaluation. In the current state of the art, the gold standard to characterize the tissue mechanical behavior is represented by experimental tests, requiring the harvesting of ex-vivo specimens. In recent years though, image-based techniques for the in vivo estimation of arterial tissue stiffness were presented. The aim of this study is to define a new approach to provide local distribution of arterial stiffness, estimated as the linearized Young's Modulus, based on the knowledge of in vivo patient-specific imaging data. In particular, the strain and stress are estimated with sectional contour length ratios and a Laplace hypothesis/inverse engineering approach, respectively, and then used to calculate the Young's Modulus. After describing the method, this was validated by using a set of Finite Element simulations as input. In particular, idealized cylinder and elbow shapes plus a single patient-specific geometry were simulated. Different stiffness distributions were tested for the simulated patient-specific case. After the validation from Finite Element data, the method was then applied to patient-specific ECG-gated Computed Tomography data by also introducing a mesh morphing approach to map the aortic surface along the cardiac phases. The validation process revealed satisfactory results. In the simulated patient-specific case, root mean square percentage errors below 10% for the homogeneous distribution and below 20% for proximal/distal distribution of stiffness. The method was then successfully used on the three ECG-gated patient-specific cases. The resulting distributions of stiffness exhibited significant heterogeneity, nevertheless the resulting Young's moduli were always contained within the 1-3 MPa range, which is in line with literature.

An Alternative Digital Image Correlation-Based Experimental Approach to Estimate Fracture Parameters in Fibrous Soft Materials

One of the main challenges in experimental fracture mechanics is to correctly estimate fracture parameters of a nonhomogeneous and nonlinear material under large deformation. The crack tip detection is strongly affected by fibers at crack tip, leading to inaccurate measures. To overcome this limitation, a novel methodology based on the Digital Image Correlation (DIC) method for crack tip detection of fibrous soft composites is proposed in this work. The unidirectional composite was manufactured using a matrix of polydimethylsiloxane reinforced with a single layer of extensible cotton knit fabric. For two different fiber orientations, the crack growth (da), Crack Tip Opening Displacement (CTOD) and Crack Tip Opening Angle (CTOA) were determined using pure shear specimens under mode I fracture. A consistent estimation of fracture parameters was obtained. The location of the crack tip position during the fracture test using the DIC-based methodology was validated against a visual inspection approach. Results indicated that the DIC-based methodology is easily replicable, precise and robust.

Critical Pressure of Intramural Delamination in Aortic Dissection

Computational models of aortic dissection can examine mechanisms by which this potentially lethal condition develops and propagates. We present results from phase-field finite element simulations that are motivated by a classical but seldom repeated experiment. Initial simulations agreed qualitatively and quantitatively with data, yet because of the complexity of the problem it was difficult to discern trends. Simplified analytical models were used to gain further insight. Together, simplified and phase-field models reveal power-law-based relationships between the pressure that initiates an intramural tear and key geometric and mechanical factors-insult surface area, wall stiffness, and tearing energy. The degree of axial stretch and luminal pressure similarly influence the pressure of tearing, which was ~88 kPa for healthy and diseased human aortas having sub-millimeter-sized initial insults, but lower for larger tear sizes. Finally, simulations show that the direction a tear propagates is influenced by focal regions of weakening or strengthening, which can drive the tear towards the lumen (dissection) or adventitia (rupture). Additional data on human aortas having different predisposing disease conditions will be needed to extend these results further, but the present findings show that physiologic pressures can propagate initial medial defects into delaminations that can serve as precursors to dissection.

Simulating progressive intramural damage leading to aortic dissection using DeepONet: an operator-regression neural network

Aortic dissection progresses mainly via delamination of the medial layer of the wall. Notwithstanding the complexity of this process, insight has been gleaned by studying in vitro and in silico the progression of dissection driven by quasi-static pressurization of the intramural space by fluid injection, which demonstrates that the differential propensity of dissection along the aorta can be affected by spatial distributions of structurally significant interlamellar struts that connect adjacent elastic lamellae. In particular, diverse histological microstructures may lead to differential mechanical behaviour during dissection, including the pressure-volume relationship of the injected fluid and the displacement field between adjacent lamellae. In this study, we develop a data-driven surrogate model of the delamination process for differential strut distributions using DeepONet, a new operator-regression neural network. This surrogate model is trained to predict the pressure-volume curve of the injected fluid and the damage progression within the wall given a spatial distribution of struts, with in silico data generated using a phase-field finite-element model. The results show that DeepONet can provide accurate predictions for diverse strut distributions, indicating that this composite branch-trunk neural network can effectively extract the underlying functional relationship between distinctive microstructures and their mechanical properties. More broadly, DeepONet can facilitate surrogate model-based analyses to quantify biological variability, improve inverse design and predict mechanical properties based on multi-modality experimental data.

A Review on Damage and Rupture Modelling for Soft Tissues

Computational modelling of damage and rupture of non-connective and connective soft tissues due to pathological and supra-physiological mechanisms is vital in the fundamental understanding of failures. Recent advancements in soft tissue damage models play an essential role in developing artificial tissues, medical devices/implants, and surgical intervention practices. The current article reviews the recently developed damage models and rupture models that considered the microstructure of the tissues. Earlier review works presented damage and rupture separately, wherein this work reviews both damage and rupture in soft tissues. Wherein the present article provides a detailed review of various models on the damage evolution and tear in soft tissues focusing on key conceptual ideas, advantages, limitations, and challenges. Some key challenges of damage and rupture models are outlined in the article, which helps extend the present damage and rupture models to various soft tissues.

On a phase-field approach to model fracture of small intestine walls

We address anisotropic elasticity and fracture in small intestine walls (SIWs) with both experimental and computational methods. Uniaxial tension experiments are performed on porcine SIW samples with varying alignments and quantify their nonlinear elastic anisotropic behavior. Fracture experiments on notched SIW strips reveal a high sensitivity of the crack propagation direction and the failure stress on the tissue orientation. From a modeling point of view, the observed anisotropic elastic response is studied with a continuum mechanical model stemming from a strain energy density with a neo-Hookean component and an anisotropic component with four families of fibers. Fracture is addressed with the phase-field approach, featuring two-fold anisotropy in the fracture toughness. Elastic and fracture model parameters are calibrated based on the experimental data, using the maximum and minimum limits of the experimental stress-stretch data set. A very good agreement between experimental data and computational results is obtained, the role of anisotropy being effectively captured by the proposed model in both the elastic and the fracture behavior. STATEMENT OF SIGNIFICANCE: This article reports a comprehensive experimental data set on the mechanical failure behavior of small intestinal tissue, and presents the corresponding protocols for preparing and testing the samples. On the one hand, the results of this study contribute to the understanding of small intestine mechanics and thus to understanding of load transfer mechanisms inside the tissue. On the other hand, these results are used as input for a phase-field modelling approach, presented in this article. The presented model can reproduce the mechanical failure behavior of the small intestine in an excellent way and is thus a promising tool for the future mechanical description of diseased small intestinal tissue.

Differential propensity of dissection along the aorta

Aortic dissections progress, in part, by delamination of the wall. Previous experiments on cut-open segments of aorta demonstrated that fluid injected within the wall delaminates the aorta in two distinct modes: stepwise progressive tearing in the abdominal aorta and a more prevalent sudden mode of tearing in the thoracic aorta that can also manifest in other regions. A microstructural understanding that delineates these two modes of tearing has remained wanting. We implemented a phase-field finite-element model of the aortic wall, motivated in part by two-photon imaging, and found correlative relations for the maximum pressure prior to tearing as a function of local geometry and material properties. Specifically, the square of the pressure of tearing relates directly to both tissue stiffness and the critical energy of tearing and inversely to the square root of the torn area; this correlation explains the sudden mode of tearing and, with the microscopy, suggests a mechanism for progressive tearing. Microscopy also confirmed that thick interlamellar radial struts are more abundant in the abdominal region of the aorta, where progressive tearing was observed previously. The computational results suggest that structurally significant radial struts increase tearing pressure by two mechanisms: confining the fluid by acting as barriers to flow and increasing tissue stiffness by holding the adjacent lamellae together. Collectively, these two phase-field models provide new insights into the mechanical factors that can influence intramural delaminations that promote aortic dissection.

A framework for incorporating 3D hyperelastic vascular wall models in 1D blood flow simulations

We present a novel framework for investigating the role of vascular structure on arterial haemodynamics in large vessels, with a special focus on the human common carotid artery (CCA). The analysis is carried out by adopting a three-dimensional (3D) derived, fibre-reinforced, hyperelastic structural model, which is coupled with an axisymmetric, reduced order model describing blood flow. The vessel transmural pressure and lumen area are related via a Holzapfel–Ogden type of law, and the residual stresses along the thickness and length of the vessel are also accounted for. After a structural characterization of the adopted hyperelastic model, we investigate the link underlying the vascular wall response and blood-flow dynamics by comparing the proposed framework results against a popular tube law. The comparison shows that the behaviour of the model can be captured by the simpler linear surrogate only if a representative value of compliance is applied. Sobol’s multi-variable sensitivity analysis is then carried out in order to identify the extent to which the structural parameters have an impact on the CCA haemodynamics. In this case, the local pulse wave velocity (PWV) is used as index for representing the arterial transmission capacity of blood pressure waveforms. The sensitivity analysis suggests that some geometrical factors, such as the stress-free inner radius and opening angle, play a major role on the system’s haemodynamics. Subsequently, we quantified the differences in haemodynamic variables obtained from different virtual CCAs, tube laws and flow conditions. Although each artery presents a distinct vascular response, the differences obtained across different flow regimes are not significant. As expected, the linear tube law is unable to accurately capture all the haemodynamic features characterizing the current model. The findings from the sensitivity analysis are further confirmed by investigating the axial stretching effect on the CCA fluid dynamics. This factor does not seem to alter the pressure and flow waveforms. On the contrary, it is shown that, for an axially stretched vessel, the vascular wall exhibits an attenuation in absolute distension and an increase in circumferential stress, corroborating the findings of previous studies. This analysis shows that the new model offers a good balance between computational complexity and physics captured, making it an ideal framework for studies aiming to investigate the profound link between vascular mechanobiology and blood flow.

A Parametric Study on Factors Influencing the Onset and Propagation of Aortic Dissection Using the Extended Finite Element Method

OBJECTIVE

Aortic dissection is a life-threatening event which starts most of the time with an intimal tear propagating along the aortic wall, while blood enters the medial layer and delaminates the medial lamellar units. Studies investigating the mechanisms underlying the initiation sequence of aortic dissection are rare in the literature, the majority of studies being focused on the propagation event. Numerical models can provide a deeper understanding of the phenomena involved during the initiation and the propagation of the initial tear, and how geometrical and mechanical parameters affect this event. In the present paper, we investigated the primary factors contributing to aortic dissection.

METHODS

A two-layer arterial model with an initial tear was developed, representing three different possible configurations depending on the initial direction of the tear. Anisotropic damage initiation criteria were developed based on uniaxial and shear experiments from the literature to predict the onset and the direction of crack propagation. We used the XFEM-based cohesive segment method to model the initiation and the early propagation of the tear along the aorta. A design of experiment was used to quantify the influence of 7 parameters reflecting crack geometry and mechanics of the wall on the critical pressure triggering the dissection and the directions of propagation of the tear.

RESULTS

The results showed that the obtained critical pressures (mean range from 206 to 251 mmHg) are in line with measurement from the literature. The medial tensile strength was found to be the most influential factor, suggesting that a medial degeneration is needed to reach a physiological critical pressure and to propagate a tear in an aortic dissection. The geometry of the tear and its location inside the aortic wall were also found to have an important role not only in the triggering of tear propagation, but also in the evolution of the tear into either aortic rupture or aortic dissection. A larger and deeper initial tear increases the risk of aortic dissection.

CONCLUSION

The numerical model was able to reproduce the behaviour of the aorta during the initiation and propagation of an aortic dissection. In addition to confirm multiple results from the literature, different types of tears were compared and the influence of several geometrical and mechanical parameters on the critical pressure and direction of propagation was evaluated with a parametric study for each tear configuration.

SIGNIFICANCE

Although these results should be experimentally validated, they allow a better understanding of the phenomena behind aortic dissection and can help in improving the diagnosis and treatment of this disease.

Quantification of true lumen helical morphology and chirality in type B aortic dissections

Chirality is a fundamental property in many biological systems. Motivated by previous observations of helical aortic blood flow, aortic tissue fibers, and propagation of aortic dissections, we introduce methods to characterize helical morphology of aortic dissections. After validation on computer-generated phantoms, the methods were applied to patients with type B dissection. For this cohort, there was a distinct bimodal distribution of helical propagation of the dissection with either achiral or exclusively right-handed chirality, with no intermediate cases or left-handed cases. This clear grouping indicates that dissection propagation favors these two modes, which is potentially due to the right-handedness of helical aortic blood flow and cell orientation. The characterization of dissection chirality and quantification of helical morphology advances our understanding of dissection pathology and lays a foundation for applications in clinical research and treatment practice. For example, the chirality and magnitude of helical metrics of dissections may indicate risk of dissection progression, help define treatment and surveillance strategies, and enable development of novel devices that account for various helical morphologies.NEW & NOTEWORTHY A novel definition of helical propagation of type B aortic dissections reveals a distinct bimodality, with the true lumen being either achiral (nonhelical) or exclusively right-handed. This right-handed chirality is consistent with anatomic and physiological phenomena such as right-handed twist during left ventricle contraction, helical blood flow, and tissue fiber direction. The helical character of aortic dissections may be useful for pathology research, diagnostics, treatment selection, therapeutic durability prediction, and aortic device design.

Multimodality Imaging-Based Characterization of Regional Material Properties in a Murine Model of Aortic Dissection

Chronic infusion of angiotensin-II in atheroprone (ApoE-/-) mice provides a reproducible model of dissection in the suprarenal abdominal aorta, often with a false lumen and intramural thrombus that thickens the wall. Such lesions exhibit complex morphologies, with different regions characterized by localized changes in wall composition, microstructure, and properties. We sought to quantify the multiaxial mechanical properties of murine dissecting aneurysm samples by combining in vitro extension-distension data with full-field multimodality measurements of wall strain and thickness to inform an inverse material characterization using the virtual fields method. A key advance is the use of a digital volume correlation approach that allows for characterization of properties not only along and around the lesion, but also across its wall. Specifically, deformations are measured at the adventitial surface by tracking motions of a speckle pattern using a custom panoramic digital image correlation technique while deformations throughout the wall and thrombus are inferred from optical coherence tomography. These measurements are registered and combined in 3D to reconstruct the reference geometry and compute the 3D finite strain fields in response to pressurization. Results reveal dramatic regional variations in material stiffness and strain energy, which reflect local changes in constituent area fractions obtained from histology but emphasize the complexity of lesion morphology and damage within the dissected wall. This is the first point-wise biomechanical characterization of such complex, heterogeneous arterial segments. Because matrix remodeling is critical to the formation and growth of these lesions, we submit that quantification of regional material properties will increase the understanding of pathological mechanical mechanisms underlying aortic dissection.

An adaptive finite element model for steerable needles

Penetration of a flexible and steerable needle into a soft target material is a complex problem to be modelled, involving several mechanical challenges. In the present paper, an adaptive finite element algorithm is developed to simulate the penetration of a steerable needle in brain-like gelatine material, where the penetration path is not predetermined. The geometry of the needle tip induces asymmetric tractions along the tool–substrate frictional interfaces, generating a bending action on the needle in addition to combined normal and shear loading in the region where fracture takes place during penetration. The fracture process is described by a cohesive zone model, and the direction of crack propagation is determined by the distribution of strain energy density in the tissue surrounding the tip. Simulation results of deep needle penetration for a programmable bevel-tip needle design, where steering can be controlled by changing the offset between interlocked needle segments, are mainly discussed in terms of penetration force versus displacement along with a detailed description of the needle tip trajectories. It is shown that such results are strongly dependent on the relative stiffness of needle and tissue and on the tip offset. The simulated relationship between programmable bevel offset and needle curvature is found to be approximately linear, confirming empirical results derived experimentally in a previous work. The proposed model enables a detailed analysis of the tool–tissue interactions during needle penetration, providing a reliable means to optimise the design of surgical catheters and aid pre-operative planning.

Modeling lamellar disruption within the aortic wall using a particle-based approach

Aortic dissections associate with medial degeneration, thus suggesting a need to understand better the biophysical interactions between the cells and matrix that constitute the middle layer of the aortic wall. Here, we use a recently extended "Smoothed Particle Hydrodynamics" formulation to examine potential mechanisms of aortic delamination arising from smooth muscle cell (SMC) dysfunction or apoptosis, degradation of or damage to elastic fibers, and pooling of glycosaminoglycans (GAGs), with associated losses of medial collagen in the region of the GAGs. First, we develop a baseline multi-layered model for the healthy aorta that delineates medial elastic lamellae and intra-lamellar constituents. Next, we examine stress fields resulting from the disruption of individual elastic lamellae, lost SMC contractility, and GAG production within an intra-lamellar space, focusing on the radial transferal of loading rather than on stresses at the tip of the delaminated tissue. Results suggest that local disruptions of elastic lamellae transfer excessive loads to nearby intra-lamellar constituents, which increases cellular vulnerability to dysfunction or death. Similarly, lost SMC function and accumulations of GAGs increase mechanical stress on nearby elastic lamellae, thereby increasing the chance of disruption. Overall these results suggest a positive feedback loop between lamellar disruption and cellular dropout with GAG production and lost medial collagen that is more pronounced at higher distending pressures. Independent of the initiating event, this feedback loop can catastrophically propagate intramural delamination.